scholarly journals On the Consistency of Scale Among Experiments, Theory, and Simulation

2016 ◽  
Author(s):  
J. McClure ◽  
A. Dye ◽  
C. Miller ◽  
W. Gray

Abstract. The career of Professor Eric F. Wood has focused on the resolution of problems of scale in hydrologic systems. Within this context, we consider an evolving approach known as the thermodynamically constrained averaging theory (TCAT), which has broad applicability to hydrology. Specifically, we consider the case of modeling of two-fluid-phase flow in porous media. Two-fluid flow processes in the subsurface are fundamentally important for a wide range of hydrologic processes, including the transport of water and air in the vadose zone and geological carbon sequestration. Mathematical models that describe these complex processes have long relied on empirical approaches that neglect important aspects of the system behavior. New data sources make it possible to access the true geometry of geologic materials and directly measure previously inaccessible quantities. This information can be exploited to support a new generation of theoretical models that are constructed based on rigorous multiscale principles for thermodynamics and continuum mechanics. The challenges to constructing a mature model are shown to involve issues of scale, consistency requirements, appropriate representation of operative physical mechanisms at the target scale of the model, and a robust structure to support model evaluation, validation, and refinement. We apply TCAT to perform physics-based data assimilation to understand how the internal behavior influences the macroscale state of two-fluid porous medium systems. Examples of a microfluidic experimental method and a lattice Boltzmann simulation method are used to examine a key deficiency associated with standard approaches. In a hydrologic process such as evaporation, the water content will ultimately be reduced below the irreducible wetting phase saturation determined from experiments. This is problematic since the derived closure relationships cannot predict the associated capillary pressures for these states. In this work, we demonstrate that the irreducible wetting-phase saturation is an artifact of the experimental design, caused by the fact that the boundary pressure difference does not approximate the true capillary pressure. Using averaging methods, we measure the true capillary pressure for fluid configurations at and below the irreducible wetting phase saturation. Results of our analysis include a state function for the capillary pressure expressed as a function of fluid saturation and interfacial area.

2017 ◽  
Vol 21 (2) ◽  
pp. 1063-1076 ◽  
Author(s):  
James E. McClure ◽  
Amanda L. Dye ◽  
Cass T. Miller ◽  
William G. Gray

Abstract. As a tool for addressing problems of scale, we consider an evolving approach known as the thermodynamically constrained averaging theory (TCAT), which has broad applicability to hydrology. We consider the case of modeling of two-fluid-phase flow in porous media, and we focus on issues of scale as they relate to various measures of pressure, capillary pressure, and state equations needed to produce solvable models. We apply TCAT to perform physics-based data assimilation to understand how the internal behavior influences the macroscale state of two-fluid porous medium systems. A microfluidic experimental method and a lattice Boltzmann simulation method are used to examine a key deficiency associated with standard approaches. In a hydrologic process such as evaporation, the water content will ultimately be reduced below the irreducible wetting-phase saturation determined from experiments. This is problematic since the derived closure relationships cannot predict the associated capillary pressures for these states. We demonstrate that the irreducible wetting-phase saturation is an artifact of the experimental design, caused by the fact that the boundary pressure difference does not approximate the true capillary pressure. Using averaging methods, we compute the true capillary pressure for fluid configurations at and below the irreducible wetting-phase saturation. Results of our analysis include a state function for the capillary pressure expressed as a function of fluid saturation and interfacial area.


The quality of a reservoir can be described in details by the application of seismo electric transfer function fractal dimension. The objective of this research is to calculate fractal dimension from the relationship among seismo electric transfer fuction, maximum seismo electric transfer function and wetting phase saturation and to confirm it by the fractal dimension derived from the relationship among capillary pressure and wetting phase saturation. In this research, porosity was measured on real collected sandstone samples and permeability was calculated theoretically from capillary pressure profile measured by mercury intrusion techniques. Two equations for calculating the fractal dimensions have been employed. The first one describes the functional relationship between wetting phase saturation, seismo electric transfer function, maximum seismo electric transfer function and fractal dimension. The second equation implies to the wetting phase saturation as a function of capillary pressure and the fractal dimension. Two procedures for obtaining the fractal dimension have been developed. The first procedure was done by plotting the logarithm of the ratio between seismo electric transfer function and maximum seismo electric transfer function versus logarithm wetting phase saturation. The slope of the first procedure = 3- Df (fractal dimension). The second procedure for obtaining the fractal dimension was completed by plotting the logarithm of capillary pressure versus the logarithm of wetting phase saturation. The slope of the second procedure = Df -3. On the basis of the obtained results of the constructed stratigraphic column and the acquired values of the fractal dimension, the sandstones of the Shajara reservoirs of the Shajara Formation were divided here into three units. The gained units from bottom to top are: Lower Shajara Seismo Electric Transfer Function Fractal Dimension Unit, Middle Shajara Seismo Electric Tranfser Function Fractal dimension Unit, and Upper Shajara Seismo Electric Transfer Function Fractal Dimension Unit. The results show similarity between seismo electric transfer tunction fractal dimension and capillary pressure fractal dimension. It was also noted that samples with wide range of pore radius were characterized by high values of fractal dimension due to an increase in their connectivity and seismo electric transfer function. In our case , and as conclusions the higher the fractal dimension, the higher the permeability, the better the shajara reservoir characteristics.


2021 ◽  
pp. 1-20
Author(s):  
A. A. Kasha ◽  
A. Sakhaee-Pour ◽  
I. A. Hussein

Summary Capillary pressure plays an essential role in controlling multiphase flow in porous media and is often difficult to be estimated at subsurface conditions. The Leverett capillary pressure function J provides a convenient tool to address this shortcoming; however, its performance remains poor where there is a large scatter in the scaled data. Our aim, therefore, was to reduce the gaps between J curves and to develop a method that allows accurate scaling of capillary pressure. We developed two mathematical expressions based on permeability and porosity values of 214 rock samples taken from North America and the Middle East. Using the values as grouping features, we used pattern-recognition algorithms in machine learning to cluster the original data into different groups. In each wetting phase saturation, we were able to quantify the gaps between the J curves by determining the ratio of the maximum J to the minimum J. Graphical maps were developed to identify the corresponding group for a new rock sample after which the capillary pressure is estimated using the average J curve of the identified group and the permeability and porosity values of the rock sample. This method also provides better performance than the flow zone indicator (FZI) approach. The proposed technique was validated on six rock types and has successfully generated average capillary pressure curves that capture the trends and values of the experimentally measured data by mercury injection. Moreover, the proposed methodology in this study provides an advanced and a machine-learning-oriented approach for rock typing. In this paper, we provide a reliable and easy-to-use method for capillary pressure estimation in the absence of experimentally measured data by mercury injection.


2014 ◽  
Vol 2014 ◽  
pp. 1-12 ◽  
Author(s):  
Olugbenga Falode ◽  
Edo Manuel

An understanding of the mechanisms by which oil is displaced from porous media requires the knowledge of the role of wettability and capillary forces in the displacement process. The determination of representative capillary pressure (Pc) data and wettability index of a reservoir rock is needed for the prediction of the fluids distribution in the reservoir: the initial water saturation and the volume of reserves. This study shows how wettability alteration of an initially water-wet reservoir rock to oil-wet affects the properties that govern multiphase flow in porous media, that is, capillary pressure, relative permeability, and irreducible saturation. Initial water-wet reservoir core samples with porosities ranging from 23 to 33%, absolute air permeability of 50 to 233 md, and initial brine saturation of 63 to 87% were first tested as water-wet samples under air-brine system. This yielded irreducible wetting phase saturation of 19 to 21%. The samples were later tested after modifying their wettability to oil-wet using a surfactant obtained from glycerophtalic paint; and the results yielded irreducible wetting phase saturation of 25 to 34%. From the results of these experiments, changing the wettability of the samples to oil-wet improved the recovery of the wetting phase.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. M35-M53 ◽  
Author(s):  
Bastien Dupuy ◽  
Stéphane Garambois ◽  
Jean Virieux

The quantitative estimation of rock physics properties is of great importance in any reservoir characterization. We have studied the sensitivity of such poroelastic rock physics properties to various seismic viscoelastic attributes (velocities, quality factors, and density). Because we considered a generalized dynamic poroelastic model, our analysis was applicable to most kinds of rocks over a wide range of frequencies. The viscoelastic attributes computed by poroelastic forward modeling were used as input to a semiglobal optimization inversion code to estimate poroelastic properties (porosity, solid frame moduli, fluid phase properties, and saturation). The sensitivity studies that we used showed that it was best to consider an inversion system with enough input data to obtain accurate estimates. However, simultaneous inversion for the whole set of poroelastic parameters was problematic due to the large number of parameters and their trade-off. Consequently, we restricted the sensitivity tests to the estimation of specific poroelastic parameters by making appropriate assumptions on the fluid content and/or solid phases. Realistic a priori assumptions were made by using well data or regional geology knowledge. We found that (1) the estimation of frame properties was accurate as long as sufficient input data were available, (2) the estimation of permeability or fluid saturation depended strongly on the use of attenuation data, and (3) the fluid bulk modulus can be accurately inverted, whereas other fluid properties have a low sensitivity. Introducing errors in a priori rock physics properties linearly shifted the estimations, but not dramatically. Finally, an uncertainty analysis on seismic input data determined that, even if the inversion was reliable, the addition of more input data may be required to obtain accurate estimations if input data were erroneous.


2018 ◽  
Vol 3 (8) ◽  
Author(s):  
James E. McClure ◽  
Ryan T. Armstrong ◽  
Mark A. Berrill ◽  
Steffen Schlüter ◽  
Steffen Berg ◽  
...  

2015 ◽  
Vol 51 (7) ◽  
pp. 5365-5381 ◽  
Author(s):  
W. G. Gray ◽  
A. L. Dye ◽  
J. E. McClure ◽  
L. J. Pyrak‐Nolte ◽  
C. T. Miller

2016 ◽  
Vol 796 ◽  
pp. 211-232 ◽  
Author(s):  
J. E. McClure ◽  
M. A. Berrill ◽  
W. G. Gray ◽  
C. T. Miller

The movements of fluid–fluid interfaces and the common curve are an important aspect of two-fluid-phase flow through porous media. The focus of this work is to develop, apply and evaluate methods to simulate two-fluid-phase flow in porous medium systems at the microscale and to demonstrate how these results can be used to support evolving macroscale models. Of particular concern is the problem of spurious velocities that confound the accurate representation of interfacial dynamics in such systems. To circumvent this problem, a combined level-set and lattice-Boltzmann method is advanced to simulate and track the dynamics of the fluid–fluid interface and of the common curve during simulations of two-fluid-phase flow in porous media. We demonstrate that the interface and common curve velocities can be determined accurately, even when spurious currents are generated in the vicinity of interfaces. Static and dynamic contact angles are computed and shown to agree with existing slip models. A resolution study is presented for dynamic drainage and imbibition in a sphere pack, demonstrating the sensitivity of averaged quantities to resolution.


2010 ◽  
Vol 13 (03) ◽  
pp. 465-472 ◽  
Author(s):  
Amund Brautaset ◽  
Geir Ersland ◽  
Arne Graue

Summary During waterfloods of six outcrop chalk core-plug samples prepared at various wettabilities, simultaneous local pressures and in-situ fluid saturations were measured. Using high-spatial-resolution magnetic-resonance imaging (MRI) to image fluid saturations and pressure taps with semipermeable disks to measure individual phase pressures allowed calculations of relative permeabilities and the dynamic capillary pressure curves for the imbibition processes. A second objective was to identify individual-fluid saturation changes caused by spontaneous imbibition and viscous displacement to determine the local recovery mechanism and to calculate local recovery factors and in-situ Amott-Harvey indices. The obtained results contribute to improved description and understanding of multiphase-fluid flow in porous media, including in situ measurements of relative permeabilities, dynamic capillary pressure curves, Amott-Harvey Indices, and local oil-recovery mechanisms.


Sign in / Sign up

Export Citation Format

Share Document