Các phân tích trường gió trung bình từ năm 1979 - 8/2015 đã cho thấy rằng trong vịnh Vũng Rô, quá trình thủy động lực chịu ảnh hưởng chính bởi dòng triều khi mà tần suất gió yếu chiếm tỷ lệ khá cao và chế độ gió chịu ảnh hưởng hoàn toàn bởi tính địa phương khu vực và ít có khả năng thay đổi đáng kể tốc độ dòng triều. Từ phân tích tác động của dòng triều, cơ chế dòng vào - ra trong vịnh khá đặc trưng.
Tạp chí Khoa học Cơng nghệ Biển; Tập 17, Số 2; 2017: 121-131 DOI: 10.15625/1859-3097/17/2/9249 http://www.vjs.ac.vn/index.php/jmst TÍNH TỐN DỊNG CHẢY TẠI VỊNH VŨNG RÔ THEO PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Trần Văn Chung, Nguyễn Hữu Huân* Viện Hải dương học, Viện Hàn lâm Khoa học Công nghệ Việt Nam * E-mail: nghhuan@gmail.com Ngày nhận bài: 23-2-2017 TÓM TẮT: Các phân tích trường gió trung bình từ năm 1979 - 8/2015 cho thấy vịnh Vũng Rơ, q trình thủy động lực chịu ảnh hưởng dịng triều mà tần suất gió yếu chiếm tỷ lệ cao chế độ gió chịu ảnh hưởng hồn tồn tính địa phương khu vực có khả thay đổi đáng kể tốc độ dịng triều Từ phân tích tác động dịng triều, chế dòng vào - vịnh đặc trưng Sự tương đồng độ lớn ngược hướng hai pha triều thể rõ ràng mơ Ảnh hưởng trường gió Đơng Bắc thể rõ pha triều xuống, có xuất vài xoáy cục nhỏ, đáng ý xốy thuận cục phía đơng bắc vịnh Có thay đổi vị trí để tốc độ dòng đạt giá trị lớn tốc độ dòng lớn gia tăng thêm khoảng 0,5 cm/s với hướng lệch 5,3o theo chiều kim đồng hồ Tác động gió mùa Tây Nam ảnh hưởng đến phân bố dòng chảy pha triều lên, phía vịnh hình thành xốy nghịch cục Các xốy hình thành có tác dụng làm suy giảm tốc độ dòng đạt cực trị, suy giảm khoảng 1,7 cm/s không làm lệch hướng dòng chảy đạt giá trị lớn (chỉ lệch 0,3o theo chiều kim đồng hồ) vị trí dịng đạt giá trị lớn Từ khóa: Thủy triều, dịng chảy, mơ hình hai chiều phi tuyến, phương pháp phần tử hữu hạn, vịnh Vũng Rô MỞ ĐẦU Vũng Rơ vịnh nhỏ thuộc xã Hịa Xn Nam, huyện Đơng Hịa, tỉnh Phú n, nằm sát rìa dãy núi Đèo Cả Vịnh ranh giới tự nhiên biển Phú n với Khánh Hịa Vũng Rơ nằm tiếp giáp với biển Đại Lãnh thuộc vịnh Vân Phong, tỉnh Khánh Hịa Vịnh Vũng Rơ có diện tích 16,4 km² mặt nước, dãy núi cao che chắn Đèo Cả, Đá Bia Hòn Bà từ phía bắc, đơng tây Phía nam vịnh đảo Hịn Nưa cao 105 m Vũng Rơ ba địa điểm có điều kiện tự nhiên tốt Việt Nam để xây dựng cảng biển lớn (hai địa điểm lại Cam Ranh Vân Phong) Hơn việc nằm cạnh cảng trung chuyển container quốc tế Vân Phong tạo cho Vũng Rô lợi lớn cho hoạt động xuất nhập hàng hóa, giao lưu với giới Nhằm khai thác lợi cảng Vũng Rô, khu kinh tế Nam Phú Yên thành lập Các hoạt động điều tra khảo sát, nghiên cứu, tính tốn yếu tố mơi trường, sinh thái, thủy văn, động lực học cho khu vực trở nên cấp bách cần thiết phục vụ cho công tác bảo vệ môi trường phát triển bền vững Các nghiên cứu trình động lực nói chung chế độ dịng chảy nói riêng cho Vũng Rô - Phú Yên phương pháp phần tử hữu hạn mẻ Việt Nam Việc sử dụng phương pháp phần tử hữu hạn khắc phục khó khăn trước mà phương pháp sai phân hữu hạn gặp phải ứng dụng nghiên cứu cho vùng có địa hình đáy phức tạp, vấn đề lưới tính, biên, dịng dọc bờ, Phương 121 Trần Văn Chung, Nguyễn Hữu Huân pháp phần tử hữu hạn giải tốn nước nơng dựa tính linh hoạt lưới hình học cao, đặc biệt nghiên cứu mạng lưới thích ứng khơng cấu trúc vùng có địa hình đáy biến đổi phức tạp Các mơ hình phần tử hữu hạn trở nên phổ biến ứng dụng nhiều nơi giới, kể tên mơ hình điển hình giới thuộc loại gồm ADCIRC, QUODDY, BELLAMY, UTBEST, SHYFEM,… Các phương trình nước nơng sử dụng để mơ q trình thủy động lực học cửa sông ven biển Các phương trình chủ đạo cho động lực hệ thống bao gồm hai phương trình liên tục phương trình động lượng Để đơn giản hai phương trình trên, sử dụng giả thiết sau: (i) mật độ không đổi; (ii) áp suất thẳng đứng thủy tĩnh; (iii) hồn lưu bình lưu lớn nhiều so với hồn lưu đối lưu MƠ HÌNH HĨA CÁC PHƯƠNG TRÌNH THỦY ĐỘNG LỰC BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Hai phương trình viết dạng sau: Phương pháp phần tử hữu hạn việc tính tốn mơ q trình Hải dương học hướng nghiên cứu mẻ giới, tiếp tục nghiên cứu hồn thiện Cụ thể, cơng trình Comblen nnk., (2009) [1] mơ mơ hình hồn lưu chung đại dương theo cách hiệu tới rời rạc hóa phương trình vi phân phần bề mặt cong phương tiện phương pháp phần tử hữu hạn mạng lưới tam giác Jones Davies (2010) [2], ứng dụng mô hình theo phương pháp phần tử hữu hạn vào nghiên cứu dòng triều Mersey Estuary Eastern Irish Sea,… Bajo nnk., (2015) [3] áp dụng mơ hình số trị dựa kỹ thuật rời rạc hóa phần tử hữu hạn (SHYFEM (Shallow Water Hydrodynamic Finite Element Model)) sử dụng nghiên cứu hoàn lưu nước thềm phía tây bắc Biển Đen, đặc biệt gần Danube Delta ven bờ Rumani Các cơng trình nghiên cứu số trị cho mơ hình thủy động lực học theo phương pháp phần tử hữu hạn, tập trung cơng trình nghiên cứu Bùi Hồng Long Trần Văn Chung (2007, 2008, 2009, 2010, 2012, 2013, 2014) [4-10], nhóm tác giả Trần Văn Chung Tống Phước Hồng Sơn (2014) [11] áp dụng nghiên cứu chế độ động lực - sinh địa hóa vùng ven biển Việt Nam bước đầu thực so sánh kết nghiên cứu mơ hình FEM với mơ hình Ecosmo với thực tế đo đạc vùng nghiên cứu Bình Cang - Nha Trang [12] so sánh với số liệu thực tế trạm mực nước nghiên cứu chế độ dòng chảy cho vịnh Bắc Bộ [13] 122 Phương trình liên tục: H t Hv (1) Phương trình động lượng nằm ngang: cd W f v v v g v v t H H v (2) Trong đó: v vận tốc lấy trung bình theo độ sâu (m/s); t thời gian (s); toán tử vi phân gradient nằm ngang; độ nâng bề mặt so với mực nước tĩnh (m); g gia tốc trọng trường (m/s2); f quay Coriolis (s-1); cd hệ số cản đáy; H tổng độ sâu cột nước (với dâng mực nước: H dz dz (h ) h với rút mực nước: H dz (h ) h h độ sâu so với mực nước tĩnh (m)) Trong trường hợp xét đến dịng chảy gió: Đưa vào giá trị vận tốc gió tồn mạng lưới tính Giá trị dùng để giải phương trình động lượng thơng qua W - ứng suất gió động học (ứng suất chia với mật độ, W τ , τ ứng suất gió, mật độ nước biển) (Pascal) Tính tốn dịng chảy vịnh Vũng Rơ… τ K aV10 V10 (3) Trong đó: V10 vận tốc gió 10 m so với bề mặt biển, K hệ số cản bề mặt biển, a mật độ khơng khí (a = 1,25 kg.m-3) Có nhiều cơng trình nghiên cứu hệ số K, mơ hình tính chúng tơi sử dụng theo tính tốn WAMDI Group (1988) dạng: 0,0012875 V10 7,5 m / s K (4) 0,0008 0,000065 V V10 7,5 m / s 10 Trong trường hợp có xét đến dịng triều: Dao động thủy triều đưa vào nút giá trị biên mở Để mà tạo tổng hợp lực triều thích hợp áp đặt xác điều kiện biên Dirichlet Phương trình sau phải lấy tổng tất thành phần triều: (t ) Z f k (t ) Ak cos k t t V (t ) U (t ) k gk 180 k Trong đó: (t) độ cao thủy triều tổng hợp vị trí biết theo thời gian (m); Z0 giá trị mực nước biển trung bình điểm cho mực nước “khơng độ sâu” (m); số k biểu thị sóng triều riêng biệt, fk V (t ) U (t ) k tham số thiên văn phụ thuộc vào thời gian quan trắc, fk gọi nhân tử biên độ, V (t ) U (t ) k gọi pha thiên văn, với V(t0) điều chỉnh thời gian vận hành U(t) thừa số thay đổi pha; Ak, gk số điều hòa biên độ (m) pha trễ Greenwich (độ) thành phần, chúng phụ thuộc vào điều kiện địa phương vị trí quan trắc; k tần số thủy triều thành phần (rad/s) Sử dụng hệ số ma sát đáy cd: Để đưa vào hệ số này, thơng thường có ba kiểu lựa chọn: đưa vào hệ số cản không đổi (CD), xấp xỉ Manning (n) xấp xỉ Chezy Các hệ số dùng cho tính tốn ứng suất đáy phương trình động lượng Ma sát đáy CD: CD hệ số cản không đổi sử dụng theo cách đơn giản để xác định ứng suất lực cản ma sát mà tạo địa hình đáy Lực cản ma sát giá trị không đổi dựa độ cao nước Chỉ giá trị mà xác định độ cản không đổi điểm liên quan độ ghồ ghề địa hình kết hợp với điểm Lựa chọn CD đưa vào trực tiếp hệ số ứng suất đáy cd Vì vậy: cd = CD (6) (5) Theo Phạm Văn Ninh [14], mơ hình chọn cd = 0,0026 Hệ số đưa vào chương trình cách tạo mảng CD cho tất nút Ma sát đáy Manning: Kiểu thứ hai ma sát đáy sử dụng công thức Manning tính tốn ứng suất đáy Manning sử dụng thực nghiệm để xác định hệ số Chezy tỉ lệ thuận với bậc sáu bán kính thủy lực tỉ lệ nghịch với hệ số Manning, n Hệ số Manning phụ thuộc tham số thực nghiệm độ sâu nước, số Reynold, mặt cắt hình học ngang, vật chất đáy thảm thực vật đáy Để mà tính cd, hệ số ứng suất đáy, cơng thức sau sử dụng theo hệ đơn vị mét: cd gn R1/3 (7) Trong đó: g gia tốc trọng trường (m/s2); n hệ số ghồ ghề Manning; R bán kính thủy lực theo mét Bán kính thủy lực eo biển với mặt cắt ngang vng góc xác định sau: R bH 2 H b (8) Trong đó: b độ rộng eo biển với mặt cắt ngang vng góc; H độ sâu eo biển 123 Trần Văn Chung, Nguyễn Hữu Huân Trong phần lớn cửa sông ven biển, độ rộng phải lớn nhiều độ sâu, bán kính thủy lực xác định sau: R=H (9) g CD Chezy = Phương trình (7) trở thành: cd Để mà chuyển CD xấp xỉ Manning vào định dạng Chezy, công thức sau sử dụng: gn H 1/3 (10) Các hệ số sử dụng chương trình tùy chọn theo kiểu đưa vào nói Thứ tự mơ hình tính đưa vào hệ số Manning ứng suất đáy Manning chuyển tới ứng suất đáy CD sau chuyển tới ứng suất đáy Chezy Bảng Các hệ số Manning cho kiểu đáy xác định Kiểu đáy N Sỏi (Gravel) Sỏi đá cuội (Gravel and boulders) Đất (Earth) Đất cỏ biển (Earth and grass) Đất, nhiều cỏ biển (Earth, very weedy) 0,025 0,040 0,030 0,026 0,080 (13) 1/ Chezy = H n (14) Hệ thống phi tuyến phương trình chủ đạo mơ hình giải theo phép lặp bước thời gian Tại điểm bắt đầu phép lặp, trạng thái hệ thống xem xét để xác định rõ sở dạng vật lý Trong q trình tính, phần tử mà thuộc nút với độ sâu 0,5 m nhỏ áp đặt cho chi phối trình động học trình bày Ip nnk., (1988) [15], tồn phần tử lại áp đặt cho chi phối q trình động lực học Cơng thức hóa hai q trình vật lý trình bày chi tiết bên Mực nước tónh z=0 Trong trường hợp nghiên cứu chế độ dòng chảy vịnh, có đáy khơng phức tạp, ta chọn hệ số Manning sau: n = 0,025 (cho vịnh), 0,075 (vùng cửa sông) Ma sát đáy Chezy: Cách tiếp cận thứ ba tới ma sát đáy cách tiếp cận Chezy Lực cản đáy Chezy phụ thuộc tỉ lệ thuận với chu vi eo biển ướt bình phương vận tốc tỉ lệ nghịch với độ dốc thủy lực diện tích mặt cắt ngang qua eo biển Lực cản (resistance) V P AS (11) Độ cao bề mặt H0 Lơ ùp x ốp Đáy biển rắn z=-h z=-(h+h0) Hình Thể dạng hình học mơ hình tính Đây biểu thức phương trình Chezy: V = C RS (12) V Trong đó: C hệ số Chezy; R lực cản Như nhìn thấy trên, kiểu ma sát đáy tương tự công thức Manning phụ thuộc vào độ sâu Vì điểm ước lượng ma sát đáy dựa điểm hệ số Chezy xác định độ cao nước 124 z= Hình Cách bố trí mạng lưới tam giác cho phần tử Tính tốn dịng chảy vịnh Vũng Rơ… Mơ hình nhóm tác giả Bùi Hồng Long Trần Văn Chung ứng dụng thành cơng vào tính dịng triều xác định số điều hòa sóng triều cụm đảo Song Tử [4], vịnh Cam Ranh [5], Vân Phong [6], Đầm Bấy (vịnh Nha Trang) [10], Bình Cang - Nha Trang [12] dịng chảy tổng hợp vịnh Bắc Bộ [13] tính tốn 0,1 m, bước thời gian 100 s, số vòng lặp bước thời gian 100, hàm trọng số =1 (sai phân theo bước thời gian hồn tồn ẩn tốn ổn định khơng điều kiện) KẾT QUẢ NGHIÊN CỨU CHẾ ĐỘ DỊNG CHẢY TẠI VŨNG RƠ Các thơng tin liệu đầu vào Thơng tin chung Hình Mạng lưới tam giác cho nghiên cứu chế độ dòng chảy Hình Trường độ sâu (m) làm khớp đồ Google Earth Để tính dịng chảy cho Vũng Rơ, chúng tơi thiết lập mạng lưới tính với kinh độ từ 109,38629oE đến 109,43760oE, vĩ độ từ 12.83287oN đến 12,88080oN (hiệu chỉnh lại đồ Google Earth thể hình 3) Mạng lưới tam giác thiết lập với góc cực tiểu 30o, diện tích cực đại 10.000 m2, trung bình 6.248 m2, cực tiểu 2.807 m2 (trên hình hiệu chỉnh lại Google Earth hình 5) Trong đó, diện tích mặt thống cho tính tốn 12,75 km2, tương ứng với 1.106 điểm nút nằm ngang 2.041 lưới tam giác Các nút gán để đưa vào điều kiện biên mở cho dao động thủy triều 32 nút, biên mở mặt cắt A (gọi Biên A) nút mặt cắt B (gọi Biên B) 23 nút, xem chi tiết hình Độ sâu cực tiểu Hình Mạng lưới tam giác làm khớp Google Earth Để cập nhật số liệu dòng chảy yếu tố khí tượng Trong khn khổ đề tài mã số VAST 06.04/14-15 thực chuyến khảo sát bổ sung vào tháng 5-6/2014 11/2014 Trên hình trạm vị khảo sát vật lý - mơi trường, trạm thực đo 125 Trần Văn Chung, Nguyễn Hữu Huân mặt rộng, riêng trạm ký hiệu “B” đo thêm liên tục ngày đêm, dòng chảy phút/số liệu liệu theo ốp giờ/số liệu theo ngày 1, 7, 13, 19 Khu vực Tuy Hịa - Phú n (109o17’E; 13o05’N) Từ hình 8, thấy ba hướng gió chiếm ưu bắc đông bắc (NNE) với tần suất xuất 19,31%, thời gian trì liên tục theo gió trung bình ngày đạt 36 ngày; đơng bắc (NE) chiếm 16,35%; hướng bắc (N) chiếm 10,57% Các số liệu thể đặc trưng trường gió mùa Đơng Bắc địa phương Các hướng lại tần suất xuất 10% Trong gió mùa Tây Nam đặc trưng hướng gió tây (W) chiếm 10,01% Phân tích gió Vũng Rơ (109,423041oE; 12,867792oN) N NNW Tốc độ gió (m/s) NNE NW NE Hình Các trạm đo cho hiệu chỉnh mơ hình WNW Điều kiện biên thủy triều: Biên mở (A) (B) (hình 2) thể theo dao động mực nước triều hình ENE W 2 - >4 - >6 - >8 - 10 >10 - 12 >12 - 14 >14 - 16 >16 - 18 >18 - 20 >20 E 0% 4% 8% 12% 16% WSW 20% ESE SW SE SSW SSE S Hình Hoa gió khu vực Tuy Hịa, Phú n N T?c d? gió (m/s) NNW Hình Giá trị dao động mực nước điển hình cho tính dịng triều cho Vũng Rơ NNE NW NE WNW Để có thơng tin chế độ gió Vũng Rơ, chúng tơi sử dụng thơng tin gió cung cấp từ NCEP CFSR từ năm 1979 đến 8/2015 theo giờ/số liệu Để đối chứng hiệu chỉnh thực tế mang tính địa phương gió Vũng Rơ, chúng tơi sử dụng số liệu gió gần khu vực Vũng Rơ nhất, trạm Tuy Hịa (6) thơng tin từ đợt khảo sát khí tượng vào tháng 5-6/2014 11/2014 Với chế độ gió trạm Tuy Hịa, chúng tơi sử dụng nguồn số liệu gió từ năm 1987 đến 2007, tần suất đo số 126 ENE W E 0% 4% 8% 12% WSW 2 - >4 - >6 - >8 - 10 >10 - 12 >12 - 14 >14 - 16 >16 - 18 >18 - 20 >20 16% ESE SW SE SSW SSE S Hình Hoa gió Vũng Rơ theo số liệu NCEP CFSR (1979 - 8/2015) Tính tốn dịng chảy vịnh Vũng Rơ… Sử dụng liệu gió cung cấp từ NCEP, CFSR từ năm 1979 đến 8/2015 theo giờ/số liệu, phân tích chúng tơi thấy tương đồng với số liệu gió Tuy Hịa, chi tiết có khác biệt đáng ý nghiên cứu đặc trưng vùng Vũng Rơ Theo kết phân tích tốc độ gió Vũng Rơ tương đối yếu, tần suất gió chiếm 36,8% tập trung tốc độ gió ≤ ws < Từ thơng tin gió thể bảng phương diện tính trung bình nhiều năm, tốc độ gió cao vào tháng 12 thấp vào tháng Tuy nhiên, theo kết phân tích tốc độ gió cao vùng đạt 22,7 m/s, hướng tây tây bắc vào tháng 11 (20 h ngày 11/11/2011) vào tháng không thấy xuất vận tốc cao 10 cm/s (bảng 2) Hướng gió có tần suất xuất nhiều hướng bắc (chiếm 15,7%) sau tới hướng NNE (chiếm 13,8%), hai hướng chịu ảnh hưởng trường gió Đơng Bắc tác động đến khu vực Tần suất xuất đứng thứ ba gió theo hướng W (9,0%) chịu ảnh hưởng trường gió Tây Nam Bảng Tốc độ gió lớn xảy Vũng Rô Tháng 10 11 12 13,7 13,6 14,3 10,3 11,0 11,5 9,6 13,6 15,7 17,1 22,7 17,7 0,0 45,0 0,0 157,5 337,5 270,0 270,0 270,0 247,5 0,0 292,5 22,5 Thời gian 13 h 26/1/ 2006 22 h 23/2/ 2013 4h 5/3/ 2005 15 h 9/4/ 1983 20 h 14/5/ 2006 0h 6/6/ 1999 13 h 29/7/ 2009 7h 16/8/ 2002 0h 29/9/ 2009 5h 30/10/ 2010 20 h 11/11/ 2011 11 h 4/12/ 2006 Vtb (m/s) 5,8 4,6 4,0 3,7 3,6 3,9 3,8 3,9 3,4 4,2 5,8 6,8 Vmax (m/s) Hướng o () Dòng chảy ảnh hưởng chế độ gió mùa Đối với dịng chảy ảnh hưởng gió triều, trung bình vùng tính, chế độ gió cụ thể sau: cho thấy, trường gió mùa Đông Bắc tác động đến vịnh cách tạo vài xốy cục vịnh, nhìn chung có tác động thay đổi tốc độ dịng (hình 10) Do ảnh hưởng trường gió Đơng Bắc: Tốc độ gió 5,8 m/s, hướng N Do ảnh hưởng trường gió Tây Nam: Tốc độ gió 3,8 m/s, hướng W Gió mùa Đơng Bắc Pha triều xuống Dưới ảnh hưởng trường gió Đơng Bắc, xuất vài xốy cục nhỏ, đáng ý xoáy thuận cục phía đơng bắc vịnh, có tâm nằm vị trí (109,42696oE, 12,87108oN), cịn vị trí khác khơng có thay đổi đáng kể chịu tác động triều Có thay đổi vị trí để tốc độ dòng đạt giá trị lớn tốc độ dòng gia tăng khoảng 0,5 cm/s với hướng lệch 5,3o theo chiều kim đồng hồ Cụ thể, tốc độ dịng đạt 44,5 cm/s, hướng 67,1o vị trí (109,38912oE; 12,83792oN), độ sâu 14,8 m Từ kết Hình 10 Phân bố dịng chảy trung bình theo độ sâu cho pha triều xuống ảnh hưởng trường gió Đơng Bắc 127 Trần Văn Chung, Nguyễn Hữu Huân Pha triều lên Trong pha triều xuống ta thấy rõ thay đổi dịng bên vịnh pha triều lên thay đổi không rõ ràng Cơ chế tạo xoáy cục cho pha triều lên gần triệt tiêu Trong pha triều lên, ảnh hưởng trường gió làm suy giảm tốc độ dòng, cụ thể cho tốc độ 31,4 cm/s (giảm 0,7 cm/s) hướng 239,4o (lệch 2,3o theo chiều ngược kim đồng hồ) Tuy vậy, tốc độ dòng đạt giá trị lớn trùng với vị trí dịng triều đạt lớn nhất, xung quanh vị trí (109,38982oE, 12,83756oN), độ sâu 2,1 m (hình 11) vịnh hình thành xoáy nghịch cục bộ, mà rõ nét hai xốy nghịch có tâm (109,42820oE, 12,87141oN) (109,41952oE, 12,86288oN) Các xốy hình thành có tác dụng làm suy giảm tốc độ dòng đạt cực trị, suy giảm khoảng 1,7 cm/s không làm lệch hướng dòng chảy đạt giá trị lớn (chỉ lệch 0,3o theo chiều kim đồng hồ) vị trí dịng đạt giá trị lớn (hình 13) Hình 12 Phân bố dịng chảy trung bình theo độ sâu cho pha triều xuống ảnh hưởng trường gió Tây Nam Hình 11 Phân bố dịng chảy trung bình theo độ sâu cho pha triều lên ảnh hưởng trường gió Đơng Bắc Gió mùa Tây Nam Pha triều xuống Tác động gió mùa Tây Nam đến vịnh cho pha triều xuống không đáng kể Ảnh hưởng thấy vị trí dịng có tốc độ yếu, nhiên tác động gió mùa Tây Nam làm tăng tốc độ dòng cực trị khoảng 1,2 cm/s, hướng lệch 0,2o theo hướng ngược kim đồng hồ không làm thay đổi vị trí dịng đạt cực trị (hình 12) Pha triều lên Đối với pha triều lên, tác động trường gió mùa Tây Nam rõ ràng, phía 128 Hình 13 Phân bố dịng chảy trung bình theo độ sâu cho pha triều lên ảnh hưởng trường gió Tây Nam Tính tốn dịng chảy vịnh Vũng Rơ… So sánh kết tính với thực tế khảo sát Để hiệu chỉnh kết tính tốn với số liệu thực đo dịng chảy, chúng tơi sử dụng số liệu chuyến khảo sát bổ sung vào tháng 5-6/2014 11/2014 Vị trí cụ thể trạm so sánh thể hình sai số mơ hình thể bảng bên Bảng So sánh kết tính số liệu khảo sát Trạm B1 B2 B1 B2 Đo đạc Tính tốn o (*) Sai số tương đối (%) o Thời điểm V (cm/s) Hg ( ) Tốc độ Hướng 129,709 5,0 128,8 0,6 0,7 h 2/6/2014 66,198 11,7 65,9 0,2 0,5 22 h 1/6/2014 V (cm/s) Hg ( ) Dòng nhỏ 5,030 Dòng lớn 11,720 Dịng trung bình 8,120 8,5 4,5 Dịng nhỏ 1,425 242,163 1,2 240,1 19,2 0,9 13 h 2/6/2014 Dòng lớn 37,925 89,895 37,0 85,4 2,5 5,3 22h 2/6/2014 Dịng trung bình 17,553 Dịng nhỏ 13,283 338,627 13,1 337,5 1,4 0,3 18 h 26/11/2014 Dòng lớn 25,192 10,960 24,2 9,4 4,1 16,6 h 27/11/2014 Dịng trung bình 20,117 Dịng nhỏ 4,703 270,598 4,1 295,0 14,7 8,3 h 28/11/2014 Dòng lớn 38,836 263,231 38,0 269,9 2,2 2,5 h 28/11/2014 Dòng trung bình 26,651 17,5 0,3 19,8 25,9 1,6 2,9 Ghi chú: V: tốc độ dòng chảy; Hg: Hướng dòng chảy; (*): Dòng chảy đo máy đo dòng COMPACT EM, Alec Electronics Co., LTD (Nhật Bản), trung bình phút/số liệu Theo kết phân tích trên, thấy tính tốn thực chịu ảnh mùa gió Tây Nam có sai số tính toán thấp so với thời điểm chịu ảnh hưởng mùa gió Đơng Bắc Số liệu tính tốn cho thấy vị trí trạm liên tục B2 cho kết sai số cao so với trạm B1 KẾT LUẬN Các phân tích trường gió trung bình từ năm 1979 - 8/2015 cho thấy vịnh Vũng Rô, q trình thủy động lực chịu ảnh hưởng dịng triều mà tần suất gió yếu chiếm tỷ lệ cao chế độ gió chịu ảnh hưởng hồn tồn tính địa phương khu vực có khả thay đổi đáng kể tốc độ dịng triều Từ phân tích tác động dịng triều, chế dòng vào - vịnh đặc trưng Sự tương đồng độ lớn ngược hướng hai pha triều thể rõ ràng mô Ảnh hưởng trường gió Đơng Bắc thể rõ pha triều xuống, có xuất vài xoáy cục nhỏ, đáng ý xốy thuận cục phía đơng bắc vịnh, có tâm nằm vị trí (109,42696oE, 12,87108oN) Có thay đổi vị trí để tốc độ dịng đạt giá trị lớn tốc độ dòng lớn gia tăng thêm khoảng 0,5 cm/s với hướng lệch 5,3o theo chiều kim đồng hồ Tác động gió mùa Tây Nam ảnh hưởng đến phân bố dịng chảy pha triều lên, phía vịnh hình thành xốy nghịch cục bộ, mà rõ nét hai xốy nghịch có tâm (109,42820oE, 12,87141oN) o o (109,41952 E, 12,86288 N) Các xốy hình thành có tác dụng làm suy giảm tốc độ dịng đạt cực trị, suy giảm khoảng 1,7 cm/s khơng làm lệch hướng dịng chảy đạt giá trị lớn (chỉ lệch 0,3o theo chiều kim đồng hồ) vị trí dịng đạt giá trị lớn Từ phân tích trên, sử dụng mơ hình theo phương pháp phần tử hữu hạn với lưới phi cấu trúc (mạng lưới tam giác), tìm vị trí có phân bố dịng chảy tương đối đặc biệt (các xốy cục bộ) Việc 129 Trần Văn Chung, Nguyễn Hữu Huân kiểm nghiệm tính đắn phương pháp phần tử với thực tế đo đạc giúp hiệu chỉnh lại thông số tính tốn cần thiết, phục vụ tốt cho mơ tốn dịng chảy biển Nếu việc chỉnh lý số liệu đo đạc thực đồng bộ, chi tiết, đáng tin cậy cung cấp tốt cho liệu đầu vào cho mơ hình thu kết tính mang tính định lượng cao vị trí cần quan tâm Kết mơ hình giúp đưa giải pháp, thơng số kỹ thuật tương đối xác để nhà quản lý có sách hoạch định, quy hoạch cơng trình - dịch vụ biển cách hợp lý, tiết kiệm tránh lãng phí khơng cần thiết, góp phần hạn chế tai biến thiên nhiên, Lời cảm ơn: Chúng xin gởi lời cảm ơn chân thành đến chủ nhiệm đề tài mã số VAST 06.04/14-15 “Đánh giá khả tự làm vịnh Vũng Rô (Phú Yên) phục vụ phát triển bền vững kinh tế biển” đồng nghiệp nhóm nghiên cứu góp ý hỗ trợ giúp chúng tơi hoàn thành báo TÀI LIỆU THAM KHẢO Comblen, R., Legrand, S., Deleersnijder, E., and Legat, V., 2009 A finite element method for solving the shallow water equations on the sphere Ocean Modelling, 28(1), 12-23 Jones, J E., and Davies, A M., 2010 Application of a finite element model to the computation of tides in the Mersey Estuary and Eastern Irish Sea Continental Shelf Research, 30(5), 491-514 Bajo, M., Ferrarin, C., Dinu, I., Umgiesser, G., and Stanica, A., 2014 The water circulation near the Danube Delta and the Romanian coast modelled with finite elements Continental Shelf Research, 78, 62-74 Bùi Hồng Long, Trần Văn Chung, 2007 Tính tốn dịng triều cụm Song Tử phương pháp phần tử hữu hạn Tuyển tập Báo cáo Hội nghị Quốc gia “Biển Đông2007”, ISSN 1859-2430 Tr 735-750 Bùi Hồng Long, Trần Văn Chung, 2008 Kết mơ chế độ dịng triều vịnh Cam Ranh phương pháp phần tử 130 hữu hạn Tạp chí Khoa học Công nghệ biển, 8(4), 19-35 Long, B H., and Chung, T V., 2009 Calculations of tidal currents in Van Phong bay using the finite element method Advances in Natural Science, 10(4), 495-478 Bui Hong Long, Tran Van Chung, 2010 Some experimental calculation for 3D currents in the strong upwelling region of southern central Vietnam using finite element method Proceedings of the International Conference marine biodiversity of east asian seas: status, challenges and sustainable development Nha Trang, Vietnam, 165-177 Bùi Hồng Long, Trần Văn Chung, 2012 Nghiên cứu chế độ dòng chảy vịnh Phan Thiết mơ hình ba chiều phi tuyến với phương pháp phần tử hữu hạn Tạp chí Khoa học Công nghệ biển, 12(4), 1-14 Bùi Hồng Long, Trần Văn Chung, 2013 Thử nghiệm tính tốn hệ thống dịng chảy khu vực biển Nam Trung Bộ mơ hình ba chiều (3D) phi tuyến Kỷ yếu Hội nghị Quốc tế “Biển Đông 2012”, Nha Trang, 12-14/09/2012, 17-28 10 Bùi Hồng Long, Trần Văn Chung, 2014 Tính tốn dịng chảy triều khu vực Đầm Bấy (vịnh Nha Trang) phương pháp phần tử hữu hạn Tạp chí Khoa học Công nghệ biển, 14(4), 332-340 11 Tran Van Chung, Tong Phuoc Hoang Son, 2014 The numerical simulations on hydrodynamic and bio-geochemistry processes in Vietnam sea waters In Proceedings of International Mini Workshop on the Western Pacific Marine Biogeochemical Environment Variability Jamstec, Tokyo, - 4, February, 2014 45-47 12 Trần Văn Chung, Bùi Hồng Long, 2014 Đặc trưng thủy động lực vực nước Bình Cang - Nha Trang qua mơ hình FEM ECOSMO Tạp chí Khoa học Công nghệ biển, 14(4), 320 - 331 13 Trần Văn Chung, Bùi Hồng Long, 2015 Tính tốn dịng chảy vịnh Vũng Rơ… Một số kết tính tốn dịng chảy vịnh Bắc Bộ mơ hình ba chiều phi tuyến Tạp chí Khoa học Cơng nghệ biển, 15(4), ISSN 1859-3097, 320 - 333 14 Chương trình điều tra nghiên cứu biển cấp nhà nước KHCN-06 (Phạm Văn Ninh (chủ biên) (1996-2000), 2003 Biển Đơng (phần khí tượng thủy văn động lực biển) Tập Nxb Đại học quốc gia Hà Nội, 565 tr 15 Ip, J T C., Lynch, D R., and Friedrichs, C T., 1998 Simulation of estuarine flooding and dewatering with application to Great Bay, New Hampshire Estuarine, Coastal and Shelf Science, 47(2), 119-141 CALCULATIONS OF CURRENT IN THE VUNG RO BAY USING THE FINITE ELEMENT METHOD Tran Van Chung, Nguyen Huu Huan Institute of Oceanography, VAST ABSTRACT: The analyses the average wind field in the period from 1979 to August 2015 have demonstrated that hydrodynamical processes are dominated by tidal currents in the Vung Ro bay In this period, the frequency of weak wind accounts for a quite high percentage Moreover, the wind field is totally influenced by local conditions, and it is less likely to change the speed of tidal currents The results from a simulation have clearly illustrated that there are similarities between speed and direction in two tidal phases by analyzing the influence of tidal current and the regime of in/out flow in the bay In addition, the influence of the northeast wind is most obvious at ebb-tide phase, a few small local vortices also appear, notably a local cyclone vortex is located at the northeast of the bay There are location changes of currents to achieve the maximum value, those velocities increase by around 0.5 cm/s with the clockwise direction of 5.3 degree The effect of southwest wind has impacted on distribution of currents in flood-tide phase, then the local anticyclonic eddies are formed inside the bay These eddies have decreased flow rate at a peak point (about 1.7 cm/s) but not prevented the direction and loction of currents from reaching the maximum value (the clockwise direction of less than 0.3 degree) Keywords: Tide, current, two-dimensional (2D) nonlinear model, finite element method (FEM), Vung Ro 131 ... với hồn lưu đối lưu MƠ HÌNH HĨA CÁC PHƯƠNG TRÌNH THỦY ĐỘNG LỰC BẰNG PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN Hai phương trình viết dạng sau: Phương pháp phần tử hữu hạn việc tính tốn mơ q trình Hải dương học... dương theo cách hiệu tới rời rạc hóa phương trình vi phân phần bề mặt cong phương tiện phương pháp phần tử hữu hạn mạng lưới tam giác Jones Davies (2010) [2], ứng dụng mơ hình theo phương pháp phần. .. chế độ dòng chảy vịnh Phan Thiết mơ hình ba chiều phi tuyến với phương pháp phần tử hữu hạn Tạp chí Khoa học Cơng nghệ biển, 12(4), 1-14 Bùi Hồng Long, Trần Văn Chung, 2013 Thử nghiệm tính tốn