CN106877363B - 一种静止同步串联补偿器抑制系统次同步振荡方法及装置 - Google Patents
一种静止同步串联补偿器抑制系统次同步振荡方法及装置 Download PDFInfo
- Publication number
- CN106877363B CN106877363B CN201710099220.8A CN201710099220A CN106877363B CN 106877363 B CN106877363 B CN 106877363B CN 201710099220 A CN201710099220 A CN 201710099220A CN 106877363 B CN106877363 B CN 106877363B
- Authority
- CN
- China
- Prior art keywords
- sssc
- damping
- parameters
- population
- control
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 46
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000010355 oscillation Effects 0.000 title claims abstract description 18
- 230000000452 restraining effect Effects 0.000 title abstract description 7
- 238000013016 damping Methods 0.000 claims abstract description 97
- 238000012546 transfer Methods 0.000 claims abstract description 37
- 210000000349 chromosome Anatomy 0.000 claims description 36
- 238000002347 injection Methods 0.000 claims description 25
- 239000007924 injection Substances 0.000 claims description 25
- 238000004422 calculation algorithm Methods 0.000 claims description 18
- 238000009434 installation Methods 0.000 claims description 18
- 230000002068 genetic effect Effects 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 10
- 239000000243 solution Substances 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 7
- 230000001427 coherent effect Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000013461 design Methods 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 49
- 230000001052 transient effect Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 239000003990 capacitor Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010845 search algorithm Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method Methods 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/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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
- H02J3/1807—Arrangements for adjusting, eliminating or compensating reactive power in networks using series compensators
-
- 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/10—Flexible AC transmission systems [FACTS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明涉及一种静止同步串联补偿器抑制系统次同步振荡方法及装置，所述方法包括：根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；将所述SSSC装置安装点的电压的频率偏差和有功潮流作为输入量，构建阻尼控制器的传递函数，其中，将所述SSSC装置安装点的电压的频率偏差和有功潮流作为阻尼控制器的输入量；优化所述阻尼控制器的传递函数的控制参数；本发明提供的技术方案，能够提高电网运行阻尼性和稳定性，同时，使得静止同步串联补偿器所在区域的电网阻尼特性达到最优，又适用于实际控制器的设计。
Description
技术领域
本发明涉及柔性交流输电技术领域，具体涉及一种静止同步串联补偿器抑制系统次同步振荡方法及装置。
背景技术
近年来电力系统快速发展，潮流需求增大，为了提高系统的输电能力串并联形式的无功补偿装置得到广泛应用。静止同步串联补偿器(SSSC)作为一种新型串联无功补偿装置，具有控制灵活、结构简单、快速有效调控系统潮流并可以一定程度抑制次同步振荡(SSR)的优点。SSSC不存在可能与交流系统电感参数发生次同步振荡的电容量，减少了次同步振荡发生机率。为确定含SSSC系统的次同步振荡问题，复转矩分析法和特征值分析法基于系统小扰动下建立微分方程，物理意义明确。虽然基于电力电子器件而建立的时域分析法具有更准确地判定依据，但是不能直接得到阻尼和频率的关系，不利于从物理意义上判断。
目前对于含SSSC系统的SSR问题和阻尼控制器的设计进行的研究方法中有的从系统角度仅考虑幅值和调制比作为输入时的单机无穷大系统模型；有的方法基于测试信号法提出SSR的多模式阻尼控制参数设计；有的方法是从经验的角度设计控制参数；这些方式将考虑SSR主要因素，实际上简化了问题，不一定能找到阻尼控制参数的最优控制量。
发明内容
本发明提供一种静止同步串联补偿器抑制系统次同步振荡方法及装置，其目的是提高电网运行阻尼性和稳定性，同时，使得静止同步串联补偿器所在区域的电网阻尼特性达到最优，又适用于实际控制器的设计。
本发明的目的是采用下述技术方案实现的：
一种静止同步串联补偿器抑制系统次同步振荡方法，其改进之处在于，包括：
根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；
根据含SSSC的发电机系统模型将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；
构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
优化所述阻尼控制器的传递函数的控制参数。
优选的，所述含SSSC的单机无穷大系统的系统线性化模型为：
上式中，δ为系统功角，ω为角速度，E′q为交轴暂态电动势，Eqe为强制空载电动势，TJ为惯性时间常数，T′d0为发电机暂态时间常数，TA和KA均为励磁系统参数，δsc为SSSC输出电压和交流线路电流的夹角，m为调制比，k1至k12均为系统线性化模型参数，按下式确定系统线性化模型参数k1至k12：
上式中，U为接入系统电压，kT为变压器的变比，xd为发电机d轴电抗，x′d为发电机d轴暂态电抗，Udc为直流电容电压，UG为端口电压，UGd为发电机d轴端口电压，UGq为发电机q轴端口电压，xq为发电机q轴电抗，xqΣ、xdΣ和x′dΣ如下式：
优选的，按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解所述含SSSC的发电机系统模型，获取阻尼的控制量为：
上式中，k为控制增益，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
优选的，按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，Ti(i＝1～4)为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量。
优选的，所述优化所述阻尼控制器的传递函数的控制参数，包括：
a.确定量子染色体的编码方式，设定所述灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤a和步骤c中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
进一步的，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，G1和G2为增益和时间常数；T1和T3作为待优化变量。
一种静止同步串联补偿器抑制系统次同步振荡装置，其改进之处在于，所述装置包括：
第一分析模块，用于根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；
第二分析模块，用于根据含SSSC的发电机系统模型，将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；
构建模块，用于构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
优化模块，用于优化所述阻尼控制器的传递函数的控制参数。
优选的，所述含SSSC的单机无穷大系统的系统线性化模型为：
上式中，δ为系统功角，ω为角速度，E′q为交轴暂态电动势，Eqe为强制空载电动势，TJ为惯性时间常数，T′d0为发电机暂态时间常数，TA和KA均为励磁系统参数，δsc为SSSC输出电压和交流线路电流的夹角，m为调制比，k1至k12均为系统线性化模型参数，按下式确定系统线性化模型参数k1至k12：
上式中，U为接入系统电压，kT为变压器的变比，xd为发电机d轴电抗，x′d为发电机d轴暂态电抗，Udc为直流电容电压，UG为端口电压，UGd为发电机d轴端口电压，UGq为发电机q轴端口电压，xq为发电机q轴电抗，xqΣ、xdΣ和x′dΣ如下式：
优选的，按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解单元，用于求解所述含SSSC的发电机系统模型，获取阻尼的控制量为：
上式中，k为控制增益，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
优选的，所述构建模块，包括：
按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，Ti(i＝1～4)为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量。
优选的，所述优化模块，包括：
a.确定量子染色体的编码方式，设定所述灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤1和步骤3中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
进一步的，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，G1和G2为增益和时间常数；T1和T3作为待优化变量。
本发明的有益效果：
本发明提供的技术方案，模型上不仅考虑交流系统，还结合SSSC的动态特性，改进模型更好的跟随电网特性，使控制设计时的稳定性水平有所提高，阻尼控制器的设计考虑了实际的采样问题，结合测量点提取阻尼控制器中的主要控制量进行优化，使用不同权重的量，更贴近电网的实控性，进一步的，采用量子遗传算法优化阻尼控制参数，其收敛速度快，全局寻优能力强，更快速的选取了阻尼控制器的最优参数。
附图说明
图1是本发明一种静止同步串联补偿器抑制系统次同步振荡方法的流程图；
图2是本发明实施例中单机无穷大系统示意图；
图3是本发明实施例中SSSC阻尼控制策略；
图4是本发明一种静止同步串联补偿器抑制系统次同步振荡装置的结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作详细说明。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
本发明提供的一种静止同步串联补偿器抑制系统次同步振荡方法，如图1所示，包括：
101.根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；
102.根据含SSSC的发电机系统模型将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；
103.构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
104.优化所述阻尼控制器的传递函数的控制参数。
具体的，现有修正模型不包含SSSC暂态特性分析，输入量简单，带来模型准确度低和系统描述不全面的缺点。本发明分析SSSC暂态特性，增加系统状态方程数，列写线性微分方程，分析小扰动下含SSSC的非线性系统对SSR的免疫能力。含SSSC的单机无穷大系统的等效示意图如图2所示。其中，发电机多质块轴系包含：励磁机质块(EXC)，发电机质块(GEN)，低压缸B质块(LPB)，低压缸A质块(LPA)，中压缸质块(MP)和高压缸质块(HP)，所述含SSSC的单机无穷大系统的系统线性化模型为：
上式中，δ为系统功角，ω为角速度，E′q为交轴暂态电动势，Eqe为强制空载电动势，TJ为惯性时间常数，T′d0为发电机暂态时间常数，TA和KA均为励磁系统参数，δsc为SSSC输出电压和交流线路电流的夹角，m为调制比，k1至k12均为系统线性化模型参数，按下式确定系统线性化模型参数k1至k12：
上式中，U为接入系统电压，kT为变压器的变比，xd为发电机d轴电抗，x′d为发电机d轴暂态电抗，Udc为直流电容电压，UG为端口电压，UGd为发电机d轴端口电压，UGq为发电机q轴端口电压，xq为发电机q轴电抗，xqΣ、xdΣ和x′dΣ如下式：
通过上述含SSSC的单机无穷大系统的系统线性化模型可以看出，SSSC影响发电机的参数为SSSC输出电压和交流线路电流的夹角和SSSC的调制比，因此，所述步骤102需根据含SSSC的发电机系统模型确定影响SSSC的调制比的参数，包括：
按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解所述含SSSC的发电机系统模型，可以根据含SSSC的发电机系统模型和二次性能指标可知控制量的计算表达式：
定义Riccati方程：
式中P为满足上述方程的解，可解得SSR的控制量；
上式中，k为控制增益，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
电压质量是电网首要指标，其对应权重越大对电压约束越强，动态调整电压时精度越高。而动态或振荡过程对频率约束越强振荡次数和幅值越小。
SSSC暂态特性取决于直流电容电压的动态特性，直流电压恒定是SSSC稳定运行的关键。现有文献详述了直流电压控制为系统提供负阻尼转矩，增益越大提供负阻尼越多，而注入电压控制提供正阻尼转矩，且增益越大提供正阻尼越多。考虑动态特性时结合了正负阻尼特性，可在协调配合下达到整体阻尼最优以及直流电容电压和系统潮流的稳定。
基于阻尼的控制量的表达式可以看出，影响SSSC的调制比的参数为SSSC装置安装点的电压的频率偏差和有功潮流，则在所述步骤103中需根据影响SSSC的调制比的参数构建阻尼控制器的传递函数，按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，Ti(i＝1～4)为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量，结合SSSC暂态特性和设计最优频率和最优潮流控制的阻尼控制策略如图3所示。
获取阻尼控制器的传递函数之后，需优化阻尼控制器的传递函数的参数，目前，常用的控制器参数确定方法有相位补偿法和智能算法。前者基于频域，控制器顺序和个数约束控制的准确性；后者对所有控制器参数统一协调，只基于目标函数和约束条件即可一次性确定所有控制器参数，相比前者更加快捷有效。国内外对于FACTS装置的SSR研究中使用较多的有遗传算法、粒子群优化和“教与学”优化算法。遗传算法在抑制SSR的阻尼控制器参数设计时考虑多目标协调，自适应能力强，收敛速度快，优化精度高，具有很大优势，研究对自适应能力要求较高的SSR问题时可优先考虑遗传算法或其改进算法。近期国际上出现了引力搜索算法用于确定大网间协调配合控制器的参数[2]，并对故障下各类算法的阻尼进行了对比，发现引力搜索算法和遗传算法两者的差别不大。本发明提出结合系数矩阵的特征根和寻优算法来确定控制参数的方法。量子遗传算法结合了量子计算和遗传算法的优势，用其来分析特征根进而优化控制参数，所述步骤104，包括：
a.确定量子染色体的编码方式，设定所述灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤1和步骤3中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
其中，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，Tf和Tp与T2和T4一样视为常数，取[5,5,0.05,0.05]，G1和G2为增益和时间常数；T1和T3作为待优化变量。
一种静止同步串联补偿器抑制系统次同步振荡装置，如图4所示，所述装置包括：
第一分析模块，用于根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；
第二分析模块，用于根据含SSSC的发电机系统模型，将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；
构建模块，用于构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
优化模块，用于优化所述阻尼控制器的传递函数的控制参数。
具体的，所述含SSSC的单机无穷大系统的系统线性化模型为：
上式中，δ为系统功角，ω为角速度，E′q为交轴暂态电动势，Eqe为强制空载电动势，TJ为惯性时间常数，T′d0为发电机暂态时间常数，TA和KA均为励磁系统参数，δsc为SSSC输出电压和交流线路电流的夹角，m为调制比，k1至k12均为系统线性化模型参数，按下式确定系统线性化模型参数k1至k12：
上式中，U为接入系统电压，kT为变压器的变比，xd为发电机d轴电抗，x′d为发电机d轴暂态电抗，Udc为直流电容电压，UG为端口电压，UGd为发电机d轴端口电压，UGq为发电机q轴端口电压，xq为发电机q轴电抗，xqΣ、xdΣ和x′dΣ如下式：
所述第二确定模块，包括：
第一确定单元，用于按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解单元，用于求解所述含SSSC的发电机系统模型，获取阻尼的控制量为：
上式中，k为控制增益，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
所述构建模块，包括：
按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，Ti(i＝1～4)为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量。
所述优化模块，包括：
a.确定量子染色体的编码方式，设定所述灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤1和步骤3中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
其中，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，Tf和Tp与T2和T4一样视为常数，取[5,5,0.05,0.05]，G1和G2为增益和时间常数；T1和T3作为待优化变量。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
Claims (8)
1.一种静止同步串联补偿器抑制系统次同步振荡方法，其特征在于，所述方法包括：
根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；根据含SSSC的发电机系统模型，将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
优化所述阻尼控制器的传递函数的控制参数；
按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解所述含SSSC的发电机系统模型，获取阻尼的控制量为：
上式中，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
2.如权利要求1所述的方法，其特征在于，按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，T1、T2、T3和T4为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量。
3.如权利要求1所述的方法，其特征在于，所述优化所述阻尼控制器的传递函数的控制参数，包括：
a.确定量子染色体的编码方式，设定灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤a和步骤c中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
4.如权利要求3所述的方法，其特征在于，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，α和β为控制目标权系数，G1和G2为控制目标增益；T1和T3作为待优化变量。
5.一种静止同步串联补偿器抑制系统次同步振荡装置，其特征在于，所述装置包括：
第一分析模块，用于根据含SSSC的单机无穷大系统的系统线性化模型，将SSSC输出电压与交流线路电流的夹角和SSSC的调制比作为SSSC影响发电机的参数；
第二分析模块，用于根据含SSSC的发电机系统模型，将SSSC装置安装点的电压的频率偏差和有功潮流作为影响SSSC的调制比的参数；
构建模块，用于构建阻尼控制器的传递函数，其中，所述阻尼控制器的输入量包括：所述SSSC装置安装点的电压的频率偏差和有功潮流；
优化模块，用于优化所述阻尼控制器的传递函数的控制参数；
按下式确定含SSSC的发电机系统模型：
上式中，Y为输出量，X为状态变量，U为控制量，A1和B1均为发电机修正模型参数矩阵，ΔPe为SSSC装置安装点的线路潮流中的有功功率，Δf为SSSC装置安装点的电压的频率偏差，ΔU为SSSC装置安装点的电压，P0为潮流平均值，f0为频率平均值；
求解单元，用于求解所述含SSSC的发电机系统模型，获取阻尼的控制量为：
上式中，kmp、kmf、kmu、kfp、kff和kfu为控制量系数。
6.如权利要求5所述的装置，其特征在于，所述构建模块，包括：
按下式构建阻尼控制器的传递函数f(s)：
上式中，α和β为控制目标权系数，α+β＝1，T1、T2、T3和T4为待定系数，Tf和Tp为隔直时间常数，G1和G2为控制目标增益，s为复变量。
7.如权利要求5所述的装置，其特征在于，所述优化模块，包括：
a.确定量子染色体的编码方式，设定灾变变速量子遗传算法的初始参数，对不包含有静止同步串联补偿器的电力系统进行阻尼量计算；
所述初始参数包括种群规模popsize、最大进化代数genmax、染色体长度lchrom、染色体变异概率Pm和最优解不变时发生灾变的代数pregenM；
b.对所述种群中各个量子染色体进行坍缩测量得到染色体二进制串；将所述染色体二进制串解码为两个十进制数，则所述两个十进制数分别为静止同步串联补偿器的安装位置和注入电压幅值；
c.依据所述安装位置和注入电压幅值对包含有所述静止同步串联补偿器的电力系统进行阻尼量计算；依据所述步骤a和步骤c中阻尼量计算的结果计算目标函数的函数值；
d.比较所述种群中每个个体目标函数的函数值，记录该种群的最优个体；
e.对所述种群依次进行量子门更新、量子相干交叉和量子变异，并判断是否对所述种群进行量子灾变；
f.判断所述种群的进化代数是否达到最大进化代数：
若达到，则步骤d得到的所述安装位置和注入电压幅值为最优的静止同步串联补偿器安装位置和注入电压幅值；
若没有达到，则返回步骤b。
8.如权利要求7所述的装置，其特征在于，按下式确定所述阻尼控制器的传递函数的控制参数的目标函数：
上式中，F为阻尼比，δi为特征根实部，ωi为特征根虚部，n为特征根数目；
按下式确定所述阻尼控制器的传递函数的控制参数的约束条件：
上式中，α和β为控制目标权系数，G1和G2为控制目标增益；T1和T3作为待优化变量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710099220.8A CN106877363B (zh) | 2017-02-23 | 2017-02-23 | 一种静止同步串联补偿器抑制系统次同步振荡方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710099220.8A CN106877363B (zh) | 2017-02-23 | 2017-02-23 | 一种静止同步串联补偿器抑制系统次同步振荡方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106877363A CN106877363A (zh) | 2017-06-20 |
| CN106877363B true CN106877363B (zh) | 2020-09-04 |
Family
ID=59168644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710099220.8A Active CN106877363B (zh) | 2017-02-23 | 2017-02-23 | 一种静止同步串联补偿器抑制系统次同步振荡方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106877363B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108011372B (zh) * | 2017-02-22 | 2021-04-16 | 全球能源互联网研究院 | 一种静止同步串联补偿器的机电暂态建模方法及装置 |
| CN108134398B (zh) * | 2018-01-04 | 2020-10-27 | 全球能源互联网研究院 | 基于电流信号反馈的抑制火电机组次同步振荡方法 |
| CN108964062B (zh) * | 2018-08-17 | 2022-03-04 | 武汉理工大学 | 一种确定分布式潮流控制器3次谐波电流取值范围的方法 |
| CN109445307A (zh) * | 2018-11-15 | 2019-03-08 | 全球能源互联网研究院有限公司 | 一种次同步振荡抑制装置控制保护系统的测试系统及方法 |
| CN109659955B (zh) * | 2018-12-25 | 2022-05-10 | 国网黑龙江省电力有限公司电力科学研究院 | 一种dfig与sssc协调电力系统广域阻尼控制方法 |
| CN110148952B (zh) * | 2019-03-15 | 2022-12-13 | 湖南工业大学 | 一种静止同步补偿器参数优化控制方法 |
| CN111262258A (zh) * | 2020-03-19 | 2020-06-09 | 国网天津市电力公司电力科学研究院 | 用于多端柔性直流系统的附加阻尼内环控制方法及控制器 |
| CN111509704B (zh) * | 2020-04-23 | 2022-07-19 | 郑州轻工业大学 | 基于sssc的含dfig-sofc多能源系统动态交互分析方法 |
| CN111965421B (zh) * | 2020-08-17 | 2021-08-24 | 广东电网有限责任公司电力科学研究院 | 一种基于频率耦合性的次超同步间谐波检测方法及装置 |
| CN113094899A (zh) * | 2021-04-07 | 2021-07-09 | 全球能源互联网研究院有限公司 | 一种随机潮流计算方法、装置、电子设备及存储介质 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5374853A (en) * | 1991-12-13 | 1994-12-20 | Electric Power Research, Inc. | Transient damping thyristor controlled series capacitor system |
| EP1980006A1 (en) * | 2006-02-01 | 2008-10-15 | ABB Technology Ltd | Controlled series compensator and method |
| US20100133915A1 (en) * | 2006-05-30 | 2010-06-03 | Abb Research Ltd | Thyristor controllied series capacitor adapted to damp sub synchronous resonances |
| CN101183790A (zh) * | 2007-12-19 | 2008-05-21 | 湖南大学 | 基于statcom与svc的电能质量调节器协调控制方法 |
| CN101834448B (zh) * | 2010-03-23 | 2012-03-28 | 浙江大学 | 一种基于静止同步串联补偿器的抑制电力系统次同步振荡的方法 |
| CN102624012B (zh) * | 2012-03-26 | 2014-10-22 | 中国电力科学研究院 | 基于力矩分解的积分算法识别负阻尼低频振荡的方法 |
| CN105449695B (zh) * | 2015-11-18 | 2018-03-06 | 华北电力科学研究院有限责任公司 | 一种发电机阻尼控制的实现方法 |
| CN105262118A (zh) * | 2015-11-20 | 2016-01-20 | 江苏省电力公司电力经济技术研究院 | 基于statcom的次同步振荡抑制方法及statcom的控制装置 |
| CN105633980B (zh) * | 2015-12-29 | 2018-01-05 | 贵州理工学院 | 一种sssc新型非线性鲁棒控制系统及其控制方法 |
| CN105811429B (zh) * | 2016-03-14 | 2021-01-15 | 全球能源互联网研究院 | 一种静止同步串联补偿器的选址定容方法 |
-
2017
- 2017-02-23 CN CN201710099220.8A patent/CN106877363B/zh active Active
Non-Patent Citations (2)
| Title |
|---|
| A Comparative Study of Damping Subsynchronous Resonance Using TCSC and IMDU;Prakash Chittora et al.;《IEEE》;20121231 * |
| Application of SSSC to Improve the Transmission Capability of Tianjin 220kV Power Grid;LU Zhen-gang et al.;《2016 China International Conference on Electricity Distribution (CICED 2016)》;20160813 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106877363A (zh) | 2017-06-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106877363B (zh) | 一种静止同步串联补偿器抑制系统次同步振荡方法及装置 | |
| CN107800146B (zh) | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 | |
| CN108631331B (zh) | 一种双馈风电场次同步振荡抑制方法及装置 | |
| CN109962495B (zh) | 一种超低频振荡扰动源定位及抑制方法 | |
| CN102130454B (zh) | 基于计算机辅助设计的电力系统动态稳定控制方法及系统 | |
| CN108631332B (zh) | 双馈风电场次同步振荡svc抑制方法及装置 | |
| CN110380429B (zh) | 一种兼顾阻尼特性的水电机组一次调频多目标优化方法 | |
| CN108964061B (zh) | 一种考虑负荷频率电压静特性的含风电交直流电力系统概率动态连续潮流计算方法 | |
| CN113883008B (zh) | 一种可抑制多扰动因素的风机模糊自适应变桨距控制方法 | |
| Chen et al. | Indirect model reference adaptive control for a class of linear fractional order systems | |
| CN110323741B (zh) | 基于录波曲线的电力系统低频振荡分析与抑制方法 | |
| CN107402515A (zh) | 一种优化同步发电机励磁系统pid参数的方法 | |
| CN110262222A (zh) | 一种区间系统最优区间pid控制器参数优化方法 | |
| CN108123434B (zh) | 一种计算pv曲线斜率以求取pv曲线运行点的方法 | |
| CN104795812A (zh) | 一种考虑定子阻尼效应并采用变量代换的慢同调分区法 | |
| CN113036769A (zh) | 一种电力系统静态电压稳定分析方法及系统 | |
| CN114997505A (zh) | 一种经变流器并网电源的稳定性分析方法及系统 | |
| Wang et al. | Design of fractional order pi controller for permanent magnet synchronous motor | |
| Sadhana et al. | Revamped Sine Cosine Algorithm Centered Optimization of System Stabilizers and Oscillation Dampers for Wind Penetrated Power System | |
| Abbasi et al. | Application of an additive self-tuning controller for static synchronous series compensator for damping of sub-synchronous resonance oscillations | |
| Kanchanaharuthai | Nonlinear controller design for hydraulic turbine regulating systems via immersion and invariance | |
| CN112968447A (zh) | 一种衡量电压测量时间常数对机组动态稳定影响的方法 | |
| CN111237116A (zh) | 非线性隐极水轮发电机水门与励磁联合调节方法和系统 | |
| CN104408330A (zh) | 水轮机调节系统高阶数学模型的降阶方法 | |
| CN116088296B (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 |