论文题目：基于FCE-SPH的海上溢油扩散模拟方法

录用期刊：计算机辅助设计与图形学学报 (CCF A)

录用时间：2024.2.21

作者列表：

1） 连远锋 中国石油大学（北京）信息科学与工程学院/人工智能学院 计算机系教师

2） 高浛钊 中国石油大学（北京）信息科学与工程学院/人工智能学院 计算机技术专业 硕21

3） 魏昊天 中国石油大学（北京）安全与海洋工程学院 博士研究生

3） 纪连恩 中国石油大学（北京）信息科学与工程学院/人工智能学院 计算机系教师

4） 董绍华 中国石油大学（北京）安全与海洋工程学院 安全工程系教师

摘要:

针对现有的海面原油泄漏模拟方法无法呈现符合物理过程的溢油风化扩散效果及模拟速度过慢问题，首先给出一种不混溶Cahn-Hilliard两相流乳化光滑粒子流体动力学扩散模拟方法FCE-SPH，结合脉动风场描述海上溢油的混合和乳化分解过程；然后设计一个两相流物理感知深度强化学习偏置图注意力网络OTF-DRL，以加速流体模拟仿真，通过异构偏置图注意力模块，增强网络获取粒子特征的表达能力。在不同的工况条件下，分别对原油在海上泄漏扩散物理过程和油品风化过程进行分析，实验结果表明，所提方法可以真实反映海上溢油扩散现象，OTF-DRL模拟计算时间缩短49.63%。

溢油两相流FCE-SPH建模

基于海面溢油扩散物理过程，本文提出结合随机脉动风场的Cahn-Hilliard两相流乳化光滑粒子动力学模型FCE-SPH。首先在SPH框架中引入Cahn-Hilliard相分离模型，然后在相分离模型基础上构建油水乳状液形成模型，最后结合随机脉动风场方法和表面粒子分类器实现海上风场模拟。

1. 基于Cahn-Hilliard方法的相分离机制

Cahn等提出了基于统计热力学描述相分离过程的Cahn-Hilliard方程，已在多相流研究中得到了广泛应用。本文将Cahn-Hillard方程与海面原油泄漏扩散模拟结合，得到离散化可以得到粒子的相分离方程为：

其中，为粒子的浓度且值域为[−1，1]；为粒子质量；是粒子的位置；为支撑半径为的核函数；为浓度；μ为化学势；为流动系数；ρ为粒子密度。

图1 Cahn-Hilliard油水两相流相分离模型

2. 乳化物生成模型

海面溢油污染通常包括扩散过程和风化过程，风化过程中乳化阶段产生的油水乳化物的流变性会发生变化，随着时间的推移直接影响其扩散过程。如图2所示，将Rasmussen模型中乳化反应视为平衡反应，即

图2 油水乳状液形成模型

设是乳化物粒子最大含水量，当时，正向反应的概率为

其中，为与温度相关的正向反应系数；v₀为黏度系数；是非压力项合力(如风场拖曳力等)。当时，乳化物分解为油和水的概率为

其中，是逆向反应系数，与原油中沥青等活性剂含量相关；是沥青百分比。

3. 随机脉动风场模型

结合拖曳力模拟空气对基于粒子的流体影响,在气液耦合方面产生了很好的效果。本文利用邻域粒子质量分类器判断表面粒子集，并结合随机脉动风场模型对粒子施加拖曳力以模拟风场作用。设表面粒子分类器的邻域粒子质量特征为，邻域粒子数量特征为，则有

其中，是粒子的位置；是邻域粒子数量；是分类器输出结果。基于Perlin噪声，使用时间作为随机数种子产生伪随机序列

其中，表示粒子受到的风场拖曳力；为噪声振幅。

OTF-DRL网络

如图3所示，FCE-SPH算法中多相流约束模型增加了算法的时间复杂度，为此本文提出了OTF-DRL网络，用于加速模拟不混溶两相流.

图3 OTF-DRL模型结构图

在两相流模拟过程中，粒子浓度与邻域粒子频繁交换变化，为了描述不同浓度粒子扩散规律，基于节点浓度的差异构建异构图，设计了基于偏置注意力的异构图注意力网络，通过计算粒子与邻域粒子之间的注意力系数，达到特征提取目的。

实验结果及分析:

图4 6种算法模拟结果对比

图4给出了6种算法的海上溢油模拟结果，假设海上风场方向为逆溢油方向，风速设置为15 m/s。可以看出，PBF算法存在粒子发散现象，第150帧至第350帧油粒子分散较为明显，泄漏体边缘不够光滑；DFSPH算法虽没有出现粒子发散现象，但原油溢出后第250帧产生粒子堆积现象，溢油边缘没有体现油水两相流耦合效果；Multiphase SPH较好地处理了两相流的交互运动效果，泄漏体边缘呈流线型，并且由于溢油喷射具有初速度，泄漏体两侧呈漩涡状分布；ICSPH泄漏体边缘不光滑且有分散的油粒子，溢油内部出现不稳定现象并产生空腔，如表1所示，由于FCE-SPH算法融入相分离、乳化和随机脉动风场模型，因此不同程度地增加了原油泄漏扩散模拟时间。

表1 6种算法海上溢油模拟用时对比

连远锋，教授，硕士生导师。研究方向为图像处理与虚拟现实、机器视觉与机器人、深度学习与数字孪生。