一种高精度的时变子波反演方法技术

本发明专利技术公开了一种高精度的时变子波反演方法，包括：一、取地震信号二阶累积量自相关作为反演初始子波的初始值；二、在地震信号时间范围内，选择至少大于四倍子波长度的子波反演时窗长度进行时窗划分；三、在每个时窗范围内，利用时窗内的井反射系数构建反射系数矩阵；四、以步骤一的反演初始子波为反演迭代初始值，采用误差下降迭代算法进行子波反演；五、按照步骤三、四顺序遍历所有时窗，完成时变子波的高精度反演。本发明专利技术对子波相位无需任何假设，依据褶积的交换率，创建反射系数矩阵，以正演记录与实际地震的最佳逼近为目标函数，在每个时窗内可反演任意相位子波，最终得到时变子波。

【技术实现步骤摘要】
一种高精度的时变子波反演方法
本专利技术涉及油气勘探
，尤其涉及一种高精度的时变子波反演方法。
技术介绍
地震子波的提取是地震高分辨率反演的重要环节，目前子波提取主要分为井震结合子波提取和地震高阶累积统计子波提取两类方法。其中井震结合子波提取方法，用地震复赛谱与井上反射系数复赛谱相减得到子波复赛谱，对子波复赛谱采用多项式谱模拟方式模拟出子波的振幅谱，而相位谱通常是人为设定；地震高阶累积量统计子波是一种全局统计子波，通常二阶统计子波只能统计出零相位子波，三阶和四阶统计量的复倒谱子波估计可以同时估算出子波的振幅谱与相位谱，子波估计的过程中需要估算反演各个相位谐波子波的系数，利用谐波子波构建混合相位统计子波，但算法对地震长度敏感，算法需要求解以谐波子波为目标的超定病态方程组，目标方程求解困难。受地层非均匀性和粘弹性的影响，实际地震子波具有时变的特征，传统的地震子波提取大多是基于时不变的假设，往往不能满足实际勘探所需的目标要求，传统的高阶统计量子波提取方法计算需要的时间范围较大，很难保证大范围地震数据的平稳性；而谱模拟方法是假定模拟类似于本文档来自技高网...

【技术保护点】
1.一种高精度的时变子波反演方法，其特征在于，包括以下步骤：/n步骤一：取地震信号二阶累积量自相关作为反演初始子波的初始值；/n步骤二：在地震信号时间范围内，选择至少大于四倍子波长度的子波反演时窗长度进行时窗划分；/n步骤三：在每个时窗范围内，利用时窗内的井反射系数构建反射系数矩阵；/n步骤四：以步骤一的反演初始子波为反演迭代初始值，采用误差下降迭代算法进行子波反演；/n步骤五：按照步骤三、四顺序遍历所有时窗，完成时变子波的高精度反演。/n

【技术特征摘要】
1.一种高精度的时变子波反演方法，其特征在于，包括以下步骤：
步骤一：取地震信号二阶累积量自相关作为反演初始子波的初始值；
步骤二：在地震信号时间范围内，选择至少大于四倍子波长度的子波反演时窗长度进行时窗划分；
步骤三：在每个时窗范围内，利用时窗内的井反射系数构建反射系数矩阵；
步骤四：以步骤一的反演初始子波为反演迭代初始值，采用误差下降迭代算法进行子波反演；
步骤五：按照步骤三、四顺序遍历所有时窗，完成时变子波的高精度反演。


2.根据权利要求1所述的一种高精度的时变子波反演方法，其特征在于，所述步骤一中的地震信号二阶累积量自相关的计算表达式为：
RR(s,t)＝E[x(s)x(t)]
式中，RR(s,t)是二阶累积量自相关，E是求和运算，x是地震振幅，s、t是振幅的序号，代表两个不同的时间点。


3.根据权利要求1所述的一种高精度的时变子波反演方法，其特征在于，所述步骤三中的反射系数矩阵G表达式如下式所示：



式中：m是子波的长度，n是时窗的长度，R是反射系数。


4.根据权利要求1...

【专利技术属性】
技术研发人员：宋明水，何文渊，毕建军，曹佳佳，
申请(专利权)人：中国石油化工股份有限公司胜利油田分公司，中国石油集团工程咨询有限责任公司，北京中恒利华石油技术研究所，
类型：发明
国别省市：山东;37