RECIPROCAL AVERAGING TECHNIQUES IN THE GEOELECTRICAL BOUNDARY ELEMENT APPROACH*

A Fredholm integral equation of the first kind with respect to the surface charge density is presented via a weighted residual formulation in the standard isotropic problem. The surface charge densities are numerically obtained by solving the well-known Fredholm integral equation of the second kind, where the new equation can be regarded as a constraint. The accuracy of boundary element solutions is examined in connection with the violation of such a constraint and the “modified reciprocal averaging techniques” are proposed.     

New Methods for Modeling Laterolog Resistivity Corrections

The paper presents methods for laterolog response modeling. In Coulomb’s charges method, Laplace’s equation is solved for the electric field distribution in rock medium with internal boundaries between different resistivity layers. There, the boundary problem is reduced to Fred-holm integral equation of the second kind. The second method uses a finite element array to model apparent resistivity from laterolog. The task is treated as DC problem and the Laplace equation is solved numerically. The presented methods were applied to borehole data covering a typical stratigraphie section of the Fore-Sudetic Monocline in southwestern Poland. Apparent resistivity was calculated using the Coulomb’s charges method and alternatively modeled using a finite element method which gave similar results. Then, a series of linear corrections for borehole, shoulder bed, and filtration effects for apparent resistivity obtained by the Coulomb’s charges method demonstrated the feasibility of calculating true resistivity of virgin and invaded zones. The proposed methods provide a flexible solution in modeling which can be adapted to other logs.     

埋地目标体矢量电磁散射的一种快速正演算法

利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。     

Constraining of the zero total surface charge in galvanic modelling

The galvanic problem is frequently solved by a Fredholm integral equation of the second kind based on a single layer source formulation. At higher conductivity contrasts between the model and its surroundings the homogeneous part of the integral equation approaches an eigenvalue equation. With infinite contrast the solution of this limiting integral equation is non-unique, but in the subspace of zero total charge the solution is unique. This mathematical property of the integral equation is reflected in its numerical solution with the result that large numerical errors may appear and convergence of the solution becomes very slow. Errors are, for the most part, related to the computed excess charge generated in the numerical solution. The effect is studied by comparing the results computed from the solution of the integral equation alone with those computed from a particular solution where the requirement of zero total charge is used as a constraint. The model examples clearly show that the use of the constraint condition significantly improves the accuracy of the results.     

Approximate Solutions to the Boundary–Volume Integral Equation for Wave Propagation in Piecewise Heterogeneous Media

Piecewise heterogeneous media that the earth presents are composed of large-scale boundary structures and small-scale volume heterogeneities. Wave propagation in such piecewise heterogeneous media can be accurately superposed through the generalized Lippmann–Schwinger integral equation (GLSIE). Two different Born series modeling schemes are formulated for the boundary–volume integral equation with 2-D antiplane motion (SH waves). Both schemes decompose the resulting boundary–volume integral equation matrix into two parts: the self-interaction operator handled with a fully implicit manner, and the extrapolation operator approximated by a Born series. The first scheme associates the self-interaction operator with each boundary itself and the volume itself, and interprets the extrapolation operator as the cross-interaction between each boundary and other boundaries/volume scatterers in a subregion. The second scheme relates the self-interaction operator to each boundary itself and its cross-interaction with the volume scatterers on both sides, and expresses the extrapolation operator as both the direct and indirect (through the volume scatterers) cross-interactions between different boundaries in a subregion. By eliminating the displacement field from the volume scatterers, the second scheme reduces the dimension of the resulting boundary-volume integral equation matrix, leading to a faster convergence than the first scheme. Both the numerical schemes are validated by dimensionless frequency responses to a heterogeneous alluvial valley with the velocity perturbed randomly in the range of ca 5–20 %. The schemes are applied to wave propagation simulation in a heterogeneous multilayered model by calculating synthetic seismograms. Numerical experiments, compared with the full-waveform numerical solution, indicate that the Born series modeling schemes significantly improve computational efficiency, especially for high frequencies.     

APPROXIMATING RAINFALL-RUNOFF MODELLING RESPONSE USING A STOCHASTIC INTEGRAL EQUATION

Rainfall-runoff modelling uncertainty can be analysed by the use of a stochastic integral formulation. The stochastic integral equation can be based on the rainfall–runoff model input of model rainfall or model rainfall excess. Similarly, the stochastic integral equation can be based on the rainfall–runoff model output of the modelled runoff hydrograph. The residual between actual measured runoff data and modelled runoff (from the rainfall–runoff model) is analysed here by the use of a stochastic integral equation. This approach is used to develop a set of convolution integral transfer function realizations that represent the chosen rainfall–runoff modelling error. The resulting stochastic integral component is a distribution of possible residual outcomes that may be directly added to the rainfall–runoff model's deterministic outcome, to develop a distribution of probable runoff hydrograph realizations from the chosen rainfall–runoff model.     

电阻率测井响应的积分方程解法

从电流守恒条件出发,得到了两个积分方程.其中第一个和在1964年提出的相同,它的未知函数是边界面上的隐电流源分布密度;第二个积分方程的未知函数则是电位函数本身. 在的文章以及本文作者的另一篇文章中,供电电极和测量电极都被当成点电极处理.用这种方法来处理电位电极系和梯度电极系已不够准确,用来处理三侧向或双侧向等具有大电极的电极系则完全无能为力. 测井中所用的电极系是包着金属外皮的绝缘心棒,金属外皮就是电极.本文提出了处理这种电极系的方法.结果表明,第二个积分方程在处理这种电极系时有明显优点.     

曲面上的位场转换（二维和三维情况）

在文献[1]中我们利用复变函数论的方法只解决了二维的曲面上的位场转换问题.在本文中,我们利用单层位势理论将它化为一个奇异积分方程求解的问题,并且提出了用奇点分离法来解此奇异积分方程,这对二维和三维情况都适用.通过理论模型的计算,结果的精度是令人满意的.     

电法勘探正演数值模拟的若干结果

本文前一部分叙述了点源二维电阻率法数值模拟正演计算的有限单元法,文中采用混合边界条件、用LL~T分解解线性方程组等优化措施,使二维有限单元法电阻率法正演计算的速度和精度都比目前国外流行的L.Rijo方法和程序有了进一步的提高。文中介绍了有关的方法并引入了若干计算结果。本文后一部分叙述了三维电阻率法数值模拟的积分方程近似解法,由于对K.Dieter等人提出的方法作了一些近似处理,并用迭代法求解积分方程表示式,实算结果说明在保证计算精度的条件下,提高了计算速度,从而提高了三维电场正演计算的有效性和实用性。所提出的方法容易推广到激发极化法的正演计算中。     

Torsional response of a rigid circular foundation on a saturated half-space to SH waves

An analytical approach is used to study the torsional vibrations of a rigid circular foundation resting on saturated soil to obliquely incident SH waves. Biot’s poroelastic dynamic theory is considered to characterize the saturated soil below the foundation, which is solved by Hankel transform later. In order to consider the scattering phenomena caused by the existence of the foundation, the total wave field in soil is classified into free-field, rigid-body scattering field and radiation scattering field. According to the classification of wave field and the mixed boundary-value conditions between the soil and the foundation, torsional vibrations of the foundation are formulated in two sets of dual integral equations. Then, the dual integral equations are reduced to Fredholm integral equation of the second kind to be solved. Combining with the dynamic equilibrium equations of the foundation, the expressions for the torsional vibrations of the foundation are obtained. Numerical results are presented to demonstrate the influence of excitation frequency, incident angle, the torsional inertia moment of the foundation and permeability of the saturated half-space on the torsional vibrations of the foundation.     

From Satellite Gradiometry Data to Subcrustal Stress Due to Mantle Convection

Subcrustal stress induced by mantle convection can be determined by the Earth’s gravitational potential. In this study, the spherical harmonic expansion of the simplified Navier–Stokes equation is developed further so satellite gradiometry data (SGD) can be used to determine the subcrustal stress. To do so, we present two methods for producing the stress components or an equivalent function thereof, the so-called S function, from which the stress components can be computed numerically. First, some integral estimators are presented to integrate the SGD and deliver the stress components and/or the S function. Second, integral equations are constructed for inversion of the SGD to the aforementioned quantities. The kernel functions of the integrals of both approaches are plotted and interpreted. The behaviour of the integral kernels is dependent on the signal and noise spectra in the first approach whilst it does not depend on extra information in the second method. It is shown that recovering the stress from the vertical–vertical gradients, using the integral estimators presented, is suitable, but when using the integral equations the vertical–vertical gradients are recommended for recovering the S function and the vertical–horizontal gradients for the stress components. This study is theoretical and numerical results using synthetic or real data are not given.     

Application of the Lorenz equations for studying thermal convection in the lithosphere

Stationary solutions including wave solutions with constant amplitudes are found for nonlinear equations of thermal convection in a layer with nonlinear rheology. The solution is based on the Fourier expansion of unknown velocities and temperatures with only the first and first two terms retained in the velocity and temperature series, respectively. This method, which can be regarded as the Lorenz method, yields the Lorenz equations that fairly well describe the thermal convection in a layer with Newtonian rheology if the Rayleigh number is not very large. The obtained generalization of the Lorenz equations to the case of an integral (having a memory) nonlinear rheology implies that only the first term is retained in the Fourier series for the stress components, i.e., the nonlinear rheological equation is harmonically linearized. However, in the Fourier series of temperature, it is essential to keep the second term: this term, which is independent of the horizontal coordinate, models the thermal boundary layer that characterizes the developed convection. We constructed the bifurcation curves that describe the stationary convection in the nonlinear integral medium simulating the rheology of the mantle, and analyzed the stability of stationary convective flows. The Lorenz method is applied to study small-scale thermal convection in the lithosphere of the Earth.     

有能量吸收的一维波动方程的阻抗反演

本文讨论了一维粘弹性波阻抗反演的迭代法。根据一维波动方程和ＫｅｌｖｉｎＶｏｉｇｔ粘弹模型导出了一个适于数值计算的第二类Ｖｏｌｔｅｒｒａ型积分方程，用迭代法求解。通过对已有方法的改进和完善，数值计算表明了该反演方法具有一定的精度、抗干扰能力和分辨率以及应用于实际水平层状地层反射地震记录反演阻抗的可能性．     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

POTENTIAL FIELD OF A STATIONARY ELECTRIC CURRENT USING FREDHOLM'S INTEGRAL EQUATIONS OF THE SECOND KIND*

An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite.     

一种新的三维大地电磁积分方程正演方法

采用规则六面体单元和并矢Green函数奇异积分等效积分技术,已有的大地电磁积分正演方法具有不能有效模拟地下复杂地质体和计算精度偏低的缺点.本文提出了一种新的三维大地电磁积分方程正演技术,即采用四面体单元、解析的并矢Green函数奇异积分表达式,达到既能模拟地下复杂异常体,又能有效提高已有积分方程法计算精度的目的.首先,采用四面体网格技术离散地下复杂异常体,获得四面体单元上的大地电磁积分方程.然后,利用针对四面体单元开发的新的奇异值积分的解析表达式,准确计算线性方程中的并矢Green函数的奇异积分,从而获得精确的线性方程.借助于PARDISO高性能并行直接求解器,实现了三维大地电磁问题的高精度求解.最后,基于国际标准3D-1模型和六棱柱模型,通过与其他方法结果的对比分析,验证了本文方法的正确性、处理高电导率对比度的能力(1000:1)和处理复杂模型的能力.     

A Neumann expansion approach to flow through heterogeneous formations

A stochastic approach is used for the study of flow through highly heterogeneous aquifers. The mathematical model is represented by a random partial differential equation in which the permeability and the porosity are considered to be random functions of position, defined by the average value, constant standard deviation and autocorrelation function characterized by the integral scale. The Laplace transform of the solution of the random partial differential equation is first written as a solution of a stochastic integral equation. This integral equation is solved using a Neumann series expansion. Conditions of convergence of this series are investigated and compared with the convergence of the perturbation series. For mean square convergence, the Neumann expansion method may converge for a larger range of variability in permeability and porosity than the classic perturbation method. Formal expressions for the average and for the correlation moments of the pressure are obtained. The influence of the variability of the permeability and porosity on pressure is analyzed for radial flow. The solutions presented for the pressure at the well, as function of the permeability coefficient of variation, may be of practical interest for evaluating the efficiency of well stimulation operations, such as hydraulic fracturing or acidizing methods, aimed at increasing the permeability around the well.     

三维地形频率域人工源电磁场的边界元模拟方法

提出了一种用边界元法计算频率域人工源三维地形电磁场的数值模拟方法.首先用矢量积分理论和电磁场边界条件，将上半空间（空气）和下半空间（地下介质）两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中，假设单元中电磁场为无限大气空间电磁场与地形影响的叠加,并假设地形影响为常项,这样既保证了计算精度又使得计算方法简便.通过分解和计算，每一个矢量面积分方程分解为对应三个坐标方向的三个常量线性方程，这些线性方程组成了对角占优的线性方程组，可用SSOR方法求解.文中给出了垂直磁偶源的垂直磁场地形影响的例子.     

Horizontal vibrations of a disk on a poroelastic half-space

This paper analytically examines the horizontal vibration of a rigid disk on a saturated poroelastic half-space. The pressure-solid displacement form of the harmonic equations of motion for asymmetric dynamic problem are developed from the form of the equations originally presented by Biot. Making use of a new method the solution of the above equations is obtained. According to the mixed boundary -value conditions, the dual integral equations of the horizontal vibration of a rigid disk on a saturated poroelastic half-space are established. By appropriate transforms, it is shown that the dual integral equations can be reduced to a pair of Fredholm integral equations of the second kind, whose solutions are then computed. Numerical results for the horizontal dynamic compliance coefficient are given at the end of this paper.     

Conditionality of inverse solutions to discretised integral equations in geoid modelling from local gravity data

The eigenvalue decomposition technique is used for analysis of conditionality of two alternative solutions for a determination of the geoid from local gravity data. The first solution is based on the standard two-step approach utilising the inverse of the Abel-Poisson integral equation (downward continuation) and consequently the Stokes/Hotine integration (gravity inversion). The second solution is based on a single integral that combines the downward continuation and the gravity inversion in one integral equation. Extreme eigenvalues and corresponding condition numbers of matrix operators are investigated to compare the stability of inverse problems of the above-mentioned computational models. To preserve a dominantly diagonal structure of the matrices for inverse solutions, the horizontal positions of the parameterised solution on the geoid and of data points are identical. The numerical experiments using real data reveal that the direct gravity inversion is numerically more stable than the downward continuation procedure in the two-step approach.