虚拟井构建方法及系统与流程

1.本发明涉及地震勘探技术领域，更具体地，涉及一种虚拟井构建方法及系统。


背景技术：

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


技术实现要素：

5.本发明提出了一种虚拟井构建方法及系统，其基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。
6.根据本发明的一方面，提出了一种虚拟井构建方法。所述方法可以包括：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过所述岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与所述目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
7.优选地，所述岩相数据包括：所述岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。
8.优选地，所述根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据所述计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函
数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
9.优选地，马尔科夫链的表达式为：
[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示出了根据本发明的虚拟井构建方法的步骤的流程图。
[0036]
在该实施例中，根据本发明的虚拟井构建方法可以包括：步骤101，确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；步骤102，根据测井数据确定每个岩相类型对应的岩相数据；步骤103，根据初始岩相类型，通过岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；步骤104，获得每一个虚拟井对应的波阻抗数据；步骤105，分别计算多个虚拟井的波阻抗数据与目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
[0037]
在一个示例中，岩相数据包括：岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。
[0038]
在一个示例中，根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
[0039]
在一个示例中，马尔科夫链的表达式为：
[0040][0041]
其中，p(π
t
)和p(π
t+1
)为时刻t和时刻t+1的岩相类型的概率，为转移概率矩阵，p(π0)为初始岩相类型的概率。
[0042]
在一个示例中，虚拟井波阻抗数据与目标波阻抗的差异的二范数为：
[0043]
j＝||i
obs-i
pse
||2ꢀꢀ
(2)
[0044]
其中，i
obs
表示目标波阻抗，i
pse
表示虚拟井的波阻抗。
[0045]
在一个示例中，岩相参数包括孔隙度、泥质含量、纵波速度、密度。
[0046]
在一个示例中，根据虚拟井的岩相参数中的纵波速度和密度计算虚拟井对应的波
阻抗数据。
[0047]
具体地，研究表明，与常规的离散时间马尔科夫链相比，连续时间马尔科夫链更符合实际地质沉积规律。根据本发明基于连续时间马尔科夫链的的虚拟井构建方法可以包括：
[0048]
确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗。
[0049]
根据测井数据确定每个岩相类型对应的岩相数据，岩相数据包括：岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。根据测井岩相解释结果确定岩相类型的计数矩阵，根据计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
[0050]
根据初始岩相类型，通过岩相数据构建连续时间马尔科夫链，马尔科夫链的表达式为公式(1)，进行多次随机模拟，在随机模拟得到岩相类型后，根据上述确定的岩相类型厚度概率分布函数随机采样可得到岩相类型的厚度值，同时根据岩相参数的概率分布函数随机采样得到岩相参数值，至此便实现连续时间马尔科夫链的构建，即一个虚拟井的构建，进而获得多个虚拟井；根据虚拟井的岩相参数计算虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与目标波阻抗的差异的二范数为公式(2)，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
[0051]
本方法基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。
[0052]
应用示例
[0053]
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0054]
图2示出了根据本发明的一个实施例的目标位置地震数据的示意图。
[0055]
图3示出了根据本发明的一个实施例的目标波阻抗的示意图。
[0056]
确定目标位置，所选地震数据长度为218个采样点，采样间隔为1ms，如图2所示，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗，如图3所示。
[0057]
图4a、图4b、图4c、图4d分别示出了根据本发明的一个实施例的最终虚拟井的孔隙度、泥质含量、纵波速度、密度曲线示意图。
[0058]
选取某实际工区的测井数据，该工区只有砂泥岩两种岩性，则岩相类型间的概率转移矩阵很简单，为分别用高斯函数拟合砂岩厚度和泥岩厚度曲线。给定转移概率矩阵和初始岩相，基于连续时间马尔科夫链随机模拟3000个虚拟井曲线，包括孔隙度、泥质含量、纵波速度、密度，进而分别对3000个虚拟井计算波阻抗曲线；计算3000个模拟波阻抗与反演波阻抗的能量差异，从中选择差异最小的1个虚拟井，作为最终虚拟井，其孔隙度、泥质含量、纵波速度、密度曲线如图4a-d所示。
[0059]
综上所述，本发明基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。
[0060]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0061]
根据本发明的实施例，提供了一种虚拟井构建系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：确定目标位置，对该位置处的地震数据进行叠后波阻抗反演，获得目标波阻抗；根据测井数据确定每个岩相类型对应的岩相数据；根据初始岩相类型，通过岩相数据构建连续时间马尔科夫链，进行多次随机模拟，获得多个虚拟井；获得每一个虚拟井对应的波阻抗数据；分别计算多个虚拟井的波阻抗数据与目标波阻抗的差异的二范数，将差异的二范数最小的虚拟井标记为地震道位置处的最终虚拟井。
[0062]
在一个示例中，岩相数据包括：岩相类型的转移概率矩阵、岩相类型厚度的概率分布函数、每个岩相参数的概率分布函数。
[0063]
在一个示例中，根据测井数据确定每个岩相类型对应的岩相数据包括：根据测井岩相解释结果确定岩相类型的计数矩阵，根据计数矩阵计算岩相类型的转移概率矩阵；统计测井上每个岩相类型的厚度数据，进而通过拟合获得每个岩相类型厚度的概率分布函数；统计测井上每个岩相类型的岩相参数，拟合每个岩相参数的概率分布函数。
[0064]
在一个示例中，马尔科夫链的表达式为：
[0065][0066]
其中，p(π
t
)和p(π
t+1
)为时刻t和时刻t+1的岩相类型的概率，为转移概率矩阵，p(π0)为初始岩相类型的概率。
[0067]
在一个示例中，虚拟井波阻抗数据与目标波阻抗的差异的二范数为：
[0068]
j＝||i
obs-i
pse
||2ꢀꢀ
(2)
[0069]
其中，i
obs
表示目标波阻抗，i
pse
表示虚拟井的波阻抗。
[0070]
在一个示例中，岩相参数包括孔隙度、泥质含量、纵波速度、密度。
[0071]
在一个示例中，根据虚拟井的岩相参数中的纵波速度和密度计算虚拟井对应的波阻抗数据。
[0072]
本系统基于叠后地震数据和工区内已有的测井数据、地质数据，构建虚拟井，适用性广，虚拟井曲线可用于后续的储层描述、高精度地震反演等环节，约束空间变化规律。
[0073]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0074]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。