本发明涉及一种用于综合能源系统运行控制的天然气气路建模方法，属于综合能源系统的运行控制技术领域。
技术背景
综合能源系统能够有效提高综合用能效率，已经成为国内外科学研究与工程实践的热点与前沿。综合能源系统的规划和运行以各个能源网络的建模和分析为基础，电力和天然气能源流紧密耦合，其中电力的分析基于从“场”到“路”的简化，已经形成了成熟的电路理论，而天然气气路的分析尚未形成与之统一的成熟理论。常规天然气气路建模方面尚存在的问题包括：缺乏直观的物理模型，可解释性不强；气-电耦合网络的分析方法无法统一，电力和天然气这两个学科之间存在知识壁垒；为保证求解精度，需要引入时空两个维度上的大量微元，面临计算复杂度高的难题。近些年来，“电路”理论的建模思想开始逐渐应用于天然气气路建模，但尚未形成完整而统一的理论框架，模型求解难度也较大，难以进一步推广到综合能源系统规划和运行的多样化应用中。因此，为实现不同能源网络研究的学科融合，以及促进综合能源系统的规划和运行工作的开展，亟需提出更加适合于综合能源系统的天然气气路模型。
技术实现要素：
本发明的目的是提出一种用于综合能源系统运行控制的天然气气路建模方法，以解决已有技术中存在的问题。基于天然气管道中质量守恒与动量守恒方程，以及天然气状态方程和流量方程，建立天然气管道中流量与压力之间的偏微分方程；利用傅里叶变换将气路映射至频域并通过二端口等值得到集总参数模型；结合天然气增压机方程，建立天然气气路一般支路模型；定义节点-支路关联矩阵和节点-流出支路关联矩阵，建立天然气气路的拓扑约束方程；结合天然气气路一般支路模型和天然气气路的拓扑约束方程，建立完整的天然气气路方程。
本发明提出的用于综合能源系统运行控制的天然气气路建模方法，包括以下步骤：
(1)建立天然气气路的管道模型，包括以下步骤：
(1-1)建立天然气在管道中一维流动过程的质量守恒方程和动量守恒方程：
式中：ρ、v和p分别为天然气的密度、流速和压力；λ、d和θ分别为管道的摩擦系数、内径和倾角，由天然气气路管理方提供，g为重力加速度，t和x分别为时间和空间；
(1-2)在步骤(1-1)的动量守恒方程中引入两个近似：一是忽略对流项，即
(1-3)将天然气状态方程p＝rtρ和管道流量方程g＝ρva代入质量守恒方程和简化后的动量方程中，得到管道中天然气流量与压力之间的时空偏微分方程：
式中：r和t分别为天然气的气体常数和温度，g为天然气质量流量，a为天然气管道的横截面积；
(1-4)建立天然气管道上一个微元的两端流量差和压降方程：
(1-5)根据步骤(1-4)中的微元的两端流量差和压降方程，定义天然气管道中气阻rg、气感lg、气容cg和受控气压源kg，rg、lg、cg和kg的计算方程如下：
rg＝λvb/(ad)
lg＝1/a
cg＝a/(rt)
从而，dx长度的管道表示为一段包括4个元件的气路，整个管道进而表示为一个分布参数气路；
(1-6)将步骤(1-5)中定义的rg、lg、cg和kg代入步骤(1-4)中的微元的两端流量差和压降方程，并通过傅里叶变换映射到频域后，获得每一个频率分量下的常微分方程如下：
并定义zg＝rg jwlg，yg＝jwcg；
(1-7)利用步骤(1-6)中的两个方程求解天然气管道末端的流量和气压为：
式中：gl和pl分别为天然气管道末端的天然气流量与压力，g0和p0分别为天然气管道首端的天然气流量和压力，l为天然气管道长度；
(1-8)定义天然气管道的传播系数为γgc＝zgyg，定义天然气管道的特征阻抗zgc＝zg/yg；
(1-9)根据步骤(1-7)中的天然气管道末端的流量和气压两个方程以及步骤(1-8)中两个定义，将天然气管道方程表示成线性二端口网络形式：
式中：a、b、c和d是网络参数，其值为：
(1-10)根据步骤(1-9)中二端口网络方程，建立π型等值气路，等值参数为：
z＝-b
k＝1-ad bc
y1＝(ad-bc-a)/b
y2＝(1-d)/b
(2)建立天然气气路的一般支路模型，包括以下步骤：
(2-1)建立天然气增压机的数学模型如下：
p1＝p2 eg
式中：p1和p2是天然气增压机两侧的压力，eg是天然气增压机提供的压力增量；
(2-2)根据步骤(1-5)中的气阻、气感、气容、受控气压源和步骤(2-1)中天然气增压机组成一般支路(如图3所示)，一般支路的方程如下：
gb＝yb(pb eb-kbpf)
式中：gb是支路中的流量，gb为未知量，pb是支路两端的气压差，pf为支路首端压力，若支路的首段为气源，则pf为已知量，若支路的首段不为气源，则pf为未知量，pt为支路末端的压力，pt为未知量，yb是气阻、气感和气容构成的支路导纳，kb和eb是支路中受控气压源和天然气增压机的元件参数，由天然气气路管理方提供；
(2-3)将天然气气路中所有支路的支路方程写成矩阵形式如下：
gb＝yb(pb eb-kbpf)
式中：gb为各个支路流量构成的向量，yb是各个支路导纳构成的对角矩阵，pb是各个支路两端气压差构成的向量，eb是各个天然气增压机的气压增量构成的向量，kb是各个受控气压源参数构成的向量，pf是各个支路首端气压构成的向量；
(3)建立天然气气路的拓扑约束方程，包括以下步骤：
(3-1)定义天然气气路中的节点-支路关联矩阵ag，该矩阵是一个n行m列的矩阵，其中n是节点数，m是支路数，用(ag)i,j表示其中第i行、第j列的元素，则(ag)i,j＝0表示支路j与节点i不相连，(ag)i,j＝1表示支路j从节点i流出，(ag)i,j＝-1表示支路j流入节点i；
(3-2)定义天然气气路中的节点-流出支路关联矩阵ag ，该矩阵保留了矩阵ag中的非负元素，即对于(ag )i,j，若支路j从节点i流出，则该元素为1，否则为0；
(3-3)建立天然气气路节点质量守恒方程：
aggb＝gn
式中：gn为每个节点上的流量注入构成的列向量，其中，天然气气路中的气负荷节点处的流量为已知量，气源节点处的流量为未知量，非气负荷和非气源的节点处的流量为0；
(3-4)建立天然气气路节点气压方程：
式中：pn为每个节点上的压力构成的列向量，其中，天然气气路中的气源节点处的压力为已知量，气负荷节点处的压力为未知量，非气负荷和非气源的节点处的压力为未知量；
(4)建立天然气气路方程，包括以下步骤：
(4-1)将步骤(3-3)和步骤(3-4)建立的方程代入步骤(2-3)建立的支路方程，得到未约简形式的天然气气路方程如下：
(4-2)定义广义节点导纳矩阵y′g和广义节点注入向量g′n如下：
g′n＝gn-agybeb
(4-3)将步骤(4-2)中定义的y′g和g′n代入(4-1)中的未约简形式的天然气气路方程，获得以下天然气气路模型方程：
y′gpn＝g′n
求解上述天然气气路模型，获得天然气气路中未知的节点压力，进而利用支路方程求取未知的支路流量，实现对综合能源系统的运行控制。
本发明提出的用于综合能源系统运行控制的天然气气路建模方法，其优点是：
本发明的用于综合能源系统运行控制的天然气气路建模方法，基于天然气管道中质量守恒与动量守恒方程，以及天然气状态方程和流量方程，建立天然气管道中流量与压力之间的偏微分方程；利用傅里叶变换将气路映射至频域并通过二端口等值得到集总参数模型；结合天然气增压机方程，建立天然气气路一般支路模型；定义节点-支路关联矩阵和节点-流出支路关联矩阵，建立天然气气路的拓扑约束方程；结合天然气气路一般支路模型和天然气气路的拓扑约束方程，建立天然气气路方程。本发明的天然气气路建模方法，与电力网络的网络矩阵和网络方程在数学形式上具有高度的统一性，从而奠定了气、电两种异质能流统一分析的基础。同时，相比传统分析方法，本发明方法具有更低的计算复杂度。
附图说明
图1是天然气管道的分布参数气路图，其中图1(a)为天然气整个管道的分布参数气路示意图，图1(b)为天然气管道微元dx的分布参数气路示意图。
图2是天然气管道的集总参数等值气路示意图。
图3是天然气气路中的一般支路示意图。
具体实施方式
本发明提出的用于综合能源系统运行控制的天然气气路建模方法，包括以下步骤：
(1)建立天然气气路的管道模型，包括以下步骤：
(1-1)建立天然气在管道中一维流动过程的质量守恒方程和动量守恒方程：
式中：ρ、v和p分别为天然气的密度、流速和压力；λ、d和θ分别为管道的摩擦系数、内径和倾角，由天然气气路管理方提供，g为重力加速度，t和x分别为时间和空间；
(1-2)在步骤(1-1)的动量守恒方程中引入两个近似：一是忽略对流项，即
(1-3)将天然气状态方程p＝rtρ和管道流量方程g＝ρva代入质量守恒方程和简化后的动量方程中，得到管道中天然气流量与压力之间的时空偏微分方程：
式中：r和t分别为天然气的气体常数和温度，g为天然气质量流量，a为天然气管道的横截面积；
(1-4)建立天然气管道上一个微元的两端流量差和压降方程：
(1-5)根据步骤(1-4)中的微元的两端流量差和压降方程，定义天然气管道中气阻rg、气感lg、气容cg和受控气压源kg，rg、lg、cg和kg的计算方程如下：
rg＝λvb/(ad)
lg＝1/a
cg＝a/(rt)
从而，dx长度的管道表示为一段包括4个元件的气路，整个管道进而表示为一个分布参数气路，天然气整个管道的分布参数气路和天然气管道微元dx的分布参数气路如图1中所示；
(1-6)将步骤(1-5)中定义的rg、lg、cg和kg代入步骤(1-4)中的微元的两端流量差和压降方程，并通过傅里叶变换映射到频域后，获得每一个频率分量下的常微分方程如下：
并定义zg＝rg jwlg，yg＝jwcg；
(1-7)利用步骤(1-6)中的两个方程求解天然气管道末端的流量和气压为：
式中：gl和pl分别为天然气管道末端的天然气流量与压力，g0和p0分别为天然气管道首端的天然气流量和压力，l为天然气管道长度；
(1-8)定义天然气管道的传播系数为γgc＝zgyg，定义天然气管道的特征阻抗zgc＝zg/yg；
(1-9)根据步骤(1-7)中的天然气管道末端的流量和气压两个方程以及步骤(1-8)中两个定义，将天然气管道方程表示成线性二端口网络形式：
式中：a、b、c和d是网络参数，其值为：
(1-10)根据步骤(1-9)中二端口网络方程，建立π型等值气路，该等值气路如图2所示，等值参数为：
z＝-b
k＝1-ad bc
y1＝(ad-bc-a)/b
y2＝(1-d)/b
(2)建立天然气气路的一般支路模型，包括以下步骤：
(2-1)建立天然气增压机的数学模型如下：
p1＝p2 eg
式中：p1和p2是天然气增压机两侧的压力，eg是天然气增压机提供的压力增量；
(2-2)根据步骤(1-5)中的气阻、气感、气容、受控气压源和步骤(2-1)中天然气增压机组成一般支路(如图3所示)，一般支路的方程如下：
gb＝yb(pb eb-kbpf)
式中：gb是支路中的流量，gb为未知量，通过求解天然气气路方程得到天然气气路中各个节点的压力后求得，pb是支路两端的气压差，pf为支路首端压力，若支路的首段为气源，则pf为已知量，若支路的首段不为气源，则pf为未知量，pt为支路末端的压力，pt为未知量，yb是气阻、气感和气容构成的支路导纳，kb和eb是支路中受控气压源和天然气增压机的元件参数，由天然气气路管理方提供；
(2-3)将天然气气路中所有支路的支路方程写成矩阵形式如下：
gb＝yb(pb eb-kbpf)
式中：gb为各个支路流量构成的向量，yb是各个支路导纳构成的对角矩阵，pb是各个支路两端气压差构成的向量，eb是各个天然气增压机的气压增量构成的向量，kb是各个受控气压源参数构成的向量，pf是各个支路首端气压构成的向量；
(3)建立天然气气路的拓扑约束方程，包括以下步骤：
(3-1)定义天然气气路中的节点-支路关联矩阵ag，该矩阵是一个n行m列的矩阵，其中n是节点数，m是支路数，用(ag)i,j表示其中第i行、第j列的元素，则(ag)i,j＝0表示支路j与节点i不相连，(ag)i,j＝1表示支路j从节点i流出，(ag)i,j＝-1表示支路j流入节点i；
(3-2)定义天然气气路中的节点-流出支路关联矩阵ag ，该矩阵保留了矩阵ag中的非负元素，即对于(ag )i,j，若支路j从节点i流出，则该元素为1，否则为0；
(3-3)建立天然气气路节点质量守恒方程：
aggb＝gn
式中：gn为每个节点上的流量注入构成的列向量，其中，天然气气路中的气负荷节点处的流量为已知量，气源节点处的流量为未知量，非气负荷和非气源的节点处的流量为0；
(3-4)建立天然气气路节点气压方程：
式中：pn为每个节点上的压力构成的列向量，其中，天然气气路中的气源节点处的压力为已知量，气负荷节点处的压力为未知量，非气负荷和非气源的节点处的压力为未知量；
(4)建立天然气气路方程，包括以下步骤：
(4-1)将步骤(3-3)和步骤(3-4)建立的方程代入步骤(2-3)建立的支路方程，得到未约简形式的天然气气路方程如下：
(4-2)定义广义节点导纳矩阵y′g和广义节点注入向量g′n如下：
g′n＝gn-agybeb
(4-3)将步骤(4-2)中定义的y′g和g′n代入(4-1)中的未约简形式的天然气气路方程，获得以下天然气气路模型方程：
y′gpn＝g′n
求解上述天然气气路模型，获得天然气气路中未知的节点压力，进而利用支路方程求取未知的支路流量，实现对综合能源系统的运行控制。
1.一种用于综合能源系统运行控制的天然气气路建模方法，其特征在于该方法包括以下步骤：
(1)建立天然气气路的管道模型，包括以下步骤：
(1-1)建立天然气在管道中一维流动过程的质量守恒方程和动量守恒方程：
式中：ρ、v和p分别为天然气的密度、流速和压力；λ、d和θ分别为管道的摩擦系数、内径和倾角，由天然气气路管理方提供，g为重力加速度，t和x分别为时间和空间；
(1-2)在步骤(1-1)的动量守恒方程中引入两个近似：一是忽略对流项，即
(1-3)将天然气状态方程p＝rtρ和管道流量方程g＝ρva代入质量守恒方程和简化后的动量方程中，得到管道中天然气流量与压力之间的时空偏微分方程：
式中：r和t分别为天然气的气体常数和温度，g为天然气质量流量，a为天然气管道的横截面积；
(1-4)建立天然气管道上一个微元的两端流量差和压降方程：
(1-5)根据步骤(1-4)中的微元的两端流量差和压降方程，定义天然气管道中气阻rg、气感lg、气容cg和受控气压源kg，rg、lg、cg和kg的计算方程如下：
rg＝λvb/(ad)
lg＝1/a
cg＝a/(rt)
从而，dx长度的管道表示为一段包括4个元件的气路，整个管道进而表示为一个分布参数气路；
(1-6)将步骤(1-5)中定义的rg、lg、cg和kg代入步骤(1-4)中的微元的两端流量差和压降方程，并通过傅里叶变换映射到频域后，获得每一个频率分量下的常微分方程如下：
并定义zg＝rg jwlg，yg＝jwcg；
(1-7)利用步骤(1-6)中的两个方程求解天然气管道末端的流量和气压为：
式中：gl和pl分别为天然气管道末端的天然气流量与压力，g0和p0分别为天然气管道首端的天然气流量和压力，l为天然气管道长度；
(1-8)定义天然气管道的传播系数为γgc＝zgyg，定义天然气管道的特征阻抗zgc＝zg/yg；
(1-9)根据步骤(1-7)中的天然气管道末端的流量和气压两个方程以及步骤(1-8)中两个定义，将天然气管道方程表示成线性二端口网络形式：
式中：a、b、c和d是网络参数，其值为：
(1-10)根据步骤(1-9)中二端口网络方程，建立π型等值气路，等值参数为：
z＝-b
k＝1-ad bc
y1＝(ad-bc-a)/b
y2＝(1-d)/b
(2)建立天然气气路的一般支路模型，包括以下步骤：
(2-1)建立天然气增压机的数学模型如下：
p1＝p2 eg
式中：p1和p2是天然气增压机两侧的压力，eg是天然气增压机提供的压力增量；
(2-2)根据步骤(1-5)中的气阻、气感、气容、受控气压源和步骤(2-1)中天然气增压机组成一般支路(如图3所示)，一般支路的方程如下：
gb＝yb(pb eb-kbpf)
式中：gb是支路中的流量，gb为未知量，pb是支路两端的气压差，pf为支路首端压力，若支路的首段为气源，则pf为已知量，若支路的首段不为气源，则pf为未知量，pt为支路末端的压力，pt为未知量，yb是气阻、气感和气容构成的支路导纳，kb和eb是支路中受控气压源和天然气增压机的元件参数，由天然气气路管理方提供；
(2-3)将天然气气路中所有支路的支路方程写成矩阵形式如下：
gb＝yb(pb eb-kbpf)
式中：gb为各个支路流量构成的向量，yb是各个支路导纳构成的对角矩阵，pb是各个支路两端气压差构成的向量，eb是各个天然气增压机的气压增量构成的向量，kb是各个受控气压源参数构成的向量，pf是各个支路首端气压构成的向量；
(3)建立天然气气路的拓扑约束方程，包括以下步骤：
(3-1)定义天然气气路中的节点-支路关联矩阵ag，该矩阵是一个n行m列的矩阵，其中n是节点数，m是支路数，用(ag)i,j表示其中第i行、第j列的元素，则(ag)i,j＝0表示支路j与节点i不相连，(ag)i,j＝1表示支路j从节点i流出，(ag)i,j＝-1表示支路j流入节点i；
(3-2)定义天然气气路中的节点-流出支路关联矩阵ag ，该矩阵保留了矩阵ag中的非负元素，即对于(ag )i,j，若支路j从节点i流出，则该元素为1，否则为0；
(3-3)建立天然气气路节点质量守恒方程：
aggb＝gn
式中：gn为每个节点上的流量注入构成的列向量，其中，天然气气路中的气负荷节点处的流量为已知量，气源节点处的流量为未知量，非气负荷和非气源的节点处的流量为0；
(3-4)建立天然气气路节点气压方程：
式中：pn为每个节点上的压力构成的列向量，其中，天然气气路中的气源节点处的压力为已知量，气负荷节点处的压力为未知量，非气负荷和非气源的节点处的压力为未知量；
(4)建立天然气气路方程，包括以下步骤：
(4-1)将步骤(3-3)和步骤(3-4)建立的方程代入步骤(2-3)建立的支路方程，得到未约简形式的天然气气路方程如下：
(4-2)定义广义节点导纳矩阵y′g和广义节点注入向量g′n如下：
g′n＝gn-agybeb
(4-3)将步骤(4-2)中定义的y′g和g′n代入(4-1)中的未约简形式的天然气气路方程，获得以下天然气气路模型方程：
y′gpn＝g′n
求解上述天然气气路模型，获得天然气气路中未知的节点压力，进而利用支路方程求取未知的支路流量，实现对综合能源系统的运行控制。
技术总结