Investigating matrix-fracture transfer via a level set method for drainage and imbibition

Maša Prodanović, Steven L. Bryant, Zuleima Karpyn

Research output: Chapter in Book/Report/Conference proceedingConference contribution

6 Citations (Scopus)

Abstract

Multiphase flow and transport phenomena within fractures are important because fractures often represent primary flow conduits in otherwise low permeability rock. Flows within the fracture, between the fracture and the adjacent matrix, and through the pore space within the matrix typically happen on different length and time scales. Capturing these scales experimentally is difficult. It is therefore useful to have a computational tool that establishes the exact position and shape of fluid/fluid interfaces in realistic fracture geometries. The level set method is such a tool. Our progressive quasistatic (PQS) algorithm based on the level set method finds detailed, pore-level fluid configurations satisfying the Young-Laplace equation at a series of prescribed capillary pressures. The fluid volumes, contact areas and interface curvatures are readily extracted from the configurations. The method automatically handles topological changes of the fluid volumes as capillary pressure varies. It also accommodates arbitrarily complicated shapes of solid confining surfaces. Here we apply the PQS method to analytically defined fracture faces and aperture distributions, to geometries of fractures obtained from high-resolution images of real rocks, and to idealized fractures connected to a porous matrix. We also explicitly model a fracture filled with proppant, using a cooperative rearrangement algorithm to construct the proppant bed and the surrounding matrix. We focus on interface movement between matrix and fracture, and snap-off of non-wetting phase into the fracture during imbibition in particular. The configuration of fluids is strongly affected by asperities in unpropped fractures and by the locally open regions at the proppant/matrix interface. The area of phase in contact with the matrix is nonlinear with phase saturation and strongly hysteretic, and thus transfer functions based on saturations should be used with caution. The effect of coupling fracture capillarity and matrix capillarity on multiphase flow properties depends on the relative sizes of typical pore throats in the matrix and typical aperture in the fracture. The simulations agree with direct obsevations of fluid configurations in fractures.

Original languageEnglish (US)
Title of host publicationSPE Annual Technical Conference and Exhibition, ATCE 2008
Pages3136-3152
Number of pages17
Volume5
StatePublished - Dec 1 2008
EventSPE Annual Technical Conference and Exhibition, ATCE 2008 - Denver, CO, United States
Duration: Sep 21 2008Sep 24 2008

Other

OtherSPE Annual Technical Conference and Exhibition, ATCE 2008
CountryUnited States
CityDenver, CO
Period9/21/089/24/08

Fingerprint

Drainage
Capillarity
Fluids
Proppants
Multiphase flow
Rocks
Confined flow
Laplace equation
Geometry
Image resolution
Contacts (fluid mechanics)
Transfer functions

All Science Journal Classification (ASJC) codes

  • Fuel Technology
  • Energy Engineering and Power Technology

Cite this

Prodanović, M., Bryant, S. L., & Karpyn, Z. (2008). Investigating matrix-fracture transfer via a level set method for drainage and imbibition. In SPE Annual Technical Conference and Exhibition, ATCE 2008 (Vol. 5, pp. 3136-3152)
Prodanović, Maša ; Bryant, Steven L. ; Karpyn, Zuleima. / Investigating matrix-fracture transfer via a level set method for drainage and imbibition. SPE Annual Technical Conference and Exhibition, ATCE 2008. Vol. 5 2008. pp. 3136-3152
@inproceedings{ae2e9fc02bd3438cab6f6a56371fdd37,
title = "Investigating matrix-fracture transfer via a level set method for drainage and imbibition",
abstract = "Multiphase flow and transport phenomena within fractures are important because fractures often represent primary flow conduits in otherwise low permeability rock. Flows within the fracture, between the fracture and the adjacent matrix, and through the pore space within the matrix typically happen on different length and time scales. Capturing these scales experimentally is difficult. It is therefore useful to have a computational tool that establishes the exact position and shape of fluid/fluid interfaces in realistic fracture geometries. The level set method is such a tool. Our progressive quasistatic (PQS) algorithm based on the level set method finds detailed, pore-level fluid configurations satisfying the Young-Laplace equation at a series of prescribed capillary pressures. The fluid volumes, contact areas and interface curvatures are readily extracted from the configurations. The method automatically handles topological changes of the fluid volumes as capillary pressure varies. It also accommodates arbitrarily complicated shapes of solid confining surfaces. Here we apply the PQS method to analytically defined fracture faces and aperture distributions, to geometries of fractures obtained from high-resolution images of real rocks, and to idealized fractures connected to a porous matrix. We also explicitly model a fracture filled with proppant, using a cooperative rearrangement algorithm to construct the proppant bed and the surrounding matrix. We focus on interface movement between matrix and fracture, and snap-off of non-wetting phase into the fracture during imbibition in particular. The configuration of fluids is strongly affected by asperities in unpropped fractures and by the locally open regions at the proppant/matrix interface. The area of phase in contact with the matrix is nonlinear with phase saturation and strongly hysteretic, and thus transfer functions based on saturations should be used with caution. The effect of coupling fracture capillarity and matrix capillarity on multiphase flow properties depends on the relative sizes of typical pore throats in the matrix and typical aperture in the fracture. The simulations agree with direct obsevations of fluid configurations in fractures.",
author = "Maša Prodanović and Bryant, {Steven L.} and Zuleima Karpyn",
year = "2008",
month = "12",
day = "1",
language = "English (US)",
isbn = "9781605604824",
volume = "5",
pages = "3136--3152",
booktitle = "SPE Annual Technical Conference and Exhibition, ATCE 2008",

}

Prodanović, M, Bryant, SL & Karpyn, Z 2008, Investigating matrix-fracture transfer via a level set method for drainage and imbibition. in SPE Annual Technical Conference and Exhibition, ATCE 2008. vol. 5, pp. 3136-3152, SPE Annual Technical Conference and Exhibition, ATCE 2008, Denver, CO, United States, 9/21/08.

Investigating matrix-fracture transfer via a level set method for drainage and imbibition. / Prodanović, Maša; Bryant, Steven L.; Karpyn, Zuleima.

SPE Annual Technical Conference and Exhibition, ATCE 2008. Vol. 5 2008. p. 3136-3152.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

TY - GEN

T1 - Investigating matrix-fracture transfer via a level set method for drainage and imbibition

AU - Prodanović, Maša

AU - Bryant, Steven L.

AU - Karpyn, Zuleima

PY - 2008/12/1

Y1 - 2008/12/1

N2 - Multiphase flow and transport phenomena within fractures are important because fractures often represent primary flow conduits in otherwise low permeability rock. Flows within the fracture, between the fracture and the adjacent matrix, and through the pore space within the matrix typically happen on different length and time scales. Capturing these scales experimentally is difficult. It is therefore useful to have a computational tool that establishes the exact position and shape of fluid/fluid interfaces in realistic fracture geometries. The level set method is such a tool. Our progressive quasistatic (PQS) algorithm based on the level set method finds detailed, pore-level fluid configurations satisfying the Young-Laplace equation at a series of prescribed capillary pressures. The fluid volumes, contact areas and interface curvatures are readily extracted from the configurations. The method automatically handles topological changes of the fluid volumes as capillary pressure varies. It also accommodates arbitrarily complicated shapes of solid confining surfaces. Here we apply the PQS method to analytically defined fracture faces and aperture distributions, to geometries of fractures obtained from high-resolution images of real rocks, and to idealized fractures connected to a porous matrix. We also explicitly model a fracture filled with proppant, using a cooperative rearrangement algorithm to construct the proppant bed and the surrounding matrix. We focus on interface movement between matrix and fracture, and snap-off of non-wetting phase into the fracture during imbibition in particular. The configuration of fluids is strongly affected by asperities in unpropped fractures and by the locally open regions at the proppant/matrix interface. The area of phase in contact with the matrix is nonlinear with phase saturation and strongly hysteretic, and thus transfer functions based on saturations should be used with caution. The effect of coupling fracture capillarity and matrix capillarity on multiphase flow properties depends on the relative sizes of typical pore throats in the matrix and typical aperture in the fracture. The simulations agree with direct obsevations of fluid configurations in fractures.

AB - Multiphase flow and transport phenomena within fractures are important because fractures often represent primary flow conduits in otherwise low permeability rock. Flows within the fracture, between the fracture and the adjacent matrix, and through the pore space within the matrix typically happen on different length and time scales. Capturing these scales experimentally is difficult. It is therefore useful to have a computational tool that establishes the exact position and shape of fluid/fluid interfaces in realistic fracture geometries. The level set method is such a tool. Our progressive quasistatic (PQS) algorithm based on the level set method finds detailed, pore-level fluid configurations satisfying the Young-Laplace equation at a series of prescribed capillary pressures. The fluid volumes, contact areas and interface curvatures are readily extracted from the configurations. The method automatically handles topological changes of the fluid volumes as capillary pressure varies. It also accommodates arbitrarily complicated shapes of solid confining surfaces. Here we apply the PQS method to analytically defined fracture faces and aperture distributions, to geometries of fractures obtained from high-resolution images of real rocks, and to idealized fractures connected to a porous matrix. We also explicitly model a fracture filled with proppant, using a cooperative rearrangement algorithm to construct the proppant bed and the surrounding matrix. We focus on interface movement between matrix and fracture, and snap-off of non-wetting phase into the fracture during imbibition in particular. The configuration of fluids is strongly affected by asperities in unpropped fractures and by the locally open regions at the proppant/matrix interface. The area of phase in contact with the matrix is nonlinear with phase saturation and strongly hysteretic, and thus transfer functions based on saturations should be used with caution. The effect of coupling fracture capillarity and matrix capillarity on multiphase flow properties depends on the relative sizes of typical pore throats in the matrix and typical aperture in the fracture. The simulations agree with direct obsevations of fluid configurations in fractures.

UR - http://www.scopus.com/inward/record.url?scp=58849134307&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=58849134307&partnerID=8YFLogxK

M3 - Conference contribution

SN - 9781605604824

VL - 5

SP - 3136

EP - 3152

BT - SPE Annual Technical Conference and Exhibition, ATCE 2008

ER -

Prodanović M, Bryant SL, Karpyn Z. Investigating matrix-fracture transfer via a level set method for drainage and imbibition. In SPE Annual Technical Conference and Exhibition, ATCE 2008. Vol. 5. 2008. p. 3136-3152