基于虚拟井随机模拟的地震岩相预测方法及系统技术方案

公开了一种基于虚拟井随机模拟的地震岩相预测方法及系统。该方法可以包括：对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；根据岩石物理模型和褶积模型正演得到每个虚拟井对应的模拟地震道数据；计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取最优相似系数最大的N个虚拟井作为优选虚拟井，以优选虚拟井中出现次数最多的岩相类型作为实际地震道的岩相类型。本发明专利技术通过模拟大量虚拟井，建立测井与地震的关系，匹配虚拟井模拟记录与实际地震记录，实现岩相预测，可避免地震反演带来的误差，适用于钻井较少的区域，减少人工干预。

Seismic facies prediction method and system based on virtual well stochastic simulation

【技术实现步骤摘要】
基于虚拟井随机模拟的地震岩相预测方法及系统
本专利技术涉及地震勘探
，更具体地，涉及一种基于虚拟井随机模拟的地震岩相预测方法及系统。
技术介绍
储层岩相分布信息对于油藏描述具有重要的参考意义。传统方法首先通过地震反演获取储层弹性参数信息，然后在岩石物理建模基础上将其进一步转换为岩相信息。这种方法忽视了地震反演结果的不确定性、井震数据尺度差异等因素对于储层地震岩相识别的影响，难以评价岩相预测结果的不确定性，为后续油藏描述、储量计算等带来一定的风险。因此，有必要开发一种基于虚拟井随机模拟的地震岩相预测方法及系统。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种基于虚拟井随机模拟的地震岩相预测方法及系统，其能够通过模拟大量虚拟井，建立测井与地震的关系，匹配虚拟井模拟记录与实际地震记录，实现岩相预测，可避免地震反演带来的误差，适用于钻井较少的区域，减少人工干预。根据本专利技术的一方面，提出了一种基于虚拟井随机模拟的地震岩相预测方法。所述方法可以包括：对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；根据岩石物理模型和褶积模型正演得到所述每个虚拟井对应的模拟地震道数据；计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取所述最优相似系数最大的N个虚拟井作为优选虚拟井，以所述优选虚拟井中出现次数最多的岩相类型作为所述实际地震道的岩相类型。优选地，所述模拟地震道数据为：dsyn(t,θ)＝w(t,θ)*Rpp(t,θ)(1)其中，dsyn(t,θ)为入射角度为θ的地震记录，w(t,θ)为地震子波，Rpp(t,θ)为反射系数。优选地，根据公式(2)计算所述最优相似系数：BSij＝BSFij·BSVij(2)其中，BSij为最优相似系数，BSFij为形相似系数，BSVij为值相似系数。优选地，所述形相似系数为：其中，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。优选地，所述值相似系数为：其中，Dij为计算参数，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。根据本专利技术的另一方面，提出了一种基于虚拟井随机模拟的地震岩相预测系统，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现以下步骤：对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；根据岩石物理模型和褶积模型正演得到所述每个虚拟井对应的模拟地震道数据；计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取所述最优相似系数最大的N个虚拟井作为优选虚拟井，以所述优选虚拟井中出现次数最多的岩相类型作为所述实际地震道的岩相类型。优选地，所述模拟地震道数据为：dsyn(t,θ)＝w(t,θ)*Rpp(t,θ)(1)其中，dsyn(t,θ)为入射角度为θ的地震记录，w(t,θ)为地震子波，Rpp(t,θ)为反射系数。优选地，根据公式(2)计算所述最优相似系数：BSij＝BSFij·BSVij(2)其中，BSij为最优相似系数，BSFij为形相似系数，BSVij为值相似系数。优选地，所述形相似系数为：其中，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。优选地，所述值相似系数为：其中，Dij为计算参数，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。本专利技术具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的基于虚拟井随机模拟的地震岩相预测方法的步骤的流程图。图2a、图2b和图2c分别示出了根据本专利技术的一个实施例的理论岩相分布及对应的孔隙度和泥质含量的示意图。图3示出了根据本专利技术的一个实施例的实际地震道的示意图。图4示出了根据本专利技术的一个实施例的预测的实际地震道的岩相分布的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的基于虚拟井随机模拟的地震岩相预测方法的步骤的流程图。在该实施例中，根据本专利技术的基于虚拟井随机模拟的地震岩相预测方法可以包括：步骤101，对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；步骤102，根据岩石物理模型和褶积模型正演得到每个虚拟井对应的模拟地震道数据；步骤103，计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取最优相似系数最大的N个虚拟井作为优选虚拟井，以优选虚拟井中出现次数最多的岩相类型作为实际地震道的岩相类型。在一个示例中，模拟地震道数据为：dsyn(t,θ)＝w(t,θ)*Rpp(t,θ)(1)其中，dsyn(t,θ)为入射角度为θ的地震记录，w(t,θ)为地震子波，Rpp(t,θ)为反射系数。在一个示例中，根据公式(2)计算最优相似系数：BSij＝BSFij·BSVij(2)其中，BSij为最优相似系数，BSFij为形相似系数，BSVij为值相似系数。在一个示例中，形相似系数为：其中，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。在一个示例中，值相似系数为：其中，Dij为计算参数，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。具体地，测井岩相划分是以岩相作为储层参数的载体，在一定的本文档来自技高网...

【技术保护点】
1.一种基于虚拟井随机模拟的地震岩相预测方法，其特征在于，包括：/n对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；/n根据岩石物理模型和褶积模型正演得到所述每个虚拟井对应的模拟地震道数据；/n计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取所述最优相似系数最大的N个虚拟井作为优选虚拟井，以所述优选虚拟井中出现次数最多的岩相类型作为所述实际地震道的岩相类型。/n

【技术特征摘要】
1.一种基于虚拟井随机模拟的地震岩相预测方法，其特征在于，包括：
对于一道实际地震道数据，基于马尔科夫先验约束生成多个虚拟井，每个虚拟井包括岩相和储层参数；
根据岩石物理模型和褶积模型正演得到所述每个虚拟井对应的模拟地震道数据；
计算每个虚拟井对应的模拟地震道数据与实际地震道数据的最优相似系数，选取所述最优相似系数最大的N个虚拟井作为优选虚拟井，以所述优选虚拟井中出现次数最多的岩相类型作为所述实际地震道的岩相类型。


2.根据权利要求1所述的基于虚拟井随机模拟的地震岩相预测方法，其中，所述模拟地震道数据为：
dsyn(t,θ)＝w(t,θ)*Rpp(t,θ)(1)
其中，dsyn(t,θ)为入射角度为θ的地震记录，w(t,θ)为地震子波，Rpp(t,θ)为反射系数。


3.根据权利要求1所述的基于虚拟井随机模拟的地震岩相预测方法，其中，根据公式(2)计算所述最优相似系数：
BSij＝BSFij·BSVij(2)
其中，BSij为最优相似系数，BSFij为形相似系数，BSVij为值相似系数。


4.根据权利要求3所述的基于虚拟井随机模拟的地震岩相预测方法，其中，所述形相似系数为：



其中，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。


5.根据权利要求3所述的基于虚拟井随机模拟的地震岩相预测方法，其中，所述值相似系数为：



其中，Dij为计算参数，xi和xj分别代表一道地震数据构成的向量，xi(k)和xj(k)分别表示向量中的第k个元素，m为每个向量的元素个数，和分别表示两个向量元素的平均值。


6.一种基于虚拟井随机模拟的地震岩相预测系统，其特...

【专利技术属性】
技术研发人员：李京南，刘婵娟，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11