首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
The procedure of polynomial fitting by least-squares technique to obtain the regional field from potential field data is commonly used in processing geophysical data. However, it has been found that the unguarded use of this procedure with measured field data may lead to erroneous regional fields. This is due to the fact that the least-squares technique pre-supposes a random variation of the residual or anomaly field. This condition is rarely fulfilled in measured field data where there is the common occurrence of ‘deviatory field points’ in the residual field data. In this paper, a procedure to overcome the problem using robust statistics is described with the straight-forward Downhill Simplex method used to obtain the optimum coefficients for the polynomial. The method is illustrated with aeromagnetic field data from the Mamfe basin of Nigeria and Cameroon. The result shows that the trend of the regional field obtained from robust statistics is consistent with the trend in the aeromagnetic field obtained by using the method of cross-correlation trend analysis.  相似文献   

2.
Gravity data, integrated with seismic refraction/reflection data, well data and geological investigations, were used to determine the location of the paleogeographic boundary between the Precambrian Saharan domain and the younger Tunisian Atlas domain. This boundary (North Saharan Flexure or NSF) has not been as clearly defined as it has been to the west in Algeria and Morocco. The gravity data analysis, which included the construction of complete Bouguer and residual gravity anomaly maps, revealed that the Atlasic domain is characterized by relative negative gravity anomalies and numerous linear gravity trends implying a thick and deformed sediment cover. The Saharan domain is characterized by relatively positive gravity anomalies with few gravity trends implying a thin and relatively undeformed sediment cover. An edge-enhancement analysis of the residual gravity anomalies revealed that the NSF is characterized by a series of discontinuous east- and northwest-trending linear anomalies south of 34°N that are not related to the well-known faults within the Gafsa and Accident de Medenine regions. Based on the continuity of the amplitudes of seismic reflection data and the trends of the residual gravity anomalies, the NSF is not an abrupt discontinuity but a series of step faults dipping toward the Atlasic domain. To obtain a more quantitative representation of the southern edge of Tunisian Atlas, a regional gravity model constrained by two wells and seismic reflection/refraction data was constructed along a north-south trending profile which confirms the presence of thicker sediments north of the NSF. Our analysis shows that the NSF has controlled the depositional environment of the sedimentary rocks within the region since at least Triassic time and has acted as a barrier to Atlasic deformation south of the NSF. The NSF is considered an important tectonic feature that has controlled the paleogeographic evolution of the southern margin of the Tethys Ocean, and it continues to be active today based on seismicity hazard studies.  相似文献   

3.
A gravimetric survey, covering a site 200 m square, was carried out in order to locate karstic cavities. After eliminating the regional trend using a polynomial fit, the residual is modelled by least-squares prediction. Correlated signals for several wavelengths are detected. The inversion of these anomalies is performed by a global 3D adjustment using spherical bodies as models. The adjustment is repeated in order to obtain a stable configuration. The results show the probable presence of a system of cavities and galleries. Data collected from boreholes and the subsequent appearance of sink-holes are consistent with the results.  相似文献   

4.
The effectiveness of a standard polynomial trend removal technique and two-dimensional wavenumber filtering in isolating gravity anomalies due to buried preglacial drainage systems was evaluated for an area in northwestern Pennsylvania. It was found that gravity surveying is a useful reconnaissance technique in this area. Both analytical methods isolated anomalies due to valleys of major hydrological interest, provided similar useful information on the valley trends, and gave reasonably accurate width estimates. Accurate depth estimates could not be inferred, although large anomalies generally coincided with thicker drift accumulations. With the station spacing used (about 0.75 miles), valleys sought must have drift accumulations in excess of 75 feet and widths of at least one mile if they are to be adequately resolved. The two-dimensional filtering method seems to offer the most complete information, while still being practical, and it avoids many of the problems of the polynomial method.  相似文献   

5.
Detrending is a widely used technique for obtaining stationary time series data in residual analysis and risk assessment. The technique is frequently applied in crop yield risk assessment and insurance ratings. Although several trend models have been proposed in the literature, whether these models achieve consistent detrending results and successfully extract the true yield trends is rarely discussed. In the present article, crop insurance pricing is evaluated by different trend models using real and historical yield data, and hypothetical yield data generated by Monte Carlo simulations. Applied to real historical data, the linear, loglinear, autoregressive integrated moving average trend models produce different risk assessment results. The differences among the model outputs are statistically significant. The largest deviation in the county crop assessment reaches 6–8 %, substantially larger than the present countrywide gross premium rate of 5–7 %. In performance tests on simulated yield trends, popular detrending methods based on smoothing techniques proved overall superior to linear, loglinear, and integrated autoregression models. The best performances were yielded by the moving average and robust locally weighted regression models.  相似文献   

6.
This paper firstly discusses the feasibility of delineating the volcanic rocks distribution by gravitational and magnetic methods on the basis of the statistical results of the magnetic susceptibility and density of the different lithologies in the study area. After the separation of gravitational and magnetic fields by bandpass filter, we determined the residual gravitational and magnetic anomalies caused by volcanic rocks. The results of potential field separation show that the residual anomalies are in beaded NE-directed distribution. In the meantime, the boundary enhancement techniques such as horizontal total gradient and vertical derivative are employed to delineate the distribution of the faults, as the result, fifteen faults including four NE-striking main faults which appear as an arc protruding to the southeast and other secondary faults are distinguished. Furthermore, two fault systems with deep and superficial attribution are revealed from the inversed section of the telluric electromagnetic sounding, and their properties, characteristics and roles in the volcanic activity are fully discussed by combining with the tectonic background. Based on the comprehensive analysis of the correlation between the distribution characteristics of the residual gravitational and magnetic anomalies and the location of the faults, three volcanic activity zones are reasonably delineated. Finally, we carry out the inversion of apparent density by taking advantage of the residual gravitational anomalies and acquire the apparent density anomalies of the top part of Paleozoic. Integrated with the apparent density anomalies, the lithologies exposed from the boreholes, the observed density data of different rocks and the residual gravitational-magnetic anomalies, we elaborately delineate the lithologies of the top of Paleozoic and further define the distribution range and the lithology combination of Kekesayi group of Middle Ordovician, Tailegula group of Lower Carboniferous and Jiamuhe group of Lower Permian.  相似文献   

7.
A method for estimating the degree of polynomial fitted to gravity anomalies to evaluate the regional anomaly is presented. The anomaly can be fitted by polynomials of different degrees with the least-squares method and the optimum degree of the polynomial evaluating this anomaly can be estimated from the point of discontinuity of the gradient on a graph of variance against the polynomial degree. The Bouguer gravity is initially separated by upward continuation to a proper height and then the degree of regional polynomial to fit the Bouguer anomaly can be estimated. Theoretical and field examples show the effectiveness of the method.  相似文献   

8.
The decrease of density contrast in sedimentary basins may be approximated by a quadratic function. A sedimentary basin may be viewed as a number of prisms placed in juxtaposition. Equations in closed form for the gravity anomalies of 3D and 2½ D prismatic models are derived. Approximate equations for these models are also derived for rapid calculations. Efficient methods are developed for anomaly calculation by an appropriate use of the exact and approximate equations, and hence, for 3D and 2½ D modelling. The depths to the basement are adjusted iteratively by comparing the calculated anomalies with the observed anomalies. These methods are applied for analysis of the residual anomaly map of the Los Angeles Basin, California.  相似文献   

9.
Magnetic anomalies of complicated 3D sources can be calculated by using a combination of analytical and numerical integration. Two surfaces and the magnetization parameters (the amplitudes of the induced and remanent components and the direction cosines) of the source can be defined by arbitrary functions or by discrete data points in a plane. When combined with a polynomial magnetization function in the direction of the third axis, 3D magnetization distribution can also be modelled. The method gives very general equations for anomaly calculation. It can be used for direct modelling of sources interpreted by seismic or other methods and also for interactive interpretation with fast computers. It is possible to calculate anomalies of, for example, intrusives or folded sedimentary beds whose surfaces are functions of horizontal coordinates and which have polynomial magnetization variations in the vertical direction due to gravitational differentiation and arbitrarily varying magnetization in the horizontal direction due to regional metamorphosis. If the distribution of magnetization parameters in the vertical direction cannot be described satisfactorily by polynomials, models can be used whose surfaces are functions of the vertical coordinate and which can then have any arbitrary magnetization distribution in the vertical direction.  相似文献   

10.
Regional–residual separation is essential in gravity and magnetic data interpretation and a variety of techniques have been proposed. Graphical determination of the regional allows geological information to be taken into account. Upward continuation can be used to obtain the regional field either empirically or using some hypothesis about the geology. In some cases, a matched filter can be designed and used to separate deep and shallow sources. Simple low pass filtering has also been used but without much success. Here we propose to use a non-linear filter approach to remove gravity and magnetic anomalies smaller than a given width. This technique attempts to mimic the graphical separation method. The results from synthetic models are presented as well as the results from a case study in eastern Canada and compared to regional gravity and magnetic anomalies obtained by other techniques. Contrary to the regional fields obtained by upward continuation, non-linear filtering does not have any physical meaning. However, its main advantage is that it gives a regional component of the gravity or magnetic field similar to the one obtained from a graphical separation.  相似文献   

11.
The methods of anomaly transformations considered are based on a system of combined analysis of the geophysical field and a priori) information on the structure of a geological object. The methods involve calculation of a transformative polynomial (describing geophysical noise) which makes it possible to separate the residual field component related to the geological characteristic under study in a correlatively optimal way. The structure of the transformative polynomial is determined by the nature of the geophysical noise that is eliminated by the field transformation. Various correlation methods of anomaly transformations arise, depending on the structure of the transformative polynomial chosen. By way of example, the correlation method employed for separating the geophysical anomalies is shown to be highly effective in investigating the local geological structure.  相似文献   

12.
—Methods and the results of estimating the anomalies characterising the density inhomo geneities in the European-Mediterranean upper mantle are described. These anomalies were obtained by subtracting the gravity effect of a crustal density model derived from seismic velocities from the observed gravity field averaging over an area of 1°× 1°. The 3-D density model of the study region comprises two regional layers of varying thickness with lateral variation of average density the sedimentary cover and the crystalline crust. The average densities for model layers were evaluated by using a velocity/density conversion function and taking into account sediment consolidation with depth. Clear correlation between residual gravity anomalies and both velocity heterogeneities and thermal regime data of the upper mantle has been revealed. An agreement of positive anomalies over the Alps, the Adriatic plate and the Calabrian Arc with high velocity domains in the upper mantle and reduced temperatures at the subcrustal layer are caused by lithospheric "roots" and thickened lithosphere below these structures. Gravity residual lows, revealed over the Western Mediterranean Basin and Pannonian Basin, are in correspondence with both low velocities and high temperatures in the upper mantle. These anomalies are the result of the presence of asthenosphere in shallow near-Moho depths below these basins.  相似文献   

13.
地磁总强度资料在地震预报中的应用   总被引:1,自引:0,他引:1       下载免费PDF全文
本文应用地磁总强度资料,用两台差消除地磁场短周期变化,同时亦削弱了长期变影响,从而使年变呈线性变化。然后采用趋势面法求得剩余长期变地区性改正量。资料处理后,显示唐山大震前可能有一年半以上的异常,变化幅度约10伽马,对台网中其它几个6级左右地震亦有较好的异常反应。  相似文献   

14.
15.
Aeromagnetic data collected in areas with severe diurnal magnetic variations (auroral zones) are difficult to level. This paper describes levelling of an aeromagnetic survey where such conditions prevail, and where sophisticated levelling techniques are needed. Corrections based on piecewise low‐order polynomial functions are often used to minimize mis‐ties in aeromagnetic data. We review this technique and describe similar mis‐tie fitting methods based on low‐pass filter levelling, tensioned B‐spline levelling and median levelling. It is demonstrated that polynomial levelling, low‐pass filter levelling and tensioned B‐spline levelling depend on the careful editing of outlying mis‐ties to avoid the introduction of false anomalies. These three techniques are equally efficient at removing level errors. Median levelling also removes level errors efficiently, but it is more robust in the sense that mis‐tie editing is not required. This is due to the inherent noise‐removal capabilities of the median filter. After mis‐tie editing, the total field anomalies of the other three techniques closely resemble the unedited median‐levelled total field anomaly.  相似文献   

16.
建立了区域重力均衡异常的物理模型,假定区域重力均衡异常源于上地幔密度分布横向不均匀以及由上地幔小尺度对流产生的边界形变.导出了区域重力均衡异常和上地幔小尺度对流之间的相关方程,讨论了该方程在地球动力学研究中的应用.应用几种模型计算上地幔小尺度对流产生的作用于岩石层底部的拖曳力场.结果表明,将这一相关方程用于利用区域重力均衡异常研究地壳及上地幔应力状态,将有广阔的前景.  相似文献   

17.
重力异常是地下不同规模、不同形态和不同埋深的不均匀地质体的综合响应,重力勘探主要通过从重力异常中提取感兴趣的局部异常,以便探测深部结构,寻找隐伏矿床.为探讨重力异常分离原则并检验方法效果,本文从各方法原理入手,加上模型试验以及在安徽省泥河矿区深部隐伏铁矿的探寻实践检验,阐明:趋势分析法是整体拟合不同于最小二乘圆滑的局部拟合,由于是多项式拟合区域场,趋势分析法不适宜范围大、地质情况复杂的测区;插值切割法以计算点场值与四点圆周平均值的插值运算为切割算子,通过连续切割,得到重力异常的切割区域场和局部场,插值切割法对于小测区单个异常的分离效果较好,切割次数选择1到2次即可;匹配滤波法通过分析实测异常功率谱曲线、选择合适的滤波段、建造适宜的低通和带通滤波器进行滤波,从而提取不同波数成分的异常场,匹配滤波更适合垂向叠加的异常分离;解析延拓是根据一个面上的一组位场数据确定另一个不同高度面上位场值,应用中要把握延拓高度;垂向二阶导数法可以起到突出浅源异常,区分水平叠加异常,确定异常体的边界,消除或削弱背景场的作用.通过安徽泥河铁矿重力异常分离实验,发现三阶趋势分析、向上延拓以及插值切割法能很好地分离出矿体异常和背景场,同时发现在泥河矿区东南部和东北部还存在剩余重力异常,可为泥河铁矿扩大规模提供新的线索.  相似文献   

18.
19.
We use the continuous wavelet transform based on complex Morlet wavelets, which has been developed to estimate the source distribution of potential fields. For magnetic anomalies of adjacent sources, they always superimpose upon each other in space and wavenumber, making the identification of magnetic sources problematic. Therefore, a scale normalization factor, a?n, is introduced on the wavelet coefficients to improve resolution in the scalogram. By theoretical modelling, we set up an approximate linear relationship between the pseudo‐wavenumber and source depth. The influences of background field, random noise and magnetization inclination on the continuous wavelet transform of magnetic anomalies are also discussed and compared with the short‐time Fourier transform results. Synthetic examples indicate that the regional trend has little effect on our method, while the influence of random noise is mainly imposed on shallower sources with higher wavenumbers. The source horizontal position will be affected by the change of magnetization direction, whereas the source depth remains unchanged. After discussing the performance of our method by showing the results of various synthetic tests, we use this method on the aeromagnetic data of the Huanghua depression in central China to define the distribution of volcanic rocks. The spectrum slices in different scales are used to determine horizontal positions of volcanic rocks and their source depths are estimated from the modulus maxima of complex coefficients, which is in good accordance with drilling results.  相似文献   

20.
In this study, a regional scale gravity data set has been inverted to infer the structure (topography) of the top of the basement underlying sub‐horizontal strata. We apply our method to this real data set for further proof of concept, validation and benchmarking against results from an earlier forward modelling done elsewhere. Our aim is to carry out implicit structural inversion, i.e., to obtain a geologically reasonable model, without specifically solving for structure. The 2.5D volume of interest is parametrized with homogeneous horizontal prisms and a two‐lithology medium is assumed. A possible regional linear trend and a general floating reference are also inverted for. Using a gridded parametrization, linear programming is used to minimize the L1 ‐norm of the data misfit, relative to a floating reference. Given a known density contrast between the lithologies, an inversion using linear programming has the intrinsic advantage that a relatively sharp image of the sub‐surface is retrieved instead of a smooth one. The model recovered is almost bi‐modal and its general features seem to be robust with respect to several parametrization scenarios investigated. The floating reference and a linear trend in the data were also retrieved simultaneously. The inversion results, indicating two depressions in the basement, are robust and agree with those obtained earlier based upon detailed 2D forward modelling using many narrow, near‐vertical prisms.  相似文献   

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

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