面向深部礦產資源勘探的井中重磁三維交叉梯度聯合反演

徐凱軍 季春暉 龐朝輝

摘要：單一井中重磁方法解決問題有限，其三維反演具有嚴重的多解性。開展面向深部礦產資源勘探的井中重磁三維交叉梯度聯合反演，采用交叉梯度函數對密度和磁化率進行耦合，綜合井中重磁多參數信息進行聯合反演，減少反解的多解性，提高深部礦體識別的準確性。通過理論模型和實際資料進行反演分析。結果表明：井中重磁交叉梯度聯合反演具有較高的縱向分辨能力，可以顯著提高深部礦體識別的準確性，特別是能有效識別深部隱伏礦體。

關鍵詞：井中重力； 井中磁力； 交叉梯度； 三維聯合反演

中圖分類號：P 631.1?? 文獻標志碼：A

文章編號：1673-5005（2024）03-0065-10?? doi：10.3969/j.issn.1673-5005.2024.03.007

Cross-gradient joint inversion of 3D borehole gravity and magnetic data for deep mineral resources exploration

XU Kaijun1， JI Chunhui2， PANG Zhaohui1

（1.School of Geosciences in China University of Petroleum（East China）， Qingdao 266580， China；2.Hainan Branch Company of CNOOC Limited， Haikou 571100， China）

Abstract：The single borehole gravity and magnetic method has its limitations， and 3D inversion techniques often suffer from non-uniqueness. To address these challenges in deep mineral resources exploration，we propose a novel approach： borehole gravity and magnetic 3D cross gradient joint inversion.Initially， we couple density and magnetic susceptibility using cross-gradient functions， followed by joint inversion based on borehole gravity and magnetic multi-parameter data. This approach aims to reduce ambiguity and enhance the identification accuracy of deep mineral resources. The theoretical model and real data inversion results demonstrate that the borehole gravity and magnetic cross-gradient joint inversion offers high vertical resolution， significantly improving the identification accuracy of deep mineral resources， particularly for effectively identifying deep concealed ore bodies.

Keywords：borehole gravity data; borehole magnetic data; cross gradient; 3D joint inversion

目前中國礦產資源勘探的目標已轉向深部礦產資源［1］，最為關鍵的是提高深部目標礦體的識別精度。井中重磁測量可以近距離感知目標礦體，顯著提高縱向分辨率且不易受淺地表干擾源的影響。井中重力測量儀器在國外已得到很好的應用［2-3］，中國目前正攻關井中高精度重磁勘探技術與裝備研發，應用后可滿足3000 m井深探測需求，提高深部礦產資源勘探能力。數據處理［4］和反演［5-6］是解釋的基礎，國內外學者在井中重磁反演理論方面已開展了相關研究，取得了較好的進展［7-12］。由于單一地球物理數據反演存在嚴重的多解性，為了提高反演的可靠性，不同地球物理數據聯合反演得到了快速發展［13-16］。針對地面重磁數據的聯合反演也得到廣泛應用［17-18］，考慮到井中重磁數據縱向分辨率高，井中重磁數據聯合反演可以更好地提高反演可靠性和深部目標體識別精度。由于密度和磁化率之間沒有明確的巖石物理關系，實現基于交叉梯度耦合的井中重磁數據聯合反演，為深部礦產資源勘探開發提供技術支持。

1 井中重磁異常正演理論

井中重磁正演通常將地下空間剖分為一系列長方體網格，采用不規則網格，礦體區域網格剖分較為細密，邊界處較為稀疏。每一個小網格可視為均勻密度體和磁性體，單元網格在井中觀測點產生的重磁三分量異常表達式如下：

gx（x，y，z）=Gσv（ε-x）r3dv，gy（x，y，z）=Gσv（η-y）r3dv，gz（x，y，z）=Gσv（ζ-z）r3dv.（1）

式中gx、gy和gz分別為x、y、z三個方向的重力異常；G為引力常量；σ為網格單元剩余密度。

Bx=μ04π（MxVxx+MyVxy+MzVxz），By=μ04π（MxVyx+MyVyy+MzVyz），Bz=μ04π（MxVzx+MyVzy+MzVzz）.（2）

其中

V=∫v1rdv.

式中，Bx、By、Bz分別為x、y、z三個方向的磁異常；μ0為真空磁導率；Mx、My、Mz分別為網格單元x、y、z三個方向的磁化強度。

基于礦體分布形態，不同單元網格可以取不同的密度和磁化強度值，依次計算地下所有剖分網格在井中觀測點產生的重磁異常進行求和即可得到井中觀測點的正演響應。

2 井中重磁聯合反演

2.1 單一井中重磁數據反演方法

井中重磁異常反演是一種求取最優化目標解的過程，反演的目標函數為

φ（m）=φd（m）+βφm（m）=Wd（Gm-d）22+βWm（m-m0）22.（3）

式中，φ（m）為目標函數；φd（m）為數據擬合目標函數；φm（m）為模型目標函數；d為觀測數據向量；G為雅克比矩陣；Wd為數據對角加權矩陣；m為模型向量；Wm為模型加權矩陣；m0為參考模型向量。反演計算就是對目標函數求解極小值的過程。因此對目標函數求偏導，并令其等于0，可得

φ（m）m=2GTWTdWd［G（m-m0）-（d-Gm0）］+2βWTmWm（m-m0）=0.（4）

基于式（4）求解，可以得到反演迭代公式為

m=m0+（GT+WTdWdG+βWTmWm）-1GTWTdWd（d-Gm0）.（5）

2.2 基于交叉梯度井中重磁數據聯合反演

2.2.1 交叉梯度函數定義

基于兩種不同物性參數的三維交叉梯度函數可定義為

t（x，y，z）=m1（x，y，z）×m2（x，y，z）.（6）

式中，為梯度運算符號;m1和m2代表著參與交叉梯度計算的不同參數，在井中重磁聯合反演中，m1和m2分別代表密度和磁化率。t值越小，說明兩種模型在該位置處的結構越為相似。t在x、y、z三個方向上的表達式分別如下：

tx=m1ym2z-m1zm2y ，ty=m1zm2x-m1xm2z ，tz=m1xm2y-m1ym2x .（7）

2.2.2 交叉梯度聯合反演目標函數

當加入交叉梯度約束項后，井中重磁聯合反演的目標函數轉化為

φ（m1，m2）=φ（m1）+φ（m2）+λφt（m1，m2），（8）

φt（m1，m2）=tTt=tTxtx+tTyty+tTztz.（9）

設B為交叉梯度t對m（m代表m1或m2）的偏導數矩陣，將t對m在m0處進行泰勒展開，忽略高階項，可得

t=t0+tm（m-m0）.（10）

記Bx=txm，By=tym，Bz=tzm，

則φtm=2BTB（m-m0）。

因此對目標函數求偏導，并令其等于0，可得

φ（m）m=2GTWTdWd［G（m-m0）-（d-Gm0）］+2βWTmWm（m-m0）+2λBTB（m-m0）=0.（11）

基于式（11）求解，可以得到聯合反演迭代公式為

m=m0+（GTWTdWdG+βWTmWm+λBTB）-1GTWTdWd（d-Gm0）.（12）

再使用共軛梯度法求解模型的更新量m（m代表密度和磁化率），再分別更新模型進行下一次反演迭代，直到得到滿足預設誤差精度的最優解。

3 理論模型反演

3.1 雙層礦體模型反演

設置了一個雙層疊置礦體模型，兩個礦體大小均為600 m×600 m×300 m，x和y方向范圍均是-300～300 m，上層礦體的深度范圍為500～800 m，下層礦體的深度范圍為1200～1500 m，模型如圖1所示。礦體的剩余密度值被設定為0.5 g/cm3，磁化率值被設定為0.01 SI，磁偏角和磁傾角為0°和90°。

井中測點從z=0起，相隔10 m設一個測點，正演計算井中測點的重磁異常。由于實際井中重磁測量主要觀測垂直分量，因此本文中反演主要基于垂直分量。圖2為井中重磁垂直分量異常曲線，圖中灰色部分顯示了礦體的深度范圍，可以看出井中重磁異常的垂直分量在礦體位置出現顯著變化，變化極值點可以確定礦體的上下邊界，表明了井中重磁測量具有很好的縱向分辨率，能夠很好地識別出縱向疊置礦體。相應地計算出地面重磁異常響應，圖3為地面重磁異常等值線圖，可以看到明顯的重磁異常平面分布，但無法識別縱向疊置的礦體。

根據正演得到的地面重磁異常數據和井中重磁異常數據，分別進行地面重磁數據獨立反演，井中重磁數據獨立反演及井中重磁數據交叉梯度聯合反演，反演結果如圖4所示。

從反演結果可以看出，對于縱向疊置的雙層礦體，僅利用地面重磁觀測數據進行反演，反演結果與實際礦體有很大差異，無法識別疊置的雙層礦體，同時反演的礦體物性與真實模型之間也存在較大的誤差。利用井中重磁數據進行反演，在縱向上可以明顯區分出兩個礦體，且礦體的位置與真實模型基本一致，與地面重磁數據反演結果相比，井中重磁數據反演結果的縱向分辨能力有了顯著的提升。反演所預測的磁化率最大值為0.009 SI，剩余密度最大值為0.31 g/cm3，與真實模型的物性也較為接近。綜合井中重磁數據進行聯合反演，在交叉梯度結構約束下，反演的礦體形態更加聚焦，礦體分布范圍更為準確，而且聯合反演結果的物性也比單獨井中重磁反演結果更加接近真實模型物性。

3.2 “W”型礦體模型反演

設置一“W”形狀礦體來模擬較為真實的地下礦脈構造。礦體的剩余密度值為0.5 g/cm3，磁化率值為0.01 SI，磁偏角和磁傾角為0°和90°。在研究區內設置5個觀測井，井位坐標分別為J1（-750，0）、J2（-350，0）、J3（0，0）、J4（350，0）、J5（750，0），井中觀測點自地面z=0 m起，相隔50 m到地下埋深z=2000 m處。模型及測井井位的分布如圖5所示。

圖6和圖7分別為“W”型礦體模型的井中重磁垂直分量異常曲線。可以看出，5口井的井中重磁垂直分量異常曲線在穿越礦體的埋深位置都出現了明顯的異常，根據井中重磁異常曲線可以大致了解礦體深度分布情況。

分別開展多井重磁異常單獨反演和基于交叉梯度約束的多井重磁異常聯合反演。圖8是多井重磁數據單獨反演結果，反演結果可以很好地顯示礦體的形態特征，但比真實模型范圍更加寬泛。圖9是多井重磁數據聯合反演結果，從圖中可以看出，基于交叉梯度約束后，聯合反演所得的結果更加聚焦，反演礦體的形狀輪廓與真實礦體模型基本一致，反演的物性值也更接近真實物性值，表明了聯合反演的有效性。

4 實際應用

4.1 研究區域背景

研究區為湖北省大冶市銅綠山銅鐵礦區，位于大冶市城區西南約3 km。銅綠山銅鐵礦是典型的矽卡巖型銅鐵礦，其次為銅礦石和鐵礦石。礦體的產出受銅綠山背斜大理巖殘留體與銅綠山石英二長閃長玢巖控制，礦體主要賦存于巖體與大理巖接觸帶［19］。

圖10為研究區地面剩余重磁異常分布等值線圖。研究區礦體的重磁數據多表現為高磁異常和局部較高重力異常，可以看出研究區內分布著多個磁異常與重力異常，特別是磁鐵礦體具有很強磁性，剩余磁異常基本反映了礦體的分布形態，而且礦區范圍內重力高異常和磁力高異常吻合性較好，表現出一定的重磁異常同源性。

目前在研究區僅有井位1鉆孔完成了井中磁測Bz分量異常的有效測量。鉆井1位于銅綠山外圍，井中磁測深度到達880 m。表1為井位1鉆孔的巖心資料統計表，表明在763～778 m深度處分布有含銅赤鐵礦磁鐵礦礦石，可以看出在礦體分布的位置，密度和磁化率都明顯偏高，具有很好的一致性，亦反映了該區域存在高磁高密度的重磁同源礦體，這為重磁聯合反演識別礦體提供了很好物性基礎。

圖11（a）為實測井中Bz分量磁異常曲線，它是由實際研究區的正常場值疊加礦體所產生的磁異常的綜合結果，需要對實測的磁異常數據進行正常場校正。當鉆井穿過磁性礦體時，鉆井井壁上產生的感應磁荷會對實測數據產生影響，需要對實測數據進行井壁磁荷校正。圖11（b）為校正以后的井中Bz分量磁異常。可以看出在760～780 m深度處顯示有一明顯磁異常，非常精確地指示了礦體的埋深位置。同時可以看到井位1處在地面剩余重磁異常圖（圖10）中均沒有表現出礦體重磁高異常特征，這主要是礦體厚度較薄，而且埋深較深，被淺部異常掩蓋，這也說明地面重磁測量對深部隱伏礦體無法有效識別，而井中測量可以對深部隱伏礦體具有很好的識別效果和探測能力。

4.2 聯合反演及解釋

首先利用地面重磁異常數據進行單獨反演，反演結果如圖12所示，然后利用地面重磁異常數據和井中磁測Bz分量數據進行交叉梯度聯合反演，結果如圖13所示。可以看出單獨反演的結果無法體現不同物性剖面上結構相似性，反演物性分布較為散亂，而聯合反演具有明顯的物性結構共同性。在重磁異常同源性較好的位置（圖13中實線黑框），聯合反演對反演結果起到了一定的約束作用，對應位置反演結果的相關性較好。在井中磁測曲線顯示異常的區域（750～800 m），聯合反演結果顯示出明顯的高磁化率，很好地反演出深部隱伏小型礦體分布（虛線黑框）。由于沒有井中重力實測數據參與聯合反演，在聯合反演的密度反演結果中沒能準確預測出深部隱伏礦體的分布，這進一步說明了地面重磁數據難以探測埋深較大的隱伏礦體，體現了井中重磁測量在深部礦體資源探測中的優勢。如果有多個井中重磁數據進行聯合反演，可以更為準確地刻畫深部礦體的空間形態。

5 結 論

（1）與地面重磁測量相比，井中重磁測量能近距離獲得深部礦體的重磁異常響應，具有很高的縱向分辨能力，可以很好地指示深部礦體的埋深范圍。

（2）利用井中重磁數據開展交叉梯度聯合反演，基于結構耦合約束，使反演的結果更加聚焦，反演的物性值更接近真實值，減少了反演的多解性，提高了反演的可靠性。

（3）井中重磁數據聯合反演獲得密度和磁化率分布，可以綜合多參數更好地對深部礦體進行評價，為深部礦產的開發提供技術支持。

參考文獻：

［1］ 胡文瑞，鮑敬偉.石油行業發展趨勢及中國對策研究［J］.中國石油大學學報（自然科學版），2018，42（4）：1-10.

HU Wenrui， BAO Jingwei. Development trends of oil industry and Chinas countermeasures［J］.Journal of China University of Petroleum（Edition of Natural Science），2018，42（4）：1-10.

［2］ NIND C， SEIGEL H， CHOUTEAU M， et al. Development of a borehole gravimeter for mining applications［J］.First Break，2007，25：71-77.

［3］ 付建偉.重力測井技術若干進展［J］.地球物理學進展，2010，25（2）：596-601.

FU Jianwei. Progress of gravity logging technology［J］. Progress in Geophysics，2010，25（2）：596-601.

［4］ 索奎，呂曉春，張貴賓，等.基于二維小波能量閾值的重磁數據去噪算法［J］.中國石油大學學報（自然科學版），2022，46（5）：36-45.

SUO Kui， L?Xiaochun， ZHANG Guibin， et al. Denoising algorithm for gravity and magnetic data based on two-dimensional wavelet energy threshold ［J］. Journal of China University of Petroleum（Edition of Natural Science）， 2022，46（5）：36-45.

［5］ 覃建華，楊琨，丁藝，等.基于KL-E的地質力學模型參數反演及應用［J］.西南石油大學學報（自然科學版），2022，44（2）：65-78.

QIN Jianhua， YANG Kun， DING Yi， et al. Inversion of geomechanical model parameters based on KL expansion and its application［J］. Journal of Southwest Petroleum University （Science & Technology Edition）， 2022，44（2）：65-78.

［6］ 王宗俊，田楠，范廷恩，等.復雜斷裂系統濁積儲層自相控反演技術研究［J］.西南石油大學學報（自然科學版），2023，45（2）：43-57.

WANG Zongjun， TIAN Nan， FAN Tingen， et al. Self-facies-control pre-stack inversion technology for turbidite sandstone reservoir with complex fault system［J］. Journal of Southwest Petroleum University （Science & Technology Edition）， 2023，45（2）：43-57.

［7］ HYOUNGREA R， LI Y G. Advantages of borehole vector gravity in density imaging［J］.Geophysics，2015，80（1）：G1-G13.

［8］ GENG M X， YANG Q J， HUANG D N.3D joint inversion of gravity-gradient and borehole gravity data［J］.Exploration Geophysics，2017，48（2）：151-165.

［9］ 李澤林，姚長利，鄭元滿.井地磁異常模量聯合反演［J］.地球物理學報，2018，61（12）：4942-4953.

LI Zelin，YAO Changli， ZHENG Yuanman. Joint inversion of surface and borehole magnetic amplitude data［J］.Chinese Journal of Geophysics，2018，61（12）：4942-4953.

［10］ 馮杰，歐洋，趙勇，等.三維井地磁測聯合約束反演［J］.地球物理學報，2019，62（10）：3686-3698.

FENG Jie， OU Yang， ZHAO Yong， et al. 3D joint constrained inversion of borehole and ground magnetic data［J］. Chinese Journal of Geophysics，2019，62（10）：3686-3698.

［11］ 徐凱軍，李猛，季春暉，等.基于模糊C均值聚類約束的井-地磁法聯合數據空間反演［J］.中國石油大學學報（自然科學版），2021，45（3）：55-64.

XU Kaijun， LI Meng， JI Chunhui， et al. Data space joint inversion of ground and borehole magnetic data with fuzzy C-means clustering constraints［J］.Journal of China University of Petroleum（Edition of Natural Science），2021，45（3）：55-64.

［12］ 魏寶君，谷靖雯，任臣，等.僅利用磁場幅度信息進行井間電磁反演成像的方法及可行性研究［J］.中國石油大學學報（自然科學版），2021，45（1）：50-59.

WEI Baojun， GU Jingwen， REN Chen， et al. Method and feasibility of cross-hole electromagnetic inversion and imaging by only using magnetic fields amplitude information［J］.Journal of China University of Petroleum（Edition of Natural Science），2021，45（1）：50-59.

［13］ 魏寶君.三維井間電磁場的聯合反演方法［J］.中國石油大學學報（自然科學版），2005，29（3）：20-26.

WEI Baojun. Combined inversion algorithm of three-dimensional cross-hole electromagnetic fields［J］.Journal of China University of Petroleum（Edition of Natural Science ），2005，29（3）：20-26.

［14］ 陳天勝，劉洋，魏修成.縱波和轉換波聯合AVO反演方法研究［J］.中國石油大學學報（自然科學版），2006，30（1）：33-37.

CHEN Tiansheng， LIU Yang， WEI Xiucheng. Joint amplitude versus offset inversion of P-P and P-SV seismic data［J］.Journal of China University of Petroleum （Edition of Natural Science ），2006，30（1）：33-37.

［15］ 王俊，孟小紅，陳召曦，等.交叉梯度理論及其在地球物理聯合反演中的應用［J］.地球物理學進展，2013，28（4）：2094-2103.

WANG Jun， MENG Xiaohong， CHEN Zhaoxi， et al. The theory of cross-gradient and its application in geophysical joint inversion［J］.Progress in Geophysics，2013，28（4）：2094-2103.

［16］ 彭國民，劉展.電磁和地震聯合反演研究現狀及發展趨勢［J］.石油地球物理勘探，2020，55（2）：465-474.

PENG Guomin， LIU Zhan. An overview of joint electromagnetic-seismic inversion and its future development［J］.Oil Geophysical Prospecting，2020，55（2）：465-474.

［17］ 丁文祥，袁峰，李曉暉，等.基于重磁聯合反演的寧蕪盆地鐘姑礦田深部地質結構解析及成礦預測［J］.地質學報，2018，92（11）：2301-2317.

DING Wenxiang， YUAN Feng， LI Xiaohui， et al.? Deep geological structure analysis and metallogenic prediction of Zhonggu ore field in the south section of Ningwu Basin based on gravity and magnetic joint inversion［J］. Acta Geologica Sinica，2018，92（11）：2301-2317.

［18］ CARRILLO J， PEREZ M A， GALLARDO L A， et al. Joint inversion of gravity and magnetic data using correspondence maps with application to geothermal fields［J］.Geophysical Journal International，2022，228（3）：1621-1636.

［19］ 吳飛，易露，尚世超，等.鄂東南礦集區銅綠山礦田深部探測物探方法技術及找礦效果探究［J］.資源環境與工程，2021，35（5）：606-610.

WU Fei， YI Lu， SHANG Shichao， et al. Deep exploration geophysical prospecting method technology and prospecting effect of Tonglüshan Orefield in Southeast Hubei Ore Concentration Area［J］.Resources Environment & Engineering，2021，35（5）：606-610.

（編輯 修榮榮）

基金項目：國家自然科學基金項目（42274181）

第一作者：徐凱軍（ 1979-），男，副教授，博士，研究方向為重磁電正反演及應用。E-mail：xkj79@upc.edu.cn。

引用格式：徐凱軍，季春暉，龐朝輝.面向深部礦產資源勘探的井中重磁三維交叉梯度聯合反演［J］.中國石油大學學報（自然科學版），2024，48（3）：65-74.

XU Kaijun， JI Chunhui， PANG Zhaohui. Cross-gradient joint inversion of 3D borehole gravity and magnetic data for deep mineral resources exploration［J］. Journal of China University of Petroleum （ Edition of Natural Science ），2024，48（3）：65-74.