Abstract
With the rapid development of urbanization and modernization, people’s requirements for electric power are also consequently increasing. It is thus essential to ensure the sustainable and stable supply of power and maintain the safety performance of the power grid. Constructing a modern safety management system (SMS) can effectively improve power enterprises’ ability to respond to emergencies. To effectively evaluate the ability of electric power enterprises to resist risk impact, and the recovery and optimization ability after impact, this paper makes a systematic investigation of the SMS of power enterprises. It evaluates the safety resilience capabilities of the SMS from four basic elements of stability, redundancy, efficiency, and adaptability. A multilevel evaluation system of electric power enterprises’ safety resilience capacity consists of four first-level and 25 second-level indices. Then, an analytic hierarchy process (AHP) and matter-element extension (MEE) joint model—the AHP-MEE evaluation model of the resilience capacity of power enterprises—is proposed. A case study of electric power enterprise safety resilience evaluation is conducted using the proposed resilience evaluation index system and assessment model. The results showed that the provided model agrees with the actual situation and verified the index system and model for electric power enterprise safety reliability and applicability. Improved enterprise safety resilience recommendations are proposed accordingly based on the evaluation results. This work proposes a robust safety resilience capacity evaluation method in the face of upcoming new risks, thus building a good working foundation to construct a new power system with new energy sources.
1. Introduction
The safety of power infrastructure is a concern in many countries. In recent decades, there have been many severe electrical safety incidents, such as the great blackouts of the United States and Canada, the blackout in London in 2003, the blackout in Western Europe in 2006, the blackout in China in 2008, and the blackout in Venezuela in 2019. These accidents always result in significant economic losses and even casualties [1]. Afterward, countries worldwide realized the importance of the power grid’s resilience and other vital infrastructure to withstand major safety risks and gradually shifted their focus from simple disaster prevention to enhancing the power grids’ ability to recover from disasters and accidents. Concurrently, along with China’s primary energy strategy of “carbon peak” and “carbon neutrality,” the development momentum of large-scale grid connection of new energy and clean energy and massive access to new energy facilities has strengthened [2]. The form of the power grid has become increasingly complex, and the pressure on risk prevention and control is increasing. Higher standard requirements are proposed to improve the power grid’s ability to deal with sudden disasters, reduce the risk and loss of large-scale power outage accidents, and guarantee the power grid’s reliability to the maximum extent [3–5].
A safety management system (SMS) can be described as a planned, documented safety program that incorporates basic management concepts and activating elements into a well-organized safety system [6], which becomes a standard tool in the safety management of large electric companies in China. Many studies have been performed from different perspectives on the design, construction, and modeling of the SMSs [5, 6]. Collier [4] proposed a ten-step plan to create a modern, intelligent, and agile power grid system based on the practical reasons for constantly changing electricity demand and frequent failures of power lines and equipment, which provides a new idea and method for power enterprises to implement intelligent management. Lu et al. [5] identified the basic concepts of power grid operation safety management, described the relationship between different standardized terms through the Venn diagram method, and finally proposed a set of standardized concept systems of power grid operation safety risk management.
For electric power enterprises, the ability to resist risk impact and recover and optimize after impact is an integral part of the safety management level of the power system. Xie et al. [7] simulated the grid bearing capacity under extreme weather, analyzed grid vulnerability and risk tolerance under different environments, and identified the nodes and lines that are highly vulnerable to extreme weather, providing theoretical support for preventing and reducing the impact of extreme weather on the grid. Bian et al. [8] dynamically analyzed the evolution process from the cause to the result of power grid trip accidents and quantitatively calculated the probability of power grid trip accidents caused by natural disasters, providing technical support for the prevention, treatment, and recovery of power grid trip accidents. Although many efforts have been made in constructing SMS, the previous studies primarily investigated accident response and theoretical surveys; few addressed emergency prevention and the recovery capacity from a disaster, which is unfavorable to the safety management of the power grid. It is vital to evaluate the ability to “cope with disasters” of power grids affected by emergencies and to study how to improve the ability of the power grid to manage emergencies and ultimately achieve the goal of maintaining power grid performance reducing accident losses. It is critical to study and construct an emergency safety management system for power grids with resilience.
In this paper, we reconstruct the power grid emergency state and SMS by introducing the concept of “resilience.” Based on the resilience theory, an optimized AHP-MEE safety resilience capability assessment model is provided to evaluate the electric power enterprise’s ability to deal with emergencies. Engaged in the former matter and afterward, the three stages of response blind spot existing in the safety management system, a corresponding power grid system construction is proposed. This process provides the theoretical basis for constructing an emergency safety management system and has fundamental practical significance for developing a modern safety management system.
2. Basics of the Resilience Theory
The resilience theory has become an essential theoretical foundation in some frontier disciplines and is expected to be integrated into more disciplines [9–12]. A power system’s recovery ability, disaster resistance ability, and postdisaster response ability after suffering various types of accidents and disasters embody its macro-resilience ability. Resilience, as an inherent capability of the system, emphasizes a system’s ability to absorb, adapt, recover, and improve external disturbances and has attributes such as situational awareness, resistance, robustness, responding capacity, recovery, timeliness, and learning ability [13–16].
2.1. Concept of System Safety Resilience
For a system, safety defense and postrecovery capabilities are essential factors that affect system resilience. Therefore, safety resilience is a valuable research subject at the intersection of resilience science and safety science. Bergström et al. [17] discussed the rationality and necessity of applying resilience theory in safety science. They believed that system safety resilience is a dynamic system attribute closely related to continuous adjustment and adaptability, developed based on vulnerability, safety capacity, and safety redundancy. System safety resilience refers to the ability of the system to absorb and adapt to external disturbances to maintain its stable safety state. It also reflects the system’s resilience, learning, and adaptability, as shown in Figure 1 [18]. System safety resilience is an essential attribute of system safety and is a critical approach to comprehensively understanding the safety attributes of complex systems [19].

Generally, system safety resilience elements can be divided into the following parts, which jointly determine the safety resilience of the system [20].(I)Stability is the ability of the system to withstand risk shocks and disturbances before failure(II)Redundancy is the ability of the system to replace parts or components affected by risk and interference(III)Efficiency is the speed at which the system recovers from the damaged state to the initial state(IV)Adaptability is the ability of the system to identify safety problems and determine priorities, and the ability to mobilize emergency and disaster relief resources can be further summarized as the availability of human, material, and financial resources in the face of an accident.
The interaction patterns between these four resilience elements and the existing organizational member system, material technology system, management system, and environmental system in the SMS of the enterprise are shown in Figure 2.

2.2. Framework of System Safety Resilience
Based on this analysis, indicators can be defined by the dynamic response curves of the system before, during, and after the event to build the basic framework of enterprise system safety resilience, as shown in Figure 3 [21].

According to Figure 3, when 0 < t < t1, the system is in the original safety state. When the disaster-bearing capacity of the system cannot absorb risks, t1 is subjected to risk shock and disturbance [21].(1)When the risk disturbance and destruction exceed the maximum total risk (safety capacity) borne by the system to maintain safe operation, accidents occur, the system is damaged, and the system enters the damage stage (t1 < t < t3). At t3, recovery measures are performed, and the system enters the recovery stage (t3 < t < t4). In the recovery phase, the system may or may not recover to presafety state point B and can also improve its safety resilience and recover to a better safety state point A.(2)The disaster-bearing capacity of the system cannot absorb risk, resulting in damage and loss of the system, and the system cannot recover to the initial safety state but a new state point C. The system still cannot withstand future risk shocks.(3)The disaster-bearing capacity of the system is sufficient to absorb the risk impact so that no damage or loss occurs. At this time, the system can maintain its presafety state at B and improve its resilience function at A through learning the ability to cope with risk shocks to manage more intense risk shocks in the future.
The quantification of system safety resilience is conducive to applying safety resilience theory. When an accident occurs, the safety state of the system will be broken, and the recovery time and efficiency of the system safety state will depend on the safety resilience capacity of the system.
3. Modern SMS for Electric Power Enterprises
With the diversified development of safety theories, accident causation theory, accident causation linkage theory, and safety resilience theory have been applied to the construction of enterprise safety management systems. We now consider the Zhejiang Electric Power Corporation in China as an example. Through continuous exploration and practice, the enterprise gradually formed a diversified and highly flexible modern SMS called the “1438” management system, where “1” represents the building mainline of the diversified and integrated high-flexibility power grid; “4” represents the four guidelines of safety policy, safety objective, safety culture, and safety theory; “3” represents the construction of safety organization system, safety work system, safety rules, and regulations system of three systems; and “8” refers to the improvement of organizational management, team quality, and other eight capabilities. These numbers are the premise of building a “three” system composition.
As shown in Figure 4, the premise of the construction and implementation of an SMS is to establish a safety policy, safety objectives, and safety culture based on safety theory. Establishing the safety organization system, safety work system, and rules and regulations clarifies and improves the responsibility, work content, and management basis of enterprise safety management [22, 23]. The details are as follows.(1)The safety organization system is supported by the safety responsibility system led by the safety production committee by creating a safety assurance system, supervision system, and assurance system to fulfill the safety subject responsibility, supervision and coordination responsibility, and support responsibility, respectively.(2)The safe working system aims to achieve continuous improvement and improvement by creating the working thread of preprevention ⟶ in-process control ⟶ postprocess treatment, supplemented by external cooperation.(3)The safety rules and regulations system is essential for constructing a safety organization system and work system. The timeliness and accuracy of this system are guaranteed by constantly supplementing the on-site regulations and timely following up on the latest laws and regulations, rules and regulations, and technical standards.

In addition, to ensure the safety organization system and a safety work system, rules and regulation system is correctly constructed and implemented. As shown in Figure 5, constructing a modern SMS of electric power enterprises improves the eight capabilities of organization management, team quality, source prevention, hidden trouble rectification, risk control, emergency response, scientific and technological innovation, and safety investment. To ensure the rationality of capacity building, the safety performance index, safety management index, rules and regulations implementation index, and scientific and technological innovation index are set in each capacity-building improvement process [24].

Although the existing SMS of electric power enterprises is performed based on the three aspects of the organization, system, and safety work, four representative indices, eight capabilities, and 15 safety evaluation indices are proposed to construct the system. However, the system’s resilience level cannot be directly determined by evaluating the above indicators in response to emergencies. Therefore, from the perspective of “safety resilience,” the safety evaluation index is organically combined with four resilience factors: stability, redundancy, efficiency, and adaptability. Additionally, the safety resilience of enterprises is quantified by safety evaluation [25].
4. Construction of the Resilience Evaluation Model
A complex system engineering is used to construct electric power enterprises’ resilience capacity assessment model. Based on the principles of science, objectivity, and systems, and combining the characteristics of electric power enterprises, this paper integrates the modern safety management system and system safety resilience theory. To ensure the rationality of the model, this paper reconstructs the existing evaluation index system from the perspective of safety resilience. It uses the analytic hierarchy process (AHP) and matter-element extension model to divide and quantify the indices, starting from the four elements of stability, redundancy, efficiency, and adaptability. An enterprise safety resilience assessment model composed of several primary and secondary indicators is constructed. The evaluation process is shown in Figure 6 [26].

4.1. Determination of Safety Resilience Evaluation Index
To effectively assess the safety resilience level of an enterprise, this paper refines the safety resilience elements of each index based on the 15 safety evaluation indicators of the existing safety management system. Four safety resilience elements are used to assess enterprise safety resilience. The 25 refined evaluation indicators are used as the second-level indicators for the first-level indicators. The corresponding relationship between the first-level indicators and the second-level indicators is established from safety resilience. Table 1 shows the refined resilience evaluation indicators [27].
The resilience evaluation indicators are considered secondary indicators and divided according to stability, redundancy, efficiency, and adaptability. The results of the classification are shown in Figure 7.

4.2. Determination of Indicator Weight
The basic principle of using AHP to divide indicators in this paper is now described. Experts construct a 4 × 4 judgment matrix for each indicator’s resilience elements, scale them, and then use the root-finding method to find its eigenvectors and conduct consistency tests. This method determines whether the experts are inconsistent when judging the importance of the index [28]. For an expert who passes the consistency test, the eigenvector of the judgment matrix constructed by the expert is the effective weight vector assigned by the expert to each resilience element in a specific index. The effective weights assigned by all experts who passed the consistency test are averaged according to the indicators. Then, the safety resilience weights determined by the expert group can be obtained. The resilience element with the most considerable weight is determined as the primary relevant element of the indicator [29].
The analytic hierarchy process is divided into the following steps to determine the weight.(1)Clarify the problem and establish a hierarchical structure(2)Construct a pairwise comparison judgment matrix (see Table 2). where is the expert’s comparison result of index factors and .(3)The weight vector is calculated from the judgment matrix, and the consistency test is performed(1)Integrate the judgment matrix by rows and then find the roots to obtain the vector (2)Normalize the vector to obtain the eigenvector of the judgment matrix (3)Calculate the maximum eigenvalue of the judgment matrix where is the -th element of the vector .(4)Calculate the consistency index of the judgment matrix (5)Calculate the random consistency ratio CR and judge the consistency of matrix ; the value of is shown in Table 3. If the consistency index and the average random index ratio are less than 0.1, the consistency test is satisfied.(4)Determine the weight Anyone who passes the consistency test is an expert who effectively weights the evaluation index. The eigenvector W = [w1, w2, …, wm] of the constructed judgment matrix M is the effective weight vector assigned to the evaluation index. Those who fail the consistency test are experts who provide the evaluation index with an invalid weight. According to the evaluation index, the arithmetic mean of the effective weight assigned by each effective weighting expert in the expert group to the corresponding index is calculated, which is the weight assigned by the expert group to each evaluation index.
4.3. Multilevel Extension Evaluation of Enterprise Safety Resilience
4.3.1. Determination of Classical Domain, Node Domain, and Matter to be Evaluated
A resilience evaluation index system composed of various indicators “ci” was developed, and U = {u1, u2, u3, u4} was used to represent the enterprise safety resilience level domain, as shown in Table 4 [30, 31].
(1) Determine the Classical Domain. We let be the classical domain of the matter element , as follows:where is the safety level in the enterprise safety resilience level domain ; is the value range of the evaluation index set C regarding the safety resilience level ; and are indicators; and ci is the lower limit and upper limit of the safety resilience level, respectively.
(2) Determine the Section Domain. We let be the node domain of the matter element , as follows:where is the value range of the evaluation index set concerning the resilience level domain , and and represent the lower limit and upper limit of the index , respectively, at all resilience levels.
(3) Determine the Evaluation Matter Element. is the matter element to be evaluated for matter element , as follows:where is the -th first-level indicator to be evaluated for the level of corporate resilience; is the set of second-level evaluation indicators for , where is the -th second-level indicator of the ci evaluation index; is the evaluation and assignment value of cij by the expert, and is the set of evaluation and assignment values of each secondary evaluation index of by the expert.
4.3.2. Evaluation of Secondary Index Safety Resilience
(1)Calculate the relevance of secondary evaluation indicators. The degree of relevance represents the closeness of the object element to be evaluated to the classical domain. The concept of classic mathematics mid-distance is introduced to establish the correlation function of the second-level evaluation index cij of the company’s safety resilience evaluation with the company’s safety resilience level . Then, the first level can be calculated. The degree of relevance of the jth secondary index of the index ci for the safety resilience level is as follows: where and are the distances from to and , respectively, as follows:(2)Determine the safety resilience level of the secondary evaluation index. Among the correlation values of a specific secondary evaluation index concerning each level, the safety resilience level can be determined according to the principle that the level corresponding to the maximum value is its safety resilience level. Thus, when , the safety resilience level of the secondary evaluation index is level .
4.3.3. Evaluation of First-Level Index Safety Resilience
(1)Calculate the relevance of the first-level evaluation index. The correlation matrix K(ci) of the first-level evaluation index of enterprise resilience to each safety resilience level is related to the weight vector of each second-level evaluation index and the correlation matrix of the second-level evaluation index to each safety resilience level. We then multiply the values as follows:(2)The method for determining the safety resilience level of the first-level evaluation index is the same as above.
4.3.4. Evaluation of the Company’s Comprehensive Safety Resilience
(1)Calculate the comprehensive correlation degree of enterprise safety resilience. The comprehensive correlation degree of enterprise safety resilience to each safety resilience level is multiplied by the weight vector of each first-level index and the correlation matrix of each first-level index to each safety resilience level is determined as follows:(2)The method for determining enterprise safety resilience is the same as above. The various evaluation indicators are divided and quantified through the above methods, and the safety resilience assessment of the modern safety management system of the enterprise is performed. The weak links of the safety resilience of the enterprise are then discovered so that the modern safety management system of the enterprise can be improved and supplemented.
5. Case Study of the Resilience Assessment Model
5.1. Basic Profile of Case Company
In this part, a branch company of State Grid Zhejiang Electric Power Corporation Ltd. (referred to as company N) was selected as a study case. It was established in 2008 and is one of the 11 municipal power supply companies in Zhejiang Province, China. It is primarily responsible for local power production, deployment, plant construction, and other services. In the past years, company N paid much more attention to safety management, continuously strengthening the risk management and control capabilities of operations, equipment, and personnel on the power generation side, power supply side, and user side. Also, some suitable measures have been taken in the past two years to promote the modernization of safety management systems and governance capabilities. In terms of safety information and digital construction, relying on the primary safety domain of the local company’s business, company N actively promotes the research and development of digital safety management and control technology, the construction of emergency command systems, and the promotion of safety tools and equipment, and promotes the construction of a safety culture based on the company’s characteristics. In addition, company N also established a research and development base for new energy sources, such as energy storage and hydrogen production technology, wind power, and solar energy. As a result of this, we made an in-depth questionnaire survey to decide the weight values and indicator domain values, with a wide range of the questionnaire personnel, such as technicians, engineers, managers, and administrative staff from company N. A total of 26 questionnaires were obtained, and 11 valid questionnaires were selected based on the objectivity of the answer sheets (abnormal results were considered invalid) and the time required to fill out the answer sheets (<2 min was regarded as invalid). The obtained quantitative value of related indicators and the basic overview of the company’s resilience is shown in Table 5 [32].
5.2. Resilience Assessment
5.2.1. Classic Domain, Knot Domain, and Determination of Matter to be Evaluated
(1) Determine the Classic Domain Rg. Based on the literature, we summarized the research content related to resilience done by the predecessors, combined the opinions of experts in the safety field, divided the enterprise’s resilience capability level boundary, and obtained the enterprise safety resilience level domain . Thus, the classical domain of each level of power, the enterprise is determined by formula (7) as follows [27]:
(2) Determine Section Domain. . The section domain is determined by formula (8) as follows:
(3) Determine the Object to be Evaluated. . This article takes the first-level index stability as an example, and the matter element to be evaluated is determined by formula (9) as follows:
5.2.2. Evaluation Index Weight Calculation
(1) First-Level Indicator Weight Calculation. According to the company’s actual situation, the ten invited experts independently compared the first-level indicators and constructed their judgment matrix to obtain the weight of the first-level indicators in the safety resilience level of the company. Considering expert (1) as an example, the method proceeds as follows.(1)Determine the target of the analysis as the “overall resilience level of the enterprise” and compare the four elements of safety resilience of the enterprise(2)Construct a judgment matrix M and determine its characteristic vector①Construct the judgment matrix M of the four resilience elements in this index, and the method is shown in formula (1). The judgment matrix M constructed by the expert is②Integrate the judgment matrix by rows and then find the roots. The method is shown in formula (2), and the matrix is described by③Normalize the matrix . The method is shown in formula (3), and the eigenvector of the judgment matrix is described by(3)The random consistency ratio of the matrix is calculated, and the consistency test is performed①Calculate the maximum eigenvalue of . The method is shown in formula (4), and the result is as follows:②The method for calculating the consistency index of the judgment matrix is shown in formula (5), and the results are as follows:③The method for calculating the random consistency ratio is shown in formula (6), and the results are as follows: From Table 1, when , is 0.89, then Because , the judgment matrix constructed by the experts passes the consistency test.(4)Determine the weight of each resilience element in the index. Because the judgment matrix constructed by the experts passes the consistency test, the weighting resilience element is effective. The eigenvector of the judgment matrix constructed is the weight assigned to its resilience element. Thus, , , , .(1)Second-level indicator weight calculation: Similarly, the eigenvectors of the judgment matrix constructed by the resilience factors in the resilience level of other experts in the expert group can be calculated. Eight experts passed the consistency test, and the weights assigned to them were calculated and averaged. The results are shown in Table 6.
The relationship between the level of safety resilience of the company and the first-level indicators is shown in Figure 8.

Similarly, each of the ten invited experts constructs a judgment matrix independently for each of the second-level evaluation indicators of the company’s safety, resilience, stability, redundancy, efficiency, and adaptability and can be calculated. The weights of the secondary indicators under each primary indicator are shown in Table 7.
5.2.3. Multilevel Extension Evaluation of Enterprise Safety Resilience
(1) Safety Resilience Evaluation of the Secondary Index. Considering the calculation method of the correlation degree of the second-level index with the second-level evaluation index to the first-level evaluation index and the calculation of the correlation degree of enterprise safety resilience level 1 as an example, we obtain the following based on formulas (10)–(12).
By analogy, the correlation degree of concerning other resilience levels can be calculated as
Similarly, the correlation value of other secondary indicators for each safety resilience level can be calculated. Then, the formula determines the safety resilience level of each secondary evaluation index, as shown in Table 5.
(2) First-Level Index Safety Resilience Evaluation. Considering the first-level index as an example to introduce the first-level index relevance calculation method from formula (13), the correlation value of each resilience grade is determined as
Similarly, the correlation values of other first-level indicators for each safety resilience level can be calculated separately.
According to the formula , the safety resilience level of each first-level evaluation index is determined, as shown in Table 5.
(3) Evaluation of Comprehensive Safety Resilience of Enterprises. According to formula (14), the total correlation value of the safety resilience level of the enterprise for each safety resilience level can be calculated as
According to , the overall safety resilience level of the enterprise can be determined to be level 2 (i.e., “strong safety resilience”).
5.3. Evaluation Results and Analysis
The power company’s safety resilience grading evaluation results are shown in Table 5. After evaluation, the power company’s comprehensive safety resilience rating is level 2 (i.e., “stronger safety resilience”). The overall resilience level is relatively high, but there is still room for improvement. Specifically, from the first-level index safety resilience evaluation results, except for the first-level index stability, the other three first-level evaluation indicators are at the second-level safety resilience level (“more robust safety resilience”), which shows that each first-level evaluation indicator has good defense and resilience to risks. However, some aspects still must be improved. From the evaluation results of the second-level index safety resilience, in addition to the nine second-level indicators such as organizational completeness and responsibility matching, the safety resilience level is level 1 (strong safety resilience), and the safety resilience levels of the remaining second-level indicators all are level 2 (more critical safety resilience). According to the evaluation results and the above analysis results, it is recommended that the electric power company comprehensively implements safety risk rectification while focusing on the following rectification work to ensure the production safety of the electric power company and to prevent gradual failures and prevent problems before they occur.
In terms of stability, it is necessary to promote the implementation of the safety production responsibility system through a detailed inventory, digital empowerment, and personnel training, strengthen the supervision of the implementation of the “two measures,” check the completion of quarterly, find problems and promptly report them, strengthen the process control and assessment of the outsourcing team, establish the credit evaluation system of the outsourcing team, improve the independent management ability of the team, and promote and popularize the enterprise “safety QR code.”
In terms of redundancy, it is necessary to improve customer-side safety management, communicate with customers quickly about existing problems, propose specific rectification measures, promote intelligent process control, accurately control the entire safe production process, and reduce communication costs.
In terms of efficiency, it is necessary to speed up the construction of emergency command centers and enhance the emergency response capabilities of enterprises to better respond to emergencies; timely feedback on the rectification of safety hazards should also be made to prevent production accidents.
In terms of adaptability, it is necessary to educate and train employees, focus on and improve the qualification rate of education and training, and improve the overall safety level of the enterprise; regularly performing training and inspecting emergency personnel and equipment can ensure the safety of emergency personnel quality and the availability of emergency equipment.
6. Conclusion
This paper investigates and evaluates the electric power enterprises’ resilience capability. The primary conclusions of this study are as follows.(1)Based on resilience theory, the evaluation indicators of electric power enterprises are reclassified on the basis of the “1438” SMS, and the resilience capacity evaluation index system of power grid enterprises is proposed, with four first-level indices, stability index, redundancy index, efficiency index, and adaptability index and 25 second-level indices.(2)An optimized AHP-MEE joint model is proposed; it provides a method to calculate the weights of evaluation indicators, precise analysis of the scope of each indicator, and exact calculation method for system safety resilience.(3)A case study is conducted to verify the reliability and applicability of the provided safety resilience capacity evaluation index system and the AHP-MEE model. It is determined as a level 2 (stronger safety resilience) for the selected sample, consistent with production reality. The evaluation results also provide targeted rectification and reform measures to strengthen safe production work.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (grant numbers 51976205, 51904282, and 51804287).