[发明专利]弹性介质主动源和被动源混采地震数据全波形反演方法

权利要求书

1.一种弹性介质主动源和被动源混采地震数据全波形反演方法，其特征在于，包括以下步骤：

a、安装MATLAB软件平台，要求采用MATLAB R2016a及以上版本，并且已配备并行工具包；

b、对被动源地震数据进行预处理，主要包括能量归一化和去噪处理，得到高质量的多分量被动源地震数据；

c、应用基于互相关的被动源地震干涉法，重构不同模式的虚拟炮记录，重构公式为：

其中，x_i和x_j代表两个检波器位置，上角标p代表检波器分量，上角标q代表震源分量，对于二维情况，均代表x分量或z分量，G^p,q(x_j,x_i,t)代表在x_i位置的q分量震源激发，在x_j位置接收的p分量格林函数，t代表时间，*表示褶积运算，S(t)表示虚拟震源子波，ρ代表密度，v_p表示纵波速度，d^p(x_j,t)表示在x_j接收的p分量地震记录，d^q(x_i,t)表示在x_i接收的q分量地震记录，<g>表示空间整体平均，公式(1)简写为：

其中，代表重构的虚拟炮记录，i和j分别表示虚拟震源和虚拟检波器的位置，q和p分别表示震源和检波器的分量；

重构后得到四种类型的虚拟炮记录，分别是和

d、构建初始速度模型，采用估计的最大速度和最小速度进行线性插值，得到由浅到深递增的梯度模型作为初始速度模型；

e、构建不依赖子波的弹性波多分量全波形反演的目标函数如下：

其中，u表示合成地震记录，上角标p代表检波器分量，上角标q代表震源分量，k表示参考道索引；

f、采用z分量虚拟震源记录和作为多分量观测记录，选取低频子波作为模拟子波，进行不依赖子波的弹性波多分量全波形反演，更新迭代后，得到一组纵波速度模型Vp和横波速度模型Vs；

g、采用x分量虚拟震源记录和作为多分量观测记录，选取低频子波作为模拟子波，进行不依赖子波的弹性波多分量全波形反演，更新迭代后，得到一组纵波速度模型Vp和横波速度模型Vs；

h、对主动源地震记录进行预处理，主要包括静校正、去噪、去多次波，得到高质量的主动源地震数据；

i、根据主动源地震数据估计震源子波；

j、以第f步的Vp模型和第g步的Vs模型作为初始模型，进行主动源弹性波全波形反演的更新迭代，达到迭代停止条件后，输出最终反演结果，即纵波速度模型和横波速度模型。