Mô hình Delft3D áp dụng vào mô phỏng quá trình vận chuyển, lan truyền của trầm tích lơ lửng và muội than tại khu vực vịnh Hạ Long cho kết quả tương đối phù hợp với một số nghiên cứu thực nghiệm trước đây. Trong mùa mưa, trầm tích lơ lửng tập trung chủ yếu ở khu vực phía tây, tây nam đảo Cát Bà và ven bờ từ Cửa Lục trở xuống phía nam với hàm lượng 50–130 g/m3.
Vietnam Journal of Marine Science and Technology; Vol 19, No 2; 2019: 233–241 DOI: https://doi.org/10.15625/1859-3097/19/2/9729 https://www.vjs.ac.vn/index.php/jmst Simulation of suspended sediment and black carbon transport in surface water layer of Ha Long bay Le Duc Cuong*, Nguyen Van Thao Institute of Marine Environment and Resources, VAST, Vietnam * E-mail: cuongld@imer.vast.vn Received: January 2018; Accepted: 26 April 2018 ©2019 Vietnam Academy of Science and Technology (VAST) Abstract Delft3D model employed to simulate the distribution and transport of suspended sediment and black carbon in Ha Long bay shows outcomes meeting with results from previous experiment studies In the rainy season, suspended matter in surface layer is mainly in waters of western and southwestern Cat Ba island regions, and from Cua Luc toward the south nearshore areas with concentration of 50–130 g/m3 The concentration of suspended setdiment in the waters from Cua Luc to the north nearshore area is from 20 g/m3 to 50 g/m3 and that of offshore areas is 2–20 g/m3 In the dry season, the average concentrations of suspended matter are lower, approximately 110–150 g/m3 compared to the rainy season In the rainy season, the total particulate carbon in surface layer is 0.0016–0.0028 kg/m3 and in the dry season, it ranges from 0.0001–0.005 kg/m3 Keywords: Black carbon, Delft3D model, Ha Long bay Citation: Le Duc Cuong, Nguyen Van Thao, 2019 Simulation of suspended sediment and black carbon transport in surface water layer of Ha Long bay Vietnam Journal of Marine Science and Technology, 19(2), 233–241 233 Tạp chí Khoa học Công nghệ Biển, Tập 19, Số 2; 2019: 233–241 DOI: https://doi.org/10.15625/1859-3097/19/2/9729 https://www.vjs.ac.vn/index.php/jmst Mơ lan truyền trầm tích lơ lửng muội than lớp nƣớc bề mặt vịnh Hạ Long Lê Đức Cƣờng*, Nguyễn Văn Thảo Viện Tài nguyên Môi trường biển, Viện Hàn lâm Khoa học Công nghệ Việt Nam, Việt Nam * E-mail: cuongld@imer.vast.vn Nhận bài: 2-1-2018; Chấp nhận đăng: 26-4-2018 Tóm tắt Mơ hình Delft3D áp dụng vào mơ q trình vận chuyển, lan truyền trầm tích lơ lửng muội than khu vực vịnh Hạ Long cho kết tương đối phù hợp với số nghiên cứu thực nghiệm trước Trong mùa mưa, trầm tích lơ lửng tập trung chủ yếu khu vực phía tây, tây nam đảo Cát Bà ven bờ từ Cửa Lục trở xuống phía nam với hàm lượng 50–130 g/m3 Ở khu vực ven bờ từ Cửa Lục ngược lên phía bắc, hàm lượng trầm tích lơ lửng khoảng 20–50 g/m3, khu vực khơi khoảng 2–20 g/m3 Trong mùa khơ, hàm lượng trầm tích lơ lửng trung bình mùa thấp khoảng 110–150 g/m3 so với mùa mưa Hàm lượng hàm lượng muội than dạng hạt tổng số (BC) lớp nước bề mặt mùa khô dao động khoảng 0,0001–0,003 kg/m3 mùa mưa khoảng 0,0016–0,005 kg/m3 Từ khóa: Muội than, mơ hình Delft3D, vịnh Hạ Long MỞ ĐẦU Khu vực vịnh Hạ Long có địa hình đường bờ biển phức tạp, biên độ thủy triều thuộc loại lớn nước Dòng chảy tổng hợp dòng triều, dòng chảy gió dòng chảy sơng, dòng triều có vai trò chủ đạo định đến xu hướng vận chuyển trầm tích lơ lửng (TTLL) Do dòng triều có độ lớn biến đổi khơng đáng kể theo chiều sâu cột nước, vận chuyển TTLL lớp nước mặt tầng bên có xu khơng khác biệt nhiều Muội than hay gọi cacbon đen dạng sol khí, phát thải từ nguồn: cháy rừng (42%), đốt nhiên liệu sinh học theo phương pháp truyền thống (18%), giao thông đường (14%), giao thông khác (10%), công nghiệp máy phát (10%), chất đốt từ dân cư kỹ thuật truyền thống (6%) Muội than khí thường lắng đọng dạng khô 234 ướt đất liền đại dương Ở môi trường nước, chúng hợp phần TTLL Muội than tồn dạng: dạng hạt tổng số, hữu dạng hạt nguyên tố dạng hạt Ở Việt Nam nay, nghiên cứu có liên quan đến muội than môi trường nước chưa nhiều tập trung vào mơi trường khơng khí Vịnh Hạ Long nơi hội tụ đầy đủ nguồn phát thải muội than điển hình vùng ven biển Việt Nam Do đó, vùng vịnh lựa chọn vùng điển hình để tiến hành nghiên cứu vận chuyển muội than dạng hạt tổng số (BC) môi trường nước Nghiên cứu đánh giá trình vận chuyển, phát tán TTLL BC không gian theo thời gian thực thông qua việc áp dụng mơ hình thủy động lực ba chiều Delft3D vào mơ Qua tạo sở để hiệu chỉnh so sánh với phương pháp phân tích khác, góp phần xây dựng sở khoa Mơ lan truyền trầm tích lơ lửng học cho nghiên cứu có liên quan cho việc bảo vệ tài nguyên môi trường biển phát triển bền vững TÀI LIỆU VÀ PHƢƠNG PHÁP Tài liệu Số liệu độ sâu đường bờ vùng biển vịnh Hạ Long số hóa từ đồ Cục Đo đạc Bản đồ xuất năm 2005 Độ sâu vùng biển phía ngồi sử dụng từ sở liệu GEBCO-1/8 Đây số liệu địa hình có độ phân dải 0,5 phút xử lý từ ảnh vệ tinh kết hợp với số liệu đo sâu [1] Số liệu khí tượng gồm số liệu gió quan trắc nhiều năm Trạm Khí tượng Hải văn Môi trường Bãi Cháy thu thập xử lý (tần suất đo h/lần) tổng hợp tài liệu từ cơng trình cơng bố nước [2–4] Số liệu mực nước dùng để hiệu chỉnh mơ hình kết đo đạc mực nước (1 h/lần) Trung tâm Khí tượng thủy văn Quốc gia Các số điều hòa thủy triều biên phía ngồi khơi trích xuất từ sở liệu FES2004 [5] Số liệu nhiệt độ độ muối nước biển vùng biển vịnh Hạ Long vùng biển phía ngồi thu thập từ kết khảo sát hai mùa nghiên cứu liên quan khu vực Số liệu nhiệt độ độ muối nước biển vùng biển xa bờ thu thập từ sở liệu WOA13 với độ phân giải 0,25o [6] Số liệu dao động mực nước, dòng chảy hàm lượng TTLL kết đo đạc nhiệm vụ : “Nghiên cứu phương pháp phân tích, đánh giá giám sát chất lượng nước ven bờ tư liệu viễn thám độ phân giải cao độ phân giải trung bình, đa thời gian; Áp dụng thử nghiệm cho ảnh vệ tinh VNREDSat-1”, 2014–2015; Các số liệu lưu lượng nước sông thu thập, xử lý để phục vụ hiệu chỉnh kiểm chứng độ tin cậy mơ hình thủy động lực Số liệu muội than đo đạc định kỳ hàng tháng từ 10/2012 đến 10/2013 Các trạm khảo sát bao gồm 16 trạm mặt rộng trạm liên tục đại diện khu vực cửa sơng, trạm liên tục phía ngồi khơi tiến hành đo đạc, phân tích chuyến khảo sát (mùa mưa mùa khô) nhiệm vụ hợp tác quốc tế „„Nghiên cứu ảnh hưởng muội than (black carbon) tới hệ sinh vật môi trường biển vịnh Hạ Long nhằm đề xuất giải pháp bảo vệ môi trường biển‟‟ [7] thực Phƣơng pháp Mơ hình Delft3D [8] sử dụng mơ với mơ đun Delft3D-Waq để tiến hành tính tốn vận chuyển TTLL Black Carbon kết mô mô đun (Delft3D-Wave) kết hợp thủy động lực (Delft3D-Flow), áp dụng kỹ thuật lưới lồng (nesting) Để tiến hành phục vụ tính tốn trường thủy động lực cho khu vực vịnh Hạ Long, trường thủy động lực vùng phía ngồi khơi mơ phỏng, kết sau trích xuất làm đầu vào cho mơ hình thủy động lực vùng vịnh Hạ Long Q trình vận chuyển muội than mơi trường nước vịnh mô cách riêng biệt mô hình Delf3D với thơng số đầu vào thiết lập dựa đặc điểm tỷ trọng hệ số nhớt sau: Đường kính hạt 500–50.000 microns, tỷ trọng 0,00313–0,004119 kg/m3 [9] độ nhớt 12.800–192.600 cps [10] Sử dụng phương pháp đánh giá sai số bình phương trung bình RMSE [12] (Root Mean Square Error) làm tiêu để đánh giá độ xác mơ hình thơng qua kiểm chứng mực nước từ kết mơ hình quan trắc Thiết lập miền tính Phạm vi tính tốn khu vực phía ngồi khơi giới hạn từ vĩ độ 19o51‟36,83”N đến 21o29‟28,14”N, từ kinh độ 106o5‟46,54”E đến 108o7‟36,68”E Số điểm lưới tính tốn vùng phía ngồi (M, N) = 162, 267 với kích thước bước lưới nhỏ 300 m lớn km Vùng tính vịnh Hạ Long giới hạn phạm vi 80 km theo phương vĩ tuyến (đông - tây) 50 km theo phương kinh tuyến 50 km (nam - bắc) Với số điểm lưới (M, N) = 625 × 435, bước lưới nhỏ 150 m lớn 500 m (hình 1) 235 Lê Đức Cường, Nguyễn Văn Thảo Hình Trường độ sâu phạm vi tính tốn mơ hình Điều kiện ban đầu điều kiện biên Tại thời điểm ban đầu, vận tốc dòng chảy mực nước gán “0” Biên mở phía sơng: Gồm lưu lượng nước sơng trung bình qua mặt cắt mùa mưa mùa khô cửa Bạch Đằng, sông Giá, sông Gianh, sông Hốt, cửa Lục Tiên Yên Sử dụng số liệu nhiệt độ, độ muối trung bình mùa từ kết thu thập số liệu khảo sát đề tài, dự án Viện Tài nguyên Môi trường biển thực Biển mở phía biển: Sử dụng sóng triều (O1, K1, Q1, P1, N2, M2, S2, K2) trích xuất từ FES2004 [5] áp dụng cho mơ hình vùng phía ngồi, sau trích xuất mực nước làm đầu vào cho vùng vịnh Hạ Long Số liệu nhiệt độ, độ muối lấy từ số liệu đo đạc, thu thập kết hợp với số liệu lấy từ đề tài, dự án Viện Tài nguyên Môi trường biển thực Thời gian tính tốn 45 ngày cho kịch bản, thời gian tính tốn để mơ đạt độ ổn định 15 ngày Hiệu chỉnh kiểm chứng mơ hình Để đánh giá hiệu chỉnh cho mơ hình khu vực vịnh Hạ Long, chúng tơi sử dụng kết tính tốn mực nước mơ hình so sánh với mực nước quan trắc cho mùa mưa mùa khô Sau lần hiệu chỉnh cuối, sai số bình phương trung bình (RMSE) kết tính tốn có giá trị 0,19 cho mùa mưa 0,21 cho mùa khơ (hình 2) Sai số chấp nhận điều kiện địa hình khu vực tính phức tạp biên độ dao động mực nước lớn khu vực vịnh Hạ Long Dao động mực nước tính tốn mơ hình với kết quan trắc vẽ hình cho thấy phù hợp tương đối pha triều độ lớn 236 Hình So sánh mực nước tính tốn từ mơ hình quan trắc Hình Hàm lượng TTLL tầng mặt mơ hình với quan trắc khu vực Cửa Lục ngày đêm, mùa mưa Mơ lan truyền trầm tích lơ lửng Hình Hàm lượng BC tầng mặt mơ hình quan trắc khu vực Cửa Lục ngày đêm, mùa khô Qua nhiều lần hiệu chỉnh, kết mơ hàm lượng trầm tích lơ lửng muội than dạng hạt tổng số vị trí Cửa Lục (hình hình 4) thể quy luật phù hợp kết tính tốn từ mơ hình số liệu đo đạc máy Compact CTD phân tích ảnh viễn thám KẾT QUẢ MƠ PHỎNG Kết mơ trƣờng dòng chảy Trong mùa khơ, khu vực biển phía ngồi vịnh Hạ Long, dòng chảy có hướng nam tây nam, nam đông nam chiếm ưu với tốc độ khoảng 50–80 cm/s Đôi mạnh hơn, khoảng 120 cm/s tần suất nhỏ, khoảng 4–5% Ngoài ra, mùa xuất dòng chảy có hướng bắc, tây bắc với tốc độ tương đương có hướng ngược lại Trong mùa dòng chảy có hướng phân tán mạnh so với mùa mưa, đặc biệt vào thời kỳ chuyển giao hai pha triều Vào mùa mưa, tác động nước lục địa sông đổ gây ảnh hưởng đến hoàn lưu nước, hướng thịnh hành dòng chảy nam tây nam tây nam hướng ngược với bắc đơng bắc đông bắc chủ yếu Ở số cửa lạch hẹp quan sát thấy triều lưu chảy mạnh với tốc độ 90–100 cm/s Ở luồng Cửa Ông, triều lưu mạnh theo trục luồng Tại Cửa Lục triều lưu chảy mạnh triều lưu chảy vào, đạt tốc độ 80–90 cm/s, pha triều xuống lớn m/s (hình 5) (a) Trường dòng chảy pha triều lên, kỳ nước kém, mùa khơ (b) Trường dòng chảy pha triều xuống, kỳ nước kém, mùa khô (c) Trường dòng chảy pha triều lên, kỳ nước cường, mùa khơ (d) Trường dòng chảy pha triều xuống, kỳ nước cường, mùa khơ Hình Kết mơ trường thủy động lực vùng vịnh Hạ Long 237 Lê Đức Cường, Nguyễn Văn Thảo TTLL nhỏ Ở khu vực khác hàm lượng TTLL có giá trị tương đối nhỏ nằm khoảng 1,5–2,5 g/m3 (hình 6) Trong mùa mưa, hàm lượng TTLL lớp nước bề mặt với hàm lượng cao (50–130 g/m3) tập trung chủ yếu khu vực phía tây, tây nam đảo Cát Bà, ven bờ từ Cửa Lục trở xuống phía nam khu vực khơi ven bờ Bãi Cháy ngược lên phía bắc Vùng ven bờ từ Cửa Lục ngược lên phía bắc, hàm lượng TTLL khoảng 20–50 g/m3, khu vực khơi khoảng 2–20 g/m3 (hình 6) Kết mơ TTLL Vào mùa khơ, hàm lượng TTLL khu vực phía tây, tây nam đảo Cát Bà ven bờ từ Cửa Lục trở xuống phía nam khoảng 60–80 g/m3 có xu hướng phát tán mở rộng so với mùa mưa Hàm lượng TTLL trung bình mùa thấp hơn, khoảng 110–150 g/m3 so với mùa mưa Trong thời kỳ nước kém, phát tán TTLL phía phía ngồi khơi khơng có khác biệt nhiều pha triều lên xuống, vùng nước có hàm lượng TTLL cao tập trung chủ yếu khu vực cửa sông Dải ven bờ khu vực Bãi Cháy ngược lên phía bắc, hàm lượng (a) Pha triều lên, mùa mưa (b) Pha triều xuống, mùa mưa (d) Pha triều xuống, mùa khơ (c) Pha triều lên, mùa khơ Hình Kết mô TTLL (g/m ) tầng mặt khu vực vịnh Hạ Long Kết mô lan truyền BC Vào mùa mưa, vùng có BC với hàm lượng cao (0,002–0,003 kg/m3) có xu hướng phát tán 238 mở rộng so với mùa khô, mùa BC có điều kiện phát tán mạnh phía ngồi khơi, hàm lượng BC trung bình thấp Mơ lan truyền trầm tích lơ lửng khoảng 0,001–0,002 kg/m3 so với mùa khô, nhiên hàm lượng BC cực tiểu (0,0016 kg/m3) lại lớn so với mùa khô (0,001 kg/m3), tức biên độ dao động BC lớp nước bề mặt nhỏ mùa khô Vùng nước ven bờ có hàm lượng BC cao tập trung chủ yếu khu vực Cửa Lục Cửa Ông Trong thời kỳ nước kém, phát tán BC phía phía ngồi khơi khơng khác biệt nhiều pha triều lên pha triều xuống Trong thời kỳ nước với biên độ thủy triều nhỏ, phát tán BC phía ngồi khơi khơng có khác biệt nhiều pha triều lên pha triều xuống, phạm vi phát tán phía ngồi khơi hẹp so với thời kỳ nước cường, kỳ nước vùng nước ven bờ có hàm lượng BC cao tập trung chủ yếu khu vực Cửa Lục Cửa Ông Trong thời kỳ nước cường, vùng nước có hàm lượng BC cao mở rộng so với thời kỳ nước Trong mùa khô, hàm lượng hàm lượng BC lớp nước bề mặt dao động khoảng 0,001–0,003 kg/m3(hình 7) (a) Phân bố muội than, pha triều lên, mùa mưa (b) Phân bố muội than, pha triều xuống, mùa mưa (c) Phân bố muội than, pha triều lên, mùa khô (d) Phân bố muội than, pha triều xuống, mùa khơ Hình Phân bố BC (kg/m3) khu vực vịnh Hạ Long Trong thời kỳ nước kém, hàm lượng BC biến đổi biên độ dao động khoảng từ 0,001–0,005 kg/m3 Tuy nhiên, hàm lượng BC khu vực ven bờ phía đơng bắc phía tây vịnh Hạ Long có biến động mạnh so với khu vực phía đơng phía nam Trong pha triều lên, ảnh hưởng khối nước biển hàm lượng BC có xu tập trung cao khu vực cửa sông dọc bờ Trong pha triều xuống, BC lớp nước bề mặt có điều kiện mở rộng phía ngồi khơi phía nam nam đông nam Hàm lượng cao 239 Lê Đức Cường, Nguyễn Văn Thảo tập trung khu vực Cửa Lục Cửa Ông, đạt 0,005 kg/m3, khu vực vịnh Hạ Long, hàm lượng BC biến đổi khoảng 0,001–0,0015 kg/m3 THẢO LUẬN Qua kết tính tốn mơ ta thấy, dòng triều định đến hàm lượng trình vận chuyển TTLL BC nước Trong mùa mưa, TTLL có hàm lượng trung bình cao so với mùa khơ, vùng nước có hàm lượng TTLL cao chủ yếu tập khu vực phía tây, tây nam đảo Cát Bà ven bờ từ khu vực Cửa Lục xuống phía nam, lưu lượng nước sông lớn nên phạm vi phân bố TTLL phía ngồi khơi mạnh so với mùa khơ, đặc biệt pha triều xuống Nhìn chung, hàm lượng trung bình TTLL mùa mưa lớn mùa khơ, vùng nước có hàm lượng TTLL cao chủ yếu tập trung vùng ven bờ Cửa Lục xuống phía nam nơi có cửa sơng nhỏ đổ biển Ở khu vực phía tây nam đảo Cát Bà, chịu tác động lưu lượng nước sông Cấm Bạch Đằng, hàm lượng TTLL khu vực đạt giá trị cao Đặc biệt, vùng nước có hàm lượng TTLL cao phát tán mạnh mẽ pha triều xuống, phần TTLL theo khối nước sông di chuyển sang vịnh Hạ Long làm cho hàm lượng TTLL xung quanh khu vực phía tây tây nam đảo Cát Bà tăng cao Theo kết nghiên cứu Trần Đức Thạnh nnk., (2007), trầm tích lơ lửng khu vực có liên quan nhiều đến dao động mực nước biến động rõ rệt theo mùa.Vào mùa mưa, hàm lượng trầm tích lơ lửng nước trung bình 27,6 g/m3, biến đổi khoảng 5–60 g/m3, thường 20–40 g/m3, đạt đến 40–60 g/m3 khu vực Bãi Cháy khu luồng tàu chạy từ Cửa Lục tới Hòn Một Khu vực có hàm lượng trầm tích lơ lửng cao với hàm lượng từ 25–40 g/m3 phía tây nam đảo Tuần Châu, khu vực cửa sơng (phía vịnh Cửa Lục) khu vực cửa sơng Mơng Dương (Cửa Ơng) Các khu vực lại có hàm lượng trầm tích lơ lửng nước hầu hết nhỏ 25 g/m3 Phân bố khơng gian khu vực có hàm lượng trầm tích lơ lửng cao vịnh Hạ Long - Bái Tử Long biến đổi rõ rệt theo pha triều Khi triều lên khối nước có hàm lượng trầm tích 240 lơ lửng cao bị thu hẹp đẩy sát lên phía pha triều xuống khối nước sơng có hàm lượng trầm tích lơ lửng cao có điều kiện mở rộng biển Vào mùa mưa, kết hợp nước từ sông chảy pha triều xuống làm cho khu vực có hàm lượng trầm tích lơ lửng cao mở rộng phía ngồi biển lớn pha triều xuống mùa khô [12] Như vậy, kết nghiên cứu tương đồng với kết nghiên cứu trên, phù hợp hàm lượng TTLL trung bình nước khoảng biến động, không gian phân bố, quy luật vận chuyển phát tán TTLL điều kiện thủy động lực Muội than dạng trầm tích lơ lửng, phát tán BC tuân theo quy luật vận chuyển phát tán trầm tích lơ lửng Trong mùa mưa biến động hàm lượng BC theo không gian lớn so với mùa khô, nhiên hàm lượng BC trung bình mùa khơ lại lớn mùa mưa Nguyên nhân mùa khô ảnh hưởng hướng gió đơng bắc áp suất khí mùa khô cao so với mùa mưa, điều làm cho muội than phát tán phía ngồi biển mùa mưa, lượng muội than từ khí vào môi trường nước diễn mạnh Trong thời kỳ nước kém, biên độ thủy triều nhỏ nên phát tán BC phía phía ngồi khơi khơng có khác biệt nhiều pha triều lên pha triều xuống, phạm vi phát tán phía ngồi khơi hẹp hợn so với thời kỳ nước cường, kỳ nước vùng nước ven bờ có hàm lượng BC cao tập trung chủ yếu khu vực Cửa Lục Cửa Ông Trong thời kỳ nước cường, vùng nước có hàm lượng BC cao mở rộng so với thời kỳ nước Tuy nhiên, theo chu kỳ ngày hàm lượng BC khu vực ven bờ phía đơng bắc phía tây vịnh Hạ Long có biến động mạnh so với khu vực phía đơng phía nam Trong pha triều lên, ảnh hưởng khối nước biển, hàm lượng BC có xu tập trung với hàm lượng cao khu vực cửa sông dọc theo ven bờ Trong pha triều xuống, BC lớp nước bề mặt có điều kiện mở rộng phía ngồi khơi phía nam nam đơng nam Mơ lan truyền trầm tích lơ lửng TÀI LIỆU THAM KHẢO [1] Jones, M T., Weatherall, P., and Cramer, R N., 2009 User guide to the Centenary Edition of the GEBCO Digital Atlas and its data sets Natural environment research council [2] Nguyễn Thế Tưởng, 2000 Sổ tay tra cứu đặc trưng khí tượng thủy văn vùng thềm lục địa Việt Nam Trung Tâm Khí tượng Thủy văn biển - Tổng Cục Khí tượng Thủy văn Nxb Nông nghiệp [3] Trạm dự báo phục vụ Khí tượng Thủy văn Quảng Ninh, 1997 Đề tài khoa học “Số liệu khí tượng thủy văn tỉnh Quảng Ninh 50 năm 1946–1995”, Quảng Ninh [4] Số liệu điều kiện tự nhiên dùng xây dựng, Phần I Quy chuẩn xây dựng Việt Nam QCXDVN 02: 2008/BXD [5] Lyard, F., Lefevre, F., Letellier, T., and Francis, O., 2006 Modelling the global ocean tides: modern insights from FES2004 Ocean Dynamics, 56(5–6), 394–415 [6] World Ocean Atlas, National Geophysical Data Center - http://www.ngdc.noaa.gov [7] Chu Văn Thuộc, 2015 Nghiên cứu ảnh hưởng muội than (black carbon) tới hệ sinh vật môi trường biển vịnh Hạ Long [8] [9] [10] [11] [12] nhằm đề xuất giải pháp bảo vệ tài nguyên môi trường biển WL|Delft Hydraulics, 1999 Delft3D-FLOW User Manual Version 3.05, Delft3D - Wave User Manual, Delft3D-Waq User Manual Version 3.01 WL| Delft Hydraulics, Delft, Netherlands, 2010 Maiwada, A D., Adamu, A A., 2015 Analysis of Black Carbon Density from 30 kva Tata and 30 kva Perkins Generator Exhausts American Journal of Mechanical Engineering and Automation, 2(6), 59–62 Current Trends in The Global Plasticcizer Industry, 2001 Tarun Khemlani, Chemical Market Resources, Inc., USA Organisedby Rapra Techology Ltd The 7th International Plastics Additives and Modifiers Conference Conference Proceedings, 8th–9th October 2001 ISBN: 14 85957 2952 Chai, T., and Draxler, R R., 2014 Root mean square error (RMSE) or mean absolute error (MAE)? Geoscientific Model Development Discussions, 7, 1525–1534 Trần Đức Thạnh nnk., 2007 Sức tải môi trường vịnh Hạ Long - Bái Tử Long Nxb Khoa học Tự nhiên Công nghệ 241 ... https://www.vjs.ac.vn/index.php/jmst Mô lan truyền trầm tích lơ lửng muội than lớp nƣớc bề mặt vịnh Hạ Long Lê Đức Cƣờng*, Nguyễn Văn Thảo Viện Tài nguyên Môi trường biển, Viện Hàn lâm Khoa học... áp dụng vào mơ q trình vận chuyển, lan truyền trầm tích lơ lửng muội than khu vực vịnh Hạ Long cho kết tương đối phù hợp với số nghiên cứu thực nghiệm trước Trong mùa mưa, trầm tích lơ lửng tập... dạng hạt tổng số (BC) lớp nước bề mặt mùa khô dao động khoảng 0,0001–0,003 kg/m3 mùa mưa khoảng 0,0016–0,005 kg/m3 Từ khóa: Muội than, mơ hình Delft3D, vịnh Hạ Long MỞ ĐẦU Khu vực vịnh Hạ Long