CN112865140B - 一种针对能源***阻尼的量化分析方法及*** - Google Patents
一种针对能源***阻尼的量化分析方法及*** Download PDFInfo
- Publication number
- CN112865140B CN112865140B CN202110260466.5A CN202110260466A CN112865140B CN 112865140 B CN112865140 B CN 112865140B CN 202110260466 A CN202110260466 A CN 202110260466A CN 112865140 B CN112865140 B CN 112865140B
- Authority
- CN
- China
- Prior art keywords
- slope
- energy system
- damping
- objective function
- value
- 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
- 238000013016 damping Methods 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000004445 quantitative analysis Methods 0.000 title claims abstract description 24
- 238000010586 diagram Methods 0.000 claims abstract description 39
- 238000011478 gradient descent method Methods 0.000 claims abstract description 22
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 41
- 238000013178 mathematical model Methods 0.000 claims description 28
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 description 8
- 230000010355 oscillation Effects 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 6
- 230000003321 amplification Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001360 synchronised effect Effects 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
-
- 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)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种针对能源***阻尼的量化分析方法及***,包括:确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;输入变量为能源***转速偏差,输出变量为能源***功率变化量;能源***为任意阶线性***或非线性***;确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼。本发明适用的能源***可以是任意阶线性或非线性***。
Description
技术领域
本发明属于能源***分析领域,更具体地,涉及一种针对能源***阻尼的量化分析方法及***。
背景技术
我国是水资源开发利用大国,已在西南地区的金沙江、雅砻江等河流建设了多座巨型水电站。随着电网多电源点并列运行、多能互补协同运行格局的形成,间歇性能源的随机波动特性给电网调度及安全带来了新的挑战。以低频振荡为代表的动稳定问题成为影响电网安全运行的重要因素。对于低频振荡的机理研究,阻尼转矩法是目前较为成熟的方法,通过分析不同机组的阻尼特性,能够对低频振荡做出合理解释。
然而,目前能源***阻尼的量化分析方法多以简单低阶模型为基础,通过经典或现代控制理论的方法计算***阻尼。尽管这样得到的阻尼是精确的,但由于模型经过了简化和线性化处理,导致所得的阻尼无法精确反映***特性,仅能作为一种定性分析的结果。此外,随着***复杂性的增加,即***阶次升高和非线性环节的引入,***阻尼的计算难度大大增加,使得传统方法不再适用于求解***阻尼。因此,需要研究一种针对复杂能源***阻尼的量化分析方法,以克服传统方法的局限性。
发明内容
针对现有技术的缺陷,本发明的目的在于提供一种针对能源***阻尼的量化分析方法及***,旨在解决目前能源***阻尼的量化分析方法多以简单低阶模型为基础,无法精确反映***特性,且随着***复杂性的增加,***阻尼的计算难度大大增加,传统方法不再适用于求解***阻尼的问题。
为实现上述目的,第一方面,本发明提供了一种针对能源***阻尼的量化分析方法,包括如下步骤:
确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼。
可选地,所述能源***的输入变量波形和对应的输出变量波形通过如下方式得到:
搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或
运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
可选地,所述目标函数的公式为:
可选地,所述利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:
Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;
Step 2:计算初始斜率k0对应的目标函数值f0;
Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:
Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;
Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step 6;
Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:
ki=kbest-di·gbest
式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;
Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step2和Step3;
Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
第二方面,本发明提供了一种针对能源***阻尼的量化分析***,包括:
变量确定单元,用于确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
直线确定单元,用于确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
阻尼分析单元,用于将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼。
可选地,所述变量确定单元,搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
可选地,所述阻尼分析单元利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;Step 2:计算初始斜率k0对应的目标函数值f0;Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step6;Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:ki=kbest-di·gbest;式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step 2和Step3;Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
本发明提供一种针对能源***阻尼的量化分析方法及***,采用一种基于相轨迹图的目标函数,利用梯度下降法的原理找到***的最优阻尼。相比于传统方法,该方法不受***复杂程度的限制,同时表达直观、计算量小、计算精度高,能够得到任意能源***阻尼的精确估计值,实现对***阻尼特性的量化分析。
附图说明
图1是本发明提供的一种针对能源***阻尼的量化分析方法流程图;
图2是本发明提供的一种简化的水电机组传递函数方框图;
图3是本发明提供的一种简化的火电机组传递函数方框图;
图4是本发明提供的相轨迹图和穿越相轨迹图的直线的示意图;
图5是本发明方法和传统方法计算的图2中水电机组的阻尼对比图;
图6是本发明方法和传统方法计算的图3中火电机组的阻尼对比图;
图7是本发明提供的一种针对能源***阻尼的量化分析***架构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明所要解决的技术问题是克服现有技术难以计算复杂能源***阻尼的缺陷,提供一种针对复杂能源***阻尼的量化分析方法,从而在量化任意能源***阻尼时,无需考虑因***过于复杂而导致阻尼难以计算的问题。本发明适用的能源***可以是任意阶线性或非线性***。
图1是本发明提供的一种针对能源***阻尼的量化分析方法流程图;如图1所示,包括如下步骤:
S101,确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
S102,确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
S103,将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼。
可选地,能源***的输入变量波形和对应的输出变量波形通过如下方式得到:
搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或
运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
可选地,所述目标函数的公式为:
可选地,所述利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:
Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;
Step 2:计算初始斜率k0对应的目标函数值f0;
Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:
Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;
Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step 6;
Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:
ki=kbest-di·gbest
式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;
Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step2和Step3;
Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
可以理解的是,本发明方法不限于能源***的阶型,适用于任意能够获取能源***输入输出变量的***,其中,输入变量为转速或频率偏差,输出变量为功率变化量。且本发明可以通过搭建能源***的数学模型,数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或直接运行要分析阻尼的能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
为说明本发明效果,因为现有阻尼分析方法仅适用于简单的能源***,下面以两种不同的简单能源***作为本发明和现有技术的实施对象对本发明方法的性能进行详细说明:一种为水电机组,其简化的传递函数方框图如图2所示;一种为火电机组,其简化的传递函数方框图如图3所示。
需要说明的是,如图2和图3所示,现有能源***阻尼分析方法需要提前知道能源***的传递函数,具有局限性。而本发明则不需要知道能源***的传递函数,相比现有技术具有极强的优势。
对于如图2所示的水电机组,利用传统方法量化其阻尼的步骤为:
步骤(1):结合***各部分传递函数,整理化简得到整体传递函数,即:
其中,ΔPm为功率变化量,Δw为转速偏差,Ka为放大倍数,Tg为液压***时间常数,Tw为水流惯性时间常数,s为拉普拉斯变量。
步骤(2):令s=jω,带入上式并经过一系列化简得:
ΔPm=DtΔw+KtΔδ
其中,Δw=sΔδ,DtΔw为阻尼转矩,KtΔδ为同步转矩,阻尼系数Dt表达式为:
其中,ω为***输入信号的震荡频率。当已知***参数和振荡频率时,***阻尼可通过上式计算得到。
对于如图2所示的水电机组,利用本申请方法量化其阻尼的步骤为:
步骤(1):在MATLAB/Simulink中搭建如图2中所示的水电机组仿真模型;
步骤(2):采用正弦信号(例如幅值为0.1Hz,频率为0.1Hz)作为转速偏差输入模型进行仿真,以模拟机组转速或频率的低频振荡现象。记录***的输入变量(即转速偏差)和输出变量(即功率变化量)波形,并绘制相应的转速-功率相轨迹图;其中,输出变量波形和输入变量波形也可通过对水电机组模型对应的实体水电机组直接进行测试得到,不需要步骤(1),也不需要进行步骤(2)的仿真步骤。
步骤(3):将上述相轨迹图的中心平移至坐标系原点,以便从所述相轨迹图的众多穿越直线中确定唯一一条具有最小目标函数值的直线,其斜率的相反数即为所述能源***的阻尼。目标函数的计算方法如下:
步骤(4):基于梯度下降法寻找步骤(3)中具有最小目标函数值的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:
Step 1:设置寻优过程中参数的初始值,即初始步长d0=0.1和初始斜率k0=0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;
Step 2:计算初始斜率对应的目标函数值f0;
Step 3:在初始斜率的某一临近点,即k0+ε,ε=10-10处,重新计算对应的目标函数值f′0。根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:
Step 4:将首次计算得到的斜率、目标函数及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;
Step 5:判断当前迭代次数i是否达到100或本次得到的最优斜率与上次之差是否小于10-10,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step 6;
Step 6:实施第i次迭代计算。基于梯度下降法的原理,计算新的斜率值ki:
ki=kbest-di·gbest
式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;
Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step2和Step3;
Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以2),将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以0.5,并返回Step 6。
对于如图3所示的火电机组,利用传统方法和本申请方法量化其阻尼的步骤与图2所示的水电机组相同。
图4本发明提供的相轨迹图和穿越相轨迹图的直线的示意图,如图4所示,对得到的转速和功率时域波形进行合成,得到转速-功率相轨迹图,具体方法是:对时域波形按时间从小到大进行取点,每一时刻的转速和功率对应相轨迹图上的一点,将所有点按时间顺序连接即可得到对应的相轨迹图。将相轨迹图中心平移至坐标原点后,其穿越直线为具有不同斜率k的正比例函数y=k·x,x和y为正比例函数的自变量和因变量。对于相轨迹图上某点的纵坐标y1i,均对应同一横坐标值下直线上一点的纵坐标y2i,它们的差值定义为Xi。
为了说明本申请方法的有效性和精度,在不同的***参数和振荡频率下(即0Hz~2.5Hz)对比本发明和传统方法的阻尼计算结果。对于如图2所示的水电机组,不变的***参数为放大倍数(Ka=20)和液压***时间常数(Tg=2.5),变化的***参数为水流惯性时间常数,即Tw=1.0、2.0和4.0,两种方法的计算结果对比如图5所示。对于如图3所示的火电机组,不变的***参数为放大倍数(Ka=20)和液压***时间常数(Tg=0.3),变化的***参数为高压蒸汽容积时间常数,即TCH=0.2、0.4和0.8,两种方法的计算结果对比如图6所示。图5和图6的对比结果表明,采用本发明所述方法计算得到的阻尼与传统方法高度一致,采用本发明所述方法计算复杂能源***阻尼是有效的,能够克服传统方法在求解复杂***阻尼时的局限性,且可以适用于任意阶线性或非线性***,并具有较高的计算精度。
图7是本发明提供的一种针对能源***阻尼的量化分析***架构图,如图7所示,包括:
变量确定单元710,用于确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
直线确定单元720,用于确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
阻尼分析单元730,用于将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼。
可选地,所述变量确定单元710,搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
可选地,阻尼分析单元730利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;Step 2:计算初始斜率k0对应的目标函数值f0;Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step6;Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:ki=kbest-di·gbest;式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step 2和Step3;Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
具体地,各个单元的功能可参见前述方法实施例的介绍,在此不做赘述。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种针对能源***阻尼的量化分析方法,其特征在于,包括如下步骤:
确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼;所述量化分析方法能够得到任意能源***阻尼的精确估计值。
2.根据权利要求1所述的针对能源***阻尼的量化分析方法,其特征在于,所述能源***的输入变量波形和对应的输出变量波形通过如下方式得到:
搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或
运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
4.根据权利要求3所述的针对能源***阻尼的量化分析方法,其特征在于,所述利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:
Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;
Step 2:计算初始斜率k0对应的目标函数值f0;
Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:
Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;
Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step 6;
Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:
ki=kbest-di·gbest
式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;
Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step2和Step3;
Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
5.一种针对能源***阻尼的量化分析***,其特征在于,包括:
变量确定单元,用于确定能源***的输入变量波形和对应的输出变量波形,并合成输入变量和输出变量之间的相轨迹图;所述输入变量为能源***转速偏差,输出变量为能源***功率变化量;所述能源***为任意阶线性***或非线性***;
直线确定单元,用于确定能够穿越所述相轨迹图的多条直线,针对每条直线求取对应的数列,其中,所述数列由相同横坐标下相轨迹图上每点的纵坐标与直线上对应点的纵坐标的差值组成;
阻尼分析单元,用于将每条直线对应数列的均方误差作为目标函数,利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼;
所述量化分析***能够得到任意能源***阻尼的精确估计值。
6.根据权利要求5所述的针对能源***阻尼的量化分析***,其特征在于,所述变量确定单元,搭建所述能源***的数学模型,所述数学模型的输入变量为转速偏差,输出变量为功率变化量;采用正弦信号作为输入对该数学模型进行仿真,将仿真得到的输入变量波形和对应的输出变量波形记录;或运行能源***,改变输入变量,记录对应的输入变量波形和相应的输出变量波形。
8.根据权利要求7所述的针对能源***阻尼的量化分析***,其特征在于,所述阻尼分析单元利用梯度下降法找到目标函数的数值最小时对应的直线,并将该直线斜率的相反数作为所述能源***的阻尼,具体包括如下步骤:Step 1:确定初始步长d0和初始斜率k0,代计数器i=0;其中,步长d0表示斜率变化的快慢程度,斜率k0表示能够穿越相轨迹图的直线的斜率,i表示当前迭代次数;Step 2:计算初始斜率k0对应的目标函数值f0;Step 3:在初始斜率k0的某一临近点k0+ε处,ε为预设邻域范围,重新计算对应的目标函数值f′0;根据两次计算结果,利用斜率公式计算初始斜率对应的目标函数一阶导数g0,计算公式为:Step 4:将首次计算得到的斜率、目标函数值及其一阶导数保存为历史最优值,即kbest=k0、fbest=f0、gbest=g0;其中,kbest表示历史最优斜率,fbest表示历史最优目标函数值,gbest表示历史最优目标函数一阶导数值;Step 5:判断当前迭代次数i是否达到最大迭代次数或本次得到的最优斜率与上次之差是否小于预设值,若满足其中一个条件,输出kbest,并将kbest的相反数作为能源***的阻尼;否则进入Step 6;Step 6:实施第i次迭代计算,基于梯度下降法的原理,计算新的斜率值ki:ki=kbest-di·gbest;式中,di为当前步长,i=0时,当前步长是初始步长d0,i>0时,当前步长则根据Step8更新;Step7:计算新斜率ki对应的目标函数fi及其一阶导数gi,计算方法同Step 2和Step3;Step8:比较新斜率ki和历史最优斜率kbest的目标函数值,若前者小于后者,即fi<fbest,则更新历史最优值,即kbest=ki、fbest=fi、gbest=gi,同时将当前步长乘以一个大于1的实数,将迭代次数i加1,并返回Step 5;否则,保持历史最优值不变,将当前步长乘以一个小于1的正实数,并返回Step 6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110260466.5A CN112865140B (zh) | 2021-03-10 | 2021-03-10 | 一种针对能源***阻尼的量化分析方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110260466.5A CN112865140B (zh) | 2021-03-10 | 2021-03-10 | 一种针对能源***阻尼的量化分析方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112865140A CN112865140A (zh) | 2021-05-28 |
CN112865140B true CN112865140B (zh) | 2022-09-30 |
Family
ID=75995213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110260466.5A Active CN112865140B (zh) | 2021-03-10 | 2021-03-10 | 一种针对能源***阻尼的量化分析方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112865140B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004038562A2 (en) * | 2002-10-23 | 2004-05-06 | Hynomics Corporation | Methods and system for generic optimization of control functions |
JP2007325359A (ja) * | 2006-05-30 | 2007-12-13 | Central Res Inst Of Electric Power Ind | 電力系統の制御系定数設定方法、装置及びプログラム |
CN102545247A (zh) * | 2012-02-27 | 2012-07-04 | 山东大学 | 基于广域测量信号的风电场侧阻尼控制器参数整定方法 |
JP2017158317A (ja) * | 2016-03-02 | 2017-09-07 | 北陸電力株式会社 | 周波数フィードバック方式のハンチング現象の解析方法、および回避方法 |
CN111654038A (zh) * | 2020-05-20 | 2020-09-11 | 中国电力科学研究院有限公司 | 一种基于轨迹特征识别电力***暂态失稳的方法及*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200074339A1 (en) * | 2013-12-11 | 2020-03-05 | U.S. Government as represented by Secretary of the Air Force | Fano-Based Information Theoretic Method (FBIT) for Design and Optimization of Nonlinear Systems |
CN104600756B (zh) * | 2015-01-29 | 2015-09-16 | 华中科技大学 | 一种中小型水力发电机组的机群等值建模方法 |
CN106451436B (zh) * | 2016-11-04 | 2019-03-12 | 广东电网有限责任公司电力科学研究院 | 一种综合量化评估发电机组动态特性的方法及*** |
CN110867903B (zh) * | 2019-11-27 | 2021-05-18 | 浙江工业大学 | 一种基于深度神经网络的智能电网***实时最优功率分配方法 |
CN111682527B (zh) * | 2020-05-28 | 2022-07-05 | 国网浙江省电力有限公司湖州供电公司 | 一种考虑负荷阻尼的在线机组频率稳定性优化解析方法 |
-
2021
- 2021-03-10 CN CN202110260466.5A patent/CN112865140B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004038562A2 (en) * | 2002-10-23 | 2004-05-06 | Hynomics Corporation | Methods and system for generic optimization of control functions |
JP2007325359A (ja) * | 2006-05-30 | 2007-12-13 | Central Res Inst Of Electric Power Ind | 電力系統の制御系定数設定方法、装置及びプログラム |
CN102545247A (zh) * | 2012-02-27 | 2012-07-04 | 山东大学 | 基于广域测量信号的风电场侧阻尼控制器参数整定方法 |
JP2017158317A (ja) * | 2016-03-02 | 2017-09-07 | 北陸電力株式会社 | 周波数フィードバック方式のハンチング現象の解析方法、および回避方法 |
CN111654038A (zh) * | 2020-05-20 | 2020-09-11 | 中国电力科学研究院有限公司 | 一种基于轨迹特征识别电力***暂态失稳的方法及*** |
Non-Patent Citations (2)
Title |
---|
"Stability Analysis of Dissipative Systems Subject to Nonlinear Damping via Lyapunov Techniques";Swann Marx 等;《 IEEE Transactions on Automatic Control 》;20190826;全文 * |
"非线性转子***稳定性量化分析方法";郑惠萍 等;《应用数学和力学》;20050915;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112865140A (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | A nonlinear generalized predictive control for pumped storage unit | |
CN103887815A (zh) | 基于运行数据的风电场参数辨识及动态等值方法 | |
US9739262B2 (en) | Static testing and calibrating method for PID link of control system of wind turbine | |
KR101474102B1 (ko) | 임펠러 형상 최적설계방법 | |
CN115313524A (zh) | 一种基于构网型变流器的光伏发电并网控制方法及*** | |
Doubrawa et al. | Optimization-based calibration of FAST. Farm parameters against large-eddy simulations | |
CN106950857A (zh) | 基于模糊逻辑控制的光伏电池mppt仿真方法 | |
CN106300338A (zh) | 基于轨迹灵敏度的受端电网动态频率安全量化评估方法 | |
CN111046327A (zh) | 适用于低频振荡与次同步振荡辨识的Prony分析方法 | |
CN111159941A (zh) | 一种用于汽车液力变矩器内流场瞬态数值模拟的方法 | |
Santamaría et al. | Novel methodology for performance characterization of vertical axis wind turbines (VAWT) prototypes through active driving mode | |
CN112865140B (zh) | 一种针对能源***阻尼的量化分析方法及*** | |
Rice et al. | Identification of the stability margin between safe operation and the onset of blade flutter | |
CN107247994B (zh) | 一种托盘塔脱硫装置脱硫效率的模糊建模方法 | |
CN104038128A (zh) | 基于nurbs的直线电机推力波动控制方法 | |
CN105116959A (zh) | 基于最优梯度变步长的最大功率点跟踪控制与预测方法 | |
CN103399489A (zh) | 一种电液加速度伺服***的波形再现控制方法 | |
CN111237116B (zh) | 非线性隐极水轮发电机水门与励磁联合调节方法和*** | |
Park et al. | A novel VLSI architecture for the real-time implementation of 2-D signal processing systems | |
Gunawan et al. | Model validation using experimental measurements from the Garfield Thomas water tunnel at the Applied Research Laboratory (ARL) at Penn State University | |
Lu et al. | Development of an integrated system design tool for helical Vertical Axis Wind Turbines (VAWT-X) | |
CN109989875B (zh) | 一种水轮机调节***参数辨识方法 | |
He et al. | Nonlinear Correction and Temperature Compensation Method of Turbine Flowmeter Based on Neural Network | |
CN116191478B (zh) | 一种多个异步电机等效惯量评估和频率响应建模方法 | |
CN110780607B (zh) | 基于adpss的水轮机调速***阻尼测试方法及装置 |
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 |