虚拟井构建方法及系统技术方案

公开了一种虚拟井构建方法及系统。该方法可以包括：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。本发明专利技术基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。化规律。化规律。

【技术实现步骤摘要】
虚拟井构建方法及系统


[0001]本专利技术涉及地震勘探
，更具体地，涉及一种虚拟井构建方法及系统。

技术介绍

[0002]虚拟井构建(预测)技术是指在没有井的位置处，利用该位置邻域内已有的地质、测井等资料获得测井曲线的一种经济可行的手段。虚拟井可用于后续储层预测的地震反演约束控制，有助于更精细地研究各类油藏参数在空间的变化规律。虚拟井构建技术被认为是一种多快好省的综合技术。
[0003]虚拟井的建立通常有两大类方法，一类方法是基于研究区域内某一空间位置的区域化变量值，从周围邻域内各点的区域化变量值进行估计，这类方法包括反距离平方加权法、残差分析法、克里金法等；另一类构建虚拟井曲线的方法主要基于遗传算法理论，即利用叠前地震入射角道集和初始层速度场数据，在所选控制点上进行基于遗传算法的叠前波形反演，求取出控制点的纵波速度、横波速度和密度曲线。然而，对于第一类方法，欲获得高精度的虚拟井，要求已知的井资料要尽量多，这难以适用于勘探程度低、钻井少的区域。对于第二类方法，在实际工作中，往往没有叠前和初始层速度场数据，叠前数据的质量通常难以得到保证，容易受噪声、处理环节等因素的影响，叠前反演结果的不确定性难以评估。因此，有必要开发一种基于连续时间马尔科夫链的虚拟井构建方法及系统。
[0004]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

技术实现思路

[0005]本专利技术提出了一种虚拟井构建方法及系统，其基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。
[0006]根据本专利技术的一方面，提出了一种虚拟井构建方法。所述方法可以包括：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过所述岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与所述目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
[0007]优选地，所述岩相数据包括：所述岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。
[0008]优选地，所述根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据所述计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函
数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
[0009]优选地，马尔科夫链的表达式为：
[0010][0011]其中，p(π
t
)和p(π
t+1
)为时刻t和时刻t+1的岩相类型的概率，为转移概率矩阵，p(π0)为初始岩相类型的概率。
[0012]优选地，虚拟井波阻抗数据与目标波阻抗的差异的二范数为：
[0013]J＝||I
obs-I
pse
||2ꢀꢀ
(2)
[0014]其中，I
obs
表示目标波阻抗，I
pse
表示虚拟井的波阻抗。
[0015]根据本专利技术的另一方面，提出了一种虚拟井构建系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过所述岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与所述目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
[0016]优选地，所述岩相数据包括：所述岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。
[0017]优选地，所述根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据所述计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
[0018]优选地，马尔科夫链的表达式为：
[0019][0020]其中，p(π
t
)和p(π
t+1
)为时刻t和时刻t+1的岩相类型的概率，为转移概率矩阵，p(π0)为初始岩相类型的概率。
[0021]优选地，虚拟井波阻抗数据与目标波阻抗的差异的二范数为：
[0022]J＝||I
obs-I
pse
||2ꢀꢀ
(2)
[0023]其中，I
obs
表示目标波阻抗，I
pse
表示虚拟井的波阻抗。
[0024]其有益效果在于：
[0025]利用叠后地震数据即可，可不必用叠前地震数据；
[0026]既适用于测井、地质资料丰富的地区，也适用于勘探程度低、测井少的工区，适用性广；
[0027]虚拟井曲线可用于后续的储层描述、地震反演、地震解释等环节；
[0028]本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。
附图说明
[0029]通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。
[0030]图1示出了根据本专利技术的虚拟井构建方法的步骤的流程图。
[0031]图2示出了根据本专利技术的一个实施例的目标位置地震数据的示意图。
[0032]图3示出了根据本专利技术的一个实施例的目标波阻抗的示意图。
[0033]图4a、图4b、图4c、图4d分别示出了根据本专利技术的一个实施例的最终虚拟井的孔隙度、泥质含量、纵波速度、密度曲线示意图。
具体实施方式
[0034]下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。
[0035]图1示出了根据本专利技术的虚拟井构建方法的步骤的流程图。本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种虚拟井构建方法，其特征在于，包括：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过所述岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与所述目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。2.根据权利要求1所述的虚拟井构建方法，其中，所述岩相数据包括：所述岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。3.根据权利要求1所述的虚拟井构建方法，其中，所述根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据所述计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。4.根据权利要求1所述的虚拟井构建方法，其中，马尔科夫链的表达式为：其中，p(π
t
)和p(π
t+1
)为时刻t和时刻t+1的岩相类型的概率，为转移概率矩阵，p(π0)为初始岩相类型的概率。5.根据权利要求1所述的虚拟井构建方法，其中，虚拟井波阻抗数据与目标波阻抗的差异的二范数为：J＝||I
obs-I
pse
||2ꢀꢀ
(2)其中，I
obs
表示目标波阻抗，I
pse
表示虚拟井的波阻抗。6.一种虚拟井构建系统，其特征在于，该系统包括：存储器，存储有计算机...

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