粘声波全波形反演方法及装置

摘要

本发明公开一种粘声波全波形反演方法及装置，该方法包括：根据由多个Maxwell体构成的GSLS模型的粘弹性波动方程，确定时间域一阶速度-应力粘声波方程；根据时间域一阶速度-应力粘声波方程和二范数目标函数，确定该方程相应的松弛模量梯度计算公式及纵波速度梯度计算公式；根据时间域一阶速度-应力粘声波方程，采用高阶交错网格有限差分法确定正向传播波场及残差逆时反传波场；根据正向传播波场、残差逆时反传波场、松弛模量梯度计算公式及纵波速度梯度计算公式，确定松弛模量梯度和纵波速度梯度；根据纵波速度梯度采用局部最优化方法进行二维粘声波全波形反演。可以更好地近似地球介质在地震频带范围内的常Q特征，达到良好的数值频散抑制能力。

主权项

一种粘声波全波形反演方法，其特征在于，包括：根据GSLS模型的粘弹性波动方程，确定时间域一阶速度‑应力粘声波方程，其中，所述GSLS模型由多个麦克斯韦Maxwell体构成；根据所述时间域一阶速度‑应力粘声波方程和二范数目标函数，确定所述时间域一阶速度‑应力粘声波方程相应的松弛模量梯度计算公式及纵波速度梯度计算公式；根据所述时间域一阶速度‑应力粘声波方程，采用高阶交错网格有限差分法确定正向传播波场及残差逆时反传波场；根据所述正向传播波场、残差逆时反传波场、松弛模量梯度计算公式及纵波速度梯度计算公式，确定松弛模量梯度和纵波速度梯度；根据所述纵波速度梯度，采用局部最优化方法进行二维粘声波全波形反演。