地层渗透率对成矿流体热场和流场的影响--以锡矿山锑矿床成矿流体为例

沉积盆地是流体活动最活跃的场所,而地层渗透率则是盆地流体运移的基础数据.沉积盆地中的热-重力驱动型流体的运移依赖于盆地的水文地质特征和成矿流体特征.为研究地层渗透率的改变对湖南锡矿山锑矿床成矿流体热场和流场的影响,首先选取一个对比地层渗透率组,计算出区域的热场和流场分布,然后改变导水地层渗透率的大小,并将其计算结果与对比渗透率组的热场结果进行对比,结果显示改变地层渗透率的大小对区域流场的影响较大,而对区域热场的影响则不大.在所讨论的地层渗透率组取值范围内,地层渗透率的改变对区域热场的影响介于5%～10%间,而地层渗透率的改变对区域流场的影响则达到2～3倍.     

Permeability semivariograms, geological structure, and flow performance

Clastic sediments may have a strong deterministic component to their permeability variation. This structure may be seen in the experimental semivariogram, but published geostatislical studies have not always exploited this feature during data analysis and covariance modeling. In this paper, we describe sedimentary organization, its importance for flow modeling, and how the semivariogram can be used for identification of structure. Clastic sedimentary structure occurs at several scales and is linked to the conditions of deposition. Lamination, bed, and bedset scales show repetitive and trend features that should be sampled carefully to assess the degree of organization and levels of heterogeneity. Interpretation of semivariograms is undertaken best with an appreciation of these geological units und how their features relate to the sampling program. Sampling at inappropriate intervals or with instruments having a large measurement volume, for example, may give misleading semivariograms. Flow simulations for models which include and ignore structure show that the repetitive features in permeability can change anisotropy and recovery performance significantly. If systematic variation is present, careful design of the permeability fields therefore is important particularly to preserve the structure effects.     

Stochastic representation for anisotropic permeability tensor random fields

In this paper, we introduce a novel stochastic model for the permeability tensor associated with stationary random porous media. In the light of recent works on mesoscale modeling of permeability, we first discuss the physical interpretation of the permeability tensor randomness. Subsequently, we propose a nonparametric prior probabilistic model for non‐Gaussian permeability tensor random fields, making use of the information theory and a maximum entropy procedure, and provide a physical interpretation of the model parameters. Finally, we demonstrate the capability of the considered class of random fields to generate higher levels of statistical fluctuations for selected stochastic principal permeabilities. This unique flexibility offered by the parameterization of the model opens up many new possibilities for both forward simulations (e.g. for uncertainty propagation in predictive simulations) and stochastic inverse problem solving. Copyright © 2011 John Wiley & Sons, Ltd.     

Review and classification of structural controls on fluid flow during regional metamorphism

Mechanisms for kilometre-scale, open-system fluid flow during regional metamorphism remain problematic. Debate also continues over the degree of fluid flow channellization during regional metamorphism, and the mechanisms for pervasive fluid flow at depth. The requirements for pervasive long-distance fluid flow are an interconnected porosity and a large regional gradient in fluid pressure and hydraulic head (thermally or structurally controlled) that dominates over local perturbations in hydraulic head due to deformation. In contrast, dynamic or transient porosity interconnection and fluid flow accompanying deformation of heterogeneous rock suites should result in moderately to strongly channellized flow at a range of scales, of which there are many examples in the literature. Classification of fluid flow types based on scale and degree of equilibration between fluid and rock, wallrock permeability, and mode of fluid transport contributes to an understanding of key factors that control fluid flow. Closed-system fluid behaviour, with restricted fluid flow in microcracks or cracks and limited fluid–rock interaction, occurs over a range of strains and crustal depths, but requires low permeabilities and/or small fluid fluxes. Long-distance, open-system fluid flow in channels is favoured in heterogeneous rocks at high strains, moderate (but variable) permeabilities, and moderate to high fluid fluxes. Long-distance, broad, pervasive fluid flow during regional metamorphism requires that the rocks are not accumulating high strains and have high permeabilities, low permeability contrasts, and high fluid fluxes. The ideal situation for such fluid flow is in situations where the rocks are undergoing stress relaxation immediately after a major deformation phase. In the mid-crust, fairly specific conditions are thus required for pervasive fluid flow. During active orogenesis, structurally controlled fluid flow (with focused open systems surrounding regions of closed-system behaviour) predominates in most, but not all, regional metamorphic situations, at a range of scales.     

Use of Border Regions for Improved Permeability Upscaling

A procedure for the improved calculation of upscaled grid block permeability tensors on Cartesian grids is described and applied. The method entails the use of a border region of fine-scale cells surrounding the coarse block for which the upscaled permeability is to be computed. The implementation allows for the use of full-tensor permeability fields on the fine and coarse scales. Either periodic or pressure–no flow boundary conditions are imposed over the extended local domain (target block plus border regions) though averaged quantities, used to compute the upscaled permeability tensor, are computed only over the target block region. Flow and transport results using this procedure are compared to those from standard methods for different types of geological and simulation models. Improvement using the new approach is consistently observed for the cases considered, though the degree of improvement varies for different models and flow quantities.     

Upscaling Permeability Measurements Within Complex Heterolithic Tidal Sandstones

We investigate numerically the effect of sample volume on the effective single-phase permeability of heterolithic tidal sandstones, using three-dimensional models reconstructed directly from large rock specimens measuring 45 × 30 × 15 cm. We find that both individual and averaged effective permeability values vary as a function of sample volume, which indicates that permeability data obtained from core-plugs will not be representative at the scale of a reservoir model grid-block regardless of the number of measurements taken. However, the error introduced by averaged data may be minimized using the appropriate averaging scheme for a given facies type and flow direction.     

Modeling permeability in imperfectly layered porous media. I. Derivation of block-scale permeability tensor for thin grid-blocks

Practical expressions are given for the nine components of the block-scale permeability tensor of a thin block. These expressions are derived from the local-scale continuity equation and Darcy's law in an anisotropic layered porous medium. The flow problem is separated in a bottom-flux problem and a top-flux problem, both of which can be solved in essentially the same way. The bottom-flux problem has been worked out in detail, and has been separated in two parts: a vertical potential difference and a horizontal potential difference part. Each is solved with a different approach specially designed for it. Depth-averaged expressions are obtained first, after which block-scale expressions are obtained by assuming a constant depth-averaged flux. In the zeroth order, this results in the well-known Dupuit approximation in geohydrology, and the vertical equilibrium (VE) approximation in petroleum reservoir engineering. The novelty of the theory presented here stems from the application of a perturbation technique to obtain first-order corrections to these well-known results. The local-scale laws are applied in the coordinate system coinciding with the principal axes of the local-scale permeability tensor. Only in this coordinate system the local-scale permeability tensor has zero off-diagonal components. However, since the porous medium is imperfectly layered, the first-order corrections show that the off-diagonal components of the block-scale permeability tensor are not zero. Furthermore, the block-scale permeability tensor is generally nonsymmetric, which implies that a coordinate system in which the off-diagonal terms disappear does not exist.     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

Modeling Hysteresis in Porous Media Flow via Relaxation

Two-phase flow in a porous medium can be modeled, using Darcy's law, in terms of the relative permeability functions of the two fluids (say, oil and water). The relative permeabilities generally depend not only on the fluid saturations but also on the direction in which the saturations are changing. During water injection, for example, the relative oil permeability k _ro falls gradually until a threshold is reached, at which stage the k _ro begins to decrease sharply. The latter stage is termed imbibition. If oil is subsequently injected, then k _ro does not recover along the imbibition path, but rather increases only gradually until another threshold is reached, whereupon it rises sharply. This second stage is called drainage, and the type of flow that occurs between the imbibition and drainage stages is called scanning flow. Changes in permeability during scanning flow are approximately reversible, whereas changes during drainage and imbibition are irreversible. Thus there is hysteresis, or memory, exhibited by the two-phase flow in the porous medium. In this work, we describe two models of permeability hysteresis. Common to both models is that the scanning flow regime is modeled with a family of curves along which the flow is reversible. In the Scanning Hysteresis Model (SHM), the scanning curves are bounded by two curves, the drainage and imbibition curves, where the flow can only occur in a specific direction. The SHM is a heuristic model consistent with experiments, but it does not have a nice mathematical specification. For instance, the algorithm for constructing solutions of Riemann problems involves several ad hoc assumptions. The Scanning Hysteresis Model with Relaxation (SHMR) augments the SHM by (a) allowing the scanning flow to extend beyond the drainage and imbibition curves and (b) treating these two curves merely as attractors of states outside the scanning region. The attraction, or relaxation, occurs on a time scale that corresponds to the redistribution of phases within the pores of the medium driven by capillary forces. By means of a formal Chapman–Enskog expansion, we show that the SHM with additional viscosity arises from the SHMR in the limit of vanishing relaxation time, provided that the diffusion associated with capillarity exceeds that induced by relaxation. Moreover, through a rigorous study of traveling waves in the SHMR, we show that the shock waves used to solve Riemann problems in the SHM are precisely those that have diffusive profiles. Thus the analysis of the SHMR justifies the SHM model. Simulations based on a simple numerical method for the simulation of flow with hysteresis confirm our analysis.     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体（分辨率3.7 μm）为研究对象，在空间节点数多达64 000 000的情况下，基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数，并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现，结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD（即10^-13m²）量级，远低于规则堆叠的圆球结构（>20 000 mD）；土壤的迂曲度为1.40~1.60，明显高于规则堆叠的圆球结构。研究结果表明，渗透率大的样品具有较小的迂曲度，这与结构的空间异质性有较强的关系；土壤的渗透率和迂曲度呈现各向异性；在水力梯度一定的前提下，渗透率较大的样品，纵向弥散系数也较大；同时，结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟，可用于土壤多孔介质在孔隙尺度下的水力学特性研究。     

利用相控模型进行井间参数预测

储层物性参数在空间的展布，一直是储层非均质性研究的难点，利用分形和蒙特卡洛方法，在精细的沉积微相及小层对比研究基础之上，预测了井间的岩石相分布，并根据取心井资料，建立了各种岩石相的渗透率相控模型，从而预测出井间渗透率的分布。该方法以地质特征为基础，将地质资料与数学方法相结合，充分体现了渗透率在井间既具有连续性，同时变化又非常剧烈的特点，为储层物性空间预测提供了新的思路。     

基于粒度控制的复杂储层渗透性建模方法

岩石颗粒大小反映了沉积岩结构的主要特征,是岩石岩性的主要评价指标,同时也是影响沉积物渗透性的主要因素。沉积环境在纵向上表现出来的差异性,使得粒序控制下的地层渗透性在纵向上也表现出非均质性。为了精细地描述由于粒度大小造成的渗透率的非均质性,本文首先结合岩心粒度分析资料,分析了研究区粒度大小对物性影响的规律,利用自然伽马、中子、密度、声波曲线,通过测井地质分析提取敏感参数,建立了粒度计算模型;然后以粒度为指标,建立了不同粒度大小条件下的渗透率计算模型。通过现场资料的处理及其与岩心资料的对比,认为这种渗透率建模方法能够较好地反映沉积物粒度对岩石渗透性的控制作用,能够更加精确地描述复杂岩性地层渗透率的非均质性特征,从而可为地质开发研究提供更为准确的渗透率剖面。     

About Darcy's law in non-Galilean frame

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupling between anisotropic damage and permeability variation in brittle rocks

In this paper, a coupled constitutive model is proposed for anisotropic damage and permeability variation in brittle rocks under deviatoric compressive stresses. The formulation of the model is based on experimental evidences and main physical mechanisms involved in the scale of microcracks are taken into account. The proposed model is expressed in the macroscopic framework and can be easily implemented for engineering application. The macroscopic free enthalpy of cracked solid is first determined by approximating crack distribution by a second‐order damage tensor. The effective elastic properties of damaged material are then derived from the free enthalpy function. The damage evolution is related to the crack growth in multiple orientations. A pragmatic approach inspired from fracture mechanics is used for the formulation of the crack propagation criterion. Compressive stress induced crack opening is taken into account and leads to macroscopic volumetric dilatancy and permeability variation. The overall permeability tensor of cracked material is determined using a micro–macro averaging procedure. Darcy's law is used for fluid flow at the macroscopic scale whereas laminar flow is assumed at the microcrack scale. Hydraulic connectivity of cracks increases with crack growth. The proposed model is applied to the Lac du Bonnet granite. Generally, good agreement is observed between numerical simulations and experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

Anisotropic Variation Law of Rock Permeability with the Burial Depth of Limestone

Permeability tensors of both macrofracture and microfracture systems were measured progressively along the depth of limestone formations at severed sites. It was found that the principal permeability values KX, Ky and Kz in these permeability tensors all decrease simultaneously and logarithmically with depth. However, the limestone aquifers are composed of an upper region where the larger permeability ellipsoid is upright or prolate and characterized by Kz>Kx and Kz>Ky, a transitional zone, and a lower zone whose smaller permeability ellipsoid is horizontal or oblate and characterized by KZ相似文献   

A model of regional fluid flow: Sand concentration factors and effective lateral and vertical permeabilities

Permeabilities of sand and shale differ by five to seven orders of magnitude. In the depositional pattern characteristic of the Gulf Coast, lenticular sand bodies occur randomly dispersed within a shale matrix, even though they can be grouped into stratigraphic sequences, contributing to net permeability of the sand-shale system and thus facilitating upward flow of pore water. In addition, a substantial lateral flow develops as fluid flow is directed toward sand bodies, thus concentrating fluid from a wide source area into sand bodies of limited lateral and vertical extent. Calculations reported here determine the effect of randomly distributed sands on both vertical and lateral permeabilities of composite sandshale systems. Effective permeabilities depend only on mean sand percent and local variance about the mean, because the effects of sand body shapes and their orientation and dip are not incorporated in the analysis. Results indicate that at a mean sand percentage of 20–50% a large lateral migration of fluid flow develops as sand bodies attract flow from a region two to seven times the radius of the sand body itself. For a mean sand percentage in excess of 50–60%, large vertical migration results and thus hydrocarbon trapping should not be expected for sand percentages in excess of 50–60%. Results of these calculations are in accord with observations that most larger oil accumulations on the Gulf Coast are found in regions having sandiness between about 20 and 50%.     

周期振荡法在低渗透测量中的应用研究

周期振荡法是一种新兴的多孔介质渗透率测量方法，具有测量周期短、稳定性好、精度高等优点。然而，由于缺乏对周期振荡法的系统研究，该方法未能在低渗透测量领域得到推广应用。采用数值模拟方法，分别研究了不同孔隙率、渗透率对周期振荡法气压测量结果的影响，以及不同形式周期波在周期振荡法渗透率测量中的适用性。结果表明：无论是孔隙率，还是渗透率都会影响气体压力传递过程；周期振荡法渗透率测量时可以不需要初始气压平衡过程，一定时间后测量结果几乎不受初始气压状态影响；与正弦波相比，方波对下游气体压力响应的影响更明显，而且方波、三角波和锯齿波具有容易加载的特点，建议采用方波取代正弦波进行周期振荡法试验。研究成果对指导周期振荡法试验和解决低渗透测量难题有一定帮助。     

Nodal and Mixed Finite Elements for the Numerical Homogenization of 3D Permeability

The aim of upscaling is to determine equivalent homogeneous parameters at a coarse-scale from a spatially oscillating fine-scale parameter distribution. To be able to use a limited number of relatively large grid-blocks in numerical oil reservoir simulators or groundwater models, upscaling of the permeability is frequently applied. The spatial fine-scale permeability distribution is generally obtained from geological and geostatistical models. After upscaling, the coarse-scale permeabilities are incorporated in the relatively large grid-blocks of the numerical model. If the porous rock may be approximated as a periodic medium, upscaling can be performed by the method of homogenization. In this paper the homogenization is performed numerically, which gives rise to an approximation error. The complementarity between two different numerical methods – the conformal-nodal finite element method and the mixed-hybrid finite element method – has been used to quantify this error. These two methods yield respectively upper and lower bounds for the eigenvalues of the coarse-scale permeability tensor. Results of 3D numerical experiments are shown, both for the far field and around wells.     

Physical properties of Mesozoic sedimentary rocks from the Perth Basin,Western Australia

The Perth Basin (PB) hosts important aquifers within the Yarragadee Formation and adjacent geological formations with potential for economic exploitation by both geothermal energy and carbon capture and sequestration. Published studies on the reservoir quality of the sedimentary units of the PB are very few. This study reports some petrophysical and lithological characteristics of the sedimentary units of interest for geothermal and geosequestration scenarios and help interpolation toward non-sampled intervals. A new fluvial-dominated lithofacies scheme was developed for the Mesozoic stratigraphy from four wells drilled in the central PB (Pinjarra-1, Cockburn-1, Gingin-1 and Gingin-2) based on grainsize, sorting, sedimentary structures and colour that relate to the environment of deposition. Systematic laboratory measurements of permeability, porosity, and thermal conductivity were conducted on core samples to investigate a variety of lithofacies and depths from these wells. Empirical correlations are established among the different physical properties, indicating encouraging relationships for full PB basin interpolation such as between porosity and permeability, when the samples are grouped into ‘hydraulic units’ defined by a ‘flow zone indicator’ parameter. The common principal controls on the PB thermal conductivity are the pore space arrangement and mineralogical content, which are strongly lithofacies-specific. Therefore, the lithofacies type could be a good first-order discriminator for describing spatial variations of thermal conductivity and then estimate their flow zone indicator.     

Modeling permeability in imperfectly layered porous media. II. A two-dimensional application of block-scale permeability

In the previous paper (Zijl and Stam, 1992), a theory has been developed to calculate the nine components of the three-dimensional intrinsic permeability tensor on the scale of a grid-block from a local-scale, predominantly layered subsurface. The resulting block-scale expressions can be written as a perturbation series of which the first term, or zeroth-order solution, coincides with the conventionally applied arithmetic and harmonic averages over the layers of the subsurface. The derived expressions permit the calculation of the diagonal and off-diagonal terms of the permeability tensor. In the present paper, these expressions will be applied in some numerical examples. Two basic two-dimensional hypothetical permeability distributions are adopted, and the various terms of the theoretical expressions are calculated. The results will be used to derive guidelines to discern the situations where higher order solutions can be neglected, and where conventional harmonic and arithmetic averages give a good estimate of the permeability on grid-block scale.