采用等效连续介质模型,引入介质渗透系数张量与介质体积应变之间的函数关系式,建立能反映介质中孔隙压力与骨架变形非线性耦合作用的数学模型。视模型中的参数为随机变量,采用Taylor展开随机有限元法对该模型进行数值模拟,研究参数具有随机性情况下耦合场中孔隙压力、有效应力、位移等场变量的随机分布特征及孔隙压力、位移等场变量对参数的灵敏度。数值算例可检验所提出方法的有效性。