孟凡興, 全旭東, 徐國蒼, 孔志召, 張 偉， 朱佳寧

(1.核工業(yè)航測遙感中心， 石家莊 050002；2.中核集團公司 鈾資源地球物理勘查技術中心(重點實驗室)， 石家莊 050002)

Occam反演在衛(wèi)境地區(qū)花崗巖型鈾礦勘查CSAMT資料處理中的應用

孟凡興1,2, 全旭東1,2, 徐國蒼1,2, 孔志召1,2, 張 偉1，2， 朱佳寧1,2

(1.核工業(yè)航測遙感中心， 石家莊 050002；2.中核集團公司 鈾資源地球物理勘查技術中心(重點實驗室)， 石家莊 050002)

花崗巖型鈾礦床一般產(chǎn)在巖體內部或巖體周邊，因此巖體發(fā)育情況的探測對鈾礦勘查工作的進一步開展具有重要意義。為了解Occam反演方法處理花崗巖型鈾礦勘查CSAMT數(shù)據(jù)時的應用效果，這里設計了隱伏巖體的理論模型，正演理論數(shù)據(jù)采用Occam反演方法進行反演。實驗表明，Occam反演對隱伏巖體的CSAMT理論數(shù)據(jù)具有穩(wěn)定收斂性，反演結果和理論模型能夠很好的吻合。結合衛(wèi)境地區(qū)的勘查實例驗證了Occam反演方法的有效性，為以后CSAMT方法間接尋找花崗巖型鈾礦提供了參考。

Occam反演； CSAMT方法； 隱伏巖體； 衛(wèi)境地區(qū)； 花崗巖型鈾礦

0 引言

可控源音頻大地電磁(CSAMT)法作為地球物理勘探方法中一種有效的電磁法勘探手段，已廣泛應用于各種深部金屬和非金屬礦產(chǎn)勘查[1-3]。總結不同地區(qū)的勘查經(jīng)驗發(fā)現(xiàn)，針對不同的勘查目的，選擇不同的反演方法處理，往往可以取得更好的效果[4-5]。

花崗巖型鈾礦床與花崗巖體有緊密空間關系和成因關系，它可產(chǎn)在巖體內部或巖體外圍不遠的一定范圍內[6-7]，因此巖體發(fā)育情況的探測對鈾礦勘查工作的進一步開展具有重要意義。

這里闡述了Occam反演的原理，并根據(jù)花崗巖型鈾礦勘查經(jīng)驗建立了隱伏巖體的正演模型，通過理論模型的反演以及在內蒙古北部衛(wèi)境地區(qū)的應用成果表明，采用Occam反演方法處理以探索隱伏巖體為目的的CSAMT數(shù)據(jù)，可以取得較好的效果。

1 Occam反演基本原理

Occam反演理論首先由Constable等[7-8]提出。地球物理反演因非唯一性，反演結果存在多個能滿足擬合數(shù)據(jù)的模型，Constable[8-9]認為反演為了獲得最優(yōu)解，反演的最終模型應盡可能的簡單、光滑，為了壓制來自非數(shù)據(jù)的模型構造，模型的粗糙度應盡可能的小。因此，Occam法是一種帶平滑約束的最小二乘法反演，這種方法反演是一種正則化的反演方法，它在尋找模型與原始數(shù)據(jù)最大擬合時，要求模型最光滑，因此受初始模型影響小，運算穩(wěn)定收斂，是一種有效的數(shù)據(jù)反演處理方法[10-11]。其反演的目標函數(shù)為：

U=R1+Q‖Rm‖2+

(1)

其中：μ-1為拉格朗日乘子；d為觀測數(shù)據(jù)向量；F為模型的響應；X*為期望擬合差水平；R=‖Rm‖2為模型粗糙度；W是m、μ駐點相關的一個泛函數(shù)。

F(m))]}

(2)

其中：m為模型參數(shù)；m0為初始模型參數(shù)；Cm為模型的協(xié)方差。

由于目標函數(shù)要同時尋求模型粗糙度和擬合差的最優(yōu)解，故在此引入拉格朗日乘子μ-1，產(chǎn)生一個無條件約束的目標泛函數(shù)U(m)來解最小化問題。

F(m)]-X*}

(3)

由于U(m)和Wμ對模型求導，駐點是相同的，則式(2)的駐點可由式(3)對一系列的μ求極小值，從而得到在擬合差允許范圍內模型粗糙度最小的迭代解[12-16]。

2 理論模型的正反演

為驗證Occam反演對花崗巖型鈾礦勘查地區(qū)CSAMT測深數(shù)據(jù)的反演效果，現(xiàn)結合衛(wèi)境地區(qū)實際工作中所采用的參數(shù)來進行理論模型的正反演。模型總深度600 m，背景電阻率為100 Ω·m，水平距離200 m～800 m的位置設計了埋深200 m、電阻率為1 000 Ω·m的高阻體對應隱伏的花崗巖體。

正演計算采用二維有限元方法，計算參數(shù)如下：計算測點距與野外實測點距50 m相同，收發(fā)距為6 km，正演網(wǎng)格為40×24。正演計算結果為1 Hz～8 192 Hz的28個頻率的TM模式的模型響應卡尼亞電阻率及阻抗相位數(shù)據(jù)。圖1為測線中心點處模型計算數(shù)據(jù)的展示。

反演模型是初始電阻率值為100 Ω·m的均勻半空間介質，圖2是理論模型和Occam反演結果圖。從反演結果與理論模型的對比中可以得出以下結論：高阻異常體反演結果的位置和理論模型設計的位置基本一致，但反演結果的高阻體有向上小幅度的延展，總體來看，反演結果的高阻體分布范圍與理論模型吻合較好，電阻率值也大致相同，低阻、高阻之間能夠漸變。說明采用Occam反演方法處理花崗巖型鈾礦勘查地區(qū)的CSAMT資料，可以大致查明隱伏花崗巖體的分布范圍，為鈾礦勘查工作的進一步開展提供參考。

圖1 500 m處模型計算數(shù)據(jù)Fig.1 Model calculation data in 500m(a)視電阻率曲線；(b)相位曲線

圖2 理論模型與反演效果對比Fig.2 The contrast of theoretical model and the inversion results(a)理論模型；(b)反演效果

3 應用實例

內蒙古衛(wèi)境地區(qū)經(jīng)過多年勘查，已發(fā)現(xiàn)數(shù)量較多的花崗巖型鈾礦化點及放射性異常點，目前工作的重點是圍繞這些礦化點和異常點進行深部找礦。采用Occam反演方法處理衛(wèi)境地區(qū)的CSAMT數(shù)據(jù)，取得了較好的效果。

3.1 研究區(qū)地質概況

研究區(qū)內出露地層主要是中新元古界青白口系艾勒格廟組(Qba)變質巖系，局部地區(qū)第四系覆蓋。

研究區(qū)內巖漿巖分布較廣，主要出露為早侏羅世灰白色中粗粒黑云母花崗巖。巖漿巖侵入于前寒武系地層中(圖3)，使本區(qū)巖性、構造復雜化，與此同時，形成了局部鈾的聚集，為鈾成礦提供了有利的條件[17]。

圖3 研究區(qū)地質及測線布置作圖Fig.3 Geological and line layout in the research area

3.2 研究區(qū)地球物理特征

該區(qū)巖石的物性特征見表1[4]，青白口系艾勒格廟組大理巖、結晶灰?guī)r電阻率平均值在400 Ω·m左右；早侏羅世中粗粒黑云母花崗巖電阻率平均值達1 000 Ω·m。依據(jù)巖石物性參數(shù)資料，總結得出艾勒格廟組變質巖(Qna)與早侏羅世花崗巖(J1γ)的電性特征，變質巖呈低阻特征，花崗巖呈高阻特征，二者接觸帶往往處于電阻率變化的梯度帶上。

3.3 野外工作參數(shù)

野外工作采用美國Zonge公司的GDP-32Ⅱ多功能電法儀，測量裝置見圖4，工作技術參數(shù)為：赤道偶極觀測裝置，標量測量；工作頻率為8 Hz～8 192 Hz，收發(fā)距不小于6 km，供電電流不小于8 A，勘探深度>600 m；視電阻率相對均方差±5 Ω·m，相位相對均方差在±20 mrad。

表1 衛(wèi)境地區(qū)巖石電阻率參數(shù)統(tǒng)計表Tab.1 The rock resistivity parameter statistics of Weijing area

圖4 CSAMT測量示意圖Fig.4 Instrumentation plan with CSAMT

3.4 應用效果分析

3.4.1 鉆孔驗證

根據(jù)勘查目的要求，在研究區(qū)布設剖面3條，其中L03線過鉆孔ZKG2。根據(jù)ZKG2鉆孔資料，淺部埋深80 m為艾勒格廟組變質巖，變質巖下部均為早侏羅世花崗巖。

圖5中清晰地反映出高阻電性層(反演電阻率一般大于500 Ω·m)分布在標高約1 020 m以下，推斷為早侏羅世花崗巖的反映；低阻電性層(反演電阻率多集中在200 Ω·m～400 Ω·m之間)分布在其上部，推斷為艾勒格廟組變質巖的反映。

圖5 實測剖面反演結果Fig.5 The inversion results of measured profile

鉆孔資料中花崗巖和變質巖接觸界線與圖2推斷解釋接觸界線存在一定的誤差(誤差范圍約50 m)，但整體上高阻電性層反映的下部早侏羅世花崗巖和低阻電性層反映的上部艾勒格廟組變質巖與鉆探揭露地質情況一致。

3.4.2 三維切片分析

圖6為海拔1 120 m和820 m反演電阻率平面分布及地質推斷解釋示意圖，從圖6中可以清晰地看出自淺部至深部(埋深約50 m～350 m)反演電阻率的分布特征及地質體的空間分布情況。

圖6 研究區(qū)反演電阻率平面切片示意圖Fig.6 The inversion of resistivity plane section sketch in the study area

海拔1 120 m反演電阻率平面圖中黑色虛線外部呈低阻電性特征，為變質巖(Qba)的反映，內部呈高阻電性特征，為花崗巖體(J1γ)的反映。海拔820 m反演電阻率平面圖整體表現(xiàn)為高阻電性特征，為花崗巖體(J1γ)的反映。

從淺部到深部，反演電阻率高阻體擴大，反映出南部的花崗巖厚度大、北部花崗巖厚度稍薄，總體上呈北薄、南厚的特征。

從圖7中可以看出：距離700 m橫切反演電阻率斷面呈上低下高的雙層電性結構，大面積花崗巖下伏在薄層變質巖之下；L02線上部低阻電性層、下部高阻電性層，高阻層為花崗巖體的反映， 低阻層為變質巖的反映，表明花崗巖體下伏在變質巖之下，變質巖厚度小于300 m；海拔820 m平面上高阻體均勻展布，為花崗巖的反映，也表明了研究區(qū)變質巖厚度不超過300 m。

圖7 研究區(qū)反演電阻率十字切片示意圖Fig.7 The inversion of resistivity cross section sketch in the study area

綜上所述：圖6和圖7清晰地顯示出衛(wèi)境地區(qū)反演電阻率高阻體和低阻體的空間分布情況，其變化特征直觀地反映出三維空間內早侏羅世花崗巖(J1γ)和艾勒格廟組變質巖(Qba)的空間展布情況。本次物探測量大致查明了隱伏花崗巖體的展布范圍，以及花崗巖和變質巖接觸帶的深部發(fā)育特征，為本地區(qū)鈾礦勘查工作的進一步開展提供了參考。

4 結論

Occam反演對初始模型的依賴程度較低，反演結果與理論模型可以很好的吻合，并且具有較高精確性和穩(wěn)定收斂特征。

通過衛(wèi)境地區(qū)勘查實例表明，采用Occam反演方法處理以探索隱伏巖體為目的的CSAMT數(shù)據(jù)，可以取得較好的效果，為以后花崗巖型鈾礦勘查區(qū)CSAMT資料的處理提供參考。

[1] 白國龍,何書躍,趙勇,等. CSAMT法在祁漫塔格隱盲礦勘查中的應用效果——以卡而卻卡礦區(qū)為例[J]. 物探化探計算技術,2015,37 (05):606-609. BAI G L ,HE SH Y ,ZHAO Y ,et al.The application effects of CSAMT in Qimantage blind ore exploration——Take Kaerqueka for example[J]. Computing techniques for Geophysical and Geochemical Exploration,2015,37(05):606-609.(In Chinese)

[2] 李寧生,馮志民,仵陽. CSAMT在固原市硝口巖鹽礦中的應用[J]. 物探與化探,2015(03):503-505. LI N SH ,FENG ZH M ,WU Y . The application of CSAMT to the Xiaokou rock salt deposit in Guyuan City[J]. Geophysical and Geochemical Exploralion,2015(03):503-505. (In Chinese)

[3] 孫博,李桐林,李鶴,等. 可控源音頻大地電磁測深法勘查深度研究[J]. 地球物理學進展, 2015,(02):836-839. SUN B ,LI T L ,LI H , et al. Study on the sounding of CSAMT[J]. Progress in Geophysics, 2015(02):836-839. (In Chinese)

[4] 孟凡興. 內蒙古四子王旗衛(wèi)境地區(qū)可控源音頻大地電磁測量報告[R].核工業(yè)航測遙感中心.2015. MENG F X . Measurement report of CSAMT in Wing jing area at Siziwang banner of Inner Mongolia[R]. Airborne Survey and Remote Sensing Center of Nuclear Industry, 2015.(In Chinese)

[5] 孟凡興,喬勇，胡鑒，等.可控源音頻大地電磁測量在花崗巖型鉑礦勘查中的應用[J].物探與化探，2016(3)：488-492+513. MENG F X,QIAO Y,HU J,et al.The applications of CSAMT to the exploration of granite uranium deposits[J].Geophysical and Geochemical Exploration,2016(03):488-492+513.(In Chinese)

[6] 余達淦，吳仁貴，陳培榮.鈾資源地質學[M] .哈爾濱:哈爾濱工程大學出版社.2005. YU D F ,WU R G ,CHEN P R , Uranium geology[M]. Harbin: engineering university press,2005.(In Chinese)

[7] 杜樂天. 中國熱液鈾礦成礦理論體系[J]. 鈾礦地質,2011(02):65-68+80. DU L T . On the theory system of hydrothermal uranium metalization in China[J]. Uranium Geology,2011(02):65-68+80. (In Chinese)

[8] CONSTABLE S C,PARKER R L, CONSTABE C G.Occam,s inversion:A practical algorithm for generating smooth modes from electromagnetic sounding data[J].Geophysics,1987,52(3):289-300.

[9] DEGROOT-HEDLIN C, CONSTABLE S.Occam's inversion to generate smooth,two-dimensional models from magnetotelluric data[J].Geo-physics,1990,55(12):1613-1624.

[10]何梅興. 可控源音頻大地電磁二維OCCAM反演研究[D].北京：中國地質大學,2009. HE M X . 2-D CSAMT OCCAM Inversion Research[D]. Beijing:China University of Geosciences,2009. (In Chinese)

[11]何梅興,胡祥云,葉益信,等. 2.5維可控源音頻大地電磁法Occam反演理論及應用[J]. 地球物理學進展,2011(06):2163-2170. HE M X ,HU X Y ,YE Y X , et al.2.5D controlled soutce audio-frequency magnetotellurics occam inversion[J]. Progress in Geophysics, 2011(06):2163-2170. (In Chinese)

[12]楊龍彬,張勝業(yè),劉文才. Occam反演在西準噶爾包古圖地區(qū)某銅礦AMT資料處理中的應用[J]. 地球物理學進展,2014(02):853-857. YANG L B ,ZHANG SH Y ,LIU W C . The application of Occam inversion in AMT data processing for a copper mine search in the west junggar Baogutu area[J]. Progress in Geophysics, 2014(02): 853-857. (In Chinese)

[13]王家映.地球物理反演理論[M].北京：中國地質大學學出版社.1998. WANG J Y . Inverse theory in geophysics[M]. Beijing: China University of Geosciences press, 1998.(In Chinese)

[14]何繼善.可控源音頻大地電磁法理論與應用[M].長沙：中南大學出版社，1990. HE J SH . Theory and application of controlled source audio-frequency magnetotelluric method[M].Changsha: Central south university press，1990.(In Chinese)

[15]王若,王妙月. 可控源音頻大地電磁數(shù)據(jù)的反演方法[J]. 地球物理學進展,2003(02):197-202. WANG R ,WANG M Y . Inversion method of controlled source audio-frequency magnetotelluric data[J]. Progress in Geophysics,2003(02):197-202.(In Chinese)

[16]吳小平,徐果明. 大地電磁數(shù)據(jù)的Occam反演改進[J]. 地球物理學報,1998(04):547-554. WU X P ,XU G M . Improvement of occam's inversion for mt data[J]. Acta geophysica Sinica, 1998(04):547-554.(In Chinese)

[17]胡鑒. 內蒙古四子王旗衛(wèi)境地區(qū)鈾資源潛力評價2013年度報告[R].核工業(yè)二〇八隊,2013. HU J . Uranium resources potential evaluation of 2013 annual report in Wing jing area at Siziwang banner of Inner Mongolia[R]. Geologic Party No.208,2013. (In Chinese)

Occam inversion in the application of the CSAMT data in granite type uranium deposit exploration in Weijing area

MENG Fanxing1,2， QUAN Xudong1,2， XU Guocang1,2，KONG Zhizhao1,2， ZHANG Wei1,2, ZHU Jianing1,2

(1.Airborne Survey and Remote Sensing Center of Nuclear Industry, Shijiazhuang 050002,China ;2.Key Laboratory for Geophysical Exploration Technology Center of Uranium Resource, Shijiazhuang 050002,China)

Granite type uranium deposit mainly locates in inner rock body or adjacent encasing body. Thus probing the developmental conditions of rock body have significance for uranium deposit exploration. In order to understand Occam inversion method processing granite type uranium deposit exploration and application effect of CSAMT data,The theory model of concealed rock mass is designed and he forward theory of Occam inversion method is adopted to improve the inversion data in this paper. Experiments show that Occam inversion of concealed rock mass possesses stable convergence theory of CSAMT data, the inversion result can match well to the theoretical model. Combined with guard position the exploration examples verify the effectiveness of Occam inversion method, for the CSAMT indirectly to find granite type uranium deposit provides the theoretical basis and practical experience.

Occam inversion; CSAMT; concealed rock mass; Weijing area; granite type uranium deposit

2016-02-22 改回日期：2016-03-31

中國核工業(yè)地質局基礎地質專項(201543)

孟凡興(1987-)，男，碩士，工程師，主要從事電磁法研究工作， E-mail：mfx911007@163.com。

1001-1749(2017)01-0032-06

P 631.4

A

10.3969/j.issn.1001-1749.2017.01.05