PhD Thesis Research FE 2 approach for the modelling of coupled - - PowerPoint PPT Presentation

phd thesis research
SMART_READER_LITE
LIVE PREVIEW

PhD Thesis Research FE 2 approach for the modelling of coupled - - PowerPoint PPT Presentation

PhD Thesis Research FE 2 approach for the modelling of coupled hydro-mechanical behaviour of partially saturated coalbeds Franois B ERTRAND Supervised by: Frdric C OLLIN University of Lige November 2015 Introduction Coalbed methane


slide-1
SLIDE 1

PhD Thesis Research

FE2 approach for the modelling of coupled hydro-mechanical behaviour

  • f partially saturated coalbeds

François BERTRAND

Supervised by: Frédéric COLLIN

University of Liège

November 2015

slide-2
SLIDE 2

Introduction

Coalbed methane

NORTH SEA GERMANY THE NETHERLANDS FRANCE LUX. Brussels Antwerpen Hasselt Mons Liege Namur

50 km

3°E 4°E 5°E 6°E 50°N 51°N

N

BELGIUM

Figure: Map of the outcropping or shallow

subsurface coal basins (shaded area) in and around

  • Belgium. Modified after [Piessens and Dusar, 2006].

6 to 9 billions tons of coal in the walloon basement 100 to 300 billions m3 of methane ? ? 10 to 20 billions m3 of recoverable methane? 3 to 10 years of gas consumption in Belgium

François BERTRAND (University of Liège) PhD Thesis Research November 2015 1 / 12

slide-3
SLIDE 3

Introduction

Coalbed methane

Coalbed methane (CBM) = unconventional resource

Source rock (= coalbeds) is also the reservoir for the methane

Coalbed methane

Conventional reservoir

Seal Source Reservoir

Coal Promising alternative to conventional gas

François BERTRAND (University of Liège) PhD Thesis Research November 2015 2 / 12

slide-4
SLIDE 4

Introduction

Coal properties Coalbeds = dual porosity systems

Micropores + Macropores

⇐ ⇒

Matrix + Cleats

Figure: From Schlumberger Oilfield Glossary Figure: From Schlumberger Oilfield Glossary

François BERTRAND (University of Liège) PhD Thesis Research November 2015 3 / 12

slide-5
SLIDE 5

Introduction

Coal properties Coalbeds = dual porosity systems

Micropores + Macropores

⇐ ⇒

Matrix + Cleats

Figure: From Schlumberger Oilfield Glossary

Cleat permeability is directly dependent

  • n the width of the cleats.

François BERTRAND (University of Liège) PhD Thesis Research November 2015 3 / 12

slide-6
SLIDE 6

Introduction

Challenge "The influence of geomechanical processes is an important issue for coalbed methane recovery and ignoring geomechanical processes may lead to errors in the evaluation of coalbed methane production." [Gu et al., 2005] But Coupled hydro-mechanical behaviour of coalbeds = complex + still defeats a comprehensive description

François BERTRAND (University of Liège) PhD Thesis Research November 2015 4 / 12

slide-7
SLIDE 7

Introduction

Challenge "The influence of geomechanical processes is an important issue for coalbed methane recovery and ignoring geomechanical processes may lead to errors in the evaluation of coalbed methane production." [Gu et al., 2005] But Coupled hydro-mechanical behaviour of coalbeds = complex + still defeats a comprehensive description

Macroscale FE computation

Current models

Phenomenological constitutive models

Laws at macroscale for microscopic phenomena

François BERTRAND (University of Liège) PhD Thesis Research November 2015 4 / 12

slide-8
SLIDE 8

Introduction

Challenge "The influence of geomechanical processes is an important issue for coalbed methane recovery and ignoring geomechanical processes may lead to errors in the evaluation of coalbed methane production." [Gu et al., 2005] But Coupled hydro-mechanical behaviour of coalbeds = complex + still defeats a comprehensive description

Macroscale FE computation

Current models FE² method

Modelling at lower scale Phenomenological constitutive models Microscale FE computation Homogenization

Laws at macroscale for microscopic phenomena

Apply a multi-scale method taking advantage of the periodical structure of coal.

François BERTRAND (University of Liège) PhD Thesis Research November 2015 4 / 12

slide-9
SLIDE 9

Contents

1

Introduction

2

Overview of the thesis Work environment Workpackage 1: Development of a multiphase flow model Workpackage 2: Development of a mechanical model Workpackage 3: Coupled hydro-mechanical model Workpackage 4: Reservoir modelling Work plan

3

Conclusion

François BERTRAND (University of Liège) PhD Thesis Research November 2015 5 / 12

slide-10
SLIDE 10

Work environment

Université de Liège

Université Joseph Fourier University of Newcastle

(France, Grenoble) (Australia)

Geomechanics and Engineering research unit Thesis supervised by F. Collin Team Pr O. Buzzi

Expertise on coal

Laboratoire 3S-R

Thesis on FE² by B. Van den Eijnden Expertise in computational geomechanics (FE code Lagamine)

François BERTRAND (University of Liège) PhD Thesis Research November 2015 6 / 12

slide-11
SLIDE 11

Workpackage 1: Multiphase flow model

FE2 method :

REV

Microscale FE computation

Constitutive equations (flow law, storage law) are applied only on the microscopic scale. Task 1.1: Extension of the FE2 method to unsaturated conditions. Task 1.2: Extension to multiphase flow conditions (liquid and gas). Task 1.3: The developed model is faced and compared with others models found in the literature:

[Pan and Connell, 2012], [Shi et al., 2014].

François BERTRAND (University of Liège) PhD Thesis Research November 2015 7 / 12

slide-12
SLIDE 12

Workpackage 2: Development of a mechanical model

The mechanical behaviour of the coalbed results from the geometry and the properties of the components. Triaxial experimental results ⇔ REV responses Task 2.1: Evaluate numerically the influence of the cleat density in the REV on the macroscopic behaviour of the material. Task 2.2: Identification of the material parameters through a back-analysis of some experimental results.

François BERTRAND (University of Liège) PhD Thesis Research November 2015 8 / 12

slide-13
SLIDE 13

Workpackage 3: Coupled hydro-mechanical model

Task 3.1: Experimental campaign Hydro-mechanical tests on an Australian coal:

Shrinkage/swelling properties and permeability measurements

Task 3.2: Coupled model implemented in the FE code LAGAMINE.

Macroscale Microscale

REV

FE computation

François BERTRAND (University of Liège) PhD Thesis Research November 2015 9 / 12

slide-14
SLIDE 14

Workpackage 4: Reservoir modelling

A synthetic reservoir model with one production well is simulated from the drilling of the well to the production period. The necessity or not of stimulation techniques is evaluated.

François BERTRAND (University of Liège) PhD Thesis Research November 2015 10 / 12

slide-15
SLIDE 15

Work plan

1st year 2nd year 3rd year 4th year

( )

: Multiphase flow model : Mechanical model : HM couplings : Reservoir modelling : Literature review : Paper writing : Thesis writing : Newcastle stay : Grenoble stay

( )

François BERTRAND (University of Liège) PhD Thesis Research November 2015 11 / 12

slide-16
SLIDE 16

Conclusion

High potential for coalbed methane in Wallonia Challenge = better understand methane migration to optimize recovery. Development of a coupled hydro-mechanical model for coalbed methane reservoirs in the framework of the finite element square method (FE2).

François BERTRAND (University of Liège) PhD Thesis Research November 2015 12 / 12

slide-17
SLIDE 17

Thank you for your attention!

FE2 approach for the modelling of coupled hydro-mechanical behaviour

  • f partially saturated coalbeds

Related work: Efficiency of shaft sealing for CO2 sequestration in coal mines,

Presentation at the Workshop on Geomechanics & Energy EU 2015 (held from 13-10-2015 to 15-10-2015 in Celle, Germany).

slide-18
SLIDE 18

References

Gu, F., Chalaturnyk, J., et al. (2005). Analysis of coalbed methane production by reservoir and geomechanical coupling simulation. Journal of Canadian Petroleum Technology, 44(10). Pan, Z. and Connell, L. D. (2012). Modelling permeability for coal reservoirs: a review of analytical models and testing data. International Journal of Coal Geology, 92:1–44. Piessens, K. and Dusar, M. (2006). Feasibility of CO2 sequestration in abandoned coal mines in belgium. Geologica Belgica. Shi, J.-Q., Pan, Z., and Durucan, S. (2014). Analytical models for coal permeability changes during coalbed methane recovery: Model comparison and performance evaluation. International Journal of Coal Geology, 136:17–24.

slide-19
SLIDE 19

SWOT Analysis

S W O T

Strengths Weaknesses Opportunities Threats

Positive Negative Internal External

  • Personal profile
  • Work environment
  • Strong modelling basis
  • Clear objectives
  • Experiments on coal required

to identify model parameters

  • Innovative approach to model

coal behaviour

  • Industrial interest
  • High CBM potential in Wallonia
  • Unfavorable political decision

about CBM ?

  • No CBM production expertise

in Wallonia

slide-20
SLIDE 20

Coalification

C a r b

  • n

c

  • n

t e n t W a t e r c

  • n

t e n t

Gas content of coals

Modified after Kentucky Geological Survey (University of Kentucky)

slide-21
SLIDE 21

Global overview

Shale gas Tight gas Coalbed methane

Country CBM resources [1012m3] CBM annual production [109m3] Russia 17.0 - 113.0 1 Canada 17.9 - 76.0 9 China 36.8 5 USA 21.5 52 Australia 8.0 - 14.0 4

From IEA (2010)

slide-22
SLIDE 22

Production

From Schlumberger (2009), Coalbed Methane: Clean Energy for the World

Methane Water Volume Time

Dewatering Stable production Decline

slide-23
SLIDE 23

CBM versus Shale gas

Coalbed methane

Conventional reservoir Shale gas Mines

Seal Source R e s e r v

  • i

r Coal

Tight sand gas

slide-24
SLIDE 24

FE2

1. discretised by finite elements Macroscopic structure 3. assigned at each macroscopic IP REV 4. : apply appropriate from the macroscopic deformation gradient tensor Localization displacements to the REV 5. : stress and deformation distributions in the REV Microscale FE computation 6. : REV returned to the macroscopic IP Homogenization averaged stress

  • 7. Macroscopic internal nodal forces

2. c from the estimation

  • f the macroscopic nodal displacements relative to the external load

Macroscopic deformation gradient tensor omputed for each IP

  • 8. Macroscopic stiffness matrix

9. between external load and internal load? Balance Next time increment evaluated

+

  • Updated estimation

displacements

  • f the nodal

required (via macroscopic stiffness matrix)