SENSITIVITY A N D U N C E R T A I N T Y ANALYSIS OF SUBSURFACE DRAINAGE DESIGN By Guangxi Wu B. A. Sc., Jilin University of Technology, 1984 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF T H E REQUIREMENTS FOR T H E DEGREE OF M A S T E R OF A P P L I E D SCIENCE in T H E FACULTY OF GRADUATE STUDIES T H E DEPARTMENT OF BIO-RESOURCE ENGINEERING We accept this thesis as conforming to the required standard T H E UNIVERSITY OF BRITISH COLUMBIA October 1988 © Guangxi Wu, 1988 In presenting degree at this the thesis in University of partial fulfilment of British Columbia, I agree freely available for reference and study. I further copying of department publication this or of thesis for by his or her representatives. requirements that the for an advanced Library shall make it agree that permission for extensive scholarly purposes may be It is granted by the understood that head of copying my or this thesis for financial gain shall not be allowed without my written permission. Department The University of British Columbia 1956 Main Mall Vancouver, Canada V6T 1Y3 DE-6(3/81) the Abstract Literature on subsurface drainage theories, determination of drainage parameters, and analysis approaches of uncertainty was reviewed. Sensitivity analysis was carried out on drain spacing equations for steady state and nonsteady state, in homogeneous soils and in layered soils. It was found that drain spacing is very sensitive to the hydraulic conductivity, the drainage coefficient, and the design midspan water table height. Spacing is not sensitive to the depth of the impermeable layer and the drain radius. In transient state, spacing is extremely sensitive to the midspan water table heights if the water table fall is relatively small. In that case steady state theory will yield more reliable results and its use is recommended. Drain spacing is usually more sensitive to the hydraulic conductivity of the soil below the drains than to that of the soil above the drains. Therefore, it is desirable to take samples from deeper soil when measuring hydraulic conductivity. A new spacing formula was developed for two-layered soils and a special case of three-layered soils with drains at the interface of the top two layers. This equation was compared with the Kirkham equation. The new formula yields spacings close to the Kirkham equation if the hydraulic conductivity of the soil above the drains is relatively small; otherwise, it tends to give more accurate results. First and second order analysis methods were employed to analyze parameter uncertainty in subsurface drainage design. It was found that conventional design methods based on a deterministic framework may result in inadequate spacing due to the uncertainty involved. Uncertainty may be incorporated into practical design by using the simple equations and graphs presented in this research; the procedure was illustrated through an example. Conclusions were drawn from the present study and recommendations were made for future research. 11 Table of Contents Abstract ii List of Tables ix List of Figures x Acknowledgement xii 1 Introduction 1 2 Literature Review 7 2.1 Subsurface Drainage Design Theory 7 2.1.1 Steady State Subsurface Drainage Theory 9 2.1.2 Non-steady State Subsurface Drainage Theory 2.2 2.3 3 14 Determination and Measurement of Drainage Parameters . 18 2.2.1 Design Criterion 18 2.2.2 Measurement of Soil Parameters 22 Subsurface Drainage Design Under Uncertain Nature . . . 25 Sensitivity Analysis of Drainage Design in Homogeneous Soil 34 3.1 Steady-state: Hooghoudt's Equation 34 3.1.1 Sensitivity of Spacing to Drainage Coefficient R 35 3.1.2 Sensitivity of Spacing to Hydraulic Conductivity 36 3.1.3 Sensitivity of Spacing to Water Table Height Midway Between Drains 37 in 3.1.4 Sensitivity of Spacing to the Equivalent Depth d 38 3.1.5 Discussion on Relationship between Equivalent Depth d and Im- e e permeable Layer depth below drains d and Drains Radius r . . . . 3.1.6 Sensitivity of Spacing to the Depth of Impermeable Layer below Drains 41 3.1.7 Sensitivity of Spacing to Radius of Drains 42 3.1.8 Estimation of Spacing Change due to Changes of All or a Portion of the Parameters 3.1.9 42 Conclusions for Sensitivity Analysis of Steady-state Drainage Design i n Homogeneous soil 3.2 43 Transient State: van Schilfgaarde Solution of Boussinesq Equation . . . . 44 3.2.1 Sensitivity of Spacing to Drainable Porosity 47 3.2.2 Sensitivity of Spacing to Hydraulic conductivity 47 3.2.3 Sensitivity of Spacing to Draw-down Time from Water Table mo to m 3.2.4 48 x Sensitivity of Spacing to the Initial Water Table Height after Recharge Ends 49 3.2.5 Sensitivity of Spacing to the Water Table Height after Recession . 49 3.2.6 Sensitivity of Sapcing to the Equivalent Depth 51 3.2.7 Sensitivity of Spacing to Actual Depth of Impermeable Layer below Drains and Drain Radius 3.3 52 Comparison of Sensitivities between Steady and Transient Drainage Design in Homogeneous Soils 4 40 53 Sensitivity Analysis of Steady State Drainage Design in Multi-layered soil 56 iv 4.1 Sensitivity Analysis of Steady State Drainage Design in Two-layered Soil 56 4.1.1 Drainage Design Using Hooghoudt's Equation 57 4.1.2 Drainage Design Using the Generalized Hooghoudt-Ernst Equation 63 4.1.3 Comparison of Spacings and Sensitivities Obtained from the Hooghoudt and the generalized Hooghoudt-Ernst Equations 4.1.4 66 Conclusions from the Sensitivity Analysis of Steady State Drainage Design in Two-Layered Soils 4.2 68 Sensitivity Analysis of Steady State Drainage Design in Three-layered Soil 4.2.1 Development of A Combined Formula for Spacing Calculation from the Hooghoudt and Ernst Equations 4.2.2 69 Comparison of Spacing between the New Spacing Equation and the Kirkham Equation 4.2.3 73 Sensitivity Analysis of Drainage Design Using the Newly Developed Spacing Formula 4.2.4 68 75 Comparison of Sensitivities of Spacing between the Newly Developed Formula and the Kirkham Equation 77 Uncertainty Analysis of Drainage Design 80 5.1 First and Second Order Analysis of Parameter Uncertainty 81 5.2 Uncertainty Analysis of Steady Drainage Design in Homogeneous Soils . 5.2.1 Drainage Coefficient R as Random Variable 5.2.2 Hydraulic Conductivity as Random Variable With Mean of K and Variance of Vax(K) 5.2.3 86 87 90 Correlated Hydraulic Conductivity K and Drainage Coefficient R with a covariance of Cov(i2, K) 90 5.2.4 Equivalent Depth as a Random Variable with Mean of d and Varie ance of Var(c? ) 93 e 5.2.5 Estimating Variance of Equivalent Depth from the Mean and Variance of the Depth to an Impermeable Layer 5.2.6 96 A l l or Some of the Parameters as Random Variables with Known Means, Variances and Covariances 5.3 97 Uncertainty Analysis of Transient Drainage Design in Homogeneous Soils 5.3.1 Drainable Porosity as Random Variable with Mean of / and Variance of Var(/) 5.3.2 100 Hydraulic Conductivity as a Random Variable with Mean of K and Variance of Var(if) 5.3.3 101 Correlated Hydraulic Conductivity and Drainable Porosity with a Known Covariance of Cov(K,f) 5.3.4 101 Equivalent Depth d as a Random Variable with Mean of d e e and Variance of Var(i ) 103 e 5.3.5 A l l or Part of the Parameters as Random Variables with Known Means, Variances and Covariances 5.4 104 Uncertainty Analysis of Steady Drainage Design in Two-layered Soils 5.4.1 . . . . 106 Hydraulic Conductivity of the Soil above Drains as Random Variable with Mean of K and Variance of V a ^ i ^ ) 108 x 5.4.3 Hydraulic Conductivity of the Soil below Drains as Random Variable with Mean of K2 and Variance of Var(i^2) 5.4.4 105 Drainage Coefficient as Random Variable with Mean of R and Variance of Vai(R) 5.4.2 99 108 Correlation of Drainage Coefficient R and Hydraulic Conductivity, Ki, of the Soil above Drains with covariance of Cov{K vi 1} R) . . . . 110 5.4.5 Correlation of Drainage Coefficient R and Hydraulic Conductivity, K of the Soil below Drains with covariance of Cov(i^2) R) 2 5.4.6 • • • • 112 Correlation of Hydraulic Conductivities of the Soil above Drains, Ki, and below Drains, K , with covariance of Cov(iCi,K ) . . . . 2 5.4.7 114 2 Equivalent Depth as Random Variable with Mean of d and Varie ance of Var( d ) 117 e 5.4.8 A l l or Part of Parameters as Random Variables with Known Means, Variances and Covariances 5.5 119 An Example Showing How Uncertainty Analysis Can Be Used in Drainage Design 6 120 5.5.1 A n Example . 121 5.5.2 Application of Uncertainty Analysis in Drainage Design 122 Conclusions and Recommendations 125 6.1 Conlusions 126 6.2 Recommendations 128 Bibliography 130 Appendices 138 A Discussion on the Simplification of the d — S Relationship 138 e B Estimation of Means and Variances of Random Parameters B.l Methods to Estimate Mean and Variance of Random Variable , 142 B.l.l The Method of Moments . . B.1.2 The Method of Maximum Likelihood 144 B.l.3 Bayesian Estimation of Mean 145 vii ., 142 143 B.1.4 Confidence Interval for Mean and Variance in Normal Distribution B.2 Estimation of the Mean and Variance of Hydraulic Conductivity 146 147 B.3 Application of Bayes Method in the Estimation of the Mean Hydraulic Conductivity 151 Graph of Geometry Factor a for Radial Flow vm 153 List of Tables 3.1 Comparison of sensitivities between steady and transient state drainage design 4.2 : . . . . Comparison of spacing and its sensitivity to individual parameters calculated from the Hooghoudt and the generalized Hooghoudt-Ernst equations. 4.3 74 Comparison of sensitivities of spacing to individual parameters based on the new spacing formula and the Kirkham equation 5.5 67 Comparison of drain spacings in two-layered soils obtained from the new spacing formula and the Kirkham equation (1971) 4.4 55 Probability that water table would not exceed certain height A.6 Comparison of sensitivities obtained from different approaches IX 78 124 141 List of Figures 2.1 Geometry of steady state drainage system 3.2 Sensitivity of spacing to midway water table height in relation with 3.3 Sensitivity of spacing to equivalent depth in relation with ^ 3.4 Symbols used to describe drainage to a ditch or a tube drain 3.5 Sensitivity of spacing to the initial water table height 3.6 Sensitivity of sapcing to water table height after recession. 51 3.7 Sensitivity of spacing to equivalent depth 52 4.8 Sensitivity of spacing to hydraulic conductivity of upper layer soil 58 4.9 Sensitivity of spacing to the hydraulic conductivity of the lower soil layer. 59 4.10 Sensitivity of spacing to water table height 11 • • 38 40 '45 ra 50 0 : . . . 61 4.11 Sensitivity of spacing to equivalent depth 62 4.12 Geometry of the subsurface drainage in three-layered soils. . . . . . . . . 71 5.13 Expected water table height due to uncertainty of drainage coefficient. . 89 5.14 Variance of water table height due to uncertainty of drainage coefficient. 89 5.15 Expected water table height due to uncertainty of hydraulic conductivity 91 5.16 Variance of water table height due to uncertainty of hydraulic conductivity 91 5.17 Expected water table height due to the correlation of the hydraulic conductivity and the drainage coefRcent . 94 5.18 Variance of water table hieght due to the correlation of the hydraulic conductivity and the drainage coefficent x 94 5.19 Expected water table height due to uncertainty of equivalent depth 5.20 Variance of water table height due to uncertainty of equivalent depth ... . . 95 95 5.21 Expected water table height due to uncertainty of the drainage coefficient. 107 5.22 Variance of water table height due to uncertainty of the drainage coefficient. 107 5.23 Expected water table height due to uncertainty of the hydraulic conductivity of the soil above the drains 109 5.24 Variance of water table height due to uncertainty of the hydraulic conductivity of the soil above the drains 109 5.25 Expected water table height due to uncertainty of hydraulic conductivity of the soil below drains. Ill 5.26 Variance of water table height due to uncertainty of hydraulic conductivity of the soil below drains. Ill 5.27 Expected water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil above drains 113 5.28 Variance of water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil above drains 113 5.29 Expected water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil below drains 115 5.30 Variance of water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil below drains 115 5.31 Expected water table height due to correlation of hydraulic conductivity of the soil above drains and that below drains 116 5.32 Variance of water table height due to correlation of hydraulic conductivity of the soil above drains and that below drains. 5.33 Expected water table height due to uncertainty of equivalent depth. . . . 116 118 5.34 Variance of water table height due to uncertainty of equivalent depth. . . 118 xi C.35 Geometry factor a for radial flow (taken from van Beers, 1976) Xll Acknowledgement I would like to express my special gratitude and appreciation to Dr. S. T. Chieng for introducing me to this topic and for his guidance, help and encouragement throughout the present study. Appreciation is also extended to Dr. W . F. Caselton and Dr. K . V . Lo for sitting on the committee, providing valuable advice and reviewing this thesis. The University of British Columbia and the Canadian people are gratefully thanked for the support they provided in the award of the University Graduate Fellowship during the second year of the program. Acknowledgement is extended to all professors, staff and students in the Department of Bio-Resource Engineering for their kindness and help that have been making study and life here easier for me. Special thanks are due to Linda for editing this thesis. I would like also to express my sincere appreciation and feeling of indebtedness to the Chinese people for the great sacrifices they have made to provide me with better education and the State Education Commission Scholarship by which I was supported to study in the University of British Columbia during the first year of the program. M y thanks are also due to my fellow Chinese students and visiting scholars in this university for their help and spiritual encouragement. I must thank my parents and my brothers for their care, love, understanding, encouragement and moral and spiritual support. I am specially indebted to my parents for the hard times they have experienced in supporting and educating their six children. In the end, I wish to express my special gratitude to my wife, Linan, for her wholehearted love, understanding and encouragement. Without her support, this thesis would not have been finished. xm Chapter 1 INTRODUCTION As a key component of the modern agriculture, drainage is playing, and will continue to play, an important role in the food production system. With the rapidly increasing world population, especially in those less developed countries, the strong demand for food and fibre has been putting greater and greater, and at times even explosive, pressures on agriculture. To increase food and fibre production, there are two basic approaches: (i) to expand cultivated area, and (ii) to increase production per unit area by ameliorating the agricultural conditions. Although the scientific and technological advances will contribute greatly to production increase, agricultural drainage will play a very important or even dominant role. By means of drainage, crop yield will be increased by more intensive and efficient use of cultivated land; more and more marshland and other waterlogged lands, estimated at over 3.5 million square kilometers (Bulavko, 1971), could be drained and put into agricultural production to meet the steadily increasing demand for food and fibre. The primary purpose of agricultural drainage is to remove excess water from the root zone for profitable food production. In humid areas, excess water mainly comes from intensive rainfall and sometimes from surplus irrigation, which may be needed as an anti-frost protection during the frost periods, or in other areas, to supplement and offset the lack of timely or regular rainfall. In irrigated areas, drainage is needed to remove the excess water for salinity control. The most important design parameters of a subsurface drainage system are drain 1 Introduction 2 depth and drain spacing. Drain depth depends mainly on the soil profile, the installation equipment, and the discharge outlet. Once the drain depth has been determined from actual field conditions, drain spacing is the only design parameter which should be determined from climatic conditions, soil properties and crop requirements. In cases where several drain depths are possible, the combination of drain depth and drain spacing should be such that the total cost of the system (including installation, material and maintenance costs) is minimized. At present, steady state theory is commonly used in agricultural drainage because of its simplicity and practicability in most situations. In cases where the water table fluctuations must be considered, more complex and relatively more realistic transient theory should be used. In both steady and transient state theories, drainage design is based on a deterministic framework, in which parameters such as soil hydraulic conductivity and drainable porosity are represented by some estimated or equivalent values. However, such soil properties, along with climatic conditions, tend to be both spatially and temporally variable. Their representation by equivalent values is also ambiguous, if not totally misleading. Such equivalent values are selected partly based on the actual nature of the system, but also largely on the experience and professional judgement of the individual designer, which adds to the variability and uncertainty inherent in these properties. Drain spacing can be calculated from certain drainage models or equations, such as Hooghoudt equation for steady state and Boussinesq equation for transient state, provided adequate information is available for various parameters. Any differences in the estimation of related parameters will eventually change the designed spacing, which in turn affects the economy of the drainage system. In the case of a large well-financed enterprise, the initial installation costs may not be the major concern, but for a farmer who has only limited funds to expend on a drainage facility and whose primary purpose is to minimize the risk of severe crop losses in the event of adverse weather, the initial Introduction 3 cost of a drainage system and the uncertainty about its performance may result in the difference of accepting or declining the drainage design. The best drainage system design is the one that minimizes cost, while maximizing the monetary return or minimizing the risk of severe crop loss. In practice this usually can not be achieved and a compromise is required. A decrease in cost may result in poorer system performance. A particular design is often dependent on the combination of the amount of money a farmer is willing to expend on drainage system and the risk level of crop failure in the event of adverse weather. Although many methods are available to determine soil hydraulic properties such as soil hydraulic conductivity and specific yield, these parameters are determined with certain amount of error, due to inherent varibility of the soil and the limitations of accuracy of the methods. For a particular field, the values of the parameters estimated are directly related to the locations and number of measurements performed, which are usually subject to the experience and professional judgement of individual designers. Even after the necessary data have been collected, different methods of interpreting the data may also result in different estimates. While a criterion to decide which estimate is better than others is difficult to set up, it is possible to estimate the effects and uncertainty of each estimate on the drain spacing and performance of a drainage system. Among the numerous drain spacing formulae (Kirkham, 1966; Mueller, 1967), the Hooghoudt equation is most commonly used in steady state drainage design owing to its simplicity compared with other more complex models. In cases of transient flow, many spacing formulae are obtained by simplifying the governing partial differential equation of groundwater flow, which is obtained by combining Darcy's law with the law of conservation of mass. The Boussinesq equation is developed by applying Dupuit- Forchheimer theory (van Schilfgaarde, 1974). Various linearizations have been applied to solve the Boussinesq equation (Dumm, 1954; Werner, 1957; Krayenhoff van de Leur, 1958; Introduction 4 Maasland, 1959; Polubarino-Kochina, 1962; Dumm, 1964; Terzidis, 1968). Most of the solutions obtained from linearizations, for example, the Glover equation, are restricted to small increment of water table drawdown and time periods. However, van Schilfgaarde (1963,1964) solved the non-linear partial differential equation without linearization, and the van Schilfgaarde solution is not restricted to relatively small increments of time and water table. In addition, the van Schilfgaarde solution is relatively simple and can still provide realistic results. For the above reasons, the van Schilfgaarde equation is used in this study for transient drainage design. With only minor modification the Hooghoudt equation can also be used in cases of two-layered soils, if the drains are at the interface of the two layers (van Beers, 1976). This is because the Hooghoudt approach primarily distinguishes the flow regions of above and below drains, and only secondarily distinguishes soil layers. However, it is claimed that the Hooghoudt equation can no longer be used in cases where a heavy clay layer of varying thickness overlies a sandy substratum (K\ <C K ), and where there are two 2 distinct pervious layers in the soil below drain level (van Beers, 1976). In these cases the use of the Ernst equation is recommended (van Beers, 1976). On the other hand, the Ernst equation is not suitable if a highly pervious layer is above a poorly pervious soil layer (Ki K ) (van Beers, 1976). A generalized equation applicable to all these 2 cases was obtained by Ernst (1975), who combined the Hooghoudt equation and the Ernst equation for radial flow. But it requires to solve a high order equation, which is sometimes inconvenient in practice. Kirkham (1971a) spacing equations for layered soils neglect the head loss in the region above drains and tend to underestimate the drain spacing if there is a highly pervious soil layer above drain level. A new spacing formula is developed in the study by combining the form of the Hooghoudt equation and the concept of equivalent depth with the radial flow component in the Ernst equation. This new spacing formula can be used in subsurface drainage design of two-layered for any Introduction 5 conditions and a special case of three-layered soils. A major objective of this study is the sensitivity and uncertainty analysis of subsurface drainage design based on commonly used models and approaches. Sensitivity analysis reveals the relative importance of various parameters in a model to the spacing of a drainage system, which can be used as a guidance in soil surveys and field geohydrologic investigations. Uncertainty analysis can be used to evaluate the performance of a drainage system in a particular field and to adjust drain spacing according to the risk level a farmer is willing to accept. The method employed in the uncertainty analysis is the first and second order analysis of parameter uncertainty. It requires only information about the first two moments (mean and variance) of a random variable, which is much easier to obtain than the complete distributions of the variables (Cornell, 1972). It also requires only basic mathematic and statistical knowledge. These two features provide a strong base for its usefulness in practical drainage design. The primary aims of this study are to apply simple but useful tools of uncertainty analysis (sensitivity analysis and first and second order analysis) to identify the dominant parameters in subsurface drainage system design, and to provide practicing engineers/designers with ready-to-use formulae and graphs to obtain an overall estimate of the uncertainty involved in the design due to the uncertainty and variability of the parameters. This thesis is organized in the following way: in Chapter 1 (this chapter), an introduction for this study is provided. In Chapter 2, various approaches of subsurface drainage design are reviewed. Methods currently used in the determination of drainage parameters including the development of drainage criteria and the measurement of hydraulic conductivity are discussed. Techniques of incorporating uncertainty into the study of groundwater flow problems and drainage-design are reviewed and discussed. In Chapter 3, the sensitivity analysis of drainage design in homogeneous soils is carried out. The Introduction results from selected steady and transient theories are compared. 6 In Chapter 4, the sensitivity of spacing to various parameters of drainage design in multi-layered soils is analyzed; several models are compared and a new spacing formula is then developed. In Chapter 5, the uncertainty analysis of subsurface drainage design in homogeneous and multi-layered soils is discussed. A n example is presented to illustrate how uncertainty analysis can be used in practical design. In Chapter 6, a summary and the conclusions from this study is presented. Recommendations are made for future research in this area. Appendix A includes a discussion on the relationship between equivalent depth and drain spacing. In Appendix B, some methods of measuring and estimating the means and variances of soil properties .are reviewed. Chapter 2 LITERATURE REVIEW Over the centuries, drainage has evolved from very primitive forms to modern state-ofthe-art technology. Extensive literature is available on this and related areas. Current research tends to consider drainage conditions more realistically, such as efforts on applying stochastic approaches to drainage design. While it is extremely difficult, if not impossible, to incorporate natural conditions exactly into design, it is relatively easier to estimate the effects of uncertainty about drainage conditions on drainage system design and its performance under specific situation. This chapter reviews the literature on major subsurface drainage design theories, uncertainty analysis approaches, and drainage parameter measurement and determination methods. 2.1 Subsurface Drainage Design Theory There are generally two ways to design a drainage system: (a) from empirical data collected through evaluation of existing drainage systems, and (b) from a theoretical analysis of the problem, applying known physical laws and tested theories. Criteria developed by empirical methods are based substantially on experience and assessment of numerous interrelated factors, and special care must be taken in transposing their use from one locality to another. Theoretical methods are more commonly used in practice: applying proved principles or laws to problems having known limiting conditions. The resulting mathematical 7 Literature Review 8 expression explains the observed action of existing drainage systems and permits the rational design of new systems as long as the site conditions are within the limits for which the equation was derived. For example, the Hooghoudt equation (1940) was originally derived for the conditions in the Netherlands, but it is also widely used in Europe and North America (Chieng, 1978). Drainage design is usually based on one of the numerous spacing formulae with its parameters evaluated from the local conditions. These spacing equations may be derived from two concepts: Darcy's Law, and the equation of continuity. For spacing equations based on the potential theory, these concepts are combined into Laplace's partial differential equation and when this equation is solved for the conditions of the particular problem, the solution yields, simultaneously, spacing equations and detailed information about the flow system. However solutions to Laplace's partial differential equation are usually very difficult to obtain for various initial and boundary conditions. On the other hand, for spacing equations by Dupuit-Forchheimer (D-F) theory (Dupuit, 1863; Forchheimer, 1930), simplifying assumptions can be made along with the use of Darcy's law and the equation of continuity to obtain useful results without solving Laplace's equation. However these results do not give the internal potential of the flow system. Drainage spacing formulae may be generally categorized as either steady state formulae or non-steady state formulae. In practical drainage design, soils are usually assumed to be homogeneous,(or homogeneous within one layer), isotropic, and only saturated flow is considered. Studies on water flow through heterogeneous soils and in the unsaturated zone are more for research purposes than direct application. In the remainder of this section, the literature on steady and non-steady subsurface drainage design will be reviewed. Literature Review 2.1.1 9 Steady State Subsurface Drainage Theory Steady state spacing formulae are based on a state of equilibrium between recharge and discharge, i.e. discharge equals recharge and the hydraulic head remains constant. In order to apply these formulae, a drainage requirement or design criterion is required, describing the quantity of water to be removed in a given period of time under conditions of a given height of the groundwater table. Steady state conditions are most likely to occur in flat areas of humid regions, where the need for drainage is determined by long duration rainfall with low intensity in winter and spring. Miiller (1967) lists 40 steady state spacing formulae and compares results of calculations with the optimal drain spacings as determined at 35 experimental fields. He concluded that only the formulae of Hooghoudt (1940), Ernst (1956), and Toksoz & Kirkham (1961) are suitable for general practical application. Hooghoudt (1940) considered only horizontal flow in the region above drains and both horizontal and radial flows in the region below drains. The radial flow caused by the convergence of flow lines near the drains is taken into account by reducing the depth of the flow layer (D) below the drains to the hypothetical depth (d ) of an "equivalent e layer", which is a function of (D), spacing (5), and drain radius (r). The Hooghoudt equation can be written as: 2 5 _ 8K d h -~~R~~ 2 e %K D h ~ R ~ l + l in which S = drain spacing (m); K 1 K 2 =• hydraulic conductivity of the soil above drains (m/day); = hydraulic conductivity of the soil below drains (m/day); h = water table height above drains midway between two drains (m); D\ = average cross-section of flow area above drains (m); Literature Review 10 R = drainage coefficient or drainage discharge rate (m/day); d =thickness of the "equivalent layer". e Figure 2.1 is a sketch to show the symbols used in the analyses. The equivalent depth d can be calculated by (Moody, 1966): e For 0 < d/S < 0.3: * l + !(!lnj!-a) a = d .d . , 3.55-1.6-+2(-) d in which 2 For d/S > 0.3: 4 -f(lnf-1.15) where d is the depth to the impermeable layer, and r is the drain radius. The Hooghoudt equation can be used for homogeneous soils (K = i^i), but also can 2 be used for a drainage situation with two layers of different permeability (K 2 ^ Ki), if the drain level is at the interface of these layers. This is because Hooghoudt's approach distinguishes primarily the flow regions of above and below drains, and only secondarily distinguishes the soil layers. Wesseling (1964) carried out extensive calculations using potential theory and compared the results with the Hooghoudt formula. The comparison shows that results obtained from using Hooghoudt's equation do not vary more than about 5% from those of potential theory for midspan water table heights. Since its appearance, Hooghoudt's equation has been favoured by practicing engineers for design in humid areas because of its accuracy, simplicity of form and ease of use. However, the Hooghoudt equation is not applicable when there is more than one layer soil below the drains and or the soil above drains has a much lower hydraulic conductivity (van Beers, 1976). In these cases, use of the Ernst equation (1956, 1962) is recommended. Literature Review Figure 2.1: Geometry of steady state drainage system. \\ Literature Review 12 Ernst (1956, 1962) divided the groundwater flow to parallel drains, and consequently the corresponding available total hydraulic head (h), into three components: a vertical (v), a horizontal (h), and a radial (r) component, or h = h + h + h = qL + qSL v h T v h + qSL r where q is the flow rate and L is the resistance. Working out the various resistance terms and changing notation of q for drainage rate to R for consistency, we can write the Ernst equation as: „D„ „ S „ S , aD h = R-^- + i ? — — - + i i - — - In K 8KD TK U 2 2 V 2 where D = v thickness of the layer over which vertical flow is considered; in most cases this component is small and may be ignored (m); K v = vertical hydraulic conductivity (m/day); KD = the sum of the product of the permeability (K) and thickness (D) of the various layers for the horizontal flow component according to the hydraulic situation: a = one pervious layer below drains: KD = K-^D-^ + two pervious layers below drains: KD = K D\ + K D X X + KD 2 2 2 2 + KzDz X 2 a = 1 2 KD = K D\ + K D = 2 geometry factor for radial flow: KD = K D u X KD 2 + KD 3 3 a can be estimated from Graph I wetted section of the drain (m); for pipe drains u = rr All other parameters are same defined as in the Hooghoudt equation. Literature Review 13 The basic difference between the Hooghoudt and Ernst equations is that Ernst took radial flow into account for the total flow, whereas Hooghoudt considered radial flow only for the flow below drains. For a drainage situation where K <C K (e.g. a clay layer on 1 2 a sandy substratum), the Ernst equation can provide accurate results. However, for a situation where K\ K (e.g. a sandy layer on a clay layer), the results obtained from 2 the Ernst equation tend to considerably underestimate drain spacing compared with the results obtained from the Hooghoudt equation. According to Ernst (1962) and Van Beers (1965), no acceptable formula has been found for the special case where K\ ^> K , and 2 the Hooghoudt equation was suggested to be used for this case. Efforts have been made to develop an equation which can be used for all K jK x 2 ratios. Such an equation was obtained by combining the approach of Hooghoudt (horizontal flow only above drains and both horizontal and radial flows below drains) and the expression for the radial flow component in the Ernst equation (Ernst, 1975). Neglecting the vertical resistance, we have: ZKiDJi S %K D h 2 2 + S 2 +^ 7T 2 2 l n ^ U By introducing an equivalent drain spacing (S ), i.e. a drain spacing that would be 0 found if horizontal flow alone is considered: _8KDh °~ R 2 b we can write the above equation as: SKDh S _ SKjDjh S 2 2 0 8K D h S + ^D ln^ 2 + 2 2 2 Rearranging the above equation gives the generalized Hooghoudt-Ernst equation: ( A)3 _8^ A _ A ( + S 0 TTS 0 )2 S Q S 0 = TTS 0 0 Literature Review 14 where s = drain spacing based on both horizontal and radial flow (m); So = drain spacing based on horizontal flow only (m); c = £ ) l n ^ , a radial resistance factor (m); B = KiDi KD 2 2 Graphs have been prepared to use this generalized equation (van Beers, 1976). Toksoz & Kirkham (1971a) used potential theory to develop an equation for twolayered soil with drains in the upper layer. The assumption used by Toksoz & Kirkham to simplify the field problem was that the loss in hydraulic head in the region lying below the arch-shaped water table and above the level of the drains is negligible compared with the head loss over the remainder of the flow region. This assumption is acceptable when the horizontal hydraulic conductivity in the soil above the drains is small and when the water table height above the drains is small compared with the depth of the soil profile below the drains. Otherwise the Kirkham equation tends to underestimate the spacing required. The Kirkham equation is also somewhat complicated and its solution requires a large number of nomographs. However, the Kirkham equation provides internal potential and stream functions which might be useful for further reserch. 2.1.2 Non-steady State Subsurface Drainage Theory Realistically speaking, steady state conditions are seldom, if ever, encountered in the field. Where high intensity rainfall is common and where the field is irrigated, non- steady state groundwater flow theory should be used in the design of subsurface drainage systems. Non-steady, or time-varying, flow problems are more difficult to solve than time invariant problems. For the viscous flow of fluids through porous media, classical theory Literature Review 15 leads to the Navier-Stokes equations (van Schilfgaarde, 1970). However these equations are very difficult to solve in many situations of practical interest and therefore have little direct application. Hence, Darcy's law is applied to replace the momentum balance to obtain more applicable equations. Darcy's law, originally stated for one-dimensional flow through a sand filter, may be generalized to q= -KVh where q is the rate of flow through the media, h is the hydraulic head, and K is the hydraulic conductivity of the media. Combined with the equation of continuity for saturated, homogeneous, isothermal and incompressible systems, Darcy's law leads to the Laplace equation: Vh 2 =0 Since time does not appears in this equation, to apply it to non-steady state problems the time variable must be introduced by means of the time-variant boundary conditions . However, the determination of such boundary conditions is very difficult, and it is assumed in the derivative of the Laplace equation that the soil is saturated, i.e. no water storage in the flow zone, which is obviously not the case in an actual field. Another equation that has been widely used for non-steady state groundwater flow problems, is based on Dupuit-Forchheimer theory. The basic D-F assumptions are that in a system of gravity flow towards a shallow sink all streamlines are horizontal and that the velocity along these streamlines is proportional to the slope of the free-water surface but independent of the depth. In addition, it is implied that the free-water surface is a streamline bounding the flow region and, in cases of an open ditch or a dam, that the water table terminates at the water level in the ditch (van Schilfgaarde, 1974). Literature Review 16 Based on D-F assumptions, if a constant specific yield or drainable porosity of the soil / is assumed, the partial differential equation governing the groundwater flow towards parallel ditches or drains, which is often referred to as the Boussinesq equation, can be written as *!<£> *-'! + where R is the net recharge rate (or drainage rate). Again, this equation is limited to saturated flow conditions because a constant drainable porosity is assumed. However, drainable porosity is far from constant, depending at least on the proximity of the water table to the soil surface and on the rate and direction of water table movement. The Boussinesq equation can be used to study the rate of fall of the water table from some initial conditions after the cessation of precipitation or irrigation, or to study fluctuation of the water table as recharge input varies with time. The Boussinesq equation is non-linear and difficult to solve. Various linearizations have been made to solve this equation (Dumm, 1954; Werner, 1957; Krayenhoff van de Leur, 1958; Maasland, 1959; Tezidis, 1968). Glover, as reported by Dumm (1954), obtained a solution to the Boussinesq equation without source of recharge in a Fourier series: , , 4m A , 1 . , .n-Kx. . h-d = —J2j-) n{ — )exp( 0 ta n iv KDt. 2 _ 2 ) starting from an initial flat water table, and the boundary condition: for 0 < x < S at t = 0 h = d for x = 0, S at t > 0 dh/dx = 0 for x = 5/2 t > 0 h = h = d+ m 0 0 If time t is not too small, all but the first term can be neglected, resulting in the Literature Review 17 simplified expression of the midspan water table height (x — 5/2) 4m , -K KDt. _exp(—^j-) 2 m i 0 = or „, Tt KDt 2 / l n ^ which is referred to as the Glover equation. Here, D = d + m /2 e 0 where d is the e equivalent depth. Tapp and Moody (Dumm, 1964) observed that the initial water table shape encountered in the field was more often parabolic than flat. They assumed a fourth degree parabolic initial water table shape of the form: h-d = ^(S x 3 - 3S x + ASx 2 2 3 2x ) A and h = 0 for x = 0, S at t = 0. The solution obtained thereafter is: , . 192m ^ i - ^ E 0 V-8, . .nvx. . (-^-)-(-^)exp( ra 5 n iv KDt. j ^ - ) 2 2 For all but the smallest time periods, all terms but the first are negligible, and the equation becomes: ~K KDt 2 n 22 S = /In 3.7m •Km i which is only a minor modification of the Glover equation. This modified Glover equation is currently recommended for use by the Bureau of Reclamation. However, it is restricted to small increments of water table recession. The Boussinesq equation has been solved without linearization for situations where the drains rest on an impermeable layer (d = 0) (Boussinesq, 1903; Dumm, 1964), and for situations where the drains do not rest on an impermeable layer (van Shilfgaarde, 1963). The van Schilfgaarde solution is: Literature Review 18 g 2 9Kd t = e r mo(2d +Tn ) •> mi(2d +m ) e 1 e 0 where d in the original equation is substituted with equivalent depth, d , to circumvent e the D - F restriction that d<^ S. If a water source is presented at the soil surface, an analytical solution becomes very difficult to obtain. In most cases, numerical techniques are employed instead to obtain approximate solutions (Moody, 1966; Skaggs, 1975). 2.2 Determination and Measurement of Drainage Parameters A successful drainage system design is often more dependent on the investigation and evaluation of local conditions than on groundwater flow theory. This includes the study and analysis of the weather conditions, crop response to water table height and soil properties. Although costs in this stage of design may be very high, it is usually worthwhile for the long term economic returns. 2.2.1 Design Criterion The design criterion is the difference between the drainage requirement, which is the total desired drainage intensity for a given region or field, and the existing natural drainage intensity of a given field. Development of a drainage design criterion is an indispensable part of the design process. Drainage criteria are developed by local experience or by analysis of soil water balance (Willardson, 1982). As far as the subsurface drainage system costs are concerned, the design drainage criterion should be no larger than necessary to remove excess water adequately and quickly enough to meet crop and trafficability requirements (Chieng, et al, 1978). The drainage criterion can be evaluated for a steady Literature Review 19 state condition, a non-steady state condition, salinity control, or soil trafficability, depending on the main purpose of the drainage system, which in turn depends on the climatic pattern and crop characteristics. Steady State Condition For steady state drainage design, the drainage requirement depends essentially on the rainfall pattern and the effects of water table position on crop yield. However, this effect is difficult to estimate because of the complex nature the crop responses to a number of factors. For areas where there are a considerable number of drainage systems installed, perhaps the most effective way to evaluate drainage criteria would be to observe both farm and experimental systems of different drainage intensities. Chieng et al (1978), on the other hand, developed a water balance computer model using a long period of weather record to estimate the minimum drainage rates that would be acceptable. For irrigated areas, the drainage criterion is evaluated from application intensity and uniformity, and evapotranspiration. Non—steady State Condition In areas where high intensity rainfall or irrigation is common, the water table may rise very rapidly into the root zone and even to the soil surface. In these cases, drainage criteria are usually expressed as the rate of fall of water table at a certain water table position, or a certain distance of fall from one water table to the next in a given period of time. The falling water table drainage criterion depends on the tolerance level of crops to temporary high water tables, which in turn depends on crop type, stage of growth, weather conditions, antecedent water table positions, soil fertility level, soil porosity and Literature Review 20 other factors (Bouwer, 1974). Observation and experience with local conditions is usually required to estimate the necessary rate of fall of the water table for various crops. A more realistic view of the situation in the field is a fluctuating water table. If the effect of the water table hydrograph on crop yield is known, and if this hydrograph for a drainage system can be predicted, undesirable patterns of water table fluctuations may be prevented. Salinity Control In arid areas, whether or not irrigated, salt may concentrate in the root zone and even at the soil surface because of the high rate of evapotranspifation. Under these conditions, drainage systems must be installed to maintain the salt balance so that the field can be profitable for a long period of time. The drainage system should be able to prevent high water tables in the growing season when leaching of salts is desired, and should keep the water tables sufficiently low between growing seasons to minimize evaporation from the water table and the resulting salt accumulation in the root zone. Drainage criteria for salinity control are based on the salt balance in the root zone. Excess irrigation water is applied to leach out the salt to maintain a salt balance suitable for crop growth. The amount of leaching water required depends on the salinity level of the irrigation water, salt content in the soil solution, evapotranspiration, irrigation uniformity, deep seepage, and other factors. Successful design of a drainage system for salinity control depends greatly on the accuracy of estimation of these factors. Trafficability Control In humid areas, heavy precipitation at inopportune times may seriously interfere with the timeliness of farming operations. Farming operations, such as planting and harvesting, may be impossible due to poor trafficability in the field. The resulting losses may vary Literature Review 21 from complete crop failure, if planting is delayed too long, to reduced crop yields if tillage, weed control, or harvesting are not performed on time (Reeve and Fausey, 1974). The design drainage criteria of drainage systems for the purpose of trafficability control are often based on providing a certain number of consecutive workable days. Skaggs (1978) defines a workable day as a day when the air volume in the soil profile exceeds some limiting value, the rainfall occurring on that day is less than a minimum value, and a minimum number of days have elapsed since that amount of rainfall occurred. The probability of a certain number of "dry" days at a given time can be predicted from the historic weather data of the region, and artificial drainage can be installed to increase the probability of a number of workable days at a given time of year. In some cases, considering the high cost of installing a subsurface drainage system, an alternative to drainage may be to more farm equipment so that planting or harvesting can be completed in shorter period of time, Soil trafficability depends on the moisture content of the top soil. Oosterbaan and Wind (1979) indicated that most farming operations can be carried out if the tension in the top soil is 100 cm of water or more. However, in practice, agricultural engineers tend to use some simple indicators for workable days, such as depth to the water table, in the design of subsurface drainage for trafficability control. From the above review of drainage criteria, it is obvious that the steady state drainage criterion is the simplest one to use. Therefore, it may be advantageous and also possible to express all drainage criteria in terms of an equivalent steady state criterion, regardless of whether the primary purpose of the drainage system is to yield a certain fall or fluctuation of the water table, salinity control or an extended workable period (Bouwer, 1974). Literature Review 2.2.2 22 Measurement of Soil Parameters Another factor that is as important as the development of the design drainage criterion, is the determination of soil properties. Physical soil properties which are of significance in subsurface drainage design or in drainage research and investigation are the hydraulic conductivity (saturated or unsaturated), the drainable porosity or specific yield, the soil-water pressure head, the volumetric water content, the water table position, etc. In practical drainage design,, only saturated hydraulic conductivity and drainable porosity are of interest; the other information is needed only for further research and detailed investigations. Measurement of Hydraulic Conductivity Hydraulic conductivity appears in both steady and non-steady state drain spacing calculations. The performance and costs of the drainage system depend greatly on the value of hydraulic conductivity. The determination of the K value of the soil is thus an important aspect of almost any drainage investigation. The value of hydraulic conductivity measured from the field or lab should represent as adequately as possible the soil and flow direction in the actual drainage system. Therefore the method, location, and number of measurements should be carefully selected so that adequate data can be collected cost- effectively. Hydraulic conductivity can be measured in the laboratory on soil samples taken from the field. These samples should be representative of the site under investigation. Due to soil variability, a large number of undisturbed core samples should be taken from representative locations with 2-3 replicates at each location. Laboratory methods are mainly used to determine the K value above the water table, to separate the vertical and horizontal hydraulic conductivities, and to predict how salt affects the hydraulic Literature Review 23 conductivity and other physical properties. Field measurements of K are generally preferred to laboratory measurements. If the water table is sufficiently high, one of the auger-hole techniques can be used. If a water table is not present, or is too low for meaningful K value measurement, the techniques for measuring K in an artificially wetted soil region may be employed. Reliable methods have been developed for measurement below the water table, whereas methods employed above the water table are generally less satisfactory (Smedema & Rycroft, 1983. Land Drainage). With the auger-hole techniques, a hole, ranging in diameter from a few centimeters to a few decimeters, is dug to some distance below the water table. After the water level in the hole has equilibrated with the water table in the soil, the water level is lowered by quickly removing water from the hole. Then the K value can be measured from the rate of rise of the water level in the hole (van Beers, 1970). The field methods for measuring hydraulic conductivity in situ above or in the absence of a water table include the shallow well pump-in technique, the cylinder permeameter method, the infiltration gradient technique, the air-entry permeameter, and the doubletube method. For all these techniques, the soil should be as close to saturation as possible and the clogging of pores where water enters the soil should be minimized for relatively reliable measurement of K value. The selection of measuring method in the field depends on the understanding of the field conditions, including the size of the soil region on which K is to be measured; the flow direction, and the presence of stones in the soil. It also depends on the manpower, equipment, time and funds available. However, even with a reasonable selection of measuring method, a certain degree of inaccuracy is inevitable due to the soil variability, which is a far greater source of inaccuracy in evaluating the K value of a field soil than the errors of inherent to the measuring techniques. It is expected that a more accurate Literature Review 24 K value can be obtained with a larger number of measurements. The number of measurements that should be taken to characterize the K of a given field depends on the soil variability, the accuracy desired, and the availability of funds and manpower. Due to other sources of uncertainty, such as entrance resistance of the drains, and water sources and sinks, errors in K of 30% or even more are usually acceptable in drainage design (Bouwer Sz Jackson, 1974). The measurement of hydraulic conductivity should be combined with soil surveys, so that the best locations for K measurement can be determined. It may also be possible to establish a correlation between hydraulic conductivity and soil texture and structure. The locations of measurements should be representative of the soil types. It is also desirable that the sampling area be divided into small subareas according to the soil types so that the K value may be better representative of the actual condition in each subarea. The "true" hydraulic conductivity can be better estimated by the geometric mean than by the harmonic mean and arithmetic mean (Bouwer, 1969). The geometric mean is defined for equally thick layers as: Kg = \JK K • • • K X where Ki, K , ..., K 2 n 2 n are the field data (Bouwer & Jackson, 1974). Measurement of Drainable Porosity The drainable porosity, or specific yield, is the volume of water that will be released from a unit volume of soil by lowering the water table. In non-steady state subsurface drainage design, drainable porosity is considered as a constant, while in reality it is far from constant, varying with time and space. It also depends on the rate and direction of fall of the water table. Literature Review 25 The technique of estimating drainable porosity from equilibrium water content profiles above two successive water table positions, using the relation between soil water content and negative pressure head, may give a reasonable estimation if the water table falls slowly enough for drainage of the pore space to "keep up" (Talsma & Haskew, 1959; French & O'Callaghan, 1966). A negative pressure head above the water table is often called tension or suction and can be measured with tensiometers. Soil water content can be measured by a gravimetric method, Neutron scattering method, or other methods reviewed by Gardner (1965) and Cope & Trickett (1965). The drainable porosity determined from the soil water content characteristics becomes less valid when the rate of fall of the water table increases. If the rate of fall of the water table is relatively high, the value of drainable porosity to be used in the nonsteady spacing equations is only a fraction of the value determined from the equilibrium water content distributions. Therefore, the value of drainable porosity used in subsurface drainage system design should be chosen cautiously. 2.3 Subsurface Drainage Design Under Uncertain Nature From the previous section, it can be seen that neither the drainage criterion nor the soil properties are known with complete certainty. It suggests that any drainage system design based on uncertain input information will result in uncertain performance, and sometimes even complete failure. Aware of this, many scientists and researchers have been looking for methods of incorporating uncertainty into the study of groundwater flow and drainage design. The falling water table criterion is a primitive step in this direction. Then a more realistic model of fluctuating water table or intermittent recharge was studied. However, they are still far from the stochastic nature of the groundwater flow in the subsurface drainage problems. Thus, more advanced and complex models and Literature Review 26 approaches have been developed. Unfortunately, at present, there are still very few, if any, techniques of satisfactory application in practical designs. Uncertainties involved in subsurface drainage design can be categorized generally into four groups: (1) forcing function uncertainty, (2) initial condition uncertainty, (3) boundary condition uncertainty, and (4) parameter uncertainty, all of which are usually presented in a specific problem (Sagar, 1978a). Forcing function uncertainty includes the unknown nature and presence of sources and sinks in the groundwater flow region. It is one of the major sources of uncertainty. Initial condition uncertainty can be minimized for the long time period. Unsaturated flow occurring above the water table is a source of boundary condition uncertainty. Parameter uncertainty is caused by the reduction of the real, but unknown, values (which are stochastic in nature), to "equivalent" values, as well as by some other sources of errors such as instrumental, reading and round off errors. The "equivalent" value of a parameter is very often ambiguous because it may not represent the "true" value of the parameter. In addition to the above uncertainties, there is also a modeling uncertainty. Because of the complexity of the interaction between soil and water, the laws or theories on which drainage and groundwater flow theories are based, are simplified and idealized. Darcy's law, for example, is only an empirical equation, not a theoretical one. Hence the derived equations inevitably involve a certain degree of modeling uncertainty. More uncertainty may enter during the process of derivation by, for example, the Dupuit-Forchheimer assumptions and linearizations. One obvious step to incorporate uncertainty into the study of groundwater flow is to consider the stochastic precipitation and evapotranspiration process, just as many researchers and scientists have done. In subsurface drainage design, net recharge (recharge minus discharge) is the forcing function and is directly related to the fluctuation of the groundwater table. Both precipitation and evapotranspiration are heavily dependent Literature Review 27 on the local meteorological conditions. Hence the uncertainty of the weather condition will result in uncertainty of the groundwater table position and therefore of the performance of a subsurface drainage system. However, analytical solutions to stochastic partial differential equations are very difficult to obtain because these equations are nonlinear. Therefore only comparatively simple problems have been analytically solved with stochastic recharge and discharge processes (Sagar and Kisiel, 1972; Bakr, 1976; Sagar, 1978a). Intermittent recharge occurs when a drainage system provides relief from precipitation with an irregular pattern and from periodic irrigation. Superposition can be used to study the intermittent recharge problems (Werner, 1957; Maasland, 1959; van Schilfgaarde, 1965a). van Schilfgaarde (1965a) reduced the continuous recharge to a bar histogram with a uniform time increment. Within each time increment the recharge was assumed to be a constant. Choosing a convenient interval, such as a day, climatological data (precipitation, temperature, wind speed, etc.) can be used first to calculate the net recharge and then to calculate the water table positions. In the calculation of net recharge, no surface runoff and deep seepage are taken into account. From the water table positions at the end of each of an arbitrarily long sequence of time increments, which may be as long as a number of years, a frequency distribution of water table heights can be developed, which can be used, in principle, in the selection of a drainage system at an acceptable level of risk. Similar techniques have been used by Chieng (1975) and Skaggs (1978). This kind of method requires a large amount of historical meteorological data and may introduce new uncertainty when historical data is used to predict future precipitation and evapotranspiration. Sagar and Preller (1980) used a different approach to study the case for drainage under uncertain weather conditions. By linearizing the Boussinesq equation and assuming Literature Review 28 different forms of stochastic functions for recharge (precipitation) and discharge (evapotranspiration), they solved the partial differential equation for the first two moments of the hydraulic head (mean and variance). The variance or standard deviation could be used as a indicator or measurement of the uncertainty about the water table position and, therefore, the performance of the drainage system. This approach also requires large amounts of data and its accuracy greatly depends on the selection of the stochastic functional forms of precipitation and evapotranspiration. By applying the Bayesian statistical analysis technique, Musy and Duckstein (1976) and Fogel et al (1979) took into account the crop response to high water table level in designing subsurface drainage systems under uncertainties. Musy and Duckstein's approach requires a probability density function of extreme precipitation and the crop loss model corresponding to the extremes. They limited the problem to situations where only one submergence event can occur in one year and such an event always occurs at the critical growth stage. "Critical" means that the crop is most likely to be damaged by high water table position. Fogel et al (1979), on the other hand, solved a similar problem without the above limitations, but still.required the probability density function of the extreme precipitation and the crop response model. The disadvantage of the Bayesian technique is due to the difficulty of obtaining the probability density function of the extreme precipitation and the crop response model, which are dependent on the local weather conditions and specific crop type. The above research only incorporate the uncertainty of climatic conditions. However, other parameters, such as soil properties, are also subject to great variability and uncertainty. Bakr (1979) and Freeze (1975) studied groundwater flow problems under climate and soil parameter uncertainty by employing analytical or numerical methods. Prasher (1982) made an effort to incorporate both climate and soil parameter uncertainty into drainage/subirrigation design. Literature Review 29 Methods used in the parameter uncertainty analysis of engineering problems can be categorized into two groups: deterministic methods and probabilistic methods (Uhl &; Sullivan, 1982). Deterministic methods assign test values to a variable within some feasible range, but without regard for the likelihood of the respective range. Probabilistic methods incorporate probabilistic information about each of the variables, which yields a probability distribution for the dependent variable. Deterministic analysis is often used to identify dominant variables so that they may be subjected to additional (probabilistic) analysis. One of the deterministic methods commonly used is sensitivity analysis. Sensitivity analysis involves testing a noncorrelated variable, or a group of correlated variables, at a number of values over a reasonably expected range while fixing other variables at their "best-guess" estimates to observe their effects on the dependent variable. Sensitivity analysis is often undertaken initially to identify dominant variables. These variables can be disaggregated and a sensitivity analysis performed on the effects of the newly introduced variables. As the analysis progresses, this method identifies sensitive variable combinations and gives quantified feedback on assumptions in the appraisal. The primary advantages of sensitivity analysis are its flexibility, versatility and simplicity. It can almost test any number of variables and in almost any desired combination subject to various correlation restrictions. It does not require extensive quantitative background information. The major disadvantages of sensitivity analysis are the problem of correlations and the proliferation of data for all but relatively simple applications. Failure to accommodate correlations can invalidate results. The output of a sensitivity analysis for a large or detailed appraisal can be voluminous. The interpretation of the output can become a time-consuming and costly task. Another disadvantage of sensitivity analysis is that it does not reflect the probability of obtaining any particular value for the dependent variable. In spite of these disadvantages, sensitivity analysis is a powerful tool Literature Review 30 to determine the behavior of a function or a model, and to identify the dominant and sensitive variables. It also may serve as a precursor to probabilistic analysis for obtaining likelihood estimates. The probabilistic methods used in the uncertainty analysis in groundwater flow problems can be generally grouped as those which require full probability distributions of all stochastic inputs and those which require only first and second moments of the stochastic parameters. A comparison and review on these two groups of methods was made by Dettinger and Wilson (1981). A summary is given below. Full distribution methods are more accurate, provided that the probability distributions of independent variables are correctly established. The derived distribution method and Monte Carlo simulation method are two of the most commonly used full distribution methods. The derived distribution technique analytically derives the probability distribution of the dependent variable which is functionally related to the random independent variables (Benjamin & Cornell, 1970). However, this analysis is very complicated and time consuming for all but simple applications, unless the problem is simplified by assuming that all variables and results are normally distributed. In such case the reliability of these techniques is dependent on the validity of these assumptions. If the function is approximately linear the results obtained from these assumptions are generally acceptable, otherwise another method (e.g. Monte Carlo simulation) should be used (Uhl and Sullivan, 1982). Sagar and Kisiel (1970) applied the full distribution method to examine the parameter uncertainty for aquifer pump tests. The Monte Carlo method, simulates the groundwater flow system with numerous replications. The simulation involves a sampling technique, usually carried out on a computer, that randomly and iteratively selects sample values for each dominant variable from its respective probability distribution to obtain sample values for the dependent variable. A large number of trials is needed for the resultant frequency distribution for Literature Review 31 the dependent variable to converge to the analytical solution. Monte Carlo simulation is a powerful method to analyse uncertainty and can be applied to most uncertainty analysis problems, such as nonlinear and relatively complex functions. However, when many variables are involved and the function is not simple and well-behaved, the Monte Carlo simulation technique can be very complicated and expensive (TJhl and Sullivan, 1982). In cases where the variables are correlated, identifying and accurately representing these interrelationships can be a complex task. Also the results obtained from Monte Carlo method are not readily transferable to new situations. The Monte Carlo method has been employed to study the uncertainty of transport processes through porous media due to the variability of soil properties (Warren and Price 1961; Freeze 1975; Smith and Freeze 1979). Massman (1987) applied the Monte Carlo method to study uncertainty in the travel time of contaminants in the groundwater flow region. The major disadvantage of full distribution methods is the fact that it is usually very difficult to obtain full probability distributions for all stochastic inputs. In addition, full distribution methods often require advanced mathematic and statistical background and intensive analyses and computations. The results obtained from the full distribution methods are highly dependent on the types of distributions selected to represent each random variable, and if correlationships among variables exist, full distribution methods can be too complicated to be practically applicable. The first and second moment method uses the mean (first moment) and the variance and/or covariance (second moment) to represent the probability properties of a random variable. The mean is a parameter representing the central tendency of a random variable, and variance is a measure of scatter from the mean. For this representation to be accurate, or in other word, the properties of a random variable can be completely characterized by its mean and variance, the distribution must be normal. In this case it Literature Review 32 has zero skewness and other higher moments can be calculated from the variance (Benjamin &; Cornell, 1970). If the function is linear or approximately linear, the dependent variable can also be approximately defined by its mean and variance. If the variables are not normally distributed or the function is highly non-linear or the coefficients of variation of random variables are not small enough, the results becomes unreliable. However, these disadvantages can usually be justified by noting that: (1) the data and physical arguments are often insufficient to establish the full probability law of a variable; (2) most engineering analyses include an important component of real, but difficult to measure, professional uncertainty (due, for example, to imperfect physical theories and to engineering approximations); and (3) the final output, namely the decision or design parameter, is often not sensitive to moments higher than the mean and variance (Cornell, 1972). In addition, first and second moment analysis employs only the same, familiar tools and procedures of algebra and calculus that are commonly used in the more deterministic analysis of the same problem. This implies that both students and practicing engineers can learn the technique and apply it more readily. Wider engineering application of stochastic methods by non-specialists will benefit everyone concerned. Also it is far better to approximately model the whole problem than to exactly model only a portion of it (Cornell, 1972). Cornell (1972) proposed the use of first and second moment analysis based on the Taylor series expansion in the uncertainty analysis of water resource system. He in- dicated that first and second moment analysis can be used for either unconditional or conditional analyses, and for either parameter or model uncertainty. Sagar (1978b) used this approach to analyze one-dimensional flow in confined aquifers by the finite element method. Prasher et al (1982) suggested that first and second moment analysis can be used in uncertainty analysis in the design of a drainage/subirrigation system. The author has employed this approach in the analysis of the uncertainty involved in steady and Literature Review 33 transient state subsurface drainage design. Simple, ready-to-use formulae and graphs have been prepared for convenient use of practical application. Chapter 3 SENSITIVITY A N A L Y S I S O F D R A I N A G E D E S I G N IN H O M O G E N E O U S SOIL In agricultural drainage design, input information of system parameters are estimated from historic climatic data, soil geohydrologic properties and crop characteristics. The estimates are more or less subjective and involve a considerable degree of uncertainty due to the complexity of the nature. It is therefore desirable to examine the relative sensitivity of system performance to each parameter, so that certain precautions can be taken in the estimation of parameter values to optimize the design of a drainage system. This chapter will carry out the sensitivity analysis of drainage design in homogeneous soils. Multi-layered soils will be the subject of the next chapter. 3.1 Steady-state: Hooghoudt's Equation. For the steady-state condition, Hooghoudt's equation (1940) has been most popularly used for drainage design (Luthin, 1978; Chieng, 1981). Hooghoudt's equation for a homogeneous soil can be described as: 4K S = ^(2d h 2 e + h) 2 where S is the drain spacing (m); R is the drainage coefficient (m/day); K is the saturated hydraulic conductivity (m/day); 34 (3.1) Sensitivity Analysis of Drainage Design in Homogeneous Soil 35 h is the water table height midway between the drains; d is the equivalent depth of impermeable layer below the drains. e d can be evaluated from tables prepared by Hooghoudt (1940) or graphs (Bouwer &; e van Schilfgaarde, 1963), or from the equations below: For 0 < f < 0.3 d e = l + f(!ln<!-a) O where \ 7T T ( 3 " 2 ) I a = 3.55 - 1.6| + 2 ( f ) 2 and for f > 0.3 * J(lnf-1.15) = ( 3 3 ) where r is the radius of the drain pipes, S is the spacing, and d is the actual depth of the impermeable layer below drains. Although d is a function of spacing S , actual depth of impermeable layer d, and drain e radius r, it is usually considered as an independent variable when sensitivity analysis is carried out. This can be justified by the fact that results show clearly the sensitivity of drain spacing to other parameters, and that the effect of considering d as an independent e variable on the sensitivity of spacing is relatively small. For detailed discussion, see appendix A . 3.1.1 Sensitivity of Spacing to Drainage Coefficient R The drainage coefficient R is the amount of water that the drainage system can remove from the soil per unit time. R is usually estimated from historic rainfall data, evaporation, consumptive water use of crops, surface runoff, deep seepage etc. It is obvious that such an estimation can not accurately represent the actual drainage requirement, and different designers make different estimations. Sensitivity Analysis of Drainage Design in Homogeneous Soil 36 From basic calculus, we know that: AS\ « dS = R ^zzAR dR 1 for a small change of R, while all other variables are held constant. Also 8S_ _ %(2d h + h ) 2 e 9R~ 2y/%(2d.h + h*) Therefore 2^(2^+ Dividing by S on both sides and rearranging gives: AS S where ^ 1AR . . (3.4) .2 represents the percentage change of drainage coefficient, and ^ is the cor- responding change of the drain spacing. The sensitivity of spacing to the drainage coefficient is approximately a constant. 3.1.2 Sensitivity of Spacing to Hydraulic Conductivity Hydraulic conductivity K is one of the main soil properties directly related to drainage design. Methods of measuring hydraulic conductivity were reviewed by Bouwer & Jackson (1974). However it is impossible to select a measuring method that can accurately evaluate the if-value because of the soil variability. That means if-values at different locations may not be the same, and even at the same location, the if-value may also change from time to time because of compaction, changes in moisture content and soil structure. Therefore, it is important to know the sensitivity of drain spacing with respect to the K- value, before choosing a certain value to be used in the design. Sensitivity Analysis of Drainage Design in Homogeneous Soil 37 Again we have: AS\ * dS = K ^-AK oK in which 8S _ dK jt(2d h + h ) 2 e ^(2d h 2 + h) 2 e The percentage change in spacing due to a small change in if-value is: A£| ±(2d h +h) S \K~ 2 /f{2d h s AK 2 e +h) S 2 e Rearranging the above equation yields: AS_ 1AK , s S where 3.1.3 is the percentage change in hydraulic conductivity. Sensitivity of Spacing to Water Table Height Midway Between Drains Water table height midway between drains is a very important parameter in agriculture drainage design. Since the depth, D, of drains below the field surface is usually determined by the soil profile, excavating equipment, or water table position, the value of h to be used in Eq. 3.1 is D — d i , where d i m n m water table midway between the drains. d i m is the minimum permissible depth of the n is dependent on the tolerance of crops to n waterlogging, and on trafficability for farm machines. Similar to the previous analysis, the change of spacing due to change in h— value is: AS\ h = ~ dS —Ah Oh where ds_ f{d e 9h ~ 2^(2d h e + h) + h) 2 38 Sensitivity Analysis of Drainage Design in Homogeneous Soil 1.0 10.0 djh Figure 3.2: Sensitivity of spacing to midway water table height in relation with ^ and the percentage change of spacing will be: _ d + h Ah AS 2d -f h h S h e e or in the dimensionless form: AS (3.6) 2% + l h where is the percentage change of h. Figure 3.2 shows that the sensitivity of spacing to water table height midway between drains decreases as ^ increases•, the minimum ratio of change of spacing and change of water table height being 0.5. ^ represents the relative importance of water flows in the soil region below and above drains. In other words, when the thickness of soil below the drains increases, the water flow in the upper layer becomes less important. 3.1.4 Sensitivity of Spacing to the Equivalent Depth d e Equivalent depth is the thickness of an imaginary layer of soil in which the water flow is equal to the water flow in the actual soil profile between drains and the impermeable Sensitivity Analysis of Drainage Design in Homogeneous Soil 39 layer. Within the equivalent layer, the water flow towards the drains is horizontal, while the actual water flow includes both horizontal and radial components. In a homogeneous soil, d is related to the spacing, the depth of the impermeable layer below the drains e and the radius of the drain pipe. In this section d is treated as a "real" parameter, and e in next section its relationship with impermeable layer depth and drain radius will be discussed. Following the same procedure as in the previous analysis, we have: AS\ de OS « dS = -KrAd e where as _ dd e ^ 2^(2d h e + h) 2 The percentage change of spacing due to change in equivalent depth is: AS ' S Ad 2d e e +h d e or in dimensionless form: AS k Ad, -T^—^r 2f + 1 d (3-7) 1 v ; e where is the percentage change in d . e Figure 3.3 shows that the sensitivity of spacing to equivalent depth increases as y increases. The ratio of percentage change in spacing to percentage change in d , ranges e from 0 to 0.5. A larger d value implies a greater importance of the equivalent layer, and e thus increases the sensitivity of spacing. Sensitivity Analysis of Drainage Design in Homogeneous Soil 40 0.3 1 ti 0.3- 0.1- 0.0 10.0 Figure 3.3: Sensitivity of spacing to equivalent depth in relation with ^ . 3.1.5 Discussion on Relationship between Equivalent Depth d and Impere meable Layer depth below drains d and Drains Radius r Equivalent depth d can be calculated from• equation 3.2 and 3.3: e ForO < | < 0.3: d, = where d l + |(jlnj-a) a = 3.55 - 1.6| + 2 ( f ) 2 and for f > 0.3: d„ = ' | ( l n | - 1-15) The partial derivatives of d with respect to d and r are, respectively: e For 0 < f < 0.3: dd ^ l - f f - 1 . 6 ( f ) e ad 2 + 4(f) [l + i ( | l n j - a ) p dd 5r e 7r Sr [l + f ( l -a)] 3 Sensitivity Analysis of Drainage Design in Homogeneous Soil 41 and for f > 0.3: dd ~dd e dd dr e f(ln?-1.15)> And for a small change in d and r, the change in d will be: e Ad Ad e (3.8) d e Ad„ Ar (3.9) where 1 / -^: /. ^ : y 1 6 ) + 2 ) 3 for.0<|<0.3 (3.10) for | > 0.3 a d for 0 < | < 0.3 h = for 4 > 0.3 (3-11) Equivalent depth, d , becomes less sensitive to d as j or d increases. When | > 0.3, e d becomes unaffected by the change of d-value. Equivalent depth becomes more sensitive e when ^ decreases or r increases, or when ^ increases. That means that in a deep soil profile, the radius of the drains is a more important factor than in a shallow soil profile. 3.1.6 Sensitivity of Spacing to the Depth of Impermeable Layer below Drains From the discussions in section 3.1.4 and 3.1.5, we can now carry out a sensitivity analysis of spacing to the depth of the impermeable layer below the drains. From section 3.1.4, we have: AS_ _ S d e w h e r e /i = Ad d e e Sensitivity Analysis of Drainage Design in Homogeneous Soil 42 And from section 3.1.5, we have: Ad ~~d Ad,_ where / can be evaluated from equation 3.10. 2 Therefore, the percentage change in spacing due to change of depth of the impermeable layer can be approximated as: AS_ — fi ~S~ d 3.1.7 /2 Ad —r d (3.12) Sensitivity of Spacing to Radius of Drains From section 3.1.4, we have: AS S d a e e And from section 3.1.5, we have: Ad, where / 3 Ar — 73 — r can be evaluated from equation 3.11. Therefore, the percentage change in spacing due to change in radius can be evaluated by: A .<? AS S 3.1.8 A t , , Ar = fifs — V (3.13) Estimation of Spacing Change due to Changes of A l l or a Portion of the Parameters The spacing change, AS, is approximately: - A c JO D S a d AS « dS = ^AR OR dS . jdS . . 8S + w^AK + — Ah + -^-Ad oK oh od A J e e Sensitivity Analysis of Drainage Design in Homogeneous Soil 43 Dividing by S on both sides: AS OS AR dS AK dS Ah dS Ad S dR S dK 8h S dd S e e S From the analyses in the previous sections, the above equation can be replaced by: AS . ^ = e where R AR —R+ e K AK Ah Ad Ar - — K + e —h + e —d + e —r h d r , (3.14) = —2\ C f l _ ~ d +h 2d +h e e ed = / i h . e = /x / P f — 3 e d J2 and f are evaluated from equation 3.10 and 3.11. 3 Let e be a five dimensional vector and equal to: = (e , SK, e , e , e ) e r h d T and let: - _ AR AK Ah Ad Ar The equation 3.14 can be writen as: AS -> ~Y=e»AX 3.1.9 (3.15) Conclusions for Sensitivity Analysis of Steady-state Drainage Design in Homogeneous soil From the analyses in the previous sections, it can be concluded that: Sensitivity Analysis of Drainage Design in Homogeneous Soil 44 1. Spacing is most sensitive to the water table height midway between drains. Therefore special care should be taken in selecting the design water table height in order that it can meet the crop and leaching requirements, while remaining as high as possible so that the spacing can be greater. 2. Spacing is least sensitive to the depth of the impervious layer below the drains and to the radius of the drains. Thus, it is not necessary to determine the depth to the impermeable layer very accurately, and when choosing drain radius, drain flow capacity should be the major concern (i.e. a drain size should be selected for the necessary flow capacity). A larger diameter for the drains will help little to increase the drain spacing. 3. Spacing is also very sensitive to the drainage coefficient and hydraulic conductivity. It is recommended that more research be done on the determination of drainage coefficients and measurement of hydraulic conductivity so that input information can be more accurate. 3.2 Transient State: van Schilfgaarde Solution of Boussinesq Equation Based on D-F assumptions, the Boussinesq equation is commonly used in nonsteady subsurface drainage design (Luthin, 1978; Chieng, 1981). It can be expressed as: where h is the height of the water table above the impermeable layer, x is the horizontal distance, t is time, / is drainable porosity, K is saturated hydraulic conductivity, and R is the rate at which water enters or leaves the saturated region (see Figure 3.4. Eq. 3.16 is nonlinear, and is frequently linearized (Polubarinova- Kochina, 1962) to Sensithity Analysis of Drainage Design in Homogeneous Soil 45 Impervious layer Figure 3.4: Symbols used to describe drainage to a ditch or a tube drain. . obtain its solution. Glover (Dumm, 1954) solved the Boussinesq equation by assuming an initially flat water table: h-d 4m 7T 0 ~ 1 . £ TITTX . —n ir KDt 2 -«»(-5-)«P( f S 2 2 ) (3.17) n=l,3,5,.. For all but the smallest time period, only the first term is considered: ir KDt (3.18) S / l n ^ ir m i where D is the average depth of the flow region, D = d -f m / 2 . m is the initial water 2 2 J e 0 0 table height; m is the water table height after fall. AD other parameters are as previously 1 defined or in figure 3.4. Equation 3.18 is often referred to as the Glover equation (Luthin, 1978). Tapp and Moody (Dumm, 1964) observed that the initial water table shape encountered in the field was more often parabolic than flat. They modified the initial conditions used by Glover accordingly, and obtained a solution that had only minor differences with Sensitivity Analysis of Drainage Design in Homogeneous Soil 46 the Glover equation: „. S = iv KDt 2 (3.19) / I n 7r m\ van Schilfgaarde (1963,1964), on the other hand, solved Eq. 3.16 without linearization: f 1 3 7 m 3 7 7 7 1 0 t = £ - M t 9Kd m^d m v \ + mo) i + J m i ' } (3.20) ' K If the equivalent depth concept is used, Eq. 3.20 will provide a practical solution to a complicated problem in simple form. Eq. 3.20 implies an initial elliptic water table shape. If the observation by Tapp and Moody (Dumm, 1964) that a flattened parabola of the fourth degree is often encountered is correct, the Glover equation would tend to underestimate the required drain spacing, and the van Schilfgaarde equation would tend to overestimate it. Putting the equation of the initial water table shape aside, Eq. 3.20 can be expected to provide more reliable results than Eq. 3.18 and Eq. 3.19 because it is based on Eq. 3.16 without linearization, and it is not restricted to relatively small increments of time and water table. However, it should be noted that, strictly speaking, Eq. 3.20 does not permit superposition because it is not linear (van Schilfgaarde; 1974). For the analysis in this section, the van Schilfgaarde equation is used, mainly because it is simple in form and is not restricted to small increments of time and water table. The equivalent depth, d , is used to substitute for the impermeable layer depth to provide e more realistic results. Eq. 3.20 can be rewriten as: S =. . °\ 2 (3.21) 9K d •> mi(2d +n»o) e with d substituting for d in Eq. 3.20. d can be evaluated from Eq. 3.2 and 3.3. As in e e Sensitivity Analysis of Drainage Design in Homogeneous Soil 47 section 3.1, d is considered as functionally related only to the impermeable layer depth e and to drain radius, to simplify the analysis and results. 3.2.1 Sensitivity of Spacing to Drainable Porosity To derive Eq. 3.16, it is assumed that drainable porosity is a constant. Drainable porosity is defined as the volume of water per unit area drained per unit change in water table height. However this quantity is far from constant, depending on the soil uniformity, on the proximity of the water table to the soil surface, and on the rate and direction of water table movement. The derivative of spacing with respect to drainable porosity is: OS df 9Kd t 1 e \ •\ ln m o ^+ m i 1 r m.\{2d -\-m.o) e This equation can be simplified by substituting: 9Kd t e and becomes: &S__15 df~~2f Therefore, the percentage change in spacing due to a small change in drainable porosity is: AS S 2 / A/ where —jis the percentage change of drainable porosity. 3.2.2 Sensitivity of Spacing to Hydraulic conductivity The derivative of spacing with respect to hydraulic conductivity is: (3.22) Sensitivity Analysis of Drainage Design in Homogeneous Soil 9d t 0S_ 1 1 e OK m (2d +mi 0 )| 2 e 48 1S_ 2K Thus the percentage change in spacing is: AS S 1 AK (3.23) K This result is the same as in steady state. 3.2.3 Sensitivity of Spacing to D r a w - d o w n T i m e from W a t e r Table m to mi 0 If the water table height is m immediately after recharge ends, it will fall to m-i after a 0 time lapse t. t is the time interval during which the water table falls from mo to mi, and is determined by the tolerance of the crop to waterlogging. A shorter t requires a higher capacity to the drainage system and costs more for the farmer. On the other hand, if t is too long, crop growth may be adversely affected due to excess water in the rootzone. This draw-down time is directly related to the cost-effectiveness of the drainage system. The partial derivative with respect to draw-down time t is: dS 1 1 9Kd„ .15 And the percentage change in spacing due to change in selected time t is: AS 5 1 Ait 2 (3.24) 7 Eq. 3.23 and 3.24 are exactly the same in form. This implies that the selection of draw-down time is as important as the determination of hydraulic conductivity to the adequate and economic design of a drainage system. Sensitivity Analysis of Drainage Design in Homogeneous Soil 3.2.4 49 Sensitivity of Spacing to the Initial Water Table Height after Recharge Ends If the recharge is natural precipitation, the initial water table height, m , is mainly 0 determined by the rainfall intensity and duration with a selected frequency of occurence. In irrigated land, m is determined by the application rate and time of irrigation. 0 The derivative of spacing with respect to initial water table height, m , is: 0 dS _ _S ^ ~ ~ 2 m 2 4 o ( 2 i e + m ) l n [ ^ g ^ ] o And the percentage change in spacing is: AS S Am-o m (2d +m )\n{Z%i Z\] + 0 e Q m o + or, in dimensionless form: AS S Am de mo m o 0 (2-i.+ 1 ) ^ ( ^ 1 1 ) v m 0 . ' ' \ m ' ° (3.25) Figure 3.5 shows that spacing becomes more sensitive to the initial water table height as ^m o increases or ra — m decreases. This indicates that the less the water table draw0 x down, the more sensitive the spacing to the initial water table. ^ little influence on the sensitivity. When ^ seems to have very is greater than 0.5, spacing becomes very sensitive to the selection of initial water table height. 3.2.5 Sensitivity of Spacing to the Water Table Height after Recession The water table height after recession within a given time period should be such that if the water table is lower than this height, recharge is needed to provide adequate water to meet crop requirement. A higher water table may induce waterlogging and a lower 50 Sensithity Analysis of Drainage Design in Homogeneous Soil O A to.o ^ 1 d8/fri0=0.2 d»/m0"2.0 V 4*M0=3.0 . m /m 1 0 • Figure 3.5: Sensitivity of spacing to the initial water table height m 0 water table may impede crop growth. The water table position after recession is directly related to drain spacing. •The derivative of spacing with respect to final water table height, m dS _ dm --2 1 S m i . (2d 2d . e + 1 } is: e m )ln[^^ } 1 ) Therefore, the percentage change in spacing is: Ami AS or in dimensionless form: AS ~S~ Amx mi (3.26) «»0 Figure 3.6 shows that spacing becomes more sensitive as ^ httle influence on the sensitivity. increases. Again ^ has 51 Sensitivity Analysis of Drainage Design in Homogeneous Soil d«/m0"0.2 <U/mO«C.O V <J«/mO«»9.0 -2.0- S -4.0<> / -6.0- -8.0- -10.0 0.0 0.2 0.4 0.6 0.8 1.0 Figure 3.6: Sensitivity of sapcing to water table height after recession. 3.2.6 Sensitivity of Sapcing to the Equivalent Depth We will follow the same procedure as in steady state analysis. First, the sensitivity of spacing to equivalent depth will be" discussed and then the results of section 3.1.5 will be used to establish the relationships between spacing and impermeable layer depth and drain radius. The derivative of spacing'with respect to equivalent depth is: i '8S dd •J 1 rmo(2ct + m ]] \-mi(2dc+mn )J e 2(Tno-mi)d x d e 1 e n e (2d +mi)(2d +m ) c e 0 [">o(2rf.+m )j 1 Lmi(2d .+mo) J 6 The percentage change in spacing is: 1 AS 2(TT1(> — m i )d c (2d,,+mo)(2ci +m ) e 1 Ad e S or, in dirhensionless form: AS S * 2(Ahi- 2 \ m i mo ,Ad P (2^ + l)(2^ + l)ln(2^ + l ) - l n ( 2 ^ + l) d J e (3.27) 52 Sensitivity Analysis of Drainage Design in Homogeneous Soil V d«^n0=0.2 tl«An 0*2.0 0.0 0.2 0.6 0.4 •mjm doAnO=9.0 l 0.8 1.0 a Figure 3.7: Sensitivity of spacing to equivalent depth Figure 3.7 shows that the sensitivity of spacing to equivalent depth does not vary much, and is less than 0.5. Neither ' 3.2.7 or 7 2 1 has significant influence on the sensitivity. mo mo ° Sensitivity of Spacing to Actual Depth of Impermeable Layer below Drains and Drain Radius From the previous section, we have: AS S de fi ,Ad e w here x r-. 2L ^ - ^ 1 mo / _ (2-*+i i1\(o ) ( 2d,^ +i 1) (2^+i 1) - li„n (/"o 2 £d- + 1) fn d 1 \ In i_ fn d i \ \ mi mo ' mi / \ mo ' 2 '•mi c And from section 3.1.5 we have: : v € c J Sensitivity Analysis of Drainage Design in Homogeneous Soil Ad e where f 2 and / 3 53 , Ar = /a— r are evaluated from Eq. 3.10 and 3.11. Therefore, the percentage changes in spacing due to the changes of the impermeable layer depth and the drain radius are: AS ~S~ = fi (3-28) , Ar = //a— r (3-29) and AS S respectively. The total percentage change in spacing due to changes in all or some of the parameters can be estimated by similar methods to those in section 3.1.8. 3.3 C o m p a r i s o n of Sensitivities between Steady and Transient Drainage D e sign i n Homogeneous Soils In this part of the chapter, a comparison of the sensitivities of spacing between steady and transient drainage design, based on the analyses of previous sections of this chapter, is made for individual parameters. The parameter values used in the comparison are as follows: For non-steady state: Initial water table height at mid-spacing: mo = 0.8 m; Water table height at mid-spacing after recession: m\ — 0.4 m; Time required for water table to recede from m to m : t = 4 days; 0 1 Sensitivity Analysis of Drainage Design in Homogeneous Soil 54 Saturated hydraulic conductivity: i f = 0.5 m/day; Drainable porosity: / = 6%; Depth to impermeable layer below drains: d = 5.0 m; Drain radius: r = 0.1 m. In order to make the comparison, the input values for steady state drainage design should be equivalent to the drainage requirement of a falling water table, which can be approximately achieved by selecting R and h values accordingly. R and h values can be chosen by: (m - )/ fi= o mi (0.8-0.4)x0.06 t 4 , mo + m.! 0.8 + 0.4 h = -— = = 0.6 m 2 2 = n ' e i f , d and r are the same as for the transient state. The results are shown in table 3.1 (columns 2 and 3). All parameters are assumed to change by 10%. The results shown in column 4 of Table 3.1 were obtined by maintaining the same falling rate of water table while adjusting m-i from 0.4 m to 0.7 m and required falling time t from 4 days to 1 days. From the results in Table 3.1, it can be seen that: 1. Transient drainage design (falling water table criterion) can be reduced to the simple steady state drainage design, by appropriately adjusting the values of R and h. Certain combinations of R-h values may result in close or the same spacing as that obtained with the transient state drainage design method. 2. The sensitivities of spacing to hydraulic conductivity, impermeable layer depth, and drain radius in steady drainage design are very close to those in transient drainage design. Sensitivity Analysis of Drainage Design in Homogeneous Soil 55 Table 3.1: Comparison of sensitivities between steady and transient state drainage design. Parameter 1 d (m) 5(m) R K h Steady state 2 2.44 33.59 -5.00 5.00 5.46 d r 1.27 0.80 e Transient state 4 3 2.58 2.81 35.24 42.32 m TTlj 0 f t 5.00 -6.94 -7.44 1.44 0.84 -5.00 5.00 5.00 -37.15 -37.73 1.70 0.75 -5.00 5.00 3. The spacing is least sensitive to impermeable layer depth and drain radius in both steady and non-steady drainage design. Therefore, in both cases it is not necessary to measure the depth to an impermeable layer very accurately, considering the additional cost of more accurate measurement. 4. The spacing is very sensitive to the initial and final water table heights, when they are close to each other. Therefore, it is recommended that for cases where the water table fall or fluctuation is small, steady state drainage theory, such as Hooghoudt's equation, be used instead of transient drainage theory, such as the van Schilfgaarde equation. In these situations, spacing obtained from steady state equations tends to be more stable and therefore reliable than that from transient state spacing equations. 5. If transient drainage theory is used, m 0 and m-i should be carefully selected, in order to avoid the over or under estimate of spacing for design. Chapter 4 SENSITIVITY A N A L Y S I S O F S T E A D Y S T A T E D R A I N A G E D E S I G N IN M U L T I - L A Y E R E D SOIL In chapter 3, the sensitivity of drainage design in homogeneous soil is discussed. However, in practice, multi-layered soil is a more common phenomenon. For example, a field under tillage may gradually form a relatively compacted layer below certain depth. Two-layered or three-layered soils are commonly encountered in real fields. In section 4.1 sensitivity analysis of drainage design in two- layered soil with drains at the interface of these layers will be carried out. In section 4.2, three-layered soil with two layers below the drains, and two-layered soil with drains in the upper layer will be analyzed. Hooghoudt's equation (1940), the generalized Hooghoudt-Ernst equation (van Beers, 1976), and the Kirkham equation (Tokoz & Kirkham, 1971a) will be used. A combined equation from the Hooghoudt and Ernst equations for layered soils will be developed. Comparison will be made accordingly. 4.1 Sensitivity Analysis of Steady State Drainage Design in Two-layered Soil Notations used in this section are: R is the drainage coefficient; K\ is the hydraulic conductivity of the soil above drains; K 2 is the hydraulic conductivity of the soil below the drains; h is the water table height midway between the drains; S is the drain spacing; 2r is the diameter of the drains; and an impermeable layer is at a depth, D , below the 2 drains. The concept of equivalent depth will again be used in the Hooghoudt equation and in the combined equation. 56 Sensitivity Analysis of Drainage Design in Multi-layered Soil 4.1.1 57 Drainage Design Using Hooghoudt's Equation The Hooghoudt equation for a two-layered soil is (van Beers, 1976): Ah S = —{2K d K + K h) 2 2 e (4.30) 1 where d can be calcualted from Eq. 3.2 and 3.3. d is again considered as independent e e of spacing to simplify the analysis and results (See Appendix A ) . Sensitivity of Spacing to Drainage Coefficient From calculus, we know that AS\ R « dS = ^zrAR OR for small AR, where dS_ _ 1 -^(2K d 9R~ 2 /%{2K d +K h) 2 y e 3 e + K h) 1 _ _ l 5 " 1 Therefore A9I - R A R The percentage change in spacing is: AS S where ^ R 1AR 2 R (4.31) is the percentage change in the drainage coefficient. The above result is the same as that in homogeneous soil. Sensitivity of Spacing to Hydraulic Conductivity of the Upper Soil Layer Similarly, we have 58 Sensitivity Analysis of Drainage Design in Multi-layered Soil K77Kt=0.2 K2/K1-1.0 K2A1«9.0 V 10.0 Figure 4.8: Sensitivity of spacing to hydraulic conductivity of upper layer soil. AS\ Kl 8S ~ dS = jr^AK, ). dK ±h 2 = —^ •2^(2K d 2 AK X + K,h) e Dividing both sides by S gives: AS S - 1 Kh AKx ~ 2 2K d + K,h K. x K l 2 e x or, in dimensionless form: AS S where ^ and K l AK 2(1 + 2 ^ ) ^ j X k (4-32) ^ are percentage change in the spacing and in the hydraulic conductivity of the upper layer of soil, respectively. Fig. 4.8 shows that the sensitivity of spacing to K decreases as ^ and ^2 increase. x Sensitivity of Spacing to the Hydraulic Conductivity of the Lower Soil Layer For the hydraulic conductivity of the lower layer of soil, we have Sensitivity Analysis of Drainage Design in Multi-layered Soil Q A 0.5n 59 KZ/K1-0.2 K2/K1-1.0 V K2/K1m9.0 3I» 10.0 Figure 4.9: Sensitivity of spacing to the hydraulic conductivity of the lower soil layer. AS\K 2 ~ dS_ * dS %2d -__ AKn • 2^{2K d + K h) AI<2 d K e = 2 2 e 1 Dividing both sides by S gives the percentage change in spacing: Kd AS 2 K S . e 2K d +K!h 2 2 AK 2 K e 2 or, in dimensionless form: AS S where K l + 2 ^ i 2 (4.33) K 2 is the percentage change in the hydraulic conductivity of the lower soil layer. Fgiure 4.9 shows that the spacing becomes more sensitive to K as | k and ^ increase. 2 Comparison of Sensitivity of Spacing to K~i and K 2 Assuming that and are equal, then if 1 KJi 2 2K d +K h 2 e l > 2"e 2K d 2 e + Kh l Sensitivity Analysis of Drainage Design in Multi-layered Soil 60 spacing is more sensitive to K\ than to K . 2 The above inequality can be simplified to: Kid* K h < x or Kj Kx 1_ 2 \h_ 2d < (4.34) e That is to say, if Equation 4.34 is true, spacing is more sensitive to K\\ otherwise spacing is more sensitive to K . In practice, condition 4.34 is not often satisfied there is 2 a highly pervious soil layer above drain level and a very poorly pervious soil layer below drain level (K 2 <C K\). Thus, the spacing is usually more sensitive to the hydraulic conductivity of the lower layer than to that of the upper layer. In other words, the hydraulic conductivity of the lower soil layer is usually of greater significance to the drainage design; in some extreme cases, the upper layer can be completely neglected without major effects on the spacing, such as in the Kirkham equation (1971a) and the modified Hooghoudt-Ernst equation for layered soils. Interestingly, the sensitivities of spacing to K\ and K add to | , which equals to the 2 sensitivity of spacing to K in homogeneous soils. Sensitivity of Spacing to Water Table Height h For the midspan water table height, we have AS\ h % dS = dS Ah = m Ah Dividing both sides by S gives: AS S Kd 2 h e 2K d 2 e + K^h Ah + Kxh h 61 Sensitivity Analysis of Drainage Design in Multi-layered Soil K2/K1sQ.2 K2/K1- 1.0 V K2/Kt«3.0 10.0 0.0 jd,/ft Figure 4.10: Sensitivity of spacing to water table height. or, in dimensionless form: AS S (4.35) Figure 4.10 shows that the sensitivity of spacing to the water table height, h, decreases as or ^ increase. Sensitivity of Spacing to Equivalent Depth The derivative of spacing with respect to the equivalent depth is: %2K dS 2 dde . 2y/%{2K d 2 e + K,h) The percentage change in spacing is: AS S %2K _ 8S_Ade _ dK ~ dd e S Ad, 2 ~ 2y/%(2K d 2 e + Kih) Rearranging the above equation gives: AS Kd Ad 2K d -f K\h de 2 d c 2 e e e s Sensitivity Analysis of Drainage Design in Multi-layered Soil S 0.5-1 K2/K1-0 K2/K1-0.2 K2/K1-1 0 62 V K2/K1-5.0 T3 * U 0.2H Figure 4.11: Sensitivity of spacing to equivalent depth or, in dimensionless form: K Ki AS 2 S d ' 1 + d h Ad„ t K\ h ° (4.36) e Eq. 4.36 is exactly the same as E q . 4.3.3, which implies that the equivalent depth is as important as the hydraulic conductivity of the lower soil in drainage design. Figure 4.11 shows that spacing becomes more sensitive to the equivalent depth as | p or increase. Sensitivity of Spacing to the Thickness of the Lower Soil Layer and the Drain Radius From the previous section, we have: „Ad AS S h d, d,. w here f" e — Kz is. Ky h TT2fff Sensitivity Analysis of Drainage Design in Multi-layered Soil 63 and from section 3.1.5, we have: AS S D D 2 AS ~S~ with f 2 „ Ar h— T and / evaluated from Eq. 3.10 and 3.11. Here, D is equivalent to d in homoge3 2 neous soil. Therefore, the percentage changes in spacing due to changes of the thickness of the lower soil layer and the drain radius are: AS S AS S D D 2 2 = f"h— r (4.38) respectively. 4.1.2 Drainage Design Using the Generalized Hooghoudt-Ernst Equation A generalized Hooghoudt-Ernst equation was obtained by considering only the horizontal flow above the drains and both horizontal and radial flows below the drains with the radial flow component being from the equation of Ernst (van Beers, 1976). The generalized Hooghoudt-Ernst equation can be written as: S> + l f l , 5 ' In 7T u + R K * D ' » S = "K>y> h TTR l n £l U (4.39) where S =drain spacing; K~i =hydrauhc conductivity of the soil above drain level; K 2 =hydraulic conductivity of the soil below drain level; h =hydraulic head: the height of the water table above drain level midway between Sensitivity Analysis of Drainage Design in Multi-layered Soil 64 drains; Di =h/2, average depth of flow region above drain level or average thickness of the soil layer through which the flow above the drains takes place; D =depth to an impervious layer below drain level; 2 R =drain discharge rate per unit surface area per unit time; u =wetted section of the drain; for pipe drains: u — irr v =drain radius. The partial derivatives of spacing with respect to other independent variables are as follows: an O 8(KiDi+K D )hS 2 . 2 R_ D dR 3S + D,Sln^ 2 dK. 35 + In Bx 2 h ^K K D )h R lDl+ 8D hS R C l O2 _ * u 2 2 2 ~ 0K 2 3 5 2 , 16 D Qg 352 = + 2 x ie D Sln 7T OS _ 1+ 1 64DiD h 8(K D +K h)S 2 8{K,D K D )h R n TT ^ 2 i_ 2 2 D 2 _ _y__ K,D,)h. 1+ R U Z - S ){1 + In ^ ) 2 + ~ 3S + ^D S\n^- 2 2 7T R U Z 8D S _ 2 2 64KiD D h 1 7TT dr D U _ U z dS dL\ n 7Til 7T dh K suhg* i D± + ^D S 2 R U It 2 l 2 K 7C dS l TTfl - *( ^+ * *) m ______ 6AK D D h f ~~ 3 5 2 , 1 6 7T o D * 2 TXRT l n L_ _ U SjK^+K.D^h R The percentage changes in spacing, due to the individual change in each variable while others are held constant, are: Sensitivity Analysis of Drainage Design in Multi-layered Soil 8(K Di+K D )h AS 1 , 64K D D h 2 1 jR 1 2 TTRS "i R ^r 2 8(K D +K D )h R 2 AS S 1 (4.40) TTRS 3 S + ^D S ln 2* R AS S 2 65 X 1 2 U TT AK u N R 2 (4.41) X it Air, tf 3S 2 + ^D S 2 In ^ IL 2 TT AS S 8(^2^2 +A"ife)fe gj^p AS 16 ir + p 1 ^ _ n 2 D c ) h 2 1 1 D 2 u 1 A L . A S i K ^ + K ^ k R l Sx) + (4.43) a h n A L > 2 (4.44) it Ar 2 7Tfl5 l (4.42) ^ K 64JCi£>i£) /i 2 16 D 2 i 2 2 U a 7T + K S l n ^ - ^ ^ 4 ^ ^ "AT 8£> S _ 352 + TTAS u D J 7T AS S l R 1 * 2 2 5 D g ^ ( M j i p _ 5)(i + 3S + ^ D 2 l . 6AD D h ~r R g 2 3 ^ n Dx _ 7T ^ D y + ^ D ^ h U (4.45) r i? The total percentage in spacing due to changes of all or some of the variables is AS „ . AR „ . AK „. A / f „, A/i „. A D „. A r -j- = =ft-g-+ BJr,-^ + tf,,-^ + K j + B i , - ^ + K X where E' , E' , E' , E' , E' R Ki K2 h D2 and 2 Dot product form can be used to write the above equation: f =K .AX where x I rrrl jpl T?I - _ AR AK\ AK ~ ^ ~Kr ^ ~ ~h''rJ ~' , . T . (4-46) can be evaluated from the coefficients of the right hand side in Eq. 4.40 through 4.45. 771/ 2 , , 2 2 771/ 2 771/ T?l\ Ah AD 2 Ar V } Sensitivity Analysis of Drainage Design in Multi-layered Soil 4.1.3 66 Comparison of Spacings and Sensitivities Obtained from the Hooghoudt and the generalized Hooghoudt-Ernst Equations The parameter values used in the calculations are: In case 1: ifi=0.5 m/day; K =0.5 2 m/day; #=0.01 m/day; _D =5.0 m; h=0.6 m; 2 r=0.1 m. In case 2: K_=0.2 m/day; K =2.0 m/day; £=0.012 m/day; £> =5.0 m; h=0.6 m; 2 2 r=0.1 m. In case 3: ^ = 2 . 0 m/day; K =0.2 m/day; #=0.004 m/day; D =3.0 m; h=0.6 m; 2 2 r=0.1 m. In case 4: #1=0.5 m/day; if =0.5 m/day; #=0.01 m/day; £> =6.0 m; h=1.0 m; 2 2 r=0.1 m. In Table 4.2, Row (1) represents the spacing obtained from the Hooghoudt equation and the sensitivities obtained from the Hooghoudt equation with simplified d -S rela- e tionship; Row (2) represents the results obtained from the generalized Hooghoudt-Ernst equation; and Row (3) represents the sensitivities obtained from the hooghoudt equation with unsimphfied d -S e relationship (see Appendix A). A l l changes in parameters are assumed to be 10% of the original values. The results calculated (Table 4.2) show that: 1. The Hooghoudt equation gives very close spacings to those obtained from the generalized Hooghoudt-Ernst equation for various conditions (K K ,K~i <C K ). 2 2 x = K ,K 2 X ^> This is in contradiction to the claim that the Hooghoudt equation can not be used in the situation where K <C K (van Beers, 1976). x 2 2. The sensitivity of spacing in the Hooghoudt equation and that in the generalized Hooghoudt-Ernst equation are not greatly different from each other. The source of difference is in part from the assumption used in the analyses that the equivalent Sensitivity Analysis of Drainage Design in Multi-layered Soil 67 Table 4.2: Comparison of spacing and its sensitivity to individual parameters calculated from the Hooghoudt and the generalized Hooghoudt-Ernst equations. Parameter Spacing (1) (m) (2) R Kx K 2 h D r 2 (1) (2) (3) (1) (2) (3) (1) (2) (3) (1) (2) (3) (1) (2) (3) (1) (2) (3) Case 1 23.96 24.24 -5.00 -6.39 -6.76 0.63 0.82 0.85 4.37 5.85 5.91 5.63 7.48 7.61 0.80 1.19 1.08 0.97 1.32 1.31 Case 2 49.25 49.48 -5.00 -6.20 -6.26 0.05 0.06 0.06 4.95 6.15 6.20 5.05 6.27 6.32 2.17 2.77 2.71 0.77 0.96 0.96 Case 3 34.88 34.92 -5.00 -4.85 -5.35 2.96 3.15 3.16 2.04 2.19 2.18 7.96 8.48 8.51 1.09 1.19 1.16 0.31 0.33 0.33 Case 4 36.08 36.37 -5.00 -6.13 -6.53 0.77 0.98 1.00 4.23 5.48 5.52 5.77 7.44 7.53 1.07 1.49 1.40 0.82 1.08 1.07 Sensitivity Analysis of Drainage Design in Multi-layered Soil 68 depth is only functionally related to the depth to the impervious layer and the drain radius, and is independent of spacing. 4.1.4 Conclusions from the Sensitivity A n a l y s i s of Steady State D r a i n a g e Design in T w o - L a y e r e d Soils From the analyses in the previous sections, it can be concluded that: 1. Drain pacing is very sensitive to the design water table height midway between the drains, the drainage coefficient, and the hydraulic conductivity of the soil below drain level. 2. Spacing is usually less sensitive to the hydraulic conductivity of the soil above the drain level than that below the drain level. The sensitivities are directly related to the ratios of and ^ , both of which represent the relative importance of the upper soil layer and lower soil layer. 3. Spacing is the least sensitive to the depth of the impermeable layer and the drain radius. Therefore, it is not necessary to measure the impermeable layer depth very accurately. 4. The generalized Hooghoudt-Ernst equation and the Hooghoudt equation yield similar spacings and sensitivities. Hooghoudt's equation can be used under most conditions for two-layered soils. 4.2 Sensitivity A n a l y s i s of Steady State Drainage Design i n Three—layered Soil Notations used in this section: K is the hydraulic conductivity of the second soil layer 2 with a thickness of D ; K$ is the hydraulic conductivity of the third soil layer with a 2 Sensitivity Analysis of Drainage Design in Multi-layered Soil 69 thickness of £> . All other parameters are the same as difined in section 4.1. K = K 3 x 2 represents a special case that there are two layers of soil and the drains is in the upper soil layer. The original Hooghoudt equation can not be used in this case (van Beers, 1976). Toksoz & Kirkham (1971a,b) developed theories and nomographs for two- layered soils with the drains in the upper layer, and with some modification, for three-layered soils with the drains at the interface of the first two layers (Wesseling, 1964). But the Kirkham equatins are difficult to use and the nomographs can not be accurately used in some situations. The generalized Hooghoudt-Ernst equation for three-layered soils also requires to solve a three order equation, which is sometimes inconvenient in practice. The trail-anderror method used by the Hooghoudt approach is easy to use at any circumstances and is commonly accepted. In this section, a combined equation using the Hooghoudt approach and the concept of equivalent depth with Ernst's component for the radial flow towards drains will be developed. The sensitivity analysis will be carried out based on the newly developed equation. Spacings and sensitivities obtained from the new equation will be compared with those from the Kirkham equation (Toksoz & Kirkham, 1971a). 4.2.1 Development of A C o m b i n e d F o r m u l a for Spacing C a l c u l a t i o n from the H o o g h o u d t and E r n s t E q u a t i o n s The general principle underlying Ernst's basic equation (1962) is that the flow of groundwater towards parallel drains, and consequently the corresponding available total hydraulic head (h), can be divided into three components: a vertical (v), a horizontal (h), and a radial (r) component as: h = h + hh + h = qL + qSLh + qSL v r v T where q is the flow rate, S is the drain spacing, and L is the resistance. Sensitivity Analysis of Drainage Design in Multi-layered Soil 70 Ernst worked out various resistance terms, and the above equation can be written as: j D S* K 8KL) S h = q-TT- + q-^jr V + q—TVK . aD In 2 u 2 where h, K ,D ,S 2 D v 2 = notation dinned in the previous section; = thickness of the layer over which vertical flow is considered; in most cases this component is small and may be ignored (m); K v = hydraulic conductivity for vertical flow (m/day); KD = the sum of the product of the hydraulic conductivity (K) and thickness (D) of the various layers for the horizontal folw component; a = geometry factor for radial flow (Appendix C); u = wetted section of the drain (m); for pipe drains u = irr. Now let's consider the flow region below the drains. The total flow towards the drains in this region, from the Ernst equation and neglecting the vertical component, can be written as: 8(K D ,+ 2 92 = g 2 K D )h 2 3 3 _i 8(K" £)2+Jir3£)3) o j wK 2 2 n a£> u 2 Borrowing the concept of equivalent depth in the Hooghoudt equation, and assuming an equivalent layer of soil exists such that the horizontal flow within this layer is equal to the total flows of the horizontal and radial components in the whole region below the drain level, with d being the thickness of this layer, then the water flow in the equivalent e layer is equal to: , 8Kd h e where K is the weighted hydraulic conductivity of the equivalent layer, - = KD + KD D + D 2 2 3 2 3 3 Sensitivity Analysis of Drainage Design in Multi-layered Soil Figure 4.12: Geometry of the subsurface drainage in three-layered soils. 71 Sensitivity Analysis of Drainage Design in Multi-layered Soil 72 Based on the assumption of an equivalent layer, we have: <?2 = <? 2 thus %{K D + K D )h 2 g i 2 2 3 SKd h 3 HK2D2+KZD3) TrK2 e a p £2 2 U Rearranging the above equation gives: d e - S + i(D 2+ f^ )ln^ ( 4 3 - 4 ? ) For the flow region above the drain level, only horizontal flow is considered, an approach used by the Hooghoudt equation (Hooghoudt, 1940). This flow is equal to: SKtDth q i = where The total flow towards the drains from both regions should be equal to the drainage coefficient. Therefore, the combined equation becomes: „ %K D h R=—-——-l S x x 2 SKd h — S e 2 or > i D i + with d evaluated from Eq. 4.47. e U_^K^ ] ( , 4 g ) Sensitivity Analysis of Drainage Design in Multi-layered Soil 4.2.2 73 Comparison of Spacing between the New Spacing Equation and the Kirkham Equation Toksoz and Kirkham (1971a) developed an equation for two-layered soils (K x , ,K , M - ^ - l K 2 „ 1 .,• 1 S-{\n—— 7T sin(-j) = e [ 1 + ^ 1 2mirD ,, . 2mivr. Y. - [ - 1 + o o t h — - ) ] [ c o s ( — =i 2 m b 2 , r b = K ): - C O S T I M T ] m Y 2mnD /S 2 " sinh(^^)gcoth(^) + coth(^) ] } ( 4 " 4 9 ) Equation 4.49 is based on the assumption that the loss in hydraulic head in the region lying below the arch-shaped water table and above the level of the drains is negligible compared with the head loss over the remainder of the flow region. In cases where K x ^ K , Wesseling (1964) indicated that Eq. 4.49 still can be used 2 by substituting the factor (K /R 2 — 1) in the left side with (K /R 2 — K /K ). 2 x The parameter values used in the comparison of spacings are as follows (adopted from van Beers, 1976): h=0.6 m; #=0.006 m/day; K = K = 1.2 m/day; r=0.1 m; D = l-6 m. x 2 2 All other parameters are as shown in Table 4.3. In Table 4.3, (1) is the spacings obtained from the new equation by assuming D\—Q, i.e., the head loss above the drains is neglected, which is the case in the Kirkham equation; (2) is the spacings obtained from the new equation considering the head loss of the region above the drains; (3) is the spacings calculated from the Kirkham equation. The spacings calculated from the Kirkham equation is slightly different from those given by Kirkham (1971a), probably because of the different converging criteria used in the calculations. The results in Table 4.3 show that the Kirkham equation is likely to underestimate the required spacing. For the new spacing formula, if the head loss above the drains is neglected (i.e. D = 0), the spacing obtained from the new spacing formula is very close x to that computed by Kirkham's equation. However, if the head loss above the drains is too great to be neglected, the spacing computed from the new equation is always larger Sensitivity Analysis of Drainage Design in Multi-layered Soil D /D D K /K 0.02 0.1 0.2 0.5 1.0 2.0 5.0 10.0 50.0 oo 3 0.25 0.4 2 3 3 2 K 0.024 0.12 0.24 0.60 1.20 . 2.40 6.00 12.0 60.0 oo 3 (1) 36.1 36.4 36.7 37.8 39.4 42.3 49.4 57.7 84.3 124.5 (2) 40.1 40.4 40.7 41.7 43.2 45.9 52.7 60.6 87.6 127.3 74 1.5 2.4 (3) 36.0 36.3 ~36.6 37.7 39.3 42.1 48.9 57.2 84.7 123.8 (1) 36.1 37.9 39.6 44.4 50 4 59.7 75.9 89.5 113.3 124.5 r (2) 40.2 41.8 43.5 47.9 53.6 62.7 78.3 91.7 115.9 127.3 4.0 6.4 (3) 36.3 37.8 39.6 44.1 50.2 59.0 74.4 87.4 111.6 123.8 (1) 35.6 39.8 43.4 51.3 60.7 73.1 91.3 103.3 119.3 124.5 (2) 39.9 43.7 47.1 54.7 63.7 75.8 93.7 105.7 121.6 127.3 (3) 36.7 39.6 42.9 50.9 60.5 72.6 90.2 102.0 118.3 123.8 Table 4.3: Comparison of drain spacings in two-layered soils obtained from the new spacing formula and the Kirkham equation (1971). than the result from Kirkham's equation, as the Kirkham equation neglects the head loss above the drains (i.e. all flow toward the drains is assumed to take place in the region below the drains, which is obviously not true). When the hydraulic conductivity of the soil above the drains is relatively small compared with that of the soil below the drains, the results given by the Kirkham equation are acceptable, otherwise it will yield much smaller spacing. The new spacing can be used in three- or two-layered soils with the simple trial-and-error approach. The calculation is much simpler compared with the Kirkham equation and tends to provide more reasonable results. Sensitivity Analysis of Drainage Design in Multi-layered Soil 4.2.3 75 S e n s i t i v i t y A n a l y s i s of Drainage Design U s i n g the N e w l y Developed Spacing Formula From Eq. 4.48, the spacing, S, is a function of R, K , h, K , D , K , D and d ; and x 2 2 3 3 from Eq. 4.47, d is in turn a function of S, K , D , K , D and r, i.e. e 2 S = 2 3 3 F(R,K ,h,K ,D ,K ,D ,d ) 1 2 2 3 3 e d =f(S,K ,D ,K ,D r) e 2 2 3 3l The partial derivative of spacing with respect to the drainage coefficient is: 8S dR dF dR dF dd„ dS + dd dS dR L P Therefore: 1 9S dR 1 - dF dd dd dS c e Similarly, we have: dS dK x "" 1 dF dKi dF -_ dd dde dS e dF dS 8h 9F 9c_ ~5K ~ 1 -. Vde~~dS~ dS dK 2 dF dde dd dK dF dde -_ dd dS _|_ 9F e ~ 1 2 e dS dD 2 dF 1 ' ~ 1- dF dde dde dD dF dde as dd e dS dK 3 dF 8K 3 1- j 2 dF dd dd dK$ dF dd dde dS e e e e Sensitivity Analysis of Drainage Design in Multi-layered Soil 8F , dF dd T dde 5£>3 dD 1 _ dF dde dde dS dS dD e 3 3 ai 1 > dr dde ' U 76 8d X dS e Various terms in the above equations are: dF_ _ _ 1 S_ dR ~ ~2R dF _2h 2 dKl~RS dF dh _i{K h+ *%X% **) RS dF ±hD d dK ~ RS(D + D ) K Dt 1 = 2 2 dF e 2 3 Ah(K - = K )D d 2 dD 3 RS(D 2 + 2 dF _ D) 3 3 e d~K~ ~ RS(D dF dD 2 dd ~ 2 2 3 3 irK (D 2 2 3 3 _ 8dl(K D 3 (D + D yS e + D )S S+ f(§-l)£> m^ 2 3 + D) 3 2 dD KD) 8K D d \n^ = dK = 3 2 e e e D) + RS(D e 2 2 2 2 dd 2 RS{D + dF_ _ 4h(K D dd 3 K )D d 3 ~ 3 + D) 2 Ah{K - = e 2 AhD d 3 3 2 + 2 irK D {D e 3 2 2 2 KD) 3 3 + D )S 3 dde_ _ ________}_______ dK 3 TK (D 2 2 + D )S 3 dde_ _ S + f ^ m ^ ( l - ^ ) dD 3 D + 2 D)S 2 3 Sensitivity Analysis of Drainage Design in Multi-layered Soil dd _ Sd (D + f D) 2 e e dr 2 Kr(D 3 dde _ d D 'dS ~ "5 3 + D )S 2 e 77 + 2 D -d 3 D + D 2 e 3 The percentage change in spacing due to the change in an individual parameter, while other parameters are held constant, can be calculated as: AS ~S~ R AS S RAR dRS R dS dK 1 s h AS AS S AS S x dh S h dS K AK 2 S AS S KxAKi S K dS h Ah AS 8K 2 S 2 K 2 OS D AD dD S D 2 D 2 2 2 2 OS K AK dK S K 3 K 3 3 dS dD Da AS S 4.2.4 dS 3 3 3 D AD S D 3 3 3 dS r Ar T dr S r Comparison of Sensitivities of Spacing between the Newly Developed Formula and the Kirkham Equation Table 4.4 is the summary of results of sensitivities of spacing to the individual parameters based on the new spacing formula, and the Kirkham equation for three-layered soils with the drains at the interface of the top two layers. The parameter values used in the calculations are: In case 1: i?=0.01 m/day; 1^=1.5 m/day; K =0.2 m/day; K =0.6 m/day; h=0.6 m; 2 D =1.0 m; D =2.0 m; r=0.1 m. 2 3 3 Sensitivity Analysis of Drainage Design in Multi-layered Soil Parameter Case 2 Case 1 (2) 12.00 -7.60 0.05 8.09 6.71 -1.35 1.18 0.56 2.55 (1) 20.98 -6.02 2.96 8.98 2.19 -0.72 0.88 0.88 0.87 Spacing (m) R K h K D Kz D r x 2 2 3 (1) 31.93 -6.25 0.18 6.42 4.53 1.05 1.53 1.53 1.00 78 Case 3 (2) 30.24 -6.22 0.30 6.22 5.02 1.62 1.16 0.84 1.03 (1) 27.16 -7.56 0.30 7.85 5.33 -1.80 1.93 1.93 2.01 (2) 25.08 -7.63 0.37 7.71 5.60 -1.17 1.87 1.12 2.12 Table 4.4: Comparison of sensitivities of spacing to individual parameters based on the new spacing formula and the Kirkham equation. In case 2: i?=0.01 m/day; Ki—Q.2 m/day; K 2=1.0 m/day; K =0.5 m/day; h—0.6 m; 3 D =2.0 m; £> =3.0 m; r=0.1 m. 2 3 In case 3: i?=0.01 m/day; K =0.2 m/day; if2=0.5 m/day; K =1.5 m/day; h=0.6 m; x 3 D =1.0 m; £> =3.0 m; r=0.1 m. 2 3 All parameters are assumed to change by 10%. In Table 4.4: (1) is the results calculated from the new spacing formula (the head loss above drains is considered) and (2) is the results obtained from the Kirkham equation. From table 4.4, it can be seen that, if a highly pervious layer is above the drains (Ki ^> K 2 and K ^> K , Case 1), the Kirkham equation is no longer reliable because x 3 the head loss above the drains is much more than negligible. The sensitivities of spacing to individual parameters from the two equations are very close to each other except to K and K when K ^> K and K 3> K . In both equations the spacing is very sensitive x 2 x 2 x 3 to the drainage coefficient and the water table height midway between the drains. In Sensitivity Analysis of Drainage Design in Multi-layered Soil 79 both equations spacing is not sensitive to the thicknesses of the lower layers, the drain radius, and the hydraulic conductivity of the third layer even when it is highly pervious. Chapter 5 UNCERTAINTY ANALYSIS OF DRAINAGE DESIGN In agricultural drainage design, varibility of soil properties and climatic conditions is an inevitable source of uncertainty in the design of a drainage system and in its performance. Such uncertainty incereases the difficulty for a farmer to decide whether to accept or reject the design. Therefore, it is desirable to estiamte the uncertainty involved in a particular design to provide the owner with necessary information, based on which the decision is made. The amount of water recharge to the soil depends on the balance of precipitation,evapotranspiration, surface runoff and deep seepage, all of which are of stochastic characteristics. Any estimate of the recharge rate is inevitably accompanied by a certain magnitude of variability. Soil properties, such as hydraulic conductivity and porosity, also exhibit a large degree of spatial and perhaps temporal variation (Freeze, 1975; de Vries, 1982). From the sensitivity analyses in previous chapters, it has been found that drain spacing is very sensitive to the soil hydraulic conductivity and the drainage rate. The current practice of drainage design is still mostly based on a deterministic framework by using equivalent climatic and soil parameters. However, the meanings of these equivalent values are not always so clear to agricultural engineers, hydrogeologists, and soil physicists. These equivalent values are chosen partly based on the actual conditions, but also largely on the experience and professional judgement of the individual designer. Much research has been done to incorporate uncertainty of parameters into stochastic and deterministic analyses (Bakr et al, 1978; Cordova and Bras, 1981; Cornell, 1972; 80 Uncertainty Analysis of Drainage Design 81 Dagan, 1979; Delhomme, 1978; Freeze, 1975; Sagar, 1978a and b; Smith and Freeze, 1979). Most of the above studies require either large amount of data, which is seldom, if ever, available in practice, or advanced mathematical and statistical knowledge which is often beyond that of many practicing designers. Some of the methods require excessive computations as prerequisites to solution of the problem. Even so, very few methods can provide an overall analysis of uncertainty. Cornell (1972) and Benjamin & Cornell (1970) introduced first and second order analysis of model and parameter uncertainty into civil engineering. Prasher (1982) applied first and second order analysis in drainage design. However, he only introduced the idea and very much work is left to be done before this method can be conveniently used by drainage practitioners. In this chapter, first and second order analysis of parameter uncertainty will be briefly illustrated and then applied in the uncertainty analysis of subsurface drainage design in homogeneous and two-layered soils. 5.1 First and Second Order Analysis of Parameter Uncertainty First and second order analysis only requires the first and second moments of the random variables and the first and second derivatives of functional relationship, which forms its bases of practical usefulness. It is usually very difficult, if not impossible, to obtain full distribution of the random variables as required by more complex and exact models, while it is relatively easier to obtain the first two moments. If the random variables are normally distributed, the first two moments can fully describe the characterizations of the random variables; otherwise this method will only provide approximate results. If coefficient of variation and/or degrees of non-linearity are relatively small, the method will provide quite accurate results (Cornell, 1972). Consider first the simplest situation. Assume that random variable, y, is functionally Uncertainty Analysis of Drainage Design 82 related only to another random variable, x, i.e. in which x has variance er . and mean fi . 2 x If the function behaves sufficiently well (i.e. it is continuous and differentiable at the mean), it can be expanded by Taylor series about the mean or expected value of x (fi )x y =/ W + /'(/*«)(* -/*») + - /O +• • • 2 where /'(/x ) and f"(fi ) are the first and second derivatives evaluated at u. , respectively. x x x If the coefficient of variation is relatively small, the above series can be truncated at the third term: V « fM + /'(/*«)(* - Mx) + ^f"(fi )(x - u. ) 2 x x Taking the expectation of both sides gives: ii = E(y) = £;[/(/*.) + /'(/*-)(*-/*-) + ^ / V - ) ( * - / ' - ) ] 2 y = f{u. ) + f'(/i )E[x x x - u. ] + ±f"ME[x x - Mx] 2 i.e. /*v = /(/*») + ^ / ' V « ) * f f To obtain the second moment of the random variable y (variance <r ) with the first 2 and second moments of x, the third term in the Taylor expansion is dropped: 2/~ Here, /*„ = /(/*«). /(/**) + /'(/**)(* - /**) Uncertainty Analysis of Drainage Design 83 Then the variance of y can be estimated as: <r = Var(y) = E[y - u. ] 2 2 y y = £[/(/i.) + / V . ) ( * - M . ) - / 0 O ] a So far,the mean, u- , and the variance, a , of y have been estimated as: 2 y (5.50) (5.51) [ / ' M M It should be noted that if the function is non-linear i.e., f"(/j, )^0, the mean value x of the dependent variable is not equivalent to deterministic solutions, and the vanriance of the dependent variable is only an approximation. The range of error is dependent on the non- linearity and the variances of the independent variables. —* More generally, if X represents the column vector of independent variables, and y the —* dependent variable functionally related to X'y = fix) then we have: (5.52) J where fJ'X'- column vector of the means of independent variables; AX'- first order difference column vector of indenpent variables: AXi = Xi- fi Xi Uncertainty Analysis of Drainage Design A X : second order difference matrix of independent variables 2 AX]: jth column of A X ; 2 B : transpose of the column vector of first order derivatives; x B\: transpose of the second order derivative matrix: df OXidXj 2 B\-\ j t h row vector of B . 2 And first and second moments of y are: 2 3 *l=B?CBi where C: covariance matrix of the vector of independent variables; dj = Cov(X ,X ) i where p(Xi,Xj) j Cov(X ,X ) i = j p(X ,X )<r * i j Xi Xi is the correlation coefficient of Xi and Xj Cj- is the jith column vector of C. For special cases where X are uncorrelated, we have: H = f(lx) + \Y.A C') T j Uncertainty Analysis of Drainage Design 85 where B T 2 is the transpose of column vector of second derivatives: df dX 2 B'2i 2 C is the variance column vector of X- As an illustration, let's consider cases where there are two independent variables, U and V, which may or may not be correlated: y = f(u,v) In equation 5.52, M VW ( U - nu AX (U-iiu) ( (U-nu)(V-p )\ 2 AX' v ) v (V - p ){U - pu) I 6_y dy v dy 2 au 2 2 \ auav In equation 5.53 and 5.54, (V-u. ) v avau \ 8y 2 av 2 / 2 Uncertainty Analysis of Drainage Design 86 Cov(U,V) >u C = { Cov(V,U) \ J Then: fly 2 \ L M JVU JU / \ ^ Cov(V,U) j c" _2 v) + -1 \fM + f Cov(v, r vu Cov(U,V) + fuv \ ft u) + f£yCov(u, v) + ft* v f(U, V) + ± \f»<r* + 2f Cov(U, V) + ft**} uv and I = (ft f ^ ft) al Cov(C/,V) \ \Cov(V,U) cr + f'VCo^U) / v ft) v ftrCov(U,V) j' + ft*} ( Ju f [ft (fuf^u 5.2 + 2ftftCov{U, V) + {ftfay Uncertainty Analysis of Steady Drainage Design in Homogeneous Soils Among the numerous steady state drainage theories available, the Hooghoudt equation is most commonly used by agricultural engineers (Luthin, 1978; Chieng, 1981). It can be written as: S* = ^(2d e + h) where d is the equivalent depth, evaluated from equation 3.2 or 3.3; K is the saturated e hydraulic conductivity; R is the drainage coefficient; h is the permissible midspan water table height above the drain level; and S is the drain spacing. Uncertainty Analysis of Drainage Design 87 For the purpose of evaluating the performance of a drainage system, the water table height h is chosen as the dependent variable. The water table height, h, is the most important parameter in drainage design. It is directly related to the soil workability and the moisture content in the root zone. For a given spacing, any uncertainty in evaluating other parameters such as K, R and d will result in uncertain performance of the drainage e system. The Hooghoudt equation (Eq. 3.1) can be rewritten as: [RS 2 h = ilK ~ +dl de ~ , ' (5 57) In the following analyses, only the drainge coefficient, R, and hydraulic conductivity, K, are to be considered as weakly correlated because in most locations where R is high, K is high and vice-versa (Prasher, 1982). A l l other parameters are to be considered as not correlated. 5.2.1 Drainage Coefficient R as Random Variable The drainge coefficient, R, is defined as the amount of water removed by a drainage system per unit surface area per unit time. It depends on the climatic conditions, the soil characteristics, the crop type and the trafficability requirement (Chieng, 1978). The proper design drainage coefficient is the one that will remove excess water before crop damage occurs. Drainage coefficients are developed by local experience or by analysis of soil water balance (Willardson, 1982). It is recognised that drainage coefficient determination are uncertain to some extent because of the complexity of nature and limitations of methods. A practical way of defining the degree of uncertainty in parameters is to estimate their first two moments by whatever means that are available. A brief review of methods to estimate means and variances from limited data can be found in Appendix B . It should Uncertainty Analysis of Drainage Design 88 be noted that there is still much work to be done in this aspect. For the analysis in this chapter, the first and second moments of random parameters are assumed to be already known. The mean and variance of the drainage coefficient are estimated, by some method, as R and Vav(R), respectively. The first and second derivatives of h with respect to R are equal to: dh _ h + 2d h dR~ 2R(h + d ) 2 e e dh _ (h + 2d h) dR ~ ~4R (h + d ) 2 2 2 e 2 2 3 e Therefore, the mean and the variance of the water table height are: E(h)\n = + h \§fiVK(R) , Var(/0|* 1. = (h + 2d h) .„ 2 2 e /r>N C^) Var(i2) 2 (h + 2d h) iR (h + d ) 2 2 e = 2 2 V e ^ R ) where \R denoted "due to R", and will be used throughout the thesis. Rearranging the above equations into dimensionless form yeilds: E(h)\ -h h _ R ( l + 2f ) V a r ( £ ) 8(1 + f ) R 2 3 Var(h)\ R h 2 _ (l + 2 2^) Var(R) 4(l + f ) (5.58) 2 2 R 2 (5.59) Figure 5.13 and 5.14 show that the effects of the randomness of the drainage coefficient, R, on the mean water table height decreases, and that on the variance of water table height increases, as % increases. 89 Uncertainty Analysis of Drainage Design 0.00 -0.04 0£ M re ft; > -0.08 1 1 -0.12 63 -0.10 -0.20 0.0 10.0 d /h t Figure 5.13: Expected water table height due to uncertainty of drainage coefficient. 1.0 5. • > IA I* 10.0 d /h e Figure 5.14: Variance of water.table height due to uncertainty of drainage coefficient. Uncertainty Analysis of Drainage Design 5.2.2 90 H y d r a u l i c C o n d u c t i v i t y as R a n d o m V a r i a b l e W i t h M e a n of K and Variance of Var(K) The expected water table height and its variance are: m \ K = h ^ Var(K) + 2 Var(/0k = ( | | ) V a r ( ; 0 2 respectively, where dh dK h + 2d h 2K(h + d ) 2 e e d h _ (h + 2d h)(3h + Qhd + U\) 9K*~ 4K (h + d ) 2 2 2 e e 2 3 e Therefore, the mean and the variance of water table height in dimensionless form are: E(h)\ .-h_ h ( l + 2t)[3 + 6 t + 4 ( t ) ] V a r ( X ) ~ 8(1 + f ) K 2 K 3 2 V a r ( f e ) | ^ ( l + 2t) Var(if) h 4(1 + f ) K 2 g 2 2 2 K > respectively. Figure 5.15 and 5.16 show that the effects of uncertainty in the hydraulic conductivity on both the mean and the variance of water table height increase as ^ increases. 5.2.3 C o r r e l a t e d H y d r a u l i c C o n d u c t i v i t y K and Drainage Coefficient R w i t h a covariance of Cox(R, K) For the convenience of writing the equations of the mean and the variance when all parameters are considered, only the additional portion of the mean and the variance of the water table height caused by the correlation between K and R will be calculated here. The total effects of K and R will be the sum of the results obtained in the previous two sections and in this section. Uncertainty Analysis of Drainage Design 91 Uncertainty Analysis of Drainage Design 92 From Eq. 5.57, we have: dh dRdK h •+ 2d h 4:KR(h + d ) 2 2 e e dh h + 2d h dR ~~ 2R(h + d ) 2 e e dh h + 2d h 2K(h + d ) 2 e 8K ' e The additional effects on the mean and the variance of the water table height are: E{h)\ = (KxR) h + J^x ^(K,R) C 4KR(h + d ) K e V™(h)\ = {KxR) _ ~ where \(KxR) ' ; 2{^){^)Cov(K,R) (h + 2d h) -2KR(h + d y 2 2 e C o v ( K e ' R ) denotes "due to the correlation of i f and R". The relationship between covariance and correlation coeffcient is: Cov(K,R) in which p(K,R) = p(K,R)a a K R is the correlation coefficient of K and R. If K and R are independent, p(K,R) = 0 Rearranging the equations of mean and variance gives: E(h)\ \ {KxR){ K x R 1 -h_ -h +2 2 ff ll + 4 Var(/Ql(*xB) _ h 2 (! + f ) Cov(K,R) K ( 5 g 2 ) R (l+2t) Cov(if,i?) 2(1 + ^ ) KR 2 2 { ' ' Uncertainty Analysis of Drainage Design 93 The negative sign in both equations indicates that R and K have opposite effects on the mean and the variance of the water table height. It can be seen in Figure 5.17 and 5.18 that effects of the correlation of R and K on both the mean and the variance of the water table height increase as * increases. 5.2.4 Equivalent Depth as a Random Variable with Mean of d and Variance e ofVar(d ) e The first and second derivatives of h with respect to d are: e dh h dd h e +d e d h _ h + 2d h 2 2 e dd ' (h + d y 2 e The mean and the variance of the water table height can be estimated as: E(h)\ = de h + ^Var(d ) e 2{h + d y K e } e Var(Mk = (J|) Var(4) 2 h (h + d ) V a r ( i ) 2 e e or, in dimensionless form: E(h)\ -h h ^ ( f ) ( l + 2t)Var( i ) 2(1 + ^ ) 3 d\ 2 de h t 2 + d e (5.64) (5.65) 2 e Figure 5.19 and 5.20 show that the uncertainty of equivalent depth becomes more important to the mean and variance of water table height as increases. 94 Uncertainty Analysis of Drainage Design -0.2 •ft! d /h e Figure 5.17: Expected water table height due to the correlation of the hydraulic conductivity and the drainage coefficent. 95 Uncertainty Analysis of Drainage Design o.e-\ d /h e Figure 5.19: Expected water table height due to uncertainty of equivalent depth 0 2 4 ; 6 8 10 djh Figure 5.20: Variance of water table height due to uncertainty of equivalent depth Uncertainty Analysis of Drainage Design 5.2.5 96 Estimating Variance of Equivalent Depth from the Mean and Variance of the Depth to an Impermeable Layer Equivalent depth, d , is dependent on the spacing, the impermeable layer depth and the e drain radius. For a given spacing and a drain radius, the variance of the "equivalent depth can be estimated from the mean and the variance of the depth to the impermeable layer. d can be calculated by equation 3.2 or 3.3 (Moody, 1966): e For 0 < f < 0.3 1 + i ( I l n * - a) where a = 3.55-1.6| + 2(|) 2 and for f > 0.3 ' |(lnf-1.15) The first derivative of d with respect to d is: e For 0 < | < 0.3 dd _ 1 - f f - 1 . 6 ( f ) e dd 2 + 4(f) [l + f ( f l n f - o : ) ] 3 2 and for f > 0.3 dd, dd=° dd Therefore, the variance of d can be estimated as: e For f > 0.3 Var(4) = 0 and for 0 < f < 0.3 VarK) = [1 _ i i _ i Q(Ay _i_ 4(iy»12 * , } l f Var(d) 1 5 S d a j J Uncertainty Analysis of Drainage Design 97 Neglecting the square of ^ and higher terms and rearranging the above equation gives; Var(ci )_ (I-!!) e d l Var(d) 2 [l + K f l n f - 3 . 5 5 ) ] ' d 2 (5.66) Let Var(4) _ Var(ti) d l ~~d?~ thus 8 d ,, d ,8. d . ( l - ; j ) = [l + 5 ( ; > « 5 - 3 . S 6 ) ] ' r , Solving the above equation yields: or - = 1.483 r i.e., when d = 0 or d = 1.483r Var(4) _ Var(fi) and for 0 < <f < 1.483r Var(4) Var(rf) d d 2 e 2 which is somewhat unreasonable. This ambiguity can be avoided by assuming that: Var(4) _ Var(ci) dl ~ d 2 for 0 < d < 1.483r, where r is the radius of the drains. 5.2.6 A l l or Some of the Parameters as Random Variables with Known Means, Variances and Covariances From section 5.1, we have L 3 =1 Uncertainty Analysis of Drainage Design 98 Var(/i) = 5 f w here dh 1 dh 2 dR C B ah 2 2 dd 8K8R 2 ah ah dRdK dK 2 2 ah ah ah 2 2 dRdd 2 dd dK8d 2 e e I dR e ah 2 2 ah \ dR ah dK B , = ah \ Var(tf) ( c = Cov(K,R) 0 V dd e I Cov(R,K) 0 Va.r{K) 0 0 Var(ci ) ) e Rearranging the above equations gives: dh dh 2 2 E{h) = h + \[^Var(R) + Var(M = ( ^ ) V a r ( t f ) + ( | ^ ) V a r ( 7 Y ) + ( ^ ) V a r ( A ) + 2 2 2[?±)(^)Cov(K,R) 2 From the results of analysis in this section, the above equation can be written in the form: E(h)-h h = E(h)\ - h R + h E{h)\ -h h Var(fe) Var(fe)| ~~hP~ ~ h 2 K fl + , Var(h)\ h 2 K E(h)\ -h h {KxR) + E(h)\ - h de h , Vax(h)\ TT~ {KxR) + , Var(fe)] ~h?~ de ( 5 ' 6 8 ) in which various terms are correspondingly evaluated from equation 5.58 through 5.66. Uncertainty Analysis of Drainage Design 5.3 99 Uncertainty Analysis of Transient Drainage Design in Homogeneous Soils In the cases where steady state drainage design is inappropriate, transient drainage theory should be used. For saturated flows and for the cases where the D-F assumptions are applicable, the water flow toward drains in the soil can be described by the Boussinesq equation (van Schilfgaarde, 1974): f d k - K (h )+B d 9h where / is the drainable porosity, h is the water table height above an impermeable layer, K is the saturated hydraulic conductivity, x is the horizontal distance from the drain center, t is the time and R is the net recharge rate. The Boussinesq equation is non-linear. Golver (Dumm, 1954) and Tapp & Moody (Dumm, 1964) solved this equation by linearization, van Schilfgaarde (1963; 1964), on the other hand, solved it without linearization and obtained the following equation for the case without recharge: fS 9Kd 2 t = m (2d +m ) m (2<f +m ) 0 e 1 e 1 e 0 in which t is the time during which the water table drops from initial height of m to m j . 0 d in the original van Schilfgaarde solution has been replaced by the equivalent depth d e to provide more practical results. For a given spacing, S, the required time, t, for the water table to recede from m 0 to m-i is determined by the drainable porosity, / , the hydraulic conductivity, K, and the equivalent depth, d . However, because of the inevitable uncertainty in the values of these e parameters,it can be expected that there will be some uncertainty in the drawdown time, which may result in either a longer period of waterlogging or water deficiency before the Uncertainty Analysis of Drainage Design 100 next recharge. It is therefore desirable to estimate the variation of the drawdown time caused by the uncertain input values of parameters. First and second order analysis of parameter uncertainty will be applied to the following cases: (1) the drainable porosity as a random variable; (2) the hydraulic conductivity as a random variable; (3) the correlation of the drainable porosity and the hydraulic conductivity; (4) the equivalent depth as a random variable; (5) all of the above. 5.3.1 Drainable Porosity as Random Variable with Mean of / and Variance of Var(/) The first and second derivatives of t with respect to / are: d t _ t df'l 8H n 8 F = ° respectively. Therefore, the mean and the variance of the required time are: E(t)\f = t + ^ V a r ( / ) = t dt t (-) Var(/)=:-Var(/) 2 Var(i)| 2 / = : or, in dimensionless form: E{t)\,-t t Var(t)|/ . t 2 = = 0 (5.69) Var(/) P (5.70) Uncertainty Analysis of Drainage Design 5.3.2 101 Hydraulic Conductivity as a Random Variable with Mean of K and Variance of Var(if) The first and second derivatives of t with respect to K are: dt _ t d~K ~~K OH dK 2 2t , K 2 Therefore, the mean and the variance of required time are: = t + i ^ V a r ( / f ) = t + ^Var(tf) Var(t)| , = ( ^ ) V a r ( i f ) = ^ V a r ( Y ) 2 i J f or, in dimensionless form: E{t)\ -t _ Var(if) t K Var(t)|jc _ Var(Jr) t K K 2 2 5.3.3 2 (5.71) (5.72) Correlated Hydraulic Conductivity and Drainable Porosity with a Known Covariance of Cov(K,f) There exists a relationship between the drainable porosity (or specific yield) and the soil hydraulic conductivity (Talsma and Haskew, 1959). Soils with a higher drainable porosity usually have a higher hydraulic conductivity. The correlation coefficient of the drainable porosity and the hydraulic conductivity can be estimated from field data for the region. The derivatives required for the calculation of the effects of the correlation on the mean and the variance of the drawdown time t are: dt dfdK 2 t_ ~ ~~Kf Uncertainty Analysis of Drainage Design 102 8t _t df~l at •_ t a~K ~~K Therefore, the mean and the variance of the drawdown time are: = * + ^Cov(ff,/) E(t)\(K*f) = Var(i)|(^x/) t-jLcov{K,f) = ( § ) ( ^ 2 Kf C o v ( ^ / ) ;Cov(if,/) or, in dimensionless form: E(t)\ -t [Kxf) = Cov(if,/) t Kf Var(*)| ( J f x / ) _ 2Cov(if,/) (5.73) (5.74) where C o v ( i f , / ) is the covariance of i f and / . If the correlation coefficient is known: where p(K,f) is the correlation of i f and / , 07 is the standard deviation of / , <TK is the standard deviation of i f , and: <r, = v/Var(/) = ^/Va^if) Then the mean and the variance of the drawdown time can be written as: m ^ n z l ^ ^ ^ ) ^ Mf^i = ^Y^}^ iKJ) , ( 75) . (5 76) Uncertainty Analysis of Drainage Design 5.3.4 103 Equivalent D e p t h d as a R a n d o m V a r i a b l e w i t h M e a n of d a n d V a r i e e ance of Var(d ) e d is dependent on the spacing, the drain radius and the impermeable layer depth. For e a given spacing and a drain radius, the variance of d can be calculated from the mean e and the variance of the depth of the impermeable layer (See section 5.2.5). The first and second derivatives of t with respect to d are: e dt t 6d d e 2(m-o — mi)t e dt {2 2t 2 ddj = d + m )(2d e 1 +rn )\n[Z%t Z]} + e 0 + 4i(mo — m i ) ( m m i + m,id + mod ) 0 d? ~ d (2d + y(2d e e e mi e e +mo)Mn[™g ^j e : Therefore: , ,2t X Var(0k = At(m - m ) ( m m I m ^ + m d ) 0 1 4 ( 2 4 + m 0 dl (^) d 2 V a r 0 2 1 ( 2 r f 0 + m e 0 ) 2 e l n [ ^ ± ^ ] / (^) (2d +m )(2d m )\n[^^ ] e e 1 e+ 0 } ) or, in dimensionless form: E(t)\ -t de ' V a r ^ 2 ( ^ - ^ ) ( l + ^ = • (2 = + ) ( t + 1 mT 2 2 1 ^ ( 4 ) + ^) ) 2 M ^ ) d { 1 (5.77) l ( 5 ( 2 t + )(2^ + ) 1 1 l n (i^Tl) . 7 8 ) ' where d can be calculated from Equation 3.2 or 3.3 and Var(d ) an be estimated from e Equation 5.66. e C Uncertainty Analysis of Drainage Design 5.3.5 104 All or Part of the Parameters as Random Variables with Known Means, Variances and Covariances From section 5.1, we have Var(fc) = S f C B x where at at dKdf dd df dt dfdK at dK at ad dK at at dfdde dKdd at ad 2 I df B\ = 2 2 e 2 2 2 2 \ at 2 2 e 2 ( 2 2 e J *^ df at dK B = x at \ ( Var(/) C = 8d Cov(/,if) 0 Vai(K) 0 Cov(KJ) 0 V ) e 0 Var(4) j Rearranging the above equations gives: l dH„ , t dt „ 2 , „ v dt dH 2 Var(t) = ( ^ ) V a r ( / ) + ( ^ ) V a r ( X ) + ( ^ ) V a r ( r f ) + 2 ( ^ ) ( ^ ) C o v ( X , / ) 2 2 2 e From the analysis results in this section, we can write the above equations in the forms: E(t) - t t = E{t)\ -t f + E(t)\ -t K t + E(t)\ - t , E(t)\ (Kxf) - t dc t t (5.79) Uncertainty Analysis of Drainage Design Var(Q - ^ r - VarCQj, - — r r 105 Vax(t)\ , Var(Q| K ~ + — — j2 + ( y x / ) , Var(Q|^ +—^i— ( - °) 5 8 in which the various terms are correspondingly evaluated from Eequation 5.69 through 5.78. 5.4 U n c e r t a i n t y A n a l y s i s of Steady Drainage Design i n T w o - l a y e r e d Soils The first and second order analysis of parameter uncertainty can also be applied to drainage design in multi-layered soils. However, more difficulties can be expected in estimating the means and the variances of various parameters and their covariances, if any exist. In this section, only uncertainty in two-layered soils with the drains lying at the interface of the layers will be analyzed, although uncertainty analysis of drainage design using the Kirkham equation, the generalized Hooghoudt-Ernst equation or the new spacing formula developed in Chapter 4 can be carried out by first and second order analysis in a similar way. The Hooghoudt equation is the most simple and commonly used theory in steady state drainage design in two-layered soils with drains at the interface of the two layers. It mainly distinguishes the flow regions above and below the drain level. It can be expressed as: Ah S = — (2K d K 2 2 e + K h) 1 where K\ and K are the saturated hydraulic conductivities of the soils above and below 2 the drains, respectively. All other parameters are as defined in section 5.2. The midspan water table height, h, is again chosen as the dependent variable for the convenience of evaluating the performance of a drainage system with a selected drain Uncertainty Analysis of Drainage Design 106 spacing. The Hooghoudt equation for two layered soils can be rewritten as: The first and second order analysis will be applied to the following cases: (1) random drainage coefficient R; (2) random hydraulic conductivity of the soil above the drains K\\ (3) random hydraulic conductivity of the soil below the drains K ; (4) correlation 2 of R and K\\ (5) correlation of R and K ; (6) correlation of K\ and K ; (7) random 2 2 equivalent depth d ; (8) all of the above. e 5.4.1 Drainage Coefficient as Random Variable with Mean of R and Variance of Var(i?) From equation 5.81, the first and second derivatives with respect to R are: dh _ K h + 2K d h ~dR ~ 2R(K h + K d ) 2 x 2 1 dh _ 8R ~~ 2 Ki{K h + 4# (if /i + 2 e K d hf Kdf 2 x 2 e 2 2 1 e 2 e Thus the expected mean of h and its variance are: W(h\\ h K\{K\h * = ~ 8R ( h + 2 E { k ) l h 2 Kl v (h\\ V a i ( / l ) | f l (Kih + = 4R*(K + 2 lh K d h) K d r 2 2 e 2 V a T { R ) e 2K d h) , Kd) ™ 2 2 e T 2W + 2 {R) e or, in dimensionless form: E{h)\ -h R h ^ [l + 2 ( f ) ( t ) ] V a r ( i ? ) 2 8[l + ( f ) ( f ) ] 3 Var(fc)|fl h 2 i? 2 fl + 2(ff)(t)] Var(i?) (5.82) 2 4[l + ( f ) ( t ) ] 2 R 2 (5.83) Figure 5.21 and 5.22 show that the effect of the randomness of R on the mean of h decreases, and that on the variance of h increases, as jj£ and/or ^ increase. 107 Uncertainty Analysis of Drainage Design O K2/K1S0.2 V K2/K1=3.0 Figure 5.21: Expected water table height due to uncertainty of the drainage coefficient. K2/K1-0.2 K2A1»1.0 V K2/K1-S.0 d /h e Figure 5.22: Variance of water table height due to uncertainty of the drainage coefficient. Uncertainty Analysis of Drainage Design 5.4.2 108 Hydraulic Conductivity of the Soil above Drains as Random Variable with Mean of K t and Variance of V a r ( i f i ) From the Hooghoudt equation, the first and second derivatives with respect to K is: x dh h W ~~ x dh _ 2 2 ~2(K h+-K d ) 1 3K h 2 + 4 x d~K\~ AK {K h 4K d h 3 2 + 2 x e e Kdf 2 e Therefore, the mean and the variance of h will be: Mh\\ m i K u , M.h* + 4K d h SK ( h K d 3 2 > = h + e w 2 Kl + 2 e y Y ^ respectively. Rearranging the above equations in dimensionless form gives: E(h)\ -h _ Ki 3 + 4(ff)(£) h 1 8[l + ( f ) ( t ) ] Var(h)\ Kl h Var(if ) _ 1 3 Var^) 4[1+ 2 K 2 if 2 (5.84) (5.85) Figure 5.23 and 5.24 show that the effect of the randomness of Ki on both the mean and the variance of h decreases as jjj^- and/or ^ increase. This is because that the flow region above the drains becomes less important when 5.4.3 and/or ^ become larger. Hydraulic Conductivity of the Soil below Drains as Random Variable with Mean of K 2 and Variance of V a r ( i f ) 2 The first and second derivatives of h with respect to K are: 2 dh d~K dh ~ Kh + K d e = 2 t 2 e 109 Uncertainty Analysis of Drainage Design ; djh Figure 5.23: Expected water table height due to uncertainty of the hydraulic conductivity of the soil above the drains. ; djh Figure 5.24: Variance of water table height due to uncertainty of the hydraulic conduc- . tivity of the soil above the drains. Uncertainty Analysis of Drainage Design dh 110 _ 2 K h d +2K hd 2 1 dK ~ {K h 2 2 3 e 2 + x e Kdf 2 e Therefore, the mean and the variance of h can be estimated as: Ff \\ , h , K h d\ + 2 x 2K hd\ 2 Rearranging the above equation in dimensionless form: E(h)\ [1 + 2 ( f )(£)](% - h Ka h 2 Var(fc)k _ ^ (t) (ff) 2 . 2 [ l + (g)(t)] + Var(K ) 2 3 Var(if ) 2 (5.87) 2 Kt)] 2 (5.86) K\ Figure 5.25 and 5.26 show that the uncertainty of K becomes more significant to 2 both the mean and the variance of h as j£ and/or ^ increase. This is because the flow region below the drains becomes more important when 5.4.4 and/or ^ becomes larger. Correlation of Drainage Coefficient R and Hydraulic Conductivity, K~i, of the Soil above Drains with covariance of C o v ( i f , i ? ) 1 From the Hooghoudt equation, the necessary derivatives needed for the calcualtion of the mean and the variance of h are: dh OR Kxh + 2K d h 2R{K h + K d ) dh dKx h 2(K h + K d ) dh dRdKx 2 2 2 1 e 2 e 2 i 2 e (Kxh + 2K d ) h AR(K h + K d f 2 2 x 2 e 2 e Therefore, the mean and the variance due to the correlation of R and K\ are: 111 Uncertainty Analysis of Drainage Design 1.<H 2 K2/K1=0.2 K2/K1-1.0 V K27K1=3.0 10.0 Figure 5.25: Expected water table height due to uncertainty of hydraulic conductivity of the soil below drains. g K2/Kl-0.2 K2/K1-0.2 Kg^KI-JiO V K2/Kt»9.0 I.On n3 > t-t > 10.0 Figure 5.26: Variance of water table.height due to uncertainty of hydraulic conductivity of the soil below drains. Uncertainty Analysis of Drainage Design v 112 ( ^ +2 ^ 4 ) ^ (h\\ or, in dimensionless form: i^faxfl)-/* ^ * [l + 2 ( f )( f)] Co (K R) 4[l + ( f ) ( t ) ] 3 d 2 V ifiiZ Varkif.xfl) l + 2(ff)(t) Cov(J<i,fl) ^ 2[1+ ( & ) ( £ ) ] ' if # 2 Figure 5.27 and 5.28 K x (5.88) u (5.89) a and R have opposite effects on both the mean and the variance of h, which generally decrease as | ^ and/or ^ increase. 5.4.5 C o r r e l a t i o n of Drainage Coefficient R a n d H y d r a u l i c C o n d u c t i v i t y , K of the Soil below D r a i n s w i t h covariance of 2 Cov(K ,R) 2 From the Hooghoudt equation, the necessary derivatives needed for the calcualtion of the mean and the variance of h are: dh^ _ K h + 2K d h d~R " 2R(K h + K d ). 2 x 2 e 1 dh 8K e d„h Kh + Kd 2 1 dh 2 {K h + 2 8R8K 2 e 2K d )K d h 2 x 2 e 2R(K h + 2 2 e Kdf 1 2 e Therefore, the mean and the variance due to the correlation of R and K are: 2 , E{h)\ » p / , (K h + =h- n {K 2 v mi Var{h)\ [K2XR) =- R { K i 2K d )K d h„ Cov(K ,R) 2 1 R) 2 h + K 2 (K h + { K l h e e 2 e f 2 2K d )d h , Cov(K ,R) 2 l R d 2 + K 2 i e e ? e r 2 113 Uncertainty Analysis of Drainage Design O A o.o a K2/K1=0.2 K2/K1-1.0 V K2/K 1=5.0 « Figure 5.27: Expected water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil above drains. d /h s Figure 5.28: Variance of water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil above drains. Uncertainty Analysis of Drainage Design 114 or, in dimensionless form: E(h)\ - h {KaxR) (fe) (^) [l + 2(^)(t)]Cov(g iZ) 2 = a (5.90) a> * 2[l + ( f ) ( t ) ] Var| « xB) _ KR 3 2 ( t ) ( t ) [ l + 2(ff)(t)]Cov(Jf ,fi) (5.91) h [l + ( t X l f ) ] ^ Figure 5.29 and 5.30 show that R and K have opposite effects on both the mean and ( 2 2 2 2 2 the variance of h, which increase as 5.4.6 and/or ^ increase. C o r r e l a t i o n of H y d r a u l i c C o n d u c t i v i t i e s of the Soil above D r a i n s , K , x and below D r a i n s , K , w i t h covariance of C o v ( i f , K ) 2 1 2 From the Hooghoudt equation, the necessary derivatives needed for the calcualtion of the mean and the variance of h are: dh _ h dK\ ~ ~ 2(K h + K d ) dh dh W ~K h + K d ~d h d h (2K h + 3K d ) dK 8K 2{K h + K d f 2 x 2 e e = 2 x 2 2 e 2 e 2 x x 2 x 2 e e Therefore, the mean and the variance due to the correlation of K and K are: x „, , E(h)\ , . d h (2K h = h+ { 2 n e {M2) 2 x Kih + 3K d ) 'Cov(gi,*») 2 + w 2 e d h 3 Var(MI( or, in dimensionless form: E(h)\ M 2 ) = - h_ {KlxKa) h Var| c x ) - l {Kih K2de)2 Cov(K ,K ) x ( f e ) ( t ) t + 3(^X^)1 2[l + ( g ) ( i ) ] » 2 Co (K ,K ) KK 2 V X x 2 (5.92) 2 Cov(X ,Jf ) (5.93) >> [l + (ft)(1f)] * Figure 5.31 and 5.32 show that effects of the correlation on both the mean and the (i lX 2 2 = 1 2 K l 2 K variance of h first increases and then decrease as jj^- and/or ^ increase. 115 Uncertainty Analysis of Drainage Design O A 0.0 K2/K1=0.2 K2/K1-1.0 V K2/K1=3.0 -0.2-1 > o V -0.4 -0.6 - 0 . 8 H -1.0 0.0 10.0 Figure 5.29: Expected water table height due to correlation of drainage coefficient and hydraulic conductivity of the soil below drains. 2 0.0 K2/K1-0.2 K2/K1-1.0 V K2/K1"5.0 -0.44 -0.8 -1.2H -1.6 -2.0 4 10.0 Figure 5.30: Variance of water table height due to correlation of drainage coefficient and hydraulic conductivity of the,soil below drains. 116 Uncertainty Analysis of Drainage Design K2/K1=0.2 V K2A 1*3.0 10.0 djh Figure 5.31: Expected water table height due to correlation of hydraulic conductivity of the soil above drains and that below drains. 2 > K2/K1-0.2 K2/K1-1.0 V K2/K1-3.0 -0.1- o U -0.2 -0.3 Figure 5.32: Variance of water table height due to correlation of hydraulic conductivity of the soil above drains and that below drains. Uncertainty Analysis of Drainage Design 5.4.7 117 E q u i v a l e n t D e p t h as R a n d o m V a r i a b l e w i t h M e a n of d and Variance e ofVar(de) The first and second derivatives of h with respect to d are: e dh _ Kh dd~ ~ ~K h + K d 2 e v 2 e d h _ K\{K.h + 2K d )h dd ~ (K h + K d f 2 2 e 2 x 2 e Therefore, the mean and the variance of h due to the randomness of d are: e KKK.h E { h ) ^ = h + 2K d )h 2 e 2(K h + K d y + 1 2 V a r K ) e or, in dimensionless form: E(h)\ de - _ ( £ ) ( t ) [ l + 2(fe)(fr)] 2 * 2 2[1+ V"(&)k ^ 2 (t) (t) 2 = V a r 4 Var(4) 2 [1 + (ff)(t)] 2 d l K ) (5.94) 2 (5.95) Figure 5.33 and-5.34 show that the uncertainty of d becomes more significant to e both the mean and the variance of h as | ^ and/or ^ increase. This is because the flow region below the drains becomes more important when ^ and/or ^ becomes larger. Equations 5.94 and 5.95 are exactly in the same form as equations 5.86 and 5.87, respectively. This again shows that the thickness of the equivalent layer is as important as the hydraulic conductivity of this layer to the drainage design. Uncertainty Analysis of Drainage Design Figure 5.33: Expected water.table height due to uncertainty of equivalent depth.. Figure 5.34: Variance of water table height due to uncertainty of equivalent depth. Uncertainty Analysis of Drainage Design 5.4.8 119 A l l or Part of Parameters as Random Variables with Known Means, Variances and Covariances From section 5.1, we have: E(h) = h + \ JZ{BlC ) 3 Vai(h) =B{•T C B x where dh dR dh B -x dh 8K 2 dh ah dR ah dKidR dh dK dR a?h dd dR dh dRBKi ah dK dh dK2&K i 2 ah ad aK ah dK ah dd dK ah dK2 dd ah ad 2 2 2 2 B 2 2 ah ah dK \dK"2 2 8RdK 2 2 2 l 2 2 8h 2 e 2 2 dK\dd e 2 2 e e Cov( K R) 0 lt Cov(i?,if ) Var(ifO 1 0 X 0 Cov( ftT ,if ) 0 Var(if ) 0 0 Var(d ) J Cov( K ,K ) 0 V e 2 Sh dRdd Var(R) C 2 2 2 , \ e 2 2 1 2 0 e } Rearranging the above equations gives: \8h dh 2 E { k ) = fe + 2 f e V dh 2 r a r ( ^ ) dd + ^ V a r 2 ( ^ l ) + 5X| V a r ( 8R8K X dh -Cov{K ,R) '8R8K 2 2 8h •Cov(K ,K ) +2 8K 8K_ 2 x 2 2 ^ 2 ) + Uncertainty Analysis of Drainage Design Var(fc) = 120 (|^) Var(#) + ( ^ - ) V a r ( i f ) + ( ^ ) V a r ( i f ) + 2 2 2 1 (^) Va K) 2(§)(^-)Cov(i, 2 r 2 + 1 ) #) + Using the results of the previous analyses in this section, we can write the above equations as: E(h)-h h = E(h)\ -h h E(h)\ R | {KlxR) E(h)\ -h h - h E(h)\ | h Var(fe) h 2 = Kl (K2xR) + V*r(h)\ h Var(h)\( h R 2 | {KlxK2) Var(h)\ Var(fe)| Var(h)\ h h h _ Var(h)\ _ Vai(h)\ h h g2 2 KlxR) | + Kl 2 2 E(fe)|jf — h £ ( ^ - 6 h h - h E(h)\ - h h~ l~ | 2 - 9 6 ) de 2 2 {K2xR) + ( 5 , [KlxK2) + 2 1 j in which various terms are evaluated from equation 5.82 through 5.95. 5.5 A n Example Showing How Uncertainty Analysis C a n Be Used in Drainage Design In this section, an example of uncertainty analysis of drainage design in homogeneous soils by using first and second order analysis of parameter uncertainty is presented. The input information of each parameter is arbitrarily made up but this does not affect the usefulness of the analyses. Appendix B briefly reviews the methods for estimating mean, variance and covariance from limited field data. Then an effort is made to show how the uncertainty analysis can be applied to practical drainage design. Uncertainty Analysis of Drainage Design 5.5.1 121 A n Example The values of parameters used as input information are as follows: Design water table height at middle point between the drains: h=0.6 m; Hydraulic conductivity: mean if=0.5 m/day, variance Var(if)=0.05 (m/day) ; 2 Drainage coefficient: mean i?=0.01 m/day, variance Var(i?)=1.5 x 10~ (m/day) ; 5 Correlation coefficient of drainage coefficient and hydraulic conductivity: p(K,R) 2 = 0.3; Depth to the impermeable layer: mean d=5 m, variance Var(d)=1.5 m ; 2 Drain radius: r=0.1 m. From the Hooghoudt equation by the deterministic design concept, the drain spacing needed is 5=24.0 m, the equivalent depth d =2.09 m. Using equations 5.58 through 5.66, e we can obtain the following results: E(h)\ -h h R = -0.013 Var(fe)| h E(h)\ -h h = 0.118 fl 2 = 0.160 K Var(fc)l* . = 0.158 E(h)\ - h KxR h Var(fe)lg h = -0.082 xfl 2 Var(rf )| e E(h)\ d = 0.002 - h de h = -0.023 = 0.001 Var(/Qk = 0.001 h 2 Uncertainty Analysis of Drainage Design 122 Then we can calculate the mean and the variance of h by Equation 5.67 and 5.68 as: E{h)-h h Var(h) -0.013 + 0.160 - 0.023 + 0.001 = 0.125 = 0.118 + 0.158 - 0.082 + 0.001 = 0.195 i.e. E[h) = 0.675 m Var(fc) = 0.070 m 2 The standard deviation of h from the mean is: SD = 0.265 m The variance or standard deviation of water table height may be considered as the risk involved in the drainage system design due to uncertainties about the nature of the problem. If we assume a normal distribution of the water table height, it implies that there is a 90% probability that h would not exceed 1.013 m (mean+1.28SD) and a 90% probability that h would exceed 0.337 m (mean— 1.2SSD), or there is an 80% probability that h is between 0.337 m and 1.013 m if the drainage system is designed based on the deterministic concept. 5.5.2 Application of Uncertainty Analysis in Drainage Design From the above example, it should be noticed that the mean water table height is actually larger than the design water table height (0.675 m and 0.600 m, respectively), which means the drain spacing based on the deterministic model (the Hooghoudt equation) is inadquate in removing the excess water before damage occurs due to the variability of soil and climatic conditions. The higher water table may result in a longer waterlogging Uncertainty Analysis of Drainage Design 123 period and thus reduce the crop yield. Through uncertainty analysis, we can correspondingly adjust the drain spacing so that the drainage system can perform better. On the other hand, even if the mean water table height meet the requirement, there is still a considerable probability that the water table will be higher than required due to the uncertainty of parametes, which may be regarded as a measure of the risk involved in selecting that design. From the results obtained previously, the increase of h is: Ah = E(h) - h = 0.675 - 0.6 = 0.075 m In order to lower the mean water table height, /inew is used in design (h ^ = h): Q h n e w =h o l d - Ah = 0.6 - 0.075 = 0.525 m Repeating the procedure in the above example, we have: spacing 5=21.64 m, equivalent depth d =1.97, and e E(h) = 0.592 m Var(fc) = 0.055 m 2 Standard deviation SD =0.234 m. This time the mean water table is close to the required height (0.6 m) (if still not the same procedure may be repeated). If a normal distribution of the water table height is assumed, the probability that the water table would not exceed certain heights are summarized in table 5.5. Table 5.5 will give the farmer an idea on approximately how much, the risk is involved in selecting the design for the particular field. The final spacing can be adjusted to be wider or narrower, depending on the risk level that the owner of the field would like to accept. Uncertainty Analysis of Drainage Design 124 Table 5.5: Probability that water table would not exceed certain height. Mean Standard Deviation Height (m) 0.592 0.710 0.826 0.892 0.943 1.060 0.592 m 0.234 m Probability 50% 69% 84% 90% 93% 98% Chapter 6 CONCLUSIONS A N D RECOMMENDATIONS This chapter summarizes the study that has been carried out on sensitivity and uncertainty analysis in the design of a subsurface drainage system. Conclusions will be drawn based on the study and recommendations will be made where appropriate for further research work in this or related aspects. The study started from the simplest subsurface drainage situation with homogeneous soils and steady state criteria. Sensitivity analysis on drain spacing to other design parameters was performed on the most popularly used steady state drainage theory, the Hooghoudt equation, for homogeneous soils to identify the dominant variable or variables among all drainage parameters. Then several transient state spacing formulae were reviewed and compared and the van Schilfgaarde solution of the Boussinesq equation was chosen for the purpose of sensitivity analysis. A comparison was made between sensitivities of spacing in steady state and transient state to corresponding parameters through an example. The Hooghoudt equation and the generalized Hooghoudt-Ernst equation for two- layered soils with the drains at the interface of the layers were analysed and their results compared. A new spacing formula was developed for the situations of two-layered soils with different hydraulic conductivities below the drains. Sensitivity analysis was carried out based on the proposed equation. The spacing and its sensitivities obtained from this new spacing formula were compared with those obtained from the Kirkham equation for multi-layered soils. The first and second moment method was applied to the analysis of parameter uncertainty involved in the subsurface drainage 125 Conclusions and Recommendations 126 system design of steady and non-steady state in homogeneous and layered soils. A n example was given to illustrate the procedure of the analysis and its application in the practical drainage design. A brief discussion on the methods of estimating the mean and the variance from limited field data of hydraulic conductivity measurements is given in Appendix B. 6.1 Conlusions 1. Sensitivity analysis can be employed to identify the dominant variables in drainage design. The results obtained can be used as guidances in the geohydrological investigations and the determination of parameter inputs in the design of a subsurface drainage system. 2. Drain spacing is very sensitive to the design midspan water table height, hydraulic conductivity, and drainage coefficient. Special care should be taken in the selection or measurement of these parameters for designing a cost- effective drainage system. 3. Drain spacing is not sensitive to the depth of the impermeable layer below drains and drain radius. Therefore high accuracy in the measurement of the depth to the impermeable layer is not desired. This is because most of the flow towards drains will take place in the region near the drains. In the selection of the drain size, the flow capacity should be the major factor to be considered. From the point of view of saving material, the drain size should be as small as possible while meets the requirement of flow capacity. 4. In the case of transient state drainage design, the sensitivities of spacing to the hydraulic conductivity, the depth of the impermeable layer below the drains, and the drain radius are very similar to those of steady state drainage theory. Conclusions and Recommendations 127 5. For transient state drainage design, the spacing is very sensitive to the initial water table height and the water table height after falling for a given period of time, particularly when these two positions are very close to each other. 6. Transient drainage design (falling water table criterion) can be reduced to the . simpler steady state drainage design by appropriately adjusting the R and h values. Certain combinations of R-h values in steady state design will result in close or the same spacing as obtained from transient state drainage design methods. 7. It is found that that when the water table fluctuation is small, steady state drainage theory , such as the Hooghoudt equation is recommended to use instead of transient drainage theory, such as the van Schilfgaarde equation. In these situations, steady state spacing formulae tend to provide more stable and therefore more reliable results than transient state spacing equations. 8. If transient drainage theory is used, the initial and final water table positions should be carefully selected, otherwise the drain system may be over or under designed. 9. Drain spacing is usually more sensitive to the hydraulic conductivity of the soil below the drains than that above the drains. Therefore, the hydraulic conductivity should be measured from soils deeper than the drain depth. 10. The Hooghoudt equation and the generalized Hooghoudt-Ernst equation give spacings and sensitivities which are very close to each other for different soil layering conditions. 11. A new spacing formula is developed for the situations where there are two soil layers below the drains, by combining the approach and the form of the Hooghoudt equation and the radial flow component of the Ernst equation. Conclusions and Recommendations 128 12. This new spacing formula yields spacings close to that obtained from the Kirkham equation of potential theory when the hydraulic conductivity of the soil above the drains is relatively small, and gives more accurate results when there is a highly pervious layer above the drains. In this case, use of the new spacing formula is recommended. 13. Drain spacing is not sensitive to the hydraulic conductivity of the lower layer if there are two soil layers below the drain level. Therefore it is not desirable to stratify the soil below drains into more than two layers. 14. First and second moment analysis can be used to incorporate the parameter uncertainty into the design of a subsurface drainage system. Effects of uncertainty in each parameter on the performance of a drainage system can be estimated from the first and second moment analysis. 15. The actual mean water table height is usually different from the design water table height due to the variability of climate and soil parameters, if the design is based on the conventional deterministic framework. 16. Simple equations and graphs are prepared to encourage the use of uncertainty analysis among the community of practicing agricultural engineers and designers. 17. Drain spacing can be adjusted according to the risk level the farmer is willing to accept and the funds available for the drainage system. The probabilities that the water table would exceed certain heights can be estimated from the uncertainty analysis. 6.2 Recommendations 1. The correlationships between various drainage parameters should be studied to Conclusions and Recommendations 129 increase the accuracy of sensitivity and uncertainty analysis. 2. The methods of estimating the means and the variances from limited historic data or field measurements should be studied in more detail. The accuracy of these estimates directly affects the design of a drainage system. 3. Other sources of uncertainty such as the effects of unsaturated flows in the region above water table should be studied and included in the analysis. 4. The correctness and limitation of the first and second moment analysis should be tested by experiments or studies of existing drainage systems. The recharge and soil properties should be better determined and other sources of uncertainty be eliminated or reduced. If the coefficients of variation of parameters are relatively small, the first and second moment analysis should provide very good approximation. 5. The effects of high water table on the crop yield/growth should be studied in more detail to relate the uncertainty in water table position with the uncertainty in the monetary return from agricultural drainage system so that farmers can have more direct criteria to accept or reject a drainage system design. Bibliography [1] Bakr, A . A . 1976. Stochastic analysis of the effects of spatial variations in hydraulic conductivity on groundwater flow. Thesis submitted to New Mexico Institute of Mining and Technology, Socorro, New Mexico, in partial fulfilment of the requirements for the degree of Doctor of Philosophy. [2] Bakr, A . A . , L. W . Gelhar, A. L. Gutjahr, and J . R. MacMillan. 1978. Stochastic analysis of spatial variability in subsurface flows. 1. Comparison of one-dimensional flows. Water Resour. Res. 14(2):263-271. [3] Benjamin, J . R., and C. A . Cornell. 1970. Probability, statistics and decision analysis for civil engineers. McGraw-Hill, New York. [4] Boussinesq, J . 1903. Sur le debit, en temps de secheresse, d'une source alimentee par une nappe d'eaux d'infiltration. Compt. Rend. Sceances Acad. Sci. 136:1511-1517. [5] Bouwer, H . 1974. Developing drainage design criteria. In: Jan. van Schilfgaarde (Ed.) Drainage for Agriculture. Agronomy Monograph 17:67-80. Amer. Soc. Agronomy. Madison, Wis. [6] Bouwer, H . 1969. Planning and interpreting soil permeability measurements. J . Irrig. & Drainage Div., A S C E 95(IR3):391-402. [7] Bouwer, H . , and R. D. Jackson. 1974. Determining soil properties. In: Jan. van Schilfgaarde (Ed.) Drainage for Agriculture. Agronomy Monograph 17:611-666. Amer. Soc. Agronomy. Madison, Wis. 130 Bibliography 131 [8] Bouwer, H., and J . van Schilfgaarde. 1963. Simplified method of predicting fall of water table in drained land. Trans, of A S A E . 6(4):288- 291,296. [9] Bulavko, A. G. 1971. The hydrology of marshes and marsh-ridden lands. U N E S C O . Nature and Resources, Vol. III(l). [10] Chieng, S. T. 1975. The effects of subsurface drain depths and drainage rates on water table levels. Thesis submitted to the McGill University, Montreal, Canada, in partial fulfilment of the requirements for the degree of Master of Science. [11] Chieng, S. T. 1986. Class notes for the course: Drainage engineering (BIOE 461), Bio-resource Engineering Department, University of British Columbia, Vancouver, B. B., Canada. [12] Chieng, S. T., R. S. Broughton, and N . Foroud. 1978. Drainage rates and water table depths. J . Irrig. & Drainage Div., A S C E 104(IR4):413- 433. [13] Cope, F., and E . S. Trickett. 1965. Measuring soil Moisture. Soils Fert. 28:201-208. [14] Cordova, J . R., and R. L. Bras. 1981. Physically based probability models of infiltration, soil moisture, and actual evapotranspiration. Water Resour. Res. 17(1):93-106. [15] Cornell, C. A . 1972. First order analysis of the model and parameter uncertainty. In: International Symposium on Uncertainties in Hydrologic and Water Resources Systems, University of Arizona, Tucson, Arizona. [16] Dagan, G. 1979. Models of groundwater flow in statistically homogeneous porous formulations. Water Resour. Res. 15(2):47-63. [17] Delhomme, J . P. 1978. l(5):251-266. Krigging in hydrosciences. Advances in Water Resources Bibliography 132 [18] Dettinger, M . D., and J. L. Wilson. 1981. First order analysis of uncertainty in numerical models of groundwater flow. 1. Mathematical development. Water REsour. Res. 17(1):149-161. [19] De Vries, J . 1981 and 1982. Personal communications. Soil Science Department, University of British Columbia, Vancouver, B. C , Canada. [20] Dumm, L . D. 1954. New formula for determining depth and spacing of subsurface drains in irrigated lands. Agr. Eng. 35:726-730. [21] Dumm, L . D. 1964. Transient-flow concept in subsurface drainage: its validity and use. Trans, of A S A E 17(2):142-145,151. [22] Dumm, L . D. 1968. Subsurface drainage by transient flow theory. J. Irrig. & Drainage Div., A S C E 94(IR4):505-519. [23] Dupuit, J . 1863. Etudes Theoretiques et Pratiques sur le Mouvement des Eaux, 2nd ed. Dunod, Paris. 304 p. [24] Ernst, L . F. 1956. Calculation of the steady flow of groundwater in vertical crosssections.Neth. J . Agr. Sci. 4:126-131. [25] Ernst, L . F. 1962. Groundwater flow in the saturated zone and its calculation parallel horizontal open conduits are present. Versl. Landb. Ond. 67.15:189 p. [26] Fogel, M . M . , L . Duckstein, and W. T. O'Brien. 1979. Optional design and operation of drainage systems. Paper No. 79-2071 presented at the summer meeting of the A S A E and C S A E at the University of Manitoba, Winnipeg, Manitoba, Canada. [27] Forchheimer, Ph. 1930. Hydraulik 3rd. ed. B. G. Teubner, Leipzig and Berlin. 595 P- , Bibliography 133 [28] Freeze, R. A . 1975. A Stochastic-conceptual analysis of one- dimensional groundwater flow in nonuniform homogeneous media. Water Resour. Res. 11(5):725—741. [29] Freeze, R. A., and J. A. Cherry. 1979. Groundwater. Prentice- Hall, Inc. Englewood Cliffs, New Jersey. [30] French, B . E., and J . R. O'Callaghan. 1966. A field test of drain spacing equations for agricultural land. J. Agr. Eng. Res. 11:282-295. [31] Gardner, W . H. 1965. Water Content. In: C. A . Black (Ed.) Methods of soil analysis. Agronomy 9:8-127. Amer. Soc. Agron., Madison, Wis. [32] Hooghoudt, S. B. 1940. Bijdragen tot de kennis van eenige natuurkundige grootheden van den ground, 7, Algemeene beschouwing van het problem van de detail ontwatering en de infiltratie door meddel van parallel loopende drains, greppels, slooten en kanalen. Verslag. Landbouwk. Onderzoek 46:515-707. [33] Krayenhoff van de Leur, D. A. 1958. A study of nonsteady groundwater flow with special reference to a reservoir coefficient. De Ingenieur 70 B:87-94. [34] Larsen, R. J., and M . L. Marx. 1981. A n introduction to mathematical statistics and its applications. Prentice-Hall, Inc. Englewood Cliffs, New Jersey. [35] Luthin, J . N . 1978. Drainage Engineering. Robert E . Grieger Publishing Company, Huntington, New York. [36] Maasland, M . 1959. Water table fluctuations induced by intermittent recharge. J . Geophys. Res. 64:549-559. [37] Massmann, J . W. 1987. Groundwater contamination from waste- management sites: the interaction between risk-based engineering design and regulatory policy. Thesis Bibliography 134 submitted to the University of British Columbia, Vancouver, Canada, in partial fulfilment of the requirements of the degree of Doctor of philosophy. [38] Miller, I., and J. E. Freund. 1977. Probability and statistics for engineers. PrenticeHall, Inc., Englewood, Cliffs, New Jersey. [39] Moody, W . T. 1966. Nonlinear differential equation of drain spacing. A S C E , proc. 92(IR2):l-9. [40] Mueller, Guenther. 1967. Der Felddurschlaessigkeitsbeiwert Kj und Seine Anwendung bie der Projektierung von Draenanlagen, Teil II. Z. Landeskultur. 803:153-170. [41] Musy, A . , and L. Duckstein. 1976. Bayesian approach to tile drain design. J . Irrig. & Drainage Div., A S C E 102(IR3):317-334. [42] Oosterbaan, R. J., and G . P. wind. 1979. Design and research. In: Jan Wesseling (Ed.) Proceedings of the International Drainage Workshop. ILRI, Wageningen, The Netherlands. [43] Prasher, S. 0. 1982. Examination of the design procedures for drainage/sub- irrigation systems in the Lower Fraser Valley, British Columbia. Thesis submitted to the University of British Columbia, Vancouver, Canada, in partial fulfilment of the requirements for the degree of Doctor of Philosophy. [44] Prasher, S. 0., S. 0 . Russell, T. H . Podmore, and S. T. Chieng. 1984. First and second order analysis of parameter uncertainty in drainage design. A S A E Paper No. 82-2564. [45] Prasher, S. 0., S. 0. Russell, S. T. Chieng, and T. H . Podmore. 1985. A probabilistic approach to agricultural drainage design. Trans, of A S A E 28(5):1494-1498. Bibliography 135 [46] Reeve, R. C , and N . R. Fausey. 1974. Drainage and timeliness of farming operations. In: Jan van Schilfgaarde (ED.) Drainage for agriculture. Agronomy Monograph 17:67-80. Amer. Soc. Agron., Madison, Wis. [47] Ross, S. M . 1987. Introduction to probability and statistics for engineers and scientists. John Wiley & Sons, Inc., New York. [48] Sagar, B. 1978a. Analysis of dynamic aquifers with stochastic forcing function. Water Resour. Res. 14(2):207-216. [49] Sagar, B. 1978b. Galerkin finite element procedure for analyzing flow through random media. Water Resour. Res. 14(6):1035-1044. [50] Sagar, B. 1979. Solution of linearized Boussinesq equation with stochastic boundaries and recharge. Water Resour. Res. 15(3):618-624. [51] Sagar, B . , and C. C. Kisiel. 1972. Limits of deterministic predictability of satu- rated flow equations. In: Proceedings of the second Symposium on Fundamentals of Transport Phenomenain Porous Media, Vol. 1:194- 205. International Association of Hydraulic Research, Guelph, Ontario, Canada. [52] Sagar, B., and C. Preller, 1980. Stochastic model of water table under drainage. J. Irrig. & Drainage Div., A S C E 106(IR3):189-202. [53] Skaggs, R. W . 1975. Drawdown solutions for simultaneous drainage and E T . J. Irrig. & Drainage Div., A S C E 91(IR3):9-22. [54]. Skaggs, R. W . 1978. A water management model for shallow watertable soils. Technical report No. 134, Water Resour. Res. Institute of the University of North Carolina, N . C , N . C. State University, Raleigh, N . C. Bibliography 136 [55] Smedema, L. K., and D. W. Rycroft. 1983. Land Drainage. Batsford Avademic and Educational Ltd., London. [56] Smith, L., and R. A . Freeze. 1979. Stochastic analysis of steady state flow in a bounded domain. 1. one-dimensional simulations. Water Resour. Res. 15(3):521528. [57] Soong, T. T. 1981. Probabilistic modeling and analysis in science and engineering. John Wiley & Sons, Inc., New York. [58] Spiegel, M . R. 1968. Mathematical Handbook. McGraw-Hill Book company. New York. [59] Talsma, T., and H . C. Haskew. 1959. Investigation of water- table response to tile drains in comparison with theory. J . Geophys. Res. 64:1933-1944. [60] Tezidis, G. 1968. Discussion: Falling water table between tile drains. ASCE 94(IR1):159. [61] Toksoz, S., and D. Kirkham. 1961. Graphical solution and interpretation of a new drain-spacing formula. J . Geophys. Res. 66(2):509-516. [62] Toksoz, S., and D. Kirkham. 1971a. Steady drainage of layered soils: I. Theory. J . Irrig. & Drainage Div., A S C E , proc. 97(IR1):1-18. [63] Toksoz, S., and D. Kirkham. 1971b. Steady drainage of layered soils: II. Nomo- graphs. J . Irrig. & Drainage Div., A S C E , proc. 97(IRl):19-37. [64] Uhl, V . W., and S. T. Sullivan. 1972. Uncertainty analysis in the appraisal of capital investment project. In: V . W . Uhl & W . E. Lowhian (Ed.) Uncertainty analysis for engineers. Amer. Inst, of Chem. Eng. Vol. 78:10-22. Bibliography 137 [65] van Beers, W. F. J . 1965. Some nomographs for the calculation of drain spacing S. Bull. No. 8. ILRI. Wageningen, the Netherlands. [66] van Beers, W . F. J. 1976. Computing drain spacing. Bull. No 15. ILRI, Wageningen, the Netherlands. [67] van Schilfgaarde, J . 1963. Design of tile drainage for falling water tables. A S C E , proc. 89(IR2):1-11. [68] van Schilfgaarde, J . 1964. Closing discussion: Design of tile drainage for falling water tables. A S C E , proc. 90(IR3):71-73. [69] van Schilfgaarde, J . 1965. Transient design of drainage systems. A S C E , proc. 91(IR3):9-22. [70] van Schilfgaarde, J . 1970. Theory of flow to drains. Advances in Hydroscience. 6:43-106. Academic Press, Inc., New York. [71] van Schilfgaarde, J . 1974. Nonsteady flow to drains. In: Jan van Schilfgaarde (Ed.) Drainage for agriculture. Agronomy Monograph 17:245-270. Amer. Soc. Agron., Madison, Wis. [72] Warren, J. E., and H . S. Price. 1961. Flow in heterogeneous porous media. Soc. Petrol. Eng. J . 1:153-169. [73] Werner, P. W. 1957. some problems in non-artesian groundwater flow. Amer. Geophys. Union, Trans. 38(4):511-518. [74] Wesseling, J . 1964. a comparison of the steady-state drain spacing formulas of Hooghoudt and Kirkham in connection with design practice. J . Hydrol. 2:25-32. [75] WiUardson, L. S. 1982. Drainage coefficient. A S A E . Paper No. 80- 2576. Appendix A DISCUSSION O N T H E SIMPLIFICATION OF T H E d - S e RELATIONSHIP Let's consider the Hooghoudt equation for homogeneous soils. It can be written as S =^ ( W + 2d h) = F{R, K, h, d ) e e where R is the drainage coefficient, K is the hydraulic conductivity, h is the water table height above drains at mid-spacing, and d is the equivalent depth. e Equivalent depth, d , is a function of the spacing, S, the depth to the impermeable e layer, d, and the drain radius, r: d = f(S,d,r) e The partial derivative of S with respect to R, in a strict mathematic sense, is: dS__dF dR ~ dR 8F dd dS dd dS dR e + e i.e. lyo u 1 = i _ 9K§jk dd 1 e as (A.98) Similarly we have * dK a K 1 138 (A.99) Discussion on the Simplification of the d — S Relationship e dF dh dS (A.100) ^ 1 - f t ds 8F dd dd dd 9d 1 - f t dS dF dd e 139 e 8d dr (A.101) e c 8d e (A.102) as However, in section 3.1, it was assumed that d is independent of 5 to simplify the e analysis. Therefore dd P as 0 and the above equations become &S _ 8F_ dR ~ dR (A.103) dS _ dF dK ~dK (A.104) dS__dF_ dh ~-~dh (A.105) dS dd 8F dde 8d dd (A.106) dS dr dF dde dd dr (A.107) e e which were used in section 3.1. It is noted that in Equation A.98 through A.102, the denominators are all the same (1 — Q^-^g), and that Equation A.103 through A.107 are obtained by letting dde dS Discussion on the Simplification of the d — S Relationship e 140 This simplification will not change their relative order of magnitude, but only their absolute values. Therefore the relative sensitivities of spacing to these parameters will not affected by such simplification. Furthermore conclusions based on the sensitivity analysis will also not affected by such simplification. Same analysis and conclusion are valid in cases of van Schilfgaarde equation for transient state drainage and Hooghoudt equation for steady state drainage in two-layered soils. For a closer estimation of the spacing change due to changes of other parameters, a coefficient C should be multiplied to the results obtained from the simplified approach. This coefficient equals - Table A.6 is a comparison of results of sensitivities obtained from the simplified approach and the unsimplified one for homogeneous soil. The parameter values used for the computation are: Case 1: R = 0.01 m/day, K = 0.5 m/day, h = 0.6 m, d = 5.0 m, r — 0.05 m, and Case 2: R = 0.007 m/day, K.= 1.0 m/day, h = 0.6 m, d = 3.0 m, r = 0.05 m. All parameters are assumed to change 5%. Column (1) is the results obtained from the simplified approach; (2) is the results obtained from unsimplified approach; and (3) is the results directly calculated from Hooghoudt equation. It can be seen that the results in column (2) are very close to those in column (3), while there are some differences between results in column (1) and column (3). however, the relative order of sensitivities of spacing is same in all three columns, which justifies the simplification because the primary purpose of sensitivity analysis is to identify the dominant parameters. cussion on the Simplification of the d — S Relationship e Spacing (m) R K h d r (1) -2.50 2.50 2.87 0.28 0.42 Case 1 21.98 (2) (3) -3.50 -3.36 3.50 3.46 4.02 3.99 0.39 0.36 0.59 0.58 (1) -2.50 2.50 2.83 1.13 0.27 Case 2 39.34 (2) (3) -2.95 -2.84 2.95 2.92 3.34 3.31 1.33 1.30 0.32 0.32 Table A.6: Comparison of sensitivities obtained from different approaches. Appendix B ESTIMATION OF MEANS A N D VARIANCES OF R A N D O M PARAMETERS In chapter 5, the means and variances of parameters to be used in the uncertainty analysis are assumed already known. However, how to estimate these values from limited data still remains a problem to be solved. While the true mean and variance of a parameter is usually unknown, a reasonable estimation can be obtained from statistical methods with given information about the parameter, although the information is usually limited. In this appendix, some methods to estimate mean and variance of a random parameter will be discussed, and then used in the estimation of the mean and variance of soil hydraulic conductivity. B.l Methods to Estimate Mean and Variance of Random Variable A best estimate of a parameter (mean or variance) should have the properties of unbiasedness, minimum variance, consistency, and sufficiency. The objective in parameter estimation is to determine a statistic e= h(x ,x ,...,x ) 1 2 n which gives a good estimate of the parameter 9. A statistic 0 is said to be an unbiased estimator for 6 if 'E\Q] = e 142 Estimation of Means and Variances of Random Parameters 143 for all 9. A statistic is said to be an efficient estimator for 9 if the variance of the statistic is minimum among all unbiased estimators. A statistic is said to be a consistent estimator for 9 if, as sample size n increases: lim P[\Q - 9\ > e] = 0 for all e > 0. The definition for sufficiency states that if 0 is a sufficient statistic for 9, all the sample information concerning 9 is contained in 0. Exact definitions for these criteria can be found in most statistics books (Miller and Freund, 1977; Larsen and Marx, 1981; Soong, 1981; Ross, 1987). B.l.l T h e M e t h o d of M o m e n t s The method of moments is simple in concept. Consider a selected probability density function f(x; 6\, .. . ,# ) whose parameters 9j,j = 1, 2 , . . . , m , are to be estimated based m on the sample X , ... ,X x n of X. Theoretical or population moments of X are oo / x f(x;9 ,...,9 )dx i 1 m 1 = 1,2,... (B.108) — OO On the other hand, sample moments of various orders can be found from the sample by: M --Y X) i t 7 1 » = 1,2,-.'.. (B.109) j=l The method.of moments suggests that, in order to determine the estimators 0 02, . . . , Q m 1 ? from the sample, we equate a sufficient number of the sample moments to the corresponding population equations (lower-order moment equations are preferred). Estimators can be obtained from solving the resulting moment equations oti = Mi i = 1, 2 , . . . , m Estimation of Means and Variances of Random Parameters 144 The method of moments is straitforward, and the moment equations are seldom difficult to solve. However, this method can not guarantee such desirable properties as unbiasedness and minimum variance (efficiency) for estimators so obtained. B.1.2 The Method of Maximum Likelihood Let f(x; 8 ) be the probability density function of population X where 9\,. .. ,9 m are the parameters to be estimated from a set of sample values x-i,x ,... 2 density function of the corresponding sample X ,X ,... 1 ,X 2 ,x . n m The joint has the form n f{x ;9 ,...,9 )f{x ;9 ,...,9 )---f{x ;9 ,...,9 ) l l m 2 l m n 1 m The likelihood function L of a set of n sample values from the population is defined by n L(x ,...,x ;e ,...,e ) n 1 1 = Y[f(x ;6 ,...,6 ) m i 1 m (B.110) i=l When the sample are given, the likelihood function L becomes a function of variables 9i,...,8 . m The concept of the method of maximum likelihood consists of choosing, as estimates of 9j,j = 1,2, ... ,m, the values of 9j,j = 1,2,... ,m, that maximizes the likelihood function L. Hence, the maximum likelihood estimates 9j of 9j,j = 1,2,... based on sample values x x ,... ly ,x 2 n can be determined from 8L —- = 0 j = 1,2,...,m 89 or — ^ = 0 .89 r=l,2,...,rn ,m, Estimation B.l.3 of Means and Variances of Random Parameters 145 Bayesian Estimation of Mean Bayesian methods of estimation have received impetus and much wider applicability in recent years through the concept of personal, or subjective, probability. The idea of Bayesian methods is to regard an unknown parameter 9 to be estimated being the value of a random variable from a given probability distribution. This usually arises when, prior to the observance of the outcomes of the data X\, X , • • • , X , we have 2 n some information about the value of 9 and this information is expressible in terms of a probability distribution, which is referred as the prior distribution of 9. Suppose that based on the knowledge and experiences, we have prior feelings about 9 that it can b regarded as being the value of a continuous random variable having probability density function /(#); and suppose that we obtained a set of sample values Xi = Xi, i = 1, 2, .. ., n, then the updated, or conditional, probability density function of 9 is as follows: }[V\x .. . ,x ) u n \D\tta\jD -—j--tt- (B.lll) J - o o / ' ^ i i • • • ,x \9)f{9)d9 n where the conditional density function f(9\x ,...,x ) x is called the posterior n density function; f(x\, . .. , x \9) is the likelihood that data values Xi = Xi,i = 1,2,...,n, when n 9 is the value of the parameter; and f(9) is the prior probability density function. The best estimate of 9, given the data values x\, x , •.. , x , is the mean of the posterior 2 n distribution f(9\xi, ... ,x ), which is called the Bayes estimator: n E oo / -oo 9f{B\x ,...,x )d9 1 As an example, let's consider a sample X\,... ,X n (B.112) n which are independent normal random variables, each having unknown mean fi and known variance a , and fi itself is a 2 random variable which is normally distributed, with known mean / i and known a\ (prior 0 Estimation of Means and Variances of Random Parameters 146 information). In order to determine the Bayes estimator, we need first to determine the conditional density function of /z given the values of Xi,i = 1, . .. ,n: f(x-L, .. . ,X \fJ,)f(/J,) f'([i\xi, ...,x ) n n ...,x ) /(aJi, n where f(x ,...,x \/i) 1 = n l EL, (-=-) e n V27Tcr 1 (t*-no) 2 /(/z) = -?=-e 27TO" -o 2 and oo / /(asi, . . . ,x \fi)f(fi)dfi n — OO Working out the above various terms, we obtained the Bayes estimator: (H = E\ji\X u ...,X }= ° X + n a n 2 2 f fi 2 0 (B.113) which is a weighted average of X, the sample mean, and fi , the a priori mean. The 0 variance of population a can be estimated from the sample variance. The variance of 2 the population mean can be calculated by aa \| nal + a 2 B.1.4 2 2 (B.114) Confidence Interval for Mean and Variance in Normal Distribution The variance of the population a is usually not known. However, it can be reasonably 2 estimated by the sample variance S . 2 For a sample size n which is large enough to approximate a with S, the random variable Y = X ~ fi is approximately a unit normal distribution. The interval Estimation of Means and Variances of Random Parameters X - Z a / 2 4 = < Ii < X + Z a / 2 147 ~ (B.115) is an approximate confidence interval for /x having the degree of confidence (1 — a). In the above equation X is the sample mean, S is the sample standard deviation. If the sample size n is not large enough, it can be proved that random variable y _ ______ has a ^-distribution with (n — 1) degrees of freedom. Therefore the interval . g g X - * _ i , / 2 - 7 = < li < X + t _ i n is a confidence interval for a n i a / 2 -= (B.116) having the degree of confidence (1 — a). As mentioned previously, an unbiased point estimator for the population variance a 2 is the sample variance S . For the construction of confidence interval for a , let's 2 2 introduce the random variable a- 2 It can be proved that D has a % distribution with (n — 1) degrees of freedom. The 2 interval <t 1 ) 3 2 X -l,a/2 n < < < r 1 ) s ' ^n-l,l-a/2 is a confidence interval for <J having the degree of confidence (1 — a). 2 B.2 E s t i m a t i o n of the M e a n and Variance of H y d r a u l i c C o n d u c t i v i t y Hydraulic conductivity is defined by Darcy's law: (B.i") Estimation of Means and Variances of Random Parameters 148 where A = area perpendicular to flow direction, Q = flux through area A, ^ = hydraulic gradient. From the above definition, hydraulic conductivity is not a point value in a strict sense. It is an average value over a volume of soil (A»dx). However, in many agricultural investigations the size of the soil sample used to measure the hydraulic conductivity is so much smaller than the scale of the flow field that the sampled value is essentially a point value, and the hydraulic conductivity can be treated as continuously distributed over the field. Strictly speaking, hydraulic conductivities do not vary in space in a purely random, unstructured manner. Neighboring values are rather correlated. The scale of correla- tion can be represented by a parameter called correlation length, which is defined as the distance over which the correlation is positive. However, the distance between two measurements in agricultural investigation is usually too large to consider any meaningful correlation between the measured values of hydraulic conductivity. That is to say that the sample values of hydraulic conductivity can be considered as independent. There is a fairly extensive literature that suggests that hydraulic conductivities are often lognormally distributed (Massmann, 1987). Suppose we have a sample of size n with n observed values, K\,... , K , then ln K n i} ..., lnK n are normally distributed. Let Xj = ln Kj . j = 1, . .. , 7i and X be the sample mean, S the sample standard deviation, the population mean (to be estimated) and o ^ the population variance. - 2 Now let's estimate the population mean fii^K by the method of moments. Xj,j 1, .. . , n are random variables from a normal distribution with mean of K a n = d variance of hiK> a which c a n be estimated from the sample variance S . 2 The probability density Estimation of Means and Variances of Random Parameters 149 function with unknown ninK is (. -l*ln x / ( s j / i i n j r ) = -j=^- K ) 2 ( and the first moment of the distribution oo too Xfixifiin^dx = / — oo fl^K On the other hand, the first sample moment is 1 71 Equating these two moments gives ^lntf = - X > j i=l (B.118) n The above result can also be obtained from the method of maximum likelihood. The likelihood function L is S" (' j-'MnK-) ! =1 = II /(^j^inic) = (-7= L{x .. lt ,•=1 The first derivative of L to dh = fi\nK 1 V^TTO-lnif is ) n E l ( a j - Wntf)^ j = 2 In i f Let dL d/i in if- = 0 then or 1 71 Minis: = - 1] asj )" 2 150 Estimation of Means and Variances of Random Parameters Recall that Xj = In Kj, (!_ K is the mean of the logarithm of hydraulic conductivity. Because of the symmetry of the normal distribution, variable X = In i f . Let K m is also the median of random fi]nK be the median of the random hydraulic conductivity, then 1 - V In Kj n fi_ = In K K m = = -\nK K ---K t n = 2 n \n^K K --K t r n Therefore the median of hydraulic conductivity is K = ^JK K ---K m X 2 n (B.119) By definition of median, there is a 50% probability that the value of hydraulic conductivity will less or larger than K , respectively. m Bouwer (1969) showed that K m (which is called geometric mean) gives the best esti- mate of hydraulic mean value of K in the studies with a resistance network analog using various randomized if— distributions in a two-dimensional system. The variance of hydraulic conductivity can be approximately calculated by: Var(tf) = ^U( i- rn) K It should be noted that K m K 2 ( B 1 2 Q ) and Var(iif) are not the mean and variance of hydraulic conductivity in a statistical sense, but rather a practical interpretation of some measurement values. From the probability density function of lognormal random variable, the mean and variance of K are p K =K m e ^ <r =&(e'**-l) 2 K (B.121) (B-122) Estimation of Means and Variances of Random Parameters 151 where o"^^- is the variance of logarithmic hydraulic conductivity, which can be estimated by: 2 rr 2 ~<? - 2 71 — i or = *" \ n — 1 = 1 k (B.123) J As an illustration, let's assume we have 8 sample measurements of hydraulic conductivity over a field under investigation. These sample values are: K\ — 0.4, K — 2 0.6, K 3 = 0.3, = 0.7, i f = 0.4, # 5 6 = 0.5, K 7 = 0.8, K 8 = 0.6 (all units are m/day). Suppose these data values are lognormally distributed. The median value of hydraulic conductivity will be K m = 0.51 m/day. From equation B.123, <r \ = 0.1079, and from l lK equations B.121 and B.122, we have = 0.54 m/day, \LK <J = 0.033 (m/day) . On the 2 2 K other hand, the variance calculated from equation B.120 is Var(if)=0.0293 (m/day) . 2 The differences between K m and fiR, u 2 K and Var(if), are fairly small, thus K m and Var(if) can be reasonably used as the mean and variance of hydraulic conductivity in the uncertainty analysis. B.3 Application of Bayes Method in the Estimation of the Mean Hydraulic Conductivity In this section, an example is given to illustrate how the Bayes method of estimation can be used in drainage design. For the purpose of simplicity, we assume that the hydraulic conductivity has a normal distribution. The mean hydraulic conductivity is unknown and can be treated as random variable which is also normally distributed. Based on the knowledge and experience on the field under investigation, the engineer feels that the hydraulic conductivity of this field can best be described by a normal distribution with a mean of /x =0-6 m/day and standard deviation of cr =0.08 m/day (prior information). o 0 Estimation of Means and Variances of Random Parameters 152 After investigation, a number of, say 8, random sample values have been obtained from measurements. The sample mean X=0.54 m/day, the sample standard deviation S = 0.17 m/day. The population standard deviation <x can be approximated with S. Then the updated mean of hydraulic conductivity after some measurements is 2 lii = — — - A 2 4- —2——zfMo = 0.56 m/day and the variance of the mean is 2 2 0.05 m/day The mean of the posterior distribution is based on the direct evidence as well as prior information. Bayesian method of estimation can be readily used to estimate the means of precipitation and evapotranspiration. The historic record can provide very good information to assume certain kind of distribution, and the recent events can be used to modify the prior distribution to yield more accurate one. Appendix C G R A P H OF G E O M E T R Y F A C T O R a F O R R A D I A L 153 FLOW Graph of Geometry Factor a for Radial Flow Figure C.35: Geometry factor a for radial flow (taken from van Beers, 1976). 154
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Sensitivity and uncertainty analysis of subsurface...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Sensitivity and uncertainty analysis of subsurface drainage design Wu, Guangxi 1988
pdf
Page Metadata
Item Metadata
Title | Sensitivity and uncertainty analysis of subsurface drainage design |
Creator |
Wu, Guangxi |
Publisher | University of British Columbia |
Date Issued | 1988 |
Description | Literature on subsurface drainage theories, determination of drainage parameters, and analysis approaches of uncertainty was reviewed. Sensitivity analysis was carried out on drain spacing equations for steady state and nonsteady state, in homogeneous soils and in layered soils. It was found that drain spacing is very sensitive to the hydraulic conductivity, the drainage coefficient, and the design midspan water table height. Spacing is not sensitive to the depth of the impermeable layer and the drain radius. In transient state, spacing is extremely sensitive to the midspan water table heights if the water table fall is relatively small. In that case steady state theory will yield more reliable results and its use is recommended. Drain spacing is usually more sensitive to the hydraulic conductivity of the soil below the drains than to that of the soil above the drains. Therefore, it is desirable to take samples from deeper soil when measuring hydraulic conductivity. A new spacing formula was developed for two-layered soils and a special case of three-layered soils with drains at the interface of the top two layers. This equation was compared with the Kirkham equation. The new formula yields spacings close to the Kirkham equation if the hydraulic conductivity of the soil above the drains is relatively small; otherwise, it tends to give more accurate results. First and second order analysis methods were employed to analyze parameter uncertainty in subsurface drainage design. It was found that conventional design methods based on a deterministic framework may result in inadequate spacing due to the uncertainty involved. Uncertainty may be incorporated into practical design by using the simple equations and graphs presented in this research; the procedure was illustrated through an example. Conclusions were drawn from the present study and recommendations were made for future research. |
Subject |
Subsurface drainage -- Mathematical models Drainage -- Mathematical models |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-09-16 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0097918 |
URI | http://hdl.handle.net/2429/28529 |
Degree |
Master of Applied Science - MASc |
Program |
Biomedical Engineering |
Affiliation |
Applied Science, Faculty of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1988_A7 W82.pdf [ 7.09MB ]
- Metadata
- JSON: 831-1.0097918.json
- JSON-LD: 831-1.0097918-ld.json
- RDF/XML (Pretty): 831-1.0097918-rdf.xml
- RDF/JSON: 831-1.0097918-rdf.json
- Turtle: 831-1.0097918-turtle.txt
- N-Triples: 831-1.0097918-rdf-ntriples.txt
- Original Record: 831-1.0097918-source.json
- Full Text
- 831-1.0097918-fulltext.txt
- Citation
- 831-1.0097918.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0097918/manifest