首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 500 毫秒
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.
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.  相似文献   

3.
欧洋  高文利  李洋  王宇航 《地球物理学报》2019,62(10):3843-3853
为了避免使用不合理初始辐射场强和方向性因子带来的误差,研究了估计辐射参数的井间电磁波层析成像技术.通过时域有限差分法模拟表明,天线长度与波长的比值、钻孔充填情况、钻孔周围介质的物性均会影响偶极天线的初始辐射场强或方向性因子;为此结合已知的分层资料,将它们设为未知参数,并设定初始辐射场强与发射点位置相关,方向性因子随射线角度而变化;采用正则化反演方法,由钻孔资料建立了模型方差目标函数,使得反演结果与钻遇的地质特征保持一致.通过理论模型试验和实例应用分析表明,相对于传统射线层析成像方法,估计辐射参数的正则化层析成像技术有助于提高反演的准确性.  相似文献   

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

5.
Time‐domain electromagnetic data are conveniently inverted by using smoothly varying 1D models with fixed vertical discretization. The vertical smoothness of the obtained models stems from the application of Occam‐type regularization constraints, which are meant to address the ill‐posedness of the problem. An important side effect of such regularization, however, is that horizontal layer boundaries can no longer be accurately reproduced as the model is required to be smooth. This issue can be overcome by inverting for fewer layers with variable thicknesses; nevertheless, to decide on a particular and constant number of layers for the parameterization of a large survey inversion can be equally problematic. Here, we present a focusing regularization technique to obtain the best of both methodologies. The new focusing approach allows for accurate reconstruction of resistivity distributions using a fixed vertical discretization while preserving the capability to reproduce horizontal boundaries. The formulation is flexible and can be coupled with traditional lateral/spatial smoothness constraints in order to resolve interfaces in stratified soils with no additional hypothesis about the number of layers. The method relies on minimizing the number of layers of non‐vanishing resistivity gradient, instead of minimizing the norm of the model variation itself. This approach ensures that the results are consistent with the measured data while favouring, at the same time, the retrieval of horizontal abrupt changes. In addition, the focusing regularization can also be applied in the horizontal direction in order to promote the reconstruction of lateral boundaries such as faults. We present the theoretical framework of our regularization methodology and illustrate its capabilities by means of both synthetic and field data sets. We further demonstrate how the concept has been integrated in our existing spatially constrained inversion formalism and show its application to large‐scale time‐domain electromagnetic data inversions.  相似文献   

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

7.
The spatial distribution of the electrical resistivity data provides useful information for investigating and modeling the fluid transport processes. 3D electrical resistivity distribution provides information about water flow and changes in electrical resistivity of the pore fluid.Therefore, to assist in understanding and modeling of the fluid transport process, 3D spatial distribution of the electrical resistivity data with the corresponded 3D geological section were mapped and interpreted in the test site located in western Germany. A process of deriving electrical resistivity values from the mechanical and radioactive parameters of cone penetration tests (CPT) and geological information of boreholes was presented. A reliable method which gives accurate resistivity values in cases of near surface sediments was introduced. Then a field test was executed where the calculated resistivity values were compared with the measured CPTe resistivity data. The CPTe (cone penetration test with electrical extension) data were also used in correlating to the ERT (electrical resistivity tomography) data. Consequently, obtained dense CPT surveys give us the possibility to determine a high resolution resistivity distribution of the investigated area.  相似文献   

8.
Model Fusion and Joint Inversion   总被引：1，自引：1，他引：0  
Inverse problems are inherently non-unique, and regularization is needed to obtain stable and reasonable solutions. The regularization adds information to the problem and determines which solution, out of the infinitely many, is obtained. In this paper, we review and discuss the case when a priori information exists in the form of either known structure or in the form of another inverse problem for a different property. The challenge is to include such information in the inversion process. To use existing known structure, we review the concept of model fusion, where we build a regularization functional that fuses the inverted model to a known one. The fusion is achieved by four different techniques. Joint inversion of two data sets is achieved by using iterative data fusion. The paper discusses four different methods for joint inversion. We discuss the use of correspondence maps or the petrophysics of the rocks, as well as structure. In particular, we suggest to further stabilize the well-known gradient cross product and suggest a new technique, Joint Total Variation, to solve the problem. The Joint Total Variation is a convex functional for joint inversion and, as such, has favorable optimization properties. We experiment with the techniques on the DC resistivity problem and the borehole tomography and show how model fusion and joint inversion can significantly improve over existing techniques.  相似文献   

9.
In this paper, evidence is presented that the combination of geospectral images and geophysical signatures (resistivity–velocity cross-plots) is a good tool to provide a natural visualization of the distribution and variations of lithological features in a test site. This was confirmed by the correlation between the electrical resistivity and seismic velocity values obtained after cross-gradient joint inversion at two profiles and geotechnical information provided by shallow boreholes in a site located in the Earth Sciences School grounds in Linares, Northeastern Mexico. The results obtained from this study show how the cross-gradient joint inversion facilitates the analysis of hydrological estimates and assists in lithological classification of subsurface materials.  相似文献   

10.
In order to gain a better understanding of the geometry of surface faults, five Controlled Source/Radio Magnetotelluric (CSRMT) profiles were measured across the Volvi basin, 45 km northeast of the city of Thessaloniki in Greece. The data were collected in two frequency ranges: a) 1–12.5 kHz using a remotely controlled double horizontal magnetic dipole transmitter (CSAMT measurements), and b) 15–250 kHz using the signal from distant radio transmitters (RMT measurements). The transition from the RMT band to the CSAMT band was smooth and continuous allowing us to combine both datasets for plane-wave modeling. The surface geology shows a predominantly 2D structure, and therefore we planned the survey into profiles perpendicular to the geological strike. We have used a 2D interpretation tool to model the data in TE, TM, TE + TM and determinant modes. Using a 4% error floor on the impedance, 2D resistivity models from inversion of the determinant data provide lower RMS data fits (4.2 and 1.2 for resistivity and phase, respectively) compared to the combined TE + TM data (4.4, 2.8, overall resistivity and phase, respectively). 2D inversion of the measured tensor data shows a sharp change in the depth to the top of resistive gneiss–schist basement that is overlain by a less resistive overburden at southern basin flanks. The change in depth to the bedrock is clearly seen in all 2D models along the measured profiles suggesting the existence of normal faults with strike directions of NE–SW to E–W. The 2D electrical resistivity models suggest that the bedrock deepens towards south-west. The resistivity models are also compared with the existing borehole information in the area and show a reasonable correlation. For example the sharp change of depth to the bedrock towards the center of the basin as seen in the resistivity models are also confirmed by the borehole data.  相似文献   

11.
无线电波透视法是常用的工作面地质构造探测方法之一,目前普遍使用的SIRT方法层析分辨率不高。本文采用约束正则化方法,推导Tikhonov正则化和全变差正则化的最小化问题表达式,讨论影响层析结果的主要因素,对典型理论模型进行了层析成像实验。结果表明:正则化方法具有比SIRT方法更好的分辨率;射线条数越多、噪声水平越低,层析分辨率越高;Tikhonov正则化在正则参数增大时层析结果更光滑,减小时则更贴近异常,全变差正则化与其相反。最后对实际坑透数据进行层析,识别出的异常构造基本吻合已知疑似构造位置,从而说明正则化方法在无线电波透视应用中的可行性。   相似文献   

12.
超高密度电法是一种新的地球物理探测技术,它通过多通道数据采集和多装置数据联合反演,极大地提高了电法勘探的成像精度.本文提出一种主成分-正则化极限学习机(PC-RELM)非线性反演方法,该方法针对超高密度电法所获取的高维勘探数据进行反演建模,通过随机设定隐层参数来简化模型的学习过程,通过主成分分析方法来进行高维数据降维,最后引入正则化因子提高反演模型的泛化能力.论文给出了超高密度电法的原理、样本构造方法和非线性反演流程,使用交叉验证方法获得了优化的隐节点数目和正则化参数,构造了优化的反演模型.通过两个经典的超高密度模型的反演结果表明,该方法能够较好地解决超高密度电法反演的高维数据非线性建模问题,能够弥补单一装置数据反演的不足,同时相较其他的非线性反演方法(ELM,BPNN和GRNN)具有更加准确的反演结果.  相似文献   

13.
The study describes a methodology used to integrate legacy resistivity data with limited geological data in order to build three-dimensional models of the near subsurface. Variogram analysis and inversion techniques more typically found in the petroleum industry are applied to a set of 1D resistivity data taken from electrical surveys conducted in the 1980s. Through careful integration with limited geological data collected from boreholes and outcrops, the resultant model can be visualized in three dimensions to depict alluvium layers as lithological and structural units within the bedrock. By tuning the variogram parameters to account for directionality, it is possible to visualize the individual lithofacies and geomorphological features in the subsurface. In this study, an electrical resistivity data set collected as part of a groundwater study in an area of the Peshawar basin in Pakistan has been re-examined. Additional lithological logs from boreholes throughout the area have been combined with local outcrop information to calibrate the data. Tectonic activity during the Himalayan orogeny has caused uplift in the area and generated significant faulting in the bedrock resulting in the formation of depressions which are identified by low resistivity values representing clays. Paleo-streams have reworked these clays which have been eroded and replaced by gravel–sand facies along paleo-channels. It is concluded that the sediments have been deposited as prograding fan-shaped bodies and lacustrine deposits with interlayered gravel–sand and clay–silt facies. The Naranji area aquifer system has thus been formed as a result of local tectonic activity with fluvial erosion and deposition and is characterized by coarse sediments with high electrical resistivities.  相似文献   

14.
根据非线性反演理论与Morozov偏差原理研究建立从双侧向测井(DLL)资料中同时重构地层原状电阻率、侵入带电阻率、侵入半径、层界面位置以及井眼泥浆电阻率的迭代正则化算法.首先利用Tikhonov正则化反演理论将双侧向测井资料的反演问题转化为含有稳定泛函的非线性目标函数的极小化问题,并利用Gauss-Newton算法确定极小化解.为得到稳定的反演结果并有效实现测井资料的最佳拟合,在迭代过程中将Morozov偏差原理和Cholesky分解技术相结合,建立了一套后验选择正则化因子的方法.最后通过理论模型和大庆油田实际测井资料的处理结果,验证了该算法能够取得更为满意的反演效果.  相似文献   

15.
Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.  相似文献   

16.
A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

17.
Magnetotelluric (MT) surveys were conducted in Pohang, Korea, for low-temperature geothermal exploration in 2002 and 2003. Pohang is located in the southeastern part of the Korean Peninsula and close to the East Sea. In the interpretation of MT data from a coastal environment, sea effects must be correctly included because seawater is a strong conductor. We first constructed a five-layered earth model with a realistic coastline and bathymetry to investigate sea effects on MT data measured in Pohang. This model clearly shows that the Pohang data are significantly influenced by sea water at frequencies blow 1 Hz at the whole measurement sites. Next, we utilized a three-dimensional inversion algorithm based on the Gauss–Newton approach to produce a reliable resistivity model. Seawater is excluded from the inversion domain to fix the resistivity, while included in the modeling domain to simulate sea effects on MT responses. Blocks for the sub-seafloor are included in unknown parameters since they are sufficiently close to the survey area to affect MT responses in Pohang. Static shifts are also considered in inversion for more accurate interpretation. The rms data misfit is smoothly reduced from 11.2 to 1.87 after 7 iterations. The resulting resistivity model shows a pattern of low–high–low resistivity with depth. The model is compatible with resistivity logs obtained from four boreholes in the survey area, and can explain major geological features in Pohang.  相似文献   

18.
Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.  相似文献   

19.
In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.  相似文献   

20.
地球物理反演是获取地球信息的重要手段,其求解具有严重的不适定性.为获得稳定的反问题结果,通常需要在目标泛函中加入正则化约束项.正确地估计正则化参数一直是地球物理反问题中的难点.目前存在的选取方法需要根据大量的试验来确定正则化参数,工作量十分巨大,并且存在很大的经验性,很难得到最优的正则化参数.针对这个问题,本文提出了一种基于广义Stein无偏风险估计的正则化参数求取方法.该方法的具体思路是通过求解模型参数均方误差的广义Stein无偏风险估计函数,在反问题求解过程中自动求取正则化参数.本文模型测试结果表明,相比于目前常用的方法,通过该方法得到的正则化参数是最优的.  相似文献   

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

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