CN105046369A - 一种基于能源中心的电气混联***建模和优化调度方法 - Google Patents
一种基于能源中心的电气混联***建模和优化调度方法 Download PDFInfo
- Publication number
- CN105046369A CN105046369A CN201510498181.XA CN201510498181A CN105046369A CN 105046369 A CN105046369 A CN 105046369A CN 201510498181 A CN201510498181 A CN 201510498181A CN 105046369 A CN105046369 A CN 105046369A
- Authority
- CN
- China
- Prior art keywords
- delta
- gas
- node
- energy
- dtri
- 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.)
- Granted
Links
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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Feedback Control In General (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公布了一种基于能源中心的电气混联***建模和优化调度方法,本发明首先建立了电力网络、天然气网络和能源中心模型,电力网络和天然气网络通过能源中心耦合形成电气混联***。然后以总的能源成本为目标函数、考虑了各种约束条件建立电气混联***优化调度数学模型。用原-对偶内点法求解,求解过程中先后引入松弛变量和障碍参数将模型变为只有等式约束的模型,再引入拉格朗日乘子得到拉格朗日函数,之后用牛顿法求解其KKT条件形成的非线性方程组。构造的算例仿真结果表明本发明对电气混联***的优化效果优于单独优化效果。
Description
技术领域
本发明涉及一种基于能源中心的电气混联***建模和优化调度方法,属于多能源综合控制利用领域。
技术背景
化石能源的广泛使用导致环境问题日益突出以及化石能源逐渐枯竭,传统的经济和社会发展模式难以持续,亟需建立更加高效、环保和可持续型的能源利用模式,能源中心(energyhub,EH)的提出使建立新的能源利用模式成为可能。能源中心被定义为由能源转化设备和储能设备构成、能够实现多种能源相互转化和存储的虚拟实体。能源中心作为多种能源形式的中间媒介,促使多种能源形式和能源***的联系更加紧密,并考虑了经济性影响。能源中心输入端连接多种能源的供应,输出端连接多种能源的消费,不同能源形式在其间相互转化,这就为多种能源形式的优化调度提供了可能性。
相比于其他一次能源,天然气对环境的影响较小、储量丰富且易于存储;燃气轮机和联合循环电厂的日渐发展,使得天然气发电投资更少、效率更高、建设周期缩短且调控更加灵活方便;此外,“页岩气革命”的不断深化,将造成天然气价格显著下降。可以预见,天然气发电的比例将逐步提高,天然气也会凭借其诸多优势在未来能源利用模式中占据重要地位。天然气网络与电力网络的联系越来越密切,可以说未来的能源利用模式将是天然气网络与电力网络高度耦合的产物。
前人对电力网络和天然气网络的独立网络都有足够的研究,如最优潮流(optimalpowerflow,OPF)等经典问题,但少有考虑电力网络与天然气网络的综合优化调度问题。近年来,由于天然气在能源利用中的地位越来越突出,电力网络与天然气网络的综合研究引起了国内外学者的广泛关注。电气混联***的研究主要分为协调规划和协调运行两个阶段,问题主要集中在不确定因素的处理、电负荷与电负荷的相互转化与影响、监管部门的协调规划与管理以及混联***的动态行为与可靠性等方面。
电气混联***优化调度本质上是非线性规划问题。本发明以总的能源成本作为目标函数,考虑了电力网络中有功、无功功率平衡约束,平衡节点相角平衡约束,发电机有功、无功出力约束,电压约束;天然气网络中流量平衡约束,压缩机状态量平衡约束,气源点输出约束和节点压力约束;能源中心的输入输出平衡约束和调度因子约束。本发明采用原-对偶内点法(primal-dualinteriorpointmethod,PDIPM)求解基于能源中心的电气混联***优化调度问题,因为PDIPM有着收敛性好,计算效率高,鲁棒性强,初值选取不敏感,没有识别约束集困难等优点。最后对构造的算例进行编程仿真,结果分析验证了本发明方法的优势。
发明内容
发明目的:本发明针对现有技术所需解决的技术问题提供一种基于能源中心的新型能源供应形式即电气混联***的建模和优化调度。
技术方案:本发明为实现上述目的,采用如下技术方案:
本发明为一种基于能源中心概念电气混联***建模和优化调度方法,其特征在于所述方法依次按以下步骤实现:
1)获取电力网络的参数信息,包括:输电线路的首端节点和末端节点编号,支路π型等效电路的电阻、电抗,对地并联电导、电纳,变压器变比和阻抗,各节点负荷以及发电机输出有功、无功约束,各节点电压约束;
对于各电力网络中节点i:
其中:Pi,Qi分别为电力网络中节点i有功功率和无功功率;ei,fi分别为节点i电压向量的实部和虚部;Gij,Bij分别为节点导纳矩阵第i行第j列元素的实部和虚部。
2)获取天然气网络的参数信息,包括:输气管道的首端节点和末端节点编号,管道的长度、内径和传输效率等物理特性,压缩机的状态量,各节点气负荷以及气源点天然气输出约束,各节点压力约束;
引入包括管道流量方程、压缩机流量消耗方程及流量平衡方程;
理想条件下,管道k从节点i到节点j的流量值可以用下面的方程表示:
其中:fkij为管道流量值; Fk为管道摩擦系数;Dk为管道内径;G为气体重力系数;Lk为管道长度;πi为节点i压力值;πj为节点j压力值;π0为标准压力值;T0为标准温度值;Tka为平均气体温度;Za为平均气体压缩系数;
对于高压网络的全紊流形式,流量方程可以进一步简化为:
其中: ε为管道效率;
理想气体条件下,压缩机的能量消耗方程可以表示为:
其中:fCk为通过压缩机的气体流量;πi为气体注入压缩机压力;πj为气体输出压缩机压力;Zki为压缩机入口的气体压缩系数;Tki为压缩机输出温度;α为热量系数;ηk为压缩机效率;
转化为消耗的流量值:
τk=αTk+βTkHkij+γTkHk 2 ij;
每一个节点的流量平衡方程可以用下面的矩阵形式表示:
(A+U)f+w-Tτ=0;
其中:A、U和T分别为支路、压缩机和压缩机消耗与节点有关的方向矩阵,f为支路流量值向量;w为各节点的气体注入向量;τ为各压缩机消耗流量值向量;
3)获取能源中心的参数信息,包括:各能源中心输入输出与电力网络及天然气网络的连接情况,各能源中心负荷以及各能源中心调度因子约束;
多输入多输出能源中心模型可以用矩阵方程描述:
其中:P和L分别为输入输出向量;矩阵C被称为耦合矩阵;
列出能源中心输入输出平衡方程:
其中:ηT表示变压器效率;ηGTe表示燃气轮机气转电的效率;ηGTh表示燃气轮机气转热的效率;ηF表示燃气锅炉的效率;
4)以总的能源成本作为目标函数f(x),选取控制变量u和状态变量根据各种约束条件列写等式约束和不等式约束,建立电气混联***非线性规划数学模型并用PDIPM求解;
建立电气混联***优化调度模型:
minf(x);
s.t.h(x)=0;
gmin≤g(x)≤gmax;
其中:h表示等式约束;g表示不等式约束;gmax,gmin分别表示不等式约束的上、下限;
5)引入广义松弛变量l和u,将模型中不等式约束转化为等式约束;
引入广义松弛变量l=[l1,…,lr]T,u=[u1,…,ur]T将广义不等式转化为等式:
g(x)+u-gmax=0;
g(x)-l-gmin=0;
应满足:u>0,l>0;
6)为了保证目标函数在可行域内取得最小值,对目标函数引入对数障碍参数μ,将模型转化为只含有等式约束的模型;
引入松弛变量和障碍参数将模型转化成只有等式约束的模型:
s.t.h(x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
其中:扰动因子(障碍参数)μ>0;
引入拉格朗日乘子得到拉格朗日函数:
7)用拉格朗日乘子法求解上述等式约束模型,引入拉格朗日乘子z、w和y,得到拉格朗日函数,根据拉格朗日极值存在必要条件(KKT条件)列写非线性方程组;
根据拉格朗日极值存在条件(KKT条件)得到一组非线性方程组,用牛顿-拉弗逊法求解,其一阶修正方程为:
其中: Lx,Ly,Lz,Lw,为KKT方程上次迭代的残差;和分别为h(x)和g(x)的海森矩阵;
8)定义对偶间隙Gap=lTz-uTw和扰动因子其中σ=[0,1]为中心参数;
9)用牛顿法求解上述非线性方程组,计算雅可比矩阵、海森矩阵和常数项等,求解修正方程组,得到各原变量和对偶变量的修正量,乘上步长对变量进行修正,直到对偶间隙Gap小于收敛精度(ε=10-6),否则不收敛;
求解修正方程得到修正量,计算原变量和对偶变量步长:
按下式修正变量:
10)输出优化调度结果,包括:电力网络发电机有功、无功出力和各节点压力,天然气网络气源点天然气输出和各节点压力值,能源中心各调度因子。
有益效果:本发明相对于现有技术而言:本发明首先建立了电力网络、天然气网络和能源中心模型,电力网络和天然气网络通过能源中心耦合形成电气混联***。然后以总的能源成本为目标函数、考虑了各种约束条件建立电气混联***优化调度数学模型。用原-对偶内点法求解,求解过程中先后引入松弛变量和障碍参数将模型变为只有等式约束的模型,再引入拉格朗日乘子得到拉格朗日函数,之后用牛顿法求解其KKT条件形成的非线性方程组。构造的算例仿真结果表明本发明对电气混联***的优化效果优于单独优化效果。
附图说明
图1为本发明方法流程图;
图2为能源中心模型示意图;
图3为能源中心示例示意图;
图4为算例结果分析的不同能源中心输入结果对比图。
具体实施方案
下面结合附图对发明的技术流程进行详细说明:
电力网络模型实施例
电力网络的拓扑结构可以用导纳矩阵Y描述。电力网络的状态变量有节点电压、节点功率及支路功率等,这里用电压的直角坐标形式给出节点有功功率和无功功率:
其中:Pi,Qi分别为电力网络中节点i有功功率和无功功率;ei,fi分别为节点i电压向量的实部和虚部;Gij,Bij分别为节点导纳矩阵第i行第j列元素的实部和虚部。
天然气网络模型实施例
典型的天然气网络包括一个或多个气源点(天然气生产点和存储点),一个或多个负荷(发电厂、其他的网络或存储点),管道,压缩机和其他设备。用于天然气网络建模的基本要素有三个:1)管道;2)压缩机;3)注入点(节点或母线)。天然气模型两个最主要的部分是管道流量模型和压缩机模型。
天然气管道流量方程描述了天然气流量值与管道两端压力值以及相关的气体性质、管道特性和运行状态之间的关系。形成流量方程需要考虑很多因素,没有一个流量方程适用于所有的情况。理想条件下,管道k从节点i到节点j的流量值可以用下面的方程表示:
其中:fkij为管道流量值; Fk为管道摩擦系数;Dk为管道内径;G为气体重力系数;Lk为管道长度;πi为节点i压力值;πj为节点j压力值;π0为标准压力值;T0为标准温度值;Tka为平均气体温度;Za为平均气体压缩系数。
对于高压网络的全紊流形式,流量方程可以进一步简化为:
其中: ε为管道效率。
压缩机安装在气体网络中为了传输气体和补偿由于摩擦阻力等造成的下游压力损失。压缩机的运行需要消耗大量的能量,在大型天然气网络中,驱动主要压缩机最经济的能量源是通过压缩机的天然气。对于压缩机的网络分析关键的方程是能量消耗,它是流过压缩机的气体流量和注入、流出压缩机的气体压力的函数。在理想气体条件下,加上经验的量化,绝热压缩机的能量消耗方程可以表示为:
其中:fCk为通过压缩机的气体流量;πi为气体注入压缩机压力;πj为气体输出压缩机压力;Zki为压缩机入口的气体压缩系数;Tki为压缩机输出温度;α为热量系数;ηk为压缩机效率。
转化为消耗的流量值:
根据基尔霍夫第一定律:所有流入与流出某网络节点的流量代数和为零。故每一个节点的流量平衡方程可以用下面的矩阵形式表示:
(A+U)f+w-Tτ=0;
其中:A、U和T分别为支路、压缩机和压缩机消耗与节点有关的方向矩阵,f为支路流量值向量;w为各节点的气体注入向量;τ为各压缩机消耗流量值向量。
能源中心模型实施例
能源中心被定义为由能源转化设备和储能设备构成、能够实现多种能源相互转化和存储的虚拟实体。能源中心之间通过能源传输设备相互连接,构成多能源网络***。在能源中心内部,能源可能被负荷消耗,或者被转化为其他形式。
能源中心有多种能源形式α,β,…∈ε,每一种能源形式可以是能源中心的输入或输出。这些输入能量(输出能量)被定义为Pα,Pβ,…,Pω(Lα,Lβ,…Lω)。如图2所示。对于多输入多输出的能源中心可以用下面的矩阵方程描述:
其中:P和L分别为输入输出向量;矩阵C被称为耦合矩阵。耦合矩阵从数学上描绘了能量从能源中心输入到输出的分配。耦合矩阵中的元素为耦合因子。每一个耦合因子表示了一个特定的输入到一个特定的输出。
图3给出了一种能源中心的示例。能源中心输入由电力和天然气提供,能源在能源中心内部通过变压器、燃气轮机和燃气锅炉转化,输出供应电负荷和热负荷。其中ν为调度因子。据此可以列出能源中心输入输出平衡方程:
Le=ηTPe+νηGTePg;
Lh=νηGThPg+(1-ν)ηFPg;
写成矩阵形式:
其中:ηT表示变压器效率;ηGTe表示燃气轮机气转电的效率;ηGTh表示燃气轮机气转热的效率;ηF表示燃气锅炉的效率。
电气混联***优化调度模型实施例
根据上述模型,电力网络与天然气网络通过能源中心耦合,形成电气混联***。以总的能源成本为目标函数,考虑电力网络、天然气网络及能源中心的各种约束,建立电气混联***优化调度模型:
minf(x);
s.t.h(x)=0;
gmin≤g(x)≤gmax;
其中:h表示等式约束;g表示不等式约束;gmax,gmin分别表示不等式约束的上、下限。
进行选取变量,控制变量:
1)电力网络控制变量:发电机有功出力PG;发电机无功出力QR;
2)天然气网络控制变量:气源点天然气供应量NG;
3)能源中心控制变量:调度因子ν。
则本发明选取的控制变量为:
u=[PG1,...,PGng,QG1,...,QGng,NG1,…NGnGS,ν1,…,νnEH]T;
其中:ng为电力网络发电机台数;nGS为天然气网络气源个数;nEH为能源中心个数。
状态变量:
1)电力网络状态变量:节点电压向量的实部e和虚部f;
2)天然气网络状态变量:节点的压力π和通过压缩机的流量fC;
3)能源中心状态变量:电能的输入量Pe和天然气的输入量Pg。
则本发明选取的状态变量为:
其中:nb为电力网络的节点个数;nGB为天然气网络节点个数;nGC为天然气网络中压缩机的个数。
总的变量定义为:
本发明以***总的能源成本作为目标函数:
其中:a2i,a1i,a0i为能量损耗曲线参数;Pi为广义能量,这里包括电力网中发电机发出的有功功率和天然气网络气源点提供气体的能量。
等式约束
1)电力网络等式约束:
ΔPi=PGi-PDi-Pei-Pi=0;
ΔQi=QRi-QDi-Qi=0;
其中:ΔPi,ΔQi为各节点有功、无功功率不平衡量;Δangle为平衡节点角度约束;tanθbal为平衡节点的角度正切值;PGi,QRi分别为发电机i的有功、无功出力;PDi,QDi分别为节点i的有功、无功负荷;Pei为节点i对能源中心的电能输入。
2)天然气网络等式约束:
ΔWi=NGi+wi-Fi-Pgi=0;
其中:ΔWi为天然气网络中各节点流量值的不平衡量;NGi为气源点对节点i的气体注入量;wi,Fi分别为节点i的气体注入量和负荷量;Pgi为节点i对能源中心的天然气输入;ΔPCk为天然气网络中压缩机k状态量的不平衡量;Rk分别为压缩机k状态量的给定值和计算值。
3)能源中心等式约束:
ΔLei=Lei-ηTPei-νiηGTePgi;
ΔLhi=Lhi-νiηGThPgi-(1-νi)ηFPgi;
其中:ΔLei,ΔLhi分别为能源中心i的电能负荷与热负荷的不平衡量;Lei,Lhi分别为能源中心i的电能负荷与热负荷;Pei,Pgi分别为能源中心i的电能与天然气输入;νi为能源中心i的调度因子。
不等式约束
1)电力网络不等式约束:
PGimin≤PGi≤PGimax;
QRimin≤QRi≤QRimax;
其中:PGimin,PGimax为发电机所发出有功功率的下限和上限;QRimin,QRimax为发电机所发无功功率的下限和上限;为节点电压幅值平方的下限和上限。
2)天然气网络不等式约束:
0≤NGi≤NGimax;
πimin≤πi≤πimax;
其中:NGimax为天然气网络中各气源点气体供应量的上限;πimin与πimax分别为各节点压力值的下限和上限。
3)能源中心不等式约束:
0≤νi≤1;
通过原-对偶内点法求解,首先引入广义松弛变量l=[l1,…,lr]T,u=[u1,…,ur]T将广义不等式转化为等式:
g(x)+u-gmax=0;
g(x)-l-gmin=0;
应满足:u>0,l>0。
然后,为了保证目标函数在可行域内取得最小值,对目标函数引入对数障碍参数μ,当u,l靠近约束边界时,目标函数趋于无穷大。此时优化调度模型为只有等式约束的模型:
s.t.h(x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
其中:扰动因子(障碍参数)μ>0。
这时可以直接通过拉格朗日乘子法求解。引入拉格朗日乘子y=[y1,...,ym]T,z=[z1,…,zr]T,w=[w1,...,wm]T得到拉格朗日函数:
根据拉格朗日函数极值存在的必要条件(KKT条件)可以得到非线性方程组:
其中:e为各元素为1的r维列向量;L=diag(l);U=diag(u);Z=diag(z);W=diag(w)均为对角矩阵。
计算可以得到:
定义互补间隙(对偶间隙)为:Gap=lTz-uTw。实践表明,目标函数中的扰动因子μ按下式计算收敛性较好:
其中:σ=[0,1]为中心参数。σ=1时,有利于提高算法的可行性,但对减少对偶间隙无作用;σ=0减少了算法的对偶间隙,有利于解的最优性。在保证可行性的基础上尽量提高解的最优性,σ一般取0.1。
可以用牛顿-拉弗逊法求解上述非线性方程组,其一阶修正方程组的矩阵形式为:
其中: Lx,Ly,Lz,Lw,为KKT方程上次迭代的残差;和分别为h(x)和g(x)的海森矩阵。
求解修正方程得到修正量,计算原变量和对偶变量步长:
按下式修正变量:
依次迭代,直到对偶间隙Gap小于收敛精度。
实施例
为验证本发明效果,将Matpower中IEEE14节点电力网络与某14节点天然气网络通过图3中的能源中心耦合构造一个算例。其中电力网络与天然气网络通过5个能源中心联系,将天然气网络中原来的5个负荷点连接到能源中心,作为能源中心的天然气输入;选取电力网络中负荷比较重的5个负荷点与能源中心连接,作为能源中心的电能输入,能源中心负荷侧考虑电能负荷和热能输出,并同时考虑能源中心调度因子的影响。各能源中心与电力网络及天然气网络具体连接情况和能源中心输出情况见表1所示。
用Matlab对构造的算例进行编程并仿真,对仿真结果进行分析。
当将天然气网络中各节点压力值限制在500~1500psia内,取天然气损耗曲线系数a2=0,a1=2,a0=0,能源中心输出取表1中数据时,混联***综合优化的最小成本为$5338.5239。电力网单独优化的最小成本为$8081.526,天然气网络单独优化的最小成本为a1·NG·GHV/baseMVA=$89.98,故当电力网与天然气单独优化的总成本为$8081.526+$89.98=$8171.506。可见混联***的综合优化成本明显低于单独优化成本,说明了电力网络与天然气网混联***的综合优化具有显著的经济效益。
为了更充分说明混联***综合优化成本低于单独优化成本,选取与上面相同的压力范围和天然气损耗曲线系数,对不同的能源中心电能输出值与热能输出值仿真,优化结果如图4所示。从对比图可知,对于不同的能源中心输出,混联***综合优化成本均低于独立优化的成本,进一步说明了混联***综合优化的优势。
表1混联***节点连接与能源中心输出
Claims (3)
1.一种基于能源中心的电气混联***建模和优化调度方法,其特征在于:包括以下步骤:
1)获取电力网络的参数信息,包括:输电线路的首端节点和末端节点编号,支路π型等效电路的电阻、电抗,对地并联电导、电纳,变压器变比和阻抗,各节点负荷以及发电机输出有功、无功约束,各节点电压约束;
2)获取天然气网络的参数信息,包括:输气管道的首端节点和末端节点编号,管道的长度、内径和传输效率等物理特性,压缩机的状态量,各节点气负荷以及气源点天然气输出约束,各节点压力约束;
3)多输入多输出能源中心模型可以用矩阵方程描述:
其中:P和L分别为输入输出向量;矩阵C被称为耦合矩阵;
列出能源中心输入输出平衡方程:
其中:ηT表示变压器效率;ηGTe表示燃气轮机气转电的效率;ηGTh表示燃气轮机气转热的效率;ηF表示燃气锅炉的效率;
4)建立电气混联***优化调度模型:
minf(x);
s.t.h(x)=0;
gmin≤g(x)≤gmax;
其中:h表示等式约束;g表示不等式约束;gmax,gmin分别表示不等式约束的上、下限;
5)引入松弛变量和障碍参数将模型转化成只有等式约束的模型:
s.t.h(x)=0;
g(x)+u-gmax=0;
g(x)-l-gmin=0;
其中:扰动因子(障碍参数)μ>0;
6)引入拉格朗日乘子得到拉格朗日函数:
7)根据拉格朗日极值存在条件(KKT条件)得到一组非线性方程组,用牛顿-拉弗逊法求解,其一阶修正方程为:
其中: Lx,Ly,Lz,Lw,为KKT方程上次迭代的残差;和分别为h(x)和g(x)的海森矩阵;
8)求解修正方程得到修正量,计算原变量和对偶变量步长:
9)按下式修正变量:
10)按上述步骤迭代,直到对偶间隙小于收敛精度,得到电气混联***优化调度的最优解。
2.根据权利要求1所述的基于能源中心的电气混联***建模和优化调度方法,其特征在于:所述步骤1中,对于电力网络中节点i:
其中:Pi,Qi分别为电力网络中节点i有功功率和无功功率;ei,fi分别为节点i电压向量的实部和虚部;Gij,Bij分别为节点导纳矩阵第i行第j列元素的实部和虚部。
3.根据权利要求1所述的基于能源中心的电气混联***建模和优化调度方法,其特征在于:所述步骤2中,引入包括管道流量方程、压缩机流量消耗方程及流量平衡方程;
理想条件下,管道k从节点i到节点j的流量值可以用下面的方程表示:
其中:fkij为管道流量值; Fk为管道摩擦系数;Dk为管道内径;G为气体重力系数;Lk为管道长度;πi为节点i压力值;πj为节点j压力值;π0为标准压力值;T0为标准温度值;Tka为平均气体温度;Za为平均气体压缩系数;
对于高压网络的全紊流形式,流量方程可以进一步简化为:
其中: ε为管道效率;
理想气体条件下,压缩机的能量消耗方程可以表示为:
其中:fCk为通过压缩机的气体流量;πi为气体注入压缩机压力;πj为气体输出压缩机压力;Zki为压缩机入口的气体压缩系数;Tki为压缩机输出温度;α为热量系数;ηk为压缩机效率;
转化为消耗的流量值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510498181.XA CN105046369B (zh) | 2015-08-13 | 2015-08-13 | 一种基于能源中心的电气混联***建模和优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510498181.XA CN105046369B (zh) | 2015-08-13 | 2015-08-13 | 一种基于能源中心的电气混联***建模和优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105046369A true CN105046369A (zh) | 2015-11-11 |
CN105046369B CN105046369B (zh) | 2019-04-30 |
Family
ID=54452896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510498181.XA Active CN105046369B (zh) | 2015-08-13 | 2015-08-13 | 一种基于能源中心的电气混联***建模和优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105046369B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958531A (zh) * | 2016-05-27 | 2016-09-21 | 清华大学 | 一种电-气耦合多能流网络状态估计方法 |
CN106056478A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合***中热网的区间潮流计算方法 |
CN106056251A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合多能流***的优化调度方法 |
CN106096777A (zh) * | 2016-06-12 | 2016-11-09 | 清华大学 | 一种电‑气耦合多能流***的优化调度方法 |
CN106447152A (zh) * | 2016-06-21 | 2017-02-22 | 中国南方电网有限责任公司电网技术研究中心 | 基于能源中心的电‑气‑热***协同调度方法和*** |
CN107292456A (zh) * | 2017-08-01 | 2017-10-24 | 重庆大学 | 基于交替方向乘子法的电‑气能量流分布式协同优化计算方法 |
CN107394828A (zh) * | 2017-08-24 | 2017-11-24 | 河海大学 | 一种基于概率潮流的电气互联综合能源***优化运行分析方法 |
CN107404118A (zh) * | 2017-09-06 | 2017-11-28 | 河海大学 | 基于随机响应面法的电气互联***概率最优潮流计算方法 |
CN108062599A (zh) * | 2017-12-12 | 2018-05-22 | 西安交通大学 | 一种电力***和天然气***协同运行的建模方法 |
CN108258679A (zh) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | 考虑供热网储热特性的电-热综合能源***优化调度方法 |
CN108510299A (zh) * | 2017-02-28 | 2018-09-07 | 腾讯科技(北京)有限公司 | 确定多媒体资源投放策略的方法及装置 |
CN108596525A (zh) * | 2018-06-29 | 2018-09-28 | 国家电网有限公司 | 冷热电气多能互补的微能源网鲁棒优化调度方法 |
CN108846507A (zh) * | 2018-05-29 | 2018-11-20 | 西安交通大学 | 基于混合整数二阶锥规划的电-气耦合***日前经济调度方法 |
CN109086474A (zh) * | 2018-06-12 | 2018-12-25 | 河海大学 | 基于二阶锥优化算法的天然气***最优潮流求解方法 |
CN109190849A (zh) * | 2018-10-18 | 2019-01-11 | 清华-伯克利深圳学院筹备办公室 | 一种能源价格的预测方法、装置、定价设备及存储介质 |
CN110428185A (zh) * | 2019-08-08 | 2019-11-08 | 河海大学 | 基于伪量测模型的电-热互联综合能源***抗差状态估计方法 |
CN112039750A (zh) * | 2020-09-01 | 2020-12-04 | 四川大学 | 一种基于综合能源***的电气联合***网络状态预估方法 |
CN112069638A (zh) * | 2020-09-08 | 2020-12-11 | 广东电网有限责任公司电力科学研究院 | 一种低压配网储能设备接入点选址方法及相关装置 |
CN112818492A (zh) * | 2021-02-02 | 2021-05-18 | 山东大学 | 计及能源集站状态变量的电-气耦合网络能量流求解方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (zh) * | 2011-12-13 | 2012-06-27 | 河海大学 | 基于原对偶内点法的vsc-hvdc交直流最优潮流方法 |
CN103366097A (zh) * | 2013-07-24 | 2013-10-23 | 国家电网公司 | 一种基于类扩展变量内点法的最优潮流计算方法 |
CN103793755A (zh) * | 2014-01-14 | 2014-05-14 | 河海大学 | 基于原对偶内点法的混合直流输电***最优潮流方法 |
CN104392285A (zh) * | 2014-11-28 | 2015-03-04 | 河海大学 | 一种含混合直流输电的电力***最优潮流获取方法 |
CN104393592A (zh) * | 2014-11-27 | 2015-03-04 | 河海大学 | 一种计及温度影响的电力***最优潮流算法 |
-
2015
- 2015-08-13 CN CN201510498181.XA patent/CN105046369B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102522746A (zh) * | 2011-12-13 | 2012-06-27 | 河海大学 | 基于原对偶内点法的vsc-hvdc交直流最优潮流方法 |
CN103366097A (zh) * | 2013-07-24 | 2013-10-23 | 国家电网公司 | 一种基于类扩展变量内点法的最优潮流计算方法 |
CN103793755A (zh) * | 2014-01-14 | 2014-05-14 | 河海大学 | 基于原对偶内点法的混合直流输电***最优潮流方法 |
CN104393592A (zh) * | 2014-11-27 | 2015-03-04 | 河海大学 | 一种计及温度影响的电力***最优潮流算法 |
CN104392285A (zh) * | 2014-11-28 | 2015-03-04 | 河海大学 | 一种含混合直流输电的电力***最优潮流获取方法 |
Non-Patent Citations (1)
Title |
---|
孙小璞: "大规模风电场并网的优化调度建模及算法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105958531A (zh) * | 2016-05-27 | 2016-09-21 | 清华大学 | 一种电-气耦合多能流网络状态估计方法 |
CN106056478A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合***中热网的区间潮流计算方法 |
CN106056251A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合多能流***的优化调度方法 |
CN106096777A (zh) * | 2016-06-12 | 2016-11-09 | 清华大学 | 一种电‑气耦合多能流***的优化调度方法 |
CN106056251B (zh) * | 2016-06-12 | 2019-06-18 | 清华大学 | 一种电-热耦合多能流***的优化调度方法 |
CN106096777B (zh) * | 2016-06-12 | 2019-05-21 | 清华大学 | 一种电-气耦合多能流***的优化调度方法 |
CN106056478B (zh) * | 2016-06-12 | 2019-05-21 | 清华大学 | 一种电-热耦合***中热网的区间潮流计算方法 |
CN106447152A (zh) * | 2016-06-21 | 2017-02-22 | 中国南方电网有限责任公司电网技术研究中心 | 基于能源中心的电‑气‑热***协同调度方法和*** |
CN108510299A (zh) * | 2017-02-28 | 2018-09-07 | 腾讯科技(北京)有限公司 | 确定多媒体资源投放策略的方法及装置 |
CN107292456B (zh) * | 2017-08-01 | 2020-11-24 | 重庆大学 | 基于交替方向乘子法的电-气能量流分布式协同优化计算方法 |
CN107292456A (zh) * | 2017-08-01 | 2017-10-24 | 重庆大学 | 基于交替方向乘子法的电‑气能量流分布式协同优化计算方法 |
CN107394828A (zh) * | 2017-08-24 | 2017-11-24 | 河海大学 | 一种基于概率潮流的电气互联综合能源***优化运行分析方法 |
CN107404118A (zh) * | 2017-09-06 | 2017-11-28 | 河海大学 | 基于随机响应面法的电气互联***概率最优潮流计算方法 |
CN108062599A (zh) * | 2017-12-12 | 2018-05-22 | 西安交通大学 | 一种电力***和天然气***协同运行的建模方法 |
CN108062599B (zh) * | 2017-12-12 | 2020-08-18 | 西安交通大学 | 一种电力***和天然气***协同运行的建模方法 |
CN108258679A (zh) * | 2017-12-25 | 2018-07-06 | 国网浙江省电力有限公司经济技术研究院 | 考虑供热网储热特性的电-热综合能源***优化调度方法 |
CN108258679B (zh) * | 2017-12-25 | 2020-05-12 | 国网浙江省电力有限公司经济技术研究院 | 考虑供热网储热特性的电-热综合能源***优化调度方法 |
CN108846507B (zh) * | 2018-05-29 | 2020-11-10 | 西安交通大学 | 基于混合整数二阶锥规划的电-气耦合***日前经济调度方法 |
CN108846507A (zh) * | 2018-05-29 | 2018-11-20 | 西安交通大学 | 基于混合整数二阶锥规划的电-气耦合***日前经济调度方法 |
CN109086474A (zh) * | 2018-06-12 | 2018-12-25 | 河海大学 | 基于二阶锥优化算法的天然气***最优潮流求解方法 |
CN108596525A (zh) * | 2018-06-29 | 2018-09-28 | 国家电网有限公司 | 冷热电气多能互补的微能源网鲁棒优化调度方法 |
CN108596525B (zh) * | 2018-06-29 | 2022-07-22 | 国家电网有限公司 | 冷热电气多能互补的微能源网鲁棒优化调度方法 |
CN109190849A (zh) * | 2018-10-18 | 2019-01-11 | 清华-伯克利深圳学院筹备办公室 | 一种能源价格的预测方法、装置、定价设备及存储介质 |
CN109190849B (zh) * | 2018-10-18 | 2022-05-31 | 清华-伯克利深圳学院筹备办公室 | 一种能源价格的预测方法、装置、定价设备及存储介质 |
CN110428185A (zh) * | 2019-08-08 | 2019-11-08 | 河海大学 | 基于伪量测模型的电-热互联综合能源***抗差状态估计方法 |
CN112039750A (zh) * | 2020-09-01 | 2020-12-04 | 四川大学 | 一种基于综合能源***的电气联合***网络状态预估方法 |
CN112069638A (zh) * | 2020-09-08 | 2020-12-11 | 广东电网有限责任公司电力科学研究院 | 一种低压配网储能设备接入点选址方法及相关装置 |
CN112818492A (zh) * | 2021-02-02 | 2021-05-18 | 山东大学 | 计及能源集站状态变量的电-气耦合网络能量流求解方法 |
CN112818492B (zh) * | 2021-02-02 | 2022-12-02 | 山东大学 | 计及能源集站状态变量的电-气耦合网络能量流求解方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105046369B (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105046369A (zh) | 一种基于能源中心的电气混联***建模和优化调度方法 | |
Lei et al. | Multi-objective stochastic expansion planning based on multi-dimensional correlation scenario generation method for regional integrated energy system integrated renewable energy | |
CN110866627B (zh) | 一种考虑阶梯气价的多区域电-气耦合综合能源***优化调度方法 | |
Jiang et al. | Multi-stage progressive optimality algorithm and its application in energy storage operation chart optimization of cascade reservoirs | |
CN109031952B (zh) | 一种电-气互联综合能源***混合控制方法 | |
CN104734155B (zh) | 一种获取电‑气互联能源***可用输电能力的方法 | |
CN108667007A (zh) | 计及电-气耦合***约束的电压稳定裕度计算方法 | |
CN105576657A (zh) | 能源互联下主动配网与输电网集成***的随机模糊连续潮流建模及算法 | |
CN107732982A (zh) | 考虑模型预测控制的综合能源***多时间尺度调度方法 | |
CN106655226A (zh) | 基于智能软开关的有源配电网不对称运行优化方法 | |
CN106655177A (zh) | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 | |
CN108023364A (zh) | 基于凸差规划的配电网分布式电源最大接入能力计算方法 | |
CN106712035A (zh) | 一种电力***经济调度方法 | |
CN106532710A (zh) | 计及电压稳定约束的微电网潮流优化方法 | |
CN105243432A (zh) | 基于nsga-ⅲ的upfc选址和容量多目标配置方法 | |
CN105046588A (zh) | 一种基于网损迭代的改进直流动态最优潮流的计算方法 | |
CN115995790A (zh) | 一种配电网故障恢复方法、***、设备及介质 | |
CN112084629A (zh) | 一种基于两阶段鲁棒优化的多能虚拟电厂聚合方法 | |
CN106505624A (zh) | 确定配电网分布式电源最优吸纳能力的调控***及方法 | |
CN107123983A (zh) | 一种基于安全域的变电站接入方案辅助评估方法 | |
Shengyu et al. | Review on interdependency modeling of integrated energy system | |
CN113901653A (zh) | 一种电-气综合能源***及其潮流求解方法 | |
An et al. | Optimal operation of multi-micro energy grids under distribution network in Southwest China | |
CN104732302A (zh) | 一种基于免疫算法的多级电磁环网优化解环方法 | |
CN105244874B (zh) | 一种计及故障风险及成本的配电网无功规划优化模型 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |