Tổng quan về luận án

Sự bùng nổ của công nghệ định vị vệ tinh toàn cầu (GNSS) và kỹ thuật đo cao vệ tinh trong kỷ nguyên trắc địa không gian hiện đại đã đặt ra yêu cầu cấp thiết về việc thiết lập một khung quy chiếu thống nhất, chính xác cao và mang tính động lực. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Trắc địa - Bản đồ (Mã số: 9.03) của tác giả Lương Thanh Thạch, dưới sự hướng dẫn khoa học của PGS. TSKH Hà Minh Hòa tại Viện Khoa học Đo đạc và Bản đồ (2018), tập trung giải quyết bài toán cốt lõi: "Nghiên cứu phương pháp xây dựng và phát triển hệ quy chiếu tọa độ không gian quốc gia".

Tại Việt Nam, Hệ quy chiếu và Hệ thống tọa độ quốc gia VN2000 được thiết lập năm 2000 dựa trên ellipsoid WGS84 định vị cục bộ và Hệ độ cao chuẩn quốc gia Hải Phòng 1972 (HP72). Mặc dù đáp ứng hiệu quả công tác đo vẽ bản đồ địa hình và quản lý đất đai trong giai đoạn đầu, VN2000 bộc lộ ba khoảng trống nghiên cứu (research gaps) mang tính bản chất:

  1. Bản chất tĩnh hai chiều (2D): Tách rời hoàn toàn hệ tọa độ mặt bằng (VN2000) và hệ độ cao chuẩn (HP72), cản trở việc xác định trực tiếp tọa độ không gian 3D $(X, Y, Z)$ và độ cao trắc địa $H$.
  2. Độ chính xác mô hình Geoid/Quasigeoid hạn chế: Mô hình VNGeo-96 trước đây dựa trên EGM96 và 367 điểm GPS/thủy chuẩn hỗn hợp có sai số dị thường độ cao lên tới $\pm 1.15\text{ m}$; mô hình Geoid2010 đạt $\pm 8.2\text{ cm}$; chưa đáp ứng chuẩn xác định độ cao trắc địa thay thế thủy chuẩn hình học.
  3. Liên kết yếu với Khung quy chiếu Trái Đất quốc tế (ITRF): Mạng lưới GPS cấp 0 gồm 69 điểm có thời gian đo ngắn (khoảng 7 giờ), xử lý bằng phần mềm thương mại TRIM NET+ và GPSurvey, độ chính xác vị trí mặt bằng chỉ đạt $\pm 0.11\text{ m}$ và chưa đạt độ chính xác tương đối $10^{-9}$ theo chuẩn IAG (Augath, 1996). 7 tham số chuyển đổi Bursa - Wolf ban hành theo Quyết định 05/2007/QĐ-BTNMT không cho phép chuyển đổi tọa độ không gian từ ITRF sang tọa độ không gian quốc gia với độ chính xác mức centimet.

Câu hỏi và giả thuyết nghiên cứu của luận án được xác định cụ thể:

  • RQ1: Mô hình quasigeoid nào tại Việt Nam đáp ứng đủ độ tin cậy để làm bề mặt quy chiếu chuyển đổi độ cao trong hệ quy chiếu không gian quốc gia (HQCTĐKGQG)?
  • RQ2: Mô hình toán học nào tối ưu hóa việc bình sai ghép nối mạng lưới GNSS quốc gia phủ trùm vào HQCTĐKGQG nhằm chính xác hóa 7 tham số liên kết ITRF?
  • RQ3: Làm thế nào để cập nhật, mở rộng các điểm cơ sở trắc địa mới vào hệ thống mà không phá vỡ tọa độ của mạng lưới gốc đã được công bố?
  • H1: Xây dựng HQCTĐKGQG thành công đòi hỏi mô hình quasigeoid có độ chính xác không thấp hơn $\pm 4\text{ cm}$ kết hợp mạng lưới GNSS phủ trùm liên kết chặt chẽ với ITRF.
  • H2: Thuật toán bình sai truy hồi với phép biến đổi xoay trực giao cho phép bổ sung điểm cơ sở không gian mới từ cơ sở dữ liệu (CSDL) trắc địa mà không cần bình sai lại toàn bộ mạng lưới gốc.

Phạm vi nghiên cứu bao quát dữ liệu mạng lưới trắc địa động lực Miền Bắc, 75 điểm độ cao hạng II, 25 điểm độ cao hạng I Nhà nước, mô hình quasigeoid VIGAC2014 và VIGAC2017. Nghiên cứu tạo bước đột phá khi hoàn thiện cơ sở phương pháp luận toán học - trắc địa vật lý, cho phép nâng độ chính xác tọa độ trắc địa lên gấp 2 lần và thiết lập nền tảng tích hợp hạ tầng dữ liệu không gian quốc gia (NSDI).

Literature Review và Positioning

Lý thuyết trắc địa vật lý và định vị không gian toàn cầu ghi nhận sự dịch chuyển mạnh mẽ từ các hệ quy chiếu cục bộ 2D sang hệ quy chiếu không gian địa tâm 3D/4D gắn với ITRS/ITRF (Boucher, 1993; Dennis, 1996). Trên thế giới, các cường quốc trắc địa đã sớm hoàn thiện hệ quy chiếu không gian quốc gia:

  • Hoa Kỳ: Cơ quan Khảo sát Trắc địa Quốc gia Hoa Kỳ (NGS) phát triển NSRS dựa trên ITRF2000, thay thế NAD83 và NAVD88; tích hợp mô hình USGG2012 ($\pm 3-4\text{ cm}$) tiến tới USGG2022 ($\pm 1-2\text{ cm}$) kết hợp mạng lưới GNSS thụ động HARN và hệ thống trạm CORS vận hành dịch vụ trực tuyến OPUS (Doyle, 2000; Snay & Soler, 2008).
  • Canada: Chuyển đổi từ CGVD28 sang CGVD2013 dựa trên mô hình geoid trọng lực CGG2013 ($\pm 3\text{ cm}$, hướng tới $\pm 1\text{ cm}$), duy trì hệ thống CSRS gắn với ITRF96 (Epoch 1997.0) thông qua 250 trạm hoạt động liên tục CACS và mạng lưới cơ sở CBN cách nhau 200 km cùng dịch vụ CSRS-PPP (Véronneau, 2013).
  • Vương quốc Anh: Cơ quan Đo đạc Quân sự và Dân sự (Ordnance Survey) tích hợp mạng lưới tam giác OSTN với ETRS89 hình thành OS Net (90 trạm CORS), sử dụng mô hình OSGM02 ($\pm 2-4\text{ cm}$) và hoàn thiện OSGM15 đạt độ chính xác $\pm 1.0\text{ cm}$ trên đất liền, $\pm 2.0\text{ cm}$ tại Orkney/Shetland và $\pm 2.4\text{ cm}$ tại Bắc Ireland (Iliffe, 2003; Greaves, 2006).
  • Australia: Triển khai GDA2020 trên nền tảng ITRF2014 và mô hình quasigeoid AUSGeoid09/AUSGeoid2020 ($\pm 3\text{ cm}$) xây dựng theo lý thuyết Molodensky từ 1.4 triệu dị thường trọng lực, mô hình DEM 9"x9" và 911 điểm GPS/thủy chuẩn, phục vụ hạ tầng ATRF và dịch vụ AUSPOS (Featherstone & Johnston, 2011; Janssen, 2017).
  • Đức: Phát triển ETRS89/DREF91/2016 kết hợp mô hình quasigeoid hỗn hợp GCG2016 đạt độ chính xác cao hơn $\pm 1\text{ cm}$.
                             [ ITRS / ITRF ]
                     (IGS Tracking Stations, Epoch t)
        [ Quốc tế / Khu vực ]               [ Việt Nam: Nghiên cứu ]
       • ETRS89 (Châu Âu / EPN)           • HQCTĐKGQG (VN2000-3D)
       • NAD83 / NSRS (Hoa Kỳ)             • Ellipsoid WGS84 định vị cục bộ
       • CSRS (Canada)                     • Mô hình Quasigeoid VIGAC2017 (±6.4cm)
       • GDA2020 / ATRF (Úc)               • 7 tham số Bursa-Wolf chuẩn hóa
                 [ Mạng lưới GNSS quốc gia phủ trùm ]
                 (Passive Control: 20-100 km; CORS)
       [ Giải pháp xử lý dữ liệu: Bình sai ghép nối khối ]
          [ Bình sai truy hồi với phép biến đổi xoay ]

Tại Việt Nam, các nghiên cứu trước đây của Hà Minh Hòa (2012a, 2013a, 2014a, 2015c) đã đặt nền móng về mặt thế trọng trường $W_0$ của mặt quasigeoid cục bộ Hòn Dấu, phần mềm ECME-GPS, GNSS Net A&C và VNAOLG-VIGAC. Tuy nhiên, tồn tại cuộc tranh luận học thuật giữa hai trường phái:

  1. Trường phái bình sai tự do gộp chung: Bình sai đồng thời toàn bộ trị đo mặt bằng, độ cao và GNSS. Phương pháp này đòi hỏi khối lượng tính toán khổng lồ, dễ gây biến dạng cục bộ các mạng lưới truyền thống ổn định và làm thay đổi toàn bộ tọa độ điểm gốc lịch sử.
  2. Trường phái bình sai riêng rẽ và ghép nối khối (Luận án kế thừa và phát triển): Tách biệt mạng lưới GNSS quốc gia trong khung ITRF và mạng lưới tọa độ không gian quốc gia, sau đó thiết lập mô hình toán bình sai ghép nối thông qua các điểm song trùng và hiệu chỉnh tọa độ.

Luận án định vị rõ đóng góp: Cung cấp giải pháp toán học hoàn chỉnh kết nối giữa mạng lưới GNSS phủ trùm với mạng lưới thiên văn - trắc địa cổ điển trong một khung quy chiếu 3D thống nhất, đồng thời đề xuất thuật toán truy hồi mở rộng mạng lưới độc lập với quy mô tính toán của toàn hệ thống.

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

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

Luận án tạo ra những bước tiến lý thuyết vững chắc trong ngành Trắc địa vật lý và Xử lý số liệu trắc địa:

  1. Mở rộng lý thuyết định vị không gian 3D trên ellipsoid quy chiếu cục bộ: Luận án chứng minh một cách tường minh mối quan hệ vật lý - hình học giữa ba bề mặt: bề mặt vật lý Trái Đất, mặt quasigeoid quốc gia và mặt ellipsoid WGS84 quy chiếu cục bộ. Tọa độ không gian $(X, Y, Z)$ liên hệ chặt chẽ với tọa độ trắc địa $(B, L, H)$, trong đó độ cao trắc địa $H$ được chuẩn hóa qua phương trình: $$H = H^\gamma + \zeta$$ với $H^\gamma$ là độ cao chuẩn quốc gia (hệ HP72) và $\zeta$ là dị thường độ cao xác định từ mô hình quasigeoid quốc gia độ chính xác cao.
  2. Chuẩn hóa mô hình chuyển đổi không gian 7 tham số Bursa - Wolf: Thiết lập hệ phương trình chuyển đổi từ tọa độ không gian $(X, Y, Z){\text{ITRF}}$ trên ellipsoid WGS84 quốc tế sang $(X, Y, Z){\text{VN2000-3D}}$ trên ellipsoid WGS84 quốc gia: $$\begin{bmatrix} X \ Y \ Z \end{bmatrix}{\text{VN2000-3D}} = \begin{bmatrix} dX_0 \ dY_0 \ dZ_0 \end{bmatrix} + (1 + m) \begin{bmatrix} 1 & \omega_Z & -\omega_Y \ -\omega_Z & 1 & \omega_X \ \omega_Y & -\omega_X & 1 \end{bmatrix} \begin{bmatrix} X \ Y \ Z \end{bmatrix}{\text{ITRF}}$$ Mô hình loại trừ sai số dịch chuyển gốc tọa độ $(dX_0, dY_0, dZ_0)$, các góc xoay vi phân $(\omega_X, \omega_Y, \omega_Z)$ và hệ số tỷ lệ $m$, bảo đảm tính tương thích chặt chẽ với khung ITRF hiện hành.
  3. Phát triển lý thuyết bình sai truy hồi với phép biến đổi xoay (Sequential Adjustment with Orthogonal Rotation Transformation): Chứng minh rằng việc bổ sung các trị đo GNSS mới hoặc các trạm trắc địa mới chỉ làm biến đổi cục bộ ma trận nghịch đảo tương quan thông qua biến đổi xoay Givens/Householder, triệt tiêu nhu cầu bình sai lại toàn bộ ma trận phương trình chuẩn cấp cao của mạng lưới gốc.
                           [ LÝ THUYẾT NỀN TẢNG ]
              Molodensky Potential Theory & ITRS/ITRF Physics
       [ Khung Chuyển Đổi 7 Tham Số ]          [ Lý Thuyết Bình Sai Khối ]
         (Mô hình Bursa - Wolf 3D)           (Tách riêng & Ghép nối qua điểm song trùng)
                  [ CƠ CHẾ TOÁN HỌC TRUY HỒI ĐỘT PHÁ ]
             (Bình sai truy hồi bằng phép biến đổi xoay trực giao)
                [ HỆ QUY CHIẾU TỌA ĐỘ KHÔNG GIAN QUỐC GIA ]

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: (1) Lý thuyết trường thế trọng trường Trái Đất (Molodensky Potential Theory), (2) Lý thuyết ma trận và tối ưu hóa bình sai lưới tự do/lưới ràng buộc (Least Squares Adjustment Theory), và (3) Mô hình kiến tạo mảng địa động lực (Plate Tectonics Kinematics).

Phương pháp tiếp cận giải quyết bài toán bình sai khối mạng lưới không gian thông qua 3 giai đoạn phân tích logic:

  • Giai đoạn 1: Bình sai độc lập mạng lưới GNSS quốc gia phủ trùm trong khung quy chiếu ITRF, thu nhận tọa độ không gian $(X, Y, Z){\text{ITRF}}$ và ma trận hiệp phương sai sai số $\mathbf{K}{\mathbf{X}_{\text{ITRF}}}$.
  • Giai đoạn 2: Thiết lập bài toán bình sai ghép nối trên tập hợp các điểm song trùng (các mốc thiên văn - trắc địa hạng I, II và độ cao hạng I, II có đo nối GNSS), xác định đồng thời 7 tham số chuyển đổi tối ưu và tọa độ bình sai trong hệ VN2000-3D.
  • Giai đoạn 3: Ứng dụng thuật toán truy hồi cập nhật trị đo để chèn các điểm cơ sở không gian mới mà không gây nhiễu loạn nghiệm tọa độ của các trạm mốc lịch sử.

Điều kiện biên (Boundary conditions) được xác định rõ ràng: Áp dụng cho các mạng lưới trắc địa quốc gia có mật độ điểm GNSS phủ trùm từ $20 - 100\text{ km/điểm}$, nền địa chất ổn định, thu nhận tối thiểu 5 vệ tinh liên tục và mô hình quasigeoid khu vực đạt độ chính xác sai số trung phương $m_\zeta \le \pm 4\text{ cm}$.

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ủ thế giới quan thực chứng khoa học (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism) trong trắc địa vật lý, sử dụng phương pháp luận định lượng nghiêm ngặt. Thiết kế nghiên cứu kết hợp giữa mô hình hóa toán học giải tích, lập trình thuật toán mô phỏng và thực nghiệm kiểm chứng trên dữ liệu đo thực tế mạng lưới quốc gia.

Thiết kế đa cấp bậc (Multi-level design) bao gồm:

  1. Cấp độ toàn cầu/khu vực: Khung quy chiếu ITRF/IGS, mô hình trọng trường toàn cầu EGM2008, dữ liệu đo cao vệ tinh altimetry.
  2. Cấp độ quốc gia: Mạng lưới trắc địa thiên văn hạng I, II; mạng lưới thủy chuẩn hình học hạng I, II Nhà nước; hệ thống trạm GNSS CORS quốc gia.
  3. Cấp độ cục bộ/địa động lực: Mạng lưới trắc địa động lực phục vụ dự báo tai biến đứt gãy khu vực Miền Bắc Việt Nam.

Tiêu chí lựa chọn điểm trắc địa trong mạng lưới GNSS quốc gia phủ trùm:

  • Thu nhận tín hiệu liên tục không ít hơn 5 vệ tinh GNSS với góc ngưỡng quan sát $\ge 15^\circ$;
  • Tọa lạc trên cấu trúc nền địa chất công trình ổn định, không bị lún sụt cục bộ;
  • Thuận tiện cho việc tiếp cận, bảo quản lâu dài và liên kết trực tiếp với mốc độ cao chuẩn hoặc mốc tọa độ hạng cao hiện hữu.
       [ THIẾT KẾ ĐA CẤP BẬC ]
       [ QUY TRÌNH THỰC NGHIỆM & PHÂN TÍCH ]

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý số liệu được chuẩn hóa thông qua hệ thống công cụ phần mềm chuyên dụng:

  • Xử lý baseline GNSS bằng phần mềm khoa học chuyên dụng kết hợp module ECME-GPS;
  • Bình sai ghép nối khối và mô phỏng thuật toán bằng hệ thống phần mềm GNSS Net A&C và VNAOLG-VIGAC do Viện Khoa học Đo đạc và Bản đồ phát triển;
  • Tính toán mô hình trọng trường và dị thường trọng lực bằng bộ chương trình GRAVSOFT tích hợp công cụ GOCE GUT;
  • Đối soát, kiểm tra chéo (triangulation) giữa 3 nguồn dữ liệu độc lập: Dữ liệu GNSS không gian, dữ liệu thủy chuẩn hình học độ chính xác cao và số liệu trọng lực mặt đất kết hợp mô hình số địa hình DTM.

Quy trình kiểm soát chất lượng dữ liệu tuân thủ nghiêm ngặt chuẩn FGCC (1984, 1988), loại trừ ảnh hưởng của độ lệch dây dọi, chiết quang khí quyển tầng đối lưu và tầng điện ly, chuẩn hóa sai số đồng hồ vệ tinh qua lịch vệ tinh chính xác của IGS.

Data và phân tích

Dữ liệu thực nghiệm của luận án bao gồm:

  • Dữ liệu đánh giá Geoid/Quasigeoid:
    • Tập mẫu 1: 75 điểm độ cao hạng II Nhà nước có đo nối GPS phân bố đồng đều tại khu vực phía Bắc phục vụ kiểm chuẩn mô hình quasigeoid hỗn hợp VIGAC2014.
    • Tập mẫu 2: 25 điểm độ cao hạng I Nhà nước thuộc mạng lưới chuẩn quốc gia phục vụ kiểm chuẩn mô hình quasigeoid quốc gia khởi đầu VIGAC2017.
  • Dữ liệu thực nghiệm mạng lưới GNSS:
    • Hệ thống các điểm song trùng thuộc Mạng lưới cơ sở các điểm địa động lực Miền Bắc (chủ trì bởi Viện Khoa học Đo đạc và Bản đồ).
    • Tọa độ trắc địa trạm gốc, các vector baseline cạnh dài và ma trận tương quan hiệp phương sai đầy đủ.

Phân tích thống kê thực nghiệm áp dụng các kiểm định sai số trung phương trọng số đơn vị ($\mu$), sai số trung phương vị trí điểm ($m_p$), hiệu chỉnh sai số hệ thống có bước nhảy và phân tích phương sai sai số dư. Thuật toán bình sai khối và bình sai truy hồi được lập trình trực tiếp bằng mã nguồn cấu trúc module hóa (ngôn ngữ Pascal/C++ tích hợp trong phụ lục luận án), đảm bảo tính tái lập (replicability) tuyệt đối của các kết quả nghiên cứu.

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

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

Quá trình phân tích thực nghiệm đem lại 4 phát hiện khoa học mang tính đột phá:

  1. Xác thực độ chính xác của các mô hình Quasigeoid tại Việt Nam:

    • Đánh giá trên 75 điểm độ cao hạng II Nhà nước, mô hình quasigeoid hỗn hợp VIGAC2014 cho thấy độ lệch có hệ thống cục bộ trước khi hiệu chỉnh, nhưng sau khi khử sai số hệ thống, sai số trung phương xác định dị thường độ cao đạt mức $\pm 4.3\text{ cm}$.
    • Đánh giá trên 25 điểm độ cao hạng I Nhà nước độc lập, mô hình quasigeoid quốc gia khởi đầu VIGAC2017 đạt sai số trung phương $\pm 6.4\text{ cm}$, vượt trội so với mô hình Geoid2010 ($\pm 8.2\text{ cm}$) và GeoidQS quân sự ($\pm 12.7\text{ cm}$). Mô hình VIGAC2017 bao trùm toàn bộ đất liền, thềm lục địa, vùng biển đảo Hoàng Sa - Trường Sa và được định vị khớp với mực nước biển trung bình trạm nghiệm triều Hòn Dấu.

  2. Nâng cao vượt bậc độ chính xác mạng lưới thông qua bình sai ghép nối:

    • Kết quả bình sai ghép nối mạng lưới GNSS quốc gia phủ trùm vào HQCTĐKGQG (VN2000-3D) chứng minh sai số trung phương vị trí điểm mặt bằng giảm từ mức $\pm 0.11\text{ m}$ (trong mạng lưới thiên văn - trắc địa cũ) xuống mức centimet ($\pm 0.01 - \pm 0.03\text{ m}$). Độ chính xác vị trí không gian được nâng lên gấp 2 lần, đồng thời xác lập bộ 7 tham số Bursa - Wolf chính xác cao có khả năng chuyển đổi trực tiếp giữa ITRF và VN2000-3D.

  3. Tính ưu việt tuyệt đối của thuật toán bình sai truy hồi với phép biến đổi xoay:

    • Thực nghiệm so sánh giữa ba phương án: (a) Bình sai tổng thể toàn mạng lưới, (b) Bình sai một phần mạng lưới, và (c) Bình sai truy hồi bổ sung trị đo làm thay đổi ẩn số. Kết quả khẳng định tọa độ của các điểm bổ sung mới tính theo thuật toán truy hồi biến đổi xoay trùng khít hoàn toàn với kết quả bình sai tổng thể (độ lệch sai số tọa độ $\Delta X, \Delta Y, \Delta Z = 0.000\text{ m}$).
    • Phương pháp này giải phóng hoàn toàn gánh nặng lưu trữ và đảo ma trận cấp lớn, cho phép cập nhật tức thời các trạm CORS hoặc điểm GNSS mới vào Hệ thống thông tin trắc địa quốc gia.

  4. Hiện tượng chênh lệch phi tuyến giữa độ cao trắc địa và độ cao chuẩn:

    • Nghiên cứu phát hiện và lý giải nguyên nhân sai số tích lũy trong các chuyển đổi trước đây: Sự không tương thích giữa hệ thống triều (tide-free system của ITRF/EGM2008 và mean-tide system của mạng lưới độ cao quốc gia) cùng sự thiếu hụt dữ liệu trọng lực chi tiết vùng núi hiểm trở gây ra độ nghiêng hệ thống trong mô hình geoid cục bộ.
Tiêu chí / Thông số VN2000 Cổ điển (Năm 2000) Kết quả Luận án / HQCTĐKGQG Đề xuất Mức độ Cải thiện
Không gian quy chiếu 2D Tĩnh (Mặt bằng $B, L$) tách rời Độ cao $H^\gamma$ 3D Không gian $(X, Y, Z)$ & $(B, L, H)$ thống nhất Tích hợp hoàn chỉnh 3D
Mô hình Geoid / Quasigeoid VNGeo-96 ($m_\zeta = \pm 1.15\text{ m}$) VIGAC2017 ($m_\zeta = \pm 6.4\text{ cm}$) Độ chính xác tăng ~18 lần
Độ chính xác mặt bằng Sai số vị trí lớn nhất $\pm 0.11\text{ m}$ Sai số vị trí bình sai ghép nối $\pm 0.01 - 0.03\text{ m}$ Nâng độ chính xác 3-5 lần
Cơ chế liên kết ITRF 7 tham số gần đúng (QĐ 05/2007) 7 tham số Bursa-Wolf chuẩn hóa qua bình sai khối Triệt tiêu sai số hệ thống
Phương thức cập nhật điểm mới Bình sai lại toàn bộ lưới hoặc gá lắp tự do Bình sai truy hồi với biến đổi xoay trực giao Không đổi tọa độ gốc, sai lệch 0.000m

Implications đa chiều

  • Về mặt lý thuyết: Luận án hoàn thiện lý luận về xây dựng hệ quy chiếu không gian động học tại các quốc gia đang phát triển; giải quyết triệt để sự mâu thuẫn giữa việc bảo tồn giá trị pháp lý của hệ tọa độ lịch sử và yêu cầu hội nhập các khung quy chiếu toàn cầu ITRS/ITRF.
  • Về mặt phương pháp luận: Cung cấp bộ công cụ thuật toán bình sai khối và truy hồi có thể chuyển giao cho các lĩnh vực liên quan như trắc địa công trình ngầm, quan trắc biến dạng công trình thủy điện, cầu lớn và mạng lưới quan trắc trượt lở đất đá.
  • Về mặt thực tiễn và chính sách:
    • Cung cấp cơ sở khoa học trực tiếp để Bộ Tài nguyên và Môi trường hoàn thiện Chiến lược phát triển ngành Đo đạc và Bản đồ Việt Nam theo Quyết định số 33/2008/QĐ-TTg;
    • Làm tiền đề kỹ thuật để triển khai mạng lưới trạm định vị vệ tinh quốc gia VNGEONET (mạng lưới CORS quốc gia) và hiện đại hóa hệ thống độ cao quốc gia;
    • Cho phép thay thế phương pháp đo thủy chuẩn hình học truyền thống tốn kém ở vùng đồi núi hiểm trở bằng công nghệ đo cao GNSS kết hợp mô hình quasigeoid chính xác cao.

Limitations và Future Research

Luận án thẳng thắn chỉ ra các giới hạn kỹ thuật nội tại:

  1. Mật độ trạm kiểm chuẩn mô hình quasigeoid: Số lượng điểm kiểm tra độc lập hạng I (25 điểm) và hạng II (75 điểm) chủ yếu tập trung tại khu vực miền Bắc, chưa bao phủ mật độ dày đặc toàn diện trên địa bàn Tây Nguyên, Nam Bộ và các đảo xa bờ.
  2. Mô hình vận tốc chuyển dịch kiến tạo mảng (Velocity Field): Luận án tập trung vào bài toán bình sai hình học không gian 3D tại thời điểm chuẩn (Epoch), chưa tích hợp đầy đủ mô hình trường vận tốc chuyển động vỏ Trái Đất 4D (với 3 góc xoay Euler đặc trưng) cho từng khối kiến tạo nhỏ tại Việt Nam.
  3. Độ chính xác mô hình Quasigeoid khởi đầu: Mức sai số $\pm 6.4\text{ cm}$ của VIGAC2017 tuy đạt bước tiến lớn nhưng vẫn chưa đạt ngưỡng lý tưởng $\pm 1 - 2\text{ cm}$ như các quốc gia phát triển (USGG2022 của Mỹ hay GCG2016 của Đức).

Chương trình nghiên cứu tiếp nối cần tập trung vào 4 hướng chiến lược:

  • Đo nối bổ sung mạng lưới trọng lực hàng không (Airborne Gravity) và trọng lực biển chi tiết để nâng độ chính xác mô hình quasigeoid quốc gia đạt ngưỡng $\le \pm 2\text{ cm}$;
  • Tích hợp dữ liệu từ mạng lưới trạm GNSS CORS hoạt động liên tục quốc gia nhằm xây dựng mô hình biến dạng hiện đại vỏ Trái Đất dạng động lực học 4D (ITRF-consistent Dynamic Datum);
  • Nghiên cứu mô hình dị thường độ cao động tích hợp ảnh hưởng triều đại dương và khí quyển thời gian thực phục vụ dẫn đường hàng hải và hàng không chính xác cao;
  • Tối ưu hóa kiến trúc CSDL trắc địa quốc gia theo mô hình phân cấp hướng đối tượng, tích hợp dịch vụ xử lý dữ liệu định vị điểm chính xác trực tuyến (VN-PPP) tương tự CSRS-PPP của Canada và AUSPOS của Australia.

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

  • Ảnh hưởng học thuật: Công trình là tài liệu tham khảo mẫu mực trong đào tạo sau đại học chuyên ngành Trắc địa - Bản đồ, Địa vật lý và Kỹ thuật Địa không gian tại Đại học Tài nguyên và Môi trường Hà Nội, Đại học Mỏ - Địa chất, Học viện Kỹ thuật Quân sự và Đại học Bách khoa TP.HCM. Các thuật toán của luận án mở ra hướng nghiên cứu mới về xử lý dữ liệu trắc địa lớn (Big Geodetic Data).
  • Chuyển đổi công nghiệp & đo đạc ứng dụng: Giảm thiểu $40 - 60%$ chi phí đo đạc ngoại nghiệp trong công tác xây dựng lưới khống chế trắc địa công trình giao thông, thủy lợi, xây dựng dân dụng thông qua việc áp dụng công nghệ GNSS/Quasigeoid thay thế thủy chuẩn hình học hạng III, IV.
  • Tác động chính sách & an ninh quốc phòng: Đảm bảo tính thống nhất của hạ tầng dữ liệu địa lý dùng chung phục vụ quy hoạch không gian biển, phân định ranh giới lãnh hải, xác định tọa độ các công trình quốc phòng, dẫn đường chính xác cho các phương tiện không người lái (UAV, USV) và giám sát tai biến địa chất đứt gãy kiến tạo phục vụ phòng chống thiên tai.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Nắm vững phương pháp luận bình sai khối hiện đại, thuật toán biến đổi xoay trong bình sai truy hồi và quy trình đánh giá mô hình geoid/quasigeoid theo chuẩn quốc tế.
  • Các nhà khoa học đầu ngành: Kế thừa khung lý thuyết để phát triển hệ quy chiếu không gian động học 4D và xây dựng mô hình trọng trường độ phân giải siêu cao.
  • Doanh nghiệp khảo sát, thiết kế & R&D trắc địa: Tiếp cận công nghệ đo cao vệ tinh GNSS độ chính xác cao, rút ngắn thời gian thi công lưới trắc địa công trình và nâng cao độ tin cậy của bản đồ địa hình số.
  • Cơ quan quản lý nhà nước (Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam; Cục Bản đồ/Bộ Tổng tham mưu): Sở hữu cơ sở khoa học và công cụ tính toán để chuẩn hóa, ban hành các quy chuẩn kỹ thuật quốc gia về hệ quy chiếu không gian và quản trị hiệu quả CSDL trắc địa 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à việc thiết lập giải pháp toán học kết hợp giữa Mô hình chuyển đổi 7 tham số Bursa - Wolf và Thuật toán bình sai truy hồi với phép biến đổi xoay trực giao. Công trình đã mở rộng lý thuyết bình sai mạng lưới trắc địa truyền thống sang không gian 3 chiều gắn với mô hình thế trọng trường Molodensky, cho phép hợp nhất hai hệ thống độc lập lịch sử (hệ tọa độ mặt bằng VN2000 và hệ độ cao chuẩn HP72) thành một hệ quy chiếu tọa độ không gian 3D thống nhất tương thích hoàn toàn với ITRF mà vẫn giữ nguyên tính ổn định của các điểm mốc gốc.

2. Đổi mới về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế?

So với hệ thống NAD83 của Mỹ (sử dụng bình sai tự do tái hiệu chỉnh toàn bộ qua HARN) và GDA94 của Úc (chuyển đổi cứng nhắc theo lưới biến dạng), luận án chọn phương pháp Bình sai ghép nối khối thông qua các điểm song trùng kết hợp xử lý truy hồi ma trận. Phương pháp này giải quyết được bài toán đặc thù của Việt Nam: nguồn lực hạn chế, mạng lưới mốc thiên văn cổ điển bị hư hao nhưng vẫn cần bảo toàn tính pháp lý của tọa độ VN2000 hiện hành, đồng thời loại bỏ nhu cầu tính toán lại toàn bộ hệ phương trình chuẩn khi mạng lưới mở rộng.

3. Phát hiện thực nghiệm đáng ngạc nhiên nhất từ dữ liệu nghiên cứu là gì?

Đó là tính chất sai số bằng không tuyệt đối ($\Delta X = 0.000\text{ m}, \Delta Y = 0.000\text{ m}, \Delta Z = 0.000\text{ m}$) khi so sánh tọa độ điểm tính bổ sung giữa phương pháp bình sai truy hồi với phép biến đổi xoay và phương pháp bình sai tổng thể toàn mạng lưới. Điều này chứng minh tính đúng đắn tuyệt đối về mặt toán học của thuật toán truy hồi, đập tan định kiến học thuật cho rằng việc bình sai cục bộ/bổ sung luôn làm suy giảm độ chính xác hình học so với bình sai gộp đồng thời.

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 trình bày chi tiết từ hệ thống công thức giải tích toán học, sơ đồ khối thuật toán, ma trận dữ liệu thực nghiệm mẫu của các vector baseline, tọa độ điểm gốc, cho đến việc công bố toàn bộ mã nguồn chương trình tính toán (viết dưới dạng các module chức năng lập trình) trong phần phụ lục, cho phép các nhà nghiên cứu độc lập kiểm chứng và chạy lại toàn bộ quy trình tính 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?

Lộ trình 10 năm định hướng chuyển đổi HQCTĐKGQG từ trạng thái 3D tĩnh sang Hệ quy chiếu không gian động học 4D (Time-dependent Geodetic Datum). Lộ trình tích hợp mạng lưới trạm GNSS CORS quốc gia thu nhận dữ liệu liên tục theo thời gian thực, xây dựng mô hình trường vận tốc vỏ Trái Đất đa mảng cục bộ, hoàn thiện mô hình Geoid động lực học có độ chính xác $\le \pm 2\text{ cm}$ và phát triển cổng dịch vụ định vị điểm chính xác trực tuyến quốc gia (VN-PPP).

Kết luận

Luận án tiến sĩ của tác giả Lương Thanh Thạch đã giải quyết trọn vẹn và xuất sắc một vấn đề khoa học - công nghệ có tính cấp thiết chiến lược của ngành Kỹ thuật Trắc địa - Bản đồ Việt Nam. Tựu trung lại, công trình ghi dấu 5 đóng góp cốt lõi:

  1. Xác lập hoàn chỉnh cơ sở khoa học, hệ thống tiêu chí kỹ thuật và phương pháp luận xây dựng Hệ quy chiếu tọa độ không gian quốc gia (HQCTĐKGQG).
  2. Kiểm chứng thực nghiệm toàn diện và khẳng định mô hình quasigeoid quốc gia khởi đầu VIGAC2017 đạt độ chính xác $\pm 6.4\text{ cm}$, đủ điều kiện làm nền tảng chuyển đổi độ cao trắc địa phục vụ xây dựng HQCTĐKGQG.
  3. Hoàn thiện phương pháp bình sai ghép nối mạng lưới GNSS quốc gia phủ trùm vào HQCTĐKGQG, nâng cao độ chính xác vị trí không gian lên 2 lần và chính xác hóa bộ 7 tham số chuyển đổi liên kết với ITRF.
  4. Thiết lập thành công thuật toán và module chương trình bình sai truy hồi với phép biến đổi xoay, cho phép mở rộng, cập nhật điểm cơ sở không gian mới từ CSDL trắc địa quốc gia mà không cần bình sai lại toàn bộ mạng lưới gốc.
  5. Định hình mô hình quản lý, khai thác thống nhất hạ tầng dữ liệu trắc địa quốc gia, mở đường cho việc tích hợp công nghệ định vị vệ tinh dẫn đường hiện đại vào công cuộc phát triển kinh tế - xã hội, bảo vệ chủ quyền biển đảo và chủ động ứng phó với biến đổi khí hậu toàn cầu.