首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 380 毫秒
1.
Airborne VLF data are routinely collected by The Geological Survey of Sweden (SGU) as part of its bedrock mapping programme. In this paper we demonstrate that the novel Tensor VLF technique developed at Uppsala University and SGU can provide useful qualitative and quantitative information about the electrical conductivity distribution in the upper few hundred meters. Single transmitter scalar VLF maps emphasize those conductive structures that have dominant strikes in the direction of the transmitter. The tensor tipper (essentially the vertical magnetic field from currents along the strike direction) calculated from multiple transmitters is dependent only upon the underlying conductivity structure. Transformation of the tipper into the peaker (the horizontal divergence) has proven to enhance the lateral resolution while the transformation to the apparent resistivity can be used to discriminate different rock types. Two case histories from the application of VLF data are presented in this study. Two dimensional structures can be quantitatively modelled by modern inversion methods developed originally for deep electromagnetic MT soundings. Direct inversion of the real and imaginary parts of the tipper provides more quantitative information about the subsurface resistivity distribution.  相似文献   

2.
The interpretation of airborne VLF data represents an important aspect of geophysical mapping of the upper few hundred meters of the Earth's crust, especially in areas with crystalline rocks. We have examined the ability of the single frequency VLF method to provide quantitative subsurface resistivity information using two generic models and standard airborne parameters with a flight altitude of 70 m and a frequency of 16 kHz. The models are long thin conductor (10 m thick, 10 Ω m resistivity and 1 km long) and a wider buried conductive dike (100 Ω m resistivity and 500 m wide). Using standard regularized inversion it turned out that for both models the conductivity of the conductors are underestimated and the vertical resolution is rather poor. The lateral positions of the minimum of the resistivity distributions coincide well with the true positions of the shallow conductors. For deeper conductors the position of the minimum resistivity moves from the edges of the conductor into the conductor. The depth to the minimum of the resistivity anomalies correlates well with the true depth to the top of the conductors although the latter is always smaller than the former.Interpretation of field airborne data collected at 70 m flight height resolved both small scale and large scale near surface conductors (conductance ∼1 S). Deeper conductors show up in the VLF data as very long wavelength anomalies that are particularly powerful in delineating the lateral boundaries of the conductors. Many of the VLF anomalies in the Stockholm area are dominated by these deep conductor responses with some near surface conductors superimposed. The deep conductors often follow topographic lows coinciding with metasediments. We interpret the frequent absence of near surface responses at 70 m flight height as a result of weak coupling between the primary VLF wave and the small scale (in all three dimensions) near-surface conductors.Radio magnetotelluric (RMT) ground measurements were carried out along a short profile coinciding with part of an airborne profile. Using data at 9 frequencies (14–250 kHz) small scale conductors in the upper few tens of meters, not identified from the airborne data, could be well resolved. Large scale deeper conductors could be identified by both methods at nearly the same positions.  相似文献   

3.
罗鸣  李予国  李刚 《地球物理学报》2016,59(11):4349-4359
本文提出了一维垂直各向异性(VTI)介质倾斜偶极源频率域海洋可控源电磁(CSEM)资料高斯-牛顿反演方法.在电阻率各向异性介质水平偶极源和垂直偶极源海洋CSEM正演算法的基础上,利用欧拉旋转方法,实现了各向异性介质倾斜偶极源海洋CSEM正演算法.海洋可控源电磁场关于地下介质横向电阻率(ρ_h)和垂向电阻率(ρ_v)的偏导数(即灵敏度矩阵)是解析计算的,结合垂直各向异性介质横向电阻率与垂向电阻率的关系,将各向异性率融入到正则化因子选择中,实现了正则化因子的自适应选择.理论模型合成数据和实测资料反演算例表明,我们提出的反演方法能够较准确的重构海底围岩和基岩的各向异性电阻率以及高阻薄层的埋藏深度、厚度和垂向电阻率.  相似文献   

4.
The applicability of three kinds of electrode configurations used to delineate a buried horizontal pipe was studied. A 3D resistivity imaging survey was carried out along eight parallel lines using pole-pole, pole-dipole, and dipole-dipole arrays with 1m minimum electrode spacings. Roll-along measurements were carried out to cover a rectangular grid. The 2D and 3D least squares algorithms based on the robust inversion method were used in the inversion of the apparent resistivity data sets. The 2D inversion of data sets could not delineate the orientation and dimension of the subsurface anomalies clearly. To obtain more accurate results, a 3D joint inversion of the pole-pole and pole-dipole data sets was performed, as well as of pole-pole and dipole-dipole data sets. In this case, both horizontal and vertical dimensions of subsurface structures were resolved. The resulting model obtained from each array was compared to those of joint inversion method. The result showed that the horizontal resolution does not improve so much as that in the vertical direction when joint inversion is applied.  相似文献   

5.
The study presents a fast imaging technique for the very low‐frequency data interpretation. First, an analytical expression was derived to compute the vertical component of the magnetic field at any point on the Earth's surface for a given current density distribution in a rectangular block on the subsurface. Current density is considered as exponentially decreasing with depth, according to the skin depth rule in a particular block. Subsequently, the vertical component of the magnetic field due to the entire subsurface was computed as the sum of the vertical component of the magnetic field due to an individual block. Since the vertical component of the magnetic field is proportional to the real part of very low‐frequency anomaly, an inversion program was developed for imaging of the subsurface conductors using the real very low‐frequency anomaly in terms of apparent current density distribution in the subsurface. Imaging results from the presented formulation were compared with other imaging techniques in terms of apparent current density and resistivity distribution using a standard numerical forward modelling and inversion technique. Efficacy of the developed approach was demonstrated for the interpretation of synthetic and field very low‐frequency data. The presented imaging technique shows improvement with respect to the filtering approaches in depicting subsurface conductors. Further, results obtained using the presented approach are closer to the results of rigorous resistivity inversion. Since the presented approach uses only the real anomaly, which is not sensitive to very small isolated near‐surface conducting features, it depicts prominent conducting features in the subsurface.  相似文献   

6.
The electromagnetic array profiling (EMAP) exploration method can be combined with a direct one-dimensional inversion process for conversion to depth to produce a subsurface resistivity cross-section. This cross-section may then be interpreted in parallel with a seismic cross-section to enhance the prediction of rock type and structure. In complex thrust environments and areas of shallow carbonate rocks, the EMAP method is often used to provide additional data either to help the seismic processor and/or to aid the seismic interpretation. In particular, the electromagnetic (EM) data can be used to build an independent seismic velocity file for depth migration. Three EMAP test areas in the western United States are used to demonstrate such a use of EMAP as an expioration tool. The first shows how a velocity file is estimated from resistivity data for seismic depth migration processing in a complex thrust environment. In the second example, the method is applied in layer-cake geology with high seismic velocity rocks at the earth's surface. The third example is another complex thrust environment, but in this case the velocity file derived from the resistivity data is used for stacking the seismic data.  相似文献   

7.
Integrated geophysical surveys using vertical electrical sounding (VES), very low frequency (VLF) EM, radiation counting, total magnetic field and self-potential (SP) measurements are carried out to characterize the geothermal area around a hot spring in the Nayagarh district, Orissa, India that lies in the East Indian geothermal province. The study was performed to delineate the fracture pattern, contaminated groundwater movement and possible heating source. VES interpretations suggest a three- to four-layer structure in the area. Resistivity survey near the hot spring suggests that weathered and fractured formations constitute the main aquifer system and extend to 60 m depth. Current flow measured at various electrode separations normalized by the applied voltage suggests that fractures extend to a greater depth. Detailed VLF study shows that fractures extend beyond 70 m depth. VLF anomaly has also very good correlation with the total magnetic field measured along the same profiles. Study results suggest that a gridded pattern of VLF survey could map the underground conductive fracture zones that can identify the movement of contaminated groundwater flow. Therefore, precautionary measures can be taken to check further contamination by delineating subsurface conducting structures. Self potential (SP) measured over the hot spring does not show a large anomaly in favor of the presence of a sulphide mineral body. A small positive (5–15mV) SP anomaly is measured which may be streaming potential due to subsurface fluid flow. A high radiation is measured about four kilometers from the hot spring, suggesting possible radiogenic heating. However, the exact nature of the heating source and its depth is not known in the area. Deep resistivity followed by a magneto-telluric survey could reveal the deeper structures.  相似文献   

8.
An analysis of the equivalence phenomenon, analogous to that encountered in resistivity sounding, has been made for electromagnetic (frequency) sounding with the following systems: horizontal coplanar coils, perpendicular coils, vertical coplanar coils, and vertical coaxial coils. Cases analyzed are three-layer H and K type resistivity distributions. The theoretical responses for the EM sounding systems have been computed by the digital linear filter method using short filters. An analysis has been made of the resolution by the four EM systems considered f the equivalence effect. It is concluded that from the equivalence point of view electromagnetic (frequency) sounding has relatively better resolution compared to resistivity sounding for the H type three-layer cases. For the K type three-layer cases electromagnetic (frequency) sounding provides a serious problem in resolution similar to that encountered in resistivity sounding. Empirical relationships have been established for the H and K type cases respectively in order to numerically quantify the equivalence involved in electromagnetic (frequency) sounding.  相似文献   

9.
In order to interpret field data from small-loop electromagnetic (EM) instruments with fixed source–receiver separation, 1D inversion method is commonly used due to its efficiency with regard to computation costs. This application of 1D inversion is based on the assumption that small-offset broadband EM signals are insensitive to lateral resistivity variation. However, this assumption can be false when isolated conductive bodies such as man-made objects are embedded in the earth. Thus, we need to clarify the applicability of the 1D inversion method for small-loop EM data. In order to systematically analyze this conventional inversion approach, we developed a 2D EM inversion algorithm and verified this algorithm with a synthetic EM data set. 1D and 2D inversions were applied to synthetic and field EM data sets. The comparison of these inversion results shows that the resistivity distribution of the subsurface constructed by the 1D inversion approach can be distorted when the earth contains man-made objects, because they induce drastic variation of the resistivity distribution. By analyzing the integrated sensitivity of the small-loop EM method, we found that this pitfall of 1D inversion may be caused by the considerable sensitivity of the small-loop EM responses to lateral resistivity variation. However, the application of our 2D inversion algorithm to synthetic and field EM data sets demonstrate that the pitfall of 1D inversion due to man-made objects can be successfully alleviated. Thus, 2D EM inversion is strongly recommended for detecting conductive isolated bodies, such as man-made objects, whereas this approach may not always be essential for interpreting the EM field data.  相似文献   

10.
Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.  相似文献   

11.
Inversion of DC resistivity data using neural networks   总被引：9，自引：0，他引：9  
The inversion of geoelectrical resistivity data is a difficult task due to its non-linear nature. In this work, the neural network (NN) approach is studied to solve both 1D and 2D resistivity inverse problems. The efficiency of a widespread, supervised training network, the back-propagation technique and its applicability to the resistivity problem, is investigated. Several NN paradigms have been tried on a basis of trial-and-error for two types of data set. In the 1D problem, the batch back-propagation paradigm was efficient while another paradigm, called resilient propagation, was used in the 2D problem. The network was trained with synthetic examples and tested on another set of synthetic data as well as on the field data. The neural network gave a result highly correlated with that of conventional serial algorithms. It proved to be a fast, accurate and objective method for depth and resistivity estimation of both 1D and 2D DC resistivity data. The main advantage of using NN for resistivity inversion is that once the network has been trained it can perform the inversion of any vertical electrical sounding data set very rapidly.  相似文献   

12.
A combination of magnetotelluric (MT) measurements on the surface and in boreholes (without metal casing) can be expected to enhance resolution and reduce the ambiguity in models of electrical resistivity derived from MT surface measurements alone. In order to quantify potential improvement in inversion models and to aid design of electromagnetic (EM) borehole sensors, we considered two synthetic 2D models containing ore bodies down to 3000 m depth (the first with two dipping conductors in resistive crystalline host rock and the second with three mineralisation zones in a sedimentary succession exhibiting only moderate resistivity contrasts). We computed 2D inversion models from the forward responses based on combinations of surface impedance measurements and borehole measurements such as (1) skin-effect transfer functions relating horizontal magnetic fields at depth to those on the surface, (2) vertical magnetic transfer functions relating vertical magnetic fields at depth to horizontal magnetic fields on the surface and (3) vertical electric transfer functions relating vertical electric fields at depth to horizontal magnetic fields on the surface. Whereas skin-effect transfer functions are sensitive to the resistivity of the background medium and 2D anomalies, the vertical magnetic and electric field transfer functions have the disadvantage that they are comparatively insensitive to the resistivity of the layered background medium. This insensitivity introduces convergence problems in the inversion of data from structures with strong 2D resistivity contrasts. Hence, we adjusted the inversion approach to a three-step procedure, where (1) an initial inversion model is computed from surface impedance measurements, (2) this inversion model from surface impedances is used as the initial model for a joint inversion of surface impedances and skin-effect transfer functions and (3) the joint inversion model derived from the surface impedances and skin-effect transfer functions is used as the initial model for the inversion of the surface impedances, skin-effect transfer functions and vertical magnetic and electric transfer functions. For both synthetic examples, the inversion models resulting from surface and borehole measurements have higher similarity to the true models than models computed exclusively from surface measurements. However, the most prominent improvements were obtained for the first example, in which a deep small-sized ore body is more easily distinguished from a shallow main ore body penetrated by a borehole and the extent of the shadow zone (a conductive artefact) underneath the main conductor is strongly reduced. Formal model error and resolution analysis demonstrated that predominantly the skin-effect transfer functions improve model resolution at depth below the sensors and at distance of \(\sim \) 300–1000 m laterally off a borehole, whereas the vertical electric and magnetic transfer functions improve resolution along the borehole and in its immediate vicinity. Furthermore, we studied the signal levels at depth and provided specifications of borehole magnetic and electric field sensors to be developed in a future project. Our results suggest that three-component SQUID and fluxgate magnetometers should be developed to facilitate borehole MT measurements at signal frequencies above and below 1 Hz, respectively.  相似文献   

13.
Surface electrical and electromagnetic methods have a limited resolution capability for determining the conductivity structure of the earth. In one-dimensional modelling a collection of many thin layers is frequently considered as one composite layer, which is then macro-anisotropic. Neither galvanic methods nor inductive methods alone can resolve the anisotropy of the ground, but a joint inversion of galvanic and inductive data may do so. The necessity of including the coefficient of anisotropy in the joint inversion of galvanic and inductive sounding data is demonstrated. An analysis is made of the combined use of geoelectrical and transient soundings to resolve the coefficient of anisotropy of a subsurface layer for varying thickness, resistivity and coefficient of anisotropy. It is found that the coefficient of anisotropy is well resolved only for layers that are many times thicker than the overburden and for coefficients of anisotropy that are not too small. The ability of the joint inversion of geoelectrical and transient sounding data to resolve macro-anisotropic layers is tested using realistic earth models determined from electrical logs.  相似文献   

14.
Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.  相似文献   

15.
A code for 3-D resistivity modelling and inversion of vertical electrical soundings has been developed based on the finite-element technique and regularisation method. Synthetic data were used to test the effectiveness of the code and to examine the resolving capability of the Schlumberger array in investigating 3-D resistivity distributions. The code was applied to experimental data set constituted by 35 Schlumberger soundings collected near the Cairo city in order to study the subsurface resistivity distribution. The results have shown that valuable imaging of the subsurface resistivity distribution can be constructed even when the vertical electrical soundings are acquired in a sparse field data set.  相似文献   

16.
Global optimization with very fast simulated annealing (VFSA) in association with joint inversion is performed for 1D earth structures. The inherent problems of equivalence and suppression in electromagnetic (EM) and direct current (DC) resistivity methods are studied. Synthetic phase data from multifrequency sounding using a horizontal coplanar coil system and synthetic apparent resistivity data from Schlumberger DC resistivity measurements are inverted individually and jointly over different types of layered earth structures. Noisy data are also inverted. The study reveals that global optimization of individual data sets cannot solve inherent equivalence or suppression problems. Joint inversion of EM and DC measurements can overcome the problem of equivalence very well. However, a suppression problem cannot be solved even after combination of data sets. This study reveals that the K-type earth structure is easiest to resolve while the A-type is the most difficult. We also conclude that the equivalence associated with a thin resistive layer can be resolved better than that for a thin conducting layer.  相似文献   

17.
Helicopter-borne frequency-domain electromagnetic (HEM) surveys are used for fast high-resolution, three-dimensional resistivity mapping. Standard interpretation tools are often based on layered earth inversion procedures which, in general, explain the HEM data sufficiently. As a HEM system is moved while measuring, noise on the data is a common problem. Generally, noisy data will be smoothed prior to inversion using appropriate low-pass filters and consequently information may be lost.For the first time the laterally constrained inversion (LCI) technique has been applied to HEM data combined with the automatic generation of dynamic starting models. The latter is important because it takes the penetration depth of the electromagnetic fields, which can heavily vary in survey areas with different geological settings, into account. The LCI technique, which has been applied to diverse airborne and ground geophysical data sets, has proven to be able to improve the HEM inversion results of layered earth structures. Although single-site 1-D inversion is generally faster and — in case of strong lateral resistivity variations — more flexible, LCI produces resistivity — depth sections which are nearly identical to those derived from noise-free data.The LCI results are compared with standard single-site Marquardt–Levenberg inversion procedures on the basis of synthetic data as well as field data. The model chosen for the generation of synthetic data represents a layered earth structure having an inhomogeneous top layer in order to study the influence of shallow resistivity variations on the resolution of deep horizontal conductors in one-dimensional inversion results. The field data example comprises a wide resistivity range in a sedimentary as well as hard-rock environment.If a sufficient resistivity contrast between air and subsurface exists, the LCI technique is also very useful in correcting for incorrect system altitude measurements by using the altitude as a constrained inversion parameter.  相似文献   

18.
Geophysical methods were applied for hydrogeological targets in many countries including Vietnam. This paper presents results of using complex geophysical techniques as well as 2D electrical resistivity imaging (ERI), vertical electrical sounding (VES), very low frequency (VLF), and seismic refraction for geological structure investigation for locating the aquifers and assessing the hydrogeological conditions for groundwater potential in industrial zones of North Hanoi, Vietnam. The locations of two aquifers are determined by their depth and thickness on the basis of resistivity and seismic velocity values which were proved by stratifications of three boreholes to 40–60 m of depth on the study area. There are connections from surface water to shallow aquifer by hydraulic windows, as follows from VLF data. The deeper aquifer can be considered as a potential groundwater supply, but the water level is descending in time, as shown by hydrological monitoring. However, with careful use and by reducing sources of pollution, groundwater can continue to be an important natural resource for future.  相似文献   

19.
Airborne electromagnetic (AEM) surveys are currently being flown over populated areas and applied to detailed problems using high flight line densities. Interpretation information is supplied through a model of the subsurface resistivity distribution. Theoretical and survey data are used here to study the character and reliability of such models. Although the survey data were obtained using a fixed-wing system, the corresponding associations with helicopter, towed-bird systems are discussed. Both Fraser half-space and 1D inversion techniques are considered in relation to their ability to distinguish geological, cultural and environmental influences on the survey data. Fraser half-space modelling provides the dual interpretation parameters of apparent resistivity and apparent depth at each operational frequency. The apparent resistivity was found to be a remarkably stable parameter and appears robust to the presence of a variety of at-surface cultural features. Such features provide both incorrect altitude data and multidimensional influences. Their influences are observed most strongly in the joint estimate of apparent depth and this accounts for the stability of the apparent resistivity. Positive apparent depths, in the example data, result from underestimated altitude measurements. It is demonstrated that increasingly negative apparent depths are associated with increasing misfits between a 1D model and the data. Centroid depth calculations, which are a transform of the Fraser half-space parameters, provide an example of the detection of non-1D influences on data obtained above a populated area. 1D inversion of both theoretical and survey data is examined. The simplest use of the 1D inversion method is in providing an estimate of a half-space resistivity. This can be undertaken prior to multilayer inversion as an initial assessment. Underestimated altitude measurements also enter the problem and, in keeping with the Fraser pseudo-layer concept, an at-surface highly resistive layer of variable thickness can be usefully introduced as a constrained parameter. It is clearly difficult to ascribe levels of significance to a ‘measure’ of misfit contained in a negative apparent depth with the dimensions of metres. The reliability of 1D models is better assessed using a formal misfit parameter. With the misfit parameter in place, the example data suggest that the 1D inversion methods provide reliable apparent resistivity values with a higher resolution than the equivalent information from the Fraser half-space estimates.  相似文献   

20.
It is proposed that the Straightforward Inversion Scheme (SIS) developed by the authors for 1D inversion of resistivity sounding and magneto-telluric sounding data can also be used in similar fashion for time-domain induced polarization sounding data. The necessary formulations based on dynamic dipole theory are presented. It is shown that by using induced polarization potential, measured at the instant when steady state current is switched off, an equation can be developed for apparent ‘chargeability–resistivity’ which is similar to the one for apparent resistivity. The two data sets of apparent resistivity and apparent chargeability–resistivity can be inverted in a combined manner, using SIS for a common uniform thickness layer earth model to estimate the respective subsurface distributions of resistivity and chargeability–resistivity. The quotient of the two profiles will give the sought after chargeability profile. A brief outline of SIS is provided for completeness. Three theoretical models are included to confirm the efficacy of SIS software by inverting only the synthetic resistivity sounding data. Then one synthetic data set based on a geological model and three field data sets (combination of resistivity and IP soundings) from diverse geological and geographical regions are included as validation of the proposal. It is hoped that the proposed scheme would complement the resistivity interpretation with special reference to shaly sand formations.  相似文献   

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

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