本发明专利技术公开了一种基于广义Morse标架的地震瞬时属性提取方法，采用广义Morse标架，将有效信号能量分布空间的确定转化成一个优化问题，利用迭代萎缩阈值算法和快速迭代萎缩阈值算法求解。在确定了有效信号能量分布空间的基础上，利用小波变换与Hilbert变换的关系，提出了含噪信号瞬时属性分析的方法。本发明专利技术的方法具有良好的抗噪性能和准确性，得到的瞬时频率剖面能够更加清楚地反映瞬时频率的变化，指示异常区域。
【技术实现步骤摘要】
本专利技术属于地球物理勘探中的信号处理领域，设及地震资料的瞬时属性提取方 法，尤其设及一种。
技术介绍
地震属性是指由地震数据经过数学变换导出的有关地震波的几何形态、运动学特 征、动力学特征和统计学特征，该些特征能够从视觉上反映储层的形态及其含油气性。20世 纪70年代，地震属性分析开始被引入地震解释中。20世纪90年代W来，由于储层描述和= 维数据体解释的需要，地震属性分析技术迅速发展，地震属性在地层构造解释、储层岩性和 物性特征描述、油气藏预测与动态监视方面得到了广泛应用。目前还没有一个公认的地震 属性分类，按照物理意义，地震属性可W分为时间、振幅、频率、相干、衰减等几大类，按照属 性拾取方法，每类又可W分为瞬时属性、单道时窗属性、多道时窗属性。其中瞬时属性是在 地震波到达的位置上拾取的属性，包括瞬时振幅、瞬时相位、瞬时频率、瞬时带宽等等。近年 来，地震属性的数量突增，然而，瞬时属性仍然是地震数据地质解释的支柱，也是目前商业 处理软件必备的技术模块。 在信号处理领域，W瞬时振幅、瞬时相位、瞬时频率为代表的信号瞬时属性的概念 由来已久，在G油or和Vilie的工作之后，大量学者在相关领域做了研究，Boashash对该 些工作做了综述。Taner等人于1979年引入复地震道分析，提出通过复地震道求取瞬时属 性的方法，并给出了该些属性的物理意义W及在地震解释中的应用。瞬时振幅与相邻层的 岩性变化及油气聚集有关。瞬时相位反映界面的不连续性、断层、不整合面和层序边界等。 瞬时频率的变化可W有效刻画地层的厚度和岩性变化，指示油气的分布等。Robedson和 Nogami将瞬时属性用于孔隙砂岩薄层厚度估计，化opra和Marfud利用瞬时属性进行不连 续性、断层和横向不连续性检测。Liu和Mar^d利用瞬时属性检测和刻画曲流河的分布并 确定其厚度。Zeng利用瞬时频率异常来指示薄层。Gao等人利用瞬时频率进行地震资料Q 值估计。 常用的提取地震瞬时属性的方法可W分为W下S类： (1化1化6的变换方法。对地震道做Hilbed变换，转化为复地震道，其中实部为 原地震道数据，虚部为其Hilbed变换。得到复地震道之后，就可W在每个采样点计算振 幅、相位和频率等属性，即瞬时属性。在Taner的工作之后，HUbed变换法被广泛用于计 算地震瞬时属性，至今大部分商业软件仍采用该方法。很多学者对该方法进行了发展和改 进，Barnes于1996年提出二维复地震道分析的概念。Luo等人于2003年提出广义化化ert 变换并给出其在地球物理方面的应用。Barnes于2007年提出加权瞬时频率的概念。Lu和 化ang将HUbed变换方法进行推广，提出基于自适应滤波器计算瞬时频率的方法。 (2)基于时频分析的方法。时频分析方法利用信号的时频分布来求取瞬时属性。 Boashash等人提出基于时频分析的自适应瞬时频率估计方法。Stankovic等人利用自适应 窗时频分布计算瞬时频率。高静怀等人提出了在相空间计算地震瞬时属性的方法。 等人提出基于窄带谱分析的地震资料瞬时属性提取方法。Stee曲s和化ijkoningen提出 基于二次时频分布的地震层序分析W及属性提取方法。化ang等人提出基于经验模式分解 (EMD)计算瞬时频率的方法。Han等人将完全总体经验模式分解（CEEMD)方法用于提取地 震资料的瞬时频率。 做基于反演的方法。Fomel等人提出利用反演的方法来获得局部属性，相比于瞬 时频率，局部频率物理意义更加明确，且应用效果明显。Liu等人提出基于反演的方法，计算 地震资料的瞬时频率。 但是，常用的基于HUbed变换估计瞬时参数的方法对噪声很敏感，且由于滤波 器的截断效应，使得计算出的瞬时属性精度低。
技术实现思路
本专利技术目的在于克服现有技术的不足，提供了一种基于广义Morse标架的地震瞬 时属性提取方法，具有良好的抗噪性能和准确性，得到的瞬时频率剖面能够更加清楚地反 映瞬时频率的变化。 为达到上述目的，本专利技术采用W下技术方案： 一种，包括W下步骤： 步骤1 ;采集地震数据体； 步骤2;计算各道含噪地震数据对应的广义Morse标架系数； 实地震道S(t)的小波变换定义为：(22) 式中基本小波(t)选用广义Morse小波； 广义Morse小波频域表达式为；(23)[001引其中U(W)为单位阶跃函数，0和丫为小波的参数，且0 >0, 丫 >0,ap，Y为 归一化常量，曰e, 丫三2(e丫/e)日/丫； 离散化后的小波族可W表示为：(24)[00过其中a。尺度因子a的离散化步长且a。> 1，t。为平移因子t的离散化步长； 和公式（1)相对应的离散化后的小波变换表示为： Cm,n= <s，IDm,n〉，侦）其中Cm，n为标架系数； 步骤3;对于各道的广义Morse标架系数迭代得到有效信号对应的系数； 定义K为从沪(护)映射到王2 (吸)的算子，将标架系数C= 映射为Z2 (吸） 上的信号，即（26)[002引算子K的伴随算子r为从I，(吸)映射到户(妒)的算子，将Z2 (吸)上的信号投影到 标架系数上： K*s= <s，Vm，n〉， (27)则s=KK*s, (28) 算子K为合成算子，r为分析算子和小波标架{> ' m,。}相对应； 将含噪信号表示为y=s+n=Kx+n, (29) 其中y表示含噪信号，s为不含噪的有效信号，n为高斯白噪声，K为公式（2。中 的合成算子，X表示变换域的系数，通过求解下列优化问题得到有效信号对应的系数宝；(30)[003引其中e和待分析信号的噪声水平有关，用Lagrange乘子A将该问题转化为W下 的无约束问题：(31) 其中A被称为正则化参数，采用下列公式所示的迭代方法对系数X进行更新x(w)=T,，k= 1，2,…，N，（32) 其中为阔值函数，定义为(33)；[00川步骤4;由公式似计算解析信号；[004引（34) 其中Wt(t，a)是有效信号对应的系数，h(t)是s(t)的化化ed变换，c(t)是s(t) 对应的解析信号； 步骤5;利用C(t)计算瞬时振幅、瞬时相位和瞬时频率：和Im分别表示c(t)的实部和虚部。 作为本专利技术进一步优选方案，所述步骤3公式（26)每一次迭代过程中，采用指数 阔值下降策略，如公式（28)所示：(38)[00川其中(39)[005引公式（29)中Am。,和Ami。分别为正则化参数的最小值和最大值，Pm。济Pmi。为最 大和最小百分比。 作为本专利技术进一步优选方案，所述步骤3中为了及时中止迭代过程，定义一个动 态停止准则(40)<当前第1页1 2 3 4 本文档来自技高网...
【技术保护点】
基于广义Morse标架的地震瞬时属性提取方法，其特征在于包括以下步骤：步骤1：采集地震数据体；步骤2：计算各道含噪地震数据对应的广义Morse标架系数；实地震道s(t)的小波变换定义为：Wψ(t,a)≡1a∫-∞+∞s(τ)ψ*(τ-ta)dτ,(1)]]>式中基本小波ψ(t)选用广义Morse小波；广义Morse小波频域表达式为：ψβ,γ(ω)=U(ω)αβ,γωβe-ωγ,(2)]]>其中U(ω)为单位阶跃函数，β和γ为小波的参数，且β＞0，γ＞0，αβ,γ为归一化常量，αβ,γ≡2(eγ/β)β/γ;]]>离散化后的小波族可以表示为：ψm,n(x)=a0-mψ[a0-m(x-nt0)]---(3)]]>其中a0尺度因子a的离散化步长且a0＞1，t0为平移因子t的离散化步长；和公式(1)相对应的离散化后的小波变换表示为：Cm,n＝<s,ψm,n>, (4)其中Cm,n为标架系数；步骤3：对于各道的广义Morse标架系数迭代得到有效信号对应的系数；定义K为从映射到的算子，将标架系数映射为上的信号，即算子K的伴随算子K*为从映射到的算子，将上的信号投影到标架系数上：K*s＝<s,ψ′m,n>, (6)则s＝KK*s, (7)算子K为合成算子，K*为分析算子和小波标架{ψ′m,n}相对应；将含噪信号表示为y＝s+n＝Kx+n, (8)其中y表示含噪信号，s为不含噪的有效信号，n为高斯白噪声，K为公式(22)中的合成算子，x表示变换域的系数，通过求解下列优化问题得到有效信号对应的系数x~=argminx||x||1s.t.||y-Kx||2≤ϵ,---(9)]]>其中ε和待分析信号的噪声水平有关，用Lagrange乘子λ将该问题转化为以下的无约束问题：x~λ=argminx12||y-Kx||22+λ||x||1,---(10)]]>其中λ被称为正则化参数，采用下列公式所示的迭代方法对系数x进行更新x(k+1)＝Tλ[x(k)+K*(y‑Kx(k))],k＝1,2,…,N, (11)其中Tλ为阈值函数，定义为Tλ(x)=x-λx|x||x|≥λ,0|x|<λ.---(12);]]>步骤4：由公式(2)计算解析信号；c(t)≡1Kψ∫0∞Wψ(t,a)daa=s(t)+ih(t),---(13)]]>其中Wψ(t,a)是有效信号对应的系数，h(t)是s(t)的Hilbert变换，c(t)是s(t)对应的解析信号；步骤5：利用c(t)计算瞬时振幅、瞬时相位和瞬时频率：e^(t)=Re2[c](t)+Im2[c](t),---(14)]]>θ^(t)=arctan[Im[c(t)]Re[c(t)]],---(15)]]>f^(t)=12πddt[arctan[Im[c(t)]Re[c(t)]]],---(16)]]>其中Re[c(t)]和Im[c(t)]分别表示c(t)的实部和虚部。...
【技术特征摘要】
【专利技术属性】
技术研发人员：高静怀，王平，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。