Regular Article
Production analysis for fractured vertical well in rectangular coal reservoirs
^{1}
College of Engineering, Peking University, Room 1010, WangKeZhen Building, Peking University, NO. 5 Yihe Road, Haidian District, Beijing, PR China
^{2}
Petroleum Engineering Department, Xi'an Shiyou University, Room 212, 18 Dianzierlu, Yanta District, Xi'an, Shannxi Province, PR China
^{*} Corresponding author: dongzz@xsyu.edu.cn
Received:
15
May
2018
Accepted:
28
August
2018
As one kind of unconventional natural gas, coalbed methane is an important energy resource that is subject to active research. Gas exists in coalbed methane reservoirs in two forms: free gas and adsorbed gas. In the course of coalbed methane production, the reservoir experiences pressure decrease, desorption, diffusion, and seepage. Previous models of coalbed methane production were mainly concerned with circular boundaries. However, field tests revealed that some fractured wells possess the characteristics of rectangular boundaries. For fractured rectangular coalbed methane reservoirs, it is necessary to deal with the four boundaries with mirror image theory, which complicates calculations. In addition, the desorption and adsorption process of coalbed methane exerts a strong effect on the seepage process. Furthermore, the complexity of the rectangular coal seam embedded with the finiteconductivity fracture results in a significant computational challenge. For the first time, this paper presented a fast analytical solution for a finiteconductivity fractured vertical well model with either rectangular closed or constantpressure boundaries in the coal seam. On the basis of the Fick diffusion law and the Darcy seepage law, a mathematical model that considers diffusion in matrix and seepage within natural fractures was established. Then, we integrated the fracture conductivity function method with the hydraulic fracture model to greatly increase computational efficiency. The analytical solutions were validated against a numerical simulation. Parameter sensitivity analysis reveals that interporosity coefficient and storage coefficient, respectively, affect the appearance time and degree of desorption and diffusion. Desorption coefficient mainly describes the capacity of desorption and diffusion. Well storage coefficient, conductivity factor, and skin factor mainly affect the early stage of production. Finally, the proposed solutions were applied to field history match. The model developed is applicable to production analysis and well testing for coalbed methane reservoirs. The new proposed model extended flow mechanism of coalbed methane, and provided a better production and pressure forecast for coalbed methane reservoirs. In addition, the analytical solutions can be used to generate type curves for fractured vertical wells with finite conductivity and in the rectangular boundary, and provide a sound theoretical basis for well tests in the coal seam. The model is also applicable to other types of unconventional gas reservoirs, such as gas shales, in which the same processes are present.
© C. Li et al., published by IFP Energies nouvelles, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Nomenclature
Symbolsr : Radial radius of gas reservoir, m
Ψ : Pseudo pressure, MPa^{2}/cp
Z : Gas deviation factor, dimensionless
c_{g} : Compressibility, MPa^{−1}
V : Gas concentration, m^{3}/m^{3}
C: Volumetric gas concentration in microspores, m^{3}/m^{3}
D: Diffusion coefficient, m^{3}/s
p_{L} : Langmuir pressure, MPa
s : Image variable of Laplace transformation
R : External radius of matrix, m
L_{f} : Fracture half length, m
DPP: Dimensionless pseudo pressure
DPPD: Dimensionless pseudo pressure derivative
Subscripts and superscripts
~: Image function of Laplace transform
Intermediate variables
Dimensionless variables
1 Introduction
Coal seams are classified as unconventional gas reservoirs, together with tight gas sands and gas shales. Coalbed methane reservoirs have a dual porosity nature with macrospores made by a natural fracture network, and microspores which exist within the coal matrix. Natural fracture systems consist of two fissure systems, namely, butt and face cleats. Butt cleats are less continuous and often terminate against the face cleat with a strike angle of 75–85 degrees (Anbarci and Ertekin, 1990).
Gas flow in a coal seam involves two mechanisms: diffusion and seepage. When the coal seam pressure drops to a certain degree, the methane in the matrix translates from the adsorbed state into the free state, diffuses into the natural fracture network, and then moves through the hydraulic fracture into the wellbore. The relationship between adsorption volume and pressure can be described by the Langmuir isothermal adsorption equation. The diffusion process follows the Fick diffusion law. Flow from the natural fracture network to the hydraulic fracture is a Darcy flow. The seepage characteristics of coalbed methane make the mathematical model of seepage more complex.
Since hydraulic fracturing technology is widely and effectively utilized to enhance production efficiency in coalbed methane reservoirs, determining how to evaluate fracturing results is becoming increasingly important. During the past decades, many scholars have contributed to analyzing production and pressure performance, as it is an essential part of reservoir evaluation. Pressure performance analysis also constitutes the basis for production analysis (Mcguire and Sikora, 1960; Raghavan et al., 1972; Albinali et al., 2016; de SwaanO, 2016).
Pressure analysis of coalbed methane was commonly performed based on the dualporosity model that was introduced by Warren and Root (1962), de Swaan O (1976) and Kazemi et al. (1992). Anbarci and Ertekin (1990) proposed the Langmuir isothermal adsorption equation and Fick diffusion law to represent the adsorption and diffusion processes in seepage flow, and to analyze pressure and production performance. Gringarten et al. (1975) developed a transient pressure performance analysis and type curve analysis for fractured wells. Their work identified three basic solutions for fractured wells: infiniteconductivity solution, uniform flux solution for vertical wells, and uniform flux solution for horizontal wells.
However, as the number of segments and fracture halflength increase, the conductivity of the fracture cannot be ignored. CincoLey et al. (1978; CincoLey and Meng, 1988) developed a semianalytical solution for finiteconductivity vertical fractured wells. Aminian and Ameri (2009) proposed a numerical coalbed methane model to predict production behavior of coalbed methane reservoirs. Clarkson et al. (2009) made a production decline curve analysis for fractured wells and horizontal wells in coalbed methane reservoirs. King et al. (1986) used a numerical method to simulate transient pressure performance with finite and infiniteconductivity vertical fractures in coal seams. Nie et al. (2012) developed an analytical solution to describe its flow characteristics in a horizontal well with multistage fracturing in coalbed methane reservoirs.
Most of the above literature modeled the transport characteristics of vertical and horizontal wells with radial boundaries. However, according to field test results, some fractured wells exhibit the characteristics of the rectangular boundary. Thus, the objective of this paper was to establish a complete semianalytical model for vertical wells with finiteconductivity fractures in a rectangularboundary coal seam.
In this paper, we established and solved a new mathematical model with rectangular boundary, adsorption, desorption, diffusion, and finite conductivity fracture. Moreover, wellbore storage effects and skin factor can be added to the solution conveniently. Based on the proposed mathematical model, we optimized the calculation process to improve computational efficiency. Secondly, the model was validated with the literature results and a numerical simulator. Thirdly, the flow characteristics were analyzed. Using parameter sensitivity analysis, the type curves of a finite conductivity fracture vertical well in a rectangular coalbed methane reservoir were established in the paper. Parameters were analyzed, including interporosity coefficient, storage coefficient, desorption coefficient, fracture conductivity factor, wellbore storage coefficient, and skin factor. Finally, we applied the newly established model and calculation process to analyze the production and pressure performance of a well from the Qinshui Basin.
2 Methodology
Some important assumptions were made for mathematic model development:

the well is in the center of a rectangular reservoir that is bounded by upper and lower impermeable layers;

the inner boundary has either constant flow rate or constant pressure (Fig. 1);
Fig. 1 Reservoir model with rectangular boundary.

adsorption in the coal matrix obeys the Langmuir isotherm adsorption equation, diffusion from the matrix to the natural fracture network obeys the Fick law, flow from the natural fracture network to the hydraulic fracture and flow to the wellbore obey the Darcy law (Fig. 2);
Fig. 2 Desorption, diffusion, and Darcy flow in the matrix and fracture.

the well is intercepted by a fully penetrating fracture, and the flow rate in the fracture is different along the fracture (Fig. 3);
Fig. 3 Flow in the finiteconductivity fracture.

the reservoir is isotropic and homogeneous with a constant height, porosity, and permeability;

the initial reservoir is uniform;

gas flow takes place only through fractures;

there is no pressure loss along the wellbore;

pressure gradients are so small that the gravity effect is negligible; and

wellbore storage effect and skin effect are nonnegligible.
2.1 Model establishment and solution
Anbarci and Ertekin (1990) obtained a series of solutions for pressure transient behavior in coal seams for various outer and innerboundary conditions. However, their solutions only work for vertical, unfractured, and radial flow. In this article, to establish a sound continuity equation for coalbed methane reservoirs, we considered the flow mechanism of free gas and adsorbed gas separated from the matrix and fractures. According to mirror image theory, the seepage condition of the one well with one fault in the reservoir is completely identical to that of two wells in the reservoir. One well is the same as the original well, and the other well is a reflected duplication of the original well that appears almost identical, but is reversed in the direction perpendicular to the fault surface. Using the Laplace transformation and the mirror image theory, we can obtain the solution of the continuity equation from Ozkan and Raghavan (1988, 1991) source solution.
2.1.1 Modeling flow in the reservoir
In the material balance equation of coalbed methane reservoirs, besides free gas flow, adsorbed gas separates from the matrix and participates in the flow, which is different from conventional gas. The continuity equation is defined as equation (1):(1)
Gas seepage from matrix to natural fracture network was usually expressed by the Fick diffusion law:(2)
The Langmuir isothermal adsorption equation describes the desorption from adsorbed gas to free gas:(3)
Matrix geometry is assumed to be composed of spherical blocks with radius R. The unsteadystate sorption equation for the coal matrix can be written as:(4)
Inner boundary condition in the dimensionless form for the microspore equation is:(5)
Outer boundary condition in the dimensionless form for the microspore equation is:(6)
The solution for equation (2):(7)
Combining with equation (4), one can linearize equation (1).
It is required for to vanish as for and at for . These conditions are expressed by:(8) (9)
Normalized pseudopressure defined by Meunier et al. (1984) is:(10)
The final continuity equation in the Laplace domain is:(11)
The closed boundary or constant pressure conditions in a rectangular coordinate system are:(12) (13) (14)
Translating a radial coordinate system to a rectangular coordinate system as:(15)
Ozkan and Raghavan (1991) presented solutions for a rectangular closed and constantpressure outer boundary in the conventional reservoir as follows.
Rectangular closed outer boundary:(16)
Rectangular constantpressure outer boundary:(17)
Different from conventional gas reservoirs, coalbed seams have strong adsorption and desorption with the seepage. So, we must reflect this phenomenon in the control equation (Eq. (1)). The new parameter γ (the specific expression shown in Intermediate Variables) will reveal this specific phenomenon in the coalbed seam. Combining with it, we obtain the solutions in the coal seam as follows.
Rectangular closed outer boundary:(18)
Rectangular constantpressure outer boundary:(19)where
2.1.2 Modeling flow in the fracture
Taking equation (18) to integral along the fracture direction yields:(20)
In this equation, is a function of the location in the fracture. To integral the equation along the fracture, we use:
Using this definition, equation (20) can be translated as:(21)
Assuming that the fracture can be divided into n segments. Discretizing equation (21) with Figure 4 yields:(22)
Fig. 4 Discretizing the equation along the fracture. 
In the Laplace domain, the sum of all of the flow rates in each part is 1/s. So, it can be expressed as:(23)
In the process of discretization, we divided the fracture into n parts. We assumed that the flow rate of each part is different, but the flow rate is constant in the same part. For each part, there is a pressure equation. Combined with the flow rate normalization equation, we can obtain an equation set with n + 1 equations and n + 1 unknown numbers as . Solving this equation set, the solution of can be determined.
Duhamel’s theory can be utilized to incorporate the wellbore storage coefficient and skin factor into well response. So, the final solution with wellbore storage and skin in the Laplace domain can be written as:(24)
2.1.3 Calculation process optimization
Three factors lead to the complex calculation of finiteconductivity fractured wells in rectangular coal reservoirs: (1) solution of the rectangular reservoir from the mirror image theory, which produces mirror wells; (2) diffuse gas flowed in the coal seam will increase complexity; and (3) finiteconductivity fracture requires solving an equation set. All of these reasons eventually result in slower computing for the normal solution that is not suitable for software development. The objective was to obtain a quick calculation method to determine the pressure or production for finiteconductivity fractured wells in a rectangular coal seam.
We introduced the fracture conductivity function to translate the infiniteconductivity fracture solution to the finiteconductivity fracture. The specific equation translates as Wang et al. (2014):(25)
From Riley et al. (1991), expresses the fracture conductivity as:(26)
2.1.4 Model verification
To verify our model and solutions, we compared it with Anbarci and Ertekin’s work (1992) and a reservoir simulator the UNConventional Oil and Gas simulator (UNCONG) (Li et al., 2015). Anbarci and Ertekin (1992) set the model with an infinite reservoir boundary and infinite conductivity fracture. Corresponding to their model, we set a large enough boundary and fracture conductivity factor. UNCONG incorporates comprehensive flowing mechanisms for unconventional reservoirs (Li et al., 2015). The same as Anbarci et al.’s work, basic data used to perform the calculation are listed in Table 1.
Basic data to verify results.
Compared the normal solution with Anbarci’s analytical solution, the solution presented in this paper fits better with the numerical solution (Fig. 5). It can thus be concluded that the new method to solve the fracture model is more accurate.
Fig. 5 Comparison between Anbarci analytical and numerical solution. 
From Figure 6, we can see that using the fracture conductivity function method can improve computation speed greatly. In fact, it is 488 times faster with a closed boundary, and 599 times faster with a constantpressure boundary. For the closed boundary, the relative error was 1.25%, and for the constant pressure boundary, the relative error was 3.23%. So, this function can be used for rectangular boundaries to dramatically increase calculation speed.
Fig. 6 Comparison between normal and conductivity function calculations. 
3 Results and discussion
3.1 Flow characteristics analysis of type curves
Because of the two different sets of borders, pressure reaches the boundary twice. Under this condition, the flow state is different from square or circle boundaries. It can be divided into seven stages (Fig. 7):

Wellbore storage stage. This stage mainly reflects the effect of wellbore storage at an early time. The linear slope of both dimensionless pseudopressure and its derivative is 1.

Transition stage of wellbore storage. In this stage, the wellbore storage influence decreases, and the gas in the reservoir flows into the wellbore gradually. The pressure is influenced by the reservoir and the wellbore storage, and the curve of the dimensionless pseudopressure derivative decreases.

Linear flow stage from natural fracture to hydraulic fracture. This stage can be divided into two substages: bilinear flow and linear flow. The slope of the dimensionless pseudopressure derivative curve changes from 1/4 to 1/2. Fracture conductivity determines the length of bilinear flow. Specifically, the larger the conductivity, the longer the duration time. The dimensionless pseudopressure decreases slightly in the reservoir in this stage; the diffusion phenomenon is not obvious.

Desorption and diffusion stage. With the pressure decreasing, desorption becomes increasingly obvious. The adsorbed gas desorbed from the matrix becomes free gas, leading to concentration increase around the matrix and diffusion to the natural fracture. This stage reflects the dimensionless pseudopressure appearing in a “V” shape; its depth and appearance time are influenced by the wellbore storage coefficient, interporosity coefficient, storage coefficient, and desorption coefficient. This stage constitutes the important flow stage that reflects the characteristics of coalbed methane.

Pseudoradial flow stage. In this stage, the gas diffused from the matrix has dynamic equilibrium with the gas flowing into the wellbore. The speed of pressure decline becomes slow.

First boundary control flow stage. At this time, the pressure wave reaches the first boundary. With the constantpressure boundary, the pressure reaches stabilization quickly. However, the closed boundary has a transition, and so the dimensionless pseudopressure and derivative reflect gradually. This is because the constantpressure boundary has a sufficient energy supply to enable the pressure to reach a balance in a short time.

Second boundary control flow stage. At this time, the constantpressure boundary has already reached balance in the first boundary control flow stage, but the closed boundary must reach the pseudosteady state again.
Fig. 7 Flow stage divided for boxed coal reservoirs. 
In these seven stages, the linear flow stage mainly reflects the characteristics of fractured vertical wells and the conductivity of a hydraulic fracture. The desorption and diffusion stage mainly reflect the characteristics of an unconventional gas reservoir that has the phenomenon of desorption and diffusion, such as coalbed methane and shale gas. The two boundary control flows mainly reflect the characteristics of a rectangular reservoir. In this picture, C represents the closed boundary type, CP represents the constant pressure boundary type, IN represents the infinite boundary type.
3.2 Sensitivity analysis
Based on the definition above, traditional factors, such as permeability, porosity and compressibility, were regrouped as new factors, such as interporosity coefficient, storage coefficient, and desorption coefficient (shown in Intermediate Variables). Combined with the fracture conductivity factor, wellbore storage coefficient and skin factor, we performed a sensitivity analysis on those six factors (Tab. 2).
Basic factors for sensitivity analysis.
Interporosity coefficient λ is the direct ratio of permeability and time constant, and the inverse ratio of porosity, viscosity, fracture halflength, and gas compressibility. Figure 8a gives the effect of λ on coal seams with closed boundary conditions. It can be seen from the curves that a greater λ will lead to a smaller time of concavedown characteristic appearance. The λ slightly influences the degree of desorption and diffusion.
Fig. 8 Influence of parameters on type curve. (a) Influence of interporosity coefficient; (b) Influence of storage coefficient; (c) Influence of desorption coefficient; (d) Influence of well storage; (e) Influence of skin factor; (f) Influence of fracture conductivity. 
Storage coefficient ω is the direct ratio of porosity; it reflects the storage ability of the fracture system. It has an inverse effect compared with the conventional dual porosity medium model reflecting the matrix storage capacity. Figure 8b gives the effect of ω on coal seams with closed boundary conditions. The larger the ω, the greater the porosity of the fracture and the gas flow rate from the fracture, the smaller the gas flow rate from the matrix diffusion. Pressure maintenance due to desorption is more obviously marked as ω becomes smaller. However, the effect is inverse at an early time due to the flow occurring in the natural fractures.
Desorption coefficient σ mainly describes the capacity of desorption and diffusion. It constitutes a direct ratio of Langmuir volume, flow rate, and original pressure. Figure 8c gives the effect of σ on coal seams with closed boundary conditions. It can be seen from the curves that a greater σ will lead to more gas that the rock can adsorb and produce. In other words, a quicker change of gas concentration leads to an enhanced desorption and diffusion effect, so as to make the concavedown phenomenon become obvious and early.
Well storage coefficient mainly affects the early stage of production. Figure 8d gives the effect of well storage coefficient on coal seams with closed boundary conditions. It can be seen from the curves that a larger coefficient leads to a more obvious well storage effect and longer well storage stage. If the well storage coefficient reaches a certain value, it covers up the linear flow phase, making the flow move directly into the matrix methane diffusion phase.
Skin factor mainly affects pressure performance for wellbore storage and its transition stage, primarily influencing the end time of this stage and dimensionless pseudopressure decrease. Figure 8e gives the effect of skin factor on coal seams with closed boundary conditions. It can be seen from the curves that a larger skin factor will prolong the time of the wellbore storage stage. This is because, at the same flow rate conditions, the higher skin factor will lead to more energy consumed and pressure decreased. After this stage, the skin factor does not influence pressure performance.
Fracture conductivity factor mainly influences the transition stage of wellbore storage and the linear flow stage; it has no influence on the wellbore storage stage or after the linear flow stage. This is because those stages have low gas flow coming from the natural fracture to the hydraulic fracture. At this time, the finiteconductivity fracture can lead these gas flows to the wellbore quickly, and can be equivalent to the infiniteconductivity fracture. In addition, with the fracture conductivity factor increasing to four times of Pi (Pi = 3.14), the finiteconductivity fracture gets close to the infiniteconductivity fracture (Fig. 8f). This will guide us to not blindly improve fracture conductivity to improve output; this will not only not influence output, but will increase fracturing cost.
3.3 Production analysis
Duhamel convolution is used to construct complex solutions (such as variable rates or pressure) with simple solutions (such as constant rates or pressure). Superposition uses the theory that a rate that changes from q_{1} at time t to a new rate q_{2} is equivalent to q_{1} continuing forever, superposed, or added on (q_{2}–q_{1}) starting at time t and continuing forever (Fig. 9).
Fig. 9 The principle diagram of Duhamel convolution. 
The steps of the Duhamel convolution are as follows: (1) determine the average pressure through the material balance equation; (2) normalize the production time; (3) calculate the dimensionless pressure or rates through the analytical model and time superposition; and (4) obtain the rate or bottomhole pressure through a lookup table.
3.4 Field application
3.4.1 Basic data
The proposed model was used to analyze production data from a fractured vertical well in a coal methane reservoir located in the Qinshui Basin in Shanxi Province of China. Daily production data, flowing wellhead pressure measurements (converted to bottomhole pressure for analysis), and some reservoir parameters were collected (Tab. 3). From the seismic data analysis, the well located in the closed rectangle boundary. We used the proposed model in this paper to match the reservoir pressure and the flow rate.
Basic data of actual field well.
3.4.2 History match result
Using the model in this paper, and combining with the principle of time superposition and the dimensional result from the model, we can obtain the match result. Figure 10a showed the pressure matching. The average error was 18%. Some additional factor may have resulted in the large error at the end of the match.
Fig. 10 History match result. (a) History match of pressure; (b) History match of flow rate; (c) History match of cumulative gas production. 
Figure 10b presented the flow rate matching. The average error was 17%. The same as the pressure match, some additional factor may have resulted in the large error at the end of the match. Figure 10c showed the cumulative gas of the calculated value and field value. There was a good match result between the calculated value and field value.
4 Conclusion
In this paper, we developed a rectangularboundary model for a finiteconductivity fractured well in a coal seam. With the mirror image theory and Laplace transformation, the solutions of this model were obtained at a constant production rate or bottomhole pressure. Through the above analysis, the following conclusions can be summarized:

By comparing the solution of the simplified model in this paper with the classical analytical model and numerical simulation model, the results show that the proposed method is more accurate and efficient to solve the fractured model and has a good match with the numerical solution.

Fracture conductivity function can dramatically improve computation speed without sacrificing accuracy for a rectangular boundary model.

Compared with the radial boundary, the new rectangular model has two new flow stages: the first and second boundary control flow stages, which results from the characteristics of rectangular reservoirs. Closed and constant pressure boundaries show the difference stabilization time in pressure performance: pressure reaches stabilization quickly with the constantpressure boundary, but the closed boundary has a transition after reaching the first boundary, and reaches stabilization after reaching the second boundary.

Interporosity coefficient influences the appearance time to desorption and diffusion, and slightly influences the degree. Storage coefficient influences the degree of desorption and diffusion, and slightly influences the appearance time. Desorption coefficient mainly describes the capacity of desorption and diffusion. Well storage coefficient and skin factor mainly affect the early stage of production. Conductivity factor mainly influences the early time, after which gas flow rates from the natural fracture decrease and can flow into the wellbore quickly – the hydraulic fracture equivalent for the infinityconductivity fracture.

The solution of this paper can solve the problem of historical matching, in which the matching speed was fast and the result was good.
Acknowledgments
This article was supported by National Science and Technology Major Project of China (Project No. 2016ZX05037003002). We thank experts from Energy Innovation Software Company (Ennosoft) who provided insight and expertise that greatly assisted the research.
References
 Albinali A., Holy R., Sarak H., Ozkan E. (2016) Modeling of 1D anomalous diffusion in fractured nanoporous media, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 71, 56. [Google Scholar]
 Aminian K., Ameri S. (2009) Predicting production performance of CBM reservoirs, J. Nat. Gas Sci. Eng. 1, 1, 25–30. [CrossRef] [Google Scholar]
 Anbarci K., Ertekin T. (1990) A comprehensive study of pressure transient analysis with sorption phenomena for singlephase gas flow in coal seams, SPE Annual Technical Conference and Exhibition of the Society of Petroleum Engineers, New Orleans, Louisiana, September 23–26. [Google Scholar]
 Anbarci K., Ertekin T. (1992) Pressure transient behavior of fractured wells in coalbed reservoirs, SPE Annual Technical Conference and Exhibition, Washington, D.C., October 4–7. [Google Scholar]
 CincoLey H., Samaniego V.F., Dominguez A.N. (1978) Transient pressure behavior for a well with a finiteconductivity vertical fracture, SPEJ 18, 4, 253–264. [CrossRef] [Google Scholar]
 CincoLey H., Meng H.Z. (1988) Pressure transient analysis of wells with finite conductivity vertical fractures in double porosity reservoirs, SPE Annual Technical Conference and Exhibition, Houston, Texas, October 2–5. [Google Scholar]
 Clarkson C.R., Jordan C.L., Ilk D., et al. (2009) Production data analysis of fractured and horizontal CBM wells, SPE Eastern Regional Meeting, Charleston, West Virginia, September 23–25. [Google Scholar]
 de Swaan O., A. (1976) Analytic solutions for determining naturally fractured reservoir properties by well testing, SPEJ 16, 3, 117–122. [Google Scholar]
 de Swaan A. (2016) Pressure transients in a fractalcluster model of porous media, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 71, 9. [Google Scholar]
 Gringarten A.C., Ramey H.J. Jr, Raghavan R. (1975) Applied pressure analysis for fractured wells, J. Pet. Technol. 27, 7, 887–892. [CrossRef] [Google Scholar]
 Kazemi H., Gilman V.J., Elsharkawy A.M. (1992) Analytical and numerical solution of oil recovery from fractured reservoirs with empirical transfer functions, SPEJ 7, 2, 219–227. [Google Scholar]
 King G.R., Ertekin T., Schwerer F.C. (1986) Numerical simulation of the transient behavior of coalseam degasification wells, SPE Form. Eval. 1, 2, 165–183. [CrossRef] [Google Scholar]
 Li X., Zhang D., Li S. (2015) A multicontinuum multiple flow mechanism simulator for unconventional oil and gas recovery, J. Nat. Gas Sci. Eng. 26, 652–669. [CrossRef] [Google Scholar]
 Mcguire W.J., Sikora V.J. (1960) The effect of vertical fractures on well productivity, J Pet Technol. 12, 10, 72–74. [CrossRef] [Google Scholar]
 Meunier D.F., Kabir C.S., Wittmann M.J. (1984) Gas well test analysis: use of normalized pressure and time functions, SPE Annual Technical Conference and Exhibition, Houston, Texas, September 16–19. [Google Scholar]
 Nie R.S., Meng Y.F., Guo J.C., et al. (2012) Modeling transient flow behavior of a horizontal well in a coal seam, Int. J. Coal Geol. 922, 54–68. [CrossRef] [Google Scholar]
 Ozkan E., Raghavan R. (1988) Some new solutions to solve problems in well test analysis problems: part 1 – analytical considerations, SPE Form. Eval. 6, 3, 359–368. [Google Scholar]
 Ozkan E., Raghavan R. (1991) New solutions for welltestanalysis problems: part 2 – computational considerations and applications, SPE Form. Eval. 6, 3, 369–378. [CrossRef] [Google Scholar]
 Raghavan R., Cady G.V., Ramey H.J. (1972) Well test analysis for vertically fractured wells, J. Pet. Technol. 24, 8, 1014–1020. [CrossRef] [Google Scholar]
 Riley M.F., Brigham W.E., Horne R.N. (1991) Analytic solutions for elliptical finiteconductivity fractures, SPE Annual Technical Conference and Exhibition, Dallas, Texas, October 6–9. [Google Scholar]
 Wang X., Luo W., Hou X., et al. (2014) Transient pressure analysis of multiplefractured horizontal wells in boxed reservoirs, Petrol. Explor. Dev. 41, 1, 74–78. [CrossRef] [Google Scholar]
 Warren J.E., Root P.J. (1962) The behavior of naturally fractured reservoirs, SPEJ 3, 3, 245–255. [Google Scholar]
All Tables
All Figures
Fig. 1 Reservoir model with rectangular boundary. 

In the text 
Fig. 2 Desorption, diffusion, and Darcy flow in the matrix and fracture. 

In the text 
Fig. 3 Flow in the finiteconductivity fracture. 

In the text 
Fig. 4 Discretizing the equation along the fracture. 

In the text 
Fig. 5 Comparison between Anbarci analytical and numerical solution. 

In the text 
Fig. 6 Comparison between normal and conductivity function calculations. 

In the text 
Fig. 7 Flow stage divided for boxed coal reservoirs. 

In the text 
Fig. 8 Influence of parameters on type curve. (a) Influence of interporosity coefficient; (b) Influence of storage coefficient; (c) Influence of desorption coefficient; (d) Influence of well storage; (e) Influence of skin factor; (f) Influence of fracture conductivity. 

In the text 
Fig. 9 The principle diagram of Duhamel convolution. 

In the text 
Fig. 10 History match result. (a) History match of pressure; (b) History match of flow rate; (c) History match of cumulative gas production. 

In the text 