Unconventional Reservoirs Require Unconventional Analysis ...

Report 10 Downloads 135 Views
Primary funding is provided by

The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow their professionals to serve as lecturers Additional support provided by AIME

Society of Petroleum Engineers Distinguished Lecturer Program www.spe.org/dl

1

Unconventional Reservoirs Require Unconventional Analysis Techniques David Anderson

Society of Petroleum Engineers Distinguished Lecturer Program www.spe.org/dl

2

This Presentation…  Introduction to rate transient analysis (RTA)  The challenge of analyzing unconventionals

 Current methodologies – how has the technology evolved?  The future of production analysis and modeling  Probabilistic approach  Field examples 3

Rate Transient Analysis (RTA) is the science (and art) of extracting useful information about the reservoir, completion and/or surface operations based on the interpretation, analysis and modeling of continuous measurements of production volumes and flowing pressures from a single well.

4

Concept of Rate Transient Analysis Company: On Stream: 03/28/2013 Field: Current Status: Flowing

Gp: 1775 MMscf Np: 0.000 Mstb Wp: 0.000 Mstb Qcond: 0.000 Mstb

W ell 01

3600

5200

- Production occurs under changing constraints - Reservoir “signal” may be in rates or pressures (or both)

3200

2800

2400

4000

3600

3200 1600

2800 1200

2400 800

2000 400 1600 0 1200

-400

800

-800

400

-1200

0 0

2

4

6

8

10

12

14

16

18

20

22

24

26

28

Time (days) Normalized Time (month)

30

32

34

36

38

40

42

44

46

48

5

Run Depth Pressure (psi(a))

Op Gas Rate (Mscfd)

4400

pwf (psia)

q (Mscfd)

2000

4800

Concept of Rate Transient Analysis Comparison View 9 . 10-3

Legend

6

Normalized Gas Rate vs. Normalized Time Normalized Gas Rate vs. MBT (2)

4 3

Normalized Gas Rate (MMscfd/psi)

q/Dp (Mscfd/psi)

2

10-3

Instantaneous normalization Superposition (Material Balance Time)

6 4 3 2

10-4

6 4 3 2

10-5 0

10

20

30

40

50

60

70

80

90

100

110

120

130

140

150

160

130

140

150

160

Normalized Time (month)

0

10

20

30

Time, Material Balance Time (months) 40

50

60

70

80

90

100

110

120

MBT

6

Type Curve Analysis – Characterize Reservoir Comparison View 4 . 100

Legend q/D  - TC Normalized Gas Rate vs. MBT

2

Log-Log Plot - Identify flow regimes Boundary Dominated Flow (connected HCPV)

1.0

4

10-1

4 2 10-2

4

Transient Flow (permeability, skin)

D

q/Dp (Mscfd/psi)

2

2 10-3

4 2 10-4

Adapted from Palacio and Blasingame: “Decline-Curve Analysis Using Type Curves” (SPE 25909) 1993

5 3 2 . 10-5 10-2

2

3 4 5 6 7 10-1

2

3 456

1.0

2

3 4 56

101

2

3 4 56

102

2

3 4 56

103

2

3 4 5 6 7 104

Material Balance Time (days)

2

3 456

105

2

3 4

7

8 . 105

Flowing Material Balance – Estimate HCPV Company: On Stream: 10/01/2002 Field: Apollo Current Status: Unknown

7000

30

6500

28

Legend Gas Rate Flowing Pressure

26

6000

24

Production Rate (Mscfd)

5500

Gp: 3409 MMscf Np: 224.268 Mstb Wp: 16.566 Mstb Qcond: 0.000 Mstb

Example 1

22

5000

Pressure (psi(a))

pwf (psia)

4000

3500

3000

2500

Operated Gas Rate (MMscfd)

20

4500

18

Measured flowing pressure

16

14

12

10

2000

8 1500

1000

4

500

2

0

Measured rate

6

0 0.00

0.50

1.00

1.50

2.00

2.50

3.00

3.50

4.00

4.50

5.00

5.50

6.00

6.50

7.00

7.50

Cumulative Production (bcf) Cumulative Gas Production (Bscf)

8.00

8.50

9.00

9.50

10.00

8

10.50

Flowing Material Balance – Estimate HCPV Company: On Stream: 10/01/2002 Field: Apollo Current Status: Unknown

7000

30

6500

28

Example 1

Legend Flowing p/Z** Gas Rate Flowing Pressure

26

6000

Calculated p/z

24

Production Rate (Mscfd)

5500

Gp: 3409 MMscf Np: 224.268 Mstb Wp: 16.566 Mstb Qcond: 0.000 Mstb

22

Pressure, p/Z** (psi(a))

4500

4000

3500

3000

2500

2000

p  p     qbpss z  z  wf

20

Operated Gas Rate (MMscfd)

pwf and p/z (psia)

5000

18

16

14

12

10

8 1500

6

1000

4

500

2

0

0 0.00

- Mattar L., Anderson, D., Dynamic Material Balance – Oil or Gas

Original Gas-In-Place

In Place Without Shut-ins - 2, CIPC 2005-113 0.50

1.00

1.50

2.00

2.50

3.00

3.50

4.00

4.50

5.00

5.50

6.00

6.50

7.00

7.50

Cumulative Production (bcf) Cumulative Gas Production (Bscf)

8.00

8.50

9.00

9.50

10.00

9

10.50

Modeling – Validate and Forecast Results 11500

9.00

11000

10000

10500

9500

8.50 8.00

10000

7.50

9500

9000

Production Forecast

8500

9000

7.00

8000

8500 7500

6.50

7000

7500

5.50

7000

6500

6500

6000

6000

5500

5500

5000

5.00 4.50 4.00 3.50

Gas Cum (MMscf)

6.00

5000

4500

4500

4000

4000 3.00

3500 3500

2.50

Pressure match

3000

3000 2500

2.00

2500

1.50

2000

2000

1500

1500

1000

1000

500

500

1.00 0.50 0.00

0 2002

2003

2004

2005

2006

Pressure (psi(a))

Cal Gas Rate (MMscfd)

8000

Benefits of RTA  Evaluation of reserves  Reliable early evaluation- choked wells  Scientific support for reserves auditors  Dynamic reservoir characterization  Estimate permeability and in-place hydrocarbons  Estimate completion effectiveness  Calibrate reservoir simulation models  Reservoir surveillance  Distinguish productivity fall-off from depletion  Identify optimization candidates

The Challenge of Analyzing Unconventionals…  Unconventional reservoirs are more complex  Complex, non-uniform fracture networks  Reservoir properties are significantly altered by completion  Low permeability – long term transient flow  Drainage area continually expands  Difficult to distinguish clear drainage boundaries

Flow Characterization – Conventional vs. Unconventional Radial Flow - Conventional

Linear Flow - Unconventional

 Fluid flows to the sandface  Pressure drawdown localized at sandface

 Fluid flows to the fracture(s)  Pressure drawdown throughout fracture(s)

High Permeability Low Contact Area

Low Permeability High Contact Area

Unconventional Analysis Methods… Square Root Time Plot- Linear Flow Simple A = 4 nf xf h h 2 xf

A k Complex A 

Dp q Skin

Dp  m t b q t

A

f

Boundaries and Drainage – Conventional vs. Unconventional a) Conventional Reservoir

b) Unconventional Reservoir Horizontal Wells

Vertical Wells

Parallel Fractures

Fracture interference Stimulated Reservoir Volume (SRV)

Geological features

Parallel and Orthogonal Fractures

Well interference 15

Unconventional Analysis Methods… Flowing Material Balance Stimulated Reservoir Volume

h

2 xf Le

p wf z In-place hydrocarbons (SRV) Cumulative Productioni

Unconventional Analysis Methods… Simplified Approachtetf = SRV

A√k

Anderson et al 2010, Analysis of Production Data from Fractured Shale Gas Wells – SPE 131787

skin

SRV

Assume – uniform fractures Calculated- xf, k, skin, SRV

Illustrating the Challenge of Analyzing Unconventionals Simulation of flow into a complex fracture network in a gas shale Company: On Stream: 25/06/2013 Field: Current Status: Flowing

Gp: 705 MMscf Np: 0.000 Mstb Wp: 0.000 Mstb Qcond: 0.000 Mstb

W ell 02

3 . 104

540 520 500 480

2

460 440 420

360

9

340 8

320 7

300 6

280 260

5

240 4

220 200 180

3

2

103

Non-uniform frac length, spacing and conductivity Ultra low matrix permeability Six months production, constant pwf

160 140 120 100 80 60 40 20 0

June

July

August

September

October 2013

November

December

Run Depth Pressure (psi(a))

Op Gas Rate (Mscfd)

380 104

pwf (psia)

q (MMscfd)

400

tetf = SRV

A√k skin

Contacted HCPV

q/Dp (MMscfd/psi)

Dp/q (psi/MMscfd)

Simplified Approach – Bulk Reservoir Properties

Gpa/ceDp (bcf)

Square Root Time

q/Dp (MMscfd/psi)

SRV = 0.75 bcf

Time (d)

Contacted HCPV = 1.1 bcf

Simplified Approach – Comparison of Analyzed Reservoir Properties with Actual

As Analyzed Stimulated reservoir width = 120 ft k (stimulated zone) = 0.011 md k (matrix) = 0.0005 md Contacted OGIP = 2 bcf

Actual

Hz fractures – 250 ft, FCD=50 Vert fractures – 500 ft, FCD = 100 k (matrix) = 0.0001 md OGIP = 46 bcf (1 section)

Simplified Approach – Comparison of Analyzed SRV with Actual

As Analyzed

Actual

SRV = 0.75 bcf

SRV ~ 0.75 bcf

Simplified Approach – Comparison of 5 year Production Forecasts Comparison View

2 . 101

Gp = 1.8 bcf

101

Gp = 1.9 bcf

8 6

q (MMscfd)

Cal Gas Rate (MMscfd) / Rate Forecast 1 (MMscfd)

5 4 3

2

1.0 8 6 5 4 3

2

10-1 2013

2014

2015

2016

Time (years)

2017

2018

2019

Field Example - Bakken Oil Company: On Stream: 06/05/2008 Field: Undefined Field Current Status: Flowing

2 . 103

Gp: 72 MMscf Np: 98.451 Mstb Wp: 19.879 Mstb Qcond: 0.000 Mstb

Bakken Oil Bakken

103

5200 5000

7

4800

103

5 4600 4

7

4400 3

5

4200

4 2

4000

3 3800 2

3600 102 3400

3

Op Gas Rate (Mscfd)

Op Water Rate (stb/d)

Op Oil Rate (stb/d)

4

5

3000

4

2800

3

2600 2400

2 2200

2 2000 1800

101 101

1600 7 1400

7 5 5

1200

4

4

1000 3

3 2

800 600

2

400 200 1.0

1.0 0 04

05

06

07

08 2008

09

10

11

12

01

02

03

04

05

06

07 2009

08

09

10

11

12

01

02

03

04

05

06

07

08

09

10

2010

23

Casing Pressure (psi(a))

5

3200

Calculated Sandface Pressure (psi(a)) Run Depth Pressure (psi(a))

7

7

Tubing Pressure (psi(a))

102

Dp/q (psi/stb/d)

tetf = SRV

A√k FCD’

Contacted HCPV = 2,800 Mstb

SRV = 850 Mstb

Np/ceDp (Mstb)

q/Dp (stb/d/psi)

Square Root Time

q/Dp (stb/d/psi)

Field Example – RTA – Simplified

Time (d)

24

Field Example – RTA - Modeling High efficiency “short” fracs

Ozkan et al. 2009, “Tri-Linear Flow”

Low efficiency “long” fracs

Stalgorova et al, 2013, “Five Region Model”

Summary of Current Unconventional RTA Technology  Provides a “bulk” reservoir interpretation  Reliable estimation of stimulated and total connected HCPV  Identification of effective system permeability and apparent skin damage  No unique interpretation of fracture properties (orientation, distribution, density, length and conductivity)  No unique interpretation of matrix permeability  Analytical models with different geometries are available

The Future of Unconventional RTA – Probabilistic Approach Rate Transient Analysis: Deterministic

Data q, pwf

Modeling – Realizations of RTA results: Probabilistic 27

Probabilistic Well Performance Analysis  Define ranges or distributions of input parameters  Completion properties  Reservoir properties

 Run the reservoir model probabilistically using Monte Carlo simulation  Keep only history matches that meet a minimum goodness of fit criteria  Report reservoir characteristics and production forecasts as distributions, not single values

28

Probabilistic Well Performance Analysis – Forecasts Rate vs Time

Expected Ultimate Recovery

Rate vs Cumulative

Original Gas in Place

29

Conclusions  RTA provides “bulk” reservoir interpretation  Ideal for establishing connected HCPV  Assists in understanding recovery mechanism  Yields reliable production forecasts

 Analyzing unconventional well production presents significant challenges  Analysis and modeling technology has evolved  Unconventional plays are statistical in nature – many wells must be analyzed to understand behavior  A probabilistic approach will help to manage and communicate uncertainty 30

Thank-you… Questions?

31

Your Feedback is Important Enter your section in the DL Evaluation Contest by completing the evaluation form for this presentation :

Click on:

Section Evaluation

Society of Petroleum Engineers Distinguished Lecturer Program www.spe.org/dl

32