SciELO - Scientific Electronic Library Online

 
vol.20 issue4The geological and geodynamic condition on the formation of the Dabashan thrust nappe structure: Based on FLAC numerical modellingAccuracy Investigation of PPP Method Versus Relative Positioning Using Different Satellite Ephemerides Products Near/Under Forest Environment author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Earth Sciences Research Journal

Print version ISSN 1794-6190

Earth Sci. Res. J. vol.20 no.4 Bogotá Oct./Dec. 2016

https://doi.org/10.15446/esrj.v20n4.61633 

Articulos Originales

The pilot appraisal of acid fracturing of coalbed methane reservoir in southeast Qinshui Basin, China

Evaluación piloto de fractura ácida en depósitos de gas metano de carbón en el suroeste de la cuenca Qinshui, China

Yu Yanga  * 

Chengwei Zhangb 

Huijun Tianb 

Wangang Chenc 

Xiaodong Pengd 

Hao Zhanga 

a.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation (Chengdu University of Technology), Chengdu 610059, China

b.College of Energy, Chengdu University of Technology, Chengdu 610059, China

c.China United Coalbed Methane Co. Ltd, Beijing 100011, China

d.CNOOC Ltd. Zhanjiang, Zhanjiang 524057, China


ABSTRACT

The reserves of Coalbed Methane (CBM) in Qinshui Basin are quite promising, but the outputs from CBM wells are quite small even after massive hydraulic fracturing. Herein the fracture system with #3 and #15 coal seams in Qinshui basin was analyzed, and it was found that both of the macro-scale fractures and micro-scale fractures are filled with clay and carbonate minerals, which explains the low productivity of CBM wells after conventional hydraulic fracturing. Acid fracturing has long been an effective method for carbonate gas reservoir to improve the gas well production. However, there were few reports about the application of acid fracturing in coal bed methane field. Based on the mineral identification and acid sensitivity test, the feasibility of acid fracturing demonstrated that the acid does more help than damage to increase the permeability of coal seams in Qinshui basin. Onsite operations have shown that acid fracturing is applicable for the CBM wells in Jincheng Mining Area. It was also observed from the microseismic survey that when applying the acid fracturing treatment, the stimulated reservoir area depends on the acid volume pumped in the first stage, which is crucial to the success of the stimulation.

Keywords: Coalbed Methane; Acid Fracturing; Acid Sensitivity Test; Stimulated Reservoir Area

RESUMEN

Las reservas de gas metano de carbón (CBM, del inglés Coalbed Methane) en la cuenca Qinshui son más que prometedoras, pero la producción en los pozos es muy pequeña, incluso después de fracturas hidráulicas masivas. En este trabajo se analizaron los sistemas de fractura de las vetas de carbón #3 y #15 de la cuenca Qinshui y se encontró que tanto las fracturas a macroescala como aquellas a microescala están cubiertas con arcillas y minerales carbonatos, lo que explica la baja productividad de los pozos de gas metano de carbón después de la fractura hidráulica convencional. La fractura ácida ha sido un método efectivo en los depósitos de gas carbonato para mejorar la producción en el pozo de gas. Sin embargo, existen pocos informes sobre la aplicación de la fractura ácida en el campo del gas metano de carbón. De acuerdo con la identificación mineral y las pruebas de sensibilidad ácida, la factibilidad de la fractura ácida demostró que el ácido es reparador en el incremento de la permeabilidad en las vetas de carbón de la cuenca Qinshui. Las operaciones in situ han demostrado que la fractura ácida es aplicable para los pozos de gas metano de carbón en el área minera de Jincheng. También se observó en el sondeo microsísmico que cuando se aplica un tratamiento de fractura ácida, el área del depósito estimulada depende del volumen de ácido bombeado en primera instancia, lo que es determinante en el éxito de la estimulación.

Palabras clave: Gas metano de carbón; fractura ácida; pruebas de sensibilidad ácida; área del reservorio estimulada

1. Introduction

Current statistical data shows that the total global CBM is more than 50 billion m3 of great exploitation potential (Al-Jubori et al., 2009). In the dual-porosity system of Coalbed methane reservoir, the methane is adsorbed in the pore of the matrix (Rogers, R. E., Ramurthy, K., Rodvelt, G., & Mullen, 2007). When the fluid pressure in nature fracture system decreases to some degree, the adsorbed methane will begin to desorb from the surface of micro-pore, diffuse to nature split, and seepage to near well bore area (Schein & Mack, 2007). Hydraulic fracturing with gel should be applied to CBM reservoir as routine, which helps to stimulate the flow channel in the near wellbore zone, to make the gas production commercially viable (McCabe et al., 1999; Cramer, 2008).

The reserves of CBM in QinShui Basin of China are quite promising (Wei et al, 2007), but most of the productions of CBM wells are quite small even after massive hydraulic fracturing (Chen et al., 2009). The measured permeability of the core plugs from Qinshui Basin is minuscule, most of which being less than one mD measured by injection tests (Wang, Tang, Wei, Xu, & Leng, 2006). The clay and calcite mineral observed in the natural fractures of coal seam are to explain the low productivity of gas wells (Puri, King, & Palmer). Therefore, increasing the permeability of natural cleat is the important target for petroleum engineers who should be stimulating the coal bed methane in Qinshui Basin.

Despite the fact that acid fracturing has been proved as an important way to improve oil and gas production for carbonate reservoir with natural fractures (Hill & Schechter, 2000; Perex, D, Huidobro, E., & Avendano, 1998), it has not been applied universally in the stimulation of CBM wells. It was reported that small volume of acid spearhead ahead of the main fracturing had been implemented to clean out the cement deeper into near well bore area and reduce the breakdown pressure of a few CBM wells (Adil, Sharma, & Bhat, 2012). The acid spotting /socking with 1000 gal of inhibited 15% HCl and 10% Acetic acid mix, which react with perforation debris and remove cement invasion, has also been proved successful in India (Pooniwala, 2012). However, it is said in the literature that the coal formations are rather clean with minimum natural deposits in the cleat system, which is entirely different from those in the QinShui Basin.

To evaluate the feasibility of acid fracturing in QinShui Basin, herein a series of tests have been designed to verify that the acid is beneficial to dissolve the carbonate fillings in the coal core of both macroscale and microscale, and the pilot acid fracturing treatments have been performed to demonstrate that the HCl acid pad, the volume of which being much larger than that in literature (Adil et al., 2012; Pooniwala, 2012; Tamayo, Lee, & Taylor, 2007) can remove the carbonate cement of the natural fracture net, and lower the breakdown pressure efficiently as expected. It is learned from the microseismic survey that the stimulated reservoir volume (SRV) is dependent on the application of acid pad before the proppant laden stage, and it is suggested that larger volume of the acid pad is crucial to the success of the stimulation operation in CBM wells with carbonate filling in cleats.

2. Natural Fracture Net and Mineral Contents

2.1 Observation of macro-scale nature fracture

The main pay zones of Qinshui Basin coal bed methane are named as 3# and 15# coal seams. From the observation of working tunnels of down-hole coal mines located in Qinshui Basin, it is found that there are two groups of cleats. The strike of the first group ranges from 20 to 50°, the density of which being unevenly distributed, with the strike of the second group being between 276 and 355°. The dip angle of cleats are more than 70° and can be taken as vertical braces (Liu, Kang, Wang, & Deng, 2008). The mineral identification with HCl acid shows that most of the fillings in cleats are carbonate type, which agrees with the observation of Wei et al.

Table 1 Statistics of coal seams in Qinshui Basin 

2.2 SEM Analysis of filling in micro-scale fracture

Scanning electron microscope (SEM) is an effective way to observe the occurrence of micro fractures, such as the dimension, filling, and connection.

Since the filling in micropores and cracks is the resistance for the gas to flow, the discrimination of mineral type, current status, and the content can be the theoretical foundation for acid fracturing in the coal seam. The Qunta250 FEG Scanning Electron Microscope made by FEI Co. and INCAx-max20 Energy Dispersive Spectrometers made by Oxford Co. of the National Key Laboratory of Oil and Gas Reservoir Geology and Exploitation are applied to describe the thin sections of the coal seam. The analysis of 3# and 15# coal in Jincheng Mine are as follows:

Clay mineral

Clay mineral is the most widely distributed and abundant mineral in China's coal. It is also found in this investigation that kaolinite, illite, chlorite, mixed layer of smectite and illite are abundant in Jincheng Mining area.

Figure 1 Flake type illite and smectite mixed layer in microfracture 

Carbonate mineral

Most of the carbonate minerals in Jincheng CBM Mine area are identified as epigenetic origin (Liu et al., 2008).

Figure 2 The calcite filling in microfracture 

Sulphide Mineral

Even if the content of sulphide minerals in coal is quite small, this kind of mineral could be the origin of insoluble precipitate of acid reaction, which can plug the fracture. The common sulphide mineral in coal is identified as pyrite (Fe2S).

Figure 3 Blackberry type pyrite aggregation in microfracture 

2.3 X-Ray Diffraction Analysis of Filling in Micro-Scale Fracture

It is hard for SEM to discriminate every clay mineral sometimes. X-ray diffraction method analyzes diffraction pattern on the photographic plate by making use of X-ray diffraction formed by crystal, so as to clarify the crystalline structures and shapes of minerals, and further analysis can tell the relative content of minerals by quantity.

According to Table 2, most of the clay minerals in coal rock are illite and kaolinite, with few illite/smectite mixed layer. It is found from more than 230 reports in Jincheng Mine that the average amount of illite is 68.7%, and that of kaolinite is 24%.

Table 2 X-ray diffraction analysis results from part of coal samples 

3 Feasibility of Acid Fracturing Stimulation in Coal Bed Methane

Coal bed methane reservoir is typical of the dual-porosity media system; the methane is adsorbed in matrix pore. When coal seam pressure decreases to a certain degree, the adsorbed methane begins to desorb from the surface of matrix pores, and diffuse to macro fracture network through pores and microfracture, and flow to well bore via hydraulic fracture.

Through downhole observations, SEM and X-ray diffraction analysis, we can see fillings in coal fractures as that most of the clay minerals are kaolinite or illite, and the carbonate fillings in both macro-scale fracture and miro-scale fracture are also abundant.

Mingshou Wang et al. measured the temperature of inclusions of fillings existing in coal seam cleat systems of Qinshui Basin and suggested that the filling sequence of minerals is that pyrite precipitates first, then quartz and clay minerals, finally carbonate minerals (calcite and Fe-dolomite). So that is to explain why the carbonate minerals are filled in the center of fractures.

From the 1990s to 2010s, more than 50 CBM wells were tested with the injection-drawdown methods to measure the permeability, and the results showed that the permeability of most of the pay zones is less than 1mD. Because of the relatively low permeability of the natural fractures in 3# and 15# coal seams of Qinshui Basin, the productivity of CBM wells are quite low even if the wells have been stimulated using massive hydraulic fracturing. Therefore, by injecting some kinds of acids into the formation, to get the carbonate mineral in the center of natural fractures resolved, the conductivity of reservoir fractures could be upgraded.

The hydrochloric acid and mud acid are screened at first to simplify the feasibility of acid applied in CBM. At the early stage of the experiment, it is found that the Calcium fluoride precipitation would follow the reaction between Hydrofluoric acid and carbonate. However, the permeability of coal sample soaked in hydrochloric acid increases significantly.

Carbonate minerals react with HCl acid as follows:

3.1 Test Preparation

The coal samples used in the experiment are taken from 3# coal seam of Jincheng mine.

Figure 4 Core samples prepared from coal seam 

3.2 Porosity Change of Coal Core

There are various methods to test the porosity of coal core, and the most widely used one is using nitrogen flooding.

Because nitrogen is inert, the molecule diameter of which is only 0.265 nm and can enter micro-pore of coals, it can cover the micro-pore and fracture. Three different concentrations of 6 %, 9 %, and 15 % hydrochloric acid are used respectively. Coal samples are taken from 3# coal seam in Jincheng mine in Qinshui Basin. Coal cores are divided into three groups, each group with ten cores, to compare the porosity changes before and after treatment with acids.

Experiment Steps

When drying the coal samples before the test, caution should be taken to avoid the destruction of coal plug due to high temperature, and the temperature is suggested as 80 °C (Dabbous, Reznik, Taber, & Fulton, 1974.)

Test the porosity of coal samples of each group before acid treatment.

Have coal samples soaked 6 hours in hydrochloric acid of 6 %, 9 %, and 15 % concentration, respectively.

Test porosity coal samples after treatment.

Porosity changes are compared in Table 3.

It is found that the average porosity increase after the acid treatment. According to the Table 3, if the acid concentration increases, the average porosity will increase accordingly.

Table 3 Porosity changes of coal samples 

3.3 Evaluation of Acid sensitivity

Acid sensitivity is the permeability decline due to the reaction between an acid and some minerals (Hsi, 1984). The minerals sensitive to HCl are iron calcite, iron dolomite, hematite, pyrite, magnesium iron ore, glauconite, hydrobiotite. The minerals sensitive to HF are calcite, dolomite, calcium feldspar, zeolites, and various types of clay minerals. Since pyrite was observed in SEM, the acid sensibility of coal seam samples is measured.

In the experiment, permeability is measured before and after injecting HCl of 15% concentration. The acid-sensitive index is calculated by formula 2.

Where: K 0 is the permeability by distilled water flooding before acidification, K1 is the permeability by distilled water flooding after acidification, is the acid sensitivity index.

Table 4 Part of acid sensitivity analysis of Coal samples 

Experiments suggest that acid treatment can improve coal seam porosity, and since there are few minerals sensitive to HCl acid in coal seams, the acid sensitive index is relatively low.

3.4 Geo-Mechanic property changes of coal samples after acid reaction

Because the cement in a fracture in Jincheng Mine can be resolved in acid, the geo-mechanic properties would change accordingly (Huang, Wang, C., Wang, T., & Zhang, Z., 2015). Before this experiment, when preparing the core sample, care should be taken to avoid turbulence which may destroy the core sample because of the fragile nature of coal. The samples listed are numbered from 1 to 10. Half of the samples are soaked in acid for 2 hours before the test, while the others are tested after acid treatment. It can be seen that after acid treatment, the cores are much easier to break down.

Figure 5 Triaxial destruction experiment of coal sample. 

Table 5 Geo-mechanic parameters of coal samples before acidizing 

Table 6 Geo-mechanic parameters of part of coal samples after acidizing 

3.5 SEM observations before and after acidizing

Through SEM experiment, the changes of minerals in the micro- fractures before and after acid treatment are studied, and the result is shown in Figure 6.

Figure 6 The change of carbonate minerals in micro fractures before and after acidizing 

It can be learned from an experiment in the laboratory that, carbonate minerals in natural fractures are dissolved almost thoroughly, the porosity and its connectivity of coal sample increase significantly, and permeability is improved obviously after the acid reaction.

4. Field Application of Acid Fracturing Engineering

4.1. Case of applying acid fracturing treatment

From 2011 to 2012, ten coal bed methane wells were stimulated with acid fracturing in Qinshui Basin ([19]. Herein the fracturing treatment of well TSD1 is given as an example.

Before the operation of TSD1 well, the hydraulic of the abnormally high breakdown pressure. The initial suggestion was to inject HCl for TSD1 well to decrease the breakdown pressure of 3# coal seam, whose measured depth was from 1031.1mracturing stimulation oftwo nearby wells, TSD5, and TDD3 had been proved as a failure because of o to 1037.1m. The operation was carried out on November 22, 2012. The total volume of proppant laden fluid was 452m3, and the pump rate was 6.5m3/min. Before the proppant laden stage, 50 m3 volume of HCl acid was performed as a pad to not only clear the debris of perforation but also remove the carbonate cement of the natural fracture net. The operation curve is shown in Figure 7.

Figure 7 Fracturing curve of well TSD1 

The operation of TS06 CBM well was compared, which was stimulated with guar based polymers, to show the advantage of acid fracturing. The target zone of TS06 well was also 3# coal seam, and the perforation interval depth was 691.6m-698.2m. Instead of the acid pad before the proppant-laden stage, 2% KCl water was pumped as a pad to break down the pay zone for the TS06 CBM well.

Figure 8 Fracturing curve of well TS06 

It can be learned from the comparison of the curves from both well TSD1 and well TS06 that, since the HCl acid was applied in the pad stage of fracturing for TSD1 well, the peak pressure in proppant-laden stage was 1.5 MPa lower than that in the pad stage, which is the evidence that the HCl acid must have effectively resolved the fillings in natural fracture system. While in the proppant-laden stage of well TS06, the peak fracturing pressure was almost 2MPa higher than that in pad stage, illustrating that the permeability of natural fracture system in well TS06 has not been improved effectively with KCl pad.

4.2. Acid fracturing monitoring analysis

To monitor the acid fracturing situation of TSD1 well, Pinnacle Company was authorized to conduct micro-seismic monitoring during the fracturing operation. Ten Three-Component Geophones, combined with the optical fiber transmission technique and mathematical analysis, were installed in monitoring well TS86 which is about 300 meters away as shown in Figure 9 to measure the position and length of the fracture.

Figure 9 Location map of TS86 and TSD1 well 

After the acid fracturing pad, it is shown in Figure 10 that the injected acid has filtrated from the hydraulic fracture to the natural cleat net, the stimulated area being approximately 140m wide and 150m long.

Figure 10 Microseismic monitoring curve at the end of the pad stage 

The variation of downhole microseismic events with time is showed in the Figure 11. The Bluepoints labeled as the early stage is representing the SRV after the acid pad stage, and red points labeled as the late stage is representing the SRV after the proppant laden stage. It is also learned that the extension of stimulated reservoir volume has almost stopped after the acid pad stage, which is showed as Figure 11. The total SRV depends on the scale of the early acid pad operation.

Figure 11 Downhole microseismic events of Well TSD1 in the early stage and the final stage 

5. Conclusions

1. Fillings in 3# coal seam, or the primary production zone of Qinshui Basin, is predominated by carbonate minerals with clay minerals. Through lab tests, it is learned that the carbonate minerals in the center of natural cleat system of coal seam can explain the feasibility of acid fracturing stimulation. If the filling in natural cleat net is clay or quartz with less carbonate, the acid will not do any good to the fracturing stimulation and should be avoided.

2. The hydraulic fracturing can create a high conductivity fracture initiated from the wellbore, but the productivity of well would be quite low if the conductivity of nature fracture has not been increased. Acid fracturing pad in CBM reservoir can increase the permeability of natural fracture net with carbonate filling, and improve the production.

3. It is learned from the microseismic survey that when fracturing the CBM well, the stimulated reservoir volume is dependent on the applied volume of the acid pad, and it is suggested that larger acid fracturing pad is crucial to the success of the whole operation.

4. Massive acid fracturing treatment may contribute to some kinds of contamination in the reservoir, changing the properties of the formation. Thus a series of tests as illustrated herein should be designed before onsite treatment to avoid these damages. It is verified that careful preparation can ensure the success application of large volume of the acid pad in stimulation treatment.

Acknowledgements

This paper is supported by the National Science and Technology Major Project of China (2011ZX05014-004-004HX).

References

Adil, F., Sharma, A., & Bhat, S. (2012). Hydraulic Fracturing of CBM Wells in India Using a Unique Fracturing-Service Technology-Operational and Technological Lessons. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/153144-MS. [ Links ]

Al-Jubori, A., Johnston, S., Boyer, C., Lambert, S. W., Bustos, O. A., Pashin, J. C., & Wray, A. (2009). Coalbed Methane: Clean Energy for the World. Oilfield Review, 21(2): 4-13. [ Links ]

Chen, Z. H., Wang, Y. B., Yang, J. S., Wang, X., Chen, Y., & Zhao, Q. (2009). Influencing factors on coal-bed methane production of single well: A case of Fanzhuang Block in the south part of Qinshui Basin. Acta Petrolei Sinica, 30(3), 409-416. [ Links ]

Cramer, D.D. 2008. Stimulating Unconventional Reservoirs: Lessons Learned, Successful Practices, Areas for Improvement. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/114172-MS. [ Links ]

Dabbous, M. K., Reznik, A. A., Taber, J. J., & Fulton, P. F. (1974). The Permeability of Coal to Gas and Water. Society of Petroleum Engineers Journal, 14(6), 563-572. DOI: http://dx.doi.org/10.2118/4711-A. [ Links ]

Hill, A. D., & Schechter, R. S. (2000). Fundamentals of Acid Stimulation. In Economides, M. J. & Nolte, K. G. (Eds), Reservoir Simulation, third. edition. John Wiley & Sons Ltd, England. [ Links ]

Hsi, C. D. (1984). Evaluation of Clay Control Additives for Matrix Acidizing Operations. Paper presented at SPE Annual Technical Conference and Exhibition, Houston, Texas, 16-19 September. DOI: 10.2118/13086-MS. [ Links ]

Huang, X. H., Wang, C. M., Wang, T. Z., & Zhang, Z. (2015). Quantification of Geological Strength Index Based on Discontinuity Volume Desity of Rock Masses. International Journal of Heat and Technology, 4(33): 255-261. DOI: 10.18280/ijht.330434. [ Links ]

Liu, H. L., Kang, Y. S., Wang, F., & Deng, Z. (2008). Coal Cleat System Characteristics and Formation Mechanisms in the Qinshui Basin. Acta Geologica Sinica, 82(10): 1377-1381. [ Links ]

Ma, Z. W., & Chen, C. (2015). Research on the Overlapping Rights of Coalbed Methane in China. Environmental and Earth Sciences Research Journal, 2(1): 21-26. DOI: http://dx.doi.org/10.18280/eesrj.020105. [ Links ]

McCabe, M. A., Robert, L. M., Blauch, M. E., Terracina, J. M., Lehman, L. V., & Bowles, B. (1999). Investigation of a New Fracturing Fluid and Conductivity Enhancement Technology on Coalbed Methane Production. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/52193-MS. [ Links ]

Perex, D, Huidobro, E., & Avendano, J. (1998). Applications of Acid Fracturing Technique to improve Gas Production in Naturally Fractured Carbonate Formations, Veracruz Field, Mexico. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/47820-MS. [ Links ]

Pooniwala, S.A. (2012). Stimulation Unlocks Coalbed Methane: Lessons Learned in India. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/149872-MS. [ Links ]

Puri, R., King, G. E., & Palmer, I. D., (1991). Damage to Coal Permeability During Hydraulic Fracturing. Society of Petroleum Engineers. DOI: http://dx.doi.org/10.2118/21813-MS. [ Links ]

Rogers, R. E., Ramurthy, K., Rodvelt, G., & Mullen, M. (2007). Coalbed Methane: Principles and Practices. Oktibbeha Publishing Company. [ Links ]

Schein, G. W., & Mack, D. J. (2007). Unconventional Gas. In: Economides, M. J. & Martin, T. (Eds), Modern Fracturing: Enhancing Natural Gas Production. Energy Tribune Publishing Inc., Houston, 387-398. [ Links ]

Tamayo, H. C., Lee, K. J., & Taylor, R. S. (2007). Enhanced Aqueous Fracturing Fluid Recovery From Tight Gas Formations: Foamed CO2 Pre-Pad Fracturing Fluid and More Effective Surfactant Systems. Journal of Canadian Petroleum Technology, 47(10). DOI: http://dx.doi.org/10.2118/08-10-33. [ Links ]

Wang, M. S., Tang, D. Z., Wei, Y. P., Xu, W., & Leng, X. (2006). Reservoir Characteristics and Enrichment Mechanism of the Coal-Bed Gas in the North of Qinshui Basin. Petroleum Geology & Experiment, 28(5), 440-444. [ Links ]

Wei, C., Qin, Y., Wang, G. G. X., Fu, X., Bo, J., & Zhang, Z. (2007). Simulation study on evolution of coalbed methane reservoir in Qinshui basin, China. International Journal of Coal Geology, 72(1): 53-69. DOI: http://dx.doi.org/10.1016/jxoaL2006.12.001. [ Links ]

How to cite item: Yang, Y., Zhang, C., Tian, H., Chen, W., Peng, X., Zhang, H. (2016). The pilot appraisal of acid fracturing of coalbed methane reservoir in southeast Qinshui Basin, China. Earth Sciences Research Journal, 20(4), C1-C6. doi: http://dx.doi.oig/10.15446/esrj.v20n4.61633

Received: September 16, 2016; Accepted: December 30, 2016

* Corresponding author: yangyu@cdut.cn

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License