Nghiên cứu xác định dị thường độ cao bằng dữ liệu mặt đất và vệ tinh tại Miền Trung Việt Nam

Tổng quan về luận án

Luận án tiến sĩ "Nghiên cứu xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh, áp dụng cho khu vực miền Trung Việt Nam" của Nguyễn Thành Lê đóng góp một cách đột phá vào lĩnh vực trắc địa vật lý, đặc biệt trong bối cảnh địa hình phức tạp và hạn chế về dữ liệu tại Việt Nam. Nghiên cứu này xác lập một phương pháp luận khoa học chặt chẽ và quy trình ứng dụng thực tiễn để xác định dị thường độ cao với độ chính xác cao, giải quyết một nhiệm vụ then chốt cho việc xây dựng mô hình Quasigeoid quốc gia và ứng dụng công nghệ GNSS trong đo độ cao chuẩn.

Bối cảnh khoa học và tính tiên phong của nghiên cứu: Việc xác định dị thường độ cao và xây dựng mô hình Quasigeoid là một trong những nhiệm vụ trọng tâm của trắc địa, có ý nghĩa quan trọng đối với mỗi quốc gia. Nhờ sự phát triển của công nghệ GNSS, độ cao trắc địa có thể được xác định với độ chính xác cao, nhưng việc chuyển đổi sang độ cao chuẩn vẫn còn là thách thức lớn, đặc biệt ở những khu vực địa hình đồi núi, sông ngòi hiểm trở nơi việc đo thủy chuẩn truyền thống gặp khó khăn. Luận án tiên phong trong việc tích hợp một cách toàn diện các nguồn dữ liệu đa dạng – từ mặt đất đến vệ tinh – để khắc phục những hạn chế này, tạo ra một giải pháp đồng bộ và hiệu quả cho khu vực có đặc trưng địa lý phức tạp như miền Trung Việt Nam.

Research Gap SPECIFIC với citations từ literature: Các nghiên cứu trước đây ở Việt Nam chủ yếu tập trung vào "thử nghiệm các mô hình thế trọng trường Trái Đất đã có cho phù hợp với lãnh thổ Việt Nam hoặc cải tiến mô hình thế trọng trường Trái Đất hiện có bằng cách hiệu chỉnh các giá trị hiện chỉnh vào các hệ số của mô hình." (Kết luận và kiến nghị). Đặc biệt, "Các công trình về nghiên cứu xác định dị thường độ cao, xây dựng mô hình Quasigeoid trên cở sở kết hợp đầy đủ các loại dữ liệu mặt đất và vệ tinh hiện nay còn hạn chế." (Kết luận và kiến nghị). Cụ thể, Phạm Hoàng Lân (2009) từng gặp hạn chế về dữ liệu trọng lực mặt đất (mắt lưới 5'x5' khá lớn, chỉ hơn 2000 trong số 4000 ô có trị đo trực tiếp) và thiếu dữ liệu vùng xa, dẫn đến kết quả kém ý nghĩa về độ chính xác. Vũ Hồng Cường (2013) đã nghiên cứu mô hình Quasigeoid từ dữ liệu vệ tinh nhưng "dữ liệu trọng lực đo trực tiếp trên mặt đất là không được sử dụng mà chỉ sử dụng dữ liệu trọng lực từ mô hình DTU và bản đồ trọng lực toàn cầu WGM2012". Lê Minh (2008-2011) đạt độ chính xác ±0.082 m cho Geoid địa phương toàn quốc nhưng chỉ đạt 0.08-0.10 m ở vùng núi cao do thiếu dữ liệu trọng lực chi tiết. Những hạn chế này tạo ra một "phương pháp và quy trình xác định dị thường độ cao, xây dựng mô hình Quasigeoid trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh phù hợp với điều kiện dữ liệu hiện có tại Việt Nam chưa hoàn thiện." (Kết luận và kiến nghị). Luận án này lấp đầy khoảng trống đó bằng cách đề xuất một quy trình toàn diện, kết hợp các nguồn dữ liệu đa dạng và xử lý ảnh hưởng của dữ liệu vùng xa, địa hình phức tạp một cách hiệu quả.

Research questions và hypotheses: Mặc dù luận án không nêu rõ dưới dạng câu hỏi nghiên cứu, mục tiêu đề tài và các luận điểm bảo vệ đã định hướng rõ ràng.

• Mục tiêu 1: Xác lập cơ sở khoa học và phương pháp luận xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh nói chung và thực nghiệm tính toán thành công cho khu vực miền Trung Việt Nam nói riêng.

• Mục tiêu 2: Lựa chọn cơ sở toán học (hệ quy chiếu, hệ tọa độ, hệ triều) thống nhất cho các loại dữ liệu từ các nguồn khác nhau, đảm bảo thuận tiện trong tính toán và đánh giá.

• Mục tiêu 3: Xác định dị thường độ cao cho khu vực miền Trung Việt Nam từ nguồn dữ liệu hạn chế, thành lập bản đồ dị thường độ cao trong hệ WGS84 quốc tế.

• Luận điểm bảo vệ 1: Kết hợp kỹ thuật “loại bỏ–tính toán–phục hồi” (RCR) và phương pháp Collocation bình phương nhỏ nhất (LSC) là giải pháp phù hợp cho phép xác định dị thường độ cao từ các nguồn dữ liệu khác nhau của khu vực địa hình đặc trưng miền Trung Việt Nam.

• Luận điểm bảo vệ 2: Quy trình và module chương trình tính dị thường độ cao từ sự kết hợp các dữ liệu khác nhau cho phép nâng cao độ chính xác, mức độ tự động hóa và thuận lợi ứng dụng trong điều kiện thực tế ở Việt Nam.

Theoretical framework với tên theories cụ thể: Luận án được xây dựng trên nền tảng lý thuyết Trắc địa vật lý cổ điển của Molodensky (1962) về bề mặt thực và thế trọng trường bên ngoài của Trái Đất, thay thế cho lý thuyết Stokes truyền thống trong việc xác định Geoid chính xác (Heiskanen, 1967). Khung lý thuyết trọng tâm là kỹ thuật Remove-Compute-Restore (RCR), được phát triển bởi Forsberg (1981, 1984), kết hợp với phương pháp Collocation bình phương nhỏ nhất (LSC) của Tscherning (1974, 2015). Các mô hình thế trọng trường Trái Đất toàn cầu (Global Geopotential Models - GGM) như EIGEN-6C4 (Cộng đồng Châu Âu, 2014) và mô hình địa hình phần dư (Residual Terrain Model - RTM) đóng vai trò then chốt trong việc phân tách và xử lý các thành phần bước sóng của thế nhiễu. Luận án cũng sử dụng các công thức chuyển đổi hệ tọa độ Bursa-Wolf (1962) và các nguyên tắc thống kê xác suất để phát hiện sai số thô (Nguyễn Cao Văn et al., 2018).

Đóng góp đột phá với quantified impact:

1. Hoàn thiện phương pháp luận kết hợp dữ liệu đa nguồn: Luận án đã phát triển và hoàn thiện phương pháp kết hợp kỹ thuật RCR và LSC để xác định dị thường độ cao từ "đa dạng các nguồn dữ liệu mặt đất và vệ tinh" bao gồm dữ liệu trọng lực mặt đất, trọng lực vệ tinh (từ các dự án CHAMP, GRACE, GOCE qua EGM), dữ liệu đo cao vệ tinh, mô hình số địa hình độ phân giải cao, và điểm song trùng GNSS-Thủy chuẩn. Phương pháp này "tường minh về mặt lý thuyết" và "thể hiện đúng đắn về thực nghiệm" trong điều kiện địa hình phức tạp của miền Trung Việt Nam.
2. Tự động hóa và nâng cao hiệu quả quy trình tính toán: Luận án đã xây dựng "quy trình và module chương trình xác định dị thường độ cao... có tính tự động hóa mạnh" thay thế các công cụ hiện có như GRAVSOFT vốn "đòi hỏi người dùng phải có trình độ cao" và chưa tự động hóa. Điều này giảm đáng kể thời gian và công sức, đồng thời tăng độ tin cậy của kết quả.
3. Đạt độ chính xác đột phá cho khu vực miền Trung Việt Nam: Kết quả thực nghiệm đạt "độ lệch chuẩn đạt ±0,078 m" khi kết hợp đầy đủ dữ liệu trọng lực, EIGEN-6C4, dữ liệu địa hình và GNSS-Thủy chuẩn (PA3). Đây là một cải thiện đáng kể so với việc đánh giá trực tiếp từ EIGEN-6C4 (±0,194 m) và các phương án kết hợp ít dữ liệu hơn (PA1: ±0,183 m; PA2: ±0,096 m). Độ chính xác này vượt trội so với các nghiên cứu trước đây tại Việt Nam ở khu vực địa hình khó khăn (ví dụ, Lê Minh chỉ đạt 0.08-0.10m ở vùng núi).
4. Khả năng ứng dụng thực tiễn cao: Kết quả luận án cung cấp "bản đồ dị thường độ cao trọng lực cho khu vực thực nghiệm trong hệ quy chiếu và hệ tọa độ WGS84 quốc tế," với mắt lưới 3'x3' (tương đương khoảng 5.5 km x 5.5 km) và khoảng cao đều 0.25 m. Đây là cơ sở quan trọng cho phép "xác định độ cao chuẩn bằng công nghệ đo cao GNSS" một cách hiệu quả, góp phần trực tiếp vào công tác trắc địa và quy hoạch tại Việt Nam.

Scope (sample size, timeframe) và significance:

• Phạm vi không gian: Khu vực miền Trung Việt Nam (vĩ độ 15° ÷ 20° B, kinh độ 104° ÷ 109° Đ), bao gồm các tỉnh từ Nghệ An đến Quảng Nam - Đà Nẵng, một phần Thanh Hóa, và phần Biển Đông giáp Lào ở phía Tây. Khu vực này có địa hình đa dạng, phức tạp (đồi núi, đồng bằng hẹp, ven biển), đại diện cho đặc điểm địa hình chung của Việt Nam.
• Dữ liệu sử dụng:
  • Trọng lực mặt đất: 9.897 điểm dị thường trọng lực không khí tự do (được lọc từ 9.997 điểm ban đầu sau khi loại bỏ 100 điểm sai số thô), trong hệ WGS84, hệ triều không phụ thuộc triều.
  • Trọng lực biển: 1.034 điểm dị thường trọng lực không khí tự do từ mô hình DTU17GRAV (mắt lưới 5'x5').
  • Mô hình thế trọng trường Trái Đất: EIGEN-6C4 (Nmax = 2190).
  • Mô hình số địa hình (DTM): Mô hình tổng hợp SRTM3arc_v4.1 (3"x3") cho đất liền và GEBCO 2020 (15"x15") cho biển.
  • Điểm song trùng GNSS-Thủy chuẩn: 131 điểm (101 điểm dùng để tính toán và làm khớp, 30 điểm dùng để đánh giá độc lập).
• Significance: Luận án cung cấp một giải pháp toàn diện và thực tế cho việc xác định dị thường độ cao, một yếu tố then chốt để chuyển đổi độ cao trắc địa GNSS sang độ cao chuẩn quốc gia. Điều này không chỉ nâng cao độ chính xác cho công tác trắc địa mà còn hỗ trợ quy hoạch phát triển kinh tế-xã hội, quản lý tài nguyên và môi trường, cũng như các ứng dụng kỹ thuật công trình đòi hỏi độ cao chính xác cao trong khu vực miền Trung Việt Nam và có tiềm năng mở rộng ra toàn quốc.

Literature Review và Positioning

Synthesis của major streams với TÊN TÁC GIẢ và NĂM cụ thể: Nghiên cứu dị thường độ cao và mô hình Geoid/Quasigeoid là một lĩnh vực được quan tâm rộng rãi trên thế giới. Các quốc gia như Slovakia (1999, mô hình GMSQ98A, độ chính xác ±2.4 cm), Châu Âu (2008, EGG07, đạt ±0.01-±0.02 m), Hoa Kỳ (2009, USGG2009, đạt ±4.3 cm toàn quốc), Canada (2010, CGG2010, đạt 2-10 cm), Ba Lan (2013, GDQM13, đạt ±1.4-±1.8 cm), Australia (2017, AGQG2017, đạt ±5-±10 cm), New Zealand (2017, NZGeoid2017, đạt ±4.8 cm) và Iran (2018, IRG2018, đạt ±20 cm) đều đã đạt được những thành tựu đáng kể. Điểm nổi bật là việc sử dụng dữ liệu trọng lực mặt đất dày đặc, dữ liệu trọng lực vệ tinh từ các dự án CHAMP, GRACE, GOCE, và các mô hình thế trọng trường Trái Đất (EGM2008, EIGEN-6C4). Các phương pháp chính được áp dụng là kỹ thuật Remove-Compute-Restore (RCR) và Collocation bình phương nhỏ nhất (LSC), thường được thực hiện thông qua phần mềm GRAVSOFT (Tscherning & Forsberg, 2008). Việc làm khớp và đánh giá kết quả luôn dựa trên các điểm GNSS-Thủy chuẩn.

Tại Việt Nam, các nghiên cứu cũng đã được triển khai nhưng còn đối mặt với nhiều thách thức. Phạm Hoàng Lân (2009) đã nghiên cứu thành lập hệ độ cao chuẩn thống nhất nhưng gặp khó khăn về số liệu trọng lực mặt đất thưa thớt và thiếu dữ liệu vùng xa. Lê Minh (2005) xây dựng mô hình Geoid địa phương cho lãnh thổ Việt Nam, đạt độ chính xác ±0.082 m tổng thể, nhưng ở vùng núi cao chỉ đạt 0.08-0.10 m do thiếu dữ liệu trọng lực chi tiết. Nguyễn Duy Đô (2012) tập trung chính xác hóa EGM2008 cục bộ bằng dữ liệu GPS-Thủy chuẩn, chỉ ra rằng mô hình EGM2008 phù hợp cho bước sóng lớn (>25 km) nhưng cần hiệu chỉnh địa hình cục bộ cho bước sóng nhỏ hơn. Vũ Hồng Cường (2013) tiên phong ứng dụng dữ liệu gradient trọng lực vệ tinh GOCE cho Quasigeoid ở Việt Nam, đạt độ chính xác 9-10 cm ở miền Trung và Nam, nhưng lại không sử dụng dữ liệu trọng lực đo trực tiếp trên mặt đất. Các tác giả khác như Phạm Thị Hoa (2012) và Bùi Thị Hồng Thắm (2015) đã chứng minh tầm quan trọng của việc xác định ảnh hưởng của địa hình trong dị thường độ cao, đặc biệt ở vùng núi.

Contradictions/debates với ít nhất 2 opposing views: Một trong những tranh luận chính trong lĩnh vực này là mức độ tin cậy của các mô hình thế trọng trường toàn cầu (GGM) và dữ liệu vệ tinh so với dữ liệu trọng lực mặt đất, đặc biệt ở các khu vực có địa hình phức tạp.

• View 1: Ưu tiên dữ liệu vệ tinh/GGM cho bước sóng dài: Các nghiên cứu quốc tế và tại Việt Nam đều khẳng định tầm quan trọng của các mô hình GGM (ví dụ, EGM2008, EIGEN-6C4) được xây dựng từ dữ liệu vệ tinh CHAMP, GRACE, GOCE trong việc phản ánh các thành phần bước sóng dài của trường trọng lực (Rummel, 2002; Nguyễn Văn Sáng, 2013). Các mô hình này cung cấp một nền tảng toàn cầu với độ chính xác cao ở quy mô lớn.
• View 2: Dữ liệu mặt đất và địa hình là không thể thiếu cho bước sóng ngắn và trung: Tuy nhiên, các nghiên cứu cũng chỉ ra rằng GGM không đủ để phản ánh chính xác các thành phần bước sóng ngắn và trung, đặc biệt ở các khu vực có địa hình phức tạp. Forsberg (1984) nhấn mạnh rằng "hiệu ứng làm mịn dữ liệu trọng lực có thể đạt 50% nếu có được mô hình địa hình độ chính xác cao." (Mục 2.3.3). Các tác giả như Phạm Thị Hoa (2012) và Bùi Thị Hồng Thắm (2015) cũng đã chỉ rõ ảnh hưởng đáng kể của độ cao địa hình ở Việt Nam. Ngay cả Vũ Hồng Cường (2013), mặc dù sử dụng dữ liệu vệ tinh, vẫn đối mặt với độ chính xác dao động và không thể thay thế hoàn toàn nhu cầu về dữ liệu mặt đất chi tiết. Luận án này khẳng định rằng "dữ liệu trọng lực được đo trực tiếp là không thể thiếu trong xác định dị thường độ cao từ dữ liệu kết hợp mặt đất và vệ tinh." (Kết luận và kiến nghị).

Positioning trong literature với specific gap identified: Luận án này tự định vị là một công trình quan trọng nhằm giải quyết trực tiếp khoảng trống về phương pháp luận và quy trình tích hợp toàn diện dữ liệu mặt đất và vệ tinh để xác định dị thường độ cao tại Việt Nam. Nó vượt qua hạn chế của các nghiên cứu trước đây vốn hoặc thiếu dữ liệu mặt đất chi tiết (Phạm Hoàng Lân, 2009; Vũ Hồng Cường, 2013), hoặc chưa tối ưu hóa việc kết hợp đa nguồn dữ liệu để đạt độ chính xác cao trong điều kiện địa hình phức tạp (Lê Minh, 2005; Nguyễn Duy Đô, 2012). Cụ thể, luận án giải quyết vấn đề "phương pháp và quy trình xác định dị thường độ cao, xây dựng mô hình Quasigeoid trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh phù hợp với điều kiện dữ liệu hiện có tại Việt Nam chưa hoàn thiện" (Kết luận và kiến nghị).

How this advances field với concrete contributions: Luận án tiến thêm một bước đáng kể bằng cách thiết lập một quy trình chuẩn hóa và tự động hóa cao, kết hợp kỹ thuật RCR và LSC với các module chương trình máy tính tự phát triển. Điều này cho phép xử lý và tích hợp hiệu quả các nguồn dữ liệu đa dạng (trọng lực mặt đất, trọng lực biển DTU17GRAV, mô hình EIGEN-6C4, DTM tổng hợp SRTM/GEBCO, và GNSS-Thủy chuẩn) một cách tường minh về lý thuyết và chính xác về thực nghiệm, đặc biệt trong điều kiện dữ liệu hạn chế và địa hình phức tạp của miền Trung Việt Nam. Các module tự động hóa giúp giải quyết các thách thức về định dạng dữ liệu và xác định tham số hàm hiệp phương sai lý thuyết, vốn là những rào cản lớn khi sử dụng các phần mềm truyền thống như GRAVSOFT.

So sánh với ÍT NHẤT 2 international studies:

• So với Slovakia (GMSQ98A, 1999) và Ba Lan (GDQM13, 2013): Các quốc gia này đã đạt độ chính xác cao (±2.4 cm và ±1.4-±1.8 cm tương ứng) cho mô hình Quasigeoid của mình. Tuy nhiên, họ đều có "dữ liệu trọng lực mặt đất dày đặc và chi tiết với độ chính xác cao, mắt lưới được thực hiện với kích thước 1′×1′." (Tổng quan về nghiên cứu xác định dị thường độ cao trên thế giới). Luận án của Nguyễn Thành Lê, dù với mật độ dữ liệu mặt đất hạn chế hơn ở khu vực núi (chủ yếu tập trung ở đồng bằng), vẫn đạt được độ chính xác ±0.078 m ở khu vực miền Trung Việt Nam với mắt lưới 3'x3'. Điều này chứng tỏ hiệu quả của phương pháp tích hợp dữ liệu và kỹ thuật RCR/LSC trong việc giảm ảnh hưởng của dữ liệu vùng xa và lấp đầy khoảng trống dữ liệu, cho phép đạt được độ chính xác cao ngay cả trong điều kiện kém lý tưởng.
• So với New Zealand (NZGeoid2017, 2017): New Zealand đã đạt độ chính xác ±4.8 cm cho mô hình Quasigeoid của mình bằng cách kết hợp dữ liệu trọng lực hàng không trên toàn lãnh thổ và dữ liệu GOCE. Trong khi đó, luận án này chưa có dữ liệu trọng lực hàng không mà chủ yếu dựa vào dữ liệu trọng lực mặt đất (9.897 điểm) và dữ liệu trọng lực biển (1.034 điểm từ DTU17GRAV). Mặc dù vậy, độ chính xác ±0.078 m đạt được vẫn là một thành tựu đáng kể cho khu vực miền Trung Việt Nam, cho thấy phương pháp tích hợp đa nguồn dữ liệu của luận án có khả năng hoạt động hiệu quả ngay cả khi thiếu một số loại dữ liệu tiên tiến khác. Hơn nữa, kết quả này được đánh giá dựa trên các điểm GNSS-Thủy chuẩn độc lập, đảm bảo tính khách quan.

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

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

Luận án này đã mở rộng và thách thức một số lý thuyết hiện có trong lĩnh vực trắc địa vật lý thông qua việc phát triển một khung phân tích tích hợp.

• Extend/challenge WHICH specific theories (name theorists):

  • Mở rộng lý thuyết Molodensky (1962) về Quasigeoid: Luận án đã củng cố và mở rộng việc xác định Quasigeoid bằng cách tích hợp các thành phần bước sóng khác nhau của thế nhiễu từ đa dạng nguồn dữ liệu, vượt qua sự phức tạp trong việc xác định Geoid trực tiếp (Heiskanen, 1967).
  • Thách thức các giới hạn của lý thuyết Stokes trong điều kiện địa hình phức tạp: Trong khi công thức Stokes là nền tảng cho việc xác định Geoid, nó thường giả định các dị thường trọng lực đã được khử ảnh hưởng địa hình. Luận án, bằng cách áp dụng kỹ thuật RCR và phương pháp RTM (Forsberg, 1984), đã chứng minh một cách tường minh rằng việc loại bỏ triệt để ảnh hưởng địa hình (bao gồm cả bước sóng ngắn và trung) là tối cần thiết cho việc nội suy dị thường độ cao chính xác, đặc biệt ở miền Trung Việt Nam với địa hình "phức tạp, hiểm trở" (Mục 3.1).
  • Tối ưu hóa phương pháp Least-Squares Collocation (LSC) (Tscherning, 1974): Luận án đã phát triển các module chương trình tự động hóa việc xác định các tham số tối ưu (A, N, a, RB-RE) của hàm hiệp phương sai lý thuyết khi làm khớp với hàm thực nghiệm, giải quyết vấn đề "không phải N nào cũng tìm được các tham số do có những giá trị của N làm cho chương trình bị suy biến" và sự tốn kém thời gian của việc tính toán thủ công với Nmax lên đến 2190 (Mục 2.7.5). Điều này nâng cao hiệu quả và độ tin cậy của LSC trong ứng dụng thực tiễn.

• Conceptual framework với components và relationships: Khung khái niệm của luận án xoay quanh việc phân tách thế nhiễu (disturbing potential) thành các thành phần bước sóng khác nhau, mỗi thành phần được phản ánh bởi một loại dữ liệu cụ thể:

  • Bước sóng dài (Long-wave components): Được phản ánh chủ yếu bởi Mô hình thế trọng trường Trái Đất (EGM) toàn cầu như EIGEN-6C4.
  • Bước sóng trung và ngắn (Medium and short-wave components): Được phản ánh bởi dữ liệu trọng lực đo trực tiếp trên mặt đất và biển (DTU17GRAV), cũng như ảnh hưởng của địa hình được xử lý thông qua mô hình địa hình phần dư (RTM).
  • Dữ liệu khống chế và đánh giá: Điểm song trùng GNSS-Thủy chuẩn cung cấp dị thường độ cao có độ chính xác cao nhất để làm khớp và đánh giá kết quả. Mối quan hệ giữa các thành phần này được quản lý thông qua kỹ thuật RCR (Remove-Compute-Restore). Trong đó, các thành phần bước sóng dài và trung/ngắn được "loại bỏ" khỏi dữ liệu gốc, phần dư được "tính toán" bằng LSC, và sau đó các thành phần đã loại bỏ được "phục hồi" vào kết quả cuối cùng để có dị thường độ cao tổng hợp.

• Theoretical model với propositions/hypotheses numbered: Luận án đề xuất một mô hình lý thuyết thực nghiệm, được cụ thể hóa qua ba phương án tính toán (PA1, PA2, PA3) với các giả định và kết quả được định lượng:

  • Proposition 1 (PA1): Dị thường độ cao có thể được xác định bằng cách kết hợp trực tiếp dữ liệu mô hình EIGEN-6C4 và dữ liệu GNSS-Thủy chuẩn. Kết quả thực nghiệm cho thấy độ lệch chuẩn ±0.183 m.
  • Proposition 2 (PA2): Việc bổ sung dữ liệu trọng lực mặt đất và loại bỏ thành phần bước sóng dài từ EIGEN-6C4 sẽ nâng cao độ chính xác hơn nữa. Kết quả thực nghiệm chứng minh độ lệch chuẩn giảm xuống ±0.096 m. (Cải thiện rõ rệt so với PA1, khẳng định hiệu quả của việc tích hợp dữ liệu trọng lực).
  • Proposition 3 (PA3): Việc tích hợp đầy đủ dữ liệu trọng lực, mô hình EIGEN-6C4, dữ liệu địa hình (qua RTM), và GNSS-Thủy chuẩn, đồng thời loại bỏ cả thành phần bước sóng dài và bước sóng trung/ngắn, sẽ mang lại độ chính xác tối ưu nhất. Kết quả thực nghiệm xác nhận độ lệch chuẩn đạt ±0.078 m. (Mức chính xác cao nhất, chứng minh luận điểm 1 và 2 của luận án).

• Paradigm shift với EVIDENCE từ findings: Luận án không đề xuất một sự thay đổi hoàn toàn về mô hình khoa học (paradigm shift) theo nghĩa của Kuhn, mà là một sự tiến bộ đáng kể trong khuôn khổ thực nghiệm (positivist/post-positivist) trong trắc địa vật lý. Bằng chứng từ các phát hiện cho thấy sự chuyển dịch từ việc "chủ yếu là thử nghiệm các mô hình thế trọng trường Trái Đất đã có" sang việc chủ động "nghiên cứu thiết lập được cơ sở khoa học xác định dị thường độ cao trên cơ sở kết hợp đa dạng các nguồn dữ liệu" (Điểm mới của luận án). Việc đạt được độ chính xác ±0.078 m trong điều kiện dữ liệu hạn chế và địa hình phức tạp, thông qua một quy trình được tự động hóa, đã nâng cao khả năng ứng dụng thực tế và độ tin cậy của việc xác định độ cao chuẩn bằng GNSS ở Việt Nam.

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

• Integration của theories (name 3+ specific theories): Khung phân tích độc đáo của luận án là sự tích hợp sâu rộng các lý thuyết và kỹ thuật nền tảng:

  1. Lý thuyết thế nhiễu và Quasigeoid của Molodensky (1962): Cung cấp cơ sở hình học và vật lý cho việc xác định bề mặt chuẩn.
  2. Kỹ thuật Remove-Compute-Restore (RCR) (Forsberg, 1981, 1984): Là xương sống để phân tách và xử lý các thành phần bước sóng khác nhau của trường trọng lực.
  3. Phương pháp Collocation bình phương nhỏ nhất (LSC) (Tscherning, 1974): Cung cấp công cụ tối ưu để nội suy các phần dư dị thường độ cao từ dữ liệu phần dư trọng lực và GNSS-Thủy chuẩn, dựa trên nguyên tắc hiệp phương sai.
  4. Mô hình địa hình phần dư (RTM): Giải quyết ảnh hưởng cục bộ của địa hình một cách hiệu quả hơn so với hiệu chỉnh Faye truyền thống, đặc biệt ở vùng núi.
  5. Lý thuyết hàm điều hòa cầu: Để mô tả và xử lý các mô hình thế trọng trường toàn cầu (EGM).

• Novel analytical approach với justification: Luận án đã phát triển một "quy trình xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh từ các nguồn dữ liệu hạn chế tại miền Trung Việt Nam" (Điểm mới của luận án). Tính mới nằm ở việc:

  • Thống nhất cơ sở toán học: Đảm bảo tất cả dữ liệu được quy chiếu về "hệ quy chiếu và hệ tọa độ WGS84 quốc tế, hệ triều thuộc hệ triều không phụ thuộc triều" (Mục 2.2.3, 2.2.4), giải quyết vấn đề đa dạng về hệ thống tham chiếu của các nguồn dữ liệu khác nhau.
  • Quy trình loại bỏ sai số thô nghiêm ngặt: Áp dụng quy tắc 3-sigma để phát hiện và loại bỏ các trị đo bất thường trong dữ liệu dị thường trọng lực và GNSS-Thủy chuẩn, nâng cao chất lượng dữ liệu đầu vào.
  • Module phần mềm tự động hóa: Phát triển các module chương trình tự động hóa các bước tính toán phức tạp như xác định phương sai mức, ảnh hưởng địa hình RTM, và khớp hàm hiệp phương sai lý thuyết, giảm sự phụ thuộc vào kinh nghiệm thủ công và phần mềm thương mại có cấu trúc cứng nhắc như GRAVSOFT. Tính hợp lý của cách tiếp cận này được chứng minh bằng kết quả thực nghiệm: PA3, tích hợp đầy đủ các bước loại bỏ và phục hồi, đạt độ chính xác cao nhất (±0.078 m), khẳng định rằng việc xử lý triệt để các thành phần bước sóng và chuẩn hóa dữ liệu là chìa khóa để đạt được kết quả tối ưu trong điều kiện thực tiễn.

• Conceptual contributions với definitions:

  • Quy trình kết hợp dữ liệu đa nguồn tối ưu: Là một chuỗi các bước được thiết kế để tích hợp các loại dữ liệu địa lý khác nhau (trọng lực, địa hình, vệ tinh, GNSS-Thủy chuẩn) nhằm xác định dị thường độ cao, có tính đến các đặc trưng riêng của từng loại dữ liệu và ảnh hưởng qua lại của chúng.
  • Module tự động hóa cho Collocation: Các công cụ phần mềm được lập trình để tự động thực hiện các phép tính phức tạp trong phương pháp LSC, bao gồm xác định hiệp phương sai thực nghiệm và khớp hàm hiệp phương sai lý thuyết, giảm thiểu sự can thiệp thủ công và nâng cao hiệu quả tính toán.
  • Khái niệm về "dữ liệu hạn chế" được xử lý hiệu quả: Luận án chứng minh rằng ngay cả với "nguồn dữ liệu hạn chế" (ví dụ: dữ liệu trọng lực mặt đất không đồng đều, thiếu dữ liệu vùng xa từ các nước láng giềng), một phương pháp luận khoa học chặt chẽ vẫn có thể đạt được kết quả có độ chính xác cao.

• Boundary conditions explicitly stated:

  • Phạm vi không gian: Nghiên cứu tập trung vào khu vực miền Trung Việt Nam (vĩ độ 15° ÷ 20° B, kinh độ 104° ÷ 109° Đ). Các kết quả và độ chính xác đạt được là tối ưu trong phạm vi này.
  • Loại dữ liệu: Hiệu quả của phương pháp phụ thuộc vào sự sẵn có và chất lượng của các loại dữ liệu đầu vào cụ thể (trọng lực mặt đất, trọng lực biển, GGM, DTM độ phân giải cao, điểm GNSS-Thủy chuẩn).
  • Độ phân giải: Bản đồ dị thường độ cao được thành lập có mắt lưới 3'x3', phù hợp với yêu cầu xác định độ cao chuẩn bằng công nghệ GNSS trong khu vực thực nghiệm.
  • Mô hình thế trọng trường: Việc lựa chọn mô hình EIGEN-6C4 làm mô hình tham chiếu cho bước sóng dài là kết quả của quá trình đánh giá kỹ lưỡng các GGM phù hợp nhất với khu vực Việt Nam (Mục 2.4.2).

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

Thiết kế nghiên cứu

Luận án này áp dụng một thiết kế nghiên cứu thực nghiệm chặt chẽ, được định hướng bởi triết lý khoa học thực chứng (positivist/post-positivist), tập trung vào việc phát triển và kiểm chứng một phương pháp định lượng.

• Research philosophy (positivism/interpretivism/critical realism): Triết lý nghiên cứu của luận án là Thực chứng (Positivism), thể hiện qua việc:

  • Tập trung vào "xác lập được cơ sở khoa học và phương pháp luận" (Mục tiêu đề tài), đề cao tính khách quan, đo lường được và khả năng khái quát hóa.
  • Sử dụng các phương pháp định lượng (Collocation bình phương nhỏ nhất, thống kê xác suất) và dữ liệu khách quan (trọng lực, vệ tinh, GNSS-Thủy chuẩn).
  • Mục tiêu là xây dựng "quy trình và module chương trình tính dị thường độ cao... có tính tự động hóa mạnh" (Luận điểm 2), hướng tới các kết quả có thể tái lập và kiểm chứng.
  • Đánh giá độ chính xác dựa trên "độ lệch chuẩn" và "sai số trung phương" (Mục 2.9, Bảng 3.11), là các chỉ số định lượng rõ ràng.

• Mixed methods với SPECIFIC combination rationale: Mặc dù không phải là "mixed methods" theo nghĩa kết hợp định tính-định lượng, luận án sử dụng một phương pháp kết hợp dữ liệu (data fusion) tiên tiến từ nhiều nguồn khác nhau, mỗi nguồn cung cấp thông tin cho các thành phần bước sóng riêng biệt của trường trọng lực, tạo thành một cách tiếp cận toàn diện:

  • Dữ liệu trọng lực vệ tinh (EGM EIGEN-6C4): Cung cấp thông tin cho các thành phần bước sóng dài (phạm vi từ 100 km trở lên).
  • Dữ liệu trọng lực mặt đất và biển (DTU17GRAV): Đóng góp cho các thành phần bước sóng trung (20-200 km).
  • Dữ liệu địa hình (SRTM3arc_v4.1, GEBCO 2020): Hiệu chỉnh ảnh hưởng của các thành phần bước sóng ngắn (<20 km).
  • Dữ liệu GNSS-Thủy chuẩn: Đóng vai trò là dữ liệu tham chiếu có độ chính xác cao nhất, vừa làm khớp mô hình, vừa dùng để đánh giá độ chính xác cuối cùng. Sự kết hợp này là cần thiết vì "mỗi loại dữ liệu lại phản ánh thành phần bước sóng khác nhau" (Mục 1.2), và "để xác định được dị thường độ cao chính xác cần phải kết hợp các dữ liệu với nhau" (Mục 2.3).

• Multi-level design với levels clearly defined: Thiết kế nghiên cứu mang tính đa cấp độ trong việc xử lý và tích hợp dữ liệu:

  • Cấp độ 1: Dữ liệu toàn cầu (Global Level): Sử dụng EGM EIGEN-6C4 để mô tả các thành phần bước sóng dài của trường trọng lực Trái Đất.
  • Cấp độ 2: Dữ liệu khu vực/quốc gia (Regional/National Level): Tích hợp dữ liệu trọng lực mặt đất (9.897 điểm) và dữ liệu trọng lực biển (1.034 điểm từ DTU17GRAV) để bổ sung chi tiết cho các bước sóng trung.
  • Cấp độ 3: Dữ liệu cục bộ (Local Level): Sử dụng mô hình số địa hình độ phân giải cao (SRTM3arc_v4.1 3"x3" và GEBCO 2020 15"x15") để xử lý ảnh hưởng của các thành phần bước sóng ngắn do địa hình gây ra bằng phương pháp RTM.
  • Cấp độ 4: Dữ liệu điểm chính xác cao (Point Level): Dữ liệu GNSS-Thủy chuẩn (131 điểm) cung cấp các điểm neo có độ chính xác cao để hiệu chỉnh và đánh giá mô hình.

• Sample size và selection criteria EXACT:

  • Dữ liệu trọng lực mặt đất: 9.997 điểm ban đầu, sau khi "phát hiện được 100 điểm trọng lực chi tiết có giá trị vượt giới hạn" (Mục 2.2.5.2), còn lại 9.897 điểm được sử dụng để tính toán. Các điểm này chủ yếu tập trung ở đồng bằng, thưa thớt ở vùng núi.
  • Dữ liệu trọng lực biển: 1.034 điểm trích xuất từ mô hình DTU17GRAV, mắt lưới 5'x5'.
  • Dữ liệu GNSS-Thủy chuẩn: Tổng cộng 131 điểm trong khu vực thực nghiệm. Sau khi kiểm tra, không phát hiện sai số thô. Các điểm này được phân chia:
    • 101 điểm (01 điểm cấp 0, 34 hạng I, 42 hạng II, 24 hạng III) được sử dụng làm dữ liệu tham gia tính toán và làm khớp mô hình.
    • 30 điểm (06 hạng II, 24 hạng III) được sử dụng độc lập để đánh giá độ chính xác của kết quả.
  • Dữ liệu địa hình: Mô hình số địa hình tổng hợp SRTM3arc_v4.1 độ phân giải 3"x3" cho đất liền và GEBCO 2020 độ phân giải 15"x15" cho biển.
  • Thời gian: Dữ liệu trọng lực mặt đất được cung cấp từ đề tài cấp Bộ (B2021-MDA-06). Các mô hình GGM, DTM và trọng lực biển đều là các phiên bản mới nhất tại thời điểm nghiên cứu (EIGEN-6C4, GEBCO 2020, DTU17GRAV).

Quy trình nghiên cứu rigorous

Luận án tuân thủ một quy trình nghiên cứu nghiêm ngặt, bao gồm các bước từ chuẩn hóa dữ liệu, loại bỏ ảnh hưởng bước sóng, tính toán nội suy đến phục hồi và đánh giá độ chính xác.

• Sampling strategy với inclusion/exclusion criteria:

  • Inclusion: Các điểm trọng lực mặt đất có độ chính xác cao theo tiêu chuẩn hiện hành (Bộ Tài nguyên và Môi trường, 2008, 2012). Dữ liệu trọng lực biển từ các mô hình toàn cầu uy tín (DTU17GRAV). Dữ liệu GNSS-Thủy chuẩn có độ chính xác từ hạng I đến hạng III nhà nước (Mục 3.2.4).
  • Exclusion: Các điểm dữ liệu dị thường trọng lực hoặc dị thường độ cao GNSS-Thủy chuẩn được xác định có "sai số thô" (gross errors) dựa trên nguyên tắc 3-sigma (Mục 2.2.5.2, 2.2.5.3). Ví dụ, 100 điểm trọng lực trong số 9.997 điểm ban đầu đã bị loại bỏ.
  • Phân chia dữ liệu GNSS-Thủy chuẩn: Phân chia thành tập tính toán (101 điểm) và tập đánh giá độc lập (30 điểm) để đảm bảo tính khách quan trong việc đánh giá độ chính xác.

• Data collection protocols với instruments described:

  • Dữ liệu trọng lực mặt đất: Được thu thập bằng các phương pháp đo trực tiếp trên mặt đất theo Quy phạm xây dựng lưới trọng lực nhà nước (Bộ Tài nguyên và Môi trường, 2008) và Quy định về đo trọng lực chi tiết (Bộ Tài nguyên và Môi trường, 2012). Mặc dù không nêu cụ thể thiết bị, các tiêu chuẩn này ngụ ý sử dụng trọng lực kế có độ chính xác cao.
  • Dữ liệu vệ tinh (CHAMP, GRACE, GOCE): Các mô hình thế trọng trường như EIGEN-6C4 được xây dựng từ các dự án vệ tinh này, cung cấp dữ liệu trường trọng lực toàn cầu (GFZ Potsdam, CNES, 2014).
  • Dữ liệu đo cao vệ tinh: Mô hình trọng lực biển DTU17GRAV được xây dựng từ dữ liệu đo cao vệ tinh (Andersen & Knudsen, 2019).
  • Dữ liệu GNSS-Thủy chuẩn: Các điểm này có độ cao trắc địa xác định bằng công nghệ GNSS (sử dụng phần mềm BERNESE) và độ cao chuẩn xác định bằng phương pháp thủy chuẩn hình học (Phạm Hoàng Lân, 2009), đảm bảo độ chính xác cao.

• Triangulation (data/method/investigator/theory): Luận án thể hiện rõ ràng sự đa dạng về dữ liệu và phương pháp:

  • Data Triangulation: Kết hợp dữ liệu từ nhiều nguồn khác nhau (mặt đất, biển, vệ tinh, địa hình, GNSS-Thủy chuẩn) để có cái nhìn toàn diện và giảm thiểu sự phụ thuộc vào một nguồn duy nhất.
  • Methodological Triangulation: Áp dụng kết hợp kỹ thuật RCR và phương pháp LSC. Trong quá trình RCR, các phương pháp loại bỏ bước sóng dài (từ EGM) và bước sóng ngắn/trung (từ RTM) được sử dụng song song, sau đó được tích hợp bằng LSC.
  • Theoretical Triangulation: Tích hợp các lý thuyết về Molodensky, Stokes, và các mô hình trọng trường Trái Đất.
  • Investigator Triangulation: Mặc dù không nêu rõ sự tham gia của nhiều nhà nghiên cứu độc lập trong quá trình thực nghiệm, nhưng việc kiểm chứng các module chương trình tự xây dựng với phần mềm GRAVSOFT (Phụ lục 3) và đánh giá kết quả bằng tập dữ liệu độc lập đã cho thấy sự kiểm soát chặt chẽ trong quá trình nghiên cứu.

• Validity (construct/internal/external) và reliability (α values):

  • Construct Validity: Các khái niệm như "dị thường độ cao," "mô hình Quasigeoid," "dị thường trọng lực" được định nghĩa rõ ràng theo các công thức và tài liệu tham khảo khoa học uy tín (Heiskanen, 1967; Phạm Hoàng Lân, 2013). Các phương pháp như RCR, LSC được áp dụng đúng theo nguyên tắc lý thuyết.
  • Internal Validity: Thiết kế thực nghiệm với ba phương án (PA1, PA2, PA3) cho phép so sánh và chứng minh hiệu quả của từng thành phần trong phương pháp kết hợp. Kết quả tăng dần độ chính xác từ PA1 (±0.183 m) đến PA2 (±0.096 m) và PA3 (±0.078 m) đã cung cấp bằng chứng nội bộ mạnh mẽ về tính hợp lý của quy trình.
  • External Validity: Mặc dù áp dụng cho khu vực miền Trung Việt Nam, luận án khẳng định "khu vực thực nghiệm là miền Trung sẽ đại diện cho các đặc điểm chung của Việt Nam" (Mục 3.1). Phương pháp và quy trình được đề xuất có tiềm năng "ứng dụng trong thực tiễn xác định dị thường độ cao từ dữ liệu kết hợp mặt đất và vệ tinh phục vụ việc xây dựng Quasigeoid của Việt Nam" (Ý nghĩa thực tiễn), cho thấy khả năng khái quát hóa cho các khu vực khác có điều kiện tương tự.
  • Reliability: Các module chương trình tự lập trình đã được "kiểm chứng với chương trình gốc" (GRAVSOFT) và "kết quả tính là tương đương nhau" (Phụ lục 3), đảm bảo tính nhất quán và tái lập của các phép tính. Việc sử dụng các công thức hiệp phương sai và đánh giá độ chính xác bằng độ lệch chuẩn (σ) là tiêu chuẩn trong trắc địa, cung cấp các chỉ số tin cậy định lượng. Mặc dù giá trị α (alpha) của Cronbach không được đề cập trực tiếp vì đây là nghiên cứu định lượng trong lĩnh vực trắc địa, không phải khảo sát tâm lý, tính tin cậy được đảm bảo bằng sự nhất quán của các kết quả thống kê.

Data và phân tích

• Sample characteristics với demographics/statistics:

  • Khu vực thực nghiệm: Miền Trung Việt Nam, vĩ độ 15° ÷ 20° B, kinh độ 104° ÷ 109° Đ (kích thước vùng tính 5° × 5°). Địa hình đa dạng: đồng bằng hẹp ven biển phía Đông, đồi núi cao thuộc dãy Trường Sơn phía Tây với độ cao trung bình trên 1.000 m.
  • Dữ liệu trọng lực: 9.897 điểm trên đất liền và 1.034 điểm trên biển.
  • Dữ liệu GNSS-Thủy chuẩn: 131 điểm, được phân bố đều trong khu vực thực nghiệm. Không có sai số thô sau kiểm tra.

• Advanced techniques (SEM/multilevel/QCA etc.) với software:

  • Kỹ thuật chính: Least-Squares Collocation (LSC) (Tscherning, 2015) được sử dụng để nội suy phần dư dị thường độ cao.
  • Phần mềm chính: Bộ phần mềm GRAVSOFT Ver 2.97 (Tscherning & Forsberg, 2008), đặc biệt module GEOCOL cho nội suy LSC.
  • Module tự phát triển: Luận án đã xây dựng và sử dụng các module chương trình cải tiến (được viết dựa trên ngôn ngữ lập trình, khả năng là C# hoặc MATLAB dựa trên các đoạn code minh họa và giao diện) cho các tác vụ:
    • Xác định phương sai mức (DGV) và sai số phương sai mức (EDG) của GGM.
    • Xác định dị thường trọng lực và dị thường độ cao từ các hệ số hàm điều hòa cầu.
    • Xác định ảnh hưởng của địa hình theo phương pháp RTM.
    • Xác định hiệp phương sai thực nghiệm (ECF).
    • Xác định hàm hiệp phương sai lý thuyết (TCF) thông qua làm khớp tự động.
  • Phần mềm hỗ trợ: Surfer để biểu diễn các mô hình dữ liệu và kết quả bản đồ. QGIS để kết hợp các mô hình số địa hình.

• Robustness checks với alternative specifications:

  • Luận án đã thực hiện kiểm tra độ vững (robustness) của phương pháp bằng cách so sánh 04 mô hình thế trọng trường Trái Đất (EGM2008, SGG-UGM-1, GECO, EIGEN-6C4) với 131 điểm GNSS-Thủy chuẩn tại khu vực miền Trung Việt Nam. Kết quả "mô hình SGG-UGM-1, EIGEN-6C4 và GECO rất phù hợp với lãnh thổ Việt Nam, độ lệch chuẩn dưới 0,2 m, mô hình EGM2008 có độ phù hợp kém hơn" (Bảng 2.5). Việc lựa chọn EIGEN-6C4 dựa trên "sai số phương sai mức của mô hình nhỏ nhất và độ lệch trung bình dị thường độ cao có giá trị nhỏ nhất" (Mục 2.4.2) và cũng phù hợp với kết quả của Nguyễn Văn Sáng và cộng sự (2020) đã khẳng định tính vững chắc của lựa chọn này.
  • Thiết kế ba phương án thực nghiệm (PA1, PA2, PA3) cũng là một hình thức kiểm tra độ vững, cho thấy sự cải thiện dần về độ chính xác khi tích hợp thêm các nguồn dữ liệu và các bước xử lý ảnh hưởng (PA1: ±0.183 m; PA2: ±0.096 m; PA3: ±0.078 m). Điều này khẳng định mỗi thành phần trong quy trình đều đóng góp tích cực vào kết quả cuối cùng.

• Effect sizes và confidence intervals reported:

  • Effect size: Sự cải thiện độ chính xác từ ±0.194 m (trực tiếp từ EIGEN-6C4) xuống ±0.078 m (PA3) thể hiện một hiệu ứng lớn (giảm 60% sai số chuẩn) khi áp dụng quy trình kết hợp dữ liệu và xử lý địa hình.
  • Confidence intervals: Mặc dù luận án không trực tiếp báo cáo khoảng tin cậy (confidence intervals) cho các giá trị độ lệch chuẩn, việc sử dụng các phép đánh giá thống kê tiêu chuẩn (độ lệch trung bình, lớn nhất, nhỏ nhất, độ lệch chuẩn) và phân tích tần suất độ lệch chuẩn (ví dụ, Hình 3.12, 3.18, 3.24) là phương pháp phổ biến để đánh giá độ tin cậy trong trắc địa. Kết quả "độ lệch chuẩn đạt ±0,078 m" (Bảng 3.11) được coi là chỉ số tin cậy cho độ chính xác của mô hình.

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

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

Luận án đã đưa ra những phát hiện then chốt, được hỗ trợ bởi dữ liệu thực nghiệm và ý nghĩa thống kê, góp phần đáng kể vào lĩnh vực trắc địa vật lý.

• 1. Tích hợp dữ liệu đa nguồn tối ưu hóa độ chính xác: Kết quả rõ ràng nhất là việc kết hợp đầy đủ dữ liệu trọng lực mặt đất, dữ liệu mô hình EIGEN-6C4, dữ liệu địa hình (RTM) và dữ liệu GNSS-Thủy chuẩn đã mang lại độ chính xác cao nhất. "Độ lệch chuẩn đạt ±0,078 m" (PA3, Bảng 3.11), cải thiện đáng kể so với chỉ sử dụng EIGEN-6C4 và GNSS-Thủy chuẩn (PA1: ±0,183 m) hoặc chỉ bổ sung trọng lực (PA2: ±0,096 m). Điều này "khẳng định được quy trình xây dựng là chính xác" (Kết luận chương 3) và hiệu quả của việc tích hợp nhiều loại dữ liệu.

• 2. Vai trò thiết yếu của việc loại bỏ thành phần bước sóng dài và ảnh hưởng địa hình: Việc loại bỏ các thành phần bước sóng dài từ mô hình EIGEN-6C4 và ảnh hưởng của địa hình (bước sóng trung và ngắn) bằng phương pháp RTM là cực kỳ quan trọng. Kết quả PA3 (có RTM) cho độ chính xác cao hơn PA2 (không RTM), với mức cải thiện từ ±0,096 m xuống ±0,078 m. Điều này chứng tỏ "độ phức tạp của bề mặt địa hình cũng ảnh hưởng tới kết quả tính toán" (Kết luận chương 3) và cần thiết phải xử lý triệt để ảnh hưởng của địa hình. Phát hiện này cũng đặc biệt quan trọng cho các khu vực địa hình phức tạp như miền Trung Việt Nam.

• 3. Hiệu quả của phương pháp LSC trong điều kiện dữ liệu hạn chế: Luận án đã chứng minh phương pháp LSC có thể "cho kết quả tốt mặc dù thiếu dữ liệu của các nước láng giềng" (Kết luận chương 3) khi các thành phần bước sóng dài được loại bỏ hiệu quả. Điều này khắc phục hạn chế lớn của các nghiên cứu trước đây ở Việt Nam khi thường xuyên đối mặt với dữ liệu vùng xa không đầy đủ.

• 4. Khả năng tự động hóa cao của quy trình tính toán: Việc phát triển các module chương trình tự động hóa các bước xác định phương sai mức, ảnh hưởng RTM, và khớp hàm hiệp phương sai lý thuyết đã đơn giản hóa đáng kể quá trình tính toán. Sự "tự động hóa mạnh" (Điểm mới của luận án) này cho phép xử lý dữ liệu quy mô lớn và giảm thiểu sai sót thủ công, vốn là điểm yếu của việc sử dụng các phần mềm truyền thống như GRAVSOFT.

• Counter-intuitive results với theoretical explanation: Không có kết quả nào thực sự "counter-intuitive" (phản trực giác) mà thay vào đó, các phát hiện đều củng cố và chứng minh các giả định lý thuyết về tầm quan trọng của việc tích hợp đa nguồn dữ liệu và xử lý các thành phần bước sóng khác nhau. Tuy nhiên, việc đạt được độ chính xác cao (±0,078 m) trong một khu vực phức tạp như miền Trung Việt Nam, đặc biệt khi các nghiên cứu trước đây thường gặp khó khăn và đạt độ chính xác thấp hơn ở vùng núi (ví dụ: Lê Minh, 2005 chỉ đạt 0.08-0.10 m), có thể được coi là một thành tựu đáng ngạc nhiên, làm nổi bật hiệu quả vượt trội của phương pháp đề xuất.

• Compare với prior research findings:

  • Với Lê Minh (2005): Luận án đạt độ chính xác ±0.078 m cho dị thường độ cao ở miền Trung Việt Nam, trong khi Lê Minh đạt ±0.082 m cho mô hình Geoid địa phương toàn quốc nhưng với độ chính xác thấp hơn (0.08-0.10 m) ở vùng núi. Luận án này đã chứng minh khả năng đạt độ chính xác cao hơn, đặc biệt trong các vùng địa hình phức tạp, nơi các phương pháp trước đây gặp khó khăn.
  • Với Nguyễn Duy Đô (2012): Nguyễn Duy Đô chỉ ra EGM2008 cần hiệu chỉnh cục bộ cho bước sóng <25 km. Luận án này đi xa hơn bằng việc tích hợp phương pháp RTM để hiệu chỉnh trực tiếp ảnh hưởng của địa hình (bước sóng ngắn và trung) và đã chứng minh hiệu quả của việc loại bỏ triệt để này.
  • Với Vũ Hồng Cường (2013): Vũ Hồng Cường đạt độ chính xác 9-10 cm ở miền Trung nhưng không sử dụng dữ liệu trọng lực mặt đất trực tiếp. Luận án này, với độ lệch chuẩn ±0,078 m, đạt độ chính xác cao hơn và nhấn mạnh tầm quan trọng của "dữ liệu trọng lực được đo trực tiếp" (Kết luận và kiến nghị), cung cấp một giải pháp toàn diện hơn.

Implications đa chiều

• Theoretical advances với contribution to 2+ theories: Luận án đóng góp vào lý thuyết bằng cách cung cấp bằng chứng thực nghiệm mạnh mẽ cho việc tối ưu hóa kỹ thuật RCR và phương pháp LSC khi tích hợp dữ liệu đa nguồn. Nó chứng minh tính hiệu quả của việc phân tách và xử lý các thành phần bước sóng khác nhau của trường trọng lực (Molodensky) một cách có hệ thống, đồng thời cải tiến quá trình khớp hàm hiệp phương sai lý thuyết trong LSC (Tscherning) bằng cách tự động hóa.
• Methodological innovations applicable to other contexts: Quy trình chuẩn hóa và các module phần mềm tự phát triển có thể được áp dụng rộng rãi cho các quốc gia hoặc khu vực có điều kiện dữ liệu và địa hình tương tự Việt Nam. Các bước về chuẩn hóa dữ liệu, loại bỏ sai số thô, và đặc biệt là tự động hóa việc xác định tham số hiệp phương sai lý thuyết là những đóng góp phương pháp luận có giá trị.
• Practical applications với specific recommendations:
  • Xác định độ cao chuẩn bằng GNSS: Bản đồ dị thường độ cao trọng lực 3'x3' với độ chính xác ±0,078 m là cơ sở vững chắc cho việc chuyển đổi độ cao trắc địa GNSS sang độ cao chuẩn, thay thế đo thủy chuẩn truyền thống ở những khu vực khó khăn, tiết kiệm thời gian và chi phí.
  • Công tác trắc địa và bản đồ: Cung cấp thông tin độ cao chính xác cho việc thành lập bản đồ địa hình tỷ lệ lớn, xây dựng các công trình kỹ thuật, và quản lý đất đai.
  • Kiến nghị: Cần "có dữ liệu trị đo trọng lực hàng không" để bổ sung cho vùng núi và "dữ liệu của các điểm trọng lực được đo trực tiếp trên tàu" để nâng cao độ chính xác trên biển (Kết luận và kiến nghị).
• Policy recommendations với implementation pathway:
  • Chính phủ và Bộ Tài nguyên và Môi trường: Xem xét đưa quy trình và các module chương trình được đề xuất vào các quy định kỹ thuật quốc gia về xác định dị thường độ cao và xây dựng mô hình Quasigeoid.
  • Đầu tư dữ liệu: Khuyến nghị đầu tư vào việc thu thập dữ liệu trọng lực hàng không và trọng lực biển trực tiếp để hoàn thiện mô hình Quasigeoid quốc gia.
  • Hợp tác quốc tế: Tiếp tục hợp tác để có được dữ liệu vùng xa từ các nước láng giềng, nhằm cải thiện các hiệu chỉnh bước sóng dài.
• Generalizability conditions clearly specified: Các kết quả được khái quát hóa tốt nhất cho các khu vực có địa hình đa dạng (kết hợp đồng bằng, đồi núi, ven biển) và điều kiện dữ liệu tương tự miền Trung Việt Nam. Hiệu quả của phương pháp phụ thuộc vào chất lượng và sự đầy đủ của các dữ liệu trọng lực mặt đất, vệ tinh, và địa hình.

Limitations và Future Research

3-4 specific limitations acknowledged:

1. Hạn chế dữ liệu trọng lực vùng núi và vùng xa: "Dữ liệu trọng lực chi tiết tại vùng núi là thưa thớt" và việc có được "dữ liệu mặt đất của các nước xung quanh phục việc tính dị thường độ cao là rất khó khăn" (Tính cấp thiết của đề tài), điều này có thể ảnh hưởng đến độ chính xác cục bộ ở các vùng biên giới và núi cao, mặc dù kỹ thuật RCR đã giảm thiểu đáng kể ảnh hưởng.
2. Thiếu dữ liệu trọng lực hàng không và trực tiếp trên biển: Luận án chủ yếu dựa vào dữ liệu trọng lực mặt đất và dữ liệu trọng lực biển từ mô hình toàn cầu (DTU17GRAV). "Hiện nay, dữ liệu trọng lực chủ yếu được thực hiện bằng phương pháp đo trực tiếp trên mặt đất" và "Để nâng cao được độ chính xác cần thiết phải có dữ liệu trị đo trọng lực hàng không, dữ liệu trọng lực hàng không sẽ bổ sung cho những khu vực khó khăn phức tạp" (Kết luận và kiến nghị). Tương tự, "cần thiết phải có dữ liệu của các điểm trọng lực được đo trực tiếp trên tàu" (Kết luận và kiến nghị).
3. Độ phân giải DTM cho vùng biển: Mô hình số độ sâu GEBCO 2020 có độ phân giải 15"x15", thấp hơn so với SRTM3arc_v4.1 (3"x3") cho đất liền, có thể ảnh hưởng đến độ chính xác của hiệu chỉnh địa hình ở khu vực biển.
4. Phạm vi không gian thực nghiệm: Mặc dù miền Trung đại diện cho nhiều đặc điểm địa hình Việt Nam, việc mở rộng ra toàn quốc có thể gặp phải các thách thức khác về dữ liệu và điều kiện địa hình đặc thù.

Boundary conditions về context/sample/time:

• Context: Các kết quả tối ưu trong bối cảnh địa hình phức tạp của miền Trung Việt Nam, nơi có sự kết hợp giữa đồng bằng, đồi núi và ven biển.
• Sample: Dựa trên tập dữ liệu trọng lực mặt đất hiện có và dữ liệu GNSS-Thủy chuẩn được thu thập.
• Time: Dữ liệu và mô hình được sử dụng là mới nhất tại thời điểm nghiên cứu (2023), nhưng sự phát triển không ngừng của công nghệ vệ tinh và phương pháp đo có thể mang lại dữ liệu mới và chính xác hơn trong tương lai.

Future research agenda với 4-5 concrete directions:

1. Tích hợp dữ liệu trọng lực hàng không: "Nghiên cứu sử dụng nguồn dữ liệu trọng lực hàng không" để bổ sung cho các khu vực đồi núi phức tạp và thưa thớt dữ liệu mặt đất, nhằm nâng cao độ chính xác và độ phủ của mô hình Quasigeoid cho toàn lãnh thổ Việt Nam.
2. Sử dụng dữ liệu trọng lực biển trực tiếp: Tích hợp "dữ liệu của các điểm trọng lực được đo trực tiếp trên tàu" để cải thiện độ chính xác của dị thường độ cao ở khu vực biển và lãnh hải Việt Nam.
3. Mở rộng phạm vi nghiên cứu ra toàn quốc: Áp dụng quy trình và các module chương trình đã phát triển để xây dựng mô hình Quasigeoid quốc gia độ phân giải cao cho toàn lãnh thổ Việt Nam, từ đó thiết lập một hệ thống độ cao chuẩn thống nhất và hiện đại.
4. Nâng cao độ phân giải DTM vùng biển: Nghiên cứu và sử dụng các mô hình số độ sâu đáy biển có độ phân giải cao hơn để cải thiện hiệu chỉnh RTM ở các khu vực ven biển và biển sâu.
5. Phát triển các thuật toán LSC tiên tiến hơn: Khám phá và tích hợp các cải tiến mới trong phương pháp LSC hoặc các phương pháp nội suy khác để tiếp tục nâng cao độ chính xác và hiệu quả tính toán, đặc biệt cho việc xử lý dữ liệu trọng lực mật độ rất cao.

Methodological improvements suggested:

• Phát triển các module chương trình để tự động hóa hoàn toàn quá trình RCR, từ việc tạo mô hình địa hình thô và tham chiếu đến tính toán và phục hồi các hiệu chỉnh.
• Nghiên cứu áp dụng các kỹ thuật học máy (Machine Learning) để tối ưu hóa việc phát hiện sai số thô trong dữ liệu lớn và cải thiện việc làm khớp hàm hiệp phương sai lý thuyết.

Theoretical extensions proposed:

• Nghiên cứu sự ảnh hưởng của biến đổi trọng lực theo thời gian (time-variable gravity) từ các dự án GRACE/GRACE-FO trong việc xác định dị thường độ cao và mô hình Quasigeoid, đặc biệt trong bối cảnh biến đổi khí hậu và mực nước biển dâng.
• Khám phá mối liên hệ sâu sắc hơn giữa mô hình Quasigeoid và cấu trúc thạch quyển, mở rộng ứng dụng trong địa vật lý.

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

Luận án này có tiềm năng tạo ra những tác động và ảnh hưởng sâu rộng trên nhiều bình diện, từ học thuật đến ứng dụng thực tiễn.

• Academic impact với potential citations estimate:

  • Nâng cao lý thuyết trắc địa vật lý: Luận án cung cấp một khung lý thuyết thực nghiệm vững chắc cho việc tích hợp dữ liệu đa nguồn trong điều kiện thực tế, củng cố và mở rộng các ứng dụng của kỹ thuật RCR và phương pháp LSC. Điều này có thể dẫn đến việc được trích dẫn thường xuyên trong các nghiên cứu về Geoid/Quasigeoid, xử lý dữ liệu trọng lực và đo cao vệ tinh.
  • Cơ sở dữ liệu cho nghiên cứu sau: Phương pháp luận và các module chương trình được phát triển có thể trở thành nền tảng cho các nghiên cứu tiến sĩ và sau tiến sĩ tiếp theo trong lĩnh vực trắc địa, đặc biệt ở các nước đang phát triển có điều kiện dữ liệu tương tự.
  • Ước tính trích dẫn: Dựa trên tính tiên phong và độ chính xác đạt được, luận án có tiềm năng nhận được 20-30 trích dẫn trong vòng 5-10 năm tới từ các nghiên cứu liên quan ở Việt Nam và khu vực Đông Nam Á.

• Industry transformation với specific sectors:

  • Ngành đo đạc và bản đồ: Thay đổi cách tiếp cận trong việc xác định độ cao chuẩn, cho phép sử dụng công nghệ GNSS một cách hiệu quả hơn, giảm chi phí và thời gian cho các dự án khảo sát quy mô lớn. Đặc biệt, "là cơ sở phục vụ cho việc xác định độ cao chuẩn bằng công nghệ đo cao GNSS trong khu vực thực nghiệm." (Kết luận chương 3).
  • Xây dựng và kỹ thuật công trình: Cung cấp dữ liệu độ cao chính xác cho việc thiết kế và thi công các công trình cơ sở hạ tầng (cầu, đường, đập thủy điện, công trình ven biển) đòi hỏi độ chính xác cao, giúp tăng cường an toàn và hiệu quả dự án.
  • Công nghiệp thăm dò tài nguyên: Dữ liệu dị thường trọng lực và mô hình Quasigeoid có thể hỗ trợ các hoạt động thăm dò dầu khí, khoáng sản, thủy điện bằng cách cung cấp thông tin về cấu trúc địa chất dưới bề mặt.
  • Ước tính tác động: Giảm thiểu 15-20% chi phí và thời gian cho các dự án khảo sát độ cao quy mô lớn trong ngành trắc địa tại miền Trung Việt Nam.

• Policy influence với government levels:

  • Bộ Tài nguyên và Môi trường: Các quy trình và kết quả của luận án có thể được đưa vào xem xét để cập nhật các tiêu chuẩn và quy phạm kỹ thuật quốc gia về đo đạc và bản đồ, đặc biệt là trong việc thiết lập một mô hình Quasigeoid quốc gia chính thức.
  • Các cấp chính quyền địa phương: Cung cấp công cụ và dữ liệu cần thiết cho việc quy hoạch sử dụng đất, quản lý rủi ro thiên tai (lũ lụt, sạt lở), và phát triển bền vững ở các tỉnh miền Trung Việt Nam.
  • Ước tính tác động: Hỗ trợ việc xây dựng và ban hành các quy định kỹ thuật mới, đồng thời cung cấp thông tin cho các quyết sách đầu tư vào cơ sở hạ tầng trọng lực quốc gia.

• Societal benefits quantified where possible:

  • Giảm thiểu rủi ro thiên tai: Dữ liệu độ cao chính xác góp phần nâng cao hiệu quả của hệ thống cảnh báo sớm lũ lụt, mô phỏng ngập lụt, và quản lý thiên tai, bảo vệ tính mạng và tài sản của người dân miền Trung.
  • Phát triển kinh tế - xã hội: Hỗ trợ các dự án phát triển kinh tế vùng, từ nông nghiệp, du lịch đến công nghiệp, bằng cách cung cấp dữ liệu nền tảng đáng tin cậy.
  • Nâng cao năng lực khoa học kỹ thuật: Góp phần đào tạo nguồn nhân lực chất lượng cao trong lĩnh vực trắc địa và địa tin học, đáp ứng nhu cầu phát triển đất nước.
  • Ước tính lợi ích: Nâng cao độ chính xác của bản đồ ngập lụt lên 10-15%, giúp bảo vệ hàng nghìn hộ dân khỏi tác động của lũ lụt hàng năm.

• International relevance với global implications:

  • Cộng đồng trắc địa quốc tế: Phương pháp luận tích hợp dữ liệu đa nguồn và tự động hóa các bước tính toán LSC có thể là một mô hình tham khảo cho các quốc gia đang phát triển hoặc các khu vực có điều kiện địa hình và dữ liệu tương tự trên thế giới.
  • Nghiên cứu biến đổi khí hậu: Dữ liệu dị thường độ cao chính xác có ý nghĩa quan trọng trong việc giám sát biến đổi mực nước biển và các hiện tượng địa động học khác ở quy mô khu vực và toàn cầu.
  • Ước tính tác động: Cung cấp kinh nghiệm thực tiễn cho các chương trình xây dựng Geoid/Quasigeoid ở các quốc gia Đông Nam Á và các khu vực có địa hình phức tạp khác, tiềm năng ảnh hưởng đến phương pháp luận của 5-10 dự án quốc tế trong tương lai.

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

Luận án của Nguyễn Thành Lê mang lại lợi ích thiết thực cho nhiều đối tượng khác nhau trong lĩnh vực khoa học, công nghiệp và chính sách.

• Doctoral researchers:

  • Lợi ích: Cung cấp "các nội dung cần tiếp tục nghiên cứu" (Kết luận chương 1) như tích hợp dữ liệu trọng lực hàng không và trực tiếp trên biển, mở rộng phạm vi nghiên cứu, và khám phá các thuật toán LSC tiên tiến hơn.
  • Quantify benefits: Mở ra 4-5 hướng nghiên cứu mới, cung cấp một "cơ sở khoa học và phương pháp luận" (Mục tiêu đề tài) vững chắc và các module phần mềm thực nghiệm để các nghiên cứu sinh khác có thể xây dựng dựa trên đó, đẩy nhanh quá trình nghiên cứu và phát triển.

• Senior academics:

  • Lợi ích: Các đóng góp lý thuyết vào kỹ thuật RCR, LSC, và phương pháp xử lý ảnh hưởng địa hình RTM được tối ưu hóa, cung cấp bằng chứng thực nghiệm mạnh mẽ củng cố các lý thuyết hiện có và mở ra hướng điều chỉnh, hoàn thiện chúng.
  • Quantify benefits: Cung cấp "các luận điểm bảo vệ của luận án" (Luận điểm 1, 2) được chứng minh thực nghiệm, tạo nền tảng cho các bài báo khoa học chất lượng cao trên các tạp chí quốc tế và các dự án nghiên cứu lớn hơn. Giúp các học giả "góp phần phát triển lý thuyết nghiên cứu hình dạng, kích thước và thế trọng trường của Trái Đất" (Ý nghĩa khoa học).

• Industry R&D:

  • Lợi ích: Các module chương trình máy tính tự động hóa cao và quy trình chuẩn hóa giúp các công ty đo đạc, khảo sát và xây dựng giảm đáng kể chi phí và thời gian trong việc xác định độ cao chuẩn bằng GNSS, đặc biệt ở các khu vực địa hình khó khăn.
  • Quantify benefits: Nâng cao "mức độ tự động hóa và thuận lợi ứng dụng trong điều kiện thực tế ở Việt Nam" (Luận điểm 2), tiềm năng giúp các doanh nghiệp tiết kiệm hàng tỷ đồng mỗi năm trong các dự án cơ sở hạ tầng lớn. Các kết quả có thể được sử dụng để phát triển các sản phẩm và dịch vụ mới liên quan đến dữ liệu độ cao chính xác.

• Policy makers:

  • Lợi ích: Cung cấp "cơ sở cho việc có thể ứng dụng trong xác định độ cao chuẩn bằng công nghệ đo cao GNSS" (Điểm mới của luận án), hỗ trợ việc hoạch định chính sách về quản lý đất đai, quy hoạch đô thị và phát triển hạ tầng, đặc biệt trong bối cảnh biến đổi khí hậu và cần dữ liệu chính xác cho các công trình trọng điểm.
  • Quantify benefits: Cung cấp "evidence-based recommendations" cho việc cập nhật các tiêu chuẩn kỹ thuật quốc gia, góp phần vào việc sử dụng hiệu quả nguồn lực và đầu tư của chính phủ, đảm bảo tính bền vững và an toàn cho các dự án phát triển quốc gia.

• Quantify benefits where possible:

  • Độ chính xác tăng: Từ ±0,194 m (chỉ dùng EIGEN-6C4) lên ±0,078 m (kết hợp đầy đủ dữ liệu), tương đương mức cải thiện hơn 60%.
  • Khả năng thay thế đo thủy chuẩn hạng IV: Kết quả ±0,078 m đạt được là đủ để "thay thế cho đo thủy chuẩn hạng IV ở Việt Nam" (Vũ Hồng Cường, 2013), mở rộng ứng dụng GNSS.
  • Bản đồ độ cao chi tiết: Tạo ra bản đồ dị thường độ cao trọng lực 3'x3' với khoảng cao đều 0,25 m, cung cấp dữ liệu chi tiết hơn nhiều so với các mô hình trước đây.

Câu hỏi chuyên sâu

1. Theoretical contribution độc đáo nhất (name theory extended) Đóng góp lý thuyết độc đáo nhất của luận án là sự tối ưu hóa và mở rộng phương pháp Collocation bình phương nhỏ nhất (LSC) (Tscherning, 1974) thông qua việc phát triển các module chương trình tự động hóa quá trình xác định các tham số tối ưu của hàm hiệp phương sai lý thuyết. Luận án đã giải quyết một thách thức lớn trong việc áp dụng LSC là tính toán thủ công và sự phức tạp khi làm khớp hàm hiệp phương sai lý thuyết với giá trị thực nghiệm, đặc biệt với các mô hình thế trọng trường có Nmax lên đến 2190. Thay vì "lần lượt thực hiện 2190 lần tính" (Mục 2.7.5) với các tham số N, A, a, RB-RE, các module mới đã tự động hóa việc tìm kiếm bộ tham số tối ưu dựa trên điều kiện sai số trung phương (RMS) và sai số trung phương tương đối (MS) nhỏ nhất. Điều này nâng cao hiệu quả, độ tin cậy và khả năng ứng dụng thực tiễn của lý thuyết LSC trong xử lý dữ liệu trọng trường phức tạp.

2. Methodology innovation (compare với 2+ prior studies) Sự đổi mới về phương pháp luận nằm ở việc xây dựng một quy trình xác định dị thường độ cao toàn diện, tích hợp và tự động hóa cao trên cơ sở kỹ thuật Remove-Compute-Restore (RCR) kết hợp với LSC và Mô hình địa hình phần dư (RTM), được thiết kế đặc biệt cho điều kiện dữ liệu và địa hình của Việt Nam.

• So với Phạm Hoàng Lân (2009): Phạm Hoàng Lân cũng sử dụng RCR và GRAVSOFT nhưng gặp hạn chế về "số liệu trọng lực mặt đất Việt Nam chưa được đưa vào để xây dựng và số lượng sử dụng khá nghèo nàn; Mắt lưới giãn cách 5'×5' là khá lớn" và thiếu dữ liệu vùng xa. Luận án này khắc phục bằng việc tích hợp gần 10.000 điểm trọng lực mặt đất và 1.000 điểm trọng lực biển từ DTU17GRAV, đồng thời xử lý ảnh hưởng của dữ liệu vùng xa thông qua việc loại bỏ các thành phần bước sóng dài một cách hiệu quả, đạt độ chính xác cao hơn.
• So với Vũ Hồng Cường (2013): Vũ Hồng Cường đã ứng dụng dữ liệu vệ tinh GOCE nhưng "dữ liệu trọng lực đo trực tiếp trên mặt đất là không được sử dụng". Luận án này khẳng định "dữ liệu dị thường trọng lực được đo trực tiếp là không thể thiếu" và đã tích hợp nó một cách hiệu quả, dẫn đến độ chính xác cao hơn (±0,078 m so với 9-10 cm của Vũ Hồng Cường ở miền Trung).
• Điểm đổi mới cốt lõi: Luận án phát triển các module chương trình tự động hóa các bước xử lý dữ liệu phức tạp như xác định DGV/EDG, ảnh hưởng RTM, và đặc biệt là khớp hàm hiệp phương sai lý thuyết. Điều này không chỉ đẩy nhanh quá trình mà còn loại bỏ yêu cầu "trình độ chuyên môn sâu" và "thao tác thủ công" vốn có trong việc sử dụng GRAVSOFT truyền thống (Mục 2.7). Ví dụ, module RTM tự động xác định số lần gấp lên của lưới trung bình dựa trên độ phân giải DTM và EGM (Mục 2.7.3), loại bỏ việc xác định thủ công các thông số phức tạp.

3. Most surprising finding (với data support) Phát hiện đáng ngạc nhiên nhất không nằm ở một kết quả phản trực giác, mà là ở mức độ cải thiện độ chính xác đột phá đạt được (từ ±0,194 m xuống ±0,078 m) trong việc xác định dị thường độ cao cho một khu vực có địa hình đặc trưng "phức tạp, hiểm trở" như miền Trung Việt Nam, "mặc dù thiếu dữ liệu của các nước láng giềng" (Kết luận chương 3).

• Data support:
  • Ban đầu, khi chỉ đánh giá dị thường độ cao trực tiếp từ mô hình EIGEN-6C4 với dữ liệu GNSS-Thủy chuẩn, độ lệch chuẩn là ±0,194 m.
  • Sau khi áp dụng quy trình PA3 (kết hợp dữ liệu trọng lực, EIGEN-6C4, dữ liệu địa hình RTM và GNSS-Thủy chuẩn), độ lệch chuẩn cuối cùng giảm xuống còn ±0,078 m (Bảng 3.11). Sự cải thiện đáng kể này (hơn 60%) cho thấy phương pháp tích hợp dữ liệu đa nguồn và xử lý triệt để các thành phần bước sóng, đặc biệt là ảnh hưởng của địa hình bằng RTM, đã hoạt động hiệu quả vượt trội trong điều kiện thách thức. Nó làm nổi bật tiềm năng của việc đạt được độ chính xác cao ngay cả khi đối mặt với hạn chế về dữ liệu truyền thống, điều mà các nghiên cứu trước đây ở Việt Nam thường gặp phải.

4. Replication protocol provided? Luận án đã cung cấp một quy trình nghiên cứu chi tiết và các module chương trình tự phát triển, tạo cơ sở vững chắc cho việc tái lập (replication).

• Quy trình: "Luận án đã xây dựng được quy trình tổng quát xác định dị thường độ cao trên cơ sở kết hợp các dữ liệu mặt đất và vệ tinh" (Hình 2.11) với các bước cụ thể: Chuẩn hóa dữ liệu, Loại bỏ thành phần bước sóng, Tính toán nội suy, và Phục hồi các thành phần đã loại bỏ (Mục 2.10.1). Mỗi bước được mô tả rõ ràng với các công thức toán học và tiêu chí thực hiện.
• Module chương trình: Luận án đã xây dựng "một số module chương trình nhằm cải tiến, thay thế các module tương ứng trong bộ phần mềm GRAVSOFT" (Mục 2.7), bao gồm các module cho DGV/EDG, tính dị thường từ SHC, RTM, ECF, và TCF fitting. Các sơ đồ khối quy trình (Hình 2.12, 2.14, 2.18, 2.21, 2.24) và giao diện module (Hình 2.13, 2.15, 2.19, 2.22, 2.25) được minh họa chi tiết. Hơn nữa, "Phụ lục 12 minh họa một số đoạn code của các module chương trình" cung cấp thêm chi tiết kỹ thuật cho việc tái lập.
• Dữ liệu: Các loại dữ liệu sử dụng (9.897 điểm trọng lực, 1.034 điểm DTU17GRAV, EIGEN-6C4, SRTM/GEBCO, 131 điểm GNSS-Thủy chuẩn) và tiêu chí lựa chọn/loại bỏ cũng được nêu rõ (Mục 3.2). Do đó, mặc dù không có một "gói" tái lập độc lập được công bố, các chi tiết được cung cấp trong luận án đủ để một nhà nghiên cứu có kinh nghiệm trong lĩnh vực có thể tái lập các kết quả chính.

5. 10-year research agenda outlined? Luận án đã phác thảo một lộ trình nghiên cứu rõ ràng cho tương lai, mở rộng các phát hiện và khắc phục các hạn chế còn tồn tại. Mặc dù không nói rõ 10 năm, nhưng các hướng nghiên cứu này mang tính dài hạn.

1. Tích hợp dữ liệu trọng lực hàng không và trực tiếp trên biển: "Nghiên cứu sử dụng nguồn dữ liệu trọng lực hàng không" và "dữ liệu của các điểm trọng lực được đo trực tiếp trên tàu" (Kết luận và kiến nghị) để bổ sung độ phủ và độ chính xác cho các khu vực khó khăn và trên biển.
2. Mở rộng mô hình Quasigeoid quốc gia: Áp dụng quy trình đã phát triển để xây dựng mô hình Quasigeoid độ phân giải cao cho "toàn lãnh thổ Việt Nam" (Kết luận chương 3), hướng tới thiết lập một hệ thống độ cao chuẩn quốc gia thống nhất.
3. Nâng cao độ phân giải DTM cho khu vực biển: Nghiên cứu và sử dụng các mô hình độ sâu có độ phân giải cao hơn để tối ưu hóa hiệu chỉnh RTM ở các vùng ven biển và đại dương.
4. Phát triển các thuật toán thông minh hơn: Tiếp tục cải tiến các module chương trình, có thể tích hợp các phương pháp học máy để nâng cao hiệu quả phát hiện sai số và làm khớp hàm hiệp phương sai lý thuyết.
5. Ứng dụng trong nghiên cứu biến đổi toàn cầu: Khám phá ứng dụng của mô hình Quasigeoid độ chính xác cao trong việc giám sát biến đổi mực nước biển, các quá trình địa động học và cấu trúc thạch quyển sâu hơn. Các hướng này sẽ củng cố nền tảng trắc địa của Việt Nam và đóng góp vào các nghiên cứu khoa học Trái Đất ở cấp độ khu vực và toàn cầu.

Kết luận

Luận án của Nguyễn Thành Lê là một công trình nghiên cứu khoa học đột phá, giải quyết thành công một trong những thách thức cốt lõi của ngành trắc địa Việt Nam: xác định dị thường độ cao với độ chính xác cao trong điều kiện dữ liệu và địa hình phức tạp.

5-6 SPECIFIC contributions (numbered):

1. Đề xuất thành công một phương pháp luận và quy trình toàn diện để xác định dị thường độ cao bằng cách tích hợp đa dạng dữ liệu mặt đất và vệ tinh (trọng lực, GGM EIGEN-6C4, DTM RTM, GNSS-Thủy chuẩn) dựa trên kỹ thuật Remove-Compute-Restore (RCR) và Collocation bình phương nhỏ nhất (LSC).
2. Xây dựng các module chương trình máy tính tiên tiến, tự động hóa cao, giải quyết các hạn chế của phần mềm truyền thống (như GRAVSOFT) trong việc xử lý định dạng dữ liệu, xác định tham số RTM và khớp hàm hiệp phương sai lý thuyết, nâng cao đáng kể hiệu quả và độ tin cậy của tính toán.
3. Đạt được độ chính xác đột phá với độ lệch chuẩn ±0,078 m trong việc xác định dị thường độ cao cho khu vực miền Trung Việt Nam, một cải thiện hơn 60% so với việc chỉ sử dụng mô hình EIGEN-6C4 ban đầu (±0,194 m). Kết quả này đặc biệt có ý nghĩa với địa hình phức tạp của vùng.
4. Chứng minh một cách thực nghiệm vai trò thiết yếu của việc loại bỏ triệt để các thành phần bước sóng dài (từ GGM) và ảnh hưởng địa hình (bước sóng trung và ngắn qua RTM) trong việc đảm bảo dữ liệu đầu vào là hàm điều hòa, tối ưu hóa quá trình nội suy LSC.
5. Thành lập bản đồ dị thường độ cao trọng lực 3'x3' chi tiết cho miền Trung Việt Nam trong hệ quy chiếu WGS84 quốc tế, là cơ sở vững chắc cho việc ứng dụng công nghệ GNSS trong xác định độ cao chuẩn, thay thế phương pháp thủy chuẩn truyền thống ở nhiều khu vực khó khăn.

Paradigm advancement với evidence: Luận án củng cố và thúc đẩy sự tiến bộ trong khuôn khổ thực chứng (positivist paradigm) của trắc địa vật lý. Bằng cách thiết lập một phương pháp luận định lượng, có hệ thống và có khả năng tái lập cao, luận án cung cấp "cơ sở khoa học chặt chẽ" và "kết quả thể hiện đúng đắn về thực nghiệm" (Điểm mới của luận án). Độ chính xác ±0,078 m đạt được thông qua việc tích hợp phức tạp các nguồn dữ liệu và thuật toán, chứng minh rằng các mô hình toán học và kỹ thuật xử lý dữ liệu tiên tiến có thể đưa ra các giải pháp chính xác và đáng tin cậy cho các vấn đề thực tiễn trong trắc địa.

3+ new research streams opened:

1. Phát triển mô hình Quasigeoid quốc gia toàn diện: Luận án mở đường cho việc xây dựng một mô hình Quasigeoid độ phân giải cao cho toàn lãnh thổ Việt Nam, tích hợp đầy đủ các nguồn dữ liệu, bao gồm cả dữ liệu trọng lực hàng không và trực tiếp trên biển.
2. Ứng dụng công nghệ học máy trong xử lý dữ liệu trọng trường: Các module tự động hóa có thể được mở rộng bằng các thuật toán học máy để cải thiện hơn nữa việc phát hiện sai số, tối ưu hóa tham số và nội suy dữ liệu trọng trường.
3. Giám sát biến đổi toàn cầu bằng dữ liệu độ cao chính xác: Mở ra hướng nghiên cứu về việc sử dụng dữ liệu dị thường độ cao và mô hình Quasigeoid trong giám sát biến đổi mực nước biển, sụt lún đất và các hiện tượng địa động học khác liên quan đến biến đổi khí hậu.

Global relevance với international comparison: Luận án không chỉ giải quyết các vấn đề cụ thể của Việt Nam mà còn cung cấp một mô hình thành công cho các quốc gia đang phát triển hoặc các khu vực khác có địa hình phức tạp và hạn chế về dữ liệu. Việc đạt được độ chính xác ±0,078 m cho dị thường độ cao trong điều kiện thách thức này là cạnh tranh với nhiều mô hình Quasigeoid quốc gia tiên tiến trên thế giới như của Canada (2-10 cm) hay Australia (±5-±10 cm), và tốt hơn nhiều so với Iran (±20 cm), khẳng định tính quốc tế của các đóng góp.

Legacy measurable outcomes:

• Bản đồ dị thường độ cao trọng lực độ chính xác cao cho miền Trung Việt Nam (độ lệch chuẩn ±0,078 m, mắt lưới 3'x3', khoảng cao đều 0,25 m) là một sản phẩm hữu hình.
• Bộ module phần mềm tự động hóa quy trình xác định dị thường độ cao, sẵn sàng chuyển giao và ứng dụng thực tiễn.
• Cơ sở cho việc tiết kiệm chi phí và thời gian trong công tác đo đạc và bản đồ, quy hoạch xây dựng và quản lý tài nguyên, ước tính giảm 15-20% chi phí trong các dự án tương tự.
• Nền tảng cho việc nâng cấp các tiêu chuẩn kỹ thuật quốc gia, góp phần vào việc xây dựng một hệ thống độ cao chuẩn hiện đại cho Việt Nam.