我国幅员辽阔,地形地貌多变,自然地质灾害种类众多。滑坡作为一种最为常见的地质灾害之一,给我国的工程建设、人民的人身安全和财产安全(尤其是西南部、西北部地区)带来了严重的危害和损失。因此,滑坡的稳定性受到了广泛的关注。由于岩土体本身具有一定的变异性,传统地按照确定性的方法计算滑坡稳定性系数很多时候无法满足需求。于是工程师们试图从概率的角度来描述和刻画滑坡的稳定性系数,从而判断滑坡的失稳概率。然而,这需要人们能够较为准确的获得计算时所需的强度参数本身的概率分布。但在实际工程中,由于工程本身规模和经费限制,尤其对于中小型工程,工程师很难得到大量的关于岩土性质试验测量数据,进而获得相关参数的统计特征。这使得这一思想地实施非常困难,而若根据这些稀少数据获得相关参数平均值和标准差,简单假设参数本身服从正态分布从而得出计算结果,又有可能带来一定的或者较大的计算误差。本文提出一种数学方法,该方法利用贝叶斯理论以及全概率公式合理利用现场稀少的试验数据以及工程经验,推断出具体岩土性质的概率分布。但由于其概率分布为隐式表达,很难直观观察出其内在的逻辑关系。本文采用马尔科夫链蒙特卡洛(MCMCS),将该分布作为平稳分布,产生大量的等效样本。这些等效样本本身的统计特征即反应了该隐式分布函数的内在特征,很好的解决了这一问题。利用该方法,对坪子山滑坡的土体性质进行了概率表征,进而得到了坪子山滑坡不同剖面在不同工况下的失稳概率。此外,本文还利用贝叶斯理论进一步确定导致该滑坡在不同工况下失稳的主要因素,从而为该滑坡的防治提供理论指导依据。