CN110729936B - 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法 - Google Patents

一种计及饱和特性的发电机实用参数辨识模型及其辨识方法 Download PDF

Info

Publication number
CN110729936B
CN110729936B CN201910679237.XA CN201910679237A CN110729936B CN 110729936 B CN110729936 B CN 110729936B CN 201910679237 A CN201910679237 A CN 201910679237A CN 110729936 B CN110729936 B CN 110729936B
Authority
CN
China
Prior art keywords
generator
axis
transient
formula
parameter identification
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
CN201910679237.XA
Other languages
English (en)
Other versions
CN110729936A (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.)
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
State Grid Fujian Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
State Grid Fujian Electric Power 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 Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd, State Grid Fujian Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Fujian Electric Power Co Ltd
Priority to CN201910679237.XA priority Critical patent/CN110729936B/zh
Publication of CN110729936A publication Critical patent/CN110729936A/zh
Application granted granted Critical
Publication of CN110729936B publication Critical patent/CN110729936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • 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
    • 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/22—Current control, e.g. using a current control loop

Abstract

本发明涉及一种计及饱和特性的发电机实用参数辨识模型及其辨识方法，首先，运用Laplace变换，把发电机六阶实用时域模型方程组转换为频域模型方程组，对频域方程组进行代数变换，消去不可观测量e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0。然后对其作Laplace反变换，转换为以积分形式表示的消去不可观测量的计及饱和特性的实用参数辨识模型。接着对积分形式的参数辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数。在此基础上，建立了适用于PMU量测数据的发电机参数辨识策略，从而有效提高发电机参数辨识的精度。

Description

一种计及饱和特性的发电机实用参数辨识模型及其辨识方法
技术领域
本发明涉及发电机实用参数辨识精度领域，特别是一种计及饱和特性的发电机实用参数辨识模型及其辨识方法。
背景技术
我国正加快推进电力系统精细化建模工作，其中励磁调速系统及PSS参数入库工作正逐年完善，但同步发电机参数仍采用铭牌参数进行仿真计算。发电机模型参数的准确性是导致电力系统时域仿真结果与实际情况偏差的重要因素，将直接影响到系统的隐患排查、事故复现、安全稳定和规划设计。近年来，利用同步相量量测装置PMU(PhasorMeasurement Unit)实现电力系统模型参数在线辨识成为重要的研究方向，其具备测量发电机实时状态变量的特点，可实现发电机自身与电网、励磁和调速系统解耦，有利于优化辨识策略，且量测到的系统扰动数据计及发电机真实运行工况，具有较高的可信度。
在实用模型不可观测量及导数项的处理方面，国内学者提出了适用于甩负荷试验的时频变换处理方法，将实用模型的微分方程组变换为积分方程，并消去中间变量，构建仅含有可观测电气量与待辨识参数的方程，避免中间变量反复迭代的误差。其中，提出的励磁感应电动势ef修正方法，计及了发电机饱和特性对参数辨识精度的影响。但该方法需使用励磁电流和电枢电流合成值If∑，由于扰动时电枢电流突变，矢量叠加合成的If∑实时值在突变瞬间存在误差，影响辨识精度。
发明内容
有鉴于此，本发明的目的是提出一种计及饱和特性的发电机实用参数辨识模型及其辨识方法，用参数辨识模型为dq解耦方程，降低参数辨识的维度，且避免对不可观测状态变量反复迭代的计算过程，提高辨识的快速性。
本发明采用以下方案实现：一种计及饱和特性的发电机实用参数辨识模型，具体为：
Figure BDA0002144044690000021
Figure BDA0002144044690000022
式中，Xq,X′q,X″q,T′q0,T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xd,X′d,X″d,T′d0,T″d0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq,X′q,X″q,T′q0,T″q0和所述Xd,X′d,X″d,T′d0,T″d0均为辨识参数；Xl为定子漏抗；id,iq分别为发电机机端电流在d、q轴分量；ef为励磁感应电动势；令
Figure BDA0002144044690000031
式中，ud0,uq0,id0,iq0为扰动发生前发电机机端稳态电压电流在d、q轴分量；令
Figure BDA0002144044690000032
k3＝1+Sd，
Figure BDA0002144044690000033
Sd,Sq分别为d、q轴饱和因数。
进一步地，本发明还提供一种基于计及饱和特性的发电机实用参数辨识模型的建立方法，包括以下步骤：
步骤S1：利用Laplace变换，将发电机六阶实用时域模型方程组转换为频域模型方程组，并对频域方程组进行代数变换，消去不可观测量e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0；
步骤S2：对步骤S1中消去不可观测量得到后的方程组进行Laplace反变换，转换为以积分形式表示的消去不可观测量的计及饱和特性的发电机实用参数辨识模型。
进一步地，所述步骤S1具体包括以下内容：
忽略定子绕组的电磁暂态过程和电枢电阻Ra，计及励磁绕组f、转子d、q轴的瞬变过程和发电机饱和特性，考虑饱和特性影响的同步发电机六阶实用模型为：
Figure BDA0002144044690000041
式中，Xq,X′q,X″q,T′q0,T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xd,X′d,X″d,T′d0,T″d0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq,X′q,X″q,T′q0,T″q0和所述Xd,X′d,X″d,T″d0,T″d0均为辨识参数；Xl为定子漏抗；ud,id,uq,iq分别为发电机机端电压电流在d、q轴分量；ef为励磁感应电动势；e′q,e″q,e′d,e″d分别为为发电机交、直轴暂态及次暂态电动势，均为状态变量；Sd,Sq分别为d、q轴饱和因数；δ,ω,TJ,Tm,Te,D依次为发电机功角、角速度、转动惯量、机械转矩、电磁转矩和阻尼系数；p为微分算子；
令式(3)中p＝0，得发电机稳态方程
Figure BDA0002144044690000042
令k1＝X′d-X″d,
Figure BDA0002144044690000043
k4＝X′q-X″q,
Figure BDA0002144044690000044
k6＝1+Sq，并将式(3)作Laplace变换，得
Figure BDA0002144044690000051
式中，e′q0,e″q0,e′d0,e″d0分别为发电机交直轴暂态和次暂态电动势初始值；
将式(5)方程组中第3式和第5式转换为e″q和e″d的表达式，为
Figure BDA0002144044690000052
将式(6)代入式(5)方程组中的第4式和第6式，得
Figure BDA0002144044690000053
将式(5)方程组中式1和式2转换为e″q和e″d的表达式，为
Figure BDA0002144044690000061
同理，给出不计阻尼回路时的发电机暂态电动势e′q和e′d的表达式，为
Figure BDA0002144044690000062
同理，将发电机交直轴暂态和次暂态电动势初始值e′q0,e″q0,e′d0,e″d0也用式(8)和式(9)表示；将e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0表达式一并代入式(7)。
进一步地，所述步骤S2具体包括以下内容：
令
Figure BDA0002144044690000063
代入式(7)，最后对公式(7)作Laplace反变换，消去不可观测量e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0，得
Figure BDA0002144044690000064
Figure BDA0002144044690000071
式中，ud0,uq0,id0,iq0分别为扰动发生前发电机机端稳态电压电流在d、q轴分量；式(10)～式(11)便是计及饱和特性的发电机实用参数辨识模型。
进一步地，本发明还提供一种计及饱和特性的发电机实用参数辨识模型的辨识方法，首先，将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；然后，在优化函数的基础上所述辨识方法包括以下步骤：
步骤SA：对同步相量量测装置PMU量测数据进行标幺化；
步骤SB：将步骤SA中处理后的量测数据代入式(4)求解稳态参数；
步骤SC：提供厂家的发电机铭牌参数作为辨识参数初值；
步骤SD：将稳态参数Xd,Xq代入式(12)～式(13)，采用遗传算法求解得到dq轴的暂态及次暂态参数；
步骤SE：若发电机参数辨识模型输出的dq轴电流值和实际电流值差值满足误差2.5％要求，则迭代结束，输出本次发电机dq轴的暂态及次暂态参数；若不满足要求，则进行下一次迭代，重复步骤SD。
进一步的，将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；具体包括以下内容：
对式(1)～式(2)离散化，以dq轴电流的拟合值和实测值差值最小作为寻优目标值构建函数，并设定辨识的边界条件，表示为：
Figure BDA0002144044690000081
Figure BDA0002144044690000082
Figure BDA0002144044690000083
Figure BDA0002144044690000084
其中，tn表示第n个采样时刻；N为总采样点数；含有下标c的参数为厂家给定参数值；边界条件限制辨识参数X′q,X″q,T′q0,T″q0,X′d,X″d,T′d0,T″d0在厂家给定参数的0.7倍到1.3倍范围内。
与现有技术相比，本发明有以下有益效果：
(1)本发明的实用参数辨识模型为dq解耦方程，降低参数辨识的维度，且避免对不可观测状态变量反复迭代的计算过程，提高辨识的快速性。
(2)本发明积分方程形式避免大扰动过程局部线性化问题，提高辨识的稳定性。
(3)本发明的模型是对发电机六阶实用饱和模型的直接推导，推导过程未增设新的假设，同时计及了发电机饱和特性，且模型避免了不可观测状态变量反复迭代引起的误差，提高辨识的精度。
附图说明
图1为本发明实施例的同步发电机参数辨识流程图。
图2为本发明实施例的单相接地故障下dq轴仿真与拟合电流对比图，其中图2(a)d轴仿真电流与拟合电流对比图，为图2(b)为q轴仿真电流与拟合电流对比图。
图3为本发明实施例的三相短路故障下dq轴仿真与拟合电流对比图，其中图3(a)为d轴仿真电流与拟合电流对比图，图3(b)为q轴仿真电流与拟合电流对比图。
图4为本发明实施例的功率波动下dq轴实测与辨识值、铭牌值计算电流对比图，其中图4(a)为d轴实测电流与辨识值、铭牌值计算电流对比图，图4(b)为q轴实测与辨识值、铭牌值计算电流对比图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
本实施例提供了一种计及饱和特性的发电机实用参数辨识模型，具体为：
Figure BDA0002144044690000101
Figure BDA0002144044690000102
式中，Xq,X′q,X″q,T′q0,T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xd,X′d,X″d,T′d0,T″d0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq,X′q,X″q,T′q0,T″q0和所述Xd,X′d,X″d,T′d0,T″d0均为辨识参数；Xl为定子漏抗；id,iq分别为发电机机端电流在d、q轴分量；ef为励磁感应电动势；令
Figure BDA0002144044690000111
式中，ud0,uq0,id0,iq0为扰动发生前发电机机端稳态电压电流在d、q轴分量；令
Figure BDA0002144044690000112
k3＝1+Sd，
Figure BDA0002144044690000113
Sd,Sq分别为d、q轴饱和因数。
较佳的，本实施例还提供一种基于计及饱和特性的发电机实用参数辨识模型的建立方法，包括以下步骤：
步骤S1：利用Laplace变换，将发电机六阶实用时域模型方程组转换为频域模型方程组，并对频域方程组进行代数变换，消去不可观测量e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0；
步骤S2：对步骤S1中消去不可观测量得到后的方程组进行Laplace反变换，转换为以积分形式表示的消去不可观测量的计及饱和特性的发电机实用参数辨识模型。
在本实施例中，所述步骤S1具体包括以下内容：
忽略定子绕组的电磁暂态过程和电枢电阻Ra，计及励磁绕组f、转子d、q轴的瞬变过程和发电机饱和特性，考虑饱和特性影响的同步发电机六阶实用模型为：
Figure BDA0002144044690000121
式中，Xq,X′q,X″q,T′q0,T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xd,X′d,X″d,T′d0,T″d0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq,X′q,X″q,T′q0,T″q0和所述Xd,X′d,X″d,T′d0,T″d0均为辨识参数；Xl为定子漏抗；ud,id,uq,iq分别为发电机机端电压电流在d、q轴分量；ef为励磁感应电动势；e′q,e″q,e′d,e″d分别为为发电机交、直轴暂态及次暂态电动势，均为状态变量；Sd,Sq分别为d、q轴饱和因数；δ,ω,TJ,Tm,Te,D依次为发电机功角、角速度、转动惯量、机械转矩、电磁转矩和阻尼系数；p为微分算子；
令式(3)中p＝0，得发电机稳态方程
Figure BDA0002144044690000122
令k1＝X′d-X″d,
Figure BDA0002144044690000123
k4＝X′q-X″q,
Figure BDA0002144044690000124
k6＝1+Sq，并将式(3)作Laplace变换，得
Figure BDA0002144044690000131
式中，e′q0,e″q0,e′d0,e″d0分别为发电机交直轴暂态和次暂态电动势初始值；
将式(5)方程组中第3式和第5式转换为e″q和e″d的表达式，为
Figure BDA0002144044690000132
将式(6)代入式(5)方程组中的第4式和第6式，得
Figure BDA0002144044690000133
将式(5)方程组中式1和式2转换为e″q和e″d的表达式，为
Figure BDA0002144044690000141
同理，给出不计阻尼回路时的发电机暂态电动势e′q和e′d的表达式，为
Figure BDA0002144044690000142
同理，将发电机交直轴暂态和次暂态电动势初始值e′q0,e″q0,e′d0,e″d0也用式(8)和式(9)表示；将e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0表达式一并代入式(7)。
在本实施例中，所述步骤S2具体包括以下内容：
令
Figure BDA0002144044690000143
并代入式(7)，最后对公式(7)作Laplace反变换，消去不可观测量e′q,e″q,e′q0,e″q0,e′d,e″d,e′d0,e″d0，得
Figure BDA0002144044690000144
Figure BDA0002144044690000151
式中，ud0,uq0,id0,iq0分别为扰动发生前发电机机端稳态电压电流在d、q轴分量；式(10)～式(11)便是计及饱和特性的发电机实用参数辨识模型。可以看出，以积分方程形式表示的发电机实用参数辨识模型仅包含发电机dq轴电气量和待辨识参数，dq轴电气量由PMU量测数据计算得到，该模型适合用于在线辨识。
在本实施例中，将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；具体内容为：
将发电机参数辨识转化为非线性优化问题，对式(1)～式(2)离散化，以dq轴电流的拟合值和实测值差值最小作为寻优目标值构建函数，并设定辨识的边界条件，表示为：
Figure BDA0002144044690000152
Figure BDA0002144044690000153
Figure BDA0002144044690000161
Figure BDA0002144044690000162
其中，式(12)-式(13)即为优化函数，tn表示第n个采样时刻；N为总采样点数；含有下标c的参数为厂家给定参数值；边界条件限制辨识参数X′q,X″q,T′q0,T″q0,X′d,X″d,T′d0,T″d0在厂家给定参数的0.7倍到1.3倍范围内，足以满足工程应用需要。
由于同步发电机的dq轴自然解耦特性，在参数辨识过程中可以将d、q轴参数分开辨识，减少了单次辨识的优化变量个数，进一步采用本实施例所述的计及饱和特性的发电机实用参数辨识模型，可以将同步发电机稳态参数与暂态/次暂态参数分开辨识，使得参数辨识过程中每次辨识优化的变量降低到4个。采用dq轴参数解耦、暂稳态参数分步的辨识策略，可以显著减小参数辨识过程中的优化维度，有利于提高参数辨识精度，同步发电机参数辨识流程如图1所示。
综合前述优化辨识模型，形成基于PMU量测数据的发电机参数辨识策略。
较佳的，本实施例还提供一种基于计及饱和特性的发电机实用参数辨识模型的辨识方法，首先，将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；然后，在优化函数的基础上所述辨识方法包括以下步骤：
步骤SA：对同步相量量测装置PMU量测数据进行标幺化预处理；
步骤SB：将步骤SA中处理后的量测数据代入式(4)求解稳态参数；
步骤SC：提供厂家的发电机铭牌参数作为辨识参数初值；
步骤SD：将稳态参数Xd,Xq代入式(12)～式(13)，采用遗传算法求解得到dq轴的暂态及次暂态参数；
步骤SE：若发电机参数辨识模型输出的dq轴电流值和实际电流值差值满足误差要求，则迭代结束，输出本次发电机dq轴的暂态及次暂态参数；若不满足要求，则进行下一次迭代，重复步骤SD。
特别的，本实施例还提供了最佳实施方式，进行算理验证。算理验证分为仿真算例和实测验证。
采用PSD-BPA搭建的某电网2018年典型运行方式，在某电厂发电机组出口处设置0.5s单相接地故障和三相短路故障，参数辨识模型分别采用计及饱和特性和非饱和特性。dq轴仿真与拟合电流曲线如图2～图3所示，参数辨识结果详见表1。
由图2～图3可知，对于典型对称或非对称故障及故障恢复后的扰动数据，10ms采样间隔的发电机实用参数辨识输出的dq轴拟合电流与仿真电流曲线高度吻合。
表1基于不同仿真故障数据的参数辨识结果对比
Figure BDA0002144044690000181
由图2～图3可知，对于典型对称或非对称故障及故障恢复后的扰动数据，采用计及饱和特性的发电机实用参数辨识模型输出的dq轴拟合电流与仿真电流曲线高度吻合。由表1可见，计及饱和特性的参数辨识模型的辨识结果误差小于2.5％，非饱和特性模型的T″q0辨识结果与设定值极为不符，造成均方根误差大，说明模型准确性对于参数辨识结果影响较大。由上分析可知，计及饱和特性的发电机实用参数辨识模型对不同故障类型具有通用性，且参数辨识结果精度高。
实测验证
为进一步验证本发明参数辨识方法的有效性，采用福建南埔电厂实测PMU扰动数据验证辨识方法的有效性。此扰动发生在机组进行单阀到顺阀方式切换过程中，受到阀切换过程中调门流量非线性影响，调速系统对机组产生负阻尼，当电网频率降低引起调速系统动作时，由于一次调频过调，造成机组功率大幅度波动。此阀切换期间，伴随电网频率多次降低，机组频差数次超过死区，调速系统频繁动作，相应PMU装置录波到多组扰动数据。用本发明计算dq轴电流、铭牌值计算dq轴电流，与本次实测电流对比如图4所示。
由图4可知，相比铭牌参数，辨识值计算的电流拟合效果优于铭牌值，证明了参数辨识值的准确性，也进一步验证了计及饱和特性的发电机实用参数辨识模型及辨识方法能够有效提高辨识精度。
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。

Claims (6)

1.一种计及饱和特性的发电机实用参数辨识模型，其特征在于：具体为：
Figure FDA0003149421430000011
Figure FDA0003149421430000012
式中，Xq，X′q，X″q，T′q0，T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xq，X′q，X″q，T′q0，T″q0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq，X′q，X″q，T′q0，T″q0和所述Xq，X′q，X″q，T′q0，T″q0均为辨识参数；Xl为定子漏抗；id，iq分别为发电机机端电流在d、q轴分量；ef为励磁感应电动势；令
Figure FDA0003149421430000013
式中，ud0，uq0，id0，iq0为扰动发生前发电机机端稳态电压电流在d、q轴分量；令
Figure FDA0003149421430000014
k3=1+Sd，
Figure FDA0003149421430000015
k6=1+Sp；Sd，Sq分别为d、q轴饱和因数。
2.一种根据权利要求1所述的计及饱和特性的发电机实用参数辨识模型的建立方法，其特征在于：包括以下步骤：
步骤S1：利用Laplace变换，将发电机六阶实用时域模型方程组转换为频域模型方程组，并对频域方程组进行代数变换，消去不可观测量e′q，e″q，e′q0，e″q0，e′d，e″d，e′d0，e″d0；
步骤S2：对步骤S1中消去不可观测量得到后的方程组进行Laplace反变换，转换为以积分形式表示的消去不可观测量的计及饱和特性的发电机实用参数辨识模型。
3.根据权利要求2所述的一种计及饱和特性的发电机实用参数辨识模型的建立方法，其特征在于：所述步骤S1具体包括以下内容：
忽略定子绕组的电磁暂态过程和电枢电阻Ra；计及励磁绕组f、转子d、q轴的瞬变过程和发电机饱和特性；考虑饱和特性影响的同步发电机六阶实用模型为：
Figure FDA0003149421430000021
式中，Xq，X′q，X″q，T′q0，T″q0分别为交轴同步电抗、交轴暂态和次暂态电抗、交轴暂态和次暂态开路时间常数；Xd，X′d，X″d，T′d0，T″d0分别为直轴同步电抗、直轴暂态和次暂态电抗、直轴暂态和次暂态开路时间常数；所述Xq，X′q，X″q，T′q0，T″q0和所述Xd，X′d，X″d，T′d0，T″d0均为辨识参数；Xl为定子漏抗；ud，id，uq，iq分别为发电机机端电压电流在d、q轴分量；ef为励磁感应电动势；e′q，e″q，e′d，e″d分别为为发电机交、直轴暂态及次暂态电动势，均为状态变量；Sd，Sq分别为d、q轴饱和因数；δ，ω，TJ，Tm，Te，D依次为发电机功角、角速度、转动惯量、机械转矩、电磁转矩和阻尼系数；p为微分算子；
令式(3)中p＝0，得发电机稳态方程
Figure FDA0003149421430000031
令k1=X′d-X″d，
Figure FDA0003149421430000032
k3=1+Sd，k4=X′q-X″q，
Figure FDA0003149421430000033
k6=1+Sq，并将式(3)作Laplace变换，得
Figure FDA0003149421430000034
式中，e′q0，e″q0，e′d0，e″d0分别为发电机交直轴暂态和次暂态电动势初始值；
将式(5)方程组中第3式和第5式转换为e″q和e″d的表达式，为
Figure FDA0003149421430000035
将式(6)代入式(5)方程组中的第4式和第6式，得
Figure FDA0003149421430000041
将式(5)方程组中式1和式2转换为e″q和e″d的表达式，为
Figure FDA0003149421430000042
同理，给出不计阻尼回路时的发电机暂态电动势e′q和e′d的表达式，为
Figure FDA0003149421430000043
同理，将发电机交直轴暂态和次暂态电动势初始值e′q0，e″q0，e′d0，e″d0也用式(8)和式(9)表示；将e′q，e″q，e′q0，e″q0，e′d，e″d，e′d0，e″d0表达式一并代入式(7)。
4.根据权利要求3所述的一种计及饱和特性的发电机实用参数辨识模型的建立方法，其特征在于：所述步骤S2具体包括以下内容：
令
Figure FDA0003149421430000044
代入式(7)，最后对公式(7)作Laplace反变换，消去了不可观测量e′q，e″q，e′q0，e″q0，e′d，e″d，e′d0，e″d0，得
Figure FDA0003149421430000045
Figure FDA0003149421430000051
式中，ud0，dq0，id0，iq0分别为扰动发生前发电机机端稳态电压电流在d、q轴分量；式(10)～式(11)便是计及饱和特性的发电机实用参数辨识模型。
5.一种根据权利要求4所述的计及饱和特性的发电机实用参数辨识模型的辨识方法，其特征在于：首先，将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；然后，在优化函数的基础上所述辨识方法包括以下步骤：
步骤SA：对同步相量量测装置PMU量测数据进行标幺化；
步骤SB：将步骤SA中处理后的量测数据代入式(4)求解稳态参数；
步骤SC：提供厂家的发电机铭牌参数作为辨识参数初值；
步骤SD：将稳态参数Xd，Xq代入式优化函数中，采用遗传算法求解得到dq轴的暂态及次暂态参数；
步骤SE：若发电机参数辨识模型输出的dq轴电流值和实际电流值差值满足误差2.5％要求，则迭代结束，输出本次发电机dq轴的暂态及次暂态参数；若不满足要求，则进行下一次迭代，重复步骤SD。
6.根据权利要求5所述的一种计及饱和特性的发电机实用参数辨识模型的辨识方法，其特征在于：将所述辨识模型离散化，构造以dq轴电流计算值与实测值差值最小为目标的优化函数；具体内容为：
对式(1)～式(2)离散化，以dq轴电流的拟合值和实测值差值最小作为寻优目标值构建函数，并设定辨识的边界条件，表示为：
Figure FDA0003149421430000061
Figure FDA0003149421430000062
Figure FDA0003149421430000063
Figure FDA0003149421430000064
其中，式(12)-式(13)即为优化函数，tn表示第n个采样时刻；N为总采样点数；含有下标c的参数为厂家给定参数值；边界条件限制辨识参数X′q，X″q，T′q0，T″q0，X′d，X″d，T′d0，T″d0在厂家给定参数的0.7倍到1.3倍范围内。
CN201910679237.XA 2019-07-25 2019-07-25 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法 Active CN110729936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910679237.XA CN110729936B (zh) 2019-07-25 2019-07-25 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910679237.XA CN110729936B (zh) 2019-07-25 2019-07-25 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法

Publications (2)

Publication Number Publication Date
CN110729936A CN110729936A (zh) 2020-01-24
CN110729936B true CN110729936B (zh) 2021-11-09

Family

ID=69217096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910679237.XA Active CN110729936B (zh) 2019-07-25 2019-07-25 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法

Country Status (1)

Country Link
CN (1) CN110729936B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103036498A (zh) * 2012-10-18 2013-04-10 中国电力科学研究院 一种基于pmu的同步发电机实用模型参数校核与辨识方法
CN107677960A (zh) * 2017-09-21 2018-02-09 国网福建省电力有限公司 基于pmu量测的同步发电机参数辨识方法
CN109274892A (zh) * 2018-11-08 2019-01-25 江苏方天电力技术有限公司 一种计及饱和效应的调相机参数分步辨识方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103036498A (zh) * 2012-10-18 2013-04-10 中国电力科学研究院 一种基于pmu的同步发电机实用模型参数校核与辨识方法
CN107677960A (zh) * 2017-09-21 2018-02-09 国网福建省电力有限公司 基于pmu量测的同步发电机参数辨识方法
CN109274892A (zh) * 2018-11-08 2019-01-25 江苏方天电力技术有限公司 一种计及饱和效应的调相机参数分步辨识方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
六阶模型同步发电机参数辨识的矩阵运算方法;陈剑等;《大电机技术》;20140731(第7期);全文 *

Also Published As

Publication number Publication date
CN110729936A (zh) 2020-01-24

Similar Documents

Publication Publication Date Title
Gagnon et al. Modeling and real-time simulation of a doubly-fed induction generator driven by a wind turbine
CN107783006B (zh) 水轮发电机转子绕组匝间短路故障检测方法
Zaker et al. Improving synchronous generator parameters estimation using $ d-q $ axes tests and considering saturation effect
CN111555312B (zh) 一种适用于电力系统超低频振荡稳定性评估的方法
CN103825521B (zh) 一种双馈风电机组中驱动系统及发电机参数的辨识方法
CN104852652B (zh) 同步风力发电机闭环矢量控制方法和系统
CN104579083A (zh) 永磁同步电机矢量控制方法及系统
CN103944481B (zh) 一种交流异步电机矢量控制系统模型参数在线修正方法
CN102510263A (zh) 基于抛载试验和数值差分的同步发电机实用参数辨识方法
CN102611380B (zh) 一种双馈电机参数在线辨识方法
CN102073012A (zh) 同步发电机参数获取方法及系统
CN103023418A (zh) 一种基于广域量测信息的同步发电机在线参数辨识方法
CN116054190A (zh) 一种双馈风电场控制参数对次同步振荡影响的分析方法
CN102983807A (zh) 异步电机转子时间常数在线识别系统及方法
CN110601175B (zh) 一种励磁电流改进标幺方法
CN110729936B (zh) 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法
CN106992732A (zh) 一种电机磁共能模型建立系统
CN105552951B (zh) 一种基于重复滑模的dfig系统控制方法
CN116722533A (zh) 一种基于随机风电预测区间的参数化暂态稳定预防控制求解方法
CN108462207B (zh) 一种计及变流器暂态调控的双馈电机短路电流计算方法
Lešić et al. Fault-tolerant control of a wind turbine with generator stator inter-turn faults
CN113783183B (zh) 弱电网下双馈风机在故障穿越期间的暂态稳定性评估方法
Hanic et al. Steady-state synchronous machine model which incorporates saturation and cross-magnetization effects
CN104809268A (zh) 一种汽轮发电机负载励磁电流计算方法
Luo et al. Synchronous generator modeling and semi-physical simulation

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