7. Cấu trúc của luận án
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 toá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). Việc tính toá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].
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)
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 ngoà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 toá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 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
(b) (a)
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 ngoà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 ngoà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 ngoà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 toá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 toá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 ngoà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.
Mô hình tính được lựa chọn kiểu liên kết đồng thời của các quá trình cơ bản TĐL 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 toá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 quá trình sóng đổ [94].
Các tham số tính toá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 toá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 toán với cách tiếp cận HLES (mô phỏng xoá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-є.
Tham số tính toán của mô hình TTLL
- 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 (). Trong quá trình tính toán, vận tốc lắng đọng sẽ tính đến cả
những ảnh hưởng do độ mặn.
- Tiêu chuẩn ứng suất cho quá trình xói của trầm tích () 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 quá trình bồi lắng của trầm tích () 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 toá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 toá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 toá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) Giá trị E tiến tới 1 thì các kết quả dự báo có hiệu quả tốt nhất, ngược lại khi E tiến 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].
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 khuôn 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 toá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 toán dự báo khác nhau đã được xây dựng. Các kịch bản tính toán dự báo gồm 2 nhóm:
+ Nhóm kịch bản tính toá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 toá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.