首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.  相似文献   

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

3.
In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.  相似文献   

4.
We present a simple and feasible approach to analyse and identify two‐dimensional effects in central loop transient electromagnetic sounding data and the correspondingly derived quasi two‐dimensional conductivity models. The proposed strategy is particularly useful in minimising interpretation errors. It is based on the calculation of a semi‐synthetic transient electromagnetic tipper at each sounding and for each observational transient time point. The semi‐synthetic transient electromagnetic tipper is derived from the measured vertical component of the induced voltage and the synthetically calculated horizontal component. The approach is computationally inexpensive and involves one two‐dimensional forward calculation of an obtained quasi two‐dimensional conductivity section. Based on a synthetic example, we demonstrate that the transient electromagnetic tipper approach is applicable in identifying which transient data points and which corresponding zones in a derived quasi two‐dimensional subsurface model are affected by two‐dimensional inhomogeneities. The one‐dimensional inversion of such data leads to false models. An application of the semi‐synthetic transient electromagnetic tipper to field data from the Azraq basin in Jordan reveals that, in total, eight of 80 investigated soundings are affected by two‐dimensional structures although the field data can be fitted optimally using one‐dimensional inversion techniques. The largest semi‐synthetic tipper response occurs in a 300 m‐wide region around a strong lateral resistivity contrast. The approach is useful for analysing structural features in derived quasi two‐dimensional sections and for qualitatively investigating how these features affect the transient response. To avoid misinterpretation, these identified zones corresponding to large tipper values are excluded from the interpretation of a quasi two‐dimensional conductivity model. Based on the semi‐synthetic study, we also demonstrate that a quantitative interpretation of the horizontal voltage response (e.g. by inversion) is usually not feasible as it requires the exact sensor position to be known. Although a tipper derived purely from field data is useful as a qualitative tool for identifying two‐dimensional distortion effects, it is only feasible if the sensor setup is sufficiently accurate. Our proposed semi‐synthetic transient electromagnetic tipper approach is particularly feasible as an a posteriori approach if no horizontal components are recorded or if the sensor setup in the field is not sufficiently accurate.  相似文献   

5.
核磁共振与瞬变电磁三维联合解释方法   总被引：3，自引：1，他引：2       下载免费PDF全文
传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用.  相似文献   

6.
The accurate estimation of sub‐seafloor resistivity features from marine controlled source electromagnetic data using inverse modelling is hindered due to the limitations of the inversion routines. The most commonly used one‐dimensional inversion techniques for resolving subsurface resistivity structures are gradient‐based methods, namely Occam and Marquardt. The first approach relies on the smoothness of the model and is recommended when there are no sharp resistivity boundaries. The Marquardt routine is relevant for many electromagnetic applications with sharp resistivity contrasts but subject to the appropriate choice of a starting model. In this paper, we explore the ability of different 1D inversion schemes to derive sub‐seafloor resistivity structures from time domain marine controlled source electromagnetic data measured along an 8‐km‐long profile in the German North Sea. Seismic reflection data reveal a dipping shallow amplitude anomaly that was the target of the controleld source electromagnetic survey. We tested four inversion schemes to find suitable starting models for the final Marquardt inversion. In this respect, as a first scenario, Occam inversion results are considered a starting model for the subsequent Marquardt inversion (Occam–Marquardt). As a second scenario, we employ a global method called Differential Evolution Adaptive Metropolis and sequentially incorporate it with Marquardt inversion. The third approach corresponds to Marquardt inversion introducing lateral constraints. Finally, we include the lateral constraints in Differential Evolution Adaptive Metropolis optimization, and the results are sequentially utilized by Marquardt inversion. Occam–Marquardt may provide accurate estimation of the subsurface features, but it is dependent on the appropriate conversion of different multi‐layered Occam model to an acceptable starting model for Marquardt inversion, which is not straightforward. Employing parameter spaces, the Differential Evolution Adaptive Metropolis approach can be pertinent to determine Marquardt a priori information; nevertheless, the uncertainties in Differential Evolution Adaptive Metropolis optimization will introduce some inaccuracies in Marquardt inversion results. Laterally constrained Marquardt may be promising to resolve sub‐seafloor features, but it is not stable if there are significant lateral changes of the sub‐seafloor structure due to the dependence of the method to the starting model. Including the lateral constraints in Differential Evolution Adaptive Metropolis approach allows for faster convergence of the routine with consistent results, furnishing more accurate estimation of a priori models for the subsequent Marquardt inversion.  相似文献   

7.
A comprehensive controlled source electromagnetic (CSEM) modelling study, based on complex resistivity structures in a deep marine geological setting, is conducted. The study demonstrates the effects of acquisition parameters and multi‐layered resistors on CSEM responses. Three‐dimensional (3D) finite difference time domain (FDTD) grid‐modelling is used for CSEM sensitivity analysis. Interpolation of normalized CSEM responses provides attributes representing relative sensitivity of the modelled structures. Modelling results show that fine grid, 1 × 1 km receiver spacing, provides good correlations between CSEM responses and the modelled structures, irrespective of source orientation. The resolution of CSEM attributes decreases for receiver spacing >2 × 2 km, when using only in‐line data. Broadside data in the grid geometry increase data density by 100 – approximately 200% by filling in in‐line responses and improve the resolution of CSEM attributes. Optimized source orientation (i.e., oblique to the strike of an elongated resistor) improves the structural definition of the CSEM anomalies for coarse‐grid geometries (receiver spacing ≥3 × 3 km). The study also shows that a multi‐resistor anomaly is not simply the summation but a cumulative response with mutual interference between constituent resistors. The combined response of constituent resistors is approximately 50% higher than the cumulative response of the multi‐resistor for 0.5 Hz at 4000 m offset. A gradual inverse variation of offset and frequency allows differentiation of CSEM anomalies for multi‐layered resistors. Similar frequency‐offset variations for laterally persistent high‐resistivity facies show visual continuity with varying geometric expressions. 3D grid‐modelling is an effective and adequate tool for CSEM survey design and sensitivity analysis.  相似文献   

8.
To advance and optimize secondary and tertiary oil recovery techniques, it is essential to know the areal propagation and distribution of the injected fluids in the subsurface. We investigate the applicability of controlled‐source electromagnetic methods to monitor fluid movements in a German oilfield (Bockstedt, onshore Northwest Germany) as injected brines (highly saline formation water) have much lower electrical resistivity than the oil within the reservoir. The main focus of this study is on controlled‐source electromagnetic simulations to test the sensitivity of various source–receiver configurations. The background model for the simulations is based on two‐dimensional inversion of magnetotelluric data gathered across the oil field and calibrated with resistivity logs. Three‐dimensional modelling results suggest that controlled‐source electromagnetic methods are sensitive to resistivity changes at reservoir depths, but the effect is difficult to resolve with surface measurements only. Resolution increases significantly if sensors or transmitters can be placed in observation wells closer to the reservoir. In particular, observation of the vertical electric field component in shallow boreholes and/or use of source configurations consisting of combinations of vertical and horizontal dipoles are promising. Preliminary results from a borehole‐to‐surface controlled‐source electromagnetic field survey carried out in spring 2014 are in good agreement with the modelling studies.  相似文献   

9.
We present a new workflow for imaging damped three‐dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short‐offset and low‐frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine‐scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite‐difference frequency‐domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full‐waveform inversion workflow is applied to the 3D Marmousi‐2 and SEG/EAGE Salt models with long‐offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one‐third of the maximum offset length at a resolution consistent with the inverted frequencies.  相似文献   

10.
In order to couple spatial data from frequency‐domain helicopter‐borne electromagnetics with electromagnetic measurements from ground geophysics (transient electromagnetics and radiomagnetotellurics), a common 1D weighted joint inversion algorithm for helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics data has been developed. The depth of investigation of helicopter‐borne electromagnetics data is rather limited compared to time‐domain electromagnetics sounding methods on the ground. In order to improve the accuracy of model parameters of shallow depth as well as of greater depth, the helicopter‐borne electromagnetics, transient electromagnetics, and radiomagnetotellurics measurements can be combined by using a joint inversion methodology. The 1D joint inversion algorithm is tested for synthetic data of helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics. The proposed concept of the joint inversion takes advantage of each method, thus providing the capability to resolve near surface (radiomagnetotellurics) and deeper electrical conductivity structures (transient electromagnetics) in combination with valuable spatial information (helicopter‐borne electromagnetics). Furthermore, the joint inversion has been applied on the field data (helicopter‐borne electromagnetics and transient electromagnetics) measured in the Cuxhaven area, Germany. In order to avoid the lessening of the resolution capacities of one data type, and thus balancing the use of inherent and ideally complementary information content, a parameter reweighting scheme that is based on the exploration depth ranges of the specific methods is proposed. A comparison of the conventional joint inversion algorithm, proposed by Jupp and Vozoff ( 1975 ), and of the newly developed algorithm is presented. The new algorithm employs the weighting on different model parameters differently. It is inferred from the synthetic and field data examples that the weighted joint inversion is more successful in explaining the subsurface than the classical joint inversion approach. In addition to this, the data fittings in weighted joint inversion are also improved.  相似文献   

11.
A detailed magnetotelluric survey was conducted in 2013 in the Sehqanat oil field, southwestern Iran to map the geoelectrical structures of the sedimentary Zagros zone, particularly the boundary between the Gachsaran Formation acting as cap rock and the Asmari Formation as the reservoir. According to the electrical well logs, a large resistivity contrast exists between the two formations. The Gachsaran Formation is formed by tens to hundreds of metres of evaporites and it is highly conductive (ca. 1 Ωm–10 Ωm), and the Asmari Formation consists of dense carbonates, which are considerably more resistive (more than 100 Ωm). Broadband magnetotelluric data were collected along five southwest–northeast directed parallel lines with more than 600 stations crossing the main geological trend. Although dimensionality and strike analysis of the magnetotelluric transfer functions showed that overall they satisfied local 2D conditions, there were also strong 3D conditions found in some of the sites. Therefore, in order to obtain a more reliable image of the resistivity distribution in the Sehqanat oil field, in addition to standard 2D inversion, we investigated to what extent 3D inversion of the data was feasible and what improvements in the resistivity image could be obtained. The 2D inversion models using the determinant average of the impedance tensor depict the main resistivity structures well, whereas the estimated 3D model shows significantly more details although problems were encountered in fitting the data with the latter. Both approaches resolved the Gachsaran–Asmari transition from high conductivity to moderate conductivity. The well‐known Sehqanat anticline could also be delineated throughout the 2D and 3D resistivity models as a resistive dome‐shaped body in the middle parts of the magnetotelluric profiles.  相似文献   

12.
This paper presents the first controlled‐source electromagnetic survey carried out in the German North Sea with a recently developed seafloor‐towed electrical dipole–dipole system, i.e., HYDRA II. Controlled‐source electromagnetic data are measured, processed, and inverted in the time domain to estimate an electrical resistivity model of the sub‐seafloor. The controlled‐source electromagnetic survey targeted a shallow, phase‐reversed, seismic reflector, which potentially indicates free gas. To compare the resistivity model to reflection seismic data and draw a combined interpretation, we apply a trans‐dimensional Bayesian inversion that estimates model parameters and uncertainties, and samples probabilistically over the number of layers of the resistivity model. The controlled‐source electromagnetic data errors show time‐varying correlations, and we therefore apply a non‐Toeplitz data covariance matrix in the inversion that is estimated from residual analysis. The geological interpretation drawn from controlled‐source electromagnetic inversion results and borehole and reflection seismic data yield resistivities of ～1 Ωm at the seafloor, which are typical for fine‐grained marine deposits, whereas resistivities below ～20 mbsf increase to 2–4 Ωm and can be related to a transition from fine‐grained (Holocene age) to unsorted, coarse‐grained, and compacted glacial sediments (Pleistocene age). Interface depths from controlled‐source electromagnetic inversion generally match the seismic reflector related to the contrast between the different depositional environments. Resistivities decrease again at greater depths to ～1 Ωm with a minimum resistivity at ～300 mbsf where a seismic reflector (that marks a major flooding surface of late Miocene age) correlates with an increased gamma‐ray count, indicating an increased amount of fine‐grained sediments. We suggest that the grain size may have a major impact on the electrical resistivity of the sediment with lower resistivities for fine‐grained sediments. Concerning the phase‐reversed seismic reflector that was targeted by the survey, controlled‐source electromagnetic inversion results yield no indication for free gas below it as resistivities are generally elevated above the reflector. We suggest that the elevated resistivities are caused by an overall decrease in porosity in the glacial sediments and that the seismic reflector could be caused by an impedance contrast at a thin low‐velocity layer. Controlled‐source electromagnetic interface depths near the reflector are quite uncertain and variable. We conclude that the seismic interface cannot be resolved with the controlled‐source electromagnetic data, but the thickness of the corresponding resistive layer follows the trend of the reflector that is inclined towards the west.  相似文献   

13.
The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section.  相似文献   

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

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

16.
We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO2 sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO2 injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO2 dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO2, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO2 plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

17.
Z‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.  相似文献   

18.
This paper discusses the asymptotic behaviour of the electromagnetic fields received on the sea‐bed (target response), as well as the fields distributed inside a thin resistive target, generated by a horizontal electric dipole above the sea‐bed in marine controlled‐source electromagnetics for hydrocarbon exploration. It is found that the guided wave supported by a thin resistive target can be expressed as a single‐mode exponential function. A simple closed‐form expression is derived to relate the single‐mode wavenumber of the guided wave to the model parameters: the resistivity and thickness of the target layer, the sea‐bed resistivity and the frequency. When the air‐wave is removed, the guided wave is dominant among the fields received on the sea‐bed at far offset. Hence the wavenumber of the guided wave can be calculated from the fields measured on the sea‐bed. The closed‐form expression can then be used to invert the target property from the calculated wavenumber and hence, can be considered as a hydrocarbon indicator.  相似文献   

19.
We present results of synthetic time‐lapse and real repeatability multi‐transient electromagnetic surveys over the North Sea Harding field. Using Archie's law to convert porosity and fluid saturation to resistivity we created 3D isotropic models of the reservoir resistivity at different stages of production from the initial state in 1996 through to complete hydrocarbon production by 2016 and, for each stage, we simulated an east‐west transient electromagnetic survey line across Harding. Unconstrained 1D full‐waveform Occam inversions of these synthetic data show that Harding should be detectable and its lateral extent reasonably well‐defined. Resistivity changes caused by hydrocarbon production from initial pre‐production state to production of the oil rim in 2011 are discernible as are significant changes from 2011–2016 during the modelled gas blowdown phase. The 2D repeatability surveys of 2007 and 2008 tied two wells: one on and the other off the structure. Between the two surveys the segment of the field under investigation produced 3.9 million barrels of oil – not enough to generate an observable time‐lapse electromagnetic anomaly with a signal‐to‐noise ratio of 40 dB. Processing of the 2007 and 2008 data included deconvolution for the measured source current and removal of spatially‐correlated noise, which increased the signal‐to‐noise ratio of the recovered impulse responses by about 20 dB and resulted in a normalized root‐mean‐square difference of 3.9% between the data sets. 1D full‐waveform Occam inversions of the real data showed that Harding was detectable and its lateral extent was also reasonably well‐defined. The results indicate that the multi‐transient electromagnetic method is suitable for exploration, appraisal and monitoring hydrocarbon production.  相似文献   

20.
In the last two decades, forward modelling for the time domain (transient) electromagnetic method has concentrated almost entirely on multi-dimensional models and algorithms. At the same time, the interpretation of real field data is still mainly one dimensional. This is caused by the lack of an efficient multi-dimensional acquisition procedure supported by sufficiently fast and reliable inversion software, on the one hand, and by the great efficiency of one-dimensional field set up and interpretation of the data on the other hand. The latter is particularly true for the short offset transient electromagnetic method, which is much less sensitive to multi-dimensional effects, compared to long offset methods. The most commonly used one-dimensional forward modelling algorithms are based on the spectral method, which requires calculating rapidly oscillating Fourier–Bessel (Hankel) integrals. Due to the very fast decay of short offset responses, the integrals become computationally unstable at late times of the transient process. Although this problem has been successfully solved for practically feasible measurement times of conventional short offset systems using transverse electric and mixed transverse electric and transverse magnetic fields, it turned out crucial for novel methods based on the use of unimodal transverse magnetic fields. These methods are much more sensitive to geoelectric parameters of the Earth in general and those of resistive targets, in particular, but they generate responses, which drop at late times significantly faster than those of conventional methods. Such behaviour of transverse magnetic fields represents severe computational problem for the spectral method, but is successfully solved by direct time domain algorithms. This article describes a generalization of the well-known Tikhonov's solution to a boundary value problem directly in time domain, which is applied to an arbitrary one-dimensional earth model excited by an arbitrary source. Contrary to existing spectral algorithms, the described method allows accurate calculations of both transverse electric and transverse magnetic transient responses at arbitrarily late times. On the other hand, it is more time efficient than finite-difference/finite element direct time domain algorithms and provides analytical late-stage asymptotic solutions.  相似文献   

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

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