CN113381453A - 交流电力系统电源发电功率调度的节点级分散方法 - Google Patents

交流电力系统电源发电功率调度的节点级分散方法 Download PDF

Info

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
Application number
CN202110599963.8A
Other languages
English (en)
Other versions
CN113381453B (zh
Inventor
江辉
彭建春
胡先昌
张凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen University
Original Assignee
Shenzhen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen University filed Critical Shenzhen University
Priority to CN202110599963.8A priority Critical patent/CN113381453B/zh
Publication of CN113381453A publication Critical patent/CN113381453A/zh
Application granted granted Critical
Publication of CN113381453B publication Critical patent/CN113381453B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
Figure BDA0003092386030000031
其中，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修正导纳；
Figure BDA0003092386030000032
是按照
Figure BDA0003092386030000033
Figure BDA0003092386030000041
确定的支路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的功率平衡的线性渐近方程：
Figure BDA0003092386030000042
Figure BDA0003092386030000043
其中，PGi为接于节点i的电源的有功功率变量，且为未知量；QGi为接于节点i的电源的无功功率变量，且为未知量；PDi为接于节点i的负荷的有功功率参数，且为已知量；QDi为接于节点i的负荷的无功功率参数，且为已知量。
上述节点功率平衡的线性渐近方程是关于节点电压幅值和相角的线性方程，且节点功率平衡的线性渐近方程随节点电压幅值和相角逼近真值而逼近按照电功率定义和Kirchhoff电流定律得到的精确的节点功率平衡方程。这正是称上述线性渐近方程为节点功率平衡的线性渐近方程的缘故。
在步骤102中，根据线性渐近方程、电源的有功供电成本以及无功消耗量，建立交流电力系统电源发电功率调度的二次规划模型。
步骤102包括：以线性渐近方程为约束、以电源的有功供电成本与电源的无功消耗量的平方和之和最小为目标函数，建立如下的交流电力系统电源发电功率调度的二次规划模型：
Figure BDA0003092386030000044
Figure BDA0003092386030000051
Figure BDA0003092386030000052
其中，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修正导纳；
Figure BDA0003092386030000053
是按照
Figure BDA0003092386030000054
Figure BDA0003092386030000055
确定的支路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包括：根据二次规划模型，按拉格朗日函数的定义建立如下的拉格朗日函数。
Figure BDA0003092386030000061
其中，
Figure BDA0003092386030000062
是拉格朗日函数；λ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修正导纳；
Figure BDA0003092386030000071
是按照
Figure BDA0003092386030000072
确定的支路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中，根据拉格朗日函数，按驻点的定义建立如下的驻点方程组：
Figure BDA0003092386030000073
其中，
Figure BDA0003092386030000074
是拉格朗日函数；λ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修正导纳；
Figure BDA0003092386030000081
是按照
Figure BDA0003092386030000082
确定的支路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中，基于驻点方程组，建立如下的节点级分散迭代公式，继而根据节点级分散迭代公式获取交流电力系统的电源发电功率调度值：
Figure BDA0003092386030000091
其中，(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
Figure FDA0003092386020000011
其中，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修正导纳；
Figure FDA0003092386020000021
是按照
Figure FDA0003092386020000022
确定的所述支路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的功率平衡的线性渐近方程：
Figure FDA0003092386020000023
Figure FDA0003092386020000024
其中，PGi为接于所述节点i的电源的有功功率变量；QGi为接于所述节点i的电源的无功功率变量；PDi为接于所述节点i的负荷的有功功率参数；QDi为接于所述节点i的负荷的无功功率参数。
3.根据权利要求1所述的交流电力系统电源发电功率调度的节点级分散方法，其特征在于，所述根据所述线性渐近方程、电源的有功供电成本以及无功消耗量，建立交流电力系统电源发电功率调度的二次规划模型包括：
以所述线性渐近方程为约束、以所述电源的有功供电成本与所述电源的无功消耗量的平方和之和最小为目标函数，建立如下的所述交流电力系统电源发电功率调度的二次规划模型：
Figure FDA0003092386020000025
Figure FDA0003092386020000031
Figure FDA0003092386020000032
其中，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修正导纳；
Figure FDA0003092386020000033
是按照
Figure FDA0003092386020000034
Figure FDA0003092386020000035
确定的所述支路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所述的交流电力系统电源发电功率调度的节点级分散方法，其特征在于，所述根据所述二次规划模型建立拉格朗日函数包括：
根据所述二次规划模型，按拉格朗日函数的定义建立如下的拉格朗日函数；
Figure FDA0003092386020000041
其中，
Figure FDA0003092386020000042
是拉格朗日函数；λ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修正导纳；
Figure FDA0003092386020000043
是按照
Figure FDA0003092386020000044
Figure FDA0003092386020000045
确定的所述支路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所述的交流电力系统电源发电功率调度的节点级分散方法，其特征在于，所述根据所述拉格朗日函数建立节点级分散迭代公式，继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值包括：
根据所述拉格朗日函数，按驻点的定义建立如下的驻点方程组：
Figure FDA0003092386020000051
其中，
Figure FDA0003092386020000052
是拉格朗日函数；λ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修正导纳；
Figure FDA0003092386020000061
是按照
Figure FDA0003092386020000062
Figure FDA0003092386020000063
确定的所述支路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的负荷的无功功率参数；
基于所述驻点方程组，建立如下的节点级分散迭代公式，继而根据所述节点级分散迭代公式获取所述交流电力系统的电源发电功率调度值：
Figure FDA0003092386020000071
其中，(t+1)表示第t+1步的迭代结果；(t)表示第t步的迭代结果；σ为大于0小于1的惯性参数；Ωi是编号为i的节点的所有邻居节点之编号集合。
6.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至5任一项所述交流电力系统电源发电功率调度的节点级分散方法的步骤。
7.一种终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述交流电力系统电源发电功率调度的节点级分散方法的步骤。
CN202110599963.8A 2021-05-31 2021-05-31 交流电力系统电源发电功率调度的节点级分散方法 Active CN113381453B (zh)

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 Active CN113381453B (zh) 2021-05-31 2021-05-31 交流电力系统电源发电功率调度的节点级分散方法

Country Status (1)

Country Link
CN (1) CN113381453B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
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 国网冀北电力有限公司承德供电公司 一种基于动态调度决策分析的电网经济运行方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
张荣权等: "基于改进萤火虫算法的冷热电联供系统多目标优化调度", 《华北电力大学学报》 *
查永星等: "基于多目标多任务进化算法的含可再生能源混合发电系统优化调度", 《华北电力大学学报》 *

Also Published As

Publication number Publication date
CN113381453B (zh) 2022-05-06

Similar Documents

Publication Publication Date Title
Yang et al. A linearized OPF model with reactive power and voltage magnitude: A pathway to improve the MW-only DC OPF
Lu et al. Distributed secondary voltage and frequency control for islanded microgrids with uncertain communication links
Watson et al. Optimized dispatch of energy storage systems in unbalanced distribution networks
Robbins et al. Optimal reactive power dispatch for voltage regulation in unbalanced distribution systems
WO2019134254A1 (zh) 一种运用分布式神经网络的实时经济调度计算方法
Moradi et al. Power flow analysis in islanded Micro-Grids via modeling different operational modes of DGs: A review and a new approach
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
Nazir et al. Grid-aware aggregation and realtime disaggregation of distributed energy resources in radial networks
Gavriluta et al. A distributed approach for OPF-based secondary control of MTDC systems
CN105406459B (zh) 利用基于一致性分布式处理估计电网中的最优潮流的方法
CN109560547A (zh) 一种考虑输配协同的主动配电网n-1安全评估方法
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
CN113381453B (zh) 交流电力系统电源发电功率调度的节点级分散方法
Han et al. Leader‐distributed follower‐decentralized control strategy for economic dispatch in cascaded‐parallel microgrids
CN113381402B (zh) 获取交流电力系统最小状态潮流的节点级分散方法
CN113381397B (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
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
CN107947153B (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