首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 312 毫秒
1.
The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data. Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined. The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.  相似文献   

2.
The Nguyen and Pinder method is one of four techniques commonly used for analysis of response data from slug tests. Limited field research has raised questions about the reliability of the parameter estimates obtained with this method. A theoretical evaluation of this technique reveals that errors were made in the derivation of the analytical solution upon which the technique is based. Simulation and field examples show that the errors result in parameter estimates that can differ from actual values by orders of magnitude. These findings indicate that the Nguyen and Pinder method should no longer be a tool in the repertoire of the field hydrogeologist. If data from a slug test performed in a partially penetrating well in a confined aquifer need to be analyzed, recent work has shown that the Hvorslev method is the best alternative among the commonly used techniques.  相似文献   

3.
2D and 3D potential-field upward continuation using splines   总被引:1,自引:0,他引:1  
The dominant upward‐continuation technique used in the potential‐field geophysics industry is the fast Fourier transform (FFT) technique. However, the spline‐based upward‐continuation technique presented in this paper has some advantages over the FFT technique. The spline technique can be used to carry out level‐to‐uneven surface 2D and 3D potential‐field upward continuation. An example of level‐to‐uneven surface upward continuation of 3D magnetic data using the spline technique is shown, and it is evident that the continued anomalies are very close to the theoretical values. The spacing can be irregular. Synthetic examples using the spline technique to continue noise‐contaminated gravity and magnetic data upward to an altitude of 15 km on irregular grids are shown. Gaussian noise with a zero mean and a standard deviation of 1% does not cause much error and can readily be tolerated. Through comparison with the FFT technique, it is found that for low‐altitude gravity and magnetic upward continuation, both the FFT technique and the spline technique are suitable; for high‐altitude upward continuation, the FFT technique is inaccurate, whereas the spline technique works very well. Also, upward continuation by the spline technique has a smaller edge effect than upward continuation by the FFT technique. The spline‐based upward continuation technique works fairly well even when the periphery of a grid is not quiet: it is rather robust in general. A real example shows that the spline technique can be employed to perform upward continuation of total‐field magnetic data and to de‐emphasize near‐surface noise.  相似文献   

4.
用最大似然法进行波场分解和震相识别   总被引:6,自引:0,他引:6       下载免费PDF全文
本文以介质中弹性波场的传播特征和多元统计分析为基础,从数学上进行严格的推导,得到了模型与资料拟合的最大似然判据,以检验地震波列中是否存在某种特殊波型的能量。利用概率滤波可得分解后的地震图,从图中可以直接得到震相到时。用此方法处理P波资料,效果较好。对于剪切波资料,有时处理结果不十分理想。这可能是由于径向分量和横向分量之间的线性相关性造成的。对于因介质各向异性而形成的分裂剪切波,用此方法得到的S波到时为快S波到时。此方法对于单台三分向记录资料的处理是十分有效的。  相似文献   

5.
魏宝君  LIU Q H 《地球物理学报》2007,50(5):1595-1605
为提高水平层状介质中三维电磁波散射和逆散射数值模拟的效率,在对角张量近似(DTA)的基础上根据不同回代方式得到了求解积分方程的DTA1和DTA2两种近似. 这两种近似可以作为计算积分方程稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法的初始猜测值和预条件因子,从而形成效率更高的混合DTA-BCGS算法. 散射实例说明了DTA2的高精度和混合DTA-BCGS算法尤其是混合DTA2-BCGS算法的高效率. 由于DTA2近似程度更高,将DTA2与变型Born迭代反演方法(DBIM)相结合形成了一种对三维异常体进行重构的快速电磁波逆散射技术. 文中的逆散射实例说明所开发的逆散射技术对重构水平层状介质中的任意三维异常体是非常有效的.  相似文献   

6.
Matrix equations are derived to transform the resistivity sounding data obtained in one type of a four-electrode array to the corresponding resistivity sounding data that would be obtained using a different four-electrode array. These expressions are based primarily on recent work in which we have established a linear relation between the apparent resistivity and the kernel function by using a powerful exponential approximation for the kernel function. It is shown that the resistivity sounding data of two different four-electrode arrays have a linear relation through an essentially non-singular matrix operator and, as such, one is derivable from the other for a one-dimensional model and it can also be extended to two-dimensions. Some numerical examples considering synthetic data are presented which demonstrates the efficiency of the method in such transformations. Two published field examples are also considered for transformation giving a reliable interpretation.  相似文献   

7.
The magnetic and gravity field produced by a given homogeneous source are related through Poisson's equation. Starting from this consideration, it is shown that some 2D interpretation tools, widely applied in the analysis of aeromagnetic data, can also be used for the interpretation of gravity gradiometric data (vertical gradient). This paper deals specifically with the Werner deconvolution, analytic signal and Euler's equation methods. After a short outline of the mathematical development, synthesized examples have been used to discuss the efficiency and limits of these interpretation methods. These tools could be applied directly to airborne gravity gradiometric data as well as ground gravity surveys after transformation of the Bouguer anomalies into vertical gradient anomalies. An example is given of the application of the Werner deconvolution and Euler's equation methods to a microgravity survey.  相似文献   

8.
An equivalent linearization technique to obtain the response of non-linear multi-degree-of-freedom dynamic systems to stationary gaussian excitations is developed. The non-linearities are assumed to be single-valued functions of accelerations, velocities and displacements. Using a property of gaussian vector processes, the closed forms of the coefficients of the equivalent linear system are obtained by the direct application of partial differentiation and expectation operators to the non-linear terms. It is shown that when the non-linearities possess potentials, the linear system has symmetric coefficient matrices. A geometrical interpretation of the linear coefficients, in connection with the original non-linearities, is presented. The accuracy is investigated by means of examples.  相似文献   

9.
Variables controlling the reliability of interpolated estimates of rainfall are identified and the significance of each is quantified. Alternative ways are investigated of expressing the gauge network and rainfall distribution patterns in a form that may be used in estimating equations for assessing error in interpolated values. A scheme is devised for obtaining equations analytically for any interpolation method used to derive information and the results of a study relevant to a specific analysis technique are shown.The purpose of the scheme is to enable confidence limits to be placed during computer processing of rainfall data on interpolated output. Comparisons may then be made between available analysis techniques, results being realisticaly assessed under operational conditions. Once a user has selected a suitable method of analysis, the factors likely to control the accuracy of the information produced are monitored during processing and likely tolerance limits are given for output values.  相似文献   

10.
本文基于地震剖面的数字化图像特征,提出用光栅算子将地震剖面按产状分离成产状相对单一的多幅数字图像剖面。对每一幅产状相对单一的数字图像剖面,用向量分解技术的向量最大特征值对应基向量投影的原理获取有效相干信息的办法实现信噪分离;在结合地质、物探分析的基础上进行最终输出剖面的滤波加权重建,获得信噪比高、波组特征保持好、分辨率相对保证的剖面。实际资料处理结果表明本文提出方法的去噪应用效果显著。  相似文献   

11.
In the year 1958 the Service for Water Management of the ?Rijkswaterstaat” started its program of geo-electrical resistivity prospecting in the western part of the Netherlands. The aim of this program was to obtain data on the salinity distribution of the ground water. The ground water regime in this part of the Netherlands is most intricate. This is due to the geological and geo-hydrological conditions and to the low elevation of the land. Many reclaimed areas are up to several metres below mean sea level. The resistivity data obtained are closely related to the salinity of the ground water. On the basis of bore hole data it was even possible to arrive at calibration curves for the salinity of the ground water in sand deposits. Under special conditions it was also possible to draw conclusions with respect to the presence of less permeable formations as e.g. clay layers. Some remarks are given on the practical performance and the interpretation of the measurements. A review is given of the work done until now. Some results are shown by means of maps of the salinity distribution of the ground water in the western and northern parts of the Netherlands. Two examples are described of the use of the data obtained during the survey in the province of Zuid-Holland. Another two examples are presented of detailed investigations for special purposes in relatively small areas.  相似文献   

12.
The standard method of surface wave tomography assumes the use of data on surface wave velocities measured on paths wholly located inside the region being investigated and oriented in different directions. If the data on group velocities obtained in the regional network of stations are used for this purpose, then the earthquake sources also should be situated within the limits of this region. If the region is small, then this requirement restricts the range of periods of surface waves and, correspondingly, the depth of research. The use of data from distant earthquakes enlarges the range of periods. However, it does not enable one to use the usual tomography method because of an inappropriate configuration of the system of paths: outside of the region in which the stations are located, the paths do not intersect, but, at the same time, the paths from a particular earthquake source toward the different stations of networks are very close over a large distance, and the mean correction to the velocity in the path sections outside of the network of stations can be accepted as identical. This assumption forms the basis of the proposed modification of the method of surface wave tomography, in which group velocities measured at stations of the local network from distant earthquakes are used. The lateral variations in the velocity are determined within the limits of the network, and the mean corrections to the velocity are determined on the paths from different earthquake sources in the sections outside of the network. As earlier, the condition of smoothness is imposed on the distribution of velocity variations, and the condition of smallness of the sum of their squares is imposed on the values of the corrections to the velocity. Testing of this method based on model examples and on real data obtained in the network of stations in Tibet showed its advantage over the standard tomography method.  相似文献   

13.
A method is presented for quantifying the uncertainty of the semivariogram of transmissivity and determining the required number of measurements. In this method, the estimated semivariogram and its 95% confidence limits are first determined from a finite number of measurements. The uncertainty of the estimated semivariogram is then quantified using the random field simulation technique. For a given value of the quantitative index of uncertainty, the required number of measured data can finally be obtained. Actual transmissivity data of an existing groundwater monitoring network are used in the application of the proposed method. The required numbers of measurements of transmissivity for four different values of the quantitative index of uncertainty are provided, from which reliable semivariograms of the transmissivity can be obtained. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

14.
This paper addresses the issue of system identification for linear structural systems using earthquake induced time histories of the structural response. The proposed methodology is based on the Eigensystem Realization Algorithm (ERA) and on the Observer/Kalman filter IDentification (OKID) approach to perform identification of structural systems using general input–output data via Markov parameters. The efficiency of the proposed technique is shown by numerical examples for the case of eight-storey building finite element models subjected to earthquake excitation and by the analysis of the data from the dynamic response of the Vincent-Thomas cable suspension bridge (Long Beach, CA) recorded during the Whittier and the Northridge earthquakes. The effects of noise in the measurements and of inadequate instrumentation are investigated. It is shown that the identified models show excellent agreement with the real systems in predicting the structural response time histories when subjected to earthquake-induced ground motion. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

15.
The main results of deep seismic sounding (DSS) are usually presented in the form of high-velocity models of the medium. Some model examples and the international DOBRE profile have shown that the informativeness of the data obtained can be significantly enhanced by the construction of wave images of the Earth’s crust, based on the migration of refracted and wide-angle reflected waves. The Donets Basin Refraction/Reflection Experiment (DOBRE) profile crosses the Dnieper-Donets paleorift zone in the Donbas region. Along the profile, refracted waves from the basement and the upper mantle and the reflections from the crust basement (the M boundary) are reliably traced. This wave migration has been used to construct a wave image of the structure of the Earth’s crust. As a result, a clear seismic image of the basement surface, whose depth changes along the profile from 0 to 20 km, was obtained. In near-slope parts of the basin, several major faults were identified that had not been identified previously during standard kinematic data processing. It is shown that the crust-upper mantle transition zone is a clearly reflective horizon only within the crystalline massifs; under a depression, it is represented by a lens-shaped highly-heterogeneous area. As shown in the model examples, the images obtained using such a migration accurately reflect the structural features of the medium, in spite of its complicated structure.  相似文献   

16.
Seismic data processing is a challenging task, especially when dealing with vector-valued datasets. These data are characterized by correlated components, where different levels of uncorrelated random noise corrupt each one of the components. Mitigating such noise while preserving the signal of interest is a primary goal in the seismic-processing workflow. The frequency-space deconvolution is a well-known linear prediction technique, which is commonly used for random noise suppression. This paper represents vector-field seismic data through quaternion arrays and shows how to mitigate random noise by proposing the extension of the frequency-space deconvolution to its hypercomplex version, the quaternion frequency-space deconvolution. It also shows how a widely linear prediction model exploits the correlation between data components of improper signals. The widely linear scheme, named widely-linear quaternion frequency-space deconvolution, produces longer prediction filters, which have enhanced signal preservation capabilities shown through synthetic and field vector-valued data examples.  相似文献   

17.
Five examples, obtained during exploration for hydrocarbons in the Pannonian Basin of Hungary, are used to show how the interpretation of seismic sections can be usefully complemented by results from MT surveys. Selection of the most appropriate MT quantities, considered to be proper ‘MT attributes’ for the purpose of visualization as well as recognition of the subsurface structures and the different inversions of MT data is essential for practical integration of seismic and MT surveys. A new technique providing a semiquantitive MT-attribute pseudosection for the purpose of visualization of the subsurface structures is proposed. The procedure utilizes derivative functions of the phase of MT impedance for visualization and derives estimated depths from the Bostick transformation of Cagniard apparent resistivities. On the basis of the MT-attribute pseudosections, constructed from the phase derivatives and transformed resistivity data, depths are estimated for interfaces between geological formations with significant resistivity contrast. In particular examples, the interface between the Tertiary sediments and the older basement rocks as well as tectonic fracture zones with decreased resistivity can be resolved.  相似文献   

18.
The paper presents 3D diffraction imaging based on the spectral decomposition of a different combination of selective or partial images. These images are obtained by the pre-stack asymmetric migration procedure, which is weighted data summation. Spectral decomposition is done in the Fourier domain with respect to spatial dip and azimuth angles. Numerical examples with the application of different workflows for the synthetic and real data examples demonstrate detailed reliable reconstruction of the fractured zones and reliable reconstruction of fracture orientation on synthetic and real 3D data examples.  相似文献   

19.
A numerical technique to compute the resistivity transform directly from the observed Wenner sounding data has been developed. In principle, the procedure is based on a decomposition method and consists of two steps: the first step determines a function that approximates the apparent resistivity data and the second step transforms this function into the corresponding kernel by an analytical operation. The proposed method is tested on some theoretical master curves. A high degree of precision is achieved with very little computer time. The applicability is shown on two field examples.  相似文献   

20.
Variations of seismic interval velocities within the cable length cause anomalies in the stacking velocity analyses. Utilizing the approximation of rectilinear ray propagation, i.e. supposing that the velocity changes cause time delays only, it is shown that the stacking velocity anomalies are linearly related to the interval velocity variations. In particular, the stacking velocity anomaly is calculated when the interval velocity of an intermediate layer undergoes a stepwise variation. The amplitude of the anomaly increases with the ratio between horizon depth and cable length. From the forward model, a program for the inversion is derived in order to identify lateral changes of interval velocities from unsmoothed stacking velocity analyses. Some examples of the application of this technique to synthetic and real data are presented.  相似文献   

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

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