一种基于对比函数的块状声波阻抗反演方法、装置及设备与流程

1.本文涉及油气资源勘探地震参数反演技术领域，尤其涉及一种基于对比函数的块状声波阻抗反演方法、装置及设备。


背景技术：

2.声波阻抗是一种非常重要的地震属性，精确的声波阻抗可以用来储层表征和岩相分类，因此声波阻抗反演成为了油气勘探中的一项常规工作流程，由于较高的计算效率，常用确定性反演方法求解声波阻抗，如何获得块状分布且具有较高频带的声波阻抗是本领域密切关注的问题。
3.现有技术中常用声波阻抗反演方法分为两类：1、稀疏脉冲反演方法，需在稀疏约束下精确地反演高分辨率的反射系数。反射系数只有用零范数约束才能得到非常精确的反演结果，因此，通常用柯西分布和拉普拉斯分布等近似稀疏的分布来约束反射系数的求解过程，反演得到的反射系数的统计分布会趋近于指定的某个近似稀疏的分布，这就使得反射系数的稀疏性变弱；或用现有的一范数和二范数的组合来逼近零范数约束，然后采用复杂的混合算法求解这种混合范数目标函数，这使得求解反射系数的过程变得非常复杂。2、基于模型的反演方法，通常需要一个光滑的先验模型约束来增强反问题求解的稳定性，减弱反问题的病态性；同时，还需要加入边缘保护正则化项来保持声波阻抗的块状特性。这个过程同样会导致在目标函数中引入多个正则化参数，从而使得反演过程变得复杂。
4.现在亟需一种高分辨率的块状声波阻抗反演方法，解决现有技术中块状声波阻抗反演精度低且计算量大的问题。


技术实现要素：

5.为解决现有技术中块状声波阻抗反演精度低且计算量大的问题，本文实施例提供了一种基于对比函数的块状声波阻抗反演方法、装置及设备，能够得到较高精度的反演结果、增强反演的稳定性并减小计算量。
6.本文提供了一种基于对比函数的块状声波阻抗反演方法，包括，
7.根据地震数据和测井数据提取地震子波；
8.根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
9.根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
10.根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
11.根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
12.根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
13.本文实施例还提供了一种基于对比函数的块状声波阻抗反演装置，包括，
14.地震子波提取单元，根据地震数据和测井数据提取地震子波；
15.初始声波阻抗模型建立单元，根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
16.相对地震记录计算单元，根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
17.正演算子构造单元，根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
18.声波阻抗对比函数计算单元，根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
19.声波阻抗计算单元，根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
20.本文实施例还提供了一种计算机设备，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述处理器执行所述计算机程序时实现上述的方法。
21.本文实施例还提供了一种计算机存储介质，其上存储有计算机程序，所述计算机程序被计算机设备的处理器运行时，执行上述的方法。
22.利用本文实施例，根据地震子波和初始声波阻抗模型计算相对地震记录，提高了块状声波阻抗反演的分辨率，根据地震数据的最大时间采样点数以及地震子波构造正演算子，根据相对地震记录和正演算子反演的声波阻抗对比函数与真实声波阻抗相比具有更强的稀疏性，因此声波阻抗对比函数的分布更加接近于稀疏分布，不需要给定很强的稀疏约束就能够得到精确的反演结果，提高了反演的精度。此外，本文实施例反演的声波阻抗对比函数与声波阻抗之间具有简单的数学关系，通过由测井数据以及地震构造解释资料建立初始声波阻抗模型和反演的声波阻抗对比函数即可求得声波阻抗，与传统的稀疏脉冲反演方法相比，不需要在反演得到高精度的反射系数序列之后再通过反演求得具有块状特征的声波阻抗，减小了反演的计算量。
附图说明
23.为了更清楚地说明本文实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本文的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
24.图1所示为本文实施例一种基于对比函数的块状声波阻抗反演方法的流程图；
25.图2所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的结构示意图；
26.图3所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的详细结构图；
27.图4所示为本文实施例基于对比函数的块状声波阻抗反演方法的逻辑示意图；
28.图5所示为本文实施例计算机设备的结构示意图；
29.图6所示为本文实施例声波阻抗、声波阻抗对比函数和反射系数对应的频数统计图、核密度估计以及标准高斯分布和标准拉普拉斯分布拟合结果。
30.【附图标记说明】：
31.201、地震子波提取单元；
32.202、初始声波阻抗模型建立单元；
33.203、相对地震记录计算单元；
34.204、正演算子构造单元；
35.205、声波阻抗对比函数计算单元；
36.206、声波阻抗计算单元；
37.301、地震子波提取单元；
38.3011、地震子波提取模块；
39.3012、振幅调整模块；
40.302、初始声波阻抗模型建立单元；
41.3021、地质层位模型构建模块；
42.3022、声波阻抗值计算模块；
43.3023、初始声波阻抗模型建立模块；
44.303、相对地震记录计算单元；
45.3031、初始地震记录计算模块；
46.3032、相对地震记录计算模块；
47.304、正演算子构造单元；
48.3041、差分算子计算模块；
49.3042、子波矩阵构建模块；
50.3043、正演算子构造模块；
51.305、声波阻抗对比函数计算单元；
52.3051、反演目标函数构建模块；
53.3052、声波阻抗对比函数计算模块；
54.306、声波阻抗计算单元；
55.502、计算机设备；
56.504、处理设备；
57.506、存储资源；
58.508、驱动机构；
59.510、输入/输出模块；
60.512、输入设备；
61.514、输出设备；
62.516、呈现设备；
63.518、图形用户接口；
64.520、网络接口；
65.522、通信链路；
66.524、通信总线。
具体实施方式
67.下面将结合本文实施例中的附图，对本文实施例中的技术方案进行清楚、完整地
描述，显然，所描述的实施例仅是本文一部分实施例，而不是全部的实施例。基于本文中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本文保护的范围。
68.如图1所示为文实施例一种基于对比函数的块状声波阻抗反演方法的流程图，在本图中描述了基于对比函数的块状声波阻抗反演的过程，该方法包括：
69.步骤101：根据地震数据和测井数据提取地震子波；
70.步骤102：根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
71.步骤103：根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
72.步骤104：根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
73.步骤105：根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
74.步骤106：根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
75.通过本文实施例的方法，根据地震子波和初始声波阻抗模型计算相对地震记录，提高了块状声波阻抗反演的分辨率，根据地震数据的最大时间采样点数以及地震子波构造正演算子，根据相对地震记录和正演算子反演的声波阻抗对比函数与真实声波阻抗相比具有更强的稀疏性，因此声波阻抗对比函数的分布更加接近于稀疏分布，不需要给定很强的稀疏约束就能够得到精确的反演结果，提高了反演的精度。此外，本文实施例反演的声波阻抗对比函数与声波阻抗之间具有简单的数学关系，通过由测井数据以及地震构造解释资料建立初始声波阻抗模型和反演的声波阻抗对比函数即可求得声波阻抗，与传统的稀疏脉冲反演方法相比，不需要在反演得到高精度的反射系数序列之后再通过反演求得具有块状特征的声波阻抗，减小了反演的计算量。
76.在本文实施例中，所述地震数据为油气厂商记录的地震历史数据，其中包括但不限于地震资料、地震记录，所述测井数据为油气厂商在勘探区域内通过测井仪器在井下测量得到的，所述地震构造解释资料为油气厂商提供的，其中包括但不限于层位标定与构造解释、露头和钻井层位标定、地层接触关系、波组反射特征、层位追踪与断层组合、构造及演化特征。根据油气厂商记录的地震数据和所述油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据采用统计的方法提取叠后的地震数据对应的地震子波。
77.根据本文的一个实施例，为了使提取的地震子波的振幅与实际的地震振幅相一致，进而提高反演的精确度，步骤101根据地震数据和测井数据提取地震子波之后，还包括，根据所述测井数据通过褶积模型进行正演模拟，得到模拟地震记录；将所述模拟地震记录与所述地震数据相对比，计算振幅缩放因子；根据所述振幅缩放因子，调整所述地震子波的振幅与所述地震数据中实际的地震振幅相一致。
78.由于实际的地震振幅是相对值，用统计的方法提取的地震子波采用褶积模型正演模拟得到的地震数据的振幅可能与实际的地震振幅存在细微差异，因此为了使提取的地震子波的振幅与实际的地震振幅相一致，提高反演的精确度，在本步骤中，对油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据通过褶积模型得到模拟地震记录，并与实际的地震数据相对比，判断模拟地震记录与实际的地震数据之间的比例关系，将所述比例关系演化为振幅缩放因子，最后通过计算的振幅缩放因子调整步骤101提取的地震子波的
振幅与地震数据中实际的地震振幅相一致，进而提高反演的精确度。
79.根据本文的一个实施例，为了简化声波阻抗的反演过程，步骤102根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型进一步包括，利用散点差值方法对所述地震构造解释资料进行差值，得到地质层位模型；利用散点差值方法对所述测井数据进行横向差值，得到地下每个点的声波阻抗值；根据所述地质层位模型和所述地下每个点的声波阻抗值建立所述初始声波阻抗模型。
80.在本步骤中，所述地震构造解释资料中包括但不限于各个地层的地质时代、埋藏深度及岩性，利用散点插值法对地震数据中各地质层位的数据分别进行插值，得到地质层位模型，所述地质层位模型中包括各层位的特征和属性，地下每个点的声波阻抗值可以将钻井结果与井的地震剖面连接起来，因此，建立的初始声波阻抗模型可以表示与地震剖面上的反射层对应的具体的地质层位、埋藏深度和岩性。
81.根据本文的一个实施例，步骤103根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录进一步包括，根据所述地震子波以及所述初始声波阻抗模型，通过褶积模型计算初始地震记录；根据所述初始地震记录和所述地震数据通过时间域的散射理论计算相对地震记录。
82.在本步骤中，所述时间域的散射理论类似于相对反射系数，所述相对地震记录可以被定义为真实地震记录(由所述地震数据得到)和初始地震记录的差，如公式(1)所示，求得的相对地震记录被用作反演的输入。
83.s
rela
＝s
real
‑
s
prior
ꢀꢀꢀ
(1)
84.其中，s
rela
为所述相对地震记录，s
real
为所述真实地震记录，s
prior
为所述初始地震记录。
85.根据本文的一个实施例，步骤104根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子的公式为公式(2)和公式(3)：
86.s
rela
＝g
χ
(χ
z
)
ꢀꢀꢀ
(2)
87.f
χ
(χ
z
)＝0.5wd ln(1 χ
z
)
ꢀꢀꢀ
(3)
88.其中，s
rela
表示所述相对地震记录，g
χ
表示所述正演算子，χ
z
表示声波阻抗对比函数，w表示由所述地震子波构成的子波矩阵其中l表示步骤101所提取的地震子波的数量，d表示由所述地震数据s
real
的最大时间采样点数构成的差分算子。
89.在本步骤中，首先需要定义声波阻抗对比函数，假设声波阻抗随时间连续变化，第i个地层界面对应的反射系数r
i
和界面两侧的声波阻抗z
i
和z
i 1
之间的关系可以如公式(4)所示：
[0090][0091]
其中，i表示地层界面序号，r
i
表示第i个地层界面对应的反射系数，z
i
、z
i 1
表示第i
个地层界面和第i 1个地层界面的声波阻抗。
[0092]
对于一个n层水平层状结构，引入差分算子，如公式(5)所示：
[0093][0094]
其中，d表示由所述地震数据的最大时间采样点数构成的差分算子，r
n
×
n
表示差分算子为一个n维实数矩阵，n表示所述地震数据的最大时间采样点数。
[0095]
反射系数和声波阻抗之间的关系用向量的形式可以表示为公式(6)：
[0096]
r＝0.5d ln z
ꢀꢀꢀ
(6)
[0097]
其中，r表示所述反射系数，d表示所述差分算子，z表示声波阻抗。
[0098]
因此，对应于真实声波阻抗和初始声波阻抗模型的反射系数如公式(7)所示：
[0099][0100]
其中，r
real
表示真实反射系数，z
real
表示真实声波阻抗，r
prior
表示初始反射系数，z
prior
表示所述初始声波阻抗模型。
[0101]
参考时间域散射理论的思想，即全波场等于初始波场和散射波场之和，在这里将r
real
视为全波场，r
prior
视为初始波场，定义相对反射系数r
rela
为散射波场，因此所述散射波场、全波场以及所述初始波场的关系为公式(8)：
[0102][0103]
其中，r
rela
表示散射波场，r
real
表示全波场，r
prior
表示初始波场，z
real
表示真实声波阻抗，z
prior
表示所述初始声波阻抗模型。
[0104]
另外，结合相对声波阻抗的定义：所述相对声波阻抗为真实声波阻抗和初始声波阻抗模型之差，如公式(9)所示：
[0105]
z
rela
＝z
real
‑
z
prior
ꢀꢀꢀ
(9)
[0106]
其中，z
rela
表示相对声波阻抗，z
real
表示真实声波阻抗，z
prior
表示所述初始声波阻抗模型。
[0107]
因此，相对声波阻抗包含了真实声波阻抗中的高频成分，由此可知，声波阻抗对比函数的概念如公式(10)所示：
[0108][0109]
其中，χ
z
表示声波阻抗对比函数，z
rela
表示相对声波阻抗，z
real
表示真实声波阻抗，z
prior
表示所述初始声波阻抗模型。
[0110]
由公式(10)可知，声波阻抗对比函数代表了相对声波阻抗相对于初始声波阻抗模型的变化，其中包含了真实声波阻抗中的高频成分，同时，声波阻抗对比函数联系了真实声波阻抗、相对声波阻抗和初始声波阻抗模型，此外，声波阻抗对比函数还表示同层中真实声波阻抗和初始声波阻抗模型的差异。
[0111]
然后推导所述正演算子。将公式(9)和公式(10)代入公式(8)中，可以得到把相对声波阻抗和声波阻抗对比函数联系在一起的相对反射系数的公式如公式(11)所示：
[0112][0113]
为了构建相对地震记录各声波阻抗对比函数之间的关系，引入褶积模型，结合公式(8)和公式(11)，可以将公式(1)表示的相对地震记录推演为公式(12)：
[0114]
s
rela
＝s
real
‑
s
prior
＝0.5wd ln(1 χ
z
)
ꢀꢀꢀ
(12)
[0115]
其中，s
rela
为所述相对地震记录，s
real
为所述真实地震记录，s
prior
为所述初始地震记录，w表示由所述地震子波构成的子波矩阵，d表示由所述地震数据的最大时间采样点数构成的差分算子，χ
z
表示声波阻抗对比函数。
[0116]
因此，公式(12)可以简写为公式(2)和公式(3)，进而得到所述正演算子，所述正演算子将由步骤103得到的相对地震记录与所述声波阻抗对比函数联系在了一起。
[0117]
根据本文的一个实施例，步骤105根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数进一步包括，基于贝叶斯原理，根据所述相对地震记录、所述正演算子以及所述声波阻抗对比函数构建反演目标函数；利用交替方向乘子法求解所述反演目标函数得到所述声波阻抗对比函数。
[0118]
在本步骤中，基于贝叶斯原理，加入对比函数更符合的拉普拉斯分布约束，构建最大后验概率意义下的反演目标函数，得到拉普拉斯约束下的高精度声波阻抗对比函数。首先假设实际地震数据中的噪声满足高斯分布，声波阻抗对比函数满足某个先验分布，似然函数和先验分布可以如公式(13)和公式(14)所示：
[0119][0120][0121]
其中，p(s
rela
|χ
z
)表示似然函数，s
rela
表示所述相对地震记录，χ
z
表示声波阻抗对比函数，为噪声和模型参数的归一化常数，exp表示以自然常数e为底的指数函数，例如exp[f(x)]是e的f(x)次方，g
χ
表示所述正演算子，p(χ
z
)表示先验分布，r(χ
z
)表示声波阻抗对比函数满足某个先验分布。
[0122]
将公式(13)和公式(14)代入到贝叶斯准则中，得到后验概率分布如公式(15)所示：
[0123][0124]
其中，p(s
rela
)表示含噪声的相对地震记录的先验边缘分布，在本文实施例中，因为地震数据是已知量，因此p(s
rela
)为常数，const是p(s
rela
)、和的组合。因此求解公式(15)所示的最大后验概率等价与求解反演目标函数。
[0125]
根据本文的一个实施例，基于贝叶斯原理，根据所述相对地震记录、所述正演算子以及所述声波阻抗对比函数构建反演目标函数的公式如公式(16)所示：
[0126][0127]
其中，c(χ
z
)表示所述反演目标函数，g
χ
表示所述正演算子，s
rela
表示所述相对地震记录，χ
z
表示所述声波阻抗对比函数，r(χ
z
)表示所述声波阻抗对比函数的先验分布，λ表示正则化参数，用于平衡所述声波阻抗对比函数χ
z
的先验分布r(χ
z
)和所述地震数据的误差。
[0128]
对于声波阻抗对比函数反演，假设相对高频的声波阻抗对比函数满足拉普拉斯分布，因此公式(16)中的所述声波阻抗对比函数的先验分布可以写成公式(17)所示的形式：
[0129]
r(χ
z
)＝||χ
z
||1ꢀꢀꢀ
(17)
[0130]
其中，r(χ
z
)表示所述声波阻抗对比函数的先验分布，χ
z
表示所述声波阻抗对比函数，||χ
z
||1表示所述声波阻抗对比函数的一范数。
[0131]
因此，公式(16)与公式(17)结合得到公式(18)：
[0132][0133]
公式(18)是一个一范数约束的反问题，可以利用交替方向乘子法进行高效地求解，得到所述声波阻抗对比函数。
[0134]
将公式(18)写为标准形式，如公式(19)所示：
[0135][0136]
其中，minimize为最小化函数，l(χ
z
)为数据误差，是一个凸函数，其公式如公式(20)所示，y
k
表示所述交替方向乘子法引入的中间变量。
[0137][0138]
其中，l(χ
z
)为数据误差，g
χ
表示如公式(3)所述的正演算子，s
rela
表示如公式(1)所述的相对地震记录。
[0139]
因此，引入拉格朗日乘子可以得到增广最小化问题，如公式(21)所示：
[0140][0141]
其中，σ表示惩罚参数，其中σ＞0，u表示拉格朗日乘子，表示所述声波阻抗对比函数与所述交替方向乘子法引入的中间变量的差值的二范数的平方。
[0142]
在进行迭代过程中，交替方向乘法的更新方式如公式(22)所示：
[0143][0144]
其中，χ
z
、χ
zk 1
表示所述声波阻抗对比函数，k表示迭代次数，表示使得目标函数取得最小值时χ
z
的取值，y
k
、y
k 1
表示所述交替方向乘子法引入的中间变量，s表示软阈值求解算子。
[0145]
基于贝叶斯理论，能够将假定的先验分布约束引入反演过程中，增强反演稳定性。
交替方向乘子法将一个1范数约束的反问题求解替代为迭代地求解两个二范数约束的最小二乘问题，这个过程增强了反演的稳定性，进一步提高了反问题的求解精度，可以为后续的储层预测和岩相分类提供高精度的数据准备。
[0146]
根据本文的一个实施例，步骤106根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗的公式如公式(23)所示：
[0147]
z
real
＝z
prior
(1 χ
z
)
ꢀꢀꢀ
(23)
[0148]
其中，z
real
为声波阻抗，z
prior
为步骤102建立的初始声波阻抗模型，χ
z
表示步骤105反演的声波阻抗对比函数。
[0149]
在本步骤中，由公式(10)的倒数可以得到公式(23)，可见，声波阻抗对比函数和声波阻抗之间具有非常简单的四则运算关系，可以很容易地由公式(8)反演得到的高精度的声波阻抗对比函数和步骤102建立的初始声波阻抗模型反求具有块状特征的声波阻抗。
[0150]
如图2所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的结构示意图，在本图中描述了基于对比函数的块状声波阻抗反演装置的结构，具体包括地震子波提取单元201、初始声波阻抗模型建立单元202、相对地震记录计算单元203、正演算子构造单元204、声波阻抗对比函数计算单元205、声波阻抗计算单元206：
[0151]
地震子波提取单元201根据地震数据和测井数据提取地震子波；
[0152]
初始声波阻抗模型建立单元202根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
[0153]
相对地震记录计算单元203根据地震子波提取单元201得到的震子波、初始声波阻抗模型建立单元202建立的初始声波阻抗模型以及所述地震数据计算相对地震记录；
[0154]
正演算子构造单元204根据所述地震数据的最大时间采样点数以及地震子波提取单元201得到地震子波构造正演算子，所述正演算子用于联系相对地震记录计算单元203计算的相对地震记录与声波阻抗对比函数；
[0155]
声波阻抗对比函数计算单元205根据相对地震记录计算单元203计算的相对地震记录以及正演算子构造单元204构造的正演算子反演所述声波阻抗对比函数；
[0156]
声波阻抗计算单元206根据初始声波阻抗模型建立单元202建立的初始声波阻抗模型和声波阻抗对比函数计算单元205反演的声波阻抗对比函数计算声波阻抗。
[0157]
如图3所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的详细结构图，在本图中具体描述了基于对比函数的块状声波阻抗反演装置的详细结构，包括地震子波提取单元301、初始声波阻抗模型建立单元302、相对地震记录计算单元303、正演算子构造单元304、声波阻抗对比函数计算单元305、声波阻抗计算单元306。
[0158]
根据本文的一个实施例，地震子波提取单元301进一步包括地震子波提取模块3011，根据油气厂商记录的地震数据和所述油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据采用统计的方法提取叠后的地震数据对应的地震子波。
[0159]
根据本文的一个实施例，地震子波提取单元301进一步包括振幅调整模块3012，对油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据通过褶积模型得到模拟地震记录，并与实际的地震数据相对比，判断模拟地震记录与实际的地震数据之间的比例关系，将所述比例关系演化为振幅缩放因子，最后通过计算的振幅缩放因子调整地震子波提取模块3011提取的地震子波的振幅与地震数据中实际的地震振幅相一致，进而提高反演
的精确度。
[0160]
根据本文的一个实施例，初始声波阻抗模型建立单元302进一步包括地质层位模型构建模块3021，利用散点插值法对地震数据中各地质层位的数据分别进行插值，得到地质层位模型，所述地质层位模型中包括各层位的特征和属性。
[0161]
根据本文的一个实施例，初始声波阻抗模型建立单元302进一步包括声波阻抗值计算模块3022，利用散点差值方法对所述测井数据进行横向差值，得到地下每个点的声波阻抗值。
[0162]
根据本文的一个实施例，初始声波阻抗模型建立单元302进一步包括初始声波阻抗模型建立模块3023，根据地质层位模型构建模块3021建立的地质层位模型和声波阻抗值计算模块3022计算的地下每个点的声波阻抗值建立初始声波阻抗模型。声波阻抗值计算模块3022计算的地下每个点的声波阻抗值可以将钻井结果与井的地震剖面连接起来，因此，初始声波阻抗模型建立模块3023建立的初始声波阻抗模型可以表示与地震剖面上的反射层对应的具体的地质层位、埋藏深度和岩性。
[0163]
根据本文的一个实施例，相对地震记录计算单元303进一步包括初始地震记录计算模块3031，根据地震子波提取单元301提取的地震子波以及初始声波阻抗模型建立单元302建立的初始声波阻抗模型，通过褶积模型计算初始地震记录。
[0164]
根据本文的一个实施例，相对地震记录计算单元303进一步包括相对地震记录计算模块3032，根据初始地震记录计算模块3031计算的初始地震记录和所述地震数据通过时间域的散射理论计算相对地震记录。所述时间域的散射理论类似于相对反射系数，相对地震记录计算模块3032计算的相对地震记录可以被定义为真实地震记录(由所述地震数据得到)和初始地震记录的差，如本说明书中公式(1)所示，求得的相对地震记录被用作反演的输入。
[0165]
根据本文的一个实施例，正演算子构造单元304进一步包括差分算子计算模块3041，首先得到所述地震数据的最大时间采样点数，然后利用本说明书中公式(5)计算差分算子。
[0166]
根据本文的一个实施例，正演算子构造单元304进一步包括子波矩阵构建模块3042，根据地震子波提取单元301提取的地震子波以及所提取的地震子波的数量构成子波矩阵，所述子波矩阵可以表示为其中l表示地震子波提取单元301所提取的地震子波的数量。
[0167]
根据本文的一个实施例，正演算子构造单元304进一步包括正演算子构造模块3043，根据差分算子计算模块3041计算的差分算子以及子波矩阵构建模块3042计算的子波矩阵通过本说明书中公式(3)构造正演算子，如本说明书中公式(12)所示，所述正演算子用于联系相对地震记录计算单元303计算的相对地震记录与声波阻抗对比函数。
[0168]
根据本文的一个实施例，声波阻抗对比函数计算单元305进一步包括反演目标函
数构建模块3051，基于贝叶斯原理，根据相对地震记录计算单元303计算的相对地震记录、所述正演算子以及所述声波阻抗对比函数，通过本说明书中公式(13)
‑
(15)构建如本说明书中公式(16)所示的反演目标函数。
[0169]
根据本文的一个实施例，声波阻抗对比函数计算单元305进一步包括声波阻抗对比函数计算模块3052，利用交替方向乘子法通过本说明书中公式(17)
‑
(22)求解反演目标函数构建模块3051构建的反演目标函数得到所述声波阻抗对比函数。
[0170]
根据本文的一个实施例，声波阻抗计算单元306根据初始声波阻抗模型建立单元302建立的初始声波阻抗模型以及声波阻抗对比函数计算单元305计算的声波阻抗对比函数通过本说明书中公式(23)计算声波阻抗。
[0171]
如图4所示为本文实施例基于对比函数的块状声波阻抗反演方法的逻辑示意图，本图中描述了利用地震数据和测井数据反演声波阻抗的逻辑。具体包括：
[0172]
步骤401：提取地震子波。
[0173]
在本步骤中，根据油气厂商记录的地震数据和所述油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据采用统计的方法提取叠后的地震数据对应的地震子波。并对油气厂商在勘探区域内通过测井仪器在井下测量得到的测井数据通过褶积模型得到模拟地震记录，并与实际的地震数据相对比，判断模拟地震记录与实际的地震数据之间的比例关系，将所述比例关系演化为振幅缩放因子，最后通过计算的振幅缩放因子调整提取的地震子波的振幅与地震数据中实际的地震振幅相一致，进而提高反演的精确度。
[0174]
步骤402：建立初始声波阻抗模型。
[0175]
在本步骤中，利用散点差值方法对所述地震构造解释资料进行差值，得到地质层位模型；利用散点差值方法对所述测井数据进行横向差值，得到地下每个点的声波阻抗值；根据所述地质层位模型和所述地下每个点的声波阻抗值建立所述初始声波阻抗模型。所述地震构造解释资料中包括各个地层的地质时代、埋藏深度及岩性，建立的初始声波阻抗模型可以表示与地震剖面上的反射层对应的具体的地质层位、埋藏深度和岩性。
[0176]
步骤403：计算初始地震记录。
[0177]
在本步骤中，根据步骤401提取的地震子波以及步骤402建立的初始声波阻抗模型通过褶积模型计算初始地震记录。
[0178]
步骤404：提取真实地震记录。
[0179]
在本步骤中，提取地震数据中的真实地震记录，本文实施例中可以通过现有技术的方法提取真实地震记录，此处不再赘述。
[0180]
步骤405：计算相对地震记录。
[0181]
在本步骤中，根据步骤403计算的初始地震记录以及步骤404提取的真实地震记录通过时间域的散射理论计算相对地震记录，所述时间域的散射理论类似于相对反射系数，因此，相对地震记录可以被定义为真实地震记录和初始地震记录的差，如本说明书中公式(1)所示。求得的相对地震记录被用作步骤407反演的输入。
[0182]
步骤406：构造正演算子。
[0183]
在本步骤中，根据本说明书中公式(4)
‑
(12)的推导，得到正演算子的公式如本说明书中公式(3)，因此，根据步骤401提取的地震子波构造本说明中公式(3)中的子波矩阵，根据油气厂商记录的地震数据通过现有技术的方法得到所述地震数据的最大时间采样点
数构造本说明书中公式(5)所示的差分算子，之后便可以通过本说明书中公式(3)构造所述正演算子，所述正演算子用于联系步骤405计算的相对地震记录和本说明书中公式(10)所示的声波阻抗对比函数。
[0184]
步骤407：基于贝叶斯原理构造反演目标函数。
[0185]
在本步骤中，基于贝叶斯原理，加入对比函数更符合的拉普拉斯分布约束，通过本说明书中公式(13)
‑
(16)构建最大后验概率意义下的反演目标函数，其中，步骤405计算的相对地震记录、步骤406构造的正演算子以及本说明书公式(10)所示额度声波阻抗对比函数的概念分别用于本说明书中公式(13)
‑
(16)的反演。
[0186]
步骤408：利用交替方向乘子法求解反演目标函数。
[0187]
在本步骤中，通过本说明书中公式(18)
‑
(22)求解步骤407构造的反演目标函数，进而得到声波阻抗对比函数的反演结果。
[0188]
步骤409：计算声波阻抗。
[0189]
在本步骤中，由本说明书中公式(10)的倒数可以得到本说明书中公式(23)，可见，声波阻抗对比函数和声波阻抗之间具有非常简单的四则运算关系，可以很容易地由本说明书中公式(8)反演得到的高精度的声波阻抗对比函数和步骤402建立的初始声波阻抗模型反求具有块状特征的声波阻抗。
[0190]
通过图4描述的方法，根据相对地震记录和正演算子反演的声波阻抗对比函数与真实声波阻抗相比具有更强的稀疏性，因此声波阻抗对比函数的分布更加接近于稀疏分布，不需要给定很强的稀疏约束就能够得到精确的反演结果，提高了反演的精度。基于贝叶斯理论，能够将假定的先验分布约束引入反演过程中，增强反演稳定性。交替方向乘子法将一个一范数约束的反问题求解替代为迭代地求解两个二范数约束的最小二乘问题，这个过程增强了反演的稳定性，进一步提高了反问题的求解精度。此外，反演的声波阻抗对比函数与声波阻抗之间具有简单的数学关系，通过由测井数据以及地震构造解释资料建立初始声波阻抗模型和反演的声波阻抗对比函数即可反求声波阻抗，与传统的稀疏脉冲反演方法相比，不需要在反演得到高精度的反射系数序列之后再通过反演求得具有块状特征的声波阻抗，减小了反演的计算量。
[0191]
如图6所示为本文实施例声波阻抗、声波阻抗对比函数和反射系数对应的频数统计图、核密度估计以及标准高斯分布和标准拉普拉斯分布拟合结果。在本图中，ai表示本文实施例中的声波阻抗，χ
z
表示本文实施例中的声波阻抗对比函数，r表示本文实施例中的反射系数，hist of ai、hist ofχ
z
、hist of r分别表示统计得到的声波阻抗、声波阻抗对比函数和反射系数的频数直方图，ked(kernel density estimation)表示核密度估计结果，sgd(standard gaussian distribution)表示标准高斯分布，sld(standard laplace distribution)表示标准拉普拉斯分布拟合结果。
[0192]
计算所述声波阻抗、声波阻抗对比函数和反射系数的频数和标准高斯分布、标准拉普拉斯分布对应的理论频数之间的卡方值如表1所示：
[0193]
表1
[0194] 标准高斯分布标准拉普拉斯分布反射系数的频数2.44e51.08e3声波阻抗对比函数的频数3.64e20.98e2
声波阻抗的频数0.97e31.37e3
[0195]
其中，卡方值越小说明变量越符合对应的分布，由表1可知，通过本文实施例的方法和装置计算的声波阻抗对比函数更符合比较稀疏的拉普拉斯分布，因此，本说明书根据相对地震记录和正演算子反演的声波阻抗对比函数与真实声波阻抗相比具有更强的稀疏性，同时与比反射系数相比稀疏性较弱，不需要给定很强的稀疏约束就能够得到精确的反演结果，提高了反演的精度。
[0196]
如图5所示为本文实施例计算机设备的结构示意图，本文中的基于对比函数的块状声波阻抗反演装置可以为本实施例中的计算机设备，执行上述本文的方法。计算机设备502可以包括一个或多个处理设备504，诸如一个或多个中央处理单元(cpu)，每个处理单元可以实现一个或多个硬件线程。计算机设备502还可以包括任何存储资源506，其用于存储诸如代码、设置、数据等之类的任何种类的信息。非限制性的，比如，存储资源506可以包括以下任一项或多种组合：任何类型的ram，任何类型的rom，闪存设备，硬盘，光盘等。更一般地，任何存储资源都可以使用任何技术来存储信息。进一步地，任何存储资源可以提供信息的易失性或非易失性保留。进一步地，任何存储资源可以表示计算机设备502的固定或可移除部件。在一种情况下，当处理设备504执行被存储在任何存储资源或存储资源的组合中的相关联的指令时，计算机设备502可以执行相关联指令的任一操作。计算机设备502还包括用于与任何存储资源交互的一个或多个驱动机构508，诸如硬盘驱动机构、光盘驱动机构等。
[0197]
计算机设备502还可以包括输入/输出模块(i/o)510，其用于接收各种输入(经由输入设备512)和用于提供各种输出(经由输出设备514)。一个具体输出机构可以包括呈现设备516和相关联的图形用户接口(gui)518。在其他实施例中，还可以不包括输入/输出模块(i/o)510、输入设备512以及输出设备514，仅作为网络中的一台计算机设备。计算机设备502还可以包括一个或多个网络接口520，其用于经由一个或多个通信链路522与其他设备交换数据。一个或多个通信总线524将上文所描述的部件耦合在一起。
[0198]
通信链路522可以以任何方式实现，例如，通过局域网、广域网(例如，因特网)、点对点连接等、或其任何组合。通信链路522可以包括由任何协议或协议组合支配的硬连线链路、无线链路、路由器、网关功能、名称服务器等的任何组合。
[0199]
本文实施例提供了的计算机设备中处理器执行所述计算机程序时实现如下步骤：
[0200]
根据地震数据和测井数据提取地震子波；
[0201]
根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
[0202]
根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
[0203]
根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
[0204]
根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
[0205]
根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
[0206]
对应于图1、图4中的方法，本文实施例还提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法的步骤。
[0207]
本文实施例还提供一种计算机可读指令，其中当处理器执行所述指令时，其中的程序使得处理器执行如图1、图4所示的方法。
[0208]
应理解，在本文的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本文实施例的实施过程构成任何限定。
[0209]
还应理解，在本文实施例中，术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系。例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
[0210]
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本文的范围。
[0211]
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0212]
在本文所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
[0213]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本文实施例方案的目的。
[0214]
另外，在本文各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0215]
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本文的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本文各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0216]
本文中应用了具体实施例对本文的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本文的方法及其核心思想；同时，对于本领域的一般技术人员，依据本
文的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本文的限制。