vực sông Dinh Ninh Hòa
* Lớp phủ thực vật
Vùng thượng lưu sông Dinh Ninh Hòa chủ yếu là rừng nguyên sinh lá rộng, xen kẽ là rừng hỗn giao tre nứa và trảng cỏ cây bụi; thảm phủ thực vật có sự phân hóa theo độ cao địa hình [16]. Theo niên giám thống kê tỉnh Khánh Hòa năm 2019, diện tích rừng của thị xã Ninh Hòa là 52147 ha, chỉ số phát triển rừng là 105,32%; rừng trồng mới tập trung là 400 ha, chỉ số phát triển trồng mới là 317,46%; rừng được giao khoán bảo vệ là 701 ha, chỉ số phát triển là 70,1%.
1.6.4. Đặc điểm lưu vực sông Cái Phan Rang
a) Vị trí địa lý
Sông bắt nguồn từ đỉnh núi vùng núi cao thuộc tỉnh Khánh Hòa và chảy theo hướng bắc - nam, khi cách cửa biển 35 km, sông đổi sang hướng tây bắc - đông nam và cuối cùng đổ ra cửa Đông Hải. Sông Cái Phan Rang có 13 phụ lưu bên hữu ngạn và 4 phụ lưu bên tả ngạn, có tổng diện tích lưu vực 3000 km2, chiều dài sông 119 km, chiều dài lưu vực 95 km, độ rộng bình quân lưu vực 31,6 km, độ dốc bình quân lưu vực 17,70/00. Hệ thống sông suối phân bố theo dạng cành cây gồm các phụ lưu hữu ngạn là Ty Cây, Ta Gu, CoRot, Hầm Leo, Đa May, Gia Nhong, Đa Cát Rum, Địa Gốc, sông Ông, sông Cha, sông Dâu,
sông Lanh Ra và sông Lu. Các phụ lưu bên tả ngạn gồm: sông Địa Gan, sông Sắt, Cho Mo, suối Sa Ra (Hình 1.3) [34].
b) Địa hình
Địa hình núi cao bao bọc gần như toàn bộ lưu vực ở các hướng: phía bắc là Núi Đào (1451 m), núi Chúa (1040 m); phía tây là phía nam dãy Trường Sơn với hàng loạt đỉnh cao hơn 1500 m, có đỉnh Bi Đúp cao nhất 2280 m; phía nam là các núi thấp hơn như Đa Khum (898 m), Cà Ná (644 m). Vùng đồng bằng trũng hạ du có cao độ trên dưới 10 m. Độ cao bình quân lưu vực 483 m. Địa hình đồi thấp xen kẽ đồng bằng phân bố ở vùng trung lưu sông Cái. Từ Tân Sơn đến Tân Mỹ, lưu vực sông mở rộng, độ dốc lòng sông còn cao, lòng sông nhiều đá tảng, một số nơi có các bãi bồi giữa sông như một sự pha trộn giữa kiểu sông miền núi và đồng bằng. Từ Tân Mỹ đến của Đông Hải, sông chảy êm trong một vùng đồi thấp và đồng bằng Phan Rang nhỏ hẹp (Hình 1.3) [34].
c) Đặc điểm địa chất, thổ nhưỡng, lớp phủ thực vật [34]
* Đặc điểm địa chất: Địa chất của tỉnh thuộc nền địa chất granit, phức
hệ Đèo Cả - Đơn Dương, giàu khoáng sản phi kim loại, nhất là nguyên vật liệu xây dựng.
* Đặc điểm thổ nhưỡng:
Thổ nhưỡng là tổ hợp 24 loại đất với phần lớn là đất đỏ, mặn, phèn mặn và nâu vàng (đất núi Feralit), có chiều sâu phong hóa dày. Trong đó, chủ yếu là đất đỏ vàng chiếm tỷ lệ lớn và phân bố ở khu vực Phan Rang, Nha Hố, Bác Ái, Ninh Sơn, Thuận Bắc; đất vàng ở Phan Rang, Ninh Phước, Thuận Nam; đất cát thành phần cơ giới nhẹ và thô, kết cấu rời rạc, phân
bố phần lớn vùng ven biển. Hình 1.3. Minh họa bản đồ lưu vực Cái Phan Rang
* Lớp phủ thực vật: Tỷ lệ che phủ rừng chiếm 67,71% tổng diện tích đất
tự nhiên, diện tích đất nông nghiệp chiếm 29,98% đất tự nhiên, chủ yếu là lúa và cây công nghiệp ngắn ngày theo thời vụ. Thảm phủ thực vật trên lưu vực nhìn chung nghèo nàn, diện tích rừng tự nhiên rất ít và thưa, chủ yếu là cỏ xen lùm bụi, cỏ lùm bụi xen cây gỗ rải rác.
1.7. NHỮNG VẤN ĐỀ NGHIÊN CỨU CỦA LUẬN ÁN
Các lưu vực sông Cái Nha Trang, sông Dinh Ninh Hòa và sông Cái Phan Rang tuy là các lưu vực có diện tích không lớn nhưng đóng vai trò cực kỳ quan trọng với sự phát triển kinh tế xã hội của khu vực Nam Trung Bộ nói chung và các tỉnh Khánh Hòa và Ninh Thuận nói riêng. Hơn nữa, các khu đô thị, trung tâm kinh tế và chính trị của các tỉnh (thành phố Nha Trang, thành phố Phan Rang, thị xã Ninh Hòa) đều nằm ở hạ lưu của các lưu vực sông này. Do vậy nhu cầu về mô phỏng và dự báo thủy văn chi tiết theo không gian là rất cấp thiết; tuy nhiên, khu vực Nam Trung Bộ chưa được áp dụng mô hình thủy văn thông số phân bố trong nghiệp vụ dự báo và ứng dụng trong mô phỏng chưa được phổ biến; do đó, vấn đề nghiên cứu của Luận án góp phần đáp ứng yêu cầu cấp thiết và hạn chế những khó khăn trên. Ngoài ra, các lưu vực sông ở khu vực Nam Trung Bộ có nhiều hồ chứa không có nhiệm vụ điều tiết, thủy triều không tác động vào sâu trong sông, khu vực thượng lưu thiếu trạm đo mưa và mặt cắt ngang sông; do đó, ứng dụng mô hình MARINE cải tiến trên các lưu vực sông ở Nam Trung Bộ là khá phù hợp.
Trong số nhiều mô hình thủy văn thông số phân bố được ứng dụng có hiệu quả ở nhiều nơi trên Thế giới và Việt Nam, mô hình MARINE được lựa chọn do có cơ sở vật lý tường minh, mã nguồn mở dễ áp dụng, đang được sử dụng tại Trung tâm Dự báo Khí tượng Thủy văn quốc gia và cho phép can thiệp, cải tiến để phù hợp với các điều kiện thực tế tại địa phương. Trên cơ sở phân tích các đặc điểm về thiếu số liệu của các lưu vực sông nghiên cứu, các hạn chế của mô hình MARINE nguyên gốc, cũng như các định hướng khắc phục để mô phỏng, dự báo dòng chảy trên lưu vực sông Cái Nha Trang, Dinh Ninh Hòa và Cái Phan Rang, Luận án đã lựa chọn MARINE và đặt ra mục tiêu cải tiến mô hình ở các khía cạnh sau:
- Cải tiến mô phỏng dòng chảy trong sông của mô hình MARINE bằng phương pháp tích hợp mô hình sóng động học một chiều và tăng khả năng ứng dụng trên lưu vực thiếu số liệu mặt cắt ngang.
- Cải tiến khả năng mô phỏng tác động của hồ chứa trong mô hình MARINE bằng phương pháp tích hợp mô đun diễn toán dòng chảy qua hồ để có thể ứng dụng cho các lưu vực đã có nhiều tác động của các hoạt động (công trình) kinh tế xã hội.
- Cải tiến tính toán mưa phân bố không gian bằng phương pháp tích hợp công cụ nội suy mưa và tăng khả năng khắc phục thiếu số liệu mưa cho các lưu vực có mạng lưới quan trắc mưa thưa thớt.
Hình 1.4. Sơ đồ nghiên cứu của Luận án
Các mô hình, mô đun, công cụ trên được kiểm tra trước và sau khi tích hợp với mô hình MARINE, đánh giá hiệu quả của quá trình cải tiến trong việc mô phỏng dòng chảy lũ trên các lưu vực sông Cái Nha Trang, Dinh Ninh Hòa
Chiều rộng và độ dốc sông Mô hình MARINE gốc Công cụ nội suy mưa Mô hình sóng
động học Mô đun diễn toán hồ chứa
Cải tiến mô hình MARINE
Mô phỏng thử nghiệm Dự báo thử nghiệm
Đánh giá
Đánh giá Đánh giá
V~F~Z Mưa trạm
DEM Thảm phủ Loại đất Mực nước ngầm
và Cái Phan Rang. Mô hình MARINE cải tiến sau đó sẽ tiếp tục được ứng dụng thử nghiệm trong nghiệp vụ dự báo hàng ngày tại Đài Khí tượng Thủy văn khu vực Nam Trung Bộ trên lưu vực sông Cái Nha Trang nhằm đánh giá tổng quát hơn cả về khả năng mô phỏng và khả năng dự báo dòng chảy lũ của mô hình.
TIỂU KẾT CHƯƠNG 1
- Yêu cầu mô phỏng và dự báo chi tiết theo không gian dẫn đến xu hướng sử dụng mô hình thủy văn thông số phân bố; tuy nhiên, ứng dụng mô hình gặp nhiều khó khăn do thiếu trạm đo và số liệu phân bố. Hiện nay có nhiều mô hình thủy văn thông số phân bố được sử dụng và mô hình MARINE là lựa chọn phù hợp do đã chứng minh được khả năng ứng dụng ở cả Việt Nam và trên Thế giới, có cơ sở khoa học tường minh và có mã nguồn mở.
- Mặc dù mô hình MARINE có khả năng mô phỏng tốt dòng chảy sườn dốc và đáp ứng yêu cầu mô phỏng, dự báo chi tiết; tuy nhiên, mô hình còn một số nhược điểm cần cải tiến và tăng khả năng mô phỏng trên lưu vực thiếu số liệu, chưa có các công cụ xử lý mưa phân bố không gian từ số liệu trạm và chưa xét đến các công trình như hồ chứa, đập dâng.
- Tích hợp mô hình sóng động học một chiều để hoàn thiện mô phỏng dòng chảy trong sông của mô hình MARINE và tăng khả năng ứng dụng trên lưu vực sông thiếu số liệu mặt cắt. Bổ sung mô phỏng tác động của hồ chứa trong mô hình MARINE bằng mô đun diễn toán dòng chảy qua hồ. Công cụ nội suy mưa được tích hợp với mô hình MARINE để khắc phục nhược điểm tính toán mưa phân bố không gian và thiếu trạm đo mưa thượng nguồn.
- Các lưu vực sông nhỏ, thiếu số liệu quan trắc ở khu vực Nam Trung Bộ (lưu vực sông Cái Nha Trang, Dinh Ninh Hòa và Cái Phan Rang) có tầm quan trọng đặc biệt với các tỉnh Khánh Hòa, Ninh Thuận và có nhu cầu cấp thiết về mô phỏng và dự báo lũ chi tiết theo không gian. Vì thế cần ứng dụng mô hình MARINE cải tiến để mô phỏng và dựa báo lũ trên các lưu vực sông không ảnh hưởng triều và thiếu số liệu mặt cắt, mưa phân bố như các lưu vực sông ở trên và qua đó chứng minh hiệu quả cải tiến của bộ công cụ mô hình.
CHƯƠNG 2: NGHIÊN CỨU CẢI TIẾN MÔ HÌNH MARINE 2.1. CƠ SỞ LÝ THUYẾT MÔ HÌNH MARINE 2.1. CƠ SỞ LÝ THUYẾT MÔ HÌNH MARINE
2.1.1. Lý thuyết dòng chảy sườn dốc của mô hình MARINE
Mô hình MARINE mô phỏng dòng chảy trên sườn dốc dựa trên phương trình xấp xỉ sóng động học và tính thấm bằng phương pháp Green - Ampt trên các ô lưới (pixel vuông). Mỗi ô lưới có khả năng trao đổi nước qua 4 ô lưới xung quanh, hướng dòng chảy là hướng có độ dốc lớn nhất; các ô lưới được liên kết lại với nhau để tính toán trao đổi nước và tập trung dòng chảy trên các ô sườn dốc đến các ô sông. Các số liệu đầu vào cần thiết cho mô hình MARINE để tính toán dòng chảy sườn dốc từ mưa gồm có:
- Bản đồ mô hình số độ cao (DEM): được sử dụng để tính độ dốc. - Bản đồ thành phần đất (thông số địa chất và thổ nhưỡng) và độ cao mực nước ngầm (tính toán độ sâu và khả năng thấm tối đa của đất).
- Bản đồ sử dụng đất và thảm phủ thực vật : được sử dụng để xác định hệ số nhám, tính tốc độ dòng chảy mặt.
- Mưa phân bố không gian các thời đoạn. - Bản đồ độ ẩm của đất trên lưu vực sông.
Phương trình mô phỏng dòng chảy trên sườn dốc trong mô hình MARINE xuất phát từ hệ phương trình Saint - Venant. Mô hình MARINE sử dụng hệ phương trình sau [116]:
- Phương trình liên tục ở trạng thái không ổn định được biểu thị cục bộ:
𝜕ℎ
𝜕𝑡 +𝜕𝑢ℎ
𝜕𝑥 = 𝑃 − 𝐼 (2.1.1)
Trong đó: h là độ cao của mặt nước (m), t là thời gian (s), u là tốc độ dòng chảy mặt (m/s), x là biến không gian (m), P thể hiện cường độ của mưa (m/s) và I là tốc độ thấm (m/s). - Phương trình động lượng: 1 𝐴 𝜕𝑄 𝜕𝑡 +1 𝐴 𝜕 𝜕𝑥(𝑄2 𝐴) + 𝑔𝜕𝑦 𝜕𝑥 − 𝑔(𝑆𝑜 − 𝑆𝑓) = 0 (2.1.2) Trong phương trình động lượng chỉ giữ lại thành phần độ dốc của đường năng lượng và độ dốc cục bộ của đất được phương trình động lượng cho sóng động học như sau:
S0 = Sf (2.1.3) Phương trình động lượng của sóng động học được rút gọn thành chế độ đồng nhất, thu được phương trình xấp xỉ sóng động học:
𝑢 = 𝑎. ℎ𝑚−1 (2.1.4)
Trong đó: a và m là các hệ số.
Khi dòng chảy này diễn ra dọc theo các rãnh, ở mức nước lớn hơn nó trở nên hỗn loạn hoàn toàn. Tốc độ có thể được biểu thị bằng công thức của Manning (Carlier, 1986) bằng cách lấy h là giá trị của bán kính thủy lực:
𝑢 = 𝑆𝑜1/2
𝑛 . ℎ2/3 (2.1.5)
Trong đó: n là hệ số Manning và S0 là độ dốc (m/m).
Thay phương trình động lượng vào phương trình liên tục (2.1.1), thu được phương trình như sau:
𝜕ℎ 𝜕𝑡 +𝑆𝑜1/2 𝑛 5 3 𝜕ℎ 𝜕𝑥 = 𝑃 − 𝐼 (2.1.6)
Đối với mỗi phần tử của lưới có cấu trúc MARINE, bốn hướng tiềm năng của dòng chảy được xem xét (các hướng chính) và hướng duy nhất được sử dụng là hướng có độ dốc lớn nhất. Đường trượt nước được giả định là phân bố đồng đều trên toàn bộ bề mặt của pixel vuông được xem xét. Các phương trình bảo toàn được giải bằng phương pháp thể tích hữu hạn. Phương trình liên tục (2.1.1) có thể được viết trên bề mặt điều khiển Ωk được giới hạn bởi Γk như sau [116]:
∫Ω 𝜕ℎ𝜕𝑡
𝑘 𝑑Ω𝑘 + ∫ 𝑑𝑖𝑣(𝑢ℎ)Ω
𝑘 𝑑Ω𝑘 = ∫ (𝑃 − 𝐼)Ω
𝑘 𝑑Ω𝑘 (2.1.7) Biến đổi phương trình (2.1.7) theo định lý Ostrogradsky, số hạng thứ hai biến thành tích phân đường bao Γk:
∫Ω 𝜕ℎ𝜕𝑡𝑘
𝑘 𝑑Ω𝑘 + ∫ 𝑢Γ 𝑘𝑑𝑘ℎ𝑘
𝑘 𝑑Γ𝑘 = ∫ (𝑃𝑘 − 𝐼𝑘)
Ω𝑘 𝑑Ω𝑘 (2.1.8)
Trong đó: dk là hướng của độ dốc lớn nhất và biến số hk không đổi trên Ωk. Áp dụng sơ đồ sai phân hiện Euler theo không gian và thời gian:
ℎ𝑘𝑡+Δ𝑡−ℎ𝑘𝑡
Δ𝑡 +𝑢𝑘𝑡ℎ𝑘𝑡
Δ𝑥 = 𝑃𝑘𝑡 − 𝐼𝑘𝑡 (2.1.9)
Trong đó: t đại diện cho thời gian và ∆t là bước thời gian. Phương trình xấp xỉ sóng động học được viết như sau:
𝑢𝑘𝑡 =𝑆𝑘 1/2
Cuối cùng, chiều cao lớp nước trên mỗi pixel được tính theo phương trình sau: ℎ𝑘𝑡+Δ𝑡 = ℎ𝑘𝑡 + Δ𝑡(𝑃𝑘𝑡− 𝐼𝑘𝑡) − Δ𝑡 Δ𝑥 𝑆𝑘1/2 𝑛𝑘 (ℎ𝑘𝑡)5/3 (2.1.11) Sơ đồ trên được thực hiện khá dễ dàng và tính toán nhanh; tuy nhiên, bước thời gian tính toán phải thấp để đảm bảo tính ổn định của dòng chảy. Có thể ước tính bước thời gian này nhờ vào điều kiện ổn định của Courant, Friedrichs, Lewy (CFL).
Hình 2.1. Sơ đồ tổng quát của mô hình MARINE [3]
2.1.2. Phương pháp tính thấm trong mô hình MARINE
Phương trình của Green và Ampt (1911) dựa trên phương trình của Darcy (1856) liên hệ tốc độ của dòng chảy I (m/s) với gradient thủy lực dh/dz thông qua độ dẫn thủy lực K (m/s), có dạng như sau [116]:
𝐼 = 𝐾𝑑ℎ
𝑑𝑧 (2.1.12)
Định luật Darcy được sử dụng cho các điều kiện dòng chảy tầng trong môi trường đồng nhất, đẳng hướng và liên tục, chất lỏng không tương tác với môi trường. Độ dẫn thủy lực K là đặc trưng của điều kiện dòng chảy trong một môi trường xốp nhất định đối với một chất lỏng nhất định.
Bản đồ DEM Bản đồ Thảm phủ Bản đồ Thổ nhưỡng Bản đồ H ngầm
Bản đồ Mưa (Đa giác) Bản đồ Độ dốc Bản đồ Hệ số nhám Bản đồ Thấm
Lưu lượng ô lưới Lưu lượng cộng
dồn các đoạn 1 sông chính
Sơ đồ tính thấm được Green và Ampt (1911) đề xuất bằng việc tích hợp mặt trước thấm ướt theo phương ngang và hướng dốc, ngăn cách đất không bão hòa với đất bão hòa. Nước có áp suất thấp hơn áp suất khí quyển, bị mao dẫn Sf hút vào môi
trường xốp (Hình 2.2) [116]. Hình 2.2. Sơ đồ thấm Green Ampt [116] Theo sơ đồ này, phương trình (2.1.12) trở thành:
𝐼 = 𝐾ℎ𝑓−ℎ0
𝑍𝑓−𝑍0 = 𝐾𝑆𝑓+𝑍𝑓
𝑍𝑓 (2.1.13)
Trong đó: Hf là cột nước từ ranh giới ẩm và h0 là cột nước bề mặt (h0 = 0, có nghĩa là không có nước bề mặt), Sf là cột nước mao dẫn.
Độ sâu đến ranh giới ẩm có liên quan đến độ sâu tích lũy của nước thấm (F), được thể hiện bằng quan hệ như sau:
F = Zt(θs - θi) (2.1.14) Trong đó: F là độ sâu tích lũy của nước thấm (m), θs là độ ẩm hoặc độ xốp bão hòa của đất (%) và θi là độ ẩm ban đầu của đất (%).
Bằng cách kết hợp các phương trình (2.1.13) và (2.1.14), Mein và Larson (1973) đề xuất công thức tính tốc độ thấm I(t) như sau:
𝐼(𝑡) = 𝑑𝐹(𝑡) 𝑑𝑡 = {
𝑃 𝑘ℎ𝑖 𝑡 ≤ 𝑡𝑝 𝐾. (1 + 𝑆𝑓.𝜃𝑠−𝜃𝑖
𝐹(𝑡) ) 𝑘ℎ𝑖 𝑡 > 𝑡𝑝 (2.1.15)
Trong đó: P là cường độ của mưa (m/s), tp là thời gian điền trũng (để hình thành vũng nước trên bề mặt đất).
Tốc độ thấm I trên thực tế bằng cường độ mưa P nhưng về sau nhỏ hơn tốc độ thấm tiềm năng I(t). Khi cường độ mưa lớn hơn tốc độ thấm, đất trở nên bão hòa và tại một thời điểm (t) hình thành lớp nước không đáng kể (Gourley và Vieux, 2006). Phương trình (2.1.15) không thể hiện rõ ràng một biến số thời
Lớp nước Lớp đất bão hòa ẩm