CN106066489A - 基于转换波地震数据构建多薄层目标函数的方法 - Google Patents

Abstract

本发明提供了一种基于转换波地震数据构建多薄层目标函数的方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面，在偏移剖面上建立稀疏约束和波阻抗约束双重约束的多薄层目标函数。本发明的方法能够基于转换波地震数据构建更加合理的多薄层目标函数，进而有利于提高目标函数求解精度和抗噪性，有利于获得高精度反射系数，能够为后续高分辨率剖面重构奠定坚实基础。

Description

基于转换波地震数据构建多薄层目标函数的方法

技术领域

本发明属于转换波地震资料处理技术领域，具体来讲，涉及一种基于转换波地震数据构建多薄层目标函数的方法。

背景技术

目前，相干体裂缝特征描述是探寻裂缝性油气藏的一项重要手段，而获取裂缝预测所需高分辨率数据体、裂缝参数的求取一直是裂缝油气藏勘探中重点研究的问题，尤其转换波低信噪比和分辨率的特点，需要提高其分辨率，便于后续解释的薄层分辨。但是，最小相位和白噪假设为前提的常规高分辨率处理方法，在主频带范围限制下，提高转换波分辨率的能力有限，无法达到后续裂缝预测分辨多薄层最小厚度的要求，所以需要采取新的高分辨率处理思路，即：转换波高精度反射系数求取与高频子波高分辨率重构的新思路，该思路不受任何假设和频带限制，能完成常规手段无法分辨多薄层最小厚度的要求，获取裂缝预测满意的高品质剖面。

通常，反射系数的求取需要构建多薄层目标函数，然而，由于常规方法构建的目标函数存在诸多不合理之处，故而，在求取过程中的反射系数存在极值和位置、大小上的误差，这种误差导致后续高分辨率剖面的错像。

综上，亟需一种能够基于转换波地震数据构建更加合理的多薄层目标函数的方法。

发明内容

本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的之一在于改善现有转换波地震数据领域中为求取反射系数所构建的多薄层目标函数的合理性。

为了实现上述目的，本发明提供一种基于转换波地震数据构建多薄层目标函数的方法。所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面，在转换波偏移剖面上取n道地震记录x(x,t)，x＝1,2,…,n；

在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数：

$O(r_e,r_o)={||\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]||}_2^2+\mathrm{\λ}{||r_e+r_o||}_1+\∂{||\frac{{(Cr-{\mathrm{\ξ}}_p)}^T(Cr-{\mathrm{\ξ}}_p)}{2}||}_1,$

其中，r表示反射系数，a_e代表反射系数的偶分量比例，a_o代表反射系数的奇分量比例，r_e为反射系数的偶分量，r_o为反射系数的奇分量，A为变换矩阵，b为地震数据和子波相关矩阵，λ为稀疏因子，为阻抗因子，C为积分算子矩阵，ξ_p为纵波波阻抗矩阵。

具体来讲，

$\{\begin{array}{c}{r}_e={\[r_e\left(t_1\right),r_e\left(t_2\right),\mathrm{.....},r_e\left(t_N\right)\]}^T\\ r_o={\[r_o\left(t_1\right),r_o\left(t_2\right),\mathrm{.....},r_o\left(t_N\right)\]}^T\end{array},$

$\left\{\begin{array}{c}{b}_e=\left[\begin{array}{c}\mathrm{Re}\[\frac{S(t,f_1)}{W(t,f_1)}{e}^{-i2{\mathrm{\πf}}_1\mathrm{\Δ}t}\]\\ \mathrm{Re}\[\frac{S(t,f_2)}{W(t,f_2)}{e}^{-i2{\mathrm{\πf}}_2\mathrm{\Δ}t}\]\\ \mathrm{...}\\ \mathrm{Re}\[\frac{S(t,f_M)}{W(t,f_M)}{e}^{-i2{\mathrm{\πf}}_M\mathrm{\Δ}t}\]\end{array}\right]\\ b_o=\left[\begin{array}{c}\mathrm{Im}\[\frac{S(t,f_1)}{W(t,f_1)}{e}^{-i2{\mathrm{\πf}}_1\mathrm{\Δ}t}\]\\ \mathrm{Im}\[\frac{S(t,f_2)}{W(t,f_2)}{e}^{-i2{\mathrm{\πf}}_2\mathrm{\Δ}t}\]\\ \mathrm{...}\\ \mathrm{Im}\[\frac{S(t,f_M)}{W(t,f_M)}{e}^{-i2{\mathrm{\πf}}_M\mathrm{\Δ}t}\]\end{array}\right],\end{array}\right.$

$\left\{\begin{array}{c}{A}_e=\left[\begin{array}{cccc}\cos \left({\mathrm{\πT}}_1{f}_1\right)& \cos \left({\mathrm{\πT}}_2{f}_1\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_1\right)\\ \cos \left({\mathrm{\πT}}_1{f}_2\right)& \cos \left({\mathrm{\πT}}_2{f}_2\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ \cos \left({\mathrm{\πT}}_1{f}_M\right)& \cos \left({\mathrm{\πT}}_2{f}_M\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_M\right)\end{array}\right]\\ A_o=\left[\begin{array}{cccc}\sin \left({\mathrm{\πT}}_1{f}_1\right)& \sin \left({\mathrm{\πT}}_2{f}_1\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_1\right)\\ \sin \left({\mathrm{\πT}}_1{f}_2\right)& \sin \left({\mathrm{\πT}}_2{f}_2\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ \sin \left({\mathrm{\πT}}_1{f}_M\right)& \sin \left({\mathrm{\πT}}_2{f}_M\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_M\right)\end{array}\right]\end{array}\right.,$

其中，t为时间，t₁、t₂、……、t_N分别表示第1、2、……、N个时间，f为频率，f₁、f₂、……、f_M分别表示第1、2、……、M个分析频率范围内的频率，T₁、T₂……、T_N/2分别表示第1、2、……、N/2个的时间厚度，Δt表示时移量，Re表示频率域下反射系数R(f)的实部，Im表示频率域下反射系数R(f)的虚部，S表示频率域下地震记录，W表示频率域下地震子波；

其中，t₀表示初始时间，t_M-1表示结束时间；

C的离散形式表示为

3、根据权利要求1所述的基于转换波地震数据构建多薄层目标函数的方法，所述λ的在0至1的范围内取值。

与现有技术相比，本发明的有益效果包括：能够基于转换波地震数据构建更加合理的多薄层目标函数，进而能有效提高目标函数求解精度和抗噪性，有利于获得高精度反射系数，能够为后续高分辨率剖面重构奠定坚实基础。

具体实施方式

在下文中，将结合示例性实施例来详细说明本发明的基于转换波地震数据构建多薄层目标函数的方法。

在本发明的一个示例性实施例中，基于转换波地震数据构建多薄层目标函数的方法可通过以下方式来实现：

(1)对采集到的转换波地震数据，经过常规预处理和叠前偏移后得到的转换波偏移剖面，取n道地震记录x(x,t)，x＝1,2,…,n；其中，x表示道数，t表示时间。

(2)在偏移剖面上建立由稀疏约束和波阻抗约束双重约束的多薄层目标函数；

步骤(2)可按下列方式完成：

(a)对于不含噪声的褶积模型，地震合成记录在频率域中可表示为：

S(f)＝W(f)·R(f) (1)

其中，S(f)表示频率域下地震记录(即，上面的x(x，t))，W(f)表示频率域下地震子波，R(f)表示频率域下反射系数；f为频率，时间t经过傅里叶变换后变换到频率f中。

利用频率域下地震记录与地震子波的丰富信息，利用合适的常规反演算法(如，共轭梯度法、匹配追踪法、模拟退火等)消除子波(即，地震子波)影响(如，子波主频、旁瓣干扰等会影响反演求取反射系数精度的内容)，从记录中得到所需的反射系数(这里，为了得到反射系数R(f))。利用反射系数奇偶分解原理，可得到频率域下目标函数表达式：

$O=\∫\[a_e\{\mathrm{Re}\[\frac{S\left(f\right)}{W\left(f\right)}\]-\mathrm{Re}\[R\left(f\right)\]\}+a_o\{\mathrm{Im}\[\frac{S\left(f\right)}{W\left(f\right)}\]-\mathrm{Im}\[R\left(f\right)\]\}\]df---\left(2\right)$

其中，a_e代表反射系数R(f)的偶分量比例，a_o代表反射系数R(f)的奇分量比例；Re表示R(f)的实部，Im表示R(f)的虚部。

(b)为建立更符合实际地层情况的反射系数模型，对步骤(a)中的多薄层反射系数模型下的目标函数进行推导或改进。

由奇偶反射系数分解原理，得多薄层反射系数奇偶分量可表示为：

$\left\{\begin{array}{c}{g}_e\left(t\right)=\{r_e\left(t\right)\[\mathrm{\δ}(t-\frac{T\left(t\right)}{2})+\mathrm{\δ}(t+\frac{T\left(t\right)}{2})\]dt\\ g_o\left(t\right)=\{r_o\left(t\right)\[\mathrm{\δ}(t-\frac{T\left(t\right)}{2})+\mathrm{\δ}(t+\frac{T\left(t\right)}{2})\]dt\end{array}\right.---\left(3\right)$

r_e(t)代表反射系数r(t)的偶分量，r_o(t)代表反射系数r(t)的奇分量，t表示时间，T表示时间厚度；r(t)为时间域反射系数，r(t)经过傅里叶变换后到上面的频率域反射系数R(f)；g(t)当前分析位置反射系数，g_e(t)代表g(t)的偶分量，g_o(t)代表g(t)的奇分量。

对式(3)进行傅里叶变换，得到其频率域实部与虚部表达式：

$\left\{\begin{array}{c}\mathrm{Re}\[R(t,f)\]=\∫r_e\left(x\right)cos\left(\mathrm{\π}fT\right(t\left)\right)dt\\ \mathrm{Im}\[R(t,f)\]=\∫r_o\left(t\right)\sin \left(\mathrm{\π}fT\right(x\left)\right)dt\end{array}\right.---\left(4\right)$

其中，π为圆周率。

利用式(2)和式(4)得到多薄层反射模型的目标函数表达式为：

$\begin{array}{c}O(r_e,r_o,T,t)=\∫\[a_e\{\mathrm{Re}\[\frac{S(t,f)}{W(t,f)}{e}^{-i2\mathrm{\π}f\mathrm{\Δ}t}\]-\∫r_e\left(t\right)\cos \left(\mathrm{\π}fT\right(t\left)\right)dt\}\\ +a_o\{\mathrm{Im}\[\frac{S(t,f)}{W(t,f)}{e}^{-i2\mathrm{\π}f\mathrm{\Δ}t}\]-\∫r_o\left(t\right)\sin \left(\mathrm{\π}fT\right(t\left)\right)dt\}\]df\end{array}---\left(5\right)$

其中，r_e代表反射系数r(t)的偶分量，r_o代表反射系数r(t)的奇分量，Δt表示时移量；i无物理意义，是指数与正弦函数变化标识符。

为了便于求解，将式(5)可写成矩阵形式：

$O(r_e,r_o)=sum\left(\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]\right)---\left(6\right)$

其中，

其中，t为时间，t₁、t₂、……、t_N分别表示第1、2、……、N个时间，f为频率，f₁、f₂、……、f_M分别表示第1、2、……、M个分析频率范围内的频率，T1、T2……、T_N/2分别表示第1、2、……、N/2个的时间厚度。

(c)由于该目标函数(6)(即，式(6))有多个全局最优解，为了降低求解问题时的多解性和不确定性，需要加上稀疏化约束，将目标函数改为范数求解形式为：

$O(r_e,r_o)={||\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]||}_2^2+\mathrm{\λ}{||r_e+r_o||}_1---\left(7\right)$

其中，λ为稀疏因子，λ一般在0-1范围内取值，||||₁表示L1范数，表示欧几里得(Euclidean)范数的平方。

(d)由于波阻抗可看作是反射系数对时间的积分，因此发明人把它引入到目标函数中作为反射系数r的先验约束。

由于转换波相对于纵波来说其信噪比和分辨率较低,为了提高反射系数精度,本发明使用分辨率更高的纵波波阻抗作为转换波的约束条件。

若I(t)为纵波波阻抗值，则纵波相对波阻抗ε_t为：

$\begin{array}{cc}{\mathrm{\ϵ}}_t=\frac{1}{2}ln\frac{I\left(t\right)}{I\left(t_0\right)}={\∫}_{t_0}^{t}r\left(\mathrm{\η}\right)d\mathrm{\η}& t=t_0,...,t_{M-1}\end{array}---\left(8\right)$

其中，I(t₀)表示纵波初始波阻抗值，r为反射系数，η表示在时间t₀与时间t中取值符合，t₀表示初始时间，t_M-1表示结束时间。

把上式(8)写成矩阵形式：

ξ_p＝Cr (9)

其中C为积分算子矩阵，其离散形式可表示为：

则由最小平方定义波阻抗约束得到新的目标函数为：

$O(r_e,r_o)={||\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]||}_2^2+\mathrm{\λ}{||r_e+r_o||}_1+\∂{||\frac{{(Cr-{\mathrm{\ξ}}_p)}^T(Cr-{\mathrm{\ξ}}_p)}{2}||}_1---\left(11\right)$

其中，r表示反射系数，a_e和a_o分别代表反射系数r(t)的偶奇分量比例，r_e是反射系数的偶分量，r_o是反射系数的奇分量；A为变换矩阵，b是地震数据和子波相关矩阵，λ是稀疏因子，是阻抗因子(其大于零即可)，C为积分算子矩阵，ξ_p是纵波波阻抗矩阵)

另外，可利用频域映射与L1范数联合优化算法来求解式(11)，从而求取高精度反射系数，进而为后续高分辨率剖面重构奠定坚实基础。

本发明的方法能够利用稀疏约束和波阻抗双重约束构建多薄层目标函数，能有效提高目标函数求解精度和抗噪性，有利于获得高精度反射系数，能够为后续高分辨率剖面重构奠定坚实基础。

尽管上面已经结合示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。

Claims (3)

1.一种基于转换波地震数据构建多薄层目标函数的方法，所述方法包括：采集转换波地震数据，进行预处理和叠前偏移后得到转换波偏移剖面，在转换波偏移剖面上取n道地震记录x(x,t)，x＝1,2,…,n；

在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数：

$O(r_e,r_o)=\left|\right|\left[\begin{array}{c}{a}_e(b_e-A_e\×r_e)\\ a_o(b_o-A_o\×r_o)\end{array}\right]|{|}_2^2+\mathrm{\λ}||r_e+r_o|{|}_1+\∂\left|\right|\frac{{(Cr-{\mathrm{\ξ}}_p)}^T(Cr-{\mathrm{\ξ}}_p)}{2}|{|}_1,$

其中，r表示反射系数，a_e代表反射系数的偶分量比例，a_o代表反射系数的奇分量比例，r_e为反射系数的偶分量，r_o为反射系数的奇分量，A为变换矩阵，b为地震数据和子波相关矩阵，λ为稀疏因子，为阻抗因子，C为积分算子矩阵，ξ_p为纵波波阻抗矩阵。

2.根据权利要求1所述的基于转换波地震数据构建多薄层目标函数的方法，其中，

$\{\begin{array}{c}{r}_e={\[r_e\left(t_1\right),r_e\left(t_2\right),\mathrm{.....},r_e\left(t_N\right)\]}^T\\ r_o={\[r_o\left(t_1\right),r_o\left(t_2\right),\mathrm{.....},r_o\left(t_N\right)\]}^T\end{array},$

$\left\{\begin{array}{c}{b}_e=\left[\begin{array}{c}\mathrm{Re}\[\frac{S(t,f_1)}{W(t,f_1)}{e}^{-i2{\mathrm{\πf}}_1\mathrm{\Δ}t}\]\\ \mathrm{Re}\[\frac{S(t,f_2)}{W(t,f_2)}{e}^{-i2{\mathrm{\πf}}_2\mathrm{\Δ}t}\]\\ \mathrm{...}\\ \mathrm{Re}\[\frac{S(t,f_M)}{W(t,f_M)}{e}^{-i2{\mathrm{\πf}}_M\mathrm{\Δ}t}\]\end{array}\right]\\ b_o=\left[\begin{array}{c}\mathrm{Im}\[\frac{S(t,f_1)}{W(t,f_1)}{e}^{-i2{\mathrm{\πf}}_1\mathrm{\Δ}t}\]\\ \mathrm{Im}\[\frac{S(t,f_2)}{W(t,f_2)}{e}^{-i2{\mathrm{\πf}}_2\mathrm{\Δ}t}\]\\ \mathrm{...}\\ \mathrm{Im}\[\frac{S(t,f_M)}{W(t,f_M)}{e}^{-i2{\mathrm{\πf}}_M\mathrm{\Δ}t}\]\end{array}\right],\end{array}\right.$

$\left\{\begin{array}{c}{A}_e=\left[\begin{array}{cccc}\cos \left({\mathrm{\πT}}_1{f}_1\right)& \cos \left({\mathrm{\πT}}_2{f}_1\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_1\right)\\ \cos \left({\mathrm{\πT}}_1{f}_2\right)& \cos \left({\mathrm{\πT}}_2{f}_2\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ \cos \left({\mathrm{\πT}}_1{f}_M\right)& \cos \left({\mathrm{\πT}}_2{f}_M\right)& \mathrm{...}& \cos \left({\mathrm{\πT}}_{N/2}{f}_M\right)\end{array}\right]\\ A_o=\left[\begin{array}{cccc}\sin \left({\mathrm{\πT}}_1{f}_1\right)& \sin \left({\mathrm{\πT}}_2{f}_1\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_1\right)\\ \sin \left({\mathrm{\πT}}_1{f}_2\right)& \sin \left({\mathrm{\πT}}_2{f}_2\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ \sin \left({\mathrm{\πT}}_1{f}_M\right)& \sin \left({\mathrm{\πT}}_2{f}_M\right)& \mathrm{...}& \sin \left({\mathrm{\πT}}_{N/2}{f}_M\right)\end{array}\right]\end{array}\right.,$

其中，t为时间，t₁、t₂、……、t_N分别表示第1、2、……、N个时间，f为频率，f₁、f₂、……、f_M分别表示第1、2、……、M个分析频率范围内的频率，T1、T2……、T_N/2分别表示第1、2、……、N/2个的时间厚度，Δt表示时移量，Re表示频率域下反射系数R(f)的实部，Im表示频率域下反射系数R(f)的虚部，S表示频率域下地震记录x(x,t)，W表示频率域下地震子波；

其中，t₀表示初始时间，t_M-1表示结束时间；

C的离散形式表示为

3.根据权利要求1所述的基于转换波地震数据构建多薄层目标函数的方法，所述λ的在0至1的范围内取值。