非稳态盲反褶积方法及装置.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN113219530A(43)申请公布日2021.08.06(21)申请号202010080341.X(22)申请日2020.02.05(71)申请人中国石油天然气股份有限公司地址100007北京市东城区东直门北大街9号(72)发明人李红兵李勇根高浩洋袁三一董世泰(74)专利代理机构北京三友知识产权代理有限公司11127代理人赵平周永君(51)Int.Cl.G01V1/30(2006.01)G01V1/28(2006.01)G01V1/36(2006.01)权利要求书2页说明书10页附图11页(54)发明名称非稳态盲反褶积方法及装置(57)摘要本发明提供了一种非稳态盲反褶积方法及装置，所述方法包含：获取待分析的地震数据，于所述地震数据中预设一初始地震子波；通过预设品质因子范围和所述初始地震子波构建衰减子波群，并对所述衰减子波群进行反射系数反演，获得反演后的反射系数；将所述反射系数代入L1-2范数函数后求二阶导数，获得求导曲线；最小化所述求导曲线中L1-2范数高阶导数，获得品质因子取值及对应的反射系数反演结果。CN113219530ACN113219530A权利要求书1/2页1.一种非稳态盲反褶积方法，其特征在于，所述方法包含：获取待分析的地震数据，于所述地震数据中预设一初始地震子波；通过预设品质因子范围和所述初始地震子波构建衰减子波群，并对所述衰减子波群进行反射系数反演，获得反演后的反射系数；将所述反射系数代入L1-2范数函数后求二阶导数，获得求导曲线；最小化所述求导曲线对应的L1-2范数函数的二阶导数，获得品质因子取值及对应的反射系数反演结果。2.根据权利要求1所述的非稳态盲反褶积方法，其特征在于，所述L1-2范数函数包含：在上式中，0＜σ＜1为权重因子，N为时间采样点数，r为反演的反射系数序列，ri为r中的第i个元素，Q为品质因子，J(Q)为所求的稀疏准则L1-2范数函数值。3.根据权利要求1所述的非稳态盲反褶积方法，其特征在于，将反演后的反射系数代入L1-2范数函数后求二阶导数，获得求导曲线包含：将反演后的反射系数代入L1-2范数函数获得扫描曲线；对所述扫描曲线求取二阶导数获得对应的二阶曲线；将所述扫描曲线和所述二阶曲线分别进行光滑处理后再叠加处理，获得求导曲线。4.根据权利要求1所述的非稳态盲反褶积方法，其特征在于，对所述衰减子波群进行反射系数反演，获得反演后的反射系数包含：利用稀疏贝叶斯学习算法对所述衰减子波群进行反射系数反演，获得反演后的反射系数。5.根据权利要求4所述的非稳态盲反褶积方法，其特征在于，所述衰减子波群包含：上式中，M为频率采样点数，FH为反傅里叶算子，W为初始子波经傅里叶变换后构成的Toeplitz矩阵，H为共轭转置，s表示非稳态地震记录，A为与品质因子有关的时频域衰减矩阵，r为反射系数序列，G(Q)为时变子波矩阵。6.一种非稳态盲反褶积装置，其特征在于，所述装置包含采集模块、反演模块、计算模块和分析模块；所述采集模块用于获取待分析的地震数据，于所述地震数据中预设一初始地震子波；所述反演模块用于通过预设品质因子范围和所述初始地震子波构建衰减子波群，并对所述衰减子波群进行反射系数反演，获得反演后的反射系数；所述计算模块用于将所述反射系数代入L1-2范数函数后求二阶导数，获得求导曲线；所述分析模块用于最小化所述求导曲线中对应的L1-2范数函数的二阶导数，获得品质因子取值及对应的反射系数反演结果。7.根据权利要求6所述的非稳态盲反褶积装置，其特征在于，所述计算模块还包含预处理单元，所述预处理单元用于将反演后的反射系数代入L1-2范数函数获得扫描曲线；对所述扫描曲线求取二阶导数获得二阶曲线；将所述扫描曲线和所述二阶曲线分别进行光滑处理后再叠加处理，获得求导曲线。2CN113219530A权利要求书2/2页8.根据权利要求6所述的非稳态盲反褶积装置，其特征在于，所述反演模块包含：利用稀疏贝叶斯学习算法对所述衰减子波群进行反射系数反演，获得反演后的反射系数。9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。3CN113219530A说明书1/10页非稳态盲反褶积方法及装置技术领域[0001]本发明涉及地球物理勘探领域，尤指一种基于L1-2范数高阶导数准则的非稳态盲反褶积方法与装置。背景技术[0002]由于地层的非完全弹性性质，地震波在地层中传播时会产生振幅能量衰减、相位畸变、频带变窄的现象。品质因子Q是描述地层弹性性质