共查询到20条相似文献，搜索用时 145 毫秒
1.
Transient wave propagation in three-dimensional unbounded domains is studied. An efficient numerical approach is proposed, which is based on using the displacement unit-impulse response matrix representing the interaction force–displacement relationship on the near field/far field interface. Spatially, an approximation is used to reduce the computational effort associated with the large size of three-dimensional problems. It is based on subdividing the fully coupled unbounded domain into multiple subdomains. The displacement unit-impulse response matrices of all subdomains are calculated separately. The error associated with this spatial decoupling can be reduced by placing the near field/far field interface further away from the domain of interest. Detailed parameter studies have been conducted using numerical examples, in order to provide guidelines for the proposed spatially local schemes, and to demonstrate the accuracy and high efficiency of the proposed method for three-dimensional soil–structure interaction problems. 相似文献
2.
An axisymmetric infinite element and a three-dimensional infinite element are developed to solve three-dimensional elastic wave propagation problems in unbounded media. The elements are capable of transmitting Rayleigh, shear and compressional waves in the frequency domain. A scheme to integrate numerically the characteristic matrices of the elements is formulated based upon Gauss—Laguerre quadrature. Finally, the axisymmetric infinite element is used to find the compliance functions of a rigid circular plate subjected to harmonic loading on a semi-infinite medium. By using infinite elements, the size of the near field may be kept small. Consequently, the system is characterized by relatively few degrees of freedom, thus providing the analyst with an inexpensive solution. 相似文献
3.
基于精细积分法的三维弹性波数值模拟(英文) 总被引：1，自引：0，他引：1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
4.
Elastic wave propagation in inhomogeneous anisotropic media 总被引：1，自引：0，他引：1
XIU CHENG WEI ) MIN YU DONG ) YUN YAI CHEN ) ) Geoscience Department University of Petroleum Beijing China ) Institute of Geophysics China Seismological Bureau Beijing China 《地震学报(英文版)》1998,11(6):655-667
ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅｄｉａｉｎｅａｒｔｈａｒｅｑｕｉｔｅｃｏｍｐｌｅｘ．Ｔｈｅｒｅｅｘｉｓｔｓｅｖｅｒａｌｕｎｃｏｎｔｉｎｕｏｕｓｐｌａｉｎｓ．Ｎｏｒｍａｌｙ,ｉｔｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｒｅａｌｍｅｄｉｕｍｗｉｔｈｌａｙ... 相似文献
5.
John P. Wolf 《地震工程与结构动力学》1986,14(4):655-673
A non-linear interaction analysis with a (generalized) non-linear structure and a linear unbounded soil is analysed in the time domain, based either on the sub-structure method, which involves global convolution integrals, or on the direct method with local boundary conditions. Alternatively, the hybrid frequency–time-domain method of analysis, which is an iterative scheme, could be used. Approximate local boundary conditions to model the wave propagation towards infinity on the artificial boundary used in the direct method of non-linear soil–structure-interaction analysis to be performed in the time domain are examined. A semi-infinite rod supported elastically, which exhibits the same properties as certain unbounded soils such as dispersion and a cut-off frequency, is used for the investigation. For a transient excitation, the superposition boundary with frequent averaging, the well-known viscous damper and the extrapolation algorithm lead to good accuracy. Moving the artificial boundary further away from the structure (or more precisely, increasing the ratio of the distance of the artificial boundary to the wave length) improves the accuracy. 相似文献
6.
Václav Vavryčuk 《Pure and Applied Geophysics》1994,143(4):563-586
Complete relations are derived for energy and energy flux of elastic waves generated by an isotropic and double-couple source in a perfectly elastic, homogeneous, isotropic, and unbounded medium. In the energy balance of elastodynamic sources near-field waves play an essential role, transforming static energy into wave energy, andvice versa. For explosive and dislocation sources, the sources surface radiates a positive wave energy that is partially distributed to the medium transforming into static energy. For implosive and antidislocation sources, the source surface generates elastic waves, but it does not necessarily imply that it also radiates a positive wave energy. The energy transported by waves can originate in gradual transformation of the static-to-wave energy during propagation of waves through a stressed medium.On leave from Geophysical Institute, Czech Academy of Sciences, Boní II/1401, 41 31, Praha 4 Czech Republic 相似文献
7.
Manuel Carbajal-Romero Esteban Flores-Mendez Norberto Flores-Guzmán Jaime Núñez-Farfán Enrique Olivera-Villaseñor Francisco José Sánchez-Sesma 《Geofísica Internacional》2013,52(1):21-30
The present work shows the propagation of Scholte interface waves at the boundary of a fluid in contact with an elastic solid, for a broad range of solid materials. It has been demonstrated that by an analysis of diffracted waves in a fluid it is possible to infer the mechanical properties of the elastic solid medium, specifically, its propagation velocities. For this purpose, the diffracted wave field of pressures and displacements, due to an initial wave of pressure in the fluid, are expressed using boundary integral representations, which satisfy the equation of motion. The source in the fluid is represented by a Hankel’s function of second kind and zero order. The solution to this wave propagation problem is obtained by means of the Indirect Boundary Element Method, which is equivalent to the well-known Somigliana representation theorem. The validation of the results is carried out by using the Discrete Wave Number Method and the Spectral Element Method. Firstly, we show spectra of pressures that illustrate the behavior of the fluid for each solid material considered, then, we apply the Fast Fourier Transform to show results in time domain. Snapshots to exemplify the emergence of Scholte’s waves are also included. 相似文献
8.
D. K. Wagh 《Pure and Applied Geophysics》1970,82(1):62-65
Summary Cauchy theory of initial stress has been applied to investigate the influence of uniform initial stress on the propagation of plane waves in an unbounded elastic medium. It is found that if the initial principal stresses are unequal, the velocity depends on the direction in which the wave propagates, whereas when the initial principal stresses are equal, the velocity is independent of the direction of wave propagation. 相似文献
9.
应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。 相似文献
10.
Myung Hoon Kim Yunseok Choi Young Ho Cha Changsoo Shin 《Pure and Applied Geophysics》2009,166(12):1967-1985
In order to correctly interpret marine exploration data, which contain many elastic signals such as S waves, surface waves and converted waves, we have developed both a frequency-domain modeling algorithm for acoustic-elastic coupled media with an irregular interface, and the corresponding waveform inversion algorithm. By applying the continuity condition between acoustic (fluid) and elastic (solid) media, wave propagation can be properly simulated throughout the coupled domain. The arbitrary interface is represented by tessellating square and triangular finite elements. Although the resulting complex impedance matrix generated by finite element methods for the acoustic-elastic coupled wave equation is asymmetric, we can exploit the usual back-propagation algorithm used in the frequency domain through modern sparse matrix technology. By running numerical experiments on a synthetic model, we demonstrate that our inversion algorithm can successfully recover P- and S-wave velocity and density models from marine exploration data (pressure data only). 相似文献
11.
12.
Based on the u–U formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A u–U formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source. 相似文献
13.
Yongming Lu Qiancheng Liu Jianfeng Zhang Kai Yang Hui Sun 《Geophysical Prospecting》2019,67(5):1296-1311
With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples. 相似文献
14.
15.
模拟弹性波在孔隙介质中传播，需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上，首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法，并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先，本文通过均匀孔隙模型的数值解与解析解的对比，验证所提出的数值方法的正确性；然后，本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能，将完全匹配层与廖氏和阻尼吸收边界进行了对比，研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明，在孔隙介质中，完全匹配层作为吸收边界能十分有效地吸收衰减外行波，无论对体波还是面波，是一种高效边界吸收算法. 相似文献
16.
A method for determination of characteristics of quasi-Rayleigh (qR) wave in a transversely isotropic homogeneous half-space with inclined axis of symmetry is outlined. The solution is obtained as a superposition of qP, qSV and qSH waves, and surface wave velocity is determined from the boundary conditions at the free surface and at infinity, as in case of Rayleigh wave in an isotropic half-space. Though the theory is simple enough, a numerical procedure for calculation of surface wave velocity presents some difficulties. The difficulty is caused by necessity to calculate complex roots of a non-linear equation, which in turn contains functions determined as roots of non-linear equations with complex coefficients. Numerical analysis shows that roots of the equation corresponding to the boundary conditions do not exist in the whole domain of azimuths and inclinations of the symmetry axis. The domain of existence of qR wave depends on the ratio of the elastic parameters: for some strongly anisotropic models the wave cannot exist at all. For some angles of inclination qR-wave velocities deviate from those calculated on the basis of the perturbation method valid for weak anisotropy, though they have the same tendency of variation with azimuth. The phase of qR wave varies with depth unlike Rayleigh wave in an isotropic half-space. Unlike Rayleigh wave in an isotropic half-space, qR wave has three components - vertical, radial and transverse. Particle motion in horizontal plane is elliptic. Direction of the major axis of the ellipsis coincides with the direction of propagation only in azimuths 0° (180°) and 90° (270°). 相似文献
17.
L. ALFANO 《Geophysical Prospecting》1973,21(3):484-525
Numerical solutions of the wave equation for a particular type of non-linear “constant Q” medium were carried out. These solutions were obtained after the transformation of the space derivatives in finite differences; power series of the time are used to express the solutions. The medium is characterized by a not single valued stress-strain relation; the stresses are greater when the absolute values of strain are increasing (passive work), and are less when they are decreasing (active work). A loss of energy follows which is constant for every cycle and independent of frequency. This model represents the simplest type of medium in agreement with the laboratory data on rock samples. For a similar medium the stress’values do not depend only on the instantaneous value of the strain, but also on the previous strain values, i.e. the history of the medium. All these characteristics are not compatible with linearity and require particular types of computing procedures similar to the one shown in this paper. The results of calculations here shown refer both to the propagation of an isolated wave and to the generation of a wave in a spherical hole by change of pressure. They refer particularly to the shape, the amplitude and the width of the isolated wave along the propagation path. The most important results for this type of medium are the following: a) for a plane single isolated wave, the displacement amplitude wave does not change along the propagated distance. The width increases linearly as function of the distance; b) the corresponding particle velocity decreases in amplitude; c) for single isolated spherical waves the displacement amplitude decreases with propagated distance only due to the geometric factor, i.e. inversely proportional to the propagated distance; its width increases in the same way as for plane waves. The comparison between these theoretical results with the field and seismological data show a sufficiently good agreement as far as the value interval of wave width and frequencies is concerned. Less satisfactory is the comparison regarding a linear dependence of the wave width on the distance. This fact happens probably because in the field often long trains of waves and not isolated waves occur. In effect, for trains of waves the behaviour is different than that of an isolated wave; particularly, for the former the frequency variations along the travelled path is less and the displacement variations greater. However, it seems likely that a further similar theoretical research for trains of waves propagating in this type of non-linear medium might be carried out to complete the present research. 相似文献
18.
The basic equation of motion to analyse the interaction of a non-linear structure and an irregular soil with the linear unbounded soil is formulated in the time domain. The contribution of the unbounded soil involves convolution integrals of the dynamic-stiffness coefficients in the time domain and the corresponding motions. Alternatively, a flexibility formulation for the contribution of the unbounded soil using the dynamic-flexibility coefficients in the time domain, together with the direct-stiffness method for the structure and the irregular soil can be applied. The dynamic-stiffness or flexibility coefficient in the time domain is calculated as the inverse Fourier transform of the corresponding value in the frequency domain. The dynamic-stiffness coefficient's asymptotic behaviour for high frequencies determines the singular part whose transformation exists only in the sense of a distribution. As the dynamic-flexibility coefficient converges to zero for the frequency approaching infinity, the corresponding coefficient in the time domain is simpler to calculate, as no singular part exists. The salient features of the dynamic-stiffness and flexibility coefficients in the time domain are illustrated using a semi-infinite rod with exponentially increasing area. The dynamic-flexibility coefficients in the time domain are calculated for a rigid circular disc resting on the surface of an elastic halfspace and of a layer built-in at its base. Material damping is also introduced using the three-parameter Kelvin and the Voigt models. 相似文献
19.
Munirdin Tohti Yibo Wang Evert Slob Yikang Zheng Xu Chang Yi Yao 《Geophysical Prospecting》2020,68(6):1927-1943
Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ. 相似文献
20.
Examples of non-linear wave propagation in an elasto-plastic building are presented for excitation by pulses of strong ground motion characteristic of the near-field shaking near earthquake faults. Conditions that lead to the occurrence of permanent deformations in the building are investigated, and the amplitudes and wavelengths of incident pulses that lead to non-linear response are shown. Because the building can fail during the first passage of the incident wave pulse up and down the building (during a period that is shorter than the first natural period of the building), it is concluded that for the analysis and the design of structures in the near-field of earthquake shaking the wave propagation method of analysis must be used in place of the response spectrum method, which is based on the vibrational solution of the same governing equations. 相似文献