نوع مقاله : پژوهشی
نویسندگان
1 استادیارگروه مهندسی عمران، واحد اردبیل، دانشگاه آزاد اسلامی، اردبیل، ایران.
2 دانشیار گروه جغرافیا، واحد اردبیل، دانشگاه آزاد اسلامی، اردبیل، ایران
چکیده
مدلسازی صحیح فرآیند بارش-رواناب به دلیل گستردگی عوامل موثر بر بارش و رواناب یکی از پیچیدگیهای علم هیدرولوژی است. هدف از این پژوهش استفاده از روشهای پیش پردازش زمانی از جمله رفـع نویز موجکی و تبدیل موجک برای پیش بینی سریهای زمانی ماهانه رواناب در دشت اردبیل میباشد. شـبیه سـازی بـارش – روانـاب بـا اسـتفاده از مدل جعبـه سـیاه شبکه عصبی مصنوعی برای سه ترکیب دادههای بارش و رواناب دشت اردبیل انجـام گردید. ترکیب اول و دوم دادهها از دادههای خود ایستگاه در زمانهای گذشته استفاده میکند و ترکیب سوم دادهها از دادههای ایستگاههای بالادست (ایستگاههای گیلانده و کوزه تپراقی) برای پیش بینی رواناب خروجی دشت (ایستگاه سامیان) استفاده میکند. نتایج نشان داد کـه اعمـال روشهای پیش پردازش زمانی رفـع نویز مـوجکی و استفاده از تبدیل موجک در شــبیه ســازی بارش-رواناب بــا مدل شبکه عصبی مصنوعی به ترتیب بطور متوسط باعــث بهبــود 4 و 39 درصدی در مرحلــه آزمایش مدل شده است.
کلیدواژهها
عنوان مقاله [English]
Predicting runoff with pre-processing approaches in Ardabil plain
نویسندگان [English]
- Farnaz Daneshvar Vousoughi 1
- Rasoul Samadzadeh 2
1 Department of Civil Engineering, Ardabil Branch, Islamic Azad University, Ardabil, Iran.
2 Associate Professor Department of Geography, Ardabil Branch, Islamic Azad University, Ardabil, Iran
چکیده [English]
1-Introduction
Nowadays water resources management is a vitally important task and is the optimum planning of irrigation projects, and the development and exploitation of water resources especially during drought and flood events are strictly dependent on the accuracy of the used rainfall-runoff modeling tool. Therefore, different models have been already developed and employed for modeling rainfall-runoff processes of the watersheds (Partovian et al., 2017).
The wavelet-based pre-processing approach in the present study was used in the modeling of runoff time series via ANN. Furthermore, the impacts of denoising (smoothing) and wavelet transform have been simultaneously investigated in the accuracy of runoff prediction for one month ahead at the outlet of Ardabil plain.
2-Methodology
2-1-Case of the Study
The plain of Ardabil (38 – 38 N and 47 – 48 E), locatedin north-western Iran, covers an area of about 990 km2 (see Fig. 1). In the present study, the trend analysis was carried out on the rainfall (P) and runoff (R) parameters for three stations including Samian (PS, RS), Gilandeh (PG, RG), and Kozatopraghi (PK, RK)) located in the Ardabil plain from 1977 to 2019. The data sampling has been reported in the one-month interval at all of the stations. Figure 2 shows the locations of the rainfall and runoff stations. In this study, five combinations of input data were consumed for runoff prediction as to the following:
Comb. 1: RS(t), RS(t-1), PS(t); Comb. 2: RS(t), RS(t-12), RS(t-24), PS(t); Comb. 3: RK(t), RG(t), RG(t-1), PG(t-12), RK(t-12)
Fig. (1): Case of the study and the position of rainfall and runoff stations.
2-2-Artificial Neural Network (ANN)
Three-layered feed-forward backpropagation, which is usually used in forecasting hydrologic time series, provides a general framework for representing the nonlinear functional mapping between a set of input and output variables.
2-3-Wavelet transform (WT)
In hydrological problems, the time series are usually in the discrete but continuous format; therefore, the discrete WT was used in the following form (Mallat, 1998):
(2)
2-4-Wavelet based de-noising
Wavelet de-noising technique is operated as follows: (1) an applicable mother wavelet and several resolution level methods are selected. An approximation subseries at the resolution level L and detailed sub-series at different resolution levels are decomposed from main time series xi (2) The absolute amounts of detailed-sub-series, which exceed the values of the fixed threshold are changed by the difference between the values of threshold and detailed sub-series.
2-5- Efficiency criteria in runoff prediction
Two different criteria were used to measure the efficiency of the proposed forecasting methods; the root means square error(RMSE) and the determination coefficient (DC).
3-Results and Discussion
Some temporal features may also exist in the runoff time series due to their highly non-stationary fluctuations. To handle such features, wavelet-based temporal pre-processed data were entered into the ANNs to improve the accuracy of runoff modeling. WT and wavelet-based de-noising approaches were used for modeling the rainfall-runoff process via the ANN model. The Daubechies-4(db4) mother wavelet, which is almost similar to the runoff signal could capture the features of the signal, especially peak values, thus, it was selected as the mother wavelet for the decomposition of the runoff time series in this study. The decomposition of runoff time series at level L yields L+1 sub-signals (one approximation sub-signal, Pa(t) and L detailed sub-signals, Pdi(t) (i=1, 2, …, L)). Decomposition level 3 was considered as the optimum decomposition level. Each of the decomposed sub-series of the runoff demonstrated a specific seasonal feature of the process. In WT-ANN (WANN) model, decomposed sub-series accompanied by the rainfall and runoff data of each compound were used in the FFNN to predict one-month-ahead runoff values at the outlet of Ardabil plain (Samian station). In the second stage, the runoff time series were denoised via WT, and the denoised runoff data were used to predict the runoff at Samian station for one month ahead. Finally, the ANN model was compared with ANN models using pre-processing inputs.
The results of three models for one-step-ahead runoff forecasting at Samian station have been presented in Table 1. Results indicated that better accuracy was comprised with another model via the WANN model in the comb. 3. WANN models via comb. 3 used the runoff data of Gilandeh and Kozatopragi that lied in the upstream and showed accurate performance. These demonstrated Gilandeh and Kozatopragi runoff time series played an important role in Samian runoff modeling. Accuracy improvement in the WANN model was 17%, 3.5%, and 35% combs. 1, 2, and 3 of inputs. The ANN model with denoised inputs showed little improvement (1, 6, and 6.2 percent in combs. 1, 2, and 3 of data) in runoff modeling at the outlet of the plain.
Table (1): The results of ANN and SVM models for one-step-ahead predictions
Input combination
Output
variable
Model Type
DC
RMSE (Normalized)
Calibration
Verification
Calibration
Verification
1
RS(t+1)
ANN
ANN with denoised data
WANN
0.592
0.594
0.791
0.434
0.438
0.587
0.065
0.065
0.047
0.051
0.052
0.044
2
RS(t+1)
ANN
ANN with denoised data
WANN
0.795
0.813
0.791
0.567
0.601
0.587
0.043
0.039
0.047
0.0434
0.042
0.044
3
RS(t+1)
ANN
ANN with denoised data
WANN
0.907
0.831
0.880
0.730
0.775
0.854
0.029
0.040
0.032
0.035
0.032
0.026
4-Conclusions
In this study, the wavelet-based denoised data and WT were employed in ANN for rainfall-runoff modeling at the outlet of Ardabil plain using data pre-processing techniques. Accordingly, first, it was sought to smooth the hydrological time series by eliminating the outliers and large noises of raw observed time series, which may be due to human or tool measurement error or systematic error. Then, different sub-series were generated by decomposing runoff time series and used to train the ANN model for rainfall-runoff modeling. Using processed and unprocessed data, the obtained results were compared; this comparison indicated the merit of applied data pre-processing approaches due to robust identification of hidden patterns in data so that the developed models could simulate and predict runoff values with lower margin of error and higher confidence and the best results were achieved by employing the decomposed runoff data via WT having different training time series with the same components of original time series. For future study, it is recommended to examine the efficiency of the proposed data pre-processing method in the rainfall-runoff modeling of other watersheds since it is expected that the merit of the method is more highlighted where the quality of the collected data is blurred due to the technical limitations. Furthermore, it is suggested to evaluate the efficiency of the proposed method in modeling the process at other time scales and for modeling other hydrological processes which may involve distinct noise levels and patterns regarding the type of process.
Keywords: Runoff modeling, Wavelet Transform, Wavelet-based de-noising, Artificial Neural Network (ANN), Ardabil plain
5-References
Donoho, D.H., 1995. Denoising by soft-thresholding. IEEE Transactions on Information Theory. 41(3):613–617.
Mallat, S.G., 1998. A Wavelet Tour of Signal Processing, second ed. Academic Press, San Diego. Partovian, A., Nourani, V., Aalami, M.T., 2016. Optimizing Neural Network for Monthly Rainfall-Runoff Modeling with Denoised-Jittered Data. Journal of Tethys. 4(4), 284–294.
کلیدواژهها [English]
- Runoff modeling
- Wavelet Transform
- Wavelet based de-noising
- Artificial Neural Network (ANN)
- Ardabil plain