幾種常用地面磁測數據濾波處理方法的編程實現

摘要：由于地面磁法勘探工作數據采集點距離地面較近，受干擾情況較多，尤其是地表不均勻磁性體的影響，導致數據不規律跳動，直接成圖不利于數據的進一步分析。該文利用VB語言對磁法勘探中常用的數據濾波處理方法進行編程實現自動化處理，提高了磁法勘探的解釋能力。

關鍵詞：VB編程 磁法勘探 數據處理 濾波

中圖分類號：P631.2文獻標識碼：A 文章編號：1674-098X（2013）04（c）-0049-02

Visual Basic（簡稱VB）是Microsoft公司推出的一種Windows應用程序開發工具。是當今世界上使用最廣泛的編程語言之一，它也被公認為是編程效率最高的一種編程方法。無論是開發功能強大、性能可靠的商務軟件，還是編寫能處理實際問題的實用小程序，VB都是最快速、最簡便的方法。從設計新型的用戶界面到利用其它應用程序的對象；從處理文字圖象到使用數據庫；從開發個人或小組使用的小工具，到大型企業應用系統，甚至通過Internet的遍及全球分布式應用程序，都可在VB提供的工具中各取所需。[1-2]

磁測工作由于工作效率高、成本低，已成為我國金屬礦產勘查最主要的物探方法之一。如何提高磁法數據處理的自動化尤為重要，通過VB編程可以輕松實現磁測數據的一維、二維分析，可以有效提高磁法勘探的分析能力。

1 磁法勘探數據處理簡介

磁法勘探獲得的原始數據為測點的總場強度T，原始數據經日變改正、高度改正、水平梯度改正和正常場改正等初步整理后得到的為磁異常ΔT數據，之后便可以利用ΔT數據初步成剖面圖、等值線圖等圖件。但以上處理工作滿足不了現有的勘探要求，往往需要對磁測數據進行進一步處理。常用的處理手段有數據圓滑、一維濾波、二維濾波等，然后進行位場分析、反演等工作[3-5]。數據濾波處理的好壞直接影響后續工作的精度，因此濾波處理是一項基礎但又十分重要的工作。

2 數據處理方法及編程實現

2.1 剖面滑動平均法

行業內廣泛使用的半定量反演方法如經驗切線法及2.5維人機交互反演對剖面數據的質量要求較高，而地面磁測的采集數據往往會受到表層不均勻磁性體的干擾產生不規律的跳動，因此在反演前需對數據進行圓滑處理，一般采用滑動平均法。設有n個測點的磁測剖面，設置一矩形窗口（-M，M），是窗口沿剖面逐點滑動，計算窗口內所有觀測值的平均值。窗口范圍則為參與平均計算的個數。

3 結語

本文介紹的幾個常用的磁測數據濾波處理方法在礦區隨機干擾、表層磁性體干擾較嚴重的情況下能夠進行濾波處理，以利于資料的進一步使用，在多個礦區的磁測數據處理中取得了較好的應用效果。此外，VB編程若能結合matlab軟件進行區域場的轉換等處理，將VB的可視化與matlab強大的數據處理能力結合，定能取到更好的應用效果。

參考文獻

[1]古俊科.VB編程技巧拾零[J].電腦知識，1996（5）：34-37.

[2]陳麗.淺談VB編程常用技巧[J].科技信息（科學教研），2008（5）：84-85.

[3]姜永蘭，姜建利.對原有地面磁測數據再處理的應用效果[J].內蒙古地質，2002（3）：25-30.

[4]裴益軒，郭民.滑動平均法的基本原理及應用[J].火炮發射與控制學報，2001（1）：21-23.

[5]馬國慶，孟令順，杜曉娟，等.磁法數據處理中的擴邊和優化中值濾波方法的研究[J].物探化探計算技術，2010，32（2）.

[6]羅孝寬，郭紹雍.應用地球物理學教程-重磁勘探[M].北京：地質出版社，1990.