改进直接包络反演与扰动分解的强散射介质反演方法技术

本发明专利技术涉及一种改进直接包络反演与扰动分解的强散射介质反演方法，包络运算为非线性运算，在波场重叠处计算包络场存在较大误差。本发明专利技术是在计算正传包络场之前对地震波场进行了上下左右行波分解，有效缓解了波场重叠对正传包络场计算的影响。将改进的包络场代入直接包络反演算法中即为改进的直接包络反演，其能对强散射体速度结构进行更好的刻画。针对强散射屏蔽区速度建模，本发明专利技术提出采用构造约束的扰动分解方法。利用逆时偏移成像结果约束扰动分解参数的选取能避免参数选取导致对强散射体的错误识别。将改进的直接包络反演方法与构造约束扰动分解方法联合，既能得到好的强散射体速度建模，也能得到较好的屏蔽区速度构造。

An improved inversion method of strong scattering medium based on direct envelope inversion and disturbance decomposition

【技术实现步骤摘要】
改进直接包络反演与扰动分解的强散射介质反演方法
：本专利技术涉及一种利用改进直接包络反演方法获得强散射体大尺度结构、利用构造约束扰动分解方法得到强散射体下边界与屏蔽区速度结构的强散射介质反演方法。
技术介绍
：目前，在地震勘探领域，速度建模精度最高的方法是全波形反演方法。全波形反演方法将速度反演问题看作以数据残差为目标函数的最优化问题，往往采用局部优化算法求解。常规全波形反演之所以可以采用局部优化算法求解，是其采用了玻恩近似这一弱散射近似。这导致常规全波形反演方法只能用来处理弱散射问题的反演。近年来，随着盐丘体构造中油气藏的发现，工业界对盐丘体的速度建模问题十分关注。盐丘体与围岩往往具有较大速度反差，是典型的强散射体。常规全波形反演方法进行这类问题的速度建模是十分困难的，除非依赖非常低频的地震数据。而低频波场的激发和接收均较为困难，因此，目前认为在常规地震采集中，低频成分往往是不充分的。因此，研究低频缺失情况下强散射介质的速度建模具有重要的理论意义和实际意义。直接包络反演方法是一种能够在低频地震信息缺失情况下进行强散射体速度建模本文档来自技高网...

【技术保护点】
1.一种改进直接包络反演与扰动分解的强散射介质反演方法，其特征在于，是在地震数据缺失低频信息且没有模型先验信息的情况下，通过基于波场分解的改进直接包络反演得到强散射体的大尺度构造，通过基于构造约束的扰动分解方法得到强散射体的下边界以及下方屏蔽区速度构造的方法；/n包括以下步骤：/na、安装MATLAB软件平台；/nb、对数据进行静校正、去噪预处理，得到高质量的观测地震数据；/nc、对地震数据进行子波估计，提取每一炮数据的震源子波；/nd、通过背景速度分析，得到大致速度范围的最大和最小值，生成一个在速度范围内，随深度逐渐递增的线性梯度模型，其中不含有强散射体的任何先验信息，作为反演的初始模型v

【技术特征摘要】
1.一种改进直接包络反演与扰动分解的强散射介质反演方法，其特征在于，是在地震数据缺失低频信息且没有模型先验信息的情况下，通过基于波场分解的改进直接包络反演得到强散射体的大尺度构造，通过基于构造约束的扰动分解方法得到强散射体的下边界以及下方屏蔽区速度构造的方法；
包括以下步骤：
a、安装MATLAB软件平台；
b、对数据进行静校正、去噪预处理，得到高质量的观测地震数据；
c、对地震数据进行子波估计，提取每一炮数据的震源子波；
d、通过背景速度分析，得到大致速度范围的最大和最小值，生成一个在速度范围内，随深度逐渐递增的线性梯度模型，其中不含有强散射体的任何先验信息，作为反演的初始模型v0；
e、在初始模型上计算模拟地震数据，对模拟数据取包络得到模拟包络数据esyn，对观测地震数据取包络得到观测包络数据eobs，通过公式(1)计算改进直接包络反演的目标函数：



式中，求和符号下角标sr表示对所有震源和检波点进行积分，T表示总的记录时间长度；
f、在初始模型上计算模拟地震波场Pf，存储所有时间点的波场快照，对模拟波场进行上下行波分解，分别得到上行波Pfup和下行波Pfdown，再分别对上行波和下行波进行左右行波分解，得到上左行波Pfupleft、上右行波Pfupright、下左行波Pfdownleft和下右行波Pfdo...

【专利技术属性】
技术研发人员：张盼，韩立国，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22