CC BY
45
УДК 004.942
МОДЕЛЬ ПРОГНОЗИРОВАНИЯ ХАРАКТЕРИСТИК НАВОДНЕНИЯ, ВЛИЯЮЩИХ НА ВЕЛИЧИНУ ПРИЧИНЯЕМОГО УЩЕРБА
М. В. Хамутова1, В. А. Кушников2
1 Хамутова Мария Васильевна, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, Саратов, Астраханская, 83, mariuka7d@rambler.ru
2Кушников Вадим Алексеевич, доктор технических наук, директор, Институт проблем точной механики и управления РАН, Россия, 410024, Саратов, Рабочая, 24, kushnikoff@yandex.ru
В соответствии с ГОСТ выбраны характеристики наводнений, влияющие на величину причиняемого им ущерба, прогнозирование которых позволит повысить эффективность ликвидации последствий наводнений. На основе формального аппарата системной динамики, который учитывает причинно-следственные связи между моделируемыми переменными, разработана математическая модель, позволяющая прогнозировать характеристики наводнений. Построен граф причинно-следственных связей, существующих между моделируемыми характеристиками наводнений. Математическая модель для прогнозирования характеристик наводнений описывается системой нелинейных дифференциальных уравнений первого порядка. Представлены функциональные зависимости правых частей системы уравнений, описывающие математическую модель. Численное решение системы уравнений получено с помощью метода Рунге-Кутты. Проведены вычислительные эксперименты, позволяющие на различных временных интервалах и с учетом изменяющихся параметров внешней среды определить моделируемые характеристики. Сравнение прогнозируемыххарактеристик, рассчитанных по модели с реальными значениями наводнения, произошедшего в Приморье в августе 2001 г., подтверждает адекватность математической модели. Полученные результаты могут быть использованы при разработке информационных систем прогнозирования последствий наводнения для оперативно-диспетчерского персонала МЧС.
Ключевые слова: математическая модель, системная динамика, прогнозирование характеристик наводнений.
001:10.18500/1816-9791-2017-17-2-231-238
ВВЕДЕНИЕ
Наводнения встречаются практически по всей территории России, причиняют материальный ущерб, наносят вред здоровью населения и приводят к гибели людей. Для повышения эффективности ликвидации последствий наводнений необходимо спрогнозировать в режиме чрезвычайной ситуации характеристики наводнения, влияющие на величину ущерба. Существующие модели прогнозирования не позволяют определить совокупность характеристик наводнения с учетом большого количества нелинейных обратных связей между ними [1,2]. Данное обстоятельство проводит к
уменьшению точности прогнозирования, что негативно влияет на процессы ликвидации последствий наводнений в условиях чрезвычайной ситуации.
В статье приведена математическая модель для прогнозирования характеристик наводнений, влияющих на величину ущерба. Предложен и обоснован ориентированный граф причинно-следственных связей, существующий между прогнозируемыми характеристиками наводнения. При разработке математической модели использован формальный аппарат системной динамики, учитывающий причинно-следственные связи между моделируемыми переменными. Результаты вычислительного эксперимента подтверждают адекватность математической модели.
1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
В соответствии с ГОСТ1 в качестве прогнозируемых характеристик наводнения были выбраны: Х1 — численность группировки сил, участвующих в аварийно-спасательных работах; Х2 — количество жилых домов, разрушенных и поврежденных в результате наводнения; Х3 — численность населения, эвакуированного из зоны затопления; Х4 — количество погибших; Х5 — протяженность железных и автомобильных дорог, оказавшихся в зоне затопления; Хб — количество промышленных предприятий в зоне наводнения; Х7 — количество транспортных средств, участвующих в аварийно-спасательных работах; Х8 — численность населения в зоне затопления; Х9 — площадь сельскохозяйственных угодий, охваченных наводнением; Хю — количество погибших сельскохозяйственных животных.
Математическая модель прогнозирования характеристик наводнения, влияющих на величину ущерба, разрабатывается на основе аппарата системной динамики, где при описании моделируемых переменных используется система нелинейных дифференциальных уравнений 1-го порядка:
= Х+ + Х-, % = 1ТП, (1)
где Х+, Х- , % = 1 , п — непрерывные или кусочно-непрерывные функции, определяющие положительную и отрицательную скорость изменения переменной Х^ [3]. В свою очередь, Х+ = , ,..., ), Х- = /-(^1, ,..., Fm) — функции от факторов Fj, у = 1,т , влияющих на скорость изменения переменной, при этом факторы могут быть функциями от прогнозируемых переменных Х^, % = 1 , п [4,5].
На основании анализа причинно-следственных связей, существующих между прогнозируемыми характеристиками наводнений, разработан ориентированный граф, представленный на рис. 1. Моделируемые характеристики Х^ являются вершинами графа, а исходящие и входящие ребра характеризуют функциональные связи между ними.
1ГОСТ 22.0.06-97/ГОСТ Р 22.0.06-95. Безопасность в чрезвычайных ситуациях. Источники природных чрезвычайных ситуаций. Поражающие факторы. Номенклатура параметров поражающих воздействий. Принят постановлением Госстандарта РФ № 308 от 20.06.1995 г.
Рис. 1. Граф причинно-следственных связей Fig. 1. The graph of cause-effect relations
Учитывая ориентированный граф причинно-следственных связей между прогнозируемыми характеристиками и параметрами внешней среды, разработанная математическая модель имеет вид
' dX1(t) dt = /i+ S (t),X (t)),
dX2(t) dt = f2+ F (t), G(t), S (t), X (t), t) - /- (Xx (t), X7 (t)),
dXs(t) dt = /з+ X8(i),X1(i),X7 (t)),
dX4(t) dt = /4+ F (t),G(t),T (t),X (t),X7 (t),Xi (t)),
dXe(t) dt = /5+ A(t), S (t)) - /— (Xi(t),X7(t)),
dXe(t) dt = /б+ S (t),X8 (t)),
dXr(t) dt = /7+ Xi(t)),
dXs(t) dt = /8+ D(t),S(t)) - /— (X4),
dXg(t) dt = /9+ I (t), S (t)) - /- (Xi(t),X7(t)),
dX io (t) dt = Ло + (F (t), G(t), T (t), S (t),Xi (t), X7 (t)),
(2)
где A(t) — плотность транспортных сетей в зоне затопления, D(t) — плотность населения, F(t), G(t), T(t) — средняя скорость течения, глубина и температура воды соответственно, I(t) — доля площади сельскохозяйственных угодий, S(t) — площадь зоны затопления.
Зависимости f— и f+ определяются на основе анализа статистических данных и аппроксимируются полиномами. При отсутствии или недостоверности статистических данных, используемых при построении модели (2), эти зависимости могут быть определены исходя из физического смысла моделируемого процесса [6,7]. Площадь зоны затопления S(t), средняя скорость течения F(t), средняя глубина воды G(t) и средняя температура воды T(t) рассчитываются экспертами по известным методикам. В частности, для определения площади зоны затопления могут использоваться методы, позволяющие моделировать динамику поверхностных вод на произвольном рельефе местности.
2. ПРОВЕРКА АДЕКВАТНОСТИ МОДЕЛИ
Проверка адекватности модели осуществляется посредством сравнения характеристик, определенных по модели (2), с реальными характеристиками наводнения, происшедшего в августе 2001 г. в Приморье. В таблице представлены функции f— и f+ системы (2), построенные для данного наводнения по методикам, представленным в [6,7].
Аналитический вид функций /— и /+ The analytic form of the functions f— and /+
/i+ ifci^S(t)X8, S(t) >e, 1 0, S(t) < e, /6+ ifc6S(t)0.5X80.i5 S(t) > e, (0, S (t) < e,
/2 + ifc2F(t)G(t)t^S(t)X8, S(t) > e, \0, S(t) < e, /7+ k7Xi
/з + k3 ^X1 /8 + ifc8D(t)S(t), S(t) > e, (0, S (t) < e,
/4 + k F(t)G(t)T(t)X8 4 X7 Xi /9 + jk91(t)S(t), S(t) > e, |0, S (t) < e,
/5 + ifc5A(t)S(t), S(t) > e, |0,S(t) < e, /io+ ikioF(t)G£™S(t), S(t) > e, (0, S (t) < e,
/2 - kii XiX7 /5 - ki2Xi X7
/8 - ki3 X4 /9 - ki4Xi X7
В частности, функция /+(£(£),Х8(£)), характеризующая скорость изменения переменной X, выбрана прямо пропорционально площади зоны затопления и численности населения в зоне затопления, что соответствует экспериментальным данным.
Для нормирования произведения факторов 5(¿) и Х8 (¿) используется возведение в степень 0.5, что обеспечивает наиболее точное совпадение с реальными результатами. Таким образом, искомая зависимость имеет вид
/,+(s<.>.*«»иS(•>
0, S(t> < е.
Коэффициенты к, г = 1,14, определяются посредством вычислительного эксперимента на этапе адаптации модели к объекту исследования. Система дифференциальных уравнений (2) с учетом данных, приведенных в таблице, была решена численным методом Рунге-Кутты 4-го порядка при начальных условиях ¿0 = 1, X (¿0) = Х^0, г = 1710.
На рис. 2 представлено решение системы уравнений (2), нормированное относительно максимальных значений моделируемых характеристик: Х^ = X/Хгтах.
Рис. 2. Результаты прогнозирования характеристик наводнения Fig. 2. Results of forecasting flood characteristics
В частности, анализ зависимости X^(t> позволяет сделать вывод, что максимальные потери среди населения приходятся на первые сутки наводнения и в дальнейшем их рост незначителен, что подтверждается экспертами. Увеличение остальных прогнозируемых характеристик обусловлено ростом площади зоны затопления, пик которого приходится на 4-й день наводнения.
На рис. 3 выполнено сравнение характеристики X^(t>, определенной по модели (2), с реальными статистическими данными, интерполированными многочленом Лагранжа Y4(t>.
Рис. 3. Сравнение значений моделируемой характеристики XN(t)
с реальными данными Fig. 3. Comparison of the values of the simulated characteristic XN(t)
with real data
Анализ полученных графиков позволяет утверждать, что значения X± (£) лишь
_ X N
незначительно отличаются от реальных данных в узлах £ = 1,4 Аср4 ^ 3%. При сравнении остальных моделируемых характеристик модели (2) с их реальными значениями, среднее значение относительных погрешностей в узлах моделирования для каждой характеристики не превышает 10%, что позволяет утверждать, что разработанная математическая модель является адекватной [8,9].
Библиографический список
1. Авакян А. Б. Наводнения. Концепция защиты // Изв. РАН. Сер. географическая. 2000. № 5. С. 40-46.
2. Воронин А. А., Васильченко А. А., Храпов С. С., Агафонникова Е. О. Анализ эффективности природовосстановительных проектов в эколого-экономической системе «Волжская ГЭС — Волго-Ахтубинская пойма» // Управление большими системами. М. : ИПУ РАН, 2014. Вып. 52. С. 133-147.
3. Садовничий В. А., Акаев А. А., Коротаев А. В., Малков С. Ю. Моделирование и прогнозирование мировой динамики / Научный совет по Программе фунд. исслед. Президиума РАН. М. : ИСПИ РАН, 2012. 360 с. Экономика и социология знания.
4. Бродский Ю. И. Лекции по математическому и имитационному моделированию. М. ; Берлин : Директ-Медиа, 2015. 240 с.
5. Белотелов Н. В., Бродский Ю. И., Павловский Ю. Н. Сложность. Математическое моделирование. Гуманитарный анализ : Исследование исторических, военных, социально-экономических и политических процессов. М. : Кн. дом «ЛИБРОКОМ», 2009. 320 с.
6. Федянин В. И., Проскурников Ю. Е. Организация и ведение аварийно-спасательных и других неотложных работ при ликвидации чрезвычайных ситуаций природного характера : учеб. пособие. Воронеж : Воронеж. гос. техн. ун-т, 2006. Ч. 1. 469 с.
7. РД 153-34.2-002-01. Временная методика оценки ущерба, возможного вследствие аварии гидротехнического сооружения (приказ Минэнерго России № 130 от 26.04.2001 г.). URL: https://znaytovar.ru/gost/2/RD_15334200201_Vremennaya_meto.html (дата обращения: 01.02.2017).
8. Хамутова М. В., Кушников В. А. Математическая модель прогнозирования последствий наводнений // Вестн. Астрахан. гос. техн. ун-та. Сер. Управление, вычислительная техника и информатика. 2016. № 3. С. 109-114.
9. Клюев В. В., Резчиков А. Ф., Кушников В. А., Богомолов А. С., Иващенко В. А., Фили-монюк Л. Ю., Хамутова М. В. Информационно-управляющая система для поддержки принятия решения по ликвидации последствий наводнения // Вестн. компьютерных и информационных технологий. 2016. № 11 (149). С. 39-45.
Образец для цитирования:
Хамутова М. В., Кушников В. А. Модель прогнозирования характеристик наводнения, влияющих на величину причиняемого ущерба // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2017. Т. 17, вып. 2. С. 231-238. DOI: 10.18500/1816-9791-2017-172-231-238.
A Model for Forecasting Characteristics of Floods Affecting the Value
of the Caused Damage
M. V. Khamutova1, V. A. Kushnikov2
1 Maria V. Khamutova, ORCID: 0000-0002-7926-1347, Saratov State University, 83, Astrakhanskaya str., 410012, Saratov, Russia, mariuka7d@rambler.ru
2Vadim A. Kushnikov, Institute of Precision Mechanics and Control, Russian Academy of Sciences, 24, Rabochayastr., 410028, Saratov, Russia, iptmuran@san.ru
According to GOST (State Standard) characteristics of floods are selected that have influence on the extent of damage. Forecasting of such characteristics will improve the efficiency of liquidation of floods effects. On the basis of formal apparatus of system dynamics that takes into account the casual relationships between modeled variables a mathematical model is developed to forecast characteristics of floods. A graph of cause-effect relations that exist between modeled characteristics is constructed. The mathematical model is described by a system of first order nonlinear differential equations. The numerical solution of the system was obtained with the aid of Runge-Kutta method. Computational experiments were conducted that allowed to determine the modeled characteristics for different time intervals. Comparison of the characteristics calculated according to the model with their real values during the Primorski Krai's flood in August 2001 confirms the adequacy of the mathematical model. The results obtained can be used in the development of information systems of forecasting the flood effects for operative and dispatching personnel of the EMERCOM.
Key words: mathematical model, system dynamics, forecasting the characteristics of floods. References
1. Avakjan A. B. Navodnenija. Koncepcija zashhity [Floods. The concept of protection]. Izvestiya RAN. Ser. Geograficheskaya [Izvestiya RAS. Ser. Geographic], 2000, no. 5, pp. 40-46 (in Russian).
2. Voronin A., Vasilchenko A., Chrapov S., Agafonnikova E. Efficiency analysis for nature restoration projects in ecologo-economic system of Volga HPP and Volga-Akhtuba Flood-plain. Large-scale Systems Control. Moscow, IPU RAN, 2014, iss. 52, pp. 133-147 (in Russian).
3. Sadovnichiy V. A., Akayev A. A., Korotayev A. V., Malkov S. Yu. Modelling and Forecasting World Dynamics. Scientific Council for ""Economics and Sociology of Knowledge" Fundamental Research Programme of the Presidium of the RAS. Moscow, RAS ISPR, 2012.
4. Brodskij Ju. I. Lekcii po matematicheskomu i imitacionnomu modelirovaniju [Lectures on mathematical and simulation modeling]. Moscow, Berlin, Direkt-Media, 2015. 240 p. (in Russian).
5. Belotelov N. V., Brodskij Ju. I., Pavlovskij Ju. N. Complexity. Mathematical modeling. Humanitarian analysis: the Study of the historical, military, socio-economic and political processes. Mosow, LIBROKOM Publ., 2009. 320 p. (in Russian).
6. Fedianin V. I., Proskurnikov Yu. E. Organization and imple-mentation of rescue and other emergency operations when liquidating the emergent situations of natural origin. Voronezh, Voronezh Techn. Univ. Press, 2006, pt. 1. 469 p. (in Russian).
7. RD 153-34.2-002-01. Temporary methods of assessment of damage, possible due to an accident waterworks (the Order of the Ministry of Energy of Russia no. 130, 26, April, 2001) (in Russian). Available at: https://znaytovar.ru/gost/2/ RD_15334200201_Vremennaya_meto.html (accessed 01, February, 2017)
8. Khamutova M. V., Kushnikov V. A. Mathematical model for forecasting of flood effects. Vestnik Astrakhan State Technical Univ. Ser. Management, Computer Sciences and Informatics, 2016, no. 3, pp. 109-114 (in Russian).
9. Klyuev V. V., Rezchikov A. F., Kushnikov V. A., Bogomolov A. S., Ivashchenko V. A., Filimonyuk L. Yu., Khamutova M. V. Informational-control system for decision-making supply during elimination of floods' consequences. Vestnik komp'iuternykh i informatsion-nykh tekhnologii [Herald of computer and information technologies], 2016, no. 11(149), pp. 39-45 (in Russian). DOI: 10.14489/vkit.2016.11.pp.039-045.
Cite this article as:
Khamutova M. V., Kushnikov V. A. A Model for Forecasting Characteristics of Floods Affecting the Value of the Caused Damage. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2017, vol. 17, iss. 2, pp. 231-238 (in Russian). DOI: 10.18500/1816-9791-2017-17-2-231-238.
