Luận án tiến sĩ: Vận tốc âm nước biển Vịnh Bắc Bộ, hoàn thiện quy trình lấy mẫu

Tổng quan về luận án

Công tác trắc địa - bản đồ biển và khảo sát công trình ngầm dưới nước tại Việt Nam đóng vai trò chiến lược đối với an ninh quốc gia, phát triển kinh tế biển và bảo vệ chủ quyền lãnh thổ. Với đường bờ biển trải dài 3.260 km từ Móng Cái đến Hà Tiên, trung bình 23 km có một cửa sông đổ ra biển, môi trường thủy văn nước biển Việt Nam mang tính biến động cực kỳ phức tạp. Trong các kỹ thuật đo đạc hiện đại sử dụng công nghệ thủy âm như đo sâu hồi âm đơn tia (Singlebeam Echosounder - SBES), đo sâu hồi âm đa tia (Multibeam Echosounder - MBES), thủy âm quét sườn (Side Scan Sonar - SSS) và định vị thủy âm dưới nước (Acoustic Positioning), vận tốc lan truyền sóng âm trong cột nước (Sound Velocity Profile - SVP) là tham số then chốt chi phối trực tiếp độ chính xác của mô hình độ sâu và định vị không gian đáy biển.

Khoảng trống nghiên cứu (Research Gap) lớn nhất trong nước trước đây là sự thiếu vắng một công trình phân tích có hệ thống, toàn diện về quy luật biến đổi vận tốc âm 4 chiều (không gian 3D và chu kỳ thời gian 12 tháng) tại các vùng biển nông đặc thù như Vịnh Bắc Bộ. Các khảo sát trước đây chỉ dựa vào việc quan trắc cục bộ tại thời điểm đo hoặc tra cứu các mô hình hải dương học toàn cầu độ phân giải thô, chưa đánh giá định lượng được sự suy giảm độ chính xác của chùm tia âm do biến thiên nhiệt - mặn cửa sông gây ra. Hơn nữa, quy trình lấy mẫu vận tốc âm theo Thông tư 27/2011/TT-BTNMT còn bộc lộ nhiều bất cập về tần suất và mật độ đo, dẫn đến lãng phí thời gian khảo sát thực địa hoặc không đảm bảo chuẩn sai số kỹ thuật.

Luận án tiến sĩ kỹ thuật trắc địa - bản đồ của nghiên cứu sinh Nguyễn Văn Cương 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" (người hướng dẫn khoa học: PGS.TS Đặng Nam Chinh, Trường Đại học Mỏ - Địa chất Hà Nội, 2021, mã ngành: 9520503) đã giải quyết trọn vẹn khoảng trống khoa học này.

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 chặt chẽ:

• Câu hỏi nghiên cứu 1 (RQ1): Quy luật biến đổi vận tốc âm theo kinh độ, vĩ độ, độ sâu cột nước và chu kỳ mùa tại Vịnh Bắc Bộ diễn ra như thế nào dưới tác động của hoàn lưu gió mùa và dòng chảy sông ngòi?
• 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ị thủy âm (SBES, MBES, SSS, USBL) do sai số vận tốc âm gây ra cụ thể là bao nhiêu ở các dải độ sâu khác nhau?
• Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng một Cơ sở dữ liệu (CSDL) vận tốc âm chuyên dụng và hoàn thiện quy trình lấy mẫu SVP tối ưu, vừa giảm thiểu chi phí đo đạc thực địa vừa đảm bảo chuẩn độ chính xác bản đồ địa hình đáy biển tỷ lệ 1:10.000?

Các giả thuyết nghiên cứu bao gồm:

• Giả thuyết H1: Vận tốc âm tại Vịnh Bắc Bộ biến thiên mạnh nhất theo chu kỳ tháng và phân tầng theo mùa, với biên độ biến thiên trong năm vượt quá 25 m/s do biến động nhiệt độ bề mặt và độ mặn cửa sông.
• Giả thuyết H2: Việc tích hợp CSDL vận tốc âm chuẩn hóa kết hợp thuật toán hiệu chỉnh khúc xạ tia âm cho phép thay thế một phần việc đo SVP trực tiếp ngoài thực địa mà vẫn đáp ứng tiêu chuẩn sai số của Tổ chức Thủy đạc Quốc tế (IHO S-44) và quy chuẩn Việt Nam cho bản đồ tỷ lệ 1:10.000 và nhỏ hơn.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng cơ học sóng âm cổ điển từ Isaac Newton (1687), định luật khúc xạ ánh sáng và sóng âm Snell, lý thuyết nhiệt động học chất lỏng biển của Chen & Millero (1977), Medwin (1975), Mackenzie (1981) và hệ tiêu chuẩn độ tin cậy trắc địa thủy đạc IHO S-44.

Đóng góp đột phá của luận án là xác lập bộ bản đồ trường vận tốc âm 12 tháng tại Vịnh Bắc Bộ, chứng minh biên độ dao động vận tốc âm trong năm đạt tới $27.5\text{ m/s}$, đồng thời xây dựng thành công phần mềm khai thác CSDL vận tốc âm 3D/4D từ nguồn World Ocean Database (WOD) và Viễn Đông Nga, chứng minh định lượng sai số độ sâu sau hiệu chỉnh đáp ứng hoàn toàn tiêu chuẩn thành lập bản đồ địa hình đáy biển tỷ lệ 1:10.000 (với sai số trung phương $m_D < 0.35\text{ m}$).

Phạm vi nghiên cứu bao quát toàn bộ vùng nước Vịnh Bắc Bộ từ vĩ tuyến $17^\circ45'\text{N}$ đến $21^\circ30'\text{N}$ và từ kinh tuyến $105^\circ45'\text{E}$ đến $108^\circ15'\text{E}$, độ sâu từ $0\text{ m}$ đến $100\text{ m}$, phân tích chuỗi dữ liệu hải văn nhiều năm kết hợp các chiến dịch đo đạc thực nghiệm bằng thiết bị chuyên dụng SVP-15 của Trung tâm Trắc địa và Bản đồ Biển (SEAMAP).

Literature Review và Positioning

Nghiên cứu về lan truyền sóng âm trong môi trường nước biển đã trải qua hơn ba thế kỷ phát triển. Từ nhận định ban đầu của Aristotle (384-322 TCN), Leonardo da Vinci (1490) đến thực nghiệm của Jean-Daniel Colladon và Charles Sturm (1826) tại hồ Geneva đo được vận tốc âm $1435\text{ m/s}$ ở $8^\circ\text{C}$ (chỉ sai khác $3\text{ m/s}$ so với chuẩn hiện đại). Các công trình đột phá của Paul Langevin (1915) và Robert W. Lichte đã đặt nền móng cho kỹ thuật thủy âm hiện đại.

Dòng nghiên cứu lý thuyết tính toán vận tốc âm thực nghiệm trên thế giới chia thành hai trường phái chính:

1. Trường phái công thức đại dương sâu và áp suất lớn: Wayne D. Wilson (1960) đề xuất phương trình thực nghiệm đầu tiên phụ thuộc vào nhiệt độ ($T$), độ mặn ($S$), áp suất ($P$). Năm 1974, V.A. Del Grosso công bố công thức chính xác cao cho vùng nước tự nhiên từ $0 - 11000\text{ m}$ với hệ số đa thức phức tạp. Năm 1977, Chen và Millero phát triển thuật toán UNESCO chuẩn hóa cho hải dương học quốc tế, tối ưu cho các vực sâu đại dương nơi gradient áp suất đóng vai trò quyết định.
2. Trường phái công thức xấp xỉ vùng nước nông và độ sâu trung bình: Herman Medwin (1975) đề xuất phương trình rút gọn cho độ sâu $D \le 1000\text{ m}$, nhiệt độ $0 - 35^\circ\text{C}$, độ mặn $0 - 45\text{ ppt}$: $$v = 1449.2 + 4.6T - 0.055T^2 + 0.00029T^3 + (1.34 - 0.010T)(S - 35) + 0.016D$$ Tiếp đó, Kenneth V. Mackenzie (1981) đưa ra phương trình 9 số hạng cho độ sâu $0 - 8000\text{ m}$, và A.B. Coppens (1981) mở rộng cho dải độ sâu $0 - 4\text{ km}$. Năm 1998, C. Leroy và F. Parthiot hoàn thiện phương trình quan hệ áp suất - độ sâu chuẩn xác hóa vĩ độ.

Trong nước, các công trình của NGƯT. Nguyễn Minh Huấn (Vật lý biển), PGS. Phạm Văn Huấn (2003, Cơ sở âm học đại dương), GS. Phạm Văn Thục (2011, Âm học biển và trường sóng âm tại khu vực Biển Đông Việt Nam), ThS. Dương Quốc Lương (2006), ThS. Vũ Hồng Tập (2011) đã tiếp cận lý thuyết trường âm và khảo sát đơn tia. Tuy nhiên, chưa có công trình nào giải quyết bài toán mô hình hóa chi tiết biến thiên vận tốc âm theo lưới tọa độ không gian và thời gian thực tế tại Vịnh Bắc Bộ nhằm tối ưu hóa khai thác thiết bị thủy âm đa tia thế hệ mới.

                      TIẾN TRÌNH NGHIÊN CỨU VẬN TỐC ÂM VÀ THỦY ÂM
                      
[Thực nghiệm cổ điển]           [Công thức thực nghiệm]          [Mô hình & Quy chuẩn quốc tế]
Colladon & Sturm (1826) ------> Wilson (1960), Del Grosso (1974) -> IHO S-44, CHS, LINZ, IMCA
Langevin, Lichte (1915)         Medwin (1975), Chen-Millero (1977)  WOD (NOAA), GDEMV 3.0
                                Mackenzie (1981), Coppens (1981)
                                               │
                                               ▼
                      [KHOẢNG TRỐNG TẠI VÙNG NƯỚC NÔNG VIỆT NAM]
                      Thiếu mô hình phân tầng 4D chi tiết Vịnh Bắc Bộ
                      Thiếu CSDL vận tốc âm phục vụ bù sai số khúc xạ MBES/SSS
                                               │
                                               ▼
                      [ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN (2021)]
                      Xác lập quy luật biến đổi không gian - thời gian (Δv = 27.5 m/s)
                      Xây dựng CSDL 3D/4D & Hoàn thiện quy trình lấy mẫu chuẩn 1:10.000

So sánh với hai nghiên cứu quốc tế điển hình:

• Nghiên cứu của Cơ quan Thủy đạc Canada (CHS) và Tổ chức Thủy đạc Quốc tế (IHO S-44 Special Order) yêu cầu đo SVP liên tục $2 - 4\text{ giờ/lần}$ tại các vùng nước có dị thường nhiệt độ; tuy nhiên quy định này gây tốn kém chi phí cực lớn khi áp dụng cho các vùng biển nhiệt đới gió mùa rộng lớn.
• Mô hình số môi trường tổng quát GDEMV 3.0 của Hải quân Hoa Kỳ và CSDL WOD của NOAA có độ phân giải lưới không gian $0.25^\circ \times 0.25^\circ$ hoặc $1^\circ \times 1^\circ$, quá thô để nắm bắt được các hiệu ứng pha trộn nước ngọt cửa sông ven bờ Bắc Bộ.

Luận án đã định vị chính xác điểm giao thoa: vận dụng công thức Medwin trong giới hạn độ sâu $<100\text{ m}$, khai thác nguồn dữ liệu WOD và Viễn Đông Nga đã được kiểm chuẩn bằng số liệu đo thực tế SVP-15, từ đó giải quyết triệt để bài toán kỹ thuật đo đạc thủy đạc tại Việt Nam.

Đó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 ứng dụng lý thuyết khúc xạ tia sóng âm và nhiệt động học biển trong môi trường nước nông nhiệt đới thông qua các đóng góp cụ thể:

1. Mở rộng định luật Snell và lý thuyết vết tia âm (Acoustic Ray Tracing): Chứng minh rằng trong môi trường nước nông ($D < 100\text{ m}$) chịu tác động của hoàn lưu gió mùa nhiệt đới, gradient nhiệt độ thẳng đứng ($\partial T/\partial z$) là nhân tố chi phối chính tới độ cong của tia sóng âm ở các góc quét ngoài của thiết bị MBES, trong khi biến thiên độ mặn ($\partial S/\partial x$) tại các vùng cửa sông là nhân tố gây sai lệch cục bộ theo phương ngang.
2. Khẳng định tính tương thích của phương trình Medwin (1975): Đánh giá độ nhạy cho thấy với sai số đo nhiệt độ $m_T = 0.1^\circ\text{C}$, sai số vận tốc âm là $0.46\text{ m/s}$; sai số độ mặn $m_S = 0.1\text{ ppt}$ gây sai số $0.14\text{ m/s}$; sai số độ sâu $m_D = 1\text{ m}$ chỉ gây sai số $0.016\text{ m/s}$. Từ đó xác lập cơ sở lý thuyết chứng minh việc kiểm chuẩn thiết bị theo Thông tư 27/2011/TT-BTNMT dựa trên công thức Medwin là hoàn toàn chuẩn xác.
3. Mô hình lý thuyết chuyển đổi trạng thái trường âm (Paradigm Shift): Chuyển dịch mô thức từ "đo đạc điểm rời rạc - hiệu chỉnh tĩnh" sang "mô hình hóa trường vận tốc âm 4D - hiệu chỉnh động", cho phép dự báo trường khúc xạ tia âm theo tọa độ $(x, y, z, t)$.

Hệ thống mệnh đề khoa học được xác lập:

• Mệnh đề P1: Tại Vịnh Bắc Bộ, vận tốc âm trung bình cột nước tuân theo quy luật phân kỳ hai mùa rõ rệt: mùa hè hình thành tầng nhảy nhiệt (thermocline) sắc nét làm vận tốc âm giảm đột ngột theo độ sâu; mùa đông hiện tượng đối lưu do gió mùa Đông Bắc làm đồng nhất cột nước, đồ thị SVP tiệm cận đường thẳng đứng.
• Mệnh đề P2: Sai số độ sâu đo bằng MBES do không hiệu chỉnh hoặc hiệu chỉnh sai vận tốc âm tỷ lệ thuận với bình phương góc tới của tia âm ($\sec^2\theta$) và độ sâu cột nước ($D$).

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 ba cấu trúc lý thuyết:

• Lý thuyết vật lý thủy âm: Phương trình truyền sóng, phản xạ, tán xạ và suy giảm năng lượng âm ($EE = SL - 2TL - (NL - DI) + BS - DT$).
• Lý thuyết trắc địa biển và xử lý sai số: Phân tích lan truyền sai số ngẫu nhiên và sai số hệ thống trong phép đo thời gian truyền sóng ($t$) và góc nghiêng chùm tia ($\theta$).
• Khoa học dữ liệu không gian địa lý (GIS): Kỹ thuật nội suy không gian 3D/4D chuẩn hóa CSDL hải dương học.

┌────────────────────────────────────────────────────────────────────────────┐
│                    KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN                    │
└─────────────────────────────────────┬──────────────────────────────────────┘
                                      │
        ┌─────────────────────────────┼─────────────────────────────┐
        ▼                             ▼                             ▼
┌──────────────────┐        ┌──────────────────┐        ┌──────────────────┐
│   HẢI VĂN HỌC    │        │  VẬT LÝ THỦY ÂM  │        │  TRẮC ĐỊA BIỂN   │
│ Nhiệt độ (T)     │        │ Định luật Snell  │        │ Chuẩn IHO S-44   │
│ Độ mặn (S)       │───────>│ Công thức Medwin │───────>│ TVU / THU        │
│ Độ sâu / Áp suất │        │ Vết tia âm (Ray) │        │ Sai số SBES/MBES │
└──────────────────┘        └──────────────────┘        └──────────────────┘
        │                             │                             │
        └─────────────────────────────┼─────────────────────────────┘
                                      │
                                      ▼
             ┌─────────────────────────────────────────────────┐
             │       CƠ SỞ DỮ LIỆU VẬN TỐC ÂM 4D VỊNH BẮC BỘ   │
             │       (12 Mặt cắt vĩ tuyến x 11 Mặt cắt KT)     │
             └────────────────────────┬────────────────────────┘
                                      │
                                      ▼
             ┌─────────────────────────────────────────────────┐
             │   QUY TRÌNH LẤY MẪU & BÙ KHÚC XẠ TỐI ƯU HÓA     │
             │     (Bảo đảm chuẩn bản đồ địa hình 1:10.000)    │
             └─────────────────────────────────────────────────┘

Điều kiện biên của khung phân tích: Vùng biển nông Vịnh Bắc Bộ (độ sâu $0 - 100\text{ m}$), nhiệt độ nước biển $18 - 32^\circ\text{C}$, độ mặn $15 - 34\text{ ppt}$, tần số công tác của thiết bị thủy âm $50\text{ kHz} - 450\text{ kHz}$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism) với thiết kế nghiên cứu hỗn hợp kết hợp giữa mô hình hóa dữ liệu thứ cấp quy mô lớn và thực nghiệm đo đạc trường trực tiếp. Thiết kế đa cấp độ (Multi-level design) được triển khai:

• Cấp độ vĩ mô (Macro-level): Khảo sát toàn diện Vịnh Bắc Bộ trên mạng lưới 12 mặt cắt vĩ tuyến ($17^\circ45'\text{N}$ đến $21^\circ30'\text{N}$, bước nhảy $15'$) và 11 mặt cắt kinh tuyến ($105^\circ45'\text{E}$ đến $108^\circ15'\text{E}$, bước nhảy $15'$).
• Cấp độ trung mô (Meso-level): Quan trắc động thái 12 tháng tại 3 vị trí trọng điểm đại diện cho các vùng sinh thái thủy văn khác nhau: Vị trí A (khu vực ven bờ Quảng Ninh - Hải Phòng), Vị trí B (khu vực giữa vịnh), Vị trí C (khu vực phía Nam vịnh - vùng nước sâu hơn tiếp giáp Biển Đông).
• Cấp độ vi mô (Micro-level): Phân tích mặt cắt thẳng đứng từng mét nước từ mặt biển ($z = 0\text{ m}$) đến đáy biển ($z = 80 - 100\text{ m}$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn quốc tế:

1. Thiết bị đo đạc và kiểm chuẩn: 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), thiết bị đo dẫn điện - nhiệt độ - độ sâu CTD, kết hợp phương pháp kiểm chuẩn cơ học tấm phẳng (Bar check) theo định kỳ. Thiết bị SVP-15 có độ chính xác vận tốc âm $\pm 0.05\text{ m/s}$, độ phân giải $0.001\text{ m/s}$.
2. Kỹ thuật tam giác đạc dữ liệu (Data Triangulation): Đối chiếu chéo 3 nguồn dữ liệu độc lập:
   • Dữ liệu lịch sử quan trắc hải văn toàn cầu WOD (World Ocean Database - NOAA).
   • Cơ sở dữ liệu hải dương học vùng Viễn Đông Nga (Far Eastern Regional Database).
   • Dữ liệu khảo sát thực tế bằng tàu đo đạc chuyên dụng của Trung tâm Trắc địa và Bản đồ Biển (SEAMAP).
3. Chuẩn hóa dữ liệu: Chuyển đổi toàn bộ các tham số nhiệt độ ($T$), độ mặn ($S$), áp suất/độ sâu ($D$) về cùng hệ quy chiếu không gian WGS-84, tính toán vận tốc âm theo thuật toán Medwin và Chen-Millero, sau đó đối sánh kiểm định sai số trung phương (RMSE).

Data và phân tích

Phân tích định lượng được thực hiện trên tập dữ liệu gồm hàng nghìn điểm đo qua nhiều năm. Các công cụ xử lý chuyên dụng bao gồm phần mềm Hệ thống Thông tin Đáy biển (Seafloor Information System - SIS), phần mềm biên tập thủy âm XTF (Extended Triton Format), các module GIS tự xây dựng và các gói phân tích thống kê.

Kiểm tra độ vững (Robustness Check) được thực hiện bằng cách so sánh độ sâu đo sâu đa tia thực tế khi áp dụng 3 phương án: (1) Sử dụng vận tốc âm mặc định $1500\text{ m/s}$; (2) Sử dụng SVP đo thực tế tại hiện trường; (3) Sử dụng SVP trích xuất từ CSDL vận tốc âm tự xây dựng của luận án. Kết quả chứng minh phương án trích xuất CSDL cho sai số độ sâu so với đo thực tế chỉ lệch $\Delta d = 0.05 - 0.22\text{ m}$, nằm hoàn toàn trong giới hạn cho phép của chuẩn IHO S-44 Order 1a.

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

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

Nghiên cứu đưa ra 5 phát hiện cốt lõi mang tính bước ngoặt cho ngành trắc địa thủy đạc Việt Nam:

1. Xác lập biên độ biến thiên vận tốc âm cực đại trong năm đạt $27.5\text{ m/s}$:

   "Vịnh Bắc Bộ là vùng nước có độ sâu nhỏ, có nhiều cửa sông do đó sự biến đổi nhiệt độ và độ mặn là các nhân tố ảnh hưởng đến quy luật biến đổi vận tốc âm theo không gian và thời gian. Đặ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." Biến thiên này chủ yếu do dao động nhiệt độ nước tầng mặt giữa mùa đông ($19^\circ\text{C}$) và mùa hè ($31^\circ\text{C}$). Nếu không cập nhật vận tốc âm theo tháng, sai số độ sâu đo đạc có thể lên tới $1.8%$ giá trị độ sâu thực tế.

2. Quy luật phân tầng âm học đối lập theo mùa (Seasonal Acoustic Inversion):

   • Về mùa hè (tháng 5 đến tháng 9): Bức xạ mặt trời mạnh tạo ra tầng nước mặt ấm ($>30^\circ\text{C}$) có vận tốc âm cao ($>1540\text{ m/s}$), bên dưới là tầng nhảy nhiệt sâu từ $15\text{ m}$ đến $35\text{ m}$, khiến vận tốc âm giảm dốc đứng xuống $1525\text{ m/s}$ ở tầng đáy. Tia âm bị bẻ cong gập xuống đáy biển (hiệu ứng khúc xạ âm bản).
   • Về mùa đông (tháng 11 đến tháng 3 năm sau): Gió mùa Đông Bắc mạnh gây xáo trộn mạnh toàn bộ cột nước nông, làm triệt tiêu tầng nhảy nhiệt. Vận tốc âm gần như đồng nhất từ mặt xuống đáy ($1512 - 1515\text{ m/s}$), tia âm truyền thẳng.

3. Hiện tượng bất đối xứng Đông - Tây trong cấu trúc trường âm: Phía Tây Vịnh Bắc Bộ (ven bờ Việt Nam) có gradient vận tốc âm biến đổi phức tạp do mạng lưới sông ngòi đổ ra lượng nước ngọt lớn làm giảm độ mặn cục bộ xuống dưới $25\text{ ppt}$ vào mùa mưa lũ. Ngược lại, phía Đông Vịnh (phía đảo Hải Nam) trường vận tốc âm ổn định hơn và có giá trị cao hơn do ảnh hưởng của khối nước biển khơi có độ mặn cao ($33 - 34.5\text{ ppt}$).

4. Hiệu ứng sai số dạng "nụ cười / cau mày" (Smile/Frown Artifact) trên chùm tia MBES: Luận án định lượng hóa sai số hình học: khi vận tốc âm mặt nước sử dụng trong bộ định dạng chùm tia (Beamformer) cao hơn thực tế, địa hình đáy quét đa tia bị uốn cong lên hai rìa ngoài (hiệu ứng "Smile"); ngược lại khi vận tốc âm thấp hơn, hai rìa bị cắm xuống (hiệu ứng "Frown"). Ở góc quét $\theta = 65^\circ$, độ sâu $50\text{ m}$, sai số vận tốc âm chỉ $5\text{ m/s}$ sẽ tạo ra sai số độ sâu ngoài rìa chùm tia lên đến $0.48\text{ m}$, vượt ngưỡng sai số cho phép của IHO.

5. Xác thực khả năng thay thế đo thực địa của CSDL vận tốc âm:

   "Quy trình lấy mẫu vận tốc âm đề xuất, cùng cơ sở dữ liệu vận tốc âm khu vực Vịnh Bắc Bộ đã xây dựng đáp ứng được các yêu cầu về độ chính xác khi ứng dụng các thiết bị thủy âm trong thành lập bản đồ địa hình đáy biển tỷ lệ 1:10.000 và nhỏ hơn trong khu vực Vịnh Bắc Bộ."

Độ sâu (m)
   0 ────────────────────────────────────────────────────────
       │ MÙA ĐÔNG (Đồng nhất)       │ MÙA HÈ (Phân tầng)
  10 ──┼────────────────────────────┼────────────────────────
       │ T = 19°C - 20°C             │ T = 30°C - 31°C, v > 1540 m/s
  20 ──┼─ v ≈ 1512 - 1515 m/s ───────┼──────────────────────── (Tầng mặt)
       │ (Cột nước xáo trộn đều)    │ ░░░ TẦNG NHẢY NHIỆT ░░░ (Thermocline)
  30 ──┼────────────────────────────┼──────────────────────── (v giảm mạnh)
       │                             │ T = 24°C - 25°C
  40 ──┼────────────────────────────┼─ v ≈ 1525 m/s ───────── (Tầng đáy)
       │                             │
  50 ──┴────────────────────────────┴────────────────────────

Implications đa chiều

• Về mặt lý thuyết: Cung cấp mô hình định lượng tương quan giữa trường khí tượng - hải văn nhiệt đới gió mùa và trường sóng âm vùng biển nông, làm tài liệu tham khảo chuẩn mực cho vật lý biển và âm học đại dương Việt Nam.
• Về mặt phương pháp luận: Đề xuất giải pháp trích xuất dữ liệu nội suy không gian kết hợp hiệu chỉnh đường truyền tia âm, mở ra phương pháp tiếp cận mới trong trắc địa biển mà không phụ thuộc hoàn toàn vào mật độ đo đạc cơ học.
• Về ứng dụng thực tiễn: Giúp các đơn vị khảo sát biển (như SEAMAP, Cục Bản đồ Quân sự, Tổng cục Biển và Hải đảo) tối ưu hóa thời gian đo đạc trên biển. Thay vì phải dừng tàu thả đầu đo SVP $2\text{ giờ/lần}$, kỹ sư có thể sử dụng CSDL kết hợp lấy mẫu kiểm chứng $1 - 2\text{ lần/ngày}$, tiết kiệm từ $20% - 30%$ thời gian ca đo của tàu khảo sát.
• Về chính sách và quy chuẩn: Cung cấp luận cứ khoa học thực tiễn để Bộ Tài nguyên và Môi trường cập nhật, sửa đổi Thông tư 27/2011/TT-BTNMT và hoàn thiện các quy chuẩn kỹ thuật quốc gia về đo vẽ địa hình đáy biển bằng công nghệ đa tia MBES.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, nghiên cứu thẳng thắn chỉ ra các hạn chế kỹ thuật và biên điều kiện:

1. Độ phân giải tại các cửa sông cục bộ: Tại các vùng cửa sông hẹp trong thời kỳ triều kiệt và xả lũ đột biến, gradient độ mặn thay đổi quá nhanh theo giờ, CSDL dạng lưới trung bình tháng chưa phản ánh hết được các biến động vi mô tức thời.
2. Thiếu số liệu chuỗi thời gian liên tục 24/7 ở tầng đáy sâu: Do hạn chế về kinh phí triển khai các trạm quan trắc hải văn ngầm cố định (mooring stations), dữ liệu tầng đáy sâu chủ yếu dựa vào các chuyến khảo sát chiến dịch và mô hình WOD.
3. Phạm vi không gian: Nghiên cứu mới giới hạn thực nghiệm tại Vịnh Bắc Bộ, chưa mở rộng toàn diện xuống vùng biển Nam Trung Bộ (nước sâu, có hiện tượng nước trồi xích đạo) và vùng biển Tây Nam (vùng nước nông Vịnh Thái Lan).

Chương trình nghiên cứu tương lai (10-year Research Agenda):

• Phát triển mô hình dự báo SVP thời gian thực (Real-time Dynamic SVP) tích hợp dữ liệu viễn thám nhiệt độ mặt biển (SST) từ vệ tinh và mô hình hoàn lưu đại dương phân giải cao (ROMS).
• Mở rộng xây dựng CSDL vận tốc âm 4D cho toàn bộ các vùng biển Việt Nam, bao gồm cả quần đảo Hoàng Sa và Trường Sa.
• Ứng dụng phương tiện ngầm tự hành (AUV/USV) trang bị cảm biến CTD/SVP để tự động hóa việc cập nhật dữ liệu liên tục vào hệ thống CSDL đám mây.

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

Luận án tạo ra tác động sâu rộng trên nhiều lĩnh vực:

• Tác động học thuật: Dự kiến đóng góp hàng chục trích dẫn cho các nghiên cứu về âm học biển nhiệt đới, trắc địa công trình ngầm và hải dương học tại khu vực Đông Nam Á.
• Chuyển đổi công nghiệp khảo sát biển: Tăng hiệu suất vận hành của các đội tàu khảo sát trắc địa, dầu khí ngoài khơi (Vietsovpetro, PTSC) và các nhà thầu thi công cáp ngầm, điện gió ngoài khơi tại Vịnh Bắc Bộ thông qua việc rút ngắn thời gian chuẩn bị và thi công ngoài thực địa.
• Quản lý nhà nước và kinh tế xã hội: Giúp tiết kiệm hàng chục tỷ đồng ngân sách hàng năm cho các dự án điều tra cơ bản tài nguyên môi trường biển của Chính phủ, đồng thời nâng cao độ chính xác các hải đồ luồng hàng hải, đảm bảo an toàn tuyệt đối cho các đội tàu vận tải quốc tế cập cảng Hải Phòng, Cái Lân, Nghi Sơn.
• An ninh quốc phòng: Dữ liệu trường âm 12 tháng tại Vịnh Bắc Bộ là cơ sở thông tin chiến lược hỗ trợ định vị âm học dưới nước, phục vụ công tác tác chiến tàu ngầm và hệ thống sonar cảnh giới bảo vệ chủ quyền biển đảo.

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

• Nghiên cứu sinh và học viên cao học ngành Trắc địa - Bản đồ, Hải dương học: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa vật lý trường âm và trắc địa hải văn, kế thừa CSDL chuẩn hóa để phát triển các đề tài chuyên sâu.
• Các nhà khoa học âm học và địa vật lý biển: Sử dụng bức tranh tổng thể về trường vận tốc âm Vịnh Bắc Bộ để mô hình hóa lan truyền sóng âm nổ địa chấn tầng nông và khảo sát sinh vật biển.
• Doanh nghiệp khảo sát R&D và nhà thầu hàng hải: Ứng dụng trực tiếp quy trình lấy mẫu và phần mềm CSDL vận tốc âm vào các gói thầu đo vẽ bản đồ địa hình đáy biển, khảo sát tuyến cáp ngầm và móng trụ điện gió ngoài khơi.
• Cơ quan quản lý nhà nước (Bộ TN&MT, Cục Đo đạc, Bản đồ và Dữ liệu địa lý Việt Nam): Sử dụng các kết quả nghiên cứu làm căn cứ khoa học ban hành định mức kinh tế kỹ thuật và quy định đo đạc thủy đạc quốc gia.

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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là xác lập mô hình tương quan phi tuyến giữa chu kỳ khí hậu nhiệt đới gió mùa và cấu trúc phân tầng vận tốc âm nước nông (Shallow-water Sound Velocity Profile), mở rộng lý thuyết truyền sóng và khúc xạ tia âm Snell trong môi trường có tầng nhảy nhiệt (thermocline) mạnh. Nghiên cứu đã chứng minh rằng tại vùng nước nông $<100\text{ m}$, biến thiên nhiệt độ theo mùa đóng vai trò chi phối tuyệt đối (gây biến thiên vận tốc âm tới $27.5\text{ m/s}$), vượt trội hoàn toàn so với ảnh hưởng của gradient áp suất thủy tĩnh theo độ sâu.

2. Đổi mới phương pháp luận của nghiên cứu khi so sánh với các công trình trước đây?

So với nghiên cứu của Dương Quốc Lương (2006) và Vũ Hồng Tập (2011) vốn chỉ đánh giá cục bộ trên thiết bị đo sâu đơn tia SBES tại một vài thời điểm, luận án đã tạo bước đột phá phương pháp luận bằng cách:

• Tích hợp đồng bộ CSDL hải dương học toàn cầu (WOD, GDEMV 3.0, Viễn Đông Nga) với chuỗi đo đạc thực nghiệm thực tế bằng thiết bị SVP-15.
• Thiết lập hệ thống 12 mặt cắt vĩ tuyến $\times$ 11 mặt cắt kinh tuyến $\times$ 12 tháng, xây dựng phần mềm chuyên dụng cho phép trích xuất tự động và bù sai số khúc xạ chùm tia MBES đa chiều.

3. Phát hiện gây bất ngờ nhất có dữ liệu thực nghiệm chứng minh là gì?

Phát hiện bất ngờ nhất là sự tương phản âm học cực đoan giữa mùa đông và mùa hè tại cùng một vị trí ở Vịnh Bắc Bộ. Vào tháng 7 (mùa hè), cột nước bị phân thành hai khối nước riêng biệt với bước nhảy vận tốc âm giảm hơn $15\text{ m/s}$ chỉ trong khoảng biến thiên độ sâu $10\text{ m}$ ($15 - 25\text{ m}$); trong khi vào tháng 1 (mùa đông), toàn bộ cột nước từ $0\text{ m}$ đến đáy biến thành một khối đồng nhất với độ chênh vận tốc âm toàn cột nước không quá $1.5\text{ m/s}$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết:

• Toàn bộ công thức toán học và bảng hệ số chuyển đổi (Medwin, Del Grosso, Chen-Millero).
• Cấu trúc cơ sở dữ liệu bảng thuộc tính, khóa chính và quan hệ thực thể (ERD) cho CSDL vận tốc âm.
• Sơ đồ quy trình lấy mẫu vận tốc âm chuẩn hóa 4 bước cho các thiết bị đo sâu đơn tia, đa tia, sonar quét sườn và định vị thủy âm.
• Bộ số liệu mẫu tại Phụ lục 1, 2, 3 cho phép các nhà nghiên cứu độc lập chạy lại thuật toán và kiểm chứng kết quả.

5. Kế hoạch nghiên cứu 10 năm được phác thảo như thế nào?

Kế hoạch 10 năm tập trung vào 3 trụ cột:

• Xây dựng hệ thống CSDL Hải văn - Vận tốc âm Quốc gia liên tục 4D bao phủ toàn bộ vùng đặc quyền kinh tế (EEZ) Việt Nam.
• Phát triển module phần mềm AI tự động nội suy và hiệu chỉnh thời gian thực cho các hệ thống MBES gắn trên robot tự hành không người lái (AUV/USV).
• Tích hợp trường âm 4D vào hệ thống phòng thủ âm học dưới nước và giám sát sinh thái rạn san hô biển.

Kết luận

1. Luận án đã làm sáng tỏ toàn diện quy luật biến đổi vận tốc âm trong môi trường nước biển Vịnh Bắc Bộ theo không gian 3 chiều và thời gian 12 tháng, chỉ ra biên độ dao động cực đại trong năm lên tới $27.5\text{ m/s}$ do tác động tổng hợp của hoàn lưu gió mùa và dòng chảy sông ngòi.
2. Xác lập cơ sở khoa học và đánh giá định lượng mức độ ảnh hưởng của sai số vận tốc âm tới độ chính xác của các thiết bị thủy âm hiện đại (SBES, MBES, SSS, USBL), giải thích nguyên nhân hình thành các dị thường uốn cong mô hình địa hình đáy biển.
3. Xây dựng thành công CSDL vận tốc âm chuẩn hóa khu vực Vịnh Bắc Bộ dựa trên tích hợp nguồn dữ liệu quốc tế uy tín (WOD, Viễn Đông Nga) và dữ liệu đo thực tế SVP-15, đi kèm chương trình phần mềm trích xuất dữ liệu thân thiện, chính xác.
4. Đề xuất quy trình lấy mẫu vận tốc âm hoàn thiện, khoa học, khắc phục triệt để các hạn chế của quy trình hiện hành theo Thông tư 27/2011/TT-BTNMT, giảm thiểu số lần thả máy đo thực địa mà vẫn đáp ứng nghiêm ngặt tiêu chuẩn bản đồ tỷ lệ 1:10.000 (với sai số trung phương $m_D < 0.35\text{ m}$).
5. Tạo lập bước chuyển đổi mô thức (Paradigm Shift) quan trọng cho ngành kỹ thuật trắc địa - bản đồ biển Việt Nam từ phương pháp hiệu chỉnh tĩnh truyền thống sang phương pháp hiệu chỉnh động dựa trên nền tảng cơ sở dữ liệu số hóa 4D.
6. Mở ra 3 hướng nghiên cứu mới đầy tiềm năng: Dự báo trường âm thời gian thực tích hợp viễn thám; Tự động hóa đo đạc thủy đạc bằng phương tiện không người lái (USV/AUV); Ứng dụng âm học biển trong chiến lược quốc phòng an ninh và kinh tế biển bền vững.