Coal and Gangue Underground Pneumatic Separation Effect Evaluation Influenced by Different Airflow Directions

Coal and gangue underground pneumatic separation is of key importance for green mining. Two kinds of arrangement schemes for high-pressure value used in pneumatic separation system are proposed in this study. Pneumatic separation effects are examined under different arrangement of high-pressure value. Here, theoretical pneumatic separation distance formulas of mineral particles affected by different airflow directions are derived and validated by a series of numerical simulations and orthogonal experiments. In the following analysis, the effects of gangue diameter (d), conveyor velocity (V 0 ), and the height difference between conveyor belt and air nozzle (h p ) are mainly considered. The numerical simulation and experimental results indicate that pneumatic separation effects under the condition of u and V 0 being in the opposite directionwill be better than that of u and V 0 being in the same direction. The pneumatic separation distance ΔS shows a decreasing trend with the increasing of the three factors. The study also shows that gangue diameter has the most significant influence on separation distance, followed by conveyor velocity V 0 and height difference


Introduction
Coal has been regarded as the most widely used fuel in many countries because of the lack of oil and gas resources, which is expected to be prolonged over decades.Initial underground separation for coal and gangue is an indispensable process of mine production and has great significance for coal preparation.Gangue from underground separation can be used for filling mined-out areas, improving the quality of raw coal, decreasing the cost of preparation, saving transport capacity, and reducing the amount of waste polluting the environment above the ground [1,2].Currently, several technologies of underground separation coal and gangue are proposed, which include an automated hydraulic separation technology and a rotating impact method of coal gangue separation underground [3][4][5][6].However, underground separation technologies mentioned above have the disadvantages of high cost, being time-consuming, and low automation.
With the development of digital image processing technology, scholars mainly probed into a general machine vision approach for online identification of coal and gangue [7][8][9][10].The influences of camera location, particle overlap, image blur degree, and conveyor velocity were investigated during the coal and gangue particles online transportation [11][12][13][14].From the above discussion, machine vision process has the advantages of real-time processing, lower cost, and higher intelligence.
Lots of researches have conducted the online identification of coal and gangue, which have proved the feasibility of image recognition of coal and gangue.Nevertheless, pneumatic separation of coal and gangue has not been discussed adequately and has not been applied in mineral separation.
Currently, pneumatic separation has been widely used in agriculture field and rare metal recycling field.A cotton seed online separation system based on machine vision was designed to realize automatic separation of red-brown and black cotton seed by compressed airflow in the air [15].Guo et al. proposed a series of methods to realize the recycling of printed circuit boards (PCBs); the whole process that 2 Advances in Materials Science and Engineering involves crushing and electrostatic and pneumatic separation has formed a closed cycle that can return material and provide salable product [16].Xu et al. described an effective mechanical process including impact crushing from printed circuit boards and investigated the pneumatic separation for metal recovery scraps [17].Kumar et al. report a simple and eco-friendly physical pneumatic separation process for the recycling of metallic values from PCBs [18,19].The effect of orifices on spherical particles was clarified by Hayashi and Oki through numerical simulations of air-solid multiphase flow in a vertical single-column pneumatic separator to realize recycle important rare metals (such as tantalum) in the PCBs of waste electronic equipment [20].Havlik et al. focused on studying mechanical-physical pretreatment of PCBs from used consumer equipment followed by extraction of copper and tin from residue fractions by leaching in hydrochloric acid solutions [21].
Moreover, researchers also have done lots of research on pneumatic separation of coal particles (particle size, i.e., mainly less than 5 mm) in fluidized conveyor.Liu et al. investigated the pneumatic classification behavior in a laboratory CCMC reactor with such a configuration by removing the coal fraction below a given size (e.g., 3.0 mm) from 0 to 20.0 mm [22,23].Yang et al. attempted to carry out a systematic process analysis of fine coal preparation in a vibrated gas-fluidized bed (VGFB) [24].
The research group has proved the feasibility of pneumatic separation technology for coal and gangue, yet pneumatic separation effects influenced by different arrangements of high-pressure value have not been presented.In this paper, theoretical models were established for the study of pneumatic separation process affected by different high-pressure airflow direction.Then, a series of numerical simulations and orthogonal experiments were conducted to evaluate the separation effect under different airflow direction.Finally, the theoretical models were corrected based on experiment data.

Principle
Underground coal and gangue pneumatic separation system, which is shown in Figure 1, contains three parts: preliminary crush separation system, machine version system, and pneumatic separation system.Preliminary crush separation system mainly contains roll-type crusher, spiral size screen, and several conveyor belts, which aims to limit the granularity of coal and gangue sent into machine vision system within a range of 50 mm∼100 mm.Machine vision system is a highly automated technology using digital image processing technology to identify the size and location of gangue based on identification algorithms; then the identified information will be sent to pneumatic separation system through image sensor.Pneumatic separation system is controlled by computer to realize coal and gangue pneumatic separation.
Two kinds of arrangement schemes for high-pressure value used in pneumatic separation system are presented in Figure 1.One is that particles' motion direction is the same as that of high-pressure airflow; the other is that particles' motion direction is opposite to that of high-pressure airflow.As shown in Figure 1, the high-pressure values 12(a) are  arranged under the conveyor belt (7) in the first arrangement scheme, mixed mineral materials are thrown from conveyor belt (7), coal materials unaffected by high-pressure airflow will fall onto conveyor belt (10), and gangue materials affected by high-pressure airflow will fall onto conveyor belt (11).In the second arrangement scheme, high-pressure values 12(b) are arranged in front of conveyor belt (7), mixed mineral materials are thrown from conveyor belt (7), coal materials unaffected by high-pressure airflow will fall on conveyor belt (11), and gangue materials affected by high-pressure airflow will fall on conveyor belt (10).

Theoretical Model.
The movement of mineral materials affected by high-pressure airflow can be divided into three stages: (1) particle's movement before falling into airflow domain; (2) movement of particle in airflow domain; (3) movement of particle after leaving the airflow domain.When the movement direction of mineral materials is the same as that of high-pressure airflow, coal particle trajectory unaffected by high-pressure airflow is shown in Figure 2(a) of curve 1 and gangue particle trajectory affected by highpressure airflow is shown in Figure 2(a) of curve 2. When the movement direction of mineral materials is opposite to that of high-pressure airflow, coal particle trajectory unaffected by high-pressure airflow is shown in Figure 2(b) of curve 1 and gangue particle trajectory affected by high-pressure airflow is shown in Figure 2(b) of curve 2.
As shown in Figure 2, Δ represents the pneumatic separation distance of coal and gangue, ℎ  represents the height difference between gangue mass center on conveyor belt and upper boundary of the airflow domain,   represents the motion time before gangue came into the airflow domain, V 1,2 is the velocity in direction  of gangue before falling  into airflow domain, and ℎ  represents the height of airflow domain.  represents the displacement in  direction of gangue before entering airflow domain. 1,2 represents the motion time of gangue in airflow domain, and   represents the motion displacement in  direction of gangue in airflow domain.V 1,2 represents the motion velocity in  direction of gangue in airflow domain, ẋ ( 1,2 ) represents the motion velocity in  direction of gangue in airflow domain,  1,2 represents the motion time of gangue after leaving airflow domain, and  1,2 represents the motion displacement in  direction of gangue after leaving airflow domain.In the following equations, subscript 1 is used to represent when airflow velocity and conveyor belt velocity are in the same direction, and subscript 2 is used to represent when airflow velocity and conveyor belt velocity are in the opposite direction.
Coal and gangue materials will be recognized by machine vision system and the recognized image information will be sent to pneumatic separation system through image sensor.Coal particles will fall from the conveyor belt freely without the effect of high-pressure airflow.When the motion direction of mineral materials is the same as that of highpressure airflow, gangue material will change its trajectory and will be thrown significantly farther than that of coal particle.When the motion direction of mineral materials is opposite to that of high-pressure airflow, gangue material will change its trajectory to an opposite direction compared with that of coal particle.The motions of coal and gangue particles are analyzed, respectively, in the following three sections.
When coal materials are recognized by machine vision system, coal particles will fall from the conveyor belt freely unaffected by high-pressure airflow. 0 and  0 can be obtained through (1) and (2).Gangue motion law at this stage is the same as coal materials.  and V  can be also shown by ( 3) and ( 4), where  0 represents the motion time of coal particle, and  0 represents the displacement of coal particle in  direction: In order to calculate the acting force of high pressure on gangue particles, multiple rectangular polyline is used to approximate gangue particle's physical shape if the dynamic pressure head of the airflow domain is known in special location.The acting force of high-pressure airflow on gangue can be presented by where  represents the acting force of the high airflow on gangue,   represents the high pressure in ring ,   represents the area of ring , and  represents the total number of rings.As shown in Figure 3, gangue is affected by airflow with different directions.In order to analyze the movement of gangue particle in high-pressure airflow domain, air resistance and horizontal momentum's increment of airflow domain are ignored.The theoretical formula of gangue displacement in direction  can be expressed as  = ∫  1,2 0 (V  +  1,2 ) = ℎ − ; then  1,2 can be obtained as shown in As shown in (6), gangue affected by different airflow direction has the same motion time in airflow domain.Gangue will be affected by two forces after falling into highpressure airflow domain.The two forces applied on gangue Advances in Materials Science and Engineering particle are gravity and high dynamic pressure.Airflow dynamic pressure will convert to static pressure on condition that the velocity of airflow  is larger than gangue's horizontal velocity V 0 and airflow must keep dynamic pressure   =  ẋ 2 /2.The pressure difference can be expressed as Δ  = ( 2 − ẋ 2 )/2.According to (5), the formula of airflow force can be expressed by When the curvature of gangue surface is not too large, it can be approximately seen as a plane, so ∀  = , (7) can be changed to According to Newton's second law, (8) can be changed to ẋ and  are in same direction, ẋ and  are in opposite direction.
As shown in ( 9), the expressions are corresponding to different airflow directions.The general solution of displacement in  direction of gangue moving in the airflow can be obtained by taking Laplace transform to (9), and the result is shown in ẋ and  are in same direction, ẋ and  are in opposite direction. ( The velocity of gangue in  direction in airflow domain can be obtained by taking derivation of (10).The general solution is shown in ẋ and  are in opposite direction. ( In order to calculate the constant coefficients  1 ,  2 ,  3 , and  4 in (10), the initial constraint conditions are as follows: the initial position of the gangue is (0) = 0 and the initial velocity of gangue is   (0) = V 0 .
When ẋ and  are in the same direction, the general solution of (  ) and   (  ) can be obtained by (10) and (11). 1 and  2 are shown in The general solutions of (  ) and   (  ) are shown in When ẋ and  are in the opposite direction, the general solution of (  ) and   (  ) can be obtained by (10) and (11). 3 and  4 are shown by The general solutions of (  ) and   (  ) are shown by ( 15) and ( 16), respectively, ) .
When gangue particles left the boundary of airflow domain, gangue particles would do flat parabolic motion.The formula of velocity V  in direction  is shown in 1,2 can be solved through (18) and the result is shown by (19): The displacement  1,2 of gangue in direction  after leaving the airflow domain can be expressed by To sum up, when ẋ and  are in the same direction, the separation distance Δ can be calculated by When ẋ and  are in the opposite direction, the separation distance Δ can be obtained by The separation distance Δ is shown by ( 21) and ( 22), which reflects the basic motion law of gangue when influenced by airflow field under different airflow directions.For the limitation of assumption, there exist big differences between theoretical model and practical model.In order to simplify the calculation and correct the difference between theoretical value and practical value, parameters   (nonlinear correction factor) and   (linear correction factor) are introduced.  reflects the convergence rate of the fitting function, and   reflects the convergent gain and is used to adjust the fitting effect of formula based on experimental value.Theoretical value will get close to experiment value by adjusting the value of   and   .The modified formulas with the correction factors are shown in Advances in Materials Science and Engineering

Air-Solid Multiphase Pneumatic Separation Simulation.
In order to evaluate the pneumatic separation effect under two kinds of arrangement schemes for high-pressure value, a "fixed coarse-grid" fluid scheme is applied in PFC3D for pneumatic separation simulations.In the fluid scheme, 550 ( In Figure 4(a), the velocity V 0 of conveyor belt 1 and airflow velocity  are in opposite directions.When gangue falls into the airflow domain, gangue particle will change its trajectory to an opposite direction and fall onto conveyor belt 2. The conveyor belt 2 has the same movement direction as airflow velocity; finally gangue particles will be transported to collecting box by conveyor belt 2. In Figure 4(b), conveyor belt V 0 and airflow velocity  are in the same direction.When gangue falls into the airflow domain, gangue particle will throw much farther and fall on conveyor belt 2. The motion of conveyor belt 2 has the same motion direction as conveyor belt 1; finally the gangue particles are transported to collecting box by conveyor belt 2.
In the fluid scheme, a pneumatic boundary is set for two available models.During the simulations, an approximation is made by specifying the velocity boundary at the right end of the model and a pressure boundary as 0.0 Pa at left end with  = 0mm.The slip boundary, in which the fluid velocity parallel to the wall surface is nonzero at the wall surface, is specified at the surrounding four walls.In simulation as shown by Figure 4(b), an approximation is made by specifying the velocity boundary at the left end of model and a pressure boundary as 0.0 Pa at the right end with  = 7000 mm.The initial setup of the slip boundary in Figure 4(b) is the same as that shown in Figure 4(a) and material properties are shown in Table 1.
As shown in Figure 4(a), airflow is injected with a velocity of 300 m/s from the negative direction of  at the starting point, while, in Figure 4(b), air is injected from the positive direction of .When the pneumatic boundaries are applied at the initial stage, the two side walls applied to confine the assembly are removed simultaneously.Figure 5 shows the front view of simulation result from the initial stage to 1 sec.
It can be obtained from Figure 5 that gangue particles have different motion trajectories under different airflow velocity.As shown in Figures 5(a) and 5(b), gangue particles will do horizontal projectile motion before entering airflow domain, as can be seen at  = 0.2 s and  = 0.4 s.As shown in Figure 5(a), velocities of conveyor belt 1 V 0 and airflow  are in the opposite direction, and gangue particles will change their motion trajectories to an opposite direction compared with that of coal without the effect of airflow, which can be seen at  = 0.6 s,  = 0.8 s, and  = 1.0 s.As shown in Figure 5(b), velocities of the conveyor belt 1 V 0 and airflow  are in the same direction, and gangue particles will be blown much farther than that of coal without being affected by airflow, as can be seen at  = 0.6 s,  = 0.8 s, and  = 1.0 s.
As can be seen in Figures 5(a     In order to study the pneumatic separation distance influenced by different airflow directions, airflow velocities ±300 m/s are chosen for the research of separation effect.Conveyor velocity V 0 and height difference ℎ  between conveyor belt and air nozzle are kept constant.The relationships between separation distances and particle diameters under different airflow velocity directions are shown in Figure 6.As can be seen in Figure 6, separation distance decreases with the increase of particle diameter.It also can be concluded that the separation effect under airflow velocity −300 m/s is better than that under airflow velocity 300 m/s.The separation effect also can be analyzed from ( 9); dynamic pressure difference Δ  can be expressed as ( 2 − ẋ 2 )/2 when  and V 0 are in the same direction, while, when  and V 0 are in the opposite direction, the dynamic pressure difference Δ  can be expressed as ( 2 + ẋ 2 )/2.Thus, it can be obtained from the above analysis that separation effect under airflow velocity −300 m/s is significantly better than that of under airflow velocity 300 m/s.

Orthogonal Experiment of Pneumatic Separation.
Digital image processing technology has been used to identify the target of various patterns of coal and gangue in underground pneumatic separation system.Before mineral materials are sent to machine vision system, coal and gangue have been crushed to 100 mm by impact crusher.The size of materials is ranging from 50 mm to 100 mm.These materials are sent by the conveyor belt for coal and gangue digital image information identification (as shown in Figure 1).In this study, separation distance Δ (m) is selected as the primary index to evaluate the pneumatic separation effect.Pneumatic separation influenced by different airflow direction can be achieved by changing the arrangement of high-pressure value.The pneumatic separation test system mainly consists of conveyor belt, queuing system, machine vision system, control system, and high-pressure air injection system.When coal and gangue materials fell down from the conveyor belt, gangue will be identified by machine vision system, and the information is transported to control system through image sensor to drive the electromagnetic value open.Thus, coal gangue pneumatic separation is realized.The air compressor used in pneumatic separation testing system is LG-6.5/10, its working pressure is 1.0 Mpa, and certified capacity is 6.5 m 3 /min.
From the analysis shown in Section 3, the conveyor velocity V 0 , height difference ℎ  , and gangue diameter  are selected as the three factors.Factors and levels are listed in Table 2.According to the identified level of factors, orthogonal table L 9 (3 4 ) is applied in the test.Orthogonal experiment arrangement and results are shown in Table 3.

Variance Analysis and Range Analysis Based on Orthogonal Test.
As can be seen from Table 2, each factor at different levels is approximate linearity, so the method of regression analysis is to obtain the linear function relationship between them appropriately [25].Thus, the regression equations of coal gangue separation distance under different airflow direction are derived.Equation ( 25) represents the regression equation when  and V 0 are in the opposite direction, and (26) represents the regression equation when  and V 0 are in the same direction, where  is the separation distance of coal and gangue,  1 represents the diameter of gangue,  2 is the height difference, and  3 is the velocity of conveyor belt: Variance analysis is carried out on the regression equation to make significance test, and the results are shown in Tables 4 and 5.In order to determine the optimal pneumatic separation solutions, comparative analyses of the ranges between various levels of each factor are carried out.The range analysis is shown in Table 6.
In Table 6, subscripts 1 and 2 used in influence factors (, , ) represent the fact that  and V 0 are in the opposite direction and the same direction, respectively.  ( = 1, 2, . . ., ) is the sum of index values corresponding to factors in column  at level .The value of   determines the optimal level and combination of factors in column .  reflects the ranges of the index with the variation of factors in column , and the influence of the factor will be more significant if the value   is greater.
According to Tables 4 and 5, all the three factors mentioned above have significant influence on pneumatic separation distance.As shown in Table 6, the pneumatic separation distances Δ are all decreasing with the increased value of the three factors.The analysis shows that gangue diameters have the most significant influence on separation distance.

Analysis of Experiment Results Based on Support Vector Machine (SVM).
From the above analysis shown in Section 4.1, the significant degrees of different factors for separation distance are obtained by variance analysis.Besides, the primary and secondary relations of the influence factors with pneumatic separation distance could be obtained according to the range analysis.However, the optimal combination is a relative definition for the limited levels and has great one-sided characteristic.Most cases occurred in the experiment; the so-called "optimal combination of the factors" is a relative optimal, not the real optimal.
For further analysis of the experiment result, SVM [26][27][28][29][30] is introduced.The detailed functional forms of SVM are given in the Appendix.Optimization settings for factors that have influences on coal gangue pneumatic separation could be divided into the following steps: (1) collect the necessary data using orthogonal experiment; (2) set SVM learning model parameters and determine SVM kernel function; (3) input learning samples, and obtain parameters; (4) establish the fitting model according to parameters obtained above; (5) determine the levels of parameters in a certain range, and then combine these levels to establish a large number of input vector samples; (6) input vector sample into the fitting model and then obtain the output sample.
The relationship of separation distance Δ and various influence factors is obtained through the above analysis.As shown in Figure 7, and -axes represent two of the three influence parameters, respectively.14 values equally spaced from the range of parameters of orthogonal experiment are taken, respectively.Thus, the comprehensive collection of the two parameters could form 196 samples, and the extreme value is selected in the third parameter in orthogonal experiment.
Figures 7(a) and 7(b) show the relationships of separation distance Δ and the three factors when  and V 0 are in the opposite direction.As can be seen from Figures 7(a) and 7(b), pneumatic separation distance Δ decreases with the increase of conveyor velocity V 0 , the height difference ℎ  , and gangue diameter .It also can be concluded that gangue diameter  has the greatest influence on separation distance Δ.
Figures 7(c) and 7(d) have shown the relationships of separation distance Δ and factors when  and V 0 are in the same direction.As can be seen from Figures 7(c) and 7(d), there is the same variation trend as that of Figures 7(a) and 7(b).Through the analysis of the two groups of figures, pneumatic separation effect will be better when  and V 0 are in the opposite direction.It can be concluded that the SVM intelligent model has important guiding significance and practical value for coal gangue pneumatic separation.

Correction of the Theoretical Formula.
Based on the least square method, (23) and (24) in Section 3.2 can be transformed into a function of   and   parameters through variable substitution; then set up equations based on the experiments data.The result of   and   can be calculated finally.
When V 0 and  are in the same direction, nonlinear correction term can be given as   = 0.0205, and linear correction term   = 1.1.Thus, the formula can be expressed as When V 0 and  are in the opposite directions, nonlinear correction term   and linear correction term   can be given as   = 1.56 × 10 5 and   = 4.4.Thus, the formula is obtained as In order to verify the effectiveness of theoretical formula, nine samples shown in Table 3 are selected as verification samples.The comparison of experiment results and calculated results is shown in Table 7.
As can be seen from Table 7, calculated value and experiment value of coal gangue pneumatic separation distance have high consistent degrees, the separation distance error between experiment value and calculated result is less for gangue with smaller diameter, and the separation distance error increases with the increase of gangue diameter.However, the separation distance error between experiment value and calculated result is still less than 16%.Thus, it can be concluded that ( 27) and ( 28) have important guiding significance and practical value for coal gangue pneumatic separation.

Conclusions
The established coal and gangue pneumatic separation model reflects the basis motion law of gangue affected by airflow and coal without being affected by airflow, which provides two feasible solutions for underground pneumatic separation  under two kinds of arrangement scheme of high-pressure value.The analysis in theory is consistent with that obtained in the experiments, which validate the established theoretical model and present the following conclusions: (1) Different high-pressure value arrangement schemes have great influence on gangue pneumatic separation, the theoretical formulas of coal gangue pneumatic separation distance affected by different airflow direction are derived, and the expressions of the two formulas are different under different airflow direction.
(2) A series of air-solid multiphase flow simulations and orthogonal experiments were conducted to clarify its effect under different airflow direction.Based on the analysis, pneumatic separation effect will be better under the condition of  and V 0 being in the opposite direction.Pneumatic separation distance Δ decreases with the increased values of the three factors (conveyor velocity V 0 , height difference ℎ  , and gangue diameter ).These analyses also show that gangue diameters have the most significant influence on separation distance, followed by conveyor velocity V 0 and height difference ℎ  .
(3) The relationship of pneumatic separation distance Δ and influence factors was obtained by SVM intelligent model; the theoretical formulas of coal gangue pneumatic separation distance are corrected based on the analysis of orthogonal experiment data.The corrected formula is suitable to serve as the theory basis of coal gangue pneumatic separation.where x  is the input vector, y  is a real constant as the output vector,  is a normal vector for fitting function,  is threshold value, and  is the fitting precision.
Based on the principle of minimum structural risk, the optimization objective could achieve better generalization ability at the minimum value of ‖‖ 2 /2.Considering the existence of approximation error  ( is a real constant) in actual application, therefore, SVM can be expressed as The optimal solution of (A.1) could be derived by calculating the saddle points of the Lagrange equation.Thus, function approximation problem can be obtained: where  SV in (A.4) is the SVM and the training sample is the support vector when (  −  *  ) is not equal to zero.

Figure 2 :
Figure 2: Pneumatic separation model of coal and gangue.

22 × 5 × 5 :
, , and  directions) fluid cells are created in a rectangular space ( = [0, 7000 mm],  = [−400 mm, 400 mm], and  = [300 mm, 400 mm]), which covers the rectangular space.A pneumatic boundary should be set for the fluid grid.Driving forces from the fluid flow are applied to the particles as body forces.These forces are also added to the fluid equations and cause change in momentum, as reflected by the change in the pressure gradient in the flow direction.As shown in Figures 4(a) and 4(b), two available models are established to reduce the computation time without loss in accuracy.The gangue hopper, containing 300 balls, is built at the top right of conveying belt 1; it is aimed at reducing the computation time.The front view of the simulation models at an initial stage under different airflow directions is shown in Figures 4(a) and 4(b), respectively.
) and 5(b), gangue particles will change their trajectories in area A and area C. Areas B and D, as shown in Figures 5(a) and 5(b), illustrated that gangue particles with smaller diameters can be blown much farther.
V 0 and  in the opposite direction V 0 and  in the same direction
Influence of conveyor belt velocity V 0 and height difference ℎ  on gangue separation distance ( and V 0 in opposite direction) Influence of conveyor belt velocity V 0 and gangue diameter  on gangue separation distance ( and V 0 in opposite direction) Influence of conveyor belt velocity V 0 and height difference ℎ  on gangue separation distance ( and V 0 in the same direction) Influence of conveyor belt velocity V 0 and gangue diameter  on gangue separation distance ( and V 0 in the same direction)

Figure 7 :
Figure 7: Influence of three factors on gangue separation distance.

Table 2 :
Levels of factors.

Table 3 :
Experimental results of separation effect.

Table 4 :
Variance analysis ( and V 0 are in opposite direction).

Table 5 :
Variance analysis ( and V 0 are in the same direction).

Table 7 :
Experimental and calculated results of separation effect.
−  −  ≤  +     +  −   ≤  +  *    ≥ 0,  *  ≥ 0, (A.2) where  is balance factor, which is used to control the degree of punishment beyond the error sample, and   and  *  are relaxation factors.  and  *  are of the same nature; in general, relaxation factor at the top of fitting curve is recorded as   ; conversely, it is recorded as  *  .Equation (28) could change into quadratic programming problem based on dual theory.Then the Lagrange equation is established:  (, ,   ,  *  ,   ,  *  ,   ,  *  )   and  *  are Lagrange multiplier,   ≫ 0,  *  ≫ 0, and   and  *  are temporary variables,   ≫ 0,  *  ≫ 0.   and  *  have the same physical significance with   and  *  .