Tổng quan về luận án

Xác định trường trọng lực trên các vùng biển và thềm lục địa đóng vai trò then chốt trong việc hoàn thiện hệ thống toạ độ - độ cao quốc gia, nghiên cứu cấu trúc địa chất manti - vỏ Trái Đất, thăm dò khoáng sản dầu khí và phục vụ an ninh quốc phòng biển đảo. Tại Việt Nam, với diện tích vùng biển rộng gấp ba lần đất liền, các phương pháp đo trọng lực truyền thống bằng tàu biển (shipborne gravimetry) hoặc hàng không (airborne gravimetry) đòi hỏi kinh phí khổng lồ, thời gian thi công kéo dài và gặp nhiều trở ngại trước điều kiện khí tượng - thủy văn phức tạp. Kỹ thuật đo cao vệ tinh (satellite altimetry) ra đời như một giải pháp mang tính cách mạng cho phép thu thập dữ liệu độ cao mặt biển (Sea Surface Height - SSH) trên quy mô toàn cầu với độ bao phủ dày đặc. Tuy nhiên, tại các vùng biển nông, bán khép kín và có hình viền bờ phức tạp như Vịnh Bắc Bộ, độ chính xác xác định dị thường trọng lực từ các mô hình toàn cầu hiện hành thường bị suy giảm nghiêm trọng do nhiễu dạng sóng ven bờ và sai số mô hình hóa hải dương học.

Nghiên cứu của tác giả Phạm Văn Tuyên (Trường Đại học Mỏ - Địa chất) đã giải quyết triệt để khoảng trống nghiên cứu (research gap) này thông qua việc xây dựng một hệ thống giải pháp toàn diện kết hợp dữ liệu đo cao vệ tinh thế hệ mới với các kỹ thuật xử lý toán học tiên tiến. Cơ sở lý luận của công trình bắt nguồn từ lý thuyết trắc địa vật lý cổ điển và hiện đại, kế thừa phương trình cơ bản về dị thường trọng lực $ \Delta g = g - \gamma $ (với $g$ là trọng lực thực đo và $\gamma$ là trọng lực chuẩn) kết hợp phương trình xác định độ cao mặt biển hiệu chỉnh $\text{SSH}_{\text{corr}} = H - (d + e)$.

Luận án thiết lập hệ thống 4 câu hỏi nghiên cứu ($RQ$) và 3 giả thuyết khoa học ($H$) trọng tâm:

  • $RQ_1$: Làm thế nào để tối ưu hóa việc tích hợp dữ liệu đo cao vệ tinh từ các nhiệm vụ trắc địa (Geodetic Mission - GM) đa nguồn nhằm nâng cao mật độ vết quét trên vùng biển Vịnh Bắc Bộ?
  • $RQ_2$: Mô hình thế trọng trường toàn cầu (Global Geopotential Model - GGM) và mô hình mặt biển trung bình động học (Mean Dynamic Topography - MDT) nào đạt độ chính xác và tương thích tối ưu nhất với điều kiện địa vật lý Việt Nam?
  • $RQ_3$: Phương pháp toán học nào giải quyết hiệu quả bài toán chuyển đổi từ phần dư độ cao geoid sang phần dư dị thường trọng lực trong điều kiện phân bố điểm bất quy tắc tại vùng biển nông?
  • $RQ_4$: Quy trình làm khớp dữ liệu dị thường trọng lực vệ tinh với số liệu đo tàu trực tiếp có khả năng triệt tiêu sai số hệ thống ở mức độ nào?
  • $H_1$: Việc kết hợp dữ liệu băng tần Ka từ SARAL/AltiKa và dữ liệu radar giao thoa từ CryoSat-2 ở chế độ GM sẽ làm tăng đáng kể độ phân giải không gian và độ chính xác của trường dị thường độ cao geoid.
  • $H_2$: Phương pháp Collocation bình phương nhỏ nhất (Least Squares Collocation - LSC) tối ưu hơn biến đổi Fourier nhanh (FFT) trong việc kiểm soát hiệu ứng bờ và xử lý số liệu đo cao bất quy tắc tại vùng nước nông.
  • $H_3$: Tích hợp đồng bộ 5 giải pháp công nghệ sẽ nâng cao độ chính xác xác định dị thường trọng lực biển tại Vịnh Bắc Bộ lên mức tốt hơn $\pm 4.0\text{ mGal}$ (đạt thực tế $\pm 2.87\text{ mGal}$), vượt trội so với các mô hình toàn cầu đương thời.

Phạm vi thực nghiệm của luận án bao quát toàn bộ vùng biển Vịnh Bắc Bộ thuộc chủ quyền Việt Nam, sử dụng chuỗi dữ liệu gồm 105 chu kỳ của vệ tinh CryoSat-2/GM, 54 chu kỳ SARAL-DP/AltiKa/GM, 818 điểm GPS - Thủy chuẩn hạng I-II-III, 31 trạm nghiệm triều ven bờ và mạng lưới 2.011 điểm đo trọng lực biển trực tiếp bằng tàu phục vụ hiệu chỉnh và kiểm chứng độc lập.

Literature Review và Positioning

Lịch sử phát triển của kỹ thuật đo cao vệ tinh phục vụ trắc địa vật lý biển đã chứng kiến những bước tiến vượt bậc từ các nghiên cứu tiên phong của Sandwell & Smith (1997), Andersen & Knudsen (1998, 2010), Hwang & Parsons (1995) đến các mô hình trường trọng lực toàn cầu độ phân giải cao như DNSC08GRA, DTU10GRA, DTU13GRA, DTU15GRA và DTU17GRA. Các học giả quốc tế đã áp dụng nhiều công cụ toán học khác nhau: Cheinway Hwang và Barry Parsons (1995) sử dụng phương pháp Collocation bình phương nhỏ nhất kết hợp dữ liệu Geosat/ERM, ERS-1 và TOPEX/Poseidon để xác định dị thường trọng lực tại sống núi Reykjanes với độ chính xác $\pm 5.0\text{ mGal}$; Sandwell và Smith (1997) ứng dụng kỹ thuật biến đổi Fourier nhanh (FFT) trên dữ liệu Geosat/GM và ERS-1/GM đạt độ chính xác toàn cầu $\pm 4.0\text{ mGal}$ đến $\pm 7.0\text{ mGal}$; Andersen và cộng sự (2010, 2017) lần lượt nâng cấp các mô hình DTU bằng cách bổ sung dữ liệu CryoSat-2 và SARAL/AltiKa, đạt độ lệch chuẩn $\pm 3.82\text{ mGal}$ tại Bắc Đại Tây Dương và $\pm 3.23\text{ mGal}$ tại Bắc Cực.

Tuy nhiên, tồn tại một tranh luận học thuật sâu sắc giữa hai trường phái tính toán: một bên ủng hộ phương pháp biến đổi Fourier nhanh (FFT/Stokes ngược) nhờ tốc độ xử lý nhanh trên lưới đều (Andersen & Knudsen, 1998; Sandwell & Smith, 2014), trong khi bên kia chứng minh phương pháp Collocation bình phương nhỏ nhất (LSC) (Moritz, 1980; Hwang, 2002) có tính ưu việt vượt trội trong việc tiếp nhận dữ liệu quan sát không đồng nhất, không phân bố trên lưới đều, đồng thời cho phép tích hợp trực tiếp ma trận hiệp phương sai sai số đo.

Tại khu vực Biển Đông và thềm lục địa Việt Nam, các nghiên cứu quốc tế lẫn trong nước trước đây bộc lộ những khoảng trống lớn:

  • Shengjun Zhang và cộng sự (2017) khi đánh giá các mô hình trọng lực vệ tinh trên Biển Đông đã chỉ ra sai số trung phương (RMSE) lên tới $\pm 5.986\text{ mGal}$; Xiaoli Deng và cộng sự (2001) ghi nhận sai số tại eo biển Đài Loan lên đến $\pm 13.9\text{ mGal}$ do ảnh hưởng của nhiễu bờ biển;
  • Ở trong nước, công trình của Nguyễn Như Trung (2003) khai thác dữ liệu Geosat và ERS-1 chỉ đạt độ chính xác khoảng $\pm 9.0\text{ mGal}$; Bùi Công Quế và cộng sự (2008) kết hợp dữ liệu tàu Gagarinsky, Attalante với dữ liệu đo cao quốc tế cho độ chính xác $\pm 8.0\text{ mGal}$;
  • Luận án của Nguyễn Văn Sáng (2012) 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 (10 chu kỳ) nhưng độ phân giải lưới còn thô ($10' \times 10'$) và độ chính xác dừng lại ở mức $\pm 6.3\text{ mGal}$;
  • Nghiên cứu của Vũ Văn Trí (2017) sử dụng CryoSat-2/GM xác định trọng lực khu vực Quần đảo Trường Sa đạt $\pm 0.67\text{ mGal}$, song khu vực nghiên cứu là vùng biển sâu ngoài khơi xa, nơi các số hiệu chỉnh khí quyển và hải dương học ổn định hơn rất nhiều so với vùng ven bờ.

Luận án của Phạm Văn Tuyên định vị chính xác vào khoảng trống học thuật: chưa có công trình nào xây dựng quy trình tối ưu hóa đa giải pháp chuyên biệt cho vùng biển nông, ven bờ và chịu tương tác thủy triều phức tạp như Vịnh Bắc Bộ. Bằng việc kết hợp các cảm biến đo cao tiên tiến nhất (băng tần Ka và SAR/SIRAL) với kỹ thuật LSC và mô hình hóa cục bộ, nghiên cứu đã tạo ra bước nhảy vọt về độ chính xác xác định trường trọng lực biển tại Việt Nam.

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

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

Nghiên cứu đóng góp quan trọng vào sự phát triển của lý thuyết Trắc địa vật lý hiện đại (Physical Geodesy), cụ thể là lý thuyết xác định trường trọng trường Trái Đất từ các bài toán giá trị biên phi tuyến tính và bán tuyến tính (Molodensky boundary-value problem) và lý thuyết Collocation bình phương nhỏ nhất do Helmut Moritz khởi xướng.

Luận án đã mở rộng mô hình thống kê ước lượng tối ưu trong trắc địa vật lý biển thông qua các luận điểm khoa học có căn cứ thực nghiệm:

  1. Xác lập mối quan hệ giải tích chặt chẽ giữa phần dư thế trọng trường $T$, phần dư độ cao geoid $N_{\text{res}}$ (hoặc $\zeta_{\text{res}}$) và phần dư dị thường trọng lực $\Delta g_{\text{res}}$ theo toán tử vi phân hướng kính: $$\Delta g_{\text{res}} = -\frac{\partial T}{\partial r} - \frac{2}{r}T$$
  2. Xây dựng hàm hiệp phương sai thực nghiệm và làm khớp tối ưu với mô hình hàm hiệp phương sai lý thuyết Hirvonen/Tscherning-Rapp cục bộ, phản ánh đúng đặc trưng thống kê của trường trọng lực khu vực Vịnh Bắc Bộ.
  3. Hoàn thiện cơ sở toán học cho kỹ thuật "Loại bỏ - Phục hồi" (Remove-Compute-Restore - RCR), trong đó trường trọng lực quan sát được phân tách thành 3 thành phần: bước sóng dài (từ GGM), bước sóng trung bình/ngắn (từ phần dư đo cao vệ tinh) và ảnh hưởng địa hình đáy biển/đo tàu trực tiếp.
       [Độ cao mặt biển vệ tinh (SSHcorr)]
[Mô hình GGM EIGEN-6C4]      [Mô hình MDT DTU15MDT]
(Triệt tiêu bước sóng dài)   (Loại bỏ hải dương động học)
         [Phần dư độ cao Geoid (Δh_res)]
      [Collocation bình phương nhỏ nhất (LSC)]
         (EMPCOV -> COVFIT -> GEOCOL)
      [Phần dư dị thường trọng lực (Δg_res)]
 [Làm khớp với 2.011 điểm đo tàu (ZLS Dynamic Meter)]
[Dị thường trọng lực biển hoàn thiện (Grid 1'x1', ±2.87 mGal)]

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

Khung phân tích của luận án tích hợp đồng bộ 5 giải pháp kỹ thuật có tính liên hoàn, giải quyết triệt để các nguồn sai số tích lũy:

  • Tích hợp cảm biến đo cao không đồng nhất: Dung hợp 105 chu kỳ dữ liệu đo cao CryoSat-2 (sử dụng radar giao thoa khẩu độ tổng hợp SIRAL) và 54 chu kỳ SARAL/AltiKa (sử dụng băng tần Ka 35.75 GHz với xung 40 Hz, đường kính dấu chân nhỏ hơn 2 lần so với băng tần Ku thông thường), tạo mạng lưới vết quét cực dày với khoảng cách giữa các tuyến dưới 2-4 km.
  • Lọc phổ thế trọng trường cấp cao: Sử dụng mô hình EIGEN-6C4 khai triển đến bậc và thứ 2190, kết hợp số liệu từ các vệ tinh chuyên dụng GOCE và GRACE, đảm bảo độ chính xác của trường chuẩn bước sóng dài tốt hơn hẳn EGM2008 và SGG-UGM-1.
  • Tách biệt mặt biển trung bình động học: Ứng dụng mô hình DTU15MDT độ phân giải $1' \times 1'$ để chuyển đổi chính xác từ mặt biển trung bình (MSS) sang mặt Geoid, loại bỏ hoàn toàn các thành phần động lực hải dương không đổi theo thời gian.
  • Xử lý số liệu ngẫu nhiên phi cách lưới bằng LSC: Sử dụng mô đun tính toán chuyên dụng trong bộ phần mềm GRAVSOFT (EMPCOV, COVFIT, GEOCOL) để thực hiện phép biến đổi tích phân không gian mà không làm suy giảm độ chính xác do nội suy lưới trung gian như phương pháp FFT.
  • Làm khớp thích nghi với dữ liệu đo trực tiếp: Thiết lập thuật toán hiệu chỉnh phần dư dị thường trọng lực vệ tinh bằng trường sai số nội suy từ 2.011 điểm đo trọng lực biển thực tế bằng máy ZLS Dynamic Meter D06.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành dựa trên thế giới quan thực chứng luận (positivism) và phương pháp luận định lượng nghiêm ngặt của khoa học trắc địa - địa vật lý. Thiết kế nghiên cứu đa cấp độ (multi-level research design) kết hợp giữa mô hình hóa vật lý không gian, phân tích chuỗi thời gian của các xung radar đo cao và kỹ thuật bình sai thống kê không gian.

Quy trình lựa chọn mẫu dữ liệu được kiểm soát chặt chẽ theo các tiêu chuẩn quốc tế:

  • Dữ liệu đo cao vệ tinh được trích xuất từ cơ sở dữ liệu phân phối bởi AVISO, CNES và RADS;
  • Mạng lưới kiểm chuẩn trên đất liền bao gồm 818 điểm GPS trùng khớp với các mốc thủy chuẩn hình học hạng I, II, III phân bố đều trên lãnh thổ Việt Nam;
  • Dữ liệu nghiệm triều gồm chuỗi quan trắc liên tục nhiều năm tại 31 trạm hải văn dọc bờ biển;
  • Mạng lưới kiểm chứng biển gồm 2.011 điểm trích xuất từ 58.989 điểm đo trọng lực trực tiếp của Liên đoàn Địa chất biển (thuộc Dự án Quyết định 47/2006/QĐ-TTg), đo bằng máy đo trọng lực biển cao cấp ZLS Dynamic Gravity Meter với độ chính xác thiết bị đạt $\pm 1.0\text{ mGal}$.

Quy trình nghiên cứu rigorous

Quy trình xử lý số liệu đo cao vệ tinh tuân thủ các bước hiệu chỉnh vật lý và hình học khắt khe nhằm thu được $\text{SSH}{\text{corr}}$ theo công thức: $$\text{SSH}{\text{corr}} = H - \left( d + \Delta d_{\text{dry}} + \Delta d_{\text{wet}} + \Delta d_{\text{iono}} + \Delta d_{\text{SSB}} + \Delta d_{\text{tide}} + \Delta d_{\text{atm}} \right)$$

  1. Hiệu chỉnh khí quyển: Số hiệu chỉnh tầng đối lưu khô ($\Delta d_{\text{dry}}$: $2.0 \div 2.4\text{ m}$) và tầng đối lưu ướt ($\Delta d_{\text{wet}}$: $5 \div 35\text{ cm}$) tính toán từ mô hình khí tượng áp suất bề mặt của Trung tâm Dự báo Thời tiết Hạn vừa Châu Âu (ECMWF); số hiệu chỉnh tầng điện ly ($\Delta d_{\text{iono}}$: $2 \div 20\text{ cm}$) tính từ mô hình điện tử toàn cầu GIM hoặc đo song tần Ku/C.
  2. Hiệu chỉnh bề mặt và động lực biển: Hiệu chỉnh độ lệch trạng thái biển ($\Delta d_{\text{SSB}}$ - Sea State Bias) dựa trên độ cao sóng hữu hiệu (SWH) và tốc độ gió biển; hiệu chỉnh thủy triều đại dương, thủy triều Trái Đất rắn và triều cực từ các mô hình triều độ chính xác cao FES2014/GOT4.8.
  3. Bình sai mạng lưới giao cắt (Crossover Adjustment): Xây dựng thuật toán phát hiện và định vị các điểm giao nhau giữa các cung thăng (ascending tracks) và cung giáng (descending tracks) của các chu kỳ vệ tinh; thiết lập phương trình số hiệu chỉnh để loại bỏ độ cao động học biến đổi theo thời gian $h_t$ và sai số quỹ đạo còn sót lại ($1 \div 2\text{ cm}$).

Data và phân tích

Phân tích thống kê và xử lý toán học được thực hiện trên các nền tảng phần mềm chuyên dụng: gói phần mềm trắc địa vật lý quốc tế GRAVSOFT (với các chương trình GEOIP, EMPCOV, COVFIT, GEOCOL) và phần mềm Geomat2015 do nhóm nghiên cứu phát triển.

Quy trình phân tích dữ liệu thực nghiệm được tổ chức theo 3 phương án đối sánh chặt chẽ:

  • Phương án 1: Chỉ sử dụng số liệu đo cao của vệ tinh CryoSat-2 ở chế độ GM (105 chu kỳ);
  • Phương án 2: Chỉ sử dụng số liệu đo cao của vệ tinh SARAL-DP/AltiKa ở chế độ GM (54 chu kỳ);
  • Phương án 3: Kết hợp đồng thời 105 chu kỳ CryoSat-2/GM và 54 chu kỳ SARAL-DP/AltiKa/GM.

Các tham số thống kê của hàm hiệp phương sai lý thuyết được xác định bằng phương pháp bình phương nhỏ nhất phi tuyến tính, bảo đảm sự phù hợp tối đa giữa đường cong hiệp phương sai giải tích và các giá trị hiệp phương sai thực nghiệm tính từ các cặp điểm khoảng cách $0 \div 110\text{ km}$. Sau bước Collocation xác định $\Delta g_{\text{res}}$, trường trọng lực được làm khớp với 2.011 điểm đo trực tiếp bằng cách mô hình hóa phần dư $\delta g = \Delta g_{\text{ship}} - \Delta g_{\text{alt}}$ cũng thông qua thuật toán LSC.

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

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

  1. Hiệu năng vượt trội của mô hình EIGEN-6C4 tại Việt Nam: Khi đánh giá sai số độ cao geoid/dị thường độ cao tại 818 điểm GPS - Thủy chuẩn trên toàn quốc, mô hình EIGEN-6C4 đạt sai số trung phương nhỏ nhất ($\pm 0.218\text{ m}$), vượt trội so với EGM2008 ($\pm 0.245\text{ m}$), GECO ($\pm 0.231\text{ m}$) và SGG-UGM-1 ($\pm 0.229\text{ m}$). Điều này chứng minh giá trị quyết định của dữ liệu GOCE trong việc tái cấu trúc phổ trọng trường bước sóng dài và trung bình tại khu vực Đông Nam Á.
  2. Độ tin cậy của mô hình DTU15MDT tại vùng biển ven bờ: Khảo sát độ lệch tại 31 trạm nghiệm triều ven biển chứng minh DTU15MDT có độ lệch chuẩn nhỏ nhất ($\pm 0.082\text{ m}$), đảm bảo khả năng tách biệt chính xác mặt Geoid khỏi mực nước biển trung bình động học mà không làm méo mó tín hiệu trọng lực.
  3. Ưu thế tuyệt đối của phương án dung hợp đa vệ tinh CryoSat-2 và SARAL/AltiKa:
    • Phương án 1 (CryoSat-2): Sai số trung phương phần dư dị thường trọng lực trước làm khớp đạt $\pm 4.32\text{ mGal}$, sau làm khớp đạt $\pm 3.12\text{ mGal}$;
    • Phương án 2 (SARAL/AltiKa): Sai số trước làm khớp đạt $\pm 4.18\text{ mGal}$, sau làm khớp đạt $\pm 3.05\text{ mGal}$;
    • Phương án 3 (CryoSat-2 + SARAL/AltiKa): Mật độ vết đo tăng gấp đôi, sai số trước làm khớp đạt $\pm 3.76\text{ mGal}$, sau khi làm khớp với 2.011 điểm đo trực tiếp, sai số giảm xuống còn $\pm 2.87\text{ mGal}$.
  4. Vượt trội so với các mô hình trọng lực biển toàn cầu: Đánh giá kiểm chứng độc lập trên vùng biển Vịnh Bắc Bộ cho thấy kết quả của luận án ($\pm 2.87\text{ mGal}$) có độ chính xác cao hơn rõ rệt so với toàn bộ các mô hình toàn cầu hiện hành: DTU10GRA ($\pm 5.12\text{ mGal}$), DTU13GRA ($\pm 4.68\text{ mGal}$), DTU15GRA ($\pm 4.21\text{ mGal}$) và DTU17GRA ($\pm 3.89\text{ mGal}$).
Mô hình / Phương án nghiên cứu Nguồn dữ liệu tích hợp Phương pháp tính toán Độ chính xác (RMSE tại Vịnh Bắc Bộ)
DTU10GRA (Andersen, 2010) Geosat, ERS-1, Envisat FFT $\pm 5.12\text{ mGal}$
DTU13GRA (Andersen et al., 2013) Bổ sung CryoSat-2, Jason-1/GM FFT + EGM2008 $\pm 4.68\text{ mGal}$
DTU15GRA (Andersen & Knudsen, 2016) CryoSat-2 (5 năm), Jason-2 FFT $\pm 4.21\text{ mGal}$
DTU17GRA (Andersen & Knudsen, 2017) CryoSat-2, Jason-1, SARAL/AltiKa FFT $\pm 3.89\text{ mGal}$
Luận án (Phương án 1) CryoSat-2/GM (105 chu kỳ) LSC + EIGEN-6C4 + Khớp đo tàu $\pm 3.12\text{ mGal}$
Luận án (Phương án 2) SARAL-DP/AltiKa/GM (54 chu kỳ) LSC + EIGEN-6C4 + Khớp đo tàu $\pm 3.05\text{ mGal}$
Luận án (Phương án 3 - Tối ưu) CryoSat-2/GM + SARAL-DP/AltiKa LSC + EIGEN-6C4 + Khớp đo tàu $\pm 2.87\text{ mGal}$

Implications đa chiều

  • Về mặt lý thuyết trắc địa: Khẳng định sự ưu việt của phương pháp Collocation bình phương nhỏ nhất trong việc tái cấu trúc các trường thế địa vật lý tại các bồn trũng trầm tích ven bờ có phân bố dữ liệu dị hướng.
  • Về mặt phương pháp luận: Thiết lập quy trình 5 bước chuẩn hóa có thể chuyển giao và áp dụng trực tiếp cho các vùng biển phức tạp khác như Vịnh Thái Lan, vùng biển quần đảo Hoàng Sa và Trường Sa.
  • Về ứng dụng thực tiễn và kinh tế: Cung cấp lưới dị thường trọng lực độ phân giải $1' \times 1'$ phục vụ giải đoán cấu trúc móng địa chất bồn trũng Sông Hồng, định vị các cấu trúc triển vọng chứa dầu khí, và mô hình hóa địa hình đáy biển sâu ở những nơi tàu đo đạc chưa thể tiếp cận.
  • Về quản lý và an ninh chủ quyền: Dữ liệu trọng lực độ chính xác cao là tham số đầu vào tối quan trọng để thiết lập mặt Geoid địa phương, chuẩn hóa hệ quy chiếu độ cao hải đảo thống nhất với đất liền, và hỗ trợ hiệu chỉnh đường đạn dẫn đường quán tính cho các phương tiện quân sự bảo vệ chủ quyền biển đảo.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, công trình vẫn tồn tại một số giới hạn học thuật khách quan:

  1. Dải dải ven bờ (Coastal buffer strip): Trong phạm vi $5 \div 10\text{ km}$ sát bờ biển, độ chính xác của số liệu đo cao vẫn bị ảnh hưởng bởi hiện tượng dạng sóng phản xạ bị méo do giao thoa đất liền (land contamination), các thuật toán retracking tiêu chuẩn chưa thể phục hồi hoàn toàn.
  2. Phân bố không gian của dữ liệu kiểm chứng tàu biển: 2.011 điểm đo trực tiếp chủ yếu tập trung theo các tuyến khảo sát địa chất sẵn có, chưa phủ trùm đồng đều 100% diện tích các vùng nước cực nông sát cửa sông.
  3. Mô hình biến động thời gian của MDT: Việc sử dụng mô hình MDT trung bình tĩnh DTU15MDT chưa phản ánh đầy đủ các biến động dòng chảy theo mùa quy mô nhỏ đặc thù của chế độ gió mùa Vịnh Bắc Bộ.

Chương trình nghiên cứu tiếp nối trong tương lai cần tập trung vào các định hướng:

  • Ứng dụng các thuật toán hiệu chỉnh dạng sóng ven bờ tiên tiến (Advanced Waveform Retracking như ALES, SAMOSA) chuyên biệt cho dữ liệu đo cao radar vùng ven bờ;
  • Tích hợp dữ liệu đo cao giao thoa dải quét rộng từ vệ tinh SWOT (Surface Water and Ocean Topography) để thu được trường độ cao mặt biển 2D với độ phân giải siêu cao (sub-kilometer);
  • Nghiên cứu mô hình hóa hàm hiệp phương sai phi tĩnh (non-stationary covariance functions) trong phương pháp LSC nhằm thích ứng hoàn hảo với các biến thiên độ sâu địa hình đáy biển đột ngột.

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

Kết quả nghiên cứu của luận án tạo ra chuỗi giá trị tác động sâu rộng:

  • Ảnh hưởng học thuật: Bài báo phối hợp với GS.TS Ole Baltazar Andersen (Đại học Kỹ thuật Đan Mạch - DTU) công bố trên tạp chí quốc tế ISI-Q1 (Advances in Space Research) đã khẳng định vị thế khoa học của trắc địa biển Việt Nam trên diễn đàn quốc tế, mở ra tiềm năng trích dẫn cao trong cộng đồng nghiên cứu trắc địa vật lý và hải dương học vệ tinh toàn cầu.
  • Thúc đẩy ngành công nghiệp dầu khí và khoáng sản biển: Cung cấp cơ sở dữ liệu trường trọng lực nền tảng cho Tập đoàn Dầu khí Việt Nam (Petrovietnam) và Liên đoàn Địa chất biển trong việc xây dựng bản đồ dị thường Bouguer tỷ lệ 1:200.000, giảm thiểu rủi ro và chi phí hàng triệu USD trong các chiến dịch khoan thăm dò đáy biển.
  • Hoàn thiện thể chế và quy phạm kỹ thuật: Cung cấp luận cứ khoa học để Bộ Tài nguyên và Môi trường ban hành các quy định, tiêu chuẩn kỹ thuật quốc gia về ứng dụng dữ liệu viễn thám và đo cao vệ tinh trong điều tra cơ bản tài nguyên - môi trường biển.
  • An ninh quốc phòng: Hoàn thiện mô hình trọng trường biển phục vụ tính toán quỹ đạo chuyển động chính xác của các khí tài quân sự tầm xa và hệ thống dẫn đường ngầm của hải quân trên vùng biển chiến lược Vịnh Bắc Bộ.

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

  • Nghiên cứu sinh và học giả chuyên ngành Trắc địa - Địa vật lý: Tiếp cận được một khung phương pháp luận hoàn chỉnh, hệ thống thuật toán và mã nguồn tính toán chuyển đổi trường thế (Geomat2015, GRAVSOFT workflows) đã được kiểm chứng thực nghiệm.
  • Các nhà khoa học Hải dương học: Sử dụng mô hình Geoid và trường trọng lực chính xác cao để nâng cao độ chính xác xác định dòng chảy hoàn lưu mặt biển, dòng hải lưu mật độ và động lực học vịnh biển.
  • Các tổ chức khảo sát công trình biển: Các đơn vị tư vấn thiết kế công trình dầu khí, điện gió ngoài khơi (offshore wind farms) tại Vịnh Bắc Bộ hưởng lợi trực tiếp từ lưới dữ liệu trọng lực chuẩn để tính toán thế trọng trường và tải trọng nền móng công trình biển.
  • Cơ quan quản lý nhà nước về biển đảo: Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam, Cục Biển và Hải đảo Việt Nam có nguồn dữ liệu tin cậy phục vụ công tác quản lý không gian biển, phân định ranh giới lãnh hải và thềm lục địa theo Luật Biển Việt Nam và UNCLOS 1982.

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?
Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng và thích ứng hóa thành công lý thuyết Collocation bình phương nhỏ nhất (Least Squares Collocation - LSC của Helmut Moritz) vào môi trường biển nông, bán khép kín. Bằng việc xây dựng ma trận hiệp phương sai thực nghiệm cục bộ phản ánh đặc tính phân bố bất đẳng hướng của vùng biển Vịnh Bắc Bộ và tích hợp thuật toán làm khớp sai số đa bậc, nghiên cứu đã giải quyết được hạn chế cố hữu của các phương pháp tích phân Stokes truyền thống vốn đòi hỏi dữ liệu phủ liên tục trên toàn cầu.

2. Điểm cải tiến phương pháp luận cốt lõi khi so sánh với các công trình quốc tế tiêu biểu?
Trả lời: So với phương pháp biến đổi Fourier nhanh (FFT) trong các mô hình toàn cầu DTU (Andersen et al., 2010, 2017) và mô hình Sandwell V23.1 vốn áp đặt lưới chữ nhật đều và dễ phát sinh hiệu ứng biên (edge effect) tại vùng ven bờ, luận án đã ứng dụng LSC trực tiếp từ tập hợp điểm đo dọc vết quét vệ tinh bất quy tắc. Đồng thời, luận án dung hợp đồng thời hai loại cảm biến có nguyên lý đo khác biệt (radar xung thông thường dải Ka của SARAL/AltiKa và radar khẩu độ tổng hợp dải Ku của CryoSat-2), tạo mật độ dữ liệu vượt trội so với các công trình chỉ khai thác một loại vệ tinh của Sandwell & Smith (2014) hay Nguyễn Văn Sáng (2012).

3. Phát hiện bất ngờ nhất trong quá trình tính toán thực nghiệm là gì?
Trả lời: Phát hiện bất ngờ nhất là mô hình thế trọng trường toàn cầu thế hệ mới SGG-UGM-1 và mô hình GECO (kết hợp dữ liệu GOCE) dù có bậc khai triển rất cao nhưng khi đánh giá đối chuẩn tại 818 điểm GPS - Thủy chuẩn và 31 trạm nghiệm triều tại Việt Nam lại cho sai số lớn hơn mô hình EIGEN-6C4 ($\pm 0.229\text{ m}$ và $\pm 0.231\text{ m}$ so với $\pm 0.218\text{ m}$). Điều này chứng minh rằng việc gia tăng bậc hàm điều hòa cầu của GGM không đồng nghĩa với việc nâng cao độ chính xác cục bộ nếu trọng số dữ liệu vệ tinh GOCE/GRACE không được tối ưu hóa chuẩn xác cho khu vực kiến tạo Đông Nam Á.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Trả lời: Luận án cung cấp chi tiết toàn bộ quy trình tái lập bao gồm: danh mục các chu kỳ vệ tinh (105 chu kỳ CryoSat-2, 54 chu kỳ SARAL/AltiKa); bảng thông số hình học ellipsoid WGS-84/GRS-80; quy thức tính toán các số hiệu chỉnh khí quyển - hải dương; sơ đồ luồng thuật toán của chương trình Geomat2015; và chuỗi tham số thiết lập cho các mô đun GEOIP, EMPCOV, COVFIT, GEOCOL thuộc bộ phần mềm GRAVSOFT. Mọi nhà nghiên cứu độc lập đều có thể tái tạo chính xác kết quả khi tiếp cận cùng nguồn dữ liệu gốc.

5. Kế hoạch nghiên cứu 10 năm tới được định hình như thế nào?
Trả lời: Luận án mở ra chương trình nghiên cứu 10 năm với lộ trình 3 giai đoạn: (1) 2024-2026: Ứng dụng kỹ thuật ALES Retracking mở rộng vùng xử lý sát bờ từ $10\text{ km}$ xuống $2\text{ km}$ trên toàn dải bờ biển miền Trung; (2) 2027-2029: Tích hợp chuỗi dữ liệu đo cao giao thoa 2D thế hệ mới từ vệ tinh SWOT và Sentinel-6 Michael Freilich nhằm thành lập bản đồ trọng lực biển toàn diện Biển Đông với độ phân giải $0.5' \times 0.5'$; (3) 2030-2034: Xây dựng hệ thống đồng hóa dữ liệu trọng lực vệ tinh - trọng lực hàng không - hải dương học phục vụ giám sát biến động độ cao mặt biển và biến động khối lượng nước biển dưới tác động của biến đổi khí hậu toàn cầu.

Kết luận

  1. Luận án đã thiết lập thành công cơ sở khoa học và quy trình công nghệ 5 bước tối ưu hóa nhằm xác định dị thường trọng lực biển từ số liệu đo cao vệ tinh, giải quyết trọn vẹn bài toán nâng cao độ chính xác tại các vùng biển ven bờ và biển nông.
  2. Đã chứng minh một cách tường minh trên cơ sở thực nghiệm tính ưu việt của việc kết hợp hai nguồn dữ liệu đo cao thế hệ mới ở chế độ Geodetic Mission: 105 chu kỳ CryoSat-2 và 54 chu kỳ SARAL-DP/AltiKa.
  3. Đã đánh giá toàn diện và lựa chọn chuẩn xác mô hình trường thế trọng trường EIGEN-6C4 (đạt độ chính xác $\pm 0.218\text{ m}$ tại 818 mốc GPS - Thủy chuẩn) và mô hình mặt biển trung bình động học DTU15MDT (đạt độ chính xác $\pm 0.082\text{ m}$ tại 31 trạm nghiệm triều) làm các mặt chuẩn tham chiếu tối ưu cho Việt Nam.
  4. Ứng dụng thành công phương pháp Collocation bình phương nhỏ nhất kết hợp thuật toán làm khớp với 2.011 điểm đo trực tiếp bằng máy ZLS Dynamic Meter, nâng cao độ chính xác xác định dị thường trọng lực biển Vịnh Bắc Bộ đạt mức $\pm 2.87\text{ mGal}$ trên lưới phân giải cao $1' \times 1'$, vượt xa mục tiêu đặt ra ($\pm 4.0\text{ mGal}$).
  5. Kết quả của luận án thể hiện bước tiến vượt bậc so với các mô hình trọng lực biển toàn cầu hiện hành (DTU10, DTU13, DTU15, DTU17) trên vùng biển thử nghiệm, đóng góp một bộ dữ liệu điều tra cơ bản vô cùng giá trị cho sự nghiệp phát triển kinh tế biển, thăm dò địa chất dầu khí và củng cố vững chắc quốc phòng - an ninh biển đảo Việt Nam.