TT Tên sông Năm đo Chiều dài (m) Số mặt cắt
1 Hồng 2000 219.064 139 2 Trà Lý 2000 64.280 35 3 Ninh Cơ 2000 53.525 25 4 Đáy 2000 231.260 100 5 Thương 2000 91.235 33 6 Lục Nam 2000 53.080 16 7 Cầu 2000 136.393 48 8 Thái Bình 2000 95.210 55 9 Văn Úc 2000 40.300 23 10 Gùa 2000 3.310 6 11 Mía 2000 3.037 4 12 Mới 2000 3.580 2 13 Lạch Tray 2000 52.000 25 14 Kinh Thày 2000 47.620 28 15 Kinh Môn 2000 40.884 19 16 Đá Bạch 2000 38.500 13 17 Cấm 2000 24.623 14 18 Lai Vu 2000 27.619 14 19 Đuống 2000 60.615 33 20 Hoá 2000 38.157 20 21 Luộc 2000 70.683 35
2.2.4 Tài liệu thuỷ văn.
Trận lũ lớn xảy ra trong năm 1996 từ ngày 9 đến ngày 28 tháng 8 được chọn để hiệu chỉnh và trận lũ từ 15 - 30/7/2002 để kiểm định mơ hình thuỷ lực trong
nghiên cứu này. Tài liệu thực đo tại các trạm thủy văn trong lưu vực được cung cấp bởi Tổng cục Khí tượng Thuỷ văn (xem Bảng 27. ).
Quá trình lưu lượng ở các biên trên và các khu giữa được tính tốn bằng mơ hình mưa - dịng chảy.
Bảng 27. Các trạm thuỷ văn dùng để mô phỏng và kiểm định mô hình.
TT Tên Trạm đo Tên sơng Yếu tố đo Vị trí (m)
1 Hà Nội Hồng Q, H 21932
2 Hưng Yên -nt- H 140750
3 Thượng Cát Đuống Q, H 2300
4 Triều Dương Luộc H 5380
5 Nam Định Đào Nam Định H 5.133
6 Trực Phương Ninh Cơ H 4.239
7 Phả Lại Thái Bình H 300
8 Phủ Lý Đáy H 131.675
9 Thái Bình Trà Lý H 168.768
10 Gián Khẩu Hoàng Long H 60.700
2.2.5 Tính tốn mơ phỏng thủy lực hệ thống sơng
a. Tính tốn mơ phỏng trận lũ tháng 8/1996
Để xác định bộ thơng số cho mơ hình thuỷ lực như: (Mặt cắt đặc trưng từng đoạn sơng, hệ số nhám lịng và bãi sơng, các thơng số của cơng trình thuỷ lợi), cần thiết phải mô phỏng và hiệu chỉnh với một năm lũ thực tế. Qua phân tích tơi đã lựa chọn trận lũ tháng 8/1996 để mô phỏng chế độ thuỷ lực cho sông Hồng và sơng Thái Bình (đây là trận lũ đã được phân tích xem xét kỹ và được lấy làm cơ sở để tính tốn. Điều này phù hợp với tài liệu địa hình đo đạc mặt cắt sơng năm 2000, mặt cắt sông không biến đổi lớn. Lũ năm 1996 cũng là năm lũ khá lớn ở thượng du sông Hồng, sơng Thái Bình có dạng lũ bất lợi cho cơng tác phòng chống lũ hạ du, mặt khác trong thời gian có lũ thì phía các cửa sơng vùng ven biển Đông do ảnh hưởng của cơn bão số 4 cấp 11, 12 đổ bộ vào đất liền từ Thái Bình đến Nghệ An, gió to và sóng lớn gây nước dâng tới 1,5 ÷ 2 ,0 m. Phía cửa sơng thuộc hệ thống sơng Hồng, nước dâng do bão vào sâu nội địa tới 20 ÷ 25 km, phía cửa sơng thuộc hệ thống sơng Thái Bình, nước dâng do bão vào sâu nội địa tới 30 ÷ 40 km. Thời đoạn tính
tốn mơ phỏng kéo dài trong 10 ngày xảy ra lũ lớn nhất từ 16/8 đến hết ngày 25/8/1996.
b. Kết quả hiệu chỉnh
Sau khi q trình tính tốn hiệu chỉnh các thơng số của mơ hình thuỷ lực đã nhận được kết quả mực nước, lưu lượng tại các vị trí trạm thuỷ văn khá phù hợp với giá trị thực đo. Mực nước lũ lớn nhất tại các vị trí kiểm tra trên tồn mạng sơng như sau:
Bảng 28. Bảng đánh giá sai số hiệu chỉnh
TT Tên trạm Yếu tố Thực đo
(m) Tính tốn (m) Sai số(m) 1 Gián Khẩu H 3.450 3.425 0.025 2 Phủ Lý H 3.830 3.994 0.164 3 Nam Định H 4.810 4.98 0.17 4 Hà Nội H 12.430 12.451 0.021 5 Hà Nội Q 14800 14300 500 6 Hưng Yên H 7.860 7.92 0.06 7 Trực Phương H 3.140 3.041 0.099 8 Thái Bình H 3.750 3.71 0.04 9 Triều Dương H 6.700 6.71 0.01 10 Thượng Cát H 11.800 11.775 0.025 11 Thượng Cát Q 5740 5200 540 12 Phả Lại H 6.520 6.445 0.075
Hình 4. Kết quả hiê ̣u chỉnh tra ̣m Hà Nơ ̣i
Hình 6. Kết quả hiê ̣u chỉnh tra ̣m Trực Phương
Hình 8. Kết quả hiê ̣u chỉnh tra ̣m Triều Dương
Hình 10. Kết quả hiê ̣u chỉnh tra ̣m Phải La ̣i
Hình 12. Kết quả hiê ̣u chỉnh tra ̣m Gián Khẩu
Hình 13. Kết quả hiê ̣u chỉnh tra ̣m Nam Đi ̣nh
Trên cơ sở bộ thông số đã xác định được trong q trình tính tốn mơ phỏng trận lũ thực tế tháng 8/1996, tiến hành tính tốn kiểm định mơ hình với trận lũ thực tế đã xảy ra năm 2002. Trận lũ năm 2002 được chọn vì đây là trận lũ lớn xảy ra gần đây nhất. Trong trận lũ này, rất nhiều bối bãi dọc sông hồng bị ngập, nhiều đoạn đê xuất hiện mạch đùn, mạch sủi. Như vậy trận lũ năm 2002 đảm bảo độ lớn về lũ, gây ngập bối bãi.
Dưới đây là mực nước lũ lớn nhất và lưu lượng lũ lớn nhất giữa kết quả tính tốn và thực đo tại một số trạm đo thủy văn:
Bảng 29. Mực nước lũ lớn nhất thực đo và kết quả tính kiểm định tại các trạm thuỷ văn chính trên hệ thống sơng Hồng - Thái Bình.
TT Trạm Sông Hmax (m) Qmax (m
3
/s)
Thực
đo Tính tốn Sai số Thực đo Tính tốn Sai số
1 Hà Nội Hồng 12.01 11.92 0.09 13100 12797 303
2 Hưng Yên Hồng 7.983 14816.3
3 Thượng Cát Đuống 11.42 11.31 0.11 6200 6100.5 99.5
4 Bến Hồ Đuống 6.43 4963.5
5 Triều Dương Luộc 6.58 6.54 0.04 2680 2570 110 6 Quyết Chiến Trà Lý 4.96 5.01 -0.05 2080 2000 80 7 Thái Bình Trà Lý 3.44 3.56 -0.12
8 Nam Định Đào 4.48 4.61 -0.13 3920 4120 200
9 Trực Phương Ninh Cơ 3.65 1489.5
10 Ba Thá Đáy 2.88 3.12 -0.24 118.3 115.2 3.1
11 Phủ Lý Đáy 3.02 125.4
12 Ninh Bình Đáy 2.51 378.5
13 Phả Lại Thái Bình 5.8 6.01 -0.21 3880.3
Qua kết quả tính tốn của mơ hình có thể thấy: Sai số mực nước tại các trạm trong phạm vi dưới 30 cm. Tại các trạm chính như Sơn Tây, Hà Nội, Thượng Cát, Hưng Yên, Triều Dương, Quyết Chiến, Thái Bình, Nam Định đều nhỏ hơn 15 cm.
Sai số về lưu lượng không lớn, lớn nhất chiếm khoảng 5% so với lưu lượng chảy qua vị trí đó.
Các kết quả so sánh cho thấy bộ thông số mô phỏng cho trận lũ năm 1996 đã kiểm định phù hợp với trận lũ năm 2002. Bộ thông số nhám này đủ tin cậy để sử dụng tính tốn mơ phỏng lũ thiết kế và quy hoạch cho sơng Trà Lý.
2.3 TÍNH TỐN LŨ THIẾT KẾ TUYẾN SƠNG TRÀ LÝ
2.3.1 Tính tốn biên mạng thủy lực hệ thống sơng
a. Biên trên:
Với mạng sơng tính tốn đã được xác định ở trên, biên trên của mơ hình thuỷ lực là quá trình lưu lượng theo thời gian Q= f(t) tại các vị trí như sau:
+ Tại trạm thuỷ văn Cầu Sơn trên sông Thương + Tại trạm thuỷ văn Chũ trên sông Lục Nam + Tại trạm thuỷ văn Thái Nguyên trên sông Cầu
+ Tại trạm Sơn Tây: sử dụng đường quá trình lưu lượng tại Sơn Tây với Q= 28.000m3/s ứng với con lũ có chu kỳ xuất hiện là 250 năm (tần suất P= 4%). b. Biên dưới:
Biên dưới của mơ hình là q trình mực nước theo thời gian Z =f(t) tại 9 cửa sông đổ ra biển của lưu vực sơng Hồng – Thái Bình, là mực nước triều thiết kế 5% của tháng VIII. Kết quả tính tốn mực nước triều thiên văn lớn nhất theo các tần suất thiết kế đã được tính tốn ở phần tính tốn ở mục 1.1.2 ở phần trước như sau:
Bảng 30. Mực nước triều thiết kế tại các biên cửa sông
STT Cửa sông Hp 5% (cm) 1 Cấm 199 2 Đá Bạch 195 3 Lạch Tray 193 4 Văn Úc 191 5 Thái Bình 189 6 Trà Lý 185 7 Hồng (cửa Ba Lạt) 177 8 Ninh Cơ 170 9 Đáy 169
Nước dâng
Biên nước dâng ứng với tần suất 20% cũng đã được tính ở mục 1.1.2 với trị số như sau:
Bảng 31. Tần suất chiều cao nước dâng tại vùng đoạn bờ từ Cửa Ông – Cửa Đáy H nd(tb)
(cm)
Cv Cs Hndp(cm)
1% 2% 5% 10% 20% 50%
78 0.67 1.01 238 212 177 149 118 70
Và như vậy nước dâng ứng với tần suất 20% là H20%=118 cm
Tổng hợp kết quả tính tốn mực nước triều thiên văn lớn nhất tháng VIII tần suất 5% tại các cửa sông và nước dâng do bão tần suất 20% là 1,18m được trình bày tại
Bảng 32. Mực nước triều thiên văn 5% tháng VIII và chiều cao nước dâng P=20% Đơn vị : cm Cửa
sông Cấm Đá Bạch Lạch Tray Văn Úc Thái Bình Trà Lý Ba Lạt Ninh Cơ Đáy Htriều 5% 199 195 193 191 189 185 177 170 169 Hnd 20% 118 118 118 118 118 118 118 118 118 Tổng hợp 317 313 311 309 307 303 295 288 287
Các tuyến sơng có đê trên địa bàn tỉnh Thái Bình nằm ở vùng hạ du của hệ thống sông Hồng – sơng Thái Bình chịu sự chi phối về chế độ thủy lực của tồn mạng sơng, nhưng do nằm sát biển có tới 3 cửa sơng ra biển (trong số 9 cửa của toàn lưu vực) do vậy chịu tác động rất mạnh của chế độ mực nước biển (thủy triều, nước dâng do bão), ngoài ra cũng chịu ảnh hưởng trực tiếp tác động của nước biển dâng do biến đổi khí hậu. Các kịch bản tính tốn sẽ là tổ hợp của lũ thượng nguồn gặp các tổ hợp mực nước ở cửa sơng. Các tính tốn lũ đối với hệ thống sông đồng bằng Bắc Bộ thường chọn các dạng lũ lớn điển hình gây nhiều bất lợi đã sảy ra trong quá khứ, đó là lũ các năm 1969, 1971 và gần đây nhất là lũ năm 1996. Lũ các năm 1969 và 1971 thì đồng bằng Bắc bộ chủ yếu có hệ thống đê điều và các khu
phân chậm lũ để chống lũ. Trong khi đó lũ năm 1996 tuy là lũ lớn và bất lợi nhưng đã có các hồ Hịa Bình và Thác Bà chống lũ nên đã hoàn toàn chủ động chống lũ, giảm thiệt hại do lũ gây ra ở mức thấp nhất. Trong tính tốn ở đây khơng tính cả ba dạng lũ trên mà chỉ lựa chọn dạng lũ bất lợi năm 1996 để tính vì:
- Dạng lũ năm 1996 là dạng lũ bất lợi nhất và lượng lũ do nhánh sơng Đà đóng góp nhiều nhất
- Điều kiện địa hình hiện trạng của hệ thống gần với thời gian xảy ra lũ 1996 nhất nên tính tốn sẽ ít sai số.
Như vậy việc tính tốn lũ thiết kế cho sơng Trà lý sẽ được chọn là một trong tổ hợp của lũ thiết kế thượng nguồn gặp tổ hợp mực nước ở cửa sông như sau:
1. Lũ thiết kế cho lưu vực sông Hồng – Thái Bình dạng lũ 8/1996 gặp triều thực năm 1996.
2. Lũ thiết kế cho lưu vực sông Hồng – Thái Bình kết hợp với triều cường tần suất 5%.
3. Lũ thiết kế cho lưu vực sông Hồng – Thái Bình kết hợp với triều cường, nước dâng do bão.
Chi tiết các trường hợp tính hợp tốn như phần nội dung tính tốn dưới đây.
2.3.2 Nội dung các trƣờng hợp tính tốn lũ thiết kế.
a. Lũ thiết kế cho lưu vực sông Hồng - Thái Bình với triều thực năm 1996 (kí hiệu LTK1)
Theo Quyết định Phê duyệt quy hoạch phòng chống lũ hệ thống sơng Hồng, sơng Thái Bình (QĐ 92/2007/QĐ-TTg), được điều chỉnh mốc thời gian quy hoạch đến năm 2020.
Tiêu chuẩn về mực nước chống lũ cho hệ thống đê hạ du là 13,1 m tại trạm thủy văn Long Biên trên sông Hồng, tại trạm thủy văn Phả lại trên sơng Thái Bình là 7,2 m.
b. Lũ thiết kế cho lưu vực sơng Hồng - Thái Bình kết hợp với triều cường thiết kế (kí hiệu LTK2)
Trường hợp tính tốn này nhằm xác định mức độ gia tăng mực nước, lưu lượng đối với trường hợp triều cường 5% xảy ra.
Triều cường được tính từ triều thiên văn cho giai đoạn từ năm 1960 - 2009. Sau đó tính tốn mực nước triều thiên văn lớn nhất với tần suất thiết kế 5% cho thời kỳ max năm và max trong thời kỳ tháng VIII là tháng xuất hiện nhiều lũ lớn trên hệ thống sông. Biên triều cường được cho ở
Bảng .
c. Lũ thiết kế cho lưu vực sơng Hồng - Thái Bình với triều thiết kế kết hợp với nước dâng do bão thiết kế. (kí hiệu LTK3)
Đây là tổ hợp giữa lũ thượng nguồn gặp bão gây nước dâng ở phía cửa sơng và trong thời gian sảy ra triều cường. Điều kiện triều cường 5% kết hợp với nước dâng do bão với tần suất 20% được xem xét tính tốn trong trường hợp này. Biên mực nước triều cho mơ hình trong trường hợp tính tốn này trình bày ở
Bảng
2.3.3 Kết quả tính tốn thủy lực mạng sơng
Tiến hành tính tốn các tổ hợp lũ thượng nguồn và các kịch bản triều, nước dâng ở cửa sông cho tồn mạng sơng Hồng sơng Thái Bình, sau đó trích kết quả tính tốn mực nước và lưu lượng cho sông Trà Lý để phân tích lựa chọn lũ thiết kế cho tuyến sơng Trà Lý như phần dưới đây.
a. Kết quả lũ thiết kế dọc tuyến sơng Trà Lý
Kết quả tính tốn cực trị về mực nước và lưu lượng trên mạng sơng của các trường hợp tính tốn ở trên được như sau:
Bảng 33. Mực nước lũ lớn nhất của các trường hợp tính tốn thiết kế lũ thiết kế
Sơng Địa danh Vị trí (m) Phƣơng án
LTK1 LTK2 LTK3
Hồng Hà Nội 72392 12.76 12.76 12.75
Hồng Hưng Yên 140750 9.07 9.07 9.10
Đuống Thượng Cát 2300 11.93 11.94 11.94
Luộc Triều Dương 5380 7.46 7.47 7.53
Sông Địa danh Vị trí (m) Phƣơng án LTK1 LTK2 LTK3 - xã Hồng Minh, H. Hưng Hà 2000 6.46 6.49 6.69 - xã Bạch Đằng, H. Đông Hưng 4650 6.01 6.05 6.26 - xã Bạch Đằng, H. Đông Hưng 5750 6.00 6.03 6.22 - xã Bạch Đằng, H. Đông Hưng 6750 5.91 5.94 6.18 - TTV Quyết Chiến 8800 5.79 5.83 6.12
- xã Hồng Giang, H. Đông Hưng 10600 5.54 5.61 5.88 - xã Hoa Nam, H. Đông Hưng 11700 5.50 5.56 5.83 - xã Hoa Nam, H. Đông Hưng 13400 5.40 5.45 5.79 - xã Đông Dương, H. Đông Hưng 15500 5.25 5.34 5.69 - xã Đông Dương, H. Đông Hưng 18000 5.03 5.13 5.47 - xã Đông Dương, H. Đông Hưng 19300 4.95 5.04 5.36 - xã Đông Dương, H. Đông Hưng 21100 4.85 4.93 5.33 - xã Đơng Hịa, TX Thái Bình 23000 4.70 4.81 5.22 - xã Hồng Diệu, TX Thái Bình 24500 4.52 4.63 5.07 - xã Hồng Diệu, TX Thái Bình 26250 4.37 4.51 4.98
- TTV Thái Bình 27850 4.31 4.46 4.91
- xã Hồng Diệu, TX Thái Bình 29000 4.05 4.19 4.71 - xã Đông Mỹ, H. Đông Hưng 30600 3.92 4.05 4.57 - xã Đơng Hồng, H. Đơng Hưng 32300 3.81 3.98 4.53 - xã Đông Á, H. Đông Hưng 33500 3.71 3.91 4.49 - xã Đông Huy, H. Đông Hưng 35050 3.63 3.83 4.46 - xã Đông Lĩnh, H. Đông Hưng 36950 3.51 3.71 4.37 - xã Đông Lĩnh, H. Đông Hưng 38250 3.44 3.64 4.33 - xã Đông Lĩnh, H. Đông Hưng 39450 3.36 3.52 4.21 - xã Thái Hà, H. Thái Thụy 41450 3.16 3.48 4.17 - xã Thái Phúc, H. Thái Thụy 43000 2.91 3.25 4.11 - xã Thái Phúc, H. Thái Thụy 45200 2.85 3.19 4.08 - xã Thái Thành, H. Thái Thụy 47600 2.81 3.14 4.06 - xã Thái Thành, H. Thái Thụy 48400 2.76 3.11 4.05 - xã Thái Thọ, H. Thái Thụy 52200 2.55 3.01 3.91 - xã Thái Thọ, H. Thái Thụy 52800 2.54 3.00 3.90 - xã Thái Thọ, H. Thái Thụy 55700 2.11 2.57 3.77
- TTV Định Cư 56600 1.90 2.48 3.71
- xã Thái Đô, H. Thái Thụy 59300 1.85 2.42 3.67