探地雷达时域有限差分法正演模拟

在以往雷达波的正演模拟中,借鉴地震波中成熟的正演模拟方法,均采用模拟声波方程的方法,精度不够高。因此从麦克思韦方程组出发,采用目前电磁场理论中最常用的时域有限差分法,对探地雷达进行了正演模拟。为了进一步说明其正确性,对正演模拟的结果又进行了偏移处理,从结果看出该方法的正确性及可行性。     

粘弹性介质中跨孔波场的交错网格有限差分法正演模拟

井间地震正演模拟技术是研究地震波在井间传播规律的重要手段之一,可以帮助认识井间地震的复杂波场.从二维井间地震波传播波动方程出发,结合初始、边界条件,推导出了交错网格任意偶阶精度差分格式,阐述了非均匀差分网格的实现方法.在此基础上,编制了跨孔波场交错网格有限差分正演程序,并应用该程序对半空间粘弹性跨孔模型进行了正演计算,得到了震源在不同位置时的波场快照及地震记录图.这对于跨孔模型波传播规律的研究,具有重要意义.     

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

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

时域有限差分法在地质雷达正演模拟中的应用

时域有限差分法是麦氏方程的差分形式的数值解法,用于求解任意复杂目标与电磁场的相互作用问题。本文把这一理论应用于地质雷达的数值模拟,得到满意的结果。     

波动方程有限差分正演模拟误差来源分析

有限差分波动方程正演数值模拟主要是解决差分格式、震源选择、边界条件处理等问题.这些问题的合理解决对得到正确和准确的波动方程数值解有着十分重要的作用.笔者通过对波动方程计算机有限差分正演模拟过程的叙述,分析了在做正演模拟时的误差来源以及如何有效地减小误差.     

可控震源粘弹性波动方程有限差分模拟

在一般的地震波计算中,地球介质可以当作各向同性的完全弹性体来对待。而对地表土壤和精细观测,粘弹性介质模型比完全弹性模型更符合实际。可控震源是近年来发展起来的精确可控、安全有效的新型震源,但国内对可控震源传播动力学的研究较少。对可控震源激发的地震波在粘弹性介质中的传播,用有限差分方法进行了数值模拟,给出了Kelvin一阶粘弹性波动方程组、吸收边界条件和震源处理。提出了检波器接收信号加速度和传播距离的三次多项式指数关系式,比线性指数关系更准确。使用该公式可以预测采集系统有效接收距离,指导测线布置。     

地震波场正演模拟中的有限单元方法

本文介绍了地震波场正演模拟中的有限单元方法。首先从粘弹性波动理论出发,推导了满足粘弹性介质的有限元方程,其次用不含阻尼矩阵的有限元方程模拟了地震记录。地表实例是自由界面的瑞雷面波。VSP记录有水平层、倾斜层;有零井源距和非零井源距;有线状和对称震源。 在程序中采用了超压缩存贮和分块计算措施,对降低有限元法要求计算机内存和运算速度的苛刻条件上做了有益的尝试。     

地震波正演模拟有限元法在微机上的实现

有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波传播过程,比较细致地再现地震图像。但由于有限元法地震模型所需网格节点数多、数据规模大,使它的计算受到限制。采用对角化集中质量矩阵、刚度矩阵变二维为一维存贮,并通过使用有限差分吸收边界条件来处理地质模型的人为边界反射问题,从而使计算所需内存降低,计算时间减少,并能比较好地实现地震声波用有限元法计算,理论模型的计算结果验证了算法的正确性和可行性。     

非均匀弹性介质地震波场伪谱法正演模拟

基于一阶速度—应力弹性波波动方程,利用傅里叶变换将空间导数变换到频率域求解,对介质波场进行正演模拟。本文结合实际勘探中可能存在的地质问题,设计了多种地质模型,包括背斜模型、向斜模型、透镜体模型和断层--侵入体模型等,进行了伪谱法地震波场正演模拟、合成地震记录的实现,并分析了其波场动力学变化特征。结果表明,应用伪谱法对复杂模型的模拟有较高的精度和可靠性,可直接用于辅助地震野外工作设计、反演解释和资料解释等工作。     

正交各向异性介质中地震记录的有限差分数值模拟

本文首先介绍了正交各向异性介质弹性波动方程的建立过程,然后讨论了高阶有限差分算法的实现原理、算法的稳定性和收敛性以及吸收边界条件的选取,并对正交各向异性介质作了正演模拟。通过波场快照、共炮点地震记录和ＶＳＰ模拟记录等正演分析,结果表明该算法在模拟各向异性介质中弹性波的传播是有效的     

基于频率域高阶有限差分法的正演模拟及并行算法

在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。     

探地雷达信号的高阶时间域有限差分模拟

探地雷达信号时域有限差分法模拟,多采用二阶精度的中心差分法近似（FDTD（2．2））,虽然计算简单,但数值色散误差较大,影响了模拟精度。在解决较复杂介质分布的探地雷达信号时,不能很好地反映信号的精细变化。而高阶时域有限差分模拟能减少数值色散带来的误差,提高了模拟的精度。采用各向异性完全匹配层（UPML）作为吸收边界条件,实现了高阶计算,有效地吸收边界电磁波反射,而且提高了计算效率。通过模拟的结果分析可知,高阶时域有限差分法能很好地提高模拟精度。     

复杂地电模型的探地雷达时域有限差分正演

目前探地雷达正演模拟。都是针对简单的层状模型、圆状空洞、正方形空洞等单一的规则模型,而对于地下弯曲的界面或“V”字形等复杂模型的正演,实现起来较为困难。然而,地下构造是复杂的,因此,对复杂地电模型的正演研究显得更为重要。这里以麦克斯韦两个旋度方程为基本出发点,运用K.S.Yee的空间网格模型理论和时域有限差分法的基本原理,推导出二维空间的探地雷达正演方程组,通过讨论数值频散关系及其产生原因,推导出符合探地雷达实际传播规律的理想频散关系。为解决正演模拟时截断边界处的超强反射,采用了Mur超吸收边界条件．并以自制的探地雷达正演模拟程序为依托,对两个复杂的探地雷达模型进行了正演模拟,得到了正演剖面图。从其后模拟实例中的两正演剖面图可以看出,边界截断处的干扰波被大大地减少,消除边界反射后的雷达剖面能更好地指导工作人员进行地质解释,从而达到了把雷达的正演研究深入到复杂地电模型中去的目的,使正演研究更符合实际的地质情况。     

基于球坐标系下有限差分的地磁测深三维正演

为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。     

均匀介质中F-K反偏移法正演模拟

F-K反偏移法正演模拟是在给定地质模型后通过构建人工偏移场,再对人工偏移场进行F-K反偏移来实现。人工偏移场的构建是在地质模型上直接放上经正确拉伸和幅度校正后的地震脉冲子波。F-K反偏移中插值映射采用两点sinc插值,其能避免假同相轴的出现且具有较快的速度。奇异点处理使用泰勒近似展开。最后对凹陷模型进行了F-K反偏移法正演模拟实现,得到了预期的效果。     

隐伏活断层地震勘探正演模拟方法

浅层高分辨率地震勘探是城市活断层探测手段中最有效、最可靠的方法之一,但是由于隐伏活断层勘探时存在着勘探深度浅、地层的属性差异小以及环境噪声大等因素,导致地震勘探数据信噪比低、反射信号弱等问题,严重影响地震数据的采集、处理和解释。采用一阶应力-速度方程的高阶交错网格有限差分方法, 结合完全匹配层(PML)的边界处理方法和通量校正(FCT)法消除频散的技术手段,对活断层模型地震勘探进行正演模拟,清晰直观地表现了断层处的细微差别,有助于地震勘探中对信号的分析与理解。结果表明,这些处理手段结合波动方程有限差分算法,能够有效实现薄层、弱反射的隐伏活断层模型的波场正演模拟,可有效识别隐伏活断层的断层倾角、断距、规模、跨越范围、断层深度以及断裂带影响范围,为分析这些参数对地震记录的影响提供有力帮助,同时结合实测地震数据的处理、解释,为防灾减灾提供具体有力的判据。     

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

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

有限单元法与有限差分法在MT一维正演模拟中的对比分析

首先从电磁场所满足的麦克斯韦方程组出发,介绍了大地电磁测深正演的基本理论,并针对一维大地电磁模型加以讨论.运用有限单元法及有限差分法分别推导了大地电磁测深一维正演算法,并运用Matlab 7.0软件编写了相应的程序.为了检验这两种一维正演算法的准确性,设计了均匀半空间模型和层状介质模型,并给出了由本文程序、解析解得到的相应结果和图件,从而对正演结果进行对比分析.结果表明:两种方法的正演结果均真实地反映了模型的地电参数.