平面应变双轴压缩岩样剪应力异常及破坏过程模拟

利用FLAC,模拟了双轴压缩岩样的破坏过程及剪应力异常.在峰值强度之前及之后,岩石的本构模型分别取为线弹性及莫尔库仑剪破坏与拉破坏复合的应变软化模型.本文仅分析了3个具有较大应力降的周期的剪应力异常及塑性区分布.在第1及第2周期,剪破坏分别发生在试样的两侧及背景空区内部,前兆明显,地震属于断错型.在第3周期,背景空区之外的新剪切破裂带引起断错型地震;空区内部的老破裂带引起走滑型地震,未观测到前兆.在塑性区边缘,剪应力梯度较高.破裂带位置剪应力值反而较低.尽管在应力-应变曲线的软化阶段之前,不同单元的剪应力表现出千姿百态的复杂形态,但尚有规律可循.若几个单元都位于某条破裂带上,则剪应力的变化可能是同步的,或演变规律类似.鉴于原地复发地震类型可能不同,前兆将有差别.     

裂纹间动态相互作用的实验观测与理论分析──以共线剪切裂纹归并为例

本文以共线剪切裂纹的归并为例,研究裂纹系的动态相互作用．实验结果表明：在张破裂的意义上,小间距共线剪切裂纹的内端部实际上是阻碍结构,通过迂回弯曲的张破裂归并．本文提出一种估计张破裂扩展途径的近似方法,利用张应力破裂判据找出了潜在的破裂扩展途径的第一近似,较好地模拟了共线裂纹的破裂途径,和现有的数值模拟方法比较,可以大大缩短机时,简化计算程序,保证精确度．本研究结出了岩石剪切破裂“隧道效应”（即破裂跳过障碍体,而障碍体本身不破坏）的机理,还给出了脆性介质在摩擦、磨损过程中碎屑产生的机理．     

在围压条件下岩体模型应力应变关系的实验研究

含切口岩体模型在差应力作用下的应力应变关系表明:岩体不同部位的应力应变曲线按斜率变化特征可以划分为四个发展阶段。根据对切口岩石破裂过程的对比,各阶段的变形机制分别与岩体内原微裂纹闭合、孔隙和颗粒弹性变形、宏观张裂剪裂以及剪切、张-剪主破裂相对应,但各阶段的发展进程在时间、空间上却不一致,这一差异均与切口分布和岩性变化相关。最后讨论了应力应变特征与发震断层孕震过程的关系,对地震预报可能有实际意义     

The effect of orthorhombic anisotropy and its implication for oil recovery and reservoir exploitation

Results of an experimental study of shear‐ and compressional‐wave propagation in an orthorhombically anisotropic medium are presented. The experiments were performed on a physical model consisting of two sets of fractures. The first set consisted of orientated rubber inclusions simulating weak material‐filled cracks. The second set consisted of a system of closely spaced parallel fractures simulated by thin plates of epoxy resin, superimposed orthogonally on the first set. Three cases of fracture orientations within the model were identified and studied. Case 1 is analogous to a jointed fracture reservoir with one vertical set of fluid‐filled cracks or fractures and one non‐filled horizontal set. This case is referred to as JFV. Case 2 is analogous to a double fracture reservoir with one horizontal set of fluid‐filled fractures or cracks and one non‐filled vertical set. This is referred to as DFH. Case 3 is analogous to a double fracture reservoir with two vertical sets of fractures or cracks, with only one fluid‐filled. Case 3 is referred to as DFV. A pulse transmission method was performed on all three modelled cases along the three principal axes. A directional variation in the compressional‐ and shear‐wave velocities, as well as distinct shear‐wave splitting, was observed. The elastic constants for each case were determined and differences between them were noted and compared with the controlled results of both layered (transverse isotropy, TI) and vertically fractured (azimuthally anisotropic models, VF) media. The differences in elastic moduli and velocities indicate the potential of recognizing the different fracture orientations and suggest an approach to designing a method of drilling to further enhance oil recovery and reservoir exploitation.     

Percolation and fracture

A statistical model of fracture, based on percolation theory, is presented which allows the quantitative evaluation of clustering of cracks in solids. Unlike models of branching processes which are more in accord with Griffith fracture, the concept of percolation lattices (finite and infinite) is used following from a physical model of multiple fracture.Some experimental results on acoustic emission, dilatancy and geophysical precursors of earthquakes can be correlated with the percolation fracture model. The model does not depend on the mechanism of crack formation, critical parameters being the number of elementary events (cracks), the dimensionality of the process and the coordination number for a network of cracks and, in finite systems, their specific size.Fracture prediction is possible from the number of elementary acts, cluster statistics and other characteristic parameters of the model. Possible applications of the percolation model for earthquake prediction are considered.     

裂缝尺度对剪切破裂能的影响

采用大理岩圆柱试样,与轴线约成30角预制不连续裂缝,再用环氧树脂胶接,以实验观测裂缝尺度及其填充物对剪切破裂能的影响.结果表明:当裂缝相对面积从0增加到90%时,破裂能从(6.92.1)103焦耳/米2逐步降到(0.60.5)1O3焦耳/米2;而使用 Knopoff方法和引入滑动弱化模式后的方法,求得的理论应力降几乎不变.它们与实测应力降的比却逐步增大并达到10.同时,分析了地震时和岩石破裂时的应力降差别很大,以及不同方法所得破裂能差别很大的原因,认为断层的蜂窝状接触对上述问题可能有重要影响.      

等吸力三轴剪切试验中非饱和黄土强度变形及水量变化特性

通过饱和原状黄土常规三轴试验和非饱和原状黄土等吸力三轴试验研究吸力和净围压对非饱和黄土强度变形的影响,并用HUANG等、胡冉等和方祥位等提出的土水特征曲线模型分析剪切过程中排水规律。研究结果表明：等吸力下原状土样固结剪切体积变形随净围压增大而增大;等净围压下原状土样固结剪切体积变形量随吸力的增大基本呈减小趋势。p-q平面内饱和土CSL线逐渐超过低吸力下非饱和土CSL线,原因在于随着p值增大,相对于非饱和土,饱和土孔隙比越来越小,较小孔隙比对抗剪强度的贡献逐渐大于非饱和土吸力对抗剪强度的贡献。应用固结稳定的数据拟合出HUANG等和胡冉等提出的模型参数,并预测剪切过程中的排水量,发现剪切速率对排水量有影响,这两种模型适用于剪切速率慢,排水充分的情况,而方祥位等提出的模型对不同剪切速率会拟合出不同的参数。     

佳木斯地震台DSQ水管倾斜仪对日本及日本海域地震的前兆反映

大地震在主破裂前往往行成地壳变形,高应力作用下进入峰值后的变形阶段,常有大量介质破裂及裂隙间的粘滑产生,反映到地倾斜记录仪上,表现为固体潮形变异常。近年来不少研究者发现,远场前兆往往出现在对应力、应变变化反应灵敏的特殊构造部位。佳木斯地震台可能处于这种特殊的构造部位,其前兆异常对1000km以外的日本及日本海域的地震有明显短临前兆反映。     

Spontaneous Complex Earthquake Rupture Propagation

— The historical development of spontaneous rupture propagation, starting from the landmark paper of Griffith in 1920, through to the late 1980s is traced, with particular emphasis on the work carried out at MIT in the 1970s by K. Aki and his co-workers. Numerical applications of Kostrov's method for planar shear cracks were developed by Hamano, Das and Aki. Simultaneously at MIT, Madariaga considered the radiated field of a dynamic shear crack. The further development of these ideas, for example, three-dimensional spontaneous planar faulting models, continued through the 1980s. Major insight into the maximum possible rupture speeds for earthquakes developed, with the acceptance of the theoretical possibility of supersonic rupture speeds for faults with cohesion and friction, the theoretical developments spurring the search for such observations for earthquake ruptures. Possible mechanisms by which faults stop were elucidated. It was shown that a propagating rupture can jump over barriers for cracks with a cohesive zone at its tip. Complex faulting models, namely the barrier and asperity models, and their associated radiated field developed. In the late 1980s, it was shown that “dynamic” or transient asperities can develop during the complex rupturing process. Even seemingly relatively simple physical situations, can lead to such complex rupturing processes that the usual idea of “rupture velocity” needs to be abandoned in those cases. Some of the work initiated by Aki and his co-workers, such as the details of the transition from sub-Rayleigh to super-shear speeds in inplane shear mode, and the behavior of the cohesive zone size as the crack extends, still remains the subject of research today.     

Experimental investigation of three-dimensional propagation process from surface fault

Experimental study of the formation and propagation of three-dimensional (3D) faults is of great sig- nificance in the understanding of the propagation process developing from initial natural faults. In the study described in this paper, experimental investigations of 3D propagation processes of a type of surface fault are carried out under biaxial compression. The strain field near the surface fault is dy- namically measured and fully analyzed with a high-density Multi-Channel Digital Strain Gauge (MCDSG) and Digital Speckle Correlation Method (DSCM) based on the white-light image analysis. Simultane- ously the micro-fracture process involved in fault formation is observed by a 3D acoustic emission (AE) location system with a set of multi-channel whole-wave record equipment. The experimental results show that the 3D propagation process of surface fault differs greatly from that of the two-dimensional (2D) state and that a new more complicated type of 3D morphological characters and deformation mechanisms are produced. The 3D propagation process of surface faults may be divided into three stages: 1) the first stage of crack propagation initiated by wing cracks; 2) the conversion stage propa- gated by petal cracks; and 3) the second stage of crack propagation formed by shell-shaped fracture surface. The primary propagation patterns of the three stages are different. The corresponding defor- mation fields and micro-fracture distributions are likewise different. The fracture activities from petal cracks especially are of vital importance during surface fault propagation. This is also a key conversion state and marks an intrinsic difference between 2D-like and the 3D state in fault development.     

基于二维三分量伪谱法模拟数据的EDA介质中横波分裂研究

为了检测定向裂隙介质中横波分裂的方位属性特征,分析地震属性随裂隙密度和方位变化,采用人工吸收边界和反周期扩展边界,用伪谱法获得不同裂隙密度和不同方位地质模型三分量地面记录;应用时频分析和剪切波偏振分析研究由于裂隙方位和密度引起的横波分裂.结果显示,裂隙密度和方位决定着横波分裂的时差和偏振.快慢横波的延迟时间随裂隙密度增大而增加,不同方位相同裂隙密度的横波分裂时差有微小的变化.在45°方位检测时间延迟时间最大.通过时频分析,可以看到不同方位的瞬时主频有显著的变化,在横波分裂处瞬时主频有明显变化.因此,瞬时主频和快横波的偏振以及延迟时间可以作为裂隙方位和密度的指示.     

地震各向异性——多组裂隙对横波偏振的影响

通过对多分量地震资料的分析,我们发现随着频率的增加横波分裂时差减小.对于深部接收的VSP数据来说快横波的偏振方向保持不变,而对于浅层接收的VSP数据来说偏振方向却存在一个最大可以达到20°的旋转.尽管多尺度随机分布微裂隙岩石物理模型已经成功地模拟并解释了横波分裂时差随频率变化的现象,却不能解释与频率相关的横波分裂.据推测,如果微裂隙的排列方向和大裂隙的排列方向不同,利用低频信息获得的偏振方向将指示裂隙主方向,而利用高频信息获得的偏振方向则指示微裂隙方向.在背景多孔隙介质中存在多组裂隙的情况下,推导出垂直入射条件下横波偏振方向的解析式,给出了系统研究横波在介质中传播的方法.研究结果表明,横波偏振方向会随着频率的变化而变化,并且在入射方位、角度一定的条件下,是裂隙方位和密度的函数,这些认识可能有助于揭示观测到的、依赖频率变化的横波偏振现象.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

强震中短期前兆异常的共性特征

强震前兆异常具有六大共性特征。地震学方面: 第一,长期阶段的平静异常,强震均发生在已发生强震震源区之间的空段上;第二,中短期阶段强震震源区之外相当大范围中小地震或震群异常活跃,震源区反而平静,形成环形或 R-t图空区。定点前兆方面;第一,异常开始时间的外早内晚特征 ;第二,异常台站空间分布的非连续性特征; 第三,短期阶段异常数量的“突增”性特征 ;第四,短临阶段异常幅度的“巨变”性特征。上述事实用已有的DD和IPE模型无法解释,可作震源硬化模型的事实基础。在强震的预报实践中,空段和空区可以作为地点判据,异常数量的增加同时伴随的巨变异常,可以作为进入短临阶段的判据。除某些强震前1~2天外,其他时段震源区内异常不明显。按传统模型指导实践,必然造成声东击西,若用已揭示的异常共性特征指导预报,则可避免上述错误。     

具有割缝的大理岩试件破坏过程的有限元模拟

本文使用应变软化的弹塑性本构模型,用有限元方法对具有中心斜割缝的岩石试件在轴压下的破坏过程进行了数值模拟。得到了与实验观察相当一致的结果:在割线端部附近,首先发生张破裂,裂纹的扩张是逐渐地弯向试件的轴线方向,并变为稳定的。随着轴向压力增加,在张裂纹的相反方向上出现共轭的剪破裂。这些现象与割缝顶端附近材料的应变软化性质有关。说明所采用的有限元方法是研究岩石材料破裂过程的一种可行的方法。     

On faulting in a medium containing an inhomogeneity inplane and antiplane strain models

Faulting in a medium with an inhomogeneity is analysed applying two-dimensional models consisting of a shear crack in a presence of a circular inclusion. The stress drop and the stress intensity factor are calculated for mode II and III cracks of various positions in relation to the inclusion. The results demonstrate that the effect of an inhomogeneity on a shear zone strongly depends on the location of a zone for either mode II or mode III shear zone. This effect is mostly due to the spatial distribution of external effective shear stress around an inhomogeneity. Depending on the position, an inhomogeneity may have either a destabilizing effect (the stress intensity factor becomes greater) or a stabilizing influence (the stress intensity factor is decreased or faulting is prohibited so the inhomogeneity acts as an asperity or a barrier). There is a substantial difference, however, between mode II and mode III shear zones approaching an inhomogeneity centrally. Namely, the effect of inhomogeneity on the mode III shear zone located in the immediate vicinity of the inhomogeneity is in this case considerably more pronounced than that for mode II shear zone and depends to a far greater extent on the rigidity contrast between the inhomogeneity and the surrounding medium. Another important conclusion is that the quantitative effect of an inhomogeneity on faulting depends essentially on the initial value of the stress drop of a shear zone approaching an inhomogeneity, being decidedly higher for a shear zone of small stress drop. It means that in specified areas in the proximity of medium inhomogeneities one should expect substantially greater faulting activity in which weak events prevail than in other regions surrounding inhomogeneities where such activity should be distinctly reduced. Such conclusions apply to both high rigidity inhomogeneities, which, in particular, may be associated with intrusions from the upper mantle, and to low rigidity inhomogeneities such as volcanos. The present model sets forth the plausible explanation regarding why earthquakes from the same region are occasionally characterized by various values of the stress drop. The model also presents the quantitative insight concerning how heterogeneity of the medium, in the sense of spatial variation of elastic constants, affects faulting.     

Growth patterns and dynamics of mud cracks at different diagenetic stages and its geological significance

This paper discusses the growth stages,spatial structures,quantitative fitting relationships among various parameters,growth patterns and influencing factors of mud cracks by field survey,core observation and SEM analysis.The study shows that:(1) Mud crack growth can go through three stages,i.e.the syndiagenetic stage,the burial diagenetic stage(including early diagenetic stage,middle-late diagenetic stage) and the epidiagenetic stage.(2) Quantitative fitting relationships among various parameters allow a great significance to describe the spatial structure,the regional distribution and the growth environment of mud cracks.(3) Mud crack growth has three models,such as the unilateral growth model including the linear growth pattern,the curvilinear growth pattern and the bifurcation growth pattern,the multilateral growth model including the inteisectional growth pattern,the join growth pattern and the dispersed growth pattern,and the mixed growth model including the combination of any patterns listed above.(4) Modern mud crack growth usually undergoes four stages.Sand beds in sand-mud rhythmic strata can play a lubricative role on crack growth and provide enough sandy deposits for filling cracks.(5) Mud crack growth usually produces bifurcation and bifurcation angles which are mostly 120° or 90° that are related to sediment heterogeneity and released energy.(6)Factors affecting mud crack growth cover many aspects:clay content and salinity can control the number of mud cracks in different areas;terrain can control mud crack morphology;and different sedimentary cycles can control the growth patterns and filling models of mud cracks.     

SKS波对地壳裂隙各向异性的响应——理论地震图研究

穿透含裂隙、裂缝地壳8s视周期的SV波的理论地震图研究表明,当地壳平均裂隙密度高于0.01即横波各向异性高于1%时,非对称面内不同方位的SKS波均发生分裂;地震图中直接的记录显示是切向T分量上出现SKS波的振动,其振幅随地壳平均裂隙密度的增大而增强,甚至能与径向R分量上的振幅相当.局限于上地壳的强裂缝各向异性同样能引起SKS分裂.长周期SKS波分裂对地壳内裂隙、裂缝的分布缺乏分辨率.直立平行排列裂隙、裂缝使得SKS分裂T分量记录特征具有方位对称性,这来自于HTI介质中快、慢波偏振和到时差随方位变化的对称性;而倾斜裂隙、裂缝使得该方位对称性丧失.对实际观测SKS分裂的偏振解释需要考虑地壳裂隙各向异性,特别是断裂附近的强裂缝各向异性.     

地震孕育过程中应力场的一个力学模型

本文根据求得的无穷远边界上作用的外力主矢量为常量的共线双半无限长裂纹的平面问题的线性解答,分析了应力场的定性特征:张压应力区的边界和裂纹将全平面分成六个区域,其分布主要由外力主矢量的方向所决定。对线性解答及其稳定性问题的讨论,使我们可从理论上分析裂纹的扩展方式,求得较“库仑破裂准则”强的破裂条件(主要依赖于外应力主矢量的方向)並给出了破裂长度的近似表达式。根据地质力学的观点,在地震孕育的整个过程中,可假定外应力沿同一方向持续加强,计算表明弹性变形使得该方向相对于闭锁段发生偏转,从而使整个孕育过程按外应力主矢量所达到的特定方向分为两个阶段:第一阶段,裂纹无扩展,通过以时间为参数的附加应力场的计算,讨论了上述微小偏转所引起的应力场的变化;第二阶段,裂纹扩展,即闭锁段两端发生纯剪切位错,通过该双位错所引起附加应力场的计算分析了该阶段的应力特征,並探讨了地震序列与大震的关系。最后讨论了平面非共线双半无限长直裂纹问题,得到了前述类似的结论。     

非穿透裂纹平板试件三维破裂的实验研究

在研究地震孕育及震源破裂的问题时,用三维破裂来模拟地震破裂过程,比二维破裂与实际更加接近。我们设计了一组实验,采用几种不同的材料（玻璃、有机玻璃和大理石等）做成平板试件,在其中心预制非穿透裂纹,并在单轴压下进行了实验研究。 实验结果表明,非穿透裂纹与穿透裂纹的扩展情况有质的差别。非穿透裂纹的扩展并不沿原裂纹面延伸,而是沿裂纹前缘扩展出许多裂纹面,每一个破裂面都是复杂的扭曲、拐折的曲面。这可能意味着实际地震断层的扩展情况也是极为复杂的。