The Study of Subsurface Land Drainage Optimal Design Model

This paper focused on choosing the best design of subsurface land drainage systems in semiarid areas. The study presented three diﬀerent soil layers with diﬀerent hydraulic conductivity and permeability, all layers are below the drain level, and the permeability is increasing with depth. A mathematical model was formulated for the horizontal and vertical drainage optimal design. The result was a nonlinear optimization problem with nonlinear constraints, which required numerical methods for its solution. The purpose of the mathematical model is to ﬁnd the best values of pipes and tubewells spacing, groundwater table drawdown, and pumps operating hours which leads to a minimum total cost of the subsurface drainage design. A computer code was developed in MATLAB environment and applied to the case study. Results show that the vertical drainage was economically better for the case study drainage network design. And the main factor aﬀecting the mathematical model for both pipe and well drainage was the distance between pipes and tubewells. In addition, considering the lifespan of vertical drainage project, the optimal design involves the minimum possible duration of pumping stations. It is hoped that the proposed optimal mathematical model will present a design methodology by which the costs of all alternative designs can be compared so that the least-cost design is selected.


Introduction
Subsurface drainage is widely used worldwide to remove excess water found below the earth's surface [1]. While surface drainage removes excess water from the soil surface before it enters the crop root zone, subsurface drainage decreases the groundwater level and provides a better environment for the crop growth [2]. Agricultural lands affected by salinity and high water levels generally require subsurface drainage [3]. Two methods of subsurface drainage are commonly used which are horizontal pipe drainage and vertical well drainage. Subsurface pipe drainage is a group of lateral pipes spread below the ground surface with specific depths, diameters, and spacing between them and connected with collector pipes that draw the surplus water out of the study area. Vertical drainage is a group of vertical wells, with specific diameters, depths, and spacing between them. ese wells are equipped with pumps that withdraw water, so the water level will decrease in the surrounding area. Land drainage helps to achieve water balance in the soil, prevent its salinization, and create favourable conditions for plant growth. Around the world, agricultural drainage plays a significant role in protecting investments in irrigation projects and agricultural production. It helps also in preserving soil resources and food production, by improving land productivity and crop yields, especially in poorly drained soils [4]. Most studies that involve subsurface drainage were about the design of the drainage network, and the task is to find the best design or the optimal design. e most important strategy in designing the subsurface drainage is choosing proper values for the spacing and depth of drains to minimize the total cost [5]. And a lot of researchers have studied the cost analysis of subsurface drainage systems for different types [6][7][8][9][10][11]. Cimorelli et al. [12] have worked on the surface drainage network optimal design and introduced a novel procedure for the best surface drainage network design. A methodology was presented by Chahar and Vadodaria [8], for determining the optimal spacing of ditches fully penetrating into an isotropic and homogeneous porous medium over an impervious layer. Another methodology was introduced by Moradi-Jalal et al. [13], for the optimal operation and design of pumping stations based on solving a nonlinear large-scale programming problem. And a new management model was presented by Moradi-Jalal et al. [14], for the optimal operation and design of water distribution systems. An objective function was constructed by Sharma and Swamee [15], for the cost structure of a pipe network system. e optimization model described by Bennett and Mays [16] was based upon dynamic programming for nonserial systems that determines the minimum cost of a drainage channel system and the detention for a watershed. A nonlinear square fitting routine was used by Wall and Miller [17] to determine the effective values of conductivity and moisture characteristics of the model of soil water drainage. And a simulation model of the optimal drainage network was presented by Howard [18], in which channels shift to minimize the total stream power within the network.
is paper presents a novel strategy for the best design of subsurface horizontal and vertical drainage in an area of different saturated soil layers with different hydraulic conductivity and permeability, and the permeability is increasing with depth. With the use of modern optimization algorithms, we can find the suitable values of groundwater table drawdown and pipes/tubewells spacing that lead to the minimum cost of the total subsurface drainage system.

Methodology
e least-cost design is that satisfying all design constraints with the minimum total cost. e objective function for pipe and well drainage was determined by considering all cost components that affect the drainage network design. en constraints were formulated depending on the hydraulic study of the study area. e result was a nonlinear objective function with nonlinear constraints. A survey of modern optimization algorisms was conducted to find the one suitable for the solution of the formulated optimization problem. It was found that the interior-point optimization algorithm was adapted to the problem and produced satisfactory results. Two computer codes for both horizontal and vertical subsurface drainage were developed in MAT-LAB environment in order to derive the optimal solution for both types of drainage systems. e solutions then were compered to know how the lifespan and type of the project will affect the total cost of the network design.

Pipe Drainage Costs
e major cost components of pipe drainage system are drainage materials, installation, and operation and maintenance.

Drainage Materials.
e costs of drainage materials and drainage pipe installation work and structures are major components of the total cost of a pipe drain project. Drainage materials include collector and lateral pipes, filters, pump set and pump house, outlet structures and manholes, and outlet pitching.
Collector and lateral pipes: the collector pipes have a transmission function to carry water to the outlet under the gravitational flow, so the UPVC (unplasticized poly vinyl chloride) corrugated nonperforated pipe is used for the collectors.
e pipe diameter of the collectors is chosen depending on the expected flow. Single Wall Perforated Flexible Corrugated UPVC Pipes (outer diameters of 80-355 mm) are widely used for lateral pipes.
Artificial filters/envelopes: the filter or envelope material around the pipes plays an important role in preventing the fine particles of the soil from entering into the pipes with the flow. Nowadays, it is preferred to use artificial filters rather than gravel filters. e artificial filters are cheaper in cost, of better quality, and easier to handle and transport than gravel filters. For medium and light soil, it is preferred to use nonwoven polypropylene fibers with a thickness of more than 2.5 mm and an opening size of more than 300 microns, and a needle punched geotextile nonwoven fabrics are used for heavy soil. Perforated pipes precoated with industrial filter are more preferred than locally coated pipes for quality installation.
Pump set and structures: pump set and structure consists of diesel pump set and pump house, manholes (junction box), and outlet structures. Other cross drain structures may also be required in larger projects but these are avoided for simplicity. e outlet structure and manholes are precast Reinforced Cement Concrete (RCC) pipes of different lengths and diameters depending on the site conditions. In the case of pump outlets, 900 mm diameter manholes and 1200 mm outlet structures are used, where 900 mm diameter manholes are used in the gravity outlets. e bottom edge of the RCC pipe for outlet structures and manholes is closed by fixed plate at the bottom. Plastic coated iron bars are provided in the walls of outlet structures and manholes to help in inspection and cleaning. In the case of gravity outlets, instead of RCC pipes, plastic manholes and outlet pipes have recently been used. Plastic manholes and outlet pipes are relatively expensive, but they are easy to carry, handle, and install at the project sites [19]. A small stone structure is constructed at the outlet to protect the pipe end from collapsing under the conditions of gravity outlet. A rodent guard is provided to prevent the pipes from getting stuck and damaged by rodents that may enter the lateral pipes.

3.2.
Installation. Installing the drainage system is the function of the installation unit according to the specifications and design. Supervisors must specify that the installation (including drainage equipment) is carried out strictly in accordance with specifications and design. Various drainage machines such as hydraulic excavators, tractor-mounted trenchers, self-propelled trenchers with laser automatic control, and self-propelled machines with L or V plough and laser control are used to install collector and lateral pipes [20,21]. Currently, laser control self-propelled trenchers are used in large-scale quality of subsurface drainage pipes installations. Other auxiliary machineries like bulldozers, excavators, tractors with trailers, backhoe, etc. are also used for the movement of manpower, drainage materials, and installation of outlets and manholes. Individual farmers can use tractor mounted trencher in small areas (1-5 ha) to install subsurface drainage pipes. In such cases, it may be desirable to use a laser control device with a tractor-mounted trencher to achieve suitable slope for collector and lateral pipes.

Operation and Maintenance.
e popular belief that subsurface drainage does not require any maintenance and operation is untenable. In the case of pump outlets, pump operation is required for at least the first few years of installation. e maintenance of subsurface drainage systems mainly involves removing sediment from outlets, manholes, and pipes, also repairing or replacing the damaged outlets, manholes, and pipes [20]. In controlled drainage systems, operations may also include closing and opening gates to reuse drained water to irrigate crops.

Pipe Drainage Cost Equation.
e most cost components that affect the subsurface pipe drainage design can be determined according to the total costs as shown in the following relationship: where i is the number of the choice (−), z i is the total costs for the choice i ($), E 0 is the interest rate (%), W i is the construction costs for the choice i ($), D is the lifespan of the drainage project (years), and U i is the annual investment costs for the choice i ($/year). e construction costs are calculated by using the following equation: where l i is the total length of all drainage pipes (m), A cut is the cross-sectional area of excavation (m 2 ), C cut is the unit cost of excavation ($/m 3 ), IC cut is the increase in unit cost of earthwork per unit depth of excavation ($/m 3 /m), h cut is the depth of excavated ditch (m), C p is the unit cost of drainage pipes and filters ($/m), C m is the unit cost of manholes per depth ($/m), and n m is the total number of manholes (−), which can be determined by where A tot is the total area of the study area (m 2 ), E is pipe spacing for subsurface pipe drainage network (m), and S is the distance between manholes along the drainage pipe (m). e cross-sectional area of excavation (A cut ) is a trapezoid section as shown in Figure 1 and it can be calculated by where b is the bottom width of excavation, h cut is the total depth of excavation, and ∅ is the lateral slop angle. Also, total length of all drainage pipes ( l i ) can be replaced by where l is the mean length of all lateral drainage pipes (m) and n i is the total number of drainage pipes (−), and it can be determined by using the following equation: where A i is the mean area served by drainage pipes (m 2 ), and it can also be determined by using the following equation: So, the construction costs can be written as Mathematical Problems in Engineering 3 Moreover, the annual investment costs are where P is the rate of depreciation, and it is taken as 15-16 years. U t i is annual cost of maintenance and service for the drainage network ($). And it can be calculated by using the following relationship: where C s is the annual cost of maintenance and service for each manhole ($). e objective is to design the least-cost pipe drainage network. us, the objective function can be stated as

Well Drainage Costs
It is recommended to establish a proper strategy for vertical drainage system design, and it should be connected with the economic factors. For example, we can choose large numbers of tubewells with a small amount of discharge and a slight decrease in the groundwater level from each tubewell, or we can choose small numbers of tubewells with a larger amount of discharge and a larger decrease in the groundwater level from each tubewell and a larger spacing between the tubewells. ere are a lot of choices, and these choices are controlled by (a) Tubewell depth (b) Tubewell spacing (c) Tubewell discharge (d) e decrease amount of the groundwater level And the most suitable choice can be determined according to the total costs as stated in (1). e construction costs for vertical drainage project are calculated by using the following equation: where n i is the total number of tubewells, l H is the total depth of the tubewell (m), and C 1 is the construction cost of each tubewell per one meter depth ($/m). Moreover, the annual investment costs are where U t i is the annual cost of maintenance and service for the tubewells ($). And it can be calculated by using the following relationship: where C 2 is the annual cost of maintenance and service for each tubewell ($) and U N i is annual cost of electricity for the withdrawal of water from each well ($). And it is given by where Q i is the discharge from each tubewell (m 3 /day), h i is the pumping depth (m), t is the pumps operating hours in the year (day), C 3 is the cost of kilowatt hours of electricity ($/kw.h), and η 0 is the pump efficiency (−). So, the objective function for the design of vertical drainage tubewells can be stated as

Case Study Description.
A total area of 500 hectares is attended to apply an irrigation network to meet the crops need and help in washing the salinity that comes from underground water. So, a drainage network is needed also along with the irrigation network to serve in releasing exceed salty water out of the study area. e study area is located in Syria, as shown in Figure 2. e net area of cultivation is 405 hectares. About 15% of the net area is covered by summer vegetables, and wheat occupies the largest proportion among the crops list 4 Mathematical Problems in Engineering which is about 44.3%, followed by cotton (35%) and barley (12.7%). Table 1 gives the general climate and crops indicators for the region. e study area contains three different soil layers (h 1 � 7 m, h 2 � 12 m, and h 3 � 15 m) with different hydraulic conductivity (k 1 � 0.8 m/day, k 2 � 3.2 m/day, and k 3 � 11.2 m/day), and there is a semi-impermeable layer under the third layer with thickness of h 4 � 2 m and hydraulic conductivity of k 4 � 0.01 m/day, as shown in Figure 3.
e groundwater table is about 3 m below the ground surface, the gradient of the groundwater surface at the beginning of the study area is J 1 � 0.006, and the gradient of the groundwater surface at the end of the study area is J 2 � 0.001. e general gradient of the study area is from north to the south, so the length of groundwater feeding line is L � 2430 m, leakage from irrigation network is φ � 6820.789 m 3 /ha/year, and the increasing amount of irrigation discharge to avoid salinity problems is g � 627.276 m 3 /ha/year. Unit cost of excavation is C cut � 50 (RMB/m 3 ), increase in unit cost of earthwork per unit depth of excavation is IC cut � 3 (RMB/m 3 /m), unit cost of drainage pipes and filters is C p � 10 (RMB/m), unit cost of manholes per depth is C m � 200 (RMB/m), annual cost of maintenance and service for each manhole is C s � 10 (RMB), ∅ � 0.78, and distance between manholes along the drainage pipe is S � 100 m. Construction cost of each tubewell per one meter depth is C 1 � 860 (RMB/m), annual cost of maintenance and service for each tubewell is C 2 � 80 (RMB/year), the cost of kilowatt hours of electricity is C 3 � 0.52 (RMB/kw·h), pump efficiency is η 0 � 0.78 , specific given factor is µ � 0.06, rate of depreciation is P � 1/15, and interest rate is E 0 � 8%. e drainage depth is m � 1 m.

Pipe Drainage Spacing Formulation.
We can determine the pipe spacing for subsurface pipe drainage network when we have three different layers and the permeability is increasing with depth by using the following formulations [22]: where As σ, σ 1 , and σ t can be calculated by using the following formulations: where H is the hydraulic head (m), which is the water table height above the drainage pipe at the midpoint between the drainage pipes as seen in Figure 4, T is the weighted transmissivity for all layers contributing to the flow (m 2 / day), q c is the drainage unit discharge that must be released by drainage pipes (m/day), K i is the hydraulic conductivity for each layer (m/day), h i is the thickness of each saturated where m is the drainage depth (m), which is the minimum depth required by each plant for better productivity. And we can obtain β 1 , β 2 , and β t values from Figure 5, after calculating the factors ψ, λ, and ε as seen in Table 2.
e drainage unit discharge (q c ) that must be released by drainage pipes can be obtained by studying the water balance in the study area.

Hydraulic Study of Drainage Pipes.
For the hydraulic study of subsurface drainage pipes, we can use Manning formula which is as follows: where Q is the discharge that must be drained by the drainage pipes (m 3 /sec), n is the roughness coefficient (−), R is the hydraulic radius (m), A is the water cross section area (m 2 ), and I is the hydraulic gradient. For better calculation of drainage pipe diameter, we can consider that the pipe is full of water, but we have to choose a pipe with an actual diameter greater than the calculated one in order to guarantee the free surface flow inside the drainage pipe. e calculated velocity inside the drainage pipes must be between these limits: And the critical pipe diameter must achieve the following formula: where d φ is the drainage ditch width (m), d cr is the critical pipe diameter (m), K m is the weighted hydraulic conductivity for all layers contributing to the flow (m/day), and q * is the drainage unit discharge (m 3 /day/m), which can be calculated by using the following equation:

Pipe Drainage Optimization
Problem. e optimization problem for the pipe drainage design can be stated as follows.
Subject to in which

Well Drainage Spacing Formulation. According to Soviet
Science Encyclopedia for calculating and design of drainage networks and land reclamation, we can determine the tubewells spacing when the permeability is increasing with depth by using the following formulation [22]: where When the well is not reaching the impermeable layer (m 2 ≠ 0), then where For calculating σ c , we put (h i � h � h * 1 + h 2 + h 3 ), for calculating σ cp , we put (h i � h p � h 2 + h 3 ), and for calculating σ c3 , we put (h i � h 3 ), as described in Figure 6.
We can obtain the Δ ε values from Table 3. And we can obtain β p , β 3 , and β 4 values from Figure 5, after calculating the factors λ, ψ, and ε as seen in Table 4.
In Table 4, r c is the radius of vertical well (m), L k is the length of the filter (m), B is the spacing between vertical wells (m), h sw is the maximum thickness of saturated layers (m), h c is the water depth inside the tubewell above the impermeable layer (m), and q c1 is the drainage unit discharge that must be released by investment wells, and it can be calculated by using the following formula: where p 1 is the discharge that comes from surface water (m/ day), p 2 is the discharge that comes from groundwater (m/ day), φ is the leakage from irrigation network (m/day), and g is the increasing amount of irrigation discharge to avoid salinity problems (m/day).   [22].  Table 4: Calculation of ψ, λ, and ε factors for the values of β 3 , β 4 , and β p [22].
Mathematical Problems in Engineering 5.6. Surrounding Wells Spacing Formula. Protection vertical wells are placed at the edge of the study area along the feeding line L, to protect the study area from the upcoming groundwater along the feeding line. is means that these wells will act as an investment and protection role; the deep groundwater component coming to the study area through the fixed-length L and for a slide of 1 m width is given by the following relationship: where G is the groundwater component (m 3 /m/day) that gathers in the study area during the day, J 1 is the gradient of the groundwater surface at the beginning of the study area (−), and J 2 is the gradient of the groundwater surface at the end of the study area (−). e spacing between surrounding wells can be calculated by using the following formula:

Determining the Number of Drainage Wells.
e total discharge that has to be released by all drainage wells can be determined by e discharge of each drainage well is us, the total number of drainage wells is e number of surrounding wells can be obtained by us, the number of investment wells is

Pump Operating
Hours. e duration of pump operating required to maintain a favourable drainage depth is given by where μ is a specific given factor (−) and R is the radius of influence of the wells, and it can be calculated according to the distribution of the wells as follows.
(i) If the tubewells are placed in a rectangular pattern, R � 0.565 * B (ii) If the tubewells are placed in a triangular pattern, R � 0.526 * B In our study, we will choose a rectangular pattern. Some researchers suggest operating the pumps only during the weeding period, but others prefer to operate the pumps in certain hours every day.

Well Drainage Optimization Problem.
e optimization problem for the well drainage design can be stated as follows. Minimize Subject to (42) Table 5 shows the general parameters calculated for both pipe and well drainage design.

Results and Discussion
By applying these values on the computer code in MATLAB environment, we derived the optimal solution for the pipe drainage design for a range of lifespan as seen in Table 6.
And the optimal solution for the vertical drainage design for a range of lifespan is shown in Table 7.
As we can see, for horizontal and vertical drainage networks, it is better to choose large distance between the lateral pipes and wells. ese distances can be calculated by applying the optimization model on the study area. And for the case study described above, the vertical drainage will be a better solution as subsurface drainage design for the whole Mathematical Problems in Engineering 9 project lifespan. e cost for operating vertical drainage pumps plays an important role in determining the optimal design when considering the lifespan of the project.

Conclusions
In order to formulate an optimization problem for the design of subsurface drainage systems, cost equations have been introduced for both horizontal and vertical drainage. e cost equations contained the most cost components that affect the subsurface drainage networks design. en the optimization problem constraints were derived from the hydraulic study of the case study. e case study contains three different soil layers with different hydraulic conductivity and permeability, and the permeability is increasing with depth. A mathematical model was formulated for the horizontal and vertical drainage optimal design in the case study. e result was a nonlinear optimization problem with nonlinear constraints, which required numerical methods for its solution. A survey of modern optimization algorisms was conducted to find the one suitable for the solution of the formulated problem. It was found that the interior-point optimization algorithm was adapted to the problem and produced satisfactory results. e results show that the proposed optimal mathematical model for both horizontal and vertical drainage networks was affected mostly by the distance between pipes and wells, and the optimal solution involved the maximum possible values of pipes and tubewells spacing. Also, for this case study, the model gave a lower cost for the designing of tubewells network compared with pipe network. And the total cost for the vertical drainage design involved minimum duration of pump operation when considering the lifespan of the subsurface drainage project. e study has shown that the pipes and tubewells spacing and the groundwater table drawdown cannot be selected randomly if we put the economic factor in consideration. Traditional pipes and tubewells design may lead to high costs compared with the optimal design. It is hoped that the proposed optimal mathematical model will present a design methodology by which the costs of all alternative designs can be compared so that the least-cost design is selected.
Data Availability e data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest
e authors declare no conflicts of interest.