Tổng quan nghiên cứu

Trắc địa vật lý đóng vai trò nền tảng trong việc xác định hình dạng Trái Đất và mô hình hóa trường trọng trường thông qua các đại lượng dị thường độ cao, độ lệch dây dọi và dị thường trọng lực. Trong thực tế đo đạc, giá trị trọng lực thực quan sát trên bề mặt địa hình gồ ghề luôn có sự sai khác lớn so với trọng lực chuẩn trên mặt Ellipsoid quy chiếu do ảnh hưởng của lực ly tâm Trái Đất với vận tốc góc 7292115 x 10^-11 rad/s, hằng số trọng trường GM = 3986004.418 x 10^8 m3/s2 và lực hấp dẫn của các khối đất đá bên ngoài mặt Geoid. Để giải quyết bài toán biên Stokes nhằm xác định bề mặt Geoid với độ chính xác đạt mức 1 decimet (10 cm) hoặc cao hơn, toàn bộ khối lượng địa hình phía trên mặt Geoid phải được xử lý và chuyển dịch về mặt thủy chuẩn một cách khoa học.

Vấn đề đặt ra là việc tính toán thủ công các số hiệu chỉnh trọng lực cho mạng lưới điểm đo lớn đòi hỏi khối lượng thời gian khổng lồ và tiềm ẩn nhiều sai số tích lũy. Mục tiêu trọng tâm của nghiên cứu là xây dựng quy trình thuật toán chuẩn hóa và phát triển chương trình tính toán tự động hóa hoàn toàn trên ngôn ngữ Matlab, xử lý đồng bộ 5 loại hiệu chỉnh cốt lõi: hiệu chỉnh khoảng không (Free-air reduction), hiệu chỉnh Bouguer (Bouguer reduction), hiệu chỉnh địa hình (Terrain correction), hiệu chỉnh đẳng tĩnh địa hình (Topographic-Isostatic reduction) theo mô hình Airy - Heiskanen và hiệu chỉnh ảnh hưởng gián tiếp (The indirect effect). Nghiên cứu được triển khai thực hiện từ tháng 01/2018 đến tháng 08/2018 tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM, mang lại giải pháp công nghệ giúp tối ưu hóa hơn 95% thời gian xử lý số liệu, đảm bảo độ tin cậy tuyệt đối phục vụ công tác xây dựng lưới tọa độ quốc gia và thăm dò địa chất khoáng sản.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên hệ thống lý thuyết trắc địa vật lý kinh điển với định luật vạn vật hấp dẫn của Newton năm 1666, sử dụng hằng số hấp dẫn chuẩn G = 6.67428 x 10^-11 m3/(kg.s2), kết hợp với lý thuyết trường thế trọng trường Trái Đất và bài toán biên xác định Geoid theo công thức Stokes năm 1849 cùng công thức Vening Meinesz. Khung lý thuyết tích hợp 2 mô hình cân bằng đẳng tĩnh vỏ Trái Đất tiêu biểu: mô hình Pratt - Hayford với giả thuyết các khối vỏ có mật độ thay đổi nằm trên mặt bù trừ đồng nhất ở độ sâu D = 100 km, và mô hình Airy - Heiskanen với giả thuyết vỏ Trái Đất có mật độ đồng nhất rho0 = 2.67 g/cm3 nổi trên lớp lót manti dầy đặc có mật độ rho1 = 3.27 g/cm3 cùng độ chênh lệch mật độ delta rho = 0.6 g/cm3 và độ dày lớp vỏ trung bình T = 30 km.

Hệ thống khái niệm chuyên ngành cốt lõi trong mô hình nghiên cứu bao gồm:

  1. Thế trọng lực và dị thường trọng lực: Hiệu số giữa trọng lực thực g và trọng lực chuẩn gamma trên Ellipsoid.
  2. Hiệu chỉnh khoảng không: Quy đổi trọng lực từ điểm đo bề mặt về Geoid trong môi trường chân không với gradient chuẩn -0.3086 mGal/m.
  3. Hiệu chỉnh Bouguer: Triệt tiêu lực hấp dẫn của đĩa phẳng vô hạn dày H với hệ số 0.1119 mGal/m khi mật độ đất đá đạt 2.67 g/cm3.
  4. Hiệu chỉnh địa hình: Bù trừ độ lệch thực tế của địa hình gồ ghề xung quanh điểm đo so với bề mặt đĩa Bouguer phẳng, luôn mang giá trị dương.
  5. Hiệu chỉnh đẳng tĩnh địa hình và ảnh hưởng gián tiếp: Phân bố lại khối lượng địa hình vào lòng đất theo nguyên lý cân bằng thủy tĩnh và hiệu chỉnh sự biến dạng thế năng làm dịch chuyển mặt Geoid thành Co-geoid.

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

Nghiên cứu sử dụng nguồn dữ liệu kép gồm mô hình số độ cao DTM (Digital Terrain Model) với độ phân giải lưới không gian từ 6 giây góc đến 1 km và tập số liệu đo trọng lực thực tế biểu diễn dưới 2 định dạng: hệ tọa độ vuông góc phẳng (x, y) và hệ tọa độ trắc địa (B, L). Cỡ mẫu thực nghiệm được thiết lập gồm tập hợp hàng trăm điểm đo trọng lực kết hợp lưới DTM đa điểm trong phạm vi không gian nghiên cứu. Phương pháp chọn mẫu theo mạng lưới ô vuông không gian phân cấp 4 vùng tính toán: vùng trong cùng có bán kính r = 200 m, vùng gần từ 1 km đến 8 km, vùng xa từ 8 km đến 16 km và vùng ngoài cùng đạt bán kính tối đa R = 45 km đến 50 km. Lý do lựa chọn phương pháp phân cấp mạng lưới không gian này nhằm tối ưu hóa độ chính xác giải tích, phản ánh đầy đủ ảnh hưởng hấp dẫn của địa hình lân cận mà vẫn kiểm soát tốt tài nguyên tính toán.

Phương pháp phân tích dựa trên thuật toán giải tích tích phân số của khối trụ tròn và vành khuyên cắt góc, được lập trình cấu trúc hóa bằng phần mềm Matlab nhằm xử lý song song ma trận số liệu lớn. Quá trình nghiên cứu được triển khai theo timeline biểu đồ từ ngày giao nhiệm vụ 26/02/2018, hoàn thiện thuật toán và phần mềm vào ngày 01/08/2018, và bảo vệ hội đồng đánh giá vào ngày 16/08/2018. Toàn bộ thuật toán được thẩm định độ tin cậy bằng phương pháp kiểm chứng chéo độc lập với bảng tính Microsoft Excel.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình xây dựng và thực nghiệm phần mềm GravityReduction trên nền tảng Matlab đã ghi nhận 4 phát hiện khoa học và kỹ thuật quan trọng:

Thứ nhất, chương trình đã tự động hóa hoàn chỉnh quy trình tính toán 5 bước phức tạp, xử lý đồng thời dữ liệu tọa độ phẳng và tọa độ trắc địa. Tốc độ xử lý dữ liệu ma trận DTM và điểm đo trọng lực đạt dưới 5 giây cho toàn bộ khu vực khảo sát, giúp cắt giảm hơn 95% thời gian thao tác so với quy trình tính toán bán tự động trên Excel.

Thứ hai, kết quả so sánh đối chứng chéo giữa chương trình GravityReduction và Microsoft Excel cho thấy độ lệch tuyệt đối giữa hai phương pháp nhỏ hơn 0.0001 mGal, tương ứng với sai số tương đối dưới 0.001%. Độ chính xác số học này khẳng định thuật toán tích phân vành khuyên và cấu trúc mã nguồn Matlab đạt độ tin cậy tiệm cận mức tuyệt đối.

Thứ ba, số liệu thực nghiệm chứng minh hiệu chỉnh địa hình luôn mang giá trị dương, dao động trung bình từ vài mGal đến hơn 15 mGal tại các khu vực núi dốc, đóng góp từ 10% đến 30% vào tổng giá trị dị thường Bouguer hoàn chỉnh. Nếu bỏ qua hiệu chỉnh địa hình, dị thường trọng lực sẽ bị sai lệch nghiêm trọng, làm sai lệch cao độ Geoid tương ứng.

Thứ tư, theo mô hình cân bằng Airy - Heiskanen, độ dày phần chân bù trừ của vỏ Trái Đất ăn sâu vào manti tại lục địa gấp 4.45 lần độ cao địa hình (t = 4.45 H), trong khi phần bù trừ dưới đại dương đạt 2.73 lần độ sâu biển (t' = 2.73 H'). Kết quả tính toán dị thường đẳng tĩnh địa hình tiệm cận giá trị 0 mGal tại các vùng có cấu trúc vỏ Trái Đất đạt trạng thái cân bằng thủy tĩnh.

Thảo luận kết quả

Nguyên nhân chính giúp chương trình đạt độ chính xác cao là việc ứng dụng công thức tích phân phân tích giải tích cho từng phần tử vành khuyên, triệt tiêu hiện tượng kỳ dị toán học tại tâm điểm đo trọng lực. Khi so sánh với công trình của Sujan Bajracharya tại Canada trên khu vực miền núi rộng lớn, kết quả nghiên cứu hoàn toàn tương đồng về mặt nhận định: các phương pháp hiệu chỉnh đẳng tĩnh địa hình theo Airy - Heiskanen kết hợp dữ liệu DTM độ phân giải cao cho độ chính xác xác định Geoid vượt trội so với phương pháp Bouguer đơn giản. Đồng thời, kết quả này củng cố nghiên cứu của Cao Đình Triều tại vùng Yên Châu (Việt Nam), khẳng định việc lựa chọn bán kính vùng trong 200 m và vùng ngoài 45 km cho phép đạt độ chi tiết hiệu chỉnh trọng lực tới 0.0174 mGal.

Trong thực tế xử lý, các tập dữ liệu có thể được trực quan hóa tối ưu thông qua bảng đối chiếu đa chỉ số (tọa độ x, y, cao độ H, hiệu chỉnh khoảng không, hiệu chỉnh Bouguer, hiệu chỉnh địa hình, dị thường đẳng tĩnh) và hệ thống biểu đồ không gian 3D. Biểu đồ 3D cho phép hiển thị lồng ghép bề mặt địa hình DTM với các vector dị thường trọng lực, giúp các nhà địa vật lý dễ dàng phân tích cấu trúc đứt gãy kiến tạo sâu trong lòng đất.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả ứng dụng kết quả nghiên cứu vào thực tiễn ngành Trắc địa - Bản đồ và Địa vật lý tại Việt Nam, 4 giải pháp cụ thể được đề xuất:

  1. Ứng dụng và chuyển giao phần mềm GravityReduction: Cục Đo đạc, Bản đồ và Thông tin địa lý Việt Nam cùng các viện nghiên cứu chuyên ngành cần đưa phần mềm vào quy trình xử lý chuẩn mạng lưới đo trọng lực quốc gia trong giai đoạn 2024-2025, hướng tới mục tiêu chuẩn hóa sai số biên tập dị thường trọng lực toàn quốc dưới 0.05 mGal.

  2. Nâng cấp độ phân giải mô hình độ cao số DTM: Bộ Tài nguyên và Môi trường cần chỉ đạo các đơn vị đo đạc thành lập lưới DTM quốc gia có độ phân giải không gian từ 3 đến 6 giây góc (tương đương kích thước ô lưới 90 m đến 180 m) trong vòng 12 tháng tới, nhằm nâng cao độ chính xác mô hình Geoid cục bộ đạt mức 2 đến 5 cm tại các vùng đồi núi phức tạp.

  3. Mở rộng mô hình mật độ vỏ Trái Đất phi đồng nhất: Nhóm nghiên cứu tại các trường đại học kỹ thuật cần tiếp tục phát triển mô-đun thuật toán trong giai đoạn 2024-2026, cho phép nạp dữ liệu mật độ biến thiên thực tế thay vì cố định giá trị chuẩn 2.67 g/cm3, giúp giảm thiểu sai số dị thường cục bộ thêm 15%.

  4. Xây dựng cơ sở dữ liệu trọng lực số hóa liên ngành: Các đơn vị khảo sát địa chất và dầu khí cần phối hợp thiết lập hệ thống cơ sở dữ liệu không gian Web-GIS đồng bộ trong thời gian 18 tháng, kết nối 100% dữ liệu từ hơn 10.000 điểm đo trọng lực trên toàn lãnh thổ phục vụ nghiên cứu biến dạng vỏ Trái Đất và phòng chống thiên tai.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và thực tiễn sâu sắc, phù hợp cho 4 nhóm đối tượng trọng tâm:

Thứ nhất, kỹ sư và chuyên viên kỹ thuật ngành Trắc địa - Bản đồ và Địa vật lý: Tài liệu cung cấp quy trình từng bước rõ ràng để tự động hóa xử lý số liệu đo trọng lực thực địa, tiết kiệm 80% thời gian biên tập bản đồ dị thường phục vụ công tác quy chuẩn lưới khống chế độ cao và thăm dò khoáng sản.

Thứ hai, giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Khoa học Trái Đất: Luận văn là tài liệu tham khảo giảng dạy mẫu mực với hệ thống hơn 20 công thức toán học tường minh về thế trọng trường, phương pháp tích phân vành khuyên và các mô hình đẳng tĩnh kinh điển.

Thứ ba, lập trình viên và chuyên gia phát triển phần mềm Trắc địa - GIS: Cung cấp thuật toán tối ưu trên Matlab để kế thừa, phát triển thành các gói phần mềm thương mại hoặc tiện ích mã nguồn mở có khả năng xử lý ma trận dữ liệu địa hình lớn hàng chục nghìn điểm trong vài giây.

Thứ tư, cơ quan quản lý nhà nước và viện nghiên cứu địa chính: Nguồn tư liệu khoa học chuẩn xác phục vụ việc hoạch định chính sách xây dựng hệ quy chiếu không gian quốc gia, xác định Geoid địa phương với độ chính xác cao tiệm cận mức 1 decimet.

Câu hỏi thường gặp

  1. Hiệu chỉnh khoảng không (Free-air reduction) có ý nghĩa vật lý như thế nào? Hiệu chỉnh khoảng không tính chuyển giá trị trọng lực từ mặt đất về mặt Geoid trong giả định không có vật chất tồn tại ở khoảng giữa. Với hệ số biến thiên chuẩn 0.3086 mGal trên mỗi mét cao độ, điểm đo ở độ cao 100 m sẽ có số hiệu chỉnh khoảng không tương ứng là +30.86 mGal.

  2. Tại sao số hiệu chỉnh địa hình (Terrain correction) luôn mang dấu dương? Khối đất đá nhô cao hơn trạm đo có lực hút hướng lên trên làm giảm giá trị trọng lực đo được, do đó khi dời đi ta phải cộng thêm vào. Ngược lại, phần thung lũng khuyết thiếu hụt khối lượng cũng làm giảm lực hút, khi bù đất vào ta tiếp tục phải cộng thêm, khiến số hiệu chỉnh địa hình luôn mang giá trị dương.

  3. Sự khác nhau cơ bản giữa mô hình đẳng tĩnh Pratt - Hayford và Airy - Heiskanen là gì? Mô hình Pratt - Hayford giả định độ sâu bù trừ cố định D = 100 km với mật độ cột đất đá thay đổi tùy độ cao. Ngược lại, mô hình Airy - Heiskanen giả định mật độ vỏ cố định 2.67 g/cm3 nhưng độ sâu chân núi thay đổi, ăn sâu vào lớp manti theo tỷ lệ t = 4.45 H.

  4. Vì sao dị thường Bouguer đơn giản không được sử dụng trực tiếp để tính độ cao Geoid theo công thức Stokes? Việc loại bỏ đĩa Bouguer làm thay đổi thế trọng trường và làm dịch chuyển vị trí mặt Geoid thực một khoảng rất lớn, gấp khoảng 10 lần cao độ Geoid ban đầu. Do đó, dị thường Bouguer đơn giản không thỏa mãn điều kiện biên bảo toàn khối lượng của bài toán Stokes.

  5. Chương trình tính toán viết bằng Matlab có ưu thế gì vượt trội so với Microsoft Excel? Chương trình Matlab xử lý tự động toàn bộ ma trận dữ liệu DTM không gian lớn, tích hợp giao diện trực quan 3D và tự động hóa 100% quy trình 5 bước. Tốc độ tính toán nhanh hơn 95% trong khi vẫn đảm bảo độ lệch kết quả so với Excel nhỏ hơn 0.0001 mGal.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện cơ sở lý thuyết về trọng lực, thế hấp dẫn và 5 phương pháp hiệu chỉnh trọng lực hiện đại phục vụ xác định Geoid.
  • Xây dựng thành công phần mềm GravityReduction trên nền tảng ngôn ngữ Matlab với giao diện đồ họa thân thiện, hỗ trợ đầy đủ hệ tọa độ phẳng và tọa độ trắc địa.
  • Kết quả thực nghiệm khẳng định thuật toán đạt độ chính xác vượt trội với sai số so với tính toán kiểm chứng Excel nhỏ hơn 0.0001 mGal, đồng thời rút ngắn hơn 95% thời gian xử lý số liệu.
  • Làm sáng tỏ cơ chế bù trừ đẳng tĩnh Airy - Heiskanen với tỷ lệ chân núi lục địa đạt 4.45 lần cao độ địa hình và tỷ lệ bù trừ biển đạt 2.73 lần độ sâu.
  • Định hướng giai đoạn 2024-2026 tập trung tích hợp mô hình mật độ biến thiên cục bộ và kết nối cơ sở dữ liệu Web-GIS phục vụ đo đạc bản đồ quốc gia.

Các đơn vị nghiên cứu, cơ quan khảo sát trắc địa và quý độc giả quan tâm hãy ứng dụng ngay quy trình thuật toán này để nâng cao độ chính xác cho các dự án đo đạc trọng lực và xây dựng Geoid thực tế!