共查询到20条相似文献，搜索用时 15 毫秒
1.
时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质. 相似文献
2.
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征. 相似文献
3.
真实地下介质具有黏弹性,地震波在传播过程中会发生耗散与频散.忽视黏弹性介质的吸收衰减效应,逆时延拓过程中地震波将会出现振幅减弱、相位失真等现象,无法准确定位震源真实位置,因此需要对黏弹性介质中传播的波场进行衰减补偿,并通过采用合适的成像算子对微地震震源进行定位与裂缝成像.本文基于耗散与频散解耦的分数阶黏弹性波动方程模拟波场,采用low rank分解近似混合域算子,分离衰减相关项并反转耗散项符号,并在补偿的衰减项波场的波数域中进行低通滤波,压制噪声的影响;使用优化后的成像算子进行微地震震源定位,并通过分离散射波场,对散射波进行逆时反传寻找裂缝.数值实验证明,本文方法通过low rank近似有效提高了计算效率,衰减补偿算子在滤波器约束下能够稳定地补偿反向延拓的波场,优化后的成像算子能够在压制随机噪声的同时进一步提高计算效率和定位分辨率. 相似文献
4.
早期的褶积微分算子都是基于正反傅立叶变换而实现的,其精度比四阶有限差分的精度稍高,本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与16阶有限差分的精度相当,远优于错格伪谱法的精确度.另外,2.5维数值模拟比二维模拟可以更真实地模拟三维介质的臬个剖面的波场,并且2.5维地震波模拟的计算量比三维模拟的计算量及计算耗时要大大减少.本文利用基于Forsyte广义正交多项式褶积微分算子法计算2.5维非均匀介质地震波场,模拟结果表明,该算法的计算速度快,计算精度高,能够直观、高效地反映复杂介质中波场的传播规律,并且2.5维波场数值模拟具有更高的计算效率,是一种非常值得深入研究并广泛应用的方法. 相似文献
5.
分数阶微分算子具有描述历史依赖性和全域相关性的特质，本文利用这种特质描述双相介质固体骨架的黏弹性特征.基于Kjartansson常Q理论将含有分数阶时间导数的黏弹固体骨架各向异性本构关系与双相介质理论有机地结合起来，并引入流变学本构关系描述孔隙流体的黏滞性力学行为，提出一种新的基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI模型.推导了相应的时间域波传播方程，然后对该方程进行了数值模拟.对整数阶导数采用高阶交错网格有限差分算法，对分数阶时间导数采用短时记忆中心差分算法，进行了不同相界、不同品质因子组及双层地质结构情况下该类介质中波场的数值模拟与特征分析.模拟结果表明：将含有分数阶时间导数的常Q黏弹固体骨架各向异性本构关系及孔隙流体的黏滞性本构关系引入双相介质理论是可行的，二者的结合能更好地反映地下介质的黏弹性特征，对于进一步认识波在黏弹各向异性孔隙介质中的传播机理具有重要意义，为反演和重构地下油气储层和结构奠定正演理论基础. 相似文献
6.
地震波在地下介质传播过程中由于非弹性衰减的存在将导致能量损失和相位变化，精确的速度与衰减参数建模对油气识别、提高强衰减介质中地震波成像的质量都起着至关重要的作用.常分数阶拉普拉斯算子黏声方程由于完全分离的速度频散项与振幅衰减项的优势，以及在强非均质衰减介质中可以高精度求解的特点，已被应用于速度与衰减参数的建模中.本文将二阶常分数阶拉普拉斯算子黏声方程拆分为等价的一阶方程组，并在此一阶方程组的基础上推导出新的梯度公式与伴随方程，建立了一种新的速度与衰减参数同时重建的全波形反演方法.相较于原二阶常分数阶拉普拉斯算子黏声方程建立的全波形反演流程，数值实验表明，新建立的反演流程可以有效避免原梯度数值计算中的噪声，尤其是可以有效提高衰减参数梯度的反演精度，从而显著提高反演的收敛速度与反演精度. 相似文献
7.
常规τ值法假设应力松弛时间与应变延迟时间近似相等, 造成了常Q模型拟合精度低. 本文利用精确的广义流变体模型Q值计算公式, 研究改进的τ值法求解常Q模型参数. 根据地震波散射理论, 推导了基于广义流变体模型的黏滞性介质一阶波恩近似方程, 结合位移-速度关系得到了含卷积完全匹配层边界条件的黏滞性介质应力-速度方程的一阶波恩近似表达式. 通过数值实验验证并对比了黏滞性介质中全波波动方程、 一阶波恩近似方程以及单程波波动方程的波场特征, 讨论了基于流变体模型的黏滞性介质一阶波恩近似方程对速度扰动和Q扰动的适应性, 以及对旅行时和振幅精度的影响. 相似文献
8.
传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性. 相似文献
9.
三维波动方程时空域混合网格有限差分数值模拟方法 总被引：1，自引：0，他引：1
常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推导出了两类非坐标轴网格点构建三维Laplace差分算子的方法,进而提出了一种利用坐标轴网格点和非坐标轴网格点共同构建三维Laplace差分算子的混合网格有限差分方法,并利用时空域频散关系和泰勒展开建立差分系数方程,推导出了差分系数的通解.相比常规高阶和时空域高阶差分格式的2阶差分精度,时空域混合网格差分离散波动方程理论上能够达到任意偶数阶差分精度,模拟精度显著提高,同时稳定性更强.频散分析表明:相比常规高阶和时空域高阶差分格式,在计算效率基本相同时,时空域混合网格差分格式能更有效地减小数值频散,减弱数值各向异性,模拟精度更高;在模拟精度基本相当时,混合网格差分格式能采用更大的时间采样间隔,计算效率更高.数值模拟实例进一步验证了混合网格差分格式在提高模拟精度和计算效率方面的先进性,也验证了其普遍适用性. 相似文献
10.
基于声波方程扩充的哈密尔顿系统,本文给出了空间精度为八阶的近似解析离散化(NAD)保辛分部Runge-Kutta方法,简称八阶NSPRK方法。该方法采用八阶精度的近似解析离散算子近似空间高阶偏微分算子,并使用二阶精度的辛分部Runge-Kutta方法进行时间离散。我们从理论和数值计算两个方面研究了八阶NSPRK方法的稳定性条件和数值频散关系,并同四阶NSPRK方法、八阶Lax-Wendroff(LWC)方法和八阶交错网格(SG)方法进行了比较。结果表明八阶NSPRK方法压制数值频散的能力显著优于传统数值计算方法。与四阶NSPRK方法和传统四阶辛格式(SPRK)方法相比,八阶NSPRK方法具有最小的数值误差和最高的计算效率:在达到同样消除数值频散的前提下,八阶NSPRK方法的计算速度约为四阶NSPRK方法的2.5倍、为四阶SPRK方法的3.4倍;八阶NSPRK方法的存储量仅为四阶NSPRK方法的47.17%、为四阶SPRK方法的49.41%。在双层介质、非均匀介质和Marmousi等复杂速度模型中,八阶NSPRK方法模拟得到的波场快照非常清晰,无可见数值频散。这些结果表明,八阶NSPRK方法在粗网格条件下能有效地压制数值频散,从而能够极大地节省计算内存,提高计算速度。总体而言,八阶NSPRK方法是一种在地震探测领域和地震学研究中有着巨大应用潜力的数值计算方法。 相似文献
11.
地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性. 相似文献
12.
早期的褶积微分算子法都是基于正反傅立叶变换而实现的,其精度比四阶有限差分稍高。本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与l6阶有限差分相当。粘弹性波动方程更真实地描述了实际地下介质中弹性波的传播规律及其波场特征。本文以二维粘弹性波动方程为例,推导了粘弹性介质波动方程的离散格式,用迭积微分算子法实现了粘弹性介质的地震波场正演模拟,并对其波传播特征进行了分析。计算结果表明该算法能正确模拟粘弹性介质中的地震波,正确地反映粘弹性介质中波场的传播规律。 相似文献
13.
附加有效阻尼比的确定是消能减震结构设计的关键。为探究布置消能器楼层数的不同对结构设计安全性的影响,以新疆地区某多层钢筋混凝土框架结构为研究对象,采用非线性时程分析方法和等效结构模型的振型分解反应谱分析方法,基于各楼层剪力计算结果,引入楼层阻尼比修正系数kmn,对等效结构模型的附加有效阻尼比进行迭代修正。结果表明:采用小震时程分析阻尼比的等效结构模型,其楼层阻尼比修正系数随着阻尼器布置楼层数的逐渐增加呈现出先增后减的趋势,当阻尼器布置层数为楼层总数的3/5时折减率最低;采用中震时程分析阻尼比的等效结构模型,当阻尼器布置楼层数不少于2/5时,可满足设防要求;对于多层框架结构,考虑综合成本,黏滞阻尼器的布置楼层数可不采用满布的形式。 相似文献
14.
稳定性问题是地震波数值模拟的一个重要问题.基于地震波传播理论,从黏弹介质本构方程出发,对矩形网格下不同黏弹模型波动方程有限差分解的稳定性进行了理论分析,导出了Kelvin-Voigt黏弹模型和Maxwell黏弹模型在任意空间差分精度下稳定性条件的表达式;给出了品质因子Q≥5时的简化式,并通过数值算例验证了理论研究结论的正确性;总结了地震波速度、频率、空间网格大小、差分系数以及品质因子与稳定性条件的关系;通过误差分析给出了近似公式的使用条件. 相似文献
15.
首先对美国喷气推进实验室（JPL）公布的2007-06-01～2007-12-31时间段内的GRACE Level 1B卫星GPS轨道位置和速度、K波段系统星间速度、加速度计非保守力以及恒星敏感器姿态实测数据对应进行了轨道拼接、粗差探测、线性内插、重新标定、坐标转换、误差分析等有效处理;其次,基于改进的能量守恒法恢复了120阶GRACE地球重力场,在120阶处累计大地水准面精度为25.313 cm;最后,检验了本文地球重力场模型IGG-GRACE的可靠性,同时分析了IGG-GRACE解算精度在低频部分略优于德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型,而在中高频部分略低的原因. 相似文献
16.
17.
波动数值模拟中透射边界的稳定实现 总被引：21，自引：0，他引：21
从波动能量在计算区内累积增大的观点出发,通过简单的一维弹性波模型,系统地阐明了在近场波动数值模拟中透射边界两类数值失稳--"高频振荡"和"零频飘移"的机理:前者源于对波动数值模拟无意义的高频波动在人工边界上的放大和波在有限计算区内多次反射产生的反复放大;后者则源于透射边界允许零频和接近零频的分量不断进入计算区. 由此提出了稳定实现透射边界的完整方案包括两项简单措施:第一,在全部计算区内按文中建议的方法注入小阻尼,以消除高频振荡;第二, 给出一种具有明确物理意义的消除零频飘移的算子算法. 最后,提供了三维波源问题和散射问题的详细数值试验结果. 相似文献
18.
给出了任意倾斜各向异性介质中二维三分量一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式及其稳定性条件,并推导出了二维任意倾斜各向异性介质完全匹配吸收层法边界条件公式和相应的交错网格任意偶数阶精度差分格式. 数值模拟结果表明,该方法模拟精度高,计算效率高,边界吸收效果好. 各向异性介质中弹性波波前面形态复杂, 且qP波波速不总是比qS波波速快. qS波波前面和同相轴的三分叉现象普遍, 且其同相轴一般不是双曲线型. 当TI介质倾斜时,3个分量上均能够观测到横波分裂现象, 而且各波形的同相轴变得不对称. 相似文献
19.
《中国科学:地球科学》2021,(8)
基于L-S热弹性理论,采用旋转交错网格伪谱法,实现了均匀各向同性介质的一阶速度-应力-温度微分方程组的数值求解和波场模拟.其中,用时间分裂法解决方程组的刚性问题,用旋转交错伪谱算子计算空间一阶导数,用中心伪谱算子计算空间二阶导数;对于热导率变化比较大的双层介质模型和参考温度随深度按梯度分布的模型,用Crank-Nicolson显式方法取代旋转交错网格伪谱法进行计算;讨论了热耦合波场性质和传播规律,对比了常规伪谱法、交错网格伪谱法和旋转交错网格伪谱法热耦合波场模拟效果.数值模拟结果表明:基于L-S热弹性理论,用时间分裂法结合旋转交错网格伪谱法对均匀各向同性介质的波场计算,能够得到稳定的、高精度的模拟结果,但是在热导率变化剧烈的情况下不能用大时间步长进行求解,而且在参考温度不均匀分布的情况下算法不稳定.本文在网格剖分方式与数值算法的优化组合应用方面进行了探索,为将这些方法推广到孔隙热弹性、热黏弹性和各向异性研究奠定了基础. 相似文献
20.
瞬变电磁响应受到发射波形的影响,正演计算中需考虑发射波形.基于位移逆Krylov子空间投影的模型降阶算法能够精确模拟三维全波形瞬变电磁正演响应,但该算法计算精确的on-time响应非常耗时.目前实际工作中多数情况下都是对off-time时间段的数据进行处理解释,因此可以将正演问题简化为精确模拟考虑发射波形的off-time时间段的响应.本文实现了一种改进的基于位移逆Krylov子空间投影的模型降阶算法,能够实现任意发射波形的瞬变电磁off-time时间段的响应的三维快速正演.基于频谱分析理论实现了一种定量化的评价on-time采样波形精度的方法,能够在保证off-time响应精度的同时,极大的减少on-time波形采样点的数量,从而显著的减少正演计算时间.典型的半正弦波、梯形波和复杂的VTEM波的模型的正演结果表明,相比密集采样的正演结果,本文算法在保证off-time响应精度的同时,能够实现5倍以上的正演加速. 相似文献