沈小慶楊炳南何 帥張德實

1.貴州省地質(zhì)礦產(chǎn)勘查開發(fā)局103地質(zhì)大隊,貴州 銅仁 554300;

2.自然資源部基巖區(qū)礦產(chǎn)資源勘查工程技術(shù)創(chuàng)新中心,貴州 貴陽550081;

3.中國地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,湖北 武漢430074

0 引言

貴州銅仁松桃錳礦整裝勘查區(qū)是中國重要的大型錳礦資源基地 (覃英等,2013;李嘯,2014)。近十年來,貴州省地礦局103地質(zhì)大隊通過長期產(chǎn)、學(xué)、研、用相結(jié)合,提出了一種新的錳礦床類型——氣液噴溢沉積型錳礦床,并建立了錳礦裂谷盆地噴溢沉積成礦理論(周琦和杜遠生,2019)。運用該成礦理論先后指導(dǎo)發(fā)現(xiàn)了整裝勘查區(qū)4個世界級超大型隱伏錳礦床和1個特大型富錳礦床(周琦等,2016a),實現(xiàn)了國內(nèi)錳礦找礦有史以來的最大突破(楊勝堂和涂飛飛,2017;李建威等,2019)。貴州松桃普覺錳礦位于錳礦整裝勘查區(qū)南部,區(qū)內(nèi)錳礦資源找礦潛力巨大(張遂等,2018)。多年來,在區(qū)內(nèi)開展了成礦地質(zhì)背景、成礦系統(tǒng)、找礦關(guān)鍵技術(shù)等大量研究工作(周琦等,2013),但對于礦體深部埋藏與特征、后期保存條件(后生構(gòu)造)的研究程度不高,研究方式局限,運用有效的物探方法在該區(qū)進行深部錳礦勘查,評價后期深部構(gòu)造與成礦地質(zhì)體顯得尤為重要。

高密度電法(楊玉蕊等,2012)、瞬變電磁(薛國強,2004;閆述等,2009)、激發(fā)極化法(李金銘,2004)等傳統(tǒng)物探方法僅適合中淺部勘查,高精度磁法在呈弱磁性的深部錳礦體勘探中收效甚微(張傳華,1985;熊光楚,1994)。隨著物探新理論、新技術(shù)、新方法的日趨成熟(底青云等,2019;米宏澤,2019),通過多年的應(yīng)用研究與實踐,具有較大探測深度的音頻大地電磁法已成為當(dāng)前應(yīng)用最廣泛的深部電磁勘探方法之一(何梅興等,2017;何帥等,2019)。音頻大地電磁法以天然的平面電磁波作為場源,通過觀測相互正交的電磁場分量來探測地下不同深度介質(zhì)的導(dǎo)電性結(jié)構(gòu) (Chave and Jones,2012;李忠等,2021)。近年來,針對錳礦整裝勘查區(qū)南華紀(jì)“大塘坡式”沉積型錳礦床(周琦等,2017;楊炳南等,2018),運用音頻大地電磁法識別礦區(qū)構(gòu)造格架、探測地下電性結(jié)構(gòu)、預(yù)測礦體展布形態(tài)及空間位置,建立了南華紀(jì)“大塘坡式”錳礦地質(zhì)-地球物理找礦模型,有效指導(dǎo)了深部找礦工作。

1 研究區(qū)概況

1.1 地質(zhì)概況

松桃普覺錳礦區(qū)位于上揚子?xùn)|緣成礦帶中南段,處于揚子地塊東南緣與華夏地塊的結(jié)合地帶,大地構(gòu)造位置位于揚子陸塊與南華活動帶過渡區(qū)的江南造山帶的西南段 (戴傳固等,2008,2010)。根據(jù)周琦等(2016b)通過對南華紀(jì)地層系統(tǒng)和地層分區(qū)的研究,研究區(qū)位于石阡-松桃-古丈Ⅲ級地塹盆地中的普覺Ⅳ級地塹盆地(圖1,圖2)。

研究區(qū)地層主要有青白口系紅子溪組(Pt3h),南華系中統(tǒng)兩界河組(Nh2l)、鐵絲坳組(Nh2t)、大塘坡組(Nh2d),南華系上統(tǒng)南沱組(Nh3n),震旦系至寒武系老堡組(Z-l)、震旦系下統(tǒng)陡山沱組(Z1d),寒武系第二統(tǒng)九門沖組(2jm)、變馬沖 組 (2b)、杷 榔 組 (2p)、清 虛 洞 組(2q),寒武系第三統(tǒng)高臺組(3g)、石冷水組(3s),寒武系第三統(tǒng)至芙蓉統(tǒng)婁山關(guān)組(3-4ls)及第四系(Q)。研究區(qū)位于梵凈山穹狀背斜北東側(cè),區(qū)內(nèi)以一系列北北東、北東向斷裂、褶曲為主的構(gòu)造格架(圖3)。區(qū)內(nèi)褶皺主要有猴子坳向斜,主要斷裂有紅石斷裂、三陽斷裂、揚立掌斷裂、木耳斷裂、水銀廠斷裂等。

1—控制Ⅲ級地塹盆地和地壘的同沉積斷層;2—控制Ⅳ級地塹盆地和地壘的同沉積斷層;3—Ⅳ級地塹盆地及所控制的錳礦床名稱;4—Ⅲ級地塹盆地范圍;5—Ⅲ級地壘范圍;6—研究區(qū)大地構(gòu)造所在位置;7—同沉積斷層編號;8—研究區(qū)在石阡-松桃-古丈(Ⅲ級)地塹盆地中的位置圖1 黔東及毗鄰區(qū)南華紀(jì)早期武陵次級裂谷盆地結(jié)構(gòu)與構(gòu)造古地理圖(周琦等,2016b)Fig.1 Structural and tectonic paleogeographic sketch of the Wuling secondary rift basin in the Early Nanhua period in eastern Guizhou and its adjacent areas (Zhou et al., 2016b)1-Synsedimentary faults that controlling 3rd-order graben basin and horst; 2-Synsedimentary faults that controlling 4th-order graben basin and horst; 3-Name of the Mn deposits controlled by 4th-order graben basin; 4-Range of 3rd-order graben basin; 5-Range of 3rdorder horst; 6-Geotectonic location of the study area; 7-Number of synsedimentary faults; 8-Location of the study area in the Shiqian-Songtao-Guzhang graben basin

1—兩界河期-鐵絲坳期砂礫巖沉積;2—大塘坡早期黑色炭質(zhì)頁巖;3—大塘坡中晚期粉砂質(zhì)頁巖;4—同沉積斷層;5—菱錳礦體圖2 黔渝湘毗鄰區(qū)南華紀(jì)兩界河-大塘坡期南華裂谷盆地復(fù)原(周琦等,2016b)Fig. 2 Rec onst r uct ion of the Nanhua rift basin in Liangjiehe-Datangpo stages of the Nanhua period in the Guizhou-Hunan-Chongqing border area, South China ( Zhou et al. , 2016b)1-Glutenite deposits in the Liangjiehe-Tiesiao stages;2 - Black carbonaceous shale in the early Datangpo stage; 3 - Silty shale in the middle and late datangpo stages; 4-Synsedimentary fault; 5-Rhodochrosite

區(qū)內(nèi)錳礦是一種新的錳礦床類型——氣液噴溢沉積型錳礦床,是典型的南華紀(jì)“大塘坡式”沉積型錳礦。錳礦體產(chǎn)于南華系下統(tǒng)大塘坡組第一段(Nh2d1)中,俗稱“含錳巖系”,其巖性主要由菱錳礦、炭質(zhì)頁巖、黏土巖組成。區(qū)內(nèi)錳礦為隱伏的錳礦床,礦體形態(tài)簡單,呈層狀、似層狀產(chǎn)出,具有品位高、厚度大、延伸展布廣等特點。錳礦體頂板為粉砂質(zhì)炭質(zhì)頁巖,局部頂板為薄層凝灰質(zhì)黏土巖。錳礦體底板為炭質(zhì)頁巖,少數(shù)礦體直接與鐵絲坳組含礫砂巖接觸,其底板為含礫砂巖。

1—走滑斷裂;2—逆斷層;3—正斷層;4—性質(zhì)不明斷層;5—AMT點位;6—研究區(qū)位置圖3 研究區(qū)主要斷裂格局及AMT測線點位布設(shè)圖Fig.3 Map showing the pattern of main faults and the layout of AMT surveyline points1-Strike-slip fault; 2-Reverse fault; 3-Normal fault; 4-Fault of unknown nature; 5-AMT point; 6-Location of the study area

1.2 物性特征

在系統(tǒng)收集研究區(qū)以往物性資料的基礎(chǔ)上,采用DZD-6A多功能電法儀補充測量了巖礦石野外露頭、鉆孔巖心物性數(shù)據(jù)。經(jīng)綜合分析,歸納了研究區(qū)物性特征(表1)。

表1 巖(礦)礦物性特征統(tǒng)計Table 1 Statistics of physical properties of rocks (minerals)

參考巖礦石電阻率特征,將巖性電阻率特征相近的相鄰地層進行合并。如寒武系第三統(tǒng)至芙蓉統(tǒng)婁山關(guān)組、寒武系第三統(tǒng)石冷水組巖性主要為電阻率特征呈高阻的白云巖,寒武系第二統(tǒng)清虛洞組巖性主要為電阻率特征呈高阻的灰?guī)r,將以上3個地層合并劃分為高阻電性層?？傮w上研究區(qū)地下電性結(jié)構(gòu)可劃分為“高阻—低阻—中高阻”3層模式(圖4)。

圖4 研究區(qū)電性結(jié)構(gòu)特征圖Fig.4 Diagram showing the resistivity structure of the study area

2 AMT野外數(shù)據(jù)采集與處理

此次AMT研究工作采用V8多功能電法工作站開展野外數(shù)據(jù)采集,投入1臺主機,型號為V8-6R,1臺輔助接收機,型號為RXU-3ER。采用“十”字交叉型布極方式進行四分量張量觀測,采集頻率范圍0.35～10400 Hz,采集參數(shù)為磁道數(shù)據(jù)和電道數(shù)據(jù)?；谘芯拷?jīng)費考慮,完成1條剖面,測線總長11 km,點距50～200 m,共計100個測點的數(shù)據(jù)采集工作。AMT剖面南東段地形起伏較大、植被茂密。北西段地勢相對平緩,觀測條件較好。測區(qū)總體電磁干擾較弱,主要干擾源為零散分布于測區(qū)的村寨民用電線。AMT剖面方位角147°,大致垂直于普覺Ⅳ級地塹盆地空間展布方向,與區(qū)內(nèi)燕山期北東向構(gòu)造垂直,與北北東向構(gòu)造斜交。通過開展不同觀測時間試驗(圖5)、不同電極距長度試驗(圖6)等系列采集參數(shù)試驗,優(yōu)選出了適合研究區(qū)的野外采集技術(shù)參數(shù),即采集參數(shù)增益設(shè)置為1,低通濾波選擇0檔,野外采集電極距長度為40 m,單點采集時間大于60 min。

a—觀測時間30 min;b—觀測時間45 min;c—觀測時間60 min圖5 不同觀測時間測點曲線對比Fig.5 Comparison of measuring point curves at different observation times(a) A 30-minute observation; (b) A 45-minute observation; (c) A 60-minute observation

野外采集的AMT原始資料記錄的是測點電場和磁場分量的時間序列,通過數(shù)據(jù)處理獲取測點的視電阻率和相位等AMT參數(shù)的響應(yīng)。首先根據(jù)測點班報記錄更改測點參數(shù)信息,采用Fourier變換法進行AMT的頻譜分析,然后使用Robust處理技術(shù)計算大地電磁的阻抗張量,獲取用于編輯的互功率譜。通過人機交互的方式對子功率譜進行篩選,剔除嚴重受干擾的功率譜,保留質(zhì)量高的功率譜參與統(tǒng)計計算。預(yù)處理后的AMT資料導(dǎo)入MTSoft2D軟件作編輯平滑、靜態(tài)校正(Jiracek,1990)、地形校正(韓騎等,2015)等進一步處理。圖7a—7c為160～164號測點經(jīng)簡單編輯平滑后的電阻率曲線。由圖可見,各曲線連續(xù)性好,形態(tài)明確,表明測點數(shù)據(jù)質(zhì)量較好。162號點由于受靜態(tài)效應(yīng)影響,xy曲線和yx曲線分離,需作靜態(tài)校正處理(圖7b)。圖7d—7f為通過觀察相鄰測點電阻率數(shù)據(jù)曲線特征,采用曲線平移法實現(xiàn)162號點靜態(tài)校正后的曲線(圖7e)。

藍色—xy曲線;紅色—yx曲線圖7 162號測點靜校正前后曲線對比Fig.7 Comparison of curves before and after static correction at measuring point 162Blue-the xy curve;Red-the yx curve

a—電極距為20 m;b—電極距為30 m;c—電極距為40 m圖6 不同電極距測點曲線對比Fig.6 Comparison of measuring point curves with different electrode distances(a) A 20-meter electrode spacing; (b) A 30-meter electrode spacing; (c) A 40-meter electrode spacing

3 理論模型反演計算

3.1 理論模型建立

根據(jù)研究區(qū)地電結(jié)構(gòu)特征,在不考慮地形因素的情況下,建立“高—低—高”3層電性結(jié)構(gòu)的簡化構(gòu)造模型(圖8)。模型長度為4 km,深度為3 km,在2.5 km處發(fā)育正斷層,斷層切穿整個模型,斷距為0.5 km。模型上部表示為寒武系第二統(tǒng)至芙蓉統(tǒng)碳酸鹽巖高阻層,其電阻率填充為1500 Ω·m;模型中部表示為寒武系第二統(tǒng)與南華系碎屑巖低阻層,忽略相對較薄的震旦系下統(tǒng)陡山沱組白云巖高阻夾層,電阻率填充為500 Ω·m,在其底部設(shè)置了以炭質(zhì)頁巖和錳礦體為主的含錳巖系低阻層,電阻率填充為50 Ω·m;模型深部表示為青白口系含礫砂巖、淺變質(zhì)巖為主的高阻層,深部電阻率填充為2500 Ω·m;正斷層電阻率填充為100 Ω·m。使用Mtsoft2D正演模塊采用50×50的矩形網(wǎng)格剖分單元對該模型進行二維有限元正演(劉向沖,2019)。按照對數(shù)等間距在0.35～10400 Hz范圍內(nèi)取與實測頻率相同的60個測量頻點,模擬計算40個測點TE模式和TM模式的視電阻率和阻抗相位。

圖8 正斷層地電模型Fig.8 Geoelectric model of the normal fault

圖9不同初始模型及參數(shù)的NLCG反演成果圖Fig.9 NLCG inversion results of different initial models and inversion methods(a)Initial model Occam-TE data+TE inversion;(b) Initial model Occam-TE data+TM inversion; (c) Initial model Occam-TE data+TE&TM inversion; (d) Initial model Bostick-TE data+TE inversion; (e) Initial model Bostick-TM data+TE inversion; (f) Initial model Bostick-(TE+TM)data+TE inversion

a—初始模型Occam-TE數(shù)據(jù)+TE反演;b—初始模型Occam-TE數(shù)據(jù)+TM反演;c—初始模型Occam-TE數(shù)據(jù)+TE&TM反演;d—初始模型Bostick-TE數(shù)據(jù)+TE反演;e—初始模型Bostick-TM數(shù)據(jù)+TE反演;f—初始模型Bostick-(TE+TM)數(shù)據(jù)+TE反演;

3.2 理論模型反演計算

目前,針對AMT資料較為成熟的反演方法有Bostick反演法(周虬,1985)、快速松弛反演法(Smith and Booker,1991)、二維Occam反演法(De Groot-Hedlin and Constable,1990)、非線性共軛梯度反演法 (Rodi and Mackie,2001)。Bostick反演法是一種擬二維的近似反演方法,可以快速直觀地顯示視電阻率隨深度的變化情況,但具有反演結(jié)果精度不高的缺點?？焖偎沙诜囱莘ㄊ峭ㄟ^減少正演的次數(shù),從而減少反演運算時間,該方法的反演結(jié)果具有多解性。二維Occam反演法相對于快速松弛反演法具有穩(wěn)定的收斂性,但其反演計算時間較長。非線性共軛梯度反演法(NLCG)相對于其他反演方法具有計算結(jié)果穩(wěn)定、效率高等優(yōu)點,但其對初始模型的依賴性較強。通過對比各反演方法的優(yōu)缺點(康敏等,2017),筆者分別選用Bostick和Occam一維反演結(jié)果作為初始模型,采用非線性共軛梯度反演法優(yōu)選適合于理論模型的二維反演模式,將該反演模式運用于研究區(qū)實測數(shù)據(jù)的反演。

理論模型的6種不同模式反演結(jié)果如圖9所示,以O(shè)ccam-TE模式一維反演為初始模型,二維反演采用TE模式、TM模式、TE和TM聯(lián)合模式的反演結(jié)果如圖9a—9c。以Bostick-TE模式一維反演為初始模型,二維反演采用TE模式、TM模式、TE和TM聯(lián)合模式的反演結(jié)果如圖9d—9f。NLCG 6種不同模式的反演結(jié)果在縱向上均呈現(xiàn)了“高—低—高”的3層電性結(jié)構(gòu),各層厚度與理論模型相當(dāng)。但對于橫向電性特征的辨識以及構(gòu)造的識別表現(xiàn)不一。在圖9a—9c中,通過上部2.5 km附近和中部低阻層的電阻率不連續(xù)性,較好地識別了斷層性質(zhì)及其發(fā)育特征,對由于斷層錯動引起的地層橫向變化有較高的靈敏度。特別是在圖9b模式下的反演結(jié)果與理論模型最為吻合。在圖9d—9f中,中部低阻層在橫向上從左至右具有抬升趨勢,這是模型橫向電性變化的綜合反映,對構(gòu)造識別不敏感。另外,理論模型反演結(jié)果可見,電性特征呈低阻的含錳巖系在反演斷面中并非位于低阻層中部電阻率值最低的位置,而是位于低阻層底部,電阻率由低向高漸變的過渡帶。這為實測剖面推測含錳巖系埋深及其空間展布提供了理論支撐。

通過反演結(jié)果對比分析,認為采用圖9b中的模式,即以O(shè)ccam-TE一維反演為初始模型的二維NLCG-TM反演模式更能客觀反映理論模型電性結(jié)構(gòu)特征,有利于斷裂構(gòu)造的識別。實測AMT數(shù)據(jù)經(jīng)處理后,采用上述反演模式進行二維帶地形反演。

4 實測AMT剖面電性結(jié)構(gòu)分析

根據(jù)巖性資料、物性資料以及電阻率等值線連續(xù)的原則,鑒于地層巖性電性差異和AMT的分辨率,參考剖面地質(zhì)模型(圖10),對反演圖(圖11)相關(guān)巖性相近的地層層位進行了電性層合并,劃分出了a、b、c、d共4個電性區(qū)塊。其中,a為高阻層,反映了寒武系第三統(tǒng)至芙蓉統(tǒng)婁山關(guān)組至寒武系第二統(tǒng)清虛洞組碳酸鹽巖地層;b為低阻層,反映了以頁巖、炭質(zhì)頁巖等碎屑巖為主的寒武系第二統(tǒng)至南華系中上統(tǒng)地層;c為高阻層,是青白口系變余砂巖、板巖電性特征的綜合反映;d為由b至c的過渡層,即是低阻向高阻過渡的相對中高阻層。通過理論模型反演分析,含錳巖系位于斷面低阻體下部,處于由低阻向高阻漸變的過渡帶。結(jié)合斷面低阻體形態(tài),推測了含錳地質(zhì)體大致空間形態(tài),其埋深由北西至南東逐漸加大,整體呈單斜狀。

1—紅子溪組;2—清水江組;3—兩界河組;4—鐵絲坳組;5—大塘坡組第二段;6—南沱組;7—震旦系;8—牛蹄塘組;9—九門沖組;10—變馬沖組;11—杷榔組;12—清虛洞組;13—寒武系第三統(tǒng)至芙蓉統(tǒng);14—板巖;15—變余砂巖;16—推測含錳巖系;17—粉砂質(zhì)頁巖;18—含礫砂巖;19—砂巖;20—炭質(zhì)頁巖;21—灰?guī)r;22—白云巖;23—第四系松散堆積體;24—推測斷層圖10 900號剖面地質(zhì)模型圖Fig.10 Geological model of the profile 9001-The Hongzixi formation; 2-The Qingshuijiang formation; 3-The Liangjiehe formation; 4-The Tiesiao formation; 5-The second section of the Datangpo formation; 6-The Nantuo formation; 7-The Sinian system; 8-The Niutitang formation; 9-The Jiumenchong formation; 10-The Bianmachong formation; 11-The Palang formation; 12-The Qingxudong formation; 13-Cambrian third system to Furong system; 14-Slate; 15-Blastopsammite; 16-Supposed manganese-bearing rock series; 17-Sand shale; 18-Pebbled sandstone;19-Sandstone; 20-Coaly shale; 21-Limestone; 22-Dolomite; 23-Quaternary loose deposits; 24-Supposed fault

斷面淺部156號點附近兩側(cè)電阻率差異明顯,表現(xiàn)為高阻與低阻的巖性分界面,推測為研究區(qū)主要的控礦斷層F0,F0斷層延伸較深,深部未對含錳巖系造成破壞。180至210號點低阻層電阻率等值線不連續(xù),推測是由已知斷層F1引起。240至280號點之間,電阻率等值線橫向不連續(xù),具有明顯的低阻上凸特征,推測是由F2、F3斷層引起。300號點下方橫向的電阻率等值線發(fā)生了較明顯的低阻上凸,在其深部伴有電性分界面,推測由F4斷層與隱伏的F5斷層組合成的“鍥形”構(gòu)造。

由于剖面地層產(chǎn)狀較陡,構(gòu)造較發(fā)育,縱向電性結(jié)構(gòu)在不同區(qū)段有所不同。在112～156號點范圍,為單一的高阻結(jié)構(gòu);在156～210號點范圍,電性結(jié)構(gòu)為“低阻—相對高阻”漸變的雙層模式;在210～230號點范圍,電性結(jié)構(gòu)為“高阻—低阻—相對高阻”的3層模式;230～310號點范圍,電性結(jié)構(gòu)為“高阻—低阻”的雙層模式。

a—寒武系婁山關(guān)組和清虛洞組巖性組合層;b—寒武系第二統(tǒng)至南華系中上統(tǒng)巖性組合層;c—青白口系地層;d—南華系與青白口系過渡帶巖性組合層;1—推測地層界線;2—推測斷層及編號;3—推測含錳地質(zhì)體圖11 900號剖面NLCG反演成果圖Fig.11 Results of the profile 900 from the NLCG inversion(a) Lithologic combination layers of the Cambrian Loushanguan formation and the Qingxudong formation; (b) Lithologic combination layers from the second Cambrian to the middle-upper Nanhua system; (c) The Qingbaikou formation; (d) Lithologic assemblage layer in the transitional zone between the Nanhua system and the Qingbaikou system 1-Supposed stratigraphic boundaries; 2-Supposed faults and numbers; 3-Supposed manganese-bearing geological bodies

5 結(jié)論

(1)通過巖礦石電性特征分析,建立了研究區(qū)總體“高阻—低阻—中高阻”3層的電性結(jié)構(gòu)分層模式。實測AMT剖面地表不同位置出露地層巖性和電性差異較大,地層產(chǎn)狀較陡,剖面電性結(jié)構(gòu)為單層、雙層、三層的多元組合類型。

(2)建立與研究區(qū)地電結(jié)構(gòu)近似的構(gòu)造理論模型,采用6種不同反演模式對理論模型進行反演計算,通過反演結(jié)果對比分析,優(yōu)選出了適合研究區(qū)的最佳反演參數(shù)組合:初始模型Occam-TE數(shù)據(jù)+NLCG-TM反演,即選用Occam反演法對TE數(shù)據(jù)做一維反演,以一維反演結(jié)果為初始模型對TM數(shù)據(jù)做非線性共軛梯度反演。

(3)AMT剖面反演成果清晰地反映了地層巖性空間展布和構(gòu)造格架。通過理論模型正反演分析,含錳巖系位于斷面低阻層底部,結(jié)合實測數(shù)據(jù)反演剖面低阻層的展布特征,間接地識別了含錳地質(zhì)體的空間展布形態(tài),查明了后期構(gòu)造F0等斷層與含錳地質(zhì)體的接觸關(guān)系,為分析研究區(qū)錳礦成礦地質(zhì)條件和深部找礦預(yù)測提供了地球物理依據(jù)。

(4)研究表明,音頻大地電磁法能夠有效的識別研究區(qū)地電結(jié)構(gòu)及深部構(gòu)造展布特征。通過理論模型的反演研究,優(yōu)選最佳反演方案,能夠更加客觀地對實測數(shù)據(jù)進行反演和解釋,為區(qū)內(nèi)錳礦深部找礦提供地球物理支撐。