深度域地震子波提取方法和装置、存储介质及电子设备.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN114200522A(43)申请公布日2022.03.18(21)申请号202010982672.2(22)申请日2020.09.17(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人李京南马灵伟胡华锋(74)专利代理机构北京聿宏知识产权代理有限公司11372代理人吴大建金淼(51)Int.Cl.G01V1/32(2006.01)G01V1/34(2006.01)G01V1/30(2006.01)权利要求书2页说明书14页附图4页(54)发明名称深度域地震子波提取方法和装置、存储介质及电子设备(57)摘要本发明涉及石油勘探开发技术领域，具体涉及一种深度域地震子波提取方法和装置、存储介质及电子设备，解决了现有技术中深度域地震子波提取方法未考虑到深度域地震数据的非稳态性、物理意义不明确、适用范围有限的问题；方法包括：获取目标地震道的深度域地震数据并进行波数谱分解，获得空间波数谱；将深度域地震数据划分为若干层段；对于每个层段，计算每个深度采样点处的波数振幅谱，并进行傅里叶逆变换，获得每个深度采样点处的深度域子波，计算所述层段的平均深度域子波；对若干层段的平均深度域子波进行融合，获得目标地震道的深度域地震子波；本发明考虑到了深度域地震数据的非稳态性，具有步骤更简便，适用范围更广的优点。CN114200522ACN114200522A权利要求书1/2页1.一种深度域地震子波提取方法，其特征在于，包括以下步骤：获取目标地震道的深度域地震数据，并对所述深度域地震数据进行波数谱分解，以获得所述深度域地震数据的空间波数谱，并确定所述空间波数谱的形态分布；根据所述空间波数谱的形态分布，将所述深度域地震数据划分为若干层段；对于每个层段，根据所述空间波数谱，计算每个深度采样点处的波数振幅谱，并对所述波数振幅谱进行傅里叶逆变换，以获得每个深度采样点处的深度域子波，计算所述层段内所有深度采样点处的深度域子波数值的平均值，以获得所述层段的平均深度域子波；对若干层段的平均深度域子波进行融合，以获得所述目标地震道的深度域地震子波。2.根据权利要求1所述的深度域地震子波提取方法，其特征在于，所述获取目标地震道的深度域地震数据，并对所述深度域地震数据进行波数谱分解，以获得所述深度域地震数据的空间波数谱，包括以下步骤：获取目标地震道的深度域地震数据；利用S变换对所述深度域地震数据进行波数谱分解，以获得所述深度域地震数据的空间波数谱。3.根据权利要求2所述的深度域地震子波提取方法，其特征在于，利用S变换对所述深度域地震数据进行波数谱分解，以获得所述深度域地震数据的空间波数谱时，采用的波数谱分解式为：其中，S(η,k)表示所述深度域地震数据的空间波数谱，η表示深度，k表示波数，u(x)表示目标地震道的深度域地震数据，x表示深度，e表示指数函数，i表示虚数单位，g(η‑x,k)表示深度域高斯窗函数，其表达式为：4.根据权利要求1所述的深度域地震子波提取方法，其特征在于，根据所述空间波数谱的形态分布，将所述深度域地震数据划分为若干层段，包括以下步骤：根据所述空间波数谱的形态分布，相邻波数变化不超过预设阈值的深度范围作为形态相似的空间波数谱，将形态相似的空间波数谱划分为一个层段，以将所述深度域地震数据划分为若干层段。5.根据权利要求1所述的深度域地震子波提取方法，其特征在于，根据所述空间波数谱，计算每个深度采样点处的波数振幅谱时，采用的波数振幅谱计算式为：其中，A(ηj,k)表示深度采样点处的波数振幅谱，ηj表示任意一个深度采样点，j＝0，1,2,...,N‑1，N表示深度域地震数据中深度采样点的个数，k表示波数，S(ηj,k)表示深度采样点ηj处的空间波数谱，其值为复数，real表示该复数的实部，imag表示该复数的虚部。6.根据权利要求1所述的深度域地震子波提取方法，其特征在于，对所述波数振幅谱进行傅里叶逆变换，以获得每个深度采样点处的深度域子波，采用的计算式为：2CN114200522A权利要求书2/2页其中，w(ηj,x)表示深度采样点处的深度域子波，A(ηj,k)表示深度采样点处的波数振幅谱，ηj表示任意一个深度采样点，j＝0，1,2,...,N‑1，N表示深度域地震数据中深度采样点的个数，x表示深度，k表示波数，e表示指数函数，i表示虚数单位。7.根据权利要求1所述的深度域地震子波提取方法，其特征在于，计算所述层段内所有深度采样点处的深度域子波数值的平均值，以获得所述层段的平均深度域子波，采用的计算式为：其中，表示第j层段的平均深度域子波，j＝0,1,2,...,M‑1，M表示所述深度域地震数