CN110796377B - Power grid service system monitoring method supporting fuzzy theory - Google Patents

Power grid service system monitoring method supporting fuzzy theory Download PDF

Info

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
Application number
CN201911049139.4A
Other languages
Chinese (zh)
Other versions
CN110796377A (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.)
Information Center of Yunnan Power Grid Co Ltd
Original Assignee
Information Center of Yunnan Power Grid Co Ltd
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 Information Center of Yunnan Power Grid Co Ltd filed Critical Information Center of Yunnan Power Grid Co Ltd
Priority to CN201911049139.4A priority Critical patent/CN110796377B/en
Publication of CN110796377A publication Critical patent/CN110796377A/en
Application granted granted Critical
Publication of CN110796377B publication Critical patent/CN110796377B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • 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
    • 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

  • 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

Power grid service system monitoring method supporting fuzzy theory
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
Figure BDA0002254849920000011
According to the obtained PUTE1Finding Ptj
Figure BDA0002254849920000021
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
Figure BDA0002254849920000022
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n)
Figure BDA0002254849920000023
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):
Figure BDA0002254849920000024
Figure BDA0002254849920000025
means all belonging to ETjSubset ET ofjMET (ET) of `j') a cumulative sum;
Figure BDA0002254849920000026
Figure BDA0002254849920000027
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''' ):
Figure BDA0002254849920000031
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:
Figure BDA0002254849920000032
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:
Figure BDA0002254849920000033
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
Figure BDA0002254849920000034
wherein, ω is1、ω2Representing the weight, ω12=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:
Figure BDA0002254849920000041
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:
Figure BDA0002254849920000042
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-mentioned
Figure BDA0002254849920000043
priIndicating that the current time uncertainty event UTE depends on the number of other time uncertainty events UTE.
The above-mentioned
Figure BDA0002254849920000044
dl 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
Figure BDA0002254849920000051
According to the obtained PUTE1Finding Ptj
Figure BDA0002254849920000052
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
Figure BDA0002254849920000061
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n)
Figure BDA0002254849920000062
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):
Figure BDA0002254849920000063
Figure BDA0002254849920000064
means all belonging to ETjSubset ET ofjMET (ET) of `j') a cumulative sum;
Figure BDA0002254849920000065
Figure BDA0002254849920000066
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''' ):
Figure BDA0002254849920000071
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:
Figure BDA0002254849920000072
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:
Figure BDA0002254849920000073
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
Figure BDA0002254849920000074
wherein, ω is1、ω2Representing the weight, ω12=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:
Figure BDA0002254849920000081
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:
Figure BDA0002254849920000082
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-mentioned
Figure BDA0002254849920000083
priIndicating 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 provided
Figure BDA0002254849920000084
dl 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
Figure BDA0002254849920000091
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
Figure BDA0002254849920000101
When in UTE5Get t1When the value is 3, P is addeddep、PdlSubstituting the value of (b) into the formula:
Figure BDA0002254849920000102
and according to the obtained PUTE1Determining Pt as 0.65001The value of (c):
Figure BDA0002254849920000103
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
Figure BDA0002254849920000104
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n)
Figure BDA0002254849920000105
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:
Figure BDA0002254849920000106
in this example, substituting each datum into the formula yields:
Figure BDA0002254849920000111
the same can be obtained:
Figure BDA0002254849920000112
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
Figure BDA0002254849920000113
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:
Figure BDA0002254849920000121
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 }:
Figure BDA0002254849920000122
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
Figure BDA0002254849920000123
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:
Figure BDA0002254849920000124
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:
Figure BDA0002254849920000125
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:
Figure BDA0002254849920000131
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
Figure BDA0002254849920000132
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
Figure FDA0002254849910000011
According to the obtained PUTE1Finding Ptj
Figure FDA0002254849910000012
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
Figure FDA0002254849910000013
According to PUTEmCalculating the number of basic probability distribution P with the n-th subset element number m(m,n)
Figure FDA0002254849910000014
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:
Figure FDA0002254849910000021
Figure FDA0002254849910000022
means all belonging to ETjSubset ET ofjMET (ET) of `j') a cumulative sum;
Figure FDA0002254849910000023
Figure FDA0002254849910000024
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''' ):
Figure FDA0002254849910000025
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:
Figure FDA0002254849910000031
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:
Figure FDA0002254849910000032
for a subset with a plurality of element numbers, the basic probability number calculation formula is as follows:
Figure FDA0002254849910000033
wherein, ω is1、ω2Representing the weight, ω12=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:
Figure FDA0002254849910000034
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:
Figure FDA0002254849910000035
wherein σ represents a standard deviation of a normal distribution function and satisfies 2 σ sigm, ω3、ω4Representing the weight, ω341 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).
5. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: the above-mentioned
Figure FDA0002254849910000041
priIndicating that the current time uncertainty event UTE depends on the number of other time uncertainty events UTE.
6. The method for monitoring the power grid service system supporting the fuzzy theory according to claim 1, wherein: the above-mentioned
Figure FDA0002254849910000042
dl represents the difference between the deadline of the current time uncertainty event UTE and the minimum deadline among all time uncertainty events.
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.
CN201911049139.4A 2019-10-31 2019-10-31 Power grid service system monitoring method supporting fuzzy theory Active CN110796377B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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