Tổng quan nghiên cứu

Lưu vực sông Trà Lý có chiều dài 64 km, giữ vai trò huyết mạch trong việc tiêu thoát lũ từ dòng chính sông Hồng ra vịnh Bắc Bộ và bảo vệ an toàn cho toàn bộ tỉnh Thái Bình với quy mô dân số hơn 1.780.954 người (thống kê năm 2009). Thực tế quan trắc thủy văn cho thấy sông Trà Lý gánh vác việc phân chuyển từ 8% đến 11% tổng lượng lũ từ trạm Sơn Tây ra biển. Điển hình trong các trận lũ lịch sử như năm 1971 và 1996, lưu lượng đỉnh lũ đổ qua trạm Quyết Chiến đã vượt mức 2.180 m³/s, gây áp lực to lớn lên hệ thống đê điều.

Tuy nhiên, sức ép từ phát triển kinh tế xã hội giai đoạn 2006–2010 với tốc độ tăng trưởng GDP đạt 12,04%/năm đã dẫn đến tình trạng khai thác bãi bồi ven sông thiếu kiểm soát. Các tuyến đê bối tự phát như bối Trà Giang, bối Hồng Thái liên tục bị tôn cao và lấn chiếm lòng dẫn, làm co hẹp nghiêm trọng hành lang thoát lũ tự nhiên. Bên cạnh đó, tác động kép từ biến đổi khí hậu toàn cầu và sự xuất hiện của các hồ chứa thủy điện lớn ở thượng nguồn (hồ Hòa Bình, Tuyên Quang, Sơn La) đã làm thay đổi căn bản chế độ thủy động lực học vùng hạ du chịu ảnh hưởng triều.

Trước tình hình đó, nghiên cứu quy hoạch phòng chống lũ chi tiết tuyến sông Trà Lý được thực hiện nhằm xác định chuẩn xác các thông số lũ thiết kế, thiết lập mực nước lũ tương ứng với các cấp tần suất (P = 0,4% đến P = 1%), và đề xuất phương án quản lý sử dụng các vùng bối bãi bền vững. Kết quả nghiên cứu cung cấp cơ sở khoa học và pháp lý trực tiếp cho việc điều chỉnh quy hoạch sử dụng đất, bảo vệ an toàn tính mạng và tài sản cho 133.693 ha diện tích tự nhiên toàn tỉnh Thái Bình.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng thủy lực học dòng chảy một chiều không ổn định trong lòng dẫn hở, biểu diễn qua hệ phương trình Saint-Venant vi phân phi tuyến gồm phương trình liên tục và phương trình động lượng. Hệ phương trình này cho phép mô phỏng chính xác sự truyền sóng lũ, tương tác sóng triều và sự biến đổi mực nước dọc theo các phân lưu hạ du.

Mô hình thủy lực MIKE 11 (do Viện Thủy động lực Đan Mạch - DHI phát triển) được ứng dụng làm công cụ tính toán mô phỏng cốt lõi. Trong đó, các khái niệm chuyên ngành then chốt được chuẩn hóa bao gồm:

  1. Lũ thiết kế: Lưu lượng và mực nước đỉnh lũ ứng với chu kỳ lặp lại xác định (tần suất thiết kế P = 1% tương ứng chu kỳ 100 năm hoặc P = 0,4% theo Quyết định 92/2007/QĐ-TTg).
  2. Hành lang thoát lũ và chỉ giới thoát lũ: Khoảng không gian lòng sông và bãi sông tối thiểu cần bảo vệ để đảm bảo chuyển tải an toàn lưu lượng lũ thiết kế mà không làm dâng mực nước đê chính vượt quá cao trình thiết kế.
  3. Đê bối và bãi bồi ven sông: Các đê bao khoanh vùng bãi trũng ngoài đê chính nhằm bảo vệ sản xuất nông nghiệp trong điều kiện lũ vừa và nhỏ, nhưng phải đảm bảo điều kiện tháo tràn hoặc ngập có kiểm soát khi lũ lớn xuất hiện.

Phương pháp nghiên cứu

Phương pháp nghiên cứu tích hợp phương pháp phân tích thống kê chuỗi thời gian thủy văn kết hợp với mô hình hóa toán học số trị thủy động lực. Nguồn dữ liệu đầu vào bao gồm:

  • Chuỗi quan trắc khí tượng thủy văn nhiều năm (1950–2008) tại các trạm thủy văn chủ lực: Quyết Chiến, Thái Bình, Định Cư (sông Trà Lý), Ba Lạt (sông Hồng), Triều Dương (sông Luộc), Thượng Cát (sông Đuống) và trạm hải văn Hòn Dấu.
  • Dữ liệu địa hình mặt cắt ngang lòng dẫn sông Hồng - Thái Bình và chi tiết tuyến sông Trà Lý từ cửa vào Hồng Minh (K0) đến cửa ra vịnh Bắc Bộ (K64) với hàng trăm mặt cắt thực đo chi tiết.

Quy trình mô phỏng thủy lực trên MIKE 11 được thiết lập qua hai giai đoạn chặt chẽ:

  • Giai đoạn hiệu chỉnh mô hình: Sử dụng số liệu thực đo của các trận lũ lịch sử (năm 1969, 1971, 1996) để điều chỉnh hệ số nhám Manning và các thông số trở lực lòng dẫn.
  • Giai đoạn kiểm định mô hình: Đối chiếu mực nước tính toán với mực nước đo thực tế tại 10 trạm thủy văn trọng điểm (Hà Nội, Hưng Yên, Trực Phương, Thái Bình, Phả Lại...), đảm bảo sai số mực nước đỉnh lũ không vượt quá 5–10 cm.

Lý do lựa chọn mô hình một chiều MIKE 11 bắt nguồn từ tính phù hợp tối ưu với mạng lưới sông phân lưu phức tạp, khả năng tính toán kết nối đa biên thủy văn thượng lưu và biên triều hạ lưu biến động liên tục theo thời gian thực.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng thủy lực và phân tích chuỗi số liệu thực tế đã mang lại nhiều phát hiện quan trọng:

Thứ nhất, vai trò chia lũ của sông Trà Lý đối với dòng chính sông Hồng đạt tỷ lệ phân lưu trung bình từ 8,6% đến 11% vào mùa lũ. Trong điều kiện lũ cực hạn, lưu lượng đổ vào sông Trà Lý có thể đạt từ 2.180 m³/s đến hơn 2.450 m³/s, trực tiếp làm giảm cao trình mực nước nguy hiểm tại trạm thủy văn Nam Định và cửa Ba Lạt.

Thứ hai, nghiên cứu đã xác định mực nước lũ thiết kế ứng với tần suất P = 1% tại các trạm then chốt dọc tuyến sông:

  • Trạm Quyết Chiến (thượng lưu sông Trà Lý): Đạt cao trình mực nước thiết kế 654 cm (vượt đỉnh lũ lịch sử năm 1971 là 645 cm khoảng 9 cm).
  • Trạm Thái Bình (trung lưu): Đạt cao trình mực nước thiết kế 444 cm (vượt mức lũ lịch sử năm 1971 là 421 cm khoảng 23 cm).
  • Trạm Định Cư (hạ lưu ven biển): Đạt cao trình mực nước thiết kế 341 cm (vượt đỉnh lũ năm 1968 là 319 cm khoảng 22 cm).

Thứ ba, sự thay đổi thủy văn sau khi hồ Hòa Bình và Tuyên Quang đi vào vận hành cho thấy tính chất phân hóa sâu sắc. Tại trạm Quyết Chiến, mực nước đỉnh lũ trung bình thời kỳ 1988–2008 giảm 7 cm so với thời kỳ 1960–1987. Tuy nhiên, ở hạ lưu (trạm Thái Bình và Định Cư), mực nước đỉnh lũ lại có xu hướng tăng từ 1 cm đến 2 cm do sự kết hợp của mực nước triều dâng do bão và hiện tượng bồi lắng lòng dẫn. Tại cửa Trà Lý, mực nước triều thiết kế P = 5% vào tháng VIII kết hợp nước dâng do bão tần suất P = 20% (chiều cao nước dâng 1,18 m) được xác định lên tới 3,03 m.

Thứ tư, việc đánh giá các vùng đê bối (bối Tịnh Thủy, bối Thái Thọ, bối Trà Giang, bối Hồng Thái) chỉ ra rằng nếu tất cả các bối bãi đều tự ý tôn cao để chống lũ triệt để, mặt cắt thoát lũ bị co hẹp cục bộ làm dâng mực nước lũ thiết kế dọc tuyến đê chính từ 15 cm đến 35 cm, tạo ra nguy cơ vỡ đê hàng loạt trên diện rộng.

Thảo luận kết quả

Diễn biến mực nước dọc tuyến sông Trà Lý là kết quả tương tác phức tạp giữa áp lực lũ thượng nguồn đổ về và lực cản của sóng triều vịnh Bắc Bộ dâng lên. Dữ liệu tính toán từ mô hình cho phép thiết lập bản đồ đường mặt nước dọc sông và mặt cắt thủy lực co hẹp, minh chứng rõ nét mối quan hệ nghịch đảo giữa diện tích bối bãi bị lấn chiếm và độ an toàn của đê chính.

Khi so sánh với các nghiên cứu trước năm 1988, hệ thống hồ chứa thượng nguồn đã cắt giảm đáng kể đỉnh lũ ở thượng lưu, nhưng không thể triệt tiêu rủi ro ngập lụt ở vùng hạ lưu ven biển tỉnh Thái Bình. Tác động của nước dâng do bão (thống kê cho thấy 31% số trận bão đổ bộ vào vùng duyên hải từ Quảng Ninh đến Ninh Bình) kết hợp kỳ triều cường nhật triều có biên độ dao động lên tới 3,94 m khiến khu vực hạ lưu luôn ở mức báo động cao. Do đó, việc quy hoạch hành lang thoát lũ và phân loại các đê bối (bối giữ lại bảo vệ dân cư, bối cho phép tràn lũ chậm, bối cần giải tỏa lòng sông) là giải pháp bắt buộc để cân bằng giữa bảo đảm an toàn công trình đê điều và duy trì không gian phát triển kinh tế bãi sông.

Đề xuất và khuyến nghị

Dựa trên kết quả tính toán thủy lực, nghiên cứu đề xuất 4 nhóm giải pháp quy hoạch công trình và phi công trình với mục tiêu và lộ trình cụ thể:

  1. Cắm mốc chỉ giới hành lang thoát lũ và thanh thải chướng ngại vật lòng sông:

    • Hành động: Tiến hành khảo sát, cắm mốc chỉ giới thoát lũ chi tiết trên toàn bộ 64 km tuyến sông Trà Lý; cưỡng chế tháo dỡ các công trình xây dựng trái phép lấn chiếm bãi bồi ngoài đê chính.
    • Mục tiêu: Đảm bảo bề rộng thoát lũ tối thiểu đạt từ 150 m đến 255 m tại các đoạn sông cong trọng yếu, giảm thiểu 100% tình trạng co hẹp dòng chảy nhân tạo.
    • Chủ thể và thời gian: Chi cục Quản lý Đê điều và Phòng chống lụt bão Thái Bình phối hợp cùng UBND các huyện (Hưng Hà, Đông Hưng, Vũ Thư, Kiến Xương, Thái Thụy) hoàn thành trong giai đoạn 2013–2015.

  2. Quy hoạch và phân loại hệ thống đê bối theo từng mức chống lũ:

    • Hành động: Phân chia rõ ràng các vùng bối nông nghiệp và bối dân sinh; thiết kế các tuyến tràn cứu hộ có kiểm soát tại bối Trà Giang, bối Hồng Thái khi mực nước sông vượt mức báo động cấp 3.
    • Mục tiêu: Duy trì an toàn tuyệt đối cho tuyến đê chính cấp II và cấp I khi xuất hiện trận lũ tần suất P = 1% và P = 0,4%.
    • Chủ thể và thời gian: Sở Nông nghiệp và Phát triển nông thôn tỉnh Thái Bình phê duyệt quy chuẩn và thực hiện cải tạo hạ tầng trong giai đoạn 2013–2018.

  3. Gia cố, nâng cấp cao trình mặt đê và xây dựng kè bảo vệ bờ sông:

    • Hành động: Nâng cấp các đoạn đê xung yếu có nền đất yếu; bê tông hóa mặt đê kết hợp giao thông nông thôn; trồng dải rừng ngập mặn phòng hộ phía ngoài đê biển và cửa sông.
    • Mục tiêu: Đảm bảo cao độ đỉnh đê luôn cao hơn mực nước lũ thiết kế tối thiểu từ 0,7 m đến 1,0 m (độ cao an toàn dự phòng); mở rộng thêm khoảng 500 ha đến 800 ha diện tích rừng phòng hộ chắn sóng.
    • Chủ thể và thời gian: Ban Quản lý Dự án Đầu tư Xây dựng các công trình Nông nghiệp tỉnh triển khai liên tục từ 2013 đến 2020.

  4. Hiện đại hóa hệ thống cảnh báo sớm và vận hành liên hồ chứa thông minh:

    • Hành động: Lắp đặt thêm 5 trạm quan trắc mực nước tự động truyền dữ liệu theo thời gian thực tại các cống đầu mối (Thiên Kiều, Sa Lung, Tam Lạc); ứng dụng phần mềm dự báo liên kết khí tượng - thủy văn - triều dâng.
    • Mục tiêu: Nâng cao thời gian dự báo cảnh báo lũ trước từ 24 đến 48 giờ, hạn chế thấp nhất thiệt hại về kinh tế và dân sinh.
    • Chủ thể và thời gian: Đài Khí tượng Thủy văn khu vực Đồng bằng Bắc Bộ phối hợp cùng tỉnh Thái Bình hoàn thiện trong giai đoạn 2014–2016.

Đối tượng nên tham khảo luận văn

Nội dung và dữ liệu thực nghiệm của luận văn mang lại giá trị thực tiễn to lớn cho nhiều nhóm độc giả chuyên môn:

  1. Cơ quan quản lý nhà nước về tài nguyên nước và đê điều (Sở NN&PTNT, Chi cục Thủy lợi và Phòng chống thiên tai):

    • Lợi ích: Cung cấp cơ sở khoa học định lượng về mực nước lũ thiết kế và hành lang thoát lũ để ban hành các quyết định quy hoạch bãi sông, quy hoạch sử dụng đất bãi bồi ven sông Trà Lý.

  2. Kỹ sư tư vấn thiết kế công trình thủy lợi, thủy điện và giao thông cầu đường:

    • Lợi ích: Sử dụng trực tiếp bảng số liệu mực nước thiết kế tần suất 1%, lưu lượng đỉnh lũ và các mặt cắt thủy lực để tính toán khẩu độ thoát nước cầu vượt sông, cao trình đáy dầm cầu, cao trình đỉnh cống tiêu ngập.

  3. Chuyên viên quy hoạch đô thị và phát triển nông thôn mới tại các địa phương ven sông:

    • Lợi ích: Nắm vững ranh giới an toàn thoát lũ để quy hoạch các cụm công nghiệp, khu dân cư mới tránh xa các vùng nguy hiểm ngập lũ hoặc vùng phân lũ chậm.

  4. Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Thủy văn học, Kỹ thuật Tài nguyên nước, Quản lý rủi ro thiên tai:

    • Lợi ích: Nguồn tài liệu tham khảo chuẩn mực về phương pháp ứng dụng mô hình số trị MIKE 11 trong bài toán thủy lực mạng sông đồng bằng chịu ảnh hưởng triều phức tạp.

Câu hỏi thường gặp

Hỏi: Vì sao cần phải lập quy hoạch phòng chống lũ chi tiết riêng cho sông Trà Lý? Trả lời: Do sông Trà Lý tiếp nhận từ 8% đến 11% lượng lũ của hệ thống sông Hồng đổ ra vịnh Bắc Bộ và trực tiếp chia đôi địa bàn tỉnh Thái Bình. Nếu không có quy hoạch chi tiết, việc lấn chiếm bãi sông bừa bãi sẽ làm nghẽn dòng chảy, gây mất an toàn đê chính bảo vệ hơn 1,7 triệu dân.

Hỏi: Mô hình MIKE 11 được ứng dụng trong nghiên cứu này có độ tin cậy như thế nào? Trả lời: Mô hình được thiết lập với dữ liệu địa hình mặt cắt đo đạc thực tế, hiệu chỉnh và kiểm định qua các chuỗi số liệu của nhiều trận lũ lịch sử (1969, 1971, 1996) tại 10 trạm thủy văn chính, đạt sai số mực nước đỉnh lũ dưới 5–10 cm, hoàn toàn đáp ứng tiêu chuẩn ngành.

Hỏi: Hồ chứa thủy điện Hòa Bình và Tuyên Quang tác động như thế nào đến đỉnh lũ sông Trà Lý? Trả lời: Các hồ chứa giúp cắt giảm trung bình từ 7 cm mực nước đỉnh lũ ở trạm Quyết Chiến (thượng nguồn sông). Tuy nhiên, tại khu vực hạ lưu ven biển, mực nước đỉnh lũ vẫn tăng từ 1 cm đến 2 cm do tác động tổng hợp của hiện tượng bồi lắng lòng dẫn và nước dâng do bão.

Hỏi: Mực nước lũ thiết kế tần suất 1% trên sông Trà Lý cao hơn đỉnh lũ lịch sử 1971 bao nhiêu? Trả lời: Theo kết quả tính toán thủy lực, mực nước thiết kế P = 1% tại trạm Quyết Chiến đạt 654 cm (cao hơn lũ năm 1971 là 9 cm), tại trạm Thái Bình đạt 444 cm (cao hơn đỉnh lũ năm 1971 khoảng 23 cm), phản ánh kịch bản an toàn tối ưu cho toàn tuyến.

Hỏi: Giải pháp nào hiệu quả nhất để xử lý mâu thuẫn giữa khai thác kinh tế bãi bồi và thoát lũ? Trả lời: Giải pháp cốt lõi là phân cấp quản lý đê bối rõ ràng kết hợp cắm mốc chỉ giới thoát lũ. Các bối sản xuất nông nghiệp chỉ duy trì mức bảo vệ ứng với báo động cấp 2 và cấp 3, sẵn sàng phương án cho tràn ngập có kiểm soát khi xuất hiện lũ thiết kế.

Kết luận

  • Luận văn đã phân tích toàn diện đặc điểm thủy văn dòng chảy sông Trà Lý, khẳng định vai trò gánh vác 8% đến 11% lượng lũ sông Hồng của tuyến sông dài 64 km này.
  • Thiết lập và kiểm định thành công mô hình thủy lực MIKE 11, tính toán chính xác mực nước lũ thiết kế P = 1% đạt 654 cm tại trạm Quyết Chiến và 444 cm tại trạm Thái Bình.
  • Làm rõ tác động trái chiều của hệ thống hồ chứa thượng lưu kết hợp mực nước triều thiết kế và nước dâng do bão lên đến 3,03 m tại vùng cửa sông.
  • Đưa ra hệ thống giải pháp phân loại đê bối khoa học, xác định hành lang thoát lũ nhằm bảo vệ 133.693 ha diện tích tự nhiên và 1,78 triệu dân tỉnh Thái Bình.
  • Đề xuất các cơ quan quản lý và địa phương nhanh chóng ứng dụng khung kết quả tính toán này vào công tác cắm mốc chỉ giới thoát lũ và lập kế hoạch đầu tư nâng cấp hệ thống đê điều giai đoạn tiếp theo.