The forward problem of electromagnetic induction: accurate finite-difference approximations for two-dimensional discrete boundaries with arbitrary geometry

We describe finite-difference approximations to the equations of 2-D electromagnetic induction that permit discrete boundaries to have arbitrary geometrical relationships to the nodes. This allows finite-difference modelling with the flexibility normally ascribed to finite-element modelling. Accuracy is demonstrated by comparison with finite-element computations. We also show that related approximations lead to substantially improved accuracy in regions of steep, but not discontinuous, conductivity gradient.     

Transient electromagnetic responses in seafloor with triaxial anisotropy

Electrical anisotropy of young oceanic crust at mid-ocean ridges is detectable by observation of the rate and geometry of the diffusion of electromagnetic fields. The anisotropy in electrical properties arises from the presence of conductive seawater in an interconnected network of mostly ridge-parallel cracks. In this paper, we first justify the choice of a triaxial model to represent young oceanic crust, with three distinct electrical conductivities in the vertical, strike and spreading directions. We then present an algorithm to calculate the transient electromagnetic responses generated by an electric dipole source over such a triaxially anisotropic seafloor. We show that if the transient passages are measured with three distinct electric dipole-dipole configurations, it is possible to discern all three unknown conductivities independently of each other.     

A general correspondence principle for time-domain electromagnetic wave and diffusion fields

A general correspondence principle is presented that relates any time-domain electromagnetic diffusion field to an electromagnetic wavefield in a 'corresponding' configuration. The principle applies to arbitrarily inhomogeneous and anisotropic media and arbitrary transmitters and receivers. For the correspondence between the two types of electromagnetic fields to hold, the electric conductivity in the diffusive case and the permittivity in the wavefield case should have the same spatial variation, while the permeability distributions in space in the two cases are to be identical. Essential steps in the derivation of the correspondence principle are the use of the time Laplace transformation of causal signals, taken at real, positive values of the transform parameter, the Schouten-Van der Pol theorem in the theory of the Laplace transformation, and the reliance upon Lerch's theorem of the uniqueness of the interrelation between causal field quantities and their time-Laplace-transform representations at real, positive values of the transform parameter. Correspondence is then established between the tensorial Green's functions in the two cases, where the Green's functions are the point-receiver responses (either electric or magnetic field) to point-transmitter excitations (either electric- or magnetic-current source).
Through the correspondence principle, all transient electromagnetic wavefields (where losses are neglected) have as a counterpart a transient diffusive electromagnetic field (where the electric displacement current is neglected). The interrelation yields the tool to compare quantitatively the potentialities of the two types of fields in transient electromagnetic geophysical prospecting.
Finally, a general medium-parameter scaling law for time-domain electromagnetic wavefields is presented.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

Genesis and distribution of mineral waters as a consequence of recent lithospheric dynamics: the Rhenish Massif,Central Europe

Three major, interdependent processes control the genesis and distribution of mineral and thermal waters in the Rhenish Massif, Central Europe: (a) Magmatic processes in the upper mantle provide most of the CO₂ to produce bicarbonate waters in shallow aquifers. (b) Extension of the brittle upper crust enables the ascent of sodium chloride waters. (c) Uplift and erosion shape the massif's relief, which determines the extent of flow systems and the distribution of thermal springs. The chemistry of mineral waters further depends on the aquifers' mineral composition. A comprehensive set of hydrological, chemical, tectonic and geophysical data on the Rhenish Massif has been compiled. It was used to classify the mineral waters and to map the spatial distribution of water properties. The composition of cuttings from several representative wells producing different water types shows that the hydrothermal alteration of the aquifer rocks consists mainly of kaolinization of chlorite and dissolution of feldspar. Numerical transport simulations favour two modes of groundwater flow: topography-driven flow and the pressure-driven ascent of basement brines along active faults. Thermal convection is less important.     

变质反应的代数识别方法

在回顾变质矿物共生组合和变质反应图解分析法的基础上，论述了录求变质矿物组合之间、同一组合内矿物之间的平衡变质反应的代数方法，即矩阵高斯变换，线性回归分析、矩阵奇值分解。     

内蒙大水清金矿带交代作用及其与金矿化关系

大水清金矿带是华北克拉通中有代表性的改造型绿岩金矿之一。它产于前东武纪绿岩地体中安家营子石英二长岩体内外接触带附近的剪切带中。金矿化交代主要受后阶段张性脆裂控制，而剪切带早阶段韧性剪切作用无明显矿化交代，后者主要为成矿交代提供有益的江游交代构造环，境即成矿交代蚀变作用和剪切带同步的两阶段演化模式。大水清金矿带的交代作用有四个阶段：即用长石化阶段，绿泥石化阶段，黄铁绢英岩化和石英－硫化物阶段以及石英－碳酸盐化阶段。第三阶段由五个亚阶段组成：即黄铁绢英岩化和强硅化－强黄铁矿化，石英－多金属硫化物交代，石英－贫扶矿充填交代和细粒黄铁矿充填交代。其中前四个亚阶段是本区金的主要成矿期．当韧性剪切转化为张性脆裂时产生的骤然扩容和伴随发生的交代体系转变为富硅体系的钾交代和石英－硫化物交代是金矿成矿的关键。本文还研究和探讨了该金矿成矿交代流体的性质及交代的物理化学环境。     

蚀变带元素地球化学成矿预测法──以金龙山微细浸染型金矿为例

金矿化作用伴围岩蚀变，且多受断裂控制，反映在地球化学上出现某些元素的异常。在已发现一害矿化体并掌握一般控制矿特征的地区，利用已有工程在蚀变带内连续捡块采样，分析有关批示元素，总结有关的地球化学指标，建立地球化学预测模式，对已知矿体及外围异常进行了地球化学评价。     

判定地下水二维水动力弥散参数的直线图解法

弥散参数的求得，是正确认识地下水污染运移及对其控制的前提。本文根据极值理论和非线性方程线性化手段，导出了浓度和时间变化的线性关系，得出了确定弥散参数的直线图解法。它可以弥补其它方法的不足。借用实际弥散试验资料对其进行了检验，表明所建数学模型与实际模型之间吻合良好，即所求得的参数结果可信。