潘勇，彭光雄，劉德波

?

基于Meta分析的不同化探異常找礦有效性定量評價

潘勇1, 2, 3，彭光雄1, 2，劉德波1, 2

(1. 中南大學(xué) 有色金屬成礦預(yù)測教育部重點(diǎn)實(shí)驗(yàn)室，湖南 長沙，410083;2. 中南大學(xué) 地球科學(xué)與信息物理學(xué)院，湖南 長沙，410083；3. 廈門大學(xué) 嘉庚學(xué)院，福建 廈門，363105)

利用青海雜多地區(qū)的1:5萬水系沉積物測量資料，分別采用常規(guī)統(tǒng)計法和多重分形法獲得區(qū)內(nèi)Pb和Zn元素的單元素地球化學(xué)異常。引入循證醫(yī)學(xué)領(lǐng)域廣泛使用的Meta分析方法對上述異常的找礦有效性進(jìn)行定量評價。通過實(shí)驗(yàn)對用于數(shù)據(jù)提取的多區(qū)域法和多級網(wǎng)格的特點(diǎn)及其適用條件進(jìn)行分析和總結(jié)。研究結(jié)果表明：利用Meta分析提取地質(zhì)數(shù)據(jù)時應(yīng)優(yōu)先采用多區(qū)域法，當(dāng)檢驗(yàn)樣本數(shù)量難以達(dá)到Meta分析的要求時方可采用多級網(wǎng)格法；基于Meta分析的找礦有效性評價不但能利用優(yōu)勢比O對找礦有效性的差異進(jìn)行定量描述，而且可利用假設(shè)檢驗(yàn)對有效性比較的結(jié)果是否具有統(tǒng)計學(xué)意義進(jìn)行判定；利用Meta分析進(jìn)行多種勘探方法的找礦有效性定量評價能為找礦預(yù)測提供關(guān)鍵、高效的找礦參數(shù)，具有較大的找礦應(yīng)用潛力。

Meta分析；化探異常；找礦有效性；定量評價；網(wǎng)格尺度

由于成礦規(guī)律及找礦預(yù)測中地質(zhì)變量的復(fù)雜性和不確定性，要準(zhǔn)確查明深部地質(zhì)體及隱伏礦體的空間位置及其規(guī)模和產(chǎn)狀往往非常困難[1?3]。目前，找礦工作已由地表淺層礦床轉(zhuǎn)向地下深部隱伏礦床，因此，對礦產(chǎn)勘查技術(shù)也提出了新的要求[4]。為適應(yīng)深部隱伏礦找礦工作的需要，地質(zhì)、物探、化探、遙感及鉆探等探礦技術(shù)和方法也得到不斷創(chuàng)新和發(fā)展。新型的礦產(chǎn)勘查技術(shù)不斷涌現(xiàn)，在深部找礦工作中發(fā)揮著重要作用。由于不同勘探方法對于不同的礦床類型、礦體特征及埋藏深度往往表現(xiàn)出不同的探測效果。如何準(zhǔn)確、高效和經(jīng)濟(jì)地獲取地下深部的有效找礦信息，并進(jìn)行地質(zhì)、物探、化探、遙感等多元找礦信息的優(yōu)選和集成已成為當(dāng)前找礦勘探中面臨的熱點(diǎn)和難點(diǎn)問題。不同找礦方法的定量評價技術(shù)是解決上述問題的重要途徑，近年來這方面的研究工作相繼得到開展，如：楊牧等[5]運(yùn)用因子分析和層次分析法對安徽銅陵鳳凰山銅礦構(gòu)造地球化學(xué)、γ能譜及可控源音頻大地電磁法等多種探測技術(shù)的找礦有效性進(jìn)行了評價；彭省臨等[2]將循證醫(yī)學(xué)領(lǐng)域廣泛使用的Meta分析方法對CSAMT，TEM和EH4等物探方法探測的二維反演剖面有效性進(jìn)行了定量對比；陳進(jìn)等[6]基于經(jīng)典的統(tǒng)計學(xué)理論提出了一種基于剖面對比的有效性定量評價方法。由于Meta分析是一種定量綜合研究的方法，在同類研究的樣本比較分析方面具有獨(dú)特的優(yōu)勢[7?9]，因而，在找礦方法有效性比較中比其他方法更加適宜。上述研究均僅對地下地質(zhì)體的二維反演剖面進(jìn)行了有效性比較，尚沒有對地表二維平面上的地質(zhì)、物探、化探、遙感異常開展研究。地表二維平面上的地質(zhì)、物探、化探、遙感異常在找礦過程中一般處于大區(qū)域掃面階段，對于不同探測方法有效性的研究和比較對后續(xù)遠(yuǎn)景靶區(qū)的圈定及地質(zhì)工程的布置均具有重要指導(dǎo)意義，因此，對不同方法獲得的地表二維平面上的找礦異常開展有效性評價顯得尤為重要。統(tǒng)計方法是化探數(shù)據(jù)處理中常用的方法，然而，嚴(yán)格地說，它并不具備刻畫異常值的功能[10]。多維分形方法既考慮了數(shù)據(jù)的頻率分布特征和空間信息屬性，同時又符合地球化學(xué)數(shù)據(jù)場具有標(biāo)度不變性和自相似性的特點(diǎn), 從而使得其在化探異常提取方面更具優(yōu)勢[11?13]。本文利用Meta分析對統(tǒng)計法和多重分形法獲得的地球化學(xué)異常開展找礦有效性研究。

1 數(shù)據(jù)來源及預(yù)處理

三江北段成礦帶位于青藏高原東北緣，是中國著名的三江銅多金屬巨型成礦帶的北延部分。研究區(qū)位于三江北段成礦帶的青海雜多地區(qū)，其TM7波段的遙感灰度影像及地球化學(xué)綜合異常如圖1所示。區(qū)內(nèi)找礦主攻的礦種為銅、金、銀、鉛、鋅，區(qū)內(nèi)的典型礦床有納日貢瑪銅鉬礦、然者涌、東莫扎抓及莫海拉亨鉛鋅礦等。地球化學(xué)數(shù)據(jù)來源于12幅1:5萬水系沉積物測量資料。水系沉積物采樣點(diǎn)空間分布比較均勻，平均采樣密度約為4點(diǎn)/km2，共測定Au，Ag，Cu，Pb，Zn，As和Sb等12種成礦元素的含量。

圖1 研究區(qū)遙感影像及地球化學(xué)綜合異常分布圖

對地球化學(xué)元素含量先取對數(shù)使其近正態(tài)分布，然后進(jìn)行統(tǒng)計。根據(jù)對數(shù)平均值和2倍對數(shù)標(biāo)準(zhǔn)離差，得到相應(yīng)元素的對數(shù)異常下限，再轉(zhuǎn)換成元素的實(shí)際異常下限值[13?14]。統(tǒng)計法得到的Pb和Zn元素的異常下限值(質(zhì)量分?jǐn)?shù)，下同)分別為74×10?6和215×10?6。多重分形采用含量?面積法分別計算大于某分組段元素對應(yīng)的面積。然后在雙對數(shù)坐標(biāo)中，采用最小二乘法求出每段曲線斜率的估計量，即得到分維數(shù)。根據(jù)多重分形特征，一般可取第二直線段(分維數(shù)為2)與第三直線段(分維數(shù)為3)的臨界點(diǎn)作為地球化學(xué)元素的異常下限值[14]。采用多重分形法得到的Pb和Zn元素的異常下限值分別為60×10?6和240×10?6。利用異常下限進(jìn)行數(shù)據(jù)分割，可分別獲得統(tǒng)計法和多重分形法提取的單元素地球化學(xué)異常。局部地區(qū)的Pb元素異常如圖2所示。利用區(qū)內(nèi)的已知礦點(diǎn)，可對上述2種方法獲得的地球化學(xué)異常的找礦有效性進(jìn)行比較分析。

(a) 統(tǒng)計法Pb異常；(b) 分形法Pb異常

圖2 統(tǒng)計法和多重分形法獲取的Pb元素地球化學(xué)異常

Fig. 2 Geochemical anomalies of Pb element derived from statistics and multifractal methods

2 Meta分析及數(shù)據(jù)提取

2.1 Meta分析

Meta分析(Meta-analysis)是對同一主題的大量獨(dú)立實(shí)驗(yàn)結(jié)果進(jìn)行綜合的統(tǒng)計方法，在同類研究的樣本比較分析方面具有獨(dú)特的優(yōu)勢[7]。從統(tǒng)計的角度，該方法起到了增大樣本量、提高檢驗(yàn)效能的作用。尤其是當(dāng)多個研究結(jié)果不一致或研究結(jié)果都沒有統(tǒng)計學(xué)意義時，采用Meta分析可得到更加接近真實(shí)情況的綜合分析結(jié)果。因此，利用Meta分析進(jìn)行找礦方法的有效性評價具有重要意義。

Meta分析是一種對若干獨(dú)立研究的統(tǒng)計結(jié)果進(jìn)行綜合、分析的統(tǒng)計方法[9]。20世紀(jì)80年代初期，Meta 分析開始在醫(yī)學(xué)界得到廣泛應(yīng)用和迅猛發(fā)展。20世紀(jì)90 年代，生態(tài)學(xué)家們又把它引入生態(tài)領(lǐng)域，引起生態(tài)學(xué)界的高度重視[15?17]。2009年彭省臨等[2]將Meta分析應(yīng)用于不同物探的測深剖面定量比較。本文采用Meta分析方法對不同方法獲得的化探異常的找礦有效性進(jìn)行定量評價。

2.2 數(shù)據(jù)提取的方法

Meta分析的統(tǒng)計學(xué)過程包括合并統(tǒng)計量的選擇、異質(zhì)性檢驗(yàn)和合并統(tǒng)計量的假設(shè)檢驗(yàn)3部分[9]。利用Meta分析進(jìn)行對比分析和評價可概況為以下6步： 1) 確定研究主題和方案；2) 廣泛收集文獻(xiàn)資料；3) 資料質(zhì)量評價；4) 數(shù)據(jù)提??；5) 統(tǒng)計學(xué)處理；6) 結(jié)果解釋及評價。其中數(shù)據(jù)提取是最關(guān)鍵的環(huán)節(jié)。Meta分析的本意主要是對以往文獻(xiàn)的研究結(jié)果進(jìn)行定量比較分析[18]。而相對于經(jīng)典的Meta分析，地質(zhì)勘探則無法像醫(yī)學(xué)研究那樣直接從文獻(xiàn)中獲取所需數(shù)據(jù)。若要利用Meta分析處理地質(zhì)勘探數(shù)據(jù)，其前提是將地質(zhì)調(diào)查獲得的各種數(shù)據(jù)經(jīng)提取轉(zhuǎn)化為Meta分析所能接收的定量信息。網(wǎng)格單元法為統(tǒng)計學(xué)進(jìn)入地質(zhì)勘探與成礦預(yù)測提供了有效的途徑[2, 6]。彭省臨等[2]利用網(wǎng)格單元法較好地解決了Meta分析中地質(zhì)數(shù)據(jù)提取的問題。針對地表二維平面上的各種地質(zhì)、物探、化探、遙感異常，基于規(guī)則網(wǎng)格可利用地質(zhì)資料或野外觀測得到的礦(化)點(diǎn)來進(jìn)行找礦有效性比較和評價，如圖3所示。實(shí)現(xiàn)地質(zhì)勘探資料Meta分析數(shù)據(jù)提取的方法可分為多區(qū)域法和多級網(wǎng)格法2種。多區(qū)域法是指對多個不同的成礦帶或同一個成礦帶的不同圖幅或者區(qū)域，分別統(tǒng)計有礦(化)點(diǎn)落入的找礦異常所在的網(wǎng)格數(shù)量和全部異常所在的網(wǎng)格總數(shù)的方法，此時，每個成礦區(qū)(區(qū)域)相當(dāng)于Meta分析中的1篇文獻(xiàn)，每個網(wǎng)格則相當(dāng)于研究文獻(xiàn)中的1個實(shí)驗(yàn)對象。多級網(wǎng)格法是指在同一個區(qū)域內(nèi)采用0.5 km×0.5 km，1.0 km×1.0 km， 2.0 km×2.0 km等多種不同網(wǎng)格尺度，分別統(tǒng)計有礦(化)點(diǎn)落入的找礦異常所在的網(wǎng)格數(shù)量和全部異常所在的網(wǎng)格總數(shù)的方法，此時，每一級網(wǎng)格相當(dāng)于Meta分析中的1篇文獻(xiàn)，每個網(wǎng)格則相當(dāng)于研究文獻(xiàn)中的1個實(shí)驗(yàn)對象。在地質(zhì)數(shù)據(jù)和檢驗(yàn)樣本數(shù)量較多、空間分布比較分散且均勻的情況下使用多區(qū)域法較適宜；反之，則不利于地質(zhì)資料的提取量化轉(zhuǎn)換，為了達(dá)到Meta分析的要求，可采用多級網(wǎng)格法。下面利用實(shí)例對多區(qū)域法和多級網(wǎng)格法的性能進(jìn)行比較分析。

圖3 基于規(guī)則網(wǎng)格不同地表找礦異常有效性定量比較示意圖

Fig. 3 Quantitative evaluation diagram of prospecting effectiveness for different surface anomalies based on grid

3 數(shù)據(jù)提取的結(jié)果

網(wǎng)格單元大小的選擇既要綜合考慮不同找礦方法的勘探精度、地質(zhì)條件及已有的檢驗(yàn)數(shù)據(jù)(如礦點(diǎn))的空間分布情況，也要保障各類地質(zhì)體或找礦異常所占據(jù)的網(wǎng)格數(shù)量在進(jìn)行有效性評價時具有統(tǒng)計學(xué)意義。區(qū)內(nèi)12個1:5萬標(biāo)準(zhǔn)圖幅中僅有8個圖幅有礦點(diǎn)資料，利用這個8個圖幅的已知礦點(diǎn)資料，分別采用0.5 km×0.5 km，1.0 km×1.0 km和2.0 km×2.0 km這3種網(wǎng)格尺度對統(tǒng)計法和多重分形法獲得的Pb和Zn異常，統(tǒng)計有礦點(diǎn)落入的異常網(wǎng)格數(shù)和有異常落入的網(wǎng)格總數(shù)，進(jìn)行找礦有效性測量指標(biāo)的提取，統(tǒng)計結(jié)果如表1和表2所示。

表1 不同方法獲得的Pb元素異常與已知礦點(diǎn)的疊合情況

表2 不同方法獲得的Zn元素異常與已知礦點(diǎn)的疊合情況

為了定量描述某種找礦方法的有效程度，通?？梢岳糜行颖竞涂倶颖镜谋戎祦砗饬?。本文定義找礦異常的網(wǎng)格有效率為：在一定網(wǎng)格尺度下，有已知礦點(diǎn)存在的異常網(wǎng)格數(shù)量與異常所占據(jù)的總網(wǎng)格數(shù)量的比值。利用如表1和表2中的數(shù)據(jù)，可以得到不同網(wǎng)格尺度下統(tǒng)計法和多重分形法獲得的Pb和Zn元素異常的網(wǎng)格有效率，如表3所示。由表3可知：在1 km×1 km的網(wǎng)格尺度時，統(tǒng)計法和多重分形法的Pb元素異常網(wǎng)格有效率分別為13.46%和21.72%，說明多重分形法的找礦有效性比統(tǒng)計法的高。同理，對于相同尺度下的Pb和Zn異常，多重分形法的找礦有效性均比統(tǒng)計法的高，但該結(jié)果無法進(jìn)行統(tǒng)計學(xué)的顯著性檢驗(yàn)。表3還顯示：隨著網(wǎng)格尺度的增加，Pb和Zn元素異常的網(wǎng)格有效率也隨之增大，說明網(wǎng)格有效率具有不穩(wěn)定性，因此，這種常規(guī)的有效性比較方法對于找礦方法的定量評價尚存在缺陷。而Meta分析為解決上述不足提供了有效的技術(shù)手段。

表3 不同網(wǎng)格尺度的Pb和Zn元素異常網(wǎng)格有效率

4 結(jié)果分析

Review Manager是國際Cochrane協(xié)作網(wǎng)開發(fā)的一套用于循證醫(yī)學(xué)Meta分析的公益軟件。本文利用Review Manager5.1版本來實(shí)現(xiàn)找礦有效性定量比較，利用表1和表2中數(shù)據(jù)可得到Meta分析森林圖和表4所示的分析結(jié)果。本文僅列出多區(qū)域法中Pb和Zn異常找礦有效性比較的部分森林圖，以及多級網(wǎng)格法中Pb+Zn元素異常的森林圖(如圖4所示)。

表4 多區(qū)域法和多級網(wǎng)格法Meta分析結(jié)果

注：2為統(tǒng)計量；為統(tǒng)計量的概率。

(a) 基于2.0 km×2.0 km網(wǎng)格的Pb異常；(b) 基于2.0 km×2.0 km網(wǎng)格的Zn異常；(c) 基于多級網(wǎng)格的化探異常

圖4 找礦有效性比較Meta分析森林圖

Fig. 4 Forest graphs of meta-analysis for prospecting effectiveness

當(dāng)異質(zhì)性檢驗(yàn)的概率＞0.05時，可認(rèn)為多個獨(dú)立研究具有同質(zhì)性，可選擇固定效應(yīng)模型(fixed effect model)計算其合并統(tǒng)計量。文中全部分析案例均具有良好的同質(zhì)性，因此，均采用固定效應(yīng)模型進(jìn)行Meta分析。從表4可見：除多區(qū)域法中2 km×2 km網(wǎng)格尺度Zn異常的總體效應(yīng)檢驗(yàn)的概率=0.17(大于0.05)，統(tǒng)計學(xué)上不具有顯著性外，其他的有效性比較顯著性明顯，都具有統(tǒng)計學(xué)意義。森林圖最下方的菱形代表多個隨機(jī)對照試驗(yàn)(randomized controlled trial, RCT)綜合結(jié)果。垂直線(代表優(yōu)勢比O=1)將圖分為左、右兩半，用于判斷結(jié)果差異有無統(tǒng)計學(xué)意義。當(dāng)菱形與垂直線相交時則表明該隨機(jī)對照試RCT 中不同方法之間差異無統(tǒng)計學(xué)意義。菱形的位置及其含義為：對于不利結(jié)局如殘廢事件及死亡等，當(dāng)菱形完全在垂直線左側(cè)時表示治療組更加有效，當(dāng)菱形完全在右側(cè)時表示對照組更有效；對于有利結(jié)局如緩解及治愈等情況，菱形的位置及其含義則恰恰相反；菱形偏離垂直線的距離越大，說明2種方案的有效性差異越明顯[18]。

Meta分析不但能有效比較2種方法的有效性，準(zhǔn)確得出“有差別”或“無差別”的結(jié)論，而且通過效應(yīng)尺度能對這種差異進(jìn)行定量描述。反映多個獨(dú)立研究綜合效應(yīng)的指標(biāo)主要包括優(yōu)勢比(O)、相對危險度(R)或率差(D)等合并統(tǒng)計量，同類型的效應(yīng)尺度之間可進(jìn)行數(shù)量對比[9]。Meta分析得到的優(yōu)勢比O越大，森林圖中菱形偏離垂直線的距離也就越大，也說明2種方案的有效性差異越明顯。由表3可知：在0.5 km×0.5 km，1.0 km×1.0 km和2.0 km×2.0 km的網(wǎng)格尺度下，多重分形法比統(tǒng)計法的Pb元素異常網(wǎng)格有效率分別高2.43%，8.26%和14.87%，相應(yīng)的Zn元素異常網(wǎng)格有效率則分別高1.44%，5.68%和7.94%?？梢婋S著網(wǎng)格尺度的變化，網(wǎng)格有效率也隨之發(fā)生較大變化，說明利用網(wǎng)格有效率作為衡量有效性的指標(biāo)具有不穩(wěn)定性。雖然在同一網(wǎng)格尺度下網(wǎng)格有效率具有比較意義，但不同的網(wǎng)格尺度之間沒有比較意義。由表4可知：基于多區(qū)域法在0.5 km×0.5 km，1.0 km×1.0 km和2.0 km×2.0 km的網(wǎng)格尺度下, Pb元素異常多重分形法與統(tǒng)計法的Meta分析優(yōu)勢比O分別為1.91，1.89和2.06，Zn元素異常相應(yīng)的優(yōu)勢比O則分別為1.29，1.38和1.56(差異不顯著)，可見網(wǎng)格尺度的變化對于優(yōu)勢比O的影響不是很明顯，也就是說優(yōu)勢比O作為衡量有效性的指標(biāo)具有良好的穩(wěn)定性。從表4可見：多級網(wǎng)格法獲得的Pb和Zn元素異常優(yōu)勢比O分別為1.81和1.36，與前述不同網(wǎng)格尺度下多區(qū)域法得到的優(yōu)勢比O差別不大。然而，無論是多區(qū)域法還是多級網(wǎng)格法，在利用網(wǎng)格單元對提取的數(shù)據(jù)進(jìn)行Meta分析時，網(wǎng)格尺度也應(yīng)當(dāng)處于一個適宜的范圍內(nèi)，超出適宜范圍將會導(dǎo)致Meta分析的結(jié)果失去統(tǒng)計學(xué)意義。多級網(wǎng)格法由于采用的網(wǎng)格尺度有限，從而造成試驗(yàn)對象的數(shù)量也受到較大限制，并且采用不同尺度的網(wǎng)格對樣本統(tǒng)計是否存在相關(guān)性(即不同網(wǎng)格尺度下的各個試驗(yàn)之間是否是相對獨(dú)立的)目前也難以確定，還需要進(jìn)行進(jìn)一步研究。因此，Meta分析中進(jìn)行數(shù)據(jù)提取時，應(yīng)優(yōu)先采用多區(qū)域法，當(dāng)檢驗(yàn)樣本數(shù)量難以達(dá)到Meta分析的要求時方可采用多級網(wǎng)格法，或者將多區(qū)域與多級網(wǎng)格混合使用。

5 結(jié)論

1) 對于常規(guī)的網(wǎng)格有效率找礦有效性評價方法，隨著網(wǎng)格尺度的增加，所得到的網(wǎng)格有效率出現(xiàn)大幅度變化，存在不穩(wěn)定缺陷，并且還無法進(jìn)行統(tǒng)計學(xué)的顯著性檢驗(yàn)。而Meta分析為解決上述不足提供了有效的技術(shù)手段。Meta分析不但能定性比較2種方法的有效性，準(zhǔn)確得出“有差別”或“無差別”的結(jié)論，而且通過優(yōu)勢比O可對這種差異進(jìn)行定量描述，優(yōu)勢比O越大，說明2種方案的有效性差異越明顯。Meta分析的優(yōu)勢比O具有不隨網(wǎng)格尺度的變化而變化的穩(wěn)定特性，并且還可利用合并統(tǒng)計量的假設(shè)檢驗(yàn)對有效性比較的結(jié)果是否具有統(tǒng)計學(xué)意義進(jìn)行判定。

2) 在Meta分析中利用網(wǎng)格單元進(jìn)行數(shù)據(jù)提取時，網(wǎng)格尺度也應(yīng)當(dāng)限定于一個適宜的范圍內(nèi)，超出適宜范圍將會導(dǎo)致Meta分析的結(jié)果失去統(tǒng)計學(xué)意義。適宜的網(wǎng)格尺度需要根據(jù)已知地質(zhì)資料的比例尺以及數(shù)據(jù)豐富程度，根據(jù)經(jīng)驗(yàn)和多次試驗(yàn)進(jìn)行選擇，以確保Meta分析的結(jié)果具有統(tǒng)計學(xué)意義。多級網(wǎng)格法由于采用的網(wǎng)格尺度比較有限，從而造成試驗(yàn)對象的數(shù)量也受到較大限制，并且不同網(wǎng)格尺度下的各個試驗(yàn)之間是否相對獨(dú)立目前也難以確定，還需要更進(jìn)一步研究。因此，采用Meta分析進(jìn)行數(shù)據(jù)提取時，應(yīng)優(yōu)先采用多區(qū)域法；當(dāng)檢驗(yàn)樣本數(shù)量難以達(dá)到Meta分析的要求時方可采用多級網(wǎng)格法，或者將多區(qū)域與多級網(wǎng)格混合使用。

3) 利用Meta分析方法進(jìn)行多種勘探方法的找礦有效性定量評價，能為找礦預(yù)測提供關(guān)鍵高效的找礦參數(shù)，為深層次挖掘地質(zhì)、物探、化探、遙感等數(shù)據(jù)信息的潛力提供新的途徑和科學(xué)手段，具有良好的找礦應(yīng)用潛力。

致謝：

感謝青海省地質(zhì)調(diào)查院為本文提供的化探及礦點(diǎn)數(shù)據(jù)。

[1] Park N W, Chi K H, Kwon B D. Accounting for spatial patterns of multiple geological data sets in geological thematic mapping using GIS-based spatial analysis[J]. Environmental Geology, 2007, 51(2): 1147?1155.

[2] 彭省臨, 王穎, 邵擁軍, 等. Meta分析在深部找礦預(yù)測研究中的應(yīng)用[J]. 高校地質(zhì)學(xué)報, 2009, 15(4): 460?469.
PENG Shenglin, WANG Ying, SHAO Yongjun, et al. Application of Meta-analysis to deep prospecting predictions: an example from the Tongshan copper deposit, Tongling, Anhui Province[J]. Geological Journal of China Universities, 2009, 15(4): 460?469.

[3] 彭光雄, 高光明. 新疆阿祖拜礦田偉晶巖型礦床遙感找礦綜合信息研究[J]. 大地構(gòu)造與成礦學(xué), 2013, 37(1): 109?118.
PENG Guangxiong, GAO Guangming. Remote sensing prospecting of pegmatite deposits in the Azubai Region, Xinjiang[J]. Geotectonica Et Metallogenia, 2013, 37(1): 109?118.

[4] 彭省臨, 邵擁軍, 張建東. 金屬礦山隱伏礦找礦預(yù)測理論與方法[J]. 地質(zhì)通報, 2011, 30(4): 538?543.
PENG Shenglin, SHAO Yongjun, ZHANG Jiandong. Prospecting prediction theory and method about concealed ore-deposit in metal mine[J]. Geological Bulletin of China, 2011, 30(4): 538?543.

[5] 楊牧, 彭省臨, 邵擁軍, 等. 隱伏礦定位預(yù)測中探查技術(shù)有效性評價: 以安徽銅陵鳳凰山銅礦南區(qū)為例[J]. 礦床地質(zhì), 2006, 25(S1): 495?499.
YANG Mu, PENG Shenglin, SHAO Yongjun, et al. The effectiveness of exploration methods in location prediction of concealed ore deposits: A case study of the Fenghuangshan Cu ore deposit in Tongling, Anhui Province[J]. Mineral Deposits, 2006, 25(S1): 495?499.

[6] 陳進(jìn), 毛先成, 彭省臨, 等. 物化探技術(shù)有效性定量評價方法研究: 以鳳凰山銅礦CSAMT 法為例[J]. 大地構(gòu)造與成礦學(xué), 2010, 34(2): 263?269.
CHEN Jin, MAO Xiancheng, PENG Shenglin, et al. Quantitative evaluation on availability through geophysical and geochemical prospecting techniques: A case study on CSAMT approach in Fenghuangshan Copper Deposit[J]. Geotectonica Et Metallogenia, 2010, 34(2): 263?269.

[7] Boris T, Peter H, Mark J, et al. Preferences for European agrarian landscapes: A meta-analysis of case studies[J]. Landscape and Urban Plannin, 2014, 132(6): 89?101.

[8] Akira S, Mori R. Retention forestry as a major paradigm for safeguarding forest biodiversity in productive landscapes: A global meta-analysis[J]. Biological Conservation, 2014, 175(7): 65?73.

[9] 劉關(guān)鍵, 吳泰相, 康德英. Meta分析中的統(tǒng)計學(xué)過程[J]. 中國臨床康復(fù), 2003, 7(4): 538?539.
LIU Guanjian, WU Taixiang, KANG Deying. The statistical process of meta-analysis[J]. Modern Rehabilitation, 2003, 7(4): 538?539.

[10] 成秋明. 多維分形理論和地球化學(xué)元素分布規(guī)律[J]. 地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報), 2000, 25(3): 311?319.
CHENG Qiuming. Multifractal theory and geochemical element distribution pattern[J]. Earth Science (Journal of China University of Geosciences), 2000, 25(3): 311?319.

[11] CHENG Qiuming. The perimeter-area fractal model and its application to geology[J]. Mathematical Geology, 1995, 27(1): 69?82.

[12] Naimark A A. Coarsely discrete fractal structure of geological medium and challenges in tectonophysical modeling[J]. Moscow University Geology Bulletin, 2009, 64(5): 273?280.

[13] 徐明鉆, 朱立新, 馬生明, 等. 多重分形模型在區(qū)域地球化學(xué)異常分析中的應(yīng)用探討[J]. 地球?qū)W報, 2010, 31(4): 611?619.
XU Mingzuan, ZHU Lixin, MA Shengming, et al. A tentative discussion on the application of multi-fractal models to the analysis of regional geochemical anomalies[J]. Acta Geoscientica Sinica, 2010, 31(4): 611?619.

[14] 鄒林, 彭省臨, 楊自安, 等. 青海阿爾茨托山地區(qū)地球化學(xué)異常場的多重分形研究[J]. 中國地質(zhì), 2004, 31(4): 436?442.
ZOU Lin, PENG Shenglin, YANG Zian, et al. Multifractal study of geochemical (anomaly) fields in the Aercituoshan area, Qinghai[J]. Chinese Geology, 2004, 31(4): 436?442.

[15] Tonhasca A, Byrne D N. The effect of crop diversification on herbivorous insects: A Meta-analysis approach[J]. Ecological Entomology, 1994, 19(1): 239?244.

[16] Curtis P. A meta-analysis of leaf gas exchange and nitrogen in trees grown under elevated carbon dioxide[J]. Plant, Cell and Environment, 1996, 19(3): 127?137.

[17] Ballantyne A P, Lavine M, Crowley T J. Meta-analysis of tropical surface temperatures during the Last Glacial Maximum[J]. Geophysical Research Letters, 2005, 26(4): 571?586.

[18] 黃魯眾, 楊騰飛, 黃晶晶, 等. 高溫作業(yè)對工人心電圖影響的Meta分析[J]. 中國循證醫(yī)學(xué)雜志, 2011, 11(2): 174?180.
HUANG Luzhong, YANG Tengfei, HUANG Jingjing, et al. ECG Changes in workers exposed to high-temperature: A Meta-analysis[J]. Chinese Journal of Evidence-based Medicine, 2011, 11(2): 174?180.

Quantitative evaluation on prospecting effectiveness of different geochemical anomalies based on Meta-analysis

PAN Yong1, 2, 3, PENG Guangxiong1, 2, LIU Debo1, 2

(1. Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education, Central South University, Changsha 410083, China;2. School of Geosciences and InfoPhysics, Central South University, Changsha 410083, China;3. College of Tan Kah Kee, Xiamen University, Xiamen 363105, China)

The survey data of stream sediment with the scale of 1:50 000 in Zaduo, Qinghai Province was used to obtain geochemical anomalies of Pb and Zn elements based on the methods of statistics and multifractal. Meta-analysis is a comprehensive statistical method used in evidence-based medicine, which was used to evaluate the prospecting effectiveness of different geochemical anomalies quantitatively. The characteristics and applicable conditions of multi-regional and multi-level grid methods used for data extraction were analyzed and summarized by experiment. The results show that when geological data are extracted based on Meta-analysis, multi-regional method should be used firstly. When the number of test samples is difficult to meet the requirements of the Meta-analysis, multi-level grid can be used as an alternative method. There are two superiorities of using Meta-analysis to evaluate the prospecting effectiveness.Odds ratio (O) of Meta-analysis can describe the differences of effectiveness quantitatively. In addition, hypothesis testing can be used to judge whether there is a statistically significant for the Meta-analysis result of effectiveness. Meta-analysis is an effective tool for effectiveness evaluation of geophysical and geochemical prospecting techniques. The evaluation results of prospecting effectiveness can provide the critical and efficient prospecting parameters for metallogenic prediction. This method has good application potential.

Meta-analysis; geochemical anomaly; prospecting effectiveness; quantitative evaluation; grid size

P632；P642

A

1672?7207(2015)01?0180?08

2014?01?13；

2014?03?24

中國博士后科學(xué)基金資助項(xiàng)目(2012T50832)；湖南省科技計劃項(xiàng)目(2012RS4047)；青海省重大科技專項(xiàng)(2010-J-A1) (Project(2012T50832) supported by China Postdoctoral Science Foundation; Project(2012RS4047) supported by Hunan Province Science and Technology Plan; Project(2010-J-A1) supported by Qinghai Province Key Project of Science and Technology)

彭光雄，博士，講師，從事遙感地質(zhì)與成礦預(yù)測研究；E-mail: pgxcsu@csu.edu.cn

10.11817/j.issn.1672?7207.2015.01.025

(編輯 陳燦華)