Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 ỨNG DỤNG MÔ HÌNH THỦY LỰC HAI CHIỀU MÔ PHỎNG ĐẶC TÍNH THỦY LỰC VÀ TÍNH TOÁN BỒI XÓI Ở VÙNG CỬA SÔNG ĐỊNH AN Nguyễn Phương Tân1, Văn Phạm Đăng Trí2 Võ Quốc Thành1 Khoa Môi trường & Tài nguyên Thiên nhiên, Trường Đại học Cần Thơ Thông tin chung: Ngày nhận: 22/10/2013 Ngày chấp nhận: 28/04/2014 Title: The application of a twodimensional hydraulic model to simulate the hydraulic properties and deposition and erosion patterns in the Dinh An estuary Từ khóa: Bồi lắng, xói lở, cửa sông, mô hình thủy lực hai chiều (CCHE2D), Đồng sông Cửu Long Keywords: Deposition, erosion, river mouth, two-dimensional hydraulic model (CCHE2D), and Vietnamese Mekong Delta ABSTRACT In the recent years, estuary deposition and erosion due to hydrodynamics changes are amongst the major problems of the Dinh An estuary, one of the two river mouths of the Hau River in the Vietnamese Mekong Delta Such natural phenomena are projected to be even greater in terms of magnitude and unpredictable in terms of time and space due to (i) the construction and operation of hydropower dams in the upstream section of the Mekong River, leading to changes of flow regime and sedimentation patterns along the river; (ii) on-going and planned hydraulic construction projects to meet different water requirements for agricultural land uses along the Mekong Basin, especially in the delta; and (iii) tidal regime changes in the East Sea In this study, a two-dimensional hydraulic model (CCHE2D) is used to study the hydrodynamics, deposition and erosion patterns in the Dinh An estuary The hydraulic component of the model is calibrated and validated based on the flows measured in August 2012 The sediment transport component is applied based on the calculated hydraulic properties and referenced sediment data (including suspended sediment concentration and transport sediment rate) The results obtained from this study set a stage for subsequent studies to understand the morphological changes of the river mouths in the Mekong Delta TÓM TẮT Trong năm gần đây, bồi lắng xói lở vấn đề cần quan tâm cửa sông Định An, hai cửa sông sông Hậu Đồng sông Cửu Long Quá trình tự nhiên dự báo nghiêm trọng khó dự đoán tương lai do: (i) trình xây dựng vận hành đập thủy điện thượng nguồn sông Mekong, dẫn đến thay đổi chế độ dòng chảy vận chuyển bùn cát dọc dòng sông; (ii) công trình thủy lợi xây dựng quy hoạch để phục vụ sản xuất nông nghiệp toàn lưu vực sông Mekong, đặc biệt vùng đồng bằng; và, (iii) thay đổi động thái thủy triều biển Đông Trong nghiên cứu này, mô hình thủy lực hai chiều (CCHE2D) sử dụng để xác định động thái dòng chảy, bồi lắng xói lở cửa sông Định An Thành phần thủy lực mô hình hiệu chỉnh kiểm định dựa số liệu đo đạc vào tháng năm 2012; thành phần vận chuyển bùn cát mô hình tính toán dựa vào đặc tính thủy lực số liệu bùn cát tham khảo từ báo (bao gồm: nồng độ bùn cát lơ lửng vận chuyển bùn cát đáy sông) Kết thu từ nghiên cứu tạo tiền đề cho nghiên cứu để tính toán thay đổi địa mạo đáy sông cho hệ thống cửa sông Đồng sông Cửu Long Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 vệ tinh từ năm 1989 – 2001, cửa Định An xói lở khoảng 4.5 m/năm bồi lắng khoảng 9.5 m/năm (Hình 1B) (Nguyễn Việt Thanh et al., 2011) Trong tương lai gần, gia tăng hoạt động phát triển nước thượng nguồn (ví dụ: việc gia tăng lấy nước phục vụ cho nông nghiệp vùng Đông Bắc Thái Lan, dự kiến xây dựng đập dòng sông Mekong thuộc Lào Campuchia) dự báo làm thay đổi đáng kể chế độ dòng chảy hệ thống nhánh sông ĐBSCL Bên cạnh đó, ảnh hưởng biến đổi khí hậu toàn cầu nên nước biển dâng lên thêm khoảng 2327 cm vào năm 2050 theo kịch trung bình (B2) (Bảng 1) (Bộ Tài nguyên Môi trường, 2012) Chính thế, việc tìm hiểu động thái dòng chảy xác định xu hướng bồi lắng, xói lở quan trọng để đánh giá tính ổn định hệ thống sông góp phần giải thích động thái nguồn tài nguyên nước mặt hệ thống dòng sông (dưới tác động thay đổi dòng chảy thượng nguồn động thái thủy triều) nhằm đảm bảo phát triển bền vững vùng GIỚI THIỆU Đồng sông Cửu Long (ĐBSCL) (Hình 1A) vùng hạ lưu cuối sông Mekong trước đổ biển Đông ĐBSCL xem vùng canh tác nông nghiệp nuôi trồng thủy sản lớn Việt Nam khu vực quan trọng đóng góp (hơn 50% sản lượng lương thực khoảng 70% sản lượng thủy sản) cho phát triển kinh tế xã hội chung Việt Nam (Trịnh Công Văn, 2009) Sông Hậu hai phân lưu sông Mekong Trước đổ biển Đông, sông Hậu chia thành hai nhánh: Định An Trần Đề (Hình 1B) Cửa sông Định An hai cửa sông giữ vai trò quan trọng phát triển kinh tế - xã hội vùng bán đảo Cà Mau (Hoa Mạnh Hùng et al., 2008) Tuy nhiên, năm gần vùng cửa sông Định An thường xuyên xảy trình bồi lắng, xói lở, xâm nhập mặn (Bùi Hồng Long & Tống Phước Hoàng Sơn, 2003; Hoa Mạnh Hùng et al., 2008) Theo kết nghiên cứu từ ảnh B A Bồi lắng Xói lở Hình 1: Vùng nghiên cứu - đoạn cửa sông Định An (cách cửa biển khoảng km) Bảng 1: Mực nước biển dâng (cm) kịch (Bộ Tài nguyên Môi trường, 2012) Các mốc thời gian kỷ 21 2040 2050 2060 2070 2080 2090 2100 Thấp (B1) 17-19 22-26 28-34 34-42 40-50 46-59 51-66 Trung bình (B2) 17-20 23-27 30-35 37-44 44-54 51-64 59-75 Cao (A1FI) 19-21 26-30 35-41 45-53 56-68 68-83 85-105 dừng lại mức độ tính toán giá trị thủy lực trung Ngày nay, với tiến khoa học máy tính bình mặt cắt đưa vào tính toán với phát triển kỹ thuật tính toán mô hình không xét đến dòng chảy ngang đại, mô hình toán thủy lực chiều (1D) mô hình thủy lực hai chiều (2D) ba xây dựng ứng dụng nhiều việc chiều (3D) chưa phản ánh đầy đủ chất tìm hiểu động thái thủy lực dòng chảy ĐBSCL động thái thủy lực dòng chảy Hơn nữa, hệ số (Trần Quốc Đạt et al., 2012; Văn Phạm Đăng Trí et nhám thủy lực áp dụng mô hình thủy al., 2013) Các mô hình thủy lực 1D lực 1D (thường) không thay đổi mặt cắt Kịch 2020 8-9 8-9 8-9 2030 11-13 12-14 13-14 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 vận chuyển bùn cát dự báo thay đổi hình thái sông như: CCHE2D, Delft2D-River, MIKE21C, TABS-MD TELEMAC (hoặc đoạn sông mô phỏng) thực tế hệ số nhám thủy lực thay đổi phân đoạn mặt cắt ngang dọc (theo chiều dòng chảy) (Nguyễn Cảnh Cầm et al., 2007; Wilson & Atkinson, 2007) Do đó, việc xác định động thái thủy lực dòng chảy thay đổi hình thái lòng sông trình phức tạp hệ thống sông/kênh tự nhiên cần có hỗ trợ mô hình toán thích hợp (Wu, 2007) PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Khảo sát thực địa Thiết bị đo đạc đáy sông kỹ thuật số ADCP (Acoustic Doppler Current Profiler), thiết bị cảm ứng (sensor) hệ thống định vị toàn cầu (Global Position System - GPS) sử dụng chuyến thực địa từ ngày 23/8/2012 đến ngày 31/8/2012 Thiết bị ADCP dùng để đo đạc cắt lớp dòng chảy sóng hồi âm theo hiệu ứng Doppler Thiết bị gắn phía bên thân tàu (Hình 2A) kết nối với máy tính để ghi nhận lại giá trị lưu lượng cao trình cột nước mặt cắt, vận tốc dòng chảy Thiết bị cảm ứng (sensor) dùng để đo giá trị mực nước đầu vào giá trị mực nước dùng để kiểm định mô hình GPS dùng để xác định vị trí điểm đo đạc điểm xói lở bờ (Hình 2B) Dữ liệu đo đạc từ ADCP xử lý từ phần mềm Winriver II với kết có bao gồm lưu lượng, độ sâu, vận tốc dòng chảy cho mặt cắt (Teledyne RD Instrument, 2007) Mô hình thủy lực 2D 3D mô chi tiết thay đổi động thái thủy lực dòng chảy đánh giá thay đổi hình thái đoạn sông (Wu, 2007) Tuy nhiên, hạn chế số liệu thu thập theo chiều thẳng đứng (từ mặt nước xuống đáy sông), mô hình thủy lực 2D thường chấp nhận nghiên cứu (thay sử dụng mô hình thủy lực 3D) (Jia & Wang, 1999) Do vậy, việc ứng dụng mô hình thủy lực 2D nhằm xác định động thái thủy lực dòng chảy, xu hướng bồi lắng xói lở vùng cửa sông ĐBSCL cần thiết, đặc biệt bối cảnh có thay đổi đáng kể điều kiện thủy văn tự nhiên, tình hình phát triển thượng nguồn lưu vực sông Mekong nước biển dâng Hiện nay, số mô hình thủy lực có khả tính toán mô Hình 2: Tàu khảo sát thiết bị đo đạc (A); xác định vị trí xói lở gần cửa sông (B) 2.2 Mô hình thuỷ lực chiều CCHE2D Mô hình thủy lực CCHE2D (có khả mô đặc tính thủy động lực học hai chiều (Langendoen, 2001)) sử dụng nghiên cứu Mô hình trường đại học công nghệ Mississippi phát triển nhằm mô trình truyền thuỷ lực, chất lượng nước, vận chuyển bùn cát biến động lòng dẫn Mô hình gồm hợp phần (module) sau (Hình 3): Hình 3: Thành phần mô hình CCHE2D 10 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 Mô hình tạo lưới (CCHE2D Mesh Generator): tạo lưới giá trị hệ tọa độ cong hình số liệu đo đạc Cao độ mực nước tự tính toán phương trình liên tục (PT1): Mô hình số (CCHE2D Numerical model): tính toán trình thủy lực Mô hình giao diện đồ họa (CCHE2D- GUI Graphical Users Intreface): nhập thông số mô (PT1) Phương trình động lượng theo phương x (Hình 4) (PT2): (PT2) Phương trình động lượng theo phương y (Hình 4) (PT3): (PT3) lớp bùn cát đáy phía dưới; đó, dòng bùn cát chiếm toàn độ sâu dòng chảy h dòng bùn cát nằm lớp đáy có bề dày nhỏ Trong : cao độ mực nước (m); t: thời gian (s); h: độ sâu cột nước (m); u,v : vận tốc trung bình theo hai phương x y (m/s); fr: hệ số ma sát đáy; , , , : ứng fCor: thông số Corilolis; suất Reynolds (N/m2); , : ứng suất tiếp tuyến đáy (N/m2) Hình 4: Hướng dòng chảy theo hai phương x, y Mô hình CCHE2D tính toán vận chuyển bùn cát hai lớp: lớp bùn cát phía (lớp mặt) Phương trình truyền tải bùn cát lơ lửng có dạng sau (PT5): (PT5) Phương trình liên tục vật chất đáy (PT6): (PT6) Trong đó: : nồng độ bùn cát (kg/m3); u,v: vận tốc trung bình theo hai phương x y (m/s); : hệ số khuếch tán rối bùn cát (m2/s); : bề dày lớp : lượng trầm tích bồi lắng; bùn cát đáy (m); : lượng trầm tích xói lở; : độ rỗng (%); : nồng độ trung bình bùn cát đáy (kg/m3); và, , : thành phần chuyển tải bùn cát theo phương x, y (kg/s/m) (Wu, 2001) 2.3 Phương pháp giải phương trình mô hình CCHE2D Các phương trình giải phương pháp phần tử hữu hiệu (Efficient Element Method); phương pháp phần tử hữu hạn đặc biệt, hệ phương trình giải hệ tọa độ cong ( , ), ứng với lưới cong tính toán miền thực A B Hình 5: Lưới cong xen kẽ tính toán u, v, z (A); vị trí tính toán thành phần vận tốc mực nước lưới tính (B) 11 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 phần hạt bùn cát đáy sông xác định Bảng (Albers & Lieberman, 2011; Nguyễn Việt Thanh et al., 2011; Wolanski et al., 1996); Trong lưới cong miền tính, thành phần vận tốc giải từ phương trình động lượng phương pháp sai phân hữu hạn sơ đồ ẩn nút theo không gian mực nước giải từ phương trình liên tục phần tử liên kết điểm (Hình A) Trong lưới cong tính toán, thành phần vận tốc bố trí nút lưới mực nước bố trí tâm ô lưới (Hình B) 2.4 Thành lập mô hình Lưới tính: lưới tính tạo thành bao gồm 9.000 nút (node) phân bố dọc miền lưới tính (bao gồm: 30 ô dọc theo hướng dòng chảy 300 ô vuông góc với hướng dòng chảy) Thành phần mẫu bùn cát độ rỗng xác định Bảng (Rijn, 1993; Zou Xue et al., 2010) Bước thời gian: Bước thời gian chọn tính toán 60 giây, thiết lập thời gian xuất kết điểm quan sát để so sánh với kết thực đo điểm quan sát Bảng 2: Thành phần bùn cát lơ lửng Kích Thước (mm) Thành phần phần trăm (%) 0.00100 25 0.00425 75 0.12500 Cao độ đáy: phương pháp nội suy tam giác (Triangulation interpolation) áp dụng nghiên cứu cao trình đáy sông nội suy xác định gần với giá trị thực đo Điều kiện biên tính toán thủy lực: biên chuỗi giá trị lưu lượng (Q) thực đo biên chuỗi giá trị mực nước (Z) thực đo (trong khoảng thời gian từ ngày 30/8/2012 đến ngày 31/8/2012) Bảng 3: Thành phần bùn cát đáy sông Kích Thước (mm) Thành phần phần trăm (%) 0.00100 0.00425 0.12500 100 Bảng 4: Thành phần mẫu bùn cát Điều kiện biên tính toán bồi lắng xói lở biên đầu vào: Kích thước (mm) 0.00100 0.00425 0.12500 Nồng độ bùn cát gáng 0.25 kg/m3 thành phần hạt bùn cát lơ lửng xác định Bảng (Wolanski et al., 1996); Bùn cát vận chuyển đáy 0.034 kg/m/s Hình 6: Lưới, biên, cao độ đáy sông, mặt cắt 12 Thành phần Độ rỗng phần trăm (%) 24.03 0.8 48.26 27.71 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 0.05 0.047 ∗ 2.5 Hiệu chỉnh kiểm định mô hình Tính toán thủy lực: mô hình hiệu chỉnh cách thay đổi hệ số nhám thủy lực Manning’s n dựa vào độ sâu cột nước Hệ số nhám thủy lực sông/kênh tự nhiên phù sa vùng đồng xác định nằm khoảng 0.020.05 (Chow, 1959) Dựa vào khoảng giá trị này, hệ số nhám thủy lực thay đổi thông qua phương pháp thử sai (Phạm Thế Bảo, 2009) tìm hệ số nhám thích hợp cho trình mô (giá trị mô gần giá trị thực đo) Hệ số nhám thủy lực tính phương pháp giảm dần theo chiều sâu cột nước (Picado et al., 2010); đó, (i) nút lưới có độ sâu cột nước nhỏ (hmin) có hệ số nhám thủy lực lớn (nmax); (ii) nút lưới có độ sâu cột nước lớn (hmax) có hệ số nhám thủy lực nhỏ (nmin); và, (iii) vị trí khác hệ số nhám thủy lực thay đổi tuyến tính bước (Wendt, 2008) theo cột nước giới hạn [nmin tới nmax] (PT 7) n n n n Trong đó: Ls: chiều dài ô lưới tính phù hợp (m); d50: kích thước trung bình hạt bùn cát đáy sông (mm); : tỷ trọng bùn cát (kg/m3); : tỷ trọng nước (kg/m3); : trọng lượng riêng bùn cát (N/m3); và, : trọng lượng riêng nước (N/m3) Wu Li (1992) đề nghị hệ số phù hợp bùn cát lơ lửng: =1 cho trường hợp xói lở mạnh; =0.25 cho trường hợp lắng động mạnh; và, =0.5 cho trường hợp bồi lắng xói lở yếu Do khu vực hạ lưu sông Hậu thường xảy trình bồi lắng (Nguyễn Văn Lập et al., 2000; Zuo Xue et al., 2012) nên hệ số =0.25 lựa chọn để tính toán trình bồi lắng xói lở Theo nghiên cứu Xue et al (2010) bùn cát lơ lửng cửa sông Định An có độ rỗng cao thành phần bùn cát chủ yếu sét, bùn cát Sự bồi lắng thành phần bùn cát có độ rỗng khoảng 0.8 (Rijn, 1993); đó, mô hình này, giá trị độ rỗng (bằng) 0.8 áp dụng Trong điều kiện giới hạn số liệu để hiệu chỉnh kiểm định mô hình, bên cạnh việc so sánh số liệu mực nước mô mực nước thực đo để kiểm chứng tính phù hợp mô hình khía cạnh thủy lực, kết mô bồi / xói bờ sông mô hình dựa vào điều kiện khảo sát thực tế đồng (PT7) Trong đó: ni: hệ số nhám thủy lực nút i (s/m1/3); nmin, nmax: hệ số nhám thủy lực lớn nhỏ (s/m1/3); (hmin) (hmax): độ sâu cột nước lớn nhỏ (m); và, hi: độ sâu cột nước nút thứ i (m) Sau hiệu chỉnh, mô hình kiểm định giá trị mực nước đo đạc từ 30 phút ngày 30/8/2012 đến 30 phút ngày 31/8/2012 Phương pháp kiểm định dựa tảng việc ứng dụng mô hình kiểm định nội (internal validation) (Shading et al., 2002) KẾT QUẢ VÀ THẢO LUẬN 3.1 Kết tính toán thủy lực Mô hình hiệu chỉnh thông qua việc thay đổi hệ số nhám thủy lực Manning’s n; với hệ số nhám thủy lực nằm khoảng 0.018-0.030 (s/m1/3), kết mực nước mô điểm quan sát phù hợp với kết thực đo (Hình A) (hệ số Nash Sutcliffe NS = 0.97) Ngoài ra, hệ số nhám áp dụng mô hình phù hợp với hệ số nhám nghiên cứu trước Trần Quốc Đạt, et al (2012) Lâm Mỹ Phụng et al (2013) Bên cạnh đó, mô hình kiểm định (Hình B) (hệ số Nash Sutcliffe NS=0.98); điều cho thấy mô hình sử dụng để mô động thái thủy lực dòng chảy nhánh sông nghiên cứu Tính toán xói lở: mô hình bồi xói cửa sông Định An, hệ số cần quan tâm bao gồm: hệ số phù hợp bùn cát lơ lửng , độ rỗng mẫu bùn cát p’ chiều dài lưới tính phù hợp bùn cát đáy sông Ls; đó, chiều dài lưới tính tính sau (PT8) (Qamar & Baig, 2012): 32.644 (PT8) ∗ Trong đó: ∗ 0.671 3,110.804 13 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 Hình 7: Mực nước thực đo mô điểm quan sát Vận tốc dòng chảy mặt cắt (Hình 6) đạt giá trị lớn thời điểm 18 (Hình A) đạt giá trị nhỏ thời điểm 22 (Hình B); kết phù hợp nghiên cứu Lane et al (1999) - vận tốc dòng chảy mặt cắt ngang tăng dần theo độ sâu Ngoài ra, đoạn sông nghiên cứu tương đối thẳng nên mực nước phân đoạn mặt cắt chênh lệch đáng kể (Hình 9) Hình 8: Vận tốc dòng chảy mặt cắt thời điểm 18 (A) 20 (B) Hình 9: Mực nước mặt cắt thời điểm 18 (A) 20 (B) 3.2 8.5 x 10-6m (Hình 10 B), bồi lắng bờ Trà Vinh 2.1 x 10-5m (Hình 10 B) Kết tính toán bồi lắng xói lở Kết tính toán bồi lắng xói lở mặt cắt (Hình 6) cho thấy sau thời gian mô có xói lở bờ Cù Lao Dung cụ thể 1.9 x 10-6 m (Hình 10 A), xói lở bờ Cù Lao Dung mặt cắt Sự xói lở bờ trái hai vị trí (mặt cắt 3) phù hợp với vị trí xác định xói lở bờ chuyến khảo sát thực địa (Hình B) Xu hướng bồi lắng 14 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 gần cửa sông mặt cắt số phù hợp với kết nghiên cứu trước Xue et al (2012) bùn cát lấy dọc bờ sông vận chuyển vùng cửa sông bồi lắng Hình 10: Thay đổi đáy sông sau thời gian mô nguyên gây xói lở bờ Cù Lao Dung thực tế Trong phần lớn thời gian mô phỏng, hướng dòng chảy mặt cắt (Hình 11 A) mặt cắt (Hình 11 B) hướng vào bờ Cù Lao Dung - Hình 11: Hướng vận tốc dòng chảy mặt cắt 1(A) (B) bờ Cù Lao Dung vận tốc dòng chảy mặt cắt số (Hình 12 B) phía bờ Trà Vinh có hướng dòng chảy vào bờ Trà Vinh – nguyên nhân gây bồi lắng vị trí Vận tốc dòng chảy mặt cắt số (Hình 12 A) hướng vào bờ Trà Vinh (TV) nhiều so với hướng sông Tuy nhiên, vận tốc hướng vào bờ nhỏ nên làm xói lở bờ Trà Vinh; Hình 12: Hướng vận tốc dòng chảy mặt cắt 1(A) (B) bờ bờ Trà Vinh 15 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 Change Studies Journal of Hydraulic Engineering, 125(9): 924-933 Lâm Mỹ Phụng, Văn Phạm Đăng Trí, & Trần Quốc Đạt (2013) Ứng dụng mô hình toán thủy lực chiều đánh giá dự báo tình hình xâm nhập mặn hệ thống sông địa bàn tỉnh Trà Vinh Tạp chí Khoa học- Trường Đại học Cần Thơ, 25 (2013): 68-75 Lane, S N., Bradbrook, K F., Richards, K S., Biron, P A., & Roy, A G (1999) The Application of Computational Fluid Dynamics to Natural River Channels: Three-Dimensional Versus TwoDimensional Approaches Geomorphology, 29: 1-20 Langendoen, E J.(2001).Evaluation of the effectiveness of selected computer models of depth-averaged free surface flow and sediment transport to predict the effects of hydraulic structures on river morphology.Project Report, USDA-ARS National Sedimentation Laboratory, Oxford M.S 10 Nguyễn Cảnh Cầm, Lưu Công Đào, Nguyễn Văn Cung, Nguyễn Như Khuê, Võ Xuân Minh, Hoàng Văn Quý, & Vũ Văn Tảo (2007) Thủy Lực Tập Hà Nội Nhà Xuất Bản Xây Dựng.323 11 Nguyễn Văn Lập, Tạ Thị Kim Oanh, & Tateishi, M (2000) Late Holocene Depositional Environments and Coastal Evolution of the Mekong River Delta, Southern Vietnam Journal of Asian Asian Earth Scienes, 18(2000): 427-439 12 Nguyễn Việt Thanh, Hai, Z J., & Hau, L P (2011) Morphological evolution of navigation channel in Dinh An estuary, Vietnam River, Coast and Estuarne Morphodynamics: RCEM2011 Tsinghua University Press, Beijing.469-482 13 Phạm Thế Bảo.(2009).Các phương pháp giải toán máy tính Khoa Toán – Tin,Trường Đại học Khoa học Tự nhiên 14 Picado, A., Dias, J M., & Fortunato, A B (2010) Tidal Changes in Estuarine Systems Induced by Local Geomorphologic Modifications Continental Sheft Reseach, 30(2010): 1854-1864 15 Qamar, M U., & Baig, F (2012) Calibration of CCHE2D for sediment simulation of Tarbela Reservoir KẾT LUẬN Bồi lắng xói lở tình trạng chung cửa sông ven biển ĐBSCL với quy mô không gian thời gian khác Kết tính toán cho thấy vùng cửa sông Định An có xu hướng xói lở dọc bờ Cù Lao Dung Kết nghiên cứu cho thấy mô hình thủy lực hai chiều (CCHE2D) mô động thái thủy lực dòng chảy đoạn sông nghiên cứu với độ tin cậy cao, có triển vọng ứng dụng nghiên cứu hệ thống nhánh sông ĐBSCL Tuy nhiên, hạn chế thời gian số liệu đo đạc, nghiên cứu thực nhằm mô đặc tính thủy lực dòng chảy (trên đoạn sông nghiên cứu) thời gian ngắn nên chưa phản ánh đầy đủ đặc tính thủy lực thay đổi hình thái lòng sông; vậy, kết nghiên cứu chứng minh biến động địa mạo vùng cửa sông (phù hợp với kết khảo sát thực địa) Hiện nay, phần lớn ứng dụng mô hình thủy lực ĐBSCL bị hạn chế mặt số liệu (số liệu đo đạc cũ số liệu thực đo) Trong bối cảnh có thay đổi lớn lưu lượng hình thái lòng sông, việc áp dụng thiết bị đo đạc đại (như thiết bị ADCP, thiết bị đo mực nước sensor) để thu thập số liệu thực tế nhằm mục tiêu phục vụ cho việc nghiên cứu động thái địa mạo thủy lực điều cần quan tâm TÀI LIỆU THAM KHẢO Albers, T., & Lieberman, N v (2011) Curent and Erosion Modelling Servey Deutsche Gesellschaft für Internationale Zusammenarbeit (GiZ) GmbH.73 Bộ Tài Nguyên Môi Trường (2012) Kịch Biến đổi Khí hậu, Nước biển dâng cho Việt Nam Nhà xuất Tài nguyênMôi trường Bản đồ Việt Nam.84 Bùi Hồng Long, & Tống Phước Hoàng Sơn.(2003).Đặc điểm địa hình biến động luồng lạch vùng cửa sông Định An.107-127 Chow, V T (1959) Open Channel Hydraulics: Blackburn Press.700 Hoa Mạnh Hùng, Nguyễn Quang Thành, & Phan Thị Thanh Hằng (2008) Động lực phát triển vùng cửa sông Hậu (cửa Định An- Tranh Đề) Các khoa học trái đất, 30: 130-135 Jia, Y., & Wang, S S Y (1999) Numerical Model for Channel Flow and Morphological 16 Tạp chí Khoa học Trường Đại học Cần Thơ Phần A: Khoa học Tự nhiên, Công nghệ Môi trường: 31 (2014): 8-17 22 Wendt, J F (2008) Computational Fluid Dynamics 72 Chaussée de Waterloo B1640 Rhode- Saint- Genèse: Belgium.299 23 Wilson, M D., & Atkinson, P M (2007) The Use of Remotely Sensed Land Cover to Derive Floodplain Friction Coefficients for Flood Inundation Modelling Hydrological Processes, 21(26): 3576-3586 24 Wolanski, E., Nguyen Ngoc Huan, Le Trong Dao, Nguyen Huu Nhan, & Nguyen Ngoc Thuy (1996) Fine- sediment Dynamic in the Mekong River Estuary, Vietnam Estuarine, Coastal and Shelt Science., 43: 565-582 25 Wu, W (2001) CCHE2D Sediment Transport Model School of Engineering The University of Mississippi.43 26 Wu, W (2007) Computational River Dynamics Taylor & Francis, London 2007 27 Wu, W., & Li, Y (1992) One and twodimensional nesting model for river flow and sedimentation Proc 5th Int.Symp on River Sedimnet 28 Xue, Z., He, R., Liu, J P., & Warner, J C (2012) Modeling transport and deposition of the Mekong River sediment Continental Shelf Research, 37(2012): 66-78 29 Xue, Z., Liu, J P., DeMaster, D., Lap, N V., & Oanh, T T K (2010) Late Holocene Evolution of the Mekong Subaqueous Delta, Southern Vietnam Marine Geology, 18: 427-439 Proceedings of the World Congress on Engineering I: 978-988 16 Rijn, L C V (1993) Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas The Netherlands: University of Utrecht Department of Physical Geography.700 17 Shading, W R., Cook, T D., & Campbell, D T (2002) Experimental and QuasiExperimental Designs for Generilized Causal Inference Boston Houghton Mifflin, Boston, New Yortk 18 Teledyne RD Instrument (2007) Winriver II User's Guide 19 Trần Quốc Đạt, Nguyễn Hiếu Trung, & Likitdecharote, K (2012) Mô Phỏng Xâm Nhập Mặn Đồng Bằng Sông Cửu Long Dưới Tác Động Mực Nước Biển Dâng Và Sự Suy Giảm Lưu Lượng Từ Thượng Nguồn Tạp Chí Khoa Học, Trường Đại Học Cần Thơ, 21b: 141-150 20 Trịnh Công Văn (2009) Identification Of Sea Level Rise Impacts on the Mekong Delta And Orientation Of Adaptation Activities VietNam National Commitee On Large Dams And Water Resousces Development.1-6 21 Văn Phạm Đăng Trí, Popescu, I., Griensven, A V., Solomatine, D., Nguyễn Hiếu Trung, & Green, A (2013) A Study of the Climate Change Impacts on Fluvial Flood Propagation in the Vietnamese Mekong Delta Hydrol Earth Syst, 9: 7227-7270 17