大型天然气管网系统在线仿真方法及软件开发

　　摘要：天然气管网系统仿真技术是进行管网系统设计方案论证、管输计划制定、运行方案优化、应急保障决策的核心技术。基于质量守恒、动量守恒、能量守恒原理及非管元件特性方程，建立了适用于任意结构形式的管网系统仿真模型;形成了基于管道数学模型泛函分析、大型稀疏矩阵压缩存储条件下的快速 LU 分解等方法的管网系统仿真模型求解算法。采用 B/S 软件架构、大数据高速缓存技术、大数据抽取实时可视化技术，开发了包括数据通讯模块、数据库存储模块、数据滤波模块、在线仿真模块、结果推送模块的大型天然气管网系统在线仿真软件。应用于总长 100 km、45 座压缩机站的管网系统和总长 31 000 km、117 座压缩机站的大型管网系统，软件适用规模、计算精度与速度均达到了国际同类商业软件水平。研究成果可为国产化仿真软件的开发应用提供借鉴。

　　关键词：天然气管网;在线仿真;数学模型;软件开发

　　天然气管网系统是由气源、管网、用户及储气库等“源 网 荷 储”各要素构成的规模庞大、组成复杂的一体化水动力系统，是重要的能源输送基础设施。截至 020 年，中国天然气干线管道总里程约 0.2 10km，基本形成了“西气东输、北气南下、海气登陆、就近供应”的供气格局 2]。

　　随着中国天然气管网系统规模不断增大、结构日趋复杂、智慧化要求逐步提高，如何实现管网系统的经济高效输送、安全稳定运行、灵活可靠调配，是管网设计运行面临的重要问题，对管网系统仿真软件的适用规模、计算精度与速度也提出了更高的要求。天然气管网系统的计算机仿真技术可基于对管网流动过程的数学表征，量化管网内天然气的流动及变化过程，从而精确预测和回溯管网系统各气源、用户、管道、设备的运行状态，是实现管网系统安全、高效设计与运行管理的核心技术，也是智慧管网背景下支撑智能管输计划制定、方案优化、应急保障的基础4]。

　　国外目前已经基本形成了较为成熟的管网系统仿真软件产品，如美国 Gregg 公司的 NextGen 软件、德国 DNV GL 的 SPS 软件、美国 Emerson 公司的 Pipeline Studio 软件、捷克 SIMONE Research Group 和德国 LIWACOM 公司联合开发的 SIMONE 软件 ，这些软件的广泛应用为中国管网系统的科学规划、建设、运行管理做出了重要贡献。在中国管网规模逐渐大型化、管网信息安全要求逐步提高、核心自主知识产权“卡脖子”问题日益凸显的背景下，国外软件的采购与维护费用高、深层次功能定制与开发、信息安全自主可控等方面的局限性也开始显现出来。特别是国外软件受其特殊版权与核心技术保护机制的影响，难以融入中国自主构建的智慧管网体系。

　　因此，掌握天然气管网系统仿真的核心理论方法，是中国在天然气干线管网和城市燃气管网规划、设计、运行、管理等方面面临的迫切需求。虽然中国相关科研院所在此方面进行了长期的研究，但相关软件产品缺乏良好的应用生态和成熟商业化发展模式，在软件适用的管网规模、计算速度与精度、稳定性方面与国外产品还存在差距。西南石油大学从 世纪 年代开始在中国较早开展管网系统仿真技术的研究，形成了大型管网系统仿真的核心理论方法 。

　　近年来，抓住国家和行业重视并推动国产仿真软件发展的契机，融合大数据、人工智能等新兴技术，实现了管网系统仿真软件架构从 C/S 架构向 /S 架构、仿真模式从离线仿真向在线仿真、适用规模从小型管网向大型复杂环状管网的快速发展。详细论述了管网系统建模及求解算法、大型天然气管网系统在线仿真软件开发及相关应用实例，以期为国产化仿真软件的开发应用提供借鉴。

　　1大型天然气管网系统仿真模型

　　对于非管道元件，从组成管网系统的实际元件出发，可建立通用阀、止回阀、调节阀、局部阻力件、换热器、过滤器、放空火炬、储罐、缓冲罐、安全阀、泄漏点、储气库、配气站、集气站、分输站、压缩机站的数学模型。考虑压缩机在实际运行过程中存在的压比控制、出口压力控制、转速控制等多种控制模式，调节阀的流量控制、压差控制、出口压力控制、入口压力控制等多种控制模式，换热器的温差控制、出口温度控制等多种模式，采用基于 PID 的控制方法，即可实现对各种控制模式的模拟。

　　为了使管网系统模型封闭，还需在管网系统的气源和用户处设置边界条件。其中，气源处一般设置压力温度边界条件、流量温度边界条件;用户处需设置压力边界条件、流量边界条件中任何一个 21 。如果是随时间变化的，还应给定这些参数随时间的变化关系。

　　2管网系统仿真模型的求解算法

　　天然气管网系统仿真模型是由一系列非线性偏微分方程组和非线性代数方程组构成，传统的求解方法首先采用隐式差分法、特征线法等方法将非线性偏微分方程展开为线性方程，然后采用 Newton Raphson迭代法、Gauss eidel 迭代法等方法进行求解，并在小型管网模型中得到了较为广泛的应用。

　　这类方法实际是用线性化方程的解来近似非线性方程的解，所以求解的准确性很大程度上取决于管道网格划分的长度，网格长度越小，计算结果越接近实际情况。对于大型天然气管网系统而言，随着元件、节点数量的增加，仿真模型规模急剧增大，如果管网网格长度划分过小，会使得数学模型规模激增，导致计算机求解速度缓慢、内存消耗急剧增加。

　　针对上述问题，结合管道流动方程的特点，基于泛函理论，创新性研究了采用大网格描述压力、温度、流量等工艺参数分布的有效基函数及其叠加形式，以不同基函数及其叠加形式来表征管道大网格内压力、流量等参数的分布，并通过研究这种分布函数来逼近偏微分方程的解 。

　　通过该方法扩大了管道网格步长，减少了仿真方程的数量，进而扩大了仿真规模。同时，针对大型管网系统仿真模型迭代求解过程中，雅克比矩阵为稀疏非对称矩阵的特点 ，将大型方程组分解成一系列小型方程组，采用稀疏矩阵的压缩存储法，利用压缩存储条件下的 LU 分解法、自适应仿真等多项技术进行快速求解，从而避免对稀疏矩阵中零元素的访问和求解。

　　测试表明，采用上述方法后，可节省计算机内存 95%以上，能够大大提升仿真软件对于管网规模的适应性。上述数学模型与求解方法对于管网系统拓扑结构、方程的分布形式没有特别要求，因此，可适用于任意结构形式的管网系统，构成了大型天然气管网系统在线仿真软件的计算核心。 大型天然气管网系统在线仿真软件开发以上述数学模型和求解算法为基础，结合天然气管网系统智能化运行的业务需求，设计了天然气管网系统运行管控平台，完成了管网在线仿真软件的开发。

　　3大型天然气管网系统在线仿真软件开发

　　基于 B/S架构，采用 Java EE 企业级开发平台开发，由图形化组态、实时数据通讯、在线数字滤波、大数据存储、在线仿真、参数评价与修正、大数据可视化模块构成;可实现大型天然气管网系统的在线仿真、运行参数预测与回溯、运行参数预警、系统能耗分析、清管作业分析、虚拟计量等功能 。

　　( 1)图形化组态模块。广泛借鉴国外管网系统仿真软件、绘图办公软件的风格，开发了全中文的图形化组态界面，具备与微软 Visio 软件类似的绘图风格，界面友好易用。

　　( 2)实时数据通讯模块。基于 HTTPS、HTTP 通讯协议及 RESTful 数据接口，实现了在线仿真软件与数字孪生体数据的集成、数字孪生体数据向仿真软件的实时传输、仿真结果向数字孪生体的实时推送。

　　( 3)在线数字滤波模块。数字孪生体采集的数据可能出现异常值或者缺失，如果将这些数据作为计算的边界条件输入在线仿真模块，将获得明显异常的仿真结果，甚至导致计算崩溃。为此，采用在线多尺度滤波方法(On line Multiscale Filtering，OLMS)将采集到的数据进行数字滤波处理，去除异常信号;采用预测参数、移动插值平均等方法进行参数补齐，保证仿真过程的稳健性。

　　(4 )大数据存储模块。针对数字孪生体与在线仿真软件数据传输量大及时序数据结构的特点，设计了管道属性数据库、天然气物性数据库、环境参数数据库、仿真结果数据库等在内的核心数据库。针对海量数据的接收与存储难题，首先采用高速缓存技术满足时序数据即时读写需求，再将数据持久化到硬盘。理论上可实现 00× 条 的数据读写要求，实现 000× 条数据分库存储，100× 条数据分表存储。

　　(5 )在线仿真模块。基于上述理论模型和求解算法，开发在线仿真软件内核，包括基于 BWRS、PR、SRK、GERG 2008 等状态方程的天然气热物性参数计算、管网水力与热力仿真、组分与热值跟踪、前景及条件预测、清管器跟踪等核心功能 29 。

　　( 6)参数评价与修正模块。为提高仿真软件的计算精度，根据仿真结果与实际运行参数之间的偏差，采用遗传算法与深度学习神经网络算法，修正管道的总传热系数、粗糙度、压缩机特性参数等可调参数，实现模型的自主学习和修正。( 7)大数据可视化模块。基于大数据可视化框架，实现实时采集数据、在线仿真数据、压缩机组能耗数据、管输计划及实际执行数据等多源大数据汇聚展示。

　　4应用实例

　　4.1西部管道天然气管网

　　西部管道天然气管网包括西气东输一线、二线、三线西段以及轮吐线，全长 100 km，包括 个主要气源， 个沿途分输点， 06 座阀室， 座压缩机站，292 条管道。西二线和西三线并行布置，站内压缩机通过管道、阀门连接，形成一个更为复杂的管网系统。

　　模型中，天然气物性计算采用 PR 状态方程，摩阻系数采用 Colebrook 公式，管道网格长度 1~10 km。气源采用流量温度边界条件，如西二线霍尔果斯站流量为 413 10 /d，入口温度为 0℃;分输用户采用压力边界或者体积流量边界条件，如西二线西段终点压力为 MP ;压缩机站采用压比或出口压力控制，如西二线霍尔果斯站压比为 1.5，站内压缩机的串并联配置按照实际压缩机组配置方式。

　　4.1.1 动态仿真

　　在稳态仿真基础上，改变南京分输站的流量，使其随时间不断变化，设置固定时间步长为 180 ，管道网格步长 0.1 10 km，进行 连续动态仿真，动态仿真完成时间为 05 ;进一步改变最小时间步长为min，最大时间步长 20 min，进行 的连续动态仿真，动态仿真完成时间为 60 。

　　通过与商业软件的对比可见，计算速度均优于 Pipeline Studio 4.0、SPS 9.7 等国外商业软件(如 SPS 软件 5 动态仿真用时 2 )，并与 SIMONE 软件( 5 动态仿真用时 20 )处于同一数量级，总体达到了国外商业软件的先进水平。采用上述软件模拟得到南京分输站及相邻扬子石化分输站 内压力、流量随时间的变化情况。可见，南京分输站压力、输量呈相反的变化趋势，而扬子石化分输站与南京分输站连接在一个阀室上，两者压力变化趋势相同，符合实际运行规律。

　　5 结论

　　开发具有自主知识产权的大型天然气管网系统在线仿真软件是掌握管网系统仿真核心技术、确保管网系统运行调度信息安全、推进天然气管网系统智能化运行的迫切需求。介绍了大型天然气管网系统仿真模型的构建方法、基于泛函分析的管网系统仿真模型高效求解算法，并基于 B/S 架构开发了由 大模块构成的大型天然气管网系统在线仿真软件，可实现管网系统的在线仿真、运行参数预测与回溯、运行参数预警、系统能耗分析、清管作业分析、虚拟计量等功能。

　　通过西部管道天然气管网及国家干线天然气管网系统的实例应用，证明该软件的计算精度符合要求，在计算速度及适用的管网系统规模方面也达到了国外同类软件的先进水平，为中国国产化高水平管网系统仿真软件的开发奠定了基础。未来可继续开展管道、站场在全生命周期内各种运行场景的在线仿真、预测、优化、预演、预警与决策，指导天然气管网的运行业务优化和应急管控。

　　参考文献：

　　[1] 刘合，梁坤，张国生，李志欣，丁麟，苏健，等 碳达峰、碳中和约束下我国天然气发展策略研究[J].中国工程科学，2021，23：33 42.LIU H, LIANG K, ZHANG G S, LI Z X, DING L, SU J, et al. China’s natural gas development strategy under the constraints ofcarbon peak and carbon neutrality[J]. Strategic Study of CAE, 2021, 23(6): 33 42.

　　[2] 李鹭光 中国天然气工业发展回顾与前景展望[J].天然气工业，2021，41： 11.LI L G. Development of natural gas industry in China: review and prospect[J]. Natural Gas Industry, 2021, 41(8): 1 11.

　　[3] MARINO A, ZIO E. A framework for the resilience analysis of complex natural gas pipeline networks from a cybe physicalsystem perspective[J]. Computers & Industrial Engineering, 2021, 162: 107727.

　　[4] 宇波，王鹏，王丽燕，向月 基于分而治之思想的天然气管网仿真方法[J].油气储运，2017，36：75 84.YU B, WANG P, WANG L Y, XIANG Y. A simulation method for natural gas pipeline networks based on the divide and conquerconcept[J]. Oil & Gas Storage and Transportation, 2017, 36(1): 75 84.

　　作者：李长俊 张员瑞 贾文龙 仇柏林 何杰 王硕

转载请注明来自发表学术论文网：http://www.fbxslw.com/dzlw/30002.html