Tổng quan về luận án

Nghiên cứu cấu trúc trường trọng lực biển đóng vai trò đặc biệt quan trọng trong trắc địa vật lý, hải dương học động lực, thăm dò tài nguyên khoáng sản biển và bảo đảm quốc phòng - an ninh. Đối với Việt Nam – quốc gia có bờ biển dài trên 3.260 km và thềm lục địa rộng lớn gấp ba lần diện tích đất liền, việc xây dựng mô hình trọng lực biển có độ chính xác và độ phân giải cao là nhiệm vụ chiến lược cấp thiết. Tuy nhiên, phương pháp đo trọng lực truyền thống bằng tàu biển hoặc máy bay bộc lộ nhiều hạn chế về chi phí vận hành khổng lồ, thời gian thực hiện kéo dài và điều kiện khí tượng thủy văn phức tạp. Trong bối cảnh đó, kỹ thuật đo cao vệ tinh radar (Satellite Altimetry) đã mở ra cuộc cách mạng trong việc xác định hình dạng Geoid và dị thường trọng lực biển ($\Delta g$) trên quy mô toàn cầu và khu vực.

Khoảng trống nghiên cứu (Research Gap) lớn nhất hiện nay nằm ở độ chính xác của số liệu đo cao tại các vùng biển nông ven bờ. Khu vực Vịnh Bắc Bộ mang đặc thù địa hình đáy biển thoải, nông, chịu ảnh hưởng phức tạp của các yếu tố động lực biển (thủy triều biên độ lớn, sóng, hoàn lưu gió mùa) và hiện tượng phản xạ radar bị nhiễu do đường bờ. Các mô hình trường trọng lực toàn cầu hiện hành như DTU10GRA, DTU13GRA, DTU15GRA, DTU17GRA của Viện Không gian Quốc gia Đan Mạch (DTU) hay VGG của Sandwell & Smith (2014, 2017) khi áp dụng tại Biển Đông và Vịnh Bắc Bộ chỉ đạt độ chính xác dao động từ $\pm 5,8\text{ mGal}$ đến $\pm 9,0\text{ mGal}$. Tại Việt Nam, các nghiên cứu trước đây (Nguyễn Như Trung, 2003; Bùi Công Quế, 2008; Nguyễn Văn Sáng, 2012) chủ yếu khai thác mô hình toàn cầu có sẵn hoặc sử dụng đơn vệ tinh ở chế độ lặp lại chính xác (Exact Repeat Mission - ERM) như ENVISAT với mật độ vết đo thưa, độ phân giải thấp ($10' \times 10'$) và sai số đạt khoảng $\pm 6,8\text{ mGal}$.

Luận án tiến sĩ kỹ thuật của tác giả Phạm Văn Tuyên (Trường Đại học Mỏ - Địa chất Hà Nội, 2023) với đề tài: "Nghiên cứu giải pháp nâng cao độ chính xác xác định dị thường trọng lực bằng số liệu đo cao vệ tinh trên vùng biển Vịnh Bắc Bộ - Việt Nam" (Mã số: 9520503) đã giải quyết triệt để bài toán này. Nghiên cứu tập trung vào 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học cốt lõi:

  • RQ1: Làm thế nào để tối ưu hóa mật độ và độ chính xác của vết đo cao vệ tinh trong vùng nước nông ven bờ Vịnh Bắc Bộ?
  • RQ2: Mô hình thế trọng trường toàn cầu (GGM) và mô hình mặt biển trung bình động học (MDT) nào tương thích tối ưu nhất với cấu trúc trọng trường lãnh thổ và vùng biển Việt Nam?
  • RQ3: Quy trình toán học nào cho phép chuyển đổi phần dư độ cao Geoid sang dị thường trọng lực và kết hợp hiệu quả với số liệu đo trọng lực tàu biển trực tiếp?
  • Giả thuyết H1: Việc tích hợp đồng thời các vệ tinh thế hệ mới ở chế độ nhiệm vụ trắc địa (Geodetic Mission - GM) bao gồm CryoSat-2/GM (băng tần Ku, công nghệ SAR/Interferometry) và SARAL-DP/AltiKa/GM (băng tần Ka 35.75 GHz) sẽ tăng mật độ vết quét lên gấp 5-10 lần, giảm triệt để nhiễu tầng điện ly và dạng sóng bờ.
  • Giả thuyết H2: Mô hình EIGEN-6C4 kết hợp với DTU15MDT cung cấp cấu trúc bước sóng dài và trung bình chính xác nhất cho kỹ thuật Loại bỏ - Phục hồi (Remove-Compute-Restore - RCR) tại Vịnh Bắc Bộ.
  • Giả thuyết H3: Thuật toán Collocation bình phương nhỏ nhất (Least Squares Collocation - LSC) kết hợp giải pháp làm khớp (fitting) với số liệu đo tàu sẽ triệt tiêu sai số hệ thống cục bộ, nâng độ chính xác dị thường trọng lực biển lên tốt hơn $\pm 3,0\text{ mGal}$.

Khung lý thuyết của nghiên cứu dựa trên lý thuyết thế trọng trường Trái Đất (Helmert, 1901; Heiskanen & Moritz, 1967), lý thuyết xác định Geoid mặt biển (Bruns, Stokes, Vening-Meinesz) và lý thuyết ước lượng thống kê tối ưu Collocation bình phương nhỏ nhất (Krarup, 1969; Moritz, 1980). Phạm vi nghiên cứu bao quát toàn bộ vùng biển Vịnh Bắc Bộ thuộc chủ quyền Việt Nam (tọa độ từ $17,8^\circ\text{N}$ đến $21,8^\circ\text{N}$; $105,5^\circ\text{E}$ đến $110,0^\circ\text{E}$). Dữ liệu thực nghiệm bao gồm 105 chu kỳ đo của CryoSat-2/GM, 54 chu kỳ đo của SARAL-DP/AltiKa/GM, 818 điểm song trùng GPS - Thủy chuẩn hạng I-II-III, 31 trạm nghiệm triều quốc gia, cùng 2.011 điểm dị thường trọng lực đo trực tiếp bằng máy đo trọng lực biển ZLS Dynamic Gravity Meter D06. Luận án đã tạo ra bước đột phá khi xây dựng lưới dị thường trọng lực độ phân giải $1' \times 1'$ với độ chính xác đạt $\pm 2,98\text{ mGal}$ (vượt mục tiêu thiết kế $\pm 4,0\text{ mGal}$).

Literature Review và Positioning

Nghiên cứu xác định trọng trường biển từ đo cao vệ tinh trên thế giới trải qua bốn thập kỷ phát triển mạnh mẽ với hai trường phái tính toán và xử lý tín hiệu chính:

  1. Trường phái biến đổi Fourier nhanh (FFT) và công thức Vening-Meinesz nghịch đảo: Sandwell & Smith (1997, 2014) thuộc Viện Hải dương học Scripps (Hoa Kỳ) đã ứng dụng biến đổi FFT trên lưới gradient độ lệch dây dọi để tính toán dị thường trọng lực toàn cầu từ dữ liệu Geosat/GM và ERS-1/GM, đạt độ chính xác $\pm 4,0\text{ mGal}$ đến $\pm 7,0\text{ mGal}$. Xiaoli Deng, Will Featherstone, Cheinway Hwang và C. Shum (2001) áp dụng công thức Vening-Meinesz nghịch đảo kết hợp xử lý lại dạng sóng radar (waveform retracking) tại eo biển Đài Loan ($20^\circ \le \text{B} \le 42^\circ\text{N}$; $110^\circ \le \text{L} \le 129^\circ\text{E}$), ghi nhận sai số $\pm 13,8\text{ mGal}$ do ảnh hưởng phức tạp của địa hình bờ. Liang Liu và cộng sự (2016) sử dụng công thức Stokes nghịch đảo với 7 năm dữ liệu Jason-1 tại Biển Đông ghi nhận sai số trung phương đạt $\pm 12,3\text{ mGal}$.
  2. Trường phái Collocation bình phương nhỏ nhất (LSC) và chuỗi mô hình DTU: Ole Baltazar Andersen và Per Knudsen (DTU Space, Đan Mạch) đã tiên phong phát triển các mô hình DNSC08GRA (Andersen et al., 2010; RMS $\pm 3,91\text{ mGal}$ tại Bắc Đại Tây Dương), DTU10GRA (Andersen, 2010; RMS $\pm 3,82\text{ mGal}$), DTU13GRA (Andersen et al., 2013; RMS $\pm 3,84\text{ mGal}$), DTU15GRA (Andersen & Knudsen, 2016; RMS $\pm 5,61\text{ mGal}$ tại Bắc Cực) và DTU17GRA (Andersen & Knudsen, 2017; RMS $\pm 3,43\text{ mGal}$). Trường phái này chứng minh LSC có ưu thế vượt trội trong việc xử lý phân bố điểm đo bất đối xứng, lọc nhiễu ngẫu nhiên và tích hợp linh hoạt các loại đại lượng trọng trường khác nhau thông qua hàm hiệp phương sai thực nghiệm.
                    TIẾN TRÌNH NGHIÊN CỨU DỊ THƯỜNG TRỌNG LỰC BIỂN

Trong nước, nghiên cứu của Bùi Công Quế (2008) kết hợp dữ liệu đo tàu Gagarinsky, L'Attalante với dữ liệu vệ tinh chỉ đạt độ chính xác $\pm 8,5\text{ mGal}$. Luận án của Nguyễn Văn Sáng (2012) tại Đại học MIIGAiK (Liên bang Nga) lần đầu tiên tự tính toán dị thường trọng lực Biển Đông từ dữ liệu ENVISAT/ERM nhưng sai số vẫn ở mức $\pm 6,8\text{ mGal}$. Đề tài cấp Bộ B2016-MDA-11ĐT (Vũ Văn Trí, Phạm Văn Tuyên, 2017) áp dụng CryoSat-2/GM cho vùng nước sâu Trường Sa đạt sai số nội bộ $\pm 0,67\text{ mGal}$ trên 625 điểm kiểm tra, song phương pháp chưa được hoàn thiện cho vùng biển nông ven bờ có độ phức tạp cao như Vịnh Bắc Bộ.

Vị thế học thuật (Academic Positioning) của luận án được xác lập rõ ràng: Đây là công trình đầu tiên tại Việt Nam thiết lập một quy trình công nghệ tích hợp 5 giải pháp đồng bộ (Đa vệ tinh GM dải sóng kép Ku/Ka, mô hình GGM tối ưu EIGEN-6C4, mô hình địa hình mặt biển động học DTU15MDT, phương pháp LSC với hàm hiệp phương sai chuyên biệt vùng và kỹ thuật làm khớp LSC với số liệu đo tàu thực tế), giải quyết triệt để "nút thắt" suy giảm độ chính xác của đo cao vệ tinh tại các vùng vịnh nông ven bờ.

Đó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à hoàn thiện lý thuyết biểu diễn trường trọng lực cục bộ theo mô hình thế nhiễu loạn $T$ trong không gian Hilbert có hạt nhân tái tạo (Reproducing Kernel Hilbert Space). Cụ thể, nghiên cứu đã phát triển sâu sắc kỹ thuật Loại bỏ - Phục hồi (Remove-Compute-Restore - RCR) trong việc phân tích phổ tần số của trường trọng lực biển nông:

$$\Delta g = \Delta g_{GGM} + \Delta g_{res} + \Delta g_{top}$$

Trong đó, phần bước sóng dài ($\Delta g_{GGM}$) được mô hình hóa chính xác từ khai triển điều hòa cầu bậc và thứ cao ($n_{max} = 2190$); phần bước sóng ngắn cục bộ ($\Delta g_{res}$) được trích xuất từ phần dư độ cao Geoid ($h_{res} = SSH_{corr} - h_{MDT} - N_{GGM}$) thông qua toán tử tích phân Collocation bình phương nhỏ nhất (LSC):

$$\Delta g_{res} = \mathbf{C}{\Delta g, h} (\mathbf{C}{h, h} + \mathbf{D})^{-1} \mathbf{h}_{res}$$

Đóng góp lý thuyết mang tính bản lề là việc chứng minh tính tương thích và sự cần thiết phải hiệu chỉnh hàm hiệp phương sai thực nghiệm (Empirical Covariance Function) phù hợp với cấu trúc địa chất Vịnh Bắc Bộ. Thay vì áp dụng trực tiếp các tham số hàm hiệp phương sai toàn cầu mặc định của Tscherning-Rapp, tác giả đã tái tham số hóa bán kính nội suy suy giảm $R_B$ (bán kính cầu Bjerhammar) và phương sai tín hiệu cục bộ $C_0$, giúp phản ánh chính xác ranh giới bất liên tục của móng kết tinh và trầm tích Đệ Tam tại bồn trũng Sông Hồng.

       KHUNG PHÂN TÍCH TỔNG HỢP 5 GIẢI PHÁP ĐỒNG BỘ CỦA LUẬN ÁN

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

Khung phân tích của luận án tích hợp đa chiều ba lý thuyết nền tảng:

  1. Lý thuyết tương tác sóng điện từ bề mặt: Đánh giá ưu thế vượt trội của bước sóng Ka-band (SARAL/AltiKa, tần số 35.75 GHz) với dấu chân radar (footprint) xung phản xạ nhỏ hơn ($\approx 8\text{ km}$ so với $15-20\text{ km}$ của băng Ku truyền thống), giúp hạn chế hiện tượng méo dạng sóng khi quét gần các đảo (đảo Cát Bà, Bạch Long Vĩ, Cô Tô) và đường bờ Vịnh Bắc Bộ.
  2. Lý thuyết phân giải đa quy mô mô hình trường trọng lực: Đánh giá định lượng 4 mô hình thế trọng trường toàn cầu tiên tiến nhất (EGM2008, GECO, EIGEN-6C4, SGG-UGM-1) dựa trên tập dữ liệu chuẩn 818 điểm GPS - Thủy chuẩn trên đất liền và 31 trạm nghiệm triều ven biển.
  3. Lý thuyết bình sai thống kê không gian: Thiết lập quy trình làm khớp 2 bước (two-step fitting) bằng LSC, cho phép tích hợp số liệu đo trọng lực tàu biển ZLS phân bố không đều vào lưới đo cao vệ tinh mà không gây ra hiện tượng méo dạng hình học bề mặt (geometric distortion).

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan thực chứng toán học (Positivism) và phương pháp luận đo đạc trắc địa vật lý thực nghiệm (Empirical Geodetic Methodology). Quy trình nghiên cứu được thiết kế dạng chuỗi khép kín với các tiêu chuẩn kiểm soát sai số nghiêm ngặt ở từng giai đoạn:

[Số liệu đo cao radar CryoSat-2 + SARAL]

Quy trình nghiên cứu rigorous

Quy trình xử lý số liệu đo cao vệ tinh được triển khai qua các bước chuẩn hóa nghiêm ngặt:

  1. Hiệu chỉnh khoảng cách và đường truyền tín hiệu: Khoảng cách đo thực tế $d$ được chuẩn hóa qua công thức: $$SSH_{corr} = H - (d + \Delta d_{inst} + \Delta d_{dry} + \Delta d_{wet} + \Delta d_{iono} + \Delta d_{ssb} + \Delta d_{tide} + \Delta d_{atm})$$

    • Hiệu chỉnh tầng điện ly ($\Delta d_{iono}$): Sử dụng mô hình GIM (Global Ionospheric Modeling) kết hợp dữ liệu song tần DORIS/PRARE, biên độ hiệu chỉnh từ $2 \div 20\text{ cm}$ ($\pm 3\text{ cm}$).
    • Hiệu chỉnh tầng đối lưu khô ($\Delta d_{dry}$): Tính toán từ mô hình áp suất khí quyển bề mặt ECMWF theo công thức Saastamoinen: $$\Delta d_{dry} = -2,277 \times 10^{-3} P_0 (1 + 0,0026 \cos 2\varphi)$$ Giá trị hiệu chỉnh đạt từ $2,2 \div 2,3\text{ m}$ với sai số cực nhỏ ($<1\text{ cm}$).
    • Hiệu chỉnh tầng đối lưu ướt ($\Delta d_{wet}$): Sử dụng mô hình tích phân hơi nước ECMWF kết hợp máy đo bức xạ vi sóng trên vệ tinh (Microwave Radiometer), biên độ hiệu chỉnh từ $5 \div 50\text{ cm}$.
    • Hiệu chỉnh địa vật lý và bề mặt phản xạ: Áp dụng mô hình thủy triều đại dương FES2014, hiệu chỉnh tải triều (Ocean Tide Loading), triều Trái Đất cứng (Solid Earth Tide), triều cực (Pole Tide) và hiệu ứng độ lệch trạng thái biển (Sea State Bias - SSB).

  2. Bình sai mạng lưới điểm giao cắt (Crossover Adjustment): Tại các giao điểm giữa cung thăng (ascending track) và cung giáng (descending track), chênh lệch độ cao $\Delta h = SSH_{asc} - SSH_{desc}$ do sai số quỹ đạo vệ tinh và biến thiên động học đại dương $h_t$ được khử triệt để thông qua mô hình hàm đa thức thời gian kết hợp bình sai điều kiện có số ẩn số.

Data và phân tích

Tác giả đã tự lập trình thành công gói phần mềm Geomat2015 (giao diện đồ họa chuyên dụng, chạy trên nền tảng Windows/MATLAB) để tính toán chính xác độ cao Geoid ($N$), dị thường độ cao ($\zeta$) và dị thường trọng lực chuẩn ($\Delta g_0$) từ hệ số hàm điều hòa cầu $\bar{C}{nm}, \bar{S}{nm}$ lên đến bậc và thứ 2190:

$$V(r, \theta, \lambda) = \frac{GM}{r} \left[ 1 + \sum_{n=2}^{n_{max}} \left(\frac{a}{r}\right)^n \sum_{m=0}^n (\bar{C}{nm} \cos m\lambda + \bar{S}{nm} \sin m\lambda) \bar{P}_{nm}(\cos\theta) \right]$$

Đồng thời, luận án sử dụng thành thạo gói phần mềm trắc địa vật lý tiêu chuẩn quốc tế GRAVSOFT của DTU/Đan Mạch với các module cốt lõi:

  • GEOIP: Nội suy độ cao địa hình mặt biển trung bình động học $h_{MDT}$ từ mô hình lưới DTU15MDT ($1' \times 1'$).
  • EMPCOV: Tính toán hàm hiệp phương sai thực nghiệm không gian giữa các phần dư độ cao Geoid $h_{res}$ theo các dải khoảng cách cung $\psi$.
  • COVFIT: Khớp hàm hiệp phương sai lý thuyết vào hàm thực nghiệm nhằm xác định chính xác các tham số trọng trường cục bộ (độ sâu Bjerhammar $R_B$, phương sai nhiễu $C_0$).
  • GEOCOL: Giải bài toán Collocation bình phương nhỏ nhất để tính chuyển ma trận $h_{res} \to \Delta g_{res}$ và thực hiện làm khớp số liệu với 2.011 điểm đo trực tiếp bằng máy ZLS Dynamic Gravity Meter D06.

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

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

      SO SÁNH SAI SỐ TRUNG PHƯƠNG (RMS) GIỮA CÁC MÔ HÌNH VÀ PHƯƠNG ÁN
  1. Xác định mô hình thế trọng trường toàn cầu tối ưu cho Việt Nam: Qua thực nghiệm so sánh dị thường độ cao tính từ 4 mô hình GGM với 818 điểm GPS - Thủy chuẩn và 31 trạm nghiệm triều, mô hình EIGEN-6C4 (kết hợp dữ liệu vệ tinh GOCE, GRACE, LAGEOS và dữ liệu trọng lực bề mặt DTU) đạt độ chính xác cao nhất với sai số trung phương $RMS = \pm 0,184\text{ m}$, vượt trội so với EGM2008 ($RMS = \pm 0,212\text{ m}$), GECO ($RMS = \pm 0,205\text{ m}$) và SGG-UGM-1 ($RMS = \pm 0,198\text{ m}$).
  2. Ưu thế vượt trội của việc kết hợp đa vệ tinh chế độ GM: Phương án kết hợp 105 chu kỳ CryoSat-2/GM với 54 chu kỳ SARAL-DP/AltiKa/GM (Phương án 2) cho mật độ điểm vết quét dày đặc và độ tin cậy vượt trội so với việc sử dụng đơn vệ tinh CryoSat-2/GM (Phương án 1). Khoảng cách giữa các vết quét giảm xuống dưới $3\text{ km}$, cho phép tái hiện chi tiết các cấu trúc dị thường dải hẹp.
  3. Hiệu quả đột phá của giải pháp làm khớp LSC với số liệu đo tàu: Khi chưa làm khớp, Phương án 2 đạt sai số so với số liệu đo trực tiếp là $\pm 3,87\text{ mGal}$. Sau khi áp dụng giải pháp 5 (làm khớp LSC với 2.011 điểm đo trực tiếp ZLS trên vùng biển Vịnh Bắc Bộ), sai số trung phương của Phương án 3 giảm mạnh xuống $\pm 2,98\text{ mGal}$, triệt tiêu hoàn toàn độ lệch hệ thống (mean bias chuyển từ $+1,24\text{ mGal}$ về xấp xỉ $0,00\text{ mGal}$).
  4. Vượt trội so với các mô hình toàn cầu tiên tiến nhất: Kết quả của luận án ($\pm 2,98\text{ mGal}$) vượt trội hoàn toàn so với các mô hình quốc tế áp dụng tại Vịnh Bắc Bộ: DTU10GRA ($\pm 5,82\text{ mGal}$), DTU13GRA ($\pm 5,45\text{ mGal}$), DTU15GRA ($\pm 4,96\text{ mGal}$) và DTU17GRA ($\pm 4,78\text{ mGal}$).

Implications đa chiều

  • Về mặt học thuật: Cung cấp luận cứ khoa học vững chắc khẳng định phương pháp LSC kết hợp dữ liệu đo cao đa tần số Ku/Ka ở chế độ GM là công cụ tối ưu để mô hình hóa trường trọng lực vùng biển cận bờ và vịnh nông có điều kiện khí tượng - hải dương phức tạp.
  • Về mặt thực tiễn và kinh tế biển: Bản đồ dị thường trọng lực Vịnh Bắc Bộ tỷ lệ lớn ($1' \times 1'$) với độ chính xác cao là cơ sở dữ liệu địa vật lý quan trọng phục vụ phân tích cấu trúc móng, phân định các đứt gãy kiến tạo sâu (đứt gãy Sông Hồng, đứt gãy Vịnh Bắc Bộ) và đánh giá tiềm năng bồn trũng trầm tích dầu khí Sông Hồng.
  • Về mặt quốc phòng - an ninh: Cung cấp các tham số thế trọng trường chuẩn xác phục vụ hiệu chỉnh hệ thống dẫn đường quán tính (INS) cho tàu ngầm, phương tiện tự hành ngầm dưới biển (UAV/UUV) và hiệu chỉnh quỹ đạo bay của các khí tài quân sự tầm xa trên biển.

Limitations và Future Research

  1. Vùng mù dữ liệu sát bờ (Coastal Blind Zone): Trong phạm vi sát bờ từ $0 \div 5\text{ km}$, tín hiệu radar đo cao vẫn bị ảnh hưởng bởi hiện tượng phản xạ đất liền (land contamination), dẫn đến số lượng xung phản xạ hợp lệ bị suy giảm.
  2. Chưa áp dụng kỹ thuật Retracking chuyên biệt cho sóng thô: Luận án sử dụng dữ liệu đo cao sản phẩm mức 2 (Level-2 GDR) đã qua tiền xử lý chuẩn của CNES/ESA/AVISO. Việc tự áp dụng các thuật toán Retracking nâng cao (như ALES - Adaptive Leading Edge Subwaveform hay SAMOSA) trực tiếp từ dữ liệu thô dạng sóng (Level-1b Sensor Waveforms) hứa hẹn sẽ cải thiện độ chính xác thêm $0,3 \div 0,5\text{ mGal}$ tại dải nước sát bờ.
  3. Mật độ số liệu đo tàu kiểm chứng phân bố chưa đồng đều: Tập 2.011 điểm đo trọng lực ZLS phân bố tập trung chủ yếu ở phần phía Tây và trung tâm vịnh, khu vực giáp ranh đường phân định phía Đông có mật độ điểm đo trực tiếp thưa hơn.

Định hướng nghiên cứu tương lai:

  • Tích hợp dữ liệu từ các thế hệ vệ tinh đo cao radar giao thoa góc quét rộng mới nhất như SWOT (Surface Water and Ocean Topography - NASA/CNES) và Sentinel-6 Michael Freilich.
  • Ứng dụng thuật toán Trí tuệ nhân tạo (AI/Machine Learning) trong việc tối ưu hóa nội suy hàm hiệp phương sai phi tham số tại các cấu trúc địa chất phức tạp.
  • Mở rộng quy trình 5 giải pháp công nghệ cho toàn bộ vùng đặc quyền kinh tế (EEZ) Biển Đông của Việt Nam, bao gồm khu vực Quần đảo Hoàng Sa và Trường Sa.

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

               MA TRẬN TÁC ĐỘNG VÀ ĐỐI TƯỢNG HƯỞNG LỢI
  • Tác động học thuật quốc tế: Công trình của tác giả đã được công bố đồng tác giả với GS.TS. Ole Baltazar Andersen (DTU Space) trên tạp chí quốc tế ISI-Q1 uy tín Advances in Space Research, khẳng định năng lực hội nhập khoa học trắc địa đỉnh cao của Việt Nam.
  • Tối ưu hóa ngân sách điều tra cơ bản: Việc sử dụng hiệu quả dữ liệu đo cao vệ tinh kết hợp đo tàu cục bộ giúp giảm thiểu hơn 70% chi phí khảo sát trọng lực biển bằng tàu chuyên dụng trên diện tích hàng chục nghìn kilomet vuông Vịnh Bắc Bộ.
  • Thực thi chủ quyền và phân định ranh giới: Dữ liệu trọng trường chính xác cao tại 21 điểm phân định lãnh hải, vùng đặc quyền kinh tế và thềm lục địa giữa Việt Nam và Trung Quốc trong Vịnh Bắc Bộ (theo Hiệp định ký năm 2000) góp phần hoàn thiện cơ sở dữ liệu pháp lý và kỹ thuật phục vụ quản lý biển đảo quốc gia.

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

  1. Nghiên cứu sinh và Giảng viên ngành Trắc địa - Địa vật lý: Tiếp cận quy trình công nghệ tính toán hiện đại, mã nguồn phần mềm Geomat2015 và phương pháp luận ước lượng tối ưu Collocation bình phương nhỏ nhất.
  2. Tập đoàn Dầu khí Việt Nam (PetroVietnam) và các Liên doanh dầu khí: Sử dụng dữ liệu dị thường trọng lực có độ phân giải $1' \times 1'$ để minh giải cấu trúc địa chất khu vực, phát hiện các bẫy dầu khí phi truyền thống tại Vịnh Bắc Bộ.
  3. Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam (Bộ Tài nguyên và Môi trường): Ứng dụng trực tiếp kết quả nghiên cứu để xây dựng quy chuẩn kỹ thuật quốc gia về xác định Geoid biển và trọng trường biển bằng công nghệ vệ tinh.
  4. Bộ Tư lệnh Hải quân và Bộ Quốc phòng: Khai thác mô hình trọng trường biển độ chính xác cao phục vụ bảo đảm hàng hải quân sự và định vị các công trình quốc phòng trên biển.

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ì?

Đóng góp độc đáo nhất là việc chuẩn hóa và tái thiết lập hàm hiệp phương sai thực nghiệm chuyên biệt cho vùng biển nông cận bờ Vịnh Bắc Bộ trong không gian Collocation bình phương nhỏ nhất. Luận án đã giải quyết thành công mâu thuẫn giữa lý thuyết trường trọng lực đại dương mở (vốn coi $h_{MDT}$ là ổn định và sóng phản xạ đồng nhất) với thực tế động lực biển phức tạp tại vịnh nông, mở rộng phạm vi áp dụng của mô hình tích phân Krarup-Moritz vào vùng biển ven bờ Việt Nam.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu quốc tế trước đây thể hiện ở điểm nào?

So với nghiên cứu của Sandwell & Smith (sử dụng thuần túy biến đổi FFT) và Andersen & Knudsen (sử dụng LSC trên quy mô toàn cầu), luận án đã thiết lập quy trình tích hợp 2 bước: (1) Kết hợp dữ liệu vệ tinh đa tần số Ku-band (CryoSat-2/GM) và Ka-band (SARAL-DP/GM) để triệt tiêu nhiễu bề mặt và tăng mật độ vết quét; (2) Sử dụng chính số liệu đo trọng lực trực tiếp bằng tàu biển (máy ZLS) làm dữ liệu chuẩn để thực hiện bài toán làm khớp ngược (fitting LSC), loại bỏ triệt để sai số hệ thống mà các mô hình toàn cầu không thể xử lý được.

3. Phát hiện bất ngờ nhất từ kết quả thực nghiệm là gì?

Phát hiện bất ngờ nhất là mô hình thế trọng trường toàn cầu mới SGG-UGM-1 (do Trung Quốc phát triển năm 2018) và GECO dù tích hợp dữ liệu trọng lực mới nhưng vẫn cho độ lệch chuẩn tại Việt Nam ($RMS = \pm 0,198\text{ m}$ và $\pm 0,205\text{ m}$) kém hơn so với mô hình EIGEN-6C4 ($RMS = \pm 0,184\text{ m}$). Điều này chứng minh sự kết hợp tối ưu giữa trường trọng lực không gian bậc cao GOCE với số liệu mặt đất trong mô hình EIGEN-6C4 có tính tương thích vượt trội với cấu trúc vỏ Trái Đất khu vực Đông Nam Á.

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

Có. Quy trình nghiên cứu được sơ đồ hóa chi tiết từ khâu tiền xử lý dữ liệu đo cao GDR (loại trừ điểm lỗi, hiệu chỉnh khí quyển GIM/ECMWF, hiệu chỉnh triều FES2014), cấu trúc tham số phần mềm Geomat2015, định dạng tệp đầu vào/đầu ra cho các module GRAVSOFT (EMPCOV, COVFIT, GEOCOL), cùng bảng tham số hàm hiệp phương sai lý thuyết, bảo đảm khả năng tái lập và kiểm chứng độc lập hoàn toàn.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình nghiên cứu dài hạn bao gồm: (1) Xây dựng mô hình Geoid biển thống nhất độ chính xác milimet cho toàn Biển Đông; (2) Tích hợp dữ liệu đo cao công nghệ SAR Altimetry và Lidar từ các vệ tinh ICESat-2, SWOT, Sentinel-6; (3) Thiết lập hệ thống giám sát biến thiên trọng trường và mực nước biển dâng thời gian thực phục vụ ứng phó biến đổi khí hậu tại dải ven biển Việt Nam.

Kết luận

  1. Xác lập cơ sở khoa học hoàn chỉnh: Luận án đã xây dựng thành công hệ thống phương pháp luận và quy trình tính toán đồng bộ, chặt chẽ về lý thuyết và khả thi tuyệt đối trong thực nghiệm để xác định dị thường trọng lực biển từ dữ liệu đo cao vệ tinh tại vùng biển nông cận bờ.
  2. Lựa chọn mô hình nền tảng tối ưu: Đã kiểm chứng khoa học và lựa chọn chính xác mô hình thế trọng trường toàn cầu EIGEN-6C4 ($n_{max}=2190$) và mô hình địa hình mặt biển trung bình động học DTU15MDT ($1' \times 1'$) tương thích tối ưu nhất với lãnh thổ và vùng biển Việt Nam.
  3. Khai thác đột phá dữ liệu đa vệ tinh Geodetic Mission: Tích hợp thành công 105 chu kỳ đo của CryoSat-2/GM và 54 chu kỳ đo của SARAL-DP/AltiKa/GM, tạo ra mạng lưới vết quét có mật độ dày đặc và độ chính xác vượt bậc.
  4. Đạt độ chính xác kỷ lục tại Vịnh Bắc Bộ: Kết quả dị thường trọng lực nhận được đạt sai số trung phương $RMS = \pm 2,98\text{ mGal}$ trên lưới phân giải $1' \times 1'$, vượt xa mục tiêu đề ra ($\pm 4,0\text{ mGal}$) và vượt trội so với các mô hình toàn cầu DTU10/13/15/17GRA.
  5. Đóng góp sản phẩm thực tiễn giá trị cao: Tác giả đã tự phát triển phần mềm chuyên dụng Geomat2015 và cung cấp bộ cơ sở dữ liệu dị thường trọng lực số hóa Vịnh Bắc Bộ, phục vụ trực tiếp cho công tác điều tra cơ bản, phát triển kinh tế biển, bảo vệ chủ quyền và an ninh quốc phòng quốc gia.