Ngày nhận bài: 5/10/2025; Ngày sửa bài: 28/10/2025; Ngày chấp nhận đăng: 19/11/2025

Tóm tắt

Nghiên cứu này mô phỏng các kịch bản tràn dầu tại khu vực Bến Gót (Hải Phòng), tập trung vào các sự cố tiềm ẩn tại cảng Lạch Huyện trong hai điều kiện gió mùa đặc trưng: Gió mùa Đông Bắc (tháng 1/2022) và gió mùa Tây Nam (tháng 7/2022). Mô hình được Đăng nhập tài khoản 28Bet dựa trên phương pháp Lagrangian do nhóm tác giả phát triển, cho phép theo dõi quỹ đạo vệt dầu với độ phân giải cao.

Kết quả cho thấy, trong gió mùa Đông Bắc, dầu từ khu vực Lạch Huyện chủ yếu trôi về hướng Nam và Tây Nam, ảnh hưởng tới luồng Lạch Huyện, kênh Hà Nam, kênh Cái Tráp, luồng Nam Triệu và lan tới vùng Đồ Sơn - Hòn Dấu. Ngược lại, trong gió mùa Tây Nam, vệt dầu có xu hướng trôi về hướng Bắc và Tây Bắc, tác động đến luồng Lạch Huyện, kênh Hà Nam, đảo Cống (Cát Hải), đảo Quả Xoài và kéo dài đến phía Nam đảo Hoàng Tân (Quảng Yên, Quảng Ninh). Mô hình đã mô tả hiệu quả quá trình phát tán dầu trong điều kiện thủy động lực phức tạp, có xét đến tương tác giữa gió và thủy triều. Kết quả nghiên cứu cung cấp cơ sở khoa học quan trọng phục vụ công tác ứng phó tràn dầu khu vực, bảo vệ các vùng sinh thái nhạy cảm và hỗ trợ quản lý rủi ro Thể thao 28bet tổng hợp ở khu vực cửa sông và ven biển.

Từ khóa: Tràn dầu, bảo vệ Thể thao 28bet biển, mô phỏng.

Abstract

This study simulates oil spill scenarios in the Ben Got area (Hai Phong), focusing on potential incidents at Lach Huyen Port under two typical monsoon conditions: The northeast monsoon (January 2022) and the southwest monsoon (July 2022). The model, developed by the authors based on the Lagrangian approach, allows high-resolution tracking of oil slick trajectories.

Results show that under northeast monsoon conditions, oil from the Lach Huyen area mainly drifts southward and southwestward, affecting the Lach Huyen Channel, Ha Nam Canal, Cai Trap Canal, Nam Trieu Channel, and extending toward the Do Son - Hon Dau region. In contrast, during the southwest monsoon, the oil slick tends to drift northward and northwestward, impacting the Lach Huyen Channel, Ha Nam Canal, Cong Island (Cat Hai), Qua Xoai Island, and extending to the southern part of Hoang Tan Island (Quang Yen, Quang Ninh). The model effectively reproduces the oil dispersion process under complex hydrodynamic conditions, accounting for the combined effects of wind and tides. The findings provide important scientific support for regional oil spill response planning, the protection of sensitive ecological zones, and the integrated management of environmental risks in estuarine and coastal areas.

Keywords: Oil spill; Marine environment protection; Simulation.

1. Đặt vấn đề

Sự cố tràn dầu trên biển và trong các vùng cửa sông ven biển là một trong những rủi ro Thể thao 28bet nghiêm trọng, có thể gây hậu quả lâu dài đối với hệ sinh thái, nguồn lợi thủy sản, hoạt động kinh tế biển và đời sống cộng đồng ven bờ [1]. Đặc biệt tại Việt Nam, với đường bờ biển dài hơn 3.260 km cùng hệ thống cảng biển phát triển nhanh chóng, nguy cơ xảy ra các sự cố tràn dầu ngày càng gia tăng. Các nghiên cứu và mô phỏng sự cố tràn dầu đóng vai trò then chốt trong công tác dự báo, ứng phó kịp thời, đồng thời là cơ sở khoa học phục vụ quản lý rủi ro Thể thao 28bet và bảo vệ các khu vực nhạy cảm.

Khu vực Bến Gót - Lạch Huyện (Hải Phòng) là trung tâm hàng hải quan trọng của miền Bắc Việt Nam, nơi tập trung nhiều hoạt động vận tải, tiếp nhiên liệu và logistics, đồng thời nằm gần hệ thống cảng lớn như cảng Lạch Huyện, Đình Vũ và khu 28 Bet 22vip nổ hũ, xổ số, casino Deep C. Đây là vùng có điều kiện thủy động lực phức tạp do chịu ảnh hưởng đồng thời của chế độ thủy triều, gió mùa và dòng chảy cửa sông [2], khiến nguy cơ lan rộng của dầu tràn khi xảy ra sự cố càng cao. Một vụ tràn dầu tại khu vực này không chỉ ảnh hưởng đến hoạt động của các cảng và cơ sở 28 Bet 22vip nổ hũ, xổ số, casino ven bờ, mà còn đe dọa các khu sinh thái nhạy cảm lân cận như Đồ Sơn, Cát Hải, Cát Bà và vịnh Hạ Long [3].

Trên thế giới, nhiều nghiên cứu đã áp dụng mô hình toán học để mô phỏng quá trình phát tán dầu trong Thể thao 28bet biển và cửa sông. Trong đó, phương pháp Lagrange - dựa trên cơ sở theo dõi chuyển động của các hạt dầu (particle tracking) - được đánh giá cao nhờ khả năng mô phỏng chi tiết quỹ đạo và sự phân tán của vệt dầu dưới tác động của gió, dòng chảy và sóng [4]. Một số nghiên cứu trên thế giới [5-9] cũng sử dụng phương phương pháp này để mô phỏng tràn dầu trên biển và tại các khu vực cửa sông. Các mô hình này dần đã được phát triển 28BET Kênh Chính Thức hóa trên toàn cầu. Ngoài ra, phương pháp Euler cũng được sử dụng để mô phỏng vệt dầu loang trên biển. Phương pháp này quan sát sự biến thiên của dòng chảy tại các điểm cố định trong không gian (ô lưới) thông qua giải hệ phương trình Navier-Stokes tại các nút hoặc ô lưới cố định trong không gian. Trên thế giới có một số nghiên cứu sử dụng phương pháp này để ứng dụng vào mô phỏng tràn dầu [10-13].

Ở Việt Nam, các nghiên cứu về tràn dầu còn hạn chế và chủ yếu tập trung ở quy mô vùng biển rộng; trong khi các nghiên cứu chi tiết cho khu vực cửa sông - cảng biển như Bến Gót - Lạch Huyện vẫn còn ít [14]. Một số Tải app 28Bet trước đây sử dụng các phần mềm 28BET Kênh Chính Thức như MIKE21 hoặc DELFT3D để mô phỏng thủy động lực và sự phát tán dầu [15-16], song việc phát triển các công cụ mô hình riêng phù hợp với đặc thù chế độ gió mùa và thủy triều của khu vực miền Bắc vẫn còn hạn chế.

Do đó, nghiên cứu này tập trung mô phỏng các kịch bản tràn dầu tại khu vực Bến Gót (Hải Phòng) dưới hai điều kiện gió mùa điển hình: gió mùa Đông Bắc (tháng 01/2022) và gió mùa Tây Nam (tháng 07/2022). Thông qua việc phát triển và ứng dụng một mô hình Lagrange do nhóm tác giả Đăng nhập tài khoản 28Bet, nghiên cứu hướng đến ba mục tiêu chính: (i) phân tích đặc điểm lan truyền và phân tán dầu trong điều kiện thủy động lực đặc trưng; (ii) xác định các khu vực ven bờ và cảng biển có nguy cơ bị ảnh hưởng; và (iii) cung cấp cơ sở khoa học cho công tác ứng phó sự cố, bảo vệ Thể thao 28bet và quản lý rủi ro tại các vùng cửa sông ven biển.

Điểm mới của nghiên cứu là việc tích hợp đồng thời tác động của gió mùa, thủy triều và dòng chảy vào mô hình dự báo quỹ đạo dầu, cho phép đánh giá chính xác hơn phạm vi và mức độ ảnh hưởng của các sự cố tràn dầu. Kết quả không chỉ mang lại cơ sở khoa học phục vụ quy hoạch bảo vệ các vùng sinh thái nhạy cảm và nâng cao năng lực ứng phó tại khu vực Hải Phòng, mà còn đóng góp quan trọng cho sự phát triển của các nghiên cứu mô hình hóa tràn dầu tại Việt Nam.

2. Tài liệu và phương pháp

2.1. Tài liệu

Dữ liệu địa hình khu vực biển Hải Phòng và lân cận bao gồm:

- Các bản đồ địa hình UTM hệ tọa độ địa lý VN 2000 tỷ lệ 1:25.000, 1:100.000 và 1:500.000 do Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam (Bộ Nông nghiệp và Môi trường) kế thừa từ dự án nghiên cứu [17]. 

- Cơ sở dữ liệu GEBCO của Trung tâm Tư liệu Hải dương học Vương quốc Anh (British Oceanigraphic Data Centre) [18].

- Cơ sở dữ liệu địa hình toàn cầu ETOPO1 của Trung tâm Tư liệu Địa vật lí Quốc gia Hoa Kỳ NGDC (National Geophysical Data Center).

Sau đó, tiến hành chuyển các dữ liệu địa hình về cùng hệ thống kinh tuyến, vĩ tuyến theo hệ tọa độ UTM 48; dữ liệu độ cao lấy theo mốc cao độ quốc gia. Các dữ liệu này được sử dụng làm đầu vào cho việc tạo lưới tính và trường độ sâu trên mô hình MIKE 21/3 Coupled FM phục vụ cho tính toán ở phần tiếp theo.

2.2. Phương pháp mô phỏng thủy động lực

MIKE 21/3 Coupled FM là một hệ thống mô hình tính toán dòng chảy hai chiều (2D - MIKE 21) và ba chiều (3D - MIKE 3) sử dụng lưới tam giác phi cấu trúc (FM - Flexible Mesh) để mô phỏng các hiện tượng thủy động lực, lan truyền chất ô nhiễm, vận chuyển trầm tích và các hiện tượng biển - ven bờ phức tạp. Điều này đặc biệt phù hợp với điều kiện thực tế tại vùng nghiên cứu ven biển Việt Nam, nên mô hình được sử dụng rộng rãi trong nghiên cứu, quy hoạch, quản lý Thể thao 28bet biển, cửa sông, hồ chứa và vùng ven biển.

Mô hình bao gồm các mô-đun thành phần sau:

- Hydrodynamic (HD): Mô phỏng dòng chảy, mực nước, vận tốc, phân bố mặn - nhiệt độ [19].

- Transport (TR): Mô phỏng sự lan truyền của các chất hòa tan hoặc huyền phù.

- Sediment Transport (ST): Tính toán quá trình vận chuyển trầm tích (bùn/cát), xói lở, bồi tụ.

- Wave (SW): Mô phỏng sóng ngắn bằng phương trình năng lượng hoặc phương trình phân bố phổ.

- ECO Lab/Oil Spill: Mô phỏng quá trình sinh thái học, sự lan truyền và phân rã của dầu trong Thể thao 28bet nước [20].

Trong nghiên cứu này, tác giả sử dụng kết hợp mô-đun thủy động lực Hydrodynamic và mô-đun tính toán tràn dầu ECO Lab/Oil spill, cho phép mô hình hóa quá trình lan truyền dầu một cách hiệu quả, đồng thời có thể mô phỏng các tình huống sự cố khác nhau để hỗ trợ công tác ứng phó và dự báo.

2.3. Phương pháp mô phỏng tràn dầu

Nghiên cứu dựa trên phương pháp mô phỏng rời rạc các chất điểm Lagrange để tính toán và dự đoán vị trí của hạt dầu trong Thể thao 28bet nước biển [21]. Một trường dòng chảy hai chiều được xem xét:

 Trong đó, X_n+1 và Y_n+1 là vị trí hạt dầu theo trục x và y tại thời điểm t+dt; X_n và Y_n là vị trí hạt dầu theo trục x và y tại thời điểm t; u_part và v_part là vận tốc trôi do dòng chảy gây ra; b là hệ số phân tán; x  là giá trị ngẫu nhiên được xác định biến đổi từ -1 đến 1; u’ và v’ là thành phần quán tính theo phương x và y của dòng chảy và gió tại thời điểm n+1; dt là bước thời gian; t là thời gian (giờ).

Phương pháp Lagrange thể hiện rõ ưu điểm trong các ứng dụng mô phỏng vật chất phi liên tục như vết dầu, rác vi thải nhựa, hay vi sinh vật. Khả năng theo dõi riêng lẻ từng hạt giúp mô phỏng chính xác sự di chuyển và phân bố không đồng đều của vật chất trong Thể thao 28bet nước. Đồng thời, khả năng loại bỏ linh hoạt các hạt không còn ý nghĩa vật lý (quá nhỏ, đã phân rã...)…) giúp tối ưu tài nguyên và tăng tốc độ tính toán, đặc biệt hữu ích trong các mô hình quy mô lớn.

3. Kết quả

3.1. Kết quả mô phỏng thủy lực thời kỳ gió mùa Đông Bắc (tháng 01/2022)

Hình 1 là biểu đồ so sánh các giá trị mực nước mô phỏng (Hsim) và giá trị mực nước quan trắc theo trạm quốc gia (Hobs_01), kết quả cho thấy các giá trị này có độ tương đồng cao về pha và độ lớn. Hình 2 thể hiện kết quả mô phỏng mực nước tại thời kỳ triều cường và triều kiệt tháng 01/2022.

3.2. Kết quả mô phỏng thủy lực thời kỳ gió mùa Tây Nam (tháng 07/2022)

Tương tự như kết quả mô phỏng ở Mục 3.1, các hình và bảng dưới đây thể hiện kết quả mô phỏng thủy lực mùa gió Tây Nam vào tháng 07/2022.

3.3. Kết quả mô phỏng tràn dầu thời kỳ gió mùa Đông Bắc (tháng 01/2022)

Tác giả tiến hành thực hiện mô phỏng cho trường hợp giả định kỳ gió mùa Đông Bắc tháng 01/2022. Trong kịch bản sự cố tại bến Gót, cửa Lạch Huyện từ 03/01/2022 đến 10/01/2022. Vệt dầu từ khu vực bến cảng Lạch Huyện, trôi dạt chủ yếu theo hướng Nam và Tây Nam ảnh hưởng tới khu vực luồng Lạch Huyện, kênh Hà Nam, kênh Cái Tráp, luồng Nam Triệu và trôi xuống khu vực Đồ Sơn, Hòn Dấu, như trong Hình 4.

3.4. Kết quả mô phỏng tràn dầu thời kỳ gió mùa Tây Nam (tháng 07/2022)

Tác giả cũng tiến hành thực hiện mô phỏng cho trường hợp giả định kỳ gió mùa Tây Nam tháng 07/2022. Trong kịch bản sự cố tại cảng Đình Vũ, cửa Bạch Đằng từ 04/07/2022 đến 11/07/2022. Vệt dầu từ khu vực bến cảng Lạch Huyện, trôi dạt chủ yếu theo hướng Bắc và Tây Bắc ảnh hưởng tới khu vực luồng Lạch Huyện, kênh Hà Nam, đảo Cống (Cát Hải), đảo Quả Xoài sau đó kéo dài lên phía nam đảo Hoàng Tân (Quảng Yên, Quảng Ninh), như trong Hình 5.

4. Kết luận

Nghiên cứu đã phát triển và ứng dụng mô hình Lagrange để mô phỏng các kịch bản tràn dầu tại bến Gót (Hải Phòng) dưới hai điều kiện gió mùa đặc trưng. Kết quả cho thấy, trong gió mùa Đông Bắc (01/2022), vệt dầu từ khu vực bến cảng Lạch Huyện, trôi dạt chủ yếu theo hướng Nam và Tây Nam ảnh hưởng tới khu vực luồng Lạch Huyện, kênh Hà Nam, kênh Cái Tráp, luồng Nam Triệu và trôi xuống khu vực Đồ Sơn, Hòn Dấu; còn trong gió mùa Tây Nam (07/2022), vệt dầu từ khu vực bến cảng Lạch Huyện, trôi dạt chủ yếu theo hướng Bắc và Tây Bắc ảnh hưởng tới khu vực luồng Lạch Huyện, kênh Hà Nam, đảo Cống (Cát Hải), đảo Quả Xoài sau đó kéo dài lên phía Nam đảo Hoàng Tân (Quảng Ninh). Mô hình đã phản ánh hiệu quả sự lan truyền dầu trong điều kiện thủy động lực phức tạp, đồng thời chứng minh tính khả thi khi so sánh với các mô hình 28BET Kênh Chính Thức.

Kết quả nghiên cứu mang lại cơ sở khoa học quan trọng cho công tác dự báo, ứng phó sự cố và bảo vệ các vùng sinh thái nhạy cảm. Đây là đóng góp mới trong việc phát triển công cụ mô hình hóa nội sinh, phù hợp với đặc điểm gió mùa khu vực miền Bắc Việt Nam. Trong thời gian tới, nhóm tác giả sẽ tích hợp thêm các quá trình vật lý – hóa học, kết hợp với dữ liệu quan trắc thực địa nhằm nâng cao độ tin cậy, mở rộng phạm vi ứng dụng cho các khu vực ven biển khác, đồng thời phát triển công cụ cảnh báo sớm hỗ trợ ứng phó nhanh các sự cố tràn dầu.

TÀI LIỆU THAM KHẢO

[1]. Beyer, Jonny, et al. “Environmental effects of the Deepwater Horizon oil spill: A review”, Marine pollution bulletin 110.1: 28-51, 2016.

[2]. Đinh Văn Huy, Đỗ Đình Chiến, Trần Đức Thạnh, Bùi Văn Vượng, “Đặc trưng hình thái, động lực và biến dạng bờ Hải Hậu, Nam Định”, Tuyển tập Tài nguyên và Môi trường biển, tập X, NXB. Khoa học và Kỹ thuật, 2003, Tr.106-125.

[3]. Tri, Doan Quang et al., “Application of environmental sensitivity index (ESI) maps of shorelines to coastal oil spills: a case study of Cat Ba Island, Vietnam”, Environmental Earth Sciences, 74.4: 3433-3451, 2005.

[4]. French, D.P.; Schuttenberg, H.Z.; Isaji, T, “Probabilities of oil exceeding thresholds of concern: Examples from an evaluation for Florida Power and Light”, In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Calgary, AB, Canada, 2 - 4 June 1999; pp. 243–270.

[5]. Brovchenko, I.; Kuschan, A.; Maderich, V.; Zheleznyak, M., “The modelling system for simulation of the oil spills in the Black Sea”, In Proceedings of the 3rd EuroGOOS Conference: Building the European Capacity in Operational Oceanography, Athens, Greece, 3 - 6 December 2002.

[6]. Carracedo, P.; Torres-Lopez, S.; Barreiro, M.; Montero, P.; Balseiro, C.F.; Penabad, E.; Leitao, P.C.; Perez-Munuzuri, V., “Improvement of pollutant drift forecast system applied to the Prestige oil spills in Galicia Coast (NW of Spain): Development of an operational system”, Mar. Pollut. Bull. 2006, 53, 350–360.

[7]. Dulière, V.; Legrand, S.; Ovidio, F., “Development of an Integrated Software for Forecasting the Impacts of Accidental Oil Pollution (OSERIT)”, Royal Belgian Institute of Natural Sciences: Brussels, Belgium, 2010.

[8]. De Dominicis, M.; Pinardi, N.; Zodiatis, G.; Lardner, R., “MEDSLIK-II, a Lagrangian marine surface oil spill model for short-term forecasting Part 1: Theory”, Geosci. Model Dev. 2013, 6, 1851–1869.

[9]. Mackay, D.; Paterson, S.; Trudel, K., “A Mathematical Model of Oil Spill Behaviour”, Report to Research and Development Division, Environment Emergency Branch, Environmental Impact Control Directorate; Environment Canada: Ottawa, ON, Canada, 1980.

[10]. Tkalich, Pavlo, Kamrul Huda, and Karina Yew Hoong Gin, “A multiphase oil spill model”, Journal of Hydraulic Research 41.2 (2003): 115-125.

[11]. Agrawal, Madhusuden, and Duraivelan Dakshinamoorthy, “Computational analysis of oil spill in shallow water due to wave and tidal motion”, Offshore Technology Conference. OTC, 2011.

[12]. Raznahan, Mohammadmehdi, et al., “Multiphase CFD simulation of the nearshore spilled oil behaviors”, Environmental pollution 288 (2021): 117730.

[13]. Sarhadi Zadeh, Ehsan, and Kourosh Hejazi, “Eulerian Oil Spills Model Using Finite‐Volume Method with Moving Boundary and Wet‐Dry Fronts”, Modelling and Simulation in Engineering 2012.1 (2012): 398387.

[14]. Vũ Duy Vĩnh, “Mô phỏng quá trình lan truyền dầu khi xảy ra sự cố tràn dầu tại khu vực cửa sông Bạch Đằng”, Tuyển tập báo cáo Hội nghị 28Bet28 link nhà cái Thể thao 28bet, tr.290-298, 2007.

[15]. Kiều, Trần Duy, “Nghiên cứu mô phỏng nguy cơ lan truyền vệt dầu trong sự cố tràn dầu trên vùng biển Phú Quốc”, Tạp chí Khoa học Kỹ thuật Thủy lợi và Môi trường, Số 52: 44-51, 2016.

[16]. Trần Anh Tú, Lê Đức Cường, “Mô phỏng một số kịch bản tràn dầu khu vực đảo Cồn Cỏ”, Tạp chí Khoa học và Công nghệ Biển, Tập 14, Số 2: 187-194, 2014.

[17]. Nguyễn Đình Dương, “Ô nhiễm dầu trên biển Việt Nam và biển Đông”, Báo cáo tổng hợp đề tài KC.09.22/ 06-10, 2008.

[18]. British Oceanigraphic Data Centre, https://www.gebco.net/data-products/gridded-bathymetry-data

[19]. MIKE 21 and MIKE 3 Flow Model FM, “Hydrodynamic and Transport Module, Scientific Documentation”, DHI, 2014.

[20]. MIKE 21/3 Oil Spill Model, “Oil Spill Template, Scientific Documentation”, DHI, 2014.

[21]. Garcia-Martinez R. and Flores Tovar H., “Computer modeling of oil spill trajectories with a high accuracy method”, Spill Sci. Technol. Bull., 5 (5