一种电-气耦合多能流系统联合静态安全分析方法技术方案

本发明专利技术涉及一种利电‑气耦合多能流系统的联合静态安全分析方法，属于含多种能源形式的运行和控制技术领域。本方法突破电力系统安全分析只考虑电网本身的局限，考虑了电力供应和天然气供应的相互影响，实现了电‑气耦合多能流系统的联合静态安全分析，从而能够同时分析电力供应和天然气供应中各种可能的预想事故对两个系统的影响，相比现有分别对力供应和天然气供应进行独立安全分析的方法，提高了安全分析的精确性，减少了误判和漏判。该方法可以应用于电‑气耦合多能流系统的在线运行中，在电‑气耦合多能流能量管理中工作，当系统存在安全问题或潜在风险时给出安全警告，有助于提高电‑气耦合多能流系统运行的安全性。

【技术实现步骤摘要】

本专利技术涉及一种电-气耦合多能流系统联合静态安全分析方法，属于含多种能源形式的电力系统运行和控制

技术介绍
能源综合利用是提高综合能源利用效率、促进可再生能源消纳的重要途径，通过打破原来电、热、冷、气、交通等能流子系统相对割裂的状态，实现多类型能源开放互联，构建多能流系统。多能流是指多种类型的能量流，表示电、热、冷、气、交通等能量流的相互耦合、转换和传输。多能流系统相比传统相互割裂的能源系统，其带来的效益包括：1)通过多类型能源的梯级开发利用和智能管理，可以降低能源消耗和浪费，提高综合能源利用效率，并有助于减少总的用能成本；2)利用不同能源的特性差异和互补、转换，有助于提高消纳间歇式可再生能源的能力；3)通过多能源的转供、互补和协调控制，有助于提高供能的可靠性，并为电网的运行提供更多可调控资源；4)通过多能流系统的协同规划和建设，可以减少基础设施的重复建设和浪费，提高资产利用率。多能流系统一方面具有可观的效益，另一方面也使原本复杂的能源系统更加复杂。多能流系统由多个能流子系统组成，这些能流子系统之间相互作用和影响，使得多能流系统复杂度显著增加，体现出许多新的特性，传统各个能流单独分析的方法已经难以适应新的要求，亟需发展出新的多能流分析方法。在我国，越来越多的燃气轮机、电制氢等耦合元件客观上增强了电-气之间的互联，促进了电-气耦合多能流系统的发展，也对电-气耦合多能流系统的运行控制和能量管理提出了新的要求。安全是多能流系统正常运行的基础，随着系统的复杂化，多能流的安全问题也更加突出。多能流系统在运行过程中会受到各种扰动的影响，包括元件故障、发电(如风电)变化、负荷变化等。为了保障多能流系统的安全运行，在绝大部分扰动下仍然可以持续传输能量，需要对系统的运行状态进行安全分析。安全分析是能量管理的核心功能之一，可以在线发现多能流系统实时运行中的薄弱环节，有助于提高多能流系统运行的安全水平。现有的安全分析技术主要集中在电力系统中，尚未扩展到多能流系统中，目前尚无电-气耦合多能
流系统静态安全分析的方法。因此，需要研究电-气耦合多能流系统静态安全分析方法。
技术实现思路
本专利技术的目的是提出一种电-气耦合多能流系统静态安全分析方法，考虑电力供应和天然气供应的相互影响和各种可能的预想事故，建立电-气耦合多能流系统静态安全分析的框架和方法，实现电-气耦合多能流系统的联合静态安全分析。本专利技术提出的电-气耦合多能流系统联合静态安全分析方法，包括以下步骤：(1)从电-气耦合多能流系统的能量管理系统中，获得当前电-气耦合多能流系统的状态估计结果Ψse，包括：电-气耦合多能流系统中电网节点的电压幅值和电网支路的有功功率，电-气耦合多能流系统中气网节点压力和气网支路流量；(2)从电-气耦合多能流系统的能量管理系统中，获得电-气耦合多能流系统中电网、气网的安全运行约束Ψsc，包括：电-气耦合多能流系统中电网节点电压幅值允许最大值、电网节点电压幅值允许最小值、电网支路有功功率允许最大值和电网支路有功功率允许最小值，电-气耦合多能流系统中气网节点压力允许最大值、气网节点压力允许最小值、气网支路流量允许最大值和气网支路流量允许最小值；(3)根据上述步骤(1)的当前状态估计结果Ψse，对上述步骤(2)的安全运行约束Ψsc进行判断：(3-1)若所有当前状态估计结果Ψse均满足安全运行约束Ψsc，转步骤(4)；(3-2)若当前状态估计结果Ψse中存在不满足安全运行约束Ψsc的状态，则判定该电-气耦合多能流系统处于紧急状态，给出安全警告，并给出不满足安全运行的具体信息，安全分析方法结束；(4)从电-气耦合多能流系统的能量管理系统中获取预想事故集Ω；(5)对上述预想事故集Ω进行判断，若预想事故集Ω为空集，则转步骤(7)，若预想事故集Ω为非空集合，则从预想事故集Ω中取出一个预想事故，形成该预想事故下的电-气耦合多能流方程组，包括：(5-1)电-气耦合多能流系统中的电网潮流方程： P i = U i Σ j ∈ i U j ( G i j c o s ( θ i - θ j ) + B i j s i n ( θ i - θ j ) ) , ( i = 1 , 2 , 3 , ... , n ) ]]> Q i = U i Σ j ∈ i U j ( G i j sin ( θ i - θ j ) - B i j cos ( θ i - θ j ) ) , ( i = 1 , 2 , 3 , ... , n ) ]]>其中本文档来自技高网...

【技术保护点】
一种电‑气耦合多能流系统联合静态安全分析方法，其特征在于该方法包括以下步骤：(1)从电‑气耦合多能流系统的能量管理系统中，获得当前电‑气耦合多能流系统的状态估计结果Ψse，包括：电‑气耦合多能流系统中电网节点的电压幅值和电网支路的有功功率，电‑气耦合多能流系统中气网节点压力和气网支路流量；(2)从电‑气耦合多能流系统的能量管理系统中，获得电‑气耦合多能流系统中电网、气网的安全运行约束Ψsc，包括：电‑气耦合多能流系统中电网节点电压幅值允许最大值、电网节点电压幅值允许最小值、电网支路有功功率允许最大值和电网支路有功功率允许最小值，电‑气耦合多能流系统中气网节点压力允许最大值、气网节点压力允许最小值、气网支路流量允许最大值和气网支路流量允许最小值；(3)根据上述步骤(1)的当前状态估计结果Ψse，对上述步骤(2)的安全运行约束Ψsc进行判断：(3‑1)若所有当前状态估计结果Ψse均满足安全运行约束Ψsc，转步骤(4)；(3‑2)若当前状态估计结果Ψse中存在不满足安全运行约束Ψsc的状态，则判定该电‑气耦合多能流系统处于紧急状态，给出安全警告，并给出不满足安全运行的具体信息，安全分析方法结束；(4)从电‑气耦合多能流系统的能量管理系统中获取预想事故集Ω；(5)对上述预想事故集Ω进行判断，若预想事故集Ω为空集，则转步骤(7)，若预想事故集Ω为非空集合，则从预想事故集Ω中取出一个预想事故，形成该预想事故下的电‑气耦合多能流方程组，包括：(5‑1)电‑气耦合多能流系统中的电网潮流方程：Pi=UiΣj∈iUj(Gijcos(θi-θj)+Bijsin(θi-θj)),(i=1,2,3,...,n)]]>Qi=UiΣj∈iUj(Gijsin(θi-θj)-Bijcos(θi-θj)),(i=1,2,3,...,n)]]>其中：Pi为电网中节点i的注入有功功率，Qi为电网中节点i的注入无功功率，θi、θj分别为节点i、节点j的电压相角，Ui、Uj分别为节点i、节点j的电压幅值，Gij为电网节点导纳矩阵Y第i行、第j列元素的实部，Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部，电网节点导纳矩阵Y从电‑气耦合多能流系统的能量管理系统中获取；(5‑2)电‑气耦合多能流系统中的天然气网水力方程：fkm=sgnp(pk,pm)×(77.54Tnpn)(Dkm)2.5(1LkmγGTakmFZg)0.5(pk2-pm2)]]>其中：fkm为天然气网中的第k个节点和第m个节点之间的管道中天然气体积流量，pk，pm分别为第k个节点和第m个节点的压强，Dkm、Lkm分别为第k个节点和第m个节点之间管道km的管道直径和管道长度，F为管道内壁的摩擦系数，F由公式计算得到，Ef为管道的效率系数，Ef取值0.92，Re为雷诺数，由公式Re＝ρvd/μ计算得到，ρ为天然气密度，v为天然气流速，μ为天然气黏性系数，d为特征长度，γG为天然气比重，取值为0‑1，Ta为天然气平均温度，Tn、pn分别为天然气标准状态下的温度和压力，Tn、pn取值分别为288K、0.1Mpa，Zg为天然气平均可压缩系数，Zg的取值为0‑1，天然气网水力方程中，当时，上式中的sgnp(pk,pm)＝1，当时，sgnp(pk,pm)＝‑1；(5‑3)电‑气耦合多能流系统中的天然气网压缩机能量消耗方程：BHP=fCinckηC(ck-1)[(pfpe)ck-1ck-1]]]>其中：pf、pe分别为天然气网中压缩机首端节点和末端节点的压强，BHP为压缩机的能耗，为该压缩机入口的体积流量，ηc为压缩机的总效率，ck为压缩机的多变系数，ηc和ck从压缩机的出厂说明书获取；(5‑4)电‑气耦合多能流系统中的通过燃气轮机耦合的电力系统与天然气网之间的耦合方程：PTur=C1fTur3+C2fTur2+C3fTur]]>其中：fTur为燃气轮机的燃气体积流量，PTur为燃气轮机的有功功率，C1、C2和C3为燃气轮机的燃气系数，从燃气轮机的出厂说明书获取；(5‑5)电‑气耦合多能流系统中天然气网的节点体积流量平衡方程：AGf＝L，其中：AG为天然气网的节点‑支路关联矩阵，从电‑气耦合多能流系统的能量管理系统获取，f为天然气网的支路体积流量，L为天然气网的节点注入体积流量；(6)利用牛顿拉夫逊方法，求解上述电‑气耦合多能流方程组，得到该预想事故下的电‑气耦合多能流系统的多能流计算结果Ψc，包括电网节点电压幅值、电网支路有功功率、气网节点压力和气网支路流量，并将该预想事故从预想事故集Ω中删除；(7)根据步骤(5)的...

【技术特征摘要】
1.一种电-气耦合多能流系统联合静态安全分析方法，其特征在于该方法包括以下步骤：(1)从电-气耦合多能流系统的能量管理系统中，获得当前电-气耦合多能流系统的状态估计结果Ψse，包括：电-气耦合多能流系统中电网节点的电压幅值和电网支路的有功功率，电-气耦合多能流系统中气网节点压力和气网支路流量；(2)从电-气耦合多能流系统的能量管理系统中，获得电-气耦合多能流系统中电网、气网的安全运行约束Ψsc，包括：电-气耦合多能流系统中电网节点电压幅值允许最大值、电网节点电压幅值允许最小值、电网支路有功功率允许最大值和电网支路有功功率允许最小值，电-气耦合多能流系统中气网节点压力允许最大值、气网节点压力允许最小值、气网支路流量允许最大值和气网支路流量允许最小值；(3)根据上述步骤(1)的当前状态估计结果Ψse，对上述步骤(2)的安全运行约束Ψsc进行判断：(3-1)若所有当前状态估计结果Ψse均满足安全运行约束Ψsc，转步骤(4)；(3-2)若当前状态估计结果Ψse中存在不满足安全运行约束Ψsc的状态，则判定该电-气耦合多能流系统处于紧急状态，给出安全警告，并给出不满足安全运行的具体信息，安全分析方法结束；(4)从电-气耦合多能流系统的能量管理系统中获取预想事故集Ω；(5)对上述预想事故集Ω进行判断，若预想事故集Ω为空集，则转步骤(7)，若预想事故集Ω为非空集合，则从预想事故集Ω中取出一个预想事故，形成该预想事故下的电-气耦合多能流方程组，包括：(5-1)电-气耦合多能流系统中的电网潮流方程： P i = U i Σ j ∈ i U j ( G i j c o s ( θ i - θ j ) + B i j s i n ( θ i - θ j ) ) , ( i = 1 , 2 , 3 , ... , n ) ]]> Q i = U i Σ j ∈ i U j ( G i j s i n ( θ i - θ j ) - B i j c o s ( θ i - θ j ) ) , ( i = 1 , 2 , 3 , ... , n ) ]]>其中：Pi为电网中节点i的注入有功功率，Qi为电网中节点i的注入无功功率，θi、θj分别为节点i、节点j的电压相角，Ui、Uj分别为节点i、节点j的电压幅值，Gij为电
\t网节点导纳矩阵Y第i行、第j列元素的实部，Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部，电网节点导纳矩阵Y从电-气耦合多能流系统的能量管理系统中获取；(5-2)电-气耦合多能流系统中的天然气网水力方程： f k m = sgn p ( p k , p m ) × ( 77.54 T n p n ) ( D k m ...

【专利技术属性】
技术研发人员：孙宏斌，郭庆来，王彬，潘昭光，陈瑜玮，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京;11