CN103632031B - A kind of rural area based on load curve decomposition load type load modeling method - Google Patents

A kind of rural area based on load curve decomposition load type load modeling method Download PDF

Info

Publication number
CN103632031B
CN103632031B CN201310500676.2A CN201310500676A CN103632031B CN 103632031 B CN103632031 B CN 103632031B CN 201310500676 A CN201310500676 A CN 201310500676A CN 103632031 B CN103632031 B CN 103632031B
Authority
CN
China
Prior art keywords
load
curve
calculate
month
sigma
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
CN201310500676.2A
Other languages
Chinese (zh)
Other versions
CN103632031A (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
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
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, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201310500676.2A priority Critical patent/CN103632031B/en
Publication of CN103632031A publication Critical patent/CN103632031A/en
Application granted granted Critical
Publication of CN103632031B publication Critical patent/CN103632031B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention provides a kind of rural area based on load curve decomposition load type load modeling method, said method comprising the steps of: decomposes rural area load curve, and calculates the ratio of total amount shared by every kind of load in the load of rural area;Calculating static load equivalent parameters, dynamic load equivalent parameters and distribution network system impedance, and combine the power supply area network topology data of rural area load bus, obtain rural area load type load model.The present invention provides a kind of rural area based on load curve decomposition load type load modeling method, the shortcoming the method overcoming conventional statistics synthesis, can be resident load Site Modeling quickly, easily and accurately, the load model generated by the present invention can improve the degree of accuracy that grid simulation calculates, and ensures power grid security, reliably, run economically.

Description

A kind of rural area based on load curve decomposition load type load modeling method
Technical field
The present invention relates to a kind of modeling method, a kind of rural area based on load curve decomposition load type load modeling Method.
Background technology
Power system digital simulation oneself become the main tool of planning and design of power system, management and running and analysis and research, power train The Mathematical Modeling of each element of uniting and the total system Mathematical Modeling being made up of it are the bases of power system digital simulation, the standard of model Whether true directly affect simulation result and decision scheme based on this.Model and parameter that emulation is used are accuracy of simulation Important decisive factor, generation current machine, excitation system, governing system, transformer, the detail mathematic model of transmission line of electricity and Modeling technique has been obtained for well developing, and electric load model is the simplest comparatively speaking, often from basic physical concept The practical model set out and use and parameter.For many years, each bulk power grid of China when Electrical power system analysis and computing, generally according to Experience selectes certain common load model (such as motor+constant-impedance model or firm power+constant impedance model) qualitatively Determine model parameter.
Power system load modeling, due to factors such as its complexity, distributivity, time variation and randomnesss, determines its mathematical modulo The difficulty that type is set up.Current load modeling method can be basically divided into 3 classes, i.e. Component Based, totally distinguish survey method and fault Fitting process.These three method is respectively arranged with its pluses and minuses.The load model obtained with Component Based has clear physics conception, is prone to The advantage understood by engineering staff, but its core is built upon on the basis of " statistics is complete, and part throttle characteristics is accurate ", This point is often difficult to, and can not often add up, thus cannot consider the time dependent characteristic of load.Always Body examination distinguishes that method avoids substantial amounts of statistical work, it is possible to obtain time dependent online Real-time Load characteristic, asking of its maximum Topic is excessively to rely on disturbance accident, and the model parameter physical significance of identification is indefinite, and another problem is to be difficult in systems All transformer stations are fitted with closing device.The advantage of fault fitting method is parameter determination process and Selection parameter during the calculating of present program Process consistent, and reproduction can be obtained under some fault.But actually it is a kind of method that examination is gathered, under some fault Load parameter whether be applicable to other faults and be difficult to ensure that, and think that total system load parameter is identical, constant and obviously do not meet The essence of load.In three of the above method, Component Based is owing to its physical model is clear, definite conception, it is simple to qualitative understanding is born Lotus characteristic, is widely used.
But there is certain shortcoming in traditional Component Based, including:
(1) to obtain load capacity not consistent with actual load power for inquiry agency, because there is the problem of simultaneity factor, and not all sets Standby is all to come into operation for 24 hours, accordingly, it would be desirable to carry out investigation statistics at times;
(2) As time goes on, actual load power, load structure and network structure all it may happen that change, if right Load is just thought once and for all after carrying out an investigation statistics modeling work, it is difficult to reach accuracy requirement;
(3) investigation work need to add up load composition and the parameter of thousands of users, and workload is huge, and is difficult to obtain standard True statistics.
The random time-dependent of synthetic load constituent is the essential reason that its part throttle characteristics has stochastic time-dependent.This time variation must The daily load curve so causing power consumer and transformer station has time variation.Therefore, the daily load curve of user and transformer station is inevitable Abundant information containing reflection load structure characteristic.
According to the convention of industry load classification, load can be divided into industrial load, Commercial Load and city dweller's load, rural area Load and this 4 class of other loads.Wherein rural area load type is one of more complicated load type, is also important bearing Lotus type, occupies larger specific gravity in each load type website.Rural area load curve has a characteristic that
Basic characteristics containing the rural area load irrigating load are as follows:
(1) rural area load changed also with season and the change of time, contained substantial amounts of irrigation load on spring and daytime in autumn, Contain certain heating load winter, and contain certain cooling load summer, but heating load and cooling load ratio are the most not Greatly.
(2) contain bigger irrigation load in spring and autumn, be the main feature of rural area load.Irrigate load and be generally present in white It the morning and afternoon, 3,4, May and 9,10, November irrigate the large percentage that load occupies;
(3) rural area load is smaller, and the regularity of load curve is not strong, and the proportion that special load curve occupies is higher.
(4) ratio that different rural area charge circuit irrigation loads occupies has certain difference.
(5) for the load curve of every day, peak load is generally present in evening, and the irrigation load on part circuit daytime is more, It is likely to occur by day.For load in evening, the load in general spring and autumn is relatively low, and the load in winter and summer is higher. For the load on daytime, the duty ratio in February is relatively low, other in month load the highest.
By to user and the signature analysis of transformer station's daily load curve, just can determine that each typical case's electricity consumption industry corresponding is used Electricity Constitution ratio of plant example and the electricity consumption industry composition of transformer station.Due to data acquisition and monitoring (SCADA) system and load Control management system (being called for short negative Ore-controlling Role) can provide the daily load curve of real-time transformer station's synthetic load and electric power to use respectively The daily load curve at family, therefore, the synthetic load composition thereby determined that will have online, the character of Real-time modeling set, thus Can fundamentally overcome the inherent shortcoming of conventional statistics synthesis load modeling.
Summary of the invention
In order to overcome above-mentioned the deficiencies in the prior art, the present invention provides a kind of rural area based on load curve decomposition load type load Modeling method, the shortcoming the method overcoming conventional statistics synthesis, can be resident load website quickly, easily and accurately Modeling, the load model generated by the present invention can improve the degree of accuracy that grid simulation calculates, and ensures power grid security, reliable, economic Ground runs.
In order to realize foregoing invention purpose, the present invention adopts the following technical scheme that:
The present invention provides a kind of rural area based on load curve decomposition load type load modeling method, and described method includes following step Rapid:
Step 1: decompose rural area load curve, and calculate the ratio of total amount shared by every kind of load in the load of rural area;
Step 2: calculate static load equivalent parameters, dynamic load equivalent parameters and distribution network system impedance, and it is negative to combine rural area The power supply area network topology data of lotus node, obtain rural area load type load model.
Described step 1 comprises the following steps:
Step 1-1: decompose rural area load curve;
Step 1-1-1: to the transformer station that supplied load type is rural area load, merges this transformer station 10kV or 6kV load outlet Daily load curve;
Step 1-1-2: calculate this transformer station's load Payload number of days monthly, if certain day load has null value, it is believed that this day nothing Effect, if certain day load has opposite sign phenomenon, be also considered as this day invalid;
Step 1-1-3: calculate this transformer station's load average daily load curve monthly and average monthly load curve;
Step 1-1-4: calculating every daily load curve period by day monthly, 6~the maximum of 17, and from February Average daily load curve selects the maximum of daytime period, and determines the load curve scope meeting reference load curve requirement, In the range from the maximum 0.8 of the daytime period calculated again to 1.15 times;
Step 1-1-5: start the load curve in all months is decomposed from February, if certain the daily load song of certain month The maximum of line period by day meets the daily load curve average value ranges that reference load curve requires, then it is assumed that this daily load is bent Line is the curve meeting reference load curve requirement, selects to decompose the curve meeting reference load curve requirement in month, calculates song The averaged curve of line, as reference load curve;
Step 1-1-6: determine the load curve mean value lower limit of satisfied big load curve requirement, this lower limit be daytime period 1.1 times of big value, if the maximum of certain daily load curve period by day is more than this lower limit, then it is assumed that this load song Line is the curve meeting big load curve requirement;
Step 1-1-7: determine the big load curve of this month, and calculate the mean value of datum curve load;
Step 1-2: calculate the ratio of total amount shared by every kind of load in the load of rural area;
Step 1-2-1: calculate each transformer station monthly average load value monthly and specify month to specify the average load value of day;
Step 1-2-2: calculate and specify month, specify day, the load value in appointment moment;
Step 1-2-3: calculate and include irrigating load, cooling load and the separate section load proportion k of heating loadp0And base Quasi-load proportion kf
In described step 1-1-3, calculate this transformer station's load average daily load curve monthly, all effective day monthly is born Lotus curve is added, the more effective number of days divided by this month;
Calculate this transformer station's load average monthly load curve monthly, by the load phase in 24 moment in monthly daily load curve every day Add again divided by the curve obtained by 24;In average monthly load curve, each effective day corresponding value is also referred to as daily load curve mean value, Daily load curve mean value is that the load in 24 moment of every day is added, then divided by the mean value obtained by 24;
In described step 1-1-5, if the number meeting the load curve of datum curve this month is less than 3, and above do not have There is the datum curve calculating any one month, then first do not calculate the curve of this month, carry out the calculating of next month curve;As Really this month meets the curve number of datum curve requirement more than or equal to 3, calculates datum curve, and all by this month meet benchmark The daily load curve of curve requirement is added, then meets the curve number of datum curve requirement, calculated datum curve divided by this month Smaller load curve as this month;Otherwise, the datum curve of month is used;
In described step 1-1-7, determine the big load curve of this month, if the load meeting big load curve requirement this month is bent Number of lines is less than 3, then ignore these load curves, do not calculate big load curve this month, does not carry out load curve and separates meter Calculate;If the load curve bar number meeting big load curve requirement this month is more than or equal to 3, calculate big load curve;
Calculate the mean value of datum curve load, big load curve, reference load curve are carried out standardization, average by three Value becomes 1.
In described step 1-2-1, calculate each transformer station monthly average load value monthly, by each in average monthly load curve monthly The value that effectively it is corresponding is added, the more effective number of days divided by this month;
Calculate the average load value specifying month to specify day, the monthly average load value in this month calculated is multiplied by the day specifying day The perunit value of load curve mean value;
In described step 1-2-2, calculate and specify month, specify day, the load value in appointment moment;
1) if only Smaller load curve, then Smaller load curve is used to calculate;
2) if there is big load curve and Smaller load curve simultaneously, the daily load curve mean value size of day is specified according to this month Determine that the big load curve of employing or Smaller load curve calculate;
Setting value is 1.1 times of the reference load mean value in certain month, then have:
If 2-1) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate;
If 2-2) daily load curve mean value is more than this definite value, using big load curve to calculate, big carry calculation result is The reference load mean value in this month calculated is multiplied by specifies little/big load curve middle finger timing of day to carve corresponding value in this month;
Also need to calculate and specify month, appointment day, the segregational load value in appointment moment, calculate big load curve and reference load curve Difference, as separate after load curve, this value equal to calculate specify month, specify day, specify the moment load value Deduct the reference load mean value in this month to be multiplied by and specify the Smaller load curve middle finger timing of day to carve corresponding value in this month, then divided by Calculate specifies month, appointment day, the load value in appointment moment;
In described step 1-2-3, calculate separate section load proportion kp0, by segregational load value divided by the appointment moon calculated Part, appointment day, the load value in appointment moment;Then reference load proportion kfFor 1-kp0
In described step 2, it is as follows that static load equivalent parameters calculates process:
By the multinomial static load model that the relationship description between load power and voltage is polynomial equation form, have
P=P0[a×(V/V0)2+b×(V/V0)+c] (1)
Q=Q0[α×(V/V0)2+β×(V/V0)+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the real-time voltage of static load, V0For The rated voltage of static load, P0And Q0It is respectively and represents at V0The specified active power of lower static load and reactive power, a, B and c is the active power coefficient of multinomial static load model, and α, β and γ are the idle of multinomial static load model Power coefficient;
The equivalence of static load is mainly a, b, c, P0With α, β, γ, Q0Equivalence, to multinomial static load The equivalence of model is based on the sensitivity to static load terminal voltage of the static load power, has
∂ P ∂ V | V = V 0 = Σ i = 1 n ∂ P i ∂ V | V = V 0 - - - ( 3 )
∂ Q ∂ V | V = V 0 = Σ i = 1 n ∂ Q i ∂ V | V = V 0 - - - ( 4 )
Wherein, n is static load number, PiAnd QiIt is respectively the meritorious of i-th static load and reactive power,With It is respectively i-th static load active power relative to the partial differential of voltage and i-th static load reactive power relative to voltage Partial differential;
Work as V=V0Time, have:
P 0 = Σ i = 1 n P 0 i - - - ( 5 )
Q 0 = Σ i = 1 n Q 0 i - - - ( 6 )
Wherein, P0iAnd Q0iIt is respectively initial active power and the reactive power of i-th static load;
a = Σ i = 1 n P 0 i a i P 0 - - - ( 7 )
b = Σ i = 1 n P 0 i b i P 0 - - - ( 8 )
c = Σ i = 1 n P 0 i c i P 0 - - - ( 9 )
Wherein, ai、biAnd ciIt is the multinomial static load model active power coefficient of i-th static load;
α = Σ i = 1 n Q 0 i α i Q 0 - - - ( 10 )
β = Σ i = 1 n Q 0 i β i Q 0 - - - ( 11 )
γ = Σ i = 1 n Q 0 i γ i Q 0 - - - ( 12 )
Wherein, αi、βiAnd γiIt is the multinomial static load model reactive power coefficient of i-th static load.
Described step 2 specifically includes following steps:
Step A: calculate all motor always absorbs active-power PΣWith reactive power QΣ, total electromagnetic power PΣem, total rotor Winding copper loss PΣcu2With total maximum electromagnetic power PΣem_max, specifically have:
P Σ = Σ j = 1 m P j
Q Σ = Σ j = 1 m Q i
P Σem = Σ j = 1 m P emj - - - ( 13 )
P Σcu 2 = Σ j = 1 m P cu 2 j
P Σem _ max = Σ j = 1 m P em _ max j
Wherein Pj、Qj、Pemj、Pcu2jAnd Pem_maxjRespectively represent the active power of jth platform equivalence motor, reactive power, Electromagnetic power, rotor windings copper loss and maximum electromagnetic power;M is equivalent motor number of units;
Step B: calculate the stator winding copper loss P of equivalent motorΣcu1And the slippage S of equivalence motor, it is expressed as:
PΣcu1=PΣ-PΣem(14)
S=PΣcu2/PΣem(15)
Initialize the maximum electromagnetic power P of equivalent motoremt_max, make Pemt_max=PΣem_max
Step C: calculate the stator winding resistance R of equivalent motors
First calculate the stator winding phase current of equivalent motorHave
I · Σ = - ( P Σ + jQ Σ 3 U · 1 ) * - - - ( 16 )
Wherein,Set end voltage for motor;
The then stator winding resistance R of equivalent motorsIt is expressed as:
R s = P Σcu 1 3 I Σ 2 - - - ( 17 )
Step D: calculate the equivalent impedance Z of equivalent machine modeldeq:
Z deq = U 1 2 P Σ - j Q Σ - - - ( 18 )
And have
Rdeq=real(Zdeq) (19)
Xdeq=imag(Zdeq) (20)
Wherein, RdeqAnd XdeqFor corresponding substitutional resistance and equivalent reactance;
Step E: calculate stator winding leakage reactance X of equivalent motorsWith rotor windings leakage reactance Xr, have
X s = X r = X 2 = - 1 2 ( U 1 2 2 P emt _ max - R s ) 2 - R s 2 - - - ( 21 )
Step F: through iteration stator winding leakage reactance X to equivalent motorsWith rotor windings leakage reactance XrIt is modified;
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor;
According to calculated Rs、Xs、XrAnd ZdeqIf, Kr=Rdeq-Rs, Kx=Xdeq-Xs, the most equivalent electronic The rotor windings resistance R of machinerReactance X excitatory with equivalencemIt is expressed as:
R r = ( K r + K x 2 / K r - ( K r + K x 2 / K r ) 2 - 4 X s 2 ) S 2 X m = K r X s + K x R r S R r S - K r - - - ( 23 )
Step H: use iterative method to calculate the maximum electromagnetic power of equivalent motor, and it is modified;
1) the maximum electromagnetic power P of the equivalent motor in calculating kth time iterationemt_maxk, it is expressed as
P emt _ max k = U 1 2 2 ( R s + R s 2 + X 2 ) - - - ( 24 )
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motordpIt is expressed as
Z dp = jX r + jX m ( R s + jX s ) R s + j ( X s + X m ) - - - ( 25 )
And have, Rdp=real(Zdp) and Xdp=imag(Zdp), RdpAnd XdpIt is respectively Thevenin's equivalence resistance and equivalence electricity Anti-;
Equivalent motor produces the condition of maximum electromagnetic power:
R pm = R r S m = R dp 2 + X dp 2 - - - ( 26 )
Wherein, RpmThe Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondencemFor critical slippage;
Thevenin's equivalence open-circuit voltageFor
V · 0 = U · 1 3 × jX m R s + j ( X s + X m ) - - - ( 27 )
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
P em _ max k = 3 U 0 2 R pm ( R dp + R pm ) 2 + X dp 2 - - - ( 28 )
2) the maximum electromagnetic power of equivalent motor is revised;
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic powermaxk, it is expressed as:
τ max k = P emt _ max k P em _ max k - - - ( 29 )
The most revised equivalent motor maximum electromagnetic power Pemt_maxIt is expressed as:
Pemt_maxmaxk×PΣem_max(30)
With Pemt_maxWith Pemt_maxkThe iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
Err P em _ max = | P emt _ max - P emt _ max k | - - - ( 31 )
WithFor iteration convergence standard, ifThen recalculate Xs、Rr、XrAnd Xm
Step I: calculate equivalent inertia time constant;
Equivalent inertia time constant H is expressed as:
H = Σ j = 1 m P nj H j Σ j = 1 m P nj - - - ( 32 )
Wherein PnjAnd HjRated power and inertia time constant for jth platform equivalence motor.
In described step 2, distribution network system impedance meter is shown as:
Z eq = [ Σ f = 1 M u f 2 ( 1 / Z f ) * ] / ( Σ l = 1 N I l ) 2 - - - ( 33 )
Wherein, ZeqFor distribution network system impedance;ufRepresent busbar voltage, ZfIndication transformer and distribution line impedance;IlTable Show load current, wherein, node total number that M and N is respectively in distribution network system and branch road total number.
Power supply area network topology data in described step 2 include distribution line data, transformer data and reactive-load compensation data.
Compared with prior art, the beneficial effects of the present invention is:
1. this method electrically-based system sales department Load Control Management System gathers the daily load curve of power consumer and number According to the daily load curve of the real-time transformer station synthetic load gathered with monitoring (SCADA) system offer, pass through load curve decomposition Obtain any time rural area Load Substation each electrical equipment type of confession and ratio, and the load utilizing traffic department to provide saves The power supply area network topology data of point, use Component Based, it is achieved the online load modeling of rural area load website, to reach It is the purpose of rural area load Site Modeling exactly, improves the degree of accuracy that grid simulation calculates, ensure power grid security, reliable, warp Ji ground runs.
2. this method daily load curve based on actual electric network data, have carried out online load modeling, have changed conventional statistics comprehensive Method relies on the pattern of off-line manual research, can carry out load modeling quick and easy, exactly, be greatly saved human and material resources, And avoiding the inaccurate problem of the load modeling caused owing to investigation result is inaccurate, this method is that tradition modeling pattern is entered one Step promotes, and provides the effect of important guidance for load modeling work.
3. (day such as the power consumer of marketing Load Control Management System collection is negative by diversified online data resource for this method The daily load curve of the real-time transformer station synthetic load that lotus curve and data acquisition and monitoring (SCADA) system provide and scheduling net Network topology information) realize onlineization of Component Based, load bus can be carried out real-time SLM modeling, model physical significance Clearly, strong adaptability, overcome the problem that traditional all load modeling methods are difficult in adapt to load time variation.
Accompanying drawing explanation
Fig. 1 is the load chart that Luohe becomes 35kV rural area load January;
Fig. 2 is the load chart that Luohe becomes 35kV rural area load May;
Fig. 3 is the load chart that Luohe becomes 35kV rural area load July;
Fig. 4 is the load chart that Luohe becomes 35kV rural area load November;
Fig. 5 is the load chart after Luohe becomes the decomposition in January of 35kV rural area load;
Fig. 6 is the load chart after Luohe becomes the decomposition in May of 35kV rural area load;
Fig. 7 is the load chart after Luohe becomes the decomposition in July of 35kV rural area load;
Fig. 8 is to become the load chart after the decomposition in November of 35kV rural area load for Luohe;
Fig. 9 is rural area based on load curve decomposition load type load modeling method flow diagram;
Figure 10 is the comparison diagram of 220kV load busbar voltage curve;
Figure 11 is the comparison diagram of 220kV load bus load active power curves;
Figure 12 is the comparison diagram of 220kV load bus load reactive capability curve.
Detailed description of the invention
Below in conjunction with the accompanying drawings the present invention is described in further detail.
The present invention provides a kind of rural area based on load curve decomposition load type load modeling method, and described method includes following step Rapid:
Step 1: decompose rural area load curve, and calculate the ratio of total amount shared by every kind of load in the load of rural area;
Step 2: calculate static load equivalent parameters, dynamic load equivalent parameters and distribution network system impedance, and it is negative to combine rural area The power supply area network topology data of lotus node, obtain rural area load type load model.
For the load curve basic characteristics of rural area load, the basic handling method of employing is as follows:
(1) smaller due to rural area load, and regular poor, therefore for the rural area load of same transformer station, Merge.
If do not merged, owing to the regularity of rural area load is poor, hardly result in preferable decomposition result.In Fig. 1-Fig. 4 Load curve be exactly that Luohe becomes the load curve after all 35kV rural areas charge circuit merges, load total after merging is relatively big, The Changing Pattern of load is the strongest.
(2) the rural area load curve decomposition containing irrigation load is base load, cooling load, heating load and irrigation load.
Cooling load mainly appears on the 6 of summer, 7, August, heating load mainly appears on the 12 of winter, 1, February, Irrigate load mainly appear on the 9 of the 3 of spring, 4, May and autumn, 10, November.Therefore, when decomposing, isolate Come load, 6,7, August be considered cooling load, 12,1, February be considered heating load, 3,4,5,9, 10, November be considered irrigate load.But the load separated in the part month of season alternation may be simultaneously present two kinds Load, owing to cannot be further continued for decomposing, is therefore approximately considered only exist that a kind of load, such as August separate negative Lotus may be simultaneously present cooling load and irrigates load, and being approximately considered the load separated August is all cooling load.
(3) basic skills of curve is separated: according to daytime period maximum, curve is decomposed, by the load curve of every month Being divided into two load curves, load curve on the basis of, another is big load curve, and the difference of two curves is as separation Load curve.
For the load curve of every month, it is classified as two suite lines according to the size of load, calculates the average bent of this two suites line Line, as typical load curve, load curve during different load level in approximate simulation every month.One therein Curve is base load curve, and this curve represents the change procedure of base load, should be less according to average load in all months Curve calculate.Urban residents' load contains irrigation load, separates the basic skills of basic load curve when being to use daytime Section maximum is decomposed, and selects the curve that daytime period maximum is less as reference load curve.
Decomposing two periods of main consideration, one is daytime period (6-17 point), and another is the period in the evening, reference load curve Should meet daytime load less, evening, load was the least.But knowable to the basic characteristics of rural area load curve, the 1. usual spring The load in season and evening in autumn is less, and the load in summer and evening in winter is relatively big, but that both differ and not very big;The whitest It load changes greatly in Various Seasonal, particularly season in spring and autumn;3. the load curve that load on daytime is less, evening load not Smaller.If selecting load on daytime and the smaller curve of load in evening as reference load curve, then meet benchmark The curve number of curve is less, can compare accurate decomposition base load curve only with daytime period maximum, therefore only with in vain It load maximum is decomposed.
Determine reference load curve daytime period maximum meet scope use two kinds of methods:
Calculate the maximum of all daily load curve daytime period, therefrom select a range of value, calculate the average of them Value;
Owing to daytime period load in February is minimum, and dispersiveness is less, and therefore calculating February, all daily load curves were white The maximum of it period, then calculates its mean value.
On the basis of above-mentioned mean value, the scope that selection reference load curve daytime period peak load meets, for every month Load, as long as daytime period peak load is within the scope of this, it is believed that be basic load curve;More than the maximum of this scope, It is considered big load curve.The average load curve calculating basic load curve is born greatly as Smaller load typical load curve, calculating The averaged curve of lotus curve is as the typical curve of big load curve, and both differences are as freezing, heating or irrigate load.
In two kinds of above-mentioned methods, first method is it needs to be determined that calculate the scope of datum curve mean value, and this scope is relatively difficult to Determine;Second method is relatively easy, and has certain adaptability.
Decompose rural area load curve process specific as follows:
Step 1-1-1: to the transformer station that supplied load type is rural area load, merges this transformer station 10kV or 6kV load outlet Daily load curve;
Step 1-1-2: calculate this transformer station's load Payload number of days monthly, if certain day load has null value, it is believed that this day nothing Effect, if certain day load has opposite sign phenomenon, be also considered as this day invalid;
Step 1-1-3: calculate this transformer station's load average daily load curve monthly and average monthly load curve;
Calculate this transformer station's load average daily load curve monthly, will be added by all effective daily load curve monthly, then remove Effective number of days with this month;
Calculate this transformer station's load average monthly load curve monthly, will the monthly load in 24 moment in daily load curve every day It is added again divided by the curve obtained by 24;In average monthly load curve, to be also referred to as daily load curve average for each effective day corresponding value Value, daily load curve mean value is the load in 24 moment of every day and is added, then divided by the mean value obtained by 24;
Step 1-1-4: calculating every daily load curve period by day monthly, i.e. 6~the maximum of 17, and from February Average daily load curve in select the maximum of daytime period, and determine the load curve scope meeting reference load curve requirement, In the range from the maximum 0.8 of the daytime period calculated again to 1.15 times;
Step 1-1-5: start the load curve in all months is decomposed from February, if certain the daily load song of certain month The maximum of line period by day meets the daily load curve average value ranges that reference load curve requires, then it is assumed that this daily load is bent Line is the curve meeting reference load curve requirement, selects to decompose the curve meeting reference load curve requirement in month, calculates song The averaged curve of line, as reference load curve;
If the number meeting the load curve of datum curve this month is less than 3, and does not the most calculate any one moon The datum curve of part, then first do not calculate the curve of this month, carry out the calculating of next month curve;If this moon meets datum curve The curve number required, more than or equal to 3, calculates datum curve, all of this month will meet the daily load that datum curve requires Curve is added, then meets the curve number of datum curve requirement divided by this month, and calculated datum curve is as the Smaller load of this month Curve;Otherwise, the datum curve of month is used;
Step 1-1-6: determine the load curve mean value lower limit of satisfied big load curve requirement, this lower limit be daytime period 1.1 times of big value, if the maximum of certain daily load curve period by day is more than this lower limit, then it is assumed that this load song Line is the curve meeting big load curve requirement;
Step 1-1-7: determine the big load curve of this month, and calculate the mean value of datum curve load;
Determine the big load curve of this month, if the load curve bar number meeting big load curve requirement this month is less than 3, then Ignore these load curves, do not calculate big load curve this month, do not carry out load curve and separate calculating;If this month meets The load curve bar number that big load curve requires, more than or equal to 3, calculates big load curve;
Calculate the mean value of datum curve load, big load curve, reference load curve are carried out standardization, will three flat Average becomes 1.
Accompanying drawing 5-8 be Luohe become 35kV rural area load 1,5,7, November decompose after load curve.
As shown in accompanying drawing 1-4 and accompanying drawing 5-8, January, load was relatively big, and dispersiveness is less, does not meet reference load curve The load curve required, uses the basic load curve in December, defines a big load curve, and the load separated is Cooling load.Irrigate the ratio that load occupies in May and November more, irrigate load and mainly appear on the morning and afternoon.7 The load separated in month is cooling load.
Calculate the ratio process of total amount shared by every kind of load in the load of rural area as follows:
Step 1-2-1: calculate each transformer station monthly average load value monthly and specify month to specify the average load value of day;
Calculate each transformer station monthly average load value monthly, will each effective day corresponding value phase in average monthly load curve monthly Add, the more effective number of days divided by this month;
Calculating the average load value specifying month to specify day, the monthly average load value in this month that will calculate is multiplied by specifies day The perunit value of daily load curve mean value;
Step 1-2-2: calculate and specify month, specify day, the load value in appointment moment;
1) if only Smaller load curve, then Smaller load curve is used to calculate;
2) if there is big load curve and Smaller load curve simultaneously, the daily load curve mean value size of day is specified according to this month Determine that the big load curve of employing or Smaller load curve calculate;
Setting value is 1.1 times of the reference load mean value in certain month, then have:
If 2-1) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate;
If 2-2) daily load curve mean value is more than this definite value, using big load curve to calculate, big carry calculation result is The reference load mean value in this month calculated is multiplied by specifies little/big load curve middle finger timing of day to carve corresponding value in this month;
Also need to calculate and specify month, appointment day, the segregational load value in appointment moment, calculate big load curve and reference load curve Difference, as separate after load curve, this value equal to calculate specify month, specify day, specify the moment load value Deduct the reference load mean value in this month to be multiplied by and specify the Smaller load curve middle finger timing of day to carve corresponding value in this month, then divided by Calculate specifies month, appointment day, the load value in appointment moment;
Step 1-2-3: calculate and include irrigating load, cooling load and the separate section load proportion k of heating loadp0And base Quasi-load proportion kf
Calculate separate section load proportion kp0, will segregational load value divided by the appointment month calculated, specify day, appointment The load value in moment;Then reference load proportion kfFor 1-kp0
In described step 2, it is as follows that static load equivalent parameters calculates process:
By the multinomial static load model that the relationship description between load power and voltage is polynomial equation form, have
P=P0[a×(V/V0)2+b×(V/V0)+c] (1)
Q=Q0[α×(V/V0)2+β×(V/V0)+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the real-time voltage of static load, V0For The rated voltage of static load, P0And Q0It is respectively and represents at V0The specified active power of lower static load and reactive power, a, B and c is the active power coefficient of multinomial static load model, and α, β and γ are the idle of multinomial static load model Power coefficient;
The equivalence of static load is mainly a, b, c, P0With α, β, γ, Q0Equivalence, to multinomial static load The equivalence of model is based on the sensitivity to static load terminal voltage of the static load power, i.e.
∂ P ∂ V | V = V 0 = Σ i = 1 n ∂ P i ∂ V | V = V 0 - - - ( 3 )
∂ Q ∂ V | V = V 0 = Σ i = 1 n ∂ Q i ∂ V | V = V 0 - - - ( 4 )
Wherein, n is static load number, PiAnd QiIt is respectively the meritorious of i-th static load and reactive power,With It is respectively i-th static load active power relative to the partial differential of voltage and i-th static load reactive power relative to voltage Partial differential;
Work as V=V0Time, have:
P 0 = Σ i = 1 n P 0 i - - - ( 5 )
Q 0 = Σ i = 1 n Q 0 i - - - ( 6 )
Wherein, P0iAnd Q0iIt is respectively initial active power and the reactive power of i-th static load;
a = Σ i = 1 n P 0 i a i P 0 - - - ( 7 )
b = Σ i = 1 n P 0 i b i P 0 - - - ( 8 )
c = Σ i = 1 n P 0 i c i P 0 - - - ( 9 )
Wherein, ai、biAnd ciIt is the multinomial static load model active power coefficient of i-th static load;
α = Σ i = 1 n Q 0 i α i Q 0 - - - ( 10 )
β = Σ i = 1 n Q 0 i β i Q 0 - - - ( 11 )
γ = Σ i = 1 n Q 0 i γ i Q 0 - - - ( 12 )
Wherein, αi、βiAnd γiIt is the multinomial static load model reactive power coefficient of i-th static load.
Described step 2 specifically includes following steps:
Step A: calculate all motor always absorbs active-power PΣWith reactive power QΣ, total electromagnetic power PΣem, total rotor Winding copper loss PΣcu2With total maximum electromagnetic power PΣem_max, specifically have:
P Σ = Σ j = 1 m P j
Q Σ = Σ j = 1 m Q i
P Σem = Σ j = 1 m P emj - - - ( 13 )
P Σcu 2 = Σ j = 1 m P cu 2 j
P Σem _ max = Σ j = 1 m P em _ max j
Wherein Pj、Qj、Pemj、Pcu2jAnd Pem_maxjRespectively represent the active power of jth platform equivalence motor, reactive power, Electromagnetic power, rotor windings copper loss and maximum electromagnetic power;M is equivalent motor number of units;
Step B: calculate the stator winding copper loss P of equivalent motorΣcu1And the slippage S of equivalence motor, it is expressed as:
PΣcu1=PΣ-PΣem(14)
S=PΣcu2/PΣem(15)
Initialize the maximum electromagnetic power P of equivalent motoremt_max, make Pemt_max=PΣem_max
Step C: calculate the stator winding resistance R of equivalent motors
First calculate the stator winding phase current of equivalent motorHave
I · Σ = - ( P Σ + jQ Σ 3 U · 1 ) * - - - ( 16 )
Wherein,Set end voltage for motor;
The then stator winding resistance R of equivalent motorsIt is expressed as:
R s = P Σcu 1 3 I Σ 2 - - - ( 17 )
Step D: calculate the equivalent impedance Z of equivalent machine modeldeq:
Z deq = U 1 2 P Σ - jQ Σ - - - ( 18 )
And have
Rdeq=real(Zdeq) (19)
Xdeq=imag(Zdeq) (20)
Wherein, RdeqAnd XdeqFor corresponding substitutional resistance and equivalent reactance;
Step E: calculate stator winding leakage reactance X of equivalent motorsWith rotor windings leakage reactance Xr, have
X s = X r = X 2 = - 1 2 ( U 1 2 2 P emt _ max - R s ) 2 - R s 2 - - - ( 21 )
Step F: through iteration stator winding leakage reactance X to equivalent motorsWith rotor windings leakage reactance XrIt is modified;
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor;
According to calculated Rs、Xs、XrAnd ZdeqIf, Kr=Rdeq-Rs, Kx=Xdeq-Xs, the most equivalent electronic The rotor windings resistance R of machinerReactance X excitatory with equivalencemIt is expressed as:
R r = ( K r + K x 2 / K r - ( K r + K x 2 / K r ) 2 - 4 X s 2 ) S 2 X m = K r X s + K x R r S R r S - K r - - - ( 23 )
Step H: use iterative method to calculate the maximum electromagnetic power of equivalent motor, and it is modified;
1) the maximum electromagnetic power P of the equivalent motor in calculating kth time iterationemt_maxk, it is expressed as
P emt _ max k = U 1 2 2 ( R s + R s 2 + X 2 ) - - - ( 24 )
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motordpIt is expressed as
Z dp = jX r + jX m ( R s + jX s ) R s + j ( X s + X m ) - - - ( 25 )
And have, Rdp=real(Zdp) and Xdp=imag(Zdp), RdpAnd XdpIt is respectively Thevenin's equivalence resistance and equivalence electricity Anti-;
Equivalent motor produces the condition of maximum electromagnetic power:
R pm = R r S m = R dp 2 + X dp 2 - - - ( 26 )
Wherein, RpmThe Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondencemFor critical slippage;
Thevenin's equivalence open-circuit voltageFor
V · 0 = U · 1 3 × jX m R s + j ( X s + X m ) - - - ( 27 )
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
P em _ max k = 3 U 0 2 R pm ( R dp + R pm ) 2 + X dp 2 - - - ( 28 )
2) the maximum electromagnetic power of equivalent motor is revised;
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic powermaxk, it is expressed as:
τ max k = P emt _ max k P em _ max k - - - ( 29 )
The most revised equivalent motor maximum electromagnetic power Pemt_maxIt is expressed as:
Pemt_maxmaxk×PΣem_max(30)
With Pemt_maxWith Pemt_maxkThe iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
Err P em _ max = | P emt _ max - P emt _ max k | - - - ( 31 )
WithFor iteration convergence standard, ifThen recalculate Xs、Rr、XrAnd Xm
Step I: calculate equivalent inertia time constant;
Equivalent inertia time constant H is expressed as:
H = Σ j = 1 m P nj H j Σ j = 1 m P nj - - - ( 32 )
Wherein PnjAnd HjRated power and inertia time constant for jth platform equivalence motor.
In described step 2, distribution network system impedance meter is shown as:
Z eq = [ Σ f = 1 M u f 2 ( 1 / Z f ) * ] / ( Σ l = 1 N I l ) 2 - - - ( 33 )
Wherein, ZeqFor distribution network system impedance;ufRepresent busbar voltage, ZfIndication transformer and distribution line impedance;IlTable Show load current, wherein, node total number that M and N is respectively in distribution network system and branch road total number.
Power supply area network topology data in described step 2 include distribution line data, transformer data and reactive-load compensation data.
Accompanying drawing 10-12 is the simulation curve comparison diagram with measured curve of the SLM Equivalent Model parameter using the present invention to draw
Accompanying drawing 10 is the comparison diagram of 220kV load busbar voltage curve;Wherein solid line is measured curve, and dotted line is for using the present invention Simulation result during the equivalent SLM model parameter drawn;
Accompanying drawing 11 is the comparison diagram of 220kV load bus load active power curves;Wherein solid line is measured curve, and dotted line is for adopting Simulation result during the equivalent SLM model parameter drawn by the present invention;
Accompanying drawing 12 is the comparison diagram of 220kV load bus load reactive capability curve;Wherein solid line is measured curve, and dotted line is for adopting Simulation result during the equivalent SLM model parameter drawn by the present invention.
Technical scheme is applied in the Henan Electric Power System stable model of Henan Electric Power Company is analyzed and researched.
Finally should be noted that: above example is only in order to illustrate that technical scheme is not intended to limit, although reference The present invention has been described in detail by above-described embodiment, those of ordinary skill in the field it is understood that still can to this Invention detailed description of the invention modify or equivalent, and without departing from spirit and scope of the invention any amendment or etc. With replacing, it all should be contained in the middle of scope of the presently claimed invention.

Claims (7)

1. rural area based on a load curve decomposition load type load modeling method, it is characterised in that: described method include with Lower step:
Step 1: decompose rural area load curve, and calculate the ratio of total amount shared by every kind of load in the load of rural area;
Step 2: calculate static load equivalent parameters, dynamic load equivalent parameters and distribution network system impedance, and it is negative to combine rural area The power supply area network topology data of lotus node, obtain rural area load type load model;
Described step 1 comprises the following steps:
Step 1-1: decompose rural area load curve;
Step 1-1-1: to the transformer station that supplied load type is rural area load, merges this transformer station 10kV or 6kV load outlet Daily load curve;
Step 1-1-2: calculate this transformer station's load Payload number of days monthly, if certain day load has null value, it is believed that this day nothing Effect, if certain day load has opposite sign phenomenon, be also considered as this day invalid;
Step 1-1-3: calculate this transformer station's load average daily load curve monthly and average monthly load curve;
Step 1-1-4: calculating every daily load curve period by day monthly, 6~the maximum of 17, and from February Average daily load curve selects the maximum of daytime period, and determines the load curve scope meeting reference load curve requirement, In the range from the maximum 0.8 of the daytime period calculated again to 1.15 times;
Step 1-1-5: start the load curve in all months is decomposed from February, if certain the daily load song of certain month The maximum of line period by day meets the daily load curve average value ranges that reference load curve requires, then it is assumed that this daily load is bent Line is the curve meeting reference load curve requirement, selects to decompose the curve meeting reference load curve requirement in month, calculates song The averaged curve of line, as reference load curve;
Step 1-1-6: determine the load curve mean value lower limit of satisfied big load curve requirement, this lower limit be daytime period 1.1 times of big value, if the maximum of certain daily load curve period by day is more than this lower limit, then it is assumed that this load song Line is the curve meeting big load curve requirement;
Step 1-1-7: determine the big load curve of this month, and calculate the mean value of datum curve load;
Step 1-2: calculate the ratio of total amount shared by every kind of load in the load of rural area;
Step 1-2-1: calculate each transformer station monthly average load value monthly and specify month to specify the average load value of day;
Step 1-2-2: calculate and specify month, specify day, the load value in appointment moment;
Step 1-2-3: calculate and include irrigating load, cooling load and the separate section load proportion k of heating loadp0And base Quasi-load proportion kf
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: in described step 1-1-3, calculate this transformer station's load average daily load curve monthly, all effective day monthly is born Lotus curve is added, the more effective number of days divided by this month;
Calculate this transformer station's load average monthly load curve monthly, by the load phase in 24 moment in monthly daily load curve every day Add again divided by the curve obtained by 24;In average monthly load curve, each effective day corresponding value is also referred to as daily load curve mean value, Daily load curve mean value is that the load in 24 moment of every day is added, then divided by the mean value obtained by 24;
In described step 1-1-5, if the number meeting the load curve of datum curve this month is less than 3, and above do not have There is the datum curve calculating any one month, then first do not calculate the curve of this month, carry out the calculating of next month curve;As Really this month meets the curve number of datum curve requirement more than or equal to 3, calculates datum curve, and all by this month meet benchmark The daily load curve of curve requirement is added, then meets the curve number of datum curve requirement, calculated datum curve divided by this month Smaller load curve as this month;Otherwise, the datum curve of month is used;
In described step 1-1-7, determine the big load curve of this month, if the load meeting big load curve requirement this month is bent Number of lines is less than 3, then ignore these load curves, do not calculate big load curve this month, does not carry out load curve and separates meter Calculate;If the load curve bar number meeting big load curve requirement this month is more than or equal to 3, calculate big load curve;
Calculate the mean value of datum curve load, big load curve, reference load curve are carried out standardization, average by three Value becomes 1.
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: in described step 1-2-1, calculate each transformer station monthly average load value monthly, by each in average monthly load curve monthly The value that effectively it is corresponding is added, the more effective number of days divided by this month;
Calculate the average load value specifying month to specify day, the monthly average load value in this month calculated is multiplied by the day specifying day The perunit value of load curve mean value;
In described step 1-2-2, calculate and specify month, specify day, the load value in appointment moment;
1) if only Smaller load curve, then Smaller load curve is used to calculate;
2) if there is big load curve and Smaller load curve simultaneously, the daily load curve mean value size of day is specified according to this month Determine that the big load curve of employing or Smaller load curve calculate;
Setting value is 1.1 times of the reference load mean value in certain month, then have:
If 2-1) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate;
If 2-2) daily load curve mean value is more than this definite value, using big load curve to calculate, big carry calculation result is The reference load mean value in this month calculated is multiplied by specifies little/big load curve middle finger timing of day to carve corresponding value in this month;
Also need to calculate and specify month, appointment day, the segregational load value in appointment moment, calculate big load curve and reference load curve Difference, as separate after load curve, this value equal to calculate specify month, specify day, specify the moment load value Deduct the reference load mean value in this month to be multiplied by and specify the Smaller load curve middle finger timing of day to carve corresponding value in this month, then divided by Calculate specifies month, appointment day, the load value in appointment moment;
In described step 1-2-3, calculate separate section load proportion kp0, by segregational load value divided by the appointment moon calculated Part, appointment day, the load value in appointment moment;Then reference load proportion kfFor 1-kp0
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: in described step 2, it is as follows that static load equivalent parameters calculates process:
By the multinomial static load model that the relationship description between load power and voltage is polynomial equation form, have
P=P0[a×(V/V0)2+b×(V/V0)+c] (1)
Q=Q0[α×(V/V0)2+β×(V/V0)+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the real-time voltage of static load, V0For The rated voltage of static load, P0And Q0It is respectively and represents at V0The specified active power of lower static load and reactive power, a, B and c is the active power coefficient of multinomial static load model, and α, β and γ are the idle of multinomial static load model Power coefficient;
The equivalence of static load is mainly a, b, c, P0With α, β, γ, Q0Equivalence, to multinomial static load The equivalence of model is based on the sensitivity to static load terminal voltage of the static load power, has
∂ P ∂ V | V = V 0 = Σ i = 1 n ∂ P i ∂ V | V = V 0 - - - ( 3 )
∂ Q ∂ V | V = V 0 = Σ i = 1 n ∂ Q i ∂ V | V = V 0 - - - ( 4 )
Wherein, n is static load number, PiAnd QiIt is respectively the meritorious of i-th static load and reactive power,With It is respectively i-th static load active power relative to the partial differential of voltage and i-th static load reactive power relative to voltage Partial differential;
Work as V=V0Time, have:
P 0 = Σ i = 1 n P 0 i - - - ( 5 )
Q 0 = Σ i = 1 n Q 0 i - - - ( 6 )
Wherein, P0iAnd Q0iIt is respectively initial active power and the reactive power of i-th static load;
a = Σ i = 1 n P 0 i a i P 0 - - - ( 7 )
b = Σ i = 1 n P 0 i b i P 0 - - - ( 8 )
c = Σ i = 1 n P 0 i c i P 0 - - - ( 9 )
Wherein, ai、biAnd ciIt is the multinomial static load model active power coefficient of i-th static load;
α = Σ i = 1 n Q 0 i α i Q 0 - - - ( 10 )
β = Σ i = 1 n Q 0 i β i Q 0 - - - ( 11 )
γ = Σ i = 1 n Q 0 i γ i Q 0 - - - ( 12 )
Wherein, αi、βiAnd γiIt is the multinomial static load model reactive power coefficient of i-th static load.
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: described step 2 specifically includes following steps:
Step A: calculate all motor always absorbs active-power PWith reactive power Q, total electromagnetic power P∑em, total rotor Winding copper loss P∑cu2With total maximum electromagnetic power P∑em_max, specifically have:
P Σ = Σ j = 1 m P j Q Σ = Σ j = 1 m Q i P Σ e m = Σ j = 1 m P e m j P Σ c u 2 = Σ j = 1 m P c u 2 j P Σ e m _ max = Σ j = 1 m P e m _ max j - - - ( 13 )
Wherein Pj、Qj、Pemj、Pcu2jAnd Pem_maxjRespectively represent the active power of jth platform equivalence motor, reactive power, Electromagnetic power, rotor windings copper loss and maximum electromagnetic power;M is equivalent motor number of units;
Step B: calculate the stator winding copper loss P of equivalent motor∑cu1And the slippage S of equivalence motor, it is expressed as:
P∑cu1=P-P∑em (14)
S=P∑cu2/P∑em (15)
Initialize the maximum electromagnetic power P of equivalent motoremt_max, make Pemt_max=P∑em_max
Step C: calculate the stator winding resistance R of equivalent motors
First calculate the stator winding phase current of equivalent motorHave
I · Σ = - ( P Σ + jQ Σ 3 U · 1 ) * - - - ( 16 )
Wherein,Set end voltage for motor;
The then stator winding resistance R of equivalent motorsIt is expressed as:
R s = P Σ c u 1 3 I Σ 2 - - - ( 17 )
Step D: calculate the equivalent impedance Z of equivalent machine modeldeq:
Z d e q = U 1 2 P Σ - jQ Σ - - - ( 18 )
And have
Rdeq=real (Zdeq) (19)
Xdeq=imag (Zdeq) (20)
Wherein, RdeqAnd XdeqFor corresponding substitutional resistance and equivalent reactance;
Step E: calculate stator winding leakage reactance X of equivalent motorsWith rotor windings leakage reactance Xr, have
X s = X r = X 2 = - 1 2 ( U 1 2 2 P e m t _ m a x - R s ) 2 - R s 2 - - - ( 21 )
Step F: through iteration stator winding leakage reactance X to equivalent motorsWith rotor windings leakage reactance XrIt is modified;
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor;
According to calculated Rs、Xs、XrAnd ZdeqIf, Kr=Rdeq-Rs, Kx=Xdeq-Xs, the most equivalent electronic The rotor windings resistance R of machinerReactance X excitatory with equivalencemIt is expressed as:
R r = ( K r + K x 2 / K r - ( K r + K x 2 / K r ) 2 - 4 X s 2 ) S 2 X m = K r X s + K x R r S R r S - K r - - - ( 23 )
Step H: use iterative method to calculate the maximum electromagnetic power of equivalent motor, and it is modified;
1) the maximum electromagnetic power P of the equivalent motor in calculating kth time iterationemt_maxk, it is expressed as
P e m t _ max k = U 1 2 2 ( R s + R s 2 + X 2 ) - - - ( 24 )
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motordpIt is expressed as
Z d p = jX r + jX m ( R s + jX s ) R s + j ( X s + X m ) - - - ( 25 )
And have, Rdp=real (Zdp) and Xdp=imag (Zdp), RdpAnd XdpIt is respectively Thevenin's equivalence resistance and equivalence electricity Anti-;
Equivalent motor produces the condition of maximum electromagnetic power:
R p m = R r S m = R d p 2 + X d p 2 - - - ( 26 )
Wherein, RpmThe Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondencemFor critical slippage;
Thevenin's equivalence open-circuit voltageFor
V · 0 = U · 1 3 × jX m R s + j ( X s + X m ) - - - ( 27 )
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
P e m _ max k = 3 V 0 2 R p m ( R d p + R p m ) 2 + X d p 2 - - - ( 28 )
2) the maximum electromagnetic power of equivalent motor is revised;
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic powermaxk, it is expressed as:
τ max k = P e m t _ max k P e m _ max k - - - ( 29 )
The most revised equivalent motor maximum electromagnetic power Pemt_maxIt is expressed as:
Pemt_maxmaxk×P∑em_max (30)
With Pemt_maxWith Pemt_maxkThe iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
Err P e m _ m a x = | P e m t _ m a x - P e m t _ max k | - - - ( 31 )
WithFor iteration convergence standard, ifThen recalculate Xs、Rr、XrAnd Xm
Step I: calculate equivalent inertia time constant;
Equivalent inertia time constant H is expressed as:
H = Σ j = 1 m P n j H j Σ j = 1 m P n j - - - ( 32 )
Wherein PnjAnd HjRated power and inertia time constant for jth platform equivalence motor.
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: in described step 2, distribution network system impedance meter is shown as:
Z e q = [ Σ f = 1 M u f 2 ( 1 / Z f ) * ] / ( Σ l = 1 N I l ) 2 - - - ( 33 )
Wherein, ZeqFor distribution network system impedance;ufRepresent busbar voltage, ZfIndication transformer and distribution line impedance;IlTable Show load current, wherein, node total number that M and N is respectively in distribution network system and branch road total number.
Rural area based on load curve decomposition the most according to claim 1 load type load modeling method, its feature exists In: the power supply area network topology data in described step 2 include distribution line data, transformer data and reactive-load compensation data.
CN201310500676.2A 2013-10-22 2013-10-22 A kind of rural area based on load curve decomposition load type load modeling method Active CN103632031B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310500676.2A CN103632031B (en) 2013-10-22 2013-10-22 A kind of rural area based on load curve decomposition load type load modeling method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310500676.2A CN103632031B (en) 2013-10-22 2013-10-22 A kind of rural area based on load curve decomposition load type load modeling method

Publications (2)

Publication Number Publication Date
CN103632031A CN103632031A (en) 2014-03-12
CN103632031B true CN103632031B (en) 2016-08-31

Family

ID=50213070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310500676.2A Active CN103632031B (en) 2013-10-22 2013-10-22 A kind of rural area based on load curve decomposition load type load modeling method

Country Status (1)

Country Link
CN (1) CN103632031B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104348154B (en) * 2014-11-22 2016-08-24 郏县供电有限责任公司 The dispatching method of a kind of power distribution network and device
CN105224812B (en) * 2015-10-21 2018-10-09 中国电力科学研究院 A kind of static load frequency factor polymerization in load model
CN107453346B (en) * 2016-05-30 2022-03-18 中国电力科学研究院 Load curve time interval division method suitable for power distribution network reconstruction
CN107069717B (en) * 2017-05-11 2020-05-15 国家电网公司 Power distribution network load modeling method
CN107392391A (en) * 2017-08-10 2017-11-24 天津求实智源科技有限公司 A kind of electric load electricity consumption decomposition method based on deep learning
CN114036762A (en) * 2021-11-15 2022-02-11 贵州电网有限责任公司 System and method for predicting load of power distribution network based on proportionality coefficient method

Citations (4)

* 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
CN102521674A (en) * 2011-12-07 2012-06-27 福建省电力有限公司 Computer generating method for short-term load forecasting in summer in power consumer electricity-utilization information acquisition system
CN103279803A (en) * 2013-04-27 2013-09-04 深圳供电局有限公司 Load modeling method and system based on comprehensive information theory and modern interior point theory

Patent Citations (4)

* 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
CN102521674A (en) * 2011-12-07 2012-06-27 福建省电力有限公司 Computer generating method for short-term load forecasting in summer in power consumer electricity-utilization information acquisition system
CN103279803A (en) * 2013-04-27 2013-09-04 深圳供电局有限公司 Load modeling method and system based on comprehensive information theory and modern interior point theory

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于日负荷曲线的负荷分类和综合建模;鞠平 等;《电力系统自动化》;20060815;第30卷(第16期);第6-9页 *
综合负荷模型参数的深化研究及适应性分析;汤涌 等;《电网技术》;20100228;第34卷(第2期);第57-63页 *
考虑配电网络的综合负荷模型;汤涌 等;《电网技术》;20070331;第31卷(第5期);第34-38页 *

Also Published As

Publication number Publication date
CN103632031A (en) 2014-03-12

Similar Documents

Publication Publication Date Title
CN103632031B (en) A kind of rural area based on load curve decomposition load type load modeling method
CN103514571B (en) A kind of business resident load type load modeling method based on load curve decomposition
CN101777765B (en) On-line load simulation method of power system
CN102116857B (en) Method for checking electric quantity at power grid gate
CN102184209A (en) Simulation data accessing method based on power grid CIM (Common Information Model) interface
Choi et al. Feasibility study: Autonomous state estimation in distribution systems
CN103887792B (en) A kind of low-voltage distribution network modeling method containing distributed power source
CN105375477A (en) Rural power distribution network line loss calculation method and system
CN108306284B (en) Online load modeling method based on local intelligent measurement
CN105427186A (en) Power distribution network line loss calculation method based on improved equivalent electric resistance method
CN106992517A (en) Electric power line loss analysis method, device and system
CN105528671A (en) Power grid multidimensional sensing and safety assessment system and method on the basis of big data
CN102651049A (en) Method for calculating loss reducing rate of newly built transformer station of electric power system
CN108258710A (en) A kind of battery energy storage system Optimal Configuration Method counted and battery capacity decays
CN102856896A (en) On-line analytical method for direct-current transmission loss
CN109494719B (en) Hierarchical impedance analysis method for medium-low voltage hybrid power distribution network
CN102738791B (en) On-line static security analysis method based on topological comparison
CN103679551A (en) Method for on-line calculation of power flow of power distribution network
CN103123715B (en) 220KV transformer substation electricity-using industry composition proportion online real-time counting method
Zhang et al. Dynamic synthesis load modeling approach based on load survey and load curves analysis
CN202190122U (en) System for analyzing theoretical line loss rate of distribution network low voltage district
CN111398859B (en) User low-voltage cause big data analysis method and system
CN107257130A (en) The low-voltage network loss computing method of decoupling is measured based on region
CN104009469B (en) A kind of low-voltage distribution network platform district network topology isopleth map collecting method
CN109286188B (en) 10kV power distribution network theoretical line loss calculation method based on multi-source data set

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model