# Evolutionary Stability Analysis of the Coal Mine Safety Management System Governed by Delay and Impulsive Differential Equations

Authors: Li Bai; Xinhua Wang; Rongwu Lu; Hao Yu; Dan Li
Journal: Complexity
 (2020)

---

## Abstract

Control of unsafe behaviors has always been a difficult point in Coal Mine Safety Management (CMSM). To study the regular pattern of behavior propagations in CMSM, a series of evolutionary dynamical system models are established by the delay and impulsive differential equations (DIDEs). The threshold constraints (TCs) and impulse responses in the models can accurately correspond to the actual inspection and supervision measures in safety management. Through the simulation analysis and comparison of the system stability, it can be concluded that some proper settings of the TCs and impulse responses bring satisfactory evolutionary stability. Based on the main results of this study, we put forward some decision-making suggestions on the control points and control strategies of unsafe behaviors in safety management.


---

## Body

## 1. Introduction

Coal mine safety responsibility accidents are the results of the continuous accumulation of unsafe behaviors, and the spread of unsafe behaviors is the main cause of the accumulation of unsafe behaviors [1, 2]. But the occurrence of coal mine safety accidents is uncertain at any moment, and employees have differences in risk perception and preferences, cognitive abilities, etc. Therefore, although some employees know that unsafe behaviors (illegal operations) may cause safety accidents, they might still choose unsafe behavior and become the implementer of unsafe behaviors under the influence of such factors as production tasks, extra payoffs, and trust to luck. China’s coal mine production is mostly arranged in district groups. There are a large amount of contacts between employees, which can communicate and transmit various information and then affect each other’s behaviors. The spread of behavior is not limited to the field of social bad behavior, but also spreads in other fields of behavior. Many scholars have also used virus propagation models or system dynamics theories and methods to study the laws of behavioral propagation from different aspects [3–5].The policies of safety management in nowadays China are always adjusted with the occurrence of accidents. The relevant supervision measures are characterized by accident types and lack of comprehensive analysis and cannot take certain effects in the long run. Actually, CMSM utilizes a combination of economic, institutional, and cultural methods to control unsafe behaviors in coal mining. Its goal is to restrain the occurrence frequencies of the unsafe behaviors and improve the coal mine safety. As a controlling method of unsafe behaviors, CMSM has already attracted many scholars’ attention [6–10]. Based on the theoretical analysis of the evolutionary stability strategy, we have already put forward some incentive measures and restraint measures to control the instability of player’s dynamical behaviors in CMSM [11]. It is found that game strategy selection is controllable under the corresponding measures, and appropriate incentive and constraint measures can mobilize the initiative of safe behaviors. However, those flexible incentive and restraint measures for the excess subsidies and severe punishments may bring some unexpected economic or social problems, which make some small and medium coal mine enterprises unable to implement those measures [12–14].The delay and impulsive differential equations (DIDEs) have been widely used in biological control, artificial intelligence, and other fields [15–22]. There are a large number of systems in various engineering fields that continuously change in a certain time interval. Due to some artificial intervention behaviors, the system state will be disturbed for a relatively short period of time, which will cause the system state to change instantly. This phenomenon is also common in CMSM. It actually takes a certain amount of time delays to check the safe behaviors and feedbacks on the inspection results, and safety management measures would be altered much more strictly when colliery safety entered a tough situation (such as a safety accident just happened in a coal mine and many other coal mines have similar production conditions with it), which corresponds exactly to the impulse response in the models. The dynamical system governed by the DIDEs is one of the right methods to analyze those sudden changes in the system state [23–27]. So in this paper, we want to set up some evolutionary game models with time delay and state feedback impulse to analyze the behavioral communications of those participants in CMSM. This microperspective analysis method can make an in-depth analysis of unsafe behavior’s communications, and we believe that the design of control mechanisms for unsafe behaviors is the key to promote the long-term effectiveness of safety measures in CMSM.The rest of this paper is organized as follows. In Section2, we build an evolutionary game model in CMSM by the DIDEs with threshold constraint (TC) on x1. Then, in Section 3, we carry out detailed and in-depth stability analysis of the DIDEs with TC on x1, which involves various possible pulse states. In Section 4, we carry out stability analysis of the DIDEs with TCs on x3 or x2. In Section 5, this paper is finally concluded and some suggestions are put forward to improve the evolutionary stability in CMSM.


## 2. Preliminaries

There are many factors that affect the game behaviors of each interested group [28]. If the employees are completely rational, the behavior selection process can be described by the expected utility function theory. But in reality the employees are not completely rational, and the behavior selection is mainly based on the psychological expectation and perceived value of the profit and loss value of the behavior choice. Behavioral communication between employees often occurs in every coal mine process. CMSM also involves many interest parties and there are complex relationships between them. Generally, there are at least three participants involved: coal miners (denoted as player I), safety inspection groups inside the coal mine (player II), and safety regulatory departments of local government (player III). In evolutionary game theory, the advantage of a population is reflected in the proportion of participants. A multiparty evolutionary game model has been established to analyze the selection strategies of safety behaviors in our previous work, and we have obtained the following replicator dynamics [29]:(1)x˙1=x11−x1x2+1−x2x3f1−c1,x˙2=x21−x21−x1f1+x3f2−c2,x˙3=x31−x31−x11−x2f1+f2−c3,where x1,x2, and x3 are the probabilities of those players choosing safety behaviors (safe operations, or safe inspections, or safe supervisions), respectively. cii=1,2,3 are, respectively, the safety costs (inputs) paid by player I, player II, and player III, and fii=1,2 are the penalty index when player I and player II make mistakes (unsafe behaviors) and are checked out, respectively. For the value of the parameter, we make an explanation here in particular. The values of the parameters are relative, i.e., the values of the parameters represent the costs and punishments paid by the participants for the benefit of a unit.According to the actual situation of China’s coal mine production and coal mine safety regulation, these parameters are assumed to meet the following conditions:(i)
c2≤c3≤c1≤1(ii)
fi≥1,i=1,2.When the implementation of the employees’ unsafe behaviors can bring them additional benefits, they may be imitated and learned by other employees so that the unsafe behaviors can be spread among the employees. Many operating activities in coal mine production are completed through the mutual cooperation of group members. Group members are regulated and restricted by group norms, thereby maintaining roughly uniform behavioral goals and cognitive structures. Some unsafe behaviors have a group nature, such as violations of regulation production. This kind of group unsafe behavior sometimes happens under the influence of higher-level illegal command, group norms, and other factors [30–33].As for the situation of CMSM, we can make further consideration on the participants’ group unsafe behaviors. If one conducts safe/unsafe behaviors in coal mine production, it always takes a period of time (inspection and supervision come one after another) to receive the consequences (wages, subsidies, or penalties). The characteristic of the delay differential equations is that the right end of the equation does not depend on the current moment, but depends on the function values of the past moment. Hence, we consider the following DIDEs:(2)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,x1>x1∗,Δx2t=α,α>0Δx3t=β,β>0x1≤x1∗,where τ1,τ2, and τ3 are, respectively, the time delays for x1,x2, and x3, which denote the times between players receiving payoffs and their behaviors been done. Because player I receives inspections and supervision from the other two players and even colleagues in the same shift, player I knows the consequences in the fastest time. Similarly, the frequency of player II’s safety inspections is faster than that of player III. So, it is assumed that(3)0<τ1≤τ2≤τ3.Coal mine workers are the most direct producers of unsafe behaviors, and they are also the ultimate managed parties in safety management. Ifx1 is found to be less than a threshold x1∗ (threshold is a certain number, which exists in some mature colliery enterprises or coal industry associations) after every fixed time T (pulse period, which is the time interval between two adjacent pulses in a periodic repeating pulse sequence), such impulses will be generated as Δx2t for x2 and Δx3t for x3, respectively. Therefore, the DIDEs (2) denote an evolutionary dynamical system with TC on x1 in CMSM.The exact solution of the evolutionary dynamical system (2) is undoubtedly difficult to obtain, so the model analysis process is mainly carried out through simulation analysis. The general process of numerical simulations is described briefly as follows:Step 1: the entire time period is divided intoN intervals of the same length tspan, and let i=1, the initial values of x1,x2,x3 are given x1,x2,x3t=0.Step 2: inith interval, the delay differential equations (DDEs) with constant delays (i.e., the first three equations in the DIDEs) are solved by MATLAB Function dde23 with given time delays and initial values of x1,x2,x3t=i−1, and some figures such as phase diagram and TSG can be drawn by plot3/plot function of MATLAB. Then, at the right end of the interval, the solutions are checked whether the pulse trigger condition (TC) is met. If TC is not satisfied, the value of the solution vector at the last moment is the initial value when DIDEs are solved in the next interval. Otherwise, the pulse value is added to the solution vector at the last moment to obtain the initial value of the next interval.Step 3:i=i+1, loop to step 2 until the DIDEs are solved and simulated in the entire time period.In the simulation process in this section, the model parameters are set as follows:c1=0.5;c2=0.2;c3=0.3;f1=1;f2=1.2;α=0.3;β=0.2. In order to illustrate the generality and consistency of this study, we first analyze the game dynamic system without time delay and impulse. Let τ1=τ2=τ3=0 and x1∗=1 in DIDEs (2). Figures 1(a) and 1(b) are time series graphs (TSGs) of players’ strategy selection probabilities x1,x2,x3 when the dynamical system (2) has no pulse and no delay. The initial values of x1,x2,x3 are P00.5,0.3,0.2 and P10.7,0.5,0.3, respectively.Figure 1
TSGs ofxi when the DIDEs (2) have no pulse and no delay. (a) x1,x2,x3t=0=P00.5,0.3,0.2 and (b) x1,x2,x3t=0=P10.7,0.5,0.3.
(a)(b)As can be seen in Figure1, there are repeated periodic fluctuations in the game process no matter what the initial value is. It illustrates that participants in the safety management system constantly change their choices based on the coal mine safety status and the behavioral selections of others. When player I performs well (x1 takes a large value), player II will relax its supervision (x2 decreases). Then, the value of x1 will gradually decrease, and then, x2 will increase, and then, x1 will increase again. The dynamical game system enters a cyclical process, where player III plays a smaller and smaller role (x3 gradually decreases towards zero). Once coal mine safety enters such a periodic fluctuation state, the efficiency of safety production management can hardly be guaranteed.Due to the sequential nature of production behaviors and supervision behaviors, the variablesx1,x2, and x3 certainly have time delays. However, the impulse response is usually artificially set, which means the supervisory behaviors suddenly change or become stronger in reality. So we should take pulse responses and TCs into account and perform detailed stability analysis on the dynamical system to demonstrate the influence of the setting of pulses and TCs.


## 3. Stability Analysis on the DIDEs with TC onx1 in CMSM

In this section, to reflect the correspondence and influence of time delays and pulses in the dynamical system (2), the evolutionary stability of CMSM is studied and simulated in a variety of cases. In the general theoretical study of the stability of dynamic systems, the stability of each variable needs to meet certain requirements. Since this study takes safety management as the application background and the main object of behavior control is player I, the stability analysis in this article is mainly on x1. The behavior control strategies of player II and player III are regarded as the adjustments of safety management decisions, and there are no special requirements for the stability of x2 or x3.During model simulation,lags=τ1,τ2,τ3 denotes the vector of the delay values of variables xi,i=1,2,3, which have been previously mentioned in Section 2. From now on, we set lags=3,5,7/30 in every process of the model simulations in this study, which are a bit small. It is based on the consideration that behavioral communication between employees in coal mining enterprises is fast and significant, where mostly a large number of workers are engaged in the same workflow and workplace, which nearly unchanged for a long time. Besides, the value of pulse period T=7, which never changes but in Section 4.3, where we will point out specially.


### 3.1. The DIDEs (2) Have Time Delays but No Pulse

As shown in Figure2, there are similar fluctuations with Figure 1 when the dynamical system (2) has time delays but no pulse. Although the fluctuation period of the variable x1 is gradually increasing, the fluctuation of x2 is becoming more severe. There is also no stability in safety management in this case, and the effectiveness of safety management will not be satisfactory.Figure 2
TSGs when the DIDEs (2) have delays but no pulse. Lags=3,5,7/30. (a) x1,x2,x3t=0=P00.5,0.3,0.2 and (b) x1,x2,x3t=0=P10.7,0.5,0.3.
(a)(b)


### 3.2. The DIDEs (2) Have Pulses

The implementers of the safety supervision measures are player II and player III, while player I is the object of its application. The pulse setting for player I has no practical meaning in CMSM. Therefore, we should not design any pulse forx1, but apply a fixed value pulse to x2 or x3. As mentioned earlier in Section 2, a pulse on x2 or x3 corresponds to the strengthening or adjustment of safety supervision measures. Specifically, if x3 has a pulse, it can correspond to strengthen the safety supervision of government authorities (such as the increase in the frequency of safety inspections, or the expansion of its scope). When x2 has a pulse, the safety inspection team of coal mine enterprises will strengthen their supervision level in proportion. As the dynamical system is found not very dependent on variables’ initial values from Figures 1 and 2, the initial value of x1,x2,x3 is given as P10.7,0.5,0.3 in this section.Figure3 shows TSGs of DIDEs (2) when x3 has a pulse separately α=0;β=0.2, where all the three variables x1,x2, and x3 involve in a cyclical fluctuation process. It illustrates that x3 becomes smaller and is even close to zero after a period of time, and then, the pulse takes its effect. In this case, player III actively participates in the safety supervision. As x1∗ changes its value from 0.5 to 0.6, the peaks of x2 and x3 are relatively reduced, but the fluctuation of x1 is more frequent.Figure 3
TSGs of DIDEs (2) when x3 has a pulse β=0.2. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Figure4 shows TSGs of DIDEs (2) when x2 has a pulse separately α=0.3;β=0. It shows x1 still has no stability, and x2 fluctuates more violently while x3 gradually decreases to zero. As x1∗ changes its value from 0.5 to 0.6, the fluctuation of x1 finally weakens while the peak of x2 intensifies. However, the system stability does not improve significantly.Figure 4
TSGs of DIDEs (2) when x2 has a pulse α=0.3. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Figure5 shows TSGs of DIDEs (2) when both x2 and x3 have pulses α=0.3;β=0.2. We can find that x1=1 lasts even longer and x2 and x3 are more actively involved in the dynamical system. And, the dynamical system (2) reduces volatility when x1∗ increases. However, all of these three cases shown in Figures 3–5 are still oscillating waveforms and there is no essential change from the previous case shown in Figure 2. It can be concluded from the simulation results in this section that there is almost no evolutionary stability in the dynamical system with TC on x1 (i.e., the DIDEs (2)).Figure 5
TSGs of DIDEs (2) when both x2 and x3 have pulses α=0.3;β=0.2. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Therefore, the probability of player I’s safety behavior selection should not be taken as a prerequisite condition for adopting safety supervision measures in the practice of safety management. Although coal miners are the ultimate object and the basic part for the implementation of safety management in CMSM, safety management of coal miners is the most difficult part due to the large number and the large-scale differences of skill and knowledge. Therefore, we analyze the evolutionary stability of the dynamical system where there are TCs imposed on other variables.


## 3.1. The DIDEs (2) Have Time Delays but No Pulse

As shown in Figure2, there are similar fluctuations with Figure 1 when the dynamical system (2) has time delays but no pulse. Although the fluctuation period of the variable x1 is gradually increasing, the fluctuation of x2 is becoming more severe. There is also no stability in safety management in this case, and the effectiveness of safety management will not be satisfactory.Figure 2
TSGs when the DIDEs (2) have delays but no pulse. Lags=3,5,7/30. (a) x1,x2,x3t=0=P00.5,0.3,0.2 and (b) x1,x2,x3t=0=P10.7,0.5,0.3.
(a)(b)


## 3.2. The DIDEs (2) Have Pulses

The implementers of the safety supervision measures are player II and player III, while player I is the object of its application. The pulse setting for player I has no practical meaning in CMSM. Therefore, we should not design any pulse forx1, but apply a fixed value pulse to x2 or x3. As mentioned earlier in Section 2, a pulse on x2 or x3 corresponds to the strengthening or adjustment of safety supervision measures. Specifically, if x3 has a pulse, it can correspond to strengthen the safety supervision of government authorities (such as the increase in the frequency of safety inspections, or the expansion of its scope). When x2 has a pulse, the safety inspection team of coal mine enterprises will strengthen their supervision level in proportion. As the dynamical system is found not very dependent on variables’ initial values from Figures 1 and 2, the initial value of x1,x2,x3 is given as P10.7,0.5,0.3 in this section.Figure3 shows TSGs of DIDEs (2) when x3 has a pulse separately α=0;β=0.2, where all the three variables x1,x2, and x3 involve in a cyclical fluctuation process. It illustrates that x3 becomes smaller and is even close to zero after a period of time, and then, the pulse takes its effect. In this case, player III actively participates in the safety supervision. As x1∗ changes its value from 0.5 to 0.6, the peaks of x2 and x3 are relatively reduced, but the fluctuation of x1 is more frequent.Figure 3
TSGs of DIDEs (2) when x3 has a pulse β=0.2. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Figure4 shows TSGs of DIDEs (2) when x2 has a pulse separately α=0.3;β=0. It shows x1 still has no stability, and x2 fluctuates more violently while x3 gradually decreases to zero. As x1∗ changes its value from 0.5 to 0.6, the fluctuation of x1 finally weakens while the peak of x2 intensifies. However, the system stability does not improve significantly.Figure 4
TSGs of DIDEs (2) when x2 has a pulse α=0.3. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Figure5 shows TSGs of DIDEs (2) when both x2 and x3 have pulses α=0.3;β=0.2. We can find that x1=1 lasts even longer and x2 and x3 are more actively involved in the dynamical system. And, the dynamical system (2) reduces volatility when x1∗ increases. However, all of these three cases shown in Figures 3–5 are still oscillating waveforms and there is no essential change from the previous case shown in Figure 2. It can be concluded from the simulation results in this section that there is almost no evolutionary stability in the dynamical system with TC on x1 (i.e., the DIDEs (2)).Figure 5
TSGs of DIDEs (2) when both x2 and x3 have pulses α=0.3;β=0.2. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x1∗=0.5 and (b) x1∗=0.6.
(a)(b)Therefore, the probability of player I’s safety behavior selection should not be taken as a prerequisite condition for adopting safety supervision measures in the practice of safety management. Although coal miners are the ultimate object and the basic part for the implementation of safety management in CMSM, safety management of coal miners is the most difficult part due to the large number and the large-scale differences of skill and knowledge. Therefore, we analyze the evolutionary stability of the dynamical system where there are TCs imposed on other variables.


## 4. Stability Analysis on the Dynamical System with TCs onx2 or x3

Compared with the coal miners, staff members in safety inspection groups of the coal mine enterprise, or in the safety supervision departments of the government, have higher knowledge and technical level and safety quality but less quantity. Their organizational integrity is much stronger, and the adjustment of their safety decisions should be more effective and faster. So in this section, we conduct stability analysis on the dynamical system when the TCs are imposed onx2 or x3.


### 4.1. Stability Analysis on the DIDEs with TC onx3

The safety supervision departments of the government (player III in the dynamical systems) play a leading role in the entire process of coal mine safety supervision [34–36]. As policy makers and direct implementers in CMSM, they are responsible for organizing, supervising, and supporting the relevant units in coal mine enterprises to perform the duties of safety inspections and timely coordinating and solving major problems in coal mine safety. The policy adjustment of the government safety departments undoubtedly can make an important impact in CMSM. Hence, we consider the following DIDEs with a TC on x3:(4)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,x3>x3∗,Δx2t=α,α≥0Δx3t=β,β≥0x3≤x3∗.To analyze the influence of the safety supervision departments adjusting their behaviors on the evolutionary stability of dynamic system (4), we consider just one case here, i.e., x3 individually imposed a pulse (β).As can be seen from Figure6, the dynamical system (4) has a satisfactory evolutionary stability when x3∗ is large enough. When x3∗ is not less than 0.3, x1 will steadily approaching 1. Even when x3∗=0.2, which is a bit small value, the volatility of the dynamical system becomes very weak and x1 gradually stabilizes to a relatively large value (see Figure 6(a)). So it is concluded that the adjustment of government departments’ safety measures plays an important role in improving the safety management level of coal mine.Figure 6
TSGs of DIDEs (4) with TC on x3 when x3 has a pulse α=0.3. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x3∗=0.2, (b) x3∗=0.3, (c) x3∗=0.5, and (d) x3∗=0.6.
(a)(b)(c)(d)


### 4.2. Stability Analysis on the DIDEs (4) with TC on x2

In accordance with the previous section, the effect of the decision adjustment of player II (safety inspection groups inside coal mine enterprises) on CMSM is studied in this section. When there is a threshold constraint (TC) onx2, we get the following DIDEs:(5)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,x2>x2∗,Δx2t=α,α≥0Δx3t=β,β≥0x2≤x2∗.Here, we also consider just one case, i.e.,x2 individually imposed a pulse (α), to analyze the evolutionary stability of the dynamic system (4). Some of the simulation results are shown in Figure 7.Figure 7
TSGs ofxi with TC on x2 when x2 has a pulse β=0.2. x1,x2,x3|t=0=P00.5,0.3,0.2. (a) x2∗=0.2, (b) x2∗=0.3, (c) x2∗=0.4, and (d) x2∗=0.5.
(a)(b)(c)(d)As can be seen from Figure7, the dynamical system (5) also has a satisfactory evolutionary stability when x2∗ is large enough (not less than 0.3). Therefore, it is also concluded that the decision adjustment of safety inspection groups can promote the safety management level of coal mining.


### 4.3. Stability Analysis on the Dynamical System with More TCs

We have studied the evolutionary stability of the dynamical system with TC imposed on the variablesx1,x2, and x3, respectively, in Section 3 and Sections 4.1 and 4.2. As a neutralization of the previous cases, we consider the following DIDEs, where there are more TCs in CMSM:(6)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,xi>xi∗,i=1,2,3,Δx2t=α,α>0Δx3t=0x1≤x1∗,x3>x3∗,Δx2t=α,α>0Δx3t=β,β>0x1≤x1∗,x3≤x3∗.To reflect the impact of TCs, we set pulse periodT=3 in the model simulations of the DIDES (6), and thus, the pulse check frequency is higher than before. By contrasting Figures 6–8, the results of the simulation show that the system stability in this case is no better than the previous two cases. The volatility of x1 may last forever (see Figure 8).Figure 8
TSGs of the DIDEs (5) with more TCs α=0.3,β=0.2. Pulse period T=3, x1,x2,x3|t=0=P00.5,0.3,0.2. (a) x1∗=0.6,x2∗=0.2,x3∗=0.1, (b) x1∗=0.8,x2∗=0.2,x3∗=0.1, (c) x1∗=0.9,x2∗=0.3,x3∗=0.1, and (d) x1∗=0.9,x2∗=0.2,x3∗=0.3.
(a)(b)(c)(d)But on the other hand, these results can also be acceptable because of the relatively large values ofx1, especially when the level of safety management technology is relatively low, which is still common in some remote areas or small- and medium-sized coal mining enterprises in China [37]. The satisfactory evolutionary stability is difficult to achieve in the practice of CMSM. As shown in Figures 8(c) and 8(d), once the pulse heuristic conditions (the TCs in the DIDEs (6)) are set properly, the dynamic system can make moderate fluctuations over a reasonable range, where the troughs of x1 are almost not less than the set value of x1∗. Meanwhile, all participants in CMSM can restrict and promote each other and actively engage in safety management. That is exactly in line with the aspirations of the Chinese government: to deal with safety management, all participators should take part in the whole process and manage it together.In this section, stability analyses are made on the dynamical system with different TCs. Table1 shows the simulation results briefly. When the TCs are imposed on x2 or x3 separately, x1 is evolutionary stable but one of the two supervisors performs inactively. In the last case, when more TCs are combined in the dynamical system, it is gratifying that x1 is progressively stable and both the other two have active performances. It is a good revelation for us. We should appropriately adjust the TCs and the thresholds of them, so that the management effect reaches the actual demand of CMSM in practice.Table 1
Evolutionary stabilities of the dynamical system with different TCs.
TCsStability ofx1Performance ofx2Performance ofx3Picturesx3≤x3∗=0.3Evolutionary stableInactiveActiveFigures6(b)–6(d)x2≤x2∗=0.3Evolutionary stableActiveInactiveFigures7(b)–7(d)More TCs combinedProgressively stableActiveActiveFigures8(c) and 8(d)


## 4.1. Stability Analysis on the DIDEs with TC onx3

The safety supervision departments of the government (player III in the dynamical systems) play a leading role in the entire process of coal mine safety supervision [34–36]. As policy makers and direct implementers in CMSM, they are responsible for organizing, supervising, and supporting the relevant units in coal mine enterprises to perform the duties of safety inspections and timely coordinating and solving major problems in coal mine safety. The policy adjustment of the government safety departments undoubtedly can make an important impact in CMSM. Hence, we consider the following DIDEs with a TC on x3:(4)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,x3>x3∗,Δx2t=α,α≥0Δx3t=β,β≥0x3≤x3∗.To analyze the influence of the safety supervision departments adjusting their behaviors on the evolutionary stability of dynamic system (4), we consider just one case here, i.e., x3 individually imposed a pulse (β).As can be seen from Figure6, the dynamical system (4) has a satisfactory evolutionary stability when x3∗ is large enough. When x3∗ is not less than 0.3, x1 will steadily approaching 1. Even when x3∗=0.2, which is a bit small value, the volatility of the dynamical system becomes very weak and x1 gradually stabilizes to a relatively large value (see Figure 6(a)). So it is concluded that the adjustment of government departments’ safety measures plays an important role in improving the safety management level of coal mine.Figure 6
TSGs of DIDEs (4) with TC on x3 when x3 has a pulse α=0.3. x1,x2,x3|t=0=P10.7,0.5,0.3. (a) x3∗=0.2, (b) x3∗=0.3, (c) x3∗=0.5, and (d) x3∗=0.6.
(a)(b)(c)(d)


## 4.2. Stability Analysis on the DIDEs (4) with TC on x2

In accordance with the previous section, the effect of the decision adjustment of player II (safety inspection groups inside coal mine enterprises) on CMSM is studied in this section. When there is a threshold constraint (TC) onx2, we get the following DIDEs:(5)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,x2>x2∗,Δx2t=α,α≥0Δx3t=β,β≥0x2≤x2∗.Here, we also consider just one case, i.e.,x2 individually imposed a pulse (α), to analyze the evolutionary stability of the dynamic system (4). Some of the simulation results are shown in Figure 7.Figure 7
TSGs ofxi with TC on x2 when x2 has a pulse β=0.2. x1,x2,x3|t=0=P00.5,0.3,0.2. (a) x2∗=0.2, (b) x2∗=0.3, (c) x2∗=0.4, and (d) x2∗=0.5.
(a)(b)(c)(d)As can be seen from Figure7, the dynamical system (5) also has a satisfactory evolutionary stability when x2∗ is large enough (not less than 0.3). Therefore, it is also concluded that the decision adjustment of safety inspection groups can promote the safety management level of coal mining.


## 4.3. Stability Analysis on the Dynamical System with More TCs

We have studied the evolutionary stability of the dynamical system with TC imposed on the variablesx1,x2, and x3, respectively, in Section 3 and Sections 4.1 and 4.2. As a neutralization of the previous cases, we consider the following DIDEs, where there are more TCs in CMSM:(6)x˙1t=x1t−τ11−x1t−τ1x2t−τ1+x3t−τ3−x2t−τ2x3t−τ3f1−c1,x˙2t=x2t−τ21−x2t−τ21−x1t−τ1f1+x3t−τ3f2−c2,x˙3t=x3t−τ31−x3t−τ31−x1t−τ11−x2t−τ2f1+f2−c3,xi>xi∗,i=1,2,3,Δx2t=α,α>0Δx3t=0x1≤x1∗,x3>x3∗,Δx2t=α,α>0Δx3t=β,β>0x1≤x1∗,x3≤x3∗.To reflect the impact of TCs, we set pulse periodT=3 in the model simulations of the DIDES (6), and thus, the pulse check frequency is higher than before. By contrasting Figures 6–8, the results of the simulation show that the system stability in this case is no better than the previous two cases. The volatility of x1 may last forever (see Figure 8).Figure 8
TSGs of the DIDEs (5) with more TCs α=0.3,β=0.2. Pulse period T=3, x1,x2,x3|t=0=P00.5,0.3,0.2. (a) x1∗=0.6,x2∗=0.2,x3∗=0.1, (b) x1∗=0.8,x2∗=0.2,x3∗=0.1, (c) x1∗=0.9,x2∗=0.3,x3∗=0.1, and (d) x1∗=0.9,x2∗=0.2,x3∗=0.3.
(a)(b)(c)(d)But on the other hand, these results can also be acceptable because of the relatively large values ofx1, especially when the level of safety management technology is relatively low, which is still common in some remote areas or small- and medium-sized coal mining enterprises in China [37]. The satisfactory evolutionary stability is difficult to achieve in the practice of CMSM. As shown in Figures 8(c) and 8(d), once the pulse heuristic conditions (the TCs in the DIDEs (6)) are set properly, the dynamic system can make moderate fluctuations over a reasonable range, where the troughs of x1 are almost not less than the set value of x1∗. Meanwhile, all participants in CMSM can restrict and promote each other and actively engage in safety management. That is exactly in line with the aspirations of the Chinese government: to deal with safety management, all participators should take part in the whole process and manage it together.In this section, stability analyses are made on the dynamical system with different TCs. Table1 shows the simulation results briefly. When the TCs are imposed on x2 or x3 separately, x1 is evolutionary stable but one of the two supervisors performs inactively. In the last case, when more TCs are combined in the dynamical system, it is gratifying that x1 is progressively stable and both the other two have active performances. It is a good revelation for us. We should appropriately adjust the TCs and the thresholds of them, so that the management effect reaches the actual demand of CMSM in practice.Table 1
Evolutionary stabilities of the dynamical system with different TCs.
TCsStability ofx1Performance ofx2Performance ofx3Picturesx3≤x3∗=0.3Evolutionary stableInactiveActiveFigures6(b)–6(d)x2≤x2∗=0.3Evolutionary stableActiveInactiveFigures7(b)–7(d)More TCs combinedProgressively stableActiveActiveFigures8(c) and 8(d)


## 5. Conclusions

In this paper, a series of evolutionary dynamical system models governed by the DIDEs are established to study the behaviors communications in CMSM, and stability analysis on the dynamical system with various TCs and pulse responses is carried out. Satisfactory evolutionary stability can be achieved when TCs and pulses are imposed onx2 or x3, which is very valuable in the practice of CMSM. The TCs and pulses in the models can accurately correspond to the actual inspection and supervision measures in safety management. Therefore, the research conclusions can provide a theoretical basis for the scientific analysis methods for the spread of unsafe behaviors and the formulation of control strategies of the supervisory measures in safety management.Based on the main results of this article, we attempt to put forward the following suggestions. To improve the management efficiency of coal mine safety, it is not enough to impose various safety regulations on coal miners separately. All the efforts of government departments, safety management departments inside enterprises, and other parties should concentrate to work together. When formulating supervision measures, it is necessary to make targeted arrangements for the safety supervision and inspection commissions of government departments and enterprise’s safety inspection groups, so as to highlight the leading role and supervision effectiveness of these safety authorities. Besides, reasonable settings of pulses and TCs can effectively improve the stability of the evolution model and enhance the application value of this research, which require continuous testing in practice of coal mine.The method of DIDEs is used to carry out the evolutionary stability analysis of behavior propagations in CMSM. It is also the research significance and main innovation of this article, yet there are still some details here for further research and analysis. In future, we can further analyze and study the time unit of the variables of the DIDEs in this study, as well as the appropriate values of pulses, the pulse periods, pulse heuristic conditions (TCs), and impulse response function,. The model research must be closely linked to the practices of specific coal mine enterprises so as to improve the accuracy of the conclusions. In addition, the settings of different model parameters may also cause the game dynamic system to enter a chaotic state or have a bifurcation [38–40]. This may cause difficulties in safety management, but it may also be an opportunity to improve the level of safety management.



---

*Source: 6210465-2020-05-28.xml*