基于地震频率信息的多波匹配方法组成比例

提供一种基于地震频率信息的多波匹配方法，所述方法包括：(a)基于频率域的转换横波和频率域的纵波，求取频率域的转换横波的扩频比例系数；(b)基于求取的扩频比例系数对所述频率域的转换横波进行扩频处理；(c)对频率域的扩频处理后的转换横波进行变换，得到时间域的扩频处理后的转换横波；(d)将时间域的扩频处理后的转换横波与时间域的纵波进行匹配。上述方法使转换横波的频带宽度得到一定程度的拓展，转换横波的高频分量和分辨率得到明显提高，从而使得转换横波和纵波的匹配效果得到明显改善。

【技术实现步骤摘要】
基于地震频率信息的多波匹配方法
本专利技术总体来说涉及多波地震勘探
，更具体地讲，涉及一种基于地震频率信息的多波匹配方法。
技术介绍
随着经济的发展，对石油、天然气等资源的需求在日益增加，因此对石油、天然气的勘探要求也越来越高。面对越来越复杂的油气勘探开发需求，常规的纵波地震勘探技术存在诸多难以解决的问题，如真假亮点识别的多解性、反演出的岩石物性参数少、储层预测精度低等问题，因此上述单一的纵波勘探技术存在明显的不足。多波地震勘探投资小，获取的地下信息比单一纵波勘探技术更加丰富，是进行岩性油气藏和裂缝性油气藏勘探的一种有效手段。但目前多波勘探还存在一些尚未解决的问题，例如，如何校准纵波、转换横波剖面中的对应层位；褶积后如何处理两个地震道的波峰波谷不一致；如何获取横波速度、生成共转换点道集等。多波匹配是多波地震勘探技术的基础，目前国内外的多波匹配方法主流是基于同相轴相似性准则在时间域内对纵波(即，PP波)和转换横波(即，PS波)进行匹配，上述匹配方法的匹配精度较低。
技术实现思路
本专利技术的目的是提出一种基于地震频率信息的多波匹配方法，以解决现有中纵波和转换横波的匹配精度低的技术问题。本专利技术的一方面提供一种基于地震频率信息的多波匹配方法，所述方法包括：(a)基于频率域的转换横波和频率域的纵波，求取所述频率域的转换横波的扩频比例系数；(b)基于求取的扩频比例系数对所述频率域的转换横波进行扩频处理；(c)对频率域的扩频处理后的转换横波进行变换，得到时间域的扩频处理后的转换横波；(d)将时间域的扩频处理后的转换横波与时间域的纵波进行匹配。可选地，所述方法可还包括：获取多道地震记录，其中，步骤(a)可包括：基于每道地震记录对应的频率域的转换横波的振幅谱和频率域的纵波的振幅谱，求取每道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数，其中，步骤(a)中，求取所述多道地震记录中的任一道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数的步骤可包括：计算所述一道地震记录对应的频率域的转换横波振幅谱与所述一道地震记录对应的频率域的纵波振幅谱的比值，将所述比值作为所述一道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数。可选地，步骤(b)可包括：基于求取的每道地震记录对应的频率域的转换横波的扩频比例系数对每道地震记录对应的频率域的转换横波的振幅谱进行扩频处理，其中，步骤(b)中，对所述多道地震记录中的任一道地震记录对应的频率域的转换横波的振幅谱进行扩频处理的步骤包括：将所述一道地震记录对应的频率域的转换横波振幅谱乘以求取的所述一道地震记录对应的频率域的转换横波的扩频比例系数，以对所述一道地震记录对应的频率域的转换横波的振幅谱进行扩频处理，其中，步骤(c)包括：对多道地震记录对应的频率域的扩频处理后的转换横波的振幅谱进行变换，得到时间域的扩频处理后的转换横波。可选地，步骤(a)可包括：基于频率域的转换横波的点谱和频率域的纵波的点谱，求取频率域的转换横波的点谱的扩频比例系数。可选地，所述方法可还包括：获取多道地震记录；针对所述多道地震记录中的一道地震记录选取多个采样点，其中，步骤(a)可包括：基于每道地震记录中的每个采样点对应的频率域的转换横波的点谱和频率域的纵波的点谱，求取每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数，其中，步骤(a)中，求取任一道地震记录中的任一个采样点对应的频率域的转换横波的点谱的扩频比例系数的步骤可包括：计算所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱与所述一道地震记录中的所述一个采样点对应的频率域的纵波的点谱的比值，将所述比值作为所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱的扩频比例系数。可选地，计算任一道地震记录中的任一个采样点对应的频率域的转换横波的点谱的步骤可包括：对所述一道地震记录中的所述一个采样点对应的时间域的转换横波进行变换，得到所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱，其中，计算任一道地震记录中的任一个采样点对应的频率域的纵波的点谱的步骤可包括：对所述一道地震记录中的所述一个采样点对应的时间域的纵波进行变换，得到所述一道地震记录中的所述一个采样点对应的频率域的纵波的点谱。可选地，每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数形成扩频比例系数矩阵，其中，所述方法在步骤(a)和步骤(b)之间可还包括：对求取的扩频比例系数矩阵进行平滑处理。可选地，对求取的扩频比例系数矩阵中的任一扩频比例系数进行平滑处理的步骤可包括：以所述扩频比例系数为中心点，选取预定大小的移动窗口；针对所述移动窗口内包含的各扩频比例系数设定对应的权重值；基于所述移动窗口内包含的各扩频比例系数以及与所述各扩频比例系数对应的设定的权重值，计算所述移动窗口内包含的各扩频比例系数的均值，并将所述均值作为处于所述中心点的所述扩频比例系数对应的平滑后的扩频比例系数。可选地，当选取的预定大小的移动窗口超出所述扩频比例系数矩阵的边界时，超出所述边界部分的移动窗口内包含的参数的值可以0值填充。可选地，步骤(b)可包括：基于求取的每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数对每道地震记录中的每个采样点对应的频率域的转换横波的点谱进行扩频处理，其中，步骤(b)中，对任一道地震记录中的任一个采样点对应的频率域的转换横波的点谱进行扩频处理的步骤包括：将所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱乘以所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱的扩频比例系数，以对所述一道地震记录中的所述一个采样点对应的频率域的转换横波的点谱进行扩频处理，其中，在步骤(b)和步骤(c)之间还包括：求取每道地震记录对应的频率域的转换横波扩频处理后的振幅谱，其中，求取任一道地震记录对应的频率域的转换横波扩频处理后的振幅谱的步骤包括：分别计算任一道地震记录中的每个采样点对应的频率域的转换横波扩频处理后的点谱的绝对值；对所述一道地震记录中的所有采样点对应的频率域的转换横波扩频处理后的点谱的绝对值进行叠加，得到所述一道地震记录对应的频率域的转换横波扩频处理后的振幅谱，其中，步骤(c)包括：对多道地震记录对应的频率域的扩频处理后的转换横波的振幅谱进行变换，得到时间域的扩频处理后的转换横波。采用本专利技术所述的基于地震频率信息的多波匹配方法，以相同地震剖面的纵波频带宽度为参考，经过上述多波匹配方法的处理，使转换横波的频带宽度得到一定程度的拓展，转换横波的高频分量和分辨率得到明显提高，从而使得转换横波和纵波的匹配效果得到明显改善。附图说明通过下面结合附图进行的详细描述，本专利技术的上述和其它目的、特点和优点将会变得更加清楚，其中：图1是示出根据本专利技术的示例性实施例的基于地震频率信息的多波匹配方法的流程图；图2是示出根据本专利技术第一示例性实施例的求取频率域的转换横波的振幅谱的扩频比例系数的步骤的流程图；图3是示出根据本专利技术第二示例性实施例的求取频率域的转换横波的点谱的扩频比例系数的步骤的流程图；图4是示出根据本专利技术的示例性实施例的对扩频比例系数矩阵进行平滑处理的步骤的流程图；图5是示出根据本专利技术的示例性实施例的选取预本文档来自技高网...

【技术保护点】
一种基于地震频率信息的多波匹配方法，所述方法包括：(a)基于频率域的转换横波和频率域的纵波，求取所述频率域的转换横波的扩频比例系数；(b)基于求取的扩频比例系数对所述频率域的转换横波进行扩频处理；(c)对频率域的扩频处理后的转换横波进行变换，得到时间域的扩频处理后的转换横波；(d)将时间域的扩频处理后的转换横波与时间域的纵波进行匹配。

【技术特征摘要】
1.一种基于地震频率信息的多波匹配方法，所述方法包括：(a)基于频率域的转换横波和频率域的纵波，求取所述频率域的转换横波的扩频比例系数；(b)基于求取的扩频比例系数对所述频率域的转换横波进行扩频处理；(c)对频率域的扩频处理后的转换横波进行变换，得到时间域的扩频处理后的转换横波；(d)将时间域的扩频处理后的转换横波与时间域的纵波进行匹配。2.根据权利要求1所述的方法，其中，所述方法还包括：获取多道地震记录；其中，步骤(a)包括：基于每道地震记录对应的频率域的转换横波的振幅谱和频率域的纵波的振幅谱，求取每道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数，其中，步骤(a)中，求取每道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数的步骤包括：计算每道地震记录对应的频率域的转换横波振幅谱与每道地震记录对应的频率域的纵波振幅谱的比值，将所述比值作为每道地震记录对应的频率域的转换横波的振幅谱的扩频比例系数。3.根据权利要求2所述的方法，其中，步骤(b)包括：基于求取的每道地震记录对应的频率域的转换横波的扩频比例系数对每道地震记录对应的频率域的转换横波的振幅谱进行扩频处理，其中，步骤(b)中，对每道地震记录对应的频率域的转换横波的振幅谱进行扩频处理的步骤包括：将每道地震记录对应的频率域的转换横波振幅谱乘以求取的每道地震记录对应的频率域的转换横波的扩频比例系数，以对每道地震记录对应的频率域的转换横波的振幅谱进行扩频处理，其中，步骤(c)包括：对多道地震记录对应的频率域的扩频处理后的转换横波的振幅谱进行变换，得到时间域的扩频处理后的转换横波。4.根据权利要求1所述的方法，其中，步骤(a)包括：基于频率域的转换横波的点谱和频率域的纵波的点谱，求取频率域的转换横波的点谱的扩频比例系数。5.根据权利要求4所述的方法，所述方法还包括：获取多道地震记录；针对所述多道地震记录中的一道地震记录选取多个采样点，其中，步骤(a)包括：基于每道地震记录中的每个采样点对应的频率域的转换横波的点谱和频率域的纵波的点谱，求取每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数，其中，步骤(a)中，求取每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数的步骤包括：计算每道地震记录中的每个采样点对应的频率域的转换横波的点谱与每道地震记录中的每个采样点对应的频率域的纵波的点谱的比值，将所述比值作为每道地震记录中的每个采样点对应的频率域的转换横波的点谱的扩频比例系数。6.根据权利要求5所述的方法，其中，计算每道地震记录中...

【专利技术属性】
技术研发人员：黄东山，邹文，朱亚东，徐敏，刘璞，唐浩，彭才，张洞君，章雄，周晶晶，郑虹，罗晶，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51