[发明专利]地震矢量波场数值模拟方法和装置

摘要

本发明提供了一种地震矢量波场数值模拟方法和装置，其中，该方法包括对原始窗函数做L次自褶积运算得到自褶积后的窗函数；对自褶积后的窗函数与原始窗函数进行加权组合运算，得到自褶积组合窗函数；通过自褶积组合窗函数截断伪谱法的空间褶积序列得到有限差分算子；根据所述有限差分算子对地震矢量波场进行数值模拟。本发明解决了现有技术中，窗函数截断逼近的有限差分算子无法兼顾谱覆盖范围和精度误差稳定性的问题，从而导致的有限差分地震矢量波场数值模拟结果数值频散较大，数值模拟精度不高的技术问题，达到了有效提高地震矢量波场数值模拟精度的技术效果。

主权项

一种地震矢量波场数值模拟方法，其特征在于，包括：对原始窗函数做L次自褶积运算得到自褶积后的窗函数，其中，L为正整数；对自褶积后的窗函数与原始窗函数进行加权运算，得到自褶积组合窗函数；通过自褶积组合窗函数截断逼近得到有限差分算子；根据所述有限差分算子对地震矢量波场进行数值模拟；其中，所述对自褶积后的窗函数与原始窗函数进行加权运算，得到自褶积组合窗函数，包括：在确定需要主瓣性能优先的情况下，原始窗函数的加权系数高于自褶积后的窗函数的加权系数；在确定需要旁瓣性能优先的情况下，自褶积后的函数的加权系数高于原始窗函数的加权系数；所述通过自褶积组合窗函数截断逼近得到有限差分算子，包括：通过得到的自褶积组合窗函数对伪谱法空间褶积序列进行截断处理，得到有限差分算子；根据得到的有限差分算子绘制逼近误差曲线；确定所述逼近误差曲线的频谱覆盖范围和逼近精度稳定性；如果频谱覆盖范围和逼近精度稳定性不满足预定要求，则改变L的值或者是改变加权运算的加权系数重新确定自褶积组合窗函数，直至确定的自褶积组合窗函数截断逼近得到的有限差分算子的频谱覆盖范围和逼近精度稳定性满足预定要求为止。