CN115117879A - Power plant parameter identification method and device and computer readable storage medium - Google Patents
Power plant parameter identification method and device and computer readable storage medium Download PDFInfo
- Publication number
- CN115117879A CN115117879A CN202210764472.9A CN202210764472A CN115117879A CN 115117879 A CN115117879 A CN 115117879A CN 202210764472 A CN202210764472 A CN 202210764472A CN 115117879 A CN115117879 A CN 115117879A
- Authority
- CN
- China
- Prior art keywords
- power
- voltage
- value
- current
- generator
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000005259 measurement Methods 0.000 claims abstract description 64
- 230000005284 excitation Effects 0.000 claims abstract description 55
- 230000006870 function Effects 0.000 claims description 33
- 230000001052 transient effect Effects 0.000 claims description 28
- 230000033228 biological regulation Effects 0.000 claims description 25
- 239000003381 stabilizer Substances 0.000 claims description 17
- 230000001360 synchronised effect Effects 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000006641 stabilisation Effects 0.000 claims 1
- 238000011105 stabilization Methods 0.000 claims 1
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0073—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source when the main path fails, e.g. transformers, busbars
-
- 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/18—Estimation of position or speed
-
- 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
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
- Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种电厂参数识别方法、装置及计算机可读存储介质。其中,该方法包括:基于电源管理单元PMU获取在预设的多个类噪声扰动场景下的目标电网的电压测量值、电流测量值和功角测量值;求解根据最小二乘法构建的目标电网的电压、电流和功角的测量值与模型预测值之间的误差函数,获取与最小误差值所对应的目标电压预测值、目标电流预测值和目标功角预测值;根据目标电压预测值、目标电流预测值和目标功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别电厂参数。本发明解决了相关技术中,存在无法对电厂参数进行准确识别的技术问题。
The invention discloses a power plant parameter identification method, a device and a computer-readable storage medium. Wherein, the method includes: acquiring, based on the power management unit PMU, the voltage measurement value, current measurement value and power angle measurement value of the target power grid under multiple preset noise-like disturbance scenarios; The error function between the measured value of voltage, current and power angle and the predicted value of the model, obtain the predicted value of target voltage, predicted value of target current and predicted value of target power angle corresponding to the minimum error value; Current predicted value and target power angle predicted value, as well as pre-built third-order generator model, speed control system model, excitation system model, identify power plant parameters. The invention solves the technical problem that the parameters of the power plant cannot be accurately identified in the related art.
Description
技术领域technical field
本发明涉及电力控制领域,具体而言,涉及一种电厂参数识别方法、装置及计算机可读存储介质。The present invention relates to the field of power control, and in particular, to a method and device for identifying parameters of a power plant, and a computer-readable storage medium.
背景技术Background technique
电厂关键参数辨识对电网合环操作、经济运行等具有重要意义。The identification of the key parameters of the power plant is of great significance to the loop closing operation and economic operation of the power grid.
相关技术中,通过傅里叶变化、拉普拉斯变化等技术将电网相关的物理量转换为频域量进行分析,基于这种转换对电厂参数进行辨识,但这种方法只能进行离线辨识。在相关技术中,存在无法对电厂参数进行准确识别的问题。In the related art, the physical quantities related to the power grid are converted into frequency domain quantities for analysis by techniques such as Fourier transformation and Laplace transformation, and power plant parameters are identified based on this transformation, but this method can only be used for offline identification. In the related art, there is a problem that the parameters of the power plant cannot be accurately identified.
针对上述的问题,目前尚未提出有效的解决方案。For the above problems, no effective solution has been proposed yet.
发明内容SUMMARY OF THE INVENTION
本发明实施例提供了一种电厂参数识别方法、装置及计算机可读存储介质,以至少解决相关技术中,存在无法对电厂参数进行准确识别的技术问题。Embodiments of the present invention provide a power plant parameter identification method, device, and computer-readable storage medium, so as to at least solve the technical problem that the power plant parameters cannot be accurately identified in the related art.
根据本发明实施例的一个方面,提供了一种电厂参数识别方法,包括:基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下电压测量值、电流测量值和功角测量值;根据所述目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值,求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,其中,所述误差函数是根据最小二乘法构建得到的,所述误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差,所述电力参数包括电压、电流和功角;根据所述与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别分别与所述发电机三阶模型、调速系统模型、励磁系统模型对应的电厂参数;其中,所述电厂参数为所述发电机三阶模型、调速系统模型、励磁系统模型中的参数。According to an aspect of the embodiments of the present invention, a method for identifying parameters of a power plant is provided, including: based on a power management unit PMU, respectively acquiring a voltage measurement value, a current measurement value, and a power measurement value of a target power grid under multiple preset noise-like disturbance scenarios. angle measurement value; according to multiple voltage measurement values, multiple current measurement values and multiple power angle measurement values of the target power grid under multiple preset noise-like disturbance scenarios, the error function is solved to obtain the difference between the minimum error value and the The corresponding voltage predicted value, current predicted value and power angle predicted value, wherein the error function is constructed according to the least squares method, and the error function is used to characterize the measured values of multiple power parameters of the power grid and the corresponding power parameters Error between predicted values, the power parameters include voltage, current and power angle; according to the voltage predicted value, current predicted value and power angle predicted value corresponding to the minimum error value, and the pre-built generator third-order model, speed regulation system model, and excitation system model, and identify the power plant parameters corresponding to the generator third-order model, speed regulation system model, and excitation system model respectively; wherein, the power plant parameters are the generator third-order model, Parameters in the speed control system model and the excitation system model.
可选地,还包括:根据所述目标电网中发电机的d轴同步电抗、d轴暂态电抗、惯性常数、d轴暂态开路时间常数,以及所述发电机的功角、加速度、q轴暂态电抗、电流实部和虚部,以及所述发电机有功功率、励磁电势,以及在噪声扰动场景下节点发电机的交轴电势、有功功率、直轴电流和交轴电流,构建所述发电机三阶模型;其中,与所述发电机三阶模型对应的电厂参数包括以下至少之一:所述发电机的d轴同步电抗、d轴暂态电抗、惯性常数、d轴暂态开路时间常数。Optionally, it also includes: according to the d-axis synchronous reactance, d-axis transient reactance, inertia constant, d-axis transient open-circuit time constant of the generator in the target power grid, and the power angle, acceleration, q Shaft transient reactance, current real part and imaginary part, as well as the generator active power, excitation potential, and quadrature axis potential, active power, direct axis current and quadrature axis current of the node generator under the noise disturbance scenario, construct the The generator third-order model; wherein, the power plant parameters corresponding to the generator third-order model include at least one of the following: d-axis synchronous reactance, d-axis transient reactance, inertia constant, d-axis transient reactance of the generator Open circuit time constant.
可选地,还包括:根据所述目标电网中调速系统的加速时间常数、所述调速系统中机组的惯性时间常数和缓冲时间常数、所述调速系统的参考角频率和运行角频率、所述调速系统的转子回路电阻、机械转矩,以及所述调速系统的机械转矩初始值,构建所述调速系统模型;其中,与所述调速系统模型所对应的电厂参数包括以下至少之一:所述调速系统中机组的惯性时间常数、所述调速系统中机组的缓冲时间常数、所述调速系统的转子回路电阻。Optionally, it also includes: according to the acceleration time constant of the speed regulation system in the target power grid, the inertia time constant and buffer time constant of the unit in the speed regulation system, the reference angular frequency and the operating angular frequency of the speed regulation system , the rotor loop resistance and mechanical torque of the speed control system, and the initial value of the mechanical torque of the speed control system, to construct the speed control system model; wherein, the power plant parameters corresponding to the speed control system model It includes at least one of the following: the inertia time constant of the unit in the speed control system, the buffer time constant of the unit in the speed control system, and the rotor loop resistance of the speed control system.
可选地,还包括:根据所述目标电网中电压调节器的实部电压、虚部电压、时间常数、增益倍数,以及励磁机、电压调节器和稳定器的输出电压,以及所述励磁机的增益倍数和延迟时间常数,以及励磁系统的延迟时间常数,以及电压稳定器的增益、衰减时间常数和延迟时间常数,构建所述励磁系统模型;与所述励磁系统模型所对应的电厂参数包括以下至少之一:所述电压调节器的时间常数、所述电压调节器的增益倍数、所述电压调节器的励磁机的增益倍数、所述励磁机的延迟时间常数、所述励磁系统的延迟时间常数、所述电压稳定器的增益、所述电压稳定器的衰减时间常数、所述电压稳定器的延迟时间常数。Optionally, it also includes: according to the real part voltage, imaginary part voltage, time constant, gain multiple of the voltage regulator in the target power grid, and the output voltage of the exciter, the voltage regulator and the stabilizer, and the exciter The gain multiple and delay time constant of the excitation system, as well as the delay time constant of the excitation system, as well as the gain, decay time constant and delay time constant of the voltage stabilizer, construct the excitation system model; the power plant parameters corresponding to the excitation system model include: At least one of the following: the time constant of the voltage regulator, the gain multiple of the voltage regulator, the gain multiple of the exciter of the voltage regulator, the delay time constant of the exciter, the delay of the excitation system time constant, gain of the voltage stabilizer, decay time constant of the voltage stabilizer, delay time constant of the voltage stabilizer.
可选地,所述根据所述目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,包括:确定目标约束条件;基于所述目标约束条件求解所述误差函数,获取与最小误差值所对应的电压预测值、电流预测值和功角预测值。Optionally, the error function is solved according to multiple voltage measurement values, multiple current measurement values and multiple power angle measurement values of the target power grid under multiple preset noise-like disturbance scenarios, to obtain a The corresponding voltage prediction value, current prediction value and power angle prediction value include: determining target constraints; solving the error function based on the target constraints, and obtaining the voltage prediction value and the current prediction value corresponding to the minimum error value and power angle predictions.
可选地,所述确定目标约束条件,包括:根据发电机在预定初始状态下运行的电压实部和电压虚部,所述发电机在所述预定初始状态下运行的电流实部和电流虚部,以及所述发电机的预定电势初始值、功角初值和q轴暂态电势初值,以及所述发电机的q轴同步电抗和d轴暂态电抗,确定所述发电机的有功功率、无功功率和电动势变化量约束条件;基于网络平衡约束,构建所述目标电网中各节点注入电流实部变化量和电流虚部变化量的约束条件;基于预定的电厂参数上限值和预定的电厂参数下限值,构建所述电厂参数的上下限约束条件。Optionally, the determining the target constraint condition includes: according to the real voltage part and the imaginary voltage part of the generator running in a predetermined initial state, the current real part and the current imaginary part of the generator running in the predetermined initial state. part, and the initial value of the predetermined potential, the initial value of the power angle and the initial value of the q-axis transient potential of the generator, as well as the q-axis synchronous reactance and the d-axis transient reactance of the generator, to determine the active power of the generator Constraints on power, reactive power and electromotive force variation; based on network balance constraints, construct constraints on the real and imaginary changes of the injected current at each node in the target grid; based on the predetermined upper limit of power plant parameters and The predetermined lower limit value of the power plant parameter is used to construct the upper and lower limit constraints of the power plant parameter.
可选地,所述多个类噪声扰动场景包括以下至少之二:对所述目标电网中多个电压控制器参考值施加的第一预定扰动、对所述目标电网中多个电厂出口处电力线路的间距施加的第二预定扰动、对多个电厂出口处的负荷施加的第三预定扰动。Optionally, the multiple noise-like disturbance scenarios include at least two of the following: a first predetermined disturbance applied to multiple voltage controller reference values in the target power grid; A second predetermined disturbance applied to the spacing of the lines, and a third predetermined disturbance applied to the loads at the plurality of power plant outlets.
根据本发明实施例的另一方面,还提供了一种电厂参数识别装置,其特征在于,包括:第一获取模块,用于基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下的电压测量值、电流测量值和功角测量值;第二获取模块,用于根据所述目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值;其中,所述误差函数是根据最小二乘法构建得到的,所述误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差,电力参数包括电压、电流和功角;识别模块,用于根据所述与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模块、调速系统模块、励磁系统模块,识别分别与所述发电机三阶模块、调速系统模块、励磁系统模块对应的电厂参数;其中,所述电厂参数为所述发电机三阶模块、调速系统模块、励磁系统模块中的参数。According to another aspect of the embodiments of the present invention, an apparatus for identifying parameters of a power plant is further provided, which is characterized by comprising: a first obtaining module, configured to obtain, based on the power management unit PMU, the target power grid in a plurality of preset categories respectively. voltage measurement value, current measurement value, and power angle measurement value under the noise disturbance scenario; the second acquisition module is configured to obtain multiple voltage measurement values, multiple The error function is solved for the current measurement value and the multiple power angle measurement values to obtain the voltage prediction value, the current prediction value and the power angle prediction value corresponding to the minimum error value; wherein, the error function is constructed and obtained according to the least square method, The error function is used to characterize the error between the measured values of a plurality of power parameters of the power grid and the corresponding predicted values of the power parameters, and the power parameters include voltage, current and power angle; the identification module is used for according to the minimum error value. The corresponding voltage prediction value, current prediction value and power angle prediction value, as well as the pre-built generator third-order module, speed regulation system module, excitation system module, are identified with the generator third-order module, speed regulation system module, Power plant parameters corresponding to the excitation system module; wherein, the power plant parameters are parameters in the generator third-order module, the speed regulation system module, and the excitation system module.
根据本发明实施例的另一方面,还提供了一种计算机可读存储介质,所述存储介质包括存储的程序,其中,在所述程序运行时控制所述存储介质所在设备执行上述任意一项所述的电厂参数识别方法。According to another aspect of the embodiments of the present invention, a computer-readable storage medium is further provided, where the storage medium includes a stored program, wherein when the program is executed, a device where the storage medium is located is controlled to execute any of the foregoing items The power plant parameter identification method.
根据本发明实施例的另一方面,还提供了一种计算机设备,包括:存储器和处理器,所述存储器存储有计算机程序;所述处理器,用于执行所述存储器中存储的计算机程序,所述计算机程序运行时使得所述处理器执行上述任一项所述的电厂参数识别方法。According to another aspect of the embodiments of the present invention, a computer device is also provided, including: a memory and a processor, where the memory stores a computer program; the processor is configured to execute the computer program stored in the memory, When the computer program runs, the processor executes the power plant parameter identification method described in any one of the above.
在本发明实施例中,基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下的电力参数,电力参数包括电压测量值、电流测量值和功角测量值;根据目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,其中,误差函数是根据最小二乘法构建得到的,误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差;根据与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别分别与发电机三阶模型、调速系统模型、励磁系统模型对应的电厂参数;其中,电厂参数为发电机三阶模型、调速系统模型、励磁系统模型中的参数。通过引入多个类噪声扰动,基于多个类噪声扰动进行电厂参数的识别,可提高电厂参数识别效率和准确度,进而解决了相关技术中,存在无法对电厂参数进行准确识别技术问题。In the embodiment of the present invention, based on the power management unit PMU, the power parameters of the target power grid under multiple preset noise-like disturbance scenarios are respectively obtained, and the power parameters include a voltage measurement value, a current measurement value, and a power angle measurement value; according to the target The error function is solved for multiple voltage measurements, multiple current measurements, and multiple power angle measurements of the power grid under multiple preset noise-like disturbance scenarios, and the voltage predicted value and current predicted value corresponding to the minimum error value are obtained. and the predicted value of power angle, in which the error function is constructed according to the least square method, and the error function is used to characterize the error between the measured values of multiple power parameters of the power grid and the corresponding predicted values of power parameters; The corresponding voltage prediction value, current prediction value and power angle prediction value, as well as the pre-built generator third-order model, speed control system model, excitation system model, identify the generator third-order model, speed control system model, excitation system model respectively. The parameters of the power plant corresponding to the model; wherein, the parameters of the power plant are the parameters in the third-order model of the generator, the model of the speed control system, and the model of the excitation system. By introducing multiple noise-like disturbances and identifying power plant parameters based on multiple noise-like disturbances, the efficiency and accuracy of power plant parameter identification can be improved, thereby solving the technical problem of inability to accurately identify power plant parameters in related technologies.
附图说明Description of drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The accompanying drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present application. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the attached image:
图1是根据本发明实施例的一种电厂参数识别方法的流程图;1 is a flowchart of a method for identifying parameters of a power plant according to an embodiment of the present invention;
图2是根据本发明实施例的一种IEEE9节点系统的结构示意图;2 is a schematic structural diagram of an IEEE9 node system according to an embodiment of the present invention;
图3是根据本发明实施例的一种可选的电厂参数识别装置的框架图。FIG. 3 is a frame diagram of an optional power plant parameter identification device according to an embodiment of the present invention.
具体实施方式Detailed ways
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。In order to make those skilled in the art better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only Embodiments are part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。It should be noted that the terms "first", "second" and the like in the description and claims of the present invention and the above drawings are used to distinguish similar objects, and are not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used may be interchanged under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having" and any variations thereof, are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those expressly listed Rather, those steps or units may include other steps or units not expressly listed or inherent to these processes, methods, products or devices.
实施例1Example 1
根据本发明实施例,提供了一种电厂参数识别方法的实施例,需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。According to an embodiment of the present invention, an embodiment of a method for identifying parameters of a power plant is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings may be executed in a computer system such as a set of computer-executable instructions, and, Although a logical order is shown in the flowcharts, in some cases steps shown or described may be performed in an order different from that herein.
图1是根据本发明实施例的电厂参数识别方法,如图1所示,该方法包括如下步骤:Fig. 1 is a power plant parameter identification method according to an embodiment of the present invention. As shown in Fig. 1 , the method includes the following steps:
步骤S102,基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下的电压测量值、电流测量值和功角测量值。Step S102 , based on the power management unit PMU, respectively acquire voltage measurement values, current measurement values and power angle measurement values of the target power grid under multiple preset noise-like disturbance scenarios.
需要明白的是,PMU(Power Management Unit,电源管理单元)是一可用于测量电力参数的测量单元。It should be understood that a PMU (Power Management Unit, power management unit) is a measurement unit that can be used to measure power parameters.
步骤S104,根据目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,其中,误差函数是根据最小二乘法构建得到的,误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差,电力参数包括电压、电流和功角。Step S104, solving an error function according to multiple voltage measurements, multiple current measurements, and multiple power angle measurements of the target power grid under multiple preset noise-like disturbance scenarios, to obtain a voltage prediction corresponding to the minimum error value value, predicted value of current and predicted value of power angle, wherein, the error function is constructed according to the least square method, and the error function is used to characterize the error between the measured values of multiple power parameters of the power grid and the corresponding predicted values of power parameters, and the power Parameters include voltage, current, and power angle.
步骤S106,根据与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别分别与发电机三阶模型、调速系统模型、励磁系统模型对应的电厂参数;其中,电厂参数为发电机三阶模型、调速系统模型、励磁系统模型中的参数。Step S106, according to the voltage prediction value, current prediction value and power angle prediction value corresponding to the minimum error value, as well as the pre-built generator third-order model, speed regulation system model, and excitation system model, identify the third-order generators respectively. The power plant parameters corresponding to the model, the speed control system model, and the excitation system model; wherein, the power plant parameters are the parameters in the generator third-order model, the speed control system model, and the excitation system model.
在本可选实施方式中,基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下的电力参数,电力参数包括电压测量值、电流测量值和功角测量值;根据目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,其中,误差函数是根据最小二乘法构建得到的,误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差;根据与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别分别与发电机三阶模型、调速系统模型、励磁系统模型对应的电厂参数;其中,电厂参数为发电机三阶模型、调速系统模型、励磁系统模型中的参数。通过引入多个不同的类噪声,施加不同类型的扰动,基于不同类型的扰动进行电厂参数的识别,可提高电厂参数识别效率和准确度,进而解决了相关技术中,存在无法对电厂参数进行准确识别技术问题。In this optional embodiment, based on the power management unit PMU, the power parameters of the target power grid under multiple preset noise-like disturbance scenarios are respectively obtained, and the power parameters include a voltage measurement value, a current measurement value, and a power angle measurement value; according to The error function is solved for multiple voltage measurements, multiple current measurements, and multiple power angle measurements of the target power grid under multiple preset noise-like disturbance scenarios, and the voltage prediction value and current prediction value corresponding to the minimum error value are obtained. value and predicted value of power angle, among which, the error function is constructed according to the least square method, and the error function is used to characterize the error between the measured value of multiple power parameters of the power grid and the corresponding predicted value of the power parameter; according to the minimum error value The corresponding voltage prediction value, current prediction value and power angle prediction value, as well as the pre-built generator third-order model, speed control system model and excitation system model, are identified with the generator third-order model, speed control system model, excitation system model respectively. Power plant parameters corresponding to the system model; wherein, the power plant parameters are parameters in the generator third-order model, the speed regulation system model, and the excitation system model. By introducing multiple different types of noise, applying different types of disturbances, and identifying power plant parameters based on different types of disturbances, the efficiency and accuracy of power plant parameter identification can be improved, thereby solving the problem of inability to accurately identify power plant parameters in related technologies. Identify technical issues.
在一些可选实时例中,根据目标电网中发电机的d轴同步电抗、d轴暂态电抗、惯性常数、d轴暂态开路时间常数,以及发电机的功角、加速度、q轴暂态电抗、电流实部和虚部,以及发电机有功功率、励磁电势,以及在噪声扰动场景下节点发电机的交轴电势、有功功率、直轴电流和交轴电流,构建发电机三阶模型;其中,与发电机三阶模型对应的电厂参数包括以下至少之一:发电机的d轴同步电抗、d轴暂态电抗、惯性常数、d轴暂态开路时间常数。由此,可以获取能够反映目标电网中多个电力参数的发电机三阶模型。In some optional real-time examples, according to the d-axis synchronous reactance, d-axis transient reactance, inertia constant, d-axis transient open-circuit time constant of the generator in the target grid, as well as the generator's power angle, acceleration, q-axis transient Reactance, real and imaginary parts of current, as well as generator active power, excitation potential, and quadrature axis potential, active power, direct axis current and quadrature axis current of node generators in the noise disturbance scenario to build a third-order generator model; The power plant parameters corresponding to the third-order generator model include at least one of the following: d-axis synchronous reactance, d-axis transient reactance, inertia constant, and d-axis transient open-circuit time constant of the generator. Thus, a third-order generator model that can reflect multiple power parameters in the target grid can be obtained.
在一些可选实时例中,根据目标电网中调速系统的加速时间常数、调速系统中机组的惯性时间常数和缓冲时间常数、调速系统的参考角频率和运行角频率、调速系统的转子回路电阻、机械转矩,以及调速系统的机械转矩初始值,构建调速系统模型;其中,与调速系统模型所对应的电厂参数包括以下至少之一:调速系统中机组的惯性时间常数、调速系统中机组的缓冲时间常数、调速系统的转子回路电阻。由此,可以获取能够反映目标电网中多个电力参数的调速系统模型。In some optional real-time examples, according to the acceleration time constant of the speed control system in the target power grid, the inertia time constant and buffer time constant of the unit in the speed control system, the reference angular frequency and the operating angular frequency of the speed control system, the speed control system The rotor loop resistance, mechanical torque, and the initial value of the mechanical torque of the speed control system are used to construct a speed control system model; wherein, the power plant parameters corresponding to the speed control system model include at least one of the following: the inertia of the unit in the speed control system Time constant, buffer time constant of the unit in the speed control system, rotor loop resistance of the speed control system. Thus, a speed regulation system model that can reflect multiple power parameters in the target power grid can be obtained.
在一些可选实时例中,根据目标电网中电压调节器的实部电压、虚部电压、时间常数、增益倍数,以及励磁机、电压调节器和稳定器的输出电压,以及励磁机的增益倍数和延迟时间常数,以及励磁系统的延迟时间常数,以及电压稳定器的增益、衰减时间常数和延迟时间常数,构建励磁系统模型;与励磁系统模型所对应的电厂参数包括以下至少之一:电压调节器的时间常数、电压调节器的增益倍数、电压调节器的励磁机的增益倍数、励磁机的延迟时间常数、励磁系统的延迟时间常数、电压稳定器的增益、电压稳定器的衰减时间常数、电压稳定器的延迟时间常数。由此,可以获取能够反映目标电网中多个电力参数的励磁系统模型。In some optional real-time examples, according to the real part voltage, imaginary part voltage, time constant, gain multiplier of the voltage regulator in the target grid, and the output voltage of the exciter, voltage regulator and stabilizer, and the gain multiplier of the exciter and delay time constant, as well as the delay time constant of the excitation system, as well as the gain, decay time constant and delay time constant of the voltage stabilizer, to construct the excitation system model; the power plant parameters corresponding to the excitation system model include at least one of the following: voltage regulation The time constant of the voltage regulator, the gain multiple of the voltage regulator, the gain multiple of the exciter of the voltage regulator, the delay time constant of the exciter, the delay time constant of the excitation system, the gain of the voltage stabilizer, the decay time constant of the voltage stabilizer, Delay time constant for the voltage stabilizer. Thus, an excitation system model that can reflect multiple power parameters in the target grid can be obtained.
在一些可选实时例中,根据目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,包括:确定目标约束条件;基于目标约束条件求解误差函数,获取与最小误差值所对应的电压预测值、电流预测值和功角预测值。In some optional real-time examples, an error function is solved according to multiple voltage measurements, multiple current measurements, and multiple power angle measurements of the target power grid under multiple preset noise-like disturbance scenarios, and the minimum error value is obtained. The corresponding voltage prediction value, current prediction value and power angle prediction value include: determining the target constraint condition; solving the error function based on the target constraint condition, and obtaining the voltage prediction value, current prediction value and power angle prediction corresponding to the minimum error value value.
在一些可选实时例中,确定目标约束条件,包括:根据发电机在预定初始状态下运行的电压实部和电压虚部,发电机在预定初始状态下运行的电流实部和电流虚部,以及发电机的预定电势初始值、功角初值和q轴暂态电势初值,以及发电机的q轴同步电抗和d轴暂态电抗,确定发电机的有功功率、无功功率和电动势变化量约束条件;基于网络平衡约束,构建目标电网中各节点注入电流实部变化量和电流虚部变化量的约束条件;基于预定的电厂参数上限值和预定的电厂参数下限值,构建电厂参数的上下限约束条件。由此,可以获取和目标电网中多个电力参数相关的目标约束条件,基于该目标约束条件,可以获取准确度高且适用性高的电压预测值、电流预测值和功角预测值。In some optional real-time examples, the target constraint conditions are determined, including: according to the voltage real part and the voltage imaginary part of the generator running in the predetermined initial state, the current real part and the current imaginary part of the generator running in the predetermined initial state, As well as the initial value of the predetermined potential, the initial value of the power angle and the initial value of the q-axis transient potential of the generator, as well as the q-axis synchronous reactance and the d-axis transient reactance of the generator, determine the active power, reactive power and electromotive force changes of the generator based on the network balance constraints, construct the constraints of the real and imaginary changes of the injected current at each node in the target power grid; The upper and lower bound constraints of the parameter. Thus, target constraints related to multiple power parameters in the target grid can be obtained, and based on the target constraints, voltage prediction values, current prediction values and power angle prediction values with high accuracy and applicability can be obtained.
多个类噪声扰动场景包括以下至少之二:对目标电网中多个电压控制器参考值施加的第一预定扰动、对目标电网中多个电厂出口处电力线路的间距施加的第二预定扰动、对多个电厂出口处的负荷施加的第三预定扰动。通过施加不同的扰动,基于施加不同扰动的场景,可以准确获取电厂参数。The multiple noise-like disturbance scenarios include at least two of the following: a first predetermined disturbance applied to a plurality of voltage controller reference values in the target power grid, a second predetermined disturbance applied to the distance between the power lines at the outlets of the multiple power plants in the target power grid, A third predetermined disturbance applied to the loads at the plurality of power plant outlets. By applying different disturbances, the power plant parameters can be accurately obtained based on the scenarios in which different disturbances are applied.
基于上述实施例及可选实施例,提供了一种可选实施方式,下面具体说明。Based on the foregoing embodiment and optional embodiments, an optional implementation manner is provided, which will be described in detail below.
电厂参数辨识对电网合环操作、经济运行等具有重要意义。在相关技术中,可以基于傅里叶变化、拉普拉斯变化等技术将物理量转换为频域变量进行分析,这种方法虽然能实现电厂参数辨识,但只能进行离线辨识。相关技术中,还可以通过对物理量进行时域分析,进而实现在线参数辨识。这些方法虽然可以实现对电厂参数的辨识,但这些辨识方法往往依赖故障扰动的发生,其本身对电力系统的稳定运行具备不利影响,且辨识结果具有一定偶然性,可能导致辨识结果的误差较大。The power plant parameter identification is of great significance to the power grid loop closing operation and economic operation. In the related art, physical quantities can be converted into frequency domain variables for analysis based on Fourier transform, Laplace transform and other technologies. Although this method can realize the identification of power plant parameters, it can only be identified offline. In the related art, online parameter identification can also be realized by performing time domain analysis on physical quantities. Although these methods can realize the identification of power plant parameters, these identification methods often rely on the occurrence of fault disturbances, which have adverse effects on the stable operation of the power system, and the identification results are contingent, which may lead to large errors in the identification results.
鉴于此,本可选实施方式引入类噪声技术,从不同角度施加扰动,并进行多扰动场景下的电厂参数辨识,由此提高电厂参数辨识结果的准确度;同时利用电厂参数辨识自由度低、多场景耦合程度低的特点,使用简约空间内点法和情景分解获取电厂参数辨识结果,进而达到显著提高参数辨识效率,实现准确、快速在线辨识电厂参数的目的。In view of this, this optional embodiment introduces noise-like technology, applies disturbances from different angles, and performs power plant parameter identification under multiple disturbance scenarios, thereby improving the accuracy of power plant parameter identification results; Due to the low degree of coupling of multiple scenarios, the simple space interior point method and scenario decomposition are used to obtain the identification results of power plant parameters, thereby significantly improving the efficiency of parameter identification and realizing accurate and rapid online identification of power plant parameters.
本公开实施方式提供的电厂参数识别方法,包括如下步骤:The power plant parameter identification method provided by the embodiment of the present disclosure includes the following steps:
步骤S1:基于PMU装置获取多个不同类噪声扰动场景下的电力数据,电力数据包括电网电压、电流和功角等数据。具体包括如下步骤:Step S1 : based on the PMU device, obtain power data under multiple different types of noise disturbance scenarios, and the power data includes data such as grid voltage, current, and power angle. Specifically include the following steps:
S11:设置不同的类噪声扰动。S11: Set different noise-like disturbances.
需要明白的是,类噪声信号是指系统在正常运行状态下的输出。对电力系统而言,类噪声信号的带宽与机电动态带宽重合,该典型带宽为0.2到2.0Hz。It should be understood that a noise-like signal refers to the output of the system under normal operating conditions. For power systems, the bandwidth of the noise-like signal coincides with the electromechanical dynamic bandwidth, which is typically 0.2 to 2.0 Hz.
在本可选实施方式中,可通过设置如下三类类噪声扰动进行电厂参数识别:分别对不同自动电压控制器参考值施加预设大小的扰动;分别对不同电厂出口处电力线路的间距施加预设大小的扰动;在预设的小范围区间内变动电厂出口处负荷变量。In this optional embodiment, the power plant parameter identification can be performed by setting the following three types of noise disturbances: respectively applying preset disturbances to different automatic voltage controller reference values; Set the size of the disturbance; change the load variable at the outlet of the power plant within a preset small range.
S12:基于PMU装置获取电厂出口电压、电流和功角数据。S12: Based on the PMU device, obtain the data of the output voltage, current and power angle of the power plant.
步骤S2:构建电厂模型,电厂模型包括用于确定发电机待辨识参数的发电机模型、调速系统模型、励磁系统模型在内的电厂模型,根据构建的电场模型确定待辨识的电厂关键参数。具体包括如下步骤:Step S2: constructing a power plant model. The power plant model includes a power plant model including a generator model, a speed regulation system model, and an excitation system model for determining the parameters of the generator to be identified, and the key parameters of the power plant to be identified are determined according to the constructed electric field model. Specifically include the following steps:
S21:构建的发电机模型为三阶模型。发电机的三阶模型如下所示:S21: The constructed generator model is a third-order model. The third-order model of the generator looks like this:
其中:in:
式中:Xdk、X′dk、Mk、T′d0k、Xqk分别发电机d轴同步电抗、d轴暂态电抗、惯性常数、d轴暂态开路时间常数、发电机q轴同步电抗,分别为发电机功角、角速度、q轴暂态电抗、电流实部和虚部;PGk、Efk分别为发电机有功功率和励磁电势;表示情景j下k节点发电机的交轴电势;ωs为同步角速度有名值;表示情景j下k节点发电机在t时刻的有功功率,表示情景j下k节点发电机在t时刻的直轴电流,表示情景j下k节点发电机在t时刻的交轴电流,nG表示节点个数,NS表示故障场景个数。下标k表征发电机节点编号;上标j为不同的类噪声扰动场景编号。In the formula: X dk , X′ dk , M k , T′ d0k , X qk are respectively the generator d-axis synchronous reactance, d-axis transient reactance, inertia constant, d-axis transient open-circuit time constant, and generator q-axis synchronous reactance , are generator power angle, angular velocity, q-axis transient reactance, real part and imaginary part of current; P Gk and E fk are generator active power and excitation potential, respectively; represents the quadrature axis potential of the k-node generator under scenario j; ω s is the named value of the synchronous angular velocity; represents the active power of the k-node generator at time t under scenario j, represents the direct-axis current of the k-node generator at time t under scenario j, represents the quadrature axis current of the k-node generator at time t under scenario j, n G represents the number of nodes, and N S represents the number of fault scenarios. The subscript k represents the generator node number; the superscript j is the number of different noise-like disturbance scenarios.
其中,Xdk、X′dk、Mk、T′d0k为待辨识的发电机参数。Among them, X dk , X' dk , M k , T' d0k are generator parameters to be identified.
S22:不计及调速系统模型中的饱和部分,构建用于确定调速系统待辨识参数的调速系统模型。调速系统模型如下所示:S22: Without taking into account the saturated part of the speed control system model, construct a speed control system model for determining the parameters to be identified for the speed control system. The model of the speed control system is as follows:
式中,tg为调速系统加速时间常数;T1和T2分别为机组惯性时间常数和缓冲时间常数;ωref、ω分别为参考角频率和运行角频率;R为转子回路电阻;Tmech为机械转矩。Tmech0为机械转矩调速前的值。其中,R、T1、T2为待辨识的调速系统参数。In the formula, t g is the acceleration time constant of the speed control system; T 1 and T 2 are the inertia time constant and buffer time constant of the unit respectively; ω ref , ω are the reference angular frequency and the operating angular frequency respectively; R is the rotor loop resistance; T mech is the mechanical torque. T mech0 is the value before mechanical torque speed regulation. Among them, R, T 1 , T 2 are the parameters of the speed control system to be identified.
S23:不计及励磁系统模型中的饱和部分,构建确定励磁系统待辨识参数的四阶励磁系统模型。四阶励磁系统模型如下所示:S23: Without taking into account the saturation part in the excitation system model, construct a fourth-order excitation system model that determines the parameters to be identified for the excitation system. The fourth-order excitation system model is shown below:
式中,Vm为测量电压;Vx、Vy、Ta、K。分别为电压调节器实部电压、虚部电压、时间常数、增益倍数;Vr1、Vr2、Vf分别为励磁机、电压调节器和稳定器的输出电压;Vref为参考电压;Kf、Tf分别为励磁机的增益倍数和延迟时间常数;Tr为整个励磁系统的延迟时间常数;Ae,Be,Te分别为电压稳定器的增益、衰减时间常数和延迟时间常数。In the formula, V m is the measurement voltage; V x , V y , Ta , K. V r1 , V r2 , V f are the output voltages of the exciter, voltage regulator and stabilizer, respectively; V ref is the reference voltage; K f , T f are the gain multiple and delay time constant of the exciter, respectively; Tr is the delay time constant of the entire excitation system; A e , Be , T e are the gain, decay time constant and delay time constant of the voltage stabilizer, respectively.
其中,Ka,Ta,Kf,Tf,Tr,Te,Ae,Be为待辨识的励磁系统参数。Among them, K a , T a , K f , T f , Tr , T e , A e , B e are the excitation system parameters to be identified.
步骤S3:以电压、电流和功角模型值和量测值误差最小为目标建立电厂参数辨识模型。其中,模型值即为根据上述发电机模型、调速系统模型和励磁系统模型得到的对应预测值;量测值为根据PMU装置采集的值。具体包括如下步骤:Step S3: Establish a power plant parameter identification model with the goal of minimizing errors between the voltage, current and power angle model values and the measured values. Wherein, the model value is the corresponding predicted value obtained according to the above generator model, speed regulation system model and excitation system model; the measured value is the value collected according to the PMU device. Specifically include the following steps:
基于ARIMA(Auto Regression and Moving Average,自回归移动平均模型)对量测电压、电流和功角进行拟合,并以三个运行变量的拟合曲线与模型预测曲线的最小二乘误差最小为目标,建立多类噪声场景下的多电厂关键参数动态辨识模型(相当于前述实施例中的误差函数):Based on ARIMA (Auto Regression and Moving Average, autoregressive moving average model), the measured voltage, current and power angle are fitted, and the goal is to minimize the least squares error between the fitting curve of the three operating variables and the model prediction curve , establish a dynamic identification model of key parameters of multiple power plants under multiple types of noise scenarios (equivalent to the error function in the foregoing embodiment):
式中,nT为时段数;ngenl为发电机节点数目;Ns为故障场景总数;分别为情景j下k节点发电机功角的预测值和测量值;分别为电流实部的预测值和测量值、电流虚部的预测值和测量值; 分别为电压实部的预测值和测量值、电压虚部的预测值和测量值。In the formula, nT is the number of time periods; n genl is the number of generator nodes; Ns is the total number of fault scenarios; are the predicted value and the measured value of the generator power angle of node k under scenario j, respectively; are the predicted and measured values of the real part of the current and the predicted and measured values of the imaginary part of the current, respectively; are the predicted and measured values of the real part of the voltage, and the predicted and measured values of the imaginary part of the voltage, respectively.
由于发电机电压、电流和功角量纲不同,且数据的数量级有所差别,故采用相对值的方式表征误差,将上述目标函数修改为:Since the dimensions of the generator voltage, current and power angle are different, and the magnitude of the data is different, the relative value is used to characterize the error, and the above objective function is modified as:
步骤S4:确定辨识模型约束条件。具体包括如下步骤:Step S4: Determine the constraints of the identification model. Specifically include the following steps:
步骤S41:假设各类噪声扰动情景下的初始状态均相同,其初值约束为:Step S41: Assuming that the initial states under various noise disturbance scenarios are the same, the initial value constraints are:
式中:ΔPk,0、ΔQk,0、ΔE′qk,0分别为节点k的有功、无功和电动势变化量;ek,fk,Iek,Ifk分别为发电机初态运行电压实部、虚部,电流实部、虚部;EQk,δk,0,E′qk,0为发电机虚拟电势初始值、功角初值和q轴暂态电势初值;Xqk,X′dk分别为发电机q轴同步电抗和d轴暂态电抗。In the formula: ΔP k,0 , ΔQ k,0 , ΔE′ qk,0 are the changes of active power, reactive power and electromotive force of node k respectively; ek , f k , I ek , I fk are the initial state operation of the generator, respectively Real and imaginary parts of voltage, real and imaginary parts of current; E Qk , δ k,0 , E′ qk,0 are the initial value of generator virtual potential, initial value of power angle and initial value of q-axis transient potential; X qk , X′ dk are the generator q-axis synchronous reactance and d-axis transient reactance, respectively.
步骤S42:各类噪声扰动情景下均满足网络平衡约束,如下:Step S42: The network balance constraint is satisfied under various noise disturbance scenarios, as follows:
其中:in:
式中:表示节点注入电流实部变化量和电流虚部变化量;分别为网络导纳矩阵实部和虚部;分别为电压实部和虚部,注入电流实部和虚部;分别为发电机在t时刻的电压实部和虚部。Xqk分别表示发电机在t时刻q轴暂态电势值、d轴暂态电抗、功角值和q轴电抗。where: Represents the change of the real part of the injected current and the change of the imaginary part of the current; are the real and imaginary parts of the network admittance matrix, respectively; are the real and imaginary parts of the voltage, and the real and imaginary parts of the injected current, respectively; are the real and imaginary parts of the generator voltage at time t, respectively. X qk respectively represent the q-axis transient potential value, d-axis transient reactance, power angle value and q-axis reactance of the generator at time t.
步骤S43:待辨识参数均满足上下限约束:Step S43: the parameters to be identified all satisfy the upper and lower limit constraints:
发电机待辨识参数上下限约束:Upper and lower limit constraints of generator parameters to be identified:
调速系统待辨识参数上下限约束为:The upper and lower limit constraints of the parameters to be identified in the speed control system are:
励磁系统待辨识参数上下限约束为:The upper and lower limit constraints of the parameters to be identified in the excitation system are:
式中,下标min表示对应参数的下限,下标max表示对应参数的上限。In the formula, the subscript min represents the lower limit of the corresponding parameter, and the subscript max represents the upper limit of the corresponding parameter.
步骤S5:基于简约空间内点法和情景分解方法对上述电厂参数辨识模型进行求解,得到最终电厂关键参数值。Step S5: Solve the above-mentioned power plant parameter identification model based on the reduced space interior point method and the scenario decomposition method, and obtain the final key parameter values of the power plant.
在本可选实施方式中,通过基于PMU装置获取不同类噪声扰动场景下电网电压、电流和功角等数据;建立包括发电机模型、调速系统模型、励磁系统模型在内的电厂模型,以确定待辨识的电厂关键参数;以电压、电流和功角模型值和量测值误差最小为目标建立电厂参数辨识模型;确定辨识模型约束条件;基于简约空间内点法和情景分解方法对上述电厂参数辨识模型进行求解,得到最终电厂关键参数值。能够从不同角度施加扰动类型,进行多扰动场景辨识,可以提高辨识的准确度;同时,利用参数辨识自由度低的特点、多场景耦合程度低,提出了简约空间内点法和情景分解求解策略,可显著提高参数辨识效率,实现准确、快速在线辨识。In this optional embodiment, data such as grid voltage, current and power angle under different types of noise disturbance scenarios are obtained based on the PMU device; power plant models including generator models, speed regulation system models, and excitation system models are established to Determine the key parameters of the power plant to be identified; establish the power plant parameter identification model with the goal of minimizing the errors of the voltage, current and power angle model values and measurement values; determine the constraints of the identification model; The parameter identification model is solved to obtain the final key parameter values of the power plant. It can apply disturbance types from different angles and perform multi-disturbance scene identification, which can improve the accuracy of identification. At the same time, taking advantage of the characteristics of low degree of freedom in parameter identification and low degree of multi-scene coupling, a simplified space interior point method and a scenario decomposition solution strategy are proposed. , which can significantly improve the efficiency of parameter identification and achieve accurate and fast online identification.
下面以IEEE9节点系统为例进行进一步说明。图2是IEEE9节点系统的示意图。图2所示,发电机连接在节点1、2、3,负荷连接在节点5、6、8。The IEEE9 node system is taken as an example for further description below. Figure 2 is a schematic diagram of an IEEE9 node system. As shown in Figure 2, the generators are connected at
其中,IEEE9节点系统中发电机的实际参数值如表1所示。Among them, the actual parameter values of the generator in the IEEE9 node system are shown in Table 1.
表1Table 1
图2是IEEE9节点系统的示意图。图2所示,发电其中,IEEE9节点系统中励磁系统的实际参数值如表2所示。Figure 2 is a schematic diagram of an IEEE9 node system. As shown in Figure 2, the actual parameter values of the excitation system in the IEEE9 node system are shown in Table 2.
表2Table 2
其中,IEEE9节点系统中调速系统的实际参数值如表3所示。Among them, the actual parameter values of the speed control system in the IEEE9 node system are shown in Table 3.
表3table 3
预设的各类类噪声扰动情景包括:对自动电压控制器参考值施加0.1的扰动;分别对线路4-6、线路9-8、线路7-5的线路间距施加一定扰动,即线路阻抗施加扰动,改变0.1;分别在负荷节点5、6、8施加小负荷扰动2MW;共计7个扰动场景。测试环境为VS2008,运行平台为CPUi5,CPU主频为3.19GHz,内存为4GB。The preset various types of noise disturbance scenarios include: applying a disturbance of 0.1 to the reference value of the automatic voltage controller; applying a certain disturbance to the line spacing of lines 4-6, 9-8, and 7-5 respectively, that is, applying line impedance. Disturbance, change 0.1; apply a small load disturbance of 2MW to load
在本可选实时方式中,优化问题规模如表4所示。In this optional real-time approach, the optimization problem size is shown in Table 4.
表4Table 4
在表4中,Ns表示情景数目,n和m分别是原对偶系统状态变量数目和等式约束数目,DIM和NNZ分别是原对偶系统维数和非零元数目,DOF是自由度。由表4可知,各算例优化规模DIM很大,但是其自由度DOF很低,非常适合简约空间下求解。In Table 4, Ns represents the number of scenarios, n and m are the number of state variables and equality constraints of the primal-dual system, respectively, DIM and NNZ are the dimension of the primal-dual system and the number of non-zero elements, and DOF is the degree of freedom. It can be seen from Table 4 that the optimization scale DIM of each calculation example is large, but its DOF is very low, which is very suitable for solving in a simple space.
基于本可选实施方式识别的发电机参数结果如表5所示。The results of the generator parameters identified based on this optional embodiment are shown in Table 5.
表5table 5
基于本可选实施方式识别的励磁参数结果如表6所示:The results of the excitation parameters identified based on this optional embodiment are shown in Table 6:
表6Table 6
基于本可选实施方式识别的调速参数结果如表7所示:The results of the speed regulation parameters identified based on this optional implementation are shown in Table 7:
表7Table 7
根据上述电厂参数的识别结果可知,在传统单扰动的情境下,由于抵抗测量误差的能力较差,辨识结果较差,当扰动情景数增为7个时,整体辨识算法辨识精度有所提高,极个别发电机参数误差较大,基本能反映系统的动态响应过程。According to the identification results of the above power plant parameters, in the traditional single disturbance scenario, due to the poor ability to resist measurement errors, the identification results are poor. When the number of disturbance scenarios increases to 7, the overall identification algorithm identification accuracy is improved. Very few generator parameters have large errors, which can basically reflect the dynamic response process of the system.
电厂关键参数辨识过程一般是以特定物理量量测值与模型值误差最小为目标的优化过程,目前主要包括频域辨识方法和时域辨识方法两大类。前者主要是基于傅里叶变化、拉普拉斯变化等技术将物理量转换为频域变量进行分析,具有计算简单、方便的优势,但只能进行离线辨识;后者直接对物理量进行时域分析,可以实现在线参数辨识,如运用最小二乘算法对三阶同步发电机进行在线辨识等。但上述辨识方法往往依赖故障扰动的发生,其本身对电力系统的稳定运行具备不利影响,且单一扰动下的参数辨识具有一定偶然性,可能导致误差较大。在本可选实施方式中,通过引入多个类噪声,从不同角度施加扰动类型,并进行多扰动场景的电厂参数辨识,可以提高电厂参数辨识的准确度,同时,利用参数辨识自由度低的特点、多场景耦合程度低,提出了简约空间内点法和情景分解求解策略,可显著提高参数辨识效率,实现准确、快速在线辨识。The identification process of key parameters of a power plant is generally an optimization process with the goal of minimizing the error between the measured value of a specific physical quantity and the model value. The former is mainly based on Fourier transformation, Laplace transformation and other technologies to convert physical quantities into frequency domain variables for analysis, which has the advantages of simple and convenient calculation, but can only be identified offline; the latter directly analyzes physical quantities in time domain , can realize online parameter identification, such as using the least squares algorithm to identify the third-order synchronous generator online. However, the above identification methods often rely on the occurrence of fault disturbance, which itself has an adverse effect on the stable operation of the power system, and the parameter identification under a single disturbance has certain contingency, which may lead to large errors. In this optional embodiment, by introducing multiple types of noise, applying disturbance types from different angles, and performing power plant parameter identification in multiple disturbance scenarios, the accuracy of power plant parameter identification can be improved. Due to its characteristics and the low degree of multi-scenario coupling, a simple space interior point method and a scenario decomposition solution strategy are proposed, which can significantly improve the efficiency of parameter identification and achieve accurate and fast online identification.
实施例2Example 2
根据本发明实施例,还提供了一种电厂参数识别装置的结构框图。参照图3所示,装置包括第一获取模块302、第二获取模块304、识别模块306。下面具体说明。According to an embodiment of the present invention, a structural block diagram of a power plant parameter identification device is also provided. Referring to FIG. 3 , the apparatus includes a first acquisition module 302 , a second acquisition module 304 , and an identification module 306 . The specific description is given below.
第一获取模块302,用于基于电源管理单元PMU,分别获取目标电网在预设的多个类噪声扰动场景下的电压测量值、电流测量值和功角测量值;第二获取模块304,连接于上述第一获取模块302,用于根据目标电网在预设的多个类噪声扰动场景下的多个电压测量值、多个电流测量值和多个功角测量值求解误差函数,得到与最小误差值所对应的电压预测值、电流预测值和功角预测值,其中,误差函数是根据最小二乘法构建得到的,误差函数用于表征电网的多个电力参数实测值和对应的电力参数预测值之间的误差,电力参数包括电压、电流和功角;识别模块306,连接于上述第二获取模块304,用于根据与最小误差值所对应的电压预测值、电流预测值和功角预测值,以及预构建的发电机三阶模型、调速系统模型、励磁系统模型,识别分别与发电机三阶模型、调速系统模型、励磁系统模型对应的电厂参数;其中,电厂参数为发电机三阶模型、调速系统模型、励磁系统模型中的参数。The first acquisition module 302 is configured to acquire, based on the power management unit PMU, the voltage measurement value, current measurement value and power angle measurement value of the target power grid under multiple preset noise-like disturbance scenarios respectively; the second acquisition module 304 is connected to In the above-mentioned first acquisition module 302, it is used to solve the error function according to the multiple voltage measurement values, multiple current measurement values and multiple power angle measurement values of the target power grid under multiple preset noise-like disturbance scenarios, to obtain the minimum value. The voltage predicted value, current predicted value and power angle predicted value corresponding to the error value, wherein, the error function is constructed according to the least square method, and the error function is used to characterize the measured values of multiple power parameters of the power grid and the corresponding predicted power parameters The error between the values, the power parameters include voltage, current and power angle; the identification module 306, connected to the above-mentioned second acquisition module 304, is used for predicting the voltage, current and power angle corresponding to the minimum error value. value, and the pre-built generator third-order model, speed regulation system model, and excitation system model, and identify the power plant parameters corresponding to the generator third-order model, speed regulation system model, and excitation system model; among them, the power plant parameter is the generator Parameters in the third-order model, the speed control system model, and the excitation system model.
此处需要说明的是,上述第一获取模块302、第二获取模块304、识别模块306对应于实施例1中的步骤S102至步骤S106,三个模块与对应的步骤所实现的实例和应用场景相同,但不限于上述实施例1所公开的内容。It should be noted here that the above-mentioned first acquisition module 302, second acquisition module 304, and identification module 306 correspond to steps S102 to S106 in
根据本发明实施例的另一方面,提供了一种计算机可读存储介质,存储介质包括存储的程序,其中,在程序运行时控制存储介质所在设备执行上述任意一项的电厂参数识别方法。According to another aspect of the embodiments of the present invention, a computer-readable storage medium is provided, the storage medium includes a stored program, wherein when the program is run, a device where the storage medium is located is controlled to execute any one of the above power plant parameter identification methods.
根据本发明实施例的另一方面,还提供了一种计算机设备,包括:存储器和处理器,存储器存储有计算机程序;处理器,用于执行存储器中存储的计算机程序,计算机程序运行时使得处理器执行上述任一项的电厂参数识别方法。According to another aspect of the embodiments of the present invention, a computer device is also provided, including: a memory and a processor, where the memory stores a computer program; the processor is configured to execute the computer program stored in the memory, and when the computer program runs, the processing The controller executes any of the power plant parameter identification methods described above.
采用本发明实施例,提供了一种图像处理的方案。通过,从而达到了目的,进而解决了相关技术中,的技术问题。With the embodiments of the present invention, an image processing solution is provided. Through, the purpose is achieved, and the technical problem in the related art is solved.
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令终端设备相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:闪存盘、只读存储器(Read-Only Memory,ROM)、随机存取器(RandomAccess Memory,RAM)、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing the hardware related to the terminal device through a program, and the program can be stored in a computer-readable storage medium, and the storage medium can Including: flash disk, read-only memory (Read-Only Memory, ROM), random access device (RandomAccess Memory, RAM), magnetic disk or optical disk, etc.
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。The above-mentioned serial numbers of the embodiments of the present invention are only for description, and do not represent the advantages or disadvantages of the embodiments.
在本发明的上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其他实施例的相关描述。In the above-mentioned embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
在本申请所提供的几个实施例中,应该理解到,所揭露的技术内容,可通过其它的方式实现。其中,以上所描述的装置实施例仅仅是示意性的,例如单元的划分,可以为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,单元或模型的间接耦合或通信连接,可以是电性或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. The device embodiments described above are only illustrative, for example, the division of units may be a logical function division, and there may be other division methods in actual implementation, for example, multiple units or components may be combined or integrated into Another system, or some features can be ignored, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be indirect coupling or communication connection through some interfaces, units or models, and may be in electrical or other forms.
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。The units described as separate components may or may not be physically separated, and components shown as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution in this embodiment.
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically alone, or two or more units may be integrated into one unit. The above-mentioned integrated units may be implemented in the form of hardware, or may be implemented in the form of software functional units.
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。The integrated unit, if implemented in the form of a software functional unit and sold or used as an independent product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium , including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes .
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above are only the preferred embodiments of the present invention. It should be pointed out that for those skilled in the art, without departing from the principles of the present invention, several improvements and modifications can be made. It should be regarded as the protection scope of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210764472.9A CN115117879A (en) | 2022-06-30 | 2022-06-30 | Power plant parameter identification method and device and computer readable storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210764472.9A CN115117879A (en) | 2022-06-30 | 2022-06-30 | Power plant parameter identification method and device and computer readable storage medium |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115117879A true CN115117879A (en) | 2022-09-27 |
Family
ID=83329285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210764472.9A Pending CN115117879A (en) | 2022-06-30 | 2022-06-30 | Power plant parameter identification method and device and computer readable storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115117879A (en) |
-
2022
- 2022-06-30 CN CN202210764472.9A patent/CN115117879A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9574511B2 (en) | System and method for a load anticipation feature and its tuning method for a generating set | |
Saha et al. | Fault detection in transmission networks of power systems | |
Mohammadi-Ivatloo et al. | Online small signal stability analysis of multi-machine systems based on synchronized phasor measurements | |
US11616365B2 (en) | Decentralized hardware-in-the-loop scheme | |
Zaker et al. | Small signal equivalent model of synchronous generator-based grid-connected microgrid using improved Heffron-Phillips model | |
CN109962495A (en) | A method for locating and suppressing ultra-low frequency oscillation disturbance source | |
CN103336882B (en) | A kind of overall process dynamic voltage stability margin appraisal procedure based on time-domain-simulation | |
Ali et al. | A trajectory piecewise-linear approach to nonlinear model order reduction of wind farms | |
WO2022231648A1 (en) | Power system model calibration using measurement data | |
Masrob et al. | Design of a neural network based power system stabilizer in reduced order power system | |
Stadler et al. | Application of dynamic REI reduction | |
CN110474322B (en) | Multi-machine system small interference stability evaluation method comprising virtual synchronous generator | |
CN115117879A (en) | Power plant parameter identification method and device and computer readable storage medium | |
CA3177640A1 (en) | System and method for automated tuning/configuring of a power system stabilizer (pss) in a digital excitation control system | |
CN116845886B (en) | Multi-port autonomous photovoltaic system network construction control method based on model prediction | |
Gong et al. | Parametric identification of DQ impedance model for three-phase voltage-source converters | |
CN117895529A (en) | Frequency response analysis method based on equivalent inertia characterization model of wind turbine | |
Ramírez-Betancour et al. | A two-time scale simulation for dynamic analysis of power systems | |
Sajjadi et al. | Participation factor-based adaptive model reduction for fast power system simulation | |
CN117175549A (en) | Method and system for analyzing frequency response of micro-grid with grid-structured converter access | |
CN116031901A (en) | Time-varying subsynchronous oscillation suppression strategy, system and device for doubly-fed wind farm | |
Kamwa et al. | Sensorless ANN-based speed estimation of synchronous generators: Improved performance through physically motivated pre-filters | |
CN107147112B (en) | A kind of power system modeling method and system based on weighting polymerization | |
CN103336864A (en) | On-line recognizing method for multistep recursion least square of dynamic equivalent model of multimachine interconnected electric power system | |
Xing et al. | On-line identification of equivalent inertia for DFIG wind turbines based on extended Kalman filters |
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 |