[发明专利]岩性油气藏弱反射地震信号的重构方法

权利要求书

1.一种岩性油气藏弱反射地震信号的重构方法，所述弱反射地震信号重构方法包括以下步骤：

S1、对原始地震信号进行小波变换的预处理去噪；

S2、将去噪后的信号进行波形匹配端点延拓处理；

所述步骤S2包括以下子步骤：

S21、计算波形匹配度：α(d₀,d_N,P_N)；

其中，d₀是信号中左端第一个最小值点到起始端的波形，其长度记为l，d_N是信号中极大值点到左边起始位置的波形，P_N是信号中的极大值点，令α(d₀,d_N0,P_N0)＝min{α(d₀,d_N,P_N),N＝0,1,2...}；

S22、若α(d₀,d_N0,P_N0)<β°l满足,其中β为一常数，则用从d_N0起左端的第一个极大值和第一个极小值的信号作为原始地震信号的左端延拓；

S23、若不满足上述条件，则取原始信号左端起始处的第一个和第二个极大值点的均值，作为延拓信号的极大值；取左端起始处的第一个和第二个极小值点的均值，作为延拓信号的极小值，进行延拓；

S3、对端点延拓处理后的地震信号进行经验模态分解，得到能够反映各种信号特征的固有模态分量；

S4、对分解得到的固有模态分量与原始地震信号进行相关性计算，将得到的相关系数与阈值相关系数进行对比，自适应选择得到有效固有模态分量；

S5、对有效固有模态分量进行希尔伯特变换得到瞬时振幅、瞬时相位、瞬时频率，将瞬时振幅进行白化处理，并对白化后的瞬时振幅和希尔伯特变换得到的瞬时相位重构信号。

2.根据权利要求1所述的一种岩性油气藏弱反射地震信号的重构方法，其特征在于：所述步骤S1包括以下子步骤：

S11、读取地震记录数据的文件头，并读取地震记录数据文件中的数据部分，得到地震信号的各项信息；

S12、将数据进行矩阵操纵，对原始地震信号进行小波去噪。

3.根据权利要求1所述的一种岩性油气藏弱反射地震信号的重构方法，其特征在于：所述步骤S3包括以下子步骤：

S31、将延拓后信号初始化,找出延拓后信号的所有的极大值和极小值点；

S32、通过三次样条插值函数求上下包络线，计算上下包络的均值,通过包络线的均值得到筛选门限的阈值限定值；

S33、判断阈值限定值是否不处于设定的阈值区间范围内，若不在阈值区间范围内，则有固有模态分量，否则，返回步骤S31,对下一个信号进行循环判断；

S34、判断余量是否为单调函数或是常量，如果是，则余量记为分解剩余项，如果不是，则返回步骤S31,对下一个信号进行循环判断。

4.根据权利要求1所述的一种岩性油气藏弱反射地震信号的重构方法，其特征在于：所述步骤S4包括以下子步骤：

S41、设置有效固有模态分量的阈值相关系数，将经验模态分解后得到的固有模态分量与分解前信号进行相关性比对；

S42、如果比对结果大于阈值相关系数，则该固有模态分量为有效固有模态分量；

S43、如果比对结果小于阈值相关系数，则该固有模态分量为伪分量。

5.根据权利要求1所述的一种岩性油气藏弱反射地震信号的重构方法，其特征在于：所述步骤S5包括以下子步骤：

S51、对筛选出的有效固有模态分量进行希尔伯特变换，变换得到的瞬时振幅、瞬时相位、瞬时频率；

S52、通过傅里叶变换得到有效固有模态分量的瞬时振幅谱，并求出振幅谱的包络，求得白化滤波器；

S53、将每个有效固有模态分量的瞬时振幅与白化滤波器相乘，实现瞬时振幅的自适应增强；

S54、通过步骤S51中的瞬时相位与S53中自适应增强后的瞬时振幅进行重构信号。