Abstract
To realize the predictive control of coal preparation quality and ensure that the quality of washing products is close to the minimum coal quality requirement of coal blending to the greatest extent is one of the important means to maximize production and maintain the interests of customers and enterprises. Therefore, the feasibility of introducing the double response surface method and multiobjective genetic algorithm to solve the aforementioned problems is further discussed. By selecting the controllable factors and noise factors affecting the output and determining their respective value levels, the product table method is used to design the robust parameter design test, and the experimental results are obtained, according to the experimental data, the second-order polynomial model of the mean and standard deviation of each response characteristic is established, and the effectiveness of the model is analyzed. Then, the double-response optimization function of each response characteristic is established according to the type of response characteristic. Finally, the corresponding parameter values are solved by multiobjective genetic algorithm. The internal and external surface method is used to design and run 60 tests. Through optimization analysis, the robust parameter settings are 150.68 kpa, 0.18143.73 kpa, and 30%, and the optimal output is ash 8.499%, which yields 69.54%, meeting the requirements of stakeholders. Moreover, compared with the traditional optimization design method, the superiority of the proposed method is verified, which shows that this method is conducive to the transformation of the coal preparation plant from fire-fighting quality management to preventive quality management and provides support for the accurate control and systematic management of the production process of the coal preparation plant.
1. Introduction
Coal preparation is the main technological process for the value-added of coal products. It is the key link for coal enterprises to identify customer needs, meet customer needs, and create enterprise value. Therefore, the qualification and stability of product quality and the washing rate in the coal preparation process directly affect the coal preparation effect. If the quality of washing products fluctuates greatly or is out of control, it will not only bring technical problems to the subsequent coal blending and affect the interests of customers but also cause enterprises to pursue too high quality and reduce the washing rate to meet the requirements of coal blending, so as to bring benefit losses to enterprises. In order to ensure the stability of coal preparation product quality and maximize the washing rate, experts and scholars at home and abroad have conducted a lot of research. As the coal preparation process is the fundamental reason affecting the washing effect, the process method has always been the focus of research. At present, the most commonly used coal preparation processes at home and abroad are heavy medium coal preparation and flotation. As the flotation technology is only suitable for coal slime with small particle size, its application scope is limited. In contrast, heavy medium coal preparation has the advantages of high separation accuracy and strong adaptability to raw coal and has gradually become the mainstream coal preparation technology at present [1]. Therefore, this study selects the heavy medium coal preparation plant as the research object and discusses the method to realize multiobjective stable output under the condition of determining the coal preparation process. In addition to process methods, some experts and scholars have studied the methods to ensure coal quality and washout rate from the perspective of management and evaluation. For example, Li [2] put forward 11 suggestions to ensure coal quality, and Gong and Xu [3] studied the performance evaluation method of quality management of coal preparation plant, both of which have played a certain role in promoting the guarantee of coal quality and the improvement of benefits. Other scholars have studied the setting of process parameters in the production process, such as feeding mode and feeding pressure [4–6], established the mapping relationship between process parameters and quality characteristics through historical data, and determined the value range of each parameter. Other scholars have studied the application of quality management tools in the coal preparation plant, such as Elevli [7] studied the application of control chart in the coal preparation plant and realized the abnormal monitoring and control of coal preparation product quality. Xi Jin et al. [8] used the nonlinear regression method to study the optimal combination of various parameters in the coal blending link, which provides a basis for the production control and optimization of coal blending. However, the existing research on coal quality and yield control methods in coal preparation process mainly focuses on postdetection. Even the research on parameter optimization is also the optimization of single-response characteristics, does not consider the common optimization problem between quality and efficiency at the same time, and ignores the influence of noise factors on the robust output of production system in the optimization process. In addition, the existing studies mainly focus on the independent optimization of single factors and lack of systematic research on the impact of coal preparation production system on coal quality and yield. As the heavy medium coal preparation is a huge and complex production system, involving hundreds of equipment and thousands of processes, if the whole system is optimized and analyzed, a large amount of data needs to be collected, even a large number of tests need to be carried out, and the research cost is too high. Therefore, the key subsystems of the whole production system are selected for analysis. The principle of subsystem optimization indirectly reflects the process of realizing the optimal output of the whole production system. As the core-processing process in the production process of heavy medium coal preparation, cyclone separation is the key control link to achieve stable quality and maximum yield of heavy medium coal preparation. It is also the key subsystem in the production system of heavy medium coal preparation. Therefore, the robust parameter design of multiobjective-oriented cyclone production system has become a problem worthy of study. Experts and scholars at home and abroad have done a lot of research on the optimal design of hydrocyclone. Sun et al. [9] and Liu et al. [10], respectively, used the grey correlation method and RSM (Reynolds stress model) to study the setting of relevant parameters such as volume concentration and flow field of cyclone in heavy medium coal preparation plant. Azizi et al. [11] studied the application of Plackett Burman method in identifying key factors of coal preparation and achieved remarkable results. Chen et al. [12] and Zhang Guo et al. [13] studied the method of parameter control database and linear regression to realize the identification and automatic control of the relationship between parameters s. In addition, many scholars have studied the online optimization and control method of heavy medium coal preparation, the optimal method of coal preparation process design, the optimization method of process yield, and the optimal realization method of cost and quality. [14–20].
However, the aforementioned optimization methods have the following shortcomings: (1) the established optimization function is a function of single-response characteristics, without considering the multiobjective optimization of quality and efficiency; (2) the fitting model established by Plackett Burman experimental design method does not consider the interaction between controllable factors; (3) the existing optimization design methods do not consider the influence of noise on the system output; (4) many existing optimization model data basically come from the historical data of field production, so the representativeness of analysis samples and the effectiveness of data seriously restrict the accuracy of the model. The shortcomings of these studies show that the existing methods cannot determine the value of the best operating parameters that can ensure the robust output of the system under the multiobjective guidance. The study of a system parameter setting is essentially a process of changing the concept of quality management, that is, from the original fire-fighting and passive quality management to preventive and autonomous management. In order to achieve this purpose, many advanced methods have studied by experts and scholars, such as Cheng and Pan [21], Lv et al. [22], and Yan [23], have studied the application of DOE method in the specific production process of machinery manufacturing industry. The significant factors affecting specific output are identified, the optimal parameter setting is confirmed, and the output performance index is greatly improved. Lecureux et al. [24] studied the influence of specific influencing factors on electromagnetic sunscreen through DOE technology, providing data support for the optimization of electromagnetic sunscreen production system. Smith et al. [25] and Moscynski [26] studied the application of DOE technology and method in aerospace field. In terms of precision analysis and control, Yoon et al. [27] applied this method to the parameter optimization analysis of the robot foot control system to realize the optimal parameter configuration of the overall system with multiple inputs and outputs. Doe has also been gradually studied and applied in the fields of chemical industry, service industry, performance evaluation, and cost management [28–30]. For the influence of noise factors on system output, Taguchi [31], a Japanese scholar, first systematically proposed the concept of robust parameter design, which was realized by internal and external table design method and signal-to-noise ratio analysis. Although this method has been widely used, it is questioned by Miller [32], Montgomery [33], Tsui [34], and other scholars because it does not consider the interaction between factors, the experiment is lack of sequence, and the SNR index is unreasonable. As an improvement, Vining and Myers [35, 36] expressed the RPD problem as a constrained optimization problem and proposed dual-response surface methodology (DRSM) to solve this problem. Its basic idea is to fit the two response surfaces of the process mean and variance respectively, and then minimize the variance with the mean target value as the constraint, so as to improve the robustness of the process. For example, Zhao et al. [37] studied the robust parameter design of double-response surface method in the process of solder paste printing, which improved the quality of solder paste thickness and reduced the volatility. He et al. [38] took the tire product design of a rubber company in Tianjin as the object, optimized an important quality feature (rubber hardness) in the tire production process using the double-response surface method and designed the test design software according to the double-response surface method. The test design results found the ratio conditions to achieve the rubber hardness target and reduce the hardness variation. By adjusting the combination of controllable factors, the double-response surface reduces the sensitivity of the system to noise factors and ensures the robustness of the system. However, the existing double-response surface analysis mainly solves the robust parameter design problem of single-response characteristics, and the optimization of production system includes two dimensions of quality and efficiency and multiple-response characteristics.
Therefore, the main contribution of this research is to explore the method of robust parameter design of hydrocyclone production system under the condition of multiresponse characteristics. First, the double-response surface method for the response characteristics of quality and efficiency is studied, and the double-response function of each response characteristic is established. Then the multiobjective genetic algorithm is used to solve the evaluation problem of multiobjective optimization function under constraints, and finally, the robust parameter values under multiobjective conditions are obtained. The superiority of the proposed method is illustrated by comparing the traditional satisfaction function method with the gradient descent method.
Based on the aforementioned discussions, the main contributions are given as follows:(1)In this paper, the dual-response surface method is applied to the robust parameter design of cyclone system for the first time, which provides a new idea for the follow-up research.(2)For the parameters that have been estimated, the optimization algorithm is used to optimize, so that more accurate and fine results can be obtained
2. Double-Response Surface Method
Suppose the response variable is Y and the controllable variable is X1, X2, and XK, repeat the measurement on the response variable Y, and calculate its mean and standard deviation. Taking the mean and standard deviation of the response as output variables, Vining and Myers fit the second-order polynomial model respectively:
As two models are used to optimize the mean and standard deviation, how to select the optimization function to optimize the mean and standard deviation at the same time is a problem that engineers and technicians must solve. For example, when Vining and Myers optimize the quality characteristics of the target type, the optimization method they select is as follows:where ω is the target value, and l and u are the boundary of x value. Formula (2) can also be transformed into
3. Multiobjective Genetic Algorithm
The multiobjective optimization problem can be described as follows:where is the objective function to be optimized; X is the variable to be optimized; lb and ub are lower bound and upper bound of variables, respectively; is the equality constraint of X; and are the inequality constraints for X.
Since the objective functions may be contradictory, that is, the improvement of one objective function may require the reduction of another objective function as a cost. As shown in Figure 1, A1 < B1 and A2 > B2 call such solutions A and B noninferior solutions, or Pareto optimal solutions. The purpose of multiobjective optimization algorithm is to find these optimal solutions.

At present, there are many multiobjective optimization algorithms. Kalyanmoy Deb’s fast nondominated sorting genetic algorithm with elite strategy (NSGA-II) is the most widely used one, and this research will adopt a multiobjective optimization algorithm based on NSGA-II.
4. The Double-Response Surface Method Proposed in This Paper
Although the double-response surface method solves the problems of lack of interaction items and lack of sequence of experiments in the process of robust design, the current research on double-response surface is mainly to analyze one quality characteristic, and now the products delivered to customers need to be comprehensively evaluated from the perspective of multiple dimensions and multiple quality characteristics. So, the robust parameter design of multiquality characteristics (response characteristics) has become an urgent problem to be solved, and the multiobjective genetic algorithm just provides the possibility for the realization of robust parameter design of multiresponse characteristics. The traditional optimization methods such as gradient descent method may fall into local extreme values, and the genetic algorithm has a strong ability in the global optimization problem and can solve the problem of local convergence.
For eye-type characteristic Y1, large-scale characteristic Y2 and small-scale characteristic Y3, the respective double-response optimization functions are as follows:
Therefore, the multi response robust parameter optimization function based on the integration of multiobjective genetic algorithm and double response surface method is shown in formula (6).
Moreover, for the acquisition of test data, the product table is used to design the test, and the response surface test design method is used for the inner table to ensure that the established model can reflect the changes of the surface and the influence of interaction terms. According to the two most unfavorable conditions, all the noise factors affecting the inner table test treatment are adopted the comprehensive error method. In other words, each treatment in the internal table runs once under each most unfavorable condition. Therefore, each test treatment will get two values, and the mean and standard deviation of these two values are used as the output response of this treatment.
5. Case Analysis
5.1. Background and Requirements of Test Design
The cyclone production system of a heavy medium coal preparation plant is a pressurized three product heavy medium cyclone production system, which is composed of two product heavy medium cyclone production systems in series. The first section is a high-density area, the products are gangue and medium clean coal, and the second section is a low-density area. Medium clean coal and some small particle size raw coal enter the system for screening, and the products are medium coal and clean coal. The schematic diagram of production process is shown in Figure 2. At present, domestic and foreign coal preparation enterprises adopt the method of “ensuring one end” for the quality control of coal preparation products, focusing on ensuring the quality of cleaned coal products. Coal quality indexes mainly include ash, calorific value, moisture, and so on. Since calorific value, moisture and other indexes can be derived from ash, clean coal ash is selected as the quality response characteristic of this study. In order to maximize the interests of the enterprise, the clean coal washing out rate, that is, the yield index, is also particularly important. Unilaterally meeting the low ash index may meet the requirements of external customers, but the clean coal yield will be greatly reduced, resulting in damage to the interests of the enterprise. Therefore, based on the recent customer demand and underground coal quality, the coal preparation plant decided to set the target of washed coal ash content of the coal preparation plant at 8.50% and the yield target at no less than 69.50% by carrying out the floating and sinking test of raw coal. The floating and sinking test results are shown in Table 1.

5.2. Test Design Scheme
In order to achieve the set response characteristic y and realize the stability of the system, it is necessary to determine the main controllable factors and noise factors affecting the response characteristics in the cyclone production system, that is, x and u in Figure 1, and reasonably set the value level of each factor in combination with the theory and practical experience of cyclone separation in heavy medium coal preparation plant. After consulting a large number of heavy medium coal preparation data [35], it is known that the cyclone production system includes three parameter types. One is the equipment structure parameters, such as the diameter of feed inlet and discharge outlet, and the other is the operating parameters, such as feed pressure and slime content. The other is fixed parameters, such as the washability of feed coal. Since the structural parameters of the equipment and the mechanism of coal are not determined by the production process, the following part will focus on the robust design of the operating parameters of the cyclone production system under the existing production conditions. Combined with relevant heavy medium coal preparation data, national and industrial standards, and with on-site management, the technicians decided to select the feed pressure (x1) in the low-density area, the solid-liquid ratio (x2) in the low-density area, the feed pressure (x3) in the high-density area and the suspended slurry content in the high-density area (x4) four parameters are used as important controllable factors, and the value level of each factor is set. Since the density of heavy medium will change in the separation process and cannot be controlled in a short time, two parameters of suspension density (U1) in low-density area and suspension density (U2) in high-density area are selected as noise factors, and the upper and lower deviations are allowed to be 5% respectively, as shown in Table 2.
The internal and external table method is used to design the test scheme. The response surface test design method is used to design the 4-factor 2-level test in the inner table. The center points in the test scheme are set as 6 and run 30 times in total. The surface is designed by the most unfavorable comprehensive error method, that is, each treatment in the internal table scheme runs once when the noise factor is the maximum value and the minimum value (one test). Therefore, the robust parameter design of the hydrocyclone production system requires a total of 60 tests. The test arrangement and test results are shown in Table 3, where 1 represents the high-level value, −1 represents the low-level value, and 0 is the central point value.
Neutralization in the table is the test results of ash characteristics under the two most unfavorable conditions, which is the mean value of the two test results, σ_Y1 is the standard deviation of the two most unfavorable cases and σ_Y2 is the test results and relevant statistical values of yield characteristics.
5.3. Robust Parameter Design of Cyclone Production System
5.3.1. Model Fitting
According to the test design plan and test results, the double-response surface model including quadratic terms for ash characteristics and yield characteristics is established using the double-response surface method, and the effectiveness of the model is analyzed.
Y1 characteristics, the variance analysis results of the mean obtained by fitting analysis are shown in Table 4, and the estimated effects and coefficients are shown in Table 5.
It can be seen from Tables 4 and 5 that the of the model is 0.02 < 0.05, so it is considered that the model is effective, and each estimated effect is significant. It is further verified that the residual error of the model meets the requirements, and there is no abnormality, as shown in Figure 3. Therefore, it can be considered that the fitting model of uy1 established is effective.

According to the aforementioned analysis, the fitting model of the mean value of Y1 characteristic is as follows:
Similarly, the fitting models obtained are as follows:
According to the double-response surface method in part 2, the optimization functions are constructed for the mean and standard deviation models of Y1 and Y2 response characteristics respectively. The optimization function of Y1 (target type) is shown in formula (9). Since the model established by the model fitting part is obtained based on coding, the value range of controllable factor is [−1, 1]. Since Y1 is the target characteristic and is strictly not allowed to exceed the set standard by 8.50%.
Similarly, the optimization function of Y2 characteristic (expected to be large) can be obtained as follows:
According to the robust parameter design method integrating double-response surface and multiobjective genetic algorithm proposed in the third part of the paper, the optimization function of multiresponse characteristics can be obtained.
The number is as follows:
Since the optimization function (11) includes nonlinear constraints, the external penalty function method is used to transfer the nonlinear inequality constraints to the objective function to obtain the final optimization function, as shown in formula (12), where m is a very large penalty coefficient, which is also the fitness function expression of multiobjective genetic algorithm.
The parameters of the multiobjective genetic algorithm are set as follows: the optimal front-end individual coefficient Pareto fraction: 0.3, population size: 100, maximum iteration generations: 200, stop algebra limit: 200, fitness function value deviation function: 1e − 100. The obtained Pareto front-end distribution and the distance between individuals in the evolution process are shown in Figure 4.

According to the Pareto front-end distribution in Figure 4, the junction distribution is relatively uniform, indicating that the obtained solution is effective and the effect is good. Meanwhile, the other three subgraphs in the figure show that the distance and propagation speed between individuals are gradually decreasing, which indicates that with the increase of evolution, the diversity of the population gradually decreases and the solution gradually converges. Since the optimal front-end individual coefficient is 0.3 and the number of populations is 100, 30 groups of optimal solutions are obtained, of which 4 groups are close to or meet the requirements, as shown in Table 6.
As the Y1 characteristic is required to be strictly less than or equal to 8.50%, only the third group solution in Table 6 strictly meets the requirements. Therefore, the robust parameter design value (code value) of the hydrocyclone production system is (X1: 0.050332, X2: −0.772998, X3: −0.295515, X4: 0.999956), and the actual values are 150.50 kpa, 0.176143.52 kpa, and 30%, respectively. The optimal response characteristics that can be realized are Y1: 8.499% and Y2: 69.59%.
5.4. Comparative Analysis of the Effectiveness of the Method
The optimization function of each output characteristic established by the double-response surface method actually adopts the method of transforming constraints into objective functions, so as to transform them into unconstrained single-objective optimization problems. At present, most of the optimization functions with multiple output characteristics are transformed by the satisfaction function method. After transforming the multiobjective function into a single-objective function, the stability point of the system is obtained by derivation, the type of stability point is determined by canonical analysis, and the confidence interval of the optimal value is estimated according to the normal hypothesis. Therefore, this part shows the effectiveness and efficiency of this method by comparing the differences between the traditional optimization method and the proposed optimization method.
First, if each output characteristic optimizes its dual-response optimization function respectively, set the upper limit of uy1 to 8.500001%, the target to 8.50%, and the lower limit to 8.40%. The type is the desired type, and set σY1 is expected to be large. The response optimizer of Minitab is used to optimize the two functions at the same time. The results are shown in Figure 5.

Similarly, the lower limit of uy2 is set to 69.5%. In order to achieve large scale, the double-response optimization result of Y2 is obtained, as shown in Figure 6.

Considering two characteristics and four response functions at the same time, the results are shown in Figure 7, It can be seen from the figure that the distribution of two characteristics and four response function is relatively uniform, thus demonstrating the effectiveness of the proposed method. The optimization results of the proposed method and the aforementioned analysis results are summarized and Table 7 is obtained. It can be seen from the table that better output results can be obtained when a single output characteristic is optimized respectively. However, the calculated parameter value (factor values) are quite different, indicating that their respective optimizations cannot guarantee the overall optimization of the system. When considering two output characteristics at the same time, the traditional method is used to optimize four response functions at the same time. Although the standard deviation of Y2 characteristics is twice as small as that of the proposed method, the performance of the traditional method is not as good as that of the proposed method in the other three indexes, especially for Y1 ash characteristics, the optimization results are not only closer to the ash requirement of 8.50%, but also have less fluctuation, so as to ensure that the coal preparation plant has high process control ability. Moreover, the maximum yield obtained by the traditional method is 68.4146%, which does not meet the design requirement of 69.50% in the floating and sinking test. In conclusion, it can be seen that the method proposed in this study is better than the traditional multiresponse robust parameter design method as a whole.

6. Conclusion
This paper shows the feasibility of integrating double-response surface method and multiobjective genetic algorithm to solve the problem of robust parameter design under the condition of multiresponse characteristics of cyclone production system in heavy medium coal preparation plant is discussed. The conclusions are as follows: (1) the optimal value of response characteristic ash is 8.499%, the optimal value of yield is 69.59%, and the optimal controllable factors are 150.50 kpa, 0.176143.52 kpa, and 30%; (2) A double-response optimization model between response characteristics and influencing factors is established; (3) The effectiveness of double-response surface method and multiobjective genetic algorithm in the optimization process of production system of heavy medium coal preparation plant is verified, and the qualified and stable coal quality and maximum yield are realized; (4) It has changed the fire-fighting and passive quality management mode of coal preparation enterprises, and gradually changed to the preventive, autonomous, and stable quality management mode. The enterprises are no longer tired of dealing with the continuous coal quality problems, but realize the expected output by controlling the stability of influencing factors. When the whole system is gradually stable, the coal quality inspection link of coal preparation plant can even be cancelled or weakened; (5) Coal preparation enterprises are more confident to realize the accurate control of the production process. In conclusion, the proposed method can explore the relationship between system input and output, reduce the influence of noise on model output and determine the robust parameter values under the condition of multiresponse characteristics, thus providing scientific support for the transformation of quality management mode.
Data Availability
The data set can be accessed upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.