共查询到20条相似文献，搜索用时 447 毫秒
1.
The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity. 相似文献
2.
Michele Cercato 《Geophysical Prospecting》2011,59(2):210-226
Practical applications of surface wave inversion demand reliable inverted shear‐wave profiles and a rigorous assessment of the uncertainty associated to the inverted parameters. As a matter of fact, the surface wave inverse problem is severely affected by solution non‐uniqueness: the degree of non‐uniqueness is closely related to the complexity of the observed dispersion pattern and to the experimental inaccuracies in dispersion measurements. Moreover, inversion pitfalls may be connected to specific problems such as inadequate model parametrization and incorrect identification of the surface wave modes. Consequently, it is essential to tune the inversion problem to the specific dataset under examination to avoid unnecessary computations and possible misinterpretations. In the heuristic inversion algorithm presented in this paper, different types of model constraints can be easily introduced to bias constructively the solution towards realistic estimates of the 1D shear‐wave profile. This approach merges the advantages of global inversion, like the extended exploration of the parameter space and a theoretically rigorous assessment of the uncertainties on the inverted parameters, with the practical approach of Lagrange multipliers, which is often used in deterministic inversion, which helps inversion to converge towards models with desired properties (e.g., ‘smooth’ or ‘minimum norm' models). In addition, two different forward kernels can be alternatively selected for direct‐problem computations: either the conventional modal inversion or, instead, the direct minimization of the secular function, which allows the interpreter to avoid mode identification. A rigorous uncertainty assessment of the model parameters is performed by posterior covariance analysis on the accepted solutions and the modal superposition associated to the inverted models is investigated by full‐waveform modelling. This way, the interpreter has several tools to address the more probable sources of inversion pitfalls within the framework of a rigorous and well‐tested global inversion algorithm. The effectiveness and the versatility of this approach, as well as the impact of the interpreter's choices on the final solution and on its posterior uncertainty, are illustrated using both synthetic and real data. In the latter case, the inverted shear velocity profiles are blind compared with borehole data. 相似文献
3.
Numerical investigations have been carried out on suitable interpretation schemes for gravity, based upon Backus-Gilbert inverse theory. A two-dimensional linear model consisting of horizontal prisms (grid shape) has been adopted for interpretation of gravity profiles. An inversion technique, which uses a linear approach suitable for rapid computation, and studies with synthetic bodies have led to improvements in the method, mainly by applying weighting factors, which improve the solution in the final stages of computational work. Gravity data from the well-known Cyprus island positive anomaly were inverted to demonstrate the feasibility of the method. 相似文献
4.
Due to its simplicity, stability, and efficiency, the use of right rectangular prisms is still widespread for potential field modelling and inversion. It is well known that modelling the subsurface with Cartesian grids has important consequences in terms of accuracy of the results. In this paper, we review the main issues that geophysicists face in day‐to‐day work when trying to use right rectangular prisms for performing gravity or full tensor gravity modelling and inversions. We demonstrate the results both theoretically and through Monte Carlo simulations, also exploiting concepts from fractal geometry. We believe that the guidelines contained in this paper may suggest a good practice for the day‐to‐day work of geophysicists dealing with gravity and full tensor gravity data. 相似文献
5.
裂陷盆地基底的起伏表现为非光滑的几何形态,传统的重力反演结果并不能很好地反映这种特点.此外,大多数情况下,重力观测面并不位于盆地上界面,应为单独的起伏观测面,盆地应为上界面和基底组成的双界面模式.基于此,本文研究了起伏观测面上裂陷盆地基底双界面模式二维重力反演方法.研究中假设沉积盆地的沉积层与基底的密度差随深度按双曲线规律变化.将沉积盆地的沉积层剖分成相邻的垂直柱体,其水平尺寸是已知的,顶面与沉积层上界面重合,底面深度代表基底的深度,即为要反演的参数.反演中引入全变差函数作为盆地模型的约束,使得反演结果呈现非光滑形态,符合裂陷盆地基底特征.为减小反演多解性,引入已知深度点作为约束.建立由重力数据拟合、已知深度约束及全变差函数组成的目标函数,采用非线性共轭梯度算法使目标函数最小化.模型试算结果表明该方法可反演裂陷盆地基底起伏,并通过调整正则化参数的值可反演坳陷盆地基底起伏.将该反演方法用于珠江口盆地惠州凹陷和运城-临汾裂陷盆地实际资料处理,其结果较好地反映了裂陷盆地基底起伏特征,为研究盆地构造、油气勘探等提供重要参考. 相似文献
6.
以全球地震面波相速度变化问题为例，使用约３００００道高质量的面波记录数据集，在数据空间和模型空间的误差服从零平均Ｇａｕｓｓ随机分布的假设下，研究反演中的阻尼和参量化问题．发现最佳分块尺度随频率增高而减小；分块参量化方法比球谐展开方法引入的人为因素少，但难以分辨甚低阶横向非均匀性，除非施加极端的空间平滑；球谐函数展开则具备恢复长波结构的优点．为此提出一种混合参量反演方法：首先用球谐函数作为全球基函数，恢复相速度的低阶球谐分量，然后用此作为进一步反演的初始模型，用分块模型迭代反演，得到最终结果． 相似文献
7.
Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models 
下载免费PDF全文
下载免费PDF全文 In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information. 相似文献
8.
Tijmen Jan Moser 《Geophysical Prospecting》2014,62(3):475-490
We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots. 相似文献
9.
We have developed a method for imaging magnetic data collected for mineral exploration to yield the following structural information: depth, model type (structural index) and susceptibility. The active nature of mineral exploration data requires we derive the structural information from a robust quantity: we propose that the first‐ or second‐order analytic‐signal amplitude is suitably stable. The procedure is to normalize the analytic‐signal amplitude by the peak value and then use non‐linear inversion to estimate the depth and the structural index for each anomaly. In our field example, different results are obtained depending on whether we inverted for the first‐ or second‐order analytic‐signal amplitude. This is probably because the two‐dimensional contact, thin sheet or horizontal cylinder models we have assumed are not appropriate. In cases such as these, when our model assumptions are not correct, the results should not be interpreted quantitatively, but they might be useful for giving a qualitative indication of how the structure might vary. With a priori information, it is possible to assume a model type (i.e. set the structural index) and generate estimates of the depth and susceptibility. These data can then be gridded and imaged. If a contact is assumed, the susceptibility contrast is estimated; for the dike model, the susceptibility‐thickness is estimated; for the horizontal cylinder, the susceptibility‐area is estimated. To emphasize that the results are dependent on our assumed model, we advocate prefixing any derived quantity by the term ‘apparent’. 相似文献
10.
In this paper a procedure for analysing the seismic response of seabed anchored floating tunnels is presented. The first step of the research was the development of an ‘ad hoc’ finite element for modelling the behaviour of anchor elements, with particular reference to the problem of transverse oscillations under time varying axial loads. The element was subsequently inserted in a step‐by‐step procedure for the numerical analysis of non‐linear response to multiple‐support seismic input; the procedure encompasses simplified modelling of fluid–structure and soil–structure interaction effects. An example of an application is given concerning two 4680 m long floating tunnels with different seabed profiles. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
11.
In this paper we propose a 3D acoustic full waveform inversion algorithm in the Laplace domain. The partial differential equation for the 3D acoustic wave equation in the Laplace domain is reformulated as a linear system of algebraic equations using the finite element method and the resulting linear system is solved by a preconditioned conjugate gradient method. The numerical solutions obtained by our modelling algorithm are verified through a comparison with the corresponding analytical solutions and the appropriate dispersion analysis. In the Laplace‐domain waveform inversion, the logarithm of the Laplace transformed wavefields mainly contains long‐wavelength information about the underlying velocity model. As a result, the algorithm smoothes a small‐scale structure but roughly identifies large‐scale features within a certain depth determined by the range of offsets and Laplace damping constants employed. Our algorithm thus provides a useful complementary process to time‐ or frequency‐domain waveform inversion, which cannot recover a large‐scale structure when low‐frequency signals are weak or absent. The algorithm is demonstrated on a synthetic example: the SEG/EAGE 3D salt‐dome model. The numerical test is limited to a Laplace‐domain synthetic data set for the inversion. In order to verify the usefulness of the inverted velocity model, we perform the 3D reverse time migration. The migration results show that our inversion results can be used as an initial model for the subsequent high‐resolution waveform inversion. Further studies are needed to perform the inversion using time‐domain synthetic data with noise or real data, thereby investigating robustness to noise. 相似文献
12.
1D joint multi‐offset inversion of time‐domain marine controlled source electromagnetic data 下载免费PDF全文
下载免费PDF全文 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. 相似文献
13.
A nonlinear method for multiparameter inversion of pre‐stack seismic data based on anisotropic Markov random field 下载免费PDF全文
下载免费PDF全文 Hongbing Zhang Qiang Guo Lifeng Liang Chenghao Cao Zuoping Shang 《Geophysical Prospecting》2018,66(3):461-477
Multiparameter inversion for pre‐stack seismic data plays a significant role in quantitative estimation of subsurface petrophysical properties. However, it remains a complicated problem due to the non‐unique results and unstable nature of the processing; the pre‐stack seismic inversion problem is ill‐posed and band‐limited. Combining the full Zoeppritz equation and additional assumptions with edge‐preserving regularisation can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function that includes edge‐preserving regularisation and soft constraints based on anisotropic Markov random fields and is intended especially for layered formations. We applied a fast simulated annealing algorithm to solve the nonlinear optimisation problem. The method directly obtains reflectivity RPP values using the full Zoeppritz equation instead of its approximations and effectively controls the stability of the multiparameter inversion by assuming a sectionally constant S‐ and P‐wave velocity ratio and using the generalised Gardner equation. We substituted the inverted parameters, i.e., the P‐wave velocity, the fitting deviation of S‐wave velocity, and the density were inverted instead of the P‐wave velocity, the S‐wave velocity, and the density, and the generalised Gardner equation was applied as a constraint. Test results on two‐dimensional synthetic data indicated that our substitution obtained improved results for multiparameter inversion. The inverted results could be improved by utilising high‐order anisotropic Markov random field neighbourhoods at early stages and low‐order anisotropic Markov random field neighbourhoods in the later stages. Moreover, for layered formations, using a large horizontal weighting coefficient can preserve the lateral continuity of layers, and using a small vertical weighting coefficient allows for large longitudinal gradients of the interlayers. The inverted results of the field data revealed more detailed information about the layers and matched the logging curves at the wells acceptably over most parts of the curves. 相似文献
14.
A parallel,scalable and memory efficient inversion code for very large‐scale airborne electromagnetics surveys 下载免费PDF全文
下载免费PDF全文 Esben Auken 《Geophysical Prospecting》2015,63(2):495-507
Over the past decade the typical size of airborne electromagnetic data sets has been growing rapidly, along with an emerging need for highly accurate modelling. One‐dimensional approximate inversions or data transform techniques have previously been employed for very large‐scale studies of quasi‐layered settings but these techniques fail to provide the consistent accuracy needed by many modern applications such as aquifer and geological mapping, uranium exploration, oil sands and integrated modelling. In these cases the use of more time‐consuming 1D forward and inverse modelling provide the only acceptable solution that is also computationally feasible. When target structures are known to be quasi layered and spatially coherent it is beneficial to incorporate this assumption directly into the inversion. This implies inverting multiple soundings at a time in larger constrained problems, which allows for resolving geological layers that are undetectable using simple independent inversions. Ideally, entire surveys should be inverted at a time in huge constrained problems but poor scaling properties of the underlying algorithms typically make this challenging. Here, we document how we optimized an inversion code for very large‐scale constrained airborne electromagnetic problems. Most importantly, we describe how we solve linear systems using an iterative method that scales linearly with the size of the data set in terms of both solution time and memory consumption. We also describe how we parallelized the core region of the code, in order to obtain almost ideal strong parallel scaling on current 4‐socket shared memory computers. We further show how model parameter uncertainty estimates can be efficiently obtained in linear time and we demonstrate the capabilities of the full implementation by inverting a 3327 line km SkyTEM survey overnight. Performance and scaling properties are discussed based on the timings of the field example and we describe the criteria that must be fulfilled in order to adapt our methodology for similar type problems. 相似文献
15.
本文对研究区(98°E～103°E, 315°N～34°N) 1∶50万航磁资料的476条南北向剖面采用磁性标量层析成像方法进行了2-D成像反演.每条剖面资料点距1 km,反演的磁性柱体剖面方向宽度1 km,共约277个柱体,深度方向共8层柱体,每层厚度4 km.反演的476个断面的磁性柱体相对磁化强度幅值组成了3-D磁性数据体.对5.5、17.5、21.5、29.5 km深的磁性水平切片和7个北东向的磁性垂直切片表现的磁性特征进行了研究,并进行了地质解释.结合研究区已经公布的地质、地球物理资料,讨论了磁性细结构能提供工程地质灾害数值模拟仿真结构模型的信息,得到了反演的磁性细结构特征比场的特征能提供更多的断层和岩体3.D细节信息的初步认识. 相似文献
16.
Bingqiang Yuan Lijun Song Li Han Shaole An Chunguan Zhang 《Geophysical Prospecting》2018,66(8):1586-1601
The Tobago Basin, which is located offshore northern Venezuela with a southern margin close to Trinidad and Tobago, has an area of approximately 59,600 km2. The Tobago Basin has relatively favourable hydrocarbon prospects, and to date, exploration work has mainly concentrated on small areas of the southwestern portion of the basin. To conduct a comprehensive study of the structural framework of the basin and the characteristics of the basement in order to identify prospective zones for hydrocarbon exploration, shipborne‐measured and satellite‐measured gravity data, shipborne‐measured magnetic data, and aeromagnetic survey data were analysed. A regularisation filtering method was used to separate and obtain regional and residual gravity and magnetic anomalies. Directional gradients of gravity and magnetic anomalies and the total horizontal gradient and vertical second derivative of gravity anomalies were employed to extract information about fault structures. Regression analysis methods were used to determine the basement depth. The geological significance of the gravity and magnetic fields was examined, the structural framework of the basin was assessed, the basement depth was estimated, and favourable hydrocarbon exploration prospects within the basin were identified. The results show that the Tobago Basin contains complex structures consisting mainly of two groups of faults trending in northeasterly and northwesterly directions and that the major northeasterly trending faults control the main structural configuration and depositional system within the basin. The basement of the Tobago Basin has deep rises and falls. It can be divided into the following four secondary tectonic units: the western sub‐basin, the central uplift area, the southern sub‐basin, and the northeastern sub‐basin. The central uplift area and northeastern sub‐basin are most likely to have developed hydrocarbon accumulations and should be targeted for further exploration. 相似文献
17.
Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration 下载免费PDF全文
下载免费PDF全文 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. 相似文献
18.
Leonardo Azevedo Rúben Nunes Amílcar Soares Guenther Schwedersky Neto Teresa S. Martins 《Geophysical Prospecting》2018,66(Z1):116-131
Seismic inversion plays an important role in reservoir modelling and characterisation due to its potential for assessing the spatial distribution of the sub‐surface petro‐elastic properties. Seismic amplitude‐versus‐angle inversion methodologies allow to retrieve P‐wave and S‐wave velocities and density individually allowing a better characterisation of existing litho‐fluid facies. We present an iterative geostatistical seismic amplitude‐versus‐angle inversion algorithm that inverts pre‐stack seismic data, sorted by angle gather, directly for: density; P‐wave; and S‐wave velocity models. The proposed iterative geostatistical inverse procedure is based on the use of stochastic sequential simulation and co‐simulation algorithms as the perturbation technique of the model parametre space; and the use of a genetic algorithm as a global optimiser to make the simulated elastic models converge from iteration to iteration. All the elastic models simulated during the iterative procedure honour the marginal prior distributions of P‐wave velocity, S‐wave velocity and density estimated from the available well‐log data, and the corresponding joint distributions between density versus P‐wave velocity and P‐wave versus S‐wave velocity. We successfully tested and implemented the proposed inversion procedure on a pre‐stack synthetic dataset, built from a real reservoir, and on a real pre‐stack seismic dataset acquired over a deep‐water gas reservoir. In both cases the results show a good convergence between real and synthetic seismic and reliable high‐resolution elastic sub‐surface Earth models. 相似文献
19.
R.E. Chvez E.L. Flores J.O. Campos M. Ladrn de Guevara M.C. Fernndez‐Puga J. Herrera 《Geophysical Prospecting》2000,48(5):835-870
A comprehensive reinterpretation of the available gravity, magnetic, geothermal, geological and borehole information has been made of the Laguna Salada Basin to establish a 3D model of the basement and sedimentary infill. According to statistical spectral analysis, the residual gravity anomaly is due to sources with a mean regional depth of 2.8 km. The topography of the basement was obtained from a three‐dimensional inversion carried out in the wavenumber domain using an iterative scheme. The maximum density contrast of ?300 kg/m3 estimated from previous studies and the mean depth of 2.5 km finally constrained this inversion. The resulting model indicated that the sedimentary infill is up to 4.2 km thick at its deepest point. According to the gravity‐derived basement topography, the basin presents an asymmetry (i.e. it is of the half‐graben type). It is deeper to the east, where it is delimited from the Sierra Cucapah by a step fault. By contrast, the limit with the Sierra de Juarez is a gently sloping fault (i.e. a listric fault). The basement is not even, but it comprises a series of structural highs and lows. N–S to NW–SE and E–W to NE–SW faults delimit these structural units. The magnetic modelling was constrained by (i) the gravity‐derived basement topography; (ii) a Curie isotherm assumed to be between 7 km and 10 km; (iii) assuming induced magnetization only; (iv) the available geological and borehole information. The magnetic anomalies were interpreted successfully using the gravity‐derived basement/sedimentary interface as the top of the magnetic bodies (i.e. the magnetic modelling supports the gravity basement topography). An elongated N–S to NW–SE trending highly magnetized body running from south to north along the basin is observed to the west of the basin. This magnetic anomaly has no gravity signature. Such a feature can be interpreted as an intrusive body emplaced along a fault running through the Laguna Salada Basin. Treatment of the gravity and magnetic information (and of their horizontal gradients) with satellite image processing techniques highlighted lineaments on the basement gravity topography correlating with mapped faults. Based on all this information, we derived detailed geological models along four selected profiles to simulate numerically the heat and fluid flow in the basin. We used a finite‐difference scheme to solve the coupled Darcy and Fourier differential equations. According to our results, we have fluid flow in the sedimentary layers and a redistribution of heat flow from the basin axis toward its rims (Sierra de Juárez and Sierra Cucapah). Our model temperatures agree within an error of 4% with the observed temperature profiles measured at boreholes. Our heat‐flow determinations agree within an error of ±15% with extrapolated observations. The numerical and chemical analyses support the hypothesis of fluid circulation between the clay–lutite layer and the fractured granitic basement. Thermal modelling shows low heat‐flow values along the Laguna Salada Basin. Deep fluid circulation patterns were observed that redistribute such flow at depth. Two patterns were distinguished. One displays the heat flow increasing from the basin axis towards its borders (temperature increase of 20°C). The second pattern shows an increasing heat flow from south to north of the basin. Such behaviour is confirmed by the temperature measurements in the thermometric boreholes. 相似文献
20.
Numerical study of the variations of magnetic resonance signals caused by surface slope 总被引：1，自引：0，他引：1
When performing forward modelling and inversion of Magnetic Resonance Sounding (MRS) data, the water-content distribution is typically assumed to be horizontal (1D case). This assumption is fully justified because MRS is often used for characterizing continuous aquifers in a nearly flat environment. However, MRS can also be used in areas with sharp topographical variations. Following a review of the standard MRS equations when using a coincident transmitter/receiver loop, the mathematical terms potentially affected by tilting of the loop are discussed. We present the results of a numerical modelling exercise, studying a case where the surface is not horizontal and the loop cannot be considered to be parallel to the top of the aquifer. This shows that maximum variations in the MRS-signal amplitude are caused mainly by north- or south-dipping slopes. Slope effects depend on the loop size (a larger loop produces a larger error) especially in the presence of shallow water. With a geomagnetic-field inclination of 65° and a slope angle ≤ 10°, the topography causes a maximum variation in amplitude of less than 10%. Near magnetic poles and equator, the slope effect is lower and undetectable in most cases. It was found that within a 10% range of variation in the amplitude, errors introduced into inversions are within the typical uncertainty for MRS inversion and hence no topographic corrections are necessary. Thus, a significant effect from non-horizontal topography might be expected only when data uncertainty is lower than the slope effect (the slope effect is lower than equivalence when data quality is poor). Today, most field data sets are inverted using the modulus of the MRS signal, but some new developments consider the complex signal (both modulus and phase). However, inversion of complex MRS signals, which would provide a higher sensitivity to groundwater distribution, may be affected by slope effect. Thus, the slope orientation and dip angle should be accurately measured in the field when the phase of MRS signals is inverted too. 相似文献