埋地管道阴极保护系统的优化分析与研究.pdf

2023年8月防腐保温DOI：10.3969/j.issn.1001-2206.2023.04.009埋地管道阴极保护系统的优化分析与研究丁淑娟1，左婧2，秦虹3，胡俊平2，郭自强2，李盼31.中国石油华北油田公司第三采油厂，河北河间 0624502.中国石油华北油田公司第五采油厂，河北辛集 0523603.中国石油华北油田公司第一采油厂，河北任丘 062552摘要：为优化埋地管道的阴极保护系统使管道全线得到更好的保护效果，在电化学实验的基础上，利用极化曲线确定待评价管道的阴极保护区间，通过COMSOL边界元仿真实验确定影响保护电位分布特征的主控因素，建立了以阳极数量、阳极位置、输出电流为决策变量的优化函数，并利用SAPSO算法进行多优化方案的求解。结果表明，可将-1.08-0.87 mV作为阴极保护电位区间；阳极数量越多、阳极与管道的径向间距越大、输出电流越大、土壤电阻率越小，管道全线的保护效果越好；在保持原有阳极地床位置不变，且将阳极2更换为深井阳极的情况下，管道全线的保护效果和电位均匀性最好，且不涉及二次征地，对其余已有阴极保护系统也不构成影响。关键词：埋地管道；阴极保护；阳极；输出电流；SAPSOOptimization analysis and research of cathodic protection system for buriedpipelineDING Shujuan1,ZUO Jing2,QIN Hong3,HU Junping2,GUO Ziqiang2,LI Pan31.No.3 Oil Production Plant of Huabei Oilfield Company,CNPC,Hejian 062450,China2.No.5 Oil Production Plant of Huabei Oilfield Company,CNPC,Xinji 052360,China3.No.1 Oil Production Plant of Huabei Oilfield Company,CNPC,Renqiu 062552,ChinaAbstract：This paper aims to optimize the cathodic protection system of buried pipelines,so as to make the whole pipeline line achievea better protection effect.On the basis of electrochemical experiments,the cathodic protection area of the pipeline to be evaluated wasdetermined by the polarization curve,and the main controlling factors affecting the distribution characteristics of protection potentialwere determined by COMSOL boundary element simulation experiment.The optimization function with anode number,anode position,and output current as decision variables was established,and the SAPSO algorithm was used to solve the multi-optimization scheme.The results show that 1.080.87 mV can be used as the cathodic protection potential range.As the number of anodes gets greater,theradial distance between the anode and pipe becomes larger.In addition,the output current increases,and the soil resistivity decreases.The whole line of the pipeline is better protected.By keeping the position of the original anode bed unchanged and replacing anode 2with the deep-well anode,the potential protection and uniformity of the whole line of the pipeline are the best.No secondary landacquisition is involved,and the remaining cathodic protection system is not affected.Keywords：buried pipeline;cathodic protection;anode;output current;SAPSO对于在役油气管道，阴极保护系统与防腐层一起构成了埋地管道的防护体系，使管道免受外腐蚀及杂散电流的影响1-2。目前，诸多学者针对阴极保护系统的运行特点和参数优化进行了研究，陈涛涛等3针对某燃气门站的阴保系统，结合馈电实验结果，总结了不同阳极地床下的优化方案；赵悦春等4通过研究电流开放系统的分布规律和保护效果，得到了适合低压输气管道的阳极地床形式；于淑珍等5通过断电实验考察了不同管道共用阴极保护系统之间的影响情况，确定了新建阳极地床的位置；冯朋鑫等6通过电连续排查和绝缘性排查测试了联合阴极保护的效果，针对阴极智能系统的升级提出了优化方案。以上研究多通过现场馈电实验及通/断电实验，测试阴极保护的保护效果，阴保系统的优化方案均基于人工经验，导致沿线难以达到较均匀的电位分布，造成现场反复施工，浪费了大量人力和物力。基于此，在电化学实验的基础上，利用极化曲线确定待评价管道的阴极保护电位区间，通过 COMSOL 边界元仿真实验确定影响保护电位分布特征的因素，利用42第49卷第4期防腐保温模拟退火粒子群算法（SAPSO）对阳极地床的参数进行优化，最后通过实例验证使电位和电流分布达到最佳。1实验部分利用PARATAT 2273电化学工作站进行实验，试片与待评价管道材质同为X60管道钢，试片尺寸10 mm 10 mm 5 mm，除工作面外均用环氧树脂密封。工作电极为试片，参比电极为饱和甘汞电极，辅助电极为铂电极，腐蚀介质为模拟溶液（按管道沿线土壤成分配置），试片在三电极体系下测试其极化曲线。动电位扫描范围为-1.20 V（SCE），扫描速率为1 mV/s。图 1 为 X60 钢的极化曲线。从阴极极化曲线上看，试片表面和腐蚀溶液界面之间发生氧化还原反应，电流密度随极化电位的增加而增加。在-1.08-0.87 mV内，腐蚀电位与电流密度呈线性关系，试片表面的阴极保护效果不断增强；当大于-1.08 mV时，电流密度突然增加，说明金属表面的负电离子较多，会引发管材发生氢脆反应。因此，将-1.08-0.87 mV作为阴极保护电位区间，该区间与GB/T 214482017埋地钢质管道阴极保护技术规范中建议的区间（-1.2-0.85 mV）相比更加精确，可为下一步系统优化提供数据支撑。-0.5-0.6-0.7-0.8-0.9-1.0-1.1腐蚀电位/V-7-6-5-4-3-2-1lg（电流密度）/（A/cm2）图1X60钢的极化曲线2模拟部分鉴于有限元法在计算迭代次数、边界条件和网格质量上的局限性，采用模型简化的边界元法进行模拟实验，得到影响阴保电位分布的主要因素。2.1控制方程和边界条件采用“电流分布、边界元”接口求解边缘区域的电解质电势和电流密度，公式如下：i1=-（1）i1=0（2）式中：i1为土壤中的电流密度，A/m2；为土壤电导率，s/m；为土壤中的电解质电势，V。在阳极地床附近，采用电解质电流密度节点确定外加电流密度。ni1=i2（3）式中：n为法向向量，指计算域以外；i2为外加电流密度，A/m2。将图1中阴极极化曲线分段插值后作为埋地管道的边界条件，埋地管道采用固定电位模型，电极电位为0 V。2.2模型建立和网格划分由于该部分只是考察影响电位分布的因素，故采用COMSOL进行简易模拟7-8。模型中管道材质为 X65 钢，长 1 600 m，规格 762 mm 12 mm，埋深2 m，防腐层电阻率为105m；辅助阳极为高硅铸铁，长5 m，直径0.1 m，埋深3 m；阳极位于管道中心处，距离管道径向距离为100 m，坐标为（0，-100，-3），土壤为无限空域。网格中最大单元格尺寸为8 m，共形成203个网格顶点、201个边单元和4个顶点单元，稳态求解精度为0.001。模型和网格划分结果见图2。-50-45000-500图2模型及网格划分结果/m2.3模型求解考察不同阳极数量对阴极保护电位的影响，结果见图3。2组阳极的坐标分别为（500m，-100m，-3 m）和（-500 m，-100 m，-3 m），3组阳极的坐标分别为（400 m，-100 m，-3 m）、（0 m，-100 m，-3 m）、（-400 m，-100 m，-3 m），阳极坐标对称分布。由于采用无限空域作为土壤，因此在土壤边界会出现电位下降的现象，但不影响分析结果。当阳极数量为 1 时，靠近阳极地床的管道电位呈漏斗状，电位从管道中心向两端逐渐增加，所有管段处于欠保护状态；当阳极数量为 2 或 3时，出现了与阳极数量相同的漏斗电位，且整丁淑娟等：埋地管道阴极保护系统的优化分析与研究432023年8月防腐保温体电位有所下降，阳极数量越多，电位分布越均匀。由于外加电流密度保持不变，因此较多的阳极数量有可能突破析氢电位，造成防腐层阴极剥离。在3组阳极的条件下，依次变化阳极坐标中的y和z（代表与管道的径向间距和埋深），考察阳极位置对阴极保护电位的影响，结果见图4。最大保护电位1 m3 m5 m7 m管道长度/m200400600800-800-600-400-2000电极电位/mV-0.82-0.90-0.97-1.05-1.12-1.20-1.27-1.35-1.43（b）不同阳极埋深（a）不同径向间距最小保护电位电极电位/mV-0.78-0.81-0.84-0.87-0.90-0.93200400600800-800-600-400-2000管道长度/m20 m40 m60 m80 m100 m图4阳极位置对阴极保护电位的影响当径向间距为20 m时，阳极地床附近的管道电位迅速负移，且移动速率随间距的减少呈指数增加；在径向间距为20、40 m时，管道电位已经超过最大保护电位。距离阳极越远，电位分布越均匀，但同时电位也会正移，导致施加在管道上的阴极电流不够，存在欠保护风险。随着阳极埋深的增加，管道电位先负移后正移，在埋深较浅时，管道全线无法达到完全保护的效果，在埋深较深时，部分管段同样存在欠保护现象。因此，对于阳极位置存在最佳的参数选择。在3组阳极的条件下，考察恒电位仪输出电流值对阴极保护电位的影响，结果见图5。随着输出电流的增加，保护电位呈不等间距负方向移动，且保护区间变大，电位分布均匀性变差，当输出电流为0.5 A时，管段全线为欠保护状态；当输出电流为2 A时，保护性最好。0.5 A1.0 A1.5 A2.0 A管道长度/m200400600800-800-600-400-2000电极电位/mV-0.85-0.87-0.89-0.90-0.92-0.93-0.95图5输出电流值对阴极保护电位的影响在3组阳极的条件下，考察土壤电阻率对阴极保护电位的影响，结果见图6。200 m100 m50 m管道长度/m200400600800-800-600-400-2000电极电位/mV-0.89-0.90-0.92-0.93-0.95-0.96图6土壤电阻率对阴极保护电位的影响随着土壤电阻率的上升，保护电位整体等间距正移，电位均匀性不发生改变，电阻率越高，土壤图3不同阳极数量对阴极保护电位的影响最小保护电位1组阳极2组阳极3组阳极管道长度/m200400600800-800-600-400-2000电极电位/mV-0.85-0.87-0.89-0.90-0.92-0.9344第49卷第4期防腐保温的IR降越大，能到达管道的保护电流越小，保护效果越差。另外，过大的土壤电阻率会使阴保电位达到析氢电位，产生氢脆风险。3基于SAPSO算法的阴极保护参数优化通过以上仿真结果可知，阳极数量、阳极位置、输出电流和土壤电阻率等参数决定着管道阴极保护系统的有效性和实用性9。考虑到土壤电阻率需根据沿线阳极地床的埋设地点，通过四电极法测量，故不将此参数作为优化变量，只作为固定参数输入。此外，管道外防腐层破损率也与阴极保护的电位分布有关，破损率越大，电位分布越不均匀，保护性越差，破损的地方互成阴阳极，加速腐蚀。当管道处于保护状态时，管道表面的阴极保护电位应在最佳保护电位附近，且不超过保护电位区间，将管道表面离散成n个点，则管道表面的平均电位-的计算公式如下：-=1nj=1nj（4）min 1Ie,De(x,y,z),ne=1n1n()j-2（5）min 2Ie,De(x,y,z),ne=()p-2（6）式中：j为j节点的电位值，mV2；1（）为管道电位与平均电位之间的差值，反映了电位分布的均匀程度，mV；Ie为输出电流，A；De（x，y，z）为阳极地床的三个位置参数，m；ne为阳极数量，个；p为最佳保护电位，mV，在此取-0.95 mV；2（）为最佳保护电位与平均电位之间的差值，反映了电位优化的效果，mV2。1、2均为目标函数，两者越小，阴极保护效果越好，反之越差。将阴极保护参数优化的式（5）、式（6）转化为多目标寻优问题，采用模拟退火（SA）和粒子群算法（PSO）的结合进行求解，算法原理见文献10-11。求解过程分3步。1）初始化粒子的各项参数，包括 Ie、De（x，y，z）和ne，设置土壤电阻率和防腐层破损率，以式（5）、式（6）为适应度函数，通过式（7）、式（8）计算粒子的速度和位置信息。vi,d（k+1）=wvi,d（k）+c1r1pi,d（k）-zi,d（k）+c2r2pg,d（k）-zi,d（k）（7）zi,d（k+1）=zi,d（k）+vi,d（k+1）（8）式中：i为粒子序号；d为维度序号；k为当前迭代次数；w、c1、c2分别为惯性权重和加速因子；v、z分别为位置和速度；r1、r2均为0到1的随机数；pi,d和pg,d分别为个体最优位置和全局最优位置。2）根据 Metropolis 法则，当适应度值的增量E 0，则新位置的解较优，选择新解；否则，按照概率 p=exp（E/T）选择较优解（式中：T为退火温度）。3）查找终止迭代条件是否满足 0.000 1 的阈值，如果满足则退出当前迭代；否则，以T为控制参数更新粒子的各项参数，重复步骤1）2）。此外，式（7）中的w、c1、c2均为定值，无法平衡全局和局部搜索的要求，故采用双正切函数构建随迭代次数变化的动值参数，对SAPSO算法进行改进，公式如下：wk=wa+tanh(-8+16(K-k)/K)ws/2c1,k=ca+tanh(-8+16(K-k)/K)cs/2c2,k=ca-tanh(-8+16(K-k)/K)cs/2（9）式中：wk、c1，k、c2，k均为对应第k次迭代的数值；K为总的迭代次数；wa、ca为该参数的最大值和最小值之和；ws、cs为该参数的最大值和最小值之差。4实例分析4.1理论计算以某输油管道为例，该管道全长143.5 km，管径D355.6 mm，壁厚8 mm，管材为X65，采用石油沥青防腐层，防腐层厚度5 mm。根据防腐层绝缘电阻，通过简易公式计算保护电流密度为30 A/m2，再利用式（10）计算外加电流阴极保护系统的单侧保护长度L12。2L=8VL/（DJsR）=49.3（km）（10）式中：VL为保护电位差值，V；D为管道外径，m；Js为保护电流密度，30 A/m2；R为单位长度管道纵向电阻，根据管道电阻率、管径和壁厚计算，20.64 10-6/m。由此得到管道沿线阴极保护站的数量N=L总/丁淑娟等：埋地管道阴极保护系统的优化分析与研究452023年8月防腐保温2L+1（L总为管道总长度，m），经计算N 4；阴极保护电流Ie=DJsL（Ie为外加阴极保护电流，A），经计算Ie=0.83 A。根据SY/T 00362016埋地钢质管道强制电流阴极保护设计规范 和 Q/SY06314.12016油气储运工程防腐绝热设计规范第1部分：管道防腐层及阴极保护的相关规定，要求设计的阳极埋深不小于1 m。最终根据经验，设计浅埋阳极的埋深为3 m，距离管道径向距离为70 m。利用管道沿线的阴极保护桩，实测管道保护电位，见图7。阳极3管道里程/km100120140204060800高程/m40200-20-40阴极保护电位/V-0.7-0.8-0.9-1.0-1.1-1.2-1.3阳极4阳极2阳极1管道高程模拟阴极保护电位实测阴极保护电位图7管道沿线阴极保护电位目前，只有阳极3所在的管段满足阴极保护电位要求，其余管段在阳极附近的电流过大，而远离阳极的管段又明显保护不足。从该管道的维抢修记录中查阅，发生外腐蚀最严重的管段为4055 km之间，与实测保护电位相符。此外，通过COMSOL建模，得到实际管道的电位模拟结果，与实测电位差距较小，证明了模拟结果的准确性。4.2优化结果4.2.1优化方案1通过 SAPSO 算法对阴极保护参数进行优化，Ie、De（x，y，z）、ne的初始值分别采用以上常规设计结果，最终得到5个阳极地床的参数信息，坐标分别为（20.6 km，-68 m，-3.7 m）、（43.8 km，-55 m，-4.1 m）、（71.2 km，-59 m，-4.2 m）、（93.1 km，-31 m，-1.9 m）、（117.7 km，-62 m，-2.1 m），保护电流均为 0.95 A。从算法迭代角度看，初次迭代时，优化结果随机分配，较高的退火温度使粒子在迭代过程中表现出随机化和多样化；随着迭代次数的增加，较低的退火温度使粒子逐渐向最优解的位置靠拢，阳极地床的坐标位置逐渐清晰，见图8。图8不同迭代次数下的收敛结果（b）第50次迭代（a）初次迭代-2-3-4-5-6-7-8-9-10z/m-20-40-60y/m-80-100140120100806040200 x/kmz/my/mx/km-2.0-2.5-3.0-3.5-4.020406080100120140-80-60-40-20根据优化方案 1 定义的阴极保护参数，采用COMSOL对其进行模拟，见图9。与实际运行方案相比，优化后的阴极保护电位均匀分布在-1.08-0.87 mV 区间内，且始终围绕在最佳保护电位-0.95 mV附近，证明了优化算法的可靠性，管道全线的保护效果良好。方案1阴极保护电位实测阴极保护电位204060801001201400管道里程/km原阳极2原阳极4新阳极5新阳极4原阳极3新阳极3新阳极2新阳极1原阳极1阴极保护电位/mV-0.8-0.9-1.0-1.1-1.2-1.3图9优化方案1的阴极保护电位46第49卷第4期防腐保温4.2.2优化方案2将已有4个阳极地床的位置固定，采用SAPSO算法优化现有阳极恒电位仪的保护电流和新建阳极地床位置。优化后原4个阳极的I1=0.75 A，I2=0.78 A，I3=0.85 A，I4=0.73 A；新建阳极地床的位置为（53.2km，-73m，-3.2m），保护电流I5=1.05A，结果见图10。与实际运行方案相比，管道全线保护效果较差，阳极2位置处接近最大保护电位。204060801001201400管道里程/km原阳极2原阳极4原阳极3新阳极2新阳极原阳极1阴极保护电位/mV方案2阴极保护电位-0.8-0.9-1.0-1.1-1.2-1.3实测阴极保护电位图10优化方案2的阴极保护电位4.2.3优化方案3将原有阳极1、3、4的位置固定，采用SAPSO算法优化现有阳极恒电位仪的保护电流，将位置2的阳极地床换为深井阳极，深井阳极由4个贵金属阳极分段构成，每段长 6 m，共 24 m，阳极直径0.2 m。优化后原有3个阳极的I1=0.71 A，I3=0.84 A，I4=0.72 A；深井阳极的I2=1.73 A，深井阳极的位置为（71.1 km，-59 m，-30.1 m），结果如图 11 所示。与实际运行方案相比，管道全线保护效果良好，管道沿线电位分布在最佳保护电位附近，且均匀性较好。204060801001201400管道里程/km原阳极2原阳极4原阳极3原阳极1阴极保护电位/mV方案3阴极保护电位-0.8-0.9-1.0-1.1-1.2-1.3实测阴极保护电位图11优化方案3的阴极保护电位4.2.4方案对比综上所述，方案 1 的优化效果较好，除增加一处阳极地床外，其余阳极地床的位置均有所变化，如对其重新调整埋深，需增加现场施工次数，造成二次开挖；方案 2 的优化效果较差，除新增阳极地床外，当阳极 2 的恒电位仪出现电流波动时，容易造成析氢风险；方案3的优化效果最好，深井阳极的位置与原位置相同，只有阳极材料和埋深发生了改变，此时只需对阳极 2 进行整改更换即可，不涉及二次征地，对其余已有阴极保护也不构成影响。因此，可选择方案 3 作为整改优化方案。5结论1）通过电化学实验，利用阴极极化曲线确定了待评价管道的最佳保护电位区间。2）针对阴极保护电位的分布特征，通过COMSOL软件进行了多因素仿真实验，阳极数量越多、阳极与保护管道的径向间距越大、输出电流越大、土壤电阻率越小，管道全线的保护效果越好；随着阳极埋深的增加，管道电位先负移后正移，保护效果有所差异。3）利用 SAPSO算法对阳极数量、阳极位置、输出电流等参数进行了优化，与实际运行方案相比，优化后的电位分布更加均匀，可减少二次开挖和阳极重新布置，具有重要意义。参考文献1 马贵利，刘海峰，王浩，等.长输管道站场区域阴极保护馈电试验结果分析J.天然气与石油，2022，40（6）：94-100.2 肖剑烽，张争伟，唐振峰.完整性管理技术在衢州压气站区域阴极保护系统的应用J.油气田地面工程，2022，41（11）：87-92.3 陈涛涛，赵茜，杜艳霞，等.某燃气门站区域阴极保护优化方案确定J.腐蚀与防护，2021，42（11）：86-91，106.4 赵悦春，庄大伟，杜艳霞，等.城市小区低压燃气管网开放式阴极保护方式及地床形式优化J.腐蚀与防护，2022，43（9）：93-98.5 于淑珍，刘明，韩鹏刚，等.某气田集输管线阴极保护检测与优化设计J.材料保护，2021，54（7）：146-149，154.（下转第54页）丁淑娟等：埋地管道阴极保护系统的优化分析与研究472023年8月防腐保温参考文献1 张维.多因素综合作用下高寒地区输油管道剩余寿命的预测J.辽宁化工，2022，51（11）：1 556-1 560，1 564.2 骆正山，李易安.基于GRA-FOA-ELM的腐蚀管道剩余寿命预测J.消防科学与技术，2020，39（4）：529-532.3 张新生，张平.不完全维修下海底腐蚀管道剩余寿命预测J.系统工程理论与实践，2019，39（11）：2 984-2 994.4 王观军，韩庆，徐树伟，等.基于相邻内腐蚀缺陷交互作用的油田注水管道剩余寿命研究J.安全、健康和环境，2021，21（7）：40-46.5 王战辉，张智芳，李瑞瑞.腐蚀缺陷对高压管道剩余强度和剩余寿命的影响J.河南科学，2020，38（11）：1 742-1 746.6 黄港港，周阳，闫骁瑾，等.基于KPCA-CS-SVM的埋地管道外腐蚀速率预测研究J.热加工工艺，2022，51（16）：38-43.7 张新生，常潆戈.基于FA-BAS-ELM的海洋油气管道外腐蚀速率预测J.中国安全科学学报，2022，32（2）：99-106.8 齐盛，帅健，顾强.基于贝叶斯和极值分析的原油储罐底板剩余寿命预测J.中国安全生产科学技术，2021，17（S1）：67-73.9 张杰，张佩颖，虞维超，等.引入失效数据的腐蚀管道结构可靠 性 评 价 方 法 J.安 全 与 环 境 工 程，2020，27（3）：185-192.10 巴振宁，王鸣铄，梁建文，等.改进的一次二阶矩计算埋地管道失效概率方法及其应用J.安全与环境学报，2021，21（6）：2 371-2 377.11 丁志林，胡风，余春波，等.基于Wiener过程的高压晶闸管寿命预测方法J.武汉大学学报（工学版），2022，55（10）：1 044-1 049.12 武文栋，施保华，郑传良，等.基于KPCA-ISSA-KELM的光伏阵列故障诊断方法J.智慧电力，2022，50（11）：69-76.13 张熙.基于Wiener过程首超概率的结构耐久性评估J.重庆建筑，2022，21（6）：50-53.14 王泽洲，陈云翔，蔡忠义，等.基于比例关系加速退化建模的设备剩余寿命在线预测J.系统工程与电子技术，2021，43（2）：584-592.15 王艺斐，苏春，谢明江.基于二元逆高斯过程的腐蚀输油管道剩余寿命预测J.东南大学学报（自然科学版），2020，50（6）：1 038-1 044.作者简介：宋正涛（1976），男，山东青岛人，工程师，2005年毕业于长江大学计算机科学与技术专业，现主要从事生产运行管理工作。Email：cy3_收稿日期：2023-05-016 冯朋鑫，宋考平，于淑珍，等.苏里格气田集输管道阴极保护运行检测及优化案例分析J.腐蚀科学与防护技术，2018，30（6）：676-682.7 高荣杰，郑斐，王焕磊，等.COMSOL软件模拟在腐蚀与防护教学中的应用J.大学教育，2020（11）：109-111，125.8 辛艳萍，梁月.基于Comsol Multiphysics的金属储罐阴极保护方案的优化J.石油化工高等学校学报，2016，29（4）：92-96.9 张奇志，雷佳，李琳.基于改进遗传算法的管道阴极保护系统辅 助 阳 极 位 置 优 化 J.材 料 保 护，2020，53（1）：84-90，105.10 刘英伟，周子杰，张洋.一种简单有效的船体外加电流阴极保护设计方法J.表面技术，2019，48（7）：309-315.11 骆正山，潘柯成.基于LASSO-ISAPSO-ELM的含蜡原油管道蜡沉积速率预测J.安全与环境工程，2022，29（6）：69-77.12 孟涛，程喜兵.城市埋地燃气钢管防腐措施研究J.管道技术与设备，2013（1）：38-40.作者简介：丁淑娟（1970），女，甘肃张掖人，工程师，2016年毕业于中国石油大学（华东）油气储运专业，现主要从事油气集输管理工作。Email：cy3_收稿日期：2023-05-05（上接第47页）54