Simulation of drainage process of subsurface pipe based on improved HYDRUS-2D model
View Fulltext  View/Add Comment  Download reader
  
DOI:10.7606/j.issn.1000-7601.2020.04.28
Key Words: boundary of subsurface pipe  soil salinization  soil salt content  drainage  salt discharge  soil water content  HYDRUS-2D model
Author NameAffiliation
LIU Yi State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
ZENG Wenzhi State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
AO Chang State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
LEI Guoqing State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
WU Jingwei State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
HUANG Jiesheng State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan, Hubei 430072, China 
Hits: 697
Download times: 620
Abstract:
      In order to improve the simulation accuracy of HYDRUS-2D in subsurface pipe drainage, this paper proposed two strategies of virtual soil layer (VSL) and actual opening area (AHA) to replace the previous seepage boundary (PSB) of HYDRUS-2D model in subsurface pipe drainage. The results showed that the simulated values of the VSL and AHA models were in good agreement with the measured values: the square root error (RMSE) of the salinity was 0.757 g·kg-1 and 0.966 g·kg-1, respectively, and the determination coefficient (R2) was 0.977 and 0.964, respectively. The RMSE of soil water content was 0.007 cm3·cm-3 and 0.008 cm3·cm-3 and the R2 was 0.885 and 0.794, respectively. In addition, the simulation accuracy of VSL in terms of drainage process (RMSE=9.48 L, R2= 0.93) and salt discharge process (RMSE=0.225 kg, R2= 0.922) were higher than that of AHA and PSB. It is obvious that VSL was more suitable for subsurface pipe boundary. At the same time, this paper further fitted the empirical formula between the ring width (RW) of VSL and its saturated hydraulic conductivity (Ks) (R2=0.99, P<0.01). The empirical formula could be used to determine the model parameters of VSL when the layout of simulated domain and subsurface pipe changes, and the simulation results could be used to guide the layout of subsurface pipe in relevant areas.