A Novel Method for Source Tracking of Chemical Gas Leakage: Outlier Mutation Optimization Algorithm

Chemical industrial parks, which act as critical infrastructures in many cities, need to be responsive to chemical gas leakage accidents. Once a chemical gas leakage accident occurs, risks of poisoning, fire, and explosion will follow. In order to meet the primary emergency response demands in chemical gas leakage accidents, source tracking technology of chemical gas leakage has been proposed and evolved. This paper proposes a novel method, Outlier Mutation Optimization (OMO) algorithm, aimed to quickly and accurately track the source of chemical gas leakage. The OMO algorithm introduces a random walk exploration mode and, based on Swarm Intelligence (SI), increases the probability of individual mutation. Compared with other optimization algorithms, the OMO algorithm has the advantages of a wider exploration range and more convergence modes. In the algorithm test session, a series of chemical gas leakage accident application examples with random parameters are first assumed based on the Gaussian plume model; next, the qualitative experiments and analysis of the OMO algorithm are conducted, based on the application example. The test results show that the OMO algorithm with default parameters has superior comprehensive performance, including the extremely high average calculation accuracy: the optimal value, which represents the error between the final objective function value obtained by the optimization algorithm and the ideal value, reaches 2.464e-15 when the number of sensors is 16; 2.356e-13 when the number of sensors is 9; and 5.694e-23 when the number of sensors is 4. There is a satisfactory calculation time: 12.743 s/50 times when the number of sensors is 16; 10.304 s/50 times when the number of sensors is 9; and 8.644 s/50 times when the number of sensors is 4. The analysis of the OMO algorithm’s characteristic parameters proves the flexibility and robustness of this method. In addition, compared with other algorithms, the OMO algorithm can obtain an excellent leakage source tracing result in the application examples of 16, 9 and 4 sensors, and the accuracy exceeds the direct search algorithm, evolutionary algorithm, and other swarm intelligence algorithms.


Introduction
The world's rapid economic development, urbanization, and industrialization make people increasingly dependent on the chemical industry, machinery industry, and electronics industry, etc. However, the storage of raw materials, the production of spare parts, the assembly of machines, and the testing of technology often need to be carried out in specific places. Among the above industrial processes, the chemical industrial park plays an important role in the storage of raw materials and chemical reaction tests [1][2][3][4]. In chemical industrial parks, there are a large number of chemicals, many of which are flammable, explosive, and hazardous gases [5]. Poisonous gas leakage accidents are challenges for various industries like the petrochemical and gas industries, which will result in financial losses and possibly human injuries [6]. In severe cases, large-scale poisoning disasters or internal domino chain explosion effects may even occur [7].
Toxic gas leakage accidents can be very dangerous to health and the environment. For instance, in December 2018, over 40 workers of a factory in Gujarat's Valsad district were hospitalized following leakage of chlorine gas; in March 2017, about 450 girl students were hospitalized in southeast Delhi's Tughlaqabad area following leakage of 'chloromethyl pyridine' chemical [8]. Inhaling an excessive amount of toxic gas will cause irreversible harm or even death. Humans have different tolerances for different toxic gases. Table 1 shows the Acute Exposure Guidelines Levels (AEGLs) standard of three kinds of typical toxic gases when supposed the exposure duration is 1 h [9][10][11][12]. The unit of AEGL standard is parts per million (ppm), NR means not recommended due to insufficient data. For the general population, AEGL-1's airborne concentration will cause notable discomfort but transient and reversible, AEGL-2 s airborne concentration will cause irreversible health effects and even disability, and AEGL-3's airborne concentration will cause threatening health effects or death [9]. It can be seen from Table 1 that different toxic gases have different concentration threshold standards; even small-scale leak accidents may have irreversible hazards.
To minimize harm to humans and avoid further accidents, it is significant to quickly and accurately get information about the leak point in the chemical leakage accident. Tracking the location of the chemical gas leak source and calculating the average leak rate of the chemical gas leak source are the most important tasks in obtaining the leak point information, which is directly related to the safety and accuracy of the fixed-point elimination work of the emergency safety department [6,13,14]. Using a gas dispersion model to match the actual gas sensor concentration distribution is the mainstream for quickly and accurately getting information about the leak point.
In recent years, many studies have focused on the problem of source tracking of chemical leakage. Singh et al. (2015) used a least-squares approach to inverse modeling, aiming to retrieve multiple point releases in the atmosphere [15,16].  used a hybrid genetic algorithm to estimate the source of hazardous material releases [17]. Mario et al. (2018) used a Computational Fluid Dynamics (CFD) model to simulate a large amount of gas leak accident data and proposed deep neural network models as a solution to the problem of leak-source tracking [18]. Omar et al. (2020) used conditional bivariate probability function and series signature to realize the spatial identification and temporal prediction of air pollution sources [19]. Huang et al. (2021) compared similarities between source profiles and ambient air profiles, aiming to improve the tracing efficiency [14]. Yet-Pole et al. (2021) combined the benefits of both the CFD model and Gaussian dispersion model, aiming to provide a model basis for the leak of hazardous materials, etc. [4]. The above-mentioned studies mainly include three aspects of such research: the development of the gas dispersion model, the optimization and improvement of the gas dispersion model, and the combination of the gas dispersion model and an external algorithm. Therefore, the gas dispersion model is the basic requirement for the traceability of chemical gas leakage.
There are five mainstream gas dispersion models: Gaussian plume model, Gaussian puff model, Britter and McQuaid (BM) model, Sutton model, and 3-D Finite Element (FEM3) model. The Gaussian plume model and Gaussian puff model are suitable for a continuous leak accident and instantaneous leak accident, respectively. The advantages of Gaussian models are their applicability to various gases, unlimited size of application scenarios, and appropriate model complexity. In addition, Gaussian models' relatively high consistency between the calculation results and the experimental values make them widely used [4,[20][21][22][23][24]. The BM model is based on a series of test data of instantaneous leak and continuous leak of heavy gas, combined by drawing calculation charts; the empirical model has a narrow application range and poor extension [25]. The Sutton model uses the statistical theory of turbulence diffusion to simulate the leak and diffusion process of gas, and there will be large errors [26]. The FEM3 model uses the finite element method to solve the gas leak process-related equations, such as the ideal gas state equation, Navier-Stokes equation, etc. The model has developed rapidly in recent years; the biggest advantage is that it can accurately solve the gas diffusion problem under complex terrain conditions. The shortcomings are detailed meshing, complex solution process, and long calculation time [4,5,[27][28][29].
The combination of a gas dispersion model and optimization algorithm has the advantage of various combinations, simple principles, and various applicable scenarios, and is thus widely applied in meeting the demands of source tracking. Table 2 shows the classification and characteristics of some commonly used methods for source tracking of gas leakage.  [18] Big Data with Probability Function [19] Sequential Monte Carlo Methods [40] K-Nearest Neighbor Classifier [41] Other methods Advanced technology, high accuracy

High technical and economic requirements
Drone-Enabled Participation [6] Random Walk Robot Participation [42,43] Hybrid Optimization Algorithm [44][45][46][47] It can be concluded from Table 2 that considering the economic level, prior technology level and calculation time, the mainstream solution to the problem of gas leakage traceability still relies on the combination of gas dispersion models and optimization algorithms. However, the accuracy of these methods has not reached a high level, which may affect the efficiency of emergency response. For instance, the accuracy of the least-squares method under the one-dimensional gas leakage model is only about 95% [15,16], the accuracy of the DNN under 25 hidden layers is only 75.43% [18], and the approximate errors of PSO, ACO and FA algorithms under the three-dimensional gas leakage models are (∆x, ∆y, ∆Q) = (43.97%, 49.43%, 19.75%) [30].
In order to trace the source of the gas leak in a satisfactory time and with as high accuracy as possible in chemical gas leakage accidents under flat terrain conditions, this article has made the following innovations:

•
Traditional gas leakage accident application examples usually take the position and number of the storage tank or environmental variables as known conditions. In order to illustrate that the OMO algorithm has the best comprehensive performance in the gas leakage traceability problem under flat terrain conditions, randomness is introduced to the various variables that will affect the calculation results of the Gaussian plume model, which can make the application examples created in this article to be strongly representative.

•
The introduction of sensor information defect conditions: the number of sensors in the three types of application examples is 16,9,4, respectively, and the intensity of the information defect increases in turn, which is a challenge to the exploration integrity of the optimization algorithm, but it is also an important part to prove the advantages of OMO algorithm.

•
The OMO algorithm is different from other swarm intelligence algorithms. The exploration phase of the OMO algorithm is dominated by outliers, and the exploration method introduces the Levi flight, which aims to avoid the local optimum with a high degree of exploration. In the exploitation phase of the OMO algorithm, two-way selection improves the convergence speed and accuracy. The complementarity of the encirclement and mutation strategies has become a major feature of the algorithm.

•
The research on the parameter law of the OMO algorithm and the introduction of precision control parameters enable the OMO algorithm to adjust parameters reasonably in practical applications and quickly adapt to various application scenarios.
The following content is divided into four sections: methodology, experimental design, results and discussion, conclusions.

Gaussian Plume Model
Last century, the Gaussian model was proposed and applied to the diffusion modeling of actual point sources, including the Gaussian plume model and Gaussian puff model [48]. The Gaussian puff model is suitable for the instantaneous gas leak diffusion from point sources, while the Gaussian plume model is suitable for the continuous gas leak diffusion from point sources [49]. In this paper, the location and intensity of the leak source will be inversely calculated based on the Gaussian plume model.
After the leaking gas is fully mixed with air, the diffusion of the leaking gas is mainly affected by air turbulence, forming passive diffusion. When the turbulent field is stable, the spatial concentration of the leaking gas presents an approximately Gaussian distribution of diffusion. The Gaussian plume model is widely used in the physical movement of gas particles [50,51], and it has been proven to have excellent robustness in the case of small-scale gas point source continuous leak accidents [52]. In contrast, some authors have pointed out that although the Gaussian plume model can better reflect the physical characteristics of gas leak diffusion particles, it cannot handle unpredictable and complex changes, such as a gas leak in extreme climates [53].
Therefore, the Gaussian plume model has the following assumptions [54]: • The concentration of leaking gas conforms to Gaussian distribution in the y-axis section and the z-axis section; • The leakage intensity of the leaking gas is continuous and uniform; • The ambient wind speed remains constant during the gas leak process, and the default direction is along the positive x-axis; • The model ignores the effects of any chemical reactions, including sedimentation and decomposition; • The leaked gas follows the ideal gas equation of state and the law of conservation of mass.
When the leak point source is at the origin of the coordinate system, in an infinite space, the expression of the concentration of the leaking gas at a certain position (x, y, z) downstream of the leak source is as follows [55]: where C(x, y, z) is the concentration of the leaking gas at the certain position (x, y, z), unit kg/m 3 ; Q is the leakage intensity of leak point source (leak rate), unit kg/s; u is the average wind speed of the environment, unit m/s; σ y and σ z respectively represent the standard deviation of the diffusion of the leaked gas along with the y and z directions, unit m.
Assuming that the coordinate of the gas leak point source is (x 0 , y 0 , h 0 ), and considering the total reflection of the ground surface, the concentration distribution calculation formula of the Gaussian plume model is finally obtained by using the mirror image method: where all coordinate units in the formula are m. According to the wind speed and solar radiation intensity, the atmospheric stability can be divided into six different grades A-F, as shown in Tables 3 and 4 [56,57]. σ y and σ z are calculated by the empirical formula of the Pasquill-Gifford model [58], which has a certain correlation with the atmospheric stability level. Their calculation expressions are shown in Table 5 (when the leak source position is at the origin of the coordinates).  Table 5. Pasquill-Gifford diffusion coefficient equation.

Outlier Mutation Optimization (OMO) Algorithm
The OMO algorithm is inspired by animal group hunting behavior and evolutionary algorithms, constructed and improved based on human computing habits. The OMO algorithm aims to overcome the shortcomings of inconsistency between the exploration phase and the exploitation phase of most SI algorithms, and reduce the probability of falling into local optimal solutions, finally improving the accuracy of the solution.
The OMO algorithm is mainly divided into two phases: exploration phase and exploitation phase. The phases are alternated according to the escape energy E of the prey: where E is the escape energy of the prey, E 0~U (−1, 1), t is the current number of iterations, T is the total number of iterations set. Figure 1 shows the relationship between the escape energy and the number of iterations in the OMO algorithm. According to the relationship between |E| and 1, the OMO algorithm divides the algorithm into the exploration phase and the exploitation phase. To appropriately increase the step length of the exploration phase and prevent excessive jumps in the exploitation phase, a parameter E para is introduced: In order to adapt to the multi-dimensional exploration interval with different proportions, the step parameters α and β of Levy flight and Brownian motion are introduced respectively: where k α and k β are constants, value 100 in general, S H is the speed control constant of predator's Levy flight, which can be adjusted according to the actual situation, default value 1, S R is the speed control constant of prey's standard Brownian motion, which can be adjusted according to the actual situation, default value 1, LB and UB represent the lower limit and upper limit of the exploration area, respectively. Standard Brownian motion and Levy flight belong to the Markov process. Figure 2 shows a schematic diagram of the trajectory of 500 steps' Levy flights and 500 steps' Brownian motions in a two-dimensional space when the step parameter is 1. Levy flight has the advantage of the wide exploration area, while Brownian motion has the advantages of small space and high density. Combining Levy flight and Brownian motion will enable the OMO algorithm to have a higher degree of exploration and more accurate convergence results. If the continuous stochastic process {B(t), t ≥ 0} defined in the time domain t satisfies the following three properties, it can be called one-dimensional standard Brownian motion (or Wiener process) [59,60]: is considered as the single-step length of standard Brownian motion. OMO algorithm takes the default time interval of 1 s, so the single step length of standard Brownian motion BM~N(0, 1).
When |E| ≥ 1, predators enter the exploration phase and randomly adopt a strategy to explore: where X m (t) represents the average coordinates of predators in the t-generation population, N is the total number of individuals in the predators' population, X i (t) represents the coordinates of the ith predator in the t-generation population, X rabbit (t) represents the position of the target prey (best fitness) in the t-generation population, LB and UB respectively represent the lower limit and upper limit of the exploration area, the symbol '•' is the Hadamard product, D is the coordinate dimension, r 1 (D) and r 2 (D) are 1 × D-dimensional vectors that obey (0, 1) uniform distribution, sign(x) is a sign function, q is a constant that obeys (0, 3) uniform distribution, LF(D) is a 1 × D-dimensional step vector of the Mantegna algorithm simulating the Levy flight, Equations (9)- (12) are the expressions of the one-dimensional Mantegna algorithm, generally, σ ν takes the value 1 and β takes the value 1.5. When q ∈ [2, 3), compared with most SI optimization algorithms, the dimension of random parameters is increased, and the dispersion degree of random distribution is improved. When q ∈ [1, 2) and q ∈ (0, 1), the opposite direction is explored based on the average position of the population and the position of the prey, and E para appropriately increases the unit step length of the predator individual for Levy flight during the exploration phase, for the purpose of jumping out of the local optimal solution and expanding the scope of exploration in the exploration phase.
When |E| < 1, suppose that r is the chance of a predator in mutating (r ≤ 0.5) or not mutating (r > 0.5) before pounce, and then predators enter the exploitation phase: where BM(D) is a 1 × D-dimensional step vector of the standard Brownian motion, r is a constant that obeys (0, 1) uniform distribution, r 1 , r 2 , r 3 are random integers between [1, N] and are not equal to each other, Mu is the coefficient of variation, values between [0.4, 1], H is a 1 × D-dimensional logical vector aimed to increase the diversity of genes in the population, Cr ∈ [0, 1], which means the crossover probability of the individual's genes, and the F function is the objective optimization function. The exploitation phase in the OMO algorithm is divided into two strategies: besiege strategy and mutant strategy. In besiege strategy, the method of combining the small-scale Brownian motion of prey and the shorter step length Levy flight of predators is adopted, which can not only jump out of the local optimal solution but also conduct sufficient small-scale exploration around the prey for high convergence accuracy of the algorithm. In the mutant strategy, predators take mutation and crossover operations; this strategy absorbs the advantages of the integration of exploration and convergence in evolutionary algorithms, guarantees the diversity of the population in the early exploration, and avoids the shortcomings of decreased exploration ability after the exploitation phase is transferred to the exploration phase.
In order to better show the logic of the OMO algorithm, Figure 3 shows a brief flow chart of the OMO algorithm. The OMO algorithm has superior performance in most standard test functions. Table 6 shows the details of the six representative standard test functions [61,62]. The population number N and maximum iterations of the OMO optimizer are set to 30 and 300, respectively. The speed control constant of the predator's Levy flight S H is set to default value 1, and the speed control constant of prey's standard Brownian motion S R is also set to default value 1, other parameters are consistent with those given in Table 6. Table 6. Details of the six representative standard test functions.

Function Name Function Equation Range Optimal Solution and Value
Sphere Figure 4 shows the test results of the OMO algorithm when the dimension is 2. The equations and parameter range of the test functions are consistent with those shown in Table 6. The population number N and the number of iterations of the OMO algorithm are set to 30 and 300, respectively.  The test results of the OMO algorithm are shown in Figure 4. The search history diagram reflects the characteristics of sufficient peripheral exploration degree, high randomness of exploration, and strong convergence ability in the OMO algorithm. The convergence curve reflects the advantage of the OMO algorithm that can effectively jump out of the local optimal solution.
The experiment of chemical gas leakage' source tracking will be detailed in Section 3 and result details of the OMO algorithm will be explained in Section 4.

Experimental Design
For the chemical gas leakage accident under flat terrain, this paper proposes the combination of the OMO algorithm and Gaussian plume model to trace the source of chemical gas leakage. The experiment is mainly to verify the accuracy and efficiency of the OMO algorithm in solving this kind of problem, the influence of different OMO algorithm parameters on the results, and the advantages compared with other optimization algorithms.

Chemical Gas Leakage Accident Application Example
This paper proposes a novel algorithm to quickly and accurately realize the source tracking of chemical gas leakage. It is, therefore, necessary to first establish a chemical gas leakage accident application example in order to verify the application effect of the OMO algorithm.
A single application example cannot prove the universality of the OMO algorithm in practical application scenarios. However, more application examples and respective tests also cannot fully cover all the possibilities. Therefore, the category characteristics of the chemical gas leakage accident scene with flat terrain and stable atmospheric conditions need to be fixed, and the variables need to increase randomness. Experiments based on this condition can strongly prove that the OMO algorithm can be efficiently applied to the chemical gas leakage accident scene with flat terrain and stable atmospheric conditions.
In the chemical gas leakage accident scene with flat terrain and stable atmospheric conditions, the main variables include scene range, the location, and intensity of the gas leakage source, average ambient wind speed, atmospheric stability, and the location distri-bution of gas sensors. Because the entire gas leakage scene can be scaled proportionally, so all application examples belonging to the scene mentioned above can be equivalent to the application examples where the scope of the scene is fixed at 100 × 100 and other variables are random. The atmospheric stability rank has a clear classification and calculation formula, which has little influence on the optimization algorithm and indirectly affects the value of the gas sensors, so this paper does not increase the randomness of the atmospheric stability. Considering the reliability of experimental results and the algorithm evaluation system, multiple groups of independent small-scale tests will be conducted in each experiment. All the obtained optimal values, location values of the gas leakage source and intensity values of the gas leakage source will be averaged respectively, and the results after averaging need to have definite comparison indexes. Therefore, each experiment will keep the intensity and location of the gas leakage source unchanged and add the randomness of the gas sensors distribution and the randomness of the average ambient wind speed to the application example. In addition, although different gas leakage source locations and intensities will form different objective functions, the formed objective functions belong to the same type and have little impact on the comprehensive performance evaluation of the optimization algorithm. Therefore, the test results in this paper only show and analyze one of the representative cases.
Suppose that in a flat terrain chemical industrial park with a range of 100 m × 100 m, a chlorine storage tank at an unknown location has a leakage accident. The average leakage rate (the leakage intensity) is Q = 80 kg/s, the location of the leakage source is (x 0 , y 0 , h 0 ) = (8, 15, 1) m, atmospheric stability level rank D, the average environmental wind speed is u = (10 × r w , 0, 0) m/s, and the plane distribution of the gas concentration sensors are as shown in Figure 5, with their height assumed to be 0 m and a measurement unit of kg/m 3 . Supplementary notes: r w obeys the uniform random distribution between (0, 1), the 16, 9, 4 sensors are randomly distributed in 4 × 4, 3 × 3, 2 × 2 squares formed by hollow circles, respectively. The accuracy of an ordinary commercial chlorine sensor (e.g., ME3-CL2) is 0.1 ppm, which is approximately equal to 3.155e-8 kg/m 3 . For the convenience of calculation, the default sensor accuracy is set to 1.0e-8 kg/m 3 in this article.

Chemical Gas Leakage's Source Tracking Model
In order to simplify the chemical gas leakage's source tracking model, this paper ideally considered that the chlorine storage tanks in the chemical industrial park are made following the Chinese industrial standard of 1 m height. Therefore, when a leakage accident occurs, the parameters that need to be back-calculated are the top view plane position (x 0 , y 0 ) of the leakage source and the average leakage rate Q of the leakage source.
The construction of the chemical gas leakage's source tracking model is based on the Gaussian plume model. Suppose the total number of deployed sensors is NS, the chlorine concentration observation value obtained by deployed sensors is C ob (i), and the theoretical value calculated by the diffusion model is C cal (i), where i is the serial number of the sensors. Then the objective function can be expressed in the form of Equation (18): Suppose the position of the ith sensor is (x i , y i , z i ), then the theoretical value calculated by the Gaussian plume diffusion model can be expressed in the form of Equation (19): According to Table 5, the Pasquill-Gifford diffusion coefficient equation of atmospheric stability level D at the position of the ith sensor is obtained: Based on the chemical gas leakage accident application example and chemical gas leakage's source tracking model, the total number NS of deployed sensors is 16, 9, 4, respectively, the lower limit of the optimization interval is LB = (0, 0, 0) and the upper limit is UB = (100, 100, 200), and the concentration observation values of monitoring points are simulated according to the Gaussian plume model.
In addition, the total number of sensors is set to 16,9,4 because the economic factors of sensor purchasing and placing in actual chemical industrial parks are taken into account, and monitoring points' data with insufficient information are more likely to lead to an increase in the number of locally optimal solutions, thus reflecting the advantages of each optimization algorithm.

Other Details and Pseudocode of OMO Algorithm
All algorithms were implemented under MathWorks MATLAB R2020a (MathWorks.Inc, Natick, MA, USA) on a computer with a Windows 10-64bit professional edition (Microsoft, Redmond, WA, USA), 16 GB RAM, and Intel ® Core™ i5-9400F. The population number N and maximum iterations of all optimizers are set to 30 and 300, respectively. The geometric mean has the advantage of good stability, especially when analyzing the value of the objective function; all results are recorded and calculated based on the geometric mean of optimizers over 50 independent runs. Each independent operation will reset the random variables in the application example to achieve a high degree of randomness in the application example, which aims to comprehensively evaluate the multi-scenario application capabilities of each algorithm. The OMO algorithm's parameters use default value and will be tested in Section 4, and the parameters of other algorithms are the same as the default or recommended value, which are reported in Table 7. Whale Optimization Algorithm (WOA) [65] Convergence constant c 1 [2,0] The pseudocode of the proposed OMO algorithm is reported in Algorithm 1.

Algorithm 1. Pseudocode of OMO algorithm.
Inputs: The population size N, maximum number of iterations T, Speed Control Constant S H , S R Outputs: Data required for drawing plotData, the location of prey preyLocation and its fitness value preyEnergy Initialize the random population X i (i = 1, 2, . . . , N), plotData, preyLocation, preyEnergy and step parameters α, β for (each iteration t) do for (each predator X i ) do Update the escape energy E and E para using Equations (3) and (4) if (|E| ≥ 1) then Exploration phase Update the location vector using Equation (8) if (|E| < 1) then Exploitation phase if (r > 0. 5

) then
Besiege strategy Update the location vector using Equations (13), (14), (16), (17) if (r ≤ 0.5) then Mutant strategy Update the location vector using Equations (13)- (17) Checkup population boundary and correct individual Calculate the fitness values of predator population Set preyLocation and preyEnergy as the location of prey and its fitness based on the minimum value between last prey and current population Update plotData and record them Return plotData, preyLocation, preyEnergy

Results
Please note that the calculation time of all algorithms in this section is obtained by summing after 50 independent runs. For experiments that are not affected by the number of sensors, this paper will only show the test results under the 16 sensor scenarios.

Qualitative Results of OMO Algorithm
Tables 8-10 reveal the qualitative results of the OMO algorithm applied to the application example constructed above. It can be seen that when the number of sensors is 16 and 9, the OMO algorithm shows extremely high accuracy in the problem of chemical gas leakage traceability, and when the number of sensors is 4, that is, when the data is severely insufficient, the OMO algorithm will have some deviation, but there is still a certain degree of accuracy. The qualitative results figures of the OMO algorithm applied to the chemical gas leakage accident application example constructed above are demonstrated in Figure 6. These results include five well-known metrics: convergence behavior, the average fitness of the population, search history, the trajectory of the first predator, and the trajectory of the prey. In addition, the escape energy of the prey is also monitored during iterations. The convergence curve reveals how the optimal value of the objective function varies during optimization. The average fitness of the population monitors how the average fitness of whole predators varies during the process of optimization. The search history distribution diagram shows the history of those positions visited by artificial predators during the calculation. The trajectory of the first predator diagram shows how the first variable of the first predator individual varies during the calculation. The map of the trajectory of the prey monitors how the variable of the prey (best solution's location) varies during the iterations. From the search history of predators, it can be concluded that predators attempt to initially boost the diversification and explore the easily overlooked areas of solution space and then exploit the vicinity of the best locations. The wide exploration area can better reduce the search omissions.
As in the trajectory in Figure 6, we see that the foremost predator's amplitude of these variations covers more than 50% of the solution space while the prey's amplitude of these variations covers about 25% of the solution space. As time passes, the explores amplitudes of the foremost predator and the prey gradually decrease. This point guarantees the transition of the OMO algorithm from the exploratory phase to the exploitative phase. Eventually, the motion patterns of the first predator and the prey become very stable which shows that the OMO algorithm is exploiting the promising regions during the concluding steps.
From the fitness and average fitness curve, it can be observed that the OMO algorithm jumped out of the local optimal solution multiple and reached the global optimal solution about the 50th iteration, but the convergence of the average fitness is reached at about 150th iteration, which means that the OMO algorithm applies the greedy algorithm which is not satisfied with the optimal solution in the previous iteration and tries to jump out of the local optimal solution. Thus, when facing the problem of location optimization with a known optimal value of the objective function, a precision control constant can be introduced to realize the function of stopping iteration after the required precision is reached, thereby greatly shortening the calculation time. The test results of the chemical gas leakage accident application example are shown in Table 11 above. Although the OMO method is mainly used to solve the traceability problem of chemical gas leakage, for the entire chemical gas leakage accident emergency and response system, the prediction of gas concentration and the marking of dangerous thresholds are also an important part of it. Figure 7 shows a preliminary schematic diagram of gas concentration prediction and the mark of risk thresholds.

Relationship between OMO Algorithm and Population Number N
In the case of a fixed number of 300 iterations, the results of the relationship between the OMO algorithm and population number N are shown in Table 12 and Figure 8a.  Taking the Briggs logarithm of the optimal value as the criterion for judging the accuracy of the OMO algorithm, the conclusions can be drawn from Table 12

Relationship between OMO Algorithm and Iterations
In the case of a fixed number of N = 30, the results of the relationship between the OMO algorithm and iterations are shown in Table 13 and Figure 8b. Taking the Briggs logarithm of the optimal value as the criterion for judging the accuracy of the OMO algorithm, the conclusions can be drawn from Table 13 and Figure 8b: • The time-requirement of the OMO algorithm is approximately proportional to the iterations; • The accuracy of the OMO algorithm gradually increases as the iterations increase; • When the iterations <100, the accuracy of the OMO algorithm is significantly improved with the increase of iterations; when the 100 < iterations < 500, the accuracy of the OMO algorithm is slightly improved with the increase of iterations; when the iterations >500, the increased speed of the OMO algorithm's accuracy with the increase of iterations gradually decreases and stabilizes.

Relationship between OMO Algorithm and Speed Control Constant
Given the high number of iterations and the high population number, the results calculated by the OMO algorithm will show high accuracy, which is not inefficient in studying the relationship between the OMO algorithm and speed control constant. Therefore, in the case of fixed numbers of r = 0, 100 iterations, and 20 population numbers, 1,000,000 sets of data are recorded based on the free combination of S H = 1:1:1000 and S R = 1:1:1000. The results of the relationship between the OMO algorithm and speed control constant are shown in Table 14, Figures 9 and 10.   These conclusions can be drawn from • Changing the values of S H and S R will have a certain impact on the accuracy of the OMO algorithm, but random values of S H and S R will have a high probability to get a more satisfactory result. • High-precision OMO algorithm results are often accompanied by suitable exploration steps (influenced by S H ), and low-precision OMO algorithm results are often accompanied by too-large exploration steps. The reason why some S H are too large but can also achieve accuracy of the algorithm is mainly due to the random allocation of super boundary data, but in the chemical gas leakage accident model established by this paper, the best value of S H is 1-60.

•
High-precision OMO algorithm results are often accompanied by suitable exploitation steps (influenced by S R ), and low-precision OMO algorithm results are often accompanied by too large exploitation steps. Logarithmic processing has been performed in the step length parameter of the standard Brownian motion to ensure that the exploitation space is small enough in the final stage and the space has been explored completely; therefore, a small-range change of S R will have little effect on the accuracy of the OMO algorithm. In the chemical gas leakage accident model established by this paper, the best value of S R is 1-100.
Therefore, the default values of the speed control constant for the chemical gas leakage accident model (16 sensors) established by this paper can be decided: S H = 4 and S R = 6. The optimal speed control constant values of different application examples are not necessarily the same and need to be verified by experiments, but most of the speed control constant settings can get a satisfactory result with a high probability.

Robustness Analysis under the Influence of White Gaussian Noise (WGN)
Considering the actual factors of the chemical gas leakage accident application example, we can add White Gaussian Noise (WGN) with different Signal-to-Noise Ratios (SNRs) to the observed value of the gas concentration to simulate the error between the calculated value of the Gaussian plume model and the actual observed value, and analyze the robustness of the OMO algorithm. In the case of fixed numbers of N = 30, iterations = 300, S H = 4, and S R = 6, the results of the relationship between the OMO algorithm's accuracy and SNR are shown in Table 15. From the data in Table 15, it can be found that as the noise increases, the error of the calculation result increases relatively. When SNR > 15, the deviation of the leakage source data calculated by the OMO algorithm is within the allowable range of error, but when SNR < 1, the results do not meet the ideal expectations.
In summary, the OMO algorithm has great robustness and can be congruently applied to chemical gas leakage accident scenarios whose SNR > 15.

Comparative Analysis among Optimization Algorithms
In the case of fixed numbers of N = 30, iterations = 300, S H = 1, and S R = 1 in the OMO algorithm, the comparison results of various optimization algorithms are shown in Tables 16-18. In addition, all optimization algorithms were implemented based on the application example established in Section 3.1. The population number N of other optimizers is set to 30, except for the PS and SA algorithms, which do not have the parameter of population number N, and the number of iterations of other optimizers is set to 300. The basic parameters of other algorithms are the same as the default or recommended value, which are reported in Section 3.3, Table 7.
The PSO algorithm has the fastest calculation speed: 2.223 s/50 times. The OMO algorithm has the lowest geometric mean error: ∆x 0 = 1.971%, ∆y 0 = 0.131%, ∆Q = 1.193%, and the stable highest accuracy, which can be reflected by the geometric mean of the optimal value: 2.464e-15.
The PSO algorithm has the fastest calculation speed: 1.638 s/50 times. The lowest x 0 error, the lowest y 0 error, and the lowest Q error are: ∆x 0 = 3.593% calculated by OMO, ∆y 0 = 0.010% calculated by DE, ∆Q = 0.077% calculated by GA. The OMO algorithm has the stable highest accuracy: 2.356e-13.
The PSO algorithm has the fastest calculation speed: 1.341 s/50 times. The lowest x 0 error, the lowest y 0 error, and the lowest Q error are: ∆x 0 = 23.427% calculated by GA, ∆y 0 = 1.121% calculated by SA, ∆Q = 0.709% calculated by GWO. The OMO algorithm has the stable highest accuracy: 5.694e-23.   Analyzing and comparing the results in Tables 11 and 16, it can be concluded that the calculation speed of PS, SA, PSO, and GA algorithms are significantly better than other algorithms, but the calculation accuracy of PS, SA, and GA algorithms is significantly lower than most algorithms, while the calculation accuracy of PSO algorithms is relatively higher. The DE algorithm has great advantages in calculation time and calculation accuracy, swarm intelligence algorithms such as HHO, GWO, SSA, WOA performance is mediocre in the calculation speed aspects, and the SSA algorithm surpasses other algorithms in accuracy, second only to the OMO and DE algorithms. The OMO algorithm has the highest calculation accuracy among all algorithms, at the expense of calculation time. However, the OMO algorithm, which introduces precision control parameters, achieves a balance between precision and calculation, when the optimal value of the function reaches 1e-10 (only second to the DE algorithm) and the calculation time is only about 6.697s (only second to PSO, GA algorithms, and PS, SA algorithms which belong to the direct search algorithm category, which is known for its computational speed).
Analyzing and comparing the results in Tables 16 to 18, it can be concluded that regardless of the severity of sensor data insufficiency in the chemical gas leakage scene, the OMO algorithm always maintains great advantages in accuracy compared with other algorithms. The geometric average of the optimal value mainly reflects the average accuracy of the algorithm; that is, the stability and accuracy of the OMO algorithm are higher than other optimization algorithms. When the number of sensors is reduced from 16 to 9, the accuracy of some optimization algorithms is severely reduced, mainly because the lack of information leads to an increase in the local optimal value. When the number of sensors is 9, optimization algorithms that perform better include OMO, PSO, DE, HHO and SSA.
When the number of sensors is only 4, the severe data insufficiency leads to dense local optimal values, and each optimization algorithm fails to obtain relatively accurate results. However, from the perspective of the relative error and the average optimal value, the OMO algorithm is more suitable than others for solving the source tracking of chemical gas leakage problems.
In summary, the OMO algorithm has shown extremely high accuracy when facing the application example of chemical gas leakage source tracking. The OMO algorithm introduces precision control parameters that can provide high-precision results while ensuring the calculation speed, thereby meeting the emergency needs in the event of chemical gas leakage accidents.

Discussion
As per results in previous sections, we can recognize that the OMO algorithm shows significantly superior results for the application example of source tracking of chemical gas leakage.
The following experiments, studies, or analyses on the OMO algorithm have been described in this Section: qualitative results of the OMO algorithm in attempting to solve the application example proposed in Section 3, relationship between OMO algorithm and population number N, relationship between OMO algorithm and iterations, relationship between OMO algorithm and speed control constant, robustness analysis under the influence of WGN and comparative analysis among optimization algorithms. The following conclusions or default values can be obtained from these experiments, studies, and analysis: the OMO algorithm has excellent exploration and exploitation capabilities; the larger the number of populations, the higher the accuracy of the OMO algorithm; the larger the iterations, the higher the accuracy of the OMO algorithm; the speed control constant has a certain degree of influence on the accuracy of the OMO algorithm, default S H = 4 and S R = 6; the OMO algorithm has relatively good robustness and can calculate reasonable results when SNR > 15 in solving the application example; the accuracy of the OMO algorithm is far superior to other algorithms, and the OMO algorithm that introduces precision control parameters is better than other algorithms in terms of overall performance.
In future work, the basic emergency safety system can be established. Real-time prediction of chemical gas dispersion will be added into this system and the Gaussian plume model will be replaced by a higher precision model with the guarantee of viable calculation time-requirements. First, we will conduct more experiments based on the CFD model and establish a database as reference data for real chemical leakage accidents. Second, the chemical gas leakage scenario experiments will be conducted to verify the established database. Third, we will use deep learning or other feasible methods aimed to establish a universal gas diffusion empirical model. Finally, we will test the accuracy of the model and apply it to the prediction of gas diffusion and replace the Gaussian plume model with the OMO algorithm.

Conclusions
This paper introduces a novel method for the purpose of chemical gas leakage accident source tracking. The OMO algorithm is mainly divided into two phases. In the exploration phase, Levi's flight and its step size vector are introduced, and some individual outlier search is adopted to improve the search range at the initial phase of the exploration. In the exploitation phase, Brownian motion and its step size vector are introduced, and the idea of individual mutation is also introduced, which enriches the path of convergence, and improves the convergence speed and accuracy. As regards the computational efficiency of the OMO algorithm, the test results show that the OMO algorithm with default parameters has superior comprehensive performance, including the extremely high average calculation accuracy-the optimal value reaches 2.464e-15 when the number of sensors is 16; 2.356e-13 when the number of sensors is 9; and 5.694e-23 when the number of sensors is 4. There is a satisfactory calculation time: 12.743 s/50 times when the number of sensors is 16; 10.304 s/50 times when the number of sensors is 9; and 8.644 s/50 times when the number of sensors is 4.
In recent years, the safety of chemical industrial parks has received great attention, especially in China. In order to actively respond to chemical gas leakage accidents, a complete set of safety emergency systems needs to be rigorously developed, including the source tracking of chemical gas leakage accidents and the real-time prediction of chemical gas dispersion. In addition to these basic emergency needs, there are challenges of source tracking of chemical gas leakage accidents with complex terrain and climatic conditions, the real-time predicting of fire and explosion caused by leaking gas, etc.
The OMO algorithm combined with a Gaussian plume model is proved in this paper to be able to quickly and accurately realize the inverse calculation of gas leakage sources under flat terrain and stable atmospheric conditions. This method can provide technical support for emergency measures in normal chemical industrial parks, and can also provide some help for the future study of chemical gas leakage accidents under complex situations.
The OMO algorithm can not only be applied to chemical gas leakage accidents in chemical industrial parks but also the three-bar truss design problem, tension/compression spring design problem, gas pipeline leakage accidents, pressure vessel design problem, etc. Conclusions can be drawn from the results of Section 4: OMO's flexible algorithm parameters can improve its effect on different application examples to a certain extent, including speed control constant, precision control coefficient, the maximum number of iterations, and population number.