Tổng quan về luận án

Nghiên cứu chuyển động hiện đại và biến dạng vỏ Trái Đất đóng vai trò nền tảng trong trắc địa địa động lực và địa vật lý hiện đại, đặc biệt tại các khu vực có cấu trúc kiến tạo phức tạp như Đông Nam Á. Luận án tiến sĩ kỹ thuật của tác giả Lại Văn Thủy với đề tài "Nghiên cứu ứng dụng lý thuyết xấp xỉ sóng nhỏ (biến đổi Wavelet) để phân tích, nội suy vận tốc chuyển dịch và biến dạng không gian từ kết quả xử lý dữ liệu đo GPS mạng lưới trắc địa địa động lực khu vực miền Bắc Việt Nam" (chuyên ngành Kỹ thuật Trắc địa - Bản đồ, Mã số: 9.03, Viện Khoa học Đo đạc và Bản đồ, 2020) là công trình tiên phong giải quyết bài toán mô hình hóa không gian 3 chiều trường vận tốc và các đại lượng biến dạng vỏ Trái Đất trên mặt cầu.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: trong hơn hai thập kỷ qua tại Việt Nam, các nghiên cứu trắc địa địa động lực (Dương Chí Công, 2000; Vy Quốc Hải, 2005, 2009; Hà Minh Hòa, 2005, 2008; Phan Trọng Trịnh, 2010; Nguyễn Văn Hướng, 2012) chủ yếu tiếp cận trên các đa giác phẳng 2D, mô hình hóa phân tán cục bộ cho từng đới đứt gãy riêng lẻ hoặc sử dụng phép nội suy Kriging/Spline trên mặt phẳng chiếu $Oxy$, hoàn toàn bỏ qua thành phần chuyển dịch đứng $V_U$ và độ cong thực tế của mặt cầu Trái Đất.

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:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để phân tích, nội suy liên tục trường vận tốc chuyển dịch 3 chiều ($V_N, V_E, V_U$) trên mặt cầu từ mạng lưới quan trắc GNSS có mật độ phân bố không đồng đều và tồn tại nhiễu đo đạc?
  • Câu hỏi nghiên cứu 2 (RQ2): Cường độ, xu thế và đặc trưng phân bố không gian của các trường đại lượng biến dạng (tốc độ xoay $\omega$, tốc độ biến dạng $\dot{\varepsilon}$, tốc độ trượt $\dot{\gamma}$, tốc độ trương nở $\Theta$) tại khu vực miền Bắc Việt Nam giai đoạn 2012–2017 diễn ra như thế nào?
  • Giả thuyết khoa học 1 (H1): Phương pháp xấp xỉ sóng nhỏ DOG (Difference of Gaussians) trên mặt cầu cho phép lọc trừ hiệu quả sai số thô, tích hợp phân tích đa phân giải (MRA) và tái tạo trường vận tốc chuyển dịch địa phương với độ ổn định số học vượt trội so với phương pháp tam giác truyền thống của Frank (1966).
  • Giả thuyết khoa học 2 (H2): Hoạt động kiến tạo hiện đại tại miền Bắc Việt Nam mang bản chất biến dạng nội mảng cường độ nhỏ, tập trung cục bộ xung quanh đới đứt gãy thạch quyển cấp II Lai Châu – Điện Biên thay vì biến dạng phân tán đồng đều.

Khung lý thuyết của nghiên cứu tích hợp Thuyết Kiến tạo mảng (Plate Tectonics), Cơ học môi trường liên tục (Continuum Mechanics) và Biến đổi sóng nhỏ đa phân giải (Multiresolution Wavelet Analysis). Phạm vi thực nghiệm bao phủ diện tích bề mặt cầu $S = 1.7358 \times 10^{11} \text{ m}^2$ (tương ứng bán kính tương đương $R = 4.7012 \times 10^5 \text{ m}$), giới hạn từ $19.5^\circ$ đến $23.0^\circ$ vĩ Bắc và $103.0^\circ$ đến $108.0^\circ$ kinh Đông, dựa trên chuỗi số liệu 6 chu kỳ đo GNSS chuẩn hóa (2012–2017) theo Quyết định số 1665/QĐ-BTNMT.

Literature Review và Positioning

Nghiên cứu biến dạng vỏ Trái Đất bằng kỹ thuật trắc địa không gian đã trải qua nhiều bước chuyển dịch mô hình (paradigm shifts). Phương pháp tam giác cổ điển do Frank (1966) đề xuất giả định trường biến dạng đồng nhất trong từng tam giác đơn lẻ, bộc lộ nhược điểm lớn là tính gián đoạn tại các cạnh biên và cực kỳ nhạy cảm với sai số phân bố trạm. Để khắc phục, Haines et al. (1998) đưa ra phương pháp nội suy hàm Spline bậc ba (bicubic splines), trong khi Spakman & Nyst (2002) phát triển kỹ thuật nghịch đảo biến thiên cạnh cơ sở (baseline-rate inversion). Tuy nhiên, các kỹ thuật này gặp hạn chế nghiêm trọng khi xử lý mạng lưới có độ phân bố trạm không đồng nhất hoặc khi có sự nhảy bậc vận tốc tại các đới đứt gãy hoạt động.

Bước đột phá phương pháp luận quốc tế xuất hiện khi Tape et al. (2009) giới thiệu lý thuyết xấp xỉ sóng nhỏ trên mặt cầu để ước tính trường vận tốc và ten-xơ biến dạng liên tục tại đới đứt gãy San Andreas (khu vực Parkfield, California, Mỹ), phát hiện tốc độ biến dạng chính đạt khoảng $1.5,\mu\text{rad/năm}$. Kỹ thuật này chứng minh khả năng tự động tách lọc tín hiệu đa tỷ lệ và phát hiện các biến động địa chấn chuyển tiếp. Tiếp đó, Cheng Pengfei et al. (2015) đã ứng dụng thành công sóng nhỏ hình cầu trên mạng lưới GPS lục địa Trung Quốc, đạt độ lệch chuẩn $RMS = 0.95,\text{mm/năm}$ (hướng Đông) và $0.97,\text{mm/năm}$ (hướng Bắc) ở vùng mạng lưới dày đặc; $1.32,\text{mm/năm}$ (hướng Đông) và $1.30,\text{mm/năm}$ (hướng Bắc) ở khu vực trạm thưa thớt. Tương tự, Kumar et al. (2003) ứng dụng sóng nhỏ phân tích chuỗi GPS tại miền Trung Nhật Bản (1996–1999), xác định tỷ lệ biến dạng trượt dao động từ $0.12$ đến $0.18,\mu\text{strain/năm}$.

- Frank (1966)                 - Haines et al. (1998)         - Tape et al. (2009) [Mỹ]
- Dương Chí Công (2000)        - Spakman & Nyst (2002)        - Cheng Pengfei et al. (2015) [TQ]
- Hà Minh Hòa (2005, 2008)     - Nguyễn Văn Hướng (2012)     - Lại Văn Thủy (2020) [VN]

Tại Việt Nam, các công trình trước đây của Dương Chí Công (2000) trên vùng Thác Bà – Yên Bái (giai đoạn 1963–1994 xác định biến dạng trượt đứt gãy Sông Hồng $< 0.6,\mu\text{rad/năm}$, tương ứng dịch ngang $6,\text{mm/năm}/10,\text{km}$), Vy Quốc Hải (2009) sử dụng GPSurvey 2.2 tại Tam Đảo – Ba Vì (điểm HUN1 dịch chuyển $33.4,\text{mm/năm}$ Đông, $12.6,\text{mm/năm}$ Nam trong ITRF-2000), và đề tài cấp Nhà nước KC.09.11/06-10 do Phan Trọng Trịnh (2010) chủ nhiệm xác định vận tốc tuyệt đối mạng lưới Biển Đông đạt $22–25,\text{mm/năm}$ hướng Đông và $5–8,\text{mm/năm}$ hướng Nam trong ITRF-08. Dù vậy, tất cả đều dừng lại ở việc tính toán vận tốc rời rạc trên mặt phẳng hoặc nội suy Kriging phẳng (Nguyễn Văn Hướng, 2012). Luận án của Lại Văn Thủy đã định vị chính xác khoảng trống khoa học này bằng việc chuyển dịch hoàn toàn bài toán sang hệ tọa độ không gian mặt cầu 3 chiều thông qua sóng nhỏ DOG.

Đó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 lý thuyết cơ học môi trường liên tục trên đa tạp mặt cầu $T^2(S)^2$ và lý thuyết biến đổi sóng nhỏ rời rạc (DWT), thiết lập mối liên kết toán học chặt chẽ giữa toán tử vi phân không gian và động học thạch quyển. Cốt lõi của mô hình là hàm sóng nhỏ hình cầu DOG (Difference of Gaussians) xác định theo biểu thức:

$$\Psi_{x,a}(\theta', \phi') = \frac{1}{\sqrt{\kappa(\alpha, a)}} \left( \frac{1}{a} e^{-\frac{1-\cos\gamma}{a^2}} - \frac{1}{\alpha a} e^{-\frac{1-\cos\gamma}{\alpha^2 a^2}} \right)$$

Trong đó $x(\theta, \phi)$ là cực trung tâm sóng nhỏ trên mặt cầu, $\gamma$ là góc ở tâm giữa điểm lưới $(\theta, \phi)$ và điểm tính toán $(\theta', \phi')$, tham số co giãn $a_q = 2^{-q}$ ứng với các bậc phân giải rời rạc hóa $q$, và $\alpha > 1$ là hằng số tỷ lệ tối ưu hóa biên độ.

Về mặt động học mảng, công trình phát triển mô hình tách vector chuyển động kiến tạo tổng thể thành hai thành phần độc lập: chuyển động xoay cứng khối mảng ổn định mô tả bởi vector quay Ơ-le ($\Omega_x, \Omega_y, \Omega_z$) và trường vận tốc chuyển dịch nội khối địa phương sinh ra biến dạng đàn hồi thuần túy.

          Trường Vận tốc Tuyệt đối V_GNSS (ITRF-2008)
   Chuyển động Mảng Ơ-le            Trường Vận tốc Địa phương
(Xoay khối Nam Trung Hoa / Sunda)     (Nội suy Wavelet DOG mặt cầu)
         Tensor Vận tốc Biến dạng D                       Tensor Tốc độ Xoay W
            (Thành phần đối xứng)                        (Thành phần phản đối xứng)
Tốc độ Trượt (γ_dot)     Tốc độ Trương nở (Θ)                 Tốc độ Xoay (ω)

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

Khung phân tích của luận án tích hợp đồng thời ba cấu phần lý thuyết: Động học mảng kiến tạo, Giải tích đa phân giải wavelet, và Tensor gradient vận tốc bậc hai. Tensor gradient vận tốc $L = \nabla v$ được phân rã tường minh thành thành phần đối xứng $D$ (tensor tốc độ biến dạng) và thành phần phản đối xứng $W$ (tensor tốc độ xoay):

$$L = D + W = \frac{1}{2}(L + L^T) + \frac{1}{2}(L - L^T)$$

Từ đó, các đại lượng vật lý biến dạng vô hướng trên mặt cầu được thiết lập chính xác:

  • Tốc độ xoay (Rotation rate): $\omega = \frac{1}{2} \text{rot}(v) = \frac{1}{2} \left( \frac{\partial v_y}{\partial x} - \frac{\partial v_x}{\partial y} \right)$
  • Tensor tốc độ biến dạng (Strain rate tensor): $D_{ij} = \frac{1}{2} \left( \frac{\partial v_i}{\partial x_j} + \frac{\partial v_j}{\partial x_i} \right)$
  • Tốc độ trương nở diện tích/thể tích (Dilatation rate): $\Theta = \text{div}(v) = \text{tr}(D) = \frac{\partial v_x}{\partial x} + \frac{\partial v_y}{\partial y} + \frac{\partial v_z}{\partial z}$

Điều kiện biên được kiểm soát chặt chẽ thông qua thuật toán chuẩn hóa tối ưu OCV (Ordinary Cross-Validation) nhằm cân bằng giữa độ trơn bề mặt và độ khớp chi tiết cục bộ, loại bỏ nguy cơ quá khớp (over-fitting) tại các trạm biên.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường thực chứng (positivism) với thiết kế định lượng nghiêm ngặt kết hợp mô hình hóa toán tin trắc địa. Thiết kế không gian đa tầng phân cấp bao gồm:

  • Tầng 1: Khung quy chiếu trắc địa quốc tế ITRF-2008 (gốc IGS).
  • Tầng 2: Khối kiến tạo khu vực Đông Nam Á (khối Nam Trung Hoa và khối Sundaland).
  • Tầng 3: Mạng lưới địa động lực 3D miền Bắc Việt Nam phủ trùm 11 đới đứt gãy thạch quyển cấp II ($F_{II}$) và 34 đới đứt gãy xuyên vỏ cấp III ($F_{III}$).
  • Tầng 4: Lưới điểm nút tam giác cầu (Spherical TIN GRID) đa tỷ lệ phân giải.

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu tuân thủ các chuẩn mực đo đạc thiên văn - vệ tinh quốc tế:

  1. Thu nhận và chuẩn hóa dữ liệu: Sử dụng máy thu GNSS 2 tần số geodetic, đo lặp 6 chu kỳ (2012–2017). Chiều cao anten được kiểm chuẩn nghiêm ngặt bằng cả phương pháp đo nghiêng và đo thẳng đứng nhằm triệt tiêu sai số tâm pha anten.
  2. Xử lý số liệu cơ sở: Tính toán bình sai lưới bằng phần mềm chuyên dụng BERNESE 5.2 (Đại học Bern, Thụy Sỹ), ràng buộc chặt chẽ với các trạm IGS quốc tế xung quanh khu vực (như HOCM, DOHO, HUES, STT1, BLV1, LANG).
  3. Phân tích đa phân giải và nội suy sóng nhỏ: Chia mặt cầu ban đầu $G_0$ gồm 20 mặt tam giác thành các cấp lưới $G_q$ với bậc rời rạc hóa $q = 7, 8, 9$. Chỉ các tam giác cầu chứa $\ge 3$ điểm trắc địa mới được chọn để tham gia tính toán hệ số sóng nhỏ nhằm đảm bảo độ tin cậy thống kê.

Data và phân tích

Đặc trưng dữ liệu không gian được kiểm soát ở mức milimet:

  • Dữ liệu GNSS vector vận tốc 3 thành phần ($V_N, V_E, V_U$) được hiệu chỉnh triệt để sai số đa đường truyền (multipath) và khúc xạ tầng đối lưu/tầng điện ly.
  • Tham số chính tắc hóa (regularization parameters) được xác định độc lập cho từng hướng $V_N, V_E, V_U$ bằng thuật toán OCV.
  • Toàn bộ thuật toán nội suy sóng nhỏ DOG trên mặt cầu, phân rã tensor gradient và tính toán trường biến dạng được nghiên cứu sinh lập trình, module hóa hoàn chỉnh trên môi trường MATLAB với giao diện GUI trực quan.

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

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

Các phát hiện chính từ kết quả xử lý dữ liệu và tính toán thực nghiệm:

[!IMPORTANT] Trích dẫn minh chứng từ luận án:

  1. "xác định trường vận tốc chuyển dịch và vận tốc biến dạng biểu diễn trên không gian mặt cầu, (2) phát hiện được trị đo thô, (3) xác định vận tốc chuyển dịch và vận tốc biến dạng trên khu vực có mật độ điểm đo trắc địa không đồng đều" (Luận án, Mục Mở đầu, tr. 3).
  2. "vận tốc chuyển dịch nhận được trong 6 chu kỳ đo (từ năm 2012-2017) đã cho thấy hoạt động của các đứt gãy trên khu vực miền Bắc Việt Nam có tốc độ nhỏ và tập trung chủ yếu ở vùng Tây Bắc (xung quanh đứt gãy Lai Châu - Điện Biên)" (Luận án, Luận điểm bảo vệ 2, tr. 4).
  3. "Kết quả tính toán tốc độ chuyển dịch bằng phần mềm BERNESE 5.0 đã xác định được vận tốc chuyển động kiến tạo hiện đại trong hệ quy chiếu trắc địa quốc tế ITRF-08 tại 12 điểm đo với tốc độ chuyển dịch về phía Đông khoảng 22-25 mm/năm và chuyển dịch về phía Nam khoảng 5-8 mm/năm" (Luận án, Chương 1, tr. 22).
Đới đứt gãy / Khu vực Kiểu động học Vận tốc dịch chuyển địa phương Đặc trưng trường biến dạng chính
Lai Châu – Điện Biên ($F_{II}.1$) Trượt bằng phải kết hợp tách giãn $3.0 – 6.0\text{ mm/năm}$ Tốc độ trượt cao, tập trung gradient biến dạng lớn nhất miền Bắc
Sông Hồng ($F_{II}.2$) Trượt bằng phải yếu $1.0 – 2.5\text{ mm/năm}$ Biến dạng cắt nhỏ, vận tốc tương đối hai cánh không đáng kể
Sông Lô ($F_{II}.3$) & Đông Bắc Trượt bằng / Ổn định $< 1.5\text{ mm/năm}$ Trường tốc độ xoay và trương nở tiệm cận $0$
Sông Mã – Sơn La ($F_{III}.4, F_{III}.5$) Trượt bằng $2.0 – 4.0\text{ mm/năm}$ Nén ép diện tích nhẹ, dị thường biến dạng cục bộ
  1. Đặc trưng trường vận tốc tuyệt đối: Toàn bộ khu vực miền Bắc Việt Nam dịch chuyển theo hướng Đông Nam với tốc độ tuyệt đối từ $22$ đến $32\text{ mm/năm}$ trong khung quy chiếu ITRF-2008, phản ánh sự trôi dạt tổng thể của khối Nam Trung Hoa dưới tác động chèn ép của khối Ấn - Úc vào mảng Á - Âu.
  2. Cường độ chuyển động nội khối nhỏ: Sau khi tách vector quay Ơ-le mảng, vận tốc chuyển dịch địa phương nội khối chỉ dao động trong khoảng $0.5 - 5.0\text{ mm/năm}$, chứng minh thạch quyển khu vực miền Bắc có độ cứng vững tương đối cao.
  3. Tập trung biến dạng cục bộ tại Tây Bắc: Hoạt động kiến tạo và dị thường biến dạng hiện đại không phân tán mà tập trung rõ rệt tại vùng Tây Bắc, đặc biệt dọc đới đứt gãy thạch quyển Lai Châu – Điện Biên ($F_{II}.1$), thể hiện cơ chế trượt bằng phải kết hợp tách dãn.
  4. Vượt trội phương pháp luận: So sánh đối chứng với phương pháp tam giác Frank (1966), phương pháp xấp xỉ sóng nhỏ DOG trên mặt cầu đã loại trừ hoàn toàn hiện tượng phân mảnh giả tạo, tái hiện trường biến dạng liên tục, mịn màng và bắt trọn các dị thường gradient tại các khu vực trạm đo phân bố thưa thớt.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở định lượng chuẩn xác bổ sung cho mô hình động học thạch quyển Đông Nam Á, khẳng định vai trò ranh giới động lực học của đới đứt gãy Lai Châu - Điện Biên trong giai đoạn hiện đại.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình xử lý số liệu GNSS địa động lực đa quy mô (BERNESE $\rightarrow$ OCV $\rightarrow$ DOG Spherical Wavelets $\rightarrow$ Deformation Tensors), có thể chuyển giao áp dụng trực tiếp cho các mạng lưới quan trắc trên toàn lãnh thổ Việt Nam và khu vực Biển Đông.
  • Ứng dụng thực tiễn và chính sách: Cung cấp bản đồ phân vùng biến dạng chính xác phục vụ đánh giá rủi ro động đất kích thích cho các công trình trọng điểm quốc gia (các bậc thang thủy điện trên lưu vực Sông Đà: Sơn La, Hòa Bình, Lai Châu), đồng thời là luận cứ khoa học để Bộ Tài nguyên và Môi trường tối ưu hóa quy hoạch mạng lưới trạm định vị vệ tinh quốc gia VNGEONET.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 3 hạn chế chính:

  1. Độ phân bố trạm chưa đồng đều: Mật độ mốc quan trắc tập trung dày ở các tuyến đứt gãy trọng điểm nhưng còn thưa ở vùng biên giới hiểm trở, khiến bậc phân giải sóng nhỏ tại một số vùng chỉ dừng ở $q=7$.
  2. Chuỗi số liệu đo định kỳ (campaign-style): Khoảng thời gian 6 chu kỳ (2012–2017) dù đủ xác định xu thế tuyến tính giữa các kỳ động đất (interseismic) nhưng chưa đủ dài để tách biệt các hiệu ứng phi tuyến tính theo mùa (seasonal variations) hoặc biến dạng chuyển tiếp chậm.
  3. Độ chính xác thành phần độ cao $V_U$: Do ảnh hưởng của khúc xạ tầng đối lưu dư thừa, độ chính xác xác định vận tốc đứng $V_U$ vẫn thấp hơn thành phần ngang $V_N, V_E$.

Hướng nghiên cứu tiếp theo:

  • Tích hợp mạng lưới trạm GNSS liên tục (CORS/VNGEONET) với dữ liệu quan sát giao thoa Radar đa thời gian (InSAR - PSInSAR/SBAS) nhằm tạo ra trường quan trắc mật độ siêu cao.
  • Mở rộng thuật toán sóng nhỏ DOG để giải bài toán nghịch đảo 3 chiều nguồn trượt đứt gãy sâu (slip deficit inversion) trong môi trường đàn - nhớt (viscoelastic modeling).

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

Công trình tạo ra bước chuyển biến rõ rệt cho nghiên cứu trắc địa địa động lực:

  • Học thuật: Xác lập phương pháp luận mới trong phân tích dữ liệu GNSS địa động lực tại Việt Nam, mở rộng tiềm năng trích dẫn trong các nghiên cứu kiến tạo khu vực Đông Nam Á.
  • Công nghiệp & Kỹ thuật: Bộ phần mềm MATLAB tự xây dựng cung cấp công cụ tính toán tự động hóa hoàn toàn cho các viện nghiên cứu chuyên ngành, doanh nghiệp khảo sát công trình ngầm và thủy điện lớn.
  • Quản lý nhà nước: Cung cấp số liệu đầu vào chuẩn hóa phục vụ Chiến lược quốc gia về phòng chống và giảm nhẹ thiên tai, hỗ trợ Viện Vật lý Địa cầu trong việc cập nhật bản đồ phân vùng vi địa chấn quốc gia.

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

  • Nghiên cứu sinh & Nhà khoa học Trắc địa - Địa chất: Tiếp cận hệ thống lý thuyết sóng nhỏ hình cầu DOG và thuật toán xử lý tensor biến dạng không gian 3 chiều tiên tiến.
  • Kỹ sư R&D Công trình thủy điện/Giao thông: Sử dụng bản đồ trường biến dạng để tính toán hệ số an toàn kháng chấn cho công trình ngầm, cầu lớn và đập thủy điện cao trình lớn.
  • Cơ quan Hoạch định Chính sách (Bộ Tài nguyên & Môi trường, Bộ Xây dựng): Có cơ sở dữ liệu tin cậy để quy hoạch phát triển không gian đô thị và hạ tầng kỹ thuật tại vùng Tây Bắc.

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ì? Đó là việc mở rộng thành công lý thuyết xấp xỉ sóng nhỏ trên mặt cầu DOG kết hợp phân tích tensor gradient vận tốc 3D để mô hình hóa trường biến dạng liên tục tại Việt Nam, thay thế hoàn toàn tư duy xấp xỉ 2D phẳng truyền thống.

  2. Phương pháp xấp xỉ sóng nhỏ vượt trội hơn các nghiên cứu quốc tế trước đây ở điểm nào? So với phương pháp nội suy Spline của Haines et al. (1998) vốn dễ làm mờ các dị thường cục bộ, và phương pháp tam giác của Frank (1966) bị đứt gãy đường biên, phương pháp sóng nhỏ DOG của luận án (kế thừa Tape et al., 2009 và Cheng Pengfei et al., 2015) cho phép phân tích đa độ phân giải thích ứng ($q=7,8,9$) theo mật độ trạm, tự động lọc nhiễu và bảo toàn biên độ gradient biến dạng tại đới đứt gãy.

  3. Phát hiện thực nghiệm gây bất ngờ nhất là gì? Đới đứt gãy Sông Hồng ($F_{II}.2$) – từng được xem là ranh giới kiến tạo chuyển động mạnh trong lịch sử – hiện nay có tốc độ chuyển dịch tương đối rất nhỏ ($1.0 - 2.5\text{ mm/năm}$), trong khi toàn bộ ứng suất và tốc độ biến dạng trượt hiện đại lại dồn tích chủ yếu về đới đứt gãy Lai Châu – Điện Biên ($F_{II}.1$).

  4. Quy trình nghiên cứu có khả năng tái lập (replication) không? Hoàn toàn có thể tái lập nhờ mô tả toán học chi tiết, các tham số sóng nhỏ tường minh ($a_q, \kappa, \alpha$), quy chuẩn xử lý BERNESE 5.2 cụ thể và bộ mã nguồn phần mềm được xây dựng hoàn chỉnh trên nền tảng MATLAB.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Định hướng phát triển mô hình đồng hóa dữ liệu đa nguồn: GNSS liên tục + InSAR + Trọng lực vệ tinh (GRACE-FO/GOCE) để quan trắc biến dạng 4D (không gian - thời gian) theo thời gian thực phục vụ cảnh báo sớm tai biến địa chất.

Kết luận

  1. Hoàn thiện khung lý thuyết: Luận án đã xây dựng thành công cơ sở khoa học và bộ công cụ toán học hoàn chỉnh ứng dụng lý thuyết xấp xỉ sóng nhỏ DOG trên mặt cầu để giải đoán trường biến dạng 3D thạch quyển.
  2. Giải quyết triệt để bài toán trắc địa thực tiễn: Xử lý và chuẩn hóa thành công chuỗi dữ liệu 6 chu kỳ đo GNSS (2012–2017) trên mạng lưới địa động lực miền Bắc, thiết lập hệ quy chiếu thống nhất trong ITRF-2008.
  3. Phát hiện quy luật động học then chốt: Làm sáng tỏ xu thế dịch chuyển Đông Nam ($22–32\text{ mm/năm}$) và khẳng định đới đứt gãy Lai Châu – Điện Biên là trung tâm dồn tích biến dạng hiện đại lớn nhất khu vực miền Bắc.
  4. Phát triển phần mềm chuyên dụng: Đóng gói trọn vẹn bộ module tính toán trên nền tảng MATLAB, cung cấp giải pháp kỹ thuật tự chủ, độ tin cậy cao cho ngành trắc địa bản đồ Việt Nam.
  5. Mở ra hướng nghiên cứu mới: Đặt nền móng vững chắc cho các nghiên cứu mô hình hóa geodynamics 4D đa tỷ lệ và phân tích nguy cơ tai biến địa chất tiệm cận trình độ khoa học quốc tế.