AN IMPROVED HOMOGENIZATION RESULT FOR IMMISCIBLE COMPRESSIBLE TWO – PHASE FLOW IN POROUS MEDIA

The paper deals with a degenerate model of immiscible compressible two-phase flow in heterogeneous porous media. We consider liquid and gas phases (water and hydrogen) flow in a porous reservoir, modeling the hydrogen migration through engineered and geological barriers for a deep repository for radioactive waste. The gas phase is supposed compressible and obeying the ideal gas law. The flow is then described by the conservation of the mass for each phase. The model is written in terms of the phase formulation, i.e. the liquid saturation phase and the gas pressure phase are primary unknowns. This formulation leads to a coupled system consisting of a nonlinear degenerate parabolic equation for the gas pressure and a nonlinear degenerate parabolic diffusion-convection equation for the liquid saturation, subject to appropriate boundary and initial conditions. The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system. The aim of this paper is to extend our previous results to the case of an ideal gas. In this case a new degeneracy appears in the pressure equation. With the help of an appropriate regularization we show the existence of a weak solution to the studied system. We also consider the corresponding nonlinear homogenization problem and provide a rigorous mathematical derivation of the upscaled model by means of the two-scale convergence. 2010 Mathematics Subject Classification. Primary: 35B27, 35K65, 74Q15, 76M50; Secondary: 76S05, 76T10.


BRAHIM AMAZIANE, LEONID PANKRATOV AND ANDREY PIATNITSKI
1. Introduction.The modeling of displacement process involving two immiscible fluids is of considerable importance in groundwater hydrology and reservoir engineering such as petroleum and environmental problems.More recently, modeling multiphase flow received an increasing attention in connection with gas migration in a nuclear waste repository and sequestration of CO 2 .
In this paper, we focus our attention on the modeling of immiscible compressible two-phase flow through heterogeneous reservoirs in the framework of the geological disposal of radioactive waste.The long-term safety of the disposal of nuclear waste is an important issue in all countries with a significant nuclear program.One of the solutions envisaged for managing waste produced by nuclear industry is to dispose the radioactive waste in deep geological formations chosen for their ability to delay and to attenuate possible releases of radionuclides in the biosphere.Repositories for the disposal of high-level and long-lived radioactive waste generally rely on a multi-barrier system to isolate the waste from the biosphere.The multibarrier system typically comprises the natural geological barrier provided by the repository host rock and its surroundings and an engineered barrier system, i.e .engineered materials placed within a repository, including the waste form, waste canisters, buffer materials, backfill and seals, for more details see for instance [43].An important task of the safety assessment process is the handling of heterogeneities of the geological formation.
In the frame of designing nuclear waste geological repositories, a problem of possible two-phase flow of water and gas, mainly hydrogen, appears, for more details see for instance [43].Multiple recent studies have established that in such installrations important amounts of gases are expected to be produced in particular due to the corrosion of metallic components used in the repository design, see e.g.[27,42] and the references therein.The French Agency for the Management of Radioactive Waste (Andra) [11] is currently investigating the feasibility of deep geological disposal of radioactive waste in an argillaceous formation.A question related to the long-term performance of the repository concerns the impact of the hydrogen gas generated in the wastes on the pressure and saturation fields in the repository and the host rock.
During recent decades mathematical analysis and numerical simulation of multiphase flows in porous media have been the subject of investigation of many researchers owing to important applications in reservoir simulation.There is an extensive literature on this subject.We will not attempt a literature review here but will merely mention a few references.Here we restrict ourselves to the mathematical analysis of such models.We refer, for instance, to the books [14,23,26,29,36,38,44] and the references therein.The mathematical analysis and the homogenization of the system describing the flow of two incompressible immiscible fluids in porous media is quite understood.Existence, uniqueness of weak solutions to these equations, and their regularity has been been shown under various assumptions on physical data; see for instance [3,14,15,21,24,25,23,29,41] and the references therein.A recent review of the mathematical homogenization methods developed for incompressible immiscible two-phase flow in porous media and compressible miscible flow in porous media can be viewed in [4,37,38].We refer for instance to [16,17,18,19,20,34,35] for more information on the homogenization of incompressible, single phase flow through heterogeneous porous media in the framework of the geological disposal of radioactive waste.
However, as reported in [9], the situation is quite different for immiscible compressible two-phase flow in porous media, where, only recently few results have been obtained.In the case of immiscible two-phase flows with one (or more) compressible fluids without any exchange between the phases, some approximate models were studied in [30,31,32].Namely, in [30] certain terms related to the compressibility are neglected, and in [31,32] the mass densities are assumed not to depend on the physical pressure, but on Chavent's global pressure.In the articles [22,33,39,40], a more general immiscible compressible two-phase flow model in porous media is considered for fields with a single rock type and [9] treated the case with several types of rocks.In [4,10] homogenization results were obtained for water-gas flow in porous media using the phase formulation, i.e.where the phase pressures and the phase saturations are primary unknowns.
Let us also mention that, recently, a new global pressure concept was introduced in [5,7] for modeling immiscible, compressible two-phase flow in porous media without any simplifying assumptions.The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic equation (the global pressure equation) and a nonlinear diffusion-convection one (the saturation equation).This new formulation is fully equivalent to the original phase equations formulation, i.e.where the phase pressures and the phase saturations are primary unknowns.For this model, an existence result is obtained in [8] and homogenization results in [6].
Let us note that all the aforementioned works are restricted to the case where the gas density is bounded from below and above, contrarily to the present work.This assumption is too restrictive for some realistic problems, such as gas migration through engineered and geological barriers for a deep repository for radioactive waste.In this case the gas obeys the ideal gas law, i.e. the equation of state is given by g (p) def = σp where g is the gas density, p g is the gas pressure and σ is a given constant.Then a new degeneracy appears in the evolution term of the gas pressure equation.In this paper we extend our previous results obtained in [4] to the more complex case of an ideal gas which is more reasonable in gas reservoir engineering.The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system.To obtain these results we elaborated a new approach based on the ideas from [4] and regularization.
The rest of the paper is organized as follows.In Section 2 we describe the physical model and formulate the corresponding mathematical problem.We also provide the assumptions on the data.
The goal of Section 3 is to prove the existence result for the corresponding system of equations.The proof is divided into a number of steps.In subsection 3.1 we consider an auxiliary δ-problem where the gas density admits a positive lower bound.The existence result for the δ-problem is given in subsection 3.2.In subsection 3.3 we obtain a number of a priori estimates for a solution of δ-problem.Then in subsection 3.4 we prove a compactness result adapted to our model.Finally, in subsection 3.5 we pass to the limit in δ-problems and complete the proof of the existence.
Section 4 is devoted to the corresponding homogenization problem.In subsection 4.1 we introduce the model with a periodic microstructure.We assume that both porosity and absolute permeability tensor are periodic rapidly oscillating functions.Then subsection 4.2 we formulate the homogenization result.This result is proved in subsection 4.3.The proof is based on the two-scale convergence technique.Our analysis relies essentially on a compactness result [4] which is rather involved due to the degeneracy and the nonlinearity of the system.The last section is followed by some concluding remarks.
2. Formulation of the problem.We consider an immiscible compressible twophase flow process in a porous reservoir Ω ⊂ R d (d = 1, 2, 3) which is a bounded Lipschitz domain.The time interval of interest is (0, T ).We focus here on the particular case of water and gas phases, but the consideration below is also valid for a general wetting phase and a non-wetting phase.Let Φ = Φ(x) be the porosity of Ω; K = K(x) be the absolute permeability tensor of Ω; w , g are the densities of water and gas, respectively; S w = S w (x, t), S g = S g (x, t) are the saturations of water and gas in Ω × (0, T ); k r,w = k r,w (S w ), k r,g = k r,g (S g ) are the relative permeabilities of water and gas; p w = p w (x, t), p g = p g (x, t) are the pressures of water and gas in Ω × (0, T ).
In what follows, for the sake of presentation simplicity we neglect the source terms.Then the conservation of mass of each phase is described by (see, e.g., [23,26,36]): where Ω T def = Ω × (0, T ) with T > 0, the velocities of water and gas q w , q g are defined by Darcy-Muskat's law: Here g, µ w , µ g are the gravity vector and the viscosities of the water and gas, respectively.
From now on we assume that the density of the water is constant, which for the sake of simplicity will be taken equal to one, i.e. w (p w ) = Const = 1, and the gas density g obeys the ideal gas law and is given by the following function: Here σ, max , and p max are positive constants.Note that g is a continuous monotone and non-negative function.
The model is completed as follows.By the definition of saturations, one has S w + S g = 1 with S w , S g 0.
We set: ) Then the curvature of the contact surface between the two fluids links the jump of pressure of two phases to the saturation by the capillary pressure law: P c (S) = p g − p w with P c (s) < 0 for all s ∈ [0, 1] and P c (1) = 0, (7) where P c (s) denotes the derivative of the function P c (s).Now due to (6) and the assumption on the water density, we rewrite the system (1) as follows: where λ g (S) := λ g (1 − S) and The system (8) have to be completed by appropriate boundary and initial conditions.
Boundary conditions.: We suppose that the boundary ∂Ω consists of two parts Γ inj and Γ imp such that Γ inj ∩ Γ imp = ∅, ∂Ω = Γ inj ∪ Γ imp .The boundary conditions are given by: where the velocities q w , q g are defined in (2), (3).Initial conditions.: The initial conditions read: p w (x, 0) = p 0 w (x) and p g (x, 0) = p 0 g (x) in Ω. (11) Notice that from (10) and (7) it follows that S = 1 on Γ inj × (0, T ).The initial condition for S is uniquely defined by the equation Then according to (9) the initial condition for Θ reads Remark 1.It is important to underline that in the earlier works (see, e.g., [4,9,10,30,31,32,33]) it was assumed that the gas density admits a strictly positive lower bound: 2.1.A fractional flow formulation.In the sequel, we use a formulation obtained after transformation using the concept of the so called global pressure.In the case of incompressible two-phase flow this concept was introduced for the first time in [12,13].Following [14,23], see also [26], we first recall the definition of the global pressure.It plays a crucial role, in particular, for compactness results.The idea of introducing the global pressure is as follows.We want to replace the water-gas flow by a flow of a fictive fluid obeying the Darcy law with a non-degenerating coefficient.Namely, we are looking for a pressure P and the coefficient γ(S) such that γ(S) > 0 for all S ∈ [0, 1], and Then the global pressure, P, is defined by: the functions G w (s) and G g (s) will be introduced later on, in (19), (20).Now it is easy to see that λ w (S)∇p w + λ g (S)∇p g = λ(S)∇P + λ g (S)∇G g (S) + λ w (S)∇G w (S) , where λ(s) Then γ(S) = λ(S).By construction, λ(S) > 0 for all S ∈ [0, 1] (see the condition (A.5) below).Thus the relation ( 15) is established.Now we specify the functions G w , G g .We define G g as follows: The functions G w are then defined by Notice that from (19), (20) we get: where Now we link the capillary pressure and the mobilities.In a standard way (see, e.g., [40] or [9] for more details) we obtain the following identity: where Ii is also convenient to introduce the following function β: where the function α is defined in (22).Notice that by the definition of the global pressure, ( 24), (25), and by the boundedness of λ w , λ g (see the condition (A.5) below) the following relations holds: In order to complete this section, let us calculate the value of the global pressure function P on Γ inj .In what follows (see condition (A.4) below) we assume that the capillary pressure function P c satisfies the condition P c (1) = 0. Then from (10) 1 we have that the saturation S equals one on Γ inj .Now the definition of the global pressure (16) implies that the function P on Γ inj is a constant which we denote by P 1 .
(A.4) The capillary pressure function The functions λ w , λ g belong to the space C([0, 1]; R + ) and satisfy the following properties: 25) is a Hölder continuous function of order θ with θ ∈ (0, 1) on the interval [0, β( ].There exists a positive constant C β such that for all s 1 , s 2 ∈ [0, β(1)] the following inequality holds: (A.8) The initial data for the pressures are such that p 0 g , p 0 w ∈ L 2 (Ω).(A.9)The initial data for the saturation is such that S 0 ∈ L ∞ (Ω) and 0 S 0 1 a.e. in Ω.The assumptions (A.1)-(A.9)are classical for two-phase flow in porous media.
3. Existence result.In order to define a weak solution of the above problem, we introduce the following Sobolev space: The space H 1 Γinj (Ω) is a Hilbert space.The norm in this space is given by u Theorem 3.1.Let assumptions (A.1)-(A.9)be fulfilled.Then there exist functions p g , p w , S such that: (I): where the function Θ is defined in ( 9); S = 1 on Γ inj .(II): the maximum principle holds: (III): For any ϕ w , ϕ g ∈ C 1 ([0, T ]; H 1 (Ω)) satisfying ϕ w = ϕ g = 0 on Γ inj × (0, T ) and ϕ w (x, T ) = ϕ g (x, T ) = 0, we have: + with Θ defined in ( 9), and P c (S) = p g − p w .(IV): The initial conditions are satisfied in a weak sense as follows: (36) Furthermore, we have and with S 0 and Θ 0 defined in (12) and (13), respectively.
The proof of Theorem 3.1 is divided into a several steps.It is based on a auxiliary existence result for the system obtained by approximation of the initial degenerate gas density g by a family of functions { δ g } δ>0 that admit a positive lower bound.For such kind of system the desired existence result is proved in [9,30,31,32,33].This result is formulated in subsection 3.1.Using the weak formulation of the regularized problem and the uniform in δ estimates for its solution, we, finally, prove Theorem 3.1.

Auxiliary δ-problem.
In this subsection we approximate the function g by a family of functions { δ g } δ>0 that admit a positive lower bound.For each δ > 0 we set: Here σ, max , p max are positive constants.
In addition to (8), consider the following family of problems: System (40) have to be completed with the corresponding boundary and initial conditions.
Boundary conditions.: The boundary conditions read where the velocities q δ w , q δ g are given by: Initial conditions.: The initial conditions read: The remaining part of the Section is organized as follows.First, in subection 3.2 we recall the existence result for the system (40).Then we obtain the uniform in δ estimates for the solution of δ-problem (40).In subection 3.4 we formulate the compactness and convergence results which we use in the proof of Theorem 3.1.

3.2.
An existence result for the δ-problem.The goal of this subsection is to recall the existence result for the δ-problem (40).First, we reformulate the condition (A.3) from subsection 2.2 in order to adapt it to our δ-problem.For this problem the condition (A.3) becomes: is given by ( 39).Now we are in position to formulate the existence result to δ-problem (40).It reads.

3.3.
A priori estimates for a solution of δ-problem (40).We start this subsection by obtaining the energy equality for δ-problem (40).The following result holds.
Lemma 3.3 (Energy equality for δ-problem).Let p δ w , p δ g , S δ be a solution to (40).Then in the sense of distributions.Here .
Notice that in the previous works (see, e.g., [4,9,10,30,31,32,33]), the function . However, in our case, with such a definition we do not have uniform in δ lower bound for the function R δ .Thus, we have to modify the definition of R g , subtracting an appropriate constant C = C(δ).The properties of the functions R δ g , R δ , and E δ are given in: , and E δ be the functions defined by ( 55), (56).Then (i) The function R δ is negative and bounded from below, that is The function E δ is bounded from below.Namely, Proof of Lemma 3.4.Using the definition of the gas density δ g given by (39), it is easy to calculate that (59) Consider now the function R δ .Due to (39) and (59), we have: The last formula immediately implies (57).Now (58) follows easily from (57) and the estimate: This completes the proof of Lemma 3.4.
In order to formulate a priori estimates for the solution to δ-problem (40), we remark first that the global pressure P δ for the problem under consideration can be introduced in a way similar to one used in subsection 2.1 above.Then the desired a priori estimates for the solution of the δ-problem can be easily derived from Lemmata 3.3, 3.4, and the equality (see subsection 2.1 for more details): where the function b(s) is defined in (24).
The following result holds.
Lemma 3.5.Let p δ w , p δ g , S δ be a solution to (40), the global pressure P δ is defined in (16), and the function β(s) is defined in (25).Then (65) Here C does not depend on δ.
Proof of Lemma 3.5.Integrating (54) over the interval (0, T ), we get: Let us estimate now the right-hand side of (66) from above.Due to the definition of the function E δ , (55), and the initial conditions (44) we have that where S 0 = S 0 (x) is the initial condition of the saturation function (see condition (A.9) in Section 2.2).Now from condition (A.1) and the maximum principle (48), we easily obtain that where |Ω| stands for the measure of the domain Ω.Therefore, where C 0 is a constant which only depends on max P c (S), and the constant φ + .
3.5.Proof of Theorem 3.1.We begin this subsection by studying the regularity properties of solution to (8).

3.5.1.
Regularity properties of a solution to system (8).Taking into account the lower semi-continuity of the norm, by Lemma 3.5, we obtain: Now we set: We also recall the relation (23): Then, taking into account (78), (80), and the last relation we obtain that the functions p w , p g defined in (81) are such that Thus properties ( 30)-( 31) are established.The maximum principle (33) follows immediately from (48) and (72).Finally, the interpretation of the initial conditions can be done as in [40] (see also [9]).
Taking into account (72), one easily gets: We then recall that Then the third term on the left-hand side of (49) takes the form: Now taking into account the convergence (72), (74), and (75), we obtain that Returning now to the water pressure function p w , we finally get: Considering (72), one can check that the fourth term of (49) satisfies the relation Thus, the saturation equation ( 34) is derived.We turn to (50) with ϕ g ∈ C 1 ([0, T ]; H 1 (Ω)), ϕ g = 0 on Γ inj × (0, T ), and ϕ g (x, T ) = 0.
Taking into account (76), one easily gets: where Considering the definition of the functions g and δ g (see ( 4) and (39), respectively) we have In order to pass to the limit in the third term of (50) we recall that (see relations (27)) Then Now taking into account the convergence results (72), (74), and (75) we obtain that Returning now to the gas pressure function p g , we finally get: (91) Thus the gas pressure equation ( 35) is obtained.Theorem 3.1 is proved.
4. Homogenization result.In this Section we consider the problem describing a reservoir with a periodic microstructure.Then in the model considered in the previous sections one has rapidly oscillating porosity function and absolute permeability tensor.Our goal is to prove the homogenization result for this model.The convergence of the homogenization process is justified by the technique of two-scale convergence [2].4.1.Formulation of the microscopic problem.In this section, we present the mathematical model describing water-gas flow in a periodically heterogeneous porous medium.As above we suppose that the gas density vanishes as the gas pressure is zero.For simplicity, we assume no source/sink terms.
We consider a bounded Lipschitz domain Ω ⊂ R d (d = 1, 2, 3) with a periodic microstructure.The microscopic length scale that represents the ratio between the cell size to the size of the whole region Ω, is denoted by ε.We assume that 0 < ε 1 is a small parameter tending to zero.We denote by Y def =(0, 1) d the periodic cell.Let Φ ε (x) = Φ(x/ε) be the porosity of Ω; K ε (x) = K(x/ε) be the absolute permeability tensor of Ω; S ε def = S ε w (x, t), is the saturations of water in Ω × (0, T ); p ε w = p ε w (x, t), p ε g = p ε g (x, t) are the pressures of water and gas in Ω × (0, T ), respectively; System (8), in the case of a periodic porous medium, takes the form (93) System (92) has to be equipped with appropriate boundary and initial conditions.Boundary conditions.: We suppose that the boundary ∂Ω consists of two parts Γ inj and Γ imp such that Γ inj ∩ Γ imp = ∅, ∂Ω = Γ inj ∪ Γ imp .The boundary conditions are given by: where the velocities q ε w , q ε g are defined as follows: Initial conditions.: The initial conditions read: Let us formulate the main assumptions on the data.First, we replace conditions (A.1), (A.2) from Section 2.2 with the following assumptions: there exist positive constants K − , K + such that We also suppose that conditions (A.3)-(A.9)from subsection 2.2 hold true.

Notational convention.
In what follows C, C 1 , .. denote generic constants that do not depend on ε.

4.2.
Statement of the homogenization result.We study the asymptotic behavior of the solution to problem (92), ( 94), (96) as ε → 0. In particular, we are going to show that the effective model reads: where S, P w , P g denote the homogenized water saturation, water pressure, and gas pressure, respectively.Φ denotes the mean value of the function Φ over the cell Y .K is the homogenized tensor with the entries K ij defined by: where the function ξ j is a Y -periodic solution of the following local problem: with e j being the j-th coordinate vector.
Remark 2. The homogenized system (100) generalizes the result obtained earlier in [4] in two ways.First, this system allows the gas density to degenerate.In addition, this system is written in terms of the homogenized phase pressures P w , P g and not in terms of the homogenized global pressure and water saturation as was done in [4] (see (3.1)).
System (100) has to be completed with the following boundary and initial conditions.
Boundary conditions.: The boundary conditions are given by: P g (x, t) = P w (x, t) = 0 on Γ inj × (0, T ); where the velocities q ε w , q ε g are defined as follows: Initial conditions.: The initial conditions read: P w (x, 0) = p 0 w (x) and P g (x, 0) = p 0 g (x) in Ω.
(105) 4.3.Proof of Theorem 4.3.The rigorous justification of the homogenization process relies on the two-scale convergence approach, see, e.g., [2].For the reader's convenience, we recall the definition of the two-scale convergence.
C, and for any test function ϕ ∈ C ∞ (Ω T ; C # (Y )) the following relation holds: This convergence is denoted by v ε (x, t) 2s v(x, y, t).
The homogenization result reads.The proof is divided into a number of steps.

4.3.1.
A priori estimates for solutions to problem (92).In this section we derive the a priori estimates for problem (92).For any ε > 0, we consider the following ε, δ-problem: where the family of functions { δ g } δ>0 is defined in (39) and The (ε, δ)-problem is completed by the boundary and initial conditions (94) and (96), respectively.The energy equality for problem (106) can be obtained as in Section 3.3.
Then following the lines of Section 3.3 one can prove the following statement which is similar to that of Lemma 3.5.Lemma 4.5.Let p ε,δ w , p ε,δ g be a solution to (106), the global pressure P ε,δ be defined in (16), and the function β(s) be defined in (25).Then Here C does not depend on ε, δ.
Now, as in Section 3.5.1,we conclude that where C is a constant that does not depend on ε, δ.The uniform estimates for the time derivatives of the functions Φ ε Θ ε and Φ ε S ε can be derived from (92) using (115).These estimates read: where C is a constant that does not depend on ε.

4.3.2.
Compactness and convergence results.First, we recall the following compactness result established in [4].
Remark 3. In the formulation of the above compactness lemma the periodicity of Φ can be replaced with the assumption that Φ ε 1 weakly in L 2 (Ω), as ε → 0. Now we turn to the compactness result for the family {Θ ε } ε>0 .
A similar result holds for the set {S ε } ε>0 .
Summarizing the above statements yields.
Lemma 4.7.There exist a function S with 0 S 1 a.e. in Ω T and a function P ∈ L 2 (0, T ; H 1 (Ω)) such that up to a subsequence: The Proof of Lemma 4.7 relies on the arguments similar to those used in the proof of Lemma 4.8 in [4].
It is easy to justify the passage to the two-scale limit in the temporal terms using the convergence results (117) and (120) from Lemma 4.7 as it was done, for example, in [4].Namely, let ϕ 0 ∈ D(Ω T ).The first two terms in (98) become: Now we pass to the limit on the right-hand side of (121).Taking into account (117), we have that lim For any η > 0, we introduce the family of functions {S ε,η } defined by: S ε,η def = min {(1 − η), max(η, S ε )} .
(126) Therefore, V η w = ξ(y) ∇ x P + ∇ x G w (S) − g (127) for all (x, t) ∈ Ω T such that S ∈ (η, 1 − η).Here ξ ∈ R d is a vector with the components ξ j that are the solutions of the auxiliary problem (102).Since η is an arbitrary positive number, representation (127) is valid for all (x, t) such that S ∈ (0, 1).In particular, V η does not depend on η: V η w = V w .This leads to the following equation: Finally, with the help of our a priori estimates we deduce in a standard way that K ε λ w (S ε ) ∇p ε w − g 2s K(y)λ w (S) I + ∇ y ξ(y) where I is the unit matrix.This allows us, with the help of (122), to obtain the weak formulation of the homogenized saturation equation (100) 2 .
The derivation of the weak formulation for the homogenized gas pressure equation can be done in a similar way.This completes the proof of Theorem 4.3.
5. Concluding remarks.We have presented new results for immiscible compressible two-phase flow in porous media.More precisely, we give a week formulation and an existence result for a degenerate system modeling water-gas flow through a porous medium.The water is assumed to be incompressible and the gas phase is supposed compressible and obeying the ideal gas law leading to a new degeneracy in the evolution term of the pressure equation.Furthermore, a homogenization result for the corresponding system is established in the case of a single rock-type model.The extension to a porous medium made of several types of rocks, i.e. the porosity, the absolute permeability, the capillary and relative permeabilities curves are different in each type of porous media, is straightforward.Let us also mention that this homogenization result has been used successfully in [1] to simulate numerically a benchmark test proposed in the framework of the European Project FORGE: Fate Of Repository Gases [28].The study still needs to be improved in several areas such as the cases of unbounded capillary pressure and double porosity media.These more complicated cases appear in the applications.Further work on these important issues is in progress.