...
| Expand |
|---|
|
Derivation of Radial VEH Aquifer Drive @model| Panel |
|---|
|
| borderColor | wheat |
|---|
| borderWidth | 10 |
|---|
|
123
Transient flow in Radial Composite Reservoir:
| LaTeX Math Block |
|---|
| \frac{\partial p_a}{\partial t} = \chi \cdot \left[ \frac{\partial^2 p_a}{\partial r^2} + \frac{1}{r}\cdot \frac{\partial p_a}{\partial r} \right] |
| | LaTeX Math Block |
|---|
| p_a(t = 0, r)= p(0) |
| | LaTeX Math Block |
|---|
| p_a(t, r=r_e) = p(t) |
| | LaTeX Math Block |
|---|
| anchor | p1_PSS |
|---|
| alignment | left |
|---|
| \frac{\partial p_a}{\partial r}
\bigg|_{(t, r=r_a)} = 0 |
|
Consider a pressure convolution:
| LaTeX Math Block |
|---|
| p_a(t, r) = p(0) + \int_0^t p_1 \left(\frac{(t-\tau) \cdot \chi}{r_e^2}, \frac{r}{r_e} \right) \dot p(\tau) d\tau |
| | LaTeX Math Block |
|---|
| \dot p(\tau) = \frac{d p}{d \tau} |
|
One can easily check that
| LaTeX Math Block Reference |
|---|
|
honors the whole set of equations | LaTeX Math Block Reference |
|---|
|
–| LaTeX Math Block Reference |
|---|
|
and as such defines a unique solution of the above problem.Water flowrate within
sector angle at interface with oil reservoir will be:| LaTeX Math Block |
|---|
|
q^{\downarrow}_{AQ}(t)= \theta \cdot r_e \cdot h \cdot u(t,r_e) |
where
is flow velocity at aquifer contact boundary, which is:| LaTeX Math Block |
|---|
|
u(t,r_e) = M \cdot \frac{\partial p_a(t,r)}{\partial r} \bigg|_{r=r_e} |
where
| LaTeX Math Inline |
|---|
| body | --uriencoded--\displaystyle M = \frac%7Bk_w%7D%7B\mu_w%7D |
|---|
|
is aquifer mobility.Water flowrate becomes:
| LaTeX Math Block |
|---|
|
q^{\downarrow}_{AQ}(t)= \theta \cdot r_e \cdot h \cdot M \cdot \frac{\partial p_a(t,r)}{\partial r} \bigg|_{r=r_e} |
Cumulative water flux:
| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \int_0^t q^{\downarrow}_{AQ}(t) dt = \theta \cdot r_e \cdot h \cdot M \cdot \int_0^t \frac{\partial p_a(t,r)}{\partial r} \bigg|_{r=r_e} dt |
Substituting
| LaTeX Math Block Reference |
|---|
|
into | LaTeX Math Block Reference |
|---|
|
leads to:| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot r_e \cdot h \cdot M \cdot \int_0^t d\xi \ \frac{\partial }{\partial r} \left[
\int_0^\xi p_1 \left( \frac{(\xi-\tau)\chi}{r_e^2}, \frac{r}{r_e} \right) \, \dot p(\tau) d\tau
\right]_{r=r_e} |
| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot h \cdot M \cdot \int_0^t d\xi \ \frac{\partial }{\partial r_D} \left[
\int_0^\xi p_1 \left( \frac{(\xi-\tau)\chi}{r_e^2}, r_D \right) \, \dot p(\tau) d\tau
\right]_{r_D=1} |
| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot h \cdot M \cdot \int_0^t d\xi \
\int_0^\xi \frac{\partial p_1}{\partial r_D} \left( \frac{(\xi-\tau)\chi}{r_e^2}, r_D \right) \Bigg|_{r_D=1} \, \dot p(\tau) d\tau
|
The above integral represents the integration over the
area in plane (see Fig. 1):| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot h \cdot M \cdot \iint_D d\xi \ d\tau \, \dot p(\tau)
\frac{\partial p_1}{\partial r_D} \left( \frac{(\xi-\tau)\chi}{r_e^2}, r_D \right) \Bigg|_{r_D=1}
|
 Image Removed
|
Fig. 1. Illustration of the integration area in plane |
Changing the integration order from
to leads to:| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot h \cdot M \cdot \int_0^t d\tau \int_\tau^t d\xi \ \dot p(\tau)
\frac{\partial p_1}{\partial r_D} \left( \frac{(\xi-\tau)\chi}{r_e^2}, r_D \right) \Bigg|_{r_D=1}
=
\theta \cdot h \cdot M \cdot \int_0^t \dot p(\tau) d\tau \int_\tau^t d\xi \
\frac{\partial p_1}{\partial r_D} \left( \frac{(\xi-\tau)\chi}{r_e^2}, r_D \right) \Bigg|_{r_D=1} |
Replacing the variable:
| LaTeX Math Block |
|---|
|
\xi = \tau + \frac{r_e^2}{\chi} \cdot t_D \rightarrow t_D = \frac{(\xi-\tau)\chi}{r_e^2} \rightarrow d\xi = \frac{r_e^2}{\chi} \cdot dt_D |
and flux becomes:
| LaTeX Math Block |
|---|
|
Q^{\downarrow}_{AQ}(t) = \theta \cdot h \cdot M \cdot \frac{r_e^2}{\chi} \cdot \int_0^t \dot p(\tau) d\tau \int_0^{(t-\tau)\chi/r_e^2}
\frac{\partial p_1( t_D, r_D)}{\partial r_D} \Bigg|_{r_D=1} dt_D = B \cdot \int_0^t \dot p(\tau) d\tau \int_0^{(t-\tau)\chi/r_e^2}
\frac{\partial p_1( t_D, r_D)}{\partial r_D} \Bigg|_{r_D=1} dt_D |
where
is water influx constant and which leads to | LaTeX Math Block Reference |
|---|
|
and | LaTeX Math Block Reference |
|---|
|
.qwe
Computational Model
...
| LaTeX Math Block |
|---|
| Q^{\downarrow}_{AQ}(t)= B \cdot \sum_\alpha W_{eD}
\left( \frac{ (t-\tau_\alpha) \chi}{r_e^2}, \frac{r_a}{r_e} \right)\Delta p_\alpha
= B \cdot W_{eD}
\left( \frac{ (t-\tau_1) \chi}{r_e^2}, \frac{r_a}{r_e} \right)\Delta p_1 +
B \cdot W_{eD}
\left( \frac{ (t-\tau_2) \chi}{r_e^2}, \frac{r_a}{r_e} \right)\Delta p_2
+ ... + B \cdot W_{eD}
\left( \frac{ (t-\tau_N) \chi}{r_e^2}, \frac{r_a}{r_e} \right)\Delta p_N |
|
...
