變異函數(shù)及結(jié)構(gòu)分析

變異函數(shù)及結(jié)構(gòu)分析第一頁(yè)，共五十九頁(yè)，編輯于2023年，星期日一、協(xié)方差函數(shù)的計(jì)算公式第一節(jié)協(xié)方差函數(shù)和變異函數(shù)的性質(zhì)設(shè)區(qū)域化變量Z(x)滿足(準(zhǔn))二階平穩(wěn)假設(shè)，h為兩樣本點(diǎn)空間分隔距離，Z(xi)與Z(xi+h)分別是Z(x)在空間位置xi和xi+h上的觀測(cè)值(i=1,2,…N(h)),則計(jì)算協(xié)方差的公式為：第二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第三頁(yè)，共五十九頁(yè)，編輯于2023年，星期日協(xié)方差函數(shù)曲線圖：以h為橫坐標(biāo)，C#(h)為縱坐標(biāo)作圖第四頁(yè)，共五十九頁(yè)，編輯于2023年，星期日二、協(xié)方差函數(shù)的性質(zhì)區(qū)域化變量Z(x)在二階平穩(wěn)假設(shè)下，其協(xié)方差函數(shù)存在且平穩(wěn)，定義為1.C(0)=Var[Z(x)]≥0，即先驗(yàn)方差不能小于零2.C(h)=C(-h)，即C(h)對(duì)h=0的直線對(duì)稱，是一個(gè)偶函數(shù)第五頁(yè)，共五十九頁(yè)，編輯于2023年，星期日證：令x-h=y，則x=y+h，帶入上式得圖形特征及含義第六頁(yè)，共五十九頁(yè)，編輯于2023年，星期日3.|C(h)|≤C(0)，即協(xié)方差函數(shù)絕對(duì)值小于等于先驗(yàn)方差證：第七頁(yè)，共五十九頁(yè)，編輯于2023年，星期日4.|h|→∞時(shí),C(h)→0，或?qū)懽鰿(∞)=0，即當(dāng)空間距離很大時(shí)，協(xié)方差函數(shù)值很小意義(空間局限性)：當(dāng)距離很大時(shí)，Z(x)和Z(x+h)之間的線性相關(guān)基本不存在5.C(h)必須是一個(gè)非負(fù)定函數(shù)，由C(xi-xj）構(gòu)成的協(xié)方差函數(shù)矩陣必須是非負(fù)定矩陣正定條件(positivedefinitecondition)區(qū)域化變量Z(x)二階平穩(wěn)，其數(shù)學(xué)期望為m，協(xié)方差為C(h),變異函數(shù)為γ(h)，令Y是該類型區(qū)域化變量的任意有限線性組合，即：第八頁(yè)，共五十九頁(yè)，編輯于2023年，星期日較難理解則由C(xi-xj）(i,j=1,2…n)構(gòu)成的協(xié)方差函數(shù)矩陣是非負(fù)定矩陣，即C(h)為非負(fù)定函數(shù)二階平穩(wěn)區(qū)域化變量的協(xié)方差函數(shù)是有條件的第九頁(yè)，共五十九頁(yè)，編輯于2023年，星期日三、實(shí)驗(yàn)(經(jīng)驗(yàn))變異函數(shù)(experimentalvariogram)的計(jì)算公式設(shè)區(qū)域化變量Z(x)滿足(準(zhǔn))二階平穩(wěn)條件或(準(zhǔn))本征假設(shè)，h為兩樣本點(diǎn)空間分隔距離，Z(xi)與Z(xi+h)分別是Z(x)在空間位置xi和xi+h上的觀測(cè)值(i=1,2,…N(h)),則計(jì)算實(shí)驗(yàn)變異函數(shù)的公式為：變異函數(shù)曲線圖：以h為橫坐標(biāo)，γ#(h)為縱坐標(biāo)作圖第十頁(yè)，共五十九頁(yè)，編輯于2023年，星期日變異函數(shù)計(jì)算實(shí)例（1）一維變異函數(shù)的計(jì)算x1x2x3x4x5x6x7x8x9x104

3

4

5797877以下為一研究對(duì)象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點(diǎn)間分隔距離h=1米，計(jì)算γ#(h)第十一頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第十二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第十三頁(yè)，共五十九頁(yè)，編輯于2023年，星期日兩方面理解：變異性的理解與相關(guān)性的理解第十四頁(yè)，共五十九頁(yè)，編輯于2023年，星期日作業(yè)：x1x2x3x4x5x6x7x8

2

4

3

15364以下為一研究對(duì)象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點(diǎn)間分隔距離h=1米，計(jì)算γ#(1),γ#(2),γ#(3)第十五頁(yè)，共五十九頁(yè)，編輯于2023年，星期日（2）二維變異函數(shù)的計(jì)算下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，*號(hào)處為無數(shù)據(jù)點(diǎn)，點(diǎn)間距離h為100米，請(qǐng)分別計(jì)算南北、東西、西北和東南方向上的變異函數(shù)值。第十六頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第十七頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第十八頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第十九頁(yè)，共五十九頁(yè)，編輯于2023年，星期日西北和東南方向上的變異函數(shù)值的計(jì)算，注意分隔距離h的確定和樣本數(shù)據(jù)對(duì)的查找第二十頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第二十一頁(yè)，共五十九頁(yè)，編輯于2023年，星期日作業(yè)：下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，網(wǎng)格交叉空白處為無數(shù)據(jù)點(diǎn)，點(diǎn)間距離h為a米，請(qǐng)分別計(jì)算南北方向γ#(a),

西北—東南方向上γ#(a)。第二十二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日四、變異函數(shù)的性質(zhì)區(qū)域化變量Z(x)滿足二階平穩(wěn)或本征假設(shè)條件，則變異函數(shù)存在且平穩(wěn)，計(jì)算公式為1.γ(0)=0，即在h=0時(shí)，變異函數(shù)為零2.γ(h)=γ(-h)，即γ(h)對(duì)h=0的直線對(duì)稱，是一個(gè)偶函數(shù)第二十三頁(yè)，共五十九頁(yè)，編輯于2023年，星期日3.γ(h)≥0,即研究現(xiàn)象的變異性只能大于或等于零4.|h|→∞時(shí),γ(h)→C(0)，或?qū)懽鳓?∞)=C(0)，即當(dāng)空間距離很大時(shí)，變異函數(shù)值接近先驗(yàn)方差5.[-γ(h)]必須是一個(gè)條件非負(fù)定函數(shù)，即由[-γ(xi-xj)]構(gòu)成的變異函數(shù)矩陣必須是條件非負(fù)定矩陣。第二十四頁(yè)，共五十九頁(yè)，編輯于2023年，星期日區(qū)域化變量Z(x)二階平穩(wěn)，其數(shù)學(xué)期望為m，協(xié)方差為C(h),變異函數(shù)為γ(h)，令Y是該類型區(qū)域化變量的任意有限線性組合，即：第二十五頁(yè)，共五十九頁(yè)，編輯于2023年，星期日區(qū)域化變量Z(x)的變異函數(shù)γ(h)是有條件的，即需滿足條件非負(fù)定條件第二十六頁(yè)，共五十九頁(yè)，編輯于2023年，星期日五、協(xié)方差函數(shù)與變異函數(shù)的關(guān)系第二十七頁(yè)，共五十九頁(yè)，編輯于2023年，星期日協(xié)方差函數(shù)和變異函數(shù)的曲線圖問題：為什么只畫出了h>0的關(guān)系圖？當(dāng)h足夠大(即存在a>0，當(dāng)h≥a)時(shí)，可以使C(h)=0,γ(h)=C(0),a稱為變程(range)1、變程a表示區(qū)域化變量從存在空間相關(guān)狀態(tài)(當(dāng)|h|<a時(shí)）轉(zhuǎn)向不存在空間相關(guān)狀態(tài)(當(dāng)|h|>a時(shí))的轉(zhuǎn)折點(diǎn)2、變程a的大小反映區(qū)域化變量影響范圍的大小，或說反映該變量自相關(guān)范圍的大小。也可說變程a是區(qū)域化變量空間變異尺度或空間自相關(guān)尺度變程a的意義：第二十八頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第二節(jié)變異函數(shù)的功能一、變異函數(shù)通過“變程”反映變量的影響范圍——變異函數(shù)的躍遷現(xiàn)象變異函數(shù)γ(h)是一個(gè)單調(diào)遞增函數(shù)，當(dāng)h超過某一數(shù)值(變程a)后，γ(h)不再繼續(xù)單調(diào)地增大，而往往穩(wěn)定在一個(gè)極限值γ(∞)附近，這種現(xiàn)象稱為“躍遷現(xiàn)象”(transitionphenomena)γ(∞)極限值稱為基臺(tái)值(sill),即C(0)【二階平穩(wěn)條件】，基臺(tái)值的大小反映變量變化幅度的大小凡具有一個(gè)變程a和一個(gè)基臺(tái)值的變異函數(shù)，稱為“躍遷型”的變異函數(shù)“變程”反映變量的影響范圍(圖示）第二十九頁(yè)，共五十九頁(yè)，編輯于2023年，星期日二、不同方向上的變異函數(shù)圖可反映區(qū)域化變量的各向異性——變異函數(shù)表示的各向異性如果在各個(gè)方向上區(qū)域化變量的變異性相同或相近，則稱區(qū)域化變量是各向同性的，反之稱為各向異性通過作出各個(gè)方向上的變異函數(shù)圖，并放到一起來比較、分析、研究，就可以確定區(qū)域化變量的各向異性（包括有無各向異性，及各向異性的類型等）第三十頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第三十一頁(yè)，共五十九頁(yè)，編輯于2023年，星期日三、塊金常數(shù)C0的大小可反映區(qū)域化變量的隨機(jī)性大小——變異函數(shù)的塊金效應(yīng)當(dāng)h=0時(shí)，變異函數(shù)γ(h)≠0，而等于一個(gè)常數(shù)C0，這種現(xiàn)象稱為“塊金效應(yīng)”(nuggeteffect),C0稱為塊金常數(shù)或塊金方差(nuggetvariance)塊金效應(yīng)的圖形表示第三十二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日“塊金效應(yīng)”

主要有兩種來源：1、區(qū)域化變量在小于抽樣尺度h時(shí)所具有的變異性2、采樣分析誤差當(dāng)樣點(diǎn)間的距離大于微域結(jié)構(gòu)的范圍，或樣點(diǎn)的大小大于微域結(jié)構(gòu)的范圍就會(huì)出現(xiàn)塊金效應(yīng)（Webester，1985）第三十三頁(yè)，共五十九頁(yè)，編輯于2023年，星期日四、變異函數(shù)在原點(diǎn)處的性狀可反映區(qū)域化變量的空間連續(xù)性變異函數(shù)在原點(diǎn)處的性狀主要有五種類型，每種類型反映了變量的不同程度的空間連續(xù)性

1、拋物線型（parabolictype）當(dāng)|h|→0時(shí),γ(h)→A|h|2(A為常數(shù))，即變異函數(shù)曲線在原點(diǎn)處趨向一條拋物線，反映區(qū)域化變量是具有高度連續(xù)性的，如礦層厚度第三十四頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

2、線性型（lineartype）當(dāng)|h|→0時(shí),γ(h)→A|h|(A為常數(shù))，即變異函數(shù)曲線在原點(diǎn)處趨向一條直線，或說在原點(diǎn)處有斜向的切線存在，反映區(qū)域化變量是具有平均的連續(xù)性，如金屬品位

3、間斷型（discontinuoustype）當(dāng)|h|→0時(shí),γ(h)→C0，即變異函數(shù)曲線在原點(diǎn)處間斷，說明塊金效應(yīng)存在，又稱“塊金效應(yīng)型”，反映區(qū)域化變量的連續(xù)性很差，但當(dāng)h增大時(shí)，γ(h)又變的較為連續(xù)了，如金品位第三十五頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

4、隨機(jī)型（randomtype）這種變異函數(shù)可看成具有基臺(tái)值C0和無窮小變程a的躍遷型變異函數(shù)，則無論h多小，h總大于a，故Z(x)與Z(x+h)總是互不相關(guān)又稱純塊金效應(yīng)型，反映了區(qū)域化變量完全不存在空間相關(guān)的情況，則本質(zhì)上此區(qū)域化變量為普通隨機(jī)變量此時(shí)，C0=C(0)第三十六頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

5、過渡型：介于拋物線型和隨機(jī)型間當(dāng)|h|→0時(shí),γ(h)→C0，即有塊金效應(yīng)；當(dāng)|h|=a時(shí),γ(a)=C(0)，即有基臺(tái)值(C0+C)和變程a,C稱為“拱高”過渡型是實(shí)際研究工作中最常遇到的一種類型第三十七頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第三節(jié)變異函數(shù)的理論模型

思考：是否有了采樣數(shù)據(jù)及變異函數(shù)計(jì)算公式就可以獲知任意距離h的區(qū)域化變量變異性？設(shè)Z(x)具有各向同性的變異函數(shù)

γ(h)，則常見的變異函數(shù)模型如下：變異函數(shù)的理論模型有基臺(tái)值模型無基臺(tái)值模型可以有或無基臺(tái)值模型：孔穴效應(yīng)模型球狀模型、指數(shù)模型高斯模型線性有基臺(tái)模型純塊金效應(yīng)模型冪函數(shù)模型對(duì)數(shù)模型線性無基臺(tái)模型第三十八頁(yè)，共五十九頁(yè)，編輯于2023年，星期日一、有基臺(tái)值模型

1、球狀模型（sphericalmodel）若模型滿足二階平穩(wěn)假設(shè)，且有有限先驗(yàn)方差，γ(h)值隨h的變大而增大，當(dāng)h達(dá)一定值(h>a)時(shí)，γ(h)達(dá)到一定值——基臺(tái)值，則稱此類模型為有基臺(tái)值模型式中：C0為塊金常數(shù)，(C0+C)為基臺(tái)值，C為拱高，a為變程第三十九頁(yè)，共五十九頁(yè)，編輯于2023年，星期日當(dāng)C0=0，C=1，稱為標(biāo)準(zhǔn)球狀模型，其圖形為：原點(diǎn)處切線的斜率為3/2a,與基臺(tái)值線交點(diǎn)的橫坐標(biāo)為2a/3球狀模型是地統(tǒng)計(jì)學(xué)應(yīng)用最廣的理論模型，許多區(qū)域化變量的理論模型都可以用球狀模型來擬合第四十頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

2、指數(shù)模型（exponentialmodel）式中：C0，C意義同前，但a不是變程當(dāng)C0=0，C=1，稱為標(biāo)準(zhǔn)指數(shù)模型，其圖形為：由于1-e-3=1-0.05=0.95≈1，則變程為3a第四十一頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

3、高斯模型（gaussianmodel）式中：C0，C意義同前，但a不是變程由于1-e-3=1-0.05=0.95≈1，則變程為√3

a當(dāng)C0=0，C=1，稱為標(biāo)準(zhǔn)高斯函數(shù)模型，其圖形為：第四十二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日模型通過原點(diǎn)切線與基臺(tái)值線交點(diǎn)的橫坐標(biāo)變程原點(diǎn)處的性狀球狀2a/3a直線指數(shù)a3a直線高斯無交點(diǎn)√3

a拋物線三種模型的比較第四十三頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

4、線性有基臺(tái)值模型（linearwithsillmodel）式中：C0，C意義同前，A為常數(shù)，表示直線的斜率，變程為a第四十四頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

5、純塊金效應(yīng)模型（purenuggeteffectmodel）此時(shí)，C0=C(0)此種模型意味著區(qū)域化變量為隨機(jī)分布，樣點(diǎn)間的協(xié)方差函數(shù)對(duì)于所有距離h均等于0，即變量不存在空間相關(guān)性第四十五頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第四十六頁(yè)，共五十九頁(yè)，編輯于2023年，星期日二、無基臺(tái)值模型

1、冪函數(shù)模型（powermodel）若與模型相應(yīng)的區(qū)域化變量不滿足二階平穩(wěn)假設(shè)，僅滿足本征假設(shè)，γ(h)值隨h的變大而增大，但不能達(dá)到一定值，即無基臺(tái)值，則稱此類模型為無基臺(tái)值模型當(dāng)改變參數(shù)θ時(shí)，可以表示原點(diǎn)處的各種性狀第四十七頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

2、線性無基臺(tái)值模型（linearwithoutsillmodel）第四十八頁(yè)，共五十九頁(yè)，編輯于2023年，星期日

3、對(duì)數(shù)模型（powermodel）第四十九頁(yè)，共五十九頁(yè)，編輯于2023年，星期日三、孔穴效應(yīng)模型(holeeffectmodel)當(dāng)變異函數(shù)γ(h)在大于一定距離后，并非單調(diào)遞增，而具有一定周期波動(dòng)，此種模型稱為孔穴效應(yīng)模型有基臺(tái)值無基臺(tái)值第五十頁(yè)，共五十九頁(yè)，編輯于2023年，星期日第四節(jié)變異函數(shù)的結(jié)構(gòu)分析一、結(jié)構(gòu)分析、套合結(jié)構(gòu)概念采樣數(shù)據(jù)計(jì)算γ#(h)試驗(yàn)變異函數(shù)曲線對(duì)區(qū)域化變量進(jìn)行分析合適的理論模型實(shí)際中區(qū)域化變量的變化性很復(fù)雜：(1)可能在不同方向上有不同的變異性；(2)在同一方向上包含不同尺度上的多層次的變異性這么復(fù)雜?。康谖迨豁?yè)，共五十九頁(yè)，編輯于2023年，星期日礦床或礦體的變異性往往由多種原因引起采樣、樣品制備及分析等過程所產(chǎn)生的誤差原因礦物成分的變化，如金礦等品位變化劇烈的礦床上尤為明顯礦層與夾層的交替變化礦床分布引起的變異≈01～ncm米至百米公里尺度不同原因引起的變異特性，其變異尺度的大小不同顯然，大尺度的變異總是包含著小尺度的變異，小尺度的變異在大尺度變異曲線上只能作為“塊金效應(yīng)”出現(xiàn)第五十二頁(yè)，共五十九頁(yè)，編輯于2023年，星期日土壤的空間變