У таблиці 1 наведені середньомісячні показнники спостережень Державної гідрометеорологічної служби в Івано-Франківській області, № 109/04-18, від 12.03.08 року, за стоком ріки Дністер на протязі 1995 – 2007 р.р. на гідрологічному посту в м. Галичі.
Таблиця 1- Середньомісячні стоки р. Дністер,
.місяцьрік | I | II | III | IV | V | VI | VII | VIII | IX | X | XI | XII |
1995 | 107 | 155 | 199 | 315 | 206 | 181 | 96,3 | 35,4 | 78,4 | 51,6 | 116 | 69,0 |
1996 | 81.0 | 51,2 | 53,2 | 607 | 265 | 95,3 | 82,8 | 164 | 364 | 205 | 127 | 119 |
1997 | 69.6 | 230 | 156 | 271 | 336 | 256 | 222 | 226 | 128 | 181 | 135 | 182 |
1998 | 187 | 213 | 156 | 417 | 315 | 544 | 507 | 206 | 139 | 248 | 357 | 106 |
1999 | 171 | 113 | 542 | 485 | 204 | 161 | 157 | 126 | 105 | 106 | 73,7 | 191 |
2000 | 70.0 | 315 | 301 | 538 | 91,9 | 52,5 | 94,1 | 69,9 | 69,2 | 39,2 | 37,3 | 56,0 |
2001 | 66.8 | 134 | 368 | 182 | 98,7 | 430 | 408 | 261 | 254 | 105 | 189 | 92,6 |
2002 | 218 | 300 | 245 | 200 | 127 | 170 | 90,8 | 155 | 127 | 217 | 150 | 64,4 |
2003 | 58.4 | 51,9 | 236 | 263 | 120 | 63,7 | 71,4 | 36,6 | 44,5 | 95,3 | 116 | 48,3 |
2004 | 147 | 180 | 257 | 133 | 108 | 68,6 | 84,2 | 415 | 115 | 135 | 216 | 154 |
2005 | 96.8 | 75,8 | 316 | 434 | 381 | 226 | 85,2 | 210 | 84,3 | 92,6 | 53,1 | 75,8 |
2006 | 82.9 | 120 | 251 | 599 | 187 | 375 | 159 | 179 | 108 | 48,9 | 128 | 52,6 |
2007 | 138 | 232 | 248 | 84,0 | 102 | 79,1 | 56,0 | 51,9 | 372 | 131 | 151 | 146 |
Для прогнозу стоку р. Дністер на 2008 рік був вибраний поліном
, (6)де
- середньорічний стік р. Дністер; - середньомісячний стік р. Дністер.Таким чином, прогнозування стоку на майбутній
, місяць здійснюється за середнім стоком за минулий місяць.Для пошуку закономірності перерозподілу середньомісячного стоку ріки за роками необхідно ідентифікувати
регресійних залежностей типу (6). Для випадку, що розглядається .Структура залежності (6) невідома, тому з використанням комбінаторного алгоритму методу групового врахуванням аргументів визначається повний набір поліномів для кожного із дванадцяти місяців за даними спостережень, які наведені у табл. 1.
Суть комбінаторного методу у тому, що здійснюється повний перебір поліномів типу (6) шляхом почергового обнуління його коефіцієнтів. Загальне число таких поліномів
.Оскільки регресійна модель (5) є функцією двох змінних, то у відповідності з формулою (2)
і число можливих часткових моделей для кожного місяця.При реалізації алгоритму прогнозування помісячного стоку табл. 1 була розбита на три частини – навчальну
, перевірочну та екзаменаційну . Були вибрані такі значення: , та . Для кожного із дванадцяти місяців генерувалось ( ) моделей, параметри яких знаходились за методом найменших квадратів на множині експериментальних точок . Із них за критерієм мінімуму зсуву було відібрано п’ять моделей. Для відібраних моделей їх параметри уточнювались на множинах і . Потім складаються всі можливі поєднання із дванадцяти поліномів так, щоб у них був присутнім лише один поліном із набору. Для всіх поєднань обчислювалось значення критерію балансу (5). Сукупність дванадцяти регресійних моделей типу (6), для яких значення критерію балансу (5) мінімальне, приймалась як модель для прогнозу стоку р. Дністер.Табл. 2 вміщує коефіцієнти моделей оптимальної складності, відібрані за критерієм балансу.
Таблиця 2 – Параметри моделей оптимальної складності для прогнозування стоку ріки Дністер
1 | 2 | 3 | 4 | 5 | 6 | ||||||||
0.1674 | 0.0239 | 0.0773 | 1.5619 | 0.3581 | 0.1747 | ||||||||
0 | 0.9063 | 0 | -3.6525 | 0 | -1.9839 | ||||||||
-1.1085 | 0.9567 | 3.8726 | -2.982 | -0.9511 | 0 | ||||||||
-0.577 | -1.2033 | 3.9814 | 5.9633 | -2.8219 | 9.7877 | ||||||||
4.522 | -1.5683 | -11.0133 | 0 | 4.455 | -1.0436 | ||||||||
0 | 0 | -2.7274 | 3.1292 | 0 | 0 | ||||||||
-0.1727 | -0.6039 | -0.4754 | -0.159 | 0.2354 | -0.2454 | ||||||||
2.205 | 6.0768 | 4.6324 | 1.4461 | -1.0732 | 1.5784 | ||||||||
0 | 0 | 0 | 0.7552 | 0 | 0.8393 | ||||||||
-4.6749 | -10.435 | -7.0557 | 0 | 0 | 0 | ||||||||
0.7235 | 0 | 0 | -3.948 | 9.4151 | -3.2964 | ||||||||
0.8397 | 0.313 | -0.104 | 0.6465 | -5.3774 | 0 |
На екзаменаційній множині
провірялась точність алгоритму прогнозування стоку р. Дністер. Результат такої перевірки відтворює рис. 2.Рисунок 2. Помісячний прогноз стоку р. Дністер
Аналіз отриманих результатів показує, що найважче прогнозувати пікові значення стоків під час весняних повеней, які припадають на березень – квітень. На цьому ж рисунку зроблений прогноз помісячного стоку р. Дністер на 2008 р. Оцінка точності методу помісячного прогнозу р. Дністер здійснювалась за допомогою коефіцієнта кореляції [6]
, де , - дійсні та прогнозовані середньомісячні значення стоку, .Граничне значення коефіцієнта кореляції, коли
дорівнює одиниці. Для випадку, який розглядається (рисунок 3) , що свідчить про задовільний прогноз помісячного стоку р. Дністер.Рисунок 3. Кореляційна залежність між дійсним і прогнозованим значенням стоку р. Дністер