共查询到20条相似文献，搜索用时 31 毫秒
1.
Kui Xiang Kjersti Solberg Eikrem Morten Jakobsen Geir Nvdal 《Geophysical Prospecting》2022,70(1):3-18
We have derived a convergent scattering series solution for the frequency-domain wave equation in acoustic media with variable density and velocity. The convergent scattering series solution is based on the homotopy analysis of a vectorial integral equation of the Lippmann–Schwinger type. By using the Green's function and partial integration, we have derived the vectorial integral equation of the Lippmann–Schwinger type that involves the pressure gradient field as well as the pressure field from the wave equation. The vectorial Lippmann–Schwinger equation can in principle be solved via matrix inversion, but the computational cost of matrix inversion scales like , where is the number of grid blocks. The computational cost can be significantly reduced if one solves the vectorial Lippmann–Schwinger equation iteratively. A simple iterative solution is the Born series, but it is only convergent when the scattering potential is sufficiently small. In this study, we have used the so-called homotopy analysis method to derive an iterative solution for the vectorial Lippmann–Schwinger equation which can be made convergent even in strongly scattering media. The computational cost of our convergent scattering series scales as . Our algorithm, which is based on the homotopy analysis method, involves a convergence control operator that we select using hierarchical matrices. We use a three-layer model and a resampled version of the SEG/EAGE salt model to show the performance of the developed convergent scattering series. 相似文献
2.
Wei-Lung Chang 《地震工程与结构动力学》1993,22(9):773-789
An alternate formulation of the ‘substructure deletion method’ suggested by Dasgupta in 19791 has been successfully implemented. The idea is to utilize simple Green's functions developed for a surface problem to replace the more complicated Green's functions required for embedded problems while still being able to generate an accurate solution. Since the exterior medium is usually represented by a continuum model, the interior medium in the present approach will also be represented by a continuum model rather than a finite element model as suggested originally, thereby eliminating the incompatibility between the solutions of the interior and exterior media. Detailed studies of the method's accuracy and limitations were performed using two-dimensional examples in wave scattering of canyons and alluvial valleys, problems which are more suitable for this method than the embedded foundation problem. The results obtained indicate that the alternate formulation gives accurate results only when the vertical dimension of the scattering object is not too large; if the aspect ratio (vertical over lateral) exceeds a certain limit, the results will not approach the known results given by boundary integral equation solutions or indirect boundary integral equations no matter what the refinement of the model may be. The greatest advantage of the present method is that the task of calculating Green's functions is reduced significantly; computational time using this new formulation is approximately five times less than for conventional boundary integral equation methods. 相似文献
3.
We study singular elastic solutions at an angular corner left by a crack that has kinked. We have in mind a geophysical context
where the faults on either side of the kink are under compression and are ready to slip, or have already slipped, under the
control of Coulomb friction. We find separable static singular solutions that are matched across the sides of the corner by
applying appropriate boundary conditions. In our more general solution we assume that one of the sides of the corner is about
to slide, i.e. it is just contained by friction, and the other may be less pressured. Our solutions display power law behaviour
with real exponents that depend continuously on the angle of the corner, the coefficient of static friction and the difference
of shear load on both sides of the corner. When friction is the same on both sides of the kink, the solutions split into a
symmetric and an antisymmetric solution. The antisymmetric solution corresponds to the simple shear case; while the symmetric
solution appears when the kink is loaded by uniaxial stress along the bisector of the kink. The antisymmetric solution is
ruled out under this model with contact since the faults cannot sustain tension. When one side of the corner is less pressured
one can also distinguish modes with contact overall from others that must open up on one side. These solutions provide an
insight into the stress distributions near fault kinks, they can also be used as tools for improving the numerical calculation
of kinks under static or dynamic loads. 相似文献
4.
M. AL-CHALABI 《Geophysical Prospecting》1973,21(4):783-795
A rigorous proof is presented to show that over a horizontally layered ground the rms velocity cannot exceed the stacking velocity. The proof helps to illustrate the difference between stacking and rms velocities in a quantitative manner. The series of Taner and Koehler (1969) is used for the purpose. Convergence of this series is tested. Including more terms will not necessarily improve the convergence. Although the series is rapidly convergent when the spread length/depth ratio is small, strong oscillations are observed when this ratio is high. 相似文献
5.
Although based on exact analytical solutions, semi‐analytical solute transport models can have significant numerical error in applications with high frequency oscillatory source terms and when parameter value combinations cause series solution approximations to converge slowly. Methods for correcting these numerical errors are presented and implemented in the AT123D code, which employs Green's functions to represent point, linear, and rectangular prismatic source zones. In order to increase its computational accuracy, a Romberg numerical integration scheme was added to AT123D with prespecified error criteria, variable time stepping, and partitioning of the integral to handle rapidly changing source terms. More rapidly converging series solution approximations for the Green's functions were also incorporated to improve both accuracy and computational efficiency for finite‐depth aquifers. AT123D also has been modified to eliminate redundant calculations at points where approximate steady‐state conditions have been reached to improve computational efficiency during numerical integration. These modifications help to decrease computer run times that can be excessive for three‐dimensional problems with large numbers of computational points, small time steps, and/or long simulation time periods. Errors in the original AT123D code also were corrected in this modified version, AT123D‐AT, in order to accurately simulate finite‐duration (pulse) source releases. 相似文献
6.
The constant-head pumping tests are usually employed to determine the aquifer parameters and they can be performed in fully or partially penetrating wells. Generally, the Dirichlet condition is prescribed along the well screen and the Neumann type no-flow condition is specified over the unscreened part of the test well. The mathematical model describing the aquifer response to a constant-head test performed in a fully penetrating well can be easily solved by the conventional integral transform technique under the uniform Dirichlet-type condition along the rim of wellbore. However, the boundary condition for a test well with partial penetration should be considered as a mixed-type condition. This mixed boundary value problem in a confined aquifer system of infinite radial extent and finite vertical extent is solved by the Laplace and finite Fourier transforms in conjunction with the triple series equations method. This approach provides analytical results for the drawdown in a partially penetrating well for arbitrary location of the well screen in a finite thickness aquifer. The semi-analytical solutions are particularly useful for the practical applications from the computational point of view. 相似文献
7.
在层状半空间精确动力刚度矩阵和斜线荷载动力格林函数的基础上建立间接边界元方法,在频域内求解无破碎带断层场地对入射平面SH波的散射。为方便求解,将总波场分解为自由波场和散射波场,自由波场由直接刚度法求得,断层两侧的散射波场通过在断层面上分别对两侧施加均布斜线荷载产生的动力响应来模拟,虚拟荷载的密度可通过引入断层表面的边界条件确定,最后叠加自由波场和散射波场求得总波场。以有落差断层和无落差断层模型为例进行数值计算,分析断层落差、断层倾角以及断层两侧介质的刚度比对散射效应的影响。研究表明,断层落差与波长相当时,断层对SH波的放大作用最大;地表位移幅值随着断层倾角的增大逐渐增大;若断层无落差且其两侧刚度不同时,一般刚度较小一侧地表位移幅值较大且振荡更为剧烈,波从刚度较小一侧入射时位移幅值放大尤为显著。 相似文献
8.
In any numerical solution of the DC resistivity experiment, care must be taken to deal with strong heterogeneity of electrical conductivity. In order to examine the importance of conductivity contrasts, we develop a scattering decomposition of the DC resistivity equation in the sparse differential domain as opposed to the traditional dense integral formulation of scattering‐type equations. We remove the singularity in the differential scattered series via separation of primary and secondary conductivity, thereby avoiding the need to address the singularity in a Green's function. The differential scattering series is observed to diverge for large conductivity contrasts and to converge for small contrasts. We derive a convergence criterion, in terms of matrix norms for the weak‐form finite‐volume equations, that accounts for both the magnitude and distribution of heterogeneity of electrical conductivity. We demonstrate the relationship between the differential scattering series and the Fréchet derivative of the electrical potential with respect to electrical conductivity, and we show how the development may be applied to the inverse problem. For linearization associated with the Fréchet derivative to be valid, the perturbation in electrical conductivity must be small as defined by the convergence of the scattered series. The differential scattering formulation also provides an efficient tool for gaining insight into charge accumulation across contrasts in electrical conductivity, and we present a derivation that equates accumulated surface charge density to the source of scattered potential. 相似文献
9.
Derived herein is the integral representation solution of a Rayleigh-damped Bernoulli–Euler beam subjected to multi-support motion, which is free from calculation of a quasi-static solution, and in which the modal participation factor for support motion is formulated as a boundary modal reaction, thus making efficient calculation feasible. Three analytical methods, including (1) the quasi-static decomposition method, (2) the integral representation with the Cesàro sum technique, and (3) the integral representation in conjunction with Stokes' transformation, are presented. Two additional numerical methods of (4) the large mass FEM simulation technique and (5) large stiffness FEM simulation technique are easily incorporated into a commercial program to solve the problem. It is found that the results obtained by using these five methods are in good agreement, and that both the Cesàro sum and Stokes' transformation regularization techniques can extract the finite part of the divergent series of the integral representation. In comparison with the Mindlin method and Cesàro sum technique, Stokes' transformation is the best way because it is not only free of calculation of the quasi-static solution, but also because it can obtain the convergence rate as rapidly as the mode acceleration method can. 相似文献
10.
Piecewise heterogeneous media that the earth presents are composed of large-scale boundary structures and small-scale volume heterogeneities. Wave propagation in such piecewise heterogeneous media can be accurately superposed through the generalized Lippmann–Schwinger integral equation (GLSIE). Two different Born series modeling schemes are formulated for the boundary–volume integral equation with 2-D antiplane motion (SH waves). Both schemes decompose the resulting boundary–volume integral equation matrix into two parts: the self-interaction operator handled with a fully implicit manner, and the extrapolation operator approximated by a Born series. The first scheme associates the self-interaction operator with each boundary itself and the volume itself, and interprets the extrapolation operator as the cross-interaction between each boundary and other boundaries/volume scatterers in a subregion. The second scheme relates the self-interaction operator to each boundary itself and its cross-interaction with the volume scatterers on both sides, and expresses the extrapolation operator as both the direct and indirect (through the volume scatterers) cross-interactions between different boundaries in a subregion. By eliminating the displacement field from the volume scatterers, the second scheme reduces the dimension of the resulting boundary-volume integral equation matrix, leading to a faster convergence than the first scheme. Both the numerical schemes are validated by dimensionless frequency responses to a heterogeneous alluvial valley with the velocity perturbed randomly in the range of ca 5–20 %. The schemes are applied to wave propagation simulation in a heterogeneous multilayered model by calculating synthetic seismograms. Numerical experiments, compared with the full-waveform numerical solution, indicate that the Born series modeling schemes significantly improve computational efficiency, especially for high frequencies. 相似文献
11.
Francisco Luzón Francisco J. Sánchez-Sesma J. Alfonso Pérez-Ruiz Leonardo Ramírez-Guzmán Andrés Pech 《Soil Dynamics and Earthquake Engineering》2009
The propagation of seismic P and SV waves within inhomogeneous alluvial valleys has been investigated using the indirect boundary element method (IBEM). An improvement on the formulae published early, for the 2D Green's functions in an inhomogeneous medium is presented in this work. A modification has been done over these functions in view of its connection to the ray theory. An accuracy analysis validates these modified Green's functions computing the relative error in frequency domain for the SH case, and with a quantitative analysis obtaining envelope and phase misfits of the solution in time domain, for the P–SV case. 相似文献
12.
A direct boundary element procedure is presented to determine the impedance matrix for a three-dimensional foundation supported on an infinitely-long canyon of uniform cross-section cut in a homogeneous half-space. The uniform cross-section of the canyon permits analytical integration along the canyon axis leading to a series of two-dimensional boundary problems involving Fourier transforms of the full-space Green's functions. Solution of these two-dimensional boundary problems leads to a dynamic flexibility influence matrix which is inverted to determine the impedance matrix. The accuracy of the procedure is demonstrated by comparison with previous solutions for a surface-supported, square foundation and results obtained by a three-dimensional boundary element method (BEM) for a foundation of finite-width supported on an infinitely-long canyon. Compared with the three-dimensional BEM, the present method requires less computer storage and is more accurate and efficient. The foundation impedance matrix determined by this procedure can be incorporated in a substructure method for earthquake analysis of arch dams. 相似文献
13.
传统的完全匹配层技术是一种能够较为有效地消除边界反射的边界条件,但是当表层为泊松比较高的自由表面时,该技术可能会产生不稳定的现象.针对传统的完全匹配层技术固有的不稳定和掠射情况下吸收效果不佳等缺陷,发展了多轴完全匹配层、卷积完全匹配层以及将两者结合的多轴卷积完全匹配层等3种边界条件.本文介绍了水平自由表面的不同处理方法以及传统、多轴、卷积和多轴卷积等4种完全匹配层条件的原理,通过二维半无限空间模型的交错网格有限差分正演模拟对比,分析了几种自由边界实施方法在这几种完全匹配层条件下的稳定性,并通过提取单道波形与解析解进行对比,定性分析了水平自由表面几种不同处理方法的准确性以及各自的适用条件. 结果表明,泊松比和水平自由表面实施方法对波场模拟效果及其稳定性有重要影响. 相似文献
14.
Apart from some special cases, calculating the dynamic stiffness matrix of foundations on a layered half-space, especially for embedded foundations, is computationally expensive. An efficient method for two-dimensional foundations in a horizontally layered soil media is presented in this paper. This method is based on indirect boundary element methods and uses discrete wave number solution methods for calculating Green's functions for displacements and analytical methods for the integrations over the boundary. For surface foundations, the present method applies at all frequencies. For embedded foundations or for constructing energy transmitting boundaries, because the free-field part is modelled by boundary elements and the excavated part is modelled by finite elements, the present method applies only at low frequencies for the spring coefficients (the real parts of the dynamic stiffness matrix) but applies at all frequencies for the damping coefficients (the imaginary part of the dynamic stiffness matrix) for undamped sites. The novelty of the method can be used for three-dimensional foundations. © 1997 by John Wiley & Sons, Ltd. 相似文献
15.
The evidence of east-west compression in northeast Japan has been reported by many investigators on the basis of geodetic,
geologic and geomorphic data, but its origin still remains far from understood. In the present study we have proposed a mechanical
model of tectonic loading at convergent plate boundary zones, and demonstrated its validity through the numerical simulation
of internal stress fields in northeast Japan with realistic 3-D geometry of plate interfaces. At convergent plate boundary
zones, in general, a part of plate convergence is consumed by steady slip along plate interfaces, and the remaining part by
inelastic deformation (seismic faulting, aseismic faulting, and active folding) of overriding plates. Such a plate boundary
process to be called ``partial collision' can be quantitatively described by introducing a collision rate defined as c = 1 − steady slip rate at plate interfaces/plate convergence rate. By this definition, we can simply represent the mechanical process of partial collision, which includes total subduction
(c = 0) and total collision (c = 1) as two extreme cases, in terms of steady slip rates at plate interfaces. On the basis of elastic dislocation theory,
first, we numerically computed the internal stress fields in northeast Japan produced by the total subduction of the Pacific
plate beneath the North American plate, however the computed stress pattern was opposite in sense to observations. Then, we
computed the internal stress fields by taking c = 0.1 on average, and succeeded in reproducing the observed east-west compression in northeast Japan. This indicates that
the concept of partial collision is essential to understand the mechanism of intraplate tectonic loading. 相似文献
16.
W.A. Mulder 《Geophysical Prospecting》2020,68(9):2857-2866
Simulations of wave propagation in the Earth usually require truncation of a larger domain to the region of interest to keep computational cost acceptable. This introduces artificial boundaries that should not generate reflected waves. Most existing boundary conditions are not able to completely suppress all the reflected energy, but suffice in practice except when modelling subtle events such as interbed multiples. Exact boundary conditions promise better performance but are usually formulated in terms of the governing wave equation and, after discretization, still may produce unwanted artefacts. Numerically exact non-reflecting boundary conditions are instead formulated in terms of the discretized wave equation. They have the property that the numerical solution computed on a given domain is the same as one on a domain enlarged to the extent that waves reflected from the boundary do not have the time to reach the original truncated domain. With a second- or higher-order finite-difference scheme for the one-dimensional wave equation, these boundary conditions follow from a recurrence relation. In its generalization to two or three dimensions, a recurrence relation was only found for a single non-reflecting boundary on one side of the domain or two of them at opposing ends. The other boundaries should then be zero Dirichlet or Neumann. If two non-reflecting boundaries meet at a corner, translation invariance is lost and a simple recurrence relation could not be found. Here, a workaround is presented that restores translation invariance by imposing classic, approximately non-reflecting boundary conditions on the other sides and numerically exact ones on the two opposing sides that otherwise would create the strongest reflected waves with the classic condition. The exact ones can also be applied independently. As a proof of principle, the method is applied to the two-dimensional acoustic wave equation, discretized on a rectangular domain with a second-order finite-difference scheme and first-order Enquist–Majda boundary conditions as approximate ones. The method is computationally costly but has the advantage that it can be reused on a sequence of problems as long as the time step and the sound speed values next to the boundary are kept fixed. 相似文献
17.
Serdar Korkmaz 《Ground water》2013,51(3):432-441
In this study, the well‐known Hantush solution procedure for groundwater mounding under infinitely long infiltration strips is extended to finite and semi‐infinite aquifer cases. Initially, the solution for infinite aquifers is presented and compared to those available in literature and to the numerical results of MODFLOW. For the finite aquifer case, the method of images, which is commonly used in well hydraulics, is used to be able to represent the constant‐head boundaries at both sides. It is shown that a finite number of images is enough to obtain the results and sustain the steady state. The effect of parameters on the growth of the mound and on the time required to reach the steady state is investigated. The semi‐infinite aquifer case is emphasized because the growth of the mound is not symmetric. As the constant‐head boundary limits the growth, the unbounded side grows continuously. For this reason, the groundwater divide shifts toward the unbounded side. An iterative solution procedure is proposed. To perform the necessary computations a code was written in Visual Basic of which the algorithm is presented. The proposed methodology has a wide range of applicability and this is demonstrated using two practical examples. The first one is mounding under a stormwater dispersion trench in an infinite aquifer and the other is infiltration from a flood control channel into a semi‐infinite aquifer. Results fit very well with those of MODFLOW. 相似文献
18.
Our aim is to introduce the Coupled Finite-Infinite Element Method (CFIEM) as a new alternative approach to the Earth’s gravity
field modelling. We show that if the computational domain is large enough in radial direction, one can obtain the qualitatively
and quantitatively comparable solution to the solution by the Finite Element Method (FEM). We study the influence of the size
of the computational domain on the final CFIEM solution as well as the successive refinement of the discretization and its
convergence to the exact solution. As an input data we use the synthetic boundary conditions computed from a Synthetic Earth
Gravity Model (SEGM) and we test the CFIEM solution by the data generated directly from SEGM and the solution by the FEM. 相似文献
19.
David P Thambiratnam 《地震工程与结构动力学》1986,14(3):475-485
Herein the propagation of transient waves in a rod due to impulsive end loading is investigated using the method of wavefront expansion. The Mindlin–Herrmann equations which include the Poisson effect are used to model the rod. Since this rod theory is restricted to lower values of Poisson's ratio (< 0.283), the present work also has the same restriction. The impulsive end loading on the rod can be prescribed in the form of stress, strain, velocity or acceleration boundary conditions. The analysis is based on the concept of a wave as a carrier of discontinuities in the field variables and their derivatives. These discontinuities are determined from a set of recurrence relations which are in turn generated by the use of time harmonic asymptotic series solutions to the equations of motion. Transients due to both velocity and strain (or stress) boundary conditions are treated. Numerical examples are presented to illustrate the method of solution. The results indicate the distortion of the transient waves due to the dispersive nature of the governing equations. Transients due to other boundary conditions at the rod end can be treated by using the results presented in the paper. 相似文献
20.
快速模拟退火算法中,要求扰动后新的模型参数必须在取值范围内,对此,人们通常采用两种方法:(1)反复扰动直到新的模型参数满足要求;(2)当新的模型参数超过最大值时,参数等于最大值;新的模型参数小于最小值时,参数等于最小值.研究发现,方法1在计算过程中增加了寻找新模型参数的时间开销,计算效率明显降低;方法2当目标函数在边界或边界的小邻域存在多个极小值时,在有限的时间内,其陷入局部极小的概率大大增加,优化效果明显下降.本文在前人研究的基础上提出了改进的模拟退火算法,同时利用改进的算法反演了地层泊松比、纵波速度和密度.结果表明,改进的算法加快了模拟退火算法的收敛速度,提高了模拟退火算法的效率和精度. 相似文献