工程结构三维疲劳裂纹最大应力强度因子计算

介绍了裂纹的类型、裂纹尖端应力场的奇异性。以一维问题为例,推导论证了奇异单元能够很好的反映裂纹尖端应力场的奇异性。应力强度因子一般表达式表明应力强度因子与载荷呈线性关系,并依赖于物体和裂纹的几何形状和尺寸。本文借助大型通用有限元软件ANSYS,采用位移外插法计算了三维表面裂纹前沿不同位置处的应力强度因子,并与《应力强度因子手册》基于实验的理论公式计算结果相比较。结果表明:有限元结果与理论解误差较小,裂纹最深处应力强度因子最大。     

Influence of electromagnetic fields on the seismotectonic strain rate; relaxation and active monitoring of elastic stresses

We study the changes in the rate of seismotectonic strains in the Earth’s crust in two of the most seismically active Central Asian areas, the Garm region in Tadjikistan and the Northern Tien Shan, which are excited by high energy pulsed electromagnetic irradiation. It is shown that in both areas the average rates of strain release increase by a factor of 10–20. The increment of the seismotectonic strains initiated during the entire period of electromagnetic treatment is 1.3 × 10⁻⁶–6.3 × 10⁻⁶, which is comparable with the value of crustal strains observed during the preparation of strong earthquakes and amounts to 3–13% of the ultimate strains for rocks. Such a contribution to the process of quasi-plastic strains of the crust results in the additional relaxation of elastic stresses by 0.4–1.67 bar, which corresponds to approximately 1–17% of the stress released in the sources of strong earthquakes. The spatial variations in the changes of the strain rate during electromagnetic treatment are studied, and their correlation with the level of stresses in the Earth’s crust is revealed. It is shown that the change in the strain rate during the electromagnetic excitation can be used as the criterion for active monitoring of the stress state of the medium for the purposes of predicting the location, time, and magnitude of strong earthquakes.     

Characterizing the effect of elastic interactions on the effective elastic properties of porous,cracked rocks

Elastic interactions between pores and cracks reflect how they are organized or spatially distributed in porous rocks. The principle goal of this paper is to understand and characterize the effect of elastic interactions on the effective elastic properties. We perform finite element modelling to quantitatively study how the spatial arrangement of inclusions affects stress distribution and the resulting overall elasticity. It is found that the stress field can be significantly altered by elastic interactions. Compared with a non‐interacting situation, stress shielding considerably stiffens the effective media, while stress amplification appreciably reduces the effective elasticity. We also demonstrate that the T‐matrix approach, which takes into account the ellipsoid distribution of pores or cracks, can successfully characterize the competing effects between stress shielding and stress amplification. Numerical results suggest that, when the concentrations of cracks increase beyond the dilute limit, the single parameter crack density is not sufficient to characterize the contribution of the cracks to the effective elasticity. In order to obtain more reliable and accurate predictions for the effective elastic responses and seismic anisotropies, the spatial distribution of pores and cracks should be included. Additionally, such elastic interaction effects are also dependent on both the pore shapes and the fluid infill.     

Dispersion and attenuation of elastic waves due to multiple scattering from inclusions

The propagation of elastic waves in a medium containing many inclusions is considered. Under the assumption that the spatial distribution of inclusions is uniform, a general equation is derived for the determination of the velocity dispersion and attenuation coefficient of the effective waves. A simple example is presented where scatterers are infinitesimally thin cracks. The calculated results show that the attenuation coefficient Q^?1 takes a peak value for the wavelength nearly equal to twice the crack length.     

Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory

Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given.     

Analytic expressions for the velocity sensitivity to the elastic moduli for the most general anisotropic media

For non‐linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density‐normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS₁, qS₂). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP‐, qS₁‐ and qS₂‐wave computations and mitigates the singularity problem for the two quasi‐shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes.     

Boundary element simulation of scattering of elastic waves by 3-D cracks

Numerical modelling techniques are now becoming common for understanding the complicated nature of seismic wave propagation in fractured rock. Here the Indirect Boundary Element Method (IBEM) is applied to study scattering of elastic waves by cracks. The problem addressed in this paper is the diffraction of P and S waves by open 3-D cracks of arbitrary shape embedded in a homogeneous isotropic medium. The IBEM yields the value of the jump of displacements between opposite surfaces of the crack, often called Crack Opening Displacement (COD). This is used to evaluate the solution away from the crack. We use a multi-regional approach which consists of splitting a surface S into two identical surfaces S⁺ and S⁻ chosen such that the crack lies at the interface. The resulting integral equations are not hyper-singular and wave propagation within media that contain open cracks can be rigorously solved. In order to validate the method, we compare results of displacements of a penny-shaped crack for a vertical incident P-wave with the classic results by Mal (1970) obtaining excellent agreement. This comparison gives us confidence to study cases where no analytic solutions exist. Some examples of incidence of P or S waves upon cracks with various shapes are depicted and the salient aspects of the method are also discussed. Both frequency and time-domain results are included.     

Effect of stress interactions on anisotropic P-SV-wave dispersion and attenuation for closely spaced cracks in saturated porous media

Generally, local stress induced by individual crack hardly disturbs their neighbours for small crack densities, which, however, could not be neglected as the crack density increases. The disturbance becomes rather complex in saturated porous rocks due to the wave-induced diffusion of fluid pressures. The problem is addressed in this study by the comparison of two solutions: the analytical solution without stress interactions and the numerical method with stress interactions. The resultant difference of effective properties can be used to estimate the effect of stress interactions quantitatively. Numerical experiments demonstrate that the spatial distribution pattern of cracks strongly affects stress interactions. For regularly distributed cracks, the resulting stress interaction (shielding or amplification) shows strong anisotropy, depending on the arrangement and density of cracks. It has an important role in the estimation of effective anisotropic parameters as well as the incident-angle-dependency of P- and SV-wave velocities. Contrarily, randomly distributed cracks with a relative small crack density generally lead to a strong cancellation of stress interactions across cracks, where both the numerical and analytical solutions show a good agreement for the estimation of effective parameters. However, for a higher crack density, the incomplete cancellation of stress interactions is expected, exhibiting an incidence-angle dependency, slightly affecting effective parameters, and differentiating the numerical and analytical solutions.     

Icequakes and glacier motion: The hans glacier,spitsbergen

The recording of glacier icequakes has been carried out on the Hans glacier, Spitsbergen. The icequakes have been connected with the formation of open cracks. The energy of the main group of icequakes has ranged from 0.01 J to 6 J, their seismic moment being from 0.5 · 10¹⁴ dyne · cm to 26 · 10¹⁴ dyne · cm. The nonelastic volume increase corresponding to each crack has ranged from 0.2 · 10⁻² m³ to 8.8 · 10⁻² m³. The stresses released due to extensional faulting ranged from 25N/m² to 500 N/m²; these values are so small that they suggest an incomplete release of the medium. The contribution into the glacier flow rate of the deformation in the surface glacier layer due to extensional faulting is lower by 2 to 3 orders of magnitude than the total plastic deformation associated with the ice creeping.     

基于方位地震数据的地应力反演方法

在页岩油气藏的开发和勘探阶段,需要对储层进行水力压裂改造,形成有利于油气聚集和运移的裂缝.地应力是进行水力压裂改造的重要参数,能够决定裂缝的大小、方向以及分布形态,影响着压裂的增产效果,且最大和最小水平应力差异比（ODHSR,Orthorhombic Differential Horizontal Stress Ratio）是评价储层是否可压裂成网的重要因子.本文探讨了基于地震数据估算地应力的方法,以指导页岩气的水力压裂开发.首先,利用叠前方位地震数据反演得到地层的弹性参数和各向异性参数;其次,基于正交各向异性水平应力差异比近似公式,利用反演得到的弹性参数和各向异性参数估算地层的ODHSR;最后,选取某工区的裂缝型页岩储层的叠前方位地震数据对该方法进行实际应用.实际工区地震数据应用表明,基于叠前方位地震数据反演得到的ODHSR能够有效的识别储层中易于压裂成网的区域.     

Predicted near-field ground motion for dynamic stress-drop models

We propose a finite difference method, using a hexagonal grid, to compute displacements (stresses, velocities, accelerations) in the near-field of a 2-D in-plane stress-drop crack, in both whole space (constant stress-drop) and half-space (depth-dependent stress-drop). To exercise the method, the stress field distribution is evaluated for both fundamental 2-D shear cracks, anti-plane. In order to test the method's reliability, the results are compared with some analytical and numerical solutions available in the literature (Kostrov, 1964;Virieux andMadariaga, 1982). For the in-plane source, the results emphasize that the method can resolve the stress concentration due to the rupture front from the stress peak associated with the shear wave propagating in front of the crack. Synthetic motions are computed on the fault, but also in an infinite medium and at the free surface. The rather complex waveforms generated in the near-field, even by simple sources, emphasize the contribution of all wave terms (near, intermediate and far-field) to the motion. The presence of near-field and the numerical procedure explain the significant low frequency content of the computed seismograms. The set of treated problems proves the method is stable and accurate.     

Effective elastic properties of randomly fractured soils: 3D numerical experiments

This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results.     

Dynamic stresses in an elastic half-space

This paper deals with the problem of time-varying point loads applied onto the surface of an elastic half-space and the stresses that such loads elicit within that medium. The emphasis is on the evaluation of the isobaric contours for all six of the stress components at various frequencies of engineering interest and for a full range of Poisson’s ratios. The extensive set of pressure bulbs presented herein may be of help in predicting the severity of dynamic effects in common practical situations in engineering—or even the lack thereof.     

Changes in geophysical properties caused by fluid injection into porous rocks: analytical models

Analytical models are provided that describe how the elastic compliance, electrical conductivity, and fluid‐flow permeability of rocks depend on stress and fluid pressure. In order to explain published laboratory data on how seismic velocities and electrical conductivity vary in sandstones and granites, the models require a population of cracks to be present in a possibly porous host phase. The central objective is to obtain a consistent mean‐field analytical model that shows how each modeled rock property depends on the nature of the crack population. The crack populations are described by a crack density, a probability distribution for the crack apertures and radii, and the averaged orientation of the cracks. The possibly anisotropic nature of the elasticity, conductivity, and permeability tensors is allowed for; however, only the isotropic limit is used when comparing to laboratory data. For the transport properties of conductivity and permeability, the percolation effect of the crack population linking up to form a connected path across a sample is modeled. However, this effect is important only in crystalline rock where the host phase has very small conductivity and permeability. In general, the importance of the crack population to the transport properties increases as the host phase becomes less conductive and less permeable.     

分析地震波估算地壳内的应力值

本文讨论了利用破裂力学理论说明地震破裂的过程, 认为地震本质上是岩石在应力作用下的低应力破裂现象.它是岩石中的裂纹不断稳态扩展、最后进入失稳扩展的结果.分析了在扩展过程中应力和位移的变化, 发现任何将要破裂的那一点的应力都要由初始应力0升高到屈服应力y 以后才破裂, 破裂后裂纹面上的点的应力降到0.在破裂前和破裂后的位移, 都可由弹性力学方程给出.在破裂的一瞬间破裂的端点产生的非弹性位移, 则不能由弹性力学方程给出.它可以由断裂力学中的裂纹滑开位移公式近似给出.根据位错模式由于计算弹性波辐射场的位错量 D(, t), 正是破裂瞬间产生的非弹性位移, 所以用弹性位移公式来计算地震位错量是错误的.我们采用了裂纹滑开位移公式来计算地震位错量, 从而导出了较合理的计算地震释放总能量的公式 ET=yDS(y 为屈服强度;D为平均位错;S 为断层面积)以及估算初始应力值0的公式:0 =[Dmax/L4y/(1-) ]1/2(L 为断层长度).用它们计算了一些地震的 Er 和0, 分别列于表1和表2.这些结果比以往的结果要更合理一些。 结果表明:(1)地震多数是在低应力作用下(即低初始应力)发生的(约100——200巴);(2)地震释放的总能量约比地震波能量大一个数量级.      

裂缝型HTI介质的纵波弹性阻抗

裂缝发育导致的地震波场各向异性,是裂缝型油气藏的一个重要特征.然而针对该类油气藏的叠前弹性阻抗技术,一直没有通过各向异性弹性阻抗公式建立与裂缝参数的直接关系.本文在回顾封闭平行硬币状裂缝模型和线性滑动模型的基础上,将裂缝填充物性质、分布密度与围岩的横纵波速度比,引入裂缝型HTI介质弹性阻抗公式,并对相应的归一化弹性阻抗响应特征进行模拟分析.分析表明,弹性阻抗受以上三参数的影响规律存在明显差异,其中含气裂隙介质随裂缝密度的变化要明显高于对应的流体裂隙介质,流体填充裂隙介质随横、纵波速度比的变化幅度要高于对应的含气裂隙介.该结论初步为裂缝型油气藏流体识别工作提供了依据.     

Cohesive Crack Analysis of Toughness Increase Due to Confining Pressure

Apparent fracture toughness in Mode I of microcracking materials such as rocks under confining pressure is analyzed based on a cohesive crack model. In rocks, the apparent fracture toughness for crack propagation varies with the confining pressure. This study provides analytical solutions for the apparent fracture toughness using a cohesive crack model, which is a model for the fracture process zone. The problem analyzed in this study is a fluid-driven fracture of a two-dimensional crack with a cohesive zone under confining pressure. The size of the cohesive zone is assumed to be negligibly small in comparison to the crack length. The analyses are performed for two types of cohesive stress distribution, namely the constant cohesive stress (Dugdale model) and the linearly decreasing cohesive stress. Furthermore, the problem for a more general cohesive stress distribution is analyzed based on the fracture energy concept. The analytical solutions are confirmed by comparing them with the results of numerical computations performed using the body force method. The analytical solution suggests a substantial increase in the apparent fracture toughness due to increased confining pressures, even if the size of the fracture process zone is small.     

Surface waves in a micropolar elastic solid containing voids

The present paper investigates the effect of voids on the propagation of surface waves in a homogeneous micropolar elastic solid medium which contains a distribution of vacuous pores (voids). The general theory for surface wave propagation in micropolar elastic media containing voids has been presented. Particular cases of surface waves (Rayleigh’s, Love’s and Stoneley’s) in micropolar media which contain vacuous pores have been deduced from the above general theory. Discussions have been made in each case to highlight the effect of voids and micropolar character of the material medium separately. Their joint effect has also been studied in details. Modulation of Rayleigh wave velocity has been studied numerically. It is observed that Love waves are not affected by the presence of voids.     

有效应力变化对致密砂岩孔隙结构及弹性波响应的影响规律

致密砂岩气藏具有裂缝发育和有效应力高的特征,研究不同有效压力下孔、裂隙介质地震波传播特征,有利于地震解释与地下储层的识别.但是前人的研究较少考虑岩石内部微观孔隙结构特征与孔隙、裂隙间流体流动的关系.本文首先通过选取四川盆地典型致密砂岩岩样,在不同有效压力下对岩石样本进行超声波实验测量.然后基于实验测得的纵、横波速度进行裂隙参数反演,得到不同有效压力下致密砂岩样本的裂隙孔隙度.再将裂隙孔隙度和样本岩石物理参数代入双重孔隙介质模型,模拟得到不同有效压力下饱水致密砂岩样本纵横波速度频散和衰减的变化规律.结果表明模型预测的速度频散曲线与纵波速度实验测量结果能够较好的吻合.最后统计分析了致密砂岩裂隙参数,得到了致密砂岩储层裂隙参数随有效压力及孔隙度变化特征.依据实际岩石物理参数建立模型,其裂隙参数三维拟合结果能够较好描述致密砂岩裂隙结构与孔隙度、应力的关联,可为实际地震勘探中预测储层裂缝性质提供基础依据.     

A Network Approach to Fracture: The Effect of Heterogeneity and Loading Conditions

—Fracture in a heterogeneous solid is simulated on a triangular network of bonds which figure the potential cracks of the medium. Heterogeneity is introduced by assuming a statistical strength distribution for the bonds. External stresses are applied to the network and the evolution of the bond population is analyzed when the stresses are increased. Bond-breaking is controlled by the crack-exten sion force which takes into account crack interactions by using an iterative procedure. Crack propaga tion leads to the coalescence of broken bonds crack clusters are formed. By using this kind of approach which combines fracture mechanics and network modelling, we are able to simulate the rupture of a rock specimen under various loading conditions without heavy computation. We discuss physical properties of the rupture process by examining the rupture stress and the geometric properties of the macroscopic fracture and their dependence on loading conditions and heterogeneity. Analysis of the geometric characteristics shows that the number of broken bonds can be fitted by a power law of the lattice size, the exponent depending on the loading conditions. Furthermore, an approximate computation of the mechanical response of the network demonstrates that the threshold secant modulus may be a more legitimate choice for a damage parameter in terms of system size independence.