Study of Cyclic Fracturing in Vertical CBM Wells

Ting Li1, 2, 3, Ji-fang Wan4, *
1 Hubei Cooperative Innovation Center of Unconventional Oil and Gas, Yangtze University, Wuhan 430100, P.R China
2 Petroleum Engineering College of Yangtze University, Wuhan 430100, P.R China
3 Key Laboratory of Drilling and Production Engineering for Oil & Gas, Hubei Province, Wuhan 430100, P.R China
4 Mechanical Engineering College of Yangtze University, Jingzhou 434023, P.R China

Article Metrics

CrossRef Citations:
Total Statistics:

Full-Text HTML Views: 1178
Abstract HTML Views: 437
PDF Downloads: 3
ePub Downloads: 6
Total Views/Downloads: 1624
Unique Statistics:

Full-Text HTML Views: 555
Abstract HTML Views: 301
PDF Downloads: 3
ePub Downloads: 6
Total Views/Downloads: 865

© 2017 Li and Wan.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Mechanical Engineering College of Yangtze University, School of Mechanical Engineering, 1 Nanhuan Rd, 434023 Jingzhou City, Hubei Province, P. R. China; Tel/Fax: +8613797507143; E-mails:,



Coal-bed methane productivity of single well is very low, and has been the bottleneck of the coal-bed methane industry in China.


Although hydraulic fracturing is the only stimulation measure to develop CBM, it cannot increase production effectively, conventional fracturing method to create opening fractures should be improved. How to make good use of natural fractures, which are plentiful in CBM reservoirs, is also an important subject for hydraulic fracturing.


In this paper, the plastic deformation of coal rock is analyzed by harnessing a pseudo-Maxwell creep phenomenon, which is normally present in rock. The Kelvin-Voigt model is utilized to describe pseudo-plastic behavior of coal rock to determine pressurization and decay cyclic time for cyclic fracturing design. The mechanical requirement for shearing natural fractures is also analyzed, and shearing distance between the faces of natural fracture can be calculated by Westergaard stress function. Ultimately, the cyclic fracturing method is proposed according to theories about stress alteration and shearing of natural fractures. This method includes such periods as fracturing, pumping shut-down and so on.


A complex fracture system can be created, which consists of opened and sheared fractures, then, large SRV(stimulated reservoir volume)and flowing drainage area can be acquired. In comparison with conventional fracturing method, this new way can make full use of the characteristics of CBM reservoirs and is more suitable to CBM. This method will lead to a significant increase of CBM production, and will achieve huge economic benefits.

Keywords: Coal-bed methane, Well production, Cyclic fracturing, Complex fractures, Method, Application.


In China, vertical wells have been applied in coalbed methane (CBM) development more than 20 years, and systematic technique has been formed, but low productivity of single well has been the restrained bottleneck for CBM industry [1-5]. Although hydraulic fracturing is the main stimulation approach for CBM reservoirs, in fact, the stimulation is not well effective. The increase of production for CBM wells depends on stimulated reservoir volume (SRV) and desorption rate of methane gas. In comparison with sand reservoirs, fractures are shorter and wider in CBM reservoirs, thus, the conventional fracturing technique should be improved to attain large SRV [6, 7]. Natural fractures and bedding planes are plenty in CBM, so how to take advantage of them is an important consideration in CBM fracturing.

Network fracturing technology applied in shale gas reservoirs has been successful, which changed the basic knowledge about conventional hydraulic fracturing and brought new visions for CBM fracturing [6-9]. Multidirectional fractures can be created in CBM, and have been testified by mine-back and other monitoring means, but this fact is hard to be accepted by many technical engineers, and they usually take measures to remove these multi-fractures to avoid sand screen-out. However, in order to increase SRV and flowing drainage area in CBM, the multi-fractures or complex fracture network should be created on purpose.

How to increase flowing drainage area as large as possible by hydraulic fracturing has been the desire for petroleum engineers. For example, Kiel proposed to create dendritic fractures using injection-hesitation method, which was applied in more than 400 wells effectively, and he also applied for patent [10-12]. But this method was not accepted by many engineers as valid in the industry, and they still believe bi-wing fracture theory with existence of multi-fractures in conventional fracturing. In the meanwhile, the quest for multifractures is still continuing. Warpinski discussed altered stress fracturing to pioneer the scientific foundation for multi-direction fracturing technology [13]. Moschovidis developed this technology and applied in Mounds-injection project to create fractures shifted from one to the next direction in the same well [14]. Surjaatmadja evolved this method from theory and experiments, and applied it together with hydra-jet fracturing in reservoirs [15-18]. Presently, the successful development of shale gas in North American has also proved the value of multidirection fracturing.

In this paper, an attempt is made to further exploit multidirectional fracturing technology to its utmost potential. Taking advantage of stress induced in the process of fracturing, and of shearing the natural fractures in CBM, complex fracture system will be created in comparison with the conventional fracturing method. And design method is also presented for increasing the productivity of CBM reservoirs.


2.1. Calculation of Induced Stress

The net pressure in hydraulic fracture affects largely the around stress, especially in the direction of fracture cross section. If the fracture is longer, and the net pressure keeps almost same, which only changes at the fracture tip, and stress induced by the opening of hydraulic fracture is shown in Fig. (1).

Fig. (1). The stress induced by fracturing in the plane of Z-Y.

According to the theory of Sneddon and Elliot, stress induced by long fracture in home-generous elastic solid can be calculated as follows [19, 20].


By use of the above equations, relationship between induced stress and distance from the fracture is demonstrated in Fig. (2). As can be seen, the induced stress at direction of minimum horizontal stress increases at the most, and at direction of maximal horizontal stress, the induced stress increases moderately. The induced stress decreases rapidly as the distance is increasing. When the distance is 3 times the half of fracture height, the induced stress changes little and keeps almost the same.

Fig. (2). The induced stress around hydraulic fracture.

2.2. Shearing of Natural Fractures

There is a large amount of butt and face cleats in CBM reservoirs, so how to take advantage of these natural fractures is also an important issue. These fractures are closed initially, but if the stress condition is changed, they will slide (Fig. 3). Usually, they will slide easier than open, the shear mechanical limit will be discussed thereafter.

Fig. (3). The shearing of natural fractures.

The shearing force acted on natural fracture can be calculated as follows [21].


Some natural fractures are bonded with minerals, such as calciteor quartz, and if they are bonded well, therefore, fracturing fluid is difficult to enter into these fractures. But comparing with the coal matrix, the strength of natural fractures is lower, so Mode II or sliding fracture will be created crossed with the major fracture, which is opened during hydraulic fracturing treatment [22-25].If the force is substantial, the surfaces of natural fracture will shear like as stick slip, which is considered as dry friction motion (Fig. 4).

Fig. (4). Stick-slip between natural fractures.

The shearing asperities are generally less than one micron in size, biting into each other’s surfaces. Since the coal rock is soft and deformed easily at a certain force, so it can be regarded as plastic solid. In this case, the use of Maxwell theory to rock viscosity can be justified, and a common Herschel-Bulkley fluid model is applied [26-29].


These parameters, such as k' and n', can be determined by laboratory tests. Also, the upper stress limit and stick-slip stress can be acquired, which are used to define whether the shearing of natural fractures will happen or not.

The above equation is suitable for Newtonian motion similar to fluid, but for soft rock, the plastic motion is common, hence, the Kelvin-Voigt model is recommended for pseudo-plastic behavior [24].


The equation describes the relationship with time t. If an initial stress is exerted to the rock, the deformation will approach the elastic model response with a difference that decays exponentially as follows.


2.3. Relaxation Rate of Rock

Assume at time t1, the exerting stress disappears, the stress or deformation in rock will decay in an exponential manner [30-33].


The stresses during the pressurization and decay (relaxation) can be plotted as shown in Fig. (5).

Fig. (5). Stress during the pressurization and decay [33].

Under different stresses, the weak plane or bedding will slide and initiate shearing fractures. The shearing fracture will follow the same stress buildup and relaxation curve as shown in Fig. (5). Therefore, composite stress chart will be observed as the stick-slip phenomenon. It means that different faces of natural fractures will shear or slide similarly, thus complex fractures will be created. Also, the hesitation between fracturing or pumping can be defined, if the pressurization and decay time is determined by experiments, thus, a cyclic time interval will be scheduled.

2.4. Shearing Length

The two faces of natural fracture will slide when shearing force is more than the upper limit, then, the shearing distance can be calculated according to the Westergaard stress function as follows.


The maximum shearing distance is as follows.


2.5. Permeability of Sheared Fractures

Assume that there are not proppants in sheared fractures, so the permeability can be calculated as follows [34, 35].


If the width of sheared fracture is 10-4 m, porosity is 1%, the permeability of sheared fracture will be 8.33μm2. The permeability of coal matrix is usually at the level of 10-3μm2, obviously, the permeability of sheared fracture is far higher than the coal permeability, even the fracture width is in the scale of millimeter. When the sheared fracture is packed or sparsely by proppants, the width will be much larger, so the permeability will raise much. Therefore, the sheared fractures are important, and if they connect with each other to form fracture network, the production rate of CBM wells will rise significantly.


The objective of cyclic fracturing design is to create complex fractures and attain large SRV and drainage area. The design method should take full advantage of the characteristics, such as low elastic module and plenty natural fractures in CBM. Based on above theory about induced stress, the induced stress and change of in-situ horizontal stress can be calculated with different fracturing net pressure. In order to create fractures with different directions, fracturing net pressure will be designed and controlled during fracturing treatment. At the same time, shearing of natural fractures in CBM can be fulfilled by cyclic fracturing. Stress alteration is used to create major fractures with different extending directions, and natural fractures along with the hydraulic fractures are sheared to form complex fracture network, ultimately, the CBM reservoir will be stimulated thoroughly.

3.1. Application of Induced Stress

Fractures always extend perpendicular to the local minimum stress direction, and induces stress around nearby areas, which can be used to creat another fracture. The hydraulic fracturing period can be divided into different parts: first injection, stop and wait, then second injection and so on. Taking advantage of the stress alteration caused by the first hydraulic fracture and after a short delay, creating another fracture will be effective. If the the process is repeated several times, the stress alteration can be used to create fractures in multidirection (Fig. 6), therefore, the SRV and drainage will be larger than the bi-wing conventional fracturing method.

3.2. Shearing Natural Fractures

From above discussion, the direction of natural fractures is important for shearing, but the shearing also can be achieved when the rock around major hydraulic fractures is pressurized plastically, no matter whether the fracturing fluid enters the natural fractures or not. The coal rock is soft, and it is easy to deform in fracturing treatment, so the natural fractures will be sheared easier contrast with the hard rock. When the limit shear stress of natural fractures is determined in CBM reservoirs, the intent of shearing them can be fulfilled by adjusting the fracturing treatment pressure and schedule. Ultimately, a complex fracture system may be created with large SRV and flowing drainage area (Fig. 7).

Fig. (6). Multidirectional fractures caused by induced stress.

Fig. (7). Multidirectional fractures in reservoir plane map of multidirectional fractures.


In Qinshui CBM basin, well No. QS 40-8 was fractured as above stimulation design theory and method. The treatment zone is shown in Fig. (8), and actual fracturing curve is Fig. (9). The fracturing process is divided into three cycles. In first cycle, after the fracture is opened, extended and connected, then pumping is shut down and the pressure is unloaded. After half an hour, fracturing treatment began again, and sand concentration was raised smoothly. When pressure rose drastically, sand pumping was stopped, and repeated several times. During the treatment, micro-seismic wave was monitoring in the nearby well No. QS 40-9, and the result (Fig. 10) showed that complex or network fractures were formed. The evaluation of the micro-seismic invasion demonstrated that the width of network fracture is 75 m, and the length is 420 m. Furthermore, the production rate testified that the fracturing result was better than the conventional method. After half five months later, the CBM production was nearly 4000 m3/d, which exceeded the average production by a factor of 4 (Fig. 11).

Fig. (8). The aimed zone of well No. QS 40-8The monitoring well QS 40-9.

Fig. (9). The treatment curve of well No. QS 40-8.

Fig. (10). The result of microseismic monitoring about well No. QS 40-8.

Fig. (11). The production curve of well No. QS 40-8.


  1. The local stress around fractures will be altered during the treatment of hydraulic fracturing. Taking advantage of this stress change, multidirectional fractures can be created on purpose, which is different from conventional hydraulic fracturing technology.
  2. Shearing the plenty natural fractures in CBM reservoirs is a creative method to form complex fractures or network fracture system. This method could take full use of the characteristics of CBM reservoirs, and has been testified by underground hydraulic fracturing for extracting CBM in coal mines.
  3. By combining altering stress with shearing natural fractures, the multidirectional fracturing method could stimulate the CBM reservoirs thoroughly, and will increase the SRV and drainage area drastically, so it is suitable to CBM reservoirs.


h  = half height of fracture,m
Pn  = net pressure in fracture,MPa
v  = Poisson’s ratio of rock
σx, ∆σy, ∆σz  = the induced stress in the direction of x, y, z respectively,MPa
τyz  = induced shear stress in the plane of y-z,MPa
r  = distance between random point in reservoir and the midpoint of fracture,m
r1  = distance between random point in reservoir and the bottom point of fracture,m
r2  = distance between random point in reservoir and upper point of fracture,m
θ  = angle from the z axis to the line of random point in reservoir and midpoint of fracture,rad
θ2  = angle from the z axis to the line of random point in reservoir and bottom point of fracture,rad
θ1  = angle from the z axis to the line of random point in reservoir and upper point of fracture,rad
τ,τ0  = shearing force between faces of natural fracture,MPa
μf  = friction factor
σn  = the net compressive pressure on the natural fracture,MPa
k', n'  = the index constant of Herschel-Bulkley model
 = shear deform rate, m/s
us  = shearing distance, m
l  = fracture length, m
x  = the point position along the fracture, m
k  = Kolosov constant
G  = shear module, MPa
umax  = the maximum shearing distance, m
S(t)  = shear stress, MPa
E  = Young’s modulus of elasticity, MPa
ε(t)  = the strain, m
μa  = the absolute viscosity of the rock
Kf  = permeability of sheared fractures, × 10-3μm2
Øf  = porosity of fractures
Wf  = the width of sheared fracture, m


The authors confirm that this article content has no conflict of interest.


The work is supported by the petroleum engineering college of Yangtze University and the Yangtze Youth Fund (2016cqn27: The research on casing damage mechanism of depleted gas storage)


[1] Hongchao W, Xiaoming W, Lianggang L. Multiple fractures in the same seam in hydraulic fracturing of CBM well. Coal Geology & Exploration 2012; 40: 20-3.
[2] Liu Z, Zhao J. Correcting hole enlargement impacts on densitylogs for coalbed methane reservoirs. Open Petroleum Engineering Journal 2015; 8: 72-7.
[3] Xilin S, Ping Y. Comparative analysis on the comprehensive benefits of unconventional gas fields in Sichuan Basin between joint development model and independent model. Open Petroleum Engineering Journal 2015; 8: 126-9.
[4] Pingya L. A discussion on how to significantly improve the single-well productivity of CBM gas wells in China. Natural Gas Industry 2013; 33: 1-6.
[5] Shuang L, Wang Y, Yanyan T. Analysis on mechanism of action of reversed transmission of the pressure on adjustment of the structure of coal production and consumption in China. Open Petroleum Engineering Journal 2015; 8: 84-9.
[6] Gu H, Siebrits E, Sabourov A. Hydraulic fracture modeling with bedding plane interfacial slip. In: SPE Eastern Regional/AAPG Eastern Section Joint Meeting. 2008; pp. 9-16.
[7] Likai Z, Youjun J, Tianhong Y, Xiaoyu L. Analysis on production of coal bed methane considering the change in permeability of coal rock. Open Petroleum Engineering Journal 2016; 9: 289-98.
[8] Aimene YE, Nairn JA. Modeling multiple hydraulic fractures interacting with natural fractures using the material point method In: SPE/EAGE European Unconventional Resources Conference and Exhibition. 2014; pp. 11-6.
[9] Flottman T, Brooke-Barnett S, Trubshaw R, Naidu SK. Influence of in-situ stresses on fracture stimulations in the Surat basin Southeast Queensland In: SPE Unconventional Resources Conference and Exhibition-Asia Pacific. 2013; pp. 3-10.
[10] Ma Z, Dong X, Lin S. Focus on the development of shale gas industrial clusters in China-based on SWOT analysis. Open Petroleum Engineering Journal 2014; 7: 142-8.
[11] Kiel O M. The Kiel process - reservoirs stimulation by dendritic fracturing Journal of Petroleum Technology SPE6984 1977; 121-5.
[12] Warpinski NR, Branagan PT. Altered-stress fracturing. Journal of Petroleum Technology 1989; 41: 990-7.
[13] Li T, Wan J. Research of network fracturing method for coalbed methane reservoirs. Open Petroleum Engineering Journal 2016; 9: 247-56.
[14] Lewis BJ, Surjaatmadja JB. Method to create multiple collocated hydraulic fractures using a temporary localized change in stress anisotropy produced during an initial stimulation treatment In: SPE Annual Technical Conference and Exhibition. 2015; pp. 38-41.
[15] Surjaatmadja JB. Placing two fractures consecutively in close proximity to significantly increase revenue-to-cost ratio In: the SPE Asia Pacific Oil and Gas Conference and Exhibition. 2008; pp. 220-32.
[16] Surjaatmadja JB, Looper DW, Coffman KR. New dynamic rock test system to assist in planning and creating multidirectional fractures In: SPE /EAGE European Unconventional Conference and Exhibition. 2014; pp. 19-22.
[17] Vincent MC. Optimizing transverse fractures in liquid-rich formations In: SPE Annual Technical Conference and Exhibition. 2011; pp. 78-83.
[18] Litchfield T, Lehman J. Inter-well interference during stimulation, flow-back and production history. In: ATW on Flowback. 2014; pp. 21-8.
[19] Meyer BR, Bazan LW. A discrete fracture network model for hydraulically induced fractures: theory, parametric and case studies. In: SPE Hydraulic Fracturing Technology Conference and Exhibition. 2011; pp. 99-105.
[20] Roman A, Ahmadi G, Issen KA, Smith DH. Permeability of fractured media under confining pressure: A simplified model. Open Petroleum Engineering Journal 2012; 5: 36-41.
[21] Soliman MY, Gamadi T. Testing tight gas and unconventional formations and determination of closure pressure In: Unconventional Resources Technology Conference. 2012; pp. 38-45.
[22] Wu K, Olson JE. Mechanics analysis of interaction between hydraulic and natural fractures in shale reservoirs In: Unconventional Resources Technology Conference. 2014; pp. 47-52.
[23] Lv D, Wu X, Li Z, et al. Sedimentary mechanism in the coal and oil shale bearing succession Graben basins: the comparision between Fushun and Huangxian basins, Eastern China. Open Petroleum Engineering Journal 2012; 8: 385-91.
[24] Deng Y, Guo R, Tian Z, Xiao C, Han H, Tan W. Productivity model for shale gas reservoir with comprehensive consideration of multi-mechanisms. Open Petroleum Engineering Journal 2015; 8: 235-47.
[25] Soliman MY, East LE, Augustine JR. Fracturingdesign aimed at enhancing fracture complexity In: SPE Europe/EAGE Annual Conference and Exhibition. 2010; pp. 33-9.
[26] Urbanic T, Baig A. Enhancing recovery in shales through stimulation of pre-existing fracture networks In: SPE Hydraulic Fracturing Technology Conference. 2014; pp. 28-34.
[27] Etminan AH, Seifi A. An improved model for geostatistical simulation of fracture parameters and their effect on static and dynamic models. Open Petroleum Engineering Journal 2008; 1: 47-57.
[28] Surjaatmadja JB. Single point of initiation, dual-fracture placement for maximizing well production In: the European Formation Damage Conference. 2007; pp. 130-9.
[29] Surjaatmadja JB. Capitalizing on stick-slip pseudo-plasticity of formation improves production of dual hydraulic fracture completions. In: SPE Asia Pacific Oil and Gas Conference and Exhibition. 2007; pp. 87-93.
[30] Meng CF, Pater C. Hydraulic fracture propagation in pre-fractured natural rocks In: 44th US Rock Mechanics Symposium. 2011; pp. 123-7.
[31] Wang WW, Olson JE. Natural and hydraulic fracture interaction study based on semicircular bending experiments In: the Unconventional Resources Technology Conference. 2013; pp. 72-6.
[32] Wang T, Li H, Wang J, Jia Y, Liu L. Volumetric source model to predict productivity of fractured horizontal well in naturally fractured reservoirs. Open Petroleum Engineering Journal 2015; 8: 373-84.
[33] Puspitasari R, Gan T. Wellbore stability modelling for horizontal and multi-branch lateral wells in CBM: practical solution to better understand the uncertainty in rock strength and coal heterogeneity In: SPE Asia Pacific Oil and Gas Conference and Exhibition. 2014; pp. 36-43.
[34] Moschovidis ZA, Gordner DC, Sund GV. Disposal of oily cuttings by downhole periodic fracturing injections, Valhall, North Sea: case study and modeling concepts In: the SPE/IADC Drilling Conference. 1993; pp. 65-70.
[35] Xu L, Zhu K, Yang X. Production forecasting of coalbed methane wells based on type-2 Fuzzy logic system. Open Petroleum Engineering Journal 2016; 9: 268-78.