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

四、克里格插值法计算步骤 (１)准备工作 首先了解地质情况，确定资料点坐标和变异函数方向，选择步长，进行不规则井网处理。 (２)计算实验变异函数 式中：Z(Xi),Z(Xi+h)――相距为h的N(h)对点Xi和 Xi+h(i=1,2,…,N(h))处的储层参数观测值； r*(h)――相隔为h的变异函数值。 (３)实验变异函数的最优拟合 选用变异函数的理论模型，对实验变异函数进行拟合，例如选用球状模型公式为 式中：C0――块金常数，反映空间变化连续性； a――变程，反映空间点相关性范围； C1,C0+C1――拱高和基台值，反映参数变化强度。 利用加权多项式回归法拟合球状模型变差图，分别指出x方向和y方向理论模型的变异函数(C0、C1和a)。 (４)不同方向的结构套合 将x和y方向理论模型进行套合，形成统一的二维同性结构。 (５)选用克立格法预测参数 利用高斯消元法求解普通克立格和泛立格方程组，得到各权系数λi和μ，代入公式： 可算出待估点X0处的估值Z(X0)； 用相应的估计方差计算公式可算出相应的估计方差。 复杂的区域化变量往往包含各种尺度上的多层次、多方向的变化性，反映在变差函数上即为多层次结构。将不同结构组合为统一结构的过程称为“结构套合” 结构套合 各层次套合 大尺度的变化性总是包含着小尺度的变化性，但却不能从大尺度的变化性中区分出小尺度的变化性。 可以用反映各种不同尺度变化性的多个变差函数之和来表示一个套合结构。 代表微观变化性的变程极小的球状模型，可近似地看作纯块金效应型 球状模型，没有块金常数，基台值为C1，变程为a1 ，反映了小规模范围的变化 球状模型，没有块金常数，基台值为C2，变程较大，为a2 ，反映了大规模范围的变化 其中 则套合结构的表达式为 ①对于几何各向异性，先根据异向比压缩距离轴，使之成为各向同性的模型； ②对于带状各向异性，运用模型叠加的方法加以处理。先用压缩距离轴的办法，使其变程变为相同，然后再把具有相同变程的两个球状模型叠加起来，构成一个新的球状模型 各方向套合 （将各向异性套合为各向同性，以便于在克里金估计时，不同方向均可用统一的结构模型计算实际的变差函数值） 例3：如图所示，油藏Ω中有一个方块，其边长为3，一个样品ω1 位于方块的中心，某油藏参数x的值为x1=1.2；另一个样品ω2 位于方块的一角，其参数x值为x2=0.5。变异函数的方程为： 假设该参数在油藏中满足平稳假设且无漂移，试用克里金法估计方块V的参数x的值。 方块与样品相对位置及变异函数示意图 解: (1)计算 r(ωi, ωj) ω1与ω2之间的距离为 ω1和ω2到自身的距离为 0根据公式(7)，有： (2)计算r(ω1,V)和r(ω2,V) 将V等分为如图所示的9个小块，根据公式(5)，有： 式中z0为ω1的位置，zi为第i个小块的中心位置。 当i=2时，样品ω1距第二小块的距离为1，因此r(z0,z2)=r(1)=0.5。 类似地可以求出任意i的r(z0,zi)： 仿照上述做法，求得： r(ω2,V)=0.882 (3)建立克里金方程组并求解： 将求得的r(ωi, ωj)和r(ωi,V)代入克里格方程组，有： 解上面的方程组得： λ1=0.673, λ2=0.327, μ=0.209 ? (4)求方块V的参数x的值 方块V的参数x的估值为： (5)计算方块参数x的克里金方差： 计算克里金方差 需要计算r(V,V) ，根据式(2) 式中zi和zj为两个小方块的位置，对于每一对zi和zj，r(zi , zj )的算法 与第(2)步中r(z0 , zi )的算法相同。 计算结果如下： r(V,V) =0.683 将有关数值代入普通克立格估计方差的计算公式得： =0.175 普通克立格法使用得比较普遍，它对孔隙度、泥质含量这些变化比较平稳的参数估值效果比较明显，因为实施普通克立格法的前提条件即是E[Z(X)]等于一个常数。也就是说，要求参数变化幅度不剧烈。但它对渗透率、含油饱和度这样一些参数估值效果不佳。 五、克立格法的适用条件 泛克立格法是一种考虑了漂移影响的克立格估值方法，即E[Z(X)]等于一个多项式(漂移)，而不是一个常数，从理论上讲，它对那些受漂移控制比较明显的参数估值效果比较好。 1.数理统计法估值没有考虑各信息点的相互关系，而只考虑信息点与待估点之间的几何特征(如距离)。而克立格法既考虑了各信息点间的相关性，以及各信息点与待估点间的相关性，同时还考虑各数据点空间分布的结构特征。 六、地质统计学