振动对充填裂隙渗透率影响的实验研究

为了研究振动和充填泥沙的共同作用对裂隙岩石渗透率的影响,通过对裂隙试件在充填和振动后条件下的渗流实验,得到了频率和振幅与充填裂隙岩石渗透率的关系曲线。结果表明,充填泥沙后的相对渗透率低于充填前,最小降幅为10.19%,最大降幅为59.74%。在相同的振幅下,试件的相对渗透率随频率的增加都基本呈下降趋势;振动后的渗透率低于振动前的渗透率;充填试件在较大频率振动下,渗透率减小的更为显著。在相同的频率下,充填试件振动后的渗透率基本上都减小了,最小降幅为8.90%,最大降幅为51.23%;而未充填试件振动后的渗透率与振幅存在一定关系;相对渗透率与振幅的变化之间存在拐点,相对渗透率并没有随着振幅的增加而一直减小。     

圆板状岩样破坏过程中的渗透特性实验研究

利用一种自行设计的附加渗透试验装置,在MTS815.02岩石力学伺服试验系统上完成了圆板状岩样破坏过程中的渗透特性试验,得到了圆板状岩样渗透率、渗透压差、渗流速度与轴向荷载的关系,并对试验结果进行了分析和讨论。研究表明:1)岩样在渗透的同时发生弯曲变形,岩样的渗透率与其受力状态密切相关,当达到峰值载荷时,渗透率急剧增大,由此验证了突水是伴随着结构破坏的渗流失稳过程;2)岩样在破坏过程中,当轴向载荷达到峰值时,岩样两端压力差急剧减小;3)岩样两端压力差保持0.2MPa,轴向荷载达到峰值时,增压器活塞位移急剧增大,即渗流速度突然变大;4)冲剪破坏为圆板状岩样的主要破坏形式。     

煤系地层岩石渗透特性试验研究

基于MTS815.02型岩石伺服渗透试验系统,应用两种实验方法测定了岩石的渗透特性,给出煤系地层中多种岩石全应力-应变过程中渗透率变化范围,对几种岩石的渗透特征进行了分析和讨论。研究表明:1)煤系地层岩石全应力-应变过程的渗透率通常在10-9到10-4Darcy量级之间;2)不同岩石在应力峰值前后其渗透性呈现不同的变化趋势;3)岩石渗透率的离散性很大,即使在相同的试验条件下,渗透率变异系数大多在0.65以上,其离散性按粗砂岩、泥岩、砂质页岩、中砂岩、细砂岩、砾岩、石灰岩顺序逐渐增大;4)在MTS815.02系统上,当所测岩样渗透率小于10-4Darcy量级时,应采用瞬态法,当渗透率量级大于10-4Darcy时,可改用稳态法测定。岩石变形过程中应力场和渗流场的耦合作用十分复杂,渗透率受孔隙压力、围压、试件尺寸、饱和程度等因素影响。     

裂隙岩体中非饱和渗流与运移的概念模型及数值模拟

探讨了裂隙岩体中非饱和地下水渗流与溶质运移的几种概念模型的构造及数值模拟问题 ,如裂隙网络模型、连续体模型、等效连续体模型、双孔隙度 (单渗透率 )模型、双渗透率模型、多组份连续体模型等。在裂隙岩体中 ,非饱和地下水的渗流可能只局限于岩体中的岩石组份、或裂隙网络 ,也可能在裂隙和岩石中同时发生 ;对前一种情形只需考虑单一连续体中的流动 ,而后一种情况则需要包括地下水在岩石和裂隙之间的交换。岩体中的裂隙网络往往是溶质运移的主要通道 ;但当溶质在裂隙与岩石之间的渗透和扩散是重要的运移机制时 ,就需要考虑岩石与裂隙界面处的溶质交换。为了模拟岩石与裂隙之间地下水和溶质的交换 ,就需要了解岩石与裂隙之间相互作用的模式和范围 ,使得这类问题的概念模型较单一连续体模型多了一层不确定性、其数值模拟也变得更为困难。因为在实际问题中不易、甚至根本不能判别非饱和渗流的实际形态 ,具体采用哪种模型主要取决于分析的目的和对现场数据的掌握程度。不论哪种模型都会受到模型及参数不确定性的影响 ,因此必须考虑与其他辅助模型的比较.     

变开度岩体裂隙多相渗流实验与有效渗透率模型

岩体裂隙的有效渗透率是描述岩体非饱和或多相渗流的关键参数,而裂隙开度是影响有效渗透率的重要因素.通过自主研发的粗糙裂隙多相渗流可视化实验平台,针对天然岩体裂隙复制而成的裂隙模型开展变开度条件下的多相渗流可视化实验,研究开度变化对多相渗流流动结构以及有效渗透率的影响.研究表明:非湿润相流体运动通道,在低流量比条件下呈现出气泡流流动结构,而在高流量比条件下呈现较为稳定的通道流流动结构.随着开度的增加,非湿润相流动通道的分支变少、等效宽度增加,两相流体的有效渗透率均增大,流动结构趋于稳定.可视化结果还阐明了柱塞流流动结构下,两相流体交替占据裂隙空间的竞争机制:当非湿润相流体通道由连续转变为不连续时,裂隙进出口压差显著增加;反之,当该通道由不连续转变为连续时,压差显著减小.最后,基于分形理论以及渗透率统计建模方法,建立了考虑开度效应的岩体裂隙多相渗流有效渗透率理论模型,并通过实验测定的有效渗透率数据验证了该模型的正确性与有效性.     

燃爆冲击作用下岩石初始破坏区形成机制与主控因素

为揭示燃爆冲击作用下井周岩石破坏区的形成机制,并分析影响初始破坏区(破碎区和初始裂隙区)的主控因素,开展了两种岩样在不同加载速率下的冲击破坏实验, 分析了岩石冲击破坏模式及岩石对加载速率的响应, 借助基于Von Mise准则建立的岩石冲击破坏的破碎区和初始裂隙区计算模型可知:加载速率低于190 GPa/s时,可依据冲击峰值压力引导的应力分布确定破碎区和初始裂隙区作用范围;燃爆压裂在近井地带主要产生破碎区和裂隙区,破碎区直径为井眼直径的1~3倍,初始裂隙区直径为井眼直径的5~7倍;冲击载荷作用下,初始破坏区与加载速率、脆性指数呈正相关,且受脆性指数影响更显著。研究结果可提高对燃爆压裂过程中岩石的破坏模式及其主控因素的认识深度,为燃爆压裂冲击条件设计提供指导。     

冲击载荷下岩石压胀对其性质的影响

采用落锤动态加载岩石实验系统,通过调节落锤质量、下落距离及活塞杆垫板材料,记录岩石瞬间形成冲击压力脉冲,研究压胀对不同岩石性质影响的规律.实验结果表明,压胀产生后岩石的性质发生了改变,岩石的渗透率有不同程度的增加,且岩石越致密,渗透率增加倍数越大;岩石的弹性模量和弹性极限随岩石体积的增加而降低;压胀产生后纵波和横波在岩...     

卸压开采损伤煤岩气体渗透率试验研究

我国煤系地层的煤岩体渗透率普遍较低,直接影响煤矿瓦斯抽采效果。卸压开采使采场围岩受采动应力作用具备了峰后力学行为,其渗透率也大幅度增加,有助于瓦斯的高效抽采。本文基于MTS815.02型电液伺服岩石渗透试验系统,对卸压开采后损伤煤岩的气体渗透率进行了测试。研究表明:(1)卸压损伤煤岩石的渗透率在10-5～10-4 m2之间,较受采动影响前要大107～108倍;(2)卸压过程中,煤样的气体渗透率与围压近似呈线性关系,与岩样呈双曲线关系;(3)卸压损伤煤岩逐渐加围压,气体渗透率近似为线性下降,卸压阶段,其渗透率增加;(4)加卸围压对煤岩的气体渗透率有很大影响。研究结果可为卸压瓦斯高效抽采提供理论依据。     

明渠中悬移质的弥散-对流方程及悬浮机理

悬移质泥沙通常构成冲积河流总输沙量的主体, 研究悬移质的悬浮机理具有重要的意义. 以双流体模型为基础, 通过引入弥散速度的概念, 建立了悬移质泥沙的输沙方程以及泥沙扩散系数的本构关系. 应用该方程分析了二维明渠均匀流中悬移质泥沙浓度垂向分布规律, 并与Einstein 和Chien 的泥沙浓度实验资料及经典扩散理论进行了对比. 以此为基础, 分析了紊动扩散、颗粒自身的紊动、颗粒碰撞应力对泥沙悬浮的影响在垂向上的变化, 以及浓度、粒径等对这些因素的影响. 结果表明, 泥沙颗粒在明渠紊流中的扩散是浑水的紊动扩散、颗粒自身的紊动、颗粒碰撞应力3 部分不同机制共同作用的结果, 把泥沙颗粒的悬浮简单归因于水流的紊动是不全面的.     

大理岩试样循环加载强化作用的试验研究

利用不同晶粒大理岩块加工成完整和含孔道的圆柱试样,在伺服实验机上分别进行单轴和常规三轴循环加载试验,其中三轴循环加载包括维持围压恒定和卸载围压两种情况.结果表明,循环加载确实可使试样的强度增加,强度大致提高5%-10%.岩石的材料强度和承载能力是两个不同的概念.大理岩晶粒之间多为裂隙,轴向压缩加载时晶粒局部接触应力远高于名义应力,在岩样未达到峰值应力之前,接触处已经产生很大的变形乃至出现局部破坏,形成的碎屑在卸载时可以脱落充填到附近的空隙,提高岩样承载能力.这种强化特征是与摩擦相关的承载能力,并非材料强度;与多次循环加载造成材料强度逐步劣化不同.从卸载点附近卸载与再加载应力应变曲线的关系,可以确认卸载是否引起岩石承载能力的增加.     

Permeability of a random array of fractures of widely varying apertures

We modelize a fractured rock by a random array of plane cracks of finite extent having a very broad distribution of apertures (or of hydraulic conductances). If the rock is permeable, the flow will essentially take place along a ‘subnetwork’ made of the less resistant cracks. Using an analogy with the treatment of variable range transport in semiconductors, we evaluate the homogenization length and the permeability of this disordered network. This evaluation makes use of the notion of the critical bonds which are the weakest cracks among the good ones necessary for percolation; the remaining weaker bonds make a negligible contribution to the permeability. The method is applicable to other examples of transport in very heterogeneous macroscopic random materials.     

Experimental Investigation on Transport Properties of Cement-Based Materials Incorporating 2D Crack Networks

This paper investigates the correlation between the geometry of crack networks and the altered transport properties of cement-based porous materials. Cracks were artificially introduced into slice specimens to obtain bidimensional (2D) crack networks, and the network was characterized by the crack density, orientation, connectivity and crack opening aperture. For the permeability, the water vapor sorption isotherms were measured and an algorithm was established to solve the intrinsic permeability of cracked specimens with the help of moisture transport modeling and the data of drying tests. The electrical conductivity of cracked specimens was measured using an alternative current method. The study on the specimens with percolated cracks shows that: (1) the pertinent geometry parameters for altered transport properties include average-based crack density, crack opening and local crack connectivity; (2) the water permeability of cracked specimens is correlated to the combination \(b^{1.7}\rho f\) and electrical conductivity to \(b^{0.45}\rho f\); (3) the different exponents on the crack opening/length ratio reflect the resistance of tortuosity of crack paths to the water and current flow and this resistance is stronger for current flow.     

P- and S-wave velocity and anisotropy in saturated rocks with aligned cracks

Rocks in the Earth’s crust contain cracks in variable amounts and on variable scales. Depending on the in situ stress level, cracks are always distributed in a common normal direction (so-called aligned cracks). In this study, we focus on the velocity and anisotropy of perfectly aligned cracks (parallel cracks). The theoretical model is introduced to calculate the elastic properties of cracked rocks based on solid matrix properties, crack parameters, and fluid properties. The theoretical results using high- and low-frequency limits and intermediate frequency for velocities and anisotropy are compared with experimental data measured in the laboratory, where rocks are saturated with water and air. We used four synthetic rocks with different crack densities (0.0, 0.0243, 0.0486, and 0.0729) in this study. The comparison shows that the results for high, low, and intermediate frequencies provide good predictions of the shear velocities and anisotropy for saturated cracked rocks. Quantitatively, the results of the P-wave velocity and anisotropy for intermediate frequency fit the measured data better than that for the high- and low-frequency limits. Compared with water saturation, dry cracks with high compliance have a strong influence on the total compliance of rock. The water saturation results in high P-velocities and a low P-anisotropy, while dry cracks induce much lower P-velocities and a high P-anisotropy. The shear wave anisotropy is fluid-independent, theoretical results fit the measured data well in both water saturation and the dry case.     

开裂孔隙材料渗透率的细观力学模型研究

采用细观力学方法对含随机裂纹网络的孔隙材料渗透性进行研究.开裂孔隙材料渗透性的影响因素包括裂纹网络的密度、连通度、裂纹的开度以及孔隙材料基体渗透性.对于不连通的裂纹网络,该文采用已有的相互作用直推法(interaction direct derivative,IDD)的理论框架,引入裂纹的密度$\rho$和裂纹开度比$b$,提出了裂纹夹杂$\!$-$\!$-$\!$基体两相复合材料渗透率的IDD理论解.对于部分连通裂纹网络,考虑局部裂纹团内部各个裂纹对有效渗透率的相互放大作用,引入裂纹网络的连通度$f$,定义与连通度相关的水平裂纹密度$\rho^{h}$,按照增量法将表征连通特征的水平裂纹嵌入有效基体中,以此方式来考虑裂纹夹杂间的相互搭接,提出了考虑裂纹连通特征的扩展IDD理论解,分别考虑了基体材料渗透率$K_{m}$、裂纹密度$\rho $、裂纹开度比$b$以及与连通度$f$相关的$\rho ^{\rm h}$.最后通过对有限区域内含随机裂纹网络孔隙材料渗透过程的有限元模拟分别验证了不连通和部分连通裂纹网络扩展IDD模型的适用性:(1)当裂纹不连通时,由于基体对流体渗透的阻隔作用,裂纹的开度对有效渗透率影响不大;(2)当裂纹部分连通时,裂纹密度分别小于1.1(无关联裂纹网络,分形维数为2.0)、1.2(关联裂纹网络,分形维数为1.75)时,扩展IDD模型能够很好地估计开裂孔隙材料的有效渗透率,但是随着裂纹进一步扩展,最大裂纹团主导作用凸显,扩展IDD模型不再适用.      

Numerical Modelling of Steady-State Flow in 2D Cracked Anisotropic Porous Media by Singular Integral Equations Method

The equations governing plane steady-state flow in heterogeneous porous media containing curved-line intersecting cracks (Pouya and Ghabezloo in Transp Porous Media 84:511?C532, 2010) and the potential solution obtained for these equations are considered here. The theoretical results are first completed for the mass balance at crack intersections points. Then, a numerical procedure based on a singular integral equations method is described concretely to derive this solution for cracked materials. Closed-form expressions of elementary integrals for special choice of collocation points lead to a very quick and easy numerical method. It is shown that this method can be applied efficiently to the study of the steady-state flow in cracked materials with anisotropic matrix permeability and a dense distribution of curved-line intersecting cracks. Some applications of this method to the permeability of cracked materials are given.     

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.     

岩体裂隙三轴应力渗透规律的试验研究

利用三轴应力渗透仪分别对煤矿底板的硬岩和软岩进行三轴应力渗透试验 ,观测渗水量、围压、轴压及水压的变化。     

裂隙型硬粘土的力学模型及其在边坡工程中的应用

硬粘土是介于软岩和硬土之间的一种介质, 其中含有规律分布的裂隙, 具有明显的方向性和各向异性特征; 本文从实际应用的角度出发, 提出一种适用于裂隙性硬粘土的损伤力学模型, 并应用于南水北调中线渠道边坡工程的数值模拟分析。     

Effective permeability of cracked unsaturated porous materials

This study focuses on a theoretical estimation of the effective permeability of unsaturated cracked porous media. The closed-form flow solution around and in a superconductive crack, embedded in an infinite porous matrix under a far-field condition, is recalled first. Then the solution of flow around a completely unsaturated (empty) crack that is considered as an obstruction against the flow is determined. The flow solution for partially saturated crack in special configurations is obtained by superposition of the two basic solutions for superconductive and empty cracks. The contribution of an unsaturated crack, with a given saturation degree, to the effective permeability is estimated by using dilute upscaling scheme. Numerical results obtained by Finite Elements Method, are in good agreement with the theoretical results for weak crack densities but show the additional effect of cracks interaction for higher densities.     

Nonlinear and parametric coupled vibrations of the rotor-shaft system as fault identification symptom using stochastic methods

In the paper several stochastic methods for detection and identification of cracks in the shafts of rotating machines are proposed. All these methods are based on the Monte Carlo simulations of the rotor-shaft lateral-torsional-longitudinal vibrations mutually coupled by transverse cracks of randomly selected depths and locations on the shaft. For this purpose there is applied a structural hybrid model of a real cracked rotor-shaft. This model is characterized by a high practical reliability and great computational efficiency, so important for hundreds of thousands numerical simulations necessary to build databases used in solving the inverse problem, i.e. crack parameter identifications. In order to ensure a good identification accuracy, for creating the Monte Carlo samples of data points there are proposed special probability density functions for locations and depths of the crack. Such an approach helps in enhancing databases corresponding to the most probable faults of the rotor-shaft system of the considered rotor machine. In the presented study six different database sizes are considered to compare identification efficiency and accuracy of considered methods. A sufficiently large database enables us to estimate almost immediately (usually in less than one second) the crack parameters with precision that is in most of the cases acceptable in practice. Then, as a next stage, one of the proposed fast improvement algorithms can be applied to refine identification results in a reasonable time. The proposed methods seem to provide very convenient diagnostic tools for industrial applications.