A Comprehensive Numerical Model for Simulating Fluid Transport in Nanopores
Yuan Zhang^{1}, Wei Yu^{2}, Kamy Sepehrnoori^{3 }& Yuan Di^{1}
^{1}College of Engineering, Peking University, Beijing, 100871, China. ^{2}Department of Petroleum Engineering, Texas A & M University, Collage Station, TX, 77843, USA. ^{3}Department of Petroleum and Geosystems Engineering, University of Texas at Austin, Austin, TX, 78712, USA.
Abstract
Since a large amount of nanopores exist in tight oil reservoirs, fluid transport in nanopores is complex due to large capillary pressure. Recent studies only focus on the effect of nanopore confinement on singlewell performance with simple planar fractures in tight oil reservoirs. Its impacts on multiwell performance with complex fracture geometries have not been reported. In this study, a numerical model was developed to investigate the effect of confined phase behavior on cumulative oil and gas production of four horizontal wells with different fracture geometries. Its pore sizes were divided into five regions based on nanopore size distribution. Then, fluid properties were evaluated under different levels of capillary pressure using PengRobinson equation of state. Afterwards, an efficient approach of Embedded Discrete Fracture Model (EDFM) was applied to explicitly model hydraulic and natural fractures in the reservoirs. Finally, three fracture geometries, i.e. nonplanar hydraulic fractures, nonplanar hydraulic fractures with one set natural fractures, and nonplanar hydraulic fractures with two sets natural fractures, are evaluated.
The multiwell performance with confined phase behavior is analyzed with permeabilities of 0.01 md and 0.1 md. This work improves the analysis of capillarity effect on multiwell performance with complex fracture geometries in tight oil reservoirs.
Introduction
Tight reservoirs make significant contributions to the oil and gas productions in recent years. These reservoirs are unique because of the low permeability, which requires multiple horizontal wells and multistage hydraulic fractures for economic production^{1}. Permeability is the intrinsic property of rocks, which measures the ability of rocks to allow liquids or gases to flow through it. Permeability is typically expressed in units of millidarcy (md), and a large value indicates the liquids or gases are highly permeable through rocks. Since the natural fractures are generally distributed in the formation, the combination of horizontal drilling and hydraulic fracturing is an efficient approach to generate a large contact area between the wellbore and low permeability formation. The actual fracturing process often creates complex fracture geometries in tight reservoirs due to the preexisting natural fractures^{2}^{,3}.
High capillary pressure due to nanopores in tight reservoirs impacts the fluid phase behavior and longterm estimated ultimate recovery (EUR) prediction in unconventional reservoirs^{4}^{,5}. Theoretical work has shown that the bubblepoint pressure decreases due to the reduction of pore radius^{6}^{,7}^{,8}^{,9}. The further the reservoir condition is from the critical point, the larger the reduction of bubblepoint pressure. Since traditional PressureVolumeTemperature (PVT) analysis failed to calculate the phase behavior with the capillary pressure effect, recent studies have focused on developing new algorithms to evaluate its effect on fluild properties and well performance in unconventional reservoirs. Du and Chu^{10} used a commercial reservoir simulator to investigate the influence of PVT variations on singlewell performance in the Bakken formation. They discussed the confined PVT properties in two permeability scenarios, however, the pore size distribution and fracture networks were not included in their work. Alternately, Wang et al.^{11} incorporated the capillary pressure effect and pore space compaction in a compositional reservoir simulator to evaluate their effects on production. They considered the change of porosity and permeability with effective stress, but no discrete fractures were considered in their simulations. On the other hand, Nojabaei et al.^{12} used a compositionallyextended black oil model to incorporate the capillarity effect on phase behavior.
The pore size distribution was considered in their model, but the hydraulic fractures were not included for simulation of singlewell performance. Rezaveisi et al.^{13} implemented the capillarity equilibrium in an inhouse compositional reservoir simulator to forecast the well performance. They compared oil recoveries of tight oil and gas condensate cases, but complex fractures geometries were not investigated. Furthermore, Melero et al.^{14} investigated the effect of fracture complexity on the well performance of CO_{2}EOR in tight oil reservoirs.
However, the capillary pressure effect was not taken into account. While Wu and Zhang^{15} discussed the gas transport in nanopores by molecular dynamic simulation; Wang et al.^{16} used the molecular dynamics model to reproduce the experimental contact angle and calibrate the parameters with nanopore confinement, but neither considered the fractures in the formations. Additionally, Siripatrachai et al.^{17} developed a compositional reservoir simulator incorporating the capillary pressure effect on phase behavior with several discrete hydraulic fractures. However, the effect of complex natural fractures in simulation of multiwell performance was not examined.
The objective of this study is to investigate the capillary pressure effect on performance of multiple wells with different fracture geometries. Newton fluid is considered and the governing equations for the oil and gas transport are applied in this work. Embedded Discrete Fracture Model (EDFM) is utilized to efficiently handle the complex fractures. First of all, EDFM was verified to approach for biwing hydraulic fractures against Local Grid Refinement (LGR) method, which can accurately simulate the oil and gas flow from shale matrix to simple fractures^{18}. Subsequently, the verified methodology was used to build a numerical reservoir model including four horizontal wells to simulate the capillary pressure effect on well performance with three complex fracture geometries. Two permeabilities of 0.01 md and 0.1 md are taken into account. This work improves the analysis of multiwell performance with complex fracture geometries in tight oil reservoirs.
Results
Base case
The Bakken formation is located in the Williston Basin. It consists of Upper and Lower Bakken Shales, Middle Bakken, and Three Forks. The Middle Bakken region, which is the main pay zone, has an estimated average oil resource of 3.65 billion barrels^{19}^{,20}. In the development of Middle Bakken formation, four horizontal wells in a reservoir size of 1 mile × 2 miles are often drilled with multistage hydraulic fractures, which permit more connected fracture networks for production^{21}.
Therefore, a base reservoir model is set up, including four horizontal wells with 30 planar hydraulic fractures for each well (Supplementary Figure 1), which is more reliable than modeling a single well. The reservoir size is 3,219 m × 1,609 m × 15 m, corresponding to length, width, and thickness, respectively. It consists of 264 grids in x direction, 132 grids in y direction, and 1 grid in z direction. The horizontal wells with lateral length of 2,829 m are also incorporated in the model. The fractures in the neighboring wells are in a zipper pattern, with fracture spacing of 97 m and fracture halflength of 79 m. The initial reservoir pressure is 51.7 MPa, and initial water saturation is 49%. The compositions and critical properties has been described previously^{19}.
Figure 1. Verification of EDFM approach. (a) Flowing bottomhole pressure for simulation of Middle Bakken tight oil reservoirs; (b) Nanopore size distribution in the reservoir model. (Color bar of 1–5 respresents five different PVT regions: less than 10 nm, 10–20 nm, 20–30 nm, 30–50 nm, and larger than 50 nm); (c) Comparison of cumulative oil production for a 30year period without the capillary pressure effect using EDFM and LGR approach; (d) Comparison of cumulative oil production for a 30year period without the capillary pressure effect using EDFM and LGR approach.
Matrix permeability is 0.071 md, and matrix porosity is 5.6%. Total compressibility is 1.5 × 10^{−15} MPa^{−1} and fracture conductivity is 15 mdm. The reservoir temperature is 116 °C. The wateroil and liquidgas relative permeability curves are from the previous study^{22}, which were obtained from history matching from a production well in the Middle Bakken formation. It is noted that the oil and gas properties, rock properties, and fracture properties are within the reasonable range of typical oil properties reported in the literature^{19}^{,23}^{,24}^{,25}. In real field production, the bottomhole pressure does not decrease rapidly, but gradually from a larger value to a smaller value within a certain period. In order to mimic the actual field observation, bottomhole pressure was set to decrease from 48.3 to 3.5 MPa around one and a half years, and remains 3.5 MPa afterwards until the end of production, as shown in Fig. 1(a). The low and high values of flowing bottomhole pressure are given based on the field observation.
Verification of Embedded Discrete Fracture Model (EDFM)
The performance of base case was evaluated to verify the EDFM approach against LGR approach with and without the capillary pressure effect. As the aforementioned discussion, EDFM approach is easier for handling the complex fracture networks and its computational speed is several times faster than the LGR approach when modeling complex fracture geometries^{26}. In order to consider the capillary pressure effect, the raw data of Mercury Injection Capillary Pressure (MICP) test and the nanopore size distribution were used for the Middle Bakken formation^{23}^{,24}^{,26}^{,27},
and the nanopore sizes were divided into five different regions: less than 10 nm, 10–20 nm, 20–30 nm, 30–50 nm, and bulk (pore sizes larger than 50 nm). It should be mentioned that different rocks need different patterns and the division mainly depends on the actual pore size distribution of the tight rock. The sensitivity studies were then performed to identify the computation time and accuracy of different patterns. Finally, distribution into five regions is found to better maintain the balance of accuracy and computation speed. Figure 1(b) displays the distribution of five different regions in the base reservoir model when the capillary pressure effect is considered.
The comparison of cumulative oil production of four horizontal wells for 30 years with and without the capillary pressure effect using LGR and EDFM approach is presented in Fig. 1(c) and (d). The cumulative oil or gas production is defined as the total amount of oil and gas production from an oil reservoir of a particular time. It is based on the mass balance equation and can be calculated by multiplying the amount of production by rate. The calculation of cumulative production can be found in the supplementary information. As shown, a good agreement is obtained, indicating that the EDFM approach can well simulate the well performance with hydraulic fractures in tight oil reservoirs.
It should be noted that in the real case studies, there are complex nonplanar hydraulic fractures with natural fractures. The LGR method is challenging to deal with such complex fracture geometries. The EDFM approach can be applied to effectively capture the complex fracture geometries and simulate their effects on oil and gas transport phenomenon in fractured porous media. Although the simulation production results were not compared with the actual field production data, which are difficult to obtain now, the simulations and prediction results are reasonable, for that the reservoir, fracture, and fluid properties used in the case studies are based on the experimental data in the Bakken tight oil reservoir^{22}^{,23}^{,24}^{,25}. In the future, great efforts will be made to further investigate the comparisons with the actual field production data.
Effect of different fracture geometries on well performance
Investigation of the effect of different fracture geometries on well performance were performed in this section. Since the nonplanar hydraulic fractured cases with natural fractures represent the realistic fracture networks in tight reservoirs, the simulation tests for four different fracture geometries for the following cases:

Case 1: Nonplanar hydraulic fractures.

Case 2: Nonplanar hydraulic fractures with 500 natural fractures (NF) in one set.

Case 3: Nonplanar hydraulic fractures with 1,000 natural fractures (NF) in one set.

Case 4: Nonplanar hydraulic fractures with 1,000 natural fractures (NF) in two sets.
Figure 2. Fracture distribution for four different models considering hydraulic and natural fractures. (a) Nonplanar hydraulic fractures; (b) Nonplanar hydraulic fractures with 500 NF in one set; (c) Nonplanar hydraulic fractures with 1,000 NF in one set; (d) Nonplanar hydraulic fractures with 1,000 NF in two sets.
Fracture geometries for four cases are depicted in Fig. 2. The fracture halflength varies from 45.7 to 137.2 m. Natural fractures is normally distributed in the models based on the assumption that their orientations are parallel to the horizontal wellbore. The angles along the x direction for one set of natural fractures are between 0–10 deg and 90–100 deg for another set. The length of natural fractures ranges from 30.5 to 91.4 m and the nature fracture conductivity is 1.5 mdm. The hydraulic fracture conductivity is 15 mdm.
Figure 3. Well performance for a 30year period for different fracture geometries. (a) Cumulative oil production; (b) Cumulative gas production.
The simulation tests for four cases investigate the impacts of different fracture geometries and the contribution of the natural fractures on well performance. The simulation results for four cases are presented in Fig. 3.
As illustrated in Fig. 3, the cumulative oil and gas production increase as the fracture geometries become complex. When there are 500 natural fractures in the model, the cumulative oil and gas production increase by 1.5% and 6.8%, respectively, compared to the nonplanar case. As 1,000 natural fractures exist in one set in the model, the cumulative oil and gas production is 2.8% and 11.4% higher, respectively. If 1,000 natural fractures set in two sets, the cumulative oil and gas production can be as high as 5.0% and 14.0%, respectively, reflecting a significant impact.
That is, complex fracture geometries have a positive influence on the well performance in tight reservoirs. When the number of natural fractures in one set increases from 500 to 1,000, more natural fractures connect to hydraulic fractures, leading to the increase in the cumulative production. When 1,000 natural fractures change from one set to two sets, it can be found that there are more natural fractures connecting to hydraulic fractures and also more natural fractures connecting each other. Figure 4 depicts the difference of pressure distribution. As shown, the more complex the fracture geometries are, the more overlaps between the neighboring wells. The pressure of Case 4 is closer to the minimum value (color is closer to purple). The difference of average pressure of four cases is plotted in Fig. 4(e) to better present the pressure drawdown in a 30year period. It can be seen that the pressure decreases faster when the fracture geometries become complex, especially at the beginning of production time (5 years).
Figure 4. Pressure distribution for four cases at the end of 30 years. (a) Nonplanar hydraulic fractures; (b) Nonplanar hydraulic fractures with 500 NF; (c) Nonplanar hydraulic fractures with 1,000 NF in one set; (d) Nonplanar hydraulic fractures with 1,000 NF in two sets; (e) Difference of average pressure of four cases.
In addition, the difference of pressure distribution for four cases in Fig. 4 shows that the drainage area becomes larger as the fracture geometries become complex. Hence, the increase in the cumulative production is more obvious. The fluid flows more easily from the formation to the wellbore due to the more connected networks.
The oil saturation of the four cases in the production period is presented in Supplementary Figure 2. The noticeable change in the oil saturation indicates that the presence of natural fractures provides more connected fracture networks. The high conductivity allows the oil and gas to transport more easily from matrix to the horizontal wells with higher transmissibility and improves the well performance. Hence, the natural fractures should be included in the simulation model in tight formations, especially with high density of natural fractures.
Effect of nanopore confinement on well performance for different fracture geometries
In this section, we evaluated the capillary pressure effect with the aforementioned fracture geometries. The cumulative oil and gas production for a 30year period with and without the capillary pressure effect are summarized in Figs 5 and and 6, respectively.
Figure 5. Cumulative oil production for a 30year period with and without the capillary pressure effect. (a) Nonplanar hydraulic fractures; (b) Nonplanar hydraulic fractures with 500 NF in one set; (c) Nonplanar hydraulic fractures with 1,000 NF in one set; (d) Nonplanar hydraulic fractures with 1,000 NF in two sets.
The cumulative oil productions in Fig. 5 increase when the capillary pressure effect is taken into account. As a result of the bubblepoint pressure decreases, the singlephase production period becomes longer and it improves the oil production. In comparison to the oil production in a period of 30 years for each fracture geometry, the cumulative oil productions with the capillary pressure effect are 15%, 17.5%, 19%, and 20% higher than the cases without the capillary pressure effect.
Figure 6. Cumulative gas production for a 30year period with and without the capillary pressure effect. (a) Nonplanar hydraulic fractures; (b) Nonplanar hydraulic fractures with 500 NF in one set; (c) Nonplanar hydraulic fractures with 1,000 NF in one set; (d) Nonplanar hydraulic fractures with 1,000 NF in two sets.
The difference of pressure distribution clearly describes the capillary pressure effect. Only nonplanar fractures case is presented for neatness in Fig. 7. As illustrated, pressure drops rapidly and drainage area becomes larger when the capillary pressure effect is taken into account, leading to the increase in the cumulative oil production at the end of 30 years. The capillary pressure effect for cumulative gas production is more complex than oil (Fig. 6). In the early period, the gas production of the case with the capillary pressure effect is higher compared to the case without the capillary pressure effect.
The reason being that when considering the capillary pressure effect, more oil is produced with solution gas.
Figure 7. Pressure distribution of the nonplanar fractures case at the end of 30year production with and without the capillary pressure effect. (a) Without the capillary pressure effect; (b) With the capillary pressure effect.
The solution gas is released and dominates gas production at the beginning. As the reservoir pressure depletes, it will reach the bubblepoint pressure; the fluid goes from singlephase to twophase region and gas production continues to increase. In this period, solution gas does not dominate the gas production any longer. It can be observed that the bubblepoint pressure reduces by 1.4 MPa at the reservoir temperature of 116 °C due to the capillary pressure effect (Supplementary Figure 3). The reduction caused a shorter twophase period for the case with the capillary pressure effect, leading to less gas production in the later period. Additionally, the intersection time between the cases with and without the capillary pressure effect will move to a later time when the fracture geometries become more complex.
Effect of matrix permeability on well performance
Matrix permeability is another key parameter affecting the well performance in tight oil reservoirs. In this section, the permeability value of 0.1 md is selected to evaluate the well performance for different fracture geometries. The capillary pressure effect is taken into account and the results were compared with the case of 0.01 md.
Figure 8. Well performance for a 30year period with the permeability of 0.1 md. (a) Cumulative oil production for a 30year period for different fracture geometries with the permeability of 0.1 md; (b) Cumulative gas production for a 30year period for different fracture geometries with the permeability of 0.1 md; (c) Cumulative oil production of the nonplanar fractures case with the permeability of 0.01 md (blue curves) and 0.1 md (red curves); (d) Cumulative gas production of the nonplanar fractures case with the permeability of 0.01 md (blue curves) and 0.1 md (red curves).
Figure 8(a) and (b) present the cumulative oil and gas production with different fracture geometries under the permeability of 0.1 md. Compared to the nonplanar hydraulic fractures case, the cumulative oil production increases 0.8% for the nonplanar hydraulic fractures with 500 NF case, 1.5% for the nonplanar hydraulic fractures with 1,000 NF in one set case, and 2.1% for the nonplanar hydraulic fractures with 1,000 NF in two sets case. The cumulative gas production increases 1.3%, 3.0%, and 4.3%, respectively. The incremental oil and gas production reflect a smaller impact compared with the results in Fig. 3. That is, the presence of natural fractures plays a less significant role in the case of high matrix permeability.
The well performance of the nonplanar fractures case is analyzed in detail, as shown in Fig. 8(c) and (d). In comparison with the results with matrix permeability of 0.01 md, the cumulative oil production increases by 24% for high permeability (blue and red solid lines). The larger incremental oil production is observed in Fig. 8(c) when the capillary pressure effect is considered (dash line and solid line), indicating that the capillary pressure effect is more significant for high permeability. For cumulative gas production in Fig. 8(d), the similar results can be observed.
The differences of the average pressure and oil saturation between low and high permeability cases without the capillary pressure are presented in Supplementary Figures 4 and 5. The pressure drops 3.1 MPa higher and oil saturation decreases 10% higher for the high permeability case than the low permeability case. The large difference indicates the increase in the production of high permeability case (0.01 md).
Although high permeability represents larger pore sizes, which is assigned to small capillary pressure, a small difference with and without capillary pressure for the high permeability is expected. However, in this study, it is assumed that the capillary pressure is independent of pore sizes and is constant for both permeability values. Consequently, a high difference with and without capillary for the high permeability case is seen. The relationship between permeability and pore sizes is beyond the scope of this work, which will be investigated in future works.
Discussion
In this study, the effects of various fracture geometries and capillary pressure on multiwell performance are investigated. The EDFM approach is applied to describe the effect of fluid flow in the complex fracture geometries. The presence of natural fractures contributes to the multiwell performance; the incremental oil and gas production of the nonplanar hydraulic fractures case with 1,000 NF in two sets is 5% and 14% higher compared to the nonplanar hydraulic fractures case.
Based on the well performance for different fracture geometries with and without the capillary pressure effect, the cumulative oil production increases when considering the capillary pressure effect.
However, the incremental oil production is similar for various fracture geometries. The gas production, in the early period, is higher compared to the case without the capillary pressure effect due to the solution gas; however, the gas production of the case without capillary pressure improves because of longer twophase period. Additionally, the intersection time between cases with and without the capillary pressure effect moves to a later time for complex fracture geometries.
For the cases with high permeability of 0.1 md, the incremental oil and gas productions for the nonplanar hydraulic fractures with 1,000 NF in two sets are 2.1% and 4.3%, respectively, which indicates the contribution of natural fractures is less significant for high permeability. Therefore, the natural fractures should be included in the numerical model of tight formations, especially for low permeability and high density of natural fracture networks.
Methodology
Blackoil fluid properties affected by nanopore confinement
Flash calculation accounting for the inequalities of liquid and vapor phases is applied to calculate the fluid properties, i.e. oil and gas viscosities, densities, solution gasoil ratio as a function of oil pressure. In this work, the oil pressure is assumed as the refrence pressure. Hence, the gas pressure is calculated using the YoungLaplace equation^{28} given below:
where P_{V} and P_{L} refer to the pressures in the vapor and liquid phase, respectively. r is the radius of capillary tube in porous media; θ is the contact angle representing the surface tension between solid and fluids; σ is the interfacial tension at the vaporliquid phase interface.
Phase equilibrium is achieved when the fugacities of components in the vapor and liquid phases are equal^{29}. The criterion is given as
where f^{i}_{L }and f^{i}_{V }represent the fugacity of component i in the liquid and vapor phases. T is the reservoir temperature; N_{c }is the number of fluid components. x_{i} and y_{i }are liquid and vapor compositions, respectively. PengRobinson equation of state and RachfordRice flash calculation are modified considering the difference between the vapor and liquid pressures. The mass balance equations are applied as well, which are given in Equations (3) through (5).
where z_{i} is the overall mole fraction of component i. F is the number of moles of original feed. L and V are the number of moles of liquid and vapor phases, respectively. K_{c}^{i} is the Kvalues considering the inequalities between liquid and vapor pressures. More details about the phase equilibrium calculation can be found in the work by Zhang et al.^{21}.
Based on phase equilibrium calculation, the blackoil fluid properties such as viscosities, solution gasoil ratio (R_{s}), and formation volume factor (B_{o}) is determined in the Middle Bakken formation for different pore sizes. In order to obtain the solution gasoil ratio, R_{s}, the equilibrium liquid composition is flashed at the standard conditions. The equation is described as
where n_{V} and n_{L} are the vapor and liquid phase molar fraction, respectively. V_{mV} and V_{mL} are the vapor and liquid molar volume, respectively. SC represents the standard condition.
The oil formation volume factor B_{o} is defined as
where RC represents the reservoir condition. The oil viscosity is calculated by^{30}
where μ is the oil viscosity at the reservior condition; a_{0} through a_{4} are 0.1023, 0.023364, 0.058533, −0.040758 and 0.0093324, respectively. ρ_{r} is the reduced density, which is defined as
where v_{ci} is identical to the critical volume of component i; the value of parameter α is equal to 1. The mixture viscosity parameter ξ is evaluated as
where M_{i} is the molecular weight of component i. The low pressure viscosities of the mixture in Equation (8) is defined as
where µ_{i}* , the low pressure viscosities for pure substances are calculated using Stiel and Thodos equation^{31}:
with ξ_{i }= Tci^{1/6} M_{i}^{1/2} Pci^{2/3} and T_{ri }=T/T_{ci}.
The calculation of the oil viscosity also can be applied in the gas viscosity. We note that the properties calculation discussed above are not limited to the fluid properties considering the capillary equilibrium. Once the capillary pressure is set to zero, it can then be applied to the bulk fluid. The tuned PVT data is input into a standard blackoil simulator to predict the well performance in the Middle Bakken formation.
Embedded Discrete Fracture Model (EDFM)
The complex fracture geometries in tight oil reservoirs are common because of the interaction of hydraulic fractures and the preexisting natural fractures. It is important to accurately model the influence of complex fractures on well performance in reservoir simulation. Traditional LGR approach can accurately model the fluid transport from matrix to fractures, but it is difficult to deal with complex nonplanar fractures^{32}. The EDFM is an efficient approach to handle the complex fracture geometries through discretizing the fractures into segments with matrix cell boundaries^{33}^{,34}^{,35}. In addition, virtual cells are added for these fracture segments. The NonNeighboring Connections (NNCs) are used for these cells to account for fluid transport associated with fractures, including the flow between matrix and fractures, flow inside an individual fracture, and flow between intersecting fractures. The volume flow rate of phase l between cells in a NNC pair is
where λ_{l} is the relative mobility of phase l, ΔP is the potential difference between the cells, and T_{NNC} is the NNC transmissibility factor, which can be calculated by
For the matrixfracture connection, k_{NNC} is the matrix permeability in the direction perpendicular to the fracture plane, A_{NNC} is the area of the fracture plane inside the matrix block, and d_{NNC} is the average normal distance from matrix block to fracture plane. For the connections between fracture segments, k_{NNC} is an average of fracture permeability, A_{NNC} is the common area between fracture segments, and d_{NNC} is the distance between centroids of the fracture segments. For the fracturewellbore intersections, an effective well index is used and calculated by
where k_{f} is the fracture permeability, w_{f} is the fracture aperture, and l and w are the fracture segment length and height, respectively. More details about the calculations of NNC transmissibility factors can be found in the work by Xu et al.^{26}. In addition, the EDFM approach can be applied in traditional reservoir simulators in a nonintrusive manner^{26}.
Acknowledgments
This work was supported by National Natural Science Foundation of China (51674010), National Science and Technology Major Project of China (2016ZX05014), and China Scholarship Council (CSC File No. 201506010205). We would like to acknowledge Computer Modeling Group Ltd. for providing the CMG software for this study. We greatly appreciate Dr. Hamid R. Lashgari for his discussion and comments on our research.
Author Contributions
Yuan Zhang and Wei Yu developed the model and conducted numerical reservoir simulation. Yuan Zhang and Wei Yu wrote the main manuscript and revised the manuscript. Yuan Di and Kamy Sepehrnoori supervised the project. All the authors reviewed the manuscript.
Additional Information
Supplementary information accompanies this paper at http://www.nature.com/srep
Competing financial interests: The authors declare no competing financial interests.
How to cite this article: Zhang, Y. et al. A Comprehensive Numerical Model for Simulating Fluid Transport in Nanopores. Sci. Rep. 7, 40507; doi: 10.1038/srep40507 (2017).
Correspondence and requests for materials should be addressed to W.Y. (email: yuwei127@gmail.com) or Y.D. (email: diyuan@mech.pku.edu.cn)
References
 Nojabaei, B., Johns, R. T. & Chu, L. Effect of capillary pressure on phase behavior in tight rocks and shales. Presented at the SPEAnnual Technical Conference and Exhibition, San Antonio, Texas, USA, October 8–10. SPE 159258MS (2012).
 Cipolla, C. L., Warpinski, N. R., Mayerhofer, M., Lolon, E. P. & Vincent, M. The relationship between fracture complexity, reservoir properties, and fracturetreatment design. Presented at the SPE Annual Technical Conference and Exhibition, Denver, Colorado, USA, September 21–24. SPE 115769MS (2008).
 Wu, K. & Olson, J. Numerical investigation of complex hydraulicfracture development in naturally fractured reservoirs. SPE Production & Operations. In press (2016).
 Nelson, P. H. Porethroat sizes in sandstones, tight sandstones, and shales. AAPG Bulletin 93, 329–340 (2009).
 Kuila, U. & Prasad, M. Specific surface area and poresize distribution in clays and shales. Geophysical Prospecting 61, 341–362 (2013).
 Brusilovsky, A. I. Mathematical simulation of phase behavior of natural multicomponent systems at high pressures with an equation of state. SPE Reservoir Engineering 7, 117–122 (1992).
 Guo, P., Sun, L., Li, S. & Sun, L. Mathematical simulation of phase behavior of natural multicomponent systems at high pressures with an equation of state. Presented at the SPE Gas Technology Symposium Calgary, Alberta, Canada, April 28May 1. SPE 35644MS (1996).
 Pang, J., Zuo, J. Y., Zhang, D. & Du, L. Impact of porous media on saturation pressures of gas and oil in tight reservoirs. Presented at the SPE Canadian Unconventional Resources Conference, Calgary, Alberta, Canada, October 30November 1. SPE 161143MS (2012).
 Nojabaei, B. & Johns R. T. Consistent extrapolation of black and volatile oil fluid properties above original saturation pressure for tight and conventional reservoirs. Presented at the SPE Reservoir Simulation Symposium Houston, Texas, USA. February 23–25. SPE173221MS (2015).
 Du, L. & Chu L. Understanding anomalous phase behavior in unconventional oil reservoirs. Presented at the SPE CanadianUnconventional Resources Conference, Calgary, Alberta, Canada, October 30November 1. SPE 161830MS (2012).
 Wang, Y., Yan, B. & Killough, J. Compositional modeling of tight oil using dynamic nanopore properties. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, USA, September 30 October 2. SPE 166267MS (2013).
 Nojabaei, B., Siripatrachai, N., Johns, R. T. & Ertekin, T. Effect of saturation dependent capillary pressure on production in tight rocks and shales: A compositionallyextended black oil formulation. Presented at the SPE Eastern Regional Meeting Charleston, WV, USA, October 21–23. SPE 171028MS (2014).
 Rezaveisi, M., Sepehrnoori, K. & Pope, G. A. Compositional simulation including effect of capillary pressure on phase behavior.Presented at the SPE Annual Technical Conference and Exhibition Houston, Texas, USA, SPE 175135MS (2015).
 ZuloagaMolero, P., Yu, W., Xu, Y., Sepehrnoori, K. & Li, B. Simulation Study of CO2EOR in Tight Oil Reservoirs with Complex Fracture Geometries. Scientific Reports. 6, 33445; 10.1038/srep33445 (2016).
 Wu, T. & Zhang D. Impact of Adsorption on Gas Transport in Nanopores. Scientific Reports. 6, 23629; 10.1038/srep23629 (2016).
 Wang, S., Javadpour, F. & Feng, Q. Confinement Correction to Mercury Intrusion Capillary Pressure of Shale Nanopores. Scientific Reports 6, 20160, doi: 10.1038/srep20160 (2016).
 Siripatrachai, N., Ertekin, T. & Johns, R. T. Compositional simulation of discrete fractures incorporating the effect of capillary pressure on phase behavior. Presented at the SPE Improved Oil Recovery Conference, Tulsa, Oklahoma, USA, April 11–13. SPE179660MS (2016).
 Rubin, B. Accurate simulation of nonDarcy flow in stimulated fractured shale reservoirs. Presented at the SPE Western Regional Meeting Anaheim, California, USA, May 27–29. SPE 132093MS (2010).
 Nojabaei, B., Johns, R. T. & Chu, L. Effect of capillary pressure on phase behavior in tight rocks and shales. SPE Reservoir Evaluation & Engineering. 16, 281–289 (2013).
 Yu, W. & Sepehrnoori, K. Optimization of well spacing for Bakken tight oil reservoirs. Presented at the SPE Unconventional Resources Technology Conference, Denver, Colorado, USA, August 25–27. SPE 20141922108 MS (2014).
 Zhang, Y., Lashgari, H. R., Di, Y. & Sepehrnoori, K. Capillary pressure effect on hydrocarbon phase behavior in unconventional reservoirs. Presented at the SPE Low Perm Symposium Denver, Colorado, USA, May 5–6. SPE 180235MS (2016).
 Yu, W., Lashgari, H. R. & Sepehrnoori, K. Simulation study of CO2 huffnpuff process in Bakken tight oil reservoirs. Presented at theWestern North American and Rocky Mountain Joint Meeting Denver, Colorado, USA, Apri 17–18. SPE 169575MS (2014).
 Jin, L. et al. Improving oil transportability using CO2 in the Bakken system–A laboratory investigation. Presented at the SPE International Conference and Exhibition on Formation Damage Control. Lafayette, Louisiana, USA, February 24–26. SPE 178948MS (2016).
 Teklu, T. W., Zhou, Z., Li, X. & Abass, H. Experimental investigation on permeability and porosity hysteresis in lowpermeability formations. Presented at the SPE Low Perm Symposium Denver, Colorado, USA, May 5–6. SPE 180226MS (2016).
 Kurtoglu, B., Sorensen, J. , Braunberger, J., Smith, S. & Kazemi, H. Geologic characterization of a Bakken reservoir for potential CO2 EOR. Presented at Unconventional Resources Technology Conference, Denver, Colorado, USA, August 12–14. SPE 168915MS (2013).
 Xu, Y., De Araujo Cavalcante Filho, J. S., Yu, W. & Sepehrnoori, K. Discrete fracture modeling of complex hydraulic fracture geometries in reservoir simulators. SPE Reservoir Evaluation & Engineering. In press (2016).
 Sorensen, J. A. et al. Characterization and evaluation of the Bakken petroleum system for CO2 storage and enhanced oil recovery.Presented at the Unconventional Resources Technology Conference, San Antonio, Texas, USA, July 20–22. SPE 178659MS (2015).
 Adamson, A. W. Physical Chemistry of Surfaces (Fifth Edition) 6–7 (John Wiley & Sons, 1990).
 Firoozabadi, A. Thermodynamics of hydrocarbon reservoirs 55–60 (New York: McGrawHill, 1999).
 Fong, D. K. S. & Nghiem, L. X. A viscosity model for reservoir fluid. Computer Modeling Group Research Report R7, (1980).
 Stiel, L. I. & Thodos, G. The viscosity of nonpolar gases at normal pressures. AIChE Journal. 7, 611–615 (1961).
 Du, S., Yoshida, N., Liang, B. & Chen, L. Application of multisegment well approach: Dynamic modeling of hydraulic fractures. Journal of Natural Gas Science and Engineering. In press (2016).
 Li, L. & Lee, S. H. Efficient fieldscale simulation of black oil in a naturally fractured reservoir through discrete fracture networks and homogenized media. SPE Reservoir Evaluation & Engineering. 11, 750–758 (2008).
 Moinfar, A., Varavei, A. Sepehrnoori, K. & Johns, R. T. Development of an efficient embedded discrete fracture model for 3D compositional reservoir simulation in fractured reservoirs. SPE Journal. 19, 289–303 (2014).
 De Araujo Cavalcante Filho, J. S., Shakiba, M., Moinfar, A. & Sepehrnoori, K. Implementation of a preprocessor for embedded discrete fracture modeling in an IMPEC compositional reservoir simulator. Presented at the SPE Reservoir Simulation Symposium Houston, Texas, USA. February 23–25. SPE 173289MS (2015).
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
© The Author(s) 2017