Next Article in Journal
Response of Plant Species Diversity to Flood Irrigation in the Tarim River Basin, Northwest China
Next Article in Special Issue
Optimal Parking Path Planning and Parking Space Selection Based on the Entropy Power Method and Bayesian Network: A Case Study in an Indoor Parking Lot
Previous Article in Journal
Heat Stress and Histopathological Changes of Vital Organs: A Novel Approach to Assess Climate Resilience in Farm Animals
Previous Article in Special Issue
Study on Freeway Congestion Propagation in Foggy Environment Based on CA-SIR Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Level Programming Model Based on Cooperative Operation Study of Stakeholders in Hazardous Chemical Storage

Business School, University of Shanghai for Science and Technology, Shanghai 200093, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(2), 1221; https://doi.org/10.3390/su15021221
Submission received: 1 December 2022 / Revised: 30 December 2022 / Accepted: 5 January 2023 / Published: 9 January 2023
(This article belongs to the Special Issue Sustainable City Planning and Development: Transport and Land Use)

Abstract

:
Due to the uncertainty of risk occurrence and the severity of accident consequences in the process of hazardous chemical storage, there are many stakeholders involved in the management and supervision of hazardous chemical storage, and their interest appeals are different. On the basis of ensuring storage safety, in order to balance the interests of stakeholders and achieve cooperative operation, a two-level programming model considering the maximization of social welfare and the interests of warehousing enterprises was proposed. First, the upper model mainly refers to the regulatory department represented by the government, including the daily supervision cost, risk loss cost, risk compensation cost, and penalty coefficient formulated by combining various indicators. In the lower model, the comprehensive risk level of the warehouse is determined by the warehouse enterprise. Based on this, the supervision coefficient is determined. Combined with the punishment coefficient, the warehousing operation cost, warehousing supervision cost, and the punishment cost when the accident occurs under different risk levels are determined. The relevant case analysis shows that, compared with the evolutionary game model, the social supervision cost of the upper level and the enterprise cost of the lower level can be reduced by 0.49% and 30.43% respectively. Compared with the traditional improved particle swarm optimization algorithm, the proposed algorithm can reduce the supervision cost of the upper society and the lower enterprise by 0.11% and 7.05%, respectively, thus achieving a better supervision effect at a relatively low cost.

1. Introduction

Due to the special physical and chemical properties of hazardous chemicals, special protective measures must be taken in the process of storage, transportation, production, and waste disposal. Otherwise, once an accident occurs, it will not only cause property loss, casualties, and environmental pollution, but also cause serious negative social impact. It is precisely because of the danger of hazardous chemicals that there are many stakeholders involved in the construction, operation, and supervision of hazardous chemicals storage facilities. Additionally, their interests are different, so there are multiple cooperative game problems. On the other hand, because the two-level programming model has the characteristics of hierarchy, independence, priority, and autonomy, it can effectively balance the interests of multiple stakeholders. Thus, the stakeholders can achieve the goal of cooperative operation under the premise of equality and mutual benefit of multiple stakeholders. This is also a kind of operation model that guarantees the normal operation of the main enterprises’ rights and obligations.
Regarding the research on the risk of hazardous chemicals, scholars have adopted different research methods according to the characteristics of the subject. Tu Yuanyuan et al. used the cause theory of trajectory cross accidents to build a quantitative risk assessment index system in which the improved risk grading index method can more objectively reflect the risk of hazardous chemicals [1]. Zhang Le and Tong Xing analyzed the serious hazardous chemical accidents of 29 cranes in 2010–2019 so as to find a solution to promote safety production management [2]. Yue Baoqiang et al. established a power system risk assessment model for hazardous chemical explosion accidents in specific scenarios [3]. Yang Li et al. used electrostatic discharge to optimize the support vector machine so as to improve the dynamic risk assessment performance of the support vector machine model [4]. Ning Zhou et al. used the fuzzy comprehensive evaluation method to build a secondary evaluation index group of fire resource demand according to different storage methods of hazardous chemicals, which scientifically and effectively improved the fire fighting rate of hazardous chemicals [5]. Xiangcui Liu et al. designed a risk assessment and decision support system for hazardous chemicals transportation, mainly based on the needs of hazardous chemicals transportation risk assessment, path planning, and emergency rescue [6]. Wei Jiang et al., on the other hand, used the HFACS model to study the human–causal relationship of hazardous chemical storage accidents. It was concluded that the divergent claims of different interests in hazardous chemical storage are important causes of accidents [7].
Regarding the research on solving the divergent claims of interest subjects in hazardous chemical storage mostly adopts a game theory approach. Firstly, Wang Wei and Wang Xiaonan built a three-way evolutionary game model to study the impact of changes in fines, regulatory costs, and regulatory success rates on the evolutionary equilibrium, puting forward more targeted policy recommendations for multi-department cooperation [8]. Under the condition of asymmetric information, Xiao Maocai built a tripartite game model of regional government, local government, and hazardous chemical logistics enterprises in hazardous chemical logistics supervision to analyze the relationship of interests of the three in the process of hazardous chemical supervision [9]. Liu Jiaguo et al. introduced the incidence of safety accidents and studied the impact of port logistics enterprises and government regulatory authorities’ strategy choices on the accident results [10]. Secondly, when the number of p-levels increases, the game strategy will change. Yingzuo Zhao et al. discussed the national supervision of hazardous chemical inspection institutions through evolutionary game on the basis of third-party chemical inspection [11]. Wang Wei et al. introduced public supervision into the supervision of hazardous chemical transport industry, proving that the degree of public participation is the decisive factor affecting the strategic choice of hazardous chemical transport enterprises [12]. It can be seen from the above literature that previous studies mostly affected the government’s regulatory strategy through the formulation of game strategies by all parties, reducing regulatory costs and accident incidence. However, there is not an eternal zero-sum game relationship between the stakeholders, and the stakeholders can also coordinate the interests of all parties through a cooperative relationship.
Regarding the research on the promotion of collaborative stakeholder operating models, a review of the existing literature shows that two-level programming model is more widely used. Since the decision behaviors between the upper and lower layers in the two-level programming model are not completely constrained by each other, but they are adjusted in time according to the changes of the other side, one can finally making the model optimal as a whole, which may be applicable to the optimization of such problems. There are many studies on the two-level programming model to seek cooperation when the main objectives are different. For example, Mingtao Ma et al. used the two-level programming model to maximize the benefits of stakeholders sharing energy storage and reduce the operating costs of the hybrid system [13]. Long Yong et al. proposed a two-level programming model for the cooperative operation of relevant stakeholders of the micro grid, which can not only reduce carbon emissions, but also increase the profits of operators [14]. Zhou Ziyu and Jiang Huiyuan established a two-level programming model for cold chain logistics distribution location and path optimization, providing a theoretical basis for multi-objective optimization of the cold chain logistics network [15]. Zheng Bin et al. based on the characteristics of the post-earthquake disaster relief network, established a two-level programming dynamic model with the goal of maximizing the satisfaction of material delivery time at the upper level and maximizing the fairness of material distribution at the lower level [16]. Jianfeng Lu et al. established a two-level programming model for distribution and scheduling of emergency materials in hazardous chemicals storage to reduce the risk of hazardous chemicals storage in the emergency network [17]. Zhou Haixia et al. built a two-level programming model with the goal of minimizing the logistics cost at the upper level and the logistics time at the lower level so as to quickly and accurately transport emergency supplies to the demand point [18].
Therefore, the main contribution of this paper are as follows:
(1)
This study focuses on the regulatory cooperation among the government, third-party regulatory agencies, the public, and hazardous chemical warehousing enterprises. It also analyzes the interest demands of each stakeholder.
(2)
This study proposes a two-level planning economy model to balance the interest demands of all parties as much as possible.
(3)
The model in this paper reduces social cost and enterprise cost to the greatest extent.
In general, game theory is generally used for interactive decision-making, where players choose strategies that benefit them. However, game theory is mostly non-cooperative games, and the two-level programming model makes up for it. Based on the above analysis, the social cost minimization pursued by the government, third-party regulatory agencies, and the public is combined with the overall cost minimization pursued by the hazardous chemicals warehousing enterprises. A two-level programming model, which is a linear programming problem, is established, and risk is introduced into the model in the form of cost [19]. The idea of the paper is divided into three steps as follows:
(1)
Firstly, the government in the upper model formulates the relevant penalty coefficient, and the hazardous chemicals warehousing enterprises in the lower model formulates the supervision cost coefficient.
(2)
Secondly, the hazardous chemicals warehousing enterprises in the lower-level model determine the comprehensive risk level of goods warehousing status according to relevant standards and estimate the probability of risk occurrence.
(3)
Then, the penalty coefficient formulated by the upper-level government will affect the penalty cost of the lower-level enterprises. Through the formulation of the penalty coefficient, the lower-level warehousing enterprises are urged to meet the safety standards in daily supervision.
Based on the above analysis, the specific model structure is shown in Figure 1. Among them, in the upper level planning model, the following can be seen:
Daily supervision cost refers to the expenses incurred by the government and third-party regulatory agencies for the daily supervision of enterprises, which is to fulfill their own supervision responsibilities.
Risk loss cost refers to the cost that the government has to spend and the expected reduction of economic benefits, which is because of the existence of risks.
Risk compensation cost is the monetary compensation for the people who may bear the risk and suffer the injury around the hazardous chemical storage enterprise.
In the lower level planning model, it can be seen:
Warehousing cost refers to the sum of the input of various factors in the warehousing business activities of an enterprise and is expressed in the form of currency.
Warehousing supervision cost refers to the additional supervision cost invested according to the comprehensive risk level of different warehouses on the basis of daily supervision.
Estimated penalty cost is the penalty that the government imposes on the company after the accident according to the risk level of the accident.

2. Method and Model

This paper is divided into an upper and lower-level model. Relevant symbols are explained in Table 1. The penalty coefficient, α r i , formulated by the government and the supervision cost coefficient, β r i , formulated by the enterprise for different comprehensive risk levels of warehousing enterprises are taken as the decision variables to urge warehousing enterprises to reach safety standards in their daily management activities [20].
In view of the complexity of the research problem, the following assumptions are made for the model establishment:
(1)
Considering the complexity of enterprise cost, this paper only studies the cost of enterprise storage, so as to replace the overall cost of the enterprise;
(2)
The government involved in this paper is the general name of all government departments;
(3)
This paper only considers one warehousing enterprise;
(4)
In order to simplify the calculation, the classification of people is not considered in risk compensation;
(5)
Assume that each costing is based on the unit area of hazardous chemical storage.

2.1. Construction of the Upper Level Model

The upper-level model mainly analyzes the daily supervision cost, risk loss cost, and risk compensation cost of the government, the third-party regulatory agencies, and the public in the process of participating in the supervision. The government estimates the risk value according to the comprehensive risk level so as to change the formulation of the penalty coefficient, as well as to affect the penalty cost of the lower level warehousing enterprises when the risk occurs. This is performed to ultimately reduce the probability of the final risk occurrence and, at the same time, minimize the social cost [21]. Specifically:
This is the upper social cost minimization objective function:
min F = r R ( Z r X r W r + U r + G r )
where: Z r is daily supervision costs per unit area of goods in warehouse, r , from the value of the government and the third-party regulatory agency; X r is storage area of goods in warehouse, r ; W r is whether to store goods in warehouse, r ; U r is the average risk loss in warehouse, r , and G r is the risk compensation cost in warehouse, r .
U r is calculated as follows:
U r 1 = E r + F r X r W r
E r = { α r 1 r X r W r , 0 < d r H D r α r 2 r X r W r , H D r < d r J D r α r 3 r X r W r , d r > J D r
U r 2 = F r X r W r
U r = q r U r 1 + ( 1 q r ) U r 2
where, U r 1 represents the loss when the estimated risk occurs in warehouse, r , E r represents the additional loss when the estimated risk occurs in warehouse, r ; F r is the unit fixed cost invested in warehouse, r , to prevent the risk from occurring; α r 1 , α r 2 , α r 3 is the penalty coefficients corresponding to the slight risk level ( 0 , H D r ] , the medium risk level ( H D r , J D r ] , and the serious risk level ( J D r , ) , which are the decision variables of the upper level model [22]. Among them, it is impossible to have a risk of 0 in the actual production activities. Therefore, the case of the comprehensive risk level d r = 0 does not exist. r indicates the additional unit loss incurred when the risk occurs in warehouse, r ; U r 2 represents the fixed cost invested to deal with the risk in warehouse, r ; and q r is used to describe the estimated risk occurrence probability in warehouse, r .
G r is calculated as follows:
G r = O r N r b
O r = | t r t ¯ | / t ¯
where O r is the risk compensation coefficient in warehouse, r ; t r is the estimated risk value of the warehouse, r ; t ¯ is the average estimated risk value of the warehousing enterprise; N r is the number of exposed population within a radius of three kilometers around the warehousing enterprises, and the number of surrounding population is obtained through GIS software; and b is the compensation for each person [23].
Constraint conditions:
W r = { 1 , Storage   of   goods   in   the   warehouse   r 0 , No   storage   of   goods   in   the   warehouse   r
q r [ 0 , 1 ]
α r i > 0 , i = 1 , 2 , 3
Equation (8) indicates the satisfaction state of the condition; Equation (9) indicates the probability range of risk occurrence; and Equation (10) indicates the value range of penalty coefficient.

2.2. Construction of the Lower Level Model

The lower-level model needs to ensure that the sum of storage cost, storage supervision cost, and penalty cost when the estimated risk occurs is minimized, that is, to minimize the overall cost of the hazardous chemicals warehousing enterprises. The objective function is as follows:
min f = r R ( C r X r W r + P r + q r M r )
where: C r is unit storage cost of goods in warehouse, r ; P r is storage supervision cost in warehouse, r ; M r is the estimated penalty cost when the risk occurs in warehouse, r . P r is determined according to the comprehensive risk level, d r , and the enterprise supervision cost coefficient, β r i ; M r is determined according to the comprehensive risk level, d r , and the penalty coefficient α r i of the corresponding level.
The storage supervision costs P r is calculated as follows:
P r = { β r 1 Y r X r W r , 0 < d r H D r β r 2 Y r X r W r , H D r < d r J D r β r 3 Y r X r W r , d r > J D r
where Y r is the unit warehousing supervision cost in warehouse, r , and β r 1 , β r 2 , β r 3 correspond to the enterprise supervision coefficients of slight, medium, and serious levels, respectively.
M r is defined as (13):
M r = { 0 , d r = 0 α r 1 L r X r W r , 0 < d r H D r α r 2 L r X r W r , H D r < d r L D r α r 3 L r X r W r , d r > L D r
Among them, L r is the unit compensation cost that the warehousing enterprise needs to pay when the estimated risk occurs in warehouse, r .
Constraint conditions:
β r i > 0 , i = 1 , 2 , 3
W r X r B r
Equation (14) indicates the value range of supervision coefficient; and Equation (15) indicates that the quantity of hazardous chemicals stored in a warehouse cannot exceed the maximum capacity of the warehouse, B r .

3. Model Solving and Case Analysis

3.1. The Solution Method of Two-Level Programming Model

Based on the analysis of some existing excellent algorithm ideas [24], this study proposes to use the improved adaptive particle swarm optimization algorithm to solve the two-level programming problem. The parameters of the model were set as follows: population size n = 50 , dimension D = 2 , learning factor c 1 = c 2 = 2 , and ω M a x and ω M i n are, respectively, 0.9 and 0.4. The termination condition of the algorithm is the maximum number of iterations, i t e r M a x = 50 , and the convergence accuracy λ α 10 10 ,   λ β 10 10 , where the convergence accuracy is solved as in Equations (16) and (17).
λ α = | α i + 1 α i | | α i |
λ β = | β i + 1 β i | | β i |
Considering the randomness of the initial setting, the program was run independently for 30 times, and the best solution was taken as the approximate optimal solution. Additionally, the main improvement strategies of the algorithm are as follows:
(1)
The disturbance factor is added to the velocity update formula to expand the population search range [25];
(2)
The adaptive weight method is used to balance the global search ability and local improvement ability of PSO;
(3)
Mutation operation is carried out on the global best solution, that is, random disturbance is added to the global optimum to improve the ability of PSO to jump out of the local best solution [26]. Assuming that random variable η follows the standard normal distribution, and its values are greater than 0 and less than 1, and the improved optimal position of particles is:
P g = P g ( 1 + η )
First of all, the upper-level decision maker formulates an appropriate accident penalty coefficient according to the comprehensive risk level of the warehouse storage state. The lower level takes this decision variable as a parameter, substitutes the result into the formulation of supervision cost, obtains the approximate optimal solution within the possible range, and feeds back the approximate optimal solution to the upper level. This cycle is repeated for many times, and the upper and lower levels of the two-level programming are optimized synchronously, and the approximate global optimal solution of the two-level programming model is finally obtained. The specific model solving process is shown in Figure 2 [27].

3.2. Case Analysis

This paper takes Shanghai Beifang Storage and Transportation Group as an example, uses the improved adaptive particle swarm optimization algorithm to solve, and gives the government’s penalty coefficient α i and the enterprise’s supervision coefficient β i for different risk levels in order to minimize the upper and lower objective functions, and then it analyzes the results.

3.2.1. Case Data Survey

Two warehouses of Beifang Group are selected as the research objects, namely, the dust sodium chlorate warehouse (warehouse 1) and the liquid hydrogen peroxide warehouse (warehouse 2), to find a relatively reasonable penalty coefficient and regulatory coefficient to minimize the social cost and the overall cost of the enterprise.
In order to make the risk assessment results accurate and credible, expert survey method and the literature review were used to collect data and information. First, the quantified values of risk impact levels and the related criteria of risk impact degree were determined, as shown in Table 2. Then, based on the probability of risk occurrence and risk impact level, the risk level classification criteria and the level quantification interval were determined, as shown in Table 3 and Table 4. From the relevant criteria in Table 2, Table 3 and Table 4, the quantitative value of risk level of each risk factor of hazardous chemical storage is calculated by using Equation (19).
Among them, seven risk indicators affecting the safety of hazardous chemical storage were selected in a generalized manner by interviewing the relevant personnel of Beifang Group and combining with the General Rules for Storage of Commonly Used Chemical Hazardous Materials. Finally, the risk weights were calculated by applying the AHP rule to each risk factor in turn to determine the comprehensive risk level of each warehouse and complete Table 5 and Table 6 [28].
Among them, since the quantization interval of the level is not a specific quantization value of the risk level, linear interpolation method is used for division, as shown in Equation (19):
Q L = Q L 1 + ( Q L 2 Q L 1 ) ( D D 1 ) ( Q Q 1 ) ( D 2 D 1 ) ( Q 2 Q 1 )
where: Q L is the quantitative value of risk level, and Q L 1 , Q L 2 are the value range of its quantization interval; D is the degree of risk impact, and D 1 , D 2 are the value range of its quantitative value; Q is the probability of risk occurrence; and Q 1 , Q 2 are the value range of its occurrence probability.
In addition, the estimated risk occurrence probability of each warehouse can be determined by calculating the weighted average of each risk factor, and the calculation model is shown in Equation (20):
R L r = q l r A r
where: R L r is the comprehensive risk level in warehouse r ; q l r is a quantitative value of the risk level of a risk factor; and A r is the risk weight value for a risk factor. Other relevant parameters in the investigation case are shown in Table 7, and the solution idea is shown in Figure 3.
On the basis of Table 5 and Table 6, the estimated risk probability of warehouse 1 is 48.73% and that of warehouse 2 is 32.43% by calculating the weighted average of all factors affecting warehouse risk. According to Equation (18), the estimated risk value of warehouse 1 is 1.904 and that of warehouse 2 is 0.998. Assuming that the average estimated risk value of warehouse 1 and warehouse 2 is taken as t ¯ , the value is 1.451. Meanwhile, according to GIS, the number of exposed people within the radius of three kilometers of Beifang Logistics Company is 562. Ignoring the distance between warehouse 1 and warehouse 2, the risk compensation coefficient of warehouse 1 is 0.904, and that of warehouse 2 is 0.312.
In addition, it is particularly emphasized that when the accident level is judged as serious, the enterprise itself is immediately rectified. Therefore, the case when the accident level is serious is not considered.

3.2.2. Case Result Analysis

MATLAB was used to program and solve the adaptive particle swarm optimization algorithm in Section 3. Considering the randomness of the initial setup, the program is run 30 times independently, and the best solution is taken as the approximate global optimal solution.
In addition, to further verify the effectiveness of the bilevel programming model for solving the cooperative operation research of the stakeholders of the hazardous chemical warehouse. Improved particle swarm optimization algorithm and classical evolutionary game method [29] are, respectively, used to compare the behaviors of regulatory entities represented by the government and warehousing enterprises in this paper. The specific numerical results are shown in Table 8, and the comparative analysis is shown in Figure 4. It can be seen that the optimal solutions of the adaptive particle swarm optimization algorithm in this paper are better than other solution methods. Not only the convergence performance of this algorithm is significantly improved, but also the standard deviation of the algorithm is 0 for 30 runs, and the results are very stable.
It can be seen that, in the case of a small degree of social cost optimization, the overall cost of the enterprise can still be significantly reduced. In addition, the method of this research can reduce the government’s penalty coefficient for slight risk level to a large extent while ensuring the reduction of social costs and overall enterprise costs. However, it also increases the penalty for medium risk level. This is because an excessively high penalty factor may not produce a better supervision effect when the risk level is slight. On the contrary, it will bring great pressure to the enterprise operation and lead to unnecessary cost increases. However, when the risk level exceeds a certain range, a low penalty coefficient may also weaken the incentive for enterprises to supervise themselves, which increases the likelihood of risk occurrence. Therefore, a relatively better value should be obtained according to the actual situation so as to achieve the approximate optimal system configuration.
In addition, we can also see that in terms of enterprise supervision coefficient, the method in this paper can make enterprises strengthen their own supervision coefficient and fundamentally reduce the probability of accidents, no matter whether the risk level is slight or medium. At the same time, with the increase in risk level, the regulatory coefficient is also increasing. In conclusion, the results obtained by the cooperative operation model of the stakeholders of hazardous chemical warehousing based on two-level programming are more practical.

4. Conclusions

Once an accident occurs in the hazardous chemical storage link, it will often cause catastrophic consequences. Moreover, there are many stakeholders, and it is difficult to balance the interests of each subject. Therefore, it is of great significance to study the way of cooperative operation among stakeholders. On the basis of summarizing the existing research results, this study analyses the cost composition of the government, third-party regulatory agencies, the public, and the hazardous chemicals warehousing enterprises in the process of production and operation under the premise that the hazardous chemicals warehousing enterprises may have an impact on society. This involves including the daily supervision cost and risk loss represented by the government and the warehousing operation cost, warehousing supervision cost, and the punishment cost when the risk occurs. A two-level programming model is established in which the upper level minimizes the social cost and the lower level minimizes the overall cost of the warehousing enterprise. Then, the improved adaptive particle swarm optimization algorithm is used to solve the model. Finally, the effectiveness of the model is verified by an example. The relevant case analysis shows that, compared with the evolutionary game model, the social supervision cost of the upper level and the enterprise cost of the lower level can be reduced by 0.49% and 30.43%, respectively. Compared with the traditional improved particle swarm optimization algorithm, the proposed algorithm can reduce the supervision cost of the upper society and the lower enterprise by 0.11% and 7.05%, respectively, thus achieving a better supervision effect at a relatively low cost.
However, only one company’s data were validated by the model, so the validation of the same type of companies needs further improvement. Moreover, the relevant cost differences caused by hazardous chemicals with different chemical properties need to be discussed in the subsequent study. This can be performed by adjusting the composition of the cost parameters in the model and the costs incurred in the actual situation. Substituting these into the formula to solve realistic problems can improve the generalizability of the model. Additionally, through the model, government penalty coefficient and enterprise supervision coefficient for practitioners in the government or the enterprise were solved to develop a reasonable regulatory strategy to ensure the safety of hazardous chemical storage. Finally, about the proposed model in this paper, since two-level programming is an NP-hard problem, it needs to design a specific algorithm for the problem studied, so it does not have wide applicability.

Author Contributions

Conceptualization, J.Y.; methodology, B.X.; validation, X.W.; formal analysis, B.X.; investigation, C.Z.; resources, J.Y.; data curation, B.X.; writing—original draft preparation, B.X. and X.W.; writing—review and editing, C.Z.; visualization, B.X.; supervision, X.W.; project administration, C.Z.; funding acquisition, J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by MOE (Ministry of Education in China) Project of Humanities and Social Sciences (22YJAZH131) “Research on mechanism of route coordinated control of commuter traffic at the urban road network in the environment of big data”. The authors, therefore, acknowledge with thanks the financial support of MOE (Ministry of Education in China).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tu, Y.Y.; Li, R.Q.; Wang, H.H.; Ciu, Y.; Liu, Y. Risk Assessment of Road Transportation of Hazardous Chemicals Based on Improved Grading Index Method. J. Hunan Univ. Sci. Technol. Nat. Sci. Ed. 2021, 36, 8–14. [Google Scholar] [CrossRef]
  2. Zhang, L.; Tong, X. Prominent Contradictions and Solutions in the Safety Production Risk Governance: Analysis Based on 29 Hazardous Chemical Accidents. J. Guangzhou Univ. Soc. Sci. Ed. 2021, 20, 54–62. [Google Scholar] [CrossRef]
  3. Yue, B.Q.; Sun, S.J.; Yang, L.C. Risk Assessment Model of Electricity System in Hazardous Chemicals Explosion Accident Based on Bayesian Networks. J. Saf. Sci. Technol. 2021, 17, 155–161. [Google Scholar] [CrossRef]
  4. Li, Y.; Wang, H.; Bai, K.; Chen, S. Dynamic Intelligent Risk Assessment of Hazardous Chemical Warehouse Fire Based on Electrostatic Discharge Method and Improved Support Vector Machine. J. Process Saf. Environ. Prot. 2021, 145, 425–434. [Google Scholar] [CrossRef]
  5. Zhou, N.; Xu, B.; Li, X.; Cui, R.; Liu, X.; Yuan, X.; Zhao, H. An Assessment Model of Fire Resources Demand for Storage of Hazardous Chemicals. J. Process Saf. Prog. 2020, 39, 1–8. [Google Scholar] [CrossRef]
  6. Liu, X.; Huang, H.; Shi, H.; Zhang, T.; Zhao, J. Design of Risk Assessment and Decision Support System for Transportation of Hazardous Chemicals. J. IOP Conf. Ser. Earth Environ. Sci. 2021, 692, 1–7. [Google Scholar] [CrossRef]
  7. Jiang, W.; Han, W.; Zhou, J.; Huang, Z. Analysis of Human Factors Relationship in Hazardous Chemical Storage Accidents. J. Int. J. Environ. Res. Public Health 2020, 17, 6217. [Google Scholar] [CrossRef]
  8. Wang, W.; Wang, X.N. Research on Evolutionary Game of Hazardous Materials Transportation and Supervision Based on Multi-sectoral Collaboration. J. China Saf. Sci. J. 2020, 30, 69–74. [Google Scholar] [CrossRef]
  9. Xiao, M.C. Research on Charging Strategy of Dangerous Goods Transportation under Multi-decision Subject. Master’s Thesis, South China University of Technology, Guangdong, China, 2019. [Google Scholar] [CrossRef]
  10. Liu, J.G.; Wang, J.J.; Zhou, H.; Zhang, Y.Q. Research on Supervision Problems of Port Hazardous Chemicals Base on Security Risk Level. J. Syst. Eng.-Theory Pract. 2018, 38, 1141–1152. [Google Scholar] [CrossRef]
  11. Zhao, Y.Z.; Yuan, W.; Zhao, Z.; Li, W.; Li, H. National Supervision of Third Party Hazardous Chemical Inspection Institution by Evolutionary Game. J. Chem. Eng. Trans. 2021, 71, 1459–1464. [Google Scholar] [CrossRef]
  12. Wang, W. Analysis on Evolution Game of Government and Hazardous Materials Transportation Enterprise Under Public Supervision. J. Chin. J. Syst. Sci. 2020, 30, 92–96. [Google Scholar]
  13. Ma, M.T.; Huang, H.J.; Song, X.; Peña-Mora, F.; Zhang, Z.; Chen, J. Optimal Sizing and Operations of Shared Energy Storage Systems in Distribution Networks: A Bi-Level Programming Approach. J. Appl. Energy 2022, 307, 118170. [Google Scholar] [CrossRef]
  14. Long, Y.; Liu, C.; Wang, Y. A Bi-level Programming Model Cooperatively Operated by the Microgrid Stakeholders. J. Manag. Sci. China 2019, 22, 69–82. [Google Scholar] [CrossRef]
  15. Zhou, Z.Y.; Jiang, H.Y. Multi-objective Optimization of Cold Chain Logistics Network Based on Bi-level Programming Model. J. Logist. Technol. 2020, 39, 65–70. [Google Scholar] [CrossRef]
  16. Zhen, B.; Ma, Z.J.; Zhou, Y.F. Bi-level Model for Dynamic Location-Transportation Problem for Post-earthquake Relief Distribution. J. Syst. Manag. 2017, 26, 326–337. [Google Scholar]
  17. Lu, J.F.; Wang, X.X.; Zhao, J.H. Optimization of Emergency Supplies Scheduling for Hazardous Chemicals Storage Considering Risk. Sustainability 2021, 13, 10718. [Google Scholar] [CrossRef]
  18. Zhou, H.X.; Mei, Y.R.; Lu, F.W. Research on the Optimization Method of Emergency Material Post Transportation Model Based on Bi-level Programming. J. Chin. J. Internet Things 2020, 4, 86–95. [Google Scholar]
  19. Tan, Z.J. Capacity and Toll Choice of An Add-on Toll Road under Various Ownership Regimes. J. Transp. Res. Part E 2012, 6, 1080–1092. [Google Scholar] [CrossRef]
  20. Wang, W.; Zhang, H.G. Transportation Network Optimization of the Hazardous Materials Considering Vehicle Speed Limit Interval. J. Oper. Res. Manag. Sci. 2021, 30, 128–134. [Google Scholar]
  21. Hu, C.Y.; Liu, G.S. The Environmental Bi-Level Programming Models to the Location Decision in Distribution Centers. J. Chin. Manag. Sci. 2007, 15, 59–62. [Google Scholar] [CrossRef]
  22. Liu, J.R. Open Vehicle Routing Optimization for Cold Chain Logistics Distribution Based on Bilevel Programming Model. J. Bull. Sci. Technol. 2019, 35, 209–212. [Google Scholar] [CrossRef]
  23. Hu, J. Research on Solid Hazardous Waste Recovery Path Optimization considering Risk Equity. Southwest Jiaotong Univ. 2013, 18, 123–126. [Google Scholar] [CrossRef]
  24. Song, Y.J.; Zhang, J.T. Multi-objective CuckooSearch Algorithm for Bi-level Programming Problems. J. Oper. Res. Manag. Sci. 2017, 26, 1–10. [Google Scholar]
  25. Zhang, Y.R. Research on Solving Bilevel Programming Model Based on Improved Particle Swarm Optimization Algorithm. Guangxi Univ. 2013, 18, 123–126. [Google Scholar] [CrossRef]
  26. Huang, W.Y. Optimal Bilevel Programming Model Based on Improved Particle Swarm Optimization Algorithm and Its Solution. J. Stat. Decis. 2018, 1, 88–91. [Google Scholar]
  27. Zhao, Z.G.; Wang, W.Q.; Huang, S.Y. Bi-level Programming Problem Based on Improved Particle Swarm Algorithm. J. Comput. Sci. 2013, 11, 115–119. [Google Scholar] [CrossRef]
  28. Ma, G.L. Study on Risk Assessment and Countermeasures of Hazardous Chemicals Storage in Petrochemical Enterprises. J. China Econ. Trade Guide Theor. Ed. 2020, 10, 128–130. [Google Scholar] [CrossRef]
  29. Zhang, T.X. Analysis of Subject Behavior in Safe Production Supervision of Hazardous Chemicals Based on Game Theory. Beijing Univ. Chem. Technol. 2020, 18, 123–126. [Google Scholar] [CrossRef]
Figure 1. Model Structure of two-level programming model for stakeholders in hazardous chemicals storage.
Figure 1. Model Structure of two-level programming model for stakeholders in hazardous chemicals storage.
Sustainability 15 01221 g001
Figure 2. Solution process of two-level programming model.
Figure 2. Solution process of two-level programming model.
Sustainability 15 01221 g002
Figure 3. Case model solution ideas.
Figure 3. Case model solution ideas.
Sustainability 15 01221 g003
Figure 4. Comparison analysis of models among this paper, evolutionary game method, and traditional particle swarm algorithm.
Figure 4. Comparison analysis of models among this paper, evolutionary game method, and traditional particle swarm algorithm.
Sustainability 15 01221 g004
Table 1. Relevant symbols.
Table 1. Relevant symbols.
SymbolsMeaningsRanges
α r i Penalty coefficient α r i > 0
β r i Supervision cost coefficient β r i > 0
O r Risk compensation coefficient in warehouse r O r > 0
W r Storage of goods
in warehouse r or not
W r = { 1 , Storage   of   goods   0 , No   storage   of   goods  
d r Comprehensive risk level { s l i g h t   r i s k   l e v e l , 0 < d r H D r m e d i u m   r i s k   l e v e l , H D r < d r J D r s e r i o u s   r i s k   l e v e l , d r > J D r
R Warehouses set R N + , r R
Z r Daily supervision costs per unit area in warehouse r Z r > 0 , r R
X r Storage area of goods
in warehouse r
X r > 0 , r R
U r Average risk loss
in warehouse r
U r > 0 , r R
G r Risk compensation cost
in warehouse r
G r > 0 , r R
U r 1 Loss when the estimated risk occurs in warehouse r U r 1 > 0 , r R
U r 2 Fixed cost invested to deal with the risk in warehouse r U r 2 > 0 , r R
E r Additional loss when the estimated risk occurs
in warehouse r
E r > 0 , r R
F r Unit fixed cost invested in warehouse r to prevent the risk from occurring F r > 0 , r R
r Additional unit loss incurred when the risk occurs
in warehouse r
r > 0 , r R
q r Estimated risk occurrence probability in warehouse r q r [ 0 , 1 ] , r R
N r Number of exposed populations within a radius
of three kilometers
N r N + , r R
b Compensation
for each person
b 0 , r R
t r Estimated risk value
of warehouse r
t r 0 , r R
t ¯ Average estimated risk value of warehousing enterprise t ¯ 0
C r Unit storage cost of goods
in warehouse r
C r > 0 , r R
P r Storage supervision cost
in warehouse r
P r > 0 , r R
M r Estimated penalty cost when risk occurs in warehouse r M r > 0 , r R
Y r Unit warehousing supervision cost
in warehouse r
Y r > 0 , r R
L r Unit compensation cost L r > 0 , r R
B r Maximum capacity
of warehouse r
B r > 0 , r R
Q L Quantitative value
of risk level
D Degree of risk impact
Q Probability of risk occurrence Q 0
R L r Comprehensive risk level
in warehouse r
q l r A quantitative value
of risk level of a risk factor
A r Risk weight value
for a risk factor
Table 2. Risk occurrence probability and risk impact level table.
Table 2. Risk occurrence probability and risk impact level table.
Risk Probability LevelProbability of Occurrence
Q 1 ,   Q 2
Influence DegreeQuantization Value of Influence Level
D 1 ,   D 2
Risk Impact Level
A(0, 40]Slight(0, 4]1
B(40, 60]Moderate(4, 6]2
C(60, 100]serious(6, 10]3
Table 3. Risk classification standard table.
Table 3. Risk classification standard table.
Risk Probability LevelRisk Level
123
A
B
C
Table 4. Risk level quantification interval table.
Table 4. Risk level quantification interval table.
Risk LevelQuantitative RangeExplain
(0, 1]The risk is small and appropriate action is required
(1, 2]The risks are high and prompt action is needed
(2, 3]The risks are enormous and require immediate action
Table 5. Risk matrix table for warehouse 1.
Table 5. Risk matrix table for warehouse 1.
Risk Factors Probability of Risk OccurrenceRisk ImpactRisk LevelRisk Weighting
A r
Comprehensive Risk Level
R L 1
Quantitative ValuesLevelQuantitative Values
q l 1
Level
Storage location is not reasonable508Serious2.250.1429
Improper temperature and humidity control909Serious2.570.2173
Product deterioration507Serious2.1250.1825
Mixed storage of chemicals of different properties109Serious1.250.2021
Lax control of ignition source507Medium2.1250.1002
Lack of awareness of personnel management407Medium1.250.1223
Improper operation303Slight0.570.0327
Table 6. Risk matrix table for warehouse 2.
Table 6. Risk matrix table for warehouse 2.
Risk Factors Probability of Risk OccurrenceRisk ImpactRisk LevelRisk Weighting
A r
Comprehensive Risk Level
R L 2
Quantitative ValuesLevelQuantitative Values
q l 2
Level
Storage location is not reasonable201Slight0.1250.0901
Improper temperature and humidity control555Medium1.380.1437
Product deterioration505Medium1.250.1638
Mixed storage of chemicals of different properties109Serious1.250.2527
Lax control of ignition source303Slight1.570.0797
Lack of awareness of personnel management405Medium0.50.1522
Improper operation303Slight0.570.1187
Table 7. Value of parameters.
Table 7. Value of parameters.
ParameterSymbolValueUnit
Daily supervision cost per unit of cargo Z r 12RMB/ton
Storage capacity of hazardous chemicals in warehouse r X r 55 (Warehouse 1)
60 (Warehouse 2)
Ton
The additional unit loss incurred when the risk occurs r 120RMB/ton
The fixed cost of dealing with the occurrence of a risk F r 50RMB/ton
Compensation per person b 1000RMB
Unit storage cost of holding goods in warehouse r C r 10 (Warehouse 1)
9 (Warehouse 2)
RMB/ton
Unit storage supervision cost Y r 18RMB/ton
The unit compensation to be paid when the risk is estimated to occur L r 300RMB/ton
Table 8. Optimization results of model in this paper, evolutionary game model, and traditional particle swarm algorithm.
Table 8. Optimization results of model in this paper, evolutionary game model, and traditional particle swarm algorithm.
PlanSocial Cost/RMBOverall Enterprise Cost
/RMB
Slight Risk Level Government Penalty Coefficient (Warehouse 2)Medium Risk Level Government Penalty Coefficient (Warehouse 1)Slight Risk Level Enterprise Supervision Coefficient of (Warehouse 2)Medium Risk Level Enterprise Supervision Coefficient of (Warehouse 1)
Evolutionary game model698,073.728,756.82.21.471.491.41
Improper temperature and humidity control555Medium1.380.1437
Traditional improved particle swarm optimization algorithm695,398.921,524.11.051.470.761.66
Mixed storage of chemicals of different properties109Serious1.250.2527
adaptive particle swarm optimization algorithm694,612.920,006.70.241.821.311.51
Lack of awareness of personnel management405Medium0.50.1522
Improper operation303Slight0.570.1187
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yao, J.; Xie, B.; Wu, X.; Zhang, C. Two-Level Programming Model Based on Cooperative Operation Study of Stakeholders in Hazardous Chemical Storage. Sustainability 2023, 15, 1221. https://doi.org/10.3390/su15021221

AMA Style

Yao J, Xie B, Wu X, Zhang C. Two-Level Programming Model Based on Cooperative Operation Study of Stakeholders in Hazardous Chemical Storage. Sustainability. 2023; 15(2):1221. https://doi.org/10.3390/su15021221

Chicago/Turabian Style

Yao, Jiao, Beibei Xie, Xiurong Wu, and Cong Zhang. 2023. "Two-Level Programming Model Based on Cooperative Operation Study of Stakeholders in Hazardous Chemical Storage" Sustainability 15, no. 2: 1221. https://doi.org/10.3390/su15021221

APA Style

Yao, J., Xie, B., Wu, X., & Zhang, C. (2023). Two-Level Programming Model Based on Cooperative Operation Study of Stakeholders in Hazardous Chemical Storage. Sustainability, 15(2), 1221. https://doi.org/10.3390/su15021221

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop