一种基于自适应瞬态提取制造技术

本发明专利技术公开了一种基于自适应瞬态提取

【技术实现步骤摘要】
一种基于自适应瞬态提取W变换的致密砂岩薄层表征方法


[0001]本专利技术涉及一种地震信号时间响应分析方法，提出了一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，该方法具有较强的时间定位能力，能够很好地描述薄层的边界结构特征
。

技术介绍

[0002]致密砂岩气是未来
2035
年远景目标内接替常规油气资源最重要和现实的来源
。
致密砂岩储层作为致密砂岩气勘探开发领域的重点和难点，具有纵向薄且窄
、
发育分散
、
物性差
、
非均质性强
、
气水分布等复杂特征，其薄层地震反射信号始终受单层厚度
、
物性差异及流体等复杂因素影响，使常规地震勘探处理解释方法的有效利用受到较大程度的制约
。
[0003]时频分析将一维地震信号映射到二维时频平面，可有效表征隐藏在地震反射信号中的频率特征，并对变化频率进行时间定位，这已广泛应用于储层表征
、
河道检测等
。W
变换是定义了具有对称频率权重高斯窗函数的传统时分析方法，它提高了低频范围内的时频谱分辨率，在近几年受到广泛关注，但在薄储层应用中仍存在分辨率不足的问题
。
瞬态提取算子
(TEO)
作为实现高度时间定位的时频后处理算子，实际应用效果中却存在明显干涉现象，这是因为在
TEO
理论中存在不收敛问题，需要进行对其进行新的界定和限制，这样才能获得较好的储层时间响应特征
。

技术实现思路

[0004]针对现有技术中的上述不足，本专利技术提出一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法
。
该方法中，我们以
W
变换为基础，引入主频估计的
chirp
调制窗口，并定义群延迟估计量下的瞬态提取算子，对其收敛性进行新的界定和限制，从而能够获得更好的薄层时间响应特征
。
[0005]为了实现上述目的，本专利技术采用的技术方案是这样的：一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，包括以下步骤：
[0006](1)
输入叠后地震数据
x(t)
，其中
t
为时间；
[0007](2)
估计
chirp
率
c(t,f)
，其中
f
为频率；
[0008](3)
构造
chirp
调制窗口
w(t,f,c(t,f))
；
[0009](4)
计算具有
chirp
调制窗口的
W
变换
WT(
τ
,f,c)
；
[0010](5)
定义瞬态提取算子其中
τ
为时延参数，为收敛的群延迟估计量；
[0011](6)
计算地震数据的自适应瞬态提取
W
变换
ATEWT(
τ
,f)
；
[0012](7)
进行频率为
f0的特征提取，得到单频剖面
ATEWT(
τ
,f0)。
[0013]作为优选：步骤
(2)
中，具体的
chirp
率
c(t,f)
估计方式为：
[0014][0015]其中
f
为频率，
f0(
τ
)≠0
表示一个波形的主频
。
[0016]作为优选：步骤
(3)
中，具体的
chirp
调制窗口
w(t,f,c(t,f))
构造为：
[0017][0018]其中
w(t,f
；
τ
)
为具有对称频率权重的高斯窗函数，形式可以为：
[0019][0020]其中
k
是一个标量因子，
f0(
τ
)≠0
表示一个波形的主频，
Δ
f(
τ
)
＝
f0(
τ
)
‑
f
主导频率与瞬时频率之差
。
[0021]作为优选：步骤
(4)
中，具有
chirp
调制窗口的
W
变换表达式为：
[0022][0023]由此结果定义瞬态提取算子
。
[0024]作为优选：步骤
(5)
中，定义的瞬态提取算子为：
[0025][0026]其中
τ
为时延参数，具体的群延迟估计量表达式为：
[0027][0028]其中表示虚部，同时，其收敛性体现在对频率的二阶导数和三阶导数的限制条件：
[0029][0030]其中表示对时间的
n
阶导数
。
[0031]作为优选：步骤
(6)
中，定义的自适应瞬态提取
W
变换
ATEWT(
τ
,f)
为：
[0032][0033]由此计算结果进行某频率的特征提取就可以得到单频剖面
。
[0034]与现有技术相比，本专利技术的优点在于：通过引入主频估计的
chirp
调制窗口并定义新的瞬态提取算子，对群延迟估计量的收敛性进行新的界定和限制，提出了一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法
。
该方法能更好地捕捉到地震资料中薄层的时间响应特征，获得更多地震数据信息
。
附图说明
[0035]图1为本专利技术流程图；
[0036]图2为中国四川西部某致密砂岩气田的地震薄层剖面图；
[0037]图3为
W
变换对致密砂岩薄层处理结果：
(a)
为
32Hz
特征提取的单频剖面图，
(b)
为
50Hz
特征提取的单频剖图；
[0038]图4为本专利技术方法对致密砂岩薄层处理结果：
(a)
为
32Hz
特征提取的单频剖面图，
(b)
为
50Hz
特征提取的单频剖图
。
具体实施方式
[0039]下面将结合附图对本专利技术作进一步说明
。
[0040]实施例1：参见图1，一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，包括以下步骤：
[0041](1)
输入叠后地震数据
x(t)
，其中
t...

【技术保护点】

【技术特征摘要】
1.
一种基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，其特征在于：包括以下步骤：
(1)
输入叠后地震数据
x(t)
，其中
t
为时间；
(2)
估计
chirp
率
c(t,f)
，其中
f
为频率；
(3)
构造
chirp
调制窗口
w(t,f,c(t,f))
；
(4)
计算具有
chirp
调制窗口的
W
变换
WT(
τ
,f,c)
；
(5)
定义瞬态提取算子其中
τ
为时延参数，为群延迟估计量；
(6)
计算地震数据的自适应瞬态提取
W
变换
ATEWT(
τ
,f)
；
(7)
进行频率为
f0的特征提取，得到单频剖面
ATEWT(
τ
,f0)。2.
根据权利要求1所述的基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，其特征在于：步骤
(2)
中，具体的
chirp
率
c(t,f)
估计方式为：其中
f
为频率，
f0(
τ
)≠0
表示一个波形的主频
。3.
根据权利要求1所述的基于自适应瞬态提取
W
变换的致密砂岩薄层表征方法，其特征在于：步骤

【专利技术属性】
技术研发人员：陈辉，谢雨桃，陈旭平，王思媛，张骥，张佳，胡英，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：