CN109120196A - 一种水电机组调速***pid参数整定方法及其装置 - Google Patents
一种水电机组调速***pid参数整定方法及其装置 Download PDFInfo
- Publication number
- CN109120196A CN109120196A CN201810956546.2A CN201810956546A CN109120196A CN 109120196 A CN109120196 A CN 109120196A CN 201810956546 A CN201810956546 A CN 201810956546A CN 109120196 A CN109120196 A CN 109120196A
- Authority
- CN
- China
- Prior art keywords
- amplification factor
- pid parameter
- proportional component
- hydropower unit
- oscillation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000001105 regulatory effect Effects 0.000 title abstract 3
- 230000003321 amplification Effects 0.000 claims abstract description 196
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 196
- 238000004088 simulation Methods 0.000 claims abstract description 104
- 230000010355 oscillation Effects 0.000 claims abstract description 94
- 230000004069 differentiation Effects 0.000 claims abstract description 56
- 230000008859 change Effects 0.000 claims description 33
- 238000013016 damping Methods 0.000 claims description 27
- 238000005457 optimization Methods 0.000 claims description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 230000010354 integration Effects 0.000 abstract 1
- 238000012795 verification Methods 0.000 abstract 1
- 230000001052 transient effect Effects 0.000 description 13
- 244000098338 Triticum aestivum Species 0.000 description 10
- 230000006872 improvement Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005404 monopole Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种水电机组调速***PID参数整定方法及其装置,其中,水电机组调速***PID参数整定方法包括步骤:采用预设的仿真模型对水电机组进行仿真,确定达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;根据达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;根据所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用预设的仿真模型进行仿真,确定达到临界稳定状态的积分环节放大倍数,作为PID参数中的积分环节放大倍数;根据PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取预设的仿真模型进行仿真验证,完成所述PID参数的整定。
Description
技术领域
本发明涉及电力***机网协调技术领域,尤其涉及一种水电机组调速***PID参数整定方法及其装置。
背景技术
水电机组调速***参数配置对电力***频率的动态特性有重要影响,不同规模的电力***,需要配置不同的水电机组调速***参数。大规模电网在实施异步联网后,原来属于大规模电网的某个地区电网将与大规模电网分离出来异步运行。异步运行后,如果地区电网与大规模电网相比,其装机和负荷等规模大幅减小,地区电网中水电机组原有的调速***参数将无法适应异步运行后的地区电网,容易导致地区电网在异步联网运行后出现长时间、大幅度、超低频率的振荡现象。
目前对大规模电网中水电机组调速***参数配置方法,以及单个水电厂带局部负荷的孤立电网中水电机组调速***参数配置方法有相应研究和成果,但对中等规模电网中如何合理配置水电机组调速***参数,优化电网频率的动态特性,避免出现超低频振荡问题还缺少相应的研究。
大规模电网中水电机组调速***配置的PID参数,容易导致机组在中等规模电网运行处响应速度过快而出现超低频率振荡阻尼不足的问题;而单个水电厂带局部负荷的孤立电网中水电机组调速***配置的PID参数,容易导致机组在中等规模电网运行处响应速度过慢而出现一次调频响应速度过慢的问题。
发明内容
本发明实施例提供一种水电机组调速***PID参数整定方法及其装置,能有效解决现有技术对中等规模电网中如何合理配置水电机组调速***参数的问题,能有效优化电网频率的动态特性,避免出现超低频振荡的问题。
本发明一实施例提供一种水电机组调速***PID参数整定方法方法,包括:
采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;
根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;
根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数;
根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
作为上述方案的改进,所述仿真模型的构建包括步骤:
将所述PID参数中的KP=0,KD=0,KI=0输入所述仿真模型进行甩负荷仿真;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数;
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于所述水电机组额定功率的2%。
作为上述方案的改进,采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期具体为:
通过所述仿真程序进行甩负荷仿真,逐步提高所述比例环节放大倍数的取值,以得到第一频率变化曲线;
获取当所述第一频率变化曲线出现正弦等幅振荡时对应的比例环节放大倍数值和正弦振荡周期,作为所述水电机组达到临界稳定状态的比例环节放大倍数和对应的正弦振荡周期。
作为上述方案的改进,所述根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数具体为:
根据公式(1)计算所述PID参数中的比例环节放大倍数
KP1=0.6×KPMAX (1)
其中,KP1为所述PID参数中的比例环节放大倍数,KPMAX为所述水电机组达到临界稳定状态时的比例环节放大倍数;
根据公式(2)计算所述PID参数中的微分环节放大倍数
KD1=min{0.075×KPMAX×T,KD0} (2)
其中,KD1为所述PID参数中的微分环节放大倍数,T为所述水电机组达到临界稳定状态时对应的振荡周期,KD0为所述水电机组正常运行时的典型微分环节放大倍数值;min{0.075×KPMAX×T,KD0}表示取两个参数中较小的数值。
作为上述方案的改进,所述根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数,之前包括:
将所述仿真模型的PID参数设置为KP=KP1,KD=KD1,KI=0;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数。
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于所述水电机组额定功率的2%。
作为上述方案的改进,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数具体为:
将所述PID参数中的比例环节放大倍数和微分环节放大倍数输入所述仿真模型进行甩负荷仿真,逐步提高所述积分环节放大倍数的取值,以得到第二频率变化曲线;
获取当所述第二频率变化曲线出现正弦等幅振荡时对应的积分环节放大倍数值,作为所述水电机组达到临界稳定状态的积分环节放大倍数。
作为上述方案的改进,所述根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定,具体为:
将所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数输入所述仿真模型进行故障仿真,以得到第三频率变化曲线;
验证所述PID参数是否优化超低频振荡阻尼,在通过验证后完成所述PID参数的整定。
作为上述方案的改进,通过以下方式验证所述PID参数是否优化超低频振荡阻尼:
对比第一频率变化曲线与所述第三频率变化曲线的差异;其中,所述第一频率变化曲线由PID参数初始为KP=0,KD=0,KI=0输入所述仿真模型进行仿真得到;
当检测到所述第三频率变化曲线的超低频振荡阻尼高于所述第一频率变化曲线的超低频振荡阻尼时,则确定所述PID参数优化超低频振荡阻尼。
作为上述方案的改进,所述方法还包括:
当检测到所述第三频率变化曲线的超低频振荡阻尼不高于所述第一频率变化曲线的超低频振荡阻尼时,则重新采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期。
本发明另一实施例对应提供了一种水电机组调速***PID参数整定装置,其特征在于,包括:
第一模块,用于采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;
第二模块,用于根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;
第三模块,用于根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数;
第四模块,用于根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
与现有技术相比,本发明实施例公开的一种水电机组调速***PID参数整定方法,其通过采用预设的仿真模型进行仿真,得到水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,根据上述得到的值,计算PID参数中的比例环节放大倍数和微分环节放大倍数,然后根据上述得到的PID参数值,采用预设的仿真模型进行仿真,确定PID参数中的积分环节放大倍数,通过仿真模型进行仿真并验证上述PID参数值,完成所述PID参数的整定,既能有效解决中等规模电网中如何合理配置水电机组调速***参数的问题,又能优化电网频率的动态特性,避免出现超低频振荡,以实现一次调频响应速度和超低频振荡阻尼间的协调优化。
附图说明
图1是本发明一实施例提供的一种水电机组调速***PID参数整定方法的流程示意图;
图2是本发明一实施例提供电网中主要的水电机组调速***PID参数整定方法的过程示意图;
图3是本发明一实施例提供的#1糯扎渡水电机组初始参数设置下达到临界稳定状态时的频率变化曲线图;
图4是本发明一实施例提供的#1糯扎渡水电机组在设置KP和KD值下达到临界稳定状态时的频率变化曲线图;
图5是本发明一实施例提供的PID参数优化前后的频率变化曲线图;
图6是本发明一实施例提供的一种水电机组调速***PID参数整定装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1,是本发明一实施例提供的一种水电机组调速***PID参数整定方法的流程示意图,包括步骤:
S11、采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期。
S12、根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数。
S13、根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数。
S14、根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
在本实施例中,利用对水电机组进行时域仿真,并计算确定PID参数中的比例环节放大倍数、微分环节放大倍数和积分环节放大倍数,然后通过对得到的PID参数进行仿真,并验证其优化电网超低频振荡的阻尼特性方面的效果,完成所述PID参数的整定,既能有效解决中等规模电网中如何合理配置水电机组调速***参数的问题,又能优化电网频率的动态特性,避免出现超低频振荡。
在另一优选实施例中,所述方法还包括仿真模型的构建,具体为:
将所述PID参数中的KP=0,KD=0,KI=0输入所述仿真模型进行甩负荷仿真;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数;
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于水电机组额定功率的2%。
优选地,该仿真模型的构建在机电暂态仿真程序中运行,构造一个单台水电机组带负荷运行的仿真模型进行仿真。
具体地,步骤S11中采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期具体为:
通过所述仿真程序进行甩负荷仿真,逐步提高所述比例环节放大倍数的取值,以得到第一频率变化曲线;
获取当所述第一频率变化曲线出现正弦等幅振荡时对应的比例环节放大倍数值和正弦振荡周期,作为所述水电机组达到临界稳定状态的比例环节放大倍数和对应的正弦振荡周期。
优选地,逐步提高所述比例环节放大倍数取值具体为取值0.5n,若频率变化曲线趋向于正弦振荡,减小取值范围,直到出现正弦等幅振荡;若频率变化曲线未出现明显变化,继续按0.5n取值递增;其中,n=0,1,2……。
例如,比例环节放大倍数首次取值0,二次取值0.5,若频率变化曲线趋向于正弦振荡,三次取值为0.51,直到出现正弦等幅振荡。
即利用机电暂态仿真程序进行甩负荷仿真,逐步提高KP的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频率出现正弦等幅振荡,记录此时的KP取值KPMAX和正弦振荡的周期T。
具体地,步骤S12中具体为
根据公式(1)计算所述PID参数中的比例环节放大倍数
KP1=0.6×KPMAX (1)
其中,KP1为所述PID参数中的比例环节放大倍数,KPMAX为所述水电机组达到临界稳定状态时的比例环节放大倍数;
根据公式(2)计算所述PID参数中的微分环节放大倍数
KD1=min{0.075×KPMAX×T,KD0} (2)
其中,KD1为所述PID参数中的微分环节放大倍数,T为所述水电机组达到临界稳定状态时对应的振荡周期,KD0为所述水电机组正常运行时的典型微分环节放大倍数值;min{0.075×KPMAX×T,KD0}表示取两个参数中较小的数值。
其中,步骤S13之前还包括对仿真模型加入KP1和KD1的值,具体为:
将所述仿真模型的PID参数设置为KP=KP1,KD=KD1,KI=0;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数。
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于水电机组额定功率的2%。
可以理解的是,采用与步骤S11中同样的单台水电机组带负荷运行的仿真模型,并将调速***PID参数设置为KP=KP1,KD=KD1,KI=0,并设置与步骤S11中同样的甩负荷故障。
具体地,步骤S13中根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数具体为:
将所述PID参数中的比例环节放大倍数和微分环节放大倍数输入所述仿真模型进行甩负荷仿真,逐步提高所述积分环节放大倍数的取值,以得到第二频率变化曲线;
获取当所述第二频率变化曲线出现正弦等幅振荡时对应的积分环节放大倍数值,作为所述水电机组达到临界稳定状态的积分环节放大倍数。
即利用机电暂态仿真程序进行甩负荷仿真,逐步提高KI的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频率出现正弦等幅振荡,记录此时的KI取值KI1。
具体地,步骤S14中根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定具体为:
将所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数输入所述仿真模型进行故障仿真,以得到第三频率变化曲线;
验证所述PID参数是否优化超低频振荡阻尼,在通过验证后完成所述PID参数的整定。
其中,所述故障包括直流闭锁、电厂机组跳闸或甩负荷;
优选地,通过以下方式验证所述PID参数是否优化超低频振荡阻尼:
对比第一频率变化曲线与所述第三频率变化曲线的差异;其中,所述第一频率变化曲线由PID参数初始为KP=0,KD=0,KI=0输入所述仿真模型进行仿真得到;
当检测到所述第三频率变化曲线的超低频振荡阻尼高于所述第一频率变化曲线的超低频振荡阻尼时,则确定所述PID参数优化超低频振荡阻尼。
其中,所述整定方法还包括:
当检测到所述第三频率变化曲线的超低频振荡阻尼不高于所述第一频率变化曲线的超低频振荡阻尼时,则重新采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期。
参见图2,是本实施例提供的电网中主要的水电机组调速***PID参数整定方法的过程示意图,将上述用于水电机组的调速***PID参数整定方法的流程应用于电网中主要的水电机组的PID参数整定过程。
优选地,主要的水电机组是指电网中单台机组装机容量较大的水电机组以及电厂总装机容量较大的水电机组。
具体地,对电网中主要的水电机组进行调速***PID整定。首先利用所述方法对这些水电机组逐一得到各机组的调速***PID参数;然后将得到的参数设置到机组中,并在电网中设置直流闭锁、电厂机组跳闸或甩负荷等故障;最后利用机电暂态仿真程序进行故障仿真,记录电网中重要变电站母线节点的频率变化曲线,通过对比机组调速***PID参数优化前后频率变化的差异,分析优化后的PID参数在优化电网超低频振荡的阻尼特性方面的效果,完成PID参数的整定。
在另一优选实施例中,在上述实施例的基础上,将上述水电机组调速***PID参数整定方法应用于实际电网中。
为优化云南电网运行后送出能力中电网超低频振荡阻尼特性,采用本发明方法对机组PSS参数进行整定,来验证本方法的有效性。
优选地,本实施例选取糯扎渡电厂的#1水电机组。
具体地,所述方法的实施过程如下:
参见上述实施例的步骤S11,在机电暂态仿真程序中,构造一个#1糯扎渡水电机组带负荷运行的仿真模型;其中负荷大小等于#1糯扎渡水电机组的额定功率650MW,调速***PID初始参数设置为KP=0,KD=0,KI=0;
设置一个甩负荷故障,负荷量等于10MW;
利用机电暂态仿真程序进行甩负荷仿真,逐步提高KP的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频出现正弦等幅振荡,记录此时的KP的取值KPMAX=3.5,正弦振荡的周期T=21.3秒。
其中,#1糯扎渡水电机组初始参数设置下达到临界稳定状态时的频率变化曲线图如图3所示,通过该频率变化曲线,得到KPMAX和对应的振荡周期T的取值。
优选地,本实施例中采用中国电科院开发的PSD-BPA机电暂态仿真程序。
参见上述实施例的步骤S12,根据上述步骤得到的KPMAX和T利用以下公式,得到PID参数中的比例环节放大倍数KP1,微分环节放大倍数KD1:
KP1=0.6×KPMAX=0.6×3.5=2.1
KD1=min{0.075×KPMAX×T,KD0}
=min{0.075×3.5×21.3,6}
=min{5.6,6}
=5.6
其中,KD0=3是机组大电网运行整定的微分环节放大倍数值。
参见上述实施例的步骤S13,设置#1糯扎渡水电机组调速***PID参数为:KP=2.1,KD=5.6,KI=0;
设置与步骤S11中同样的甩负荷故障;
利用机电暂态仿真程序进行甩负荷仿真,逐步提高KI的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频率出现正弦等幅振荡,记录此时的KI的取值KI1=0.45。
其中,#1糯扎渡水电机组在设置KP和KD值下达到临界稳定状态时的频率变化曲线如图4所示,通过该频率变化曲线,得到KI的取值。
参见上述实施例的步骤S14,对云南电网内主要水电机组利用步骤S11-S13的方法对逐一进行调速***PID参数优化。
设置楚穗直流单极闭锁作为故障;
利用机电暂态仿真程序对云南电网进行楚穗直流单极闭锁故障仿真,记录电网中重要变电站母线节点的频率变化曲线,对比机组调速***PID参数优化前后频率变化的差异,可以看出,优化后的PID参数明显提高了电网超低频振荡的阻尼。
例如,对其中的糯扎渡、小湾和溪洛渡电厂水电机组进行了调速***PID参数整定,主要整定结果如下:
糯扎渡:KP1=2.1,KI1=0.45,KD1=5.6
小湾:KP1=2.4,KI1=0.67,KD1=5.05
溪洛渡:KP1=3.0,KI1=1.10,KD1=4.1
参见图5,本实施例提供的PID参数优化前后的频率变化曲线图。
参见图6,为本发明一实施例提供的一种水电机组调速***PID参数整定装置的结构示意图,包括:
第一模块1,用于采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;
第二模块2,用于根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;
第三模块3,用于根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数;
第四模块4,用于根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
在一种可选的实施例中,第一模块1包括:
仿真模型初始设置单元,用于将所述PID参数中的KP=0,KD=0,KI=0输入所述仿真模型进行甩负荷仿真;
其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数;
甩负荷故障设置单元,用于以预设的负荷量设置甩负荷故障;
其中,所述负荷量不大于水电机组额定功率的2%;
第一频率变化曲线仿真单元,用于通过所述仿真程序进行甩负荷仿真,逐步提高所述比例环节放大倍数的取值,以得到第一频率变化曲线;
临界稳定状态仿真结果获取单元,用于获取当所述第一频率变化曲线出现正弦等幅振荡时对应的比例环节放大倍数值和正弦振荡周期,作为所述水电机组达到临界稳定状态的比例环节放大倍数和对应的正弦振荡周期。
优选地,该步骤在机电暂态仿真程序中运行,构造一个单台水电机组带负荷运行的仿真模型进行仿真。
即利用机电暂态仿真程序进行甩负荷仿真,逐步提高KP的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频率出现正弦等幅振荡,记录此时的KP取值KPMAX和正弦振荡的周期T。
在一种可选的实施例中,第二模块2包括:
比例环节放大倍数计算单元,用于根据公式KP1=0.6×KPMAX,计算PID参数中的比例环节放大倍数;
微分环节放大倍数计算单元,用于根据公式KD1=min{0.075×KPMAX×T,KD0},计算PID参数中的微分环节放大倍数;
其中,KP1为所述PID参数中的比例环节放大倍数,KD1为所述PID参数中的微分环节放大倍数;KPMAX为所述水电机组达到临界稳定状态时的比例环节放大倍数,T为所述水电机组达到临界稳定状态时对应的振荡周期,KD0为所述水电机组正常运行时的典型微分环节放大倍数值;min{0.075×KPMAX×T,KD0}表示取两个参数中较小的数值。
在一种可选的实施例中,第三模块3包括:
仿真模型参数设置单元,用于将所述仿真模型的PID参数设置为KP=KP1,KD=KD1,KI=0;
其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数;
甩负荷故障设置单元,用于以预设的负荷量设置甩负荷故障;
其中,所述负荷量不大于水电机组额定功率的2%;
可以理解的是,采用与第一模块1中同样的单台水电机组带负荷运行的仿真模型,并将调速***PID参数设置为KP=KP1,KD=KD1,KI=0,并设置与第一模块1中同样的甩负荷故障。
第二频率变化曲线仿真单元,用于将所述PID参数中的比例环节放大倍数和微分环节放大倍数输入所述仿真模型进行甩负荷仿真,逐步提高所述积分环节放大倍数的取值,以得到第二频率变化曲线;
积分环节放大倍数获取单元,用于获取当所述第二频率变化曲线出现正弦等幅振荡时对应的积分环节放大倍数值,作为所述水电机组达到临界稳定状态的积分环节放大倍数。
即利用机电暂态仿真程序进行甩负荷仿真,逐步提高KI的取值,观察每次甩负荷仿真时水电机组频率变化曲线,直到频率出现正弦等幅振荡,记录此时的KI取值KI1。
在一种可选的实施例中,第四模块4包括:
第三频率变化曲线仿真单元,用于将所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数输入所述仿真模型进行故障仿真,以得到第三频率变化曲线;
其中,所述故障包括直流闭锁、电厂机组跳闸或甩负荷;
PID参数验证模块,用于验证所述PID参数是否优化超低频振荡阻尼,在通过验证后完成所述PID参数的整定。
优选地,通过以下方式验证所述PID参数是否优化超低频振荡阻尼:
对比第一频率变化曲线与所述第三频率变化曲线的差异;当检测到所述第三频率变化曲线的超低频振荡阻尼高于所述第一频率变化曲线的超低频振荡阻尼时,则确定所述PID参数优化超低频振荡阻尼。
其中,所述第一频率变化曲线由PID参数初始为KP=0,KD=0,KI=0输入所述仿真模型进行仿真得到;
其中,所述方法还包括:
当检测到所述第三频率变化曲线的超低频振荡阻尼不高于所述第一频率变化曲线的超低频振荡阻尼时,则重新采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期。
例如,对主要的水电机组进行调速***PID整定;其中,主要的水电机组指电网中单台机组装机容量较大的水电机组以及电厂总装机容量较大的水电机组。
具体地,对电网中主要的水电机组进行调速***PID整定,首先利用所述方法对这些水电机组逐一得到各机组的调速***PID参数;然后将得到的参数设置到机组中,并在电网中设置直流闭锁、电厂机组跳闸或甩负荷等故障;最后利用机电暂态仿真程序进行故障仿真,记录电网中重要变电站母线节点的频率变化曲线,通过对比机组调速***PID参数修改前后频率变化的差异,分析修改后的PID参数在优化电网超低频振荡的阻尼特性方面的效果。
需说明的是,以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。另外,本发明提供的装置实施例附图中,模块之间的连接关系表示它们之间具有通信连接,具体可以实现为一条或多条通信总线或信号线。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (10)
1.一种水电机组调速***PID参数整定方法,其特征在于,包括步骤:
采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;
根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;
根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数;
根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
2.如权利要求1所述的水电机组调速***PID参数整定方法,其特征在于,所述仿真模型的构建包括步骤:
将所述PID参数中的KP=0,KD=0,KI=0输入所述仿真模型进行甩负荷仿真;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数;
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于所述水电机组额定功率的2%。
3.如权利要求2所述的水电机组调速***PID参数整定方法,其特征在于,所述采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期具体为:
通过所述仿真程序进行甩负荷仿真,逐步提高所述比例环节放大倍数的取值,以得到第一频率变化曲线;
获取当所述第一频率变化曲线出现正弦等幅振荡时对应的比例环节放大倍数值和正弦振荡周期,作为所述水电机组达到临界稳定状态的比例环节放大倍数和对应的正弦振荡周期。
4.如权利要求1所述的水电机组调速***PID参数整定方法,其特征在于,所述根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数具体为:
根据公式(1)计算所述PID参数中的比例环节放大倍数
KP1=0.6×KPMAX (1)
其中,KP1为所述PID参数中的比例环节放大倍数,KPMAX为所述水电机组达到临界稳定状态时的比例环节放大倍数;
根据公式(2)计算所述PID参数中的微分环节放大倍数
KD1=min{0.075×KPMAX×T,KD0} (2)
其中,KD1为所述PID参数中的微分环节放大倍数,T为所述水电机组达到临界稳定状态时对应的振荡周期,KD0为所述水电机组正常运行时的典型微分环节放大倍数值;min{0.075×KPMAX×T,KD0}表示取两个参数中较小的数值。
5.如权利要求1所述的水电机组调速***PID参数整定方法,其特征在于,所述根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数,之前包括:
将所述仿真模型的PID参数设置为KP=KP1,KD=KD1,KI=0;其中,KP为所述比例环节放大倍数,KD为所述微分环节放大倍数,KI为所述积分环节放大倍数。
以预设的负荷量设置甩负荷故障;其中,所述负荷量不大于所述水电机组额定功率的2%。
6.如权利要求5所述的水电机组调速***PID参数整定方法,其特征在于,所述采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数具体为:
将所述PID参数中的比例环节放大倍数和微分环节放大倍数输入所述仿真模型进行甩负荷仿真,逐步提高所述积分环节放大倍数的取值,以得到第二频率变化曲线;
获取当所述第二频率变化曲线出现正弦等幅振荡时对应的积分环节放大倍数值,作为所述水电机组达到临界稳定状态的积分环节放大倍数。
7.如权利要求1所述的水电机组调速***PID参数整定方法,其特征在于,所述根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定,具体为:
将所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数输入所述仿真模型进行故障仿真,以得到第三频率变化曲线;
验证所述PID参数是否优化超低频振荡阻尼,在通过验证后完成所述PID参数的整定。
8.如权利要求7所述的水电机组调速***PID参数整定方法,其特征在于,通过以下方式验证所述PID参数是否优化超低频振荡阻尼:
对比第一频率变化曲线与所述第三频率变化曲线的差异;其中,所述第一频率变化曲线由PID参数初始为KP=0,KD=0,KI=0输入所述仿真模型进行仿真得到;
当检测到所述第三频率变化曲线的超低频振荡阻尼高于所述第一频率变化曲线的超低频振荡阻尼时,则确定所述PID参数优化超低频振荡阻尼。
9.如权利要求8所述的水电机组调速***PID参数整定方法,其特征在于,所述方法还包括:
当检测到所述第三频率变化曲线的超低频振荡阻尼不高于所述第一频率变化曲线的超低频振荡阻尼时,则重新采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期。
10.一种水电机组调速***PID参数整定装置,其特征在于,包括:
第一模块,用于采用预设的仿真模型对水电机组进行时域仿真,确定所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期;
第二模块,用于根据所述水电机组达到临界稳定状态时的比例环节放大倍数和对应的振荡周期,计算PID参数中的比例环节放大倍数和微分环节放大倍数;
第三模块,用于根据计算所得的PID参数中的比例环节放大倍数和微分环节放大倍数,采用所述预设的仿真模型进行时域仿真,确定所述水电机组达到临界稳定状态的积分环节放大倍数,作为所述PID参数中的积分环节放大倍数;
第四模块,用于根据所述PID参数中的比例环节放大倍数、微分环节放大倍数以及积分环节放大倍数,采取所述预设的仿真模型进行仿真验证,完成所述PID参数的整定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810956546.2A CN109120196B (zh) | 2018-08-21 | 2018-08-21 | 一种水电机组调速***pid参数整定方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810956546.2A CN109120196B (zh) | 2018-08-21 | 2018-08-21 | 一种水电机组调速***pid参数整定方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109120196A true CN109120196A (zh) | 2019-01-01 |
CN109120196B CN109120196B (zh) | 2020-05-26 |
Family
ID=64853432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810956546.2A Active CN109120196B (zh) | 2018-08-21 | 2018-08-21 | 一种水电机组调速***pid参数整定方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109120196B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112366733A (zh) * | 2020-11-30 | 2021-02-12 | 南方电网科学研究院有限责任公司 | 水电机组调速器pid参数优化方法、装置、***及介质 |
CN117155171A (zh) * | 2023-10-31 | 2023-12-01 | 苏州泰科贝尔直驱电机有限公司 | 直驱电机的启动控制参数优化方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5229699A (en) * | 1991-10-15 | 1993-07-20 | Industrial Technology Research Institute | Method and an apparatus for PID controller tuning |
CN104601082A (zh) * | 2014-12-18 | 2015-05-06 | 东南大学 | 一种直线电机控制器 |
CN105487375A (zh) * | 2015-12-31 | 2016-04-13 | 中国船舶重工集团公司第七一九研究所 | 一种离散pid控制器参数整定方法 |
CN107086589A (zh) * | 2017-06-23 | 2017-08-22 | 云南电网有限责任公司 | 一种基于临界参数抑制水电机组超低频振荡的调速器参数优化方法 |
KR20180024092A (ko) * | 2016-08-26 | 2018-03-08 | 한국로봇융합연구원 | Pid 게인 자동 튜닝 방법 |
CN107800146A (zh) * | 2017-11-16 | 2018-03-13 | 国网四川省电力公司电力科学研究院 | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 |
CN107942648A (zh) * | 2017-11-29 | 2018-04-20 | 中国飞机强度研究所 | 一种超大空间温度场pid控制器参数整定方法 |
-
2018
- 2018-08-21 CN CN201810956546.2A patent/CN109120196B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5229699A (en) * | 1991-10-15 | 1993-07-20 | Industrial Technology Research Institute | Method and an apparatus for PID controller tuning |
CN104601082A (zh) * | 2014-12-18 | 2015-05-06 | 东南大学 | 一种直线电机控制器 |
CN105487375A (zh) * | 2015-12-31 | 2016-04-13 | 中国船舶重工集团公司第七一九研究所 | 一种离散pid控制器参数整定方法 |
KR20180024092A (ko) * | 2016-08-26 | 2018-03-08 | 한국로봇융합연구원 | Pid 게인 자동 튜닝 방법 |
CN107086589A (zh) * | 2017-06-23 | 2017-08-22 | 云南电网有限责任公司 | 一种基于临界参数抑制水电机组超低频振荡的调速器参数优化方法 |
CN107800146A (zh) * | 2017-11-16 | 2018-03-13 | 国网四川省电力公司电力科学研究院 | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 |
CN107942648A (zh) * | 2017-11-29 | 2018-04-20 | 中国飞机强度研究所 | 一种超大空间温度场pid控制器参数整定方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112366733A (zh) * | 2020-11-30 | 2021-02-12 | 南方电网科学研究院有限责任公司 | 水电机组调速器pid参数优化方法、装置、***及介质 |
CN117155171A (zh) * | 2023-10-31 | 2023-12-01 | 苏州泰科贝尔直驱电机有限公司 | 直驱电机的启动控制参数优化方法及*** |
CN117155171B (zh) * | 2023-10-31 | 2024-01-30 | 苏州泰科贝尔直驱电机有限公司 | 直驱电机的启动控制参数优化方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109120196B (zh) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shoults et al. | Multi-area adaptive LFC developed for a comprehensive AGC simulator | |
Zhang et al. | Multi-objective virtual inertia control of renewable power generator for transient stability improvement in interconnected power system | |
CN107086589B (zh) | 一种基于临界参数抑制水电机组超低频振荡的调速器参数优化方法 | |
CN105703382B (zh) | 一种风电场次同步振荡抑制方法和*** | |
CN109301815A (zh) | 水电机组调速***附加阻尼控制器参数整定方法与装置 | |
CN102624013B (zh) | 基于相位补偿原理的储能阻尼控制器的设计方法 | |
Khayyatzadeh et al. | Sub-synchronous resonance damping using high penetration PV plant | |
CN102904266B (zh) | 一种确立风电场无功补偿容量适网性的方法 | |
CN109120196A (zh) | 一种水电机组调速***pid参数整定方法及其装置 | |
CN110212513A (zh) | 一种稳定直流微网母线电压波动的灵活虚拟电容控制方法 | |
CN109301814A (zh) | 一种接入电网风电容量分析方法和*** | |
CN109103927A (zh) | 提高一次调频动态响应特性的调速***pid控制器参数整定方法 | |
CN109149602A (zh) | 一种水电机组调速***附加阻尼控制器模型及其构建方法 | |
CN109412180A (zh) | 水电机组异步联网场景下的超低频振荡控制方法 | |
CN109217289A (zh) | 一种基于灰狼算法的光伏附加阻尼控制器参数优化方法 | |
CN109742773A (zh) | 一种自适应广域阻尼控制器 | |
CN110137948B (zh) | 一种双馈风电场次同步振荡风险的评估方法及装置 | |
Wang et al. | Control Method for Additional Damper in Hydro-turbine Speed Governor of Hydro-dominant Power Systems | |
Hinners et al. | Model order reduction of active distribution networks with TSO-DSO interconnection power flow control | |
CN107171344B (zh) | 一种基于极点配置抑制水电机组超低频振荡的调速器参数优化方法 | |
CN109638871A (zh) | 考虑风电接入的大规模交直流混联***主网划分方法 | |
CN109390972A (zh) | 水电为主电网异步互联后调速器参数调整方法及*** | |
CN105914752A (zh) | 基于快速搜索与密度峰值聚类的主导节点选择方法 | |
CN109066765A (zh) | 基于一致性策略的孤岛微电网自适应同步频率控制方法 | |
CN108879688A (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 |