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

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

3

4

5797877以下為一研究對象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點間分隔距離h=1米，計算γ#(h)變異函數(shù)計算實例（1）一維變異函數(shù)的計算x1x211變異函數(shù)及結(jié)構(gòu)分析課件12變異函數(shù)及結(jié)構(gòu)分析課件13兩方面理解：變異性的理解與相關(guān)性的理解兩方面理解：變異性的理解與相關(guān)性的理解14作業(yè)：x1x2x3x4x5x6x7x8

2

4

3

15364以下為一研究對象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點間分隔距離h=1米，計算γ#(1),γ#(2),γ#(3)作業(yè)：x1x2x3x4x515（2）二維變異函數(shù)的計算下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，*號處為無數(shù)據(jù)點，點間距離h為100米，請分別計算南北、東西、西北和東南方向上的變異函數(shù)值。（2）二維變異函數(shù)的計算下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，*號處16變異函數(shù)及結(jié)構(gòu)分析課件17變異函數(shù)及結(jié)構(gòu)分析課件18變異函數(shù)及結(jié)構(gòu)分析課件19西北和東南方向上的變異函數(shù)值的計算，注意分隔距離h的確定和樣本數(shù)據(jù)對的查找西北和東南方向上的變異函數(shù)值的計算，注意分隔距離h的確定和樣20變異函數(shù)及結(jié)構(gòu)分析課件21作業(yè)：下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，網(wǎng)格交叉空白處為無數(shù)據(jù)點，點間距離h為a米，請分別計算南北方向γ#(a),

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

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

1、拋物線型（parabolictype）當(dāng)|h|→0時,γ(h)→A|h|2(A為常數(shù))，即變異函數(shù)曲線在原點處趨向一條拋物線，反映區(qū)域化變量是具有高度連續(xù)性的，如礦層厚度四、變異函數(shù)在原點處的性狀可反映區(qū)域化變量的空間連續(xù)性變34

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

3、間斷型（discontinuoustype）當(dāng)|h|→0時,γ(h)→C0，即變異函數(shù)曲線在原點處間斷，說明塊金效應(yīng)存在，又稱“塊金效應(yīng)型”，反映區(qū)域化變量的連續(xù)性很差，但當(dāng)h增大時，γ(h)又變的較為連續(xù)了，如金品位2、線性型（lineartype）當(dāng)|h|→0時,γ35

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

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

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

γ(h)，則常見的變異函數(shù)模型如下：變異函數(shù)的理論模型有基臺值模型無基臺值模型可以有或無基臺值模型：孔穴效應(yīng)模型球狀模型、指數(shù)模型高斯模型線性有基臺模型純塊金效應(yīng)模型冪函數(shù)模型對數(shù)模型線性無基臺模型第三節(jié)變異函數(shù)的理論模型思考：是否有了采樣數(shù)據(jù)及變異函38一、有基臺值模型

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

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，則變程為3a2、指數(shù)模型（exponentialmodel）式中：41

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ù)模型，其圖形為：3、高斯模型（gaussianmodel）式中：C0，42模型通過原點切線與基臺值線交點的橫坐標(biāo)變程原點處的性狀球狀2a/3a直線指數(shù)a3a直線高斯無交點√3

a拋物線三種模型的比較模型通過原點切線與基臺值線交點的橫坐標(biāo)變程原點處的性狀球狀243

4、線性有基臺值模型（linearwithsillmodel）式中：C0，C意義同前，A為常數(shù)，表示直線的斜率，變程為a4、線性有基臺值模型（linearwithsillm44

5、純塊金效應(yīng)模型（purenuggeteffectmodel）此時，C0=C(0)此種模型意味著區(qū)域化變量為隨機(jī)分布，樣點間的協(xié)方差函數(shù)對于所有距離h均等于0，即變量不存在空間相關(guān)性5、純塊金效應(yīng)模型（purenuggeteffect45變異函數(shù)及結(jié)構(gòu)分析課件46二、無基臺值模型

1、冪函數(shù)模型（powermodel）若與模型相應(yīng)的區(qū)域化變量不滿足二階平穩(wěn)假設(shè)，僅滿足本征假設(shè)，γ(h)值隨h的變大而增大，但不能達(dá)到一定值，即無基臺值，則稱此類模型為無基臺值模型當(dāng)改變參數(shù)θ時，可以表示原點處的各種性狀二、無基臺值模型1、冪函數(shù)模型（powermodel）47

2、線性無基臺值模型（linearwithoutsillmodel）2、線性無基臺值模型（linearwithoutsil48

3、對數(shù)模型（powermodel）3、對數(shù)模型（powermodel）49三、孔穴效應(yīng)模型(holeeffectmodel)當(dāng)變異函數(shù)γ(h)在大于一定距離后，并非單調(diào)遞增，而具有一定周期波動，此種模型稱為孔穴效應(yīng)模型有基臺值無基臺值三、孔穴效應(yīng)模型(holeeffectmodel)當(dāng)變50第四節(jié)變異函數(shù)的結(jié)構(gòu)分析一、結(jié)構(gòu)分析、套合結(jié)構(gòu)概念采樣數(shù)據(jù)計算γ#(h)試驗變異函數(shù)曲線對區(qū)域化變量進(jìn)行分析合適的理論模型實際中區(qū)域化變量的變化性很復(fù)雜：(1)可能在不同方向上有不同的變異性；(2)在同一方向上包含不同尺度上的多層次的變異性這么復(fù)雜！？第四節(jié)變異函數(shù)的結(jié)構(gòu)分析一、結(jié)構(gòu)分析、套合結(jié)構(gòu)概念采樣數(shù)51礦床或礦體的變異性往往由多種原因引起采樣、樣品制備及分析等過程所產(chǎn)生的誤差原因礦物成分的變化，如金礦等品位變化劇烈的礦床上尤為明顯礦層與夾層的交替變化礦床分布引起的變異≈01～ncm米至百米公里尺度不同原因引起的變異特性，其變異尺度的大小不同顯然，大尺度的變異總是包含著小尺度的變異，小尺度的變異在大尺度變異曲線上只能作為“塊金效應(yīng)”出現(xiàn)礦床或礦體的變異性往往由多種原因引起采樣、樣品制備及分析等過52土壤的空間變異性與土壤母質(zhì)、氣候、水文、地形和生物等因素相關(guān)Z(x)Z(x+h）hh≈0h≤1mh≤100m取樣和測定誤差+其它因素，如水分+地形影響合適的理論模型??？結(jié)構(gòu)分析土壤的空間變異性與土壤母質(zhì)、氣候、水文、地形和生物等因素相關(guān)53結(jié)構(gòu)分析:就是構(gòu)造一個變異函數(shù)模型對于全部有效結(jié)構(gòu)信息作定量化的概括，以表征區(qū)域化變量的主要特征。結(jié)構(gòu)分析的主要方法:套合結(jié)構(gòu)套合結(jié)構(gòu)(nestedstructure):就是把分別出現(xiàn)在不同距離h上和(或）不同方向α上同時起作用的變異性組合起來套合結(jié)構(gòu)表達(dá)式:套和結(jié)構(gòu)可以表示為多個變異函數(shù)之和，每一個變異函數(shù)代表一種特定尺度上的變異性，表達(dá)式為：γi(h)可以是相同或不同的理論模型結(jié)構(gòu)分析:就是構(gòu)造一個變異函數(shù)模型對于全部有效結(jié)構(gòu)信息作定量54二、單一方向上的套合

設(shè)區(qū)域化變量Z(x)在某一方向上的變異性是由γ0(h)、γ1(h)和γ2(h)組成微觀尺度，純塊金效應(yīng)模型變程為a1的球狀模型變程為a2的球狀模型(a2>a1)二、單一方向上的套合設(shè)區(qū)域化變量Z(x)在某一方向上的55變異函數(shù)及結(jié)構(gòu)分析課件56變異函數(shù)及結(jié)構(gòu)分析課件57作業(yè)：

設(shè)區(qū)域化變量Z(x)在某一方向上的變異性是由γ0(h)、γ1(h)和γ2(h)組成，請求出其套合結(jié)構(gòu)表達(dá)式并寫出計算過程微觀尺度，純塊金效應(yīng)模型變程為a的球狀模型變程為3a的指數(shù)模型作業(yè)：設(shè)區(qū)域化變量Z(x)在某一方向上的變異性是由γ058變異函數(shù)及結(jié)構(gòu)分析課件59——地統(tǒng)計學(xué)的工具第四章變異函數(shù)及結(jié)構(gòu)分析——地統(tǒng)計學(xué)的工具第四章變異函數(shù)及結(jié)構(gòu)分析60一、協(xié)方差函數(shù)的計算公式第一節(jié)協(xié)方差函數(shù)和變異函數(shù)的性質(zhì)設(shè)區(qū)域化變量Z(x)滿足(準(zhǔn))二階平穩(wěn)假設(shè)，h為兩樣本點空間分隔距離，Z(xi)與Z(xi+h)分別是Z(x)在空間位置xi和xi+h上的觀測值(i=1,2,…N(h)),則計算協(xié)方差的公式為：一、協(xié)方差函數(shù)的計算公式第一節(jié)協(xié)方差函數(shù)和變異函數(shù)的性質(zhì)61變異函數(shù)及結(jié)構(gòu)分析課件62協(xié)方差函數(shù)曲線圖：以h為橫坐標(biāo)，C#(h)為縱坐標(biāo)作圖協(xié)方差函數(shù)曲線圖：以h為橫坐標(biāo)，C#(h)為縱坐標(biāo)作圖63二、協(xié)方差函數(shù)的性質(zhì)區(qū)域化變量Z(x)在二階平穩(wěn)假設(shè)下，其協(xié)方差函數(shù)存在且平穩(wěn)，定義為1.C(0)=Var[Z(x)]≥0，即先驗方差不能小于零2.C(h)=C(-h)，即C(h)對h=0的直線對稱，是一個偶函數(shù)二、協(xié)方差函數(shù)的性質(zhì)區(qū)域化變量Z(x)在二階平穩(wěn)假設(shè)下，其協(xié)64證：令x-h=y，則x=y+h，帶入上式得圖形特征及含義證：令x-h=y，則x=y+h，帶入上式得圖形特征及含義653.|C(h)|≤C(0)，即協(xié)方差函數(shù)絕對值小于等于先驗方差證：3.|C(h)|≤C(0)，即協(xié)方差函數(shù)絕對值小于等于先664.|h|→∞時,C(h)→0，或?qū)懽鰿(∞)=0，即當(dāng)空間距離很大時，協(xié)方差函數(shù)值很小意義(空間局限性)：當(dāng)距離很大時，Z(x)和Z(x+h)之間的線性相關(guān)基本不存在5.C(h)必須是一個非負(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ū)域化變量的任意有限線性組合，即：4.|h|→∞時,C(h)→0，或?qū)懽鰿(∞)=0，即當(dāng)67較難理解則由C(xi-xj）(i,j=1,2…n)構(gòu)成的協(xié)方差函數(shù)矩陣是非負(fù)定矩陣，即C(h)為非負(fù)定函數(shù)二階平穩(wěn)區(qū)域化變量的協(xié)方差函數(shù)是有條件的較難理解則由C(xi-xj）(i,j=1,2…n)構(gòu)成的協(xié)方68三、實驗(經(jīng)驗)變異函數(shù)(experimentalvariogram)的計算公式設(shè)區(qū)域化變量Z(x)滿足(準(zhǔn))二階平穩(wěn)條件或(準(zhǔn))本征假設(shè)，h為兩樣本點空間分隔距離，Z(xi)與Z(xi+h)分別是Z(x)在空間位置xi和xi+h上的觀測值(i=1,2,…N(h)),則計算實驗變異函數(shù)的公式為：變異函數(shù)曲線圖：以h為橫坐標(biāo)，γ#(h)為縱坐標(biāo)作圖三、實驗(經(jīng)驗)變異函數(shù)(experimentalvari69變異函數(shù)計算實例（1）一維變異函數(shù)的計算x1x2x3x4x5x6x7x8x9x104

3

4

5797877以下為一研究對象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點間分隔距離h=1米，計算γ#(h)變異函數(shù)計算實例（1）一維變異函數(shù)的計算x1x270變異函數(shù)及結(jié)構(gòu)分析課件71變異函數(shù)及結(jié)構(gòu)分析課件72兩方面理解：變異性的理解與相關(guān)性的理解兩方面理解：變異性的理解與相關(guān)性的理解73作業(yè)：x1x2x3x4x5x6x7x8

2

4

3

15364以下為一研究對象在水平方向上的采樣數(shù)據(jù)，滿足二階平穩(wěn)或本征假設(shè)，采樣值如圖所示，點間分隔距離h=1米，計算γ#(1),γ#(2),γ#(3)作業(yè)：x1x2x3x4x574（2）二維變異函數(shù)的計算下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，*號處為無數(shù)據(jù)點，點間距離h為100米，請分別計算南北、東西、西北和東南方向上的變異函數(shù)值。（2）二維變異函數(shù)的計算下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，*號處75變異函數(shù)及結(jié)構(gòu)分析課件76變異函數(shù)及結(jié)構(gòu)分析課件77變異函數(shù)及結(jié)構(gòu)分析課件78西北和東南方向上的變異函數(shù)值的計算，注意分隔距離h的確定和樣本數(shù)據(jù)對的查找西北和東南方向上的變異函數(shù)值的計算，注意分隔距離h的確定和樣79變異函數(shù)及結(jié)構(gòu)分析課件80作業(yè)：下圖為正方形網(wǎng)格狀的采樣數(shù)據(jù)，網(wǎng)格交叉空白處為無數(shù)據(jù)點，點間距離h為a米，請分別計算南北方向γ#(a),

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

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

1、拋物線型（parabolictype）當(dāng)|h|→0時,γ(h)→A|h|2(A為常數(shù))，即變異函數(shù)曲線在原點處趨向一條拋物線，反映區(qū)域化變量是具有高度連續(xù)性的，如礦層厚度四、變異函數(shù)在原點處的性狀可反映區(qū)域化變量的空間連續(xù)性變93

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

3、間斷型（discontinuoustype）當(dāng)|h|→0時,γ(h)→C0，即變異函數(shù)曲線在原點處間斷，說明塊金效應(yīng)存在，又稱“塊金效應(yīng)型”，反映區(qū)域化變量的連續(xù)性很差，但當(dāng)h增大時，γ(h)又變的較為連續(xù)了，如金品位2、線性型（lineartype）當(dāng)|h|→0時,γ94

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

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

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

γ(h)，則常見的變異函數(shù)模型如下：變異函數(shù)的理論模型有基臺值模型無基臺值模型可以有或無基臺值模型：孔穴效應(yīng)模型球狀模型、指數(shù)模型高斯模型線性有基臺模型純塊金效應(yīng)模型冪函數(shù)模型對數(shù)模型線性無基臺模型第三節(jié)變異函數(shù)的理論模型思考：是否有了采樣數(shù)據(jù)及變異函97一、有基臺值模型

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

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，則變程為3a2、指數(shù)模型（exponentialmodel）式中：100

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ù)模型，其圖形為：3、高斯模型（gaussianmodel）式中：C0，101模型通過原點切線與基臺值線交點的橫坐標(biāo)變程原點處的性狀球狀2a/3a直線指數(shù)a3a直線高斯無交點√3

a拋物線三種模型的比較模型通過原點切線與基臺值線交點的橫坐標(biāo)變程原點處的性狀球狀2102

4、線性有基臺值模型（linearwithsillmodel）式中：C0，C意義同前，A為常數(shù)，表示直線的斜率，變程為a4、線性有基臺值模型（linearwithsillm103

5、純塊金效應(yīng)模型（purenuggeteffectmodel）此時，C0=C(0)此種模型意味著區(qū)域化變量為隨機(jī)分布，樣點間的協(xié)方差函數(shù)對于所有距離h均等于0，即變量不存在空間相關(guān)性5、純塊金效應(yīng)模型（purenuggeteffect104變異函數(shù)及結(jié)構(gòu)分析課件105二、無基臺值模型

1、冪函數(shù)模型（powermodel）若與模型相應(yīng)的區(qū)域化變量不滿