CN107171354A - 柔性直流电网换流站下垂系数计算方法 - Google Patents
柔性直流电网换流站下垂系数计算方法 Download PDFInfo
- Publication number
- CN107171354A CN107171354A CN201710466233.4A CN201710466233A CN107171354A CN 107171354 A CN107171354 A CN 107171354A CN 201710466233 A CN201710466233 A CN 201710466233A CN 107171354 A CN107171354 A CN 107171354A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- power
- vsc
- station
- 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
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
-
- 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]
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Direct Current Feeding And Distribution (AREA)
Abstract
本发明涉及一种柔性直流电网换流站下垂系数计算方法:定义广义功率共享系数,选取设定公共参考端,直流电压作为下垂控制站的公共电压参考信号;将换流站VSCi端有功功率平衡关系在直流网络扰动前的稳态工作点按一阶泰勒级数展开,并写出矩阵形式;根据克莱姆法则求得系统中任一换流站端直流电压的变化量与功率阶跃量之间的线性解析关系;求解广义功率共享系数向量T的最小二乘问题得到广义功率共享系数向量T的最优解,再计算各下垂控制站的下垂系数。
Description
技术领域
本发明涉及高压柔性直流输电技术领域,特别涉及针对换流站N-1故障或发生大的功率阶跃扰动,柔性直流电网直流电压的控制策略。
背景技术
柔性直流输电网络内部某个换流站发生大的功率阶跃扰动或因故障退出运行,导致整个直流系统内部功率不平衡,引起各换流站直流电压升高或降低发生较大的变化,因而控制系统内部直流电压的稳定有极其重要的意义[1]。下垂控制因具备直流电压波动量小、动态响应速度快、稳定度高等技术优势[2,3],在现有直流网络电压控制方法中受到更多关注,其控制思想为:系统内部某个换流站因故障退出运行或发生功率阶跃扰动后,系统中具有功率调节能力的换流站均采用下垂控制,按照下垂特性共同分担系统内的不平衡功率。对于下垂控制而言,下垂系数是分析和研究的关键,其大小直接决定了下垂控制的控制性能。目前,柔性直流电网采用的下垂控制主要包括三类:一是固定下垂系数的下垂控制[3-5],下垂系数的取值由换流站额定容量来确定,由于未考虑各换流站实际可用传输容量的差异,按固定下垂系数进行功率调节,可能造成有的换流站满载运行,而有的换流站轻载运行。二是变下垂系数的下垂控制[6-8],下垂系数根据换流站功率裕度或直流电压的大小自适应调节,由于下垂系数的取值依赖于系统的当前运行状态,导致下垂系数的计算繁琐复杂,并且下垂系数的实时调整将频繁地改变直流系统功率的分配方式,引起直流电压的频繁波动。三是具有比例分担特性的下垂控制[9],选取某一端直流电压作为下垂控制站的公共电压参考信号,在忽略直流网络损耗变化的前提下,换流站下垂系数由负荷分担系数、系统功率阶跃量以及最大允许直流电压偏差来确定,然而下垂系数的确定仍面临需要人为给定最大允许直流电压偏差的问题。鉴于以上分析,下垂控制站下垂系数的快速、直接、通用的计算方法尚未形成,针对换流站N-1故障或发生大的功率阶跃扰动,有必要提出换流站下垂系数的快速直接计算新方法,在保证直流电网电压波动尽可能小的目标下实现柔性直流电网直流电压的快速稳定。
参考文献
[1]赵成勇.柔性直流输电建模和仿真技术[M].北京:中国电力出版社,2014.
[2]唐庚,徐政,刘昇,顾益磊,陆羿,裘鹏.适用于多端柔性直流输电系统的新型直流电压控制策略[J].电力系统自动化,2013,37(15):125-132.
[3]JefBeerten,Stijn Cole,Ronnie Belmans,Modeling of multi-terminalVSC HVDC systems with distributed DC voltage control[J],IEEE Trans.on PowerSystem,2014,29(1):34-42.
[4]Xu L,Williams B W,Yao L.Multi-terminal DC transmission systems forconnecting large offshore wind farms[C]//Power and Energy Society GeneralMeeting-Conversion and Delivery of Electrical Energy in the,Century.IEEE,2008:1-7.
[5]阎发友,汤广福,贺之渊,孔明.基于MMC的多端柔性直流输电系统改进下垂控制策略[J].中国电机工程学报,2014,34(03):397-404.
[6]Nilanjan Ray Chaudhuri,Balarko Chaudhuri,Adaptive Droop Controlfor Effective Power Sharing in Multi-Terminal DC(MTDC)Grids[J],IEEE Trans.onPower Systems,2013,28(1):21-29.
[7]刘瑜超,武健,刘怀远,等.基于自适应下垂调节的VSC-MTDC功率协调控制[J].中国电机工程学报,2016,36(1):40-48.
[8]陈朋,李梅航,严兵,等.适用于多端柔性直流输电系统的灵活下垂控制策略[J].电网技术,2016,40(11):3433-3440.
[9]Berggren B,Linden K,Majumder R.DC grid control through the pilotvoltage droop concept:methodology for establishing droop constants[J].IEEETransactions on Power Systems,2015,30(5):2312-2320.
发明内容
本发明的目的是克服现有技术的上述不足,针对换流站N-1故障或发生大的功率阶跃扰动,以整个直流电网电压变化尽可能小为目标,提出一种下垂控制站下垂系数的快速、直接计算新方法,将所计算的下垂系数引入到换流站有功外环控制器,通过执行换流站间有功功率的快速协调控制,可以实现直流电网电压的快速稳定。技术方案如下:
一种柔性直流电网换流站下垂系数计算方法,设定多端柔性直流系统有n个换流站,其中换流站VSC1~VSCm为下垂控制站,VSCm+1~VSCn为定有功功率站,VSCn端发生功率阶跃,功率阶跃量为ΔPsn,计算步骤如下:
(1)定义广义功率共享系数Ti,i=1,2,…,m,在换流站N-1故障或发生大的功率阶跃扰动后,系统中各下垂控制站的传输功率变化量ΔPsi=-TiΔPsn,
(2)选取VSCg端为公共参考端,设定此端直流电压作为下垂控制站的公共电压参考信号,对于下垂控制站VSCi,i=1,2,…,m,由其下垂特性可得
其中,Kdroop,i为VSCi端的下垂系数,ΔUdcg为VSCg端的直流电压变化量;
(3)将换流站VSCi端有功功率平衡关系在直流网络扰动前的稳态工作点处按一阶泰勒级数展开:
其中,Psi为PCC点传输的有功功率,Udci为换流站直流侧极对极电压,ΔPsi、ΔUdci分别表示Psi、Udci的变化量,Idci为换流站直流侧输出的电流,Yii、Yij分别表示直流网络节点导纳矩阵中的对角和非对角元素,上标带*的量表示直流系统扰动前变量的稳态值。
将式(2)写成矩阵形式:
ΔPs=CΔUdc (3)
其中,ΔPs、ΔUdc分别为ΔPsi、ΔUdcii=1,2,…,n,构成的向量,C为n阶系数方阵;
(4)对于具有n个换流站的多端柔性直流系统,换流站功率变化量因换流站类型的不同分为以下三种情况:
1)对于下垂控制站VSCii=1,2,…,m,功率变化量ΔPsi=-TiΔPsn;
2)对于定有功功率站VSCi,i=m+1,m+2,…,n-1,功率变化量ΔPsi=0;
3)对于发生功率阶跃的换流站VSCn,功率变化量为ΔPsn;
将式(3)ΔPs中的元素按上述三种情况对应代入,根据克莱姆法则,由式(3)求得系统中任一换流站VSCl端直流电压的变化量ΔUdcl与功率阶跃量ΔPsn之间的线性解析关系如下:
其中,l=1,2,…,n,det(C)表示方阵C的行列式,Mnl、Mkl分别是矩阵C的元素Cnl、Ckl的代数余子式;
取式(4)中下标l=g,并代入式(1)可得
将式(4)写成紧凑形式:
其中,L为由Mnl(l=1,2,…,n)构成的列向量;T为由Ti,i=1,2,…,m,构成的列向量;A为n×m的系数矩阵,Mkl为其第l行k列的元素。
(5)求解广义功率共享系数向量T的最小二乘问题:
其中Kn0为正实数,E为n维列向量,且E=[1 1... 1]T,
得到广义功率共享系数向量T的最优解,再计算各下垂控制站的下垂系数Kdroop,i,i=1,2,...,m。
本发明的有益效果如下:
1.快速性。本发明建立基于广义功率共享系数的最小二乘求解模型,通过方便地求解该模型,实现下垂系数的快速、直接计算。
2.通用性。本发明所提出的下垂控制站下垂系数的计算仅依赖于系统扰动前的运行状态和拓扑结构,且无需人为给定最大允许直流电压偏差,下垂系数的计算方法具有通用性。
3.适用性。所提下垂控制站下垂系数的计算不受扰动换流站位置及功率阶跃量大小的限制,对于系统内任一端换流站发生功率阶跃扰动或退出运行的情况均适用。
附图说明
图1多端直流系统换流站VSCi端模型
图2基于下垂控制的直流电网快速功率协调控制
具体实施方式
假定多端柔性直流系统有n个换流站,其中换流站VSC1~VSCm为下垂控制站,VSCm+1~VSCn为定有功功率站,VSCn端发生功率阶跃,功率阶跃量为ΔPsn。
为快速计算下垂控制站的下垂系数,本发明提出“广义功率共享系数”的概念,广义功率共享系数用Ti(i=1,2,…,m)表示,换流站N-1故障或发生大的功率阶跃扰动后,系统中各下垂控制站的传输功率变化量ΔPsi=-TiΔPsn。选取某一端直流电压作为下垂控制站的公共电压参考信号(假定VSCg端为公共参考端),对于下垂控制站VSCi(i=1,2,…,m),由其下垂特性可得
其中,Kdroop,i为VSCi端的下垂系数,ΔUdcg为VSCg端的直流电压变化量。
换流站VSCi端模型如图1所示,忽略Ri上的有功损耗的变化,将换流站VSCi端有功功率平衡关系在直流网络扰动前的稳态工作点处按一阶泰勒级数展开:
其中,Psi为PCC点传输的有功功率,Udci为换流站直流侧极对极电压,ΔPsi、ΔUdci分别表示Psi、Udci的变化量,Idci为换流站直流侧输出的电流,Yii、Yij分别表示直流网络节点导纳矩阵中的对角和非对角元素。上标带*的量表示直流系统扰动前变量的稳态值。
将式(2)写成矩阵形式:
ΔPs=CΔUdc (3)
其中,ΔPs、ΔUdc分别为ΔPsi、ΔUdci(i=1,2,...,n)构成的向量,C为n阶系数方阵,由式(2)可见,C矩阵元素值取决于扰动前直流网络的拓扑及运行状态,从而C为常数方阵。
对于具有n个换流站的多端柔性直流系统,换流站功率变化量因换流站类型的不同分为以下三种情况:
1)对于下垂控制站VSCi(i=1,2,…,m),功率变化量ΔPsi=-TiΔPsn;
2)对于定有功功率站VSCi(i=m+1,m+2,…,n-1),功率变化量ΔPsi=0;
3)对于发生功率阶跃的换流站VSCn,功率变化量为ΔPsn。
将式(3)ΔPs中的元素按上述三种情况对应代入,根据克莱姆法则,可由式(3)求得系统中任一换流站VSCl端直流电压的变化量ΔUdcl与功率阶跃量ΔPsn之间的线性解析关系如下:
其中,l=1,2,…,n,det(C)表示方阵C的行列式,Mnl、Mkl分别是矩阵C的元素Cnl、Ckl的代数余子式。由于C为常数方阵,故上述代数余子式均为常数。
取式(4)中下标l=g,并代入式(1)可得
由式(5)可知,下垂控制站的下垂系数仅和广义功率共享系数、扰动前直流网络的拓扑和运行状态、扰动换流站的位置有关,而与功率阶跃量ΔPsn无关。式(5)是本发明建立的巧妙的转换关系,为求取各下垂控制站的下垂系数,只需确定广义功率共享系数Ti(i=1,2,...,m)的取值。
将式(4)写成紧凑形式:
其中,L为由Mnl(l=1,2,…,n)构成的列向量;T为由Ti(i=1,2,...,m)构成的列向量;A为n×m的系数矩阵,Mkl为其第l行k列的元素。
由式(6)可知,对于确定的向量L和矩阵A,要使直流电网电压变化量ΔUdc尽可能小,应满足:
min||L-AT|| (7)
可得无约束条件下广义共享功率系数向量T的最小二乘解如下:
在系统的实际运行中,还需要考虑系统的运行约束条件。对于式(4),令
其中,Knl为描述功率阶跃量ΔPsn与直流电压变化量ΔUdcl数值关系的系数。
由系统运行特性可知,当系统内发生正向功率阶跃(即ΔPsn>0)时,注入直流电网的有功功率增加,换流站直流电压升高(即ΔUdcl>0);反之,当ΔPsn<0时,ΔUdcl<0。
因此,为满足系统运行特性,Knl必为正数,取Knl0为Knl的下限取值。为简单起见,将n个下限值取为同一个值,即Kn10=...=Knn0=Kn0。
将式(9)写成紧凑形式,结合上述分析及式(6),得到如下不等式约束条件:
其中,E为n维列向量,且E=[1 1 ... 1]T。式(10)可写成形如βT≤b的标准形式,其中
考虑式(7)与式(10),则求解广义功率共享系数向量T的最小二乘问题可用下述模型表示:
为便于求解式(11),将该式等价转化为二次规划问题,其数学模型如式(12)所示。
其中,H=2ATA;
求解式(12),得到广义功率共享系数向量T的最优解
将式(13)代入式(5),即可完成下垂控制站下垂系数的计算。将所确定的下垂系数加入下垂控制站外环控制器,即可实现系统内直流功率的快速协调控制,控制结构如图2所示。
概括而言,本发明的技术方案如下:
1.依照式(2)、(3)关系,由直流网络导纳矩阵元素及直流系统扰动前各换流站直流电压、电流形成系数矩阵C。
2.根据扰动换流站及下垂控制站位置,计算矩阵C的元素Cnl、Ckl的代数余子式Mnl(l=1,2,…,n)、Mkl(k=1,2,…,m,l=1,2,…,n),依照式(6)形成系数矩阵A和向量L。
3.求解式(7),由式(8)确定无约束条件下广义共享功率系数向量T的最小二乘解
4.考虑线性不等式约束条件式(10),建立含线性不等式约束的最小二乘求解模型式(11),将其等价转化为二次规划模型式(12)。
5.求解二次规划模型式(12),得到广义功率共享系数的最优解代入式(5)计算各下垂控制站的下垂系数Kdroop,i(i=1,2,...,m)。
6.将下垂系数Kdroop,i(i=1,2,...,m)加入图2所示的换流站有功外环控制器,执行直流系统的快速功率平衡控制。
Claims (1)
1.一种柔性直流电网换流站下垂系数计算方法,设定多端柔性直流系统有n个换流站,其中换流站VSC1~VSCm为下垂控制站,VSCm+1~VSCn为定有功功率站,VSCn端发生功率阶跃,功率阶跃量为ΔPsn,计算步骤如下:
(1)定义广义功率共享系数Ti,i=1,2,…,m,在换流站N-1故障或发生大的功率阶跃扰动后,系统中各下垂控制站的传输功率变化量ΔPsi=-TiΔPsn,
(2)选取VSCg端为公共参考端,设定此端直流电压作为下垂控制站的公共电压参考信号,对于下垂控制站VSCi,i=1,2,…,m,由其下垂特性可得
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>o</mi>
<mi>o</mi>
<mi>p</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>g</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>g</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Kdroop,i为VSCi端的下垂系数,ΔUdcg为VSCg端的直流电压变化量;
(3)将换流站VSCi端有功功率平衡关系在直流网络扰动前的稳态工作点i=1,2,…,n,处按一阶泰勒级数展开:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>j</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>i</mi>
</mrow>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
</munder>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>i</mi>
</mrow>
<mi>n</mi>
</munderover>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>i</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Psi为PCC点传输的有功功率,Udci为换流站直流侧极对极电压,ΔPsi、ΔUdci分别表示Psi、Udci的变化量,Idci为换流站直流侧输出的电流,Yii、Yij分别表示直流网络节点导纳矩阵中的对角和非对角元素,上标带*的量表示直流系统扰动前变量的稳态值。
将式(2)写成矩阵形式:
ΔPs=CΔUdc (3)
其中,ΔPs、ΔUdc分别为ΔPsi、ΔUdcii=1,2,…,n,构成的向量,C为n阶系数方阵;
(4)对于具有n个换流站的多端柔性直流系统,换流站功率变化量因换流站类型的不同分为以下三种情况:
1)对于下垂控制站VSCii=1,2,…,m,功率变化量ΔPsi=-TiΔPsn;
2)对于定有功功率站VSCi,i=m+1,m+2,…,n-1,功率变化量ΔPsi=0;
3)对于发生功率阶跃的换流站VSCn,功率变化量为ΔPsn;
将式(3)ΔPs中的元素按上述三种情况对应代入,根据克莱姆法则,由式(3)求得系统中任一换流站VSCl端直流电压的变化量ΔUdcl与功率阶跃量ΔPsn之间的线性解析关系如下:
<mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>k</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>det</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,l=1,2,…,n,det(C)表示方阵C的行列式,Mnl、Mkl分别是矩阵C的元素Cnl、Ckl的代数余子式;
取式(4)中下标l=g,并代入式(1)可得
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
<mi>o</mi>
<mi>o</mi>
<mi>p</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>k</mi>
<mi>g</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mi>det</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
1
将式(4)写成紧凑形式:
<mrow>
<msub>
<mi>&Delta;U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>-</mo>
<mi>A</mi>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>s</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>det</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,L为由Mnl(l=1,2,…,n)构成的列向量;T为由Ti,i=1,2,…,m,构成的列向量;A为n×m的系数矩阵,Mkl为其第l行k列的元素。
(5)求解广义功率共享系数向量T的最小二乘问题:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>L</mi>
<mo>-</mo>
<mi>A</mi>
<mi>T</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&beta;</mi>
<mi>T</mi>
<mo>&le;</mo>
<mi>b</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Kn0为正实数,E为n维列向量,且E=[1 1 ...1]T,得到广义功率共享系数向量T的最优解,再计算各下垂控制站的下垂系数Kdroop,i,i=1,2,...,m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710466233.4A CN107171354B (zh) | 2017-06-19 | 2017-06-19 | 柔性直流电网换流站下垂系数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710466233.4A CN107171354B (zh) | 2017-06-19 | 2017-06-19 | 柔性直流电网换流站下垂系数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107171354A true CN107171354A (zh) | 2017-09-15 |
CN107171354B CN107171354B (zh) | 2020-01-14 |
Family
ID=59818904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710466233.4A Active CN107171354B (zh) | 2017-06-19 | 2017-06-19 | 柔性直流电网换流站下垂系数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107171354B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109120005A (zh) * | 2018-06-22 | 2019-01-01 | 华北电力大学(保定) | 一种多端柔性直流输电系统功率协调控制方法 |
CN110635485A (zh) * | 2019-10-09 | 2019-12-31 | 湖南大学 | 针对多端直流系统自适应下垂控制的潮流模型计算方法 |
CN110675276A (zh) * | 2019-09-06 | 2020-01-10 | 中国电力科学研究院有限公司 | 一种用于直流输电系统反演下垂控制的方法及系统 |
CN113036798A (zh) * | 2021-03-11 | 2021-06-25 | 天津大学 | 一种柔性直流配电网n-1安全域求解方法 |
EP3711134B1 (en) | 2017-11-14 | 2022-01-05 | Hitachi Energy Switzerland AG | Voltage droop-based method in a power transmission system |
CN115001044A (zh) * | 2022-07-19 | 2022-09-02 | 东南大学溧阳研究院 | 一种换流器自主功率共享控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130021829A1 (en) * | 2011-07-22 | 2013-01-24 | Sao Charles | Arrangement and a method for supplying electric power |
CN105680465A (zh) * | 2016-01-26 | 2016-06-15 | 上海交通大学 | 用于多端柔直输电系统的改进协调控制方法 |
CN105896517A (zh) * | 2014-12-31 | 2016-08-24 | 国家电网公司 | 一种直流电网的电压下垂控制方法 |
CN105978016A (zh) * | 2016-06-30 | 2016-09-28 | 东北电力大学 | 一种基于最优潮流的多端柔性直流输电系统优化控制方法 |
CN106786715A (zh) * | 2016-12-08 | 2017-05-31 | 四川大学 | 一种多端vsc‑hvdc系统下垂控制系数确定方法 |
-
2017
- 2017-06-19 CN CN201710466233.4A patent/CN107171354B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130021829A1 (en) * | 2011-07-22 | 2013-01-24 | Sao Charles | Arrangement and a method for supplying electric power |
CN105896517A (zh) * | 2014-12-31 | 2016-08-24 | 国家电网公司 | 一种直流电网的电压下垂控制方法 |
CN105680465A (zh) * | 2016-01-26 | 2016-06-15 | 上海交通大学 | 用于多端柔直输电系统的改进协调控制方法 |
CN105978016A (zh) * | 2016-06-30 | 2016-09-28 | 东北电力大学 | 一种基于最优潮流的多端柔性直流输电系统优化控制方法 |
CN106786715A (zh) * | 2016-12-08 | 2017-05-31 | 四川大学 | 一种多端vsc‑hvdc系统下垂控制系数确定方法 |
Non-Patent Citations (1)
Title |
---|
刘瑜超等: "基于自适应下垂系数调节的VSC-MTDC功率协调控制", 《中国电机工程学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3711134B1 (en) | 2017-11-14 | 2022-01-05 | Hitachi Energy Switzerland AG | Voltage droop-based method in a power transmission system |
CN109120005A (zh) * | 2018-06-22 | 2019-01-01 | 华北电力大学(保定) | 一种多端柔性直流输电系统功率协调控制方法 |
CN109120005B (zh) * | 2018-06-22 | 2022-04-01 | 华北电力大学(保定) | 一种多端柔性直流输电系统功率协调控制方法 |
CN110675276A (zh) * | 2019-09-06 | 2020-01-10 | 中国电力科学研究院有限公司 | 一种用于直流输电系统反演下垂控制的方法及系统 |
CN110635485A (zh) * | 2019-10-09 | 2019-12-31 | 湖南大学 | 针对多端直流系统自适应下垂控制的潮流模型计算方法 |
CN113036798A (zh) * | 2021-03-11 | 2021-06-25 | 天津大学 | 一种柔性直流配电网n-1安全域求解方法 |
CN113036798B (zh) * | 2021-03-11 | 2023-11-07 | 天津大学 | 一种柔性直流配电网n-1安全域求解方法 |
CN115001044A (zh) * | 2022-07-19 | 2022-09-02 | 东南大学溧阳研究院 | 一种换流器自主功率共享控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107171354B (zh) | 2020-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107171354A (zh) | 柔性直流电网换流站下垂系数计算方法 | |
AU2018101070A4 (en) | Automatic voltage control method, device and system for wind farm | |
CN109390962B (zh) | 一种多端柔直系统的不平衡功率自适应优化分配方法 | |
Hui et al. | A robust control strategy to improve transient stability for AC-DC interconnected power system with wind farms | |
CN109638870A (zh) | 一种特高压直流送端电网的调相机配置方法 | |
CN109787297A (zh) | 一种考虑暂态特性的特高压直流送端电网无功优化方法 | |
CN114819693A (zh) | 一种多直流馈入电网的新能源承载能力评估方法和装置 | |
CN116365506A (zh) | 含分布式光伏主动配电网节能降损优化方法及系统 | |
CN105406488B (zh) | 一种基于光伏逆变器无功调节的过电压抑制方法 | |
CN105470978B (zh) | 一种静止无功补偿装置成组协调控制方法 | |
CN106786629A (zh) | 一种风场内部无功电压协调控制方法 | |
Li et al. | Research on capacity planning of renewable energy grid integration based on effective short circuit ratio | |
CN111614110A (zh) | 一种基于改进多目标粒子群算法的受端电网储能优化配置方法 | |
CN108306300B (zh) | 一种用于确定主网断面上统一潮流控制器容量的方法及系统 | |
CN110323779B (zh) | 一种分布式发电与储能装置的功率动态聚合的方法和系统 | |
CN115313397B (zh) | 一种基于智能软开关的电压就地控制方法 | |
CN109245090A (zh) | 一种大功率缺失下频率最低点预测的解析模型的建模方法 | |
Ye et al. | Improved droop control strategy for an MMC-MTDC connected to offshore wind farms with dynamic correction of the actual operating point | |
CN112199822B (zh) | 一种用于寻求最优系统阻抗值的外点罚函数方法 | |
CN111682584B (zh) | 一种基于电压调节器和多光伏系统的配电网电压管理方法 | |
CN105207249A (zh) | 一种风电场恒定功率因数控制系统及其控制方法 | |
Mengxi et al. | Reactive power coordinated control strategy based on PSO for wind farms cluster | |
Zhan et al. | Coordinating Operation and Control Strategy for the Integrated Wind-Storage VSC-MTDC System | |
CN117709689B (zh) | 一种计及整体效能与能量阻抗的风电场功率分配优化方法 | |
CN114204537B (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 |