1自然气管道泄漏速率

　　自然气扩散是一种非常复杂的大气污染物扩散过程，在对其模拟评价过程中存在大量不确定性因素，这些不确定性因素直接影响了模拟结果的准确性，泄漏速率或泄漏量即为其中之一。

　　自然气泄漏速率是模拟评价自然气管线泄漏事故过程中的重要因素之一，对其研究将有利于进步模拟评价结果的准确性。应用热力学和气体动力学理论，结公道想气体状态方程、绝热方程和能量守恒方程，研究分析了自然气管道的泄漏过程，分别给出了在临界泄漏阶段以及亚临界泄漏阶段的泄漏速率计算公式，为自然气管输安全工程的相关计算提供了可靠的理论依据。

　　2基于HYSYS的动态自然气集输管道泄漏模型

　　自然气管道泄漏模型与合用性分析：

　　1）稳态泄漏模型。这类模型*常见的是一般泄漏速率模型，利用气体在泄漏点活动状态有关的特性，以为泄漏发生时泄漏处管内压力不变，通常根据各种守恒方程可得出解析解的泄漏速率方程。该方法合用于泄漏时管道内气体压力恒定工况。也有良多研究将小孔泄漏近似为稳态泄漏进行计算。

　　2）动态非不乱泄漏模型。目前动态泄漏模型均主要对大孔泄漏和管道泄漏情景进行计算和研究，因难以确定实际泄漏量，基本都假设管道发生泄漏后ESD阀同时封闭，以为泄漏量即为两ESD阀间气体量，并根据能量守恒方程等得到泄漏速率方程。实际上，自然气管道泄漏后，因为受泄漏检测技术、职员干涉干与、执行机构响应等因素的制约，泄漏处上下游ESD阀关断动作并不能保证同步，即存在不同程度的关断延迟。在此鉴戒加拿大阿尔伯特能源与公共事业局（AEUB）推荐的自然气管道有效泄漏长度概念，建立管道泄漏的等效启齿系统模型，结合过程模拟技术，研究ESD阀关断延迟对泄漏过程的影响。

　　3ESD阀关断动作对泄漏率的影响

　　3.1计算算例前提。

　　因为高压自然气管道泄漏速度较高，为清楚对比泄漏后的动态平衡过程和ESD阀关断延迟影响，基于前述理论和模型，建立20000m长、内径为254mm水平自然气管道HYSYS泄漏等效模型。计算管道泄漏常见的小孔、大孔范围内的泄漏口面积，小孔取内管道内横截面积的1%（50615mm2），大孔取10%（5065mm2）。对于全尺寸断裂，泄漏视为管道泄漏模式，泄漏速率为上下游管道单独泄漏速率之和，2个ESD阀延迟相互影响很小，篇幅所限，这里不再分析。为求不同泄漏孔径的气体泄漏速度系数，将等效模型内介质指定为纯甲烷并与美国环保署和国家海洋大气治理局推荐的ALOHA计算数据进行对比验证，固然ALOHA只能对纯物质提供ESD阀封闭后泄漏量的分钟内均值，但是可以近似得到不同泄漏孔径的泄漏系数，进步模型计算泄漏速率正确度。在此小孔和大孔速度系数分别取0165和0185.设定的前提为：1）因管道处于集输管网的背景压力环境，指定管段的进口压力（本文取9512kPa），下游出口压力（稳态模式求解后指定）；2）进口吻体温度为45℃；3）指定实际的气体成分，相关参数如流经上下游ESD阀的流量，泄漏点处压力、温度等由HYSYS计算。基于前述构建的泄漏等效模型，在验证其有效性后，本文利用HYSYS的过程计算上风，计算管道摩擦压力损失变化和沿程热量损失，将管道气流改为实际自然气成分（本文为相对分子质量为19141的高含硫自然气），分析管道泄漏段上下游ESD阀封闭时机对泄漏的影响。为分析和对比ESD阀关断顺序对小孔及大孔泄漏速率的影响，定义以下动作时间点：1）50min为管道由正常平稳运行状态转为泄漏事故状态的开始时刻；2）150min时上游（或下游）ESD阀立刻执行关断；3）250min时剩余的下游（或上游）ESD阀执行关断。

　　需说明的是，所设各动作时间点只为清楚对比ESD阀关断影响而设定，并不代表管道ESD阀关断动作的实际执行时间。

　　3.2不同关断前提下的管道流量变化。

　　泄漏的发生使管道出进口流量不再相等，为计算有效泄漏质量，在分析泄漏事故状态时，了解自然气管道流量变化有助于分析泄漏速率的变化影响关系。对于小孔泄漏，泄漏发生后（50～150min）泄漏段管道进口流量增加，出口流量降低，泄漏量由泄漏段的上游管道所提供，即在泄漏段管道出口无回流。

　　当上游ESD阀先关断后，下游管道出口将立刻泛起逆流，增补泄漏量。当下游ESD阀先行封闭时，上游管道进口处流量降低，但持续增补流量显著大于上游先关断段时的逆向增补流量，将造成较多的气体开释。当泄漏为大孔泄漏时，泄漏量则由ESD阀上游管道和ESD阀出口的下游管道共同提供，上游进口流量明显增加，下游管道出口泛起大的逆向回流，造成更多的气体开释。当上游ESD阀先行封闭后，泄漏管段下游出口的逆向回流量大幅增加，增补泄漏量。当下游ESD阀先行封闭后，管内上游流量大幅增加，但比拟低于上游ESD阀先关断后的出口逆向增补流量。不同孔径和ESD阀不同关断顺序下管内流量存在差异，流量变化存在动态平衡过程，这些因素都决定实际的气体开释质量和速率。显然，如不能保证ESD阀同时关断，应优先关断管道进口ESD阀，此时不管大小孔增补流量都相对较低。

　　3.3不同关断前提下的泄漏速率变化。

　　泄漏段管道出进口流量在不同关断顺序和时差前提下的变化决定着泄漏率的变化，泄漏发生后（50min时）因为管内压力存在动态平衡过程，小孔因为泄漏速率低，平衡时间很短，初始时刻流量稍高，十几秒内达到压力平衡后的稳态泄漏流量（由6001215kg.min-1降至585159kg.min-1）。而大孔泄漏后初始时刻流量较高，需经由10余分钟系统方达到压力平衡后的稳态泄漏，降幅较高（由5565121kg.min-1降至4377133kg.min-1）。实际过程中当上游ESD阀先关断时，大孔和小孔泄漏速率都会降低，但是大孔泄漏速率降低程度显著高于小孔泄漏，泄漏速率在上游ESD阀关断达到不乱后泄漏速率降低34188%小孔泄漏速率仅仅降低3145%.当下游ESD阀先行封闭时，大孔泄漏速率也会大幅下降27153%而小孔泄漏速率并未降低反而有所增加。

　　4结论

　　利用本文方法可对不同间隔和梯度变化的管道，定义不同的泄漏位置和尺寸，结合管道过程控制参数建立泄漏速率计算模型，为自然气管道泄漏控制提供有益参考。通过本文研究可得到以下结论。

　　1）跟着泄漏孔径的增加，管道内将泛起回流现象，管道ESD阀关断延迟顺序对小孔泄漏和大孔泄漏速率的控制效果存在不同影响。先关断上游ESD阀时，不管小孔仍是大孔泄漏，都使下游出口ESD阀处泛起逆向流，增补泄漏。而下游先关断时，小孔泄漏使上游进口处管道流量有所降低，而大孔泄漏却使进口处流量大幅增加。2）ESD阀的关断延迟和关断顺序，主要表现在对泄漏速率的影响，即影响单位时间内开释到大气的自然气量，而相同距离不同顺序的关断对开释持续时间的影响较小。3）在现有自然气管道关断多采用压力、流量监控等自动关断的现状下，为*大限度地控制事故自然气泄漏速率和开释时间，管道的紧急关断应优先采用联动的同时关断方式，其次为上游先关断方式。