A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach
Modeling flow and transport in fractured porous media has been a topic of intensive research for a number of energy- and environment-related industries. The presence of multiscale fractures makes it an extremely challenging task to resolve accurately and efficiently the flow dynamics at both the loc...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2020-12-01
|
| Series: | Energies |
| Subjects: | |
| Online Access: | https://www.mdpi.com/1996-1073/13/24/6667 |
| _version_ | 1797544423390707712 |
|---|---|
| author | Wenjuan Zhang Waleed Diab Hadi Hajibeygi Mohammed Al Kobaisi |
| author_facet | Wenjuan Zhang Waleed Diab Hadi Hajibeygi Mohammed Al Kobaisi |
| author_sort | Wenjuan Zhang |
| collection | DOAJ |
| description | Modeling flow and transport in fractured porous media has been a topic of intensive research for a number of energy- and environment-related industries. The presence of multiscale fractures makes it an extremely challenging task to resolve accurately and efficiently the flow dynamics at both the local and global scales. To tackle this challenge, we developed a computational workflow that adopts a two-level hierarchical strategy based on fracture length partitioning. This was achieved by specifying a partition length to split the discrete fracture network (DFN) into small-scale fractures and large-scale fractures. Flow-based numerical upscaling was then employed to homogenize the small-scale fractures and the porous matrix into an equivalent/effective single medium, whereas the large-scale fractures were modeled explicitly. As the effective medium properties can be fully tensorial, the developed hierarchical framework constructed the discrete systems for the explicit fracture–matrix sub-domains using the nonlinear two-point flux approximation (NTPFA) scheme. This led to a significant reduction of grid orientation effects, thus developing a robust, applicable, and field-relevant framework. To assess the efficacy of the proposed hierarchical workflow, several numerical simulations were carried out to systematically analyze the effects of the homogenized explicit cutoff length scale, as well as the fracture length and orientation distributions. The effect of different boundary conditions, namely, the constant pressure drop boundary condition and the linear pressure boundary condition, for the numerical upscaling on the accuracy of the workflow was investigated. The results show that when the partition length is much larger than the characteristic length of the grid block, and when the DFN has a predominant orientation that is often the case in practical simulations, the workflow employing linear pressure boundary conditions for numerical upscaling give closer results to the full-model reference solutions. Our findings shed new light on the development of meaningful computational frameworks for highly fractured, heterogeneous geological media where fractures are present at multiple scales. |
| first_indexed | 2024-03-10T14:00:18Z |
| format | Article |
| id | doaj.art-a54a3d35b3cc40e5b64520abedbe10ba |
| institution | Directory Open Access Journal |
| issn | 1996-1073 |
| language | English |
| last_indexed | 2024-03-10T14:00:18Z |
| publishDate | 2020-12-01 |
| publisher | MDPI AG |
| record_format | Article |
| series | Energies |
| spelling | doaj.art-a54a3d35b3cc40e5b64520abedbe10ba2023-11-21T01:15:48ZengMDPI AGEnergies1996-10732020-12-011324666710.3390/en13246667A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling ApproachWenjuan Zhang0Waleed Diab1Hadi Hajibeygi2Mohammed Al Kobaisi3Department of Petroleum Engineering, Khalifa University of Science and Technology, Abu Dhabi 127788, UAEDepartment of Petroleum Engineering, Khalifa University of Science and Technology, Abu Dhabi 127788, UAEDepartment of Geoscience and Engineering, Delft University of Technology, 5048 Delft, The NetherlandsDepartment of Petroleum Engineering, Khalifa University of Science and Technology, Abu Dhabi 127788, UAEModeling flow and transport in fractured porous media has been a topic of intensive research for a number of energy- and environment-related industries. The presence of multiscale fractures makes it an extremely challenging task to resolve accurately and efficiently the flow dynamics at both the local and global scales. To tackle this challenge, we developed a computational workflow that adopts a two-level hierarchical strategy based on fracture length partitioning. This was achieved by specifying a partition length to split the discrete fracture network (DFN) into small-scale fractures and large-scale fractures. Flow-based numerical upscaling was then employed to homogenize the small-scale fractures and the porous matrix into an equivalent/effective single medium, whereas the large-scale fractures were modeled explicitly. As the effective medium properties can be fully tensorial, the developed hierarchical framework constructed the discrete systems for the explicit fracture–matrix sub-domains using the nonlinear two-point flux approximation (NTPFA) scheme. This led to a significant reduction of grid orientation effects, thus developing a robust, applicable, and field-relevant framework. To assess the efficacy of the proposed hierarchical workflow, several numerical simulations were carried out to systematically analyze the effects of the homogenized explicit cutoff length scale, as well as the fracture length and orientation distributions. The effect of different boundary conditions, namely, the constant pressure drop boundary condition and the linear pressure boundary condition, for the numerical upscaling on the accuracy of the workflow was investigated. The results show that when the partition length is much larger than the characteristic length of the grid block, and when the DFN has a predominant orientation that is often the case in practical simulations, the workflow employing linear pressure boundary conditions for numerical upscaling give closer results to the full-model reference solutions. Our findings shed new light on the development of meaningful computational frameworks for highly fractured, heterogeneous geological media where fractures are present at multiple scales.https://www.mdpi.com/1996-1073/13/24/6667simulationfracturesdiscrete fracture network (DFN)embedded discrete fracture–matrix (EDFM)upscalingnonlinear two-point flux approximation (NTPFA) |
| spellingShingle | Wenjuan Zhang Waleed Diab Hadi Hajibeygi Mohammed Al Kobaisi A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach Energies simulation fractures discrete fracture network (DFN) embedded discrete fracture–matrix (EDFM) upscaling nonlinear two-point flux approximation (NTPFA) |
| title | A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach |
| title_full | A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach |
| title_fullStr | A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach |
| title_full_unstemmed | A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach |
| title_short | A Computational Workflow for Flow and Transport in Fractured Porous Media Based on a Hierarchical Nonlinear Discrete Fracture Modeling Approach |
| title_sort | computational workflow for flow and transport in fractured porous media based on a hierarchical nonlinear discrete fracture modeling approach |
| topic | simulation fractures discrete fracture network (DFN) embedded discrete fracture–matrix (EDFM) upscaling nonlinear two-point flux approximation (NTPFA) |
| url | https://www.mdpi.com/1996-1073/13/24/6667 |
| work_keys_str_mv | AT wenjuanzhang acomputationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT waleeddiab acomputationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT hadihajibeygi acomputationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT mohammedalkobaisi acomputationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT wenjuanzhang computationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT waleeddiab computationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT hadihajibeygi computationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach AT mohammedalkobaisi computationalworkflowforflowandtransportinfracturedporousmediabasedonahierarchicalnonlineardiscretefracturemodelingapproach |