布格重力异常是由自由空气重力异常作地形校正和中间层校正而得,近地表密度是这两个校正的重要参数。近地表密度的直接获取方法是采集岩石标本并进行密度测定,但在某些区域（比如沉积覆盖区）标本采集条件并不充分,这就需要间接估计近地表密度的方法。此外,当研究区面积较大时候,近地表密度往往不再是常密度,此时若按照常规办法采用常密度进行地形和中间层校正,将对得到的布格重力异常产生不可忽略的偏差。本文研究一种基于重力贝叶斯分析的近地表密度估计方法,用于获取区域近地表横向密度分布和光滑布格重力异常。本文方法由自由空气重力异常观测数据和近地表密度、布格重力异常模型参数构建反演目标函数,同时基于布格重力异常相对地表地形起伏光滑分布的先验假设,对布格重力异常模型进行三次B样条光滑约束,采用基于ABIC准则的贝叶斯最优后验估计算法迭代求解。理论一维、二维模型数据试验验证了本文方法的有效性。相关参数影响评价表明,足够多的样条系数可以提高估计结果的准确性;当异常的噪声水平变化或模型粗糙程度变化时,ABIC准则可有效地自动调整先验光滑约束的权重。云南地区云县-会东、思茅-七甸两条重力剖面应用结果表明,剖面沿线的近地表密度值起伏变化明显,达2.45².8g·cm^-3,前寒武纪和古生代地层密度相对较高,主要为2.53².75g·cm^-3,而中生代密度较低,为2.45².73g·cm^-3;红河断裂两侧近地表密度差异较大,可达0.4g·cm^-3。云南地区的二维区域应用结果表明,研究区近地表密度值大致分布在2.3².75g·cm^-3之间,整体分布特征为北部区域和东部区域密度值略高于西部和南部区域;前寒武纪密度分布在2.55².65g·cm^-3,古生代密度分布在2.5².65g·cm^-3,部分区域有达到2.7g·cm^-3以上的密度值,中生代密度分布较广,在2.45².7g·cm^-3之间。本文估计的近地表密度结果与区域物性资料及地表地质特征大致吻合,估计的布格重力异常具有良好的光滑性。本文方法为区域近地表密度结构估计提供技术支撑,为区域深部结构与构造研究提供更可靠的布格重力异常数据。