一种叠后地震波相对波阻抗反演方法及系统.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN113885075A(43)申请公布日2022.01.04(21)申请号202111096121.7(22)申请日2021.09.15(71)申请人成都理工大学地址610000四川省成都市成华区二仙桥东三路1号(72)发明人王兴建薛雅娟曹俊兴廖万平李卿武王崇名(74)专利代理机构济南圣达知识产权代理有限公司37221代理人闫圣娟(51)Int.Cl.G01V1/28(2006.01)G01V1/30(2006.01)G01V1/36(2006.01)权利要求书2页说明书6页附图3页(54)发明名称一种叠后地震波相对波阻抗反演方法及系统(57)摘要本发明公开一种叠后地震波相对波阻抗反演方法及系统，包括：获取目标层位的叠后地震波数据；对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；对重构的叠后地震波数据进行带通积分得到相对波阻抗。无需井位信息，避免了传统波阻抗反演技术受限于地震波频带宽度及需要大量井位信息的缺陷，可适用于低信噪比地震数据。CN113885075ACN113885075A权利要求书1/2页1.一种叠后地震波相对波阻抗反演方法，其特征在于，包括：获取目标层位的叠后地震波数据；对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；对重构的叠后地震波数据进行带通积分得到相对波阻抗。2.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据逐道进行同步挤压小波包变换得到叠后地震波数据的时频能量分布，在时频能量分布中提取预设阈值范围内的模极大值曲线。3.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据x(t)进行同步挤压小波包变换的过程为：2Tx(v,b)＝∫|Wx(p,b)|δ(vx(p,b)‑v)dp；其中，Wx(p,b)是叠后地震波数据x(t)在点(p,b)处的小波包变换，vx(p,b)为叠后地震波数据x(t)在点(p,b)处的局域小波向量估计。4.如权利要求3所述的一种叠后地震波相对波阻抗反演方法，其特征在于，叠后地震波数据x(t)的小波包变换Wx(p,b)为：ss2πi(x‑b)p其中，为母小波的共轭函数，ωpb(t)＝|p|ω(|p|(x‑b))e，参数s∈(1/2,1)，|p|≥1；p,b∈R2；叠后地震波数据x(t)的局域小波向量估计vx(p,b)为：其中，Wx(p,b)≠0。5.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，在时频域内提取预设阈值范围内的模极大值曲线的过程为：对经同步挤压小波包变换得到的时频能量j分布按列优先原则进行矩阵元素的扫描，得到尺度因子f下，区间t∈(t0‑δ,t0+δ)内，满足|jjjjW(f,t)|≤|W(f,t0)|的点(f,t0)，则点(f,t0)为同步挤压小波包变换的局部模极大值点，其中，δ为常数；在预设阈值范围内，将局部模极大值点连接后，得到模极大值曲线。6.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，重构的叠后地震波数据为：其中，为ωpb(t)的双边框架；Uk为第k个聚类子集。7.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，获取目标层位的叠后地震波数据的过程为：综合利用地质信息、测井以及合成地震记录，标定叠后地震波数据的目标层位，确定叠后地震波数据的分析时间范围。2CN113885075A权利要求书2/2页8.一种叠后地震波相对波阻抗反演系统，其特征在于，包括：数据获取模块，被配置为获取目标层位的叠后地震波数据；同步挤压小波包变换模块，被配置为对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；数据重构模块，被配置为根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；反演模块，被配置为对重构的叠后地震波数据进行带通积分得到相对波阻抗。9.一种电子设备，其特征在于，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成权利要求1‑7任一项所述的方法。10.一种计算机可读存储介质，其特征在于，用于存储计算机指令，所述计算机指令被处理器执行时，完成权利要求1‑7任一项所述的方法。3CN113885075A说明书1/6页一种叠后地震波相对波阻抗反演方法及系统技术领域[0001]本发明涉及油气勘探地球物理处理方法