Modeling and Simulation of Inductive Coupling Interference from Suburban Railways to Buried Pipelines
-
摘要:
随着市域铁路和油气管道的快速发展,铁路与埋地管道平行或交叉铺设的情况已不可避免. 为分析市域铁路牵引供电系统对油气管道安全运行的影响,首先,建立铁路交流牵引供电系统与邻近埋地管道的数学模型,并利用CDEGS软件进行仿真验证;然后,采用等效距离方法将斜接近和平行接近2种情况进行统一,探究牵引供电系统对邻近埋地管道的影响机理,分析土壤电阻率、导线与埋地管道间接近距离、机车负荷电流、铁路与埋地管道间并行长度、管道涂层电阻率、动车组电流谐波次数等因素对管道沿线感应耦合电压分布的影响;最后,结合贯通地线的设置,提出4种抑制感应耦合电压的方案进行对比分析. 研究结果表明:当市域铁路牵引供电系统与埋地管道的最大间距与最小间距比值小于4.5时,采用等效距离方法的偏差在5%以内;管道上的感应耦合电压极大值随着土壤电阻率的增加而增大;接近距离由50 m增大到250 m时,感应耦合电压降幅为50.6%;机车负荷电流在200~
1000 A变化时,感应耦合电压增幅较大;并行长度在2~10 km变化时,感应耦合电压从22.6 V增加到170.7 V;谐波含有率和谐波次数对感应耦合电压的影响不可忽视;在贯通地线基础上增设一根回流线的方案对抑制感应耦合电压的效果最佳.Abstract:With the rapid development of suburban railways and oil and gas pipelines, parallel laying or cross-laying is inevitable for suburban railways and buried pipelines. In order to assess the impact of the traction power supply system of suburban railways on the safe operation of oil and gas pipelines, firstly, the mathematical model of the alternating current (AC) traction power supply system of the suburban railways and the adjacent buried pipelines was established, and the simulation based on CDEGS software was conducted. Next, the two cases of the oblique approach and parallel approach were unified by the equivalent distance method. The influence mechanism of the traction power supply system of the suburban railways on the adjacent buried pipeline was investigated, and the influence of factors including soil resistivity, distance between the conductor and the buried pipeline, locomotive load current, parallel length of railway and buried pipeline, pipeline coating resistivity, and number of current harmonics of electric multiple units on the inductive coupling voltage distribution along the pipeline was investigated. Finally, combined with the setting of a through ground line, four schemes to suppress the inductive coupling voltage were put forward for comparative analysis. The results show that the error of the equivalent distance method is within 5% when the ratio of the maximum distance to the minimum distance between the traction power supply system of the suburban railway and the buried pipeline is less than 4.5. The maximum value of the inductive coupling voltage of the pipeline increases with increasing soil resistivity. The decrease in the inductive coupling voltage is 50.6% when the distance between the suburban railway and the buried pipeline varies from 50 m to 250 m. The increase in the inductive coupling voltage rises significantly when the locomotive load current varies from 200 A to 1 000 A. The inductive coupling voltage increases from 22.6 V to 170.7 V when the parallel length varies from 2 km to 10 km. The harmonic content and the number of harmonics have a significant influence on the inductive coupling voltage. The best inductive coupling voltage suppression effect is achieved by adding a return line on the basis of the through ground line.
-
Key words:
- suburban railway /
- buried pipelines /
- inductive coupling voltage /
- oblique approach /
- harmonics
-
光伏(photovoltaic,PV)阵列因建筑物、树木、云层等遮挡而处于局部阴影条件(partially shaded conditions,PSC),产生热斑效应[1]. 为避免热斑效应的产生,光伏组件两端通常并联旁路二极管,这将导致局部阴影条件下光伏阵列的功率-电压特性曲线呈现多峰特性. 扰动观察法、电导增量法等传统最大功率点跟踪(maximum power point tracking,MPPT)算法[2-3] 容易陷入局部最大功率点(local maximum power point,LMPP),造成能量损失,不适用于局部阴影条件. 为此,改进的MPPT技术被陆续提出,并促进了光伏发电技术的快速发展[4-17].
针对局部阴影条件下的MPPT可以分为硬件方法与软件方法两类. 硬件方法通过改变光伏发电系统的结构以减少局部阴影导致的能量损失. 文献[6-7]在光伏阵列与变换器之间引入可控的开关矩阵,实时改变阵列中光伏组件的组合方式以适应阴影条件;文献[8]通过增加直流变换器使光伏组件独立运行,且均工作在MPPT状态. 尽管硬件方法可有效解决传统MPPT算法在局部阴影条件下的误跟踪问题,但也增加了电路结构的复杂度和系统成本. 软件方法,又称全局最大功率点跟踪(global maximum power point tracking,GMPPT)算法,不需要额外增加硬件电路,即可识别多峰曲线的全局最大功率点(global maximum power point, GMPP) ,相较于硬件方法展现出明显的成本优势. GMPPT算法主要分为智能算法、快速功率线算法、基于传统搜索算法的改进型GMPPT算法等. 智能算法,如粒子群算法和模糊逻辑控制[9-10],能够精确地从多个极值中找到最大值,实现优化目标,但算法较为复杂. 快速功率线算法则以Direct算法[11]和Fibonacci算法[12]较为常见,通过利用光伏阵列输出曲线的特征缩小搜索范围,快速定位GMPP,然而,若参数初始值设置不当,可能导致算法收敛到局部功率点. 基于传统搜索算法的改进型GMPPT算法原理简单,且易于实现. 文献[13]在已知光伏组件开路电压的情况下,控制算法扫描特定的电压点附近,有效缩短了扫描区间,但算法可靠性低;文献[14]分析光伏阵列中分段处的短路电流与峰值点电流的近似比例关系,提出基于光伏阵列电流特性的GMPPT算法,但不适用于曲线峰值点从左向右依次升高的情况;最大功率梯形算法[15] (maximum power trapezium,MPT)通过限制扫描电压的上下限,减少搜索范围,提高跟踪速度,但其扫描范围仍然偏大.
尽管改进型GMPPT算法在性能和成本上取得了较好的平衡,但其跟踪速度受限于扰动步长和扰动时间. 扰动时间一定时,扰动步长越小,算法的跟踪精度越高,但在扫描区间内的扰动次数增加,跟踪时间越长,跟踪速度越慢. 扰动步长一定时,扰动时间越小,跟踪时间越短,跟踪速度越快. 在电压闭环控制中,对于相同扰动步长,参考电压在电流-电压特性曲线恒流区所需的调整时间比恒压区要长[16-17]. 因此,现有算法多以恒流区为基准设计扰动时间.
为进一步提高局部阴影条件下光伏发电系统的跟踪速度,本文提出一种改进型快速GMPPT算法. 该算法通过在跟踪过程中跳过调整时间较长的恒流区,提高系统的跟踪速度. 本文通过仿真和实验验证所提算法的有效性,并将其与现有的MPT算法和全局扫描算法的跟踪效果进行对比.
1. 局部阴影条件下光伏组串的输出特性
光伏阵列由多个光伏组件组合而成,分为并联型光伏阵列、串联型光伏阵列与串并联型光伏阵列. 其中,并联型光伏阵列由于输出电压较低,不适用于光伏发电系统. 相比之下,串并联型光伏阵列可以实现较大的输出电流和较高的输出电压,但由于电压和电流之间相互影响,在局部阴影条件下其输出功率会显著降低. 而光伏组串作为串联型光伏阵列,输出电压可根据需求进行设置,在建筑集成光伏系统中应用广泛[18]. 图1为光伏组串的电路结构,图中,PV1、PV2和PVn为光伏组件,Ipv、Vpv分别为组串的电流和电压.
图2为光伏组件的输出特性. Ipv-Vpv 曲线从膝点(knee point,KP)处被划分为恒流区和恒压区. Ppv (组串功率)-Vpv 曲线表现出单峰特性,最大功率点(maximum power point,MPP)与KP相对应. MPP对应的电压VMPP和电流IMPP近似为[19]
{VMPP=0.76Voc,IMPP=0.9Isc, (1) 式中:Voc和Isc分别为同一条件下光伏组件的开路电压和短路电流.
光伏组件的短路电流近似与光照强度成正比,开路电压近似与光照强度的自然对数成正比,即光照强度改变,短路电流随之改变,而开路电压近似不变[20]. 因此,在任意光照条件下,VMPP近似为标准测试条件(光照强度S =
1000 W/m2,温度25 ℃)下组件开路电压Voc,STC,module的0.76倍,即恒流区的长度为0.76Voc,STC,module.图3为局部阴影条件下光伏组串的输出特性. 在Ipv-Vpv 曲线中,随着KP个数的增加,恒流区和恒压区也相应增加. Ppv-Vpv 曲线表现出多峰特性,峰值点中的最大值为GMPP,其他峰值点为LMPP,并与KP相对应. 在局部阴影条件下,第i个峰值点的电压VMPPi和电流IMPPi分别近似为[19]
{VMPPi=i−0.24nVoc,array,IMPPi=0.9Isci, (2) 式中:Voc,array为光伏阵列的开路电压,Isci为第i个组件的短路电流,i = 1, 2, …, n.
由于峰值点与膝点相对应,则第i个膝点电压满足式(2). 在光伏组串中,每个组件对应的KP电压不变,恒流区长度近似不变.
2. 全局最大功率点跟踪
2.1 MPT算法
为缩短搜索范围,MPT算法根据GMPP的分布特点定义扫描电压的上限Vrange,max和下限Vrange,min分别为
{Vrange,max=0.9Voc,STC,array,Vrange,min=PGMPPIGMPP,STC,aray, (3) 式中:Voc,STC,array和IGMPP,STC,array分别为标准测试条件下光伏阵列的开路电压和GMPP处对应的输出电流,PGMPP为GMPP处的功率,即最大输出功率.
式(3)中,Voc,STC,array与IGMPP,STC,array均为常数. 因此,在搜索过程中,MPT算法的电压上限Vrange,max固定,而电压下限Vrange,min随着最大输出功率PGMPP的更新而增大,扫描区间相应减小,从而缩短跟踪时间.
图4为MPT算法的跟踪示例,其光伏组串由3个受光照强度不同(条件1~3)的组件串联构成. 图中,绿线为标准条件下MPT算法的功率边界,红线为扫描路径,阴影区域为恒流区,余图同;横坐标上变量为对应点的电压,点e对应的电压Ve = 0.9Voc,STC,array. 光照强度改变时,光伏组串的工作点由点b跳变为点c. MPT算法启动搜索过程,初始化点c为GMPP,且根据式(3)更新Vrange,min为点k对应的电压Vk;然后,MPT算法控制组串运行至Vrange,max工作在点e;再调节系统向左扫描至Vrange,min. 随着扫描过程中GMPP的更新,Vrange,min不断增大,扫描范围相应缩小,最终Vrange,min更新为点f对应的电压Vf. MPT算法的扫描路径为:c→d→e→f.
MPT算法通过限制搜索电压的动态上、下限,减少搜索区间的范围,有效提高跟踪速度. 但MPT算法仅是在开路电压范围内缩短了搜索区间,未在搜索区间内进一步减少搜索过程,搜索区间仍能进一步改善. 同时,光伏组串的GMPP所对应的输出电压较低或光照强度低导致其全局最大输出功率PGMPP较小时,MPT算法的电压下限Vrange,min小,搜索区间仍较大,跟踪时间长.
2.2 改进型快速GMPPT算法
基于传统搜索算法的改进型GMPPT算法通过缩短扰动时间和搜索范围来减少跟踪时间,提高能量的利用率. MPT算法正是通过限定扫描电压的上下限动态地减少搜索区间,缩短跟踪时间. 为避免误扰动,GMPPT算法扰动时间应大于电压扰动后的调整时间,以保证新一轮扰动在上一轮扰动稳定后进行. 但调整时间受光伏阵列的工作区域影响. 当电压扰动工作在恒流区时,其调整时间比电压扰动在恒压区时的调整时间长. 调整时间越长,扰动时间越长,跟踪时间相应增加,能量损失更严重. 本文提出基于MPT的改进型快速GMPPT算法,在MPT算法的动态搜索区间内进一步跳过调整时间较长的恒流区,缩短扫描区间.
图5为改进型快速GMPPT算法流程,主程序包含阴影检测和计数器溢出检测2个条件. 图中,Vold、Iold、Pold分别为初始化电压、电流、功率,V、I、P分别为当前周期检测电压、电流、功率,ΔV为电压步长,ΔP为当前周期与上一个周期的功率差,ILMPP为光伏阵列在峰值点的输出电流. 阴影检测利用功率差∆PGMPP与功率阈值Pthreshold的大小判断光照条件是否发生改变. 计数器溢出检测利用固定时间计数器定期执行GMPPT子程序,避免算法因未检测到光照变化而长时间运行在GMPP外,造成能量损失.
改进型快速GMPPT算法在MPT的电压上下限内扫描,在扫描过程中定位到峰值点时,跨越当前组件相对应的恒流区,其参考电压Vref为
Vref=VLMPP−VCCR, (4) 式中:VLMPP为光伏阵列在峰值点的输出电压;VCCR为恒流区的长度,近似为0.76Voc,STC,module.
受环境因素的影响,恒流区的范围有所波动. 为避免算法在跨越恒流区过程中错过GMPP,本文所取恒流区的长度为0.7Voc,STC,module. 另外,由式(1)和式(2)可知,峰值点电流近似为0.9倍的短路电流,且峰值点的左侧范围电流近似不变,视为恒流区. 因此,算法利用跨越前后的电流大小判断其是否可能错过GMPP. 若跨越恒流区后的当前电流I满足
I<ILMPP0.9, (5) 则算法仍处于恒流区;反之,算法可能错过峰值点,则在当前电压的基础上叠加VCCR/2,缩短跨越区间,并利用式(5)重新判断,以避免错过GMPP.
图6为改进型快速GMPPT算法的跟踪示例. 初始时,光伏阵列稳定工作于点b,光照强度突然改变,工作点由点b跳变为点c. 算法检测到光照条件变化,启动搜索过程. 其中,蓝色区间[0,Vf]、[Vj,Vh]和[Vg,Vd]为光照强度改变后光伏阵列输出特性曲线的恒流区.
改进型快速GMPPT算法的搜索过程分为以下3个阶段.
1) 初始化:初始化GMPP为点c,并更新电压下限Vrange,min为Vk;设置初始参考电压Vref为Vrange,max,调节光伏组串工作在点e,并根据点c、e的功率大小Pc、Pe,判断是否更新GMPP与Vrange,min.
2) 扫描与判断:减小参考电压Vref,调节光伏组串向左扫描至LMPP (点d);比较点d与当前GMPP的功率大小,更新GMPP为点d,且电压下限Vrange,min更新为Vf.
3) 跨越:在扫描到LMPP后,跳过该组件的恒流区;即在Vd的基础上减去恒流区的电压范围,调节系统工作在点g.
在此基础上,不断重复上述过程,直到光伏阵列的输出电压小于最终的电压下限Vrange,min. 最后,设置Vref为最终的GMPP所对应电压,调节系统稳定运行于点d.
综上所述,改进型快速GMPPT算法的扫描路径为:e→d,g→h,j→f.
3. 仿真与实验
3.1 仿真分析
为验证算法的有效性,在MATLAB/Simulink中搭建了仿真模型. 光伏MPPT系统的主电路由光伏组串、变换器、负载电阻组成,如图7所示. 图中:D、D1~D3为二极管,Ci为输入电容,Co为输出电容,L为电感,S为开关管,R为电阻.
为与传统MPT算法[15]进行对比,本文选用了文献[15]中型号为Kyocera KC200GT的3个光伏组件进行串联,构成光伏组串,光照条件如表1所示. 组件参数为:最大功率点电压Vm=26.30 V,最大功率点电流Im=7.61 A,最大输出功率Pm=200.14 W,开路电压Voc=32.90 V,短路电流Isc=8.21 A. 变换器的参数为:Ci = 330 μF,L = 220 μH,Co = 330 μF,开关频率fs = 50 kHz,R = 200 Ω. 表1中列出的不同局部阴影条件下的光伏阵列输出特性如图8所示.
表 1 光伏阵列的输出参数Table 1. Output parameters of the PV array条件 光照强度/(W·m−2) 全局最大功率点 PV1 PV2 PV3 功率/W 电压/V 电流/A PSC1 120 210 400 90.13 54.60 1.65 PSC2 200 300 560 134.60 84.80 1.59 PSC3 150 270 750 143.60 25.20 5.70 图9为改进型快速GMPPT算法(算法1)、全局扫描算法(算法2)和MPT算法(算法3)的仿真波形. 图中,t为时间. 为对比算法的跟踪效果,3种算法的电压步长均为0.2 V,扰动时间为0.8 ms. 仿真跟踪性能对比结果如表2所示.
表 2 阴影条件变化下3种算法的仿真跟踪性能比较Table 2. Simulation tracking performance comparison of the three algorithms under shading variation光照条件 扫描时间/s 总路径的扫描电压/V 能量损失/J 算法 1 算法 2 算法 3 算法 1 算法 2 算法 3 算法 1 算法 2 算法 3 PSC1→PSC2 0.110 0.330 0.340 17.34 69.09 71.14 2.94 7.81 8.25 PSC2→PSC3 0.095 0.295 0.300 17.57 69.09 69.96 4.76 13.66 13.69 PSC3→PSC1 0.125 0.320 0.355 17.34 69.09 76.99 2.63 6.13 7.41 扫描时间为光照条件未突变前的稳定状态到突变后的稳定状态的跟踪时间. 光照条件在0.500 s时由PSC1变化为PSC2,改进型快速GMPPT算法在0.500 s时检测到光照变化,启动搜索过程,约在0.660 s时稳定,跟踪至突变后的GMPP,扫描时间为0.110 s;全局扫描算法约在0.830 s稳定,扫描时间为0.330 s; MPT算法则约在0.840 s稳定,扫描时间为0.340 s. 总路径的扫描电压是各算法通过电压步长进行扰动的电压区间总和. 其中,改进型快速GMPPT算法包含3个扫描区间,分别为从0.9Voc,STC,array (88.83 V)到最右侧峰值点(84.80 V)、从第一次跨越恒流区后的电压(61.77 V)到中间的峰值点(55.84 V)、从第二次跨越恒流区后的电压(32.81 V)到最左侧峰值点(25.43 V),总路径的扫描电压约为17.34 V. 全局扫描算法的扫描区间为[0.6Voc,STC,module, 0.9Voc,STC,array],总路径的扫描电压为69.09 V. MPT算法的扫描区间为[PGMPP/IGMPP,STC,array, 0.9Voc,STC,array],总路径的扫描电压为71.14 V. 综上,与全局扫描算法相比,所提算法在扫描时间、总路径的扫描电压、能量损失方面分别减少66.67%、74.86%和58.19%;与MPT算法相比,这些指标则分别减少67.65%、75.63%和62.31%.
当光照条件由PSC2改变为PSC3时,与全局扫描算法相比,所提算法在扫描时间、总路径的扫描电压、能量损失方面分别减少67.80%、74.57%和66.13%;与MPT算法相比,这些指标则分别减少68.33%、74.89%和66.14%.
当光照条件由PSC3突变为PSC1时,所提算法与全局扫描算法相比,扫描时间、总路径扫描电压、能量损失分别减少60.94%、74.57%和72.96%;与MPT算法相比,这些指标则分别减少64.79%、77.48%和76.91%.
3.2 实验结果
为验证改进型快速GMPPT算法的理论分析和仿真结果,本文搭建了由光伏模拟器、Boost变换器、电子负载、上位机、GMPPT数字控制器和示波器组成的实验平台,如图10所示. 其中,上位机控制光伏模拟器Chroma 63206A-600-420模拟表1的局部阴影条件. 图11为3种算法的实验波形,实验中各项参数保持一致,电压步长为1 V,扰动时间为0.4 s,其余硬件参数与仿真设置相同.
由图11可知:当光照条件由PSC1改变为PSC2时,相较于全局扫描算法,所提算法的扫描时间、总路径扫描电压、能量损失分别减少68.25%、74.86%和58.19%;与MPT算法相比,分别减少68.00%、75.63%和62.31%. 当光照条件由PSC2改变为PSC3时,所提算法相比全局扫描算法的上述3项指标分别减少65.85%、74.57%和66.13%;与MPT算法相比,分别减少66.13%、74.89%和66.14%. 当光照条件由PSC3改变为PSC1时,所提算法与全局扫描算法相比,3项指标分别减少64.29%、74.57%和72.96%;与MPT算法相比,分别减少67.39%、77.48%和76.91%.
表3为3种算法的实验跟踪性能对比. 可以看出,这3种算法均能准确实现GMPPT;实验结果与仿真结果一致,进一步验证改进型快速GMPPT算法的有效性.
表 3 阴影变化下3种算法的实验跟踪性能比较Table 3. Experiment tracking performance comparison of the three algorithms under shading variation光照条件 扫描时间/s 总路径的扫描电压/V 能量损失/J 算法 1 算法 2 算法 3 算法 1 算法 2 算法 3 算法 1 算法 2 算法 3 PSC1→PSC2 4.0 12.6 12.5 17.34 69.09 71.14 98.19 234.86 260.49 PSC2→PSC3 4.2 12.3 12.4 17.57 69.09 69.96 158.23 467.10 467.35 PSC3→PSC1 4.5 12.6 13.8 17.34 69.09 76.99 72.00 266.25 311.85 4. 结 论
本文提出了一种适用于局部阴影条件下的改进型快速GMPPT算法,分析局部阴影条件下光伏阵列的输出特性,讨论了膝点与开路电压的关系;详细阐述了现有MPT算法与改进型快速GMPPT算法的工作原理,对比全局扫描算法、MPT算法以及本文算法在扫描时间、电压扫描路径与能量损失方面的表现. 通过理论分析、仿真与实验的验证,可得以下结论:
1) 改进型快速GMPPT算法无须增加额外硬件电路,节省了成本.
2) 与软件方法中的智能算法相比,本文算法更为简单,实现方便.
3) 相比于现有MPT算法和全局扫描算法,本文算法能够更快、更准确地跟踪到GMPP,缩短电压扫描路径,减少扫描时间,有效提高跟踪速度,降低能量损失,同时不增加算法复杂度,保持了算法的精确度和稳定性.
-
表 1 直供加回流供电方式牵引网导线规格参数
Table 1. Specifications of traction network conductors for direct power supply with return wire
导线类型 导线规格 计算半
径/mm电阻/
(Ω·km−1)等效半
径/mm接触线 CTA-120 6.45 0.148 5.03 承力索 JTMH-95 5.90 0.303 5.30 回流线 LBGLJ-185 7.85 0.146 6.11 表 2 钢轨规格参数
Table 2. Specifications of rail
参数 规格/(kg·m−1) 截面面积/cm2 计算半径/mm 电阻/(Ω·km−1) 等效半径/mm 取值 60 77.08 109.10 0.135 12.79 表 3 埋地管道参数
Table 3. Parameters of buried pipeline
参数 规格 管道金属相对磁导率 300 管道金属电阻率/(× 10−7 Ω·m) 1.7 管道涂层相对介电常数 5 管道涂层电阻率/(× 105 Ω·m2) 1.0 管道半径/m 0.3 管道涂层厚度/m 0.003 管道埋深/m 1.5 表 4 负荷电流变化时埋地管道感应耦合电压最大值
Table 4. Maximum values of inductive coupling voltage of buried pipeline during load current variations
负荷电流幅值/A 数学模型/V 仿真模型/V 偏差/% 200 32.5 30.7 5.5 400 64.9 61.5 5.2 600 97.1 92.2 5.0 800 129.5 123.2 4.9 1000 162.2 154.4 4.8 表 5 间距比值变化时埋地管道感应耦合电压最大值
Table 5. Maximum values of inductive coupling voltage of buried pipeline with varying distance ratios
Ri 数学模型/V 仿真模型/V 偏差/% 1.5 75.9 73.5 3.2 2.0 72.2 70.2 2.8 2.5 69.7 67.8 2.7 3.0 67.3 65.2 3.1 3.5 65.6 63.2 3.7 4.0 63.8 61.1 4.2 4.5 62.4 59.4 4.8 表 6 某动车组谐波参数
Table 6. Harmonic parameters of electric multiple unit
主要谐波次数/次 含有率/% 3 1.38 5 2.84 21 3.00 23 4.07 25 2.40 -
[1] 曹建猷. 电气化铁道供电系统[M]. 北京:中国铁道出版社,1983. [2] 李群湛,贺建闽. 牵引供电系统分析[M]. 3版. 成都:西南交通大学出版社,2012. [3] 国家发展改革委,交通运输部,国家铁路局,等. 关于推动都市圈市域(郊)铁路加快发展的意见[Z]. 2020. [4] 万红霞,李婷婷,宋东东,等. 杂散电流对埋地管道的腐蚀及排流方式的研究进展[J]. 表面技术,2021,50(4): 125-134.WAN Hongxia, LI Tingting, SONG Dongdong, et al. Research progress of stray current on corrosion and drainage method of buried pipeline[J]. Surface Technology, 2021, 50(4): 125-134. [5] 梁毅,杜艳霞. 交流干扰和阴极保护协同作用下的腐蚀评判标准与机理研究进展[J]. 中国腐蚀与防护学报,2020,40(3): 215-222. doi: 10.11902/1005.4537.2019.053LIANG Yi, DU Yanxia. Research progress on evaluation criteria and mechanism of corrosion under cathodic protection and AC interference[J]. Journal of Chinese Society for Corrosion and Protection, 2020, 40(3): 215-222. doi: 10.11902/1005.4537.2019.053 [6] 朱久国. 交流电气化铁路对埋地油气管道电磁干扰特性研究[D]. 成都: 西南交通大学,2018. [7] 胡士信,路民旭,杜艳霞,等. 管道交流腐蚀的新观点[J]. 腐蚀与防护,2010,31(6): 419-424. [8] 吴荫顺,曹备. 阴极保护和阳极保护:原理、技术及工程应用[M]. 北京: 中国石化出版社,2007: 223. [9] WAKELIN R G, SHELDON C. Investigation and mitigation of AC corrosion on a 300 mm diameter natural gas pipeline[C]//The 59th NACE Annual Conference. New Orleans: Corrosion, 2004: 04205.1-04205.9. [10] FLOYD R. Testing and mitigation of AC corrosion on 8 line: a field study[C]//The 59th NACE Annual Conference. New Orleans: Corrosion, 2004: 04210.1-04210.9. [11] LINHARDT P, BALL G. AC corrosion: results from laboratory investigations and from a failure analysis[C]//The 61th NACE Annual Conference. Houston: Corrosion, 2006: 06160.1-06160.9. [12] 任增珺. 日东管道杂散电流干扰检测与防护[J]. 油气储运,2015,34(1): 111-114. doi: 10.6047/j.issn.1000-8241.2015.01.024REN Zengjun. Detection and protection of stray current interference of Rizhao–Dongming oil pipeline[J]. Oil & Gas Storage and Transportation, 2015, 34(1): 111-114. doi: 10.6047/j.issn.1000-8241.2015.01.024 [13] CARSON J R. Wave propagation in overhead wires with ground return[J]. Bell System Technical Journal, 1926, 5(4): 539-554 doi: 10.1002/j.1538-7305.1926.tb00122.x [14] 高攸纲,沈远茂,石丹. 交流电气化铁道对周围电气及电子系统的阻性耦合影响[J]. 邮电设计技术,2007(3): 57-60. doi: 10.3969/j.issn.1007-3043.2007.03.013GAO Youzang, SHEN Yuanmao, SHI Dan. Resistive coupling effects of AC electric railway to surrounding electric and electronic system[J]. Designing Techniques of Posts and Telecommunications, 2007(3): 57-60. doi: 10.3969/j.issn.1007-3043.2007.03.013 [15] 蒋俊. 交流线路正常运行时对平行敷设油气管道的电磁影响[J]. 电网技术,2008,32(2): 78-80,92.JIANG Jun. Electromagnetic influence of normally operating AC power transmission line on gas/oil pipeline parallel to transmission line[J]. Power System Technology, 2008, 32(2): 78-80,92. [16] 汪可. 电气化铁路对油气管道的影响及防护措施[D]. 成都: 西南交通大学,2013. [17] YONG J, XIA B, YONG H, et al. Harmonic voltage induction on pipelines: measurement results and methods of assessment[C]//IEEE Transactions on Power Delivery. [S.l.]: IEEE, 2018: 2170-2179. [18] 陈民武,朱久国,解绍锋,等. 牵引供电系统对埋地管道阻性耦合交流干扰建模及仿真[J]. 中国铁道科学,2018,39(2): 80-86. doi: 10.3969/j.issn.1001-4632.2018.02.10CHEN Minwu, ZHU Jiuguo, XIE Shaofeng, et al. Modeling and simulation of resistive coupling AC interference of traction power supply system to buried pipeline[J]. China Railway Science, 2018, 39(2): 80-86. doi: 10.3969/j.issn.1001-4632.2018.02.10 [19] 盛望群. 基于CDEGS的交流电气化铁路对沿线油气管道电磁干扰影响研究[J]. 铁道科学与工程学报,2020,17(8): 2101-2108.SHENG Wangqun. Study on influence of AC electrified railway on electromagnetic interference of oil and gas pipelines along the line based on CDEGS[J]. Journal of Railway Science and Engineering, 2020, 17(8): 2101-2108. [20] CHARALAMBOUS C A, DEMETRIOU A, LAZARI A L, et al. Effects of electromagnetic interference on underground pipelines caused by the operation of high voltage AC traction systems: the impact of harmonics[J]. IEEE Transactions on Power Delivery, 2018, 33(6): 2664-2672. doi: 10.1109/TPWRD.2018.2803080 [21] 齐磊,崔翔,郭剑,等. 特高压交流输电线路正常运行时对输油输气管道的感性耦合计算模型[J]. 中国电机工程学报,2010,30(21): 121-126.QI Lei, CUI Xiang, GUO Jian, et al. Inductive coupling modelling of normally operating UHV AC transmission line to adjacent oil/gas pipeline[J]. Proceedings of the CSEE, 2010, 30(21): 121-126. [22] 国家质量监督检验检疫总局. 接地系统的土壤电阻率、接地阻抗和地面电位测量导则 第1部分:常规测量:GB/T 17949.1—2000[S]. 北京: 中国标准出版社,2000. [23] 解绍锋,孙镜堤,骆冰祥,等. 高速铁路对邻近普速铁路电力电缆的干扰机理[J]. 西南交通大学学报,2021,56(1): 206-213.XIE Shaofeng, SUN Jingdi, LUO Bingxiang, et al. Mechanism of high-speed railway interference on power cables of adjacent normal-speed railway[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 206-213. [24] 周胜军,谈萌. 基于监测数据的高铁动车组谐波特性分析[J]. 电力科学与技术学报,2018,33(3): 128-133. doi: 10.3969/j.issn.1673-9140.2018.03.020ZHOU Shengjun, TAN Meng. Harmonic characteristics analysis of electric multiple units in high-speed railway based on the monitoring data[J]. Journal of Electric Power Science and Technology, 2018, 33(3): 128-133. doi: 10.3969/j.issn.1673-9140.2018.03.020 [25] 周秀荣. 电气化铁路对通信线路的干扰影响及防护措施[J]. 中国铁路,2007(6): 66-68. doi: 10.3969/j.issn.1001-683X.2007.06.017 -