共查询到20条相似文献，搜索用时 593 毫秒
1.
The decrease of density contrast in sedimentary strata may be approximated by a quadratic function. A sedimentary basin may be considered as a number of horizontal polygonal prisms of finite thickness placed one over the other. Equations for the gravity anomaly of ann-sided polygon prism have been derived using a quadratic density function. A method has been developed for inversion of the gravity anomalies using a polygonal prismatic model. An analytical method has been employed for evaluating the necessary derivatives as the computing time by this method is much slower than by a numerical method. Approximate equations have been derived for rapid calculation of the anomalies and derivatives. Efficient computer programs have been developed for calculation of the anomalies and derivatives by an appropriate use of the exact and approximate equations. The depths to the polygon prisms constituting the basin are iteratively adjusted by minimizing the sum of the squares of the differences between the observed and calculated anomalies. These methods have been applied to analyze the Bouguer anomaly map of the Los Angeles basin, California. 相似文献
2.
V. Chakravarthi 《Pure and Applied Geophysics》1995,145(2):327-335
The decrease in density contrast of sedimentary rocks with depth in many sedimentary basins can be approximated by a parabolic density function. Analytical gravity expression of an outcropping two-dimensional vertical step along which the density contrast decreases parabolically with depth is derived in the space domain. A modification ofBott's (1960) method of gravity interpretation is proposed by considering two outcropping vertical steps on either side of the first and last observation points in addition toN outcropping vertical prisms in order to interpret the gravity anomalies of nonoutcropping basins. The thicknesses of the two outcropping vertical steps are made equal to the thicknesses of the two outcropping vertical prisms placed below the first and last observation points. The initial depth estimates of the sedimentary basin are calculated by the infinite slab formula ofVisweswara Rao
et al. (1993). The gravity effects of theN outcropping prisms and the two outcropping vertical steps are calculated at each anomaly point and the depth to the floor of the basin are adjusted based on the differences between the observed and calculated anomalies. A gravity anomaly profile of Los Angeles basin, California is interpreted. 相似文献
3.
In this study, a Markov Random Field (MRF) approach is used to locate source boundary positions which are difficult to identify
from Bouguer gravity and magnetic maps. As a generalized form of Markov Chains, the MRF approach is an unsupervised statistical
model based algorithm and is applied to the analysis of images, particularly in the detection of visual patterns or textures.
Here, we present a dynamic programming based on the MRF approach for boundary detection of noisy and super-positioned potential
anomalies, which are produced by various geological structures. In the MRF method, gravity and magnetic maps are considered
as two-dimensional (2-D) images with a matrix composed of N1 × N2 pixels. Each pixel value of the matrix is optimized in real time with no a priori processing by using two parameter sets; average steering vector (θ) and quantization level (M). They carry information about the correlation of neighboring pixels and the locality of their connections. We have chosen
MRF as a processing approach for geophysical data since it is an unsupervised, efficient model for image enhancement, border
detection and separation of 2-D potential anomalies. The main benefit of MRF is that an average steering vector and a quantization
level are enough in evaluation of the potential anomaly maps. We have compared the MRF method to noise implemented synthetic
potential field anomalies. After satisfactory results were found, the method has been applied to gravity and magnetic anomaly
maps of Gelibolu Peninsula in Western Turkey. Here, we have observed Anafartalar thrust fault and another parallel fault northwest
of Anafartalar thrust fault. We have modeled a geological structure including a lateral fault, which results in a higher susceptibility
and anomaly amplitude increment. We have shown that the MRF method is effective to detect the broad-scale geological structures
in the Gelibolu Peninsula, and thus to delineate the complex tectonic structure of Gelibolu Peninsula. 相似文献
4.
MAHBOUB AL-CHALABI 《Geophysical Prospecting》1972,20(1):1-16
The interpretation of a gravity anomaly in terms of the shape of the anomalous body is a non-linear problem and may, therefore, be carried out using non-linear optimisation techniques. The formulation is extended to include cases where the density contrast and the regional background are also unknown. For a given model the objective function is provided by the discrepancy between the observed anomaly and the calculated anomaly due to the model. Given an initial model, the optimisation procedure searches for a minimum of the objective function by an iterative adjustment of the parameters. A number of suitable objective functions is given. The behaviour of these functions in the parameter hyperspace is quite complex. Accordingly, direct search methods should be employed at the early stages of the search, changing to gradient methods at later stages. The use of constraints is also necessary to ensure the geological feasibility of the model. The required computer time may be largely reduced by careful programming. Two examples of interpretation by optimisation methods are given. 相似文献
5.
A constrained 3D density model of the upper crust along a part of the Deccan Syneclise is carried out based on the complete Bouguer anomaly data. Spectral analysis of the complete Bouguer gravity anomaly map of the study region suggests two major sources: short wavelength anomalies (<100 km) caused primarily due to the density inhomogeneities at shallow crustal level and long wavelength anomalies (>100 km) produced due to the sources deeper than the upper crust. A residual map of the short wavelength anomalies is prepared from the complete Bouguer anomaly using Butterworth high‐pass filter (100 km cut‐off wavelength). Utilizing the constraints from deep resistivity sounding, magnetotellurics and deep seismic sounding studies, 2.5D density models have been generated along 39 profiles of this region. The mismatch between the calculated response of the a priori 2.5D model with the residual (short wavelength) gravity anomalies is minimized by introducing high‐density intrusive bodies (≥2.81 g/cm3) in the basement. With these 2.5D density models, the initial geometry of our 3D density model, which includes alluvium, Deccan trap, Mesozoic sediment and high‐density intrusive bodies in the basement up to a depth of 7 km (upper crust), is generated. In the final 3D model, Deccan trap extends from 200 m to nearly 1700 m below the 90–150 m thick Quaternary sediment. Further down, the sub‐trappean Mesozoic sediment is present at a depth range of 600–3000 m followed by the basement. The derived 3D density model also indicates six intrusive bodies of density 2.83 g/cm3 in the basement at an average depth of about 4–7 km that best fits the residual gravity anomaly of the study area. 相似文献
6.
In this paper, we attempt to use satellite gravity data and a new inversion method to study the lateral density anomaly distribution
in the mantle.
First, density difference Δρ(τ,θ,φ) is expanded in terms of a three—dimensional orthogonal function system, the coefficients of the expansion are to be determined.
Then, a set of observation equations is established from the relationship between density anomaly and disturbing geopotential.
In the equations the unknown vector contains the coefficients of density anomaly expansion, the observational vector is obtained
by computing geopotential perturbations using the potential coefficients of GEM10B, and a filtering process is done for the
observational values by properly selecting the harmonic degrees of geopotentical. Finally, the lateral density variations
in the lower mantle (670 km toCM boundary) are investigated. In this case, the degrees of disturbing geopotential are selected as 2–11, the truncated degrees
of density anomaly expansion are taken asL=6 andK=4, and the damping least squares method is used to solve the observation equations.
The resulting model shows the high level of density perturbations at 670 km discontinuity and core — mantle boundary, a high
— density zone circumscribing the Pacific and a lower—density region under the center of Pacific. These features are in agreement
with the three—dimensional seismic velocity variation features by Dziewonski (1984). In the Antarctic region and some parts
of Atlantic and Indian Ocean, however, the resulting density anomalies are negatively correlated with the seismic velocity
anomalies, the cause resulting in these phenomena is preliminarily analysed in this paper.
The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 53–65, 1991.
The principle and method represented in this paper can also be suitable to study the lateral density anomaly distribution
in the earth’s crust and the upper mantle. 相似文献
7.
I. V. Radhakrishna Murthy P. Rama Rao P. Ramakrishna 《Pure and Applied Geophysics》1989,130(4):711-719
The exponential decrease in the density contrast of anomalous bodies is equated to a quadratic function. The three-dimensional body is divided into a series of parallel vertical polygonal sections of unit thickness each. The gravity effect of each side of the vertical polygonal cross-section with a quadratic density variation is derived in closed form. Gravity effects of a few equispaced and representative cross-sections are integrated to obtain the gravity anomaly of the entire body. When used in modelling gravity data, the method of considering vertical cross-sections has the unique advantage of easily allowing changes in theZ-coordinates of polygonal sections. 相似文献
8.
An equivalent layer magnetization model obtained from inversion of long-wavelength satellite magnetic anomaly data indicates a very magnetic source region centered in south central Kentucky. The magnetization maximum nearly coincides with a gravity high elongated north-south and extending into Tennessee. Previous refraction profiles suggest that the source of the gravity anomaly is a large mass of rock occupying much of the crustal thickness. The outline of the source delineated by gravity contours is also discernible in aeromagnetic anomaly patterns. Taken together, the geophysical data suggest a large, localized mass of intracrustal rock which is both dense and very magnetic. A simple magnetization/density model is given which accounts for the gravity and long-wavelength aeromagnetic anomalies due to the body. We interpret it as a mafic plutonic complex, and several lines of evidence are consistent with a rift association. The body is, however, clearly related to the inferred position of the Grenville Front. It is bounded on the north by the fault zones of the 38th Parallel Lineament. The inferred mean magnetization (4 A/m) of the body is large, but not inconsistent with values reported by others for deep crustal bodies associated with long-wavelength magnetic anomalies. Such magnetization levels can be achieved with magnetic mineralogies produced by normal oxidation and metamorphic processes and enhanced by viscous build-up, especially in mafic rocks of alkaline character. 相似文献
9.
南海位于太平洋板块、印澳板块和欧亚板块交汇处,自晚中生代以来历经张裂作用、海底扩张以及印藏碰撞、菲律宾海板块西向运动等构造事件的叠加改造,不仅形成了复杂多样的构造格局,而且堆积了厚薄不均的沉积层.为了考察沉积层密度改正对利用重力资料分析南海不同尺度构造特征的影响,本文利用南海各区域不同深度沉积层的地震波速度及钻孔密度等数据,建立了沉积层与沉积基底密度差随深度变化的二次函数关系式,并基于该关系式,计算了南海沉积层相对基底密度低而产生的重力异常值.结果显示,南海沉积层的重力异常值在海盆区介于-40~-60 mGal,而在堆积巨厚沉积物的莺歌海盆地可达到-135 mGal;相对于空间重力异常、布格重力异常,经沉积层重力异常改正后的地壳布格重力异常更能突出深部不同尺度的密度结构和莫霍面的起伏特征,其总水平导数模更突显了南海西北部红河断裂带的海上延伸;利用谱分析技术估算岩石圈强度时,经沉积层重力异常改正的地壳布格重力异常数据获得的岩石圈有效弹性厚度值更为符合地质实际,特别是在长条形的巨厚沉积区如莺歌海盆地和马来盆地.分析表明,重力异常的沉积层密度改正对揭示南海构造特征具有重要的意义. 相似文献
10.
利用宁夏地区2011—2013年期间多期重力观测资料,采用中国地震局推广的LG-ADJ程序对观测资料进行平差处理,给出地震前后3年的区域重力场变化和点值变化图像。结合构造活动分析,对2012永宁4.6级地震前后宁夏地区重力异常变化进行研究。认为：2012年永宁地震前,测区重力场出现较好的中期前兆性变化图像,临震前短期异常突增,重力场在空间和量值上均出现大幅的上升变化,重力点值震前趋势下降-转折上升-加速上升-发震的变化,是非常值得关注的前兆异常。 相似文献
11.
利用哈佛大学GCMT数据中心和前人积累的历史地震资料(1962~2016年M W>4.0地震)以及Crust2.0地壳结构统计分析了喜马拉雅地区、天山地区的地壳区域构造与地震活动间的相关性。此外,利用GFZ地学研究中心提供的静态卫星重力模型GGM03S/EGM2008和地形模型Topo计算了2个地区的各类重力异常场,同时还模拟了不同地壳弹性参数下的重力异常场,结果表明喜马拉雅地区重力异常场在水平、垂直方向的梯度特征远大于天山地区的异常特征,且喜马拉雅地区的有效弹性板厚度Te(6~15km)小于天山地区的有效弹性板厚度Te(20~30km)。最后,利用喜马拉雅地区与天山地区的GPS震间三维形变场约束了断层运动模型,结果显示两者主前缘断裂的断层闭锁深度及应力积累状态存在较大的差异。因此认为,造成青藏高原及邻区的边界地壳区域地震活动性差异的动力学因素,与地壳有效弹性板厚度、孕震断层参数及区域应力积累状态等密切相关。 相似文献
12.
An inversion technique using the Marquardt optimization is developed to interpret the gravity anomalies due to anticlinal and synclinal structures with density contrast varying continuously with depth. The algorithm simultaneously estimates the parameters of the respective models, in addition to the regional gravity background that is invariably associated with the residual gravity anomaly. Forward modelling is realized through analytically derived gravity expressions for the respective models in the space domain. The efficacy of the inversion is demonstrated with the gravity anomaly due to a theoretical model, in each case with and without the regional background. In addition, the applicability is illustrated using the gravity anomalies of the Pays De Bray anticline, situated north‐west of Paris, France. The interpreted depth of the Pays De Bray anticline using the present inversion compares well with the drilling depth. 相似文献
13.
本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征, 以鲁甸和景谷地震为例, 认识其深部构造环境和动力学过程, 为该区域的构造运动和地震孕育环境研究提供依据。 结果表明, 研究区布格重力异常和均衡重力异常与地质构造格局相关性较好, 川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形、 断裂及地震活动密切相关。 强震多分布在断裂带两侧重力异常的过渡地带和高梯度带, 断裂带两侧横向和垂向的显著介质密度差异是强震孕育的深部构造背景。 布格重力异常和均衡重力异常揭示的鲁甸、 景谷震源区深浅差异性的重力异常特征, 暗示鲁甸和景谷地震孕震环境的不同。 相似文献
14.
基于EGM2008重力场模型计算获得了渭河盆地及邻区布格重力异常。采用小波多尺度分解方法对布格重力异常进行了4阶小波逼近和小波细节分解,同时基于平均径向对数功率谱方法定量化地计算出1~4阶小波细节和小波逼近所对应的场源平均埋深。结合区域地质和地震资料,对获得的重力场结果进行分析,得到如下结论:①鄂尔多斯地块、渭河盆地、秦岭造山带3个一级构造单元的布格重力异常之间存在明显差异;构造区内部重力异常也存在横向的显著差异。布格重力异常的走向、规模、分布特征与二级构造区及主要的断裂具有一定的对应关系。②渭河盆地及邻区布格重力异常1~4阶细节对应4~23 km不同深度的场源信息,鄂尔多斯地块南缘东、西部的地壳结构存在明显的差异;渭河盆地凹陷、凸起构造区边界清晰,断裂边界与重力异常边界具有较好的一致性;秦岭造山带重力异常连贯性不好,东、西部重力异常变化特征表现出明显的差异。③渭河盆地及邻区布格重力异常分布与莫霍面埋深具有非常明显的镜像关系。渭河盆地及邻区地震主要分布在六盘山—陇县—宝鸡断裂带、渭河断裂与渭南塬前断裂交汇处、韩城断裂与双泉—临猗断裂交汇处。渭河盆地及邻区重力异常主要由中上地壳剩余密度体所影响,这可能是该区地震以浅源地震为主的主要原因。 相似文献
15.
2017年九寨沟地震前重力场异常特征提取 总被引：2，自引：2，他引：0
2017年8月8日九寨沟发生7.0级地震,此次地震发生在塔藏断裂、岷江断裂和虎牙断裂附近,震源机制为走滑型。重力场静态异常和动态异常变化可为研究深地壳结构性质提供重要的物理场信息。本文利用WGM计算的布格重力异常数据及实测流动重力数据进行小波多尺度分解,分离得到不同深度场源特征的地壳介质横向不均匀性。结果表明,WGM布格重力异常数据和实测重力数据的小波多尺度结果较为一致,通过对比分析发现,九寨沟地震发生在区域重力场的四象限分布中心位置,不同时空尺度的重力场变化对于深入认识潜在地震危险性趋势具有一定的科学意义。 相似文献
16.
The method of downward continuation is well known to those working in gravity, magnetic, SP and low-frequency electromagnetic exploration. It is demonstrated that the method of continuation can also be usefully employed in the interpretation of induced polarization gradient profiling using point electrodes to determine target depth. The apparent resistance Ra and chargeability Ma measurements obtained with point electrode excitation of the ground have been used to compute the values of (Ra)l and (Ma)l that would be obtained with a linear array. Continuation of the apparent polarizability profile thus obtained with the linear array gives a value for the depth of the target which agrees closely with that obtained by the continuation of the SP profile. On the other hand, continuation of the profile of the secondary transient signal (VS)L alone, yields a depth of the target which is in agreement with the borehole information. However, it is seen that the secondary transient voltage profiling response splits into two anomalies which fall on either side of the SP and/or (Ma)l anomaly centre, and does not coincide with that of the latter. 相似文献
17.
The MS7.0 Jiuzhaigou earthquake occurred on August 8, 2017. The earthquake occurred in the vicinity of the Tazang fault, the Minjiang fault and the Huya fault, where the focal mechanism is of the strike slip type. The static and dynamic anomalies of the gravity field can provide important physical field information for studying the structural properties of deep crust. Multi-scale decomposition techniques are used to separate Bouguer gravity at different depths and give some explanation to gravity variations at different time space scales.The results indicate that the wavelet multi-scale results of the EGM2008 model and the measured gravity data are consistent. Through comparative analysis, it is found that the Jiuzhaigou earthquake occurred in the stress enhanced region. The variation of gravity field at different time scales has a certain scientific significance for further understanding potential earthquake risk trend. 相似文献
18.
19.
青藏高原东北缘重力异常多尺度横向构造分析 总被引：8，自引：6，他引：2
本文研究了青藏高原东北缘地区布格重力异常特征,采用优化滤波法和归一化总水平梯度垂直导数法对研究区重力异常进行多尺度分离和横向构造分析.分离出的多尺度重力异常特征表明:1) 青藏高原东北缘地区大致以东经106°线为界,有一条醒目的重力异常梯级带,即贺兰山-六盘山-川滇南北构造带的北段,其东西两侧布格重力异常特征在形态和走向上截然不同,意味着两侧密度结构和构造特征存在明显差异. 2) 鄂尔多斯地块内部定边以北,重力异常高带走向由北东向转为近南北向,推测定边附近存在一个密度或构造界面,其两侧物质组成和构造特征具有差异,对比大尺度重力异常和中尺度重力异常,表明异常特征的这种差异主要是由上地幔深部结构引起的. 3) 青藏高原东北部各块体深部边界位置与地表构造分布不同,反映出该区构造复杂,深浅构造差异大. 4) 由于印度-欧亚板块碰撞及随后印度板块持续向北的挤压作用,造成青藏高原东北缘中、下地壳物质在巨大的北东向推挤力和鄂尔多斯刚性块体阻挡的共同作用下,沿着相对软弱的秦岭造山带方向蠕动.依据多尺度重力异常及其横向构造特征,综合推断出研究区内五条断裂带,即秦岭地轴北缘断裂带、海原-六盘山断裂带、香山-天景山断裂带、烟筒山断裂带和青铜峡-固原断裂带,并分析了它们在地壳深部的可能展布特征. 相似文献
20.
The magnetic anomaly due to a long tabular body usually consists of a maximum and a minimum. The distances and the amplitudes of the maximum and the minimum, when defined in dimensionless quantities, may be used as characteristics of the source. In this paper, a method based on the positions of the maximum and the minimum on the magnetic anomaly due to a long tabular body has been presented. Characteristic ratios,D andA involving the distances and amplitudes of the maximum and the minimum points on the anomaly curve are defined. Nomograms showing the variations ofD andA with the parameters of (1) the dike and (2) the vertical fault models are presented. The parameters of the causative source are evaluated from the two ratiosD andA and the nomograms, using some simple analytical relations presented here. From the nomograms, it is observed that (a) for a thick dike,A is always greater thanD, (b)A=D for a thin sheet and (c) for a vertical fault,A is always less thanD. Thus from the characteristic ratiosD andA it is possible to evaluate the source parameters and also to distinguish whether the source is a dike, sheet or a vertical fault. The method is fast and is applicable for the magnetic anomalies either in total, vertical or horizontal component. The method has been applied on two field examples and the results are found to be in close agreement with those obtained by using other methods. A simple method of locating the origin on the anomaly curve is included. The limitations of the method are also discussed. 相似文献