一种突出主频的振幅谱平滑方法技术

本发明专利技术提供了一种突出主频的振幅谱平滑方法，属于油气、页岩油气及煤层气地震勘探与开发领域。本方法包括：(1)对振幅谱平均划分数据间隔，并提取极值；(2)对步骤(1)提取的数据进行线性插值得到已插值数据；(3)对步骤(2)得到的已插值数据进行平滑处理。本发明专利技术通过划分区域寻找主频，在肯定不会漏掉整个数据最大值的情况下，也保留了其他区间的最大值，并使得结果数据大致体现了原数据的整体趋势。

【技术实现步骤摘要】

本专利技术属于油气、页岩油气及煤层气地震勘探与开发领域，具体涉及一种突出主 频的振幅谱平滑方法。
技术介绍
在地震数据的处理与解释中时常会进行地震数据的时频转换，计算其频谱（时间 域的信号在频率域下的表示方式）。所得的结果会是以分别以振幅及相位为纵轴，频率为横 轴的两张图。振幅谱有时也称"振幅频谱"，表示振幅随频率变化的情形。 数据平滑是有效去除数据噪音以提高信噪比的技术手段。数据平滑的技术手段很 多，如移动窗口平均法、移动窗口拟合多项式平滑、回归平滑等等。 地震资料的主频含义是指，已转换到频率域的信号中其极值主要集中的频段范 围。 目前并没有突出主频的振幅谱平滑的相关技术。
技术实现思路
本专利技术的目的在于解决上述现有技术中存在的难题，提供一种突出主频的振幅谱 平滑方法，利用一定的技术手段对振幅谱进行平滑处理、得到突出体现振幅谱的主频信息 平滑数据。利用此类数据绘制的振幅图有助于清晰的显示出相应地震数据的主频范围，并 有一定的美观效果，为常规油气、页岩油气及煤层气的勘探、开发服务。 本专利技术是通过以下技术方案实现的： -种突出主频的振幅谱平滑方法，包括： (1)对振幅谱平均划分数据间隔，并提取极值； (2)对步骤（1)提取的数据进行线性插值得到已插值数据； (3)对步骤（2)得到的已插值数据进行平滑处理。 所述步骤（1)是这样实现的： 首先统计振幅谱的样点数据及个数，然后给定数据划分间隔值，把数据划分为若 干个区间，提取每个区间的样点数据的最大值，并记录原样点坐标。 所述间隔值为把振幅谱的样点数据平均划分为5~10个区间的值。 所述步骤（2)是这样实现的： 利用线性插值法对步骤（1)提取的数据进行插值处理，插值的个数为原振幅谱的 样点个数。 与现有技术相比，本专利技术的有益效果是：本专利技术通过划分区域寻找主频，在肯定不 会漏掉整个数据最大值的情况下，也保留了其他区间的最大值，并使得结果数据大致体现 了原数据的整体趋势。【附图说明】 图1是原振幅谱样点数据（共25个样点，等分为5区间，每个区间样点为5)。 图2是每个区间都提取一个最大值，并记住其在原样点数据中的位置。 图3是已知坐标（x0, y0)与（xl，yl)，要得到区间内某一位置X在直线 上的值的示意图。 图4是对已提取的数据进行线性插值（插值后的样点个数为原振幅谱的样点个 数） 图5是对已插值的数据进行平滑处理。【具体实施方式】 下面结合附图对本专利技术作进一步详细描述： 本专利技术主要分为以下三个步骤： (1)平均划分数据间隔、并提取极值： 首先统计振幅谱的样点数据及个数，然后给定一个数据划分间隔值（以把振幅谱 数据平均划分为5~10个区间的间隔值为适宜）。把数据划分为若干个区间，提取每个区 间的样点数据的最大值，并记录原样点坐标，如图1、图2所示。振幅谱数据体现的就是频率 信息，每一个区间提取一个最大值，把全部区间得到的最大值做个统计，那么其中最大的一 个值，可以理解为这个数据的主频。 (2)对已提取数据进行线性插值： 对第一步提取的数据进行插值处理，插值的个数为原振幅谱样点个数。利用简单 的线性插值即可。线性插值法是指使用连接两个已知量的直线来确定在这两个已知量之间 的一个未知量的值的方法。 假设已知坐标（x0, y0)与（xl，yl)，要得到区间内某一位置X在直线上 的值。枏据图3中所示,得至1丨 假设方程两边的值为α，那么这个值就是插值系数-从x〇到X的距离与从x〇到 Xl距离的比值。由于X值已知，所以可以从公式得到α的值 这样，在代数上就可以表示成为： y = (I- a ) y〇+ a yl 或者， y = y0+ a (yl-y〇) 这样通过（！就可以直接得到y。实际上，即使x不在x0到xl之间并且α也不是 介于0到1之间,这个公式也是成立的。在这种情况下,这种方法叫作线性外插。已知y求 X的过程与以上过程相同，只是X与y要进行交换。如图4所示。 (3)对已插值数据进行平滑处理： 对第二步得到的插值后的数据进行平滑处理。利用常规平滑算法即可，例如可用 五点三次平滑，如下： 设已知η个等距点X。< X1 <~< xn 1上的观测数据（或实验数据）为…， yn 1，则可以在每个数据点的前后各取两个相邻的点，用三次多项式y = ao+aj+aj^+a；^3进 行逼近。 根据最小二乘原理确定出系数a。，ap a2, a3,最后可得到五点三次平滑公式如下： 具中:1?表不Y1的平滑值，如图5所不。 本专利技术涉及的方法技术是一项适用于常规油气的勘探开发和非常规（煤层气和 页岩油气）油气的勘探开发，其方法技术的关键点包括将振幅谱样点按一定的间隔提取区 间最大值并记录原样点坐标，然后对此组样点进行插值处理（插满原始样点个数），最后进 行平滑处理。利用此类数据绘制的振幅图有助于清晰的显示出相应地震数据的主频范围， 并有一定的美观效果，为常规油气、页岩油气及煤层气的勘探、开发服务。 上述技术方案只是本专利技术的一种实施方式，对于本领域内的技术人员而言，在本 专利技术公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 专利技术上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。【主权项】1. ，其特征在于：所述方法包括： (1) 对振幅谱平均划分数据间隔，并提取极值； (2) 对步骤（1)提取的数据进行线性插值得到已插值数据； (3) 对步骤（2)得到的已插值数据进行平滑处理。2. 根据权利要求1所述的突出主频的振幅谱平滑方法，其特征在于：所述步骤（1)是 这样实现的： 首先统计振幅谱的样点数据及个数，然后给定数据划分间隔值，把数据划分为若干个 区间，提取每个区间的样点数据的最大值，并记录原样点坐标。3. 根据权利要求2所述的突出主频的振幅谱平滑方法，其特征在于：所述间隔值为把 振幅谱的样点数据平均划分为5~10个区间的值。4. 根据权利要求1所述的突出主频的振幅谱平滑方法，其特征在于：所述步骤（2)是 这样实现的： 利用线性插值法对步骤（1)提取的数据进行插值处理，插值的个数为原振幅谱的样点 个数。【专利摘要】本专利技术提供了，属于油气、页岩油气及煤层气地震勘探与开发领域。本方法包括：(1)对振幅谱平均划分数据间隔，并提取极值；(2)对步骤(1)提取的数据进行线性插值得到已插值数据；(3)对步骤(2)得到的已插值数据进行平滑处理。本专利技术通过划分区域寻找主频，在肯定不会漏掉整个数据最大值的情况下，也保留了其他区间的最大值，并使得结果数据大致体现了原数据的整体趋势。【IPC分类】G01V1/28【公开号】CN105093284【申请号】CN201410208297【专利技术人】宋志翔, 刘百红, 岳承琪, 杨强 【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院【公开日】2015年11月25日【申请日】2014年5月16日本文档来自技高网...

【技术保护点】
一种突出主频的振幅谱平滑方法，其特征在于：所述方法包括：(1)对振幅谱平均划分数据间隔，并提取极值；(2)对步骤(1)提取的数据进行线性插值得到已插值数据；(3)对步骤(2)得到的已插值数据进行平滑处理。

【技术特征摘要】

【专利技术属性】
技术研发人员：宋志翔，刘百红，岳承琪，杨强，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11