3D image-based characterization of fluid displacement in a Berea core

Improved network flow models require the incorporation of increasingly accurate geometrical characterization of the microscale pore structure as well as greater information on fluid–fluid interaction (interfaces) at pore scales. We report on three dimensional (3D) pore scale medium characterization, absolute permeability computations for throat structures, and pore scale residual fluid distribution in a Berea core. X-ray computed microtomography combined with X-ray attenuating dopants is used to obtain 3D images of the pore network and to resolve phase distributions in the pore space.     

Approximation of interfacial properties in multiphase porous medium systems

We investigate the estimation of interfacial areas, curvatures, and common curve lengths in multiphase porous medium systems. Algorithms are developed to obtain estimates of these quantities based upon a variety of potential data sources and estimation approaches. The accuracy of the derived approximations are evaluated as a function of the data type and resolution of the data. The methods advanced improve upon standard approaches now in use and show excellent accuracy at resolutions on the order of five lattice points per minimum radius of curvature of the object being resolved. Finally, we suggest a promising class of extensions that could lead to further improvements in the accuracy of such methods.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

Quantification of salt weathering in porous stones using an experimental continuous partial immersion method

In this study, an experimental salt weathering simulation and porous stone durability classification are proposed. There are many laboratory tests that quantify durability against salt crystallisation weathering action. These are usually based on the total immersion of samples into a saline solution, which is not representative of the salt weathering mechanism. An experimental test based on partial immersion is suggested. This is a comparable study of weight loss and degradation of visual appearance due to salt crystallisation using, on the one hand, a standard durability test (UNE), and, on the other, the proposed durability test. The weight loss and visual appearance in our test is comparable to the degradation of building stone. The differences between weight loss data in both tests depend on the petrophysical properties: porous media and degree of coherence.

From this testing, a new durability classification as a function of dry weight loss in the partial immersion test is proposed. Four divisions of different types of materials can be made in this classification, which quantifies salt weathering action mainly in environments and mild climatic conditions.     

声波在两种多孔介质界面上的反射和透射

本文导出了声波在两种多孔介质界面上反射、透射的一般计算公式.作为例子,数值计算了P₁波入射于界面时,P₁、P₂和S波的反射、透射系数与声波频率、入射角等量之间的关系.结果表明,各种模式波的反射、透射系数与入射角、多孔介质性质有关,在Biot特征频率附近与频率有关,并用界面两侧的法向能流相等验证了结果的正确性.若把多孔介质当作均匀固体处理,将会得到显著不同的结果.     

Another look at the conceptual fundamentals of porous media upscaling

We suggest a critical look at the epistemic foundations of the porous media upscaling problem that focuses on conceptual processes at work and not merely on form manipulations. We explore the way in which critical aspects of scientific methodology make their appearance in the upscaling context, thus generating useful effective parameters in practice. The fons et origo of our approach is a conceptual blending of knowledge states that requires the revision of the traditional method of scientific argument underlying most upscaling techniques. By contrast to previous techniques, the scientific reasoning of the proposed upscaling approach is based on a stochastic model that involves teleologic solutions and stochastic logic integration principles. The syllogistic form of the approach has important advantages over the traditional reasoning scheme of porous media upscaling, such as: it allows the rigorous derivation of the joint probability distributions of hydraulic gradients and conductivities across space; it imposes no restriction on the functional form of the effective parameters or the shape of the probability laws governing the random media (non-Gaussian distributions, multiple-point statistics and non-linear models are automatically incorporated); it relies on sound methodological principles rather than being ad hoc; and it offers the rational means for integrating the multifarious core knowledge bases and uncertain site-specific information sources about the subsurface system. Previous upscaling results are derived as special cases of the proposed upscaling approach under limited conditions of porous media flow, a fact that further demonstrates the generalization power of the approach. Our hope is that looking at the upscaling problem in this novel way will direct further attention to the methodological exploration of the problem at the length and the detail that it deserves.I would like to thank Drs. A. Kolovos and D.T. Hristopulos for their valuable comments. The work was supported by grants from the Army Research Office (Grant no. DAAG55–98–1-0289) and the National Institute of Environmental Health Sciences (P42-ES05948 & P30-ES10126).     

Nutrients adsorption from seawater by new porous carrier made from zeolitized fly ash and slag

This study investigated the adsorption and precipitation of phosphate by blast furnace slag (BFS) separately.

In order to evaluate the adsorption capacity of BFS, BFS was treated before its use by acid. The authors aim to develop a new porous carrier to adsorb simultaneously ammonium and phosphate from seawater under eutrophic conditions. The current paper deals with a promising new approach to improve the utilization of some industrial solid wastes such as BFS and zeolite synthesized from fly ash [ZFA(Fe)] by their solidification to cylindrical porous carriers using a hydrothermal hot-pressing (HHP) method.

Attempts to produce porous carriers using an arranged HHP method with different porosities (24%, 40% and 52% (v/v)) were carried out. Physical properties of carriers such as porosity, compressive strength and height have been investigated. Laboratory studies showed strong evidence that the porous carrier was very selective towards phosphate and ammonium. The results demonstrated the role of porosity in enhancing phosphate and ammonium adsorption by the increase of the surface area per weight. The estimates of the parameters and the correlation coefficients according to the Freundlich equations revealed that adsorption was related to the porosity of carriers and phosphate and ammonium were adsorbed well on the carriers having large porosity.

The results suggested that developing carrier with high porosity was a promising way to enhance nutrients adsorption.     

Convective stability analysis of the long-term storage of carbon dioxide in deep saline aquifers

Deep saline aquifers are one of the most suitable geologic formations for carbon sequestration. The linear and global stability analysis of the time-dependent density-driven convection in deep saline aquifers is presented for long-term storage of carbon dioxide (CO₂). The convective mixing that can greatly accelerate the CO₂ dissolution into saline aquifers arises because the density of brine increases upon the dissolution of CO₂ and such a density difference may induce instability. The effects of anisotropic permeability on the stability criteria, such as the critical time for the appearance of convective phenomena and the critical wavelength of the most unstable perturbation, are investigated with linear and global stability analysis. The linear stability analysis provides a sufficient condition for instability while the global stability analysis yields a sufficient condition for stability. The results obtained from these two approaches are not exactly the same but show a consistent trend, both indicating that the anisotropic system becomes more unstable when either the vertical or horizontal permeability increases.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Rock-magnetic properties of Neogene sediments, northern flank of Tianshan Mountains

Analyses of rock-magnetic properties of Neogene sediments of the Taxihe section, northern Tianshan Mountains, show that the section can be classified into three categories including lacustrine facies, fluvial facies and alluvial facies, which correspond to the lower, middle and upper of the Taxihe section respectively. The magnetic minerals of the lacustrine facies may be affected by the process of weath- ering, lithogenesis and biolithogenesis besides the source of the sediments. The natural remanence intensities are between 10-3 A/m and 10-2 A/m. The minerals are dominated by magnetite and the high coercive magnetic mineral may be goethite. The magnetic grains are the mixture of PSD SD or SD SP. The natural remanence intensities of the strata of fluvial facies are between 10-2 A/m and 10-1 A/m, about ten times that of the lacustrine facies. The magnetic minerals are mainly magnetite and hematite, and the magnetic grains are mainly PSD. The characteristic remanence (ChRM) carriers are magnetites. In the alluvial facies, the natural remanence intensities are mostly less than 1×10-2 A/m. The magnetic minerals of the series are dominated by magnetite and hematite, almost the same as the fluvial facies. But the difference is that most of the stepwise demagnetization can reveal two components and the ChRM carriers are hematites. The magnetic grains are PSD in terms of the hysteresis parameters.