[发明专利]垂变网格地震波场外推方法及装置

权利要求书

1.一种垂变网格地震波场外推方法，其特征在于，所述方法包括：

根据预设的规则网格的速度模型，确定垂变网格的比例参数，并利用所述比例参数，通过坐标变换确定垂变网格的速度模型与规则网格的速度模型的坐标变换关系；

根据多元函数微分法则，利用所述坐标变换关系确定垂变网格与规则网格的偏导关系；

根据所述偏导关系进行有限差分运算，得到差分算子，利用所述差分算子对垂变网格进行地震波场数值模拟；

利用拉格朗日插值算法，对得到地震波场数值模拟结果进行插值恢复，恢复至规则网格。

2.根据权利要求1所述的方法，其特征在于，所述利用所述比例参数，通过坐标变换确定垂变网格的速度模型与规则网格的速度模型的坐标变换关系包括：根据坐标变换，利用拉格朗日插值算法，将比例参数插值至垂变网格，得到垂变网格的速度模型与规则网格的速度模型的坐标变换关系。

3.根据权利要求1所述的方法，其特征在于，所述根据多元函数微分法则，利用所述坐标变换关系确定垂变网格与规则网格的偏导关系包括：根据多元函数微分法则，利用所述坐标变换关系确定垂变网格与规则网格的偏导关系，并根据所述偏导关系，将规则网格的偏导变换成垂变网格的偏导。

4.根据权利要求1所述的方法，其特征在于，所述利用所述差分算子对垂变网格进行地震波场数值模拟包括：将所述差分算子离散化，利用有限差分方法对垂变网格进行地震波场数值模拟。

5.一种垂变网格地震波场外推装置，其特征在于，所述装置包括：

坐标变换模块，用于根据预设的规则网格的速度模型，确定垂变网格的比例参数，并利用所述比例参数，通过坐标变换确定垂变网格的速度模型与规则网格的速度模型的坐标变换关系；

偏导关系模块，用于根据多元函数微分法则，利用所述坐标变换关系确定垂变网格与规则网格的偏导关系；

数值模拟模块，用于根据所述偏导关系进行有限差分运算，得到差分算子，利用所述差分算子对垂变网格进行地震波场数值模拟；

插值恢复模块，用于利用拉格朗日插值算法，对得到地震波场数值模拟结果进行插值恢复，恢复至规则网格。

6.根据权利要求5所述的装置，其特征在于，所述坐标变换模块还具体用于根据坐标变换，利用拉格朗日插值算法，将比例参数插值至垂变网格，得到垂变网格的速度模型与规则网格的速度模型的坐标变换关系。

7.根据权利要求5所述的装置，其特征在于，所述偏导关系模块还具体用于根据多元函数微分法则，利用所述坐标变换关系确定垂变网格与规则网格的偏导关系，并根据所述偏导关系，将规则网格的偏导变换成垂变网格的偏导。

8.根据权利要求5所述的装置，其特征在于，所述数值模拟模块还具体用于将所述差分算子离散化，利用有限差分方法对垂变网格进行地震波场数值模拟。

9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一项所述方法的计算机程序。