OPTIMIZATION OF SHORT DIGITAL LINEAR FILTERS FOR INCREASED ACCURACY*

The accuracy of short length digital linear filter operators can be substantially increased if the sampling interval as well as the abscissa shift are properly adjusted. This may be done by a trial and error process of adjustment of these parameters until the error made by the filter operator, applied to a suitably chosen test function, is smallest. As an illustration of the application of this method, 7-, 11- and 19-point filters for the calculation of Schlumberger apparent resistivity from a known resistivity transform are designed. Errors with the new 7-point filter are seen to be less than those with a 19-point filter of conventional design. The errors with the new 19-point filter are two to three orders of magnitude smaller than those made by the conventional 19-point filter. The new method should provide digital linear operators that allow significant improvements in accuracy for comparable computation efforts, or substantial reduction in computation for comparable accuracy of results, or something of both.     

Hankel transform filters for dipole antenna radiation in a conductive medium

We discuss the Hankel transforms related to a particular application, i.e. the dipole antenna radiation in conductive media, such as the antenna radiation in sea-bed electromagnetic applications. In this application, the electromagnetic wavefields decay very rapidly with distance. A good filter means that it can be used to evaluate weak fields. Exponential sampling transforms a Hankel transform into a convolution equation, which must be solved to obtain the filter coefficients. Here, we use a direct matrix inversion method to solve the convolution equation in the sample domain, instead of the Fourier transform method and the Wiener–Hopf method, previously used to solve the convolution equation. This direct method is conceptually simple and is suitable for our optimization process: by using the Sommerfeld identity, we search for the optimum sampling interval, which corresponds to the minimum wavefield, evaluated for a given length filter. The performances of the new filters obtained are compared with some well-known filters. We find that our filters perform better for our application; that is, for the same length filters, our filters are able to calculate weaker fields. For users working in similar applications, three sets of filters with lengths 61, 121, 241 are available from the author.     

OPTIMIZED FAST HANKEL TRANSFORM FILTERS1

In the linear digital filter theory for calculation of Hankel transforms it is possible to find explicit series expansions for the filter coefficients. A method is presented for optimizing the Hankel filters calculated in this way. For a certain desired accuracy of computation, the sampling density and filter length are minimized by choosing the parameters determining the filter characteristics according to the analytical properties of the input function. A new approach to the calculation of the filter coefficients has been developed for these optimized filters. The length of the filters may be further reduced by introducing a shift in the sampling scheme.     

The Green￠s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green's functions ofhorizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous con-ditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green's functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.     

半空间饱和土在内部简谐水平力作用下的Ge函数

根据积分变换方法得出了半空间内部作用简谐水平力时的Gree函数．首先，利用Hankel积分变换方法，直接对频域内的Biot波动方程进行求解，得出Biot波动方程的通解；利用通解和半空间内部作用水平力时边界上的边界条件，以及力作用面上的连续性条件，可以得出上述边值问题的解；对于边值问题在变换域内的解进行相应的逆变换，就可以得出频域内的Gree函数．本文得到的线弹性退化解与文献中的结果吻合．最后，文中给出了两个算例．      

The Green''''s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

IntroductionThe wave propagation problems in saturated soil are very important for the civil engineering, geophysics and seismology. Biot (1956,1962) established the theory of wave propagation in saturated soil firstly, and hereafter many researchers have used Biot theory to study wave propagation problems in saturated soil. By using integral transform and potential function method, Philippacopoulos (1988) studied the Lamb(s problem of a vertical point force applied to the surface of saturate…     

The Green’s function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green’s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above-mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green’s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper. Foundation item: State Natural Science Foundation (59879012) and Doctoral Foundation from State Education Commission (98024832).     

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.     

Transient foundation solution of saturated soil to impulsive concentrated loading

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.     

Torsional vibration of a finite cylindrical cavity in a two-layer transversely isotropic half-space

The aim of this paper is to present a rigorous investigation for a two-layered transversely isotropic linear elastic half-space containing a circular cylindrical cavity of length equal to the top layer undergoing mono-harmonic ring shape shear stress applied either on the vertical cylindrical surface or on the base of the cavity. To this end, a combination of Fourier cosine integral transform for depth and Hankel integral transform for radial distance are used, which translate the boundary value problem to a singular integral equation for the shear stress comes out from the continuity of two layers. The integral equation is solved for some collocation points with a smoothed variable of distance, which is adapted with the use of a free parameter. It is shown that, although the shear stress is highly singular, it does not highly depend on this free parameter. Both the analytical and numerical results are verified with both the static isotropic and dynamic transversely isotropic homogeneous cases. In addition, some new graphical results are presented for more understanding in engineering point of view.     

水平层状介质视电阻率的高精度计算公式及其误差分析

为适应地电阻率法预报地震深入研究的需要,本文介绍了几种高精度计算水平层状介质视电阻率的公式及方法。对于MN→0的视电阻率计算,本文给出了J_1(x)滤波系数法的计算步骤,此外对二层、三层介质的传统级数求和法,给出了估算精度的公式;对于MN≠O的视电阻率计算,本文讨论了数值积分法和J_0(x)滤波系数法的计算步骤以及为提高精度所采取的技术措施,并且给出了一批用两种方法所得到的计算结果。     

GRAVITY INTERPRETATION USING THE HANKEL TRANSFORM*

The Hankel transform theorem can be applied to the inversion of gravity data for the buried sphere, the horizontal cylinder, and the vertical rod. This new approach leads to exact solutions of the transforms for the assumed bodies. A comparison with the classical procedure by Fourier transform reveals that for the sphere and the vertical rod, the Hankel transform is preferable.     

OPTIMUM DIGITAL FILTERS FOR SIGNAL TO NOISE RATIO ENHANCEMENT*

One of the main objectives of seismic digital processing is the improvement of the signal-to-noise ratio in the recorded data. Wiener filters have been successfully applied in this capacity, but alternate filtering devices also merit our attention. Two such systems are the matched filter and the output energy filter. The former is better known to geophysicists as the crosscorrelation filter, and has seen widespread use for the processing of vibratory source data, while the latter is. much less familiar in seismic work. The matched filter is designed such that ideally the presence of a given signal is indicated by a single large deflection in the output. The output energy filter ideally reveals the presence of such a signal by producing a longer burst of energy in the time interval where the signal occurs. The received seismic trace is assumed to be an additive mixture of signal and noise. The shape of the signal must be known in order to design the matched filter, but only the autocorrelation function of this signal need be known to obtain the output energy filter. The derivation of these filters differs according to whether the noise is white or colored. In the former case the noise autocorrelation function consists of only a single spike at lag zero, while in the latter the shape of this noise autocorrelation function is arbitrary. We propose a novel version of the matched filter. Its memory function is given by the minimum-delay wavelet whose autocorrelation function is computed from selected gates of an actual seismic trace. For this reason explicit knowledge of the signal shape is not required for its design; nevertheless, its performance level is not much below that achievable with ordinary matched filters. We call this new filter the “mini-matched” filter. With digital computation in mind, the design criteria are formulated and optimized with time as a discrete variable. We illustrate the techniques with simple numerical examples, and discuss many of the interesting properties that these filters exhibit.     

ON THE CONSTRUCTION OF THE 3D BAND-LIMITED EXTRAPOLATION OPERATOR IN THE SPACE-FREQUENCY DOMAIN1

Based on an expansion of the band-limited 3D extrapolation operator in terms of orthogonal Chebyshev polynomials, a closed form expression of the space-frequency response is presented. A key step is an evaluation of the (inverse) 2D Fourier transform of circularly symmetric functions, which is related to the (zero-order) Hankel transform. Hankel transforms of individual members of the orthogonal set of polynomials are available from tables and summation of series; hence, the real and the imaginary parts of the space-frequency response can be found in terms of cylindrical and spherical Bessel functions, respectively. The procedure permits an efficient and accurate evaluation of the space-frequency response.     

The 3-D Truncation Filtering Methodology Defined for Planar and Spherical Models: Interpreting Gravity Data Generated by Point Masses

This paper focuses on one particular way of linear filtering the gravity data to facilitate gravity inversion or interpretation. With the use of integral transforms the gravity anomalies are transformed into new quantities that allow an easier interpretation with the help of pattern recognition. As the integral transforms are in fact filters, and as the regions of integration are caps with a variable radius, which can be systematically changed as a free parameter, we refer to such methodology as the truncation filtering. Such filters may be understood as weighted spherical windows moving over the surface, on which the gravity anomaly is defined, the kernel of the transform being the weight function. The objective of this paper is to define and deploy the truncation filtering for a planar model, i.e. for a homogenous horizontally infinite layer with embedded anomalous masses, and for a spherical model, i.e., for a homogenous massive sphere with embedded anomalous masses. Instead of the original gravity anomaly, the quantities resulting from the truncation filtering are interpreted/inverted. As we shall see, this approach has certain benefits. The fundamental concept of the truncation filtering methodology is demonstrated here in terms of the model consisting of one point mass anomaly.The relationship between the depth of the point mass and the instant of the onset of the dimple pattern observed in sequences produced by truncation filtering the synthetic gravity data generated by point masses is, for both the planar and spherical models, compiled by computer simulations, as well as derived analytically. It is shown, that the dimple pattern is a consequence of truncating the domain of the filter and is free of the choice of the kernel of the filter. It is shown, that in terms of the mean earth and depths of point masses no greater than some 100 km the spherical model may be replaced by a planar model from the perspective of the truncation filtering methodology. It is also shown, that from the viewpoint of the truncation filtering methodology the rigorous gravity anomaly may be approximated by the vertical component of the gravity disturbance. The relationship between the instant of the dimple onset and the depth of the point mass thus becomes linear and independent of the magnitude (mass) of the point mass.     

层状介质中水平磁偶极子电磁场空间分布特征

针对岩层中介电常数和磁导率随深度变化的情况,给出水平磁偶极子源电磁场的水平或垂直空间分布特征及相应的变化规律。采用Kong给出的汉克尔（Hankel）J0变换线性滤波器（241点）和汉克尔J1变换线性滤波器（241点）算法,选用均匀大地和两层地电模型,讨论介质的介电常数和磁导率变化时的水平磁偶极子电磁场空间分布特征。计算结果表明,磁场的峨分量幅值与磁导率反相关;磁场H2和电场Ey分量随着磁导率的增大而增大;低频电磁测深中介电常数变化时对电磁场基本没有影响。     

An approach to handling non-Gaussianity of parameters and state variables in ensemble Kalman filtering

The ensemble Kalman filter (EnKF) is a commonly used real-time data assimilation algorithm in various disciplines. Here, the EnKF is applied, in a hydrogeological context, to condition log-conductivity realizations on log-conductivity and transient piezometric head data. In this case, the state vector is made up of log-conductivities and piezometric heads over a discretized aquifer domain, the forecast model is a groundwater flow numerical model, and the transient piezometric head data are sequentially assimilated to update the state vector. It is well known that all Kalman filters perform optimally for linear forecast models and a multiGaussian-distributed state vector. Of the different Kalman filters, the EnKF provides a robust solution to address non-linearities; however, it does not handle well non-Gaussian state-vector distributions. In the standard EnKF, as time passes and more state observations are assimilated, the distributions become closer to Gaussian, even if the initial ones are clearly non-Gaussian. A new method is proposed that transforms the original state vector into a new vector that is univariate Gaussian at all times. Back transforming the vector after the filtering ensures that the initial non-Gaussian univariate distributions of the state-vector components are preserved throughout. The proposed method is based in normal-score transforming each variable for all locations and all time steps. This new method, termed the normal-score ensemble Kalman filter (NS-EnKF), is demonstrated in a synthetic bimodal aquifer resembling a fluvial deposit, and it is compared to the standard EnKF. The proposed method performs better than the standard EnKF in all aspects analyzed (log-conductivity characterization and flow and transport predictions).     

THE USE OF FILTERED BESSEL FUNCTIONS IN DIRECT INTERPRETATION OF GEOELECTRICAL SOUNDINGS *

We start from the Hankel transform of Stefanescu's integral written in the convolutionintegral form suggested by Ghosh (1971). In this way it is possible to obtain the kernel function by the linear electric filter theory. Ghosh worked out the sets of filter coefficients in frequency domain and showed the very low content of high frequencies of apparent resistivity curves. Vertical soundings in the field measure a series of apparent resistivity values at a constant increment Δx of the logarithm of electrode spacing. Without loss of information we obtain the filter coefficient series by digital convolution of the Bessel function of exponential argument with sine function of the appropriate argument. With a series of forty-one values we obtain the kernel functions from the resistivity curves to an accuracy of better than 0.5%. With the digital method it is possible to calculate easily the filter coefficients for any electrode arrangement and any cut-off frequency.     

电阻率测深的数字解释

本文主要介绍了应用积分变换的方法和采样定理将视电阻率ρ_s曲线作线性滤波,得出一新的电阻率转换函数T′曲线,然后,以层参数（各层的电阻率和厚度）算得的T用最优化数值方法在DJS-6型电子计算机上与其进行自动拟合,以达到解释电测深曲线的目的。 文中简述了戈什（Ghosh）提出的对ρ_s作线性滤波的原理,介绍了与国外不同的取样间距和滤波系数的确定以及阻尼最小二乘法和变尺度最优化法的计算框图和应用,最后附有实例和简要的讨论。