非饱和土中热弹性波的传播特性分析

基于非饱和多孔介质的研究成果,考虑热效应和孔隙流体迂曲度的影响,研究了非饱和土中热弹性波的传播特性。利用非饱和土中耦合热的固-液-气三相介质的质量平衡方程、渗流连续方程、动量平衡方程和广义非Fourier热传导定律,建立了问题的热弹性波动方程。通过引入势函数,经过理论推导给出了非饱和土中热弹性波的弥散特征方程。结合数值算例,分析了几类热弹性波的波速和衰减系数随迂曲度、热膨胀系数和介质温度等热物理参数的变化规律。结果表明:孔隙水迂曲度的增大将引起P₁波、P₃波和S波的波速增大,而孔隙气体迂曲度的增大仅使得P₂波的波速增大;热膨胀系数的增大将造成P₁波波速的增大和热(T)波波速的减小;介质温度的升高将引起各类热弹性波波速的增大;频率、热膨胀系数和介质温度的变化对各类热弹性波的衰减系数均有较大影响,不可忽视。     

三相非饱和土参数对波的传播的影响研究

基于已建立的非饱和土中的波动方程,导出4种体波的波速以及衰减的解析表达式。数值分析了4种体波的波速、衰减与孔隙率以及Lame常数 、 之间的关系。结果表明,随孔隙率的增大,P1波和S波的波速有所提高,P2波和P3波的波速随孔隙率的增大而有所减小;3种P波的衰减有所增大,其中以P3波变化最为显著。P1波、P3波和S波的波速基本上不随Lame常数 变化,P2波的波速随Lame常数 的增大而稍有增大。P2及S波的衰减基本不随Lame常数 变化,P1、P3波的衰减随Lame常数 的增大而有所减小。3种P波的波速基本上不随Lame常数 变化,S波的波速随Lame常数 的增大而有所增大。P2及S波的衰减基本不随Lame常数 变化,P1、P3波的衰减随Lame常数 的增大而有所减小,其中P3波的变化最为显著。     

应变软化土坡渐进破坏的可靠度分析

在考虑土性参数的不确定性和土体应变软化特性的基础上,以可靠度方法建立了边坡渐进破坏的极限状态方程。通过计算边坡局部破坏概率,确定边坡中最危险的土条,并推导出破坏起始于该土条,逐渐向两侧扩展的渐进破坏概率计算公式。该方法能很好地反映边坡局部破坏的产生、扩展及对边坡整体稳定性的影响,找到最贴近实际情况的破坏传播路径,为边坡稳定分析提供了理论依据。通过算例分析可知,滑动面上各个土条的局部破坏概率均不相等,破坏起始于中部第8土条,破坏先向上传递,后向下传递至坡脚,并非从坡脚向坡顶传递的牵引式破坏,或从坡顶向坡脚传递的推移式破坏。     

Hydro‐mechanical modeling of cohesive crack propagation in multiphase porous media using the extended finite element method

In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.     

基于水文模型的雷达监测降雨量误差传递研究

基于多普勒雷达和相应的雨量计资料,利用卡尔曼滤波校准法对分组Z~I关系估算的雷达降雨进行同化,结合新安江模型,提出采用增长模繁殖法对雷达监测降雨量资料进行扰动,定量分析了模型输入误差对径流模拟的影响。经过湖北省白莲河流域的实际应用研究表明,卡尔曼滤波校准法估算降水量的精度比Z~I关系法有了明显的提高,其相应的洪水预报效果也都优于Z~I关系;提出的增长模繁殖法对雷达监测降雨量资料产生的扰动输入误差在经过模型的传递后有增大的趋势。     

Effects of medium and pore‐fluid compressibility on chemical‐dissolution front instability in fluid‐saturated porous media

Theoretical analysis and computational simulations have been carried out to investigate how medium and pore‐fluid compressibility affects the chemical‐dissolution front propagation, which is associated with a fully‐coupled nonlinear problem between porosity, pore‐fluid pressure, pore‐fluid density and reactive chemical‐species transport within a deformable and fluid‐saturated porous medium. When the fully‐coupled nonlinear system is in a subcritical state, some analytical solutions have been derived for a special case, in which the ratio of the equilibrium concentration to the solid molar density of the chemical species is approaching zero. To investigate the effect of either medium compressibility or pore‐fluid compressibility on the evolutions of chemical dissolution fronts in supercritical chemical dissolution systems, numerical algorithms and procedures have been also proposed. The related theoretical and numerical results have demonstrated that: (i) not only can pore‐fluid compressibility affect the propagating speeds of chemical dissolution fronts in both subcritical and supercritical systems, but also it can affect the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; (ii) medium compressibility may have a little influence on the propagating speeds of chemical dissolution fronts, but it can have significant effects on the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; and (iii) both medium and pore‐fluid compressibility may stabilize irregular chemical‐dissolution‐fronts in supercritical chemical dissolution systems. Copyright © 2011 John Wiley & Sons, Ltd.     

冲击载荷作用下岩石动态断裂试验研究

以自行研制的可调速落锤冲击试验机进行加载,通过高速相机搭建数据采集系统,采用数字散斑相关方法作为试验的观测方法,开展5组不同预制裂纹长度的花岗岩试件在冲击载荷作用下的动态断裂试验。并对岩石I型裂纹在冲击载荷作用下的位移场演化、裂纹动态断裂的裂尖张开位移、裂纹尖端的扩展历史、动态断裂的应力强度因子进行了定量研究。试验结果表明：对于单一试件而言,在落锤开始与岩石试件的边界接触到试件的预制裂纹开始扩展整个过程中,裂纹尖端的应力强度因子KI呈递增趋势,KII与KI相比,相差1～2个数量级。对于不同预制裂纹试件,在落锤冲击加载作用下,总体表现为一定的预制裂纹长度范围内,随着预制裂纹的长度增加,裂纹尖端应力强度因子呈增加的趋势。     

祁连山北缘玉门砾岩的磁性地层年代与褶皱过程*

在我国西部致使玉门砾岩强烈褶皱变形的构造运动是喜马拉雅运动中重要的一幕。祁连山北缘老君庙逆断裂-褶皱带前翼发育的生长地层为认识玉门砾岩褶皱的过程提供了重要证据。通过对这套生长地层及其内部渐进不整合几何形态和结构的详细填图以及磁性地层年代学研究,认为玉门砾岩的底界具有穿时性,其年龄可能由东南向西北变小,其底界年龄在牛胳套剖面约4Ma,在青草湾西剖面约3.55Ma。玉门砾岩的褶皱变形在牛胳套剖面起始于约3Ma,沉积速率由约260m/Ma变为170m/Ma,在青草湾西剖面起始于约1.2Ma。结合构造变形的样式、强度及侵蚀程度推断老君庙逆断裂-褶皱带是以约7.8km/Ma的速率由东南向西北侧向扩展生长的。发育在背斜前翼的多个退覆-超覆生长楔表明玉门砾岩的褶皱变形由3~4个活跃期和平静期组成,这反映了祁连山及其北部前陆盆地内褶皱逆冲席体向前陆方向的不断扩展生长过程, 但并不能简单地以此作为整个祁连山或青藏高原多次垂直隆升的证据。     

Shear wave anisotropy in the crust around the San Andreas fault near Parkfield: spatial and temporal analysis

We systematically analysed shear wave splitting (SWS) for seismic data observed at a temporary array and two permanent networks around the San Andreas Fault (SAF) Observatory at Depth. The purpose was to investigate the spatial distribution of crustal shear wave anisotropy around the SAF in this segment and its temporal behaviour in relation to the occurrence of the 2004 Parkfield M 6.0 earthquake. The dense coverage of the networks, the accurate locations of earthquakes and the high-resolution velocity model provide a unique opportunity to investigate anisotropy in detail around the SAF zone. The results show that the primary fast polarization directions (PDs) in the region including the SAF zone and the northeast side of the fault are NW–SE, nearly parallel or subparallel to the SAF strike. Some measurements on the southwest side of the fault are oriented to the NNE–SSW direction, approximately parallel to the direction of local maximum horizontal compressive stress. There are also a few areas in which the observed fast PDs do not fit into this general pattern. The strong spatial variations in both the measured fast PDs and time delays reveal the extreme complexity of shear wave anisotropy in the area. The top 2–3 km of the crust appears to contribute the most to the observed time delays; however substantial anisotropy could extend to as deep as 7–8 km in the region. The average time delay in the region is about 0.06 s. We also analysed temporal patterns of SWS parameters in a nearly 4-yr period around the 2004 Parkfield main shock based on similar events. The results show that there are no appreciable precursory, coseismic, or post-seismic temporal changes of SWS in a region near the rupture of an M 6.0 earthquake, about 15 km away from its epicentre.