Dispersion of a filtration stream in media with random inhomogeneities

The transport of a dynamically neutral impurity by a stream into a porous medium with random inhomogeneities is considered. In contrast to [1, 2], in which a Markov random-walk process of the impurity particles was postulated substantially (taking a hypothesis about Markov random walk processes contradicts to a definite degree the representation of particle motion along a streamline, finiteness of the velocity, and smoothness of the trajectory), the complete system of equations for the filtration concentration and velocities is investigated here by a perturbation method, which results in a non-local equation for the mean concentrations after taking the average. It is shown that the local equation (parabolic or hyperbolic) is the limit case in the scheme considered. The effect of a regular drift of saturation, analogous to the effect of directed transport in the theory of inhomogeneous turbulent diffusion [3], is established. One-dimensional, plane, and three-dimensional flows are considered. The fundamental relationships contain moments of the random velocity field. The relationship between these moments and the characteristics of the random permeability and porosity fields has been established in [1, 2].     

Conditional averaging of the equations of flow in random composite porous media

The problem of conditional averaging of the transport equations is solved for a neutral impurity in a composite medium with random porosity and impurity diffusion tensor. An unclosed system of conditionally averaged equations is constructed and closed using the globally averaged equations. The average impurity concentration fields for the individual phases of the composite medium and the phase-continuum interaction characteristics are calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 1, pp. 75–81, January–February, 1987.     

Averaging of stochastic evolution equations of transport in porous media

A study is made of the problem of averaging the simplest one-dimensional evolution equations of stochastic transport in a porous medium. A number of exact functional equations corresponding to distributions of the random parameters of a special form is obtained. In some cases, the functional equations can be localized and reduced to differential equations of fairly high order. The first part of the paper (Secs. 1–6) considers the process of transport of a neutral admixture in porous media. The functional approach and technique for decoupling the correlations explained by Klyatskin [4] is used. The second part of the paper studies the process of transport in porous media of two immiscible incompressible fluids in the framework of the Buckley—Leverett model. A linear equation is obtained for the joint probability density of the solution of the stochastic quasilinear transport equation and its derivative. An infinite chain of equations for the moments of the solution is obtained. A scheme of approximate closure is proposed, and the solution of the approximate equations for the mean concentration is compared with the exactly averaged concentration.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 127–136, September–October, 1985.We are grateful to A. I. Shnirel'man for pointing out the possibility of obtaining an averaged equation in the case of a velocity distribution in accordance with a Cauchy law.     

Percolation with a limiting gradient in a medium with random inhomogeneities

The coverage of a medium by percolation and the effective permeability of a medium with stagnant zones are determined. It is shown that effective permeability is a function of external conditions, particularly the average pressure gradient. Three-, two-, and one-dimensional flows are discussed. The theory of overshoots of random functions and fields beyond a prescribed level [1, 2] is used for the investigation. Overshoots of elements of the percolation field in media with random inhomogeneities are studied. Overshoots of energy being dissipated in a volume are discussed in particular; this permits an approximate determination of the coverage of an inhomogeneous porous medium by migration during percolation with a limiting gradient, i.e., in the case of formation of stagnant zones chaotically disseminated in the flow region.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 159–165, September–October, 1970.The authors thank V. M. Entov for discussing the article and useful comments.     

Some approximate solutions of filtration problems in a fissured-porous medium

The equations for the filtration of a fluid in a fissured-porous medium [1] under the assumption that the permeability of the porous blocks is negligible in comparison with the permeability of the cracks and that the porosity of the cracks is negligible in comparison with the porosity of the blocks may be written in the form Here p₁ is the pressure in the cracks, p₂ is the pressure in the porous blocks, is the characteristic lag time, , is the piezoconductivity coefficient. We shall consider the approximate solutions of this system of equations in the case of filtration to a well which penetrates a fissured-porous stratum of thickness h and begins to operate at the moment t=0 with the flow rate Q.The author wishes to tank V. N. Nikolaevskii for discussions of the study.     

Macroscale Properties of Porous Media from a Network Model of Biofilm Processes

We demonstrate how a network model can predict porosity and permeability changes in a porous medium as a result of biofilm buildup in the pore spaces. A biofilm consists of bacteria and extracellular polymeric substances (EPS) bonded together and attached to a surface. In this case, the surface consists of the walls of the porous medium, which we model as a random network of pipes.Our model contains five species. Four of these are bacteria and EPS in both fluid and adsorbed phases. The fifth species is nutrient, which we assume to reside in the fluid phase only. Bacteria and EPS transfer between the adsorbed and fluid phases through adsorption and erosion or sloughing. The adsorbed species influence the effective radii of the pipes in the network, which affect the porosity and permeability.We develop a technique for integrating the coupled system of ordinary and partial differential equations that govern transport of these species in the network. We examine ensemble averages of simulations using different arrays of pipe radii having identical statistics. These averages show how different rate parameters in the biofilm transport processes affect the concentration and permeability profiles.     

Vorticity of the velocity field for flow in a porous medium with random inhomogeneities

The vorticity field of the flow velocity in a porous medium with random inhomogeneities is considered in the correlation approximation of perturbation theory. The correlation tensor of the vorticity, the correlation between the vorticity and the permeability field, and the circulation of the velocity are calculated for three- and two-dimensional flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 157–160, July–August, 1982.     

高温高压气体在低含水率介质中迁移的数值模拟

针对气液两相非等温渗流模型高度非线性的特点,发展了适宜的数值离散方法。根据相态转换准则和控制方程的性质,采用最低饱和度法简化算法。空间离散方面,使用有限体积法;时间离散方面,设计了一套包含合理求解顺序的Picard迭代法,解决了方程组强耦合的问题。利用上述数值方法对高温高压气体的迁移行为进行数值模拟,证明了气体在低含水率介质和等效孔隙度的干燥介质内的运动基本一致,并分析了空腔内的气液相态转变过程。在此基础上,研究了多孔介质孔隙度和渗透率对气体压强演化和示踪气体迁移的影响。研究表明,孔隙度越小（相同渗透率）、渗透率越高（相同孔隙度）,示踪气体的迁移距离越远,并给出了估算不同孔隙度和渗透率下迁移距离的半经验公式。     

Transport equations for elastic and other waves in random media

We derive and analyze transport equations for the energy density of waves of any kind in a random medium. The equations take account of nonuniformities of the background medium, scattering by random inhomogeneities, polarization effects, coupling of different types of waves, etc. We also show that diffusive behavior occurs on long time and distance scales and we determine the diffusion coefficients. The results are specialized to acoustic, electromagnetic, and elastic waves. The analysis is based on the governing equations of motion and uses the Wigner distribution.     

Elasticoplastic conditions of the filtration of a liquid in porous media

The article discusses questions in the theory of filtration in porous media, taking account of elastic, elasticoplastic, and plastic deformations. Parameters are introduced to evaluate irreversible effects in petroleum- and water-bearing strata, i.e., coefficients of the change in the porosity and the permeability. Equations are derived for filtration under unsteady-state and steady-state working conditions of wells and galleries. Two limiting cases, which allow analytical solutions, are separated out. In the general case, the equations of elasticoplastic filtration conditions are solved on an electronic computer. The numerial calculations show that the predominating effect results from taking account of the irreversible change in the permeability, depending on the change of the pressure in the stratum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 84–90, September–October, 1973.The author is grateful to V. N. Nikolaevskii for his evaluation of the work.     

Statistical theory of radiative transfer in layered media

An equation for the probability density of the wave intensity which takes into account absorption, is obtained with a help of the invariant imbedding method. The limiting case when the medium occupies a half-space, is considered. The field intensity is found for the case of a source inside the medium. The conditions of applicability of the linear theory or radiative transfer are obtained. Numerical solutions of the equations corresponding to the statistical theory of radiative transfer in a layered medium with random inhomogeneities are discussed.     

Modeling the Short-Term Damageability of a Transversely Isotropic Piezoelectric Material

A short-term microdamage theory for porous transversely isotropic piezoelectric materials is set forth. Microdamages are modeled by pores. The fracture criterion for a microvolume of a transversely isotropic medium is assumed to have the Huber–Mises form. The ultimate strength is a random function of coordinates with an exponential or Weibull distribution. The stress–strain distribution and effective properties of the material are determined from the stochastic electroelastic equations. The deformation and microdamage equations are closed by the porosity balance equations. For various values of electric intensity, the microdamage–macrodeformation relationships and deformation curves are plotted. The effect of electric intensity on the microdamage of piezoelectric materials is studied     

Internal rotation in the hydrodynamics of weakly conducting dielectric suspensions

Hydrodynamic phenomena in weakly conducting single-phase media due to interphase electric stresses are reviewed in [1]. In the present paper, a model is constructed of a dielectric suspension with body couples due to the field acting on free charges distributed on the surface of the particles of the suspension. Averaging of the microscopic fields yields macroscopic equations for the field and the polarization of the dielectric suspension with allowance for the finite relaxation time of the distribution of the free charge on the phase interface. The developed model is used to consider the occurrence of spontaneous rotation of a dielectric cylinder in a weakly conducting suspension in the presence of an electric field; compared with the case of single-phase media [2], this is characterized by a significant reduction in the threshold intensity of the electric field with increasing concentration of the particles [3]. In the present model of a dielectric suspension, the destabilization of the cylinder is due to the occurrence of rotations of the particles of the suspension due to the interaction between the polarization and the motion of the medium. The relaxation equation for the polarization for the given model is analogous to the corresponding equation for media which can be magnetized [4–6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 86–93, March–April, 1980.     

Instability of two streaming conducting and dielectric bounded fluids in porous medium under time-varying electric field

In this paper, we consider the instability of the interface between two superposed streaming conducting and dielectric fluids of finite depths through porous medium in a vertical electric field varying periodically with time. A damped Mathieu equation with complex coefficients is obtained. The method of multiple scales is used to obtain an approximate solution of this equation, and then to analyze the stability criteria of the system. We distinguish between the non-resonance case, and the resonance case, respectively. It is found, in the first case, that both the porosity of porous medium, and the kinematic viscosities have stabilizing effects, and the medium permeability has a destabilizing effect on the system. While in the second case, it is found that each of the frequency of the electric field, and the fluid velocities, as well as the medium permeability, has a stabilizing effect, and decreases the value of the resonance point, while each of the porosity of the porous medium, and the kinematic viscosities has a destabilizing effect, and increases the value of the resonance point. In the absence of both streaming velocities and porous medium, we obtain the canonical form of the Mathieu equation. It is found that the fluid depth and the surface tension have a destabilizing effect on the system. This instability sets in for any value of the fluid depth, and by increasing the depth, the instability holds for higher values of the electric potential; while the surface tension has no effect on the instability region for small wavenumber values. Finally, the case of a steady electric field in the presence of a porous medium is also investigated, and the stability conditions show that each of the fluid depths and the porosity of the porous medium ɛ has a destabilizing effect, while the fluid velocities have stabilizing effect. The stability conditions for two limiting cases of interest, the case of purely fluids), and the case of absence of streaming, are also obtained and discussed in detail.     

A stochastic model for filtration of particulate suspensions with incomplete pore plugging

A population balance model for particulate suspension transport with capture of particles by porous medium accounting for complete and incomplete plugging of pores by retained particles is derived. The model accounts for pore space accessibility, due to restriction on finite size particle movement through the overall pore space, and for particle flux reduction, due to transport of particles by the fraction of the overall flux. The novel feature of the model is the residual pore conductivity after the particle retention in the pore and the possibility of one pore to capture several particles. A closed system of governing stochastic equations determines the evolution of size distributions for suspended particles and pores. Its averaging results in the closed system of hydrodynamic equations accounting for permeability and porosity reduction due to plugging. The problem of deep bed filtration of a single particle size suspension through a single pore size medium where a pore can be completely plugged by two particles allows for an exact analytical solution. The phenomenological deep bed filtration model follows from the analytical solution.     

The spherical problem of the filtration of a liquid and a gas with a change in the permeability of rocks in accordance with an exponential law

The article considers the problem of the filtration of liquids (or gases), pumped through a borehole at a constant rate with elastic filtration conditions. The permeability of the stratum is assumed to be an exponential function of the coordinates. The viscosities of the injected and displaced liquids are assumed to be different. To increase the capacity of strata, i.e., of collectors used for the burial of industrial waste flows and gases, various methods are employed to increase the fracturing and the permeability of the rocks (hydro-pulse techniques, explosions, and other methods). As a result of this, a spherical region is formed in the rocks, in which the permeability varies along the radius. The character of this change is well described by an exponential function. The pumping of waste flows or industrial gases into such a cavity leads to the displacement of the stratum liquid (or gas). The problem of the displacement of one liquid by another liquid not miscible with it under rigid filtration conditions was first discussed in [1–5]. Here a study was made of a region of finite dimensions, bounded by two boundaries, with given pressures or mass flow rates (the linear and axisymmetric flow problems). The permeability of the stratum was assumed to be independent of the coordinates. A special characteristic of these problems is the fact that it is impossible to consider unbounded or semi-bounded filtration conditions in them since, under rigid filtration conditions, the condition of bounded character of the pressure (the head) is not satisfied at infinity. Elastic filtration conditions for two immiscible liquids were first discussed in [6], and later in [7, 8] and other reports. Here an investigation was made of the linear and axisymmetric problems for an unbounded region. In [9, 10] solutions are given to some problems with spherical symmetry for an unbounded region, with rigid filtration conditions and a jumpwise change of the permeability along the radius. In the problems of [6–10] the condition of the bounded character of the pressure is satisfied. In [11] the case of a hyperbolic change in the permeability of the rocks is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 42–51, November–December, 1974.     

Statistical characteristics of the motion of a fluid particle in a medium with random porosity

In the study of flow of a neutral admixture in a porous medium, it is most often assumed in the stochastic formulation that the porosity is constant and a determinate quantity, and the velocity is a random function [1–4]. The velocity distribution is usually regarded as known. Flow in a porous medium with random porosity has been studied to a far lesser extent. We note [5], which studies the averaged equations obtained within the framework of the correlation approximation. We consider the model problem of one-dimensional motion of a fluid particle (position of the front for flow of a neutral admixture in a porous medium) in a medium with random porosity. For a particular form of random porosity field, expressions are obtained for the one- and two-point densities of the distribution of the position of the particle. A study is made of the dependences of the first four moments and the correlation function of the position of the particle as functions of the time. It is shown that for large values of the time the motion of the particle is asymptotically similar to Brownian motion. It is shown by means of numerical modeling that the results obtained transfer to the case of an arbitrary random porosity field. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 59–65, November–December, 1986.     

Multi-fluid model of suspension filtration in a porous medium

In the framework of a three-fluid approach, a new model of suspension filtration in a porous medium is constructed with account for the formation of a dense packing of trapped particles with finite permeability and porosity. The following three continua are considered: the carrier fluid, the suspended particles, and the deposited particles. For a one-dimensional transient flow of suspension, a system of equations for the concentrations of the suspended and deposited particles, the suspension velocity, and the pressure is constructed. Two cases of the flow in a porous medium are considered: plane and radial. Numerical solution is found using a finite-difference method. Numerical calculations are shown to be in agreement with an analytical solution for the simplest case of filtration with a constant velocity and constant porosity and permeability. A comparison is performed with the classic filtration models for five sets of experimental data on the contamination of a porous sample. It is shown that near the inlet boundary, where an intense deposition of particles takes place, the new model describes the concentration profile of the deposited particles more accurately than the classical model.     

Flow through isotropic granular porous media

A generalization of the Navier-Stokes equation is developed to include laminar flow through a rigid isotropic granular porous medium of spatially varying permeability. The model is based on a theory of interspersed continua and the mean geometrical properties of an idealized granular porous microstructure. The derived momentum transport equations are applicable to granular porous media over the entire porosity range from zero through unity. No restriction with respect to flow velocity is imposed, except for the assumption of laminar flow within the pores. The results provide useful and versatile equations and substantiate many of the empirical equations currently in use. One of the major advantages of the generalized momentum equation is its adaptability to numerical simulation.     

Porosity–Permeability Relations for Evolving Pore Space: A Review with a Focus on (Bio-)geochemically Altered Porous Media

Reactive transport processes in a porous medium will often both cause changes to the pore structure, via precipitation and dissolution of biomass or minerals, and be affected by these changes, via changes to the material’s porosity and permeability. An understanding of the pore structure morphology and the changes to flow parameters during these processes is critical when modeling reactive transport. Commonly applied porosity–permeability relations in simulation models on the REV scale use a power-law relation, often with slight modifications, to describe such features; they are often used for modeling the effects of mineral precipitation and/or dissolution on permeability. To predict the reduction in permeability due to biomass growth, many different and often rather complex relations have been developed and published by a variety of authors. Some authors use exponential or simplified Kozeny–Carman relations. However, many of these relations do not lead to fundamentally different predictions of permeability alteration when compared to a simple power-law relation with a suitable exponent. Exceptions to this general trend are only few of the porosity–permeability relations developed for biomass clogging; these consider a residual permeability even when the pore space is completely filled with biomass. Other exceptions are relations that consider a critical porosity at which the porous medium becomes impermeable; this is often used when modeling the effect of mineral precipitation. This review first defines the scale on which porosity–permeability relations are typically used and aims at explaining why these relations are not unique. It shows the variety of existing approaches and concludes with their essential features.