Research Article  Open Access
Yan Zeng, Qing Wang, Zhengfu Ning, Hongliang Sun, "A Mathematical Pressure Transient Analysis Model for Multiple Fractured Horizontal Wells in Shale Gas Reservoirs", Geofluids, vol. 2018, Article ID 8065949, 16 pages, 2018. https://doi.org/10.1155/2018/8065949
A Mathematical Pressure Transient Analysis Model for Multiple Fractured Horizontal Wells in Shale Gas Reservoirs
Abstract
Multistage fractured horizontal wells (MFHWs) have become the main technology for shale gas exploration. However, the existing models have neglected the percolation mechanism in nanopores of organic matter and failed to consider the differences among the reservoir properties in different areas. On that account, in this study, a modified apparent permeability model was proposed describing gas flow in shale gas reservoirs by integrating bulk gas flow in nanopores and gas desorption from nanopores. The apparent permeability was introduced into the macroseepage model to establish a dynamic pressure analysis model for MFHWs dualporosity formations. The Laplace transformation and the regular perturbation method were used to obtain an analytical solution. The influences of fracture halflength, fracture permeability, Langmuir volume, matrix radius, matrix permeability, and induced fracture permeability on pressure and production were discussed. Results show that fracture halflength, fracture permeability, and induced fracture permeability exert a significant influence on production. A larger Langmuir volume results in a smaller pressure and pressure derivative. An increase in matrix permeability increases the production rate. Besides, this model fits the actual field data relatively well. It has a reliable theoretical foundation and can preferably describe the dynamic changes of pressure in the exploration process.
1. Introduction
Shale gas is known as a key resource to meet the increasing world energy demand because of its rich reserves and extensive distribution [1]. To be more specific, many pieces of research have shown that shale gas reservoirs are characterized by extremely low permeability values. Consequently, hydraulic fracturing simulations and horizontal drilling are developed for commercial exploitation. In shale gas reservoirs, the complex pore structure, macro/microfracture network distribution, and microscale reservoir properties make the gas flow characteristics become more complicated than those in conventional reservoirs. Gas flow in organic matter, natural fractures, secondary fractures, and hydraulic fractures is controlled by different mechanisms [2]. To be more specific, 20%–80% of shale gas exists in the adsorbed state on the surface of organic particulates. Desorption of adsorbed gas in organic matter affects gas transport mechanisms significantly [3, 4]. Understanding the pressure and ratetransient behaviors of multistage fractured horizontal wells in shale gas reservoirs is of great importance for production forecasting, well placement, and configuration optimization [5].
In recent studies [6–8], gas transport mechanisms in organic matter can be described by continuum flow, slip flow viscous flow, the Knudsen diffusion, and transition flow. Wang and Li [9] made a comparison of collision frequency in gasgas and gassolid interactions. They concluded that the collisions between gas molecules and the solid wall account for a large percentage of the total collisions. Therefore, the Knudsen diffusion should be considered when studying gas transport in nanopores. Darabi et al. [10] also reported that, under typical shale gas reservoir conditions, the Knudsen diffusion dominates gas transport and the contribution to cumulative production can reach 20%. Adsorbed gas desorption from the surface of organic matter also plays a significant role in gas transport [11, 12]. Swami et al. [13] developed a modified model to investigate the contribution of adsorbed gas to the final recovery. They found that desorption of adsorbed gas can increase the pore diameter, reduce tortuosity, and cause extra slippage at the solid boundary. Yu and Sepehrnoori [14] analyzed the shale gas development history for a period of 30 years in North America and concluded that different shale samples demonstrate quite different adsorption capacities.
Mayerhofer et al. [15] studied the difference between conventional double wings symmetric fractures and hydraulic fractures of Barnet shale and were the first to propose the simulated reservoir volume (SRV) conception. In SRV, multistage hydraulic fracturing connects existing natural fractures, which generates a large fracture network [16]. Kucuk and Sawyer [17] firstly developed an analytical model to investigate transient pressure in shale gas reservoirs. However, the model ignored diffusion flow and the desorption effect. Lee and Brockenbrough [18] analyzed production data coming from vertical wells and developed an analytical solution to describe transient pressure. Ozkan et al. [19] proposed a dualmechanism, that is to say, a dualporosity horizontal well model for shale gas reservoirs, which includes hydraulic fractures, as well as the inner region and the outer region. The model can adopt the diffusive flow mechanism but fail to take adsorbed gas into consideration. Stalgorova and Mattar [20–22] extended the trilinear flow model to a fiveregion model for MFHWs in homogeneous shale gas reservoirs. In order to investigate the effect of natural fractures on matrix permeability, Apaydin et al. [23] combined composite blocks into an analytical trilinear flow model. Zhao et al. [24, 25] analyzed the pressure response and production performance of MFHWs in shale gas reservoirs and they introduced the source function theory to characterize SRV. Liu et al. [26] discussed the desorption and diffusion effects in the analytical model and studied the productivitydecline law. Wang [27] considered the stress sensitive effect of natural fractures and hydraulic fracture angles in the semianalytical model. His research implied that the stress sensitive effect has a significant impact on bottom hole flowing pressure.
In recent years, many general semianalytical models, considering microscopic seepage mechanisms, have been proposed with the rapid development of computer technology [24, 25, 28–30]. These models can describe predominant flow regimes for shale well but fail to consider all of the typical characteristics and possible situations in shale gas reservoirs. In the present study, a novel apparent permeability model was put forward by taking into consideration viscous flow, slip flow, transition flow, gas adsorption/desorption [31], and the poromechanical response [32, 33]. The novelty of this study concentrates on the improvement of the semianalytical model to conduct pressure and ratetransient analysis in SRV, by considering microscale percolation mechanisms and heterogeneity. The rest of this study can be described in the following way. In Section 2, the apparent permeability model for gas transport in shale nanopores is introduced. The semianalytical model for MFHWs in a heterogeneous shale gas reservoir, by considering microscale percolation mechanisms in a dualporosity formation, is introduced in Section 3. Loglog dimension pressure, dimension pressure derivative, and production type curves are plotted firstly, and then a sensitivity analysis of formation parameters is conducted in Section 4. An actual field case is studied in Section 5. Finally, some conclusions are drawn in Section 6.
2. Apparent Permeability Model for Gas Transport in Shale Nanopores
2.1. Viscous Flow
The total shale gas flux is composed of a viscous flow flux and a Knudsen diffusion flux component. Viscous flow is strongly prominent when the short molecule free path is far less than the pore diameter. In other words, the transport mechanism is governed by viscous flow in macropores (where the pore diameter > 50 nm). The HagenPoiseuille equation can describe the molar flux [34] as stated below:where is the viscous flow flux, mol/(m^{2}s); is porosity, dimensionless; is tortuosity, dimensionless; is the pore radius, m; is viscosity, Pas; is the universal gas constant, J/(molK); is temperature, K; is pressure, Pa; and is the shale gas transport distance, m.
Gas viscosity gradually deviates from the traditional viscosity definition for high Knudsen numbers. Karniadakis et al. [35] modified the viscosity definition in the Knudsen layer by considering the rarefaction effect.where is the viscosity considering rarefaction effect, Pas; is the rarefaction coefficient, dimensionless; and Kn is the Knudsen number, dimensionless.
According to the Darcy equation, the apparent permeability of viscous flow can be expressed as [36]:where is the apparent permeability, taking into account the effect of rarefaction, m^{2}; is the gas mole volume at standard conditions, 22.414 10^{−3} m^{3}/mol; and is the permeability of viscous flow, m^{2}.
2.2. The Knudsen Diffusion
When there is shale gas transport through micropores under low pressure, the molecule free path is long and equal to the pore diameter. The Knudsen diffusion is prominent because collisions between molecules and the micropores wall are more frequent than the intermolecular collisions. The flux of the Knudsen diffusion can be expressed as [34]where is the Knudsen diffusion flux, mol/(m^{2}s); is the Knudsen diffusion coefficient, m^{2}/s; and is the gas concentration, m^{3}/mol. The Knudsen diffusion coefficient is obtained aswhere is the gas molar mass, kg/mol. The rough wall has a significant effect on the Knudsen diffusion. Darabi et al. [10] proposed a function to describe the roughness effect as follows:where is the Knudsen diffusion coefficient which takes into consideration the roughness effect, m^{2}/s; is the mean free path of gas molecules, m; and is the fractal dimension of the rough wall, dimensionless.
By combining (4) and (6), the apparent permeability of the Knudsen diffusion can be obtained and written as follows:
2.3. Weight Factor
In the actual shale gas reservoirs, different transport mechanisms exist at the same time. Therefore, the total shale gas flux should be a weighted summation of the viscous flow flux and the Knudsen diffusion flux, based on their different contributions. The weight factors of the viscous flow and the Knudsen diffusion flow are defined as the ratios of intermolecular collisions and molecular/porewall collisions to total collisions, respectively. Based on the definitions for a gas molecule free path [11] and collision numbers [37], the weight factors of viscous flow and Knudsen diffusion can be approximated as follows:
Therefore, the total apparent permeability in nanopores is
With a decrease in pressure, the gas molecule free path increases immediately. The weight factor of the Knudsen diffusion also increases. The piece of research carried out by Wu et al. [38] showed that, under the 10 nm radius condition, viscous flow dominates when the pressure is larger than 4.02 10^{5} Pa; otherwise, the Knudsen diffusion dominates.
2.4. SorptionInduced Swelling Response
Organic matter is described by weak strength and strong sensitivity to stress change. With an increase in reservoir pressure, adsorbed gas begins to desorb. Gas desorption results in shrinkage of the organic matrix and to an increased effective hydraulic diameter [39]. According to the solid deformation theory [40] and the Langmuir isotherm equation [41], the relationship between the degree of solid deformation and the reservoir pressure can be written as follows:where is the solid deformation degree, dimensionless; is the Langmuir volume, m^{3}/kg; is the density of the shale matrix, kg/m^{3}; is the shale matrix Young modulus, Pa; is the Langmuir pressure, Pa; is the initial pressure, Pa; and is the rock compressibility, Pa^{−1}.
Meanwhile, based on Seidle’s model [42], the approximation relation between the porosity and the solid deformation degree is given by
Assume that the pore volumes are proportional to the gas flow channels. In accordance with the capillary model, the effective hydraulic diameter can be obtained as
2.5. Gas Adsorption/Desorption
Because of the large surface area and oilwet characteristic, nanopores have a strong adsorbed gas capacity. The mass balance equation, which considers the adsorbed gas, can be expressed aswhere is the shale gas density, kg/m^{3}; is the gas compressibility, Pa^{−1}; and is the mass of adsorbed gas per unit volume of shale, kg/m^{3}. The apparent permeability, which considers desorption of absorbed gas, can be expressed as
Cui et al. [12] defined the effective adsorption porosity as
According to the Langmuir isotherm equation, the effective adsorption porosity can be expressed as
Finally, the apparent permeability model for gas transport in shale nanopores can be expressed as
3. Productivity Prediction Model for MFHWs
3.1. Physical Model
Figure 1 shows the diagram of MFHWs in a shale gas reservoir. The reservoir is divided into seven contiguous regions: two upperreservoir regions (regions 6 and 7), two outerreservoir regions (regions 4 and 5), two innerreservoir regions (regions 2 and 3), and a hydraulic region (region 1). Not all reservoir zones between primary hydraulic fractures are stimulated by induced fractures. 1D linear flow is assumed within unstimulated region and the flow direction depends on the location. As shown in Figure 1, a region of higher permeability around each fracture, the socalled SRV area, is introduced to represent fracture branching (region 2). The SRV region occupies part of the space between fractures and the flow in the other part (region 3) towards region 2, parallel to the wellbore. Regions 4 and 5 are in the outer reservoir beyond the tips of the hydraulic fractures and they connect with regions 2 and 3 in direction. Because of the higher permeability in region 2, the pressure decreases in this area first. Therefore, the flow direction in regions 4 and 5 is towards regions 2 and 3, perpendicular to the wellbore. Although the main extending direction of the hydraulic fractures is horizontal ( direction in Figure 1), we also need to consider fractures height in vertical direction ( direction in Figure 1). Based on this consideration, regions 6 and 7 are defined in the upper reservoir beyond the tips of the hydraulic (connecting with regions 2 and 3 in direction) and the flow direction is vertical to the wellbore, towards the lower reservoir.
The main assumptions considered for this model are listed below:(1)The shale reservoir has been fractured and has a constant reservoir thickness. All hydraulic fractures are symmetric and arranged uniformly along the horizontal well. Length and conductivity of the hydraulic fractures are the same.(2)Natural fracture permeability is stressdependent and shale gas flow in the matrix is driven by concentration difference.(3)The initial pressure throughout the reservoir is uniform.(4)Hydraulic fractures penetrate the entire pay formation. Fracture interference is ignored. After SRV, large fracture networks can be generated in region 2 and it has higher natural fracture permeability compared to the other regions.(5)The shale gas flow in the reservoir is considered to be an isothermal flow and the effect of gravity and capillary pressure are negligible.(6)Nonsimulated areas are described by a dualporosity system consisting of a shale matrix and natural fractures.(7)The only route for shale gas to reach the horizontal wellbore is through hydraulic fractures. The flow rate of MFHW is the sum of the flow rate from every hydraulic fracture.
3.2. Mathematical Model
According to the mass conservation equation, the motion equation, and the state equation, the partial differential formulas for shale gas flow in a matrix system and a fracture system can be obtained as follows [19, 27]:where is a characteristic dimensionless parameter of stress sensitivity and relates to the permeability modulus . is defined as , where .
Pedrosa’s substitution [43] is used to weaken the nonlinearity of (19). By setting , (19) can be simplified into
To solve (18) and (20), these two equations were converted to the Laplace domain. And then the definite boundary continuity is introduced to develop an integral seepage flow differential equation. The final equation can be written aswhere .
A series of variables and dimensionless variables are defined in Table 1.

3.2.1. Region 7
In this upperreservoir region, the governing diffusivity equation for shale gas flow in the direction becomeswhere
At the top of the reservoir , the boundary condition is noflow:
On the interfaces of regions 7 and 3 and regions 7 and 5, pressure continuity condition can be obtained as
Equation (22) is solved simultaneously with (24) and (25). The final solution for (22) can be given as
3.2.2. Region 6
Consistent with region 7, the governing diffusivity equation for shale gas flow in the direction becomeswhere
The boundary condition at the top of the reservoir is
The boundary conditions on the interfaces of regions 6 and 2 and regions 6 and 4 are
The solution for (27) can be given as
3.2.3. Region 5
In this region, shale gas flow takes place in the direction and the direction. So the governing diffusivity equation for region 5 iswhere
Along the horizontal well , there is a noflow boundary:
Flux across the interface of regions 5 and 7 is continuous; therefore,
Substituting (34) and (35) into (32), each of its terms from 0 to is integrated with respect to . The governing diffusivity can be simplified as
In this region, shale gas flow in the direction and with a closed boundary condition is expressed as
Pressure continuity condition between regions 5 and 3 (at ) is expressed as
By solving (36), the solution can be obtained aswhere
3.2.4. Region 4
In a similar way, shale gas flow in this region takes place in the direction and the direction. The governing diffusivity equation for region 4 iswhere
The same method is also used in region 5, where the governing diffusivity equation can be simplified as
Noflow condition at the outerreservoir boundary is given as
On the interface of regions 4 and 2, the boundary condition is
The solution for (41) is given aswhere
3.2.5. Region 3
This region is adjacent to regions 5 and 7, so shale gas flow direction takes place along the direction, direction, and direction. The governing diffusivity equation becomeswhere
For the direction, the boundary condition used is a noflow condition for a horizontal well and a flux continuity condition is used on the interface of regions 5 and 3 :
For the direction, (51) can be obtained:
By substituting (50) and (51) into (48), the final governing diffusivity equation becomes
Between two fractures, there is a noflow boundary (at ):
On the interface of regions 3 and 2, the pressure continuity condition is written as
The solution for (48) is given aswhere
3.2.6. Region 2
Similarly, as region 2 is adjacent to hydraulic fracture, the governing diffusivity equation is
Consistent with region 3, the simplified equation can be written aswhere
At the interface of regions 2 and 3, the boundary condition is characterized by flux continuity:
And between the hydraulic fracture level and region 2, the boundary condition is
3.2.7. Region 1
In this region, the diffusivity equation is given bywhere . The boundary condition can be stated by the expressions:
The solution can be obtained:where
Assuming that pressure distribution along the horizontal well is uniform, the dimensionless downhole pressure can be solved by setting [44]:
Shale gas flow in hydraulic fractures is not completely linear. The nonlinear streamlines result in an additional pressure drop. Mukherjee and Economides [45] presented a formulation to take this additional pressure drop into consideration.
Based on the Duhamel principle and the superposition principle, the dimensionless flowing bottom hole pressure (FHBP), which takes into consideration the wellbore storage effect and the skin effect, can be expressed as [30]
The FHBP for MFHWs, when considering the stress sensitivity effect, can be obtained as
The relationship between the dimensionless rate under the condition of constant pressure and dimensionless pressure caused by the constant rate is described as [46]
4. Discussion and Analysis
4.1. Apparent Permeability in the Shale Matrix
Figure 2 shows the variation of the apparent permeabilities against the average pressure. The input parameters are listed in Table 2. Based on Figure 2, one can observe that, in the early period, with a decrease of the average pressure, the bulk permeability decreases dramatically, but the apparent permeability caused by gas desorption increases greatly. In the later period, increases slightly with an increase of the pressure. also increases in a similar way with an increase in pressure. The explanation is that, at high pressure, the Knudsen diffusion increases significantly. As can be seen in Figure 2, viscous flow dominates at high pressure and the Knudsen diffusion dominates at low pressure. Figure 2 also indicates that the pore diameter has a great impact on the apparent permeability. With a decrease of the pressure, nanopores become larger and larger under the influence of matrix shrinkage, so the apparent permeability increases gradually. However, there is an inflection point between the apparent permeability and the average pressure. Furthermore, pore diameter is bigger and the time of the inflection point for apparent permeability occurring appears later. The apparent permeability decreases faster after the inflection point. The results presented in Figure 2 agree well with those presented by Wu et al. [38].

(a)
(b)
(c)
4.2. Effects of Hydraulic Fracture HalfLength
Figures 3 and 4 show the effect of hydraulic fracture halflength on pseudopressure, pseudopressure derivative, and dimensionless production. The input parameters are listed in Table 3. According to the results, fracture halflength affects the whole production process. The longer the hydraulic fracture halflength, the more significant the pressure drop and the higher the production rate. The contact area between the hydraulic fracture and the shale gas reservoir is larger when the fracture halflength is longer. Therefore, when gas flow resistance is reduced, more shale gas can flow to the fracture. Furthermore, larger flow rate results in significant pressure drop and a larger production.

4.3. Effect of Hydraulic Fracture Permeability
Figures 5 and 6 present the curve of the relationships of dimensionless time and pseudopressure and dimensionless production rate under different hydraulic fracture permeability. The figures show that hydraulic fracture permeability dominates all of the production period except for when we have early pure wellbore storage regime and boundary flow regime. The results are quite different from those observed in a single porosity tight reservoir where hydraulic fracture permeability only exerts an influence on early time pressure behavior [47]. The figures also indicate that, with an increase in fracture permeability, the pseudopressure drop and the production rate growth tend to slow down. Therefore, there is an optimal fracture permeability for a certain reservoir.
4.4. Effect of the Langmuir Volume.
The effects of the Langmuir volume on pressure and production are plotted in Figures 7 and 8. The Langmuir volume is defined as the maximum volume of shale gas that can be absorbed to a unit weight organic matrix at infinite pressure. The Langmuir volume is a significant parameter for gas adsorption and desorption. As illustrated in Figures 7 and 8, the larger the Langmuir volume is, the lower the pressure drop and the production will be. This is because a larger Langmuir volume results in a greater amount of absorbed gas which provides an extra contribution for well reserves. However, the Langmuir volume does not affect the linear flow stage as only shale gas in hydraulic fractures contributes to production at that period. It can also be seen that, with an increase of the Langmuir volume, the dimensionless rate decreases slowly. This results in longer production time.
4.5. Effect of Matrix Radius
Figures 9 and 10 show the effect of matrix block radius on pseudopressure, pseudopressure derivative, and dimensionless production rate. The shape and size of the matrix blocks exert a very significant effect on the interporosity flow stage. In this model, the value of the matrix block radius determines the magnitudes of the interporosity flow coefficient and the storativity ratio. According to the results, the matrix block radius mainly affects the transitional stage and the stage of bilinear flow. That is to say, the larger the matrix block radius, the more significant the pressure drop of the transitional stage. Consequently, the stage of bilinear flow occurs later. When the matrix block radius is large enough, bilinear flow disappears. Meanwhile, the magnitude of the matrix block radius also affects the production in the early and middle stages. With an increasing matrix block radius, the production rate gradually decreases, and the production decline rate decreases as well. In the late stage of production, the production rates under different matrix block radii are approximately the same [30].
4.6. Effect of Matrix Permeability
Figures 11 and 12 demonstrate how matrix permeability affects pseudopressure and dimensionless production rate. As can be seen in these two figures, as the matrix permeability decreases, pressure drop in the transient period becomes more significant and the fracture linear flow period becomes longer. Dimensionless production rate increases slowly with an increase of matrix permeability. The results indicate that the rock has more capacity to transfer gas molecules from the matrix to the fracture system. Under ultralow matrix permeability conditions, the pressure variation is more sensitive.
4.7. Effect of Induced Fracture Permeability
Figures 13 and 14 represent the effect of induced fracture permeability on pressure and production. The induced fracture permeability of shale reservoirs is an extremely important reservoir parameter, as shale gas transport to hydraulic fracture should take place through induced fracture systems. Therefore, the induced fracture permeability determines the scale of gas well productivity. The results show that the higher the induced fracture permeability, the more significant the pressure drop during the transitional period. With a continuous increase in induced fracture permeability, the duration of the stage of bilinear flow becomes shorter and shorter until it completely disappears. Figure 13 indicates that the induced fracture permeability has a more significant influence on pressure than matrix radius and matrix permeability. Meanwhile, the induced fracture permeability also affects the scale of production rate; that is, the higher the permeability, the larger the scale. However, as far as production rate is concerned, the secondary fracture permeability only affects its scale, but not its decline rate.
5. Case Study
The MFHW is located in the transitional zone between the Sichuan Basin and the YunnanGuizhou Plateau (400–1,200 m in altitude) and mainly developed a marine and continental sedimentary facies. The specific stratum parameters are listed in Table 4.

As seen in Figure 15, the production data show that there is a relatively high degree of correlation between pressure and production; the production data present an obvious response to the pressure fluctuations caused by the test carried out in the early stage; therefore, the data have a certain degree of reliability. As shown by a comparative analysis on the actual production data and production curve, the measured data are deviated from the typical curve by a relatively significant extent. This is mainly because of the skin effect and wellbore storage effect on the early stage seepage. The fitting curve shows that the data points are concentrated in the stage of linear flow of matrix (see Figure 16), mainly because the extremely low permeability of matrixfractures has extended the duration of the stage of nonsteady state linear flow.
6. Conclusion
In this study, an analytical model was presented to simulate shale gas flow through MFHW in a shale reservoir. Based on microscale percolation mechanisms in nanopores, a sevenlinearflow model for dualporosity formations was established. This model takes viscous flow, the Knudsen diffusion, the sorptioninduced swelling response, the adsorption/desorption, the matrix shrinkage and the stress sensitivity effect into consideration. The pressure transient analysis was carried out under conditions of constant rate. The following conclusions were deduced:
The shale gas transport mechanism in nanopores mainly includes bulk gas transport and gas adsorption/desorption. Bulk gas apparent permeability is a weighted summation of viscous flow permeability and the Knudsen diffusion permeability. Gas desorption can increase the total apparent permeability.
A semianalytical model was established for the pressure behavior of MFHWs in shale gas reservoirs. Based on the proposed model, flow regimes were identified into two top reservoir regions, two innerreservoir regions, two outerreservoir regions, and one artificial fracture region.
Larger fracture halflength results in a more significant pressure drop and a higher production rate. Fracture permeability dominates throughout the production period except during the early pure wellbore storage stage and the boundary flow stage. The effect of desorption cannot be neglected when simulating, as the Langmuir volume has a significant influence on the production rate. Compared with the matrix radius and matrix permeability, pressure and production are more sensitive to induced fracture permeability.
The semianalytical model fits relatively well with the actual production data, has a reliable theoretical foundation, and can preferably describe the dynamic changes of pressure in the exploration process of shale gas wells.
Nomenclature
:  Temperature, K 
:  Shale tortuosity 
:  Gas molar mass, kg/mol 
:  Universal gas constant, J/(molK) 
:  Boltzmann constant, J/K 
:  Shale density, t/m^{3} 
:  Gas density, t/m^{3} 
:  Shale compressibility, MPa^{−1} 
:  Gas compressibility, MPa^{−1} 
:  Mole volume of gas at standard constant, m^{3}/mol 
:  Langmuir volume, m^{3}/t 
:  Langmuir pressure, MPa 
:  Shale matrix Young modulus, MPa 
:  Rarefaction coefficient 
:  Pseudopressure, psi^{2}/cp 
:  Production rate, Mscf/D 
:  Laplace transform parameter 
:  Time, hours 
:  Permeability modulus, psi^{−1} 
:  Coordinate, ft 
:  Coordinate, ft 
:  Coordinate, ft 
:  Hydraulic fracture width, ft 
:  Hydraulic fracture halflength, ft 
:  Reservoir size in direction, ft 
:  SRV width, ft 
:  Halfwidth between two fractures, ft 
:  Half fracture height, ft 
:  Half reservoir thickness, ft 
:  Radius, ft 
:  Shale viscosity, cp 
:  Porosity 
:  Shape factor 
:  Diffusivity, ft/hour 
:  Flow capacity ratio 
:  Desorption coefficient at average pressure 
:  Storativity ratio 
:  Skin factor 
:  Dimensionless wellbore storage coefficient. 
:  Dimensionless 
:  Natural fracture system 
:  Hydraulic fracture 
:  The th region (i = 2, 3, 4, 5, 6, 7) 
:  Matrix system 
sc:  Standard condition 
:  Dimensionless wellbore 
:  Dimensionless length of external boundary. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to acknowledge the support from MOE Key Laboratory of Petroleum Engineering in China University of Petroleum (Beijing). The authors would also like to acknowledge the financial support from the National Nature Science Foundation of China (no. 51504265) and Science Foundation of China University of Petroleum, Beijing (no. 2462015YQ0223).
References
 R. Heller, J. Vermylen, and M. Zoback, “Experimental investigation of matrix permeability of gas shales,” AAPG Bulletin, vol. 98, no. 5, pp. 975–995, 2014. View at: Publisher Site  Google Scholar
 J. B. Curtis, “Fractured shalegas systems,” AAPG Bulletin, vol. 86, no. 11, pp. 1921–1938, 2002. View at: Google Scholar
 F. Javadpour, “Nanopores and apparent permeability of gas flow in mudrocks (shales and siltstone),” Journal of Canadian Petroleum Technology, vol. 48, no. 8, pp. 16–21, 2009. View at: Google Scholar
 F. Yang, Z. Ning, Q. Wang, R. Zhang, and B. M. Krooss, “Pore structure characteristics of lower Silurian shales in the southern Sichuan Basin, China: Insights to pore development and gas storage mechanism,” International Journal of Coal Geology, vol. 156, pp. 12–24, 2016. View at: Publisher Site  Google Scholar
 F. Civan, D. Devegowda, and R. F. Sigal, “Critical evaluation and improvement of methods for determination of matrix permeability of shale,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 4620–4629, Society of Petroleum Engineers, October 2013. View at: Google Scholar
 W. N. Yuan, Z. J. Pan, X. Li et al., “Experimental study and modelling of methane adsorption and diffusion in shale,” Fuel, vol. 117, pp. 509–519, 2014. View at: Publisher Site  Google Scholar
 R. Aguilera, “Flow units: from conventional to tight gas to shale gas reservoirs,” in Proceedings of the Trinidad and Tobago Energy Resources Conference, pp. 572–587, Society of Petroleum Engineers, June 2010. View at: Google Scholar
 M. R. Rahmanian, R. Aguilera, and A. Kantzas, “A new unified diffusionviscousflow model based on porelevel studies of tight gas formations,” SPE Journal, vol. 18, no. 1, pp. 38–49, 2012. View at: Publisher Site  Google Scholar
 M. Wang and Z. Li, “Nonideal gas flow and heat transfer in micro and nanochannels using the direct simulation Monte Carlo method,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 68, no. 4, Article ID 046704, 2003. View at: Publisher Site  Google Scholar
 H. Darabi, A. Ettehad, F. Javadpour, and K. Sepehrnoori, “Gas flow in ultratight shale strata,” Journal of Fluid Mechanics, vol. 710, pp. 641–658, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 F. Javadpour, D. Fisher, and M. Unsworth, “Nanoscale gas flow in shale gas sediments,” Journal of Canadian Petroleum Technology, vol. 46, no. 10, pp. 55–61, 2007. View at: Google Scholar
 X. Cui, A. M. M. Bustin, and R. M. Bustin, “Measurements of gas permeability and diffusivity of tight reservoir rocks: Different approaches and their applications,” Geofluids, vol. 9, no. 3, pp. 208–223, 2009. View at: Publisher Site  Google Scholar
 V. Swami, A. T. Settari, and F. Javadpour, “A numerical model for multimechanism flow in shale gas reservoirs with application to laboratory scale testing,” in Proceedings of the SPE EAGE Annual Conference & Exhibition incorporating (EUROPEC '13), pp. 1274–1293, Society of Petroleum Engineers, June 2013. View at: Google Scholar
 W. Yu and K. Sepehrnoori, “Simulation of gas desorption and geomechanics effects for unconventional gas reservoirs,” Fuel, vol. 116, pp. 455–464, 2014. View at: Publisher Site  Google Scholar
 M. J. Mayerhofer, E. P. Lolon, J. E. Youngblood, and J. R. Heinze, “Integration of microseismic fracture mapping results with numerical fracture network production modeling in the Barnett shale,” in Proceedings of the SPE Annual Technical Conference and Exhibition (ATCE '06), pp. 976–983, Society of Petroleum Engineers, San Antonio, Tex, USA, September 2006. View at: Google Scholar
 C. R. Clarkson, “Production data analysis of unconventional gas wells: review of theory and best practices,” International Journal of Coal Geology, vol. 109, pp. 101–146, 2013. View at: Publisher Site  Google Scholar
 F. Kucuk and W. K. Sawyer, “Transient flow in naturally fractured reservoirs and its application to devonian gas shales,” in Proceedings of the SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, 1980. View at: Publisher Site  Google Scholar
 S.T. Lee and J. R. Brockenbrough, “A new approximate analytic solution for finiteconductivity vertical fractures,” SPE Formation Evaluation, vol. 1, no. 1, pp. 75–88, 1986. View at: Publisher Site  Google Scholar
 E. Ozkan, R. S. Raghavan, and O. G. Apaydin, “Modeling of fluid transfer from shale matrix to fracture network,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 3314–3331, Society of Petroleum Engineers, September 2010. View at: Google Scholar
 E. Stalgorova and L. Mattar, “Analytical model for unconventional multifractured composite systems,” SPE Reservoir Evaluation & Engineering, vol. 16, no. 3, pp. 246–256, 2013. View at: Publisher Site  Google Scholar
 E. Stalgorova and L. Mattar, “Analytical model for history matching and forecasting production in multifrac composite systems,” in Proceedings of the SPE Canadian Tight Resources Conference, pp. 450–466, Society of Petroleum Engineers, November 2012. View at: Google Scholar
 E. Stalgorova and L. Mattar, “Practical analytical model to simulate production of horizontal wells with branch fractures,” in Proceedings of the SPE Canadian Tight Resources Conference, pp. 433–449, Society of Petroleum Engineers, November 2012. View at: Google Scholar
 O. G. Apaydin, E. Ozkan, and R. Raghavan, “Effect of discontinuous microfractures on ultratight matrix permeability of a dualporosity medium,” SPE Reservoir Evaluation & Engineering, vol. 15, no. 4, pp. 473–485, 2012. View at: Publisher Site  Google Scholar
 Y.L. Zhao, L.H. Zhang, J.Z. Zhao, J.X. Luo, and B.N. Zhang, “’Triple porosity‘ modeling of transient well test and rate decline analysis for multifractured horizontal well in shale gas reservoirs,” Journal of Petroleum Science and Engineering, vol. 110, pp. 253–262, 2013. View at: Publisher Site  Google Scholar
 Y.L. Zhao, L.H. Zhang, J.X. Luo, and B.N. Zhang, “Performance of fractured horizontal well with stimulated reservoir volume in unconventional gas reservoir,” Journal of Hydrology, vol. 512, pp. 447–456, 2014. View at: Publisher Site  Google Scholar
 Q. Liu, K. Li, W. Wang, X. Hu, and H. Liu, “Production behavior of fractured horizontal well in closed rectangular shale gas reservoirs,” Mathematical Problems in Engineering, vol. 2016, Article ID 4260148, 9 pages, 2016. View at: Publisher Site  Google Scholar
 H.T. Wang, “Performance of multiple fractured horizontal wells in shale gas reservoirs with consideration of multiple mechanisms,” Journal of Hydrology, vol. 510, pp. 299–312, 2014. View at: Publisher Site  Google Scholar
 R.H. Zhang, L.H. Zhang, R.H. Wang, Y.L. Zhao, and D.L. Zhang, “Research on transient flow theory of a multiple fractured horizontal well in a composite shale gas reservoir based on the finiteelement method,” Journal of Natural Gas Science and Engineering, vol. 33, Article ID 66e80, pp. 587–598, 2016. View at: Publisher Site  Google Scholar
 R. H. Zhang, L. H. Zhang, R. H. Wang et al., “Simulation of a multistage fractured horizontal well with finite conductivity in composite Shale gas reservoir through finiteelement method,” Energy & Fuels, vol. 30, no. 11, pp. 9036–9049, 2016. View at: Google Scholar
 Y. Zeng, Z. Ning, Y. Lei, L. Huang, C. Lv, and Y. Hou, “Analytical model for shale gas transportation from matrix to fracture network,” in Proceedings of the SPE Europec featured at 79th EAGE Conference and Exhibition, Society of Petroleum Engineers, Paris, France, 2017. View at: Publisher Site  Google Scholar
 J. Wang, L. Chen, Q. Kang, and S. S. Rahman, “Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect,” Fuel, vol. 181, pp. 478–490, 2016. View at: Publisher Site  Google Scholar
 A. Tinni, E. Fathi, R. Agarwal, C. Sondergeld, Y. Akkutlu, and C. Rai, “Shale permeability measurements on plugs and crushed samples,” in Proceedings of the SPE Canadian Unconventional Resources Conference, pp. 342–355, Society of Petroleum Engineers, November 2012. View at: Google Scholar
 S. Wang, D. Elsworth, and J. Liu, “A mechanistic model for permeability evolution in fractured sorbing media,” Journal of Geophysical Research: Solid Earth, vol. 117, no. 6, 2012. View at: Publisher Site  Google Scholar
 J.G. Choi, D. D. Do, and H. D. Do, “Surface diffusion of adsorbed molecules in porous media: Monolayer, multilayer, and capillary condensation regimes,” Industrial & Engineering Chemistry Research, vol. 40, no. 19, pp. 4005–4031, 2001. View at: Publisher Site  Google Scholar
 G. Karniadakis, A. Beskok, and N. Aluru, Simple Fluids in Nanochannels, Springer, New York, NY, USA, 2005.
 S. Yuhong and W. K. Chan, “Analytical modeling of rarefied Poiseuille flow in microchannels,” Journal of Vacuum Science & Technology A: Vacuum, Surfaces, and Films, vol. 22, no. 2, pp. 383–394, 2004. View at: Publisher Site  Google Scholar
 F. Reif, Fundamentals of Statistical and Thermal Physics, Waveland Press, 2009.
 K. Wu, X. Li, C. Wang et al., “Apparent permeability for gas flow in shale reservoirs coupling effects of gas diffusion and desorption,” in Proceedings of the Unconventional Resources Technology Conference (URTEC '14), August 2014. View at: Publisher Site  Google Scholar
 K. R. Alnoaimi and A. R. Kovscek, “Experimental and numerical analysis of gas transport in shale including the role of sorption,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 3527–3542, Society of Petroleum Engineers, October 2013. View at: Google Scholar
 A. D. Bangham, “Membrane models with phospholipids,” Progress in Biophysics and Molecular Biology, vol. 18, pp. 29–95, 1968. View at: Publisher Site  Google Scholar
 I. Langmuir and V. J. Schaefer, “The effect of dissolved salts on insoluble monolayers,” Journal of the American Chemical Society, vol. 59, no. 11, pp. 2400–2414, 1937. View at: Publisher Site  Google Scholar
 J. P. Seidle and L. G. Huitt, “Experimental measurement of coal matrix shrinkage due to gas desorption and implications for cleat permeability increases,” in Proceedings of the International Meeting on Petroleum Engineering, pp. 575–582, November 1995. View at: Google Scholar
 O. A. Pedrosa, “Pressure transient response in stresssensitive formations,” in Proceedings of the SPE California Regional Meeting, Society of Petroleum Engineers, 1986. View at: Publisher Site  Google Scholar
 Q. Deng, R.S. Nie, Y.L. Jia, X.Y. Huang, J.M. Li, and H.K. Li, “A new analytical model for nonuniformly distributed multifractured system in shale gas reservoirs,” Journal of Natural Gas Science and Engineering, vol. 27, pp. 719–737, 2015. View at: Publisher Site  Google Scholar
 H. Mukherjee and M. J. Economides, “Parametric comparison of horizontal and vertical well performance,” SPE Formation Evaluation, vol. 6, no. 2, pp. 209–216, 1991. View at: Publisher Site  Google Scholar
 A. F. Van Everdingen and W. Hurst, “The application of the Laplace transformation to flow problems in reservoirs,” Journal of Petroleum Technology, vol. 1, no. 12, pp. 305–324, 1949. View at: Publisher Site  Google Scholar
 S. Yao, F. Zeng, H. Liu, and G. Zhao, “A semianalytical model for multistage fractured horizontal wells,” Journal of Hydrology, vol. 507, pp. 201–212, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yan Zeng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.