High production indexes and key factors influencing coalbed methane (CBM) horizontal well productivity

Similar documents
The Effect of Coal Seam Gas saturation on CBM Well Productivity-A Case Study of Central Region of Hedong Area

Study on Gas Content Critical Value for Regional Prediction in Hebi Area

The scope of fracture zone measure and high drill gas drainage technology research

Journal of Chemical and Pharmaceutical Research, 2014, 6(3): Research Article

Numerical simulation on down-hole cone bit seals

SOIL IMPROVEMENT BY VACUUM PRELOADING FOR A POWER PLANT PROJECT IN VIETNAM

Basketball field goal percentage prediction model research and application based on BP neural network

The Study on the Influence of Gust Wind on Vehicle Stability Chen Wang a, Haibo Huang b*, Shaofang Xu c

Coupling and Analysis of 981 Deep Water Semi-submersible. Drilling Platform and the Mooring System

Handlebar Bi-Directional Portable Pump on Vehicle Yingying Tang1, Yang Jiang1, Yugang Xin1 and Baida Yang1

A Correlation Study of Nadal's and Federer's Technical Indicators in Different Tennis Arenas

Gas Emission Modelling of Gate Road Development

Investigations into the identification and control of outburst risk in Australian underground coal mines

Traffic safety analysis of intersections between the residential entrance and urban road

Coal Bed Methane (CBM) Permeability Testing

Review of Gas Emission Prediction and Control Methods for Multi-seam Mining in Chinese Coal Mines

67. Sectional normalization and recognization on the PV-Diagram of reciprocating compressor

Study on the shape parameters of bulbous bow of. tuna longline fishing vessel

Study on Mixing Property of Circulation/Shear Coupling Mechanism-Double Impeller Configuration

The Application of Pedestrian Microscopic Simulation Technology in Researching the Influenced Realm around Urban Rail Transit Station

Research on Small Wind Power System Based on H-type Vertical Wind Turbine Rong-Qiang GUAN a, Jing YU b

Feasibility Analysis of China s Traffic Congestion Charge Legislation

ARTICLE. Overcome casing integrity issues in your plug-and-perf operations

The Mechanism Study of Vortex Tools Drainage Gas Recovery of Gas Well

Yuele Jia Computer Science School, Southwest Petroleum University, Chengdu , Sichuan, China

Investigation of Winning Factors of Miami Heat in NBA Playoff Season

CFD Analysis and Experiment Study of the Rotary Two-Stage Inverter Compressor with Vapor Injection

Numerical simulation of three-dimensional flow field in three-line rollers and four-line rollers compact spinning systems using finite element method

Decompression Method For Massive Compressed Files In Mobile Rich Media Applications

Elastic-Plastic Finite Element Analysis for Zhongwei--guiyang Natural Gas Pipeline Endangered by Collapse

Influence of Ventilation Tube Rupture from Fires on Gas Distribution in Tunnel

2. Introduction of nitrogen injection displacement for fire extinguishment

BioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.

Method of Determining the Threshold Pressure Gradient

Keywords: contra-rotating fan, air pressure, deformation of blades.

Enhanced Gas Drainage from Undersaturated Coalbed Methane Reservoirs

NUMERICAL SIMULATION OF STATIC INTERFERENCE EFFECTS FOR SINGLE BUILDINGS GROUP

PREDICTING the outcomes of sporting events

Study on Water Hammer Prevention in Pumping Water Supply Systems by Multi-valves

Vibration Analysis and Test of Backup Roll in Temper Mill

EXPERIMENTAL STUDY ON BEESWAX USING WATERJET DRILLING

Aerodynamic Performance of Trains with Different Longitudinal Section Lines under Crosswind

Broadly speaking, there are four different types of structures, each with its own particular function:

Optimized Gas Injection Rate for Underground Gas Storage; Sensitivity Analysis of Reservoir and Well Properties

A computational study on structural barrier to vertical spread of window spill plume along building exterior façade under external wind

u = Open Access Reliability Analysis and Optimization of the Ship Ballast Water System Tang Ming 1, Zhu Fa-xin 2,* and Li Yu-le 2 " ) x # m," > 0

SCIENTIFIC COMMITTEE SEVENTH REGULAR SESSION August 2011 Pohnpei, Federated States of Micronesia

Natural Gas Properties Analysis of Bangladesh: A Case Study of Titas Gas Field

Large olympic stadium expansion and development planning of surplus value

BioTechnology. BioTechnology. An Indian Journal FULL PAPER KEYWORDS ABSTRACT. Based on TOPSIS CBA basketball team strength statistical analysis

4 RESERVOIR ENGINEERING

Wind Characteristic Analysis of 100-Kilometer Wind Area in Lanzhou-Xinjiang High-Speed Railway

RESEARCH OF BLOCKAGE SEGMENT DETECTION IN WATER SUPPLY PIPELINE BASED ON FLUID TRANSIENT ANALYSIS ABSTRACT

Study on Fire Plume in Large Spaces Using Ground Heating

A COMPARATIVE STUDY OF PARAFFIN WAX

How to plan a safe and successful permeability test program in coal seams

EVALUATION OF GAS BUBBLE DURING FOAMY OIL DEPLETION EXPERIMENT USING CT SCANNING

Theory of multi-bench retaining for large area foundation pit and its engineering application

SIMULATION OF THE FLOW FIELD CHARACTERISTICS OF TRANSIENT FLOW

Research on Profile Control and Water Shut-off Performance of Pre-crosslinked Gel Particles and Matching Relationship between Particle and Pore Size 1

The Effect Analysis of Rudder between X-Form and Cross-Form

Chapter 8: Reservoir Mechanics

New power in production logging

Study on the Influencing Factors of Gas Mixing Length in Nitrogen Displacement of Gas Pipeline Kun Huang 1,a Yan Xian 2,b Kunrong Shen 3,c

Stress Analysis of The West -East gas pipeline with Defects Under Thermal Load

Accurate Measurement of Steam Flow Properties

W I L D W E L L C O N T R O L PRESSURE BASICS AND CONCEPTS

Procedia Engineering Procedia Engineering 2 (2010)

Effect of Gas-wetness on Gas-water Two-phase Seepage in Visual Microscopic Pore Models

Performance Analysis of Centrifugal Compressor under Multiple Working Conditions Based on Time-weighted Average

The Aerodynamic Design and Investigation of Loading Distribution of a Mixed Flow Compressor

A07 Surfactant Induced Solubilization and Transfer Resistance in Gas-Water and Gas-Oil Systems

SPE Cyclic Shut in Eliminates Liquid Loading in Gas Wells

Combined impacts of configurational and compositional properties of street network on vehicular flow

Developments in Longwall Ventilation

Numerical Simulation of the Basketball Flight Trajectory based on FLUENT Fluid Solid Coupling Mechanics Yanhong Pan

Using Experimental Procedure to Improve the Efficiency of the Two Stand Reversing Cold Mill

Wind Directional Effect on a Single Storey House Using Educational Wind Tunnel

Situated 250km from Muscat in

FIELD MEASUREMENTS OF SURFACE SUSPENDED SEDIMENT CONCENTRATION IN THE YANGTZE ESTUARY, CHINA

Analysis and Research of Mooring System. Jiahui Fan*

Experimental Research on Oxygen-enriched Air Supply System for Gasoline Engine

COAL SEAM METHANE - EXPLORATION AND PRODUCTION

EFFECTS OF IMPORT AND INVENTORY AMOUNTS ON CHANGES IN WHOLESALE PRICES OF SALMON IN JAPAN

Analysis of Steady-state Characteristics of Aerodynamic Bearing Based on FLUENT JIA Chenhui 1,a, Du Caifeng 2,b and QIU Ming 3,c

DRILLSCENE PROACTIVE DRILLING DECISIONS

Quantifying Bullwhip Effect and reducing its Impact. Roshan Shaikh and Mudasser Ali Khan * ABSTRACT

Structural Design and Analysis of the New Mobile Refuge Chamber

Availability, Economics and Potential of North American Unconventional Gas

CALIBRATION OF THE PLATOON DISPERSION MODEL BY CONSIDERING THE IMPACT OF THE PERCENTAGE OF BUSES AT SIGNALIZED INTERSECTIONS

Aalborg Universitet. Published in: Proceedings of Offshore Wind 2007 Conference & Exhibition. Publication date: 2007

CAN MUDLOGGING SHOWS BE QUANTIFIED?

Aerodynamic Performance Comparison of Head Shapes for High-Speed Train at 500KPH

CBM Geology & Well Design 5TH ANNUAL CBM & UNCONVENTIONAL GAS WEDNESDAY 27 JUNE 2012 GEOFF BARKER, PARTNER

Keywords: multiple linear regression; pedestrian crossing delay; right-turn car flow; the number of pedestrians;

Dynamics of the bubble near a triangular prism array

International Journal of Petroleum and Geoscience Engineering Volume 03, Issue 02, Pages , 2015

NUMERICAL SIMULATION OF WIND INTERFERENCE EFFECT

CSU_Yunlu 2D Soccer Simulation Team Description Paper 2015

by Michael Young Human Performance Consulting

Transcription:

International Conference on Energy and Environmental Protection (ICEEP 2016) High production indexes and key factors influencing coalbed methane (CBM) horizontal well productivity Yanhui Yang1, a, Shu Tao2, b,yanlei Yang1, c, Wei Tian1, Yongping Zhang1, Ning Wei1 1 2 PetroChina Huabei Oilfield Company, Renqiu Hebei, 062552, China School of energy resources, China University of Geosciences (Beijing), Beijing, 100083, China a yjy_yyh@petrochina.com.cn, b cugbst@126.com, c yjzx_yyl@petrochina.com.cn, Keywords: Fanzhuang block; horizontal well; geological factors; high production index Abstract: To evaluate the development effects systematically, the coalbed methane (CBM) production performance of 58 horizontal wells in the Fanzhuang block of the Qinshui Basin was tracked and correlations between the gas rate per 100 m and the influencing factors were discussed by single and multiple factor analysis. The grey correlation analysis (GRA) was then applied to quantitatively evaluate the key geological factors influencing CBM horizontal well productivity. The results show that the order of geological factors influencing horizontal well productivity is burial depth, gas production capacity (C), gas production potential (P), ratio of critical desorption pressure to initial reservoir pressure (pc/pi), gas content, and thickness, where the gas production potential and capacity can reflect the comprehensive impacts of thickness, gas content, and permeability. The larger the gas production potential or capacity is, the better the gas productivity is. The high production indexes are the depth of 300~500m, pc/pi >0.6, P>130m4/t, C>60 md m4/t, and far away from the faults. Introduction Coalbed methane (CBM) is an important source of clean energy which has been a strategic supplement to conventional natural gas [1]. As the largest consumer and producer of coal in the world, China has the third largest CBM reserves, behind those of the Russia and Canada [2-3]. With the rapid development of CBM industrialization, lower costs and higher gas production per well in CBM commercial development have gradually become important topics in CBM research.contrast with the vertical wells, the horizontal well is one of the important drilling technologies for the development of coalbed methane, and is also the technical development trend of the development of coalbed methane, which integrates drilling, well completion and stimulation on production [4]. In view of the low permeability, strong heterogeneity of coal reservoirs in china, Huabei Oilfield introduce horizontal wells draw foreign experience of low rank coal mining during "The 11th Five" periods, and take the lead in the implementation of horizontal well development in in Fanzhuang block, southern Qinshui Basin.However, the difference of gas and water production is gradually emerging with the development of CBM horizontal wells. Compared with the vertical well, the horizontal well has a good effect of increasing production while the higest production well even up to 15 times around the vertical wells.but there are still some questions: High variation in the production of horizontal wells,and there are a lot of inefficient wells; The key factors that affect the development of horizontal wells are not clear; Low permeability and strong heterogeneity of high rank coal reservoir is difficult to optimize the high yield region of the reservoir. In order to evaluate the development effects systematically, the coalbed methane (CBM) production performance of 58 horizontal wells in the Fanzhuang block of the Qinshui Basin was tracked to point out the problems existing in the development of horizontal wells. Besides, correlations between the gas rate per 100 m and the influencing factors were discussed by single and multiple factor analysis. The grey correlation analysis (GRA) was then applied to quantitatively evaluate the key geological factors influencing CBM horizontal well productivity with the goal of providing theory and engineering support for the development of CBM. 2016. The authors - Published by Atlantis Press 78

The main problems in the development of multi branch horizontal wells (1) Coal reservoir prone to sensitivity damage: the coal seam is an organic reservoir with low permeability and easy to be polluted [5-6]. Compared with vertical wells, horizontal wells is very easy to appear reservoir sensitivity damage in the drilling process, which resulting in effective flow channel blocked and the permeability reduction.these questions lead to the production of water in unit length is below the adjacent vertical wells, and desorption pressure is difficult to reach quickly. Statistical the water production of 58 horizontal wells and the adjacent vertical wells in Fanzhuang block, we found that the water production per unit length only 0.25 m 3, which is far lower than surrounding straight wells (average of 3.34 m 3 ). The desorption pressure of horizontal well is also much lower than that of the vertical well, as shown in Table 1. Table 1 Water production of horizontal well and adjacent vertical well Desorption Water Water production Well type N pressure production per unit length (Well) (MPa) (m 3 ) (m 3 ) Horizontal 58 1.61 1089 0.29 Vertical 130 1.82 564 3.53 (2) Borehole wall collapse caused by the instability of wellbore in the process of drilling: due to the attitude, coal structure, roof and high rank inherent mechanical properties of coal reservoirs, the wellbore are are easy to collapse during drilling process, which result in the diffierence between the actual horizontal trajectory and the design trajectory. These wells is is difficult to realize area reduction with low resource utilization. (3) A large number of inefficient wells with no improvement measures: the horizontal well always use barefoot well completion, which is unable to carry out the production increase. The borehole wall will collapse in a few years after the coal seam gas exploitation, which greatly shortens the gas recovery life and reduces the efficiency of the exploration and development of CBM. It can be seen that the production of 38 horizontal wells is lower than 3000m 3 /d, accounted for 66% of the total number of wells (Fig.1). Fig. 1 Pie chart of horizontal well productivity Effects of the geological factors on the gas production There are some inefficient horizontal wells (<3000m 3 /d) in the Fanzhuang block, souther Qinshui basin. The geological condition is a decisive component for the high production of horizontal wells. In this study, five reservoir parameters (including the burial depth, the coal thickness, the ratio of the critical desorption pressure to the initial reservoir pressure, the gas content and the permeability) of 35 wells (away from faults with no engineering collapse) were analyzed with correlation analysis to reveal the factor that has the greatest effect on gas production of horizontal wells in the Fanzhuang block. Moreover, in order to avoid the influence of coal seam footage on gas production, the average gas rate per 100m was calculated to represent the development effect of horizontal wells. 79

Single factors Burial depth In the Fanzhuang Block, the middle depth of the target coal seams was considered the burial depth, which was used to analyze the correlation of the burial depth and gas production. Fig.2a shows that the burial depth of the target coal seam ranges from 273 to 812 m based on 35 data points, and there is a negative relationship between the burial depth and average gas rate per 100m. Overall, the average gas rate per 100m decreases with the increase of the buried depth and high gas production (>150m 3 /(d/100m)) are mainly concentrated in the depth of 250~600m.Besides, the burial depths more than 630m is of poor gas recovery (<50m 3 /(d/100m)). Therefore, a coal seam with a depth of from 250 to 600 m tends to have a much greater potential for producing CBM with horizontal wells. Fig. 2 Relation between average gas rate per 100m and burial depth(a), thickness(b). Thickness Coal seam is the fundamental for CBM accumulation [7]. Under the same conditions, the greater the thickness of the coal, the more abundant gas source with the stronger the ability of gas supply, which is favorable to the higher gas production. The thickness of coal seam can be determined according to the well logging interpretation. Fig.2b indicate that the thickness of the target coal seam ranges from 4.5 to 7 m. However, the gas production of the horizontal well has a trend of decrease with the increase of the thickness of the coal seam, which is obviously contrary to the above understanding. This phenomenon indicates that the productivity of horizontal well is affected by many factors but not a single factor. Besides, the vertical heterogeneity of permeability may be an important reason for the decrease of gas rate. Ratio of the critical desorption pressure to the initial reservoir pressure In general, the greater the depth, the greater the reservoir pressures with the more powerful the reservoir compaction, which result in lower permeability and lower productivity of CBM wells [8]. The average pressure coefficient of coal reservoirs (0.67MPa/100m) in the study area can be used to calculate the initial reservoir pressure(p c ). The results show that the p c in the study area is 2 to 6MPa. Fig.3a indicates that the average gas rate per 100m has a negative relationship with the p c. When the p c is more than 4MPa, the wells show a poor gas recovery. For development wells, when successive gas appears, the flowing bottom hole pressure obtained from real-time monitoring system can be regarded as the critical desorption pressure(p i ). It can be seen that the greater the p i, the higher production of CBM wells (Fig.3b). When the p i is in the range of 2 to 2.5MPa, the gas rate of horizontal wells is basically more than 150m 3 /(d/100m)). Actually, the horizontal well is more complex than the vertical well due to the horizontal well has some branches and veins branch, the ratio of the p c to p i is more significance for the development of CBM horizontal wells [9]. It can be seen that p c /p i has a positive relationship with the gas rate per 100 m. When the p c /p i is lower than 0.6, the horizontal well performed a poor gas recovery; Contrary, the gas rate per 100 m is basically greater than 150m 3 /(d/100m) while the p c /p i is more than 0.6(Fig.3c). 80

Gas content Fig. 3 Scatter map of the average gas rate per 100m and initial p c (a), p i (b) and p c / p i(c). Gas content is the material base for the development of CBM wells, while the CBM wells will get high and stable production in the late production process due to adequate gas supply [10]. In this work, the gas content was calculated according to Langmuir parameters and the petrophysical properties of the reservoir. The gas content in the study area is generally high (10 to 25 m 3 /t). Fig.4a indicate that the gas content has a significant influence on the gas production of CBM horizontal well, there are almost no high-yield wells in lower gas content area (<15 m 3 /t). When the gas content is more than 20m 3 /t, the proportion of high yield wells increased significantly. Fig. 4 Scatter map of the average gas rate per 100m and the gas content(a),permeability(b). Permeability Although the branch and vein branch of horizontal well can solve the difficult problem of low permeability reservoir, the drainage pressure drop velocity and the expansion of the pressure drop funnel are still affected by the development of fracture. According to the relationship between well testing permeability and depth, the permeability of coal reservoirs in actual CBM wells can be predicted. The results show that the coal reservoirs in Fanzhaung block is a typical extremely-low permeability reservoir with the permeability lower than 1mD. There is a slightly positive relation between the gas rate per 100m and permeability, but it is worth noted that the horizontal well is of poor production when the permeability lower than 0.3 md(fig.4b). Multi-factor A coal seam with a larger thickness contributes to more controlled reserve of the CBM wells, which results in stronger gas supply capability and higher gas production. In the Fanzhuang Block, less thickness correlates with greater gas production (Fig. 2b), showing that CBM production is the result of multiple factors, and a single factor (e.g., thickness) does not directly affect CBM production. Thus, two parameters, the gas production potential (P=the thickness the gas content) and the gas production capacity (C= P the permeability), characterize the comprehensive effect of the thickness, the gas content and the permeability. Fig. 5 indicate that P and C have a good positive correlation with the gas rate per 100m of the horizontal wells. When C is more than 60mD m 4 /t or P is more than 120 81

m 4 /t, the production of horizontal well is basically greater than 150m 3 / (d/100m), which indicate that high gas production potential and capacity are favorable for gas productivity. Fig.5 Scatter map of the average gas rate per 100m and the gas production potential(a), the gas production capacity(b). Gray correlation analysis and reservoir classification Gray correlation analysis is an important part of the gray system theory, which is also be frequently used to determine whether a close relationship exists between array curves, based on their similarity. The purpose of gray correlation analysis is to quantify the correlation degrees of various factors and present the key factors that affect the system development [11]. In this work, the relationship between the factors and the gas rate per 100m was determined by the method of gray correlation analysis. The results show that the order of geological factors influencing horizontal well productivity is burial depth, C, P, p c /p i,, gas content, and thickness (Table 2). In the process of CBM development, the burial depth determine the permeability pressure and the development of fracture, while the ratio of the critical desorption pressure to the initial reservoir pressure frequently determines the degree of difficulty of the water drainage and decompression in the CBM reservoir. Besides, it can be seen that P and C can reflect the comprehensive impacts of thickness, gas content, and permeability. The key geological factors affecting CBM horizontal wells productivity in the Fanzhuang Block can be characterized as burial depth, p c /p i, P and C. According to the relationship between these four factors and gas production per 100m of actual CBM wells, the coal reservoir is divided into three categories(table 3): Table 2 The results from the grey correction analysis. Factors Correlation degree Ran k Burial depth 0.7941 1 C 0.6865 2 P 0.6271 3 p c/p i 0.5244 4 Gas content 0.5041 5 Thickness 0.5017 6 Table 3 classification of the reservoirs of horizontal wells Type Average gas rate per 100m (m 3 /(d/100m)) Depth(m) p c/p i P(m 4 /t) C(mD m 4 /t) I >150 <500 >0.6 >120 >60 II 50~150 500~700 0.2~0. 100~120 40~60 6 III <50 >700 <0.2 <100 <40 (I): Burial depth < 500m, p c /p i >0.6, P>120m 4 /t,c>60md m 4 /t. This kind of coal reservoir has good development effect and high economic benefit with the high gas recovery (>150 m 3 / (d/100m)). 82

(II): Burial depth is 500 to 700m, 0.2<p c /p i <0.6, 100<P <120mD m 4 /t, 40<C <60mD m 4 /t.the development of this kind of coal reservoir basically has no economic benefit, while the gas production per 100m is about 50~150 m 3 / (d/100m). (III): Burial depth>700m, p c /p i <0.2, P<100mD m 4 /t, C <40mD m 4 /t.the development effect of this kind of reservoir is poor with no economic benefit, which should be avoided in the process of development as soon as possible. Conclusions a. The coalbed methane (CBM) production performance of 58 horizontal wells in the Fanzhuang block of the Qinshui Basin was tracked, and the main problems in the development of horizontal wells was addressed: Coal reservoir is sensitive to reservoir sensitivity damage; Borehole wall collapse caused by the instability of wellbore in the process of drilling; There are a large number of inefficient wells with no improvement measures. b. The order of geological factors influencing horizontal well productivity is burial depth, gas production capacity (C), gas production potential (P), ratio of critical desorption pressure to initial reservoir pressure (p c /p i ), gas content, and thickness, where the gas production potential and capacity can reflect the comprehensive impacts of thickness, gas content, and permeability. c. Three types of reservoir model related to the development of CBM horizontal wells was quantitative characterized according to the combination of four parameters (including burial depth, ratio of critical desorption pressure to initial reservoir pressure, the gas production potential and capacity). The high production indexes are the depth of 300~500m, p c /p i >0.6, P>130m 4 /t, C>60 md m 4 /t, and far away from the faults. Acknowledgments This work was financially supported by the Key Project of the National Science and Technology (Grant No. 2011ZX05061), the Project of the University-enterprise Cooperation of PetroChina Huabei Oilfield Company Development optimization design and matching technology optimization in the coalbed methane (CBM) favorable areas, southern Qinshui Basin (Grant No. HBYT-YJY-2015-JS-300) and the Project of the PetroChina Huabei Oilfield Company Research on High-Performance Development Technology in Eastern Mabi Block (Grant No. 2016-HB-M06). We appreciate the constructive comments and suggestions from the reviewers and the editors. References [1] Y. Lv, D. Tang, H. Xu, and H. Luo: Int. J. Coal Geol. Vol. 96-97(2012), p. 93-108 [2] Z. Chen, Y. Wang, J. Yang, X. Wang, Y. Chen, and Q. Zhao: Acta Petrol. Sin. Vol. 30(2009), p. 409-416 [3] H. Xu, D. Tang, S. Tang, J. Zhao, Y. Meng, and S. Tao: Int. J. Coal Geol. Vol. 121(2014), p. 44-52 [4] Sinayuc. C, Shi, JQ, Imrie. C, Syed. S., Korre. A, and Durucan. S: Energy Procedia, Vol. 22(2011), p.2150-2156. [5] Y. Li, D. Tang, H. Xu, and T. Yu: J. Petrol. Sci. Eng. Vol. 122(2014), p. 488-496 [6] Y. Yang, D. Tang, Y. Yang, L. Chen, and S. Tao: Coal Sci. Technol. Vol. 43(2015), p. 96-99 [7] Pashin. J. C: Int. J. Coal Geol. Vol. 82(2010), p. 135-146. [8] Shu. T, Y. Wang, D. Tang, H. Xu, Y. Lv, W. He, Y. Li: Int. J. Coal Geol. Vol. 93(2012), p. 16-22 [9] McKee. C. R., Bumb. A. C: SPE Form Evaluation. Vol. 2(1987), p. 599-608. 83

[10] Ross. H. E, Hagin. P, Zoback. M. D: International Journal of Greenhouse Gas Control. Vol.3 (2009), p.774-786. [11] Deng. J. L: Systems & Control Letters. Vol.1 (1982).p.288-294 84