CN110796377B - Power grid service system monitoring method supporting fuzzy theory - Google Patents
Power grid service system monitoring method supporting fuzzy theory Download PDFInfo
- Publication number
- CN110796377B CN110796377B CN201911049139.4A CN201911049139A CN110796377B CN 110796377 B CN110796377 B CN 110796377B CN 201911049139 A CN201911049139 A CN 201911049139A CN 110796377 B CN110796377 B CN 110796377B
- Authority
- CN
- China
- Prior art keywords
- subset
- time
- ute
- met
- fuzzy
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012544 monitoring process Methods 0.000 title claims abstract description 15
- 238000009826 distribution Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000005315 distribution function Methods 0.000 claims description 9
- 230000001174 ascending effect Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000007669 thermal treatment Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 description 5
- 230000010485 coping Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 230000002411 adverse Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000008186 active pharmaceutical agent Substances 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 238000011418 maintenance treatment Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a power grid service system monitoring method supporting a fuzzy theory, which comprises the steps of constructing basic probability data of a corresponding subset of events by considering the urgency degree of the events, and establishing another group of basic probability data by considering the dependency relationship and the deadline time among the events; and combining basic probability data obtained by multiple attributes, fusing the multiple attributes by adopting D-S to enable the multiple attributes to become the characteristics of the multiple attributes, deducing the membership degree and the non-membership degree of the fuzzy ending moment through the relation between the D-S and the IFS, further reversely deducing the membership degree and the non-membership degree of the fuzzy starting moment according to the operating time reasoning characteristic, wherein the starting time represents the moment most relevant to the scheduling sequence, so that scoring is further performed through the obtained membership degree and the non-membership degree of the fuzzy starting moment, and the events are sequentially loaded into the time scheduling sequence according to the fuzzy starting moment time sequence and the score, so that the scheduling of the parallel risk events for responding to the time sequence is realized.
Description
Technical Field
The invention relates to a power grid service system monitoring method supporting a fuzzy theory, and belongs to the field of power grid system monitoring.
Background
The power grid is a main constituent part of a power system and is one of important infrastructures for urban modern construction, and the advantages and disadvantages of the modern power grid not only directly influence the safe operation of a power department, but also relate to the development of national economy and the life quality of people. The comprehensive evaluation criteria for grid services are not absolute, but it is important to reduce the losses in the grid supply and to restore operation as soon as possible. The method and the device can accurately determine the coping time sequence of the parallel risk events occurring when the power grid system is abnormally operated, so that the power grid system can be recovered to normally operate as soon as possible, the energy consumption in the work of detecting the power grid system is reduced, and the adverse effects and the loss caused by the abnormal operation of the power grid system are reduced.
Disclosure of Invention
The invention provides a power grid service system monitoring method supporting a fuzzy theory, which is used for coping when a power grid system is abnormally operated and scientifically arranging scheduling of coping time sequences of parallel risk events.
The technical scheme of the invention is as follows: a power grid service system monitoring method supporting a fuzzy theory comprises the following steps:
step1, a plurality of events which have influences on the normal rotation of the power grid system and occur at the same time when the power grid service system is monitored are called as time uncertain events; obtaining any subset ET of the general urgent time uncertainty event UTE by adopting an average distribution methodiMET of1(ETi) (ii) a Obtaining any subset ET of urgent time uncertainty event UTE and untightened time uncertainty event UTE by adopting normal distribution methodiMET of1(ETi);
Step2, according to ETiAnd (3) respectively calculating the basic probability distribution number of each subset according to different numbers of the subset elements:
if the number of elements in the subset is 1, the number of elementary probability distributions for the subset is represented as: pt1,…,Ptj,…,Ptsigm(ii) a All subsets with the number of elements in the subset 1 in the time uncertainty event UTE are assigned a constraint probability PUTE1(ii) a Satisfies Pt1+…+Ptj+…+Ptsigm=PUTE1According to Pdep、PdlPush out PUTE1:
According to the obtained PUTE1Finding Ptj:
If the number of elements in the subset is m, the subset with the number of elements being m has C (sigm, m), and the basic probability distribution number for constructing the subset with the number of elements being m is expressed as: p(m,1),…,P(m,n),…,P(m,C(sigm,m))(ii) a All subsets with m elements in the subset of time uncertainty event UTE are assigned a constraint probability PUTEmAccording to PUTE1Calculating PUTEm:
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n):
The obtained Pt is subjected to heat treatmentj、P(m,n)MET as corresponding respective subsets2(ETi);
Wherein the fuzzy end time t of the time uncertainty event UTE is determinedjAs sample space DE, ETiRepresenting the ith subset of the sample space DE, and sigm representing the total number of elements of the sample space DE of the time uncertainty event UTE; t is tjRepresents the jth blurring end time and also represents the jth element of the sample space DE; pdep、PdlIs a constant number 0<Pdep、Pdl<1;j∈[1,sigm]C (sigm, m) represents a permutation combination;
step3, MET obtained from Step11(ETi) MET with results obtained at Step22(ETi) Obtaining MET (E) according to the basic probability orthogonality of the D-S evidence theoryTi);
Step4, according to MET (ET)j) To ETjDegree of membership mu (ET)j) Non-membership gamma (ET)j) The calculation of (2):
represents MET (ET)j'') of a plurality of different combinations, wherein ETj'' indicates all subsets in the event of uncertainty and the currently sought subset ETjSubsets whose intersections are not null;
wherein, ETjRefers to a subset with 1 element number;
step5, calculating the membership degree mu 1 (ET) of the matching fuzzy starting time according to the result obtained in the Step4k' ' ') and non-membership degree gamma 1(ETk''' ):
Wherein the fuzzy start time represents a difference between the fuzzy end time and the event execution time; ETk'' denotes ETjA subset formed by fuzzy start time corresponding to the middle element;
the hesitation degree of the blurring start timing is:
π(ETk''' )=1-(μ1(ETk''' )+γ1(ETk''' ))
mixing mu 1 (ET)k''' )、γ1(ETk''' )、π(ETk' ' ') into the scoring function:
step6, loading the events into the time scheduling sequence in turn according to the fuzzy start time sequence and the score.
For the general urgent time uncertainty event UTE, any subset ET of the general urgent time uncertainty event UTE is obtained by adopting an average distribution methodiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein, ω is1、ω2Representing the weight, ω1+ω2=1;ETjRefers to a subset of 1 element, ETpRefers to a subset of a plurality of elements.
The method for normally distributing the urgent time uncertainty event UTE and the untightened time uncertainty event UTE is adopted to obtain any subset ET of the urgent time uncertainty event UTE and the untightened time uncertainty event UTEiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
subjecting the obtained MET to a thermal treatment1(tj) T as a corresponding subset elementjSubset ET ofjMET of2(ETj);
For a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein σ represents a standard deviation of a normal distribution function and satisfies 2 σ sigm, ω3、ω4Representing the weight, ω3+ω 41 is ═ 1; mu get t1Means for expressing normal distribution function under non-urgent condition, mu is tsigmExpressing the mean of normal distribution functions in a pressing state; ETpRefers to a subset of a plurality of elements.
In Step2, for the case that the number of the subset elements is C (signm, m), the ordering of the subsets is 1,2,. n,. C (signm, m); the sequencing mode is as follows:
the subsets are sorted in an ascending order according to the size of the first element of the subsets, and then in an ascending order on the basis of the last sorting according to the sizes of the second element to the C (sigm, m).
The above-mentionedpriIndicating that the current time uncertainty event UTE depends on the number of other time uncertainty events UTE.
The above-mentioneddl represents the difference between the deadline of the current time uncertainty event UTE and the minimum deadline among all time uncertainty events.
The Step6 is specifically as follows:
if only a single event which is not loaded with the time scheduling sequence exists at the same time, loading the event into the time scheduling sequence;
if there are multiple events not loaded in the time schedule sequence at the same time, the high scoring events are loaded in the time schedule sequence.
The invention has the beneficial effects that: step1 of the invention considers the urgency degree of the events to construct basic probability data of the corresponding subset of the events, and step2 considers the dependency relationship and the deadline time among the events to establish another group of basic probability data; combining basic probability data obtained by multiple attributes, fusing the multiple attributes by using D-S to enable the multiple attributes to become the characteristics of the multiple attributes, deducing the membership and the non-membership of a fuzzy ending moment through the relation between the D-S and the IFS, further reversely deducing the membership and the non-membership of the fuzzy starting moment according to the reasoning characteristic of working time, wherein the starting time represents the moment most relevant to a scheduling sequence, so that the obtained membership and the non-membership of the fuzzy starting moment are used for scoring, and events are sequentially loaded into the time scheduling sequence according to the time sequence and the score of the fuzzy starting moment, so that the scheduling of parallel risk event handling time sequence is realized, the working efficiency of a power grid, the running stability of the power system are enhanced, the energy loss and the manpower resource waste in detection are reduced, and the speed of abnormal work restoration of the power grid system is accelerated through the scheduling, the adverse effects of business service stop and the like and adverse consequences of incapability of timely recovering normal operation due to unscientific coping of power grid system risk events are reduced.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a time uncertainty event constraint topology diagram;
FIG. 3 is a time interval diagram of a time uncertainty event;
fig. 4 is a time indeterminate event scheduling sequence diagram.
Detailed Description
Example 1: as shown in fig. 1 to 4, a method for monitoring a power grid service system supporting a fuzzy theory includes the following steps:
step1, normally rotating a plurality of pairs of power grid systems at the same time when the power grid service systems are monitoredThe influential events are called time indeterminate events; obtaining any subset ET of the general urgent time uncertainty event UTE by adopting an average distribution methodiMET of1(ETi) (ii) a Obtaining any subset ET of urgent time uncertainty event UTE and untightened time uncertainty event UTE by adopting normal distribution methodiMET of1(ETi);
Step2, according to ETiAnd (3) respectively calculating the basic probability distribution number of each subset according to different numbers of the subset elements:
if the number of elements in the subset is 1, the number of elementary probability distributions for the subset is represented as: pt1,…,Ptj,…,Ptsigm(ii) a All subsets with the number of elements in the subset 1 in the time uncertainty event UTE are assigned a constraint probability PUTE1(ii) a Satisfies Pt1+…+Ptj+…+Ptsigm=PUTE1According to Pdep、PdlPush out PUTE1:
According to the obtained PUTE1Finding Ptj:
If the number of elements in the subset is m, the subset with the number of elements being m has C (sigm, m), and the basic probability distribution number for constructing the subset with the number of elements being m is expressed as: p(m,1),…,P(m,n),…,P(m,C(sigm,m))(ii) a All subsets with m elements in the subset of time uncertainty event UTE are assigned a constraint probability PUTEmAccording to PUTE1Calculating PUTEm:
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n):
The obtained Pt is subjected to heat treatmentj、P(m,n)MET as corresponding respective subsets2(ETi);
Wherein the fuzzy end time t of the time uncertainty event UTE is determinedjAs sample space DE, ETiRepresenting the ith subset of the sample space DE, and sigm representing the total number of elements of the sample space DE of the time uncertainty event UTE; t is tjRepresents the jth blurring end time and also represents the jth element of the sample space DE; pdep、PdlIs a constant number 0<Pdep、Pdl<1;j∈[1,sigm]C (sigm, m) represents a permutation combination;
step3, MET obtained from Step11(ETi) MET with results obtained at Step22(ETi) Obtaining MET (ET) according to the basic probability orthogonality of the D-S evidence theoryi);
Step4, according to MET (ET)j) To ETjDegree of membership mu (ET)j) Non-membership gamma (ET)j) The calculation of (2):
represents MET (ET)j'') of a plurality of different combinations, wherein ETj'' indicates all subsets in the event of uncertainty and the currently sought subset ETjSubsets whose intersections are not null;
wherein, ETjRefers to a subset with 1 element number;
step5, calculating the membership degree mu 1 (ET) of the matching fuzzy starting time according to the result obtained in the Step4k' ' ') and non-membership degree gamma 1(ETk''' ):
Wherein the fuzzy start time represents a difference between the fuzzy end time and the event execution time; ETk'' denotes ETjA subset formed by fuzzy start time corresponding to the middle element;
the hesitation degree of the blurring start timing is:
π(ETk''' )=1-(μ1(ETk''' )+γ1(ETk''' ))
mixing mu 1 (ET)k''' )、γ1(ETk''' )、π(ETk' ' ') into the scoring function:
step6, loading the events into the time scheduling sequence in turn according to the fuzzy start time sequence and the score.
Further, the general urgent time uncertainty event UTE can be set, and any subset ET thereof can be obtained by adopting an average distribution methodiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein, ω is1、ω2Representing the weight, ω1+ω2=1;ETjRefers to a subset of 1 element, ETpRefers to a subset of a plurality of elements.
Further, it may be set that the normal distribution method is adopted for both the urgent time uncertainty event UTE and the urgent time uncertainty event UTE to obtain any subset ET thereofiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
subjecting the obtained MET to a thermal treatment1(tj) T as a corresponding subset elementjSubset ET ofjMET of2(ETj);
For a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein σ represents a standard deviation of a normal distribution function and satisfies 2 σ sigm, ω3、ω4Representing the weight, ω3+ω 41 is ═ 1; mu get t1Means for expressing normal distribution function under non-urgent condition, mu is tsigmTo representThe mean of normal distribution functions in a pressing state; ETpRefers to a subset of a plurality of elements.
Further, in Step2, for the case that the number of the subset elements is C (sigm, m), the ordering of the subsets may be 1,2,. n.,. C (sigm, m); the sequencing mode is as follows:
the subsets are sorted in an ascending order according to the size of the first element of the subsets, and then in an ascending order on the basis of the last sorting according to the sizes of the second element to the C (sigm, m).
The above-mentionedpriIndicating that the current time uncertainty event UTE depends on the number of other time uncertainty events UTE. (As shown in FIG. 2, the number is the number of current time uncertainty events counted by the arrow)
Further, the above may be provideddl represents the difference between the deadline of the current time uncertainty event UTE and the minimum deadline among all time uncertainty events. The cutoff time represents the last fuzzy end time in the sample space of the current time uncertainty event UTE;
further, the Step6 may be set to specifically be:
if only a single event which is not loaded with the time scheduling sequence exists at the same time, loading the event into the time scheduling sequence;
if there are multiple events not loaded in the time schedule sequence at the same time, the high scoring events are loaded in the time schedule sequence.
Further, the following example is made for the steps in the present application:
supposing that when the abnormal condition of the power grid system occurs, a plurality of risk events are monitored to occur simultaneously, and the risk events are called as time uncertainty events UTE. How to get a response to multiple risk events when these risk events occur simultaneouslyIt is crucial that the sequence is such that an optimal healing effect is achieved. The monitored time uncertainty event UTE in the power grid system comprises the following steps: UTE1、UTE2、UTE3、UTE4、UTE5(e.g., the events may be, in order, a system response timeout event, a system access slow event, a system unavailable event, a network failure event, a service down event). A time uncertainty event constraint (UTEC) topological graph is shown in fig. 2 and a time uncertainty event time interval is shown in fig. 3.
The steps in this example are all UTE5For example, the following steps are carried out:
step1, the time uncertainty event UTE adopts different methods according to the respective urgency: UTE2The degree of influence on the operation of the power grid does not change with the lapse of time, so that the average distribution is adopted; UTE1、UTE3、UTE4、 UTE5The influence on the normal operation of the power grid is more and more serious along with the time, so normal distribution is adopted. And the UTE can be better distinguished according to different urgency degrees by adopting different modes for classification, and the UTE is more suitable for actual conditions. The preset information of each event of the UTEF is shown in Table 1 (the data given in the embodiment are all generated at the integral point time, but when the system is actually used for collecting event information, half an hour can be taken as a step length, such as 0.5h and 1h, and for the actual time between 0h and 0.5h, 0.5h is selected as the fuzzy end time, other similar reasons, further, the fuzzy end time interval is designed to be the same event in all the participating events, and the maximum number of the events is 2):
TABLE 1 Preset information for time uncertain event streams
Known UTE5End of blur time tjHas a sample space DE of {3,4}, and ETiAs any subset of DE, i.e. ET1{3} or ET2{4} or ET3The basic probability number is not considered because the empty set can not be calculated, and the following processes are not considered because of the same reasonAnd (4) empty collection. UTE5The blur end time of (2) is time 3 and time 4, and is obtained by an average distribution according to the degree of urgency: MET1{3}=0.3333、MET1{4}=0.3333、 MET1{3,4}=0.3334。
Step2, let UTE5The fuzzy end interval of (2) has 2 primitive elements: t is t1,t2And t is1=3,t2Further, let Pt be 41+Pt2=PUTE1The resulting probabilities are shown in table 2 below:
TABLE 2
When in UTE5Get t1When the value is 3, P is addeddep、PdlSubstituting the value of (b) into the formula:
and according to the obtained PUTE1Determining Pt as 0.65001The value of (c):
by analogy, the basic probability distribution number Pt of all elementary elements can be obtainedjWhen t is2When equal to 4, P40.2167, i.e. MET2(ET1)=MET2{3}=0.4333,MET2(ET2)=MET2{4} ═ 0.2167. After the formula is improved, the event can be scheduled as early as possible in the scheduling process, and the event scheduling efficiency is improved.
For UTE5When the number of elements in the subset is m is 2, that is, {3,4}, there is 1 in the subset having the number of elements m, that is, C (sigm, m) ═ C (2, 2).
According to PUTE1Calculating PUTEm:
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n):
Step3, in UTE5When ET1When the measured value is {3}, MET obtained by Step1 is obtained according to the orthogonality of basic probabilities of D-S evidence theory1{3} MET obtained with Step22{3} orthogonal combination result MET {3 }.
When ET is a non-empty set, the following orthogonal formula is calculated:
in this example, substituting each datum into the formula yields:
the same can be obtained:
MET (ET) result orthogonally found from fundamental probability of D-S evidence theoryi) The characteristics of the two calculation modes are integrated, the accuracy of the obtained result is improved, the result is more comprehensive, and the multi-attribute characteristic is expressed by using the result of the single attribute.
In other UTEs, repeating the above process can obtain MET calculated in two ways from Step1-Step3 in all UTEs1(ETi)、MET2(ETi) And finally orthogonally combined MET (ET)i) As a result, the results are asTable 3 below shows:
TABLE 3
Step4, according to MET (ET)j) To ETjDegree of membership mu (ET)j) Non-membership gamma (ET)j) The calculation of (2):
in UTE5In when taking ET1When not equal to 3, MET (ET)1) MET {3} ═ 0.5177, substituting the following equation:
because of ET1Is only 3 regardless of the empty set, so ET1' {3}, the cumulative sum of the elementary probability numbers is MET {3} itself:
all subsets ET1={3}、ET2={4}、ET3In {3,4}, ET1={3}、ET3The intersection of {3,4} and {3} is not null, and the condition is satisfied. Since only one {4} of the total subset satisfies the condition, the cumulative sum of the elementary probability numbers is MET {4 }:
by analogy, since the results from the D-S evidence theory do not have the advantage of time derivation, Step4 prepares for the next application of the negative-going time inference theory of IFS, the results are shown in table 4 below:
TABLE 4
Step5, in UTE5In when getting tjWhen the time is 3, the corresponding module is obtained according to the negative time reasoning theory of IFSPaste start time tkThe intuitive fuzzy set membership of the fuzzy set 2 on the sample space DS is:
in the same way, find tjWhen 4, the corresponding blurring start time tkThe intuitive fuzzy set membership of the fuzzy set of 3 on the sample space DS is:
solving UTE through an intuitionistic fuzzy set negative time reasoning theory5When t is injT corresponding to 3iDegree of membership μ 1 (ET) at the start of blurring of 2k' ' ') and non-membership degree gamma 1(ETk' ' ') and further determining a blurring start time tiHesitation degree pi (ET) of 2k''' ):
π(ETk”')=1-(μ1(ETk''' )+γ1(ETk''' ))=1-(0.7589+0.1667)=0.0744
And finally, utilizing a scoring function formula:
by analogy, all the membership degree, the non-membership degree and the hesitation degree of the fuzzy starting time and the scores of all the fuzzy starting times are obtained, the characteristics and the calculation process of fusing the various attributes in the front can be embodied by using the intuitive digital scores, the scheduling sequence of the UTE is clearly embodied, and the obtained result is shown in the following table 5:
TABLE 5
Score each UTE high in time series and verticallyThe low-level, time scheduling queues are loaded in sequence, and the finally obtained scheduling sequence is shown in fig. 4, and a brief analysis shows that: UTE at time 25Loading a time scheduling queue; at time 3 due to UTE5Has already been performed, UTE4Loading a time scheduling queue; UTE is compared at time 43And UTE2Fuzzy start time score, UTE since 0.8784 > 0.68143Loading a time scheduling queue; UTE is compared at time 52And UTE1Fuzzy start time score, UTE since 0.6459 < 0.66501Loading a time scheduling queue; at time 6, due to UTE1、UTE3、UTE4、UTE5Has all been loaded into the time scheduling queue, so UTE will be2The time scheduling queue is loaded.
In summary, the final response sequence is: UTE5, UTE4, UTE3, UTE1, UTE 2. When abnormal power supply condition occurs in a certain area, the system carries out maintenance treatment on the conditions of service downtime, network fault, system unavailability, system response timeout, system access slow and the like in sequence.
The working principle of the invention is as follows: according to the method, when the power grid system is abnormal in operation, the risk event response time sequence generated in the power grid system is scheduled according to the fuzzy theory. First, MET is obtained by two different ways1(ETi)、 MET2(ETi) For fusing multi-attribute features, MET is mapped according to D-S evidence theory1(ETi)、MET2(ETi) Orthogonal combination is carried out to obtain MET (ET) fusing the information obtained in the calculation process of the two waysi) Substituting the result into a formula to calculate the IFS membership and the non-membership of each UTE fuzzy end moment, providing for utilizing the advantages in the aspect of time derivation of an intuitive fuzzy set, deducing the corresponding membership and the non-membership of the fuzzy start moment and the hesitation degree of the fuzzy start moment from the IFS membership and the non-membership of each UTE fuzzy end moment and the hesitation degree of the fuzzy start moment through an intuitive fuzzy set negative time reasoning theory, and finally obtaining the probability score of each UTE fuzzy start moment according to a scoring function formula, thereby sequentially loading the UTEs into a time scheduling queue and establishing an accurate and effective risk event coping partAnd (4) processing corresponding time sequence.
While the present invention has been described in detail with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, and various changes can be made without departing from the spirit of the present invention within the knowledge of those skilled in the art.
Claims (7)
1. A power grid service system monitoring method supporting a fuzzy theory is characterized in that: the method comprises the following steps:
step1, a plurality of events which have influences on the normal rotation of the power grid system and occur at the same time when the power grid service system is monitored are called as time uncertain events; obtaining any subset ET of the general urgent time uncertainty event UTE by adopting an average distribution methodiMET of1(ETi) (ii) a Obtaining any subset ET of urgent time uncertainty event UTE and untightened time uncertainty event UTE by adopting normal distribution methodiMET of1(ETi);
Step2, according to ETiAnd (3) respectively calculating the basic probability distribution number of each subset according to different numbers of the subset elements:
if the number of elements in the subset is 1, the number of elementary probability distributions for the subset is represented as: pt1,…,Ptj,…,Ptsigm(ii) a All subsets with the number of elements in the subset 1 in the time uncertainty event UTE are assigned a constraint probability PUTE1(ii) a Satisfies Pt1+…+Ptj+…+Ptsigm=PUTE1According to Pdep、PdlPush out PUTE1:
According to the obtained PUTE1Finding Ptj:
If the number of elements in the subset is m, the subset with the number of elements being m has C (sigm, m), and the basic probability distribution number for constructing the subset with the number of elements being m is expressed as: p(m,1),…,P(m,n),…,P(m,C(sigm,m))(ii) a All subsets with m elements in the subset of time uncertainty event UTE are assigned a constraint probability PUTEmAccording to PUTE1Calculating PUTEm:
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n):
The obtained Pt is subjected to heat treatmentj、P(m,n)MET as corresponding respective subsets2(ETi);
Wherein the fuzzy end time t of the time uncertainty event UTE is determinedjAs sample space DE, ETiRepresenting the ith subset of the sample space DE, and sigm representing the total number of elements of the sample space DE of the time uncertainty event UTE; t is tjRepresents the jth blurring end time and also represents the jth element of the sample space DE; pdep、PdlIs a constant number 0<Pdep、Pdl<1;j∈[1,sigm]C (sigm, m) represents a permutation combination;
step3, MET obtained from Step11(ETi) MET with results obtained at Step22(ETi) Obtaining MET (ET) according to the basic probability orthogonality of the D-S evidence theoryi);
Step4, according to MET (ET)j) To ETjDegree of membership mu (ET)j) Non-membership gamma (ET)j) Is/are as followsAnd (3) calculating:
represents MET (ET)j'') of a plurality of different combinations, wherein ETj'' indicates all subsets in the event of uncertainty and the currently sought subset ETjSubsets whose intersections are not null;
wherein, ETjRefers to a subset with 1 element number;
step5, calculating the membership degree mu 1 (ET) of the matching fuzzy starting time according to the result obtained in the Step4k' ' ') and non-membership degree gamma 1(ETk''' ):
Wherein the fuzzy start time represents a difference between the fuzzy end time and the event execution time; ETk'' denotes ETjA subset formed by fuzzy start time corresponding to the middle element;
the hesitation degree of the blurring start timing is:
π(ETk''' )=1-(μ1(ETk''' )+γ1(ETk''' ))
mixing mu 1 (ET)k- ''' )、γ1(ETk''' )、π(ETk' ' ') into the scoring function:
step6, loading the events into the time scheduling sequence in turn according to the fuzzy start time sequence and the score.
2. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: for the general urgent time uncertainty event UTE, any subset ET of the general urgent time uncertainty event UTE is obtained by adopting an average distribution methodiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein, ω is1、ω2Representing the weight, ω1+ω2=1;ETjRefers to a subset of 1 element, ETpRefers to a subset of a plurality of elements.
3. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: the method for normally distributing the urgent time uncertainty event UTE and the untightened time uncertainty event UTE is adopted to obtain any subset ET of the urgent time uncertainty event UTE and the untightened time uncertainty event UTEiMET of1(ETi):
For a subset of one element number, the basic probability number calculation formula is as follows:
subjecting the obtained MET to a thermal treatment1(tj) T as a corresponding subset elementjSubset ET ofjMET of2(ETj);
For a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
wherein σ represents a standard deviation of a normal distribution function and satisfies 2 σ sigm, ω3、ω4Representing the weight, ω3+ω41 is ═ 1; mu get t1Means for expressing normal distribution function under non-urgent condition, mu is tsigmExpressing the mean of normal distribution functions in a pressing state; ETpRefers to a subset of a plurality of elements.
4. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: in Step2, for the case that the number of the subset elements is C (signm, m), the ordering of the subsets is 1,2,. n,. C (signm, m); the sequencing mode is as follows:
the subsets are sorted in an ascending order according to the size of the first element of the subsets, and then in an ascending order on the basis of the last sorting according to the sizes of the second element to the C (sigm, m).
7. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: the Step6 is specifically as follows:
if only a single event which is not loaded with the time scheduling sequence exists at the same time, loading the event into the time scheduling sequence;
if there are multiple events not loaded in the time schedule sequence at the same time, the high scoring events are loaded in the time schedule sequence.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911049139.4A CN110796377B (en) | 2019-10-31 | 2019-10-31 | Power grid service system monitoring method supporting fuzzy theory |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911049139.4A CN110796377B (en) | 2019-10-31 | 2019-10-31 | Power grid service system monitoring method supporting fuzzy theory |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110796377A CN110796377A (en) | 2020-02-14 |
CN110796377B true CN110796377B (en) | 2022-03-29 |
Family
ID=69442195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911049139.4A Active CN110796377B (en) | 2019-10-31 | 2019-10-31 | Power grid service system monitoring method supporting fuzzy theory |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110796377B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116976682B (en) * | 2023-09-22 | 2023-12-26 | 安徽融兆智能有限公司 | Fuzzy algorithm-based operation state evaluation method for electricity consumption information acquisition system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624782A (en) * | 2011-10-31 | 2012-08-01 | 李宗诚 | Internal concentrated harmonization system/information and communication technology (ICH/ICT) information fusion basis of internet |
CN102802158A (en) * | 2012-08-07 | 2012-11-28 | 湖南大学 | Method for detecting network anomaly of wireless sensor based on trust evaluation |
CN105809287A (en) * | 2016-03-10 | 2016-07-27 | 云南大学 | High-voltage transmission line icing process integrated prediction method |
CN106932670A (en) * | 2017-02-08 | 2017-07-07 | 国家电网公司 | A kind of distribution power automation terminal method for diagnosing status based on D S evidence theories |
CN108877218A (en) * | 2018-07-04 | 2018-11-23 | 西北工业大学 | Method for controlling traffic signal lights based on D-S evidence theory |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7742972B2 (en) * | 1999-07-21 | 2010-06-22 | Longitude Llc | Enhanced parimutuel wagering |
US9680855B2 (en) * | 2014-06-30 | 2017-06-13 | Neo Prime, LLC | Probabilistic model for cyber risk forecasting |
-
2019
- 2019-10-31 CN CN201911049139.4A patent/CN110796377B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624782A (en) * | 2011-10-31 | 2012-08-01 | 李宗诚 | Internal concentrated harmonization system/information and communication technology (ICH/ICT) information fusion basis of internet |
CN102802158A (en) * | 2012-08-07 | 2012-11-28 | 湖南大学 | Method for detecting network anomaly of wireless sensor based on trust evaluation |
CN105809287A (en) * | 2016-03-10 | 2016-07-27 | 云南大学 | High-voltage transmission line icing process integrated prediction method |
CN106932670A (en) * | 2017-02-08 | 2017-07-07 | 国家电网公司 | A kind of distribution power automation terminal method for diagnosing status based on D S evidence theories |
CN108877218A (en) * | 2018-07-04 | 2018-11-23 | 西北工业大学 | Method for controlling traffic signal lights based on D-S evidence theory |
Non-Patent Citations (1)
Title |
---|
"基于云-模糊融合技术的电机故障诊断";王惠中;《计算机测量与控制》;20150825(第8期);第2614-2616、2621页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110796377A (en) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111475804B (en) | Alarm prediction method and system | |
US7472097B1 (en) | Employee selection via multiple neural networks | |
WO2019218475A1 (en) | Method and device for identifying abnormally-behaving subject, terminal device, and medium | |
CN112016905B (en) | Information display method and device based on approval process, electronic equipment and medium | |
CN112650580B (en) | Industrial big data monitoring system based on edge calculation | |
CN110555477A (en) | municipal facility fault prediction method and device | |
CN116307928A (en) | Object supervision management system | |
CN112416562B (en) | Method and device for distributed task scheduling engine | |
CN110796377B (en) | Power grid service system monitoring method supporting fuzzy theory | |
CN109412838A (en) | Server cluster host node selection method based on hash calculating and Performance Evaluation | |
CN113190426B (en) | Stability monitoring method for big data scoring system | |
CN112906738A (en) | Water quality detection and treatment method | |
CN114647684A (en) | Traffic prediction method and device based on stacking algorithm and related equipment | |
US6978222B2 (en) | Method of determining level of effect of system entity on system performance, by use of active time of same entity | |
CN114401158A (en) | Flow charging method and device, electronic equipment and storage medium | |
Baykasoğlu et al. | Genetic programming based data mining approach to dispatching rule selection in a simulated job shop | |
CN105897498A (en) | Business monitoring method and device | |
CN109740750B (en) | Data collection method and device | |
CN114817408B (en) | Scheduling resource identification method and device, electronic equipment and storage medium | |
Chandio et al. | Towards adaptable and tunable cloud-based map-matching strategy for GPS trajectories | |
CN115080215A (en) | Method and system for performing task scheduling among computing nodes by state monitoring chip | |
CN115130913A (en) | New energy project risk assessment and quantification method and system | |
Garcia-Almanza et al. | Evolving decision rules to predict investment opportunities | |
CN113888318A (en) | Risk detection method and system | |
CN113673957A (en) | Offline crowdsourcing labeling method for text data |
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 |