断层内破碎岩石渗透性的实验研究

断层破碎带岩石的渗透性是影响地下工程中由渗流失稳而引发灾害事故的重要影响因素之一。利用一种专门的破碎岩石压实渗透仪,在MTS815.02岩石力学伺服试验系统上测试了取自断层内的破碎灰岩、破碎砂岩、破碎泥岩的渗透系数,得到并分析了轴向压力、岩样粒径、水流速度与渗透系数的关系曲线。研究表明:破碎岩石渗透系数与其压实状态密切相关,随着轴压的增加,渗透系数都相应降低;在相同的轴向压力状态下,岩样粒径越小,其渗透系数也越小,因此,挤压程度和充填程度都较低的小断层导水性好,容易导致突水;破碎岩石的渗透系数随着水流速度的增大呈减小趋势,但水流速度对破碎岩石渗透系数的影响程度有限;在相同的压力、粒径情况下,破碎泥岩的渗透系数比破碎灰岩、破碎砂岩的要小1～2个量级,说明当断层上下盘为坚硬岩石时,断层带渗透性好、易导水。研究结果可为采动过程中评价断层的导水性提供参考。     

破碎岩石气体渗透性的试验测定方法

为了测试采矿工程中松散破碎岩体的气体渗透特性,设计了与MTS815.02岩石力学伺服机配套的气体渗透仪和测试系统。依靠MTS伺服机改变压力和位移条件,渗透仪用以稳定破碎岩石并连接气体测试系统,接入气源、组成回路和测试气透特性等由测试系统完成。通过对几组破碎砂岩岩样氮气渗透性能的测试,结合考虑状态方程的渗透率换算公式,得出破碎砂岩气体渗透率随压力和粒径条件的变化规律。该测试方法和测试结果在研究采动破碎岩体瓦斯等气体流动规律方面具有重要意义。     

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

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

Experimental Study of Seepage Properties of Broken Sandstone Under Different Porosities

In coal mining the water flow in broken rock is a very common phenomenon. Study of seepage properties of broken rock is one of the basic subjects required in order to understand the stability of rock surrounding roadways, preventing disasters such as water inrush and gas outbursts and developing underground resources. So far, quantitative studies on the nonlinear seepage properties of broken sandstone under different porosities are not extensive in the research literature. In this article, by means of an electro-hydraulic servo-controlled test system (MTS815.02) and a patent seepage device, the seepage properties under different conditions of porosity were tested on broken sandstone of five different grain sizes. Based on the loading method of controlling the axial compression displacement and steady permeating method, we obtained curves of the relation of pore pressure with time, as well as the relation curves between the pore pressure gradient for steady seepage and velocity. Furthermore, we calculated the permeability k and non-Darcy coefficient β corresponding to different porosities by fitting these curves with the binomial expression. This study indicates that: (1) the seepage properties of broken sandstone are closely related to grain size, load levels, and porosity structure; (2) the permeability k decreases, while the coefficient β increases with a decrease in porosity φ, but both the k − φ and the β − φ curves show some local fluctuations; (3) the permeability k of the broken sandstone has a magnitude of 10⁻¹⁴–10⁻¹² m², while the coefficient β ranges from 10¹⁰ to 10¹² m⁻¹. The results obtained provide some information for further study of the nonlinear seepage behavior of broken rock theoretically.     

软硬相间层状复杂岩体综合变形模量原位试验研究

某拟建大型电站坝基岩体为一套河湖相沉积岩,包含了砾岩、砂岩、粉砂岩和泥岩等多种岩性岩体。不同岩体变形参数各异,常规变形试验只能获得单一岩性岩体的变形模量,而这种软硬相间层状复杂岩体的综合变形模量是与单一岩性岩体变形模量、岩层倾角、岩层厚度多个因素相关的,如何通过试验获得这种复杂组合岩体的综合变形模量是一个重要问题。本文首先对坝址区不同岩性所占的比例进行了详细的统计,获得砂岩、粉砂岩和泥岩的厚度比为6：2.5：1.5,根据这一统计结果在坝轴线附近选择含有砂岩、粉砂岩和泥岩三种岩体的大型变形试验区,将承压板尺寸设定为200cm×37cm,承压板就能同时覆盖砂岩、粉砂岩和泥岩三种岩体,且将其厚度比确定为6：2.5：1.5。通过在承压板上施加大荷载（250余吨）将应力传至下部多层岩体,从而获得软硬相间层状复杂岩体的综合变形模量。试验结果表明坝基弱下风化岩体的综合变形模量可达到7.0GPa,已经满足筑坝的要求。通过原位试验成功地论证了坝基岩体综合变形模量量值的问题。
     

An Analytical Model of Porosity–Permeability for Porous and Fractured Media

The classic Kozeny–Carman equation (KC) uses parameters that are empirically based or not readily measureable for predicting the permeability of unfractured porous media. Numerous published KC modifications share this disadvantage, which potentially limits the range of conditions under which the equations are applicable. It is not straightforward to formulate non-empirical general approaches due to the challenges of representing complex pore and fracture networks. Fractal-based expressions are increasingly popular in this regard, but have not yet been applied accurately and without empirical constants to estimating rock permeability. This study introduces a general non-empirical analytical KC-type expression for predicting matrix and fracture permeability during single-phase flow. It uses fractal methods to characterize geometric factors such as pore connectivity, non-uniform grain or crystal size distribution, pore arrangement, and fracture distribution in relation to pore distribution. Advances include (i) modification of the fractal approach used by Yu and coworkers for industrial applications to formulate KC-type expressions that are consistent with pore size observations on rocks. (ii) Consideration of cross-flow between pores that adhere to a fractal size distribution. (iii) Extension of the classic KC equation to fractured media absent empirical constants, a particular contribution of the study. Predictions based on the novel expression correspond well to measured matrix and fracture permeability data from natural sandstone and carbonate rocks, although the currently available dataset for fractures is sparse. The correspondence between model calculation results and matrix data is better than for existing models.     

三种红层岩石常规三轴压缩下的强度与变形特性研究

利用MTS815Teststar程控伺服岩石力学试验系统研究了川东地区一红层边坡中的砂岩、粉砂岩和泥岩围压为0³MPa的应力-应变全过程曲线,建立了峰值强度、峰值强度前弹性模量以及峰值强度后的弹性模量和围压的关系。将低围压下红层的全应力-应变曲线概化成5个阶段,分别为压密段、弹性段、屈服段、应变软化段和塑性流动阶段。试验结果得出,红层弹性模量随围压的增加而提高且变化明显,砂岩和粉砂岩在此围压内为脆性破坏,泥岩为塑性破坏的规律。     

贝雷砂岩射孔孔道流动效应正交计算分析

对射孔孔道流动效应进行了量纲分析,建立了轴向和径向流量与影响流量因素间的工程表达式。综合考虑影响流动效应的主要因素和取值范围,设计了数值模拟计算工况的正交计算方案,并利用ANSYS/FLOTRAN流体程序进行了带损伤区贝雷砂岩流量数值模拟研究,获得了各个计算方案的径向流量和轴向流量,用以确定出流量公式的系数。最终确定的流量公式表明,各要素对流量影响的重要顺序依次为固有渗透率、损害区厚度、孔深、渗透率剩余百分数和孔径。     

A New Method for Calculating Two-Phase Relative Permeability from Resistivity Data in Porous Media

Many resistivity data from laboratory measurements and well logging are available. Papers on the relationship between resistivity and relative permeability have been few. To this end, a new method was developed to infer two-phase relative permeability from the resistivity data in a consolidated porous medium. It was found that the wetting phase relative permeability is inversely proportional to the resistivity index of a porous medium. The proposed model was verified using the experimental data in different rocks (Berea, Boise sandstone, and limestone) at different temperatures up to 300°F. The results demonstrated that the oil and water relative permeabilities calculated from the experimental resistivity data by using the model proposed in this article were close to those calculated from the capillary pressure data in the rock samples with different porosities and permeabilities. The results demonstrated that the proposed approach to calculating two-phase relative permeability from resistivity data works satisfactorily in the cases studied.     

Attrition and particle breakage under monotonic and cyclic loadingAttrition et rupture de particules sous chargements monotones et cycliques

The morphological evolution of gravel grains under monotonic and cyclic loading, typically present in road structures, is presented. Two types of mineralogy are considered: a limestone and a sandstone. Loading in uniaxial strain conditions are applied, with a maximum mean pressure of 5 MPa and a maximum number of cycles of 250?000. Variations in grain size distribution and grain shape, measured with image analysis techniques, are discussed in detail. In the case of uniform grain size distribution, particle breakage does appear for stresses around 500 kPa for the limestone grains and 100 kPa for the sandstone grains. For well-graded materials under low cyclic pressures (<1 MPa), changes in grain size distribution can hardly be detected. To cite this article: F. Mayoraz et al., C. R. Mecanique 334 (2006).     

Stress-Dependent Porosity and Permeability of Porous Rocks Represented by a Mechanistic Elastic Cylindrical Pore-Shell Model

The stress dependency of the porosity and permeability of porous rocks is described theoretically by representing the preferential flow paths in heterogeneous porous rocks by a bundle of tortuous cylindrical elastic tubes. A Lamé-type equation is applied to relate the radial displacement of the internal wall of the cylindrical elastic tubes and the porosity to the variation of the pore fluid pressure. The variation of the permeability of porous rocks by effective stress is determined by incorporating the radial displacement of the internal wall of the cylindrical elastic tubes into the Kozeny–Carman relationship. The fully analytical solutions of the mechanistic elastic pore-shell model developed by combining the Lamé and Kozeny–Carman equations are shown to lead to very accurate correlations of the stress dependency of both the porosity and the permeability of porous rocks.

     

A Fractal Model for Gas–Water Relative Permeability in Inorganic Shale with Nanoscale Pores

A reliable gas–water relative permeability model in shale is extremely important for the accurate numerical simulation of gas–water two-phase flow (e.g., fracturing fluid flowback) in gas-shale reservoirs, which has important implication for the economic development of gas-shale reservoir. A gas–water relative permeability model in inorganic shale with nanoscale pores at laboratory condition and reservoir condition was proposed based on the fractal scaling theory and modified non-slip boundary of continuity equation in the nanotube. The model not only considers the gas slippage in the entire Knudsen regime, multilayer sticking (near-wall high-viscosity water) and the quantified thickness of water film, but also combines the real gas effect and stress dependence effect. The presented model has been validated by various experiments data of sandstone with microscale pores and bulk shale with nanoscale pores. The results show that: (1) The Knudsen diffusion and slippage effects enhance the gas relative permeability dramatically; however, it is not obviously affected at high pressure. (2) The multilayer sticking effect and water film should not be neglected: the multilayer sticking would reduce the water relative permeability as well as slightly decrease gas relative permeability, and the film flow has a negative impact on both of the gas and water relative permeability. (3) The increased fractal dimension for pore size distribution or tortuosity would increase gas relative permeability but decrease the water relative permeability for a given saturation; however, the effect on relative permeability is not that notable. (4) The real gas effect is beneficial for the gas relative permeability, and the influence is considerable when the pressure is high enough and when the nanopores of bulk shale are mostly with smaller size. For the stress dependence, not like the intrinsic permeability, none of the gas or water relative permeability is sensitive to the net pressure and it can be ignored completely.     

不同尺寸砂岩动态力学性质和应力平衡性的试验研究

采用大直径分离式霍普金森压杆系统获得的不同尺寸试样的实验冲击动态力学参数有差异,因此在直径100 mm压杆上进行了3种直径（50、75和100 mm）和5种长径比（0.4、0.5、0.6、0.8和1.0）的砂岩试样冲击试验,分析了不同尺寸试样应力-应变曲线和应变率曲线随尺寸的变化,提出了用于比较波形对齐重合度的波形叠加系数,并与应力平衡因子共同构建了动态应力平衡性研究体系,由此确定大直径霍普金森压杆试验的试样建议尺寸。同时,利用高速摄影机监测试样的动态破坏状况。结果表明:当长径比相同时,直径75与100 mm岩石试样的动态抗压强度测试值相近,直径50 mm试样具有更明显的长度效应;随着试样直径的增大,应变率曲线从单峰变为双峰;小尺寸试样更易发生轴向劈裂破坏,大尺寸试样受内部应力波叠加影响产生了较大的拉应力,易发生层裂拉伸和轴向劈裂的复合型破坏;对直径75 mm且长径比0.3~0.4的试样,波形对齐后重合度较高,在起始破坏前拥有充足的应力平衡时间,应变率加载效果较好。获得了砂岩试样冲击压缩试验的尺寸效应,可为大直径岩石试样的尺寸选择提供参考。     

Efficiency of diffusion controlled miscible displacement in fractured porous media

Experiments were performed to study the diffusion process between matrix and fracture while there is flow in fracture. 2-inch diameter and 6-inch length Berea sandstone and Indiana limestone samples were cut cylindrically. An artificial fracture spanning between injection and production ends was created and the sample was coated with heat-shrinkable teflon tube. A miscible solvent (heptane) was injected from one end of the core saturated with oil at a constant rate. The effects of (a) oil type (mineral oil and kerosene), (b) injection rates, (c) orientation of the core, (d) matrix wettability, (e) core type (a sandstone and a limestone), and (f) amount of water in matrix on the oil recovery performance were examined. The process efficiency in terms of the time required for the recovery as well as the amount of solvent injected was also investigated. It is expected that the experimental results will be useful in deriving the matrix–fracture transfer function by diffusion that is controlled by the flow rate, matrix and fluid properties.     

基于高速3D-DIC技术的砂岩动力特性粒径效应研究

利用分离式霍普金森压杆（SHPB）,对粗砂岩、中等粒径砂岩和细砂岩进行了应变率为69～83 s–1的动态单轴抗压实验,研究了粒径尺寸效应对砂岩动力特性的影响。通过三维数字图像相关（3D-DIC）技术分析高速摄像图像,获得了砂岩的实时应变场,据此分析了动态荷载下3种粒径砂岩的动力变形特性和裂纹开展行为。结果表明,砂岩弹性应变储能可逆释放的临界应变率随着粒径的减小而升高,动态压缩强度随着粒径减小而增大,动态强度应变率敏感度则与强度规律相反。相较于静态条件下,中等粒径砂岩和细砂岩的动态弹性模量增长了2～3倍,粗砂岩的动态弹性模量增长达5倍以上。细砂岩的动态泊松比相较于静态提高了约25%,中等粒径砂岩的动态泊松比约为静态时的70%。动态裂纹首先出现于试件内部,然后传播至表面,呈现出应变局部化,动态荷载下岩石裂纹的孕育和扩展相比静态条件下均有所提前,其中细砂岩在动态荷载条件下的归一化裂纹起裂阈值仅为峰值强度的10%。微观分析表明,矿物粒径大小和黏土矿物含量分别在砂岩的动力力学性质和裂纹开展行为方面发挥主要作用。     

A Novel Method for Gas–Water Relative Permeability Measurement of Coal Using NMR Relaxation

Using the conventional volumetric method in unsteady-state relative permeability measurements for unconventional gas reservoirs, such as coal and gas shale, is a significant challenge because the movable water volume in coal or shale is too small to be detected. Moreover, the dead volume in the measurement system adds extra inaccuracy to the displaced water determination. In this study, a low-field nuclear magnetic resonance (NMR) spectrometer was introduced into a custom-built relative permeability measurement apparatus, and a new method was developed to accurately quantify the displaced water, avoiding the drawback of the dead volume. The changes of water in the coal matrix and cleats were monitored during the unsteady-state displacement experiments. Relative permeability curves for two Chinese anthracite and bituminous coals were obtained, matching the existing research results from the Chinese coalbed methane area. Moreover, the influences of confining pressure on the shape of the relative permeability curve were evaluated. Although uncertainties and limits exist, the NMR-based method is a practical and applicable method to evaluate the gas/water relative permeability of ultra-low permeability rocks.     

Numerical Analysis of Imbibition Front Evolution in Fractured Sandstone under Capillary-Dominated Conditions

Fractures serve as primary conduits having a great impact on the migration of injected fluid into fractured permeable media. Appropriate transport properties such as relative permeability and capillary pressure are essential for successful simulation and prediction of multi-phase flow in such systems. However, the lack of a thorough understanding of the dynamics governing immiscible displacement in fractured media, limits our ability to properly represent their macroscopic transport properties. Previous experimental observations of imbibition front evolution in fractured rocks are examined in the present study using an automated history-matching approach to obtain representative relative permeability and capillary pressure curves. Predicted imbibition front evolution under different flow conditions resulted in an excellent agreement with experimental observations. Sensitivity analyses, in combination with direct experimental observation, allowed exploring the competing effects of relative permeability and capillary pressure on the development of saturation distribution and imbibing front evolution in fractured porous media. Results show that residual saturations are most sensitive to matrix relative permeability to oil, while the ratio of oil and water relative permeability, rock heterogeneity, boundary condition, and matrix–fracture capillary pressure contrast, affect displacement shape, speed, and geometry of the imbibing front.     

Attenuation of strain waves in core samples of three types of rock

Dynamic photoelasticity was employed to determine the velocity of longitudinal stress waves, dynamic modulus of elasticity and attenuation coefficients in rockcore samples 1 in. (25 mm) in diameter, 18 in. (0.46 m) long. Birefringent strips bonded to the core samples of Salem limestone, Charcoal granite and Berea sandstone provided all the data needed for the dynamic characterization of these rock types. The rods were dynamically loaded at one end with a lead-azide charge. A multiple-spark-gap camer was used to record the dynamic isochromatic-fringe patterns occurring in the birefringent strip. Of the three rock types investigated, the Berea sandstone exhibited the largest energy losses as characterized by an attenuation coefficient of 0.0910. Salem limestone and Charcoal granite exhibited much smaller losses with attenuation coefficients of 0.0196 and 0.0024, respectively. The extremely low-energy loss associated with Charcoal granite indicates that this material transmits stress waves as well as most metals.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

贵州关岭大寨崩滑碎屑流灾害初步研究

2010年6月28日,贵州省关岭县岗乌镇大寨村发生特大型崩滑碎屑(石)流灾害,造成99人死亡或失踪。通过现场考察崩滑区的地质环境与斜坡岩体结构,认为斜坡体由似"干砌块石结构"的裂隙化岩体组成是发生崩溃式破坏的主要内在原因。超常暴雨(过程雨量237mm)条件下斜坡岩体后缘裂缝充水形成持续的"水楔作用"是斜坡岩体松动、倾倒垮塌的主要外部引发因素。碎屑(石)流块度的空间分布具有从源头向沟口逐次减小,碎屑(石)流运动冲击高度逐步降低,冲击速度逐步减小,并显示4个能级4个冲程的特点。根据动势能守恒定律,计算了每个冲程的最大速度,得出第1冲程为高速崩滑,其它冲程属于碎屑(石)流动冲击。未发现区域天然地震、光照水库诱发地震与外围历史采矿活动与本次事件相关的直接证据。由于滑坡后缘仍存在不稳定岩体,碎屑(石)流堆积体上多处分布直径3~5m的堰塞塘,说明碎屑(石)空隙的排泄能力不足,在未来暴雨条件下引发新的崩滑或形成沟谷型泥石流的可能性是存在的。