petrel建模算法优选与适用，附建模入门教程.doc

3建模途径 对于建模方法目前各种书籍文献论述的比较清楚了，这里只对目前通过长期建模验证的比较适用的建模方法进行一个简单概括。 3.1确定性建模 3.1.1确定性建模概述 确定性建模是对井间未知区给出确定性的预测结果，即试图从具有确定性资料的控制点(如井点)出发，推测出井间确定的、唯一的储层地质特征。 3.1.2确定性建模方法 井间储层插值，即对井间储层参数进行内插和井外推测，是确定性建模的重要方法，包括数理统计插值方法和地质统计学克里金插值方法。 （一）数理统计插值 常规数理统计方法很多，如三角网法、距离反比加权法、径向基函数法、多重网格逼近法、离散光滑插值法、样条插值方法、最近邻点法、移动最小二乘法等。一般插值方法可分为局部插值与整体插值两大类。局部插值法的特点是每个插值点只影响其周围的局部区域，如距离反比、B样条插值等；整体插值法是基于整体插值点，一般要求解一个线性方程组，变动或改变一个插值点，就会改变整个插值曲面，如薄板样条插值等。 在储层建模中，应用该方法的前提是地质参数在井间具有数理统计关系，即某种数学函数关系，如三角网方法的前提是井间参数值是井孔参数与井间距离的函数。应用这一函数关系，即可对多井井间（2维剖面、2维平面及3维空间）进行储层参数插值，并建立储层地质模型。当然，这类方法也可整合地震信息进行插值。 不同插值方法的内涵及应用范畴各有所差别（表3-1）。在实际应用中，应根据地质参数的空间分布特征、原始数据类以及插值网格节点规模，选择合适的插值方法。 表3-1 主要数理统计插值方法简表 算法算法描述算法应用范围三角网线形插值Triangulation with linear interpolation，算法首先基于已知点连接三角形网络，并以此为基础，将处于各个三角形内的未知点，通过各个三角形边作线性插值。算法稳定性好，效率高；插值结果数值介于已知点数值范围；在已知点少的情况下会出现明显的受三角形控制趋势；同时三角网络外的未知点将无法进行插值计算。算法在已知点数目较多且分布均匀的情况下插值效果较好，适合于构造或是平面参数分布插值；在已知点稀疏且分布不均的情况下，插值效果表现出明显的受已知点三角形控制的趋势；距离反比 加权插值Inverse distance to a power，算法将未知点值表示为其周围已知点的加权平均，权系数与到各已知点的距离成反比。算法简单、效率高；相比三角网法，可对每个未知点进行插值计算；插值结果数值介于已知点数值范围，但容易在井点处出现“牛眼”现象。距离反比法是应用最广的插值方法，其可单独使用，也常与其他插值方法进行综合。基本上适合于任意类型的网格面或三维网格体属性插值。径向基函数插值Radio basic function，算法在径向基函数空间中，以各已知点为中心，距离r为变量的基函数的线性组合构建插值函数。采用不同特性的径向基函数，可适应不同建模模型的需要。该方法在求解基函数权系数时，需要解一个同已知点规模的线性方程组，因此已知点个数最好限定在一定规模。径向基函数法可适用于各种类型的网格化面插值，如构造层面、平面储层参数分布等；已知点个数不能太多，最好在200个左右，且分布均匀；只适用于井点数据，不能综合地震、等值线及断层信息。多重网格逼近法Convergent grid，算法的基本思想，是将整个插值过程分解为网格大小粗到细的多次迭代计算。算法效率非常高、稳定性好；对没有原始点分布位置处，也能保持很好的趋势。可综合地震、等值线及井点数据插值。例如，地震层位解释测线数据的网格化插值；按等值线或密集散点方式给出的构造面或平面储层属性参数的网格化插值。离散光滑插值DSI—Discrete Smooth Interpolation，算法基于对目标体的离散化，通过设立目标准则及约束条件，将各种地质模型特征加入到算法中，并最终通过迭代方法求模型的最优解。DSI算法非常灵活，可根据不同情况，设立相适应的目标准则与约束条件。由于算法基于整体的迭代求解，当插值网格节点较大时，算法效率会变得很低。目前DSI方法已成为地质体几何建模的主流技术，特别是对于复杂地质体几何外形建模，如岩浆侵入岩体、地层复杂断裂以及挠曲褶皱，等等。同时，可将地质目标体边界作为约束条件，可实现平面相控储层参数插值等。样条插值法Spline，算法使用某种数学函数，对一些限定的点值，通过控制估计方差，利用一些特征节点，用多项式拟合的方法来产生平滑的插值曲线。样条插值法包括多种类型：算法适用于构造层面插值。最邻近点法Closest point，算法对每一个未知点，从所有已知点中找到与之最近的一个，然后将此已知点的值赋给这个未知点。计算结果是一个阶梯函数。算法适用于大规模已知点