CN112464436A - Step length adjusting method for parameter identification of power simulation model - Google Patents

Step length adjusting method for parameter identification of power simulation model Download PDF

Info

Publication number
CN112464436A
CN112464436A CN202011190349.8A CN202011190349A CN112464436A CN 112464436 A CN112464436 A CN 112464436A CN 202011190349 A CN202011190349 A CN 202011190349A CN 112464436 A CN112464436 A CN 112464436A
Authority
CN
China
Prior art keywords
value
simulation
delta
parameter
simulation model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011190349.8A
Other languages
Chinese (zh)
Other versions
CN112464436B (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.)
Electric Power Research Institute of Hainan Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Hainan Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electric Power Research Institute of Hainan Power Grid Co Ltd filed Critical Electric Power Research Institute of Hainan Power Grid Co Ltd
Priority to CN202011190349.8A priority Critical patent/CN112464436B/en
Publication of CN112464436A publication Critical patent/CN112464436A/en
Application granted granted Critical
Publication of CN112464436B publication Critical patent/CN112464436B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a step length adjusting method for parameter identification of an electric power simulation model, which comprises the following steps: s1, establishing a power system simulation model, and setting an initial parameter value of any index to be measured in the simulation model
Figure DDA0002752597430000011
And sets a simulation threshold deltaopt(ii) a S2, establishing a parameter identification submodel: s3, setting an optimal solution of f (T), calculating control coefficients a and b based on the optimal solution, and obtaining a step function expression; and S4, iteratively selecting an adjusting coefficient delta, and iteratively adjusting the step function expression continuously until the result of the parameter identification submodel meets the requirement.

Description

Step length adjusting method for parameter identification of power simulation model
Technical Field
The invention relates to the technical field of power simulation, in particular to a step length adjusting method for parameter identification of a power simulation model.
Background
With the increasing development of power systems, the types of equipment parameters and operation data used for calculation and analysis of the power systems are more and more, the data volume is also increased continuously, and accurate equipment parameters are the basis for protection setting, fault analysis, load flow calculation and loss analysis of the power systems, and have very important significance for safe and stable operation of the power systems. Because the required parameters are too large to be combed, all the parameters do not have actual measurement data, and the actual measurement results of the parameters have deviation from the actual measurement results under the influence of factors such as a test method, environmental conditions, operation conditions and the like. Therefore, the method for researching the simulation parameter identification of the power system is used for identifying and checking parameters required by simulation, provides simulation accuracy and plays an important role in improving the safety and stability of the power grid.
Because the simulation calculation process of the power system is a solving process of a nonlinear equation, the related parameter quantity is huge, mutual coupling relations exist among parameters, and independent identification for each parameter is difficult, a model-free identification method is needed to identify the parameters of the power system.
Disclosure of Invention
The present invention is directed to a step size adjustment method for parameter identification of an electrical simulation model, so as to solve the problems in the background art.
The invention is realized by the following technical scheme: a step size adjusting method for parameter identification of a power simulation model comprises the following steps:
s1, establishing a power system simulation model, and setting an initial parameter value of any index to be measured in the simulation model
Figure BDA0002752597410000021
And sets a simulation threshold deltaopt；
S2, establishing a parameter identification submodel:
Figure BDA0002752597410000022
wherein, C(j+1)For the parameter values after iteration, μ is the step size, Δ C(j)For parametric disturbance values, Δ δiFor simulating the difference, a and b are control coefficients, T is an iteration progress coefficient, and T is deltaoptDelta, delta is an adjusting coefficient;
s3, setting an optimal solution of f (T), calculating control coefficients a and b based on the optimal solution, and obtaining a step function expression;
and S4, iteratively selecting an adjusting coefficient delta, and iteratively adjusting the step function expression continuously until the result of the parameter identification submodel meets the requirement.
Preferably, the indexes to be measured of the line include resistance, reactance and admittance of the line.
Preferably, the theoretical value or the first measured value of the line index is set as an initial parameter value
Figure BDA0002752597410000023
Preferably, any normal number is set as the parameter disturbance value Δ C(j)And solving the simulation difference deltai：
Calculating the perturbed parameter value
Figure BDA0002752597410000024
The perturbed parameter value is used for
Figure BDA0002752597410000025
Putting the simulation model into the simulation model to obtain a first simulation value
Figure BDA0002752597410000026
The perturbed parameter value is used for
Figure BDA0002752597410000027
Putting the simulation model into the simulation model to obtain a second simulation value
Figure BDA0002752597410000028
The first simulation value
Figure BDA0002752597410000029
And a second simulation value
Figure BDA00027525974100000210
The difference is the simulation difference deltai。
Preferably, when T is 0.683, there is an optimal solution for f (T).
Preferably, the control coefficient a is 1.47 to 4.63 f (0.683) calculated by the following formula, and the control coefficient b is 4.63 f (0.683) to 2.47 calculated by the following formula.
Preferably, the parameter identification submodel is solved to obtain the iterated parameter value C after iteration(j+1)The iterated parameter value C(j+1)Putting the simulated model into the simulation model for calculation to obtain a third simulation value delta(j+1)；
If delta(j+1)When the value is more than or equal to delta, the third simulation value delta(j+1)If the parameter is not in accordance with the requirement, the parameter disturbance value delta C is generated iteratively(j)And adjusting the coefficient delta, and solving the parameter identification submodel again.
Compared with the prior art, the invention has the following beneficial effects: according to the step length adjusting method for parameter identification of the power simulation model, provided by the invention, through self-adaptive adjustment of the step length, the convergence and convergence speed of parameter identification can be effectively improved, so that the parameter identification of a power system can be realized more quickly, the practicability of the model parameter identification of the wide area power system is further improved, and positive effects are generated on improving the simulation analysis precision of the power system and correctly making a power grid construction plan and an operation mode. In addition, the invention can be suitable for various power system simulation software used in the power industry of China at present, thereby having good popularization and application prospects.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only preferred embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without inventive exercise.
FIG. 1 is a flowchart of a step size adjustment method for parameter identification of a power simulation model according to the present invention.
Detailed Description
In order to better understand the technical content of the invention, specific embodiments are provided below, and the invention is further described with reference to the accompanying drawings.
In the parameter identification method of the power system, the convergence and the convergence speed of the parameter identification are greatly dependent on the step length coefficient, the larger the step length coefficient is, the faster the convergence speed is, but the poorer the convergence is, the smaller the step length is, the smaller the convergence speed is, the better the convergence is, but the parameter identification method is easy to fall into a local extreme value. In the initial stage of algorithm iteration, the algorithm convergence speed can be improved due to the large step length, the situation that the algorithm falls into a local extreme value is avoided, and a better parameter is found more quickly. In the middle and later stages of algorithm iteration, a smaller step length is adopted, the fineness of searching extreme values is improved, and the algorithm convergence is improved.
In order to improve the convergence and the convergence speed of the method, referring to fig. 1, the invention discloses a step length adjustment method for parameter identification of a power simulation model, which comprises the following steps:
s1, establishing a power system simulation model, and setting an initial parameter value of any index to be measured in the simulation model
Figure BDA0002752597410000041
And sets a simulation threshold deltaopt；
Specifically, the line index to be measured includes resistance, reactance and admittance of the line, and in the embodiment of the present invention, a theoretical value or a first measurement value of the line index is set as an initial parameter value
Figure BDA0002752597410000042
S2, establishing a parameter identification submodel:
Figure BDA0002752597410000043
wherein, C(j+1)For the parameter values after iteration, μ is the step size, Δ C(j)For parametric disturbance values, Δ δiFor simulation differenceThe values a and b are control coefficients, T is an iteration progress coefficient, and T is deltaoptDelta, delta is an adjusting coefficient;
parameter disturbance value Δ C(j)Calculating the value of the parameter after disturbance for any normal number
Figure BDA0002752597410000044
Figure BDA0002752597410000045
The perturbed parameter value is used for
Figure BDA0002752597410000046
Putting the simulation model into the simulation model to obtain a first simulation value
Figure BDA0002752597410000047
The perturbed parameter value is used for
Figure BDA0002752597410000048
Putting the simulation model into the simulation model to obtain a second simulation value
Figure BDA0002752597410000049
The first simulation value
Figure BDA00027525974100000410
And a second simulation value
Figure BDA00027525974100000411
The difference is the simulation difference deltai
S3, setting an optimal solution of f (T), calculating control coefficients a and b based on the optimal solution, and obtaining a step function expression;
in the present embodiment, when T is 0.683, f (T) has an optimum solution, the control coefficient a is calculated by the following equation to be 1.47 to 4.63 f (0.683), the control coefficient b is calculated by the following equation to be 4.63 f (0.683) to 2.47, and the function expression of the step μ can be obtained based on the calculated control coefficients a and b.
And S4, iteratively selecting an adjusting coefficient delta, and iteratively adjusting the step function expression continuously until the result of the parameter identification submodel meets the requirement.
Based on parameter disturbance value Delta C(j)Function expression of step size mu, initial parameter value
Figure BDA0002752597410000051
Simulated difference deltaiThe parameter identification submodel may be solved,
preferably, any normal number is set as the parameter disturbance value Δ C(j)And solving the simulation difference deltaiObtaining the iterated parameter value C after iteration(j+1)The iterated parameter value C(j+1)Putting the simulated model into the simulation model for calculation to obtain a third simulation value delta(j+1)；
If delta(j+1)When the value is more than or equal to delta, the third simulation value delta(j+1)If the parameter is not in accordance with the requirement, the parameter disturbance value delta C is generated iteratively(j)And adjusting the coefficient delta, and solving the parameter identification submodel again when the coefficient delta is equal tooptWhen δ, the iteration ends.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (7)

1. A step length adjusting method for parameter identification of a power simulation model is characterized by comprising the following steps:
s1, establishing a power system simulation model, and setting an initial parameter value of any index to be measured in the simulation model
Figure FDA0002752597400000011
And sets a simulation threshold deltaopt；
S2, establishing a parameter identification submodel:
Figure FDA0002752597400000012
wherein, C(j+1)For the parameter values after iteration, μ is the step size, Δ C(j)For parametric disturbance values, Δ δiFor simulating the difference, a and b are control coefficients, T is an iteration progress coefficient, and T is deltaoptDelta, delta is an adjusting coefficient;
s3, setting an optimal solution of f (T), calculating control coefficients a and b based on the optimal solution, and obtaining a step function expression;
and S4, iteratively selecting an adjusting coefficient delta, and iteratively adjusting the step function expression continuously until the result of the parameter identification submodel meets the requirement.
2. The step size adjustment method for parameter identification of the power simulation model according to claim 1, wherein the line to-be-tested indicators include resistance, reactance, and admittance of the line.
3. The method as claimed in claim 2, wherein the theoretical value or the first measured value of the line index is set as the initial parameter value
Figure FDA0002752597400000018
4. The method as claimed in claim 2, wherein any normal number is set as the parameter disturbance value Δ C(j)And solving the simulation difference deltai：
Calculating the perturbed parameter value
Figure FDA0002752597400000013
The perturbed parameter value is used for
Figure FDA0002752597400000014
Putting the simulation model into the simulation model to obtain a first simulation value
Figure FDA0002752597400000015
The perturbed parameter value is used for
Figure FDA0002752597400000016
Putting the simulation model into the simulation model to obtain a second simulation value
Figure FDA0002752597400000017
The first simulation value
Figure FDA0002752597400000021
And a second simulation value
Figure FDA0002752597400000022
The difference is the simulation difference deltai。
5. The method of claim 4, wherein when T is 0.683, f (T) has an optimal solution.
6. The method of claim 5, wherein the control coefficient a is calculated as 1.47-4.63 f (0.683), and the control coefficient b is calculated as 4.63 f (0.683) -2.47.
7. The method of claim 6, wherein the sub-model is solved to obtain the iterated parameter value C(j+1)The iterated parameter value C(j +1)Putting the simulated model into the simulation model for calculation to obtain a third simulation value delta(j+1)；
If delta(j+1)When the value is more than or equal to delta, the third simulation value delta(j+1)If not, generating parameters iterativelyNumber of perturbation values Δ C(j)And adjusting the coefficient delta, and solving the parameter identification submodel again.
CN202011190349.8A 2020-10-30 2020-10-30 Step length adjusting method for parameter identification of power simulation model Active CN112464436B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011190349.8A CN112464436B (en) 2020-10-30 2020-10-30 Step length adjusting method for parameter identification of power simulation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011190349.8A CN112464436B (en) 2020-10-30 2020-10-30 Step length adjusting method for parameter identification of power simulation model

Publications (2)

Publication Number Publication Date
CN112464436A true CN112464436A (en) 2021-03-09
CN112464436B CN112464436B (en) 2022-11-08

Family

ID=74835715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011190349.8A Active CN112464436B (en) 2020-10-30 2020-10-30 Step length adjusting method for parameter identification of power simulation model

Country Status (1)

Country Link
CN (1) CN112464436B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255268A (en) * 2021-05-21 2021-08-13 北京华大九天科技股份有限公司 Method for detecting and repairing transient analysis non-convergence in circuit simulation

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200293627A1 (en) * 2019-03-13 2020-09-17 General Electric Company Method and apparatus for composite load calibration for a power system

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200293627A1 (en) * 2019-03-13 2020-09-17 General Electric Company Method and apparatus for composite load calibration for a power system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘道伟等: "基于轨迹灵敏度的戴维南等效参数迭代优化辨识", 《中国电机工程学报》 *
梁钰等: "基于动态模式分解的电力***主导振荡模式及参数识别方法研究", 《电子设计工程》 *
申中一等: "基于变步长自适应线性神经网络的PMSM参数辨识", 《电子测量技术》 *
赵菲等: "基于WAMS的输电线路参数在线辨识的研究", 《电力科学与工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113255268A (en) * 2021-05-21 2021-08-13 北京华大九天科技股份有限公司 Method for detecting and repairing transient analysis non-convergence in circuit simulation
CN113255268B (en) * 2021-05-21 2022-05-24 北京华大九天科技股份有限公司 Method for detecting and repairing transient analysis non-convergence in circuit simulation

Also Published As

Publication number Publication date
CN112464436B (en) 2022-11-08

Similar Documents

Publication Publication Date Title
CN107016489A (en) A kind of electric power system robust state estimation method and device
CN103617816B (en) The measuring method of reactor core power distribution
CN111224404B (en) Power flow rapid control method for electric power system with controllable phase shifter
CN108336739A (en) A kind of Probabilistic Load Flow on-line calculation method based on RBF neural
CN112288326A (en) Fault scene set reduction method suitable for toughness evaluation of power transmission system
CN112464436B (en) Step length adjusting method for parameter identification of power simulation model
CN114583767B (en) Data-driven wind power plant frequency modulation response characteristic modeling method and system
CN117748507A (en) Distribution network harmonic access uncertainty assessment method based on Gaussian regression model
CN111814284A (en) On-line voltage stability evaluation method based on correlation detection and improved random forest
CN106532712A (en) Rectangular coordinate Newton method load flow calculation method for small-impedance-branch-containing power grid based on compensation method
CN107832959B (en) Voltage stability evaluation method considering load characteristics and power supply constraints
CN112464437B (en) Parameter identification method of electric power simulation model
CN105701265A (en) Double-fed wind generator modeling method and apparatus
CN110943473A (en) Generator coherence identification method based on wide area measurement system and clustering theory
CN111175608A (en) Power distribution network harmonic responsibility quantitative division method based on accelerated independent component analysis
CN110568260A (en) Power transmission line harmonic parameter estimation method for power grid harmonic analysis
CN116400266A (en) Transformer fault detection method, device and medium based on digital twin model
CN115542236A (en) Method and device for estimating running error of electric energy meter
CN108804843A (en) A kind of cutting load execution station emulation mode, apparatus and system
CN114744631A (en) Data driving voltage estimation method based on non-PMU power distribution network
CN113794198A (en) Method, device, terminal and storage medium for suppressing broadband oscillation
CN113946973A (en) Power supply reliability related index analysis method based on grey correlation algorithm
CN113537821A (en) Rapid assessment method and system for state of power system
CN113537338A (en) Robust line parameter identification method based on LSTM neural network and improved SCADA data
Fan et al. Dielectric loss angle data processing based on adaptive weighted data fusion algorithm of the aging mine cable

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