衡山花岗岩体西缘韧性剪切带的特征,成因及岩体定位机制

衡山花岗岩体由南岳岩体(早期)和白石峰岩体(晚期)组成,衡山岩体西缘形成有分带明显的韧性剪切带,向西倾斜,在XZ面上呈倾向滑动剪切运动,矿物显微变形构造显示为倾向简单剪切特征。经测定,韧性剪切带形成时温度下限为343.6℃,古差应力下限为52～62MPa,应变速率下限为2.05×10~(-14)/s-3.13×10~(-14)/s,推算总位移量在13.5km左右。在花岗岩体以外没有韧性剪切变形现象,说明该韧性剪切带是衡山岩体定位过程中由白石峰岩体沿南岳岩体边缘侵位的简单剪切作用所形成。     

1962—1986年JYD1968.0地极坐标的均匀化及其与BIH的比较

本文研究了１９６２－１９８６年ＪＹＤ１９６８．０地极坐标均匀中星表系统差，章动系列的变动和光行星常数变动的影响，并将经过改正后的ＪＹＤ１９６８．０地极坐标系统和ＢＩＨ系统作了比较。     

偏微分方程数值解法在地学应用中的对比分析

在地学中广泛应用的偏微分方程数值解法有两种:有限差分法及有限单元法。对于定常态问题两种方法完全等价,对于非定常态问题有限单元法形成的代数方程归根到底仍是有限差分方程,但在一定条件下会引起反常问题,原因是与代数组相容的不是原来的热传导方程,而是反热传导方程。     

北京大海坨4.9级地震前白家疃井地热地形变前兆异常

北京大海坨ＭＬ＝４．９地震前，白家疃地热，地形变两种观测项目，同时观测记录到明显的异常变化，具有短期和短促，突发的临震异常特点，采取针对性的分析方法，结果较好，为单台项（前兆方法）综合预报提供了震兆观测实例。     

有限元离散模型中的出平面运动形式

波动在有限元离散模型中传播时,会出现许多不同于连续介质中波动的现象.本文引用经典物理学中晶格动力学的分析方法,由有限元离散形式的出平面波动方程得到频散关系,在此基础上,讨论了集中质量有限元模型中各种可能的出平面运动形式,并分析了考虑时域离散化影响下各种形式可能存在的频率范围.     

南亚高压强度与邻近地区垂直速度的相互依赖关系

利用ERA-Interim 1979—2018年6—8月的再分析资料,通过相关分析、信息流、合成分析等方法研究了南亚高压强度与其邻近地区垂直速度的相互依赖关系.结果表明：南亚高压在上对流层-下平流层区域具有上冷下暖的热力结构,冷、暖中心分别在70 hPa和250 hPa,以这两层的温度异常建立的温差指数可以反映南亚高压的强度.在不同时间尺度上,南亚高压强度与其邻近地区垂直速度的相互依赖关系是不同的.在月时间尺度上,南亚高压强度通过动力作用影响邻近地区的垂直速度,南亚高压增强（减弱）时,其东部地区的上升运动和西部地区的下沉运动同步增强（减弱）;在日时间尺度上,南亚高压中部的垂直速度通过热力强迫影响南亚高压强度,南亚高压中部地区上升运动增强（减弱）时,南亚高压增强（减弱）且位置偏西（偏东）.     

天山造山带及邻区孕震应力场及断裂带活动性特征的数值实验分析

天山造山带构造环境复杂,活动断裂带和强震分布广泛,且主要分布于阿尔泰山、天山、西昆仑—帕米尔弧形构造带上,尤以天山地区最为集中.迄今为止,天山造山带地区的主要断裂带的活动特征与孕震应力场特征之间的动力学机理尚未有清晰的认识.本文以GPS等实际观测数据为约束,建立有限元数值模型,计算了研究区域地壳形变、应力/应变积累速率、弹性应变能密度以及库仑应力变化率等关键因素.模拟计算结果显示地表速度场与研究区域实际GPS观测值基本一致,且主要断裂带上弹性应变能密度分布与实际地震活动性也基本吻合,验证了数值模型和结果的可靠性.结合最新的观测和数值模拟结果分析发现,研究区的断层和地震活动性主要受控于近南北向的主压应力,与主要观测特征相一致.同时,帕米尔高原北部边界带—塔什库尔干断裂（TKF）、天山造山带南边界的东侧—迈丹断裂（MDF）、兴地断裂（XDF）库仑应力增大明显,在未来强震发生的可能性较高,应予密切关注.     

Finite difference modeling of contaminant transport using analytic element flow solutions

The two-dimensional implementation of the analytic element method (AEM) is commonly used to simulate steady-state saturated groundwater flow phenomena at regional and local scales. However, unlike alternative groundwater flow simulation methods, AEM results are not ordinarily used as the basis for simulation of reactive solute transport. The use of AEM-simulated flow fields is impeded by the discrepancy between a continuous representation of flow and a typically discrete representation of transport, and requires translation of the flow solution to a discrete analog. This paper presents a variety of methods for analytically calculating conservative discrete water fluxes and integrated components of the dispersion tensor across cell interfaces. An Eulerian finite difference method based on these AEM-derived parameters is implemented for use in simulation of 2D (vertically averaged) solute transport. This implementation is first benchmarked against existing methods that use standard finite difference flow solutions, then used to investigate the effects of an inaccurate discrete water balance. It is shown that improper translation of AEM fluxes leads to significant water balance errors and inaccurate simulation of contaminant transport.     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

井地直流电法三维数值模拟中若干问题研究

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。