汶川地震断层岩气体和液体渗透率实验研究

以N2和水为孔隙流体在20-180 MPa围压范围内详细测量了汶川地震断裂带断层岩的渗透率,并同时测量了作为参照的砂岩的渗透率.实验结果表明,气体和液体渗透率均随围压增加而幂次衰减.在相同围压和孔隙压条件下,N2渗透率高于水渗透率接近1个数量级.将Klinkenberg效应校正后的气体渗透率（绝对渗透率）与水渗透率进行了比较.对比结果显示,砂岩的绝对渗透率与水渗透率基本一致;但断层岩绝对渗透率显著高于水渗透率,表明断层岩气体渗透率与水渗透率之差别不能完全由Klinkenberg效应解释.对实验数据的拟合分析表明,滑脱因子b值与绝对渗透率kl存在幂律关系（b=λkld）.断层岩符合b=0.2×10–3kl–0.557关系（R2=0.998）,显示绝对渗透率在10^–16-10^–20 m2范围内,控制断层岩渗透率的因素是一致的.另一方面,断层岩的d值明显小于砂岩的d值,暗示二者的滑脱效应存在差异.分析表明,断层岩中的粘土矿物颗粒表面吸附水及粘土矿物吸水膨胀导致有效孔隙尺寸减小是造成断层岩渗透率显著低于绝对渗透率的主要原因.研究结果显示,对于富含粘土矿物的断层岩,经Klinkenberg校正后的气体渗透率与水渗透率并不一致,因此采用相应的液体作为测量介质才能够更为准确和真实地揭示地下流体的渗流状况.     

汶川地震断层带结构及渗透率

对汶川地震断层带进行了跨断层的渗透率测量.结果显示汶川地震断层由低渗的核部(2.4×10-19～3.8×10-16m2)、高渗的破碎带(3.7×10-16～3.0×10-15m2)以及含裂隙原岩(6.0×10-18～4.3×10-13 m2)组成(有效压力40 MPa),其中新鲜断层泥具有最低的渗透率.断层泥和两侧原岩...     

低渗砂岩储层渗透率各向异性规律的实验研究

应用1摩尔/升的盐水作为孔隙介质,采用液体压力脉冲法对采自鄂尔多斯盆地某油田三叠系延长组低渗砂岩样品进行了渗透率随有效应力的变化规律实验研究.在围压0～100MPa,孔隙压力0～12MPa范围内,样品的渗透率变化范围在0～60×10-18 m2之间,实验结果表明样品的渗透率随有效应力的增加而减少.通过拟合实验结果,得到渗透率随有效应力的变化规律较好的符合幂函数关系,相关系数介于0.903～0.984之间.同时对同一样品的X,Y,Z相互垂直三个方向的渗透率随有效应力的变化进行了比较,结果表明渗透率的各向异性同样是压力的函数,且随着有效应力的增加,不同平面内的渗透率各向异性表现出了不同的变化规律.但在实验压力范围内,渗透率各向异性皆为正值,表明在有效应力为100MPa范围内孔隙内流体的流动方向未发生改变.研究结果为鄂尔多斯盆地低渗(超低渗)油气田开发,特别是对深部油气田开发过程中开发方式的选择与设计,充分利用渗透率各向异性的特点,提高采收率提供了新的岩石物性资料.     

超低渗透率测量仪的测试标定及初步测量结果

低渗-超低渗透率测量是流体渗流力学和岩石物理学研究领域的热点问题.断层带流体输运的定量化研究对于了解断层的力学性质和预测地下流体的流动具有重要意义.文中介绍了新研制的超低渗透率测量系统和工作原理.这套设备使用稳态法和目前逐渐流行的孔隙压力振荡法进行测量,其孔隙流体为蒸馏水和氮气,最高围压可达200M Pa,孔隙压上限为...     

水库蓄水与断层带流体孔隙压时-空分布的数值模拟——以紫坪铺水库为例

通过对汶川地震破裂带上的断层岩及围岩的渗透率测量,获得了研究区岩石的渗透率随深度变化的规律。依据所得到的渗透率实验数据,数值模拟了紫坪铺水库蓄水所产生的流体孔隙压的时空分布。模拟结果显示,渗透率是否随深度变化,其流体孔隙压的模拟结果也相差很大。由于断层带具有高渗的特点,因此与水库底部存在水力联系的断层带及其宽度对流体渗透的模拟结果有明显的影响。断层带越宽,下渗作用越明显,流体孔隙压的分布与断层带的产状越接近。若断层带宽度为100~300m,2008年汶川MS8.0地震发生时,其震源区的流体孔隙压为0.1~0.15MPa,与库体载荷在断层面上产生的附加正应力相当。该模拟结果可以作为进一步分析汶川地震的发生是否与紫坪铺水库蓄水存在关联的重要约束。     

汶川地震断裂带水岩相互作用及其对断裂带演化影响

汶川地震断层岩的矿物学和地球化学特征揭示出地震断层经历了漫长时间演化和复杂的水岩相互作用.间震期水岩相互作用导致断层岩中的破碎矿物蚀变,尤其是长石等矿物含量渐渐减少甚至消失,而黏土矿物（蒙脱石、伊利石、伊/蒙混层、绿泥石等）含量逐渐增高,以及如黄铁矿、石膏、重晶石、坡缕石等热液系统中常见的矿物大量出现;Mg、P、Ti、Mn、Fe等元素倾向富集在断层带中,而Si、K和Na等出现明显的亏损;元素的大量迁移导致断层带的体积巨量亏损.实验结果表明,黏土矿物的亲水性引起水渗透率比干燥气体渗透率明显偏低,并且二者偏差无法通过Klinkenberg校正消除.蒙脱石吸水膨胀和黏土矿物颗粒表面吸附孔隙流体造成孔隙度降低是导致水渗透率偏低的重要原因.断层岩碎屑结构使得其中的孔隙可能在600 MPa围压下得以保存,从而有助于流体沿断层带下渗,并在断层带深部形成高流体孔隙压.地震断层的主要矿物学及粒度分布特征并非在地震破裂过程中形成,因此利用断层岩粒度分布资料估算地震破裂能并不合适.     

汶川地震断层带传输性质研究——对地震同震弱化作用的启示

正断层带的流体压力因关系着地震断层滑动强度以及应力释放过程一直倍受重视。断层带渗透率结构对于了解断层带内流体活动状况、断层强度和摩擦稳定性以及同震热压作用等至关重要。本文全面总结了断层带的渗透性,尤其是大的地震断层带渗透性的研究成果。为了开展断层岩传输性质的测量工作,开发了超低渗岩石渗流测量系统。以汶川地震发震断层映秀—北川断层带上两个典型露头剖面为对象,分析了断层带的内部结构,对断层岩进行了成     

汶川地震破裂带断层岩纵波速度与孔隙度关系的实验研究

分布于地震破裂带上的断层岩具有高孔隙度的特征.该特点造成了其弹性波速度与结晶岩石和沉积岩存在明显的差异.确定断层岩的弹性波速度与孔隙度和矿物组成的关系对于利用地震资料探测深部断层和测井资料的解释至关重要.在10~600 MPa条件下,本文对地震断层岩的纵波波速（V_p）和总孔隙度（φ_t）进行了测量,并深入分析了V_p与孔隙度的关系.结果表明在10~600 MPa的压力范围内,V_p（p）随着压力的增高呈现对数增加,其增长率随着压力的上升而逐渐减小,遵从∂V_p（p）/∂p=a_v/p的变化规律.断层岩中的孔隙度随着压力的增高呈对数减小.与传统的认识不同,实验发现在压力高达600 MPa以上,大多数断层岩中仍然可以残留可观的孔隙量.分析显示V_p与总孔隙度及总粘土含量呈负线性相关.该发现有助于认识深部流体的活动通道特征,有助于理解断层带中存在大量粘土矿物、断层带内的物质可被大量带出、围陷波的形成等地质和地球物理现象.     

不同流体与黏土矿物对岩石渗透率影响实验研究

储层的渗透率是描述流体在岩石中渗透性的固有属性,在实验室测量的样品液体渗透率往往小于克式校正法修正后的Klinkenberg渗透率.本研究分别以氮气和盐水为介质,分别测量了 30块不含黏土矿物的枫丹白露砂岩试块与30块含黏土矿物的Y砂岩试块的渗透率,用以研究不同流体与黏土矿物对渗透率的影响以及气体测量渗透率时滑脱因子b...     

四川安岳地区须家河组二段低渗储层不同孔隙结构岩电参数变化规律

a、b、m、n、F、I是阿尔奇公式中重要的岩电参数,在测井解释中具有重要作用.近年的研究发现,在低渗储层中,孔隙结构直接影响储层产能评价的准确性.然而在目前的低渗储层岩电实验研究中,仅仅对几个岩电参数进行分析,而没有从孔隙结构的角度全面和系统的研究a、b、m、n、F、I六个岩电参数及其变化规律.为此,在充分利用高压隔板岩电实验的基础上,结合压汞测试、核磁共振等资料,对安岳地区须二段低渗储层不同孔隙结构的岩电参数变化规律及其与产能的关系进行分析.结果表明,研究区Ⅰ~Ⅳ四类孔隙结构对应的a、b、m、n、F、I各岩电参数的变化范围不同;Ⅰ类到Ⅳ类孔隙结构的a、n、F变化范围值递增,b、m变化范围值递减.低渗储层的产能与孔隙结构类型密切相关,且与岩电参数的变化范围具有对应性.以上不同孔隙结构岩电参数变化规律的分析为准确理解和应用实际岩电参数提供了依据.     

Evidence for extreme variations in the permeability of laterite from a detailed analysis of well behaviour in Nigeria

Laterite soils are widespread in tropical Africa and have a large impact on the hydrology of the areas they cover. The permeability of laterite helps determine the partitioning of runoff and interflow and regulates groundwater recharge to underlying bedrock. Groundwater within laterite also forms a widespread source of drinking water, typically from unimproved hang‐dug‐wells. Despite its importance, there is little published information on laterite aquifer properties. In this study, data from a 6 m deep well in Nigeria have been analysed to characterise the hydraulic conductivity of the laterite from repeated pumping tests. Transmissivity measurements from 40 tests spread out across a hydrological year varied from 0.1 to 1000 m²/d. Further interpretation of the data demonstrate a strong non‐linear decrease in horizontal hydraulic conductivity with depth, characterised by an upper horizon of extreme permeability (400 m/d), and a much lower permeability profile beneath (<0.1 m/d). These data are substantiated with observations from other wells throughout the area. This non‐linear permeability structure has several implications: the upper laterite can facilitate rapid lateral throughflow in the wet season, enabling contaminants to be transported significant distances (up to 1 km); natural groundwater levels are restricted to a narrow range for much of the year; and, in the dry season, the lower permeability of the deeper laterite restricts the amount of water which can be abstracted from shallow wells, leading to well failure. The work highlights the need for a wider study to better understand laterite soils and the role they play in regional hydrology. © 2013 Natural Environment Research Council. Hydrological Processes published by John Wiley & Sons Ltd.     

Comparison of the effective pressure law for permeability and resistivity formation factor in Chelmsford granite

Permeability, resistivity formation factor, and pore volume change were simultaneously measured on samples of Chelmsford granite subjected to confining pressure and pore pressure cycles. Using a technique described in a previous paper, the tangent coefficients of the effective pressure law for permeability _k and for formation factor _F were determined. _k and _F did not differ significantly from one another. They showed a strong stress history dependence as has already been observed for _k in several crystalline rocks. According to the definition of the effective pressure law used here, two physical properties with identical 's must be related through a one-to-one functional relationship. Hence, the observation above suggests that such a relationship may be empirically found between permeability and formation factor. Indeed, analysis of the data revealed that, to a good approximation, permeability was inversely proportional to the formation factor. The same relation has previously been observed in other crystalline rocks. This relationship was included in a recent version of the so-called equivalent channel model. Using this model, the specific surface area of the cracksA _c/V_S, the standard deviation of the distribution of asperities heightsh and the hydraulic radiusm _o were evaluated. The following values were respectively found: 850 cm^–1, 0.008 m and 0.14 m. The specific surface area of the cracks was independently estimated on micrographs of polished sections using a standard quantitative stereology method. The result was in good agreement with the values estimated from the transport properties data.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

Scales of rock permeability

Permeability is a transport property which is currently measured in Darcy units. Although this unit is very convenient for most purposes, its use prevents from recognizing that permeability has units of length squared. Physically, the square root of permeability can thus be seen as a characteristic length or a characteristic pore size. At the laboratory scale, the identification of this characteristic length is a good example of how experimental measurements and theoretical modelling can be integrated. Three distinct identifications are of current use, relying on three different techniques: image analysis of thin sections, mercury porosimetry and nitrogen adsorption. In each case, one or several theoretical models allow us to derive permeability from the experimental data (equivalent channel models, statistical models, effective media models, percolation and network models). Permeability varies with pressure and temperature and this is a decisive point for any extrapolation to crustal conditions. As far as pressure is concerned, most of the effect is due to cracks and a model which does not incorporate this fact will miss its goal. Temperature induced modifications can be the result of several processes: thermal cracking (due to thermal expansion mismatch and anisotropy, or to fluid pressure build up), and pressure solution are the two main ones. Experimental data on pressure and temperature effects are difficult to obtain but they are urgently needed. Finally, an important issue is: up to which point are these small scale data and models relevant when considering formations at the oil reservoir scale, or at the crust scale? At larger scales the identification of the characteristic scale is also a major goal which is examined.     

Seismoelectric measurements in a porous quartz‐sand sample with anisotropic permeability

Seismoelectric coupling coefficients are difficult to predict theoretically because they depend on a large numbers of rock properties, including porosity, permeability, tortuosity, etc. The dependence of the coupling coefficient on rock properties such as permeability requires experimental data. In this study, we carry out a set of laboratory measurements to determine the dependence of seismoelectric coupling coefficient on permeability. We use both an artificial porous “sandstone” sample, with cracks, built using quartz‐sand and Berea sandstone samples. The artificial sample is a cube with 39% porosity. Its permeability levels are anisotropic: 14.7 D, 13.8 D, and 8.3 D in the x‐, y‐, and z‐directions, respectively. Seismoelectric measurements are performed in a water tank in the frequency range of 20 kHz–90 kHz. A piezoelectric P‐wave source is used to generate an acoustic wave that propagates through the sample from the three different (x, y, and z) directions. The amplitudes of the seismoelectric signal induced by the acoustic waves vary with the direction. The highest signal is in the direction of the highest permeability, and the lowest signal is in the direction of the lowest permeability. Since the porosity of the sample is constant, the results directly show the dependence of seismoelectric coefficients on permeability. Seismoelectric measurements with natural rocks are performed using Berea sandstone 500 and 100 samples. Because the Berea samples are nearly isotropic in permeability, the amplitudes of the seismoelectric signals induced in the different directions are the same within the measurement error. Because the permeability of Berea 500 is higher than that of Berea 100, the amplitude of the seismoelectric signals induced in Berea 500 is higher than those in Berea 100. To determine the relative contributions of porosity and permeability on seismoelectric conversion, we carried out an analysis, using Pride (1994) formulation and Kozeny–Carman relationship; the normalized amplitudes of seismoelectric coupling coefficients in three directions are calculated and compared with the experimental results. The results show that the seismoelectric conversion is related to permeability in the frequency range of measurements. This is an encouraging result since it opens the possibility of determining the permeability of a formation from seismoelectric measurements.     

Heat transfer coefficients of natural volcanic clasts

Heat transfer coefficients used in numerical simulations of volcanic eruptions are typically borrowed from industrial settings where the coefficients are well determined for non-permeable, machined (spherical) materials. Volcanic clasts, in contrast, are permeable and have irregular shapes. We performed a series of laboratory experiments to determine heat transfer coefficients for natural volcanic particles. We measured the surface and interior temperatures during cooling at wind speeds ranging from 0 to 10 m/s. We also measured the permeability and density of the particles. We find that the permeability of the particles has little effect on clast cooling. In the absence of any wind, heat loss occurs by free convection, and we find no relationship between the heat transfer coefficient and particle density. However, for non-zero Reynolds numbers (finite wind speed), the heat transfer coefficient decreases with increasing porosity. We obtain a correlation for the dimensionless heat loss, or Nusselt number, of the form Nu = 2 + aRe^1/2Pr^1/3 where a is a density dependent coefficient given by a = 0.00022ρ + 0.31, with ρ in kg/m³, and Re and Pr are the Reynolds number and Prandtl number, respectively. Compared with non-porous particles, heat transfer coefficients for natural pumice clasts are reduced by a factor of 2–3 for particles with similar Re. Numerical simulations show that this leads to an increase in depositional temperature by 50–90 °C.     

Threshold pressure gradient of fluid flow through multi-porous media in low and extra-low permeability reservoirs

After analyzing many studies of fluid flow theory of multi-porous media in low and extra-low permeability reservoirs and the numerical simulation of non-Darcy flow, we found that a negative flow rate occurs in the existing non-Darcy flow equation, which is unreasonable. We believe that the existing equation can only be considered as a discriminant to judging Darcy flow or non-Darcy flow, and cannot be taken as a fluid flow governing equation of multi-porous media. Our analysis of the experimental results shows that the threshold pressure gradient (TPG) of low and extra-low permeability reservoirs is excessively high, and does not conform to fluid flow through multi-porous media in the actual reservoir situation. Therefore, we present a reasonable TPG ranging from 0.006 to 0.04 MPa/m at the well depth of 1500 m and oil drainage distance of 500 m. The results of our study also indicate that the non-Darcy flow phenomenon will disappear when the TPG reaches a certain value. In addition, the TPG or non-Darcy flow in low and extra-low permeability reservoirs does not need to be considered in the productivity prediction and reservoir numerical simulation. At present, the black oil model or dual-porous media is suitable for simulating low and extra-low permeability reservoirs.     

Influence of geostatic stresses on permeability of jointed rock masses

Ｉｎｆｌｕｅｎｃｅｏｆｇｅｏｓｔａｔｉｃｓｔｒｅｓｅｓｏｎｐｅｒｍｅａｂｉｌｉ┐ｔｙｏｆｊｏｉｎｔｅｄｒｏｃｋｍａｓｅｓＣＵＡＮＧ－ＢＢＩＮＧＺＨＯＵ（周创兵）ａｎｄＷＥＮ－ＬＩＮＸＩＯＮＧ（熊文林）ＷｕｈａｎＵｎｉｕｅｒｓｉｔｙｏｆＨｙｄｒａｕｌ...     

Fluid transport by solitary waves along growing faults: A field example from the South Eugene Island Basin, Gulf of Mexico

The Red Fault system is one of the main growth faults found in the South Eugene Island Basin, a salt withdrawal minibasin located offshore Louisiana, in the Gulf of Mexico. This fault system corresponds to a lateral boundary between fluid overpressured compartments. In addition, there is a set of observations indicating that the Red Fault system exhibits rapid episodic migration of fluids. This fault represents an example of preferential pathway for the upward episodic migration of overpressured hydrocarbons from deep, heavily pressured, compartments on time scales of years. The migrations of fluids into active growing faults could take the form of propagating surges (solitary waves) that propagate upward along the fault planes in a wave-like manner at km/yr. Solitary waves represent a very efficient mechanism for the upward transport of fluids along growth faults in sedimentary basins generating its own permeability. In addition, this mechanism is compatible with the fact that the fault plane is observed to sustain a static pore fluid pressure difference between its two sides. The propagation of solitary waves in active growth faults appears as a fundamental mechanism to understand the nature of upward fast migration of fluids along active growth faults in compartimentalized sedimentary basins.     

Permeability measurements of natural and experimental volcanic materials with a simple permeameter: Toward an understanding of magmatic degassing processes

Permeability measurement of quenched volcanic porous materials is an important approach to understand permeability development and degassing of vesicular silicic magmas. In this study, we developed a gas permeameter to measure permeability of natural samples and experimental products. The permeameter has broad measurement ranges of pressure difference (10¹–10⁵ Pa) and gas-flow rate (10^− 9–10^− 5 m³/s). These ranges enable us to measure viscous permeability in the range of 10^− 17–10^− 9 m² for 1 centimeter-scale samples, using the Forchheimer equation, which includes the inertial effect of gas flow permeating through the samples. In addition, we improved the procedure for performing permeability measurements of mm-sized products of decompression experiments. Although a previous study reported the first permeability data for vesicular silicic glass products of decompression experiments, we found an overestimation in their permeability data due to problems in sample preparation, especially for very low permeability samples. Our improved measurements give lower permeability values than those of Takeuchi et al. (2005)(Takeuchi, S., Nakashima, S., Tomiya, A., Shinohara, H., 2005. Experimental constraints on the low gas permeability of vesicular magma during decompression. Geophys. Res. Lett., 32, L10312 doi:10.1029/2005GL022491).