首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 218 毫秒
1.
We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.  相似文献   

2.
We have developed a least-squares method to determine simultaneously the depth and the width of a buried thick dipping dike from residualized magnetic data using filters of successive window lengths. The method involves using a relationship between the depth and the half-width of the source and a combination of windowed observations. The relationship represents a family of curves (window curves). For a fixed window length, the depth is determined for each half-width value by solving one nonlinear equation of the form f (z) = 0 using the least-squares method. The computed depths are plotted against the width values representing a continuous curve. The solution for the depth and the width of the buried dike is read at the common intersection of the window curves. The method involves using a dike model convolved with the same moving average filter as applied to the observed data. As a result, this method can be applied to residuals as well as to measured magnetic data. Procedures are also formulated to estimate the amplitude coefficient and the index parameter. The method is applied to theoretical data with and without random errors. The validity of the method is tested on airborne magnetic data from Canada and on a vertical component magnetic anomaly from Turkey. In all cases examined, the model parameters obtained are in good agreement with the actual ones and with those given in the published literature.  相似文献   

3.
Ground water flow associated with pumping and injection tests generates self-potential signals that can be measured at the ground surface and used to estimate the pattern of ground water flow at depth. We propose an inversion of the self-potential signals that accounts for the heterogeneous nature of the aquifer and a relationship between the electrical resistivity and the streaming current coupling coefficient. We recast the inversion of the self-potential data into a Bayesian framework. Synthetic tests are performed showing the advantage in using self-potential signals in addition to in situ measurements of the potentiometric levels to reconstruct the shape of the water table. This methodology is applied to a new data set from a series of coordinated hydraulic tomography, self-potential, and electrical resistivity tomography experiments performed at the Boise Hydrogeophysical Research Site, Idaho. In particular, we examine one of the dipole hydraulic tests and its reciprocal to show the sensitivity of the self-potential signals to variations of the potentiometric levels under steady-state conditions. However, because of the high pumping rate, the response was also influenced by the Reynolds number , especially near the pumping well for a given test. Ground water flow in the inertial laminar flow regime is responsible for nonlinearity that is not yet accounted for in self-potential tomography. Numerical modeling addresses the sensitivity of the self-potential response to this problem.  相似文献   

4.
In the present paper a new method is proposed for the quantitative interpretation of self-potential anomalies which are produced by a vertical dipole. First the mathematical expression of the wavenumber spectrum of the self-potential anomaly is deduced. It is pointed out that at relatively high wavenumbers the behavior of the amplitude spectrum is controlled by the closer to the surface pole at depth h. On the other hand, the “width” of the amplitude spectrum depends on the depth h and the dipole length L.Making a proper mathematical transformation of the amplitude spectrum, and applying the least squares method, it is possible to calculate the depth to the upper pole. The dipole length may then be calculated, by solving numerically a characteristic algebraic equation, as long as the “width” of the amplitude spectrum has been previously defined.The proposed method is applied on a well known self-potential profile from Greece. The calculated parameters of the polarized body are in good agreement with real data. Experimentation with synthetic models in which random noise was introduced, showed that this method gives reliable results if the noise amplitude is not more than 20% of the signal amplitude. It is clearly more efficient than the methods which are based on the model of the point pole or the dipole with a small length. It can also give good results if the horizontal extensions of the polarized body are not more than a few tenths of the depth of the upper pole. If the polarized body is tilted, the depth of the upper pole can be calculated with satisfactory accuracy.The direct interpretation method which is proposed in the present paper, may be useful in mineral exploration, and particularly if the target of interest is the detection of massive sulfide mineralization.  相似文献   

5.
A new approach is proposed in order to interpret field self-potential (SP) anomalies related to simple geometric-shaped models such as sphere, horizontal cylinder, and vertical cylinder. This approach is mainly based on solving a set of algebraic linear equations, and directed towards the best estimate of the three model parameters, e.g., electric dipole moment, depth, and polarization angle. Its utility and validity are demonstrated through studying and analyzing synthetic self-potential anomalies obtained by using simulated data generated from a known model and a statistical distribution with different random errors components. Being theoretically tested and proven, this approach has been consequently applied on two real field self-potential anomalies taken from Colorado and Turkey. A comparable and acceptable agreement is obtained between the results derived by the new proposed method and those deduced by other interpretation methods. Moreover, the depth obtained by such an approach is found to be very close to that obtained by drilling information.  相似文献   

6.
The self-potential anomaly due to a two-dimensional inclined sheets of finite depth extent has been analysed in the frequency domain using the Fourier transform. Expression for the Fourier amplitude and phase spectra are derived. The Fourier amplitude and phase spectra are analysed so as to evaluate the parameters of the sheet. Application of this method on two anomalies (synthetic and field data) has given good results.  相似文献   

7.
自然电场法常用于环境与工程等领域的监测作业,但各时刻观测数据往往单独反演解释.为了充分利用时序数据间的关联信息,提高监测数据的反演解释可靠性,提出基于卡尔曼滤波的自然电场监测数据时序反演方法.根据达西定律和阿尔奇公式建立污染物在孔隙介质中的运动扩散的动态地电模型,作为用于构建卡尔曼滤波的状态模型.而卡尔曼滤波的观测模型则通过常规的自然电场法正演获得.在建立状态模型和观测模型的基础上,构建起卡尔曼滤波递归,将地电模型演化信息与自然电场观测数据进行信息融合,实现自然电场监测数据的时序反演.加入噪声的自然电场模拟数据测试表明时序反演算法具有较好的鲁棒性,对噪声不敏感.沙槽物理实验监测数据的计算测试也同样证明时序反演能有效处理监测数据,实现对动态模型的准确重构.  相似文献   

8.
A quantitative interpretation method of self-potential field anomalies has been proposed. The method is designed and implemented for the estimation of center depth, electric dipole moment or magnitude of polarization, polarization angle, and geometric shape factor of a buried body from SP field data, related to simple geometric structures such as cylinders, spheres and sheet-like bodies. The proposed method is based on Fair function minimization and also on stochastic optimization modeling. This new technique was first tested on theoretical synthetic data randomly generated by a chosen statistical distribution from a known model with different random noise components. Such mathematical simulation shows a very close agreement between assumed and estimated model parameters. Being theoretically proven, it has been applied and tested on self-potential field data taken from the United States, Germany, India and Turkey. The agreement between results obtained by the suggested method and those obtained by other previous methods is good and comparable. Moreover, the depth obtained by this method is found to be in high accordance with that obtained from drilling information.  相似文献   

9.
Summary A quantitative method of interpreting self-potential anomaly caused by a spherical ore body using downward continuation method is presented. Master curves to determine the depth, radius and angle of polarization have been prepared.  相似文献   

10.
Iodine geyser, located in the Waimangu Geothermal Valley (New Zealand), has been studied by both passive electrical and seismic methods. The activity of the geyser was monitored at various distances from the vent using self-potential method. The self-potential signals display cyclic negative variations with respect to a baseline drawn when the geyser is quiet. The minimum in the self-potential signals coincides with the maximum overflow. We provide a numerical model able to explain both the polarity and magnitude of the observed signal. This model is based on the fluctuations of the hydraulic head in the conduit of the geyser; the divergence of the streaming current density is created at the interface between the pipe and the surrounding rock. Passive seismic experiments were used to localize ambient noise sources. These signals have been processed with the so-called Matched-Field Processing technique (MFP); a dominant source emerged from this processing, that we characterized in range and depth with a good accuracy.  相似文献   

11.
自然电位、自然伽马测井曲线形态和特征与沉积相带及其储集砂体关系密切,它们对不同岩性地层特征响应十分敏感.利用自然伽马、自然电位同步减小的较大幅度评价识别渗砂层,指示划分水下分流河道主体沉积微相带;利用自然伽马、自然电位减小的幅度差评价识别低渗砂、致密砂层,指示划分水下分流河道堤泛(侧翼)沉积微相带;并以自然电位比自然伽马减小幅度的相对细小差异识别油水层.在安塞油田沿河湾地区长6自然伽马、自然电位曲线划分沉积相带及储层的应用中,建立了自然伽马、自然电位划分储层的下限标准,有效地评价了特低渗透储层沉积微相带及储层类型,提高和扩大了测井曲线的应用及效果.  相似文献   

12.
A nomogram has been devised for situations, in which the source of a self-potential anomaly can be approximated by an obliquely polarized sphere or horizontal cylinder embedded in a homogeneous half space. The nomogram can be used for rapid determination of three parameters of the target: (1) depth to the centre, (2) angle between the axis of polarization and the horizontal, (3) shift of the point vertically above the centre of the body from zero potential value. The nomogram has been tested and the parameters determined for SP results obtained over ore bodies Weiss and Süleymanköy in the Ergani Copper district, Turkey. The curves computed for the estimated parameters match the field curves well.  相似文献   

13.
The assessment of hydraulic conductivity of heterogeneous aquifers is a difficult task using traditional hydrogeological methods (e.g., steady state or transient pumping tests) due to their low spatial resolution. Geophysical measurements performed at the ground surface and in boreholes provide additional information for increasing the resolution and accuracy of the inverted hydraulic conductivity field. We used a stochastic joint inversion of Direct Current (DC) resistivity and self-potential (SP) data plus in situ measurement of the salinity in a downstream well during a synthetic salt tracer experiment to reconstruct the hydraulic conductivity field between two wells. The pilot point parameterization was used to avoid over-parameterization of the inverse problem. Bounds on the model parameters were used to promote a consistent Markov chain Monte Carlo sampling of the model parameters. To evaluate the effectiveness of the joint inversion process, we compared eight cases in which the geophysical data are coupled or not to the in situ sampling of the salinity to map the hydraulic conductivity. We first tested the effectiveness of the inversion of each type of data alone (concentration sampling, self-potential, and DC resistivity), and then we combined the data two by two. We finally combined all the data together to show the value of each type of geophysical data in the joint inversion process because of their different sensitivity map. We also investigated a case in which the data were contaminated with noise and the variogram unknown and inverted stochastically. The results of the inversion revealed that incorporating the self-potential data improves the estimate of hydraulic conductivity field especially when the self-potential data were combined to the salt concentration measurement in the second well or to the time-lapse cross-well electrical resistivity data. Various tests were also performed to quantify the uncertainty in the inverted hydraulic conductivity field.  相似文献   

14.
Application of particle swarm optimization on self-potential data   总被引：1，自引：0，他引：1  
Particle swarm optimization (PSO) is a global search method, which can be used for quantitative interpretation of self-potential data in geophysics. At the result of this process, parameters of a source model, e.g., the electrical dipole moment, the depth of the source, the distance from the origin, the polarization angle and regional coefficients are estimated. This study investigates the results and interpretation of a detailed numerical data of some simple body responses, contaminated and field data. The method is applied to three field examples from Turkey and the results are compared with the previous works. The statistics of particle swarm optimization and the corresponding model parameters are analyzed with respect to the number of generation. We also present the oscillations of the model parameters at the vicinity of the low misfit area. Further, we show how the model parameters and absolute frequencies are related to the total number of PSO iterations. Gaussian noise shifts the low misfit area region from the correct parameter values proportional to the level of errors, which directly affects the result of the PSO method. These effects also give some ambiguity of the model parameters. However, the statistical analyses help to decrease these ambiguities in order to find the correct values. Thus, the findings suggest that PSO can be used for quantitative interpretation of self-potential data.  相似文献   

15.
We have developed a new numerical method to determine the shape (shape factor), depth, polarization angle, and electric dipole moment of a buried structure from residual self-potential (SP) anomalies. The method is based on defining the anomaly value at the origin and four characteristic points and their corresponding distances on the anomaly profile. The problem of shape determination from residual SP anomaly has been transformed into the problem of finding a solution to a nonlinear equation of the form q = f (q). Knowing the shape, the depth, polarization angle and the electric dipole moment are determined individually using three linear equations. Formulas have been derived for spheres and cylinders. By using all possible combinations of the four characteristic points and their corresponding distances, a procedure is developed for automated determination of the best-fit-model parameters of the buried structure from SP anomalies. The method was applied to synthetic data with 5% random errors and tested on a field example from Colorado. In both cases, the model parameters obtained by the present method, particularly the shape and depth of the buried structures are found in good agreement with the actual ones. The present method has the capability of avoiding highly noisy data points and enforcing the incorporation of points of the least random errors to enhance the interpretation results.  相似文献   

16.
—We have developed a least-squares minimization approach to determine the shape (shape-factor) of a buried polarized body from a residual self-potential anomaly profile. By defining the zero anomaly distance and the anomaly value at the origin on the profile, the problem of the shape-factor determination is transformed into the problem of finding a solution of a nonlinear equation of the form f(q) = 0. Procedures are also formulated to estimate the depth of polarization angle, and the electric dipole moment. The method is applied to synthetic data with and without random noise. The obtained shape-factor agrees very well with the model shape-factor when using synthetic data. After adding ± 2 percent random error in the synthetic data, the shape factor obtained is within ± 4 percent. Finally the validity of the method is tested on a field example from the Ergani copper district, Turkey.  相似文献   

17.
The processing of multicomponent seismic data, carried out individually on the different wavetypes (P-, S- and converted waves), should result in an improved image of the subsurface. We examine the wavefield‐separation method proposed by Cho and Spencer . We discuss practical aspects related to the separation of interfering waves on two-component surface seismic data and illustrate these using synthetic data. A sliding spatial window is used for analysis. The choice of its width represents a trade-off between stabilizing the method in the presence of random noise and ensuring a good spatial resolution. No a priori knowledge of the subsurface is required, but locally, the characteristic parameters of the waves, i.e. horizontal slowness and polarization, are assumed to be constant inside the analysis window. These parameters are estimated at each frequency, but a statistical analysis provides a more robust estimate, especially in the presence of random noise. This approach also solves the problem of eigenvalue sharing and switching. Additional smoothing of the estimates according to a model may further improve the results. The width of the analysis window may be decreased if the waves inside the data window differ significantly in amplitude. The dominant wave in each case is separated from the lower-amplitude waves and subtracted from the data. This novel iterative approach thereby allows for the isolation of low-amplitude events.  相似文献   

18.
Salt and related structures have played important roles in controlling hydrocarbon accumulations in the Gulf of Mexico. Using a two-dimensional fluid flow/compaction model, which allows for both conduction and convection of heat, an examination is given of the effects on thermal patterns of the combined influence of multiple salt features, including diapirs, pillows, sheets and wedges. The focusing and defocusing of heat due to the higher thermal conductivity of salt are accounted for in the modeling. The results show that there could be as much as a 30°C anomaly above multi-salt bodies due to the focusing of heat by salt, and as much as 50°C temperature contrast between internal salt positions and sediments external to the salt in the deep part of a section. The magnitude of the thermal anomaly depends on the size (or width) of the salt and on the depth of the rooted salt. The modeled results provide estimates of the influence of salt in expanding the oil generation window by approximately half of the salt thickness.  相似文献   

19.
Common studies on the static electric field distribution over a conductivity anomaly use the self-potential method. However, this method is time consuming and requires nonpolarizable electrodes to be placed in the ground. Moreover, the information gained by this method is restricted to the horizontal variations of the electric field. To overcome the limitation in the self-potential technique, we conducted a field experiment using a non conventional technique to assess the static electric field over a conductivity anomaly. We use two metallic potential probes arranged on an insulated boom with a separation of 126 cm. When placed into the electric field of the free air, a surface charge will be induced on each probe trying to equalize with the potential of the surrounding atmosphere. The use of a plasma source at both probes facilitated continuous and quicker measurement of the electric field in the air. The present study shows first experimental measurements with a modified potential probe technique (MPP) along a 600-meter-long transect to demonstrate the general feasibility of this method for studying the static electric field distribution over shallow conductivity anomalies.Field measurements were carried out on a test site on top of the Bramsche Massif near Osnabrück (Northwest Germany) to benefit from a variety of available near surface data over an almost vertical conductivity anomaly. High resolution self-potential data served in a numerical analysis to estimate the expected individual components of the electric field vector.During the experiment we found more anomalies in the vertical and horizontal components of the electric field than self-potential anomalies. These contrasting findings are successfully cross-validated with conventional near surface geophysical methods. Among these methods, we used self-potential, radiomagnetotelluric, electric resistivity tomography and induced polarization data to derive 2D conductivity models of the subsurface in order to infer the geometrical properties and the origin of the conductivity anomaly in the survey area. The presented study demonstrates the feasibility of electric field measurements in free air to detect and study near surface conductivity anomalies. Variations in Ez correlate well with the conductivity distribution obtained from resistivity methods. Compared to the self-potential technique, continuously free air measurements of the electric field are more rapid and of better lateral resolution combined with the unique ability to analyze vertical components of the electric field which are of particular importance to detect lateral conductivity contrasts. Mapping Ez in free air is a good tool to precisely map lateral changes of the electric field distribution in areas where SP generation fails. MPP offers interesting application in other geophysical techniques e.g. in time domain electromagnetics, DC and IP.With this method we were able to reveal a ca. 150 m broad zone of enhanced electric field strength.  相似文献   

20.
v--vUsing Frank and Wolfe's algorithm, a new interesting nonlinear programming technique has been developed in an attempt to estimate the geometric shape factor of a buried polarized body from a residual self-potential anomaly. Furthermore, the depth, the polarization angle and the electrical dipole moment have also been derived. This algorithm is noted to be robust and its application to SP data converges rapidly towards the optimal solution. The developed technique is tested through studying synthetic data with and without random noise. As a result, the near agreement between the model geometric shape factor and the evaluated one is well recognized. The validity of this proposed technique is tested on a field example from the Ergani Copper district, Turkey. The superiority of the nonlinear programming technique over other recently published methods is shown.  相似文献   

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

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