流动保障、多相流模拟、水合物与蜡沉积防治、管道设计与操作优化的知识全景图
海底管道连接井口与处理设施，跨越几公里到几百公里。流动保障（Flow Assurance）的核心任务：确保油气安全、连续地从海底输送到平台或岸上。
| 风险 | 成因 | 后果 |
|---|---|---|
| 水合物堵塞 | 低温高压下天然气+水形成冰状结晶 | 管道堵塞，停输，清管代价极高 |
| 蜡沉积 | 温度低于析蜡点，石蜡析出附着管壁 | 流通截面减小，压降增大，堵塞 |
| 段塞流 | 多相流态下液相间歇性聚集成塞 | 分离器过载、压力波动、设备损坏 |
| 腐蚀 | CO₂/H₂S 溶于水形成酸性环境 | 管壁减薄、穿孔、泄漏 |
| 砂沉积 | 地层出砂随流体进入管道 | 磨蚀、堵塞、设备损坏 |
| 沥青质沉积 | 压力温度变化导致沥青质析出 | 地层/管道堵塞，影响生产 |
油气水三相在管道中的分布形态决定了压降、持液率和段塞特性。流态随流量、管径、倾角和流体物性变化。
| 流态 | 特征 | 风险等级 | 典型场景 |
|---|---|---|---|
| 泡状流（Bubble） | 气相以小气泡分散在液相中 | 低 | 高液量、低气量 |
| 段塞流（Slug） | 液塞与气塞交替，液塞高速通过 | 高 | 中低流量、地形起伏 |
| 分层流（Stratified） | 气液分层，液相在底部流动 | 低 | 低流量、水平管道 |
| 环状流（Annular） | 液相沿管壁成环，气芯在中央 | 中 | 高气量 |
| 分散泡状流 | 气泡极细密地分散在液相中 | 低 | 极高液量 |
以表观气速（Vsg）和表观液速（Vsl）为坐标轴，绘制流态分区图。常用的有 Taitel-Dukler 图、Baker 图和 Mandhane 图。
实际工程中，流态判定由模拟软件自动完成（基于 Beggs-Brill、OLGA 等模型）。
摩擦压降：流体与管壁摩擦引起的压力损失，取决于流速、管径、粗糙度和流态
重力压降：管道高程变化引起的静水压力差，海底管道地形起伏大时尤其重要
加速压降：流体密度变化引起的动量变化，多相流中通常较小
| 参数 | 含义 | 单位 | 软件视角 |
|---|---|---|---|
| 持液率（Hold-up） | 管道截面上液相占据的面积分数 | 无因次（0~1） | 时序数据，每个网格点每个时间步 |
| 表观速度 | 假设该相单独占据全部截面时的流速 | m/s | 输入条件或计算结果 |
| 混合密度 | 气液混合物的等效密度 | kg/m³ | 重力压降计算的基础 |
| 混合黏度 | 气液混合物的等效黏度 | Pa·s | 摩擦压降和雷诺数计算 |
| 段塞长度 | 单个液塞的长度 | m | 分离器容积设计依据 |
| 段塞频率 | 单位时间内通过某点的液塞数 | 1/h | 设备疲劳和载荷分析 |
多相流模拟的本质是数值求解流体力学方程组。算法按重要性分为四层：地基层（压降+传热）→ 瞬态层（双流体模型）→ 风险层（水合物/蜡）→ 数值层（收敛求解）。如果地基层算不对，上面全部白费。
第一层（地基）：给定流量和管道参数，计算沿线压力和温度分布。常用模型有 Beggs-Brill（倾斜管）、Dukler（水平管）。稳态求解需要 P-V-T 三重闭环迭代。
第二层（传热）：多层热阻模型计算管道热损失。保温层是核心控制项，设计得好能减少 80%+ 热损失。冷却时间是停输应急的关键指标。
第三层（瞬态）：双流体模型（OLGA 核心）分别求解气相和液相的质量/动量守恒方程。段塞追踪在此基础上追踪每个液塞的前沿和后沿位置。
第四层（风险）：van der Waals-Platteeuw 模型预测水合物形成条件，Hammerschmidt 公式计算抑制剂用量。蜡沉积速率由分子扩散和剪切剥离的平衡决定。
给定恒定的输入条件，计算管道在稳定运行状态下的全线压力和温度分布。
管径选择：给定输量，多大管径能让终点压力满足处理设施入口要求？
保温设计：多厚保温层能让终点温度高于水合物形成温度？
最大输量：给定管道规格和入口压力，最大能输多少？
注醇量：需要注入多少甲醇/MEG 来抑制水合物？
模拟管道在启动、停输、流量变化、清管等动态过程中的行为。这是流动保障最核心的模拟类型。
| 瞬态工况 | 场景描述 | 核心关注 |
|---|---|---|
| 启动（Startup） | 管道从停输状态恢复生产 | 冷却后的管道内流体状态、注醇时机、升温速率 |
| 停输（Shutdown） | 计划停输或紧急关井 | 冷却时间（Cool-down Time）、水合物形成风险窗口 |
| 冷却（Cool-down） | 停输后管道温度随时间下降 | 温度降到水合物形成温度所需时间 = 应急响应窗口 |
| 清管（Pigging） | 发送清管器清除积液/蜡 | 清管前后的段塞量、接收端容量 |
| 流量变化 | 产量递减或增产调整 | 流态转换、段塞风险变化 |
| 泄漏模拟 | 管道破裂或穿孔 | 泄漏检测、压力传播、泄漏量估算 |
停输后管道温度持续下降。当温度降到水合物形成温度以下时，管道面临堵塞风险。
冷却时间 = 停输时刻温度 降到 水合物形成温度所需时间。这个时间窗口决定了：
- 操作人员有多长时间进行应急处理（注醇、排空等）
- 保温层厚度设计的依据（保温越好，冷却越慢）
- 是否需要设置死管（Deadleg）保温或电伴热
天然气水合物在特定温度-压力范围内形成，可用相平衡图判定：
典型条件：0~25°C，压力 > 30 bar。
预测方法：Hammerschmidt 公式、Makogon 方法、PVT 软件（PVTSim、Multiflash）。
抑制剂：甲醇（MeOH）、乙二醇（MEG）、动力学抑制剂（KHI）、防聚剂（AA）。
| 防治策略 | 原理 | 适用场景 | 软件模拟要点 |
|---|---|---|---|
| 热力学抑制（注醇） | 降低水的活度，移动水合物相平衡线 | 短距离管道、停输应急 | 计算所需注醇量 vs 水合物裕度 |
| 保温 | 减少热损失，维持管内温度高于水合物温度 | 长距离管道 | 保温层厚度 vs 终点温度的灵敏度分析 |
| 电伴热 | 主动加热管道，维持温度 | 关键管段（立管、死管） | 功率需求 vs 冷却时间延长 |
| 降压 | 降低系统压力，提高水合物形成温度门槛 | 停输应急 | 降压操作的安全窗口 |
| 动力学抑制（KHI） | 延缓水合物晶体生长速率 | 深水长距离 | 过冷度（Subcooling）在 KHI 有效范围内 |
原油温度低于析蜡温度（WAT）时，石蜡分子结晶析出，在管壁上沉积。沉积速率取决于：
- 管壁与流体的温差（分子扩散驱动）
- 剪应力（高流速可冲刷部分蜡层）
- 原油组成（高含蜡原油风险大）
WAT（Wax Appearance Temperature）：析蜡温度，蜡开始析出的温度
Pour Point：倾点，原油失去流动性的温度
Cloud Point：浊点，原油出现浑浊的温度（通常 ≥ WAT）
| 蜡防治策略 | 原理 | 软件模拟 |
|---|---|---|
| 保温 | 维持管内温度高于 WAT | 沿线温度曲线 vs WAT 线 |
| 清管 | 定期发送清管器刮除蜡层 | 清管频率优化：沉积量 vs 清管风险 |
| 化学防蜡剂 | 抑制蜡结晶或改变蜡的结构 | 药剂注入量模拟 |
| 热洗 | 注入热流体融化蜡沉积 | 热洗温度和持续时间 |
水动力段塞（Hydrodynamic Slug）：管道内气液相间的自然流动不稳定性产生的液塞。与流量、管径、倾角相关。
地形段塞（Terrain Slug）：管道沿低点积水，积到一定程度被气流推动形成大段塞。在起伏地形尤其明显。
操作段塞（Operational Slug）：启动、清管、流量变化等操作引起的瞬态段塞。通常液塞最大。
| 段塞管理策略 | 说明 | 软件模拟输入/输出 |
|---|---|---|
| 段塞捕集器 | 在管道出口设置缓冲容器，吸收液塞 | 输入：段塞体积和频率 → 设计段塞捕集器容积 |
| 流量控制 | 调整气液比，避开段塞流区域 | 输入：流量范围 → 流态图上的操作区域 |
| 顶部节流 | 在出口设置节流阀，背压抑制段塞 | 输入：节流开度 → 段塞频率/幅度变化 |
| 主动控制 | 基于实时数据反馈控制阀门 | 输入：实时 P/T → 控制策略优化 |
多相流模拟的准确性高度依赖流体物性（PVT）模型的精度。流体物性描述了油、气、水各相在不同温度压力下的物理性质。
| 参数 | 含义 | 对模拟的影响 |
|---|---|---|
| 泡点压力（Pb） | 液相开始析出气泡的压力 | 决定管道中气相出现的条件 |
| 气油比（GOR） | 单位体积原油中溶解气的体积 | 决定管内气液比和流态 |
| 原油体积系数（Bo） | 地层油体积与地面油体积之比 | 储量计算和流量转换 |
| 黏度曲线 | 各相黏度随 T/P 变化的关系 | 直接影响压降计算 |
| 密度曲线 | 各相密度随 T/P 变化的关系 | 影响重力压降和持液率 |
| 溶解度 | 气相在液相中的溶解量随 P/T 变化 | 组分模型的核心输入 |
| 工具 | 类型 | 核心能力 | 厂商 | 典型场景 |
|---|---|---|---|---|
| OLGA | 瞬态多相流 | 1D 瞬态多相流模拟，段塞、清管、启动停输 | Schlumberger | 操作窗口定义、瞬态分析 |
| PIPESIM | 稳态管网 | 稳态多相流管网模拟，管径优化 | Schlumberger | 概念设计、管径选择 |
| LedaFlow | 瞬态多相流 | OLGA 的主要替代品，4x4 模型更精细 | Kongsberg | 瞬态分析、段塞预测 |
| PVTSim / Multiflash | PVT 物性 | 流体相行为、水合物预测、注醇量计算 | Calsep / KBC | PVT 建模、水合物分析 |
| HYSYS | 流程模拟 | 全流程工艺模拟，设备设计 | AspenTech | 处理设施设计、能量平衡 |
| PROSPER | 井筒-管道耦合 | 单井与管道耦合的节点分析 | Petroleum Experts | 生产系统优化 |
| GAP | 管网优化 | 多井多管道网络优化 | Petroleum Experts | 管网调度优化 |
概念设计阶段：PIPESIM（稳态快速筛选方案）→ 确定管径和保温方案
详细设计阶段：OLGA/LedaFlow（瞬态验证设计方案的动态行为）
操作阶段：OLGA 实时模型（在线监测和预测）+ 在线 PVT 更新
PVT 建模：PVTSim/Multiflash 建立流体模型 → 导入 OLGA/PIPESIM
从软件开发角度，海底管道混输模拟系统的核心数据实体和关系。
| 实体 | 核心字段 | 说明 |
|---|---|---|
| 管道（Pipeline） |
管道编号、名称、类型（出口/进口/回注） 总长度、管段数 起终点坐标、水深 |
管道的顶层实体，由多个管段组成 |
| 管段（Segment） |
内径、外径、壁厚 保温层类型、厚度、导热系数 粗糙度、高程起点/终点 海水温度、埋深、环境传热系数 |
管道按特征（管径、保温、地形）划分为多个管段 |
| 高程剖面（Profile） |
距离（m）、高程（m） 离散点数（通常 50~500 个点） |
管道沿线的地形起伏，是重力压降和段塞分析的关键输入 |
| 流体模型（Fluid） |
模型类型（黑油/组分） 黑油参数：API、GOR、含水率 组分参数：C1~C30+ 摩尔分数 PVT 表或 EOS 参数 |
定义流体物性。一个流体模型可被多个工况引用 |
| 边界条件（Boundary） |
入口：流量（油/气/水）、温度、压力 出口：压力设定 或 流量设定 环境条件：海水温度剖面 |
模拟计算的输入条件 |
| 模拟工况（Scenario） |
工况名称、类型（稳态/瞬态） 关联：管道 + 流体 + 边界条件 瞬态控制：时间步长、总模拟时间 操作事件序列（启动/停输/清管/注醇） |
一次完整的模拟任务定义 |
| 模拟结果（Result） |
时序数据：时间 × 管段网格点 × 参数 参数：压力、温度、持液率、流速、流态 段塞统计：液塞体积、频率 水合物裕度：沿线温度 vs 水合物温度 |
三维数据矩阵（时间-空间-参数），数据量大 |
| 风险报告（Risk） |
水合物风险区域（管段 + 时间段） 蜡沉积厚度预测 段塞频率和体积 冷却时间 |
基于模拟结果的风险判定 |
稳态模拟：一条管道一次计算，输出 ~100~500 个空间点 × ~10 个参数 = 数千个数据点
瞬态模拟：时间步长 1~60 秒，模拟 24~72 小时 = 数千个时间步 × 100~500 空间点 × 10 参数 = 百万~千万级数据点
多工况对比：一个项目通常包含 10~50 个工况，总数据量可达 GB 级
- 结果数据用时序数据库或 HDF5 存储，不用关系数据库存原始网格数据
- 提供聚合接口：最大值、最小值、时间平均、空间平均
- 支持结果对比：同管道不同工况的 P/T 曲线叠加
| 缩写/术语 | 全称 | 含义 |
|---|---|---|
| Flow Assurance | 流动保障 | 确保油气从井口到处理设施安全、连续输送的技术体系 |
| Multiphase Flow | 多相流 | 油、气、水（有时含砂）在同一管道中混合流动 |
| Hold-up | 持液率 | 管道截面上液相占据的面积分数 |
| Slug | 段塞 | 液相间歇性聚集形成的液塞 |
| Cool-down Time | 冷却时间 | 停输后温度降到水合物形成温度的时间窗口 |
| WAT | Wax Appearance Temperature | 析蜡温度，蜡开始结晶析出的温度 |
| GOR | Gas-Oil Ratio | 气油比，单位体积原油中溶解气的体积 |
| MEG | Monoethylene Glycol | 乙二醇，常用的水合物热力学抑制剂 |
| KHI | Kinetic Hydrate Inhibitor | 动力学水合物抑制剂，延缓水合物晶体生长 |
| FEED | Front End Engineering Design | 前端工程设计，详细设计之前的关键阶段 |
| EOS | Equation of State | 状态方程，描述流体相行为的数学模型（PR、SRK） |
| PVT | Pressure-Volume-Temperature | 描述流体物性随压力、温度变化的关系 |
| SCADA | Supervisory Control and Data Acquisition | 数据采集与监控系统，管道运营的实时数据源 |
| DCS | Distributed Control System | 分布式控制系统，平台/炼厂的控制系统 |
| Pig / Pigging | 清管器 / 清管 | 发送管道内检测器或清管器清除杂物/积液 |
| FID | Final Investment Decision | 最终投资决策 |
| Subcooling | 过冷度 | 实际温度低于水合物形成温度的差值 |
数字孪生: 在线实时模型从"锦上添花"变为标准配置，SCADA 驱动的实时模拟成为管道运营的核心工具。
云计算模拟: 多工况并行模拟从本地工作站迁移到云端，大幅缩短设计周期。
数据标准化: WITSML/PRODML 等标准在油气行业加速普及，模拟数据与生产数据的互操作性提升。
机器学习辅助: 用 ML 模型近似 OLGA 等物理模型，实现毫秒级响应的实时预警。
自动化工况分析: 自动生成数百个工况的参数扫描，替代人工逐个调试。
与结构/腐蚀耦合: 流动保障-结构-腐蚀一体化分析，减少各专业之间的数据断档。
模型迷信: 再精确的模型也依赖输入数据质量。错误的 PVT 数据会让最精密的模拟产生谬误。
工具替代理解: 会用 OLGA 不等于理解多相流。工具使用必须建立在理论基础之上。