[发明专利]基于角度域直接包络反演的盐丘速度密度估计方法

摘要

本发明公开了基于角度域直接包络反演的盐丘速度密度估计方法，具体为：S1、设置炮点和检波点，采集得到观测数据；S2、构建矩形网格地质模型；S3、给定直接包络反演阶段待优化的目标函数Je(v(x),ρ(x))；S4、给定用于速度长波长速度分量反演的散射角范围；S5、通过直接包络反演和迭代优化算法对目标函数Je(v(x),ρ(x))进行优化，得到速度模型的长波长分量vl(x)；S6、给定下一阶段新的待优化目标函数J(v(x),ρ(x))；S7、通过最速下降法对目标函数J(v(x),ρ(x))进行优化，得到精细的速度结构vf(x)；S8、给定用于密度反演的散射角范围；S9、通过最速下降法对目标函数进行进一步优化，得到精细的密度结构ρf(x)。本方法能够使用直接包络中新的Fréchet导数进行能量传播。

主权项

1.基于角度域直接包络反演的盐丘速度密度估计方法，其特征在于，具体按照以下步骤实施：步骤1、设置炮点和检波点，采集得到观测数据u_obs(t,x_r；x_s)，其中t表示时间变量；x_r,x_s分别表示检波点和炮点的位置；步骤2、构建矩形网格地质模型，根据技术指标及工区要求设定模型大小Nx×Nz，并设定正演模拟的空间采样间隔、时间采样间隔dt、最大采样时间Nt，所述空间采样间隔包括横向采样间隔dx和纵向采样间隔dz；步骤3、给定速度模型v(x)和密度模型ρ(x)的初始速度模型v₀(x)和初始密度模型ρ₀(x)，给定直接包络反演阶段待优化的目标函数J_e(v(x),ρ(x))；步骤4、给定用于速度长波长速度分量反演的散射角范围；步骤5、固定密度模型为初始密度模型ρ₀(x)，通过直接包络反演和迭代优化算法对目标函数J_e(v(x),ρ(x))进行优化，得到速度模型的长波长分量v_l(x)；步骤6、通过Gardner公式得到密度模型的长波长分量ρ_l(x)，并给定下一阶段新的待优化目标函数J(v(x),ρ(x))；步骤7、固定密度模型为长波长分量ρ_l(x)，通过最速下降法对目标函数J(v(x),ρ(x))进行优化，得到精细的速度结构v_f(x)；步骤8、给定用于密度反演的散射角范围；步骤9、固定速度模型v_f(x)，通过最速下降法对目标函数进行进一步优化，得到精细的密度结构ρ_f(x)。