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阀关断动作的实际执行时间。

　　2.不同关断条件下的泄漏速率变化。

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

　　3.不同关断条件下的管道流量变化。

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

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