Ảnh hưởng của điều tiết các hồ chưa thương lưu đến dòng chảy kiệt cấp nước vùng hạ du sông Cả

Để phục vụ phát triển kinh tế xã hội cho các địa phương trên lưu vực, một số hồ chứa trên thượng lưu đã và đang được xây dựng với mục đích cấp nước trong mùa khô, cắt lũ trong mùa m ưa và phát điện. Kết quả tính toán thủy lực mô tả dòng chảy kiệt trong mạng sông ứng với các phương án điều tiết và các tần suất khác nhau để làm cơ sở cho việc đề xuất giải pháp phù hợp, đảm bảo cấp nước cho việc phát triển nông nghiệp, thủy sản, các ngành kinh tế và môi trường.

KHOA HỌC CÔNG NGHỆ ẢNH HƯỞNG CỦA ĐIỀU TIẾT CÁC HỒ CHƯA THƯƠNG LƯU ĐẾN DÒNG CHẢY KIỆT CẤP NƯỚC VÙNG HẠ DU SÔNG CẢ PGS.TS Nguyễn Q uang Trung ThS Lương Ngọc Chung Viện Nước, Tưới tiêu Môi trường Tóm tắt: Sơng Cả đóng vai trị quan trọng phát triển kinh tế xã hội hai tỉnh Nghệ An Hà Tĩnh Tổng lượng nước trung bình hàng năm vào khoảng 21 đến 23 tỷ m ét khối Tuy nhiên mùa khô kéo dài tháng, lượng nước chiếm khoảng 25 đến 30%, đặc biệt tháng III IV mực nước sông hạ xuống thấp ảnh hưởng lớn đến việc cấp nước cho ngành kinh tế Để phục vụ phát triển kinh tế xã hội cho địa phương lưu vực, số hồ chứa thượng lưu xây dựng với mục đích cấp nước mùa khô, cắt lũ mùa m ưa phát điện Kết tính tốn thủy lực m tả dịng chảy kiệt mạng sơng ứng với phương án điều tiết tần suất khác để làm sở cho việc đề xuất giải pháp phù hợp, đảm bảo cấp nước cho việc phát triển nông nghiệp, thủy sản, ngành kinh tế môi trường Summary: Ca River plays an important role in socio-econom ic development of Nghe An and Ha Tinh provinces Total annual discharge is about of 21 to 23 m illion m However, in m onths of dry season, water discharge occupies of 25 to 30% only, especially in March and April, water level drops down very low which impact to water supply for econom ic sectors In order to develop socio-econom ic in the basin, reservoirs in upstream have been built for flood control in rainy season, water supply in dry season and electricity generation Results of hydraulic sim ulation with the operation options and frequencies to determine proper m easures for ensuring water supply in dry season in lower river basin I ĐẶT VẤN ĐỀ1 Sông Cả m ột hệ thống sơng liên quốc gia, có tổng diện tích tự nhiên tồn lưu vực 27.200km 2, phần diện tích nằm đất Cộng hoà Dân chủ Nhân dân Lào 9.470km2, đất Việt Nam 17.730km2 Sông Cả gồm nhiều nhánh sông nhỏ nhập lưu sông Hiếu, sông Giăng, sông La (bao gồm Ngàn Sâu, Ngàn Phố) Vùng hưởng lợi từ hệ thống sông Cả vùng chịu tác động nguồn nước sông Cả nằm chủ yếu hạ du sông thuộc địa bàn hai tỉnh Nghệ An Hà Tĩnh Mặc dù diện tích lưu vực sơng Cả lớn, nguồn nước dồi dào, trung bình m ột năm tải biển tổng lượng nước từ 21-23 tỷ m3 Người phản biện: Ngày nhận bài: Ngày thông qua phản biện: Ngày duyệt đăng: phân bố không theo thời gian Trong ba tháng mùa m ưa tổng lượng dòng chảy chiếm tới 15-16 tỷ m3, mùa khơ tháng tổng lượng dịng chảy từ 6-7 tỷ m3 làm cho m ực nước sông bị hạ thấp, mặn xâm nhập sâu ảnh hưởng lớn đến khả khai thác cơng trình thủy lợi dọc sơng gây nên tình trạng hạn hán, mùa khó khăn việc cấp nước cho ngành kinh tế Mục đích nghiên cứu tính tốn thuỷ lực nhằm xác định mực nước, lưu lượng mùa kiệt cho hệ thống sơng Cả nhằm mơ tả chế độ dịng chảy mùa kiệt mạng sông ứng với tần suất thiết kế 75%; 85% 90% cho giai đoạn năm 2020 trường hợp có tham gia điều tiết hồ chứa thượng lưu để làm sở cho việc đề xuất giải pháp phù hợp, đảm bảo cấp nước cho việc phát triển nông nghiệp, thủy sản, ngành TẠP CH Í KH OA H ỌC VÀ CÔN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 43 KHOA HỌC CÔNG NGHỆ Choang, sông Trai, suối Rộ, sông Rào Gang, sông Cầu Nậy… kinh tế môi trường II PHẠM VI NGHIÊN CỨU Phạm vi nghiên cứu bao gồm: Sơng Cả (tồn III ĐẶC ĐIỂM DỊNG CHÍNH SƠNG CẢ dịng thuộc Việt Nam ), tồn lưu Đặc trưng hình thái lưu vực sông Cả thể vực sông Hiếu, sông Giăng, Sông La (bao gồm qua tiêu đánh giá lưu vực nhánh Sông Ngàn sâu, Ngàn Phố) Bảng nhánh sông suối nhỏ: Huổi Nguyên, Khe Bảng 1: Đặc trưng hình thái m ột số lưu vực sông TT Luu vực F (km2) Lsôn g (km ) Độ cao bq(m) Sông Cả S Nậm Mô S Giăng Sông Hiếu Sông La 27.200 3.970 1.050 5.340 3.210 531 173 77 228 135 294 960 492 303 362 Độ dốc bqlv (%o) 1,83 2,57 1,72 1,30 2,82 Bbq km /k m2 89 38,2 15,8 32,5 46,6 Hệ số Mật số lưới không sông đối km /km xứng 0,60 -0,14 0,22 -0,09 0,71 0,02 0,87 0,53 Hệ số hình dạng lưu vực 0,29 0,27 0,24 0,20 0,68 IV MÔ HÌNH THỦY LỰC HỆ THỐNG SƠNG trí bắt đầu có đê đến nhập lưu vào sông Cả 4.1 Sơ đồ mạng sông: + Sông Ngàn Phố: Từ trạm thuỷ văn Sơn Diệm đến ngã ba Linh Cảm Sông Cả m ột hệ thống sơng hồn chỉnh, nghiên cứu toán thủy lực m ùa kiệt hạ lưu cần phải xem xét tính tốn tồn m ạng sơng Mạng tính tốn thuỷ lực bao gồm tồn dịng phụ lưu vùng trung, hạ du lưu vực sông Cả, cụ thể sau: + Sông Ngàn Sâu: Từ trạm thuỷ văn Hoà Duyệt đến ngã ba Linh Cảm + Sông La: Từ Linh Cảm đến nhập lưu vào sông Cả ngã ba Chợ Tràng + Dịng sơng Cả: Từ ngã ba Cửa Rào đến cửa sông Cửa Hội + Sông Hiếu: Từ trạm thuỷ văn Nghĩa Khánh đến nhập lưu vào sông Cả ngã ba Cây Chanh + Sông Giăng: Từ tuyến Thác Muối đến nhập lưu vào sông Cả + Sông Gang: Tại vị 4.2 C ác biên tính tốn 44 Hình Sơ đồcủa mạng a Biên mơsơng hình: TẠP C H Í KH OA HỌC VÀ CÔN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 KHOA HỌC CÔNG NGHỆ Với m ạng sơng tính tốn xác định trên, biên mơ hình thuỷ lực trình lưu lượng theo thời gian Q=f(t) cụ thể sau: + Trên sông Ngàn Phố trạm thủy văn Sơn Diệm (Flv= 790 km 2) + Biên lưu lượng trạm thủy văn Cửa Rào (Flv= 12.800 km 2) dịng sơng Cả b Biên m hình: + Trên sơng Hiếu biên lưu lượng Nghĩa Khánh (Flv= 4.020 km2) + Trên sông Giăng biên lưu lượng Thác Muối (Flv= 785 km2) + Trên sơng Rào Gang vị trí đầu đê (Flv= 478 km2) + Trên sông Ngàn Sâu trạm thủy văn Hoà Duyệt (Flv= 1.880 km2) Biên mơ hình thuỷ lực q trình m ực nước theo thời gian Z=f(t) Cửa Hội c Biên dọc sơng m hình: Các biên gia nhập vào trung lưu sông Cả sông khác: Huổi Nguyên, Khe Choang, Hói Quất, khu khác sơng Cả, sông Hiếu, sông Giăng, sông Ngàn Sâu, Ngàn Phố Bảng 2: Chỉ tiêu biên gia nhập khu sông Cả TT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 Biên gia nhập Tả sông Cả Huổi Nguyên Hữu sông Cả Khe Choang Hữu sông Cả Tả sông Cả Tả sơng Cả Hói Quất Tả sơng Cả Hữu sông Sông Trai Sông Rỗ Khe Lung Hữu sông Cả Tả sông Hiếu Tả sông Hiếu Tả sông Hiếu Hữu sông Hiếu Hữu sông Hiếu Khu sông Giăng Khu sông Rào Gang Khu sông Ngàn Phố Khu sông Ngàn Sâu Flv (km 2) 130 800 156 431 147 351 150 152 180 51,5 201,5 125 219 79 177 150 134 284 569 265 183 270 188 Thuộc sông Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Hiếu Sông Hiếu Sông Hiếu Sông Hiếu Sông Hiếu Sông Giăng Sông Rào Gang Sông Ngàn Phố Sông Ngàn Sâu d Các vị trí lấy nước khu giữa: TẠP CH Í KH OA H ỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 47 KHOA HỌC CÔNG NGHỆ Bảng 3: Chỉ tiêu vị trí lấy nước dọc sơng Các cụm dùng nước Các trạm bơm Anh Sơn Các trạm bơm Thanh Chương Trên Rào Gang Trên sông Lam -Thượng sông Giăng - Hạ sông Giăng - Gần cửa sơng Gang - TB sơng Giăng Các vị trí lấy nước Đô Lương - Trạm bơm - Hệ thống Đơ Lương Các vị trí lấy nước Nam Đàn - Cống Nam Đàn - Hệ thống trạm bơm - Hệ thống trạm bơm Các trạm bơm Hưng Nguyên - Tại Yên Xuân - Tại TV Chợ Tràng TP Vinh Hệ thống lấy nước Đức Thọ - TB Linh Cảm - TB nhỏ - Cống Đức Xá - Cống Trung Lương TB Sông Ngàn Phố TB Sông Ngàn Sâu Tổng SƠ ĐỒ C ÁC CƠN G TRÌN H C HÍN H LẤ Y NƯỚC DỌC SƠN G C Ả Trên sông Cả Lưu lượng (m 3/s) Tương lai 1,739 2,635 Rào Gang 2,148 2,460 Cả Cả Cả Giăng 1,509 2,266 0,540 1,486 2,182 3,324 1,932 1,958 Cả Cả 5,854 30,417 7,245 44,000 Cả Cả Cả 23,000 1,395 2,693 44,040 2,082 3,245 Cả Lam Lam La La La La Hào Ngàn phố Ngàn Sâu 1,345 1,402 0,114 1,541 1,990 0,214 10,158 0,617 3,600 6,500 2,464 0,735 99,981 Các trạm kiểm tra dọc sông: 0,000 0,750 4,970 8,970 3,881 1,446 138,880 Chợ Tràng, Bến Thủy, Hưng Nhân, Hưng Hoà, Rào Đừng, (sông Lam ); Linh Cảm (sông La); Trung Lương (sơng Hào) 4.3 Khai thác m hình MIKE11 tính tốn thủy lực m ùa kiệt 4.3.1 Tính tốn m phỏng: Nhằm xác định thơng số dịng chảy xâm nhập mặn m ùa kiệt cho vùng hạ du sông Cả, nghiên cứu chúng tơi chọn thời kỳ từ ngày 07-20/04/1989 để tính tốn m ô Đây thời đoạn triều thường bị kiệt nhất, có nhu cầu nước 48 cao năm lưu vực, đồng thời có số TẠP C H Í KH OA HỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 KHOA HỌC liệu quan trắc đầy đủ, đồng CƠNG NGHỆ số trạm có quan trắc m ực nước thời kỳ kiệt sông Cả sau: Kết tính tốn mơ thực đo Bảng 4: Kết mực nước thực đo tính tốn m TT Vị Trí Sơng TV n Thượng TV Chợ Tràng TV Bến Thủy Hưng Nhân Hưng Hoà Rào Đừng TV Linh Cảm Trung Lương Cả Lam Lam Lam Lam Lam La Hào Hmax (m) Thực đo Tính toán 1,97 1,98 1,15 1,05 1,08 1,01 1,08 1,03 1,16 1,03 1,06 1,05 1,18 1,08 1,04 1,04 Đường trình m ực nước tính tốn mơ kiệt năm 1989 thực đo m ột số vị trí sơng Cả 4.3.2 Tính tốn kiểm định m hình: Sai số 0,01 0,10 0,07 0,05 0,13 0,01 0,10 0,00 Thực đo 1,76 -1,07 -1,34 -1,38 -1,35 -1,25 -0,55 -1,34 Hmin (m) Tính tốn 1,79 -1,10 -1,14 -1,27 -1,27 -1,33 -0,34 -1,10 Sai số 0,03 0,03 0,20 0,11 0,08 0,08 0,21 0,24 Nghiên cứu chọn thời kỳ kiệt có số liệu thực đo từ 25-30/07/2011 để tính tốn kiểm định Kết tính tốn m thực đo m ột số trạm thủy văn sông Cả sau: Bảng 5: Kết mực nước thực đo tính tốn m TT Vị Trí TV Chợ Tràng TV Bến Thủy Hưng Nhân Rào Đừng C Đức Xá Sông Lam Lam Lam Lam La Hmax (m) Thực đo 1,16 1,05 1,11 1,10 1,25 Tính tốn 1,09 1,03 1,07 0,99 1,22 Kết mô kiệt từ 07-20/04/1989 kiểm định từ 25-30/07/2011 cho thấy: - Mực nước, nồng độ m ặn tính tốn giá trị thực đo chênh khơng đáng kể; đường q trình diễn biến m ực nước thời kỳ kiệt năm 1989 kiểm định 2011 q trình tính tốn q trình thực đo tương đối bám sát - Kết tính tốn mơ kiểm định phù hợp với thực tế khảo sát, thông số sử dụng mơ hình thuỷ lực phản ánh xác chế độ thủy lực m ùa kiệt mạng sông Cả, đủ độ tin cậy để tiến hành tính tốn thuỷ lực kiệt, m ặn cho trường hợp nghiên cứu hệ thống sông Cả Sai số 0,07 0,02 0,04 0,11 0,03 Hmin (m) Thực đo 0,89 -1,01 -0,85 -1,11 -0,30 Tính tốn 0,81 -0,89 -0,85 -1,08 -0,31 Sai số 0,08 0,08 0,00 0,03 0,01 4.3.3 Các phương án tính tốn: a Tính tốn khơng có hồ chứa bổ sung lưu lượng mùa kiệt cho hạ du: - Trường hợp khai thác nguồn nước sông Cả phục vụ nhu cầu nước khơng có cơng trình hồ chứa thượng nguồn bổ sung lưu lượng cho hạ du m ùa kiệt ứng với tần suất (Ký hiệu HT-75%; HT 85% HT 90%); - Trường hợp khai thác nguồn nước sông Cả phục vụ nhu cầu nước tương lai 2020, khơng có cơng trình hồ chứa thượng nguồn bổ sung lưu lượng cho hạ du mùa kiệt với tần suất (ký hiệu: TL2020-75%; TL202085% TL2020-90%) b Tính tốn có hồ chứa bổ sung lưu lượng TẠP CH Í KH OA H ỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 49 KHOA HỌC CÔNG NGHỆ mùa kiệt cho hạ du: THC90%-3) với trường hợp: Trường hợp 1: - Như trường hợp Có hồ Bản Vẽ: Nhiệm vụ cơng trình phát điện, tham gia cắt lũ cho hạ du bổ sung nguồn nước tham gia đẩy m ặn cho hạ du sông Cả mùa kiệt với lưu lượng 80m 3/s ứng với tần suất (Ký hiệu THC75%-1, THC85%-1 THC90%-1) - Xây dựng hồ Thác Muối nhánh sông Giăng tưới bổ sung nguồn nước cho hạ du mùa kiệt Q=18,5m 3/s Trường hợp (Ký hiệu THC75%-2; 85%-2 THC90%-2 ), với trường hợp: - Hồ Bản Vẽ bổ sung xuống hạ du sông Cả mùa kiệt Q=80m3/s - Xây dựng hồ Đá Gân, Khe Chè, Nước Sốt, Rào Mắc nhánh Ngàn Phố tưới bổ sung nguồn nước cho hạ du m ùa kiệt Q=11,5m 3/s - Hồ Bản Mồng khởi công xây dựng năm 2010 nhánh sông Hiếu, với nhiệm vụ bổ sung nước cho hạ du sông Cả m ùa kiệt 22m3/s Tóm lại: Trên dịng sơng Cả có hồ bổ sung cho hạ du với lưu lượng 120,5m3/s; Trên nhánh sông La có hồ bổ sung với lưu lượng 24,9m 3/s - Hồ Ngàn Trươi khởi công xây dựng năm 2009 nhánh Ngàn Trươi, với nhiệm vụ bổ sung cho hạ du m ùa kiệt với lưu lượng 4m 3/s 4.3.4 Kết tính tốn kiệt với tần suất P=75% - Xây dựng hồ Chúc A, Trại Dơi nhánh Ngàn Sâu tưới bổ sung nguồn nước cho hạ du mùa kiệt Q=9,4m3/s Kết tính tốn mực nước, lưu lượng vị trí hệ thống sông Cả sau: Trường hợp (Ký hiệu THC75%-3;THC85%Bảng 6: Mực nước max, số vị trí sơng - P=75% Đơn vị: m Vị trí Sơng Cả Thủy văn Dừa TL Bara Đô Lương HL Bara Đô Lương TV Yên Thượng TV Nam Đàn Cống Nam Đàn T.bơm Hưng Lĩnh T.bơm Hưng Long Cầu Yên Xuân Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh H Châu-H Nhân TV Bến Thủy Cống Hưng Hoà Cống Rào Đừng TV Cửa Hội HT-75% TL2020-75% Min Max Tương lai 2020 THC75%-1 THC75%-2 Min Max Min Max THC75%-3 Min Max Min Max 13.54 9.88 6.13 1.34 0.95 0.93 0.07 -0.48 -0.54 13.66 9.92 6.33 1.53 1.12 1.10 0.80 0.78 0.76 13.54 9.86 6.04 1.08 0.29 0.29 -0.22 -0.62 -0.64 13.66 9.91 6.26 1.31 0.64 0.64 0.73 0.75 0.75 13.99 10.06 6.81 2.01 1.36 1.33 0.30 -0.33 -0.43 14.07 10.10 6.93 2.12 1.47 1.44 0.87 0.83 0.81 14.08 10.11 6.96 2.16 1.51 1.47 0.41 -0.24 -0.36 14.16 10.15 7.07 2.25 1.60 1.56 0.91 0.85 0.83 14.08 10.11 6.96 2.27 1.62 1.57 0.49 -0.17 -0.31 14.16 10.15 7.08 2.36 1.71 1.66 0.95 0.87 0.85 -0.77 -0.76 -0.79 -0.82 -0.81 -0.86 -0.91 0.72 0.71 0.73 0.74 0.76 0.79 0.80 -0.80 -0.78 -0.80 -0.83 -0.82 -0.87 -0.91 0.73 0.73 0.75 0.75 0.79 0.80 0.80 -0.74 -0.75 -0.77 -0.82 -0.81 -0.86 -0.91 0.78 0.76 0.77 0.76 0.79 0.81 0.80 -0.72 -0.74 -0.76 -0.81 -0.80 -0.86 -0.91 0.80 0.78 0.77 0.76 0.80 0.81 0.80 -0.69 -0.72 -0.74 -0.80 -0.80 -0.86 -0.91 0.82 0.79 0.78 0.77 0.80 0.81 0.80 TẠP C H Í KH OA HỌC VÀ CÔN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 49 KHOA HỌC CÔNG NGHỆ HT-75% Vị trí Min Sơng La Linh Cảm CầuChợ Thượng Cống Đức Xá Liên Minh-Đ.Tùng Ngã ba Chợ Tràng Sông Hào Cống Trung Lương TL2020-75% Min Max Max Tương lai 2020 THC75%-1 THC75%-2 Min Max Min Max THC75%-3 Min Max -0.41 -0.42 -0.45 -0.47 -0.77 0.71 0.70 0.70 0.71 0.72 -0.38 -0.41 -0.45 -0.48 -0.80 0.73 0.73 0.73 0.72 0.73 -0.36 -0.38 -0.43 -0.46 -0.74 0.78 0.78 0.78 0.77 0.78 -0.34 -0.36 -0.41 -0.44 -0.72 0.80 0.80 0.79 0.78 0.80 -0.22 -0.26 -0.34 -0.38 -0.69 0.85 0.84 0.82 0.81 0.82 -0.72 0.70 -0.79 0.72 -0.74 0.76 -0.73 0.78 -0.69 0.80 Bảng 7: Lưu lượng m ax, số vị trí sơng - P=75% Đơn vị: m3/s Vị trí Sơng Cả Thủy văn Dừa TL Para Đơ Lương HL Para Đô Lương TV Yên Thượng Cống Nam Đàn Trạm bơm Hưng Lĩnh Trạm bơm Hưng Long Cầu Yên Xuân Ngã ba Chợ Tràng Sông Giăng Thác Muối Cửa sông Giăng Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh Hưng Châu-H Nhân Thủy văn Bến Thủy Hưng Hồ Cống Rào Đừng Sơng Ngàn Phố Sơn Diệm Cửa sông Ngàn Phố Sông Ngàn Sâu Hồ Duyệt Cửa sơng Ngàn Sâu Sơng La TV Linh Cảm Cống Đức Xá Cửa sông La Sông Hào Cống Trung Lương Min Max TL2020-75% Min Max Tương lai 2020 THC75%-1 THC75%-2 Min Max Min Max 74 77 50 53 36 -21 -76 -166 -419 91 94 67 69 56 70 92 138 275 74 77 43 46 -17 -67 -110 -185 -413 91 94 61 62 18 39 61 116 255 154 157 123 127 82 43 -23 -98 -382 171 174 141 143 98 108 121 170 301 176 179 145 149 104 72 22 -75 -370 193 196 163 165 121 129 140 186 315 176 179 145 167 122 93 51 -56 -364 193 196 163 183 139 146 157 200 329 6 6 6 6 24 24 25 26 -714 -940 -1006 -1540 -1445 -1831 545 709 819 1023 1234 1510 -645 -839 -943 -1283 -1411 -1768 524 683 794 999 1208 1484 -610 -812 -898 -1256 -1365 -1724 572 727 836 1036 1245 1520 -596 -800 -880 -1244 -1347 -1707 587 740 849 1048 1255 1530 -573 -780 -851 -1223 -1320 -1685 609 762 870 1067 1275 1550 -12 11 20 -10 11 20 -10 11 21 -10 11 21 20 22 31 16 -67 24 68 16 -62 24 66 16 -68 24 68 20 -66 28 71 29 -59 37 78 -89 -188 -244 94 155 208 -78 -180 -238 92 155 208 -85 -193 -254 95 158 210 -85 -192 -255 98 160 211 -66 -175 -239 114 173 221 -25 19 -24 19 -24 20 -24 20 -22 21 HT-75% 50 THC75%-3 Min Max TẠP CH Í KH OA H ỌC VÀ C ÔN G N GHỆ THỦ Y LỢI SỐ 20 - 2014 KHOA HỌC CÔNG NGHỆ 4.3.5 Kết tính tốn kiệt với tần suất P=85% : Kết tính tốn mực nước, lưu lượng vị trí hệ thống sông Cả sau: Bảng 8: Mực nước max, số vị trí sơng - P=85% Đơn vị: m HT-85% Vị trí Min Sơng Cả Thủy văn Dừa TL Bara Đô Lương HL Bara Đô Lương TV Yên Thượng TV Nam Đàn C Nam Đàn T.bơm Hưng Lĩnh T.bơm Hưng Long Cầu Yên Xuân Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh H Châu-H Nhân TV Bến Thủy Cống Hưng Hoà Cống Rào Đừng TV Cửa Hội Sông La Linh Cảm Cầu Chợ Thượng Cống Đức Xá Liên Minh-Đ.Tùng Ngã ba Chợ Tràng Sông Hào Cống Trung Lương TL2020-85% Min Max Max Tương lai 2020 THC85%-1 THC85%-2 Min Max Min Max THC85%-3 Min Max 13.49 13.61 13.49 9.86 9.90 9.85 6.04 6.23 5.94 1.21 1.39 0.95 0.83 1.02 -0.21 0.82 1.00 -0.21 0.01 0.77 -0.30 -0.52 0.76 -0.65 -0.57 0.74 -0.67 13.61 9.88 6.15 1.17 0.47 0.47 0.70 0.73 0.73 13.95 10.04 6.75 1.94 1.29 1.26 0.25 -0.37 -0.44 14.03 10.08 6.87 2.06 1.43 1.39 0.87 0.83 0.81 14.05 10.09 6.91 2.10 1.46 1.42 0.37 -0.28 -0.38 14.12 10.13 7.01 2.20 1.56 1.52 0.90 0.85 0.83 14.05 10.09 6.91 2.22 1.57 1.53 0.45 -0.20 -0.32 14.12 10.13 7.02 2.31 1.67 1.62 0.94 0.87 0.85 -0.78 -0.77 -0.79 -0.83 -0.82 -0.87 -0.91 0.71 0.71 0.73 0.74 0.75 0.79 0.80 -0.81 -0.78 -0.81 -0.84 -0.82 -0.87 -0.91 0.72 0.73 0.75 0.75 0.78 0.80 0.80 -0.75 -0.75 -0.78 -0.82 -0.81 -0.86 -0.91 0.78 0.76 0.76 0.76 0.79 0.81 0.80 -0.72 -0.74 -0.76 -0.81 -0.80 -0.86 -0.91 0.80 0.78 0.77 0.76 0.80 0.81 0.80 -0.69 -0.72 -0.75 -0.80 -0.80 -0.86 -0.91 0.82 0.79 0.78 0.77 0.80 0.81 0.80 -0.43 -0.44 -0.46 -0.48 -0.78 0.70 0.68 0.69 0.70 0.71 -0.46 -0.47 -0.50 -0.51 -0.81 0.70 0.71 0.71 0.71 0.72 -0.44 -0.44 -0.47 -0.49 -0.75 0.76 0.76 0.77 0.76 0.78 -0.36 -0.38 -0.42 -0.45 -0.72 0.80 0.80 0.79 0.78 0.80 -0.24 -0.28 -0.36 -0.40 -0.69 0.85 0.84 0.82 0.81 0.82 -0.72 0.70 -0.80 0.71 -0.76 0.76 -0.73 0.78 -0.70 0.80 Bảng 9: Lưu lượng m ax, số vị trí sơng - P=85% Vị trí Sơng Cả Thủy văn Dừa TL Para Đô Lương HL Para Đô Lương TV Yên Thượng Cống Nam Đàn Trạm bơm Hưng Lĩnh Trạm bơm Hưng Long Cầu Yên Xuân Ngã ba Chợ Tràng HT-85% Min Max TL2020-85% Min Max 67 70 43 45 25 -34 -85 -176 -422 83 86 59 59 46 62 87 134 270 67 70 36 38 -22 -76 -119 -189 -413 83 86 52 52 31 53 108 248 TẠP C H Í KH OA HỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 Đơn vị: m3/s Tương lai 2020 THC85%-1 THC85%-2 Min Max Min Max THC85%-3 Min Max 147 150 116 118 73 42 -22 -96 -380 169 172 138 159 113 92 51 -56 -364 51 163 166 132 133 90 104 119 162 295 169 172 138 140 95 71 22 -75 -370 185 188 154 155 111 123 136 179 309 185 188 154 174 130 141 152 193 322 KHOA HỌC CƠNG NGHỆ HT-85% Vị trí Sơng Giăng Thác Muối Cửa sông Giăng Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh Hưng Châu-H Nhân Thủy văn Bến Thủy Hưng Hồ Cống Rào Đừng Sơng Ngàn Phố Sơn Diệm Cửa sông Ngàn Phố Sông Ngàn Sâu Hồ Duyệt Cửa sơng Ngàn Sâu Sơng La TV Linh Cảm Cống Đức Xá Cửa sông La Sông Hào Cống Trung Lương TL2020-85% Min Max Tương lai 2020 THC85%-1 THC85%-2 Min Max Min Max THC85%-3 Min Max Min Max 5 6 5 6 5 6 5 6 23 23 24 25 -719 -944 -1011 -1544 -1451 -1836 537 701 811 1016 1226 1502 -659 -852 -953 -1301 -1422 -1777 511 669 780 986 1193 1468 -612 -814 -907 -1256 -1375 -1731 559 714 823 1025 1234 1509 -596 -800 -884 -1244 -1351 -1710 579 733 841 1041 1248 1523 -573 -780 -856 -1223 -1324 -1688 601 754 862 1060 1268 1543 -14 19 -13 18 -13 20 -14 20 19 21 30 14 -68 20 65 14 -62 20 65 14 -68 20 67 18 -68 24 69 27 -61 34 75 -92 -192 -246 90 151 205 -81 -193 -247 89 146 200 -87 -205 -264 92 149 202 -90 -197 -258 94 157 209 -68 -180 -243 110 170 218 -25 18 -25 19 -25 19 -24 20 -22 21 4.3.6 Kết tính tốn kiệt với tần suất P=90% : Kết tính tốn mực nước, lưu lượng vị trí hệ thống sơng Cả sau: Bảng 10: Mực nước m ax, số vị trí sơng - P=90% Đơn vị: m Vị trí HT-90% Min Sơng Cả Thủy văn Dừa TL Bara Đô Lương HL Bara Đô Lương TV Yên Thượng TV Nam Đàn C Nam Đàn T.bơm Hưng Lĩnh T.bơm Hưng Long Cầu Yên Xuân Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh H Châu-H Nhân TV Bến Thủy Cống Hưng Hoà Cống Rào Đừng Max TL2020-90% Min Max Tương lai 2020 THC90%-1 THC90%-2 Min Max Min Max THC90%-3 Min Max 13.46 13.57 13.46 9.85 9.89 9.84 5.97 6.17 5.88 1.13 1.31 0.85 0.75 0.95 -0.72 0.74 0.94 -0.69 -0.02 0.71 -0.31 -0.55 0.67 -0.68 -0.59 0.67 -0.69 13.57 9.87 6.08 1.05 0.28 0.28 0.58 0.62 0.63 13.94 10.03 6.72 1.89 1.25 1.22 0.22 -0.39 -0.47 14.00 10.06 6.83 2.00 1.35 1.32 0.77 0.74 0.72 14.04 10.08 6.88 2.07 1.42 1.39 0.34 -0.30 -0.40 14.10 10.12 6.98 2.16 1.51 1.47 0.81 0.76 0.74 14.04 10.08 6.88 2.07 1.42 1.39 0.34 -0.30 -0.39 14.10 10.12 6.98 2.16 1.51 1.47 0.82 0.77 0.75 -0.79 -0.77 -0.79 -0.83 -0.82 -0.87 0.64 0.64 0.68 0.69 0.71 0.76 -0.76 -0.76 -0.78 -0.82 -0.81 -0.86 0.69 0.67 0.68 0.68 0.69 0.74 -0.74 -0.75 -0.77 -0.82 -0.81 -0.86 0.71 0.68 0.69 0.68 0.70 0.74 -0.73 -0.74 -0.77 -0.81 -0.81 -0.86 0.72 0.69 0.70 0.69 0.70 0.75 0.65 0.64 0.67 0.67 0.69 0.74 -0.82 -0.79 -0.81 -0.84 -0.83 -0.87 52 TẠP CH Í KH OA H ỌC VÀ C ÔN G N GHỆ THỦ Y LỢI SỐ 20 - 2014 KHOA HỌC HT-90% Vị trí TV Cửa Hội Sông La Linh Cảm Cầu Chợ Thượng Cống Đức Xá Liên Minh-Đ.Tùng Ngã ba Chợ Tràng Sơng Hào Cống Trung Lương CƠNG NGHỆ Tương lai 2020 THC90%-1 THC90%-2 Min Max Min Max -0.91 0.80 -0.91 0.80 THC90%-3 Min Max -0.91 0.80 Min -0.91 Max 0.80 TL2020-90% Min Max -0.91 0.80 -0.44 -0.45 -0.47 -0.49 -0.79 0.60 0.60 0.60 0.63 0.65 -0.47 -0.48 -0.51 -0.52 -0.82 0.58 0.59 0.59 0.62 0.64 -0.45 -0.46 -0.48 -0.50 -0.76 0.62 0.62 0.62 0.66 0.69 -0.43 -0.43 -0.46 -0.48 -0.74 0.66 0.65 0.64 0.69 0.71 -0.32 -0.34 -0.40 -0.44 -0.73 0.72 0.70 0.69 0.71 0.72 -0.73 0.64 -0.81 0.63 -0.77 0.66 -0.75 0.68 -0.73 0.70 Bảng 11: Lưu lượng m ax, m ột số vị trí sơng - P=90% Đơn vị: m3/s Vị trí Sơng Cả Thủy văn Dừa TL Para Đô Lương HL Para Đô Lương TV Yên Thượng Cống Nam Đàn Trạm bơm Hưng Lĩnh Trạm bơm Hưng Long Cầu Yên Xuân Ngã ba Chợ Tràng Sông Giăng Thác Muối Cửa sông Giăng Sông Lam Ngã ba Chợ Tràng TB Hưng Khánh Hưng Châu-H Nhân Thủy văn Bến Thủy Hưng Hoà Cống Rào Đừng Sông Ngàn Phố Sơn Diệm Cửa sông Ngàn Phố Sơng Ngàn Sâu Hồ Duyệt Cửa sơng Ngàn Sâu Sông La TV Linh Cảm Cống Đức Xá Cửa sông La Sông Hào Cống Trung Lương Min Max Tương lai 2020 TL2020-90% THC90%-1 Min Max Min Max 63 65 38 40 17 -49 -99 -183 -446 78 80 53 53 40 61 86 134 267 63 65 31 34 -16 -78 -121 -195 -452 78 80 46 47 27 51 106 247 143 145 111 113 68 17 -51 -124 -417 158 160 126 127 82 95 112 158 292 165 167 133 135 90 51 -7 -80 -405 180 182 148 149 104 114 128 174 305 165 167 133 135 90 49 -10 -84 -413 180 182 148 149 104 114 128 175 306 5 6 5 6 5 6 5 6 5 6 -793 -1115 -1253 -1965 -2141 -2573 532 696 806 1012 1222 1497 -772 -1089 -1200 -1866 -1962 -2365 507 664 777 983 1188 1462 -769 -1108 -1241 -1956 -2129 -2563 554 709 818 1021 1229 1504 -754 -1101 -1229 -1947 -2116 -2553 568 723 831 1032 1240 1514 -745 -1094 -1220 -1939 -2105 -2544 580 735 844 1043 1252 1528 -13 17 -14 17 -15 18 -16 18 18 19 27 12 -67 18 64 12 -63 18 63 12 -70 18 65 16 -69 22 68 26 -65 32 74 -90 -198 -259 88 149 203 -85 -197 -251 87 143 198 -95 -211 -268 89 147 201 -96 -212 -270 92 149 202 -72 -195 -258 107 161 212 -27 18 -25 19 -27 19 -27 19 -25 20 HT-90% THC90%-2 Min Max THC90%-3 Min Max Ghi chú: Lưu lượng âm thể dòng chảy ngược TẠP C H Í KH OA HỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 53 KHOA HỌC CƠNG NGHỆ V PHÂN TÍCH KẾT Q UẢ TÍNH TỐN CÁC TRƯỜNG HỢP KIỆT LƯU VỰC SƠNG CẢ VÀ ĐÁNH GIÁ KHẢ NĂNG LẤY NƯỚC CỦA CÁC CƠNG TRÌNH HẠ DU Qua kết tính tốn thủy lực kiệt tần suất kiệt (75%, 85% 90%) lưu vực sơng trường hợp tính tốn cho thấy: - Vùng chịu ảnh hưởng mạnh chế độ thủy triều: Trên sông Lam từ ngã ba Chợ Tràng đến cửa Hội, sông La từ ngã ba Linh Cảm đến ngã ba Chợ Tràng khơng có biến động m ạnh vùng chịu tác động m ạnh thủy triều: + Các trường hợp TL2020-75%, THC75%-1, THC75%-2, THC75%-3 ngã ba Chợ Tràng m ực nước nhỏ biến động từ -0,80÷0,69m , m ực nước lớn biến động từ 0,73÷0,82m; cống Đức Xá sơng La m ực nước nhỏ biến động từ -0,45÷0,34m , m ực nước lớn biến động từ 0,73÷0,82m Về lưu lượng: Đây vùng chịu ảnh hưởng chế độ thủy triều mạnh Khi lưu lượng thượng nguồn mùa khô nhỏ nên hoạt động thủy triều thắng sinh dòng chảy ngược (lưu lượng âm) thủy triều lên Khi thủy triều rút, lưu lượng theo chảy biển (lưu lượng dương) bao gồm lưu lượng sông phần lưu lượng đẩy vào triều lên Vì vùng hạ du từ Linh Cảm , ngã ba Chợ Tràng trở xuống lưu lượng tương đối dồi Tuy nhiên, chưa có hồ chứa thượng nguồn bổ sung nước vùng bị mặn uy hiếp mạnh gây khó khăn cho việc lấy nước (diễn biến mặn hạ du sơng Cả trình bày m ột viết khác) + Các trường hợp TL2020-85%, PA75%-1, PA85%-2, PA85%-3 ngã ba Chợ Tràng m ực nước nhỏ biến động từ -0,81÷0,69m , m ực nước lớn biến động từ 0,72÷0,82m; cống Đức Xá sông La 54 m ực nước nhỏ biến động từ -0,50÷0,36m, mực nước lớn biến động từ 0,71÷0,82m + Các trường hợp TL2020-90%, PA75%-1, PA90%-2, PA90%-3 ngã ba Chợ Tràng m ực nước nhỏ biến động từ -0,82÷0,73m, mực nước lớn biến động từ 0,68÷0,72m ; cống Đức Xá sơng La m ực nước nhỏ biến động từ -0,49÷0,40m, mực nước lớn biến động từ 0,65÷0,69m Tương tự trường hợp 75% lưu lượng từ Linh Cảm, ngã ba Chợ Tràng trở xuống hạ du tương đối dồi dào, nhiên chưa có hồ chứa thượng nguồn bổ sung nguồn vùng bị m ặn uy hiếp m ạnh so với tần suất 75% nên việc lấy nước diễn khó khăn Từ phân tích cho thấy: m ực nước, lưu lượng hạ du sông Lam , sông La không bị tác động nhiều từ biến động lưu lượng thượng nguồn, mà chịu tác động m ạnh chế độ thủy triều cửa Hội Đây vùng chịu tác động mạnh xâm nhập mặn, ảnh hưởng lớn đến hoạt động cơng trình thủy lợi - Vùng khơng bị ảnh hưởng thủy triều: Về m ùa kiệt nhu cầu lấy nước sông tăng (qua cống Nam Đàn, cống Đô Lương, Trung Lương, Đức Xá, trạm bơm dọc sông ) làm cho m ực nước sông giảm mạnh có chênh lệch tần suất tính tốn: + So sánh trường hợp trạng tương lai 2020 chưa có hồ điều tiết với tần suất tính tốn hạ lưu Bara Đơ Lương m ực nước tương lai bị hạ thấp 7÷10cm Do phía thượng nguồn có m ột số trạm bơm nhỏ lấy nước tưới, nên lưu lượng đến đập Đơ Lương cịn tương đối lớn từ 77-94m3/s (trường hợp chưa có hồ thượng nguồn điều tiết), đảm bảo lấy nước cấp cho vùng Bắc Nghệ An với nhu cầu 30,4m3/s tương lai 44m 3/s TẠP C H Í KH OA HỌC VÀ CƠN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 KHOA HỌC Tại trạm thủy văn Yên Thượng mực nước đến 2020 bị hạ thấp so với 17÷26cm; lưu lượng Yên Thượng bổ sung thêm nguồn từ sông Giăng, nhiên từ hạ lưu Bara Đô Lương đến Nam Đàn có nhiều trạm bơm lấy nước Lưu lượng Yên Thượng chưa có hồ điều tiết tính cho giai đoạn 2020 cịn 46÷62m3/s (tần suất 75%), 38÷52m3 (P=85%) 34÷47m3/s (P=90%) lưu lượng sông Cả trước lấy vào cống Nam Đàn Qua cho thấy tương lai, lượng nước cần phải lấy vào cống Nam Đàn cấp cho vùng Nam Hưng Nghi 4,04m 3/s Khi chưa có thượng nguồn điều tiết, lưu lượng sông Cả khơng đảm bảo cấp cho vùng này, thiếu khoảng 6÷10m3/s, mặn xâm nhập sâu + Khi có hồ Bản Vẽ bổ sung lưu lượng m ùa kiệt thêm 80m 3/s làm m ực nước vùng không ảnh hưởng triều tăng lên nhiều, lưu lượng bổ sung thêm đảm bảo lượng nước cấp cho cơng trình lấy nước dọc sơng, cụ thể: Mực nước hạ lưu Bara Đô Lương tăng thêm so với chưa có hồ điều tiết 67÷77cm ứngvới tần suất 75%, 72÷81cm ứng với tần suất 85% 75÷84cm ứng với tần suất 90% Lưu lượng tần suất 75%, 85%, 90% thượng lưu Bara Đơ Lương 157÷196m3s thoả mãn lượng nước cấp theo yêu cầu tương lai hệ thống Đô Lương Mực nước m ax cống Nam Đàn 1,44m m in 1,33m tăng thêm so với chưa có hồ điều tiết 80÷107cm (ứng với tần suất 75%); tần suất 85% có mực nước max 1,39m m in 1,26m tăng thêm 145÷147cm; tần suất 90% có m ực nước max 1,32m 1,22m tăng thêm 1,04÷191cm; Lưu lượng thượng lưu cơng trình (n Thượng) từ 113÷143m3/s đảm bảo lượng nước cấp cho cơng trình, có cống Nam Đàn (44,04m3/s) CƠNG NGHỆ Tóm lại: Khi bổ sung nước từ hồ Bản Vẽ mực nước, lưu lượng triền sông Cả tăng cao mùa kiệt, đảm bảo cho việc khai thác cơng trình thủy lợi dọc sông Mặt khác điều tiết làm giảm nồng độ m ặn hạ du sông Cả + Khi có hồ Bản Vẽ bổ sung lưu lượng nước mùa kiệt thêm 80m3/s, Bản Mồng bổ sung 22m 3/s hồ Ngàn Trươi bổ sung 4m 3/s (Phía sơng Cả tăng cường thêm lưu lượng m ùa kiệt thêm 102m3/s, sông La thêm 4m 3/s): Mực nước hạ lưu Bara Đô Lương tăng thêm so với có hồ Bản Vẽ điều tiết với tần suất từ 14÷16cm; Mực nước max cống Nam Đàn m ới tăng thêm so với trường hợp có hồ Bản Vẽ từ 12÷17cm Lưu lượng thượng lưu Bara Đơ Lương từ 167÷196m3/s, Yên Thượng (thượng lưu cống Nam Đàn) từ 135÷165m3/s đảm bảo lưu lượng cho cơng trình lấy nước (hệ thống Đô Lương, cống Nam Đàn trạm bơm dọc sơng) Tóm lại: Khi có bổ sung lưu lượng từ hồ Bản Vẽ, Bản Mồng, Ngàn Trươi m ực nước, lưu lượng dọc sông Cả cải thiện tốt, đảm bảo cao trình m ực nước, lưu lượng cho việc khai thác cơng trình thủy lợi dọc sông VI KẾT LUẬN Lưu vực sông Cả có vị trí địa lý thuận lợi để phát triển kinh tế, sản xuất hàng hóa đa dạng với với ngành m ũi nhọn đặc thù Đây vùng m ạnh cơng nghiệp, nông nghiệp, ngư nghiệp gia công chế biến hàng xuất Cơng tác thủy lợi đóng m ột vai trò quan trọng nhằm khai thác nguồn nước nước m ột cách bền vững, đáp ứng nhu cầu phát triển kinh tế - xã hội địa phương lưu vực Mặc dù nguồn nước sông Cả lớn phân bổ không theo không gian thời gian, mùa kiệt chiếm khoảng 25-30% tổng lượng dòng chảy năm gây khó khăn nhiều cho việc lấy nước cấp cho sản xuất nơng TẠP CH Í KH OA H ỌC VÀ CÔN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 55 KHOA HỌC CÔNG NGHỆ nghiệp ngành kinh tế lưu vực, đặc biệt vùng hạ du bị m ặn xâm nhập sâu, đặc biệt năm cạn kiệt năm 1998, 2002, 2005, 2010…gây cản trở nhiều cho sản xuất nông nghiệp, m ặn xâm nhập vào sâu nội dịa hạn chế việc lấy nước trạm bơm , cửa cống lấy nước Chun đề tính tốn thủy lực kiệt nghiên cứu tính tốn cho trường hợp trạng, tương lai khơng có hồ có hồ điều tiết thượng nguồn với tần suất thiết kế 75%, 85% 90% từ đánh giá, phân tích khó khăn, thuận lợi việc khai thác nguồn nước hạ du sông Cả vị trí lấy nước dọc sơng dịng sơng Cả cống Nam Đàn, trạm bơm…; sông La cống Đức Xá, cống Trung Lương… làm sở để đề xuất giải pháp phù hợp với vùng nhằm đáp ứng phát triển kinh tế xã hội vùng hạ du khai thác bền vững nguồn nước lưu vực sông Cả TÀI LIỆU THAM KHẢO [1] Chuyên đề “Đánh giá tác động dịng chảy kiệt đến vùng hạ lưu sơng Cả", thuộc đề tài Nghiên cứu đề xuất giải pháp giảm thiếu ảnh hưởng dòng chảy kiệt phục vụ sản xuất nông nghiệp thủy sản vùng hạ du sông Cả sông Mã 2013 [2] Chuyên đề thủy lực dự án “Quy hoạch tổng thể thủy lợi khu vực Miền Trung điều kiện Biến đổi khí hậu - nước biển dâng", Viện Quy hoạch Thủy lợi 2012 [3] Chuyên Thủy lực dự án “Rà soát quy hoạch sử dụng tổng hợp nguồn nước lưu vực sông Cả", Viện Quy hoạch Thủy lợi, 2012 [4] Bảng đặc trưng hình thái lưu vực sông Việt Nam [5] Báo cáo Tổng kết tình hình hạn hán, xâm nhập m ặn năm 2010 Sở Nông nghiệp Phát triển Nông thôn tỉnh Nghệ An, Hà Tĩnh [6] Báo cáo trạng cơng trình khai thác nguồn nước hệ thống sông Cả Sở Nông nghiệp Phát triển Nông thôn tỉnh Nghệ An, Hà Tĩnh 56 TẠP C H Í KH OA HỌC VÀ CÔN G N GH Ệ TH ỦY LỢI SỐ 20 - 2014 ... Thuộc sông Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Cả Sông Hiếu Sông Hiếu Sông Hiếu Sông Hiếu Sông Hiếu Sông Giăng Sông Rào... gồm tồn dịng phụ lưu vùng trung, hạ du lưu vực sông Cả, cụ thể sau: + Sông Ngàn Sâu: Từ trạm thuỷ văn Hoà Duyệt đến ngã ba Linh Cảm + Sông La: Từ Linh Cảm đến nhập lưu vào sông Cả ngã ba Chợ Tràng... phát triển kinh tế xã hội vùng hạ du khai thác bền vững nguồn nước lưu vực sông Cả TÀI LIỆU THAM KHẢO [1] Chuyên đề “Đánh giá tác động dòng chảy kiệt đến vùng hạ lưu sông Cả" , thuộc đề tài Nghiên