首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 359 毫秒
1.
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.  相似文献   

2.
常规的重力梯度法对地质体位置划分很模糊,当深度越大,地质体越小时,对Vzz和Vzx求导的误差也越大,导致结果与实际偏差很大。本文提出状态判别因子对角点位置进行优化确定,比常规高阶重力导数法更精确,使重力梯度法在对小地质体和断裂上有更好的分辨率,理论模型试算和实际资料处理效果很好,为找矿或划分地下地质体产状提供了更有利的依据。  相似文献   

3.
A method is presented for determining bounds of the properties of axial symmetric bodies from a finite number of gravity and magnetic observations based on Parker's theory of ideal bodies. Bounds on the density contrast and the intensity of magnetization are calculated as a function of depth to the top of the anomalous source, restricting the range of smallest possible solutions to fit the data. The model studied is approximated by an array of vertical annuli cylinders, each of uniform density and magnetization. Linear programming algorithms based on the ideal body theory were used to calculate the distribution of these parameters within the body. Simultaneous inversion of gravity and magnetic data is performed assuming a constant ratio between the density contrast and the intensity of magnetization and that a common body is responsible for both observed fields. The parameter k(|J|/δp) provides information about the rock type of the structure. Interpretation of gravity and aeromagnetic data from Darnley Bay, NWT, Canada, indicated the presence of a shallow ultrabasic intrusion.  相似文献   

4.
Underground gravity observations in deep coal mines using the conventional gravity meters Worden (type Master) and LaCoste-Romberg (model D), both of which have been adapted to the fire damp regulations, can be accurate to ± 10 and ± 3 μgal, respectively. For underground determination of the vertical gradient of gravity the LaCoste-Romberg meter is used together with a specially designed measuring tower. Using this euipment an accuracy in tower gradient observations of ± 30 E was obtained. To apply the equipment to precision gravity observations in underground situations an additional correction, i.e. a gallery correction, is needed. High accuracy in correction is achieved by a new method of three-dimensional modelling. The gravity effect is computed for bodies with a surface approximated by triangular elements, which are generated from corner points of the body. The combination of gallery correction with tower gradient data leads to a new method for in situ density determination. It offers the possibility of horizontal instead of vertical density profiling. To demonstrate the effectiveness of the developments in underground observations the localization of a pump room is presented. Microgravity and tower gradient observations were carried out to detect the cavity. The horizontal gradient was also calculated to give a more reliable location.  相似文献   

5.
Summary A solution of the direct gravity problem for a finite body with variable density is given. The method is based on Green's formula and is applicable when a particular solution of Poisson's equation is known. The attraction due to the body is expressed by integrals over its surface The exact solution of the direct gravity problem, as known from the theory of two-dimensional fields [1–3], is closely connected with the problem of the analytic continuation of the exterior field of the attracting mass system into its interior. In the first place, this is a problem of determining the singularities of the exterior field, their distribution within the system and their nature. This approach to the solution of the direct problem is also meaningful from the point of view of determining the characteristics of the attracting system and, therefore, also of solving the inverse problem. In the case of two-dimensional fields the methods of analytical continuation were widely developed in a series of well-known papers by V. N. Strakhov, and they are mainly based on the methods of the theory of the functions of the complex variable. These methods were also successfully applied by Tsirulskii and Golizdra [1, 2] in treating the homogeneous and inhomogeneous, two-dimensional direct problem by means of Cauchy's integrals. However, as regards three-dimensional fields a number of fundamental problems has not been solved in this respect.Dedicated to 90th Birthday of Professor Frantiek Fiala  相似文献   

6.
7.
Summary A summation method of upward continuation of gravity data has been considered under the assumption that observations are available at regular intervals. The upward continued value has been obtained as the sum of products of individual gravity values and corresponding theoretical coefficients. Besides the usual parameter involving horizontal and vertical distances, the theoretical coefficients have been generalized to be dependent also on i) the order of a low order polynomial assumed to represent the gravity variation around a grid point and ii) the weights assigned to the gravity values at the nearest four grid points used for least-squares determination of the polynomial. Since the observations in practical cases are available over a finite area only, the effect of truncation of the area of the integration has also been discussed separately. The method has been programmed and tested on a three-dimensional model, whose true gravity effects were computed at various levels over a finite area. Upward continued values have been computed under various assumptions about the gravity field in the outside region. Comparisons of these results with the true values indicate that the truncation effect becomes increasingly important for larger values of the ratio of elevation to grid separation and/or when the gravity field is not a local one. It has also been found that the greater is the above ratio, the less important is the effect of weights on the theoretical coefficients and practically vanishes(<10 –4 ) when the ratio is greater than5.0.  相似文献   

8.
Summary As an aid to the interpretation of gravity data for three-dimensional causative bodies, a method of evaluation of the gravity effect of a homogeneous polyhedron at an external point has been worked out. It has been shown that the gravity effect is expressible in terms of some prinitive integrals over the basic triangular faces of the polyhedron. The formula for evaluation of these integrals has also been derived.Contribution from the Earth Physics Branch No. 501.  相似文献   

9.
Repetitive gravity surveys at Pacaya Volcano from 1975 to 1979 revealed time-dependent changes in the gravity field, which although related to volcanic activity, could not be uniquely attributed to elevation changes or mass redistributions because elevation control was lacking. Elevation control was established in July 1979 using precision leveling. Relative elevation and gravity measurements in June and July of 1979, January 1980 and June 1980 indicate concurrent gravity and elevation changes contemporaneous with variations in eruptive activity. From June 1979 to January 1980, while fumarolic activity was dominant, relative to the most remote station, the volcano deflated by at least 195 mm and the gravity field increased by up to 221 μgal. From January 1980 to June 1980, preceding a Strombolian eruption beginning about June 1980, the volcano inflated by at least 19 mm and the gravity field decreased by up to 231 μgal. Gravity change maps for the intervals of January 1978 to June 1979, June 1979 to January 1980, and January 1980 to June 1980 show areas subject to repeated positive and negative gravity change. Some of those areas coincide with areas of maximum elevation change observed in the June 1979–January 1980 and January 1980–June 1980 intervals; however, gravity changes were observed in areas lacking elevation changes. Adjusting observed gravity changes for elevation changes using a free-air value of −3.086 μgal/cm does not substantially alter the pattern, position, or amplitude of the gravity changes. The relationship between gravity changes, elevation changes, and volcanic activity requires a mechanism producing gravity decreases with little inflation during times of increased eruptive activity, and producing gravity increases with subsidence during times of declining eruptive activity. Such a pattern of changes could be produced by a near-surface magma body in which high-density degassed magma is displaced volume for volume by low-density vesiculated magma during time of increased eruptive activity, and in which loss of gasses by fumarolic activity produces a density increase and a reduction in volume of the magma body during periods of declining eruptive activity. Such a pattern of changes could also be induced by a low-density, vesiculated magma body moving upward in the volcanic pile by piecemeal stoping where the high-density rocks of the volcano are replaced on a volume for volume basis by low-density magma during periods of increasing eruptive activity; and by later density increases and magma body volume reductions accompanying devolatilization and devesiculation during periods of declining eruptive activity. Simple density change and density contrast models involving shallow magma bodies at depths of 100 to 200 m indicate density changes or contrasts of about 0.4 g/cm3 could produce the gravity changes.  相似文献   

10.
Summary The effect of an additional homogeneous magnetic field with an intensity of 0–4.5 Oe on the Worden quartz gravity meter No. 961 and on Sharpe quartz gravity meters Nos 173 and 174 was tested. Whereas no effect was observed with the Worden gravity meter, the magnetic field had a measurable effect on both the Sharpe gravity meters. The largest deviation of the reading beam is caused by the horizontal component of the magnetic field which acts in the plane of oscillation of the gravity-meter arm. The Sharpe gravity meter No. 173 is considerably sensitive; a field of 0.2 Oe intensity, corresponding to the magnitude of the horizontal component of the geomagnetic field in mid-latitudes, causes an error in the measurement of gravity of as much as 0.08 mGal. With a view to the different behaviours of the individual quartz gravity meters of the same type in a magnetic field, it should prove expedient to carry out check measurements with all gravity meters and, with regard to the sensitivity of the gravity meter to the magnetic field and the required accuracy of the gravity determination, take into account this perturbing factor in field measurements, as well as laboratory tests of gravity meters.  相似文献   

11.
The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.  相似文献   

12.
In this paper, we present a case study on the use of the normalized source strength (NSS) for interpretation of magnetic and gravity gradient tensors data. This application arises in exploration of nickel, copper and platinum group element (Ni‐Cu‐PGE) deposits in the McFaulds Lake area, Northern Ontario, Canada. In this study, we have used the normalized source strength function derived from recent high resolution aeromagnetic and gravity gradiometry data for locating geological bodies. In our algorithm, we use maxima of the normalized source strength for estimating the horizontal location of the causative body. Then we estimate depth to the source and structural index at that point using the ratio between the normalized source strength and its vertical derivative calculated at two levels; the measurement level and a height h above the measurement level. To discriminate more reliable solutions from spurious ones, we reject solutions with unreasonable estimated structural indices. This method uses an upward continuation filter which reduces the effect of high frequency noise. In the magnetic case, the advantage is that, in general, the normalized magnetic source strength is relatively insensitive to magnetization direction, thus it provides more reliable information than standard techniques when geologic bodies carry remanent magnetization. For dipping gravity sources, the calculated normalized source strength yields a reliable estimate of the source location by peaking right above the top surface. Application of the method on aeromagnetic and gravity gradient tensor data sets from McFaulds Lake area indicates that most of the gravity and magnetic sources are located just beneath a 20 m thick (on average) overburden and delineated magnetic and gravity sources which can be probably approximated by geological contacts and thin dikes, come up to the overburden.  相似文献   

13.
蒋甫玉  高丽坤 《应用地球物理》2012,9(2):119-130,233
针对利用重磁资料增强地质体边界在图像中的视觉效果和提高地质解译的准确性问题,提出应用改进的小子域滤波方法对重力异常及重力梯度张量数据进行增强处理。根据滑动窗口子域平均选择原理,探讨了改进的小子域滤波方法在位场异常数据含有高斯白噪声时,不同窗口大小对地质体边界的识别效果及其在具有不同边界延伸方向地质体中的应用效果。模型试验结果表明,利用改进的小子域滤波对重力梯度张量数据进行增强处理,得到的地质体边界形态失真更小,且受滤波窗口大小、噪声以及地质体边界方向的影响较小;对深部场源体,通过增大滤波窗口的方式,可以较好地反映深部场源体的边界。结合黑龙江省虎林盆地重力异常以及计算的重力梯度张量的处理实例表明改进的小子域滤波法较传统的小子域滤波法增强了对断裂水平位置信息的识别。  相似文献   

14.
The gravimetric parameters of the gravity pole tide are the amplitude factor δ, which is the ratio of gravity variations induced by polar motion for a real Earth to variations computed for a rigid one, and the phase difference κ between the observed and the rigid gravity pole tide. They can be estimated from the records of superconducting gravimeters (SGs). However, they are affected by the loading effect of the ocean pole tide. Recent results from TOPEX/Poseidon (TP) altimeter confirm that the ocean pole tide has a self-consistent equilibrium response. Accordingly, we calculate the gravity loading effects as well as their influence on the gravimetric parameters of gravity pole tide at all the 26 SG stations in the world on the assumption of a self-consistent equilibrium ocean pole tide model. The gravity loading effect is evaluated between 1 January 1997 and 31 December 2006. Numerical results show that the amplitude of the gravity loading effect reaches 10−9 m s−2, which is larger than the accuracy (10−10 m s−2) of a SG. The gravimetric factor δ is 1% larger at all SG stations. Then, the contribution of a self-consistent ocean pole tide to the pole tide gravimetric parameters cannot be ignored as it exceeds the current accuracy of the estimation of the pole tide gravity factors. For the nine stations studied in Ducarme et al. [Ducarme, B., Venedikov, A.P., Arnoso, J., et al., 2006. Global analysis of the GGP superconducting gravimeters network for the estimation of the pole tide gravimetric amplitude factor. J. Geodyn. 41, 334–344.], the mean of the modeled tidal factors δm = 1.1813 agrees very well with the result of a global analysis δCH = 1.1816 ± 0.0047 in that paper. On the other hand, the modeled phase difference κm varies from −0.273° to 0.351°. Comparing to the two main periods of the gravity pole tide, annual period and Chandler period, κm is too small to be considered. Therefore, The computed time difference κL induced by a self-consistent ocean pole tide produces a negligible effect on κm. It confirms the results of Ducarme et al., 2006, where no convincing time difference was found in the SG records.  相似文献   

15.
Abstract

In the context of ageostrophic theory in a homogeneous ocean, a nondimensional number is determined which corresponds to the Ursell number for long gravity waves. It is defined as Q = NL 2/h 3, where N is the amplitude of the wave travelling along the long length-scale direction, L is its length and h (which for gravity waves is the water depth) is given by h=(l 4 f 2/g)1/3. where l is the short length-scale, f the Coriolis parameter and g the acceleration due to gravity. The physical meaning of Q is as follows: if Q? O(1) the free evolution of the wave is linear and weakly dispersive, if Q = O(1) nonlinear and dispersive effects balance out and finally if Q ?O(1) the evolution is nonlinear and non-dispersive. Expressions for the time scales for the development of dispersive and nonlinear effects are also determined. These results apply to topographically trapped waves, namely barotropic continental shelf and double Kelvin waves travelling along a rectilinear topographic variation.  相似文献   

16.
Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, we investigate the possibility of using observations of the glacial isostatic adjustment (GIA) process to constrain the thermal contribution to lateral variations in mantle viscosity. In particular, global historic relative sea level, GPS in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide are used. The lateral viscosity perturbations are inferred from the seismic tomography model S20A by inserting the scaling factor β to determine the contribution of thermal effects versus compositional heterogeneity and non-isotropic pre-stress effects on lateral heterogeneity in mantle viscosity. When β = 1, lateral velocity variations are caused by thermal effects alone. With β < 1, the contribution of thermal effect decreases, so that for β = 0, there is no lateral viscosity variation and the Earth is laterally homogeneous. These lateral viscosity variations are superposed on four different reference models which differ significantly in the lower mantle viscosity. The Coupled Laplace Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic Earth with self-gravitating oceans, induced by the ICE-4G deglaciation model.Results show that the effect of β on uplift rates and gravity rate-of-change is not simple and involves the trade-off between the contribution of lateral viscosity variations in the transition zone and in the lower mantle. Models with small viscosity contrast in the lower mantle cannot explain the observed uplift rates in Laurentide and Fennoscandia. However, the RF3S20 model with a reference viscosity profile simplified from Peltier's VM2 with the value of β around 0.2–0.4 is found to explain most of the global RSL data, the uplift rates in Laurentide and Fennoscandia and the BIFROST horizontal velocity data. In addition, the changes in GIA signals caused by changes in the value of β are large enough to be detected by the data, although uncertainty in other parameters in the GIA models still exists. This may encourage us to further utilize GIA observations to constrain the thermal effect on mantle lateral heterogeneity as geodetic and satellite gravity measurements are improved.  相似文献   

17.
The spectrum of a magnetic or a gravity anomaly due to a body of a given shape with either homogeneous magnetization or uniform density distribution can be expressed as a product of the Fourier transforms of the source geometry and the Green's function. The transform of the source geometry for any irregularly-shaped body can be accurately determined by representing the body as closely as possible by a number of prismatic bodies. The Green's function is not dependent upon the source geometry. So the analytical expression for its transform remains the same for all causative bodies. It is, therefore, not difficult to obtain the spectrum of an anomaly by multiplying the transform of the source geometry by that of the Green's function. Then the inverse of this spectrum, which yields the anomaly in the space domain, is calculated by using the Fast Fourier Transform algorithm. Many examples show the reliability and accuracy of the method for calculating potential field anomalies.  相似文献   

18.
利用北天山地区2016~2019年观测的4期流动重力观测资料,分析研究一年尺度的重力场动态变化特征,并利用小波分析方法,将不同场源深度的重力异常进行分离。通过功率谱分析,获取各阶小波重力细节对应的场源深度。研究结果表明,2017年8月9日精河MS6.6地震前,震中位于负值集中区,四阶小波重力细节显示震中附近出现明显的四象限分布;2020年1月16日库车MS5.6地震前,震中位于负值区,小波重力细节整体量值较小;功率谱估算的场源近似深度与2次地震的震源深度相近。  相似文献   

19.
An underground gravity study was carried out under extreme conditions of the Alpine regions. The lead–zinc mine Bleiberg, Carinthia, was selected as an example to show the possibilities and limitations of the subsurface gravity method. For in situ density determinations, gravity measurements were made in two vertical mine shafts passing through Triassic sedimentary rocks of the Bleiberg Unit. The main prblem in gravity data reduction in extremely rugged topography is the accurate calculation of the terrain effect on underground stations. A general discussion of the various corrections required for the gravity measurements in the mine is presented. The mean interval densities in the two shafts, in limestone, dolomite, and schists formations, were determined as 2.76 and 2.77 g/cm3, respectively, with an accuracy of better than 0.01 g/cm3 for a depth interval of 50 m. The interval densities provide valuable information about the lithological and structural changes in the shaft surroundings and also agree well with the representative hand sample densities. In the second part, the applications of subsurface horizontal gravity surveys in exploration are discussed. Since the influence of topography is less underground because of the greater distance to the surface, subsurface surveys have definite advantages over surface surveys and can be very helpful in locating anomalous density zones in the mines. An example of gravity survey with a station spacing of 10 m at a depth of about 540 m is presented.  相似文献   

20.
Summary The general solution of the direct gravimetric problem for homogeneous and inhomogeneous bodies was derived in[1] in the form of Green-type integrals. In the present paper it is proved that this solution simplifies considerably, if the solution of the external Dirichlet problem, defined by conditions(2), is known for a closed surface bounding the attracting body.Dedicated to RNDr Jan Pícha, CSc., on his 60th Birthday  相似文献   

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

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