### Abstract

Darcian steady 2-D flow to a point sink (vertical well) placed eccentrically with respect to two circles demarcating zones of contrasting permeability is studied by the methods of complex analysis and numerically by MODFLOW package. In the analytical approach, two conjugated Laplace equations for a characteristic flow function are solved by the method of images, i.e. the original sink is mirrored about two circles that generates an infinite system of fictitious sinks and source. The internal circle of the annulus models formation damage (gravel pack) near the well and the ring-shaped zone represents a pristine porous medium. On the external circle the head (pressure) is fixed and on the internal circle streamlines are refracted. The latter is equivalent to continuity of pressure and normal component of specific discharge that is satisfied by the choice of the intensity and loci of fictitious sinks. Flow net and dependence of the well discharge on eccentricity are obtained for different annulus radii and permeability ratios. A non-trivial minimum of the discharge is discovered for the case of the ring domain permeability higher than that of the internal circle. In the numerical solution, a finite difference code is implemented and compared with the analytical results for the two-conductivity zone. Numerical solution is also obtained for an aquifer with a three-conductivity zonation. The case of permeability exponentially varying with one Cartesian coordinate within a circular feeding contour is studied analytically by series expansions of a characteristic function obeying a modified Helmholtz equation with a point singularity located eccentrically inside the feeding contour. The coefficients of the modified Bessel function series are obtained by the Sommerfeld addition theorem. A trivial minimum of the flow rate into a small-radius well signifies the trade-off between permeability variation and short-cutting between the well and feeding contour.

Original language | English |
---|---|

Pages (from-to) | 1252-1260 |

Number of pages | 9 |

Journal | Computers and Geosciences |

Volume | 36 |

Issue number | 10 |

DOIs | |

Publication status | Published - Oct 2010 |

### Fingerprint

### Keywords

- Analytical and FDM solution
- Complex potential
- Helmholtz equation
- Laplace equation
- Refraction
- Specific discharge

### ASJC Scopus subject areas

- Information Systems
- Computers in Earth Sciences

### Cite this

*Computers and Geosciences*,

*36*(10), 1252-1260. https://doi.org/10.1016/j.cageo.2010.01.014

**Can heterogeneity of the near-wellbore rock cause extrema of the Darcian fluid inflow rate from the formation (the Polubarinova-Kochina problem revisited)?** / Obnosov, Yurii; Kasimova, Rouzalia; Al-Maktoumi, Ali; Kacimov, Anvar.

Research output: Contribution to journal › Article

*Computers and Geosciences*, vol. 36, no. 10, pp. 1252-1260. https://doi.org/10.1016/j.cageo.2010.01.014

}

TY - JOUR

T1 - Can heterogeneity of the near-wellbore rock cause extrema of the Darcian fluid inflow rate from the formation (the Polubarinova-Kochina problem revisited)?

AU - Obnosov, Yurii

AU - Kasimova, Rouzalia

AU - Al-Maktoumi, Ali

AU - Kacimov, Anvar

PY - 2010/10

Y1 - 2010/10

N2 - Darcian steady 2-D flow to a point sink (vertical well) placed eccentrically with respect to two circles demarcating zones of contrasting permeability is studied by the methods of complex analysis and numerically by MODFLOW package. In the analytical approach, two conjugated Laplace equations for a characteristic flow function are solved by the method of images, i.e. the original sink is mirrored about two circles that generates an infinite system of fictitious sinks and source. The internal circle of the annulus models formation damage (gravel pack) near the well and the ring-shaped zone represents a pristine porous medium. On the external circle the head (pressure) is fixed and on the internal circle streamlines are refracted. The latter is equivalent to continuity of pressure and normal component of specific discharge that is satisfied by the choice of the intensity and loci of fictitious sinks. Flow net and dependence of the well discharge on eccentricity are obtained for different annulus radii and permeability ratios. A non-trivial minimum of the discharge is discovered for the case of the ring domain permeability higher than that of the internal circle. In the numerical solution, a finite difference code is implemented and compared with the analytical results for the two-conductivity zone. Numerical solution is also obtained for an aquifer with a three-conductivity zonation. The case of permeability exponentially varying with one Cartesian coordinate within a circular feeding contour is studied analytically by series expansions of a characteristic function obeying a modified Helmholtz equation with a point singularity located eccentrically inside the feeding contour. The coefficients of the modified Bessel function series are obtained by the Sommerfeld addition theorem. A trivial minimum of the flow rate into a small-radius well signifies the trade-off between permeability variation and short-cutting between the well and feeding contour.

AB - Darcian steady 2-D flow to a point sink (vertical well) placed eccentrically with respect to two circles demarcating zones of contrasting permeability is studied by the methods of complex analysis and numerically by MODFLOW package. In the analytical approach, two conjugated Laplace equations for a characteristic flow function are solved by the method of images, i.e. the original sink is mirrored about two circles that generates an infinite system of fictitious sinks and source. The internal circle of the annulus models formation damage (gravel pack) near the well and the ring-shaped zone represents a pristine porous medium. On the external circle the head (pressure) is fixed and on the internal circle streamlines are refracted. The latter is equivalent to continuity of pressure and normal component of specific discharge that is satisfied by the choice of the intensity and loci of fictitious sinks. Flow net and dependence of the well discharge on eccentricity are obtained for different annulus radii and permeability ratios. A non-trivial minimum of the discharge is discovered for the case of the ring domain permeability higher than that of the internal circle. In the numerical solution, a finite difference code is implemented and compared with the analytical results for the two-conductivity zone. Numerical solution is also obtained for an aquifer with a three-conductivity zonation. The case of permeability exponentially varying with one Cartesian coordinate within a circular feeding contour is studied analytically by series expansions of a characteristic function obeying a modified Helmholtz equation with a point singularity located eccentrically inside the feeding contour. The coefficients of the modified Bessel function series are obtained by the Sommerfeld addition theorem. A trivial minimum of the flow rate into a small-radius well signifies the trade-off between permeability variation and short-cutting between the well and feeding contour.

KW - Analytical and FDM solution

KW - Complex potential

KW - Helmholtz equation

KW - Laplace equation

KW - Refraction

KW - Specific discharge

UR - http://www.scopus.com/inward/record.url?scp=77956738367&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77956738367&partnerID=8YFLogxK

U2 - 10.1016/j.cageo.2010.01.014

DO - 10.1016/j.cageo.2010.01.014

M3 - Article

VL - 36

SP - 1252

EP - 1260

JO - Computers and Geosciences

JF - Computers and Geosciences

SN - 0098-3004

IS - 10

ER -