共查询到20条相似文献，搜索用时 62 毫秒
1.
C. C. Tscherning 《Geophysical Journal International》1979,59(1):147-153
Summary. The anomalous (gravitational) potential of the Earth, T , is split in two parts, T= T C + T M . Here T M is a harmonic function generated by known mass density anomalies and T C =T-T M . This function will also be a harmonic function, which therefore may be approximated using the method of collocation, based on known gravity anomalies or altimeter derived geoid undulations, for example. Gravity anomalies can then be predicted using the known linearized relationship between T and Δ g . This procedure may give a 40–50 per cent increase in the precision of the prediction results as compared to a procedure where mass density anomalies are not taken into account. 相似文献
2.
一种基于超曲面样条函数的三维空间插值方法 总被引：1，自引：0，他引：1
将二维曲面样条函数插值法进一步引申,提出一种基于超曲面样条函数进行三维空间插值的新方法,介绍超曲面样条函数的构造方法,给出利用该方法进行三维空间插值的数值实现过程。最后以插值形成三维空间规则数据场为例,证明利用超曲面样条函数法进行三维空间插值的可行性和准确性。 相似文献
3.
Tim van Zon Pavel Ditmar Kabir Roy Chowdhury Jaap Mondt 《Geophysical Journal International》2007,170(2):481-491
In this paper we present a method which allows delineation of geologic structures in a bi-modal lithotype setting. We propose to use gravity data in combination with a priori information about the density contrast between the two lithotypes. The iterative method uses an objective function with five tunable parameters which need to be set. Using an efficient parameter search, suitable ranges of these are investigated to determine their optimal values, respectively, which in turn, ensures good inversion results.
The approach produces structural images of the subsurface, without the need of an a priori density model; the depth to the top of the inhomogeneity is also retrieved.
Besides synthetic simulations, the methodology has also been applied to a small gravity data set, acquired by the industry over a basinal structure. A consistent, bi-modal image of the bedrock depression is obtained from the data, which, in this case, was the goal. Other potential areas of application include delineation of salt structures and ore deposits. 相似文献
The approach produces structural images of the subsurface, without the need of an a priori density model; the depth to the top of the inhomogeneity is also retrieved.
Besides synthetic simulations, the methodology has also been applied to a small gravity data set, acquired by the industry over a basinal structure. A consistent, bi-modal image of the bedrock depression is obtained from the data, which, in this case, was the goal. Other potential areas of application include delineation of salt structures and ore deposits. 相似文献
4.
Rosaria Tondi Roberto de Franco & Riccardo Barzaghi' 《Geophysical Journal International》2000,141(3):679-698
A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area. 相似文献
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area. 相似文献
5.
S. Takasaki N. Sato A. Kadokura H. Yamagishi H. Kawano Y. Ebihara Y.-M. Tanaka 《Polar Science》2008,2(2):73-86
In the austral summer of 2006–2007, the 48th Japanese Antarctic Research Expedition (JARE-48) installed two unmanned low-power magnetometers to form a closely spaced magnetometer network in combination with the permanent sites at Japan's Syowa Station in Antarctica. To identify field line resonances (FLRs), gradient methods are applied to the data from three adjacent sites in Antarctica and data from conjugate points in Antarctica and Iceland. By analyzing the data from the Antarctic and Icelandic sites individually, the structure of FLRs with high coherence is clearly identified. However, by analyzing the data from closely spaced Antarctic sites, it is more difficult to identify the signature of FLRs because of the inclusion of multiple signals related to the local geomagnetic pulsations over a broad frequency range. The frequency and resonance width of FLRs are determined by applying the amplitude phase gradient method (APGM) to the data from Antarctic sites. This yields the eigenfrequency as a continuous function of ground latitudes in the area surrounding Syowa Station. The mass density in the equatorial region at the L of the auroral zones is estimated from the obtained FLR frequency by numerically solving the standing Alfvén wave equation. The mass density thus obtained is consistent with observational results from previous in situ measurements by spacecraft. The results of the present study demonstrate that data from geomagnetic conjugate points are helpful in identifying FLR in cases in which the magnetometers are too close to each other to enable identification. Once FLR is identified, APGM can be applied to the identified FLR, yielding the FLR frequency as a continuous function of ground latitudes. Therefore, the magnetospheric equatorial mass density is readily estimated with high spatial resolution. 相似文献
6.
M. A. Brodsky 《Geophysical Journal International》1981,66(3):727-732
Summary. The inverse gravity potential problem consists in the determination of the form and the density of the body by its exterior gravity potential. We describe two similar classes of bodies for which this problem has a unique constructive solution.
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ1 ( z )ρ2 ( R , ø) density distribution, where z is the cylindrical coordinate; R , ø are the polar coordinates in a section plane. This class is important for prospecting geophysics in that it allows us to determine in a unique and constructive way, the function ρ1 ( R , ø), the length, form and orientation of the cylinder if we know the function ρ1 ( z ) and the exterior potential. The classical moment problem of functions is the basis for the solution of this problem.
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ1 (τ)ρ2 (θ, ø) where τ, θ, ø are the spherical coordinates. The function ρ1 (θ, ø), length and form of spherical sections can be uniquely determined by exterior potential if we know the function ρ1 (τ). We propose a new constructive method for harmonic continuation of the gravity potential into the region containing the perturbing masses for the solution of the problem. 相似文献
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ
7.
The correspondence between gravitational attraction and surface displacement due to volume expansion
Summary. It is shown that the formulas for the calculation of surface displacement due to a volume source in a half-space may be derived from a potential of the same form as the gravitational potential. The similarity between these expressions allows the use of methods developed for gravity in the study of surface displacement. As an example of the applicability of the correspondence, an efficient algorithm is developed for the calculation of uplift due to volume expansion within a homogeneous half-space. The method may be generalized to three-dimensional volumes undergoing expansion in a half-space. 相似文献
8.
张力样条是在三次样条函数中加入张力系数而构成。根据相同的粒度数据用两种方法在计算机上所作的粒度分布曲线图形的对比,显示出张力样条方法比三次样条方法有明显的优越性,它消除了插值曲线的不必要的摆动,在保持曲线平滑的前提下,能更准确地表达沉积物的粒度特征。用张力样条结合将离散点连成光滑曲线的阿克玛方法进行插值,使复杂的粒度分布曲线的插值同时保持准确性和平滑性,在沉积物粒度的数据的计算机处理方面有明显的意义。 相似文献
9.
Choice of norm for the density distribution of the Earth 总被引：1，自引：0，他引：1
Summary. The determination of the density distribution of the Earth from gravity data is called the inverse gravimetric problem. A unique solution to this problem may be obtained by introducing a priori data concerning the covariance of density anomalies. This is equivalent to requiring the density to fulfil a minimum norm condition. The generally used norm is the one equal to the integral of the square of the density distribution ( L2 -norm), the use of which implies that blocks of constant density are uncorrelated. It is shown that for harmonic anomalous density distributions this leads to an external gravity field with a power spectrum (degree-variances) which tends too slowly to zero, i.e. implying gravity anomalies much less correlated than actually observed. It is proposed to use a stronger norm, equal to the integral of the square sum of the derivatives of the density distribution. As a consequence of this, base functions which are constant within blocks, are no longer a natural choice when solving the inverse gravimetric problem. Instead a block with a linearly varying density may be used. A formula for the potential of such a block is derived. 相似文献
10.
地球重力场的变化是导致陆地水储量变化的重要因素之一,利用GRACE(Gravity Recovery and Climate Experiment)重力场恢复与气候实验重力卫星数据,结合GLDAS(Global Land Data Assimilation Systems)全球陆面数据同化系统和实测地下水位数据,反演和田地区克里雅河流域11年间四季和田地区的陆地水储量动态变化,模拟计算地下水等效水高变化趋势,构建了地下水水位估算模型。研究结果表明：和田地区春、夏两季的陆地水储量呈现出增加趋势,而秋、冬两季出现亏损状态;GRACE地球重力卫星所反演的陆地水储量比GLDAS同化系统所模拟的水资源变化更为剧烈,但2类数据的动态变化拟合度很高;GLDAS水资源等效水高二阶微分、GLDAS水资源变化倒数一阶微分、GRACE陆地水储量变化倒数变化、地下水储量变化一阶微分的敏感程度最高,构建的多元逐步回归模型明显优于线性函数,且水位深度越浅,该估算模型的适用性越高。 相似文献
11.
Elastic finite element models are applied to investigate the effects of topography and medium heterogeneities on the surface deformation and the gravity field produced by volcanic pressure sources. Changes in the gravity field cannot be interpreted only in terms of gain of mass disregarding the ground deformation of the rocks surrounding the source. Contributions to gravity changes depend also on surface and subsurface mass redistribution driven by dilation of the volcanic source. Both ground deformation and gravity changes were firstly evaluated by solving a coupled axisymmetric problem to estimate the effects of topography and medium heterogeneities. Numerical results show significant discrepancies in the ground deformation and gravity field compared to those predicted by analytical solutions, which disregard topography, elastic heterogeneities and density subsurface structures. With this in mind, we reviewed the expected gravity changes accompanying the 1993–1997 inflation phase on Mt Etna by setting up a fully 3-D finite element model in which we used the real topography, to include the geometry, and seismic tomography, to infer the crustal heterogeneities. The inflation phase was clearly detected by different geodetic techniques (EDM, GPS, SAR and levelling data) that showed a uniform expansion of the overall volcano edifice. When the gravity data are integrated with ground deformation data and a coupled FEM modelling was solved, a mass intrusion could have occurred at depth to justify both ground deformation and gravity observations. 相似文献
12.
This paper summarizes principles and measures of comprehensive division of mountainous areas, as well as clarifies meaning, structure, function and path established for the map and file information visibility system (MFIVS). And then, taking Huaihua City of Hunan Province as an example, and based on the MFIVS means, concrete attempts on regionalization are carried out. The result is relatively objective and accurate, which illuminates that the method is a comprehensive one with the characteristics of concision, applicability and effectiveness. 相似文献
13.
A new method for comprehensive division of mountainous area: taking Huaihua City as an example 总被引：2，自引：0，他引：2
1 IntroductionAs human society has entered the 21st century, knowledge-based economy and "Digital Earth" project are paid much attention to. Studies on regional sustainable development come into the manipulating stage. Thus, traditional division needs to be reconsidered[1]. In new era, division work should be both in conformity with rationality of eco-environment, and be beneficial to socio-economic development, i.e., eco-economic comprehensive division should be taken into account[2,3]. As … 相似文献
14.
Abel Amirbekyan Volker Michel Frederik J. Simons † 《Geophysical Journal International》2008,174(2):617-628
We present a mathematical framework and a new methodology for the parametrization of surface wave phase-speed models, based on traveltime data. Our method is neither purely local, like block-based approaches, nor is it purely global, like those based on spherical harmonic basis functions. Rather, it combines the well-known theory and practical utility of the spherical harmonics with the spatial localization properties of spline basis functions. We derive the theoretical foundations for the application of harmonic spherical splines to surface wave tomography and summarize the results of numerous numerical tests illustrating the performance of a practical inversion scheme based upon them. Our presentation is based on the notion of reproducing-kernel Hilbert spaces, which lends itself to the parametrization of fully 3-D tomographic earth models that include body waves as well. 相似文献
15.
当前气温插值方法多将高程、经纬度等作为影响因素,无法解决风向、风速对气温空间分布的影响问题。该文提出一种顾及风向和风速的气温空间插值方法:1)根据气象观测站点的风向、风速数据进行插值,生成连续的风场表面;2)基于该风场数据利用高斯扩散模型构建成本表面;3)利用Dijkstra最短路径算法计算观测点与待求点的累计移动成本最短路径(SPOCMC);4)将SPOCMC作为协变量进行薄盘光滑样条插值以实现气温插值。为验证该方法的有效性和可靠性,选取山东省109个气象观测点样本数据,以SPOCMC、DEM和SPOCMC-DEM分别作为协变量对气温数据进行薄盘光滑样条插值,结果显示:SPOCMC-DEM法的MAE和RMSE均值(分别为0.517、0.779)略低于SPOCMC法(0.583、1.016),显著低于DEM法(0.809、1.231),表明添加SPOCMC作为协变量可有效提高气温空间插值的准确性。在江苏省与贵州省的插值实验结果证明了该方法的普适性。 相似文献
16.
A self-adaptive automated parametrization approach is suggested for the sequential inversion of controlled-source seismic tomography and gravity data. The velocities and interfaces are parametrized by their Haar wavelet expansion coefficients. Only those coefficients that are well constrained by the data, as measured by the number of rays that cross the corresponding wavelet function support area and their angular coverage, are inverted for, others are set to zero. This approach results in a reasonable distribution of resolution throughout the model even in cases of irregular ray coverage and does overcome the trade-off between different types of model parameters. A modified sequential inversion approach is suggested to join the traveltimes and gravity anomalies inversion. An algorithm is developed that inverts for smooth velocity and density variations inside the seismic layer, the position of its bottom interface as well as for optimal values of the velocity-to-density regression coefficients. The algorithm makes use of direct (diving), reflected and head (critically refracted) wave traveltimes. The algorithm workflow is demonstrated on a synthetic data example. 相似文献
17.
Error covariance estimates are necessary information for the combination of solutions resulting from different kinds of data or methods, or for the assimilation of new results in already existing solutions. Such a combination or assimilation process demands proper weighting of the data, in order for the combination to be optimal and the error estimates of the results realistic. One flexible method for the gravity field approximation is least-squares collocation leading to optimal solutions for the predicted quantities and their error covariance estimates. The drawback of this method is related to the current ability of computers in handling very large systems of linear equations produced by an equally large amount of available input data. This problem becomes more serious when error covariance estimates have to be simultaneously computed. Using numerical experiments aiming at revealing dependencies between error covariance estimates and given features of the input data we investigate the possibility of a straightforward estimation of error covariance functions exploiting known characteristics of the observations. The experiments using gravity anomalies for the computation of geoid heights and the associated error covariance functions were conducted in the Arctic region north of 64° latitude. The correlation between the known features of the data and the parameters variance and correlation length of the computed error covariance functions was estimated using multiple regression analysis. The results showed that a satisfactory a priori estimation of these parameters was not possible, at least in the area considered. 相似文献
18.
Interpolation using a generalized Green's function for a spherical surface spline in tension 总被引：1，自引：0，他引：1
A variety of methods exist for interpolating Cartesian or spherical surface data onto an equidistant lattice in a procedure known as gridding. Methods based on Green's functions are particularly simple to implement. In such methods, the Green's function for the gridding operator is determined and the resulting gridding solution is composed of the superposition of contributions from each data constraint, weighted by the Green's function evaluated for all output–input point separations. The Green's function method allows for considerable flexibility, such as complete freedom in specifying where the solution will be evaluated (it does not have to be on a lattice) and the ability to include both surface heights and surface gradients as data constraints. Green's function solutions for Cartesian data in 1-, 2- and 3-D spaces are well known, as is the dilogarithm solution for minimum curvature spline on a spherical surface. Here, the spherical surface case is extended to include tension and the new generalized Green's function is derived. It is shown that the new function reduces to the dilogarithm solution in the limit of zero tension. Properties of the new function are examined and the new gridding method is implemented in Matlab® and demonstrated on three geophysical data sets. 相似文献
19.
J.M. Abril 《Journal of Paleolimnology》2003,30(4):363-370
Classical treatment of mass conservation of solids in growing sediments states an advection-diffusion equation for the bulk sediment density. Thus, the diffusion coefficient has to account for elemental processes of exchange of solid particles by pore water or reciprocally. Nevertheless, in a gravity field, these exchanges are forced and cannot be treated as diffusion but as mass flow. A compaction potential energy is defined so that its spatial gradients force a mass flow involving a conductivity function. This leads to a more consistent definition of mass sedimentation rates and to a writing of the continuity equation for density involving only an advection term. Typical bulk density profiles show an asymptotic increase with depth. With the present formulation, this can be obtained as a steady-state solution under constant sedimentation rate and constant conductivity, while the classical formulation fails to do it. Alternatively, it can be found that under these conditions, the compaction potential is a linear function of the bulk density. The mass flow due to compaction and the compaction potential are found for several sediment cores from literature data. From this basis, the advection-diffusion equation for a particle-associated tracer is rewritten. In particular, when the mass flow term due to compaction is considered, the resulting sub-grid scale processes lead to a different formulation of the diffusive fluxes: they are proportional to the gradients of specific concentration in solids instead of the concentration per unit bulk volume. This new formulation is most suitable to find out analytical solutions for radiometric dating models involving mixing and compaction. Numerical solutions are found for the new and the classical treatment for some particular cases to illustrate differences. 相似文献