地电台址大范围介质条件的数值模拟研究

本文取两类震源模型,即“闭锁”模型和“膨胀”模型进行数值模拟,计算了大范围层状介质在多种不同组合情况下,其应力、应变场的分布,并进而计算了由附加应力、应变产生的视电阻率变化量。结果表明,不同介质组合对地电阻率前兆异常观测的影响不同,有利于发现前兆异常的大范围介质条件是:(1)坚硬层孕震;(2)孕育地震的高强度坚硬层要薄;(3)地表疏松覆盖层不能太厚。此外,计算结果还表明,地电台站探测范围内的介质力学,电性条件对前兆异常观测也有重大影响。     

西藏南部地堑构造成因的数值模拟

青藏高原上地壳在印度洋板块的挤压作用下发育出大量断层，因此上地壳的破碎已经从整体上降低了其有效粘滞系数。其所承受的水平向压应力逐渐转移到中、下地壳中。随着应力转移的进行，上地壳和上地幔内σ１可能已转为垂直方向，处于拉张减薄状态，并在藏南地区形成了一系列张性地堑构造。而目前在中、下地壳内σ１可能仍为近南北的水平方向。在印度洋板块的强烈挤压下，该深度范围内柔性的中、下地壳物质仍在垂直方向继续膨胀、增厚，并导致青藏高原的继续抬升。本文利用粘弹性蠕变本构关系的有限元方法模拟了上述设想的可能性。     

长周期地震动三维有限元数值模拟方法

采用数值方法对长周期地震动进行模拟已成为当前强地面运动预测的研究方向之一。本文提出了一种有限元并行计算的高效地震动数值模拟方法,可以同时考虑震源、传播途径和场地条件对地震动的影响,并根据昆明地区未来可能发生的设定地震对昆明市的影响地震动进行了预测。这一方法可以应用于国内已探测的发震活断层进行地震危害性评价工作。     

有限元等数值模拟方法在我国岩石圈构造研究中的应用与发展

本文综述了二十年来数值模拟方法在我国岩石圈构造研究中的应用与发展状况,主要涉及造山带数值模拟、俯冲带数值模拟、盆地数值模拟、上地幔小尺度对流数值模拟等几方面内容．为了研究中国大陆及其相邻地区的岩石圈三维结构．已经取得许多区域性地质调查和地球物理、地球化学等探测资料．同时以这些探测资料为基础,国内外科学家做了许多以该地区岩石圈构造与演化为研究目的数值模拟工作．可以看出,数值模拟是岩石圈构造与演化研究的有力工具．并已经取得了卓著的成就．这些成果为继续进行岩石圈构造的研究奠定了基础．     

基于双二次插值的探地雷达有限元数值模拟

从探地雷达(GPR)满足的波动方程出发,详细介绍了二维GPR模型单元剖分、二次插值、数值积分和有限元刚度矩阵总体合成的GPR有限元求解过程.为解决数值模拟时截断边界处的超强反射,采用Clay Bout透射边界条件对雷达波进行衰减,进而压制了来自截断边界处的反射波.在满足时间步长与空间网格差分稳定性前提下,采用中心差分法对GPR有限元方程进行离散,并用不完全LU分解预处理的BICGSTAB算法求解系数方程组,然后编制了基于双二次插值的GPR有限元正演模拟matlab程序.运用该程序分别对矩形和"V"字形两个典型地电模型进行正演计算,得到了正演剖面图,将该正演剖面图与基于线性插值的FEM算法的正演剖面图做了对比分析.结果表明基于双二次插值FEM算法相比基于双线性插值FEM算法异常响应更明显,具有更高的模拟精度,更有利于指导雷达剖面的数据解译.     

喜马拉雅西构造结及邻区岩石圈演化三维有限元数值模拟

用数值模拟的方法探讨了自10 Ma以来喜马拉雅弧形造山带的西构造结及邻区岩石圈的形变和位移变化特征,并采用不同的流变参数和流变结构模拟了前2 Ma内的形变和应力场的演化过程,通过对不同黏性系数、本构关系和初始地壳厚度的改变在演化过程中所起的作用研究表明：（1）采用黏弹性模型时,西构造结区黏性参数的选取对西构造结及邻区隆升高度和隆升范围有重要影响,计算结果给出西构造结区地壳黏性系数应小于1023Pa·s;（2）相对刚性的塔里木盆地采用弹性结构时,对整个西构造结区域的应力和应变场均能引致明显变化;（3）力学性质稳定块体的存在并没有使应力场的传播出现明显的解耦现象.     

东亚地区应力场的三维数值模拟

本文利用三维线弹性有限元方法模拟欧亚板块东部(主要是中国大陆地区)的三维构造应力场。计算中着重考虑东亚地区地形对该区地应力场的影响。同时考虑了对边界条件的限制,以及岩石在漫长地质时期中体现的流变性质,并在线弹性解上叠加了静岩压力等因素。 对一系列的计算结果进行数值分析,结果表明: 1.除各板块间的相互作用外,现代地形特征对东亚地区应力场分布情况有重要影响。 2.计算结果中,水平向最大主压应力轴走向在中国大陆地区呈辐射状分布,其中心大体上在青藏高原。随着与西部地区距离的增加,水平向最大主压应力值呈降低的趋势。应力方向和数值上的上述特征与震源机制解及原地应力测量结果是一致的。     

青藏高原现代构造应力状态及构造运动的三维弹粘性数值模拟

根据现有的地球物理和地质资料，建立了青藏高原三维弹－粘性有限单元模型，模拟了由地震震源机制解和地质调查资料推断的高原岩石圈的构造应力状态，并模拟岩石圈物质的运动特征，着重考察高原物质向东南方向运动的可能性，并把模拟结果与实际资料结果进行比较。施加载荷时考虑了处于Ａｉｒｙ均衡条件下的高原地壳的重力势能和底面浮力，以及水平边界的挤压作用，根据试错法的试验结果，能够符合观测结果的模型边界条件表明：①印度     

三维地震动场数值模拟并行计算系统

并行计算系统是进行三维复杂场地地震动场数值模拟的基础。本文作者利用10台微机构建了一套并行计算集群(常称为Beowulf系统),节点微机上运行Windows2000操作系统,采用MPI标准的MPICH1．2．4版本作为并行计算的支撑环境,并以100M bps高速交换式以太网作为互连网络。利用编制的并行计算程序对该并行计算系统性能进行了测试。结果表明：该系统能够达到非常高的并行加速比和并行效率。该系统将用于三维地震动场数值模拟研究。     

钻孔施工对新丰江和平地电台深井地电阻率观测影响的数值模拟

1 观测背景深井地电阻率观测在消除地表干扰、季节天气变化影响等方面效果良好(刘昌谋等, 1994),是当前地电阻率观测的发展趋势.广东新丰江和平地电台于1992年起开展深孔布极地电阻率观测,并在几次显著地震前观测到明显的异常变化.2015年3—5月,该台地电阻率测值出现显著异常变化,怀疑为测区钻孔施工对电阻率测值的变...     

大地电磁二维对称各向异性介质的有限元数值模拟

假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验     

电位双二次变化二维地电断面电阻率测深有限元数值模拟

以异常电位为研究对象,给出了一种新的二维地电断面电阻率测深的有限元数值模拟方法.该方法与以往这类方法的主要区别在于,网格单元中的电导率采用双线性插值,电位则用二次函数进行双二次插值,从而使得电位的正演结果以及由电位导出的视电阻率更加符合实际地电场的变化规律;另外,有限元网格单元数也大为减少.对几例模型的试算结果表明,文中提出的计算方法是行之有效的.     

电位双二次变化二维地电断面电阻率测深有限元数值模拟

以异常电位为研究对象,给出了一种新的二维地电断面电阻率测深的有限元数值模拟方法.该方法与以往这类方法的主要区别在于,网格单元中的电导率采用双线性插值,电位则用二次函数进行双二次插值,从而使得电位的正演结果以及由电位导出的视电阻率更加符合实际地电场的变化规律;另外,有限元网格单元数也大为减少.对几例模型的试算结果表明,文中提出的计算方法是行之有效的.     

三维有限元数值模拟在宁夏大坝电厂地震安全性评价中的应用

本文综合考虑了宁夏大坝电厂周围２５０ｋｍ范围内的地质、地球物理方面的有关资料，建立了三维有限元模型。利用该模型模拟了该区域的构造应力场，进而对大坝电厂厂址区及附近进行了地震危险性分析。在分析工作中结合安全度的定义提出了危险系数的概念及潜在震源区划分的定量－定性标志。经过分析，作者认为大坝电厂附近构造稳定性较好。     

地震孕育、发生、发展动态过程的三维有限元数值模拟

为探讨用数值方法模拟地震孕育、发生、发展动态过程的途径,从线性流变体介质内制约质点运动的运动方程出发,导出了模拟地震孕育、发生、发展动态过程的三维有限元方程及程序。还给出了模拟进程中地震孕育、发生、发展过程的约束条件,使得可以用同一个程序完整地模拟地震孕育、发生、发展的全过程,为地震过程本质的认识以及物理预报的实现提供了一个潜在的新手段。     

三维交流电测井的有限元方法

为了模拟三维地层中交流电测井的响应,本文考虑位移电流的影响,导出了用位函数A-Ф描述此问题的矩阵算子方程;证明了该矩阵算子是自伴算子;并由自伴算子的具势性得到了矩阵算子方程的泛函.最后用基于A-Ф的有限元方法计算了一维、二维及三维地层模型的交流电测井的响应.     

地震孕育、发生、发展动态过程的三维有限元数值模拟

为探讨用数值方法模拟地震孕育、发生、发展动态过程的途径,从线性流变体介质内制约质点运动的运动方程出发,导出了模拟地震孕育、发生、发展动态过程的三维有限元方程及程序。还给出了模拟进程中地震孕育、发生、发展过程的约束条件,使得可以用同一个程序完整地模拟地震孕育、发生、发展的全过程,为地震过程本质的认识以及物理预报的实现提供了一个潜在的新手段。     

连云港─银川岩石圈壳─幔组合结构剖面构造应力场的有限元数值模拟

从岩石圈壳－幔组合结构出发，选择了大华北岩石围块体内一条NWW向连云港－银川的岩石圈壳－幔组合结构剖面，试图先从一个简单的力学模型开始，用二维有限元法得到岩石围（地壳和上地幔顶盖层）中的构造应力场的分布特征，从而解释板内地震与　岩石围亮－幔结构比值（R）相关的力学原因。     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.