CN112183995A - Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting - Google Patents

Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting Download PDF

Info

Publication number
CN112183995A
CN112183995A CN202011008701.1A CN202011008701A CN112183995A CN 112183995 A CN112183995 A CN 112183995A CN 202011008701 A CN202011008701 A CN 202011008701A CN 112183995 A CN112183995 A CN 112183995A
Authority
CN
China
Prior art keywords
interval
protection system
score
index
defect
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011008701.1A
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hubei Electric Power 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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202011008701.1A priority Critical patent/CN112183995A/en
Publication of CN112183995A publication Critical patent/CN112183995A/en
Pending legal-status Critical Current

Links

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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F1/00—Details not covered by groups G06F3/00 - G06F13/00 and G06F21/00
    • G06F1/26—Power supply means, e.g. regulation thereof
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • General Engineering & Computer Science (AREA)
  • Development Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Evolutionary Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提供一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，包括：根据直流保护系统可靠性综合评估指标体系组成，构建指标体系层次；依据直流保护系统的各个层次以及缺陷所在部位不同，对直流保护系统进行分类，统计以往各换流站的全部缺陷数据，根据缺陷数据设置直流保护系统各个层次的分值；根据直流保护系统缺陷的历史数据得出直流保护系统缺陷的指标数据；对指标量进行熵权法，得出三个指标的权重；运用熵权法对直流保护系统各类设备及整个系统进行可靠性评价；引入权重因子，将区间层次分析法和区间熵权法相结合，得到指标评价综合区间权重，由此对直流保护系统进行区间层次法和区间熵权法混合赋权法可靠性评价。The invention provides a reliability evaluation method of a DC protection system based on the interval analysis hierarchy process and the interval entropy weight mixed weighting, including: constructing the index system hierarchy according to the composition of the DC protection system reliability comprehensive evaluation index system; Each level and the location of the defect are different, classify the DC protection system, count all the defect data of each converter station in the past, and set the score of each level of the DC protection system according to the defect data; according to the historical data of the DC protection system defect, get the DC The index data of the defects of the protection system; the entropy weight method is applied to the index quantity to obtain the weights of the three indicators; the entropy weight method is used to evaluate the reliability of various equipment and the entire system of the DC protection system; the weight factor is introduced to analyze the interval level. Combined with the interval entropy weight method, the comprehensive interval weight of the index evaluation is obtained, and the reliability of the DC protection system is evaluated by the interval hierarchy method and the interval entropy weight method.

Description

Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting
Technical Field
The invention relates to the technical field of direct current protection system reliability research, in particular to a direct current protection system reliability evaluation method based on mixed weighting of an interval analytic hierarchy process and an interval entropy weight.
Background
The direct current protection is the core of the secondary part of the high-voltage direct current transmission system, is the first defense line for safe and reliable operation of the direct current transmission system, and the reliability directly determines the safe operation of direct current transmission equipment and the stable operation of a power grid, so that the reliability of the direct current protection system needs to be evaluated. For the reliability evaluation of the alternating current power grid, a large amount of research has been carried out at home and abroad, and currently, widely applied evaluation methods include a markov model method, a fault tree analysis method, a GO probability analysis method and the like. The reliability evaluation of the direct current protection system is less, and the evaluation method is not intuitive and has poor applicability.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a direct current protection system reliability evaluation method based on the mixed weighting of the interval analytic hierarchy process and the interval entropy weight.
The invention aims to be realized by the following technical scheme:
a direct current protection system reliability evaluation method based on mixed weighting of an inter-zone analytic method and an inter-zone entropy weight comprises the following steps:
the method comprises the following steps: according to the DC protection system reliability comprehensive evaluation index system composition, constructing an index system hierarchy: the method comprises the steps that a target layer, a criterion layer B and a scheme layer C are compared and scored pairwise, the index layers of the power user equipment are obtained, the interval weight of each layer is obtained through an inter-layer analytic method, a characteristic vector is obtained according to a judgment matrix of the layer B when the scheme layer C is aligned, a total ranking weight of the layers is established, and the total ranking weight of the scheme layer C is obtained;
step two: classifying the direct current protection system according to different levels and different positions of defects of the direct current protection system, counting all defect data of each converter station in the past, setting scores of each level of the direct current protection system according to the defect data, setting the scores as the bottom scores of each element of the direct current protection system, and obtaining the score of the interval hierarchy method of the direct current protection system by combining with the total sorting weight of the scheme layer C;
step three: classifying the direct current protection system according to different levels and different positions of the defects of the direct current protection system, collecting historical data of the defects of the direct current protection system, and obtaining three index data of the defects of the direct current protection system according to the historical data: the component defect-free time rate MTTF, the component mean repair time rate MTTR and the component defect rate EDR;
step four: performing an entropy weight method on the three indexes in the third step to obtain weights of the three indexes, and setting two weight values in the entropy weight method, wherein the numerical value is unchanged after the elements without the defect data in the elements are converted into a standard matrix, and the numerical value is set to be 0.5 after the elements without the defect data in the elements are converted into the standard matrix, so that the score of the entropy weight method of the whole direct current protection system is obtained;
step five: introducing a weight factor theta, combining an interval hierarchy analysis method and an interval entropy weight method to obtain an index evaluation comprehensive interval score w, wherein the comprehensive interval score w changes along with the change of theta, and the calculation formula of the index evaluation comprehensive interval score w is as follows:
w＝θwIAHP+(1-θ)wIEN,0≤θ≤1 (19)
wherein wIAHPGrading by interval hierarchy method obtained in step two, wIENAnd step four obtained in step four, when theta is 1 or 0, the interval hierarchy method and the interval entropy weight method are respectively corresponded.
Further, the specific implementation process of the step one is as follows:
(a) establishing each layer of the DC protection system according to the composition of the DC protection system, and establishing a judgment matrix of an analytic hierarchy process, wherein each layer of the DC protection system is shown in the following table 1, wherein the DC protection system (A)1) Devices, apparatus (B) for the target layerj) The interface board, CT, PT (C) are criterion layersi) Is a scheme layer;
TABLE 1 respective hierarchy of DC protection system
Figure RE-GDA0002783794870000021
Figure RE-GDA0002783794870000031
(b) According to each level of the direct current protection system, selecting experts to carry out pairwise comparison and scoring on index levels of the power user equipment, carrying out pairwise comparison and scoring on the power user indexes according to a 'interval mark table' in a table 2, wherein scoring results are generally carried out according to a in an intervalijPoint and width mu, wherein the middle point a of the intervalijThe value of (1) is a base number of judgment or an average value of a random variable in a judgment interval scale; width mu is a base a given by the expert according to uncertainty and ambiguity of the judgment after determining a good baseijThe interval variation range of (a);
Figure RE-GDA0002783794870000041
TABLE 2 interval ratio Scale Table
Figure RE-GDA0002783794870000042
Forming interval number according to the score of the interval proportion scale table
Figure RE-GDA0002783794870000043
When index i is more important than index j, i.e. aijWhen not less than 1, i is not equal to j
Figure RE-GDA0002783794870000044
When the index j is more important than the index i, i.e. ajiWhen j is not equal to i and is not equal to 1
Figure RE-GDA0002783794870000045
(c) Constructing interval judgment matrix A
The selected experts compare and score each index of the same layer with respect to the importance of the index of the previous layer according to an interval proportional scaling method to construct a judgment matrix A, and the corresponding element of the judgment matrix A is aijAs shown in formula (4):
Figure RE-GDA0002783794870000051
wherein i is 1, …, n; j is 1, …, n; n is the total number of the single-layer indexes, and the maximum eigenvalue lambda of the judgment matrix A is calculated by utilizing the power methodmaxAnd a feature vector ξ;
(d) consistency check
If the relative consistency CR of the judgment matrix A is less than 0.1, the judgment matrix A is considered to be feasible, and the smaller the CR value is, the better the consistency inspection is; if the CR of the judgment matrix A is not less than 0.1 and does not pass the consistency test, returning to the step (b) to compare two by two again, and reconstructing a qualified judgment matrix A, wherein the calculation of the CR is as follows:
Figure RE-GDA0002783794870000052
TABLE 3 degree of freedom index RI
Figure RE-GDA0002783794870000053
From A xi ═ lambdamaxwIAHPDetermining weights w of indexes in a hierarchical modelIAHP；
(e) According to tables 4 to 11, each level of the DC protection system is divided into two according to an interval level method, then, experts are selected to compare and score the index levels of the power user equipment, interval weights of each level are obtained by adopting an interval level analysis method, characteristic vectors are obtained according to judgment matrixes of a layer B when a scheme layer C is aligned, a total level sorting weight is established, and a total sorting interval weight of the scheme layer C is obtained; the algorithm of the total sort interval weight of the scheme layer C is shown in table 12:
TABLE 4 DC PROTECTION SYSTEM A1Hierarchy
Figure RE-GDA0002783794870000054
Figure RE-GDA0002783794870000061
TABLE 5 measurement device B1Hierarchy
Figure RE-GDA0002783794870000062
TABLE 6 measurement interface B2Hierarchy
Figure RE-GDA0002783794870000063
TABLE 7 hostForm B3Hierarchy
Figure RE-GDA0002783794870000064
TABLE 8 device model B4Hierarchy
Figure RE-GDA0002783794870000065
TABLE 9 independent form B5Hierarchy
Figure RE-GDA0002783794870000066
Watch 10 two-out-of-three device B6Hierarchy
Figure RE-GDA0002783794870000071
Meter 11 trip outlet and secondary loop B7Hierarchy
Figure RE-GDA0002783794870000072
Table 12 Total sorting weight Algorithm
Figure RE-GDA0002783794870000073
Further, the specific implementation process of the step two is as follows:
classifying the direct current protection system according to different levels and different positions of defects of the direct current protection system, counting all defect data of each converter station in the past, setting the score of each type of equipment as 100, and subtracting the data obtained by the ratio as the type of element CiIf this element C is a fraction ofiDefect-free, then the element C is determinediIs 100, thereby obtaining each of the DC protection systemsThe value of each layer is set as the bottom of each element of the direct current protection system, and the solution layer C is combinediObtaining the interval value of the interval hierarchical method of the direct current protection system; when a certain converter station is analyzed, according to the bottom value of the direct current protection system and the defect frequency of the converter station in the previous year, the fault frequency is subjected to division reduction on the basis of bottom division; and similarly, the interval value of each criterion layer is obtained according to the weight of each criterion layer B and the value of each layer of the direct-current protection system, so that the operation quality of various devices can be judged.
Further, the specific implementation process of the step three is as follows:
classifying the direct current protection system according to different levels and different positions of the defects of the direct current protection system in the step one, collecting historical data of the defects of the direct current protection system, and obtaining three index data of the defects of the direct current protection system according to the historical data: the component defect-free time rate MTTF, the component mean repair time rate MTTR and the component defect rate EDR are calculated by the following method:
1) MTTF (maximum time transfer frequency) without defect of element
Figure RE-GDA0002783794870000081
Figure RE-GDA0002783794870000082
In the formula, TND、TTNORespectively, the non-defective running time and the normal running total time of a certain type of element; n is a radical ofDThe defect times of the type of element in a given time period are taken as the defect times; sigma TNDThe defect-free time sum of various elements is obtained;
2) mean Time To Repair (MTTR) of components
Figure RE-GDA0002783794870000083
Figure RE-GDA0002783794870000084
In the formula, TAC、TTDMRespectively the average maintenance time and the total maintenance time of certain types of elements; n is a radical ofRThe number of times of repairing the defects of the elements in a given time period; sigma TACThe average maintenance time sum of various components;
3) element defect rate EDR
Figure RE-GDA0002783794870000091
In the formula, ND、NTThe defect times of certain type of elements in a given time period and the total defect times of the direct current transmission protection system are respectively.
Further, the specific implementation process of the step four is as follows:
(a) based on the three index quantities in step three, form (x)ij)m×nThe relation matrix (a) is obtained by preprocessing the data to obtain a standard matrixij)m×nM is the number of types of the defect elements, n is the number of the defect indexes, and the index data in the matrix is divided into a forward index, a reverse index and an interval index, wherein the forward index is the better the evaluation value is, the reverse index is the better the evaluation value is, the interval index is the better the evaluation value is closer to the middle part of the interval, and the processing of the index data is as follows:
the forward direction index is consistent with the formula
Figure RE-GDA0002783794870000092
The consistency processing formula of the reverse indexes is
Figure RE-GDA0002783794870000093
The consistent processing formula of the interval index is
Figure RE-GDA0002783794870000094
In the formula, xmaxThe maximum value in the same index data; x is the number ofminThe minimum value in the same index data; x is the number ofmidThe intermediate values in the same type of index data;
(b) and (3) solving the characteristic weight, and solving the characteristic weight of the jth index of the ith equipment according to the standard matrix:
Figure RE-GDA0002783794870000101
when P is presentijWhen 0, it makes no sense to calculate the entropy value, so for PijTo correct it[10]Is defined as
Figure RE-GDA0002783794870000102
Calculating entropy of jth index
Figure RE-GDA0002783794870000103
In the formula: j is 1,2, …, n;
(c) calculating the entropy weight of the jth index
Figure RE-GDA0002783794870000104
In the formula: j is 1,2, …, n;
setting two weight values in an entropy weight method, wherein firstly, the numerical value of an element without defect data in the element is unchanged after the element is converted into a standard matrix; secondly, the numerical value of the elements without defect data in the elements is set to be 0.5 after the elements are converted into the standard matrix, because some elements do not collect defect data in the past, the elements have no data, and therefore the numerical value of the elements after the defect conditions of the elements are converted into the standard matrix is set to be 0.5;
(f) weight coefficient of each type of element
According to the two sets of weight values, two sets of weight coefficients are calculated according to the formula (18), and the weight coefficient Z of each type of element of the entropy weight direct current power transmission protection system in the defect dataiThe calculation method is as follows:
Figure RE-GDA0002783794870000105
in the formula: i is 1,2, …, m;
(g) calculating a score for an entropy weight method of a DC protection system
And (4) obtaining the weight coefficients of various elements of the two groups of direct current power transmission protection systems in the defect data according to the formula (18), and combining the base scores of the various elements of the direct current protection systems obtained in the step two to obtain the interval score of the entropy weight method of the direct current protection systems.
Furthermore, in the second, fourth and fifth steps, the score is over 85 points, which shows that the system state is excellent; a score of between 70 and 85 indicates a good system status; a score of between 60 and 70, indicating that the system is in condition; a score of 60 or less indicates poor system status.
Furthermore, in the fifth step, the value of theta is 0.5, and the analytic hierarchy process and the entropy weight process respectively account for half of the specific gravity.
The method comprehensively summarizes the protection defect conditions of high-voltage direct-current transmission in recent years, based on the defect types, data and evaluation indexes of the converter stations, and by combining the actual conditions of the converter stations, the direct-current protection system of the converter stations and various protection devices are evaluated by adopting an interval analytic hierarchy process and an interval entropy weight mixed weighting method, and according to the scores, which device has a good running state can be seen, so that the method plays a decisive role in improving the reliability of the relay protection device; the method not only can find weak links needing to be improved of the direct current protection system to obtain qualitative and quantitative conclusions, but also can scientifically and reasonably evaluate and analyze the operation conditions of the direct current protection system and various devices according to the method, has good application effect on solving the actual engineering problem, and can ensure that the relay protection capability is fully exerted in actual application.
Detailed Description
The present invention will be further described with reference to specific examples to fully understand the objects, features and effects of the present invention.
In this embodiment, the data is derived from the evaluation of the defect data from the north Hu converter station 2014 1 to 2019 6, and includes 5 converter stations such as Longquan converter station, Jianglin converter station, Tuyan converter station, Yidu converter station and Kudzuvin dam converter station.
The embodiment provides a method for evaluating reliability of a direct current protection system based on mixed weighting of an interval analytic hierarchy process and an interval entropy weight, which comprises the following steps of:
the method comprises the following steps: according to the DC protection system reliability comprehensive evaluation index system composition, constructing an index system hierarchy: and selecting experts to compare and score the index levels of the power user equipment in pairs by using a target level, a criterion level (first level index) and a scheme level (second and third level indexes), obtaining interval weights of all levels by adopting an inter-level analysis method, obtaining a characteristic vector according to a judgment matrix of a level B aligned by a scheme level C, establishing a total ranking weight of the level, and obtaining the total ranking weight of the scheme level C.
(1) Establishing each hierarchy of a DC protection system
Establishing each level of the direct current protection system according to the composition of the direct current protection system, and establishing a judgment matrix of an analytic hierarchy process, wherein each level of the direct current protection system is shown as the following table:
Figure RE-GDA0002783794870000121
Figure RE-GDA0002783794870000131
(2) DC protection system A1
Figure RE-GDA0002783794870000132
Ranking the importance degree: device type independent three-taking-two device measuring interface trip outlet and secondary circuit
The discrimination matrix (interval width of the previous column in the column labeled yellow, the same as below):
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
2.0000 0.2000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
4.0000 0.1000 3.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
4.0000 0.1000 3.0000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
4.0000 0.1000 3.0000 0.1000 1.0000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
3.0000 0.3000 2.0000 0.3000 0.5000 0.3000 0.3333 0.2000 0.3333 0.2000 1.0000 0.0000 0.0000 0.0000
0.2000 0.3000 0.3000 0.3000 0.2000 0.2500 0.2000 0.2500 0.2000 0.2000 0.3333 0.2000 1.0000 0.0000
the results met the consistency test as follows:
DC protection system A1 Interval low value Interval high value
Measuring apparatus B1 0.0661 0.0687
Measurement interface B2 0.0850 0.0851
Host type B3 0.2207 0.2226
Device type B4 0.2360 0.2387
Independent type B5 0.2355 0.2386
Two-out-of-three device B6 0.1136 0.1208
Trip outlet and secondary circuit B7 0.0339 0.0347
(3) Measuring apparatus B1
Figure RE-GDA0002783794870000141
Ranking the importance degree: photoelectric CT (pure light CT) and photoelectric PT (zero flux transformer) are conventional transformers
Discrimination matrix
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
1.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
1.0000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
1.0000 0.1000 1.0000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000
1.0000 0.1000 1.0000 0.1000 1.0000 0.1000 1.0000 0.1000 1.0000 0.0000
The results met the consistency test as follows:
measuring apparatus B1 Interval low value Interval high value
Photoelectric CTC1 0.2008 0.2008
Pure light type CTC2 0.2004 0.2004
Photoelectric PTC3 0.2000 0.2000
Zero flux transformer C4 0.1996 0.1996
Conventional mutual inductor C5 0.1992 0.1992
(4) Measurement interface B2
Figure RE-GDA0002783794870000142
Ranking the importance degree: processor board, power supply module, on-off interface board, analog interface board, communication board and case back board
Discrimination matrix
Figure RE-GDA0002783794870000143
Figure RE-GDA0002783794870000151
The results met the consistency test as follows:
measurement interface B2 Interval low value Interval high value
Processor board C6 0.4182 0.4394
Switching value interface board C7 0.1092 0.1162
Analog quantity interface board C8 0.1086 0.1153
Communication board C9 0.1079 0.1145
Power supply module C10 0.1854 0.1889
Chassis backboard C11 0.0469 0.0496
(5) Host type DC protection system B3
Figure RE-GDA0002783794870000152
Ranking the importance degree: CPU board, power module, PCI board and software
Discrimination matrix
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.5000 0.2000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
1.0000 0.1000 2.0000 0.1000 1.0000 0.0000 0.0000 0.0000
0.3333 0.2000 0.3330 0.1000 0.3330 0.1000 1.0000 0.0000
The results met the consistency test as follows:
host type B3 Interval low value Interval high value
CPU board C12 0.3402 0.3518
PCI board C13 0.2074 0.2123
Power supply module C14 0.3441 0.3477
Software C15 0.0967 0.0998
(6) Device type DC protection system B4
Figure RE-GDA0002783794870000153
Ranking the importance degree: processor board, power supply module, on-off interface board, analog interface board, communication board and case back board
Discrimination matrix
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.2500 0.2000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.2500 0.2000 1.0000 0.2000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.2500 0.2000 1.0000 0.2000 1.0000 0.2000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.3333 0.2000 2.0000 0.2000 2.0000 0.2000 2.0000 0.2000 1.0000 0.0000 0.0000 0.0000
0.1667 0.2000 0.3333 0.2000 0.3333 0.2000 0.3333 0.2000 0.3333 0.2000 1.0000 0.0000
The results met the consistency test as follows:
device type B4 Interval low value Interval high value
Processor plug-in C16 0.4182 0.4394
Communication plug-in component C17 0.1092 0.1162
Opening-in and opening-out plug-in component C18 0.1086 0.1153
Power supply plug-in component C19 0.1079 0.1145
Chassis backboard C20 0.1854 0.1889
Software C21 0.0469 0.0496
(7) Independent type direct current protection system B5
Figure RE-GDA0002783794870000161
Ranking the importance degree: processor board, power module, MMI plug-in, A/D (VFC) plug-in, communication plug-in and AC plug-in software
Discrimination matrix
Figure RE-GDA0002783794870000162
Figure RE-GDA0002783794870000171
The results met the consistency test as follows:
independent type B5 Interval low value Interval high value
Processor plug-in C22 0.2791 0.3044
Opening-in and opening-out plug-in component C23 0.0954 0.0980
MMI plug-in component C24 0.0937 0.1002
A/D (VFC) plug-in component C25 0.0933 0.0996
Communication plug-in component C26 0.0929 0.0991
AC plug-in C27 0.0924 0.0985
Power supply plug-in component C28 0.1716 0.1722
Software C29 0.0539 0.0557
(8) Two-out-of-three device B6
Figure RE-GDA0002783794870000172
Ranking the importance degree: processor board card (independent module) power board card (communication board)
Discrimination matrix
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.3333 0.1500 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.3333 0.1500 2.0000 0.1500 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.5000 0.2000 1.0000 0.1000 0.5000 0.2000 1.0000 0.0000 0.0000 0.0000
1.0000 0.1000 3.0000 0.1000 2.0000 0.1000 2.0000 0.2000 1.0000 0.0000
The results met the consistency test as follows:
two-out-of-three device B6 Interval low value Interval high value
Processor board card C30 0.3222 0.3245
Open-in and open-out board card C31 0.0989 0.1000
Power supply board card C32 0.1672 0.1681
Communication board C33 0.1170 0.1220
Independent module C34 0.2889 0.2913
(9) Trip outlet and secondary circuit B7
Figure RE-GDA0002783794870000173
Figure RE-GDA0002783794870000181
Ranking the importance degree: switch operation box cable loop or optical fiber loop or communication channel accessory
Discrimination matrix
1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.5000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.5000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.5000 0.1000 1.0000 0.1000 1.0000 0.1000 1.0000 0.0000 0.0000 0.0000
0.3333 0.3000 0.5000 0.2000 0.5000 0.2000 0.5000 0.2000 1.0000 0.0000
The results met the consistency test as follows:
Figure RE-GDA0002783794870000182
and obtaining a characteristic vector according to the judgment matrix of the layer B when the scheme layer C is aligned, establishing a hierarchical total sorting weight, and obtaining a total sorting interval weight of the scheme layer C. The algorithm of the total sorting weight of the scheme layer C is as follows:
Figure RE-GDA0002783794870000183
Figure RE-GDA0002783794870000191
the total sort weight of the scheme layer C is obtained according to the table and the related data as follows:
Figure RE-GDA0002783794870000192
Figure RE-GDA0002783794870000201
step two: classifying the direct current protection system according to different levels and different positions of the defects of the direct current protection system, counting all defect data of each converter station in the past, setting scores of all levels of the direct current protection system according to the defect data, and combining with a scheme layer CiAnd obtaining the interval value of the interval hierarchical method of the direct current protection system. And obtaining the interval value of each criterion layer according to the weight of each criterion layer B and the value of each layer of the direct current protection system, thereby judging the operation quality of various devices.
The invention is based on the defect data of Hubei converter stations from 1 month to 2019 months in 2014 to carry out evaluation, and the defect data comprises 5 converter stations such as a Longquan converter station, a Jiangling converter station, a Tuoling converter station, a Yidu converter station, a Kuzhou dam converter station and the like, and specific defect information is shown in tables 1 to 6.
Table 1 Longquan convertor station defect information
Figure RE-GDA0002783794870000202
TABLE 2 Jiangling converter station Defect information
Figure RE-GDA0002783794870000203
Figure RE-GDA0002783794870000211
TABLE 3 Tuolin converter station Defect information
Figure RE-GDA0002783794870000212
Table 4 defect information for preferably all converter stations
Figure RE-GDA0002783794870000213
TABLE 5. Kudzuvine river dam converter station Defect information
Figure RE-GDA0002783794870000214
Figure RE-GDA0002783794870000221
TABLE 6 Defect information of five converter stations
Figure RE-GDA0002783794870000222
And (4) counting all defect data of five convertor stations in Hubei in nearly five years, and calculating the proportion of the defects of each type of equipment to the total defects. As shown in the table below.
Figure RE-GDA0002783794870000223
Figure RE-GDA0002783794870000231
And then, setting the score of each type of equipment as 100, subtracting the data obtained by percentage as the score of the element, and if the element is not defective, determining the score of the element as 100. Specific scores are shown in the following table.
Figure RE-GDA0002783794870000232
Figure RE-GDA0002783794870000241
According to the scores of all the elements in the table, the total sorting weight of the sub-criterion layer C is combined to obtain the direct current protection system A1Has a section value of [93.9370,99.1927 ]]。
Similarly, the interval values of various devices are calculated as follows, and the following table shows that the host dc protection system device and the measuring device have lower reliability than other devices and need to pay attention.
Figure RE-GDA0002783794870000242
Step three: classifying the direct current protection system according to different levels of the direct current protection system and different positions of the defects, collecting historical data of the defects of the direct current protection system, and obtaining index data of the defects of the direct current protection system according to the historical data, wherein the index data are respectively as follows: a component defect free time rate MTTF, a component mean repair time rate MTTR, and a component defect rate EDR.
Counting all defect data of five convertor stations in Hubei in the last five years, and calculating various index values of each type of equipment by combining the actual operation conditions of the five convertor stations as follows:
Figure RE-GDA0002783794870000243
Figure RE-GDA0002783794870000251
step four: and carrying out an entropy weight method on the three index quantities to obtain the weights of the three indexes. Since the first two of the three defect indexes are reverse indexes, the smaller the index value is, the better the index value is, the larger the latter index is, the better the index value is. In order to obtain the interval value of the direct current protection system, two weighted values of the entropy weight method are calculated, wherein the value of the zero-defect data element in the element is unchanged after being converted into a standard matrix, and the value of the zero-defect data element in the element is set to be 0.5 after being converted into the standard matrix, so that the interval value of the entropy weight method of the whole direct current protection system is obtained.
The above table data was preprocessed according to the entropy weight method to obtain the standard matrix as shown below.
Figure RE-GDA0002783794870000261
Figure RE-GDA0002783794870000271
Since the first two of the three defect indexes are reverse indexes, the smaller the index value is, the better the index value is, the larger the latter index is, the better the index value is. In order to obtain the interval value of the direct current protection system, two weighted values of the entropy weight method are calculated, wherein the value of the zero-defect data element in the element is unchanged after being converted into a standard matrix, and the value of the zero-defect data element in the element is set to be 0.5 after being converted into the standard matrix, so that the interval value of the entropy weight method of the whole direct current protection system is obtained.
According to the calculation method of the entropy weight method, the weight of each index is obtained as follows
Figure RE-GDA0002783794870000272
The calculated weights for each element are as follows:
Figure RE-GDA0002783794870000273
Figure RE-GDA0002783794870000281
finally, the interval value of the entropy weight method of the whole direct current protection system is obtained to be [86.1057, 90.8727 ].
For the evaluation of various devices of the direct current protection system, the same method is adopted, the targeted data only contains the respective data of various devices, the data of other devices is not considered, if the devices have no defect data, the operation condition of the devices is determined to be good, the evaluation value is 100, if the devices have only one element, the defect value of other elements of the devices is determined to be 1 percent of the element value, the average maintenance time of the elements is calculated according to 0.5, the operation time of the elements is calculated according to 5 years, the evaluation value of the devices is obtained according to the value, and the interval value of the entropy method of various devices of the direct current protection system is obtained by the calculation of the method is shown in the following table.
Figure RE-GDA0002783794870000282
Figure RE-GDA0002783794870000291
Step five: in order to take the subjective advantage of the interval hierarchy analysis method and the objective advantage of the interval entropy weight method into consideration, the advantages of the interval hierarchy analysis method and the interval entropy weight method are taken, the weight factor theta is introduced, and the interval hierarchy analysis method and the interval entropy weight method are combined to obtain the index evaluation comprehensive interval weight w. And when theta is 1 or 0, the comprehensive interval weight w respectively corresponds to an interval AHP method and an interval entropy weight method, so that the evaluation of the direct current protection system is obtained, and the direct current protection system is evaluated.
DC protection system A1Has a section value of [93.9370,99.1927 ]]The interval value of the entropy weight method of the direct current protection system is [86.105 ]7,90.8727]Setting theta to be 0.5, and obtaining interval values of [90.0215,95.0327 ] of the direct current protection system]。
The scores of the various devices of the dc protection system are shown in the following table.
Figure RE-GDA0002783794870000292
As can be seen from the above table, the overall reliability of the dc protection system is excellent, but the reliability of the measurement device is lower than that of other devices, and the maintenance is important.
The invention focuses on counting the faults and defect conditions of equipment such as a protection device, a board card, a secondary circuit, an optical CT and the like in a direct current control protection system, comprehensively summarizes the protection defect conditions of high-voltage direct current transmission in recent years, and based on the defect types, data and evaluation indexes of converter stations, the method of mixed weighting of interval analytic hierarchy process and interval entropy weight is adopted to evaluate the direct current protection system of the converter stations and various protection equipment by combining the actual conditions of the converter stations.
The method can not only find weak links needing to be improved of the direct current protection system to obtain qualitative and quantitative conclusions, but also scientifically and reasonably evaluate and analyze the operation conditions of the direct current protection system and various devices according to the method, has good application effect on solving the practical engineering problem, and can ensure that the relay protection capability is fully exerted in practical application. The use effect of the method can be obtained through a specific implementation scheme, and the method provided by the embodiment finds weak links needing to be improved and obtains qualitative and quantitative conclusions, so that the method has the advantages of effectiveness and feasibility. The invention is also suitable for the reliability analysis of the alternating current system.

Claims (7)

1.一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于包括以下步骤：1. a DC protection system reliability evaluation method based on interval AHP and interval entropy weight mixed weighting, is characterized in that comprising the following steps: 步骤一：根据直流保护系统可靠性综合评估指标体系组成，构建指标体系层次：目标层、准则层B以及方案层C，然后对电力用户设备的指标层次进行两两比较打分，采用区间层次分析法得出各个层次区间权重，根据方案层C对准则层B的判断矩阵得出特征向量，建立层次总排序权值，得出方案层C的总排序权值；Step 1: According to the composition of the comprehensive evaluation index system for the reliability of the DC protection system, construct the index system layers: the target layer, the criterion layer B and the scheme layer C, and then compare and score the indicator layers of the power user equipment, using the interval analysis hierarchy process. Obtain the interval weights of each level, obtain the eigenvectors according to the judgment matrix of the scheme layer C to the criterion layer B, establish the total ranking weights of the layers, and obtain the total sorting weights of the scheme layer C; 步骤二：根据直流保护系统的各个层次以及缺陷所在部位不同，对直流保护系统进行分类，统计以往各换流站的全部缺陷数据，根据缺陷数据设置直流保护系统各个层次的分值，此分值设置为直流保护系统各个元件的底分，结合方案层C的总排序权值，得到直流保护系统的区间层次法的评分；Step 2: According to the different levels of the DC protection system and the location of the defect, classify the DC protection system, count all the defect data of each converter station in the past, and set the score of each level of the DC protection system according to the defect data. It is set as the bottom score of each component of the DC protection system, and combined with the total ranking weight of the scheme layer C, the score of the interval level method of the DC protection system is obtained; 步骤三：根据直流保护系统的各个层次以及缺陷所在部位不同，对直流保护系统进行分类，收集直流保护系统缺陷的历史数据，根据历史数据得出直流保护系统缺陷的三个指标数据：元件无缺陷时间率MTTF、元件平均维修时间率MTTR以及元件缺陷率EDR；Step 3: According to the different levels of the DC protection system and the location of the defect, classify the DC protection system, collect the historical data of the DC protection system defect, and obtain three index data of the DC protection system defect according to the historical data: no defects in the components Time rate MTTF, component mean repair time rate MTTR and component defect rate EDR; 步骤四：对步骤三中的三个指标进行熵权法，得出三个指标的权重，在熵权法中设置两个权重值，一是元件中无缺陷数据的元件换算成标准矩阵后数值不变，二是元件中无缺陷数据的元件换算成标准矩阵后数值设置为0.5，由此得出整个直流保护系统熵权法的评分；Step 4: Perform the entropy weight method on the three indicators in step 3 to obtain the weight of the three indicators, and set two weight values in the entropy weight method. First, the components without defect data in the components are converted into standard matrix values. The second is that the component with no defect data in the component is converted into a standard matrix and the value is set to 0.5, from which the entropy weight method score of the entire DC protection system is obtained; 步骤五：引入权重因子θ，将区间层次分析法和区间熵权法相结合，得到指标评价综合区间评分w，综合区间评分w随着θ的改变而改变，指标评价综合区间评分w的计算公式如下：Step 5: Introduce the weight factor θ, combine the interval AHP and the interval entropy weight method to obtain the comprehensive interval score w of the index evaluation. The comprehensive interval score w changes with the change of θ. The calculation formula of the comprehensive interval score w of the index evaluation is as follows : w＝θwIAHP+(1-θ)wIEN,0≤θ≤1 (19)w=θw IAHP +(1-θ)w IEN , 0≤θ≤1 (19) 其中wIAHP为步骤二所得区间层次法的评分，wIEN为步骤四所得步骤四，当θ＝1或者0时，分别对应于区间层次法和区间熵权法。Wherein w IAHP is the score of the interval hierarchy method obtained in the second step, and w IEN is the step 4 obtained in the step 4. When θ=1 or 0, it corresponds to the interval hierarchy method and the interval entropy weight method respectively. 2.根据权利要求1所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤一的具体实现过程如下：2. a kind of DC protection system reliability evaluation method based on interval analytic hierarchy process and interval entropy weight mixed weighting according to claim 1, is characterized in that: the concrete realization process of step 1 is as follows: (a)根据直流保护系统的组成建立直流保护系统的各个层次，构建层次分析法的判断矩阵，直流保护系统的各个层次如下表1所示，其中直流保护系统(A1)为目标层，各装置、设备(Bj)为准则层，各接口板、CT、PT(Ci)为方案层；(a) According to the composition of the DC protection system, establish the various levels of the DC protection system, and construct the judgment matrix of the AHP. The various levels of the DC protection system are shown in Table 1 below, in which the DC protection system (A 1 ) is the target layer, and each The device and equipment (B j ) are the criterion layer, and each interface board, CT and PT (C i ) are the solution layer; 表1直流保护系统的各层次Table 1 Various levels of DC protection system
Figure RE-FDA0002783794860000021
Figure RE-FDA0002783794860000021
Figure RE-FDA0002783794860000031
Figure RE-FDA0002783794860000031
 (b)根据直流保护系统的各层次，选择专家对电力用户设备的指标层次进行两两比较打分，电力用户指标两两比较打分按照表2“区间标度表”进行，打分结果一般由区间中aij点和宽度μ两部分组成，其中，区间中点aij的取值，是判断的一个基数，或随机变量在判断区间标度内的均值；宽度μ是在确定好基数之后，专家根据判断的不确定性以及模糊性而给定的基数aij的区间变化范围；(b) According to the various levels of the DC protection system, select experts to compare and score the index levels of the power user equipment. The a ij point and the width μ are composed of two parts, among which, the value of the point a ij in the interval is a basis for judgment, or the mean value of the random variable within the scale of the judgment interval; the width μ is the basis after the basis is determined. The range of interval variation of the given base number a ij due to the uncertainty and ambiguity of judgment;
Figure RE-FDA0002783794860000032
Figure RE-FDA0002783794860000032
 表2区间比例标度表Table 2 Interval proportional scale table
Figure RE-FDA0002783794860000033
Figure RE-FDA0002783794860000033
 根据区间比例标度表打分，形成区间数
Figure RE-FDA0002783794860000034
当指标i比指标j重要时，即aij≥1，i≠j时Score according to the interval proportional scale table to form the interval number
Figure RE-FDA0002783794860000034
When the index i is more important than the index j, that is, when a ij ≥ 1, i≠j
Figure RE-FDA0002783794860000041
Figure RE-FDA0002783794860000041
 当指标j比指标i重要时，即aji≥1，j≠i时When the index j is more important than the index i, that is, a ji ≥ 1, j≠i
Figure RE-FDA0002783794860000042
Figure RE-FDA0002783794860000042
 (c)构造区间判断矩阵A(c) Constructing the interval judgment matrix A 选定的专家根据“区间比例标度法”对同层的每一个指标相对于上一层指标重要性进行两两比较打分，构建判断矩阵A，其对应元素为aij，如式(4)所示：According to the "interval proportional scaling method", the selected experts compare and score the importance of each index of the same layer relative to the index of the previous layer, and construct a judgment matrix A, whose corresponding element is a ij , as shown in formula (4) shown:
Figure RE-FDA0002783794860000043
Figure RE-FDA0002783794860000043
 其中i＝1,…,n；j＝1,…,n；n是指单层指标的个数总数，利用幂法，计算判断矩阵A的最大特征值λmax及特征向量ξ；where i=1,...,n; j=1,...,n; n refers to the total number of single-layer indicators, and the power method is used to calculate the maximum eigenvalue λmax and eigenvector ξ of the judgment matrix A; (d)一致性检验(d) Consistency check 若判断矩阵A的相对一致性CR小于0.1，则认为判断矩阵A可行，通过一致性检验，通常CR值越小越好；若判断矩阵A的CR不满足小于0.1，没有通过一致性检验，则退回步骤(b)重新两两比较，重新构造合格的判断矩阵A，CR的计算如公式如下：If the relative consistency CR of the judgment matrix A is less than 0.1, then the judgment matrix A is considered feasible, and the consistency test is passed. Generally, the smaller the CR value, the better; Return to step (b) to compare the two pairs again, and reconstruct the qualified judgment matrix A. The calculation of CR is as follows:
Figure RE-FDA0002783794860000044
Figure RE-FDA0002783794860000044
 表3自由度指标RITable 3 Degree of freedom index RI
Figure RE-FDA0002783794860000045
Figure RE-FDA0002783794860000045
 由Aξ＝λmaxwIAHP确定层次结构模型中各指标的权重wIAHP；Determine the weight w IAHP of each index in the hierarchical structure model by Aξ=λ max w IAHP ; (e)按照表4到表11对直流保护系统的各个层次按区间层次法，然后选择专家对电力用户设备的指标层次进行两两比较打分，采用区间层次分析法得出各个层次区间权重，根据方案层C对准则层B的判断矩阵得出特征向量，建立层次总排序权值，得出方案层C的总排序区间权值；方案层C的总排序区间权值的算法如表12所示：(e) According to Table 4 to Table 11, use the interval level method for each level of the DC protection system, and then select experts to compare and score the indicator levels of the power user equipment, and use the interval analysis hierarchy process to obtain the interval weights of each level. Scheme layer C obtains eigenvectors from the judgment matrix of criterion layer B, establishes the total ranking weight of the hierarchy, and obtains the total sorting interval weight of scheme layer C; the algorithm of the total sorting interval weight of scheme layer C is shown in Table 12 : 表4直流保护系统A1层次Table 4 DC protection system A 1 level
Figure RE-FDA0002783794860000051
Figure RE-FDA0002783794860000051
 表5测量设备B1层次Table 5 Measuring Equipment B Level 1
Figure RE-FDA0002783794860000052
Figure RE-FDA0002783794860000052
 表6测量接口B2层次Table 6 Measurement Interface B Level 2
Figure RE-FDA0002783794860000053
Figure RE-FDA0002783794860000053
 表7主机型B3层次Table 7 Mainframe type B 3 levels
Figure RE-FDA0002783794860000054
Figure RE-FDA0002783794860000054
 表8装置型B4层次Table 8 Device Type B 4 Tiers
Figure RE-FDA0002783794860000061
Figure RE-FDA0002783794860000061
 表9独立型B5层次Table 9 Independent Type B 5 Tiers
Figure RE-FDA0002783794860000062
Figure RE-FDA0002783794860000062
 表10三取二装置B6层次Table 10 Three out of two device B 6 levels
Figure RE-FDA0002783794860000063
Figure RE-FDA0002783794860000063
 表11跳闸出口及二次回路B7层次Table 11 Trip outlet and secondary circuit B 7 levels
Figure RE-FDA0002783794860000064
Figure RE-FDA0002783794860000064
 表12总排序权值的算法Table 12 Algorithms for total ranking weights
Figure RE-FDA0002783794860000065
Figure RE-FDA0002783794860000065
Figure RE-FDA0002783794860000071
Figure RE-FDA0002783794860000071
 3.根据权利要求1或2所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤二的具体实现过程如下：3. a kind of DC protection system reliability evaluation method based on interval analytic hierarchy process and interval entropy weight mixed weighting according to claim 1 and 2, is characterized in that: the concrete realization process of step 2 is as follows: 根据直流保护系统的各个层次以及缺陷所在部位不同，对直流保护系统进行分类，统计以往各换流站的全部缺陷数据，设每类设备分数是100分，减去占比后得出的数据为这类元件Ci的分数，若这个元件Ci无缺陷，则定这个元件Ci的分数为100，由此得到直流保护系统各个层次的分值，这个分值设置为直流保护系统各个元件的底分，结合方案层Ci的总排序权值，得到直流保护系统的区间层次法的区间值；分析某个换流站时，根据直流保护系统的底值，再根据这个换流站前一年的缺陷次数，故障次数在底分的基础上再进行减分制；同理根据各准则层B的权重以及直流保护系统各个层次的分值得出各准则层的区间分值，由此可以判断出各类设备运行的好坏。According to the different levels of the DC protection system and the location of the defect, the DC protection system is classified, and all the defect data of each converter station in the past are counted. The score of each type of equipment is set to 100 points, and the data obtained after deducting the proportion is The score of this type of component C i , if this component C i has no defects, the score of this component C i is set as 100, and the score of each level of the DC protection system is obtained, and this score is set as the value of each component of the DC protection system. The bottom score, combined with the total sorting weight of the scheme layer C i , can obtain the interval value of the interval hierarchy method of the DC protection system; when analyzing a certain converter station, according to the bottom value of the DC protection system, and then according to the previous value of the converter station. The number of defects in a year, and the number of faults will be deducted based on the bottom score. Similarly, the interval score of each criterion layer can be determined according to the weight of each criterion layer B and the score of each layer of the DC protection system. Check out the performance of various equipments. 4.根据权利要求1或2所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤三的具体实现过程如下：4. a kind of DC protection system reliability evaluation method based on interval analytic hierarchy process and interval entropy weight mixed weighting according to claim 1 and 2, is characterized in that: the concrete realization process of step 3 is as follows: 按照步骤一中直流保护系统的各个层次以及缺陷所在部位不同，对直流保护系统进行分类，收集直流保护系统缺陷的历史数据，根据历史数据得出直流保护系统缺陷的三个指标数据：元件无缺陷时间率MTTF、元件平均维修时间率MTTR以及元件缺陷率EDR，计算方法如下：According to the different levels of the DC protection system and the location of the defects in step 1, classify the DC protection system, collect the historical data of the DC protection system defects, and obtain three index data of the DC protection system defects according to the historical data: no defects in the components Time rate MTTF, component mean repair time rate MTTR and component defect rate EDR are calculated as follows: 1)元件无缺陷时间率MTTF1) Component defect-free time rate MTTF
Figure RE-FDA0002783794860000081
Figure RE-FDA0002783794860000081
Figure RE-FDA0002783794860000082
Figure RE-FDA0002783794860000082
 式中，TND、TTNO分别为某类元件的无缺陷运行时间和正常运行总时间；ND为给定时间段内该类元件的缺陷次数；∑TND为各类元件无缺陷时间总和；In the formula, T ND and T TNO are the defect-free running time and the total normal running time of a certain type of components, respectively; ND is the number of defects of this type of components in a given time period; ∑T ND is the sum of the defect-free time of various components ; 2)元件平均维修时间率MTTR2) Component mean maintenance time rate MTTR
Figure RE-FDA0002783794860000083
Figure RE-FDA0002783794860000083
Figure RE-FDA0002783794860000084
Figure RE-FDA0002783794860000084
 式中，TAC、TTDM分别为某类元件的平均维修时间和维修总时间；NR为给定时间段内该类元件缺陷的修复次数；∑TAC为各类元件平均维修时间总和；In the formula, T AC and T TDM are the average maintenance time and total maintenance time of a certain type of components, respectively; NR is the number of repairs for defects of this type of component in a given time period; ∑T AC is the sum of the average maintenance time of various components; 3)元件缺陷率EDR3) Component defect rate EDR
Figure RE-FDA0002783794860000085
Figure RE-FDA0002783794860000085
 式中，ND、NT分别为给定时间段内的某类元件缺陷次数和直流输电保护系统缺陷总次数。In the formula, N D and N T are the number of defects of a certain type of components and the total number of defects of the DC transmission protection system in a given time period, respectively. 5.根据权利要求1所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤四的具体实现过程如下：5. a kind of DC protection system reliability evaluation method based on interval AHP and interval entropy weight mixed weighting according to claim 1, is characterized in that: the concrete realization process of step 4 is as follows: (a)基于步骤三中的三个指标量，形成(xij)m×n的关系矩阵，首先对其进行数据预处理，求出标准矩阵(aij)m×n，m为缺陷元件类别数，n为缺陷指标数，矩阵中的指标数据分为正向指标、逆向指标和区间指标，其中正向指标是评价值越大越好，逆向指标是评价值越小越好，区间指标是评价值越靠近区间中部越好，指标数据的处理如下：(a) Based on the three index quantities in step 3, a relationship matrix of (x ij ) m×n is formed, and data preprocessing is performed first to obtain a standard matrix (a ij ) m×n , where m is the defect component category Number, n is the number of defect indicators, the indicator data in the matrix is divided into positive indicators, reverse indicators and interval indicators, where the positive indicator is that the larger the evaluation value, the better, the reverse indicator is that the smaller the evaluation value, the better, and the interval indicator is the evaluation value. The closer the value is to the middle of the interval, the better. The indicator data is processed as follows: 正向指标的一致化处理公式为The consistent processing formula of the positive index is:
Figure RE-FDA0002783794860000091
Figure RE-FDA0002783794860000091
 逆向指标的一致化处理公式为The consistent processing formula of the inverse index is:
Figure RE-FDA0002783794860000092
Figure RE-FDA0002783794860000092
 区间指标的一致化处理公式为The consistent processing formula of the interval index is:
Figure RE-FDA0002783794860000093
Figure RE-FDA0002783794860000093
 式中，xmax为同一类指标数据中的最大值；xmin为同一类指标数据中的最小值；xmid为同一类指标数据中的中间值；In the formula, x max is the maximum value in the same type of indicator data; x min is the minimum value in the same type of indicator data; x mid is the middle value in the same type of indicator data; (b)求取特征权重，依据标准矩阵，求出第i个设备的第j个指标的特征权重：(b) Find the feature weight, according to the standard matrix, find the feature weight of the j-th index of the i-th device:
Figure RE-FDA0002783794860000094
Figure RE-FDA0002783794860000094
 当Pij＝0时，计算熵值无意义，因此对Pij加以修正[10]，定义为When P ij = 0, it is meaningless to calculate the entropy value, so P ij is modified [10] , which is defined as
Figure RE-FDA0002783794860000095
Figure RE-FDA0002783794860000095
 计算第j个指标的熵值Calculate the entropy value of the jth index
Figure RE-FDA0002783794860000096
Figure RE-FDA0002783794860000096
 式中：j＝1,2,…,n；In the formula: j=1,2,…,n; (c)计算第j个指标的熵权重(c) Calculate the entropy weight of the jth indicator
Figure RE-FDA0002783794860000101
Figure RE-FDA0002783794860000101
 式中：j＝1,2,…,n；In the formula: j=1,2,…,n; 在熵权法中设置两个权重值，一是元件中无缺陷数据的元件换算成标准矩阵后数值不变；二是元件中无缺陷数据的元件换算成标准矩阵后数值设置为0.5，这是由于有些元件以往未收集到缺陷数据，导致这些元件无数据，因此将这些元件的缺陷情况换算成标准矩阵后数值设置为0.5；Two weight values are set in the entropy weight method. One is that the components with no defect data in the components are converted into standard matrices and their values remain unchanged; Since some components have not collected defect data in the past, resulting in no data for these components, the value of these components is set to 0.5 after converting the defect status of these components into a standard matrix; (d)各类元件的权重系数(d) Weight coefficients of various components 根据上述设置的两组权重值，分别按式(18)计算出两组权重系数，熵权重直流输电保护系统各类元件在缺陷数据中的权重系数Zi计算方法如下：According to the two sets of weight values set above, the two sets of weight coefficients are calculated according to formula (18). The calculation method of the weight coefficient Z i of various components of the entropy-weighted DC transmission protection system in the defect data is as follows:
Figure RE-FDA0002783794860000102
Figure RE-FDA0002783794860000102
 式中：i＝1,2,…,m；In the formula: i=1,2,...,m; (e)计算直流保护系统熵权法的评分(e) Calculate the score of the entropy weight method of the DC protection system 根据式(18)得到的两组直流输电保护系统各类元件在缺陷数据中的权重系数，结合步骤二中得到的直流保护系统各个元件的底分，得出直流保护系统熵权法的区间评分。According to the weight coefficients of various components in the defect data of the two groups of HVDC transmission protection system obtained according to formula (18), combined with the bottom score of each component of the DC protection system obtained in step 2, the interval score of the DC protection system entropy weight method is obtained . 6.根据权利要求1所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤二、四、五中，评分在85分及以上，表明系统状态优；评分在70到85之间，表明系统状态良；评分在60到70之间，表明系统状态中；评分小于等于60，表明系统状态差。6. a kind of DC protection system reliability evaluation method based on interval analytic hierarchy process and interval entropy weight mixed weighting according to claim 1, is characterized in that: in step 2, 4, 5, the score is 85 points and above , indicating that the system state is excellent; the score is between 70 and 85, indicating that the system state is good; the score is between 60 and 70, indicating that the system state is in good condition; the score is less than or equal to 60, indicating that the system state is poor. 7.根据权利要求1所述的一种基于区间层次分析法和区间熵权混合赋权的直流保护系统可靠性评价方法，其特征在于：步骤五中θ取值为0.5，此时层次分析法和熵权法各占一半的比重。7. a kind of DC protection system reliability evaluation method based on interval analytic hierarchy process and interval entropy weight mixed weighting according to claim 1, is characterized in that: in step 5, the value of θ is 0.5, this time analytic hierarchy process and the entropy weight method each account for half of the proportion.
CN202011008701.1A 2020-09-23 2020-09-23 Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting Pending CN112183995A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011008701.1A CN112183995A (en) 2020-09-23 2020-09-23 Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011008701.1A CN112183995A (en) 2020-09-23 2020-09-23 Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting

Publications (1)

Publication Number Publication Date
CN112183995A true CN112183995A (en) 2021-01-05

Family

ID=73955411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011008701.1A Pending CN112183995A (en) 2020-09-23 2020-09-23 Direct-current protection system reliability evaluation method based on inter-zone analytic method and inter-zone entropy weight mixed weighting

Country Status (1)

Country Link
CN (1) CN112183995A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118094571A (en) * 2024-04-25 2024-05-28 南昌理工学院 Enterprise information security processing method, system, readable storage medium and computer

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933629A (en) * 2015-05-21 2015-09-23 天津大学 Power user equipment evaluation method based on interval level analysis and interval entropy combination
CN107563680A (en) * 2017-10-20 2018-01-09 广东电网有限责任公司电力科学研究院 A kind of distribution network reliability evaluation method based on AHP and entropy assessment
CN110580387A (en) * 2019-08-27 2019-12-17 国网湖北省电力有限公司电力科学研究院 Hybrid Weibull Reliability Evaluation Method for DC Protection System Based on Entropy Weight Method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104933629A (en) * 2015-05-21 2015-09-23 天津大学 Power user equipment evaluation method based on interval level analysis and interval entropy combination
CN107563680A (en) * 2017-10-20 2018-01-09 广东电网有限责任公司电力科学研究院 A kind of distribution network reliability evaluation method based on AHP and entropy assessment
CN110580387A (en) * 2019-08-27 2019-12-17 国网湖北省电力有限公司电力科学研究院 Hybrid Weibull Reliability Evaluation Method for DC Protection System Based on Entropy Weight Method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118094571A (en) * 2024-04-25 2024-05-28 南昌理工学院 Enterprise information security processing method, system, readable storage medium and computer

Similar Documents

Publication Publication Date Title
CN103971171B (en) A Method for State Evaluation of Power Transmission Equipment
CN110555596B (en) Sampling inspection strategy making method and system based on power distribution material quality evaluation
CN105740975B (en) A Method for Evaluation and Prediction of Equipment Defects Based on Data Correlation
CN104408562B (en) A kind of photovoltaic system generating efficiency comprehensive estimation method based on BP neural network
CN106651189B (en) Transformer state evaluation method based on multilayer composite rule
CN105956779A (en) Power transformer running state evaluation method and power transformer running state evaluation device
CN105719048A (en) Intermediate-voltage distribution operation state fuzzy integrated evaluation method based on principle component analysis method and entropy weight method
CN106384210A (en) Power transmission and transformation equipment maintenance priority ordering method based on maintenance risk premium
CN105023195A (en) Distributed photovoltaic power distribution network reliability evaluating method
CN102063657A (en) Operating level and power supplying capability evaluation method for urban electric distribution network
CN109934447A (en) A Fuzzy Comprehensive Evaluation Method for the Efficiency of Secondary Equipment in Smart Substations
CN110175789A (en) Quan Fangfa is assigned based on fuzzy hierarchy and the power distribution network diagnostic assessment index weights of CRITIC method
CN104376413A (en) Power grid planning scheme evaluation system based on analytic hierarchy process and data envelopment analysis
CN104462718A (en) Method for evaluating economic operation year range of transformer substation
CN110705887A (en) A comprehensive evaluation method of low-voltage station operating state based on neural network model
CN106485333A (en) A kind of transmission line of electricity running status appraisal procedure
CN114354783A (en) Health degree evaluation method of extra-high voltage oil chromatography monitoring device based on-operation data
CN104750976A (en) Establishment method of transmission line state evaluation parameter system
CN105225021A (en) The optimum choice method of power distribution network project yet to be built
CN109283293B (en) Power Transformer Fault Diagnosis Method Based on Coefficient of Variation and TOPSIS Method
CN106208043A (en) A kind of low and medium voltage distribution network key power consumption link discrimination method
CN110209996A (en) It is a kind of for evaluating the method and system of the general level of the health of protective relaying device
CN110930057A (en) Quantitative evaluation method of reliability of distribution transformer inspection results based on LOF algorithm
CN111667161B (en) Method for determining high-speed railway operation safety index
CN104616209B (en) A kind of power cable connector information fusion evaluation method based on on-line monitoring

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210105

RJ01 Rejection of invention patent application after publication