CN106684866A - 一种静态开环综合增益的计算方法及装置 - Google Patents
一种静态开环综合增益的计算方法及装置 Download PDFInfo
- Publication number
- CN106684866A CN106684866A CN201710014788.5A CN201710014788A CN106684866A CN 106684866 A CN106684866 A CN 106684866A CN 201710014788 A CN201710014788 A CN 201710014788A CN 106684866 A CN106684866 A CN 106684866A
- Authority
- CN
- China
- Prior art keywords
- excitation
- gain
- voltage
- exciting current
- open loop
- 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.)
- Pending
Links
- 230000003068 static effect Effects 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000005284 excitation Effects 0.000 claims abstract description 161
- 238000005259 measurement Methods 0.000 claims abstract description 54
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 230000010354 integration Effects 0.000 claims description 69
- 230000010363 phase shift Effects 0.000 claims description 26
- 238000000205 computational method Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 10
- 230000005611 electricity Effects 0.000 claims description 8
- 230000009467 reduction Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 17
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000010248 power generation Methods 0.000 abstract 1
- 238000004804 winding Methods 0.000 description 31
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 12
- 229910052710 silicon Inorganic materials 0.000 description 12
- 239000010703 silicon Substances 0.000 description 12
- 238000004088 simulation Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 9
- 210000004027 cell Anatomy 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000003786 synthesis reaction Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000003321 amplification Effects 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000036772 blood pressure Effects 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 206010068052 Mosaicism Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Tests Of Circuit Breakers, Generators, And Electric Motors (AREA)
Abstract
本发明的实施例提供一种静态开环综合增益的计算方法及装置，涉及发电及电力系统建模技术领域，能够提高发电机自并励静止励磁系统的静态开环综合增益实测计算结果的准确性。具体方案包括：通过实测发电机空载特性确定励磁电流基准值；采用外部恒压源供电进行空载升压试验测得机端参考电压、机端电压、励磁电流；对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统静态开环综合增益。本发明用于计算静态开环综合增益。
Description
技术领域
本发明的实施例涉及发电及电力系统建模技术领域，尤其涉及一种静态开环综合增益的计算方法及装置。
背景技术
我国电网大规模互联以及交直流联合运行过程中大扰动功角稳定以及区域弱阻尼功率振荡小扰动稳定性问题凸显，促使电网公司要求建立更为精确和有效的发电机励磁和调速系统模型，通过实测确定准确的模型参数，为开展电网动态分析提供准确的基础数据。
自1980年后IEEE先后多次发布用于电力系统稳定分析的典型励磁系统标准模型，包括直流励磁机励磁系统、交流励磁系统、静止可控硅励磁系统等。这些励磁模型都是基于对典型励磁系统物理结构各环节的详细分析，然后经过合理的等值和简化，模型的详细程度能够满足电力系统稳定性问题(0～3Hz)研究。
IEEE推荐的励磁模型中，大多数参数可以明确对应到实际的物理环节，从而可以通过实测直接得到，比如励磁机时间常数，励磁机饱和特性，各类校正控制器的增益和时间常数及其限幅值等。而有的参数则是对多个环节特性进行综合等值得到，单独对某一个物理环节测试不能确定该参数，励磁系统的静态开环综合增益就是这种类型的参数，其数值与励磁控制器增益整定、励磁功率单元放大倍数以及励磁系统基准值选取相关，实测建模时需要对相关环节进行仔细分析，制定合适的实测方案。
目前实测建模中存在以下两方面问题，导致静态开环综合增益计算结果的准确性欠佳，具体说明如下：
1)励磁系统建模导则中定义励磁系统放大倍数作为静态开环综合增益，该放大倍数具体是励磁功率整流输出相对励磁控制单元输入即电压偏差的放大倍数。
目前实测计算中一般采用励磁电压、励磁电流、换相电抗来计算整流器空载输出电压，从而得到静态开环综合增益，该方法忽略了整流换相压降的影响；
2)采用额定励磁电压除以额定励磁电流作为励磁电阻基准值，认为大多数需要模拟的工况对应的励磁绕组电阻与发电机额定运行状态下的励磁绕组电阻相近，但实测状态下励磁绕组电阻很可能偏离额定状况较多，将会导致较大的励磁电压基准值偏差，导致静态开环综合增益计算结果出现偏差。
发明内容
本发明的实施例提供一种静态开环综合增益的计算方法及装置，能够提高发电机自并励静止励磁系统的静态开环综合增益实测计算结果的准确性。为了达成上述目的，本申请的实施例采用如下技术方案：
第一方面，提供一种静态开环综合增益的计算方法，包括：
通过实测发电机空载特性确定励磁电流基准值；
采用外部恒压源供电进行空载升压试验测得机端参考电压、机端电压、励磁电流；
对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统的静态开环综合增益。
第二方面，提供一种静态开环综合增益的计算装置，用于执行第一方面所提供的方法。
静态开环综合增益大小不仅影响到励磁系统电压调整动态特性，也决定了发电机静态电压调节精度(即电压调差率)，一般要求其标幺值在200～400，以满足相关标准中电压静差率为1％的要求。另外，该增益是影响空载电压阶跃响应的重要参数，对电力系统动态品质有直接影响。为了准确地确定励磁系统开环综合增益，首先必须从标准模型出发准确地理解综合增益的涵义，明确与其直接相关的励磁系统物理环节；另一方面，不同的励磁装置厂家在励磁控制参数整定方面采用了各自不同原则和方法，而制定实测方案时，实测条件可能与实际运行的情况不一样，必须仔细了解相关励磁系统的控制参数整定原则，从而根据实测方案提出正确的修正方法。本发明的实施例所提供的静态开环综合增益的计算方法及装置，既能考虑换相压降的影响，又避免励磁绕组基准值选取困难问题，从而提高了静态开环综合增益实测计算结果的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的实施例中静态开环综合增益的计算方法流程示意图；
图2(a)为空载发电机自并励励磁调节系统结构的示意图；
图2(b)为空载发电机自并励励磁调节系统数学模型的示意图；
图3为洪家渡#3机组励磁系统结构和相关参数图；
图4为洪家渡#3机空载特性和励磁电流基准确定示意图；
图5为洪家渡#3机励磁系统综合开环增益求取的直线拟合结果；
图6为洪家渡#3机空载5％阶跃仿真曲线与实测曲线对比；
图7为昌江核电#2机组励磁系统结构和相关参数；
图8为昌江核电#2机励磁电流基准确定示意图；
图9为昌江核电#2机励磁系统综合开环增益求取的直线拟合结果；
图10为昌江核电#2机空载5％阶跃仿真曲线与实测曲线对比；
图11为本发明的实施例中静态开环综合增益的计算装置结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
本发明的实施例提供一种静态开环综合增益的计算方法，结合图1所示，包括以下步骤：
101、通过实测发电机空载特性确定励磁电流基准值。
具体的，发电机空载运行，机端电压从0.1p.u.逐渐升压然后再降压到0.1p.u，升压最大值应该不低于1.2p.u.。
记录升压和降压过程中各点的机端电压、励磁电流、励磁电压。对三相机端电压进行平均得到平均机端电压，绘制平均机端电压对应于励磁电流的上升和下降二条曲线，对二条曲线再进行机端电压值的平均化处理，得到空载曲线，然后平移到过零点，从零点作切线得到空载气隙线。
在气隙线上对应于额定机端电压的励磁电流，即为励磁电流基准值。
102、采用外部恒压源供电进行空载升压试验。
通过升压试验可测得机端参考电压、机端电压、励磁电流等。
103、对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统的静态开环综合增益。
正常运行时励磁系统的静态开环综合增益以K表示，则K＝KR×KA。其中，KR表示正常运行条件下励磁控制单元增益设定值，KA表示励磁功率单元增益。
采用完全余弦移相触发方式时，K'表示实测条件下励磁系统的静态开环综合增益，K'R表示实测时励磁控制单元增益的设定值。
当采用非完全余弦移相触发方式时，U2N、U2分别为励磁变二次侧额定励磁电压、实测时励磁变二次电压。
对K'的计算过程说明如下：
设置励磁控制单元为纯比例环节，采用外部恒电压源为励磁变提供电源，发电机空载运行，机端电压从0.4p.u.逐渐升压然后再降压，升压最大值可以在1.0p.u.上下。
记录升压和降压过程中各点的参考电压、机端电压、励磁电流。根据励磁电流基准值计算实测励磁电流的标幺值IFD，用额定机端电压分别计算实测记录参考电压的标幺值VREF和机端电压(可求三相平均值)的标幺值Vt，作出IFD与差值电压ΔV＝VREF-Vt的特性曲线，并采用最小二乘拟合法进行直线拟合，拟合直线的斜率即为实测条件下的静态开环综合增益K'。
结合上述步骤，确定实测条件下的静态开环综合增益K'后，实际计算静态开环综合增益的过程可用以下公式描述：
采用完全余弦移相触发方式时，有：
采用非完全余弦移相触发方式时，有：
正常运行时励磁系统静态开环综合增益为：
K＝KR×KA (3)
为使得本发明的发明目的、特征、优点能够更加的明显和易懂，下面将结合本发明实施例中的附图，本实施例后文对本申请方案的原理做说明，并提供具体实例。原理部分主要介绍本发明相关的技术理论以及依据；实例部分基于实际发电机实测数据，采用本发明方法和常规方法分别计算静态开环综合增益，并基于实测阶跃曲线比较两种方法的仿真效果。
第一部分 基本原理和方法
1静止自并励励磁系统开环综合增益
1.1静止自并励励磁系统结构与数学模型
图2(a)为采用静止自并励励磁系统的发电机励磁调节系统结构，包括励磁控制器、励磁功率单元以及空载发电机三个部分，图2(b)为分别对应上述三部分元件的系统数学模型。注意图2励磁系统模型省去了低励、过励以及限幅等非线性环节，这些环节与本发明方法无关。在图2中：
1)励磁控制器由超前滞后校正单元(TC1/TB1,TC2/TB2)和增益放大环节KR组成，VREF、Vt分别为机端电压给定值和测量值，VC为移相触发控制电平；
2)励磁功率单元包括移相触发控制单元和可控硅整流电路，移相触发控制单元根据移相触发控制电平VC计算移相触发角α，Vfd、EFD分别为可控硅整流器直流励磁电压输出的有名值与标幺值，Ifd为励磁电流有名值。励磁功率单元采用具有增益KA和时间常数TA的一阶惯性环节模拟；
3)T’d0为励磁绕组开路时间常数，对应的一阶惯性环节代表了不计饱和特性的发电机空载状态数学模型；
励磁系统模型参数一般都采用标幺值表示，因此建立模型首先需要确定励磁系统各变量的基准值。
1.2励磁系统标幺值基准选取方法
理论上来说，励磁系统基准值可以有不同基准值选取方式，只要模型参数归算正确，模拟得到的系统动态特性基本相同，但是不同基准值对应的模型参数却可以差别很大。
电力系统稳定计算程序中，由于励磁系统模型中励磁电压、励磁电流和励磁电阻的基准值与发电机模型中相应变量的基准值不一样，为此需要在两种不同基准值系统下对励磁电压、励磁电流和励磁电阻进行换算，称为励磁系统模型与发电机模型接口。励磁模型的基准值应该与所要应用计算程序中的模型接口规定一致。
目前国内外主流的电力系统稳定计算程序中，均按照IEEE推荐励磁系统模型的基准值定义以实现励磁模型与发电机模型接口。IEEE推荐的励磁系统标幺值基准选取方法如下：
1)选取发电机定子额定电压和额定电流分别为机端电压和电流的基准值，电压给定的基准值与机端电压基准值相同；
2)选取发电机空载特性气隙线上对应机端额定电压的励磁电流为励磁电流基准值IFDB，相应的励磁电压为励磁电压基准值EFDB，即
EFDB＝IFDBRFDB (4)
上式中，RFDB表示励磁绕组电阻基准值，应选取为指定工作温度t°c下的绕组电阻值Rfd@t℃，即对应于待建模工况励磁绕组实际电阻，有：
RFDB＝Rfd@t℃ (5)
IFDB可由发电机空载特性气隙线确定，RFDB可根据指定绕组工作温度和已知温度下的励磁绕组电阻(制造厂家通常给出25℃和125℃的电阻值)修正计算得到。
1.3励磁系统模型的开环综合增益定义
励磁系统模型的开环综合增益定义为励磁系统模型输出励磁电压标幺值与模型输入电压偏差标幺值之比值，即
另一方面，从图2(b)可见，励磁系统开环综合增益等于励磁控制器整定增益与励磁功率单元增益的乘积，即：
K＝KRKA (7)
励磁系统开环综合增益可以根据式(6)采用实测方法确定，也可以首先由励磁控制器整定参数得到KR，根据移相触发方式和励磁电压基准值计算出KA，进而根据式(7)计算开环综合增益。
无论采用上述哪种方式确定开环综合增益，必须注意模型输出励磁电压的基准值必须与稳定计算中励磁模型与发电机模型接口的相关规定一致，否则将影响到开环综合增益结果的正确性。
稳定计算中励磁模型输出通常都是采用IEEE推荐的基准值定义方法，励磁绕组电阻标幺值RFD≡1.0p.u.，当IFD＝1.0p.u.时，相应地在发电机空载特性气隙线上对应的机端电压为额定值；即在稳态情况下，EFD≡IFD。
IEEE也有推荐选取指定工作温度下(汽轮发电机组100℃，水轮发电机组75℃)的励磁绕组电阻为其基准值，是因为实际运行工况下励磁绕组温度与之比较接近，从而近似认为RFD≈1.0p.u.。如果励磁绕组实际工作温度偏离较多导致其电阻变化较大时，需要对RFDB进行修正，并按式(4)重新计算EFDB。
国内开展励磁实测建模按照额定励磁电压VfN与额定励磁电流IfN的比值确定RFDB，即
也是认为大多数需要模拟的工况对应的励磁绕组电阻与发电机额定运行状态下的励磁绕组电阻相近，近似认为RFD≈1.0p.u.。但是进行开环综合增益实测时，实测状态下的Rfd@t℃很可能偏离RFDB较多，则需要根据Rfd@t℃重新确定EFDB。
2励磁功率单元增益的理论计算分析
静止自并励励磁系统采用静止可控硅三相全波可控整流方式。励磁功率单元增益与移相触发控制方式以及移相触发控制电平定标方式有关。目前典型的移相触发控制方式有余弦移相、非完全余弦移相二种类型。不同的移相触发控制方式，可控硅整流器的增益特性各不相同。以下按移相触发方式分别阐述功率单元增益的理论计算方法。
2.1采用余弦移相触发的励磁功率单元增益计算
2.1.1不计及换相电抗压降情况
三相全波可控整流器输出的空载直流电压为：
Ud0＝1.35U2cosα (9)
式中，U2为励磁变压器二次侧线电压。当可控硅整流采用余弦移相控制时，移相控制的同步电压幅值USYN为整流桥交流侧电压标幺值，即：
而移相触发角α由下式求得:
其中VCB为控制电平的基准值。将式(11)代入式(9)，可得:
由式(12)可见，在整流器直流输出达到限幅值之前，输出Ud0正比于VC，与交流侧电压U2大小无关。选取发电机励磁电压基准值为EFDB，则励磁功率单元增益为：
由(13)可见，KA2取决于VC定标方式和励磁电压基准值EFDB，与励磁变二次侧电压无关。
2.1.2计及换相电抗压降情况
当计及可控整流器换相电抗压降，励磁电压和励磁电流有名值之间关系为：
Vfd＝Ud0-Ifd@t℃XC (14)
式中XC为可控整流器换相电抗有名值，稳态下励磁电压和励磁电流之间满足如下关系：
Vfd＝Ifd@t℃Rfd@t℃ (15)
将(15)代入(14)消去Ifd@t℃，考虑式(11)的定标方式，则有:
选取励磁绕组电阻基准值RFDB＝Rfd@t℃，因此有：
则得计及换相电抗压降后可控硅整流器增益系数为:
式中，KC＝XC/Rfd@t℃为换相压降系数(即换相电抗的标幺值，一般KC≈1)。对比(13)和(18)可知，计及换相压降系数后的增益略有减小，是不计及换相电抗压降增益的1/(1+KC)。
2.2采用非完全余弦移相触发的励磁功率单元增益
非完全余弦移相触发控制方式下，与余弦移相方式不同，其移相触发角采用下式确定：
式中UC0为按照控制电平基准值VC＝VCmax定标常数，则可控硅整流器空载输出电压为：
可控硅整流器单元增益为：
可见，可控硅整流器单元增益与励磁变二次侧交流电压大小有关。
如果整定控制电平达最大值(1.0p.u.)时，移相触发角α＝0°，则UC0＝1.0。可控硅整流器放大倍数为：
比较式(13)和(22)可见，非完全余弦移相控制方式下，只有当可控硅整流器交流电源电压恒定时，整流输出直流电压才与控制电平之间成比例。
类似地，计及换相电抗压降，非完全余弦移相触发控制方式下励磁功率单元增益为:
2.3励磁绕组电阻基准值和换弧压降对功率单元增益的影响分析
以余弦移相触发方式为例，假定选取的励磁绕组电阻基准值RFDB，励磁电压基准值VFDB＝RFDBIFDB，采用式(13)得到采用可控硅空载输出计算的功率单元增益为：
先假设励磁绕组电阻实际值应该为Rfd@t℃，且假设RFDB≠Rfd@t℃，励磁电压基准值为EFDB＝Rfd@t℃IFDB。在考虑换弧压降和励磁电阻基准值为Rfd@t℃情况下，由式(17)可得功率单元的准确增益应为：
由(25)并考虑式(24)可得：
式中，RFD表示当前励磁绕组电阻Rfd@t℃的标幺值，可由下式求取：
式(26)定量描述了换弧压降以及电阻基准值差异对增益计算结果的影响，误差大小取决于(KC+RFD)的数值大小。
3励磁系统开环增益系数的实测方法
从式(18)、(23)可见，理论计算方法涉及到励磁电压基准值(由励磁绕组电阻基准值确定)选择准确性问题，在实际励磁系统建模时通常对开环综合增益通过实测进行确定和校核，包括励磁控制单元增益KR校核和功率单元增益KA的确定。KR校核可采用频域响应测试或时域阶跃测试，已有的实测工作经验表明控制器增益整定值KR与校核值具有很好的一致性，可直接采用；
功率单元增益KA需要采用动态试验实测予以确定，试验主要步骤包括：励磁系统直接连接发电机励磁绕组，励磁功率整流单元由恒定电源供电，发电机维持额定转速运行，调整励磁系统电压给定值VREF，使得发电机电压逐渐升高至额定电压左右为止，期间分别记录多组发电机机端电压Vt、励磁电压Vfd、励磁电流Ifd数据。然后基于实测数据计算，常规方法和本发明方法阐述如下。
3.1常规方法
第一步：按照空载试验数据作出空载特性曲线，从空载特性曲线求取气隙线，从气隙线得到额定电压下的励磁电流作为励磁电流基准；一般用额定励磁电压除以额定励磁电流得到的励磁电阻作为励磁电阻基准，或者采用厂家给出某温度下的励磁绕组电阻作为基准；最后励磁电压基准值等于励磁电流基准值乘以励磁电阻基准值。
第二步：对每组记录数据中的励磁电压和励磁电流计算整流器空载输出电压：
第二步：采用选取的励磁系统基准值，对测量得到的机端电压偏差、励磁电压、励磁电流以及整流器输出开路电压数据进行标幺化处理；
第三步：根据每一步调整电压给定值产生的电压偏差变化量和整流器输出开路电压变化量，计算励磁系统实测静态综合放大倍数(可采用平均值或取线性化后的斜率)：
第四步：根据励磁系统实测静态综合放大倍数K'、励磁控制器实测整定增益K'R、励磁控制器运行整定增益KR，可计算励磁功率单元的增益KA。当采用完全余弦移相触发方式时，KA由式(1)得到；当采用非完全余弦移相触发方式时，KA由式(2)计算得到；则正常运行时励磁系统静态开环综合增益K由式(3)计算得到。
3.2本发明方法
严格来说，励磁系统开环综合增益应根据式(6)按照整流器负载输出励磁电压来计算，而上述常规法(见式(28)、(29))采用了整流器空载输出电压计算开环综合增益的方法，忽略了整流换相压降的影响。另一方面，从1.3节分析可知，按照IEEE推荐的励磁系统基准值定义，在稳态情况下，模型输出的励磁电压与励磁电流的标幺值应该是相等的，即EFD≡IFD，即由(6)式可得：
按照式(30a)、(30b)计算开环综合增益时，也考虑了整流换相压降的影响。当励磁电流基准值可以由空载特性曲线准确获取，则励磁电流标幺值就是准确的。另外，式(30a)、(30b)也并没有涉及到励磁绕组电阻基准值选择问题，可无需考虑励磁绕组电阻基准值与实际值不一致的问题，巧妙地解决了励磁绕组电阻测量难的问题。
综上所述，本发明方法实施步骤如下：
第一步：确定励磁电流基准值。进行发电机空载特性试验，按照空载特性试验记录的各点机端电压、励磁电流、励磁电压数据；为进一步保证空载特性气隙线的准确性，对三相机端电压进行平均得到平均机端电压，绘制平均机端电压-励磁电流的上升和下降二条曲线，对二条曲线进再行机端电压值的平均化处理，得到一条空载曲线，然后平移到过零点，从零点作切线得到空载气隙线；然后在气隙线上计算对应于额定机端电压的励磁电流，得到励磁电流的基准值。
第二步：采用第一步得到励磁电流基准值以及额定机端电压，对对动态试验测量得到的励磁电流、机端电压偏差进行标幺化处理，得到对应的标幺值IFD、ΔV＝VREF-Vt；
第三步：采用式(30a)、(30b)计算静态开环综合增益。针对多个数据一般作出ΔV＝VREF-Vt与励磁电流IFD的关系曲线，采用最小二乘拟合法进行直线拟合，拟合直线的斜率即为实测条件下的静态开环综合增益K'。
第三步：根据励磁系统实测静态综合放大倍数K'、励磁控制器实测整定增益K'R、励磁控制器运行整定增益KR，可计算励磁功率单元的增益KA。当采用完全余弦移相触发方式时，KA由式(1)得到；当采用非完全余弦移相触发方式时，KA由式(2)计算得到；则正常运行时励磁系统静态开环综合增益K由式(3)计算得到。
实施本发明的上述技术方案的有益效果包括：既能考虑换相压降的影响，又避免励磁绕组基准值选取困难问题，提高静态开环综合增益计算结果的准确性。
第二部分 实例验证
4实例分析和验证
4.1实例1
洪家渡3#机组采用静止自并励励磁系统，使用广州擎天实业有限公司EXC-9200微机励磁控制器，它是双通道励磁控制器，控制方式是串联型PID+PSS控制方式，采用完全余弦移相触发方式，相关结构和参数如图3所示。
如图4所示，首先按照空载特性曲线确定有效准确的气隙线，然后取定子额定电压对应励磁电流作为励磁电流的基准值：IFDB＝670.2A。
选取发电机铭牌额定励磁电压与额定励磁电流之比为发电机励磁绕组电阻的基准值，即RFDB＝UfN/IfN＝300/1408＝0.213Ω，则励磁电压基准为UFDB＝RFDBIFDB＝142.75V。相关数据见表1。
表1洪家渡3#机组励磁系统综合增益实测数据
对上表数据进行曲线拟合如图5所示，拟合直线的斜率即为静态综合增益K'。另外实测时K'R＝22，实际运行时KR＝36，按照式(1)计算得到KA，计算结果如表2。其中理论值采用式(18)考虑换弧压降公式求得，其中的换弧压降系数KC＝0.0619。
表2洪家渡3#机组励磁系统实测增益参数结果
从上表可看出，三种方法求出的综合增益在数值上有较大差异。为比较准确性，采用实测空载机端电压5％实测曲线进行核对，图6为BPA仿真曲线与实测曲线对比。
仿真曲线与实测曲线的指标对比见表3，结合图6可见，本发明方法、理论方法的指标均全部达标，而本发明方法在上升时间、峰值时间的偏差为0，效果更好；而常规方法效果最差，其上升时间、峰值时间指标均超标。充分验证本发明方法更为精确和有效。
表3洪家渡3#机组实测仿真曲线指标对比
5.2实例2
昌江核电2#机组静止自并励励磁系统开环综合增益实测为例，采用的是上海ABB工程有限公司产品(型号为T6S-O/U251-S6000)，采用了串联型的PID控制，为非完全移相触发方式，其测试系统结构及方案如图7所示。
如图8所示，首先按照空载特性曲线确定有效的气隙线，然后取定子额定电压对应励磁电流作为励磁电流的基准值为IFDB＝1204.8A。
选取发电机铭牌额定励磁电压与额定励磁电流之比为励磁绕组电阻基准值，即RFDB＝0.1026Ω，则励磁电压基准则为UFDB＝123.57V。计算增益的相关数据见表4。
表4昌江核电2#机组励磁系统综合增益实测数据
同样对上表数据进行曲线拟合如图9所示，可得静态增益K'。注意到二种方法得到的拟合直线斜率基本相等，直线几乎重合了。原因解释如下。因：
VD0＝VFD+KCIFD
即：
当时，则两种方法得到的结果就会基本相同。本例中实测数据表可取平均值：
即因此即二种方法结果恰好相同。
由于采用非完全移相触发方式，实测增益实验时，励磁变二次侧电压294.3V，控制器增益设置为60；而励磁系统正常运行时，励磁变二次侧额定电压为1000V，控制器增益整定为500，励磁功率增益需按照式(2)进行计算，计算结果如表5所示：
表5昌江核电#2机组励磁系统增益实测结果
图5中采用常规方法拟合得到直线的解析式为Y＝17.05x-0.018，采用发明方法拟合得到直线的解析式为Y＝17.09x-0.021。
采用BPA程序进行机端5％电压阶跃仿真，仿真曲线与实测曲线对比如图10所示，可见曲线吻合程度较好。
仿真曲线与实测曲线的指标对比见表6，可见偏差均满足指标要求。本实例也验证本发明方法的有效性。
表6昌江核电2#机组实测仿真曲线指标对比
根据在上述实施例的一种静态开环综合增益的计算方法，解决了常规方法存在的问题，既能考虑换相压降的影响，又避免励磁绕组基准值选取困难问题，提高静态开环综合增益计算结果的准确性。
需要说明的是，上述一种静态开环综合增益的计算方法的实施方式中，空载特性试验仅是针对某台机组的情况，实际应用中可以根据实际机组的具体情况，确定试验的最低电压和最高电压，也可以采用单相电压代替三相平均电压以完成空载特性曲线制作。
此外，上述一种静态开环综合增益的计算方法的实施方式中，励磁控制器类型仅仅是针对某具体机组实际配置情况，实际应用中可以根据实际机组配置励磁控制单元类型和方式，如串联型改为并型或者串-并联混合型，同样适用于本发明方法。
此外，上述一种静态开环综合增益的计算方法的实施方式中，励磁电流基准值除采用本发明方法获得外，如能有其他方式得到更为精确的数值，也可在本发明中直接采用。
本实施例还提供一种静态开环综合增益的计算装置，用于执行上述实施例中所描述的方法，结合图11所示，计算装置包括：
测试单元111，用于通过实测发电机空载特性确定励磁电流基准值。采用外部恒压源供电进行空载升压试验测得机端参考电压、机端电压、励磁电流。
数据处理单元112，用于对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统的静态开环综合增益。
可选的，数据处理单元具体用于确定功率环节部分增益KA。计算KR与KA之积，作为正常运行时励磁系统的静态开环综合增益K，即：K＝KR×KA。
其中，KR为正常运行条件下励磁控制单元的增益设定值。
可选的，数据处理单元112用于确定实测条件下励磁系统的静态开环综合增益K'。采用完全余弦移相触发方式时，其中，K'R为实测时励磁控制单元静态增益的设定值。
或者，数据处理单元112用于确定实测条件下励磁系统的静态开环综合增益K'。当采用非完全余弦移相触发方式时，
其中，K'R为实测时励磁控制单元静态增益的设定值，U2N、U2分别为励磁变二次侧额定励磁电压、实测时励磁变二次电压。
可选的，测试单元111，具体用于设置励磁控制单元为纯比例环节，采用外部恒电压源为励磁变提供电源，发电机空载运行，机端电压从0.4p.u.逐渐升压然后再降压，机端电压最高升压到1.0p.u.。
数据处理单元112，具体用于记录升压和降压过程中各点的参考电压、机端电压、励磁电流。根据励磁电流基准值计算实测励磁电流的标幺值IFD，用额定机端电压分别计算实测记录参考电压的标幺值VREF和机端电压的标幺值Vt，作出IFD与差值电压ΔV＝VREF-Vt的关系曲线，并采用最小二乘拟合法进行直线拟合，拟合直线的斜率即为实测条件下的静态开环综合增益K'。
以上为对本发明所提供的一种静态开环综合增益的计算方法及装置的描述，对于本领域的一般技术人员，依据本发明实施例的思想，在具体实施方式及应用范围上均会有改变之处，综上，本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种静态开环综合增益的计算方法，其特征在于，包括：
通过实测发电机空载特性确定励磁电流基准值；
采用外部恒压源供电进行空载升压试验测得机端参考电压、机端电压、励磁电流；
对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统的静态开环综合增益。
2.根据权利要求1所述的计算方法，其特征在于，所述通过实测发电机空载特性确定励磁电流基准值，包括：
发电机空载运行，机端电压从0.1p.u.逐渐升压到1.2p.u.，然后再降压到0.1p.u.，记录升压和降压过程中各点的机端电压、励磁电流；
对三相机端电压进行平均得到平均机端电压，绘制平均机端电压对应于励磁电流的上升和下降二条曲线，对二条曲线再进行机端电压值的平均化处理，得到空载曲线，然后平移到过零点，从零点作切线得到空载气隙线；
在气隙线上对应于额定机端电压的励磁电流，即为所述励磁电流基准值。
3.根据权利要求1所述的计算方法，其特征在于，所述计算得到正常运行时励磁系统的静态开环综合增益，包括：
确定励磁功率单元增益KA；
计算KR与KA之积，作为正常运行时励磁系统的静态开环综合增益K，即K＝KR×KA；
其中，KR为正常运行条件下励磁控制单元增益的设定值。
4.根据权利要求3所述的计算方法，其特征在于，所述确定功率环节部分增益KA，包括：
确定实测条件下励磁系统的静态开环综合增益K'；
采用完全余弦移相触发方式时，
其中，K'R为实测时励磁控制单元静态增益的设定值。
5.根据权利要求3所述的计算方法，其特征在于，所述确定功率环节部分增益KA，包括：
确定实测条件下励磁系统的静态开环综合增益K'；
当采用非完全余弦移相触发方式时，
其中，K'R为实测时励磁控制单元静态增益的设定值，U2N、U2分别为励磁变二次侧额定励磁电压、实测时励磁变二次电压。
6.根据权利要求4或5所述的计算方法，其特征在于，所述确定实测条件下励磁系统的静态开环综合增益K'，包括：
设置励磁控制单元为纯比例环节，发电机空载运行，机端电压从0.4p.u.逐渐升压到1.0p.u.，记录升压和降压过程中各点的参考电压、机端电压、励磁电流；
根据励磁电流基准值计算实测励磁电流的标幺值IFD，用额定机端电压分别计算实测记录参考电压的标幺值VREF和机端电压的标幺值Vt，作出IFD与差值电压ΔV＝VREF-Vt的特性曲线，并采用最小二乘拟合法进行直线拟合，拟合直线的斜率即为实测条件下的静态开环综合增益K'。
7.一种静态开环综合增益的计算装置，其特征在于，包括：
测试单元，用于通过实测发电机空载特性确定励磁电流基准值；采用外部恒压源供电进行空载升压试验测得机端参考电压、机端电压、励磁电流；
数据处理单元，用于对测得的数值进行标幺化后，采用励磁电流标幺值计算得到正常运行时励磁系统的静态开环综合增益。
8.根据权利要求7所述的计算装置，其特征在于，
所述数据处理单元具体用于确定励磁功率单元增益KA；计算KR与KA之积，作为正常运行时励磁系统的静态开环综合增益K，即K＝KR×KA；
其中，KR为正常运行条件下励磁控制单元增益的设定值。
9.根据权利要求8所述的计算装置，其特征在于，
所述数据处理单元用于确定实测条件下励磁系统的静态开环综合增益K'；
采用完全余弦移相触发方式时，
其中，K'R为实测时励磁控制单元静态增益的设定值；
或者，所述数据处理单元用于确定实测条件下励磁系统的静态开环综合增益K'；当采用非完全余弦移相触发方式时，
其中，K'R为实测时励磁控制单元静态增益的设定值，U2N为励磁变二次侧额定励磁电压，U2为实测时励磁变二次电压。
10.根据权利要求9所述的计算装置，其特征在于，
所述测试单元，具体用于设置励磁控制单元为纯比例环节，采用外部恒电压源为励磁变提供电源，发电机空载运行，参考电压从0.4p.u.逐渐升压到1.0p.u.所述数据处理单元，具体用于记录升压过程中各点的参考电压、机端电压、励磁电流；根据励磁电流基准值计算实测励磁电流的标幺值IFD，用额定机端电压分别计算实测记录参考电压的标幺值VREF和机端电压的标幺值Vt，作出IFD与差值电压ΔV＝VREF-Vt的特性曲线，并采用最小二乘拟合法进行直线拟合，拟合直线的斜率即为实测条件下的静态开环综合增益K'。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710014788.5A CN106684866A (zh) | 2017-01-09 | 2017-01-09 | 一种静态开环综合增益的计算方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710014788.5A CN106684866A (zh) | 2017-01-09 | 2017-01-09 | 一种静态开环综合增益的计算方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106684866A true CN106684866A (zh) | 2017-05-17 |
Family
ID=58850494
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710014788.5A Pending CN106684866A (zh) | 2017-01-09 | 2017-01-09 | 一种静态开环综合增益的计算方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106684866A (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107966615A (zh) * | 2017-12-05 | 2018-04-27 | 华北电力大学 | 一种集成运放开环幅频特性测试系统 |
| CN109633443A (zh) * | 2018-12-28 | 2019-04-16 | 华中科技大学 | 一种适合于自励磁大型同步电机空载特性测量方法 |
| CN110138295A (zh) * | 2019-04-29 | 2019-08-16 | 中国电力科学研究院有限公司 | 一种测量自并励励磁系统电压控制环静态放大倍数的方法及系统 |
| CN112383252A (zh) * | 2020-10-30 | 2021-02-19 | 华北电力科学研究院有限责任公司 | 双馈发电机组励磁控制系统标幺方法及装置 |
| CN112769366A (zh) * | 2020-12-31 | 2021-05-07 | 徐州中矿大传动与自动化有限公司 | 一种电励磁同步电机励磁变流器控制方法及装置、系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103823183A (zh) * | 2014-02-28 | 2014-05-28 | 国家电网公司 | 一种同步发电机定子开路瞬变时间常数的测量方法 |
| CN104730936A (zh) * | 2015-03-05 | 2015-06-24 | 国家电网公司 | 核电机组39相自励无刷励磁系统机电暂态模型实测建模方法 |
-
2017
- 2017-01-09 CN CN201710014788.5A patent/CN106684866A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103823183A (zh) * | 2014-02-28 | 2014-05-28 | 国家电网公司 | 一种同步发电机定子开路瞬变时间常数的测量方法 |
| CN104730936A (zh) * | 2015-03-05 | 2015-06-24 | 国家电网公司 | 核电机组39相自励无刷励磁系统机电暂态模型实测建模方法 |
Non-Patent Citations (2)
| Title |
|---|
| 单文培主编: "《电气设备试验及故障处理实例》", 28 February 2006, 中国水利水电出版社 * |
| 汪盛波: "发电机励磁系统精确建模与参数辨识研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107966615A (zh) * | 2017-12-05 | 2018-04-27 | 华北电力大学 | 一种集成运放开环幅频特性测试系统 |
| CN109633443A (zh) * | 2018-12-28 | 2019-04-16 | 华中科技大学 | 一种适合于自励磁大型同步电机空载特性测量方法 |
| CN110138295A (zh) * | 2019-04-29 | 2019-08-16 | 中国电力科学研究院有限公司 | 一种测量自并励励磁系统电压控制环静态放大倍数的方法及系统 |
| CN112383252A (zh) * | 2020-10-30 | 2021-02-19 | 华北电力科学研究院有限责任公司 | 双馈发电机组励磁控制系统标幺方法及装置 |
| CN112383252B (zh) * | 2020-10-30 | 2022-05-06 | 华北电力科学研究院有限责任公司 | 双馈发电机组励磁控制系统标幺方法及装置 |
| CN112769366A (zh) * | 2020-12-31 | 2021-05-07 | 徐州中矿大传动与自动化有限公司 | 一种电励磁同步电机励磁变流器控制方法及装置、系统 |
| CN112769366B (zh) * | 2020-12-31 | 2024-01-09 | 江苏国传电气有限公司 | 一种电励磁同步电机励磁变流器控制方法及装置、系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106684866A (zh) | 一种静态开环综合增益的计算方法及装置 | |
| CN104221241B (zh) | 用于控制变频器的方法和装置 | |
| CN106020168B (zh) | 双馈风机次同步谐振硬件在环测试系统及方法 | |
| Lenwari et al. | The use of genetic algorithms for the design of resonant compensators for active filters | |
| CN106680714A (zh) | 一种励磁系统正常运行时功率单元增益的计算方法及装置 | |
| CN101364735B (zh) | 燃气轮机发电机组调速系统负荷频率调节测试的试验方法 | |
| CN101719753B (zh) | 抽水蓄能电站静止变频启动电机转子初始位置检测方法 | |
| Cao et al. | Research on PID parameters optimization of synchronous generator excitation control system | |
| Buja et al. | Steady-state performance degradation of a DTC IM drive under parameter and transduction errors | |
| CN109066793A (zh) | 一种微网逆变器的虚拟柴油发电机组控制方法 | |
| CN111130135A (zh) | 一种适用于高比例新能源接入下的电力系统惯量计算方法 | |
| US9910076B2 (en) | Method and system for analysis of energy savings in an electrical network | |
| CN108988387A (zh) | 基于转子电流移相平均的双馈风电机组次同步谐振抑制方法 | |
| CN109936148B (zh) | 一种发电机组无功电流补偿率的测量方法及在线监测系统 | |
| Bossa et al. | A field test to determine PSS effectiveness at multigenerator power plants | |
| CN110829465B (zh) | 计及多调速器死区的电力系统超低频振荡风险评估方法 | |
| CN108363007A (zh) | 一种光伏虚拟同步发电机性能测试装置及方法 | |
| CN110474315B (zh) | 一种用于柔性直流小干扰稳定性判定的方法及系统 | |
| CN112968447B (zh) | 一种衡量电压测量时间常数对机组动态稳定影响的方法 | |
| CN113078646A (zh) | 一种计及逆变电源控制器饱和特性的故障暂态电流解析方法 | |
| CN114997505A (zh) | 一种经变流器并网电源的稳定性分析方法及系统 | |
| CN112736935A (zh) | 一种利用pss补偿角在线校核电力系统稳定器模型参数的方法 | |
| Skamyin et al. | Influence of nonlinear load on the measurement of harmonic impedance of the power supply system | |
| CN110879357A (zh) | 一种基于同步相量测量的强迫功率振荡扰动源定位方法 | |
| CN110514892B (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170517 |
|
| RJ01 | Rejection of invention patent application after publication |