基于模拟仿真的石化油气输送系统数字孪生体技术
来源:用户上传
作者:周波
摘要：数字孪生是真实物体在虚拟世界的一种建构过程，形成物理维度上的实体和信息的虚拟世界同生共存、虚实交融的数字发展新格局。研究通过对天然气管道电子仪表系统及数字孪生体的搭建，运用数据融合与数字模拟仿真大数据技术体系，对数字孪生体的仿真度M行验证分析，验证结果证明了数字孪生体算法仿真能够有效的对异常数据预警和主动推送，能够提高系统应用的安全性。
关键词：数字孪生体;天然气;模拟仿真;输送管道
中图分类号：TE83
文献标志码：A文章编号：1001-5922（2022）09-0102-04
Applicationtechnologyofdigitaltwinofnaturalgastransmissionpipelinebasedonsimulation
ZHOUBo
（CPECCSouthwestCompany，Chengdu610041，China）
Abstract：Digitaltwinisaconstructionprocessofrealobjectsinthevirtualworld，formthephysicaldimensionofentityandinformationvirtualworldcoexistence，integrationofdigitaldevelopmentpattern，throughthegaspipelineelectronicinstrumentsystemanddigitaltwin，usingdatafusionanddigitalsimulationsystem，verificationanalysisofdigitaltwinsimulation，toprovethedigitaltwinalgorithmsimulationcaneffectiveabnormaldatawarningandactivepush，canimprovethesecurityofthesystemapplication.
Keywords：digitaltwin;naturalgas;simulation;transmissionpipeline
天然气能否准确计量直接影响着开发、运行、配给及用户等各方的经济利益，而管路结构的安全性是影响天然气流量计量准确性的主要因素之一。天然气管道事故常处于突发情况，需要建立能即使发现天然气输送管道的应急响应技术[1]。
数字孪生是通过物理模型、传感器更新、运行历史等数据，包括多学科、多物理量和时间尺度及多概率的实时仿真过程，在虚拟空间中进行映射，创造一个能表达其相对应的实体设备全部周期过程的数字模型[2]。数字孪生体的工程意义是在实际运行系统之外，构建一个与实际运行系统运行环境基本一致的虚拟化系统，且使虚拟化系统的运行时间前推一定周期，从而判断实际运行系统在未来短周期内可能发生的故障，并作出提前预警[3]。
受制于工程逻辑学的局限性，实际搭建液化天然气（（LiquefiedNaturalGas，简称LNG）输送管道数字孪生体时，存在海量的隐藏传导机制和不完备数据，所以，只能通过机器学习系统搭建仿真系统，而无法实现完全一致且无需后续干预的数字孪生体[4]。建筑信息模型（BuildingInformationModeling，简称BIM）系统中内置的数字孪生体的时间前推能力难以满足实际运行需求。[5]
该研究设计了一套相对完善的基于机器学习算法的数字孪生体仿真运行算法，并对其进行了验证。[6]
1天然气管道电子仪表系统及数字孪生体的搭建
今天，大部分LNG天然气管道已经实现了压力、流量、温度等仪表的1s周期远程实时抄表工作且形成了物联网大数据体系。随着量子力学发展，使用高精度三轴激光陀螺仪捕捉管道的震动，可以有效反映出管道内的气液混合情况、涡流湍流情况等[7]。数字孪生体建设中，这些仪表采集的物联网大数据称为其主要数据来源。[8]
该研究设计的数字孪生体，对各探头数据在BIM系统中构建数据LOD（LevelsofDetail，简称LOD）技术图，且在授时系统驱动下形成图时序序列，进而使用模糊神经网络对所有图序列分别执行卷积压缩，使每个图序列形成1个双精度浮点型特征变量，再将这些特征变量输入到多列神经网络中，重新组建各个压力表、流量表、温度表在前推时间点的曲线估计前推结果数据[23]。该数字孪生体构建策略如图1所示：
图1中，使用了图论统计法、空间卷积法、时序图序列同构化方法等数据前置治理算法，模糊神经网络、对数神经网络等机器学习算法，多列神经网络数据重构法等数据解模糊算法等，形成数字孪生体中前推一定时间周期的数据曲线估计前推模型。相关算法将在后续章节中展开论述[911]。
2数据融合与数字模拟仿真大数据技术体系
2.1数据的前置治理算法
2.1.1图论统计学方法
使用LOD法将仪表数据定义到BIM的三维模型中，形成数据投影。数据表达为{A，x，y，z，b1，b2，n，t}，其中，{A}为仪表序列号，{x，y，z}为仪表测量点的BIM三维空间坐标，{b1，b2}为仪表管道逻辑关系前后连接点，{n}为仪表读数;{t}为仪表读数时间戳;将相同{A，x，y，z，b1，b2}标识的序列在时间轴{t}上展开，形成图序列原始统计结果;
2.1.2四维时空空间卷积方法
四维时空空间卷积方法是将上述{A，x，y，z，b1，b2，n，t}序列，以{x，y，z，t}四个坐标维度作为卷积维度，对{n}值执行空间卷积，强化特征数据，其基函数如公式（1）：
nlc202210121127
Gx，y，z，t=∫+
SymboleB@
-
SymboleB@
gbxm-bym-bzm-btm-bdb（1）
式中：Gx，y，z，t为卷积后输出Y果;gb为卷积核，为一个小规模四维阵列数据，不同的卷积核赋值方式会实现不同的卷积效果，如边缘强化、离群捕捉、峰值强化等;xm-bym-bzm-btm-b为被卷积序列分量函数;
2.2基于机器学习的数据深度挖掘算法
2.2.1基于六阶多项式深度迭代回归算法的模糊神经网络
六阶多项式深度迭代回归算法的统计学本质是将大宗数据卷积到1个双精度浮点型变量内，使其在最小的信息损失量前提下使用较小规模数据代表大宗数据的数据特征。六阶多项式深度迭代回归算法的节点基函数如公式（2）：
y=∑ni=1∑5j=0Ajxji（2）
式中：j为多项式阶数;i为遍历指针;n为前一层神经网络的节点总量;xji为上一层神经网络中第i个节点输入值的j阶次方;y为该神经网络节点的输出值;Aj为第j阶多项式的回归系数，即该公式内每个节点包含A0至A5共6个待回归系数。
2.2.2基于对数深度迭代回归算法的对数神经网络
对数深度迭代回归算法的统计学本质是利用对数曲线的非线性特征放大数据在自变量坐标轴0值附近的细节，该数字孪生体算法中，将其作为多列神经网络的节点函数，其基函数如公式（3）：
y=∑ni=1A・logxi+B（3）
式中：A、B为待回归系数;其他数学符号含义同式（2）。
3数字孪生体的仿真度验证分析
3.1仿真验证方案
在Matlab中加载Simulink控件构建仿真环境，数据训练中，采用某单位LNG天然气输送管道中压力、流量、温度、震动加速度仪表数据作为原始数据作为原始数据，使用中间数据作为输入数据，对应的前推数据作为输出数据执行数据训练，后续类似数据组织方式组织收敛效果验证。每个特定分列的算法模块内部逻辑结构相同，根据不同的大数据训练需求对其分别训练，使其向特定的应用方向数据收敛。
验证试验中，使用仿真数据与原始数据之间的决定系数、标准偏差率和平均离群率作为验证方法。其中：
决定系数（R2）的算法如式（4）：
R2=∑i（xi-x）∑i（xi-xi），x=1n∑ni=1xi（4）
式中：x为考察样本序列算数平均值;xi为序列中第i个回归值;xi为序列中第i个输入值;n为考察样本数量。
标准偏差率的算法如式（5）：
σ=1n-1∑ni=1xi-μ2，μ=1n∑ni=1xi（5）
式中：σ为输入序列x的标准偏差率;xi为输入序列x中的第i个输入项;μ为输入序列x的算数平均值;n为统计样本数量。
平均离群率的算法如式（6）：
B=1n∑ni=1axi+byi+cA2+B2・1max（B）-min（B）
（6）
式中：B为标准离群率;a、b、c为离群率线性回归预测直线的系数;（xi，yi）为第i个散点的二维空间坐标;max（B）;min（B）为离群率估计的最大值和最小值。
对数据孪生体前推时间周期曲线估计预测数据的算法效能验证方法使用验证敏感度的方式进行计算，敏感度计算方法为为数据预测值与数据实测值之差与数据实测值的比值均值，其算法如式（7）：
M=1n∑ni=1ri-firi（7）
式中：n为统计数据节点的数量;ri为第i个数据实测值;fi为第i个数据预测值。
3.2流量分布数据耦合度
对天然气流量分布数据流在定时域窗口进行值域区间的数据频率分布统计和频域区间的数据分散度运算，可以建立动态的流量分布耦合数据，比较分析流量分布数据耦合度的统计如表1所示。
表1中，决定系数一般指可决系数，表示一个随机变量与多个随机变量关系的数字特征;标准偏差率时一种度量数据分布的分散程度标准，用以衡量实际值与理论值或者估计值的偏差程度;平均离群率是多点数据集中偏离大部分数据的一种算法。
3.3压力分布数据耦合度
天然气输送管道的渗透率和孔隙度等相关参数会随着输送压力的变换而发生动态变化，同时该变化又会影响到压力的分布状况，针对压力分布数据耦合度如表2所示。
由表2可知，数字孪生体算法仿真在压力分布数据耦合度统计中的决定系数验证结果略高于BIM自带的结果，在标准偏差率和平均离群率的验证结果却明显低于BIM自带的数据结果。
3.4数据前推估计敏感度
对数字孪生体算法仿真和BIM自带数据系统依据上述的系统数据耦合度验证结果数据曲线，进行系统的数据前推估计敏感度统计，如表3所示。
由表3可知，显示数字孪生体算法仿真在5s和曲线估计前推数据结果略高于BIM自带的结果，在15、30、60s时的估计前推数据均明显优于BIM自带的敏感度结果，且比较结果为t<10.000，P<0.05，说明具有明显统计学意义。
4结语
通过在MATLAB中加载Simulink控件构建仿真环境，对流量分布和压力分布的数据耦合度，以及数据前推估计敏感度进行了测试，结果证明了数字孪生体算法仿真能够有效的对异常数据预警和主动推送，能够提高系统应用的安全性，具有很好的应用价值。
【参考文献】
[1]余宏伟，高小霞.天然气输送管道事故应急响应现存问题及解决方案[J].化工管理，2021（25）：125126.
nlc202210121127
[2]曹玉法，王洪申，郑成凯.基于数字孪生概念的实验探究[J].机械制造与自动化，2021，50（4）：3234.
[3]张卓雷.数字孪生助力未来智慧城市新基建[J].信息化建设，2021（9）：3032.
[4]杨毅，周志斌，李长俊，等.天然气管输调节控制仿真模型[J].天然气工业，2008（10）：98100.
[5]曹学文，曹恒广，赵湘阳，等.直板清管器与管壁间摩擦阻力的模拟研究[J].石油机械，2021，49（5）：138146.
[6]崔蓝予，刘贵军，黄凯良，等.高层住宅厨房集中排烟道弯头阻力模拟分析[J].建筑与预算，2021（4）：5355.
[7]叶建君，吴凡，林俊光，等.燃气蒸汽联合循环电厂天然气系统阻力计算及设计[J].节能，2021，40（4）：1920.
[8]朱杰，匡海鹰，王海明，等.次氯酸钠在长距离原水管中辅助减阻的工程效应[J].净水技术，2021，40（4）：4348.
[9]慕青松，汪儒生，牛永红，等.地下管道受冻土冻胀影响和轴力作用的力学分析[J].兰州大学学报（自然科学版），2021，57（2）：278284.
[10]张加蓉，朱桥，高嵩，等.饱和地热蒸汽管道长距离水力计算研究[J].电力勘测设计，2021（3）：6166.
[11]戎宗义.用Excel计算烟道和空煤气管道的阻力[J].工业炉，2021，43（2）：4553.
[12]赵峰辉，王Y，关旭辉，等.天然气长输管道节能设计探讨[J].轻工科技，2021，37（7）：1819.
[13]戚伟，郭沫川，谭玉叶，等.大冶铁矿充填管道自流输送参数分析及管网磨损规律研究[J].矿业研究与开发，2021，41（6）：164169.
[14]张白云，王俊新，唐泽润，等.长距离输水泵站中空气罐进出口阻力系数对其水锤防护效果影响的研究[J].中国农村水利水电，2021（6）：197201.
[15]刘一鸣，姜晗琳，李振山.黄河下游泥沙长距离管道输送目标优化方法[J].人民黄河，2021，43（5）：3034.
[16]曹学文，曹恒广，赵湘阳，等.直板清管器与管壁间摩擦阻力的模拟研究[J].石油机械，2021，49（5）：138146.
[17]崔蓝予，刘贵军，黄凯良，等.高层住宅厨房集中排烟道弯头阻力模拟分析[J].建筑与预算，2021（4）：5355.
[18]叶建君，吴凡，林俊光，等.燃气-蒸汽联合循环电厂天然气系统阻力计算及设计[J].节能，2021，40（4）：1920.
[19]朱杰，匡海鹰，王海明，等.次氯酸钠在长距离原水管中辅助减阻的工程效应[J].净水技术，2021，40（4）：4348.
[20]慕青松，汪儒生，牛永红，等.地下管道受冻土冻胀影响和轴力作用的力学分析[J].兰州大学学报（自然科学版），2021，57（2）：278284.
[21]张加蓉，朱桥，高嵩，等.饱和地热蒸汽管道长距离水力计算研究[J].电力勘测设计，2021（3）：6166.
[22]戎宗义.用Excel计算烟道和空煤气管道的阻力[J].工业炉，2021，43（2）：4553.
[23]章书远，魏新江，王霄，等.顶管技术顶进阻力的研究现状与展望[J].低温建筑技术，2021，43（3）：105110.
[24]孙志明.探析横管式煤气初冷器阻力大的原因及应对措施[J].中小企业管理与科技（下旬刊），2021（3）：158159.
[25]杜东明，谭灿觯张翔宇.1000MW机组烟风煤粉管道节能降耗探讨[J].南方能源建设，2021，8（1）：8792.
nlc202210121127
转载注明来源:https://www.xzbu.com/1/view-15440885.htm
