本发明提供了一种球坐标系下重震联合反演密度界面分布的方法，该方法包括以下步骤：估算初始的密度界面深度，计算修正后的密度界面深度#imgabs0#迭代之后重力异常矩阵值#imgabs1#计算，计算已知地震约束点与反演计算点的差值Δhi，均方根误差分析，多次迭代计算及迭代终点判断。本发明的有益效果在于：该方法针对于区域大尺度的密度界面反演研究，将地球近似为球体，用球单元体替代常规的棱柱体，克服了大尺度研究忽略地球曲率的问题，反演结果更接近真实情况，同时反演过程中引入了深度加权因子和地震约束，提高了反演精度，减少了密度界面反演的多解性和不确定性。