雷明選,魏明非,高 虔,孫 青,惠偉光
(1.陜西省地礦局 第六地質隊,臨潼 710611;2.中國地質大學,北京 100083)
馬角寺鉬礦位于小秦嶺西段南坡,行政區劃屬商洛市商州區腰市鎮管轄。該地區褶皺及斷裂構造復雜,成礦地質條件良好。上世紀八十年代在小秦嶺地區開展的1∶50 000 水系沉積物測量工作,曾在本區圈定出綜合異常一處,異常分布面積較大,組分復雜,主要元素組合為Au、Pb、Zn、Mo,顯示該區有較好的Mo、Au多金屬找礦前景。之前曾有多家地勘單位在本區開展過異常查證及地質找礦工作,但由于該區地形惡劣,地表覆蓋嚴重,且工作中使用的勘查手段較為單一,一直未取得重大突破。
2010年陜西地礦局第六地質隊開展了涵蓋本區的“陜西省洛南縣洛源多金屬礦普查”工作,綜合運用大比例尺地質測量、土壤地球化學測量、巖石地化剖面測量及槽探、鉆探工程驗證等工作手段和方法,發現了該礦床。尤其是利用土壤地球化學測量成果對該區土壤地球化學特征進行總結和研究,在縮小找礦靶區直至最終發現該礦床過程中起到了至關重要的作用,成為運用地球化學方法成功找礦的一個典型實例,在今后其他地區的找礦工作中值得借鑒。
目前該礦區共圈出鉬礦體四個。礦床中鉬為主成礦元素,伴生釩礦。經研究該礦區外圍及鉬礦體深部找礦潛力較大,有進一步找礦前景。
馬角寺鉬礦區位于華北陸塊南緣的豫皖陸塊區西南邊緣,成礦區劃屬古亞洲成礦域華北陸塊成礦省小秦嶺—豫西太古代、元古代、古生代、中生代金鉬鋁土礦鉛鋅成礦區(Ⅲ-17),劃歸為金堆城元古—早古生代上疊盆地燕山期鉬鐵銅鉛黃鐵礦成礦帶(Ⅲ-172)金堆城—木龍溝鉬鐵銅鉛鋅黃鐵礦成礦遠景區(Ⅲ-172-1)。
礦區及其外圍出露地層由老到新依次為:元古界、古生界、新生界等(圖1)。
(1)元古界分布于礦區北部,出露地層包括:①長城系熊耳群斜長玢巖、細碧巖等淺變質火山巖系;②長城系高山河組石英巖、變石英砂巖、變粉砂巖等淺變質濱海—淺海相碎屑巖;③薊縣系龍家園組硅質灰巖、燧石條帶灰巖等淺海相鎂質碳酸鹽巖;④薊縣系巡檢組海相沉積的燧石條帶狀白云巖,其下部為土黃色粘土質含礫白云巖,底部有一層厚十幾厘米的褐鐵礦風化殼;⑤震旦系砂質板巖、長石粉砂巖,下部為冰磧礫巖。
(2)古生界包括為寒武系(∈)和奧陶-志留系陶灣群((O—S)tw)地層。寒武系分布于礦區中部,為一套淺變質的淺海相碳酸鹽巖層夾碎屑巖建造,按巖性組合特征分為上、中、下三個統。陶灣群((O—S)tw)出露于礦區南部,為一套淺變質的淺海相碎屑巖、碳酸鹽巖夾中基性火山巖建造。
高山河組與下覆地層熊耳群、上伏地層龍家園組均為角度不整合接觸,震旦系與下覆巡檢司組為角度不整合接觸,陶灣群與下覆的寒武系為角度不整合接觸,而震旦系與寒武系之間為斷層接觸。
含鉬礦化構造帶及鉬礦體賦存于震旦系與寒武系地層的斷層接觸帶。礦區構造表現為在南北向主壓應力持續作用下,形成東西向緊閉線性褶皺,相伴而生東西向逆沖斷層以及與其相配套的北東向、北西向平移斷層。
褶皺構造:主要表現為東西向輞峪緊閉倒轉向斜,其軸線走向為300°,向東傾伏,傾伏角40余度。北翼出露元古界及寒武系地層,傾向北北東,傾角40°~70°,地層倒轉,指示倒轉的劈理隨處可見。南翼出露寒武系地層,層序正常,傾向北北東,傾角20°~45°,多處地段發育有舒緩的次一級褶曲。核部出露地層為陶灣群廟灣組。馬角寺鉬礦床位于該緊閉倒轉向斜的北翼。
礦區斷裂較為發育,以東西向或近東西為主,其次為北東向及北西向。東西向斷裂:傾向北或南,斷層帶內劈理發育,可見擠壓構造片巖、擠壓構造圍巖透鏡體、圍巖角礫、斷層泥等,反映斷裂具多期活動的特征,為區內主要控礦、儲礦構造,帶內常見石英細脈群充填,構成含鉬礦化構造帶,馬角寺鉬礦床賦存于其中。北東向斷層為左行平移斷層,北西向斷層為右行平移斷層。區內巖漿巖不發育,脈巖主要為石英脈、長英質脈等,多分布于東西向斷裂帶內。礦區周邊地區金屬礦產分布較多,主要有金堆城特大型鉬礦床、鐵爐子鉛鋅礦床、鐵源鉛鋅礦床、火石溝金礦床、潘河鉬礦床等。

圖1 馬角寺鉬礦區地質概略圖Fig.1 Geological sketch map of mining area in Majiaosi molybdenum deposit area
1∶50 000水系沉積物測量在區內圈定出綜合異常一處,面積約44 km2,分布于東西向、北東向及北西向斷裂帶附近及其交匯部位(圖2及表1)。綜合異常分布面積較大,組分復雜,主要元素組合為Au、Pb、Zn、Mo,其次為Cu。各元素異常位置大致吻合,元素間套合性好。異常區褶皺及斷裂構造復雜,顯示有較好的Mo、Au多金屬找礦前景。

圖2 馬角寺鉬礦區1∶50 000水系沉積物地球化學異常Fig.2 Geochemical anomalies(1∶50 000) of drainage deposits in Majiaosi molybdenum deposit area

元素面積異常下限極大值均值濃度分帶備注Mo6.1526.05.83為主要成礦元素Au6.215.62.33Pb15.8100520150.62為主要伴生元素Zn5.91503151562Cu202040.726.81為高背景元素
注:Au元素含量單位為10-9,其余元素含量單位為10-6
1∶10 000土壤測量范圍位于1∶50 000水系沉積物異常區內,東西向長約10 km,南北向寬約4.3 km,面積約40 km2。區內地層及構造線近東西向展布,樣品布置以南北方向布置采樣線,采樣網度100 m×40 m。土壤樣品采取殘坡積層B層,由一點多坑等份量采集組合成一個樣品,采樣深度10 cm~30 cm。共采集土壤樣品10 306件,所有樣品均晾曬干燥后過60#尼籠篩,送樣重量大于80 g,送國土資源部西安礦產資源監督檢測中心分析,分析項目為Au、Ag、Cu、Pb、Zn、Mo等六種元素。運用化學光譜法分析Au;運用原子吸收法分析Ag、Cu、Pb、Zn、Mo。
3.2.1 元素含量特征
對各元素平均值、最大值、最小值、標準離差、變化系數及濃度克拉克值等地球化學參數進行了統計(表2)。
從各元素地球化學參數總體特征來看,區內Mo、Pb、 Au等元素含量相對較高,變化系數較大,具有成礦的良好地球化學特征。
3.2.2 元素含量分布型式
圖3為區內各元素含量對數—頻數直方圖。從圖上可以看出,區內Au、Ag、Cu、Pb、Zn、Mo等六種元素含量分布均為單峰形態,分配型式均近似服從對數正態分布,其中Mo、Au、Ag、Pb、Zn等元素含量的分布型式均向右邊高含量區域傾斜,分布范圍較廣,離散程度較大,表明這些元素參與了次生富集成暈作用及過程。Cu元素含量的分布型式均向左邊低含量區域傾斜,表明其離散度相對較小,該元素在次生富集成暈作用及過程貢獻較小。
3.3.1 元素相關性分析
對區內土壤樣品Au、Ag、Cu、Pb、Zn、Mo等元素進行了相關性分析(表3)。可以看出,Mo、Ag、Pb、Zn間為呈較強的正相關,顯示Mo與Ag、Pb、Zn可能來自同一成礦階段;Au與Ag、Pb、Zn、Mo相關性不明顯,顯示Au的成礦與Ag、Pb、Zn、Mo不屬于同一成礦階段;Cu與其他元素極弱相關或不相關,說明Cu與上述各元素的成礦關系不強。

表2 馬角寺鉬礦區各元素地球化學參數統計一覽表
注:樣品數量n=10 396;Au、Ag元素含量單位為10-9,其余元素含量單位為10-6;測試單位為國土資源部西安礦產資源監督檢測中心。

圖3 馬角寺鉬礦區元素含量對數—頻數直方圖Fig.3 Logarithm-frequency histogram of element content in Majiaosi molybdenum deposit area

表3 馬角寺鉬礦區元素相關系數矩陣
3.3.2 R型聚類分析
圖4為本區土壤樣品各元素相關性聚類分析譜系圖,從中可以清楚地看出元素之間的組合關系。在相關系數為0.5的水平下,Mo、Ag、Pb、Zn組成一類,而Au、Cu各自單獨成為一類。由此進一步說明了Mo、Ag、Pb、Zn系同一成礦階段的產物。

圖4 馬角寺鉬礦區元素聚類分析譜系圖Fig.4 The pedigree chart of element cluster analysis in Majiaosi molybdenum deposit area
分形理論創立于上世紀七十年代中期,由法國數學家曼德爾布羅特(Mandelbrot)提出,其研究對象為自然界和社會活動中廣泛存在的無序(無規則)而具有自相似性的系統。該理論自問世以來,已被廣泛應用于地球科學的眾多領域,并取得了眾多研究成果。在礦床學領域,分形已被應用于確定區域化探異常、多期次礦化疊加[5-7]、礦床遠景儲量預測[8-10]和指導勘探網度布置[11]等諸多領域。
在這里我們采用含量-頻數法對本區土壤元素含量數據進行分形分析,探討分維值與地球化學異常的相關性,并嘗試利用分維值計算元素背景值,從而確定異常下限。
分數維的定義:D=-lgN(c)/lgc,其中:c表示元素含量;N(c)表示元素含量大于c的樣品數;D為分形的維數。將c和N(c)標繪在雙對數圖上,元素含量各點大致在一條線上分布,利用直線的斜率可求出D值[17]。若元素含量各點大致在一條線上分布,則反映了元素含量之間具有自相似性[3]。
我們對本區10 306件土壤樣品中Au、Ag、Cu、Pb、Zn、Mo等六種元素含量數據進行分形分析,將c和N(c)標繪在雙對數圖上,繪制成元素含量分形曲線圖(圖5)。
從圖5可以看出,Cu、Pb、Zn等三種元素含量分形曲線特征基本相似,均表現為明顯的兩段式分布,其低含量構成的分維值D1反映各元素的背景值分布,高含量構成的分維值D2反映各元素的土壤地球化學異常;Au、Ag、Mo等三種元素含量分形曲線呈明顯的分段式分布,具多標度分形特征,低含量構成的分維值D1反映各元素的背景值分布,D2、D3、D4反映了本區Mo、Au、Ag元素疊加有不同程度和作用的地球化學異常分布。
分析元素含量分形曲線特征,我們可以得出:①Mo、Au、Ag元素,尤其是Mo元素,其含量呈具多標度分形特征,表明其經歷兩次或兩次以上的次生富集成暈作用和過程,指示本區極有可能存在鉬礦(化)體;②Pb、Zn元素的分維值D2相對較小,分布范圍較廣,容易經過次生作用而富集形成實質性的地球化學異常,這對于在本區尋找金、鉛、鋅礦提供了一定的有用線索;③Cu元素的分維值D2相對較大,表明該元素含量的空間差異性較小,分布比較均勻,不易經過次生作用而富集形成實質性的地球化學異常。
3.5.1 用算術“迭代法”計算背景值及異常下限
對本區土壤樣品Au、Ag、Cu、Pb、Zn、Mo等六種元素分別采用算術計算法、對數計算法計算背景值及異常下限,用“迭代法”對高于背景值加上或低于背景值減去三倍標準差的值予以剔除(表4)。
3.5.2 利用分形特征曲線確定背景值及異常下限
在元素含量分形分析中,分維值D1可以反映沒有受到礦化干擾的背景區的元素含量分布,分維值D2、D3、D4可以代表受到礦化影響的異常區元素含量異常特征。在元素含量分形特征圖上,分維值D1、D2的交點代表了由背景區向異常區過渡的拐點,該交點對應的含量值可以認為是元素的背景值。
基于此,在本區元素含量分形特征圖(圖5)上,可以確定出元素的背景值c0(表5)。再根據公式:T=c0+ 2σ計算各元素的異常下限(表4)。

圖5 元素含量分形特征圖Fig.5 Element content fractal graph in Majiaosi molybdenum deposit area

元素AuMoAgPbZnCu初始數據計算結果樣品數103061030610306103061030610306最大值175.0198.04355.01400.95911.0159.8最小值0.20.0919.05.822.63.7算術法平均值2.4741.33087.30954.006121.20124.446標準差4.0673.64087.51156.407124.2068.614變異系數1.6442.7371.0021.0441.0250.352對數法平均值0.2890.0211.8981.6492.0181.363標準差0.2580.2150.1650.2350.2080.150變異系數0.89210.0770.0870.1430.1030.110剔除野值后計算結果迭代次數554555樣品數101431009110134100381005410201最大值9.243.29206.0168.3335.060.7最小值0.390.3229.010.830.08.8算術法背景值2.1911.09581.71347.690109.31024.346標準差1.3130.45828.98625.69949.5717.766變異系數0.5990.4180.3550.5390.4530.319異常下限T4.1892.012148.6699.089208.45139.879對數法背景值0.2780.0051.8881.6292.0021.364標準差0.2300.1710.1420.1990.1750.141變異系數0.82631.2310.0750.1220.0870.104異常下限T5.4582.223149.371106.225224.36344.348分形法背景值C0 1.8620.83279.43341.68783.17625.119標準差σ1.3130.45828.98625.69949.5717.766異常下限T4.4881.748137.40593.085182.31840.651綜合確定異常下限T5215010020040
綜合上述三種方法計算的結果,結合本區土壤測量實際情況,取整后最終確定出本區各單元素異常下限(表4)。
依據確定出的異常下限,在本區共圈出單元素異常122個(表5、圖6),其中金異常18個、銀異常25個、鉬異常23個、鉛異常18個、鋅異常21個、銅異常17個。
3.7.1 Au元素土壤地球化學異常特征
金異常區Au元素含量為5.0×10-9~175.0×10-9,最高值為175.0×10-9,為背景值(2.191×10-9)的79.87倍,是異常下限(5×10-9)的35倍,顯示區內有較好的金礦找礦前景。

表5 單元素異常特征一覽表
注:注:Au、Ag元素含量單位為10-9,其余元素含量單位為10-6

圖6 馬角寺鉬礦區單元素土壤地球化學異常Fig.6 Soil geochemical anomalies (1∶10,000) of single elements in Majiaosi molybdenum deposit area
金異常集中連片分布于礦區東南,異常多為不規則的條帶狀、橢圓狀,長軸多呈近東西向展布,產于輞峪緊閉倒轉向斜南翼的中、上統寒武系地層中。
3.7.2 Mo、Ag、Pb、Zn元素土壤地球化學異常 特征
(1)鉬異常區Mo元素含量為2.0×10-6~198.0×10-6,最高值為198.0×10-6,是背景值(1.095×10-6)的180.82倍,是異常下限(2×10-6)的99倍。多數鉬異常面積大、強度高、分帶性良好、濃集中心較明顯,指示區內極可能有鉬的成礦作用。
(2)銀異常區Ag元素含量為150.0×10-9~4 355.0×10-9,最高值為4 355.0×10-9,是背景值(81.713×10-9)的53.30倍,是異常下限(150×10-9)的29.03倍。
(3)鉛異常區Pb元素含量100.0×10-6~1 400.9×10-6,最高值為1 400.9×10-6,是背景值(47.690×10-6)的29.38倍,是異常下限(100×10-6)的14倍。
(4)鋅異常區Zn元素含量為200.0×10-6~5 911.0×10-6,最高值為5 911.0×10-6,是背景值(109.310×10-6)的54.08倍,是異常下限(200×10-6)的29.56倍。
Mo、Ag、Pb、Zn各元素異常在空間位置上接近,套合性較好,多集中分布于本區中部震旦系與寒武系地層接觸帶附近。
3.7.3 Cu元素土壤地球化學異常特征
銅異常區Cu元素含量為40.0×10-6~159.8×10-6,最高值為159.8×10-6,僅是背景值(24.346×10-6)的6.56倍,只有異常下限(40×10-6)的4.0倍。本區銅異常不發育,規模均較小,分布比較零散,多分布于寒武系地層中,其余地層分布較少。
區內共圈出六個綜合異常,編號分別為Ht-1-甲、Ht-2-甲、Ht-3-甲、Ht-4-丙、Ht-5-丙、Ht-6-丙(圖7)。
3.8.1 Ht-1-甲、Ht-2-甲、Ht-3-甲等綜合異常特征
Ht-1-甲綜合異常位于本區中西部,呈北東東向近似帶狀展布,長約2.1 km,寬0.4 km~1.3 km。由3個鉬異常、4個銀異常、3個鉛異常、3個鋅鋅異常、1個金異常組成,元素含量最高值:鉬含量為133.0×10-6,銀含量為4 355.0×10-9、鉛含量為459×10-6、鋅含量為2 900×10-6。

圖7 馬角寺鉬礦區綜合異常及成礦預測圖Fig.7 Integrated anomaly and metallogenic prediction of molybdenum ore zone in Majiaosi molybdenum deposit area
Ht-2-甲綜合異常位于礦區中部,近似帶狀分布,長約3.5 km,寬0.4 km~1.3 km。由6個鉬異常、6個銀異常、2個鉛異常、4個鋅異常、1個金異常、2個銅異常組成,元素含量最高值:鉬含量為198.0×10-6、銀含量為2 112.0×10-9、鉛含量為893.3×10-6、鋅含量為5 911×10-6、銅含量為92.5×10-6。
Ht-3-甲綜合異常位于本區中東部,呈東西向帶狀分布,長約2.4 km,寬約0.3 km~1.2 km。由1個鉬異常、3個銀異常、2個鉛異常、2個鋅異常、1個金異常組成,元素含量最高值:鉬含量為51.7×10-6,銀含量為41 921×10-9,鉛含量為283.2×10-6,鋅含量為1 208×10-6,金含量為133×10-9。
Ht-1-甲、Ht-2-甲、Ht-3-甲等3個綜合異常,集中分布于本區中部全堂溝腦-屈家-張家-寺通溝口-上岸溝腦一帶,異常強度較高,規模較大,其中鉬元素異常遠較其他元素異常的強度高、規模大。元素組合復雜,主要為Mo、Ag、Pb、Zn,其次為Au、Cu,其中Mo、Ag、Pb、Zn各單元素異常套合性好,異常疊加現象比較明顯,顯示具有尋找鉬多金屬礦產的地球化學條件。異常沿震旦系與寒武系地層接觸帶呈近東西向帶狀分布,局部呈北西向。
3.8.2 Ht-4-丙、Ht-5-丙、Ht-6-丙等綜合異常特征
Ht-4-丙、Ht-5-丙、Ht-6-丙等3個異常集中分布于本區南東部扁擔川—羊山底南部-輞峪溝口一帶,異常面積較大,但強度不高。元素組合較為復雜,主要為Au、Mo,其次為、Pb、Zn、Ag、Cu,主要元素異常套合性較差。異常產于中、上寒武系地層中,與北西向斷裂關系密切。
工作中對Ht-1-甲、Ht-2-甲等2個主要綜合異常利用巖石地球化學剖面測量開展查證,共布設5條地化剖面。
地化剖面比例尺為1∶2 000,布置在Ht-1-甲、Ht-2-甲兩個重點異常,基本垂直異常長軸方向布設剖面。按不同的地質體分別采取巖石樣,樣品間距一般為20 m,含礦帶及礦化體樣品間距加密為1 m~2 m。連續打塊取樣,分析項目:Au、Ag、Cu、Mo、 Pb、Zn。依據巖石元素含量值編制地化綜合剖面圖。
在Ht-1-甲異常區布設兩條地化剖面,經查證發現震旦系與寒武系地層接觸帶附近及含碳質片巖中鉬含量較高(圖8)。利用槽探工程揭露,表明震旦系與寒武系地層之間為斷層接觸,斷層帶為含鉬礦化構造帶,編號Q1(圖7)。通過槽探初步控制,并結合鉆孔(ZK06)深部驗證,在含鉬礦化構造帶中圈出K2、K3等兩個鉬礦體。
在Ht-2-甲異常區布設三條地化剖面,經查證發現含碳質片巖中鉬含量較高(圖9)。利用槽探工程揭露,發現震旦系與寒武系地層之間為斷層接觸,斷層帶為含鉬礦化構造帶,即Q1(圖7)。經槽探初步控制,并結合鉆孔(ZK03、ZK04、ZK07)深部驗證,在該含鉬礦化構造帶中圈出K1、K4兩個鉬礦體。

圖8 馬角寺鉬礦區南溝2號巖石地化剖面圖Fig.8 Geochemistry profile of nangou No.2 rock in Majiaosi molybdenum deposit area

圖9 馬角寺鉬礦區竹溝對面梁4號巖石地化剖面圖Fig.9 Geochemistry profile of No.4 rock at beam of opposite side in Zhugou in Majiaosi molybdenum deposit
Ht-3-甲綜合異常區地表覆蓋較厚,未能進行查證,推斷其應為Q1含鉬礦化構造帶引起的異常。分布于本區南東部的Ht-4-丙、Ht-5-丙、Ht-6-丙等3個綜合異常未進行查證,性質不明。
依據所圈綜合異常特征,結合本區成礦地質背景、礦化信息,在本區圈出二級找礦靶區及三級成礦靶區各1個(圖7)。
該找礦靶區位于本區中部木子溝及以東地區,呈規則的四邊形,面積約1.2 km2。包含Ht-3-甲綜合異常,Mo、Ag、Pb、Zn元素異常面積較大,空間位置上套合較好,元素相關性較強,Au、Cu元素異常亦有分布。
找礦靶區位于東西向輞峪緊閉倒轉向斜的北翼,多數地段地層倒轉,東西向逆沖斷層及北東、北西向平移斷裂極為發育,且多期活動。區內雖無大規模巖漿巖出露,但斷裂帶中常見多期次的石英脈等脈巖產出,表明靶區內曾有較為強烈的熱液活動。震旦系和寒武系地層為一近東西向逆沖斷層接觸,在靶區推斷Q1含鉬礦化構造帶賦存于該斷層接觸帶中,是尋找Mo礦床的有利靶區。
該找礦靶區位于本區扁擔川—輞峪溝一帶,呈規則四邊形,面積約6.5 km2。包含Ht-4-丙、Ht-5-丙、Ht-6-丙等3個綜合異常,Au、Mo、Pb、Zn元素異常面積較大、強度高,其他單元素異常亦有分布,Au、Pb、Zn元素異常套合性較好。
靶區內主要出露寒武系灰巖、白云質灰巖、白云巖,次為奧陶—志留系陶灣群廟灣組千枚巖。靶區位于東西向輞峪緊閉倒轉向斜的南翼,北西向斷裂較為發育。斷裂帶內常見石英脈等脈巖,是尋找Au、Mo礦產的較有利靶區。
(1)區內1∶50 000水系沉積物異常分布面積較大,組分復雜,主要元素組合為Au、Pb、Zn、Mo,其次為Cu,成礦元素為Au、Mo,元素間套合性好,區域成礦地質條件良好,表明該區具有找到鉬、金多金屬礦產的可能。
(2)1∶10 000土壤測量在區內圈出5個綜合異常,已在其中兩個綜合異常中找到鉬礦床,其余綜合異常屬重要的找礦靶區,可以安排進一步的查證工作。
(3)本次工作說明1∶10 000土壤測量在小秦嶺地區是尋找鉬、金多金屬礦的可行方法,應大力推廣使用。
致謝
作者在編寫過程中得到陜西省地礦局第六地質隊白和總工程師、夏元鵬高級工程師大力支持和幫助,在此表示感謝;在這里還要感謝項目組所有的工作人員。
參考文獻:
[1] 潘桂棠,肖慶輝,陸松年,等.中國大地構造單元劃分[J].中國地質, 2009,36(1):1-15.
[2] 宋小文,侯滿堂,陳如意.陜西省成礦區(帶)的劃分[J].西北地質,2004,37(3):29-42.
[3] 蒙勇,文件生,呂宇明,等.廣西貴港銀山嶺鉛鋅礦土壤地球化學找礦研究[J].現代地質,2011, 25(5):1013-1020.
[4] 樸壽成,劉樹田.吉林小石人金礦地球化學異常特征及成礦預測[J].地質與勘探,2003,39(2):26-29.
[5] 孟憲偉,張曉華.多標度分形與地球化學場分解[J].地質與勘探,1996,32(4):47-49.
[6] 成秋明.多維分形理論和地球化學元素分布規律[J].地球科學,2000,25(3):311-318.
[7] 成秋明.多重分形與地質統計學方法用于勘查地球化學異常空間結構和奇異性分析[J].地球科學,2001,26(2):161-166.
[8] 申維.分形混沌與礦產預測[M].北京:地質出版社, 2002.
[9] 魏民,趙鵬大,劉紅光,等.中國巖金礦床品位-噸位模型研究[J].地球科學-中國地質大學學報,2001,26(2):176-178.
[10] 王慶飛,萬麗,劉學飛.典型構造蝕變巖型金礦遠景資源量數學模型與預測[J].礦床地質,2007,26(3):341-345.
[11] 程小久,翟裕生.不同類型容礦構造中礦體鉛鋅品位分形特征[J].礦床地質,1994,13(S1):113-115.
[12] 萬麗,王慶飛.成礦元素品位有序數據集自仿射分形方法應用性評價[J].地球科學進展,2007,22(4):357-360.
[13] 鄧軍,王慶飛,楊立強,等.膠東西北部金熱液成礦系統內部結構解析[J].地球科學,2005,30(1):102-108.
[14] 謝和平.分形-巖石力學導論[M].北京:科學出版社,1996.
[15] 韓東昱,龔慶杰,向運川.區域化探數據處理的幾種分形方法[J].地質通報,2004,23(7): 714-719.
[16] 曾慶棟,沈遠超,劉鐵兵,等.山東平邑卓家莊金礦金品味的分形結構特征及意義[J].地質與勘探,2000,36(3):43-46.