华东成品油管道泄漏量及压降的计算(导师:董荣国编写:张一飞杨荣丁孙韬吴进)摘要:在十二五规划期间,我国又将建立多条成品油管道,这给我国成品油管道的发展带来了契机,同时,对成品油管道的安全运行也提出了更高的要求。特别是“11.22”青岛东黄复线输油管线泄漏爆炸事故的发生,引发了人们对于成品油管道泄漏这一安全问题的高度重视和深度思考。华东分公司,做为华东区最大的成品油物流中心,现有三条在运成品油管线(浙苏管线、苏南管线、苏北管道),华东成品油管网已初具规模,为了更好的应对管道泄漏事故,保护国有资产安全,迎合公司倡导的“平安管道”建设,本文主要对成品油管线泄漏量,以及泄漏所引起的漏点压降这一课题展开深入的研究。关键词:压力波漏点压降泄漏量1、华东管网压力波传播速度的分析1.1实际波速采集在实际生产过程中,通过各种工况发生变化或调节阀调节下载流量而产生相应管段的压力波,利用泄漏检测软件对管线各站进出站压力数据的高精度采集,查找对应的压力特征突变点,来测算总结压力波传递速度。1.1.1浙苏管线1)柴油状态下:表1.1序号陈山站嘉兴站间距(km)时间差(s)实际波速(km/s)19:07:019:07:4250.591411.2229:09:549:10:3550.591411.2239:27:499:28:3150.591421.2049:32:219:34:0350.591421.20序号嘉兴站苏州站间距(km)时间差(s)实际波速(km/s)112:58:5012:59:4770.245571.23213:00:1813:01:1570.245571.23313:02:0713:03:0470.245571.23413:04:0313:05:0070.245571.23序号嘉兴站桐乡站间距(km)时间差(s)实际波速(km/s)14:13:194:12:3751.613421.2324:15:024:14:2051.613421.2334:17:294:16:4751.613421.2344:30:084:29:2651.613421.23序号桐乡站湖州站间距(km)时间差(s)实际波速(km/s)14:14:194:14:5949.36401.2324:16:464:17:2649.36401.2334:29:254:30:0549.36401.2344:33:194:33:5949.36401.232)汽油状态下:表1.2序号陈山站嘉兴站间距(km)时间差(s)实际波速(km/s)120:12:2020:13:0750.591471.07220:12:4520:13:3250.591471.07320:13:4820:14:3550.591471.07420:13:3220:14:1850.591471.07序号嘉兴站苏州站间距(km)时间差(s)实际波速(km/s)14:59:264:58:2170.245651.0825:00:545:01:5970.245651.0835:02:435:03:4870.245651.0845:04:395:05:4470.245651.08序号嘉兴站桐乡站间距(km)时间差(s)实际波速(km/s)123:26:4423:25:5651.613481.08223:26:5823:26:1051.613481.08323:30:4723:29:5951.613481.08423:31:1423:30:2651.613481.08序号桐乡站湖州站间距(km)时间差(s)实际波速(km/s)121:20:2321:21:0949.36451.09221:22:0321:22:4949.36461.07321:49:0521:49:5149.36461.07422:00:0122:00:4749.36451.091.1.2苏南管线1)柴油状态下表1.3序号南京站镇江站间距(km)时间差(s)实际波速(km/s)19:37:229:38:3183.54691.2129:41:509:42:5983.54691.2139:44:399:45:4883.54691.2149:53:359:54:4483.54691.21序号镇江站常州站间距(km)时间差(s)实际波速(km/s)122:32:2322:33:1562.59521.20222:38:5122:39:4362.59521.20322:41:4022:42:3162.59511.22422:51:3622:52:2762.59511.22序号常州站江阴站间距(km)时间差(s)实际波速(km/s)117:42:0617:43:1078.29641.22217:48:3417:49:3978.29651.20317:51:2317:52:2878.29651.20418:01:1918:02:2378.29641.22序号江阴站无锡站间距(km)时间差(s)实际波速(km/s)117:19:3517:20:2865.17531.23217:23:0317:23:5665.17531.23317:30:5217:31:4565.17531.23417:36:4817:37:4165.17531.23序号无锡站苏州站间距(km)时间差(s)实际波速(km/s)12:42:072:43:1077.18631.2322:46:352:47:3877.18631.2332:50:242:51:2777.18631.2342:59:203:00:2377.18631.232)汽油状态下表1.4序号南京站镇江站间距(km)时间差(s)实际波速(km/s)116:07:2316:08:4183.54781.07216:21:5116:23:0983.54781.07316:44:4016:45:5883.54781.07416:53:3616:54:5483.54781.07序号镇江站常州站间距(km)时间差(s)实际波速(km/s)115:18:4515:19:4462.59591.06215:33:1315:34:1162.59581.08315:56:0215:57:0162.59591.06416:04:5816:05:5662.59581.08序号常州站江阴站间距(km)时间差(s)实际波速(km/s)17:55:297:56:4278.29731.0728:09:578:11:1078.29731.0738:32:468:33:5978.29731.0748:41:428:42:5578.29731.07序号江阴站无锡站间距(km)时间差(s)实际波速(km/s)17:47:387:48:3865.17601.0827:56:067:57:0665.17601.0838:16:158:17:1565.17601.0848:21:088:22:0865.17601.08序号无锡站苏州站间距(km)时间差(s)实际波速(km/s)17:20:457:21:5677.18711.0927:26:137:27:2477.18711.0937:41:227:42:3477.18721.0747:56:157:57:2777.18721.07根据华东管网浙苏和苏南管道的实际数据,代入压力波速理论公式11kakDceE进行计算,可以得到浙苏管线与苏南管线中不同管径管道的理论波速数据。式中:a—瞬变压力波速,m/s;k—液体体积弹性系数,Pa;ρ—液体密度,kg/m3;E—管材弹性模量,Pa;D—管径,m;e—管壁厚度,m;c1—管道约束系数;g—重力加速度,m/s2。406mm管径管道中输送汽油时的压力波速度:20℃时,汽油K=9.16×108Pa,ρ=740kg/m³,D=(406.6-7.1-7.1)/1000=0.3922m,E(钢)=206.9×109Pa,e=0.0071m,c1=1-0.32。波速a约为1006.23m/s。输送柴油时的压力波速度:20℃时,柴油K=12×108Pa,ρ=850kg/m³,D=(406.6-7.1-7.1)/1000=0.3922m,E(钢)=206.9×109Pa,e=0.0071m,c1=1-0.32。波速a约为1045.50m/s。323mm管径管道中输送汽油时的压力波速度:20℃时,汽油K=9.16×108Pa,ρ=740kg/m³,D=(323.9-6.4-6.4)/1000=0.3111m,E(钢)=206.9×109Pa,e=0.0064m,c1=1-0.32。波速a约为1017.41m/s。输送柴油时的压力波速度:20℃时,柴油K=12×108Pa,ρ=850kg/m³,D=(323.9-6.4-6.4)/1000=0.3111m,E(钢)=206.9×109Pa,e=0.0064m,c1=1-0.32。波速a约为1059.96m/s。说明:以上计算结果仅为根据理论公式算出的理论值,华东管网陈山至嘉兴段、南京至江阴均为406mm管径的管道,嘉兴至湖州、嘉兴至苏州,江阴至苏州则为323mm管径的管道,由于温度、油品密度、体积弹性系数等均为变化数据,因而这里选取的数据均为标准状态下理论数据值。在实际生产过程中,通过各种工况发生变化或调节阀调节下载流量而产生相应管段的压力波,利用泄漏检测软件对管线各站进出站压力数据的高精度采集,查找对应的压力特征突变点,来测算总结压力波传递速度。经分析浙苏管线与苏南管线的波速基本相同,区别在与管径大小及管道内输送的油品介质,归纳总结出华东管网管段波速数据:表1.5油品种类406管径波速(km/s)323管径波速(km/s)汽油1.071.08柴油1.211.23将实际数据与理论计算数据进行比较,发现压力波的传播速度存在一定的差距,其主要原因是由于理论计算中所取的汽柴油密度与实际管道中油品的密度并不一致,两者密度之间的偏差是造成实际数据和理论计算结果有所差异的原因之一,同时温度对于液体的体积弹性系数也有影响,温度越高,液体的体积弹性系数减小,液体的密度减小,液体的可压缩性增大,压力波的传播速度减小,反之温度越低,压力波的传播速度则增大。理论计算中所取的为标准环境20℃的油品体积弹性系数等相关系数,因而也是造成差异的原因之一。同时压力仪表的精确度及时间计算的时间精度也是造成差异的原因之一。但是通过大量的数据采集与整理,可以证明理论公式的普遍理论,即406mm管径管道的压力波速在相同介质的情况下略小于323mm管径管道的压力波速;在相同管径的管道中,柴油介质的压力波速较为明显的大于汽油介质的压力波速。由于符合理论计算的对比结果,证明了华东管网压力波速经验值的真实可用性。利用大量的数据统计,可确认华东管网的压力波速如表1.5所示,并可用于日后实际生产中的相关计算。2、华东管网压力波衰减率的分析2.1压力波衰减数据采集方式针对压力波衰减的程度,目前在专业领域也还未有较为详细明确的理论公式,主要原因为影响压力波衰减的因素较多,相互影响的方式也较为复杂,且由于长输成品油管道无法进行高密度多频次的检测,同时也无法进行精确的多向实验,因而无法确定压力波在长输成品油管道中衰减程度的具体情况。因而我们只能选取较为有代表性的管段进行大量的数据采集比对,来逐步判别影响压力波衰减的各项主要因素,并得出相应管段的压力波衰减率经验值。利用华东管网上线的泄漏检测软件进行相应的测试数据采集,因为泄漏检测软件在进出站压力的数据采集上精度达到了0.0001MPa,高于SCADA系统采集的数据精度。同样由于条件限制,只能采集到沿线各站的进出站压力,无法采集管道中压力的变化。因此我们将调节某站调节阀开大下载来作为压力波发生的初始点,该站的进出站压力突变值视为压力波产生的自变初始值;在该站的上游站场或下游站场因该压力波传递站场而产生的进出站压力突变值视为压力波的因变值。通过自变初始值与因变值的比对,并结合站场之间的距离得出该管段的压力波衰减率。即:自变初始值=调节前压力值-调节后压力值因变值=(上