实际中的土体大多处于非饱和状态,目前针对非饱和土渗流问题的研究,较多采用的是流-固或气-液两相耦合模型,较少涉及同时考虑固、液、气的三相耦合模型。另外,在研究中,为了方便,主要针对的是均质非饱和土,而事实上由于岩土体的形成过程十分复杂,它的性质往往表现出明显的空间变异性,而非饱和土在考虑空间变异性条件下的多场耦合分析则更为复杂。研究非饱和土的参数空间变异性对渗流-变形及稳定性问题的影响,对实际工程具有重要理论意义和指导作用。本文根据严格的理论推导,得出了固-液-气三相耦合控制方程组,并将三相耦合模型与考虑岩土参数空间变异性的非高斯随机场结合,利用蒙特卡洛法对非饱和土的降雨入渗以及稳定性问题进行了分析,主要工作内容和成果如下:1)建立了非饱和土的固-液-气三相耦合的控制方程组,并引入COMSOL中求解计算。利用两个著名的非饱和土固结试验验证了所推导三相耦合模型的正确性和可靠性。对比了三相耦合模型和两相耦合模型之间的差异,并详细解释了固结过程中,气体的迁移、液体的渗流、土骨架的变形以及它们之间的相互作用。2)采用基于Cholesky分解的中点法对岩土参数随机场进行离散,并使其服从非高斯分布。将岩土参数的空间变异性理论与固-液-气三相耦合模型相结合,对一维非饱和土柱的降雨入渗问题进行分析。对土柱中的水分渗流、气体迁移、土骨架变形进行了研究,利用蒙特卡洛法,详细分析了弹性模量(E)、土-水特征曲线参数(a、n)和饱和渗透系数(K_s)的空间变异性对结果的影响。结果表明,它们的相关距离和变异系数对竖向位移、孔隙水压力、孔隙气压力的均值影响不大,但对标准差影响显著;a和K_s、n和K_s的互相关系数对竖向位移、孔隙水压力、孔隙气压力的均值基本无影响,但对标准差却有着显著影响。3)将一维三相耦合随机模型扩展至二维,采用强度折减法对堤坝稳定性问题进行研究,重点分析了粘聚力(c)、内摩擦角(?)、饱和渗透系数(K_s)、弹性模量(E)的空间变异性对堤坝稳定性系数均值和可靠指标的影响;另外还对降雨边坡模型试验进行了研究,探讨了K_s的空间变异性对边坡降雨过程的影响。结果表明,c和?的水平相关距离对稳定性系数均值基本无影响,但竖向相关距离对其则有较大影响,尤其在变异系数取值高时显著;K_s的变异系数和相关距离对可靠指标有着重要影响;透水性较大夹层的分布对稳定性系数有着重要的影响,当透水性较大夹层位于中下部时,不利于堤坝稳定,而当透水性较大夹层位于中上部时,有利于堤坝稳定;在二维问题中,K_s的相关距离和变异系数对孔隙水压力、孔隙气压力、水平位移和竖向位移的均值和标准差的影响与一维问题类似,但也有不同之处。