劉漢棟，王巧云，李秀章，陳建平，祝培剛，朱繼托，劉 麗

（1．山東省地質(zhì)調(diào)查院，山東 濟南250013；2．中國地質(zhì)大學(xué) （北京）國土資源與高新技術(shù)研究中心，北京100083；3．北京市國土資源信息開發(fā)研究重點實驗室，北京100083）

成礦預(yù)測是應(yīng)用地質(zhì)理論和科學(xué)方法，綜合地質(zhì)、地球物理、地球化學(xué)等基礎(chǔ)資料獲得找礦信息，總結(jié)成礦條件和礦體賦存規(guī)律，建立礦床模型，圈定不同級別成礦遠景區(qū)［1］。隨著科學(xué)技術(shù)的發(fā)展，成礦預(yù)測評價由定性向定量、由二維向三維發(fā)展，20世紀(jì)80年代以來，地理信息系統(tǒng)（GIS）技術(shù)被成功引入成礦預(yù)測中，應(yīng)用GIS技術(shù)，提高了礦產(chǎn)資源預(yù)測的效率，縮短了資源評價周期，并能定量圈定各級有利成礦區(qū)段［2－6］。證據(jù)權(quán)模型是一種立足于數(shù)據(jù)的空間位置關(guān)系，結(jié)合GIS技術(shù)對各種有利成礦因素（證據(jù)層）進行有效綜合的數(shù)學(xué)預(yù)測評價模型［7］。目前國內(nèi)不同專家、學(xué)者對一些主要成礦區(qū)帶、礦床進行了基于證據(jù)權(quán)法的礦產(chǎn)資源預(yù)測工作，取得了良好的找礦成果［4－5，7－17］。

本文在研究焦家成礦帶地質(zhì)背景和成礦規(guī)律的基礎(chǔ)上，利用GIS技術(shù)進行成礦有利信息分析與提取，建立評價證據(jù)指標(biāo)模型，使用證據(jù)權(quán)法對金礦進行了成礦預(yù)測，指出該區(qū)具有進一步的找礦潛力。

1 研究區(qū)地質(zhì)背景

研究區(qū)位于膠東半島西北部，西靠沂沭斷裂帶，南接膠萊拗陷，北鄰龍口斷陷盆地和渤海拗陷，東接牟平－即墨構(gòu)造混雜巖帶［18］。古老的前寒武紀(jì)變形變質(zhì)巖系、多期多成因的巖漿活動和以NE－NNE向斷裂為主的構(gòu)造格架，構(gòu)成了本區(qū)金礦的區(qū)域成礦地質(zhì)背景（圖1）。

圖1 膠東西北部地區(qū)地質(zhì)礦產(chǎn)簡圖

區(qū)域主要控礦構(gòu)造為招平斷裂、焦家斷裂和三山島斷裂，斷裂總體呈NE走向，與三條主干斷裂平行展布的一組近NE向次級斷裂發(fā)育，少量NW向和近EW向斷裂。三山島斷裂傾向SE，焦家斷裂傾向NW，招平斷裂傾向SE，三山島斷裂和焦家斷裂之間及招平斷裂東側(cè)的斷裂上盤主要地質(zhì)體為早前寒武紀(jì)變質(zhì)雜巖，由新太古代TTG花崗質(zhì)片麻巖、變輝長巖（斜長角閃巖）和少量呈殘留包體出現(xiàn)的膠東巖群及中太古代唐家莊巖群組成；焦家斷裂和招平斷裂之間的斷裂下盤主要分布侏羅紀(jì)玲瓏花崗巖，侵入有白堊紀(jì)郭家?guī)X花崗閃長巖、偉德山花崗巖、嶗山花崗巖和相關(guān)的脈巖（煌斑巖、輝綠玢巖、閃長玢巖、花崗質(zhì)脈巖等）。主干斷裂切割了玲瓏花崗巖體［19］。

焦家主干斷裂及次級低序次斷裂（望兒山斷裂、河西斷裂等）構(gòu)成控礦構(gòu)造系統(tǒng)，焦家斷裂總體走向30°，傾向北西，大致沿新太古代變質(zhì)巖系與侏羅紀(jì)玲瓏花崗巖的接觸帶分布，但有的地段穿切玲瓏花崗巖和郭家?guī)X花崗巖，破碎蝕變帶寬80～200m不等［20］。

在焦家金成礦帶上，先后發(fā)現(xiàn)焦家、新城、河西、馬塘、望兒山、上莊等大型金礦床，近幾年先后在成礦帶的深部發(fā)現(xiàn)寺莊深部、馬塘深部、焦家深部等大型－超大型金礦床，礦床成因類型主要為破碎帶蝕變巖型，少量石英脈型。

2 成礦有利信息提取

本次研究收集編制了多元地質(zhì)數(shù)據(jù)庫，包括地質(zhì)、地球物理、地球化學(xué)、礦產(chǎn)勘查數(shù)據(jù)等。根據(jù)成礦條件，基于GIS技術(shù)分析了研究區(qū)的成礦地質(zhì)背景和已查明的十余個金礦床的實測數(shù)據(jù)，從建立的數(shù)據(jù)庫中提取了6個成礦有利證據(jù)層。

2．1 構(gòu)造信息

從構(gòu)造緩沖區(qū)、構(gòu)造展布特征、構(gòu)造交匯特征、構(gòu)造中心對稱度等對構(gòu)造特進行成礦有利信息的提取與分析。

構(gòu)造緩沖區(qū)：焦家主干斷裂破碎帶發(fā)育，是金礦體產(chǎn)出的主要位置，結(jié)合研究區(qū)實際情況，選取斷裂面兩側(cè)200m為緩沖區(qū)，經(jīng)統(tǒng)計，87．07%的礦塊落在緩沖區(qū)內(nèi)，斷裂緩沖區(qū)是一個非常重要的預(yù)測要素。

構(gòu)造展布特征：主要包括斷裂等密度、斷裂頻數(shù)、主干斷裂分析等，這些變量從不同的角度反映線性構(gòu)造的特征，從中發(fā)掘與成礦有關(guān)的信息，提取成礦信息是成礦預(yù)測的要求。構(gòu)造等密度反映線性構(gòu)造發(fā)育程度，是斷裂長度的加和；構(gòu)造頻數(shù)是單元中斷裂構(gòu)造產(chǎn)出的條數(shù)，是斷裂條數(shù)的加和，反映了區(qū)域線形構(gòu)造的復(fù)雜程度，體現(xiàn)了區(qū)域構(gòu)造格架的主體特征；主干斷裂是指傾向上延深大，平面上延伸長的斷裂構(gòu)造。采用斷裂等密度與斷裂頻數(shù)之比來定量化分析主干斷裂發(fā)育區(qū)。其比值大表現(xiàn)為主干斷裂的特征。以焦家金礦帶為背景，北東向斷裂（焦家斷裂、望兒山斷裂和河西斷裂）為焦家金礦帶內(nèi)的主干斷裂。圖2為研究區(qū)內(nèi)斷裂的構(gòu)造長度、頻數(shù)玫瑰花圖。經(jīng)統(tǒng)計選取主干斷裂值（0．00917，0．01167）的范圍作為成礦有利區(qū)間（圖3）。

圖2 研究區(qū)構(gòu)造長度、頻數(shù)玫瑰花圖

圖3 主干斷裂（長度／頻度）分布直方圖

構(gòu)造交匯特征：斷裂交匯部位是成礦有利位置，構(gòu)造交點數(shù)代表著單元面積內(nèi)交點數(shù)多少的一個定量指標(biāo)，經(jīng)統(tǒng)計選取構(gòu)造交匯點值為（0．005，1．250）的范圍作為成礦有利區(qū)間。

構(gòu)造中心對稱度：代表了構(gòu)造對稱的特征。經(jīng)統(tǒng)計，選?。?．001，0．020）為成礦有利因子，研究區(qū)北部有利區(qū)間塊體周圍區(qū)域較好的成礦事實預(yù)示著南部相似部位具有很好的找礦前景。

2．2 蝕變巖信息

研究區(qū)金礦體均賦存于斷裂構(gòu)造所控制的蝕變巖內(nèi)，蝕變巖的發(fā)育空間、形態(tài)、產(chǎn)狀與斷裂帶基本一致。本文通過統(tǒng)計已知礦體分布于不同蝕變巖的比例來確定不同蝕變巖對成礦的貢獻大小。經(jīng)統(tǒng)計，研究區(qū)內(nèi)礦體41．38%落在絹英巖化花崗質(zhì)碎裂巖內(nèi)，26．88%在黃鐵絹英巖化碎裂巖內(nèi)，16．06%在絹英巖化花崗巖內(nèi)，5．08%落在變輝長質(zhì)碎裂巖帶內(nèi)（圖4）?？梢?，絹英巖化花崗質(zhì)碎裂巖和黃鐵絹英巖化碎裂巖是礦體產(chǎn)出的主要蝕變巖。蝕變巖的分布圖層可作為成礦有利證據(jù)層。

圖4 不同蝕變巖中含礦單元圖

2．3 成礦有利巖體

研究區(qū)內(nèi)侵入巖主要為古元古代馬連莊序列變輝長巖、中生代侏羅紀(jì)玲瓏序列花崗巖和白堊紀(jì)郭家?guī)X序列花崗閃長巖。經(jīng)統(tǒng)計，研究區(qū)內(nèi)62．26%的金礦體落在玲瓏序列內(nèi)，31．59%礦體落在馬連莊序列，6．14%在郭家?guī)X序列（圖5），顯示玲瓏序列和馬連莊序列中含礦性較好。

2．4 等間距控礦信息

構(gòu)造帶和構(gòu)造形跡空間展布的韻律定向性和間距的倍數(shù)性是構(gòu)造距離的兩個特點［21］。區(qū)內(nèi)控礦構(gòu)造的等距性和礦液流動方向及疊加程度導(dǎo)致了礦床呈等距性展布，對指導(dǎo)盲區(qū)找礦具有指示性作用［22］。

同產(chǎn)于焦家主干斷裂中的新城、東季、焦家、馬塘、寺莊等金礦床，雖然它們的規(guī)模大小不一，但其礦床成因類型、蝕變巖、礦體形態(tài)及產(chǎn)狀、礦體側(cè)伏及礦化分帶規(guī)律都是相同的。盧作祥等［23］在對本區(qū)研究后，得出了同一控礦構(gòu)造帶內(nèi)的礦床（點）具有等距性分布的規(guī)律。近年來的研究及新城、望兒山、焦家、馬塘、寺莊等金礦床第二富集帶礦體的發(fā)現(xiàn)表明，這種等距性不只是在走向上，也存在于傾向上。經(jīng)統(tǒng)計，研究區(qū)內(nèi)93．91%的金礦床落在等間距控礦模型中（圖6）。

圖5 不同巖體中含礦單元圖

圖6 焦家金礦帶主要礦床分布圖

2．5 地球物理成礦信息

將研究區(qū)可控源大地音頻電磁測量資料與區(qū)內(nèi)已知礦體進行疊加統(tǒng)計分析，發(fā)現(xiàn)59．7%的金礦體落在視電阻率值為1000～20000Ω·m的區(qū)間，顯示視電阻率值1000～20000Ω·m的區(qū)間為金礦床有利地球物理異常信息。

2．6 地球化學(xué)成礦信息

化探信息是一種直接的找礦信息，是找礦預(yù)測中的直接變量［24］。利用研究區(qū)鉆孔中Au元素分析數(shù)據(jù)和原生暈化探數(shù)據(jù)分析單元塊的元素異常分布，使用距離反比加權(quán)對未知單元塊進行Au元素含量插值處理。依據(jù)鉆孔原生暈化探數(shù)據(jù)信息，建立 Au、Ag、As、Bi、Cu、Hg、Pb、Sb、Zn等微量元素異常模型，依據(jù)元素縱向分帶規(guī)律確定As、Hg、Sb為礦體頭部異常元素；Ag、Au、Zn、Pb、Cu為近礦異常元素；Bi、Mo為礦體尾部異常元素。經(jīng)統(tǒng)計，鉆孔數(shù)據(jù)覆蓋范圍內(nèi)金礦床78．49%落在Au元素異常內(nèi)，37．77%的金礦床落在元素組合異常內(nèi)。

3 基于證據(jù)權(quán)模型的多元綜合信息成礦預(yù)測

3．1 預(yù)測模型的建立

依據(jù)成礦地質(zhì)條件和找礦標(biāo)志，基于GIS技術(shù)的找礦模型和成礦有利信息處理的基礎(chǔ)上，對研究區(qū)金礦成礦的控礦因素進行了綜合分析，總結(jié)出與研究區(qū)成礦有關(guān)的預(yù)測變量有14個，包括賦礦巖體（玲瓏序列、馬連莊序列）、蝕變帶（黃鐵絹英巖化碎裂巖、絹英巖化花崗質(zhì)碎裂巖、變輝長質(zhì)碎裂巖），構(gòu)造（斷裂緩沖區(qū)、構(gòu)造交點數(shù)、構(gòu)造中心對稱度、主干斷裂）、Au品位及元素組合異常、等間距控礦、視電阻率異常，研究區(qū)定量預(yù)測模型見表1。

3．2 成礦證據(jù)因子分析

本次研究采用證據(jù)權(quán)法對有利成礦信息和定量預(yù)測模型進行綜合分析。證據(jù)權(quán)法是加拿大數(shù)學(xué)地質(zhì)學(xué)家Agterberg提出的一種地學(xué)統(tǒng)計方法，是通過對一些與礦產(chǎn)形成相關(guān)的地學(xué)信息的疊加復(fù)合分析來進行礦產(chǎn)遠景區(qū)的預(yù)測。其中的每一種地學(xué)信息都被視為成礦遠景區(qū)預(yù)測的一個證據(jù)因子，而每一個證據(jù)因子對成礦預(yù)測的貢獻是由這個因子的權(quán)重值來確定的［17］。

依據(jù)成礦規(guī)律，根據(jù)對地質(zhì)、地球物理和地球化學(xué)成礦信息的分析，對研究區(qū)的各控礦因素做了綜合評定，提取出14個有利證據(jù)項：Au元素異常、元素組合異常、變輝長雜巖質(zhì)碎裂巖、黃鐵絹英巖化碎裂巖、等間距控礦、斷裂緩沖區(qū)、絹英巖化花崗質(zhì)碎裂巖、斷裂、主干斷裂、玲瓏序列、斷裂中心對稱度、馬連莊序列、視電阻率異常、構(gòu)造交匯點。針對以上14項，利用MRAS軟件進行了證據(jù)權(quán)法計算。各證據(jù)因子的權(quán)重值見表2。

對表2數(shù)據(jù)分析，得出以下4點認識：①研究區(qū)Au元素異常和 Au、Ag、Cu、Pb、Zn、As、Sb、Bi、Hg化探元素組合異常與成礦關(guān)系密切，其C值皆大于5；②蝕變碎裂巖與金礦關(guān)系密切，變輝長質(zhì)碎裂巖、黃鐵絹英巖化碎裂巖、絹英巖化花崗質(zhì)碎裂巖其C值皆大于3；③斷裂構(gòu)造作為賦礦部位，與成礦關(guān)系密切，等間距控礦、斷裂緩沖區(qū)的C值大于4，主干斷裂的C值大于3，斷裂中心對稱度C大于2；④賦礦圍巖與成礦有一定的關(guān)系，其C值大于2。

對14個證據(jù)項進行條件獨立性檢驗，在顯著性水平為0．05的情況下，上述14個因素基本上滿足條件獨立性。

表1 焦家金成礦帶預(yù)測模型

表2 焦家金礦帶各證據(jù)因子權(quán)重值

4 預(yù)測區(qū)的圈定及驗證

依據(jù)證據(jù)權(quán)重分析，計算出各單元含礦預(yù)測的后驗概率，結(jié)合實際地質(zhì)情況和礦體分布規(guī)律，在已有礦區(qū)、礦體范圍外圈定找礦靶區(qū)。研究區(qū)內(nèi)圈定6個找礦預(yù)測區(qū)（圖7）。

1）后坡預(yù)測區(qū)：位于焦家斷裂、望兒山斷裂深部延伸部位，淺部為曲家超大型金礦床和焦家深部超大型金礦床，二個礦床主礦體深部尚未封閉。預(yù)測模型主要落在玲瓏巖體、黃鐵絹英巖化碎裂巖以及斷裂的200m緩沖區(qū)內(nèi)，預(yù)測模型的礦體產(chǎn)狀與焦家主干斷裂一致，為成礦有利部位。

圖7 焦家金成礦帶成礦預(yù)測圖

2）任家預(yù)測區(qū)：位于焦家斷裂深部產(chǎn)狀變緩部位，紗嶺超大型金礦床的深部，預(yù)測模型主要落在馬連莊巖體和玲瓏巖體的接觸帶附近，焦家主干斷裂的200m緩沖區(qū)、黃鐵絹英巖化碎裂巖帶和絹英巖化花崗質(zhì)碎裂巖帶內(nèi)，預(yù)測模型的礦體產(chǎn)狀與焦家主干斷裂斷裂一致，成礦條件有利。

3）臥龍預(yù)測區(qū)：位于焦家主干斷裂深部，淺部為前陳中型金礦床，預(yù)測模型主要落在斷裂的200m緩沖區(qū)和黃鐵絹英巖化碎裂巖內(nèi)，位于玲瓏巖體與馬連莊巖體的接觸帶附近，預(yù)測模型的礦體產(chǎn)狀與焦家主干斷裂一致，有較好的成礦地質(zhì)條件。

4）綦家預(yù)測區(qū)：位于焦家斷裂南延部位，預(yù)測區(qū)及附近尚未發(fā)現(xiàn)金礦床（點），預(yù)測立方塊主要落在斷裂的200m緩沖區(qū)內(nèi)，位于玲瓏巖體與馬連莊巖體接觸帶附近，有視電阻率異常存在，成礦條件較好。

5）龍埠預(yù)測區(qū)：位于望兒山斷裂帶的深部，其南側(cè)淺部為望兒山中型金礦床，預(yù)測模型主要落在斷裂的200m緩沖區(qū)和絹英巖化花崗質(zhì)碎裂巖、玲瓏巖體內(nèi)，預(yù)測模型的礦體產(chǎn)狀與望兒山斷裂大體一致。

6）城子預(yù)測區(qū)：位于望兒山斷裂南延部位，其北部為望兒山金礦床和城北小型金礦床，預(yù)測模型主要落在斷裂的200m緩沖區(qū)、玲瓏巖體內(nèi)，預(yù)測模型的礦體產(chǎn)狀與望兒山斷裂大體一致。

在圈定的后坡預(yù)測區(qū)內(nèi)施工1個驗證鉆孔，在1328．50～1355．30m見到品位大于0．5×10－6的礦化體，視厚度26．80m，圈定礦體3個，主礦體位于孔深1344．61～1349．11m 處，視厚度4．50m，平均品位2．63×10－6，驗證結(jié)果與預(yù)測基本一致。

5 結(jié) 論

1）利用GIS技術(shù)對研究區(qū)開展多源地質(zhì)找礦信息綜合集成分析，提取斷裂構(gòu)造、巖漿巖、物化探異常等與成礦的關(guān)系，結(jié)合區(qū)域成礦地質(zhì)條件和礦體分布規(guī)律，總結(jié)出研究區(qū)內(nèi)有利于金礦產(chǎn)出的成礦條件組合，提出了研究區(qū)控礦要素為巖漿巖、構(gòu)造、蝕變帶、物探異常、化探異常和等間距控礦等。

2）結(jié)合研究區(qū)地質(zhì)特征、區(qū)域成礦規(guī)律和成礦地質(zhì)條件，應(yīng)用基于GIS的證據(jù)權(quán)法進行多源成礦信息綜合預(yù)測評價，提出了6個金礦找礦預(yù)測區(qū)，為研究區(qū)勘查工作部署提供了依據(jù)。

3）對預(yù)測區(qū)擇優(yōu)進行鉆探驗證，在焦家主干斷裂的深部發(fā)現(xiàn)了主礦體，驗證了預(yù)測結(jié)果。

［1］ 肖克炎，張曉華，陳鄭輝，等．成礦預(yù)測中證據(jù)權(quán)重法與信息量法及其比較［J］．物探化探計算技術(shù)，1999，21（3）：223－226．

［2］ 向運川，任天祥，楊竹溪．開發(fā)利用地理信息系統(tǒng)（GIS）綜合分析地學(xué)信息進行礦產(chǎn)預(yù)測［J］．物探與化探，1996，20（1）：1－13．

［3］ 池順都，趙鵬大．應(yīng)用GIS圈定找礦可行地段和有利地段［J］．地球科學(xué)，1998，23（2）：125－128．

［4］ 劉星，胡光道．應(yīng)用MORPAS系統(tǒng)證據(jù)權(quán)重法進行多源信息成礦預(yù)測——以瀾滄江南段地區(qū)為例［J］．地質(zhì)與勘探，2003，39（4）：44－68．

［5］ 蘇紅旗，葛艷，劉冬林，等．基于GIS的證據(jù)權(quán)重法礦產(chǎn)成礦預(yù)測系統(tǒng)（EWM）［J］．地質(zhì)與勘探，1999，35（1）：44－46．

［6］ 曹瑜，胡光道．圈定“5P”找礦地段的GIS成礦預(yù)測空間模型及其應(yīng)用［J］．地球科學(xué)，1999，24（4）：409－412．

［7］ 呂鵬，朱鵬飛，畢志偉，等．基于GIS和證據(jù)權(quán)模型的克什克騰旗有色金屬成礦預(yù)測與評價［J］．地質(zhì)與勘探，2011，47（5）：909－917．

［8］ 陳建平，董慶吉，郝金華，等．基于GIS的證據(jù)權(quán)重法青?！叭北倍伟邘r型鉬銅礦產(chǎn)資源成礦預(yù)測［J］．巖石礦物學(xué)雜志，2011，30（3）：519－529．

［9］ 徐善法，陳建平，葉繼華．證據(jù)權(quán)法在三江北段銅金礦床成礦預(yù)測中的應(yīng)用研究［J］．地質(zhì)與勘探，2006，42（2）：54－59．

［10］ 申維，堵海燕，依力哈木，等．基于GIS的證據(jù)權(quán)法在浙西銅礦資源成礦預(yù)測中的應(yīng)用［J］．地質(zhì)學(xué)刊，2014，38（2）：245－253．

［11］ 孫巖，王訓(xùn)練，陳建平，等．基于證據(jù)權(quán)重法的沱沱河地區(qū)鉛鋅銀礦成礦預(yù)測［J］．地質(zhì)通報，2010，29（4）：556－564．

［12］ 鄧勇，邱瑞山，羅鑫．基于證據(jù)權(quán)重法的成礦預(yù)測——以廣東省鎢錫礦的成礦預(yù)測為例［J］．地質(zhì)通報，2007，26（9）：1228－1234．

［13］ 劉世翔，薛林福，郄瑞卿，等．基于GIS的證據(jù)權(quán)重法在黑龍江省西北部金礦成礦預(yù)測中的應(yīng)用［J］．吉林大學(xué)學(xué)報（地球科學(xué)版），2007，37（5）：889－894．

［14］ 劉曉玲，陳建平．基于GIS的證據(jù)權(quán)重法在內(nèi)蒙古阿魯科爾沁旗地區(qū)成礦預(yù)測中的應(yīng)用［J］．地質(zhì)通報，2010，29（4）：571－580．

［15］ 張啟明，陳建平，齊先茂．基于GIS的證據(jù)權(quán)法在三江北段銅多金屬成礦預(yù)測中的應(yīng)用［J］．沉積與特提斯地質(zhì)，2007，27（2）：106－112．

［16］ 成秋明，陳志軍，Khaled A．模糊證據(jù)權(quán)方法在鎮(zhèn)沅（老王寨）地區(qū)金礦資源評價中的應(yīng)用［J］．地球科學(xué)：中國地質(zhì)大學(xué)學(xué)報，2007，32（2）：175－184．

［17］ 陳建平，陳勇，王全明．基于GIS的多元信息成礦預(yù)測研究——以赤峰地區(qū)為例［J］．地學(xué)前緣，2008，15（4）：18－26．

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

［19］ 宋明春，伊丕厚，徐軍祥，等．膠西北金礦階梯式成礦模式［J］．中國科學(xué)：地球科學(xué)，2012，42（7）：992－1000．

［20］ 宋明春，宋英昕，沈昆，等．膠東焦家深部金礦礦床地球化學(xué)特征及有關(guān)問題討論［J］．地球化學(xué)，2013，42（3）：274－289．

［21］ 翟裕生．關(guān)于礦田構(gòu)造研究的若干問題［J］．地質(zhì)論評，1984，30（1）：19－25．

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

［23］ 盧作祥．山東招掖金礦帶焦家式金礦床的礦床分帶［J］．地球科學(xué)，1998，13（2）：54－56．

［24］ 董英君，張德全，徐文藝，等．礦產(chǎn)資源評價系統(tǒng)及其在東昆侖的應(yīng)用［J］．地球?qū)W報，2005，26（1）：83－88．