共查询到20条相似文献，搜索用时 31 毫秒
1.
流体饱和多孔介质黏弹性动力人工边界 总被引：1，自引：0，他引：1
基于Biot流体饱和多孔介质本构方程，采用平面波和远场散射波经验叠加来反映外行波传播，以经验参数反映人工边界外行波动的衰减和多角度透射特性。在人工边界处分别施加反映固相和液相介质传播效应的弹簧及阻尼来模拟人工边界以外的无限域介质对来自有限域的外行波的能量的吸收作用。从而形成一种流体饱和多孔介质的黏弹性动力人工边界。数值算例表明：边界的精度和稳定性高于现有的黏性边界、黏弹性人工边界及一阶透射边界。 相似文献
2.
T. AkiyoshiK. FuchidaH. L. Fang 《Soil Dynamics and Earthquake Engineering》1994,13(6):387-397
Based on Biot's two-phase mixture theory and the paraxial approximation, the absorbing boundary conditions in the time domain for u-w, u-U and u-p formulations are presented for the dynamic analysis of fluid-saturated porous media. These absorbing boundaries are equivalent to the viscous boundaries in the fundamental mode. The expressions for the energy ratio and reflection coefficient between the reflected and incident waves along the absorbing boundary are given. The numerical results show that the proposed absorbing boundaries can greatly suppress the spuriously reflected wave and model the far field of the system. These results also dynamic analysis of infinite fluid-saturated porous media. for the transient dynamic analysis of infinite fluid-saturated porous media. 相似文献
3.
4.
针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题。以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型。通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱和多孔热弹性介质边界上的四种反射波的振幅反射率的理论表达式。在此基础上进行数值计算,分别讨论平面S波的入射频率、入射角和热膨胀系数等参数对四种反射波的振幅反射率的影响情况。结果表明:各反射波的振幅反射率分别随频率和热膨胀系数的增大而增大,同时也受到平面S波入射角变化的影响。该结论对于土动力学的理论研究及其相关的工程勘探具有一定的指导意义。 相似文献
5.
Jianwen Liang Hongbing You Vincent W. Lee 《Soil Dynamics and Earthquake Engineering》2006,26(6-7):611-625
The scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space is modeled using the indirect boundary element method in the frequency domain. The free-field responses are calculated to determine the displacements and stresses at the surface of the canyon, and fictitious distributed loads are then applied at the surface of the canyon in the free field to calculate the Green's functions for displacements and stresses. The amplitudes of the fictitious distributed loads are determined from the boundary conditions, and the displacements arising from the waves in the free field and from the fictitious distributed loads are summed to obtain the solution. The effects of fluid saturation, boundary conditions, porosity, and soil layers on the surface displacement amplitudes and phase shifts are discussed, and some useful conclusions are obtained. It is shown that the surface displacement amplitudes due to saturation and boundary conditions, different porosities, or the presence of a soil layer can be very dissimilar, and large phase shifts can be observed. The resulting wavelengths for an undrained saturated poroelastic medium are slightly longer than those for a drained saturated poroelastic medium; and are longer for a drained saturated poroelastic medium than those for a dry poroelastic medium. As porosity increases, the wavelengths become longer; and a layered half-space produces longer wavelengths than a homogeneous half-space. 相似文献
6.
Diffraction of plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space 总被引：7，自引：0，他引：7
Jianwen Liang Zhenning Ba Vincent W. Lee 《Soil Dynamics and Earthquake Engineering》2006,26(6-7):582-610
An analytical solution for the scattering and diffraction of incident plane SV waves by a shallow circular-arc canyon in a saturated poroelastic half-space is derived by the wave function expansion method. The solution is utilized to analyze the dependence of the computed surface motions on the incident frequencies, incident angles, porosity, boundary drainage and Poisson's ratio. It is shown that, depending on the incident angles, the surface displacement amplitudes around a canyon in a dry poroelastic half-space and saturated poroelastic half-space can be very different. The surface displacement amplitudes of an undrained saturated poroelastic half-space are close to those of a drained saturated poroelastic half-space. For low porosity, the surface displacement amplitudes of a saturated poroelastic half-space are almost identical to those of a dry poroelastic half-space, and drainage condition has little influence on the surface displacement amplitudes. But for high porosity, the effect of drainage condition becomes significant, and for the same porosity, the displacement amplitudes of an undrained saturated half-space will be larger than those of a drained saturated half-space. Poisson's ratio is also an important factor affecting the surface displacement amplitudes around the canyon, both in drained and undrained conditions, but leads to larger effects for an undrained saturated half-space than for a drained saturated half-space. Large pore pressures are found around the canyon and their amplitudes depend on the incident angles and frequencies. Below the surface, the amplitudes of pore pressures are less than they are at the surface, especially for high frequencies. 相似文献
7.
Amongst different existing techniques developed for absorbing boundaries in transient dynamic analysis, the paraxial approximation represents an elegant framework. This approach is extended to the case of saturated porous media, which is of vital interest in many earthquake engineering problems. With respect to the frequency content of the dynamic process, several formulations exist for modelling the two-phase behaviou of soil or rock materials. The existence of a shear wave and two dilatational waves, which are standard features of porous media, is verified for the adopted formulation. The paraxial approximation is presented in the case of the u – p Biot two-phase dynamic formulation, suitable for seismic analyses and implemented by using a finite element approach. The efficiency of these boundaries is further verified through same numerical examples. 相似文献
8.
The strong coupling of applied stress and pore fluid pressure, known as poroelasticity, is relevant to a number of applied problems arising in hydrogeology and reservoir engineering. The standard theory of poroelastic behavior in a homogeneous, isotropic, elastic porous medium saturated by a viscous, compressible fluid is due to Biot, who derived a pair of coupled partial differential equations that accurately predict the existence of two independent dilatational (compressional) wave motions, corresponding to in-phase and out-of-phase displacements of the solid and fluid phases, respectively. The Biot equations can be decoupled exactly after Fourier transformation to the frequency domain, but the resulting pair of Helmholtz equations cannot be converted to partial differential equations in the time domain and, therefore, closed-form analytical solutions of these equations in space and time variables cannot be obtained. In this paper we show that the decoupled Helmholtz equations can in fact be transformed to two independent partial differential equations in the time domain if the wave excitation frequency is very small as compared to a critical frequency equal to the kinematic viscosity of the pore fluid divided by the permeability of the porous medium. The partial differential equations found are a propagating wave equation and a dissipative wave equation, for which closed-form solutions are known under a variety of initial and boundary conditions. Numerical calculations indicate that the magnitude of the critical frequency for representative sedimentary materials containing either water or a nonaqueous phase liquid is in the kHz–MHz range, which is generally above the seismic band of frequencies. Therefore, the two partial differential equations obtained should be accurate for modeling elastic wave phenomena in fluid-saturated porous media under typical low-frequency conditions applicable to hydrogeological problems. 相似文献
9.
基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波－慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高. 相似文献
10.
Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with small-scale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference (FD) method of Biot’s poroelastic equations, with an arbitrary even-order (2L) accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer (CPML) absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultrasonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thicknesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Q sc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments. 相似文献
11.
The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly. 相似文献
12.
The phenomenon of reflection and transmission of plane harmonic waves at the plane interface between two dissimilar poroelastic solids saturated with two immiscible viscous fluids is investigated. Both porous media are considered dissipative due to the presence of viscosity in pore‐fluids. Four attenuated (three dilatational and one shear) waves propagate in such a dissipative porous medium. A finite non‐dimensional parameter is used to define the effective connections between the surface‐pores of two media at their common interface. Another finite parameter represents the gas‐share in the saturation of pores. An attenuated wave in a dissipative medium is described through the specification of directions of propagation and maximum attenuation. A general representation of an attenuated wave is defined through its inhomogeneous propagation, i.e., different directions for propagation and attenuation. Incidence of an inhomogeneous wave is considered at the interface between two dissipative porous solids. This results in four reflected and four transmitted inhomogeneous waves. Expressions are derived for the partition of incident energy among the reflected and transmitted waves. Numerical examples are studied to determine the effects of saturating pore fluid, frequency, surface‐pore connections and wave inhomogeneity on the strengths of reflected and transmitted waves. Interaction energy due to the interference of different (inhomogeneous) waves is calculated in both the dissipative porous media to verify the conservation of incident energy. 相似文献
13.
模拟弹性波在孔隙介质中传播，需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上，首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法，并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先，本文通过均匀孔隙模型的数值解与解析解的对比，验证所提出的数值方法的正确性；然后，本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能，将完全匹配层与廖氏和阻尼吸收边界进行了对比，研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明，在孔隙介质中，完全匹配层作为吸收边界能十分有效地吸收衰减外行波，无论对体波还是面波，是一种高效边界吸收算法. 相似文献
14.
An analytical model for describing the propagation and attenuation of Rayleigh waves along the free surface of an elastic porous medium containing two immiscible, viscous, compressible fluids is developed in the present study based on the poroelastic equations formulated by Lo et al. [Lo WC, Sposito G, Majer E. Wave propagation through elastic porous media containing two immiscible fluids. Water Resour Res 2005;41:W02025]. The dispersion equation obtained is complex-valued due to viscous dissipation resulting from the relative motion of the solid to the pore fluids. As an excitation frequency is stipulated, the dispersion equation that is a cubic polynomial is numerically solved to determine the phase speed and attenuation coefficient of Rayleigh waves in Columbia fine sandy loam permeated by an air–water mixture. Our numerical results show that, corresponding to three dilatational waves, there is also the existence of three different modes of Rayleigh wave in an unsaturated porous medium, which are designated as the R1, R2, and R3 waves in descending order of phase speed, respectively. The phase speed of the R1 wave is non-dispersive (frequency-independent) in the frequency range we examined (10 Hz–10 kHz) and decreases as water saturation increases, whose magnitude ranges from 20% to 49% of that of the first dilatational wave with respect to water content. However, it is revealed numerically that the R2 and R3 waves are functions of excitation frequency. Given the same water saturation and excitation frequency, the phase speeds of the R2 and R3 waves are found to be approximately 90% of those of the second and third dilatational waves, respectively. The R1 wave has the lowest attenuation coefficient whereas the R3 wave attenuates highest. 相似文献
15.
Multi—transmitting formula for attenuating waves 总被引：1，自引：0，他引：1
The MTF is extended to case of attenuating incident wave by introducing an attenuation coefficient.The reflection coefficients of this modified MTF and MTF are evaluated and compared when an attenuating wave impinges on the boundary,and the results demonstrate that MTF can be used to absorb slightly attenuating waves and the modi-fied MTF is more capable of absorbing heavily attenuating waves than MTF.The accuracy of modified MTF is also tested by numerical examples of fluid saturated porous media. 相似文献
16.
A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media. 相似文献
17.
本文以饱水两相介质的土力学模型为研究对象,在假定两相介质为弹性介质条件下,采用了显式有限元法和透射边界进行了饱和弹性半空间动力响应问题的研究。为避免谐波输入初始间断的影响,文中提出了一个处理函数,并以弹性半空间为算例,对饱水介质和单相介质分别进行了在底边界P波垂直入射时的动力响应分析,验证了该处理函数的有效性和实用性。 相似文献
18.
Local fluid flow (LFF) at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band. LFF is easily influenced by the structure and boundary conditions of the porous media, which leads to different behaviors of the peak frequency of attenuation. The associated transition frequency can provide detailed information about the trend of LFF; therefore, research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation (i.e., inverse of quality factor) facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media. In this study, we firstly obtain the transition frequency of fluid flux based on Biot’s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell (RUC). We then analyze changes of these two frequencies in porous media with different porous properties. Finally, we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs. This setup can facilitate the understanding of the effect from the undrained boundary condition. Results demonstrate that these two frequencies have the same trend at low water saturation, but amplitude variations differ between the frequencies as the amount of saturation increases. However, for cases of high water saturation, both the trend and the amplitude variation of these two frequencies fit well with each other. 相似文献
19.
平面P波在饱和半空间中洞室周围的散射(I):解析解 总被引：1，自引：0，他引：1
利用波函数展开法给出了入射平面P波在饱和半空间中圆形洞室周围散射问题的一个解析解。半空间假定为无粘性流体饱和介质,满足Biot理论。采用一种基于实验数据的孔隙率和模量之间的线性关系来确定Biot模型中的介质参数。解答考虑了透水边界和非透水边界两种情况。对边界条件进行了数值检验,结果表明,随着级数截断项数的增大,边界残量衰减很快。解答为进一步研究入射波频率和角度、边界渗透条件、孔隙率、泊松比等参数对散射的影响奠定了基础。 相似文献
20.
This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated (either drained or undrained) cases. There are large phase... 相似文献