油藏描述储层参数空间分布预测.ppt

3.趋势面模型参数的计算 4.趋势面分析应用实例 (3) 变异函数的性质 设Z(x)是区域化变量，在满足二阶平稳假设条件下，变异函数式具有如下性质： ① r(0) =0，即在h=0处，变异函数为0； ② r(h) = r(-h) ，即 r(h) 关于直线h=0是对称的，它是一个偶函数； ③ r(h)≥0，即 r(h)只能大于或等于0； ④ |h|→∞时，r(h)→C(0)，或r(∞)=C(0)，即当空间距离增大时，变异函数接近先验方差 5.变异函数的参数 变异函数有四个非常重要的参数，即基台值、变程、块金值和分维数。 前3个参数可以直接从变异函数图中得到。它们决定变异函数的形状与结构。 变异函数的形状反映自然现象空间分布结构或空间相关的类型，同时还能给出这种空间相关的范围。 当变异函数随着间隔距离h的增大，从非零值达到一个相对稳定的常数时，该常数称为基台值C0+C， 当间隔距离h=0时，r(0)= C0，该值称为块金值。 基台值是系统或系统属性中最大的变异，变异函数达到基台值时的间隔距离a称为变程。变程表示在h≥a以后，区域化变量Z(x)空间相关性消失。 块金值表示区域化变量在小于抽样尺度时非连续变异，由区域化变量的属性或测量误差决定。 上述三个参数可从变异函数曲线图直接得到，或通过估计曲线回归参数得到。 第4个参数，即分维数用于表示变异函数的特性，由变异函数r(h)和间隔距离h之间的关系确定： 2 r(h) = h(4-2D) 分维数D为双对数直线回归方程中的斜率，它是一个无量纲数。 分维数D的大小，表示变异函数曲线的曲率，可以作为随机变异的量度。 ③ 高斯模型 (2) 无基台值的模型 ① 幂函数模型 (5) 变异函数图（变差图） 　　r(x，h)一般是依赖于x和h两个自变量的。如果又知道r(x，h) 与X取何值无关，只依赖于h(基本步长),则可把变异函数r(x，h)写为r(h)。此时，以h为横坐标，以r(h)值为纵坐标作出的图形就叫变异函数图（变差图） 。 设Z(x)为一区域化随机变量，并满足二阶平稳假设，h为两样本点空间分隔距离，Z(xi)、Z(xi+h)分别是区域化变量Z(x)在空间位置xi和xi+h处的实测值[i=1,2,…,N(h)]，那么，变异函数r*(h)的离散计算公式为: 3.变异函数的计算 例1：Z(x)是一维区域化变量，Z(1)=2, Z(2)=4, Z(3)=3, Z(4)=1, Z(5)=5, Z(6)=3, Z(7)=6, Z(8)=4。计算r*(1), r*(2), r*(3) 的值并绘制变差图。 解：由公式 得 ={[Z(1)-Z(2)]2+ [Z(2)-Z(3)]2 +…+ [Z(7)-Z(8)]2}/(2×7) =(22 +12 +22 +42 +22 +32 +22)/14 =42/14=3.00 同理 r*(2) = (12 +32 +22 +22 +12 +12)/2×6 = 20/12 = 1.67 r*(3) = (12 +12 +02 +52 +12)/2×5 = 28/10 = 2.80 以h为横坐标，r*(h)为纵坐标绘制变差图。 h r*(h) 例2.设Z(x)为二维区域化变量。下图所示的正方形网格节点处有一批已知数据，小正方形边长为a。试计算右图所示的a1,a2,a3,a4四个方向的实验变异函数值（只算前三个h值），并作出a1,a3两个方向的变差图。 解：由公式 计算各个方向的变异函数值。 (1)a1方向：r*(a)=197/(2*24)=4.10 r*(2a)=336/(2*20)=8.40 r*(3a)=435/(2*18)=12.08 (2)a3方向：r*(a)=187/(2*22)=4.25 r*(2a)=296/(2*18)=8.22 r*(3a)=327/(2*15)=10.90 (3)a2方向：r*(√2 a)=233/(2*18)=6.47 r*(2√2 a)=315/(2*14)=11.25 r*(3√2 a)=247/(2*8)=15.44 (4)a4方向：r*(√2 a)=191/(2*19)=6.47 r*(2√2 a)=381/(2*16)=11.91 r*(3√2 a)=345/(2*10)=17.25 0 a √2a 2a 2√2a 3a 3√2a h 5 10 15 r*(h) a3方向变差图 a2方向变差图 a1方向变差图 a4