对球近似的大地测量边值问题,从扰动位的Poisson积分公式出发,推导了大地测量Stokes边值问题的重力异常球谐函数级数表达式,再基于虚功原理,将重力异常球谐函数级数表达式写为变分方程,即大地测量边值问题的自然积分方程,然后再用离散化技术,将所得到的自然积分方程转化为代数线性方程组,通过对代数线性方程组求解,即可得大地测量Stokes边值问题的数值解.