基于Oasis Montaj平臺的航空物探數據處理軟件系統

胥值禮，孟慶敏，李文杰，崔志強

(中國地質科學院 地球物理地球化學勘查研究所，廊坊 065000)

0 前言

航空物探具有效率高、成本較低、適用于大面積快速勘查等優點，它可為礦產資源、油氣資源、海洋地質、地下水資源、工程地質與環境，基礎地質、軍事與國防等方面提供大量信息。近年來在國民經濟建設中發揮了重要的作用，得到了快速發展和廣泛應用。隨著硬件和軟件技術的發展，高分辨率、超高分辨率航空物探技術也得到快速發展。因此相應的數據處理軟件面臨著大數據量、不同采樣率、多參數、不同存儲方式等問題，這就對數據處理的效率和質量提出了更高要求。

針對上述問題，本次研究對現有的軟件(多為使用 FORTRAN 語言編寫的獨立程序，且不具備直觀的軟件操作界面)統一使用C/C++語言進行了移植、優化，并開發了新的功能，集成為一套能處理航空電磁法、航空磁法、航空放射性數據的航空物探數據處理軟件系統，以滿足越來越多、越來越高的社會需求。

本系統選用的開發平臺為加拿大 Geosoft公司的Oasis Montaj(簡稱 OM)基本軟件平臺。它是目前在地球科學領域中應用最廣泛的軟件平臺，是最著名的礦產勘查和地球科學調查領域的領先軟件，得到了各著名礦業公司、油氣公司、勘探公司、咨詢公司、科研院所和政府部門的一致認可和長期使用。它具有高性能的數據庫系統，可高效地存儲和處理多達 64T 字節的數據，能為多種地球物理勘查海量數據提供強有力的存儲管理支撐。它提供的 GX Developer 二次開發軟件包，具有數據庫操作、數據濾波、數據成圖等諸多基本功能，可進行應用功能的快速定制開發[6]。因此采用這種開發模式，可將開發重點集中在對專業數據處理方法的研究上，適合于多種數據處理軟件功能的開發，具有較強的靈活性，能大大地提高開發效率和開發質量、節約開發成本。

1 系統設計

1.1 開發環境

本系統的開發環境為：Windows XP 及以上版本的操作系統、加拿大 Geosoft 公司的 Oasis montaj V7.2 基本軟件平臺、Visual C++ 6.0 軟件編程平臺、C/C++ 及 GXC(Geosoft 的 GX 編程語言)程序語言。

1.2 系統架構

本系統采用三層結構：第一層為Windows XP操作系統；第二層為 OM 平臺(包括數據庫系統)；第三層為航空物探數據處理應用軟件系統(圖1)。

圖1 系統架構圖Fig．1 The diagram of system architecture

1.3 系統功能

本系統的主要功能分為：數據管理、航磁數據處理、航電數據處理、航放數據處理、數據調平處理五大類(圖2)。

(1)數據管理功能。它包括數據輸入、數據合并、數據輸出、數據裁剪模塊。該功能適用于航電、航磁、航放數據。

(2)航磁數據處理功能[7-10]。它包括日變數據預處理、日變校正、滯后校正、方向差校正、正常場校正、航磁總精度統計模塊。

(3)航電數據處理功能，包括視電阻率轉換模塊。

(4)航放數據處理功能[8-9]，包括氣壓高度轉換、數據配置、數據預處理、背景消除、含量計算、比率計算模塊。

(5)數據調平處理功能。它包括基于測線數據的調平(一維人機交互調平[1]、二維自動調平[2-3])、基于網格數據的二維自動調平模塊，適用于航電、航磁、航放數據。航磁實測切割線數據調平模塊，適用于航磁數據。

1.4 模塊開發

如圖3所示，應用功能模塊自上而下的調用流程為：OM 平臺菜單-> GX 接口模塊->航空物探應用軟件模塊-> OM 開發庫模塊(GX API)。在本系統中，開發的模塊為：GX 接口模塊、航空物探應用軟件模塊。

(1)OM 平臺菜單。在 OM 平臺上動態加載的菜單，由它調用 GX 接口模塊。

(2)GX 接口模塊。在 OM 平臺上開發的菜單調用接口模塊，由它調用航空物探應用軟件模塊。

(3)航空物探應用軟件模塊。在Visual C++ 6.0 平臺上開發的動態鏈接庫(DLL)或可執行模塊(EXE)，由它們調用 OM 平臺的 GX API 開發庫模塊，完成航空物探數據處理功能。它是本系統的核心模塊，提供參數設置用戶界面和數據處理功能。

(4)OM開發庫模塊(GX API)。OM 平臺提供的二次開發庫，由它完成對數據庫(簡稱 GDB)的操作和一些基本處理功能(如數據濾波、圖形生成與顯示等)。

1.5 軟件編譯

(1)GX 接口模塊由 GX Developer Toolkit 7.3.0 I8 提供的編譯器進行編譯。

(2)航空物探應用軟件模塊由 Visual C++ 6.0 提供的編譯器進行編譯。

圖2 系統功能結構圖Fig．2 The diagram of system function structure

圖3 軟件功能調用流程圖Fig．3 The flow chart of system function

1.6 系統集成

圖4為系統集成主界面示意圖。圖2中的功能模塊可以“插件”方式動態加載到 OM 平臺上，加載的主菜單分別為AGS_Data(數據管理)、AGS_Mag(航磁數據處理)、AGS_EM(航電數據處理)、AGS_Rad(航放數據處理)和AGS_Lev(數據調平處理)。各個主菜單分別由其相應的功能子菜單組成。

2 系統特點及主要優勢

本系統具有以下幾個主要技術特色。

2.1 剖面平面圖方式的數據顯示與處理功能

利用地理信息系統(GIS)的圖層管理模式[4] [5]，以剖面平面圖方式(剖面平面曲線圖、剖面平面顏色漸變圖，分別見圖4的左下、右下界面)動態直觀地顯示和處理測線數據，克服了單一剖面顯示方式不能觀察測區整體情況、不能即時觀察處理效果以及切換測線的重復操作而導致處理效率低下等問題。此“所見即所得”方式便于數據的選擇性處理和人機交互分析，可極大地提高數據處理的質量和效率。

2.2 數據調平處理功能

(1)基于測線數據的自動二維調平功能。采用了二維空間移動平均濾波和二分法搜索等算法，來消除零漂、基值不一致對數據造成的影響，以改善成圖質量，突出目標異常。

(2)基于網格數據的自動二維調平功能。采用了規則網格數據運算和二維空間褶積濾波等算法，對數據進行精細調平處理，消除局部不水平和隨機噪聲的影響，以改善成圖質量，突出目標異常。

2.3 數據批處理功能

本系統開發的宗旨是面向生產應用。在設計時，充分考慮了航空物探的數據特點和生產中用戶提出的需求，實現了按飛行時間段或測量分區對數據進行自動批處理功能。可分別保存不同采集系統、不同時間段的處理參數，克服了處理參數輸入的繁雜性和重復性，大大提高了數據處理的效率與質量，具有較強的實用性。

3 應用實例

3.1 航空電磁數據處理

圖5為北方某區域航電數據的處理結果。圖5(a)為使用視電阻率轉換功能對原始數據(固定翼飛機采集)轉換得到的視電阻率，水平不一致和噪聲影響較為嚴重；圖5(b)為二維自動調平處理后的視電阻率，較好地消除了水平不一致和噪聲的影響，從而使視電阻率的分布特征更加清晰。

圖4 系統集成主界面示意圖Fig．4 The main interface of integrated system

圖5 北方某區域航電視電阻率處理圖(中頻1563 Hz)Fig．5 The processed map of apparent resistivity from an airborne electromagnetic survey in northern China (1563Hz)(a)原始數據轉換的視電阻率；(b)二維自動調平后的視電阻率

圖6 南方某區域航磁數據處理圖Fig．6 The processed map of △T from an airborne magnetic survey in southern China(a)原始數據；(b)常規處理結果；(c)2D自動調平結果

圖7 南方某區域航放鈾元素數據處理圖Fig．7 The processed map of uranium element from an airborne radioactivity survey in southern China(a)原始鈾計數率；(b)常規處理得到的鈾含量；(c)基于網格數據二維自動調平的鈾含量

3.2 航空磁測數據處理

圖6為南方某區域航磁數據(直升飛機采集)的處理結果。圖6(a)為原始數據(總場△T)；圖6(b)為常規處理(日變數據預處理、日變校正、滯后校正、方向差校正和正常場校正)的結果，從圖中可以看出，數據的水平不一致現象和隨機噪聲干擾較明顯；圖6(c)為先經基于測線數據的二維自動調平，再經基于網格數據的二維自動調平的結果，較好地消除了圖6(b)中的水平不一致及隨機噪聲的影響，清晰地顯示出了異常的特征。

3.3 航空放射性數據處理

圖7為南方某區域航放鈾元素數據(直升飛機采集)的處理結果。圖7(a)為鈾元素的原始計數率；圖7(b)為常規處理(飛機本底消除、宇宙本底消除、大氣氡本底消除、康普頓譜剝離、高度衰減影響消除和放射性元素含量計算)的鈾含量；圖7(c)為基于網格數據二維自動調平的鈾含量，較好地消除了水平不一致及噪聲的影響，使異常形態更加清晰。

4 結束語

本系統可同時處理航空電磁法、磁法、放射性數據，具有集成度高、功能強、實用性強、可視化程度高、運行效率高、運行穩定可靠、操作簡便、容錯性強等諸多優勢，尤其是剖面平面圖方式的數據顯示與處理、數據調平處理、數據批處理功能，可極大地提高數據處理的效率與質量。對大量實測數據的處理結果表明，所實現的功能能很好地滿足航空物探的數據處理要求。

在本系統基礎上，可對功能進行逐漸增加和完善，使之成為一套高效適用的集測線規劃、質量分析、數據處理和綜合解釋為一體的航空物探軟件系統，以完成航空物探的一系列工作流程，更好地服務于生產應用。

參考文獻：

[1] 胥值禮，劉還林，丁志強，等.頻率域航空電磁法數據調平軟件設計與實現[J].物探化探計算技術，2009,31(4):333-338.

[2] 胥值禮，孟慶敏，崔志強.基于GIS的航磁實測切割線網自動調平軟件[J].物探化探計算技術，2010，32(4):423-428.

[3] 李文杰.用于頻率域航空電磁數據的二維自動調平[J].成都理工大學學報:自然科學版，2007,34(4):447-451.

[4] 于長春，郭志宏,眭素文，等.航空物探領域的GIS開發與應用[J].物探化探計算技術，2003,25(1):39-44.

[5] 陳建春. Visual C++開發GIS系統[M]. 北京：電子工業出版社，2000.

[6] Geosoft Inc. GX Developer Custom Application Development Toolkit for extending Oasis montaj and standalone applications USER GUIDE and REFERENCE MANUAL (V7.0)[S], 2009.

[7] 中國國土資源部發布.航空磁測技術規范[S].2010.

[8] 核工業航測遙感中心譯.IAEA航空伽瑪能譜測量323規定[S].1991.

[9] 中國國防科學技術委員會發布.航空伽瑪能譜測量規范[S].2005.

[10] 穆石敏，申寧華，孫運生.區域地球物理數據處理方法及其應用[M].吉林:吉林科學技術出版社，1990.