共查询到20条相似文献，搜索用时 125 毫秒
1.
A Stanley Raj Y SRINIVAS D Hudson Oliver D Muthuraj 《Journal of Earth System Science》2014,123(2):395-411
The non-linear apparent resistivity problem in the subsurface study of the earth takes into account the model parameters in terms of resistivity and thickness of individual subsurface layers using the trained synthetic data by means of Artificial Neural Networks (ANN). Here we used a single layer feed-forward neural network with fast back propagation learning algorithm. So on proper training of back propagation networks it tends to give the resistivity and thickness of the subsurface layer model of the field resistivity data with reference to the synthetic data trained in the appropriate network. During training, the weights and biases of the network are iteratively adjusted to make network performance function level more efficient. On adequate training, errors are minimized and the best result is obtained using the artificial neural networks. The network is trained with more number of VES data and this trained network is demonstrated by the field data. The accuracy of inversion depends upon the number of data trained. In this novel and specially designed algorithm, the interpretation of the vertical electrical sounding has been done successfully with the more accurate layer model. 相似文献
2.
Effects of galvanic distortions on magnetotelluric data: Interpretation and its correction using deep electrical data 总被引：1，自引：0，他引：1
The non-inductive galvanic disturbances due to surficial bodies, lying smaller than high frequency skin depth, cause serious
interpretational errors in magnetotelluric data. These frequency independent distortions result in a quasi-static shift between
the apparent resistivity curves known as static shift. Two-dimensional modelling studies, for the effects of surficial bodies
on magnetotelluric interpretation, show that the transverse electric (TE) mode apparent resistivity curves are hardly affected
compared to the transverse magnetic (TM) mode curves, facilitating the correction by using a curve shifting method to match
low frequency asymptotes. But in the case of field data the problem is rather complicated because of the random distribution
of geometry and conductivity of near surface inhomogeneities. Here we present the use of deep resistivity sounding (DRS) data
to constrain MT static shift. Direct current sensitivity studies show that the behaviour of MT static shift can be estimated
using DC resistivity measurements close to the MT sounding station to appreciable depths. The distorted data set is corrected
using the MT response for DRS model and further subject to joint inversion with DRS data. Joint inversion leads to better
estimation of MT parameters compared to the separate inversion of data sets. 相似文献
3.
Much information on the regional lithospheric structure may come from MTS data acquired by synchronous 2D arrays and processed with regard to the nonlocal response of a laterally inhomogeneous subsurface. We suggest to invert the nonlocal MT responses applying correlation of all surface horizontal and vertical components of the geomagnetic field recorded simultaneously at all stations. The inversion algorithm has been applied to 2004–2005 European observatory data of diurnal Sq variations for first five harmonics and yielded lateral conductivity patterns for different periods. The maps show spatial correlation between conductivity maxima and lithospheric thickness minima and, specifically, highlight the contours of the Pannonian basin, where lithosphere is as thin as ～50 km, from seismic data. 相似文献
4.
The applications of intelligent techniques have increased exponentially in recent days to study most of the non-linear parameters.In particular,the behavior of earth resembles the non-linearity applications.An efficient tool is needed for the interpretation of geophysical parameters to study the subsurface of the earth.Artificial Neural Networks(ANN) perform certain tasks if the structure of the network is modified accordingly for the purpose it has been used.The three most robust networks were taken and comparatively analyzed for their performance to choose the appropriate network.The single-layer feed-forward neural network with the back propagation algorithm is chosen as one of the well-suited networks after comparing the results.Initially,certain synthetic data sets of all three-layer curves have been taken for training the network,and the network is validated by the Held datasets collected from Tuticorin Coastal Region(78°7′30″E and 8°48′45″N),Tamil Nadu.India.The interpretation has been done successfully using the corresponding learning algorithm in the present study.With proper training of back propagation networks,it tends to give the resistivity and thickness of the subsurface layer model of the field resistivity data concerning the synthetic data trained earlier in the appropriate network.The network is trained with more Vertical Electrical Sounding(VES) data,and this trained network is demonstrated by the field data.Groundwater table depth also has been modeled. 相似文献
5.
微机技术的发展和石油重磁勘探的需要促使了石油重磁资料微机处理和解释系统的建立。系统由输入、输出辅助子系统、重磁场分离和转换子系统及五个反演解释子系统八部分组成。五个反演解释子系统针对重磁力方法在油气勘探中所能解决的地质问题的类型不同而设计。计算模块的选择注意到方法的有效和适应，地球物理综合解释的需要及平面与剖面、定量正演与反演相结合的重要性。实际应用例子说明系统建立的成功和有效 相似文献
6.
Joint inversion of 2D resistivity and seismic travel time data to image saltwater intrusion over karstic areas 总被引：1，自引：0，他引：1
The simultaneous inversion of multiple geophysical data types has been proven to be a powerful tool to both improve subsurface imaging and help in the interpretation process. The main goal of this paper was to develop joint inversion strategies to provide improved resistivity and seismic velocity images for delineating saline water zones in karstic geological formations. The cross-gradient constraint approach was adopted to jointly invert resistivity and seismic first arrival data. The basic idea of this approach is to quantitatively estimate the structural similarity between resistivity and seismic velocity models, using the cross product of their gradients and to achieve a unified geological model which satisfies both data sets. Initially, synthetic data were employed to help develop a joint inversion strategy to be used over such complex geological structures. The proposed strategy uses a weighting factor for the cross-gradient constraints and separate damping factors for the resistivity and seismic data. This strategy was applied successfully on field data from the karstic region of Stilos, Crete, Greece. 相似文献
7.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
8.
The problem of multiphase phase flow in heterogeneous subsurface porous media is one involving many uncertainties. In particular, the permeability of the medium is an important aspect of the model that is inherently uncertain. Properly quantifying these uncertainties is essential in order to make reliable probabilistic-based predictions and future decisions. In this work, a measure-theoretic framework is employed to quantify uncertainties in a two-phase subsurface flow model in high-contrast media. Given uncertain saturation data from observation wells, the stochastic inverse problem is solved numerically in order to obtain a probability measure on the space of unknown permeability parameters characterizing the two-phase flow. As solving the stochastic inverse problem requires a number of forward model solves, we also incorporate the use of a conservative version of the generalized multiscale finite element method for added efficiency. The parameter-space probability measure is used in order to make predictions of saturation values where measurements are not available, and to validate the effectiveness of the proposed approach in the context of fine and coarse model solves. A number of numerical examples are offered to illustrate the measure-theoretic methodology for solving the stochastic inverse problem using both fine and coarse solution schemes. 相似文献
9.
This paper reviews the validity of earlier models obtained after quantitative interpretation of GDS data and presents a fresh
model using the inversion scheme EM2INV. The 2-D inversion of data is more objective than the earlier interpretation performed
by using trial and error method. The inversion results indicate that the present model differs from the earlier ones. The
reason could be that available GDS data are sufficient only for deriving the horizontal variation of subsurface resistivity.
In order to study the vertical resistivity variation additional MT sounding data would be required. It would therefore be
desirable to carry out MT survey in the specified area. A more comprehensive/appropriate model could be derived from joint
inversion of GDS and MT data. 相似文献
10.
The paper deals with modeling the response of pavement systems to subsurface deformations. An analytic formulation is developed in the context of layered elastic theory wherein the bottom boundary is distorted in a prescribed manner – forcing the medium to deform. In the basic case considered, an axisymmetric displacement field is imposed, consisting of vertical and radial components. A version of the model, dealing with a single elastic layer deformed only in the vertical direction (with frictionless bottom), is parametrically interrogated to study the influence of layer thickness and modulus on the resulting surface profile and vertical contact stresses. Furthermore, it is shown that by spatially superposing several axisymmetric solutions, general displacement shapes can be reproduced and analyzed. Such a scheme is demonstrated in an inverse type of analysis to a case study dealing with a sagging runway; using measured surface elevations, an attempt is made to backcalculate subsurface deformations and stress changes at the bottom interface of the structure. 相似文献
11.
Seawater intrusion (SWI) is a complex process, where 3D modeling is often necessary in order to monitor and manage the affected aquifers. Here, we present a synthetic study to test a joint hydrogeophysical inversion approach aimed at solving the inverse problem of estimating initial and current saltwater distribution. First, we use a 3D groundwater model for variable density flow based on discretized flow and solute mass balance equations. In addition to the groundwater model, a 3D geophysical model was developed for direct current resistivity imaging and inversion. The objective function of the coupled problem consists of data misfit and regularization terms as well as a coupling term that relates groundwater and geophysical states. We present a novel approach to solve the inverse problem using an alternating direction method of multipliers (ADMM) to minimize this coupled objective function. ADMM enables to treat the groundwater and geophysical part separately and thus use the existing software with minor changes. To further reduce the computational cost, the sensitivities are derived analytically for the discretized system of equations, which allows us to efficiently compute the gradients in the minimization procedure. The method was tested on different synthetic scenarios with groundwater and geophysical data represented by solute mass fraction data and direct current resistivity data. With the ADMM approach, we were able to obtain better estimates for the solute distribution compared to just considering each data separately, solving the problem with a simple coupled approach or by a direct substitution of the coupling constraint. 相似文献
12.
Faisal Rehman Sherif M. El-Hady Ali H. Atef Hussein M. Harbi 《Arabian Journal of Geosciences》2016,9(8):519
The present study deals with the seismic site classification of Bahrah area, Wadi Fatima, to characterize the local site conditions. The dynamic behavior of sediments was studied by the application of surface wave inversion. The multichannel analysis of surface waves (MASW) shallow geophysical technique was utilized for site classification. MASW survey was carried out at 66 sites along with 13 seismic refraction profiles at suitable localities. MASW and seismic refraction profiles were processed and compared with the available borehole data. The integration of MASW techniques with seismic refraction and borehole data progressively enhanced the subsurface visualization and reliability of the shear wave velocity estimation in the subsurface in the study area. The subsurface shear-wave velocity model was achieved by the solution of an inverse problem-based dispersion of surface waves and propagation in a vertically heterogeneous medium. The 2D genetic algorithm was employed for the inversion of dispersion curves to obtain velocity and thickness of subsurface layers. The depth to engineering bedrock and velocity of shear waves in the first 30 m was deciphered and mapped. The depth of bedrock in study area varies from 4 to 30 m, and V S 30 ranges from 320 to 800 m/s. The most of study area falls in B and C class categories in addition to few sites of D class according to the NEHRP guidelines. 相似文献
13.
Geostatistical Seismic Inversion with Direct Sequential Simulation and Co-simulation with Multi-local Distribution Functions 总被引：1，自引：0，他引：1
Rúben Nunes Amílcar Soares Leonardo Azevedo Pedro Pereira 《Mathematical Geosciences》2017,49(5):583-601
Stochastic sequential simulation is a common modelling technique used in Earth sciences and an integral part of iterative geostatistical seismic inversion methodologies. Traditional stochastic sequential simulation techniques based on bi-point statistics assume, for the entire study area, stationarity of the spatial continuity pattern and a single probability distribution function, as revealed by a single variogram model and inferred from the available experimental data, respectively. In this paper, the traditional direct sequential simulation algorithm is extended to handle non-stationary natural phenomena. The proposed stochastic sequential simulation algorithm can take into consideration multiple regionalized spatial continuity patterns and probability distribution functions, depending on the spatial location of the grid node to be simulated. This work shows the application and discusses the benefits of the proposed stochastic sequential simulation as part of an iterative geostatistical seismic inversion methodology in two distinct geological environments in which non-stationarity behaviour can be assessed by the simultaneous interpretation of the available well-log and seismic reflection data. The results show that the elastic models generated by the proposed stochastic sequential simulation are able to reproduce simultaneously the regional and global variogram models and target distribution functions relative to the average volume of each sub-region. When used as part of a geostatistical seismic inversion procedure, the retrieved inverse models are more geologically realistic, since they incorporate the knowledge of the subsurface geology as provided, for example, by seismic and well-log data interpretation. 相似文献
14.
《Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy》2000,25(2):195-200
In the present paper we study the possibility of using joint quantitative interpretation of experimental electromagnetic and seismic data for lithospheric studies and to establish the situations in which the joint inversion of seismic and EM experimental data is possible. For this purpose a theoretical model of a solid-liquid rock mixture with fractal structure of grain-to-grain contact (Bahr, 1997; Spangenberg, 1998) was used to investigate dependence of elastic and electrical properties on rock microstructure. The electrical conductivity and seismic wave velocity were estimated by Monte-Carlo modeling for low porosity water-saturated crystalline rock and olivine-partially molten basalt system. It was shown that both elastic properties and conductivity of solid-liquid rock mixture depend non-linearly upon two main common factors, i.e. amount of liquid phase and fractal rock microstructure that is controlled by the degree of contact between grains. The dependence of rock physical macroproperties (seismic velocities, electrical conductivity and seismic and electrical anisotropy) upon these two main factors makes it possible to construct algorithms for joint quantitative interpretation of experimental seismic and electromagnetic data. 相似文献
15.
地震波速层析成像方法研究进展 总被引：3，自引：0，他引：3
本文回顾了利用天然地震观测获取地下速度结构的方法。尽管有不同的新方法涌现,天然地震波速层析成像方法,尤其是多震相联合反演的格点层析成像方法,是当今使用广泛使用的层析成像方法之一,是对地球内部成像的最有效方式。波速层析成像方法的未来发展首先是提高第一手的观测资料,即增加接收地震波信息的地震台站分布密度;同时,通过多种地球物理方法联合反演相互约束可以给出较为严格的地球物理模型,并来降低了地球物理反演和解译的多解性,这是地球物理探测研究的趋势,也是天然地震波速层析成像方法的研究趋势。 相似文献
16.
This paper presents the combination of two complementary methods, magnetic and joint inversion of resistivity/TEM data, as
an effective approach to characterize groundwater reservoirs. Twenty stations of transient electromagnetic (TEM) and vertical
electrical soundings (VES) were measured and interpreted using a joint inversion technique to evaluate the subsurface low
resistivity zones connected to the groundwater reservoir. A complementary survey including 871 land magnetic stations was
carried out at the same area to detect the upper surface of the basaltic sheet, which represents the bottom of the Miocene
aquifer in the study area. The geological log from one borehole drilled in the zone was used to partially calibrate the calculated
models. The results revealed that the study area consists of five different geological units with ages ranging from Paleogene
(Oligocene) to Quaternary. The methodology provides good results at a very low cost when compared with drilling boreholes. 相似文献
17.
This study concerns the identification of parameters of soil constitutive models from geotechnical measurements by inverse analysis. To deal with the non‐uniqueness of the solution, the inverse analysis is based on a genetic algorithm (GA) optimization process. For a given uncertainty on the measurements, the GA identifies a set of solutions. A statistical method based on a principal component analysis (PCA) is, then, proposed to evaluate the representativeness of this set. It is shown that this representativeness is controlled by the GA population size for which an optimal value can be defined. The PCA also gives a first‐order approximation of the solution set of the inverse problem as an ellipsoid. These developments are first made on a synthetic excavation problem and on a pressuremeter test. Some experimental applications are, then, studied in a companion paper, to show the reliability of the method. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
18.
本文在大地电磁测深松弛（ｒａｐｉｄ）反演理论方法的基础上，探讨了电偶源频率电磁测深（赤道装置）阻抗视电阻率的二维反演问题。由二维层状模型反演目标函数的近似简化出发。将全剖面的二维反演约化为沿剖面各测点的逐一测点一维反演，从而构制了一个简易可行的二维层状模型参数化反演方案。反演试验表明，这一反演方法是实用和有效的 相似文献
19.
We present inversion results for a 100 site, broadband magnetotelluric (MT) survey in the Penola Trough, Otway Basin, South Australia. The Penola Trough is host to several petroleum reservoirs and has more recently been a target for unconventional geothermal exploration. We present two interpretations of the MT data. A 1D anisotropic interpretation, where anisotropy is determined within the Otway Basin sequence and basement in the northeastern Penola Trough, fits the impedance tensor well. However, the anisotropy strike is inconsistent with the known orientation of electrically conductive fractures in the Penola Trough. On the other hand, a 3D interpretation, which incorporates lateral variations in resistivity, requires no anisotropy yet it matches the data equally well. Both the 1D and 3D inversions resolve several layers within the Otway Basin sequence, which correspond to stratigraphic units defined in wells and in the coincident Haselgrove–Balnaves 3D seismic survey. These include the Eumeralla and Dilwyn formations, which are poorly resolved in the seismic data. The basin architecture, defined in the 3D inversion, in particular the depth to basement, is consistent with previous interpretations based on seismic reflection data that show that the Otway Basin thins in the northeastern Penola Trough. This does not occur in the anisotropic model. We therefore conclude that the subsurface resistivity appears to be isotropic in the Penola Trough. This contrasts with the anisotropic resistivity structure determined in a previous study in the Koroit region, eastern Otway Basin. The difference in the MT responses between the two regions is supported by resistivity and permeability information from well logs and may reflect differences in the orientation of subsurface fractures, or differences in the present-day stress field, between the two regions. 相似文献
20.
A Distance-based Prior Model Parameterization for Constraining Solutions of Spatial Inverse Problems 总被引：4，自引：3，他引：1
Spatial inverse problems in the Earth Sciences are often ill-posed, requiring the specification of a prior model to constrain
the nature of the inverse solutions. Otherwise, inverted model realizations lack geological realism. In spatial modeling,
such prior model determines the spatial variability of the inverse solution, for example as constrained by a variogram, a
Boolean model, or a training image-based model. In many cases, particularly in subsurface modeling, one lacks the amount of
data to fully determine the nature of the spatial variability. For example, many different training images could be proposed
for a given study area. Such alternative training images or scenarios relate to the different possible geological concepts
each exhibiting a distinctive geological architecture. Many inverse methods rely on priors that represent a single subjectively
chosen geological concept (a single variogram within a multi-Gaussian model or a single training image). This paper proposes
a novel and practical parameterization of the prior model allowing several discrete choices of geological architectures within
the prior. This method does not attempt to parameterize the possibly complex architectures by a set of model parameters. Instead,
a large set of prior model realizations is provided in advance, by means of Monte Carlo simulation, where the training image
is randomized. The parameterization is achieved by defining a metric space which accommodates this large set of model realizations.
This metric space is equipped with a “similarity distance” function or a distance function that measures the similarity of
geometry between any two model realizations relevant to the problem at hand. Through examples, inverse solutions can be efficiently
found in this metric space using a simple stochastic search method. 相似文献