共查询到20条相似文献，搜索用时 15 毫秒
1.
Gary Barnes 《Geophysical Prospecting》2014,62(3):646-657
When anomalous gravity gradient signals provide a large signal‐to‐noise ratio, airborne and marine surveys can be considered with wide line spacing. In these cases, spatial resolution and sampling requirements become the limiting factors for specifying the line spacing, rather than anomaly detectability. This situation is analysed by generating known signals from a geological model and then sub‐sampling them using a simulated airborne gravity gradient survey with a line spacing much wider than the characteristic anomaly size. The data are processed using an equivalent source inversion, which is used subsequently to predict and grid the field in‐between the survey lines by means of forward calculations. Spatial and spectral error analysis is used to quantify the accuracy and resolution of the processed data and the advantages of acquiring multiple gravity gradient components are demonstrated. With measurements of the full tensor along survey lines spaced at 4 × 4 km, it is shown that the vertical gravity gradient can be reconstructed accurately over a bandwidth of 2 km with spatial root‐mean square errors less than 30%. A real airborne full‐tensor gravity gradient survey is presented to confirm the synthetic analysis in a practical situation. 相似文献
2.
Site characterization activities at potential unexploded ordnance (UXO) sites rely on sparse sampling collected as geophysical
surveys along strip transects. From these samples, the locations of target areas, those regions on the site where the geophysical
anomaly density is significantly above the background density, must be identified. A target area detection approach using
a hidden Markov model (HMM) is developed here. HMM’s use stationary transition probabilities from one state to another for
steps between adjacent locations as well as the probability of any particular observation occurring given each possible underlying
state. The approach developed here identifies the transition probabilities directly from the conceptual site model (CSM) created
as part of the UXO site characterization process. A series of simulations examine the ability of the HMM approach to simultaneously
determine the target area locations within each transect and to estimate the unknown anomaly intensity within the identified
target area. The HMM results are compared to those obtained using a simpler target detection approach that considers the background
anomaly density to be defined by a Poisson distribution and each location to be independent of any adjacent location. Results
show that the HMM approach is capable of accurately identifying the target locations with limited false positive identifications
when both the background and target are intensities are known. The HMM approach is relatively robust to changes in the initial
estimate of the target anomaly intensity and is capable of identifying target locations and the corresponding target anomaly
intensity when this intensity is approximately 60% higher than the background intensity at intensities that are representative
of actual field sites. Application to data collected from a wide area assessment field site show that the HMM approach identifies
the area of the site with elevated anomaly intensity with few false positives. This field site application also shows that
the HMM results are relatively robust to changes in the transect width. 相似文献
3.
依据航空重力测量基本原理,构建了航空重力异常解算的卡尔曼滤波模型,将新息自适应卡尔曼滤波器(IAE,Innovation based Adaptive Estimation)应用于量测噪声未知的航空重力异常解算.针对IAE滤波器滑动窗口宽度难以准确确定的问题,通过对多个不同滑动窗口新息协方差估计的加权平均,获得改进的IAE滤波器,该IAE滤波器不仅具有量测噪声自适应估计能力,还能实现滑动采样窗口的优化选取.试验结果表明,IAE滤波器可以降低因量测噪声统计信息不明引起的解算误差,改进IAE解算的重力异常误差约为1mGal. 相似文献
4.
介绍了“中国大陆构造环境监测网络”（简称“陆态网络”）重力观测体系。从分形几何学角度研究了陆态网络重力测网的分形特征,计算得到测网分维数Df=1.5598、最佳网格化间距r=139km。综合测网分维数与重力时变距、重力异常范围、测网范围等量化指标,对陆态网络重力测网的地震监测能力进行分析,结果表明,测网最佳网格化间距等特征值与MS6.0地震引发的重力场变化区域的特征异常区半径等统计量化指标接近,具备监测我国除边界及藏北无人区以外大陆区域的MS6.0及以上地震的能力。 相似文献
5.
Pierre Keating 《Geophysical Prospecting》1995,43(4):569-580
Gravity data are often acquired over long periods of time using different instruments and various survey techniques, resulting in data sets of non-uniform accuracy. As station locations are inhomogeneously distributed, gravity values are interpolated on to a regular grid to allow further processing, such as computing horizontal or vertical gradients. Some interpolation techniques can estimate the interpolation error. Although estimation of the error due to interpolation is of importance, it is more useful to estimate the maximum gravity anomaly that may have gone undetected by a survey. This is equivalent to the determination of the maximum mass whose gravity anomaly will be undetected at any station location, given the data accuracy at each station. Assuming that the maximum density contrast present in the survey area is known or can be reasonably assumed from a knowledge of the geology, the proposed procedure is as follows: at every grid node, the maximum mass whose gravity anomaly does not disturb any of the surrounding observed gravity values by more than their accuracies is determined. A finite vertical cylinder is used as the mass model in the computations. The resulting map gives the maximum detection error and, as such, it is a worst-case scenario. Moreover, the map can be used to optimize future gravity surveys: new stations should be located at, or near, map maxima. The technique is applied to a set of gravity observations obtained from different surveys made over a period of more than 40 years in the Abitibi Greenstone Belt in eastern Canada. 相似文献
6.
Summary Conditions are given under which two thick plates, differing in dip, apparent susceptibility, and remanence, will produce similar magnetic anomalies. From these conditions correction formulae are developed. Using these formulae the dip and susceptibility of a plate with remanent magnetization can be obtained from those of non-remanent plate. An interpretation procedure is suggested where the magnetic anomaly is first interpreted by means of a plate without remanence, dip and apparent susceptibility are then estimated by using the correction formulae developed. Thickness, position and depth of the plate are unaffected by the remanence correction procedure. The procedure is independent of the field component measured. 相似文献
7.
Following the probability tomography principles previously introduced to image the sources of electric and electromagnetic anomalies, we demonstrate that a similar approach can be used to analyse gravity data. First, we give a coherent derivation of the Bouguer anomaly concept as a Newtonian-type integral for an arbitrary mass distribution buried below a non-flat topography. A discretized solution of this integral is then derived as a sum of elementary contributions, which are cross-correlated with the gravity data function in the expression for the total power associated with the Bouguer anomaly. To image the mass distribution underground we introduce a mass contrast occurrence probability function using the cross-correlation product of the observed Bouguer anomaly and the synthetic field due to an elementary mass contrast source. The tomographic procedure consists of scanning the subsurface with the elementary source and calculating the occurrence probability function at the nodes of a regular grid. The complete set of grid values is used to highlight the zones of highest probability of mass contrast concentrations. Some synthetic and field examples demonstrate the reliability and resolution of the new gravity tomographic approach. 相似文献
8.
Bayes estimate of the probability of exceedance of annual floods 总被引：1，自引：1，他引：1
In this paper Lindley's Bayesian approximation procedure is used to obtain the Bayes estimate of the probability of exceedence of a flood discharge. The Bayes estimates of the probability of exceedence has been shown by S.K. Sinha to be equivalent to the estimate of the probability of exceedence from the predictive or Bayesian disribution, of a future flood discharge. The evaluation of complex ratios of multiple integrals common in a Bayesian analysis is not necessary using Lindley's procedure. The Bayes estimates are compared to those obtained by the method of maximum likelihood and the method of moments. The results show that Bayes estimates of the probability of exceedence are larger as expected, but have smaller posterior standard deviations. 相似文献
9.
Sean E. Walker Leonard R. Pasion Douglas W. Oldenburg Stephen D. Billings 《Journal of Applied Geophysics》2007,61(3-4):254-278
Using field data and numerical simulations we investigate the effect of data quality on time domain electromagnetic discrimination. Data quality decreases when measurements contain responses not accounted for by our mathematical modelling. This can include instrument noise, inaccurately reported position and orientation information, geologic contributions to the signal, and loss of validity of the forward modelling. Survey design is critical to data quality in order to have sufficient sampling of data anomalies, and also to ensure that each target is illuminated such that both the axial and transverse components of the polarization can be excited and measured. For dipole model based discrimination algorithms, success is contingent upon the accuracy with which the components of the polarization tensor can be estimated. Field data from different survey modes are analysed to identify noise sources and provide quantitative estimates of the noise in each survey. Inversion results show that increased noise levels lead to greater spread in recovered parameters. Monte Carlo simulations are performed in order to investigate the importance of other data quality factors. Analysis of inversion results from the simulations show that anomaly size, signal to noise ratio, positioning error, line spacing and station spacing all play a role in the spread of recovered parameters. Through the analysis of our simulation results we propose a figure of merit as a means of quantifying different data quality factors with a single number and relate this number to the accuracy with which parameters can be estimated. 相似文献
10.
This paper presents a new inversion method for the interpretation of 2D magnetic anomaly data, which uses the combination of the analytic signal and its total gradient to estimate the depth and the nature (structural index) of an isolated magnetic source. However, our proposed method is sensitive to noise. In order to lower the effect of noise, we apply upward continuation technique to smooth the anomaly. Tests on synthetic noise-free and noise corrupted magnetic data show that the new method can successfully estimate the depth and the nature of the causative source. The practical application of the technique is applied to measured magnetic anomaly data from Jurh area, northeast China, and the inversion results are in agreement with the inversion results from Euler deconvolution of the analytic signal. 相似文献
11.
We present an inversion technique based on the Marquardt algorithm to estimate the depth of a 2.5-D sedimentary basin in addition
to the regional gravity anomaly that is associated with the residual gravity anomaly, wherein the density contrast varies
parabolically with depth. Forward modeling is carried out through a derived analytical gravity expression of a 2.5-D vertical
prism. Inversion of a theoretical gravity anomaly with and without a regional gravity anomaly illustrates the procedure that
it is found to be insensitive to the regional gravity effect. Furthermore, the algorithm is exemplified with the gravity anomalies
of the derived density-depth model of the Godavari subbasin, India with a parabolic density profile resulting in a more consistent
geological model rather than a constant density profile. The main advantage of this method is that it works well even when
the profile of interpretation does not bisect the strike length of the sedimentary basin. 相似文献
12.
E. M. Abdelrahman T. M. El-Araby A. A. Ammar H. I. Hassanein 《Pure and Applied Geophysics》1997,150(1):121-128
—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. 相似文献
13.
Field testing of modular borehole monitoring with simultaneous distributed acoustic sensing and geophone vertical seismic profiles at Citronelle,Alabama 
下载免费PDF全文
下载免费PDF全文 A modular borehole monitoring concept has been implemented to provide a suite of well‐based monitoring tools that can be deployed cost effectively in a flexible and robust package. The initial modular borehole monitoring system was deployed as part of a CO2 injection test operated by the Southeast Regional Carbon Sequestration Partnership near Citronelle, Alabama. The Citronelle modular monitoring system transmits electrical power and signals, fibre‐optic light pulses, and fluids between the surface and a reservoir. Additionally, a separate multi‐conductor tubing‐encapsulated line was used for borehole geophones, including a specialized clamp for casing clamping with tubing deployment. The deployment of geophones and fibre‐optic cables allowed comparison testing of distributed acoustic sensing. We designed a large source effort (>64 sweeps per source point) to test fibre‐optic vertical seismic profile and acquired data in 2013. The native measurement in the specific distributed acoustic sensing unit used (an iDAS from Silixa Ltd) is described as a localized strain rate. Following a processing flow of adaptive noise reduction and rebalancing the signal to dimensionless strain, improvement from repeated stacking of the source was observed. Conversion of the rebalanced strain signal to equivalent velocity units, via a scaling by local apparent velocity, allows quantitative comparison of distributed acoustic sensing and geophone data in units of velocity. We see a very good match of uncorrelated time series in both amplitude and phase, demonstrating that velocity‐converted distributed acoustic sensing data can be analyzed equivalent to vertical geophones. We show that distributed acoustic sensing data, when averaged over an interval comparable to typical geophone spacing, can obtain signal‐to‐noise ratios of 18 dB to 24 dB below clamped geophones, a result that is variable with noise spectral amplitude because the noise characteristics are not identical. With vertical seismic profile processing, we demonstrate the effectiveness of downgoing deconvolution from the large spatial sampling of distributed acoustic sensing data, along with improved upgoing reflection quality. We conclude that the extra source effort currently needed for tubing‐deployed distributed acoustic sensing vertical seismic profile, as part of a modular monitoring system, is well compensated by the extra spatial sampling and lower deployment cost as compared with conventional borehole geophones. 相似文献
14.
15.
Field tests were conducted to determine representative total-intensity magnetic anomalies due to the presence of underground storage tanks and 55-gallon steel drums. Three different drums were suspended from a non-magnetic tripod and the underlying field surveyed with each drum in an upright and a flipped plus rotated orientation. At drum-to-sensor separations of 11 feet, the anomalies had peak values of around 50 gammas and half-widths about equal to the drum-to-sensor separation. Remanent and induced magnetizations were comparable; crushing one of the drums significantly reduced both. A profile over a single underground storage tank had a 1000-gamma anomaly, which was similar to the modeled anomaly due to an infinitely long cylinder horizontally magnetized perpendicular to its axis. A profile over two adjacent tanks had a smooth 350-gamma single-peak anomaly even though models of two tanks produced dual-peaked anomalies. Demagnetization could explain why crushing a drum reduced its induced magnetization and why two adjacent tanks produced a single-peak anomaly.
A 40-acre abandoned landfill was surveyed on a 50- by 100-foot rectangular grid and along several detailed profiles. The observed field had broad positive and negative anomalies that were similar to modeled anomalies due to thickness variations in a layer of uniformly magnetized material. It was not comparable to the anomalies due to induced magnetization in multiple, randomly located, randomly sized, independent spheres, suggesting that demagnetization may have limited the effective susceptibility of the landfill material. A different 6-acre site survey conducted on a 10- by 10-foot grid was analyzed to determine the maximum station spacing and line separation that could have been used. Essentially, all of the anomalies at this site would have been resolved by a survey conducted on a 20- by 20-foot grid and the larger anomalies would have been detected by a 50- by 50-foot grid. 相似文献
A 40-acre abandoned landfill was surveyed on a 50- by 100-foot rectangular grid and along several detailed profiles. The observed field had broad positive and negative anomalies that were similar to modeled anomalies due to thickness variations in a layer of uniformly magnetized material. It was not comparable to the anomalies due to induced magnetization in multiple, randomly located, randomly sized, independent spheres, suggesting that demagnetization may have limited the effective susceptibility of the landfill material. A different 6-acre site survey conducted on a 10- by 10-foot grid was analyzed to determine the maximum station spacing and line separation that could have been used. Essentially, all of the anomalies at this site would have been resolved by a survey conducted on a 20- by 20-foot grid and the larger anomalies would have been detected by a 50- by 50-foot grid. 相似文献
16.
The depth of causative bodies may be derived from the power spectrum of their magnetic anomalies. After a short review of the theoretical basis of the method, several examples of its application to synthetic cases are shown. Disturbing effects due to improper choice of the sampling interval and to anomalies only partially contained in the examined segment of the profile are studied. The spectral method is then applied to real cases; the tapering of the anomalies, the application of non-linear filters and the effects of anomaly superpositions are investigated. In conclusion, an appraisal of the method and of its possible practical impact has been given. 相似文献
17.
CHRISTINA KRAJEWSKI LOTHAR DRESEN CHRISTOPH GELBKE HORST RÜTER 《Geophysical Prospecting》1989,37(7):717-751
The possibilities for reconstructing seismic velocity distributions containing low-velocity anomalies by iterative tomographic methods are examined studying numerical and analogue 2D model data. The geometrical conditions of the model series were designed to generalize the geometrical characteristics of a typical cross-hole tomographic field case. Models with high (30%) and low (8%) velocity contrasts were realized. Traveltimes of 2D ultrasonic P-waves, determined for a dense net of raypaths across each model, form the analogue data set. The numerical data consists of traveltimes calculated along straight raypaths. Additionally, a set of curved-ray traveltimes was calculated for a smoothed version of the high-contrast model. The Simultaneous Iterative Reconstruction Technique (SIRT) was chosen from the various tomographic inversion methods. The abilities of this standard procedure are studied using the low-contrast model data. The investigations concentrate on the resolving power concerning geometry and velocity, and on the effects caused by erroneous data due to noise or a finite time precision. The grid spacing and the source and receiver patterns are modified. Smoothing and slowness constraints were tested. The inversion of high-contrast analogue model data shows that curved raypaths have to be considered. Hence, a ray-tracing algorithm using velocity gradients was developed, based on the grid structure of the tomographic inversion. This algorithm is included in the SIRT-process and the improvements concerning anomaly localization, resolution and velocity reconstruction are demonstrated. Since curved-ray tomography is time-consuming compared with straight-ray SIRT, it is necessary to consider the effects of grid spacing, ray density, slowness constraints and the 相似文献
18.
S. E. HJELT 《Geophysical Prospecting》1973,21(2):243-265
Various experiments are described in designing two-dimensional magnetic interpretation algorithms using computer curve fitting techniques. For a single anomaly the position of the anomaly maximum and the half-width of the anomaly give good initial estimates of the plate position and thickness. A nomogram and formulae for improving these estimates is given. Curves and estimates for the effects of finite depth extent of a plate show the limits, when the lower surface of the plate can be neglected in curve fitting. The combined anomaly of two parallel plates can be separated into partial anomalies with no common points using the horizontal derivative of the anomaly. The changes of the anomaly maxima and changes in anomaly half-widths are studied as a function of plate separation. The position of the maxima and the half-widths can be corrected before applying the one-plate procedure for obtaining initial estimates of plate positions and thicknesses. The performance of standard optimization methods of Powell, Davidon, and Marquardt in improving the values of the plate parameters are compared. The Powell method seems to be the most reliable for both single and multi-plate anomalies. All methods become unacceptably slow when the number of plates is greater than 2 or 3. In these cases feasible interpretation times are obtained using the partial anomalies and sequential parabolic search of the parameter values as tailored specially to the thick plate model. Experiments with three different error norms, the classical least squares, weighted least squares and minimax, show that the first norm gives the best overall performance in automatic interpretation. The behaviour of the classical least squares norm as a function of the plate parameters is also briefly described. 相似文献
19.
以诸城—宜川深地震测深速度剖面为约束, 对沿该剖面得到的高精度重力数据进行拟合, 并对临汾强震区平面布格重力异常进行处理, 得到了该地区地壳密度结构及平面重力异常分布． 利用上述结果分析了临汾强震区的地壳结构及构造环境, 结合前人相关研究成果, 认为临汾强震区地壳中存在塑性相对较强的介质, 洪洞和临汾两次历史地震皆发生在其与周边弹性介质的转换边界上． 另外, 临汾凹陷南北两侧局部构造环境存在差异, 在区域应力场作用下, 导致了洪洞地震和临汾地震的发生. 两次地震在发震时间、 地点和震级等地震要素上有所不同． 相似文献
20.
A new unbiased plotting position formula for Gumbel distribution 总被引：1，自引：0，他引：1
The probability plots (graphical approach) are used to fit the probability distribution to given series, to identify the
outliers and to assess goodness of fit. The graphical approach requires probability of exceedence or non exceedence of various
events. This is obtained through the use of plotting position formula. In literature many plotting position formulae have
been reported. All of the many existing formulae provide different results particularly at the tails of the distribution and
hence there is need of unbiased plotting position formulae for different distributions. Expression for the largest expected
order statistics is found in a simple form. Using exact plotting position from Gumbel order statistics a new unbiased plotting
position formula has been developed for the Gumbel distribution. The developed formula better approximates the exact plotting
positions as compared to other existing formulae. 相似文献