毛 賀,盧昕婷
(1.吉林省地質調查院,吉林 長春 130102;2.長春工程學院,吉林 長春 130012)
傳統的磁異常處理手段包括:磁異常化極、解析延拓、磁源重力異常、磁異常分量換算、加深區域內礦石磁性研究等。這些方法將復雜的磁異常簡單化,或者突出磁異常的某一特點,使得異常更加容易解釋。這些方法都有自己的特點和應用條件,并且使用都受到一定的限制。比如化極需要地磁場和地質體磁化強度矢量的偏角和傾角信息。一般來說,地質體都具有顯著的剩余磁化強度,其和感應磁化強度往往還不是同一方向。在這種情況下,地質體磁化強度矢量的方向可能是未知和不準確的。低磁緯度礦區化極計算同樣也不穩定,其會產生平滑但部分錯誤的異常[1]。解析延拓使用比測量磁場高一個數量級的磁位導數,其異常圖分析可能會遺漏一些深度磁源異常引起的低頻信息。解析延拓產生的異常還取決于地質體磁化矢量和周圍地磁場的方向。磁源重力異常的問題和類似化極處理。另外,選擇磁位而不是測量磁場進行解釋,會降低分辨率,特別是針對淺部位置礦體的異常,因為計算磁位會模糊高頻信息。
在文獻[2]中,提出了級磁變換MMTs。MMTs包括TMA本身Ta,以及基于TMA的四種變換:TMA梯度的模R、TMALaplacian算子結果L、E為TMA平方的Laplacian算子結果再開方取二分之一、Q為TMA和TMALaplacian算子結果的乘積的平方根。上述四種變換結果具有非負的分布,對磁化強度矢量方向的敏感程度比TMA自身低得多,并且它們對磁異常源能顯示出比測量的磁異常場更高的中心性。本文使用MMTs對吉林省某礦區的磁異常數據進行處理分析,MMTs得到的異常更接近磁異常源的真實水平位置。

本文使用MMTs對吉林省某礦區的磁異常數據進行處理,測區ΔT平面圖如圖1a所示。
由圖1a可以看出,磁異常主要集中在該研究礦區的南部,異常走向主要有北西向和北東向,個別異常沿近東西向展布,多呈團塊狀或帶狀分布;測區的中部偏南主要是低緩的負異常,變化梯度小;測區的北部主要以大面積低緩正異常為主,局部疊加低緩負異常。
MMTs處理結果如圖1b至圖1d所示。在Ta等值線平面圖中,正異常向北側移動,異常的線性特征相對突出、明顯。Q、E、R等值線圖的異常分布與Ta相似,異常范圍邊界呈遞進收斂關系,位置更加清晰、明確,也更接近磁異常源的真實水平位置。L等值線圖中的異常解比較分散,難以準確的確定異常礦體的邊界位置,但其對淺部異常礦體有一點的指示作用。


圖1
本文使用MMTs對吉林省某礦區磁測數據進行處理,取得了較好的處理效果。
MMTs和傳統磁異常處理方法相比較具有以下優勢:MMTs對某些噪聲的放大作用小于解析延拓;在低磁緯度礦區比化極穩定。此外,MMTs算法簡單,程序實現方便,具有很高的計算效率。