一种利用薛定谔方程进行地震波衰减梯度估计的方法技术

本发明专利技术属于油气勘探地球物理处理方法领域。本发明专利技术公开了一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法（附图1）。该方法通过将地震数据利用时频分析方法转换到时频域，对时频域的地震数据采用量子力学薛定谔方程求解自适应基波函数，在获得自适应基波函数的映射系数序列的基础上，结合最小二乘法实现衰减梯度的高分辨率估计。本发明专利技术提供了一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法，使用了薛定谔方程进行地震数据的量子特征表征，有利于深度挖掘多频率成分地震信号的内蕴变化关系，将传统衰减梯度估计算法拓展到势能

【技术实现步骤摘要】
一种利用薛定谔方程进行地震波衰减梯度估计的方法


[0001]本专利技术涉及油气勘探地球物理处理方法领域，是一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法。

技术介绍

[0002]地震波属性衰减技术是目前油气勘探领域广泛应用的技术之一，被广泛应用于储层描述和油气检测。地震波衰减梯度估计是地震波属性衰减技术中的一种方法。在传统的地震波动学领域，地震波衰减梯度估计是利用时频分析方法对地震波衰减导致的高频信息损耗进行估计。传统的衰减梯度分析方法采用两点斜率或线性拟合的方法对地震波能量的高频衰减进行估计。为了提高衰减梯度估计的精确性和有效性，高分辨率的时频分析方法是实现这类衰减异常估计方法的关键。常规的地震信号时频分析方法如短时傅里叶变换，S变换，连续小波变换，Wigner
‑
Ville分布等各具特色，但由于受测不准原理的限制，时频分辨率不能同时达到最优，在时空定位上不够准确。如窗函数的固定不变导致短时傅里叶变换在时、频域不能都有足够高的分辨率；S变换由于小波基函数的固定不能满足实际数据处理的需求；小波变换分析方法的有效性依赖于小波函数的选取；而交叉项的存在也限制了Wigner
‑
Ville分布的应用。此外，常规衰减梯度估计方法仅仅适合于地震波频谱较为平滑且地震资料信噪比较高的情况。实际上，由于地震信号是一个多频率成分信号，它的频谱宽带较宽，不同频率成分衰减不同，对整个地震波频谱进行高频段拟合效果并不能够真实反映不同频率成分的高频衰减情况，最终计算结果准确度不够高。
[0003]量子力学是微观世界中最好的物理理论，它的基本假设就是所有事物的不连续能量和薛定谔波动方程。在量子力学中，微观物体的演化是通过薛定谔方程描述的。当我们正确调制薛定谔波动方程的潜在场，我们就可以准确地描述微观物体的演化，然后我们可以准确地描述信号。
[0004]本专利技术的目的在于解决上述现有技术中存在的难题，从现有量子力学薛定谔方程出发，结合地震波传播特性提供一种新的地震波衰减梯度估计方法。采用的是将地震信号通过薛定谔方程转换到势能波函数域，利用自适应基波函数的映射系数挖掘隐藏在地震数据中的物理表现，对地震能量高频衰减进行估计，给出高精度的衰减梯度估计值。利用量子力学薛定谔方程进行地震波衰减梯度估计的方法更适合多频率成分地震波衰减梯度估计，较常规方法具有更强的识别性和准确性，能够提高油气指示的准确性和有效性，计算精度更高、速度快，适合海量地震数据的处理。

技术实现思路

[0005]利用量子力学薛定谔方程进行地震波衰减梯度估计的方法是一种自适应的高分辨率地震能量高频衰减估计技术。利用量子力学薛定谔方程进行地震波衰减梯度估计的方法，包括以下步骤：
[0006](1)综合利用测井、地质和合成地震记录等资料准确标定地震数据的目标层。
[0007](2)针对目标层位地震数据，利用时频分析方法将地震数据转换到时频域。
[0008](3)对时频域的地震数据，逐道逐点对地震数据的频率幅度谱采用薛定谔方程进行求解，获得自适应基波函数，计算波函数在势能
‑
波函数域的映射系数序列。
[0009](4)对势能
‑
波函数域的映射系数序列，求取对数序列，结合最小二乘法计算衰减幅度信息，对计算获得的衰减幅度剖面进行高斯平滑，归一化，获得最终的基于量子力学薛定谔方程的衰减梯度剖面。
[0010]本专利技术的一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法，其核心问题是将地震数据通过时频分析方法转换到时频域，对时频域的地震数据采用量子力学薛定谔方程求解自适应基波函数，在获得自适应基波函数的映射系数序列的基础上，结合最小二乘法实现衰减梯度的高分辨率估计。
[0011]本专利技术的具体实现原理如下：
[0012]1.综合利用测井、地质和合成地震记录等资料准确标定地震数据的目标层。
[0013]2.针对目标层位地震数据，利用时频分析方法将地震数据转换到时频域。记地震信号为x(u)，u为时间变量。经过时频变换后，地震信号的时频谱记为X(f,t)。这里的时频变换方法可以为短时傅里叶变换，S变换，连续小波变换等。其中，
[0014]短时傅里叶变换为：
[0015]X(f,t)＝∫x(u)h
*
(u
‑
t)e
‑
j2πfu
du,
ꢀꢀꢀ
(1)
[0016]式中，h(t)是短时分析窗函数，可以为矩形窗、汉宁窗等。h
*
(
·
)运算表示对窗函数求取共轭函数。t为时间变量，f为频率。
[0017]S变换为：
[0018][0019]连续小波变换为：
[0020]X(a,t)＝∫x(u)ψ
t,a*
(u)du，
ꢀꢀꢀ
(3)
[0021]式中，ψ
*
(
·
)运算表示对母小波ψ(
·
)函数求取共轭函数。a为尺度因子。小波尺度因子a与对应的信号频率具有下列关系式：
[0022][0023]式中，f
s
为采样频率，f
c
为小波中心频率。
[0024]3.对时频域的地震数据，逐道逐点对地震数据的频率幅度谱采用薛定谔方程进行求解，获得自适应基波函数，计算波函数在势能
‑
波函数域的映射系数序列。
[0025]3.1对时频域的地震数据X(f,t)，逐道逐点对地震数据的频率幅度谱采用薛定谔方程进行求解，获得自适应基波函数。
[0026]这里，设地震数据包含m条地震道，每条地震道包含n个采样点。将一道地震信号在i(i＝1,2,L n)点处的频率幅度谱X(f,t
i
)设为一个量子系统的势能V(f)，利用下式对频率幅度谱H(f,t
i
)进行薛定谔方程分解：
[0027][0028]其中，h为普朗克常数，m是量子粒子的质量，E为能量，
▽
表示梯度算子，波函数ψ(f)属于具有可积功能的希尔伯特空间，满足
[0029]∫|ψ(f)|2df＝1。
ꢀꢀꢀ
(6)
[0030]在周期性边界条件下，即ψ(f+L)＝ψ(f)，其中，L为周期长度，方程(5)的所有解可以表示为：
[0031][0032]其中，A表示幅度，每个解即波函数ψ(f)与一个确定的能量E关联，且E＞V。
[0033]构造哈密顿算子则方程(5)的解等同于对角度化哈密顿矩阵H，这里,H是一个n
×
n矩阵。从而有：
[0034][0035]其中，k＝1,2,L,n。H(k,l)表示哈密顿矩阵中第(k,l)个元素。哈密顿矩阵具有如下形式：
[0036][0037]计算哈密顿矩阵H的特征向量Psi，则特征向量Psi的集合构成地震信号在希尔伯特空间上的一组自适应基波函数本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法，其特征在于采用以下步骤：（1）综合利用测井、地质和合成地震记录等资料准确标定地震数据的目标层；（2）针对目标层位地震数据，利用时频分析方法将地震数据转换到时频域；（3）对时频域的地震数据，逐道逐点对地震数据的频率幅度谱采用薛定谔方程进行求解，获得自适应基波函数，计算波函数在势能
‑
波函数域的映射系数序列；（4）对势能
‑
波函数域的映射系数序列，求取对数序列，结合最小二乘法计算衰减幅度信息，对计算获得的衰减幅度剖面进行高斯平滑，归一化，获得最终的基于量子力学薛定谔方程的衰减梯度剖面。2.根据权利要求1所述的一种利用量子力学薛定谔方程进行地震波衰减梯度估计的方法，其特征在于：步骤1中准确标定地...

【专利技术属性】
技术研发人员：王兴建，薛雅娟，曹俊兴，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：