一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

三维生物礁滩建模及地震数值模拟     

黄继伟  熊晓军  贺振华 《物探化探计算技术》2014,(3):326-330

基于塔礁、点礁、生物滩等地质体的几何特征及地质特征,采用高精度的三维块状建模方法,建立了适合于复杂地区的三维礁滩复合体模型。该模型不但立体地反映了生物礁、生物滩及地层间的空间展布特征,还可以直观地体现了生物礁滩储层内部的物性差异。此外采用FK域三维波动方程数值模拟方法对其进行了地震波场数值模拟,得到了一系列具有代表性时间剖面上的地震波场响应特征,为正确认识实际生物礁滩的复杂结构和特性提供了依据。  相似文献