CN116299677A - 岩性圈闭边界的确定方法、装置、电子设备及存储介质 - Google Patents

Abstract

本发明公开了一种岩性圈闭边界的确定方法、装置、电子设备及存储介质。该方法包括：基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；基于岩相约束的弹性参数蒙特卡洛随机模拟，通过所述不同岩性组合的伪井岩相曲线确定目标弹性参数；基于传播矩阵的波动方程正演模拟方法，通过所述目标弹性参数确定不同岩性组合的伪井模拟地震数据；依据所述不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。采用本申请技术方案，避免了传统岩性圈闭边界的确定存在主观性强、不同人员圈定的边界存在较大偏差、客观依据不充分等问题，从而避免了人为因素造成的圈定边界存在偏差的情况。

Description

岩性圈闭边界的确定方法、装置、电子设备及存储介质

技术领域

本发明涉及油气勘探技术领域，尤其涉及一种岩性圈闭边界的确定方法、装置、电子设备及存储介质。

背景技术

深水环境沉积的浊积砂岩储层与周围泥岩形成岩性圈闭，圈闭边界的确定不仅是岩性圈闭有效性判定的主要问题，而且直接影响圈闭的规模。

但通常圈闭边界附近砂体厚度远低于地震分辨率，无法通过地震解释直接获得岩性尖灭点。深水沉积砂体广泛发育，横向相变快，岩性组合更加复杂多样，而深水区勘探钻井数量少，井控范围大，刻画圈闭边界的振幅阈值难以确定，地震响应复杂，手工圈定边界多解性强。因此需要一种能够准确实现岩性圈闭边界定量刻画的技术。

发明内容

本发明提供了一种岩性圈闭边界的确定方法、装置、电子设备及存储介质，以解决圈闭边界识别不确定性强的难题。

根据本发明的一方面，提供了一种岩性圈闭边界的确定方法，该方法包括：

基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；

基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数；

基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据；

依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

根据本发明的另一方面，提供了一种岩性圈闭边界的确定装置，该装置包括：

岩相曲线确定模块，用于基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；

弹性参数确定模块，用于基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数；

地震数据确定模块，用于基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据；

圈闭边界确定模块，用于依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

根据本发明的另一方面，提供了一种电子设备，电子设备包括：

至少一个处理器；以及

与至少一个处理器通信连接的存储器；其中，

存储器存储有可被至少一个处理器执行的计算机程序，计算机程序被至少一个处理器执行，以使至少一个处理器能够执行本发明任一实施例的岩性圈闭边界的确定方法。

根据本发明的另一方面，提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机指令，计算机指令用于使处理器执行时实现本发明任一实施例的岩性圈闭边界的确定方法。

根据本发明实施例的技术方案，通过基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线，使得在对不同岩性组合的伪井进行模拟时，能够得到较为准确的模拟结果，进而减少模拟结果的误差。通过基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数，基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据，使得在计算伪井模拟地震数据时能够便于计算，从而减少计算设备的运算压力。通过依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界，使得最终岩性圈闭边界的确认结果能够尽可能的满足真实情况，提高系统确定结果的准确性。采用上述方法，将圈闭边界刻画从定性阶段发展到定量阶段，解决了圈闭边界识别不确定性强的难题，避免了传统岩性圈闭边界的确定存在主观性强、不同人员圈定的边界存在较大偏差、客观依据不充分等问题，从而避免了人为因素造成的圈定边界存在偏差的情况。

应当理解，本部分所描述的内容并非旨在标识本发明的实施例的关键或重要特征，也不用于限制本发明的范围。本发明的其它特征将通过以下的说明书而变得容易理解。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例一提供了一种岩性圈闭边界的确定方法的流程图；

图2为本发明实施例所使用的不同岩性组合中各层岩性厚度模拟结果的示意图；

图3为本发明实施例二提供的另一种岩性圈闭边界的确定方法的流程图；

图4为本发明实施例所适用的一种模拟岩相的示意图；

图5为本发明实施例所适用的纵波速度的示意图；

图6为本发明实施例所适用的横波速度的示意图；

图7为本发明实施例所适用的密度的示意图；

图8为本发明实施例二提供的另外一种岩性圈闭边界的确定方法的流程图；

图9为本发明实施例所适用的目标岩性的地震振幅的示意图；

图10为本发明实施例所适用的岩性圈闭边界的示意图；

图11为本发明实施例四提供的一种岩性圈闭边界的确定装置的结构示意图；

图12是实现本发明实施例的岩性圈闭边界的确定方法的电子设备的结构示意图。

具体实施方式

为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

实施例一

图1为本发明实施例一提供了一种岩性圈闭边界的确定方法的流程图，本实施例可适用于在缺乏钻井时，实现准确的定量确定岩性圈闭的边界情况，该方法可以由岩性圈闭边界的确定装置来执行，该岩性圈闭边界的确定装置可以采用硬件和/或软件的形式实现，该岩性圈闭边界的确定装置可配置于具有数据处理能力的电子设备中。如图1所示，该方法包括：

S110、基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线。

连续时间马尔科夫链模型可以是用以表示相邻时刻的岩性类别的转移概率，可以利用以下公式(1)进行表示。

p(π_t|π_t+1,π_t+2,…,π_T)＝p(π_t|π_t+1)(1)

式中，π_t为t时刻的岩性，p为不同岩性之间的转移概率矩阵。

岩性类别可以是包括但不限于泥岩、胶结砂岩、灰岩以及砂岩等。伪井垂向岩性组合可以是用以模拟钻井在垂向方向上的岩性组合，其中岩性组合可以是用以表述不同岩石在分布范围上的组合。伪井岩相曲线可以是用以描述模拟钻井的岩性组合。伪井可以是依据连续时间马尔科夫链模型进行模拟得到的模拟钻井。

由于地质过程中把地层剖面视为不同岩石组合的序列，而对某一区域地质剖面往往均符合马尔科夫性，而马尔科夫链具有无后效性，在某一时刻t的岩性的概率分布只与前一时刻t+1的岩性有关，而与其它时刻无关，即

p(π_t|π_t+1,π_t+2,…,π_T)＝p(π_t|π_t+1)

因此可以通过马尔科夫链对岩性在垂向的发育情况进行模拟，依据连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，得到不同岩性组合的伪井的伪井岩相曲线。

在一种可选方案中，根据模拟结果确定不同岩性组合的伪井岩相曲线之前，还包括：

基于指数概率模型对不同岩性组合中各层岩性厚度进行模拟。

图2为本发明实施例所使用的不同岩性组合中各层岩性厚度模拟结果的示意图。参见图2，由于钻井中除了存在各岩性组合的不同之外，相同岩性组合的不同岩性厚度也不相同，因此需要在基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟之后，依据指数概率模型对不同岩性组合中各层岩性厚度进行模拟，从而确定出不同岩性组合的不同岩性厚度的伪井岩相曲线。

其中，指数概率模型为：

λe^-λx

式中，λ参数可以调节概率的均值。

S120、基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数。

目标弹性参数可以是用以进行不同岩性确定的弹性参数。

在得到不同岩性组合的不同岩性厚度的伪井岩相曲线后，能够通过岩相约束的弹性参数蒙特卡洛随机模拟方法进行模拟，从而得到不同岩性组合的不同岩性厚度的目标弹性参数。

S130、基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据。

伪井模拟地震数据可以是通过给定实际地震子波模拟得到的伪井中不同岩性的地震振幅。实际地震子波可以是用以模拟地震状态时使用地震子波。

在得到目标弹性参数后，通过传播矩阵的波动方程正演模拟方法进行模拟，从而确定伪井岩性剖面与地震振幅的定量响应关系，进而得到不同岩性组合的伪井模拟地震数据。

S140、依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

岩性圈闭边界可以是用以描述监测区域内目标岩性所围成的区域的边界。监测区域可以是需要确定岩性圈闭边界的区域。

在得到不同岩性组合的伪井模拟地震数据后，基于伪井模拟地震数据，分别对不同岩性的地震振幅进行统计，从而确定不同岩性的分布，进而确定出岩性圈闭边界。

根据本发明实施例的技术方案，通过基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线，使得在对不同岩性组合的伪井进行模拟时，能够得到较为准确的模拟结果，进而减少模拟结果的误差。通过基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数，基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据，使得在计算伪井模拟地震数据时能够便于计算，从而减少计算设备的运算压力。通过依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界，使得最终岩性圈闭边界的确认结果能够尽可能的满足真实情况，提高系统确定结果的准确性。采用上述方法，将圈闭边界刻画从定性阶段发展到定量阶段，解决了圈闭边界识别不确定性强的难题，避免了传统岩性圈闭边界的确定存在主观性强、不同人员圈定的边界存在较大偏差、客观依据不充分等问题，从而避免了人为因素造成的圈定边界存在偏差的情况。

实施例二

图3为本发明实施例二提供的另一种岩性圈闭边界的确定方法的流程图，本实施例在上述实施例的基础上对前述实施例中基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数的过程进行进一步优化，本实施例可以与上述一个或多个实施例中各个可选方案进行结合。如图3所示，该方法包括：

S210、基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线。

S220、基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定不同岩性组合的候选弹性参数。

其中，候选弹性参数包括纵波速度、横波速度以及密度。

候选弹性参数可以是用以描述不同岩性组合的特征，包括但不限于纵波速度、横波速度以及密度等。

在得到不同岩性组合的伪井岩相曲线后，依据岩相约束的弹性参数蒙特卡洛随机模拟，开展各岩相弹性参数蒙特卡洛随机模拟，从而得到不同岩性的候选弹性参数。

S230、根据不同岩性组合的候选弹性参数确定纵波速度和密度的第一线性回归关系和纵波速度和横波速度的第二线性回归关系。

S240、依据第一线性回归关系和第二线性回归关系对纵波速度进行随机抽样。

S250、根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到目标弹性参数。

图4为本发明实施例所适用的一种模拟岩相的示意图。图5为本发明实施例所适用的纵波速度的示意图。图6为本发明实施例所适用的横波速度的示意图。图7为本发明实施例所适用的密度的示意图。参见图4、5、6以及7，由于得到的候选弹性参数之间具备一定的相关性，因此需要对候选弹性参数之间的关系进行确定。因此将不同岩性组合的候选弹性参数确定纵波速度和密度的第一线性回归关系与建立纵波速度和横波速度的第二线性回归关系，并对纵波速度进行随机抽样，再根据纵波速度的抽样值和第一线性与第二线性回归关系对横波速度和密度进行抽样。

并依据纵波速度、横波速度和密度的抽样值确定目标弹性参数。

在一种可选方案中，根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到目标弹性参数，可包括步骤A1-A2：

步骤A1、根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到抽样弹性参数。

步骤A2、根据抽样弹性参数确定目标弹性参数。

在一种可选方案中，根据如下公式确定目标弹性参数：

其中，EEI(χ)为目标弹性参数，

GI为梯度阻抗，AI₀＝ρ₀V_p0，AI₀为参考纵波阻抗；AI＝ρV_p，AI为纵波阻抗，式中，V_p、V_s和ρ为目标弹性参数中的目标纵波速度、目标横波速度和目标密度；V_p0、V_s0和ρ₀为参考纵波速度、参考横波速度和参考密度，K为横纵波速度比常数，χ为CHI角，角度范围为-90°<χ<90°。

在得到纵波速度、横波速度和密度的抽样值后，可以利用以下公式(2)计算得到目标弹性参数。

其中，EEI(χ)为目标弹性参数，

GI为梯度阻抗，AI₀＝ρ₀V_p0，AI₀为参考纵波阻抗；AI＝ρV_p，AI为纵波阻抗，式中，V_p、V_s和ρ为目标弹性参数中的目标纵波速度、目标横波速度和目标密度；V_p0、V_s0和ρ₀为参考纵波速度、参考横波速度和参考密度，K为横纵波速度比常数，χ为CHI角，角度范围为-90°<χ<90°。

S260、基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据。

S270、依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

根据本发明实施例的技术方案，通过基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定不同岩性组合的候选弹性参数根据不同岩性组合的候选弹性参数确定纵波速度和密度的第一线性回归关系和纵波速度和横波速度的第二线性回归关系，依据第一线性回归关系和第二线性回归关系对纵波速度进行随机抽样，根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到目标弹性参数，使得目标弹性参数的确定能够更为准确，避免了简单抽样时出现不符合候选弹性参数相关性的数据，从而影响系统的运算结果的准确性。

实施例三

图8为本发明实施例二提供的另外一种岩性圈闭边界的确定方法的流程图，本实施例在上述实施例的基础上对前述实施例中依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界的过程进行进一步优化，本实施例可以与上述一个或多个实施例中各个可选方案进行结合。如图8所示，该方法包括：

S310、基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线。

S320、基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数。

S330、基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据。

S340、依据蒙特卡洛模拟确定不同岩性组合的伪井特征弹性曲线。

S350、基于传播矩阵的波动方程正演模拟方法，依据预先确定的实际地震子波进行地震正演模拟，确定目标岩性的地震振幅响应数据；目标岩性至少包括砂岩和围岩。

地震振幅响应数据可以是不同岩性组合的地震振幅数据。目标岩性至少包括砂岩和围岩。

通过蒙特卡洛模拟获得了不同岩性组合的伪井特征弹性曲线，给定实际地震子波，借助基于传播矩阵的波动方程正演模拟方法对伪井开展地震正演模拟，建立伪井岩性剖面与地震振幅的定量响应关系，从而确定目标岩性的地震振幅响应数据。

S360、基于目标岩性的地震振幅响应数据，确定目标岩性的地震振幅概率分布。

S370、依据目标岩性的地震振幅概率分布，确定岩性圈闭边界。

地震振幅概率分布可以是用以描述不同振幅的概率分布情况。

图9为本发明实施例所适用的目标岩性的地震振幅的示意图。参见图9，基于地震振幅响应数据，分别对目标岩性的地震振幅进行统计，获取砂岩段的地震振幅的概率分布以及围岩的地震振幅概率分布。

图10为本发明实施例所适用的岩性圈闭边界的示意图。参见图10，对两个概率分布进行概率函数曲线的拟合，两个概率分布函数曲线的交点处对应的地震振幅值，即为定量区分岩性圈闭尖灭线的地震振幅阈值。以岩性圈闭阈值为截止值，将振幅值低于截止值的区域过滤掉，剩余的振幅异常区即为圈闭范围，振幅异常区的边界即为岩性圈闭边界。

根据本发明实施例的技术方案，通过基于传播矩阵的波动方程正演模拟方法，依据预先确定的实际地震子波进行地震正演模拟，确定目标岩性的地震振幅响应数据，基于目标岩性的地震振幅响应数据，确定目标岩性的地震振幅概率分布，依据目标岩性的地震振幅概率分布，确定岩性圈闭边界，使得岩性圈闭边界的确定结果更为准确，提高了整个系统运算结果的稳定性。

实施例四

图11为本发明实施例四提供的一种岩性圈闭边界的确定装置的结构示意图。本实施例可适用于在缺乏钻井时，实现准确的定量确定岩性圈闭的边界的情况。该岩性圈闭边界的确定装置可以采用硬件和/或软件的形式实现，该岩性圈闭边界的确定装置可配置于具有数据处理能力的电子设备中。如图11所示，本实施例的岩性圈闭边界的确定装置，可包括：岩相曲线确定模块410、弹性参数确定模块420、地震数据确定模块430以及圈闭边界确定模块440。其中：

岩相曲线确定模块410，用于基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；

弹性参数确定模块420，用于基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定目标弹性参数；

地震数据确定模块430，用于基于传播矩阵的波动方程正演模拟方法，通过目标弹性参数确定不同岩性组合的伪井模拟地震数据；

圈闭边界确定模块440，用于依据不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

在上述实施例的基础上，可选的，岩相曲线确定模块410之前，还包括：

厚度模拟确定模块，用于基于指数概率模型对不同岩性组合中各层岩性厚度进行模拟。

在上述实施例的基础上，可选的，弹性参数包括纵波速度、横波速度、密度。

在上述实施例的基础上，可选的，弹性参数确定模块420，包括：

候选弹性确定单元，用于基于岩相约束的弹性参数蒙特卡洛随机模拟，通过不同岩性组合的伪井岩相曲线确定不同岩性组合的候选弹性参数；

回归关系确定单元，用于根据不同岩性组合的候选弹性参数确定纵波速度和密度的第一线性回归关系和纵波速度和横波速度的第二线性回归关系；

随机抽样单元，用于依据第一线性回归关系和第二线性回归关系对纵波速度进行随机抽样；

目标弹性确定单元，用于根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到目标弹性参数。

在上述实施例的基础上，可选的，目标弹性确定单元，包括：

抽样弹性确定子单元，用于根据纵波速度的抽样值和线性回归关系对横波速度和密度进行抽样，得到抽样弹性参数；

目标弹性确定子单元，用于根据抽样弹性参数确定目标弹性参数。

在上述实施例的基础上，可选的，根据如下公式确定目标弹性参数：

其中，EEI(χ)为目标弹性参数，

GI为梯度阻抗，AI₀＝ρ₀V_p0，AI₀为参考纵波阻抗；AI＝ρV_p，AI为纵波阻抗，式中，V_p、V_s和ρ为目标弹性参数中的目标纵波速度、目标横波速度和目标密度；V_p0、V_s0和ρ₀为参考纵波速度、参考横波速度和参考密度，K为横纵波速度比常数，χ为CHI角，角度范围为-90°<χ<90°。

在上述实施例的基础上，可选的，圈闭边界确定模块440，包括：

弹性曲线确定单元，用于依据蒙特卡洛模拟确定不同岩性组合的伪井特征弹性曲线；

响应数据确定单元，用于基于传播矩阵的波动方程正演模拟方法，依据预先确定的实际地震子波进行地震正演模拟，确定目标岩性的地震振幅响应数据；目标岩性至少包括砂岩和围岩；

概率分布确定单元，用于基于目标岩性的地震振幅响应数据，确定目标岩性的地震振幅概率分布；

岩性边界确定单元，用于依据目标岩性的地震振幅概率分布，确定岩性圈闭边界。

本发明实施例所提供的岩性圈闭边界的确定装置可执行本发明任意实施例所提供的岩性圈闭边界的确定方法，具备执行方法相应的功能模块和有益效果。

实施例五

图12示出了可以用来实施本发明的实施例的电子设备10的结构示意图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备(如头盔、眼镜、手表等)和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。

如图12所示，电子设备10包括至少一个处理器11，以及与至少一个处理器11通信连接的存储器，如只读存储器(ROM)12、随机访问存储器(RAM)13等，其中，存储器存储有可被至少一个处理器执行的计算机程序，处理器11可以根据存储在只读存储器(ROM)12中的计算机程序或者从存储单元18加载到随机访问存储器(RAM)13中的计算机程序，来执行各种适当的动作和处理。在RAM 13中，还可存储电子设备10操作所需的各种程序和数据。处理器11、ROM 12以及RAM 13通过总线14彼此相连。输入/输出(I/O)接口15也连接至总线14。

电子设备10中的多个部件连接至I/O接口15，包括：输入单元16，例如键盘、鼠标等；输出单元17，例如各种类型的显示器、扬声器等；存储单元18，例如磁盘、光盘等；以及通信单元19，例如网卡、调制解调器、无线通信收发机等。通信单元19允许电子设备10通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。

处理器11可以是各种具有处理和计算能力的通用和/或专用处理组件。处理器11的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的处理器、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。处理器11执行上文所描述的各个方法和处理，例如岩性圈闭边界的确定方法。

在一些实施例中，岩性圈闭边界的确定方法可被实现为计算机程序，其被有形地包含于计算机可读存储介质，例如存储单元18。在一些实施例中，计算机程序的部分或者全部可以经由ROM 12和/或通信单元19而被载入和/或安装到电子设备10上。当计算机程序加载到RAM 13并由处理器11执行时，可以执行上文描述的岩性圈闭边界的确定方法的一个或多个步骤。备选地，在其他实施例中，处理器11可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行岩性圈闭边界的确定方法。

本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、负载可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。

用于实施本发明的方法的计算机程序可以采用一个或多个编程语言的任何组合来编写。这些计算机程序可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器，使得计算机程序当由处理器执行时使流程图和/或框图中所规定的功能/操作被实施。计算机程序可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。

在本发明的上下文中，计算机可读存储介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的计算机程序。计算机可读存储介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。备选地，计算机可读存储介质可以是机器可读信号介质。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。

为了提供与用户的交互，可以在电子设备上实施此处描述的系统和技术，该电子设备具有：用于向用户显示信息的显示装置(例如，CRT(阴极射线管)或者LCD(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给电子设备。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。

可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(LAN)、广域网(WAN)、区块链网络和互联网。

计算系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，又称为云计算服务器或云主机，是云计算服务体系中的一项主机产品，以解决了传统物理主机与VPS服务中，存在的管理难度大，业务扩展性弱的缺陷。

应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。

上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。

Claims (10)

1.一种岩性圈闭边界的确定方法，其特征在于，包括：

基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；

基于岩相约束的弹性参数蒙特卡洛随机模拟，通过所述不同岩性组合的伪井岩相曲线确定目标弹性参数；

基于传播矩阵的波动方程正演模拟方法，通过所述目标弹性参数确定不同岩性组合的伪井模拟地震数据；

依据所述不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

2.根据权利要求1所述的方法，其特征在于，根据模拟结果确定不同岩性组合的伪井岩相曲线之前，还包括：

基于指数概率模型对不同岩性组合中各层岩性厚度进行模拟。

3.根据权利要求1所述的方法，其特征在于，所述弹性参数包括纵波速度、横波速度、密度。

4.根据权利要求3所述的方法，其特征在于，基于岩相约束的弹性参数蒙特卡洛随机模拟，通过所述不同岩性组合的伪井岩相曲线确定目标弹性参数，包括：

基于岩相约束的弹性参数蒙特卡洛随机模拟，通过所述不同岩性组合的伪井岩相曲线确定不同岩性组合的候选弹性参数；

根据所述不同岩性组合的候选弹性参数确定所述纵波速度和所述密度的第一线性回归关系和所述纵波速度和所述横波速度的第二线性回归关系；

依据所述第一线性回归关系和所述第二线性回归关系对纵波速度进行随机抽样；

根据纵波速度的抽样值和所述线性回归关系对横波速度和密度进行抽样，得到所述目标弹性参数。

5.根据权利要求4所述的方法，其特征在于，根据所述纵波速度的抽样值和所述线性回归关系对所述横波速度和所述密度进行抽样，得到所述目标弹性参数，包括：

根据所述纵波速度的抽样值和所述线性回归关系对所述横波速度和所述密度进行抽样，得到抽样弹性参数；

根据所述抽样弹性参数确定目标弹性参数。

6.根据权利要求5所述的方法，其特征在于，根据如下公式确定所述目标弹性参数：

其中，EEI(χ)为所述目标弹性参数，

GI为梯度阻抗，AI₀＝ρ₀V_p0，AI₀为参考纵波阻抗；AI＝ρV_p，AI为纵波阻抗，式中，V_p、V_s和ρ为所述目标弹性参数中的目标纵波速度、目标横波速度和目标密度；V_p0、V_s0和ρ₀为参考纵波速度、参考横波速度和参考密度，K为横纵波速度比常数，χ为CHI角，角度范围为-90°<χ<90°。

7.根据权利要求1所述的方法，其特征在于，依据所述不同岩性组合的伪井模拟地震数据确定岩性圈闭边界，包括:

依据蒙特卡洛模拟确定所述不同岩性组合的伪井特征弹性曲线；

基于传播矩阵的波动方程正演模拟方法，依据预先确定的实际地震子波进行地震正演模拟，确定目标岩性的地震振幅响应数据；所述目标岩性至少包括砂岩和围岩；

基于所述目标岩性的地震振幅响应数据，确定目标岩性的地震振幅概率分布；

依据所述目标岩性的地震振幅概率分布，确定所述岩性圈闭边界。

8.一种岩性圈闭边界的确定装置，其特征在于，包括：

岩相曲线确定模块，用于基于连续时间马尔科夫链模型对伪井垂向岩性组合进行模拟，根据模拟结果确定不同岩性组合的伪井岩相曲线；

弹性参数确定模块，用于基于岩相约束的弹性参数蒙特卡洛随机模拟，通过所述不同岩性组合的伪井岩相曲线确定目标弹性参数；

地震数据确定模块，用于基于传播矩阵的波动方程正演模拟方法，通过所述目标弹性参数确定不同岩性组合的伪井模拟地震数据；

圈闭边界确定模块，用于依据所述不同岩性组合的伪井模拟地震数据确定岩性圈闭边界。

9.一种电子设备，其特征在于，所述电子设备包括：

至少一个处理器；以及

与所述至少一个处理器通信连接的存储器；其中，

所述存储器存储有可被所述至少一个处理器执行的计算机程序，所述计算机程序被所述至少一个处理器执行，以使所述至少一个处理器能够执行权利要求1-7中任一项所述的岩性圈闭边界的确定方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使处理器执行时实现权利要求1-7中任一项所述的岩性圈闭边界的确定方法。

Images