Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks
Determination of the representative elementary volume (REV) of fractured rock masses based on equivalent permeability (K) is significantly dependent on the geometric characteristics of fractures. In this work, a series of numerical simulations were performed to analyze the relationship between geome...
Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Wiley
2022-01-01
|
| Series: | Geofluids |
| Online Access: | http://dx.doi.org/10.1155/2022/8073013 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1832562184926789632 |
|---|---|
| author | Ding Liu Hai Pu Shiru Guo Ziheng Sha Chong Li |
| author_facet | Ding Liu Hai Pu Shiru Guo Ziheng Sha Chong Li |
| author_sort | Ding Liu |
| collection | DOAJ |
| description | Determination of the representative elementary volume (REV) of fractured rock masses based on equivalent permeability (K) is significantly dependent on the geometric characteristics of fractures. In this work, a series of numerical simulations were performed to analyze the relationship between geometric characteristics of fractures and the REV size, in which fracture length follows a power-law distribution. A method to evaluate the K of a three-dimensional (3D) discrete fracture network (DFN) by extracting the equivalent pipe network (EPN) model from the DFN model was utilized and verified. The results show that K of the 3D DFN model has an exponential relationship with the power exponent (a) of fracture length distribution and the evaluation of K agrees well with that reported in previous studies, confirming the reliability of the EPN model for calculating seepage properties of complex 3D DFN models. When the side length of submodels (Ln) is small, the K varies significantly due to the influence of random number seeds used to generate fracture length, location, and orientation. The K holds a constant value after Ln exceeds some specific value. The critical model scale is determined as the REV size, and the corresponding volume of the 3D DFN model is represented by VREV. The VREV varies within a narrow range when a≤4.0. When a=4.5, the VREV rapidly increases to more than 3.4 times than that when a=4.0. The fluid flow becomes more inhomogeneous due to the small nonpersistent fractures that dominate the preferential flow paths when a exceeds a certain value (i.e., 4.5). The K at the REV size decreases exponentially with the increment of a. This tendency can be explained by the decrease of the average intersection length (Li) with the increment of a, which is a geometric parameter for reflecting the connectivity of the fracture network. |
| format | Article |
| id | doaj-art-5bd7182f33ca4a72af94bba3191eacbd |
| institution | Kabale University |
| issn | 1468-8123 |
| language | English |
| publishDate | 2022-01-01 |
| publisher | Wiley |
| record_format | Article |
| series | Geofluids |
| spelling | doaj-art-5bd7182f33ca4a72af94bba3191eacbd2025-02-03T01:23:15ZengWileyGeofluids1468-81232022-01-01202210.1155/2022/8073013Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture NetworksDing Liu0Hai Pu1Shiru Guo2Ziheng Sha3Chong Li4State Key Laboratory for Geomechanics and Deep Underground EngineeringState Key Laboratory for Geomechanics and Deep Underground EngineeringState Key Laboratory for Geomechanics and Deep Underground EngineeringState Key Laboratory for Geomechanics and Deep Underground EngineeringSchool of MinesDetermination of the representative elementary volume (REV) of fractured rock masses based on equivalent permeability (K) is significantly dependent on the geometric characteristics of fractures. In this work, a series of numerical simulations were performed to analyze the relationship between geometric characteristics of fractures and the REV size, in which fracture length follows a power-law distribution. A method to evaluate the K of a three-dimensional (3D) discrete fracture network (DFN) by extracting the equivalent pipe network (EPN) model from the DFN model was utilized and verified. The results show that K of the 3D DFN model has an exponential relationship with the power exponent (a) of fracture length distribution and the evaluation of K agrees well with that reported in previous studies, confirming the reliability of the EPN model for calculating seepage properties of complex 3D DFN models. When the side length of submodels (Ln) is small, the K varies significantly due to the influence of random number seeds used to generate fracture length, location, and orientation. The K holds a constant value after Ln exceeds some specific value. The critical model scale is determined as the REV size, and the corresponding volume of the 3D DFN model is represented by VREV. The VREV varies within a narrow range when a≤4.0. When a=4.5, the VREV rapidly increases to more than 3.4 times than that when a=4.0. The fluid flow becomes more inhomogeneous due to the small nonpersistent fractures that dominate the preferential flow paths when a exceeds a certain value (i.e., 4.5). The K at the REV size decreases exponentially with the increment of a. This tendency can be explained by the decrease of the average intersection length (Li) with the increment of a, which is a geometric parameter for reflecting the connectivity of the fracture network.http://dx.doi.org/10.1155/2022/8073013 |
| spellingShingle | Ding Liu Hai Pu Shiru Guo Ziheng Sha Chong Li Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks Geofluids |
| title | Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks |
| title_full | Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks |
| title_fullStr | Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks |
| title_full_unstemmed | Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks |
| title_short | Numerical Investigations on the Effect of Fracture Length Distribution on the Representative Elementary Volume of 3D Discrete Fracture Networks |
| title_sort | numerical investigations on the effect of fracture length distribution on the representative elementary volume of 3d discrete fracture networks |
| url | http://dx.doi.org/10.1155/2022/8073013 |
| work_keys_str_mv | AT dingliu numericalinvestigationsontheeffectoffracturelengthdistributionontherepresentativeelementaryvolumeof3ddiscretefracturenetworks AT haipu numericalinvestigationsontheeffectoffracturelengthdistributionontherepresentativeelementaryvolumeof3ddiscretefracturenetworks AT shiruguo numericalinvestigationsontheeffectoffracturelengthdistributionontherepresentativeelementaryvolumeof3ddiscretefracturenetworks AT zihengsha numericalinvestigationsontheeffectoffracturelengthdistributionontherepresentativeelementaryvolumeof3ddiscretefracturenetworks AT chongli numericalinvestigationsontheeffectoffracturelengthdistributionontherepresentativeelementaryvolumeof3ddiscretefracturenetworks |