基于局部各向异性的非平稳多点地质统计学算法.PDFVIP

第41卷第2 期 物 探 与 化 探 Vol .41，No .2 2017 年 4 月 GEOPHYSICAL GEOCHEMICAL EXPLORATION Apr .，2017 doi ： 10.11720/wtyht .2017.2.11 喻思羽，李少华,段太忠，等.基于局部各向异性的非平稳多点地质统计学算法[J ].物探与化探,2017,41(2) :262-269.http ://d〇i .〇rg/10.11720/ wtyht .2017.2.11 Yu S Y,Li S H,Duan T Z,et al .Non-stationaiy multiple-point geostatistics algorithm base on local anisotropy [J ] .Geophysical and Geochemical Explora - tion ,2017,41(2) :262-269.http ://doi .org/10.11720/wtyht .2017.2.11 基于局部各向异性的非平稳多点 地质统计学算法 喻 思 羽 1 , 李 少 华 1 ,段 太 忠 2 , 王 鸣 川 2 (1.长江大学地球科学学院，武 汉 430100; 2.中石化石油勘探开发研究院，北 京 100083) 摘 要 ：基于平稳性假设的传统多点地质统计学在非平稳建模中存在不足。提出一种基于局部各向异性分区的多 点地质统计建模方法，利用不同方向变差函数变程组合的玫瑰花图定量表征训练图像的局部各向异性，结合经典 多维尺度分析和K 均值聚类分析对非平稳训练图像进行自动分区，各子区域相互独立，分区内运用传统平稳性多 点建模算法模拟，最终实现非平稳性建模。以裂缝网络建模为例，新方法较好地模拟并再现训练图像的非平稳性， 对多点地质统计建模算法的非平稳建模具有很好借鉴意义。 关键词 ：各向异性;玫瑰花图；多维尺度分析;多点地质统计学;非平稳性 中图分类号 ：P613.4 文献标识码： A 文章编号 ：1000-8918(2017)02-0262-08 作为空间独立不相关的多点统计信息库,Zhang[7]和 0 引言 Arpat [8]把这种思路分别运用于SNESIM和 SIMPAT 算法上，取得一定效果。但建立若干旋转、缩放变换 多点地质统计建模方法（multiple-point geosta- 的训练图像的多点统计信息库，不仅增加内存开销 tistics ,MPS )是储层地质建模研究领域的前沿[1，2]。 同时效率较低，此外旋转、缩放参数的选取较为困 基于空间多点相关性,多点地质统计建模方法预先 难。第 2 种方法是基于分区方法，2011年 , Mehrdad 用数据样板扫描提取训练图像的数据样式，即重复 Honarkhah通过形态学特征，例 如 Gabor纹理过滤 观察训练图像所有位置的样本点，得到训练图像内 器[9]，把非平稳的训练图像分割为多个子区域，然 在的几何形态及非均质特征，完成未知区域的预 后提取各个子区域内部的空间多点相关性，再用于 测[3]。平稳性假设是经典地质统计学的理论基 相应子区域的模拟。由于纹理过滤器Gabor是二维 础[4]，即区域化变量均值、方差必须平稳且与位置 图像分析工具，其难以用于三维建模。第 3 种方法 无关。因此训练图像必须具备统计平稳性[5]。然 是基于距离函数[10]建立储层模式与位置的关联性， 而自然界没有绝对