CN107834538B - Multimode electric system redundancy optimization method based on sequence optimization and markov chain - Google Patents

Multimode electric system redundancy optimization method based on sequence optimization and markov chain Download PDF

Info

Publication number
CN107834538B
CN107834538B CN201710944284.3A CN201710944284A CN107834538B CN 107834538 B CN107834538 B CN 107834538B CN 201710944284 A CN201710944284 A CN 201710944284A CN 107834538 B CN107834538 B CN 107834538B
Authority
CN
China
Prior art keywords
reliability
fundamental
parallel
current
electric system
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710944284.3A
Other languages
Chinese (zh)
Other versions
CN107834538A (en
Inventor
胡怡霜
丁一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710944284.3A priority Critical patent/CN107834538B/en
Publication of CN107834538A publication Critical patent/CN107834538A/en
Application granted granted Critical
Publication of CN107834538B publication Critical patent/CN107834538B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70Smart grids as climate change mitigation technology in the energy generation sector
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Power Engineering (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of multimode electric system redundancy optimization methods based on sequence optimization and markov chain.Classification is carried out to all systems and finds out fundamental system, the reliability for finding out fundamental system is calculated using markov chain and its general generating function;The rough reliability for obtaining remaining system is calculated by the first modification method, is carried out processing for the rough reliability of all multimode electric system using sequence optimization algorithm BP and is carried out first time screening;Reliability based on fundamental system calculates the accurate reliability for obtaining remaining system by the second modification method, carries out programmed screening using sequence optimization algorithm BP;Optimal multimode electric system is finally acquired using redundancy optimization algorithm.It is greatly to reduce that the present invention, which seeks the time used, and there is screening process more accurate accuracy to be more applicable for reality suitable for the redundancy optimization algorithm of system.

Description

Multimode electric system redundancy optimization method based on sequence optimization and markov chain
Technical field
The present invention relates to a kind of electric power system optimization processing methods, are related in multimode Power System Reliability Analysis Based on sequence optimization and markov chain multimode electric system redundancy optimization method.
Background technique
Reliability engineering is grown up first from aerospace industry and electronics industry after World War II, at present Aerospace, electronics, chemical industry, many industrial departments such as machinery are penetrated into.Reliability engineering penetrates into power industry and electrical equipment Manufacturing industry starts from middle 1960s, and what is developed later is very fast.The function of electric system as far as possible may be used to user Qualified electric energy is economically provided by ground, its reliability may be defined as providing a user up-to-standard, the energy of continuous electric energy Power, this ability are usually indicated with probability.It is so-called up-to-standard, just refer to that the frequency of electric energy and voltage preferably must be held in regulation model Within enclosing.
Model in Reliability Evaluation of Power Systems is to calculate analysis possible breakdown shape probability of state and consequence, show that reflection system is reliable Property a series of horizontal indexs.However, at one in several hundred or even thousands of a elements real systems, it may occur however that event The enormous amount of barrier state.It, can not be to all possible in actual assessment due to calculating the limitation of time and computing resource Malfunction is assessed.Therefore, State enumeration method, which only screens, contributes big malfunction to assess system reliability.Most Common selection method is off failure tuple, that is, selects 2 weights or the 3 following malfunctions of weight, ignore the malfunction of more Gao Chong. The advantages of this method is the sum of selected shape probability of state close to 1, and negligible amounts.But in systems in practice, due to element Outage probability it is different, some high heavy failures can be bigger than the probability of happening of low heavy failure.With the IEEE-RTS with 71 elements For system, when element uses 2 state model, consider that the system mode quantity of N-3 is 57226, the sum of probability is 0.95110503.Actually, biggish preceding 57226 states of probability include 16786 0~3 weight malfunctions and 40440 4 The heavy malfunction in weight~6, the sum of probability are 0.98976138.Probability of these height weight malfunctions are big and consequence is serious, to being The reliability effect of system is very big, and these maximum probability height weight failure can be neglected by carrying out system states filter by cut-off failure tuple.
From the foregoing, it will be observed that in Power System Reliability Analysis, the quantity of state selection mostly on, retain or delete, all can Very big influence is caused to final result, this mean that the Power System Reliability Analysis in multistate model be one very It is necessary to research direction.
In existing multimode Power System Reliability redundancy optimization calculation method, most common is exactly GA algorithm, if Wish to seek having the system i.e. reliability of optimum structure to meet the requirements in a series of multimode electric system, economy is most Excellent system calculates the systematic reliability of institute, when this method can consume very long calculating often using the method for exhaustion Between, a few houres even several days.
Summary of the invention
In order to solve the problems, such as background technique, the invention proposes a kind of more shapes based on sequence optimization and markov chain State electric system redundancy optimization method.The method of the present invention improves conventional method from time and precision, applies to a series of In multimode electric system, it can seek obtaining the best multistate system of the cost lowest economic met under reliability requires, And the used time is far smaller than pervious traditional algorithm.
As shown in Figure 1, technical scheme is as follows:
Step 1: classify to all multimode electric system, found out in every one kind one it is most representational Fundamental system, and the reliability for finding out fundamental system is calculated using markov chain and its general generating function;
Step 2: the difference of Series Parallel Elements element in multimode electric system is considered, based on the reliable of fundamental system Degree calculates the rough reliability for obtaining each multimode electric system in addition to fundamental system, benefit by the first modification method Processing is carried out for the rough reliability of all multimode electric system with sequence optimization algorithm BP and carries out first time screening, passes through sieve Choosing again reduces calculation amount;
Step 3: the difference of inside and outside Series Parallel Elements element in multimode electric system is considered, based on fundamental system Reliability is calculated by the second modification method and obtains the accurate reliable of each multimode electric system in addition to fundamental system Degree carries out programmed screening using sequence optimization algorithm BP;
Step 4: being acquired most for all multimode electric system that third step is screened using redundancy optimization algorithm Excellent multimode electric system.
The system quantity that the present invention greatly reduces calculating by carrying out category filter the step of second step and third step, and And improve the accuracy that reliability screening calculates.
Multistate system of the present invention is defined as: system and its element may all show multiple performance levels, this germline It is referred to as multistate system.
The first step specifically:
1.1) all multimode electric system are classified, to be parallel with add ons and phase in multimode electric system The sum of concatenated major component is identical more by the sum for the major component for being parallel with add ons and being connected in series as classification foundation State electric system is classified as a classification;
1.2) be classified as same classification from fundamental system is successively extracted, fundamental system refers to the element (packet in same classification Include add ons and major component) the least system of sum;
It is subsequent to be handled again for the other systems in same classification in addition to fundamental system, in addition to fundamental system The series-parallel element of other systems and the difference of fundamental system be all carried out on the basis of fundamental system it is in parallel or altogether because The change of failure structure.
1.3) reliability for finding out fundamental system in each classification is calculated with markov chain and its general generating function.
The second step specifically:
In front on the basis of the fundamental system and reliability of the affiliated class of the first step, for the residue system in same classification System, remaining system refer to other multimode electric system in addition to fundamental system, successively carry out reliability in the following ways It calculates:
2.1) firstly, the reliability for carrying out additional parallel element is updated and calculated:
2.1.a) if remaining system is compared with fundamental system, increase has a parallel element, then reliability calculating after updating Formula are as follows:
A '=1- (1-A) (1-R)
Wherein, A is current reliability, and R is the reliability of parallel element, and A ' is reliability after updating;
2.1.b) if remaining system is compared with fundamental system, increase the add ons for there are multiple parallel connections, repeat step 2.1.a) successively all increased parallel elements are iterated to calculate in the same manner, when initial calculation, current reliability A was basic The reliability of system, reliability is as next increased after the update obtained later using current increased parallel element alignment processing Current reliability A when parallel element alignment processing, to obtain the first intermediate reliability of remaining system;
2.2) then, (fundamental system thinks not containing in the present invention for the reliability update calculating of progress common cause failure structure Common cause failure structure):
2.2.a) if remaining system is compared with fundamental system, wherein there is one group of common cause failure structure more, then it can after updating By spending the product for being current reliability B multiplied by all element reliabilitys of existing common cause failure structure;
2.2.b) if remaining system is compared with fundamental system, there are multiple groups common cause failure structures, repeat step 2.2.a) according to It is secondary that all common cause failure structures are iterated to calculate in the same manner,
Current reliability B is the first intermediate reliability that step 2.1) obtains when initial calculation, later currently altogether because losing Current reliability B when reliability is as next common cause failure structure alignment processing after the update that effect structure alignment processing obtains, To obtain the rough reliability of remaining system;
Common cause failure structure refers to causes two or more units due to common failure in a system While fail.
2.3) by the rough reliability of the reliability of the fundamental system of all classification and the remaining system of all classification according to Descending is arranged, and the first collating sequence is obtained, and descending alignment curve needed for drawing out sequence optimization algorithm;
2.4) the first collating sequence is handled using the blind choosing method BP in sequence optimization algorithm, specifically from first row G system finds out " selection with blind choosing method BP as " the subset G good enough " in blind choosing method BP before selecting in sequence sequence Collect the number s of S ".
In order to improve the accuracy of final result, the present invention is carried out another Calculation of Reliability and is resequenced with third step Specifically which system of the number s of " selection subset S " sought.
The third step specifically:
Third step is similar with second step, will be reliable on the basis of the fundamental system and reliability that select in second step The variable quantity of degree makes corresponding changes according to the inside and outside level of change.
3.1) firstly, the reliability for carrying out additional parallel element is updated and calculated:
3.1.a) if remaining system is compared with fundamental system, increase has an internal parallel element, and internal parallel element is Refer to the other attachment element being in series with again in the add ons that major component is parallel with, then formula of reliability after updating Are as follows:
A '=A+ (R '-AR ') R '
Wherein, A is current reliability, and R ' is the reliability of internal parallel element, and A ' is reliability after updating;
3.1.b) if remaining system is compared with fundamental system, increase has multiple internal parallel elements, repeats step 3.1.a) Successively all increased internal parallel elements are iterated to calculate in the same manner, when initial calculation, current reliability A was fundamental system The reliability of system, reliability is first as next internal parallel after the update obtained later using current internal parallel element alignment processing Current reliability A when part alignment processing, to obtain the second intermediate reliability of remaining system;
3.2) reliability for then, carrying out additional parallel element, which updates, to be calculated:
3.2.a) if remaining system is compared with fundamental system, increase has a parallel connection outside element, and parallel connection outside element is Refer to the other attachment element being parallel on major component or add ons, then formula of reliability after updating are as follows:
A '=1- (1-A) (1-R ")
Wherein, A is current reliability, and R " is the reliability of parallel connection outside element, and A ' is reliability after updating;
The reliability formula of parallel connection outside element and the step 2.1) formula of nonrated parallel element are consistent.
3.2.b) if remaining system is compared with fundamental system, increase has multiple parallel connection outside elements, repeats step 3.2.a) Successively all increased parallel connection outside elements are iterated to calculate in the same manner, when initial calculation, current reliability A was step 3.1) obtain the second intermediate reliability, later using current external parallel element alignment processing obtain update after reliability as Current reliability A when next parallel connection outside element alignment processing, to obtain reliability among the third of remaining system;
3.3) then, (fundamental system thinks not containing in the present invention for the reliability update calculating of progress common cause failure structure Common cause failure structure):
3.3.a) if remaining system is compared with fundamental system, wherein there is one group of common cause failure structure more, then it can after updating By spending the product for being current reliability B multiplied by all element reliabilitys of existing common cause failure structure;
3.3.b) if remaining system is compared with fundamental system, there are multiple groups common cause failure structures, repeat step 3.3.a) according to Secondary to iterate to calculate in the same manner to all common cause failure structures, when initial calculation, current reliability B was what step 3.2) obtained Reliability among third, reliability is as lower altogether because losing after the update obtained later using current common cause failure structure alignment processing Current reliability B when structure alignment processing is imitated, to obtain the accurate reliability of remaining system;
3.4) for the preceding g system selected in step 2.4), arranged with accurate reliability according to descending, and again from For s system as the element in " selection subset S ", s is the number of " selection subset S " that step 2.4) obtains before middle selection.
The s system that step 3.4) obtains herein is obtained after being handled by preceding screening technique twice with high-reliability System, i.e. s system is to meet the system of reliability requirement in all systems, and subsequent step uses in this s system Redundancy optimization algorithm is handled.
Second modification method of the invention not only allows in parallel and common cause failure structure compared to the first modification method It influences, it is also contemplated that the influence of internal parallel and internal common cause failure structure, parallel connection outside and external common cause failure structure, So as to obtain the reliability more accurate compared to the first.
The system quantity that the present invention greatly reduces calculating by carrying out category filter the step of second step and third step, and And improve the accuracy that reliability screening calculates.
4th step specifically:
For all multimode electric system that third step is screened, carried out using the redundancy optimization algorithm of following formula Processing, obtains the system of wherein cost minimization as optimal system by min function:
Min C=∑ ci*ni
Wherein, c is the cost of each element in single multimode electric system, and n is the total quantity of element, and C is single more The cost of state electric system.
The beneficial effects of the present invention are:
Scheme compared with the existing technology, for optimal system, (meet that reliability requires minimum cost is the present invention System) time used of seeking reduce significantly.
Scheme compared with the existing technology, screening process proposed by the present invention have more accurate accuracy, ensure that height The accuracy of reliability screening system.
Scheme compared with the existing technology, Monte Carlo are more applicable for system compared to markov chain, and the time used compares horse Family name's chain is few, and the present invention is suitable for the redundancy optimization algorithm of system, is more applicable for the actual conditions of reality.
Detailed description of the invention
Fig. 1 is the logic diagram of the method for the present invention.
Fig. 2 is the multimode electric system that the structure type number of the embodiment of the present invention is 1-1.
Fig. 3 is the multimode electric system that the structure type number of the embodiment of the present invention is 1-2.
Fig. 4 is the multimode electric system that the structure type number of the embodiment of the present invention is 1-3.
Fig. 5 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-1.
Fig. 6 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-2.
Fig. 7 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-3.
Fig. 8 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-4.
Fig. 9 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-5.
Figure 10 is the multimode electric system that the structure type number of the embodiment of the present invention is 2-6.
Figure 11 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-1.
Figure 12 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-2.
Figure 13 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-3.
Figure 14 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-4.
Figure 15 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-5.
Figure 16 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-6.
Figure 17 is the multimode electric system that the structure type number of the embodiment of the present invention is 3-7.
Figure 18 is the multimode electric system that the structure type number of the embodiment of the present invention is 4-1.
Figure 19 is the multimode electric system that the structure type number of the embodiment of the present invention is 4-2.
Figure 20 is the multimode electric system that the structure type number of the embodiment of the present invention is 4-3.
Figure 21 is the multimode electric system that the structure type number of the embodiment of the present invention is 4-4.
Figure 22 is the multimode electric system that the structure type number of the embodiment of the present invention is 4-5.
Figure 23 is the multimode electric system that the structure type number of the embodiment of the present invention is 5-1.
Figure 24 is the multimode electric system that the structure type number of the embodiment of the present invention is 5-2.
Figure 25 is the multimode electric system that the structure type number of the embodiment of the present invention is 5-3.
Figure 26 is the multimode electric system that the structure type number of the embodiment of the present invention is 5-4.
Figure 27 is the multimode electric system that the structure type number of the embodiment of the present invention is 6-1.
Figure 28 is the multimode electric system that the structure type number of the embodiment of the present invention is 6-2.
Figure 29 is the multimode electric system that the structure type number of the embodiment of the present invention is 6-3.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
The embodiment of the present invention is as follows:
Step 1: classify to all multimode electric system, found out in every one kind one it is most representational Fundamental system, and the reliability for finding out fundamental system is calculated using markov chain and its general generating function;
Assuming that there are 28 kinds of system structures, as shown in Fig. 2~Figure 29,6 classes are divided into according to classification one.Assuming that element is identical, And cost is 1.
In each system structure, box represents each force device, and two box series connections represent force device Series connection, two box parallel connections represent the parallel connection of force device, and the box between two two elements represents to be deposited between the two elements In common cause failure structure.
Each system structure just represents a kind of NETWORK STRUCTURE PRESERVING POWER SYSTEM, from power generation, power transformation, and transmission of electricity, distribution and electricity consumption side Topological structure between each class component.The i.e. structure number of series element number of add ons in parallel is had in each structure Preceding digital, each structure are number behind the structure number in the chart sequence number in affiliated class.
The structure type number of Fig. 2~Fig. 4 is respectively 1-1,1-2,1-3, and the structure type number of Fig. 5~Figure 10 is respectively The structure type number of 2-1,2-2,2-3,2-4,2-5,2-6, Figure 11~Figure 17 are respectively 3-1,3-2,3-3,3-4,3-5,3- 6, the structure type number of 3-7, Figure 18~Figure 22 are respectively 4-1,4-2,4-3,4-4,4-5, the structure type of Figure 13~Figure 26 Number is respectively 5-1,5-2,5-3,5-4, and the structure type number of Figure 27~Figure 29 is respectively 6-1,6-2,6-3.
The reliability for defining each element is 0.98, finds out every a kind of fundamental system using markov chain and its general generating function The accurate reliability of system, selects fundamental system for 1-1,2-1,3-1,4-1,5-1,6-1, accurate reliability 0.9036, 0.9216,0.9401,0.9589,0.9780,0.9975.
Step 2: the reliability based on fundamental system, is calculated by the first modification method and is obtained in addition to fundamental system The rough reliability of each multimode electric system, can for all the rough of multimode electric system using sequence optimization algorithm BP Processing, which is carried out, by degree carries out first time screening;
Using the calculation formula of definition, the rough reliability of the remaining system in every one kind is sought, preceding 23 kinds of substitutions sequence is selected In optimization.Corresponding reliability is shown in Table 1.
The rough reliability of 1 system of table
Step 3: the reliability based on fundamental system, is calculated by the second modification method and is obtained in addition to fundamental system The accurate reliability of each multimode electric system carries out programmed screening using sequence optimization algorithm BP;
By the BP rule of sequence optimization algorithm, acquires and select the element number of subset for 15.
Using the formula of definition, the systematic more accurate reliability of institute in 23 systems is sought, the reliable of preceding 15 systems is asked Property is as shown in the table
The more accurate reliability of 2 15 systems of table
Step 4: being acquired most for all multimode electric system that third step is screened using redundancy optimization algorithm Excellent multimode electric system, wherein following table is the cost of 15 systems.
The cost of 3 15 systems of table
The reliability of system 6-2 and 6-3 is high as can be seen from Table 4 and cost is minimum, so system 6-2 and 6-3 are exactly Optimal system.
Embodiment verifying
A) reliability accuracy validation
The accurate reliability of all 28 systems is solved using markov chain
The accurate reliability of the whole systems of table 4
System 6-2 and 6-3 has high-reliability and least cost as can be seen from the above table, meets calculated result.
B) compare and put in order
In terms of accurate reliability, top 10 system is 6-2,6-3,6-1,5-4,5-3,5-2,5-1,4-2,4-5,4-4.
It is 6-2,6-3,6-1 according to the top 10 of rough reliability and the system of more accurate reliability filtered out twice, 5-4,5-3,5-2,4-2,4-5,4-3,3-5,3-4, two arrangements generally remain identical.This demonstrate rough reliabilitys and more smart The accuracy of true reliability screened twice.
C) with only with the time superiority compared with markov chain and its traditional algorithm of general generating function
As can be seen from Table 5, the substantially traditional algorithm volume half of time used in new method, greatly reduces the time
Table 5 calculates the time
It can be seen that the present invention has more accurate accuracy, method high reliablity, and it is big to seek the time used Big reduction, is more suitable for the actual conditions of reality, and obvious technical effects are prominent.

Claims (3)

1. a kind of multimode electric system redundancy optimization method based on sequence optimization and markov chain, it is characterised in that including walking as follows It is rapid:
Step 1: classify to all multimode electric system, found out in every one kind one it is most representational basic System, and the reliability for finding out fundamental system is calculated using markov chain and its general generating function;
Step 2: the reliability based on fundamental system, is calculated by the first modification method and is obtained each of in addition to fundamental system The rough reliability of multimode electric system is directed to the rough reliability of all multimode electric system using sequence optimization algorithm BP It carries out processing and carries out first time screening;
Step 3: the reliability based on fundamental system, is calculated by the second modification method and is obtained each of in addition to fundamental system The accurate reliability of multimode electric system carries out programmed screening using sequence optimization algorithm BP;
Step 4: being acquired using redundancy optimization algorithm optimal for all multimode electric system that third step is screened Multimode electric system;
The second step specifically:
For the remaining system in same classification, remaining system refers to other multimode electric system in addition to fundamental system, Reliability calculating is successively carried out in the following ways:
2.1) firstly, the reliability for carrying out additional parallel element is updated and calculated:
2.1.a) if remaining system is compared with fundamental system, increase has a parallel element, then formula of reliability after updating Are as follows:
A '=1- (1-A) (1-R)
Wherein, A is current reliability, and R is the reliability of parallel element, and A ' is reliability after updating;
2.1.b) if remaining system is compared with fundamental system, increase the add ons for having multiple parallel connections, repeat step 2.1.a) according to It is secondary that all increased parallel elements are iterated to calculate in the same manner, when initial calculation current reliability A be fundamental system can By degree, later using reliability after the update of current increased parallel element alignment processing acquisition as next increased parallel element Current reliability A when alignment processing, to obtain the first intermediate reliability of remaining system;
2.2) reliability for then, carrying out common cause failure structure, which updates, to be calculated:
2.2.a) if remaining system is compared with fundamental system, wherein there is one group of common cause failure structure more, then reliability after updating It is current reliability B multiplied by the product of all element reliabilitys of existing common cause failure structure;
2.2.b) if remaining system is compared with fundamental system, there are multiple groups common cause failure structures, repeat step 2.2.a) it is successively right All common cause failure structures iterate to calculate in the same manner,
Current reliability B is the first intermediate reliability that step 2.1) obtains when initial calculation, later with current common cause failure knot Current reliability B when reliability is as next common cause failure structure alignment processing after the update that structure alignment processing obtains, thus Obtain the rough reliability of remaining system;
2.3) by the rough reliability of the reliability of the fundamental system of all classification and the remaining system of all classification according to descending It is arranged, obtains the first collating sequence, and draw out descending alignment curve;
2.4) the first collating sequence is handled using the blind choosing method BP in sequence optimization algorithm, specifically from the first sequence sequence G system finds out " selection subset S " with blind choosing method BP as " the subset G good enough " in blind choosing method BP before selecting in column Number s;
The third step specifically:
3.1) firstly, the reliability for carrying out additional parallel element is updated and calculated:
3.1.a) if remaining system is compared with fundamental system, increase has an internal parallel element, and internal parallel element refers to The other attachment element being in series with again in the add ons that major component is parallel with, then formula of reliability after updating are as follows:
A '=A+ (R '-AR ') R '
Wherein, A is current reliability, and R ' is the reliability of internal parallel element, and A ' is reliability after updating;
3.1.b) if remaining system is compared with fundamental system, increase has multiple internal parallel elements, repeats step 3.1.a) successively All increased internal parallel elements are iterated to calculate in the same manner, when initial calculation, current reliability A was fundamental system Reliability, later using reliability after the update of current internal parallel element alignment processing acquisition as next internal parallel element pair Current reliability A when should handle, to obtain the second intermediate reliability of remaining system;
3.2) reliability for then, carrying out additional parallel element, which updates, to be calculated:
3.2.a) if remaining system is compared with fundamental system, increase has a parallel connection outside element, and parallel connection outside element refers to The other attachment element being parallel on major component or add ons, then formula of reliability after updating are as follows:
A '=1- (1-A) (1-R ")
Wherein, A is current reliability, and R " is the reliability of parallel connection outside element, and A ' is reliability after updating;
3.2.b) if remaining system is compared with fundamental system, increase has multiple parallel connection outside elements, repeats step 3.2.a) successively All increased parallel connection outside elements are iterated to calculate in the same manner, when initial calculation, current reliability A was that step 3.1) obtains The second intermediate reliability, reliability is as next outer after the update obtained later using current external parallel element alignment processing Current reliability A when portion's parallel element alignment processing, to obtain reliability among the third of remaining system;
3.3) reliability for then, carrying out common cause failure structure, which updates, to be calculated:
3.3.a) if remaining system is compared with fundamental system, wherein there is one group of common cause failure structure more, then reliability after updating It is current reliability B multiplied by the product of all element reliabilitys of existing common cause failure structure;
3.3.b) if remaining system is compared with fundamental system, there are multiple groups common cause failure structures, repeat step 3.3.a) it is successively right All common cause failure structures iterate to calculate in the same manner, and when initial calculation, current reliability B was the third that step 3.2) obtains Intermediate reliability, later using reliability after the update of current common cause failure structure alignment processing acquisition as next common cause failure knot Current reliability B when structure alignment processing, to obtain the accurate reliability of remaining system;
3.4) it for the preceding g system selected in step 2.4), is arranged with accurate reliability according to descending, and therefrom select again For s system as the element in " selection subset S ", s is the number of " selection subset S " that step 2.4) obtains before taking.
2. a kind of multimode electric system redundancy optimization method based on sequence optimization and markov chain according to claim 1, It is characterized by: the first step specifically:
1.1) all multimode electric system are classified, to be parallel with add ons in multimode electric system and be connected in series Major component sum as classification foundation, by the identical multimode of sum for the major component for being parallel with add ons and being connected in series Electric system is classified as a classification;
1.2) be classified as same classification from fundamental system is successively extracted, fundamental system refers to that component population is minimum in same classification System, element includes add ons and major component in same classification;
1.3) reliability for finding out fundamental system in each classification is calculated with markov chain and its general generating function.
3. a kind of multimode electric system redundancy optimization method based on sequence optimization and markov chain according to claim 1, It is characterized by: the 4th step specifically:
For all multimode electric system that third step is screened, handled using the redundancy optimization algorithm of following formula Obtain optimal system:
MinC=∑ ci*ni
Wherein, c is the cost of each element in single multimode electric system, and n is the total quantity of element, and C is single multimode The cost of electric system.
CN201710944284.3A 2017-10-12 2017-10-12 Multimode electric system redundancy optimization method based on sequence optimization and markov chain Active CN107834538B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710944284.3A CN107834538B (en) 2017-10-12 2017-10-12 Multimode electric system redundancy optimization method based on sequence optimization and markov chain

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710944284.3A CN107834538B (en) 2017-10-12 2017-10-12 Multimode electric system redundancy optimization method based on sequence optimization and markov chain

Publications (2)

Publication Number Publication Date
CN107834538A CN107834538A (en) 2018-03-23
CN107834538B true CN107834538B (en) 2019-04-05

Family

ID=61647897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710944284.3A Active CN107834538B (en) 2017-10-12 2017-10-12 Multimode electric system redundancy optimization method based on sequence optimization and markov chain

Country Status (1)

Country Link
CN (1) CN107834538B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110601187B (en) * 2019-09-18 2021-02-05 浙江大学 Multi-state power system optimization construction method based on continuous discrete function

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9460478B2 (en) * 2012-12-17 2016-10-04 Arizona Board Of Regents On Behalf Of Arizona State University System and method for wind generation forecasting
CN105205502B (en) * 2015-10-30 2019-01-01 山东大学 A kind of Load time series classification method based on markov Monte Carlo
CN105762826B (en) * 2016-04-15 2018-12-18 华北电力大学 A kind of total transfer capability calculation method of ac and dc systems containing VSC-HVDC

Also Published As

Publication number Publication date
CN107834538A (en) 2018-03-23

Similar Documents

Publication Publication Date Title
CN104037776B (en) The electric network reactive-load capacity collocation method of random inertial factor particle swarm optimization algorithm
CN111416359B (en) Distribution network reconstruction method considering weighted power flow entropy
CN107871216A (en) A kind of recognition methods of power distribution network fragility node
CN109636009B (en) Method and system for establishing neural network model for determining line loss of power grid
CN108564592A (en) Based on a variety of image partition methods for being clustered to differential evolution algorithm of dynamic
CN107909185B (en) Multimode Power System Reliability Analysis method based on sequence optimization and markov chain
CN105868863B (en) Electromechanical system multi-level reliability prediction method based on fuzzy judgment
CN114692505A (en) Error prediction method and device for capacitor voltage transformer
CN107909194B (en) System-level testability design multi-objective optimization method
CN107834538B (en) Multimode electric system redundancy optimization method based on sequence optimization and markov chain
CN108446214A (en) Test case evolution generation method based on DBN
CN110456217A (en) A kind of MMC Fault Locating Method based on WPD-FOA-LSSVM dual model
CN111652478B (en) Umbrella algorithm-based power system voltage stability evaluation misclassification constraint method
CN107944097A (en) Plane enclosed cuts with scissors the numbering matrix description method of kinematic chain again
CN107038306A (en) A kind of optimal Latin hypercube experimental design method based on self-adapted genetic algorithm
CN104133766A (en) Software system maintainability evaluation and promotion method based on multi-objective community detection
CN107863771B (en) Multimode electric system redundancy optimization method based on sequence optimization and Monte Carlo
CN107909184B (en) Multimode Power System Reliability Analysis method based on sequence optimization and Monte Carlo
CN116306769A (en) Bayesian optimization method and system containing countermeasure generation network
CN105896545A (en) Accident chain mode based distributed power flow controller location method
Kumar et al. Comparative analysis of particle swarm optimization variants on distributed generation allocation for network loss minimization
CN117175548B (en) Sequential fault emergency control method based on random power flow
Samaan et al. Assessment of the annual frequency and duration indices in composite system reliability using genetic algorithms
CN110581554A (en) Power grid N-k fault analysis and screening method and device based on influence increment
CN110601187B (en) Multi-state power system optimization construction method based on continuous discrete function

Legal Events

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