徐 彬 陳建平 相 軒

(1.中國地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京100083;2.北京市國土資源信息研究開發(fā)重點實驗室，北京100083)

目前，山東焦家金成礦帶內(nèi)已發(fā)現(xiàn)了新城金礦、焦家金礦、河西金礦等3 座特大型金礦床和河?xùn)|金礦、靈山溝金礦、上莊金礦、東季金礦、望兒山金礦、馬塘金礦等6 座大中型金礦床［1］。但隨著地質(zhì)勘查程度的不斷提高，區(qū)內(nèi)的平均勘查深度已接近－500 m，預(yù)示著地表、淺部礦產(chǎn)資源越來越少，金礦找尋難度越來越大，加之礦產(chǎn)資源開發(fā)利用程度的不斷提高，使得許多老礦山已經(jīng)出現(xiàn)了資源危機［2-9］。因此，為了實現(xiàn)礦山資源接替，拓展找礦空間，攻深找盲已成為該地區(qū)金礦的重點找礦方向。

1 研究區(qū)地質(zhì)概況

研究區(qū)大地構(gòu)造位置處于華北板塊( Ⅰ) 、膠北地塊( Ⅱ) 、膠北隆起及凹陷區(qū)( Ⅲ) 、膠北隆起( Ⅳ)之膠北凸起( Ⅴ) 西緣，我國著名的膠西北金成礦區(qū)內(nèi)［10］。區(qū)內(nèi)地層比較簡單，主要出露荊山群祿格莊組( Pt1jL) 、疙瘩狀石榴矽線黑云片巖和第四系( Q)含礫黏土質(zhì)砂土、輕亞黏土、礫石和含礫石的混粒級砂以及粉砂沉積物。區(qū)內(nèi)斷裂構(gòu)造發(fā)育，主要發(fā)育脆性斷裂構(gòu)造，以焦家斷裂構(gòu)造為主干，同時包括其上下盤伴生以及派生的低序次斷裂，其中，NE 向的斷裂基本控制了該區(qū)帶內(nèi)金礦床的分布。區(qū)內(nèi)巖漿巖廣布，形成時代自中太古代至中生代均有出露，其中以中生代燕山早期( 晚侏羅世) 玲瓏序列為主體，中生代以來的侵入巖受近EW，NE，NNE 向斷裂制約。區(qū)內(nèi)已發(fā)現(xiàn)評價的金礦床( 點) 共計60 余處，已發(fā)現(xiàn)的金礦床均受控于NE 向的焦家主干斷裂及下盤伴生、派生的次級斷裂及裂隙群［11-13］。

2 三維成礦條件分析與成礦預(yù)測

三維成礦預(yù)測流程見圖1。

2.1 找礦地質(zhì)模型

圖1 三維成礦預(yù)測流程Fig.1 Three-dimensional metallogenic prediction process

在分析焦家金成礦帶大中型典型礦床地質(zhì)特征和控礦要素的基礎(chǔ)上，建立了區(qū)內(nèi)地質(zhì)找礦模型，見表1。

表1 研究區(qū)金礦床地質(zhì)找礦模型Table 1 Geological prospecting model of gold deposit in research area

2.2 三維地質(zhì)建模及成礦條件

采用Surpac［14-16］軟件對區(qū)內(nèi)的巖體、已知礦體、物化探異常等進行三維實體建模，區(qū)內(nèi)巖體的三維實體模型見圖2。根據(jù)現(xiàn)有地質(zhì)資料，特別是勘探線的網(wǎng)度、礦體的大小、礦體邊界的復(fù)雜度、采礦設(shè)計的要求及實際情況，確定單元塊體大小為100 m ×100 m×50 m，共劃分單元塊1 992 672 個，其中1 193 148個單元塊位于研究區(qū)內(nèi)部。采用立方體預(yù)測模型法對各個變量進行成礦有利條件的分析與提取，在此基礎(chǔ)上采用三維信息量法對各預(yù)測變量進行評價。

2.2.1 賦礦巖體信息

巖漿巖為熱液礦床成礦提供熱源及成礦物質(zhì)，對于隱伏巖體來說，礦體一般發(fā)現(xiàn)于巖體周邊及表面一定區(qū)域內(nèi)，因此，巖體緩沖區(qū)是不可或缺的預(yù)測變量。分別使用巖體實體模型和已知礦體實體模型對立方體塊模型進行約束，將不同巖體所包含的塊體單元和不同巖體所包含的礦體( 塊) 單元作為礦床預(yù)測中的巖性變量和先驗條件。經(jīng)統(tǒng)計，區(qū)內(nèi)已知礦塊數(shù)量為4 036 個，其中大部分礦塊位于玲瓏序列和馬連莊序列內(nèi)，可見金礦床在玲瓏序列和馬連莊序列內(nèi)含礦性較好。

2.2.2 蝕變帶信息

研究區(qū)內(nèi)金礦床均賦存于斷裂帶所控制的蝕變巖帶內(nèi)，蝕變巖帶的發(fā)育空間、形狀產(chǎn)狀與對應(yīng)的斷裂帶基本一致。因此，可以通過統(tǒng)計已知礦體( 塊)中各蝕變帶含量來確定各蝕變帶對成礦的影響。經(jīng)統(tǒng)計，區(qū)內(nèi)89.40% 的金礦床產(chǎn)于蝕變帶內(nèi)，其中41.38%的金礦床位于絹英巖化花崗質(zhì)碎裂巖帶內(nèi)，26.88%的金礦床位于黃鐵絹英巖化碎裂巖帶內(nèi)，16.06%的金礦床位于絹英巖化花崗巖內(nèi)，5.08%的金礦床位于變輝長質(zhì)碎裂巖帶內(nèi)?？梢姡佊r化花崗質(zhì)碎裂巖和黃鐵絹英巖化碎裂巖是區(qū)內(nèi)礦體產(chǎn)出的主要蝕變帶類型。

2.2.3 構(gòu)造信息

(1) 構(gòu)造帶特征。強烈的構(gòu)造活動區(qū)域是成礦流體運移的通道，而礦質(zhì)沉積需要一個相對平靜的環(huán)境，在主干斷裂的旁側(cè)，構(gòu)造發(fā)育相對較弱，該類區(qū)域通常為礦體形成的相對有利區(qū)。因此，結(jié)合研究區(qū)實際情況，取200 m 對斷裂進行緩沖區(qū)分析，采用斷裂及其緩沖區(qū)實體模型對三維立方體模型進行限定，劃分出斷裂及其緩沖區(qū)所包含的單元塊體。經(jīng)統(tǒng)計，有87.07%的礦體( 塊) 落在該緩沖區(qū)內(nèi)，證明緩沖區(qū)是一個非常重要的預(yù)測變量。

(2) 構(gòu)造展布特征。采用構(gòu)造等密度( 反映線性構(gòu)造發(fā)育程度) 與構(gòu)造頻數(shù)( 反映線性構(gòu)造復(fù)雜程度) 的比值來定量分析區(qū)域主干斷裂的特征［17］。以焦家金礦帶為背景，NE 向斷裂( 其中，焦家斷裂、望兒山斷裂和河西斷裂為區(qū)內(nèi)的重要構(gòu)造) 為焦家金礦帶區(qū)域內(nèi)的主干斷裂，計算研究區(qū)內(nèi)斷裂等密度和斷裂頻數(shù)，并計算兩者之比。經(jīng)統(tǒng)計，選取主干斷裂值為0.009 17 ～0.011 67，作為成礦有利區(qū)間。

(3) 構(gòu)造交匯特征。斷裂交點數(shù)指單元面積內(nèi)斷裂交點的總數(shù)，其高值區(qū)代表多組構(gòu)造交匯的部位，是巖漿熱液上移最重要的特征，其旁側(cè)為成礦的有利空間。計算每個單元塊體內(nèi)的斷裂交點數(shù)，經(jīng)統(tǒng)計，選取構(gòu)造交匯點值為0.005 ～1.250，作為成礦有利區(qū)間。

(4) 構(gòu)造巖漿活動特征。巖漿熱液以斷裂為通道上涌，形成了礦體賦存的模式。構(gòu)造中心對稱度與構(gòu)造的疊加分析可反映構(gòu)造巖漿活動的特征，構(gòu)造中心對稱度與巖體的疊加分析可反映巖體的分異中心特征。經(jīng)統(tǒng)計，選取中心對稱度值為0.001 ～0.020，作為成礦有利因子。

2.2.4 等間距控礦找礦信息

構(gòu)造帶和構(gòu)造形跡空間展布的韻律定向性和間距的倍數(shù)性是構(gòu)造距離的2 個特點，研究礦床的等距性展布特征對指導(dǎo)盲區(qū)找礦具有指示性作用。研究區(qū)內(nèi)同一控礦構(gòu)造帶內(nèi)的礦床( 點) 具有等距性分布規(guī)律，該等距性不僅在走向上，而且存在于傾向上，據(jù)此可認(rèn)為區(qū)內(nèi)深部可能存在第二梯段盲礦體［18-19］。焦家、河西以及新城等礦床沿斷裂帶具有一定的等距性，礦體往往賦存于斷裂構(gòu)造的交叉、復(fù)合部位或構(gòu)造發(fā)生突然轉(zhuǎn)折的部位，由此可確定沿焦家主干斷裂方向的間距為1 500 ～1 800 m。焦家斷裂及其低序次構(gòu)造控制礦體的產(chǎn)狀具有走向NE，傾向NW，傾角約45°，在剖面上普遍具有向SW 側(cè)伏的規(guī)律，據(jù)此，可沿低級序斷裂方向確定間距為1 200 ～1 500 m。參考距焦家金礦約3 km，具有相同成因類型及控礦條件的新城金礦和望兒山金礦的成礦深度，將第一控礦空間延深至－1 100 m，第二控礦空間延深至－1 500 m。蝕變圍巖平行主斷裂以斷層泥為中心成對稱帶狀分布，蝕變帶一般寬70 ～250 m，最大370 m，平均150 m，主要圍巖蝕變類型有綠泥石化、鉀長石化、紅化，據(jù)此可確定等間距帶帶寬為250 m。

基于以上分析，可建立一個理想的等間距控礦菱形模型，將該模型與焦家礦區(qū)疊加，結(jié)合焦家實際地質(zhì)情況得出焦家金成礦帶等間距控礦模型，經(jīng)統(tǒng)計，研究區(qū)內(nèi)93.91%的礦體落在該模型中。

2.2.5 地球物理信息

研究區(qū)東部焦家斷裂帶視電阻率等值線呈舒緩波狀向NW 緩傾，反映出焦家斷裂帶向深部平緩延伸。電阻率等值線的緩陡反映了斷裂相對平緩和陡傾的部位，其陡變帶和拐點一般是找礦有利層位。對41 條可控源大地音頻電磁剖面進行觀察統(tǒng)計分析，發(fā)現(xiàn)視電阻率值為0，1 000，2 000，26 000，39 000 Ω·m的等值線陡變帶和拐點穿過的已知礦體較多，因此可提取等值線的陡變帶以及拐點作為找礦的指示變量。根據(jù)可控源大地音頻電磁剖面上的視電祖率等值線，建立視電阻率等值線曲面模型，將其與區(qū)內(nèi)已知礦體模型進行疊加分析，選取已知礦體賦存較好的2 層( 視電阻率值為1 000，20 000 Ω·m) ，在立方體模型中賦值，取位于異常值區(qū)間內(nèi)的礦體提取地球物理異常信息。經(jīng)統(tǒng)計，可控源大地音頻電磁剖面數(shù)據(jù)覆蓋范圍內(nèi)59.7%的礦體落在地球物理異常區(qū)域內(nèi)，充分說明視電阻率值為1 000 ～20 000 Ω·m的等值線曲面為金礦床的主要賦存層位，夾持著構(gòu)造蝕變帶，是區(qū)內(nèi)成礦的有利地球物理信息。

2.2.6 地球化學(xué)信息

對鉆孔樣品的各元素含量進行克里金插值［20］，依據(jù)井中元素取樣化驗分析結(jié)果，建立了Au、Ag、As、Bi、Cu、Hg、Pb、Sb、Zn 等9 類巖石樣品微量元素異常塊體模型，并依據(jù)元素在縱向上的分帶規(guī)律，確定了以Hg、Ag、Sb、Pb 為礦體頭部異常元素，以As、Au、Zn 為近礦異常元素，以Bi、Cu 為礦體尾部異常元素。經(jīng)統(tǒng)計，鉆孔數(shù)據(jù)覆蓋范圍內(nèi)金礦床78.49%落在Au 異常內(nèi)，井中元素剖面覆蓋范圍內(nèi)37.77%的礦體落在上述各元素的組合異常內(nèi)。

2.3 成礦預(yù)測及靶區(qū)圈定

將賦礦巖體( 玲瓏序列、馬連莊序列) 、蝕變帶( 黃鐵絹英巖化碎裂巖、絹英巖化花崗質(zhì)碎裂巖、變輝長質(zhì)碎裂巖) 、斷裂緩沖區(qū)、構(gòu)造交點數(shù)、構(gòu)造中心對稱度、主干斷裂、Au 品位值、元素組合異常、等間距控礦、視電阻率異常等13 個預(yù)測變量的成礦有利區(qū)間作為三維信息法［21］的計算因子，計算各地質(zhì)因素、找礦標(biāo)志所提供的找礦信息量，并采用統(tǒng)計分析法將信息量值分為3 個級別，結(jié)合區(qū)內(nèi)實際地質(zhì)情況、已有見礦工程分布等因素圈定了6 處找礦靶區(qū)，即A1－1，B3－1，B3－2，B3－3，C2－1，C2－2。A1－1區(qū)共包括預(yù)測含礦立方塊6 704 個，主要位于焦家主干斷裂、望兒山主斷裂和河西斷裂的深部產(chǎn)狀較緩處的巖體( 玲瓏巖體、馬連莊巖體、郭家?guī)X巖體) 接觸部位，黃鐵絹英巖化碎裂巖帶和斷裂的200 m 緩沖區(qū)內(nèi)及其兩側(cè)，預(yù)測礦體模型的產(chǎn)狀與焦家主干斷裂和河西斷裂大體一致，屬成礦有利部位。B3－1區(qū)共包括預(yù)測含礦立方塊2 119 個，主要位于焦家主干斷裂深部產(chǎn)狀較緩處的200 m 緩沖區(qū)內(nèi)及其兩側(cè)、黃鐵絹英巖化碎裂巖帶、馬連莊巖體、玲瓏巖體邊緣及接觸帶，預(yù)測礦體模型的產(chǎn)狀與焦家主干斷裂大體一致，屬成礦有利部位。B3－2區(qū)共包括預(yù)測含礦立方塊2 406 個，主要位于焦家主干斷裂產(chǎn)狀由陡變緩處的200 m 緩沖區(qū)內(nèi)及其兩側(cè)、黃鐵絹英巖化碎裂巖帶、玲瓏巖體與馬連莊巖體的接觸帶附近，預(yù)測礦體模型的產(chǎn)狀與焦家主干斷裂大體一致。B3－3區(qū)共包括預(yù)測含礦立方塊1 558 個，主要位于焦家主干斷裂南部與苗家斷裂帶交匯處的200 m 緩沖區(qū)內(nèi)及其兩側(cè)、玲瓏巖體接觸帶與鐵絹英巖化碎裂巖帶。C2－1區(qū)共包括預(yù)測含礦立方塊1 011 個，主要位于望兒山斷裂南部拐彎部位較淺處的200 m 緩沖區(qū)內(nèi)及其兩側(cè)、絹英巖化花崗質(zhì)閃長巖帶，預(yù)測礦體模型的產(chǎn)狀與望兒山斷裂大體一致。C2－2區(qū)共包括預(yù)測含礦立方塊260 個，主要位于望兒山主斷裂南部產(chǎn)狀由陡變緩處的200 m 緩沖區(qū)內(nèi)及其兩側(cè)，預(yù)測礦體模型的產(chǎn)狀與望兒山斷裂大體一致。采用體積估計法計算得到上述6 處找礦靶區(qū)的資源總量為382.482 8 t。

3 結(jié) 語

通過對山東焦家金礦帶地質(zhì)背景和成礦條件進行分析，總結(jié)出地質(zhì)找礦模型，在此基礎(chǔ)上建立了該區(qū)的定量預(yù)測模型，對區(qū)內(nèi)13 個預(yù)測變量進行了三維成礦有利條件分析與提取，采用三維信息法計算各地質(zhì)因素、找礦標(biāo)志所提供的找礦信息量，從而圈定了找礦靶區(qū)并估算了資源儲量，反應(yīng)出區(qū)內(nèi)具有較好的找礦潛力。

［1］ 石玉臣. 山東省焦家成礦帶深部金礦成礦預(yù)測研究及其應(yīng)用［D］. 長春:吉林大學(xué)，2005.

Shi Yuchen. Research on Deep Space Metallogenic Predication and its Application in Jiaojia Gold Belt in Shandong Province ［D］.Changchun:Jilin University，2005.

［2］ 宋明春，崔書學(xué)，楊之利. 山東焦家金礦帶深部找礦的重大突破及其意義［J］. 地質(zhì)與勘探，2008，44(1) :1-8.

Song Mingchun，Cui Shuxue，Yang Zhili. Great progress and far-reaching significance of deep exploration in Jiaojia metallogenic belt，Shandong Province［J］. Geology and Exploration，2008，44(1) :1-8.

［3］ 毛先成，鄒艷紅，陳 進，等. 危機礦山深部、邊部隱伏礦體的三維可視化預(yù)測——以安徽銅陵鳳凰山礦田為例［J］. 地質(zhì)通報，2010，29(2) :401-413.

Mao Xiancheng，Zou Yanhong，Cheng Jin，et al. Three-dimensional visual prediction of concealed ore bodies in the deep and marginal parts of crisis mines:a case study of the Fenghuangshan ore field in Tongling，Anhui Province［J］. Geological bulletin of China，2010，29(2) :401-413.

［4］ 張新宇. 地學(xué)空間三維可視化儲量計算輔助分析系統(tǒng)關(guān)鍵技術(shù)的研究［D］. 長春:吉林大學(xué)，2006.

Zhang Xinyu. Research on the Computer-aided System of Mineral Reserve Visualized Computation in Three-dimensional Geosciences Space［D］. Jilin:Jilin University，2006.

［5］ 陳建平，陳 勇，曾 敏，等. 基于數(shù)字礦床模型的新疆可可托海3 號脈三維定位定量研究［J］. 地質(zhì)通報，2008，27( 4) :552-559.

Chen Jianping，Chen Yong，Zeng Min，et al. 3D positioning and quantitative prediction of the Koktokay No.3 Dike，Xinjiang based on the digital mineral deposit model［J］. Geological Bulletin of China，2008，27(4) :552-559.

［6］ 王麗梅，陳建平，唐菊興. 基于數(shù)字礦床模型的西藏玉龍斑巖型銅礦三維定位定量預(yù)測［J］. 地質(zhì)通報，2010，29(4) :565-570.

Wang Limei，Chen Jianping，Tang Juxing.3D positioning and quantitative prediction of Yulong porphyry copper deposit，Tibet based on digital mineral deposit model［J］. Geological Bulletin of China，2010，29(4) :565-570.

［7］ 陳建平，嚴(yán) 瓊，尚北川，等. 湖南黃沙坪地區(qū)鉛鋅礦床三維預(yù)測研究［J］. 地質(zhì)學(xué)刊，2012，36(3) :243-249.

Chen Jianping，Yan Qiong，Shang Beichuan，et al. Research on 3D prediction of lead-zinc deposit in Huangshaping District in Hunan Province［J］. Journal of Geology，2012，36(3) :243-249.

［8］ 戎景會，陳建平，尚北川.基于找礦模型的云南個舊某深部隱伏礦體三維預(yù)測［J］. 地質(zhì)與勘探，2012，48(1) :191-198.

Rong Jinghui，Chen Jianping，Shang Beichuan. Three-dimensional prediction of blind ore bodies in Gejiu City，Yunnan Province based on the ore-search model［J］. Geology and Exploration，2012，48(1) :191-198.

［9］ 陳建平，于萍萍，史 蕊，等. 區(qū)域隱伏礦體三維定量預(yù)測評價方法研究［J］. 地學(xué)前緣，2014，21(5) :211-220.

Chen Jianping，Yu Pingping，Shi Rui. Research on three-dimensional quantitative prediction and evaluation methods of regional concealed ore bodies［J］. Earth Science Frontiers，2014，21(5) :211-220.

［10］ 趙鵬沄，顧雪祥，鄧小華. 山東焦家金礦礦床成因及成礦模式［J］. 地質(zhì)與勘探，2007，43(4) :29-35.

Zhao Pengyun，Gu Xuexiang，Deng Xiaohua. Genesis and metallogenic model of Jiaojia gold deposit，Shandong Province［J］. Geology and Exploration，2007，43(4) :29-35.

［11］ 劉玉強，楊東來，黃太嶺. 山東膠萊盆地金礦床地質(zhì)特征及找礦方向［J］. 礦床地質(zhì)，1999，18(3) :195-207.

Liu Yuqiang，Yang Donglai，Huang Tailing. Geological characteristics and ore-prospecting targets in gold deposits of Jiaolai Basin，Shandong Province［J］. Mineral Deposits，1999，18(3) :195-207.

［12］ 陳 勇. 云南個舊錫礦東區(qū)隱伏礦體定位與定量預(yù)測研究［D］. 北京:中國地質(zhì)大學(xué)( 北京) ，2009.

Chen Yong. Study on Localized and Quantitative Prediction of Concealed Deposits in the Eastern Gejiu Deposit［D］. Beijing: China University of Geosciences( Beijing) ，2009.

［13］ 呂 鵬. 基于立方體預(yù)測模型的隱伏礦體三維預(yù)測和系統(tǒng)開發(fā)［D］. 北京:中國地質(zhì)大學(xué)，2007.

Lu Peng. Development of the 3D Prediction System Based on Cube Prediction Model［D］. Beijing: China University of Geosciences( Beijing) ，2007.

［14］ 鄧金燦，甘文志. 基于Surpac 平臺設(shè)計的廣西高峰錫礦數(shù)字化信息系統(tǒng)三維空間分析及找礦預(yù)測［J］. 礦產(chǎn)與地質(zhì)，2009，23(3) :287-291.

Deng Jincan，Gan Wenzhi. Three-dimensional space analysis and exploration forecast with digital information system based on Surpac platform［J］. Mineral Resources and Geology，2009，23( 3) : 287-291.

［15］ 白 蕓. Surpac 在內(nèi)蒙古阿榮旗鉬礦中的應(yīng)用［D］. 北京: 中國地質(zhì)大學(xué)( 北京) ，2008.

Bai Yun. Application of Surpac in Molybdenum Deposit in Arongqi，Inaer Mgaolia［D］. Beijing: China University of Geosciences( Beijing) ，2008.

［16］ 楊曉坤，秦德先，馮美麗，等. 基于Surpac 的礦山三維地學(xué)模型及綜合信息成礦預(yù)測研究［J］. 地質(zhì)與資源，2008，17( 1) :61-64.

Yang Xiaokun，Qin Dexian，F(xiàn)eng Meili，et al.The Surpac-based 3D geologic model for mine and metalligenic prognosis by comprehensive information［J］. Geology and Resources，2008，17(1) :61-64.

［17］ 董慶吉，肖克炎，陳建平，等. 西南“三江”成礦帶北段區(qū)域成礦斷裂信息定量化分析［J］. 地質(zhì)通報，2010，29( 10) :1479-1485.

Dong Qingji，Xiao Keyan，Chen Jianping，et al. The quantitative analysis of regional metallogenic fault in the northern segment of the Sanjiang Metallogenic Belt，Southwestern China［J］. Geological Bulletin of China，2010，29(10) :1479-1485.

［18］ 龐緒成，顧雪祥，唐菊興，等. 山東焦家金礦深部盲礦體定位預(yù)測研究［J］. 成都理工大學(xué)學(xué)報: 自然科學(xué)版，2005，32( 6) :551-555.

Pang Xucheng，Gu Xuexiang，Tang Juxing，et al. Prediction of the deep blind ore bodies in the Jiaojia gold deposit，Shandong Province［J］. Journal of Chengdu University of Technology:Science ＆ Technology Edition，2005，32(6) :551-555.

［19］ 董 敏. 新城金礦構(gòu)造控礦及礦體三維定位預(yù)測［D］. 西安:長安大學(xué)，2010.

Dong Min. Structure of Ore-controlling and Steroprediction of Orebody Location of Xinchen Gold Deposit［D］. Xi’an:Chang'an University，2010.

［20］ 張 雋. 基于克里金方法的地球物理三維空間數(shù)據(jù)網(wǎng)格化［D］. 北京:中國地質(zhì)大學(xué)( 北京) ，2012.

Zhang Jun. Kring-based Geophysical Data Mesh of Three-dimensional Space［D］. Beijing: China University of Geosciences ( Beijing) ，2012.

［21］ 毛先成，戴塔根，吳湘濱，等. 危機礦山深邊部隱伏礦體立體定量預(yù)測研究——以廣西大廠錫多金屬礦床為例［J］. 中國地質(zhì)，2009，36(2) :424-435.

Mao Xiancheng，Dai Tagen，Wu Xiangbin，et al. Stereoscopic quantitative prediction of concealed ore bodies in the deep and marginal parts of crisis mines: a case study of the Dachang tin polymetallic ore deposit in Guangxi Province［J］. Geology in China，2009，36(2) :424-435.