GIS技術(shù)在勘查中應(yīng)用第四章預(yù)測(cè)模型及其應(yīng)用

第四章礦產(chǎn)資源預(yù)測(cè)模型及其應(yīng)用中國(guó)地質(zhì)大學(xué)資源學(xué)院數(shù)據(jù)輸入數(shù)據(jù)管理模型建立空間分析數(shù)據(jù)輸出4.1GIS應(yīng)用于成礦預(yù)測(cè)的一般流程

▲GIS應(yīng)用于成礦預(yù)測(cè)，除了需要滿足一般GIS應(yīng)用的數(shù)據(jù)流程要求外，又需要滿足本學(xué)科特有的數(shù)據(jù)要求，其流程為：（1）地礦信息獲取和空間數(shù)據(jù)庫(kù)的建立（2）預(yù)測(cè)信息的提取----變量（3）成礦預(yù)測(cè)模型的建立（4）空間分析（2、3、4相互交叉驗(yàn)證）（5）解釋最終結(jié)果

4.1.1地礦信息獲取和空間數(shù)據(jù)庫(kù)的建立利用GIS提供的數(shù)據(jù)輸入平臺(tái)將地學(xué)資料分類(lèi)輸入。金屬礦產(chǎn)預(yù)測(cè)的多源地礦信息源一般包括：研究區(qū)地質(zhì)礦產(chǎn)圖，研究區(qū)地球物理數(shù)據(jù)或物探異常圖，研究區(qū)地球化學(xué)數(shù)據(jù)或化探異常圖，研究區(qū)遙感數(shù)據(jù)或遙感圖像及遙感解譯圖等。輸入的方法可采用數(shù)字化儀、掃描儀和鍵盤(pán)輸入。對(duì)輸入的各類(lèi)數(shù)據(jù)應(yīng)標(biāo)注其屬性及數(shù)據(jù)間的拓?fù)潢P(guān)系，并將輸入數(shù)據(jù)轉(zhuǎn)換為統(tǒng)一的某一種坐標(biāo)投影系統(tǒng)和有利于空間分析的數(shù)據(jù)格式（矢量或柵格格式），進(jìn)行誤差校正。從而建立起可供查詢(xún)、檢索及進(jìn)行各種分析的空間數(shù)據(jù)庫(kù)。4.1.2預(yù)測(cè)信息的提取根據(jù)所研究礦產(chǎn)模型的需要，從原始數(shù)據(jù)庫(kù)中提取有用的證據(jù)單元，或重新組合成基本的證據(jù)層。例如：利用物、化探資料產(chǎn)生地球物理和地球化學(xué)異常圖；提供有利于成礦的線元或面元形成單獨(dú)或組合的證據(jù)圖層；對(duì)礦點(diǎn)、礦床的產(chǎn)出位置進(jìn)行提取得到礦點(diǎn)、礦床分布圖等。4.1.3預(yù)測(cè)模型的建立證據(jù)權(quán)法信息量法BP神經(jīng)網(wǎng)絡(luò)ART神經(jīng)網(wǎng)絡(luò)特征分析法邏輯回歸模糊邏輯4.1.4空間分析以重新組合的證據(jù)層為基本單元，利用GIS強(qiáng)大的圖形分析功能，按照專(zhuān)家系統(tǒng)提供的知識(shí)和經(jīng)驗(yàn)，對(duì)各個(gè)圖形進(jìn)行操作分析。通過(guò)GIS的圖形顯示功能觀察中間結(jié)果，對(duì)比、分析、評(píng)價(jià)、修改初始模型和圖形分析方法，產(chǎn)生新的中間圖層，直到滿足礦產(chǎn)資源預(yù)測(cè)的要求為止。4.1.5解釋最終結(jié)果對(duì)多層地學(xué)證據(jù)進(jìn)行疊加，產(chǎn)生研究區(qū)礦產(chǎn)預(yù)測(cè)圖。即產(chǎn)生預(yù)測(cè)的最終成果。利用GIS產(chǎn)生的礦產(chǎn)預(yù)測(cè)圖是以不同色級(jí)標(biāo)示的一幅礦產(chǎn)資源空間可能產(chǎn)出的概率圖。根據(jù)其概率的不同可以圈定各級(jí)遠(yuǎn)景區(qū)。同時(shí)可以通過(guò)組織專(zhuān)家分析成果應(yīng)用的有效性，并可以修改或重新構(gòu)成基本證據(jù)的模型參數(shù)，產(chǎn)生新的礦產(chǎn)資源預(yù)測(cè)圖件。4.1.6礦產(chǎn)資源GIS評(píng)價(jià)系統(tǒng)基于GIS的礦產(chǎn)資源評(píng)價(jià)工作目前國(guó)內(nèi)成果：金屬礦產(chǎn)資源評(píng)價(jià)分析系統(tǒng)（MORPAS）（趙鵬大、胡光道等）礦產(chǎn)資源綜合信息礦產(chǎn)資源預(yù)測(cè)評(píng)價(jià)系統(tǒng)（王世稱(chēng)等）礦產(chǎn)資源綜合評(píng)價(jià)系統(tǒng)（MRAS）

（肖可炎等）礦產(chǎn)資源評(píng)價(jià)與預(yù)測(cè)工作一般流程地質(zhì)物探化探遙感礦產(chǎn)數(shù)據(jù)處理與分析成礦信息提取經(jīng)驗(yàn)?zāi)Ｐ捅O(jiān)督模型非監(jiān)督模型找礦有利度與靶區(qū)分類(lèi)踏勘驗(yàn)證GIS找礦標(biāo)志提取地質(zhì)背景，地質(zhì)異常分析與找礦標(biāo)志提取物探資料處理分析與找礦標(biāo)志提取化探資料處理分析與找礦標(biāo)志提取遙感資料處理分析與找礦標(biāo)志提取成礦信息綜合多源信息空間數(shù)據(jù)庫(kù)成果表達(dá)成果表達(dá)成果表達(dá)兩個(gè)核心概念

一、單元?jiǎng)澐謫卧獎(jiǎng)澐质沁M(jìn)行地質(zhì)異常提取、礦產(chǎn)資源綜合評(píng)價(jià)最基礎(chǔ)的環(huán)節(jié)，單元?jiǎng)澐值煤侠砼c否將對(duì)礦產(chǎn)預(yù)測(cè)結(jié)果產(chǎn)生直接的影響。它不僅是統(tǒng)計(jì)觀測(cè)的基本單位，也是提取與資源特征有密切關(guān)系的地質(zhì)變量的基礎(chǔ)，而且是統(tǒng)計(jì)對(duì)比和綜合評(píng)價(jià)的對(duì)象；評(píng)價(jià)的結(jié)果是以單元為基礎(chǔ)或通過(guò)單元表達(dá)的。單元?jiǎng)澐值哪康臑榱舜_定地質(zhì)變量觀察尺度和取值范圍，提高評(píng)價(jià)結(jié)果的準(zhǔn)確性，而單元類(lèi)型和大小，猶如樣品采集和分析那樣，其取樣的方法及大小不同，獲得的結(jié)果對(duì)地質(zhì)現(xiàn)象描述的精確程度不同，從而直接影響地質(zhì)研究與資源評(píng)價(jià)的效果。預(yù)測(cè)單元?jiǎng)澐痔?，造成同一地質(zhì)體分布于多個(gè)單元，人為割裂地質(zhì)現(xiàn)象，而且明顯地?cái)U(kuò)大了無(wú)礦單元和單一控礦單元的數(shù)目，增加了預(yù)測(cè)工作量，不利于地質(zhì)模型的建立；而網(wǎng)格單元?jiǎng)澐痔?，則歪曲了有礦單元的分布形態(tài)，使誤判有礦的面積增大，不利于找礦工作的進(jìn)行，并使預(yù)測(cè)靶區(qū)的信度降低，因此如何確定最佳網(wǎng)格單元大小并非易事，它必須結(jié)合實(shí)際資料水平和采用的評(píng)價(jià)模型，選擇合理的單元?jiǎng)澐址椒āＤ壳坝幸?guī)則網(wǎng)格單元?jiǎng)澐趾筒灰?guī)則單元?jiǎng)澐謨煞N方法。規(guī)則網(wǎng)格單元規(guī)則網(wǎng)格化單元是指采用一定規(guī)則的網(wǎng)格作為統(tǒng)計(jì)預(yù)測(cè)樣品的方法總稱(chēng)。自20世紀(jì)60年代D．P．Harsis,首先采用網(wǎng)格法劃分單元建立資源統(tǒng)計(jì)預(yù)測(cè)模型以來(lái)，該方法一直被廣泛地應(yīng)用于資源預(yù)測(cè)、統(tǒng)計(jì)性地質(zhì)數(shù)據(jù)處理等領(lǐng)域。該方法母體構(gòu)成是以“相似類(lèi)比”原理為基礎(chǔ)的，即相似的地質(zhì)構(gòu)造環(huán)境應(yīng)含有相似的礦產(chǎn)。研究區(qū)域內(nèi)所有地質(zhì)構(gòu)造環(huán)境相似的地段都是母體的組成部分，因此，采樣必須沿一定網(wǎng)格連續(xù)進(jìn)行，丟掉任何一個(gè)小的局部地段都有可能漏掉礦化單元。網(wǎng)格化單元有如下優(yōu)點(diǎn)：①網(wǎng)格化單元的劃分原則統(tǒng)一，預(yù)測(cè)區(qū)和模型區(qū)網(wǎng)格單元規(guī)格易統(tǒng)一，它不受地質(zhì)工作程度的不平衡性所限制；②網(wǎng)格化單元是獨(dú)立的，單元內(nèi)已知礦床是隨機(jī)分布的，有利于隨機(jī)抽樣，形成簡(jiǎn)單樣本，同時(shí)還可以研究礦床在單元內(nèi)的分布規(guī)模；③某些礦種（同生礦床系列、石油等）采用網(wǎng)格化單元有利于擴(kuò)大樣品數(shù)，擴(kuò)大樣本容量，以提高統(tǒng)計(jì)的準(zhǔn)確性；④網(wǎng)格化單元?jiǎng)澐址椒ê?jiǎn)便，劃分方案統(tǒng)一，易于坐標(biāo)化，便于大范圍信息的儲(chǔ)存和計(jì)算機(jī)成圖工作。網(wǎng)格單元?jiǎng)澐值年P(guān)鍵問(wèn)題是選擇單元的位置、大小。如何確定最佳網(wǎng)格單元大小，并非易事，也缺乏明確的準(zhǔn)則（趙鵬大，1995）通常對(duì)單元的劃分應(yīng)該考慮以下因素：①預(yù)測(cè)比例尺和精度要求；②預(yù)測(cè)區(qū)地質(zhì)條件復(fù)雜程度，礦點(diǎn)數(shù)及空間分布特征；③研究區(qū)范圍大小及保證統(tǒng)計(jì)分析所需的單元數(shù)；④地質(zhì)特征的空間變異性；③應(yīng)遵循抽樣的隨機(jī)性與樣品的代表性原則。常用的經(jīng)驗(yàn)算法有1經(jīng)驗(yàn)性最佳面積：根據(jù)區(qū)內(nèi)礦點(diǎn)數(shù)量以及預(yù)測(cè)范圍的大小，經(jīng)驗(yàn)性地確定最優(yōu)的單元面積，即S=研究區(qū)總面積／礦點(diǎn)數(shù)*L式中：L為給出的期望礦點(diǎn)平均數(shù)，一般取1一3）之間的值。2趙鵬大等提出一種采用概率統(tǒng)計(jì)劃分方法，即以落入單元內(nèi)的期望礦點(diǎn)數(shù)等于或小于實(shí)際落人單元礦點(diǎn)數(shù)標(biāo)準(zhǔn)差的三倍為準(zhǔn)則，即有

δ／E＝Sqrt（（l－S）／nS）＜＝l／3式中：E為落人單元內(nèi)的礦點(diǎn)數(shù)的數(shù)學(xué)期望；δ為實(shí)際落人單元內(nèi)的礦點(diǎn)數(shù)標(biāo)準(zhǔn)差；S為單元面積與研究區(qū)面積之比；n為礦點(diǎn)總數(shù)。3）根據(jù)圖件比例尺確定單元大?。簩?duì)于1：50萬(wàn)比例尺的地質(zhì)圖，用25－100km2的面積為基本單元的大?。粚?duì)于1：20萬(wàn)比例尺的地質(zhì)圖，用4－16km2的面積為基本單元的大小；對(duì)于I：5萬(wàn)地質(zhì)圖單元大小為0．25-1km2的面積比較適宜。地質(zhì)體單元法也稱(chēng)自然單元法前述單元?jiǎng)澐质且詽M足抽樣的隨機(jī)性或代表性為前提而進(jìn)行的等大小規(guī)則單元?jiǎng)澐?，此法給成礦預(yù)測(cè)帶來(lái)了方便，但可能不能很好反映地層體內(nèi)容的變化特征。不規(guī)則單元——地質(zhì)體、地質(zhì)異常體、資源體等，則可以在不割裂地質(zhì)體內(nèi)部變化性的情況下，開(kāi)展預(yù)測(cè)研究。地質(zhì)單元法是指應(yīng)用對(duì)預(yù)測(cè)礦種具有明顯控制作用的地質(zhì)條件和找礦意義明確的標(biāo)志圈定地質(zhì)統(tǒng)計(jì)單元的方法。其優(yōu)點(diǎn)是能保證直接建立預(yù)測(cè)地質(zhì)體與礦產(chǎn)的關(guān)系。其關(guān)鍵是確定劃分地質(zhì)單元邊界的條件。礦產(chǎn)預(yù)測(cè)的地質(zhì)體單元?jiǎng)澐址椒ㄖ饕Q于綜合信息找礦模型的特點(diǎn)。在綜合信息找礦模型中，有兩種找礦標(biāo)志：一種是成礦的必要條件；另一種是成礦有利（或不利）標(biāo)志。地質(zhì)統(tǒng)計(jì)單元的劃分以成礦的必要條件為基礎(chǔ)，并以成礦有利或不利標(biāo)志為補(bǔ)充。例如：對(duì)于巖漿巖型礦床，多以巖體專(zhuān)屬性特征為主要控礦因素。若礦體的圍巖即是成礦母巖而這類(lèi)礦床的礦體很少產(chǎn)于巖休以外的圍巖中時(shí)，一般選擇巖體為單元。1）將各種自然地質(zhì)體作為單元，這是常用的一種劃分方法。自然地質(zhì)體有地層、巖體和蝕變帶等。例如地層，可將表示地層的那些GIS圖元作為數(shù)據(jù)提取單元。此時(shí)，單元?jiǎng)澐值降牡貙訂挝患?jí)別（例如，是劃分到“群”還是劃分到“組”）應(yīng)由用戶結(jié)合實(shí)際情況確定。另外，如果研究區(qū)域中若干個(gè)地層圖元都帶有用戶指定的某個(gè)屬性項(xiàng)，那么在數(shù)據(jù)提取時(shí)應(yīng)將這若干圖元作為同一個(gè)單元看待，雖然這些圖元分散在GIS圖的不同位置。譬如，某研究區(qū)出露新生代的第四紀(jì)。新近紀(jì)和古近紀(jì)地層及中生代的白堊紀(jì)、侏羅紀(jì)和三疊紀(jì)地層，單元?jiǎng)澐值降募?jí)別是“代”?，F(xiàn)在要考察銅礦（化）點(diǎn)在新生代和中生代地層中的分布情況，那就應(yīng)把圖中那些表示第四紀(jì)、新近紀(jì)和古近紀(jì)地層的圖元作為一個(gè)數(shù)據(jù)提取單元，因?yàn)檫@些地層都屬于新生代；把表示白堊紀(jì)、侏羅紀(jì)和三疊紀(jì)地層的圖元作為另一個(gè)數(shù)據(jù)提取單元，因?yàn)檫@些地層都屬于中生代。如果單元?jiǎng)澐值降募?jí)別是“組”，那么就應(yīng)將同一組的地層作為一個(gè)單元，而不管這個(gè)組在圖上的分布情況。這實(shí)際上是利用了GIS的屬性統(tǒng)計(jì)功能。也可以選擇表示巖體的那些圖元作為數(shù)據(jù)提取單元以查明對(duì)找礦有利的巖體。地層、巖體在GIS中就是以區(qū)圖元實(shí)體的形式存在的，在實(shí)際操作時(shí)就是根據(jù)用戶指定的地層單位級(jí)別依據(jù)屬性將它們聯(lián)系起來(lái)形成一個(gè)組合，同一個(gè)組合中的所有圖元被視為一個(gè)數(shù)據(jù)提取單元。2）根據(jù)研究需要也可把多個(gè)（種）地質(zhì)體或構(gòu)造組合形成的區(qū)域作為單元稱(chēng)之為構(gòu)造單元。有的成礦區(qū)（帶）已經(jīng)根據(jù)本區(qū)的地層、大地構(gòu)造及巖體分布情況完成了構(gòu)造分區(qū)，并生成了相應(yīng)的GIS圖件，我們可以直接把這樣的圖件作為構(gòu)造單元圖。如果要在GIS中完成這樣的工作，顯然是一個(gè)人機(jī)交互過(guò)程，它要求用戶充分研究工作區(qū)的地質(zhì)、地理及構(gòu)造背景，憑借自己的地質(zhì)認(rèn)識(shí)和經(jīng)驗(yàn)劃分預(yù)測(cè)單元，這樣提取出的信息會(huì)更真實(shí)地反映地質(zhì)事實(shí)。構(gòu)造單元邊界通常是地質(zhì)界線，如斷層線、地層界線、巖體界線、不整合面之間的分界線等，在GIS的輔助下它的形成過(guò)程可以概括為：選取地質(zhì)底圖上具有不同地學(xué)意義的原始地質(zhì)界線作為單元的邊界輪廓，適當(dāng)添加輔助線圖元，用一個(gè)新的圖層存放這些線元（選中的原始地質(zhì)界線和輔助線）

經(jīng)過(guò)線編輯、線轉(zhuǎn)弧段、弧段編輯、弧段拓?fù)溴e(cuò)誤檢查、拓?fù)渲亟?，最后生成一個(gè)區(qū)圖層，并將它保存起來(lái)以備后用。例如以褶皺軸為界，將兩翼的地層（實(shí)際上是地層組合）各劃作一個(gè)單元；以主干斷裂、次級(jí)斷裂為界，形成具有不同構(gòu)造背景的構(gòu)造單元。3）以現(xiàn)代地貌為依據(jù)劃分單元，如將匯水盆地（一級(jí)、二級(jí)、三級(jí)匯水域）作為單元。在地球化學(xué)找礦中，經(jīng)常要用匯水盆地來(lái)圈定和研究區(qū)域化探異常。4）以人機(jī)交互的形式將用戶指定的區(qū)域作為單元。如用戶可用鼠標(biāo)圈定某一特定（這一區(qū)域可以是矩形也可以是任意多邊形），將其作為數(shù)據(jù)取值單元，以考察其中的礦點(diǎn)、斷裂。地層和巖體分布情況。

從GIS角度來(lái)講，不管采用哪種單元?jiǎng)澐诸?lèi)型，最終得到的總是一個(gè)區(qū)文件，在MapGIS中是一個(gè)＊．Wp區(qū)文件。在進(jìn)行小比例尺的礦產(chǎn)預(yù)測(cè)時(shí)，因?yàn)轭A(yù)測(cè)工作區(qū)的面積很大，單元?jiǎng)澐謶?yīng)使成礦預(yù)測(cè)的實(shí)際研究范圍大大縮小，使研究具有明確的對(duì)象。因此，單元?jiǎng)澐忠谡业V模型的指導(dǎo)下進(jìn)行，以最大限度地反映成礦信息和預(yù)測(cè)單元面積最小為原則，通過(guò)對(duì)工作區(qū)內(nèi)的地質(zhì)、地球物理、地球化學(xué)、遙感資料的綜合分析，對(duì)可能的成礦地段劃分預(yù)測(cè)單元。二、地質(zhì)信息的定量化提取不同的地學(xué)信息有不同的表達(dá)方式（數(shù)據(jù)。圖形、圖像及文字描述），尤其是地質(zhì)信息主要以地質(zhì)構(gòu)造圖件表達(dá)。如何將圖形、圖像及文字描述數(shù)字化和定量化，如何保證這種轉(zhuǎn)化過(guò)程盡量減少信息的損失和失真，而盡量通過(guò)這種轉(zhuǎn)換增加信息量并減少問(wèn)題的多解性，是實(shí)現(xiàn)定量化成礦預(yù)測(cè)的一個(gè)關(guān)鍵。地質(zhì)信息的定量化提取地質(zhì)變量的取值或地質(zhì)異常的提取以往是一項(xiàng)比較費(fèi)時(shí)也是很復(fù)雜的工作。有些地質(zhì)信息是顯式的，而更多的信息是隱蔽的，需要通過(guò)一定的分析、處理等才能獲得。傳統(tǒng)上，一般是在地質(zhì)圖上打好網(wǎng)格或按不同地質(zhì)體用人工讀圖的方式記錄有關(guān)地層、巖體、構(gòu)造、蝕變及礦化等信息，再填入相關(guān)的表格，不僅費(fèi)工、費(fèi)時(shí)，而且采集的信息精度也不高，如地層出露面積、斷層長(zhǎng)度等是靠讀方格紙估計(jì)的，此外還有一些信息根本就不可能采集，如單元網(wǎng)格內(nèi)的地層熵、斷裂的累計(jì)長(zhǎng)度、優(yōu)異度、成礦多樣度等。GIS系統(tǒng)提供的地質(zhì)變量自動(dòng)取值模塊不僅能有效地解決上述問(wèn)題，而且提高了數(shù)據(jù)采集精度，可靠程度，對(duì)工作效率的提高更是無(wú)法估算的。在實(shí)際工作中，有以下兩種定量取值方法。（1）用空間疊加與查詢(xún)的方式將圖形信息定量化）用空間疊加與查詢(xún)的方法將圖形信息定量化：通過(guò)將研究區(qū)地質(zhì)點(diǎn)、線、面圖層與統(tǒng)計(jì)單元作空間疊加分析，然后運(yùn)用GIS的空間交互查詢(xún)功能，實(shí)現(xiàn)各統(tǒng)計(jì)單元內(nèi)圖形信息定量提取的任務(wù)。如單元內(nèi)有無(wú)某一類(lèi)型的斷層出現(xiàn)。單元內(nèi)有無(wú)某一類(lèi)型的地層出現(xiàn)、單元內(nèi)某一類(lèi)型的地層出露的面積等。（2）構(gòu)造新的地質(zhì)變量通過(guò)熵計(jì)算、地殼升降系數(shù)、線性體統(tǒng)計(jì)、地質(zhì)復(fù)雜度等一系列新變量的定義，提取深層次的地質(zhì)異常（變量）1）熵（Hr）的計(jì)算公式為式中：N為研究區(qū)地質(zhì)對(duì)象種類(lèi)的多少，如地質(zhì)圖中地層、巖漿巖的總種類(lèi)數(shù)；P；為計(jì)算單元中第i類(lèi)地質(zhì)體面積廣(i=l，2，…N）與單元面積之比；Hm為最大可能嫡，又稱(chēng)最大不確定性值。

三證據(jù)權(quán)法證據(jù)權(quán)法（WeightsofEvidence）是一種定量評(píng)價(jià)方法，它最初用于醫(yī)療診斷；20世紀(jì)80年代末，加拿大地調(diào)所FritsAgterberg和GraemeBonham-Carter將該方法引入到GIS支持下的礦產(chǎn)資源潛力評(píng)價(jià)中。其基本過(guò)程是：將每一種地學(xué)信息視為成礦預(yù)測(cè)的一個(gè)證據(jù)因子；通過(guò)分析，計(jì)算出每一個(gè)證據(jù)因子對(duì)成礦貢獻(xiàn)的權(quán)重值；最終對(duì)各證據(jù)因子加權(quán)求和，得出成礦有利度值，從而對(duì)礦產(chǎn)遠(yuǎn)景區(qū)進(jìn)行定位預(yù)測(cè)。

證據(jù)權(quán)法假設(shè)研究區(qū)被劃分成面積相等的Ｔ個(gè)單元，其中有Ｄ個(gè)單元為有礦單元。則隨機(jī)選取一個(gè)單元有礦的概率是：

Ｐ（Ｄ）＝Ｄ/T

則先驗(yàn)概率：

O（D）=Ｐ（Ｄ）/（1-Ｐ（Ｄ））=D/（T-D）對(duì)于任意一個(gè)證據(jù)因子，其權(quán)重定義為：Ｗ＋＝ｌｎ[Ｐ（Ｂ／Ｄ）/Ｐ（Ｂ／Ｄ－）]Ｗ－＝ｌｎ[Ｐ（Ｂ－／Ｄ）/Ｐ（Ｂ－／Ｄ－）]

Ｂ為因子存在區(qū)的單元數(shù)；Ｂ－為因子不存在區(qū)的單元數(shù)，Ｄ表示有礦，Ｄ－表示無(wú)礦Ｗ＋、Ｗ－分別為證據(jù)因子存在區(qū)和不存在區(qū)的權(quán)重值（即成礦關(guān)聯(lián)度），對(duì)于原始數(shù)據(jù)缺失區(qū)域權(quán)重值為０。定義關(guān)聯(lián)度顯著性指標(biāo)C為：

C=W+-W-

證據(jù)權(quán)法證據(jù)權(quán)重法要求各證據(jù)因子之間相對(duì)于礦點(diǎn)分布滿足條件獨(dú)立。對(duì)于Ｎ個(gè)證據(jù)因子，若它們都關(guān)于礦點(diǎn)條件獨(dú)立，則研究區(qū)內(nèi)任一ｋ單元為有礦的后驗(yàn)概率為：由此可得出研究區(qū)內(nèi)任一ｋ單元為有礦單元的后驗(yàn)概率為：

Ｐ＝Ｏ／（１＋Ｏ）最后根據(jù)后驗(yàn)概率圈出找礦遠(yuǎn)景區(qū)。

證據(jù)權(quán)法

工作步驟各類(lèi)圖件證據(jù)層面積計(jì)算點(diǎn)-區(qū)分析點(diǎn)-線分析礦點(diǎn)密度分析權(quán)重計(jì)算W+&W-根據(jù)C值進(jìn)行優(yōu)化

證據(jù)加權(quán)求和先驗(yàn)概率成礦有利度后驗(yàn)概率

證據(jù)權(quán)法工作步驟---1、權(quán)重計(jì)算研究區(qū)某類(lèi)地層落在該地層上的的礦點(diǎn)數(shù)不落在該地層的礦點(diǎn)數(shù)出露該地層的面積未出露該地層的面積

證據(jù)權(quán)法工作步驟---1、權(quán)重計(jì)算

1 50 7 0.7/0.5=1.4 ln(1.4)=+0.332 50 3 0.3/0.5=0.6 ln(0.6)=-0.51Total 100 10=礦床點(diǎn)21某類(lèi)地層不存在某類(lèi)地層存在類(lèi)別 面積 礦點(diǎn)數(shù) 礦點(diǎn)／面積 權(quán)重

證據(jù)權(quán)法工作步驟---1、權(quán)重計(jì)算

1MissingData(25)

3 - 0.02 60 3 0.3/0.60=0.5 ln(0.5)= - 0.693 15 4 0.4/0.15=2.7 ln(2.7)= + 0.98Total 75 10123=礦床點(diǎn)類(lèi)別 面積 礦點(diǎn)數(shù) 礦點(diǎn)／面積 權(quán)重

證據(jù)權(quán)法工作步驟---2、權(quán)重優(yōu)化0.00.51.01.52.02.53.0012345678910顯著度曲線C證據(jù)存在證據(jù)不存在斷層BUFFER距離或物化探異常范圍

證據(jù)權(quán)法工作步驟---3、證據(jù)加權(quán)求和w+w+w+w-w-w-圖層加權(quán)有利度圖最小有利度 < < < < 最大有利度

證據(jù)權(quán)法工作步驟---3、證據(jù)加權(quán)求和

證據(jù)層1 證據(jù)層2 加權(quán)和

A 1 1 0.33 + 0.0= + 0.33B 1 2 0.33 - 0.69= - 0.36C 1 3 0.33 + 0.98= + 1.31E 2 2 -0.51 - 0.69= - 1.20ABCE21+13=證據(jù)層1證據(jù)層2圖層綜合22023/2/4413.“證據(jù)權(quán)”方法的優(yōu)點(diǎn)1.“證據(jù)權(quán)”技術(shù)最終結(jié)果是以權(quán)的形式或以后驗(yàn)概率圖的形式表達(dá)的異常組合圖。2.“證據(jù)權(quán)”技術(shù)的優(yōu)點(diǎn)在于權(quán)的解釋是相對(duì)直觀的，并能夠獨(dú)立的確定，易于產(chǎn)生重現(xiàn)性。該方法亦適用于獲取局部特征和區(qū)域模型的信息（譬如地球化學(xué)和地球物理異常）。2023/2/4424.基于證據(jù)權(quán)法信息集成與資源潛力評(píng)價(jià)（據(jù)陳永清）建立了評(píng)價(jià)模型，借助MORPAS3.0軟件，應(yīng)用證據(jù)權(quán)法對(duì)“三江”南段有色金屬等礦產(chǎn)資源的潛力進(jìn)行綜合信息提取及評(píng)價(jià)；在“三江”南段共圈定Cu-Pb-Zn-Ag找礦遠(yuǎn)景區(qū)20個(gè)，其中A級(jí)遠(yuǎn)景區(qū)6個(gè)、B級(jí)遠(yuǎn)景區(qū)7個(gè)、C級(jí)遠(yuǎn)景區(qū)7個(gè)。2023/2/443（MORPAS3.0）礦產(chǎn)資源評(píng)價(jià)分析系統(tǒng)2023/2/4442023/2/445地層含礦性分析2023/2/446地層W+W-CE0.448-0.0810.529K-0.1750.057-0.232J-0.1300.075-0.205T0.230-0.2460.476P0.191-0.0940.285C-P-0.3650.017-0.382C0.201-0.0230.224D-C0.265-0.0250.290D0.356-0.0480.404S-D0.876-0.0620.938S-0.6930.022-0.715O-S0.311-0.0150.326O0.694-0.0360.730∈0.876-0.0350.911Z-∈-0.4320.006-0.438Pt0.102-0.0490.1512023/2/447巖性數(shù)W+W-C1~2-1.1940.106-1.3003~4-0.1690.047-0.2165~6-0.0810.031-0.1127~80.395-0.1300.5259~100.539-0.0720.61111~120.224-0.0050.22913~140.945-0.0100.955巖性多樣性2023/2/448線性構(gòu)造控礦性分析線性構(gòu)造W+W-C線性構(gòu)造0.5KmBuffer0.322-0.7121.034線性構(gòu)造1KmBuffer0.158-1.5751.733線性構(gòu)造1.5KmBuffer0.123-2.4782.601線性構(gòu)造2KmBuffer0.093-2.2422.335線性構(gòu)造2.5KmBuffer0.121-1.8421.963線性構(gòu)造——包括實(shí)測(cè)的斷裂構(gòu)造和用遙感和物探方法推測(cè)的隱伏斷裂構(gòu)造2023/2/449環(huán)形構(gòu)造W+W-C環(huán)形構(gòu)造1KmBuffer0.301-0.5860.887環(huán)形構(gòu)造2KmBuffer0.260-0.6230.883環(huán)形構(gòu)造3KmBuffer0.236-0.7150.951環(huán)形構(gòu)造4KmBuffer0.174-0.6370.811環(huán)形構(gòu)造5KmBuffer0.133-0.5440.677遙感解譯環(huán)形構(gòu)造控礦性分析2023/2/450中酸性巖體W+W-C中酸性巖體1KmBuffer0.343-0.2480.591中酸性巖體2KmBuffer0.364-0.3110.675中酸性巖體3KmBuffer0.357-0.3420.699中酸性巖體4KmBuffer0.317-0.3260.643中酸性巖體5KmBuffer0.276-0.3050.581中酸性巖體含礦性分析2023/2/451基性超基性巖體W+W-C基性超基性巖體1KmBuffer0.508-0.0410.549基性超基性巖體2KmBuffer0.617-0.0610.678基性超基性巖體3KmBuffer0.539-0.0640.603基性超基性巖體4KmBuffer0.558-0.0800.638基性超基性巖體5KmBuffer0.531-0.0830.614基性超基性巖體含礦性分析2023/2/452

化探因子W+W-CF1Co-Cr-Cu-As-V0.278-0.2230.501F2Pb-Zn-Ag0.510-0.2830.793F5Au0.038-0.0430.081地球化學(xué)控礦性分析2023/2/453（后驗(yàn)概率）單元總數(shù)：2484含礦單元數(shù)：191先驗(yàn)概率：0.082023/2/454序號(hào)證據(jù)圖層W+W-C1線性構(gòu)造Buffer1.5Km0.123-2.4782.6012環(huán)形構(gòu)造Buffer3Km0.236-0.7150.9513含礦地層組合熵0.205-0.5180.7234巖性多樣性（巖性種類(lèi)數(shù)7~14）0.441-0.2970.7385中酸性巖體Buffer3Km0.357-0.3420.6996基性超基性巖體Buffer4Km0.558-0.0800.6387Cu異常(SA)0.537-0.4791.0168Pb-Zn-Ag異常0.510-0.2830.7939Au異常(SA)0.316-0.2080.525西南“三江”南段資源潛力評(píng)價(jià)證據(jù)權(quán)模型2023/2/455Cu-Pb-Zn-Ag找礦靶區(qū)(Ⅰ類(lèi))A級(jí)6個(gè)

B級(jí)7個(gè)

C級(jí)7個(gè)Au-PGE找礦靶區(qū)（Ⅱ類(lèi)）

A級(jí)4個(gè)有色金屬貴金屬綜合找礦靶區(qū)（Ⅲ類(lèi)）

A級(jí)4個(gè)

B級(jí)1個(gè)

C級(jí)4個(gè)2023/2/456西南“三江”南段Cu、Pb、Zn、Ag、Au等找礦靶區(qū)分類(lèi)表類(lèi)型級(jí)別編號(hào)面積（km2）典型礦床Ⅰ類(lèi)找礦靶區(qū)AI-A-1812.36德欽羊拉銅礦I-A-21229.92騰沖老廠坪子銅鉛鋅多金屬礦I-A-31041.53保山核桃坪鉛鋅礦I-A-4523.40思茅大坪掌銅礦BI-B-1203.21I-B-2605.45I-B-3651.00I-B-4542.27I-B-5273.60I-B-6169.69I-B-7139.24I-B-8396.91I-B-9141.282023/2/457類(lèi)型級(jí)別編號(hào)面積（km2）典型礦床I類(lèi)CI-C-1131.31I-C-2604.84Ⅱ類(lèi)靶區(qū)AⅡ-A-196.24金寶山鉑鈀礦Ⅱ-A-2475.90老王寨金礦BⅡ-B-1213.50Ⅲ類(lèi)靶區(qū)AⅢ-A-11843.96中甸普朗銅礦、雪雞坪銅礦Ⅲ-A-23898.59金頂鉛鋅礦、白秧坪銅礦、北衙金礦Ⅲ-A-31712.35鎮(zhèn)康蘆子園鉛鋅礦、勐興鉛鋅礦BⅢ-B-1187.35Ⅲ-B-2474.64Ⅲ-B-31403.872023/2/458實(shí)例--川西北地區(qū)證據(jù)權(quán)法的應(yīng)用

——第一步，數(shù)據(jù)庫(kù)建立預(yù)測(cè)的礦種為金礦，類(lèi)型為微細(xì)浸染型