Recent increases in climate-induced natural disasters have amplified the risk of Natech (natural hazard-triggered technological) accidents, particularly in the chemical industry. These emergencies, characterized by their urgency and resource constraints, pose significant challenges for emergency planning. The Functional Resonance Analysis Method (FRAM) offers a systematic approach to enhance emergency response strategies. This study introduces a FRAM-based methodology specifically designed for fuel storage tank farms, structured into four critical stages: understanding, designing, analyzing, and enhancing the response process. It promotes a cycle of continuous improvement. A case study on a seismic Natech incident at a fuel storage facility demonstrates the methodology's effectiveness and its potential to boost the resilience of emergency response systems against Natech challenges.@en

Since the development of the best–worst method (BWM) in 2015, it has become a popular research focus in multi-criteria decision-making. The original optimization problem of the BWM is a nonlinear min–max model that can lead to multiple optimal solutions, while the linear model of the BWM produces a unique solution. The two models need to be solved by optimization software packages. In addition, although the linear model of the BWM can obtain a unique solution, it produces different feasible regions than the nonlinear model of the BWM, and it changes the objective function. This study aims to solve the nonlinear model of the BWM mathematically to obtain the analytical forms of the optimal solutions. First, we transform the original nonlinear model of BWM into an equivalent optimization model driven by the optimally modified comparison vectors. The equivalent BWM provides a solid basis for computing the analytical solutions. Second, for not-fully consistent pairwise comparison systems, we strictly prove that there is only one unique optimal solution with three criteria, and there might be multiple optimal solutions with more than three criteria. We further develop the analytical forms of these unique and multiple optimal solutions and the optimal interval weights. Third, we develop a secondary objective function to select a unique solution for the BWM. The secondary objective function retains all the characteristics of the original nonlinear model of the BWM, and we find the unique solution analytically. Finally, some numerical examples are examined, and a comparative analysis is performed to demonstrate the effectiveness of our analytical solution approach.@en