排序方式： 共有7条查询结果，搜索用时 218 毫秒
1
1.
Spectral analysis of the full gravity tensor 总被引：9，自引：0，他引：9
2.
Gravity observations made over the Groningen gasfield from 1978 onwards are available. The existing gravity data were reanalysed and a new survey was carried out. The observed gravity changes obtained during four surveys, spanning an 18-year period, were compared with the gravity effect due to natural gas extraction computed from the reservoir model and the production data. The random error in the gravity values is small enough to detect the effect of gas extraction after a few years (up to 2 μgal/year). The trends obtained from observed and calculated gravity changes agree well within their expected error margins after statistical data snooping. Due to the inadequate measurement set-up of the previous surveys, systematic errors present in the gravity data hampered the use of this data for refinement of the reservoir and the production models. However, with improved gravity monitoring, in particular a very systematic network set-up and well-calibrated instruments, this should now be feasible. By using the network analysis from the 1996 measurements for the planning of a network in space and time, most of the current problems can be avoided. 相似文献
3.
When standard boundary element methods (BEM) are used in order to solve the linearized vector Molodensky problem we are confronted with
two problems: (1) the absence of O(|x|−2) terms in the decay condition is not taken into account, since the single-layer ansatz, which is commonly used as representation
of the disturbing potential, is of the order O(|x|−1) as x→∞. This implies that the standard theory of Galerkin BEM is not applicable since the injectivity of the integral operator
fails; (2) the N×N stiffness matrix is dense, with N typically of the order 105. Without fast algorithms, which provide suitable approximations to the stiffness matrix by a sparse one with O(N(logN)
s
), s≥0, non-zero elements, high-resolution global gravity field recovery is not feasible. Solutions to both problems are proposed.
(1) A proper variational formulation taking the decay condition into account is based on some closed subspace of co-dimension
3 of the space of square integrable functions on the boundary surface. Instead of imposing the constraints directly on the
boundary element trial space, they are incorporated into a variational formulation by penalization with a Lagrange multiplier.
The conforming discretization yields an augmented linear system of equations of dimension N+3×N+3. The penalty term guarantees the well-posedness of the problem, and gives precise information about the incompatibility
of the data. (2) Since the upper left submatrix of dimension N×N of the augmented system is the stiffness matrix of the standard BEM, the approach allows all techniques to be used to generate
sparse approximations to the stiffness matrix, such as wavelets, fast multipole methods, panel clustering etc., without any
modification. A combination of panel clustering and fast multipole method is used in order to solve the augmented linear system
of equations in O(N) operations. The method is based on an approximation of the kernel function of the integral operator by a degenerate kernel
in the far field, which is provided by a multipole expansion of the kernel function. Numerical experiments show that the fast
algorithm is superior to the standard BEM algorithm in terms of CPU time by about three orders of magnitude for N=65 538 unknowns. Similar holds for the storage requirements. About 30 iterations are necessary in order to solve the linear
system of equations using the generalized minimum residual method (GMRES). The number of iterations is almost independent
of the number of unknowns, which indicates good conditioning of the system matrix.
Received: 16 October 1999 / Accepted: 28 February 2001 相似文献
4.
5.
A general scheme is given for the solution in a least-squares sense of the geodetic boundary value problem in a spherical,
constant-radius approximation, both uniquely and overdetermined, for a large class of observations. The only conditions are
that the relation of the observations to the disturbing potential is such that a diagonalization in the spectrum can be found
and that the error-covariance function of the observations is isotropic and homogeneous. Most types of observations used in
physical geodesy can be adjusted to fit into this approach. Examples are gravity anomalies, deflections of the vertical and
the second derivatives of the gravity potential.
Received: 3 November 1999 / Accepted: 25 September 2000 相似文献
6.
Least-squares by observation equations is applied to the solution of geodetic boundary value problems (g.b.v.p.). The procedure is explained solving the vectorial Stokes problem in spherical and constant radius approximation. The results
are Stokes and Vening-Meinesz integrals and, in addition, the respective a posteriori variance-covariances.
Employing the same procedure the overdeterminedg.b.v.p. has been solved for observable functions potential, scalar gravity, astronomical latitude and longitude, gravity gradients
Гxz, Гyz, and Гzz and three-dimensional geocentric positions. The solutions of a large variety of uniquely and overdeterminedg.b.v.p.'s can be obtained from it by specializing weights. Interesting is that the anomalous potential can be determined—up to a
constant—from astronomical latitude and longitude in combination with either {Гxz,Гyz} or horizontal coordinate corrections Δx and Δy, or both. Dual to the formulation in terms of observation equations the overdeterminedg.b.v.p.'s can as well be solved by condition equations.
Constant radius approximation can be overcome in an iterative approach. For the Stokes problem this results in the solution
of the “simple” Molodenskii problem. Finally defining an error covariance model with a Krarup-type kernel first results were
obtained for a posteriori variance-covariance and reliability analysis. 相似文献
7.
The use of degree variances in satellite gradiometry 总被引：1，自引：0，他引：1
1