High-Voltage Cable Condition Assessment Method Based on Multi-Source Data Analysis

: In view of the problem that the weight value given by the previous state evaluation method is ﬁxed and single and cannot analyze the inﬂuence of the weight vector deviation on the evaluation result, a method based on the weight space Markov chain and Monte Carlo method (Markov chains Monte Carlo, MCMC) is proposed. The sampling method is used for evaluating the condition of high-voltage cables. The weight vector set obtained by MCMC sampling and the comprehensive degradation degree of the high-voltage cable sample are weighted and summed then compared in pairs to obtain the comprehensive degradation degree result. The status probability value and overall priority ranking probability of the object to be evaluated are obtained based on probability statistics, and the order of maintenance is determined according to the status probability value and the ranking result. It is realized that the cable line that needs to be identiﬁed in the follow-up defect is clariﬁed according to the evaluation result. This is helpful for operational and maintenance personnel to more accurately implement the maintenance plan for the cable and improve the operational and maintenance efﬁciency. Author Contributions: Conceptualization, X.-K.M. and Y.-B.J.; methodology, Z.-H.L. and Z.-Q.Y.; val-idation, P.-J.H.; formal analysis, Z.-M.L.; investigation, T.J.; resources, X.-K.M.; data curation, Z.-H.L.; writing—original draft preparation, Z.-Q.Y.; writing—review and editing, Y.-B.J.; funding acquisition,


Introduction
We have made great progress in high-level voltage and large-capacity power transmission with XLPE cables in recent years. High-voltage cables of 110 kV and above have been widely used in cities in China [1]. Factors such as the operating environment, working conditions and external forces of the cable may cause defects in operation. Nowadays, a variety of online monitoring systems for cable temperature, grounding circulation, partial discharge, etc. have been commonly installed on high-voltage cables of 110 kV and above [2]. However, the lack of integration of each cable monitoring system makes the monitoring data of various status variables relatively independent and unable to be comprehensively utilized.
Domestic and foreign research results on power cables are mainly focused on analyzing and judging the operating status of power cables, such as on-line monitoring of cable temperature [3], grounding circulation and partial discharge [4], on-line insulation measurement and infrared thermal imaging technology [5], etc. The methods are relatively simple for accurately judging the operating status of the cable, and they have not formed a comprehensive and systematic power cable operating status evaluation model. The research on the identification of high-voltage cable defects mainly focuses on the identification of partial discharge patterns (creeping discharge and corona discharge, etc.) [6], using the detection signal changes caused by cable defects (partial discharge signals) and directly using the detection signal threshold for fault warning and judgement [7]. There are many other types of defects in cables, except for partial discharge in actual situations, such as abnormal cable joint temperature and outer sheath damage [8]. Therefore, the cable-related index data obtained from the offline inspection of the cable and the online inspection can reflect the insulation condition of the cable to a certain extent. However, the research methods are relatively simple, and the multi-source information of the cable is not comprehensively considered, so it is impossible to establish a systematic cable condition evaluation model to comprehensively evaluate the cable condition.
Nowadays, scholars have begun to study comprehensive state evaluation methods of cables and have achieved some results. This has established a 10 kV cable operating state evaluation model based on the test results obtained by conducting macroscopic and microscopic tests on the test cable, and has fully considered the cable operating life, operating environment and other ancillary factors, and evaluated the new and old cables [9,10]. It applies D-S evidence theory to the state assessment of cables [11], uses multiple monitoring indicators, and proposes a comprehensive evaluation model for cable insulation that integrates evidence theory and fuzzy theory. It also established a comprehensive indicator system for cables [12], and obtained the weight values of various indicators based on the improved analytic hierarchy process [13]. On this basis, a comprehensive evaluation model of cable operating status with improved radar chart theory was established [14]. An analysis method combining fuzzy theory and analytic hierarchy process has been introduced to evaluate the cable status based on the particularity and difference of different influencing factors [15], and on the basis of the evaluation results, a maintenance strategy considering economy and reliability is formulated [16]. The insulation resistance, dielectric loss angle and partial discharge of the disabled cables and their accessories were comprehensively measured. A cable condition assessment model was established based on the measurement results and the national standard GB12706 [17]. Compared with the cable body, the terminals and intermediate joints of the cable accessories are more prone to failure due to their complicated design and production defects [18,19]. The aforementioned cable state evaluation methods all give fixed weight values, but do not consider the difference and uncertainty of the group evaluation weight vector. It is impossible to analyze the influence of the weight value on the evaluation result when the weight value changes within a given interval, nor can it give the probability value of each state of the object to be evaluated. A thermal state evaluation method based on the maximum temperature difference on the terminal surface is proposed by analyzing the thermal field distribution characteristics of the test cable terminal surface [20]. The multi-state index evaluation of high-voltage cables is a multi-attribute comprehensive evaluation problem, and the weight of each index has a great impact on the evaluation result [21]. The advantages of Markov chain Monte Carlo (MCMC) in the evaluation of multi-state indicators are followed, which can solve the problem of the influence of the weights of different indicators on the evaluation result. MCMC requires to evaluate the full data likelihood at different parameter values iteratively [22]. MCMC users often decide convergence by applying empirical diagnostic tools in practice [23]. It is applied to multi-source data analysis of cable PD ultrasonic signals based on the engineering value of MCMC [24]. The high-voltage cables work in a very strong electric field, and the existing online and offline detection technologies, test methods, condition evaluation methods, and defect identification methods are not complete. How to comprehensively utilize and quantify the multi-source data of cables as the basis for evaluating the status of high-voltage cables is a difficult problem for improving the reliability and safety of the power supply operation system.
In order to solve the problem that the weight value of the high-voltage cable status evaluation is fixed and single and cannot analyze the influence of the weight deviation on the evaluation result; the cable status evaluation system is established based on the weight space of MCMC sampling, and the multi-source status indicator system of the cable is built. The proposed method was verified by an example based on the relevant data provided by a certain electric power bureau.

The Weight Space of Multi-Attribute Comprehensive Evaluation
In order to fully dig the multi-source data such as online monitoring and offline testing of high-voltage cables, comprehensive evaluation of the operating status of high-voltage cables is carried out in consideration of the difference and uncertainty of the weight value. In actual situations, the weight samples determined by a small number of expert groups are often limited, and the comprehensive evaluation results obtained with only a few weight vectors are not necessarily accurate, and it is impossible to analyze the influence of the weight vector deviation on the comprehensive evaluation results. According to the limited weight vector data, the random distribution that the weight vector of the expert group may satisfy is estimated, and the amount of available weight sample data is artificially expanded, which can make the evaluation result more reasonable.

Definition of the Weight Space
The determination of the weight value is closely related to the expert's work experience, knowledge structure system and the degree of grasp of the evaluated object. Scientific and reasonable processing of the evaluation results of the expert group on the weight is an important link in the reasonable determination of the weight of each attribute.
Assuming that the n-dimensional of the weight vectors given by a finite group of different experts obey a certain distribution, the relevant parameter quantities of the distribution are determined according to the weight sample information that has been obtained, and then random weight samples that meet the requirements are drawn from the known weight distribution through the sampling method. In this way, the required weight sample data amount is obtained. The concept of the weight space is introduced. The weight space can be defined as a set of all possible weight combinations, which is expressed as: There, the value range of the weight is 0-100%, and the sum of all weight is 100%.

Weight Space Constraints
More weight vector data can be obtained for comprehensive evaluation by sampling the high-dimensional weight space of multiple attributes, and more reasonable comprehensive evaluation results can be obtained. Because the importance of each attribute is different from each other, two points need to be determined when sampling the weight space: (1) the constraint condition that the weight of each attribute satisfies; (2) the probability that the weight vector given by different weighting methods or expert groups is satisfactorily distributed. The constraints on the corresponding weight of each attribute of the weight vector include the weight interval of a single attribute, the relative relationship between the weight of two attributes, and so on. The constraints on the corresponding weight of each attribute of the weight vector include the weight interval of a single attribute, the relative relationship between the weight of two attributes, and so on. If the i-th attribute interval of a given weight vector W k is w ik [w ikmin , w ikmax ], then the sum of all attribute weights is 1, and the interval constraint can be converted to, According to the above formula, the interval constraint of a single attribute can be transformed into a relative relation constraint with the sum of the weights of the remaining attributes. The relative relationship constraint between the i-th and j-th attribute weights of the weight vector W k is expressed as: where, η w ij is a constrained weight space of order n. Taking the third-order weight vector as an example(W K = [w 1k ,w 2k ,w 3k ] T ), the three-dimensional weight space constrained by two conditions in Figure 1 (w 1k ∈ [0.1, 0.8],w 1k /w 3k ≤ 3). The three-dimensional weight space is represented by the equilateral triangles. The area of the upper triangle of the parallel lines L 1 and L 2 represents the size of the weight w 1k . The area of the lower trapezoid represents the sum of weights w 2k and w 3k . The middle trapezoidal region surrounded by L 1 and L 2 is the feasible region of the weight space determined by the w 1k interval constraint. One end of L 3 is vertex w 2k , and the ratio of the area of the upper and lower triangles is w 1k /w 3k . The feasible region formed by L 1 , L 2 , and L 3 is the constrained weight space.
[ ] 1 , , 1 j n j ≠ ∈ According to the above formula, the interval constraint of a single attribute can be transformed into a relative relation constraint with the sum of the weights of the remaining attributes.
The relative relationship constraint between the i-th and j-th attribute weights of the weight vector Wk is expressed as: where, ij w η is a constrained weight space of order n. Taking the third-order weight vector as an example(WK = [w1k,w2k,w3k] T ), the three-dimensional weight space constrained by two conditions in Figure 1 (w1k ∈ [0.1,0.8],w1k/w3k ≤ 3). The three-dimensional weight space is represented by the equilateral triangles. The area of the upper triangle of the parallel lines L1 and L2 represents the size of the weight w1k. The area of the lower trapezoid represents the sum of weights w2k and w3k. The middle trapezoidal region surrounded by L1 and L2 is the feasible region of the weight space determined by the w1k interval constraint. One end of L3 is vertex w2k, and the ratio of the area of the upper and lower triangles is w1k/w3k. The feasible region formed by L1, L2, and L3 is the constrained weight space. We can extract the weight vector of the known probability distribution from the constrained weight space, and perform multi-attribute comprehensive evaluation on the object to be evaluated, based on extracting a given number of weighted sample sets, and finally obtain the status result of each object to be evaluated.

MCMC Sampling
(1) Markov Chain Markov Chain (Markov Chain, MC) is a discrete-time random process that satisfies Markov properties in mathematics. According to the transition probability between each state, we can transit from the current state to the rest of the state or transition to its own state. The mathematical definition of a Markov chain is [22]: There, {xt, t ∈ T} is set to a random process. For any positive integer n and t1 < t2 < … < tn, there will be: We can extract the weight vector of the known probability distribution from the constrained weight space, and perform multi-attribute comprehensive evaluation on the object to be evaluated, based on extracting a given number of weighted sample sets, and finally obtain the status result of each object to be evaluated.

MCMC Sampling
(1) Markov Chain Markov Chain (Markov Chain, MC) is a discrete-time random process that satisfies Markov properties in mathematics. According to the transition probability between each state, we can transit from the current state to the rest of the state or transition to its own state. The mathematical definition of a Markov chain is [22]: There, {x t , t ∈ T} is set to a random process. For any positive integer n and t 1 < t 2 < . . . < t n , there will be: and the conditional distribution is: {x t , t ∈ T} is called the Markov process. In the random process, the conditional probability can satisfy the following formula within any integer n ∈ T and any i 0 , i 1 , . . . , i n+1 ∈ I.
{x n , n ∈ T} is called the Markov chain. In the Markov process, the matrix formed by the transition probabilities of each state is called the transition probability matrix, which is denoted as P = [P ij ] nxn , where P ij represents the probability that the system is in state i at time t and state j at time t + 1. The state transition matrix P and the probability distribution π(x) satisfy the following expressions for all i and j: where, π(x) is the stationary distribution of the state transition matrix P.
(2) Monte Carlo The Monte Carlo method (MC) is a numerical method based on random experiments and statistical calculations. The problem to be solved is abstracted as a mathematical problem and solved by establishing a probability model. The principle is to take the mathematical expectation of random variables or the probability of occurrence of random events as the solution to the problem to be sought.
First of all, a probabilistic process was constructed artificially, and some of its parameters are just the solutions to the problems to be sought. Essentially, the process is to transform non-random problems into random problems. Secondly, the constructed various probability models can be considered to be composed of different probability distributions, so simulation experiments can be realized by generating random vectors (variables) with known probability distributions. Finally, when the simulation experiment is completed, the solution of the problem is obtained based on the construction of various relevant estimators for the simulation experiment results that need to be investigated and registered.
(3) MCMC sampling The Markov chains with Monte Carlo (Markov chains Monte Carlo, MCMC) sampling methods were combined to sample the high-dimensional weight space, and enough weighting schemes can be obtained. Gibbs sampling in MCMC is used to obtain a series of weight samples approximately equal to a given multi-dimensional probability distribution.
The weight space sampling based on MCMC mainly includes three steps: (1) A suitable Markov chain is constructed in the weight space W, which can converge to a stable distribution π(x); (2) A sequence of sample points x 1 , x 2 , x 3 , . . . x N is generated based on the Markov chain in step 1), which starts from a certain point x 0 in W; (3) When the Markov chain is iterated by M samples that starts from x 0 , the marginal probability distribution of x M , at each moment is a stable distribution π(x) and reaches convergence.
The distribution that the weight vector satisfies is unknown in the actual situation. It is necessary to estimate the random distribution that the weight vector of the expert group may satisfy based on the limited weight vector data. The n-dimensional weight vectors given by different expert groups obey the n-dimensional normal distribution supposed, which is W~N n (µ, ∑), where µ is the expectation vector and ∑ is the covariance matrix. When q elements in W are unknown (q < n) and other elements are known, there W, µ, and ∑ can be divided into block matrices.
where, W (1) and µ (1) are the weight and expectation vector of q × 1 order. ∑ 11 is the q × q order covariance sub-matrix. The posterior conditional distribution of each element of different n-dimensional weight vectors at the given time of other elements can be obtained when the different elements form W (1) is chose and q = 1. The conditional distribution of weight vector elements w 1 and w 2 can be obtained at n = 2.
The MCMC sampling output based on the software Matlab is shown in Figure 2. The Markov chain in MCMC sampling moves forward in the direction of x 1 , and then moves in the direction of x 2 to perform alternate sampling of the coordinate axis, which can be seen from Figure 2 in each iteration.  (8) where, W (1) and μ (1) are the weight and expectation vector of q × 1 order. ∑11 is the q × q order covariance sub-matrix. The posterior conditional distribution of each element of different n-dimensional weight vectors at the given time of other elements can be obtained when the different elements form W (1) is chose and q = 1. The conditional distribution of weight vector elements w1 and w2 can be obtained at n = 2.
The MCMC sampling output based on the software Matlab is shown in Figure 2. The Markov chain in MCMC sampling moves forward in the direction of x1, and then moves in the direction of x2 to perform alternate sampling of the coordinate axis, which can be seen from Figure 2 in each iteration.

Comprehensive Evaluation Process
The multi-attribute comprehensive evaluation problem contains m objects Si (i = 1, 2, …, n) to be evaluated, and each object has n attributes. It can be expressed as Si = [xi1, xi2, …, xin]. Each object attribute should be standardized or degraded as F(Si) = [f(xi1), f(xi2),…, f(xin)] in multi-attribute decision-making to obtain a standardized data matrix F(S) as shown in Table 1. The comprehensive degree of degradation is used as a measure of the final state of the cable, since the degree of degradation can characterize the health of power equipment. The comprehensive evaluation results of cables are divided into four intervals: (0,0.25), (0.25,0.5), (0.5,0.75), (0.75,1), corresponding to the normal Q/GDW11262-2014 power cable and channel maintenance procedures status, attention status, abnormal status, and severe status.

Comprehensive Evaluation Process
The multi-attribute comprehensive evaluation problem contains m objects S i (i = 1, 2, . . . , n) to be evaluated, and each object has n attributes. It can be expressed as . . , f (x in )] in multi-attribute decision-making to obtain a standardized data matrix F(S) as shown in Table 1. The comprehensive degree of degradation is used as a measure of the final state of the cable, since the degree of degradation can characterize the health of power equipment. The comprehensive evaluation results of cables are divided into four intervals: (0, 0.25), (0.25, 0.5), (0.5, 0.75), (0.75, 1), corresponding to the normal Q/GDW11262-2014 power cable and channel maintenance procedures status, attention status, abnormal status, and severe status.
The weights of each attribute given by the expert group or the weights obtained by the weight space sampling method will be inconsistent in practical applications. There are l types of empowerment schemes is supposed that expressed as W k = [w 1k , w 2k , . . . , w nk ] T , Energies 2022, 15, 1369 7 of 16 k = 1, 2, . . . , l. The comprehensive evaluation result vector of each evaluated object under different weights can be described as: The probability matrix of S i belonging to different states can be obtained based on the 1 × l-order comprehensive evaluation result vector V(S i , W) of the evaluation object S i under l-weighting scheme, which can be expressed as follows: where, a m1 , a m2 , a m3 , a m4 are representing channel maintenance procedures status, attention status, abnormal status, and severe status. For any two evaluated objects S i1 and S i2 , the probability that the evaluation result of S i1 is better than S i2 under a given weight vector distribution is: where, P(S i1 , S i2 ) is called the priority ranking probability of S i1 to S i2 under the set of weight vectors, and W k is the random variable in the weight space. V(S i , W) is the 1 × l-order comprehensive evaluation result vector of the evaluated object S i under l-weighting scheme. F(S i ) is the normalized data vector of 1 × n-order of the object S i . W k is the n × 1-order weight vector of the k-th weighting scheme.
Pairs of m evaluated objects can be combined to get m 2 -m combinations. Each combination can be performed MCMC simulation by using constraint weight space. Finally, the priority probability P between the two objects is calculated. MCMC simulation is performed on the combination of m 2 -m evaluated objects, and the priority ranking probability matrix of m evaluated objects can be described as: where, P is the m × m-order priority ranking probability matrix, and the diagonal elements are P mm = 0.5. The off-diagonal elements are P ij = P(S i ≥ S j ), and there exist P ij + P ji =1 theoretically., there is |P ij + P ji − 1| ≤ α(the simulation error α = 0.001), and the precise values of P ij and P ji (i, j = 1, 2, . . . , m) can be simulated independently, which shows that the credibility of the simulation results can be supported by the number of simulations, otherwise the number of simulations l needs to be increased. According to the non-diagonal elements of each row of the priority ranking probability matrix P, the overall priority ranking probability H i of the i-th evaluated object can be described as: The i represents the average priority ranking probability of the evaluation object S i relative to the subsequent evaluation object in the i-th place. In summary, the specific steps for comprehensive evaluation of high-voltage cable status are as follows: (1) A single abnormality judgment is performed on the index with high priority based on the collected cable sample data. These cables that show obvious abnormal data in a single index are directly included in the scope of maintenance and defect identification. The remaining cable samples are comprehensively evaluated. (2) The degree of degradation of the cable data to be evaluated is obtained based on the above formula. (3) There l-weight vector samples (given by the relevant decision maker) of the highdimensional weight space are inferred by MCMC sampling. The probability matrix A and priority ranking probability matrix P of each evaluation object belonging to different states are obtained after weighted summation of the sample degradation degree and pairwise comparison. (4) The overall priority ranking probability H i and the credibility coefficient σ i of the i-th evaluated object are obtained based on the priority ranking probability matrix. Finally, the inspection scope and the priority of the inspection are determined according to the sorting results.

Comprehensive State Evaluation of High-Voltage Cables
The data sources for comprehensive state evaluation and defect identification of highvoltage cables can be divided into five aspects: online monitoring data, offline test data, operation and maintenance information, equipment information and environmental and external force damage data, which is shown in Figure 3. According to the actual situation of 220 kV high-voltage cables in a certain power supply bureau, this paper verifies the highvoltage cable condition evaluation method proposed based on the multi-source information of high-voltage cables.
The partial discharge ultrasonic sensor and grounding circulation sensor are mainly installed on the grounding wire of each terminal and intermediate joint of the cable, and the cable temperature data is measured through distributed optical fibers. There are differences of monitoring and test data in the sampling period, data structure, measurement location, etc. For example, each cable usually has multiple cable terminals and intermediate joints with the grounding circulation and partial discharge data of cables. Multiple ground circulation and partial discharge monitoring data are received at every acquisition time. The monitoring data of the entire line length can be obtained at each acquisition moment for the temperature and current carrying capacity of high-voltage cables. It constitutes a data set of the cable real-time data based on online monitoring and other offline data.  The partial discharge ultrasonic sensor and grounding circulation sensor are mainly installed on the grounding wire of each terminal and intermediate joint of the cable, and the cable temperature data is measured through distributed optical fibers. There are differences of monitoring and test data in the sampling period, data structure, measurement location, etc. For example, each cable usually has multiple cable terminals and intermediate joints with the grounding circulation and partial discharge data of cables. Multiple ground circulation and partial discharge monitoring data are received at every acquisition time. The monitoring data of the entire line length can be obtained at each acquisition moment for the temperature and current carrying capacity of high-voltage cables. It constitutes a data set of the cable real-time data based on online monitoring and other offline data.

Selection of Comprehensive Evaluation Indexes for High-Voltage Cable State
The health of the cables will gradually change during the operation of high-voltage cables, due to factors such as operating time and working environment. In order to ensure the comprehensive evaluation result of the high-voltage cable status is accurate, the factors that have a certain impact on the cable operating status should be selected as comprehensively as possible. The comprehensive evaluation indicators of cable status selected in this paper are shown in Table 2.

Selection of Comprehensive Evaluation Indexes for High-Voltage Cable State
The health of the cables will gradually change during the operation of high-voltage cables, due to factors such as operating time and working environment. In order to ensure the comprehensive evaluation result of the high-voltage cable status is accurate, the factors that have a certain impact on the cable operating status should be selected as comprehensively as possible. The comprehensive evaluation indicators of cable status selected in this paper are shown in Table 2. When the degradation degree is set to 0, it indicates that the device is in a healthy state. when the degradation degree is set to a value of 1, it indicates that the device is in a fault state. The final state of the cable is measured by the comprehensive degradation degree in this paper. It is generally hoped that the number of failures or overhauls is as small as possible for discrete data such as the number of failures and the number of overhauls. The calculation formula for defining their degradation degree is as follows: where, X and X min are the historical faults or overhaul times since the cable was put into operation and the corresponding limit of times.

Case Analysis
There are 20 circuits of 220 kV cables in the jurisdiction of a power supply bureau, of which 12 circuits are equipped with optical fiber temperature monitoring devices, and 16 are equipped with cable outer sheath circulation monitoring devices; there are also 7 cables with partial discharge monitoring devices installed. A total of 15 cables with five cycles of monitoring devices for temperature, partial discharge, and grounding circulation are installed at the same time as examples to test the effectiveness of the comprehensive evaluation method in this paper. The basic information of each evaluated cable is shown in Table 3. The cable model is 220 kV-XLPE-2500 mm copper core single core aluminum sheath. The worst value of cable online monitoring, offline test, operation and maintenance and environmental historical data before 20 July 2020 was chosen to comprehensively evaluate the operating status of each cable, assuming that all collected data correspond to the true status of the cable. Table 3. Basic information of cables.

Length
(m) Firstly, whether each cable has a single indicator abnormal state need to be judged. It can be seen that the conductor temperature of Phase A near the GIS terminal of the substation 750 m exceeds the maximum allowable temperature of 130 • C based on the temperature data of Line 1 in Figure 4, and it is judged as S 1A abnormal. According to the ratio of the circulating current to the load current at each terminal and intermediate joint of Line 2 in Figure 5, it can be seen that the circulating current at the B-phase directly grounded joint #9 exceeds 20% of the load current, and it is judged that the cable S 2B is abnormal.
The insulation resistance value of the outer sheath of the abnormal section of Line 3 and Line 4 in Table 4. It can be seen from the table that the insulation resistance of the outer sheath of each phase in the #4-#5 sections of Line 3 does not meet the requirement of less than 0.5 MΩ/km. The insulation resistance of the outer sheath of the #9-#10 sections and #12-#13 sections of Line 4 are also unqualified, so it is determined that the cables S 3A , S 3B , S 3C and S 4A , S 4B , S 4C are abnormal.  Where, #0 represents the GIS terminal of the substation. There #1,2,4,5,7,8,10,11 represent the cross interconnection connector of cable. There #3,6,9,12 represent the directly grounded connector of cable.
The insulation resistance value of the outer sheath of the abnormal section of Line 3 and Line 4 in Table 4. It can be seen from the table that the insulation resistance of the outer sheath of each phase in the #4-#5 sections of Line 3 does not meet the requirement of less than 0.5 MΩ/km. The insulation resistance of the outer sheath of the #9-#10 sections and #12-#13 sections of Line 4 are also unqualified, so it is determined that the cables S3A, S3B, S3C and S4A, S4B, S4C are abnormal.
It can be determined that the eight cable lines S1A, S2B, S3 and S4 all have abnormalities in different situations, which is judged by a single index, and corresponding maintenance and repair work need to be prioritized. The abnormal cables are eliminated, and the validity of the evaluation method in this paper is tested based on the remaining seven cables. The test result is compared with the imported reference sample, which corresponds to the best value of each status indicator in the historical data set.    Where, #0 represents the GIS terminal of the substation. There #1,2,4,5,7,8,10,11 represent the cross interconnection connector of cable. There #3,6,9,12 represent the directly grounded connector of cable.
The insulation resistance value of the outer sheath of the abnormal section of Line 3 and Line 4 in Table 4. It can be seen from the table that the insulation resistance of the outer sheath of each phase in the #4-#5 sections of Line 3 does not meet the requirement of less than 0.5 MΩ/km. The insulation resistance of the outer sheath of the #9-#10 sections and #12-#13 sections of Line 4 are also unqualified, so it is determined that the cables S3A, S3B, S3C and S4A, S4B, S4C are abnormal.
It can be determined that the eight cable lines S1A, S2B, S3 and S4 all have abnormalities in different situations, which is judged by a single index, and corresponding maintenance and repair work need to be prioritized. The abnormal cables are eliminated, and the validity of the evaluation method in this paper is tested based on the remaining seven cables. The test result is compared with the imported reference sample, which corresponds to the best value of each status indicator in the historical data set.   It can be determined that the eight cable lines S 1A , S 2B , S 3 and S 4 all have abnormalities in different situations, which is judged by a single index, and corresponding maintenance and repair work need to be prioritized. The abnormal cables are eliminated, and the validity of the evaluation method in this paper is tested based on the remaining seven cables. The test result is compared with the imported reference sample, which corresponds to the best value of each status indicator in the historical data set.
Actually, the distribution that the weight vector satisfies is not known, so it is necessary to estimate the random distribution that the weight vector of the expert group may satisfy based on the limited weight vector data. According to the estimated probability distribution, by constructing the transition conditional probability, enough weight sample data are drawn upon. The transition conditional probability is constructed according to the estimated probability distribution, and extract enough weighted sample data. The amount of available weight sample data is artificially expanded based on this method. Since the uncertainty and difference of the weights are considered, the evaluation results are more reasonable. This is assuming that the weight vectors given by multiple experts satisfy the n-dimensional normal distribution.
The importance and reliability of cable data from different sources are different. The data obtained through offline testing are the most accurate, and the monitoring data may have errors in the monitored data due to sensor or communication abnormalities, so its accuracy is lower than offline testing. In practice, cable failures occur from time to time due to external force damage such as construction, but the probability of occurrence is lower than the accuracy of offline testing and online monitoring data. The environment has the weakest influence on the evaluation of cable status. Therefore, the reliability of the cable status data in Table 2 from high to low is offline test data, online monitoring data, external force damage data, operation and maintenance data, and environmental data. The relative relationship of the weight vectors of various indicators can be obtained as: 18 , w 19 , w 20 , w 21 } ≥ {w 9 , w 10 , w 11 , w 12 , w 13 , w 14 } ≥ {w 15 , w 16 , w 17 } (17) where, w 5 > (w 1 , w 2 , w 3 , w 4 ) is equivalent to w 5 /w 1 > 1, w 5 /w 2 > 1, w 5 /w 3 > 1 and w 5 /w 4 > 1.
By analogy, 74 conditional constraints of the weight vector can be obtained. Therefore, the 500 weight vector obtained by MCMC sampling for the constrained weight space is shown in Figure 6.
satisfy the n-dimensional normal distribution. The importance and reliability of cable data from different sources are different. The data obtained through offline testing are the most accurate, and the monitoring data may have errors in the monitored data due to sensor or communication abnormalities, so its accuracy is lower than offline testing. In practice, cable failures occur from time to time due to external force damage such as construction, but the probability of occurrence is lower than the accuracy of offline testing and online monitoring data. The environment has the weakest influence on the evaluation of cable status. Therefore, the reliability of the cable status data in Table 2 from high to low is offline test data, online monitoring data, external force damage data, operation and maintenance data, and environmental data. The relative relationship of the weight vectors of various indicators can be obtained as: The sampling results of the first 100 sets of weight vectors are shown in Figure 5, where the maximum weight is distributed in the interval w5~w8, and the maximum value does not exceed 0.2. The weights of cable operation and maintenance and environmental data w9~w14 and w15~w17 are relatively small. The operating status of the cable is mainly reflected by offline testing, online monitoring and external force damage data, but the weights of different status indicators show random distribution characteristics under given constraints. According to the weight vector sampling results in Figure 5 and the The sampling results of the first 100 sets of weight vectors are shown in Figure 5, where the maximum weight is distributed in the interval w 5 -w 8 , and the maximum value does not exceed 0.2. The weights of cable operation and maintenance and environmental data w 9 -w 14 and w 15 -w 17 are relatively small. The operating status of the cable is mainly reflected by offline testing, online monitoring and external force damage data, but the weights of different status indicators show random distribution characteristics under given constraints. According to the weight vector sampling results in Figure 5 and the evaluated cable data of the reference sample, the evaluation results and their probability distributions of seven evaluated cables can be obtained, as shown in Figures 7 and 8.
The probability of the cable S 1B is 1, which has priority over S 1C , S 2A , S 2C , and S 5C , and the probability of priority over S 5A (0.994) and S 5B (0.812), there the cable S 1B is in the best state of all the cables evaluated. It can be seen from Figure 6 that the evaluation results of S 1B are relatively smallest under most of the weight vectors, but in the 500 groups of weight vectors, there are still 94 groups and 3 groups of weight vectors. The evaluation results of S 1B are worse than S 5A and S 5B . According to Formula (12), the overall priority ranking probability vector can be calculated as H = [0.8295, 0.44, 0.0065, 0.1396, 0.5670, 0.737, 0.2915] T , from which the order of the evaluated cable status is S 1B > S 5B > S 5A > S 1C > S 5C > S 2C > S 2A . According to Formula (15), the credibility coefficient of this ranking is calculated as θ = 0.981. ping area; the evaluation results of S1C and S2A have a large distribution interval, which shows that the dispersion of the weight vector has a great influence on the evaluation of the state of S1C and S2A. When MCMC is not used to sample the weight space, the weight vector obtained is a fixed value. If the weight vectors obtained at this time are W(104) and W(413) as shown in Figure 7, the sorting results of each cable status are respectively S5B > S1B > S5A > S5C > S1C > S2C > S2A and S5B > S5A > S1B > S5C > S1C > S2C > S2A, and the evaluation results of the corresponding cables are also marked in Figure 7.  The method in this paper performs MCMC sampling on the weight space, which can obtain the overall probability distribution of the relative relationship of the evaluated objects from the perspective of the randomness of the group weighting; compared with the existing weighting methods, the method in this paper is helpful for analyzing the weights in a given interval. The impact of internal changes on the evaluation results is According to the probability distribution of the status evaluation results of each cable in Figure 7, it can be seen that the relative relationship between the expected values of the evaluation results of different cables is consistent with the aforementioned status ranking results. The evaluation results of S 1B , S 5B , and S 5A are distributed in the overlapping area; the evaluation results of S 1C and S 2A have a large distribution interval, which shows that the dispersion of the weight vector has a great influence on the evaluation of the state of S 1C and S 2A . When MCMC is not used to sample the weight space, the weight vector obtained is a fixed value. If the weight vectors obtained at this time are W(104) and W(413) as shown in Figure 7, the sorting results of each cable status are respectively S 5B > S 1B > S 5A > S 5C > S 1C > S 2C > S 2A and S 5B > S 5A > S 1B > S 5C > S 1C > S 2C > S 2A , and the evaluation results of the corresponding cables are also marked in Figure 7.
The method in this paper performs MCMC sampling on the weight space, which can obtain the overall probability distribution of the relative relationship of the evaluated objects from the perspective of the randomness of the group weighting; compared with the existing weighting methods, the method in this paper is helpful for analyzing the weights in a given interval. The impact of internal changes on the evaluation results is shown.
According to Formula (11), the probability matrix of each evaluation cable belonging to different states is:  Figure 7 and the calculation results of the probability matrix, it can be seen that S 1B , S 5A , and S 5B are in the normal state, and S 1C and S 5C are in the normal state with a probability of 97.7% and 85.8%, respectively. S 2A is in the interval of the attention state, and S 2C belongs to the attention state with a probability of 97.5%. It can be seen that the temperature values of S 2A and S 2C exceed the allowable long-term operation temperature of 90 • C but do not exceed the maximum allowable temperature of 130 • C based on the cable data, and the amount of partial discharge is close to the extreme value of the normal aging standard. In practice, cable failure events rarely occur, and it can be seen that the result is consistent with the actual situation. In actual situations, cable failure events occur rarely, and it can be seen that the result is consistent with the actual situation. For this reason, it can be determined that S 2A and S 2C cables need to be prioritized for routine testing and periodic maintenance.

Conclusions
The weight vector set obtained by MCMC sampling proposed in this paper takes into account the difference and uncertainty of the weight, and can reflect the influence of the weight value on the evaluation result of the high-voltage cable state when the weight value changes within a given interval. The overall priority probability of each cable operating state can be obtained by constructing the priority sequence probability matrix of the evaluated cables, which can guide the operation and maintenance personnel to formulate the sequence of cable maintenance. The superiority of this method is even more obvious when hundreds or thousands of cables are being overhauled at the same time. The scope of maintenance can be shortened, and the cable lines to be inspected can be clarified, which greatly saves manpower and financial resources. Funding: This work was supported in part by State Grid Shanxi Electric Power Corporation headquarter's science and technology project (52053020000V), the Hebei University "One Province One University" Special Fund (801260201232)", and the president fund project of Hebei University (XZJJ201908), and postdoctoral programme of State Grid Shanxi Electric Power Corporation and State Grid Shanxi Electric Power Corporation headquarter's.