基于數(shù)學(xué)理論的成礦預(yù)測(cè)方法研究

基于數(shù)學(xué)理論的成礦預(yù)測(cè)方法研究

0成礦預(yù)測(cè)中統(tǒng)計(jì)理論的提出及應(yīng)用必要性建立數(shù)學(xué)模型，預(yù)測(cè)成礦是探索礦位的重要方法。利用統(tǒng)計(jì)方法對(duì)已有數(shù)據(jù)進(jìn)行處理,進(jìn)而推斷其內(nèi)部規(guī)律,并據(jù)此進(jìn)行隱伏礦體定位仍是成礦預(yù)測(cè)的主流。常用的統(tǒng)計(jì)方法包括多元統(tǒng)計(jì)、地質(zhì)統(tǒng)計(jì)學(xué)等。近年來(lái),既有傳統(tǒng)方法的更新,如新的統(tǒng)計(jì)分布模型的應(yīng)用、克立格方法的擴(kuò)展等,也有新理論和新方法的引入,如統(tǒng)計(jì)分形、隨機(jī)過(guò)程等。本文從數(shù)學(xué)理論分析角度入手,討論了成礦預(yù)測(cè)中統(tǒng)計(jì)理論的數(shù)學(xué)表述及其物理意義,并探討了它們之間的內(nèi)在聯(lián)系。同時(shí),提出了Weibull分布模型是成礦元素分布的另一種更為普遍的分布形式,并將Weibull模型與分形相結(jié)合對(duì)膠東礦集區(qū)大尹格莊金礦的元素分布進(jìn)行了應(yīng)用研究。1概率統(tǒng)計(jì)1.1概率分布特征數(shù)據(jù)概率分布特征分析是通過(guò)有限或局部的地質(zhì)數(shù)據(jù)推斷整體的分布特征,進(jìn)而開(kāi)展成礦預(yù)測(cè)研究的重要步驟。對(duì)數(shù)據(jù)概率分布的分析是通過(guò)分布函數(shù)來(lái)描述的,分布函數(shù)既能說(shuō)明隨機(jī)變量的可能取值,又能指出將以多大的概率取這些值,能較完整地描述隨機(jī)變量的變化特征。因此,分析觀測(cè)數(shù)據(jù)并確定研究對(duì)象的概率分布函數(shù)是一項(xiàng)十分重要的工作。分析觀測(cè)數(shù)據(jù)分布特征的最常用方法是作直方圖和統(tǒng)計(jì)分布曲線圖。在實(shí)際應(yīng)用中多使用正態(tài)分布、對(duì)數(shù)正態(tài)分布或指數(shù)分布來(lái)描述成礦元素的品位分布,但是上述數(shù)學(xué)模型的擬合效果與實(shí)際地質(zhì)數(shù)據(jù)分布之間仍存在一定差距,并且在描述不同元素分布的過(guò)程中,常需變換不同的數(shù)學(xué)模型。Weibull模型的分布曲線對(duì)單調(diào)減、單峰對(duì)稱和單峰偏斜等分布具有靈活的適應(yīng)性,由此推測(cè)Weibull模型可能更適于描述成礦元素分布。在某礦區(qū)成礦元素的應(yīng)用中發(fā)現(xiàn),Weibull分布模型更接近于元素分布,該模型的擬合曲線與原始數(shù)據(jù)直方圖的對(duì)應(yīng)效果較好。Weibull分布模型的概率密度函數(shù)可表示為:f(x)=αβ(xβ)α-1exp[-(xβ)α]?x＞0(1)f(x)=αβ(xβ)α?1exp[?(xβ)α]?x＞0(1)其中:形狀參數(shù)α>0;尺度參數(shù)β>0。它們控制密度函數(shù)曲線的二維形態(tài),可以描述多種對(duì)象。1.2多元統(tǒng)計(jì)分析1.2.1利用內(nèi)在聯(lián)系生成主要成分?jǐn)?shù)據(jù)屬性分析是要判別各組數(shù)據(jù)的關(guān)系,把具有相似性的數(shù)據(jù)進(jìn)行分類;將存在復(fù)雜關(guān)系的較多變量依據(jù)某種內(nèi)在聯(lián)系作最佳組合,生成幾個(gè)主要的成分等。數(shù)據(jù)屬性的分析方法主要有相關(guān)分析與回歸分析、點(diǎn)群分析、判別分析、典型相關(guān)分析、因子分析、非線性映射分析和對(duì)應(yīng)分析等。1.2.2進(jìn)行趨勢(shì)分析由于成礦預(yù)測(cè)觀測(cè)值可分解為趨勢(shì)部分和隨機(jī)部分,對(duì)觀測(cè)數(shù)據(jù)進(jìn)行趨勢(shì)分析是要將趨勢(shì)部分分離出來(lái)。采用的方法是用空間或時(shí)間的某種函數(shù)來(lái)逼近趨勢(shì),從而分解出趨勢(shì)部分。主要方法包括趨勢(shì)面分析、調(diào)和趨勢(shì)分析、典型趨勢(shì)分析和單位向量場(chǎng)分析。2可變異性函數(shù)地質(zhì)統(tǒng)計(jì)學(xué)是以區(qū)域化變量理論為基礎(chǔ),以變異函數(shù)為基本工具,研究在時(shí)空分布上既具有隨機(jī)性又具有結(jié)構(gòu)性的自然現(xiàn)象的科學(xué)。變異函數(shù)是研究時(shí)空分布數(shù)據(jù)的結(jié)構(gòu)性、變異性的有效工具。在成礦預(yù)測(cè)中的許多變量,如礦石品位、礦體厚度、巖礦的各種物化參數(shù)、地形的標(biāo)高等都和觀測(cè)點(diǎn)的坐標(biāo)有關(guān),因而都是區(qū)域化變量,都可以用變異函數(shù)研究。地質(zhì)統(tǒng)計(jì)分析在成礦數(shù)據(jù)分析中的應(yīng)用非常廣泛,而且歷史悠久,它已經(jīng)成為地質(zhì)研究工作中必不可少的應(yīng)用工具,其理論和方法都將隨著地質(zhì)研究中新問(wèn)題的出現(xiàn)而不斷發(fā)展和完善。2.1a點(diǎn)x的變異函數(shù)區(qū)域化變量是空間上的一個(gè)數(shù)值函數(shù),它在空間的每一個(gè)點(diǎn)取一個(gè)確定值。由于區(qū)域化變量具有結(jié)構(gòu)性和隨機(jī)性,因此不能用通常的數(shù)學(xué)分析方法對(duì)這種數(shù)值空間函數(shù)進(jìn)行直接的研究,為此需作平穩(wěn)假設(shè):(1)假定的局部范圍內(nèi),變量的數(shù)學(xué)期望值為一常數(shù),并且不依賴于支撐點(diǎn)x,即E[Ζ(x)]=m?x(2)E[Z(x)]=m?x(2)(2)任意一對(duì)隨機(jī)變量值{Z(x),Z(x+h)}的協(xié)方差存在,并且只依賴于相隔的距離h,即C(h)=E{Ζ(x+h)?Ζ(x)}-m2(3)C(h)=E{Z(x+h)?Z(x)}?m2(3)變異函數(shù)是任一方向α、相距為h的兩個(gè)區(qū)域變量值Z(x)和Z(x+h)增量的方差的一半,記為2γ(h,α)=Var{Ζ(x)-Ζ(x+h)}(4)2γ(h,α)=Var{Z(x)?Z(x+h)}(4)由實(shí)測(cè)數(shù)據(jù)獲得的變異函數(shù)稱為先驗(yàn)變異函數(shù),記為γ*(h)=12Ν(h)Ν(h)∑i=1[Ζ(xi)-Ζ(xi+h)]2(5)其中要求Z(x)存在且平穩(wěn),N(h)是被矢量分隔的數(shù)據(jù)對(duì)的數(shù)目。利用變異函數(shù)能有效分析成礦預(yù)測(cè)變量的不同連續(xù)性、不同類型的各向異性和不同程度的可遷性。2.2金融學(xué)變量范圍的設(shè)定克立格法是一種局部估計(jì)方法,它以最小的估計(jì)方差給出對(duì)于塊段平均品位的無(wú)偏估計(jì)量。實(shí)際上,它也是一種廣義的最小二乘回歸算法,而其最優(yōu)目標(biāo)是誤差的期望值為零,方差達(dá)到最小。根據(jù)不同的研究目的和不同的約束條件,有不同的克立格法,但都是通過(guò)先確立克立格方程組,然后再求解的。根據(jù)所研究對(duì)象的不同趨勢(shì),克里格法可有3種不同類型:(1)簡(jiǎn)單克立格法,E[Z(x)]=m為已知,且Z(x)滿足二階平穩(wěn)(或內(nèi)蘊(yùn))假設(shè);(2)普通克立格法,E[Z(x)]=m為未知,且Z(x)滿足二階平穩(wěn)(或內(nèi)蘊(yùn))假設(shè);(3)泛克立格法,E[Z(x)]=m(x),Z(x)非平穩(wěn),其中Z(x)=Y(x)+m(x),Y(x)為漂移,且E{Y(x)}=0。依據(jù)研究目的不同,克立格法還可以包括:(1)因子克立格法,用以理解各種不同尺度影響的根源;(2)協(xié)同克立格法,用于多個(gè)變量且具有協(xié)同區(qū)域性的情況;(3)塊狀克立格法,用特定的尺度“塊”——稱為支集的研究;(4)對(duì)數(shù)克立格法,用于區(qū)域變量Z(x)服從對(duì)數(shù)正態(tài)分布情況;(5)指示克立格法,對(duì)有特異值的數(shù)據(jù)進(jìn)行分析。綜上所述,地質(zhì)統(tǒng)計(jì)學(xué)的核心思想是通過(guò)變異函數(shù),把所有的點(diǎn)對(duì)按照間隔距離的大小、方向進(jìn)行分組,在每一個(gè)組內(nèi),計(jì)算每個(gè)點(diǎn)對(duì)屬性值的差異,取其平均作為該組屬性值的差異(變異值)。這樣,將整個(gè)空間分為不同大小和方向的組,并有相對(duì)應(yīng)的屬性差值。然后,再以線性、無(wú)偏和最小估計(jì)方差為前提,利用克立格方程組的解,對(duì)每一樣品值分別賦給一定的權(quán)系數(shù),最后進(jìn)行加權(quán)平均來(lái)進(jìn)行估計(jì)。隨著研究的不斷深入,數(shù)據(jù)空間結(jié)構(gòu)進(jìn)一步細(xì)化,對(duì)估計(jì)量的精度要求會(huì)越來(lái)越高,因此,還會(huì)有新的克立格法不斷產(chǎn)生。3分形的內(nèi)涵、測(cè)量分形是其組成部分以某種方式與整體相似的形,是研究自然界空間結(jié)構(gòu)復(fù)雜性的一門(mén)學(xué)科。分形理論中最核心的概念是分維數(shù),它是分形的定量表征。研究表明,通過(guò)地質(zhì)體的分維值可以挖掘潛在的成礦信息,建立分形模型。自然界中的分形有兩種基本類型:統(tǒng)計(jì)自相似分形和統(tǒng)計(jì)自仿射分形。3.1特征尺度及模型統(tǒng)計(jì)自相似分形具有不隨標(biāo)度變化的自相似特征,當(dāng)其局部被放大時(shí),它表現(xiàn)為與整體統(tǒng)計(jì)相似。也可以說(shuō)自相似分形具有各向同性的特點(diǎn),即在由x和y坐標(biāo)所確定的二維情況下,結(jié)構(gòu)與x軸和y軸的幾何取向無(wú)關(guān)。常用的統(tǒng)計(jì)分形模型為:Ν(r)=Cr-Dr＞0(6)其中:r表示特征尺度;C>0,稱為比例常數(shù);D>0,稱為分維數(shù);N(r)=Cr-D,表示尺度大于等于r的數(shù)目,記為N(≥r)。分維數(shù)的求值有兩種方法:一元線性回歸方法和非線性回歸方法。一元線性回歸法是將觀測(cè)數(shù)據(jù)(N(r1),N(r2),…,N(rn))和(r1,r2,…,rn)代入線性回歸方程:lgΝ(r)=-Dlgr+lgC(7)用最小二乘法求出斜率D,即分維值的估計(jì)值。非線性回歸法是以(6)式為非線性回歸模型,用非線性回歸的最小二乘法直接求出參數(shù)的估計(jì)量,也就是分維數(shù)。此方法求出的分維數(shù)比用一元線性回歸模型求出的分維數(shù)更精確,即誤差更小,而且不必受取對(duì)數(shù)的條件所限制。3.2lq1h及tt-lt特征與統(tǒng)計(jì)自相似分形不同的是,統(tǒng)計(jì)自仿射分形自相似性不單一,隨著尺度的不同而變化,具有各向異性的特點(diǎn)。其描述整體行為的分形維數(shù)與描述局部行為的分形維數(shù)不相等,也就是說(shuō)它具有與自相似分形完全不同的特點(diǎn)。如分?jǐn)?shù)Brown運(yùn)動(dòng),記為BH(t),是時(shí)間段[0,T]上的函數(shù),H是Hausdorff測(cè)度,且H∈,BH(t)具有統(tǒng)計(jì)特征:σ(BΗ(t))∝ΤΗ(8)σ(BΗ(bt))∝bΗΤΗ(9)其中b是伸縮常數(shù)。根據(jù)(8)、(9)兩式可得σ(BH(bt))∝bHσ(BH(t)),這說(shuō)明分?jǐn)?shù)Brown運(yùn)動(dòng)是自仿射統(tǒng)計(jì)分形。利用數(shù)盒子方法,可求得分?jǐn)?shù)Brown運(yùn)動(dòng)的局部維是D=2-H,而它的整體維當(dāng)b→∞時(shí)是D=1。實(shí)際上,現(xiàn)實(shí)中地質(zhì)問(wèn)題往往表現(xiàn)為自仿射分形,如多數(shù)地形、地下深部構(gòu)造、Brown運(yùn)動(dòng)曲線等都是自仿射統(tǒng)計(jì)分形。可以預(yù)測(cè),自仿射分形理論在地質(zhì)學(xué)研究中有廣泛的應(yīng)用前景。自仿射分形在地質(zhì)研究中的應(yīng)用還不像自相似分形的應(yīng)用那么普及,主要原因是由于自仿射分形的分維數(shù)不易求出,因?yàn)樽苑律浞中尉S數(shù)是整體維還是局部維與所選尺度單位有關(guān)。Mandelbrot的研究認(rèn)為:對(duì)于自仿射分形,存在一個(gè)新的量,稱為跨越長(zhǎng)度tc(crossoverlength),在它兩側(cè)存在著局域和整體兩個(gè)不同的分維值。而合理地選擇跨越長(zhǎng)度tc是十分困難的。4受系統(tǒng)本身存在的概率有隨機(jī)過(guò)程所描述的是自變量和因變量的關(guān)系不確定的物體的運(yùn)動(dòng)過(guò)程,自變量和因變量間既存在某種相關(guān)關(guān)系,又受系統(tǒng)本身存在的特定概率的支配。在數(shù)學(xué)上,將依賴于一個(gè)參數(shù)t的隨機(jī)變量集合{xt}稱為隨機(jī)過(guò)程。隨機(jī)過(guò)程的種類很多,按照過(guò)程的概率結(jié)構(gòu)可分為二階矩過(guò)程、平穩(wěn)過(guò)程、馬爾可夫過(guò)程、鞅過(guò)程和隨機(jī)點(diǎn)過(guò)程等。下面只介紹在地學(xué)數(shù)據(jù)處理中應(yīng)用較多和有應(yīng)用潛力的馬爾可夫過(guò)程和隨機(jī)點(diǎn)過(guò)程兩種方法。4.1成礦過(guò)程的狀態(tài)轉(zhuǎn)移馬爾可夫過(guò)程是研究事物狀態(tài)及狀態(tài)轉(zhuǎn)移的理論,該過(guò)程考慮了以前事件對(duì)后來(lái)事件的影響,即從一種狀態(tài)轉(zhuǎn)移到另一種狀態(tài),隨時(shí)間變化所作的狀態(tài)轉(zhuǎn)移,且狀態(tài)轉(zhuǎn)移具有概率性質(zhì)。它是以時(shí)間序列來(lái)處理觀測(cè)數(shù)據(jù)的。若馬爾可夫過(guò)程的狀態(tài)是離散的,則稱此過(guò)程為馬爾可夫鏈。在成礦預(yù)測(cè)中,元素豐度變化可看成是一個(gè)隨機(jī)運(yùn)動(dòng)著的狀態(tài)變化。由于狀態(tài)的轉(zhuǎn)移是一隨機(jī)事件,在研究中雖然不能準(zhǔn)確地預(yù)言下一時(shí)刻的豐度值,但可以根據(jù)馬爾可夫鏈理論和方法,事先給出轉(zhuǎn)移到某一狀態(tài)的概率,即轉(zhuǎn)移概率,從而達(dá)到研究和預(yù)測(cè)的目的。4.2預(yù)測(cè)的過(guò)程屬于序列的疊加過(guò)程隨機(jī)點(diǎn)過(guò)程是隨機(jī)過(guò)程理論的一個(gè)重要分支,它是依賴于循環(huán)時(shí)間的隨機(jī)過(guò)程。用數(shù)學(xué)術(shù)語(yǔ)表達(dá)為:設(shè){ξn}n=0,1,2,…是一個(gè)離散參數(shù)隨機(jī)過(guò)程,對(duì)于所有n,存在正的常數(shù)M,并且有M≥ξn≥0,若取tn=n∑i=1ξi,則稱{tn}是一個(gè)隨機(jī)點(diǎn)過(guò)程。實(shí)際上,它是一個(gè)序列的累加過(guò)程。在成礦預(yù)測(cè)中,可用其分析沉積作用,因?yàn)槌练e物的厚度可看成是一定重復(fù)時(shí)間的累加過(guò)程。在地學(xué)中,沉積地層厚度和剖面元素的分布等系列數(shù)據(jù)均可視為隨機(jī)點(diǎn)過(guò)程。Huang等從隨機(jī)動(dòng)力學(xué)角度,在厘米級(jí)的巖性變化尺度上,定量分析了鄂東二疊系大隆組頂部和三疊系大冶組底部的沉積行為及動(dòng)力學(xué)特征,詳細(xì)研究了沉積過(guò)程中的馬爾可夫過(guò)程,為該區(qū)的P—T之交的沉積環(huán)境突變,給出了動(dòng)力學(xué)解釋。5多種方法相互配合的方法本文從數(shù)學(xué)理論分析角度出發(fā),討論了多元統(tǒng)計(jì)、地質(zhì)統(tǒng)計(jì)分析、統(tǒng)計(jì)分形和隨機(jī)過(guò)程等多種方法的數(shù)學(xué)表述及其物理內(nèi)涵;這些方法與其分支之間存在內(nèi)在的聯(lián)系(圖1)。從圖1中可以看出,概率統(tǒng)計(jì)是多元統(tǒng)計(jì)、地質(zhì)統(tǒng)計(jì)分析、統(tǒng)計(jì)分形和隨機(jī)過(guò)程的公共基礎(chǔ)。多元統(tǒng)計(jì)、地質(zhì)統(tǒng)計(jì)分析、統(tǒng)計(jì)分形和隨機(jī)過(guò)程從不同的角度描述了數(shù)據(jù)集的物理內(nèi)涵,如多元統(tǒng)計(jì)學(xué)的十幾種方法,可概括為相依性、趨勢(shì)性與依賴型。多種傳統(tǒng)方法相互配合使用,可以進(jìn)行數(shù)據(jù)集的深度挖掘。以分形為例,分形規(guī)律與上述分析方法有諸多聯(lián)系之處。如在Weibull模型中,當(dāng)α?β時(shí),表現(xiàn)為自相似分形,其形狀參數(shù)α對(duì)應(yīng)于數(shù)據(jù)集分維值;區(qū)域化變量Z(x)滿足二階平穩(wěn)假設(shè)時(shí),變異函數(shù)r(h)與間距h存在冪律關(guān)系,即r(h)=ch4-2D,其中D是分維值;隨機(jī)行走是自仿射分形的典型實(shí)例。筆者將Weibull分布和分形模型有機(jī)地結(jié)合,分析了膠東礦集區(qū)大尹格莊金礦中S、Cu、Pb、Au4種元素的空間分布特征。首先將4種元素的空間分布進(jìn)行分形統(tǒng)計(jì),得出盒子維數(shù);再以盒子維數(shù)、均值作為Weibull模型的形狀參數(shù)α及尺度參數(shù)β,做出統(tǒng)計(jì)對(duì)象的分布曲線;進(jìn)一步將Weibull曲線與元素空間分布的直方圖進(jìn)行對(duì)比,判斷兩者是否吻合。S、Cu、Pb、Au4種元素的盒子維數(shù)各不相同,直方圖也互有差異(圖2至圖5),將盒子維數(shù)作為Weibull模型的形狀參數(shù),做出Weibull分布曲線,其與元素含量直方圖擬合較好。Weibull模型的成功應(yīng)用說(shuō)