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.     

变质反应的代数识别方法

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

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

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

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

Regionalization of unit hydrograph Parameters: 2. Uncertainty analysis

Hydrologic model parameters obtained from regional regression equations are subject to uncertainty. Consequently, hydrologic model outputs based on the stochastic parameters are random. This paper presents a systematic analysis of uncertainty associated with the two parameters, N and K, in Nash's IUH model from different regional regression equations. The uncertainty features associated with N and K are further incorporated to assess the uncertainty of the resulting IUH. Numerical results indicate that uncertainty of N and K from the regional regression equations are too significant to be ignored.     

Regionalization of unit hydrograph Parameters: 2. Uncertainty analysis

Hydrologic model parameters obtained from regional regression equations are subject to uncertainty. Consequently, hydrologic model outputs based on the stochastic parameters are random. This paper presents a systematic analysis of uncertainty associated with the two parameters, N and K, in Nash's IUH model from different regional regression equations. The uncertainty features associated with N and K are further incorporated to assess the uncertainty of the resulting IUH. Numerical results indicate that uncertainty of N and K from the regional regression equations are too significant to be ignored.     

电法勘探的发展和展望

电法勘探在经历了近一个世纪的发展后，其方法理论、仪器设备、野外数据采集、处理和解释等方面都经历了一系列重大变化．本文以方法理论的进展为主线，回顾、展望了目前电法勘探中几个重要而令人关注的研究焦点．这些问题的研究进展将会对21世纪的电法勘探产生深远的影响．