[发明专利]地震矢量波场有限差分数值模拟方法、设备以及系统

说明书

本发明提供了一种地震矢量波场有限差分数值模拟方法、系统、计算机设备以及计算机可读存储介质，涉及地震勘探技术领域。该方法包括：根据多种不同参数的窗函数以及余弦组合窗函数确定混合褶积优化窗函数；根据所述混合褶积优化窗函数对伪谱法的空间褶积序列进行截断，得到优化的有限差分算子；根据所述优化的有限差分算子对地震矢量波场进行数值模拟。本发明构建一种基于最小二乘组合的混合褶积优化窗函数，利用该窗函数去截断伪谱法的空间褶积序列，得到优化的有限差分算子，再利用该有限差分算子进行地震矢量波场数值模拟，提高地震矢量波场数值模拟的精度与效率。

技术领域

本发明关于地震勘探技术领域，特别是关于地震领域的数值模拟技术，具体的讲是一种地震矢量波场有限差分数值模拟方法、系统、计算机设备以及计算机可读存储介质。

背景技术

地震波场数值模拟是基于波动理论，通过数值计算，模拟地震波在地球介质中的传播过程。波场数值模拟是地球物理反演以及成像技术的基础，可以用地震波场数值模拟的手段在弹性参数和地震响应之间建立起一个联系，通过波场数值模拟的方式知道在已知的一个储层条件下能得到一个什么样的地震响应。相反，假设有了一个地震响应，可以应用反演算法来得到弹性参数，有了这些弹性参数，可以对反演结果做一个基于岩石物理模型的解释，得到相应的储层参数。微分方程法基于连续介质微分体积元弹性动力学原理，常用的微分方程法包括有限差分和有限元法，是地震波场数值模拟的一个重要工具和方法。因此，有了精度和效率更高的有限差分算法，就可以获得精度更高的地震波场数值模拟结果，更好的服务于成像和反演的需要。

有限差分算法就是用差分算子近似微分算子，这种近似势必会导致引入数值频散，也就是误差。针对数值频散，现行主流有两种方法去优化有限差分算法，提高算法精度，其一是最优化方法，其二是窗函数法。两种方法的本质都是以期获得最大谱覆盖范围的前提下，降低误差限，只是最优化方法是通过搜索最优解提高算法精度，窗函数法是通过交互设计最优解，以提高算法精度。近年来，国内外地球物理学界针对有限差分算法的优化以提高地震波场数值模拟的精度，做了许多研究工作，取得了较好的效果。

应用最优化方法计算得到优化的有限差分系数，实际上是一个多参数优化的问题，目前被广泛使用优化有限差分算子的最优化方法有最小二乘法，模拟退火法以及Remez算法等，针对不同的优化方法，需要使用不同的目标函数。通过诸如最小二乘法、Remez算法去搜索解空间内的最优解，以满足误差限的要求，通常面临收敛慢，收敛不到全局最优解或者不收敛的问题。优选窗函数法采用不同的窗函数去截断伪谱法的空间褶积序列推导出优化的有限差分算子。

因此，如何提供一种新的有限差分算法，将其应用于地震矢量波场数值模拟，以实现地震矢量波场的高精度高效率的数值模拟是本领域亟待解决的技术难题。

发明内容

有鉴于此，本发明实施例提供了一种地震矢量波场有限差分数值模拟方法、系统、计算机设备以及计算机可读存储介质，首先对多种窗函数进行混合褶积，其次是对混合褶积之后的多个混合褶积窗函数以及余弦组合窗函数进行基于最小二乘方法的优化组合，得到一种基于最小二乘组合的混合褶积优化窗函数，最后利用该窗函数去截断伪谱法的空间褶积序列，得到优化的有限差分算子，再利用该有限差分算子进行地震矢量波场数值模拟，提高了地震矢量波场数值模拟的精度与效率。

本发明的目的之一是，提供一种地震矢量波场有限差分数值模拟方法，包括：

根据多种不同参数的窗函数以及余弦组合窗确定混合褶积优化窗函数；

根据所述混合褶积优化窗函数截断伪谱法的空间褶积序列，得到优化的有限差分算子；

根据所述优化的有限差分算子对地震矢量波场进行数值模拟。

优选的，所述根据多种不同参数的窗函数以及余弦组合窗确定混合褶积优化窗函数包括：

选取多种窗函数；