共查询到20条相似文献，搜索用时 31 毫秒
1.
Piecewise 1D laterally constrained inversion of resistivity data 总被引：1，自引：0，他引：1
Esben Auken ers V. Christiansen Bo H. Jacobsen Nikolaj Foged Kurt I. Sørensen 《Geophysical Prospecting》2005,53(4):497-506
In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non‐geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning. We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well‐resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result. Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example. 相似文献
2.
Yves Albouy Pierre Andrieux Gérard Rakotondrasoa Michel Ritz Marc Descloitres Jean-Lambert Join Eddy Rasolomanana 《Ground water》2001,39(1):87-97
Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently. 相似文献
3.
水平井钻进过程中,地层边界的准确预测对地质导向决策至关重要.但水平井测井环境复杂,随钻方位电磁波测井响应不直观,储层性质、地层界面等信息的精确、定量评价困难.本文提出了一种适用于大斜度井/水平井随钻方位电磁波测井资料的电阻率及地层界面信息的实时提取方法.该方法主要基于以下几个技术:(1)引入滑动开窗策略将高维随钻方位电磁波测井资料的反演问题转化为多个连续窗口的一维反演问题;(2)对反演参数进行多初值初始化,结合正则化Levenberg-Marquardt算法保证解的全局最优性;(3)充分利用随钻方位电磁波测井多种探测深度信息,采用分级反演策略提高反演速度.数值和实测资料反演结果表明:降维反演方法可实现对地层界面的实时、准确提取,且适用于多数复杂地层结构和任意井眼轨迹;利用浅探测地质信号可实现随钻方位电磁波测井资料的分级反演,最终反演速度可提升10倍以上. 相似文献
4.
An iterative refinement method for determining a layered resistivity model from a Schlumberger or Wenner sounding curve is adapted to determine a layered resistivity model by using apparent resistivity and phase derived from the magnetotelluric impedance. Magnetotelluric observations presented as a function of period are first converted to an approximate resistivity–depth profile using Schmucker's transformation and this is used to construct an initial guess (starting) model. A two-stage procedure is then invoked. Keeping resistivities constant, layer boundaries are first adjusted to give a minimum misfit between measured data and responses and this is followed by resistivity adjustments with fixed layer boundaries to reduce the misfit further. The method is illustrated by application to some synthetic data both exact and with added noise, to a real field data set and to some magnetotelluric profile data obtained in a survey over the Carnmenellis granites in south Cornwall. The method is validated by recovering conductivity models from the exact and noisy 1D synthetic data. For complicated three-dimensional data at a single site and along a profile of stations, the method is shown to produce acceptable solutions which may be used as starting models in further two- or three-dimensional studies. 相似文献
5.
A finite number of layers must normally be used to represent the electrical properties of a horizontally-layered region. It is well known that these representative layers must include anisotropy. Applying the concept of joint inversion to DC resistivity and magnetotelluric (electromagnetic) data from a structure consisting of many fine layers, we show how the resistivities in the horizontal and vertical directions can be separately determined in each layer. The results are related to the contrasts existing in the interbeds which comprise the composite layers. The method has applications in both engineering and exploration. 相似文献
6.
Introduction Jiashi-Artux area in southwest Xinjiang is one of the most active earthquake provinces at pre-sent in Chinese mainland. In the last century, about 3/4 strong earthquakes in Chinese mainland hit this area, and especially from January 21 to April 16 in 1997, 7 earthquakes with the magnituderanging from 6.0 to 6.9 occurred in a very small area of 9 km×18 km near Jiashi (ZHU et al, 1998). It has never taken place before in Chinese mainland that a series of strong earthquakes shoo… 相似文献
7.
Seismic impedance inversion is a well-known method used to obtain the image of subsurface geological structures. Utilizing the spatial coherence among seismic traces, the laterally constrained multitrace impedance inversion (LCI) is superior to trace-by-trace inversion and can produce a more realistic image of the subsurface structures. However, when the traces are numerous, it will take great computational cost and a lot of memory to solve the large-scale matrix in the multitrace inversion, which restricts the efficiency and applicability of the existing multitrace inversion algorithm. In addition, the multitrace inversion methods are not only needed to consider the lateral correlation but also should take the constraints in temporal dimension into account. As usual, these vertical constraints represent the stratigraphic characteristics of the reservoir. For instance, total-variation regularization is adopted to obtain the blocky structure. However, it still limits the magnitude of model parameter variation and therefore somewhat distorts the real image. In this paper, we propose two schemes to solve these issues. Firstly, we introduce a fast algorithm called blocky coordinate descent (BCD) to derive a new framework of laterally constrained multitrace impedance inversion. This new BCD-based inversion approach is fast and spends fewer memories. Next, we introduce a minimum gradient support regularization into the BCD-based laterally constrained inversion. This new approach can adapt to sharp layer boundaries and keep the spatial coherence. The feasibility of the proposed method is illustrated by numerical tests for both synthetic data and field seismic data. 相似文献
8.
The problem of equivalence in direct current (DC) resistivity and electromagnetic methods for a thin resistive and conducting layer is well‐known. Attempts have been made in the past to resolve this problem through joint inversion. However, equivalence still remains an unresolved problem. In the present study, an effort is made to reduce non‐uniqueness due to equivalence using global optimization and joint inversion by successive refinement of the model space. A number of solutions derived for DC resistivity data using very fast simulated annealing global inversion that fits the observations equally well, follow the equivalence principle and show a definite trend. For a thin conductive layer, the quotient between resistivity and thickness is constant, while for a resistive one, the product between these magnitudes is constant. Three approaches to obtain very fast simulated annealing solutions are tested. In the first one, layer resistivities and thicknesses are optimized in a linear domain. In the second, layer resistivities are optimized in the logarithmic domain and thicknesses in the linear domain. Lastly, both layer resistivities and thicknesses are optimized in the logarithmic domain. Only model data from the mean models, corresponding to very fast simulated annealing solutions obtained for approach three, always fit the observations. The mean model defined by multiple very fast simulated annealing solutions shows extremely large uncertainty (almost 100%) in the final solution after inversion of individual DC resistivity or electromagnetic (EM) data sets. Uncertainty associated with the intermediate resistive and conducting layers after global optimization and joint inversion is still large. In order to reduce the large uncertainty associated with the intermediate layer, global optimization is performed over several iterations by reducing and redefining the search limits of model parameters according to the uncertainty in the solution. The new minimum and maximum limits are obtained from the uncertainty in the previous iteration. Though the misfit error reduces in the solution after successive refinement of the model space in individual inversion, it is observed that the mean model drifts away from the actual model. However, successive refinement of the model space using global optimization and joint inversion reduces uncertainty to a very low level in 4–5 iterations. This approach works very well in resolving the problem of equivalence for resistive as well as for conducting layers. The efficacy of the approach has been demonstrated using DC resistivity and EM data, however, it can be applied to any geophysical data to solve the inherent ambiguities in the interpretations. 相似文献
9.
Mohammad Israil 《Acta Geophysica》2006,54(2):126-141
Two techniques have been presented for the delineation of boundaries from smooth models obtained by smooth inversion techniques
of geoelectrical sounding data, such as straightforward inversion scheme, Occam’s and Zohdy’s methods. The smooth model consists
of a large number of equally spaced layers, wherein the real geological boundaries are missing. The techniques proposed here
suppress the geologically irrelevant boundaries and support the real structural boundaries present in the geoelectrical data.
In the first technique, solution of linear inverse problem is improved iteratively through weighted minimum norm inverse,
the weight being taken from the current solution. The technique is referred as Iterative Straightforward Inversion Scheme.
The second method is analytical, based on the application of smoothing filter, referred in the literature as edge-preserving
smoothing. A few examples of theoretical magnetotelluric, dc resistivity and field sounding data have been presented to demonstrate
the capabilities of the techniques. The methodologies also reduce the conspicuous oscillations in the smooth solutions caused
by the conversion of sharp boundaries to the smooth ones. 相似文献
10.
Gylfi Pll Hersir Axel Bjrnsson Laust Brsting Pedersen 《Journal of Volcanology and Geothermal Research》1984,20(3-4)
The results of 10 magnetotelluric soundings, performed along a 110-km-long profile crossing the constructive plate boundary in southwest Iceland, are presented. Apparent resistivities are interpreted by a horizontally stratified earth model to yield a pseudo cross-section along the profile. The crust-mantle interface contains a well conductive layer. The depth to the good conductor increases with age of the crust and the distance from the axial zone. This layer is interpreted as partially molten basalt, at a temperature about 1100°C and a volume fraction of the melt phase in the range 10–20%.The high-conductivity layer probably disappears west of the Borgarnes anticlinal axis, which separates the older (to the west) and younger (to the east) flood basalts in western Iceland, indicating that the temperature below the oldest part of the profile lies below the solidus curve of basalt. Recent seismic crustal investigations in the same area indicate a state of partial melting or a magma chamber, which agrees with the results of the magnetotelluric soundings. 相似文献
11.
L.P. Vinnik R.A. Avetisjan N.G. Mikhailova 《Physics of the Earth and Planetary Interiors》1983,33(3):149-163
The method of detection of P-to-SV converted waves from distant earthquakes (Vinnik, 1977) was applied to sets of long-period records from a few seismograph stations in Europe and the west of North America. The results obtained suggest that the converted phases related to the major boundaries in the mantle can be reliably detected and the depths of conversion evaluated with an accuracy of a few kilometres. The depth of the olivine-spinel transition is close to 400 km and no difference between the estimates for the north of Europe and the west of North America is found. The depth of the boundary separating the upper and lower mantle is close to 640 km, which is 30 km less than in the recent Earth-reference models. Fine S velocity stratification of this transition changes laterally from a high-gradient layer 50 km thick, terminated at the bottom by a sharp discontinuity, to a gradient layer 100 km or more thick without the discontinuity. A striking anomaly of the mantle transition zone is found in the Rio Grande rift area where a well pronounced boundary is found at 510 km depth. 相似文献
12.
A 3-D crustal geoelectric model of the Ukrainian Shield (USh) is constructed from magnetic variation data. High electrical conductivity anomalies with resistivities of 1–100 Ω m are located at depths of up to 30 km from the basement surface. A high conductivity layer with ρ = 25 Ω m and with its upper boundary at a depth of 70 km is supposed to exist in the upper mantle of the southwestern USh. 相似文献
13.
Niels B. Christensen 《Geophysical Prospecting》2000,48(1):1-19
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.
Three separate conditions controlling the structure of the inversion are discussed. The rate of change of boundary layer temperature and the vertical gradient of potential temperature in the overlying air set the maximum limit to possible rise rates. In transient conditions, the rise of the inversion layer is sometimes dominated by turbulent transport and erosion, however this cannot continue indefinitely. The horizontal advection above the inversion can also change the relative temperature increase at the inversion, and so may partly determine inversion history.In this paper the theory of the plume-connective boundary layer has been extended to allow the plumes to erode material from the overlying stable air, and so permit study of how the boundary layer could deepen by eroding upward into the capping inversion. Only clear-air processes without condensation are included in the model described here.When the mechanism of convection is examined, it is demonstrated that, because of energetic considerations, the rate of erosion, by buoyancy generated turbulence of the overlying air down across the inversion, will not usually determine the rate of rise of the inversion, except briefly. Other effects often dominate in determining the average erosion rate. Furthermore, because of the limitations set by the buoyancy of the entrained air, sharp inversions, exceeding about 1°C, cannot be formed without the action of latent heat from evaporating droplets. An intermediate buffer layer forms in the inversion zone whenever the erosion rate is high enough to produce a temperature increase above the boundary layer, sufficient to inhibit direct entrainment. 相似文献
15.
16.
Many joint inversion schemes use 1D forward modelling in the integrated interpretation of various geophysical data. In extending the joint inversion approach to the investigation of 2D structures, the discretization of the model parameters and the appropriate choice of the forward‐modelling procedure play a very important role. In this paper, a hybrid seismic–geoelectric joint inversion method is proposed for the investigation of 2D near‐surface geological structures. The electric and seismic models are coupled together through the use of common boundaries between the adjacent layers. Assuming a 2D model composed of homogeneous layers with curved boundaries, a fast ray‐tracing algorithm is used for the calculation of refraction seismic traveltime data. In the geoelectric forward modelling, a locally 1D approximation is used. The boundary surfaces are written in the form of series expansion; the inversion algorithms are formulated for the expansion coefficients and the petrophysical parameters as unknowns. Two versions of the inversion method are proposed: in versions A and B, interval‐wise constant functions and Chebyshev polynomials are, respectively, used as basis functions of the series expansion. The versions are tested by means of synthetic and in situ measured data. The tests show that both methods are stable and accurate. 相似文献
17.
K. MALLICK 《Geophysical Prospecting》1970,18(Z1):738-757
This work deals with a theoretical investigation of the magnetotelluric sounding where the horizontal layers of a geological section do not have a sharp but a transitional boundary. The ratio of the horizontal electric to the horizontal magnetic components gives a measure of “Cagniard impedance”, which together with the phase relationship between them, enables one to determine the subsurface geometry and the true resistivities. Utilising the principle of similitude a set of Resistivity Master Curves and the corresponding phase relationships are presented here for the case where the lower half-space is infinitely resistive. The relation between the relative transition thickness and the cut-off period is illustrated. 相似文献
18.
The electrical potential due to a point source of current placed on the ground surface is studied for a multi-layered earth consisting of homogeneous overburden of constant conductivity over a stack of transition layers where conductivity varies with depth according to power and exponential laws in even and odd layers, respectively. The general recursion relations are derived and are used to obtain expressions for the apparent resistivities for Schlumberger and Wenner electrode arrays. Their asymptotic behavior has been studied. The solutions for some particular cases are given: (i) odd layers with conductivities exponentially varying with depth while all even layers (and the first) have constant conductivities; (ii) even layer conductivities varying as a power law while odd number layers are of constant conductivity; and (iii) any two successive layers as transition layers and all others having constant conductivities. Further it is shown that Patella's theory is a particular case of the present study. It is concluded, therefore, that the present treatment is more general as all earlier models consisting of trasition layers can be derived from this study. 相似文献
19.
A radiomagnetotelluric survey on an oil-contaminated area near the Brazi Refinery, Romania 总被引：2，自引：0，他引：2
Scalar radiomagnetotelluric measurements were carried out on a contaminated test area close to the Brazi Refinery in Romania in order to detect and to monitor a 1 m thick oil layer expected at 5 m depth. Radio transmitters broadcasting in a frequency range from 10 kHz to 300 kHz were selected to observe the apparent resistivity and the phase data associated with the E‐ and B‐polarizations. They were located parallel and perpendicular to the assumed strike direction of the contamination plume. The data were interpreted by a 2D inversion technique from which the conductivity structure of the area was derived. The 2D inversion models of all profiles on the contaminated area show a poor‐conductivity zone above the groundwater table which could be associated with the oil contamination. A first attempt was also made to monitor the contaminated layer: the radiomagnetotelluric measurements were repeated on the same profiles a year later, but this time in a dry period, not in a rainy one. The 2D inversion results of the measurements in the dry period indicate that the high‐resistivity layer moved closer to the surface. Additional reference measurements were then carried out on a non‐contaminated area situated at a distance from the refinery, in the opposite direction to the flow of the groundwater. These reference measurements were used for the derivation of the unperturbed geology and they were also compared with the measurements of the contaminated test area. There is a significant difference in the frequency dependences of the apparent resistivities of the reference and contaminated areas, which could indicate a contamination at shallow depth. The 2D inversion results show the increase of resistivity at a depth of about 5 m beneath the contaminated area where the oil contamination is expected according to the information from the boreholes. 相似文献
20.
Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field. 相似文献