To calculate the risk of an electrical fire resulting from violations of electrical codes and related risk factors, a specialist system as suggested in [6], was applied, considering the electrical infrastructure of the National Museum. This methodology offered a thorough risk evaluation in the period before the unfortunate fire.
The chosen methodology for risk indexing employs a fuzzy Petri net to represent the sequence of states and transitions associated with electrical failure mechanisms capable of causing fires. Fuzzy Petri nets are used in the construction of expert systems due to their ability to represent reasoning and knowledge through fuzzy production rules [23].
Transitions, on the other hand, represent events that could cause a change in the state of the electrical system. They may denote occurrences such as the overheating of a component due to an excessive electrical load or the onset of a short-circuit due to worn-out insulation. Arcs illustrate the causal relationships between places and transitions, providing a roadmap for the potential sequence of events that could lead to an electrical fire. This comprehensive representation helps in formulating a clear picture of the complex interplay of factors that may lead to potential electrical fires.
The fuzzy transition rules encapsulate the dynamics and transition logic of a fuzzy Petri network. There are essentially three types of rule productions with which any complex system can be represented: Type 1 (Simple); Type 2 (Logical AND) and Type 3 (Logical OR) [25].
To represent these relationships, FPRs use fuzzy operators, which allow for the representation of linguistic (fuzzy) variables. The fuzzy simulation theory adopts ‘‘MIN” operator to manage ‘‘AND” problems and ‘‘Max” operator to manage ‘‘OR” problems. [25], as shown Table 1.
The minimum operator is used to represent the logical AND operator in fuzzy logic. It is used to determine the degree of support for a rule based on the weakest premise. In FPRs, the minimum operator is used to calculate the degree of membership of a fuzzy set. The maximum operator is used to represent the logical OR operator in fuzzy logic. It is used to determine the degree of support for a rule based on the strongest premise. In FPRs, the maximum operator is used to calculate the degree of membership of a fuzzy set [26].
Together, these rules allow for the comprehensive representation and simulation of intricate complex systems within a fuzzy Petri network. This approach, using expert systems and fuzzy Petri nets, provides a comprehensive understanding of the potential chain of events leading to a fire. It thereby helps identify potential weaknesses and risks in the electrical system, allowing for a prioritization of mitigation efforts based on the assessed risks. The complete failure chain is shown in Fig. 21.
3.1. Matrix representation
The matrix representation of WFPNs provides an effective way to perform computational simulations of the system. By representing the system as a two-dimensional matrix, it is possible to perform various operations on the system, such as fire sequence determination of transitions and the behavior analysis of the system under different conditions.
One advantage of the matrix representation is that it allows for efficient computation of the system's behavior. The matrix can be updated based on the firing of transitions, allowing for real-time simulations, leading for rapid prototyping and testing of complex systems before their implemented in the real-life [27].
In addition, the matrix representation of WFPNs allows for the use of various mathematical and computational tools to analyze the system. For example: it is possible to use linear algebra techniques to analyze the matrix and determine the steady-state behavior of the system.
One of these approaches is presented in Liu H. et al. [28], where the authors introduce a matrix representation that enables computational simulation of complex systems. The Dynamic Adaptative Fuzzy Petri net (DAFPN) can be defined as an 11-tuple [26]
\(\:\text{D}\text{A}\text{F}\text{P}\text{N}=\:\left(\text{P};\:\text{T};\:\text{I};\:\text{O};\:\text{D};\:{\alpha\:};\:{\beta\:};\:\text{W};\:\text{U};\:\text{T}\text{h};\:\text{M}\right)\) , where:
\(\:\text{P}=\left\{{\text{p}}_{1},\:{\text{p}}_{2},\dots\:,\:{\text{p}}_{\text{m}}\right\}\) denotes a finite nonempty set of places; \(\:\text{T}=\left\{{\text{t}}_{1},\:{\text{t}}_{2},\dots\:,\:{\text{t}}_{\text{n}}\right\}\) denotes a finite nonempty set of transitions.
\(\:\text{I}\::\text{P}\times\:\text{T}\to\:\left\{0,\:1\right\}\) is an \(\:\text{m}\:\times\:\text{n}\) input incidence matrix defining the directed arcs from place to transitions.
\(\:{\text{I}}_{\text{i}\text{j}}=1\) , if there is a directed arc from pi to tj, and \(\:{\text{I}}_{\text{i}\text{j}}=0\), Iij if there is no directed arcs from pi to \(\:{\text{t}}_{\text{j}}\), for pi to \(\:\text{i}=1,\:2,\:\dots\:,\:\text{m}\) and \(\:\text{j}=1,\:2,\:\dots\:,\:\text{n}\).
\(\:\text{O}\::\text{T}\times\:\text{P}\to\:0,\:1\) is an \(\:\text{m}\:\times\:\text{n}\) output incidence matrix defining the directed arcs from transitions to places. \(\:{\text{O}}_{\text{i}\text{j}}=1\), if there is a directed arc from pi to \(\:{\text{t}}_{\text{j}}\), and \(\:{\text{O}}_{\text{i}\text{j}}=0\) if there is no directed arcs from pi to \(\:{\text{t}}_{\text{j}}\), for \(\:\text{i}=1,\:2,\:\dots\:,\:\text{m}\) and \(\:\text{j}=1,\:2,\:\dots\:,\:\text{n}\), \(\:\text{D}=\left\{{\text{d}}_{1},\:{\:\text{d}}_{2},\dots\:,\:{\text{d}}_{\text{m}}\right\}\) denotes a finite set of propositions.
\(\:\text{P}\cap\:\text{T}\cap\:\text{D}=\:{\varnothing}\) , \(\:\left|\text{P}\right|=\left|\text{D}\right|\). \(\:{\alpha\:}\::\text{P}\to\:\left[0,\:1\right]\) is an association function which maps from places to real values between 0 and 1.
\(\:{\beta\:}\::\text{P}\to\:\text{D}\) is an association function representing a bijective mapping from places to propositions.
\(\:\text{W}\::\text{I}\to\:\left[0,\:1\right]\) is an input function and it can be expressed as a \(\:\text{m}\:\times\:\text{n}\)-dimensional matrix. The value of an element in \(\:\text{W}\), \(\:{\text{w}}_{\text{i}\text{j}}\in\:\left[\text{0,1}\right]\), is the weight of the input place, which indicates how much the place pi impacts its following transition \(\:{\text{t}}_{\text{j}}\) connected by \(\:{\text{I}}_{\text{i}\text{j}}\).
\(\:\text{U}\::\text{O}\to\:\left[0,\:1\right]\) is an output function and can be expressed as a \(\:\text{m}\:\times\:\text{n}\)-dimensional matrix. The value of an element in \(\:\text{U}\), \(\:{{\mu\:}}_{\text{i}\text{j}}\in\:\left[\text{0,1}\right]\), is the value of certainty factor, which indicates how much a transition tj impacts its output places \(\:{\text{p}}_{\text{i}}\), if the transition fires.
\(\:\text{T}\text{h}\::\text{O}\to\:\left[\text{0,1}\right]\) is an output function which assigns a certainty value between 0 and 1 to each output place of a transition, \(\:\text{T}\text{h}=\:{\left({{\tau\:}}_{\text{i}\text{j}}\right)}_{\text{m}\text{x}\text{n}},\:\text{i}=1,\:2,\:...,\:\text{m};\text{j}=1,\:2,\:\dots\:,\:\text{n}\), denoting the output threshold of this place. \(\:{{\tau\:}}_{\text{i}\text{j}}\in\:\left[\text{0,1}\right]\), if there is a direct arc from \(\:{\text{t}}_{\text{j}}\) to pi, and \(\:{{\tau\:}}_{\text{i}\text{j}}=+{\infty\:}\), if there are no direct arcs from \(\:{\text{t}}_{\text{j}}\) to pi, for \(\:\text{i}=1,\:2,\:...,\:\text{m}\:\text{a}\text{n}\text{d}\:\text{j}=1,\:2,\:...,\:\text{n}\). \(\:\text{M}\) denotes a marking of the Petri Net, \(\:\text{M}=\:{\left({\alpha\:}\left({\text{p}}_{1}\right),\:{\alpha\:}\left({\text{p}}_{2}\right),\:\dots\:,\:{\alpha\:}\left({\text{p}}_{\text{m}}\right),\right)}^{\text{T}}\:\), where \(\:{\alpha\:}\left({\text{p}}_{1}\right)\) is the truth value of place pj. The initial marking is denoted by \(\:{\text{M}}_{0}\). Its determination should be based on an actual status. Hence, its value could be understood as the dynamic input and directly influence the dynamic behavior of DAFPN.
The evolution of failures chain depends on the execution of production rules, related to a fire of transitions mechanisms.
The notation \(\:\text{I}\left(\text{t}\right)=\left\{{\text{p}}_{\text{I}1},\:{\text{p}}_{\text{I}2},\dots\:,{\text{p}}_{\text{m}},\:\right\}\) represents the input with corresponding weights \(\:{\text{w}}_{\text{I}1}\), \(\:{\text{w}}_{\text{I}2},\dots\:,\:{\text{w}}_{\text{I}\text{m}}\), Similarly, \(\:\text{O}\left(\text{t}\right)=\left\{{\text{p}}_{\text{O}1},\:{\text{p}}_{\text{O}2},\dots\:,{\text{p}}_{\text{n}},\:\right\}\) represents the output with corresponding output thresholds \(\:{{\tau\:}}_{\text{O}1},\:{{\tau\:}}_{\text{O}2},\:\dots\:,\:{{\tau\:}}_{\text{O}\text{n}}\) and certainty factors \(\:{{\mu\:}}_{\text{O}1},\:{{\mu\:}}_{\text{O}2},\:\dots\:,\:{{\mu\:}}_{\text{O}\text{n}}\). The enabling and firing rules are specified as follows.
Enabling rule: \(\:\forall\:\text{t}\in\:\text{T}\), \(\:\text{t}\) is enabled and fired if \(\:\forall\:{\text{p}}_{\text{I}\text{j}}\in\:\text{I}\left(\text{t}\right)\),
$$\:\left\{{\alpha\:}\left({\text{p}}_{\text{i}\text{j}}\right)>0\right\}\wedge\:\left\{{\mu\:}\left(\text{t}\right)\ge\:\text{min}\left({{\tau\:}}_{\text{O}\text{K}}\right)\right\},\:\text{j}=1,\:2,\:\dots\:,\:\text{m};\text{k}=1,\:2,\:\dots\:,\:\text{n}.$$
1
Where \(\:{\alpha\:}\left({\text{p}}_{\text{i}\text{j}}\right)\), \(\:\:{\alpha\:}\left({\text{p}}_{\text{i}\text{j}}\right)\in\:\left[\text{0,1}\right]\), is the fuzzy truth value in place \(\:{\text{p}}_{\text{I}\text{j}},\:{\text{p}}_{\text{I}\text{j}}\in\:\text{I}\left(\text{t}\right)\:\), wich indicates the truth degree of proposition \(\:{\text{d}}_{\text{I}\text{j}},\) if \(\:{\beta\:}\left({\text{p}}_{\text{I}\text{j}}\right)=\:{\text{d}}_{\text{I}\text{j}};\:{\mu\:}\left(\text{t}\right)=\:\sum\:_{\text{j}=1}^{\text{m}}{\alpha\:}\left({\text{p}}_{\text{i}\text{j}}\right){\text{w}}_{\text{I}\text{j}},\:\text{j}=1,\:2,\:\dots\:,\:\text{m}\) is the equivalente fuzzy truth of input places at transitions when \(\:\text{t}\) is enabled.
After \(\:\text{t}\) is fired, the tokens in input places are copied, and tokens with fuzzy truth are put into each output place whose output thresholds are lesser than the equivalent fuzzy truth. The new fuzzy truth values of output places are defined as:
$$\:{\alpha\:}\left({\text{p}}_{\text{O}\text{i}}\right)=\left\{\begin{array}{c}{{\mu\:}}_{\text{O}\text{i}}\mu\:\left(\text{t}\right),\:\mu\:\left(t\right)\ge\:{{\tau\:}}_{\text{O}\text{i}}\\\:0,\:\:\:\:\:\:\:\:\:\:\:\:\:\:\mu\:\left(t\right)<{{\tau\:}}_{\text{O}\text{i}}\:\end{array}\right.\text{i}=1,\:2,\:\dots\:,\:\text{n}.$$
2
If a place has more than one input transition and more than one of its input transitions fires, then the transition produces the new fuzzy truth value of the output place with the maximum fuzzy truth.
3.3. Model parametrization
This comprehensive fuzzy Petri net simulation provides a probabilistic risk index for various electrical fault events that can lead to a fire. Each proposition represents events or states, which in this context can be physical and chemical mechanisms, risk factors, electrical malfunctions, and fires. The value assigned to each proposition represents the likelihood of occurrence of each of these events or states.
The model allows for prioritizing corrections based on the obtained results. The ability to discern which factors pose the greatest risk and are more readily correctable can guide effective intervention strategies.
3.3.1. General fire risk indexing before the fire
Figure 21 illustrates the influence of each risk factor on others, specifically concerning the mechanisms involved in electrical fires. Based on fuzzy arithmetic, events with feedback don't result in an "overflow" of the proposition's value. This facilitates an effective depiction of the intricate interplay between risk factors.
The circles in black represent the issues within the museum's electrical installation prior to the fire, in accordance with the investigation report, serving as inputs to the developed mathematical model. The green color circles represent events/states with a zero probability, i.e., implying unlikely. The color scale, ranging from yellow to red, denotes the occurrence possibility for each event, where yellow and red suggest the lowest and highest likelihood respectively. The onset risk of a fire is signified by the proposition p60, displayed in orange, highlighted by the black arrow in Fig. 19. This signifies a tangible fire possibility, given the reported conditions and characteristics of the electrical installation.
This outcome clearly signifies the urgent need for corrective interventions in the installation, which were unfortunately neglected, potentially resulting in the fire incident. Despite financial constraints cited by the museum staff, it seems there was a misjudgment in evaluating the severity of existing issues, given that addressing many of the non-conformities would not necessarily have required substantial investment.
3.3.2. Hazard ranking through the application of the methodology
This study enables the hazardous hierarchization by applying the proposed methodology, simulating the individual impact of each risk factor on the probability of fire ignition (p60) as depicted in Fig. 19. The hierarchy is derived from calculating the value of each proposition through computational simulation, representing the likelihood of occurrence for each event. This method provides a more strategic approach to risk management. By identifying, quantifying, and prioritizing risks, it aids in making better-informed and efficient decisions, ultimately enhancing the overall safety of electrical installations. Considering a range from 0 to 1, Table 2 showcases the impact on safety for each electrical issue found in the National Museum of Rio de Janeiro, Brazil.
Table 3
Impact on the safety of each electrical issues.
Relevance | Place | Electrical issues found | Likelihood of fire p(60) |
1 | p32 | Are there signs of oxidation on electrical connections? | 0.500 |
2 | p22 | Are there loose connections? | 0.375 |
3 | p03 | Are there inconsistencies in the design of circuit breakers and/or RCCB? | 0.281 |
4 | p29 | Is there excessive dust or moisture on electrical connections or terminals? | 0.250 |
5 | p04 | Are there reports of frequent breaker tripping? | 0.211 |
6 | p20 | Are there poorly executed cable splices, lack of signaling? | 0.141 |
7 | p23 | Are there connections outside junction boxes or sealing? | 0.125 |
8 | p45 | Are there issues with the grounding and bonding system? | 0.125 |
9 | p01 | Has the installation more than 30 years old? | 0.094 |
10 | p02 | Are there lack of periodic maintenance and inspections? | 0.094 |
11 | p24 | Are there conductors exposed to mechanical damage? | 0.094 |
12 | p36 | Are there wires exposed to sun showing color fading? | 0.079 |
13 | p43 | Are there signs of rodent existence? | 0.047 |
14 | p46 | Is there no SPD in the main electrical panel? | 0.031 |
15 | p34 | Are there non-protected PVC cables exposed to direct sunlight? | 0.020 |
Based on the data in Table 2, the following rationale can be provided for the varying criticality of different failures:
-
Signs of Oxidation on Electrical Connections (weight = 0.500): Oxidation increases the resistance at electrical connections. This, in turn, can lead to excessive heating during current flow, which has the potential to ignite nearby materials, hence the high fire risk.
-
Loose Connections in Outlets, Switches, Breaker Terminals (weight = 0.375): Loose connections can cause electrical arcing. This phenomenon generates a significant amount of heat that can potentially start a fire, explaining its high fire risk.
-
Inconsistencies in Circuit Breakers/RCCB Design (weight = 0.281): Circuit breakers and RCCBs are crucial protective devices in an electrical system. Any inconsistency in their design can prevent them from functioning properly during overloads or short-circuits, thereby escalating the risk of an electrical fire.
-
Excessive Dust or Moisture on Electrical Connections/Terminals (weight = 0.250): Dust and moisture could lead to short-circuits or provide fuel for a potential fire. Moisture can also cause corrosion, which increases resistance and generates heat.
-
Frequent Breaker Tripping (weight = 0.211): This often indicates an underlying issue in the electrical system, such as overloading or short-circuits, both of which can cause excessive heat or sparks and potentially lead to fires.
-
Poorly Executed Cable Splices, Lack of Signaling (weight = 0.141): Poorly executed splices increase resistance, generate heat, and can also cause arcing. The lack of signaling may lead to using the system in unsafe conditions, unknowingly.
-
Connections Outside Junction Boxes or Sealing Problems in Electrical Cabinets (weight = 0.125): These issues expose electrical connections to environmental factors that could lead to short circuits or fires.
The other risk factors from the middle to the bottom of the list, such as the existence of rodents or non-protected PVC cables exposed to sunlight, while still posing a risk, are less direct or immediate in their potential to cause a fire. Rodents could cause damage over time, and unprotected cables could degrade under sunlight, but these risks are slower acting and thus less likely to cause a fire within a short time frame.
The prioritization of electrical issues can serve as an indicative guide for initiating corrections. However, the individual contribution of each risk factor may vary due to its interaction with other present factors, akin to what occurs in real life scenarios. While an issue might initially seem less critical, its significance can increase substantially when combined with other factors within a particular environment or situation. Consequently, understanding the interplay between these elements is just as crucial for effective and comprehensive risk management as identifying the individual issues. This complexity mirrors real-life situations where multiple variables coexist and interact in ways that can escalate or mitigate potential risks.