泥岩是地下工程建设中比较常见的一种软弱岩体,具有明显的流变特性,随着资源开采以及地下工程建设向深部的转移,泥岩在多场耦合作用下的破坏机制及其对工程长期稳定性的影响已成为工程界急需解决的关键技术。泥岩由于其低渗透性、良好的蠕变性和遇水损伤自修复特征被认为是一种良好的核废料储存介质,法国、比利时、瑞士等国相继建成地下试验平台,开展泥岩现场多场耦合试验研究。结合比利时地下大型试验室正开展的泥岩研究课题,本文从泥岩的力学特性、渗透性、开挖扰动区、地下水-围岩相互作用机理以及围岩流变性等方面完成如下工作:(1)建立了考虑最大拉应力准则的修正Mohr-Coulomb模型,采用向后欧拉隐式应力积分算法编制了UMAT本构程序;为了描述泥岩的硬化和软化行为,将损伤引入到修正的Mohr-Coulomb准则中,基于损伤势函数的概念建立了泥岩弹塑性损伤本构模型,导出了泥岩的损伤演化方程,编制了泥岩弹塑性损伤本构程序。(2)基于拉丁超立方抽样,提出了采用非参数统计方法中的秩相关系数来评价多因素敏感性的方法;分别建立了基于侧压力系数的三维地应场反演模型、位移量测反演模型和渗流场反演模型,提出Nelder-Mead法与有限元联合反演法,将有限元程序作为一个单独模块嵌入到Nelder-Mead算法程序中,以测点的实测值与计算值建立目标函数,采用精确罚函数法实现有约束的反演问题转化成无约束的反演问题,编制了优化反演分析程序;(3)针对本构模型的参数辨识问题,编制了本构模型参数反演程序,并根据非排水条件下泥岩三轴试验结果,采用多目标函数优化反分析法获得了泥岩本构模型参数。(4)在分析软岩与水相互作用机理、不同应力-应变阶段渗透性演化规律以及隧洞围岩渗透性分区的基础上,建立了渗透系数、孔隙度等参数的动态演化方程,导出了岩土介质的渗流.应力动态全耦合模型;基于应力应变.渗透率全过程试验和渗透性工程试验,通过引入损伤的概念,建立了描述岩体破坏过程的渗流-应力-损伤耦合模型;对隧洞围岩裂隙自愈合机理进行了分析,通过引入愈合应力和水化学愈合因子的概念,建立了描述泥岩裂隙自愈合特性的渗透性自愈合模型;基于现场监测的孔隙压力资料,采用有限元优化法反演了泥岩的渗透系数和渗透性演化方程中的待定参数。(5)根据泥岩的非线性蠕变变形特点,建立了蠕变损伤与蠕应变的关系式,构造了基于Mohr-Coulomb准则的蠕变势,导出了考虑渗流-应力-损伤耦合的非线性蠕变损伤本构模型,采用了显式积分算法编制了uMAT蠕变本构程序;基于近20年的现场监测成果,采用位移反分析法获得了本构模型中的待定参数。(6)在泥岩力学特性、渗透性、本构模型及长期强度准则等研究成果的基础上,采用数值模拟方法,研究储库围岩孔隙压力分布规律、开挖扰动区的损伤演化过程及渗透性演化规律,对围岩稳定性进行评价,预测试验巷道长期稳定性。