Axial and Radial Permeability Evolutions of Compressed Sandstones: End Effects and Shear-band Induced Permeability Anisotropy

The influence of hydrostatic and uniaxial stress states on the porosity and permeability of sandstones has been investigated. The experimental procedure uses a special triaxial cell which allows permeability measurements in the axial and radial directions. The core sleeve is equipped with two pressure samplers placed distant from the ends. They provide mid-length axial permeability measure as opposed to the overall permeability measure, which is based on the flow imposed through the pistons of the triaxial cell. The core sleeve is also equipped to perform flows in two directions transverse to the axis of the sample. Two independent measures of axial and complementary radial permeability are thus obtained. Both Fontainebleau sandstone specimens with a porosity of about 5.8% to 8% and low permeability ranging from 2.5 mD to 30 mD and Bentheimer sandstone with a porosity of 24% and a high permeability of 3D have been tested. The initial axial permeability values obtained by each method are in good agreement for the Fontainebleau sandstone. The Bentheimer sandstone samples present an axial mid-length permeability 1.6 times higher than the overall permeability. A similar discrepancy is also observed in the radial direction, also it relates essentially to the shape of flow lines induced by the radial flow. All the tested samples have shown a higher stress dependency of overall and radial permeability than mid-length permeability. The effect of compaction damage at the pistons/sample and radial ports/sample interfaces is discussed. The relevance of directional permeability measurements during continuous uniaxial compression loadings has been shown on the Bentheimer sandstone until the failure of the sample. We can efficiently measure the influence of brittle failure associated to dilatant regime on the permeability: It tends to increase in the failure propagation direction and to decrease strongly in the transverse direction.     

Instream flow requirements for sediment transport in the lower Weihe River

Estimation of instream flow requirements for sediment transport (IFRST) in the downstream Weihe River is important for the maintenance of a healthy ecosystem of the river. In this study, the IFRST for the lower Weihe River is estimated on the basis of the observed hydrological data during the period 1960–2001. The results showed that the mean annual IFRSTs for the reach between Lintong and Huaxian and the reach below Huaxian were 6·85 and 7·62 billion m³/year, respectively, and the standard errors were 0·50 and 0·76 billion m³/year, respectively. The results also showed that the Sanmenxia Reservoir is an important driver for the changes in channel morphology and hydraulics, as well as the IFRST. Furthermore, according to the hydrological frequency (p value) estimated from annual instream flow data during the period 1960–2001, four typical years (p = 25%, 50%, 75% and 90%) for the two reaches of the Weihe River were determined. The analysis showed that the IFRST has a negative power functional relationship with the sediment concentration. On the basis of the efficiency coefficients (R²) during three sub‐periods 1960–1973, 1974–1990 and 1991–2001, the annual IFRSTs in the above two reaches for the four typical years were estimated under different deposition and erosion–deposition conditions. The results provide useful references for restoration and water resource management of the Weihe River. Copyright © 2010 John Wiley & Sons, Ltd.     

Upper mantle P wave velocity structure of the northern part of China and Mongolia

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｅｔｅｒｍｉｎａｔｉｏｎｏｆｆｉｎｅｒａｄｉａｌｖｅｌｏｃｉｔｙｓｔｒｕｃｔｕｒｅｏｆｕｐｐｅｒｍａｎｔｌｅｐｌａｙｓａｎｉｍｐｏｒｔａｎｔｒｏｌｅｉｎｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｍａｎｔｌｅｃｏｍｐｏｓｉｔｉｏｎａ...     

Permeability Evolution During Non-linear Viscous Creep of Calcite Rocks

Permeability, storage capacity and volumetric strain were measured in situ during deformation of hot-pressed calcite aggregates containing 10, 20, and 30 wt% quartz. Both isostatic and conventional triaxial loading conditions were used. The tests were performed at confining pressure of 300 MPa, pore pressures between 50 to 290 MPa, temperatures from 673 to 873 K and strain rates of 3 × 10⁻⁵ s⁻¹. Argon gas was used as the pore fluid. The initial porosities of the starting samples varied from 5% to 9%, with higher porosity correlated to higher quartz content. Microstructural observations after the experiment indicate two kinds of pores are present: 1) Angular, crack-like pores along boundaries between quartz grains or between quartz and calcite grains and 2) equant and tubular voids within the calcite matrix. Under isostatic loading conditions, the compaction rate covaries with porosity and increases with increasing effective pressure. Most of the permeability reduction induced during compaction is irreversible and probably owes to plastic processes. As has been found in previous studies on hot-pressed calcite aggregates, permeability, k, is nonlinearly related to porosity, ϕ. Over small changes in porosity, the two parameters are approximately related as k ∝ ϕⁿ. The exponent n strongly increases as porosity decreases to a finite value (from about 4 to 6% depending on quartz content), suggesting a porosity percolation threshold. When subjected to triaxial deformation, the calcite-quartz aggregates exhibit shear-enhanced compaction, but permeability does not decrease as rapidly as it does under isostatic conditions. During triaxial compaction the exponent n only varies between 2 and 3. Non-isostatic deformation seems to reduce the percolation threshold, and, in fact, enhances the permeability relative to that at the same porosity during isostatic compaction. Our data provide constraints on the governing parameters of the compaction theory which describes fluid flow through a viscous matrix, and may have important implications for expulsion of sedimentary fluids, for fluid flow during deformation and metamorphism, and melt extraction from partially molten rocks.     

Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

Application of <Superscript>137</Superscript>Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau,China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of reservoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq kg-1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq kg-1 and 0.92 Bq kg-1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t km-2 a-1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t km-2 a-1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t km-2 a-1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low‐shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional‐wave velocity is strongly correlated with the shear‐wave velocity in this suite of rocks. The ratio V_P/V_S varies significantly with change of both pore‐fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot‐flow and squirt‐flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt‐flow of fluid between the pores of various geometries. The contribution of the Biot‐flow loss mechanism to the total loss is negligible. The compressional‐wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low‐shale, high‐porosity sandstone reservoirs.     

Shear-enhanced compaction during non-linear viscous creep of porous calcite-quartz aggregates

In this paper, we extend the previous studies of semi-brittle flow of synthetic calcite-quartz aggregates to a range of temperatures and effective pressures where viscous creep occurs. Triaxial deformation experiments were performed on hot-pressed calcite-quartz aggregates containing 10, 20 and 30 wt% quartz at confining pressure of 300 MPa, pore pressures of 50-290 MPa, temperatures of 673-1073 K and strain rates of 3.0×10⁻⁵/s, 8.3×10⁻⁵/s and 3.0×10⁻⁴/s. Starting porosity varied from 5 to 9%. We made axial and volumetric strain measurements during the mechanical tests. Pore volume change was measured by monitoring the volume of pore fluid that flows out of or into the specimen at constant pore pressure. Yield stress increased with decreasing porosity and showed a dependence on effective pressure. Thus, the yield stress versus effective pressure can be described as a yield surface with negative slope that expands with decreasing porosity and increasing strain hardening, gradually approaching the envelope of strength at 10% strain, which has a positive slope. Creep of porous rock can be modeled to first order as an isolated equivalent void in an incompressible nonlinear viscous matrix. An incremental method is used to calculate the stress-strain curve of the porous material under a constant external strain rate. The numerical simulations reproduce general trends of the deformation behavior of the porous rock, such as the yield stress decreasing with increasing effective pressure and significant strain hardening at high effective pressure. The drop of yield stress with increasing porosity is modeled well, and so is the volumetric strain rate, which increases with increasing porosity.     

Transect study on solute transport in a macroporous soil

Solute transport experiments using a non-reactive tracer were conducted on short, undisturbed, saturated columns of a sandy loam soil. All columns, 20 cm in diameter and 20 cm long, were collected along a transect of 35 m. Most of the soil columns had pre-existing macropores. The columns were leached at a steady flow-rate under ponding conditions. The resulting breakthrough curves (BTCs) showed a large heterogeneity. Several of the BTCs displayed early breakthrough and long tailing. All the data were interpreted in terms of dimensional time moments, the classical convection-dispersion equation (CDE) and the mobile-immobile transport model (MIM). Experimental time moments were found to vary significantly among the different BTCs. Analysis of the time moments also revealed that the variance of the field-scale BTC was several times larger than the average of the local-scale variance. The pore water velocity v and dispersion coefficient D were obtained by fitting the CDE to the local-scale BTCs, resulting in an average dispersivity of 7·4 cm. Frequency distributions for the CDE parameters v and D were equally well described by a normal or log-normal probability density function (pdf). When a log-normal pdf for D is considered, the variance of the log_e transformed D values (σ_{ln D}²) was found to be 2·1. For the MIM model, two additional parameters were fitted: the fraction of mobile water, θ_m/θ, and the first-order mass transfer coefficient, α. The MIM was more successful in describing the data than the CDE transport model. For the MIM model, the average dispersivity was about 2 cm. The MIM parameters v, D and θ_m/θ were best described by a log-normal pdf rather than a normal pdf. Only the parameter α was better described by a normal pdf. Mobile water fractions, θ_m/θ ranged from 0·01 to 0·98, with a mean of 0·43 (based on a log-normal pdf). When the CDE and MIM were applied to the data, the fitted pore water velocities, v, compared favourably with the effective pore water velocities, v_eff, obtained from moment analysis.     

Network flow modeling via lattice-Boltzmann based channel conductance

Lattice-Boltzmann (LB) computations of single phase, pore-to-pore conductance are compared to models in which such conductances are computed via standard pore body–channel–pore body series resistance (SR), with the conductance of each individual element (pore body, channel) based on geometric shape factor measurements. The LB computations, based upon actual channel geometry derived from X-ray computed tomographic imagery, reveal that the variation in conductance for channels having similar shape factor is much larger than is adequately captured by the geometric models. Fits to the dependence of median value of conductance versus shape factor from the LB-based computations show a power law dependence of higher power than that predicted by the geometric models. We introduce two network flow models based upon the LB conductance computations: one model is based upon LB computations for each pore-to-pore connection; the second is based upon a power law fit to the relationship between computed conductance and throat shape factor. Bulk absolute permeabilities for Fontainebleau sandstone images are computed using the SR-based network models and the two LB-based models. Both LB-based network models produce bulk absolute permeability values that fit published data more accurately than the SR-based models.     

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

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

Permeability of Triaxially Compressed Sandstone: Influence of Deformation and Strain-rate on Permeability

— The influence of differential stress on the permeability of a Lower Permian sandstone was investigated. Rock cylinders of 50 mm in diameter and 100 mm length of a fine-grained (mean grain size 0.2 mm), low-porosity (6–9%) sandstone were used to study the relation between differential stress, rock deformation, rock failure and hydraulic properties, with a focus on the changes of hydraulic properties in the pre-failure and failure region of triaxial rock deformation. The experiments were conducted at confining pressures up to 20 MPa, and axial force was controlled by lateral strain with a rate ranging from 10^?6 to 10^?7 sec^?1. While deforming the samples, permeability was determined by steady-state technique with a pressure gradient of 1 MPa over the specimen length and a fluid pressure level between 40 and 90% of the confining pressure. The results show that permeability of low-porosity sandstones under increasing triaxial stress firstly decreases due to compaction and starts to increase after the onset of dilatancy. This kind of permeability evolution is similar to that of crystalline rocks. A significant dependence of permeability evolution on strain rate was found. Comparison of permeability to volumetric strain demonstrates that the permeability increase after the onset of dilatancy is not sufficient to regain the initial permeability up to failure of the specimen. The initial permeability, which was determined in advance of the experiments, usually was regained in the post-failure region. After the onset of dilatancy, the permeability increase displays a linear dependence on volumetric strain.     

Estimates of effective permeability for non-Newtonian fluid flow in randomly heterogeneous porous media

An understanding of the interplay between non-Newtonian effects in porous media flow and field-scale domain heterogeneity is of great importance in several engineering and geological applications. Here we present a simplified approach to the derivation of an effective permeability for flow of a purely viscous power–law fluid with flow behavior index n in a randomly heterogeneous porous domain subject to a uniform pressure gradient. A standard form of the flow law generalizing the Darcy’s law to non-Newtonian fluids is adopted, with the permeability coefficient being the only source of randomness. The natural logarithm of the permeability is considered a spatially homogeneous and correlated Gaussian random field. Under the ergodic hypothesis, an effective permeability is first derived for two limit 1-D flow geometries: flow parallel to permeability variation (serial-type layers), and flow transverse to permeability variation (parallel-type layers). The effective permeability of a 2-D or 3-D isotropic domain is conjectured to be a power average of 1-D results, generalizing results valid for Newtonian fluids under the validity of Darcy’s law; the conjecture is validated comparing our results with previous literature findings. The conjecture is then extended, allowing the exponents of the power averaging to be functions of the flow behavior index. For Newtonian flow, novel expressions for the effective permeability reduce to those derived in the past. The effective permeability is shown to be a function of flow dimensionality, domain heterogeneity, and flow behavior index. The impact of heterogeneity is significant, especially for shear-thinning fluids with a low flow behavior index, which tend to exhibit channeling behavior.     

Laboratory Characterization of Permeability and Its Anisotropy of Chelungpu Fault Rocks

In Taiwan an international project to drill into the Chelungpu fault (TCDP) was initiated after the M _w 7.6 Chi-Chi earthquake in 1999. At Takeng, two vertical holes (A and B) to depths of about 2 km have been drilled through the northern portion of the Chelungpu fault system. In this study, we conducted systematic hydromechanical tests on TCDP drillcores collected from Hole-A at various depths above and below the major slip zone of the Chelungpu fault. We focus on the measurements of permeability as function of pressure and the brittle failure behavior. Evolution of permeability as a function of pressure and porosity was measured using either steady-state flow or a pulse transient technique. When subjected to an effective pressure reaching 100 MPa, permeability values of shaly siltstone samples range from 10^?16 to 10^?19 m². In comparison, permeability values of porous sandstones are at least an order of magnitude higher, ranging from 10^?14 to 10^?18 m². To characterize permeability anisotropy associated with the bedding structure of the rocks of the Chelungpu fault, cylindrical samples were taken from the TCDP drillcores along three orthogonal directions, denoted X, Y and Z respectively. Direction Z is parallel to the TCDP core axis, and the other two directions are perpendicular to the core axis, with X (N105°E) perpendicular and Y (N15°E) parallel to the strike of the bedding. In shaly siltstones, permeability values of samples cored along the strike of bedding (direction Y) can be up to 1 order of magnitude higher than those cored perpendicular to the strike of bedding (direction X). These observations indicate that permeability anisotropy is controlled by the spatial distribution of bedding in Chelungpu fault host rocks. Permeability evolution of fault rocks plays an important role in dynamic weakening processes, which are particularly pertinent to large earthquakes such as the Chi-Chi earthquake. Our experimental data on permeability and its anisotropy of TCDP core samples provide necessary constraints on fault models and proposed weakening mechanisms.     

Hydromechanical Heterogeneities of a Mature Fault Zone: Impacts on Fluid Flow

In this paper, fluid flow is examined for a mature strike‐slip fault zone with anisotropic permeability and internal heterogeneity. The hydraulic properties of the fault zone were first characterized in situ by microgeophysical (V_P and σ_c) and rock‐quality measurements (Q‐value) performed along a 50‐m long profile perpendicular to the fault zone. Then, the local hydrogeological context of the fault was modified to conduct a water‐injection test. The resulting fluid pressures and flow rates through the different fault‐zone compartments were then analyzed with a two‐phase fluid‐flow numerical simulation. Fault hydraulic properties estimated from the injection test signals were compared to the properties estimated from the multiscale geological approach. We found that (1) the microgeophysical measurements that we made yield valuable information on the porosity and the specific storage coefficient within the fault zone and (2) the Q‐value method highlights significant contrasts in permeability. Fault hydrodynamic behavior can be modeled by a permeability tensor rotation across the fault zone and by a storativity increase. The permeability tensor rotation is linked to the modification of the preexisting fracture properties and to the development of new fractures during the faulting process, whereas the storativity increase results from the development of micro‐ and macrofractures that lower the fault‐zone stiffness and allows an increased extension of the pore space within the fault damage zone. Finally, heterogeneities internal to the fault zones create complex patterns of fluid flow that reflect the connections of paths with contrasting properties.     

Belowground carbon balance and carbon accumulation rate in the successional series of monsoon evergreen broad-leaved forest

The belowground part of terrestrial ecosystem is a huge carbon pool. It is believed that of the total 2500Gt carbon stored in global terrestrial ecosystem, soil carbon storage within the 1 m surface layer ac- counts for 2000Gt, which is 4-fold of vegetation car- bon storage[1,2]. Compared with the carbon in the vegetation, carbon in the deep soil layers is much more stable, and it will stay in soil profile permanentlyunless geological vicissitude occurs. Essentially, forest restoration is the…     

A multiple-point statistics algorithm for 3D pore space reconstruction from 2D images

Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of long-range, curvilinear features of geological structures, serves as the simulation engine. Various validating criteria such as marginal distributions of pore and grain, directional variograms, multiple-point connectivity curves, single phase effective permeability and two phase relative permeability calculations are used to analyze the results. The method is tested on a sample of Berea sandstone for which a 3D micro-CT scanning image is available. The results confirm that the equi-probable 3D realizations obtained preserve the typical patterns of the pore space that exist in thin sections, reproduce the long-range connectivities, capture the characteristics of anisotropy in both horizontal and vertical directions and have single and two phase flow characteristics consistent with those of the measured 3D micro-CT image.     

The porosity trend and pore sizes of the rocks in the continental crust of the earth: Evidence from experimental data on permeability

The depth trends of permeability are constructed from the measurements of the tight rocks typical of the basement of the continental crust at temperatures up to 600°C and pressures up to 200 MPa. It is established that the permeability decreases with depth. The statistical processing of the experimental data yielded the generalized dependence logk = ?12.6-3.23H ^0.223. The method is suggested and, based on the experimental data on permeability, the estimates are obtained for the effective diameters of the pore channels and effective porosity at the PT parameters corresponding to the in situ deep zones of the continental crust. It is found that porosity decreases with depth, while distinct depth dependence of the pore sizes is not observed. The dependence of porosity on the depth is approximated by the relationship logφ = ?0.65–0.1H + 0.0019H ². The porosity is estimated at a few percent for a depth of 10 km with a decline to 0.01–0.1% at 35 km. The estimates of porosity retrieved from the experimental data agree with the theoretical calculations based on the present-day ideas of the structure of the discrete media and with the results of magnetotelluric sounding. Thus, according to three independent estimates, the porosity of the rocks of the continental crust decreases with depth. At the same time, in both the intermediate and lower crust there are intervals where the porosity values suggest the presence of fluid-saturated horizons at these depths.     

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra‐fragment porosity provides water storage and the inter‐fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0·3, 1·0, and 1·5 cm h^?1. Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel‐shaped cells each 3·5 cm in diameter. The total porosity of the columns was 0·47 ± 0·02 m³ m^?3, approximately 13% of pores were ≥ 15 µm diameter, and the saturated hydraulic conductivity was 12·66 ± 1·31 m day^?1. Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0·3, 1·0, 1·5 cm h^?1 application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C_o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two‐region convection–dispersion model was used to describe the BTCs and fitted well (r² = 0·97–0·99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0·95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright © 2007 John Wiley & Sons, Ltd.     

含不同孔隙流体的砂岩地震波速度随压力变化的实验研究

岩石的地震波性质是区域构造研究和浅部地震勘探应用的基础.延长油田是我国重要的油气生产基地之一,但目前仍缺乏地震波性质方面的基础资料.作者利用Autolab2000多功能岩石物性自动测试设备,在0~180MPa及饱含不同孔隙流体(干燥、饱水及饱油)条件下,研究了三种来自延长油田砂岩岩芯的纵波、横波速度.结果表明:三种砂岩的V_p、V_s1和V_s2均随压力增加(或降低)而基本呈对数曲线增大(或减小);干燥、饱水和饱油三种波速间的关系因砂岩类型不同而不同,这主要取决于岩石的有效弹性模量、孔隙流体性质以及岩石的内部结构等;含同种孔隙流体的不同类型砂岩,其V_p、V_s1和V_s2随压力变化的规律主要受岩石孔隙度和粒度的影响;而含不同孔隙流体的同种砂岩,其V_p、V_s1和V_s2随压力变化的规律则主要受控于岩石的有效弹性模量和流体密度.另外,含水或含油饱和度的变化对V_s1和V_s2基本没有影响.实验结果可以为该地区地震资料的解释及与声波测井之间的对比提供重要的基础数据和参考依据.