Random Violation Risk Degree Based Service Channel Routing Mechanism in Smart Grid

Smart gird, integrated power network with communication network, has brought an innovation of traditional power for future green energy. Optical fiber technology and synchronous digital hierarchy (SDH) technology is widely used in smart grid communication transmission network. It is a challenge to reduce impact of the availability of smart grid communication services caused by random failures and random time to repair. Firstly, we create a service channel violation risk degree (SCVRD) model to precisely track the violation risk change of communication service channel. It is denoted by the probability of service channel cumulative failure duration exceeding the prescribed duration. Secondly, a service channel violation risk degree routing mechanism is proposed to improve the availability of communication service. At last, the simulation is implemented with MATLAB and network data in one province are used as data instance. The simulation results show that the average service channel failure rate of availability-aware routing based on statistics (AAR-OS) algorithm and risk-aware provisioning algorithm are reduced by 15% and 6%, respectively.


Introduction
Smart gird, integrated power network with communication network, incorporates the latest innovative technologies to bring a revolutionary change and innovation of traditional power for future green energy [1][2][3][4][5].Although smart grid has lots of promising features, such as intelligent de-centralized control, resilience, flexibility, sustainability, digitalization, intelligence, consumer empowerment, renewable energy, smart infrastructure and so on, a number of critical challenges and open issues like need to be further discussed [6][7][8][9].One of these critical challenges is risk of smart grid communication which is always the critical constraint in ultra high voltage (UHV) Grid [7], Smart Home [8], Microgrids [9] and other smart grid applications.
Due to the complexity of smart grid communication environment, different communication technologies are used for the realization of smart grid, such as Optical fiber technology, power line communication (PLC), 4G/5G, wireless mesh network and so on [10][11][12].Among these communication technologies, optical fiber technology and the corresponding synchronous digital hierarchy (SDH) technology is widely used in the smart grid communication transmission network for the high bandwidth, high anti-interference, small signal attenuation and long transmission distance.However, there still exists communication violation risk.
With SDH technology, one optical fiber can carry multiple service channels.Each single service channel may have a great transmission capacity, such as STM-16 (10 Gbit/s).In this case, even a short interruption of the fiber can still cause a large amount of data loss.Therefore, the occurrence of interrupted service channel must be limited.To achieve this purpose, the statistical path availability A SPA model is normally adopted to model the service channel and the basic rules to guarantee the availability of service channels are also specified [12].According to the rules, in the process of service channel planning, different levels of electric power communication services should be allocated to the different channels.For example, to guarantee the service quality, a service with statistical availability greater than 99.9% should be allocated to a channel with the predetermined statistical availability higher than 99.9%.In addition, to ensure the requirements of high availability and high real-time, services are planned with the primary channel and backup channel.Thus, when some electric power communication failure events occur and some channels are interrupted, the services can be safeguarded.
However, in practice, the faults of transmission equipment and optical cables carrying the service channels occur randomly.Therefore, during a period of time, the actual path availability A APA (t) of service channels may be significantly higher than 99.9%, or may be less than 99.9% due to a sudden failure.The backup channel strategy cannot ensure that the actual channel availability will not violate the rules.Thus, in practice, there exists a risk that service channel may violate the availability requirements, namely there exists a service channel violation risk (SCVR).In this case, the challenge that electric power communication network channel planning faces is how to effectively control the availability decline because of the channel random failures, thus to reduce the probability of SCVR.
To solve this problem, the availability-aware routing mechanism should be considered.Through scientific and rational route planning, when a channel fails, the service carried by the channel can be conveyed by another channel and will not be affected, thus the violation risk caused by channel random failure can be avoided.To achieve the above goal, firstly, a probability distribution model of SCVR should be studied and designed.Then, a SCVR based routing mechanism should be proposed to reduce the failure number (FN) and failure duration of the electric power communication service.The goal is to minimize the violation risk caused by the random failures of transmission equipment and optical cable and thus to improve the availability of electric power communication service.
In this paper, the main contributions include: (1) A probability distribution model of SCVR, named service channel violation risk degree (SCVRD) model, is proposed, which is denoted by the probability of service channel cumulative failure duration exceeding the prescribed duration.(2) Based on SCVRD, a service channel violation risk degree routing (SCVRD-R) algorithm is proposed to improve the availability of electric power communication service.
The remainder of the paper is organized as follows.Section 2 reviews the related work and analyzes the limitation of the current works.In Section 3, the differences between A SPA and A APA (t) are analyzed and the SCVRD model is proposed.Section 4 gives the approximate transformation of violation risk distribution and proposes the SCVRD-R algorithm.Section 5 discusses the simulation results.Finally, Section 6 concludes the paper.

Related Work
Currently, the body of work related to smart grid communication robustness is rapidly increasing.For the realization of smart grid, many studies have looked at the communication challenges in smart grid.
Papers [6][7][8][9] pointed the critical communication challenge of smart grid communications and also gave the feasible solutions and future directions from the overall perspective.One way is to improve the communication reliability by eliminating the defects of the technology itself.Paper [10] proposed an orthogonal poly-phase-based multicarrier code division multiple access (OPP-MC-CDMA) system and implemented with a minimum mean square error equalizer and nonlinear preprocessing to overcome the effects of noise and multipath frequency-selective fading commonly experienced in PLC channels.This way is the most effective but it depends on the update of corresponding communication technology and it is difficult to make great progress.
The other way is to optimize the communication routing of communication service.Some studies adopt the method of optimizing the routing protocol of the corresponding communication network technology.Paper [11] presented a QoS-aware wireless mesh network (WMN) routing technique that employed multiple metrics in optimized link state routing (OLSR) for AMI applications in a smart grid neighbor area network based wireless mesh network.They indicate to guarantee the optimized communication routing.Other studies turn to solve the service channel failure risk problem to guarantee the effective communication routing.The main way is to control the SCVR by routing method.There are two kinds of routing methods: availability-aware routing based on statistics (AAR-OS) and availability-aware routing based on uncertainty (AAR-OU).
In respect of AAR-OS, the difference between A SPA and A APA (t) and how availability changes over time and geographical locations are pointed out in Reference [13].Based on the new availability calculation method, the 3W-availability aware routing (3WAR) algorithm was proposed in that paper which effectively narrowed the gap between the actual availability and target availability.In paper [14], the definitions of min cross layer cut (MCLC) and min cross layer spanning tree (MCLST) were given and the availability routing algorithm under different failure probability conditions was proposed to maximize the MCLC and minimize the MCLST.Paper [15] adopted the log information of path state as the basis for routing and considered the path with highest statistics availability as the service channel.Papers [16][17][18] proposed A SPA based multipath routing mechanism by increasing the redundancy of resources to enhance the A SPA of the channel.Paper [19] proposed a primary-backup sharing routing mechanism in optical networks to improve the resource utilization in premise of ensuring the statistical availability.In paper [20], the cost of routing was taken into account in routing algorithm to minimize the cost in premise of ensuring the statistics availability.
The methods mentioned above are all using the A SPA as the decision indicator in routing mechanism.The advantages of those methods are simple, easy to reflect the availability in the overall trend and clear in physical meaning.But A SPA only has statistical meaning, which can reflect the availability variation trend on the whole but cannot reflect the actual availability fluctuations of channel.Thus, the threat to the grid due to network random failures cannot be effectively reduced by those methods.
In the aspect of AAR-OU, papers [21][22][23][24][25] all studied the uncertainty of the path availability during a short time period.Paper [21] proposed a dynamic availability-aware survivable routing architecture to provide the service path protection based on the partial restorability.Papers [22,23] defined the concept of availability border and proposed the path availability evaluation method and the routing algorithm on the assumption that failure arrival rate was dynamic and corresponding repair time was fixed.Paper [24] replaced the statistical availability with service continuity and proposed the probability Equation of service uninterrupted, which effectively promoted the actual availability of the service.By statistical methods, paper [25] obtained the accurate probability of service channel failure time exceeding the specified time according to a lot of simulation based a given network environment.But when the network environment changed, the simulation needed to be restarted which reduced the universality of this method.
The above methods mentioned are considering the actual availability which is more accurate than A SPA in routing.However, because the time to repair (TTR) of channel is changing with the environment and the geographical location, thus the assumption of a fixed TTR will limit the application scope of the above methods.On the other aspect, if the correlation among TTRs of different channels is not considered when selecting the primary and backup routing, the risk of primary and backup channels simultaneously being in failure will be increasing.Therefore, random TTR and its impact should be further considered for universal service channel failure risk routing mechanism.
Further studies showed that the occurrence of an electric power communication failure event and the fluctuations of service channel A APA (t) are interrelated and the fluctuations of A APA (t) is one of the root cause of electric power communication failure events [26].Because the fluctuations of A APA (t) has the random nature according to the interaction of many factors, the occurrence of the electric power communication failure event is a complex random process.Moreover, the fluctuations of A APA (t) makes the occurrence of the events that violate the availability rules (i.e., service channel failure events (SFE)) inevitable and hard to be precisely tracked.Therefore, it is necessary to precisely quantify the occurrence rule of SFE and SCVR.The preliminary work in paper [27] adopted the influence factors of service channel availability (FN and TTR) instead of A APA (t) to equivalently quantify SCVR.It just simply mentioned the idea without mathematical proof.However, it is proved to be an effective way to start with analysis of the influence factors of service channel availability and their relationship.
Thus, the study of distribution models of FN, failure arrival rate and TTR of transmission equipment and fiber cable and their internal relationship is the fundamental way to precisely quantify SCVR and control the influence degree of electric power communication failure events.According to the analysis above, the innovativeness of this paper are as follows.
(1) To precisely track the violation risk change of service channel under the condition that all of the FN, failure arrival rate and TTR are random, we deduce SCVRD model from the service channel violation risk model which is denoted by A APA (t) and denote SCVRD model by the probability of service channel cumulative failure duration exceeding the prescribed duration.We prove the deduction and simplify the SCVRD model with mathematical method.(2) Based on SCVRD model, the SCVRD-R mechanism is proposed to reduce the FN and failure duration of the electric power communication service.The goal is to minimize the violation risk caused by the random failures and TTR of transmission equipment and optical cable and thus to improve the availability of electric power communication service.

SCVRD Model
In this section, the differences between A SPA and A APA (t) are compared firstly and the existing problem of AAR-OS algorithm is analyzed.Subsequently, the SCVRD model is established according to the joint distribution of failure arrival rate and repair time of the transmission equipment and optical cable.

Differences Between A SPA and A APA (t)
Before establishing the SCVRD model, the varying rule of A APA (t) should be analyzed to identify the influence factors which cause the differences between A SPA and A APA (t).
A SPA is defined as the ratio of service un-interrupted duration in a statistical period of time [22].Let the graph G(V,E) represent the network topology, V denotes the node set and E denotes the edge set.A SPA can be calculated by Equations ( 1)-( 3) based on the statistical data at the end of each time period.
) and e ij ∈ E denotes the edge between node i and node j in G(V,E).A(v i ) denotes the statistical availability of the ith node.MTBF(v i ) and MTTR(v i ) respectively denote the mean time between failures (MTBF) and mean TTR (MTTR) of the i-th node.A(e ij ) denotes the statistical availability of the edge e ij .MTBF(e ij ) and MTTR(e ij ) respectively denote the MTBF and MTTR of the edge.MTBF and MTTR in Equations ( 2) and ( 3) are calculated by time between failures (TBF) and TTR in statistical period respectively.A APA (t) can be expressed as follows: Next, we will discuss the situations of A APA (t).Because the state of v i and e ij are always varying between normal state and failure state, the varying rule of A APA (t) can be analyzed by investigating the conversion process of TBF i and TTR i .This process can be divided into three phases for analysis: (1) During the period from initial time t 0 to the first failure occurred time t 1 , A APA (t) can be expressed as Equation ( 5): (2) When the kth (k ≥ 1) fault occurs, the items in Equation ( 4) are all constants except for TTR k , so A APA (t) varies only with the change of TTR k , expressed by Equation ( 6).a, b respectively, denote the cumulative availability time and the cumulative time when the kth fault occurs and both of them are constants in this phase.
The varying curve of A APA (t) in this phase is shown with the solid line in Figure 1.(3) When the channel returns to normal state from the kth failure, the items in Equation ( 4) are all constants except for TBF k , so A APA (t) varies only with the change of TBF k , expressed by Equation (7).c, d respectively denote the cumulative repair time and the cumulative time when the channel returns to normal state from the kth failure and both of them are constants in this phase.
The varying curve of A APA (t) in this phase is shown as the solid line in Figure 2.  According to the analysis above, the varying curve of A APA (t) changes with the channel's conversion process between normal states and failure states, as is shown in Figure 3.A Thr in Figure 3 denotes the availability threshold.As shown in Figure 3, the differences between A SPA , A APA (t) and A Thr are clearly illustrated.A SPA is the statistical mean value of A APA (t), which must be always greater than A Thr and is expressed as a horizontal dashed line in Figure 3.A SPA can be used to describe the overall availability of service channel from a global perspective but it cannot reflect the influence of each failure and repair time to the service.
A APA (t) is the actual availability of the service channel and its value can be calculated by (4) after each failure.However, since TBF i , TTR i , failure arrival rate and FN of each channel all have the nature of randomness, A APA (t) may be significantly higher than A Thr sometimes or may be lower than A Thr sometimes because of the sudden random failure.Therefore, the occurrence of SFE is inevitable.
If the degree of SCVR can be quantified, then the impact caused by SFE can be controlled, thereby the number and the duration of the impacted services can be reduced.To achieve the purpose, the SCVRD model is discussed and established in the following paragraphs.

SCVRD Model
As shown in Figure 3, SCVRD which represents the probability of service channel violation risk P SCVRD can be expressed by Equation (8): However, SCVRD is hard to be accurately quantified and tracked due to the fluctuations of A APA (t).Thus, we introduce cumulative repair time TR and corresponding repair time threshold T Thr in the statistical period to describe the service channel violation risk.
Thus, the probability of service channel violation risk based on availability is converted to the probability of cumulative repair time exceeding the corresponding repair time threshold in the statistical period.The equivalence relationship is shown in Figure 4. Next, the characteristics of SCVRD will be further analyzed.By further analyzing Equation (4), we can find that A APA (t) is co-determined by FN of the service channel and each failure TTR i .According to [23], TTR of v i and e ij is independent and they all obey the log-normal distribution as expressed by Equation (12).Therefore, TR subjects to the joint distribution of ∑ ∞ i=1 TTR i , which means that SFE is a random process and P SCVRD is a random probability model.Assume that FN of service channel obeys the Poisson distribution with average arrival rate λ, as expressed by Equation ( 13).Thus, the service channel violation risk caused by the kth failure P SCVRD (X = k, TR k > T Thr ) during the statistical period can be expressed by Equation ( 14): Because the failure occurrences of v i and e ij are independent, according to the conditional probability Equation and the n-fold convolution Equation, Equation ( 14) can be expanded to Equation (15): According to the total probability Equation, the Equation ( 15) can be expanded to (16).
Figure 5 is an example to illustrate the difference between P SCVRD (TR > T Thr ) and A SPA in the process of SCVRD quantization.During the calculation, we can find that P R (k) < 10 −6 when k > 10.As a result, when k ≤ 10, the precision requirement can be satisfied and this result is consistent with the fact that the repair time cannot be infinitely small.
Through this example, it can be found that although the A SPA of both SC N1−N2−N3 and SC N1−N4−N3 are all higher than A Thr , there is still a failure risk of about 0.389~0.426for service channel.In addition, although the A SPA of SC N1−N2−N3 and SC N1−N4−N3 are the same, their failure risks are still different: P SCVRD (TR Link1 > 0.72) = 0.426 and P SCVRD (TR Link2 > 0.72) = 0.389.Compared to A SPA , P SCVRD is more accurate to distinguish the differences of A APA (t) among the channels.Thus P SCVRD is more suitable in the quantization process of SFE.

Approximate Transformation of Violation Risk Distribution
In the electric power communication network, v i and e ij have the relatively high A SPA .It means that there are only several hours that the channel is in failure state during the thousands of hours' execution.Then it can be assumed that TTR of v i and e ij will not overlapped.In addition, the failure of v i and e ij is independent.Here the node violation risk P v i (TR > T Thr ) is analyzed and the violation risk model of edge P e ij (TR > T Thr ) can be got this way.Here the node violation risk model is introduced.As to node v i , the repair time TTR j , j = 1, .. obey the log-normal distribution with parameters µ i , δ i , because TTR j is a random variable, there is possibility that the TTR 1 of first failure may be bigger than T Thr and there also may be k failure occurs before k ∑ j=2 TTR j > T Thr .Therefore, from the total probability Equation we have: Equation ( 16) gives the probability distribution of violation risk P SCVRD but the distribution function is too complex to be a decision parameter.We will simplify the distribution function of ∑ k i=1 TTR i > T Thr according to the FN probability generating function, TTR i moment generating function and their distribution functions.So that P SCVRD can be used as a decision parameter in the channel routing process.The detailed derivation process is described as follows.Assume that the probability generating function of FN is expressed by: The moment generating function of TTR is expressed by: f (x) is the probability density function of TTR, then the moment generating function of down time (DT) is expressed by: Since TTR i is independent and obeys the same log-normal distribution, making use of the nature that the mathematical expectation of the product of the independent random variables equals to the product of the mathematical expectation of each independent random variable, the Equation ( 20) can be converted into Equation ( 21): From Equations ( 18) and ( 21), we get: Equation ( 18) is simplified to: From Equations ( 22) and ( 23), we get: From (24), the relationship between the moment generating function of total channel violation time and the moment generating function M TTR (t) of single violation time is explicit.Then we use the k-order moments of TTR to expand the moment generating function M DT (t).The k-order origin moment of TTR is expressed as Equation ( 25) and f (t) is probability distribution function of TTR.
From Equations ( 20) and ( 25), we get: Then: , m k represents the k-order origin moment of TTR and substituting Y into , then the moment generating function of Y is expressed as Equation ( 28): From Equation ( 26), we get: Put Equation (29) into Equation (28), we get: Since the moment generating function of the standard normal distribution is M N(0,1) (t) = e t 2 /2 , Equation (31) means that the moment generating function of Y converges to the moment generating function of standard normal distribution.Because the probability generating function or moment generating function of random variable is corresponding to its distribution function, the distribution function of Y converges to the standard normal distribution N(0, 1).Thus, the service channel violation risk P R is approximately calculated by converting Y = DT−λ i •m 1 λ i •m 2 and standard normal distribution function Φ(Y), which is expressed by Equation (32) By Equation (32), the expected value and variance of log-normal distribution can be respectively calculated as: µ i and σ i are the expected value and variance of normal distribution.According to Var(X) = E(X 2 ) − [E(X)] 2 , we get: Substituting Equations ( 35) and (36) into Equation (32), we get: Equation (37) gives the approximate method to calculate the violation risk of v i , which can be used as the decision-making parameter to select the route with minimum violation risk for service channel.

SCVRD Based Routing (SCVRD-R) Algorithm
Since the service channel is composed of v i and e ij , v i and e ij is independent of each other, the service channel violation risk minimization problem Min(P R (DT > T Thr )) can be converted into the shortest path problem with SCVRD(v i ) and SCVRD(e ij ) as the weight.To solve this problem, SCVRD-R algorithm based on Dijkstra is proposed in this paper.Here the nodes and edges that service channel passes are all called as channel section (CS), then in order to simplify the description of the algorithm, SCVRD(v i ) and SCVRD(e ij ) will be denoted as SCVRD(CS i ).The steps of SCVRD-R algorithm are described as follows.Figure 6 is the flowchart of SCVRD-R Algorithm 1.
(1) Initialize the graph G(N,E) and the aTTRibute vector of service (s, d, B, A, r).s, d, B, A and r(s,d) denote the source node, destination node, bandwidth requirements, availability threshold and the violation risk degree between s and d respectively.r 0 (s, d) = 0.When there is no channel between s and d, r(s,d) +∞.Set S = {s}, D = G-s.
(2) If the bandwidth of any path in the graph G(N,E) does not meet the bandwidth requirements B, then the path is deleted from G(N,E).Thus, a new sub-graph G'(N,E) is generated.
(3) For each CS i ∈ G(N, E), Calculate the probability of violation risk P CS i (TR > T Thr ) according to Equation (37).End for each.
(4) Choose the CS i with minimum violation risk probability P CS i from D, then D = D − CS i , S = S + CS i .
(5) Take CS i as the middle path, if the violation risk from s to CS k shrinks, then modify r(s, CS k ).

Simulation and Results
In this paper, two ways have been used to analyze the performance of the algorithm.Firstly, the Monte Carlo analysis method is adopted to simulate the random failures of service channel.The result of this method is compared to AAR-OS algorithm.Secondly, based on the backbone topology of one provincial electric power company, we randomly generate the distribution parameters of failure arrival rate and repair time for v i and e ij .Then the performance of SCVRD-R algorithm, AAR-OS algorithm and risk-aware provisioning (RAP) algorithm [23] is compared under the intensity of 10 kinds of services.
To reduce the impact on experiment results which is brought by our algorithm's random nature, we get the probability of service channel violation risk through repeating experiment 10,000 times and the method that the channel monthly statistical average failure times are divided by the experiment times.The statistical average data and the routing results of two algorithms are shown in Table 1.In the Table 1, the thick arrows denote the route selected by SCVRD-R algorithm, while the thin arrows denote the route selected by AAR-OS algorithm.From Table 1, it can be seen that the violation risk of the two channels are fluctuating.If AAR-OS algorithm is adopted, even when the violation risk of SC N1−N2−N3 is lower than SC N1−N4−N3 (for example on month 6 and 7), AAR-OS will always choose the channel SC N1−N4−N3 for the service S i because the channel SC N1−N4−N3 has slightly higher statistical availability than SC N1−N2−N3 .If SCVRD-R algorithm proposed in this paper is adopted, then the lower risk route will be selected for S i according to the change of the violation risk value.Thus, the number of the service violation will be reduced and its A APA (t) will be enhanced.

Network Simulation
The network simulation scenario is set up based on the backbone network topology of one provincial electric power company, as shown in Figure 7.The raw data about topology, business and actual routing are all provided by this electric power company because of the research cooperation between this company and our research team.
The network simulation topology contains 74 nodes and 104 edges.The purpose of the simulation experiment is to compare the failure rate of the service channel routed by SCVRD-R algorithm, AAR-OS algorithm and RAP algorithm in 12 months in Table 2.  Here, we divide 30 services into ten kinds of service intensities as an instance to evaluate the algorithm performance in the different ten service intensities.However, the partition of service intensities is not only this one way.Assume that the service is generated between nodes in 2 hops of 14 aggregation nodes (marked with circle in Figure 7) and the services have the same availability threshold of 99.5%.The failure arrival rates of nodes and edges are randomly generated in the interval of [0.006, 1.5] and their repair time are randomly generated in the interval of [0.4,0.003].The service channel failure rate is calculated using the ratio between the number of failure service channels and the number of total service channels.The bandwidth of each optical cable is 5GB.At present the bandwidth in the electric power communication transmission network is redundant.Besides, the SDH electric power communication service channel uses fixed bandwidth.So, we temporarily do not consider the variable bandwidth of cable.
The simulation results are shown in Figure 8, which shows that the average service channel failure rate of SCVRD-R algorithm is respectively about 15% and 6% lower than that of AAR-OS algorithm and RAP algorithm under the different service intensities.Other service intensities are similar.The failure rate in the various service intensities has nothing to do with the change of month, just randomly changes within months.The temporal resolution of data is not necessarily fixed for month, which means that week or other durations can also be set as the temporal resolution of data.

Service Channel Failure Rate Comparison
Figure 8a shows the service channel failure rates of SCVRD-R algorithm, AAR-OS algorithm and RAP algorithm in 12 months under service intensity 1.Because SCVRD-R algorithm finds route according to SCVRD, the service channel can adjust the route based on the distribution of failure arrival rate and TTR to avoid the fault paths, which can reduce the service channel failure rate.The results in Figure 8b,c are similar to the results in Figure 8a.

Failure Rate Comparison Under the Saturation Status
Figure 8d shows the service-channel failure rate of SCVRD-R algorithm, AAR-OS algorithm and RAP algorithm in 12 months under service intensity 4.Under this service strength, the traffic carried by parts of the optic cables may reach the bandwidth limit and the network will begin to be the saturation status.In this case, the service channel with low SCVR cannot carry the service which is transferred from the service channel with relatively high SCVR.Then, parts of the services cannot select the route according to the lowest violation risk, resulting in the ascending service channel failure rate.In addition, with the growing of number of services, the ability of these three algorithms to control A APA (t) and SCVRD of channels is gradually weakened.So, the results of the failure rate turn to convergence and will finally converge to the failure rate without routing risk control as the number of services grows.

Conclusions
On account of the problems in electric power communication service route planning, a probability model of service channel violation risk, named SCVRD model, is proposed in this paper.Generally, the A SPA is calculated by mean value and A APA (t) is calculated based on some assumption, both of them cannot precisely track the violation risk change of service channel which is caused by the random failure under random TTR condition.To solve this problem, we deduce SCVRD model from the service channel violation risk model which is usually denoted by A APA (t) and denote SCVRD model by the probability of service channel cumulative failure duration exceeding the prescribed duration.The deduction is proved and then SCVRD model is simplified using mathematical method.Based on SCVRD, a service channel violation risk degree routing algorithm, named SCVRD-R algorithm, is proposed to reduce the risk caused by random failure of transmission equipment and optical cable and improve the availability of electric power communication service.Finally, the simulation results show that the average service channel failure rate of AAR-OS algorithm and RAP algorithm are respectively reduced by 15% and 6%.
In future work, we plan to further investigate how the failure rate convergence changes when optic cables reach the bandwidth limit.Furthermore, we intend to figure out the accurate numerical relationship among service intensity, service channel bandwidth and failure rate convergence to actually guide the service routing planning and optimization in electric power communication network.In many cases, the success of smart grid depends on their ability to support real-time decision-making.Therefore, we would make further efforts on the exploration of our proposed methodology in medium and short-term planning scenarios.

Figure 1 .
Figure 1.The varying curve of A APA (t) with TTR k .

Figure 2 .
Figure 2. The varying curve of A APA (t) with TBF k .

Figure 3 .
Figure 3.The integral varying curve of A APA (t).

( 6 )
Repeat step 4 and step 5 until CS i = d.(7) If the new route is different from the original one, then a new service channel is created.Repeat Step 2 to Step 7.

Table 1 .
Comparison Table of Violation Risk and Routing Result.

Table 2 .
Service Intensity Distribution Table.