Tổng quan về luận án

Nghiên cứu động lực học vỏ Trái Đất và biến dạng kiến tạo hiện đại đóng vai trò cốt lõi trong việc đánh giá tai biến địa chất, dự báo động đất và hiện đại hóa hạ tầng trắc địa quốc gia. Luận án tiến sĩ kỹ thuật trắc địa bản đồ của Nghiên cứu sinh Lại Văn Thủy (2020) với tiêu đề "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", thực hiện tại Viện Khoa học Đo đạc và Bản đồ dưới sự hướng dẫn khoa học của GS. Đặng Hùng Võ và TS. Dương Chí Công, đại diện cho bước chuyển biến mang tính bước ngoặt từ các phương pháp xử lý biến dạng hình học phẳng hai chiều (2D) truyền thống sang mô hình phân tích đa phân giải ba chiều (3D) trên mặt cầu thực thể.

                                  KHUNG NGHIÊN CỨU TỔNG THỂ
[DỮ LIỆU GNSS ĐỊA ĐỘNG LỰC]                                             [LÝ THUYẾT SÓNG NHỎ CẦU DOG]
- 6 chu kỳ đo (2012 - 2017)                                              - Không gian hàm T²(S²)
- Khung quy chiếu ITRF-2008                                              - Đa phân giải: Bậc q = 7, 8, 9
- 22 đới đứt gãy Miền Bắc                                               - Tối ưu hóa chính tắc OCV
                             [KHUNG PHÂN TÍCH KHÔNG GIAN 3D]
                       Véc-tơ chuyển dịch = Vận tốc quay Ơ-le + Biến dạng địa phương
                           [TENSOR GRADIENT VẬN TỐC (L = D + W)]
[TENSOR VẬN TỐC BIẾN DẠNG (D)]                                         [TENSOR TỐC ĐỘ XOAY (W)]
- Tốc độ biến dạng trượt (Shear rate)                                   - Tốc độ xoay (Rotation rate)
- Tốc độ trương nở (Dilatation rate Θ)                                  - Phân tích vặn xoắn kiến tạo

Bối cảnh khoa học và khoảng trống nghiên cứu (Research Gap)

Trước khi công trình này được công bố, các phân tích biến dạng vỏ Trái Đất tại Việt Nam chủ yếu dựa trên:

  1. Phương pháp mạng lưới tam giác phẳng của Frank (1966).
  2. Các thuật toán nội suy đa thức, hàm Spline bậc ba (Haines et al., 1998).
  3. Mô hình khối cứng rời rạc hai chiều chỉ xét đến hai thành phần vận tốc mặt bằng ($V_N, V_E$) và bỏ qua thành phần chuyển dịch thẳng đứng ($V_U$).

Các phương pháp cổ điển này bộc lộ hạn chế nghiêm trọng:

  • Giả định tính đồng nhất đàn hồi trên toàn bộ mạng lưới tam giác.
  • Gây ra sai số biên lớn tại những khu vực có mật độ trạm đo thưa thớt hoặc phân bố bất đối xứng.
  • Bất khả thi trong việc tách lọc các dị thường cục bộ bậc cao khỏi xu thế chuyển động mảng kiến tạo khu vực.

Khoảng trống học thuật chính là sự thiếu vắng một công cụ toán học đa tỷ lệ (multiscale) có khả năng định xứ không gian - tần số đồng thời trên mặt cầu $S^2$, nhằm xử lý mạng lưới trắc địa phân bố không đều và đánh giá đồng bộ tensor biến dạng không gian 3 chiều.

Câu hỏi nghiên cứu và Giả thuyết khoa học

Luận án tập trung giải quyết các câu hỏi và giả thuyết nghiên cứu cốt lõi:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một khung toán học xấp xỉ sóng nhỏ mặt cầu (Spherical Wavelets) nhằm nội suy chính xác trường vận tốc chuyển dịch 3D từ các chuỗi số liệu GPS phân bố không đồng đều?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ tích tụ ứng suất, tốc độ trượt cắt và tốc độ trương nở thực tế tại các đới đứt gãy chính ở miền Bắc Việt Nam (đặc biệt là đứt gãy Sông Hồng, Sông Chảy, Sơn La, Lai Châu - Điện Biên) diễn ra như thế nào trong giai đoạn 2012–2017?
  • Giả thuyết 1 (H1): Hàm sóng nhỏ dạng hiệu của hai hàm Gauss (Difference of Gaussians - DOG) trên mặt cầu $S^2$ kết hợp chuẩn hóa Tikhonov cho phép trích xuất tín hiệu biến dạng kiến tạo cục bộ với độ tin cậy vượt trội so với phương pháp Frank truyền thống.
  • Giả thuyết 2 (H2): Hoạt động kiến tạo hiện đại tại miền Bắc Việt Nam mang tính phân dị sâu sắc, trong đó hoạt động trượt và tích tụ biến dạng tập trung chủ yếu ở phía Tây Bắc (đới Lai Châu - Điện Biên), trong khi đới đứt gãy Sông Hồng duy trì trạng thái tích tụ biến dạng ở mức thấp đến trung bình.

Phạm vi và Quy mô nghiên cứu

  • Phạm vi không gian: Tọa độ từ $19.5^\circ\text{N}$ đến $23.0^\circ\text{N}$ vĩ độ Bắc và từ $103.0^\circ\text{E}$ đến $108.0^\circ\text{E}$ kinh độ Đông. Diện tích mặt cầu phân tích $S = 1.7358 \times 10^{11} \text{ m}^2$, tương ứng bán kính vùng bao phủ $R = 4.7012 \times 10^5 \text{ m}$.
  • Phạm vi thời gian và Dữ liệu: Chuỗi dữ liệu quan trắc 6 chu kỳ liên tục từ năm 2012 đến 2017 thuộc Đề án xây dựng mạng lưới trắc địa địa động lực theo Quyết định số 1665/QĐ-BTNMT của Bộ Tài nguyên và Môi trường, kết hợp mạng lưới trạm IGS toàn cầu và khu vực.

Literature Review và Định vị học thuật

Tiến trình nghiên cứu địa động lực bằng trắc địa tại Việt Nam

Nghiên cứu chuyển dịch vỏ Trái Đất bằng kỹ thuật trắc địa tại Việt Nam khởi đầu từ năm 1980 với mạng lưới tam giác khe nứt Chí Linh của Cục Đo đạc và Bản đồ Nhà nước (Đặng Hùng Võ, 1980), tiếp nối bằng các tuyến đo thủy chuẩn hình học quan trắc chuyển dịch đứng (Nguyễn Khắc Vinh, 1986). Từ năm 1993, việc ứng dụng công nghệ GPS đã mở rộng khảo sát tại đới đứt gãy Sông Hồng và Sông Chảy (Dương Chí Công, 1994, 2005). Các nghiên cứu hợp tác quốc tế Việt - Đức, Việt - Pháp (1994–2000) và các công trình của Vi Quốc Hải (2007), Nguyễn Tuấn Anh (2010), Hà Minh Hòa (2004, 2008), Phan Trọng Trịnh (2010) thuộc Đề tài trọng điểm cấp Nhà nước KC.08/06-10 đã từng bước xác định vận tốc tuyệt đối của mảng Sundaland và khối Nam Trung Hoa.

Tuy nhiên, như tác giả Lại Văn Thủy nhận định:

"Các kết quả tính toán, phân tích chuyển dịch và biến dạng mới chỉ thực hiện riêng lẻ cho từng đới đứt gãy hoặc trên khu vực nhỏ theo các đa giác phẳng; việc nghiên cứu đánh giá chuyển dịch của vỏ Trái Đất theo quy mô khu vực trong mô hình không gian hiện chưa được nghiên cứu ở Việt Nam."

Tổng quan quốc tế về Lý thuyết sóng nhỏ trong Trắc địa

Biến đổi sóng nhỏ (Wavelet Transform) xuất hiện nhằm khắc phục triệt để nhược điểm thiếu tính định xứ thời gian - không gian của biến đổi Fourier truyền thống (Chui, 1992; Daubechies, 1992). Trong trắc địa địa vật lý quốc tế:

  • Sideris (2001) ứng dụng xấp xỉ sóng nhỏ rời rạc (DWT) để lọc nhiễu trọng lực hàng không.
  • Isabelle Panet, Yuki Kuroishi, Matthias Holschneider (2006) kết hợp dữ liệu GRACE và trọng lực mặt đất để xây dựng mô hình trọng trường Nhật Bản bằng sóng nhỏ cầu.
  • Kaloop & Kim (2014), Gao et al. (2017) áp dụng Wavelet để lọc đa đường dẫn (multipath) trong chuỗi thời gian GPS-RTK.
                    TIẾN TRÌNH PHÁT TRIỂN PHƯƠNG PHÁP LUẬN
                    
  1966: Frank                1998: Haines et al.        2009: Tape et al.          2020: Lại Văn Thủy
  - Đàn hồi đồng nhất        - Cải tiến độ trơn         - Tích hợp đa phân giải     - 6 chu kỳ GPS Miền Bắc
  - Lỗi biên mạng lưới       - Vẫn phụ thuộc phẳng      - Khảo sát đứt gãy San      - Tách véc-tơ Euler 3D
                                                          Andreas (Parkfield)       - Minh giải 4 đại lượng

Đặc biệt, bước đột phá về mô hình hóa trường biến dạng mặt cầu được khởi xướng bởi Tape et al. (2009) khi ứng dụng sóng nhỏ DOG phân tích mạng lưới GPS tại Parkfield (đứt gãy San Andreas, Hoa Kỳ), phát hiện tốc độ biến dạng chính đạt xấp xỉ $1.5\ \mu\text{rad/năm}$. Kỹ thuật này sau đó được Cheng Pengfei et al. (2015) triển khai trên mạng lưới 1.000 trạm GPS lục địa Trung Quốc, đạt độ lệch bình phương trung bình giữa mô hình và thực đo là $0.95\text{ mm/năm}$ (hướng Đông) và $0.97\text{ mm/năm}$ (hướng Bắc) ở vùng dày đặc trạm, và $1.32\text{ mm/năm}$ ở vùng trạm thưa.

Tranh luận học thuật và Sự định vị của Luận án

Tiêu chí so sánh Phương pháp Tam giác Frank (1966) Phương pháp Spline cầu (Haines et al., 1998) Luận án Lại Văn Thủy (2020) [Dựa trên Tape et al., 2009]
Không gian toán học Hình học phẳng Euclidean 2D ($R^2$) Bề mặt cầu toán học ($S^2$) Không gian hàm khả tích mặt cầu $T^2(S^2)$
Bản chất phép nội suy Cục bộ từng tam giác đơn lẻ Làm trơn toàn cục (Global Smoothing) Đa phân giải đa tỷ lệ (Multiresolution Analysis)
Khả năng lọc nhiễu Rất nhạy cảm với sai số vị trí điểm Phụ thuộc trọng số làm trơn Khử nhiễu tự nhiên qua ngưỡng sóng nhỏ và OCV
Xử lý mật độ trạm thưa Gây biến dạng méo mó ở biên Dễ tạo ra các cực trị giả mạo Tự động thích ứng theo cấp phân giải $q$
Thành phần vận tốc Chỉ xử lý độc lập ($V_N, V_E$) Thường giới hạn ($V_N, V_E$) Đồng thời 3 thành phần không gian ($V_N, V_E, V_U$)

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

Đóng góp cho lý thuyết trắc địa động lực

Luận án đã mở rộng lý thuyết động học mảng kiến tạo hiện đại bằng cách tích hợp mô hình quay khối cứng Euler với lý thuyết biến dạng môi trường liên tục thông qua cơ sở hàm sóng nhỏ mặt cầu Difference of Gaussians (DOG).

                               KHUNG TOÁN HỌC SÓNG NHỎ MẶT CẦU
          [KHÔNG GIAN HÀM T²(S²)]                         [HÀM SÓNG NHỎ MẶT CẦU DOG]
          Phân rã đa cấp độ phân giải:                   Ψ_{a,x}(x') = (1/2πα²) [2 exp(-d²/2a²) 
          V_0 ⊂ V_1 ⊂ ... ⊂ V_{q_max}                                  - (1/2) exp(-d²/4a²)]
                             [LƯỚI ĐA PHÂN GIẢI TAM GIÁC CẦU G_q]
                       - Bậc q = 0: G_0 = 20 mặt (Icosahedron)
                       - Bậc q bất kỳ: N_q = 10 × 4^q + 2 đỉnh tam giác
                       - Tỷ lệ rời rạc hóa: a_q = 2^{-q}
                             [TENSOR GRADIENT VẬN TỐC TOÀN PHẦN]
                                           L = D + W
          [TENSOR VẬN TỐC BIẾN DẠNG: D]                     [TENSOR TỐC ĐỘ XOAY: W]
          D = (1/2)(L + L^T)                                W = (1/2)(L - L^T)
          - Tốc độ mở rộng / trượt                          - Tốc độ xoay vi phân
          - Tốc độ trương nở: Θ = tr(D)                     - Chiều hướng vặn xoắn

Cơ sở giải tích sóng nhỏ DOG trên mặt cầu

Không gian hàm bình phương khả tích trên mặt cầu $T^2(S^2)$ được phân rã thành các tầng tỷ lệ thông qua hàm cơ bản DOG: $$\Psi_{\alpha, x}(x') = \frac{1}{2\pi \alpha^2} \left[ 2\exp\left(-\frac{\Delta(\theta, \lambda; \theta', \lambda')^2}{2a^2}\right) - \frac{1}{2}\exp\left(-\frac{\Delta(\theta, \lambda; \theta', \lambda')^2}{4a^2}\right) \right]$$ Trong đó:

  • $\Delta(\theta, \lambda; \theta', \lambda')$ là khoảng cách góc cung lớn (góc ở tâm) giữa điểm cực trung tâm $x(\theta, \lambda)$ và điểm quan sát $x'(\theta', \lambda')$.
  • $a_q = 2^{-q}$ là tỷ lệ rời rạc hóa ở bậc $q$ ($q = 0, 1, 2, \dots, q_{\max}$).
  • $\alpha$ là hệ số tỷ lệ hình học ($\alpha > 1$). Khi $q$ không đổi, biên độ hàm sóng nhỏ tăng tỷ lệ thuận với $\alpha$.

Lưới phân chia tam giác cầu $G_q$ bắt đầu từ khối đa diện đều 20 mặt (Icosahedron, $q=0$, $N_0 = 12$ đỉnh). Ở bậc $q$, tổng số điểm nút lưới cầu là: $$N_q = 10 \cdot 4^q + 2$$ Đối với miền Bắc Việt Nam, tác giả xác định vùng không gian nghiên cứu có $S = 1.7358 \times 10^{11} \text{ m}^2$, lựa chọn bậc $q_{\min} = 7$ (khoảng cách nút lưới $\sim 40\text{ km}$), $q = 8$ ($\sim 20\text{ km}$) và $q_{\max} = 9$ ($\sim 10\text{ km}$) để nắm bắt đầy đủ các biến động kiến tạo từ cấp khu vực đến cấp đới đứt gãy hẹp.

Khung tensor vi phân xác định các đại lượng biến dạng

Từ trường vận tốc chuyển dịch địa phương $\vec{v}(x) = (v_\theta, v_\lambda, v_r)^T$ được xấp xỉ liên tục qua Wavelet, tensor gradient vận tốc $L = \nabla \vec{v}$ được phân tích thành hai thành phần trực giao: $$L = D + W$$

  1. Tensor tốc độ biến dạng (Đối xứng - Symmetric Strain Rate Tensor): $$D = \frac{1}{2}(L + L^T), \quad D_{ij} = \frac{1}{2}\left(\frac{\partial v_i}{\partial x_j} + \frac{\partial v_j}{\partial x_i}\right)$$ Bao gồm tốc độ biến dạng cắt cực đại (maximum shear strain rate $\dot{\gamma}_{\max}$) và tốc độ nén ép/giãn nở chính.

  2. Tensor tốc độ xoay (Phản đối xứng - Anti-symmetric Rotation Rate Tensor): $$W = \frac{1}{2}(L - L^T), \quad W_{ij} = \frac{1}{2}\left(\frac{\partial v_i}{\partial x_j} - \frac{\partial v_j}{\partial x_i}\right)$$ Đặc trưng cho vận tốc quay vi phân cứng của vi khối mà không làm biến đổi hình dạng thể tích.

  3. Tốc độ trương nở thể tích (Dilatation Rate $\Theta$): $$\Theta = \operatorname{tr}(D) = \operatorname{div}(\vec{v}) = \frac{\partial v_x}{\partial x} + \frac{\partial v_y}{\partial y} + \frac{\partial v_z}{\partial z}$$ Giá trị $\Theta > 0$ đặc trưng cho chế độ tách giãn địa chất (extension/dilation), $\Theta < 0$ biểu thị chế độ nén ép (compression/contraction), và $\Theta = 0$ xác lập trạng thái bảo toàn diện tích/thể tích không nén.

Phương pháp nghiên cứu thực nghiệm

                      QUY TRÌNH XỬ LÝ DỮ LIỆU ĐỊA ĐỘNG LỰC
                  [THU THẬP VÀ XỬ LÝ DỮ LIỆU GNSS (2012 - 2017)]
                  - 6 chu kỳ đo đồng bộ trên 22 đới đứt gãy
                  - Xử lý mạng lưới bằng BERNESE 5.2 / GAMIT
                  - Khung quy chiếu toàn cầu ITRF-2008 qua các trạm IGS
                  [TÍNH VẬN TỐC CHUYỂN DỊCH KHU VỰC VÀ ĐỊA PHƯƠNG]
                  - Xác định cực quay Ơ-le mảng nội khối: V_kv = Ω × R
                  - Tách chuyển dịch địa phương: V_dp = V_tuyetdoi - V_kv
                  [THIẾT LẬP LƯỚI SÓNG NHỎ CẦU VÀ NỘI SUY WAVELET]
                  - Xác định bán kính R = 4.7012 × 10^5 m
                  - Lập lưới tam giác cầu G_q (q = 7, 8, 9)
                  - Chuẩn hóa chính tắc Tikhonov & Đánh giá chéo OCV
                  [TÍNH TOÁN TENSOR BIẾN DẠNG VÀ THÀNH LẬP BẢN ĐỒ]
                  - Tensor gradient L -> D (Strain rate), W (Rotation rate)
                  - Tốc độ trượt (Shear), Tốc độ trương nở (Dilatation Θ)
                  - So sánh kiểm chứng chéo với phương pháp Frank (1966)

Thiết kế nghiên cứu và Chuẩn hóa dữ liệu GNSS

Dữ liệu quan trắc trắc địa địa động lực 6 chu kỳ (2012–2017) được xử lý bình sai nghiêm ngặt bằng phần mềm khoa học BERNESE 5.2 và GAMIT:

  • Sử dụng mô hình quỹ đạo vệ tinh chính xác cao của IGS (International GNSS Service).
  • Hiệu chỉnh pha sóng mang, tầng đối lưu (Troposphere ZTD) bằng mô hình VMF1/GMF kết hợp lập mô hình tải trọng đại dương (Ocean Tide Loading).
  • Tọa độ và véc-tơ vận tốc tuyệt đối được quy chiếu về hệ quy chiếu quốc tế ITRF-2008.

Thuật toán tách vận tốc quay Euler và chuyển dịch địa phương

Chuyển động tuyệt đối của một điểm trắc địa $\vec{V}{\text{tuyệt đối}}$ được phân tích thành: $$\vec{V}{\text{tuyệt đối}} = \vec{V}{\text{khu vực}} + \vec{V}{\text{địa phương}}$$ Trong đó, chuyển dịch khu vực $\vec{V}{\text{khu vực}}$ thể hiện sự quay cứng của mảng thạch quyển, xác định bởi véc-tơ quay Euler $\vec{\Omega} = (\Omega_X, \Omega_Y, \Omega_Z)^T$: $$\vec{V}{\text{khu vực}} = \vec{\Omega} \times \vec{R}$$ Nghiên cứu sinh đã lập trình module tính toán các tham số góc quay Euler và tọa độ cực quay cho khối kiến tạo miền Bắc Việt Nam, từ đó loại trừ chuyển động nền toàn cục để thu được trường vận tốc chuyển dịch địa phương thuần túy $\vec{V}_{\text{địa phương}}$ tại từng trạm đo.

                  CƠ CHẾ PHÂN RÃ CHUYỂN DỊCH KHÔNG GIAN

Xây dựng phần mềm tích hợp trên nền tảng Matlab

Tác giả đã tự thiết kế và phát triển gói phần mềm chuyên dụng trên nền tảng Matlab bao gồm 3 phân hệ chức năng:

  1. Module 1 - Phân tích động học Euler: Tính toán tham số quay mảng kiến tạo, tọa độ cực quay và trích xuất véc-tơ vận tốc địa phương.
  2. Module 2 - Nội suy Wavelet đa phân giải: Tạo lưới tam giác cầu $G_q$, tính ma trận thiết kế sóng nhỏ DOG, tối ưu hóa hệ số chính tắc hóa Tikhonov bằng thuật toán Ordinary Cross-Validation (OCV) nhằm triệt tiêu nhiễu đo ngẫu nhiên.
  3. Module 3 - Tính toán tensor biến dạng 3D: Khai triển đạo hàm riêng giải tích của các hệ số sóng nhỏ để xuất ra trường tốc độ biến dạng, tốc độ xoay, tốc độ trượt và tốc độ trương nở trên mạng lưới ô vuông (GRID).

Phát hiện đột phá và Ý nghĩa khoa học

Những phát hiện định lượng then chốt

                     PHÂN HÓA BIẾN DẠNG MIỀN BẮC VIỆT NAM
                     
  TÂY BẮC (Đới Lai Châu - Điện Biên)       ĐÔNG BẮC - ĐỒNG BẰNG BẮC BỘ (Đới Sông Hồng)
  - Tốc độ biến dạng: 1.5 - 3.2 × 10⁻⁸/yr  - Tốc độ biến dạng: 0.3 - 0.8 × 10⁻⁸/yr
  - Trượt bằng trái + Tách dãn             - Trượt phải chậm, ổn định kiến tạo
  - Trương nở diện tích: Dương (Dilation)  - Trương nở: Trung hòa / Co ngót nhẹ
  - Tích tụ ứng suất địa chấn cao          - Nguy cơ địa chấn tức thời thấp

Tác giả đã công bố luận điểm bảo vệ cốt lõi:

"Kết quả tính toán thực nghiệm với 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)."

  1. Phân hóa trường chuyển dịch ngang và đứng:

    • Trường vận tốc tuyệt đối ITRF-2008 cho thấy toàn bộ lãnh thổ miền Bắc chuyển dịch chủ đạo theo hướng Đông Nam với tốc độ từ $22\text{ mm/năm}$ đến $32\text{ mm/năm}$.
    • Chuyển dịch địa phương tương đối giữa hai bờ các hệ thống đứt gãy chính có độ lớn nhỏ, dao động trung bình từ $1.0\text{ mm/năm}$ đến $3.5\text{ mm/năm}$.
    • Chuyển dịch thẳng đứng ($V_U$) phân hóa thành các vệt trồi lún cục bộ xen kẽ với biên độ nâng hạ dao động từ $-4.0\text{ mm/năm}$ đến $+4.5\text{ mm/năm}$.

  2. Đặc trưng biến dạng đới đứt gãy Lai Châu - Điện Biên (FII.1):

    • Đạt tốc độ biến dạng trượt và biến dạng chính lớn nhất khu vực miền Bắc, dao động trong khoảng $(1.5 \div 3.2) \times 10^{-8}/\text{năm}$.
    • Tốc độ trương nở mang giá trị dương ($\Theta > 0$), phản ánh cơ chế trượt bằng kết hợp tách dãn tích cực. Đây là tâm điểm tích tụ năng lượng biến dạng kiến tạo hiện đại mạnh nhất toàn miền Bắc.

  3. Đặc trưng biến dạng đới đứt gãy Sông Hồng (FII.2) và Sông Lô (FII.3):

    • Duy trì trạng thái biến dạng thấp, tốc độ trượt cắt không vượt quá $0.8 \times 10^{-8}/\text{năm}$ (tương ứng vận tốc chuyển dịch trượt ngang dưới $1.5\text{ mm/năm}$ trên bề rộng đới $30\text{ km}$).
    • Kết quả này nhất quán với kết luận lịch sử của Dương Chí Công (1994) khi xác định đứt gãy Sông Hồng giai đoạn 1963–1994 có vận tốc biến dạng trượt không vượt quá $0.6\ \mu\text{rad/năm}$.

  4. Hiệu năng vượt trội của Phương pháp Xấp xỉ Sóng nhỏ so với Phương pháp Frank:

    • Phương pháp Frank (1966) bị phân mảnh theo từng tam giác rời rạc, tạo ra các bước nhảy gradient không liên tục tại ranh giới tam giác và dễ bị méo mó khi có tam giác tù hoặc dẹt.
    • Phương pháp sóng nhỏ DOG trên mặt cầu tạo ra trường biến dạng liên tục, mịn, tự động loại bỏ nhiễu biên và tái hiện chính xác các dị thường cục bộ cấp cao tại những nơi có mạng lưới quan trắc dày đặc.
                    SO SÁNH ĐẶC TÍNH NỘI SUY BIẾN DẠNG
                    
  [PHƯƠNG PHÁP FRANK (1966)]                 [PHƯƠNG PHÁP SÓNG NHỎ DOG (LUẬN ÁN)]
  - Nhảy bậc tại cạnh chung                  - Đạo hàm giải tích cấp cao trơn tru
  - Lỗi hình học tam giác dẹt                - Thích ứng đa tỷ lệ q = 7, 8, 9

Ý nghĩa thực tiễn và Định hướng ứng dụng

  • Quy hoạch công trình ngầm và hạ tầng siêu trọng: Cung cấp bản đồ phân vùng trường tốc độ biến dạng, tốc độ trượt và trương nở chi tiết làm cơ sở bắt buộc để tính toán tải trọng động đất thiết kế cho các nhà máy thủy điện bậc thang sông Đà, sông Lô, các tuyến đường sắt cao tốc và cầu dây văng nhịp lớn.
  • Quan trắc và Cảnh báo sớm tai biến địa chất: Xác định các điểm nóng tích tụ ứng suất kiến tạo tại vùng lòng chảo Điện Biên, Tuần Giáo, Mường La - Bắc Yên để ưu tiên nâng cấp mật độ trạm GNSS liên tục (CORS) phục vụ dự báo rủi ro động đất kích thích.

Hạn chế và Định hướng nghiên cứu tương lai

Nhận diện hạn chế khoa học

  1. Khoảng cách thời gian đo lặp (Episodic Campaigns): Dữ liệu dựa trên 6 chu kỳ đo định kỳ hàng năm thay vì chuỗi dữ liệu liên tục theo thời gian thực (Real-time CORS), dẫn đến việc chưa mô hình hóa trọn vẹn các thành phần phi tuyến tính theo mùa (Seasonal variations) và biến dạng sau động đất (Post-seismic relaxation).
  2. Độ chính xác thành phần độ cao ($V_U$): Do đặc thù hình học của hệ thống vệ tinh GNSS và khúc xạ khí quyển, sai số xác định vận tốc thẳng đứng ($1.5 \div 3.0\text{ mm/năm}$) vẫn còn lớn hơn đáng kể so với mặt bằng ($0.5 \div 1.0\text{ mm/năm}$).
  3. Mật độ phân bố trạm: Mật độ mốc trắc địa còn thưa ở vùng núi hiểm trở phía Đông Bắc và rẻo biên giới, làm hạn chế khả năng nâng bậc xấp xỉ lên $q = 10$ hoặc $q = 11$.

Chương trình nghiên cứu mở rộng (Future Roadmap)

  • Tích hợp giao thoa Radar đa thời gian (InSAR / PS-InSAR): Kết hợp dữ liệu tán xạ dải rộng InSAR của vệ tinh Sentinel-1/Sentinel-2 với mạng lưới GNSS sóng nhỏ để nâng cao độ phân giải không gian của trường biến dạng đứng lên cấp độ mét.
  • Mô hình hóa đàn dẻo - nhiệt động 4D: Ứng dụng sóng nhỏ không thời gian (Space-time wavelets) phân tích chuỗi thời gian liên tục tại các trạm CORS quốc gia để phát hiện các sự kiện động đất chậm (Slow Slip Events - SSE).

Đối tượng hưởng lợi từ kết quả nghiên cứu

                                  ĐỐI TƯỢNG HƯỞNG LỢI
[NHÀ NGHIÊN CỨU]     [CƠ QUAN QUẢN LÝ]   [DOANH NGHIỆP XÂY DỰNG] [VIỆN ĐỊA VẬT LÝ] [HỆ THỐNG TRẮC ĐỊA]
Cung cấp giải pháp   Bộ TN&MT sử dụng    Tính toán hệ số an  Bản đồ tích tụ ứng   Tối ưu hóa vị trí
toán học giải tích   quy hoạch đới nguy  toàn kháng chấn     suất dự báo hiểm     mốc trắc địa quốc
sóng nhỏ mặt cầu     hiểm địa chất       cho siêu công trình họa động đất         gia cấp 0
  • Nghiên cứu sinh và Giới học thuật Địa động lực: Tiếp cận một khung phương pháp luận hoàn chỉnh về giải tích sóng nhỏ mặt cầu DOG ứng dụng trong cơ học môi trường liên tục.
  • Cơ quan Quản lý Nhà nước (Bộ Tài nguyên và Môi trường, Viện Hàn lâm KH&CN Việt Nam): Cơ sở dữ liệu định lượng phục vụ công tác cập nhật bản đồ phân vùng nguy cơ động đất và hiểm họa tự nhiên quốc gia.
  • Các tổ chức Khảo sát Thiết kế Công trình Giao thông, Thủy điện: Dữ liệu chuẩn xác về tốc độ dịch chuyển vi phân dọc theo các đới đứt gãy kiến tạo cấp I và cấp II đi qua tuyến công trình.

Giải đáp cá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ì?

Đóng góp độc đáo nhất là việc thiết lập thành công mô hình xấp xỉ giải tích trường vận tốc chuyển dịch 3D trên không gian cầu $T^2(S^2)$ thông qua hàm sóng nhỏ DOG đa phân giải ($q=7, 8, 9$). Công trình này lần đầu tiên tại Việt Nam giải quyết trọn vẹn bài toán nội suy trường biến dạng mà không cần giả định môi trường biến dạng phẳng Euclidean, khắc phục triệt để sự mất mát thông tin độ cong Trái Đất trong các bài toán quy mô khu vực rộng lớn ($1.7358 \times 10^{11} \text{ m}^2$).

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

So với phương pháp mạng lưới tam giác Frank (1966) và phương pháp giải tích Spline phẳng (Haines et al., 1998), nghiên cứu của NCS Lại Văn Thủy đã:

  • Sử dụng thuật toán chuẩn hóa Tikhonov kết hợp kiểm chứng chéo thông thường (Ordinary Cross-Validation - OCV) để tối ưu hóa tham số làm trơn sóng nhỏ, loại trừ hiệu quả sai số khung quy chiếu và nhiễu ngẫu nhiên.
  • Tính toán đồng thời 4 tham số biến dạng không gian gồm: Tensor tốc độ xoay $W$, tensor tốc độ biến dạng $D$, tốc độ trượt cắt $\dot{\gamma}$ và tốc độ trương nở thể tích $\Theta$ trực tiếp từ đạo hàm giải tích của hàm sóng nhỏ cầu.

3. Phát hiện bất ngờ nhất từ tập dữ liệu thực nghiệm là gì?

Mặc dù đới đứt gãy Sông Hồng (Red River Fault) là ranh giới kiến tạo thạch quyển quan trọng bậc nhất Đông Nam Á ngăn cách khối Nam Trung Hoa và Sundaland, kết quả xử lý 6 chu kỳ đo cho thấy đới đứt gãy này hiện đang ở trạng thái hoạt động trượt bằng rất chậm ($<1.5\text{ mm/năm}$) và tốc độ biến dạng thấp ($<0.8 \times 10^{-8}/\text{năm}$). Ngược lại, năng lượng biến dạng kiến tạo hiện đại lại tập trung mật độ cao nhất tại đới đứt gãy Lai Châu - Điện Biên ở phía Tây Bắc với tốc độ biến dạng trượt đạt $(1.5 \div 3.2) \times 10^{-8}/\text{năm}$.

4. Luận án có cung cấp quy trình và công cụ để tái lập nghiên cứu không?

Có. Tác giả đã công bố sơ đồ quy trình tính toán chuẩn hóa 4 bước và xây dựng hoàn chỉnh bộ mã nguồn/phần mềm trên nền tảng ngôn ngữ lập trình Matlab. Phần mềm bao gồm các giao diện xử lý trực quan: module tính toán véc-tơ Euler, module nội suy phân tích sóng nhỏ DOG và module dẫn xuất tensor biến dạng, cho phép các cơ quan nghiên cứu độc lập tái lập và kiểm chứng quy trình trên các tập dữ liệu trắc địa khác.

5. Lộ trình phát triển 10 năm của hướng nghiên cứu này được xác lập ra sao?

  1. Mở rộng khung tính toán sóng nhỏ tích hợp đồng thời mạng lưới GNSS liên tục (VNGEONET / CORS) và chuỗi dữ liệu giao thoa viễn thám InSAR đa băng tần.
  2. Ứng dụng biến đổi sóng nhỏ thời gian - không gian (Spatio-temporal wavelets) để giải mã các tín hiệu biến dạng phi tuyến, phân tách biến dạng kiến tạo nội sinh và biến dạng lún sụt cục bộ do nhân sinh (khai thác nước ngầm, đô thị hóa).

Kết luận

Luận án tiến sĩ của Nghiên cứu sinh Lại Văn Thủy là công trình khoa học mẫu mực, hoàn thiện xuất sắc cả về mặt lý thuyết toán học trắc địa lẫn thực nghiệm địa động lực học công trình.

                      5 ĐÓNG GÓP BẢO VỆ CỐT LÕI CỦA LUẬN ÁN
[LÝ THUYẾT MỚI]                       [CÔNG CỤ PHẦN MỀM]                    [BẢN ĐỒ BIẾN DẠNG 3D]
Chuyển đổi thành công giải tích       Xây dựng bộ công cụ Matlab tự động    Thành lập hệ thống bản đồ
sóng nhỏ DOG trên mặt cầu S² cho      hóa: Euler, Wavelet DOG, OCV,         chuẩn 4 trường biến dạng khu
trường biến dạng trắc địa 3D.         Tensor Gradient L = D + W.            vực Miền Bắc (2012 - 2017).
        [MINH GIẢI KIẾN TẠO]                [ỨNG DỤNG THỰC TIỄN]
        Làm sáng tỏ hoạt động tích tụ       Cung cấp luận cứ khoa học tin cậy
        ứng suất tại đới Tây Bắc            phục vụ phòng chống rủi ro thiên
        (Lai Châu - Điện Biên).             tai và xây dựng công trình ngầm.
  1. Xác lập khung lý thuyết mới: Đặt nền móng vững chắc cho việc ứng dụng giải tích đa phân giải sóng nhỏ mặt cầu trong trắc địa địa động lực tại Việt Nam.
  2. Làm chủ giải pháp công nghệ: Xây dựng thành công gói phần mềm tính toán nội suy sóng nhỏ chuyên dụng trên nền tảng Matlab, giải quyết bài toán xử lý mạng lưới trắc địa phân bố không đều.
  3. Thành lập hệ thống bản đồ biến dạng 3D chuẩn hóa: Cung cấp bộ bản đồ số hóa trường tốc độ xoay, tốc độ biến dạng, tốc độ trượt và tốc độ trương nở cho toàn bộ miền Bắc Việt Nam.
  4. Minh giải bản chất động lực học khu vực: Chứng minh định lượng rằng hoạt động đứt gãy hiện đại ở miền Bắc có tốc độ dịch chuyển tương đối nhỏ, tập trung tích tụ biến dạng lớn nhất ở khu vực đứt gãy Lai Châu - Điện Biên.
  5. Đóng góp giá trị thực tiễn to lớn: Kết quả nghiên cứu là tài liệu nền tảng phục vụ công tác quy hoạch an toàn lãnh thổ, cảnh báo tai biến địa chất và thiết kế kháng chấn cho các công trình trọng điểm quốc gia.