3.3.3.1.Hiệu chỉnh mơ hình
Số liệu mực nƣớc sóng, dịng chảy thực đo tại trạm C và mực nƣớc tại trạm D thời gian đo đạc từ ngày 15/11 - 28/11/2015 sẽ đƣợc dùng để hiệu chỉnh mơ hình.
Nhằm đơn giản hóa, luận văn sử dụng các số liệu trích xuất từ ECMWF tại điểm gần nhất với khu vực tính tốn: Điểm 1 (13.25o, 109.375o), điểm 2 (15.125o,
Hình 3.27a. Lưới Mesh II cho miền tính lớn lớn
Hình 3.27b. Lưới tính Mesh I cho khu vực nghiên cứu nghiên cứu
109.5o), điểm 3 (13o, 109.5o) làm điều kiện biên cho các tính tốn phục vụ hiệu chỉnh mơ hình. Tuy nhiên, do điểm trích xuất cịn ở khá xa so với khu vực tính tốn, đồng thời các giá trị từ ECMWF có xu hƣớng thiên bé so với các kết quả phân tích và thực đo, do vậy nếu sử dụng trực tiếp các số liệu này sẽ dẫn đến những sai số hệ thống trong q trình mơ phỏng đặc biệt ở các hiện tƣợng thời tiết dị thƣờng (bão, áp thấp nhiệt đới, gió mùa Đơng Bắc...). Vì thế, về nguyên tắc cần tiến hành mơ hình hóa cho khu vực rộng hơn hoặc tồn bộ Biển Đơng và sử dụng kỹ thuật lƣới lồng cho khu vực nghiên cứu.
Trong khuôn khổ luận văn chỉ mô phỏng trƣờng thủy động lực khu vực cửa Đà Diễn với lƣới Mesh I và các điều kiện biên ngoài biển sẽ sử dụng các kết quả tính tốn của lƣới Mesh II đƣợc thực hiện bởi nhóm nghiên cứu thuộc đề tài ĐTĐL.CN 15/15.
Qua q trình mơ phỏng chế độ thủy thạch động lực học cửa Đà Diễn, dựa vào các điều kiện biên và bộ thơng số mơ hình đƣợc hiệu chỉnh. Kết quả hiệu chỉnh mơ hình đƣợc trích tại các vị trí trạm C và trạm D so sánh với số liệu thực đo cho thấy, mực nƣớc đã đƣợc mô phỏng khá tốt, chỉ số Nash-Sutcliffe ở trạm C là 0,89, đƣờng quá trình mực nƣớc tính tốn và thực đo phù hợp về pha dao động, đỉnh và chân triều khá sát nhau (Hình 3.28, Hình 3.31).
Các mơ phỏng về dịng chảy tại trạm C cũng cho thấy có sự phù hợp cả về xu hƣớng và pha dao động tuy còn một số sai khác nhất định có thể sự phân chia tỷ lệ lƣu lƣợng tại biên dịng chảy trong sơng vào các lạch chƣa sát với điều kiện thực tiễn.
Hình 3.28. Mực nước trạm C và mực nước tính tốn được từ mơ hình
Hình 3.29. So sánh dịng chảy tính tốn với thực đo ở trạm C
Hình 3.31. So sánh mực nước tính tốn với thực đo ở trạm D
Các mơ phỏng về sóng nhìn chung có sự tƣơng đồng khá cao, tuy nhiên có một số đỉnh sóng thực đo mà trong tính tốn chƣa thể hiện đƣợc có thể là do việc bố trí trạm C ở gần cửa biển nơi thƣờng xuyên có tàu bè qua lại nên xuất hiện một số các đỉnh nghi là do sóng tàu gây nên. Xu thế tƣơng tự cũng có thể quan sát thấy ở các trạm đo cịn lại, điều đó chứng tỏ bộ thơng số đã lựa chọn và các dữ liệu địa hình đã bám sát thực đo và có thể sử dụng cho bƣớc đánh giá kiểm định tiếp theo.
3.3.3.2.Kiểm định mơ hình
Sử dụng lƣới, địa hình của đợt tháng 3 năm 2016 và bộ thông số mô đƣợc giữ ngun sau q trình hiệu chỉnh để mơ phỏng lại trƣờng thủy động lực khu vực nghiên cứu giai đoạn từ 18/5 -31/5/2016 nhằm mục đích kiểm định lại mơ hình. Cụ thể, mơ hình đƣợc kiểm định với độ cao sóng, dịng chảy thực đo tại trạm E và mực nƣớc, độ cao sóng tại trạm F (Hình 3.32 đến Hình 3.35).
Hình 3.32. Dịng chảy tính tốn với thực đo ở trạm E
Hình 3.33. Độ cao sóng tính tốn với thực đo ở trạm E
Hình 3.35. Độ cao sóng tính tốn với thực đo ở trạm F
Kết quả cho thấy, tƣơng tự nhƣ khi hiệu chỉnh, mực nƣớc tính tốn và thực đo tại trạm F tƣơng đối bám sát nhau, hệ số Nash-Sutcliffe tính tốn đƣợc là 0,84, đạt kết quả tốt. Tuy nhiên, ở các yếu tố về dịng chảy và sóng, vẫn cịn sự khác biệt khá rõ nét giữa tính tốn và thực đo tại các trạm E, F. Tuy nhiên, điều đó có thể giải thích đƣợc là do hạn chế về thời gian, trong khuôn khổ luận văn. Các kết quả nêu trên có thể chấp nhận đƣợc để sử dụng cho các mô phỏng trƣờng thủy động lực vùng cửa Đà Diễn theo các kịch bản tính tốn.
3.3.4. Mơ phỏng chế độ động lực và diễn biến bồi xói khu vực từ cầu Đà Rằng mới ra đến cửa Đà Diễn mới ra đến cửa Đà Diễn