共查询到20条相似文献，搜索用时 31 毫秒
1.
近年来围绕四川盆地侏罗系陆相致密砂岩已取得了勘探突破,其中川中—川西过渡带具备形成大气田的地质条件,但对该套致密砂岩弹性性质变化规律的研究还较少,致使利用地震方法进行"甜点"储层预测的精度不高.本文利用四川盆地侏罗系沙溪庙组32块样品开展了系统的声学测量,在此基础上,分析了样品弹性性质的变化规律.结合X射线衍射矿物组分分析、扫描电镜、铸体薄片和岩石薄片特征确定了不同成岩作用对岩石储集性能的影响.研究结果表明,研究区致密砂岩储层表现为孔隙型储层,受差异性成岩作用影响,黏土含量、钙质含量和硅质含量的差异以及它们分布特征之间的差异对岩石弹性性质造成了很大的影响.在研究区对岩石物性及弹性性质有明显影响的成岩作用包括早期的钙质胶结作用、压实作用和溶蚀作用,因此针对不同时期的成岩作用对岩石弹性及物性的影响,利用接触-胶结模型、微分等效模量模型和临界孔隙度校正的Hashin-Shtrikman上限模型建立了研究区致密砂岩的岩石物理模型. 相似文献
2.
在页岩油气藏的开发和勘探阶段,需要对储层进行水力压裂改造,形成有利于油气聚集和运移的裂缝.地应力是进行水力压裂改造的重要参数,能够决定裂缝的大小、方向以及分布形态,影响着压裂的增产效果,且最大和最小水平应力差异比（ODHSR,Orthorhombic Differential Horizontal Stress Ratio）是评价储层是否可压裂成网的重要因子.本文探讨了基于地震数据估算地应力的方法,以指导页岩气的水力压裂开发.首先,利用叠前方位地震数据反演得到地层的弹性参数和各向异性参数;其次,基于正交各向异性水平应力差异比近似公式,利用反演得到的弹性参数和各向异性参数估算地层的ODHSR;最后,选取某工区的裂缝型页岩储层的叠前方位地震数据对该方法进行实际应用.实际工区地震数据应用表明,基于叠前方位地震数据反演得到的ODHSR能够有效的识别储层中易于压裂成网的区域. 相似文献
3.
Akinori Hashima Yukitoshi Fukahata Chihiro Hashimoto Mitsuhiro Matsu’ura 《Pure and Applied Geophysics》2014,171(8):1669-1693
We derived explicit expressions in the time domain for 3-D quasi-static strain and stress fields, due to a point moment tensor source in an elastic surface layer overlying viscoelastic half-space under gravity. The expressions of strain in the elastic surface layer were directly obtained from the expressions of displacement in our previous paper. The conversion of strain into stress is easy, because the stress–strain relation of elastic material is linear. In the viscoelastic substratum, the expressions of strain were obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. The strain–stress conversion is not straightforward, as the stress–strain relation of viscoelastic material is usually given in a differential form. To convert strain into stress, we used an integral form of the stress–strain relation instead of the usual differential form. The expressions give the responses of elastic half-space at \( t = 0 \) , and the responses of an elastic plate floating on non-viscous liquid at \( t = \infty \) . The moment tensor is rationally decomposed into the three independent force systems, corresponding to isotropic expansion, shear faulting and crack opening, and so the expressions include the strain and stress fields for these force systems as special cases. As the first numerical example, we computed the temporal changes in strain and stress fields after the sudden opening of an infinitely long vertical crack cutting the elastic surface layer. Here, we observe that the stress changes caused by the sudden crack opening gradually decay with time and vanish at \( t = \infty \) everywhere. After the completion of stress relaxation, a characteristic pattern of shear strain remains in the viscoelastic substratum. Since the strain and stress fields at \( t = \infty \) can be read as the strain- and stress-rate fields caused by steady crack opening, respectively, this numerical example demonstrates the realization of a steady stress state supported by steady viscous flow in the asthenosphere, associated with steady seafloor spreading at mid-ocean ridges. For the second numerical example, we computed the temporal changes in strain and stress fields after the 2011 Tohoku-oki mega-thrust earthquake, which occurred at the North American-Pacific plate interface. In this numerical example, the stress changes caused by coseismic fault slip vanish at \( t = \infty \) in the viscoelastic substratum, but remain in the elastic surface layer. The coseismic stress changes (and also strain changes) in the elastic surface layer diffuse away from the source region with time, due to gradual stress relaxation in the viscoelastic substratum. 相似文献
4.
采用规则网格有限差分方法对二维平面弹性波动方程进行差分离散,得到相应的弹性波动方程的有限差分方程,再将弹性波动方程的差分格式与吸收边界、自由边界的离散形式结合形成弹性波动方程有限差分方程解决问题的主体,将其应用于含方形凹陷半无限非均匀介质的模型中进行数值模拟,得到此离散化模型中不同时刻不同节点的位移值。针对具体算例,运用上述方法结合科学计算软件MATLAB和结果后处理软件DIFEM ISOLINE PLOTER得到不同时刻的水平方向位移等值线图与接收器测量点处的合成位移记录,讨论非均匀介质、吸收边界、方形凹陷等对波动特性的影响。 相似文献
5.
This paper presents a time-dependent semi-analytical artificial boundary for numerically simulating elastic wave propagation problems in a two-dimensional homogeneous half space. A polygonal boundary is considered in the half space to truncate the semi-infinite domain, with an appropriate boundary condition imposed. Using the concept of the scaled boundary finite element method, the wave equation of the truncated semi-infinite domain is represented by the partial differential equation of non-constant coefficients. The resulting partial differential equation has only one spatial coordinate variable and time variable. Through introducing a few auxiliary functions at the truncated boundary, the resulting partial differential equations are further transformed into linear time-dependent equations. This allows an artificial boundary to be derived from the time-dependent equations. The proposed artificial boundary is local in time, global at the truncated boundary and semi-analytical in the finite element sense. Compared with the scaled boundary finite element method, the main advantage in using the proposed artificial boundary is that the requirement for solving a matrix form of Lyapunov equation to obtain the unit-impulse response matrix is avoided, so that computer efforts are significantly reduced. The related numerical results from some typical examples have demonstrated that the proposed artificial boundary is of high accuracy in dealing with time-dependent elastic wave propagation in two-dimensional homogeneous semi-infinite domains. 相似文献
6.
Non-linear buffeting response analysis of long-span suspension bridges with central buckle 总被引：1，自引：1，他引：0
The rigid central buckle employed in the Runyang Suspension Bridge (RSB) was the first time it was used in a suspension bridge in China. By using a spectral representation method and FFT technique combined with measured data,a 3D fluctuating wind field considering the tower wind effect is simulated. A novel FE model for buffeting analysis is then presented,in which a specific user-defined Matrix27 element in ANSYS is employed to simulate the aeroelastic forces and its stiffness or damping matrices are parameterized by wind velocity and vibration frequency. A nonlinear time history analysis is carried out to study the influence of the rigid central buckle on the wind-induced buffeting response of a long-span suspension bridge. The results can be used as a reference for wind resistance design of long-span suspension bridges with a rigid central buckle in the future. 相似文献
7.
Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation. 相似文献
8.
Thermo-elastic interactions in an infinite elastic solid due to a concentrated transient heat source
Snehanshu Kumar Roy Choudhuri 《Pure and Applied Geophysics》1972,97(1):57-69
Summary The present paper is concerned with the determination of thermo-elastic stress and temperature distribution in an infinite elastic solid when it is subjected to a concentrated transient heat source taking into account the effect of coupling. The fundamental partial differential equation for the thermo-elastic potential is solved by means of the operational method. The solutions are obtained for small values of time in terms of known functions such as complementary error functions and associated complementary error functions and have the form of power series with respect to the parameter of coupling . Since this parameter is generally very small for most of the metals, the expressions for stress and temperature distribution have been obtained by considering terms upto linear . 相似文献
9.
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. 相似文献
10.
A nearly analytic exponential time difference method for solving 2D seismic wave equations 
下载免费PDF全文
下载免费PDF全文 In this paper, we propose a nearly analytic exponential time difference (NETD) method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approximate the high-order spatial differential operators and transform the seismic wave equations into semi-discrete ordinary differential equations (ODEs). Then, the converted ODE system is solved by the exponential time difference (ETD) method. We investigate the properties of NETD in detail, including the stability condition for 1-D and 2-D cases, the theoretical and relative errors, the numerical dispersion relation for the 2-D acoustic case, and the computational efficiency. In order to further validate the method, we apply it to simulating acoustic/elastic wave propagation in multilayer models which have strong contrasts and complex heterogeneous media, e.g., the SEG model and the Marmousi model. From our theoretical analyses and numerical results, the NETD can suppress numerical dispersion effectively by using the displacement and gradient to approximate the high-order spatial derivatives. In addition, because NETD is based on the structure of the Lie group method which preserves the quantitative properties of differential equations, it can achieve more accurate results than the classical methods. 相似文献
11.
Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix 下载免费PDF全文
下载免费PDF全文 To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems. 相似文献
12.
V. M. Pavlov 《Izvestiya Physics of the Solid Earth》2013,49(1):24-33
A new algorithm is proposed for calculating the complete synthetic seismograms from a point source in the form of the sum of a single force and a dipole with an arbitrary seismic moment tensor in a plane layered medium composed of homogenous elastic isotropic layers. Following the idea of (Alekseev and Mikhailenko, 1978), an artificial cylindrical boundary is introduced, on which the boundary conditions are specified. For this modified problem, the exact solution (in terms of the displacements and stresses on the horizontal plane areal element) in the frequency domain is derived and substantiated. The unknown depth-dependent coefficients form the motion-stress vector, whose components satisfy the known system of ordinary differential equations. This system is solved by the method that involves the matrix impedance and propagator for the vector of motion, as previously suggested by the author in (Pavlov, 2009). In relation to the initial problem, the reflections from the artificial boundary are noise, which, to a certain degree, can be suppressed by selecting a long enough distance to this boundary and owing to the presence of a purely imaginary addition to the frequency. The algorithm is not constrained by the thickness of the layers, is applicable for any frequency range, and is suitable for computing the static offset. 相似文献
13.
The objective of this research is to study the dynamic response characteristics of a three-beam system with intermediate elastic connections under a moving load/mass-spring. In this study, the finite Sine-Fourier transform was performed for the dynamic partial differential equations of a simply supported three-beam system (SSTBS) under a moving load and a moving mass-spring, respectively. The dynamic partial differential equations were transformed into dynamic ordinary differential equations relative to the time coordinates, and the equations were solved and the displacement Fourier amplitude spectral expressions were obtained. Finally, based on finite Sine-Fourier inverse transform, the expressions for dynamic response of SSTBS under the moving load and moving mass-spring were obtained. The proposed method, along with ANSYS, was used to calculate the dynamic response of the SSTBS under a moving load/mass-spring at different speeds. The results obtained herein were consistent with the ANSYS numerical calculation results, verifying the accuracy of the proposed method. The influence of the load/mass-spring’s moving speed on the dynamic deflections of SSTBS were analyzed. SSTBS has several critical speeds under a moving load/mass-spring. The vertical acceleration incurred by a change in the vertical speed of SSTBS due to the movement of mass-spring and the centrifugal acceleration produced by the movement of massspring on the vertical curve generated by SSTBS vibration could not be neglected. 相似文献
14.
Deviatoric stress in a diamond anvil cell with gold as a pressure and stress indicator is measured by two complementary techiques using synchrotron radiation. The first method employs a white X-ray beam using energy disperisive X-ray diffraction. The incident X-ray beam is parallel to the load axis and the diffraction pattern is recorded at a low two-theta angle. Using powder diffraction patterns of polycrystallin gold, we measured the elastic strain of two crystal planes oriented normal to the diffraction vector. Stresses nearly parallel and perpendicular to the load axis can be calculated by stress-strain tensor relatonship. The other method uses a monochromatic wiggler X-ray beam. In this case, the diamond cell is oriented so that the incident beam is perpendicular to the load axis. The diffraction pattern is recorded on an image plate area detector. Elastic strains responding to stresses perpendicular and parallel to the load axis can be measured and stresses of the same orientations can be calculated from the strain data. These measurements provide a lower bound of the actual differential stress in a diamond cell. With these techniques, we can measure stress distribution in a less deviatoric gasketted sample and determine yield strength of mantle materials at high pressures and temperatures. 相似文献
15.
A general solution is deduced of the differential equations describing the propagation of elastic waves in a dissipative liquid-filled viscoelastic porous solid. The velocities of three existing waves have been expressed in convenient form using the moduli of the solid phase and by introducing the frequency-dependent equivalent mass densities. The solution is then used to examine some of the phenomena which arise when each of the three-body waves, in turn, are incident on a traction-free plane boundary. Analytic expressions for the reflection coefficients are obtained. Numerical calculations have been made, for a particular model, in case of incidentP
I wave. Effect of viscoelasticity and viscosity on the reflection coefficients has also been exhibited. 相似文献
16.
用布希涅斯克定义的弹性半空间内的垂直位移包括两项积分,除了积分号前面系数的差别之外,第一项积分是单层位势而第二项积分为双层位势。若扁壳基础是正高斯曲率的几何曲面,则壳底与半空间表面间的挤压强度就是半空间表面作用的分布垂直荷载。当越过边界时,双层势位的函数值和单层势位的法向导数值发生跳跃。利用这些性质,本文得出布希涅斯克积分的反演公式,从而避开要求解偏微分—积分方程组的巨大数学困难而易于得出解析解。以椭园抛物面扁壳为例说明本文方法的应用。 相似文献
17.
微分求积法（DQM）是1种求解微分方程初（边）值问题的数值方法,通常以较小的计算工作量即可获得较高的数值精度。这种方法应用于工程领域时多用来解决梁、板等结构的静力分析或结构特征值分析等问题,即对边值问题的微分方程的求解。结构动力分析属于初值问题,荷载和结构反应都具有特殊性,直接套用DQM求解边值问题并不能获得问题的解。本文尝试利用微分求积原理建立求解结构动力反应的具体方法。借鉴单元法的思想,将荷载持时划分为若干个时步,在每个时步内对动态荷载和结构反应进行离散,然后用DQM对时步逐个进行求解,得到体系在整个时域内的反应过程。通过对3种不同自振周期的线弹性单自由度体系在不同频率简谐激励下反应的计算,阐释了本文方法的可行性以及高精度、高效率的特点,通过数值试验确定了时步内相对较优的节点数,并为时步长度的选取提供了建议。 相似文献
18.
To resolve the issue regarding inaccurate prediction of the hysteretic behavior by micro-based numerical analysis for partially-restrained(PR)steel frames with solid reinforced concrete(RC)infill walls,an innovative simplified model of composite compression struts is proposed on the basis of experimental observation on the cracking distribution,load transferring mechanism,and failure modes of RC infill walls filled in PR steel frame.The proposed composite compression struts model for the solid RC infill walls is composed ofαinclined struts and main diagonal struts.Theαinclined struts are used to reflect the part of the lateral force resisted by shear connectors along the frame-wall interface,while the main diagonal struts are introduced to take into account the rest of the lateral force transferred along the diagonal direction due to the complicated interaction between the steel frame and RC infill walls.This study derives appropriate formulas for the effective widths of theαinclined strut and main diagonal strut,respectively.An example of PR steel frame with RC infill walls simulating simulated by the composite inclined compression struts model is illustrated.The maximum lateral strength and the hysteresis curve shape obtained from the proposed composite strut model are in good agreement with those from the test results,and the backbone curve of a PR steel frame with RC infill walls can be predicted precisely when the inter-story drift is within 1%.This simplified model can also predict the structural stiffness and the equivalent viscous damping ratio well when the inter-story drift ratio exceeds 0.5%. 相似文献
19.
完全匹配层(PML)作为一种稳定高效的吸收边界条件,广泛应用于基于一阶电磁波动方程的探地雷达(GPR)数值模拟中.为解决基于二阶电磁波动方程的GPR数值模拟的吸收边界问题,本文借鉴二阶弹性波动方程的PML边界条件构建思想,提出了一种适合二阶电磁波动方程GPR时域有限元模拟的PML边界条件.从二阶电磁波动方程出发,基于复拉伸坐标变换,推导了PML算法的频域表达式;通过合理构造辅助微分方程,得到了PML算法的时域表达式,并以变分形式(弱形式)加载到GPR时域有限元方程中,实现了PML边界条件在二阶电磁波动方程GPR时域有限元模拟中的应用.在此基础上,对比了无边界条件、Sarma边界条件和PML边界条件下均匀模型的波场快照、单道波形、时域反射误差和能量衰减曲线,结果表明:PML边界条件的吸收效果要远优于Sarma边界条件,具有近似零反射系数.一个复杂介质模型的正演模拟验证了PML边界条件在非均匀地电结构中电磁波传播模拟的良好吸收效果. 相似文献
20.
本以两承压岩体相互作用系统作为“弹性回跳”地震模型，用损伤力学的方法研究了系统失稳即地震前位移速率、声发射等前兆规律与岩石破裂过程的关系，用差分法在计算机上求得了地震发生过程中的变形与地声序列曲线，研究发现在某些情况下出现地声相对平静是由于两岩体介质中某一介质体破裂而另一个介质体弹性恢复所致。 相似文献