Đánh giá đặc trưng trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng Trần Anh Tú Trường Đại học Khoa học Tự nhiên Luận văn ThS Chuyên ngành: Hải dương học; Mã số: 60 44 97 Người hướng dẫn: GS TS Đinh Văn Ưu Năm bảo vệ: 2012 Abstract: Tổng quan vấn đề nghiên cứu điều kiện tự nhiên khu vực Hải Phòng Nghiên cứu tài liệu phương pháp qua địa hình, khí tượng, thủy hải văn, trầm tích lơ lửng, mơ hình thủy động lực, mơ hình lan truyền trầm tích lơ lửng Đánh giá trạng trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng: phân bố TTLL (Trầm tích lơ lửng) theo thời gian, đặc điểm TTLL khu vực sơng Hải Phịng, đặc điểm TTLL khu vực xa bờ Hải Phịng Mơ trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng mơ hình delft3d Keywords: Hải dương học; Trầm tích; Cửa sơng ven biển; Hải Phòng Content Các tượng sa bồi luồng cảng, cửa sơng, xói lở-bồi tụ bờ biển, độ đục nước gia tăng làm ảnh hưởng đến chất lượng nước bãi tắm, khu nuôi trồng thủy sản liên quan đến trầm tích lơ lửng (TTLL) Ngồi ra, khu vực có giá trị hàm lượng TTLL cao làm ảnh hưởng tới tầm nhìn xuyên suốt khối nước, quang hợp thực vật sống lồi sinh vật mơi trường nước Thành phố cảng Hải Phịng năm có đóng góp quan trọng hai ngành kinh tế đặc trưng dịch vụ cảng biển du lịch Tuy nhiên, đặc thù địa lý vùng cửa sông, khu vực ven biển thành phố Hải Phòng chịu ảnh hưởng nặng nề dòng vật chất từ lục địa đưa qua hệ thống sơng Thái Bình, sơng Hồng Trong dịng vật chất đó, dịng trầm tích lơ lửng có chế phức tạp nguyên nhân tự nhiên (dịng chảy, sóng, xói lở bờ) người (nạo vét luồng, khai hoang lấn biển, phá rừng ngập mặn nuôi trồng thủy sản) gây Sa bồi luồng vào cảng Hải Phịng có xu hướng gia tăng kèm việc chi phí cho việc nạo vét luồng lạch tốn Theo thống kê Cảng vụ Hàng hải Hải Phòng (2005), khối lượng nạo vét luồng vào cảng Hải Phòng năm 2003 2004 tương ứng 2.394.000m3 2.854.000m3 [29] Mặt khác, trình nạo vét luồng vào cảng diễn thường xuyên khiến cho bùn cát vật chất ô nhiễm lắng xuống lại bị đưa lên, hòa tan nước làm gia tăng nguy gây ô nhiễm đến môi trường nước hệ sinh thái xung quanh [2] Mặt khác, dòng vật chất làm ảnh hưởng đến chất lượng bãi tắm Đồ Sơn khu nuôi trồng hải sản đảo Cát Bà làm giảm hiệu đáng kể mặt kinh tế Ngồi thành phố Hải Phịng có kế hoạch thực dự án đê quai lấn biển phục vụ xây dựng Sân bay Quốc tế vùng ven bờ Tiên Lãng Việc nhiều làm thay đổi chế dịng chảy, vận chuyển trầm tích lơ lửng sơng Văn Úc Thái Bình nói riêng vùng cửa sơng ven bờ Hải Phịng nói chung [5] Bởi vậy, việc đánh giá TTLL vùng cửa sơng ven biển Hải Phịng điều cần thiết Với mục tiêu luận văn mô trạng trầm tích lơ lửng theo mùa khu vực cửa sơng ven biển Hải Phòng Học viên đánh giá lựa chọn số liệu thu thập từ tài liệu nghiên cứu từ trước đến trầm tích lơ lửng làm số liệu đầu vào cho mơ hình tính Nội dung luận văn trình bày thành 04 chương: Chƣơng 1: Tổng quan vấn đề nghiên cứu điều kiện tự nhiên khu vực Phòng Hải Chƣơng 2: Tài liệu phương pháp Chƣơng 3: Đánh giá trạng trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng Chƣơng 4: Mơ trầm tích lơ lửng khu vực cửa sơng ven biển Phịng mơ hình delft3d Hải CHƢƠNG I TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU VÀ ĐIỀU KIỆN TỰ NHIÊN KHU VỰC HẢI PHÒNG I.1 Tổng quan vấn đề nghiên cứu I.1.1 Tình hình nghiên cứu ngồi nước Nghiên cứu phân bố trầm tích lơ lửng vùng cửa sơng ven biển nhà khoa học nước quan tâm từ hàng trăm năm đạt nhiều thành tựu quan trọng Các kết nghiên cứu áp dụng phục vụ cho việc bảo vệ cơng trình ven bờ phát triển mơi trường bền vững Những nghiên cứu lý thuyết trầm tích (vận chuyển) đáng kể cơng trình H.A Einstein (1950), Krone Partheniades (1962, 1968), E.W BijJker (1967, 1971), Leo C Van Rijn (1993), J.W Vander Meer (1990), Richard Soulsby (1997) Theo hướng mơ hình hóa có loại mơ hình: Mơ hình vật lý mơ hình toán Những trung tâm, viện hàng đầu nghiên cứu, tính tốn, dự báo q trình thủythạch động lực, kể đến là: Trung tâm Thủy lực Hà Lan (Delft Hydraulics) với phần mềm DELFT3D, UNIBEST; Viện Thủy lực Đan Mạch (Danish Hydraulic Instiute-DHI) tiếng với phần mềm: MIKE 21, MIKE 3, hay Trung tâm Nghiên cứu Cơng trình Ven bờ thuộc Qn đội Mỹ (Coastal Engineering Research CenterCERC) có mơ hình GENESIS, SBEACH; mơ hình TELEMAC (Pháp), ứng dụng tốt cho tính tốn dịng chảy, sóng, vận chuyển bùn cát, biến động địa hình đáy biển, đường bờ, bồi lấp cửa sơng Ngồi ra, số mơ hình có mã nguồn mở COHERENS (Bỉ), SHYFEM (Italia) phát triển khơng phần mềm nói I.1.2 Tình hình nghiên cứu nước Vấn đề nghiên cứu trầm tích lơ lửng q trình động lực mơ hình hóa Việt Nam bắt đầu phát triển từ khoảng đầu năm 1980 [10] tiếp tục trọng Chương trình Biển KT.03 (1991-1995); KHCN.06 (1996-2000); TTLL liên quan đến xói lở bờ biển cịn đặt khn khổ đề tài độc lập cấp nhà nước chương trình biển giai đoạn 2001-2005 Ngoài nhiều đề tài, dự án liên quan đến TTLL thực cấp Các tác giả Trần Hồng Thái [18], Đinh Văn Ưu [28], Nguyễn Thọ Sáo [16] ứng dụng phát triển số phương pháp mơ hình tính tốn động lực vận chuyển trầm tích cho số vùng cửa sơng ven biển (Cửa ƠngQuảng Ninh, Cửa Tùng-Quảng Trị, Hải Phòng) Việt Nam Các tác giả thuộc Viện Hải dương học, Nha Trang, tiêu biểu Bùi Hồng Long nghiên cứu vùng Phan Rí, Hàm Tiến, Phước Thể với mục tiêu cung cấp thông số kỹ thuật, đưa phương án thiết kế thi cơng đê, kè chống xói lở [14] Trong năm gần Nhà nước cho triển khai nhiều đề tài nghiên cứu, ứng dụng trình thủy-thạch động lực xói lở, bồi tụ chẳng hạn như: Ứng dụng mơ hình WAM, STWAVE để dự báo sóng đề tài KC.09.04; Các đề tài KT.03.14, KHCN.06.08 (1996-2000), KC.09.05 (2001-2005) tiến hành nghiên cứu, dự báo trình xói lở-bồi tụ bờ biển cửa sơng Việt Nam; dự án Việt Nam-Thụy Điển (2004-2007): nghiên cứu xói lở bờ biển Hải Hậu, Nam Định, dự án ứng dụng nhiều mơ hình sóng, vận chuyển bồi tích; Đề tài cấp nhà nước KHCN-0610 ”Cơ sở khoa học đặc trưng đới bờ phục vụ yêu cầu xây dựng cơng trình biển ven bờ” Viện Cơ học chủ trì Các đề tài ngồi việc đo đạc thực địa xây dựng áp dụng mơ hình nhằm tính tốn q trình sóng, dịng chảy, vận chuyển trầm tích, biến đổi địa hình bãi, đường bờ, nhằm lý giải nguyên nhân gây tác động môi trường Các tác giả thuộc Viện TN&MTB ứng dụng mơ hình DELFT3D để nghiên cứu vấn đề có liên quan đến TTLL khu vực khác Quảng Ninh [20], Hải Phịng [1], Thái Bình Nam Định [23] Các nghiên cứu giúp cho nhà quản lý địa phương nói có cách nhìn cách tổng thể mối quan hệ phát triển kinh tế - xã hội, quản lý tổng hợp dải ven bờ bảo vệ môi trường biển Vấn đề nghiên cứu trầm tích lơ lửng vùng cửa sơng Hải Phịng có số kết định [21, 26, 28] Tuy nhiên kết nhận định chuyên gia đặc trưng trầm tích lơ lửng phạm vi hẹp mà chưa có cách nhìn cách tổng quan phạm vi khơng gian biến đổi theo thời gian (hàng chục năm) I.2 Điều kiện tự nhiên khu vực nghiên cứu Thành phố Hải Phòng nằm bờ tây vịnh Bắc Bộ thuộc phía đơng vùng dun hải Bắc Bộ Chế độ khí hậu khu vực chia thành hai mùa: mùa đơng có đặc điểm lạnh, khơ mưa mùa hè có đặc điểm nóng, ẩm mưa nhiều Nhiệt độ khơng khí trung bình năm khu vực dao động khoảng từ 22,5-30,00C Lượng mưa trung bình nhiều năm vùng ven biển Hải Phòng lớn với giá trị từ 1.600 2.000mm Tuy nhiên, lượng mưa phân bố không mà chủ yếu tập trung vào tháng mùa hè, cao vào tháng đạt 200mm.Chế độ gió khu vực ven biển Hải Phịng thể rõ rệt ảnh hưởng hồn lưu chung khí thay đổi theo mùa Về mùa đông thịnh hành gió hướng bắc đơng bắc, vận tốc gió trung bình dao động khoảng 3,2-3,7m/s, mạnh đạt tới 25-30m/s Vào mùa hè, chế độ gió khu vực chịu chi phối hệ thống gió mùa tây nam, hướng gió chủ yếu đơng nam nam, tốc độ gió trung bình khoảng 3,5-4,0 m/s, cực đại đạt 20 - 25m/s Khu vực Hải Phịng giai đoạn 1945-2007, có 53 bão ảnh hưởng trực tiếp tỉnh/thành lân cận Khu vực nghiên cứu chịu tác động trực tiếp từ nguồn cung cấp nước trầm tích sơng chảy vào Các sông phần hạ lưu cuối trước đổ biển hệ thống sông Thái Bình, gồm sơng Bạch Đằng, Cấm, Lạch Tray, Văn Úc Thái Bình Thủy triều vùng ven biển Hải Phòng nhật triều nhất, với biên độ dao động mực nước từ - m Dòng chảy ven bờ Hải Phòng dòng chảy tổng hợp, gờ m thành phần dịng chảy triều, gió, sóng, dịng chảy sơng tính chất dịng chảy vịnh Bắc Bộ Dòng chảy khu vực cửa sông thường chịu ảnh hưởng mạnh sông, phía ngồi biển dịng chảy tầng mặt chủ yếu gió gây Trong trường gió đơng bắc, dịng chảy tầng mặt dọc bờ có hướng thống trị từ bắc xuống nam cịn trường gió đơng nam, dòng chảy tầng mặt dọc bờ thể nhiều hướng khác tùy thuộc vào cường đô ̣ dòng chảy sơng đ ịa hình bờ biển Sóng ven biển Hải Phịng chủ yếu sóng truyền từ khơi bị khúc xạ phân tán lượng ma sát đáy Độ cao sóng chiếm tần suất lớn khoảng 0,25 - 1,0m, chủ yếu sóng hướng đơng đơng bắc Sóng có độ cao từ 3-4m chiếm tần suất nhỏ 0,15% [1, 2, 12] Trầm tích lơ lửng khu vực cửa sông ven biển khu vực nghiên cứu nhiều nguồn cung cấp khác nguồn chủ yếu từ sơng đưa Hàm lượng trầm tích lơ lửng sông biến thiên khoảng rộng, từ 10 đến 1.000g/m3 năm CHƢƠNG II TÀI LIỆU VÀ PHƢƠNG PHÁP II.1 Tài liệu - Địa hình: số liệu độ sâu đường bờ khu vực ven biển thành phố Hải Phịng số hố từ đồ địa hình UTM tỷ lệ 1: 50.000 Cục Đo đạc Bản đồ xuất bản, đồ xuất năm gần với hệ tọa độ nhà nước VN2000 - Khí tượng: q trình tính tốn, sử dụng số liệu gió tháng năm 2010 tháng năm 2010 quan trắc trạm Khí tượng Thủy văn Hịn Dấu - Thủy hải văn: số liệu lưu lượng sông lấy từ chuỗi quan trắc từ trạm cố định sơng Hải Phịng, ngồi số liệu dòng chảy tham khảo từ kết khảo sát đánh giá đề tài [4, 23, 24, 25] Học viên sử dụng chương trình dự báo thuỷ triều Đài Thiên văn Quốc gia Nhật Bản NAO (National Astronomical Observatory) để tính cho vùng vịnh Bắc Bộ, sau trích xuất mực nước biên lỏng khu vực nghiên cứu để tính trường dịng chảy - Trầm tích lơ lửng: số liệu quan trắc phân tích trầm tích lơ lửng đề tài [1, 3, 4, 6, 11, 12] sử dụng để đánh giá trạng làm tư liệu tham khảo cho mô trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng II.2 Phƣơng pháp II.2.1 Mơ hình thủy động lực Cơ sở tốn học mơ hình thuỷ động lực giải phương trình Navier Stokes với chất lỏng không nén nước nông phương pháp xấp xỉ Boussinesq Sự biến đổi thành phần vận tốc thẳng đứng phương trình động lượng bỏ qua Với mơ hình chiều, thành phần vận tốc thẳng đứng tính tốn từ phương trình liên tục [30] Theo phương nằm ngang, tuỳ thuộc vào điều kiện cụ thể khu vực tính, hệ toạ độ chủ yếu sau sử dụng: - Hệ toạ độ Đề (Cartesian): (x,y) Dùng cho vùng tính có địa hình đường bờ đơn giản - Hệ toạ độ cong trực giao: (  , ) Sử dụng cho khu vực có địa hình đường bờ phức tạp vùng cửa sông, ven biển, vũng vịnh v.v - Hệ toạ độ cầu: (  ,  ) Áp dụng cho khu vực rộng lớn, trải dài nhiều kinh - vĩ độ khác Phương trình liên tục (viết hệ toạ độ cong trực giao):      (d   )U G  (d   )V G  1   Q t   G G G G (1) Với Q thể thêm vào hay nguồn nước, bốc mưa đơn vị diện tích: Q  H  (qin  qout )d  P  E 1 đó:  , lµ hệ số hệ toạ độ cong trực giao G , G hệ số chuyển đổi từ hệ toạ độ Đề Các sang hệ tạo độ cong trực giao d độ sâu điểm tính (độ sâu nước đường chuẩn (0 hải đồ))  mực nước điểm tính (mực nước đường chuẩn) U, V thành phần vận tốc theo hướng  , qin qout nguồn nước đưa vào đơn vị thể tích H độ sâu điểm tính (H=d +  ) P, E lượng mưa bốc Phương trình bảo tồn động lượng theo hướng   ( toạ độ cong trực giao):  G  G u u u v u  u uv v2       fv t  G  G  d    G G  G G  1   u  P  F   vv   M (d   )      G (2)  G  G v u v v v  v uv u2       fu t G  G  d    G G  G G   1   v  P  F   vv   M  (d   )      G (3) Phương trình viết cho thành phần vận tốc theo phương thẳng đứng:      (d   )u G  (d   )v G  1      H (qin  qout ) (4) t    G G G G II.2.2 Mơ hình lan truyền trầm tích lơ lửng Cơ sở tốn học mơ hình vận chuyển trầm tích phương trình lan truyền khuyếch tán vật chất:    C C   C    C    Dx  uxC    Dy  u yC    Dz  uzC    z t x  x y  y   z   Nếu tính nguồn đưa từ ngồi vào thì:    C   C C C       Dx  u xC    Dy  u y C    Dz  u z C   F (C , t )   z t x  x y z  y     (6) (7) Trong phương trình (6, 7) trên: Dx, Dy, Dz hệ số khuyếch tán theo phương x, y, z ; F(C, t) nguồn vật chất thêm vào đi; C: hàm lượng vật chất Đối với mơ hình lan truyền trầm tích lơ lửng, giả thiết hàm lượng trầm tích lơ lửng cột nước giảm xảy q trình lắng đọng trầm tích Mặt khác q trình xói (tái lơ lửng-resuspension) xảy hàm lượng trầm tích lơ lửng cột nước tăng lên Q trình lắng đọng trầm tích phụ thuộc vào ứng suất xung quanh bề mặt (ambient shear stress-Tau) ứng suất tới hạn cho trình lắng đọng (Taucrsed) Nếu ứng suất xung quanh thấp ứng suất tới hạn, diễn q trình lắng đọng trầm tích Dịng trầm tích lắng đọng (Sedimentation flux)= Psed x Vsed x (IM1) (g/m2/ngày) Trong : Vsed : vận tốc lắng đọng trầm tích IM1 : nồng độ vật chất vô Psed : khả lắng đọng tính theo cơng thức :  Tau Psed = max  0,1  sed  Taucr      (8) Ngược lại, q trình xói xảy ứng suất xung quanh cao ứng suất tới hạn cho q trình tái lơ lửng (Taucrres): Dịng tái lơ lửng (Ressuspension flux) = Pres x Zres (g.m2/ngày) Trong đó: Zres là: tỷ lệ tái lơ lửng ban đầu Pres là: khả tái lơ lửng tính theo cơng thức: res  Taucr  Pres = max  0,  Tau 1    (9) Trong đó: IM1 hàm lượng trầm tích lơ lửng; Psed- xác suất xảy trình lắng đọng trầm tích; Vsed- Vận tốc lắng đọng; Tau-ứng suất xung quanh; Taucrsed-ứng suất tới hạn cho trình lắng đọng trầm tích; Taucrres-ứng suất tới hạn cho q trình tái lơ lửng; Pres-xác suất xảy trình tái lơ lửng; Zres- tốc độ tái lơ lửng từ bề mặt đáy Ứng suất xung quanh (Tau) phụ thuộc vào trình động lực sóng, gió, dịng chảy, mực nước, độ nhám đáy CHƢƠNG III ĐÁNH GIÁ HIỆN TRẠNG TRẦM TÍCH LƠ LỬNG VỰC CỬA SƠNG VEN BIỂN HẢI PHỊNG KHU III.1 Phân bố TTLL theo thời gian Dựa chuỗi số liệu quan trắc giai đoạn 1996-2010 (15 năm), học viên phân tích đánh giá biến đổi TTLL khu vực cửa sông ven biển Hải Phịng: giá trị TTLL trung bình năm năm 1996, 1997, 2001, 2006, 2007, 2008, 2009 vượt từ 1,04 đến 2,54 lần cho phép so với QCVN10: 2008/BTNMT Trong 10 năm trở lại nồng độ TTLL trung bình khu vực Đồ Sơn cho thấy: tầng mặt có giá trị xấp xỉ GHCP tầng đáy vượt GHCP so với QCVN10: 2008/BTNMT So với khu vực khác ven bờ biển Bắc Bộ Trà Cổ, Cửa Lục, Sầm Sơn, Cửa Lị, khu vực nghiên cứu có giá trị TTLL trung bình lớn hơn, thấp khu vực cửa Ba Lạt Về giá trị trung bình tầng đáy lớn tầng mặt, giống xu chung khu vực khác III.2 Đặc điểm TTLL khu vực sơng Hải Phịng Theo [4] cho thấy TTLL trung bình ngày đêm sơng có giá trị tầng mặt lớn tầng đáy hai mùa (mưa khô) Tuy nhiên, khu vực cửa sông Thái Bình có xu ngược lại, tầng đáy có giá trị TTLL lớn tầng mặt Tại cửa sông Văn Úc giá trị TTLL trung bình ngày đêm tầng đáy tầng mặt chênh không đáng kể Giá trị TTLL trung bình ngày lớn 133 mg/l quan trắc tầng đáy khu vực sông Cấm nhỏ 18 mg/l tầng mặt khu vực sông Bạch Đằng Các sông (Bạch Đằng, Cấm, Lạch Tray) phía bắc bán đảo Đồ Sơn giá trị TTLL tầng đáy lớn tầng mặt, sơng khu vực phía nam bán đảo Đồ Sơn có đặc điểm ngược lại tầng mặt lớn tầng đáy Vào mùa mưa cho thấy giá trị trung bình tầng mặt khu vực sông Bạch Đằng cửa sơng Thái Bình lớn tầng đáy Tuy nhiên khu vực cửa sông Cấm, Lạch Tray Văn Úc có xu ngược lại, giá trị TTLL tầng đáy lớn tầng mặt Giá trị TTLL trung bình ngày lớn 427 mg/l quan trắc tầng đáy khu vực sông Cấm nhỏ 36 mg/l tầng đáy khu vực sơng Thái Bình III.3 Đặc điểm TTLL khu vực xa bờ Hải Phòng Trong khuôn khổ Dự án [12] hợp tác thành phố Hải Phịng thành phố BrestPháp có tiến hành thu mẫu phân tích TTLL khu vực xa cửa sơng ven bờ biển Hải Phịng Theo số liệu cho thấy hàm lượng TTLL nước thấp GHCP so với QCVN10: 2008/BTNMT Hàm lượng TTLL trung bình nước ven biển Hải Phịng đợt khảo sát năm 2009 dao động từ 28,05 - 99,1 mg/l mùa mưa từ 20,26 - 68,34mg/l mùa khô Hàm lượng TTLL nước biển năm 2009 khu vực Hải Phịng có xu hướng tăng cao nước tầng đáy, mùa mưa lớn mùa khô Xét theo mặt cắt, nhận thấy phân bố hàm lượng TTLL theo mặt cắt phức tạp, không tuân theo quy luật chung giảm từ bờ khơi Điều cho thấy vận chuyển trầm tích lơ lửng vùng ven biển Hải Phòng phức tạp, ngồi nguồn sơng đưa cịn chịu ảnh hưởng chế độ triều, chế độ dịng chảy sơng chế độ thuỷ động lực khác Theo bảng 3.2 cho thấy vào mùa mưa mặt cắt có giá trị TTLL trung bình lớn (62,14 mg/l) sau đến mặt cắt 5, mặt cắt 4, mặt cắt 1, mặt cắt (34,78 mg/l) Vào mùa khô mặt cắt có giá trị TTLL trung bình lớn (43,08 mg/l) sau đến mặt cắt 1, mặt cắt 4, mặt cắt 3, mặt cắt (26,41 mg/l) Khu vực phía bắc quần đảo Long Châu có giá trị TSS trung bình tầng mặt 26,15 mg/l trung bình tầng đáy 27,25 mg/l mùa mưa Khu vực bãi tắm Cát Cò I (đảo Cát Bà) có giá trị TSS tầng mặt tầng đáy 27,8 30,2 mg/l Phía tây nam đảo Cát Bà có giá trị TSS tầng mặt tầng đáy 37,6 38,3 mg/l Nhìn chung, giá trị hàm lượng TTLL khu vực xa bờ Hải Phòng thấp GHCP so với QCVN10: 2008/BTNMT Dựa số liệu quan trắc trung bình tầng đáy tầng mặt, tranh phân bố TTLL theo khơng gian mùa khơ khu vực trung tâm có giá trị 30 mg/l, phía đơng bắc có giá trị khoảng 35-40 mg/l, vùng cửa sông (Cấm Bạch Đằng) phía bắc mũi Đồ Sơn có giá trị khoảng 35-80 mg/l phía nam (cửa sơng Văn Úc, Thái Bình) mũi Đồ Sơn có giá trị khoảng 35-60 mg/l Trong mùa mưa, khu vực trung tâm có giá trị 30-35 mg/l, phía đơng bắc có giá trị khoảng 40-45 mg/l, khu vực cửa sông Cấm-Bạch Đằng có giá trị khoảng 45-120 mg/l (Hình 3.10) Đây giá trị làm sở cho việc tiến hành mô TTLL (chương IV) khu vực nghiên cứu theo đặc trưng mùa CHƢƠNG IV MƠ PHỎNG TRẦM TÍCH LƠ LỬNG KHU VỰC SƠNG VEN BIỂN HẢI PHỊNG BẰNG MƠ HÌNH DELFT3D CỬA IV Triển khai mơ hình thủy động lực Phạm vi miền tính mơ hình  Mô dao động mực nước khu vực vịnh Bắc Bộ: - Số liệu đường bờ sử dụng từ đồ địa hình UTM tỷ lệ 1: 25000 (Cục Đo đạc Bản đồ), đồ với hệ tọa độ nhà nước VN-2000 Những đồ qt vào máy tính, số hố xử lý phần mềm Acview, MapInfo - Số liệu độ sâu lấy từ nguồn sở liệu địa hình ETOPO5 (Earth Topography - Minute) Trung tâm Tư liệu Địa vật lí Quốc gia Mỹ NGDC (National Geophysical Data Center) GEBCO -1 (General Bathymetric Chart of the Ocean (GEBCO) one minute) Trung tâm tư liệu hải dương học vương quốc Anh (British Oceanographic Data Centre-BODC) - Số liệu đầu vào tính tốn cho khu vực vịnh Bắc Bộ lấy từ số điều hịa tồn cầu FES2004 dự án Topex/ Poseidon với độ phân giải 1/8 độ tham khảo tài liệu Takashi TAKANEZAWA, 2000 [17]  Mô trường dịng chảy khu vực nghiên cứu: Mơ hình thuỷ động lực cho khu vực cửa sông ven biển Hải Phòng sử dụng hệ lưới cong trực giao Phạm vi vùng tính mơ hình bao gồm vùng nước cửa sông Bạch Đằng, Cấm, Lạch Tray phần phía ngồi cửa sơng mở rộng phía ngồi biển đến độ sâu khoảng 11 đến 12 mét Miền tính chia thành 185 x 356 lưới (M = 185, N = 356), trải dài từ 106o68E - 106o98 E, 20o65N - 20o948N (Hình 4.1b) Kích thước ô lưới biến đổi từ 48,24 đến 158,3 mét Thời gian tính tốn - Mùa khơ (từ ngày 01 đến ngày 31 tháng năm 2010) - Mùa mưa (từ ngày 01 đến ngày 31 tháng năm 2010) Bước thời gian chạy mơ hình thủy động lực 0,5 phút Hiệu chỉnh kiểm nghiệm kết mơ hình Sau lần hiệu chỉnh cuối, kết tính tốn cho thấy sai số bình phương trung bình mực nước tính từ mơ hình bảng dự báo thủy triều mùa mưa khô 0,29 m 0,37 m Dao động mực nước tính từ mơ hình với mực nước thực đo cho thấy có phù hợp định pha triều độ lớn (Hình 4.1) Các kết sau lần hiệu chỉnh cuối cho thấy phù hợp tương đối tính tốn số liệu quan trắc thực tế mực nước (m) mơ hình thực đo 3.5 2.5 1.5 0.5 thời gian 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 a) Tháng năm 2010 mơ hình mực nước (m) thực đo 3.5 2.5 1.5 0.5 thời gian 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 b) Tháng năm 2010 Hình 4.1 Đƣờng trình mực nƣớc thực đo kết tính từ mơ hình Trạm Hịn Dấu Vận tốc dịng chảy tính tốn quan trắc có khác biệt, nhiên số liệu quan trắc tiến hành đo tầng mặt (cách mặt 0,5 mét), kết từ mơ hình tính cho trung bình cột nước Đây hạn chế q trình tính tốn nên nhiều ảnh hưởng đến kết nghiên cứu IV.2 Triển khai mơ hình lan truyền trầm tích lơ lửng Phạm vi miền tính mơ hình Mơ hình mơ lan truyền trầm tích lơ lửng dùng kết mơ hình thủy động lực làm điều kiện cho việc tính tốn Vì vậy, tất kết mơ hình thủy động lực chuyển đổi sang định dạng đầu vào để sử dụng cho tính tốn mơ hình lan truyền TTLL Thời gian tính tốn - Mùa khô (từ ngày 01 đến ngày 31 tháng năm 2010) - Mùa mưa (từ ngày 01 đến ngày 31 tháng năm 2010) Điều kiện biên tham số tính Trên biển lỏng cửa sơng biển sử dụng số liệu TTLL trung bình mùa kết tính tốn từ mơ hình so sánh với số liệu quan trắc (Hình 4.1) So sánh giá trị TTLL kết tính từ mơ hình số liệu quan trắc có xu nhau, giá trị có khác biệt đáng kể Tuy nhiên số liệu quan trắc tiến hành đo tầng mặt (cách mặt 0,5 mét), cịn kết từ mơ hình tính cho trung bình cột nước Đây hạn chế q trình tính tốn nên nhiều ảnh hưởng đến kết nghiên cứu IV.3 Kết tính tốn IV.3.1 Dịng chảy  Mùa khô: vào mùa khô, biến đổi mùa trường gió suy giảm đáng kể lưu lượng nước từ sông đưa tạo khác biệt tương đối trường dòng chảy so với mùa mưa Ở thời kỳ cuối pha triều lên nước lớn, vận tốc dòng chảy nhỏ, ảnh hưởng khối nước sông vào thời điểm hạn chế nên khối nước biển xâm nhập sâu vào phía cửa sơng Thời gian chuyển pha nước lớn thời điểm triều xuống nhỏ, khoảng Trong thời điểm nước ròng, trường dòng chảy khu vực ven biển Hải Phịng có vận tốc nhỏ - khoảng 0,1 đến 0,2 m/s phân tán mạnh hướng chảy Cũng tải lượng nước sông nhỏ nên thời gian dừng chảy vào thời điểm nước ròng ngắn hơn, trường dòng chảy nhanh chóng chuyển trạng thái từ dừng chảy thành chảy lên Khu vực sơng Bạch Đằng có vận tốc 0,7-0,8 m/s giảm dần đến Đình Vũ với vận tốc dịng chảy 0,2-0,4 m/s Phía Đơng Bắc khu vực nghiên cứu lượng nước từ Quảng Ninh đưa sang lưu lượng từ sông Chanh chảy thấp nên tốc độ dịng chảy mùa khơ đạt từ 0,1-0,2 cm/s Kết mơ trường dịng chảy pha triều lên vào mùa khô cho thấy ảnh hưởng sâu vào lục địa khối nước biển, đặc biệt khu vực cửa Bạch Đằng Hướng dịng chảy vùng phía ngồi cửa sơng đồng định hướng phía sơng Tuy nhiên đến khu vực Đình Vũ với vận tốc 0,2-0,3 m/s, vào sâu ảnh hưởng dịng chảy sơng Bạch Đằng với vận tốc 0,5-0,7 m/s Với hình dạng đường bờ, mùa khơ xu hướng di chuyển khối nước khu vực nghiên cứu phía nam tây nam bán đảo Đồ Sơn nhiều Mặc dù tải lượng nước từ sông đưa mùa khô giảm mạnh tăng cường gió đơng bắc làm cho xu hướng mùa khô thể rõ rệt  Mùa mưa: vào mùa mưa lưu lượng nước sông lớn nên pha triều lên, vận tốc dịng chảy từ biển hướng vào phía cửa sơng có giá trị nhỏ Trong pha triều này, hướng dịng chảy chủ yếu nam - đơng nam với giá trị vận tốc biến đổi từ 0,2-0,7m/s Ở khu vực cửa sông Cấm-Bạch Đằng, nơi lưu lượng nước từ sông lớn sông đưa vùng ven biển Hải Phịng khơng có dịng chảy ngược từ biển vào Ở thời điểm nước lớn, hướng dòng chảy khu vực ven biển Hải Phòng phân tán mạnh mẽ với giá trị vận tốc nhỏ, đặc biệt vùng nước Hòn Dấu, Cát Bà Cát Hải Một số nơi lòng dẫn hẹp khu vực cửa Lạch Huyện, cửa Nam Triệu vận tốc dịng chảy đạt đến giá trị 1,0m/s Ở thời điểm nước ròng, khối nước từ sơng có điều kiện phát triển mạnh mẽ phía biển, nhiên giới hạn vốn có lưu lượng nước sơng nên dịng chảy có hướng phía biển tồn phạm vi khoảng 10-20 km từ bờ phía ngồi Khi trường dịng chảy khu vực sông Bạch Đằng lớn, biến thiên từ 0,9-1,4 m/s Các kết tính tốn cho thấy xu hướng dịch chuyển phía nam tây nam bán đảo Đồ Sơn khối nước trạng thái biến đổi khác pha triều (trừ pha triều lên) Nguyên nhân tượng dồn ép khối nước trao đổi nước nên phía bắc hạn chế địa hình phía đơng bắc khu vực nghiên cứu nơng IV.3.2 Trầm tích lơ lửng  Mùa khô: đặc điểm vận chuyển lan truyền trầm tích lơ lửng (TTLL) có liên quan chặt chẽ đến chế độ thủy động lực nguồn cung cấp trầm tích Vào mùa khơ có đặc điểm suy giảm mạnh mẽ nguồn cung cấp trầm tích từ cửa sơng khiến cho vào thời điểm nước lớn, dòng TTLL phát tán hạn chế vùng ven biển khu vực nghiên cứu so với mùa mưa Chỉ vùng nước nhỏ phía ngồi cửa Nam Triệu có hàm lượng TTLL tương đối cao (từ 20-40mg/l) cịn lại khu vực khác có hàm lượng TTLL nhỏ Vào pha triều xuống dịng trầm tích lơ lửng từ lục địa có điều kiện thuâ ̣n lới để phát tán phía ngồi Vùng nước có hàm lượng TTLL tương đối cao (khoảng 3050mg/l) tập trung sát cửa Nam Triệu, Lạch Tray Cơ chế lan truyền biến đổi TTLL vùng cửa sông ven biển khu vực nghiên cứu trường hợp triều lên mùa khô tương tự mùa mưa Mặc dù suy giảm nguồn trầm tích lơ lửng tải lượng nước t sông đưa làm cho ảnh hưởng khối nước biển có hàm lượng TTLL thấp trở lên mạnh mẽ, vùng nước có hàm lượng TTLL thấp tiến sâu vào lịng sơng khu vực nghiên cứu, sơng có tải lượng nước thấp Lạch Tray Ở khu vực phía tây nam đảo Cát Bà, hàm lượng TTLL khu vực nhỏ với giá trị dao động 10mg/l Điều giải thích sau, nơi chịu ảnh hưởng nguồn trầm tích từ lục địa đưa so với khu vực cịn lại phạm vi tính tốn hai mùa giá trị TTLL ổn định Ở khu vực phía nam tây nam đảo Cát Hải, vị trí gần cửa sơng phía bắc bán đảo Đồ Sơn nên hàm lượng TTLL nước giảm dần từ cửa sơng khu vực dao động khoảng 5-30mg/l Biến thiên theo thời gian hàm lượng TTLL phụ thuộc chặt chẽ vào dao động mực nước thể vai trị dịng trầm tích từ lục địa Giá trị hàm lượng TTLL thường đạt cực đại nước ròng giảm dần thủy triều tăng lên đến cực tiểu hàm lượng vào gần thời điểm nước lớn Biến thiên hàm lượng TTLL mạnh vào ngày triều cường Vào ngày tính tốn cuối khu vực phía nam đảo Cát Hải TTLL có giá trị khoảng 40-45 mg/l, khu vực vụng Cát Bà (các khu ni trồng thủy sản, bãi tắm) có giá trị 25-30 mg/l  Mùa mưa: vào mùa mưa hàm lượng TTLL sông khu vực nghiên cứu có giá trị lớn 100 mg/l Khu vực sơng Cấm có hàm lượng trầm tích cao sơng khác Với đặc điểm kết mô phân bố TTLL mùa mưa 2010 cho thấy vùng nước có hàm lượng TTLL cao chủ yếu xuất khu vực cửa Nam Triệu sông Cấm với ảnh hưởng từ nguồn trầm tích từ phía thượng nguồn Do ảnh hưởng trường dòng chảy nên phân bố biến động TTLL vùng cửa sông ven biển khu vực nghiên cứu chủ yếu theo pha dao động mực nước triều Sự phát tán tán TTLL từ sông vùng ven biển bị hạn chế tập trung sát cửa sông Trong khu vực cịn lại bị khối nước biển với hàm lượng TTLL nhỏ (< 30mg/l) Biến đổi theo thời gian hàm lượng TTLL: khu vực phía tây nam đảo Cát Bà, biến động hàm lượng TTLL theo thời gian cho thấy khu vực không chịu ảnh hưởng nguồn trầm tích từ lục địa Hàm lượng TTLL khu vực nhỏ với giá trị dao động nhỏ 20mg/l Ở khu vực phía nam tây nam đảo Cát Hải, vị trí gần cửa sơng phía bắc bán đảo Đồ Sơn nên hàm lượng TTLL nước giảm dần từ cửa sông khu vực cịn dao động khoảng 15 60mg/l Biến thiên theo thời gian hàm lượng TTLL phụ thuộc chặt chẽ vào dao động mực nước thể vai trị dịng trầm tích từ lục địa Giá trị hàm lượng TTLL thường đạt cực đại thủy triều xuống giảm dần thủy triều lên đến cực tiểu hàm lượng vào gần thời điểm nước lớn Biến thiên hàm lượng TTLL mạnh vào ngày triều cường Vào ngày tính tốn cuối cho thấy ven bờ phía tây khu vực nghiên cứu hàm lượng TTLL có giá trị khoảng 100 mg/l (từ cửa sông Cấm-Bạch Đằng đến bán đảo Đồ Sơn) Khu vực vực vụng Cát Bà (các khu nuôi trồng thủy sản, bãi tắm) có giá trị khoảng 40 mg/l Do nước sơng mùa mưa từ lục địa đưa lớn nhiều so với mùa khô nên phạm vi không gian lan truyền TTLL khối nước có giá trị > 100 mg/l vượt qua đảo Hòn Dấu KẾT LUẬN Luận văn thu thập số liệu trầm tích lơ lửng thay đổi theo thời gian (19962010), không gian (các cửa sơng, đảo) vùng cửa sơng ven biển Hải Phịng cách có hệ thống từ trước đến  Kết phân tích, đánh giá số liệu quan trắc cho thấy tranh tổng thể đặc điểm trầm tích lơ lửng có đặc trưng mùa rõ rệt phân bố chúng chịu ảnh hưởng mạnh q trình tương tác lục địa-biển  Mơ đun thủy động lực chất lượng nước mơ hình delft3d sử dụng để tính tốn trường thủy động lực, lan truyền trầm tích lơ lửng từ cửa sơng Hải Phịng biển Kết tính tốn cho thấy, xu lan truyền, hàm lượng TTLL có biến đổi theo mùa, theo pha triều rõ Vào mùa khơ, hướng lan truyền phía đơng nam (qua đảo Cát Bà) với hàm lượng đạt 42 mg/l (pha triều xuống), đơi phạm vi lan truyền TTLL có giá trị khoảng 40 mg/l đến tận mũi Đồ Sơn-sát với khu vực đảo Hịn Dấu Vào mùa mưa, lưu lượng sơng đưa lớn nhiều so với mùa khô, hướng lan truyền hướng đơng nam với hàm lượng đạt 120 mg/l (pha triều xuống) Phạm vi lan truyền TTLL có giá trị > 100 mg/l qua khu vực đảo Hòn Dấu, nhiên suốt trình tính tốn, khu vực Bến Gót khơng có tượng xảy đặc biệt khu vực phía tây đảo Cát Bà giá trị lớn đạt 30 mg/l  Trong thời gian tiếp theo, học viên tiến hành nghiên cứu phân tầng trường dịng chảy khu vực cửa sơng ven biển Hải Phịng tính tốn trầm tích lơ lửng có tính đến hoạt động kinh tế-xã hội (nguồn thải)  References Đỗ Trọng Bình, Trần Anh Tú, Vũ Duy Vĩnh (2010) Nghiên cứu đánh giá lan truyền chất gây ô nhiễm khu vực cửa sông ven biển Hải Phịng mơ hình tốn học Báo cáo tổng hợp Đề tài cấp Thành phố Hải Phòng Mã số: ĐT.MT.2008.500 Đỗ Đình Chiến, Vũ Duy Vĩnh, Trần Anh Tú (2005) Mơ q trình vận chuyển phân bố trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng mơ hình delft-3D Báo cáo tổng kết Nhiệm vụ sở năm 2005 Lưu trữ Viện TN&MTB Cục Môi trường, 1999-2008 Báo cáo quan trắc môi trường biển hàng năm Trạm Quan trắc phân tích mơi trường biển miền Bắc, từ năm 1999 đến 2008 Lưu trữ Viện tài nguyên Môi trường biển Nguyễn Đức Cự (2011) Nghiên cứu, đánh giá tác động cơng trình hồ chứa thượng nguồn đến diễn biến hình thái tài ngun - mơi trường vùng cửa sông ven biển đồng Bắc Bộ Báo cáo tổng hợp Đề tài độc lập cấp Nhà nước (Mã số: ĐTĐL 2009T/05) Nguyễn Đức Cự, Nguyễn Văn Thảo, Vũ Duy Vĩnh (2011) Nghiên cứu đánh giá tác động thủy thạch - động lực hệ thống đê quai lấn biển phục vụ xây dựng Sân bay quốc tế khu vực ven bờ Tiên Lãng - Hải Phòng Báo cáo tổng hợp Nhiệm vụ cấp thành phố Hải Phòng Lưu Văn Diệu, Nguyễn Chu Hồi, Nguyễn T P Hoa (2001) Đánh giá mức độ ô nhiễm nguồn thải từ lục địa, đề xuất giải pháp kiểm sốt, quản lý nhiễm nguồn lục địa đưa số khu vực cửa sông ven biển phía bắc (từ Quảng Ninh đến Thanh Hóa) Báo cáo tổng hợp Nhiệm vụ cấp Viện KH&CNVN Lưu Văn Diệu (2010) Đánh giá biến động thông số độ đục, nồng độ chất rắn lơ lửng (TTLL), nhu cầu oxy sinh hóa (BOD5) nhu cầu oxy hóa học (COD) nước biển ven bờ phía Bắc qua kết quan trắc trạm quan trắc phân tích mơi trường miền Bắc Việt Nam từ năm 1998-2008 Tuyển tập Tài nguyên môi trường biển, tập XV; Nhà xuất Khoa học Tự nhiên Công nghệ; tr 321-334 Nguyễn Minh Hải (2010) Nghiên cứu tượng nước dâng vùng ven biển Hải Phòng Báo cáo tập Hội đồng lịch sử thành phố Hải Phịng (1990) Địa chí Hải Phịng 10 Phạm Sỹ Hồn (2009) Nghiên cứu vận chuyển trầm tích từ cửa sơng biển vịnh Bình Cang-Nha Trang mơ hình tốn Luận văn Thạc sỹ 11 Trần Đình Lân, Lê Xuân Sinh (2008) Dự báo nguy ô nhiễm đề xuất giải pháp bảo vệ môi trường khu cơng nghiệp Bến Rừng, Thủy Ngun, Hải Phịng Đề tài cấp thành phố Hải Phòng (Mã số: ĐT.MT.2006.446) 12 Trần Đình Lân, Nguyễn Văn Thảo, Nguyễn T T Hà (2010) Đánh giá trạng môi trường xác định vấn đề ưu tiên phục vụ quản lý tổng hợp vùng bờ biển Hải Phòng Báo cáo Tổng hợp Đề tài cấp thành phố Hải Phòng Mã số: ĐT.MT.2008.498 13 Leo C VanRijn (1993) Principles of sediment transport in rivers, estuaries and coastal seas AQUA Publacations 14 Bùi Hồng Long (2004) Một số kết khảo sát, nghiên cứu tượng xói lở bồi tụ khu vực ven biển Bình Thuận Tuyển tập nghiên cứu biển, tập XIV, Nxb KHKT 15 Mulla D J and Addiscott T M (2005) Calibration and Validation of Water shedScale Models 16 Nguyễn Thọ Sáo, Nguyễn Minh Huấn, Ngơ Chí Tuấn, Đặng Đình Khá (2010) Biến động trầm tích diễn biến hình thái khu vực cửa sông ven bờ Cửa Tùng, Quảng Trị Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên Cơng nghệ 26, Số 3S (2010); tr 427-434 17 Takashi TAKANEZAWA (2000) NAOTIDE User manual 18 Trần Hồng Thái, Lê Vũ Việt Phong, Nguyễn Thanh Tùng, Phạm Văn Hải (2010) Mô phỏng, dự báo trình vận chuyển bùn cát lơ lửng khu vực Cửa Ông Tuyển tập báo cáo Hội thảo khoa học lần thứ 10, Viện Khoa học KTTV&MT; tr 332-341 19 Trần Đức Thạnh, Nguyễn Đức Cự, Nguyễn Hữu Cử, Đỗ Đình Chiến (2001) Nghiên cứu dự báo, phòng chống sạt lở biển Bắc Bộ từ Quảng Ninh tới Thanh Hoá Báo cáo tổng hợp Đề tài cấp Nhà nước (Mã số: KHCN.5A) 20 Trần Đức Thạnh, Đỗ Đình Chiến, Trần Anh Tú, Nguyễn T Kim Anh (2007) Xây dựng mơ hình lan truyền chất nhiễm cho Vịnh Hạ Long - Vịnh Bái Tử Long Báo cáo tổng hợp Dự án cấp tỉnh Quảng Ninh 21 Trần Đức Thạnh, Vũ Duy Vĩnh, Yoshiki Saito, Đỗ Đình Chiến, Trần Anh Tú (2008) Bước đầu đánh giá ảnh hưởng đập Hịa Bình đến mơi trường trầm tích ven bờ châu thổ sơng Hồng Tạp chí Khoa học Công nghệ Biển (Phụ trương (T.8)/2008); tr 01-17 22 Trần Đức Thạnh, Nguyễn T P Hoa, Cao T T Trang (2008) Đánh giá tình trạng nhiễm suy thối mơi trường khu vực cửa sơng Cấm - Bạch Đằng đề xuất giải pháp bảo vệ Báo cáo tổng hợp Nhiệm vụ cấp Viện KH&CNVN 23 Cao T T Trang, Vũ Thị Lựu (2007) Đánh giá khả tích tụ phân tán chất nhiễm vùng cửa sông ven biển Việt Nam Báo cáo tổng hợp Đề tài cấp Viện KH&CNVN 24 Cao T T Trang, Nguyễn T P Hoa, Dương Thanh Nghị (2009) Đánh giá sức tải môi trường đảo Cát Bà đề xuất giải pháp phát triển bền vững Báo cáo tổng hợp Đề tài cấp Thành phố Hải Phòng (Mã số: ĐT.MT.2006.442) 25 Cao T T Trang, Nguyễn Mạnh Thắng, Lê Xuân Sinh (2009) Thử nghiệm đánh giá sức tải môi trường sông Bạch Đằng đề xuất giải pháp bảo vệ môi trường, phát triển bền vững Báo cáo tổng hợp Đề tài cấp Viện KH&CNVN 26 Trần Anh Tú, Vũ Duy Vĩnh (2011) Mô lan truyền chất gây ô nhiễm khu vực cửa sơng ven biển Hải Phịng Hội nghị Khoa học cơng nghệ biển tồn quốc lần thứ V, Hà Nội Tập Quyển 2, tr 163-170 27 U.S Army Corps of Engineers (2002) Coastal Engineering Manual EM 1110-21100, Part III 28 Đinh Văn Ưu (2009) Mơ hình vận chuyển trầm tích biến động địa hình đáy áp dụng cho vùng biển cửa sơng cảng Hải Phịng Tạp chí Khoa học ĐH QG Hà Nội, Khoa học Tự nhiên Công nghệ 25, số 1S (2009); tr 133-139 29 Bùi Văn Vượng, Trần Đức Thạnh, Đỗ T T Hương, Cao T T Trang (2007) Nạo vét cảng Hải Phòng số ảnh hưởng đến mơi trường hệ sinh thái biển Tuyển tập báo khoa học Hội nghị KHCN&MT (Kỷ niệm 10 năm thành lập Trung tâm đào tạo TT KHCN Bảo vệ Môi trường thủy, Trường ĐH Hàng hải Việt Nam); tr 202-208 30 WL/Delft Hydraulics (1999) Delft3D-FLOW User Manual Version 3.05, Delft3DWaq User Manual Version 3.01 Delft, Netherlands  ... nhiên khu vực Phòng Hải Chƣơng 2: Tài liệu phương pháp Chƣơng 3: Đánh giá trạng trầm tích lơ lửng khu vực cửa sơng ven biển Hải Phịng Chƣơng 4: Mơ trầm tích lơ lửng khu vực cửa sơng ven biển Phịng... vùng cửa sơng ven bờ Hải Phịng nói chung [5] Bởi vậy, việc đánh giá TTLL vùng cửa sông ven biển Hải Phòng điều cần thiết Với mục tiêu luận văn mô trạng trầm tích lơ lửng theo mùa khu vực cửa sơng... suất nhỏ 0,15% [1, 2, 12] Trầm tích lơ lửng khu vực cửa sơng ven biển khu vực nghiên cứu nhiều nguồn cung cấp khác nguồn chủ yếu từ sơng đưa Hàm lượng trầm tích lơ lửng sông biến thiên khoảng rộng,