Production performance analysis for horizontal wells in composite coal bed methane reservoir

Transcription

1 Original Article Production performance analysis for horizontal wells in composite coal bed methane reservoir Energy Exploration & Exploitation 2017, Vol. 35(2) ! The Author(s) 2017 DOI: / journals.sagepub.com/home/eea Zhang Wei, Jiang Ruizhong, Xu Jianchun, Gao Yihua and Yang Yibo Abstract In this paper, the mathematical model of production performance analysis for horizontal wells in composite coal bed methane reservoir is introduced. In this model, two regions with different formation parameters are distinguished, and multiple mechanisms are considered including desorption, diffusion, and viscous flow. Then the solution of horizontal well performance analysis model is obtained by using point source function method, Laplace transform, and Stehfest algorithm comprehensively. The solution of the proposed model is verified with previous work thoroughly. The pressure transient analysis for horizontal well when producing at a constant rate is obtained and discussed. At last, different flow regimes are divided based on pressure transient analysis curves. They are early wellbore storage period, skin factor period, first radial flow regime, transition regime, second radial flow regime, transfer regime, and late pseudo-radial flow regime. The effects of related parameters such as storativity ratio, transfer coefficient, adsorption coefficient, ratio of vertical permeability to horizontal permeability, skin factor, horizontal well position in vertical direction, and inner region radius are analyzed as well according to pressure transient analysis and rate transient analysis curves. The presented work in this paper can give a better understanding of coal bed methane production performance in composite reservoir. Keywords Horizontal well, coal bed methane composite reservoir, pseudo-steady diffusion, pressure transient analysis, rate transient analysis China University of Petroleum (East China), Qingdao, China Corresponding authors: Xu Jianchun, Engineering Building B255, China University of Petroleum (East China), Qingdao , China. illeyupc@gmail.com Zhang Wei, Engineering Building B255, China University of Petroleum (East China), Qingdao , China @s.upc.edu.cn Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (

2 Wei et al. 195 Introduction Coal bed methane (CBM) which is an unconventional gas resource (Law and Curtis, 2002) has emerged as an important energy source worldwide over the last several decades and will play an important role in the future energy industry (Flores, 1998; Letcher, 2013; Moore, 2012; Rightmire, 1984). Taking the U.S. as an example, by the end of 2013, the proved reserve of CBM was 12,392 billion cubic feet which took up about 3.5% of U.S. s domestic natural gas reserves. The estimated production was 1466 billion cubic feet which accounted for about 4.75% of U.S. s domestic natural gas production (Source: dnav/ng/ng_enr_coalbed_dcu_nus_a.htm). The pore structure, gas storage, and flow mechanisms for CBM reservoir are different from that of conventional gas reservoirs. CBM reservoir is composed of the matrix, micropores, cleat system, and gas (Anbarci and Ertekin, 1991; Gong et al., 2014; Nie et al., 2012; Remner et al., 1986). CBM reservoir is always taken as a dual porosity media (Warren and Root, 1963) the macropores in the cleat system and those in the coal matrix (Pillalamarry et al., 2011; Yu-long et al., 2014). Most gas stores as adsorption gas in matrix (Gray, 1987; Pillalamarry et al., 2011). To describe the flow mechanisms in CBM reservoir, three processes are distinguished: (1) Desorption. The decrease of reservoir pressure will result in the desorption of gas from coal grains. (2) Diffusion. The gas will diffuse in micropores of matrix. (3) Viscous flow. The transport of gas in the cleat network is Darcy flow because the diameter of the cleat network is large enough (Nie et al., 2012; Yu-long et al., 2014). To analyze the flow characteristics and well performances of CBM reservoir, production performances such as pressure transient analysis (PTA) and rate transient analysis (RTA) are always used as effective tools (Doublet et al., 1994; Nie and Ding, 2010; Yu-long et al., 2014). Besides vertical wells, horizontal wells are also commonly employed in CBM development (Palmer, 2010; Palmer et al., 1993). Many scholars have studied the PTA or RTA for horizontal wells. Sarkar and Rajtar (1994) presented an analytical solution for pressure transient behavior of horizontal wells in a coal bed reservoir considering wellbore storage and skin. A semi-infinite, anisotropic CBM reservoir was assumed. Engler and Rajtar (1992) studied the pressure response for horizontal wells in CBM gas reservoirs, and both analytical solution and numerical simulation results were given in their work. The effects of the diffusion coefficient and desorption coefficient were also analyzed in detail. Clarkson et al. (2009) studied the production performance analysis workflows for both single- and two-phase flow and for both vertical well and horizontal well. King and Ertekin (1988) developed the CBM numerical simulator to study the flow characteristics of vertical fractured well or horizontal well and it can be used for two-phase flow (water and gas). Nie et al. (2012) studied horizontal well production performance analysis method in coal bed seam. Their model was solved using Laplace transform method. The rate decline and pressure transient curves were used to analyze single- and two-phase gas flow. All above works assumed the CBM reservoir to be homogenous and none of them considered composite CBM reservoirs. For composite reservoirs, Shi et al. (2012) demonstrated the equation for horizontal well pressure response in composite oil reservoirs considering the whole horizontal well located in the inner region. However, the type curves were not presented. Zhang et al. (2010) introduced a new well testing model in composite oil reservoirs and the reservoirs were assumed to be linear and have two zones. Formation properties were assumed to be different in different zones. Raghavan (2012) analyzed the pressure transient behaviors in

3 196 Energy Exploration & Exploitation 35(2) a two-region composite system using the Green s function method and the method of sources and sinks. Jiang et al. (2014) presented the rate transient behavior of fractured horizontal well in tight oil reservoir and Zhao et al. (2014) gave the performance of multistage fractured horizontal wells considering stimulated reservoir volume (SRV) for both tight gas and shale gas. The PTA and RTA were presented considering multiple effect parameters. Zhao et al. (2014) illustrated the PTA and PTA behaviors for vertical fractured wells considering multiple transport flow mechanisms including adsorption/desorption, diffusion in matrix, and Darcy flow in cleat fractures. Chen and Raghavan (1995) presented some solving methods for fractured well and horizontal well in composite oil reservoirs. As can be seen, though many scholars have done great jobs on PTA and RTA in composite system, little is related to horizontal well of CBM. This article studies the production performance behavior for single-phase CBM horizontal wells in composite CBM reservoir. The goal is to demonstrate the PTA and RTA results for the better understanding of the production performances in CBM development. The rest of the article is organized as follows: In the next section, the physical model including the transport mechanisms and the production scenario of the CBM reservoir are introduced. In the following section, the mathematical model is established, and the solution of the mathematical model is described in Model solution section. Result and discussion section will demonstrate the type curves and flow regimes based on the solution. Later, the difference among the composite CBM model, the homogenous CBM model, and the composite noncoal reservoir model will be presented. At last, sensitivity analysis will be made on several important parameters affecting CBM production. Physical model The CBM is always taken as a dual-porosity media which contains the matrix and cleat, and they generally have different properties. In our model, four flow stages are considered in CBM production process (see Figures 1 and 2). They are as follows: (1) Desorption. The pressure drop in the cleats will result in the gas desorption from the matrix grain. This phenomenon can be described using Langmuir isotherm theory (Gong et al., 2014; Nie et al., 2012; Zhao et al., 2014). (2) Diffusion in matrix micropores. After desorption from the matrix grain, gas will diffuse into the cleat. In this paper, the pseudo-steady diffusion is used. (3) Flow in cleat. Flow in the cleat network is viscous flow which obeys Figure 1. Gas storage and transport in the coal seam (modified from Remner et al. (1986)).

4 Wei et al. 197 Figure 2. Flow process in CBM reservoir. CBM: coal bed methane. Darcy s law. (4) Flow between the cleat and wellbore. The gas will flow into the wellbore from the cleat due to the decline of wellbore pressure; the flow process among different systems is illustrated in Figure 2. Because the well is always drilled in a region with best reservoir properties, the reservoir properties near the production well are always better than that far away from the well. In order to characterize the heterogeneity of the CBM reservoir, the reservoir in this work is set to be a composite system containing an inner region and an outer region. Also composite models can be applied to describe hydraulic fracturing wells with SRV; however, the proposed model in this work is only developed for characterizing the heterogeneity of CBM reservoir but not for fractured reservoirs. The permeability of inner region is assumed to be larger than that of outer region, and the horizontal well is in the inner region completely (see Figure 3). Mathematical model Mathematical model The production scenario is a radial cylindrical composite reservoir and the cylindrical coordinates are used in model description. Some assumptions are made as follows: 1. Well produces at a constant rate or at a constant wellbore pressure. The outer boundary is assumed to be infinite for the outer region and the top and bottom boundary to be closed. The inner region radius is r 1 (see Figure 3). 2. The CBM flux along the horizontal wellbore is assumed to be uniform. This assumption makes it possible to use point source method. 3. Wellbore storage phenomenon and skin factor are considered. 4. The reservoir is horizontal with uniform thickness of h and initial pressure p i. 5. For the inner region, the horizontal permeability is k h1, the vertical permeability is k v1, the compressibility C t1, and the porosity is 1 ; while for the outer region, they are k h2, k v2, C t2, 2.

5 198 Energy Exploration & Exploitation 35(2) Figure 3. Production scenario of CBM used for mathematical model. CBM: coal bed methane. For the fracture system, according to the assumptions of the physical model, the diffusivity equation for flow of CBM in fracture system can be derived by coupling mass conservation equation, equation of state, and equation of motion k h1 r p 1 þ 1 k h1 p 1 @z k v1 p 1 p 1 ¼ 1 c t1 þ p m T To linearize gas-flow equations, the pseudo-pressure approach is applied. The pseudopressure is defined as follows Z pf p ðpþ ¼ p o Z dp ð2þ ð1þ The mathematical model in the radial cylindrical coordinate system is established. Flow equation in the inner region 1 rk h1 þ r þ ¼ 1C þ p m T Flow equation in the outer region 2 rk h2 þ r þ ¼ 2C þ p m T ð3þ ð4þ The Langmuir isotherm adsorption theory is used to characterize the adsorption characteristics for CBM (Langmuir, 1916) V E ¼ V L L þ ð5þ

6 Wei et al. 199 The diffusion can be described using the pseudo-steady model as follows Outer boundary Interface ¼ DðV m V E Þ 2 ¼ z ¼ 0 or h ð7þ 2 ¼ i, r ¼1 ð8þ 1 ¼ 2, r ¼ r 1 ð9þ k ¼ k r ¼ r 1 ð10þ To describe the point source, the spherical coordinate system is utilized and the radial coordinate is R, the inner boundary condition can be described as p p 1 T sc 4 ffiffiffiffiffiffiffiffiffiffiffiffi k h1 k v1 lim R 1 "!0 p sc ¼ q point source in inner region ð11þ R¼" p p 1 T sc 4 ffiffiffiffiffiffiffiffiffiffiffiffi k h2 k v2 lim R 2 "!0 p sc ¼ q point source in outer region ð12þ R¼" The conversion relationships of all coordinate systems p r ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x 2 þ y 2 p R ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x 2 þ y 2 þ z 2 ð13þ ð14þ Dimensionless mathematical model The dimensionless variables are presented in Appendix 1. According to Appendix 1, the dimensionless equations can be obtained. For inner region, with equations (3), (5), and (6), the following equations can be D r D r D þ 2 1D r D md þð1! D ð1! Þ¼ ðv md 1D Þ D For outer region, with equations (4) to (6), the following equations can be 2D r D þ 2 2D r D r D þm1 ð! md ð17þ md ð15þ ð1! Þ¼ ðv md 2D Þ D

7 200 Energy Exploration & Exploitation 35(2) Taking the Laplace transformation with respect to t D, the following model can be obtained for equations (15) to D r D r r 1D r 2D D þ 2 1D r ¼!s 1D þð1! ÞsV md ð19þ 1D þ 2 2D r D sv md ð1! sv md ð1! ¼!s 2D þm1 ð! ÞsV md ð20þ Þ¼ V md 1D Þ¼ V md 2D Substituting equation (21) into (19) and equation (22) into (20), we D r D r r 1D r 2D D þ 2 1D r D ¼!s 1D þð1! f 1 ðsþ ¼!sþð1! Þs sð1! Þþ þ 2 2D r D ¼!s 2D þm1 ð! ð21þ ð22þ Þs sð1! Þþ 1D ð23þ ð24þ Þs sð1! Þþ 2D ð25þ f 2 ðsþ ¼!sþMð1! Þs sð1! Þþ For equations (7) to (10), the dimensionless model is 1D 2D ¼ 0, z D ¼ 0 or h D D D2 ¼ 0, r D ¼1 ð28þ D1 ¼ D2, r D ¼ r 1D ð29þ ¼ D2, r D ¼ r 1D ð30þ D M is the mobility ratio of inner region and outer region. Note that the conventional mobility ratio is defined between two different fluids and here, we use the definition to characterize the CBM mobility ability of inner region and outer region. Equations (11) and (12) and (23) to (30) are the mathematical model for composite CBM reservoir. Model solution To solve the model, the point function method is used. The detailed derivation is shown in Appendix 2. By integrating the point source regarding r 0 from L h =2toL h =2, we can obtain the pressure distribution in the reservoir for horizontal well which can be taken as the line source in the center of composite reservoir. When the horizontal well is completely in the

8 Wei et al. 201 inner region, the pressure distribution when the observation point is in the inner and outer region is shown, respectively, as follows Z 1 HDN1 ¼ 1 Z 1 HDN2 ¼ 1 DN1 0 ¼0 dr0 D r D 5 r 1D ð31þ DN2 0 ¼0 dr0 D r 1D 5 r D 5 1 ð32þ According to the above analysis, the wellbore pressure distribution can be obtained. The integral average method is applied wd ¼ 1 Z 1 HDN1 2 ¼0 dr D ð33þ 1 With the method of Duhamel s principle (Van Everdingen and Hurst, 1949), the solution considering well storage effect and skin effect can be obtained s wdðs, C D Þ ¼ wd þ S ð34þ s þ C D s 2 s wd þ S Using Stehfest numerical inversion algorithm (Stehfest, 1970), the pressure distribution in the real space can be solved and production performance type curves can be plotted. When the well is producing at a constant bottom-hole pressure, the dimensionless rate can be calculated as follows Q D ¼ 1 s 2 wd ð35þ Result and discussion The code of the framework is programed by Matlab2013a. Model verification In this paper, the CBM reservoir is assumed as a composite one which contains two regions with different formation parameters. If some parameters are set to satisfy certain conditions, the new model proposed in this paper can be converted into some other models. For instance, if f 1 ¼ f 2 ¼ f, this new model can be simplified as triple-porosity model in a homogenous reservoir; this simplified model is similar with the model proposed by Zhang et al. (2010). f can be described as equation (A-28) in the paper by Wang et al. (2013). Figure 4 depicts the bottom-hole pressure and bottom-hole pressure derivative curves. As can be seen from Figure 4, this new model can match perfectly with the existing model, which indicates that this proposed model is reliable. Type curves of horizontal wells in CBM reservoir In this section, the production performance curves are analyzed in detail. The dimensionless wellbore pseudo-pressure curve ( wd ) and the pseudo-pressure derivative curve ( 0 wd * t D/ C D ) with respect to t D /C D are plotted as a tool to recognize different flow regimes. Figures 5 and 6 illustrate the type curves of wellbore pseudo-pressure responses for horizontal well

9 202 Energy Exploration & Exploitation 35(2) Figure 4. Comparison between this study and the result proposed by Wang et al. (2013), and the basic parameters are as follows: (S ¼ 0,2,4; C D ¼ 10; L D ¼ 10; r wd ¼ ; z wd ¼ 0.5; x f ¼ 0.01; x m ¼ 0.85; x v ¼ 0.1; vf ¼ 10; mf ¼ 1; mv ¼ 0.2) Note: Definitions of parameters S, C D,L D,r wd,z wd, x f, x m, x v, vf, mf, mv can be found in the paper of Wang et al. (2013). Figure 5. Pressure-transient-type curves of horizontal well in the inner region. model with pseudo-steady diffusion in a CBM reservoir. When producing at a constant production rate, seven main transport regimes can be divided as shown in Figure 5: Regime I: The early wellbore storage period (Guo et al., 2012; Liu et al., 2015; Wang, 2014). The slope is 1 in both the pressure and pressure derivative curve in log-log coordinates.

10 Wei et al. 203 Regime II: The skin factor period (Zhang et al., 2014). This period is between wellbore storage and the early bilinear flow. Regime III: The first radial flow regime. As can be seen, this regime follows regime II. The value of pseudo-pressure derivative curves is constant which is h 1D /4. This regime may not be obvious if the thickness of the reservoir is relatively small. Regime IV: Transition regime. When the pressure spreads to the top or bottom boundaries, the first radial flow regime ends and the transition regime happens. For a homogeneous formation, this regime is also called linear flow regime (Nie et al., 2012). Regime V: Second radial flow regime. The value of pseudo-pressure derivative curves in this regime is a constant which is 0.5. Regime VI: Transfer regime. In this regime, the gas will flow from matrix to cleat because the pressure in cleat decreases significantly and the diffusion/desorption gas becomes an important supplement flowing into the cleats. There is a V shape concave that appears in pressure derivative curves in this regime, which is a typical characteristic for dual porosity media. Regime VII: Late pseudo-radial flow regime. Following regime VI, this regime starts when the gas flow from matrix to cleat and that from cleat to wellbore reach a dynamic balance state. The value in pseudo-pressure derivative curves is constant which is 0.5*M. Type curves for noncoal composite reservoir Until now, the PTA for composite CBM reservoir has been obtained. In addition, when M ¼ ¼ 1, the PTA can also be applied to homogeneous reservoir. Actually, little work has been done on PTA and RTA curves in detail for horizontal well in circular composite reservoir. In our model, if is set to be zero, the PTA or RTA curves can be calculated for the horizontal well in circular composite reservoir without considering the desorption and diffusion flow mechanisms, i.e. the PTA or RTA curves for single porosity reservoir. Figures 5 and 6 are type curves for noncoal gas without desorption and diffusion mechanisms. As can be seen, for the noncoal curves, the regime VI disappears as shown in the brown ellipse. Sensitivity analysis In this section, the influences of seven sensitive parameters are analyzed. Besides the dimensionless wellbore pseudo-pressure curve and the pseudo-pressure derivative curve, the RTAtype curves of the Q D and Q 0 D * t D/C D with respect to t D /C D are also shown. 1. Effect of storage coefficient x. Figure 7 shows the effect of storage coefficient of cleat x on PTA-type curves. According to the analysis in Type curves of horizontal wells in CBM reservoir section, there are seven flow regimes for CBM PTA curves. From Figure 7, we can see the storage coefficient mainly affects regime II, regime IV, and regime VI. A smaller storage coefficient indicates relatively smaller reserves in cleat. When producing at a constant production rate, the pressure spreads quicker than that with a bigger storage coefficient. This leads to earlier appearance of regime IV and regime VI. On the other hand, a smaller storage coefficient causes a deeper V shape concave in regime VI. Figure 8 shows the effect of storage coefficient of cleat x on RTA-type curves. When the storage coefficient increases, the production rate grows bigger, which reflects a high supply capacity for the CBM to flow from cleat to wellbore with the

11 204 Energy Exploration & Exploitation 35(2) Figure 6. Pressure-transient-type curves for noncoal reservoir of horizontal well in the inner region. Figure 7. The effect of storage coefficient on pressure-transient-type curves. constant wellbore pressure. Also, similar to the PTA curves, a smaller storage coefficient causes a deeper V shape concave in regime VI in RTA-type curves. 2. Effect of transfer coefficient. Figure 9 shows the effect of transfer coefficient on PTA-type curves. From Figure 9, we can notice that the transfer coefficient mainly affects regime VI. A smaller transition

12 Wei et al. 205 Figure 8. The effect of storage coefficient on rate-transient-type curves. Figure 9. The effect of transfer coefficient on pressure-transient-type curves. coefficient means relatively smaller diffusion coefficient for CBM. When the diffusion coefficient is small, the CBM is difficult to flow into the cleat from CBM matrix, which makes the beginning time of regime VI to be late. Figure 10 shows the effect of transfer coefficient on RTA-type curves. When the transfer coefficient increases, the production rate will increase accordingly in regime VI. This reflects a high flow capacity from matrix to cleat with the constant wellbore pressure.

13 206 Energy Exploration & Exploitation 35(2) Figure 10. The effect of transfer coefficient on rate-transient-type curves. Figure 11. The effect of adsorption coefficient on pressure-transient-type curves. 3. Effect of adsorption coefficient. Figure 11 shows the effect of adsorption coefficient on PTA-type curves. From Figure 11, we can see the adsorption coefficient mainly affects regime VI. When the adsorption coefficient is small, the CBM is difficult to flow into the cleat from CBM matrix. This results in the late appearance of regime VI and shallow V shape concave in the transfer regime. The effect of adsorption coefficient on RTA-type curves is shown in Figure 12. When the

14 Wei et al. 207 Figure 12. The effect of adsorption coefficient on rate-transient-type curves. Figure 13. The effect of ratio of vertical permeability to horizontal permeability on pressure-transienttype curves. adsorption coefficient is big, the production rate is big in regime VI. This manifests that the CBM has a high flow capacity from matrix to cleat with the constant wellbore pressure. 4. Effect of ratio of vertical permeability to horizontal permeability k v /k h. Figure 13 presents the effect of ratio of vertical permeability to horizontal permeability on PTA-type curves. From Figure 13, we can see the ratio of vertical permeability to horizontal

15 208 Energy Exploration & Exploitation 35(2) Figure 14. The effect of ratio of vertical permeability to horizontal permeability on rate-transient-type curves. Figure 15. The effect of skin factor on pressure-transient-type curves. permeability mainly affects regime II, regime III, and regime IV. When the ratio of vertical permeability to horizontal permeability is small, the horizontal permeability is much larger than the vertical permeability. This results in great difficulty for gas to flow in vertical direction and consumes more energy when the well is producing at a constant gas rate. Consequently, the value of first radial flow in pressure derivative curve becomes big, and

16 Wei et al. 209 Figure 16. The effect of skin factor on rate-transient-type curves. Figure 17. The effect of inner region radius on pressure-transient-type curves. this value equals approximately to (h 1D )/4. Figure 14 shows the effect of ratio of vertical permeability to horizontal permeability on RTA-type curves. When the vertical permeability to horizontal permeability ratio increases, the production rate will increase accordingly in the whole production lifetime. 5. Effect of skin factor. Figure 15 shows the effect of skin factor on PTA-type curves. In the proposed model, the negative skin factor which is difficult to deal with for horizontal well testing is considered.

17 210 Energy Exploration & Exploitation 35(2) Figure 18. The effect of inner region radius on rate-transient-type curves. Figure 19. The effect of horizontal well position in vertical direction on pressure-transient-type curves. This problem can be solved by bringing in effective wellbore radius (Wang et al., 2013). Skin factor mainly affects regime I and regime II as shown in Figure 15. When the skin factor increases, the resistance around wellbore will increase, engendering greater difficulty for gas stream to flow into the wellbore. Figure 16 shows the effect of skin factor on RTAtype curves. When the skin factor increases, the production rate will decrease accordingly in all time for well production lifetime.

18 Wei et al Effect of inner region radius. Figure 17 shows the effect of inner region radius on PTA-type curves. It can be seen that inner region radius mainly affects regimes V and VI. If the inner radius is small, e.g. r 1D ¼ 2, the second radial flow will last short and there will be a deep V shape in regime VI. A smaller inner region radius will decrease the duration of second radial flow. Besides, seen from Figure 18, smaller inner region radius will lead to high production rate in regimes V and VI. When M ¼ 0.5, the outer region permeability is bigger than that of inner region. A smaller inner region radius will lead to smaller flow resistance for gas after the pressure spreads to the inner region boundary, therefore resulting in higher production rate at the late time. 7. Effect of horizontal well position in vertical direction. Figure 19 shows the effect of horizontal well position in vertical direction on PTAtype curves. Three values (0.1, 0.3, 0.5) are set to study the effect of the position. Seen from Figure 19, the horizontal well position in vertical direction mainly affects regime III and regime IV. When the horizontal well position in vertical direction z wd decreases, i.e. the distance between the horizontal well and the upper/bottom boundary decreases, the first radial flow regime lasts shorter and an additional radial flow regime appears. As can be seen from the curves that when z wd ¼ 0.1, the pseudo-pressure derivative curve is a constant. When z wd increases, for example when z wd ¼ 0.3, the additional radial flow regime becomes less obvious. When z wd ¼ 0.5, the additional radial flow regime disappears. This regime is caused by the horizontal well being located away from the center position in z direction. Figure 20 shows the effect of horizontal well position in vertical direction on RTA-type curves. When the horizontal well position in vertical direction z wd decreases, i.e. the distance Figure 20. The effect of horizontal well position in vertical direction on rate-transient-type curves.

19 212 Energy Exploration & Exploitation 35(2) between the horizontal and the upper/bottom boundary decreases, the production rate will decrease in all time for well production lifetime, but the difference is not apparent. Conclusions The main goal of this paper is to present the model and results of an extended study on the production performance behaviors of horizontal wells in composite CBM reservoir. The motivation of this study comes from the fact that the CBM reservoirs are generally heterogeneous. In this paper, the performance analysis model is established considering the transport mechanisms for CBM and solved using semianalytical method. Different flow regimes are divided based on the solution of this model, and sensitivity analysis of uncertain parameters is given. According to the work and results presented, the following conclusions are obtained: 1. The mathematical model is derived for production performance analysis of horizontal wells in composite CBM reservoirs. This model can be used to analyze PTA and RTA for CBM. 2. Seven flow regimes exist in CBM production process including wellbore storage period, skin effect regime, first radial flow regime, transition regime, second radial flow regime, transfer regime, and late-pseudo radial flow regime. 3. The flow mobility ratio of inner and outer region decides the shape of the production performance curves. The value of pseudo-pressure derivative curve for late pseudoradial flow regime is 0.5* M. The value of pseudo-pressure derivative curve for second radial flow regime is There is a big difference between the coal PTA and that of noncoal reservoir: there is no transfer regime for noncoal reservoir. 5. Storage coefficient, transfer coefficient, and adsorption coefficient all have significant effects on pressure transient behavior and production performance behavior, and these parameters mainly affect the transfer regime. Bigger storage coefficient, transfer coefficient, and adsorption coefficient lead to bigger production rate. 6. Considering ratio of vertical permeability to horizontal permeability, the value of pseudopressure derivative curve for early radial flow regime equals approximately to h 1D /4. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is supported by National Natural Science Foundation of China (Grant Nos and ). References Anbarci K and Ertekin T (1991) A simplified approach for in-situ characterization of desorption properties of coal seams. SPE low permeability reservoirs symposium, Society of Petroleum Engineers, Denver, Colorado, USA, April. Chen CC and Raghavan R (1995) Supplement to SPE 28393: Modelling a fractured well in a composite reservoir. SPE Formation Evaluation 10(4):

20 Wei et al. 213 Clarkson CR, Jordan CL, Ilk D, et al. (2009) Production data analysis of fractured and horizontal CBM wells. SPE eastern regional meeting, Society of Petroleum Engineers, Charleston, West Virginia, USA, September. Doublet LE, Pande PK, McCollum TJ, et al. (1994) Decline curve analysis using type curves Analysis of oil well production data using material balance time: Application to field cases. International petroleum conference and exhibition of Mexico, Society of Petroleum Engineers, Veracruz, Mexico, October. Engler TW and Rajtar JM (1992) Pressure transient testing of horizontal wells in coalbed reservoirs. SPE rocky mountain regional meeting, Society of Petroleum Engineers, Casper, Wyoming, USA, May. Flores RM (1998) Coalbed methane: From hazard to resource. International Journal of Coal Geology 35(1): Gong B, Zhang Y, Fan Y, et al. (2014) A novel approach to model enhanced coal bed methane recovery with discrete fracture characterizations in a geochemical simulator. Journal of Petroleum Science and Engineering 124: Gray I (1987) Reservoir engineering in coal seams: Part 1 The physical process of gas storage and movement in coal seams. SPE Reservoir Engineering 2(01): Guo J, Zhang L, Wang H, et al. (2012) Pressure transient analysis for multi-stage fractured horizontal wells in shale gas reservoirs. Transport in Porous Media 93(3): Jiang R, Xu J, Sun Z, et al. (2014) Rate transient analysis for multistage fractured horizontal well in tight oil reservoirs considering stimulated reservoir volume. Mathematical Problems in Engineering 2014: King GR and Ertekin T (1988) Comparative evaluation of vertical and horizontal drainage wells for the degasification of coal seams. SPE Reservoir Engineering 3(2): Langmuir I (1916) The constitution and fundamental properties of solids and liquids. Part I. Solids. Journal of the American Chemical Society 38(11): Law BE and Curtis JB (2002) Introduction to unconventional petroleum systems. AAPG Bulletin 86(11): Letcher T (ed.) (2013) Future Energy: Improved, Sustainable and Clean Options for Our Planet. Amsterdam, The Netherlands: Elsevier. Liu M, Xiao C, Wang Y, et al. (2015) Sensitivity analysis of geometry for multi-stage fractured horizontal wells with consideration of finite-conductivity fractures in shale gas reservoirs. Journal of Natural Gas Science and Engineering 22: Moore TA (2012) Coalbed methane: A review. International Journal of Coal Geology 101: Nie RS and Ding Y (2010) Research on the nonlinear spherical percolation model with quadratic pressure gradient and its percolation characteristics. Natural Science 2(02): 98. Nie RS, Meng YF, Guo JC, et al. (2012) Modeling transient flow behavior of a horizontal well in a coal seam. International Journal of Coal Geology 92: Ozkan E (1994) New solutions for well-test-analysis problems: Part III Additional algorithms. SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, New Orleans, Louisiana, September. Ozkan E and Raghavan R (1991) New solutions for well-test-analysis problems: Part 1 Analytical considerations (includes associated papers and 29213). SPE Formation Evaluation 6(03): Palmer I (2010) Coalbed methane completions: A world view. International Journal of Coal Geology 82(3): Palmer ID, Lambert SW and Spitler JL (1993) Coalbed methane well completions and stimulations. Hydrocarbons from Coal. Boulder, USA: AAPG Publications, pp Pillalamarry M, Harpalani S and Liu S (2011) Gas diffusion behavior of coal and its impact on production from coalbed methane reservoirs. International Journal of Coal Geology 86(4):

21 214 Energy Exploration & Exploitation 35(2) Raghavan R (2012) A horizontal well in a composite system with planar interfaces. Advances in Water Resources 38: Remner DJ, Ertekin T, Sung W, et al. (1986) A parametric study of the effects of coal seam properties on gas drainage efficiency. SPE Reservoir Engineering 1(06): Rightmire CT (1984) Coalbed methane resource. American Association of Petroleum Geologists Studies in Geology 138: Sarkar PS and Rajtar JM (1994) Transient well testing of coalbed methane reservoirs with horizontal wells. Permian basin oil and gas recovery conference, Society of Petroleum Engineers, Midland, Texas, March. Shi G, Nie R, Lu J, et al. (2012) Well test model of horizontal well in 2-Zoned composite reservoir and example interpretation. Journal of Southwest Petroleum University (Science & Technology Edition) 34(5): Stehfest H (1970) Algorithm 368: Numerical inversion of Laplace transforms (D5). Communications of the ACM 13(1): Van Everdingen AF and Hurst W (1949) The application of the Laplace transformation to flow problems in reservoirs. Transactions of AIME 186(305): Wang H, Zhang L and Guo J (2013) A new rod source model for pressure transient analysis of horizontal wells with positive/negative skin in triple-porosity reservoirs. Journal of Petroleum Science and Engineering 108: Wang HT (2014) Performance of multiple fractured horizontal wells in shale gas reservoirs with consideration of multiple mechanisms. Journal of Hydrology 510: Warren JE and Root PJ (1963) The behavior of naturally fractured reservoirs. SPE Journal 3(3): Yu-long Z, Lie-Hui Z, Guo-Qing F, et al. (2014) Performance analysis of fractured wells with stimulated reservoir volume in coal seam reservoirs. Oil & Gas Science and Technology Revue d IFP Energies nouvelles 71(1): Zhang LH, Guo JJ and Liu QG (2010) A new well test model for a two-zone linear composite reservoir with varied thicknesses. Journal of Hydrodynamics, Series B 22(6): Zhang LH, Zhao YL and Liu QG (2014) Well-test-analysis and applications of source functions in a bi-zonal composite gas reservoir. Petroleum Science and Technology 32(8): Zhao YL, Zhang LH, Luo JX, et al. (2014) Performance of fractured horizontal well with stimulated reservoir volume in unconventional gas reservoir. Journal of Hydrology 512: Appendix Notation C wellbore storage factor, m 3 /Pa C t total compressibility, Pa 1 D diffusion coefficient, m 2 /s h reservoir thickness, m I k (x) the modified Bessel function of first kind k permeability, m 2 K k (x) the modified Bessel function of second kind L reference length, m L h the length of horizontal well, m M mobility ratio of inner region and outer region, fraction

22 Wei et al. 215 P pressure, Pa q point source production rate, m 3 /s Q well production rate, m 3 /s r,, z cylindrical coordinates s variable of Laplace space with respect to t D S skin factor, dimensionless t time, s T temperature, K x,y,z Cartesian coordinates V E equilibrium volumetric gas concentration, m 3 /m 3 V m volumetric gas concentration, m 3 /m 3 V L Langmuir volume, m 3 /m 3 Z compressibility factor, fraction a shape factor, m 2 adsorption coefficient for matrix, dimensionless Z diffusivity ratio, fraction transfer coefficient, dimensionless l fluid viscosity, Pa s porosity, fraction pseudo pressure, Pa/s L Langmuir pseudo pressure, Pa/s x storativity ratio, fraction Subscripts and superscripts D dimensionless i initial h horizontal direction sc standard condition v vertical direction - Laplace 1 inner region 2 outer region Appendix 1. Dimensionless definition D1 ¼ 2k h1ht sc ð i 1 Þ p sc QT D2 ¼ 2k h1ht sc ð i 2 Þ p sc QT ¼ 1 C t1 þ 2k h1hv Ei Q

23 216 Energy Exploration & Exploitation 35(2)! ¼ 1C t1 t D ¼ k h1t L 2 M ¼ k h kh 1 2 k h kh ¼ C t C t 1 r D ¼ r L sffiffiffiffiffiffi z 1D ¼ z k h1 L k v1 sffiffiffiffiffiffiffiffiffiffiffiffiffi k v1 k h2 z 2D ¼ z 1D k h1 k v2 V md ¼ V Ei V m V Ei 2 V ED ¼ V Ei V E V Ei ¼ p scqt L 2k h1 ht sc i ð L þ ¼ 2L2 h Q DV Ei sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi n ¼ f 1 þ n2 2 n ¼ h 2 1D sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi f 2 þ n2 2 h 2 2D Þ C D ¼ C 2hC t1 L 2 z wd ¼ z h Appendix 2. Solving the mathematic model The point function method is used (Jiang et al., 2014; Ozkan, 1994; Ozkan and Raghavan, 1991; Wang et al., 2013) to solve this model. As the point function may be in the inner region or in the outer region shown in Figure 21.

24 Wei et al. 217 Figure 21. Point sources. Using the same method, when the point source is in the inner region (Ozkan, 1994; Zhang et al., 2014) shown in Figure 21, the pressure distribution in the inner and outer region is " DN1 ¼ 1 X þ1 pffiffiffiffi pffiffiffiffi I k f 1r 0 2u D K k f 1rD cosð 0 Þþ 2 Xþ1 cos n z 1D cos n z0 1D h 1D h 1D Xþ1 k¼ 1 k¼ 1 þ 2 Xþ1 n¼1 DN2 ¼ 1 2u I k n r 0 D Kk ð n r D Þcosð 0 cos n z 1D h 1D cos n z0 1D h 1D Xþ1 " Xþ1 k¼ 1 X þ1 k¼ 1 k¼ 1 Þþ Xþ1 k¼ 1 c k0 I k n¼1 p rd ffiffiffiffi f 1 c kn I k ð n r D Þcosð 0 Þ cosð 0 Þ # ð36þ pffiffiffiffi d k0 K k f 2rD cosð 0 Þþ 2 Xþ1 cos n z 2D cos n z0 2D h n¼1 2D h 2D # ð37þ d kn K k n r D cos 0 ð Þ c kn ¼ C kn I k n r 0 D ð38þ d kn ¼ M D kn I k n r 0 D r 1D ð39þ C kn ¼ M nk 0 k ð nr 1D ÞK k n r 1D n K k ð n r 1D ÞK 0 k nr 1D M n I 0 k ð nr 1D ÞK k n r 1D n I k ð n r 1D ÞK 0 k ð40þ nr 1D D kn ¼ A kn ð41þ