首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We investigate a novel way to introduce resistivity models deriving from airborne electromagnetic surveys into regional geological modelling. Standard geometrical geological modelling can be strengthened using geophysical data. Here, we propose to extract information contained in a resistivity model in the form of local slopes that constrain the modelling of geological interfaces. The proposed method is illustrated on an airborne electromagnetic survey conducted in the region of Courtenay in France. First, a resistivity contrast corresponding to the clay/chalk interface was interpreted confronting the electromagnetic soundings to boreholes. Slopes were then sampled on this geophysical model and jointly interpolated with the clay/chalk interface documented in boreholes using an implicit 3D potential‐field method. In order to evaluate this new joint geophysical–geological model, its accuracy was compared with that of both pure geological and pure geophysical models for various borehole configurations. The proposed joint modelling yields the most accurate clay/chalk interface whatever the number and location of boreholes taken into account for modelling and validation. Compared with standard geological modelling, the approach introduces in between boreholes geometrical information derived from geophysical results. Compared with conventional resistivity interpretation of the geophysical model, it reduces drift effects and honours the boreholes. The method therefore improves what is commonly obtained with geological or geophysical data separately, making it very attractive for robust 3D geological modelling of the subsurface.  相似文献   

2.
The possibility of a time‐domain electromagnetic sounding method using excitation and measurement of vertical electric fields to search for and identify deeply buried reservoirs of hydrocarbons offshore is investigated. The method operates on source–receiver offsets, which are several times less than the depth of the reservoir. Geoelectric information is obtained from the transient responses recorded in the pauses between the pulses of electric current in the absence of the source field. The basics of the method, as well as its sensitivity, resolution, and the highest accessible depth of soundings for various geological conditions in a wide range of sea depths, are analyzed. For the analysis, 1D and 3D geoelectric models of hydrocarbon reservoirs are used. It is shown that under existing technologies of excitation and measurement of vertical electric fields, the highest accessible depth of soundings can be up to 4 km. Technology for the inversion and interpretation of transient responses is demonstrated on experimental data.  相似文献   

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

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

5.
The results of ten tensor magnetotelluric soundings recorded within the Charlevoix crater are presented. The geological complexity of the region (the crater is located at the junction of two major geological provinces, i.e. the crystalline basement of the Canadian shield and the more conductive sediments of the Appalachian system) and the amount of cultural noise present in our data prevent a highly definitive interpretation of the shallow structure. A subhorizontal conductive fault related to Rondot's (1971) structural sketch of the original impact structure may exist within the upper 5 km, but it is not completely constrained by the data. Pseudo 3D modelling using stacked thin sheets indicates that superficial patches of high conductivity may explain the large anisotropy observed at the long periods and suggests that our deeper structures fit well with the 2D model previously proposed by Kurtz (1982) for the region.  相似文献   

6.
One of the major aspects of rock-physics forward modelling is to predict seismic behaviour at an undrilled location using drilled well data. It is important to model the rock and fluid properties away from drilled wells to characterize the reservoir and investigate the root causes of different seismic responses. Using the forward modelling technique, it is possible to explain the amplitude responses of present seismic data in terms of probable rock and reservoir properties. In this context, rock-physics modelling adds significant values in the prospect maturation process by reducing the risk of reservoir presence in exploration and appraisal phases. The synthetic amplitude variation with offset gathers from the forward model is compared with real seismic gathers to ensure the fidelity of the existing geological model. ‘Prospect A’ in the study area has been identified from seismic interpretation, which was deposited as slope fan sediments in Mahanadi basin, East Coast of India. The mapped prospect has shown class-I amplitude variation with offset response in seismic without any direct hydrocarbon indicator support. The existing geological model suggests the presence of an excellent gas reservoir with proven charge access from the fetch area, moderate porosity and type of lithology within this fan prospect. But, whether the seismic response from this geological model will exhibit a class-I amplitude variation with offset behaviour or ‘dim spot’ will be visible; the objective of the present study is to investigate these queries. A rock-physics depth trend analysis has been done to envisage the possibilities of class-I reservoir in ‘Prospect A’. Forward modelling, using a combination of mechanical and chemical compaction, shows the synthetic gas gathers at ‘Prospect A’, which are class I in nature. The study has also depicted 2D forward modelling using lithology and fluid properties of discovery well within similar stratigraphy to predict whether ‘dim spot’ will be seen in seismic. The estimated change in synthetic amplitude response has been observed as ∼5% at contact, which suggests that the changes will not be visible in seismic. The study connects the existing geological model with a top-down seismic interpretation using rock-physics forward modelling technique to mature a deep-water exploratory prospect.  相似文献   

7.
When electric soundings are made over an irregular terrain, topographic effects can influence the values of apparent resistivity and lead to erroneous 1D interpretation. A 3D finite-element method has been applied to study the topographical effect of a slope on Schlumberger soundings parallel to the strike. When the resistivity survey is performed at the top of the slope, the apparent resistivity values can be two times higher than in the flat-earth case, depending on the angle (α) and height (H) of the slope, and on the distance (X) between the sounding and the slope top. The results are presented as nondimensional curves which can be used for evaluating topographic anomalies for any value of the parameters α, H and X. It is numerically shown that the topographic effects can be removed from measurements on horizontally layered structures with an irregular earth surface. Real measurements were performed in different geological conditions over an irregular terrain. The correction method based on the nondimensional curves has been applied to the data and has enabled the determination of the correct layered ground configuration using 1D interpretation.  相似文献   

8.
South America presents several unique geomagnetic features, one of which is the Equatorial Electrojet (EEJ), a current system which extends itself east to west in Northeastern Brazil for almost 3500 km. Considering the fact that the influence of this phenomenon at low magnetic latitudes can be detected at great distances from its central axis, it is important to simulate its effect in magnetotelluric (MT) exploration. To accomplish this and by the use of an integral equation scheme, we have modeled the MT response of confined three-dimensional (3D) bodies (dykes in an homogeneous half-space) and deep 3D structures (horst and graben located at Marajo Basin in Northern Brazil). As the inductive source, we replace the classical plane wave source by a common line of current, besides gaussian and undulatory electrojets. The results of the modelling indicate that the studied effect is prominent in both one-dimensional (1D) and 3D media for periods (T) greater than 10 s. It decreases with distance, but it is detectable as far as 3000 km from the center of the EEJ. It is also observed that for T greater than 10 s, the computed soundings can be strongly distorted, mainly by EEJ effects in the host medium which, in turn, cause changes in the final branch of the 3D soundings. For structures in the Marajo Basin, our results reveal that the 3D galvanic effect prevails in the interval 0.1<T<10 s, while source influence occurs mainly for T>10 s. On the other hand, automatic inversion of these 3D data shows that, in equatorial regions having complex geology, some errors arise in conventional 1D interpretation of the MT soundings. This is due to the superposition of the host medium response, the galvanic effect of 3D structures, and the EEJ influence.  相似文献   

9.
当地表存在三维非均匀电导率分布时,区域大地电磁响应发生畸变. 以往对这种畸变研究多假设近地表为三维,区域构造为一维或二维. 对于更一般的三维/三维构造,为了分析并消除这种畸变影响,真实反映地下三维区域构造信息,本文实现了三维大地电磁相位张量积分方程数值算法,并研究在不同地质模型下相位张量响应. 结果表明,相位张量不仅可以反映一般三维构造信息,亦可有效反映复杂近地表构造下三维区域构造信息,而无须假设区域构造为一维或二维,证明相位张量具有较强抗近地表局部非均匀构造干扰能力,能够保持更为一般的三维区域构造信息. 为了加快正演计算,同时保持一定精度,算法采用了积分方程多网格法.  相似文献   

10.
Magnetotelluric (MT) soundings and gravity methods were employed to study the deep freshwater aquifer in the area north of Abo Zenema city on the eastern side of the Gulf of Suez, Egypt. Seven MT sites and 48 gravity stations were surveyed along northeast–southwest profiles as close as possible to a line perpendicular to the coast of the Gulf of Suez. The MT survey was conducted using high and low frequencies to investigate shallow and deep areas, respectively. One-dimensional inversion was conducted using a heuristic inversion scheme of the Bostick algorithm. The MT data were also inverted with a 2-D smooth model inversion routine using the nonlinear conjugate gradient method to infer variation in vertical and lateral resistivity inside the Earth. A 100-Ohm-m homogeneous half-space initial model was used to invert the TE mode data only. Then, the inverted model obtained from the TE mode data was used as an initial model for inversion of the TM mode data. The inverted model thus obtained from the TM mode data inversion was used as an initial model for the inversion of the joint TE and TM responses. Two-dimensional (2-D) forward modeling of the gravity data was conducted using the 2-D polygon method of Talwani’s algorithm for an arbitrarily shaped body and was based on the subsurface information from the MT survey and the available information about the geological structure of the study area. This method enabled us to obtain the basement structure of the coastal aquifer in the study area. The results from the analysis and the interpretation of MT and gravity data were used to detect and delineate the groundwater coastal aquifer in the study area.  相似文献   

11.
The efficacy of the magnetotelluric and audiomagnetotelluric (MT/AMT) methods for detailing the structure of a hypothetical geological section is investigated by using the singular value decomposition (SVD) technique. The section is representative of southeastern Turkey, which is mostly covered by basalt and is a prime area for oil exploration. One of the geological units, the Germav shale at a depth of 600 m, is a problem layer for electromagnetic surveys because of its very low resistivity (on average 3 Ωm) and highly variable thickness across the area (200–900 m). In the MT frequency range (0.0004–40 Hz) its total conductance—or, since its resistivity is known from resistivity log information, its thickness—is the best resolved model parameter. The total depth to the Germav shale and the resistivity of the Cambrian/Precambrian basement are the marginally resolved parameters. In the AMT frequency range (4–10000 Hz) the resistivity of the surface basalt layer strongly affects the resolution of the other, less important, model parameters which are the total depth to the Germav shale and the total conductance of the Germav shale. The errors in the measurements determine the number of model parameters resolvable, and are also important for interpretation of the geological model parameters to within a desired accuracy. It is shown that statistical evaluation of the MT and/or AMT interpretations by using an SVD factorization of the sensitivity matrix can be helpful to define the importance of some particular stage of the interpretation, and also provides a priori knowledge to plan a proposed survey. Arrangements of MT and AMT observations, together with some Schlumberger resistivity soundings, on a large grid will certainly provide three-dimensional detailed information of the deep geoelectric structure of the area.  相似文献   

12.
Resistivity measurements were carried out in a survey area in the south of Germany. This area is characterized by complicated subsurface geology. Schlumberger full-arrays and their respective half-arrays were recorded simultaneously. The results obtained by the one-dimensional (1D) interpretation of the full-array measurements were incorrect because of a resistivity discontinuity. This discontinuity, under a relatively thick overburden, could only be located by the half-array soundings. Its exact location and the resistivity distribution in the subsurface were ascertained by comparing the sounding curves with 2D model curves, which are calculated by a finite-difference method.  相似文献   

13.
The theory of electrical dipole soundings proved that this method can produce resistivity measurements, which are comparable with those obtained by electrical soundings of the Wenner or Schlumberger type. Their main advantage is the use of short cable lengths, which is important if the depth of penetration should be large. A considerable disadvantage of the dipole method is the great sensitivity to lateral discontinuities. Though these have an influence on the Schlumberger arrangement as well, they can disturb a dipole sounding to such an extent than an interpretation based on a horizontal layer case is no more possible. There are six different dipole arrays, which differ from each other with respect to the angle enclosed by the two dipole orientations-the current dipole AB and the measuring dipole MN. The theoretical comparison of the dipole arrays with the Schlumberger array concerning their sensitivity to lateral discontinuities is a useful basis for the choice of the most suitable configuration. Considering geological subsurface conditions the right choice of a dipole array can give an optimal result, i.e. a dipole sounding for which the sensitivity to lateral discontinuities is as small as possible under the given circumstances.  相似文献   

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

15.
We invert 2D surface gravity data constrained both by geological and seismic information. We use a number of pre-processing tools in order to reduce the general multi-body inversion into several single-body inversions, whereby we can reduce the overall complexity of the inversion task. This is done with as few assumptions as possible. Furthermore, for a single-body inversion we uncouple the determination of the shape of the causative sources from the determination of their mass density contrast to the surroundings. The inversion for the geometrical shape of the source body is done in steps. Firstly, a rough 3D shape of the source is modelled—a model consisting of the vertical mass columns of equal height. The horizontal extension is implied by the surface gravity signal. Subsequently, the shape of each source body is modified to obtain a better fit to the surface gravity data. In each modification step, the overall change of the shape of the source body is followed by an update of the mass density contrast to the surroundings. The technique was applied to a set of gravity data from the Eastern Goldfield area in Western Australia. The area has been widely studied in the past. In 1999, two seismic profiles that cross-sect the area were measured. Furthermore, an extensive geological modelling for the area has been conducted. The practical goal of this work was to verify the geological interpretation using the potential field data (mainly the gravity data although magnetic data were also available) and only weakly constrained by the seismic information. The result was the reconstruction of the ‘rough’ 3D geometry of the source bodies and the estimation of a constant mass density contrast to the surroundings. A possible extension of this technique for detailed studies of the geological model is briefly discussed.  相似文献   

16.
The possibilities of reconstructing two-dimensional (2D) cross sections based on the data of the profile soundings by the transient electromagnetic method (TEM) with a single ungrounded loop are illustrated on three-dimensional (3D) models. The process of reconstruction includes three main steps: transformation of the responses in the depth dependence of resistivity ρ(h) measured along the profile, with their subsequent stitching into the 2D pseudo section; point-by-point one-dimensional (1D) inversion of the responses with the starting model constructed based on the transformations; and correction of the 2D cross section with the use of 2.5-dimensional (2.5D) block inversion. It is shown that single-loop TEM soundings allow studying the geological media within a local domain the lateral dimensions of which are commensurate with the depth of the investigation. The structure of the medium beyond this domain insignificantly affects the sounding results. This locality enables the TEM to reconstruct the geoelectrical structure of the medium from the 2D cross sections with the minimal distortions caused by the lack of information beyond the profile of the transient response measurements.  相似文献   

17.
Adverse geologies are often encountered during tunnel construction, which could seriously endanger the construction. To ensure the safety, it is essential to detect adverse geologies and their water‐bearing situation ahead the tunnel face. Ground‐penetrating radar is a suitable instrument, but the accurate interpretation of its detection results is difficult. In this paper, at first, an improved back projection imaging algorithm is proposed, which can make reflection waves closer to the real geological boundaries with few artificial clutters. And then, forward modelling of ground‐penetrating radar is carried out for typical adverse geologies, such as karst caves, faults, fractured rock masses, fracture network, and water‐bearing body. Their corresponding response features are obtained, accumulating experience for geological interpretation. The above two methods provide the basis for target identification and geological interpretation. In the last part, the application of the above two methods in several engineering cases are given, and their effectiveness is verified.  相似文献   

18.
A magnetotelluric acquisition and modelling campaign, comprising 400 full tensor soundings, was carried out in the Sub-Andean Foothills in the area located south of Santa Cruz de la Sierra, Bolivia. The objective of the survey was to improve the imaging of the Paleozoic section and provide insight to the overall structure of the Devonian Silurian complex. The acquired data were inverted assuming both two-dimensional and three-dimensional dimensionality of the subsurface structure, adopting a multi-step iterative workflow, which began applying unconstrained inversion and introduced increasing constraints determined at each following step based on the interpretation of the latest inversion results. The geological interpretation of the final inversion results allowed to better image the expected complexity in the structures, particularly highlighting several different levels of detachment, mainly in the area of Kirusillas Formation. Furthermore, an oblique sense was identified of the main tendency of N-S deformation, related to a dextral transgressional movement that generates discrete dislocations, at a length scale of the order of a few kilometres.  相似文献   

19.
The Belvedere Spinello salt mine is located in the Catanzaro Province of Calabria in Southern Italy. An extensive mining program has caused the development of Underground cavities filled with brine and the migration of this brine has been of great environmental concern to the mine owners. This paper presents the results of a multidimensional interpretation of a two-phase resistivity and magnetotelluric (MT) survey that was performed in an attempt to determine the complex conductivity structure of the mine area and to gain information on brine development and migration pathways. Key resistivity soundings were interpreted using a 2.5D algorithm based on the Polozhii decomposition method. The MT data were interpreted using a 2D finite-element code. A conductivity model was developed, integrating available geological and drill-hole information. The interpretation of the MT data, collected five years after the acquisition of the resistivity data, shows a conductive feature of depth that is not resolved in the resistivity interpretation. This feature has been interpreted as a thick brine zone that has developed as a result of mining during the interval between the resistivity and the MT measurements.  相似文献   

20.
The oil shale exploration program in Jordan is undertaking great activity in the domain of applied geophysical methods to evaluate bitumen‐bearing rock. In the study area, the bituminous marl or oil shale exhibits a rock type dominated by lithofacies layers composed of chalky limestone, marls, clayey marls, and phosphatic marls. The study aims to present enhancements for oil shale seam detection using progressive interpretation from a one‐dimensional inversion to a three‐dimensional modelling and inversion of ground‐based transient electromagnetic data at an area of stressed geological layers. The geophysical survey combined 58 transient electromagnetic sites to produce geoelectrical structures at different depth slices, and cross sections were used to characterise the horizon of the most likely sites for mining oil shale. The results show valuable information on the thickness of the oil shale seam at 3.7 Ωm, which is correlated to the geoelectrical layer between 2‐ and 4 ms transient time delays, and at depths ranging between 85 and 105 m. The 300 m penetrated depth of the transient electromagnetic soundings allows the resolution of the main geological units at narrow resistivity contrast and the distinction of the main geological structures that constrain the detection of the oil shale seam. This geoelectrical layer at different depth slices illustrates a localised oil shale setting and can be spatially correlated with an area bounded by fold and fault systems. Also, three‐dimensional modelling and inversion for synthetic and experimental data are introduced at the faulted area. The results show the limitations of oil shale imaging at a depth exceeding 130 m, which depends on the near‐surface resistivity layer, the low resistivity contrast of the main lithological units, and the degree of geological detail achieved at a suitable model's misfit value.  相似文献   

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

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