岩心孔隙介质中流体的核磁共振弛豫

弛豫时间是核磁共振研究中的一个重要参数,岩心孔隙介质流体的弛豫过程是自由流体弛豫机制、表面弛豫机制和流体的扩散弛豫机制共同作用的结果,它包含了丰富的孔隙和流体本身的信息. 弛豫时间和自扩散系数的测量及对弛豫时间的分析是核磁共振技术应用于岩心分析和石油勘测的重要内容.     

The internal magnetic field distribution, and single exponential magnetic resonance free induction decay, in rocks

When fluid saturated porous media are subjected to an applied uniform magnetic field, an internal magnetic field, inside the pore space, is induced due to magnetic susceptibility differences between the pore-filling fluid and the solid matrix. The microscopic distribution of the internal magnetic field, and its gradients, was simulated based on the thin-section pore structure of a sedimentary rock. The simulation results were verified experimentally. We show that the 'decay due to diffusion in internal field' magnetic resonance technique may be applied to measure the pore size distribution in partially saturated porous media. For the first time, we have observed that the internal magnetic field and its gradients in porous rocks have a Lorentzian distribution, with an average gradient value of zero. The Lorentzian distribution of internal magnetic field arises from the large susceptibility contrast and an intrinsic disordered pore structure in these porous media. We confirm that the single exponential magnetic resonance free induction decay commonly observed in fluid saturated porous media arises from a Lorentzian internal field distribution. A linear relationship between the magnetic resonance linewidth, and the product of the susceptibility difference in the porous media and the applied magnetic field, is observed through simulation and experiment.     

多孔介质中卡森流体的分形分析

研究了非牛顿流体中的卡森流体在多孔介质中的流动特性.基于服从分形分布的弯曲毛细管束模型,运用分形几何理论推导出了该流体在多孔介质中流动的流量、流速、启动压力梯度和有效渗透率的分形解析解.模型中的每一个参数都有明确的物理意义,它将卡森流体在多孔介质中的流动特性与多孔介质的微结构参数有机联系起来.文中给出了卡森流体的流速、启动压力梯度和有效渗透率随着各影响因素的变化趋势,并进行了讨论.所得分形模型可以更深刻地理解卡森流体在多孔介质中流动的内在物理机理. 关键词： 多孔介质 卡森流体 分形     

Direct measurement of porous media local hydrodynamical permeability using gas MRI.

The concept of hydraulic permeability is at the core of modeling single phase or multi-phase flow in heterogeneous porous media, as it is the spatial distribution of the permeability that primarily governs the behavior of fluid flow in the medium. To date, the modeling of fluid flow in porous media has been hampered by poor estimates of local permeability. Magnetic Resonance Imaging is well known for its ability to measure non-invasively the local density and flow rate of different fluids saturating porous media [1,2]. In this paper we demonstrate the first non-invasive method for the direct measurement of a single projection of the local permeability tensor of a porous medium using gas-phase MRI. The potential for three-dimensional imaging of the medium permeability is also discussed. The limitations of the method are listed and results are presented in a model porous medium as well as in a real oil reservoir rock.     

Correlation functions for inhomogeneous magnetic field in random media with application to a dense random pack of spheres

In porous media subject to applied magnetic field, the internal field arises out of susceptibility contrast of the constituents. We have examined the spatial inhomogeneity of the internal fields in a random pack of spheres using numerical computation. We find that the pair-correlation function of the internal field (K2) is a close approximation to the structure factor of the material, thus K2 can be used to characterize pore geometry. The magnetic length scale LambdaM exhibited in K2 is shown to be related to the fluid transport in the medium.     

非均一多孔介质中的水热迁移研究

孔隙裂隙非均一多孔介质中水热迁移的研究，国际上只是近年来才开始取得明显的进展。在发育裂隙的孔隙岩层中同时存在着两种渗流系统：孔隙总体积较大、渗透性相对弱的多孔岩块和总体积较小、渗透性却相对较强的分割多孔块体的裂隙。从而提出了“孔隙一裂隙二重性”假定，即地下水主要贮存在孔隙中，而水的运动主要在裂隙中进行，用一个一阶量描述孔隙裂隙间水流（热流）的传递耦合项。本文导出描述孔隙─裂隙岩层中的水流和热迁移的基本微分方程，建立起相应的数学模型，并成功地用于我国西藏羊八井热田分布参数模型的水热迁移研究。     

气孔分布特性对含水多孔介质气体扩散的影响

本文基于多孔介质的气孔分布特性,计算了多孔介质在含水状态下的扩散性能,并且比较了采用两种方式计算相对渗透率时的相对扩散性能。其结果表明,基于气孔分布的计算结果低于与气孔分布无关的计算结果。另外,疏水性含水多孔介质的扩散性能低于亲水性含水多孔介质的扩散性能,基于气孔分布计算含水多孔介质的气体扩散性能时,Wyllie公式并不适用。     

Pores resolving simulation of Darcy flows

A theoretical formulation and corresponding numerical solutions are presented for microscopic fluid flows in porous media with the domain sufficiently large to reproduce integral Darcy scale effects. Pore space geometry and topology influence flow through media, but the difficulty of observing the configurations of real pore spaces limits understanding of their effects. A rigorous direct numerical simulation (DNS) of percolating flows is a formidable task due to intricacies of internal boundaries of the pore space. Representing the grain size distribution by means of repelling body forces in the equations governing fluid motion greatly simplifies computational efforts. An accurate representation of pore-scale geometry requires that within the solid the repelling forces attenuate flow to stagnation in a short time compared to the characteristic time scale of the pore-scale flow. In the computational model this is achieved by adopting an implicit immersed-boundary method with the attenuation time scale smaller than the time step of an explicit fluid model. A series of numerical simulations of the flow through randomly generated media of different porosities show that computational experiments can be equivalent to physical experiments with the added advantage of nearly complete observability. Besides obtaining macroscopic measures of permeability and tortuosity, numerical experiments can shed light on the effect of the pore space structure on bulk properties of Darcy scale flows.     

XFEM modeling of hydraulic fracture in porous rocks with natural fractures

Hydraulic fracture (HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method (XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures (NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk (KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.     

Laboratory observations of altered porous fluid flow behavior in Berea sandstone induced by low-frequency dynamic stress stimulation

It has been observed repeatedly that low-frequency (1–500 Hz) seismic stress waves can enhance oil production from depleted reservoirs and contaminant extraction from groundwater aquifers. The physics coupling stress waves to fluid flow behavior in porous media is not understood, although numerous physical mechanisms have been proposed to explain the observations. To quantify the effects of low-frequency, dynamic-stress stimulation on multiphase fluid flow and in situ particle behavior in porous media, laboratory experiments were conducted with a core flow stimulation apparatus that allows for precise control and measurement of applied stress and strain, static confinement, and fluid flow parameters. Results are reported for experiments on stimulated single-phase and two-phase fluid flow behavior in 2.54-cm-diameter Berea sandstone cores. For all experiments, stimulation was applied to the cores in the form of sinusoidal, axial, mechanical stress coupled to the solid porous matrix at frequencies of 25 to 75 Hz. Applied stress RMS amplitudes ranged from 300 to 1200 kPa and, at these levels, produced coupled, pore-pressure fluctuations of much less than 1.2 to 4.8 kPa, respectively. During single-phase brine flow, stimulation increased the absolute permeability of the rock by 10–20%. This was caused by mobilizing in situ clay particles that were partially plugging the pore throats. During two-phase, steady-state, constant-rate flow of oil-brine and decane-brine mixtures, stimulation caused significant changes in the bulk fluid pressure drop across the core. The pressure changes showed a strong dependence on the viscosity of the nonwetting fluid phase (oil or decane) relative to the wetting phase (brine). This may indicate that relative changes in the mobility of wetting versus nonwetting fluid phases were induced by the dynamic stress. Under the specific experimental conditions used, pore-scale particle perturbation and altered wettability are possible physical mechanisms that can explain the results.     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

多孔介质内溶解与沉淀过程的格子Boltzmann方法模拟

本文采用格子Boltzmann方法模拟了多孔介质内的溶解和沉淀现象, 并分析了雷诺数、施密特数、达姆科勒数对多孔介质孔隙结构及浓度分布的影响. 结果表明: 对于多孔介质内的溶解(沉淀)过程, 当雷诺数越大时, 孔隙率越大(小), 平均浓度值越小(大); 当达姆科勒数或施密特数较小时, 溶解和沉淀过程均受反应控制, 此时反应在多孔介质的固体表面较为均匀的发生; 当达姆科勒数或施密特数较大时, 溶解和沉淀过程均受扩散控制, 此时反应主要发生在上游及大孔隙区域.     

Elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media

This paper studies the elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media with the nonzero boundary slip velocity for pore size distribution. The coefficient bF_m(ω) that measures the deviation from Poiseuille flow friction in such media is presented. Based on this coefficient, we investigate the properties of elastic waves by calculating their phase velocities and attenuation coefficients as functions of frequency and the behaviour of the dynamic permeability. The study shows that the pore size distribution removes oscillations in all physical quantities in the non-Newtonian regime. Consideration of the nonzero boundary slip effect in non-Newtonian (Maxwell) fluid-saturated porous media results in (a) an overall increase of the dynamic permeability, (b) an increase of phase velocities of fast Biot waves and shear waves except in the low frequency domain and an overall increase of phase velocity of slow Biot waves and (c) an overall increase of the attenuation of three Biot waves in the intermediate frequency domain except in the deeply non-Newtonian regime. The study also shows that the attenuation coefficient of slow Biot waves is small in the deeply non-Newtonian regime at higher frequency, which encourages us to detect slow Biot waves in oil-saturated porous rock.     

Sodium NMR relaxation in porous materials

The NMR relaxation of hydrogen nuclei of a fluid in a porous material is generally interpreted in terms of the Brownstein-Tarr model, in which the relaxation rate of the signal is inversely proportional to the pore size. We have investigated whether this model can be applied to the relaxation of Na nuclei in a NaCl solution in a porous material. The results indicate that the ion distribution over the pores can be obtained from an analysis of the Na NMR signal decay, if the pore sizes are roughly below 1 microm. This information is very useful for studies of combined moisture and ion transport in porous building materials.     

An analysis of the radial flow in the heterogeneous porous media based on fractal and constructal tree networks

In this paper, an analysis of the radial flow in the heterogeneous porous media based on fractal and constructal tree networks is presented. A dual-domain model is applied to simulate the heterogeneous porous media embedded with a constructal tree network based on the fractal distribution of pore space and tortuosity nature of flow paths. The analytical expressions for seepage velocity, pressure drop, local and global permeability of the network and binary system are derived, and the transport properties for the optimal branching structure are discussed. Notable is that the global permeability (K_n) of the network and the volume fraction (f_n) occupied by the network exhibit linear scaling law with the fractal dimension (D_p) of channel diameter bylogK_n∼0.46D_p and logf_n∼1.03D_p, respectively. Our analytical results are in good agreement with the available numerical results for steady-state soil vapor extraction and indicate that the fractal dimension for pore space has significant effect on the permeable properties of the media. The proposed dual-domain model may capture the characteristics of heterogeneous porous media and help understanding the transport mechanisms of the radial flow in the media.     

Effective Gradients in Porous Media Due to Susceptibility Differences

In porous media, magnetic susceptibility differences between the solid phase and the fluid filling the pore space lead to field inhomogeneities inside the pore space. In many cases, diffusion of the spins in the fluid phase through these internal inhomogeneities controls the transverse decay rate of the NMR signal. In disordered porous media such as sedimentary rocks, a detailed evaluation of this process is in practice not possible because the field inhomogeneities depend not only on the susceptibility difference but also on the details of the pore geometry. In this report, the major features of diffusion in internal gradients are analyzed with the concept of effective gradients. Effective gradients are related to the field inhomogeneities over the dephasing length, the typical length over which the spins diffuse before they dephase. For the CPMG sequence, the dependence of relaxation rate on echo spacing can be described to first order by a distribution of effective gradients. It is argued that for a given susceptibility difference, there is a maximum value for these effective gradients,g_max, that depends on only the diffusion coefficient, the Larmor frequency, and the susceptibility difference. This analysis is applied to the case of water-saturated sedimentary rocks. From a set of NMR measurements and a compilation of a large number of susceptibility measurements, we conclude that the effective gradients in carbonates are typically smaller than gradients of current NMR well logging tools, whereas in many sandstones, internal gradients can be comparable to or larger than tool gradients.