排序方式： 共有4条查询结果，搜索用时 0 毫秒
1
1.
三维波动方程时空域混合网格有限差分数值模拟方法 总被引：1，自引：0，他引：1
常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推导出了两类非坐标轴网格点构建三维Laplace差分算子的方法,进而提出了一种利用坐标轴网格点和非坐标轴网格点共同构建三维Laplace差分算子的混合网格有限差分方法,并利用时空域频散关系和泰勒展开建立差分系数方程,推导出了差分系数的通解.相比常规高阶和时空域高阶差分格式的2阶差分精度,时空域混合网格差分离散波动方程理论上能够达到任意偶数阶差分精度,模拟精度显著提高,同时稳定性更强.频散分析表明:相比常规高阶和时空域高阶差分格式,在计算效率基本相同时,时空域混合网格差分格式能更有效地减小数值频散,减弱数值各向异性,模拟精度更高;在模拟精度基本相当时,混合网格差分格式能采用更大的时间采样间隔,计算效率更高.数值模拟实例进一步验证了混合网格差分格式在提高模拟精度和计算效率方面的先进性,也验证了其普遍适用性. 相似文献
2.
常规虚源点Marchenko自聚焦多次波预测方法只适用于预测不含自由表面的多次波模型,局限于压制层间多次波,该方法在构建上下行格林函数场前,必须从反射响应中去除所有与表面相关的多次波.本文对构建上下行Marchenko格林函数方程进行改进,得到了包含一次波、层间多次波和自由表面多次波的格林函数,利用改进的Marchenko自聚焦预测方法预测自由表面多次波.本文利用水平层状模型数据及SMARRT模型数据证明,改进后的Marchenko法预测海底相关的自由表面多次波效果较为理想,该方法避免了常规SRME自由表面多次波预测方法需要近道重构的缺陷,能够有效提高地震资料的信噪比和分辨率. 相似文献
3.
真振幅成像是一种代表性的定量估计模型参数扰动高波数部分的地震波成像方法.经典的真振幅成像方法在高频近似和理想照明假设条件下求取显式对角Hessian逆矩阵作为偏移振幅加权算子,用以校正波传播过程中的几何扩散效应,得到模型参数扰动的带限估计.真振幅保真成像方法在利用逆时偏移（RTM）框架实现时会产生低波数噪声,影响对高波数参数估计的精度.本文给出了一种新的基于RTM框架的真振幅保真成像条件,该成像条件针对反射波数据,在高频近似下散射模式对应正问题及Bayes反问题框架下导出.与传统基于高频渐进反演的波动方程成像方法类似,利用本文提出RTM成像条件能够保证计算结果与高频近似下反演结果的一致性.同时,利用本文提出RTM真振幅成像条件能够在成像过程中自动保真的消除传统真振幅RTM算法中存在低波数噪声,模型数值实验结果验证了本文方法的正确性和有效性. 相似文献
4.
各向异性研究对地下介质精确成像有着重要的意义,在当前计算机硬件迅速发展及宽方位地震数据采集日益普遍的情况下,成像必须考虑介质的各向异性.逆时偏移是基于双程波动方程的较为精确的数值解的成像方法,所以相对于其他地震成像方法,它具有很大的优势,譬如不受反射界面的倾角限制、偏移速度结构合适时能够使回转波及多次波正确成像.在各向同性介质中,可使用标量波方程来模拟波场.而在各向异性介质中,P波和SV波是相互耦合的,即不存在单纯的标量波传播,通常利用能代表耦合波场中P波分量运动学特征的拟声波(qP波)进行偏移成像.本文中,我们推导出了TTI介质下qP波控制方程.该方程可采用显式有限差分格式进行求解.通过声学近似,若沿对称轴方向的剪切波速度为零,对于对称轴方向不变且ε≥δ的模型来说,可得到稳定的数值解.但对于TTI介质来说,由于沿对称轴方向各向异性参数是变化的,声学近似会引起波场传播及数值计算的不稳定.因此,我们提出了正则化有限横波的方法,很好地解决了这一问题.最后,给出了Foothill模型的测试结果及某探区实际资料试算结果,展示了采用这个方程进行复杂TTI模型正演和高质量逆时偏移成像结果,证实了该方法的正确性和实际资料应用中的有效性. 相似文献
1