碳酸盐岩古暗河充填性预测方法技术

本发明专利技术公开了一种碳酸盐岩古暗河充填性预测方法，该方法包括：分析测井资料，建立溶洞的充填参数的先验分布，以及建立表征所述充填参数与岩石弹性参数的先验关系的岩石物理模型；根据所述岩石物理模型计算溶洞的岩石弹性参数；以岩石弹性参数作为已知参数，充填参数作为反演目标，并结合充填参数的先验分布，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数。由此可见，利用本发明专利技术可以定量获得古暗河充填物性质参数，例如泥质含量、砂质含量、灰质含量及孔隙度，实现暗河内充填物性质及充填程度的定量预测。

Prediction method of carbonate rock filling

The present invention discloses a method for predicting the fillability of carbonate ancient underflow, which includes: analyzing logging data, establishing a prior distribution of filling parameters of karst caves, and establishing a petrophysical model to characterize the prior relationship between filling parameters and rock elastic parameters, and calculating the rock of karst caves according to the said petrophysical model. The rock elastic parameters are taken as known parameters, the filling parameters are taken as inversion targets, and combined with the prior distribution of filling parameters, the filling parameters of karst cave based on Bayesian are inverted to obtain the predicted filling parameters of karst cave. Thus it can be seen that the method can be used to quantitatively obtain the property parameters of the fillings in the ancient dark river, such as the shale content, the sand content, the ash content and the porosity, so as to realize the quantitative prediction of the filling property and the filling degree in the dark river.

【技术实现步骤摘要】
碳酸盐岩古暗河充填性预测方法
本专利技术涉及碳酸盐岩缝洞型储集体预测领域，具体涉及一种碳酸盐岩古暗河充填性预测方法。
技术介绍
塔河碳酸盐岩缝洞型油藏储集空间以溶洞为主，溶洞储集性能主要受溶洞规模、溶洞内充填物类型及充填程度的影响，不同充填特征的溶洞型储集体产量差异大。古岩溶期形成的岩溶暗河体系是大型溶洞最为发育的区域，但暗河在形成过程中及埋藏后经历了广泛的充填、垮塌破坏，致使原洞穴空间部分或全部被砂泥质或角砾岩等充填，目前常用的溶洞充填预测方法有频率差异分析技术、叠后波阻抗反演技术、正演模拟技术。频率差异分析技术是利用溶洞充填不同介质时出现的反射强度特征差异，而反射波强弱的变化可以由溶洞壁反射系数及调谐复合反射系数来反映，通过模拟试验溶洞中不同充填物时的调谐反射强度，进而实现溶洞充填物性质的识别。叠后波阻抗反演技术是利用溶洞不同充填条件下的波阻抗差异对充填物性质进行判别，当溶洞内充填砂泥质或未被充填时均表现为低于碳酸盐岩基质的低阻抗特征，通过实钻标定可实现不同充填条件的定性识别。正演模拟技术是采用正演模拟方法对溶洞充填特征进行分析,总结溶洞不同充填特征对应的波形响应特征,在此基础上建立溶洞不同充填情况下的地震反射特征识别方法。在实现本专利技术的过程中，专利技术人发现现有技术至少存在如下问题：现有的频率差异分析、叠后波阻抗反演和正演模拟技术等溶洞充填预测方法，从地震反射强度、波阻抗大小及地震反射特征等方面对溶洞内是否有砂泥质充填进行了预测，形成了溶洞充填的定性识别方法，但现有方法缺乏针对充填特征参数(充填程度、充填物性质)定量表征的有效技术。
技术实现思路
鉴于上述问题，提出了本专利技术以便提供一种克服上述问题或者至少部分地解决上述问题的碳酸盐岩古暗河充填性预测方法。根据本专利技术的一个方面，提供了一种碳酸盐岩古暗河充填性预测方法，该方法包括：分析测井资料，建立溶洞的充填参数的先验分布，以及建立表征充填参数与岩石弹性参数的先验关系的岩石物理模型；根据岩石物理模型计算溶洞的岩石弹性参数；以岩石弹性参数作为已知参数，充填参数作为反演目标，并结合充填参数的先验分布，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数。可选的，岩石弹性参数包括：纵波速度、横波速度和/或密度。可选的，泥质含量、砂质含量、灰质含量以及孔隙度。可选的，求取溶洞的预测充填参数具体为：求取使充填参数的后验概率分布最大的充填参数作为预测充填参数。可选的，在建立表征充填参数与岩石弹性参数的先验关系的岩石物理模型之后，方法还包括：采用蒙特卡洛仿真模拟技术获取先验关系的全区分布，并对充填参数进行抽样；根据岩石物理模型计算溶洞的岩石弹性参数具体为：根据抽样得到的抽样充填参数，并结合岩石物理模型模拟出抽样岩石弹性参数。可选的，充填参数的先验分布具体为：高斯分布、均匀分布或者指数分布。利用本专利技术提供的碳酸盐岩古暗河充填性预测方法，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数，可以获得定量的古暗河充填物性质参数，例如泥质含量、砂质含量、灰质含量及孔隙度，实现暗河内充填物性质及充填程度的定量预测以及准确识别，进而能够指导缝洞型油藏的高效开发。上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，而可依照说明书的内容予以实施，并且为了让本专利技术的上述和其它目的、特征和优点能够更明显易懂，以下特举本专利技术的具体实施方式。附图说明通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本专利技术的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：图1示出了根据本专利技术一个实施例的碳酸盐岩古暗河充填性预测方法的流程图；图2示出了根据本专利技术另一个实施例的碳酸盐岩古暗河充填性预测方法的流程图；图3示出了本专利技术一个具体实施例的纵波速度与实测纵波速度的变化规律图。具体实施方式下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。图1示出了根据本专利技术一个实施例的碳酸盐岩古暗河充填性预测方法的流程图。如图1所示，该方法包括以下步骤：步骤S110：分析测井资料，建立溶洞的充填参数的先验分布，以及建立表征充填参数与岩石弹性参数的先验关系的岩石物理模型。在抽取样本之前，对所要估计的未知参数所了解的信息称为先验信息，而先验分布，即：未知参数的先验信息的分布形式表示。在统计学中，先验信息可以更好的帮助解决统计决策问题。在本实施例中，充填参数的先验分布是指根据测井资料，例如三维地震资料、地质信息以及钻井信息，了解到的关于充填物性质及充填程度的信息的分布形式表示，也即关于溶洞的充填参数的信息的分布形式表示，其中，充填参数包括但不限于以下参数中的至少一项：泥质含量、砂质含量、灰质含量以及孔隙度。同时，根据广泛调研的结果以及现有的测井资料，采集与溶洞的充填参数具有较强的相关关系的岩石弹性参数数据，进而确定能够对充填参数进行预测的岩石弹性参数，并根据溶洞的充填参数与确定的岩石弹性参数的相关关系建立能够表征古暗河的充填参数同岩石弹性参数的先验关系的岩石物理模型。步骤S120：根据岩石物理模型计算溶洞的岩石弹性参数。利用该岩石物理模型，输入充填参数，则可以拟合出对应的岩石弹性参数。具体地，输入抽样充填参数，则可以拟合出抽样岩石弹性参数。其中，采样充填参数可以结合充填参数的先验分布，并以等概率方式接受新的随机抽样点的方法来获取，或者通过随机生成方式来获取。具体地，接受概率分布通常选择均匀分布，以等概率方式接受新的随机抽样点。步骤S130：以岩石弹性参数作为已知参数，充填参数作为反演目标，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数。在贝叶斯反演框架下，后验概率函数与先验分布及似然函数的关系为：其中，P(m|d)为后验概率函数，P(m)为先验分布，P(d|m)为似然函数，P(d)为模型的全模型空间概率，可视为一个常数。对应的，在本实施例中，m为充填参数，d为岩石弹性参数；P(d|m)为在给定充填参数的条件下的条件概率，也称为似然函数，其表达的是充填参数m和岩石弹性参数d之间的相对概率，在贝叶斯反演计算中，它的大小反映了充填参数m和岩石弹性参数d的失配程度；P(m)为对应岩石物理模型的充填参数的先验概率，其涵盖了各种各样的先验信息，包括人们对充填参数的经验认识及其他相关数据对充填参数的约束，其值与岩石弹性参数d无关；P(d)为岩石物理模型的全模型空间的概率，可视为一个常数；P(m|d)为在给定岩石弹性参数的条件下的条件概率，也称为后验概率函数。具体地，将通过岩石物理模型拟合出的抽样岩石弹性参数作为贝叶斯反演中输入的已知参数，而将充填参数作为反演目标，其中，反演目标函数可根据经验选取，本专利技术中，选取最大的后验概率函数为反演目标函数来进行说明，但是并不以此限定反演目标函数的选择。通过反演目标函数，求取溶洞的预测充填参数。由此可见，根据本实施例提供的碳酸盐岩古暗河充填性预本文档来自技高网...

【技术保护点】
1.一种碳酸盐岩古暗河充填性预测方法，其特征在于，包括：分析测井资料，建立溶洞的充填参数的先验分布，以及建立表征所述充填参数与岩石弹性参数的先验关系的岩石物理模型；根据所述岩石物理模型计算溶洞的岩石弹性参数；以岩石弹性参数作为已知参数，充填参数作为反演目标，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数。

【技术特征摘要】
1.一种碳酸盐岩古暗河充填性预测方法，其特征在于，包括：分析测井资料，建立溶洞的充填参数的先验分布，以及建立表征所述充填参数与岩石弹性参数的先验关系的岩石物理模型；根据所述岩石物理模型计算溶洞的岩石弹性参数；以岩石弹性参数作为已知参数，充填参数作为反演目标，进行基于贝叶斯的溶洞的充填参数的反演，求取溶洞的预测充填参数。2.根据权利要求1所述的方法，其特征在于，所述岩石弹性参数包括：纵波速度、横波速度和/或密度。3.根据权利要求1或2所述的方法，所述充填参数包括：泥质含量、砂质含量、灰质含量以及孔隙度。4.根据权利要求1-3任一项所述的方法，...

【专利技术属性】
技术研发人员：鲍典，姜应兵，张慧涛，李正民，赵青，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：北京,11