Luận án: Biến đổi vận tốc âm Vịnh Bắc Bộ và hoàn thiện quy trình lấy mẫu thủy âm

Tổng quan về luận án

Nghiên cứu của tác giả Nguyễn Văn Cương (2021) với đề tài "Nghiên cứu quy luật biến đổi vận tốc âm và hoàn thiện quy trình lấy mẫu xác định vận tốc âm trong môi trường nước biển khu vực Vịnh Bắc Bộ Việt Nam phục vụ khai thác hiệu quả các thiết bị thủy âm" thuộc chuyên ngành Kỹ thuật trắc địa - bản đồ (Mã số: 9520503), dưới sự hướng dẫn khoa học của PGS.TS Đặng Nam Chinh tại Trường Đại học Mỏ - Địa chất, đặt nền móng tiên phong cho việc chuẩn hóa âm học thủy đạc tại vùng biển đặc thù Việt Nam. Với bờ biển dài hơn 3.260 km và trung bình cứ 23 km lại có một cửa sông đổ ra biển, môi trường nước biển Việt Nam – đặc biệt là Vịnh Bắc Bộ – có cấu trúc hải văn phân tầng vô cùng phức tạp. Tuy nhiên, trước nghiên cứu này, công tác đo đạc địa hình đáy biển bằng thiết bị thủy âm tại Việt Nam chủ yếu áp dụng quy trình lấy mẫu rời rạc hoặc sử dụng các công thức tham chiếu toàn cầu chưa qua hiệu chuẩn theo đặc thù thủy văn địa phương, dẫn tới sai số khúc xạ âm học nghiêm trọng khi đo sâu.

Khoảng trống học thuật (Research Gap) cốt lõi được xác định rõ ràng: Chưa từng có một công trình nào lượng hóa chi tiết quy luật biến thiên không - thời gian của trường vận tốc âm (Sound Velocity - $v$) trên toàn diện tích Vịnh Bắc Bộ, đồng thời thiếu hụt một cơ sở dữ liệu (CSDL) chuẩn hóa tích hợp thuật toán xử lý để tối ưu hóa tần suất lấy mẫu trực địa. Các nghiên cứu trước đây (như Dương Quốc Lương, 2006; Phạm Văn Thục, 2011) chỉ dừng lại ở khảo sát điểm hoặc nghiên cứu vật lý biển tổng quan, chưa giải quyết bài toán cốt lõi trong đo đạc kỹ thuật: Làm thế nào để giảm thiểu số lần thả đầu đo SVP/CTD ngoài hiện trường mà vẫn kiểm soát sai số độ sâu đạt chuẩn quốc tế IHO S-44 và Tiêu chuẩn ngành của Bộ Tài nguyên và Môi trường (BTNMT)?

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

• Câu hỏi nghiên cứu 1 (RQ1): Các tham số hải văn (nhiệt độ $T$, độ mặn $S$, độ sâu/áp suất $P$) chi phối quy luật biến thiên của vận tốc âm tại Vịnh Bắc Bộ theo không gian 3 chiều và thời gian 12 tháng trong năm như thế nào?
• Câu hỏi nghiên cứu 2 (RQ2): Mức độ suy giảm độ chính xác của các thiết bị đo sâu hồi âm đơn tia (SBES), đa tia (MBES), quét sườn (SSS) và định vị thủy âm dưới tác động của biến thiên trường vận tốc âm tại vùng nước nông Vịnh Bắc Bộ là bao nhiêu?
• Câu hỏi nghiên cứu 3 (RQ3): Liệu việc thiết lập CSDL vận tốc âm chuyên dụng kết hợp mô hình nội suy không - thời gian có đáp ứng được độ chính xác thành lập bản đồ địa hình đáy biển tỷ lệ 1:10.000 và nhỏ hơn theo chuẩn IHO S-44 mà không cần lấy mẫu liên tục?
• Giả thuyết khoa học 1 (H1): Sự phân tầng nhiệt - muối tại Vịnh Bắc Bộ do tác động đan xen của chế độ hoàn lưu gió mùa và dòng xả lục địa tạo ra độ biến thiên vận tốc âm trong năm vượt ngưỡng $\Delta v > 25\text{ m/s}$, tập trung chủ yếu ở tầng nước mặt đến độ sâu 30 m.
• Giả thuyết khoa học 2 (H2): Áp dụng CSDL vận tốc âm chuẩn hóa kết hợp thuật toán truy xuất theo tọa độ $(\varphi, \lambda, z, t)$ cho phép hiệu chỉnh khúc xạ tia âm với sai số $RMSE < 0{,}25\text{ m}$, đáp ứng hoàn toàn tiêu chuẩn đo vẽ bản đồ địa hình đáy biển tỷ lệ $\le 1:10.000$, giúp cắt giảm trên 50% khối lượng lấy mẫu trực tiếp bằng thiết bị đo sâu vận tốc âm (SVP).

Khung lý thuyết của luận án được xây dựng vững chắc trên nền tảng cơ học sóng âm đàn hồi xuất phát từ định luật Isaac Newton (1687, Principia), định luật khúc xạ Snell ($c_1^{-1}\sin\theta_1 = c_2^{-1}\sin\theta_2$), lý thuyết lan truyền âm trong chất lỏng phân tầng của Medwin (1975), Del Grosso (1974), Chen & Millero (1977), Mackenzie (1981) và mô hình chuyển đổi áp suất - độ sâu của Leroy & Parthiot (1998). Phạm vi nghiên cứu bao quát toàn bộ vùng biển Vịnh Bắc Bộ trong tọa độ $17^\circ 45'\text{N} - 21^\circ 30'\text{N}$ và $105^\circ 45'\text{E} - 108^\circ 15'\text{E}$, xử lý chuỗi dữ liệu hải dương học lịch sử nhiều năm kết hợp dữ liệu đo thực nghiệm thực tế bằng thiết bị chuyên dụng SVP-15 và máy đo sâu đa tia đa tần số.

Literature Review và Positioning

Lịch sử nghiên cứu âm học dưới nước bắt đầu từ các quan sát định tính của Aristotle (384–322 TCN), Leonardo da Vinci (1490) cho đến thí nghiệm đo vận tốc âm đầu tiên trong nước tại Hồ Geneva năm 1826 của Jean-Daniel Colladon và Charles Sturm (đo được $v = 1435\text{ m/s}$ ở $8^\circ\text{C}$, sai khác chỉ $3\text{ m/s}$ so với giá trị hiện đại). Sự phát triển của các công thức giải tích xác định vận tốc âm thực nghiệm chứng kiến hai luồng quan điểm và trường phái mô hình hóa đối lập trong tài liệu quốc tế:

1. Trường phái phương trình phức hợp đa biến đại dương sâu: Khởi xướng bởi Wayne D. Wilson (1960), tiếp nối bởi V.A. Del Grosso (1974) với phương trình 19 tham số bậc cao và Chen & Millero (1977) (thuật toán được UNESCO công nhận chính thức). Nhóm nghiên cứu này nhấn mạnh vào tính chính xác tuyệt đối ở áp suất cực cao (độ sâu $0 - 11.000\text{ m}$, áp suất $0 - 1000\text{ bar}$), ưu tiên khảo sát hải dương học toàn cầu và rãnh sâu. Tuy nhiên, điểm hạn chế khi áp dụng vào vùng biển ven bờ nước nông là độ phức tạp tính toán thừa thãi và độ nhạy kém với các biến thiên cục bộ vùng cửa sông có độ mặn thấp ($S < 25\text{ ppt}$).
2. Trường phái phương trình xấp xỉ rút gọn vùng thềm lục địa/nước nông: Tiêu biểu là H. Medwin (1975) với công thức 5 tham số cho độ sâu $D \le 1000\text{ m}$, K.V. Mackenzie (1981) với phương trình 9 số hạng, và A.B. Coppens (1981) cho vùng nước có chiều sâu trung bình ($D \le 4000\text{ m}$). Medwin đề xuất phương trình tối ưu: $$v = 1449{,}2 + 4{,}6T - 0{,}055T^2 + 0{,}00029T^3 + (1{,}34 - 0{,}010T)(S - 35) + 0{,}016D$$ Công thức này chứng minh tính hiệu quả vượt trội trong thực thi thuật toán thời gian thực trên các hệ thống đo sâu thương mại tại vùng nước nông.

                    ┌─────────────────────────────────────────────────────────┐
                    │               LỊCH SỬ PHÁT TRIỂN LÝ THUYẾT              │
                    │               VẬN TỐC ÂM TRONG NƯỚC BIỂN                │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                  ┌──────────────────────────────┴──────────────────────────────┐
                  ▼                                                             ▼
┌───────────────────────────────────┐                         ┌───────────────────────────────────┐
│     TRƯỜNG PHÁI BIỂN SÂU/ÁP CAO   │                         │    TRƯỜNG PHÁI VÙNG NƯỚC NÔNG     │
│  - Wilson (1960)                  │                         │  - Medwin (1975): D ≤ 1000m       │
│  - Del Grosso (1974): 19 tham số  │                         │  - Mackenzie (1981): 9 số hạng    │
│  - Chen & Millero (1977): UNESCO  │                         │  - Coppens (1981): D ≤ 4000m      │
└─────────────────┬─────────────────┘                         └─────────────────┬─────────────────┘
                  │                                                             │
                  └──────────────────────────────┬──────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │       RESEARCH GAP TẠI VỊNH BẮC BỘ & VIỆT NAM           │
                    │  - Chưa có CSDL không-thời gian chuẩn hóa               │
                    │  - Chưa liên kết lý thuyết khúc xạ với tiêu chuẩn IHO   │
                    │  - Thiếu quy trình lấy mẫu SVP tối ưu cho MBES/SSS      │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │              ĐÓNG GÓP CỦA LUẬN ÁN (2021)                │
                    │  - CSDL chuẩn hóa từ WOD & Viễn Đông Nga                │
                    │  - Mô hình biến thiên Δv lên tới 27,5 m/s               │
                    │  - Quy trình lấy mẫu tối ưu cho bản đồ tỷ lệ 1:10.000   │
                    └─────────────────────────────────────────────────────────┘

Tại Việt Nam, các công trình của Phạm Văn Thục (2011) trong chuyên khảo "Âm học biển và trường sóng âm tại khu vực Biển Đông Việt Nam" hay dự án KC.09 (2007) đã phân tích trường âm quy mô lớn trên toàn Biển Đông, nhưng thiếu độ phân giải lưới không gian chi tiết cho Vịnh Bắc Bộ. Luận văn của Vũ Hồng Tập (2011) chỉ dừng lại ở việc khảo sát sai số vận tốc âm trên máy đo sâu đơn tia (SBES), trong khi nghiên cứu của Phạm Vũ Vinh Quang (2013) khảo sát thiết bị quét sườn (SSS) mà chưa đi sâu vào bản chất vật lý phân tầng cột nước.

Luận án của Nguyễn Văn Cương định vị chính xác điểm giao thoa giữa âm học đại dương vật lý và kỹ thuật trắc địa công trình biển. Bằng việc so sánh đối chiếu chuẩn mực với Tiêu chuẩn S-44 của Tổ chức Thủy đạc Quốc tế (IHO Standards for Hydrographic Surveys - Khóa ấn bản thứ 5 và 6), tiêu chuẩn của Cơ quan Thủy đạc Canada (CHS), Cơ quan Thủy đạc New Zealand (LINZ) và Hiệp hội Nhà thầu Hàng hải Quốc tế (IMCA), tác giả đã chứng minh rằng các mô hình toàn cầu (như GDEMV 3.0 của Hải quân Hoa Kỳ) tồn tại độ lệch cục bộ đáng kể tại dải ven bờ Việt Nam. Luận án bù đắp khoảng trống này bằng việc xây dựng bộ dữ liệu tích hợp từ World Ocean Database (WOD), CSDL Viễn Đông Nga và số liệu quan trắc thực địa, thiết lập tiêu chuẩn kỹ thuật lấy mẫu thích ứng hoàn toàn với điều kiện nhiệt đới gió mùa.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và kiểm chứng thực nghiệm giới hạn áp dụng của các mô hình lý thuyết âm học kinh điển trong môi trường nước nông nhiệt đới chịu ảnh hưởng mạnh của chế độ thủy triều và xả lũ cửa sông:

1. Lượng hóa độ nhạy tham số môi trường: Xác định chính xác trọng số ảnh hưởng của các thành phần vi phân đối với gradient vận tốc âm tại Vịnh Bắc Bộ. Trong điều kiện tiêu chuẩn ($T = 0^\circ\text{C}, S = 35\text{ ppt}, D = 0\text{ m}$), vận tốc âm cơ sở là $1449{,}14\text{ m/s}$. Tại Vịnh Bắc Bộ, mức độ ảnh hưởng thực nghiệm được lượng hóa:
   • Khi nhiệt độ $T$ biến thiên $1^\circ\text{C}$, $\Delta v \approx 4{,}6\text{ m/s}$ (tương ứng gradient $\partial v/\partial T$).
   • Khi độ mặn $S$ biến thiên $1\text{ ppt}$, $\Delta v \approx 1{,}34 - 1{,}4\text{ m/s}$ ($\partial v/\partial S$).
   • Khi độ sâu $D$ tăng $100\text{ m}$, $\Delta v \approx 1{,}6\text{ m/s}$ do tác động của áp suất thủy tĩnh ($\partial v/\partial D$).
2. Chứng minh sự dịch chuyển vị thế chi phối giữa nhiệt độ và độ mặn: Luận án phát hiện rằng tại vùng nước nông ven bờ ($D < 20\text{ m}$) vào mùa mưa lũ, gradient độ mặn $\partial v/\partial S$ tại các cửa sông Bạch Đằng, Ba Lạt, Đáy trở thành yếu tố chi phối đột biến cục bộ, lấn át gradient nhiệt độ, tạo ra dị thường khúc xạ âm dạng nêm muối (salt-wedge refraction anomaly).
3. Mô hình hóa đường truyền tia âm bằng phương pháp dò tia (Ray Tracing): Tác giả đã thiết lập hệ phương trình vi phân mô tả quỹ đạo chùm tia âm phát ra từ đầu phát biến (Transducer) qua $n$ lớp nước đồng nhất giả định theo định luật Snell: $$z_i = z_{i-1} + \frac{\cos\theta_{i-1} - \cos\theta_i}{g_i} \cdot v_{i-1}, \quad x_i = x_{i-1} + \frac{\sin\theta_i - \sin\theta_{i-1}}{g_i} \cdot v_{i-1}$$ trong đó $g_i = (v_i - v_{i-1})/(z_i - z_{i-1})$ là gradient vận tốc âm lớp thứ $i$. Điều này cho phép lượng hóa chính xác sai số vị trí điểm đáy biển $(\Delta x, \Delta y, \Delta z)$ khi góc tới $\theta$ của chùm tia ngoài (outer beams) trong hệ thống MBES mở rộng đến $60^\circ - 75^\circ$.

             Transducer (Bộ phát biến)
                    │ \ θ0
                    │  \
     Lớp 1 (v1)     │   \
    ────────────────┴────\────────────────── Mặt phân cách 1 (Refraction)
                          \ θ1
     Lớp 2 (v2 > v1)       \
    ────────────────────────\─────────────── Mặt phân cách 2
                             \
     Lớp 3 (v3)               \ θ2
                               └───► Đáy biển (Seafloor)
    Quỹ đạo thực: Uốn cong liên tục theo định luật Snell (c / sin θ = const)

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên ngành giữa 3 cấu trúc lý thuyết: Lý thuyết lan truyền sóng đàn hồi trong môi trường liên tục, Lý thuyết phân tích sai số trắc địa công trình biển và Lý thuyết tối ưu hóa quá trình thu thập thông tin không gian (Spatial Sampling Optimization).

Cách tiếp cận mới nằm ở giải pháp liên kết trực tiếp giá trị chênh lệch vận tốc âm $\Delta v = |v_{\text{thực tế}} - v_{\text{CSDL}}|$ với Tổng sai số độ sâu cho phép (Total Vertical Uncertainty - TVU) theo chuẩn IHO S-44: $$TVU = \sqrt{a^2 + (b \cdot d)^2}$$ trong đó $a$ là sai số độc lập với độ sâu (hằng số sai số thiết bị), $b$ là hệ số sai số phụ thuộc độ sâu $d$.

Điều kiện biên được xác lập tường minh: Khung phân tích áp dụng tối ưu cho vùng biển thềm lục địa có độ sâu $D \le 120\text{ m}$, biên độ dao động nhiệt độ $T \in [15^\circ\text{C}, 31^\circ\text{C}]$, độ mặn $S \in [18\text{ ppt}, 34{,}5\text{ ppt}]$, không áp dụng cho các vùng biển hở đại dương có độ sâu vượt quá $1000\text{ m}$ hoặc vùng cửa sông nước ngọt hoàn toàn ($S = 0\text{ ppt}$).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án vận dụng thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism) trong thiết kế nghiên cứu thực nghiệm trắc địa biển. Phương pháp nghiên cứu đa tầng (Multi-level Design) được cấu trúc chặt chẽ qua 3 cấp độ:

• Cấp độ vĩ mô (Macro-level): Khảo sát cấu trúc hoàn lưu và trường nhiệt - muối trung bình nhiều năm trên toàn diện tích Vịnh Bắc Bộ dựa trên lưới tọa độ địa lý chuẩn hóa $15' \times 15'$ (khoảng $27 \times 27\text{ km}$).
• Cấp độ trung mô (Meso-level): Thiết lập các mặt cắt âm học kinh tuyến ($105^\circ 45'\text{E}$ đến $108^\circ 15'\text{E}$) và vĩ tuyến ($17^\circ 45'\text{N}$ đến $21^\circ 30'\text{N}$) để bóc tách quy luật biến đổi mặt cắt vận tốc âm (SVP) theo 4 mùa và 12 tháng.
• Cấp độ vi mô (Micro-level): Đánh giá sai số tức thời trên từng tia đo của đầu phát biến MBES và SSS tại các khu vực thực nghiệm trọng điểm (Hải Phòng, Quảng Ninh, vùng biển khơi Vịnh Bắc Bộ).

┌─────────────────────────────────────────────────────────────────────────────┐
│                       THIẾT KẾ NGHIÊN CỨU ĐA TẦNG                           │
├─────────────────────────────────────────────────────────────────────────────┤
│ 1. CẤP VĨ MÔ (Macro): Trường nhiệt-muối toàn Vịnh Bắc Bộ (Lưới 15' x 15')   │
│    Data Sources: World Ocean Database (WOD) + Viễn Đông Nga                 │
├─────────────────────────────────────────────────────────────────────────────┤
│ 2. CẤP TRUNG MÔ (Meso): 28 Mặt cắt SVP kinh/vĩ tuyến theo 12 tháng/4 mùa    │
│    Phân tích: Gradient phân tầng cột nước (Water Column Stratification)     │
├─────────────────────────────────────────────────────────────────────────────┤
│ 3. CẤP VI MÔ (Micro): Kiểm chứng thực địa thiết bị SVP-15 + MBES đa tia     │
│    Xác thực: Tính toán RMSE độ sâu, so sánh chuẩn IHO S-44 & Thông tư 27    │
└─────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được tiến hành với độ tin cậy cơ sở dữ liệu cao:

1. Thu thập dữ liệu nguồn: Tích hợp $1.200+$ trạm đo hải văn lịch sử từ CSDL Hải dương Quốc tế (WOD - NOAA), CSDL Hải dương học Viễn Đông Nga và các trạm đo đạc của Trung tâm Trắc địa và Bản đồ Biển (SEAMAP). Dữ liệu bao gồm chuỗi đo nhiệt độ và độ mặn phân tầng từ bề mặt ($0\text{ m}$) đến đáy ($100\text{ m}$).
2. Giao thức đo đạc thực nghiệm: Sử dụng thiết bị đo biên dạng vận tốc âm chuyên dụng SVP-15 (hãng Valeport/Reson) với độ chính xác đo vận tốc âm $\pm 0{,}05\text{ m/s}$, cảm biến áp suất độ chính xác $0{,}01%$, kết hợp máy đo CTD đo điện dẫn, nhiệt độ, độ sâu chuẩn hóa theo Thông tư 27/2011/TT-BTNMT.
3. Tam giác đạc dữ liệu (Data Triangulation): Kiểm tra chéo giữa 3 nguồn dữ liệu: (1) Giá trị $v$ tính từ công thức thực nghiệm Medwin và Del Grosso qua số liệu CTD; (2) Giá trị $v$ đo trực tiếp bằng cảm biến thời gian truyền âm thực (Time-of-Flight) của thiết bị SVP-15; (3) Giá trị $v$ truy xuất từ mô hình số CSDL xây dựng. Độ tin cậy nhất quán nội tại đạt hệ số tương quan $R^2 > 0{,}98$.

Data và phân tích

Phân tích số liệu sử dụng các công cụ tin học chuyên dụng và thuật toán địa thống kê tiên tiến:

• Chuẩn hóa không gian CSDL: Chuyển đổi dữ liệu hải văn về định dạng bảng cấu trúc quan hệ (RDBMS), phân bố khóa chính theo kinh độ $(\lambda)$, vĩ độ $(\varphi)$, độ sâu $(z)$ và tháng quan trắc $(t)$.
• Phần mềm và thuật toán: Phát triển module phần mềm chuyên dụng trên nền tảng GIS và ngôn ngữ lập trình cấu trúc để tự động trích xuất mặt cắt SVP từ tệp định dạng chuẩn trao đổi quốc tế (như .CSV, định dạng nhị phân Seafloor Information System .SIS, Extended Triton Format .XTF).
• Thống kê sai số và Robustness Checks: Tính toán sai số trung phương (Root Mean Square Error - RMSE) giữa vận tốc âm CSDL và vận tốc âm đo thực tế: $$RMSE_v = \sqrt{\frac{1}{N}\sum_{i=1}^N (v_{\text{CSDL}, i} - v_{\text{Thực tế}, i})^2}$$ Kết quả so sánh trên hàng trăm điểm kiểm tra độc lập cho thấy $RMSE_v$ dao động trong khoảng $0{,}82 - 1{,}45\text{ m/s}$. Sai số này khi đưa vào mô hình khúc xạ tia âm chỉ gây ra sai số độ sâu $\Delta D < 0{,}08\text{ m}$ ở độ sâu $50\text{ m}$ (hoàn toàn nằm trong ngưỡng sai số cho phép $TVU = 0{,}25\text{ m}$ của IHO S-44 Order 1a và Thông tư BTNMT cho bản đồ tỷ lệ 1:10.000).

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Quy luật biến thiên vận tốc âm cực đại theo mùa lên tới $27{,}5\text{ m/s}$: Trích dẫn dữ liệu gốc từ luận án: "Đặc điểm đó đã ảnh hưởng tới sự biến thiên vận tốc âm trong môi trường nước theo thời gian các tháng trong một năm lên tới 27,5 m/s". Vận tốc âm trung bình tại tầng mặt Vịnh Bắc Bộ đạt giá trị thấp nhất vào mùa đông (tháng 1 - tháng 2) với $v \approx 1502 - 1508\text{ m/s}$ do nhiệt độ nước hạ thấp ($18 - 20^\circ\text{C}$), và đạt cực đại vào mùa hè (tháng 7 - tháng 8) lên tới $v \approx 1530 - 1535{,}5\text{ m/s}$ khi nhiệt độ bề mặt vượt $30^\circ\text{C}$.
2. Hiệu ứng phân tầng âm học ba lớp rõ rệt trong mùa hè: Vào mùa hè (tháng 5 đến tháng 9), cột nước Vịnh Bắc Bộ phân hóa thành 3 tầng đặc trưng:
   • Tầng đồng nhất mặt (Surface mixed layer, $0 - 10\text{ m}$): Vận tốc âm cao ($> 1530\text{ m/s}$), gradient thẳng đứng gần như bằng 0.
   • Tầng bước nhảy âm/nhiệt (Thermocline/Sonic layer, $10 - 35\text{ m}$): Vận tốc âm sụt giảm đột ngột với gradient âm rất lớn ($\partial v/\partial z < -0{,}8\text{ s}^{-1}$), làm bẻ cong mạnh các tia âm chiếu xiên của máy đo đa tia MBES về phía đáy.
   • Tầng đáy ($> 35\text{ m}$ đến sát đáy): Vận tốc âm ổn định trở lại ($1515 - 1520\text{ m/s}$).
3. Tính ổn định âm học vào mùa đông: Trong các tháng mùa khô (tháng 11 đến tháng 3 năm sau), do ảnh hưởng của gió mùa Đông Bắc gây xáo trộn đối lưu mạnh, cột nước gần như đồng nhất từ mặt đến đáy, gradient vận tốc âm rất bé ($\Delta v < 2\text{ m/s}$ trên toàn bộ cột nước), tia âm truyền gần như thẳng.
4. Phát hiện dị thường âm học vùng cửa sông ven bờ: Tại các khu vực chịu ảnh hưởng trực tiếp của hệ thống sông Hồng và sông Thái Bình, vào mùa lũ xuất hiện hiện tượng "nghịch nhiệt và nghịch mặn" cục bộ, tạo ra mặt cắt vận tốc âm dạng dị thường (vận tốc âm tầng mặt nhỏ hơn tầng giữa), gây biến dạng dữ liệu quét sườn Side Scan Sonar nếu chỉ dùng giá trị vận tốc âm bề mặt đơn lẻ.

Độ sâu (m) 0 ──────────────────────────────────────────  1535 m/s (Tầng đồng nhất mặt)
            │                                         │
           10 ─────────────────────────────────────────  Thermocline (Bước nhảy âm)
            │    \ Gradient âm rất mạnh               │  Tia âm MBES bị bẻ cong
            │     \ (∂v/∂z < -0.8 s⁻¹)                │  xuống dưới đáy biển
           35 ─────\───────────────────────────────────
            │       │                                 │
           50 ──────┴──────────────────────────────────  1518 m/s (Tầng nước đáy ổn định)
              1510   1520   1530   1540   v (m/s)
                    MẶT CẮT SVP ĐẶC TRƯNG MÙA HÈ VỊNH BẮC BỘ

Implications đa chiều

• Về mặt lý thuyết: Bổ sung vào kho tàng âm học biển thế giới bức tranh chi tiết về động lực học trường âm vùng vịnh nhiệt đới nửa kín; cung cấp bộ tham số thực nghiệm chính xác để hiệu chỉnh thuật toán UNESCO/Medwin trong vùng biển thềm lục địa Việt Nam.
• Về mặt phương pháp luận: Thiết lập quy trình công nghệ 2 bước đột phá: Thay vì bắt buộc phải đo SVP trực tiếp trước mỗi ca đo (gây lãng phí thời gian neo tàu), kỹ sư thủy đạc có thể trích xuất trực tiếp profile từ CSDL vận tốc âm của luận án để nạp vào phần mềm thu thập số liệu (như QPS QINSy, HYPACK, Teledyne PDS, Kongsberg SIS). Việc lấy mẫu trực tiếp chỉ thực hiện 1 lần duy nhất khi có cảnh báo dị thường thời tiết.
• Về mặt ứng dụng thực tiễn: Tối ưu hóa hiệu suất vận hành của các tàu khảo sát đo đạc biển; giảm thời gian dừng chờ đo hải văn từ 20-30% tổng thời gian hành trình khảo sát, tiết kiệm hàng tỷ đồng chi phí ca máy và nhiên liệu mỗi năm cho các đơn vị đo đạc trắc địa biển.
• Về chính sách và tiêu chuẩn hóa: Cung cấp cơ sở khoa học trực tiếp để Tổng cục Biển và Hải đảo Việt Nam, Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam soát xét, cập nhật Thông tư quy định kỹ thuật đo vẽ địa hình đáy biển bằng công nghệ đo sâu đa tia.

Limitations và Future Research

Mặc dù đạt được những kết quả có tính đột phá cao, luận án vẫn thẳng thắn thừa nhận 3 hạn chế khoa học do điều kiện khách quan:

1. Mật độ trạm đo chưa đồng đều ở vùng biển giáp ranh: Dữ liệu lịch sử tại khu vực đường phân định Vịnh Bắc Bộ và vùng biển xa bờ còn thưa hơn so với dải ven bờ Hải Phòng - Quảng Ninh, dẫn tới độ tin cậy nội suy không gian ở một số ô lưới $15' \times 15'$ phía ngoài vịnh chưa tuyệt đối đồng nhất.
2. Độ phân giải thời gian giới hạn theo bước tháng: CSDL hiện tại được chuẩn hóa theo giá trị trung bình 12 tháng và 4 mùa, chưa mô hình hóa được biến thiên vận tốc âm tức thời theo chu kỳ nhật triều/bán nhật triều (chu kỳ 12h và 24h) tại các luồng hàng hải ven bờ hẹp.
3. Thiếu vắng dữ liệu trường âm cực đoan: Chưa tích hợp các kịch bản dị thường âm học trong điều kiện thời tiết cực đoan như bão biển cấp $11 - 12$ hoặc hiện tượng nước dâng do bão khi các trạm đo không thể vận hành.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng nghiên cứu chuyên sâu:

• Hướng 1: Tích hợp trí tuệ nhân tạo (AI/Machine Learning) và mạng nơ-ron tích chập (CNN/LSTM) để dự báo mặt cắt vận tốc âm thời gian thực 4D $(\varphi, \lambda, z, t)$ dựa trên dữ liệu viễn thám nhiệt độ mặt biển (SST) từ vệ tinh kết hợp mô hình hoàn lưu đại dương ROMS/FVCOM.
• Hướng 2: Mở rộng xây dựng CSDL vận tốc âm chi tiết cho toàn bộ các vùng biển trọng điểm khác của Việt Nam (vùng biển Miền Trung, Quần đảo Hoàng Sa, Trường Sa, và vùng biển Tây Nam Bộ).
• Hướng 3: Nghiên cứu ảnh hưởng của trường âm học phân tầng đến công nghệ định vị thủy âm tầm xa (Long Baseline - LBL, Ultra-short Baseline - USBL) phục vụ lắp đặt công trình dầu khí và robot ngầm tự hành (AUV/ROV).
• Hướng 4: Phát triển thiết bị đo sâu hồi âm thông minh tự động hiệu chỉnh khúc xạ thời gian thực tích hợp chip xử lý nhúng dữ liệu CSDL âm học.

Tác động và ảnh hưởng

• Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo nền tảng cho các viện nghiên cứu biển, trường đại học đào tạo chuyên ngành Trắc địa - Bản đồ, Địa vật lý biển, Kỹ thuật Hàng hải. Ước tính công trình tạo tiền đề trích dẫn cho hàng chục bài báo khoa học trên các tạp chí Scopus/ISI chuyên ngành âm học và đo đạc biển (như Applied Acoustics, IEEE Journal of Oceanic Engineering, Continental Shelf Research).
• Chuyển đổi ngành công nghiệp đo đạc biển: Chuẩn hóa quy trình đo đạc thủy đạc cho các đơn vị khảo sát như Tổng công ty Bảo đảm An toàn Hàng hải miền Bắc/miền Nam, Viện Hải dương học, các doanh nghiệp khảo sát công trình dầu khí (PVEP, PTSC), giúp nâng cao độ chính xác đo quét địa hình đáy biển phục vụ nạo vét luồng hàng hải, thiết kế cáp ngầm và quy hoạch điện gió ngoài khơi.
• Tác động quản lý nhà nước và an ninh quốc phòng: Cung cấp cơ sở dữ liệu nền tảng phục vụ phân định ranh giới biển, quản lý tổng hợp tài nguyên môi trường biển và hải đảo. Đặc biệt, dữ liệu phân tầng âm học là thông tin chiến lược tối mật phục vụ công tác tác chiến thủy âm, giám sát tàu ngầm và bố trí hệ thống cảm biến âm thanh cảnh giới ngầm bảo vệ chủ quyền biển đảo Tổ quốc.

                    ┌─────────────────────────────────────────────────────────┐
                    │               HỆ SINH THÁI TÁC ĐỘNG                     │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
         ┌───────────────────────┬───────────────┴───────────────┬───────────────────────┐
         ▼                       ▼                               ▼                       ▼
┌──────────────────┐    ┌──────────────────┐           ┌──────────────────┐    ┌──────────────────┐
│  HỌC THUẬT &     │    │  CÔNG NGHIỆP &   │           │ QUẢN LÝ NHÀ NƯỚC │    │ QUỐC PHÒNG &     │
│   NGHIÊN CỨU     │    │    KHẢO SÁT      │           │  & QUY HOẠCH     │    │   AN NINH        │
│ - Chuẩn hóa âm   │    │ - Cắt giảm 50%   │           │ - Cập nhật Thông │    │ - Tác chiến thủy │
│   học nhiệt đới  │    │   thời gian SVP  │           │   tư kỹ thuật    │    │   âm tàu ngầm    │
│ - Đào tạo kỹ sư  │    │ - Nâng cao độ tin│           │ - Bản đồ chuẩn   │    │ - Bố trí sonar   │
│   thủy đạc cao   │    │   cậy đo quét    │           │   IHO S-44       │    │   cảnh giới ngầm │
└──────────────────┘    └──────────────────┘           └──────────────────┘    └──────────────────┘

Đối tượng hưởng lợi

• Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận mẫu mực về xử lý dữ liệu trắc địa biển, các công thức giải tích âm học và phương pháp thiết lập CSDL chuyên ngành.
• Các nhà khoa học âm học đại dương: Khai thác bộ số liệu thực nghiệm chuẩn hóa để kiểm chứng các mô hình truyền sóng âm phức tạp trong môi trường nước nông.
• Kỹ sư khảo sát thủy đạc và các công ty đo đạc biển: Sở hữu quy trình lấy mẫu tinh gọn, công cụ trích xuất vận tốc âm tin cậy để tối ưu hóa tiến độ thi công các dự án đo vẽ bản đồ đáy biển tỷ lệ lớn.
• Cơ quan quản lý nhà nước (Bộ TN&MT, Bộ GTVT, Bộ Quốc phòng): Có luận cứ khoa học chính xác để ban hành các quy chuẩn kỹ thuật quốc gia, quy hoạch luồng tàu và hoạch định chiến lược kinh tế - quốc phòng biển đảo.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
   Trả lời: Luận án đã mở rộng mô hình truyền sóng đàn hồi trong môi trường liên tục của Medwin (1975) và Del Grosso (1974) vào không gian nước nông nhiệt đới Vịnh Bắc Bộ, xác lập mối tương quan định lượng giữa gradient nhiệt - muối phức hợp với sự uốn cong tia âm trong kỹ thuật đo sâu đa tia MBES, chứng minh rằng biến thiên độ mặn ven bờ có thể tạo ra nghịch chuyển vận tốc âm cục bộ lấn át tác động của nhiệt độ.
2. Đột phá phương pháp luận so với các nghiên cứu trước đây thể hiện thế nào?
   Trả lời: So với nghiên cứu của Dương Quốc Lương (2006) chỉ đo điểm đơn lẻ và Vũ Hồng Tập (2011) chỉ xét máy đơn tia SBES, luận án đã tích hợp thành công CSDL hải dương toàn cầu (WOD) với dữ liệu Viễn Đông Nga thành CSDL 4D chuẩn hóa lưới $15' \times 15'$, liên kết trực tiếp giá trị $\Delta v$ với chỉ số Total Vertical Uncertainty (TVU) theo chuẩn quốc tế IHO S-44 Khóa ấn bản mới nhất.
3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm?
   Trả lời: Đó là biên độ dao động vận tốc âm theo mùa tại Vịnh Bắc Bộ đạt tới $27{,}5\text{ m/s}$ (từ $1508\text{ m/s}$ lên $1535{,}5\text{ m/s}$), vượt xa các dự báo lý thuyết đại dương thông thường ($\approx 10 - 15\text{ m/s}$). Sự biến thiên dữ dội này khẳng định việc áp dụng một giá trị vận tốc âm cố định ($1500\text{ m/s}$) cho công tác đo đạc tại Vịnh Bắc Bộ sẽ gây ra sai số thô bạo về độ sâu, hoàn toàn vi phạm chuẩn đo vẽ trắc địa.
4. Luận án có cung cấp giao thức tái lặp nghiên cứu (Replication Protocol) không?
   Trả lời: Có. Luận án mô tả chi tiết cấu trúc bảng CSDL, thuật toán trích xuất dữ liệu, sơ đồ khối phần mềm xử lý, bảng hệ số hiệu chỉnh Del Grosso/Medwin, và quy trình 4 bước kiểm chuẩn thiết bị SVP-15 theo Thông tư 27/2011/TT-BTNMT, cho phép bất kỳ nhóm nghiên cứu độc lập nào cũng có thể tái lập và mở rộng quy trình.
5. Chương trình nghiên cứu 10 năm được phác thảo ra sao?
   Trả lời: Tác giả định hướng lộ trình 10 năm phát triển hệ thống CSDL âm học 4D thời gian thực cho toàn bộ vùng đặc quyền kinh tế (EEZ) của Việt Nam, tích hợp công nghệ trí tuệ nhân tạo (AI/Machine Learning) dự báo trường vận tốc âm từ viễn thám và thiết lập mạng lưới phao quan trắc âm học tự động (Acoustic Ocean Buoy Network).

Kết luận

Luận án tiến sĩ của NCS Nguyễn Văn Cương đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, tạo nên dấu ấn học thuật quan trọng với 6 đóng góp cốt lõi:

1. Xác lập toàn diện quy luật biến thiên không gian 3 chiều và thời gian 12 tháng của trường vận tốc âm tại Vịnh Bắc Bộ với biên độ dao động năm đạt ngưỡng kỷ lục $27{,}5\text{ m/s}$.
2. Xây dựng hoàn chỉnh Cơ sở dữ liệu vận tốc âm số hóa chuẩn hóa lưới $15' \times 15'$ từ nguồn dữ liệu tích hợp quốc tế và thực địa, phục vụ trực tiếp công tác biên tập trắc địa biển.
3. Lượng hóa chính xác tác động của sai số vận tốc âm đến độ suy giảm vị trí không gian của các hệ thống thủy âm hiện đại (MBES, SBES, SSS, Định vị thủy âm).
4. Đề xuất quy trình công nghệ lấy mẫu vận tốc âm cải tiến, cho phép khai thác CSDL để thay thế phần lớn các ca thả máy trực tiếp ngoài hiện trường.
5. Chứng minh thực nghiệm CSDL và quy trình mới đảm bảo độ chính xác tuyệt đối theo chuẩn IHO S-44 cho thành lập bản đồ địa hình đáy biển tỷ lệ 1:10.000 và nhỏ hơn.
6. Cung cấp luận cứ kỹ thuật chuẩn mực phục vụ xây dựng văn bản quy phạm pháp luật và chiến lược quản lý tài nguyên, bảo vệ chủ quyền biển đảo quốc gia.