共查询到20条相似文献，搜索用时 750 毫秒
1.
SH波在表面多层介质中传播的精确模拟 总被引：4，自引：2，他引：4
针对地震横波在地表低速层内的振幅放大效应问题，提出了一种模拟SH波在地表层状介质中传播的递推算法，并用它模拟了新西兰Alfredton盆地A10场址的SH波地震动响应特性。这个方法适用于具线性吸收性质的粘弹性介质。由于方法不受介质层厚薄制约，层厚可以无限薄化，实践上可以用许多薄层逼近的办法来模拟纵向上任意变化的连续介质。通过求取不同频率不同波数平面简谐波解并按实际问题的加权迭加可求解具特定波形和传播方向组合的任意SH波场。此方法在计算上具有解析解特有的精确性，稳定性和方便性，特别适用于模拟薄层介层，次波长现象及需要进行大量而又精确模拟计算的情形。 相似文献
2.
介质的弹性常数为三维四阶张量的分量,共有８１个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为２１个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有５个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的（即对称轴可以不平等于铅直方向）,在此情况下,需要进行坐标变换,如果已知介质在某一坐标系（其坐标轴平行或垂直于介质的对称轴）中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。 相似文献
3.
相速度和偏振方向是研究地震波传播规律和描述介质特性的重要参数,在理论研究和实际应用中有重要作用.本文假定倾斜横向各向同性(TTI)介质对称轴位于观测坐标系XOZ面内,在此观测坐标系下直接推导了TTI介质弹性波相速度和偏振方向的解析表达式,再进一步利用Thomsen弱各向异性理论,推导了弱各向异性近似条件下弹性波相速度以及qP波和qSV波偏振方向表达式.理论分析和数值试例表明,在相速度方面,随着各向异性介质参数改变,qP波和qSH波速度变化较为平缓,qSV波速度变化较为剧烈.弹性波相速度近似式误差均较小,能较好地近似精确相速度.在偏振方向方面,SH波偏振方向只是传播方向和对称轴倾角的函数,而与各向异性参数无关,SH波偏振方向既垂直于传播方向,又垂直于TTI介质对称轴方向.除特定方向外,qP波和qSV波的偏振方向与传播方向均成一定角度,并且随TTI介质对称轴倾角的改变而改变;在精确和近似情况下,qP波和qSV波的偏振方向始终垂直;在精度允许范围内,偏振方向的弱各向异性近似式与理论解析式吻合较好. 相似文献
4.
传统的Kirchhoff传播算子结构简洁,适用于描述横向均匀介质中波的传播.Ray-Kirchhoff传播算子较为精确地描述了波在非均匀介质中传播的运动学特征,其理论上的先天不足依赖于介质的复杂性.本文通过Born序列逼近波在非均匀介质中传播的大角度波分量,提出一种Born-Kirchhoff传播算子,将传统Kirchhoff传播算子的适用范围扩展至非均匀介质,同时描述波的运动学和动力学特征,其精度取决于Born序列逼近的阶数.利用Born序列频散方程,可以精确分析各阶Born-Kirchhoff传播算子对波长、传播角和非均质性的尺度依赖特征,其中,一阶Born-Kirchhoff传播算子的精度高于传统的相屏传播算子.波数域的Born-Kirchhoff传播算子对于高波数波是奇异的,导致波数域数值计算发散,但其空间域版本是非奇异的,无条件数值稳定,可通过Kirchhoff求和数值实施.本文给出各阶Born-Kirchhoff传播算子及其频散方程,可用于不同程度非均匀介质中的波传播模拟,复杂构造地震成像和速度估计.本文利用零阶和一阶Born-Kirchhoff传播算子计算简单二维模型的合成地震图,并与边界元法进行了比较. 相似文献
5.
由于构造运动等作用,TI介质对称轴往往沿空间任意方向分布,具有任意空间取向对称轴的TI（ATI）介质更符合实际地质情况.VTI介质与ATI介质的相速度在形式上具有一致性,VTI介质中地震波的相角对应ATI介质对称轴与地震波传播方向的夹角.本文基于Tsvankin的VTI介质精确相速度公式,利用TI介质对称轴和地震波传播方向上单位向量的数量积和向量积来计算ATI介质的精确相速度.根据弱各向异性假设,导出qP波和qSV波的近似相速度,分析了近似公式的误差,讨论总结了ATI介质qP波和qSV波的相速度特征.本文中的单位向量采用观测坐标系表示,通过相角关系,可以较为方便地由ATI介质近似相速度导出频散关系,然后借助傅里叶逆变换推导出时间-波数域qP波和qSV波解耦的波动方程.数值算例表明本文的波动方程是qP波和qSV波解耦的,波场计算结果稳定,未出现明显的数值频散,验证了本文方法的有效性. 相似文献
6.
我们发展了一种模拟复杂地表下含裂缝介质地震波场的方法,这对于解释山地地区的地震资料具有重要意义。基于Coates-Schoenberg方法,把裂缝引入到有限差分法(FD)中,从而使包含裂缝的单元里的弹性介质就具有了局部的各向异性。为了模拟起伏的地表地形,我们借助于贴体网格,将笛卡尔坐标系的具有水平对称轴的横向各向同性介质(HTI)的弹性波方程和自由边界条件变换到曲线坐标系中,采用一种稳定的、显式的二阶精度的有限差分方法离散(曲线坐标系)HTI介质中的弹性波方程。数值实例充分地展现了在不规则地球表面的影响下裂缝介质中地震波传播的复杂性。合成地震记录和波场快照表明裂缝端点产生的散射波在地表处会受不规则地表地形的作用,再次被散射;同理,地表地形产生的散射波,经过裂缝端点时也会被再次散射,尤其是瑞利面波产生的散射波,因其能量很强,严重污染了地震记录,使得识别地下裂缝等产生的有效信息变得异常困难。这对山地地震勘探中资料的解释具有重要意义。 相似文献
7.
结合振动波在土介质中的传播特性,在传统矩形空沟的基础上,设计6种不同截面形状的新型空沟.通过有限元建立各种新型空沟在不同面积、不同空沟排数及不同截面形式等工况下的隔振模型,研究各新型空沟在均匀、各向同性的弹性半空间下对地表稳态垂直简谐激励的隔振效果,从而探索如何通过改变波的传播方向、削弱波在原方向上的传播能力等措施,提... 相似文献
8.
实验室可以产生两种振动模式的横波,一种是剪切振动,另一种是扭转振动.本文利用扭转换能器观测了扭转横波在各向异性样品中的传播特性.实验表明,无偏振方向的扭转波在各向异性介质中会以两种不同速度传播,并出现波的分裂现象,其快慢波的速度值与剪切横波的快慢横波速度一致.用扭转波换能器接收时,快慢扭转波的波形振幅不受各向异性方位影响,表现为无方向性.通过两块不同来源的均匀各向异性样品,用实验观测方式初步揭示了各向异性介质中扭转波的一些传播特征. 相似文献
9.
本文给出了重力波在中层大气温度波导中的导制传播模型,并在此模型的基础上详细讨论了重力波部分导制传播下的对称模式与非对称模式,导出了不同模式下相应的特征函数和色散方程,进一步用离散的方法对两类色散方程进行了求解;同时还利用二维全隐欧拉格式（FICE）对重力波在温度波导中的传播进行了模拟,模拟的结果也成功地展现了对称与非对称两种传播模式.研究表明,下边界的扰动能量在向上传播进入波导区域后被俘获,形成导制传播.不同周期的初始扰动,在波导内均会形成对称与非对称形式两种模式的导制传播,由于两者的行进速度不一致,最终会引起两种不同模式的分离.数值模拟中重力波的水平行进速度与线性模型预测值非常接近.波导中不同模式下重力波的水平波长与初始扰动的水平波长非常一致,然而波导区域内重力波的频率与初始扰动的频率无关,频率不同的初始扰动会激发出相同频率的重力波对称与非对称导制传播模式.这表明在确定的温度波导中,水平波数才是决定重力波传播特性的决定因素.进一步的分析显示,初始扰动的水平波数-频率分布越接近完全导制传播的色散关系时,温度波导中更易于生成以该种模式部分导制传播的重力波. 相似文献
10.
本文介绍实际地球介质（地球物理介质、地质介质）中的波的激发与传播问题的某些新的研究方向与进展。首先简要说明实际地球介质的一般特性。其次分别介绍随时间变化的三维非均匀地球介质中的地震波和非线性地球介质中的地震波的传播问题，也就是四维地震学和非线性地震学的研究进展。本文最后研究了三种类型（广义瑞利型、广义勒夫型和广义导波型）的构造活动波的激发与传播特性，并用它们来解释过去已发现的几种不同类型的地震活动波。 相似文献
11.
基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段. 相似文献
12.
当地下介质存在各向异性时,在观测坐标系下的弹性参数与自然坐标系下的弹性参数不一定相同．首先,根据势能密度和耗散能密度与坐标轴无关的原理,推导出了双相各向异性介质中观测坐标系下弹性参数与自然坐标系下弹性参数之间的关系;然后,从任意双相各向异性中弹性波波动方程出发,得出了该方程的伪谱法数值解法;最后,通过数值模拟,观测到了存在于双相各向异性介质中的４类波,即快纵波、慢纵波、快横波和慢横波．在双相各向异性介质中,ＳＶ波传播的波前面上仍然存在波面尖角,这些尖角在界面上要发生反射和透射．另外,数值模拟结果中可见转换慢纵波和慢纵波的转换波． 相似文献
13.
将各向异性介质分成各向同性背景介质和各向异性异常介质,并提出一种空间-波数混合域方法实现了各向异性介质下直流电阻率法的三维数值模拟.不同于传统直流电阻率数值模拟方法,本文算法直接对空间域异常电位满足的偏微分方程沿水平方向进行二维傅里叶变换,使水平方向转换成波数域,保留垂直方向为空间域,可根据地下介质电流密度变化的快慢灵活剖分.这样可把空间域异常电位满足的三维偏微分方程转化成不同波数满足的一维常微分方程,把一个大规模三维数值模拟问题分解为多个一维数值模拟问题,利用一维有限单元法求解方程组,并通过采用压缩算子迭代计算,最终获得较为精确的数值解.与自适应有限单元法对比验证了本文算法的正确性;测试了算法的收敛性,结果表明在满足精度要求的情况下,算法的收敛性只与异常体和围岩之间的电导率差异相关,而与异常体大小和埋深无关;分析了算法计算效率,结果表明算法的计算效率与剖分网格节点成线性关系,算法可在微型计算机中较快计算出剖分节点总数超过千万的各向异性模型的结果;设计简单观测系统并验证其具有反映地下各向异性结构特性的能力;最后模拟异常体沿着不同方向旋转不同角度时的响应特征,对比分析可知异常体为各向异性... 相似文献
14.
《地震研究》2017,(4)
准确模拟地震波在复杂介质中的传播过程,能够为分析主动源探测中获得的数据资料以及反演区域内部介质结构提供理论支持和依据。基于二维介质模型,利用有限差分法对地震波场传播过程进行数值模拟。结果表明,震源在不同介质中被激发后所产生的波场特征不同,有限水体的存在对波场传播有显著影响,随着水体水位的加深,激发产生的波场能量先增强后减弱,存在一个最佳水位深度,此时激发产生的波场能量最强。在有限水体一侧激发震源时,两侧等震中距台站的记录会有所差异,靠近陆地的一侧的波场能量更大,传播速度更快,且水平分量的差异相比垂直方向更明显。因此,可以通过数值模拟判断水库水位的最佳深度与震源的最佳沉放位置,为水库激发气枪震源实验提供理论依据。 相似文献
15.
Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利用同时处理这两种力学机制的BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程中出现的各种波动现象． 相似文献
16.
三角网格有限元法能够准确模拟复杂构造和复杂介质条件下的地震波场,数值频散和稳定性条件是地震波数值模拟中参数选择的主要依据.基于均匀的线性三角网格单元,根据结构刚度矩阵的组装原理以及平面波理论,推导了集中质量矩阵下两种网格结构的声波频散函数以及稳定性条件,并对数值频散特性以及稳定性进行了详细研究:三角网格单元中波动的数值频散除了受到空间采样间隔、单元网格纵横比和波传播方向等常规因素的影响外,还受到网格布局的影响,过锐或过钝的三角单元会对波动数值频散产生不良的影响,不同类型的单元网格、单元纵横比对应着不同的稳定性条件,正三角单元中的波动具有较好的数值频散特性,其数值各向异性(频散随波传播方向的变化)效应最弱,稳定性条件也较为宽松.最后通过数值模拟直观地验证了以上分析结果,为有限元正演三角网格的剖分和参数的设置提供一定的理论依据. 相似文献
17.
《地球物理学进展》2015,(1)
BISQ模型同时包含了合流体孔隙介质中Biot流动和喷射流动两种重要的力学机制.基于BISQ模型的双相各向同性介质弹性波波动方程是一个复变系数偏微分方程组.本文率先建立了该方程的25点频率-空间域有限差分数值解法,在理想相界和黏滞相界情况下,对Biot流动和喷射流动共同作用下的双相各向同性介质中的波传播进行了数值模拟,通过与仅受Biot流动作用下的波场模拟结果的对比,分析了两种机制耦合作用对弹性波传播的影响.同时,本文也研究了波在双层双相各向同性介质分界面处的反射和透射特征.数值模拟结果表明:在Biot流动和喷射流动耦合作用下,双相介质中传播的快P波的速度和振幅都小于仅考虑Biot流动的双相介质中的快P波,且慢P波的衰减也更为强烈,而S波的波速和振幅则无明显差异.这表明局部喷射流动对P波的衰减和频散具有重要影响,而对S波的影响较小;慢P波的强烈衰减使得其在波场快照中无法被观测到,双层双相介质中的波传播现象类似于单相介质的情况.同时本文的研究结果也表明,频率-空间域有限差分法在基于BISQ模型的双相介质中波传播数值模拟中的正确性和有效性,为开展孔隙弹性介质全波形反演问题的研究提供了研究基础. 相似文献
18.
为了实现交错沉积等复杂环境中的电磁场数值模拟,本文在常规横向同性模型的基础上引入了电导率主轴坐标系相对地层坐标系的层理方位角和倾角,建立了交错地层模型.并利用传播矩阵法建立了一维层状交错地层模型中的多分量感应测井仪器响应的正演模拟算法.首先将频率-波数域中的电磁场分解为上行和下行模式波,给出了任意朝向的磁偶极子在无限大地层中模式波的解析解.进一步通过引入地层界面上的透射、局部反射以及广义反射系数矩阵,推导了一维层状地层中的模式波表达式.在此基础上,利用二维Gauss-Legendre积分实现了Fourier逆变换,得到了可用于多分量感应测井模拟的频率-空间域磁场并矢格林函数.最后,通过多个数值模拟结果考察了井眼倾角、层理方位角和倾角变化对多分量感应测井响应的影响. 相似文献
19.
旋转轴对称介质中的qP波射线追踪 总被引：1，自引：1，他引：0
本文使用qP波一阶射线追踪方程(FORT)计算光滑、非均匀旋转轴对称弱各向异性介质中qP波传播的路径和走时.此FORT方程只依赖于15个弱各向异性参数,而非标准射线方程中的21个弹性参数.通常弹性参数模型是在局部坐标系中给定的,而在实际中需要的是全局坐标系下的弹性参数,因此为了解决两个坐标系下弹,性参数的变换问题,本文... 相似文献
20.
地震学百科知识(五)——地震各向异性 总被引：1，自引：0，他引：1
1基本概念地震波在地球的各向异性介质中传播时,其传播速度与质点偏振方向等特性随波的传播方向而变化的现象,称为地震各向异性。地震各向异性通常表现为三个方面:①波的传播速度随传播方向而变化;②波的传播速度随波动的质点偏振方向不同而发生改变;例如,S波经过各向异性介质后会分 相似文献