Building fires are characterized by high uncertainty, so their fire risk assessment is a very challenging task. Many indexes and parameters related to building fires are ambiguous and uncertain; as a result, a flexible and robust method is needed to process quantitative or qualitative data and update existing information when new data are available. This paper presents a novel model to deal with the uncertainty of the residential building fire risk and systematically optimize its performance effectiveness. The model includes fuzzy theory, evidence reasoning theory, and expected utility methods. Fuzzy analysis hierarchy process is applied to analyze the residential building fire risk index system and determine the weights of the risk indexes, while the evidence reasoning operator is used to synthesize them. Three buildings were selected as a case study to illustrate the proposed fire risk model. The results show that the fire risk level of three buildings corresponds to “moderate” or below which is consistent with the previous study. These results also truly reflect the actual situation of fire safety in these residential buildings. The application of this model provides a powerful mathematical framework for cooperative modeling of the fire risk assessment system and allows data to be analyzed step by step in a systematic manner. It is expected that the proposed model could provide managers and researchers with flexible and transparent tools to effectively reduce the fire risk in the system.

1. Introduction

With the acceleration of industrialization, urbanization, and marketization in China, building construction industry has developed rapidly. Particularly, the structure and function of buildings are becoming more complex, and various new technologies and techniques are emerging constantly, which have led to the increasingly severe situation of building fires. According to the statistics provided by the Ministry of Public Security in 2013, a total of 388,821 fires were recorded in China, in which 52% (202,299) of fires occurred in buildings, resulting in 3410 civilian deaths or injuries and 3760 million Chinese yuan (CNY) direct property losses. Nowadays, building fire is considered to be an enormous threat to people’s life and production in China, and a growing concern is how to take appropriate measures to reduce the fire risk, minimize the damage and loss caused by fire in buildings, and guarantee building fire safety. Therefore, it is urgent to establish a suitable fire risk assessment model, and it provides information through quantitative or qualitative analysis results to make decisions on whether to take steps to reduce the risk [1, 2].

There are mainly four conventional types of fire risk analysis methods: checklist, description, index, and probability method [3]. However, most of these approaches have prescriptive drawbacks which make them difficult to quantitative fire risk analysis due to the inability to deal with the uncertainties associated with the fire risk factors of the system. With the improvement of performance-based fire-protection design, some fire risk analysis models and corresponding software have emerged, such as FiRECAM™ (Fire Risk Evaluation and Cost Assessment Model) [4, 5], FIERAsystem (Fire Evaluation and Risk Assessment system) [6], CESARE-RISK (Centre for Environment Safety and Risk Engineering, RISK) [7, 8], and Crisp II (Computation of Risk Indices by Simulation Procedures) [9]. However, these models should depend on some strict constraints, such as a large number of input data, specific fire scenarios, and the large amount of calculations. Consequently, researchers are concentrating on developing flexible fire risk analysis tools based on systematic safety theory. For example, Ibrahim et al. presented a fire risk method based on the analytical hierarchy method (AHP) for heritage buildings [10]. Lo developed a fire risk ranking system for existing buildings using the fuzzy set approach [11]. Liu et al. built a fire risk analysis system for commercial buildings by using the structure entropy weight method [12]. Xin and Huang proposed scenario cluster methods in the process of the fire risk analysis model for residential buildings [2]. Briefly speaking, these methods reveal two main challenges in an uncertain environment associated with the fire risk factors of the system. The first challenge faced by these methods is the lack of the ability to process a variety of data suitable for fire risk reasoning mechanisms, and the second is the lack of the ability to analyze the interdependence of risk factors. In this paper, a fire risk analysis model integrated fuzzy theory, and evidential reasoning (ER) theory is presented for residential buildings. Compared with the traditional fuzzy reasoning approach, ER has the advantage of avoiding losing useful information; therefore, it can be applied to model complex systems. The framework of this model is organised as follows. Section 2 illustrates the methodology of the research. Section 3 presents a case study to verify the feasibility of the methodology. Sections 4 and 5 discuss the empirical results and conclusions.

2. Methodology

Quantitative risk assessment (QRA) techniques are usually used for assessing uncertainties in building fires. However, due to the lack of fire accident statistics, an effective solution is to integrate expert judgments into the QRA process. QRA consists of four main procedures: hazard identification, occurrence probability calculation, consequence severity assessment, and risk quantification [13, 14]. In order to process the complex system structure and promote a flexible implementation method, different decision-making techniques can be used, such as fuzzy analytic hierarchy process, fuzzy set theory, and evidence reasoning method. Due to the fact that fuzzy logic could provide a flexibility way to represent the vague information resulting from the lack of data or knowledge. Therefore, the fuzzy set theory has a wide application in different fields such as reliability engineering, system safety, and risk assessment [15].

The proposed framework, shown in Figure 1, allows step-by-step analysis of the utility tunnel fire risk in a transparent way, as described as follows:(1)Identifying fire risk factors and establishing the hierarchical structure of the index system(2)Using fuzzy analytic hierarchy process (FAHP) to calculate the weights of indexes(3)Applying the belief degree structure based on the fuzzy set theory to measure the fire risk(4)Aggregating the result of the fire risk using the evidence reasoning (ER) algorithm(5)Using the expected utility method to obtain a clear result of the fire risk(6)Sensitivity analysis

2.1. Identifying Fire Risk Factors and Establishing the Hierarchical Structure of the Index System

In order to make better decisions on fire control protection and emergency evacuation measures, a structured and systematic approach is needed. It is better to describe the fire risk problem in a hierarchical structure so that decision makers could have a thorough understanding of the system, especially when it is a complex system with multilevel structural indexes.

According to NFPA550 Guidelines, to achieve fire safety, reducing the fire risk mainly starts from two aspects: one is to prevent the occurrence of fire, and the other is to control the impact of fire [16]. In this paper, fire risk factors of these two aspects are, respectively, defined as disaster-causing factors and loss-controlling factors. Disaster-causing factors may cause the fire risk to be transformed into disaster before fire occurs, while loss control factors signified various fire protection and management measures to control the development process of fire and mainly involved four aspects: passive measures, active measures, fire management, and fire brigade fighting.

Based on the characteristics of residential building fire and the literature review [6, 1720], the factors influencing the risk of building fires are analyzed from the two aspects of disaster-causing factors and loss-controlling factors. A general hierarchical structure (presented in Figure 2) is finally established after theoretical preparation, the initial construction of the index system, the optimization of the index system, and the determination of the index system.

2.2. Fuzzy Analytic Hierarchy Process (FAHP)

The traditional analytic hierarchy process (AHP) constructs the judgment matrix by comparing the two factors with the 1–9 scale method. However, due to the subjectivity of human judgment, different people will get different conclusions. Using the triangular fuzzy number to scale the two-pair comparison of factors can consider the uncertainty of experts in analysis and judgment, which give the range of expert’s judgment in the form of intervals to reduce subjectivity. In 1996, Chang [21] applied triangular fuzzy numbers to construct judgment matrices and combined with the extent analysis method to calculate the weights of each index in the hierarchical structure. Finally, the traditional AHP is transformed into the FAHP in the fuzzy environment, which can provide more practical results [22].

2.2.1. Triangular Fuzzy Number

Suppose the triangular fuzzy number is , and its membership function is equal to

Herein, , and represent the lower and upper boundary value of triangular fuzzy number , respectively, and represents the median value of triangular fuzzy number . Generally, triangular fuzzy number can be abbreviated as . Let and be triangular fuzzy numbers; then, the possibility degree of is defined as follows:

2.2.2. Fuzzy Synthetic Extent

Consider as a set of analytic objects and as a target set; we can get the extent value of the -th object satisfying the -th goal, in which the sign is . Then, the value of synthetic extent of the -th object is defined as [21, 23]

2.2.3. The Procedure of the FAHP

In the evaluation of the fire risk, the determination of the weight of each fire risk factor is particularly important. The weight represents the relative importance of each factor in the overall evaluation. Only when the weight of each factor is obtained, the fire risk assessment can be carried out. The steps of determining the weight by the FAHP method are as follows:(1)According to the objective of fire risk assessment, the hierarchical system structure is established, which is composed of fire risk factors.(2)The judgment matrix is constructed by triangular fuzzy numbers (according to Table 1) through a pairwise comparison of the index system by experts [24, 25].(3)According to equation (3), the value Si of synthetic extent of each factor is obtained.(4)The possibility degree is calculated such that factor is more important than others:

Then, the weight vector is obtained:

Finally, the normalized weight vector is obtained.

2.3. Application of the Belief Structure for the Fire Risk Calculation

After identifying fire risk factors and establishing the hierarchical structure of the index system, another important task of risk management is to assess the risk, which is an effective way to prevent or reduce the effect of the fire [2]. In this paper, the fire risk of residential buildings is defined as the result of comprehensive measurement associated to the occurrence likelihood and the consequence severity of the fire. The formula is as follows:where is the magnitude of the fire risk presented by various potential fire hazards, refers to the occurrence likelihood of potential fire hazards or fire risk factors, implies the consequence severity of potential fire hazards or fire risk factors, and represents the interconnection relationship between and .

2.3.1. Fuzzy Linguistic Variables for the Fire Risk

After defining the fire risk, it is necessary to transform the factors into the same form of fuzzy evaluation grade. Due to the uncertainty, analysts tend to use linguistic variable terms rather than precise numerical values to evaluate the fire risk. Therefore, this paper uses a ranking form of fuzzy linguistic variables to represent the fire risk profile of each factor.

A belief degree is generally used to describe the level of expectations for trust events, and it must be less than or equal to 1 to express the degree of which answer is considered true. Individual differences of belief degree depend on assessor’s expertise and the knowledge to understand the assessment system. A belief structure can solve the problems of fuzziness, uncertainty, and imprecision in human decision-making. Therefore, this paper presents a model that combines fuzzy linguistic variables and a belief degree to construct a belief structure with the same set of assessment grades [26]. These sets’ form of each factor could be expressed as follows:

Among them, , , and represent the evaluation grade variables of the occurrence likelihood of fire, consequence severity of fire, and fire risk, respectively.

2.3.2. Fire Risk Level Based on a Belief Structure

Because of the complexity and uncertainty of the system, the type of membership function is not the dominant factor in the risk assessment analysis of the system [27]. Therefore, as listed in Table 2 and Figure 2, this paper applies the triangular membership function which is the most commonly used one to describe the subjective linguistic variables [15] and adopts the five-phase method, adjusted and modified from Ngai and Wat [28] to represent the occurrence likelihood of fire (L) and the consequence severity of building fire (S), respectively. Suppose that the occurrence likelihood of building fire (L) and the consequence severity of building fire (S) for each factor are independent of each other; they are denoted by triangular fuzzy numbers and . Then, the desired fire risk of each factor can be expressed as

Fuzzy risk P with a belief structure can be obtained through the following steps:(1)According to formula (8), calculate of each factor(2)Map the calculated to the membership curve, and obtain the intersection points of each fuzzy language level variable (note: if there is more than one intersection point on a certain fuzzy language level variable, take the intersection point with the largest longitudinal coordinate value), as shown in Figures 3 and 4(3)Obtain a set of intersection values (), which denote five nonstandardized linguistic variable levels of risk P in the form of fuzzy sets(4)Normalize , and obtain the basic belief degree of each factor related to its fire risk

As listed in Table 2, if a single factor judged by experts’ knowledge and experience takes a fire risk value of that the occurrence likelihood of building fire corresponds to (0.5,0.75,1), the consequence severity of building fire corresponds to (0.25,0.50,0.75). The corresponding value of will be (0.125,0.375,0.75). Then, map to to get the set of intersection values (), shown in Figure 4. Finally, the basic belief degree is obtained after the normalization of , which denotes that five nonstandardized linguistic variables of very low, low, general, high, and very high correspond to 0.25, 0.75, 0.8, 0.4, and 0, respectively.

It is noteworthy that the triangular fuzzy numbers for the occurrence likelihood (L) and the consequence severity (S) of building fire judged by experts cannot be used directly as input data for the synthesis of fire risk results by the evidential reasoning algorithm. They need to convert to five standardized linguistic variable terms before synthesizing the fire risk of each factor [29].

2.4. Synthesizing Assessment Result Using the Evidence Reasoning Algorithm

The theory of evidential reasoning was first proposed by Dempsterin 1967 [30]. Then, in 1976, Shafer further expanded and improved Dempster’s work to form a complete and systematic theory [31]. Subsequently, in commemoration of Dempster and Shafer’s contribution to the theory of evidence reasoning, the theory was often called Dempster–Shafer theory or D-S theory for abbreviation. D-S theory can be used to deal with uncertain, imprecise, and or inaccurate information. It was originally used as an approximate reasoning tool for information synthesis in expert systems [32]. Later, it was applied to the decision-making judgment of uncertain problems [33]. Due to the uncertainty of the changing system environment and qualitative descriptive information and to consider the influence of the weight in the synthesis of evidence, evidence reasoning algorithm (ER algorithm) was proposed [34].

After knowing the basic belief degree and the weight of each factor, suppose is a basic probability mass, denoting the degree to which the -th basic factor supports the general factor to be evaluated as the -th grade:

The unassigned probability mass is composed of two parts, which represent the unassigned mass function due to the weight and the unassigned mass function due to the lack of information and incompleteness:

Suppose represent the combined masses of basic factors synthesized on the -th evaluation grade. Suppose represent the unassigned probability mass to the first basic factors. The formula is as follows:where represents the normalizing factor, which reflects the degree of conflict between the indicators (evidence). Suppose that there is a total of basic factors for evaluation objectives; then, , , and are obtained by iteration calculation. After that, the combined belief degree can be obtained by the following normalization process:where represents the unassigned belief degree to the general factor after aggregation. and represent the comprehensive belief degree to the evaluation object.

2.5. Obtaining a Clear Result Using the Expected Utility Method

In fact, the belief degree vector obtained in the former evaluation is the trust distribution of risk under the identification framework, and the result cannot be shown clearly. For example, the identification framework of a building fire risk (i.e., assessment set) is recorded as “very low,” “low,” “general,” “high,” and “very high.” Suppose that the combined degree of the belief vector is (0,0.45,0.5,0.05,0), calculated by the above formula, which means that the construction risk level corresponds to a “low” level of 45%, a “general” level of 50%, and a “high” level of 5%. However, this information cannot clearly indicate the magnitude of the fire risk. Therefore, the concept of utility value is introduced in [35] as follows:where represents the utility of the evaluation grade . In order to further clarify the level of the fire risk corresponding to the utility value, it is necessary to classify the grade of the fire risk. This paper presents the classification as shown in Table 3.

Quantitative evaluation results (utility values) can be obtained by processing the above methods. However, if the basic attribute (factor) information is incomplete or the expert’s information about the factor is uncertain, the result obtained by the ER algorithm is also uncertain. [34, 3638] refer to the concept of utility interval and conquer this problem through minimum utility , maximum utility , and average utility :

2.6. Verification of the Model Using Sensitivity Analysis

Due to the influence of external factors, input values obtained from different experts or the same experts in different periods are different. Consequently, the uncertainty is inherent in fire risk assessment. In this paper, a sensitivity analysis method is introduced for studying and predicting the disturbance degree of the model output value (risk magnitude) caused by the change of the input value of each index. Sensitivity analysis is a systematic analysis method, which identifies weak points or areas in the system with the insight of managers in quantitative evaluation and continuously improves the design of the system and improves the stability of the system [39].

If the validated model is reliable and its reasoning process is logically feasible, then the sensitivity analysis of the model at least satisfies the following three theorems:(1)A slight increase/decrease in the degrees of belief at any linguistic variables of the lowest-level factors will result in increase/decrease in the fire risk level of the output of the model(2)If the belief degree at the lowest preference linguistic variable of the lowest-level factors increases by and (meanwhile, the belief degree at the highest preference linguistic variable decreases by and (1 >  > )) and the utility values of the model output are and , then should be greater than (3)In the lowest-level factors, the total influence of factors on the output of the model is always greater than that of factor sets

3. Case Study

Three residential buildings marked from BUILDING-1 to BUILDING-3 were selected as a case study to illustrate the proposed fire risk model. This paper takes BUILDING-1, for example, to describe the calculation process of the model step by step. Based on the hierarchical structure of the fire risk model in Figure 2 and the available information in [40], the fire risk of BUILDING-1can be assessed through the following steps.

3.1. Develop a Generic Fire Risk Model for BUILDING-1

At this phase, the identified fire risk factors and a generic fire risk model are presented in Figure 1. The index system of fire risk assessment mainly consists of three levels, including the total target risk, the first-level factor set, and the second-level factor set. According to Wang et al. [41], fuzzy linguistic terms for risk expression are used for effective information processing in the range of 4 to 7. Therefore, this study uses five linguistic terms to denote the assessment of fire risk based on the viewpoint of experts in the field.

3.2. Determine the Weights of Each Factor

Given the hierarchical structure of fire risk in Figure 2, the weight calculations for fire risk factors are conducted. The weight calculations of factors U1, U2, U3, U4, and U5 are taken as an example. Firstly, the judgment matrix is constructed through the pairwise comparison of these five factors by experts (according to Table 1) and presented in Table 4. Then, according to equations (3) and (4), the value of synthetic extent and the possibility degree of each factor are obtained, respectively. Finally, the normalized weight vector for five factors is obtained. Using a similar way, the weights of all factors can be calculated and listed in Table 5.

3.3. Application of the Belief Structure for Fire Risk in BUILDING-1

According to the actual situation of BUILDING-1 fire safety, the occurrence likelihood of fire (L) and the consequence severity of fire (S) for each bottom index should be scored, and the scoring standards are mainly based on the code for fire-protection design of buildings (GB50016-2014) [42], code for fire prevention in design of interior decoration of buildings (GB 50222-2017) [43], guidance on building fire risk assessment for property insurance, and CIB W14 Workshop Report [44]. For example, the detailed scoring rules of index U15 (building service life) and index U4 (property fire management) are shown in Table 6. According to these rules, it is easy to obtain the value of and of each bottom index. Accordingly, by utilising equation (8), the fire risk of each bottom index is presented in Table 7 in the form of . Then, is mapped to for obtaining the intersection point. Finally, the basic belief degree is obtained after the normalization of , and the results are shown in Table 8.

3.4. Synthesizing Assessment Result Using the Evidence Reasoning Algorithm

On the premise that the weight of each index was obtained, the aggregation calculations for U11, U12, U13, U14, U15, and U16 were implemented according to the D-S operator (equations (9)–(20)); then, the aggregation result of disaster-causing factor U1 is obtained. Similarly, the aggregation results of passive measures U2, active measures U3, property fire management U4, the rescue capability of fire brigade U5 and the objective fire risk R can also be obtained, and the results of the first-level index are presented in Table 9.

3.5. The Target Fire Risk Assessment Using the Expected Utility Method

From Table 9, the objective fire risk R corresponding to five-level linguistic terms can be expressed as R = {VL (0.36), L (0.43), M (0.16) H (0.04), VH (0.005)}, which cannot reveal the magnitude of the target fire risk in a clear way. Thus, the final fire risk (FR) is evaluated using equation (21), and the result is 0.2205, shown in Table 10. According to Table 3, it could be observed that the objective fire risk R is acceptable, but if cost-effectiveness is reasonable, measures can be taken in this building to reduce its fire risks.

3.6. Sensitivity Analysis

In order to verify the model, the degrees of belief at the lowest preference linguistic variable of the lowest-level factors should increase by 10%, 20%, and 30% (meanwhile, the degrees of belief at the highest preference linguistic variable decrease by 10%, 20%, and 30%). The model output data are tabulated in Table 11, and the graphic display results are listed in Figure 5. It is obvious that all the results are consistent with theorems 1 and 2, respectively. According to theorem 3, if the model is logically reasonable and feasible, the belief degree at the lowest level of the hierarchy structure associated with factors will always be smaller than the one associated with factors. This can be illustrated by comparing the results of different input data, such as if the belief degree at the lowest preference linguistic variable associated with all the lowest-level factors increases by 10% (simultaneously, the one at the highest preference linguistic variable decreases by 10%), the output utility value is 0.1717. However, if the belief degree at the lowest preference linguistic variable associated with U11, U12, U13, U14, U15, U16, U21, U22, U23, U24, U31, U32, and U33 factors increases by 10% (simultaneously, the one at the highest preference linguistic variable decreases by 10%), the output utility value is 0.1826. Considering that 0.1717 is smaller than 0.1826, it can be concluded that the verified model satisfies theorem 3.

4. Results and Discussion

Based on the results of the case study in Table 12 and Figure 6, it can be observed that the fire risk level of three buildings corresponds to “moderate” or below. However, it is noteworthy to mention that some aspects should be paid attention to.

In the aspect of disaster-causing factor U1: since the service life of the three residential buildings is less than 10 years and there is no dangerous disaster-causing factor in the internal and external environment of these buildings, the fire risk corresponding to U1 of these buildings is all acceptable.

In the aspect of passive measures U2: U2 of BUILDING-2 and BUILDING-3 was higher than that of BUILDING-1. This was mainly due to obstruction of safe evacuation in the stairwell of BUILDING-2 and BUILDING-3, such as some evacuation passageways are littered with debris and some safety evacuation signs are missing, which mean that the residents may fail to evacuate from these buildings in case of a fire.

In the aspect of active measures U3 and property fire management U4: U3 and U4 of BUILDING-3 were higher than those of BUILDING-1 and BUILDING-2. This was mainly due to the lack of regular maintenance of fire-fighting equipment in BUILDING-3. It can be assured that if there is no regular maintenance and inspection, the reliability of fire-fighting equipment will be reduced. In BUILDING-3, it was found that some fire-fighting equipment were rusty or even abandoned, such as safety monitoring device was out of use, and the fire extinguisher was out of the service date range. In addition, U4 of BUILDING-2 was higher than that of BUILDING-1. This is mainly because that, in BUILDING-2, there is no prejob training of safety management personnel, and daily fire hazard investigation is not carried out.

In the aspect of the rescue capability of fire brigade U5: BUILDING-2 and BUILDING-3 are all located in CBD of the city. It means that the traffic around the buildings is congested, and the nearby fire brigade may not be able to arrive in time. In particular, BUILDING-2 is further away from the fire brigade than BUILDING-3.

Furthermore, the fire risk of residential buildings is determined by many factors in the complex external environment. It is noteworthy from the analysis that a slight change will lead to the corresponding change in the output value of the model. According to Figure 5, it is obvious that the fire risk model is more sensitive to occupant density U12, electrical equipment U11, property fire management U4, portable fire-extinguisher apparatus U32, and indoor hydrant water supply system U31 than other factors. In other words, the uncertainty of these factors has a relatively large influence on the disturbance of the model system. Therefore, the most effective way to reduce the fire risk of residential buildings is to control these five indicators at first. The analysis results are also consistent with the actual fire prevention measures.

In the previous studies [40, 45], grey correlation method and fuzzy clustering method were applied for fire risk assessment in these buildings of China, and the results of these studies are in accordance with the results of our research, which indicated that the presented model is logically feasible and can still maintain its specific function under turbulence or uncertainty conditions.

5. Conclusions

This study proposes a novel model which combines evidence theory, fuzzy theory, and sensitivity analysis technique for assessing the building fire risk using inaccurate input data in order to optimize system operating efficiency by a standardized fuzzy linguistic term. This model is different from the traditional risk assessment model and characterized with flexible data acquisition capability and unified input and output modes. Therefore, it is easy to deal with the uncertainty of the fire risk problem in the complex system.

Furthermore, the model adopts a series of processes, such as weight calculation based on the FAHP, two-dimensional measurement of the fire risk based on triangular fuzzy numbers, construction of the belief structure, factor aggregation via the evidential reasoning algorithm, and assessment results using the expected utility method, to effectively address uncertainties of subjective estimation. In summary, the proposed model has the following advantages for fire risk analysis on the complex system: (1) this model presents a managerial view to analysts in a reasonable, reliable, and transparent way so that they can collaborate with experts’ suggestion or on-site investigation to model complex systems under external uncertainties. (2) The model provides an effective tool for researchers to make full use of limited information to assess the fire risk of the whole system and improve its operational flexibility. (3) The model has strong flexibility, has high robustness, and is easy to program. It can be used as a computer tool for fire risk assessment of complex systems under high uncertainty.

This research proposes a quantitative fire risk assessment model which could provide building fire managers and researchers with flexible and transparent tools to effectively reduce the fire risk under the disturbance of fire risk uncertainty of the system. It should be noted that, in our study, the index scoring rules are mainly based on codes and standards, which lead to conservative results. Therefore, the acceptable level of fire risk based on performance-based codes needs to be determined in the future [46].

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by the National Natural Science Foundation of China (NSFC, 51704054 and 51874255), Key Technologies for Prevention and Control of Serious and Extraordinary Accidents of Ministry of Emergency Management (no. Chongqing-0001-2018AQ), the Natural Science Foundation of Chongqing (cstc2019jcyj-msxmX0462), and Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJQN201801531 and KJQN201801519).