Abstract
In this paper, the sparrow search algorithm is used to predict the discharge coefficient (Cd) of the triangular side orifice for the first time. Dimensionless parameters influencing the size of Cd of side orifices are obtained as input values and discharge coefficient as output values of the model. The results show that the determination coefficient R2 is 0.973, the root means square error RMSE is 0.0122, and the average absolute percentage error is 0.010% in the testing phase. The model has high forecast accuracy, strong generalization ability and higher accuracy than other models and traditional empirical formulas. Quantitative analysis by Sobol's method shows that the ratio W/H of top orifice height to side orifice height, Fr of upstream Froude number, and ratio B/L of channel width to a bottom edge length of side orifice are the main factors influencing the discharge capacity of triangular side orifice. The first-order sensitivity coefficient and global sensitivity coefficient are 0.23, 0.11, 0.17 and 0.41, 0.39, 0.35 respectively.
HIGHLIGHTS
In this paper, the sparrow search algorithm is used to predict the discharge coefficient of the triangular side orifice for the first time.
The importance of dimensionless parameters on discharge coefficient is quantified using the Sobol's method.
The flow characteristics of the side orifice are analyzed.
Graphical Abstract
INTRODUCTION
With the increasing shortage of water resources, it is imperative to implement water-saving irrigation. It is particularly important to study measuring equipment with high precision, strong adaptability and good convenience (Wang et al. 2016). As a hydraulic control structure widely used in hydraulic and irrigation works, side orifices in the open channel are used to divert water discharge from the main channel to other channels (Hussain et al. 2011a). Because of the importance of this diversion structure, many experiments, analyses and numerical studies have been carried out on its hydraulic characteristics. The discharge coefficient is the most important hydraulic parameter of the side orifice. Therefore, it is of great significance to calculate the discharge coefficient of the side orifice accurately in the irrigation system and hydraulic engineering.
In recent decades, to measure water accurately, different scholars have carried out in-depth research on hydraulic characteristics of triangular side orifices through model tests and given the calculation formula of discharge coefficient as shown in Table 1. Vatankhah & Mirnia (2018) analyzed hydraulic characteristics of triangular side orifices through model tests and obtained a discharge coefficient relationship based on the classical regression method. Jamei et al. (2021) used multiple linear regression to give the discharge coefficient calculation formulas for both known and unknown upstream channel discharges. Vatankhah (2019) developed the calculation equation of discharge coefficient of side orifice considering the change of orifice centerline. Considering that the discharge coefficient equation fitted by different scholars is difficult to be unified, it is not convenient for people to use. Therefore, it is necessary to develop a systematic and highly accurate discharge coefficient calculation method to solve the problem of discharge measurement in small channels.
Calculation formula of discharge coefficient fitted by different scholars
Serial number . | Equations . | Authors . |
---|---|---|
1 | ![]() | Vatankhah & Mirnia (2018) |
2 | ![]() | Jamei et al. (2021) |
3 | ![]() | Jamei et al. (2021) |
4 | ![]() | Vatankhah (2019) |
Serial number . | Equations . | Authors . |
---|---|---|
1 | ![]() | Vatankhah & Mirnia (2018) |
2 | ![]() | Jamei et al. (2021) |
3 | ![]() | Jamei et al. (2021) |
4 | ![]() | Vatankhah (2019) |
Where Cd is the discharge coefficient of the side orifice; L is the orifice length (m); H is the orifice height (m); B is the main channel width (m); W is the orifice crest height (m); y1 is the flow depth upstream of the channel (m); Fr1 is the upstream Froude number.
In recent years, with the development of artificial intelligence technology, many scholars have successfully applied this technology to solve complex hydraulic engineering problems. Azimi et al. (2017) used the adaptive neuro-fuzzy inference system (ANFIS) and a hybrid of ANFIS and a genetic algorithm (ANFIS-GA) to predict the discharge coefficient of rectangular side orifices. At the same time, the discharge rate of side orifices was simulated by FLOW-3D software. The results show that ANFIS-GA has the highest accuracy. Moghadam et al. (2019) used ANFIS and Firefly algorithm (FA) to predict the discharge coefficient of side orifices. The results show that ANFIS-FA has stable performance. Eghbalzadeh et al. (2016) used an artificial neural network to predict the discharge coefficient of rectangular side orifices. The results show that the predicted results of all neural networks are better than those of non-linear regression, and the results of radial basis neural network (RBNN) are better than those of generalized neural network (GRNN) and feed-forward neural network (FFNN). Jamei et al. (2021) and Qian et al. (2019) carried out prediction research on the discharge coefficient of triangular side orifices through the intelligent model. To explore the sensitivity of dimensionless parameters of the model, they arranged and combined the dimensionless parameters affecting Cd into 26 different input combinations. This method not only increases the calculation amount, but also tends to ignore the interaction of different parameters on the discharge coefficient. Based on the current literature research, most scholars focus on comparing the accuracy of the rectangular side-orifice intelligent model but have little research on triangular side-orifice. In addition, the current research lacks quantitative analysis of model input parameters and ignores the specific contribution of dimensionless parameters to Cd.
Therefore, this study aimed at the prediction of the side orifice discharge coefficient. Through Buckingham-π theorem, the dimensionless parameters that affect the triangular side orifice discharge coefficient were obtained as the input parameters of each model, and the discharge coefficient was taken as the output parameter of each model. Considering that back propagation neural network (BPNN) was easy to fall into local optimum, a sparrow search algorithm was used to optimize its weights and thresholds, and three statistical indexes were used to get the performance and accuracy of the model, and the results were compared with those of the literature. Finally, the input parameters of the model were quantitatively analyzed by Sobol's method to determine the sensitivity of each dimensionless parameter of the Cd and the variation rule after the interaction. This study explored the application of artificial intelligence technology in hydraulic engineering and provided new ideas for the design and application of side orifices for the optimization of water use and farmland water-saving irrigation.
DATA AND MODELS
Experimental datas
The range of parameters in the experimental data set
Serial number . | Height H/cm . | Length L/cm . | Bottom height W/cm . | Channel discharge Qu/(L/s) . | Side-orifice discharge Qs/(L/s) . | Upstream water depth y1/cm . | Water depth in the center of the orifice yc/cm . | Downstream water depth y2/cm . |
---|---|---|---|---|---|---|---|---|
Min | 4 | 30 | 5 | 13.33 | 1.77 | 9.41 | 10.48 | 10.82 |
Max | 10 | 40 | 10 | 34.64 | 17.58 | 28.57 | 28.86 | 82.80 |
Serial number . | Height H/cm . | Length L/cm . | Bottom height W/cm . | Channel discharge Qu/(L/s) . | Side-orifice discharge Qs/(L/s) . | Upstream water depth y1/cm . | Water depth in the center of the orifice yc/cm . | Downstream water depth y2/cm . |
---|---|---|---|---|---|---|---|---|
Min | 4 | 30 | 5 | 13.33 | 1.77 | 9.41 | 10.48 | 10.82 |
Max | 10 | 40 | 10 | 34.64 | 17.58 | 28.57 | 28.86 | 82.80 |
Dimensional analysis
Data preprocessing
Model building
BPNN
Back propagation neural network, as a traditional multi-layer FFNN, consists of an input layer, a hidden layer and an output layer. It mainly includes forward propagation of signal and reverses transmission of error. For error back-propagation, the output errors of each neuron layer are calculated step by step through the output layer, and then the weights and thresholds of each layer are adjusted according to the error gradient descent method so the final output of the modified network can approach the expected value. Due to the small amount of calculation and strong parallelism, it has been widely used in the engineering field (Shi et al. 2019; An et al. 2021; Yang et al. 2021).
SSA-BPNN

ELM







Therefore, , where
is Moore-Penrose generalized inverse matrix of matrix H. Studies have shown that ELM is a general-purpose approximator, and its simplification contributes to rapid calculation and high generalization ability, and is suitable for various situations (Sun et al. 2008; Luo & Zhang 2014; Roushangar et al. 2018).
Sobol's sensitivity analysis
Evaluation index


RESULTS AND DISCUSSION
Model comparison
Scatter fit plots of two models in the testing and training phases. (a) Training stage. (b) Testing stage.
Scatter fit plots of two models in the testing and training phases. (a) Training stage. (b) Testing stage.
Table 3 shows the performance indicators of the two models at different stages. The values of RMSE and MAPE are closer to 0, the smaller the model error and average deviation. Table 3 shows the RSME of SSA-BPNN decreases by 45.16 and 37.11% respectively, based on BPNN in the training and testing stages, indicating that the measured value of the Cd of the SSA-BPNN deviates less from the predicted value. Compared with BPNN, MAPE of SSA-BPNN decreases by 42.86 and 37.5% respectively, which indicates that SSA-BPNN has a better prediction effect and higher accuracy. Therefore, after using the SSA algorithm, the performance of BPNN is improved effectively, and it is enhanced for the global search ability of the neural network.
Performance indicators of BP and SSA-BP models at different stages
Model . | Training stage . | Testing stage . | ||||
---|---|---|---|---|---|---|
RMSE . | MAPE (%) . | R2 . | RMSE . | MAPE . | R2 . | |
BPNN | 0.0186 | 0.007 | 0.873 | 0.0194 | 0.016 | 0.881 |
SSA-BPNN | 0.0102 | 0.004 | 0.962 | 0.0122 | 0.010 | 0.973 |
Model . | Training stage . | Testing stage . | ||||
---|---|---|---|---|---|---|
RMSE . | MAPE (%) . | R2 . | RMSE . | MAPE . | R2 . | |
BPNN | 0.0186 | 0.007 | 0.873 | 0.0194 | 0.016 | 0.881 |
SSA-BPNN | 0.0102 | 0.004 | 0.962 | 0.0122 | 0.010 | 0.973 |
Table 4 shows the performance index of ELM model in the training and testing stages. When the number of neurons in the hidden layer is 20, the model achieves the best performance (shown in bold in Table 4). In the training stage, RMSE = 0.011, MAPE = 0.004, R2 = 0.955. In the testing stage, RMSE = 0.012, MAPE = 0.010, R2 = 0.953. Obviously, the predictive and fitting ability of SSA-BPNN is better than ELM. Therefore, SSA-BP is the best intelligent model to predict the discharge capacity of triangular side orifices.
Performance indicators of ELM models at different stages
Number of neurons . | Training stage . | Testing stage . | ||||
---|---|---|---|---|---|---|
RMSE . | MAPE (%) . | R2 . | RMSE . | MAPE (%) . | R2 . | |
5 | 0.019 | 0.007 | 0.869 | 0.018 | 0.016 | 0.894 |
10 | 0.016 | 0.006 | 0.911 | 0.016 | 0.014 | 0.910 |
15 | 0.012 | 0.004 | 0.945 | 0.014 | 0.012 | 0.940 |
20 | 0.011 | 0.004 | 0.955 | 0.012 | 0.010 | 0.953 |
25 | 0.011 | 0.004 | 0.950 | 0.012 | 0.010 | 0.945 |
Number of neurons . | Training stage . | Testing stage . | ||||
---|---|---|---|---|---|---|
RMSE . | MAPE (%) . | R2 . | RMSE . | MAPE (%) . | R2 . | |
5 | 0.019 | 0.007 | 0.869 | 0.018 | 0.016 | 0.894 |
10 | 0.016 | 0.006 | 0.911 | 0.016 | 0.014 | 0.910 |
15 | 0.012 | 0.004 | 0.945 | 0.014 | 0.012 | 0.940 |
20 | 0.011 | 0.004 | 0.955 | 0.012 | 0.010 | 0.953 |
25 | 0.011 | 0.004 | 0.950 | 0.012 | 0.010 | 0.945 |
Comparison with literature research
Comparison with previous models
Model . | RMSE . | MAPE(%) . | R2 . |
---|---|---|---|
SSA-BP | 0.0122 | 0.010 | 0.9730 |
GA-SVR (Hussain et al. 2010) | 0.0135 | / | 0.9205 |
Empirical (Vatankhah & Mirnia 2018) | 0.0150 | 2.370 | 0.9216 |
Model . | RMSE . | MAPE(%) . | R2 . |
---|---|---|---|
SSA-BP | 0.0122 | 0.010 | 0.9730 |
GA-SVR (Hussain et al. 2010) | 0.0135 | / | 0.9205 |
Empirical (Vatankhah & Mirnia 2018) | 0.0150 | 2.370 | 0.9216 |
Sensitivity analysis
First-order sensitivity coefficient S1 and global versus sensitivity coefficient Si of different parameters.
First-order sensitivity coefficient S1 and global versus sensitivity coefficient Si of different parameters.
Variation of the most influential inputs predicted discharge coefficient.
CONCLUSION
In this study, artificial intelligence technology was applied to solve the complex hydraulic characteristics of triangular side orifices. Firstly, dimensionless parameters influencing the Cd of side orifices were obtained by dimension analysis of the Buckingham-π theorem as the input value and Cd as the output value of the model. The predicted value of the Cd was obtained through the algorithm learning principle. To further obtain the variation rule of influence of dimensionless parameters on Cd, Sobol's method was used to quantify the parameters, and the following conclusions were obtained:
- (1)
SSA-BPNN can be used to predict the Cd of triangular side orifices, where maximum error and average error are 6.56 and 1.73%, respectively. The model has high accuracy, strong generalization ability and higher accuracy than other intelligent models and empirical equations.
- (2)
The first-order sensitivity of W/H, B/L, B/H, y1/H and Fr is 0.23, 0.11, 0.03, 0.05 and 0.17 respectively, and the global sensitivity is 0.41, 0.39, 0.21, 0.31 and 0.35 respectively. To optimize the discharge capacity of the side orifices, W/H, Fr and B/L should be considered as the primary factors in the design of similar projects.
- (3)
When 0.63 < B/L < 0.68 and 0.50 < W/H < 0.65, the discharge capacity of the side orifice is the largest; when 0.80 < B/L < 0.83 and 2.40 < W/H < 2.50, the discharge capacity of side orifices is the smallest.
- (4)
Considering the complexity of hydraulic characteristics of side orifices, an intelligent method not only provides a new method to solve similar engineering problems but also further improves the accuracy and efficiency of the Cd of triangular side orifices, which provides an important reference basis for accurate water consumption for farmland irrigation.
Finally, the flow regime of the upstream of the side orifice is measured by the size of Fr, and the flow characteristics of the side orifice are further analyzed in the future research.
DATA AVAILABILITY STATEMENT
All relevant data are available from an online repository or repositories (https://doi.org/10.1061/(ASCE)IR.1943-4774.0001343).
CONFLICT OF INTEREST
The authors declare there is no conflict.