地球表面负荷理论在研究海潮负荷形变、冰川负荷形变以及地下水负荷形变等诸多领域有着广泛的应用。传统理论在求解负荷形变边值问题时中往往采用龙格库塔数值积分方法，这种方法在求解过程中存在“球心奇异性”以及“高阶不稳定性”等缺陷。本研究旨在提出一种基于谱元法的数值新方法，用于解决龙格库塔方法的诸多缺陷，实现对地球表面负荷形变的高效高精度求解。推导了真实地球形变边值问题的弱形式解，发展了基于球谐函数以及高阶样条函数为基函数的谱元法，计算了表面负荷形变问题的格林函数。相关计算结果与解析解比较，显示出较好的收敛性和准确性。