Often in oil reservoirs a layer of water lies under the layer of oil.
The suction pressure due to a distribution of oil wells will cause the
oil-water interface to rise up towards the wells. A three-dimensional
boundary integral formulation is presented for calculating the steady
interface shape when the oil wells are represented by point sinks. So
phisticated integration techniques are implemented in an effort to obt
ain accurate results. In particular, the efficiency of various integra
tion methods are compared for this problem, including QUADPACK routine
s, adaptive methods based on the IMT rule, the Kronrod rule, the metho
d of degenerate quadrilaterals, and the Gauss-Rational rule for infini
te integrals. Numerical results for various general multi-sink distrib
utions are discussed, as are some further results for the axisymmetric
single well problem. (C) 1996 American Institute of Physics.