共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
N. I. Berezina V. I. Dmitriev N. A. Merschikova 《Izvestiya Physics of the Solid Earth》2013,49(3):350-355
We consider the iterative numerical method for solving two-dimensional (2D) inverse problems of magnetotelluric sounding, which significantly reduces the computational burden of the inverse problem solution in the class of quasi-layered models. The idea of the method is to replace the operator of the direct 2D problem of calculating the low-frequency electromagnetic field in a quasi-layered medium by a quasi-one dimensional operator at each observation point. The method is applicable for solving the inverse problems of magnetotellurics with either the E- and H-polarized fields and in the case when the inverse problem is simultaneously solved using the impedance values for the fields with both polarizations. We describe the numerical method and present the examples of its application to the numerical solution of a number of model inverse problems of magnetotelluric sounding. 相似文献
3.
Determination of spherical harmonic coefficients of the Earth’s gravity field is often an ill-posed problem and leads to solving an ill-conditioned system of equations. Inversion of such a system is critical, as small errors of data will yield large variations in the result. Regularization is a method to solve such an unstable system of equations. In this study, direct methods of Tikhonov, truncated and damped singular value decomposition and iterative methods of ν, algebraic reconstruction technique, range restricted generalized minimum residual and conjugate gradient are used to solve the normal equations constructed based on range rate data of the gravity field and climate experiment (GRACE) for specific periods. Numerical studies show that the Tikhonov regularization and damped singular value decomposition methods for which the regularization parameter is estimated using quasioptimal criterion deliver the smoothest solutions. Each regularized solution is compared to the global land data assimilation system (GLDAS) hydrological model. The Tikhonov regularization with L-curve delivers a solution with high correlation with this model and a relatively small standard deviation over oceans. Among iterative methods, conjugate gradient is the most suited one for the same reasons and it has the shortest computation time. 相似文献
4.
Linearized inversion methods such as Gauss‐Newton and multiple re‐weighted least‐squares are iterative processes in which an update in the current model is computed as a function of data misfit and the gradient of data with respect to model parameters. The main advantage of those methods is their ability to refine the model parameters although they have a high computational cost for seismic inversion. In the Gauss‐Newton method a system of equations, corresponding to the sensitivity matrix, is solved in the least‐squares sense at each iteration, while in the multiple re‐weighted least‐squares method many systems are solved using the same sensitivity matrix. The sensitivity matrix arising from these methods is usually not sparse, thus limiting the use of standard preconditioners in the solution of the linearized systems. For reduction of the computational cost of the linearized inversion methods, we propose the use of preconditioners based on a partial orthogonalization of the columns of the sensitivity matrix. The new approach collapses a band of co‐diagonals of the normal equations matrix into the main diagonal, being equivalent to computing the least‐squares solution starting from a partial solution of the linear system. The preconditioning is driven by a bandwidth L which can be interpreted as the distance for which the correlation between model parameters is relevant. To illustrate the benefit of the proposed approach to the reduction of the computational cost of the inversion we apply the multiple re‐weighted least‐squares method to the 2D acoustic seismic waveform inversion problem. We verify the reduction in the number of iterations in the conjugate'gradient algorithm as the bandwidth of the preconditioners increases. This effect reduces the total computational cost of inversion as well. 相似文献
5.
Nonlinear seismic response analysis of earth dams 总被引：1，自引：0，他引：1
The objective of this paper is to propose a general and efficient numerical procedure for analysing the dynamic response of geotechnical structures, which are considered as both nonlinear and two phase systems. In Section 2, the appropriate coupled dynamic field equations for the response of a two-phase soil system are briefly reviewed. The finite element spatial discretization of the field equations is described and time integration for the resulting nonlinear semi-discrete finite element equations is discussed. In Section 3, iterative techniques are examined for the solution of the global nonlinear system of finite element equations. A large amount of computational effort is expended in the iterative phase of the solution and so the iterative procedure used must be both reliable and efficient. The performance of three iterative procedure is examined: Newton Raphson, Modified Newton Raphson and Quasi-Newton methods, including BGFS and Broyden updates. Finally, in Section 4, the elasto-plastic earthquake response analysis of a two phase nonhomogeneous earth dam is presented. Extensive documentation exists1 for the particular problem selected including recorded earthquake motions at the base and crest of the dam. The results of the numerical calculations are compared to the recorded response of the dam. 相似文献
6.
David E. Dougherty 《Advances in water resources》1985,8(4):201-209
An unconditionally stable explicit time integrator has recently been developed for parabolic systems of equations. This rational Runge Kutta (RRK) method, proposed by Wambecq1 and Hairer2, has been applied by Liu et al.3 to linear heat conduction problems in a time-partitioned solution context. An important practical question is whether the method has application for the solution of (nearly) hyperbolic equations as well.In this paper the RRK method is applied to a nonlinear heat conduction problem, the advection-diffusion equation, and the hyperbolic Buckley-Leverett problem. The method is, indeed, found to be unconditionally stable for the linear heat conduction problem and performs satisfactorily for the nonlinear heat flow case. A heuristic limitation on the utility of RRK for the advection-diffusion equation arises in the Courant number; for the second-order accurate one-step two-stage RRK method, a limiting Courant number of 2 applies. First order upwinding is not as effective when used with RRK as with Euler one-step methods. The method is found to perform poorly for the Buckley-Leverett problem. 相似文献
7.
R. W. James 《Pure and Applied Geophysics》1967,68(1):83-89
Summary General recurrence relations between the coefficients in thenth and (n+1)th order spherical harmonic multipole expansions are derived. The particular application presented here is the derivation of the equations concerned with representing the geomagnetic field by magnetic multipoles. The equations up to the 3rd order multipole are given as an example of the method. The main advantage in using these recurrence relations rather than other methods is that the mathematics is reduced to merely a matter of successive substitutions and this allows a fast step by step generation of the required equations, in a form for which there is a simple numerical program for solution. 相似文献
8.
Complex aquifer systems are often modeled with quasi-three-dimensional models, which consider two-dimensional horizontal flow in the aquifers and one-dimensional vertical flow through aquitards. When the aquifer system consists of a phreatic aquifer and one or more semiconfined aquifers connected by aquitards, the discrete model consists of a nonlinear system of algebraic equations, because the transmissivity of the phreatic aquifer depends on the phreatic head. If the water extraction is very high, the phreatic aquifer can be depleted and the equations of the model must be modified accordingly. There are not simple and general criteria to state if the phreatic aquifer is depleted before solving the system of equations. Therefore, the iterative procedures (e.g., relaxation methods), used to find the solution to the forward problem, must handle these particular conditions and can suffer several problems of convergence. These problems can be caused by the choice of the initial head values or of the relaxation coefficient of the iterative algorithms; however, they can also be caused by the nonexistence or nonuniqueness of the solution to the system of nonlinear equations. The study of existence and uniqueness of the general problem is very difficult and, therefore, we consider a simplified problem, for which the discrete model can be handled analytically. The results of the numerical experiments show that the solution to the forward problem can be nonunique. Only for some cases it is possible to invoke physical arguments to eliminate tentative solutions. 相似文献
9.
I. V. Yegorov 《Izvestiya Physics of the Solid Earth》2009,45(9):812-821
The algorithm for numerically solving the direct 3-D problem of calculating the electromagnetic field varying harmonically
in an arbitrary inhomogeneous 3-D media is developed on the base of the Trefftz method, which has not been used previously
in geoelectrics. The corresponding system of algebraic equations has been solved with the use of a modification of the well-known
Kaczmarz iterative method. A cyclic method of equalization is used as a procedure of preconditioning the matrix of the system. 相似文献
10.
利用伴随算子L*,直接的偏移方法通常导致一个低分辨率或模糊的地震成像.线性化偏移反演方法需求解一个最小二乘问题.但直接的最小二乘方法的数值不稳定,为目视解译带来困难.本文建立约束正则化数学模型,研究了地震偏移反演成像问题的迭代正则化求解方法.首先对最小二乘问题施加正则化约束,接着利用梯度迭代法求解反演成像问题,特别是提出了共轭梯度方法的混合实现技巧.为了表征该方法的可实际利用性,分别对一维,二维和三维地震模型进行了数值模拟.结果表明该正则偏移反演成像方法是有效的,对于实际的地震成像问题有着良好的应用前景. 相似文献
11.
在三维频率域电磁法的正演模拟方法中,有限元方法具有计算精度高、适应性强的优点,近年来来得到了越来越多的关注.在正演过程中,主要的计算量集中在求解由偏微分方程组离散得到的线性方程组上,因此求解线性方程组关系着正演计算速度以及模拟精度.由于由有限元方法离散得到的复系数线性方程组条件数非常大,使用常规的迭代法和预条件很难收敛.目前大多数的研究工作采用直接解法,需要大量的计算机内存,限制了可求解问题的规模.本文研究了线性方程组的迭代解法,通过将复系数线性方程组转化为其实对称形式,构造分块对角预条件.在应用预条件的过程中,需要求解两个较小的实数方程,通过辅助空间解法求解.本文的算法适用于可控源电磁法和大地电磁法,对一系列的数值算例的模拟结果证明了迭代算法的效率,结果表明迭代算法可以在小于20次迭代内收敛,同时迭代次数与模型电阻率、问题规模和频率无关. 相似文献
12.
Numerical Studies on the Harmonic Downward Continuation of Band-Limited Airborne Gravity 总被引：1，自引：0，他引：1
In this paper, the numerical stability and efficiency of methods of harmonic downward continuation from flying altitudes are treated for sampled gravity field data. The problem is first formulated in its continuous form, i.e. as the inverse solution of the spherical Dirichlet problem, and is then approximated by Gaussian quadrature to yield a finite system of linear equations. The numerical stability of this system is investigated for both error-free gravity data and for the noisy and band-limited gravity measurements usually obtained from airborne gravity surveys. It can be shown that the system becomes ill-conditioned, once the ratio between flying altitude and data sampling rate exceeds a certain limit. It can also be shown that noisy measurements tend to generate a solution that is practically useless, long before the system becomes ill-conditioned. Therefore, instead of treating the general solution of the discrete downward continuation problem, the more modest question is studied, for which range of flying altitudes and sampling rates, the numerical solution of the discrete linear system can be considered as practically useful. Practically useful will be defined heuristically as of sufficient accuracy and stability to satisfy the requirements of the user. The question will be investigated for the specific application of geoid computation from gravity data sampled at flying altitudes. In this case, a stable solution with a standard deviation of a few centimeters is required. Typical flight parameters are heights of 2–6 km, a minimum half-wavelength resolution of 2 km, and data noise between 0.5 and 1.5 mGal. Different methods of geoid determination, different solution techniques for the resulting systems of linear equations, and different minimization principles will be compared. As a result operational parameters will be defined which, for a given noise level, will result in a geoid accuracy of a few centimeters for the estimated band-limited gravity field spectrum. 相似文献
13.
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 相似文献
14.
Abstract In dealing with the transient sediment transport problem, the commonly used uncoupled model may not be suitable. The uncoupling technique is intended to separate the physical coupling phenomenon of water flow and sediment transport into two independent processes. Very often, as a result, severe numerical oscillation and solution instability problems appear in the simulation of transient sediment transport in alluvial channels. The coupled model, which simultaneously solves water flow continuity, momentum and sediment continuity equations, gives fewer numerical oscillation and solution instability problems. In this article, a coupled model using a matrix double-sweep method to solve the system of nonlinear algebraic equations has been developed. Several test runs designed on the basis of a schematic model have been performed. The numerical oscillation and solution instability problems have been investigated through a comparison with those obtained from an uncoupled model. Based on the proposed case studies, it can be concluded that, for transient bed evolution, the performance of the coupled model is much better than that of the uncoupled model. The numerical oscillation is reduced and the solution is more stable. This newly developed coupled model was also applied to the Cho-Shui River in Taiwan. This application study implied that the effect of the peaky flood wave propagation on the bed evolution could be simulated better by the coupled model than by the uncoupled model. 相似文献
15.
《Advances in water resources》1998,21(5):327-337
A p finite element scheme and parallel iterative solver are introduced for a modified form of the shallow water equations. The governing equations are the three-dimensional shallow water equations. After a harmonic decomposition in time and rearrangement, the resulting equations are a complex Helmholz problem for surface elevation, and a complex momentum equation for the horizontal velocity. Both equations are nonlinear and the resulting system is solved using the Picard iteration combined with a preconditioned biconjugate gradient (PBCG) method for the linearized subproblems. A subdomain-based parallel preconditioner is developed which uses incomplete LU factorization with thresholding (ILUT) methods within subdomains, overlapping ILUT factorizations for subdomain boundaries and under-relaxed iteration for the resulting block system. The method builds on techniques successfully applied to linear elements by introducing ordering and condensation techniques to handle uniform p refinement. The combined methods show good performance for a range of p (element order), h (element size), and N (number of processors). Performance and scalability results are presented for a field scale problem where up to 512 processors are used. 相似文献
16.
Xiwen Wang 《地震学报(英文版)》1992,5(1):133-143
In this paper, the application of Backus—Gilbert’s inversion method to the potential field anomalies for evaluating gravity
and magnetic inversion solutions is discussed. Errors in data and singularity of kernels in the equations result in difficulties
in solving equations. The application of regularization method similar to spectral expansion method makes calculation fast
and easily. To make solution stable, constraints are used, which make the spread of solutions become narrow, standard deviation
become small and iterative computations of inversion become fast. Finally, the author analyses specifically two profiles of
Yunchen basin and calculated the Moho interface and the Curie isotherm of these two profiles.
The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 212–221, 1991. 相似文献
17.
Summary This paper is a continuation of[1]. It is mainly devoted to problems connected with the application of the method of determination of geometrical spreading in laterally inhomogeneous media with curved interfaces based on the solution of eight (in a three-dimensional medium) or two (in a two-dimensional medium) linear ordinary differential equations of the first order. The method of determination of the partial derivatives of velocity with respect to the special coordinates, connected with the ray under investigation, and the methods of determination of the initial values for the system of differential equations at the source and at the interfaces are proposed. 相似文献
18.
Summary Green's theorem on harmonic functions makes it possible to determine the integral relationship between the harmonic function and its derivative with respect to the normal on a closed Lyapunov surface. The conditions of solvability are given by Fredholm's theory of integral equations. The solution for a sphere was presented by Molodenskii[3] and the general solution with the help of Molodenskii's parameter k by Ostach[4]. The present paper indicates a possibility of solving this problem with the help of a system of linear algebraic equations, a simplified modification of the Ostach-Molodenskii solution and, finally, a method, based on Eremeev's solution of the fundamental integral equation[5]. 相似文献
19.
Haitjema HM 《Ground water》2006,44(1):102-105
The analytic element method, like the boundary integral equation method, gives rise to a system of equations with a fully populated coefficient matrix. For simple problems, these systems of equations are linear, and a direct solution method, such as Gauss elimination, offers the most efficient solution strategy. However, more realistic models of regional ground water flow involve nonlinear equations, particularly when including surface water and ground water interactions. The problem may still be solved by use of Gauss elimination, but it requires an iterative procedure with a reconstruction and decomposition of the coefficient matrix at every iteration step. The nonlinearities manifest themselves as changes in individual matrix coefficients and the elimination (or reintroduction) of several equations between one iteration and the other. The repeated matrix reconstruction and decomposition is computationally intense and may be avoided by use of the Sherman-Morrison formula, which can be used to modify the original solution in accordance with (small) changes in the coefficient matrix. The computational efficiency of the Sherman-Morrison formula decreases with increasing numbers of equations to be modified. In view of this, the Sherman-Morrison formula is only used to remove equations from the original set of equations, while treating all other nonlinearities by use of an iterative refinement procedure. 相似文献
20.
Computational seismic modelling (CSM) plays an important role in the geophysical industry as an established aid to seismic interpreters. Numerical solution of the elastic wave equations has proved to be a very important tool for geophysicists in both forward modelling and migration. Among the techniques generally used in CSM, we consider the finite-element method (FEM) and investigate its computational and visualization requirements. The CSMFEM program, designed for this purpose and developed on an IBM 3090 computer with vector facility, is described in detail. It constitutes a numerical laboratory for performing computer experiments. Two Newmark type algorithms for time integration are compared with other time integration schemes, and both direct and iterative methods for solving the corresponding large sparse system of linear algebraic equations are analysed. Several numerical experiments to simulate seismic energy propagation through heterogeneous media are performed. Synthetics in the form of common shot gathers, vertical seismic profiles and snapshots are suitably displayed, since with the large amounts of data obtained from CSM research, methods for visualization of the computed results must be developed. The FEM is compared with other numerical tools, such as finite-difference and pseudo-spectral methods. 相似文献