天然气管道放空时间的计算(精).docVIP

xi ? xL xL ? xi ?1 ? ? pL ? ?x p ? i ? 1, j ? ? ?x p ? i, j ? ? xi ? xL x ?x ? ? ? i ? 1, j ? ? L i ?1 ? ? i, j ? ? ?L ? ?x ?x ? xi ? xL xL ? xi ?1 ? ? TL ? ?x T ? i ? 1, j ? ? ?x T ? i, j ? ? 3.2.3 显式差分方程组的建立 将全微分方程组改写为如下形式： (22 dp ? B1 d? ? ? 其中： z d? ? C1 d? y z (23 dp ? B2 d? dp ? B3 d? d? ? C2 d? d? ? C3 d? y (24 C C (25 B1 ? as ? C1 ? as2 ? T ? ?Z ? ? ? ? ?? 3 ? 1 ??? 3 as ?1 ? ? ?? ? ?? c pT ? ? Z ? ?T ? p ? ?? 2 d ? 2 d （26） （27） B2 ? as ? C2 ? ? as2 c pT ? T ? ?Z ? ? ? ? ?? 3 ? 1 ??? 3 as ?1 ? ? ?? ? ?? 2 2 Z T d d ? ? ? ? ? p ? ? ? ? （28） （29） ? T ? ?Z ? ? B3 ? ? c p ?1 ? ? ? ? ? Z ? ?T ? p ? C3 ? ? as2 pc p ? T ? ?Z ? ? ? ? ?? 3 ? ? p ? ?Z ? ? ?1 ? ? ? ?1 ? ? ? ?? ? ? 2 Z T d ? ? ? ? ? p ? ?? ? ? ? Z ? ?T ? p ? ? （30） （31） 将全微分方程组（32）?（34）离散成差分方程组： p ? i, j ? 1? ? pL ? ? i, j ? 1? ? ? L ? B1 z ? C1 ?? ?? ? p ? i, j ? 1? ? pR ? B2 ?? y z (32 ? ? i, j ? 1? ? ? R ? C2 ?? y (33 6 ? p ? i, j ? 1? ? pC ? B3 ?? C T ? i, j ? 1? ? TC ? C3 ?? C (34 式（32）?（34）共有 3 个未知量： p ? i, j ? 1? 、 ? ? i, j ? 1? 、 T ? i, j ? 1? ，解方程组可以得出 3 个未知量 的值。 3.2.4 边界条件的处理 对计算模型的左右边界进行处理，可得模型的左边界条件，随着放空的不断进行，干线压力随时间 不断下降，i+1 时刻的压力值为 i 时刻压力值减去由于 ?? 时段的放空而产生的压降。而对于干线的温度 值，则认为在整个放空过程中与环境温度相同。天然气从较大管径的干线流入较小管径的放空管，可视 为等熵绝热膨胀，放空管入口的压力、温度可由下式计算： p1 ? p0 1 ?1 ? 0.5 ?? ? 1? Ma ? ? ?1 ? ? （35） 2 ? ?1 1 T1 ? p1 ? ?? ? T0 ? p0 ? （36） 式中， p1 、 T1 为放空管入口处的压力、温度值， p0 、 T0 为干线的压力、温度值， ? 为气体的绝热指 数，Ma 为放空管入口处气体的马赫数。将右边界条件简化计算可得放空火炬出口处的压力为大气压。而 对于出口处温度，也做定值处理，等于环境温度。 3.3 实例计算 选取西气东输二线干线某段输气管道的实际放空数据，对所建计算模型进行对比验证。该段干线管 径 1 219×18.4 mm，管长 10.02 km，放空前管内压力 4.5 MPa，放空结束时管内压力 0.05 MPa（均为表 压）。放空管管径 406×12.5 mm，放空管长 450 m。将特征线法计算的放空参数和实际数据进行对比。 使用显示特征差分方法计算所得放空管入口压力随放空时间的变化曲线以及气体平均流速随时间变化 的曲线与实际管道放空曲线基本吻合（图 4，图 5）。实际放空时间为 16.5 h，计算所得放空时间约为 15 h，计算误差为 9.6%，在可接受范围内，说明该计算模型可以较为准确地预测管道的放空时间。 在放空前期，管道压力高，气体流速较大，可达 270 m/s，接近声速，放空速率大，管道内压降速率 也很大。随着放空的进行，气体流速随着压力的快速下降大幅减小，管道内的压降速率随之减小（图 4， 图 5）。在放空后期，管道内压力接近大气压时，气体速率缓慢，仅 15 m/s 左右。此次放空由 4.5 MPa 放 空至 2.0 MPa 仅用约 3.5 h，而从 0.5 MPa 放空至 0.05 MPa 则用时 6.5 h（图 4）。低压力下放空时间过 长