关于《基于顺序耦合方法的多年冻土融化对埋地管道力学响应研究》的学术报告
本研究的主要作者为来自江苏大学土木工程与力学学院的王飞、中国科学院西北生态环境资源研究院冻土工程国家重点实验室的吴刚、陈敦（通讯作者）、李国玉、钱玉龙、席飞龙和王玲。该研究于2023年3月24日发表在学术期刊《Atmosphere》（第14卷，第620期）上。
一、 研究的学术背景
本研究属于寒区工程与岩土工程交叉领域，具体聚焦于冻土力学与管道结构工程。研究的开展源于一个迫切的现实问题：在气候变化背景下，多年冻土（Permafrost）的持续退化和融化，对穿越寒区、特别是富含冰多年冻土区的埋地温输油管道（如中俄原油管道、阿拉斯加管道等）的结构安全构成了长期且日益严重的威胁。管道运营散发的热量会扰动其下方及周围的冻土热平衡，导致冻土融化、地基发生不均匀融沉（Thaw Settlement），进而引发管道产生过大的附加应力和变形，可能导致管道失效，引发严重的环境与安全事故。
尽管前人已在冻土-管道热相互作用、以及管道在冻胀（Frost Heave）条件下的响应方面开展了大量研究，但对于管道在多年冻土融化沉降下的力学行为，特别是考虑热-力（Thermo-Mechanical）顺序耦合过程的深入研究相对有限。现有研究多将冻土融化与管道力学响应作为两个独立过程处理，或采用简化的模型，未能充分反映冻土融化固结（Thaw Consolidation）这一复杂的热-水-力（Thermo-Hydro-Mechanical, THM）耦合过程对管道受力的真实影响。
因此，本研究旨在：1）验证一种顺序热-力耦合数值方法用于预测冻土融化沉降的适用性与准确性；2）基于中俄原油管道沿线实际监测站点的数据，建立埋地管道-冻土相互作用的三维热-力耦合模型，模拟分析在差异融沉（Differential Thaw Settlement）地基条件下，管道的力学响应（变形与应力）随时间的发展规律；3）评估管道在30年运营期内的安全性，并为类似寒区管道工程的设计、维护和风险评估提供理论依据和参考。
二、 详细研究流程
本研究主要包含两大核心流程：顺序耦合方法的验证与管道-冻土相互作用的全尺寸三维模拟。
流程一：顺序热-力耦合方法的验证 为了确保后续管道-冻土相互作用分析的可靠性，研究首先对一个一维冻土融化固结试验进行了数值模拟复现与验证。 1. 研究对象与模型建立：研究对象为 Huang (2021) 报告中所述的粉质粘土试样（直径101毫米，高100毫米，干密度1.6克/立方厘米）。研究团队根据试验条件，建立了一个简化的轴对称数值模型。 2. 边界条件与参数设置：模型严格复现了试验的周期性热边界（试样顶部温度按正弦函数变化，模拟冻融循环）和力学边界（顶部施加100 kPa恒定荷载，底部固定温度并约束位移）。土壤的热物理参数（导热系数λ、体积热容c）和力学参数（弹性模量E、泊松比μ、渗透系数k等）均采用先前已发表的试验数据（Xu等，2001；Yao等，2012），具体数值在文中列表给出。 3. 数值方法：采用顺序热-力耦合方法。第一步，对整个模型进行纯热传导计算，求解考虑冰水相变的瞬态温度场（控制方程如原文公式(1)）。第二步，将第一步计算得到的温度场结果作为输入，仅在已融化的区域进行大应变固结计算（控制方程包括平衡方程(2)、几何方程(3)、本构方程(4)、有效应力原理(5)和连续性方程(6)）。这种方法将复杂的热-水-力耦合问题解耦为依次进行的热分析和力学分析，显著降低了计算复杂度。 4. 验证流程：运行模型，计算试样在不同深度处的温度变化以及顶部的竖向变形（融沉量） 随时间的发展过程。将模拟结果与Huang (2021)的实测数据进行对比。
流程二：埋地管道-冻土相互作用的三维热-力耦合模拟 在验证了方法可靠性后，研究建立了基于中俄原油管道实际工况的全尺寸三维模型。 1. 研究对象与计算域：研究对象为中俄原油管道沿线一个温暖、富冰多年冻土区监测断面的管道-地基系统。管道采用X65钢材，直径和壁厚依据工程实际设定（文中提及壁厚16mm的案例）。地基土体根据地质资料简化为三层：粉质粘土、砾石土和风化花岗岩。计算域总长120米（沿管道轴向），宽60米（从管道中心线向两侧各30米），深30米。为了模拟差异融沉，模型沿管道轴向（Z轴）被等分为三段，中间40米为强融沉区（赋予粉质粘土和砾石土较低的力学参数和较高的渗透系数），两端各40米为弱融沉区（赋予相同的土层但较高的力学参数和较低的渗透系数）。 2. 模型构建与网格：管道采用四节点减缩积分壳单元模拟，周围土体采用八节点减缩积分实体单元模拟。这种单元选择兼顾了计算精度与效率。 3. 边界条件设定： * 热边界：模型侧面假设为绝热；底部施加地热通量（0.01 °C/m）；顶部地表采用对流换热边界，大气温度基于加格达奇气象站的长期观测数据，拟合为考虑气候变暖趋势（未来50年升温2.4°C）的正弦函数（原文公式(7)）。管道内壁与流动原油之间也采用对流换热边界，原油温度依据泵站监测数据拟合为正弦函数（原文公式(8)）。初始地温场通过不考虑气候变暖的长期瞬态解获得。 * 力学边界：模型顶部自由位移且孔隙水压力为零（自由排水）；底部边界固定；四个侧面约束法向位移。 4. 材料参数：所有土层的热物理参数（冻/融状态下的λ, c, 相变潜热L）均基于现场测量和室内试验结果。冻土和融土的力学行为采用Mohr-Coulomb本构模型描述，其参数（E, μ, 粘聚力c, 内摩擦角φ）引自先前文献（Li等，2009；Wen等，2010）。管道钢材的弹塑性行为采用Ramberg-Osgood本构模型描述（原文公式(9)），参数包括弹性模量（207 GPa）、屈服强度（450 MPa）和材料常数（α, n）。 5. 模拟流程：首先，进行长达30年的瞬态热分析，模拟在气候变暖和管道散热的共同作用下，管道周围地基土温度场的演变以及融化盘（Thaw Bulb）的形成与发展。然后，将每年（或关键时间点）的温度场结果顺序导入力学分析步，计算由冻土融化和融土固结引起的地基沉降，进而求解管道在此沉降作用下的应力与变形响应。重点分析了管道运营第2、5、10、20、30年的结果。
三、 主要研究结果
流程一的结果：如图2所示，模拟得到的试样不同深度处温度变化趋势与实测数据一致，但在振幅和相位上存在一定差异，作者认为这可能是由于对相变潜热处理方式的简化所致。更重要的是，模拟得到的试样竖向变形（融沉）随时间呈阶梯式增长的趋势与试验观测结果高度吻合（图2b）。这有力地证明了所采用的顺序热-力耦合方法能够有效地预测周期性边界条件下冻土的融化沉降行为，为后续大规模三维模拟提供了方法学信心。
流程二的结果： 1. 管道周围地温场演变：模拟结果（图5）清晰显示，管道运营显著扰动了原始平直的地温等值线，使其在管道周围呈“凹陷”状分布。随着运营时间增长，温度梯度减小，地温场趋于均匀。管道下方的多年冻土上限（Permafrost Table） 持续下降，其深度随时间的变化可用幂函数（y=3.45x^0.34）描述（图7）。在前5年，冻土上限下降速率高达0.9米/年，随后逐渐减缓。在运营10-15年后，融化锋面已延伸至风化花岗岩层（埋深约8米以下），该层冰含量低、融沉性弱。 2. 融化盘形成：即使在冬季，管道上方的季节冻结深度也小于天然场地，导致在多年冻土上限与季节冻结层底之间形成残留融化层（即融化盘），且该融化盘随运营时间不断向周围和深处扩展（图5, 6）。 3. 管道沉降发展：位于强融沉区中部的管道竖向沉降在运营前10年快速发展，随后20年趋于稳定，最终沉降量在33-36厘米之间（图8）。沉降趋于稳定的原因正是融化锋面进入了沉降量很小的风化花岗岩层。这突显了地基土融沉特性对管道沉降的显著影响。 4. 管道应力响应：差异融沉在管道中引起了显著的应力集中。椭圆形的高应力区出现在强、弱融沉区交界处附近管道的顶部和底部（图9）。随着地基沉降发展，这些高应力区范围扩大。 * 具体数据：管道运营第2、5、10年，最大拉应力分别为36 MPa、79 MPa和145 MPa；最大压应力分别为35 MPa、73 MPa和133 MPa。 * 安全评估：这些应力值均小于X65钢材的许用应力（324 MPa，基于450 MPa屈服强度和0.72安全系数计算得出）。在模拟条件下（强、弱融沉区长度比为1:1，管壁厚16 mm），管道在30年内的最大应力约为其许用应力的45%，因此被认为在此时间段内是安全的（图10）。 5. 高风险位置：应力分析明确指出，管道在强、弱融沉区交界面的两侧应力最大，是潜在的失效高风险位置。
四、 研究结论与意义
本研究得出以下主要结论： 1. 通过一维融沉固结试验对比验证，顺序热-力耦合数值方法能够可靠地预测周期性边界条件下退化冻土的融沉固结行为。 2. 埋地温油管道会显著扰动周围地温场，导致多年冻土上限持续下降、融化盘形成并扩展。管道下方冻土上限深度随运营时间呈幂函数形式增加。 3. 管道沉降随地基冻土融化而增加，当融化锋面进入弱融沉的风化花岗岩层后趋于稳定。地基差异融沉会在管道强、弱融沉区交界面附近产生椭圆形的高应力集中区。 4. 在本研究设定的参数和30年模拟期内，管道应力水平低于材料许用应力，结构安全。但作者强调，若考虑模型中未包含的积水热侵蚀、管道持续下移以及热-力完全耦合效应（即沉降对热传导的反饋），实际冻土退化速度和地基位移可能被低估，从而导致管道应力被低估。此外，管道材料的老化也会增加长期风险。
本研究的价值体现在： * 科学价值：发展并验证了一种用于分析冻土-管道相互作用复杂问题的顺序热-力耦合数值框架，深化了对多年冻土区埋地管道在热扰动下长期力学行为演化规律的理解，特别是量化了差异融沉导致的管道应力集中效应。 * 应用价值：为中俄原油管道等寒区在役管道的安全评估、风险预测和维护决策提供了重要的分析工具和理论依据。研究指出的高风险位置（差异沉降交界区）可为管道重点监测和加固段落的选取提供指导。同时，研究也暗示了在管道规划阶段进行详细冻土工程地质勘察的重要性。
五、 研究亮点
六、 其他有价值的补充
文章在讨论部分坦诚指出了本研究的两个主要局限性：一是未考虑积水热侵蚀和管道下沉对热边界条件的反馈，可能低估长期融化；二是采用了解耦的热-力分析，若考虑融沉固结对热传导过程的影响，可能会产生更大的融化盘。这些局限性可能导致对管道应力的估计偏于保守（即低估）。作者因此提出，未来的研究应致力于建立考虑冻土热-水-力学完全耦合过程以及管-土接触特性的更全面模型，并进行参数化分析，以更精确地评估寒区管道的结构完整性。这一讨论为后续研究指明了方向。