CN109800968B - 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 - Google Patents
考虑天然气系统热力过程的电-气互联系统概率能流分析方法 Download PDFInfo
- Publication number
- CN109800968B CN109800968B CN201811633617.1A CN201811633617A CN109800968B CN 109800968 B CN109800968 B CN 109800968B CN 201811633617 A CN201811633617 A CN 201811633617A CN 109800968 B CN109800968 B CN 109800968B
- Authority
- CN
- China
- Prior art keywords
- node
- natural gas
- model
- gas
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了考虑天然气系统热力过程的电‑气互联系统概率能流分析方法,主要步骤为:1)对电‑气互联系统进行抽样。2)建立考虑天然气系统热力过程的设备模型。3)建立电‑气互联系统的概率能流模型,电‑气互联系统的概率能流模型包括节点平衡方程和不确定因素模型。4)建立风险指标模型。5)将样本数据输入到电‑气互联系统的概率能流模型,计算电‑气互联系统的能流,并将结果输入到风险指标模型中,计算电‑气互联系统的风险指标。6)收敛判断。本发明提出了一种较全面地考虑了天然气系统内存在的各种设备,以及各设备内热力过程的概率能流计算分析方法,可以更加准确、有效的分析系统状态并评估风险情况。
Description
技术领域
本发明涉及电-气综合能源系统领域,具体是考虑天然气系统热力过程的电-气互联系统概率能流分析方法。
背景技术
近年来,随着燃气轮机装机容量的日益攀升,电力系统和天然气系统之间的耦合日益紧密,将电-气互联系统看作一个统一整体来进行计算分析的做法日益普遍。在电-气互联系统中,存在着大量的不确定因素,这使得概率能流计算方法称为了一种有效的分析工具,以反映电压、气压等状态变量的越限风险等。
目前,在进行电-气互联系统的统一计算分析时,往往只考虑天然气压力和流量是变量,而天然气温度被当成和环境温度相等的常数,天然气系统中各设备的热力过程没有被计及。热力过程指的是当气体流过某设备时,热、功和各种能量间的传递与转换关系,这种关系可以通过气温、气压、流量等状态变量来刻画。
基于两个原因,在电-气互联系统的概率能流的计算中,天然气温度应该被当做变量。
其一,天然气的温度和压力间存在耦合关系,两变量会相互影响,且同一系统中不同部位的气温可能差异巨大(例如,天然气温度在中亚-中国天然气管道的首端为45℃,而在该管道的末端,天然气温度接近环境温度,其在冬天低至5℃以下),若在计算时将天然气温度看成常数,将给气压等计算结果都带来显著的误差。
其二,天然气系统存在一种名为“水合物生成”的特殊现象,该现象是否出现是由气温和气压共同决定的。若系统中生成了水合物,会带来一系列的不良后果,如局部气压过高、设备损坏和供气中断等,由于电力系统与天然气系统紧密耦合,因此供电的可靠性也会受到影响。此外,调压器作为天然气系统中的一种重要设备,通常设置在输气系统的末端,通过固定出口节点的气压来满足下游用户的要求。当天然气流经调压器时,气温和气压同时都会降低,这很有可能加重水合物生成的风险。因此,热力过程在准确分析气压等系统状态变量和评估水合物生成概率等系统风险时具有重要意义。
电-气互联系统的概率能流分析大体可以分为两步:能流模型的建立和风险指标的计算。对于能流模型,天然气设备模型是重要组成部分。首先,已有文献考虑了管道的热力过程并建立了节点热力平衡方程,但矛盾的地方在于,该研究使用Weymouth方程作为管道的流量模型,而事实上Weymouth方程是基于整个管道温度不变的假设推导的,因此在考虑热力过程的背景下使用该模型的计算精度有限。其次,对于调压器,目前在电-气互联系统的研究中,仅有一些文献在优化计算时将调压比上下限纳入约束,或根据其一端气压固定的特性,将调压器与压缩机一并归类为主动支路,而尚未见到对调压器的详细建模。再次,对于压缩机,其热力过程将导致出口温度的显著上升,而这一特性在现有的电-气互联系统研究中都尚未得到有效考虑。总而言之,现有研究对天然气系统设备及其热力过程的考虑都不够全面和精细。
对于指标,现有研究可以计算节点电压、气压的越限概率以及支路潮流的过载概率等,但仍缺乏可以量化水合物生成(同时受天然气温度和压力两个状态变量的影响)概率的指标。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,考虑天然气系统热力过程的电-气互联系统概率能流分析方法,主要包括以下步骤:
1)对电-气互联系统进行抽样,得到样本数据。所述样本数据包括电负荷、气负荷和风速。
2)建立考虑天然气系统热力过程的设备模型。所述设备模型包括天然气管道模型、压缩机模型和调压器模型。
所述天然气管道模型包括天然气管道热力模型和天然气管道流量模型。
所述压缩机模型包括压缩机热力模型和压缩机流量模型。
所述调压器模型包括调压器热力模型和调压器流量模型。
建立天然气管道模型的主要步骤如下:
2.1)天然气水平管道m1n1的气体流动守恒方程如下所示:
2.2)与管道起点位置距离x处的天然气温度T(x)如下所示:
式中,为计算系数。Ts为环境温度。Tm1为从天然气系统节点m1流出的天然气的温度。ηJT为焦耳-汤姆逊系数。pm1和pn1分别为天然气管道支路中首节点m1与末节点n1的气压。为天然气管道m1n1的长度。
2.3)将天然气温度T作为状态变量,则公式1改写为下式:
2.4)将公式2和公式3带入公式4中,得到:
2.5)标准状态下体积流量与质量流量的转换关系如下所示:
气体常数R如下所示:
2.6)将公式6和公式7带入公式5中,得到适用于非等温条件的代数形式的管道流量模型,即:
式中,T0为标准状态下天然气温度。p0为标准状态下天然气压力。λ为摩擦系数。
2.7)根据公式5,热力模型,即天然气管道m1n1出口处天然气温度的计算式如下所示:
压缩机热力模型如下所示:
调压器热力模型如下所示:
式中,Tm3为从天然气系统调压器节点m3流出的天然气的温度;为调压器m3n3出口处气体温度;pm3和pn3分别为天然气调压器支路中首节点m3与末节点n3的气压;为调压器孔径和管道内径的比值;为调压器系数。
调压器流量模型如下所示:
压缩机流量模型分别如公式13和公式14所示:
3)建立电-气互联系统的概率能流模型。电-气互联系统的概率能流模型包括节点平衡方程和不确定因素模型。节点平衡方程包括天然气系统模型、电力系统模型和耦合元件模型。不确定因素包括电力系统节点i的普通有功负荷PD,i、电力系统节点i的普通无功负荷QD,i、天然气系统节点n的普通气负荷FD,n和电力系统节点i的风电出力PW,i。不确定因素中电/气负荷被建模为正态分布,风速被建模为威布尔分布。
天然气系统模型包括天然气系统的节点热力平衡模型和天然气系统的节点流量平衡模型。
天然气系统节点的热力平衡如下所示:
式中,m∈n表示节点m与节点n相邻。FG,n和TG,n分别为天然气系统节点n的气源注入流量和该气源天然气的温度。m=m1,m2,m3;n=n1,n2,n3;
符号函数sgn1(m,n)的取值如下所示:
天然气系统节点的流量平衡模型如下所示:
式中,FD,n和FGAS,n分别为天然气系统节点n的普通气负荷和燃气轮机消耗的气流量。Nm为天然气系统节点的总数。
符号函数sgn2(m,n)的取值如下所示:
电力系统模型包括电力系统节点有功功率平衡方程和无功功率平衡方程。
电力系统节点有功功率平衡方程如下所示:
式中,PG,i和PGAS,i为电力系统节点i的常规机组和燃气轮机的有功出力。PD,i为电力系统节点i的普通有功负荷。PW,i为电力系统节点i的风电出力。
电力系统节点无功功率平衡方程如下所示:
式中,QG,i和QGAS,i分别为电力系统节点i的常规机组和燃气轮机的无功出力。QD,i为电力系统节点i的普通无功负荷。QC,i为电力系统节点i的并联无功补偿装置的出力。Vi和θi分别为电力系统节点i的电压幅值和相角。Gij和Bij为节点导纳矩阵的第i行第j列元素的实部和虚部。Ne为电力系统节点总数。θij为电力系统节点i和电力系统节点j的相角差。
耦合元件燃气轮机的能量转换模型如下所示:
式中,GHV为天然气高热值。αi,n、βi,n和γi,n为连接电力系统节点i和天然气系统节点n的燃气轮机的能量转换参数。
4)建立风险指标模型。所述风险指标模型主要包括水合物生成概率风险指标模型、气压越上限概率风险指标模型、气压越下限概率风险指标模型和气温越上限概率风险指标模型。
建立风险指标模型的主要步骤如下:
4.1)计算节点n的水合物生成概率PHFn,即:
式中,Pn(x)为节点n处于系统状态x的概率。H1,n为节点n生成水合物的状态的集合。X为节点n生成水合物的状态。
4.2)计算节点n的气压越上限概率PNHPVn,即:
式中,H2,n为节点n气压越上限的状态集合。
4.3)计算节点n的气压越下限概率PNLPVn,即:
式中,H3,n为节点n气压越下限的状态集合。
4.4)计算节点n的气温越上限概率PNHTVn,即:
式中,H4,n为节点n气温越上限的状态集合。
5)将样本数据输入到电-气互联系统的概率能流模型,计算电-气互联系统的能流,并将结果输入到风险指标模型中,计算电-气互联系统的风险指标。
6)判断是否达到最大抽样次数k,若是,则结束,若否,则返回步骤1。
本发明的技术效果是毋庸置疑的。本发明提出了一种较全面地考虑了天然气系统内存在的各种设备,以及各设备内热力过程的概率能流计算分析方法,可以更加准确、有效的分析系统状态并评估风险情况。
附图说明
图1为水合物生成的临界条件曲线图;
图2为电-气互联系统概率能流分析流程图;
图3为13节点天然气系统示意图;
图4为天然气系统各节点的气压越上限概率图;
图5为天然气系统各节点的气压越下限概率图;
图6为天然气系统各节点气压的标准差示意图;
图7为PHFn与风电渗透率的关系示意图;
图8为PNHPVn与风电渗透率的关系示意图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图2,考虑天然气系统热力过程的电-气互联系统概率能流分析方法,主要包括以下步骤:
1)利用非时序蒙特卡洛模拟法对电-气互联系统进行随机抽样,得到样本数据。所述样本数据包括电负荷、气负荷和风速。本实施例中涉及的所有随机样本都是连续型的。
2)建立考虑天然气系统热力过程的设备模型。所述设备模型包括天然气管道模型、压缩机模型和调压器模型。
所述天然气管道模型包括天然气管道热力模型和天然气管道流量模型。
所述压缩机模型包括压缩机热力模型和压缩机流量模型。
所述调压器模型包括调压器热力模型和调压器流量模型。
建立天然气管道模型的主要步骤如下:
2.1)天然气水平管道m1n1的气体流动守恒方程如下所示:
2.2)与管道起点位置距离x处的天然气温度T(x)如下所示:
式中,为计算系数。Ts为环境温度。Tm1为从天然气系统节点m1流出的天然气的温度。ηJT为焦耳-汤姆逊系数。pm1和pn1分别为天然气管道支路中首节点m1与末节点n1的气压。为天然气管道m1n1的长度。
在计算公式2时,由于右边的第三项的值远小于前面两项,
可以忽略右边的第三项。
2.3)将天然气温度T作为状态变量,则公式1改写为下式:
2.4)将公式2和公式3带入公式4中,得到:
2.5)标准状态下体积流量与质量流量的转换关系如下所示:
气体常数R如下所示:
2.6)将公式6和公式7带入公式5中,得到适用于非等温条件的代数形式的管道流量模型,即:
式中,T0为标准状态下天然气温度。P0为标准状态下天然气压力。λ为摩擦系数。
2.7)根据公式5,热力模型,即天然气管道m1n1出口处天然气温度的计算式如下所示:
压缩机热力模型如下所示:
式中,Tm2为从天然气系统压缩机节点m2流出的天然气的温度;为压缩机m2n2出口处气体温度;pm2和pn2分别为天然气压缩机支路中首节点m2与末节点n2的气压;为多变指数,其表明压缩机的实际热力过程是介于两种理想状况(绝热压缩和等温压缩)之间的。
调压器热力模型如下所示:
式中,Tm3为从天然气系统调压器节点m3流出的天然气的温度;为调压器m3n3出口处气体温度;pm3和pn3分别为天然气调压器支路中首节点m3与末节点n3的气压;为调压器孔径和管道内径的比值;为调压器系数。
调压器流量模型如下所示:
压缩机流量模型分别如公式13和公式14所示:
3)建立电-气互联系统的概率能流模型。电-气互联系统的概率能流模型包括节点平衡方程和不确定因素模型。节点平衡方程包括天然气系统模型、电力系统模型和耦合元件模型。不确定因素包括电力系统节点i的普通有功负荷PD,i、电力系统节点i的普通无功负荷QD,i、天然气系统节点n的普通气负荷FD,n和电力系统节点i的风电出力PW,i。不确定因素中电/气负荷被建模为正态分布,风速被建模为威布尔分布。
天然气系统模型包括天然气系统的节点热力平衡模型和天然气系统的节点流量平衡模型。
热力平衡反映的是,从某节点流出的天然气的温度,是流入该节点的各支路末端气体温度的加权平均。流量平衡表明,流入某节点和流出它的天然气流量应该是相等的。
天然气系统节点的热力平衡如下所示:
式中,m∈n表示节点m与节点n相邻。FG,n和TG,n分别为天然气系统节点n的气源注入流量和该气源天然气的温度。m=m1,m2,m3;n=n1,n2,n3;m1和n1表示天然气管道节点,m2和n2表示天然气压缩机节点。m3和n3表示天然气调压器节点。
m和n的取值包括以下三组:1)m=m1,且n=n1。2)m=m2,且n=n2。3)m=m3,且n=n3。
当m=m1,且n=n1时,建立的是电-气综合能源系统的天然气管道节点求解模型。当m=m2,且n=n2时,建立的是电-气综合能源系统的天然气压缩机节点求解模型。当m=m3,且n=n3时,建立的是电-气综合能源系统的天然气调压器节点求解模型。
符号函数sgn1(m,n)的取值如下所示:
天然气系统节点的流量平衡模型如下所示:
式中,FD,n和FGAS,n分别为天然气系统节点n的普通气负荷和燃气轮机消耗的气流量。Nm为天然气系统节点的总数。
符号函数sgn2(m,n)的取值如下所示:
电力系统模型包括电力系统节点有功功率平衡方程和无功功率平衡方程。
电力系统节点有功功率平衡方程如下所示:
式中,PG,i和PGAS,i为电力系统节点i的常规机组和燃气轮机的有功出力。PD,i为电力系统节点i的普通有功负荷。PW,i为电力系统节点i的风电出力。
电力系统节点无功功率平衡方程如下所示:
式中,QG,i和QGAS,i分别为电力系统节点i的常规机组和燃气轮机的无功出力。QD,i为电力系统节点i的普通无功负荷。QC,i为电力系统节点i的并联无功补偿装置的出力。Vi和θi分别为电力系统节点i的电压幅值和相角。Gij和Bij为节点导纳矩阵的第i行第j列元素的实部和虚部。Ne为电力系统节点总数。θij为电力系统节点i和电力系统节点j的相角差。
耦合元件燃气轮机的能量转换模型如下所示:
式中,GHV为天然气高热值。αi,n、βi,n和γi,n为连接电力系统节点i和天然气系统节点n的燃气轮机的能量转换参数。
所建的电-气综合能源系统统一能流模型是一组以X=[θi,Vi,πm,,Tm]T作为系统状态变量的非线性方程。使用牛顿-拉夫逊法可以实现该模型的有效求解。
4)所述风险指标模型主要包括水合物生成概率风险指标模型、气压越上限概率风险指标模型、气压越下限概率风险指标模型和气温越上限概率风险指标模型。
建立风险指标模型的主要步骤如下:
4.1)计算节点n的水合物生成概率PHFn,即:
式中,Pn(x)为节点n处于系统状态x的概率。H1,n为节点n生成水合物的状态的集合。X为节点n生成水合物的状态。
4.2)气压过高可能造成管道破裂,为避免这种情况,节点气压需设置一个上限,并统计计算节点n的气压越上限概率(PNHPVn)。当气压过低时会对气流的正常传输造成影响,因此节点气压还需设置一个下限,并统计计算节点n的气压越下限概率(PNLPVn)。此外,过高的气温还会损坏系统设备中防腐材料的性能,因此节点气温需设置一个上限,并统计计算节点n的气温越上限概率(PNHTVn)。
计算节点n的气压越上限概率PNHPVn,即:
式中,H2,n为节点n气压越上限的状态集合。
计算节点n的气压越下限概率PNLPVn,即:
式中,H3,n为节点n气压越下限的状态集合。
计算节点n的气温越上限概率PNHTVn,即:
式中,H4,n为节点n气温越上限的状态集合。
5)将样本数据输入到电-气互联系统的概率能流模型,计算电-气互联系统的能流,并将结果输入到风险指标模型中,计算电-气互联系统的风险指标。
6)判断是否达到最大抽样次数k,若是,则结束,若否,则返回步骤1。
实施例2:
一种验证考虑天然气系统热力过程的电-气互联系统概率能流分析方法的实验,主要包括以下步骤:
1)建立测试系统
以IEEE14-NGS13系统,即IEEE14节点系统和13节点天然气系统为例,测试本发明所提出的考虑天然气系统中不同设备热力过程的电-气综合能源系统概率能流分析方法,参见图3。
13节点的天然气系统包括2个气源,3个普通负荷,2台燃气轮机,8条普通管道,3台压缩机和一条管道+调压器的组合支路。其中,所有压缩机均由天然气驱动。调压器被设置在连接节点12和13的管道的末端。2台燃气轮机分别连接的是电力系统的节点1(平衡节点)和天然气系统的节点8,以及电力系统的节点3和天然气系统的节点6。其示意图如图2所示。电力系统节点6与一个风电场相连。
2)不同比较模型
为验证本发明所提考虑热力过程的电-气互联系统概率能流计算分析方法的有效性,采用如下3个模型进行比较:
M1:现有的考虑了管道热力过程的电-气互联系统稳态能流模型。
M2:区别于M1,使用本发明所推导的管道流量模型和首次考虑的压缩机热力模型。
M3:本发明提出的考虑天然气系统热力过程的电-气互联系统概率能流计算模型,在M2的基础上考虑调压器的存在。
3)管道流量模型和压缩机热力模型的有效性验证
本实施例仿真将模型M1和M2的计算结果进行了对比,在仿真时保持电力系统的风电渗透率为30%。表1给出的是采用不同的模型时,能够实现计算收敛的样本的比例。
表1M1和M2的计算收敛比例
从表1可以看出,采用M1进行计算时有3.21%的样本无法收敛,而当采用M2时所有样本都能够成功计算。造成这种差异的原因是,在M1所用的管道模型中存在一项ln((Tm-Ts)/(Tn-Ts)),若Tm和Tn其中一项比环境温度Ts大而另一项比Ts小,对数函数的自变量小于零,从而无法实现计算。
图4、图5和表2给出了采用不同的模型时,各节点水合物的生成概率以及节点气压的越限概率。表2中没有给出的各节点的水合物生成概率均为0。
表2M1和M2的水合物生成概率
从图4和图5可以看出,采用不同的模型进行计算,得到的节点气压越限概率具有显著差别,造成这种差异的原因是有二。其一,M1和M2对压缩机首末两端的气温的处理方式不同,对于M1,认为压缩机首末两端节点的温度为与气源温度相等的常量;而对于M2,压缩机首末两端节点温度是根据相应的模型计算得到的,其值可能与气源温度相差很大。其二,各节点的气温、气压等状态变量间存在紧密的耦合关系,任何变量计算结果的不同都会给其他变量的计算结果带来影响。表2表明,采用M1计算得到的PHFn和采用M2的PHFn具有显著差别,在天然气系统节点13,这种差异达到了18.02%。根据图1得知,当天然气温度高于275K时,水合物生成的临界压力随温度的上升而急剧上升,因此,哪怕分别采用M1和M2的气温计算结果差异不是很大,对水合物的生成概率也会产生很大的影响。
总而言之,采用不同的模型会得到不同的分析结果,因此,合理、准确的建模具有重要意义,本发明建立了考虑天然气系统热力过程的管道和压缩机模型,更接近设备的实际运行情况,具有较高的实用价值。
4)调压器模型的有效性验证
本实施例的仿真将模型M2和M3的计算结果进行了对比,在仿真时保持电力系统的风电渗透率为30%。图6给出了当调压器出口节点气压被固定为620psi时,系统各节点气压的标准差。
标准差代表了节点气压的离散程度。从图6可以看出,受负荷随机性和状态变量间耦合关系的影响,采用M2计算时,系统下游节点的气压标准差较大,而系统上游节点的气压受平衡节点的束缚作用而离散程度较小。然而,调压器的出口节点气压是固定的,它可以起到类似平衡节点的作用,因此,采用M3计算时调压器附近节点的气压标准差较小,如节点11和12。此外,采用M3计算时各节点气压的标准差都比采用M2时显著减小,这证明调压器对全网的节点气压具有重要的稳定作用。
表3M2和M3的水合物生成概率
表3给出了分别采用M2和M3计算时,节点水合物的生成概率。可以看出,考虑调压器的作用时,节点13的水合物生成概率升高了约28%。这是因为当天然气通过调压器时,其出口节点的气温和气压会同时降低,综合考虑这两个变量的下降程度和水合物生成临界条件的关系,该位置的水合物生成概率往往是会升高的。
5)电力系统对天然气系统的影响分析
本实施例仿真以M3作为仿真对象,对比了电力系统风电渗透率从0%到35%变化时,对天然气系统的影响。图7表示的是不同风电渗透率下,天然气系统节点8的水合物生成概率,图中的μ和δ表示的是风速的均值和标准差。可以看出,天然气系统节点的水合物生成概率随电力系统风电渗透率的增加而增加,且随风电场风速大小和波动性的增加而增加。
图8表示的是不同风电渗透率下,天然气系统节点5的气压越上限概率。可以看出,天然气系统节点的气压越限概率随电力系统风电渗透率的增加而增加,且增加的速度还与风电场风速的大小和波动性有关。这些结果体现了电力系统和天然气系统的紧密耦合关系,因此,把电-气互联系统看作一个统一整体进行分析是十分有必要的。
Claims (1)
1.考虑天然气系统热力过程的电-气互联系统概率能流分析方法,其特征在于,主要包括以下步骤:
1)建立电-气互联系统的概率能流模型;电-气互联系统的概率能流模型包括节点平衡方程和不确定因素模型;节点平衡方程包括天然气系统模型、电力系统模型和耦合元件燃气轮机的能量转换模型;
不确定因素包括电力系统节点i的普通有功负荷PD,i、电力系统节点i的普通无功负荷QD,i、天然气系统节点n的普通气负荷FD,n和电力系统节点i的风电出力PW,i;不确定因素中电/气负荷被建模为正态分布,风速被建模为威布尔分布;
利用非时序蒙特卡洛模拟法对电-气互联系统进行抽样,得到所述样本数据;所述样本数据包括电负荷、气负荷和风速;
2)建立考虑天然气系统热力过程的设备模型;所述设备模型包括天然气管道模型、压缩机模型和调压器模型;
所述天然气管道模型包括天然气管道热力模型和天然气管道流量模型;
所述压缩机模型包括压缩机热力模型和压缩机流量模型;
所述调压器模型包括调压器热力模型和调压器流量模型;
建立天然气管道模型的主要步骤如下:
2.1)天然气水平管道m1n1的气体流动守恒方程如下所示:
2.2)与管道起点位置距离x处的天然气温度T(x)如下所示:
式中,为计算系数;Ts为环境温度;Tm1为从天然气系统节点m1流出的天然气的温度;ηJT为焦耳-汤姆逊系数;pm1和pn1分别为天然气管道支路中首节点m1与末节点n1的气压;为天然气管道m1n1的长度;
2.3)将天然气温度T作为状态变量,则公式(1)改写为下式:
2.4)将公式(2)和公式(3)带入公式(4)中,得到:
2.5)标准状态下体积流量与质量流量的转换关系如下所示:
气体常数R如下所示:
2.6)将公式(6)和公式(7)带入公式(5)中,得到适用于非等温条件的代数形式的管道流量模型,即:
式中,T0为标准状态下天然气温度;p0为标准状态下天然气压力;λ为摩擦系数;
压缩机热力模型如下所示:
调压器热力模型如下所示:
式中,Tm3为从天然气系统调压器节点m3流出的天然气的温度;为调压器m3n3出口处气体温度;pm3和pn3分别为天然气调压器支路中首节点m3与末节点n3的气压;为调压器孔径和管道内径的比值;为调压器系数;
调压器流量模型如下所示:
压缩机流量模型如公式(13)和公式(14)所示:
3)建立天然气系统模型,包括天然气系统的节点热力平衡模型和天然气系统的节点流量平衡模型;
天然气系统的节点热力平衡如下所示:
式中,m∈n表示节点m与节点n相邻;FG,n和TG,n分别为天然气系统节点n的气源注入流量和该气源天然气的温度;m=m1,m2,m3;n=n1,n2,n3;
符号函数sgn1(m,n)的取值如下所示:
天然气系统节点的流量平衡模型如下所示:
式中,FD,n和FGAS,n分别为天然气系统节点n的普通气负荷和燃气轮机消耗的气流量;Nm为天然气系统节点的总数;
符号函数sgn2(m,n)的取值如下所示:
电力系统模型包括电力系统节点有功功率平衡方程和无功功率平衡方程;
电力系统节点有功功率平衡方程如下所示:
式中,PG,i和PGAS,i为电力系统节点i的常规机组和燃气轮机的有功出力;PD,i为电力系统节点i的普通有功负荷;PW,i为电力系统节点i的风电出力;
电力系统节点无功功率平衡方程如下所示:
式中,QG,i和QGAS,i分别为电力系统节点i的常规机组和燃气轮机的无功出力;QD,i为电力系统节点i的普通无功负荷;QC,i为电力系统节点i的并联无功补偿装置的出力;Vi和θi分别为电力系统节点i的电压幅值和相角;Gij和Bij为节点导纳矩阵的第i行第j列元素的实部和虚部;Ne为电力系统节点总数;θij为电力系统节点i和电力系统节点j的相角差;
耦合元件燃气轮机的能量转换模型如下所示:
式中,GHV为天然气高热值;αi,n、βi,n和γi,n为连接电力系统节点i和天然气系统节点n的燃气轮机的能量转换参数;
4)将样本数据输入到电-气互联系统的概率能流模型中,计算电-气互联系统的能流;
5)建立风险指标模型;所述风险指标模型包括水合物生成概率、气压越上限概率、气压越下限概率和气温越上限概率;
建立风险指标模型的主要步骤如下:
5.1)计算节点n的水合物生成概率PHFn,即:
式中,Pn(x)为节点n处于系统状态x的概率;H1,n为节点n生成水合物的状态的集合;X为节点n生成水合物的状态;
5.2)计算节点n的气压越上限概率PNHPVn,即:
式中,H2,n为节点n气压越上限的状态集合;
5.3)计算节点n的气压越下限概率PNLPVn,即:
式中,H3,n为节点n气压越下限的状态集合;
5.4)计算节点n的气温越上限概率PNHTVn,即:
式中,H4,n为节点n气温越上限的状态集合;
6)将电-气互联系统的能流输入到风险指标模型中,计算电-气互联系统的风险指标;电-气综合能源系统统一能流模型是一组以X=[θi,Vi,pm,Tm]T作为系统状态变量的非线性方程;
7)收敛判断,即判断是否达到最大抽样次数k,若是,则结束,若否,则返回步骤1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633617.1A CN109800968B (zh) | 2018-12-29 | 2018-12-29 | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811633617.1A CN109800968B (zh) | 2018-12-29 | 2018-12-29 | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109800968A CN109800968A (zh) | 2019-05-24 |
CN109800968B true CN109800968B (zh) | 2020-11-06 |
Family
ID=66558025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811633617.1A Active CN109800968B (zh) | 2018-12-29 | 2018-12-29 | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109800968B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110532642B (zh) * | 2019-08-07 | 2021-09-07 | 西安交通大学 | 一种综合能源系统概率能流的计算方法 |
CN110717643B (zh) * | 2019-08-07 | 2023-04-07 | 重庆大学 | 一种基于全局灵敏度分析的天然气网络储气配置方法 |
CN110929405B (zh) * | 2019-11-28 | 2023-08-04 | 国网辽宁省电力有限公司经济技术研究院 | 一种计及风电机组及燃气轮机组的电-气动态分析方法 |
CN111682520B (zh) * | 2020-03-11 | 2023-08-22 | 国网辽宁省电力有限公司大连供电公司 | 用网络流对电-气互联系统优化分析的方法 |
CN111950135B (zh) * | 2020-07-28 | 2022-09-23 | 重庆大学 | 一种基于网络解耦的电-气互联系统概率能流计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764689A (zh) * | 2015-04-07 | 2015-07-08 | 中国科学院广州能源研究所 | 一种天然气水合物合成和力学测试一体化压力室 |
CN106447113A (zh) * | 2016-10-08 | 2017-02-22 | 东南大学 | 一种基于运行优化模型的多区域综合能源系统运行方法 |
CN107194055A (zh) * | 2017-05-17 | 2017-09-22 | 重庆大学 | 考虑电转气装置的电‑气互联系统可靠性建模及其评估方法 |
CN107730129A (zh) * | 2017-10-24 | 2018-02-23 | 重庆大学 | 考虑光热热电联产与电锅炉的电‑气‑热互联系统风险评估方法 |
CN108734391A (zh) * | 2018-05-08 | 2018-11-02 | 重庆大学 | 基于堆栈降噪自动编码器的电-气综合能源系统概率能流量计算方法 |
CN109063379A (zh) * | 2018-09-11 | 2018-12-21 | 东北电力大学 | 基于最大熵原理的电-气综合能源系统概率能流计算方法 |
-
2018
- 2018-12-29 CN CN201811633617.1A patent/CN109800968B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104764689A (zh) * | 2015-04-07 | 2015-07-08 | 中国科学院广州能源研究所 | 一种天然气水合物合成和力学测试一体化压力室 |
CN106447113A (zh) * | 2016-10-08 | 2017-02-22 | 东南大学 | 一种基于运行优化模型的多区域综合能源系统运行方法 |
CN107194055A (zh) * | 2017-05-17 | 2017-09-22 | 重庆大学 | 考虑电转气装置的电‑气互联系统可靠性建模及其评估方法 |
CN107730129A (zh) * | 2017-10-24 | 2018-02-23 | 重庆大学 | 考虑光热热电联产与电锅炉的电‑气‑热互联系统风险评估方法 |
CN108734391A (zh) * | 2018-05-08 | 2018-11-02 | 重庆大学 | 基于堆栈降噪自动编码器的电-气综合能源系统概率能流量计算方法 |
CN109063379A (zh) * | 2018-09-11 | 2018-12-21 | 东北电力大学 | 基于最大熵原理的电-气综合能源系统概率能流计算方法 |
Non-Patent Citations (6)
Title |
---|
Analysis of probabilistic energy flow for integrated electricity-gas energy system with P2G based on cumulant method;Qiang Hu;《2017 IEEE Conference on Energy Internet and Energy System Integration》;20180104;全文 * |
Equivalent Optimal Power Flow Method Considering Natural Gas Network Constraints;Wei Dai;《2018 5th International Conference on Electric Power and Energy Conversion Systems》;20180923;全文 * |
Risk assessment of integrated electrical, natural gas and district heating;Juan Yu;《Electrical Power and Energy Systems》;20180615;第277~287页 * |
基于可靠性的天然气管道水合物形成概率;阮超宇;《油气储运》;20170531;第36卷(第5期);第487~494页 * |
湿气管道水合物形成概率计算;邓道明等;《化工学报》;20140226;第65卷(第6期);第2270~2275页 * |
电 - 气互联综合能源系统多时段暂态能量流仿真;卫志农;《电力自动化设备》;20170630;第37卷(第6期);第41~47页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109800968A (zh) | 2019-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109800968B (zh) | 考虑天然气系统热力过程的电-气互联系统概率能流分析方法 | |
CN108494021B (zh) | 电-热-气综合能源系统的稳定评估与静态控制方法 | |
CN109242365B (zh) | 一种电-热互联综合能源系统的区间潮流计算方法 | |
Zhu et al. | An integrated multi-energy flow calculation method for electricity-gas-thermal integrated energy systems | |
CN109740242B (zh) | 考虑天然气热力过程的电-气综合能源系统统一能流计算方法 | |
CN107730129B (zh) | 考虑光热热电联产与电锅炉的电-气-热互联系统风险评估方法 | |
Sundar et al. | State and parameter estimation for natural gas pipeline networks using transient state data | |
CN107194055B (zh) | 考虑电转气装置的电-气互联系统可靠性建模及其评估方法 | |
Fu et al. | Failure probability estimation of gas supply using the central moment method in an integrated energy system | |
US11181874B2 (en) | N-1 static security analysis method for integrated energy system | |
CN110851960B (zh) | 一种计及天然气网络影响的配电网运行可靠性评估方法 | |
CN109376428B (zh) | 综合能源系统的可靠性评估方法、装置、设备及存储介质 | |
CN111563315B (zh) | 一种基于拓扑分析的电-气综合能源系统稳态能量流计算方法 | |
Kou et al. | Many-objective optimization for coordinated operation of integrated electricity and gas network | |
CN112069692B (zh) | 一种天然气管网输差计算的优化求解方法 | |
CN110532642B (zh) | 一种综合能源系统概率能流的计算方法 | |
CN110866213A (zh) | 电-气综合能源系统多网络稳态能量流分析方法及装置 | |
CN111241479B (zh) | 基于交叉熵及客观熵权法的电-热互联综合能源系统风险评估方法 | |
CN108389136A (zh) | 一种考虑多重不确定性的天然气概率-模糊能流分析方法 | |
CN113283077B (zh) | 考虑相关性的电-热-氢综合能源系统区间能流计算方法 | |
Bao et al. | Probabilistic energy flow and risk assessment of electricity–gas systems considering the thermodynamic process | |
CN111950135B (zh) | 一种基于网络解耦的电-气互联系统概率能流计算方法 | |
Wang et al. | Finite-time stabilization of Port-controlled Hamiltonian systems with application to nonlinear affine systems | |
CN112257355A (zh) | 一种掺氢天然气的中低压配气管网建模方法 | |
CN111414675A (zh) | 一种电热综合能源系统的双层抗差状态估计方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |