Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Pore fluid pressure at DSDP Leg 84 Holes

Evidence of considerable overpressuring of pore fluids in the sediment drilled during Leg 84 was obtained from direct measurement of pressure by two methods. The first involved measurement of back pressure when the annulus of the drill hole became constricted with unremoved drill cuttings or constriction was caused by plastic inflow of the drill hole walls. The second involved measurement of pressure ahead of the bit in conjunction with in situ water samples and heat flow. All measurements indicated abnormally high pore pressure even in slope deposits of the Middle America Trench off Guatemala.

Sediment and pore fluid chemistry for IODP Sites 311-U1325 and 311-U1329

Analytical challenges in obtaining high quality measurements of rare earth elements (REEs) from small pore fluid volumes have limited the application of REEs as deep fluid geochemical tracers. Using a recently developed analytical technique, we analyzed REEs from pore fluids collected from Sites U1325 and U1329, drilled on the northern Cascadia margin during the Integrated Ocean Drilling Program (IODP) Expedition 311, to investigate the REE behavior during diagenesis and their utility as tracers of deep fluid migration. These sites were selected because they represent contrasting settings on an accretionary margin: a ponded basin at the toe of the margin, and the landward Tofino Basin near the shelf's edge. REE concentrations of pore fluid in the methanogenic zone at Sites U1325 and U1329 correlate positively with concentrations of dissolved organic carbon (DOC) and alkalinity. Fractionations across the REE series are driven by preferential complexation of the heavy REEs. Simultaneous enrichment of diagenetic indicators (DOC and alkalinity) and of REEs (in particular the heavy elements Ho to Lu), suggests that the heavy REEs are released during particulate organic carbon (POC) degradation and are subsequently chelated by DOC. REE concentrations are greater at Site U1325, a site where shorter residence times of POC in sulfate-bearing redox zones may enhance REE burial efficiency within sulfidic and methanogenic sediment zones where REE release ensues. Cross-plots of La concentrations versus Cl, Li and Sr delineate a distinct field for the deep fluids (z > 75 mbsf) at Site U1329, and indicate the presence of a fluid not observed at the other sites drilled on the Cascadia margin. Changes in REE patterns, the presence of a positive Eu anomaly, and other available geochemical data for this site suggest a complex hydrology and possible interaction with the igneous Crescent Terrane, located east of the drilled transect.

Ba, Al and SO4 concentrations in sediments and pore fluid of ODP Leg 127/128 samples

The barium distribution in sediments and pore fluids from five sites drilled in the Japan Sea have been used to illustrate the geochemical behavior of this element as it pertains paleoproductivity reconstructions, diagenetic remobilization, and barite precipitation in authigenic fronts. Sites where sulfate is depleted in the pore fluids also show high concentrations of dissolved barium, reflecting dissolution of biogenic barite. The high rate of sedimentation at Sites 798 and 799 results in a rapid sulfate depletion, which in turn leads to barite dissolution and reprecipitation in diagenetic fronts. The dissolved barium distribution at these sites has been used to quantify the rate of barite dissolution; we estimate a first-order rate constant for barite dissolution to be 2*10**-6/s at Site 799 and 2*10**-7/s at Site 798. Authigenic barite has been documented in sediments from Site 799 at 323 meters below seafloor by scanning electron microscopy and X-ray fluorescence analysis. These results indicate barite precipitation in a diagenetic front near the zone of sulfate depletion by upward migration of dissolved barium and downward diffusion of sulfate. Barite precipitation has also been inferred at Sites 796 and 798 based on sedimentary and dissolved barium distributions. Sulfate is not depleted in the pore fluids of Site 794. The lack of diagenetic remobilization of biogenic barium at this site preserves the high barium signal associated with the high-productivity sequences deposited during the late Miocene to Pliocene. Significantly, the organic carbon distribution does not indicate high accumulation rates during the periods of high opal and barium deposition. Instead, higher organic carbon accumulations are recorded in the Quaternary and middle Miocene sequences; intervals that are also characterized by deposition of siliciclastic turbidites. The presence of a terrestrial component in the organic carbon record renders barium a more useful indicator than organic carbon for paleoproductivity reconstructions in this marginal sea.