The general form of objective function for production targets optimisation is given by:

G(t) = \sum_{p=1}^{N^{\uparrow}_P} \left[ R_O \cdot q^{\uparrow}_{O, p} + R_G \cdot  q^{\uparrow}_{G, p} \right] 
- \sum_{p=1}^{N^{\uparrow}_P} C^{\uparrow}_{L,p} \cdot q^{\uparrow}_{L, p}
- \sum_{p=1}^{N^{\uparrow}_P} C^{\uparrow}_{O,p} \cdot q^{\uparrow}_{O, p} 
- \sum_{p=1}^{N^{\uparrow}_P} C^{\uparrow}_{G,p} \cdot q^{\uparrow}_{G, p} 
- \sum_{p=1}^{N^{\uparrow}_P} C^{\uparrow}_{W,p} \cdot q^{\uparrow}_{W, p}
- \sum_{i=1}^{N^{\downarrow}_W} C^{\downarrow}_{W,j} \cdot q^{\downarrow}_{W, i} 
- \sum_{j=1}^{N^{\downarrow}_G} C^{\downarrow}_{G,j} \cdot q^{\downarrow}_{G, j} \rightarrow \rm max

where

volume/day	oil production rate for -th producer,	cash/volume	cost of produced oil treatment and transportation from wellhead to CMS	cash/volume	oil selling price
volume/day	gas production rate for -th producer,	cash/volume	cost of produced gas treatment and transportation from wellhead to CMS	cash/volume	gas selling price
volume/day	water production rate for -th producer	cash/volume	cost of produced water treatment and transportation from wellhead to CMS	counts	number of producers at
volume/day	liquid production rate for -th producer	cash/volume	cost of fluid lift to the wellhead, cash/volume	counts	number of water injectors at
volume/day	water injection rate for -th water injector	cash/volume	cost of water injection, including purchase, treatment, transportation and pumping	counts	number of gas injectors at
volume/day	gas injection rate for -th gas injector	cash/volume	cost of gas injection, including purchase, treatment, transportation and pumping	months	time

This can be rewritten in terms of sandface flowrates:

G = \sum_{p=1}^{N^{\uparrow}_P} G^{\uparrow}_{t,p} \cdot q^{\uparrow}_{t, p}
- \sum_{i=1}^{N^{\downarrow}_W} G^{\downarrow}_w  \cdot 
q^{\downarrow}_{w, i}
- \sum_{j=1}^{N^{\downarrow}_G} G^{\downarrow}_g  \cdot q^{\downarrow}_{g, j} \rightarrow \rm max

where

G(t) = \sum_{p=1}^{N^{\uparrow}_P} \left[ 
(R_O -  C^{\uparrow}_{O,p}) \cdot q^{\uparrow}_{O, p} + (R_G - C^{\uparrow}_{G,p}) \cdot  q^{\uparrow}_{G, p} 
- C^{\uparrow}_{L,p} \cdot q^{\uparrow}_{L, p} - C^{\uparrow}_{W,p} \cdot q^{\uparrow}_{W, p}
\right]  
- \sum_{i=1}^{N^{\downarrow}_W} C^{\downarrow}_{W,j} \cdot q^{\downarrow}_{W, i} 
- \sum_{j=1}^{N^{\downarrow}_G} C^{\downarrow}_{G,j} \cdot q^{\downarrow}_{G, j}

G(t) = \sum_{p=1}^{N^{\uparrow}_P} \left[ 
\left[  (R_O -  C^{\uparrow}_{O,p}) + (R_G - C^{\uparrow}_{G,p}) \cdot  Y_{g,p} \right]  \cdot q^{\uparrow}_{O, p} 
- C^{\uparrow}_{L,p} \cdot q^{\uparrow}_{L, p} - C^{\uparrow}_{W,p} \cdot Y_{w,p} \cdot q^{\uparrow}_{L, p}
\right]  
- \sum_{i=1}^{N^{\downarrow}_W} C^{\downarrow}_{W,j} \cdot q^{\downarrow}_{W, i} 
- \sum_{j=1}^{N^{\downarrow}_G} C^{\downarrow}_{G,j} \cdot q^{\downarrow}_{G, j}

G(t) = \sum_{p=1}^{N^{\uparrow}_P} \left[ 
\left[  (R_O -  C^{\uparrow}_{O,p}) + (R_G - C^{\uparrow}_{G,p}) \cdot  Y_{g,p} \right]  \cdot (1- Y_{w,p}) 
- C^{\uparrow}_{L,p} \cdot q^{\uparrow}_{L, p} - C^{\uparrow}_{W,p} \cdot Y_{w,p} 
\right]  \cdot q^{\uparrow}_{L, p}
- \sum_{i=1}^{N^{\downarrow}_W} C^{\downarrow}_{W,j} \cdot q^{\downarrow}_{W, i} 
- \sum_{j=1}^{N^{\downarrow}_G} C^{\downarrow}_{G,j} \cdot q^{\downarrow}_{G, j}

Substituting liquid production rate from one arrives to:

G(t) = \sum_{p=1}^{N^{\uparrow}_P}  \frac{\left[  (R_O -  C^{\uparrow}_{O,p}) + (R_G - C^{\uparrow}_{G,p}) \cdot  Y_{g,p} \right]  \cdot (1- Y_{w,p}) 
- C^{\uparrow}_{L,p} \cdot q^{\uparrow}_{L, p} - C^{\uparrow}_{W,p} \cdot Y_{w,p} }
{B_w Y_{w,p} + \left[ (B_o - R_s B_g] + (B_g - R_v B_o) \, Y_{g,p} \right] \cdot (1-Y_{w,p})}
 
 \cdot q^{\uparrow}_{t, p}
- \sum_{i=1}^{N^{\downarrow}_W} C^{\downarrow}_{W,j} \cdot B_w \cdot q^{\downarrow}_{w, i} 
- \sum_{j=1}^{N^{\downarrow}_G} C^{\downarrow}_{G,j} \cdot B_g \cdot q^{\downarrow}_{g, j}

See Also