王春女，陳建平，韓振哲，李永勝，賈文彬

(1.中國(guó)地質(zhì)調(diào)查局發(fā)展研究中心，北京 100037；2.國(guó)土資源部礦產(chǎn)勘查技術(shù)指導(dǎo)中心，北京 100120；3.中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083)

礦產(chǎn)預(yù)測(cè)中的不確定性評(píng)估方法與應(yīng)用:以贛北冷水坑礦區(qū)為例

王春女1,2，陳建平3，韓振哲1,2，李永勝1,2，賈文彬3

(1.中國(guó)地質(zhì)調(diào)查局發(fā)展研究中心，北京 100037；2.國(guó)土資源部礦產(chǎn)勘查技術(shù)指導(dǎo)中心，北京 100120；3.中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院，北京 100083)

礦產(chǎn)預(yù)測(cè)過(guò)程一方面應(yīng)用地質(zhì)成礦基礎(chǔ)理論，另一方面充分結(jié)合地質(zhì)、地球化學(xué)、地球物理、遙感等信息數(shù)據(jù)所反映的成礦規(guī)律及蘊(yùn)含信息，而信息處理及預(yù)測(cè)過(guò)程存在一定程度的不確定性，因此不確定性的正確評(píng)估對(duì)找礦預(yù)測(cè)工作和礦業(yè)經(jīng)濟(jì)具有重要意義。本文分析了礦產(chǎn)預(yù)測(cè)不確定性的來(lái)源和一般評(píng)價(jià)方法，從預(yù)測(cè)過(guò)程的內(nèi)部影響因素和地質(zhì)影響因素兩個(gè)方面出發(fā)，提出了基于模糊集理論的礦產(chǎn)預(yù)測(cè)不確定性評(píng)價(jià)方法，并將其在贛北冷水坑礦區(qū)三維礦產(chǎn)預(yù)測(cè)中進(jìn)行了應(yīng)用。該方法對(duì)預(yù)測(cè)的五個(gè)遠(yuǎn)景區(qū)不確定性評(píng)估效果良好，與實(shí)際地質(zhì)情況一致，從而驗(yàn)證了評(píng)價(jià)方法的有效性和三維預(yù)測(cè)方法的可行性。

礦產(chǎn)預(yù)測(cè);不確定性;冷水坑;方法

隨著找礦難度的加大，礦產(chǎn)預(yù)測(cè)的不確定性研究日益得到廣泛重視。目前，礦產(chǎn)預(yù)測(cè)不確定性的研究主要集中在兩方面，一方面是不確定性來(lái)源的研究(趙鵬大，1996；Donald A.Singer，2001；Bardossy and Fodor，2004)[1-3]；另一方面是評(píng)價(jià)不確定性的方法(桂寶林，1984；Cheng and Agterberg，1999；Bardossy and Fodor，2004；Zuo,R,etal.，2009；張善明，2012)[4-7]。此外，一些學(xué)者已將不確定性評(píng)價(jià)進(jìn)行了一些實(shí)例應(yīng)用(Laurent Menut，2008；方維萱，2009；Frits Agterberg，2011)[8-10]。因此，結(jié)合地質(zhì)信息固有屬性和礦產(chǎn)預(yù)測(cè)過(guò)程，研究礦產(chǎn)預(yù)測(cè)中的不確定性分析方法，對(duì)開展礦產(chǎn)勘查，降低區(qū)域找礦風(fēng)險(xiǎn)具有重要意義。

1 礦產(chǎn)預(yù)測(cè)不確定性來(lái)源

礦產(chǎn)預(yù)測(cè)是在充分收集預(yù)測(cè)區(qū)的地質(zhì)、地球物理、地球化學(xué)、遙感等信息的基礎(chǔ)上，以成礦理論為指導(dǎo)，梳理研究成礦規(guī)律，分析提取成礦要素，應(yīng)用數(shù)學(xué)統(tǒng)計(jì)方法和計(jì)算機(jī)技術(shù)相結(jié)合，從而確定成礦靶區(qū)并對(duì)礦床數(shù)量和資源量進(jìn)行評(píng)估。礦產(chǎn)預(yù)測(cè)是一個(gè)復(fù)雜多階段的工作過(guò)程，它在進(jìn)行的每一階段都可能產(chǎn)生不確定性，按性質(zhì)不確定性來(lái)源可分為地質(zhì)事物的固有屬性和人為工作的誤差屬性，按階段可分為數(shù)據(jù)收集不充分性、理論模型模糊性、操作方法不精確性、專家意見不準(zhǔn)確性[11]。綜合來(lái)講，不確定性的來(lái)源主要有以下五個(gè)方面。

1)地質(zhì)自然屬性。由于地質(zhì)作用本身具有多期性、復(fù)雜性，其對(duì)礦床形成的影響也是多因素和多階段的，既體現(xiàn)在成礦期的物質(zhì)富集、成礦流體、成礦動(dòng)力系統(tǒng)等過(guò)程因素，又體現(xiàn)在礦床形成后礦床的保存環(huán)境等因素，從而導(dǎo)致了成礦過(guò)程和地質(zhì)特征的不確定性。

2)工作手段。一方面，地質(zhì)工作中所獲得的地質(zhì)信息一般是來(lái)自一定范圍內(nèi)具有代表性的樣品數(shù)據(jù)，而樣品數(shù)據(jù)代表性的選擇存在一定的不確定性。另一方面，地質(zhì)工作中所選擇的工作方法和工作手段的準(zhǔn)確性和適應(yīng)性也存在一定程度的不確定性。

3)數(shù)據(jù)質(zhì)量。首先，數(shù)據(jù)在獲取過(guò)程中受儀器精度、人為操作誤差和環(huán)境干擾等因素的限制和影響，所以收集到得數(shù)據(jù)本身就具有不確定性。其次資料的充分性，比如詳細(xì)程度、比例尺都將對(duì)數(shù)據(jù)質(zhì)量產(chǎn)生一定影響。

4)信息處理。信息的轉(zhuǎn)化、加工一般是建立在一定數(shù)學(xué)模型和數(shù)學(xué)算法的基礎(chǔ)進(jìn)行的，而建立數(shù)學(xué)模型的完善性，算法選擇的適用性和精確性都存在一定程度的不確定性。此外，預(yù)測(cè)過(guò)程中的信息提取、信息綜合分析等都是依托GIS軟件的空間分析功能完成，也會(huì)產(chǎn)生不確定性，而且這些不確定性還可能被放大、分散、傳遞或積累，從而導(dǎo)致相當(dāng)數(shù)量的不確定性[12-13]。

5)主觀因素。在礦產(chǎn)預(yù)測(cè)定位過(guò)程中找礦模型的建立、有利成礦信息的選擇都是依據(jù)人類對(duì)地質(zhì)現(xiàn)象的認(rèn)識(shí)程度和知識(shí)的完備程度，這些過(guò)程都會(huì)由于人類認(rèn)識(shí)的不確定性產(chǎn)生不同程度的信息失真或不確定性。在礦產(chǎn)預(yù)測(cè)資源量的估計(jì)階段，人類對(duì)資源量估算方法的選擇導(dǎo)致了預(yù)測(cè)資源量的不確定性。

2 礦產(chǎn)預(yù)測(cè)不確定性評(píng)價(jià)方法

傳統(tǒng)的評(píng)價(jià)礦產(chǎn)預(yù)測(cè)不確定性的方法有概率理論和地質(zhì)統(tǒng)計(jì)學(xué)法。概率理論研究隨機(jī)現(xiàn)象數(shù)量規(guī)律的方法。在礦產(chǎn)預(yù)測(cè)過(guò)程中，通過(guò)對(duì)比發(fā)現(xiàn)有些變量在每個(gè)模型單元都會(huì)出現(xiàn)，這實(shí)際上是一個(gè)確定性的事件，稱之為確定性地質(zhì)事件，確定性地質(zhì)事件對(duì)礦產(chǎn)資源量預(yù)測(cè)不起作用，但對(duì)礦產(chǎn)預(yù)測(cè)的靶區(qū)確定和單元邊界條件的確定有重要作用；而有些地質(zhì)變量只在某些模型單元中出現(xiàn)，出現(xiàn)概率介于0～1之間，稱之為統(tǒng)計(jì)地質(zhì)變量，它是礦產(chǎn)資源量預(yù)測(cè)的關(guān)鍵性地質(zhì)變量。概率論方法評(píng)價(jià)不確定性要做的就是確定這些統(tǒng)計(jì)地質(zhì)變量。地質(zhì)統(tǒng)計(jì)學(xué)法是以區(qū)域化變量為基礎(chǔ)，借助變異函數(shù)，研究既具有隨機(jī)性又具有結(jié)構(gòu)性，或空間相關(guān)性和依賴性的自然現(xiàn)象的一門科學(xué)。目前，模糊集是礦產(chǎn)預(yù)測(cè)不確定性評(píng)價(jià)較優(yōu)越的一種評(píng)價(jià)方法。

模糊集概念是美國(guó)控制論專家L.A.Zaden于1965年提出，他把普通集合論中μ要么屬于A、要么不屬于A的兩種絕對(duì)概念靈活化，用隸屬度來(lái)代替絕對(duì)的“屬于”或“不屬于”關(guān)系，表示屬于程度是多少，取值可取0～1之間的任意實(shí)數(shù)值。也就是說(shuō)，把隸屬度關(guān)系從普遍集合論中只取“0”或“1”兩個(gè)值推廣到[0，1]閉區(qū)間。

設(shè)A是集合X到[0，1]的一個(gè)映射，A：X→[0，1]，X→A(x) 則稱A是X上的模糊集，A(x)稱為模糊集A的隸屬函數(shù)，或稱A(x)為x對(duì)模糊集A的隸屬度。在礦產(chǎn)預(yù)測(cè)中，隸屬度的大小表示地質(zhì)變量、預(yù)測(cè)結(jié)果的可靠性及誤差的大小。

考慮一般情況，礦產(chǎn)預(yù)測(cè)評(píng)價(jià)中使用了不確定性和影響因素的信息，根據(jù)模糊集理論，定義如式(1)所示模型結(jié)構(gòu)[14]。

(1)

式中：X=(X1,X2,…,Xm)，代表不確定性的隨機(jī)輸入變量；Y=(Y1,Y2,…,Yn)，代表確定性輸入；μ代表礦產(chǎn)預(yù)測(cè)評(píng)價(jià)的輸出結(jié)果。

3 贛北冷水坑礦區(qū)三維成礦預(yù)測(cè)不確定性分析

目前，我國(guó)主要的礦產(chǎn)資源預(yù)測(cè)理論與方法有地質(zhì)異常(趙鵬大，等，1991)[15]、證據(jù)權(quán)重法(Agterberg，etal.，1989)[16，17]、地質(zhì)統(tǒng)計(jì)學(xué)(王仁鐸，1985；孫英君，等，2004)[18，19]、綜合信息礦產(chǎn)預(yù)測(cè)法(王世稱，2010)[20]等。傳統(tǒng)的以平面圖和剖面圖為主的地質(zhì)信息模擬與表達(dá)已難以滿足現(xiàn)代礦山信息化發(fā)展的迫切需要。本文提出的三維成礦預(yù)測(cè)方法是一種基于真三維地質(zhì)建模的隱伏礦體預(yù)測(cè)方法，是在傳統(tǒng)二維預(yù)測(cè)的基礎(chǔ)上，將二維成礦預(yù)測(cè)中被認(rèn)可的變量分析拓展到三維空間內(nèi)，從而為三維成礦預(yù)測(cè)提供新的變量，在實(shí)現(xiàn)可視化的同時(shí)實(shí)現(xiàn)定位定量的精確預(yù)測(cè)。

3.1 研究區(qū)三維成礦預(yù)測(cè)

研究區(qū)的三維成礦預(yù)測(cè)是借助于“立方體模型”來(lái)實(shí)現(xiàn)傳統(tǒng)的二維找礦向三維找礦的新突破。以三維數(shù)字模型為基礎(chǔ)，地質(zhì)統(tǒng)計(jì)方法為手段，基于GIS技術(shù)對(duì)成礦信息進(jìn)行多元提取與分析，確定有利成礦條件。采用信息量的方法計(jì)算成礦有利度，根據(jù)成礦有利度進(jìn)行分級(jí)，圈定成礦遠(yuǎn)景區(qū)。研究具體流程如下所示。

1)對(duì)已有地質(zhì)研究資料和礦區(qū)生產(chǎn)資料進(jìn)行收集與整理，建立矢量化數(shù)據(jù)庫(kù)。

2)根據(jù)研究區(qū)域地質(zhì)背景及前人研究成果的梳理，總結(jié)出該地區(qū)的成礦規(guī)律，建立找礦模型。

3)數(shù)字區(qū)域三維建模。選擇基于平行剖面的建模方法，即通過(guò)相鄰剖面之間對(duì)應(yīng)的輪廓線連接來(lái)模擬地質(zhì)體的形態(tài)，分別建立了贛北冷水坑礦田及外圍地形、地層、構(gòu)造、巖體以及已知礦體的三維實(shí)體模型。

4)成礦有利信息提取。根據(jù)搜集到的資料，確定找礦標(biāo)志，分別將找礦標(biāo)志變量與已知礦體進(jìn)行對(duì)比，提取出有利成礦的三維預(yù)測(cè)變量并將進(jìn)行立方體預(yù)測(cè)單元的提取。

5)三維定位預(yù)測(cè)。根據(jù)對(duì)研究區(qū)礦床的成礦條件分析及有利成礦信息提取，建立研究區(qū)礦產(chǎn)預(yù)測(cè)評(píng)價(jià)模型，確定研究區(qū)的主要控礦因素。將找礦數(shù)字模型所確定的控礦因素作為屬性賦給每一個(gè)單元塊，計(jì)算每個(gè)地質(zhì)單元所包含的綜合找礦信息。根據(jù)提取的成礦有利信息，結(jié)合最后我們根據(jù)每種找礦標(biāo)志的相對(duì)權(quán)值計(jì)算出該立方塊單元的綜合信息量，在研究區(qū)內(nèi)的三維空間范圍內(nèi)圈定成礦遠(yuǎn)景區(qū)(圖1)[21-23]。

6)預(yù)測(cè)分析對(duì)預(yù)測(cè)結(jié)果及控礦要素進(jìn)行進(jìn)一步分析總結(jié)，實(shí)現(xiàn)研究區(qū)定性、定量、定位及定概率預(yù)測(cè)。

3.2 不確定信息分析

綜合本次研究的實(shí)際情況和預(yù)測(cè)過(guò)程，我們從預(yù)測(cè)過(guò)程的內(nèi)部因素和地質(zhì)因素對(duì)成礦的作用兩方面進(jìn)行分析，最后兩部分因素進(jìn)行疊加綜合處理得出本次礦產(chǎn)預(yù)測(cè)的不確定性。

3.2.1 內(nèi)部因素不確定信息分析

針對(duì)本次研究的實(shí)際情況對(duì)式(1)中X進(jìn)行定量分析，確定X和μ之間的關(guān)系。假設(shè)理想狀態(tài)下，影響因素對(duì)礦產(chǎn)預(yù)測(cè)結(jié)果無(wú)影響，將其值賦為1，依據(jù)其內(nèi)部精度的預(yù)測(cè)結(jié)果的影響依次對(duì)影響因素進(jìn)行賦值。綜合考慮本次研究的各個(gè)方面，我們確定了下面幾個(gè)影響不確定性的因素。

1)工程密度。工程密度涉及到勘探線間距和鉆孔間距?？碧骄€間距是指在地質(zhì)體建模時(shí)所用到的各個(gè)剖面之間的距離，間距越小，所建的實(shí)體模型越逼真，其精度也越高。將工程間距劃分不同的等級(jí)，如表1所示。本次建模所使用的勘探線剖面有兩部分組成，一部分是已知礦區(qū)的實(shí)測(cè)勘探線剖面，間距為200m，一部分是用來(lái)區(qū)域建模的圖切剖面，間距是2000m。

表1 工程間距賦值表

2)數(shù)據(jù)質(zhì)量。主要涉及資料精度和數(shù)據(jù)完整性兩個(gè)方面。剖面精度涉及到剖面的比例尺問(wèn)題，從某礦區(qū)到區(qū)域上比例尺可以是1∶500、1∶1000、1∶2000、1∶5000、1∶1萬(wàn)、1∶5萬(wàn)、1∶10萬(wàn)和1∶20萬(wàn)，顯然，比例尺越大，圖面表現(xiàn)的地質(zhì)信息越詳細(xì)，從而預(yù)測(cè)結(jié)果的精度越高。將以上比例尺劃分不同的等級(jí)，根據(jù)相對(duì)精度進(jìn)行賦值，見表2。本次建模最重要的是勘探線剖面，比例尺越大，模型精度越高，比例尺主要有1∶5萬(wàn)和1∶2000兩種，其中區(qū)域地質(zhì)模型建模主要應(yīng)用的是比例尺為1∶5萬(wàn)的勘探線剖面，已知礦體建模應(yīng)用的是比例尺為1∶2000的勘探線剖面。本次研究收集到的數(shù)據(jù)以1∶5萬(wàn)的區(qū)域地質(zhì)調(diào)查和礦產(chǎn)地質(zhì)調(diào)查數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)，收集的數(shù)據(jù)內(nèi)容包括實(shí)際材料圖、地質(zhì)礦產(chǎn)圖、物探數(shù)據(jù)、化探數(shù)據(jù)和文字資料，缺乏遙感資料，比較全面，數(shù)據(jù)完整性賦值為0.9。

圖1 技術(shù)路線圖

表2 比例尺賦值表

3)礦塊大小。根據(jù)建模經(jīng)驗(yàn)及軟件公司的建議，對(duì)于礦塊大小的確定，可根據(jù)勘探線的網(wǎng)度、礦體的大小，礦體邊界的復(fù)雜度及采礦設(shè)計(jì)的要求。一般情況下，礦塊大小可取勘探線間距的1/5～1/10，或者礦塊大小可以設(shè)置為采礦時(shí)的一個(gè)礦塊大小(如一個(gè)臺(tái)階)等。因此，礦塊劃分的太大或者太小對(duì)于模型精度和最終的預(yù)測(cè)結(jié)果均有影響，選擇適當(dāng)大小的塊體極為重要，使它即能保證精度又能保證塊體模型順利為預(yù)測(cè)服務(wù)。本研究采用的塊體大小為100m×100m×50m，雖然按一般礦塊劃分原則稍大，但由于本研究中區(qū)域面積較大，導(dǎo)致范圍塊體量較多，數(shù)據(jù)量會(huì)較大，因此這個(gè)規(guī)格符合實(shí)際研究需要。將這一指標(biāo)精度賦值為0.7。

4)找礦模型。找礦模型是工作人員對(duì)于該區(qū)域地質(zhì)找礦研究的認(rèn)識(shí)程度，能否準(zhǔn)確全面地概括出研究區(qū)域的找礦模型，在一定程度上影響著預(yù)測(cè)結(jié)果的精度。本研究找礦模型的提出是建立在專家對(duì)研究區(qū)域地質(zhì)背景和成礦條件充分認(rèn)知總結(jié)基礎(chǔ)上的，精度相對(duì)較高，故賦值為0.7。

為了定量化計(jì)算μ的精度指標(biāo)，可用式(2)先確定f(Xm,Yn)的值。

(2)

式中，Yn為影響因素的重要程度，Y值的確定依據(jù)影響不確定性各個(gè)因素對(duì)預(yù)測(cè)結(jié)果影響的大小進(jìn)行賦值，對(duì)影響因素按照重要程度進(jìn)行分級(jí)，我們采用專家打分法確定影響因素的權(quán)重(表3)。從表4看出，影響最為重要的是工程密度，賦權(quán)重值為3.5；其次是找礦模型和剖面精度，賦權(quán)重值分別為3和2；次重要的是數(shù)據(jù)完整性和礦塊大小，賦權(quán)重值為1.25和1.5。

表3 專家打分權(quán)重表

表4 不確定性因素賦值及權(quán)重值表

根據(jù)不確定性因素賦值及權(quán)重值表，見式(3)。

(3)

3.2.2 地質(zhì)因素的不確定性分析

由于成礦地質(zhì)過(guò)程的復(fù)雜性，每個(gè)地質(zhì)因素都對(duì)礦床形成具有不同的貢獻(xiàn)值，并具有不確定性，因此我們用成礦有利度來(lái)表示每個(gè)地質(zhì)因素對(duì)成礦作用的貢獻(xiàn)值。根據(jù)對(duì)大量資料及前人研究成果的梳理總結(jié)，確定贛北冷水坑鉛鋅銀礦的主要地質(zhì)控礦因素包括地層(早白堊世下渡組、新元古代南華紀(jì)下峰巖組)、北東向斷裂及其與火山斷裂的重迭處；巖漿巖(燕山早期第二階段正長(zhǎng)花崗巖、燕山晚期次火山旋回花崗斑巖、加里東期混合花崗巖)[24]。控礦因素成礦有利度賦值見表5。

表5 控礦因素成礦有利度賦值表

根據(jù)遠(yuǎn)景區(qū)內(nèi)控礦地質(zhì)因素的情況，地質(zhì)因素的可靠性可用式(4)表達(dá)。

(4)

則冷水坑三維預(yù)測(cè)的可靠性可表達(dá)為預(yù)測(cè)過(guò)程的內(nèi)部因素和地質(zhì)因素的綜合評(píng)價(jià)，可用式(5)進(jìn)行表達(dá)。

(5)

研究所使用的資料比例尺統(tǒng)一、礦塊大小一致、單一找礦模型指導(dǎo)找礦，因此從表達(dá)式可以看出，幾個(gè)成礦遠(yuǎn)景區(qū)的內(nèi)部因素的不確定性基本相同，主要影響不確定性變化的是遠(yuǎn)景區(qū)內(nèi)出露的地質(zhì)體。以遠(yuǎn)景區(qū)A-1為例，該遠(yuǎn)景區(qū)內(nèi)存在已知礦體，因此工程間距涉及兩部分，一部分是已知礦區(qū)的勘探線剖面，間距為200m，一部分是用來(lái)區(qū)域建模的圖切剖面，間距是2000m，因此，工程密度賦值取0.8和0.5的中間值0.65。同理，遠(yuǎn)景區(qū)內(nèi)的資料比例尺也涉及兩部分，區(qū)域建模的比例尺是1∶50000，礦體建模的比例尺是1∶2000，所以資料的綜合精度賦值為0.75。遠(yuǎn)景區(qū)內(nèi)出露的地質(zhì)體有早白堊系下渡組、燕山晚期花崗斑巖、燕山早期第二階段正長(zhǎng)花崗巖、加里東期混合花崗巖、并有區(qū)域性大斷裂湖石斷裂、冷水坑斷裂貫徹區(qū)內(nèi)。

則遠(yuǎn)景區(qū)A-1的可靠性如式(6)所示。

(6)

因此，在遠(yuǎn)景區(qū)A-1找礦的最大風(fēng)險(xiǎn)為24.27%。按照此方法，我們可以依次計(jì)算出其他六個(gè)遠(yuǎn)景區(qū)的三維成礦預(yù)測(cè)的不確定性，分別為：46.77%、40.38%、46.78%、40.66%、55.86%、46.77%。遠(yuǎn)景區(qū)A-1的成礦的可靠性最高，該區(qū)的東北部已發(fā)現(xiàn)冷水坑大型鉛鋅銀礦礦床，區(qū)內(nèi)的控礦地層與巖體占遠(yuǎn)景區(qū)面積的一半以上，區(qū)域性大斷裂貫穿整個(gè)遠(yuǎn)景區(qū)，成礦地質(zhì)條件良好；兩個(gè)C級(jí)遠(yuǎn)景區(qū)成礦的風(fēng)險(xiǎn)性最高，主要由自身的地質(zhì)成礦條件決定的，區(qū)內(nèi)只有少數(shù)小斷裂和小部分控礦巖體。研究區(qū)成礦遠(yuǎn)景區(qū)預(yù)測(cè)結(jié)果如圖2所示。由此可見，該方法能夠?qū)Τ傻V預(yù)測(cè)的不確定性和風(fēng)險(xiǎn)性進(jìn)行較好的定量評(píng)價(jià)。但是該方法也存在一些缺點(diǎn)，如：內(nèi)在因素的權(quán)重值判定雖然是根據(jù)專家打分的平均值，但是仍然存在一定的主觀因素，屬于以主觀判定主觀，評(píng)估方法還有進(jìn)一步研究改善的空間。

圖2 研究區(qū)成礦遠(yuǎn)景區(qū)預(yù)測(cè)結(jié)果

4 結(jié) 論

1)礦產(chǎn)預(yù)測(cè)的不確定性隨著人們對(duì)成礦系統(tǒng)的知識(shí)和信息的掌握程度變化而改變。因此，減少礦產(chǎn)預(yù)測(cè)的不確定性可以從以下幾方面進(jìn)行著手：加強(qiáng)科技研究，提高儀器的精度和可靠程度；深入研究研究區(qū)的成礦地質(zhì)背景和礦床成因，增強(qiáng)找礦模型的準(zhǔn)確性；提高數(shù)據(jù)質(zhì)量，建立全面的空間數(shù)據(jù)庫(kù)，應(yīng)用GIS技術(shù)提高信息處理的精度水平，減少信息處理過(guò)程中的不確定性。

2)本文從數(shù)據(jù)質(zhì)量、主觀因素兩個(gè)不確定性來(lái)源的內(nèi)部因素精度對(duì)礦產(chǎn)預(yù)測(cè)的不確定性進(jìn)行了定量估計(jì)，方法簡(jiǎn)單易操作，且符合模糊集理論。該方法雖然是以贛北冷水坑礦區(qū)為例，但是其計(jì)算理論和方法對(duì)其他礦區(qū)也試用。

3)礦產(chǎn)預(yù)測(cè)的不確定性評(píng)價(jià)除了從來(lái)源方向進(jìn)行評(píng)估，還應(yīng)考慮經(jīng)濟(jì)價(jià)值、環(huán)境、社會(huì)等宏觀因素的影響。

4)本次礦產(chǎn)預(yù)測(cè)研究工作，除了兩個(gè)C級(jí)遠(yuǎn)景區(qū)，其他遠(yuǎn)景區(qū)的成礦可靠性都在60%左右，特別是遠(yuǎn)景區(qū)A-1的可靠性高達(dá)70%以上，并與區(qū)內(nèi)已知礦體進(jìn)行了很好的互相驗(yàn)證，說(shuō)明基于立方塊的三維預(yù)測(cè)方法的可靠性較高，既很好的呈現(xiàn)了地質(zhì)體的形態(tài)，實(shí)現(xiàn)了可視化，又準(zhǔn)確的對(duì)研究區(qū)進(jìn)行了定位定量預(yù)測(cè)，并且預(yù)測(cè)成果可靠性較高，是一種值得深入研究與推廣的礦產(chǎn)預(yù)測(cè)方法。

[1] 趙鵬大，池順都，陳永清.查明地質(zhì)異常:成礦預(yù)測(cè)的基礎(chǔ)[J].高校地質(zhì)學(xué)報(bào),1996，2(4):361-372.

[2] Donald A.Singer.資源定量評(píng)價(jià)發(fā)展方向展望[J].地球科學(xué),2001,26(2):152-156.

[3] G.Bardossy，J.Fodor.Evaluation of uncertainty and risksin geology[M].Berlin:Springer Group,2004:3-37.

[4] 桂寶林.礦產(chǎn)資源預(yù)測(cè)評(píng)價(jià)精確性與可靠性的估計(jì)原則[J].地質(zhì)論評(píng)，1984,30(6):536-542.

[5] Cheng,Q.，Agterberg,F.P.Fuzzy weights of evidence method and its application in mineral potential mapping[J].Natural resources researeh,1999,8(1):27-35.

[6] Zuo,R.Cheng,Q.Agterber,F.P.Application of a hybrid method combining multilevel fuzzy comprehensive evaluation with asymmetric fuzzy relation analysis to mapping prospectivity[J].Ore Geology Reviews,2009，35:101-108.

[7] 張善明,呂新彪,唐小春，等.礦產(chǎn)預(yù)測(cè)中的不確定性評(píng)估方法與應(yīng)用-以廣西大瑤山礦區(qū)為例[J].地球信息科學(xué)學(xué)報(bào)，2012,12(3):342-347.

[8] Laurent Menut，Isabelle Chiapello，Cyril Moulin.Predictability of mineral dust concentrations:The African Monsoon Multidisciplinary Analysis first short observation period forecasted with CHIMERE-DUST[J].JOURNAL OF GEOPHYSICAL RESEARCH，2009,114.

[9] 方維萱，郭玉乾.基于風(fēng)險(xiǎn)分析的商業(yè)性找礦預(yù)測(cè)新方法與應(yīng)用[J].地學(xué)前緣，2009，16(2)：209-226.

[10] Frits Agterberg.Principles of Probabilistic Regional Mineral Resource Estimation[J].Earth Science-Journal of China University of Geosciences，2011,36(2):189-200.

[11] 左仁廣.基于地質(zhì)異常的礦產(chǎn)資源定量化預(yù)測(cè)與不確定性評(píng)價(jià)[D].武漢：中國(guó)地質(zhì)大學(xué)(武漢)，2009.

[10] Dominy,S.C.and Johansen,G.F.,.Reducing grade unertainty in high-nugget effect gold veins:application of geological and geochemical proxies[J].Publicaton Series-Australasian Institute of Mining and Metallurgy 2004:291-302.[12] D.A.Singer，W.D.Menzie,.Statistical guides to estimating the number of undiscovered mineral deposit：an example with porphyry copper deposits[J].Proceedings of IAMG'05:GIS and Spatial Analysis，2005.

[13] Alok Porwal，E.J.M.Carranza,M.Hal.Extended Weights-of-Evidence Modelling for Predictive Mapping of Base Metal Deposit Potential in Aravalli Province,Western India，Explor[J].Mining Geol.，2011，10(4):273-287.

[14] L.A.Zadeh.FuzZysets[J].IEEE Information and Control，1965(8):338-353.

[15] 趙鵬大，池順都.初論地質(zhì)異常[J].地球科學(xué)，1991,16(3):241-248.

[16] Agterberg F.P.，Computer programs for mineral exploration[J].Science，1989(245):76-81.

[17] Bonham-Carter GF，Agterberg F.P.，Wright D.F.，Weights of evidence modelling：a new approach to mapping mineral potential.Statistical applications in the earth science[J].geological survey of Canada，1989:171-183.

[18] 王仁鐸.線性地質(zhì)統(tǒng)計(jì)學(xué)[M].北京:中國(guó)地質(zhì)大學(xué)出版社，1985.

[19] 孫英君，王勁峰，柏延臣.地統(tǒng)計(jì)學(xué)方法進(jìn)展研究[J].地球科學(xué)進(jìn)展，2004，19(2):268-274.

[20] 王世稱.綜合信息礦產(chǎn)預(yù)測(cè)理論與方法體系新進(jìn)展[J].地質(zhì)通報(bào)，2010,29(10)：1399-1403.

[21] 陳建平，呂鵬，吳文，等.基于三維可視化技術(shù)的隱伏礦體預(yù)測(cè)[J].地學(xué)前緣，2007，14(5)：55-62.

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

[23] 呂鵬，陳建平，張路鎖，等.基于礦床規(guī)模模型的西南三江北段區(qū)域資源潛力定量預(yù)測(cè)與評(píng)價(jià)[J].地質(zhì)與勘探，2006，42(5):66-71.

[24] 趙志剛，萬(wàn)浩章，董光裕，等.江西貴溪冷水坑銀鉛鋅礦田及外圍構(gòu)造-巖漿-成礦系統(tǒng)解析[J].中國(guó)西部科技，2008,7(30)：4-6.

[25] M.Z.Abzalov，M.Humphreys.Resource Estimation of Structurally Complex and Discontinuous Mineralization Using Non-linear Geostatistics：Case Study of a Mesothermal Gold Deposit in Northern Canada，Explor[J].Mining Geol，2002，11(1-4)：19-29.

[26] T.Claquin，M.Schulz，Y.Balkanski,et al.Uncertainties in assessing radiative forcing by mineral dust[J].Tellus,1998(6):491-505.

[27] 郝天錄.地質(zhì)異常的礦產(chǎn)資源定量化預(yù)測(cè)與不確定性評(píng)價(jià)探討[J].企業(yè)導(dǎo)報(bào)，2011(12):111-112.

[28] 史長(zhǎng)義，趙永平，劉莉.礦產(chǎn)資源地球化學(xué)評(píng)價(jià)和預(yù)測(cè)的幾個(gè)問(wèn)題[J].地質(zhì)與勘探，2003,39(6):14-17.

[29] 左仁廣，夏慶霖.礦產(chǎn)勘查與評(píng)價(jià)中的風(fēng)險(xiǎn)與不確定性[J].礦床地質(zhì)，2010,29(S1):769-770.

[30] 翟裕生，彭潤(rùn)民，向運(yùn)川，等.區(qū)域成礦研究法[M].北京:中國(guó)大地出版社,2004:100-125.

[31] 朱裕生，肖克炎，宋國(guó)耀，等.成礦預(yù)測(cè)方法[M].北京:地質(zhì)出版社,1997.

[32] 張壽庭，趙鵬大，陳建平，等.多目標(biāo)礦產(chǎn)預(yù)測(cè)評(píng)價(jià)及其研究意義[J].成都理工大學(xué)學(xué)報(bào)：自然科學(xué)版,2003,30(5)：441-446.

[33] 趙鵬大，陳建平,張壽庭.“三聯(lián)式”成礦預(yù)測(cè)新進(jìn)展[J].地學(xué)前緣，2003,10(2)：455-463.

Evaluation of information uncertainty in mineral resource prediction and assessment-Take Gan North Lengshuikeng ore field as example

WANG Chun-nv1,2,CHEN Jian-ping2,HAN Zhen-zhe1,2，LI Yong-sheng1，2，JIA Wen-bin2

(1.Development and Research Center of China Geological Survey,Beijing 100037,China; 2.National Exploration &Development Planning Technical Guidance Center in Ministry of Land and Resources，Beijing 100120，China;3.School of Earth Sciences and Resources，China University of Geosciences(Beijing)，Beijing 100083，China)

Mineral prediction uses the geologic metallogenic theory，on the other hand,and combines fully geological，geochemical，geophysical and remote sensing data,and so on.There is a certain degree of uncertainty in the processing and the prediction of information.Therefore，it is important for ore prediction and mining industry economics that estimates correctly uncertainty.This paper analyzes the source of uncertainty in mineral resource prediction and general evaluation method，and proposes a general method of uncertainty evaluation for mineral resources prediction taking Gan North Lengshuikeng ore field as an example，which starts from two aspects of the internal factors and the forecast of geological factor.The methoed get a good result which consistent with the actual geological conditions for five prospective area,therefore，verify the validity of the evaluation method and the feasible of three dimensional prediction method.

mineral resource prediction；uncertainty；Lengshuikeng ore field；method

2014-09-05

重要礦產(chǎn)遠(yuǎn)景調(diào)查區(qū)找礦前景綜合分析與評(píng)價(jià)項(xiàng)目資助(編號(hào)：1212011120997)

王春女(1986-)，女，碩士，地球探測(cè)與信息技術(shù)專業(yè)，從事成礦預(yù)測(cè)方面研究。E-mail：wangchunnv009@163.com。

P612;F407.1

A

1004-4051(2015)05-0142-07