地震波场模拟中的褶积微分算子法

褶积微分算子是将偏导数算子的带限傅氏谱进行反傅氏变换得到的,具有对称或反对称和衰减速度快的特点。根据频谱分析理论,选用序列在结束时能光滑趋于零的截断窗口将褶积算子截断,结合考虑微分算子守恒定律,将褶积算子中心值进行修正,得到适合于波场模拟的最优权系数。给出用褶积微分算子模拟声波波场的实现途径与实例。理论分析与模型试算证明了褶积微分算子法的正确性与有效性。     

地震波传播的褶积微分算子法数值模拟

为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息.      

地震波场的高阶交错网格有限差分模拟

应用高阶交错网格有限差分算法,并加入吸收边界条件和衰减带,对弹性波方程进行模拟,分析了其稳定性和收敛性。通过对各向同性和各向异性介质模型的模拟表明,高阶差分波动方程模拟的网格频散较小、精度较高、效果较好,可为地震勘探及其资料解译提供技术手段。     

地震波混合阶褶积算法模拟

基于正反傅立叶变换,提出了地震波模拟的混合阶褶积算法。该方法原理简单、易于实现,结合了有限差分与伪谱法的优点,具有较高的精度和计算效率,适用于地震波场的正演计算,同时还给出了二维地震波场的理论计算实例。计算结果表明,此算法模拟结果正确、精度高、速度快、能适应较为复杂地质模型,并且易于推广到各向异性介质中去。     

褶积微分算子法数值模拟的边界问题研究

褶积微分算子法是一种全新的数值模拟方法,已被广泛应用于复杂介质的地震波场数值模拟,但是其边界反射问题一直没有解决。这里将最佳匹配层(PML)吸收边界条件引入到褶积微分算子法中,此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波,随传播距离按指数规律衰减,不产生任何反射,以达到消除边界反射的目的。构造不同的模型,通过对比分析证明,PML吸收边界条件能比较好地解决褶积微分算子法的边界问题,从而验证了完全匹配层吸收效果的优越性。     

交错网格地震波场模拟及频散校正策略

交错网格高阶差分解法是地震波场模拟的一种有效方法。较大的时间延拓步长可以提高效率,但可能引起算法不稳定,使模拟无法进行。因此,对于空间网格大小和时间步长的选择要适中。虽然震源子波的主频对模拟精度的影响不大,但适当提高震源子波的主频,可以提高地震波对于薄层的分辨率。FCT方法能够有效地压制数值频散,将FCT方法与交错网格有限差分法相结合,可以提高波场模拟的精度和运算效率。     

煤储层各向异性地震波场模拟

对中国煤层的裂隙系统进行了系统分类;阐述了各类裂隙的主要成因。基于二维二分量各向异性弹性波方程交错网格高阶有限差分法数值模拟,对各向同性介质模型和高裂隙密度(含干裂隙和饱和水裂隙)煤层的二维层状模型进行了弹性波场数值模拟。结果表明:煤层含裂隙后与不含裂隙时相比,反射时差的变化已不足以识别,但由于波速的改变会使煤层顶底、界面的波阻抗发生变化,从而引起波的动力学特征的改变,AVO现象尤其明显。含气和含水裂隙煤层中在x分量和z分量表现为极性相反,而且含水裂隙煤层中的波速比含气裂隙煤层的波速要快。因此,应用各向异性理论研究煤层裂隙的特征是可能的。     

二阶弹性波动方程高精度交错网格波场分离数值模拟

给出了一种等价的二阶弹性波动方程,以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波波场问题.应用高阶交错网格有限差分法求解该波动方程,并使用通量校正技术(FCT)进一步压制频散,采用均匀介质模型和层状介质模型进行波场分离数值试验,精确得到了混合波场、完全分离的纯纵波及纯横波波场.数值结果分析表明,本文方法在均匀介质情况下准确可靠,在分离后的纯纵、横波波场中可观察到较为丰富的能量转换信息,这对认识复杂弹性波的传播规律及弹性波理论具有重要意义.     

碳酸盐岩礁相地层地震波场特征分析

碳酸盐岩礁作为特殊的碳酸盐岩隆,是石油和天然气良好的储存场所。为了揭示其地震波场响应特征,采用高精度相移加插值波场延拓算子,对根据某地区生物礁相地层的地质特征所设计的模型,进行地震波场正演和偏移处理。结果表明,生物礁的特殊地质条件决定了其在地震剖面上必然出现异常响应,对响应特征的分析解释能够在油气勘探中发挥作用,所做的地震数值模拟工作,对碳酸盐岩礁体油气藏勘探具有指导性意义。     

利用显示微分算子对松散网格上2D波传播的数值模拟

我们提供了一种在松散2D网格上模拟地震波传播的数值方法。该算法以弹性波方程的速度——应力公式为基础，因此利用一种交错排列网格法。不像有限单元或频谱单元法，这种方法也能处理柔性松散网格，我们利用显函数微分算子计算每一时间步的空间导数。