4. Cấu trúc luận văn
3.3.1 Xây dựng phươngtrình dự báo mưa
Các tháng mùa lũ tại tỉnh Gia Lai được xác định như sau. Mùa lũ là những tháng có dòng chảy bình quân tháng lớn hơn hoặc bằng 8.3% dòng
chảy bình quân năm và có tần suất xuất hiện (P%) lớn hơn hoặc bằng 50% liên tục theo thời gian trong năm và còn lại là các tháng mùa cạn. theo tiêu chí này, mùa lũ trên các sông suối vùng phía Tây và trung tâm tỉnh bắt đầu từ tháng 7 kết thúc vào tháng 10. Mùa lũ trên các sông suối ở vùng phía Đông và
Đông Nam Tỉnh bắt đầu từ tháng 9 kết thúc vào tháng 12⌊6⌋.
Như vậy theo tính toán, trạm khí tượng Pleiku và Yaly các tháng mùa
mưa lũ là từ tháng 7 đến tháng 11. Trạm khí tượng An Khê và Ajun Pa các
tháng mùa mưa lũ là từ tháng 9 đến tháng 12. Tuy nhiên trong tháng 11 tại khu vực phía Tây và trung tâm tỉnh Gia Lai phổ biến thường không mưa, cá biệt có một số năm có mưamuộn gây lũ. Trong tháng 12 tại khu vực phía Đông vàĐông Nam tỉnh Gia Lai phổ biến là kết thúc mùa mưa vào thời kỳ đầu tháng, cá biệt có một số năm có mưa muộn xảy ra vào thời kỳ giữa tháng gây lũ. Do vậy tại luận văn này sẽ xây dựng phương trình hồi quy dự báo định lượng mưa cho trạm khí tượng Pleiku và Yaly ( đại diện cho khu vực phía tây và trung tâm tỉnh) trong các tháng 7, 8, 9, 10, trạm khí tượng An Khê và
AjunPa ( đại diện cho khu vực Đông, Đông Nam tỉnh) trong các tháng 9, 10,
11.
1, Tập các nhân tố dự báo
10 nhân tố dự báo ban đầu được lựa chọn là số liệu trung bình của các khu vực nói trên tại mực 850 mb và 700 mb, gồm các yếu tố là:
- Pv850: Giá trị xoáy thế trung bình tại mực 850 mb
-Z 850: Giá trị độ cao địa thế vịtrung bình tại mực 850 mb - U850: Giá trị gió vĩ hướngtrung bình tại mực 850 mb
- V850: Giá trị gió kinh hướngtrung bình tại mực 850 mb
- W850: Giá trị dòng thăngtrung bình tại mực 850 mb
- Pv700: Giá trị xoáy thế trung bìnhtrung bình tại mực 700 mb
-Z 700: Giá trị độ ao địa thế vịtrung bình tại mực 700 mb
- V700: Giá trị gió kinh hướngtrung bình tại mực 700 mb
- W700: Giá trị dòng thăngtrung bình tại mực 700 mb
Các số liệu được chiết suất trong thời gian 2012-2017 (06 năm, trong đó 05 năm dùng để tính toán xây dựng phương trình, năm 2017 dùng để kiểm nghiệmdự báo độc lập).
2, Kết quả xây dựng phương trình hồi quy dự báo mưa
Kết quả xây dựng phương trình hồi quy dự báo mưa 24h
Qua tính toán dùng các nhân tố dự báo là 10 yếu tố synop động lực chiết suất từ mực 850 mb và 700 mb của thời điểm 00Z của 24h trước để xây dựng phương trình dự báo mưa tại trạm lúc hiện tại, theo hồi quy từng bước,
kết quả tại các trạm như sau:
Trạm PleiKu: Tháng 7: R = 234,731 + 1,096 U850 - 0,156 Z850 – 41,135 W850 Tháng 8: R = 640,701 + 0,273 U700 - 0,201 Z700 + 45,792 W850 Tháng 9: R = 13,856 + 2,03 U700 - 1,314 U850 – 22,695 W700 Tháng 10: R = 0,125 + 33,322 Pv700 + 0,699 V850 – 15,483 W700 Trạm YaLy: Tháng 7: R = 836,258 - 0,264 Z700 + 21,116 Pv700 – 0,301 U850 Tháng 8: R = 906,011- 0,286Z700 - 35,898W850 + 0,371V700 Tháng 9: R = 426,287 + 0,946 U700 - 39,33W850 - 0,133Z700 Tháng 10: R = -4,176 + 41,53 Pv850-28,958W700+0,485V850 Trạm An Khê: Tháng 9: R = 14,457 - 35,257 W700 - 21,51 Pv700 - 43,155 W850 Tháng 10: R = 358,849-30,918W700 - 1,336U700 - 0,235Z850 Tháng 11: R = -11,565 - 4,336 U700 - 105,18 W700 - 1,192 V850 Trạm AjunPa: Tháng 9: R=6,84-0,448V700+1,01 U700-0,907U850 Tháng 10: R=0,502–22,328W850+16,71Pv700+0,239V850
Tháng 11: R=395,393–1,799U700–35,343W850–0,127Z700
Kết quả xây dựng phương trình hồi quy dự báo mưa 48h
Qua tính toán dùng các nhân tố dự báo là 10 yếu tố synop động lực chiết suất từ mực 850 mb và 700 mb của thời điểm 00Z của 48h trước để xây dựng phương trình dự báo mưa tại trạm lúc hiện tại, theo hồi quy từng bước,
kết quả tại các trạm như sau:
Trạm Pleiku: Tháng 7: R=468,77+0,816U850-0,306Z850-62,234Pv850 Tháng 8: R=26,355-124,55W700-136,37W850+1,934V700 Tháng 9: R=683,405+0,655U700-0,214Z700+38,254W850 Tháng 10: R=279,326-45,617W850-0,086Z700+0,354V850 Trạm Yaly: Tháng 7: R=585,939-0,388Z850+1,178V700-0,535V850 Tháng 8: R = 489,566 - 0,32Z850 - 38,138W700 - 0,545U850 Tháng 9: R = 504,016 - 0,328Z850 - 0,661U700 - 26,645Pv850 Tháng 10: R=-0,56-27,491W850-33,728W700+18,094Pv850 Trạm An Khê: Tháng 9: R=344,666-34,69Pv850-0,217Z850-0,498U850 Tháng 10: R = 998,373 - 0,314Z700 - 1,451U850 - 22,57Pv850 Tháng 11:R= -582,2 - 4,107U850 - 106,35W850 + 0,18Z700 Trạm Ajunpa: Tháng 9: R=8,209+0,266U700-17,889Pv700-24,856W700 Tháng 10: R = 513,99 - 0,163Z700 - 0,954U850 + 0,246V850 Tháng 11: R = -3,575 - 0,834U700 - 0,841U850 - 26,929W700