CN104200115B - 一种基于地质统计学模拟的全地层速度建模方法 - Google Patents
一种基于地质统计学模拟的全地层速度建模方法 Download PDFInfo
- Publication number
- CN104200115B CN104200115B CN201410465201.9A CN201410465201A CN104200115B CN 104200115 B CN104200115 B CN 104200115B CN 201410465201 A CN201410465201 A CN 201410465201A CN 104200115 B CN104200115 B CN 104200115B
- Authority
- CN
- China
- Prior art keywords
- velocity
- control layer
- speed
- speeds control
- formation velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004088 simulation Methods 0.000 title claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 43
- 238000013508 migration Methods 0.000 claims abstract description 19
- 230000005012 migration Effects 0.000 claims abstract description 19
- 238000005315 distribution function Methods 0.000 claims abstract description 4
- 239000010410 layer Substances 0.000 claims description 64
- 238000011160 research Methods 0.000 claims description 26
- 239000011229 interlayer Substances 0.000 claims description 10
- 238000003325 tomography Methods 0.000 claims description 8
- 208000035126 Facies Diseases 0.000 claims description 7
- 238000004062 sedimentation Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 abstract 1
- 238000005070 sampling Methods 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 102000000584 Calmodulin Human genes 0.000 description 1
- 108010041952 Calmodulin Proteins 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于地质统计学模拟的全地层速度建模方法，所述方法包括以下步骤：确定预设数量的速度控制层；对叠后时间偏移地震数据进行井震标定并对各速度控制层的顶、底界层位进行构造解释，对比追踪获得全区的各速度控制层的顶、底界层位和断层；建立全区的地层框架模型；建立各速度控制层的速度概率密度分布函数和速度纵向变差函数；建立各速度控制层的速度横向变差函数；在所述全区的地层框架模型的控制下，利用所述各速度控制层的速度概率密度分布函数、速度纵向变差函数和速度横向变差函数进行地质统计学模拟，得到多个速度模型实现；对所述多个速度模型实现进行平均计算，得到全地层速度模型。
Description
技术领域
本发明属于石油天然气地震勘探技术领域，更具体地讲，涉及一种应用于石油地震勘探的野外地震资料处理解释的基于地质统计学模拟的全地层速度建模方法。
背景技术
准确建立全地层速度模型是地震资料处理及解释工作的核心任务之一，是获取合理优良地质成果的关键工作内容和重要前提条件。
在地震资料处理工作中，由于各类处理方法理论的假设条件的存在，使得处理人员为了获得高质量的成像而开展速度分析，借此获得研究区的最初叠加速度模型。若将该叠加速度模型直接用作后续叠前时间或叠前深度偏移的初始速度模型，将带来较多的偏移速度迭代次数，这不仅带来繁复的工作量，而且容易造成局部速度迭代收敛方向的偏差。如果初始速度模型存在较大局部偏差或者错误，速度迭代分析也可能无法消除趋势误差，从而降低局部区域偏移成像的质量和精度，甚至在局部区域无法做到成果资料的保真保幅。
在地震资料构造解释工作中，全地层速度模型主要用于时深转换。较为成熟的方法是以处理叠加或偏移速度为背景，利用研究区的测井资料进行约束校正，并以较为大套的地震层位框架进行控制，从而建立全区的层速度模型。该方法能基本满足构造解释的需要，而且也不易发生构造畸变。但是在岩性油气藏及其它各类隐蔽油气藏的勘探任务中，需要确定储层的准确顶底深度，对于水平井的开发部署，储层的准确位置及深度则显得更为重要，粗糙的层速度模型的精度远远无法满足此类勘探任务的需要，而在储层预测的地震反演工作中，对速度模型的精度要求更加不言而喻。随着叠前深度偏移技术的进步与日趋成熟，时深转换工作已经逐渐淡化，速度研究逐渐集中于偏移速度分析和地震反演的速度建模工作之中。
地层速度的本质包含地层框架的控制因素，在偏移速度模型的生成过程 中，多数情况下未加入地层框架控制，使得速度分析的约束控制条件单一化，仅从成像质量和能量聚焦的因素进行修整。尤其对于低信噪比地区，即使加入早期解释地层框架，也极易造成纵向速度样点求取的盲目性和速度横向变化规律的凌乱。
鉴于上述难题，需要从另一个角度开展速度研究并建立全地层速度模型，以在地层框架控制下模拟地层速度并能较为准确地反映地层速度的纵横向变化规律。
发明内容
针对现有技术中存在的问题，本发明的目的是解决上述技术问题中的一个或多个。
本发明的目的在于提供一种引入地层框架的约束以及各速度控制层地质统计学特征并更加合理可靠的基于地质统计学模拟的全地层速度建模方法。
为了实现上述目的，本发明提供了一种基于地质统计学模拟的全地层速度建模方法，所述方法包括以下步骤：利用研究区内的声波测井及VSP测井的速度数据确定预设数量的速度控制层；对叠后时间偏移地震数据进行井震标定并对各速度控制层的顶、底界层位进行构造解释，对比追踪获得全区的各速度控制层的顶、底界层位和断层；利用所述全区的各速度控制层的顶、底界层位和断层建立全区的地层框架模型；利用所述研究区内的声波测井及VSP测井的速度数据建立各速度控制层的速度概率密度分布函数和速度纵向变差函数；在所述构造解释的基础上对各速度控制层进行地震属性研究并获得各速度控制层的地震属性，利用所述各速度控制层的地震属性以及各速度控制层的沉积相研究成果和岩相研究成果建立各速度控制层的速度横向变差函数；在所述全区的地层框架模型的控制下，利用所述各速度控制层的速度概率密度分布函数、速度纵向变差函数和速度横向变差函数进行地质统计学模拟，得到多个速度模型实现；对所述多个速度模型实现进行平均计算，得到全地层速度模型。
根据本发明的基于地质统计学模拟的全地层速度建模方法的一个实施例，所述叠后时间偏移地震数据是对所述研究区的地震数据进行常规预处理及叠后时间偏移处理得到的。
根据本发明的基于地质统计学模拟的全地层速度建模方法的一个实 施例，利用对所述研究区的目的层进行地震反演获得的目的层反演速度对所述全地层速度模型进行内嵌替换，得到目的层速度更新的全地层速度模型。
根据本发明的基于地质统计学模拟的全地层速度建模方法的一个实施例，确定速度控制层的原则是同一个速度控制层的背景速度基本一致并且同一个速度控制层具有基本相同的沉积背景。
本发明的基于地质统计学模拟的全地层速度建模方法能够适用于各类复杂构造及低信噪比的地震工区，所获得的全地层速度模型较为真实合理，符合研究区各地层的地质特征，可用于地震时深转换、确定性反演以及叠前偏移速度迭代的初始速度模型，尤其对隐蔽油气藏的精细刻画与评价，水平井的开发部署等工作具有重要的应用意义。
附图说明
图1是示出根据本发明的示例性实施例的基于地质统计学模拟的全地层速度建模方法的流程图。
具体实施方式
以下，参照附图来详细说明本发明的实施例。
可通过计算机软件模块实现本发明的基于地质统计学模拟的全地层速度建模方法。在本申请中描述的各个步骤不限制为上述步骤，其中的一些步骤可被进一步拆分为更多的步骤，并且一些步骤可合并为较少的步骤。
对于某一研究区的某一地层，在其沉积环境及岩相变化不大的情况下，地层的弹性性质具有较好的地质统计学规律。根据地质统计学理论，结合地层速度的概率密度分布函数和地层框架控制下的纵横向变差函数，可以对地震采样率下的地层样点速度进行统计学模拟，从而建立全区的全地层速度模型。因此，本发明基于以上思路提出了一种基于地质统计学模拟的全地层速度建模方法，对于地层框架内速度样点的插值采用地质统计学模拟算法取代以往的线性插值算法，更加符合地质统计学规律，其横向变差函数将地震数据以及其它地质资料所表现的同一套地层的沉积学规律引入到该速度模型之中，纵向变差函数则将同一套地层的岩性组合的变化规律引入到速度模型中，使获得的速度模型更加符合研究区的地质规律。
图1是示出根据本发明的示例性实施例的基于地质统计学模拟的全地层 速度建模方法的流程图。
参照图1，在步骤S101中，确定速度控制层。根据本发明的示例性实施例，利用研究区内的声波测井及VSP测井的速度数据确定预设数量的速度控制层。其中，速度控制层对应的是这样一套地层，即其地层的背景速度基本一致并且具有基本相同的沉积背景。具体地，利用研究区内钻井的声波测井及VSP测井的速度数据来分析地层速度结构及速度背景，并通过结合研究区的区域地层资料，确定纵向上合适数量的速度控制层。其中，上述的预设数量实际上就是所确定的合适数量。在确定速度控制层时，其基本原则是同一个速度控制层的背景速度基本一致并且同一个速度控制层具有基本相同的沉积背景，从而避免同一个速度控制层内有较大的岩相变化，其中的“基本”的涵义是清楚的，因为速度控制层的背景速度和沉积背景不可能是完全相同的，只要基本一致或基本相同就可以划分为同一个速度控制层。
其中，声波测井数据是指利用声学物理原理制造的测井仪器，由测井电缆下入井内，使地面电测仪沿着井筒连续记录的随深度变化的声学参数，声波测井能测量地层的声波速度。VSP(即垂直地震剖面法)测井数据(或称为VSP地震记录、VSP地震资料、VSP数据)为地面激发井中接收的数据，即炮点设置在地面而检波点设置在井中，其下行直达波受噪声的影响较小，较地面地震记录能更精确地反映地下介质对振幅的衰减作用，其作用之一是能测量地层的地震波速度。
在确定速度控制层之后，在步骤S102中，利用上述声波测井及VSP测井的速度数据对叠后时间偏移地震数据进行井震标定，在步骤S103中，对各速度控制层的顶、底界层位进行构造解释，对比追踪获得全区的各速度控制层的顶、底界层位和断层。其中，叠后时间偏移地震数据是对研究区的地震数据进行常规预处理及叠后时间偏移处理得到的。预处理是采用本领域的常规方法，例如将原始采集的地震数据经过静校正、去噪、振幅补偿、动校正、叠加以及偏移等处理，最终形成叠后或偏移后的地震数据，在此不再赘述。
进行井震标定是利用声波测井制作的合成地震记录和VSP测井处理得到的走廊叠加剖面对地震处理成果剖面进行对比标定，确定各速度控制层顶、底界对应的地震层位位置及其地震响应特性，进行构造解释是根据井震标定的结果分析并获得各速度控制层顶、底界的地震层位位置和其地震响应特征，结合区域地层及构造研究成果或认识，对比追踪获得全区的各速度控制层的 顶、底界层位和断层，井震标定和构造解释均可以采用本领域公知的方法进行。
在步骤S104中，利用全区的各速度控制层的顶、底界层位和断层建立全区的地层框架模型。建立地层框架模型是利用构造解释获得的地震层位和断层，在地震网格控制下进行网格化得到的一种框架模型，其中地震网格是指原始采集的地震数据经过地震处理之后，由地震测网和地震采样率构成的三维网格化空间，建立地层框架模型可以采用本领域公知的方法进行。
在步骤S105和S106中，利用研究区内的声波测井及VSP测井的速度数据建立各速度控制层的速度概率密度分布函数和速度纵向变差函数。
其中，概率密度分布函数是一个描述某个随机变量的输出值在某个确定的取值点附近的可能性的函数，例如，对于一维实随机变量X，设它的累积分布函数是FX(x)，如果存在可测函数fX(x)满足式1：
——式1
那么X是一个连续型随机变量，并且fX(x)就是它的概率密度分布函数。
具体地，概率密度分布函数的含义可以理解为，fX(x)是X落在x处“单位宽度”内的概率，速度的概率密度分布函数描述了速度样点值落在某一速度值附近的概率，某一速度控制层的速度概率密度分布函数的建立方法是利用该速度控制层内的声波测井及VSP测井的速度数据样点进行先验概率密度的计算，再利用某种理论概率密度分布函数对先验概率密度进行匹配，并通过调整理论函数的参数使其与先验概率密度样点达到较好的拟合度。
变差函数是Motheron(1965)年提出的一种矩估计方法，为区域化变量的增量平方的数学期望，也就是区域化变量的增量的方差，它反映了区域化变量在某个方向上某一距离范围内的变化程度，某一速度控制层的速度纵向变差函数描述了速度在某一速度控制层内纵向上的变化程度，速度纵向变差函数建立是利用该速度控制层内的声波测井及VSP测井的速度数据样点进行先验变差值的计算，再利用某种理论变差函数对先验变差值进行匹配，通过调整变差函数的参数使其与先验变差值达到较好的拟合度，其中，纵向变差函数的参数包括函数组分及纵向变程。
在步骤S107中，在构造解释的基础上对各速度控制层进行地震属性研究并获得各速度控制层的地震属性，利用所述各速度控制层的地震属性以及各速度控制层的沉积相研究成果和岩相研究成果建立各速度控制层的速度横向 变差函数。
其中，地震属性主要是以能够反映岩性变化的属性为主，例如振幅属性，它能够反映某一速度控制层速度的横向变化，因此某一速度控制层的速度横向变差函数可以用该速度控制层的地震属性的横向变差函数来表征，根据上述变差函数的涵义，某一速度控制层的速度横向变差函数描述了速度在某一速度控制层内横向上的变化程度，其横向变差函数的建立是利用该速度控制层的地震属性数据样点进行先验变差函数值的计算，再利用某种理论变差函数对先验变差值进行匹配，通过调整变差函数的参数使其与先验变差值达到较好的拟合度，其中，横向变差函数的参数包括函数组分及横向变程，横向变程的确定需要兼顾该速度控制层的沉积相研究成果和岩相研究成果。
其中，在建立各速度控制层的速度横向变差函数时，所利用的地震属性主要是以能够反映岩性变化的属性为主，例如振幅属性等；而各速度控制层的沉积相研究成果和岩相研究成果则是已有的并且可以采用本领域的现有方法获得。利用上述地震属性以及各速度控制层的沉积相研究成果和岩相研究成果就可以确定各速度控制层的速度横向变差函数的参数，如横向变程、函数组成类型等，从而建立适用于研究区各速度控制层的速度横向变差函数。
总之，上述各速度控制层的速度概率密度分布函数、速度纵向变差函数的建立是本领域技术人员可以在各速度控制层的声波测井及VSP测井的速度数据的基础上采用公知的理论和方法获得；各速度控制层的速度横向变差函数的建立亦是本领域技术人员可以在各速度控制层的地震属性研究结果及其沉积相、岩相研究成果的基础上采用公知的理论和方法获得。
在步骤S108中，在全区的地层框架模型的控制下，利用各速度控制层的速度概率密度分布函数、速度纵向变差函数和速度横向变差函数进行地质统计学模拟，得到多个速度模型实现。其中，多个速度模型实现是等概率的并且是通过定义模拟得到。
其中，地质统计学模拟主要基于序贯指示模拟类型算法或马尔科夫链-蒙特卡罗算法，后者更为先进。对于某一速度控制层速度的地质统计学模拟，马尔科夫链-蒙特卡罗算法根据概率密度分布函数获得统计意义上正确的样点集，即根据概率密度分布函数能够得到多个等概率速度组合的结果，纵向变差函数约束马尔科夫链-蒙特卡罗模拟算法的纵向分辨率和概率权重，横向变差函数控制横向未知样点处速度的模拟结果，其控制方式可以是依据已 知样点到未知样点处速度的横向变差，对未知样点处速度的概率密度函数和纵向变差函数的进行统计学定义。地质统计学模拟的实现方式要繁杂许多，但基本原理相同，许多地质统计学著作和论文中已有详细论述，也已产生多个商业化计算机软件模块，在此不再赘述。
在步骤S109中，对多个速度模型实现进行平均计算，得到全地层速度模型。根据上述地质统计学模拟算法的基本原理，模拟算法根据概率密度分布函数能够得到多个等概率速度组合结果，一个等概率速度组合结果即为一个速度模型实现，对多个等概率速度模型实现进行平均计算，可以降低地质统计学模拟的多解性，更加接近真实结果，一般地，定义越多个速度模型实现，其平均化结果越接近真实结果。
如果已对研究区的目的层进行了地震反演工作，则可以利用对研究区的目的层进行地震反演获得的目的层反演速度对上述全地层速度模型进行内嵌替换，从而可以得到目的层速度更新的全地层速度模型。
其中，内嵌替换操作的实现方法可以是，将所述目的层的顶底层位定义为时窗以控制需要进行内嵌替换的数据范围，在时窗范围之内Vnew＝Vinversion，在时窗范围之外Vnew＝Voriginal，其中Voriginal是多个速度模型实现进行平均计算得到的原始的全地层速度模型，Vinversion是目的层地震反演速度体，Vnew即为目的层速度更新的全地层速度模型。
以本发明所阐述的方法或思路建立的全地层速度模型能够较好地体现研究区各地层的地质学特征，保证速度模型趋势的真实性，该全地层速度模型不仅能为各类确定性反演提供合理的初始速度模型，也能作为叠前偏移速度分析的初始输入速度，使偏移速度分析迭代次数减少并向正确的方向迭代收敛，从而提高偏移处理成果资料的保真度。
因此，本发明基于地质统计学模拟的全地层速度建模方法以地质统计学模拟算法为基础建立全地层速度模型，建模过程中地层框架的约束以及各控制层地质统计学特征的引入，使获得的速度模型更加符合研究区的地质规律；目的层反演速度的内嵌替换，能够提高目的层储层深度及形态的刻画精度，本发明比传统的基于线性插值的全地层层速度建模方法更加合理可靠。
尽管已参照优选实施例表示和描述了本发明，但本领域技术人员应该理解，在不脱离由权利要求限定的本发明的精神和范围的情况下，可以对这些实施例进行各种修改和变换。
Claims (4)
1.一种基于地质统计学模拟的全地层速度建模方法，其特征在于，所述方法包括以下步骤：
利用研究区内的声波测井及VSP测井的速度数据来分析地层速度结构及速度背景，并通过结合研究区的区域地层资料，确定预设数量的速度控制层；
对叠后时间偏移地震数据进行井震标定并对各速度控制层的顶、底界层位进行构造解释，对比追踪获得全区的各速度控制层的顶、底界层位和断层；
利用所述全区的各速度控制层的顶、底界层位和断层建立全区的地层框架模型；
利用所述研究区内的声波测井及VSP测井的速度数据建立各速度控制层的速度概率密度分布函数和速度纵向变差函数；
在所述构造解释的基础上对各速度控制层进行地震属性研究并获得各速度控制层的地震属性，利用所述各速度控制层的地震属性以及各速度控制层的沉积相研究成果和岩相研究成果建立各速度控制层的速度横向变差函数；
在所述全区的地层框架模型的控制下，利用所述各速度控制层的速度概率密度分布函数、速度纵向变差函数和速度横向变差函数进行地质统计学模拟，得到多个速度模型实现；
对所述多个速度模型实现进行平均计算，得到全地层速度模型。
2.根据权利要求1所述的基于地质统计学模拟的全地层速度建模方法，其特征在于，所述叠后时间偏移地震数据是对所述研究区的地震数据进行预处理及叠后时间偏移处理得到的。
3.根据权利要求1所述的基于地质统计学模拟的全地层速度建模方法，其特征在于，利用对所述研究区的目的层进行地震反演获得的目的层反演速度对所述全地层速度模型进行内嵌替换，得到目的层速度更新的全地层速度模型。
4.根据权利要求1所述的基于地质统计学模拟的全地层速度建模方法，其特征在于，确定速度控制层的原则是同一个速度控制层的背景速度基本一致并且同一个速度控制层具有基本相同的沉积背景。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410465201.9A CN104200115B (zh) | 2014-09-12 | 2014-09-12 | 一种基于地质统计学模拟的全地层速度建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410465201.9A CN104200115B (zh) | 2014-09-12 | 2014-09-12 | 一种基于地质统计学模拟的全地层速度建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104200115A CN104200115A (zh) | 2014-12-10 |
| CN104200115B true CN104200115B (zh) | 2017-05-17 |
Family
ID=52085408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410465201.9A Active CN104200115B (zh) | 2014-09-12 | 2014-09-12 | 一种基于地质统计学模拟的全地层速度建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104200115B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105093293A (zh) * | 2014-05-14 | 2015-11-25 | 中国石油天然气股份有限公司 | 低频补偿提高洞穴型碳酸盐岩储层地震量化预测的方法 |
| CN104612660B (zh) * | 2015-01-23 | 2017-04-26 | 中南大学 | 一种基于mrf的油气储层岩相随机模拟及实现方法 |
| CN105353412B (zh) * | 2015-12-14 | 2017-08-29 | 中国石油大学(华东) | 一种井震联合平均速度场的计算方法及系统 |
| CN105842736B (zh) * | 2016-05-27 | 2018-07-10 | 中国石油集团东方地球物理勘探有限责任公司 | 一种构建地层速度模型的方法 |
| CN106094027B (zh) * | 2016-06-01 | 2019-03-12 | 中国海洋石油集团有限公司 | 一种垂直地震剖面vsp钻前压力预测方法和系统 |
| CN107065013B (zh) * | 2016-12-27 | 2019-01-25 | 中国石油大学(华东) | 一种地震尺度下的层速度确定方法及装置 |
| CN108318921B (zh) * | 2018-01-11 | 2021-08-10 | 中国石油大学(华东) | 一种基于横向约束的快速地震随机反演方法 |
| CN108646294B (zh) * | 2018-06-29 | 2019-09-06 | 中海石油(中国)有限公司 | 一种复杂断裂系统下深水浊积砂岩储层高精度反演方法 |
| CN109212610B (zh) * | 2018-10-31 | 2021-06-22 | 中国石油化工股份有限公司 | 基于欧式距离近似匹配的井震数据提取方法 |
| CN110579802B (zh) * | 2019-10-09 | 2021-02-26 | 中国科学院海洋研究所 | 一种天然气水合物储层物性参数的高精度反演方法 |
| CN115113266B (zh) * | 2021-03-19 | 2024-08-27 | 中国石油化工股份有限公司 | 火山岩深度域相控速度建模方法、存储介质和计算机设备 |
| CN114114399B (zh) * | 2021-11-26 | 2023-06-23 | 同济大学 | 一种基于Bayes统计决策的智能速度谱解释与建模方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101796507A (zh) * | 2007-08-27 | 2010-08-04 | 兰德马克绘图国际公司，哈里伯顿公司 | 计算变差函数模型的系统和方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ATE533075T1 (de) * | 2007-12-21 | 2011-11-15 | Prad Res & Dev Ltd | Mehrpunkt-geostatistikverfahren mit abzweigungs- lauflängen-kompression und lokaler rastertransformation |
-
2014
- 2014-09-12 CN CN201410465201.9A patent/CN104200115B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101796507A (zh) * | 2007-08-27 | 2010-08-04 | 兰德马克绘图国际公司，哈里伯顿公司 | 计算变差函数模型的系统和方法 |
Non-Patent Citations (2)
| Title |
|---|
| 地质统计学反演在三道桥工区储层预测中的应用;宁松华 等;《石油天然气学报（江汉石油学院学报）》;20140731;第36卷(第7期);第52-54，IV页 * |
| 巴楚夏河工区的三维速度建模方法;武丽 等;《石油地球物理勘探》;20060228;第41卷(第1期);第87-92页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104200115A (zh) | 2014-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104200115B (zh) | 一种基于地质统计学模拟的全地层速度建模方法 | |
| EP0891562B1 (en) | 3-d geologic modelling | |
| US9696442B2 (en) | Method for estimating subsurface properties from geophysical survey data using physics-based inversion | |
| US7062072B2 (en) | Methods of producing images of underground formations surrounding a borehole | |
| Høyer et al. | Multiple-point statistical simulation for hydrogeological models: 3-D training image development and conditioning strategies | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| US11740376B2 (en) | Reflection seismology multiple imaging | |
| CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
| CN108680951A (zh) | 一种基于地震信息判断煤层气富集沉积控制作用的方法 | |
| Yasin et al. | Estimation of petrophysical parameters from seismic inversion by combining particle swarm optimization and multilayer linear calculator | |
| WO2017035104A1 (en) | Velocity model seismic static correction | |
| Borgomano et al. | Stratigraphic well correlations for 3-D static modeling of carbonate reservoirs | |
| US20220155480A1 (en) | Secondary recovery surveillance using validated streamline-based simulation | |
| US10451769B2 (en) | Method for determining petrophysical properties from logging measurements | |
| Ayodele et al. | Static reservoir modeling using stochastic method: A case study of the cretaceous sequence of Gamtoos Basin, Offshore, South Africa | |
| Moret et al. | Investigating the stratigraphy of an alluvial aquifer using crosswell seismic traveltime tomography | |
| CN111077578B (zh) | 岩层分布预测方法和装置 | |
| Liu et al. | Geostatistical modeling for fine reservoir description of Wei2 block of Weicheng oilfield, Dongpu depression, China | |
| Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
| Mehdipour et al. | The Best Scenario for Geostatistical Modeling of Porosity in the Sarvak Reservoir in an Iranian Oil Field, Using Electrofacies, Seismic Facies, and Seismic Attributes | |
| Volpi et al. | Vertical proportion curves: a qualitative and quantitative tool for reservoir characterization | |
| CN113806674A (zh) | 古河道纵向尺度的量化方法、装置、电子设备及存储介质 | |
| Masoud et al. | Reservoir characterization and geostatistical model of the Cretaceous and Cambrian-Ordovician reservoir intervals, Meghil field, Sirte basin, Libya | |
| CN110941029B (zh) | 一种涉及地质戴帽的速度建模方法 | |
| Correia | Integration of reservoir characterization with history matching guided by pilot wells: Application to the Norne Field |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180202 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200916 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Co-patentee after: BGP Inc., China National Petroleum Corp. Patentee after: CHINA NATIONAL PETROLEUM Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |
|
| TR01 | Transfer of patent right |