Nghiên cứu trình bày chế độ thủy động lực đầm Đề Gi chủ yếu là sự kết hợp giữa dòng chảy do thủy triều, dòng chảy do sông đổ ra và chịu ảnh hưởng của yếu tố địa hình. Mô hình 2D được xây dựng dựa vào hệ phương trình nước nông phi tuyến đã mô phỏng được bức tranh dòng chảy trong đầm.
Tuyển Tập Nghiên Cứu Biển, 2015, tập 21, số 1: 11-20 ĐẶC ĐIỂM DỊNG CHẢY ĐẦM ĐỀ GI, BÌNH ĐỊNH Phạm Sỹ Hoàn Viện Hải dương học, Viện Hàn lâm Khoa học & Cơng nghệ Việt Nam Tóm tắt Chế độ thủy động lực đầm Đề Gi chủ yếu kết hợp dòng chảy thủy triều, dòng chảy sông đổ chịu ảnh hưởng yếu tố địa hình Mơ hình 2D xây dựng dựa vào hệ phương trình nước nơng phi tuyến mơ tranh dịng chảy đầm Dịng chảy đầm có hai hướng ngược chiều triều lên triều xuống, mạnh mùa mưa, đạt 100 cm/s khu vực trao đổi nước sông- đầm đầmbiển, khu vực khác chủ yếu nhỏ 40 cm/s Dịng chảy mùa khơ chủ yếu nhỏ 20 cm/s, khu vực trao đổi nước đầm- biển đầm- sơng đạt 60 cm/s Có thể tồn xốy nước cục khu vực cửa sơng- đầm khu vực phía cửa đầm- biển tương tác thủy triều, nước sơng yếu tố địa hình SIMULATION OF THE CURRENT SYSTEM IN DE GI LAGOON (BINH DINH PROVINCE) Pham Sy Hoan Institute of Oceanography, Vietnam Academy of Science & Technology Abstract The hydrodynamical regime in the De Gi lagoon is mainly as a result of interaction between tides, rivers and terrain The 2D model build based on nonlinear shallow water equations was successfully simulated the current field in the De Gi lagoon The current in lagoon is in the opposite directions in the ebb and flood tides The current in the rainy season is stronger than that in dry season In the rainy season, the current speed is mostly less than 40 cm/s It can reach over 100 cm/s in the areas of water exchange between lagoon and sea, lagoon and river In the dry season, the current speed is mostly less than 20 cm/s and may over 60 cm/s in the areas of water exchange between lagoon and sea, lagoon and river The interaction of the tides, rivers and terrain can form the local eddies in the river mouth- lagoon and lagoon- sea areas I MỞ ĐẦU bắt thủy sản… Đặc biệt nguồn sa khoáng titan dồi dào, mỏ sa khống Đề Gi mỏ lớn tỉnh Bình Định lớn thứ hai nước Đầm Đề Gi lưu vực tương đối đơn giản theo quan điểm điều kiện địa lý tự nhiên Địa hình đầm tương đối nơng, chỗ sâu cửa đầm khoảng m Đầm có cửa thơng với biển, có hai nhánh sông thuộc loại nhỏ đổ Vào mùa khô, hai nhánh sơng khơng Đầm Đề Gi nằm phía bắc cách Quy Nhơn (Bình Định) khoảng 40 km, có tổng diện tích 2.000 ha, thuộc địa bàn xã Cát Khánh, Cát Minh (Phù Cát) Mỹ Cát, Mỹ Chánh, Mỹ Thành (Phù Mỹ) Theo đánh giá nhà khoa học, đầm Đề Gi có tiềm đa dạng tài nguyên sinh vật, phi sinh vật, nuôi trồng đánh 11 thể rõ vai trò cấp nước vào đầm Vào mùa mưa, hai nhánh sông không gây lũ lớn, khốc liệt sơng khác (lưu lượng nước trung bình ngày mùa lũ trạm An Xuyên khoảng 335 m3/s, theo trung tâm dự báo KTTV Bình Định) Một yếu tố thủy động lực quan trọng góp phần trực tiếp hay gián tiếp tác động vào chất lượng mơi trường, sinh thái đầm Đề Gi nói riêng, vùng biển khác nói chung dịng chảy Để nắm rõ đặc trưng dịng chảy, ngồi việc phải tiến hành điều tra, khảo sát, mơ hình hóa q trình cơng cụ đắc lực để nghiên cứu vấn đề Hiện nay, có nhiều mơ hình tốn tính tốn thủy- thạch động lực khu vực vũng vịnh ven bờ, có dịng chảy Các mơ hình tốn có uy tín nâng cấp liên tục giới thường mô hình thương mại, điển MIKE (DHI Group, 2011), EFDC (USA Environmental Protection Agency), TELEMAC (EDF-DER, Paris), HYDRO2D, ADCIRC (CHL, USA),… Ứng dụng thành cơng mơ hình MIKE, Phạm Thu Hương cs (2011) tính tốn dịng chảy sóng khu vực cửa sơng Đà Rằng; Nguyễn Chí Cơng cs (2012) tính dịng chảy lan truyền chất gây ô nhiễm lên bãi tắm ven bờ tây vịnh Nha Trang; Phạm Sỹ Hoàn cs (2013) tính dịng chảy tổng hợp vịnh Quy Nhơn… Ngồi mơ hình thương mại, số mơ hình mã nguồn mở tính tốn dịng chảy khu vực ven bờ, vũng vịnh nói riêng, yếu tố thủy- thạch động lực khác nói chung ECOMSED (HydroQual, Inc, 2004), HAMSOM, ECOSMO, ROM… Bùi Hồng Long Phạm Xn Dương (2010) sử dụng mơ hình ROM để tính tốn dịng chảy cho vịnh Bình Cang- Nha Trang Phạm Sỹ Hồn Lê Đình Mầu (2011) khai thác ứng dụng mơ hình ECOMSED tính tốn dịng chảy lan truyền vật chất ô nhiễm khu vực ven bờ cửa sông Mê Kông Cũng ứng dụng ECOMSED, tác giả tập thể nghiên cứu tính tốn dòng chảy, vận chuyển vật chất biến đổi địa hình đáy cho khu vực cửa sơng Tắc, Nha Trang Trong nghiên cứu này, tác giả sử dụng cải tiến mơ hình xây dựng mơ thành cơng dịng chảy lan truyền vật chất cho vịnh Bình Cang - Nha Trang vào mùa mưa để mơ dịng chảy đầm Đề Gi tỉnh Bình Định Mơ hình kiểm chứng với số liệu đo mực nước, số liệu đo dòng chảy 02 mùa (mùa mưa: 10/2009 mùa khô: 04/2010) Kết kiểm chứng cho thấy mơ hình mơ tốt dịng chảy đầm Đề Gi II TÀI LIỆU VÀ PHƯƠNG PHÁP Tài liệu Số liệu đo dòng chảy tầng (mặt, đáy) 16 trạm mặt rộng 01 trạm liên tục 01 ngày đêm hai mùa (mùa mưa: tháng 10/2009, mùa khô: tháng 4/2010) Trạm đo liên tục mặt cắt cửa Đề Gi- biển thực 03 thủy trực (Hình 1); Số liệu đo mực nước thời gian với đo dòng chảy trạm liên tục cảng cá Đề Gi; Địa hình khu vực nghiên cứu; Các số liệu thực dự án “Điều tra, đánh giá trạng đề xuất giải pháp quản lý tổng hợp đầm Đề Gi theo hướng phát triển bền vững” Viện Hải dương học thực hiện; Ngoài ra, số liệu đo lưu lượng nước sông La Tinh trạm An Xuyên (từ 25/10 đến 08/11/2007) trạm An Mỹ (từ 25/10 đến 08/11/2007) Trung tâm dự báo KTTV Bình Định thực (2009) Vị trí hai trạm thể hình Phương pháp Đã ứng dụng mơ hình dịng chảy 2D tác giả xây dựng mơ thành cơng dịng chảy vận chuyển vật chất lơ lửng cho vực nước Bình Cang - Nha Trang vào mùa mưa (Phạm Sỹ Hoàn Nguyễn Thọ Sáo, 2009; Báo cáo tổng kết đề tài cấp Cơ sở, Viện Hải dương học Phạm Sỹ Hồn cs thực năm 2010) để mơ dòng chảy đầm Đề Gi ảnh hưởng thủy triều nước sơng 12 14.18 Xã Mỹ Thành BIỂN ĐÔNG 21 14.17 N Hồng Sơn Xã Mỹ Chánh 14.16 Bờ bắc M1 M3 Bờ nam M2 Hưng Lạc ĐẦM ĐỀ GI 10 Xã Mỹ Cát 14.15 Vónh Lợi 14.14 Hòn Giữa 12 Xã Cát Minh 13 14.13 14 11 15 T1 T2 16 Xã Cát Khánh 14.12 109.15 109.16 109.17 109.18 109.19 Cảng Cá An Quang 109.2 109.21 Hình Sơ đồ địa hình khu vực nghiên cứu vị trí trạm đo Fig Bathymetry of the study area and positions of surveyed stations Mơ hình dùng hệ phương trình nước nơng phi tuyến mơ trường dịng chảy trung bình theo độ sâu, có tính đến ma sát đáy, nhớt rối ngang, bỏ qua ma sát mặt vùng nghiên cứu có địa hình núi bao quanh, bề mặt nước có diện tích khơng lớn (khoảng 2.000 ha) nên tác động gió bỏ qua Hệ phương trình bao gồm: Phương trình bảo tồn khối lượng (1): ∂ζ ∂ (h + ζ )u ∂ (h + ζ )v + + =0 ∂t ∂x ∂y (1) Phương trình bảo tồn động lượng (2) (3): τ xb ∂u ∂u ∂u ∂ 2u ∂ 2u ∂ζ + u + v − ε − ε − fv = − g + ∂t ∂x ∂y ∂x ∂y ∂x ρ (h + ζ ) (2) b τy ∂v ∂v ∂v ∂ v ∂ v ∂ζ + u + v − ε − ε + fu = − g + ∂t ∂x ∂y ∂x ∂y ∂y ρ (h + ζ ) 2 (3) τ x,y = ρ.Cb.(u + v ) (u,v) (4) Trong đó: u, v – thành phần tốc độ dịng chảy trung bình toàn cột nước theo phương Ox Oy; h – độ sâu cột nước b 2 1/2 13 trạng thái mặt biển yên tĩnh; ζ – dao động mặt biển so với mặt biển yên tĩnh; f – tham số Coriolis; g- gia tốc trọng trường; τbx,y - ứng suất ma sát đáy tương ứng với phương Ox, Oy; ρ – mật độ nước biển; Cb- hệ số ma sát đáy Hệ số ma sát đáy: Theo thực nghiệm thơng thường mơ hình tốn chiều, hệ số nhận giá trị từ 0,0026 – 0,0030; Để tiếp cận với điều kiện địa phương thực tế, nghiên cứu sử dụng hệ số ma sát đáy phụ thuộc vào đặc tính đáy thể qua cấp hạt đáy độ sâu thể qua hệ số ma sát đáy Chezy theo công thức: Cb = g/C2; với C = 18 log(12h /r) – hệ số Chezy; r= 2,5 d50/30 – độ ghồ ghề (độ nhám) đáy biển; d50 – đường kính hạt vật chất đáy Đường kính hạt d50 tính tốn từ kết phân tích học trầm tích khu vực đầm Đề Gi lưu trữ phịng Địa chất biển (Viện Hải dương học) cơng bố Phạm Bá Trung (2012) Các phương trình rời rạc theo sơ đồ Arakawa- C (Mesinger Arakawa, 1976), đó, biến ζ đặt tâm ô lưới, biến u, v đặt cạnh ô lưới Các phương trình giải theo phương pháp sai phân hữu hạn Điều kiện ổn định số lấy theo Courant- Friedrichs- Levy Bước thời gian dt = 1(s); bước không gian dx = dy = 1(s) Sơ đồ sử dụng nhiều toán thủy động lực đơn giản tính tốn nhanh, điển Kiyoshi Horikawa (1987), Christopher G Koutitas (1988) Tại biên lỏng phía sơng, cho lưu lượng nước trạm An Xun, An Mỹ (Hình 2) Do có lưu lượng vào mùa mưa, nên lưu lượng nước sông vào mùa khơ xác định theo đặc điểm dịng chảy mùa khô chiếm khoảng 30% năm (Đặc điểm Khí hậu Thủy văn Bình Định, 2004) Lưu lượng trung bình ngày thời kỳ mưa lũ điển hình trạm An Mỹ khoảng 42,5 m3/s; trạm An Xuyên khoảng 162,6 m3/s Hình Lưu lượng trung bình ngày (m3/s) sông La Tinh trạm An Xuyên An Mỹ từ ngày 25/10 đến 01/11/2007 Fig The daily average water discharge of La Tinh river at An Xuyen and An My stations during 25/10 – 01/11/2007 Tại biên lỏng phía biển cho dao động mực nước dựa vào số liệu mực nước trạm Quy Nhơn cách vùng nghiên cứu phía nam khoảng 40 km tính theo đường thẳng Kết tính tốn so sánh với số liệu dao động mực nước (Hình 3) dịng chảy thực đo (Hình 4, hình 5) cho thấy kết mơ hình tính tốn hợp lý Cụ thể: -Về dao động mực nước, sai số tuyệt đối dao động mực nước vào thời kỳ mùa khơ lớn 14,16 cm, trung bình 5,91 cm; sai số tuyệt đối lớn vào mùa mưa 14,56 cm, trung bình 5,51 cm Các sai số lớn chủ yếu tập trung vào lúc chuyển pha triều Điều quan trọng pha dao động mực nước kết tính tốn mơ hình số liệu thực đo 14 -Về dòng chảy, tính tốn tài liệu thực đo có dao động gần pha (Hình 4, hình 5) Các thành phần dịng chảy thực đo có phần lớn dịng chảy tính tốn từ mơ hình Sai số tuyệt đối mùa khô (04/2010) lớn khoảng 16,3 cm/s, trung bình khoảng cm/s; mùa mưa (10/2009) tương ứng 20 cm/s 8,8 cm/s Các sai số lớn gặp phải thường xuất chuyển pha triều Trong mơ hình tính tốn, u cầu lưu lượng nước sông làm đầu vào liên tục biến đổi theo bước thời gian tính (1 s), tài liệu lưu lượng thu dạng trung bình ngày Đây có lẽ nguyên nhân gây sai số lớn dao động mực nước dòng chảy chuyển pha triều (a) (b) Hình Dao động mực nước theo tính tốn mơ hình theo số liệu đo (a- mùa khô: 04/2010; b- mùa mưa: 10/2009) Fig Comparison between model and measured data of sea water level (a- dry season: 04/2010; b- rainy season: 10/2009) (a) (b) Hình Các thành phần tốc độ dòng chảy (a- thành phần u; b- thành phần v) theo tính tốn mơ hình số liệu thực đo vào mùa khô 4/2010 Fig Comparison between calculated model and measured data of velocity components (a- for u component; b- for v component) in dry season (4/2010) (a) (b) Hình Các thành phần tốc độ dịng chảy (a- thành phần u; b- thành phần v) theo tính tốn mơ hình số liệu thực đo vào mùa mưa 10/2009 Fig Comparison between calculated model and measured data of velocity components (a- for u component; b- for v component) in rainy season (10/2009) 15 III KẾT QUẢ Kết tính tốn lấy đại diện cho 02 mùa: mùa mưa (10/2009) mùa khô (4/2010) Mỗi mùa tương ứng, tính tốn cho pha triều đặc trưng (triều lên triều xuống) Mùa mưa Các kết thống kê từ tài liệu thực đo kết tính tốn cho bảng Tài liệu thực đo thống kê tức thời cho tồn vùng, cịn kết tính tốn thống kê cho chu kỳ dài tính tốn (5 ngày) Có thể thấy, dịng chảy có tốc độ lớn tập trung khu vực có địa hình biến đổi nhanh, diện tích mặt cắt ướt thu hẹp nhanh (vùng cửa sông vùng cửa đầm) Pha triều xuống, tốc độ dòng chảy lớn pha triều lên đáng kể Tốc độ dòng chảy trung bình thực đo pha triều lên khoảng 10,32 cm/s, pha triều xuống từ 20,53 – 22,37 cm/s Theo kết tính tốn, pha triều lên, tốc độ dịng chảy trung bình dao động từ 6,03 – 35,76 cm/s, pha triều xuống từ 6,45 – 69,89 cm/s Bảng Một số đặc trưng dòng chảy đầm Đề Gi vào mùa mưa (10/2009) Table Some features of the current in De Gi lagoon in rainy season (10/2009) Thời gian Pha triều lên Pha triều xuống Đặc trưng Cực đại Cực tiểu Trung bình Độ lệch chuẩn Cực đại Cực tiểu Trung bình Độ lệch chuẩn Tốc độ dịng chảy thực đo toàn khu vực (cm/s) Tầng mặt Tầng đáy 15,90 4,74 10,32 7,89 82,00 42,40 6,72 8,50 20,53 22,37 19,15 13,86 Tốc độ dịng chảy theo tính tốn mơ hình (cm/s) Cửa sông Cửa biển Trong đầm 65,45 17,61 8,24 11,45 3,26 5,82 35,76 8,06 6,03 14,15 3,46 1,59 77,71 26,19 45,43 42,93 8,95 1,79 69,89 20,52 6,45 5,96 4,58 6,17 Hướng dòng chảy mùa mưa đầm Đề Gi biến đổi đơn giản toàn vùng nghiên cứu theo biến đổi thủy triều (dòng chảy có hướng chảy thủy triều xuống ngược lại) Hướng dòng biến đổi phức tạp chuyển tiếp pha triều chảy vùng nghiên cứu có hướng chảy ra, tốc độ dịng chảy tăng dần lên theo thủy triều xuống (Hình 6b) Tốc độ dòng chảy đầm thường nhỏ 30 cm/s, tốc độ dịng cửa sơng cửa đầm đạt 100 cm/s 1.1 Pha triều xuống Khi triều bắt đầu lên, dịng chảy phía ngồi cửa đầm có hướng chảy vào, trễ pha, dịng chảy phía cửa sơng có hướng chảy Sự nghịch hướng gây nên xoáy nghịch cục chiều kim đồng hồ phần đỉnh đầm – cửa sơng Xốy tồn thời gian ngắn thủy triều lên đủ mạnh để đẩy khối nước sông trở vào Khi thủy triều lên mạnh, hướng dịng chảy tồn khu vực có hướng chảy vào Khu vực Khi triều bắt đầu xuống, dịng chảy phía đỉnh đầm có hướng chảy vào, dịng chảy khu vực cửa đầm lại có hướng chảy Trong thời điểm này, dịng chảy cửa sơng đỉnh đầm có tốc độ lớn vùng cịn lại Điều thể trễ pha dòng chảy từ khu vực cửa sông đỉnh đầm Theo tính tốn mơ hình, thời gian trễ pha khoảng 20 phút Khi triều xuống mạnh, hướng dòng 16 1.2 Pha triều lên phía đỉnh đầm, dịng chảy phía bờ tây có phần mạnh phần đỉnh phía bờ đơng Khu vực phía cửa đầm, tốc độ dịng chảy phía bờ đơng lại có phần vượt trội phía bờ tây (Hình 6a) Những đặc điểm có lẽ đặc điểm hình dạng đầm, với 02 cửa hẹp trao đổi nước 02 đầu (đỉnh 109.17 109.18 109.19 109.20 109.21 109.16 109.19 N Hồng Sơn BIỂN ĐÔNG 14.16 14.16 14.16 14.16 14.15 14.14 Hòn Giữa Hòn Giữa 14.13 14.13 14.12 14.12 xã Cát Minh 14.13 14.13 xã Cát Minh Vónh Lợi 14.14 14.14 xã Mỹ Caùt An Quang An Quang 109.16 109.17 109.18 109.19 109.20 109.21 109.15 (a) xã Cát Khánh : 100 cm/s 109.16 14.12 xã Cát Khánh : 50 cm/s 14.15 Vónh Lợi 14.14 xã Mỹ Cát Hưng Lạc 14.15 14.15 BIỂN ĐÔNG xã Mỹ Chánh Hưng Lạc 14.12 109.21 14.17 14.17 N Hồng Sơn xã Mỹ Chánh 109.15 109.20 xã Mỹ Thành 14.17 xã Mỹ Thành 14.17 109.18 14.18 14.18 14.18 109.17 14.18 109.16 cửa) Phần phía cửa biển, diện tích mở rộng phía tây nam, phần phía cửa sơng, diện tích đầm mở rộng phía đơng bắc Thêm vào đó, nước chảy vào, theo qn tính nó, với địa nêu hình thành nên đặc điểm dịng chảy lệch phía bờ nêu 109.17 109.18 109.19 109.20 109.21 (b) Hình Dịng chảy đầm Đề Gi tính tốn mơ hình mùa mưa 10/2009 (a- triều lên; b- triều xuống) Fig Simulated current field in De Gi lagoon in rainy season 10/2009 (a- flood tides; b- ebb tides) Mùa khô Cấu trúc hệ dịng chảy đầm Đề Gi vào mùa khơ tương tự mùa mưa với đổi đồng hướng dịng theo khơng gian, đổi hướng dịng theo lên xuống pha triều (dịng chảy có hướng chảy vào triều lên, hướng chảy triều xuống) Tốc độ dòng chảy lớn tập trung khu vực cửa sơng (phía đỉnh đầm) khu vực cửa đầm (bắt đầu đổ biển) Do đặc điểm hình dạng qn tính dịng chảy mà dịng chảy có đặc điểm lệch phía bờ tây khu vực đỉnh đầm lệch phía bờ đơng khu vực phía cửa đầm (đặc biệt pha triều lên) Tuy nhiên, khơng có đóng góp lượng nước sơng nên tốc độ dịng chảy lớn tốc độ dịng chảy vào khơng đáng kể (khoảng 10 cm/s khu vực cửa sông, cm/s khu vực cửa biển (cửa đầm), gần cm/s khu vực đầm) (Bảng 2) Trên bảng thấy được, tốc độ dòng chảy trung bình thực đo tính tốn mơ hình khơng lớn Xét cho tồn vùng nghiên cứu, tốc độ dịng chảy trung bình theo thực đo pha triều lên dao động từ 14,01 – 16,99 cm/s, pha triều xuống 13,53 – 13,68 cm/s Còn theo kết tính tốn mơ hình, pha triều lên 5,93 – 21,26 cm/s, pha triều xuống 8,65 – 32,6 cm/s 17 Bảng Một số đặc trưng dịng chảy đầm Đề Gi vào mùa khơ (04/2010) Table Some features of the current in De Gi lagoon in dry season (04/2010) Thời gian Pha triều lên Pha triều xuống Cực đại Cực tiểu Trung bình Độ lệch chuẩn Cực đại Cực tiểu Trung bình Độ lệch chuẩn 2.1 Pha triều xuống Khi thủy triều bắt đầu xuống, lượng nước sông đổ mùa khơng đáng kể nên dịng chảy tồn khu vực có hướng chảy ra, ngoại trừ phận dịng chảy chảy vào cửa đầm, lệch phía bờ bắc, kết hợp với dòng chảy cửa, gần bờ phía nam dịng chảy từ đầm tạo nên xốy thuận tức thời phía khu vực cửa đầm (Hình 8a) Xốy thuận yếu, quy mô nhỏ tồn thủy triều xuống mạnh Khi thủy triều xuống mạnh, toàn khu vực đầm có dịng chảy đồng theo khơng gian có tốc độ nhỏ đáng kể so với mùa mưa 109.16 109.17 109.18 Tốc độ dịng chảy theo tính tốn mơ hình (cm/s) Cửa sông Cửa biển Đầm 18,88 49,52 9,64 5,51 4,41 3,91 14,49 21,26 5,93 3,63 11,77 1,63 28,46 59,35 13,47 14,53 10,83 5,45 21,57 32,60 8,65 3,78 49,52 2,21 Tốc độ dòng chảy thực đo (cm/s) Tầng mặt Tầng đáy 22,18 20,68 9,13 5,67 16,99 14,01 6,05 5,84 38,62 19,43 8,37 9,87 13,53 13,68 4,60 3,66 Đặc trưng 109.19 109.20 Tốc độ dòng lớn chỗ trao đổi nước với biển, đạt gần 60 cm/s Tốc độ dòng vị trí khác chủ yếu nhỏ 30 cm/s (Hình 7b) 2.2 Pha triều lên Khi triều bắt đầu lên, đặc điểm địa hình mà xốy nghịch cục phía đỉnh đầm- cửa sơng hình thành mùa mưa (Hình 8b) Xốy biến triều xuống mạnh (Hình 7a) Tốc độ dịng chảy vào pha triều lên đạt khoảng 50 cm/s (vị trí trao đổi nước với biển) Tốc độ dịng vị trí khác chủ yếu nhỏ 20 cm/s 109.21 109.16 109.19 14.18 N Hoàng Sơn BIỂN ĐÔNG 14.16 14.16 14.16 14.16 14.15 14.14 Hòn Giữa Hòn Giữa 14.13 14.13 14.12 14.12 xã Cát Minh 14.13 14.13 xã Cát Minh Vónh Lợi 14.14 14.14 xã Mỹ Cát An Quang An Quang 109.17 109.18 109.19 109.20 109.21 109.15 xã Cát Khánh : 30 cm/s 109.16 14.12 xã Cát Khánh : 30 cm/s 14.15 Vónh Lợi 14.14 14.15 Hưng Lạc 14.15 xã Mỹ Cát 14.12 BIỂN ĐÔNG xã Mỹ Chánh Hưng Lạc 109.16 109.21 14.17 14.17 N Hồng Sơn xã Mỹ Chánh 109.15 109.20 xã Mỹ Thành 14.17 xã Mỹ Thành 14.17 109.18 14.18 14.18 14.18 109.17 109.17 109.18 109.19 109.20 109.21 (a) (b) Hình Dịng chảy đầm Đề Gi tính tốn mơ hình mùa khô 04/2010 (a- pha triều lên; b- pha triều xuống) Fig Simulated current field in De Gi lagoon in dry season 04/2010 (a- flood tides; b- ebb tides) 18 (a) (b) Hình Dịng chảy đầm Đề Gi lúc đỉnh triều (a- hình thành xốy thuận cục phía cửa đầm) lúc chân triều (b- hình thành xốy nghịch cục phía đỉnh đầm ngồi cửa sông) Fig Simulated current field in De Gi lagoon at high tide (a- existance of local cyclone in lagoon mouth) and low tide (b- existance of local anticyclonic in river mouth area) Nhìn chung, hệ dịng chảy đầm Đề Gi đơn giản có đặc điểm biến đổi theo mùa theo pha triều rõ rệt Kết mơ cho thấy tồn xốy nghịch cục khu vực cửa sơng phía đỉnh đầm đổi pha triều từ triều xuống sang triều lên (chân triều) Xoáy nghịch tác nhân làm giảm khả truyền tải – khuếch tán vật chất từ cửa sơng đầm Cịn đổi pha triều từ triều lên sang triều xuống (đỉnh triều), hình thành xốy thuận cục khu vực phía cửa trao đổi nước với biển làm cho vật chất khuếch tán từ đầm biển bị tồn đọng lại đầm với thời gian lâu Đến dịng chảy mạnh lên, tồn khu vực đầm có dịng chảy với tốc độ gần nhau, tức trình truyền tải – khuếch tán chất đầm biển diễn mạnh mẽ tính tốn mơ hình kiểm chứng với số liệu thực tế cho kết phù hợp Hệ dòng chảy đầm đơn giản, với 02 hướng chủ đạo chảy phía biển pha triều xuống chảy lên phía đỉnh đầm vào cửa sơng pha triều lên Tốc độ dịng chảy lớn mùa mưa có đóng góp đáng kể lưu lượng nước từ sơng đổ Trong mùa này, dịng chảy đạt tới 100 cm/s khu vực tiếp giáp với biển (pha triều xuống) đạt tới 70 cm/s khu vực cửa sơng Trong mùa khơ, dịng chảy khu vực cửa biển (chảy ra) đạt 60 cm/s, chảy vào đạt tới 50 cm/s, khu vực cịn lại, dịng chảy có tốc độ nhỏ 20 cm/s Dịng chảy đầm hình thành nên xốy cục cửa sơng phía đỉnh đẩm phía cửa đầm ảnh hưởng điều kiện địa hình, tương tác nước từ sông đổ đảo pha triều Với đặc thù đầm bán kín có trao đổi chất với lục địa với biển, thời gian dòng triều lên lớn thời gian dịng triều rút, tồn xốy cục đầm, với hoạt động kinh tế biển IV KẾT LUẬN VÀ THẢO LUẬN Đối với đầm có độ sâu nơng biến đổi đầm Đề Gi, mơ hình mơ dịng chảy 2D dùng để mơ q trình dịng chảy hồn tồn hợp lý Điều chắn kết 19 khu vực, vật chất có nguồn gốc sơng tải bị ứ đọng lại khu vực đầm phần, gây nên tượng bồi lấp khu vực cửa biển cửa sơng đầm Cần có điều tra, tính tốn cụ thể để khẳng định thêm luận Lời cảm ơn Tác giả xin trân trọng cảm ơn Ban chủ nhiệm dự án “Điều tra, đánh giá trạng đề xuất giải pháp quản lý tổng hợp đầm Đề Gi theo hướng phát triển bền vững” cho phép tác giả sử dụng tài liệu số liệu để hoàn thành nghiên cứu TÀI LIỆU THAM KHẢO Bùi Hồng Long Phạm Xuân Dương, 2010 Một số kết tính tốn dịng chảy theo mùa vùng vịnh Bình Cang – Nha Trang mơ hình ROMS Tuyển tập nghiên cứu biển, tập XVII: 30-42, Nhà xuất Khoa học Tự nhiên Công nghệ Christopher G Koutitas, 1988 Mathematical models in coastal engineering Pentech Press, London, 155 p Kiyoshi Horikawa, 1987 Nearshore dynamics and coastal processes University of Tokyo Press, 489 p Mesinger F and A Arakawa, 1976 Numerical methods used in atmospheric models Volume 1, Garp Publications series No 17, 64 p Nguyễn Chí Cơng, Nguyễn Minh Huấn Phan Thành Bắc, 2012 Mô lan truyền vật chất gây ô nhiễm từ cửa sông tới khu vực bãi tắm vịnh Nha Trang Kỷ yếu Hội nghị Quốc tế “Biển Đông 2012”, tập 2: 37-44, Nhà xuất Khoa học Tự nhiên Công nghệ Phạm Bá Trung, 2012 Đặc điểm địa hình trầm tích tầng mặt đầm Đề Gi Kỷ yếu Hội nghị Quốc tế “Biển Đông 2012”, tập 2: 232-242, Nhà xuất Khoa học Tự nhiên Cơng nghệ Phạm Sỹ Hồn Lê Đình Mầu, 2011 Tính tốn vận chuyển vật chất lơ lửng dải ven biển cửa sơng Mê Kơng mơ hình tốn Kỷ yếu Hội nghị KH&CN biển tồn quốc lần thứ V Quyển 2- Khí tượng, thủy văn động lực học biển, tr 98-105 Nhà xuất Khoa học Tự nhiên Cơng nghệ Phạm Sỹ Hồn Nguyễn Thọ Sáo, 2009 Mơ dịng chảy vận chuyển trầm tích lơ lửng cửa sơng vịnh Bình Cang- Nha Trang mơ hình tốn chiều Tạp chí Khí tượng Thủy văn, số 584: 42-50, 8/2009 Phạm Sỹ Hồn, Nguyễn Chí Cơng, Lê Đình Mầu, 2013 Đặc điểm khí tượng, thủy văn động lực vùng biển vịnh Quy Nhơn Tạp chí khoa học cơng nghệ biển, 1(T.13): 1- 11 Phạm Thu Hương, Nguyễn Bá Quỳ Ngơ Lê Long, 2011 Ứng dụng mơ hình MIKE 21 FM nghiên cứu ảnh hưởng sóng dịng chảy đến cửa sơng Đà Rằng tỉnh Phú n Tạp chí Khoa học Cơng nghệ Hàng hải, số 27: 42- 46 8/2011 Sở Khoa học Công nghệ tỉnh Bình Định, 2004 Đặc điểm Khí hậu - Thủy văn tỉnh Bình Định 20 ... dịng chảy tồn khu vực có hướng chảy vào Khu vực Khi triều bắt đầu xuống, dịng chảy phía đỉnh đầm có hướng chảy vào, dịng chảy khu vực cửa đầm lại có hướng chảy Trong thời điểm này, dòng chảy. .. độ dòng chảy lớn tập trung khu vực cửa sơng (phía đỉnh đầm) khu vực cửa đầm (bắt đầu đổ biển) Do đặc điểm hình dạng qn tính dịng chảy mà dịng chảy có đặc điểm lệch phía bờ tây khu vực đỉnh đầm. .. lưu lượng nước sông vào mùa khô xác định theo đặc điểm dịng chảy mùa khơ chiếm khoảng 30% năm (Đặc điểm Khí hậu Thủy văn Bình Định, 2004) Lưu lượng trung bình ngày thời kỳ mưa lũ điển hình trạm