[发明专利]一种多点地质统计学沉积相建模方法

摘要

本发明公开了一种多点地质统计学沉积相建模方法。该方法在原有的多点地质统计学沉积相建模方法的基础上，提出了几何因子的概念，采用几何因子限定空间求取反映沉积相带规模的条件概率，应用校正比例恒等式的方法对扫描训练图像得到的条件概率分布曲线进行修正，对原有的多点地质统计学方法进行了改进。应用改进后的算法进行建模实现表明，引入几何因子的多点地质统计学方法继承了传统多点地质统计学再现沉积相几何形态的优点，还弥足了其在处理相带规模和连续性等方面的不足。

主权项

一种多点地质统计学沉积相建模方法，包括以下步骤：S100，以地质概念模型为指导对研究区进行模式拟合，建立原型地质模型；S200，根据原型地质模型建立符合沉积模式的训练图像并创建搜索树，对原型地质模型进行变差函数分析得到各沉积相主方向的变程，以获取几何因子；S300，判断区域空间中是否存在条件数据：若存在条件数据，则建立网格，将条件数据分配到最近的网格节点，引入几何因子，在每个访问节点的周围构建几何因子空间，通过指示克里金方法确定几何因子空间内待估点处的沉积相规模概率作为几何因子概率P(A|C)；S400，指定能够访问所有网格节点的随机访问路径；S500，对于每个访问节点，判断是否为条件数据：若否，执行步骤S510；若是，确定该节点的沉积相类型，执行步骤S520；S510，判断该节点处条件数据事件在搜索树中是否有记录：若有记录，执行步骤S520；若无记录，执行步骤S540；S520，根据以该节点为中心的条件数据事件在训练图像中出现的重复次数得到沉积相估计概率作为该点的局部条件概率P(A|B)；S530，通过校正比例恒等式将局部条件概率P(A|B)与几何因子概率P(A|C)相结合，得到条件累积概率P(A|B,C)，由此获得条件概率分布函数，执行步骤S550；S540，令该节点的条件概率为几何因子概率P(A|C)，由此获得条件概率分布函数；S550，在条件概率分布函数中提取一个模拟值，将其当作实际的观测数据加入到原来的条件数据集中，作为随后其他节点的模拟条件；S600，在访问完所有节点并确定了所有节点的模拟结果的基础上，根据沉积相带的空间结构和几何形态，得到所需的沉积相模型。