The relative biologic effectiveness versus linear energy transfer curve as an output-input relation for linear cellular systems

  • Received: 01 September 2008 Accepted: 29 June 2018 Published: 01 June 2009
  • MSC : Primary: 92C99, 93A30; Secondary: 92B05.

  • Experiments have established that different radiation types have different magnitudes of biological response. When biological response is defined in terms of the Relative Biologic Effectiveness (RBE) and different radiation type is characterized by Linear Energy Transfer (LET), the plot of the RBE versus LET (RBE-LET) curve shows RBE to increase with increasing LET, to reach a maximum, and to decrease with further increasing LET. Perhaps due to the descriptive nature of biology, most quantitative models for the RBE-LET curve ignore the reality of the underlying molecular biology. On the other hand, the molecular basis for the RBE-LET curve is not completely known despite recent efforts.
       Here we introduce a differential equation formulation for a signal-and-system model that sees cells as systems, different radiation types as input, and cellular responses as output. Because of scant knowledge of the underlying biochemical network, the current version is necessarily a work in progress. It explains the RBE-LET curve using not just input parameters but also systems internal state parameters. These systems internal state parameters represent parts of a biochemical network within a cell. Although multiple biochemical parts may well be involved, the shape of the RBE-LET curve is reproduced when only three system parameters are related to three biochemical parts: the molecular machinery for DNA double strand break repair; the molecular pathways for handling oxidative stress; and the radiolytic products of the cellular water.
       Despite being a simplified ''toy model,'' changes in the systems state parameters lead to model curves that are refutable in a modern molecular biology laboratory. As the parts in the biochemical network of the radiation response are being further elucidated, this model can incorporate new systems state parameters to allow a more accurate fit.

    Citation: Quoc T. Luu, Paul DuChateau. The relative biologic effectiveness versus linear energy transfercurve as an output-input relation for linear cellular systems[J]. Mathematical Biosciences and Engineering, 2009, 6(3): 591-602. doi: 10.3934/mbe.2009.6.591

    Related Papers:

    [1] Tung-Sheng Zhan, Yih-Der Lee, Jheng-Lun Jiang . Optimal OCR coordination in a high penetration distribution power system using a refined immune algorithm with an auto-tuning reproductive mechanism. AIMS Energy, 2024, 12(6): 1225-1263. doi: 10.3934/energy.2024056
    [2] Pooja Khandare, Sanjay Deokar, Arati Dixit . Relay Coordination and Optimization techniques using DWT-Differentiation Algorithms for Fault Detection in Microgrid. AIMS Energy, 2020, 8(4): 563-579. doi: 10.3934/energy.2020.4.563
    [3] Amandeep Gill, Pushpendra Singh, Jalpa H. Jobanputra, Mohan Lal Kolhe . Placement analysis of combined renewable and conventional distributed energy resources within a radial distribution network. AIMS Energy, 2022, 10(6): 1216-1229. doi: 10.3934/energy.2022057
    [4] Sunimerjit Kaur, Yadwinder Singh Brar, Jaspreet Singh Dhillon . Multi-objective real-time integrated solar-wind-thermal power dispatch by using meta-heuristic technique. AIMS Energy, 2022, 10(4): 943-971. doi: 10.3934/energy.2022043
    [5] Minh Y Nguyen . Optimal voltage controls of distribution systems with OLTC and shunt capacitors by modified particle swarm optimization: A case study. AIMS Energy, 2019, 7(6): 883-900. doi: 10.3934/energy.2019.6.883
    [6] W. E. P. Sampath Ediriweera, N. W. A. Lidula . Design and protection of microgrid clusters: A comprehensive review. AIMS Energy, 2022, 10(3): 375-411. doi: 10.3934/energy.2022020
    [7] Yan Li, Yaheng Su, Qixin Zhao, Bala Wuda, Kaibo Qu, Lei Tang . An electricity price optimization model considering time-of-use and active distribution network efficiency improvements. AIMS Energy, 2025, 13(1): 13-34. doi: 10.3934/energy.2025002
    [8] Shafiuzzaman Khan Khadem, Malabika Basu, Michael F. Conlon . Capacity enhancement and flexible operation of unified power quality conditioner in smart and microgrid network. AIMS Energy, 2018, 6(1): 49-69. doi: 10.3934/energy.2018.1.49
    [9] Sabir Rustemli, Zeki İlcihan, Gökhan Sahin, Wilfried G. J. H. M. van Sark . A novel design and simulation of a mechanical coordinate based photovoltaic solar tracking system. AIMS Energy, 2023, 11(5): 753-773. doi: 10.3934/energy.2023037
    [10] Steven B. Sherman, Zachary P. Cano, Michael Fowler, Zhongwei Chen . Range-extending Zinc-air battery for electric vehicle. AIMS Energy, 2018, 6(1): 121-145. doi: 10.3934/energy.2018.1.121
  • Experiments have established that different radiation types have different magnitudes of biological response. When biological response is defined in terms of the Relative Biologic Effectiveness (RBE) and different radiation type is characterized by Linear Energy Transfer (LET), the plot of the RBE versus LET (RBE-LET) curve shows RBE to increase with increasing LET, to reach a maximum, and to decrease with further increasing LET. Perhaps due to the descriptive nature of biology, most quantitative models for the RBE-LET curve ignore the reality of the underlying molecular biology. On the other hand, the molecular basis for the RBE-LET curve is not completely known despite recent efforts.
       Here we introduce a differential equation formulation for a signal-and-system model that sees cells as systems, different radiation types as input, and cellular responses as output. Because of scant knowledge of the underlying biochemical network, the current version is necessarily a work in progress. It explains the RBE-LET curve using not just input parameters but also systems internal state parameters. These systems internal state parameters represent parts of a biochemical network within a cell. Although multiple biochemical parts may well be involved, the shape of the RBE-LET curve is reproduced when only three system parameters are related to three biochemical parts: the molecular machinery for DNA double strand break repair; the molecular pathways for handling oxidative stress; and the radiolytic products of the cellular water.
       Despite being a simplified ''toy model,'' changes in the systems state parameters lead to model curves that are refutable in a modern molecular biology laboratory. As the parts in the biochemical network of the radiation response are being further elucidated, this model can incorporate new systems state parameters to allow a more accurate fit.


    Integrating distributed generators into power systems has emerged as a transformative force in today's dynamic and ever-evolving energy landscape. Distributed generators, encompassing a wide range of renewable and non-renewable sources such as solar panels, wind turbines, diesel and microturbines, are increasingly prevalent in modern power networks. While these DGs promise cleaner and more sustainable energy generation, their integration presents unique challenges, particularly concerning the coordination of power system relays. Power relays play a critical role in safeguarding power systems by detecting and isolating faults, thereby ensuring the continuous delivery of electricity to consumers. However, traditional relay coordination practices face new complexities as DG penetration grows. This article delves into the intricate issue of power system relay coordination and explores the profound influence of DGs on this essential aspect of grid operation.

    Power system relay coordination is a well-established discipline aimed at designing relay protection schemes that ensure the selective operation of protective devices during fault conditions. Traditionally, this coordination has focused on centralized generation sources, allowing for precise adjustment of relay settings based on a predetermined hierarchy. In this classical approach, generators are typically large, centralized units whose behavior can be predicted with a high degree of accuracy. Relay coordination is achieved by adjusting settings such as time-current characteristics, fault-clearing times, and coordination margins to ensure that the nearest relay responds to a fault while isolating it, leaving the rest of the system intact.

    The integration of DGs challenges the traditional relay coordination paradigm in several ways:

    Intermittency and variability: DGs, especially renewable sources like solar and wind, are inherently variable and intermittent. Their power output depends on weather conditions, time of day, and other factors, making it difficult to predict and coordinate relay settings accurately.

    Bi-directional power flow: Unlike conventional generators, DGs can inject power into the grid and draw power from it. This bi-directional power flow introduces the potential for reverse power flows during fault conditions, which can complicate relay coordination strategies [1].

    Islanding operation: DGs can create isolated sections or "islands" within the power grid during disturbances or outages. These islands pose unique challenges for relay coordination, as traditional settings may lead to unnecessary tripping of distributed generators.

    Protection coordination time scales: DGs often operate on shorter time scales than centralized generators, responding rapidly to changes in grid conditions. Relay coordination must now consider these fast-acting devices while protecting against grid faults.

    To deal with these challenges, power system engineers are developing innovative solutions considering the generators' distributed nature. Adaptive relay coordination schemes, which use real-time data and advanced algorithms, are gaining prominence. These schemes allow for dynamic adjustments of relay settings based on the actual operating conditions of the grid, including the presence and behavior of DGs. Adaptive relay coordination promises enhanced reliability and resilience in power systems, ensuring that protective devices respond appropriately to faults while minimizing unnecessary tripping of distributed generators. Moreover, these strategies can improve the overall efficiency and stability of the grid, making it more accommodating to the variability of renewable energy sources.

    As mentioned, DG can impact the protection coordination of over-current relays on distribution systems. [2] proposed a "protection coordination index" (PCI), which can serve as an effective measure when planning the protection of meshed distribution systems with DGs. Then, a two-phase non-linear programming (NLP) optimization method is applied to determine the PCI by calculating variations in the maximum DG penetration level with changes in the protection coordination time interval. [3,4] resents the multi-agent theory to deal with protection coordination. In [3], the effects of DG grid connection have been considered, and the Java Agent Development Framework platform open software is used to adjust the protection coordination curve. [4] considers the issue of feeder reorganization when DG is connected to the grid and uses the signals measured by the phasor measuring unit (PMU) to dynamically adjust the relay settings by the central control system in the power station. The protection coordination problem was regarded as an optimization model in traditional mathematics, and the mixed integer nonlinear programming (MINLP) and mixed integer linear programming (MILP) were used to decide the TMS setting for fault current variation of the DGs integration in [5,6]. [7] uses the interior-point (IP) algorithm to solve the optimal adaptive protection coordination strategy under high penetration of green generation and network reconfiguration. Many contributions to the application of metaheuristic intelligence-based random search techniques applied to the adaptive coordination decision are represented here by document [8], which uses firefly optimization theory (FA) and evolutionary programming method (EP) to simultaneously consider feeder reorganization, DGs integration, DGs capacity planning and optimal adaptive protection coordination strategy. [9] utilizes the genetic algorithm (GA) to deal with the optimal DG placement problem to maximize the penetration level of DG in MV distribution networks without changing the original relay protection parameter. GA is used to discover the optimal sizes and locations of DG. Coincidentally, [10] also worked on the optimal placement issue of DGs and using optimization methodology to classify the cases, either coordination holds or coordination lost, by considering DG location changing and DC capacity variation. [11] presents an approach for protection coordination in microgrids that incorporates non-standard characteristic features of directional over-current relays (DOCRs), and the coordination model corresponds to an MINLP problem. There are four famous metaheuristic techniques, i.e., particle swarm optimization (PSO), GA, teaching-learning based optimization (TLBO) algorithm and shuffled frog leaping algorithm (SFLA), implemented for solving the optimal coordination problem. [12] proposes the use of adaptive protection, using voltage and current measurement of the system to overcome the challenges of overcurrent protection in distribution systems with DGs. The trip characteristics of OCRs are updated by detecting system operating states (grid-connected or island) and the faulted section.

    References [13,14,15,16] apply artificial neural networks (ANN), using known training samples to train ANN weight values and construct a decision-making neural network system. When the system status changes–such as DG ON/OFF–the line after switching or line reorganization, the change in line current is used as the input value of the system and the corresponding relay setting adjustment value is output through the network decision-making process. In [13], an adaptive distance relaying method was presented. The multi-layered perceptron (MLP), also known as the error backpropagation network, was applied to estimate the actual power system condition and to calculate the appropriate tripping impedance under varying power system conditions. [14] and [15] apply ANN based on the Levenberg-Marquardt algorithm to model the DOCR characteristics where the trained ANN model can compute TMS, the pickup current setting (PCS) and the operation time of each DOCR in terms of changes caused by DG integration in the radial distribution system and looped distribution system respectively.

    In this paper, the relay protection coordination issue was regarded as an optimization model by minimizing the objective function of the accumulated sum of all relays' operating times in the distribution power system subject to the TMS limits digital OCR and CTI restriction of relay operation time. The optimal protection coordination model is applied to a power system by considering several conditions of DG incorporation. According to DGs' various operating conditions, each OCR's TMS setting value is investigated. First, the author proposes the concept of tracking paths. Tracking paths can be divided into two categories: the main feeder path and the branch line path. The trace path will be represented as a set of relay numbers in order according to the topology of feeder connectivity. Then, the TMS of each relay can be used to calculate the OT of the relay, and the cumulative OT of the relay on each tracking path will be used to determine the fitness of the TMS setting combination. Furthermore, the relay protection coordination problem can be modeled to minimize the accumulated summation of all primary and backup relay OT subject to the CTI limitation.

    Second, as the author's earlier contribution, an adaptive evolution mechanism with a self-adjusting scheme was proposed in [17]. The self-adjusting scheme was combined successfully with the PSO algorithm to improve the shortcomings of premature convergence. Then the modified particle swarm optimization algorithm with adaptive self-cognition and society operation scheme was proposed in this paper. The MPSO-ASSOS determines TMS for each relay on the tracking path, subject to CTI and OT limitation. A 16-bus test MV system with five DGs will be applied to test the adaptive protection coordination approach proposed in this paper. The simulation results show that the proposed method can quickly and stably obtain the solution with maximum fitness (i.e., minimum overall OT) and overcome the impact on protection coordination settings of the system with DG incorporation.

    As mentioned in the previous section, the power injection of RESs can affect the protection coordination of over-current relays (OCR) for the RDNs; there is a simple radial feeder, with six buses, five feeders, four DGs and five OCRs with paired circuit breakers (Figure 1).

    Figure 1.  Parts of single line diagram of an RDN with DGs.

    The circuit breakers, CB 18, CB 19, CB 20, and CB 21, switch and simulate whether the DG is connected to the power distribution system. The fault analysis and protection coordination calculation of the demonstrated radial feeder was performed by using ETAP software. The parameter setting of each relay of the mentioned RDN in Figure 1 includes PCS and TMS, and these are assumed to carry out protection coordination calculation according to formal design steps in advance, as shown in Table 1. From Table 1, it can be found that the TMS of the relay is designed to be arranged in a sequential order of value, and the larger the value of the TMS, the longer the operation time of the relay. Therefore, it can think of an OCR with a larger TMS as the backup protection of an OCR with a smaller TMS, which is called primary protection; for instance, because Relay 8 is the primary protection of the line section 7–8, Relay 7 will be the backup protection of the Relay 8. Similarly, Relay 6 is the backup protection of Relay 7, and Relay 5 serves as the backup relay of Relay 6.

    Table 1.  Parameter setting of each relay for the RDN.
    Line section Corresponding relay TMS PCS (A)
    4–5 Relay 5 0.37 800
    5–6 Relay 6 0.26 800
    6–7 Relay 7 0.17 800
    7–8 Relay 8 0.08 800

     | Show Table
    DownLoad: CSV

    Table 2 shows the fault current distribution of the three-phase grounding fault at Bus 6 for the condition without DG injection and with DG1 integration. It can be seen from Figure 1 that only DG1 is incorporated into the end of this feeder. Table 2 shows that the source of the fault current flowing through Relay 5 and Relay 6 is the power source of this feeder, and due to the integration of DG1, it will generate an opposite-direction load flow from the power source. For this reason, the value of fault current through Relay 7 and Relay 8 when DG1 is integrated will be higher than if DG1 is not incorporated. Then, the higher fault current will affect the protection coordination mechanism of the original design. The protection coordination impact of DG1 integration for the fault that occurred at Bus 6 is shown in Figure 2. The operation time of Relay 5 and Relay 6 will not change due to the integration of DG1, which are 0.894 seconds and 0.646 seconds, respectively, as shown at points B and A. However, after DG1 is incorporated, the fault current value flowing through Relay 7 and Relay 8 becomes higher, so that the operation time of Relay 7 and Relay 8 can be calculated as 1.08 seconds and 0.557 seconds respectively, as shown at points D and C., according to the computed operation time, the order of operation sequence of Relay 5–8 will be changed to Relay 8, Relay 6, Relay 7 and Relay 5, which has caused the original protection coordination out of order. Simply put, when Bus 6 fails, Relay 6 should trip the circuit breaker, but due to the incorporation of DG1, Relay 8 will act instead and the protection coordination is lost.

    Table 2.  The fault current distribution of fault occurs at Bus 6.
    Relay No. Fault current without
    DGs injection (kA)
    Fault current with
    DG1 injection (kA)
    5 12.44 12.44
    6 12.44 12.44
    7 0.791 2.474
    8 0.483 2.167

     | Show Table
    DownLoad: CSV
    Figure 2.  Protection coordination impact of DG1 integration when a fault occurs at Bus 6.

    A further example illustrates the situation when the power connection point of DG is not located at the end of the feeder. Considering the feeder in Figure 1, DG1 and DG2 are connected at Bus 8 and Bus 6, respectively, to the supply power at the same time. The resulting fault current distribution when the three-phase grounding fault occurred at Bus 5 is shown in Table 3. The fault current flowing to Bus 5 comes from DG1, DG2 and the power company, and it is not difficult to find that only the fault current flowing through Relay 5 is provided by the power company, regardless of whether the feeder is incorporated into DG or not. When DG1 and DG2 are combined, the fault current of Relay 6, Relay 7 and Relay 8 is increasing due to the contribution of DG1 and DG2. The protection coordination impact of two DGs integration for the fault that occurred at Bus 5 is shown in Figure 3. The operation time of Relay 5 remains the same value with the integration of DG1 and DG2, which are 0.854 seconds, as shown at point F. in Figure 3, however, after two DGs are incorporated, the fault current value flowing through Relay 6 to Relay 8 becomes higher, so that the operation time of three relays can be found as 1.03, 1.11 and 0.537 seconds, respectively, and are marked at points E, H and G in Figure 3. The order of operation sequence of Relay 5 to Relay 8 will be changed to Relay 8, Relay 5, Relay 6 and Relay 7, which has also caused the original protection coordination to be out of order. From the results of the above two case studies, it can be concluded that the existing relay settings need to be reviewed and updated to deal with the effect of DER integration to reduce the occurrence of relay miscoordination.

    Table 3.  The fault current distribution of fault occurs at Bus 5.
    Relay Fault current without DGs injection (kA) Fault current with DG1 and DG2 injection (kA)
    5 14.125 14.125
    6 1.302 4.582
    7 0.784 2.403
    8 0.479 2.105

     | Show Table
    DownLoad: CSV
    Figure 3.  Protection coordination impact of DG1 and DG2 integration when a fault occurs at Bus 5.

    The relay protection coordination problem can be formulated as a nonlinear programming problem for minimizing all relay operating time summation. A relay's operation time depends on its PCS and TMS for a given fault current through the relay. Obtaining suitable PCS and TMS parameter values for all the relays installed in the distribution system and appropriately coordinating all the relays' operations is challenging [18,19]. To simplify the protection coordination optimization problem, this paper will assume that the PCS was given and has undergone a reasonably rigorous design and calculation process. Then, the setting of TMS can be obtained by formulating the OCRs coordination issue as an optimization problem with the objective function expressed as

    MinizeObjF=TNk=1FNj=1SNi=1Wki,j×TMSi (1)
    Wki,j=βi(Ifki,jPCSi)αi1, (2)

    where

    ObjF : Object function of calculating the total summation of operation time for each specified relay.

    Wki,j : The weighting factor of relay i when the jth fault occurs in the kth coordination search path condition.

    i : The relay identifier.

    j : The fault location identifier.

    k : The coordination tracking path identifier (relay connectivity).

    FN : The total number of fault locations is considered.

    TN : The total number of coordination tracking paths.

    SN : The total relay number of the set of tracking paths.

    TMSi : TMS parameter of relay i.

    PCSi : The given PCS parameter of relay i.

    Ifki,j : Fault current through the relay i when the jth fault occurs in the kth coordination tracking path.

    αi,βi : Parameters of the standard normal time-current characteristic function [4] for relay i.

    Equation 1 is an objective function to minimize the summation of the operation times of the primary and backup relays when faults occur at different pre-setting positions. According to different fault locations, TN different protection coordination tracking paths can be determined and these paths may partially overlap, so the relay operation time on the overlapping paths will not be accumulated repeatedly. The normal inverse time-current characteristic function of the IEC standard [20], as shown in Eq 2, was only considered in this paper for evaluating the relay's operation time, and the fault current in Eq 2 can be simulated previously with several power engineering software, such as ETAP, PSS/E, etc. The parameters, αi and βi, of the standard normal time-current curve are set to the value of 0.14 and 0.02 in this paper, respectively.

    Coordination time interval (CTI) is the operation time difference between the primary and backup relays. When the CTI of all the primary and backup paired relays in the power system is within a limit value or an allowable value range, it represents that protection relays are well-coordinated. The CTI evaluation and its limitations are clearly described in Eqs 3 and 4. Then, the selection of the TMS value will depend on the brand and model of the selected relay, and there are also restrictions on its upper limit and lower limit, as shown in Eq 5.

    CTIki,j=tkib,jtkip,j (3)
    CTIminCTIki,jCTImax (4)
    TMSminTMSkiTMSmax (5)

    where

    CTIki,j : Object function of calculating the total summation of operation time for each specified relay.

    tkip,j : Operation time of primary relay ip when the jth fault occurs in the kth coordination tracking path.

    tkib,j : Operation time of backup relay ip when the jth fault occurs in the kth coordination tracking path.

    CTImin : The lower CTI bound is set to 0.2 seconds in this paper.

    CTImax : The upper CTI bound is set to 0.24 seconds in this paper.

    TMSmin : The lower bound of TMS is set to 0.05 in this paper.

    TMSmax : The upper bound of TMS is set to 1.0 in this paper.

    In the procedure of the PSO, a given objective function is optimized by simulating the natural behavior of birds (particles) flocking. Each particle (pbest) knows its best value and position so far; this information is the analogy to the personal experiences of each particle. Moreover, each particle knows the group's best value (gbest) among pbests so far. This information is an analogy to the knowledge of how the other particles around them have performed [21,22,23].

    In the original model of PSO, the velocity v will control the direction and distance of particle movement and can be classified into three operating modes as follows [21,22,23], and for the sake of distinction, it is called simple self-cognition and society operation scheme (SSSOS) here

    1) Self-cognition operation mode

    v(g+1)i=v(g)i+C1rand1(pbest(g)ip(g)i) (6)

    2) Society operation mode

    v(g+1)i=v(g)i+C2rand2(gbest(g)ip(g)i) (7)

    3) Hybrid operation mode

    v(g+1)i=v(g)i+C1rand1(pbest(g)ip(g)i)+C2rand2(gbest(g)ip(g)i) (8)

    where

    rand1 : The uniform random number in (0, 1) for the self-cognition process.

    rand2 : The uniform random number in (0, 1) for the society process.

    g : The current iteration numbers.

    v(g)i : The velocity of particle i at iteration g.

    C1, C2 : Weighting factors.

    p(g)i : Current position of particle i at iteration g.

    pbest(g)i : pbest of particle i at iteration g.

    gbest(g)i : gbest of particle i at iteration g.

    A certain velocity, which gradually gets close to pbest and gbest, can be calculated using the above equation. The particle position of the next iteration (searching point within the solution space) with a unit time can be modified by the following equation:

    p(g+1)i=p(g)i+v(g+1)i (9)

    The PSO algorithm prefers to use the hybrid operation to generate the next generation of new solution combinations in the SSSOS. A higher weighting factor, C1, allows exploring solution territory around the local best solution and controlling the probability of introducing new particles. The weighting factor C2 controls the velocity toward the area of the global best solution space. Then, the hybrid operation has the functions of the above two operations to generate the next generation of new solution combinations. If C1 is too high, the solution might settle at a local optimum. On the contrary, a lower rate could generate too many possibilities of uncertainty. The offspring particle loses its resemblance to the parent particles, and the algorithm will not learn from the past and could become unstable. Choosing a suitable weighting factor of self-cognition and society operation for PSO is a dilemma. This article proposes an adaptive self-cognition and society operation scheme below to avoid such difficulty. The ASSOS scheme was illustrated as follows,

    1) Assume that the operation mechanism of each particle generated in the last iteration process is recorded. There are three operation mechanisms: self-cognition, society and hybrid.

    2) Select one particle orderly for producing a next-generation new particle according to

    a) If rand1 < P1(g) and rand2 < P2(g): neither self-cognition nor society process be executed;

    b) If rand1 ≥ P1(g) and rand2 < P2(g): execute self-cognition process only;

    c) If rand1 < P1(g) and rand2 ≥ P2(g): execute society process only;

    d) If rand1 ≥ P1(g) and rand2 ≥ P2(g): hybrid process be executed.

    where

    P1(g) : control parameter of self-cognition with initial value P1(0) = 0.5 and 0 P1 1.

    P2(g) : control parameter of society process with initial value P2(0) = 0.5 and 0 P2 1.

    The next-generation new particles will be generated until all parent particles are processed. Figure 4 shows the initial relationship between self-cognition, society, and hybrid operation in PSO, which can be performed to generate next-generation particles in equal initial probability. If the random number of self-cognition(rand1) or society(rand2) processes is less than the corresponding control parameter, the related procedures will not be implemented. Self-cognition operation will play a more important role than that in PSO since the mutation can explore new regions. Suppose the search is very close to the local or global optimum. In that case, self-cognition may need to become dominant, especially in the absence of the critical suitable particles in a generation. Since all the procedures are randomly operated, there is no telling which is better.

    Figure 4.  Initial probability map of self-cognition, society and hybrid operation in ASSOS.

    3) A competition mechanism is thus implemented in the search process according to the fitness score. For instance, if the best current particle comes from a hybrid process, there is more likelihood for this procedure to generate a better solution for the next iteration. The area of the hybrid procedure must get bigger by reducing P1(g) and P2(g) to expand the probability, as shown in Figure 5.

    Figure 5.  Variation for decreasing hybrid operation probability.

    If the best fitness of the iteration g-1 is greater than the generation g, i.e., Fitness(g1)min>Fitness(g)min comes from hybrid procedure, both of control parameters will increase. We have

    P(g+1)1=P(g)1+D1=P(g)1+(K1gmax) (10)
    P(g+1)2=P(g)2+D2=P(g)2+(K2gmax) (11)

    K1 and K2 are the regulating factors, and in general, K1 < K2. gmax is the maximum iteration number. Figure 5 shows the decrease in the probability of hybrid operation area. On the contrary, there is a greater likelihood for the other two procedures to generate better particles. Both control parameters must decrease to subjoin the probability, as shown in Figure 6.

    Figure 6.  Variation for increasing hybrid operation probability.

    If the best solution remains the same, hybrid operation needs to be held back to recover the related area. If the best fitness of generation g-1 is less than that of generation g, i.e., Fitness(g1)minFitness(g)min comes from the hybrid operation, the control parameters will decrease as the following manner,

    P(g+1)1=P(g)1D1=P(g)1(K1gmax) (12)
    P(g+1)2=P(g)2D2=P(g)2(K2gmax) (13)

    4) If the best fitness of the iteration g-1 is greater than the iteration g, (i.e., Fitness(g1)min>Fitness(g)min) comes from only self-cognition procedure, then the control parameter will increase by using Eq 10. Conversely, if Fitness(g1)minFitness(g)min comes from self-cognition, then the control parameters will decrease by using Eq 12. In this situation, the control parameter P2 is fixed. The probability variation of self-cognition is shown illustrated in Figures 7 and 8.

    Figure 7.  Variation for decreasing self-cognition operation probability.
    Figure 8.  Variation for increasing self-cognition operation probability.

    5) If the best fitness of the iteration g-1 is greater than the iteration g, i.e., Fitness(g1)min>Fitness(g)min, comes from only society procedure, then the control parameter will increase by using Eq 11. Conversely, if Fitness(g1)minFitness(g)min comes from society operation, then the control parameters will decrease by using Eq 13. In this situation, the control parameter P1 is fixed. The probability variation is shown illustrated in Figures 9 and 10.

    Figure 9.  Variation for decreasing society operation probability.
    Figure 10.  Variation for increasing society operation probability.

    Particle generation is the first procedure of the PSO, and it has been mentioned in the previous section that the relay protection coordination problem can be regarded as a nonlinear programming problem. Therefore, in this paper, the MPSO-ASSOS will be used to combine the concept of nonlinear programming and apply it to the TMS parameters setting for each relay. First, the protection coordination tracking path in the system must be determined. The relay numbers on the path must be arranged in sequence according to the order of the main relay and the backup relay, and the set of relay numbers can be obtained; for instance, two tracking paths can be determined in Figure 1 as Path1: { 5 6 7 8 } and Path2: { 7 9 }. Path 1 can be called as the main feeder, and Path 2 is the branch separated from the main feeder. The TMS of relays on the main feeder and branch lines need to be designed in different manuals as explained in the following sections.

    As a main feeder path set example, Path: { x y z }, the last relay number z can be chosen as the starting point for protection coordination and its operation time can be set to the minimum value allowed by power regulations, such as 0.2 seconds. This relay can be regarded as the primary relay due to the fault current and PCS are given. Then, its TMS can be calculated using Eq 2 in the following way,

    TMSz=tkz/Wkz,j (14)

    tkz is the presetting operation time of the starting point relay. Then, the TMSy of the up-stream (backup) relay y will be estimated as

    TMSy=random(UPBy,LOWBy) (15)
    UPBy=tkz+CTImaxWky,j (16)
    LOWBy=tkz+CTIminWky,j, (17)

    where UPBy is the upper bound of TMSy, it can be obtained by adding the maximum CTI allowable value to the operation time of the previous relay. Similarly, LOWBy is the lower bound can be calculated by adding the minimum CTI allowable value to the operating time of the previous relay. The operation time tky can be found by

    tky=Wky,j×TMSy (18)

    Since relay x is the next up-stream (backup) protection of relay y, TMSx can be calculated repeat according to Eqs 15–18, as long as the subscript y in these equations are changed to x.

    For instance, the tracking path mentioned earlier in this section, Path1: { 5 6 7 8 }, if the operation time of relay 8, tk8, is set to 0.2 sec., the TMS8 can be obtained by using Eq 2 and Eq 14. Then, the TMS7 could be generated randomly within upper and lower bound by

    TMS7=random(UPB7,LOWB7)

    with

    UPB7=0.2+CTImaxWk7,j
    LOWB7=0.2+CTIminWk7,j

    The operation time tk7 can be found by

    tk7=Wk7,j×TMS7

    As mentioned above, the TMS of remaining rely in the path can be obtained by repeating Eqs 15–18.

    After designing the relay parameters of the main feeder, the relay TMS of branch lines that branched from the main feeder can estimated next. Let us consider the continuation of the tracking path in section 4.2.1. and assume that the set of branch path is { y w }. This path means relay w is located on the branch line that is branching from main feeder after relay y. The TMS and operation time of relay w can be appraised by using Eqs 15 and 18 respectively, as long as the subscript y change to w and the upper and lower bound should be estimated as

    UPBw=tkyCTIminWkw,j (19)
    LOWBw=tkyCTImaxWkw,j (20)

    After all the TMS setting values of all relays in the system be estimated, they can be used as a combination of possible solutions (a particle) and used as the basis for calculating the next-generation solution in PSO and MPSO-ASSOS.

    The fitness score of each particle is obtained by calculating the objective function in Eq 1 and considering the inequivalent constraints mentioned in Eq 4 and Eq 5. If one or more variables of the particle violate their limits, the corresponding particle will be punished by multiply a punishing factor leading to a lower fitness value. The fitness function of the coordination problem is defined as

    Fitness=1ρOBJF+εPUNFTMS+μPUNFCTI (21)
    PUNFTMSi={0,TMSminTMSkiTMSmaxP1,otherwise} (22)
    PUNFCTIi={0,CTIminCTIki,jCTImaxP2,otherwise}, (23)

    where

    PUNFTMS : The total summation of TMS punishing factor of all relays of each particle.

    PUNFCTI : The total summation of CTI punishing factor of all relays of each particle.

    ρ,ε,μ : Proportional constants.

    PUNFTMSi : TMS punishing factor of relay i.

    PUNFCTIi : CTI punishing factor of relay i.

    P1 : TMS punishing value of relay i is set to more than 103 when the TMS value is against the boundary.

    P2 : CTI punishing value of relay i is also set to more than 103.

    The proposed MPSO-ASSOS algorithm is simulated on 16-bus radial distribution systems. The proposed MPSO-ASSOS algorithm was applied to determine the optimal protection coordination issue on a 16-bus radial distribution system with DGs. The test system structure is shown in Figure 11; it has one main transformer, sixteen load buses and five DGs integration. The main substation consists of a 161 kV/22.8 kV, 50 MVA, Wye-Delta connection transformer. The DGs and load data are demonstrated in Table 4 and Table 5, respectively. The optimal protection coordination problem solved by MPSO-ASSOS was coded by MATLAB 2023 software. All programs were executed on a personal computer with Intel Core i7-10510U 2.3 GHz CPU and 16.0 GB RAM.

    Figure 11.  The topology of a 16-bus radial distribution system.
    Table 4.  Installed DG data of the system.
    DG Installed location bus Installed capacity (MW) Generation power (MW)
    1 8 10 6
    2 15 10 5
    3 11 10 4
    4 12 10 5
    5 5 10 5

     | Show Table
    DownLoad: CSV
    Table 5.  Lumped load data of the system.
    Bus Percentage of load type Lumped load
    Constant kVA (%) Constant Z (%) Active power (kW) Reactive power (kVAR)
    3 100 0 2244 1391
    6 100 0 2652 1130
    8 0 100 2412 1168
    9 0 100 1680 300
    11 100 0 1870 1159
    12 100 0 2771 1717
    13 0 100 2295 1422
    14 100 0 2592 1607
    15 0 100 3395 1446
    16 100 0 2125 1317
    Total 24036 12657

     | Show Table
    DownLoad: CSV

    The corresponding fault current, which flow through the specified relay under different fault, can be calculated from the ETAP software. In the simulation test of this paper, when no DG is integrated into the system, the ETAP is first used to calculate the fault current value flowing through each relay under conditions of different fault locations as the initial value of the relay's TMS setting. This test system will consider five DG-integrated operation scenarios, so there are eight operating conditions, including the above-mentioned initial operating states, as shown in Table 6. The fault current variation when the fault occurred at Bus 8, Bus 10 and Bus 14 for each DG operation condition are demonstrated in Figures 12, 13 and 14, respectively.

    Table 6.  Possible operation statuses of DGs under study.
    Statues (Condition) DGs ON DGs OFF DGs power injection (MW)
    1 None 1–5 0
    2 1–5 None 25
    3 2–5 1 19
    4 1, 2, 3, 5 4 20
    5 2–4 1, 5 19
    6 1, 3, 5 2, 4 15
    7 1, 3 2, 4, 5 10
    8 2, 5 1, 3, 4 10

     | Show Table
    DownLoad: CSV
    Figure 12.  Fault current passing through each relay when the fault occurred at Bus 8.
    Figure 13.  Fault current passing through each relay when the fault occurred at Bus 10.
    Figure 14.  Fault current passing through each relay when the fault occurred at Bus 14.

    According to the structure of the test distribution system shown in Figure 11, the coordination tracking path can be investigated and arranged as shown in Table 7. The Table lists 15 tracking paths. In addition to being used as the initial TMS solutions for estimating each relay, any two adjacent relays listed in each path also describe the connectivity relationship between the primary and backup relay. Therefore, the difference in operating time between the two relays is the CTI value.

    Table 7.  Coordination tracking path of the test system.
    Path List of relay No. (in order of circuit connectivity)
    1 1, 2
    2 1, 2, 3
    3 1, 2, 3, 4
    4 1, 2, 3, 4, 5
    5 1, 2, 3, 4, 5, 6
    6 1, 2, 3, 4, 5, 6, 7
    7 1, 2, 3, 4, 5, 6, 7, 8
    8 1, 2, 3, 4, 5, 6, 7, 9
    9 4, 10
    10 4, 10, 11
    11 4, 10, 12
    12 5, 13
    13 5, 13, 14
    14 14, 15
    15 14, 16

     | Show Table
    DownLoad: CSV

    The proposed MPSO-ASSOS and PSO original models will be used to solve the protection coordinating optimization problem considering the parallel operation of DGs. These two approaches were applied to investigate the required adjustments to the TMS parameters of existing protection relays with the same 16-bus MV system. In this section, the total relay operation time mentioned in Eq 1 is the objective function that calculates the total operation time summation according to the list of relay No. as shown in Table 7. Table 8 shows the total operation time of relays in the coordination tracking paths, which are evaluated from the optimal TMS setting value optimized by PSO and MPSO-ASSOS under different numbers of DGs integrated conditions, as shown in Table 6. It can be found from the data that as the number of DGs incorporated increases, the distribution of fault current changes. The reverse fault current causes the original fault current to decrease. The TMS of the power relay must be revised downwards, which shortens the operation time of the power relay so that the total operation time is reduced. Two heuristic search methods are tested using the same fitness function and constraint sets. Then, the convergence solutions were shown from the average value of each method after executing the corresponding program 30 times, with each time having 160,000 particles for 100 iterations. It is obvious that the proposed MPSO-ASSOS possesses a better search ability to obtain the suitable TMS parameter combination for the shorter total relay operating time than the PSO.

    Table 8.  Total relay operating time comparisons under different DG quantities.
    Quantity of DGs integrated Total DG injection power (MW) PSO (sec.) MPSO-ASSOS (sec.)
    0 0 30.75 30.31
    2 10 29.43 28.65
    3 15/19 33.15/32.95 29.65/29.10
    4 19/20 29.32/29.51 28.99/29.10
    5 25 29.57 28.74

     | Show Table
    DownLoad: CSV

    In this study, all the relays are assumed to be numerical OCR, which considers continuous values of TMS and PCS. As mentioned earlier, the PCS of each relay was prespecified in the work. The optimum TMS settings of OCRs solved by the proposed MPSO-ASSOS are used to protect the network in each given condition, as shown in Table 9. These settings can provide proper decisions for TMS setting of protection coordination in any possible DGs integrated conditions.

    Table 9.  The optimum TMS settings of OCRs are solved by proposed by MPSO-ASSOS.
    Relay Condition
    1 2 3 5 8
    1 0.1098 0.3542 0.3219 0.3116 0.2669
    2 0.5473 0.5414 0.5326 0.5708 0.5774
    3 0.5856 0.5636 0.5601 0.6099 0.6244
    4 0.4797 0.4539 0.4539 0.4965 0.5070
    5 0.3751 0.3736 0.3769 0.4206 0.4167
    6 0.2752 0.2802 0.2801 0.3181 0.3165
    7 0.1918 0.1815 0.1815 0.2222 0.2222
    8 0.0990 0.0817 0.0731 0.1081 0.1222
    9 0.0954 0.0729 0.0703 0.1080 0.1272
    10 0.3771 0.3726 0.3721 0.3538 0.4069
    11 0.2697 0.2670 0.2676 0.2624 0.3086
    12 0.2755 0.2696 0.2727 0.2603 0.3126
    13 0.2673 0.2804 0.2857 0.2776 0.2771
    14 0.1652 0.1664 0.1771 0.1716 0.1690
    15 0.0667 0.0582 0.0767 0.0755 0.0694
    16 0.0658 0.0591 0.0763 0.0739 0.0768

     | Show Table
    DownLoad: CSV

    The operating times of primary-backup relay pairs with their CTI are shown in Figures 15, 16, 17, 18 and 19 for different DG quantities mentioned in Table 6. Similar results of other DG combinations are not demonstrated because of the brevity. The symbols, tb and tp, inside the figure are the operation time of the backup relay and primary relay, respectively. CTI can be evaluated by subtracting tb by tp. It is to be noted that relay pairs 9 and 12 present the miscoordination situation in Figures 18 and 19 because of the high value of fault current flow through the backup relay, resulting in the CTI greater than 0.24 except for the above situations, and the remainder of these figure show that the CTI is between 0.2 and 0.24 sec. for all the relay pairs.

    Figure 15.  Operation times of the primary-backup relay pairs and the CTIs for Condition 1 with no DG incorporation.
    Figure 16.  Operation times of the primary-backup relay pairs and the CTIs for Condition 2 with all DG incorporation.
    Figure 17.  Operation times of the primary-backup relay pairs and the CTIs for Condition 3 with four DGs integration.
    Figure 18.  Operation times of the primary-backup relay pairs and the CTIs for Condition 5 with three DGs integration.
    Figure 19.  Operation times of the primary-backup relay pairs and the CTIs for Condition 8 with two DGs integration.

    The convergence of the MPSO-ASSOS, PSO and GA of searching for the best solution for optimal TMS setting is presented in Figure 20.

    Figure 20.  Convergence comparison in 100 generations(iterations) for Condition 2 with five DGs integration.

    These metaheuristic search methods were tested based on the same fitness function of Eq 21 and constraint set, and then the convergence diagrams were drawn from the average value of each method after executing the corresponding program 50 times, with each time having total 500 chromosomes or particles populations for 100 iterations. Figure 20 demonstrates that the proposed MPSO-ASSOS can obtain the best and stable solution about in the 36th generation (or iteration) and that it performs better than the other two algorithms in finding the maximum fitness value. Moreover, it obtains the best TMS combination of each relay and the approximate lowest total accumulated operation time.

    This article proposes an approach to the adaptive protection coordination problem considering the MV distribution system with high DG incorporation. The protection coordination model was formulated as the optimization problem model. This paper included two main contributions: First, a concept of tracking paths is proposed, the distribution system has divided into main feeder path and branch line path. Each tracking path set contains the relay number installed on the transmission line topology. Then, the TMS of each relay can be used to calculate the OT of the relay to determine the fitness of the TMS setting combination. Second, an MPSO- ASSOS algorithm with adaptive self-cognition and society operation scheme was proposed, and it is utilized to determine TMS for each relay on the tracking path subject to CTI and TMS limitation. A protection coordination problem of a 16-bus MV system with five DGs and sixteen OCRs is applied to test the proposed method. The simulation results show that the proposed MPSO-ASSOS reduces the overall OT and overcomes the impact on protection coordination settings of the system with different conditions of the DGs integration.

    Although the proposed algorithm can successfully handle the problem of relay protection coordination, according to the test of this article, there may still be situations in the system that cannot be coordinated due to the inherent design nature. This study recommends replacing all the existing electromechanical OCRs with modern numerical types. A central protection agent can apply the optimum TMS setting of the numerical OCRs to coordinate relays through telecommunication or wireless technology links.

    The authors declare they have not used Artificial Intelligence (AI) tools in creating this article.

    The author would like to thank the Financial support given to this work by the National Atomic Research Institute, Taiwan, R.O.C. under contract number NL1120262 is appreciated.

    The authors declare no conflict of interest.

    T.-S.Z. performed system model, data analysis, model simulations, article writing, and editing. Y.-D. L. and J.-L. Jiang provided the research funding support and related experiences, and C.-L. Su provided project ideas, conceptualization, and related experiences. J.-T. Yu performed the data curation and also conducted simulations. All authors have read and agreed to the published version of the manuscript.

  • This article has been cited by:

    1. Tung-Sheng Zhan, Yih-Der Lee, Jheng-Lun Jiang, Optimal OCR coordination in a high penetration distribution power system using a refined immune algorithm with an auto-tuning reproductive mechanism, 2024, 12, 2333-8334, 1225, 10.3934/energy.2024056
  • Reader Comments
  • © 2009 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3123) PDF downloads(645) Cited by(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog