Tổng quan luận án

Luận án Tiến sĩ Khoa học ngành Địa lý – Địa chất, chuyên ngành Hải dương học (mã số: 010707) với đề tài "Dao động tự do và dao động mùa của mực nước Biển Đông" được tác giả Phạm Văn Huấn thực hiện tại Trường Đại học Tổng hợp Hà Nội (Bộ Giáo dục và Đào tạo), dưới sự hướng dẫn khoa học của GS.TS. Nguyễn Ngọc Thụy và bảo vệ năm 1993.

Về tính cấp thiết, biến động thời gian và không gian của mực nước biển là hiện tượng tự nhiên quy mô lớn, tác động trực tiếp tới các hoạt động kinh tế kỹ thuật biển như vận tải hàng hải, xây dựng công trình ven bờ và ngoài khơi, hệ thống đê điều bảo vệ bờ, tưới tiêu nông nghiệp ven biển, cấp thoát nước đô thị và phòng chống thiên tai ngập lụt do bão, nước dâng. Về mặt khoa học, dao động mực nước phản ánh tương tác động lực học phức tạp giữa khí quyển, thủy vực và địa hình đáy biển.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ thực trạng các công trình trước đây tại Biển Đông chủ yếu tập trung vào thủy triều và nước dâng do bão, trong khi các dao động ở dải tần số khác (dao động mùa, dao động chu kỳ xi nốp, dao động tự do/dao động riêng nhiều giờ) chưa được khảo sát đầy đủ. Đặc biệt, bản chất cơ chế cộng hưởng động lực học giải thích sự ngự trị của chế độ nhật triều và tính phân hóa không gian của mực nước trên toàn Biển Đông chưa có công trình nào giải quyết bằng mô hình số trị dao động riêng trên địa hình thực.

Luận án đặt ra và giải quyết bốn nhiệm vụ nghiên cứu cụ thể:

  1. Khảo sát chi tiết về dao động tự do của Biển Đông;
  2. Giải thích cơ chế hình thành của hiện tượng thủy triều rất phức tạp và độc đáo ở Biển Đông;
  3. Nghiên cứu cấu trúc của những dao động mực nước ở các vùng bờ khác nhau dọc bờ biển nước ta;
  4. Nghiên cứu chế độ dao động mùa của mực nước ở Biển Đông và thử nghiệm mô hình tính mực nước dâng trong gió theo dữ liệu về trường gió và áp suất trên biển.

Đối tượng nghiên cứu là trường độ cao mực nước, các dao động tự do và dao động mùa. Phạm vi không gian bao trùm toàn bộ thủy vực Biển Đông (từ $99^\circ 30'$ Đ đến $121^\circ$ Đ và từ $3^\circ$ N đến $25^\circ$ B với diện tích khoảng $3,5 \cdot 10^6 \text{ km}^2$), bao gồm vịnh Bắc Bộ, vịnh Thái Lan và dải bờ biển phía Tây. Phạm vi thời gian và dữ liệu gồm chuỗi số liệu thực đo mực nước từng giờ kéo dài 1 năm tại các trạm Hòn Dấu, Đà Nẵng, Quy Nhơn, Vũng Tàu, Rạch Giá; trường gió mùa đông bắc và tây nam điển hình; cùng chuỗi mô phỏng số trị 1000 giờ cho dao động tự do.

Tổng quan tài liệu và vị trí của luận án

Tổng quan các công trình nghiên cứu về dao động mực nước Biển Đông được tác giả phân chia theo các nhóm chủ đề:

  • Nhóm nghiên cứu thủy triều: Các nghiên cứu sơ lược đầu thế kỷ XX của các tác giả phương Tây [32, 61]; các nghiên cứu áp dụng phương pháp hiện đại từ các công trình [55, 41]; các bản đồ triều giải tích và số trị của các tác giả trong nước và quốc tế [24, 30, 12, 15, 25, 45, 13]. Công trình của Nguyễn Ngọc Thụy [15] chỉ ra rằng nguyên nhân hình thành chế độ thủy triều Biển Đông là do sóng triều Thái Bình Dương truyền qua các eo biển phía bắc/đông bắc và chịu sự biến đổi căn bản bởi điều kiện địa hình địa phương, chuyển tính chất nhật triều từ thứ yếu thành chủ yếu. Chương trình nghiên cứu cấp nhà nước KT-03-03 tiếp tục hoàn thiện mô hình số tính triều và dòng triều.
  • Nhóm nghiên cứu nước dâng do bão: Nghiên cứu đặc trưng chế độ bão, tần suất và phân bố xác suất nước dâng nguy hiểm [14, 11]; phương pháp hồi quy tương quan phục vụ dự báo nghiệp vụ [4]; các mô hình số trị giải hệ phương trình sóng dài nước nông từ những năm 1980 của Đỗ Ngọc Quỳnh [33], cùng các công trình [18, 1, 10, 2].
  • Nhóm nghiên cứu dao động mùa và chu kỳ dài: Công trình của các tác giả Liên Xô [26] phân tích sóng năm theo phương pháp đẳng độ cao. Công trình [42] đã phê phán bản đồ trong [26] vì sự tồn tại bất hợp lý của đỉnh sóng năm ở trung tâm biển, đồng thời chứng minh dao động mùa chủ yếu do gió mùa kích động các sóng dài tiến lan truyền trên địa hình nước nông. Các nghiên cứu phân tích phổ và phổ tương hỗ giữa khí quyển và mực nước [14, 17] ghi nhận sự xuất hiện của các sóng bội nước nông bậc ba, bậc bốn, bậc năm trong vùng nhật triều mạnh mà sơ đồ Darwin dễ bỏ sót.
  • Nghiên cứu dao động tự do trên thế giới: Garrett [62] nghiên cứu cộng hưởng hệ thống vịnh Fundy-Maine (chu kỳ $13,3 \pm 0,4$ giờ); nghiên cứu dao động riêng Bắc Đại Tây Dương với chu kỳ 9,3 và 14,8 giờ [67]; các nghiên cứu của các tác giả quốc tế [66, 29, 31, 38, 40, 52] dùng phương pháp tích phân số trị để tìm tần số riêng và mốt riêng.

Vị trí của luận án là công trình nghiên cứu hệ thống đầu tiên áp dụng mô hình số trị thủy động phi ma sát trên lưới địa hình thực để xác định toàn bộ phổ dao động tự do của Biển Đông, lấy đó làm cơ sở vật lý giải thích cơ chế cộng hưởng thủy triều, kết hợp phân tích phổ thực đo chuỗi năm tại bờ Tây và mô phỏng trường nước dâng do gió mùa.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án kết hợp phương pháp lý thuyết giải tích, phương pháp phân tích thống kê xử lý chuỗi thời gian và phương pháp mô hình số trị thủy động:

1. Phương pháp phân tích điều hòa mực nước

Xuất phát từ công thức tổng quát của độ cao mực nước thủy triều: $$H_t = A_0 + \sum_{i=1}^r F_i H_i \cos[Q_i t + (V_0 + u)_i - G_i]$$ Trong đó $A_0$ là mực nước trung bình; $Q_i$ là tốc độ thay đổi pha; $F_i, (V_0+u)_i$ là hệ số suy giảm và pha ban đầu thiên văn; $H_i, G_i$ là biên độ và độ lệch pha cần xác định.

Tác giả đã cải tiến dạng tính toán truyền thống (1.2) sang công thức biến đổi (1.3): $$A_i = F_i \cos[Q_i t + (V_0 + u)_i], \quad B_i = F_i \sin[Q_i t + (V_0 + u)_i]$$ $$X_i = H_i \cos G_i, \quad Y_i = H_i \sin G_i$$ Hệ số $A_i, B_i$ được tính trước cho từng sóng $i$ tại từng thời điểm $t$, còn $X_i, Y_i$ là ẩn số được giải bằng phương pháp bình phương nhỏ nhất. Cải tiến này khắc phục sai số do lấy gần đúng thông số thiên văn ở giữa kỳ quan trắc và cho phép xử lý chuỗi số liệu không liên tục (chuỗi dòng chảy ngắt quãng).

2. Phương pháp phân tích phổ trong hải dương học

  • Áp dụng cặp biến đổi Fourier thuận - nghịch và thuật toán biến đổi nhanh (Fast Fourier Transform) [57], công thức 12 và 24 tọa độ.
  • Đánh giá mật độ phổ $S(\omega)$ thông qua biến đổi Fourier đối với hàm tương quan $R(\tau)$ (định lý Wiener-Khinchin): $$S(\omega) = \frac{1}{2\pi} \int_{-\infty}^{\infty} R(\tau) e^{-i\omega\tau} d\tau$$
  • Sử dụng ước lượng phổ thực tế có trọng số làm trơn $\lambda(\tau)$ với điểm cắt $T_m$: $$S^(\omega) = \frac{1}{\pi} \int_0^{T_m} \lambda(\tau) R_x^(\tau) \cos \omega\tau d\tau$$ Các hàm làm trơn được khảo sát gồm: Bartlett, Bartlett cải biến, Tukey ($a = 0,25$), Hanning, Parzen, Hamming.
  • Thực nghiệm loại trừ 30 sóng triều đã biết theo phương pháp lọc kiểu Doodson nhằm tránh hiện tượng rò rỉ phổ do đỉnh nhật triều quá mạnh làm lu mờ các tần số khác.

3. Phương pháp mô hình số trị thủy động

Hệ phương trình tuyến tính chuyển động sóng dài trong nước nông tích phân theo độ sâu với xấp xỉ thủy tĩnh trên Trái Đất quay: $$\frac{\partial U}{\partial t} - \lambda V + \frac{1}{\rho_w}\frac{\partial P_a}{\partial x} - \frac{T_x}{\rho_w h} + g \frac{\partial \zeta}{\partial x} + \frac{r U \sqrt{U^2 + V^2}}{h} = 0$$ $$\frac{\partial V}{\partial t} + \lambda U + \frac{1}{\rho_w}\frac{\partial P_a}{\partial y} - \frac{T_y}{\rho_w h} + g \frac{\partial \zeta}{\partial y} + \frac{r V \sqrt{U^2 + V^2}}{h} = 0$$ $$\frac{\partial \zeta}{\partial t} + \frac{\partial (h U)}{\partial x} + \frac{\partial (h V)}{\partial y} = 0$$ Trong đó: $\zeta$ là độ dâng mực nước; $U, V$ là các thành phần vận tốc trung bình theo chiều sâu; $\lambda = 2\omega \sin \phi$ là thông số Coriolis; $P_a$ là áp suất khí quyển; $T_x, T_y$ là ứng suất gió mặt biển; $r$ là hệ số ma sát đáy; $h$ là độ sâu đáy biển; $g$ là gia tốc trọng lực; $\rho_w$ là mật độ nước biển.

Mô hình được rời rạc hóa trên lưới sai phân hiện Richardson (staggered grid). Biên cứng áp dụng điều kiện không xuyên qua ($U \cos \alpha + V \sin \alpha = 0$). Biên lỏng áp dụng mực nước bằng 0, mực nước triều dự tính hoặc điều kiện phát xạ sóng.

Nội dung chính theo từng chương

Chương 1 - Khái quát về điều kiện tự nhiên và đặc điểm dao động mực nước Biển Đông. Cơ sở phương pháp nghiên cứu

Chương 1 tổng hợp các nét đặc trưng về hình thái, khí tượng thủy văn và cơ sở toán học:

  • Đặc điểm hình thái và độ sâu: Biển Đông trải dài từ $99^\circ 30'$ Đ đến $121^\circ$ Đ, $3^\circ$ N đến $25^\circ$ B, diện tích khoảng $3,5 \cdot 10^6 \text{ km}^2$. Địa hình đáy bất đối xứng rõ rệt: bồn trũng sâu cực đại 5560 m nằm lệch về phía đông; thềm lục địa sâu dưới 200 m chiếm hơn một nửa diện tích tập trung ở phía nam, tây nam và hai vịnh lớn (vịnh Bắc Bộ và vịnh Thái Lan có độ sâu từ vài mét đến dưới 100 m). Biển thông với Thái Bình Dương qua eo Basi (sâu trên 3000 m, rộng 350 km), eo Đài Loan (sâu 50 m, rộng 200 km) và thông với biển Giava qua eo Malacca (sâu 30–40 m, rộng 450 km).
  • Chế độ gió: Gió mùa đông bắc (tháng 10 đến tháng 4 năm sau) có cường độ mạnh và độ ổn định cao; gió mùa tây nam (mùa hè) có tính bất đồng nhất không gian lớn.
  • Chế độ mực nước: Tồn tại cả 4 loại thủy triều, trong đó nhật triều chiếm ưu thế tại vùng khơi và hai vịnh (độ lớn triều tại đỉnh vịnh Bắc Bộ đạt tới 6 m); bán nhật triều chỉ phân bố ở một số khu vực hẹp (eo Đài Loan, Thuận An). Nước dâng do bão ven bờ Việt Nam có thể đạt tới 250 cm; dao động mùa do gió mùa đạt biên độ 30–40 cm tại bờ Tây.
  • Hệ thống hóa các phương pháp giải tích sóng đứng trong kênh (Merian, Dubois), thuật toán phân tích Fourier, phương pháp hàm tương quan – mật độ phổ và thiết lập hệ phương trình vi phân sai phân thủy động nước nông.

Chương 2 - Khảo sát dao động tự do của Biển Đông

Chương 2 là nội dung trung tâm của luận án, tập trung giải bài toán dao động riêng:

  • Mô hình số dao động tự do: Giải hệ phương trình sóng dài tuyến tính không ma sát ($r = 0$, $T_x = T_y = 0$, $\nabla P_a = 0$). Điều kiện ban đầu đặt trường vận tốc triệt tiêu ($u = v = 0$) và kích động một mô nước hình dạng hàm sin bậc hai ở trung tâm biển với biên độ $A = 1,5 \text{ m}$.
  • Thông số lưới tính: Bước lưới không gian $\Delta x = \Delta y = 0,5^\circ \approx 55,569 \text{ km}$; bước thời gian $\Delta t = 150 \text{ giây}$ (thỏa mãn tiêu chuẩn ổn định Courant-Friedrichs-Lewy). Thời gian tính 1000 giờ, lưu trữ chuỗi mực nước tại 16 điểm tính đại diện trên toàn biển để phân tích phổ và khôi phục cấu trúc không gian các mốt dao động.
  • Kết quả 5 nhóm chu kỳ dao động tự do:
    1. Nhóm chu kỳ trên một ngày: 33,3 giờ và 55,6 giờ.
    2. Nhóm chu kỳ toàn nhật: 23,8; 24,8 và 25,0 giờ.
    3. Nhóm chu kỳ trung gian: 14,3; 17,2 và 19,2 giờ.
    4. Nhóm chu kỳ bán nhật: 10,6; 11,6 và 13,2 giờ.
    5. Nhóm chu kỳ nước nông: Từ 10,6 giờ trở xuống (9,8; 9,6; 9,4; 9,1; 8,6; 8,2; 7,9; 7,6; 7,1; 6,9; 6,1 và 4,7–4,2 giờ).
Chu kỳ mốt riêng (giờ) Trị số bình phương trung bình toàn biển $\Delta$ (mm)
60,0 16
24,0 10
19,2 64
17,2 14
11,6 18
10,6 15
9,7 8

Chu kỳ 19,2 giờ là mốt dao động riêng chủ đạo toàn biển mang năng lượng lớn nhất, có hệ thống amphidromy xoay trái với điểm vô triều ảo ở đông bắc biển.

  • Giải thích cơ chế cộng hưởng thủy triều dọc bờ biển Việt Nam: Ở phía Bắc và trung tâm biển, chu kỳ tự do nhóm toàn nhật là 25,0 giờ gần với chu kỳ sóng $O_1$ (25,82 giờ) hơn sóng $K_1$ (23,93 giờ), nên sóng $O_1$ cộng hưởng mạnh hơn. Càng xuống phía Nam, chu kỳ tự do giảm xuống 23,8–24,8 giờ, gần như trùng với chu kỳ sóng $K_1$, làm sóng $K_1$ tăng vọt biên độ. Đối với sóng bán nhật, chu kỳ tự do 11,6 giờ ở phía Nam gần với sóng $M_2$ (12,42 giờ) hơn sóng $S_2$ (12,00 giờ), khiến biên độ sóng $M_2$ tăng mạnh.
Trạm quan trắc Biên độ $K_1$ (cm) Biên độ $O_1$ (cm) Tỷ số $K_1 / O_1$ Biên độ $S_2$ (cm) Biên độ $M_2$ (cm) Tỷ số $S_2 / M_2$
Hòn Dấu 65,16 74,71 0,87 5,03 9,29 0,54
Đà Nẵng 19,44 12,93 1,50 5,75 17,23 0,33
Quy Nhơn 30,88 26,49 1,16 6,65 16,06 0,41
Vũng Tàu 59,48 45,22 1,32 28,64 74,83 0,38
Rạch Giá 20,46 11,82 1,73 3,04 16,12 0,18
  • Phân loại 6 vùng cộng hưởng của Biển Đông:
    • Vùng 1 (Trung tâm biển và bờ sâu Trung Nam Bộ): Cộng hưởng chu kỳ ngày và trung gian $\rightarrow$ Nhật triều không đều, biên độ nhỏ.
    • Vùng 2 (Vịnh Bắc Bộ trừ tây nam cửa vịnh): Cộng hưởng trung gian, ngày và nước nông (không cộng hưởng nửa ngày) $\rightarrow$ Nhật triều đều, biên độ lớn.
    • Vùng 3 (Tây nam cửa Vịnh Bắc Bộ và tây bắc Kalimantan): Chỉ cộng hưởng chu kỳ trung gian $\rightarrow$ Bán nhật triều đều, biên độ rất nhỏ.
    • Vùng 4 (Dải hẹp ven bờ Đông Nam Bộ và tây nam Kalimantan): Cộng hưởng cả 5 nhóm, nhóm bán nhật mạnh hơn $\rightarrow$ Thủy triều hỗn hợp, bán nhật triều không đều.
    • Vùng 5 (Bắc vịnh Thái Lan): Cộng hưởng cả 5 nhóm, chu kỳ toàn nhật mạnh nhất $\rightarrow$ Nhật triều, sóng nước nông phong phú, biên độ lớn.
    • Vùng 6 (Ven bờ nam Trung Quốc từ ĐB Hải Nam đến eo Đài Loan): Không cộng hưởng chu kỳ ngày, cộng hưởng bán nhật và nước nông $\rightarrow$ Bán nhật triều kết hợp sóng nước nông.

Chương 3 – Phổ mực nước ở ven bờ Tây Biển Đông

  • Đặt vấn đề phân tích phổ thực nghiệm từ chuỗi số liệu quan trắc từng giờ độ dài 1 năm tại 5 trạm bờ Tây: Hòn Dấu, Đà Nẵng, Quy Nhơn, Vũng Tàu, Rạch Giá.
  • Phân tích khó khăn trong xử lý chuỗi: Hàm tương quan của các yếu tố hải dương suy giảm rất chậm và mang tính chu kỳ; năng lượng sóng nhật triều quá lớn gây hiệu ứng "thấm" làm lu mờ các đỉnh phổ lân cận.
  • Đánh giá thử nghiệm phương pháp lọc triều kiểu Doodson loại trừ 30 sóng điều hòa và lựa chọn hàm cửa sổ làm trơn.
  • Kết quả phổ thực đo khẳng định sự tồn tại của các đỉnh phổ trùng khớp với các chu kỳ dao động tự do đã tính toán ở Chương 2, chứng minh tính thực tế của các mốt dao động riêng trong hệ thống thủy động Biển Đông.

Chương 4 – Tính mực nước trong trường gió mùa

  • Khảo sát biến trình năm và dao động mùa của mực nước biển do gió mùa gây ra (biên độ thực đo đạt 30–40 cm ở dải bờ Tây).
  • Xây dựng mô hình số trị thủy động tuyến tính với trường gió xuất phát gồm hai trường gió mùa đặc trưng: gió mùa đông bắc và gió mùa tây nam.
  • Kết quả tính toán chỉ ra trường mực nước dâng/rút tương ứng với sự phân bố ứng suất gió và địa hình độ sâu: gió mùa đông bắc gây dâng nước mạnh ở thềm lục địa phía nam và tây nam; gió mùa tây nam gây dâng nước ở phía bắc và đông bắc biển.
  • Xác nhận nhận định lý thuyết rằng dao động mùa của Biển Đông là quá trình sóng dài tiến dưới tác dụng cưỡng bức của trường gió mùa, chịu sự biến dạng bởi đường bờ và độ sâu nước nông.

Kết quả và những đóng góp mới

  1. Đóng góp mới về lý luận và phương pháp:

    • Xác định toàn diện phổ dao động tự do của Biển Đông bằng phương pháp số trị trên địa hình thực, làm rõ 5 nhóm chu kỳ dao động riêng và chứng minh mốt 19,2 giờ là dao động chủ đạo toàn biển.
    • Làm sáng tỏ cơ sở động lực học của hiện tượng nhật triều chiếm ưu thế tại Biển Đông: tính chất nhật triều ngự trị là kết quả của sự cộng hưởng giữa sóng triều Thái Bình Dương với các mốt dao động riêng toàn nhật (23,8–25,0 giờ) của thủy vực.
    • Giải thích cơ chế biến đổi tỷ số biên độ $K_1/O_1$ tăng dần từ Bắc vào Nam và sự phát triển của sóng bán nhật $M_2$ ở phía Nam dựa trên sự tương thích tần số với các chu kỳ dao động tự do cục bộ.
    • Hoàn thiện phương pháp phân tích điều hòa mực nước theo sơ đồ bình phương nhỏ nhất (công thức 1.3), cho phép xử lý chính xác chuỗi số liệu dài, chuỗi gián đoạn mà không phụ thuộc vào phép lấy gần đúng thông số thiên văn.
    • Thiết lập sơ đồ phân vùng Biển Đông thành 6 vùng cộng hưởng động lực học có giá trị dự báo cấu trúc thủy triều.

  2. Đóng góp thực tiễn:

    • Cung cấp bức tranh định lượng về trường mực nước dâng do gió mùa (đạt 30–40 cm), hỗ trợ tính toán thiết kế đê điều, luồng tàu, công trình cảng biển và hệ thống thủy lợi ven bờ.
    • Kết quả phân tích phổ và hằng số điều hòa tại các trạm Hòn Dấu, Đà Nẵng, Quy Nhơn, Vũng Tàu, Rạch Giá cung cấp thông số đầu vào tin cậy cho công tác hiệu chỉnh bảng thủy triều và tính toán mực nước cực trị.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế:

    • Mô hình số trị sử dụng hệ phương trình sóng dài nước nông tuyến tính, chưa tính đến các số hạng phi tuyến đầy đủ trong tương tác giữa triều, nước dâng và sóng nước nông.
    • Thiếu nguồn dữ liệu quan trắc khí tượng chi tiết ngoài khơi trong bão để kiểm chứng khả năng cộng hưởng của ngoại lực với mốt dao động riêng 19,2 giờ.
    • Việc lọc tách 30 sóng triều để xử lý phổ thực nghiệm chuỗi mực nước còn gặp khó khăn do tính dao động kéo dài của hàm tương quan.

  • Hướng nghiên cứu tiếp:

    • Phát triển mô hình số trị phi tuyến với lưới tính chi tiết hóa và điều kiện biên thực tế hơn.
    • Nghiên cứu bài toán tương tác tổng cộng giữa thủy triều, nước dâng do gió và sóng thềm lục địa.
    • Hoàn thiện mô hình số trị trường gió phục vụ công tác dự báo nghiệp vụ các dao động mực nước phi triều.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và giảng viên chuyên ngành Hải dương học, Khí tượng Thủy văn biển: Luận án là tài liệu tham khảo về phương pháp tính dao động riêng của thủy vực phức tạp, kỹ thuật phân tích phổ chuỗi thời gian thực tế và thuật toán bình phương nhỏ nhất trong phân tích điều hòa.
  • Đối với kỹ sư công trình thủy, giao thông hàng hải và quy hoạch ven biển: Cung cấp số liệu định lượng về biên độ sóng triều thành phần, dao động mùa, mức độ dâng nước do gió mùa và phân bố năng lượng sóng nước nông tại các vùng ven bờ Việt Nam.

Câu hỏi thường gặp

1. Chu kỳ dao động tự do nào đóng vai trò chủ đạo trên toàn Biển Đông và có đặc điểm cấu trúc ra sao?
Chu kỳ 19,2 giờ là mốt dao động tự do chủ đạo trên toàn Biển Đông, mang phương sai lớn nhất ở hầu hết 16 điểm tính với trị số bình phương trung bình toàn biển đạt $\Delta = 64\text{ mm}$ (cao nhất trong các mốt). Cấu trúc không gian của mốt này tạo thành hệ thống amphidromy quay trái với điểm vô triều ảo nằm ở phía đông bắc biển, kéo theo toàn bộ khối nước trung tâm và các vịnh cùng tham gia dao động.

2. Luận án giải thích thế nào về quy luật tỷ số biên độ $K_1 / O_1$ tăng dần từ Bắc vào Nam dọc bờ biển Việt Nam?
Ở phần phía bắc và trung tâm Biển Đông, chu kỳ dao động tự do nhóm toàn nhật là 25,0 giờ, gần với chu kỳ sóng $O_1$ (25,82 giờ) hơn sóng $K_1$ (23,93 giờ), nên sóng $O_1$ được cộng hưởng mạnh hơn. Càng đi xuống phía Nam, chu kỳ tự do địa phương giảm xuống 23,8–24,8 giờ, trùng khớp với chu kỳ sóng $K_1$, làm biên độ sóng $K_1$ khuếch đại vượt trội. Do đó, tỷ số $K_1/O_1$ tăng dần từ 0,87 tại Hòn Dấu lên 1,16 tại Quy Nhơn, 1,32 tại Vũng Tàu và đạt 1,73 tại Rạch Giá.

3. Tiêu chí phân chia 6 vùng cộng hưởng của Biển Đông trong luận án là gì?
Sự phân vùng được thiết lập thuần túy dựa trên phổ các chu kỳ dao động tự do xuất hiện tại từng khu vực địa hình. Chẳng hạn, Vùng 2 (Vịnh Bắc Bộ trừ tây nam) cộng hưởng với các chu kỳ trung gian, toàn nhật và nước nông nhưng không cộng hưởng chu kỳ bán nhật, dẫn đến chế độ nhật triều đều biên độ lớn; trong khi Vùng 3 (tây nam cửa Vịnh Bắc Bộ) chỉ cộng hưởng chu kỳ trung gian, không cộng hưởng chu kỳ triều, nên có chế độ bán nhật triều đều nhưng biên độ rất nhỏ.

4. Cải tiến phương pháp phân tích điều hòa của tác giả mang lại ưu thế kỹ thuật gì?
Bằng cách chuyển đổi công thức điều hòa sang dạng biến đổi (1.3) với các ẩn số $X_i = H_i \cos G_i$ và $Y_i = H_i \sin G_i$ giải bằng phương pháp bình phương nhỏ nhất, tác giả không phải lấy gần đúng các thông số thiên văn ở giữa kỳ quan trắc, loại bỏ sai số tích lũy khi xử lý chuỗi số liệu dài và cho phép phân tích các chuỗi số liệu không liên tục.

Kết luận

Luận án của tác giả Phạm Văn Huấn (1993) đã giải quyết hệ thống bài toán dao động tự do của Biển Đông bằng mô hình số trị thủy động, chỉ ra 5 nhóm chu kỳ riêng với mốt dao động chủ đạo 19,2 giờ. Kết quả nghiên cứu đã chứng minh bản chất hiện tượng nhật triều chiếm ưu thế tại Biển Đông là do hiệu ứng cộng hưởng động lực học giữa sóng triều đại dương và các mốt dao động riêng của thủy vực. Công trình đóng góp cải tiến phương pháp phân tích điều hòa theo nguyên lý bình phương nhỏ nhất, đồng thời xác lập mô hình tính toán dao động mùa của mực nước dưới tác động của trường gió mùa.