CHƯƠNG 2 PHƯƠNG PHÁP NGHIÊN CỨU VÀ TÀI LIỆU
2.2. Phương pháp nghiên cứu
2.2.2. Phương pháp mơ hình
Các yếu tố thủy động lực (TĐL), vận chuyển TTLL được mơ hình hóa trên cơ sở mơ hình Delft3D do Viện Thủy lực Delft (Hà Lan) nghiên cứu phát triển với các module cơ bản như TĐL (Delft3d-Flow), sóng (Delft3d-Wave), vận chuyển bùn cát (Delft3d- Sed) và các module khác. Mơ hình này có thể mơ phỏng tốt điều kiện TĐL - sóng, vận chuyển bùn cát, chất lượng nước ở vùng cửa sơng ven bờ [94].
Hình 2. 1. Tương tác sóng- dịng chảy và vận chuyển trầm tích trong mơ hình Delft3D
Module TĐL (Delft3d-Flow) có thể tính tốn kết hợp đồng thời (online coupling) với các module khác như sóng (Delft3d-Wave), vận chuyển trầm tích (Delft3d-Sed).
Vận chuyển trầm tích
(Delft3d-Sed) Địa hình đáy
(Bottom) Thủy động lực (Delft3d-Flow) Sóng (Delft3d-Wave) Tương tác sóng- dịng chảy
Việc tính tốn kết hợp đồng thời có thể cho thấy được sự tương tác giữa các quá trình TĐL - sóng và vận chuyển trầm tích (hình 2.1) tại mỗi thời điểm tính của mơ hình [94].
Hình 2. 2. Các lưới tính của mơ hình thủy động lực
(a-lưới thơ phía ngồi; b- lưới chi tiết cho vùng ven biển Hải Phòng)
Trong nghiên cứu này, kiểu tính tốn kết hợp đồng thời (online –coupling) của các module chính là TĐL, sóng và vận chuyển trầm tích đã được sử dụng để đánh giá đặc điểm vận chuyển TTLL.
Lưới tính cho vùng bờ biển Hải Phịng
Mơ hình thuỷ động lực cho khu vực cửa sơng ven bờ Hải Phịng sử dụng hệ lưới (a)
cong trực giao. Phạm vi vùng tính của mơ hình bao gồm các vùng nước của các cửa sơng ven biển trải dài từ vùng phía bắc khu vực vịnh Hạ Long đến phía nam cửa Trà Lý. Miền tính có kích thước khoảng 106 km theo chiều Đơng Bắc-Tây Nam và 64 km theo chiều Tây Bắc-Đơng Nam, với diện tích mặt nước khoảng 5085km2 được chia thành 628 x 488 điểm tính và các ơ lưới có kích thước biến đổi từ 8,3 đến 340m (hình 2.2-b). Các ơ lưới tính theo chiều thẳng đứng sử dụng hệ toạ độ . Trong mơ hình thuỷ động lực cho vùng cửa sơng ven biển Hải Phịng, độ sâu của các điểm tính được chia thành 4 lớp nước với tỷ lệ đều nhau từ mặt xuống đáy là 25% độ sâu cho mỗi lớp nước.
Lưới độ sâu cho mơ hình tính ở khu vực này là file số liệu địa hình đã được xử lý, gắn với lưới tính của mơ hình. Lưới tính của mơ hình thơ phía ngồi dùng để NESTHD các điều kiện biên mở phía biển cho mơ hình lưới chi tiết khu vực nghiên cứu cũng là hệ lưới cong trực giao. Phạm vi vùng tính của mơ hình này mở rộng ra phía ngồi gần với lưới tính thơ hơn (hình 2.2-a). Miền tính này có kích thước khoảng 129km theo phương Đoong-Tây và 122km theo phương Bắc-Nam, diện tích mặt nước khoảng 15738km2 được chia thành 608 x 605 điểm tính với các ơ lưới có kích thước biến đổi từ 9,3 đến 1800,39m. Lưới tính theo chiều thẳng của mơ hình này cũng được chia thành 4 lớp nước với tỷ lệ đều nhau từ mặt xuống đáy là 25%.
Lưới độ sâu của mơ hình thuỷ động lực khu vực nghiên cứu (detailed model) và mơ hình phía ngồi dùng để NESTHD (overall model) được tạo ra bởi trường độ sâu, kết hợp với lưới tính và chương trình Delft3D- QUICKIN trong bộ phần mềm Delft3d.
Điều kiện biên và điều kiện ban đầu
Trong mơ hình Delft3D, điều kiện ban đầu của mơ hình có thể sử dụng từ kết quả tính tốn của các lần chạy trước đó thơng qua các restart file. Đối với trường hợp áp dụng cho vùng ven bờ Hải Phòng, điều kiện ban đầu của các kịch bản hiện trạng là các kết quả tính tốn trong khoảng thời gian trước đó 1 tháng. Đối với các biên sông, các giá trị lưu lượng nước, độ mặn, nhiệt độ, hàm lượng TTLL trung bình theo mùa được sử dụng. Đối với các biên mở phía biển, các kết quả tính dao động mực nước, độ muối, nhiệt độ của mơ hình phía ngồi (phương pháp NESTHD) được áp dụng. Số liệu nhiệt muối của mơ hình lấy từ cơ sở dữ liệu WOA13 với độ phân giải 0,25 độ cho khu vực Biển Đơng.
lực- sóng- vận chuyển bùn cát [94, 95]. Trong đó các yếu tố chính được tính đến bao gồm: độ muối; nhiệt độ; ảnh hưởng của gió bề mặt; tương tác với sóng (tính đồng thời kết hợp sóng-online coupling) và TTLL (tính đồng thời). Trong luận án, chỉ xét tới ảnh hưởng của gió, các yếu tố khí tượng khác như độ ẩm, lượng mưa, bức xạ, nhiệt độ khơng khí khơng tính đến trong mơ hình. Số liệu gió đưa vào mơ hình tính là các số liệu quan trắc tại Bạch Long Vỹ và Hòn Dáu trong năm 2016 và 2017 với tần suất 6h/lần.
Mơ hình sóng
Mơ hình sóng được thiết lập chạy đồng thời (online coupling) với mơ hình TĐL và mơ hình vận chuyển trầm tích. Tại mỗi thời điểm tính tốn (1h), mơ hình sóng sử dụng lưới tính, trường gió, các kết quả tính độ sâu, mực nước, dịng chảy của mơ hình TĐL.
- Điều kiện biên mở của mơ hình sóng cịn được tham khảo từ kết quả dự báo sóng Wave Climate ở vùng vịnh Bắc Bộ [96] trong thời gian 2014-2019.
- Kiểu ma sát đáy trong mơ hình sóng được lựa chọn là phổ JONSWAP với hệ số ma sát đáy có giá trị 0,067 [94]. Mơ hình B&J [97] được lựa chọn để tính ảnh hưởng của nước nơng nơi diễn ra q trình sóng đổ [94].
Các tham số tính tốn khác của mơ hình
- Tham số nhám đáy (bottom roughness) được lựa chọn sử dụng các hệ số Manning (n) biến đổi theo không gian với giá trị 0,018-0,023 m-1/3s. Các hệ số manning lớn hơn ở điều kiện trầm tích đáy là vật liệu thơ và nhỏ hơn ở điều kiện trầm tích đáy là hạt mịn. Căn cứ để tính tốn các hệ số này dựa vào phân bố trầm tích tầng mặt và Tài liệu hướng dẫn lựa chọn các hệ số Manning [98-99].
- Các giá trị liên quan đến điều kiện rối có thể được xác định như là một hằng số, hoặc tham số biến đổi theo khơng gian hoặc tính tốn với cách tiếp cận HLES (mơ phỏng xốy lớn bình lưu - Horizontal Large Eddy Simulation). Mơ hình HLES gần đây đã được tích hợp trong hệ thống mơ hình Delft3D theo lý thuyết của Uittenbogaard [100] và được thảo luận trong nghiên cứu của Van Vossen [101]. Trong nghiên cứu này, hệ số khuyếch tán rối và nhớt rối nền theo phương ngang được lựa chọn là 10m2/s. Các hệ số này theo phương thẳng đứng là 10-5m2/s. Mơ hình khép kín rối 2 chiều là mơ hình HLES trong Delft3D. Mơ hình khép kín rối 3 chiều là mơ hình k-є.
- Vận tốc lắng đọng của TTLL được chọn là 0,1mm/s. Đây là giá trị vận tốc lắng đọng trong nước ngọt (ws,f). Trong q trình tính tốn, vận tốc lắng đọng wssẽ tính đến
cả những ảnh hưởng do độ mặn.
- Tiêu chuẩn ứng suất cho q trình xói của trầm tích (c,e) biến đổi trong khoảng
từ 0,1-1,0 N/m2 [102]. Sau các lần hiệu chỉnh, tiêu chuẩn xói được lựa chọn là 0,25 N/m2. - Tiêu chuẩn ứng suất cho q trình bồi lắng của trầm tích (c,d) biến đổi trong khoảng từ 0,005-0,25 N/m2 [102]. Sau các lần hiệu chỉnh, tiêu chuẩn xói được lựa chọn là 0,1N/m2.
- Tốc độ xói trong tự nhiên đo đạc được biến đổi trong khoảng 10-5-10-3kg/m2s. Với tỷ trọng bùn cát đáy là 2650kg/m3, tỷ trọng TTLL gần lớp biên đáy là 500kg/m3, tốc độ xói ban đầu được giả thiết là 10-3 kg/m2.s.
Để đánh giá mức độ tin cậy trong các tính tốn, các chỉ số được sử dụng là hệ số tương quan Bravais-Pearson và chỉ số Nash và Sutcliffe.
Hệ số tương quan Bravais-Pearson được xác định theo cơng thức sau:
(1)
Trong đó: Oi – giá trị quan trắc, đo đạc;𝑂̅ - giá trị quan trắc trung bình Pi - giá trị tính tốn, dự báo; 𝑃̅ - giá trị trung bình dự báo.
r nằm trong khoảng từ 0-1. Các kết quả tính tốn dự báo tốt nhất khi r tiến tới 1 và độ tin cậy thấp nhất khi r tiến tới 0 [103].
Một chỉ số khác để đánh giá độ tin cậy của mơ hình là chỉ số hiệu quả của dự báo- chỉ số Nash và Sutcliffe [104]. Chỉ số này đánh giá số lượng dự báo đảm bảo độ tin cậy cho phép.
(2)
tới 0 thì các dự báo khơng đáng tin cậy. Khi E mang dấu âm (-), các đặc trưng trung bình tính từ chuỗi quan trắc cho kết quả dự báo tốt hơn từ mơ hình [103].
Số liệu mực nước dùng để hiệu chỉnh, kiểm chứng kết quả tính của mơ hình trong là các giá trị đo đạc tại trạm Hải văn Hòn Dáu trong năm 2015, 2016 với tần suất 1 giờ/lần. Ngoài số liệu mực nước, các chuỗi số liệu đo dòng chảy, hàm lượng TTLL được sử dụng của đề tài cấp thành phố “Nghiên cứu xây dựng luận cứ phục vụ lập qui hoạch
các bãi đổ bùn cát do nạo vét trên địa bàn Hải Phòng, MT.2015.721” và các kết quả
khảo sát về dịng chảy, hàm lượng TTLL trong khn khổ thực hiện đề tài này để kiểm chứng, hiệu chỉnh mơ hình.
Các kịch bản tính tốn dự báo
Dựa trên các căn cứ về điều kiện tự nhiên, quy hoạch khai thác cát, hiện trạng hoạt động khai thác cát ở vùng cửa sông ven biển Hải Phịng, các nhóm kịch bản tính tốn dự báo khác nhau đã được xây dựng. Các kịch bản tính tốn dự báo gồm 2 nhóm:
+ Nhóm kịch bản tính tốn mơ phỏng ảnh hưởng đến thay đổi điều kiện dịng chảy, sóng, vận chuyển bùn cát và bồi xói do thay đổi địa hình ở các vị trí khai thác cát.
+ Nhóm kịch bản mơ phỏng, dự báo phát tán chất gây ô nhiễm trong quá trình khai thác cát: khai thác với điều kiện hiện trạng (giả thiết 30% công suất) của các dự án đã cấp phép (khoảng 2.53 triệu m3/năm); dự báo khai thác hết 100% công suất của các dự án đã cấp phép (khoảng 8.429 triệu m3/năm); dự báo khai thác 100% công suất theo quy hoạch khai thác cát.
Thời gian mơ phỏng tính tốn cho mỗi kịch bản là 15 ngày bao gồm kỳ triều cường và triều kém. Đối với mỗi nhóm kịch bản tính cho 3 hướng sóng gió khác nhau: NE, SE và SW. Mỗi hướng gió sẽ được tính làm 2 nhóm trường hợp điều kiện gió bình thường và gió cực đoan.
Nhằm đánh giá trao đổi nước qua một số mặt cắt trong khu vực dưới ảnh hưởng của khai thác cát. Với mỗi kịch bản thời gian chạy là khoảng 15 ngày để tính đến tất cả các ảnh hưởng của kỳ triều cường và kỳ triều kém, bước thời gian tính (t): 60 giây.