Đang tải... (xem toàn văn)
Để thấy rõ tác động của hoạt động khai thác cát tới chế độ dòng chảy, biến hình lòng dẫn của đoạn sông nghiên cứu, nhóm tác giả đã ứng dụng mô hình Mike 21C để mô phỏng chế độ động lực và chuyển động bùn cát của một số kịch bản khai thác cát theo quy hoạch được duyệt của các địa phương, sau đó so sánh kết quả nhận được với phương án hiện trạng.
ĐÁNH GIÁ ẢNH HƯ ỞNG CỦA CÁC KỊCH BẢN KHAI THÁC CÁT ĐẾN DIẾN BIẾN LỊNG DẪN SƠNG TIỀN ĐOẠN HẠ LƯ U CẦU MỸ THUẬN PGS.TS Đinh Công Sản, PGS Lê Mạnh Hùng, KS Nguyễn Ngọc Thành Viện Khoa học Thủy lợi Việt Nam Tóm tắt: Sơng Tiền đoạn hạ lưu cầu Mỹ Thuận đoạn sông phân nhánh, nhánh trái chảy địa phận tỉnh Tiền Giang sau đổ biển Đơng qua cửa Đại cửa Tiểu, nhánh phải chảy qua địa phận tỉnh Vĩnh Long đổ biển Đông theo cửa Cổ Chiên, Cung Hầu Ba lai Đây đoạn sơng có chế độ dòng chảy phức tạp lũ, triều, phân lưu dòng chảy nhánh, đồng thời với xáo trộn dòng chảy tác động người khai thác thường xuyên liên tục đoạn sông như: nuôi trồng thủy sản, khai thác cát, giao thông thủy v.v… Để thấy rõ tác động hoạt động khai thác cát tới chế độ dịng chảy, biến hình lịng dẫn đoạn sơng nghiên cứu, nhóm tác giả ứng dụng mơ hình Mike 21C để mơ chế độ động lực chuyển động bùn cát số kịch khai thác cát theo quy hoạch duyệt địa phương, sau so sánh kết nhận với phương án trạng Summary: Tien river at downstream of My Thuan bridge is a branching river, left brach runs through Tien Giang province before going out to the sea via Tieu and Dai river mouths, right branch runs through Vinh Long province prior to flowing to the ocean via Co Chien, Cung Hau and Balai mouths This part of the river with complicated flow regime due to flood, tide, bifurcation and along with turbulent flow due to the human impacts such as: sand mining, aquaculture activities, navigation and so on…To assess the impact of sand mining to hydrodynamics and morphology of this part of river, the authors used MIKE21C model to simulate hydrodynamics and sediment transport of some sand mining scenarios based on provinces’s approved sand mining plans Subsequently the results were compared to the existing scenario I ĐẶT VẤN ĐỀ Sông Tiền đoạn hạ lưu cầu Mỹ Thuận, đoạn sông phân nhánh có chế độ dịng chảy phức tạp, chịu ảnh hưởng khơng chế độ dịng chảy thượng nguồn, chế độ thủy triều biển Đơng, hình thái sơng phân nhánh, mà cịn chịu tác động khơng nhỏ từ hoạt động khai thác dịng sơng người như: giao thơng thủy, nuôi trồng thủy sản, khai thác cát v.v… Tác động nhiều yếu tố khác vào đoạn sông gây nên diễn biến phức tạp, khó kiểm sốt xói lở bờ nhánh phải chảy qua thành phố Vĩnh Long, bồi lắng lòng dẫn nhánh trái chảy qua địa phận tỉnh Tiền Giang, tỷ lệ phân lưu dịng chảy hai nhánh ln thay đổi Hệ lụy việc nhiều nhà cửa, ruộng vườn, sở hạ tầng bị dịng nước trơi nhiều năm nhánh phải Ngược lại lịng sơng nhánh trái bị bồi lắng gây ảnh hưởng không nhỏ tới giao thơng thủy lũ sơng Tiền Để hạn chế diễn biến bất lợi xảy cho đoạn sông này, cần phải xem xét cách đầy đủ ảnh hưởng yếu tố khách quan chủ quan tác động tới chế độ thủy động lực bùn cát đoạn sông Trong khuôn khổ báo xin trình bày phương pháp tính tốn, đánh giá thay đổi chế độ động lực, bùn cát đoạn sông nghiên cứu tác động hoạt động khai thác cát, mơ hình tốn Mike 21C II PHƯ ƠNG PHÁP NGHIÊN CỨU Phương pháp nghiên cứu Đánh giá tác động hoạt động khai thác cát, sở đề xuất giải pháp khai thác cát hợp lý cho đoạn sơng, sử dụng mơ hình vật lý, mơ hình tốn, đo đạc theo dõi thực tế chi tiết thời gian khai thác theo kịch khác kết hợp phương pháp với Trong khuôn khổ báo chúng tơi sử dụng phương pháp mơ hình tốn Mike 21C để nghiên cứu, có tiến hành đo đạc thực tế số liệu thủy văn, dòng chảy bùn cát đoạn sông nghiên cứu để hiệu chỉnh kiểm định mơ hình Mơ hình tốn Mike 21C Viện kỹ thuật tài nguyên nước môi trường Đan M ạch (DHI) thiết lập, ứng dụng phổ biến nhiều nước giới, để nghiên cứu chế độ động lực diễn biến lòng sơng Nhìn chung kết nghiên cứu mơ hình tốn Mike 21C cho kết xác, chấp nhận [0] Thiết lập mơ hình phục vụ nghiên cứu Mô đoạn sông nghiên cứu mơ hình tốn M IKE 21C hệ lưới cong (Curvilinear Grid Generator) vng góc với phương dịng chảy vị trí (thể hình 1), dựa tài địa hình đo đạc vào năm 2010 Viện Khoa học thủy lợi miền Nam thực Hình Phạm vi nghiên cứu sông Tiền đoạn hạ lưu cầu Mỹ Thuận hệ lưới cong Dữ liệu thủy văn, gồm ba biên hở, biên thượng lưu trạm M ỹ Thuận, sử dụng đường trình lưu lượng, hai biên mực nước phía hạ lưu khu vực phà Đình Khao (nhánh Cổ Chiên) Đồng Phú (nhánh sơng Tiền) dùng đường q trình mực nước Biên lưu lượng hai biên mực nước, ba biên hở cung cấp cho mơ hình M IK E 21C trích xuất từ kết tính tốn từ mơ hình MIKE11 cho tồn hệ thống sơng Mê Kơng, sau kiểm định độ xác Dữ liệu bùn cát [0], lấy từ mơ hình M IKE 11 chạy cho tồn hệ thống đồng sơng Cửu Long, đường kính hạt cát đáy, cấp phối hạt xác định từ mẫu bùn cát đáy Viện Khoa học Thủy lợi miền Nam lấy phân tích mẫu vào tháng 12/2010 (lấy d=0,15mm cho khu vực tính tốn) Khi tính tốn biến hình lịng dẫn đoạn sơng nghiên cứu M ike 21C, lựa chọn cơng thức Engelund and Fredsore [0,2] để tính vận chuyển bùn cát đáy công thức Lane-Kalinske để tính bùn cát lơ lửng Kết cho thấy việc lựa chọn cơng thức để tính tốn hồn toàn phù hợp với số liệu thực đo Kết hiệu chỉnh kiểm định mơ hình Hiệu chỉnh kiểm định mơ hình tốn quan trọng từ việc làm đảm bảo độ xác kết mơ tính tốn dự báo theo trường hợp, theo kịch quan tâm Để tiến hành hiệu chỉnh kiểm định mô hình, cần có số liệu thực đo chi tiết với độ xác cao, để làm sở cho việc hiệu chỉnh cơng thức sử dụng, thơng số tính tốn … cho kết tính tốn nhận sai khác so với số liệu thực đo Điều thể qua trị số hiệu Nash-Sutcliffe (NSE) tỷ số độ lệch quan trắc tiêu chuẩn (RSR) Chỉ số NSE thông số xác định phương sai dư so với phương sai chuỗi số liệu thực đo, tính cơng thức [0]: NSE = 1- (1) Trong đó: - NSE hệ số Nash-Sutcliffe; giá trị thực đo thứ i; - giá trị mơ thứ i mơ hình; Giá trị thực đo trung bình; - n số lần thực đo Trị số RSR dùng để tính độ lệch cho phép kết từ mơ hình với số liệu thực đo Giá trị tối ưu RSR 0, với giá trị thấp biểu thị mô xác cao RSR tính theo cơng thức: RSR = = (2) Trong tham số cơng (2) nêu công thức (1) Căn để đánh giá mức độ xác mơ hình thực tế theo hai số trình bày Bảng Bảng Đánh giá độ xác mơ hình theo số RSR, NSE [0] Đánh giá độ xác mơ hình Rất tốt Tốt Đạt Không đạt RSR NSE 0,00 RSR 0,50 0,50 < RSR 0,60 0,60 < RSR 0,70 RSR > 0,70 0,75 < NSE 1,00 0,65 < NSE 0,75 0,50 < NSE 0,65 NSE 0,50 Để làm sở cho việc hiệu chỉnh kiểm định mơ hình, chúng tơi tiến hành mơ q trình động lực biến đổi lịng dẫn đoạn sơng nghiên cứu khoảng thời gian từ tháng năm 2009 đến tháng 12 năm 2011, sau trích xuất số liệu để so sánh kết với số liệu thực đo ngày 22, 23 24/12/2010 Hình thể phù hợp kết tính tốn với giá trị lưu lượng thực đo hai nhánh (a) (b) Hình Kết kiểm định phân chia lưu lượng nhánh trái (a) nhánh phải (b) Hình thể phân bố vận tốc mặt trị số vận tốc mặt cắt trích xuất từ mơ hình Mike 21C thực tế đo đạc thiết bị đo vận tốc ADCP, lúc ngày 22/12/2010. Hình So sánh phân bố lưu tốc mặt cắt MT-2 MT-3 lúc 9h00' ngày 22-12-2010 kết mô MIKE 21C tài liệu thực đo thiết bị ADCP Hình thể trình diễn biến lòng dẫn hai mặt cắt ngang nhánh phải Trường An Tp Vĩnh Long Hình Diễn biến hình thái mặt cắt nhánh phải Trường An (a) Tp.Vĩnh Long (b) Số liệu tính tốn cuối mơ hình tốn so sánh với giá trị thực đo, sau tính trị số hiệu Nash-Sutcliffe (NSE) theo công thức (1) tỷ số độ lệch quan trắc tiêu chuẩn (RSR), theo công thức Kết nhận NSE = 0,975 RSR = 0,156 Từ kết kiểm định cho t hấy yếu tố thủy lực biến hình lịng dẫn mơ hình so với thực đo có sai khác khơng nhiều Chính dùng mơ hình với tham số hiệu chỉnh để nghiên cứu tác động kịch khai thác cát III KẾT Q UẢ NGHIÊN CỨU VÀ THẢO LUẬN Q ui hoạch khai thác cát kịch 3.1 Quy hoạch khai thác cát phê duyệt M ỏ cát khu vực nghiên cứu thuộc địa phận tỉnh Vĩnh Long Tiền Giang Theo quy hoạch khai thác cát sông Tiền tỉnh Vĩnh Long Tiền Giang đến năm 2020, khu vực nghiên cứu có thân cát ST-TQ.1 (thuộc Vĩnh Long) thân cát II (thuộc Tiền Giang) với quy mô: chiều dài 12,4 km, chiều rộng trung bình 600 m, chiều sâu thân cát 3,85 m Vĩnh Long 5,5 m Tiền Giang [0,0] Hình Bản đồ Quy hoạch khai thác cát khu vực hạ lưu cầu Mỹ Thuận 3.2 Các kịch khai thác cát Các kịch khai thác cát xây dựng dựa yêu cầu thực tế địa phương quản lý, cấp phép kỹ thuật khai thác Trong phạm vi nghiên cứu khu vực khai thác chia thành vùng khai thác kí hiệu KV1, KV2 KV3 tương ứng, hình Hình Vị trí khu vực khai thác cát đoạn sông nghiên cứu vị trí theo dõi biến đổi lưu tốc dịng chảy hai nhánh sơng Để thấy rõ ảnh hưởng kịch khai thác cát đoạn sông đến tỷ lệ phân lưu lưu lượng hai nhánh thay đổi độ lớn vận tốc dòng chảy số điểm đặc trưng (các vị trí VT1, VT2, VT3 nhánh phải vị trí VT4, VT5 VT6 nhánh trái), điểm nhạy cảm với xói lở bồi lắng lịng dẫn Trong báo tiến hành xem xét kịch khai thác cát khác nhau: - Kịch 1, khai thác cát khu vực KV1 với quy mô khai thác chiều rộng B=600m chiều dài L=5km độ sâu H= -15m, nhánh trái (xem hình 6); - Kịch 2, khai thác cát khu vực KV2 với quy mô khai thác chiều rộng B=500m chiều dài L=3km độ sâu H= -18m nhánh phải (hình 6); - Kịch 3, khai thác đồng thời hai khu vực KV1và KV3 (khu vực KV3 khai thác với quy mô B=750m, L=1,4km H=-18m) - Kịch 4, khai thác cát đồng thời khu vực KV1, KV2 KV3 với quy mô kịch KB1, KB2 KB3 Kết tính tốn thảo luận 3.3 Ảnh hưởng kịch khai thác cát tới tỷ lệ phân lưu dòng chảy hai nhánh Xem xét khả biến đổi phân chia lưu lượng có lợi kịch KTC KV1 KV3 xét mức độ giảm thiểu xói lở nhánh sơng Cổ Chiên, giảm lưu lượng khoảng 5% cho nhánh (ứng với cấp lưu lượng tạo lòng) Đồng thời, lưu lượng nhánh sông Tiền gia tăng đạt tới 5% tương ứng, làm giảm thiểu bồi lắng Kết cụ thể xem bảng Bảng Tỷ lệ phân lưu trung bình hai nhánh sơng Cổ Chiên Sông Tiền theo kịch khai thác cát ứng với cấp lưu lượng thượng lưu Cấp Q (m3/s) 14,000 ÷20,500 Tăng giảm so với khơng khai thác cát 9,000 ÷13,900 Tăng giảm so với khơng khai thác cát 4,000÷ 8,900 Tăng giảm so với không khai thác cát Hiện trạng 2009 2010 s Cổ Sông s Cổ Sông s Cổ Sông Chiên Tiền Chiên Tiền Chiên Tiền 69.56 30.44 65.41 34.59 64.61 35.39 -4.16 4.16 -4.96 4.96 2011 s Cổ Sông Chiên Tiền 64.60 35.40 -4.96 4.96 70.83 29.17 67.29 -3.54 32.71 3.54 65.23 34.77 -5.60 5.60 65.87 34.13 -4.96 4.96 71.86 28.14 67.96 -3.90 32.04 3.90 67.79 32.21 -4.07 4.07 66.67 33.33 -5.19 5.19 3.4 Ảnh hưởng kịch khai thác cát tới độ lớn vận tốc dòng chảy vị trí đặc trưng Xem xét đường q trình lưu tốc vị trí VT1,VT2,VT3 nhánh phải (sơng Cổ Chiên) VT4,VT5, VT6 nhánh trái (sông Tiền) có nhận xét: - Nhìn chung, tác động kịch khai thác cát lưu tốc dòng chảy không lớn, sông rộng sâu, mặt cắt thay đổi khai thác cát không nhiều - Vị trí VT1, VT2 VT3 khai thác cát theo KB2 lưu tốc tăng so với trạng, dịng chảy nhánh sơng Cổ Chiên tăng lên Các kịch khác lưu tốc vị trí giảm so với trạng, giảm nhiều KB3 VT3, với mức giảm khoảng 0,05m/s tương đương với khoảng 9% lưu tốc phương án trạng Tác động kịch làm giảm khả xói lở khu vực Tp.Vĩnh Long - Tại vị trí VT4, VT5 VT6 lưu tốc có thay đổi so với trạng tất kịch bản, mức độ không nhiều Hie n tra ng Kha i thac KV Kha i thac KV Kha i thac KV 1+K V3 Kha i thac KV 1+K V2+K V3 [m/s] [m/s] [m/s] [m/s] [m/s] Hien tra ng Kha i thac KV1 Kha i thac KV2 Kha i thac KV1+KV3 Kha i thac KV1+KV2+ KV3 [m/s] [m/s] [m/s] [m/s] [m/s] 0.65 0.55 Van toc (m/s ) Van toc ( m/s ) 0.60 0.55 0.50 0.45 0.50 0.45 0.40 0.40 12:00 2011-08- 21 00:00 08- 22 12:00 00:00 08-23 T hoi gian (gio) 12 :00 00:00 08-24 12: 00 2:00 011-09-20 00:00 09- 21 12:00 Vị trí Hie n tra ng Kha i thac KV Kha i thac KV Kha i thac KV 1+K V3 Kha i thac KV 1+K V2+K V3 0:00 9-22 T hoi gian (g io) 12:00 00:00 09-23 12: 00 Vị trí [m/s] [m/s] [m/s] [m/s] [m/s] H ien tra ng K i thac KV K i thac KV K i thac KV 1+KV K i thac KV 1+KV 2+K V3 [m/s] [m/s] [m/s] [m/s] [m/s] Van toc (m/s ) Van toc ( m/s) 0.60 0.55 0.50 0.45 0.30 0.25 0.20 0.40 12:00 2011-08- 21 00:00 08-22 12:00 00:00 08-23 T hoi gian ( gio) 12:00 00 :00 08 -24 00:00 2011- 08- 23 12: 00 12:00 00: 00 08- 25 12:00 Vị trí Vị trí Hie n tra ng Kha i thac KV Kha i thac KV Kha i thac KV 1+K V3 Kha i thac KV 1+K V2+K V3 00:00 12:00 08-24 Thoi gian (gio) [m/s] [m/s] [m/s] [m/s] [m/s] H ien tra ng [ m/s ] K hai thac KV [ m/s ] K hai thac KV [ m/s ] K hai thac KV [ m/s ] K hai thac KV 1+KV 2+KV [ m/s ] Van toc (m/s ) Van toc ( m/s) 0.180 0.25 0.20 0.160 0.140 0.120 0.100 00:00 2011-09- 19 12:00 00:00 09-20 12:00 00:00 09- 21 T hoi gian (gio) 12 :00 00:00 09-22 12: 00 2011- 09-05 00:00 09-06 12:00 Vị trí 00:00 12:00 09-07 T hoi gian ( gio) 00:00 09-08 12:00 00:00 09-09 Vị trí Hình Đường trình thay đổi lưu tốc kịch khai thác cát so với trạng (không khai thác cát) 3.5 Ảnh hưởng kịch khai thác cát tới biến hình lịng dẫn Kết biến hình lịng dẫn thời gian năm, với liệt thủy văn, bùn cát tương tự năm 2009, 2010 2011 (đại diện cho năm lũ nhỏ, trung bình lũ lớn) theo kịch trạng (không khai thác cát) kịch khai thác cát trình bày hình hình MzRe sultView1 M zRe sultView1 1142 25 500 00 1142500 1142 20 000 00 1142000 1141 15 500 00 1141500 1141 10 000 00 1141000 1140 05 500 00 1140500 1140 00 000 00 1140000 1139 95 500 00 1139500 1139000 1139 90 000 00 1138500 1138 85 500 00 Da ta [-] Abo ve 0- -2.5 -5 - -2 -7.5 -10 - -7 -12 - -15 - -12 -17 - -20 - -17 -22 - -25 - -22 -27 - -30 - -27 -32 - Bel o w -32 Un de fi n ed 1138 80 000 00 1137 75 500 00 1137 70 000 00 1136 65 500 00 1136 60 000 00 1135 55 500 00 1135 50 000 00 1134 45 500 00 1134 40 000 00 60 00 000 00 60 020 00 00 01/01/02 00:00:00, Time step0 of 40 604 000 00 06 60 600 00 60 080 00 00 Bed e l ve l [-] Abo ve - -2.5 -5 - -2 -7.5 -5 -10 - -7 -12 - -1 -15 - -1 2.5 -17 - -1 -20 - -1 7.5 -22 - -2 -25 - -2 2.5 -27 - -2 -30 - -2 7.5 -32 - -3 Bel o w -3 2.5 Un de fi n ed Va 1138000 1137500 1137000 1136500 1136000 1135500 1135000 1134500 1134000 600000 00 610 000 00 602000 604000 606000 608000 610000 12/3 1/10 15:40:00, Time s e t p 1652 of 1652 Hình Địa hình trạng địa hình sau năm khơng khai thác cát Mz Res ultVie w6 MzResultV iew4 11 425 00 11 11 25 25 00 00 11 420 00 11 11 20 20 00 00 11 415 00 11 11 15 15 00 00 11 410 00 11 11 10 10 00 00 11 405 00 11 11 05 05 00 00 11 400 00 11 11 00 00 00 00 11 395 00 11 11 95 95 00 00 11 390 00 11 11 90 90 00 00 11 385 00 11 11 85 85 00 00 11 380 00 B ed l ev el [ -] 11 11 80 80 00 00 Ab ov e - -2 -5 - -2 -7 -5 -10 - -7 -12 - -10 -15 - -12 -17 - -15 -20 - -17 -22 - -20 B edl evel [ - ] 11 375 00 A bove 0- -2 -5 - - -7 - - -1 0- - -1 2.5 - - 10 -1 5- - 12 -1 7.5 - - 15 -2 0- - 17 -2 2.5 - - 20 -2 5- - 22 -2 7.5 - - 25 -3 0- - 27 -3 2.5 - - 30 Bel ow - 32 U ndef ined V alu e 11 370 00 11 365 00 11 360 00 11 355 00 11 350 00 11 345 00 11 340 00 0000 60 2000 6040 00 0600 60 8000 11 11 75 75 00 00 11 11 70 70 00 00 11 11 65 65 00 00 11 11 60 60 00 00 11 11 55 55 00 00 -25 - -22 -27 - -25 -30 - -27 -32 - -30 Be o l w -32 Unde f n i ed V al ue 11 11 50 50 00 00 11 11 45 45 00 00 11 11 40 40 00 00 610 000 60 60 00 000 60 02 200 00 0 6/0 6/1 1:0 0:0 0, Tim e s te p of 22 10/ 30/ 11 3: 40: 00, Ti mest ep 42 of 42 60 60 40 000 60 06 600 00 60 60 80 000 61 10 000 00 Hình Địa hình sau năm sau khai thác cát KV1 +KV3 (bên trái) khai thác KV1+KV2+KV3 (bên phải) Để xem xét biến hình lịng dẫn sau năm mơ cách rõ ràng hơn, xem xét số mặt cắt ngang khu vực nghiên cứu, thể hình 10 Kết biến hình lịng dẫn mặt cắt hình 11 đến hình 15 Hình 10 Vị trí trích xuất mặt cắt ngang lịng dẫn để đánh giá diễn biến sau năm Hình 11 Biến hình mặt cắt MC1 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 12 Biến hình mặt cắt MC2 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 13 Biến hình mặt cắt MC3 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 14 Biến hình mặt cắt MC4 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 15 Biến hình mặt cắt MC5 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Qua kết tính tốn cho thấy: - Khi khơng khai thác cát, biến hình lịng dẫn sau năm nhánh phải (sơng Cổ Chiên) có xu bị xói, nhánh trái (sơng Tiền) có xu bồi tụ mạnh đặc biệt khu vực ngang với đầu cù lao An Bình Vấn đề khơng có lợi cho việc phịng chống xói lở bồi lắng khu vực - Khai thác cát khu vực quy hoạch khơng ảnh hưởng nhiều đến xói lở bồi lắng lịng sơng Xu biến đổi lịng dẫn có lợi có khai thác cát, M C1, M C2 MC3 - Khai thác cát theo kịch (KB3) cho kết có lợi biến hình lịng dẫn cho đoạn sơng phương diện chống xói lở, bồi lắng: nhánh phải (sông Cổ Chiên - Tp.Vĩnh Long) giảm xói lở nhánh trái (sơng Tiền) giảm bồi lắng, giảm ách tắc giao thông thủy IV KẾT LUẬN VÀ KIẾN NGHỊ Khai thác cát việc làm cần thiết, giá trị cao nguồn tài nguyên cát đem lại, khai thác vị trí, kỹ thuật góp phần điều khiển dịng chảy, giảm thiểu xói bồi Ngược lại gây tác động xấu phương diện gia tăng xói bồi khu vực khơng mong muốn Từ kết nghiên cứu khu vực sông Tiền hạ lưu cầu Mỹ Thuận, diễn biến vận tốc dòng chảy, thay đổi tỷ lệ phân lưu diễn biến xói bồi lịng dẫn sông Tiền (đoạn Tiền Giang) sông Cổ Chiên (đoạn thị xã Vĩnh Long) đưa số ý kiến sau: - Khai thác cát theo kịch KB1, KB3 BK4 giúp khai thơng dịng chảy bên nhánh sông Tiền, làm cho lưu lượng, lưu tốc nhánh sơng Cổ Chiên giảm Xói bồi lịng sơng khu vực nghiên cứu phân bố lại theo xu có lợi hai nhánh sơng Trong kịch kịch KB3 có hiệu (tức không khai thác cát nhánh sông Cổ Chiên) - Khai thác cát theo kịch KB2 cho kết ngược lại với kịch trên, làm cho lưu lượng, lưu tốc nhánh sông Cổ Chiên tăng lên, mặt khác khu vực khai thác cát theo kịch KB2 có kè Cổ Chiên thi cơng hồn thành giai đoạn tới, khơng có lợi chống sạt lở Chúng tơi kiến nghị hạn chế dần tiến tới dừng khai thác khu vực Tình hình khai thác cát thực tế khu vực diễn phức tạp, ngồi nạn khai thác cát lậu khai thác cát khơng tn thủ theo vị trí quy hoạch, thời gian quy trình khai thác diễn phổ biến địa phương Kết việc gây sạt lở bờ sông khu vực lân cận Chúng tơi kiến nghị quyền địa phương phải theo dõi khắt khe trình khai thác cát doanh nghiệp để đảm bảo khai thác cát theo quy hoạch duyệt, để tránh gia tăng xói lở khu vực khơng làm thất tài ngun TÀI LIỆU THAM KHẢO [1] Đinh Công Sản (2007), “Một số vấn đề động lực học dịng chảy hình thái sơng Cửu Long”, Luận án Tiến sĩ kỹ thuật, Viện Khoa học Thủy lợi miền Nam 10 [2] Dinh Cong San, (2005) “Research on river bed erosion and sedimentation prediction by MIKE21C model at Tan Chau-Hong Ngu area, in the Mekong River”, Proceedings of the International Symposium on Sustainable Development in the Mekong River basin, pp 188-195 [3] Dinh Cong San, (2008) “River bed scour and deposition, the causes and suitable measure to stabilize river bank and protection work at Lasan – M aithon reach, Thanh Da peninsula – Sai Gon River – Hochiminh City”, 1st International Conference for Environment and Natural Resources March 17th - 18th, 2008, Environmental Protection for Urban and Industrial Zones to International Integration, HCM C.Vietnam [4] Lê Mạnh Hùng (2004), Báo cáo tổng kết dự án Khoa học công nghệ cấp nhà nước, “Nghiên cứu dự báo xói lở, bồi lắng lịng dẫn đề xuất biện pháp phịng chống cho hệ thống sơng Đồng sông Cửu Long”, Viện KHTL miền Nam [5] Lê Mạnh Hùng nnk (2010), Báo cáo khảo sát thủy văn, “Nghiên cứu ảnh hưởng hoạt động khai thác cát đến thay đổi lịng dẫn sơng cửu long (sơng Tiền, sông Hậu) đề suất giải pháp quản lý, quy hoạch khai thác hợp lý” – Viện Khoa học Thủy lợi miền Nam [6] Moriasi D.N et al (2007), Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations Transactions of the ASABE 50, 885 [7] Sở TN&M T tỉnh Tiền Giang (2006), Quy hoạch thăm dò, khai thác sử dụng tài ngun khống sản cát lịng sông Tiền, tỉnh Tiền Giang đến năm 2020 [8] UBND Tỉnh Vĩnh Long 29/12/2009 Quyết định việc phê duyệt Quy hoạch thăm dị, khai thác, sử dụng khống sản cát lịng sơng địa bàn tỉnh Vĩnh Long đến năm 2020 số 29/2009/QĐ-UBND 11 ... Quy hoạch khai thác cát khu vực hạ lưu cầu Mỹ Thuận 3.2 Các kịch khai thác cát Các kịch khai thác cát xây dựng dựa yêu cầu thực tế địa phương quản lý, cấp phép kỹ thuật khai thác Trong phạm vi nghiên... 2011 theo kịch khai thác cát khác Hình 12 Biến hình mặt cắt MC2 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 13 Biến hình mặt cắt MC3 hạ lưu cầu Mỹ Thuận. .. 2011 theo kịch khai thác cát khác Hình 14 Biến hình mặt cắt MC4 hạ lưu cầu Mỹ Thuận sau năm thủy văn 2009-2010 2011 theo kịch khai thác cát khác Hình 15 Biến hình mặt cắt MC5 hạ lưu cầu Mỹ Thuận
