CN106300337A - Power distribution network partition method based on dynamic Decomposition method - Google Patents
Power distribution network partition method based on dynamic Decomposition method Download PDFInfo
- Publication number
- CN106300337A CN106300337A CN201610645531.5A CN201610645531A CN106300337A CN 106300337 A CN106300337 A CN 106300337A CN 201610645531 A CN201610645531 A CN 201610645531A CN 106300337 A CN106300337 A CN 106300337A
- Authority
- CN
- China
- Prior art keywords
- matrix
- distribution network
- power distribution
- threshold value
- predetermined threshold
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005192 partition Methods 0.000 title claims abstract description 20
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 82
- 230000001174 ascending effect Effects 0.000 claims description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
The invention provides power distribution network partition method based on dynamic Decomposition method, belong to field of electrical equipment.Including obtaining the relative distance of power supply unit in power distribution network, build distance matrix based on relative distance, and then matrix of adjusting the distance carries out a series of process and obtains matrix.Choose non-zero diagonal element the most successively, and carry out multidirectional subsequent searches for the non-zero diagonal element chosen, thus obtain the system subdomain matrix corresponding with non-zero diagonal element, and then complete the subregion to power distribution network according to whole system subdomain matrixes.By introducing electrical distance during subregion and dynamic Decomposition method, improve the accuracy of the division result of power distribution network.
Description
Technical field
The invention belongs to field of electrical equipment, particularly to power distribution network partition method based on dynamic Decomposition method.
Background technology
Being faced with under the dual-pressure of current energy source situation and living environment, people start constantly to seek energy and environment and ask
The solution of topic, wherein distributed new generating is owing to using clean reproducible energy as its one-level energy, and having can
Regeneration, the advantage such as pollution-free, thus, the development of distributed generation technology by the extensive concern of countries in the world and is supported energetically.
Along with being incorporated into the power networks of increasing distributed power source (Distributed Generation, DG), owing to joining
The impact of distributed power source in electricity system so that power distribution network flow of power characteristic is become new two-way by traditional one-way flow
Flowing, thus bring many new challenges to the stable operation of distribution system.Especially build the distribution system planning containing DG to join
When putting model, the most accurately the power distribution network containing DG is carried out accurate subregion, have become as problem demanding prompt solution.
Summary of the invention
In order to solve shortcoming and defect present in prior art, present invention provide for improving joining of subregion accuracy
Sub area division method.
In order to reach above-mentioned technical purpose, the invention provides power distribution network partition method based on dynamic Decomposition method, described
Power distribution network partition method, including:
Step one, obtains power distribution network and includes distributed power source relative distance between interior power supply unit, according to phase
Adjust the distance structure distance matrix r, and the element adjusted the distance in matrix r is normalized, and obtains normalized cumulant matrix D；
Step 2, carries out ascending order arrangement by the element in normalized cumulant matrix D, obtains finite sequence N (k), to limited
Element in sequence N (k) carries out initialization process, obtains predetermined threshold value α0, build and distance matrix rijSize identical complete zero
Matrix A；
Step 3, by the element D in normalized cumulant matrix DijWith predetermined threshold value α0Carry out numerical value contrast, if Dij≥
α0, then the elements A in order matrix Aij=1, if Dij＜ α0, then the elements A in order matrix Aij=0, by normalized cumulant matrix D
In whole elements contrast after, obtain matrix A；
Step 4, chooses non-zero diagonal element A from matrix Amn, along first AmnThe direction that row increase scans for, choosing
Go out the elements A being not zeromk, respectively along elements AmkThe direction that row increase and row increases scan for, until at the end of Sou Suo,
Using the element that searches as the first system subdomain matrix；
Step 5, chooses second non-zero diagonal element from matrix A, determines second system by the mode in step 4
Domain matrix, after N whole non-zero diagonal element has been chosen in matrix A, obtains N number of system subdomain matrix, according to the N got
Individual system subdomain matrix, obtains based on predetermined threshold value α0Power distribution network division result.
Optionally, described predetermined threshold value α0Obtaining step be:
According to number of elements k in finite sequence N (k), determine variable middle=ceil [(low+high)/2], its
In, determining initial value low=1, high=k, ceil represent that return is more than or equal to the computing of smallest positive integral in follow-up expression formula；
According to the variable middle got, from finite sequence N (k), choose N (middle) as predetermined threshold value α0's
Value.
Optionally, if the quantity of the system subdomain matrix got is less than preset requirement, predetermined threshold value α is shown0It is bigger than normal,
Make high=middle, and make αi+1=αi-l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
Wherein, dynamic variable d=[N (high)-N (low)]/k.
Optionally, if the quantity of the system subdomain matrix got is higher than preset requirement, predetermined threshold value α is shown0It is less than normal,
Make low=middle, and make αi+1=αi+ l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
Wherein, dynamic variable d=[N (high)-N (low)]/k.
Optionally, the span of i, j, m, n, k, N is natural number.
The technical scheme that the present invention provides has the benefit that
By introducing electrical distance and dynamic Decomposition method during subregion, improve division result accurate of power distribution network
Property.
Accompanying drawing explanation
In order to be illustrated more clearly that technical scheme, the required accompanying drawing used in embodiment being described below
It is briefly described, it should be apparent that, the accompanying drawing in describing below is only some embodiments of the present invention, general for this area
From the point of view of logical technical staff, on the premise of not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the schematic flow sheet of the power distribution network partition method based on dynamic Decomposition method that the present invention provides；
Fig. 2 is the structural representation of the power distribution network comprising distributed power supply unit that the present invention provides.
Detailed description of the invention
Structure and advantage for making the present invention are clearer, make the structure of the present invention further below in conjunction with accompanying drawing
Describe.
Embodiment one
The invention provides power distribution network partition method based on dynamic Decomposition method, described power distribution network partition method, including:
Step one, obtains power distribution network and includes distributed power source relative distance between interior power supply unit, according to phase
Adjust the distance structure distance matrix r, and the element adjusted the distance in matrix r is normalized, and obtains normalized cumulant matrix D；
Step 2, carries out ascending order arrangement by the element in normalized cumulant matrix D, obtains finite sequence N (k), to limited
Element in sequence N (k) carries out initialization process, obtains predetermined threshold value α0, build and distance matrix rijSize identical complete zero
Matrix A；
Step 3, by the element D in normalized cumulant matrix DijWith predetermined threshold value α0Carry out numerical value contrast, if Dij≥
α0, then the elements A in order matrix Aij=1, if Dij＜ α0, then the elements A in order matrix Aij=0, by normalized cumulant matrix D
In whole elements contrast after, obtain matrix A；
Step 4, chooses non-zero diagonal element A from matrix Amn, along first AmnThe direction that row increase scans for, choosing
Go out the elements A being not zeromk, respectively along elements AmkThe direction that row increase and row increases scan for, until at the end of Sou Suo,
Using the element that searches as the first system subdomain matrix；
Step 5, chooses second non-zero diagonal element from matrix A, determines second system by the mode in step 4
Domain matrix, after N whole non-zero diagonal element has been chosen in matrix A, obtains N number of system subdomain matrix, according to the N got
Individual system subdomain matrix, obtains based on predetermined threshold value α0Power distribution network division result.
In force, in order to the power distribution network comprising distributed power source (Distributed Generation, DG) is carried out
Accurately subregion, reduces the impact that power distribution network is brought by DG, and the present invention proposes power distribution network partition method based on dynamic Decomposition method,
Including entering step one to the step shown in step 7.
Relative distance between the DG arrived according to reality measurement in step one, builds distance matrix r, owing to saving in space
The electrical distance of point is symmetrical, it is possible to D matrix changes into corresponding upper triangular matrix, obtains new D matrix, system
Meter D in nonzero element number and press ascending order arrange after, obtain finite sequence N (k).In step 2, first framework one
Individual size is identical with distance matrix r, element be zero matrix A.And then according to the manner of comparison in step 3, matrix A is carried out
Assignment processes.
In matrix A after treatment, choose non-zero diagonal element.First first element in non-zero diagonal element is chosen
Amn, along first AmnThe direction (the most to the right) that row increase scans for, and selects the elements A being not zeromk, respectively along elements Amk's
The direction that row increase and row increases scans for, until at the end of Sou Suo, using the element that searches as the first system subdomain square
Battle array.Choose second element in non-zero diagonal element the most again, obtain second system subdomain matrix in the same manner, with this
Analogize, based on last element in nonzero element, obtain n-th system subdomain matrix, thus utilize electrical distance and move
State decomposition method related notion carries out dynamic partition to power distribution network.
Optionally, described predetermined threshold value α0Obtaining step be:
According to number of elements k in finite sequence N (k), determine variable middle=ceil [(low+high)/2], its
In, determining initial value low=1, high=k, ceil represent that return is more than or equal to the computing of smallest positive integral in follow-up expression formula；
According to the variable middle got, from finite sequence N (k), choose N (middle) as predetermined threshold value α0's
Value.
In force, said process needs use predetermined threshold value α0, obtaining this threshold value needs to combine finite sequence N (k).
And it also requires determine variable middle and dynamic variable d.Variable middle is obtained according to initial value low, the high built
Numerical value, and then finite sequence N (k) will be chosen N (middle) as predetermined threshold value α0Value.
Optionally, if the quantity of the system subdomain matrix got is less than preset requirement, predetermined threshold value α is shown0It is bigger than normal,
Make high=middle, and make αi+1=αi-l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
If the quantity of the system subdomain matrix got is higher than preset requirement, show predetermined threshold value α0Less than normal, make low=
Middle, and make αi+1=αi+ l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
Wherein, dynamic variable d=[N (high)-N (low)]/k.
In force, if according to predetermined threshold value α0The quantity of the system subdomain matrix got is less than preset requirement, then
Show predetermined threshold value α0Setting value there is deviation.
When the quantity of system subdomain matrix is less than preset requirement, show predetermined threshold value α0Bigger than normal, make high=middle, and
And make αi+1=αi-l d, to predetermined threshold value α0It is updated；
If when the quantity of the system subdomain matrix got is higher than preset requirement, showing predetermined threshold value α0Less than normal, make low
=middle, and make αi+1=αi+ l d, to predetermined threshold value α0It is updated.
To predetermined threshold value α0After being updated, repeat the step 2 in preceding method to step 5, in order to obtain and meet
The system subdomain matrix of predetermined number, thus obtain based on predetermined threshold value α0Power distribution network division result.
In unified statement herein, in the present embodiment, the span of i, j, m, n, k, N is natural number, and whole non-negative is whole
Number.
The invention provides power distribution network partition method based on dynamic Decomposition method, including obtaining power supply unit in power distribution network
Relative distance, builds distance matrix based on relative distance, and then matrix of adjusting the distance carries out a series of process and obtains matrix.At matrix
In choose non-zero diagonal element successively, and carry out multidirectional subsequent searches for the non-zero diagonal element chosen, thus obtain
The system subdomain matrix corresponding with non-zero diagonal element, and then complete power distribution network is divided according to whole system subdomain matrixes
District.By utilizing electrical distance and dynamic Decomposition method during subregion, improve the accuracy of the division result of power distribution network.
Emulation experiment
There is distribution net work structure as shown in Figure 2 in practice, DG therein represents the distributed electrical being located in power distribution network
Source.The stain being scattered between DG represents other kinds of power supply unit, and the sequence number of lower section is for marking the sequence number of each unit.
The power distribution network partition method provided according to the present invention, carries out result such as table 1 institute of subregion to the power supply unit in Fig. 2
Show.
| Partition number | Division result |
| Subregion 1 | 1,2,3,18,19,20,21,22,23 |
| Subregion 2 | 4,5,6,7,8,24,25,26,27,28,29,30 |
| Subregion 3 | 9,10,11,12,13,14,15,16,17,31,32 |
Table 1 carries out the result of subregion to power supply unit
Each sequence number in above-described embodiment just to describe, do not represent each parts assemble or use during elder generation
Rear order.
The foregoing is only embodiments of the invention, not in order to limit the present invention, all in the spirit and principles in the present invention
Within, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.
Claims (5)
1. power distribution network partition method based on dynamic Decomposition method, it is characterised in that described power distribution network partition method, including:
Step one, obtains power distribution network and includes distributed power source relative distance between interior power supply unit, according to relatively away from
From building distance matrix r, the element adjusted the distance in matrix r is normalized, and obtains normalized cumulant matrix D；
Step 2, carries out ascending order arrangement by the element in normalized cumulant matrix D, obtains finite sequence N (k), to finite sequence N
K the element in () carries out initialization process, obtain predetermined threshold value α0, build and distance matrix rijFull null matrix A that size is identical；
Step 3, by the element D in normalized cumulant matrix DijWith predetermined threshold value α0Carry out numerical value contrast, if Dij≥α0, then
Elements A in order matrix Aij=1, if Dij＜ α0, then the elements A in order matrix Aij=0, by normalized cumulant matrix D
After all element contrasts, obtain matrix A；
Step 4, chooses non-zero diagonal element A from matrix Amn, along first AmnThe direction that row increase scans for, and selects not
It it is the elements A of zeromk, respectively along elements AmkRow increase and row increase direction scan for, until at the end of Sou Suo, will search
The element that rope arrives is as the first system subdomain matrix；
Step 5, chooses second non-zero diagonal element from matrix A, determines second system subdomain square by the mode in step 4
Battle array, after N whole non-zero diagonal element has been chosen in matrix A, obtains N number of system subdomain matrix, according to the N number of system got
System subdomain matrix, obtains based on predetermined threshold value α0Power distribution network division result.
Power distribution network partition method based on dynamic Decomposition method the most according to claim 1, it is characterised in that described default threshold
Value α0Obtaining step be:
According to number of elements k in finite sequence N (k), determine variable middle=ceil [(low+high)/2], wherein, really
Determining initial value low=1, high=k, ceil represent that return is more than or equal to the computing of smallest positive integral in follow-up expression formula；
According to the variable middle got, from finite sequence N (k), choose N (middle) as predetermined threshold value α0Value.
Power distribution network partition method based on dynamic Decomposition method the most according to claim 2, it is characterised in that:
If the quantity of the system subdomain matrix got is less than preset requirement, show predetermined threshold value α0Bigger than normal, make high=
Middle, and make αi+1=αi-l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
Wherein, dynamic variable d=[N (high)-N (low)]/k.
Power distribution network partition method based on dynamic Decomposition method the most according to claim 2, it is characterised in that:
If the quantity of the system subdomain matrix got is higher than preset requirement, show predetermined threshold value α0Less than normal, make low=
Middle, and make αi+1=αi+ l d, to predetermined threshold value α0It is updated；
Based on predetermined threshold value α after updating0, repeat step 2 to step 5；
Wherein, dynamic variable d=[N (high)-N (low)]/k.
5. according to the power distribution network partition method based on dynamic Decomposition method described in any one of Claims 1-4, it is characterised in that:
The span of i, j, m, n, k, N is natural number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610645531.5A CN106300337A (en) | 2016-08-09 | 2016-08-09 | Power distribution network partition method based on dynamic Decomposition method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610645531.5A CN106300337A (en) | 2016-08-09 | 2016-08-09 | Power distribution network partition method based on dynamic Decomposition method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106300337A true CN106300337A (en) | 2017-01-04 |
Family
ID=57667000
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610645531.5A Pending CN106300337A (en) | 2016-08-09 | 2016-08-09 | Power distribution network partition method based on dynamic Decomposition method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106300337A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130015710A1 (en) * | 2007-05-17 | 2013-01-17 | Enphase Energy, Inc. | Distributed inverter and intelligent gateway |
| CN103595039A (en) * | 2013-10-17 | 2014-02-19 | 广东电网公司茂名供电局 | Distributed loop closing risk analysis method for power grid system |
| CN104917203A (en) * | 2015-06-26 | 2015-09-16 | 国家电网公司 | Distributed power supply combined power generation system based on geographical division and coordination control method |
-
2016
- 2016-08-09 CN CN201610645531.5A patent/CN106300337A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130015710A1 (en) * | 2007-05-17 | 2013-01-17 | Enphase Energy, Inc. | Distributed inverter and intelligent gateway |
| CN103595039A (en) * | 2013-10-17 | 2014-02-19 | 广东电网公司茂名供电局 | Distributed loop closing risk analysis method for power grid system |
| CN104917203A (en) * | 2015-06-26 | 2015-09-16 | 国家电网公司 | Distributed power supply combined power generation system based on geographical division and coordination control method |
Non-Patent Citations (2)
| Title |
|---|
| 刘增杰，胡同夫: "《MariaDB从零开始学（视频教学版）》", 31 October 2015 * |
| 叶丽红: ""计及分布式电源的配电***综合优化配置方法的研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tolbert et al. | Real-time selective harmonic minimization for multilevel inverters connected to solar panels using artificial neural network angle generation | |
| CN105790279B (en) | Reactive voltage partition method based on spectral clustering | |
| CN105260952B (en) | Photovoltaic plant reliability estimation method based on Markov chain Monte-Carlo method | |
| CN110099017A (en) | The channel estimation methods of mixing quantization system based on deep neural network | |
| CN109873501A (en) | A kind of low-voltage network topology automatic identification method | |
| CN104242306B (en) | Power system self-adaption partition method based on principal component analysis method | |
| CN104638672B (en) | Determining method of photovoltaic transmission power limit considering variable correlation | |
| He et al. | Energy-efficient power allocation with individual and sum power constraints | |
| CN109255514A (en) | A kind of independently-powered capability assessment method of intelligent distribution network subregion | |
| CN105912812A (en) | Method and device for determining shed parameters of post insulator | |
| CN108183754B (en) | Electromagnetic radiation prediction method for GSM base station | |
| CN109615246A (en) | A kind of active distribution network economical operation state determines method | |
| CN106300337A (en) | Power distribution network partition method based on dynamic Decomposition method | |
| CN106356867B (en) | Method for determining node reactive voltage sensitivity | |
| CN107817684A (en) | A kind of carrier rocket quick fault testing policy optimization method | |
| CN106953331A (en) | A kind of Harmonic Power Flow of Power Systems algorithm for considering uncertain and three-phase imbalance | |
| CN116404642A (en) | Distributed power supply cluster division method and device, electronic equipment and storage medium | |
| CN106155177A (en) | The method realizing SPWM waveform based on digital processing unit | |
| CN105117814A (en) | Method for actively power distribution network bi-layer wind-power planning based on improved Cuckoo search algorithm | |
| CN105277852A (en) | Classification and identification method of line conditions of power distribution network | |
| CN105447251B (en) | A kind of verification method based on transaction types excitation | |
| Barari et al. | Analytical evaluation of beam deformation problem using approximate methods | |
| CN105789676B (en) | A kind of electrokinetic cell method for group matching propagated based on message | |
| CN104008305B (en) | For ten million kilowatt of wind power base can power generating wind resource distribution method of estimation | |
| CN106327079A (en) | Method for evaluating reactive optimization control of power distributing network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170125 Address after: 313000 Huzhou economic and Technological Development Zone, Zhejiang, Phoenix Road, No. 777, No. Applicant after: Huzhou Power Supply Company of State Grid Zhejiang Electric Power Company Applicant after: State Grid Corporation of China Applicant after: State net Zhejiang Deqing County electric company Applicant after: Zhejiang University Address before: 313200 Zhejiang city of Huzhou province Deqing County ZTE Wukang Road No. 9 Applicant before: State net Zhejiang Deqing County electric company Applicant before: Zhejiang University |
|
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170104 |
|
| RJ01 | Rejection of invention patent application after publication |