CN108763831B - 一种水泵水轮机线性模型参数辨识方法 - Google Patents

一种水泵水轮机线性模型参数辨识方法 Download PDF

Info

Publication number
CN108763831B
CN108763831B CN201810662132.9A CN201810662132A CN108763831B CN 108763831 B CN108763831 B CN 108763831B CN 201810662132 A CN201810662132 A CN 201810662132A CN 108763831 B CN108763831 B CN 108763831B
Authority
CN
China
Prior art keywords
model
storage unit
pumped storage
coefficient
linear model
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
Application number
CN201810662132.9A
Other languages
English (en)
Other versions
CN108763831A (zh
Inventor
彭绪意
杨文�
聂赛
杨雄
洪云来
常国庆
莫旭晶
刘泽
胥千鑫
汤凯
秦程
章志平
温锦红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIANGXI HONGPING PUMPED STORAGE CO LTD
State Grid Corp of China SGCC
State Grid Xinyuan Co Ltd
Original Assignee
JIANGXI HONGPING PUMPED STORAGE CO LTD
State Grid Corp of China SGCC
State Grid Xinyuan Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by JIANGXI HONGPING PUMPED STORAGE CO LTD, State Grid Corp of China SGCC, State Grid Xinyuan Co Ltd filed Critical JIANGXI HONGPING PUMPED STORAGE CO LTD
Priority to CN201810662132.9A priority Critical patent/CN108763831B/zh
Publication of CN108763831A publication Critical patent/CN108763831A/zh
Application granted granted Critical
Publication of CN108763831B publication Critical patent/CN108763831B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/044—Recurrent networks, e.g. Hopfield networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • 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/16—Mechanical energy storage, e.g. flywheels or pressurised fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Geometry (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Evolutionary Biology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Control Of Water Turbines (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种水泵水轮机线性模型参数辨识方法，它包括以下步骤：由引水***近似弹性水击模型、水泵水轮机IEEE六参数模型、发电电动机***模型和PID调速器模型耦合得到发电工况下的抽水蓄能机组调速***线性模型，构造抽水蓄能机组调速***的线性映射模型；确定抽水蓄能机组调速***待辨识参数的上下边界，并获取参数辨识样本数据；采用BP神经网络方法训练参数辨识样本数据，建立BP神经网络参数辨识模型；以抽水蓄能机组调速***的运行实测数据作为BP神经网络参数辨识模型的输入，求解获得抽水蓄能机组调速***的参数辨识结果。本发明不仅较容易确定线性模型参数、数据的精确度较理想，还具有样本的获取难度较小、选取样本的灵活性较好的优点。

Description

一种水泵水轮机线性模型参数辨识方法
技术领域
本发明属于抽水蓄能机组精确建模领域，具体涉及一种水泵水轮机线性模型参数辨识方法。
背景技术
水泵水轮机是抽水蓄能机组的核心设备，是机组调速***中的控制对象，针对抽水蓄能机组调速***动力学响应机理、控制优化以及故障诊断等相关研究的基础是***模型的精确描述。模型参数辨识是解决机组调速***模型精确表达的有效途径，也是相关领域的研究热点，方法研究也逐步趋于成熟。国内外学术与工程界对抽水蓄能机组控制***参数辨识方法可以分为：***线性模型辨识和非线性模型辨识。
抽水蓄能机组调速***线性模型具有模型结构简单、计算实时性好等特点，其一直是基于模型的仿真装置开发、复杂电力***稳定性研究的首选。线性模型参数的确定一直是困扰工程界的难题，传统的方法是根据水轮机综合特性曲线或者水泵水轮机全特性曲线，在某一稳定工况点进行线性化处理计算出能线性模型的六个系数，线性模型参数较难确定。但是，随着机组的长时间运行或在发生重大检修后其工况运行边界会发生迁移，依据特性曲线计算的线性模型参数对机组实时运行状态的精确描述具有一定的局限，数据的精确度不理想。因此，现有的对抽水蓄能机组控制***参数辨识方法中存在着线性模型参数较难确定和数据的精确度不理想的问题。
发明内容
本发明的目的在于，提供一种水泵水轮机线性模型参数辨识方法。本发明不仅能较容易确定线性模型参数，还具有数据的精确度较理想的优点。
本发明的技术方案：一种水泵水轮机线性模型参数辨识方法，包括以下步骤：
a、由引水***近似弹性水击模型、水泵水轮机IEEE六参数模型、发电电动机***模型和PID调速器模型耦合得到发电工况下的抽水蓄能机组调速***线性模型，再基于模型映射理论，构造抽水蓄能机组调速***的线性映射模型，线性映射模型的输出即为抽水蓄能机组调速***待辨识参数；
b、确定抽水蓄能机组调速***待辨识参数的上下边界，并获取参数辨识样本数据；
c、采用BP神经网络方法训练参数辨识样本数据，建立基于调速***线性模型映射的BP神经网络参数辨识模型；
d、以抽水蓄能机组调速***的运行实测数据作为BP神经网络参数辨识模型的输入，求解获得抽水蓄能机组调速***的参数辨识结果。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤a中的PID调速器模型包括并联PID控制器和辅助接力器型执行机构；所述并联PID控制器传递函数表达式为
Figure BDA0001706643910000021
其中s为拉普拉斯算子，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，xc为机组转速给定值，x为机组转速，yc为接力器开度给定值，bp为永态差值系数，σ为调速器电气控制输出；
所述辅助接力器型执行机构传递函数表达式为
Figure BDA0001706643910000022
其中Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，k0为放大系数，y为执行机构输出，s为拉普拉斯算子，σ为调速器电气控制输出。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤a中的水泵水轮机IEEE六参数模型采用相对值表达式为
Figure BDA0001706643910000031
和
Figure BDA0001706643910000032
其中mt为力矩偏差相对值，q为流量偏差相对值，y为导叶开度偏差相对值，x为转速偏差相对值，h为水头偏差相对值，ey为力矩对导叶开度相对系数，ex为力矩对转速传递系数，eh为力矩对工作水头传递系数，eqy为流量对导叶开度相对系数，eqx为流量对转速相对系数，eqh为流量对工作水头相对系数，f为摩阻系数。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤a中的引水***近似弹性水击模型的传递函数表达式为
Figure BDA0001706643910000033
其中hw为管道特性系数，Tr为水击相长，f为摩阻系数，s为拉普拉斯算子。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤a中的发电电动机***模型采用转子运动方程描述发电电动机动态特性，发电电动机***模型的传递函数表达式为
Figure BDA0001706643910000034
其中Ta为机组惯性时间常数，eg为发电机自调节系数，s为拉普拉斯算子。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤a中的抽水蓄能机组调速***线性模型中，选取X＝[x,σi,x1,ya,y,x2,x3]T为状态变量，U＝[mg0,xc,yc]T为输入变量，Y＝[x]T为输出变量，其中x为转速偏差相对值，x1为微分环节的状态变量，x2、x3均为引水***环节的状态变量，则调速***线性模型的状态空间表达式为
Figure BDA0001706643910000041
其中
Figure BDA0001706643910000042
Figure BDA0001706643910000043
a22＝bp·Ki，a23＝1，
Figure BDA0001706643910000044
b23＝-Ki·bp；
Figure BDA0001706643910000045
Figure BDA0001706643910000046
Figure BDA0001706643910000047
a67＝1，a71＝eqx，a75＝eqy，
Figure BDA0001706643910000048
Figure BDA0001706643910000049
ex，eqx，ey，eqy，eh，eqh均为水泵水轮机线性模型待辨识参数，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，bp为永态差值系数，Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，Tr为水击相长，Ta为机组惯性时间常数；
所述抽水蓄能机组调速***线性模型的***输出表达式为
Figure BDA0001706643910000051
其中
Figure BDA0001706643910000052
ex、eh均为水泵水轮机线性模型待辨识参数，Ta为机组惯性时间常数，ey为力矩对导叶开度相对系数；
对于抽水蓄能机组调速***线性模型的***输出表达式，当
Figure BDA0001706643910000053
中***的参数选定时，抽水蓄能机组调速***能够作为Y与U的单变量***，一组模型输入变量U会对应一个***响应输出向量Y，提取模型的输入变量U，则模型表示为
Figure BDA0001706643910000054
将
Figure BDA0001706643910000055
定义为源模型；
基于映射理论象与原象直接的对应关系，当***输入U不变时，定义由Y到
Figure BDA0001706643910000056
的变维映射过程为映射模型，其表示式为
Figure BDA0001706643910000057
即映射模型的输出为抽水蓄能机组调速***待辨识参数。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤b中的确定抽水蓄能机组调速***待辨识参数的上下边界，对于水泵水轮机线性模型待辨识参数[ex,eqx,ey,eqy,eh,eqh]，由机组运行实测数据得其对应的机组运行工作水头Hw，当Hw为定值时，机组的流量和力矩变化的因变量则为导叶开度y，保持机组单位转速为额定单位转速n11r，将导叶开度在零开度到额定开度区间内，以5％的偏移量计算20组水泵水轮机线性***模型六系数，从中选取相应的最大值和最小值作为抽水蓄能机组调速***待辨识参数的上下边界值。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤b中的获取参数辨识样本数据为，先基于建立的***源模型，保持源模型中的输入信号U不变，以相同的偏移量增加或减小
Figure BDA0001706643910000058
中的待辨识参数，并将待辨识参数向量xi代入源模型进行仿真，获得模型输出Yi；其次，重复上述操作，直到获得参数辨识所需的样本数量；接下来，为了增强神经网络的学习能力，计算调速***输出Yi的4个指标量包括：上升时间、调节时间、超调量和波峰次数，将指标量与Yi构成新的Yi，并将xi与Yi进行组合，得到BP神经网络辨识所需的样本集{(x1,Y1),(x2,Y2),...,(xi,Yi),...,(xN,YN)}，N为样本集中的样本数；最后，从获得的样本集中随机抽取部分样本作为训练样本，其余为测试样本。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤c中基于训练样本为，采用BP神经网络方法训练逼近抽水蓄能机组线性模型的映射模型，并用测试样本集验证逼近的精度，待精度符合辨识要求时，即建立了基于调速***线性模型映射的BP神经网络参数辨识模型。
前述的一种水泵水轮机线性模型参数辨识方法中，所述步骤d中以抽水蓄能机组调速***在空载开机工况的原始运行实测数据作为步骤c中建立的BP神经网络参数辨识模型的输入，即得到一组参数向量x，参数向量x为***的参数辨识结果。
与现有技术相比，本发明改进了现有的对抽水蓄能机组控制***参数辨识方法，通过利用BP神经网络训练样本集，并将测试样本输入训练完成后的网络，得到的网络训练截断误差为10-7时的误差正态检验结果以及测试样本的平均参数误差，BP神经网络的输出值与样本值之间的线性拟合曲线的方差值R大于0.99，且与正态分布的检验直线吻合度高。进一步，由训练样本的输出与期望输出的平均参数误差分布图可以看出，神经网络训练结果的分布误差均小于4％，网络训练结果稳定、可靠，线性模型参数较容易确定，且数据的精确度较理想；将机组空载开机工况的实际运行数据作为输入，依据建立的辨识模型，得到参数辨识结果。基于参数辨识结果的机组转速仿真输出与实测数据对比所示，结果表明，依照参数辨识结果的调速***机组转速仿真曲线与实测数据的吻合度高，且调速***的负荷变化受电网调度部门的约束和限制较小，对机组任意施加激励获取足够的样本的获取难度较小。本发明提供的基于抽水蓄能机组调速***线性映射模型的BP神经网络水泵水轮机参数辨识方法，能够克服神经网络方法在调速***参数辨识时训练样本不足的难题，且该方法改善了神经网络参数辨识对样本选取的灵活性，避免了对大量试验数据的依赖，实现了机组实际运行数据微样本集的调速***线性模型的高精度参数辨识，选取样本的灵活性较好。因此，本发明不仅较容易确定线性模型参数、数据的精确度较理想，还具有样本的获取难度较小、选取样本的灵活性较好的优点。
附图说明
图1是抽水蓄能机组调速***线性模型的框架图；
图2是空载开机工况正态误差分布示意图；
图3是空载开机工况测试样本误差示意图；
图4是基于参数辨识结果的机组转速仿真输出与实测数据对比示意图。
附图中的标记为：1-PID调速器模型，2-并联PID控制器，3-辅助接力器型执行机构，4-水泵水轮机IEEE六参数模型，5-发电电动机***模型。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明，但并不作为对本发明限制的依据。
实施例一。一种水泵水轮机线性模型参数辨识方法，包括以下步骤：
a、由引水***近似弹性水击模型、水泵水轮机IEEE六参数模型、发电电动机***模型和PID调速器模型耦合得到发电工况下的抽水蓄能机组调速***线性模型，如图1所示，再基于模型映射理论，构造抽水蓄能机组调速***的线性映射模型，线性映射模型的输出即为抽水蓄能机组调速***待辨识参数；
b、确定抽水蓄能机组调速***待辨识参数的上下边界，并获取参数辨识样本数据；
c、采用BP神经网络方法训练参数辨识样本数据，建立基于调速***线性模型映射的BP神经网络参数辨识模型；
d、以抽水蓄能机组调速***的运行实测数据作为BP神经网络参数辨识模型的输入，求解获得抽水蓄能机组调速***的参数辨识结果。
所述步骤a中的PID调速器模型包括并联PID控制器和辅助接力器型执行机构；所述并联PID控制器传递函数表达式为
Figure BDA0001706643910000081
其中s为拉普拉斯算子，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，xc为机组转速给定值，x为机组转速，yc为接力器开度给定值，bp为永态差值系数，σ为调速器电气控制输出；
所述辅助接力器型执行机构传递函数表达式为
Figure BDA0001706643910000082
其中Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，k0为放大系数，y为执行机构输出，s为拉普拉斯算子，σ为调速器电气控制输出。
所述步骤a中的水泵水轮机IEEE六参数模型采用相对值表达式为
Figure BDA0001706643910000083
和
Figure BDA0001706643910000091
其中mt为力矩偏差相对值，q为流量偏差相对值，y为导叶开度偏差相对值，x为转速偏差相对值，h为水头偏差相对值，ey为力矩对导叶开度相对系数，ex为力矩对转速传递系数，eh为力矩对工作水头传递系数，eqy为流量对导叶开度相对系数，eqx为流量对转速相对系数，eqh为流量对工作水头相对系数，f为摩阻系数。
所述步骤a中的引水***近似弹性水击模型的传递函数表达式为
Figure BDA0001706643910000092
其中hw为管道特性系数，Tr为水击相长，f为摩阻系数，s为拉普拉斯算子。
所述步骤a中的发电电动机***模型采用转子运动方程描述发电电动机动态特性，发电电动机***模型的传递函数表达式为
Figure BDA0001706643910000093
其中Ta为机组惯性时间常数，eg为发电机自调节系数，s为拉普拉斯算子。
所述步骤a中的抽水蓄能机组调速***线性模型中，选取X＝[x,σi,x1,ya,y,x2,x3]T为状态变量，U＝[mg0,xc,yc]T为输入变量，Y＝[x]T为输出变量，其中x为转速偏差相对值，x1为微分环节的状态变量，x2、x3均为引水***环节的状态变量，则调速***线性模型的状态空间表达式为
Figure BDA0001706643910000101
其中
Figure BDA0001706643910000102
Figure BDA0001706643910000103
a22＝bp·Ki，a23＝1，
Figure BDA0001706643910000104
Figure BDA0001706643910000105
Figure BDA0001706643910000106
a67＝1，a71＝eqx，a75＝eqy，
Figure BDA0001706643910000107
Figure BDA0001706643910000108
ex，eqx，ey，eqy，eh，eqh均为水泵水轮机线性模型待辨识参数，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，bp为永态差值系数，Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，Tr为水击相长，Ta为机组惯性时间常数；
所述抽水蓄能机组调速***线性模型的***输出表达式为
Figure BDA0001706643910000109
其中
Figure BDA00017066439100001010
ex、eh均为水泵水轮机线性模型待辨识参数，Ta为机组惯性时间常数，ey为力矩对导叶开度相对系数；
对于抽水蓄能机组调速***线性模型的***输出表达式，当
Figure BDA0001706643910000111
中***的参数选定时，抽水蓄能机组调速***能够作为Y与U的单变量***，一组模型输入变量U会对应一个***响应输出向量Y，提取模型的输入变量U，则模型表示为
Figure BDA0001706643910000112
将
Figure BDA0001706643910000113
定义为源模型；
基于映射理论象与原象直接的对应关系，当***输入U不变时，定义由Y到
Figure BDA0001706643910000114
的变维映射过程为映射模型，其表示式为
Figure BDA0001706643910000115
即映射模型的输出为抽水蓄能机组调速***待辨识参数。
所述步骤b中的确定抽水蓄能机组调速***待辨识参数的上下边界，对于水泵水轮机线性模型待辨识参数[ex,eqx,ey,eqy,eh,eqh]，由机组运行实测数据得其对应的机组运行工作水头Hw，当Hw为定值时，机组的流量和力矩变化的因变量则为导叶开度y，保持机组单位转速为额定单位转速n11r，将导叶开度在零开度到额定开度区间内，以5％的偏移量计算20组水泵水轮机线性***模型六系数，从中选取相应的最大值和最小值作为抽水蓄能机组调速***待辨识参数的上下边界值。
所述步骤b中的获取参数辨识样本数据为，先基于建立的***源模型，保持源模型中的输入信号U不变，以相同的偏移量增加或减小
Figure BDA0001706643910000116
中的待辨识参数，并将待辨识参数向量xi代入源模型进行仿真，获得模型输出Yi；其次，重复上述操作，直到获得参数辨识所需的样本数量；接下来，为了增强神经网络的学习能力，计算调速***输出Yi的4个指标量包括：上升时间、调节时间、超调量和波峰次数，将指标量与Yi构成新的Yi，并将xi与Yi进行组合，得到BP神经网络辨识所需的样本集{(x1,Y1),(x2,Y2),...,(xi,Yi),...,(xN,YN)}，N为样本集中的样本数；最后，从获得的样本集中随机抽取部分样本作为训练样本，其余为测试样本。
所述步骤c中基于训练样本为，采用BP神经网络方法训练逼近抽水蓄能机组线性模型的映射模型，并用测试样本集验证逼近的精度，待精度符合辨识要求时，即建立了基于调速***线性模型映射的BP神经网络参数辨识模型。
所述步骤d中以抽水蓄能机组调速***在空载开机工况的原始运行实测数据作为步骤c中建立的BP神经网络参数辨识模型的输入，即得到一组参数向量x，参数向量x为***的参数辨识结果。
实施例二。以我国某单机300MW可逆式抽水蓄能机组的空载开机工况的实测数据作为辨识数据样本，辨识如图1所示的调速***线性模型中的水泵水轮机参数[ex,ey,eh,eqx,eqy,eqh]。
a、由引水***近似弹性水击模型、水泵水轮机IEEE六参数模型、发电电动机***模型和PID调速器模型耦合得到发电工况下的抽水蓄能机组调速***线性模型，再基于模型映射理论，构造抽水蓄能机组调速***的线性映射模型，线性映射模型的输出即为抽水蓄能机组调速***待辨识参数；
所述样本数据的工况条件：
1)上库水位735.45m,下库水位181.0m；
2)机组空载开机，PID参数为Kp＝0.6、Ki＝0.09、Kd＝0.5。
b、确定抽水蓄能机组调速***待辨识参数的上下边界，并获取参数辨识样本数据；
由机组运行实测数据计算得到参数边界如表1所示。
表1调速***线性模型水泵水轮机参数边界
Figure BDA0001706643910000121
c、采用BP神经网络方法训练参数辨识样本数据，建立基于调速***线性模型映射的BP神经网络参数辨识模型；
为保证样本的多样性和随机性，在确定的参数[ex,eqx,ey,eqy,eh,eqh]的边界域值内，随机生成100组参数xi。以实测数据对应的空载开机工况实际输入分别作为模型输入源模型中的输入信号U，保持U不变，设置仿真时间为70s，仿真时间间隔0.02s，将100组参数xi输入如图1所示的源模型，得到机组转速仿真输出Yi及其性能指标量。依据建立的映射模型，至此获得该工况条件下以Yi为输入、xi为目标输出、容量为100的样本集{(Y1,x1),(Y2,x2),...,(Yi,xi),...,(Y200,x200)}。依据随机抽样原则抽取80％的样本作为神经网络的训练样本(80个样本)，其余为测试样本(20个样本)。
d、以抽水蓄能机组调速***的运行实测数据作为BP神经网络参数辨识模型的输入，求解获得抽水蓄能机组调速***的参数辨识结果。
利用BP神经网络训练样本集，并将测试样本输入训练完成后的网络，得到的网络训练截断误差为10-7时的误差正态检验结果以及测试样本的平均参数误差如图3和图4所示。由图3和图4可知，BP神经网络的输出值与样本值之间的线性拟合曲线的方差值R大于0.99，且与正态分布的检验直线吻合度高。进一步，由训练样本的输出与期望输出的平均参数误差分布图可以看出，神经网络训练结果的分布误差均小于4％，网络训练结果稳定、可靠。
将机组空载开机工况的实际运行数据作为输入，依据建立的辨识模型，参数辨识结果如表2所示。
表2水泵水轮机的参数辨识结果
Figure BDA0001706643910000131
所述步骤a中的PID调速器模型包括并联PID控制器和辅助接力器型执行机构；所述并联PID控制器传递函数表达式为
Figure BDA0001706643910000141
其中s为拉普拉斯算子，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，xc为机组转速给定值，x为机组转速，yc为接力器开度给定值，bp为永态差值系数，σ为调速器电气控制输出；
所述辅助接力器型执行机构传递函数表达式为
Figure BDA0001706643910000142
其中Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，k0为放大系数，y为执行机构输出，s为拉普拉斯算子，σ为调速器电气控制输出。
所述步骤a中的水泵水轮机IEEE六参数模型采用相对值表达式为
Figure BDA0001706643910000143
和
Figure BDA0001706643910000144
其中mt为力矩偏差相对值，q为流量偏差相对值，y为导叶开度偏差相对值，x为转速偏差相对值，h为水头偏差相对值，ey为力矩对导叶开度相对系数，ex为力矩对转速传递系数，eh为力矩对工作水头传递系数，eqy为流量对导叶开度相对系数，eqx为流量对转速相对系数，eqh为流量对工作水头相对系数，f为摩阻系数。
所述步骤a中的引水***近似弹性水击模型的传递函数表达式为
Figure BDA0001706643910000151
其中hw为管道特性系数，Tr为水击相长，f为摩阻系数，s为拉普拉斯算子。
所述步骤a中的发电电动机***模型采用转子运动方程描述发电电动机动态特性，发电电动机***模型的传递函数表达式为
Figure BDA0001706643910000152
其中Ta为机组惯性时间常数，eg为发电机自调节系数，s为拉普拉斯算子。
所述步骤a中的抽水蓄能机组调速***线性模型中，选取X＝[x,σi,x1,ya,y,x2,x3]T为状态变量，U＝[mg0,xc,yc]T为输入变量，Y＝[x]T为输出变量，其中x为转速偏差相对值，x1为微分环节的状态变量，x2、x3均为引水***环节的状态变量，则调速***线性模型的状态空间表达式为
Figure BDA0001706643910000153
其中
Figure BDA0001706643910000154
Figure BDA0001706643910000155
a22＝bp·Ki，a23＝1，
Figure BDA0001706643910000156
b23＝-Ki·bp；
Figure BDA0001706643910000157
Figure BDA0001706643910000161
Figure BDA0001706643910000162
a67＝1，a71＝eqx，a75＝eqy，
Figure BDA0001706643910000163
Figure BDA0001706643910000164
ex，eqx，ey，eqy，eh，eqh均为水泵水轮机线性模型待辨识参数，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，bp为永态差值系数，Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，Tr为水击相长，Ta为机组惯性时间常数；
所述抽水蓄能机组调速***线性模型的***输出表达式为
Figure BDA0001706643910000165
其中
Figure BDA0001706643910000166
ex、eh均为水泵水轮机线性模型待辨识参数，Ta为机组惯性时间常数，ey为力矩对导叶开度相对系数；
对于抽水蓄能机组调速***线性模型的***输出表达式，当
Figure BDA0001706643910000167
中***的参数选定时，抽水蓄能机组调速***能够作为Y与U的单变量***，一组模型输入变量U会对应一个***响应输出向量Y，提取模型的输入变量U，则模型表示为
Figure BDA0001706643910000168
将
Figure BDA0001706643910000169
定义为源模型；
基于映射理论象与原象直接的对应关系，当***输入U不变时，定义由Y到
Figure BDA00017066439100001610
的变维映射过程为映射模型，其表示式为
Figure BDA00017066439100001611
即映射模型的输出为抽水蓄能机组调速***待辨识参数。
所述步骤b中的确定抽水蓄能机组调速***待辨识参数的上下边界，对于水泵水轮机线性模型待辨识参数[ex,eqx,ey,eqy,eh,eqh]，由机组运行实测数据得其对应的机组运行工作水头Hw，当Hw为定值时，机组的流量和力矩变化的因变量则为导叶开度y，保持机组单位转速为额定单位转速n11r，将导叶开度在零开度到额定开度区间内，以5％的偏移量计算20组水泵水轮机线性***模型六系数，从中选取相应的最大值和最小值作为抽水蓄能机组调速***待辨识参数的上下边界值。
所述步骤b中的获取参数辨识样本数据为，先基于建立的***源模型，保持源模型中的输入信号U不变，以相同的偏移量增加或减小
Figure BDA0001706643910000171
中的待辨识参数，并将待辨识参数向量xi代入源模型进行仿真，获得模型输出Yi；其次，重复上述操作，直到获得参数辨识所需的样本数量；接下来，为了增强神经网络的学习能力，计算调速***输出Yi的4个指标量包括：上升时间、调节时间、超调量和波峰次数，将指标量与Yi构成新的Yi，并将xi与Yi进行组合，得到BP神经网络辨识所需的样本集{(x1,Y1),(x2,Y2),...,(xi,Yi),...,(xN,YN)}，N为样本集中的样本数；最后，从获得的样本集中随机抽取部分样本作为训练样本，其余为测试样本。
所述步骤c中基于训练样本为，采用BP神经网络方法训练逼近抽水蓄能机组线性模型的映射模型，并用测试样本集验证逼近的精度，待精度符合辨识要求时，即建立了基于调速***线性模型映射的BP神经网络参数辨识模型。
所述步骤d中以抽水蓄能机组调速***在空载开机工况的原始运行实测数据作为步骤c中建立的BP神经网络参数辨识模型的输入，即得到一组参数向量x，参数向量x为***的参数辨识结果。
基于参数辨识结果的机组转速仿真输出与实测数据对比如图4所示，结果表明，依照参数辨识结果的调速***机组转速仿真曲线与实测数据的吻合度高。

Claims (7)

1.一种水泵水轮机线性模型参数辨识方法，其特征在于，包括以下步骤：
a、由引水***近似弹性水击模型、水泵水轮机IEEE六参数模型、发电电动机***模型和PID调速器模型耦合得到发电工况下的抽水蓄能机组调速***线性模型，再基于模型映射理论，构造抽水蓄能机组调速***的线性映射模型，线性映射模型的输出即为抽水蓄能机组调速***待辨识参数；
b、确定抽水蓄能机组调速***待辨识参数的上下边界，并获取参数辨识样本数据；
c、采用BP神经网络方法训练参数辨识样本数据，建立基于调速***线性模型映射的BP神经网络参数辨识模型；
d、以抽水蓄能机组调速***的运行实测数据作为BP神经网络参数辨识模型的输入，求解获得抽水蓄能机组调速***的参数辨识结果；
所述步骤a中的抽水蓄能机组调速***线性模型中，选取X＝[x，σi，x1，ya，y，x2，x3]T为状态变量，U＝[mg0，xc，yc]T为输入变量，Y＝[x]T为输出变量，其中x为转速偏差相对值，x1为微分环节的状态变量，x2、x3均为引水***环节的状态变量，则调速***线性模型的状态空间表达式为
Figure FDA0003533977270000011
其中
Figure FDA0003533977270000012
a22＝bp·Ki，a23＝1，
Figure FDA0003533977270000013
b23＝-Ki·bp；
Figure FDA0003533977270000014
Figure FDA0003533977270000015
Figure FDA0003533977270000016
a67＝1，a71＝eqx，a75＝eqy，
Figure FDA0003533977270000017
ex，eqx，ey，eqy，eh，eqh均为水泵水轮机线性模型待辨识参数，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，bp为永态差值系数，Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，Tr为水击相长，Ta为机组惯性时间常数；
所述抽水蓄能机组调速***线性模型的***输出表达式为
Figure FDA0003533977270000021
其中
Figure FDA0003533977270000022
ex、eh均为水泵水轮机线性模型待辨识参数，Ta为机组惯性时间常数，ey为力矩对导叶开度相对系数；
对于抽水蓄能机组调速***线性模型的***输出表达式，当
Figure FDA0003533977270000023
中***的参数选定时，抽水蓄能机组调速***能够作为Y与U的单变量***，一组模型输入变量U会对应一个***响应输出向量Y，提取模型的输入变量U，则模型表示为
Figure FDA0003533977270000024
将
Figure FDA0003533977270000025
定义为源模型；
基于映射理论象与原象直接的对应关系，当***输入U不变时，定义由Y到X的变维映射过程为映射模型，其表示式为
Figure FDA0003533977270000026
即映射模型的输出为抽水蓄能机组调速***待辨识参数；所述步骤b中的确定抽水蓄能机组调速***待辨识参数的上下边界，对于水泵水轮机线性模型待辨识参数[ex，eqx，ey，eqy，eh，eqh]，由机组运行实测数据得其对应的机组运行工作水头Hw，当Hw为定值时，机组的流量和力矩变化的因变量则为导叶开度y，保持机组单位转速为额定单位转速n11r，将导叶开度在零开度到额定开度区间内，以5％的偏移量计算20组水泵水轮机线性***模型六系数，从中选取相应的最大值和最小值作为抽水蓄能机组调速***待辨识参数的上下边界值；所述步骤b中的获取参数辨识样本数据为，先基于建立的***线性模型，保持线性模型中的输入信号U不变，以相同的偏移量增加或减小X中的待辨识参数，并将待辨识参数向量xi代入线性模型进行仿真，获得模型输出Yi；其次，重复上述操作，直到获得参数辨识所需的样本数量；接下来，为了增强神经网络的学习能力，计算调速***输出Yi的4个指标量包括：上升时间、调节时间、超调量和波峰次数，将指标量与Yi构成新的Yi，并将xi与Yi进行组合，得到BP神经网络辨识所需的样本集{(x1，Y1)，(x2，Y2)，...，(xi，Yi)，...，(xN，YN)}，N为样本集中的样本数；最后，从获得的样本集中随机抽取部分样本作为训练样本，其余为测试样本。
2.根据权利要求1所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤a中的PID调速器模型包括并联PID控制器和辅助接力器型执行机构；所述并联PID控制器传递函数表达式为
Figure FDA0003533977270000031
其中s为拉普拉斯算子，Kp为比例系数，Ki为积分系数，Kd为系数，Td为微分时间常数，xc为机组转速给定值，x为机组转速，yc为接力器开度给定值，bp为永态差值系数，σ为调速器电气控制输出；
所述辅助接力器型执行机构传递函数表达式为
Figure FDA0003533977270000032
其中Ty为主接力器响应时间常数，TyB为辅助接力器响应时间常数，k0为放大系数，y为执行机构输出，s为拉普拉斯算子，σ为调速器电气控制输出。
3.根据权利要求1所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤a中的水泵水轮机IEEE六参数模型采用相对值表达式为
Figure FDA0003533977270000033
和
Figure FDA0003533977270000034
其中mt为力矩偏差相对值，q为流量偏差相对值，y为导叶开度偏差相对值，x为转速偏差相对值，h为水头偏差相对值，ey为力矩对导叶开度相对系数，ex为力矩对转速传递系数，eh为力矩对工作水头传递系数，eqy为流量对导叶开度相对系数，eqx为流量对转速相对系数，eqh为流量对工作水头相对系数，f为摩阻系数。
4.根据权利要求1所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤a中的引水***近似弹性水击模型的传递函数表达式为
Figure FDA0003533977270000035
其中hw为管道特性系数，Tr为水击相长，f为摩阻系数，s为拉普拉斯算子。
5.根据权利要求1所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤a中的发电电动机***模型采用转子运动方程描述发电电动机动态特性，发电电动机***模型的传递函数表达式为
Figure FDA0003533977270000041
其中Ta为机组惯性时间常数，eg为发电机自调节系数，s为拉普拉斯算子。
6.根据权利要求1所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤c中基于训练样本为，采用BP神经网络方法训练逼近抽水蓄能机组线性模型的映射模型，并用测试样本集验证逼近的精度，待精度符合辨识要求时，即建立了基于调速***线性模型映射的BP神经网络参数辨识模型。
7.根据权利要求1至6中任一权利要求所述的一种水泵水轮机线性模型参数辨识方法，其特征在于：所述步骤d中以抽水蓄能机组调速***在空载开机工况的原始运行实测数据作为步骤c中建立的BP神经网络参数辨识模型的输入，即得到一组参数向量x，参数向量x为***的参数辨识结果。
CN201810662132.9A 2018-06-25 2018-06-25 一种水泵水轮机线性模型参数辨识方法 Active CN108763831B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810662132.9A CN108763831B (zh) 2018-06-25 2018-06-25 一种水泵水轮机线性模型参数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810662132.9A CN108763831B (zh) 2018-06-25 2018-06-25 一种水泵水轮机线性模型参数辨识方法

Publications (2)

Publication Number Publication Date
CN108763831A CN108763831A (zh) 2018-11-06
CN108763831B true CN108763831B (zh) 2022-05-31

Family

ID=63976862

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810662132.9A Active CN108763831B (zh) 2018-06-25 2018-06-25 一种水泵水轮机线性模型参数辨识方法

Country Status (1)

Country Link
CN (1) CN108763831B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109654266B (zh) * 2018-11-15 2020-01-21 西安建筑科技大学 一种基于多相流水击模型的泄压阀压力精度设定方法
CN109840382B (zh) * 2019-02-15 2023-10-24 东莞维科电池有限公司 一种挤压涂布机泵速计算方法、装置和存储介质
CN110361974B (zh) * 2019-07-16 2022-10-04 华北水利水电大学 基于bp-foa混合算法的水轮机调速***优化方法
CN110889218B (zh) * 2019-11-20 2023-09-01 天生桥二级水力发电有限公司天生桥水力发电总厂 基于神经网络的水轮机非线性建模方法
CN112270088B (zh) * 2020-10-26 2023-07-18 中国长江电力股份有限公司 集群水轮发电机组仿真模型及其建模方法
CN112507480B (zh) * 2020-11-25 2024-06-07 浙江同善人工智能技术有限公司 一种惯性参数辨识方法
CN113642133B (zh) * 2021-09-01 2023-05-09 中国电建集团成都勘测设计研究院有限公司 水斗式水轮机转轮内宽计算方法
CN113987802B (zh) * 2021-10-28 2023-04-11 华中科技大学 一种水泵水轮机特性系数的计算方法及***
CN114352469B (zh) * 2021-12-08 2024-02-09 国网江西省电力有限公司电力科学研究院 一种高阶水轮机调节***稳定域计算方法及***
CN114677532B (zh) * 2022-04-20 2024-05-28 重庆邮电大学 一种利用出厂数据和特性曲线的感应电动机参数辨识方法
CN114997073B (zh) * 2022-07-29 2022-12-09 浙江大学 冲击式水轮机喷嘴结构参数优化方法及***
CN117471905B (zh) * 2023-09-21 2024-04-12 昆明理工大学 一种改进樽海鞘群的抽水蓄能机组调节***的控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104281130A (zh) * 2014-09-22 2015-01-14 国家电网公司 一种基于大数据技术的水电设备监测和故障诊断***
CN105956350A (zh) * 2016-07-05 2016-09-21 华中科技大学 一种抽水蓄能机组过水***的建模方法
CN106125552A (zh) * 2016-08-08 2016-11-16 国家电网公司 抽水蓄能机组调速***模糊分数阶pid控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7454270B2 (en) * 2003-05-13 2008-11-18 Siemens Power Transmission & Distribution, Inc. Dynamic economic dispatch for the management of a power distribution system
CN103529698B (zh) * 2013-10-17 2016-09-14 广东电网公司电力科学研究院 发电机调速***参数辨识方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104281130A (zh) * 2014-09-22 2015-01-14 国家电网公司 一种基于大数据技术的水电设备监测和故障诊断***
CN105956350A (zh) * 2016-07-05 2016-09-21 华中科技大学 一种抽水蓄能机组过水***的建模方法
CN106125552A (zh) * 2016-08-08 2016-11-16 国家电网公司 抽水蓄能机组调速***模糊分数阶pid控制方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A real-time accurate model and its predictive fuzzy PID controller for pumped storage unit via error compensation;Jianzhong Zhou等;《energies》;20171224;1-24 *
An adaptively fast fuzzy fractional order PID control for pumped storage hydro unit using improved gravitational search algorithm;YanheXu等;《Energy Conversion and Management》;20160301;第111卷;67-78 *
Real Time Implementation of Grid Connected Wind Energy Conversion System: Fractional Order Controller;A. Beddar等;《International Conference on Technological Advances in Electrical Engineering (ICTAEE’16.)》;20161031;1-8 *
基于ANN的抽水蓄能电站建模与过渡过程优化;曾洪涛等;《排灌机械工程学报》;20141031;第32卷(第10期);864-870+876 *
抽水蓄能机组原动机及调速***参数实测、建模与仿真;彭天波等;《第十八次中国水电设备学术讨论会论文集》;20161027;450-457 *

Also Published As

Publication number Publication date
CN108763831A (zh) 2018-11-06

Similar Documents

Publication Publication Date Title
CN108763831B (zh) 一种水泵水轮机线性模型参数辨识方法
CN107066712B (zh) 基于导叶开度—功率特性的水轮机模型参数辨识方法及装置
CN103926875B (zh) 一种滚珠丝杠进给***摩擦补偿方法
CN103529698B (zh) 发电机调速***参数辨识方法
CN110206685B (zh) 风电场中的风力发电机组的前馈控制方法和设备
Saarinen et al. Field measurements and system identification of three frequency controlling hydropower plants
CN107168103B (zh) 水电机组调速控制仿真方法和***
CN103219725A (zh) 一种基于实时运行数据的风电场等值建模方法
CN109787251A (zh) 一种集群温控负荷聚合模型、***参数辨识及反推控制方法
CN106452247A (zh) 一种永磁同步电机转动惯量辨识的方法及装置
CN104038128B (zh) 基于nurbs的直线电机推力波动控制方法
Jasniewicz et al. Wind turbine modelling and identification for control system applications
GB2555010B (en) Determining loads on a wind turbine
CN105024612A (zh) 一种基于参数辨识的电机电流控制方法及***
CN106059412B (zh) 基于信度规则库推理的他励直流电动机转速控制方法
CN111123698A (zh) 一种水力发电机组调节***的无模型自适应pid控制方法
Burian et al. Investigation of the pump unit control system with the neural network productivity estimator
CN109657380A (zh) 一种基于扩展卡尔曼滤波的双馈风电场动态等值方法
CN116255493B (zh) 一种伺服阀的调节控制方法及***
CN109254569A (zh) 一种面向快速频率响应下的水电机组动态模型建立方法
CN103532463A (zh) 基于状态观测器的直流马达转速控制方法及***
CN104123442B (zh) 一种风电机组有功功率在线等值模型建模方法
CN111637008B (zh) 共享机制下风电场风力发电机变桨***的控制方法及***
CN112270088B (zh) 集群水轮发电机组仿真模型及其建模方法
Olexa et al. Simulation and Laboratory Model of Small Hydropower Plant

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