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

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

Pore-Size Distributions and Tortuosity in Heterogeneous Porous Media

The diffusion coefficient, measured at long observation times by pulsed-held-gradient NMR, can in principle be used to estimate the tortuosity of a porous medium. This method is useful for glass-sphere packs, but we find that it does not generally work for porous sedimentary rock. Natural sedimentary rocks are characterized by complex microgeometries and broad distributions of pore sizes, which cannot be adequately sampled by diffusing molecules in experimentally accessible observation times. The time-dependent diffusion coefficient D(t) can be distinctly irregular for rocks with very large pores. In heterogeneous porous media, determination of pore-size distribution by relaxation-time measurements and tortuosity by PFG diffusion measurements are mutually exclusive.     

A model for strong attenuation and dispersion of seismic P-waves in a partially saturated fractured reservoir

In this work we interpret the data showing unusually strong velocity dispersion of P-waves (up to 30%) and attenuation in a relatively narrow frequency range. The cross-hole and VSP data were measured in a reservoir, which is in the porous zone of the Silurian Kankakee Limestone Formation formed by vertical fractures within a porous matrix saturated by oil, and gas patches. Such a medium exhibits significant attenuation due to wave-induced fluid flow across the interfaces between different types of inclusions (fractures, fluid patches) and background. Other models of intrinsic attenuation (in particular squirt flow models) cannot explain the amount of observed dispersion when using realistic rock properties. In order to interpret data in a satisfactory way we develop a superposition model for fractured porous rocks accounting also for the patchy saturation effect.     

Stress-driven phase transformation and the roughening of solid-solid interfaces

The application of stress to multiphase solid-liquid systems often results in morphological instabilities. Here we propose a solid-solid phase transformation model for roughening instability in the interface between two porous materials with different porosities under normal compression stresses. This instability is triggered by a finite jump in the free energy density across the interface, and it leads to the formation of fingerlike structures aligned with the principal direction of compaction. The model is proposed as an explanation for the roughening of stylolites-irregular interfaces associated with the compaction of sedimentary rocks that fluctuate about a plane perpendicular to the principal direction of compaction.     

低本底伽马能谱测量在沉积岩分类判别中的应用及其影响因素探讨

沉积岩中放射性元素铀、钍、钾的含量主要取决于岩石中泥质的含量,而泥质含量的高低是判别沉积岩类别的主要依据,因此,可以将铀、钍、钾的比活度或含量值作为判别沉积岩岩类的量化指标。利用低本底伽马能谱法测量岩石样品中铀、钍、钾的比活度或含量,可以快速、准确地进行岩类判别。由于放射性测量受样品几何形状、样品质量及其所含水分含量变化的影响,因此岩类判别的准确性会受到影响。通过理论探讨和实验验证,发现各待测核素特征峰计数率值与样品质量成线性关系, 且不同的能量区域、不同的岩性,线性关系系数与趋势拟合程度都是不同的;当样品的水分含量不超过10%(取样时,岩屑水分不超过10%)时,对测量结果的影响很小(变化值均在一倍均方差以内),因此,水分含量对低本底伽马能谱法判别沉积岩岩类的准确性不产生显著的影响,一般可以不予考虑。对陕西定边某油气勘查区钻孔岩屑样品进行判别试验,对实测数据进行质量校正后,利用铀、钍、钾道的“归一化”计数率值(即单位质量的计数率)只能粗略地划分不同沉积岩的大类,再以铀、钍、钾特征峰的合峰计数率建立综合判别模式,能够对岩类进行进一步细分,且判别准确度达到75%以上。     

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.     

Grain segregation mechanism in aeolian sand ripples

Many sedimentary rocks are formed by migration of sand ripples. Thin layers of coarse and fine sand are present in these rocks, and understanding how layers in sandstone are created has been a longstanding question. Here, we propose a mechanism for the origin of the most common layered sedimentary structures such as inverse graded climbing ripple lamination and cross-stratification patterns. The mechanism involves a competition between three segregation processes: (i) size-segregation and (ii) shape-segregation during transport and rolling, and (iii) size segregation due to different hopping lengths of the small and large grains. We develop a discrete model of grain dynamics which incorporates the coupling between moving grains and the static sand surface, as well as the different properties of grains, such as size and roughness, in order to test the plausibility of this physical mechanism. Received 19 July 1999 and Received in final form 4 August 1999     

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.     

The use of apatite fission track thermochronology to constrain fault movements and sedimentary basin evolution in northeastern Brazil

An apatite fission track study of crystalline rocks underlying sedimentary basins in northeastern Brazil indicate that crustal blocks that occur on opposite sides of a geological fault experienced different thermal histories. Samples collected on the West block yielded corrected fission-track ages from 140 to 375 Ma, whereas samples collected on the East block yielded ages between 90 and 125 Ma. The thermal models suggest that each block experienced two cooling events separated by a heating event at different times. We concluded that the West block moved downward relative to the East block ca. 140 Ma ago, when sediments eroded from the East side were deposited on the West side. This process represents the early stage of sedimentary basin formation and the opening of the South Atlantic Ocean in the region. Downward and upward movements related to heating and cooling events of these crustal blocks at different periods until recent times are proposed.     

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.     

Nd-isotopic composition of Phanerozoic sediments in the Inner Zone of Southwest Japan Arc: implications on provenance characteristics and contribution to formation of mature island arc system

The Nd-isotopic data on sedimentary and metamorphic rocks of SW Japan Arc allow their discrimination into five different depleted mantle (Tdm) model age clusters, 2.6-2.45 Ga, 2.3-2.05 Ga, 1.9-1.55 Ga, 1.45-1.25 Ga, 1.2-0.85 Ga. The 2.6-2.45 Ga and 1.9-1.55 Ga model ages are also coincident with U-Pb inherited zircon ages of the above two epochs as well as the major magmatic activity in the Sino-Korean Craton (SKC). The 2.3-2.05 Ga model ages can be considered as the initial formation ages for the precursors of sedimentary rocks. The Nd-isotopic data suggest that the Hida Belt was most likely formed as a part of the SKC. The mantle underlying the Ryoke Belt had continental lithospheric signature during Triassic-Jurassic period. The 1.9-1.55 Ga model ages, especially 1.8 Ga~, can be associated with the formation of this belt. The source material for the sedimentary rocks occurring in the accretionary terrane of northeastern areas in the SW Japan Arc was probably in and around the SKC of the Ryoke Belt itself. The sedimentary rocks occurring in southwestern areas of the Arc were mainly composed of materials derived from a relatively younger source (1.45-0.85 Ga).     

On the material constants measurement method of a fluid-saturated transversely isotropic poroelastic medium

The fluid-saturated transversely isotropic poroelastic medium could be widely found in nature, e.g., the sedimentary rocks underground. To determine the eight independent material constants for the transversely isotropic poroelastic medium, a series of tests are discussed. Two undrained tests and one drained test are suggested as a set of tests of the least amount. For the verification purpose, two additional drained tests are also introduced as an option. The atmospheric dried test is discussed as a replacement of the traditional infiltrated drained test to save the time waiting for an equilibrium state. Some microscopic material constants, i.e.,the unjacketed bulk coefficients, the porosity, and the compressibility of porous fluid, are measurable but unnecessary to determine the independent material constants of a poroelastic medium.     

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.     

Probing the analogy of the proton relaxation in biological tissues and porous media

Summary Experimental data are reported that show the analogy of longitudinal and transverse proton relaxation in heterogeneous systems as different as biological tissues and water-saturated rocks. Published data on the τ-dependence of the transverse-relaxation rate for biological tissues with magnetite grains, used as a contrast agent in MRI, are discussed in the light of our recent results on water-saturated porous media, and give for the liver another case of a behaviour parallel to that in rocks. There are enough similarities between NMR relaxation in tissues and in other porous media that, for work in either area, attention to the other is likely to be fruitful. Work supported by Italian CNR and MURST grants     

Flow through a porous membrane simulated by cellular automata and by finite elements

Computational results concerning incompressible viscous flow through two channels connected by a porous membrane are presented. The example is extraordinary for its four different types of boundary conditions that are necessary to make the problem complete. The solution is accomplished by two methods: by cellular automata and by finite elements. The numerical means to satisfy the boundary conditions are given for both methods. Overall agreement is achieved, but significant differences show up in details.     

The effect of external mean flow on sound transmission through double-walled cylindrical shells lined with poroelastic material

Sound transmission through a system of double shells, lined with poroelastic material in the presence of external mean flow, is studied. The porous material is modeled as an equivalent fluid because shear wave contributions are known to be insignificant. This is achieved by accounting for the energetically most dominant wave types in the calculations. The transmission characteristics of the sandwich construction are presented for different incidence angles and Mach numbers over a wide frequency range. It is noted that the transmission loss exhibits three dips on the frequency axis as opposed to flat panels where there are only two such frequencies—results are discussed in the light of these observations. Flow is shown to decrease the transmission loss below the ring frequency, but increase this above the ring frequency due to the negative stiffness and the damping effect added by the flow. In the absence of external mean flow, porous material provides superior insulation for most part of the frequency band of interest. However, in the presence of external flow, this is true only below the ring frequency—above this frequency, the presence of air gap in sandwich constructions is the dominant factor that determines the acoustic performance. In the absence of external flow, an air gap always improves sound insulation.     

兴城地区岩石光谱特征与物性参数、金属元素间响应关系研究

岩石物性参数、元素含量和光谱特征三者之间的相互关系并非独立存在, 研究三者之间的相互关系,对于探索通过遥感信息进行岩石矿物成分、物性参数的定量反演方法奠定基础。采集了兴城地区590块岩样本,分析其物性参数（密度、介电常数、电阻率、磁化率）、金属成分含量（Fe, Ti, V, Mn, Zr, Co, Zn, Nb, Bi, Pb）以及岩石光谱之间的关系。将各物化参数与原始光谱、光谱吸收深度、光谱小波包分解后的高低频做相关性分析,找出各物化参数影响光谱吸收和反射的特征波段,探究关系密切的参数。该研究为岩石岩性分类、某种金属元素和物理特性的反演以及用某种参数来预测关系密切的参数奠定了基础,取得如下结论：（1）获得了火成岩中Fe,Ti,Mn,V,Zn,Bi和Pb等元素的特征谱带,其中Fe元素含量更高,与光谱的相关性更显著。在0.4和0.54 μm波段附近存在Fe的特征反射波段,1.0~1.5 μm波段范围内存在Fe元素的特征吸收波段;在0.4~0.55和0.6~0.65 μm波段范围内,Ti元素与光谱的相关性更显著,在2.28 μm波段附近存在Ti的特征吸收波段;在0.41 μm波段附近存在Mn的特征反射波段;火成岩、沉积岩V元素与光谱相关性差异性较大,在0.76,0.81,0.89,0.95 μm波段附近可能存在火成岩的特征吸收波段和沉积岩的特征反射波段;沉积岩的Zn元素含量与光谱的相关性比火成岩要显著,在0.41,1.36,1.59 μm波段附近可能存在火成岩Zn的特征反射波段,在2.34 μm波段附近可能存在沉积岩Zn的特征波段;2.14 μm波段附近,Bi元素对沉积岩光谱的吸收有一定的作用;Pb的特征谱带可能存在于0.45,0.54,2.29 μm附近;（2）在岩石各物理特性和光谱的关系研究中,在0.57~0.85 μm波段范围内,密度对光谱有很好的吸收反射特征,在0.53 μm波段附近,磁化率使光谱有较强的反射,在1.08 μm波段附近,磁化率使光谱产生较强的吸收;电阻率和光谱的相关关系与密度和光谱的相关关系极为相似;（3）在岩石各物性参数间的关系中发现,密度和电阻率呈显著正相关的关系;（4）在岩石各物性参数与金属元素的关系中研究发现,密度与各金属元素相关性较弱;磁化率与Fe和Ti元素呈显著的正相关关系;电阻率与各金属元素间的相关关系较弱;介电常数与各金属元素含量呈正相关关系,其中与V,Zn,Bi元素的相关关系最显著;（5）在金属元素间,Fe与Ti有较显著的正相关关系,Ti与Fe,V元素间有较强的正相关关系,Zn和Pb存在较强的正相关。     

Magnetic resonance imaging: applications of novel methods in studies of porous media.

NMR imaging is finding broad applications in nonbiological areas including the study of fluid flow and fluid ingress in porous media. The porous media include at the one end mineral rocks and various building materials through various solid plastic materials to foodstuffs at the other end of the spectrum. The fluids within these various media range from crude oil and water mixtures, and water itself, to a range of organic solvents in plastic materials. This paper is concerned with the flow and ingress of water through Bentheimer sandstone and Ninian reservoir specimens, and also in solid nylon blocks.     

Computation of dynamic seismic responses to viscous fluid of digitized three-dimensional Berea sandstones with a coupled finite-difference method

The seismic response of saturated porous rocks is studied numerically using microtomographic images of three-dimensional digitized Berea sandstones. A stress-strain calculation is employed to compute the velocities and attenuations of rock samples whose sizes are much smaller than the seismic wavelength of interest. To compensate for the contributions of small cracks lost in the imaging process to the total velocity and attenuation, a hybrid method is developed to recover the crack distribution, in which the differential effective medium theory, the Kuster-Tokso?z model, and a modified squirt-flow model are utilized in a two-step Monte Carlo inversion. In the inversion, the velocities of P- and S-waves measured for the dry and water-saturated cases, and the measured attenuation of P-waves for different fluids are used. By using such a hybrid method, both the velocities of saturated porous rocks and the attenuations are predicted accurately when compared to laboratory data. The hybrid method is a practical way to model numerically the seismic properties of saturated porous rocks until very high resolution digital data are available. Cracks lost in the imaging process are critical for accurately predicting velocities and attenuations of saturated porous rocks.