CN113381453A - 交流电力系统电源发电功率调度的节点级分散方法 - Google Patents
交流电力系统电源发电功率调度的节点级分散方法 Download PDFInfo
- Publication number
- CN113381453A CN113381453A CN202110599963.8A CN202110599963A CN113381453A CN 113381453 A CN113381453 A CN 113381453A CN 202110599963 A CN202110599963 A CN 202110599963A CN 113381453 A CN113381453 A CN 113381453A
- Authority
- CN
- China
- Prior art keywords
- node
- power
- branch
- power supply
- alternating current
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000010248 power generation Methods 0.000 title claims abstract description 30
- 239000006185 dispersion Substances 0.000 title claims abstract description 24
- 230000006870 function Effects 0.000 claims description 31
- 238000004590 computer program Methods 0.000 claims description 14
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Power Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种交流电力系统电源发电功率调度的节点级分散方法,属于电力工程领域,首先根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程;根据线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型;根据二次规划模型建立拉格朗日函数;根据拉格朗日函数建立节点级分散迭代公式,继而根据节点级分散迭代公式获取交流电力系统的电源发电功率调度值。这使交流电力系统电源发电功率调度的求解结果唯一和全局最优,避免了传统方法中电源发电功率调度之解的全局最优性无保障的缺陷;同时,这种交流电力系统的电源发电功率调度的求解既节点级分散又无源荷的功率私密信息泄露。
Description
技术领域
本申请涉及电力工程领域,尤其涉及一种交流电力系统电源发电功率调度的节点级分散方法。
背景技术
交流电力系统电源发电功率调度是保障其运行高经济性的关键环节。目前,它要么集中、要么分片构建和求解以非线性节点功率平衡方程组为约束的优化模型实现。但约束的非线性导致调度之解的全局最优性无保障缺陷、求解的集中或分片性导致节点级源荷的功率私密数据泄露的缺陷。同时,已有调度方法因难以计入网损而存在劣化调度结果经济性的缺陷。
发明内容
本申请实施例提供一种交流电力系统电源发电功率调度的节点级分散方法,能够解决传统的交流电力系统电源发电功率调度方法所存在的调度之解的全局最优性无保障、源荷的功率私密信息泄露以及难以计入网损的问题。
本申请实施例的第一方面提供了一种交流电力系统电源发电功率调度的节点级分散方法,包括:
根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程;
根据所述线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型;
根据所述二次规划模型建立拉格朗日函数;
根据所述拉格朗日函数建立节点级分散迭代公式,继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值。
本申请实施例的第二方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述交流电力系统电源发电功率调度的节点级分散方法的步骤。
本申请实施例的第三方面提供了一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述交流电力系统电源发电功率调度的节点级分散方法的步骤。
本申请实施例与现有技术相比存在的有益效果是:由于采用线性渐近方程建立交流电力系统电源发电功率调度的二次规划模型,故交流电力系统电源发电功率调度的求解结果唯一和全局最优,避免了传统方法中调度之解的全局最优性无保障的缺陷;又由于根据电源的有功供电成本以及无功消耗量建立交流电力系统电源发电功率调度的二次规划模型,故解决了网损难以计入的问题;同时,由于建立节点级分散迭代公式,故交流电力系统的电源发电功率调度的求解既节点级分散又无源荷的功率私密信息泄露。
附图说明
为了更清楚地说明本发明实施例技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种交流电力系统电源发电功率调度的节点级分散方法的实现流程图;
图2是本发明实施例提供的交流电力系统通用模型的结构示意图;
图3是本发明实施例提供的一种终端设备的结构示意图。
具体实施方式
为了使本申请所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
参见图1,图1是本发明实施例提供的一种交流电力系统电源发电功率调度的节点级分散方法的实现流程图。如图所示的交流电力系统电源发电功率调度的节点级分散方法可包括以下步骤:
在步骤101中,根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程。
具体实施中,步骤101可以包括步骤A1和步骤B1。
在步骤A1中,根据交流电力系统的支路导纳参数、支路两端的电压幅值以及支路两端的电压相角,运用电功率定义式并撤分其中的耦合项,建立如下的支路传输功率的线性渐近表达式:
Pij=αijVi+βijVj+γijθi+δijθj
其中,Pij为支路ij传输的有功功率;Qij为支路ij传输的无功功率;i和j均为交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为交流电力系统中节点的总个数;gij为支路ij的电导参数;bij为支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为节点j的电压相角变量;ξij是按照ξij=sinθij-θi+θj确定的支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的支路ij的第4修正导纳;是按照 确定的支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的三十三支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的支路ij的第8修正导纳;Vi和Vj均为标幺值电压。gij和bij均为已知的电力系统参数。
通过将非线性的电功率的定义式变换为线性渐近表达式,避免了以非线性方程为约束的优化规划模型难以求解的问题。
在步骤B1中,根据支路传输功率的线性渐近表达式和交流电力系统的支路连接结构,按Kirchhoff电流定律建立如下的节点i的功率平衡的线性渐近方程:
其中,PGi为接于节点i的电源的有功功率变量,且为未知量;QGi为接于节点i的电源的无功功率变量,且为未知量;PDi为接于节点i的负荷的有功功率参数,且为已知量;QDi为接于节点i的负荷的无功功率参数,且为已知量。
上述节点功率平衡的线性渐近方程是关于节点电压幅值和相角的线性方程,且节点功率平衡的线性渐近方程随节点电压幅值和相角逼近真值而逼近按照电功率定义和Kirchhoff电流定律得到的精确的节点功率平衡方程。这正是称上述线性渐近方程为节点功率平衡的线性渐近方程的缘故。
在步骤102中,根据线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型。
步骤102包括:以线性渐近方程为约束、以电源的有功供电成本与电源的无功消耗量的平方和之和最小为目标函数,建立如下的交流电力系统电源发电功率调度的二次规划模型:
其中,aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点;i和j均为交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为交流电力系统中节点的总个数;gij为支路ij的电导参数;bij为支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为节点j的电压相角变量;Vn是编号为n的节点的电压幅值,且Vn=1是已知量;θn是编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的支路ij的第4修正导纳;是按照 确定的支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的所述支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的支路ij的第8修正导纳;Vi和Vj均为标幺值电压。PGi为接于节点i的电源的有功功率变量,且为未知量;QGi为接于节点i的电源的无功功率变量,且为未知量;PDi为接于节点i的负荷的有功功率参数;QDi为接于节点i的负荷的无功功率参数。gij、bij、QDi和PDi都为已知的电力系统参数。
通过上述二次规划模型,在满足节点功率平衡的线性渐近方程约束下,实现发电成本最小。
上述二次规划模型中目标函数的二次项系数都大于零,故是凸函数,约束条件是线性等式,因此是凸二次规划。按最优化理论,它的局部最优解唯一且是全局最优解。因此,二次规划模型的拉格朗日函数的驻点就是唯一的全局最优解。
在步骤103中,根据二次规划模型建立拉格朗日函数。
步骤103包括:根据二次规划模型,按拉格朗日函数的定义建立如下的拉格朗日函数。
其中,是拉格朗日函数;λi为对应节点i的有功功率平衡方程的拉格朗日乘子;ξi为对应节点i的无功功率平衡方程的拉格朗日乘子;aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点。i和j均为交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为交流电力系统中节点的总个数;gij为支路ij的电导参数;bij为支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为节点j的电压相角变量;Vn是编号为n的节点的电压幅值,且Vn=1是已知量;θn是编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的支路ij的第4修正导纳;是按照确定的支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的支路ij的第8修正导纳;Vi和Vj均为标幺值电压。PGi为接于节点i的电源的有功功率变量,且为未知量;QGi为接于节点i的电源的无功功率变量,且为未知量;PDi为接于节点i的负荷的有功功率参数;QDi为接于节点i的负荷的无功功率参数。gij、bij、QDi和PDi都为已知的电力系统参数。
在步骤104中,根据拉格朗日函数建立节点级分散迭代公式,继而根据节点级分散迭代公式获取交流电力系统的电源发电功率调度值。
具体实施中,步骤104可以包括步骤A2和步骤B2。
在步骤A2中,根据拉格朗日函数,按驻点的定义建立如下的驻点方程组:
其中,是拉格朗日函数;λi为对应节点i的有功功率平衡方程的拉格朗日乘子;ξi为对应节点i的无功功率平衡方程的拉格朗日乘子;aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点。i和j均为交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为交流电力系统中节点的总个数;gij为支路ij的电导参数;bij为支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为节点j的电压相角变量;Vn是编号为n的节点的电压幅值,且Vn=1是已知量;θn是编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的支路ij的第4修正导纳;是按照确定的支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的三十三支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的支路ij的第8修正导纳;Vi和Vj均为标幺值电压。PGi为接于节点i的电源的有功功率变量,且为未知量;QGi为接于节点i的电源的无功功率变量,且为未知量;PDi为接于节点i的负荷的有功功率参数;QDi为接于节点i的负荷的无功功率参数。gij、bij、QDi和PDi都为已知的电力系统参数。
通过求驻点方程组的解,从而获取到目标函数取极小值时各个变量的值。
在步骤B2中,基于驻点方程组,建立如下的节点级分散迭代公式,继而根据节点级分散迭代公式获取交流电力系统的电源发电功率调度值:
其中,(t+1)表示第t+1步的迭代结果;(t)表示第t步的迭代结果;σ为大于0小于1的惯性参数;Ωi是编号为i的节点的所有邻居节点之编号集合。
按上述节点级分散迭代公式迭代计算,直至收敛,所得各节点的电源的有功功率和无功功率最终解,就是交流电力系统的电源发电功率调度结果。从而实现交流电力系统电源发电功率的节点级分散调度。
步骤B2按照控制理论,将连续方程组(驻点方程组)转化为离散迭代表达式(节点级分散迭代公式)。按上述节点级分散迭代公式计算编号为i的节点的θi、Vi、λi、ξi、PGi和QGi时,只需要编号属于集合Ωi的节点(也就是只需要邻居节点)的电压幅值和相角及拉格朗日乘子,不需要邻居节点的源荷功率私密数据。因此,上述迭代公式是节点级分散的,且邻居节点的源荷功率私密信息无泄露。这正是称本发明给出的方法为交流电力系统电源发电功率调度的节点级分散方法的缘故。
本申请实施例首先根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程;根据线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型;根据二次规划模型建立拉格朗日函数;根据拉格朗日函数建立节点级分散迭代公式,继而根据节点级分散迭代公式获取交流电力系统的电源发电功率调度值。由于采用线性渐近方程建立交流电力系统电源发电功率调度的二次规划模型,故交流电力系统电源发电功率调度的求解结果唯一和全局最优,避免了传统方法中调度之解的全局最优性无保障的缺陷;又由于根据电源的有功供电成本以及无功消耗量建立交流电力系统电源发电功率调度的二次规划模型,故解决了网损难以计入的问题;同时,由于建立节点级分散迭代公式,故交流电力系统的电源发电功率调度的求解既节点级分散又无源荷的功率私密信息泄露。
本申请实施例第二方面提供的是一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述交流电力系统电源发电功率调度的节点级分散方法的步骤。
图3是本申请实施例第三方面提供的终端设备的示意图。该实施例的终端设备3包括:处理器30、存储器31以及存储在所述存储器31中并可在所述处理器30上运行的计算机程序32,所述处理器30执行所述计算机程序32时实现上述交流电力系统电源发电功率调度的节点级分散方法实施例中的各个步骤,例如图1所示的步骤101至104。本领域技术人员应当理解,图3仅仅是终端设备3的示例,并不构成对终端设备3的限定。所述终端设备3包括但不限于处理器30、存储器31以及存储在所述存储器31中并可在所述处理器30上运行的计算机程序32,例如所述终端设备是本身存储器上或外接可移动存储器上存储有所述计算机程序32的服务器、计算机、掌上电脑及其与输入输出设备和网络接入设备的组合。
以上所述实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围,均应包含在本申请的保护范围之内。
Claims (7)
1.一种交流电力系统电源发电功率调度的节点级分散方法,其特征在于,所述交流电力系统电源发电功率调度的节点级分散方法包括:
根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程;
根据所述线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型;
根据所述二次规划模型建立拉格朗日函数;
根据所述拉格朗日函数建立节点级分散迭代公式,继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值。
2.根据权利要求1所述的交流电力系统电源发电功率调度的节点级分散方法,其特征在于,所述根据已知的交流电力系统的结构和参数,建立节点功率平衡的线性渐近方程包括:
根据所述交流电力系统的支路导纳参数、支路两端的电压幅值以及所述支路两端的电压相角,运用电功率定义式并撤分其中的耦合项,建立如下的支路传输功率的线性渐近表达式:
Pij=αijVi+βijVj+γijθi+δijθj
其中,Pij为支路ij传输的有功功率;Qij为所述支路ij传输的无功功率;i和j均为所述交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力系统中所述节点的总个数;gij为所述支路ij的电导参数;bij为所述支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为所述节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为所述节点j的电压相角变量;ξij是按照ξij=sinθij-θi+θj确定的所述支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的所述支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的所述支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的所述支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的所述支路ij的第4修正导纳;是按照确定的所述支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的所述支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的所述支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的所述支路ij的第8修正导纳;所述Vi和所述Vj均为标幺值电压;
根据所述支路传输功率的线性渐近表达式和所述交流电力系统的支路连接结构,按Kirchhoff电流定律建立如下的所述节点i的功率平衡的线性渐近方程:
其中,PGi为接于所述节点i的电源的有功功率变量;QGi为接于所述节点i的电源的无功功率变量;PDi为接于所述节点i的负荷的有功功率参数;QDi为接于所述节点i的负荷的无功功率参数。
3.根据权利要求1所述的交流电力系统电源发电功率调度的节点级分散方法,其特征在于,所述根据所述线性渐近方程、电源的有功供电成本以及无功消耗量,建立交流电力系统电源发电功率调度的二次规划模型包括:
以所述线性渐近方程为约束、以所述电源的有功供电成本与所述电源的无功消耗量的平方和之和最小为目标函数,建立如下的所述交流电力系统电源发电功率调度的二次规划模型:
其中,aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点;i和j均为所述交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力系统中所述节点的总个数;gij为支路ij的电导参数;bij为所述支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为所述节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为所述节点j的电压相角变量;Vn是所述编号为n的节点的电压幅值,且Vn=1是已知量;θn是所述编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的所述支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的所述支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的所述支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的所述支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的所述支路ij的第4修正导纳;是按照 确定的所述支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的所述支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的所述支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的所述支路ij的第8修正导纳;所述Vi和所述Vj均为标幺值电压;PGi为接于所述节点i的电源的有功功率变量;QGi为接于所述节点i的电源的无功功率变量;PDi为接于所述节点i的负荷的有功功率参数;QDi为接于所述节点i的负荷的无功功率参数。
4.根据权利要求1所述的交流电力系统电源发电功率调度的节点级分散方法,其特征在于,所述根据所述二次规划模型建立拉格朗日函数包括:
根据所述二次规划模型,按拉格朗日函数的定义建立如下的拉格朗日函数;
其中,是拉格朗日函数;λi为对应节点i的有功功率平衡方程的拉格朗日乘子;ξi为对应所述节点i的无功功率平衡方程的拉格朗日乘子;aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点;i和j均为所述交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力系统中所述节点的总个数;gij为支路ij的电导参数;bij为所述支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为所述节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为所述节点j的电压相角变量;Vn是所述编号为n的节点的电压幅值,且Vn=1是已知量;θn是所述编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的所述支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的所述支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的所述支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的所述支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的所述支路ij的第4修正导纳;是按照 确定的所述支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的所述支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的所述支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的所述支路ij的第8修正导纳;所述Vi和所述Vj均为标幺值电压;PGi为接于所述节点i的电源的有功功率变量;QGi为接于所述节点i的电源的无功功率变量;PDi为接于所述节点i的负荷的有功功率参数;QDi为接于所述节点i的负荷的无功功率参数。
5.根据权利要求1所述的交流电力系统电源发电功率调度的节点级分散方法,其特征在于,所述根据所述拉格朗日函数建立节点级分散迭代公式,继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值包括:
根据所述拉格朗日函数,按驻点的定义建立如下的驻点方程组:
其中,是拉格朗日函数;λi为对应节点i的有功功率平衡方程的拉格朗日乘子;ξi为对应所述节点i的无功功率平衡方程的拉格朗日乘子;aGi和cGi及dGi为节点i的电源的有功供电成本函数中二次项、一次项的系数及常数项;hGq为电源的无功消耗量的权重参数;编号为n的节点是交流电力系统电源发电功率调度的电压参考节点;i和j均为所述交流电力系统中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力系统中所述节点的总个数;gij为支路ij的电导参数;bij为所述支路ij的电纳参数;Vi为节点i的电压幅值变量;θi为所述节点i的电压相角变量;Vj为节点j的电压幅值变量;θj为所述节点j的电压相角变量;Vn是所述编号为n的节点的电压幅值,且Vn=1是已知量;θn是所述编号为n的节点的电压相角,且θn=0是已知量;ξij是按照ξij=sinθij-θi+θj确定的所述支路ij的导纳修正系数;αij是按照αij=gijVi-gijcosθijVj/2-bijVjθi/3+bijVjθj/3-bijVjξij/2确定的所述支路ij的第1修正导纳;βij是按照βij=-bijViθi/3-Vigijcosθij/2-bijViθj/3+bijViξij/2确定的所述支路ij的第2修正导纳;γij是按照γij=-bijViVj/3确定的所述支路ij的第3修正导纳;δij是按照δij=bijViVj/3确定的所述支路ij的第4修正导纳;是按照 确定的所述支路ij的第5修正导纳;φij是按照φij=bijcosθijVi/2-gijViθi/3+gijViθj/3-ξijgijVi/3确定的所述支路ij的第6修正导纳;χij是按照χij=-gijViVj/3确定的所述支路ij的第7修正导纳;ψij是按照ψij=gijViVj/3确定的所述支路ij的第8修正导纳;所述Vi和所述Vj均为标幺值电压;PGi为接于所述节点i的电源的有功功率变量;QGi为接于所述节点i的电源的无功功率变量;PDi为接于所述节点i的负荷的有功功率参数;QDi为接于所述节点i的负荷的无功功率参数;
基于所述驻点方程组,建立如下的节点级分散迭代公式,继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值:
其中,(t+1)表示第t+1步的迭代结果;(t)表示第t步的迭代结果;σ为大于0小于1的惯性参数;Ωi是编号为i的节点的所有邻居节点之编号集合。
6.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至5任一项所述交流电力系统电源发电功率调度的节点级分散方法的步骤。
7.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述交流电力系统电源发电功率调度的节点级分散方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110599963.8A CN113381453B (zh) | 2021-05-31 | 2021-05-31 | 交流电力系统电源发电功率调度的节点级分散方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110599963.8A CN113381453B (zh) | 2021-05-31 | 2021-05-31 | 交流电力系统电源发电功率调度的节点级分散方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113381453A true CN113381453A (zh) | 2021-09-10 |
CN113381453B CN113381453B (zh) | 2022-05-06 |
Family
ID=77575046
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110599963.8A Expired - Fee Related CN113381453B (zh) | 2021-05-31 | 2021-05-31 | 交流电力系统电源发电功率调度的节点级分散方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113381453B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104158191A (zh) * | 2014-09-09 | 2014-11-19 | 武汉大学 | 一种多机电力系统稳定运行的分散协调控制方法 |
US20150051744A1 (en) * | 2013-08-19 | 2015-02-19 | Board Of Trustees Of Michigan State University | Linear Optimal Power Flow System and Method |
CN104659812A (zh) * | 2015-01-31 | 2015-05-27 | 华南理工大学 | 一种基于预测控制的多微网协调控制方法 |
US20170264094A1 (en) * | 2016-03-09 | 2017-09-14 | Mitsubishi Electric Research Laboratories, Inc. | Decentralized Control of Electricity Passing through Electrical Grid |
CN107529678A (zh) * | 2017-10-11 | 2018-01-02 | 中国南方电网有限责任公司 | 可扩展安全稳定约束的交直流电网多目标优化控制方法 |
CN109256810A (zh) * | 2018-11-14 | 2019-01-22 | 南京邮电大学 | 考虑风机出力不确定成本的多目标优化方法 |
CN109314390A (zh) * | 2017-05-15 | 2019-02-05 | 深圳大学 | 获取直流电力网潮流的等量电导补偿型全局线性对称方法 |
CN111130118A (zh) * | 2020-01-09 | 2020-05-08 | 清华大学 | 一种基于分段线性化的电力系统最优潮流计算方法 |
CN112200498A (zh) * | 2020-09-02 | 2021-01-08 | 国网冀北电力有限公司承德供电公司 | 一种基于动态调度决策分析的电网经济运行方法 |
-
2021
- 2021-05-31 CN CN202110599963.8A patent/CN113381453B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150051744A1 (en) * | 2013-08-19 | 2015-02-19 | Board Of Trustees Of Michigan State University | Linear Optimal Power Flow System and Method |
CN104158191A (zh) * | 2014-09-09 | 2014-11-19 | 武汉大学 | 一种多机电力系统稳定运行的分散协调控制方法 |
CN104659812A (zh) * | 2015-01-31 | 2015-05-27 | 华南理工大学 | 一种基于预测控制的多微网协调控制方法 |
US20170264094A1 (en) * | 2016-03-09 | 2017-09-14 | Mitsubishi Electric Research Laboratories, Inc. | Decentralized Control of Electricity Passing through Electrical Grid |
CN109314390A (zh) * | 2017-05-15 | 2019-02-05 | 深圳大学 | 获取直流电力网潮流的等量电导补偿型全局线性对称方法 |
CN107529678A (zh) * | 2017-10-11 | 2018-01-02 | 中国南方电网有限责任公司 | 可扩展安全稳定约束的交直流电网多目标优化控制方法 |
CN109256810A (zh) * | 2018-11-14 | 2019-01-22 | 南京邮电大学 | 考虑风机出力不确定成本的多目标优化方法 |
CN111130118A (zh) * | 2020-01-09 | 2020-05-08 | 清华大学 | 一种基于分段线性化的电力系统最优潮流计算方法 |
CN112200498A (zh) * | 2020-09-02 | 2021-01-08 | 国网冀北电力有限公司承德供电公司 | 一种基于动态调度决策分析的电网经济运行方法 |
Non-Patent Citations (2)
Title |
---|
张荣权等: "基于改进萤火虫算法的冷热电联供系统多目标优化调度", 《华北电力大学学报》 * |
查永星等: "基于多目标多任务进化算法的含可再生能源混合发电系统优化调度", 《华北电力大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113381453B (zh) | 2022-05-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Moradi et al. | Power flow analysis in islanded Micro-Grids via modeling different operational modes of DGs: A review and a new approach | |
WO2019134254A1 (zh) | 一种运用分布式神经网络的实时经济调度计算方法 | |
Watson et al. | Optimized dispatch of energy storage systems in unbalanced distribution networks | |
Zamzam et al. | Beyond relaxation and Newton–Raphson: Solving AC OPF for multi-phase systems with renewables | |
Nazir et al. | Grid-aware aggregation and realtime disaggregation of distributed energy resources in radial networks | |
Zheng et al. | Fully distributed multi‐area economic dispatch method for active distribution networks | |
US20180158152A1 (en) | Methods and Software for Calculating Optimal Power Flow in an Electrical Power Grid and Utilizations of Same | |
Lai et al. | Distributed power control for DERs based on networked multiagent systems with communication delays | |
Spanias et al. | A system reference frame approach for stability analysis and control of power grids | |
CN105406459B (zh) | 利用基于一致性分布式处理估计电网中的最优潮流的方法 | |
Gavriluta et al. | A distributed approach for OPF-based secondary control of MTDC systems | |
Chen et al. | Distributed finite-step iterative algorithm for economic dispatch of generation | |
Yan et al. | Distributed continuous-time algorithm for economic dispatch problem over switching communication topology | |
Babayomi et al. | Distributed consensus-based reactive power sharing in microgrids: A predictive virtual capacitance control technique | |
Xia et al. | Optimal secondary control of islanded AC microgrids with communication time-delay based on multi-agent deep reinforcement learning | |
CN113381453B (zh) | 交流电力系统电源发电功率调度的节点级分散方法 | |
Han et al. | Leader‐distributed follower‐decentralized control strategy for economic dispatch in cascaded‐parallel microgrids | |
CN113381402B (zh) | 获取交流电力系统最小状态潮流的节点级分散方法 | |
Mallick et al. | Optimal decentralized microgrid coordination via the schur’s complement and s-procedure | |
Velasco et al. | Complex power sharing is not complex | |
CN113381397B (zh) | 获取直流电力系统最小状态潮流的节点级分散方法 | |
CN113224788B (zh) | 含有风电接入的交直流混联电网静态电压稳定评估方法 | |
Shen et al. | Coordinated decentralized reactive power optimization of multi-area power systems with discrete variables based on the alternating direction method of multipliers and Ward equivalence | |
CN108920879A (zh) | 移频建模仿真方法及装置 | |
Farhangi et al. | An analytical approach for the stability analysis of power networks through Kuramoto oscillators model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220506 |
|
CF01 | Termination of patent right due to non-payment of annual fee |