共查询到20条相似文献，搜索用时 843 毫秒
1.
两相饱和多孔介质的基本方程组和计算参数的选取在一定程度上是混乱的.本文利用连续介质力学理论,结合空间平均化方法,根据应力-应变关系、运动学关系、连续性方程及广义Darcy定律,建立了横观各向同性液体饱和多孔介质的基本方程组;通过将该基本方程组与Biot理论的基本方程组进行比较,确定了Biot方程组中的弹性常数Bi与弹性常数cij的关系,并得到了确定Biot参数mi和ri的计算公式——Biot参数与渗透率、孔隙率及黏性系数的关系;最后对几种特殊情形进行了讨论,给出了简化方程组. 相似文献
2.
本文针对CHAMP型卫星建立了顾及非线性改正的轨道扰动方程定轨理论与方法.首先从卫星运动的二阶微分方程出发,引入了正常引力位以及相应的参考轨道,然后分别推导了线性化轨道扰动方程与顾及非线性改正的轨道扰动方程,同时说明了建立的线性化轨道扰动方程与目前处理CHAMP卫星数据的动力学定轨方法是等价的.其次分别对线性化轨道扰动方程与顾及非线性改正的轨道扰动方程的精度进行了估计,在卫星定位精度为3cm与非惯性力测量精度为3×10~(-10)m·s~(-2)的前提下证明了下列结论:当参考轨道与实际轨道之间的距离ρ≤4.7m时线性化轨道扰动方程的精度能达到非惯性力的测量精度以及当ρ≤4.14×10~3m时顾及非线性改正的轨道扰动方程能达到非惯性力的测量精度.由此便可得出结论:相对于线性化轨道扰动方程,顾及非线性改正的轨道扰动方程具有更高的精度,且适合在更长的时间弧段上建立关于引力场位系数的法方程组,特别是针对CHAMP卫星计划进行的模拟计算也完全验证了该结论.最后利用叠加原理,给出了顾及非线性改正的轨道扰动方程的求解方法.此外,还针对GRACE卫星计划利用顾及非线性改正的轨道扰动方程进行了恢复引力场的模拟计算,结果表明:分段建立位系数的法方程组时子弧段分别取值2h、1d、6 d对恢复引力场的结果几乎不产生影响,这表明在处理GRACE数据时能够以6d的弧长来建立法方程组. 相似文献
3.
Two linear-hysteretic-damping models that provide energy dissipation independent of the deformation frequency, are studied in this paper: a hysteretic Kelvin element and a hysteretic Maxwell element. Both models use the Hilbert transform and yield integro–differential equations for the equations of motion of structures when real-valued signals are utilized in the formulation. It is shown that the use of analytic (complex-valued) signals allows the transformation of these integro–differential equations into differential equations with analytic input signals and complex-valued coefficients. These differential equations show both stable and unstable poles. A technique for the solution of these differential equations is presented; it consists of a conventional modal decomposition of the state-space equations and the integration of the differential equations forward in time for the modal co-ordinates associated with stable poles, and backwards in time for the modal co-ordinates associated with unstable poles. Some numerical examples are presented to illustrate the characteristics of the models and the proposed analysis technique. 相似文献
4.
Jinghuai GAO Weimin HAN Yanbin HE Haixia ZHAO Hui LI Yijie ZHANG Zongben XU 《中国科学:地球科学(英文版)》2021,(3):377-387
Tight oil/gas medium is a special porous medium, which plays a significant role in oil and gas exploration. This paper is devoted to the derivation of wave equations in such a media, which take a much simpler form compared to the general equations in the poroelasticity theory and can be employed for parameter inversion from seismic data. We start with the fluid and solid motion equations at a pore scale, and deduce the complete Biot's equations by applying the volume averaging technique.The underlying assumptions are carefully clarified. Moreover, time dependence of the permeability in tight oil/gas media is discussed based on available results from rock physical experiments. Leveraging the Kozeny-Carman equation, time dependence of the porosity is theoretically investigated. We derive the wave equations in tight oil/gas media based on the complete Biot's equations under some reasonable assumptions on the media. The derived wave equations have the similar form as the diffusiveviscous wave equations. A comparison of the two sets of wave equations reveals explicit relations between the coefficients in diffusive-viscous wave equations and the measurable parameters for the tight oil/gas media. The derived equations are validated by numerical results. Based on the derived equations, reflection and transmission properties for a single tight interlayer are investigated. The numerical results demonstrate that the reflection and transmission of the seismic waves are affected by the thickness and attenuation of the interlayer, which is of great significance for the exploration of oil and gas. 相似文献
5.
6.
This paper analytically examines the horizontal vibration of a rigid disk on a saturated poroelastic half-space. The pressure-solid displacement form of the harmonic equations of motion for asymmetric dynamic problem are developed from the form of the equations originally presented by Biot. Making use of a new method the solution of the above equations is obtained. According to the mixed boundary -value conditions, the dual integral equations of the horizontal vibration of a rigid disk on a saturated poroelastic half-space are established. By appropriate transforms, it is shown that the dual integral equations can be reduced to a pair of Fredholm integral equations of the second kind, whose solutions are then computed. Numerical results for the horizontal dynamic compliance coefficient are given at the end of this paper. 相似文献
7.
The specific energy and specific force equations have many applications in open-channel flow problems. At present, these equations have analytical solution only for rectangular channels. Trial and error procedure also graphical solutions are the existing methods of solving these equations. No analytical solutions are available in the technical literature for these equations in trapezoidal and triangular channels because it is presumed that these equations are quintic equations. The inversion of such equations consists of finding the roots of quintic equations. In the current study for a given channel geometry and discharge, the subcritical (supercritical) depth is analytically found in terms of the other supercritical (subcritical) depth. For this purpose, by considering physically meaningful domains, a quintic equation has been reduced to a quartic equation. In the next step, this quartic equation has been converted to a resolvent cubic equation and two quadratic equations. This research shows these steps clearly to reach an acceptable physical analytic solution for water depth in trapezoidal and triangular channels. 相似文献
8.
用Biot介质参数说明了波达波夫震电波方程组中弹性动力学 参数的含义，解释了第一类和第二类震电效应的意义，在忽略第一类震电效应条件下将该方 程组与Pride方程组进行比较，说明了二者在描述第二类震电效应方面的异同点. 同时指出 ：波达波夫方程组忽略了流体与固体的耦合质量；方程中的黏性耗散项丢掉了一个孔隙度因 子，依据该方程组计算出的弹性波和转换电场的幅度将偏大；边界条件之一存在错误，会影 响对波在界面上的反射透射规律的描述. 相似文献
9.
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equations. In this study, we compare two kinds of such wave equations: the first-order (velocity–stress) and the second-order (displacement–stress) separate elastic wave equations, with the first-order (velocity–stress) and the second-order (displacement–stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-difference method. Comparisons are given of wavefield snapshots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corresponding first-order or second-order full elastic wave equations. These mixed equations are computationally slightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-component processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements. 相似文献
10.
Heterogeneous wave equations are more complicated numerically than homogeneous wave equations, but are necessary for physical validity. A wide variety of numerical solutions of seismic wave equations is available, but most produce strong numerical artefacts and local instabilities where model parameters change rapidly. Accuracy and stability of heterogeneous equations is achieved through staggered-grid formulations. A new pseudospectral staggered-grid algorithm is developed for the poroelastic (Biot) equations. The algorithm may be reduced to handle the elastic and acoustic limits of the Biot equations. Comparisons of results from poroelastic, elastic, acoustic and scalar computations for a 2D model show that porous medium parameters may affect amplitudes significantly. The use of homogeneous wave equations for modelling of a heterogeneous medium, or of a centred rather than a staggered grid, or of simplified (e.g. acoustic) wave equations when elastic or poroelastic media are synthesized, may produce erroneous or ambiguous interpretations. 相似文献
11.
Backus and Crampin derived analytical equations for estimating approximate phase-velocity variations in symmetry planes in weakly anisotropic media, where the coefficients of the equations are linear combinations of the elastic constants. We examine the application of similar equations to group-velocity variations in off-symmetry planes, where the coefficients of the equations are derived numerically. We estimate the accuracy of these equations over a range of anisotropic materials with transverse isotropy with both vertical and horizontal symmetry axes, and with combinations of transverse isotropy yielding orthorhombic symmetry. These modified equations are good approximations for up to 17% shear-wave anisotropy for propagations in symmetry planes for all waves in all symmetry systems examined, but are valid only for lower shear-wave anisotropy (up to 11%) in off-symmetry planes. We also obtain analytical moveout equations for the reflection of qP-, qSH-, and qSV- waves from a single interface for off-symmetry planes in anisotropic symmetry. The moveout equation consists of two terms: a hyperbolic moveout and a residual moveout, where the residual moveout is proportional to the degree of anisotropy and the spread length of the acquisition geometry. Numerical moveout curves are computed for a range of anisotropic materials to verify the analytical moveout equations. 相似文献
12.
13.
This paper considers the vertical dynamic response of a disk on a saturated poroelastic half-space. Firstly the pressure-solid displacement form of the harmonic equations of motion for a poroelastic solid are developed from the form of the equations originally presented by Biot. These equations are solved by a new method. Then the mixed boundary value problem for the vertical harmonic vibration of a disk on a poroelastic half-space is studied. The two types of drainage conditions at the surface of the poroelastic half-space are considered: (a) the surface of the poroelastic half-space is assumed to be completely pervious both within and exterior to the plate; (b) The interface between the plate and the poroelastic half-space is assumed to be impervious and the exterior region is assumed to be pervious. By using the Hankel transform techniques, the paper develops the governing dual integral equations. These governing integral equations are further reduced to systems of standard Fredholm integral equations of the second kind by Abel transform. 相似文献
14.
本文针对声-弹耦合介质,为尽可能的减少频率域正演模拟的计算内存,提高计算效率,在一阶非均质位移-应力波动方程的基础上,借助等效交错网格思想并充分考虑密度参数空间变化对地震波传播的影响,推导了声-弹耦合地震波波动方程.在流相介质和固相介质中分别采用非均质情况频率域二阶声压标量波、二阶纯位移控制方程,为保证流、固相介质间地震波能量的稳定传输和有效交换,提出了声-弹耦合界面转换过渡层方法,并详细阐述了过渡层与上下介质空间差分具体耦合方法.在与非均质纯位移波动方程正演结果对比分析的基础上,首先采用各向同性单层流相介质模型进行正演模拟验证了声-弹耦合方程数值模拟中过渡层策略的有效性和准确性,随后又数值模拟了地震波在声-弹耦合介质简单模型和复杂Marmousi2模型中的传播,验证了本文方法稳定性和准确性,同时该方法可以简单的推广到三维情况. 相似文献
15.
《Advances in water resources》1999,22(5):461-478
Fractional flow formulations of the multi-phase flow equations exhibit several attractive attributes for numerical simulations. The governing equations are a saturation equation having an advection diffusion form, for which characteristic methods are suited, and a global pressure equation whose form is elliptic. The fractional flow approach to the governing equations is compared with other approaches and the implication of equation form for numerical methods discussed. The fractional flow equations are solved with a modified method of characteristics for the saturation equation and a finite element method for the pressure equation. An iterative algorithm for determination of the general boundary conditions is implemented. Comparisons are made with a numerical method based on the two-pressure formulation of the governing equations. While the fractional flow approach is attractive for model problems, the performance of numerical methods based on these equations is relatively poor when the method is applied to general boundary conditions. We expect similar difficulties with the fractional flow approach for more general problems involving heterogenous material properties and multiple spatial dimensions. 相似文献
16.
17.
The accuracy of the new equations for long frame-type structures, derived by Kerr and Zarembski, and extended recently by Kerr and Accorsi also to dynamic analyses, is investigated. At first, the natural frequencies of a tall frame-type structure, with twenty floors, are determined using the new equations, and then they are compared with the corresponding results calculated using the finite element method and the common shear-building analysis. The found agreement with the finite element analysis is close. To check further the accuracy of the new equations, small scale models are tested on a shake table. The natural frequencies recorded in the tests are then compared with those determined from the new equations, for a wide range of geometrical parameters. The agreement with the test data is good. The presented study indicates that the new equations are well suited for analyzing the dynamic response of tall frame-type structures. 相似文献
18.
This study is aimed at developing statistical equations to estimate the inelastic displacement ratio of singledegree- of-freedom systems subjected to far fault repeated earthquakes. In the study, peak ground motion parameters are used to define the scatter of the original data. The ratio of peak ground acceleration to peak ground velocity, and peak ground velocity of the ground motion records and structural parameters such as period of vibration and lateral strength ratio are used in the proposed equations. For the development of the equations, nonlinear time history analyses of single-degree-offreedom systems are conducted. Then, the results are used in a multivariate regression procedure. The equations are verified by comparing the estimated results with the calculated results. The average error and coefficient of variation of the proposed equations are presented. The analyses results revealed that the direct use of peak ground motion parameters for the estimation of inelastic displacement ratio significantly reduced the scatter in the original data and yielded accurate results. From the comparative results it is also observed that results obtained using equations specific to peak ground velocity or peak ground acceleration to peak ground velocity ratio are similar. 相似文献
19.
20.
Jan Hrabě 《Studia Geophysica et Geodaetica》1994,38(2):157-167
Summary Principles of an alternative approach to the ray theory and paraxial ray approximation (PRA) theory are discussed. Invariant
equations for the ray, the eikonal equation and basic equations of the PRA-theory are derived on the basis of Riemannian geometry.
Paraxial rays and paraxial time field equations in general curvilinear coordinates are shown as an example of application. 相似文献