由于地震动存在空间变异性，长大隧道沿纵向的地震反应将有较为显著的行波效应，进行隧道纵向地震反应分析是评价其抗震性能的重要方面。而饱和地基为固液两相体，其动力反应较之单相体更为复杂。本文对饱和地基中隧道在非一致地震作用下纵向反应分析的数值方法进行了较为系统的研究，包括饱和多孔介质波动方程及反应特性的梳理、传输边界的研发、非一致地震输入的三维有限元实现、隧道-地基系统动力反应的数值模拟等。取得的主要研究成果如下：1、对饱和土动力方程及波动传播特点进行了系统的梳理。重点讨论了流体运动方程的建立，指出两相体动力分析时土骨架和孔隙水之间的相互作用包含渗透力和惯性耦合力两项。并进一步澄清，由于惯性耦合力的存在，即使渗透系数无穷大，两类压缩波波速也不会分别等于土骨架和孔隙水单相体中的压缩波波速。2、采用渗透系数为0时柱面波的u-p方程，推导得出了饱和两相介质无限域时域动力有限元分析的高阶弹簧-阻尼-质量传输边界。数值算例表明只要渗透系数的量级满足u-p方程的适用条件，该高阶传输边界均能给出足够精确的结果。将其直接运用于一般二维平面应变内源波动问题时，只要将传输边界布置在离散射波源稍远的位置，所得计算结果仍足以满足工程问题的精度要求。3、通过引入精确的人工边界改进了二维平面内自由波场的一维化时域有限元算法。将一维算法得到的自由场反应扩展至三维，基于一致质量有限元平衡方程，给出了通过边界单元节点的自由场位移求取自由波场产生的等效荷载的公式。采用粘弹性边界，通过在人工边界上施加等效荷载实现波动输入。编制FORTRAN辅助程序并借助有限元软件ANSYS，实现了瞬态平面内体波斜入射至成层半空间场地的三维数值模拟，数值算例验证了其精度和有效性。4、针对饱和地基中长大隧道纵向地震反应数值模拟开展研究。论证了按照非一致地震输入方法计算边界等效荷载，对隧道结构的截断边界不作特殊处理，只需将模型沿隧道纵向取得较长，即可得到满意的精度。研究了体波入射角、隧道刚度对纵向反应的影响，讨论了成层土基中隧道地震反应行波效应的分析方法。5、对非一致地震作用下港珠澳沉管隧道的动力反应进行了较为系统的数值模拟，重点讨论了沉管隧道采用不同接头类型时接头处的内力及变形规律，对实际工程的抗震性能做出评价。