CN107133749B - Power information physical coupling modeling method considering demand response information - Google Patents

Power information physical coupling modeling method considering demand response information Download PDF

Info

Publication number
CN107133749B
CN107133749B CN201710537021.0A CN201710537021A CN107133749B CN 107133749 B CN107133749 B CN 107133749B CN 201710537021 A CN201710537021 A CN 201710537021A CN 107133749 B CN107133749 B CN 107133749B
Authority
CN
China
Prior art keywords
demand response
power
information
physical coupling
frequency
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
CN201710537021.0A
Other languages
Chinese (zh)
Other versions
CN107133749A (en
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.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Nari Technology Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Nari Technology 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 State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Nari Technology Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201710537021.0A priority Critical patent/CN107133749B/en
Publication of CN107133749A publication Critical patent/CN107133749A/en
Application granted granted Critical
Publication of CN107133749B publication Critical patent/CN107133749B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067—Enterprise or organisation modelling
    • 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
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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
    • 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]
    • 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
    • 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
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention discloses a power information physical coupling modeling method considering demand response information, and belongs to the technical field of power information physical system analysis. Aiming at the problem of frequency stability control in the power information physical coupling system, the invention provides a method for participating in the system frequency control process by using demand response resources, a demand response feedback control loop is added into a traditional frequency response model, and the feedback loop establishes a power information physical coupling model by considering the characteristics of the demand response participating in the system frequency modulation capability, the information system communication delay and the like. The invention realizes the close connection between the information side and the physical side in the power information physical coupling system, and greatly improves the frequency stability of the power system.

Description

Power information physical coupling modeling method considering demand response information
Technical Field
The invention belongs to the technical field of analysis of power information physical systems, and particularly relates to a power information physical coupling modeling method considering demand response information.
Background
The Cyber-Physical Systems (CPS), also called Cyber-Physical Systems (Cyber-Physical Systems) and Cyber-Physical coupling Systems (Cyber-Physical Systems), refer to a multidimensional heterogeneous complex system which integrates a computing system, a communication network and a Physical environment into a whole by a 3C technology to form real-time perception, dynamic control and information service fusion. With the rapid development of computing technology, communication technology and intelligent control technology, once the information physical coupling system is proposed, it draws great attention to and keeps developing rapidly in academic and industrial fields.
In recent years, with the continuous development of smart grid construction, the number of grid sensors, the scale of an information network and the number of decision units are all rapidly increased, and the automation degree of a power system is greatly improved. The powerful functions of the information system provide technical support for the operation of the power grid, and the requirement response technology can participate in the frequency adjustment of the system. Therefore, how to consider the physical interaction characteristics of the power information and mine the effective information of a large amount of data in the power information physical coupling system, a power information physical coupling modeling method considering the demand response information is researched, and a problem to be researched is urgently needed. However, no mature protocol has been proposed for this study.
Disclosure of Invention
The invention aims to: aiming at the defect that demand response information is not considered in the process of power information physical coupling modeling in the prior art, a power information physical coupling modeling method considering the demand response information is provided, and a new thought is provided for frequency stability control of an information physical coupling system.
Specifically, the invention is realized by adopting the following technical scheme: the method comprises the following steps of incorporating demand response information into a power system physical model, and establishing a power information physical coupling model as shown in the following formula:
ΔPT(s)-ΔPL(s)+ΔPDR(s)＝2H·s·Δω(s)+D·Δω(s)
in the formula,. DELTA.PT(s) representing the power increment of the thermal power generating unit; delta PL(s) represents the load power increment；ΔPDR(s) represents a demand response resource power increment; H. d represents the system inertia time constant and the damping coefficient respectively; Δ ω(s) represents the system frequency increment; s represents the Laplace operator;
considering that the power change of the demand response resource participating in the frequency control of the power information physical coupling system is instantly completed, the delta P is adjustedDR(s) is characterized in that the transfer function form of a proportional control delay link is as follows:
Figure GDA0002551243640000021
in the formula, KPThe regulation relation of the demand response power and the frequency variation is represented;
Figure GDA0002551243640000022
representing a delay link; t isdRepresenting a demand response resource response delay time constant.
The technical scheme is further characterized in that a Pade approximation (Pade approximation) response delay link to the demand is adopted
Figure GDA0002551243640000023
Linearization was performed as follows:
Figure GDA0002551243640000024
wherein G(s) represents
Figure GDA0002551243640000025
And (5) the transfer function of the link after linearization.
The invention has the following beneficial effects: according to the invention, a demand response feedback control loop is added in the traditional frequency response model, and the feedback loop considers the capacity of demand response participating in system frequency modulation and the characteristics of information system communication delay and the like, so that a power information physical coupling model is established, the close relation between an information side and a physical side in a power information physical coupling system is realized, and the frequency stability of the power system is greatly improved.
Drawings
FIG. 1 is a block diagram of a power information physical coupling model of the present invention that considers demand response information;
FIG. 2 is a diagram of a simulation model according to an embodiment.
FIG. 3 is a diagram of simulation results of an embodiment.
Detailed Description
The present invention will be described in further detail with reference to the following examples and the accompanying drawings.
Example 1:
the invention discloses a power information physical coupling modeling method considering demand response information, and particularly relates to a power information physical coupling modeling method which includes the following main processes of incorporating the demand response information into a power system physical model and establishing a power information physical coupling model as shown in the following formula.
ΔPT(s)-ΔPL(s)+ΔPDR(s)＝2H·s·Δω(s)+D·Δω(s)
Wherein, Δ PT(s) representing the power increment of the thermal power generating unit; delta PL(s) represents a load power increment; delta PDR(s) represents a demand response resource power increment; H. d represents the system inertia time constant and the damping coefficient respectively; Δ ω(s) represents the system frequency increment; s represents the laplace operator.
Considering that the power change of the demand response resource participating in the frequency control of the power information physical coupling system is instantly completed, the delta P is adjustedDRAnd(s) is characterized in a transfer function form of a proportional control delay element, and a specific formula is as follows.
Figure GDA0002551243640000031
Wherein, KPThe regulation relation of the demand response power and the frequency variation is represented;
Figure GDA0002551243640000032
representing a delay link; t isdIndicating demand response resource response latencyA time constant.
Delay link for responding to demand by using Pade approximation (Pade approximation)
Figure GDA0002551243640000033
Linearization was performed as follows.
Figure GDA0002551243640000034
Wherein:
Figure GDA0002551243640000035
is a polynomial of order p,
Figure GDA0002551243640000036
Figure GDA0002551243640000037
is a polynomial of order q,
Figure GDA0002551243640000038
p and q are polynomial numbers, generally 5-10, in the embodiment, 5; k is a non-negative integer.
By adopting five-stage simplification to obtain
Figure GDA0002551243640000039
Wherein G(s) represents
Figure GDA0002551243640000041
And (5) the transfer function of the link after linearization.
Through the modeling process, the present embodiment finally forms a complete model block diagram as shown in fig. 1. As shown in FIG. 1, Δ PL(s) represents a load power increment; kPThe regulation relation of the demand response power and the frequency variation is represented;
Figure GDA0002551243640000042
representing a delay link; t isdA response delay time constant representing the demand response resource; H. d represents the system inertia time constant and the damping coefficient respectively; fHPRepresents a reheat coefficient; t isRHRepresenting a reheat generator set time constant; t isGRepresents the governor time constant; t ischRepresenting the turbine time constant; Δ ω(s) represents the system frequency increment; r represents a primary frequency modulation difference adjustment coefficient; s represents the laplace operator.
Practical applications of the above scheme are given below. Assume that the system parameters are as shown in the table below.
TABLE 1 simulation parameters of electric power systems
TG Tch TRH R FHP 2H D Td △PL(s) Kp
0.2s 0.4s 7s 3.0Hz/p.u. 0.3 0.1667pu.s 0.015p.u./Hz 0.1s 0.01p.u. 0.2
Wherein, TGRepresents the governor time constant; t ischRepresenting the turbine time constant; t isRHRepresenting the inertia time constant of the reheating unit; r represents a primary frequency modulation difference adjustment coefficient; fHPRepresents a reheat coefficient; H. d represents the system inertia time constant and the damping coefficient respectively; t isdA response delay time constant representing the demand response resource; delta PL(s) represents a load power increment; kPIndicating the regulation relationship of the demand response power and the frequency variation.
Substituting the parameters in the table into the complete model block diagram in fig. 1, and building a simulation model in Simulink/Matlab, wherein the structure of the simulation model is shown in fig. 2, outputting a frequency change signal from a frequency change module to represent that the frequency of the power grid falls, and outputting a frequency recovery curve through a frequency modulation module of the power grid and a proportion and delay module added with demand response.
The simulation results are shown in fig. 3. According to simulation results, the effect of considering the participation of demand response to the frequency primary frequency modulation is better than that of not considering the demand response, the lowest point of frequency drop is improved, and the final recovery stable value of the frequency is higher than that of not considering the demand response; considering the effect of communication delay on the participation of demand response in primary frequency modulation, it can be seen that in the case of communication delay, the lowest frequency point is lower than that in the case of no consideration of communication delay, but the final frequency recovery effect is similar.
Although the present invention has been described in terms of the preferred embodiment, it is not intended that the invention be limited to the embodiment. Any equivalent changes or modifications made without departing from the spirit and scope of the present invention also belong to the protection scope of the present invention. The scope of the invention should therefore be determined with reference to the appended claims.

Claims (1)

1. A power information physical coupling modeling method considering demand response information is characterized in that the demand response information is incorporated into a power system physical model, and the power information physical coupling model is established as shown in the following formula:
ΔPT(s)-ΔPL(s)+ΔPDR(s)＝2H·s·Δω(s)+D·Δω(s)
in the formula,. DELTA.PT(s) representing the power increment of the thermal power generating unit; delta PL(s) represents a load power increment; delta PDR(s) represents a demand response resource power increment; H. d represents the system inertia time constant and the damping coefficient respectively; Δ ω(s) represents the system frequency increment; s represents the Laplace operator;
wherein, Δ PDR(s) is characterized in that the transfer function form of a proportional control delay link is as follows:
Figure FDA0002551243630000011
in the formula, KPThe regulation relation of the demand response power and the frequency variation is represented;
Figure FDA0002551243630000012
representing a delay link; t isdRepresenting the response delay time constant of the demand response resource and adopting Pade approximation to delay the demand response
Figure FDA0002551243630000013
Linearization was performed as follows:
Figure FDA0002551243630000014
wherein G(s) represents
Figure FDA0002551243630000015
And (5) the transfer function of the link after linearization.
CN201710537021.0A 2017-07-04 2017-07-04 Power information physical coupling modeling method considering demand response information Active CN107133749B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710537021.0A CN107133749B (en) 2017-07-04 2017-07-04 Power information physical coupling modeling method considering demand response information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710537021.0A CN107133749B (en) 2017-07-04 2017-07-04 Power information physical coupling modeling method considering demand response information

Publications (2)

Publication Number Publication Date
CN107133749A CN107133749A (en) 2017-09-05
CN107133749B true CN107133749B (en) 2020-11-03

Family

ID=59735975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710537021.0A Active CN107133749B (en) 2017-07-04 2017-07-04 Power information physical coupling modeling method considering demand response information

Country Status (1)

Country Link
CN (1) CN107133749B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108021029A (en) * 2017-11-17 2018-05-11 北京航空航天大学 A kind of intelligent domestic electricity demanding response platform
CN107947169B (en) * 2017-12-07 2020-01-10 清华大学 Information flow modeling method for power grid energy management system
CN108489015B (en) * 2018-03-05 2020-07-14 苏州科技大学 Air conditioning system temperature control method based on pole allocation and Pade approximation
CN109345115B (en) * 2018-09-30 2021-08-06 国网浙江省电力有限公司 Evaluation method for evolution driving metric value of power information physical system
CN109755937B (en) * 2019-01-18 2021-08-31 内蒙古电力(集团)有限责任公司电力调度控制分公司 Regional power grid inertia calculation method and device based on measurement
CN111368443A (en) * 2020-03-06 2020-07-03 北京旷视机器人技术有限公司 Equipment simulation method and device and electronic equipment
CN112531791B (en) * 2020-11-27 2023-03-24 国网宁夏电力有限公司电力科学研究院 Method and device for controlling resource coordination frequency of demand side
CN112818545B (en) * 2021-02-03 2023-03-21 东北电力大学 Power information physical joint simulation system based on OPC
CN113078650B (en) * 2021-04-25 2023-03-31 上海交通大学 Active power distribution network information physical fusion control method of hybrid system model

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549869A (en) * 2009-07-27 2012-07-04 响应负载有限公司 Responsive load monitoring system and method
CN103527407A (en) * 2012-06-28 2014-01-22 西门子公司 Method and arrangement for responding to a grid event, such as fast frequency drop, by combining demand response, inertial response and spinning reserve
JP2015154595A (en) * 2014-02-14 2015-08-24 アズビル株式会社 Energy management device and energy management method
EP2937962A1 (en) * 2012-12-18 2015-10-28 Panasonic Intellectual Property Management Co., Ltd. Supply and demand control device and supply and demand control method
CN106063068A (en) * 2014-11-13 2016-10-26 松下知识产权经营株式会社 Apparatus control device and demand response method
WO2017016617A1 (en) * 2015-07-30 2017-02-02 Siemens Aktiengesellschaft Method for computer-assisted parameterisation of a converter in an electrical grid

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549869A (en) * 2009-07-27 2012-07-04 响应负载有限公司 Responsive load monitoring system and method
CN103527407A (en) * 2012-06-28 2014-01-22 西门子公司 Method and arrangement for responding to a grid event, such as fast frequency drop, by combining demand response, inertial response and spinning reserve
EP2937962A1 (en) * 2012-12-18 2015-10-28 Panasonic Intellectual Property Management Co., Ltd. Supply and demand control device and supply and demand control method
JP2015154595A (en) * 2014-02-14 2015-08-24 アズビル株式会社 Energy management device and energy management method
CN106063068A (en) * 2014-11-13 2016-10-26 松下知识产权经营株式会社 Apparatus control device and demand response method
WO2017016617A1 (en) * 2015-07-30 2017-02-02 Siemens Aktiengesellschaft Method for computer-assisted parameterisation of a converter in an electrical grid

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
电力需求响应技术标准化研究;张晶 等;《中国电机工程学报》;20140805;第34卷(第22期);3623-3629 *

Also Published As

Publication number Publication date
CN107133749A (en) 2017-09-05

Similar Documents

Publication Publication Date Title
CN107133749B (en) Power information physical coupling modeling method considering demand response information
CN102801158B (en) Method for calculating time-lag electric power system eigenvalue and discriminating stability based on Pade approximation
CN110829487B (en) Dynamic frequency prediction method for power system
CN105978725A (en) Non-fragile distributed fault estimation method based on sensor network
CN106354924A (en) FMI-based equipment cooperation simulation system and construction method
CN106452202A (en) Multi-shaft motor synchronous control system and method thereof
CN103104931A (en) Primary frequency modulation method and primary frequency modulation system based on boiler primary wind pressure dynamic compensation
CN109787251A (en) A kind of cluster temperature control Load aggregation model, Parameter identification and counter push away control method
CN102870309A (en) Method and system for facilitating design of a high voltage (HVDC) control system, an HVDC system and a method for optimising an HVDC system
CN106505591B (en) A kind of multi-machine power system ultra-low frequency oscillation analysis method
CN104156886B (en) Method of evaluating flexibility of power supply with renewable energy power system
CN102377190B (en) Grid frequency rate limiting system
CN103942091A (en) MATLAB user-defined model and PSASP joint simulation excitation system simulation method and system
CN110649596B (en) Frequency full-response analytic model considering system initial state
Debbağ et al. Internet based monitoring and control of a wind turbine via PLC
CN106571637B (en) A kind of power system frequency response analysis method considering frequency space division cloth
CN103886419A (en) Power grid operation mode planning, checking and analyzing system
CN109408884A (en) Information processing method for central air-conditioning system Automated Design
CN105406517A (en) Finite time average consistency algorithm-based economic dispatching method for power system
CN104767205B (en) Method for establishing automatic generation control system of electric power system based on wind power plant access
CN109245090B (en) Modeling method of analytic model for predicting frequency lowest point under high-power deficiency
CN105122165B (en) Plant control system and factory control method
CN104035336B (en) Method for designing non-linear controller of MIMO non-minimum phase CSTR
CN106342297B (en) A kind of satellite modeling method with stand-alone device modeling
Makhmalbaf et al. Co-simulation of detailed whole building with the power system to study smart grid applications

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