高压天然气管网动态模拟与壅塞流动特性研究

【摘要】：天然气作为优质清洁能源得到世界各国的高度重视,随着我国粗放型经济发展方式的转变,能源政策相应调整,天然气在能源结构中的比重不断增加,输送天然气的长输管线和城市高压管网数量越来越大,储存天然气的设施规模不断增大。对高压天然气管道的科学计算和设计,以及在实际运行中正确的判断和调节,直接决定着高压天然气管网系统的安全稳定运行和保障城市正常供气。在设计过程中,由于缺乏计算高压天然气的有效方法和工具,设计人员在水力计算过程中,往往采用稳定流动状态代替不稳定流动状态进行计算,所设计的管网压力、流量等参数与实际运行数据存在较大偏差。在实际生产运行中,缺乏与实际管网相匹配的模拟工具,不能准确评估出输气管线在各种情况下的工况,不能准确控制管线的安全稳定运行。因此有必要针对我国高压天然气管网的实际情况,深入研究高压天然气管网中的流动状况及影响因素,通过动态模拟找出运行参数的变化规律,并提出相应的控制策略。 鉴于此,本文采用建立高压天然气管网数学模型、数值求解、动态仿真模拟、实例对比分析等方法,对高压天然气管网的水力、热力动态特性、高压储气、以及壅塞流动等方面进行了研究。主要结论如下： (1)对不稳定流数学模型的建立作了大量简化：①天然气沿程温度按等温考虑；②忽略质量流量随时间的变化；③忽略水平输气过程中沿程线速度变化的影响；④忽略管道位置高差的影响；⑤分别采用恰尔内(И.А.Чарный)线性化法和最小二乘法线性化方法对数学模型进行线性化处理等。本文针对目前天然气管道不稳定流动的计算过程中忽略对流项、惯性项和中立项等问题,首次在没有忽略任何项的连续性方程、动量方程、能量方程、状态方程和焓方程所建立的数学模型的情况下,分别采用了特征线法和隐式中心有限差分法对天然气管网不稳定流进行了求数值解。采用数值解法对高压天然气管网不稳定流动进行研究,不用对模型方程组进行简化和复杂的数学变换,求解出来的结果与管道实际运行数据更为接近。求解时要将待求的时间层次上所有的未知量联立起来同时进行求解,求解的非线性方程组很庞大,方程组过大时求解需要的时间长。为了保证数值解的稳定性,时间步长、距离步长都不能够取得过小。数值解考虑了气体与外界环境的热交换,考虑了天然气沿程温度变化,能够计算出天然气沿线的温度分布。由于供气与用气的不均匀性,管网中天然气的压力、流量与温度也随时间波动。如果管网存在有分气点,则在分气点处气体的压力、温度、密度的变化最为剧烈。 (2)与隐式法相比,特征线法主要优点在于它们能对各种边界条件进行显式处理,且计算占用计算机内存较小,但其时间步长和距离步长的比值受到一定稳定条件的限制,这使得时间步长只能取得很小如果分析时间过长时,就需要消耗很长的计算时间。在采用特征线法对天然气管道进行数值求解时,利用逆步进法差分网格建立起特征差分方程,对常微分方程组结合初始条件和边界条件进行求数值解,由于边界节点的参数随时间的变化比较明显,因此该方法能提高计算精度,求解过程容易实现。本文还采用了适用于慢瞬变流的守恒型求解数学模型的隐式中心有限差分法,并利用Newton-Raphson法对形成的高阶非线性方程组进行了迭代计算。隐式中心有限差分法是隐式格式的有限差分方法,表现出了二阶精度高、计算速度快的优点。 (3)对天然气管网进行模拟时,可以单独采用特征线法或隐式有限差分法。本文综合考虑这两种方法的优点,采用综合法求解。结果表明：在模拟高压天然气管网动态运行工况时,在长管元件的边界点上采用特征线法计算,在其它情况下采用隐式中心有限差分法。根据本文所建模型及提出的解法,用Visual C++编程开发出的高压天然气管网动态模拟软件,利用实例计算验证,计算结果可以满足实际工程需要。 (4)采用SHBWR状态方程替代常用的对比状态方程,采用Newton迭代法求解,大幅度提高了天然气状态参数的计算速度与精度。由于高压天然气管道中天然气的流态一般在紊流区,因此在确定管道摩阻系数时,选择计算较为准确的柯列勃洛克公式,采用牛顿迭代方式可以计算出管道摩阻系数。 (5)对输气管道进行水力分析,从输气管道的天然气流动基本方程出发,结合实际气体状态方程与焓方程建立了动态模拟的数学模型；区别于以往的试算法,采用弦截法进行摩阻系数的求解,提高了计算精度。将摩阻系数公式进行变形改写,使其与管道流量进行关联,随着输气管道流量改变而改变,减少了由于摩阻系数的定值输入而带来的误差。 (6)对输气管道的温度变化进行热力学分析,研究了输气管道内天然气温度的变化规律。通过计算比较表明：单独求解输气管道温度与利用含有能量方程的稳定流模型求解出的管道温度基本吻合,因而单独求解管道温度是可行的,为多点接气不稳定流动数学模型模型中温度初始值的求解提供了有效的方法。 (7)多点接气和单气源高压天然气管网数学模型所依据的连续性方程、动量方程以及能量方程是相同的,区别主要在边界条件的处理上,其初始值的求解方法和所建立的方程个数不同。本文提出了多点接高压天然气管网动态模拟初始值计算方法,区别于单气源动态模型,可以从边界条件和约束条件入手对多点接高压天然气管网动态方程组封闭性进行分析。其封闭性分析来看,每增加一个气源点,就要补充两个关于该气源的方程,一个是流量方程或压力方程,另一个是温度方程。从而构成封闭的方程组,达到可以求解条件。 (8)通过对多点接高压天然气管网的计算,得到各个气源点、各个分气点以及各个管段在不同剖面上不同时刻的压力和流量的变化情况。可根据各点压力要求,对各气源压力进行有效控制,对流量进行合理的调度分配。利用SCADA系统采集的数据作为边界条件,对管网进行动态模拟计算,同时将计算结果与SCADA采集的节点流量数据对比分析,可以作为检测数据采集系统的准确性与可靠性 (9)对天然气流动中的壅塞现象进行研究,建立高压天然气管网壅塞流动的数学模型,分别采用四阶Rung-Kutta法和Newton-Raphson迭代法对模型进行求解,并通过工程实例进行对比验证。结果表明计算结果完全可用于天然气管网的设计、运行及调度管理。 (10)以气体动力学为基础,对径管和等径管天然气管道的壅塞现象进行分析研究,采用四阶Rung-Kutta法对等径管道壅塞模型进行求解。结合工程实例,讨论了壅塞状态下管道内气体参数的变化规律。结果再从压缩机室向CNG槽车内加气以及通过CNG槽车向储罐卸气过程中,开始充气阶段由于加(卸)气管两端的压差较大并且管径较小,天然气在低压端急剧膨胀易发生壅塞流动。产生壅塞流动的判定式为：根据实例计算当管道始末端压力比在0.32-0.41之间时,将发生流动壅塞。在工程设计时,管道流量不能超过壅塞流量,否则正常供气会受到影响。等径输气管道天然气流动发生壅塞时,其压力、温度等参数是沿程变化的；压力、温度沿程降低,马赫数及流速沿程增大,流速趋向于声速；各截面的当地声速与密度和温度有关,沿线降低。长输管道末端门站所连接的城市输气管网的管径不应太小,否则在用气高峰时刻在变径处会产生流动壅塞,导致某个时段内输入城市管网的天然气量小于城市实际用气需求量,造成供气不足,引发供气安全事故。 (11)分别对高压管道的稳态储气和非稳态储气进行研究,通过拟合管道内压力变化规律,推导出稳态储气量的计算公式,并分析了管长、压力等参数对稳态储气的影响；结合非稳态储气的定义与变径管道壅塞流动模型,对长输管道末段流动壅塞及其对管道储气能力的影响进行分析。结果表明：长输管线末段管道储气量计算可以按稳态,也可按非稳态来考虑。稳态储气计算结果偏于保守,计算结果与实际结果有误差。按稳态工况计算管段储气量的时候,通常根据起点最大压力、末端最小压力以及管道的实际长度来计算储气量,未考虑其最大储气量和最优长度。如果管段过长,由于流动阻力加大,不但不能多储气,反而会影响管道的储气。在利用长输管线末段储气的设计中应该注意这一特点。按照非稳态计算结果德更准确。 (12)根据本文提出的动态模型和解法,对第一版天然气管网动态模拟软件进行改进,使其具有多点接气高压天然气管网动态模拟计算功能,改进了历史数据的输入、查询、新数据计算与存储等功能,有利于实现负荷预测与实时模拟的结合。通过对多个实际或设计的高压天然气管网在多种不同运行工况下的计算,并将计算结果与国外同类模拟软件(PIPELINE、TGNET、SPS等)计算结果、城市管网SCADA系统采集的实际运行数据进行对比,结果表明,各个调压站计算压力的变化趋势与采集压力的变化趋势基本一致,与国外软件计算结果的相对偏差在6%以内,与SCADA系统采集的实际运行数据的相对误差在3%以内。改进型软件可对实际管网进行动态模拟计算,模拟结果准确反映出管道压力的变化规律,满足高压管网运行的要求,能够指导管网的智能化运行。 从以上结论来看,研究建立的高压天然气动态模拟数学模型和提出的解法计算结果精度高、计算速度快。揭示了城市高压天然气管网不稳定流动和储、放气的变化规律。结合SCADA系统采集的实际运行数据,可对高压管线进行动态模拟,计算结果可指导高压天然气管线的设计、运行工况模拟和储放气的分析。该研究为高压天然气管网动态研究提供了新的思路,为高压天然气管网的设计、运行、调度调控提供了理论依据,为分析壅塞流动和城市天然气高压储气能力提供了新的方法,为高压天然气安全输送和储存提供了技术支持。 【关键词】：高压天然气管网 动态模拟 综合数值法 储气能力 壅塞流动
【学位授予单位】：山东大学
【学位级别】：博士
【学位授予年份】：2010
【分类号】：TE832
【目录】：

• 摘要13-18
• ABSTRACT18-25
• 主要符号说明25-28
• 第1章 绪论28-36
• 1.1 课题研究的背景28-31
• 1.2 世界高压天然气管线的发展状况31-32
• 1.2.1 国外高压天然气管线发展31
• 1.2.2 国内高压天然气管线发展31-32
• 1.3 国内外高压天然气管网动态模拟研究进展32-34
• 1.3.1 高压天然气管网动态模拟理论的发展32-33
• 1.3.2 高压天然气管网动态模拟工具的开发33-34
• 1.4 本文主要内容34-36
• 第2章 高压天然气管道流动特性与热力特性36-52
• 2.1 引言36
• 2.2 高压天然气管道的流动特性36-42
• 2.2.1 高压天然气在管道中流动基本方程36-40
• 2.2.2 柯列勃洛克公式中摩阻系数确定40-42
• 2.3 高压天然气管道热力特性42-51
• 2.3.1 高压天然气管道流动主要热力状态参数的求解42-48
• 2.3.2 高压天然气管道沿程温度变化规律的分析48-51
• 2.4 本章小结51-52
• 第3章 高压天然气管网动态模型的建立及求解52-93
• 3.1 利用特征线法建立动态模型及求解52-72
• 3.1.1 引言52-53
• 3.1.2 特征线方程的建立53-59
• 3.1.3 特征差分方程的建立59-63
• 3.1.4 边界条件的确定63-67
• 3.1.5 迭代计算过程67-71
• 3.1.6 节点方程的建立71-72
• 3.2 利用隐式中心差分法建立动态模型及求解72-87
• 3.2.1 引言72-73
• 3.2.2 差分格式方程的建立73-76
• 3.2.3 差分格式方程的边界条件和初始条件76-78
• 3.2.4 对差分方程组封闭性的分析78-79
• 3.2.5 非线性方程组的Newton-Raphson迭代法求解79-85
• 3.2.6 Jacobi矩阵的形成过程分析85-87
• 3.3 利用特征线法与有限差分法相结合的综合法求解87-90
• 3.3.1 引言87-88
• 3.3.2 综合法的适用性88-90
• 3.3.3 综合法计算的程序流程图90
• 3.4 本章小结90-93
• 第4章 多点接气高压天然气管网动态模型的赋值 #66及其封闭性分析93-104
• 4.1 引言93
• 4.2 多点接气高压天然气管网系统动态模拟初值的确定93-96
• 4.2.1 管段剖分面上流量初始值的确定94
• 4.2.2 管段剖分面上压力初始值的确定94-95
• 4.2.3 管段剖分面温度初始值的确定95-96
• 4.3 多点接气高压天然气管网动态模拟方程组封闭性分析96-100
• 4.3.1 从管段剖分面划分上进行分析96-97
• 4.3.2 从约束条件和边界条件上进行分析97-100
• 4.4 实例计算验证100-103
• 4.5 本章小结103-104
• 第5章 高压天然气管道壅塞流动模型的建立与特性分析104-120
• 5.1 引言104
• 5.2 壅塞发生的原理104-105
• 5.3 壅塞流动基本参数105-107
• 5.3.1 壅塞流动中声速的计算105-107
• 5.3.2 壅塞流动中马赫数的计算107
• 5.4 壅塞流动数学模型建立与分析107-118
• 5.4.1 等径管道壅塞流动特性107-110
• 5.4.2 变径管道壅塞流动特性110-114
• 5.4.3 实例验证壅塞流动模型与分析114-118
• 5.5 本章小结118-120
• 第6章 高压天然气管道储气能力确定及壅塞影响分析120-132
• 6.1 引言120
• 6.2 流动状态为稳定流时储气量的计算分析120-124
• 6.2.1 稳定流高压天然气管道储气能力的确定120-124
• 6.2.2 稳定流动状态下管道储气能力影响因素分析124
• 6.3 不稳定流动状态下储气能力计算分析124-130
• 6.3.1 不稳定流动状态下储气量的确定125
• 6.3.2 壅塞流动对管道储气能力的影响125-130
• 6.4 本章小结130-132
• 第7章 高压天然气管网动态模拟软件的开发与应用132-154
• 7.1 引言132
• 7.2 高压天然气管网动态模拟软件功能的改进132-137
• 7.2.1 存储与查询功能的改进132-136
• 7.2.2 数据库功能的改进136-137
• 7.3 模拟软件可靠性对比验证与模拟计算应用137-152
• 7.3.1 模拟结果与国外同类软件计算结果对比137-144
• 7.3.2 模拟结果与SCADA系统采集的实际运行数据对比144-150
• 7.3.3 改进版软件的模拟计算应用150-152
• 7.4 本章小结152-154
• 第8章 全文总结和展望154-159
• 8.1 本文主要结论154-157
• 8.2 本文创新点157-158
• 8.3 工作展望158-159
• 附录一至附录五159-163
• 参考文献163-170
• 致谢170-171
• 攻读博士学位期间发表的学术论文及科研成果171-173
• ENGLISH DISSERTATION173-205
• 学位论文评阅及答辩情况表205

您可以在本站搜索以下学术论文文献来了解更多相关内容

枝状高压天然气管网动态分析及软件开发    黄明

我国天然气资源的可持续发展与利用    李宁,马庆元,郭继平

微型喷管内气体流动的流量壅塞现象    丁英涛;谢君堂;仲顺安;

等截面摩擦管中气粉流的壅塞现象    韩珽，彭一川

城市高压燃气管网动态模拟计算    马江平;李帆;管延文;

对加快发展我国天然气供气调峰设施的思考    王莉,庄建远,檀建超

2008年中国天然气行业发展综述    杨建红;李伟;王峰;

2009年我国油气管道新进展    蒲明;

沿途有分输的长输燃气管道不稳定流计算    隋元春,薛世达

压缩天然气供应工艺及规模的探讨    申粤,乔珩,徐正康

长输管道末段储气量的计算与分析    李猷嘉

中国能源与经济之间关系的模型及实证分析    侯建朝

天然气管网数值模拟    李越

天然气管网的水力瞬变分析    于恩禄

基于设计选型及现场实测风机单位风量耗功率的分析    余鹏;唐辉强;

求解扩散方程的一类显示交替分组方法    王晨;徐安农;蒋心学;

风管流量测量实验的设计    方武;

空气悬架弹性元件特性的理论研究    杨启耀;张文娜;王存宝;周孔亢;

样条最小二乘法解泊松方程    王造奇

建平县大管径城市供水管路水锤防治措施    刘淑珍;苏建伟;黄志刚;

高压天然气管道孔口泄漏扩散浓度与范围仿真探讨    朱彦凝;吴赟城;王强;

移动进样口液相制备色谱分离过程参数的回归分析    李彤,林炳昌

用预测－校正方法解非线性含扩散的色谱方程    王际达，张雄文

ICP刻蚀机反应腔室气流仿真研究    程嘉;朱煜;

燃气管道瞬态仿真的研究与应用    高兰周;

油气管道维抢修技术进展    张仕民;谭桂斌;朱霄霄;串俊刚;陈迎春;

FDTD结合Mur二阶吸收边界条件分析三维散射问题时的稳定性研究    张玉胜;汪文秉;

MEI结合有限差分法二维应用的几个技术    张新军;梁昌洪;金谋平;

长输气管道采用内壁覆盖层减阻试验研究    翟建军;尹崇禄;林发布;魏中磊;张伯寅;林竹;张丽萍;秦延龙;杨晓鸿;

极超音速强激波作用下真实气体效应分析    沈钦灿;陶钢;

封闭式聚能装药爆炸噪声控制填充材料影响研究    徐全军;陈威;秦军;李裕春;孔福利;

可持续发展概念的广谱分析    张玉祥;

夏热冬冷地区地源热泵技术应用问题的几点思考    姜晨光;林长胜;钱保国;何跃平;朱烨昕;顾持真;

焊接车间污染源放散与控制的研究    刘秋新;张琳琳;

产业集群与循环经济的关系研究    张伟

中国能源消费与经济增长关系研究    吴明明

热解、气化过程中燃料-N的形态转化及迁移规律研究    王宗华

基于无网格方法的声学问题数值模拟研究    李坤

不锈钢冲压焊接离心泵能量转换特性与设计方法    张翔

水润滑轴承数值仿真及其材料摩擦学性能研究    段海涛

重型汽车动力传动系热管理系统研究    谭罡风

曲线公路隧道营运通风关键参数研究    王峰

致密气藏的储层特征与流体识别    倪根生

基于管网计算理论的通用锅炉水动力计算系统的研究与应用    谢金芳

水下超声速准二维垂直气体射流研究    郭强

增压器噪声控制与进气消声器设计研究    刘丽媛

基于TRNSYS的空调系统自动控制仿真器的开发    迟光亮

矿井风机监测系统的研究与设计    郑宝鑫

基于PLC的清洗机控制研究    陈永超

架空导线覆冰时捕获空气中过冷水滴规律研究    郭华

西部地区可再生能源的立法保障研究    杨艳艳

居住建筑室内设计中的建筑节能策略研究    丁翠

燃气管网的稳态分析与模拟    张勇

碳载Au基低温燃料电池阳极催化剂电化学性能研究    李红彬

不稳定流计算的探讨    瞿宏寿,黄一苓,朱谊章,刘壁澄

沿途有分输的长输燃气管道不稳定流计算    隋元春,薛世达

燃气管网节点方程牛顿法数学模拟    管延文,张鸣,李波,李帆

不同模型燃气管道动态模拟的对比分析    李军;玉建军;

燃气管网水力计算程序的开发    张宁;田贯三;李成乐;解东来;

浅析软件开发项目中的需求分析    朱寅非;

燃气管网拓扑结构的计算机生成方法及评价    彭继军,田贯三,张增刚,李兴泉,陈洪涛

天然气输配管网的矩阵描述    吕光宇,陈晓鸣

天然气在复杂枝状管网系统内不稳定流动及其计算机模拟    江茂泽;曾自强;刘良坚;

燃气管网动态仿真的研究及应用    杨昭;张甫仁;朱强;

天然气及西气东输工程有关情况介绍    李玉发

中国的能源效率及国际比较    王庆一

中国工业化中的能源消耗强度变动及因素分析——基于分解模型的实证分析    吴巧生;成金华;

中国能源强度变动的区域因素分解分析——基于LMDI分解方法    李国璋;王双;

我国工业CO_2排放影响因素差异性研究——基于高耗能行业与中低耗能行业    查冬兰;周德群;

中国能源消费因果关系分析    张志柏;

世界天然气发展前景探讨    孙国华,栾振华

中国能源强度变化的结构与效率因素贡献——基于AWD的实证分析    周勇;李廉水;

等截面管道中气固两相流的颗粒速度的理论计算    彭一川,肖泽强

GPRS技术在SCADA系统中的应用    韩冰,李芬华

中国工业能源消耗与CO_2排放影响效应研究    查冬兰;周德群;

天然气中油含量的测定    徐梅;邹林虎;

SCADA系统在天然气管网中的应用    郭强,莫德举

苏州:收拾停当迎“西气”——访苏州天然气管网有限公司总经理顾金明    岳勘

长三角地区天然气利用前景    李斌杰;

海南省环岛天然气管网规划(2008～2020年)启动    苏欣;

中国石油与山东省合资建设天然气管网    本刊通讯员;

热烈庆祝大型天然气管网国产化仿真软件产品正式发布    

天然气调压站工艺方案设计实践    卢祥林;

浙江省人民政府办公厅转发省计委关于加强天然气利用工作意见的通知    

天然气管网合理运作模式研究    罗佐县;

天然气管网压力能用于废旧橡胶粉碎的制冷装置    熊永强;华贲;罗东晓;徐文东;

利用天然气管网压力能的轻烃分离方法    申安云;余祖珊;王保全;

大型天然气管网系统抗震可靠性分析    刘武;谷雪琴;赫德明;

大型天然气管网系统抗震可靠性分析    刘武;谷雪琴;赫德明;

埋地管道防腐检测技术在燃气管网中的应用探析    楚洁;韩鹏辉;

西部开发的东进序曲——西气东输管道工程的力学论证    庄茁;由小川;冯耀荣;霍春勇;庄传晶;

21世纪城市燃气管网安全检测技术的回顾与展望    Tim G.Browm;李英杰;

多管齐下 突出事故预防 控制风险 提高本质安全    

钢骨架塑料复合管研制与应用    王俊良;赵勤宽;

燃气热水锅炉的工程设计初探    宋红伟;张娟;

江西省天然气管网工程开工建设    姜小毛周益民

永济天然气管网开工建设    皮海峰　行勇安

云南首个城市天然气管网开建    首席记者 吴晓松

海西天然气管网二期工程开工    黄西雷　郭永生

市城区天然气管网覆盖8万余人    本报记者 万亮发

三镇天然气管网主体工程竣工    记者 冯劲松　通讯员 王思敏

江西省天然气管网工程开建    记者 桂榕何宝庆

临汾一市四县(区)天然气管网建设项目开工    李玟

天然气管网年底通至平谷    记者 王东亮

我市高压天然气管网发展规划通过专家评审    记者 郑代玉

高压天然气管网动态模拟与壅塞流动特性研究    刘庆堂

大型天然气管网仿真计算引擎的研究与实现    郑建国

城市天然气输配管网水力模拟研究与实践    王兴畏

城市地下燃气管网泄漏监测关键技术研究    高华

城镇天然气管网优化设计研究    王鑫

城镇天然气管网优化系统的研究    陈晓妮

苏州天然气管网股份有限公司发展战略分析    熊非

天然气管网运行优化技术研究与软件开发    于进

基于GIS的城市天然气管网预警与应急系统研究    马磊

山东济华天然气管网远程监控管理系统的设计与实现    李亚松

多气源天然气管网水力工况研究    周晓静

转输天然气管网用气负荷研究与预测    薛晋

天然气管网数值模拟    李越

燃气蒸汽联合循环电站与天然气管网耦合特性研究    杨文正