利用边界积分方程方法可以高效计算界面Stokes流动。在三维情况下，目标点处的流体速度由对所有界面的二维表面积分给出，从而降低了问题的维度。一个核心挑战是，当目标点靠近但不在界面时，表面积分是近奇异的，标准的数值积分方法会失去精度。本文提出了一种方法，能够准确计算三维椭圆边值问题中出现的近奇异积分。该方法基于通过目标点在表面上的正交投影得到的基点处对被积函数的局部级数展开式。由此得到的级数展开式中的初等函数可以在基点附近的区域内利用递归算法实现高精度积分。剩余的积分则使用标准的数值积分规则进行计算，这里选择的是4阶梯形法则。该方法可简化为标准积分加上修正项，并且在整个范围内保持4阶精度。该方法被应用于解决几个椭球刚体周围的Stokes流动。我们比较了在刚体附近的速度误差，以及在有无修正的情况下计算的粒子围绕刚体运动的时间和位移误差。

Interfacial Stokes flow can be efficiently computed using the Boundary Integral Equation method. In 3D, the fluid velocity at a target point is given by a 2D surface integral over all interfaces, thus reducing the dimension of the problem. A core challenge is that for target points near, but not on, an interface, the surface integral is near-singular and standard quadratures lose accuracy. This paper presents a method to accurately compute the near-singular integrals arising in elliptic boundary value problems in 3D. It is based on a local series approximation of the integrand about a base point on the surface, obtained by orthogonal projection of the target point onto the surface. The elementary functions in the resulting series approximation can be integrated to high accuracy in a neighborhood of the base point using a recursive algorithm. The remaining integral is evaluated numerically using a standard quadrature rule, chosen here to be the 4th order Trapezoidal rule. The method is reduced to the standard quadrature plus a correction, and is uniformly of 4th order. The method is applied to resolve Stokes flow past several ellipsoidal rigid bodies. We compare the error in the velocity near the bodies, and in the time and displacement of particles traveling around the bodies, computed with and without the corrections.