固体潮应变各向异性研究

钻孔应变仪的观测过程中发现潮汐应变响应存在方向各向异性, 推测是受到台站附近断层影响导致的。 我们对此展开数值模拟研究, 讨论软弱断层对潮汐变形传播的影响。 我们的研究结果证实了软弱断层的存在确实能够产生能够观测到的各向异性现象。 同时我们还发现其各向异性变化对断层弹性参数较为敏感。 观测到的部分地震震前反映固体潮应变各向异性的潮汐因子"玫瑰图"变化有可能是由于台站附近断层在震前的性质变化引起的。     

巴颜喀拉块体周缘几次强震对甘肃地区地震形势的影响

以青藏高原为目标,采用黏弹性体本构关系建立3D有限元模型,模拟了在印度板块持续向北推进、挤压下,巴颜喀拉块体周缘强震对邻近地区的影响.通过降低强震发生位置处单元弹性模量的方法模拟了2001年昆仑山口西MS8.1、2008年汶川MS8.0、2010年玉树MS7.1地震,计算了发生这3次MS7.0以上强震后50年内每5年的黏弹性调整过程,讨论了3次强震对祁连地震带及甘东南地区地震形势的影响,尝试开展依据区域应力调整特征判定未来强震危险区的探索试验.     

井孔水温对远场巨震同震响应及其机制的数模研究

热量的流动必然伴随着温度的改变,地震活动期间的井孔水温同震响应也服从热力学的基本规律。本文以汶川8.0级地震和日本福岛9.0级地震为例,以海口ZK26井水温同震响应的观测数据为基础,运用热力学传导方程的数值模拟方法,研究了同震响应过程中井孔水温度变化与热量传导之间的关系及相关的热力学机制问题。通过正演的方法,得出沿井深方向不同时刻水温同震响应变化的数值模拟曲线,水温数值模拟曲线与实际测量数据一致。研究结果表明,井孔中水温同震响应方式（上升、下降或不变）,与水温传感器的位置、热源的位置和分布、传感器与热源之间的相对空间等因素有关。     

不同初始扰动对电离层扩展F影响的数值模拟

重力波、中性风场、电场是激发电离层扩展F的主要影响因子,本文基于中低纬电离层扩展F发展的物理模型,通过电场强度、背景风场对扩展F影响作用的分析和经验对比,首先验证了模型的有效性,后借助该模型数值模拟了给定背景环境下三种尺度初始电子密度扰动条件下扩展F的发展情况,同时研究了利用化学物质释放实现一定尺度扰动,进而激发扩展F的过程.结果表明,较强的背景电场、东向风场有利于扩展F的形成和抬升,与经验结论相吻合;电离层从被作用初始扰动到激发扩展F的过程中存在拐点效应,拐点之后扩展F被激发形成并且抬升迅速,同时短波长扰动相对于长波长扰动更有利于扩展F的激发和发展;化学物质H₂O释放通过耗散电子密度,形成了一定尺度扰动并诱发了扩展F的形成,该方法可作为一种人工激发扩展F的探索手段.

     

鲜水河断裂带北西段跨断层基线变化机理探讨

2013年芦山M_S7.0地震发生后,鲜水河断裂北西段的侏倭、虚墟等跨断层测点基线测项均出现不同程度的异常变化,异常核实结果认为基线异常属实.该异常变化是否与芦山地震的发生存在一定关联性,值得深入研究.本文首先分析芦山地震发生前后鲜水河断裂北西段不同基线长度的变化特征.基于鲜水河断裂带及邻区三维非线性黏弹性有限元模型,以GPS观测资料、同震静态滑移量作为约束,通过开展多组数值模拟实验,探讨下地壳不同流变特征下,芦山地震发生前后不同阶段,不同跨断层基线长度的时序动态变化特征.初步研究结果表明,1）侏倭、虚墟基线测点原始观测资料均表现出地震发生时基线长度减小,发生后短时间内反向快速恢复及增加的协同变化;2）在以GPS观测数据为约束的块体间差异性运动的动力学边界条件下,鲜水河断裂表现为左旋走滑运动;3）模拟结果显示,地震发生时测点间距离迅速减小,表现为断层右旋的特征,这与实际观测资料反映的结果一致;4）鲜水河断裂带两侧地块的下地壳黏滞系数分别取10¹⁸Pa·s、10¹⁹ Pa·s时,由于黏弹性松弛效应的影响,造成基线长度在地震发生后短时间尺度内快速增加.对比分析认为,地震后短时间内实测资料反映的基线长度快速增加的特征可能是黏弹性松弛效应与构造作用共同作用的结果.

     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

Influence of spatial variations of microtopography and infiltration on surface runoff and field scale hydrological connectivity

Surface runoff on agricultural fields arises when rainfall exceeds infiltration. Excess water ponding in and flowing through local microtopography increases the hydrological connectivity of fields. In turn, an increased level of hydrological connectivity leads to a higher surface runoff flux at the field boundaries. We investigated the functional hydrological connectivity of synthetical elevation fields with varying statistical properties. For this purpose, we developed an object-oriented ponding and redistribution model to which Philip’s infiltration model was coupled. The connectivity behaviour is determined by the presence of depressions with a large area and spatial organization of microtopography in rills or channels. The presence of microdepressions suppresses the effect of the spatial variation of infiltration properties. Connectivity behaviour of a field with a varying spatial distribution of infiltration properties can be predicted by transforming the unique connectivity function that was defined for a designated microtopography.     

A 3D gravity data interpretation of the Matagami mining camp, Abitibi Subprovince, Superior Province, Québec, Canada: Application to VMS deposit exploration

Two inversions, unconstrained and constrained, of a gravity survey of the Matagami mining camp (Abitibi Archaean Subprovince, Canada) have been performed in order to identify the downward extension of a rhyolitic horizon hosting VMS-type base metals deposit and the morphologies of the major felsic plutons. A comparative study exhibits the similarities between measured and calculated densities from chemical compositions of the Matagami lithologies. This allows building an initial 3D geodensity model which integrates densities and available structural and geological surface mapping data. This model is integrated during the iteration process of the constrained inversion in the objective function. The resulting true density model and two derived cross-sections upgrade the 3D imaging of this area. Also, the model gives new insight for regional geological interpretation exposing possible shapes of the main geological units at depth and suggests the potential existence of deep fertile geological bodies.     

Numerical simulation of deep-hole resistivity anomaly caused by drilling construction in Xinfengjiang geoelectric station

The Xinfengjiang reservoir in Guangdong Province, is one of the large reservoirs that have triggered earthquakes of magnitude greater than 6. Numerous earthquakes have occurred since the impoundment of the reservoir, making it one of the most active seismic zones in south China. In 2015, a set of deep-hole resistivity anomalies was observed in the Heping geoelectric station in Dongyuan county, located near the Xinfengjiang reservoir. After a field investigation, we found that a planned well drilling construction of new measuring channels was being carried out during that corresponding period of time. After careful comparison and analysis on the basis of the collected raw data, we had a reason to believe that drilling construction, rather than the inducement of the Xinfengjiang reservoir, was the main culprit for those unusual georesistivity values. So as to verify the above conjecture, we constructed a series of 3D finite element models based on the geological and hydrological information around Heping station and analyzed the drilling disturbances, respectively. Some significant conclusions were finally drawn according to the precise numerical simulation. This study gives a good example by combining numerical simulation with engineering practice as a way to understand the root cause of georesistivity anomalies in reality.