Xử lý và minh giải dữ liệu địa từ từ tuyến bay đến bản đồ

Xử lý dữ liệu địa từ biến chuỗi đo theo tuyến bay thành các lớp bản đồ có thể kiểm tra và diễn giải: dữ liệu trên không được đồng bộ với GNSS và trạm từ nền, hiệu chỉnh biến thiên theo thời gian, loại trường tham chiếu IGRF, cân bằng giữa các tuyến, nội suy thành lưới rồi mới tạo TMI, dị thường và những sản phẩm dẫn xuất phù hợp. Quy trình tốt phải giữ đường dẫn ngược về dữ liệu nguồn; nó không tạo lại tín hiệu chưa từng được đo và không tự động phát hiện quặng.

Trong phạm vi tiêu chuẩn hiện hành, BlueCap bàn giao dữ liệu BlueCapBird® và trạm nền đã đồng bộ, quyền truy cập đọc dữ liệu tuyến nguồn, lưới TMI sẵn sàng cho GIS sau hiệu chỉnh biến thiên, loại IGRF và cân bằng tiêu chuẩn, cùng mô tả cấu trúc và hệ tọa độ. Cân bằng chính xác, microlevelling, RTP, analytic signal, tilt, 1VD/2VD, lọc, nghịch đảo và minh giải địa chất là xử lý tùy chọn hoặc theo yêu cầu dự án, trừ khi chúng được ghi rõ trong phạm vi đã thống nhất.

Nếu cần thiết kế và thu nhận dữ liệu mới, xem dịch vụ khảo sát địa từ bằng drone. Nếu cần hiểu cảm biến, trường từ toàn phần và ý nghĩa của một dị thường, đọc hướng dẫn từ kế drone và dị thường từ.

1. Dữ liệu địa từ nguồn cần có trước khi bắt đầu xử lý

Một tệp giá trị nanotesla đơn lẻ chưa đủ để tạo bản đồ địa từ có thể bảo vệ về kỹ thuật. Mỗi phép đo phải được giữ cùng vị trí, thời gian và trạng thái của hệ thống thu nhận. Trên BlueCapBird®, dữ liệu từ kế QuSpin QTFM Gen2 được ghép với các kênh hỗ trợ để người xử lý có thể phân biệt đáp ứng địa chất với lỗi thời gian, thay đổi độ cao, chuyển động nền tảng hoặc đoạn đo không hợp lệ.

Dữ liệu nguồn nên được lưu ở trạng thái chỉ đọc hoặc có phiên bản bất biến. Việc xóa spike, loại đoạn quay đầu hay sửa vị trí phải tạo kênh xử lý mới và ghi lại lý do, không ghi đè lên phép đo ban đầu. Quyền truy cập PostgreSQL chỉ đọc trong gói BlueCap giúp nhóm GIS, Oasis montaj hoặc pipeline riêng kiểm tra lại dữ liệu tuyến chưa lọc và từ điển trường.

BlueCapBird® Fast phục vụ cấu hình thu nhận tốc độ cao và địa hình phức tạp; BlueCapBird® Light phù hợp cấu hình tốc độ thấp hơn trên địa hình đơn giản. Dù dùng cấu hình nào, xử lý không thể bù hoàn toàn cho nhiễu nền tảng, độ cao không ổn định hoặc khoảng trống đã xảy ra khi thu nhận.

2. Đồng bộ thời gian, hiệu chỉnh trạm nền, biến thiên từ và IGRF

Thứ tự xử lý cụ thể phụ thuộc thiết bị và dự án, nhưng các bước phải có thể lặp lại từ dữ liệu gốc. Hiệu chỉnh nhật biến là cách gọi thường gặp của hiệu chỉnh biến thiên từ theo thời gian; nó không phải cùng một phép tính với loại trường chính IGRF.

Trạm từ nền theo dõi thay đổi trường địa từ tại một vị trí tham chiếu trong thời gian khảo sát. Nó hỗ trợ loại phần biến thiên theo thời gian, nhưng không loại được trường riêng của drone, nhiễu đường điện, sai số AGL hoặc tín hiệu từ một vật thể kim loại trên mặt đất. Nếu trạm nền quá xa, bị nhiễu cục bộ, thiếu thời gian hoặc trường biến động mạnh, người xử lý phải ghi rõ giới hạn thay vì coi phép hiệu chỉnh là hoàn hảo.

IGRF là mô tả toán học chuẩn của trường từ chính Trái Đất và biến thiên thế kỷ, không phải bản đồ địa chất địa phương. Khi loại IGRF, hồ sơ xử lý cần ghi thế hệ mô hình, ngày giờ, tọa độ và độ cao dùng để tính. IGRF-14 là thế hệ hiện hành cho giai đoạn này; tên phiên bản vẫn phải được lưu để dữ liệu có thể tái xử lý trong tương lai.

Dự thảo quy định kỹ thuật bay đo từ và phổ gamma của cơ quan quản lý Việt Nam đề cập đồng bộ trường từ với tọa độ và thời gian, theo dõi biến thiên, kiểm tra tuyến và tránh thu nhận trong bão từ. Tài liệu này được dùng làm nguồn tham khảo thuật ngữ và nguyên tắc kỹ thuật, không được trình bày như quy định cuối cùng hoặc tư vấn pháp lý cho một chuyến bay cụ thể.

3. Cân bằng tuyến, crossover levelling và microlevelling khác nhau thế nào?

Sau hiệu chỉnh thời gian và IGRF, hai tuyến đi qua cùng khu vực vẫn có thể không khớp hoàn toàn. Nguyên nhân gồm sai khác độ cao, độ trễ còn dư, nhiễu theo hướng, mức nền giữa các chuyến hoặc biến thiên thời gian chưa được mô tả hết. Cân bằng tuyến (levelling) dùng các giao điểm giữa tuyến chính và tuyến khống chế để ước tính và giảm sai khác có hệ thống.

Microlevelling không đồng nghĩa với làm mượt toàn bộ lưới. Nó cần một kênh hiệu chỉnh riêng, giới hạn biên độ, độ dài lọc và hình ảnh phần tín hiệu đã bị lấy đi. USGS mô tả xử lý dữ liệu từ hàng không tại Santa Cruz cho thấy đạo hàm có thể khuếch đại nhiễu tuyến còn lại; vì thế phải kiểm tra corrugation trước khi dùng các phép tăng cường biên và ước tính độ sâu.

Một bản đồ trông “sạch” hơn chưa chắc trung thực hơn. QC cần so sánh lưới trước và sau cân bằng, histogram mis-tie, profile trên các tuyến đại diện, tương quan với độ cao và lưới hiệu chỉnh đã trừ. Tín hiệu bị xóa phải được xem xét nghiêm túc như tín hiệu được giữ lại.

4. Từ dữ liệu tuyến đến TMI, dị thường từ và lưới GIS

TMI (Total Magnetic Intensity) về nghĩa vật lý là cường độ trường từ toàn phần. Trong thực hành, tên “bản đồ TMI” đôi khi được dùng cho trường đã qua hiệu chỉnh và đôi khi cho phần dư sau khi loại IGRF. Vì vậy, tên tệp không đủ để hiểu lớp dữ liệu: metadata phải nói rõ trạm nền, IGRF, mức tham chiếu, levelling và mọi phép lọc đã áp dụng.

Nội suy chuyển các điểm hoặc profile không đều thành raster có ô lưới đều. Thuật toán, kích thước ô, khoảng tìm kiếm, cách xử lý rìa và vùng không có dữ liệu phải phù hợp với khoảng cách tuyến, khoảng cách điểm dọc tuyến, AGL và bước sóng mục tiêu. Chọn ô lưới rất nhỏ chỉ tạo nhiều pixel hơn; nó không tạo thêm độ phân giải không gian mà hệ thống chưa đo được.

Trước khi phát hành lưới, người xử lý nên kiểm tra tối thiểu:

• profile tuyến nguồn có khớp với raster tại cùng vị trí hay không;
• rìa lưới, lỗ dữ liệu và khu vực ngoại suy có được mask rõ ràng hay không;
• sọc theo hướng tuyến, bullseye quanh điểm lỗi hoặc bước nhảy giữa block có còn xuất hiện hay không;
• kích thước ô và phép nội suy có phù hợp hình học khảo sát hay không;
• giá trị cực trị có phải phép đo hợp lệ, nhiễu văn hóa hay artefact xử lý;
• CRS, datum, đơn vị, nodata và quy ước tên lớp có được ghi trong metadata.

Ví dụ quy trình công khai của USGS cho Royal Society Range đi từ kiểm tra dữ liệu nguồn, hiệu chỉnh trạm nền và loại IGRF tới levelling tại giao điểm rồi mới lập lưới. Đây là một ví dụ kỹ thuật có truy vết, không phải bộ tham số mặc định áp dụng cho mọi dự án BlueCap.

5. RTP, analytic signal, tilt, 1VD, 2VD và các bộ lọc địa từ

Sản phẩm dẫn xuất chỉ có giá trị khi trả lời một câu hỏi địa chất cụ thể. Tất cả các lớp dưới đây là tùy chọn hoặc theo phạm vi dự án; BlueCap không mặc định tạo một chồng bản đồ giống nhau cho mọi khu đo.

RTP không phải luôn là lựa chọn tốt nhất cho dữ liệu ở Việt Nam. Ở vĩ độ từ thấp, phép biến đổi có thể kém ổn định; từ dư mạnh cũng làm giả định cảm ứng đơn giản không còn phù hợp. Khi đó analytic signal, tilt, gradient, mô hình forward hoặc một cách RTP ổn định có thể hữu ích hơn, nhưng quyết định phải dựa trên dữ liệu và mô hình địa chất chứ không phải tên sản phẩm phổ biến.

Đạo hàm và bộ lọc không phát hiện vật chất mới. Chúng chỉ biểu diễn lại nội dung phổ và gradient đã có trong lưới đầu vào. Nếu dữ liệu nguồn thiếu tín hiệu vì cảm biến bay quá cao, khoảng cách tuyến quá rộng hoặc đã bị làm mượt trước đó, 1VD và 2VD không thể phục hồi thông tin đã mất.

6. Nghịch đảo từ 3D và minh giải địa chất không phải phát hiện quặng tự động

Nghịch đảo dữ liệu từ tìm một mô hình phân bố độ từ cảm hoặc tham số nguồn có thể tái tạo trường quan sát trong giới hạn sai số. Bài toán này không duy nhất: nhiều hình học, độ sâu, độ từ cảm và từ dư khác nhau có thể tạo đáp ứng gần giống nhau. Kết quả phụ thuộc lưới đầu vào, topography, kích thước ô mô hình, regularisation, mô hình nền, ràng buộc địa chất và giả định về hướng từ hóa.

• 01

  Quan sát đã đo

  Profile trường từ, vị trí, độ cao, trạm nền, tuyến và chỉ số QC là bằng chứng gốc cần giữ nguyên.

• 02

  Sản phẩm đã xử lý

  TMI, dị thường, RTP, gradient và đạo hàm là phép biến đổi có tham số, không phải một quan sát độc lập mới.

• 03

  Mô hình địa vật lý

  Forward model hoặc nghịch đảo 3D thể hiện một hay nhiều phân bố thuộc tính vật lý phù hợp dữ liệu và ràng buộc.

• 04

  Giả thuyết địa chất

  Đơn vị đá, tiếp xúc, đứt gãy, xâm nhập hoặc hệ khoáng hóa được đề xuất bằng cách tích hợp nhiều lớp bằng chứng.

• 05

  Kiểm chứng hiện trường

  Lập bản đồ địa chất, địa hóa, phương pháp địa vật lý khác, trench, mẫu và khoan kiểm tra giả thuyết.

• 06

  Quyết định tài nguyên

  Hàm lượng, tính liên tục, tài nguyên và trữ lượng cần dữ liệu khoan, assay, QA/QC và quy trình chuyên môn phù hợp.

Minh giải có thể phân chia miền trường, theo dấu lineament, ước tính vị trí tiếp xúc, so sánh dị thường với bản đồ địa chất và xếp hạng mục tiêu khảo sát tiếp. Với từng loại khoáng sản, vai trò của địa từ khác nhau; ma trận tại địa vật lý thăm dò khoáng sản theo hệ quặng giải thích riêng cho sắt, vàng, đồng, niken, lithium và đất hiếm.

Không nên chuyển thẳng một màu nóng trên bản đồ thành “mục tiêu quặng”. Bản đồ từ không đo trực tiếp vàng, đồng hoặc hàm lượng khoáng sản; mô hình nghịch đảo không phải mô hình tài nguyên; AI hoặc pipeline tự động không thay thế chuyên gia chịu trách nhiệm, khảo sát địa chất, lấy mẫu, khoan và assay. Tự động hóa hữu ích cho các bước lặp lại và kiểm tra nhất quán, còn giả thuyết địa chất và độ bất định phải được chuyên gia xem xét.

7. Bàn giao GIS, hồ sơ QC và phạm vi xử lý BlueCap

Phạm vi xử lý phải được chốt trước khi bắt đầu để nhóm dự án biết lớp nào là sản phẩm tiêu chuẩn, lớp nào là dịch vụ chuyên môn bổ sung và lớp nào chỉ là bản xem trước hiện trường.

Các định dạng có thể thống nhất gồm GeoTIFF cho raster GIS, GeoPackage cho vector và lớp dự án, grid Geosoft cho quy trình địa vật lý, cùng PDF/PNG cho bản đồ trình bày. Định dạng cuối, hệ tọa độ, độ phân giải, nodata, bảng màu, tên lớp và bộ sản phẩm phải được ghi trong phạm vi; việc một định dạng có thể xuất không có nghĩa nó mặc định có trong mọi hợp đồng.

Checklist QC trước khi phát hành dữ liệu địa từ

1. Dữ liệu nguồn được khóa phiên bản, có checksum hoặc lịch sử phiên bản và không bị ghi đè.
2. Tất cả kênh dùng cùng quy ước UTC, đơn vị, CRS và mốc độ cao đã công bố.
3. Trạm nền, khoảng thời gian hợp lệ, biến thiên mạnh, gap và phép nội suy theo thời gian đã được kiểm tra.
4. Các kênh hiệu chỉnh — diurnal, IGRF, levelling và tùy chọn microlevelling — được lưu hoặc mô tả riêng.
5. Mis-tie tại crossover, đoạn bay lại, mask, nhiễu văn hóa và phụ thuộc AGL được ghi trong hồ sơ QC.
6. Raster được đối chiếu lại với profile; cell size, thuật toán nội suy, rìa và vùng ngoại suy được công bố.
7. Mỗi dẫn xuất và mô hình có tham số, phiên bản phần mềm, giả định và giới hạn diễn giải.
8. Sản phẩm cuối được chuyên gia xem xét theo câu hỏi dự án trước khi dùng cho xếp hạng mục tiêu.

Nguồn kỹ thuật chính

• Quy định kỹ thuật đo từ của Việt Nam trên Cơ sở dữ liệu quốc gia về văn bản pháp luật — thuật ngữ trường từ toàn phần, trường bình thường và trường dị thường.
• Dự thảo kỹ thuật bay đo từ và phổ gamma bằng phương tiện bay — nguồn tham khảo về tuyến, biến thiên từ, đồng bộ, kiểm tra và xử lý sơ bộ; đây là dự thảo, không phải tư vấn pháp lý.
• IAGA / NCEI — International Geomagnetic Reference Field — mô hình trường từ chính, hệ số IGRF-14 và yêu cầu nhận diện đúng thế hệ mô hình.
• USGS — Royal Society Range Aeromagnetic Data Processing — ví dụ công khai về QC, trạm nền, IGRF, levelling tại giao điểm và gridding.
• USGS — Processing and Interpretation of Aeromagnetic Data, Santa Cruz Basin — ảnh hưởng của nhiễu tuyến tới đạo hàm, RTP, analytic signal và ước tính nguồn.
• QuSpin QTFM specifications — thông số nguồn của từ kế nguyên tử rubidi đo trường từ toàn phần dùng trong BlueCapBird®.

Các nguồn trên mô tả nguyên lý và ví dụ kỹ thuật. Tham số áp dụng cho một dự án BlueCap được quyết định từ hình học khảo sát, dữ liệu thực tế, mục tiêu địa chất và phạm vi bàn giao đã thống nhất.

Để chuẩn bị đầu vào đúng từ đầu, xem hướng dẫn lập kế hoạch tuyến bay khảo sát drone và cấu hình BlueCapBird® Fast. Công cụ dự toán chi phí khảo sát drone hỗ trợ ước tính line-km, thời gian và ngân sách ban đầu.