Dynamic response of a multi-layered poroelastic medium

An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed.     

Efficiency of trenches along railways for trains moving at sub- or supersonic speeds

A full 3D analytical approach is adopted to account for trenches on one or both sides close to a railroad. Low-frequency ground vibrations are investigated due to the passing of trains, and open trenches are used as wave barriers. The modelling technique is based on Fourier transforms and Fourier series. The ground is modelled as a layered semi-infinite domain and the embankment with finite layers. The trenches are obtained by simulating the upper surface layer with two or three finite rectangular regions with appropriate widths. A particular boundary condition is adopted at the vertical sides of all finite regions to enable the solution procedure. Rails and sleepers are accounted for with Euler–Bernoulli beams and an anisotropic Kirchhoff plate with transversal isotropy. The wheel loads from the boogie wheel pairs of the train are simulated as moving forces. Hence, no irregularities in rails or wheels are accounted for.     

Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays

In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result.     

各向异性层状介质中视电阻率与磁场响应研究

针对任意各向异性地层,利用极向型和环向型标量位函数,导出相应的直流电视电阻率和磁电阻率的磁场响应关系.计算了各向异性地层的直流电视电阻率和磁电阻率响应特征,重点分析了电阻率测深方法与磁电阻率探测方法对地下各向异性介质的探测能力.文中采用状态矩阵的分析方法,首先采用极向型和环向型标量位构造了各向异性层状介质电场与磁场的通解,利用各层界面电场、磁场的连续性及地面激励源的耦合条件,推导了不同层之间电磁场的状态矩阵,建立了空间电场与磁场的递归计算关系.其次,针对递归计算中指数项数值计算的不稳定性,借用状态矩阵的性质,导出了将不稳定指数计算项转化为稳定的指数项的转换关系.针对横向各向同性(TI)介质中极向位与环向位解耦的特点,导出了电磁场的直接积分解.最后,采用解析解验证了算法的正确性,给出了多层各向异性地层模型的视电阻率和磁场响应曲线,分析了直流电法探测裂缝性地层、估计裂缝分布性状的可能性.     

Solution of forward 2-D geoelectric problems using integrodifferential equations

The solution of forward 2-D geoelectric problems is considered for anisotropic and bianisotropic media. The problems are reduced to the solution of integrodifferential equations obtained with the use of integral Fourier and Fourier-Bessel transforms for horizontally heterogeneous layered media and radially heterogeneous layered media. This approach is particularly effective for studying the influence of small-scale geological structures (such as faults, borehole wall mud cakes, and the like) on results of geoelectric studies.     

TRANSFORMATION OF RESISTIVITY SOUNDING MEASUREMENTS OBTAINED IN ONE ELECTRODE CONFIGURATION TO ANOTHER CONFIGURATION BY MEANS OF DIGITAL LINEAR FILTERING*

The technique of digital linear filtering is used for transformation of apparent resistivity data from one electrode configuration into another. Usually filter spectra are determined via the discrete Fourier transforms of input and output functions: the filter characteristic is the quotient of the spectra of the output function and input function. In this paper, the transformation of the apparent resistivities is presented for four electrode configurations (Wenner, the two-electrode, Schlumberger, and dipole configurations). In our method, there is no need to use the discrete Fourier transform of the input and output functions in order to determine the filter spectrum for converting apparent resistivity in one electrode configuration to any other configuration. Sine responses for determination of the derivative of apparent resistivities are given in analytical form. If the filter spectrum for converting the apparent resistivity to the resistivity transform for one electrode configuration is known, the filter spectra for transforming the apparent resistivity to the resistivity transform for any electrode configurations can be calculated by using newly derived expressions.     

井下地电阻率观测影响系数分析——以江宁地震台为例

采用水平层状均匀介质中点电源位于任意深度时的电位解析表达式,以江宁台3层电性结构为例,分析了井下对称四极地电阻率观测时各层影响系数随深度、极距的变化,并结合探测深度探讨了实施井下观测时影响系数在选择供电极距和电极埋深时的作用。结果表明,对于"K"型电性结构,江宁台井下观测对地表、浅层干扰有较强的抑制作用,其短极距观测对地表、浅层干扰的抑制能力显著优于长极距观测;长极距观测在电极埋深H小于100m时对地表介质季节性的干扰具有放大作用;浅层影响系数一定时,电极埋深和供电极距需同时增加;江宁台井下观测供电极距AB/2取100～150m、电极埋深H为250m较为合理。     

重力异常场空间-波数混合域三维数值模拟

重力勘探中复杂条件下的三维正演计算量大存储要求高,使得这种条件下重力勘探高效、精细正反演变得困难.针对这一问题,提出一种空间-波数混合域数值模拟方法,该方法将空间域引力位积分进行水平方向二维傅里叶变换,将三维空间域卷积问题转换为多个不同波数之间相互独立的空间垂向一维积分问题,一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征密度变化,可得出单元积分的解析表达式.该方法计算量和存储需求少,算法高度并行;保留垂向为空间域,优势之一在于可根据实际情况合理调整单元疏密程度,准确模拟任意复杂地形和密度异常体的重力异常,兼顾计算精度与计算效率;优势之二在于用形函数拟合求得积分的解析解,计算精度和效率高;充分利用一维形函数积分的高效和高精度,不同波数之间一维积分高度并行性及快速傅里叶变换的高效性,实现重力异常场三维数值模拟.设计棱柱体模型,通过数值解和解析解对比验证了该方法的正确性、适用性和高效性.针对任意复杂地形条件下的重力场及其张量的模拟问题,提出一种快速算法,对其有效性进行了验证.探究标准FFT法的截断效应对计算精度的影响,对比分析Gauss-FFT法和标准FFT扩边法两种方法的计算精度和效率,总结了二者的选取策略,结果表明选用标准FFT扩边法计算效率更高.实际地形的数值模拟表明本文算法适用于任意复杂地形的高效计算.     

A FAST FOURIER TRANSFORM METHOD FOR RAPID COMPUTATION OF GRAVITY AND MAGNETIC ANOMALIES DUE TO ARBITRARY BODIES*

The spectrum of a magnetic or a gravity anomaly due to a body of a given shape with either homogeneous magnetization or uniform density distribution can be expressed as a product of the Fourier transforms of the source geometry and the Green's function. The transform of the source geometry for any irregularly-shaped body can be accurately determined by representing the body as closely as possible by a number of prismatic bodies. The Green's function is not dependent upon the source geometry. So the analytical expression for its transform remains the same for all causative bodies. It is, therefore, not difficult to obtain the spectrum of an anomaly by multiplying the transform of the source geometry by that of the Green's function. Then the inverse of this spectrum, which yields the anomaly in the space domain, is calculated by using the Fast Fourier Transform algorithm. Many examples show the reliability and accuracy of the method for calculating potential field anomalies.     

3D modelling and sensitivity in DC resistivity using charge density

A three‐dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite‐difference codes for two models: the two‐layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole‐pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression.     

Cross‐hole resistivity tomography using different electrode configurations

This paper investigates the relative merits and effectiveness of cross‐hole resistivity tomography using different electrode configurations for four popular electrode arrays: pole–pole, pole–bipole, bipole–pole and bipole–bipole. By examination of two synthetic models (a dipping conductive strip and a dislocated fault), it is shown that besides the popular pole–pole array, some specified three‐ and four‐electrode configurations, such as pole–bipole AM–N, bipole–pole AM–B and bipole–bipole AM–BN with their multispacing cross‐hole profiling and scanning surveys, are useful for cross‐hole resistivity tomography. These configurations, compared with the pole–pole array, may reduce or eliminate the effect of remote electrodes (systematic error) and yield satisfactory images with 20% noise‐contaminated data. It is also shown that the configurations which have either both current electrodes or both potential electrodes in the same borehole, i.e. pole–bipole A–MN, bipole–pole AB–M and bipole–bipole AB–MN, have a singularity problem in data acquisition, namely low readings of the potential or potential difference in cross‐hole surveying, so that the data are easily obscured by background noise and yield images inferior to those from other configurations.     

Pilot‐scale field validation of the long electrode electrical resistivity tomography method

A validation experiment, carried out in a scaled field setting, was attempted for the long electrode electrical resistivity tomography method in order to demonstrate the performance of the technique in imaging a simple buried target. The experiment was an approximately 1/17 scale mock‐up of a region encompassing a buried nuclear waste tank on the Hanford site. The target of focus was constructed by manually forming a simulated plume within the vadose zone using a tank waste simulant. The long electrode results were compared to results from conventional point electrodes on the surface and buried within the survey domain. Using a pole‐pole array, both point and long electrode imaging techniques identified the lateral extents of the pre‐formed plume with reasonable fidelity but the long electrode method was handicapped in reconstructing vertical boundaries. The pole‐dipole and dipole‐dipole arrays were also tested with the long electrode method and were shown to have the least favourable target properties, including the position of the reconstructed plume relative to the known plume and the intensity of false positive targets. The poor performance of the pole‐dipole and dipole‐dipole arrays was attributed to an inexhaustive and non‐optimal coverage of data at key electrodes, as well as an increased noise for electrode combinations with high geometric factors. However, when comparing the model resolution matrix among the different acquisition strategies, the pole‐dipole and dipole‐dipole arrays using long electrodes were shown to have significantly higher average and maximum values within the matrix than any pole‐pole array. The model resolution describes how well the inversion model resolves the subsurface. Given the model resolution performance of the pole‐dipole and dipole‐dipole arrays, it may be worth investing in tools to understand the optimum subset of randomly distributed electrode pairs to produce maximum performance from the inversion model.     

Near surface geophysical surveys with a high frequency mutual impedance measuring system

Relationships for the mutual impedance of horizontal coplanar, vertical coplanar, perpendicular and vertical coaxial loop antenna arrays operating at a high frequency, located over a layered conducting medium, have been derived. The relationships take into account the occurrence of displacement currents in both the free space and the conducting medium. The relations can be used for calculating the earth parameters in geophysical prospecting. The use of high frequency antenna array makes it possible to determine distributions of apparent resistivity as well as apparent electrical permittivity. The potential application of the measurements of mutual impedance of loop antennas operating at a high frequency to the geophysical examinations is presented. The measuring system developed for this purpose is described. A simple method of interpreting the mutual impedance measurements to obtain the traces of apparent resistivity and apparent permittivity of the investigated earth is given.     

The electromagnetic response of a horizontal electric dipole buried in a multi‐layered earth

The electromagnetic response of a horizontal electric dipole transmitter in the presence of a conductive, layered earth is important in a number of geophysical applications, ranging from controlled‐source audio‐frequency magnetotellurics to borehole geophysics to marine electromagnetics. The problem has been thoroughly studied for more than a century, starting from a dipole resting on the surface of a half‐space and subsequently advancing all the way to a transmitter buried within a stack of anisotropic layers. The solution is still relevant today. For example, it is useful for one‐dimensional modelling and interpretation, as well as to provide background fields for two‐ and three‐dimensional modelling methods such as integral equation or primary–secondary field formulations. This tutorial borrows elements from the many texts and papers on the topic and combines them into what we believe is a helpful guide to performing layered earth electromagnetic field calculations. It is not intended to replace any of the existing work on the subject. However, we have found that this combination of elements is particularly effective in teaching electromagnetic theory and providing a basis for algorithmic development. Readers will be able to calculate electric and magnetic fields at any point in or above the earth, produced by a transmitter at any location. As an illustrative example, we calculate the fields of a dipole buried in a multi‐layered anisotropic earth to demonstrate how the theory that developed in this tutorial can be implemented in practice; we then use the example to examine the diffusion of volume charge density within anisotropic media—a rarely visualised process. The algorithm is internally validated by comparing the response of many thin layers with alternating high and low conductivity values to the theoretically equivalent (yet algorithmically simpler) anisotropic solution, as well as externally validated against an independent algorithm.     

二维层状介质点电源边界单元解

本文用付氏变换将点源二维地电断面的三维边值问题变成二维边值问题,然后讨论了边界元法在解决层状介质中存在二维柱体时点源场视电阻率计算问题中的应用。由于选取的基本解避免了对水平层状界面的积分计算,从而大大减少了边界剖分节点数。在计算中设法略去了对边界上位场法向导数的求解,使未知量大为减少,从而便利了求解二维边值问题,并由付氏反变换求得三维空间中的位。本方法所需计算机内存少,原始数据的输入简单,且有较高的计算精度。     

The effect of the electrical anisotropy on the response of helicopter-borne frequency-domain electromagnetic systems

Helicopter electromagnetic (HEM) systems are commonly used for conductivity mapping and the data are often interpreted using an isotropic horizontally layered earth model. However, in regions with distinct dipping stratification, it is useful to extend the model to a layered earth with general anisotropy by assigning each layer a symmetrical 3 × 3 resistivity tensor. The electromagnetic (EM) field is represented by two scalar potentials, which describe the poloidal and toroidal parts of the magnetic field. Via a 2D Fourier transform, we obtain two coupled ordinary differential equations in the vertical coordinate. To stabilize the numerical calculation, the wavenumber domain is divided into two parts associated with small and large wavenumbers. The EM field for small wavenumbers is continued from layer to layer with the continuity conditions. For large wavenumbers, the EM field behaves like a DC field and therefore cannot be sensed by airborne EM systems. Thus, the contribution from the large wavenumbers is simply ignored. The magnetic fields are calculated for the vertical coaxial (VCX), horizontal coplanar (HCP) and vertical coplanar (VCP) coil configurations for a helicopter EM system. The apparent resistivities defined from the VCX, VCP and HCP coil responses, when plotted in polar coordinates, clearly identify the principal anisotropic axes of an anisotropic earth. The field example from the Edwards Aquifer recharge area in Texas confirms that the polar plots of the apparent resistivities identify the principal anisotropic axes that coincide well with the direction of the underground structures.     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

A DIGITAL LINEAR FILTER FOR RESISTIVITY SOUNDING WITH A GENERALIZED ELECTRODE ARRAY*

This paper presents a digital linear filter which maps composite resistivity transforms to apparent resistivities for any four—electrode array over a horizontally layered earth. A filter is provided for each of three sampling rates; the choice of filter will depend on resistivity contrasts and computational facilities. Two methods of filter design are compared. The Wiener-Hopf least-squares method is preferable for low sampling rate filters. The Fourier transform method is more successful in producing a filter with a high sampling rate which can handle resistivity contrasts of 100 000: 1.     

任意各向异性介质中三维可控源音频大地电磁正演模拟

在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向.