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 PDFInfo
 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
Links
 238000000354 decomposition reaction Methods 0.000 title claims abstract description 24
 230000003068 static Effects 0.000 claims abstract description 67
 230000002354 daily Effects 0.000 claims description 61
 238000004804 winding Methods 0.000 claims description 17
 230000000875 corresponding Effects 0.000 claims description 14
 238000001816 cooling Methods 0.000 claims description 12
 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 10
 240000002853 Nelumbo nucifera Species 0.000 claims description 10
 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 10
 RYGMFSIKBFXOCRUHFFFAOYSAN copper Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwLjAnIGhlaWdodD0nMzAwLjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PScxMzguMCcgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5DPC90ZXh0Pgo8dGV4dCB4PScxNjUuNicgeT0nMTcwLjAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID51PC90ZXh0Pgo8cGF0aCBkPSdNIDE4OC40LDE1MC4wIEwgMTg4LjQsMTQ5LjggTCAxODguNCwxNDkuNyBMIDE4OC40LDE0OS41IEwgMTg4LjMsMTQ5LjMgTCAxODguMywxNDkuMiBMIDE4OC4yLDE0OS4wIEwgMTg4LjEsMTQ4LjkgTCAxODguMCwxNDguNyBMIDE4Ny45LDE0OC42IEwgMTg3LjcsMTQ4LjUgTCAxODcuNiwxNDguNCBMIDE4Ny41LDE0OC4zIEwgMTg3LjMsMTQ4LjIgTCAxODcuMiwxNDguMSBMIDE4Ny4wLDE0OC4xIEwgMTg2LjgsMTQ4LjAgTCAxODYuNywxNDguMCBMIDE4Ni41LDE0OC4wIEwgMTg2LjMsMTQ4LjAgTCAxODYuMSwxNDguMCBMIDE4Ni4wLDE0OC4xIEwgMTg1LjgsMTQ4LjEgTCAxODUuNiwxNDguMiBMIDE4NS41LDE0OC4yIEwgMTg1LjMsMTQ4LjMgTCAxODUuMiwxNDguNCBMIDE4NS4xLDE0OC41IEwgMTg0LjksMTQ4LjcgTCAxODQuOCwxNDguOCBMIDE4NC43LDE0OC45IEwgMTg0LjcsMTQ5LjEgTCAxODQuNiwxNDkuMiBMIDE4NC41LDE0OS40IEwgMTg0LjUsMTQ5LjYgTCAxODQuNSwxNDkuNyBMIDE4NC40LDE0OS45IEwgMTg0LjQsMTUwLjEgTCAxODQuNSwxNTAuMyBMIDE4NC41LDE1MC40IEwgMTg0LjUsMTUwLjYgTCAxODQuNiwxNTAuOCBMIDE4NC43LDE1MC45IEwgMTg0LjcsMTUxLjEgTCAxODQuOCwxNTEuMiBMIDE4NC45LDE1MS4zIEwgMTg1LjEsMTUxLjUgTCAxODUuMiwxNTEuNiBMIDE4NS4zLDE1MS43IEwgMTg1LjUsMTUxLjggTCAxODUuNiwxNTEuOCBMIDE4NS44LDE1MS45IEwgMTg2LjAsMTUxLjkgTCAxODYuMSwxNTIuMCBMIDE4Ni4zLDE1Mi4wIEwgMTg2LjUsMTUyLjAgTCAxODYuNywxNTIuMCBMIDE4Ni44LDE1Mi4wIEwgMTg3LjAsMTUxLjkgTCAxODcuMiwxNTEuOSBMIDE4Ny4zLDE1MS44IEwgMTg3LjUsMTUxLjcgTCAxODcuNiwxNTEuNiBMIDE4Ny43LDE1MS41IEwgMTg3LjksMTUxLjQgTCAxODguMCwxNTEuMyBMIDE4OC4xLDE1MS4xIEwgMTg4LjIsMTUxLjAgTCAxODguMywxNTAuOCBMIDE4OC4zLDE1MC43IEwgMTg4LjQsMTUwLjUgTCAxODguNCwxNTAuMyBMIDE4OC40LDE1MC4yIEwgMTg4LjQsMTUwLjAgTCAxODYuNCwxNTAuMCBaJyBzdHlsZT0nZmlsbDojMDAwMDAwO2ZpbGwtcnVsZTpldmVub2RkO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTojMDAwMDAwO3N0cm9rZS13aWR0aDowLjBweDtzdHJva2UtbGluZWNhcDpidXR0O3N0cm9rZS1saW5lam9pbjptaXRlcjtzdHJva2Utb3BhY2l0eToxOycgLz4KPC9zdmc+Cg== data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NS4wJyBoZWlnaHQ9Jzg1LjAnIHg9JzAuMCcgeT0nMC4wJz4gPC9yZWN0Pgo8dGV4dCB4PSczNS4wJyB5PSc1My42JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+QzwvdGV4dD4KPHRleHQgeD0nNTEuMCcgeT0nNTMuNicgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPnU8L3RleHQ+CjxwYXRoIGQ9J00gNjYuNCw0Mi4wIEwgNjYuNCw0MS45IEwgNjYuNCw0MS44IEwgNjYuNCw0MS43IEwgNjYuMyw0MS42IEwgNjYuMyw0MS41IEwgNjYuMiw0MS40IEwgNjYuMiw0MS4zIEwgNjYuMSw0MS4zIEwgNjYuMSw0MS4yIEwgNjYuMCw0MS4xIEwgNjUuOSw0MS4xIEwgNjUuOCw0MS4wIEwgNjUuNyw0MS4wIEwgNjUuNiw0MC45IEwgNjUuNSw0MC45IEwgNjUuNSw0MC45IEwgNjUuNCw0MC44IEwgNjUuMyw0MC44IEwgNjUuMiw0MC44IEwgNjUuMSw0MC45IEwgNjUuMCw0MC45IEwgNjQuOSw0MC45IEwgNjQuOCw0MC45IEwgNjQuNyw0MS4wIEwgNjQuNiw0MS4wIEwgNjQuNSw0MS4xIEwgNjQuNCw0MS4yIEwgNjQuNCw0MS4yIEwgNjQuMyw0MS4zIEwgNjQuMiw0MS40IEwgNjQuMiw0MS41IEwgNjQuMiw0MS42IEwgNjQuMSw0MS43IEwgNjQuMSw0MS44IEwgNjQuMSw0MS45IEwgNjQuMSw0Mi4wIEwgNjQuMSw0Mi4wIEwgNjQuMSw0Mi4xIEwgNjQuMSw0Mi4yIEwgNjQuMSw0Mi4zIEwgNjQuMiw0Mi40IEwgNjQuMiw0Mi41IEwgNjQuMiw0Mi42IEwgNjQuMyw0Mi43IEwgNjQuNCw0Mi44IEwgNjQuNCw0Mi44IEwgNjQuNSw0Mi45IEwgNjQuNiw0My4wIEwgNjQuNyw0My4wIEwgNjQuOCw0My4xIEwgNjQuOSw0My4xIEwgNjUuMCw0My4xIEwgNjUuMSw0My4xIEwgNjUuMiw0My4yIEwgNjUuMyw0My4yIEwgNjUuNCw0My4yIEwgNjUuNSw0My4xIEwgNjUuNSw0My4xIEwgNjUuNiw0My4xIEwgNjUuNyw0My4wIEwgNjUuOCw0My4wIEwgNjUuOSw0Mi45IEwgNjYuMCw0Mi45IEwgNjYuMSw0Mi44IEwgNjYuMSw0Mi43IEwgNjYuMiw0Mi43IEwgNjYuMiw0Mi42IEwgNjYuMyw0Mi41IEwgNjYuMyw0Mi40IEwgNjYuNCw0Mi4zIEwgNjYuNCw0Mi4yIEwgNjYuNCw0Mi4xIEwgNjYuNCw0Mi4wIEwgNjUuMiw0Mi4wIFonIHN0eWxlPSdmaWxsOiMwMDAwMDA7ZmlsbC1ydWxlOmV2ZW5vZGQ7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOiMwMDAwMDA7c3Ryb2tlLXdpZHRoOjAuMHB4O3N0cm9rZS1saW5lY2FwOmJ1dHQ7c3Ryb2tlLWxpbmVqb2luOm1pdGVyO3N0cm9rZS1vcGFjaXR5OjE7JyAvPgo8L3N2Zz4K [Cu] RYGMFSIKBFXOCRUHFFFAOYSAN 0.000 claims description 9
 229910052802 copper Inorganic materials 0.000 claims description 9
 239000010949 copper Substances 0.000 claims description 9
 238000010438 heat treatment Methods 0.000 claims description 9
 238000000034 method Methods 0.000 claims description 9
 230000005611 electricity Effects 0.000 claims description 7
 230000003203 everyday Effects 0.000 claims description 7
 230000000694 effects Effects 0.000 claims description 5
 101710031899 moon Proteins 0.000 claims description 4
 238000004364 calculation method Methods 0.000 claims description 3
 230000035945 sensitivity Effects 0.000 claims description 3
 230000002964 excitative Effects 0.000 claims 2
 238000004088 simulation Methods 0.000 abstract description 12
 230000015572 biosynthetic process Effects 0.000 abstract description 3
 238000003786 synthesis reaction Methods 0.000 abstract description 3
 230000002194 synthesizing Effects 0.000 abstract description 3
 238000010586 diagram Methods 0.000 description 8
 230000002262 irrigation Effects 0.000 description 7
 238000003973 irrigation Methods 0.000 description 7
 238000009114 investigational therapy Methods 0.000 description 4
 230000036962 time dependent Effects 0.000 description 4
 238000004458 analytical method Methods 0.000 description 3
 239000007787 solid Substances 0.000 description 3
 241000196324 Embryophyta Species 0.000 description 1
 230000005540 biological transmission Effects 0.000 description 1
 239000000470 constituent Substances 0.000 description 1
 238000005516 engineering process Methods 0.000 description 1
 230000005284 excitation Effects 0.000 description 1
 238000007710 freezing Methods 0.000 description 1
 230000005484 gravity Effects 0.000 description 1
 239000000463 material Substances 0.000 description 1
 230000001932 seasonal Effects 0.000 description 1
 238000000926 separation method Methods 0.000 description 1
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
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+constantimpedance 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 Realtime 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 timedependent of synthetic load constituent is the essential reason that its part throttle characteristics has stochastic timedependent.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 Orecontrolling Role) can provide the daily load curve of realtime 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 Realtime modeling set, thus
Can fundamentally overcome the inherent shortcoming of conventional statistics synthesis load modeling.
Summary of the invention
In order to overcome abovementioned 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 11: decompose rural area load curve；
Step 111: 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 112: 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 113: calculate this transformer station's load average daily load curve monthly and average monthly load curve；
Step 114: 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 115: 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 116: 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 117: determine the big load curve of this month, and calculate the mean value of datum curve load；
Step 12: calculate the ratio of total amount shared by every kind of load in the load of rural area；
Step 121: calculate each transformer station monthly average load value monthly and specify month to specify the average load value of day；
Step 122: calculate and specify month, specify day, the load value in appointment moment；
Step 123: calculate and include irrigating load, cooling load and the separate section load proportion k of heating load_{p0}And base
Quasiload proportion k_{f}。
In described step 113, 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 115, 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 117, 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 121, 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 122, 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 21) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate；
If 22) 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 123, calculate separate section load proportion k_{p0}, by segregational load value divided by the appointment moon calculated
Part, appointment day, the load value in appointment moment；Then reference load proportion k_{f}For 1k_{p0}。
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=P_{0}[a×(V/V_{0})^{2}+b×(V/V_{0})+c] (1)
Q=Q_{0}[α×(V/V_{0})^{2}+β×(V/V_{0})+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the realtime voltage of static load, V_{0}For
The rated voltage of static load, P_{0}And Q_{0}It is respectively and represents at V_{0}The 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, P_{0}With α, β, γ, Q_{0}Equivalence, to multinomial static load
The equivalence of model is based on the sensitivity to static load terminal voltage of the static load power, has
Wherein, n is static load number, P_{i}And Q_{i}It is respectively the meritorious of ith static load and reactive power,With
It is respectively ith static load active power relative to the partial differential of voltage and ith static load reactive power relative to voltage
Partial differential；
Work as V=V_{0}Time, have:
Wherein, P_{0i}And Q_{0i}It is respectively initial active power and the reactive power of ith static load；
Wherein, a_{i}、b_{i}And c_{i}It is the multinomial static load model active power coefficient of ith static load；
Wherein, α_{i}、β_{i}And γ_{i}It is the multinomial static load model reactive power coefficient of ith static load.
Described step 2 specifically includes following steps:
Step A: calculate all motor always absorbs activepower P_{Σ}With reactive power Q_{Σ}, total electromagnetic power P_{Σem}, total rotor
Winding copper loss P_{Σcu2}With total maximum electromagnetic power P_{Σem_max}, specifically have:
Wherein P_{j}、Q_{j}、P_{emj}、P_{cu2j}And P_{em_maxj}Respectively 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_{Σcu1}And 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 motor_{emt_max}, make P_{emt_max}=P_{Σem_max}；
Step C: calculate the stator winding resistance R of equivalent motor_{s}；
First calculate the stator winding phase current of equivalent motorHave
Wherein,Set end voltage for motor；
The then stator winding resistance R of equivalent motor_{s}It is expressed as:
Step D: calculate the equivalent impedance Z of equivalent machine model_{deq}:
And have
R_{deq}=real(Z_{deq}) (19)
X_{deq}=imag(Z_{deq}) (20)
Wherein, R_{deq}And X_{deq}For corresponding substitutional resistance and equivalent reactance；
Step E: calculate stator winding leakage reactance X of equivalent motor_{s}With rotor windings leakage reactance X_{r}, have
Step F: through iteration stator winding leakage reactance X to equivalent motor_{s}With rotor windings leakage reactance X_{r}It is modified；
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor；
According to calculated R_{s}、X_{s}、X_{r}And Z_{deq}If, K_{r}=R_{deq}R_{s}, K_{x}=X_{deq}X_{s}, the most equivalent electronic
The rotor windings resistance R of machine_{r}Reactance X excitatory with equivalence_{m}It is expressed as:
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 iteration_{emt_maxk}, it is expressed as
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motor_{dp}It is expressed as
And have, R_{dp}=real(Z_{dp}) and X_{dp}=imag(Z_{dp}), R_{dp}And X_{dp}It is respectively Thevenin's equivalence resistance and equivalence electricity
Anti；
Equivalent motor produces the condition of maximum electromagnetic power:
Wherein, R_{pm}The Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondence_{m}For critical slippage；
Thevenin's equivalence opencircuit voltageFor
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
2) the maximum electromagnetic power of equivalent motor is revised；
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic power_{maxk}, it is expressed as:
The most revised equivalent motor maximum electromagnetic power P_{emt_max}It is expressed as:
P_{emt_max}=τ_{maxk}×P_{Σem_max}(30)
With P_{emt_max}With P_{emt_maxk}The iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
WithFor iteration convergence standard, ifThen recalculate X_{s}、R_{r}、X_{r}And X_{m}；
Step I: calculate equivalent inertia time constant；
Equivalent inertia time constant H is expressed as:
Wherein P_{nj}And H_{j}Rated power and inertia time constant for jth platform equivalence motor.
In described step 2, distribution network system impedance meter is shown as:
Wherein, Z_{eq}For distribution network system impedance；u_{f}Represent busbar voltage, Z_{f}Indication transformer and distribution line impedance；I_{l}Table
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 reactiveload compensation data.
Compared with prior art, the beneficial effects of the present invention is:
1. this method electricallybased 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 realtime 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 offline 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 realtime 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 realtime 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. 1Fig. 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 (617 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 abovementioned 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 abovementioned 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 111: 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 112: 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 113: 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 114: 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 115: 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 116: 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 117: 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 58 be Luohe become 35kV rural area load 1,5,7, November decompose after load curve.
As shown in accompanying drawing 14 and accompanying drawing 58, 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 121: 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 122: 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 21) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate；
If 22) 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 123: calculate and include irrigating load, cooling load and the separate section load proportion k of heating load_{p0}And base
Quasiload proportion k_{f}。
Calculate separate section load proportion k_{p0}, will segregational load value divided by the appointment month calculated, specify day, appointment
The load value in moment；Then reference load proportion k_{f}For 1k_{p0}。
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=P_{0}[a×(V/V_{0})^{2}+b×(V/V_{0})+c] (1)
Q=Q_{0}[α×(V/V_{0})^{2}+β×(V/V_{0})+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the realtime voltage of static load, V_{0}For
The rated voltage of static load, P_{0}And Q_{0}It is respectively and represents at V_{0}The 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, P_{0}With α, β, γ, Q_{0}Equivalence, 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.
Wherein, n is static load number, P_{i}And Q_{i}It is respectively the meritorious of ith static load and reactive power,With
It is respectively ith static load active power relative to the partial differential of voltage and ith static load reactive power relative to voltage
Partial differential；
Work as V=V_{0}Time, have:
Wherein, P_{0i}And Q_{0i}It is respectively initial active power and the reactive power of ith static load；
Wherein, a_{i}、b_{i}And c_{i}It is the multinomial static load model active power coefficient of ith static load；
Wherein, α_{i}、β_{i}And γ_{i}It is the multinomial static load model reactive power coefficient of ith static load.
Described step 2 specifically includes following steps:
Step A: calculate all motor always absorbs activepower P_{Σ}With reactive power Q_{Σ}, total electromagnetic power P_{Σem}, total rotor
Winding copper loss P_{Σcu2}With total maximum electromagnetic power P_{Σem_max}, specifically have:
Wherein P_{j}、Q_{j}、P_{emj}、P_{cu2j}And P_{em_maxj}Respectively 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_{Σcu1}And 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 motor_{emt_max}, make P_{emt_max}=P_{Σem_max}；
Step C: calculate the stator winding resistance R of equivalent motor_{s}；
First calculate the stator winding phase current of equivalent motorHave
Wherein,Set end voltage for motor；
The then stator winding resistance R of equivalent motor_{s}It is expressed as:
Step D: calculate the equivalent impedance Z of equivalent machine model_{deq}:
And have
R_{deq}=real(Z_{deq}) (19)
X_{deq}=imag(Z_{deq}) (20)
Wherein, R_{deq}And X_{deq}For corresponding substitutional resistance and equivalent reactance；
Step E: calculate stator winding leakage reactance X of equivalent motor_{s}With rotor windings leakage reactance X_{r}, have
Step F: through iteration stator winding leakage reactance X to equivalent motor_{s}With rotor windings leakage reactance X_{r}It is modified；
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor；
According to calculated R_{s}、X_{s}、X_{r}And Z_{deq}If, K_{r}=R_{deq}R_{s}, K_{x}=X_{deq}X_{s}, the most equivalent electronic
The rotor windings resistance R of machine_{r}Reactance X excitatory with equivalence_{m}It is expressed as:
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 iteration_{emt_maxk}, it is expressed as
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motor_{dp}It is expressed as
And have, R_{dp}=real(Z_{dp}) and X_{dp}=imag(Z_{dp}), R_{dp}And X_{dp}It is respectively Thevenin's equivalence resistance and equivalence electricity
Anti；
Equivalent motor produces the condition of maximum electromagnetic power:
Wherein, R_{pm}The Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondence_{m}For critical slippage；
Thevenin's equivalence opencircuit voltageFor
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
2) the maximum electromagnetic power of equivalent motor is revised；
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic power_{maxk}, it is expressed as:
The most revised equivalent motor maximum electromagnetic power P_{emt_max}It is expressed as:
P_{emt_max}=τ_{maxk}×P_{Σem_max}(30)
With P_{emt_max}With P_{emt_maxk}The iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
WithFor iteration convergence standard, ifThen recalculate X_{s}、R_{r}、X_{r}And X_{m}；
Step I: calculate equivalent inertia time constant；
Equivalent inertia time constant H is expressed as:
Wherein P_{nj}And H_{j}Rated power and inertia time constant for jth platform equivalence motor.
In described step 2, distribution network system impedance meter is shown as:
Wherein, Z_{eq}For distribution network system impedance；u_{f}Represent busbar voltage, Z_{f}Indication transformer and distribution line impedance；I_{l}Table
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 reactiveload compensation data.
Accompanying drawing 1012 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 abovedescribed 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 11: decompose rural area load curve；
Step 111: 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 112: 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 113: calculate this transformer station's load average daily load curve monthly and average monthly load curve；
Step 114: 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 115: 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 116: 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 117: determine the big load curve of this month, and calculate the mean value of datum curve load；
Step 12: calculate the ratio of total amount shared by every kind of load in the load of rural area；
Step 121: calculate each transformer station monthly average load value monthly and specify month to specify the average load value of day；
Step 122: calculate and specify month, specify day, the load value in appointment moment；
Step 123: calculate and include irrigating load, cooling load and the separate section load proportion k of heating load_{p0}And base
Quasiload proportion k_{f}。
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 113, 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 115, 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 117, 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 121, 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 122, 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 21) the daily load curve mean value of this month is less than this definite value, Smaller load curve is used to calculate；
If 22) 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 123, calculate separate section load proportion k_{p0}, by segregational load value divided by the appointment moon calculated
Part, appointment day, the load value in appointment moment；Then reference load proportion k_{f}For 1k_{p0}。
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=P_{0}[a×(V/V_{0})^{2}+b×(V/V_{0})+c] (1)
Q=Q_{0}[α×(V/V_{0})^{2}+β×(V/V_{0})+γ] (2)
Wherein, P and Q is respectively active power and the reactive power of static load, and V is the realtime voltage of static load, V_{0}For
The rated voltage of static load, P_{0}And Q_{0}It is respectively and represents at V_{0}The 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, P_{0}With α, β, γ, Q_{0}Equivalence, to multinomial static load
The equivalence of model is based on the sensitivity to static load terminal voltage of the static load power, has
Wherein, n is static load number, P_{i}And Q_{i}It is respectively the meritorious of ith static load and reactive power,With
It is respectively ith static load active power relative to the partial differential of voltage and ith static load reactive power relative to voltage
Partial differential；
Work as V=V_{0}Time, have:
Wherein, P_{0i}And Q_{0i}It is respectively initial active power and the reactive power of ith static load；
Wherein, a_{i}、b_{i}And c_{i}It is the multinomial static load model active power coefficient of ith static load；
Wherein, α_{i}、β_{i}And γ_{i}It is the multinomial static load model reactive power coefficient of ith 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 activepower P_{∑}With reactive power Q_{∑}, total electromagnetic power P_{∑em}, total rotor
Winding copper loss P_{∑cu2}With total maximum electromagnetic power P_{∑em_max}, specifically have:
Wherein P_{j}、Q_{j}、P_{emj}、P_{cu2j}And P_{em}__{maxj}Respectively 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_{∑cu1}And 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 motor_{emt_max}, make P_{emt_max}=P_{∑em_max}；
Step C: calculate the stator winding resistance R of equivalent motor_{s}；
First calculate the stator winding phase current of equivalent motorHave
Wherein,Set end voltage for motor；
The then stator winding resistance R of equivalent motor_{s}It is expressed as:
Step D: calculate the equivalent impedance Z of equivalent machine model_{deq}:
And have
R_{deq}=real (Z_{deq}) (19)
X_{deq}=imag (Z_{deq}) (20)
Wherein, R_{deq}And X_{deq}For corresponding substitutional resistance and equivalent reactance；
Step E: calculate stator winding leakage reactance X of equivalent motor_{s}With rotor windings leakage reactance X_{r}, have
Step F: through iteration stator winding leakage reactance X to equivalent motor_{s}With rotor windings leakage reactance X_{r}It is modified；
Step G: calculate rotor windings resistance and the excitatory reactance of equivalence of equivalent motor；
According to calculated R_{s}、X_{s}、X_{r}And Z_{deq}If, K_{r}=R_{deq}R_{s}, K_{x}=X_{deq}X_{s}, the most equivalent electronic
The rotor windings resistance R of machine_{r}Reactance X excitatory with equivalence_{m}It is expressed as:
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 iteration_{emt_maxk}, it is expressed as
Based on Thevenin's theorem, the Thevenin's equivalence impedance Z of equivalent motor_{dp}It is expressed as
And have, R_{dp}=real (Z_{dp}) and X_{dp}=imag (Z_{dp}), R_{dp}And X_{dp}It is respectively Thevenin's equivalence resistance and equivalence electricity
Anti；
Equivalent motor produces the condition of maximum electromagnetic power:
Wherein, R_{pm}The Thevenin's equivalence impedance magnitude of maximum electromagnetic power, S is produced for correspondence_{m}For critical slippage；
Thevenin's equivalence opencircuit voltageFor
Then, the maximum electromagnetic power of the Thevenin's equivalence motor in kth time iteration is represented by
2) the maximum electromagnetic power of equivalent motor is revised；
Calculate the correction factor τ of kth time iteration medium value motor maximum electromagnetic power_{maxk}, it is expressed as:
The most revised equivalent motor maximum electromagnetic power P_{emt_max}It is expressed as:
P_{emt_max}=τ_{maxk}×P_{∑em_max} (30)
With P_{emt_max}With P_{emt_maxk}The iteration error that absolute error is equivalent motor maximum electromagnetic powerHave
WithFor iteration convergence standard, ifThen recalculate X_{s}、R_{r}、X_{r}And X_{m}；
Step I: calculate equivalent inertia time constant；
Equivalent inertia time constant H is expressed as:
Wherein P_{nj}And H_{j}Rated 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:
Wherein, Z_{eq}For distribution network system impedance；u_{f}Represent busbar voltage, Z_{f}Indication transformer and distribution line impedance；I_{l}Table
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 reactiveload compensation data.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201310500676.2A CN103632031B (en)  20131022  20131022  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)  20131022  20131022  A kind of rural area based on load curve decomposition load type load modeling method 
Publications (2)
Publication Number  Publication Date 

CN103632031A CN103632031A (en)  20140312 
CN103632031B true CN103632031B (en)  20160831 
Family
ID=50213070
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201310500676.2A Active CN103632031B (en)  20131022  20131022  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)
Publication number  Priority date  Publication date  Assignee  Title 

CN104348154B (en) *  20141122  20160824  郏县供电有限责任公司  The dispatching method of a kind of power distribution network and device 
CN105224812B (en) *  20151021  20181009  中国电力科学研究院  A kind of static load frequency factor polymerization in load model 
CN107453346B (en) *  20160530  20220318  中国电力科学研究院  Load curve time interval division method suitable for power distribution network reconstruction 
CN107069717B (en) *  20170511  20200515  国家电网公司  Power distribution network load modeling method 
CN107392391A (en) *  20170810  20171124  天津求实智源科技有限公司  A kind of electric load electricity consumption decomposition method based on deep learning 
CN114036762A (en) *  20211115  20220211  贵州电网有限责任公司  System and method for predicting load of power distribution network based on proportionality coefficient method 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN101777765A (en) *  20100127  20100714  中国电力科学研究院  Online load simulation method of power system 
CN101789598A (en) *  20100305  20100728  湖北省电力试验研究院  Power system load modelling method 
CN102521674A (en) *  20111207  20120627  福建省电力有限公司  Computer generating method for shortterm load forecasting in summer in power consumer electricityutilization information acquisition system 
CN103279803A (en) *  20130427  20130904  深圳供电局有限公司  Load modeling method and system based on comprehensive information theory and modern interior point theory 

2013
 20131022 CN CN201310500676.2A patent/CN103632031B/en active Active
Patent Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

CN101777765A (en) *  20100127  20100714  中国电力科学研究院  Online load simulation method of power system 
CN101789598A (en) *  20100305  20100728  湖北省电力试验研究院  Power system load modelling method 
CN102521674A (en) *  20111207  20120627  福建省电力有限公司  Computer generating method for shortterm load forecasting in summer in power consumer electricityutilization information acquisition system 
CN103279803A (en) *  20130427  20130904  深圳供电局有限公司  Load modeling method and system based on comprehensive information theory and modern interior point theory 
NonPatent Citations (3)
Title 

基于日负荷曲线的负荷分类和综合建模;鞠平 等;《电力系统自动化》;20060815;第30卷(第16期);第69页 * 
综合负荷模型参数的深化研究及适应性分析;汤涌 等;《电网技术》;20100228;第34卷(第2期);第5763页 * 
考虑配电网络的综合负荷模型;汤涌 等;《电网技术》;20070331;第31卷(第5期);第3438页 * 
Also Published As
Publication number  Publication date 

CN103632031A (en)  20140312 
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)  Online 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 lowvoltage 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)  Online analytical method for directcurrent transmission loss  
CN109494719B (en)  Hierarchical impedance analysis method for mediumlow voltage hybrid power distribution network  
CN102738791B (en)  Online static security analysis method based on topological comparison  
CN103679551A (en)  Method for online calculation of power flow of power distribution network  
CN103123715B (en)  220KV transformer substation electricityusing industry composition proportion online realtime 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 lowvoltage cause big data analysis method and system  
CN107257130A (en)  The lowvoltage network loss computing method of decoupling is measured based on region  
CN104009469B (en)  A kind of lowvoltage distribution network platform district network topology isopleth map collecting method  
CN109286188B (en)  10kV power distribution network theoretical line loss calculation method based on multisource 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 