智能海底维护：蜡与法兰管理的一体化方法

《Chemical Thermodynamics and Thermal Analysis》：Smart subsea maintenance: An integrated approach to wax and flange management

　　

当原油温度低于蜡出现温度（WAT）时，原油输送管道中会发生蜡结晶。这给石油和天然气行业带来了关键的流动保障风险，常见后果包括有效管径受限、泵送成本增加以及通过压力尖峰对结构完整性造成威胁，最坏情况下会导致管道完全堵塞。研究人员开发了一个集成数学模型，利用Bonny Light原油数据模拟40公里管道中的蜡积累。该框架耦合了热力学溶解度理论（Van't Hoff）、沉淀动力学、成核动力学和Fick传质，以及改进的Darcy-Weisbach压降分析和符合ASME B16.5 Class 900标准的Von Mises准则法兰完整性评估。开发了一个三级自动预警系统（GREEN、AMBER、RED），持续将计算得到的Von Mises利用率映射到操作风险水平，并带有用于早期预警的升级率图。该模型针对Shell Nigeria Bonny Light现场数据进行了验证，显示出强大的预测精度（R² = 0.905，r = 0.988）。模拟结果表明，WAT（318 K）在11.43公里处达到，使得71.5%的管道容易发生沉积。30天后，峰值沉积在距离入口25至27公里之间达到10.99 mm（20.5 kg/m³）。这种积累导致总压降增加了6.6%（从0.9745 MPa到1.0382 MPa）。Von Mises法兰应力在入口处达到峰值66 MPa，对应的利用率为0.40，确认在整个30天模拟期间所有空间节点仍处于GREEN安全分类中。为了维持管道健康，清管间隔应从最初的28天缩短至14–21天，机械清管集中在20–32公里段。未来的工作应包括瞬态流动效应、利用Bonny Light样品进行实验室流动回路验证以完善文献中得出的蜡动力学常数，以及多段几何形状。

在石油和天然气行业中，海底管道输送原油时，蜡沉积是一项关键的流动保障挑战。当流体温度低于蜡出现温度（WAT）时，原油中的石蜡结晶并沉积在管壁，导致有效管径减小、泵送成本增加，严重时可能引发压力尖峰甚至完全堵塞。现有研究多集中于蜡沉积的预测模型，但忽略了蜡沉积引起的压降变化对法兰连接件结构完整性的影响。法兰作为管道与仪表之间的关键机械连接件，承受轴向和弯曲载荷，当蜡沉积导致有效内径减小、压力升高时，可能超出原始设计规格。因此，有必要开展一体化研究，同时评估蜡沉积对流动保障和结构完整性的影响，以优化管道维护策略。

研究人员在《Chemical Thermodynamics and Thermal Analysis》上发表了本研究，开发了一个集成数学模型，用于模拟海底原油管道中的蜡沉积并自动评估法兰完整性。该模型基于Bonny Light原油的现场数据，耦合了热力学溶解度理论、沉淀动力学、成核动力学、Fick传质，以及改进的Darcy-Weisbach压降分析和符合ASME B16.5 Class 900标准的Von Mises准则法兰应力评估。模型在MATLAB中实现，采用有限差分空间离散和显式欧拉时间步进，模拟了40公里管道在30天内的蜡沉积过程。通过Shell Nigeria Bonny Light现场数据验证，模型预测精度高（R²=0.905，r=0.988）。研究结论表明：WAT（318 K）在11.43公里处达到，71.5%的管道长度处于沉积风险区；30天后蜡沉积峰值厚度为10.99 mm，浓度20.5 kg/m³，位于入口25–27公里处；总压降从0.9745 MPa增至1.0382 MPa，增幅6.6%；Von Mises法兰应力在入口处最大为65.80 MPa，利用率0.39，所有空间节点均处于GREEN安全分类；建议将清管间隔从28天缩短至14–21天，并集中清管20–32公里段。该研究为优化清管计划和提高管道可靠性提供了定量依据，填补了尼日利亚原油管道数据空白，并与尼日利亚石油工业法案的目标一致。

为开展研究，研究人员采用了以下主要关键技术方法：①基于Van't Hoff方程的热力学溶解度模型，忽略压力效应，仅考虑温度对溶解度的影响；②沉淀动力学模型，以过冷度（ΔT）、壁面剪切速率（γ˙_wall）和蜡含量（wt%）为驱动因素，引入对数时间硬化项以模拟沉积老化；③经典成核理论，结合过饱和度（σ）和剪切修正，计算临界成核半径（R_critical），并考虑剪切对成核的破坏作用；④Fick扩散定律计算质量传递通量，考虑浓度梯度；⑤Darcy-Weisbach方程结合有效直径计算局部压降，使用Haaland公式计算摩擦因子；⑥Von Mises准则计算法兰等效应力，包括轴向、环向和径向应力分量，并与ASME B16.5 Class 900标准中允许应力比较得到利用率；⑦三级自动预警系统（GREEN/AMBER/RED），基于利用率阈值（0.60和0.85）进行风险分类，并定义升级率以提供早期预警。样本来源为Bonny Light原油，现场数据来自Shell Nigeria。

研究结果涵盖以下内容：
- **3.1.1 温度分布分析**：通过模拟得到管道内温度沿程指数衰减，热衰减常数α=1.47×10^-4 m^-1，WAT（318 K）在11.43公里处达到，剩余71.5%管道处于低于WAT的沉积风险区。
- **3.1.2 蜡沉积特性**：30天后蜡沉积呈非均匀分布，峰值沉积浓度20.5 kg/m³出现在距入口约27公里处，沉积速率受过冷度、剪切诱导成核和分子扩散共同控制。
- **3.1.3 蜡厚度分析**：蜡沉积厚度从WAT起始点逐渐增加，在27公里处达到最大约10.99 mm，之后略有下降，反映了热驱动力和流体动力学条件的竞争。
- **3.1.4 蜡厚度积累分析**：时间上，蜡厚度呈对数增长，前10天快速积累，之后趋近准稳态；空间上，20天内积累显著，之后生长速率因热平衡和剪切去除增强而降低。
- **3.1.5 蜡浓度梯度分析**：浓度梯度沿管道单调增大，从WAT起始点的约-2000 kg/m⁴增加到出口的约-9000 kg/m⁴，表明扩散驱动力增强。
- **3.1.6 与现场数据对比**：模型预测与Shell Nigeria Bonny Light现场厚度数据高度相关（Pearson r=0.988，R²=0.905），误差分布呈正弦模式，平均绝对误差低于8%，所有数据点位于±10%置信区间内。
- **3.1.7 沿管道压降梯度**：在WAT之前，压降梯度约0.0235 kPa/m；30天后，在25公里处峰值梯度约0.0275 kPa/m，较清洁管道基线增加17%。
- **3.1.8 总压降分析**：30天内总压降从0.9745 MPa增至1.0382 MPa，增幅6.6%，中间值呈对数增长，与沉积动力学一致。
- **3.1.9 内部压力分布**：压力沿程下降，30天后在27公里处压力从14.37 MPa降至14.28 MPa，出口处从14.08 MPa降至13.82 MPa，仍高于最低输送压力。
- **3.1.10 法兰应力分析（解析法与Von Mises准则）**：解析法总应力范围22–23 MPa，利用率0.13–0.14；Von Mises应力在入口处最大65.80 MPa（利用率0.39），沿程递减至57.02 MPa（利用率0.34）。Von Mises准则因包含环向和径向分量而更为保守，推荐作为主要评估准则。
- **3.1.11 自动预警分类**：所有1000个空间节点在30天内均处于GREEN分类（利用率0.34–0.39），升级率均低于阈值0.01/天，未触发预警。表明蜡沉积的主要后果是流动保障退化而非结构风险。
- **3.2 敏感性分析**：对峰值厚度和总压降的敏感性分析表明，热传递系数（U¨）和蜡含量（wt%）是最高敏感参数（±13–19%变化每±20%变化），剪切去除常数（Cs）中等敏感，熔化焓（ΔH_fusion）敏感度低，确认压降预测的稳健性。

讨论部分总结了研究结论，研究人员翻译如下：
本研究开发并验证了一个集成数学框架，用于海底原油管道中的蜡沉积预测和法兰完整性自动评估，并应用于尼日利亚Bonny Light原油系统。得出以下结论：
1. 开发了耦合热力学蜡溶解度（Van't Hoff）、沉淀动力学、经典成核理论、Fick传质和Darcy-Weisbach压降公式的综合数学模型，模型通过Shell Nigeria Bonny Light现场数据验证，R²=0.905，r=0.988，为尼日利亚原油管道基础设施提供了本地校准工具。
2. 模型预测蜡出现温度（318 K）在入口11.43公里处达到，71.5%的40公里管道易受沉积影响。峰值沉积厚度10.99 mm和峰值蜡浓度20.5 kg/m³出现在25–27公里处，与Bonny Light原油沉积模式的现场观测一致。
3. 30天内蜡积累导致总压降增加6.6%（0.9745 MPa至1.0382 MPa），最陡压降梯度位于约25公里处，局部阻力较清洁管道基线高17%。
4. 集成法兰完整性框架将蜡引起的压力变化转化为符合ASME B16.5 Class 900标准的Von Mises等效应力估计。Von Mises应力在入口处峰值65.80 MPa（利用率0.39），因捕获完整三轴应力状态而推荐作为主要评估准则。
5. 三级自动预警系统（GREEN/AMBER/RED）阈值校准至ASME B16.5、API RP 579-1和DNV-ST-F101，所有1000个空间节点在30天模拟中保持GREEN分类，利用率在0.34–0.39之间，未触发升级标志。这确认在模拟操作条件下，蜡沉积的主要后果是流动保障退化而非即时结构风险。
6. 敏感性分析识别出整体热传递系数U¨和蜡含量wt%为最高敏感参数（±20%变化引起峰值厚度±13–19%变化），强调了现场准确表征这些参数对于特定地点部署的重要性。
7. 操作上，模型支持缩短清管间隔从28天至14–21天，集中清管活动于20–32公里段，并将生产速度提高至2.0–2.5 m/s以增强剪切去除。应在蜡形成区域每5公里安装温度传感器。
8. 未来工作应包括：使用Bonny Light原油样品进行实验室流动回路实验以从第一性原理优化动力学常数Cs和k；考虑瞬态流动效应（如段塞流和启停热循环）；多段几何形状（变壁厚和绝缘）；以及与机器学习模型集成实现实时参数更新。