一種地質(zhì)儲量估算方法——點(diǎn)克里格估值法編程實(shí)現(xiàn)探究(圖文)

1、一種地質(zhì)儲量估算方法一種地質(zhì)儲量估算方法——點(diǎn)克里格估值法編程實(shí)現(xiàn)點(diǎn)克里格估值法編程實(shí)現(xiàn)探究探究(圖文圖文)克立格法是由法國G馬特隆教授以南非礦山地質(zhì)工程師DG克立格的名字命名的，是一種對空間分布數(shù)據(jù)求最優(yōu)、線性、無偏內(nèi)插估計(jì)量的方法?？死锔穹ㄊ堑刭|(zhì)統(tǒng)計(jì)學(xué)（空間信息統(tǒng)計(jì)學(xué)）的核心，最早應(yīng)用于礦山開采。目前地質(zhì)統(tǒng)計(jì)學(xué)在理論和軟件方面發(fā)展迅速。理論方面：Journel將克立格法的估值問題，從一般矢量空間擴(kuò)展到n個(gè)原始數(shù)據(jù)的全部可測度函數(shù)所形

2、成的矢量空間（希爾伯特空間）進(jìn)行考察。為地統(tǒng)計(jì)學(xué)由低維向高維動態(tài)發(fā)展奠定了基礎(chǔ)。Matheron為首的參數(shù)地質(zhì)統(tǒng)計(jì)學(xué)派以正態(tài)假設(shè)為前提，在協(xié)同區(qū)域化理論的基礎(chǔ)上，提出多元地質(zhì)統(tǒng)計(jì)學(xué)的基本思想。Journel發(fā)展了無須對數(shù)據(jù)分布作任何假設(shè)的非參數(shù)地質(zhì)統(tǒng)計(jì)學(xué)，提出了一些非參數(shù)地質(zhì)統(tǒng)計(jì)學(xué)克立格方法。同時(shí)，多學(xué)科相互滲透，拓展了地質(zhì)統(tǒng)計(jì)學(xué)的理論體系，使地統(tǒng)計(jì)學(xué)呈現(xiàn)出以下特點(diǎn)：地質(zhì)統(tǒng)計(jì)學(xué)與數(shù)學(xué)的結(jié)合更加緊密；非線性科學(xué)思想向地質(zhì)統(tǒng)計(jì)學(xué)滲透；穩(wěn)健統(tǒng)

3、計(jì)學(xué)與地質(zhì)統(tǒng)計(jì)學(xué)的結(jié)合更加密切。軟件方面：國外的空間信息統(tǒng)計(jì)學(xué)軟件主要有：法國巴黎高等礦院地質(zhì)統(tǒng)計(jì)學(xué)研究中心研制出兩種大型軟件系統(tǒng)：ISATIS系統(tǒng)及HERESIM系統(tǒng)。這兩個(gè)系統(tǒng)對油氣田地質(zhì)及開采工程等方面考慮得尤為周到。美國斯坦福大學(xué)地球科學(xué)應(yīng)用系研制的GEOEAS、GEOSTATTOOLBOX程序包及1922年、1996年又推出了GSLIB相應(yīng)的理論變差函數(shù)模型參與Kriging估值。地統(tǒng)計(jì)學(xué)將變異函數(shù)理論模型分為三大類：第一類是

4、有基臺值模型，包括球狀模型、指數(shù)模型、高斯模型、線性有基臺值模型和純塊金效應(yīng)模型；第二類是無基臺值模型，包括冪函數(shù)模型、線性無基臺值模型、拋物線模型；第三類是孔穴效應(yīng)模型。球狀模型是地統(tǒng)計(jì)分析中應(yīng)用最廣泛的理論模型，本文采用球狀模型進(jìn)行變差函數(shù)擬合。3軟件設(shè)計(jì)本軟件系統(tǒng)將總的功能實(shí)現(xiàn)分解為原始數(shù)據(jù)輸入、理論變差函數(shù)計(jì)算、變差函數(shù)擬合、擬合參數(shù)調(diào)整、克里格估值、軟件操作幫助六個(gè)模塊，如圖1所示。在理論變差函數(shù)計(jì)算模塊中，將存入預(yù)定義數(shù)組中