毛 賀,盧昕婷
(1.吉林省地質(zhì)調(diào)查院,吉林 長春 130102;2.長春工程學(xué)院,吉林 長春 130012)
傳統(tǒng)的磁異常處理手段包括:磁異?;瘶O、解析延拓、磁源重力異常、磁異常分量換算、加深區(qū)域內(nèi)礦石磁性研究等。這些方法將復(fù)雜的磁異常簡單化,或者突出磁異常的某一特點,使得異常更加容易解釋。這些方法都有自己的特點和應(yīng)用條件,并且使用都受到一定的限制。比如化極需要地磁場和地質(zhì)體磁化強度矢量的偏角和傾角信息。一般來說,地質(zhì)體都具有顯著的剩余磁化強度,其和感應(yīng)磁化強度往往還不是同一方向。在這種情況下,地質(zhì)體磁化強度矢量的方向可能是未知和不準(zhǔn)確的。低磁緯度礦區(qū)化極計算同樣也不穩(wěn)定,其會產(chǎn)生平滑但部分錯誤的異常[1]。解析延拓使用比測量磁場高一個數(shù)量級的磁位導(dǎo)數(shù),其異常圖分析可能會遺漏一些深度磁源異常引起的低頻信息。解析延拓產(chǎn)生的異常還取決于地質(zhì)體磁化矢量和周圍地磁場的方向。磁源重力異常的問題和類似化極處理。另外,選擇磁位而不是測量磁場進(jìn)行解釋,會降低分辨率,特別是針對淺部位置礦體的異常,因為計算磁位會模糊高頻信息。
在文獻(xiàn)[2]中,提出了級磁變換MMTs。MMTs包括TMA本身Ta,以及基于TMA的四種變換:TMA梯度的模R、TMALaplacian算子結(jié)果L、E為TMA平方的Laplacian算子結(jié)果再開方取二分之一、Q為TMA和TMALaplacian算子結(jié)果的乘積的平方根。上述四種變換結(jié)果具有非負(fù)的分布,對磁化強度矢量方向的敏感程度比TMA自身低得多,并且它們對磁異常源能顯示出比測量的磁異常場更高的中心性。本文使用MMTs對吉林省某礦區(qū)的磁異常數(shù)據(jù)進(jìn)行處理分析,MMTs得到的異常更接近磁異常源的真實水平位置。
本文使用MMTs對吉林省某礦區(qū)的磁異常數(shù)據(jù)進(jìn)行處理,測區(qū)ΔT平面圖如圖1a所示。
由圖1a可以看出,磁異常主要集中在該研究礦區(qū)的南部,異常走向主要有北西向和北東向,個別異常沿近東西向展布,多呈團(tuán)塊狀或帶狀分布;測區(qū)的中部偏南主要是低緩的負(fù)異常,變化梯度小;測區(qū)的北部主要以大面積低緩正異常為主,局部疊加低緩負(fù)異常。
MMTs處理結(jié)果如圖1b至圖1d所示。在Ta等值線平面圖中,正異常向北側(cè)移動,異常的線性特征相對突出、明顯。Q、E、R等值線圖的異常分布與Ta相似,異常范圍邊界呈遞進(jìn)收斂關(guān)系,位置更加清晰、明確,也更接近磁異常源的真實水平位置。L等值線圖中的異常解比較分散,難以準(zhǔn)確的確定異常礦體的邊界位置,但其對淺部異常礦體有一點的指示作用。
圖1
本文使用MMTs對吉林省某礦區(qū)磁測數(shù)據(jù)進(jìn)行處理,取得了較好的處理效果。
MMTs和傳統(tǒng)磁異常處理方法相比較具有以下優(yōu)勢:MMTs對某些噪聲的放大作用小于解析延拓;在低磁緯度礦區(qū)比化極穩(wěn)定。此外,MMTs算法簡單,程序?qū)崿F(xiàn)方便,具有很高的計算效率。