Page tree

You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 23 Next »

Given:

  • a function  q^*(t, {\bf p}) of real-value argument  t \in \R and set of model parameters  {\bf p} = \{ p_m\}_{m = 1..M} = \{p_1, p_2, ... p_M\}
  • a training data set:  \{ (t_k, q_k)\}_{k = 1..N} = \{ (t_0, a_0), (t_1, q_1), ..., (t_N, q_N) \}

the matching procedure assumes searching for thee specific set of model parameters  {\bf p} to minimize the goal function:

F({\bf p}) = \sum_{k=1}^N \, \Psi \left( q^*(t_k) - q_k \right) \rightarrow \textrm{min}

where  \Psi(z) is the discrepancy distance function.

Most popular choices are \Psi(z) = z^2 and \Psi(z) = |z|.


There are few approaches to match the Arps decline to the historical data:



Unconstrained matching

All three model parameters  \{ q^*_0, \, D_0, \, b \} are being varied to achieve the best fit to the training dataset.


Exponential Production DeclineHyperbolic Production DeclineHarmonic Production Decline

b=0

0<b<1

b=1

(1) q(t)=q^*_0 \exp \left( -D_0 \, t \right)
(2) q(t) = \frac{q^*_0}{ \left( 1+b \cdot D_0 \cdot t \right)^{1/b} }
(3) q(t)=\frac{q^*_0}{1+D_0 \, t}


The best-fit model may not match: 

  • the initial production rate  q^*(t=0) \neq q_0
  • the current production rate  q^*(t=t_N) \neq q_N
  • the current cumulative production  Q^*(t=t_N) \neq Q_N


Match the value of the initial rate  q^*(t=0) = q_0

The value of the model rate at the initial time moment is set to training dataset:  q^*(t=0) = q_0 and the two other model properties \{ D_0, \, b \} are being varied to achieve the best fit to the training dataset.

Exponential Production DeclineHyperbolic Production DeclineHarmonic Production Decline

b=0

0<b<1

b=1

(4) q(t)=q_0 \exp \left( -D_0 \, t \right)
(5) q(t) = \frac{q_0}{ \left( 1+b \cdot D_0 \cdot t \right)^{1/b} }
(6) q(t)=\frac{q_0}{1+D_0 \, t}


The best-fit model may not match: 

  • the current production rate  q^*(t=t_N) \neq q_N
  • the current cumulative production  Q^*(t=t_N) \neq Q_N


Match the value of the current rate  q^*(t=t_N) = q_N


To ensure the smooth transition from historical data [(t_1,q_1)... (t_N, q_N)] to the production forecasts in future time moments [(t_{N+1},q_{N+1}), ...] one may wish to constrain the model by firm matching the production at the last historical moment (t_N, q_N) which leads to the following form of Arp's model:

Exponential Production DeclineHyperbolic Production DeclineHarmonic Production Decline

b=0

0<b<1

b=1

(7) q(t)=q_N \cdot \exp \big[ -D_0 \cdot (t-t_N) \big]
(8) q(t) = q_N \cdot \left[ \frac{1+b \cdot D_0 \cdot t_N } { 1+b \cdot D_0 \cdot t } \right]^{1/b}
(9) q(t) = q_N \cdot \left[ \frac{1+D_0 \cdot t_N } { 1+ D_0 \cdot t } \right]
(10) Q(t) - Q_N = [ q_N - q(t)] \, \tau_0
(11) Q(t) - Q_N = \frac{q_N^b \, (\tau_0 + b \, t_N)}{1-b} \left[ q_N^{1-b} - q^{1-b}(t) \right]
(12) Q(t) - Q_N = q_N \, (\tau_0 + t_N) \cdot \ln \frac{q_N}{q(t)}


The best-fit model may not match: 

  • the initial production rate  q^*(t=0) \neq q_0
  • the current cumulative production  Q^*(t=t_N) \neq Q_N

Match the value of the current cumulative  Q^*(t=t_N) = Q_N

To ensure the smooth transition from historical data [(t_1,q_1)... (t_N, q_N)] to the production forecasts in future time moments [(t_{N+1},q_{N+1}), ...] one may wish to constrain the model by firm matching the production at the last historical moment (t_N, q_N) which leads to the following form of Arp's model:

Exponential Production DeclineHyperbolic Production DeclineHarmonic Production Decline

b=0

0<b<1

b=1

(13) q(t)=q_N \cdot \exp \big[ -D_0 \cdot (t-t_N) \big]
(14) q(t) = q_N \cdot \left[ \frac{1+b \cdot D_0 \cdot t_N } { 1+b \cdot D_0 \cdot t } \right]^{1/b}
(15) q(t) = q_N \cdot \left[ \frac{1+D_0 \cdot t_N } { 1+ D_0 \cdot t } \right]
(16) Q(t) - Q_N = [ q_N - q(t)] \, \tau_0
(17) Q(t) - Q_N = \frac{q_N^b \, (\tau_0 + b \, t_N)}{1-b} \left[ q_N^{1-b} - q^{1-b}(t) \right]
(18) Q(t) - Q_N = q_N \, (\tau_0 + t_N) \cdot \ln \frac{q_N}{q(t)}


The best-fit model may not match: 

  • the initial production rate  q^*(t=0) \neq q_0
  • the current production rate  q^*(t=t_N) \neq q_N

Match the value of the current rate and cumulative  q^*(t=t_N) = q_N Q^*(t=t_N) = Q_N


To ensure the smooth transition from historical data [(t_1,q_1)... (t_N, q_N)] to the production forecasts in future time moments [(t_{N+1},q_{N+1}), ...] one may wish to constrain the model by firm matching the production at the last historical moment (t_N, q_N) which leads to the following form of Arp's model:

Exponential Production DeclineHyperbolic Production DeclineHarmonic Production Decline

b=0

0<b<1

b=1

(19) q(t)=q_N \cdot \exp \big[ -D_0 \cdot (t-t_N) \big]
(20) q(t) = q_N \cdot \left[ \frac{1+b \cdot D_0 \cdot t_N } { 1+b \cdot D_0 \cdot t } \right]^{1/b}
(21) q(t) = q_N \cdot \left[ \frac{1+D_0 \cdot t_N } { 1+ D_0 \cdot t } \right]
(22) Q(t) - Q_N = [ q_N - q(t)] \, \tau_0
(23) Q(t) - Q_N = \frac{q_N^b \, (\tau_0 + b \, t_N)}{1-b} \left[ q_N^{1-b} - q^{1-b}(t) \right]
(24) Q(t) - Q_N = q_N \, (\tau_0 + t_N) \cdot \ln \frac{q_N}{q(t)}


The best-fit model may not match: 

  • the initial production rate  q^*(t=0) \neq q_0


See Also


Petroleum Industry / Upstream /  Production / Subsurface Production / Field Study & Modelling / Production Analysis / Decline Curve Analysis / DCA Arps @model




  • No labels