CN105356454B - A kind of power system load modeling method based on typical load group - Google Patents

A kind of power system load modeling method based on typical load group Download PDF

Info

Publication number
CN105356454B
CN105356454B CN201510746371.9A CN201510746371A CN105356454B CN 105356454 B CN105356454 B CN 105356454B CN 201510746371 A CN201510746371 A CN 201510746371A CN 105356454 B CN105356454 B CN 105356454B
Authority
CN
China
Prior art keywords
load
msub
typical
model
load group
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
CN201510746371.9A
Other languages
Chinese (zh)
Other versions
CN105356454A (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.)
Shanghai University of Electric Power
State Grid Shanghai Electric Power Co Ltd
Original Assignee
Shanghai University of Electric Power
State Grid Shanghai Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power, State Grid Shanghai Electric Power Co Ltd filed Critical Shanghai University of Electric Power
Priority to CN201510746371.9A priority Critical patent/CN105356454B/en
Publication of CN105356454A publication Critical patent/CN105356454A/en
Application granted granted Critical
Publication of CN105356454B publication Critical patent/CN105356454B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The present invention relates to a kind of power system load modeling method based on typical load group, including step：1) typical load group model is established：Typical load group model is generated according to node load information；2) node load model is established：According to typical load group model parameter and the accounting of each typical load group, weighting polymerization obtains node load model parameter, and establishes node load model based on this.Relative to the typical load parameter used in current emulation, the present invention obtains the load parameter of node load model parameter using the weighting polymerization of a variety of typical load groups, therefore node active voltage response characteristic and reactive voltage response characteristic are closer to actual measurement situation.

Description

A kind of power system load modeling method based on typical load group
Technical field
The present invention relates to Power System and its Automation field, more particularly, to a kind of power train based on typical load group System load modeling method.
Background technology
Substantial amounts of calculate shows with result of the test：Load model influences very on the result of calculation of Electrical Power System Dynamic behavior Greatly, also there is certain influence to Load flow calculation, short circuit calculation, safety analysis, voltage stability etc..Under critical condition, also may be used Qualitatively conclusion can fundamentally be changed.
Mechanism dynamic model is applied widely because of its clear concept, strong applicability in Practical Project.China at present Dynamic load model in power networks at different levels uses the motor parallel connection static load of fixed proportion, and wherein motor uses mostly Canonical parameter.There are two classes in the source of the canonical parameter typically used, and one kind is the typical motor that IEEE load operations group is recommended Parameter, another kind of is the recommendation of electrical engineering association of China.In actual applications, may influence to imitate using unified canonical parameter True accuracy：Typical motor only includes limited several classes, is 7 classes as IEEE divides, but actual node load structure may be more complicated, Therefore, typical motor parameter can not reflect the dynamic behaviour of some given area node load well.Can be just so finding Really the asynchronous motor parameter of reflection given area node dynamics is significant.
Conventional modeling method has Component Based and Measurement-based approach.The advantages of Component Based is without the scene of carrying out Actual measurement and experiment, physical model is clear, definite conception, with strong applicability.The common practices of Component Based is to yet to be built at present The network of mould node subordinate, load structure carry out Data Collection, then by work such as load classification, synthesis, finally obtain The load model and parameter of the node, workload is huge and modeling accuracy is general.How the condition of certain model precision is being ensured Lower simplified modeling work amount is practical to load modeling work to be of great importance.
Measurement-based approach has the characteristics of simple and practical, but it is needed in node installation independent of user's statistics Corresponding measurement apparatus, the component parameters for recognizing to obtain can not explain actual conditions well, be not easy to be received by Field Force. Have ready conditions and obtain the node of measuring equipment data, how with reference to two kinds of modeling methods, learn from other's strong points to offset one's weaknesses, establish a kind of both facilitated and grasp Make, the modeling method that and can obtains accurate practical node load model is also a job highly significant.
The content of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide one kind is based on typical load The power system load modeling method of group.
The purpose of the present invention can be achieved through the following technical solutions：
A kind of power system load modeling method based on typical load group, including step：
1) typical load group model is established：Typical load group model is generated according to node load information；
2) node load model is established：It is poly- according to typical load group model parameter and the accounting of each typical load group, weighting Conjunction obtains node load model parameter, and establishes node load model based on this.
The typical load group includes industrial load group, resident load group and Commercial Load group, and step 1) specifically includes Step：
11) node load information is loaded into, divides typical load group, and determines that typical load group is formed；
12) each typical load group is emulated using ETAP, obtains the injecting power voltage responsive characteristic of node；
13) typical load group model is established, and determines its model parameter.
The typical load group model is made up of ZIP static loads submodel and asynchronous motor submodel,
The ZIP static loads submodel is：
Wherein：PsAnd QsThe respectively active power and reactive power of static load, Ps0,Qs0To be respectively that static load has The initial value of work(power and reactive power, U are node voltage, U0For the initial value of node voltage, ap,bp,cp,aq,bq,cqFor allusion quotation The static load parameter of type load group model parameter, and ap+bp+cp=1, aq+bq+cq=1,
The asynchronous motor submodel uses third-order model.
The typical load group model parameter also includes：
Induction-motor load parameter：RsFor stator resistance, XsFor stator reactance, XmFor excitation reactance, XrFor rotor reactance, RrFor Rotor resistance, H motor inertia time constants, A are the accounting of part directly proportional to rotating speed square in mechanical load torque, and B is The accounting of part directly proportional to rotating speed in mechanical load torque,
Sound ratio coefficient kPWith induction-motor load rate coefficient kL。
A in the typical load group model parameterp,bp,cp,aq,bq,cq,Rs,Xm,Xr, A, B use empirical value, Xs,Rr,H, kP,kLObtained using Identification of Genetic Algorithm.
The step 2) specifically includes step：
21) it is loaded into and the accounting of each typical load group and power distribution network equivalence pressure drop；
22) weighting polymerization obtains node load model parameter, is specially：
Wherein：SmN,k=ρkS0kP,k/kL,k, Sm0,k=ρkS0kP,k, S0For node injecting power, kP,kIt is negative for k-th of typical case The sound proportionality coefficient of lotus group model, ρkFor the accounting of k-th of typical load group, kL,kFor the electricity of k-th of typical load group model Engine load rate coefficient, HkFor the motor inertia time constant of k-th of typical load group model, Rs,kFor k-th of typical load The stator resistance of group model, Xs,kFor the stator reactance of k-th of typical load group model, Xm,kFor k-th of typical load group model Excitation reactance, Xr,kFor the rotor reactance of k-th of typical load group model, Rr,kFor rotor resistance, AkFor k-th of typical load The accounting of part directly proportional to rotating speed square, B in the mechanical load torque of group modelkFor in mechanical load torque with rotating speed into just Than the accounting of part, SmN,kFor ρkS0kP,k/kL,k；
23) the node load model parameter based on acquisition establishes node load model.
The accounting of each typical load group and the pressure drop of power distribution network equivalence are obtained by fault oscillograph measurement.
Compared with prior art, the present invention has advantages below：
1) relative to the typical load parameter used in current emulation, it is polymerize using the weighting of a variety of typical load groups and is obtained The load parameter of node load model parameter, therefore node active voltage response characteristic and reactive voltage response characteristic are closer in fact Survey situation.
2) typical load group is designed as three kinds of industrial load group, resident load group and Commercial Load group, can ensure essence While spending, reduce amount of calculation to greatest extent.
3) typical load group model is by ZIP static loads submodel and asynchronous motor submodel, such a modeling pattern skill Art is ripe, simple and convenient.
4) in typical load group model parameter, only the parameter of high sensitivity is recognized, remaining parameter uses representative value, Further modeling efficiency can be improved on the premise of precision is not influenceed.
5) accounting of each typical load group and the pressure drop of power distribution network equivalence are obtained by fault oscillograph measurement, and power voltage rings Answer the goodness of fit of curve and measured curve higher, when having ready conditions acquisition node measured data, it should using measured data to system The node load model that meter obtains is modified.
Brief description of the drawings
Fig. 1 is the key step schematic flow sheet of the present invention；
Fig. 2 is asynchronous induction motor model circuit diagram；
Fig. 3 is improved adaptive GA-IAGA parameter identification steps flow chart schematic diagram；
Fig. 4 is node load precircuit schematic diagram；
Fig. 5 (a) is the typical resident load group power response characteristic in area；
Fig. 5 (b) is regional typical commercial load group power response characteristic；
Fig. 5 (c) is regional typical industry load group power response characteristic；
Fig. 6 is the application method compared to the prior art node load model power response characteristic schematic diagram；
Fig. 7 is situation 1 and the node load model power response characteristic contrast schematic diagram of situation 2 in the present embodiment.
Embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with technical solution of the present invention Premised on implemented, give detailed embodiment and specific operating process, but protection scope of the present invention is not limited to Following embodiments.
A kind of power system load modeling method based on typical load group, including step：
1) typical load group model is established：Typical load group model is generated according to node load information；
In general, regional load group can be divided into industrial load, Commercial Load and resident load and the class of specific load four, Its typical case forms as shown in table 1.Industrial load refers to as Factory Electricity Supply load, wherein more than 50% is induction-motor load, with asynchronous Based on motor.Commercial Load refers to the load for commercial building or business block, and wherein air-conditioning and lighting load accounting is big, air-conditioning Up to 40%-50%.Resident load refers to the load of resident residential area, and lotus is mainly made up of illumination and household electrical appliance, with household electric Based on device load, summer, winter was based on air-conditioning and warmer based on air-conditioning (monomer).Specific load is primarily referred to as actively negative Lotus, such as electric automobile, and impact load, such as smelting and electric railway.Each department can be according to regional load character and modeling Need to divide subgroup, such as industrial load, can further be divided by load character.Specific load typically individually models, the application For general node modeling, therefore do not include specific load in model.
Table 1
Step 1) specifically includes step：
11) node load information is loaded into, divides typical load group, and determines that typical load group is formed；
12) each typical load group is emulated using ETAP, obtains the injecting power voltage responsive characteristic of node；
The typical load group obtained in ETAP softwares to investigation emulates, and the power voltage for obtaining load group node rings Answer characteristic P (U) and Q (U).
13) typical load group model is established, and determines its model parameter.Industrial load group, Commercial Load group and resident bear Lotus group can be expressed as being combined by one group of ZIP static load and asynchronous motor (third-order model).
Static load is represented with ZIP models, as shown in formula (1)：
Wherein：PsAnd QsThe respectively active power and reactive power of static load, Ps0,Qs0To be respectively that static load has The initial value of work(power and reactive power, U are node voltage, U0For the initial value of node voltage, ap,bp,cp,aq,bq,cqFor allusion quotation The static load parameter of type load group model parameter, wherein ap,bp,cpFor active voltage coefficient correlation, aq,bq,cqFor idle electricity Press coefficient correlation, and ap+bp+cp=1, aq+bq+cq=1,
Ps0,Qs0Calculating formula such as formula (2),
Wherein：P0,Q0The respectively total active power and reactive power of node injection, kP' it is motor in node total load Proportion.
Asynchronous motor is formed as shown in Fig. 2 using third-order model, during xy coordinate systems, model can use formula (3)~ (5) represent.
Output equation：
In above-mentioned formula, except ωB=2 π fBOutside=100 π, time t and inertia time constant H, take specified with motor Perunit value on the basis of capacity：
Wherein：Sm0=S0×kPMotor injecting power when being stable state, kLFor rate of load condensate.
To sum up, the parameter in typical load group model includes：
Static load parameter：ap,bp,cpAnd aq,bq,cq,
Induction-motor load parameter：Rs,Xs,Xm,Xr,Rr,H,A,B
Sound ratio coefficient kPWith induction-motor load rate coefficient kL
In above-mentioned parameter, the respective implication of induction-motor load parameter is respectively：RsFor stator resistance, XsFor stator reactance, Xm For excitation reactance, XrFor rotor reactance, RrFor rotor resistance, H motor inertia time constants, A is with turning in mechanical load torque The accounting of fast square directly proportional part, B are the accounting of part directly proportional to rotating speed in mechanical load torque.
In transient state load model, static load typically can be represented directly with constant-impedance, be calculated using steady-state value, i.e., Formula (6) can be used to solve：
To equivalent induction conductivity, by carrying out sensitivity analysis to parameter, stator reactance, rotor resistance are electronic Machine inertia time constant sensitivity is higher.In addition, the influence of sound proportionality coefficient and induction-motor load rate coefficient to transient stability It is more significant.Therefore, only the parameter of high sensitivity is recognized in the application, remaining parameter is using representative value (in the present embodiment Using IEEE typical motors).To sum up, the parameter that simplify needs to recognize in identification afterload group model is 5, is specially：Xs, Rr,H,kP,kL, its identification process carries out typical load swarm parameter identification using improved adaptive GA-IAGA, as shown in figure 3, due to specific Identification details belong to prior art, therefore the application is not reinflated.
2) node load model is established：It is poly- according to typical load group model parameter and the accounting of each typical load group, weighting Conjunction obtains node load model parameter, and establishes node load model based on this.
Load configuration that the application interior joint load model uses considers to match somebody with somebody as shown in figure 4, referred to as SLM models in model Network equivalence reactance, equivalent static load and equivalent motor.Equivalent motor uses third-order model, and specific pattern is referring to formula (3)~(5), equivalent static load uses constant-impedance model, as shown in formula (6).Model parameter includes motor parameter (8 It is individual), and sound ratio coefficient kPWith induction-motor load rate coefficient kL(2) outside totally 10 parameters, plus power distribution network equivalent reactance XD。
Step 2) specifically includes step：
21) it is loaded into and the accounting of each typical load group and power distribution network equivalence pressure drop, loading shares two kinds of situations.
Situation 1：Each typical load group ratio only in known node load, the accounting of each typical load group are obtained by investigating, Design load or estimated value are selected in the pressure drop of power distribution network equivalence；
Situation 2：The accounting of each typical load group and the pressure drop of power distribution network equivalence are obtained by fault oscillograph measurement.
22) weighting polymerization obtains node load model parameter, for situation 1：
A) power distribution network equivalent impedance
Estimated according to the capacity of the transformer station and line end voltage condition：
Known UiAnd UL, then Δ U=Ui-UL(when node voltage situation is unknown, then it can be passed on the line according to code, power It is 2%~5% that transporting, which causes terminal voltage landing limit value, i.e. Δ U=(2%~5%) UN, UL=Ui-ΔU)
(as Δ U=(2%~5%) UNWhen,
B) equivalent motor parameter, sound proportionality coefficient and induction-motor load rate coefficient
Taking each typical load group, proportion coefficient is in node：
Equivalent motor parameter, sound proportionality coefficient and induction-motor load rate coefficient are obtained using weighting polymerization methodses, are gathered Conjunction method is described below：
Node injecting power under known limitIt can estimate that load direct access point is steady according to formula (7)~(9) State injecting powerFor statement it is convenient for the sake of, below in text withRepresent, respective formula in each parameter aggregation method such as table 2 It is shown：
Table 2
※：If parameter is all 0, equivalent parameters take 0 and (for example work as B1=B2=B3When=0, B=0)
C) equivalent static load
The parameter k obtained using polymerizationP, substituting into formula (6) to ask for.
For situation 2, in such cases, by power distribution network equivalent reactance, each typical load group accounting coefficient of node is used as and needed The parameter to be recognized, identification use improved adaptive GA-IAGA, parameter [X to be identifiedD,ρ1,ρ2], calculation process and the parameter shown in Fig. 3 It is similar to recognize flow.After power distribution network equivalent reactance and each typical load group accounting coefficient is obtained, calculated using formula shown in table 2 Load model parameters.
Specifically, obtaining node load model parameter and being：
Wherein：SmN,k=ρkS0kP,k/kL,k, Sm0,k=ρkS0kP,k, S0For node injecting power, kP,kIt is negative for k-th of typical case The sound proportionality coefficient of lotus group model, ρkFor the accounting of k-th of typical load group, kL,kFor the electricity of k-th of typical load group model Engine load rate coefficient, HkFor the motor inertia time constant of k-th of typical load group model, Rs,kFor k-th of typical load The stator resistance of group model, Xs,kFor the stator reactance of k-th of typical load group model, Xm,kFor k-th of typical load group model Excitation reactance, Xr,kFor the rotor reactance of k-th of typical load group model, Rr,kFor rotor resistance, AkFor k-th of typical load The accounting of part directly proportional to rotating speed square, B in the mechanical load torque of group modelkFor in mechanical load torque with rotating speed into just Than the accounting of part, SmN,kFor ρkS0kP,k/kL,k；
23) the node load model parameter based on acquisition establishes node load model.
Application example：
The inspection for the load modeling method that the application is illustrated is carried out by taking Anshun Prefecture's 110kV nodes as an example.
11) node load information is loaded into, divides typical load group, and determines that typical load group is formed；
First, investigate 110kV node subordinate's load structures, formed typical industry load group, typical commercial load group and Typical resident load group, the specific composition of all kinds of typical load groups are following as shown in 3~table of table 5：
The typical industry load group of table 3 is formed
The typical commercial load group of table 4 is formed
Load type Commercial air-conditioner (%) Illuminate (%) Heat-pump-type heater (%) Other equipment (%)
Load accounting (%) 56.5 26.1 8.7 8.7
The typical resident load group of table 5 is formed
Load type Air-conditioning (%) Heater Illumination Refrigerator Other
Load accounting (%) 20.4 13.5 16.4 28.6 21.1
12) typical load group is emulated using ETAP, obtains the injecting power voltage responsive characteristic of load bus, point Not as Fig. 5 (a)~5 (c) are shown.
13) typical load group model is established, and determines its model parameter.Letter is carried out to typical load group using genetic algorithm Change parameter identification, obtain each typical load swarm parameter.Identification result is as shown in table 6.
Table 6
2) node load model is established
Situation 1：Each typical load group proportionality coefficient of known node.
Investigation knows, the transformer station shares two 110kV transformers, allows voltage landing to be 4% during nominal load；Load Based on industrial load, comprising a part of resident load, Commercial Load is smaller, and its accounting is：Industrial load 60%, resident load 30%, Commercial Load 10%. calculates the node synthetic load parameter obtained as shown in situation 1 in table 7.
Situation 2：X is obtained by node failure oscillograph data identificationD,ρ1,ρ2, the weighting polymerization of equivalent motor parameter obtains , obtained node synthetic load parameter is as shown in situation 2 in table 7.
Table 7
In figure 6 to actual measurement, the typical load parameter used at present, the work(for the load parameter (situation 2) that the application obtains Rate voltage characteristic compares.It can be seen that relative to the typical load parameter used in current emulation, the application obtains Load parameter node active voltage response characteristic and reactive voltage response characteristic closer to actual measurement situation, illustrate in the application It is proposed it is preferable using the modeling method effect of typical load group.
Fig. 7 lists the actual active reactive response curve that fault oscillograph captures, situation 1 in the present invention is respectively adopted With the power voltage resonse characteristic of the node load model of the parameter of situation 2.
From figure 7, it is seen that situation 1 (only polymerize with investigational data and obtain load model parameters) and situation 2 are (using actual measurement number Obtain load model parameters according to polymerizeing again after being recognized) under the node load model parameter that obtains preferably describe node The power voltage response characteristic of load, active similar with measured curve shape with idle dynamic characteristic, active power is coincide Degree is better than reactive power, and error mainly appears on voltage descending branch.Two groups of parameters are compared, the power voltage response curve of situation 2 It is higher with the goodness of fit of measured curve, illustrate when having ready conditions acquisition node measured data, it should using measured data to statistics The node load model of acquisition is modified.

Claims (2)

  1. A kind of 1. power system load modeling method based on typical load group, it is characterised in that including step：
    1) typical load group model is established：Typical load group model is generated according to node load information,
    2) node load model is established：According to typical load group model parameter and the accounting of each typical load group, weighting polymerize Node load model is established to node load model parameter, and based on this；
    The typical load group includes industrial load group, resident load group and Commercial Load group, and step 1) specifically includes step：
    11) node load information is loaded into, divides typical load group, and determines that typical load group is formed,
    12) each typical load group is emulated using ETAP, obtains the injecting power voltage responsive characteristic of node,
    13) typical load group model is established, and determines its model parameter；
    The typical load group model is made up of ZIP static loads submodel and asynchronous motor submodel,
    The ZIP static loads submodel is：
    <mrow> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>P</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>p</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>p</mi> </msub> <mo>(</mo> <mfrac> <mi>U</mi> <msub> <mi>U</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>c</mi> <mi>p</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>U</mi> <msub> <mi>U</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow>
    <mrow> <msub> <mi>Q</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>Q</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>q</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mi>q</mi> </msub> <mo>(</mo> <mfrac> <mi>U</mi> <msub> <mi>U</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>c</mi> <mi>q</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>U</mi> <msub> <mi>U</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow>
    Wherein：PsAnd QsThe respectively active power and reactive power of static load, Ps0,Qs0To be respectively static load wattful power The initial value of rate and reactive power, U are node voltage, U0For the initial value of node voltage, ap,bp,cp,aq,bq,cqIt is negative for typical case The static load parameter of lotus group model parameter, and ap+bp+cp=1, aq+bq+cq=1,
    The asynchronous motor submodel uses third-order model；
    The typical load group model parameter also includes：
    Induction-motor load parameter：RsFor stator resistance, XsFor stator reactance, XmFor excitation reactance, XrFor rotor reactance, RrFor rotor Resistance, H are motor inertia time constant, and A is the accounting of part directly proportional to rotating speed square in mechanical load torque, and B is machine The accounting of part directly proportional to rotating speed in tool load torque,
    Sound ratio coefficient kPWith induction-motor load rate coefficient kL,
    And a in the typical load group model parameterp,bp,cp,aq,bq,cq,Rs,Xm,Xr, A, B use empirical value, Xs,Rr,H, kP,kLObtained using Identification of Genetic Algorithm；
    Step 2) specifically includes step：
    21) accounting and the power distribution network equivalence pressure drop of each typical load group are loaded into,
    22) weighting polymerization obtains node load model parameter, is specially：
    Wherein：SmN,k=ρkS0kP,k/kL,k, Sm0,k=ρkS0kP,k, S0For node injecting power, kP,kFor k-th of typical load group's mould The sound proportionality coefficient of type, ρkFor the accounting of k-th of typical load group, kL,kMotor for k-th of typical load group model is born Lotus rate coefficient, HkFor the motor inertia time constant of k-th of typical load group model, Rs,kFor k-th of typical load group model Stator resistance, Xs,kFor the stator reactance of k-th of typical load group model, Xm,kFor the excitation of k-th of typical load group model Reactance, Xr,kFor the rotor reactance of k-th of typical load group model, Rr,kFor rotor resistance, AkFor k-th of typical load group model Mechanical load torque in part directly proportional to rotating speed square accounting, BkFor part directly proportional to rotating speed in mechanical load torque Accounting,
    23) the node load model parameter based on acquisition establishes node load model.
  2. A kind of 2. power system load modeling method based on typical load group according to claim 1, it is characterised in that The accounting of each typical load group and the pressure drop of power distribution network equivalence are obtained by fault oscillograph measurement.
CN201510746371.9A 2015-11-05 2015-11-05 A kind of power system load modeling method based on typical load group Active CN105356454B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510746371.9A CN105356454B (en) 2015-11-05 2015-11-05 A kind of power system load modeling method based on typical load group

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510746371.9A CN105356454B (en) 2015-11-05 2015-11-05 A kind of power system load modeling method based on typical load group

Publications (2)

Publication Number Publication Date
CN105356454A CN105356454A (en) 2016-02-24
CN105356454B true CN105356454B (en) 2017-11-21

Family

ID=55332363

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510746371.9A Active CN105356454B (en) 2015-11-05 2015-11-05 A kind of power system load modeling method based on typical load group

Country Status (1)

Country Link
CN (1) CN105356454B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018144009A1 (en) * 2017-02-03 2018-08-09 Duke Energy Corporation Power management methods for a circuit of a substation, and related apparatuses and computer program products
CN106921159A (en) * 2017-02-21 2017-07-04 国网山东省电力公司电力科学研究院 A kind of active power modulator approach of many air conditioner load groups
CN106972497B (en) * 2017-04-24 2019-12-10 国网山东省电力公司烟台供电公司 probability model establishing and solving method for air conditioner load group dispersion modulation strategy
CN107069717B (en) * 2017-05-11 2020-05-15 国家电网公司 Power distribution network load modeling method
CN107147112B (en) * 2017-05-24 2019-09-24 河海大学 A kind of power system modeling method and system based on weighting polymerization
CN109245100B (en) * 2018-11-07 2020-09-01 国网浙江省电力有限公司经济技术研究院 Load dynamic modeling method considering time-varying property of load composition of alternating-current and direct-current power distribution network
CN109687438B (en) * 2018-12-19 2022-06-24 西南交通大学 Power grid fragile line identification method considering high-speed rail impact load effect
CN109764539A (en) * 2019-01-04 2019-05-17 北京华建网源电力设计研究院有限公司 Dynamic system and control model for load group of electric water heater
CN110472339B (en) * 2019-08-16 2023-05-16 东莞理工学院 Room lighting load calculation method based on Monte Carlo
CN111415059B (en) * 2020-01-08 2022-04-29 天津大学 Practical model machine construction and online application method
CN111796143B (en) * 2020-09-10 2020-12-15 深圳华工能源技术有限公司 Energy-saving metering method for energy-saving equipment of power distribution and utilization system
CN112290538A (en) * 2020-09-30 2021-01-29 天津大学 Load model parameter online correction method based on aggregation-identification double-layer framework

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101777765A (en) * 2010-01-27 2010-07-14 中国电力科学研究院 On-line load simulation method of power system
CN101789598A (en) * 2010-03-05 2010-07-28 湖北省电力试验研究院 Power system load modelling method
CN104992017A (en) * 2015-07-01 2015-10-21 河海大学 Wind power random fluctuation based generalized load modeling method and apparatus

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002034150A (en) * 2000-07-14 2002-01-31 Tokyo Electric Power Co Inc:The Synthesizing method for load model in power system simulator

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101777765A (en) * 2010-01-27 2010-07-14 中国电力科学研究院 On-line load simulation method of power system
CN101789598A (en) * 2010-03-05 2010-07-28 湖北省电力试验研究院 Power system load modelling method
CN104992017A (en) * 2015-07-01 2015-10-21 河海大学 Wind power random fluctuation based generalized load modeling method and apparatus

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
针对电压慢恢复现象的电机负荷建模方法;华嘉成等;《电力自动化设备》;20150930;第35卷(第9期);第94页-第98页 *

Also Published As

Publication number Publication date
CN105356454A (en) 2016-02-24

Similar Documents

Publication Publication Date Title
CN105356454B (en) A kind of power system load modeling method based on typical load group
CN105071381B (en) State enumeration reliability evaluation method and device based on influence increment
Morison et al. Practical issues in load modeling for voltage stability studies
CN107451364B (en) Method for identifying equivalent parameters of DFIG wind power plant
CN103412994B (en) A kind of method determining high-speed wind tunnel large aircraft model contracting ratio
CN104062618B (en) A kind of capacitive apparatus on-Line Monitor Device method of calibration of dual signal source
CN103033716B (en) Calculation method of proportion of each lode component in power grid comprehensive load model
CN109004652A (en) A kind of direct current multi-infeed systems receiving end grid strength preparation method containing asynchronous machine
CN106443246A (en) Online identification method of small-disturbance stability parameters based on PMU measurement data
CN108614190B (en) Power distribution network T-connection line parameter estimation method based on mu PMU and AMI
CN103793556B (en) Comprehensively go out transformer station&#39;s load modeling method of specific electric load
CN103595040A (en) Power distribution network comprising direct-current distributed power supply and analysis method of comprehensive load characteristics of power distribution network
CN109116156A (en) A kind of method and apparatus that transmission line of electricity line loss is determined based on mutual inductor output signal
CN106207284A (en) A kind of method promoting ferric phosphate lithium cell partial volume Capacity uniformity
CN103441493B (en) Method for automatically selecting key sections on load side of electrical partition of power grid
CN112907075A (en) Method for identifying parameters of comprehensive load model of power system
CN108254712A (en) A kind of electric energy meter electric source loop protection Simulation Analysis method and device
CN103632031B (en) A kind of rural area based on load curve decomposition load type load modeling method
CN103514571A (en) Load curve decomposition based load modeling method for commercial load and resident load
CN102508086A (en) Parameter estimation method for electric power circuit in T-type connection
CN107516115B (en) A kind of load model canonical parameter extracting method based on k- central point algorithm
CN104701839A (en) Air conditioner load modeling method based on least squares parameter identification
CN114970154A (en) Transient voltage supporting capability quantitative evaluation index construction method
Ma et al. Impact analysis of wind generation on voltage stability and system load margin
CN105866585A (en) Method for harmonic source identification and liability sharing based on distortion power

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant