Abstract
Keywords
Introduction
As a clean energy source and an important strategic complement to conventional hydrocarbon resources, coalbed methane (CBM) has achieved large-scale commercial development in the United States, Australia, and Canada (Al-Jabouri et al., 2009; American Association of Petroleum Geologists, 2015; Moore, 2012; Islam, 2015). China has abundant CBM resources (Yun et al., 2012). Since the beginning of this century, with major investments from the state and related enterprises, the CBM industry in China has developed rapidly. In 2000, the number of CBM wells drilled in China was less than 500, and it increased to 12,574 by the end of 2012 (Ye et al., 2013). However, China still faces some urgent problems in the development of CBM, for example, poor economic benefits due to low average productivity of individual wells (Fu, 2015; Zhu et al., 2015).
Many literatures have reported the CBM production characteristics and its influencing factors in the basin-scale. Kaiser et al. (1994) analyzed and contrasted the CBM gas systems, resources and production between the Sand Wash Basin and the San Juan Basin. Based on statistical analysis, Ayers (2002) revealed the key factors determining CBM resources and output in the Powder River Basin and the San Juan Basin. Su et al. (2005) concluded eight main parameters affecting CBM production in the Qinshui Basin. Qin et al. (2008) conducted a comprehensive analysis of production influencing factors of CBM wells in the Fuxin Basin in China. Gentzis et al. (2008) dissected the geological and engineering factors affecting CBM well performance in Alberta.
With the continuous development of CBM in China, however, it is found that CBM productivity varies significantly in area-scale or even in wells adjacent to (Qin et al., 2012; Ye and Lu, 2016). The research results from the basin-scale have important guiding significance in the exploration period of CBM, but it is becoming more and more difficult to meet the actual demand in the production stage. Therefore, detailed studies about CBM productivity on finer scale (i.e., block-scale or area-scale) are urgently needed. In recent years, minority researchers have conducted fine studies about CBM production and the main influencing factors in block-scale (Lv et al., 2012; Tao et al., 2014; Zhao et al., 2015). However, these studies predominantly focused on early (primarily 2–4 years) production characteristics and its influencing factors, while the producibility performance of CBM wells in different production stages vary greatly. In this study, a further detailed study of CBM producibility was carried out, using the Hancheng pilot test area (HPTA), a precursor trial area in Hancheng block with mature, well-characterized CBM reservoirs and long-term production data (primarily 7–8 years, which can reflect the regional production potential more accurately). The paper begins with a review of the geological setting of the studied area, and continues with a detailed analysis of CBM productivity subarea in HPTA, and then a spectrum of geological and engineering factors was extracted and discussed, followed by a quantitative evaluation of these affecting factors based on the rough set theory. The paper concludes with a comprehensive analysis of the key factors on long-term production performance in HPTA and the implications that are helpful in optimizing the reservoir in proven and emerging CBM plays around the world.
Geological setting
The Hancheng block is located in the Southeastern Ordos Basin, China (Figure 1(b)), with an area of 1120 km2 (Chen et al., 2006). The total CBM reserves of the Hancheng block is 1.7 × 1012 m3 (Ma and Yin, 2002). HPTA is located in the northeast of the Hancheng Block. The basic structural form in HPTA is a monoclinic, inclining to the Southwest, with some small fold growth (Figure 1(c)). There are several small-scale and high-angle faults in the study area, which could cut through the surface (Wang, 2002; Zhang, 2008).

(a) Location map of the Ordos Basin. (b) Map of the Ordos Basin and the Hancheng block, and (c) Structural map of the coal seams in Hancheng pilot test area.
The Permian strata bearing with coal spread widely in the study area. The Taiyuan and Shanxi Formations in the lower Permian System are the two main coal-bearing strata (Sun et al., 2017), where the No. 3, No. 5, and No. 11 coal seams are the main targets for CBM development (Figure 2). The thickness of the three coal seams is 1–3.8, 1–6, and 2∼10.8 m, respectively. The coals in HPTA are mainly low volatile bituminous (

Stratigraphic column in the Hancheng block.
Data and methods
Data
The drilling data, logging data, core samples test data, fracturing data, and production data were collected from 70 CBM wells in HPTA, which belongs to the PetroChina Coalbed Methane Company Limited. These data were used to analyze the production characteristics and its key controlling factors of CBM wells.
A spectrum of structural, sedimentologic, geothermal, reservoir petrologic, hydrogeological, and engineering variables have been understood to affect CBM production (Ayers, 2002; Gentzis et al., 2008; Kaiser et al., 1994; Lv et al., 2012; Qin et al., 2008; Su et al., 2005; Tao et al., 2014; Zhao et al., 2015). In this study, the processes of influencing factors selection are as follows: (1) Listing all factors extensively that have potential impact on gas production. In HPTA, these factors include the structural conditions, the depositional setting, the hydrological conditions, the coal thickness, the gas content, the depth, the coal rank, the initial reservoir pressure, the critical desorption pressure, the porosity and permeability, the well selection, the well completion, the hydraulic fracturing, and dynamic liquid level drop rate. (2) Excluding the unavailable and obviously invalid parameters. In HPTA, the porosity and permeability tests are unavailable for most of 70 CBM wells, so they have to be excluded. Furthermore, the collected data show that the depositional and hydrological conditions, the coal rank change very little, and the well selection as well as the well completion are the same. Subsequently, nine geological and engineering factors (structural curvature of the coal seams (SCCS), effect of faults, burial depth, thickness, gas content, critical reservoir ratio (CRR), volume of the fracturing liquids per meter (VFLPM), volume of the fracturing sand per meter (VFSPM), and dynamic liquid level drop rate (DLLDR)) were analyzed in this study.
The rough set and its application
The rough set theory is a new mathematical method to describe the vague, imprecise, and incomplete information system (Pawlak, 1982, 1992; Zhang et al., 2001). Based on the classical set theory, the rough set theory defines an information system:
The CBM production information system (examples of 20 CBM wells in all 70 CBM wells in HPTA).
CBM: Coalbed methane; CRR: critical reservoir ratio; DLLDR: dynamic liquid level drop rate; DWF: distance between the well and the fault; HPTA: Hancheng pilot test area; SCCS: structure curvature of the coal seams; VFLPM: volume of the fracturing liquids per meter; VFSPM: volume of the fracturing sand per meter.
Each condition attribute
The result from the rough set theory analysis.
CRR: critical reservoir ratio; DLLDR: dynamic liquid level drop rate; DWF: distance between the well and the fault; SCCS: structure curvature of the coal seams; VFLPM: volume of the fracturing liquids per meter; VFSPM: volume of the fracturing sand per meter.
Results and discussions
CBM productivity subarea in HPTA
The CBM production process is thought of as having three stages: (1) dewatering stage, (2) stable stage, and (3) decline stage, and these different stages have different production characteristics (Ayers, 2002; Clarkson, 2013). In more than 7 years of production, CBM wells have experienced a long period of the stable stage, and thus the production data can reflect the production potential more accurately.
The production performance of CBM wells in HPTA is obviously different. Of the 70 wells, 21 wells produce more than 1000 m3/day, 33 wells less than 1000 m3/day, while the other 16 wells cannot produce gas. Additionally, the wells can be divided into three categories: type I (≥1000 m3/d), type II (<1000 and >0 m3/d), type III (no gas production). Furthermore, HPTA can be partitioned into four zones (Figure 3): zone I, zone II, zone III1, and zone III2. Zone I is distributed in the central and eastern regions of HPTA, all wells in zone I have gas production rate exceeding 1000 m3/d, while the water production rate is between 1.1 and 5.3 m3/d. Zone II, which surrounds zone I, is an approximate annular region where majority of wells are stripper wells with an average gas production of 541 m3/d, while the water production rate is between 1.4 and 18.2 m3/d. Zone III1, located around the northern faults, has no gas production. But the water production is relatively large, varying from 9.5 to 29.1 m3/d. Surrounding the southern small-scale fault, the wells in Zone III2 also have a large water production, ranging from 4.5 to 24.9 m3/d.

CBM productivity subarea in Hancheng pilot test area.
Geological factors affecting CBM well productivity
Structural curvature of the coal seams
Curvature is a parameter that reflects the curving degree of a curve or surface. In the viewpoint of mathematical mechanics, the distribution of fractures, which is affected by the stress field, can be predicted through the curvature of the structural planes (Bergbauer and Pollard, 2003; Lisle, 1994; Mandujano et al., 2005; Ouahed et al., 2005; Watkins et al., 2015). The larger the absolute value of the structural curvature, the greater the curving degree. Among all regions where the absolute value of structural curvature is large, the positive curvature region above the neutral plane and the negative curvature region under the neutral plane are under the environment of extension, and thus fractures generate and expand more easily in these regions (Figure 4; Shen et al., 2010; Suo et al., 2012). For the coal reservoir, the influence of the fracture (or cleat) system on permeability is significant (Close, 1993; Levine, 1996), so it is reliable to predict coal seam permeability by using the structure curvature method. In HPTA, most of the production wells lack well testing data, and thus SCCS (instead of permeability) was used to analyze the influence on gas production.

Formation deformation and fracture dominant position (left), sketch of relation between fold and stress strain (right).
According to the calculation method in Appendix 2, SCCS was calculated and the structural curvature map of the coal seams in HPTA was plotted (Figure 5). SCCS in the study area ranges from −0.001635 to 0.0011 m−1. A prominent area where the gas production rate is more than 1000 m3/d in the central of Figure 3 corresponds, in part, with a natural anomaly in Figure 5 where SCCS is a negative number with larger absolute value. The possible reason for this phenomenon is that the natural anomaly is also the region of extensional deformation, so fractures (or cleat) generate and expand easier in coal seams, guaranteeing larger permeability, and thus contributing to a better production performance.

Structural curvature map of the coal seams in Hancheng pilot test area.
Fault
Faults have an important influence on hydrocarbon accumulation (Markowski, 1998; Xu et al., 2012). During the development of CBM, the impact of fault is clearly worth considering (Gentzis et al., 2008; Sang et al., 2009; Tao et al., 2014). In HPTA, the CBM wells in zone III1 and zone III2 always produce a large amount of water and no gas (Figure 3). The following reasons may account for this phenomenon: (1) as a channel, the fault connects the coal seams with the other aquifers. Water analysis data show that the Cl− content of the coalbed water in the study area is generally 800–1300 mg/L, while the Cl− content in the water produced from wells near faults is less than 500 mg/L, the lowest is only approximately 300 mg/L, which indicates that the shallow groundwater flows to the wells through faults (Shao et al., 2014). During the dewatering stage, the depressurization of the coal reservoir in the fault-affecting regions is difficult to realize, thus, the gas could not desorb from the coal matrix, let alone move by Darcy flow to the wellbore. (2) Water is relatively active in the fault area, and thus, gas in the coal seam escapes easily, resulting in a lower gas content, which also has a negative impact on gas production.
Depth
The burial depth of the coal seams (i.e., the middle buried depth) in HPTA ranges from 349 to 771 m. Statistics show that, among all the type I wells, 68% have coal seams buried shallower than 500 m, while only 32% of the wells are deeper than 500 m. Overall, the burial depth has no significant impact on gas production, because the distribution of points in Figure 6 is relatively scattered.

Scatterplot of daily gas production and burial depth.
Thickness
Coal seam thickness is one of the important parameters for CBM resources evaluation. When other conditions keep the same, thicker coal seams are more productive (Pashin, 1991, 1997). According to the data of 70 wells, total thickness of the coal seams ranges from 2.5 to 15.9 m. Figure 7 show that the thickness of coal seams have no obvious effect on gas production. But it is to be observed that the wells with coal seams thicker than 10 m do not show high yield. This phenomenon, which also exists in the CBM exploitation of Qinshui Basin in China (Lv et al., 2012), is not consistent with the general perception. Compared to a single thin coal seam, the overlay of multiple coal seams with larger total thickness, enhances vertical reservoir heterogeneity (Jin et al., 2004), leading to a serious interlayer interference problem, and thus restricting the production potential of CBM wells (Huang et al., 2014; Qin et al., 2014; Yao et al., 2014).

Scatterplot of daily gas production and thickness.
Gas content
The gas-bearing volume (gas content) is a material base and essential condition of desired gas production (Moore, 2012). The coalbed gas content in HPTA varies from 2 to 14 m3/t, and are much lower than those in Qinshui Basin (22–25 m3/t; Lv et al., 2012). Figure 8 shows a positive correlation between the gas content and gas production. All of the type I wells have coalbed gas content greater than 5 m3/t. Generally, in a certain productive section, higher gas content makes the CBM resources more enriched, which ensures high gas production.

Scatterplot of daily gas production and gas content.
Critical reservoir ratio
Gas saturation is one of the key sources of coal reservoir heterogeneity (Bustin, 1997; Cui and Bustin, 2005; Pashin, 2010). To obtain gas saturation, the desorption testing determining gas content and the adsorption isotherm testing determining gas capacity should be taken first. However, most of the production wells in HPTA do not have the above testing data, hence, another related parameter—the critical desorption pressure to original formation pressure ratio, namely the critical reservoir ratio (CRR)—is used to analyze its influence on gas production. Figure 9 shows that CRR ranges from 0 to 0.79 based on the 70 data points, with 0 representing that wells do not produce gas. There is a positive correlation between CRR and gas production rates. Among all the CBM wells with CRR less than 0.4, none of them have a gas production rate greater than 1000 m3/d. However, 62.8% of the CBM wells with CRR larger than 0.4 have a gas production rate greater than 1000 m3/d.

Scatterplot of daily gas production and critical reservoir ratio.
Engineering factors affecting CBM well productivity
Hydraulic fracturing
The coal reservoir in the Hancheng Block is characterized as low permeability (Yao et al., 2009). To gain commercial gas flow, reservoir stimulation is necessary (Johnson et al., 2002). Presently, hydraulic fracturing stimulation technology is widely used in HPTA. During the process of fracturing, enhancing the discharge capacity and increasing fracturing fluid volume are important measures to enlarge the length as well as the width of the fractures. Similarly, a larger quantity of sand is more conducive to expanding the scale of sand-packed fractures and their conductivity (Zhao et al., 2015). Therefore, VFLPM and VFSPM are used to characterize the fracture scale for lack of fracture monitoring data. Figures 10 and 11 are scatter maps of gas production rates and VFLPM, VFSPM. The maps show that VFLPM and VFSPM of type I wells are >80 and >7 m3/m, respectively. Overall, the effective reservoir stimulation is necessary for CBM wells in HPTA to achieve an adequate production rate.

Scatterplot of daily gas production and volume of the fracturing liquids per meter.

Scatterplot of daily gas production and volume of the fracturing sand per meter.
Dynamic liquid level drop rate during the dewatering stage with single-phase water flow
Compared with the conventional sandstone reservoir, the coal reservoir has a significant stress sensitivity (Li et al., 2013; Palmer and Mansoori, 1998; Somerton et al., 1975; Seidle et al., 1992). During the CBM production process (i.e., dewatering stage, stable stage, and decline stage), with the reduction of reservoir pressure, the damage of stress sensitivity to coal reservoir is inevitable, especially in the dewatering stage with single-phase water flow (Clarkson, 2013; Moore, 2012; Tang et al., 2015).
For the CBM wells in HPTA, the bottom-hole flowing pressure, which controls the depressurization of the coal reservoir, is controlled by the dynamic liquid level in the dewatering stage with a single-phase water flow. Hence, the dynamic liquid level drop rate (DLLDR) can reflect the decline rate in the coal seam pressure. The average DLLDR of the 70 wells in the dewatering stage with a single-phase water flow were collected (Figure 12, points (0, 0) representative a well only produce water all the time). Based on current data, DLLDR of type I wells in the dewatering stage are mostly less than 4 m/d, while the wells with unreasonable DLLDR always have low gas production rates. The reason for this phenomenon is that an unreasonable DLLDR in the dewatering stage could cause severe stress sensitivity of coal reservoir, reducing permeability, and has long-term bad effects on production performance.

Scatterplot of daily gas production and the dynamic liquid level drop rate in the dewatering stage with single-phase water flow.
Comprehensive analysis
The results from the rough set theory analysis show that the degressive order of the influencing degree of these nine factors is DWF, SCCS, gas content, CRR, VFLPM, VFSPM, DLLDR, depth, and thickness. Geological factors play a major role in the control of long-term gas production, especially the effect of fault, structural curvature of the coal seams, and the gas content (
Faults play an obvious passive role for CBM accumulation and exploitation in HPTA. They create passages for gas to escape during the accumulation of CBM causing the coal seams to slide and rupture. Moreover, the coal reservoir adjacent to faults cannot be depressurized effectively, due to the continuous water supplement from other aquifer through faults. Therefore, the wells adjacent to faults, making up to 22.9% of the 70 selected wells, produced excessive water but low or no gas. Accordingly, the attribute dependency degree of DWF (0.5643) is the highest among the nine single attribute. SCCS is most likely to affect the heterogeneity of fracture development in coal seams (i.e., permeability), and thus has vital effects on gas production. Generally, gas content is the fundamental material in gas production, therefore, it is also an important factor.
In the early (primarily 2–4 years) production, CRR and hydraulic fracturing are the most significant factors to CBM production (Lv et al., 2012; Zhao et al., 2015). For long-term production performance in HPTA, however, the dependency degrees of CRR, VFLPM and VFSPM to gas production rate only rank 4–6 in the nine factors. The following reasons can account for this phenomenon: (1) CRR greatly affects the difficulty of the water drainage, because even a small degree of under-saturation can necessitate prolonged dewatering before a large reservoir volume can reach the critical desorption pressure (Lv et al., 2012). In addition, during the early production period, depressurization of coal seams mainly occurs in the fracturing affected areas near the wellbore. Consequently, the CRR and the hydraulic stimulation are the most prominent factors to early production behavior of CBM wells. (2) After a long period of production, the gas in the coal seams near the wellbore is increasingly exhausted, and depressurization extends from the fracturing affected coal seams to the original coal seams. Therefore, original permeability (i.e., the fracture (or cleat) development in original coal seams) and gas content is more important to long-term gas production than the CRR and the hydraulic fracturing. Moreover, the comprehensive dependency degree of the nine factors that influence gas production is 0.9481, which indicates that they cannot completely determine gas productivity due to other unconsidered factors and that the result obtained by the method of the rough set theory analysis are reasonable.
Conclusions
The long-term production data of 70 CBM wells in the HPTA demonstrate that production performance varies significantly in area-scale. There are certain distribution characteristics of productivity, and HPTA can be partitioned into four zones: zone I, zone II, zone III1, and zone III2. The results from the rough set theory analysis show that the degressive order of the influencing degree of these nine factors on CBM production is (i) DWF, (ii) SCCS, (iii) the gas content, (iv) CRR, (v) VFLPM, (vi) VFSPM, (vii) DLLDR, (viii) the depth, and (ix) the thickness. Geological factors, especially faults, SCCS and gas content, play a major role in controlling long-term gas production. While engineering factors, such as effective reservoir fracturing and reasonable DLLDR, are indispensable in the development of CBM. The CBM wells in HPTA have experienced long-term production, and thus the regional productivity could be reflected more accurately. The high-yield CBM wells (>1000 m3/d) should meet the following conditions: keeping away from the faults; drilling wells in predominant vadose zones (especially in the zones where SCCS are negative numbers with larger absolute values); gas contents greater than 5 m3/t; CRR larger than 0.4; VFLPM of more than 80 m3/m; VFSPM of more than 7 m3/m; DLLDR lower than 4 m/d. While the buried depth and the thickness of coal seams have no obvious effect on gas production.
