首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Koefoed has given practical procedures of obtaining the layer parameters directly from the apparent resistivity sounding measurements by using the raised kernel function H(λ) as the intermediate step. However, it is felt that the first step of his method—namely the derivation of the H curve from the apparent resistivity curve—is relatively lengthy. In this paper a method is proposed of determining the resistivity transform T(λ), a function directly related to H(λ), from the resistivity field curve. It is shown that the apparent resistivity and the resistivity transform functions are linearily related to each other such that the principle of linear electric filter theory could be applied to obtain the latter from the former. Separate sets of filter coefficients have been worked out for the Schlumberger and the Wenner form of field procedures. The practical process of deriving the T curve simply amounts to running a weighted average of the sampled apparent resistivity field data with the pre-determined coefficients. The whole process could be graphically performed within an quarter of an hour with an accuracy of about 2%.  相似文献   

2.
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.
It has been found that the Wiener-Hopf least-squares method is a very successful tool for the determination of resistivity sounding filters. The values of the individual filter coefficients differ quite appreciably from those obtained by the Ghosh procedure. These differences in the filter coefficients, however, have only a negligible effect on the output of the filter. It seems that these differences in the coefficients correspond to a filter function of a rather narrow frequency band around the Nyquist frequency, which is only very weakly present in the input and output functions.  相似文献   

4.
Different sets of filter coefficients for the linear filter technique for the computations of resistivity and EM sounding curves are evaluated for several electrode and coil configurations. Instead of this procedure, the two-electrode filter can be used for computations of Wenner, Schlumberger, and dipole—dipole apparent resistivity model curves by defining convolutional expressions which contain the new input functions in terms of the resistivity transform function. Similarly, the Schlumberger filter performs the computations of dipole—dipole apparent resistivity model curves. The Wenner, Schlumberger, and dipole—dipole filter functions are defined in terms of the two-electrode filter using the new convolutional expressions. A relationship between the Schlumberger and dipole—dipole filter functions is given. The above arguments are adopted for the computations of EM sounding curves. It is shown that the EM filter for the horizontal coplanar loop system (which is identical to the two-electrode filter) performs the computations of the mutual coupling ratios for perpendicular, vertical coplanar, and vertical coaxial loop systems. In the same way, the Schlumberger filter can be used to compute vertical coaxial sounding curves. The corresponding input functions are defined in terms of the EM kernel for all convolutional expressions presented. After these considerations, integral expressions of the mutual coupling ratios involving zero-order Bessel function are derived. The mutual coupling ratio for the vertical coaxial loop system is given in the same form as the mutual coupling ratio for the vertical coplanar loop system.  相似文献   

5.
A simple unified equation of apparent resistivity for a general four-electrode array is developed. The main idea is the analytical integration of the Stefanescu expression for potential over a layered earth by writing an exponential approximation for the kernel function. Finally a matrix equation is developed to estimate the kernel function from observed apparent resistivity values. The general equation automatically reduces to the particular configuration once the electrode separations are modified suitably. Examples for Schlumberger and Wenner configurations are numerically calculated to estimate the precision of the method. Good results in a short execution time are obtained, irrespective of the shape of the apparent resistivity curve. Finally, the full interpretation of one theoretical resistivity curve and two field resistivity curves is demonstrated. The more stable ridge-regression estimation method is used in the identification of layer parameters from the kernel function.  相似文献   

6.
This paper describes certain procedures for deriving from the apparent resistivity data as measured by the Wenner electrode configuration two functions, known as the kernel and the associated kernel respectively, both of which are functions dependent on the layer resistivities and thicknesses. It is shown that the solution of the integral equation for the Wenner electrode configuration leads directly to the associated kernel, from which an integral expression expressing the kernel explicitly in terms of the apparent resistivity function can be derived. The kernel is related to the associated kernel by a simple functional equation where K1(λ) is the kernel and B1(λ) the associated kernel. Composite numerical quadrature formulas and also integration formulas based on partial approximation of the integrand by a parabolic arc within a small interval are developed for the calculation of the kernel and the associated kernel from apparent resistivity data. Both techniques of integration require knowledge of the values of the apparent resistivity function at points lying between the input data points. It is shown that such unknown values of the apparent resistivity function can satisfactorily be obtained by interpolation using the least-squares method. The least-squares method involves the approximation of the observed set of apparent resistivity data by orthogonal polynomials generated by Forsythe's method (Forsythe 1956). Values of the kernel and of the associated kernel obtained by numerical integration compare favourably with the corresponding theoretical values of these functions.  相似文献   

7.
The technique of linear digital filtering developed for the computation of standard curves for conventional resistivity and electromagnetic depth soundings is applied to the determination of filter coefficients for the computation of dipole curves from the resistivity transform function by convolution. In designing the filter function from which the coefficients are derived, a sampling interval shorter than the one used in the earlier work on resistivity sounding is found to be necessary. The performance of the filter sets is tested and found to be highly accurate. The method is also simple and very fast in application.  相似文献   

8.
In the linear filter method of interpreting resistivity sounding data, as developed by Ghosh (1971), it appears that the filter function in the x-domain approaches an oscillating function for both large positive and large negative abscissa values. In the present note the reason for this oscillating behaviour is derived, and a possible practical application is indicated.  相似文献   

9.
In this paper a fast method is developed for computing apparent resistivity curves for known layer configurations. The method is based on the application of a linear filter to determine the apparent resistivity curve from, the kernel function.  相似文献   

10.
The values of the filter coefficients used for the computation of electromagnetic sounding curves are studied in conjunction with the values of the input function to the filter, or the range of values which the input function may assume, and the filters are broken off at such a place that the error in the sum of the products of filter coefficient and input function does not exceed a prescribed value. This analysis is carried out for the horizontal coils system, the perpendicular coils system, and the vertical coplanar coils system. The lengths of the filters so derived depend on the layer parameters, the frequency and the coil spacing. Even in the most unfavourable cases the filters are shorter than the filters published by Koefoed, Ghosh, and Polman (1972).  相似文献   

11.
In a previous paper by the present author a method was developed for direct interpretation of resistivity observations made with a Schlumberger electrode configuration. This method consisted of two steps. The first of these was to derive the kernel function in the integral expression for the apparent resistivity from the observed data; the second step was to derive the resistivity stratification from this kernel function. The first of these two steps depends on the electrode configuration that has been used. In the present paper the above mentioned method is modified so as to make it apply to a Wenner electrode configuration. The procedure is indicated by which the method may be adapted to any other electrode configuration in which the distances between the electrodes are finite. The second step in the interpretation, i.e. the derivation of the resistivity stratification from the kernel function, is independent of the electrode configuration used, and therefore needs no further discussion in the present context.  相似文献   

12.
The induced polarization response of a horizontally multilayered earth with no resistivity contrast can rapidly be calculated on a desk calculator or minicomputer for any electrode array. The formulation is a simple series summation of the products of weighting coefficients and the true induced polarization responses for each of the layers. The coefficients are directly derivable from the corresponding resistivity model. This series approach to IP formulation was originally described by Seigel but has not been treated extensively in the present-day geophysical literature. This method can be applied to either time or frequency domain induced polarization measurements. Once the coefficients are known, apparent induced polarization response can readily be obtained by judicious substitution of known, suspected, or assumed values of the true induced polarization of each layer. Basic formulation is presented for the IP potential coefficients (pole-pole or two array) with no resistivity contrast between the layers. From these coefficients, response of any number of layers for any electrode array can be obtained by suitable differentiation. Some examples of Wenner array for a three-layered earth and dipole-dipole array for a four-layered earth are used to illustrate the application. The results of this technique are valid for many natural situations of modest resistivity contrast. However, they definitely cannot be used if there are highly contrasting resistivity layers present. Such an approach is conceptually simple and is useful for survey planning, checking or setting the “depth-of-penetration”of a given array. For field induced polarization data that fits reasonably well to the no-resistivity-contrast model, this simple approach facilitates quantitative interpretation.  相似文献   

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

14.
A mathematical analysis is given of the phenomenon of equivalence in resistivity sounding, which is based upon the properties of the raised kernel function. Analysis of this function instead of the apparent resistivity function is justified because, as has been shown in a previous publication, variations in the apparent resistivity function lead to variations in the raised kernel function with relative values of the same order of magnitude The expression for the raised kernel function is expanded into a Mac Laurin series. Equivalence can occur only if the second order term of this series is negligible. The coefficient of the first order term depends on the resistivity and the thickness of the layer under consideration. There is an infinite set of combinations of values for these two quantities, for which the coefficient of the first order term has the same value. All these combinations represent equivalent layer distributions.  相似文献   

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

16.
本文采用数字滤波技术,研究了电阻率测井的理论曲线计算和Ra曲线的环境影响校正,并对影响数字滤波效果的因素进行了分析。结果表明,用数字滤波新技术进行环境校正,提高了校正质量,比传统的图版校正方法有明显优点。文中采用数字滤波方法绘制的Ra理论曲线,证明所推导的数学模型是正确的。  相似文献   

17.
The computation method described in this paper is based on the existence of a linear relationship between the mutual coupling ratio and the kernel function in the integral expression for it. Accordingly, the mutual coupling ratio can be determined by first computing sample values of the kernel function and then subjecting these to a digital linear filter. In the present paper the appropriate sampling distance is determined and the values of the digital filter coefficients are computed, both for electromagnetic sounding with horizontal coils and for electromagnetic sounding with perpendicular coils.  相似文献   

18.
This paper describes the procedure for interpreting the apparent resistivity data measured with the two-electrode array directly with the help of kernel function. The calculation of kernel function from the observed resistivity curve is done by the method of decomposition. In the method of decomposition the resistivity curve is approximated by a sum of certain functions, whose choice is only restricted by the requirement that the contribution to the kernel function corresponding to them should be easily computable. A few such functions are classified. These, and the standard curves for corresponding kernel functions obtained by utilising an integral expression for two-electrode array expressing the kernel explicitly in terms of the apparent resistivity functions, are plotted on log-log scale. The determination of layer parameters, that is, the layer resistivities and thicknesses from the kernel function can be carried out by a method proposed by Pekeris (1940).  相似文献   

19.
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.  相似文献   

20.
It was found in Part I of this paper that approximating the sharp cut-off frequency characteristic best in a mean square sense by an impulse response of finite length M produced a characteristic whose slope on a linear frequency scale was proportional to the length of impulse response, but whose maximum overshoot of ±9% was independent of this length (Gibbs' phenomenon). Weighting functions, based on frequency tapering or arbitrarily chosen, were used in Part II to modify the truncated impulse response of the sharp cut-off frequency characteristic, and thereby obtain a trade-off between the value of maximum overshoot and the sharpness of the resulting characteristic. These weighting functions, known as apodising functions, were dependent on the time-bandwidth product , where , corresponded to the tapering range of frequencies. Part III now deals with digital filters where the number 2N–1 of coefficients is directly related to the finite length M of the continuous impulse response. The values of the filter coefficients are taken from the continuous impulse response at the sampling instants, and the resulting characteristic is approximately the same as that derived in Part II for the continuous finite length impulse response. Corresponding to known types of frequency tapering, we now specify a filter characteristic which is undefined in the tapering range, and determine the filter coefficients according to a mean square criterion over the rest of the frequency spectrum. The resulting characteristic is dependent on the time bandwidth product = (N–1/2)ξ up to a maximum value of 2, beyond which undesirable effects occur. This optimum partially specified characteristic is an improvement on the previous digital filters in terms of the trade-off ratio for values of maximum overshoot less than 1%. Similar to the previous optimum characteristic is the optimum partially specified weighted digital filter, where greater “emphasis is placed on reducing the value of maximum overshoot than of maximum undershoot”. Such characteristics are capable of providing better trade-off ratios than the other filters for maximum overshoots greater than 1/2%. However these filters have critical maximum numbers 2.NC–1 of coefficients, beyond which the resulting characteristics have unsuitable shapes. This type of characteristic differs from the others in not being a biassed odd function about its cut-off frequency.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号