共查询到20条相似文献，搜索用时 46 毫秒
1.
本文对董良国、马在田等人发表在<地球物理学报>上的一阶弹性波方程交错网格高阶差分解法一文进行了差分格式的改进.发现原文对一阶空间导数采用2N阶差分精度,而对高阶导数却采用简单的中心差分,即对同一方程中的各阶导数的处理存在有的采用高阶格式、有的采用低阶格式的不一致性问题,本文的改进点正是消除了这种不一致性,建立了时间和空间均为四阶精度的差分格式.该方法具有编程简便、易于计算机实现和精度较高的特点. 相似文献
2.
3.
在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高. 相似文献
4.
给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式. 数值模拟结果表明,该方法模拟精度高,计算效率高,边界吸收效果好. 各向异性介质中弹性波波前面形态复杂, 且qP波波速不总是比qS波波速快. qS波波前面和同相轴的三分叉现象普遍, 且其同相轴一般不是双曲线型. 当TI介质倾斜时,3个分量上均能够观测到横波分裂现象, 而且各波形的同相轴变得不对称. 相似文献
5.
6.
《地球物理学进展》2016,(1)
为克服各向同性和VTI介质逆时偏移方法对复杂地质构造成像的局限,研究了TTI介质拟声波逆时偏移方法.首先从精确的TTI介质频散关系出发,引入一个各向异性控制参数σ,推导了新的二阶耦合TTI介质拟声波方程,以保证波场延拓的稳定性;然后引入波场的伪速度分量,推导了等价的一阶拟声波方程.相比于规则网格有限差分法,交错网格有限差分(SGFD)法能够有效地压制数值频散,模拟精度更高;因此利用高阶SGFD法求解TTI介质一阶拟声波方程,构建逆时偏移所需的正向和逆时波场延拓算子,并应用归一化互相关成像条件实现精确的TTI介质逆时偏移成像.最后,简单讨论了伪横波的产生机制,并给出了伪横波的联合压制策略.模型试验结果验证了方法的有效性和稳定性. 相似文献
7.
一阶弹性波方程交错网格高阶差分解法 总被引：150，自引：34，他引：116
提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。 相似文献
8.
有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。 相似文献
9.
地震波场数值模拟是理解地震波在地下介质中的传播特点，帮助解释观测数据的有效手段，而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法（Flux|corrected transport method，FCT）相结合, 对横向各向同性介质（Transverse isotropic medium，TI）一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明，本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比，在压制网格数值频散方面有明显的优势，计算精度提高，而且可以利用较大的空间步长，提高计算效率. 相似文献
10.
密度是岩石物性参数的重要组成部分,对储层评价、岩性解释和油藏描述等具有重要作用,因此实际资料应用中含有密度信息的正演建模是必不可少的.作为波动理论应用最广泛的实际资料处理技术——逆时偏移和全波形反演,其正演建模一般采用经典的二阶标量波方程,且不考虑密度的空间变化.而一阶速度-应力控制方程通过交错网格有限差分法正演建模,其精度高于二阶方程但计算成本过高,很难在三维实际资料中应用.鉴于非均质正演在实际生产中所面临问题,本文借用交错网格的思想,充分考虑空间变密度对波传播的影响,提出基于等效交错网格的高阶有限差分建模算法并应用于非均质标量波正演.文中从数学上给出一阶方程和二阶方程在数值模拟中的等价性证明,并分析了震源、边界条件和稳定性.最后通过简单层状模型进行数值测试,对比不同建模方法以验证新方法的准确性;利用高阶精度方法对Sigsbee-2a模型做正演测试,验证本文提出的方法的稳定性. 相似文献
11.
12.
柯坪塔格前缘断裂东段是柯坪推覆构造系前缘的一条活动断裂,野外调查获得了其晚第四纪错断洪积扇、冲沟阶地面的证据,实测了变形地貌面上的断层陡坎,分析了断层的形变量,通过采样测年估算了断层的缩短速率。由7个观测点的断层陡坎剖面测量,计算了观测点处断层的水平缩短速率,结果表明,断裂弧顶部位的五道班、三间房以西及其大山口道班附近,断层错断了Ⅰ级和Ⅱ级洪积扇(阶地)。断层在这些地点最新活动强烈,水平缩短平均速率全新世以来为0·35~0·44mm/a,更新世晚期末以来为0·16~0·30mm/a,而在非弧顶部位的巴楚磷矿、三岔口以北及大山口北断层只错断了更新世晚期Ⅲ级洪积扇,且水平缩短速率较小,断层水平缩短平均速率更新世晚期以来为0·05~0·07mm/a 相似文献
13.
大震以后,跟随着一系列的余震,其震源机制有的和主震一致,有的则大不一样,这表现在P波初动的反向,反应了震源应力场的变化。 余震震源机制变化的原因,既不是岩块的错动过头,也不是液体流出的下陷,而是发生大震的滑动岩块对其前后邻接岩块产生的转换应力场所引起的剪切破裂。 相似文献
14.
15.
16.
17.
18.
在不同含水条件下研究了光面、含脆性物质夹层和含延性物质夹层 3类断层的摩擦滑动行为 ,认为水对断层滑动稳定性和摩擦滑动的速度依赖性有重要影响。建立了含有岩性、围压 (深度 )和含水量 3因素的断层活动稳定性的模式 ,认为脆性断层更多地表现为速度弱化 ,少量水的存在易于出现黏滑 ,高孔隙压时会表现为稳滑 ;半脆性或半延性的断层在少量水或较低孔隙压的情况下会表现出速度弱化 ,也就可能出现黏滑 ,但孔隙压较高后 ,会表现出速度强化 ,滑动也就会是稳定的 ;延性断层多为速度强化 ,不出现黏滑。提出水的进入降低了系统的稳定性 ,水的存在增加了系统的稳定性。这是考虑水或一般流体对岩体或断层活动稳定性影响时需要考虑的两个方面 相似文献
19.
20.
采用大理岩圆柱试样,与轴线约成30角预制不连续裂缝,再用环氧树脂胶接,以实验观测裂缝尺度及其填充物对剪切破裂能的影响.结果表明:当裂缝相对面积从0增加到90%时,破裂能从(6.92.1)103焦耳/米2逐步降到(0.60.5)1O3焦耳/米2;而使用 Knopoff方法和引入滑动弱化模式后的方法,求得的理论应力降几乎不变.它们与实测应力降的比却逐步增大并达到10.同时,分析了地震时和岩石破裂时的应力降差别很大,以及不同方法所得破裂能差别很大的原因,认为断层的蜂窝状接触对上述问题可能有重要影响. 相似文献