CN115036995B - 一种多级时间复用的无功优化控制方法 - Google Patents
一种多级时间复用的无功优化控制方法 Download PDFInfo
- Publication number
- CN115036995B CN115036995B CN202210640964.7A CN202210640964A CN115036995B CN 115036995 B CN115036995 B CN 115036995B CN 202210640964 A CN202210640964 A CN 202210640964A CN 115036995 B CN115036995 B CN 115036995B
- Authority
- CN
- China
- Prior art keywords
- reactive power
- matrix
- power source
- time
- switching
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 126
- 238000004146 energy storage Methods 0.000 claims abstract description 23
- 230000035945 sensitivity Effects 0.000 claims abstract description 23
- 239000003990 capacitor Substances 0.000 claims abstract description 22
- 230000004044 response Effects 0.000 claims description 20
- 238000013507 mapping Methods 0.000 claims description 8
- 230000006798 recombination Effects 0.000 claims description 5
- 238000005215 recombination Methods 0.000 claims description 5
- 238000012163 sequencing technique Methods 0.000 claims description 5
- 238000004891 communication Methods 0.000 claims description 4
- 238000011161 development Methods 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000008093 supporting effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
Abstract
本发明公开了一种多级时间复用的无功优化控制方法,建立包括地调及新能源场站各无功源的投切时序的矩阵模型;建立包括地调及新能源场站各无功源的投切增量的矩阵模型;对时序矩阵和增量矩阵进行求解;对时序矩阵和增量矩阵进行灵敏度优化。本发明综合考虑电容器、电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的调压能力和特性,在各无功源的投切时序和投切增量上进行快速建模,形成多时间复用的无功控制方法,并通过不断的更新模型矩阵达到适应快速的网源变化的能力。
Description
技术领域
本发明涉及无功控制技术,具体地,涉及一种多级时间复用的无功优化控制方法。
背景技术
在智能电网背景下,在“碳达峰”、“碳中和”的目标驱动下,新能源的发展引起了广泛的关注、研究和应用,尤其风电、光伏的普及以及储能的应用多样化将占据几十年的发展势头已成共识。传统以火电、水电主导的电网的无功控制方案在大比例新能源接入后还能起到多大的支撑作用受到严峻的挑战,如何充分挖掘在大量新能源接入后电网的无功控制能力已经引起了广泛的关注。
早期无功电压优化控制的理论研究大都基于新能源具有统一的标准形式,而现状是新能源厂家的控制方案各具特色、形式各异,在电网侧、场站侧、厂家侧难以达成一个让各方满意的技术方案。在大量新能源接入传统能源电网后,为解决各无功源之间的协调配合问题,专家学者在优化算法方面做了大量研究,虽取得了一定成效,但也越来越不适应网源的快速发展,电压的日波幅和越限问题愈加频发。
因此本发明提出一种多级时间复用的无功优化控制方法。
发明内容
为解决现有技术中存在的不足,本发明的目的在于,提供一种多级时间复用的无功优化控制方法,综合考虑电容器、电抗器、FC、变压器分接头、风机、光伏逆变器、储能装置的调压能力和特性,在各无功源的投切时序和投切增量上进行快速建模,形成多时间复用的无功控制方法,并通过不断的更新模型矩阵达到适应快速的网源变化的能力。
本发明采用如下的技术方案。
一种多级时间复用的无功优化控制方法,所述无功优化控制方法包括以下步骤:
(1)建立包括地调及新能源场站各无功源的投切时序的矩阵模型;
(2)建立包括地调及新能源场站各无功源的投切增量的矩阵模型;
(3)对时序矩阵和增量矩阵进行求解;
(4)对时序矩阵和增量矩阵进行灵敏度优化。
进一步地,步骤(1)中,根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的调压特性,建立各无功源的投切时序矩阵,包括状态矩阵和时间响应矩阵。
进一步地,对于各无功源的状态矩阵,考虑因素包括是否闭锁、是否故障、是否有剩余可调容量;状态矩阵A:
式中,元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如aw1即代表风电场1的可调状态,awn即代表风电场n的可调状态,其他元素代表涵义以此类推;其状态值通过设为0或1来描述,若为0则表示当前状态不可参与无功调节,若为1则表示当前状态可参与无功调节。
进一步地,对于各无功源的时间响应矩阵,考虑因素包括通讯时长和执行时长;时间响应矩阵B:
式中,各元素代表各无功源的响应速度,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头;如tw1即代表风电场1的响应时间,twn即代表风电场n的响应时间,其他元素以此类推。
进一步地,步骤(2)中,根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的容量,建立各无功源的投切增量矩阵,包括状态矩阵和容量矩阵。
进一步地,对于容量矩阵,考虑因素包括各无功源的无功总容量、需对外提供的无功总容量、所需保留的无功裕度容量;
容量矩阵D:
式中,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如dw1即代表风电场1的可调容量,dwn即代表风电场n的可调状态,其他元素代表涵义以此类推。
进一步地,步骤(3)中,对时序矩阵和增量矩阵求解具体步骤为,
对各无功源投切时序矩阵中的状态矩阵和时间响应矩阵进行一一映射并进行元素排序重组,形成地调和新能源场站各自无功源的投切时序;对各无功源的投切增量矩阵中的状态矩阵和容量矩阵进行一一映射并求取映射后矩阵的特征值,以特征值大小作为对应无功源无功出力大小的权重,得到地调和新能源场站各自无功源的投切增量。
进一步地,步骤(4)中,以各无功源对某节点的电压支撑程度建立各无功源对某节点电压的调节灵敏度,以此灵敏度优化已形成的各无功源的投切时序与增量;并对该灵敏度持续在线更新;对于优化结果中存在同等无功源优先级的用分时加权矩阵进一步优化。
进一步地,灵敏度描述如下:
式中,下标j表示不同的无功源,下标i表示电网中的不同节点,Sj→i为某无功源j对某节点i的无功调节灵敏度;ΔQj为某无功源j参与的无功调节量,ΔVj→i为对应单一无功源j参与调节后对某节点i电压影响的幅值大小;
分时加权矩阵[Wji]j×n描述如下:
式中,j为不同种类的无功源,n为不同种类无功源中最大的列数,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如ww1即代表风电场1的加权值,wwn即代表风电场n的加权值,其他元素代表涵义以此类推。
一种多级时间复用的无功优化控制系统,所述系统用于实现所述的多级时间复用的无功优化控制方法。
本发明的有益效果在于,与现有技术相比,本发明综合考虑电容器、电抗器、FC、变压器分接头、风机、光伏逆变器、储能装置的调压能力和特性,在各无功源的投切时序和投切增量上进行快速建模,形成多时间复用的无功控制方法,并通过不断的更新模型矩阵达到适应快速的网源变化的能力。
附图说明
图1为本发明所述的多级时间复用的无功优化控制方法流程图。
具体实施方式
下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本申请的保护范围。
如图1所示,本发明所述的一种多级时间复用的无功优化控制方法,包括如下步骤:
(1)建立包括地调及新能源场站各无功源的投切时序的矩阵模型;
根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的调压特性,建立各无功源的投切时序矩阵:状态矩阵和时间响应矩阵。
对于各无功源的状态矩阵,考虑的因素包括:是否闭锁、是否故障、是否有剩余可调容量;其状态值通过设为0或1来描述,若为0则表示当前状态不可参与无功调节,若为1则表示当前状态可参与无功调节。
状态矩阵描述如下:
式中:矩阵A即为各无功源的状态矩阵,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如aw1即代表风电场1的可调状态,awn即代表风电场n的可调状态,其他元素代表涵义以此类推;对于长时间参与无功电压调节的情况,优先考虑矩阵中aci/afci/ali/ati;对于短时间参与无功电压调节的情况,优先考虑矩阵中的awi/apvi/asi/asvgi,此处不同无功源的排列无优先级区别。
对于各无功源的时间响应矩阵,考虑的因素包括:通讯时长和执行时长。通讯时长包括从命令下达到场站接受到命令这段时间,设为t1;各无功源执行完毕命令的时间设为t2;t1代表命令的传输速度,与网络和通道方式相关,可用两地距离描述t1=ks,s即为主站至某场站的距离,k为系数(可由通道方式和装置时延得出),也可由现场调试得出;t2为各无功源的命令执行时间,不同类型的无功源有所差异,可由实际调试得出。
时间响应矩阵描述如下:
式中:矩阵B为各无功源的响应时间矩阵,其中各元素代表各无功源的响应速度;其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头;如tw1即代表风电场1的响应时间,twn即代表风电场n的响应时间,其他元素以此类推。
(2)建立包括地调及新能源场站各无功源的投切增量的矩阵模型;
根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的容量,建立各无功源的投切增量矩阵:状态矩阵和容量矩阵。
对于状态矩阵同步骤(1)中所述。
对于容量矩阵,考虑因素包括:各无功源的无功总容量、需对外提供的无功总容量、所需保留的无功裕度容量;容量矩阵中的元素以最终能够向外提供的无功容量为基值。
容量矩阵描述如下:
式中:矩阵D即为各无功源的容量矩阵,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如dw1即代表风电场1的可调容量,dwn即代表风电场n的可调状态,其他元素代表涵义以此类推。
(3)对时序矩阵和增量矩阵进行求解;
对各无功源投切时序矩阵中的状态矩阵和时间响应矩阵进行一一映射并进行元素排序重组,形成地调和新能源场站各自无功源的投切时序;对各无功源的投切增量矩阵中的状态矩阵和容量矩阵进行一一映射并求取映射后矩阵的特征值,以特征值大小作为对应无功源无功出力大小的权重,得到地调和新能源场站各自无功源的投切增量。
先将矩阵A的元素一一映射到B矩阵中得到B'矩阵,将不参与无功调节的元素过滤,即B'=M(A→B);然后将B'矩阵进行排序重组得到矩阵S,也即S=sort(B');则矩阵S即为各无功源的执行时序矩阵。
同样的,将A的元素一一映射到D矩阵中得到D'矩阵,即D'=M(A→D);然后将D'矩阵进行排序重组得到矩阵S',也即S'=sort(D');则矩阵S'即为各无功源的执行增量矩阵。
(4)对时序矩阵和增量矩阵进行灵敏度优化;
以应对实际复杂的轻重缓急的无功需求为目的,直接以各无功源对某节点的电压支撑程度建立各无功源对某节点电压的调节灵敏度,以此灵敏度优化已形成的各无功源的投切时序与增量;并对该灵敏度持续在线更新;对于优化结果中存在同等无功源优先级的用分时加权矩阵进一步优化。
该步骤在调试、试运行阶段形成初级优化样本,并在实际运行阶段持续更新样本以适应动态的网络变化和无功需求的变化。其中灵敏度描述如下:
式中:下标j表示不同的无功源,下标i表示电网中的不同节点,Sj→i为某无功源j对某节点i的无功调节灵敏度;ΔQj为某无功源j参与的无功调节量,ΔVj→i为对应单一无功源j参与调节后对某节点i电压影响的幅值大小。
在调试期间,形成原始灵敏度矩阵Sr(j→i),试运行期间和正式运行期间持续更新矩阵Sr(j→i),为下一次无功源的投切时序和增量做优化参考,以适应实际复杂的实际工况。
对形成的灵敏度矩阵与无功源状态矩阵一一映射,即Sr'=M(A→Sr),过滤掉不参与本次无功调压的无功源,并对矩阵元素排序,得到地调和新能源场站各自可参与到本次无功调压的灵敏度矩阵。
最后对时序矩阵S、增量矩阵S'、灵敏度矩阵Sr'中的各类无功源求取交集,得到最终的各无功源的投切时序和增量。若均无交集或两两无交集则以满足无功增量优先;对其中存在同等优先级的无功源用分时加权矩阵进一步优化评估,即在设定时间范围内,每调节成功一次即加一分,形成分时加权矩阵[Wji]j×n,超过设定时间范围,加权值清零,重新计数。
分时加权矩阵描述如下:
式中:j为不同种类的无功源,n为不同种类无功源中最大的列数,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头,如ww1即代表风电场1的加权值,wwn即代表风电场n的加权值,其他元素代表涵义以此类推。
本发明的有益效果在于,与现有技术相比,本发明综合考虑电容器、电抗器、FC、变压器分接头、风机、光伏逆变器、储能装置的调压能力和特性,在各无功源的投切时序和投切增量上进行快速建模,形成多时间复用的无功控制方法,并通过不断的更新模型矩阵达到适应快速的网源变化的能力。
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施示例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。
Claims (5)
1.一种多级时间复用的无功优化控制方法,其特征在于,所述无功优化控制方法包括以下步骤:
(1)建立包括地调及新能源场站各无功源的投切时序的矩阵模型;
根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的调压特性,建立各无功源的投切时序矩阵,包括状态矩阵和时间响应矩阵;
对于各无功源的状态矩阵,考虑因素包括是否闭锁、是否故障、是否有剩余可调容量;状态矩阵A:
式中,元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头;其状态值通过设为0或1来描述,若为0则表示当前状态不可参与无功调节,若为1则表示当前状态可参与无功调节;
对于各无功源的时间响应矩阵,考虑因素包括通讯时长和执行时长;时间响应矩阵B:
式中,各元素代表各无功源的响应速度,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头;
(2)建立包括地调及新能源场站各无功源的投切增量的矩阵模型;
根据电容器电抗器、FC、变压器分接头、风机、光伏逆变器、SVG及储能装置的容量,建立各无功源的投切增量矩阵,包括状态矩阵和容量矩阵;
对于容量矩阵,考虑因素包括各无功源的无功总容量、需对外提供的无功总容量、所需保留的无功裕度容量;
容量矩阵D:
式中,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头;
(3)对时序矩阵和增量矩阵进行求解;
(4)对时序矩阵和增量矩阵进行灵敏度优化。
2.根据权利要求1所述的多级时间复用的无功优化控制方法,其特征在于,
步骤(3)中,对时序矩阵和增量矩阵求解具体步骤为,
对各无功源投切时序矩阵中的状态矩阵和时间响应矩阵进行一一映射并进行元素排序重组,形成地调和新能源场站各自无功源的投切时序;对各无功源的投切增量矩阵中的状态矩阵和容量矩阵进行一一映射并求取映射后矩阵的特征值,以特征值大小作为对应无功源无功出力大小的权重,得到地调和新能源场站各自无功源的投切增量。
3.根据权利要求1所述的多级时间复用的无功优化控制方法,其特征在于,
步骤(4)中,以各无功源对某节点的电压支撑程度建立各无功源对某节点电压的调节灵敏度,以此灵敏度优化已形成的各无功源的投切时序与增量;并对该灵敏度持续在线更新;对于优化结果中存在同等无功源优先级的用分时加权矩阵进一步优化。
4.根据权利要求3所述的多级时间复用的无功优化控制方法,其特征在于,
灵敏度描述如下:
式中,下标j表示不同的无功源,下标i表示电网中的不同节点,Sj→i为某无功源j对某节点i的无功调节灵敏度;ΔQj为某无功源j参与的无功调节量,ΔVj→i为对应单一无功源j参与调节后对某节点i电压影响的幅值大小;
分时加权矩阵[Wji]j×n描述如下:
式中,j为不同种类的无功源,n为不同种类无功源中最大的列数,其中元素下标w/pv/s/svg/c/fc/l/t分别代表风机、光伏逆变器、储能、SVG、电容器、FC、电抗器、变压器分接头。
5.一种多级时间复用的无功优化控制系统,其特征在于,所述系统用于实现权利要求1-4任意一项所述的多级时间复用的无功优化控制方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210640964.7A CN115036995B (zh) | 2022-06-08 | 一种多级时间复用的无功优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210640964.7A CN115036995B (zh) | 2022-06-08 | 一种多级时间复用的无功优化控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115036995A CN115036995A (zh) | 2022-09-09 |
CN115036995B true CN115036995B (zh) | 2024-07-02 |
Family
ID=
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106655204A (zh) * | 2016-11-10 | 2017-05-10 | 华北电力大学 | 基于多无功源互动的风电场/群无功电压实时控制方法 |
CN108964023A (zh) * | 2018-06-29 | 2018-12-07 | 国网上海市电力公司 | 一种用于电网的母线电压态势短期预测方法及系统 |
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106655204A (zh) * | 2016-11-10 | 2017-05-10 | 华北电力大学 | 基于多无功源互动的风电场/群无功电压实时控制方法 |
CN108964023A (zh) * | 2018-06-29 | 2018-12-07 | 国网上海市电力公司 | 一种用于电网的母线电压态势短期预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110247438B (zh) | 基于天牛须算法的主动配电网资源优化配置 | |
CN107622324A (zh) | 一种考虑多微网能量交互的鲁棒环境经济调度方法 | |
CN102684201B (zh) | 一种基于电压越限概率的含风电场电网无功优化方法 | |
CN103219748A (zh) | 一种电网风电消纳能力的确定方法及系统 | |
CN103944175A (zh) | 风光储联合发电系统出力特性优化方法 | |
CN111614110B (zh) | 一种基于改进多目标粒子群算法的受端电网储能优化配置方法 | |
CN113555908B (zh) | 一种智能配电网储能优化配置方法 | |
CN107611970B (zh) | 分布式光伏和电动汽车的不确定性配网的优化方法 | |
CN115036995B (zh) | 一种多级时间复用的无功优化控制方法 | |
CN107294105B (zh) | 分布式光伏集群无通信条件下的动态调压控制方法 | |
CN112332421A (zh) | 基于自适应下垂控制的光伏电站参与电网电压调节方法 | |
CN108416132B (zh) | 一种分布式光伏集群的自动变步长仿真加速方法 | |
CN111478335A (zh) | 一种计及分布式光伏的配电网潮流计算方法及系统 | |
CN109508062B (zh) | 一种基于模糊电导的光伏发电控制方法及系统 | |
CN114530848B (zh) | 一种光储虚拟电厂多时间尺度动态划分方法 | |
CN115036995A (zh) | 一种多级时间复用的无功优化控制方法 | |
CN106410867B (zh) | 一种电能路由器功率流优化分析方法 | |
CN110829440A (zh) | 一种三相中低压一体化配电网电压控制方法及系统 | |
CN117526299B (zh) | 一种微电网有功无功功率协调控制系统及方法 | |
CN116470564A (zh) | 多无功源协同的海上风电场无功电压自动控制系统及方法 | |
CN117910836B (zh) | 一种提升大电网灵活性的储能电站规划方法 | |
CN114221385B (zh) | 一种考虑电动汽车用户出行需求的柔性充电方法及系统 | |
CN108054758A (zh) | 新能源电站电压平衡优化方法及存储介质 | |
CN115659779B (zh) | 一种多直流馈入受端电网的新能源接入优化策略 | |
CN203243117U (zh) | 风光储综合控制系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |