.5, #3- E3?! ‘ .a ~ "' a- ‘ ,1; ‘ ‘. H“, «I ‘ fi ‘1. ‘ n»‘ I“,." 97“" 3&3)“ Al . 'J‘ 5333' _':"x‘. "\lf f gig-fix “a . I l tf ' ’ ‘\ ‘. 35:31.2“ ‘ . s- . {"4 w ‘l G 2 'V‘ ' II M bug: ;) *f-l Ii -‘.‘ .niV, n .1" ’ ‘ 'Q‘ Ia“ ““4? “‘2‘0'3 "y $91;- .131 7. jw‘ I“: 'r‘: é ‘fi‘! 1' 7“ 'SJ-h .5» mg 0 ICHtGAN STAT v as ml l 1|le illmilllllllllll 3 1 93 01399 2130 This is to certify that the dissertation entitled Use of Scaling Theory in Predictive Methods to Efficiently Construct Two- and Three—Phase Constitutive Relationships presented by Herr Soeryantono has been accepted towards fulfillment of the requirements for Ph . D . degree in Civil Engineering Date April, 1995 MSU is an Affirmative Action/Equal Opportunity Institution 0-12771 t.“ fig”..- __ v LlBRARY Michigan $tate Unlverslty PLACE II RETURN BOXto mallow: momma your mead. TO AVOID FINES Mum on or “on dd. duo. Use of Scaling Theory in Predictive Methods to Efficiently Construct Two- and Three-Phase Constitutive Relationships by Herr Soeryantono A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Civil and Environmental Engineering 1995 Roger B. Wallace ABSTRACT Use of Scaling Theory in Predictive Methods to Efl‘iciently Construct Two- and Three-Phase Constitutive Relationships by Herr Soeryantono Scaling theory (Parker et al., 1987) combined with Leverette’s (1941) predictive method permits prediction of the relative permeability-saturation-pressure (kt-S-P) relationship of any two- or three-phase system from P-S data of a particular two-phase system and interfacial tension data. Implementation of the theory is studied by comparing the two- and three-phase kr-S-P curves constructed with Leverette’s predictive method with and without the use of scaling theory. The observations include predictions which are based on air-water or air-oil P-S data using two soils: uniform graded Ottawa sand and loamy sand Metea soil. Comparison of the kr-S—P curves obtained with and without the scaling theory showed that, at least for the investigated cases, the use of air-oil system as the basis of prediction is only suggested for a well-graded soil. When involving a uniform grain type soil the use of air-water system as the basis of prediction is preferable. The required P-S data were measured in a pressure cell using the pressure equilibrium method. This method reverses the conventional technique (Klute, 1986) by imposing rapid finite desaturation and subsequently allowing hydrostatic equilibrium under a zero flux boundary conditions. Theoretical verification showed that this method is more efficient than the conventional technique. Laboratory experiments confirmed the theoretical result and showed that the P-S relation obtained from the data, statistically, indistinguishable from that obtained by conventional method. To Ibu, Christina, Titon and Nino iv ACKNOWLEDGMENT I gratefully acknowledge the contributions of those who have made this dissertation possible. My deepest gratitude goes to Dr. Roger B. Wallace, Associate Professor, Department of Civil and Environmental Engineering, Michigan State University, my dissertation director and advisor during my course work. I appreciate his fatherly advice, patience, and continual encouragement throughout my entire graduate program and his stamina in reviewing a great many drafis. He was instrumental to the successful completion of this dissertation. I would also like to thank Dr. Thomas C. Voice, Associate Professor Department of Civil and Environmental Engineering; Dr. David C. Wiggert, Professor, Department of Civil and Environmental Engineering; and Dr. Raymond J. Kunze, Professor, Department of Crops and Soil Science, Michigan State University for serving in my dissertation committee and for their extremely valuable comments in this study. For the further development of pressure saturation measurement apparatus which was essential in this study I owe Dr. Arthur T. Corey, Profesor Emeritus, Department of Agricultural Engineering, Colorado State University. Dr. Mokma, Professor, Department of Crop and Soil Sciences helped to select the soil types I used. Research funding was obtained from several sources. I wish to thank the Second Higher Education Development Project, Ministry of Education and Culture Republic of Indonesia that provided support for the first three years of my doctoral program. This research could not have been completed without support from the National Institute for ' Environmental Health Sciences (Grant ESO4911) and the Michigan State University Institute for Environmental Toxicology. Ruth Dukesbury and Kahlil Rowter helped in proofreading. For that I am very grateful. To Mike Annable, Lizzette Chevalier, Reza Rakhshandehroo, Munjed Maraqa, Michel Hsu, Yao J .1, Linda Steinmann and Brenda Minott, my colleagues in the Department of Civil and Environmental Engineering, Michigan State University, I would like to express my deep appreciation to for all their support and encouragement Finally, I want to thank my wife and children: Christina, Titon and Nino for their love, care, support and understanding. vi TABLE OF CONTENT LIST OF TABLES ....................................................................................................................................... ix LIST OF FIGURES ....................................................................................................................................... x LIST OF SYMBOLS - - __ - ..................................................................................................... xii CHAPTER I: INTRODUCTION ................................................................................................................. 1 1.1. BACKGROUND .............. .................................................................. I 1.2. GENERAL OBJECTIVES, APPROACHES AND PRESENTATION ORGANIZATION ............................................. 7 1.3. DEFINITIONS AND GENERAL OVERVIEW. - ................................................ 9 1.3.1. Multiphase System ................................................................................................... 9 1.3.2. Multiphase Flow Equation and the Predictive Method ................................................................... 11 1.3.3. Methods to Scale the Pressure-Saturation Relationship .................................................................. 14 1.3.3.1. General Formulation of the Pressure-Saturation Scaling ........................................................................ 14 1.3.3.2. Leverette (1941) .1 Function . ........................................................ 18 1.3.3.3. Similar Media Concept by Miller and Miller (1956) ............................. 19 1.3.3.4. Brooks and Corey (1964) Pressure-Saturation Function as A Scaling Method ...................................... 21 1.3.3.5. Lenhard and Parker (1987) Scaling Theory. ............................................ 22 1.3.3.6. The Issue of Contact Angle. ....................................................................................... 23 1.3.3.7. Dependency of Scaling Method on Soil Texture .- ....... 24 CHAPTER II: THE ACCEPTABILITY OF THE PRESSURE EQUILIBRIUM METHOD TO MEASURE DRAINAGE PATH OF SOIL PRESSURE-SATURATION ........... 25 2.1. BACKGROUND ....................................... - -. - - ............................ 25 2.2. MEASUREMENT METHOD AND MATERIAL .............................................................................................. 29 2.3 . RESULT AND DISCUSSION ........................................................................................................................ 36 2.3.1. Pressure-Saturation Measurement Simulation... -. ......................................... 36 2.3.2. Hysteretic Issue in the P-S Measurement Employing the Pressure Equilibrium Method ............... 41 2.3.3. Pressure-Saturation Measurement employing the Pressure Equilibrium and the Saturation Equilibrium methods. ...................................................................................................................... 44 2.4. SUMMARY AND CONCLUSION .................................................................................................................. 55 CHAPTER III: THE ACCEPT ABILITY OF SCALING THEORY IMPLEMENTATION IN PREDICT ING THREE-PHASE CONSTITUTIVE EQUATIONS FROM TWO-PHASE PRESSURE-SATURATION DATA. .................................................... 58 3.1. BACKGROUND- _- _ .................................................................................................... 58 3 .2. THEORETICAL BACKGROUND. ................................................................................................................. 62 3 .3. MATERIAL AND METHOD. ....................................................................................................................... 67 vii 3.4.3.111 3.5. Sum APPENDlC A. Pressu M M :12. .\l A3 P: B. Pressu Bl A; 8.3 711 C REgreg Cl. a; C2. A1 C3, A: C4. Ai D. [men’a DI. A. DZ. 0: D3. A: Di A: 0.5. 0. D 6. A. E. Press” 51. A. El» 0: E3 A] 3.4. RESULT AND DISCUSSION ........................................................................................................................ 69 3.4.1. Comparison Across the Soil Types ................................................................................................. 83 3.4.2. Comparison Across the Methods ................................................ 85 3.4.3. The Effect of the Propagated Error on the Constitutive Variables ................................................. 91 3.5. SUMMARY AND CONCLUSION .................................................................................................................. 95 APPENDICES .................................................................................................................... 102 A. Pressure-Saturation Measurement Simulation ................................................................................... 102 A.1. Model Conceptualization. ........................................................................................ 103 A2. Model Flowchart .................................................................................................... 104 A.3. Program Listing in Excel Worksheet ..................................... 105 B. Pressure-Saturation Data Obtained Using Saturation Equilibrium Method. ...................................... 138 8.1. Air-Oil Data. .................................... 139 3.2. Air-Water Data. ...................................... . ....................................................................... 142 C. Regression Output. ...... . . . . ...................................... 145 CI. Air-Water System, Based on Data Measured Using Saturation Equilibrium Method. ............................... 146 C2. Air-Water System, Based on Data Measured Using Pressure Equilibrium Method. .. 148 C.3. Air-Oil System, Based on Data Measured Using Saturation Equilibrium Method ..................................... 152 C4. Air-Water System, Based on Data Measured Using Pressure Equilibrium Method. .................................. 145 D. Interfacial Tension Data - -- - - -- ............................................... 157 DJ. Air-Water System - Metea Soil .................................................................................................................. 158 D2. Oil-Water System - Metea Soil .................................................................................................................. 160 D3. Air-Oil System - Metea Soil ........................................................................................... 161 0.4. Air-Water System - Ottawa Sand ........ .. ................................................ 162 D5. Oil-Water System - Ottawa Sand . .............................................................................. 164 0.6. Air-Oil System - Ottawa Sand .............................................................. 165 E. Pressure-Saturation Data Obtained Using Pressure Equilibrium Method .......................................... 166 El. Air-Water Data - Metea Soil. ............................................. 167 13.2. Oil-Water Data - Metea Soil. ...................................................... 177 E.3. Air-Oil Data - Metea Soil. .......................................................................................................................... 185 8.4. Air-Water Data - Ottawa Sand. ................................................................................................................... 188 15.5. Oil-Water Data - Ottawa Sand .................................................................................................................... 198 13.6. Air-Oil Data - Ottawa Sand ................................................................................................ 208 viii Tabl Iabl Tabl- LIST OF TABLES Table 2-1: The best estimated air-water and air-oil P-S parameters based on the pooled data measured using the PE and SE methods .................... 54 Table 3-1: The obtained interfacial tensions, P-S parameters and scaling factors ....................................................................................................... 70 Table 3-2: Summary of the agreement between prediction outcomes by methods 2 or 3 to method 1 ...................................................................... 82 ix ll Figure Figur: Ilgun Figun Figure Fight Figure LIST OF FIGURES Figure 2-1: Pressure cell detail ................................................................................... 30 Figure 2-2: The apparatus set-up ................................................................................ 31 Figure 2-3: Pressure head at F015 cm and accumulated outflow with respect to time obtained from the PE and SE method simulations using hp= -50 cmHzO and hp= -37 cmHZO respectively. ............................................................................................. 39 Figure 2—4: Saturation profiles obtained from the saturation equilibrium and the pressure equilibrium simulations at initial time (heavy solid lines), during desaturation stage (light solid lines), at the end of desaturation stage (heavy dash-dotted lines), during redistribution stage (light dashed line) and at simulation end time (heavy dashed line). ................................................ 42 Figure 2-5: (a)-(d) Water pressure head changes with respect to time as recorded by the transducer during redistribution time in a measurement employing the PE method and (e) the obtained P-S data from that particular measurement superimposed with the fitted P—S function by van Genuchten 1981) .............................. 46 Figure 2-6: P—S curves obtained by the PE method (dashed lines) superimposed with those obtained by the SE method (solid lines). The confidence limits (shaded area) were based on regression upon the SE data (diamond marker). ...................................... 49 Figure 2-7: P-S data obtained from measurement using the PE method (plus marker) and the SE method (diamond marker).The P-S curves (solid lines) and the confidence limits (dashed lines) were based on regression upon the PE data (plus marker) ....................... 50 Figure 3-1: Two-phase oil-water and air-oil pressure-saturation curves obtained by method 1 (solid lines) with their 95% confidence limits (shaded area) superimposed with those obtained by method 2 (heavy dashed lines) with their uncertainty limits (light dashed lines) and the measured data. ............................................. 71 Figure 3-2: Two-phase air-water and oil-water pressure-saturation curves obtained by method 1 (solid lines) with their 95% confidence limits (shaded area) superimposed with those obtained by method 3 (heavy dashed lines) with their uncertainty limits (light dashed lines) and the measured data. ............................................. 72 9.3 2 A... 2.3:...“ w: Figure 3-3: Two-phase air-water, oil water and air-oil relative permeability curves obtained by method 1 (solid lines) with their uncertainty limits (shaded area) superimposed with those obtained by method 2 (heavy dashed lines) with their uncertainty limits (light dashed lines) and the measured data. ................ 73 Figure 3-4: Two-phase air-water, oil-water and air-oil relative permeability curves obtained by method 1 (solid lines) with their uncertainty limits (shaded area) superimposed with those obtained by method 3 (heavy dashed lines) with their uncertainty limits (light dashed lines). ..................................................... 74 Figure 3-5: Three-phase oil saturation curves obtained by method 1 (solid lines) with their uncertainty limits (shaded area) superimposed with those obtained by method 2 (heavy _ dashed lines) with their uncertainty limits (light dashed ' lines) ......................................................................................................... 75 Figure 3-6: Three-phase oil saturation curves obtained by method 1 (solid lines) with their uncertainty limits (shaded area) superimposed with those obtained by method 3 (heavy dashed lines) with their uncertainty limits (light dashed lines) ......................................................................................................... 76 Figure 3-7: Three-phase oil relative permeability obtained by method 1 (solid lines) and method 2 (dashed lines). ............................................... 77 Figure 3-8: Three-phase oil relative permeability obtained by method 1 (solid lines) with their uncertainty limits (shaded area) and those obtained by method 2 (heavy dashed lines) with their uncertainty limits (light dashed lines). ..................................................... 78 Figure 3-9: Three-phase oil relative permeability obtained by method 1 (solid lines) with their uncertainty limits (shaded area) and those obtained by method 3 (heavy dashed lines) with their uncertainty limits (light dashed lines). ..................................................... 79 Figure 3-10: Parameters Sr and n estimated by the individual regression with their 95% confidence limits ............................................................. 87 Figure 3-11: Two-phase pressure-saturation curves obtained by the pooled regression (heavy dashed lines) with their 95% confidence limits (light dashed lines) and those obtained by the individual regression (solid lines) with their 95% confidence limits (shaded area). ................................................................................. 94 xi LIST OF SYMBOLS P-S parameter of van Genuchten (1982) P-S function in two phase system pq. or of air-water P-S curve. or of oil-water P-S curve. or of air-oil P-S curve. scaling factor. air-oil scaling factor. air-water scaling factor. oil-water scaling factor. uncertainty limits of parameter or. uncertainty limits of parameter n. uncertainty limits of saturation. uncertainty limits of residual saturation. Brooks and Carey (1964) pore size index. porosity. specific density of the reference phase. specific density of phase p. interfacial tension. air-water interfacial tension. oil-water interfacial tension. air-oil interfacial tension. oil-water interfacial tension. accumulated outflow volume. target outflow volume in simulation. contact angle. covariance matrix [C]=(FT.F)'1 xii E Wgfi. 3* h- 3'- m an o 2’ O 3‘3‘3‘3‘3“ “‘2’ E ' 3‘ S m, "9%) ’1' "M as(h,-) BS(h,-) BS(h,-)] 3Sr 8a an grain size scale factor in similar media concept. jacobian matrix [F]: |: degree of fieedom. gravitational acceleration pressure head. air phase pressure head. capillary pressure head of air-water system. capillary pressure head of air-water system. oil phase pressure head. capillary pressure head of air-water system. prescribed pressure head lower boundary condition. capillary pressure head in two-phase system pq. water phase pressure head. summation index, i,j=1, 2, 3, . . . Leverette J function. saturated relative permeability. saturated permeability tensor. relative permeability. oil relative permeability in air-oil system. oil relative permeability in oil-water system. oil relative permeability in three-phase system. relative permeability of phase p. water relative permeability in three-phase system. water relative permeability in air-water system. water relative permeability in oil-water system. P-S parameter of van Genuchten (1982) P-S function in two phase system pq. P-S parameter of van Genuchten (1982) P-S function in two phase system pq. capillary pressure. xiii S58 S 80W — 30W St Saow W Saw W — 30W Sw phase type, p = air, oil, water. displacement capillary pressure. wetting-phase pressure. capillary pressure in two-phase system pq. nonwetting-phase pressure. two-phase system p==oil, water. q=air, oil. wetting phase flux. macroscopic radius of curvature. microscopic radius of curvature, radius of curvature interfaces in capillary tube. degree of saturation. degree of saturation as a function of capillary pressure in two phase system pq. saturation as a function of macroscopic radius of curvature, pore size distribution function. water saturation profile at td. water saturation profile at tc. water saturation profile at to. effective saturation S,,E(S-Sr)/(1-Sr ). oil phase saturation in air-oil system. oil phase saturation in three-phase system. degree of saturation of phase p. residual saturation, residual saturation of wetting phase p in two phase system P9- sum of squared error total liquid saturation, saturation of oil phase plus water phase in three-phase system. total liquid effective saturation in three-phase system. water phase saturation in three phase system. water phase saturation in air-water system. water phase effective saturation in three-phase system. xiv water phase saturation in oil-water system. time. time required to drain Vd, drainage time. student t value at 95% probability and degree of freedom df. time required to reach equilibrium. initial time. direction in Cartesian coordinates. vertical distance. XV CHAPTER I INTRODUCTION 1.1. BACKGROUND Volatile organic contaminants (V OCs) are extensively present in waste dumps. Some of the V005 exist in the form of non-aqueous phase liquid (NAPL). When NAPL is allowed to infiltrate the soil, some fraction will partition into and contaminate the water phase that exist in soil pores. Due to its immiscible nature with water and gas, the remaining fraction will stay as a separate liquid phase. This fraction will migrate further leaving a trail of residual NAPL in the soil pore. This residual fraction will become a long-term source of ground water contamination and further NAPL migration will extend the contaminated region. As a special case in immiscible flow problems, the coexistence of gas-NAPL- water in porous medium is referred to as a three-phase system. The phases involved in such a system are commonly referred to as air, oil and water phases with the VOC being the oil phase. According to Gee et a1. (1991) the development of a mathematical model to accurately assess the spatial and temporal quantification of these contaminants is impeded by the lack of three-phase constitutive information relating relative permeability (kt), saturation (S) and capillary pressure (P). As pointed out by Abriola (1989), the primary difficulty is in measuring the capillary pressure of an individual phase in a three-phase system. Therefore, most analyses of constitutive properties governing three-phase flow 2 (Corey, et al., 1956; Stone, 1973; Azis and Settari, 1979) are based on predictive methods which extend the information obtained from two-phase measurements. The predictive method conceptualizes the three-phase system as two separate interdependent two-phase systems; air-total liquid and oil-water. The degree of total liquid saturation (ie. the saturation of the oil and water altogether) is assumed to be dictated only by the air-oil capillary pressure and that the degree of water saturation is is dictated only by the oil-water capillary pressure. Thus, a complete set of two- and three phase kr-S-P relationships can be determined based only on the two-phase P-S information. Parker et al. (1987) proposed a scaling theory which simplifies the implementation of the predictive method even further. Within the same porous medium, the scaling theory permits prediction of the PS relationship of any two-phase system from the P-S relationship of a particular two-phase system provided that the pertaining scaling factor is known. The scaling theory determines the scaling factors from interfacial tension data. Therefore, when the P-S data of one two-phase system is available the kr-S- P information involving any fluids in any combinations can be determined. The additional data required are only the interfacial tension data. The scaling theory was incorporated in a versatile numerical three-phase flow model (Kaluarachi et al., 1990), which was further improved and has been shown valid for both desaturation and imbibition scanning paths (Parker and Lenhard, 1987b) and to accept multi component organic chemical transport (Parker and Kaluarachchi, 1990). This model is one among the few formulated to predict DNAPL movement (Huling and 3 Weaver, 1991). Incorporation of scaling theory into the model eliminates the necessity to provide air-oil and oil-water P-S information. It allows the use of air-water PS and interfacial tension data only. NAPL contamination assessment in the unsaturated zone most likely involves unconsolidated soils. Air-water P-S measurement on an unconsolidated soil sample is a common practice among soil scientist and the measurement protocol for this case is well established. Thus, air-water P-S data are commonly available. Interfacial tension measurement can also be considered standard measurements with an established protocol as described by Adamson (1976). Considering that the data requirements are less, scaling theory greatly improves the applicability of the multi-phase and multi-component model by Kaluarachi et al. (1990). Beside the advantages stated above, however, the use of scaling theory in the predictive method may contribute to additional loss in accuracy of the obtained kr-S relationship as Lenhard etal. (1987) has suggested. Scaling theory regards the soil matrix as a rigid porous medium. Accordingly, the pore size distribution is also assumed constant with respect to time. This assumption is widely adopted in the arcs of ground water engineering. However, as noted in Lenhard et al. (1987), the k,-S relationship prediction via scaling theory might be sensitive to any disturbance in the pore size distribution. Small changes in pore size distributions, such as that due to swelling, which is not clearly evidenced in the P-S plot might afi'ect the prediction of k,-S curves significantly. On the other hand, the pore size distribution is inherent in P-S relationship 4 and this relationship is represented only by the P-S parameters. Therefore, any error or bias in these parameters might also affect the prediction of kr-S curves significantly such as the effect of pore size distribution changes. Note that, kr-S-P relationship prediction via scaling theory is based on P-S parameters measured in one two-phase system only. These parameters are referred to as base parameters. The above discussion suggests the importance of a study to investigate the effect of error or bias embodied in the base parameters to the prediction of the k,-S-P curves via the scaling theory. The magnitude of error and the potential of obtaining bias parameters can be deduced from the precision of those parameters. Accuracy and precision in determining the P-S parameters might be affected by soil texture and the type of the immiscible fluid involved. A study by McCuen et al. (1981), based on 1,085 samples measured by Rawls et al. (1976) and Holtan et al. (1968), observed distinct P-S parameter variance between soil texture classes. As observed in the measurement result reported by Lenhard and Parker (1987), the P-S parameters derived from air-water data possessed higher variation than those determined from either air-oil or oil-water data . If the accuracy and precision of the P-S parameters differs from one two-phase or from one porous medium to the other, so does the performance of the prediction of the k,-S-P relationship. In turn, the scaling theory might perform differently when it is implemented on different soil texture using base parameters from different two-phase system. 5 In providing the three-phase kr-S-P relationship information the predictive method relies on the availability of two-phase P-S information. When involving unconsolidated soil samples, which is typical for subsurface soil samples, P-S measurement of the air-water system is more common than the air-oil or oil-water system. Therefore, air-water P-S data in general is widely available. When information of air-oil or oil-water P-S relationship are required, in most cases, measurements should be conducted. On the other hand, established P—S measurement, such as that by Klute (1986), are primarily intended to measure air-water system. The characteristics of an air-water measurement protocol are: suitable only for a low wetting phase volatility, limited to gaseous type of non-wetting phase, and takes a considerable amount of time. In parallel with current growing interest in the multiphase flow problem, it is desirable to have an efficient P-S measurement technique that measures not only air-water systems but also air-oil and oil-water systems. Conventionally, P—S data are measured by allowing saturation to equilibrate under a series of fixed pressures imposed at the boundaries. While this technique is considered to be standard, a considerable amount of time is required to measure the entire P-S relationship. Typically, time required to obtain sufficient data to adequately represent the desaturation scanning path is approximately two weeks (Su and Brooks, 1980; Bear, 1982). Such a disadvantage significantly constrains the gathering a large amount of data, sample replication, and to measure data near the region of residual saturation. Yet, it is desirable to include these measurements. Replication and a large number of data insure a 6 representative sample mean as well as knowledge about the variance. Inclusion of data from the region near residual saturation reduces the risk of bias. The time required to measure P-S data using the conventional method is primarily dictated by the time required for the system to reach hydrostatic equilibrium. Evidently, any measurement readings in other than the hydrostatic equilibrium state would not represent the primary drainage or imbibition P—S relationships. (Topp et al., 1967; Smiles et al., 1971; Carey and Brooks, 1975; Parker et al., 1985). A P-S measurement method devised by A. T. Corey, as reported by White et al. (1970) is capable of reducing the measurement time while obtaining data that are representative of primary drainage. This technique reverses the conventional method by requiring capillary pressure to equilibrate under a series of fixed externally imposed saturation. Henceforth in this dissertation (based on the phenomenon to deduce the equilibrium state of the sample), the conventional method is referred to as the saturation equilibrium method. The method devised by Corey is referred to as the pressure equilibrium method. The pressure equilibrium method has been employed by White et al. (1970) and White et al. (1972) and was improved upon by Su and Brooks (1980). Typical total time required to measure 9 data points in air-water desaturation and imbibition scanning paths was less than 24 hours (Su and Brooks, 1980). While the pressure equilibrium technique has been shown capable of reducing the measurement time significantly, the obtained data have never been formally compared to that obtained using the conventional method. As such, the time reduction achieved may only be the result of a premature reading which might be biased relative to the 7 conventional measurement. If the converse was true and experimental time was indeed shortened, this technique offers an efficient method which might be firrther improved to measure any type of immiscible fluid. 1.2. GENERAL OBJEcrrVEs, APPROACHESAND PRESENTATION ORGANIZATION This study focuses on the effects of incorporating scaling theory in the predictive method to construct the three-phase kr-S-P relationship. Without using scaling theory, the P-S parameters required by the predictive method should be determined from P-S desaturation data of air-water, oil-water and air-oil systems. With scaling theory, however, only one of those systems is required to be measured. This study investigates the alternative of using either the air-water or air-oil systems as the base system. The air-water system is not directly involved in a three phase system. However, air-water data are commonly available. Conversely, the air-oil system involves in a three-phase system but P-S data of this system are not as widely available as that of air-water. The performances of these alternatives are investigated in two soil textures. The textures chosen are sand with uniform grain size distribution and loamy sand. Lenhard and Parker (1987) proposed the scaling theory as a way of improving the efficiency of the predictive method. Accordingly, the discussion throughout this dissertation is in the context of testing the hypothesis that the predictive methods with and without the scaling theory are not different regardless of the usage of the two-phase system to obtain the base parameters and the involvement of soil texture. 8 Additionally, this study verifies the pressure equilibrium measurement technique against the conventional measurement method as part of an effort of developing an efficient P-S measurement technique. The Su and Brooks (1980) apparatus is modified and improved to make it suitable for air-water, oil-water and air-oil P-S desaturation measurement. The concept of the modification is verified theoretically. An unsaturated mathematical model based on Richard’s equation is developed to simulate P-S measurement employing the pressure equilibrium method as well as the saturation equilibrium method. The modified apparatus is used to measure P-S data of air-water and air-oil systems. The pressure equilibrium is proposed to provide an efficient two-phase P-S measurement technique substituting the conventional saturation equilibrium method. As such, discussion in this dissertation is focused to observe whether the pressure equilibrium method is capable of reducing the time required to measure two-phase P-S data while conserving the characteristic of the primary desaturation. The two major issues addressed pertain to the implementation of the scaling theory and the verification of the pressure equilibrium measurement method. These two issues are discussed in Chapter II and Chapter III. Each of these chapters is presented as a self-contained section. Each one contains a concise discussion of the background, description of the material and method used and subsequently discusses the results and findings specific to the particular issue. Experimental protocols, measurement result, data reduction and statistical calculations are presented in the appendices. Background of the issue is presented in three separate sections. First, the Definition and General Overview section in this chapter clarifies the terminology used 9 and summarizes the major existing concepts pertaining to the pressure-saturation scaling theory. It is intended to draw a frame work from a broad point of view. Additional background, more specific to the separate issues is presented concisely in the introductory sections of Chapters 11 and 111. These sections are intended to put the addressed issues in perspective. 1.3. DEFHVITIONSAND GENERAL OVER VIEW. As mentioned above, this dissertation is concerned with the scaling method in the context of its implementation in the predictive method. The predictive method is a concept to provide an alternative means for solving multiphase flow equations. Prior to discussing the existing methods to scale the P-S relationship, the subsequent sections describes some terminology pertaining to multiphase flow used in this dissertation. Following that, the involvement of the predictive method in solving the multiphase flow is presented. The last section reviews the existing methods to scale the P-S relationship and highlights the porous medium hydraulic properties that are postulated by some authors as being the factors in the scaling method. The notations used in the presentation, generally follow the symbols used in Parker (1989). 1.3.1. Multiphase System When two or more fluids coexist in a soil pore and one or more of those fluids are in motion, the types of flow possible are either miscible or immiscible displacements (Bear, 1972). Miscible displacement refers to the case where the fluids are completely 10 soluble in each other. A distinct fluid-fluid interface which separates those fluids does not exist. In the immiscible displacement, each fluid for practical purpose exists separately. As such, a distinct interface between fluids exists and the interfacial tension between the two fluids is not zero. Fluid motion accompanying immiscible displacement is commonly referred to as multiphase flow. In this dissertation the term phase is used specifically to refer to the fluids involved and not the soil grains. A multiphase system refers to the coexistence of two or more phases in soil pore in conjunction with the multiphase flow concept. A multiphase system in general may involve two or more fluid types. However, a conceptualization of multiphase system which considers only three types of fluid is commonly regarded as sufficient and representative as the most general form of a multiphase system. This dissertation uses the term three-phase system in this context. In conjunction with an interest in subsurface contamination assessment, the three fluids involved in a three-phase system are ofien classified as water, NAPL and gaseous phases. The NAPL and gaseous phase are commonly referred to as oil and air phases respectively. Henceforth, the terms air, oil and water are used to refer to the phases involved in the three-phase system. A two-phase system is a multiphase system involving only two phases. The phases in a two-phase system are referred to as wetting-and nonwetting. The phase that tends to have direct contact with the solid is referred to as the wetting phase. The two- phase system consisting of air and water is a special case. This system describes condition in the unsaturated zone prior to contamination. The water and the air phases are 11 the wetting— and nonwetting-phases, respectively. Unsaturated flow is a particular case of two-phase flow. In this case, the interest is focused on the water phase motion only and the air phase is assumed stagnant at atmospheric pressure. 1.3.2. Multiphase Flow Equation and the Predictive Method. The equation governing the multiphase flow is derived fiom the generalization of Darcy’s equation combined with the continuity equation. (Bear, 1972; Corey, 1982; Parker, 1989). According to Parker (1989), the governing equation of the multiphase flow in a non-deformable porous medium, under the assumption of no sources or sinks of the phases involved, is given by: ¢a(pp ___p_S=) 3 k kn[3_:p a_z_] 3"! — x,'Pp tr 11 ""pp ax The subscript p is reserved to indicate the fluid phase (p=a, o, w) while i and j indicates directions (i, j=1 ,2,3) with repeating values indicating summation in tensor notation, x is the Cartesian coordinate, z is the vertical elevation, 4) denotes the porosity of the medium, I is time, pp is the phase density, km is the phase relative permeability, k,- is the intrinsic permeability tensor, Sp is the degree of saturation by the phase which is the fraction of the soil pore occupied by phase p, and hp is the pressure head of phase p. Note that h, S and k1. are interrelated. When capillary pressure h increases the saturation is decreased. Saturation decrement reduce the pore volume that participate in providing the channel for the fluid in question to flow. The relationship between kr-S-h is 12 formulated in the constitutive equations. The constitutive equation goveming multiphase flow consists of two functional relationships. The first relationship is a function relating the saturation of a particular phase and the pressure heads of the phases involved. The second is the relationship between the relative permeability and the saturation. The first function is referred to as pressure-saturation (P-S) function. In a three-phase system, this function is formulated as: Sp=f(ha’ ho, kW) [1.1] The P-S function in a two-phase system is a special case that involves only the pressures of the wetting and nonwetting phases. Interaction between these pressures can be represented only by the capillary pressure which is defined as: PM 5 Pp — Pq Here PM is the capillary pressure, Pp is the pressure of the nonwetting phase p and Pq is the pressure of the wetting phase q. The capillary pressure Pm is commonly expressed in terms of pressure head relative to some reference phase as: m=mem4 where hm is the capillary pressure of the pertaining immiscible phases pq in term of pressure head, g is the gravitational force, and p“ is the density of the reference phase. The two-phase P-S fiinction will then take a simple form as: s= 11h”). [1-2] The procedure to solve the flow equation commonly involves a step to reduce the nonlinearity of that function by applying the chain rule to the term of the derivative of 13 saturation with respect to time. The form of the chain rule for the case of three-phase flow is written as: as, as, ah, as, ah, as, ahw = 1-3 at ah, at + ah, at + ahw 6t [ 1 This means that the constitutive equation that describe the relation of oil saturation and the capillary pressure heads should be formulated in the form of Eqn 1-1. Construction of such a relationship requires direct measurement of the saturation and the pressure head of the air, oil and water phases simultaneously. This requires three phase measurements which are not commonly made. Furthermore an established protocol and technique to conduct three-phase P-S measurement is not yet available. Lenhard and Parker (1988) are the only investigators to report three-phase measurements for P-S relationships. The common method for solving the three-phase flow equation is by employing the predictive method proposed by Leverett (1941). The predictive method assumes that the soil matrix is a hydrophilic porous medium. Under this assumption, the water phase has a direct contact with the soil grain while the oil phase always resides between the water and the air phases. Provided that the oil volume is high enough to cover the water surface, the air phase would have no direct contact with the water phase. The predictive method further assumes that the total liquid saturation (ie. oil plus water) is dictated only by the air-oil capillary pressure and that the water saturation is dictated only by the oil-water capillary pressure. Employing the above point of view, a three-phase system consisting of air, oil and water can be conceptualized as two interdependent two-phase systems. The first system is 14 a two-phase system consisting of a total liquid phase and an air phase. Since the saturation of the total liquid is a function of the air-oil capillary pressure only, the total liquid-air system is assumed to be identical to an air-oil system. Secondly, within the total liquid phase, the oil-water two.phase system exists. Employing this approach, the above chain rule would only involve the derivative of total saturation with respect to air-oil capillary pressure and derivatives of water saturation with respect to oil-water capillary pressure. Therefore, it only requires knowledge of the air-oil and oil-water P-S relationship. These relationships are easier to measure compared to a true three-phase P-S relationship. 1.3.3. Methods to Scale the Pressum-Saturation Relationship. ”“2 [E l' Elf. -S 'Sl' P-S scaling theory was developed to relate a P-S curve across porous media hydraulic system. This knowledge allows transformation of the P-S relationship from certain system to another system. Prerequisite to this knowledge is an understanding of any determinants that dictate the P-S relationship. According to Bear (1972), because of the dependency of capillary pressure P on interfacial tension 0 and radius of curvature R, on a microscopic level, macroscopically it depends on (a) the geometry of the void space, (b) the nature of the soil-fluid interaction and (c) the degree of saturation S. On a macroscopic level, the porous medium property associated with the characteristic of the void space geometry is the pore size distribution. The factor of soil-fluid interaction includes contact angle (.0 as suggested by Morrow (1976) and Demond and Roberts 15 (1990) among others and wettability as suggested by Scheideger (1969) and Adamson (1967). The study by Anderson (1987), however, concluded that the PS relationship was insensitive to the wettability at large saturation ranges as it is overwhelmed by the effect of surface roughness. This dissertation assumes that the effect of wettability is insignificant. The third factor which is the dependency of capillary pressure to the saturation, is formulated in the form of P—S relationship curve. Therefore, it is obvious that inherent in the P-S relationship there exist interactions between P, R,, 0', co and S. The following paragraphs discuss further the conceptualization of these interaction. Capillary phenomena in a soil matrix are commonly explained based on Laplace’s equation. Under this concept, the capillary pressure Pp, is related to the radius of curvature Rc of the interface by: 2 a Note that, a P-S relationship describes the relationship between the capillary pressure and saturation at a macroscopic level. Hence, a given capillary pressure does not necessarily refer to a particular radius of curvature in a particular soil pore. On a macroscopic level, R, is rather the radius of a conceptual capillary interface in the soil- fluid system that exits as the response of the system when subjected to that capillary pressure. To distinguish between the microscopic and the macroscopic radius of curvature let the macroscopically measure be denoted as R. Carey (1982) defined R as the hydraulic radius equivalent to that commonly used in hydraulic engineering. Bear (1982) put this R in the context of “statistical average taken over the void space in the vicinity of a 16 considered point in the porous medium”. The most important characteristic of R is that it decreases as the capillary pressure Pp, increase and when R decreases S decreases proportionally. The conceptual relation between S and R is commonly associated with pore size distribution and is taken to be a porous medium dependent parameter. In this dissertation the term of pore size distribution function is used specifically to refer a function that relates S and R in the form of R(S). Based on the conception above, a P-S relationship in the form of saturation as a function of capillary pressure 5(qu) can be factored into two relationships. The first relationship relates the capillary pressure and the radius of curvature Pp,(R). This relationship is formulated following the concept behind Equation 1-4. The second relates the macroscopic radius of curvature and a function of saturation which is conceptualized as the pore size distribution function R(S). Therefore, the function of S(Pp,) can be viewed as the result of substituting the R(S) into Pp,(R). Following this conceptualization, an S‘( P3,) curve of soil-fluids system A can be transformed to become 53(P;q) of soil-fluids system B by scaling their capillary pressures such that at any given S there exist a relation of P3, (R) = [3(5) Pg, (R). The factor B is referred to as scaling factor and hence is defined as: 13(8); Pym/Pia). [1-51 The general formulation of the scaling factor can be obtained by substituting Equation 1- 4 into Equation 1-5 and is written as: 17 a” R‘(S) s = — 1-6 [OM . 1 Within a soil-fluid system the interfacial tension is commonly regarded as a constant. The macroscopic radius of curvature R, however, varies as the saturation changes. When the saturation in the system A and B is changed with the same amount, the R in systems A and B is not necessarily changed with the same amount. Therefore, the ratio RA/RB varies when the saturation is varied. Equation 1-6 accommodates this characteristic by expressing the macroscopic radius of curvature R as the pore size distribution function R(S). Therefore, to obtain a scaling factor that transfonn P-S curves across two different two-phase systems and porous media requires knowledge of the interfacial tensions and pore size distribution functions. The scaling factor as formulated in Equation 1-6 might obscure its practical value. To date, most common technique to deduce the pore size distribution function requires P-S data When P—S data are measured and available for both systems A and B, scaling procedures are not required. Some investigators have proposed relating the pore size distribution function to the grain size distribution. For example, Smith (1933) and White et al. (19703) simplify the shape of the soil grain into sphere shapes to formulate the void structure. Even after such simplification, the results were still too complex for practical use and yet too specific to serve any general purpose. Apparently, procedures for scaling the P—S curve are only feasible within identical soils or at least within similar porous media as defined in the similar media concept by Miller and Miller (1956). 18 WW Leverette (1941) proposed a function, referred to as .1 function, to represent a general form of the P-S relations. This function is given as: Mara Here k is the saturated permeability and (p is the porosity. In his report, Laverette (1941) showed that Equation 1-7 coalesced the P-S curves measured in several unconsolidated sands. Based on this notion, the J function is commonly regarded as function that can be used to scale P-S curves not only across the two-phase system but also across porous media. To scale P-S relationship, Equation 1-6 utilizes a factor which transforms P-S relationship of a particular soil-fluids system to another. The J function follows different approach. This method transforms any P-S curve into its stande form, which is the J function as written in Equation 1-7. In order to compare the scaling concept behind Equation 1—6 and Equation 1-7, the scaling factor is reformulated based on Equation 1-7. The factor which scales the capillary pressure in system A into system B at a given J(S"') is given by: "7 B = 9;. 8A [1-8] a k3 A. ' According to Bear (1972) and Corey (1982), based on the concept of relative permeability proposed by Burdine (1956), the term of (k,/ OB R _____._fi 1-9 B a, R r 1 B ax Be: When comparing Equation 1-9 with Equation 1-6, it shows that the J function scales the P-S relationship across porous media using a constant (the ratio of Rm“) rather than using the pore size distribution function. It is expected that this particular scaling method will not always hold. In a later study using air-water system in various consolidated and unconsolidated porous media (Rose and Bruce, 1949), shows that plots of the J function are unique to the porous medium type. Based on this, further studies (Rapoport 1955; Richardson, 1961) that advocated the use of this function to scale the P-S curves, restricted the use of their method for identical porous medium. Under this restriction the ratio of Rm, becomes unity. In turn P-S scaling by J function accommodates the scaling across two-phase systems only and excludes the feature of across porous media. liii S' 'l M l' E 1 Mil lll'll 1222! Miller and Miller (1956) proposed the similar media concept as a method to scale the unit of a physical model to simulate unsaturated flow in porous media. Here two porous media are considered as similar when the grain size axis of their grain size frequency distribution curves can be scaled to each other by a constant factor d. 20 According to this method, all P-S curves within similar media will coalesce into a standard curve when their capillary pressures are scaled by: P=[&9—]d {1.10} where P is the scaled capillary pressure. Using superscript A and B to denote the oil- fluids systems, the scaling factor equation derived based on Equation 1-10 will take form of: oB dA “-077 [1-11] This formulation is equivalent to that formulated in Equation 1—9 where the scaling across porous media is accommodated by including the term of the ratio of d. It should be noted, however, to satisfy the similar media condition the two porous media should have similar grain size frequency distribution. Inherently, this condition requires that the two porous media have similar shape pore size distribution curve. When this condition is applied in Equation 1-6, the term of the ratio of the pore size distribution functions will be constant. As the result the scaling factor B in Equation 1-6 will take form of constant rather than function. Thus, the scaling factor as formulated in Equation 1-6 is in agreement with that in Equation 1-11. Study by Klute and Wilkinson (1958) confirmed that when the criterion of similar media was satisfied P-S curves of several different-but-similar sands were coalesced. The study by Elrick et al. (1959) using soil containing not only sand fractions but also silt and clay, however, concluded that when the theory was checked against condition far from 21 ideal a deviation was observed. Moreover, this study demonstrated the sensitivity of the scaling procedure to any pore size distribution variation. .v Hi, 41.! I a A special case of interest in scaling the PS curve is when this curve is expressed following the equation proposed by Brooks and Carey (1964). Here the P-S relationship is formulated as: P A s, = [—4-] [1-12] PP? where S, is the effective saturation which is defined as (S-S,.)/(l-Sr ), Sr is the residual saturation, Pd is the displacement pressure and A is the pore size index. All the required parameters to scale the P-S curve are already inherent in this formula since all parameters involved have a clear physical interpretation. When S is expressed in term of S, the differences between the two P-S curves due to their differences in Sr are factored out. Thus P-S curves from any two-phase system will have common range of S, from 1 at full saturation to zero at residual saturation. The displacement pressure Pd is the capillary pressure associated with the largest pore in the soil-fluid system in question. When the capillary pressure in a P-S plot is factored by Pd, a standard capillary pressure equivalent toPused in Miller and Miller (1956) is obtained. The parameter It. represents the pore size distribution characteristic 22 and was postulated to be a porous medium dependent parameter. In turn, any two P-S curves that have identical 2. will satisfy the requirement for similar media. Therefore, if those two curves are expressed in term of effective saturation and reduced capillary pressure they will coalesce into a single curve. Based on this particular characteristic, the parameters Pd and A are often used as parameters upon which the characteristic of the scaling method is observed such as in studies by Demand and Roberts (1990), Lenhard and Parker (1987) among others. Even though P-S formulation given by Eqn 1-12 provides all properties required for scaling procedure, the use of this function to conduct scaling across porous media is not practical. When a scaling equation comparable to Equation 1-6 is derived based on Equation 1-12, the resulting equation will be a function of A. This means that the knowledge of pore size indexes of the porous media involved is required. As mention previously, to date, these indexes can only be obtained by means of P-S measurement. manmmammsmzmm As part of their effort to develop an efficient method to construct a three-phase constitutive relationship, Lenhard and Parker (1987) proposed their scaling theory. This method was intended to scale the PS curves across two-phase systems within a single porous medium. Therefore, their scaling factor took the form of Equation 1-6 in which the term of RA(S)/RB(S) is set to unity since the soil in systems A and B are identical. Further detail description of their theory is presented as part of Chapter III. 23 Waugh. In its original form, the contact angle a) is involved in the Laplace equation. When (D is included, Equation 1-1 takes the form: a” case)” )[R‘(S)] [1-13] [3(5) = [0* -cosa)A R"(S) Among investigators, the inclusion of contact angle a.) is still a point of controversy. Several investigators who worked with artificial unconsolidated porous medium (Bethel and Calhoun, 1953; Morrow, 1976; Demand, 1988; and Demand and Roberts, 1990) observed the effect of a) on the P-S relationship. Subsequent work (Demand and Roberts, 1991) on a sand-organic liquid-water system showed the discrepancy between the observed and scaled displacement pressure when (u was excluded. Exclusion of a) (Leverette, 1941; Rapoport, 1955; Parker, 1987) was based on the argument that conceptualization of soil pore geometry as a tube was not appropriate. Hence, the radius of tube should be expressed in term of radius of curvature. Arithmetically, when the radius of tube is expressed in term of radius of curvature the contact angle a.) will be canceled out. By not including contact angle to, Dumore and Schols (1974) were able to successfully scale their data of heptane-water and toluane- water in consolidated sample. Anderson (1987) observed that both in strongly water-wet and oil-wet consolidated rock sample, the capillary pressure is insensitive to a) at large range values of (0, especially in the drainage path. 24 It is interesting that most of the investigations that reported an observable effect of a), used an artificial unconsolidated porous medium or clean sand. On the connary, those who used natural rock core or unconsolidated soil concluded that due to the existence of surface roughness, the apparent effect of contact angle is almost non-existent. In this dissertation the effect of contact angle co is assumed negligible. [SIZE 1 CSI' Mil $11.1 Several studies indicate the dependency of the parameters Pd and A on soil texture. This dependency is based on the assumption that in an air-water system, the pore size distribution characteristic is strongly affected by the grain size frequency distribution. Several studies even advocated deduction of the P-S curve based on the parameters associated with the grain size distribution (McCueen et al., 1981; Campbell, 1985; Mishra et al., 1989 among others).These approaches were based on regressing an empirical function onto grain size distribution data. Brakensiek et al. (1981) used the Brooks-Corey (1964) equation as their regression function. In this study the parameters Pd and A were classified into the twelve Soil Conservation Service textural classes. A study by Cosby et al., (1984) confirmed that, statistically, the magnitudes of Pd and A were strongly dictated by the soil texture. Similarly, the study by McCueen et al., (1981) found that these parameters are unique for each soil texture class. CHAPTER II THE A CCEPTABILITY OF THE PRESSURE EQUILIBRIUM METHOD To MEASURE DRAINA GE PA TH 01“ San PRESSURE-SA TURA TION 2.1. BACKGROUND The pressure-saturation (P-S) relationship is one of the constitutive relationships governing multiphase flow through unsaturated porous media. It expresses the functional relationship of the equilibrium wetting phase saturation at a given capillary pressure. Petroleum engineers employ the P-S relationship (Collins, 1961), while soil scientists present the same information in a water retention or moisture characteristic curve (Child, 1940). The soil scientist’s term, however, is commonly restricted to air-water systems, while the more general P-S relationship is employed with a variety of immiscible fluid systems. In parallel with current growing interest in a multiphase flow problem, it is desirable to have an efficient P-S measurement method that measures not only air-water systems but also air-oil and oil-water systems. Conventionally, the P-S data are measured by allowing saturation to equilibrate under a series of fixed pressure gradients imposed at a boundary. While this technique is considered to be a standard, it is a very time,consuming method. The typical time period required to obtain sufficient data to adequately represent the desaturation scanning path is approximately two weeks (Su and Brooks, 1980). Such a disadvantage places significant limitations in collecting larger numbers of data, replicating the sample, and measuring the data near the residual saturation region. Replication and the large numbers of data are 25 26 desirable to obtain representative mean and to increase precision reliability. An inclusion of data from the region near residual saturation eliminates the risk of bias. The length of experimental time in the conventional P-S measurement is primarily dictated by the length of time required for the system to reach equilibrium. Up to this time, studies that investigated the feasibility of conducting measurements during the transient state have had qualified success. At a given saturation, capillary pressure measured under dynamic flow tended to be higher than those in a static condition (Topp et al., 1967; Smiles et al., 1971; Carey and Brooks, 1975; Parker et al., 1985). Apparently, in order to obtain P-S data that share similar characteristic to the conventional method, any alternative method that measures the P-S relationship under an equilibrium state is more sound. An alternative method that measures the P-S relationship under the equilibrium state was devised by A. T. Carey (White et al., 1970). This method had been used in White et al. (1970) and White et al. (1972) and was improved by Su and Brooks (1980). In contrast to the conventional method, suction was applied in two stages. In the first stage suction was set to drive a reasonably quick desaturation. While in the second stage, after a volume of desaturation had been obtained, the suction was adjusted to balance the equilibrium wetting phase pressure. Implementing this method to measure air-water data using 50 different soils covering a wide range of properties, Su and Brooks (1980), observed that typical experimental time to complete both desaturation and imbibition scanning paths was less than 24 hours. The transient time for the pressure to reach ll 27 equilibrium after the desaturation step took only a few minutes at high saturation and a few hours at low saturation region. While this technique significantly shortens the experimental time, the obtained data have never been formally compared to data using the conventional method. As such, the time reduction achieved may only be the result of a premature reading which could introduce bias relative to the conventional measurement. If the converse were true and experimental time was indeed shortened, this technique offers an efficient method which might be further improved to measure any type of immiscible fluids. Note that in the Su and Brooks (1980) apparatus, a wetting phase evaporation was allowed to accelerate the desaturation rate. If the wetting phase is a volatile liquid, the desaturation might be dominated by evaporation. This condition may introduce a problem in adjusting the second stage suction. Furthermore, this approach prohibits the use of a liquid nonwetting phase. However, applying a high positive pressure on the nonwetting phase to desaturate the sample may resolve those problems. This modification, however, is only justified providing that it will not affect the characteristics of the method for reaching equilibrium quickly. According to Su and Brooks (1980), the equilibrium state was quickly reached due to the adjustment made in applying the second stage suction. Unfortunately, no further explanation was given. In their measurement technique this adjustment was conducted in the following manner. After a finite volume had been drained from the sample, the outflow was directed toward a vertical capillary glass tube containing the meniscus of the air-wetting phase interface. At this instance, the air pressure on that 28 interface was adjusted to maintain the meniscus level at a prescribed datum until this meniscus became stationary. Leveling the meniscus during the transient state will permit determination of the wetting phase pressure. In our opinion leveling the meniscus changes the condition at the external face of the plate from a constant head to a zero flux condition. This suggests that the important behavior underlying the Su and Brooks (1980) method is that the transient state within the soil sample which occurs under a no- flow boundary condition will be terminated more quickly than it would under a constant head boundary condition as in the conventional method. Based on this premise, applying a positive pressure to desaturate the sample should give no effect to the basic principle of the method. Furthermore, replacing the capillary tube with a zero displacement pressure transducer seems to be more effective in providing a zero-flux boundary condition while monitoring the wetting phase pressure. F urthermore, the use of a transducer eliminates the laborious second stage suction adjustment. If the transducer used is capable of measuring the pressure difference between the wetting and nonwetting pressure, a more efficient measurement protocol can be obtained. In order to distinguish between the Carey method and the conventional method based on the phenomenon used to indicate the equilibrium state of the system, the Carey method will be referred to as the pressure equilibrium (PE) method and the conventional one as the saturation equilibrium (SE) method. The current paper studies the implementa- tion of the PE method by developing a P-S measurement apparatus based on the concept referred to in the previous paragraph. The apparatus developed was tested to measure 29 air-water, air-oil and oil-water desatmation P-S relationships in a natural loamy sand soil. To justify the acceptability of the PE method, the air-water and air-oil data obtained were compared to those obtained by the SE method. In addition, a 1-D unsaturated flow mathematical model was used to verify the method theoretically. 2.2. MEASUREMENT METHoD AND MATERIAL The PE method measures the P-S desaturation data by imposing a finite desaturation on a presaturated sample and subsequently allowing the fluids in the sample to equilibrate under a no-flow boundary conditions prior to recording the equilibrium data. The equilibrium state in this case is indicated by a stationary wetting-nonwetting pressure difference. With the SE method, desaturation is allowed to occur under a fixed wetting-nonwetting pressure difference imposed at the boundaries until equilibrium is reached as indicated by an outflow dirninishment. In order to compare the performance of these two methods, a P-S measurement device capable of performing either PE or SE methods in measuring air-water, oil-water or air-oil desaturation data was developed. This device was developed utilizing a standard pressure cell by Soil Moisture (Model No: 1400) as shown in detail in Figure 2-1. The connection between the pressure cell and the other parts of the device is shown in Figure 2-2. Three pressurized fluid sources are available: air, oil and water. The air phase source is obtained from the building line through a cascade of coarse and fine regulators. The oil and water phase pressure chambers provide the source of the oil and water phases respectively. The pressures of these liquids are governed by the air pressure in the chamber. This pressure 30 : 1. Top assembly 4; 2. Inflow port 3. Filter paper 4. Wing nut 5. Wm O-ring ———> 6. Soil sample -——> 7. Brass retaining cylinder ———> 8.Filterpaper .——’ 5. Viton O-ring 9. Viton O-ring 12. l-baggiggrlme % 10. Bottom assembly plate ; 11. Outflow port Figure 2-1: Pressure cell detail 31 an “on magma—u 2E. ”N1N 85mg an; “Banana «BB—.8 880:00 Avg: Danna coca? ”firmed .833 m5 EMT—MW Av Iv kw N» 9 .5 g .505me ................. OBmmOhnH Nm m m Hausa—95E c> D5802 AN 2, a, T 8: can :0 AT 9 D kw _ AV 8: SSE 5.. W m v a a N; E: §58> 09:5... financed “EV/luv ~85 h< SEER If 32 can be altered utilizing the building pressurized air. The zero pressure of these liquids is defined as the pressure which level the air-liquids interface in burettes B1 and BZ (Figure 2-2) at the center of the cell. The applied building air pressure at any instance is monitored via a mercury manometer. The connections between the fluid sources and the cell are designed such that the air and water phases can be used as the nonwetting and the wetting phases respectively, while the oil phase may function as either one. The nonwetting phase line is connected to the top of the cell. The positive pressure applied on the nonwetting phase drains the wetting phase from the bottom of the cell. The outflow is directed into a vial collector in which the pressure is atmospheric. The outflow line to the vial is controlled by a valve, referred to as the collector valve. The collector valve is opened only when desaturation is permitted. A differential pressure transducer is installed to measure the pressure gradient between the wetting and the nonwetting phase pressures. Typical procedures in conducting measurement are: packing the sample, saturating the sample, calibrating the transducer, conducting P-S measurement based on PE or SE method, and determining the residual saturation by an independent means. Note that in the SE method the nonwetting phase pressure is determined using a mercury manometer. Therefore, transducer calibration is not required. The loose soil sample with approximately 2-3% water content was packed in a brass retaining cylinder (2 ‘/4"OD, 0.065" wall thickness, and 3 cm height) using a compactor with a 1.6 kg weight falling 10 times from a 7 cm height. This procedure proved to be effective in maintaining consistent bulk density between samples. 33 Subsequently, the packed soil in the retaining cylinder was assembled with the other pressure cell parts. The pressure plate was oven dried prior to being installed in the cell. The final cell assembly (Figure 2-1) was placed and connected to the other part of the measurement device and ready to be saturated. To saturate the sample, the top of the cell was connected to a vacuum source. The bottom of the cell was connected to the wetting phase line and the wetting phase pressure was adjusted to be atmospheric at the center of the cell. Subsequently, in order to evacuate the soil’s pore gas, a 900 mbar vacuum was applied to the nonwetting phase without allowing wetting phase to flow into the cell (the valves V10 and V11 were closed, see Figure 2-2). After several hours, wetting phase absorption was allowed, by opening either V0 or V11, until the top surface of the soil sample became slightly oversaturated. At this instance, the suction was dropped to atmospheric pressure. Using this technique, air bubbles trapped during the advancement of the wetting phase were confined and collapsed. Next, the cell was removed and weighed in order to obtain what is refened as the initial weight of the cell. Subsequently, the cell was reconnected to the wetting and non-wetting phase lines. These lines, including the voids in the transducer, were filled with the pertinent fluid and the pressure of each fluid was set to zero. The cell was then left to equilibrate at this setting for approximately 14 hours after which the transducer was calibrated. The transducer is calibrated against a mercury manometer in the following manner: The wetting phase pressure, either water or oil phase pressure, is fixed at zero. Once this pressure is set, the cell is removed and the wetting and nonwetting phase lines 34 to the cell are sealed. Subsequently, the non-wetting phase pressure is varied stepwise from zero to maximum and back to zero. The transducer reading then, is calibrated against the applied nonwetting pressure as indicated by the manometer. This calibration is conducted prior to any experimental run. The PE method is employed to determine the P-S relationship in the following manner. After the transducer is calibrated, the nonwetting phase pressure is increased to approximately twice the displacement pressure without permitting any desaturation; this aspect of the procedure will be discussed in the Result and Discussion section. Since no desaturation has occurred, the increment in the nonwetting phase pressure should not , affect the transducer reading even though a small disturbance may be registered. Once the fluctuation due to this disturbance has diminished, the first capillary pressure measurement is taken. The initial capillary pressure and saturation are therefore known. The subsequent desaturation is initiated by opening the collector valve. After a finite volume has been collected the valve is closed. The collected wetting phase is weighed to determine the current state of saturation in the sample. The capillary pressure reading is taken after the transducer reading stabilizes. The procedure is repeated again starting at the step of opening the collector valve. When the drainage rate becomes too low, the non- wetting phase pressure is increased. The desaturation process is continued until saturation in the sample is near residual saturation. With the SE method, the saturated cell is connected to the wetting phase burette and the air phase lines are left open to the atmosphere. The pressure in the wetting phase fluid is adjusted to place the interface in the burette at the center of the cell and the fluids 35 are allowed to equilibrate overnight. The wetting phase pressure at the center of the cell being atmospheric results in a calculated pressure difference, relative to the air pressure, of zero. This what is recorded as the capillary pressure at initial saturation. The next desaturation point is obtained by applying an incremental step of positive pressure to the non-wetting phase while the outflow is directed to the collector vial and the line to the burette is closed. Drainage is permitted for at least 22 hours. Either the mercury manometer or the transducer can be used to determined the nonwetting phase pressure. However, it is preferable to use the manometer, since it does not require calibration. The wetting phase pressure is the distance of the wetting phase ‘ fluid surface in the vial to the center of the cell. The saturation in the sample is determined from the weight of the drained wetting phase. The procedure is repeated for the subsequent desaturation points. After the last measurement point is obtained, the cell is weighed. The difference between this weight and the initial weight of the cell reflects the total amount of the wetting phase drained. The agreement between this amount and the sum of the wetting phase collected in the vial is used to determined the mass balance error. The cell is then taken apart and all the soil in the retaining cylinder is weighed prior to being oven dried. The difference in weight before and after oven drying is the residual saturation and the oven dried weight divided by the brass retainer volume is the bulk density. This apparatus was used to examine the performance of the pressure equilibrium method in measuring air-water and air-oil pressure-saturation data. Likewise, the same apparatus was used to conduct the saturation equilibrium method measurements of the 36 air-water and air-oil systems. The porous medium used was Metea soil with a texture of 6% non-expansive clay, 11% silt and 83% sand. It contained 1.1% organic matter after dry sieving to pass 850 micron to exclude larger gravel and other organic matter. In order to have a stable packed sample, the water-phase used was a deaerated 0.005 M CaSO4 solution as suggested by Klute (1986). The oil-phase was tetrachloroethylene (PCE) at 99.99% purity with a density of 1.623 g/cm3. Building air was used as the air phase. 2. 3. RESULTAND DISCUSSION 2. 3. 1. Pressure-Saturation Measurement Simulation. The PE method is based on the premise that transient saturation profiles within the soil sample will reach equilibrium relatively quickly when occuning under a zero flux boundary condition. To verify this notion theoretically, a numerical model was developed and used to compare the time required by the PE and SE methods to collect P-S data. This model was a l-D finite difference model, in which the air mass density and viscosity are assumed to be negligible. The air pressure and water mass density are assumed constant with respect to time and space. The appropriate governing equation for this case is given by Richard equation: US 3 8h earn-air, -k,(-5;-l)] {2.1} where Q is the porosity (cm3/cm3), ks is the saturated hydraulic conductivity (cm/sec), t is the time (sec), h is the pressure head of water (cm) and z is the vertical distance (cm). The 37 bottom of the soil sample was set at z=0 cm and the top was at F 3 cm. This height was divided into 10 uniform elements. The water saturation (S) and relative permeability (k,) are assumed to follow van Genuchten (1980) closed-form relations: 1-1/ s = s, + (1 + SJ + (ah)"] " [2-2] 2 -1 / k. = S,"2[1-é— 53/04))" ) "] [2-3] where Se is the effective saturation defined as (S- S,)/(l - r). Here Sr, (1 (cm") and n are the van Genuchten P-S parameters. Note that the model was simplified by omitting hysteresis in the P-S relationship and Equation 2-2 was assumed to hold regardless of the saturation history. A simulation of the transient conditions within the cell when the PE method is employed begins at time t, when drainage is initiated by opening the collector valve. The initial pressure head distribution h(z, ttd, the S(td) profile changed to becoming the redistribution transient profile and formed the final profile at the end of simulation time. The relative position of these transient profiles relative to the final profile suggest that, in the region of z bellow approximately 1 cm, saturation was increasing. As a result, the data obtained in the PE case might not represent the true primary drainage and a bias relative to the SE data would then be observed. The hysteresis issue may not be the only determinant which could bias the data obtained by the PE method. The dynamic state of the fluids at the time when the measurement readings were taken might also be an important factor. T opp (1967), Smiles et al. (1971) and Corey and Brooks (1975) observed that P-S measurements taken under unsteady, steady and hydrostatic states yielded different data sets. They found that at any given saturation, capillary pressure measured in a steady or unsteady state was higher than that measured in hydrostatic equilibrium. The model developed cannot be used to address any discrepancy associated with these hysteretic and dynamic state issues. The 44 hysteresis was omitted in the model so that similarity in the final saturation profile S(tc) obtained by the PE and SE simulations was artificially imposed. With the model, the occurrence of the hydrostatic equilibrium state can be precisely determined by monitoring the change in the pressure distribution profile. During an actual measurement, however, such a profile might not be practical to be acquired since the height of the pressure cell is only 3 cm. Thus, any means used to infer the state of hydrostatic equilibrium during the actual measurement could be less accurate than that employed in the model. In turn, the P-S point obtained in the actual measurement employing the PE method might be unintentionally taken while the fluids were still at an unsteady state condition. As the . model cannot be used to address theses issues, it is necessary to conduct an inspection on data obtained from experimental measurements. 2. 3. 3. Pressure-Saturation Measurement employing the Pressure Equilibrium and the Saturation Equilibrium methods. To verify the PE method experimentally, the method was implemented to measure P-S data for air-water, and air-oil systems using Metea soil. The PE method would be justified if it was unbiased relative to the SE method. In order to achieve this condition, the following criteria should be satisfied. First, the measurement should be taken as close as possible to the hydrostatic equilibrium state. Second, any discrepancies observed should still be within tolerable limits such that the P-S information obtained is statistically comparable. This information would include the P—S curve derived from the measured data and the P-S parameters describing that curve. Subsequent discussions will be presented from this point of view. 45 232lE'l'i'Zi"!!! 511.5“! In order to verify the sufficiency of the time given for the system to reach hydrostatic equilibrium for each measurement with the PE method, the transient behavior of the capillary pressure as indicated by the transducer during redistribution was recorded approximately every 5 minutes. Typical plots of such data are shown in Figure 2-5. The open squares in Figures 2-5a through 2-5d denote the capillary pressure readings. The series of points belonging to a particular desaturation step end with a solid point for which hydrostatic equilibrium was assumed to be achieved. The criterion for identifying these end points was that the rate of capillary pressure change be less than 0.25 mbar per minute. If this rate was achieved at the second point, additional readings were taken. The label next to each solid point shows the degree of saturation at that step. At S = 0.6T and 0.25, the pressures were equilibrated overnight (785 minutes). The horizontal lines indicate the pressure level at which the non-wetting phase pressure was set (labeled as hp in the Figure 2-5). Having knowledge of the trend of the pressure head profile in the vicinity of the pressure plate as shown in Figure 2-3, it was possible to reconstruct the redistribution curve of each desaturation step. Each redistribution curve in Figure 2-5 (solid line) was traced manually through the open points that belong to a particular step ending at the solid point. The beginning and the end points of this solid curve were then manually extrapolated (shown as dotted lines) so that at the beginning of each desaturation step it crossed hp and at the end point it continued the trend of the solid line. It is assumed that the quasi-equilibrium state was already reached at the points where the redistribution 46 9...... .23: emioanmb =5 .3 :2qu Mi Rafi 2: 53> woman—toga EoEoSmmoE 5:3th 35 82m 3% Mi 35890 05 A8 98 peace um 05 flame—95 EoEuBmmoE a E can covenmummvfi meta—e 32.355 05 .3 @0388 we eat 8 “8&2 5T5 woman". e8: 2529:“ .533 €13 ”9N 053m GQSEEV 05:. m - gage 8e Inuit I one. 8.: ; ca... I one to»... I 8.0.. 39. In...» -31.... ------IIIIEMoum *8? F mud m.o mad 0 A3 _I I seem w. 0 mm o uuuuuuuuuuu mud _ # o8.W m 1., fl svwm . , - i llllllllllllllllllll .l. L I com. on .5 8e. us 80— 00°F Over ONE. 08? SO 8m 03 O E o Scum . 8- oo. 2 o 18.0”..- ---- orq .. c ........ II IIIHJ .. e .I. 8 O IIIII _I llllll II n 8. $ m a a u a 0,. Eomhfl.‘ eramdm 09 d IIIIIIII a ,I IIIII ._ out M 80 mN—dafl ov—u m com m 82 o8 8a SN 8m 8m 8.. 8m 8m 8. o 'lt o — Owl "- U: 9 omN \an ADV om- WW «gum -) 852: mm m s m m 8% 533;? A8 7. AF “Hf com I I I a. IIIIIIIIIIIIIIII I I I I - m as.-- .. 8m 8m 8. 8. 8 o o . I 8a 3.3 IIIIIIIII 9- fl. . 3 I - - . .3... 3W“ cow «3 as «3 o; E. bnpao. tad/t 8am” Iloddddz calm/If: I j uh} I II It? I I Ir j - OVI( Eulmanlnfi‘l 'lrIILIIILILlL-llrllpl emu 47 curves became horizontal. Beyond these points, the rate of the capillary pressure change was assumed to be insignificant. It was observed that the solid points were laid in the horizontal region of the curve. This indicates that the criteria of 0.25 mbar per minute was sufficient to assure the quasi-equilibrium state of the system. Inspection of the constructed curve indicates that, with the air pressure setting hp as shown in Figure 2-5, in an air-water system with high saturation (980%), 20 minutes was sufficient to reach quasi-equilibrium. At moderate saturation (30% h b h p n p h P r- c— .. 8 .t 8 .t 3 1 t a t 8 , a 2. e - 8 .. .. a m. on .. E38 m n a h g men use: E e: E a SN 38% S 2 a... I. 2 2. a .“ t 238 v "on: team 2230 A3 “ I no or r r I I! r ass 8v use: new use: 3 e3... 84 an .2 .2 .8 .2 .2 u .5: 88288338 (mm)°‘q Emggav 8— sane. app ope a... I v..o . N..o oc _ _ _ _ e , , , , / / / I I / In Ill/l / Ill/l I l I / r2 / / u. x m / 12 \WI , < w .... Ta g88vdufi1§d1m6| & p 52553 ...a 77 (b) Metea Soil Figure 3-7: Three-phase oil relative permeability obtained by method 1 (solid lines) and method 2 (dashed lines). 13.12.34“ 80W km 0.01 kgw : (a) Ottawa Sand ha0 (mbar) 10 how= 100 mbar Sw= 0.00 how (mbar) 1 0 (c) Ottawa Sand h = 4 mbar s?- 0.95 0.01 78 kg“ ~ 'n - ~ " u (b) Metea Soil = 100 mbar 0W sw= 0.05 0.01 ‘ how (mbar) _ .100 .-- U... ,- a ,- 4 (d) Metea Soil h = 4 mbar 5? 0.99 Figure 3-8: Three-phase oil relative permeability obtained by method 1 (solid lines) with their uncertainty limits (shaded area) and those obtained by method 2 (heavy dashed lines) with their uncertainty limits (light dashed lines). I342” 1.3:“ 0.01 (a) Ottawa Sand h = 100 mbar 0W sw= 0.00 10 0.01 (c) Ottawa Sand f hm: 4 mbar s,= 0.95 ....... ..... ,- _o ,- .0 . o .0 r 79 0.1 1 k2,?“ (b) Metea Soil hm (mbar) 10 ...... ...... ..- ..... .- . c u m u s m a - h = 100 mbar sf} 0.04 0.01 0.1 1 it : 0.01 " Figure 3-9: Three-phase oil relative permeability obtained by method 1 (solid lines) with their uncertainty limits (shaded area) and those obtained by method 3 (heavy dashed lines) with their uncertainty limits (light dashed lines). til-13nd: 80 In method 2, the base parameters were used to construct the curves of water saturation Saw against capillary pressure h for the air-water system using Equation 3-2. The 95% confidence limits obtained during regression were used as the uncertainty limits of the curve prediction. The same base parameters plus the measured interfacial tension were used to plot the air-oil Sf,lo and oil-water Saw curves using Equation 3-4. The uncertainty limits of these curves were determined based on the uncertainty in the base parameters and the scaling factors. The uncertainty limits of the base parameters were assumed to be identical to their confidence limits obtained in the parameter estimation procedure. The uncertainty limits of the scaling factors were determined by propagating the confidence limits of the interfacial tension values. The two-phase relative permeability curves k3: k3: and k3,},v were plotted employing the two-phase form of Equation 3-10. With this two-phase form, the terms of Stww and S3," in Equation 3-10 were substituted by SEW , Saw or 53° as necessary. Note that, the saturation 53‘” , Saw or 83° were substituted in the form of Equation 3-2 and Equation 3-4. As such the relative permeability functions were expressed in term of the capillary pressure, base parameters and interfacial tension. This expanded function was used to determined the relative permeability at a given capillary pressure. This function also served as the base function to propagate the error of the base parameters and the scaling factors in order to determine the uncertainty limits. Similar procedures were employed to construct the rest of the constitutive functions and to determine their uncertainty limits. 81 In general, method 3 and method 1 followed similar procedures to construct the constitutive curves. The base parameters used in method 3, however, were the air-oil parameters obtained from the individual regression while those used in method 1 were from the pooled regression. Each of the saturation functions Saw , S§° or Saw in method 1 has its own base parameter set. Therefore, these curves were constructed in a similar manner to those used in the case of 53‘” in method 2. As shown in the previous section, the three-phase oil saturation (SSW) and relative permeability (kfém) functions are formulated as fimctions of the oil-water and air-oil capillary pressures (how and hm). When how and ha0 are varied, a family of 53°“ or k3,?“ curves are obtained (Figures 3-5a, 3-5d, 3-6a, 3-6d and 3-7). Among these curves, two which allow the widest range of oil saturation change were selected for further analysis. The selected curves were referred to as vadose zone and saturated zone curves. The vadose zone curves were plotted by increasing hao while holding h0W constant at high capillary pressure (Figures 3-5b, 3-5e, 3-6b and 3-6e for SS” and Figures 3-8a, 3-8b, 3-9a, and 3-9b for 1:33“). These plots represent oil desaturation as the air saturation increased at constant low water saturation. In the field, this case might occur in the vadose zone when the water phase is at residual saturation and the oil phase percolates as it is driven by gravitational force. The saturated zone curves were plotted reversing the vadose zone situation, ie., how was increased but h,lo was held constant at 82 low capillary pressure (3-5c, 3-5f, 3-6c and 3-6f for SS” and Figures 3-8c, 3-8d, 3-9c, and 3-9d for k3,?" ). In this setting the oil saturation was increased as water desaturation occurs at constant high total liquid saturation. Such a situation may occur when the DNAPL penetrates below the capillary fringe zone and replaces the water volume in the soil pore. Table 3-2: Summary of the agreement between prediction outcomes by methods 2 or 3 to method 1 Ottawa Sand : Metea Soil Method 2 I Method 3 1 Method 2 } Method 3 Saw - i O i - i 0 53° 0 i- ' i- O i- ' SW 0 i- 0 i- O i- O ...‘l .......... -1 ........ 1. ....... 1. ....... 1. ........ kg? 0 | . I O | O ---------------------- ,__-___-,____-__,_____-_. km 0 1 . l O 1 O _JSL ___________________ +— _______ 1. _______ +— ........ kg: 0 . o . o . e “an; ------------------ ‘l" ------- 1" ------- '1‘ -------- Adamant“: ...... 0. ___,'____'___-I--_0____.',___9_--. S3” 1n sat’d zone 0 1|_ (9 r- 0 1I” @ £on9329532112--____°____,'___.'___-I_-_C3___,'_-__C?.___. kaow in sat’d zone 0 ' . ' Q ' G) Note: 0 poor agreement, no overlapping uncertainty limits in most of the range. G) fair agreement, overlapping uncertainty limits. 0 good agreement, the mean laid within the uncertainty limit of method 1. - no prediction required Table 3-2 summarizes the subjective interpretation of the agreement between method 1 and 2 as well as between method 1 and 3 as observed on the obtained curves. The agreement levels in this table were classified into three categories. A good agreement score was given for the cases where the mean responses of the prediction curves (method 2 or 3) were laid within the uncertainty limits of the reference curve (method 1). A fair 83 agreement score was given for the cases in which the mean of the prediction curve was slightly off of the uncertainty limits of the reference curve, however, their uncertainty limits overlapped each other in all of the observed ranges. If those uncertainty limits were off at some parts of the observed ranges, a poor agreement score was given. The Ottawa sand and Metea soil columns in Table 3-2 contrast the effect of soil texture. The columns labeled as Method 2 and Method 3 show the effects of using air-water and air-oil as the system basis for prediction, respectively. The arrangement of the constitutive variables listed in the rows of Table 3-2, from top to bottom, were sorted according to the order of the procedure followed in constructing the two- and three-phase k,-S-P relationships. The variables which were potentially affected more by the error in the base parameters propagated along the prediction were placed in the lower rows. The subsequent sections discuss the issue associated with soil texture, and are followed by discussion concerning the methods. The characteristic of the error propagation along the prediction procedure is discussed in the last sections. 3.4.1. Comparison Across the Soil Types Table 3-2 shows that, when the Metea soil was used as the porous medium, all of the obtained k,-S-P curves were in the categories of good or fair. When the Ottawa sand was used, however, the outcomes range from good to poor. It appears that the soil texture affects the prediction outcomes. The Metea soil has loamy sand texture while the size of the Ottawa sand grains are practically uniform. Study by McCuen et al. (1981), based on 84 1,085 samples measured by Rawls et al., (1976) and Holtan et al., (1968), found that the standard error of the squared root of the Brooks and Corey’s (1964) A. (which is equivalent to van Genuchten ’s parameter n) for sand was 0.039. Compared to 0.013 for loamy sand, that for the sand was 300% higher. The pooled regressions regressed the air-water, oil-water and air-oil simultaneously. As such, any factor associated with a two-phase type involved was lumped. Accordingly, any differences in the estimated parameters including their statistical properties should be related to the soil type. The entries under sub-heading Method 1 in Table 3-1 list the confidence limits of the P-S parameters estimated by the pooled regression for the investigated soil. It is observed in the table that the confidence limits of all parameters obtained in the Ottawa sand were in almost in all cases wider than those measured in the Metea soil. A wider confidence limit suggests a higher possibility of obtaining various outcomes when a number of data sets are used to estimate the parameters. This explaines the pattern exhibited by the scoring outcomes across the soil type in Table 3-2. In Ottawa sand cases, where the confidence limits of the base parameters used were relatively wide, the outcomes were varied in a wider range of outcomes: from good to poor. A favorable outcome may be obtained (Method 2-Ottawa sand in Table 3-1). Yet, there was also a risk of obtaining unfavorable outcome (Method 3-Ottawa sand in Table 3-1). In contrast, the outcomes of the Metea soil cases, where the confidence limits of the base parameters used were relatively narrow, were consistently favorable. This evidence suggests that when prediction via scaling theory was implemented in a uniform graded porous medium the obtained result may vary. At the 85 extreme case, an unfavorable result such as those observed in the case of Ottawa sand might be obtained. 3. 4. 2. Comparison Across the Methods Across the methods, as observed in Table 3-2, method 2 received a greater number of good scores compared to method 3. This tendency was consistently observed in both soils. This scoring outcome might suggest that method 2 is a better substitute for method 1 than is method 3. The major difference between method 2 and 3 was in the two- phase type from which the base parameters were determined. If it is true that the superiority of method 2 over method 3 is due to the use of the two-phase type, then a parallel tendency should also be observed in their base parameters. The following discussion inspects any possible correlation between the base parameters used in methods 2 and 3 and the superiority of method 2 over method 3. Except the scaling factor B, the other parameters S,, n and Claw (method 2) or Ciao (method 3) were inspected. At this stage, each of these parameters is assumed to have potential for influencing the difference in the scores across the two methods. As shown in Equation 3-5 and Equation 3-6 the scaling factors B were involved in the scaling procedure. These factors, however, were determined using the interfacial tension data and the same interfacial tension data was used by the two methods. Based on this, it is unlikely that B has potential to influence the difference between methods. 86 Table 3-3: Discrepancies between the base parameters used in method 2 and 3 relative to those used in method 1. Base Ottawa sand Metea soil Parameters Method 2 Method 3 Method 2 Method 3 S, 27 % 20 % 5 % 2 % a“, 1 % - 4 % - a“, — 26 % - 8 % n 5 % 31 % 9 % 7 % Table 3-3 shows the discrepancies of the base parameters S,, n and (law or 0:30 employed by the two methods relative to the base parameters used in method 1. The magnitude of these descrepancies are regarded as the level of the error of those parameters. As observed in this table, the discrepancies in the parameters S, and n were not consistent with the scoring outcome observed in Table 3-2. For example, while method 2 received a better score than method 3, the errors of the S, employed in method 2 was higher than those employed in method 3. That S, and n cannot be associated to the two-phase type involved appears to be in agreement with the assumption which states that these parameters are porous medium dependent parameters rather than the two-phase type (Lenhard and Parker, 1987). Note that, methods 2 and 3 estimate the P-S parameters using the individual regression procedure. This procedure did not pool the data according to the soil type. The S, and n used in method 2 were obtained from air-water measurements while those used in method 3 were from air-oil measurements. To permit more general inspection, S, and n estimated by the individual regression upon each of the air-water, air-oil and oil-water data were 87 l 1- n- r H ... 1 i .... l l' J 1 l7 ... l L r7 5 " . 1 i fi u r c i *- I I \O V ("I N ll 1 n l .. I u 7 1 2 ~ —u p l n I _ l H l I I I I In ”I "1 N. "I Q d c o c c o Metea Soil Air-Oil Metea Soil Oil-Water Metea Soil Air-Water Metea Soil Air-Oil Metea Soil Oil-Water Metea Soil Air-Water Ottawa Sand Air-Oil Ottawa Sand Oil-Water Ottawa Sand Air-Water Ottawa Sand Air-Oil Ottawa Sand Oil-Water Ottawa Sand Air-Water Figure 3-10: Parameters S, and n estimated by the individual regression with their 95% confidence limits I-rtedr 88 plotted in Figure 3-10. With the individual regression procedure, neither their independence to the two-phase type nor their correlation to the soil type was artificially imposed. However, as observed in Figure 3-10, except in the case of "are of the Ottawa sand, S, and n appear to be more consistent across the soil types than across the methods. Based on all of the mentioned evidence, it is concluded that S, and n are independent of the two-phase type used. As such, the outcome across the two methods observed in Table 3-2 was unlikely to be the signature of the errors of S, and n. Unfavorable scores in the cases of method 2 using Ottawa sand migh be strongly influenced by the error in nao. ‘ However, it is viewed as more characteristic of the soil type rather than the methods. The discrepancies in the parameter at in Table 3-3, on the other hand, appear to be correlated to the differences in the scores obtained by the two methods. Method 2, which received overall good results had or discrepancies of only 1% and 4%. Method 3, however, received overall poorer results and had at discrepancies of 26% and 8%. It appears that good result might be correlated to low discrepancies, whereas poorer results to high discrepancies. As concluded in the previous section, S, and n appeared to be porous medium dependent parameters. The pooled regression accommodates this condition. Thus, oraw cow and am obtained via pooled regression are regarded as the most representative on. As such, any discrepancies relative to these on can be used to measure the level of errors 89 embodied in or. Additionally, the width of the confidence limits of the or obtained by the pooled regression indicates the range in which different outcomes from different data sets might be obtained. The wider these limits the more the risk of obtaining higher discrepancies. In turn it also indicates higher risk of obtaining stronger errors. The parameters listed under sub-heading Method 1 in Table 3-1 were obtained via pooled regression. An inspection on a values in this list, shows that the order of the magnitude of a was aao>aow>aaw and the width of the confidence limits of those obtained by the pooled regression followed the same order. These patterns were consistently observed in both soils. The observed pattern suggest that the higher the magnitude of a the lower its precision. A study by van Genuchten and Nielsen (1985) postulated that or is inversely proportional to displacement capillary pressure(ade '1). Utilizing this relationship, the observed pattern indicated that low Pd value possesses low accuracy. Lenhard and Parker (1987) observed similar problems. Their P-S measurement using benzyl alcohol as the oil phase covering the lowest range of capillary pressures among the measured two-phase systems, suffered the highest uncertainties. In the investigated cases, the aaw in method 2 was equivalent to Pd= 43.5 mbars for the Metea soil and 17.1 for the Ottawa sand. In method 3, 030 was equivalent to 18.7 (Metea soil) and 7.5 mbars (Ottawa sand). In order to anticipate 200 mbar maximum capillary pressure measurement in the air-oil system of the Ottawa sand, a 5 psi (=345 mbar) rated pressure transducer membrane was used. While it gave a good resolution for 90 low pressure measurement, it lacked accuracy due to a zero shift problem. On the average, the calibration in the air-oil measurement setting indicated a $0.7 mbar standard error of estimate. When this standard error of estimate is expressed as 95% confidence limits it is equivalent to approximately 20% of the Pd of the air-oil system in the case of Ottawa sand. For the air-water system, on the other hand, it equivalent to only 3.5% of Pd. Based on this inspection, it is concluded that am was less reliable compared to aaw. Method 3 used aao while method 2 used oraw. When identical S, and n were employed in the two methods and since (1a,, was less reliable than oraw, the chance of obtaining bias in the prediction outcome is higher in method 3 than in method 2. As observed in Table 3-2, method 2 performed reasonably well not only when a well-graded soil such as Metea soil was used but also when the poorly graded Ottawa sand was used. Method 3 appears to be acceptable only when Metea soil was used. The source of discrepancies between methods 1 and 2 or methods 1 and 3 was due to the discrepancies in the P-S parameters used. When method 3 was applied using Ottawa sand, the base parameters used suffered the combination of two problems; less precision in the parameters S, and n due to the texture of the soil and low accuracy in aao due to difficulty in measuring low capillary pressure. As a result, the parameters determined solely from the air-oil data measured in a uniform grain soil might not be representative enough to be used as the scaling reference compared to those derived from the air-water data. 91 As presented in the previous section the two- and three-phase oil saturation S,, and relative permeability k,0 are functions of the two-phase oil saturation in an air-oil system (S,,) and water saturation in an oil-water system (Sow). In method 2 both Sa0 and SW are prediction curves while in method 3 only S0,, is a prediction curve. Thus, in obtaining So and k,,,, method 2 would propagate any error embodied in the base parameters more than would method 3. Based on this, prediction via method 3 should not be inferior when compared to the prediction via method 2. The scoring outcome on Table 3-2, however, shows that method 2 was better than method 3. As concluded previously, method 3 has a problem in determining the parameter aw. Apparently, the difficulty of measuring low capillary pressure has a stronger effect upon the outcome than when this additional source of error is induced. 3. 4.3. Efl'ect of Propagated Error on the Constitutive Variables Parker et al. (1987) postulated that prediction of the relative permeability function might be sensitive to any disturbance in the base P-S parameters. The arrangement of the rows in Table 3-2, as mentioned previously, roughly reflects the sequential steps in constructing the two- and three-phase constitutive relationship curves. In turn, it also reflects the way the errors in the base parameters and scaling factors were propagated. The constitutive variables in the lower rows were functions of the variables in the upper rows. Therefore, the lower row variables possessed the uncertainty of the variables in the upper rows. Thus, the lower row variables possessed greater uncertainties 92 as the error that embodied in the base parameters and the scaling factors were propagated further. The columns under the Metea Soil heading in Table 3-2 show the above trend. The upper rows, up to variable k3: , were scored in the category good. Below this row, some scores of fair were observed. This pattern indicates that the errors propagated up to variable hf: produced only insignificant bias. Further propagation slightly increased the bias. The bias resulting from this increment at the end point (the three-phase oil relative permeability), however, is still within tolerable limits. This was possible as the error embodied in the original bias was relatively small. Therefore, it follows that the magnitude of error of the base parameters as tabulated in the pertaining columns in Table 3-3 can be regarded as small. It is interesting to observe that 27% error in the S,, as shown in the column Ottawa Sand - Method 2 in Table 3-3, yields insignificant bias in all of the prediction curves. This is only possible, however, if the accuracy of the other parameter was high. Nevertheless, this particular case suggests the insensitivity of the employed prediction methods employed to the parameter S,. It is speculated that the bias in the prediction curve associated with error in this parameter is localized in the vicinity of the low saturation region. On the other hand, even with parameter estimation based on measured data the uncertainty in this region is high due to the non-linearity of the constitutive relation firnction. The errors of the base parameters in the case of Ottawa sand - Method 3 (Table 3-2) were high. The predictions employing these parameters, as shown in Table 3-2, were 93 unfavorable. Despite the significant discrepancy observed in the prediction curves the errors of the base parameters in this particular case was not clearly evident in the P-S curves plot. Figure 3-11 was plotted to demonstrate this issue. This figure plots the measured data and the fitted P-S curves on these data via the pooled regression (dashed line) superimposed by the curves obtained via the individual regressions (solid line). The parameters used by the pooled regression curves were the parameters employed in method 1. Those used by the individual regression curves were the base parameters used in methods 2 or 3. Thus, the parameter discrepancies listed in Table 3-3 were the discrepancies of the parameters used by the dashed and shaded areas. In Figure 3-11 the discrepancy between the solid and dashed curves in the case of the air-oil system in Ottawa sand appears to be comparable, for example, to the case of the air-water system in Metea soil. As a P-S plot, error as high as 20%, 26% and 31% in parameters S,, n and a appears similar to errors of 5%, 4% and 9% (Table 3-3). If these errors were included in the scaling procedure, however, a significant difference in outcomes would be clearly observed. In the same example, the air-oil and air-water systems were scaled to predict the oil-water P—S curve. The resulting curves were the plots in Figure 3-2a and Figure 3-lc, respectively. It is observed that the plot in Figure 3-2a was significantly worse than that in Figure 3-lc. Note that, the errors of the scaling factors involved in constructing Figure 3-2a was only 0.02% while in Figure 3-1c it was 9.5%. Therefore, it is believed that the errors observed in Figure 3-2a and Figure 3-1c were strongly affected by the errors in their base parameters rather than to the errors in their scaling factors. 94 it... see 885 aE_ 8838... $3 hoe a? $8: Boa anemone nausea 2: B engage ones one and: Bean awe aha: convenes $3 hon“ 5S, Anon: @053 .0605 eofimawoc ©2er 05 .3 8583 morn—o nougaméuamoa omega; ”:-m 0.53m a... 3?. saw a... ... no a... 2 S 3 o... z. 2 2.. 2 ... . 2... . 2 , . .2 . mu .8 .a a .3 .3 .8 . a . 8 i a . 8. .. on u: r s "I n: u... ) ) ) .8 m .2. m m . a . .... ...... ,8 ( .R (\ EI\ . e. _ - 8. , e8 . 8. gag-=0 " . a: . a “ new $802 “a. w . on a . 8. -an 8a .5 tea he. sum 3...“. 2 2 S no no I. 2 .... no. I n. 2 .... no 4... 2 o... «9 ... q. 3 2. a. .... 2 a... «so .2 .3 . a ”I u- u- ) ) I 3 ) m m m a. a. a m. m. m ( . 8. 93531.": Us“ 3% r 8a 80 .9 95 In the actual implementation of the scaling theory, only one system is required to be measured. When data is available from only one system, the goodness of the obtained base parameters can only be verified through an inspection of the agreement upon the best estimated P-S curve against the data points. As demonstrated above, this technique might not be sufficient to reveal even errors as high as 20%, 26% and 31% in parameters S,, n and on. Based on this, in order to increase the possibility of obtaining more representative parameters, it is necessary to replicate the sample. 3. 5. SUMMAR YAND CONCLUSION Parker et al. (1987) proposed a predictive method employing scaling theory. Using their theory, a complete set of two- and three-phase k,-S-P relationships can be constructed based only on one set of P-S parameters, referred to as the base parameters, measured in one two-phase system plus the pertaining interfacial tension data. Two alternatives to implement this theory were investigated, referred to as method 2 and 3 respectively. Method 2 uses the air-water system to determined the base parameters while Method 3 uses the air-oil system. Method 2 was regarded as the most efficient alternative in the case where air-water P-S data are available. The base parameters S,, n, (law were determined from the available data and the parameters (low and cm were predicted via interfacial tension measurements. Method 3 was similar to method 2, but the S,, n, Ciao were determined 96 from air-oil P-S data while the (raw and aow were predicted. This alternative might be suitable in cases where sufficient data are not available, yet more efficient procedures than those used in method 2 are desired. These methods were tested to predict the two- and three phase k,-S-P relationships using two contrasting soils: uniform grain Ottawa sand and loamy sand Metea. As the reference to justify the performance of methods 2 and 3, the obtained prediction curves were compared to those obtained via method 1 which is a predictive method without employing the scaling theory. When the Metea soil was used as the porous medium, the prediction of the k,-S-P curves yielded a favorable outcome. For the cases where Ottawa sand was used the outcomes vary. At the extreme case, an unfavorable result such as those observed in the case of Ottawa sand with method 3 might be obtained. Inspection of the statistical properties of the base parameters, suggests that determination of these parameters in a poorly graded soil, such as Ottawa sand, might result in a high variation. In turn, the risk of obtaining less representative parameters is also high. The investigated cases suggest that the parameters S, and n were independent of the two-phase system used. The deference between the outcomes obtained via methods 2 and 3 was associated with the errors in their values of the parameter 0:. Determination of or in the air-water system was more reliable than in the air-oil system. This was due to the 97 difficulty of measuring low capillary pressure. As a result, prediction by method 2 which is based on air-water data appears to be superior over method 3. When method 3 was applied using Ottawa sand, the base parameters used suffered the combination of two problems; less precision in parameters S, and n due to the texture of the soil and low accuracy in am due to the difficulty in measuring low capillary pressure. As a result, the parameters determined solely from the air-oil data measured in a uniform grain soil might not be representative enough to be used as the scaling reference compared to those derived from the air-water data. A less representative parameters set may not deteriorate in agreement between the data points and the regressed P-S curve, however, it may cause a problem if used as the basis of prediction. Hence, the use of method 3 should be limited to a well-graded soil. For a uniform grain type soil, it is more rigorous to employ method 2 which gives better assurances in providing representative parameters. Furthermore, the issue of difficulty in measuring low capillary pressure inherent in method 3 has a more significant effect than the error propagated by the additional step required in method 2. BIBLIOGRAPHY Anderson, W.G. 1987. Wettability literature survey Part 4: Effects of wettability on capilary pressure. J. of Petrol. Technol, 39, 10:1283-1300. Azis, K., and A. Settari. 1979 Petroleum Reservoir simulation. Applied Science Publisher. London Bear, J. 1982. 1972. Dynamics of Porous Media. American Elsevier Publishing Co. Inc. Beck, J .V., and KJ. Arnold. 1976. Parameter Estimation in Engineering and Science. 501 pp. John Wiley and Son Inc. New York. Bethel, FT, and J .C. Calhoun. 1953. Capilary desaturation in unconsolidated beads. Petroleum Transaction, [98:197-202 Brakensiek, D.L., RL. Engleman, W.J. Rawls. 1981. Variation within texture classes of soil water parameters, Trans, ASAE, 24,2:335-339. Brooks, RH, and AT Corey. 1964. Hydraulic properties of porous media. Hydrology. papers no. 3., Colorado State University., Forth Collins, Colorado. Burdine, NT. 1953. Relative permeability calculations from pore size distribution data. Petroleum Transaction, AIME, 198:71-78 Corey, AT 1986. Mechanics of immiscible fluids in porous media Water Resource Publ.Chelsea. Michigan. Corey, A.T., and RH. Brooks. 1975. Drainage characteristics of soils. Soil Sci. Soc. Amer. Proc., 392251-255 Corey, A.T., C.H. Rathjens, and J.H. Anderson. 1956. Technical note: Three-phase relative permeability. Petroleum Transaction, AIME, 207:349-351 Cosby, B.J., G.M. Homberger, R.B. Clapp, and TR. Ginn. 1984. A statistical exploration of the relationships of soil moisture characteritics to the physical properties of soils. Water Resour. Res, 20, 6:682-690. Demond, A.H., and P.V. Roberts. 1991. Effect of interfacial forces on two-phases capillary pressure-saturation relationships. Water Resour. Res, 27, 3 :423-437 Dumore, J.M., and RS. Schols. 1974. Drainage capillary presure functions and the influence of connate water. Soc. of Petrol. Eng. J., 14, 5 :437-444 Elrick, D.E. J .H. Scandrett, and EB. Miller. 1959. Test of capillary flow scaling. Soil Sci. Soc. Amer. Proc., 23:329-332 Gee, G.W., C.T. Kincaid, R.J. Lenhard, and CS. Simmons.1992 Recent studies of flow and transport in the vadose zone. Kaluarachchi, J .J ., and J .C. Parker. 1989. An efficient finite element method for modeling multiphase flow. Water Resour. Res, 25, 1:43-54 98 99 Kia, SR, and A. Abdul. 1990. Retention of diesel fuel in aquifer material. J. Hydraul. Eng, 116, 7:881-894 Klute, A. 1986. Water retention laboratory methods; Methods of soil analysis, Part I: Physical and mineralogical methods, Amer. Soc. of Agron, Madison, 26:635-661 Klute, A, and J .E. Wilkinson. 1958. Some tests of the similar media concept of capillary flow: 1. Reduced capillary conductivity and moisture characteristic data. Soil Sci. Soc. Amer. Proc., 22:278-281 Kool, J .B., J .C. Parker, and M.T. van Genuchten. 1985. Determining soil hydraulic properties from one-step outflow experiments by parameter estimation: 1. Theory and numerical studies. Soil Sci. Soc. Am. J., 49:1348-1354. Lenhard, R.J., and J .C. Parker. 1987. Measurement and prediction of saturation-pressure relationships in three-phase porous media systems. J. Contam. Hydrol, 1:407-424. Lenhard, R.J., J .C. Parker, and 1.1. Kaluarachchi. 1989. A model for hysteretic constitutive relations governing multiphase flow. 3. Refinements and numerical simulations. Water Resour. Res. J., 25, 7:1727-1736. Lenhard, R.J., J.C. Parker. 1988. Experimental validation of the theory of extending two- phase saturation-pressure relations to three-phase fluid phase systems for monotonic drainage paths. Water Resour. Res. J., 24, 3:373-380. Lenhard, R1, J.C. Parker. 1988a. Experimental validation of the theory of extending two-phase saturation-pressure relations to three-phase fluid phase systems for monotonic drainage paths. Water Resour. Res. J., 24, 3 :373-380. Lenhard, R.J., J.H. Dane, J.C. Parker, and JJ. Kaluarachchi. 1988b. Measurement and simulation of one dimensional transient three-phase flow for monotonic liquid drainage. Water Resour. Res. J., 24, 6:853-863. Leverrett, MC. 1942. Capillary behavior in porous solids. Trans. Amer. Ins. of Mining Eng, 142:152-169 Luckner, L., M.Th. van Genuchten, and DR. Nielssen. 1989. A consisten set of parametric models for the two-phase flow of immicible fluids in the suburface. Water Resour. Res. J., 25, 10:2187-2193. McCueen R.H., W.J. Rawls, and TL. Brakensiek. 1981. Statistical analysis of the Brooks-Corey and the Green-Ampt parameters across soil textures. Water Resour. Res. J., 17, 4:1005-1013. Mercer, J .W., and RM. Cohen.l990. A riview of immiscible fluids in the subsurface: properties, models, characteristization and remediation. J. of Contam. Hydrol, 6: 107-1 63 Miller, BE, and RD. Miller. 1956. Physical theory for capillary flow phenomena. J. of Applied Physics, 27, 4:324-332 Mishra, S., and J.C. Parker. 1989. Parameter estimation for coupled unsaturated flow and transport. Water Resour. Res, 25, 3: 385-396 100 Morrow , N. R. 1976. Capillary pressure correlation for unifomrly wetted porous media. J. Can. Pet. T echnol., 15, 4:49-69 Mualem, Y. 1976. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. J., 12, 3 :513-521. Parker, J .C. 1986. Hydrostatic of porous media. Soil physical chemistry, Ed. D.L. Sparks, CBC Press. Parker, J .C., and RJ. Lenhard. 1987. A model for hysteretic constitutive relations governing multiphase flow. I. Saturation-pressure relations. Water Resour. Res. J., 23, 12:2187-2196 Parker, J.C., J.B. K001, and M.T. van Genuchten. 1985. Determining soil hydraulic properties from one-step outflow experiments by parameter estimation: 11. Experimental studies. Soil Sci. Soc. Am. J., 49:1354-1359. Parker, J .C., R.J. Lenhard, and T. Kuppusamy. 1987. A parametric model for constitutive properties governing multiphase flow in porous media. Water Resour. Res. J., 23, 4:618-624 Parker, J .C., R.J. Lenhard, and T. Kuppusamy. 1987a. A parametric model for constitutive properties governing multiphase flow in porous media. Water Resour. Res. J., 23, 4:618—624 Pinder, G.F. 1982. Groundwater Hydrology - Research Needs for the Next Decade. Fundamental Research Needs for Water and Waste Water Systems, AEEP/N SF Conference, M.S. Switzenbaum, ed., Arlington, Virginia Purcell, W.R. 1949. Capillary pressure - their measurement using mercury and the calculation of permeability therefrom. Petroleum Trans, TP 25 44, 186 :39-48 Rawls, W.J., D.L. Brakensiek, and K.E. Saxton. 1981 Soil water characteristics. Paper no. 81-2510, ASCE., St. Joseph, MI. Rapoport, LA. 1955. Scaling laws for use in design and operation of water-oil flow models. Petroleum Trans, AIME 204: 143-150 Russo, D. 1988. Determining soil hydraulic properties by parameter estimation : On the selection of a model for the hydraulic properties. Water Resour. Res. J., 24, 3:453- 459. Schwille, F. 1967. Petroleum Contamination of the Subsoil - A Hydrological Problem. Joint Problem of Oil and Water Industries, P. Hepple, ed., Elsevier Publ. Co., New York. Seber, G.A.F, and C.J. Wild. 1988. Nonlinear Regression. 768 pp. John Weley & Sons Inc. New York. Smiles, D.E., G. Vachaud, M. Vauclin. 1971. A test of the uniqueness of the soil moisture characteristic during transient, nonhysteretic flow of water in a rigid soil. Soil Sci. Soc. Amer. Proc., 35 :534-539 101 Smith, W.O. 1933. Minimum capillary rise in an ideal uniform soil. Physics, 4: 1 84-193 Smith, W.O. 1933. The final distribution of retained liquid in an ideal uniform soil. Physics, 4 :425-438 Stone, H.L. 1973. Probability model for estimating three-phase relative permeability. J. Petrol. T echnol., 20:214-218 Su, C., and RH. Brooks. 1980. Water retention measurement for soils. J. Irrig. Drain. Div., 105-112. Topp, G.C., A. Klute, and D. B. Peters. 1967. Comparison of water content-pressure head data obtained by equilibrium, steady-state, and unsteady-state methods. Soil Sci. Soc. Amer. Proc., 31:312-314. van Genuchten, M.T. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J., 44 :892-898. van Genuchten, M.T., and DR. Nielsen. 1985. On describing and predicting the hydraulic properties of unsaturated soils. Annales Geophysicae., 3, 5 :615-628 Warrick, A.W., J .J . Mullen, and DR Nielsen. 1977. Scaling field-measured soil hydraulic properties uing a similar media concept. Water Resour. Res. J., 13, 2:355- 362 White, N.F., D.K. Sunada, H.R. Duke, and A.T. Corey. 1970b. Boundary effects in desaturation of porous media. Soil Science, 113:7-12. White, N.F., H.R. Duke, D.K. Sunada, and A.T. Corey. 1970a. Physics of desaturation in porous materials. J. Irrig. Drain. Div., 165-190. Wilkinson, J .E., and A. Klute. 1959. Some tests of the similar media concept of ccapillary flow: 11. Flow system data. Soil Science, 23 :434-437 APPENDICES A. PRESSURE-SATURATION MEASUREMENT SINIULATION A. 1. Model Conceptualization. A.2. Simulation Flowchart A.3. Program Listing in Excel Worksheet. 102 103 A1 MODEL CONCEPTUALIZATION x=3.0 r 81.315 “r 2.45 L 13am: 7:. 2.95-- o n/f/l, 175 a? / 50,533,... 1) """"""""""""" 4 105 M 0" .. ,0. f, ........................... ............ ............... ............ .................... ................... ................... Boundary condition: a) SE method b)PEnl:thod Initial oonditiom Hydrostatic profile with hxgo==h (20 —) qx=3=0 ’20 -’ hx=0=hprescnbed I20 -) qx=3=0 0tdelat "" 9x=0=0 presa'ibed 104 A2, SMULATION Specify Simulation Parameters FLOWCHART Construct Initial h Profile / Specify Initial Guess of At 1 - L 80) k 0) do) / l e / t— - t+ At ] I At 1 AS— - — 21x q(t) A”? / S(t+At)= S(t+At) t h(t + At): f(S(t + N» ha) = 110+ At) V(t) = V(t + At) AI=AIX2 Yes EITOI'inhU‘l‘A!) A , Switch N0 "'— ac. J .3 Yes 11 We +111): V(:)+§(x= 0)-At/ or Irregular h profile h‘+1(t+At)= h'(t+At)- it? at S(t + At) = f(hi+l(t + 111)) k, (t + At) = f(S(t + At)) q(t + At) = f(S(t + At), k,(t + 111)) N End 0 of simulation ? owes 105 an 8 65392:? anoint: «a 8 «Ego 1.250 5 3383309390.» none.» on 8. nee—gen as as one 89688.th "258.com cu o8 3888 5. 9: nos: .58 nu . E3 2%” . . 8» BEImIo.§I use: 83m on 3. man on man mg no new mu 0 I I I I I + o Isseom...8aee§o..?oam§u nouns. 3 rz...>..o~en_n an \ «a md i 1 md :82: .35ch o «.83. 3 58%. ...-80.9 I «88. on \ 8 Lone}. . . om "8.83 . . e i - P 3.2092» Louis I Ran queoI§< . Z 2 . Me I- 1 mm. 95 Bates. : .1 ( can. ..ro No "Item " . as». .re no «Elam “ 923.. as I N 1- -- N gonad no :. so: Ended .- . «88.88983 to. I on t .- ....w 8.6.38.8 I. 92808.» 58.95 c .=Om-o_:._n_ m. M w n n v M II: x 1 I ... ceases. 3 a T E95203 III: 11 ii _qmoxm E 022.24. 509E .m.< [III _ r 106 n no 3 no 3 no an we 3 :58. rsgifiohmeén ozhsaévohgaé. on : Ears—328983;." aghafighamocu 8 :58- rsgrpmokmelw E§§§9§83r B Izaarsérfiohmeén Q§§§9§327 8 camiégrfimokmeén Shama+~v9§39=n an : agrsgfifiohfircu $.5§lvo.§337 3 : amaisamrfiohm? :u ox.“ 5863..“ .393" no cami<6§2umobame+sn Bu.»3.80..~.33r «a Agavggrcmogeiu ShamQSohamvsu E gunfirsamiomoamaéu Bhaazhohzmvcn 3 car» cart at 3 eon hv I8m9§8mw$en I.5m3<§32998m9+5.633788» «2833388 E9??? 0 3 6395:3199" 233.6.329653?5.639983" 5588355.” $33? no 3 ARmWIXSmgor c53:6.3930937fag—$83., 3905885..“ §¥<+9$u 0.0 3 §m¢§8m§7 c53.5.3233?5.639563? 358.585 §9<+uvfim-wg.fin8§m$u .133...$o$o§§0§§+g.633:.. «N .3....miagnbgozame—u :.avm-.3§mo$o:.¢m0§3?50.8.33? r4. .vvmgvxougéugfiozggu :-._Eg.>vmo$o....80b;m?90.§»m9:7 on .Qmfiexmrzégozagn .19mg.xfio.umov...~mob.3?80.§§Eu a. .8089§E.t.¢.88»§m$n ....flmgvxgfiovz08939896839" 3 .83851903028»... .1838:.830:..80.n_3»£8.§83:u 2 2 same can 3 V v £39mvxnc.ur.r8»8§m$u 8...§v0¢o:.$0.§mnz.u u. 83¢9§I¢E§8an$u ..Amvmgvxovog...90u “. gig—”3.5.8828“... :..mvmgivoevosevou Y 3m¢miofi$vgo§3$u ..Avvmgivvoéougn 2. .g-~vmx§$r88»kno$u ..§.gi99uwo:§0u .. .g-_vm.>fi.t...m~»8§m$n .igénxgéofvou = :vmgzufivgonkamsu :.cvmdvngvodvozgu, .. Sggumvgkamgu .Iovmfimxovogzfion o fingifii.§8&8$u .153§80$o=.80u .— .Bm.§xI&r§o€3$u .....Bmsnmxsqaefmon : o A0158 . ea. 8.. . 0 u. lll :.....Ectt q. :..... . _c 112 _M E: 8- on 3 3 8 ma 3 8% 9.6.32.3.3%....983553Emm. 5.56.9898?“ no .8...&2..8..&E.~. «....ogacéamer... . £13.. EFFIEQ§§+E§33u «a gutabaanréhn ~<.6.8&._¢5mer.... . . rs Inn .Y...m9§6m.8mzbam»n 3 3.3.5.9“? 9.633253»: 5. , , r}... H-.. 3.50.98.898508? 8 .8u.+8".5$.8“..~n N..o.30._¢_umer..... , . .l .. - -.. ......m.9.~..8w+8m:bam»uumfl. .8“.+8u..auu.8u.~. 2.9822983... . ii........uF..IEEEBmiomxbamnn on .Bmkouiotfihu ~<..o.»me..¢.3er.... . , . . .... ..Irméficbméomzamou no ...Eéogvkougmhu 9.6.82.358?...3.890.3aiusomfiomz.simiusomémvrsau on 88.80.08??? 9.6.30.388?.....ofieaamarfifiméomfv:.m.9.~..mom§om=.8$u no .Buéoufifiguhu 9.6.32.3.32:...Sameaamaiuimémz.FIVm9§$m+8mzbnSu 3 o o no 6:31.35... 63.5.... no 8 8“.» 9.6.82.7939......33335355-5.5.9.2983» on .8uéuvaomhuuu 9.6.82.gamer.....o.33.o.32<.u.2m+fimz-.I..m.9.u.2mém=.8mnn an .Ruébatkuhu «....oama...oam§<..... . , 3.35:3...:....m.9.~..2m+kmsb»non as .Rufititkthn 929833.331:Bamavaamuiufmémx.Erméfifmémsaan k Euétvfimbthn 2.930.3.33:....¢_3ea$m?.§m.mm¥:....m.9.u.2.m+m..mzbanan on fiufiuitbnuhn 9.6.82.3.83:...3.335.»..nrfifimxfimvv_Y...m.9.e.m..mi.m:b8¢n on .Ruiufiufituu 2.932.Samar...3.823303%???va2.5.09.0.2méu39m? I .vtéuflvthu.~.u .2.. E $5 . in? «5.5. .883"... .E.E.E..E.~. “.633.8.33:...3.395.3«E5ém....Y5.09.55.20.38" 2. no 0 _ N N :..“... 5.53122: __< EEuEE. Sc Bali! gaggicnu .§«¢§2.2C_P¢uutlnp n. 113 0 SF 2: an 3 .3 3 an 6.293380605323890." Soh.»n?$a8..~.83.u 80b.3»63096.83.n 3 ....ggfigegag.» roofinmnéxorfiégr 80658.9»9633..." 8 1‘ 2-6333...8.-$...&%u Boa—303386.33.“ .8b.3»$809.983.n no ....v3§<»...$_.5....§.u 80h.8»§8..~.33.u 805333.933." 3 111 ..-.§§$.>S_8_..bmxu Sofiaaéaonhgfr 80932309533.“ on I 2.3.2333...8.$....qu .8&.§+.38..~.83.u 80939309658...» 8 223535...8..8....8:« 80h.8»6309.§33.u 50989989583." 8 2-6333...8..§..Lo:u Bohgafiuahfiz.“ 80.n.8?8»0968mn...u .3 ....§§S.EQ$..§u v8§3$§3§83r m8h§9§09§837 on ...gbflfiggflsmxu 80h.»a?§8..~.33.u «8933809....33.» on o 3 6:63;... .o.%+....d>a..E 6E?m.fl9§n.. no no a .o..$§.80-§...§..on88».u__u Refinmnévoxusmfr 80939906289...» 3 ....93:8:F8Eu wa—oéamaaxoahgwsu 865.8»6392933." v2 2§$§<§82§§§Eu 8_o&.3n§oa.a3m3r 36.namuggn8mvsu a: 11 2.§$¢»5:§_.5:F§:u «u—ohrwmagofinamivn 820bfinior§8ain «3 2-3a.3§29.8§85u §o.~_3a§ou.§337 S.oh3$§.§895n 3. 2.2332332228282382? 8653??»82833.» 66b§$§ofiaamnéu 92 22359392827822.8277 8553338259.? 86982332535“ a2 2-29ogfiiszsozfmiu 86.2»3é395893u 8693aghamficu o3 2-28.35.29.82335u 8553958238.? 8209893096835" n3 2.A§<».~.8ulx Eta! hnN 0 n. 128 6; _3figgvgfvuxdu83nrv :.sxgéévuévgrmvnxu ~v~o&.3a§39§33r 98h.3n§.§3o.cu 3 :.ggéxmxiusbnxu Qwofiomnaxofifiamvcu $85.3»3gn337 3 Aligfirfigiusgxa 9853388339? 9~ob.3»§.§39:u ow . :-gggfiivvgévusfiuzn 3~o.~_3?9»8.§337 9~ob.3»§§.§33r 3 . :-gifixgwgsfiuxu 9853233539:- vvuoéamugkflamaéu 3 :.A.3 ..3 2: 3 8 no 8 an 3 «a §§$.§F§.ar Rag—oaoaiam? 632:2... maggoainau, «a §S§S§NF3~§u 3369093358» 633-353878? fiaofigooinnau 8 .§<§35~.785u 502.35.39.53- ?»figfixgv 6...- 33.893358", 3 .3533.>8_7-§u 83385630233»- A3§§i§r83u 33.89.039.53- 3 .333359293- magma—939.53- §<23<§87§7 33.80.38.58- 3 saggiflvgau Rauncfioénoekm? :..»gfiigr §§B.Fv»oa.haoan B AS§S.>8.3:2- 882853359» Aggioolo? Raggoekmau on .agfixat3zau 883338..»{39 $335969..- Rnoefioaosknan no .8»?§3§:7t5u 833953.353- .Baggafigovmos.» goggaoaiam? 3 CRS¢RSE:3:JH 8853.53.53? 983383.373...- oaaoaaoogfloninmnu no «a giant-$3 §8uémoh 58591.8 goat-$3... 3 8 exgfivxnzazau 832.95.333.33? 63.935078...- gflogvfiau an 633$...«:3:2u 833.53.33.58- .Ogsgau §§Ro¢oi3an 2 .§5.>:.3:_..u gonzo—Egonfnmou 333579.... $8503.58- k .§§g$.>o:.,8:.r «803853353- Svgvsxfivkan 5032030953? 2 335.35.558.3- 882.85.39E3au .95v<.>k723- 502653.53- an :Xgfiis$§3u 53:20-39Z3? 335.653... 302393.33»- I 6338.359783- gnlofigoeknm? 3.3%??? acknogfamnu H ill- .3 mg. r-Ebll' 3 (:..-1|! .898... £379... TE ghb. .§$~8$.>8$35. 839.35.839.53»- .Qhfizfivfi? 833.589.53- E j- ; $978... .- 8» 85.5.9 .3». ngmavxsr gaseofaufi .a. 9.2.. 9.3.. no. Emmi tar?“ z .2 4 _ z 130 131 132 IterNum=O TRUE SET. SET. IF ENDJF SET. IF SET. ENDJF E sht'! sht'! , sht‘! sht‘! sht‘!l41 sht'! sht'! sht'! sht‘! sht’! 133 sht'! sht' sht'!C sht'EC oulTime” . sht'!$C$26 TomlTimc>tPfim+l 0.0000001 1 5 sht'!$C$26 sht’! SCSB l 5 5111'! sht'!$C $76 sht’!$D$76 sht'!$C$74 sht'! SC$73 sht'! SD$73 sht'ESDS72 sh1'2SCS71 134 A FORML . 5. CELLXL ’rk sht'!$D$71 CELLXL 19.1.1 FORMI . 5, CELLXI. rk sht'!$C$70).()FFSET([CELL.XLWIRauhESAS6.RecCoum.20.l.1 ET. 5, CELL.XLW]er sht'!SDS70).OFFSET([CELLXL SSA56.RecCounL21.l.1 FORMLMGETCELUSJCELLXLWler sht'!SCSJO).OFFSET([CEU.XLW]Result!$AS6.RecCount,22,l,l FORNILWGETflI-ZLIAS, CELLXL rk sht'ESCSBZ CELLXL“ 'SA$6,RecCounL23,l,1 SETN ,RecCounHl SETN ° +1 Ime*0.0000001 1 =END RecordlnitCond =F =F =F =F =F =F GET. sht‘!$C$31 sht' sht' sht‘! sht'! =FORMU . sht' =F =F =F =F =F =F =F =F =F =F =F =F =F =F =F =F =F sht' sht' sht’ sht‘ 135 sht’ !$C$1 OFFSET 00000 RecCount=0 136 a _' 137 SIMULATION SCENARIO Initial Condition Continue current state to with Pw z=0 : cm suction: cm Redistribution Criteria Switch to redistribution at time: Outflow second cm3 0K Canoe! Reset Tune/omFlow 138 B. PRESSURE-SATURATION DATA OBTAINED USING SATURATION EQUILIBRIUM METHOD B. l. Air-Water Data 3.2. Air-Oil Data 139 PROJECT ME—109a SOIL TYPE Metea <850u SYSTEM Air 0“ PROC. NOTE Reff. for air-oil DATE 09-Jul-90 rho 1.622 @ 20 deg PREPARATION DATA 1 W by filter Paper SAMPLE WEIGHT OF Oil NO Air Dry Saturated absorbed 1 0.162 0.6832 0.5212 2 0.184 0.7808 0.5968 0.5g PREPARATION DATA 2 Initial Moisture Content SAMPLE Tin Air dry Oven dry lnit. NO + Tin + Tin Moisture 1 1.379 5.304 5.250 0.014 2 1.387 3.551 3.510 0.019 3 1.378 2.398 2.382 0.016 0.016399478 PREPARATION DATA 3 (in gram) CAPS: 3.994 ram CELL RING ALL SATURATED 8 AT THE LABEL CELL CELL,CAPS END OF PART BOTTOM FILLED EXPERIMENT <1) <2> <3> <4> 1 70.43 409.45 582.44 549.23 2 70.31 406.71 582.05 #NIA 3 70.59 407.70 582.95 549.71 4 71.96 407.59 583.03 550.16 5 70.54 406.52 583.13 550.95 M NOTE: Col<1>: annm+wmu+ Rimg+Upp¢part+ Lwrpan+0vn¢yplu+2ADFiham Col<3>z<2>+batomcqu+amuodooil+wbellawfillod FINAL PROC DATA For Residual Moisture Calculation CELL TIN TIN, BRAS TIN, BRAS LABEL RING, MO- RING, OVN IST SMPL DRY SMPL <5> <6> <7> 1 6.97 186.89 176.37 2 7.71 #N/A #NIA 3 6.97 187.03 176.45 4 6.82 188.61 177.84 5 6.9L 187.84 176.64 140 PRESSURE DATA Pressure ' 7.5 7.5 15 15 30 30 Date 12-Jul 13-Jul 13-Jul 14-Jul 14—Jul 16-Jul Time 01:14 PM 11:28 AM 11:29 AM 12:25 PM 02:47 PM 10:20 AM Temperature : 22 23 23 23 23 _23 Step: —START- --END—- -START- -—END— -START- —-END—- CELL <8a> <8a> <8b> <8c> <8c> 1 9.74 11.87 11.87 21.01 9.47 21.40 2 9.55 11.14 11.14 19.61 9.48 19.65 3 10.65 15.60 15.60 24.12 9.18 17.70 4 10.43 12.11 12.11 20.96 9.30 19.13 5 9.12 11.10 11.10 19.21 10.31 20.72 Pressure 40 40 50 50 75 75 Date 16-Jul 17-JuI 17-Jul 18-Jul 18-Jul 19-Jul Time 12:16 PM 10:20 AM 10:27 AM 11:50 AM 11:50 AM 09:16 AM Temperature : 23 23 23 23 23 23 Step: ~START— --END-— —START- -END-—- -START— -END— CELL <8d> <8d> <8a> <8e> <8f> <8f> 1 9.54 12.45 12.45 16.76 16.76 17.62 2 9.34 #N/A #N/A #N/A #N/A #N/A 3 9.66 13.16 13.16 15.82 15.82 17.65 4 9.07 13.88 13.88 17.14 17.14 18.09 5 10.50 16.31 16.31 18.49 18.49 19.14 Pressure 100 100 238 238 Date 19-Jul 20-Jul 20-Jul 21-Jul Time 10:39 AM 10:33 AM 10:39 AM 10:33 AM Temperature : 23 23 23 23 Step: -START- -END—- -START- —END—- CELL <8g> (89> <8h> <8b> 1 9.62 11.52 1 1.52 12.29 2 #NIA #NIA #NIA WA 3 9.65 12.17 12.17 12.41 4 9.09 10.67 10.67 11.97 5 10.14 11.85 11.85 12.82 DATA REDUCTION 1 Bulk Densil OVEN DRY BULK CELL SOIL WGT DENSITY <10> 1 1.453 2 WA 3 1.452 4 1.454 5 1.456 Note: Ring Vol: 68.119 cm3 141 DATA REDUCTION 2 ResidualMorsture Content _ RESIDUAL INITIAL WATER WATER MASS CELL 140151qu MOISTURE RELEASED COLLECTED BAIJANCE comm comm <11>=6-7 <12>=11+13 <13>=34 =sum(d8) «9:04.133: 1 10.52 43.73 33.21 33.95 1.69 2 #NIA #NIA #N/A #N/A #N/A 3 10.58 43.82 3324 32.74 -1.14 4 10.77 43.64 32.87 32.26 -1.40 5 11.20 43.38 32.18 31.82 -0.83 WEIGHT OF MOISTURE LOSS PRE: 7.5 15 30 40 50 75 100 238 CELL AT EACH STEP 1 2.13 9.14 11.93 2.91 4.31 0.86 1.90 0.77 2 1.59 8.47 10.17 #NIA #NIA WA #NIA #NIA 3 4.95 8.52 8.52 3.50 2.66 1.83 2.52 0.24 4 1.68 8.85 9.83 4.81 3.26 0.95 1.58 1.30 5 1.98 8.11 10.41 5.81 2.18 0.65 1.71 0.97 ACCUMULATIVE 1 2.13 9.14 11.93 2.91 4.31 0.86 1.90 0.77 2 3.72 17.61 22.10 #NIA #NIA #NIA #NIA WA 3 7.08 17.66 20.45 6.41 6.97 2.69 4.42 1.01 4 3.81 17.99 21.76 7.72 7.57 1.81 3.48 2.07 5 4.11 17.25 22.34 8.72 6.49 1.51 3.61 1.74 1 0.402 0.383 0.300 2 #N/A #NIA WA 3 0.392 0.347 0.270 4 0.389 0.374 0.294 5 0.389 0.371 0.298 AVG #NIA #NIA #NIA STD‘ NA NA #NIA DEGREE OF SATURATION 1 1.00 0.95 0.75 2 NA NA #NIA 3 1.00 0.89 0.89 4 1.00 0.98 0.78 5 1.00 0.95 0.77 0.193 #NIA 0.193 0.205 0.204 WA #NIA 0.188 NA 0.181 0.182 0.151 #NIA WA 0.41 #N/A 0.41 0.42 _ 0.39_ 0.127 #N/A 0.137 0.132 0.132 #N/A NA .32 #N/A 0.35 0.34 0.119 #NIA 0.121 0.124 0.128 NA NA 0.30 #NIA 0.31 0.32 0.25 #NIA 0.25 0.28 0.095 #NIA 0.098 0.097 0.101 #NIA #NIA 0.24 WA 0.24 0.25 142 PROJECT ME-109 SOIL TYPE Metea SYSTEM Air Water PROC. NOTE Replica #1 of 108 reff. for aw DATE 29-Sep-90 rho 1.021 @ 20 deg PREPARATION DATA 1 Wei ht of water absorbed b filter r SAMPLE WEIGHT OF Water NO Air Dry Saturated absorbed 1 0.171 0.427 0.256 2 0.194 0.488 0.294 SAMPLE Tin Moist Oven dry lnit. NO soiI+Tin + Tin Moisture 1 0.920 7.210 6.860 0.059 2 0.920 3.280 3.150 0.058 3 0.910 4.190 4.040 0.048 0.055047282 PREPARATION DATA 3 (in gram) CAPS: 8.201 Tram __ “ALL _“ATTRATE' "AT T~ ’ CELL. RING CELL CELL,CAPS END OF LABEL PART DTTOM FILLE EXPERIMENT <1> <2> <3> <4> 1 70.39 399.76 552.99 532.29 2 73.27 398.88 555.03 532.23 3 69.43 397.74 555.70 534.13 4 72.81 397.08 553.85 533.78 5 73.50 393.05 552.38 531 .33 M NOTEICOI <1>: Wingnuts+Washa+ijg+Uppaput+Lwrpu1+0vn¢yplatc+2ADFflwpaper Col<3>:<2>+bottomcaps+aaunatcdsoi1+apaocbeflow FINAL PROC DATA For Residual Moisture Calculation w CELL . TIN TIN, BRAS TIN, BRAS LABEL RING, MO- RING, OVN IST SMPL DRY SMPL <5> <6> <7> 1 3.83 177.51 172.51 2 3.83 180.44 174.43 3 3.83 176.40 171.74 4 3.81 180.52 175.99 5 3.81 181.36 176.19 143 DATA REDUCTION 2 _R_e§°dual Moisture Content _ _ _ . RESIDUAL INITIAL ¥ WATER 7 WATER MASS CELL MOISTURE MOISTURE RELEASED COLLECTED BALLANCE CONTENT CONTENT <11>=6-7 =l1+l3 <13>=34 =sum(d8) <15>=(I4—13)/12 1 5.00 25.70 20.70 21.13 1.88 2 8.01 28.81 22.80 22.53 -0.94 3 4.88 28.23 21.57 21.12 -1.72 4 4.53 24.80 20.07 20.19 0.49 5 5.17 28.22 21.05 21.70 2.48 WEIGHT OF MOISTURE LOSS PRESS 17 . 49 88 89 117 187 317 717 CELL AT EACH STEP 1 0.73 5.80 8.97 1.97 1.42 2.07 2.04 0.33 2 0.92 1.72 8.80 5.03 2.74 1.18 2.08 0.08 3 0.80 8.40 8.40 1 .53 2.89 0.93 1 .80 0.77 4 0.98 5.99 8.84 1.43 2.01 1.13 1.74 0.29 5 1.09 1.02 9.82 4.03 1.72 0.95 2.58 0.71 ACCUMULATIVE 1 0.73 8.33 13.30 15.27 18.89 18.78 20.80 21.13 2 0.92 2.84 11.44 18.47 19.21 20.39 22.45 22.53 3 0.80 7.00 13.40 14.93 17.82 18.55 20.35 21.12 4 0.98 8.95 13.59 15.02 17.03 18.18 19.90 20.19 5 1.09 2.11 11.73 15.78 17.48 18.43 20.99 21.70 VOLUMETRIC MOISTURE CONTENT PRESS 0 17 49 88 89 117 187 317 717 CELL 1 0.384 0.373 0.291 0.188 0.159 0.139 0.108 0.078 0.073 2 0.419 0.408 0.380 0.251 0.177 0.137 0.120 0.089 0.088 3 0.378 0.370 0.278 0.182 0.159 0.120 0.108 0.080 0.088 4 0.383 0.349 0.281 0.183 0.142 0.113 0.098 0.071 0.087 5 0.394 0.378 0.383 0.222 0.183 0.138 0.124 0.088 0.078 AVG 0.388 0.375 0.314 0.201 0.180 0.129 0.111 0.081 0.075 STDV 0.019 0.018 0.048 0.031 0.011 0.011 0.010 0.007 0.008 DEGREE OF SATURATION 1 1.00 0.97 0.78 0.49 0.42 0.38 0.28 0.20 0.19 2 1.00 0.97 0.91 0.80 0.42 0.33 0.29 0.21 0.21 3 1.00 0.98 0.73 0.48 0.42 0.32 0.28 0.21 0.18 4 1.00 0.98 0.72 0.45 0.39 0.31 0.27 0.20 0.18 5 1.00 0.98 0.92 0.58 0.41 0.35 0.31 0.22 0.19 144 PRESSURE DATA Pressure : 1 7 49 66 Date : 01-Aug 02-Aug 02-Aug 03-Aug 03-Aug 05-Aug Time : 01:40 PM 11:28AM 11:28AM 11:17AM 11:17AM 01:30AM Temperature : 22 23 23 23 23 23 Step: «START— --END-— -START- --END-- -START-— --END-—- CELL 10.40 11.13 11.13 18.73 9.35 18.32 9.37 10.29 10.29 12.01 9.88 18.68 9.65 10.25 10.25 16.85 9.45 15.85 10.01 10.97 10.97 18.98 10.12 18.76 9.32 10.41 10.41 11.43 9.18 18.80 I IUI&00Nd T>ressure : 89 117'? 167 Date : 05-Aug 06-Aug 06-Aug 07-Aug 07-Aug 08-Aug Time : 01:30 AM 11:34 AM 11:34 AM 01:57 PM 01 :57 PM 01:57 PM Temperature : 23 23 23 23 23 23 Step: -START- --END- -START- --END— -START—- --END— CELL 1 9.65 11.62 11.62 13.04 13.04 15.11 2 9.12 14.15 14.15 16.89 16.89 18.07 3 9.56 11.09 11.09 13.78 13.78 14.71 4 9.88 11.31 11.31 13.32 13.32 14.45 5 9.21 13.24 13.24 14.96 14.96 15.91 _Pressure : 317 717 Date : 08-Aug 09-Aug 09-Aug 10—Aug Time : 01:57 PM 12:50 PM 12:50 PM 12:50 PM Temperature : 23 23 23 23 Step: -START- --END-—- -START- --END—- CELL 1 15.11 17.15 9.43 9.76 2 18.07 20.13 9.45 9.53 3 14.71 16.51 10.05 10.82 4 14.45 16.19 10.11 10.40 5 15.91 18.47 9.32 10.03 DATA REDUCTION 1 Bulk Densit OVEN DR‘ BULK CELL SOIL WG' DENSITY <9> <10> 1 98.29 1.443 2 97.33 1.429 3 98.48 1.446 4 99.37 1.459 5 98.88 1.452 m Note: Ring Vol: 68.119 145 C. REGRESSION OUTPUT C. l. Air-Water System, Based on Data Measured Using Saturation Equilibrium C.2. 145133;.“ System, Based on Data Measured Using Pressure Equilibrium C.3. lAiiret-IOIIOdSystem, Based on Data Measured Using Saturation Equilibrium C.4. $31an System, Based on Data Measured Using Pressure Equilibrium 146 SPSS Result Air-Hater System A” Saturation Equilibrium. Ami. Cell 1 Nonlinear Regression Summary Statistics Dependent Variable SE_Sl Source DF Sum of Squares Mean Square Regression 3 3.21827 1.07276 Residual 6 2.095629E-03 3.492714E-04 Uncorrected Total 9 3.22037 (Corrected Total) 8 .79303 R squared = 1 - Residual SS / Corrected SS = .99736 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .185867309 .016738882 .144908741 .226825877 A .023725148 .001067714 .021112546 .026337751 N 2.608953612 .153587725 2.233137989 2.984769236 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 -.1204 .7100 A -.1204 1.0000 -.6550 N .7100 -.6550 1.0000 1.2. Cell 2 Nonlinear Regression Summary Statistics Dependent Variable SE_SZ Source DF Sum of Squares Mean Square Regression 3 3.59751 1.19917 Residual 6 1.829279E-03 3.048798E-04 Uncorrected Total 9 3.59934 (Corrected Total) 8 .85753 R squared = 1 - Residual SS / Corrected SS = .99787 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .223956252 .012046083 .194480548 .253431956 A. .017169866 .000465797 .016030102 .018309629 N 3.759266303 .232879374 3.189431003 4.329101604 .Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 .0872 .5541 .A .0872 1.0000 -.5472 N .5541 -.5472 1.0000 Am3. Cell 3 Nonlinear Regression Summary Statistics Dependent variable SE_53 Source DF Sum of Squares Mean Square Regression 3 3.13448 1.04483 Residual 6 2.192198E-03 3.653664E-04 Uncorrected Total 9 3.13667 (Corrected Total) 8 .80941 R squared = 1 - Residual SS / Corrected SS = .99729 147 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .179457012 .016845170 .138238366 .220675657 A .024686491 .001141628 .021893028 .027479954 N 2.611190050 .156043455 2.229365469 2.993014630 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 -.1260 .7096 A -.1260 1.0000 -.6586 N .7096 -.6586 1.0000 A.4. Cell 4 Nonlinear Regression Summary Statistics Dependent variable SE_S4 Source DF Sum of Squares Mean Square Regression 3 3.01430 1.00477 Residual 6 1.613228E-O3 2.688714E-04 Uncorrected Total 9 3.01592 (Corrected Total) 8 .82072 R squared = 1 - Residual SS / Corrected SS = .99803 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .175272285 .013755495 .141613801 .208930769 A .025508022 .000986397 .023094396 .027921649 N 2.681842856 .137600831 2.345145753 3.018539960 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 -.1170 .6982 A. -.1170 1.0000 -.6552 N .6982 -.6552 1.0000 .A.5. Cell 5 Nonlinear Regression Summary Statistics Dependent variable SE_SS Source DF Sum of Squares Mean Square Regression 3 3.53873 1.17958 Residual 6 7.839523E-03 1.306587E-03 Uncorrected Total 9 3.54657 (Corrected Total) 8 .86994 R squared = 1 - Residual 38 / Corrected SS = .99099 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .219392728 .025020222 .158170451 .280615005 A. .017738378 .001020385 .015241585 .020235171 N 3.674507015 .461516478 2.545216874 4.803797155 Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 .0830 .5612 .A .0830 1.0000 -.5479 N .5612 -.5479 1.0000 148 ASS. 83 all Nonlinear Regression Summary Statistics Dependent Variable SE_S Source DF Sum of Squares Mean Square Regression 3 16.45007 5.48336 Residual 42 .06880 1.637981E-03 Uncorrected Total 45 16.51887 (Corrected Total) 44 4.17228 R squared = 1 - Residual SS / Corrected SS = .98351 Asymptotic 95 5 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .193063491 .014628308 .163542370 .222584612 A .021273045 .000791200 .019676340 .022869751 N 2.907024543 .168339440 2.567301800 3.246747286 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 -.0372 .6603 A. -.0372 1.0000 -.6124 N .6603 -.6124 1.0000 8. Preaure‘lquilibriumi 8.1. Cell 1 Nonlinear Regression Summary Statistics Dependent Variable PE_Sl Source DF Sum of Squares Mean Square Regression 3 8.72360 2.90787 Residual 15 .01589 1.059621E-03 Uncorrected Total 18 8.73949 (Corrected Total) 17 1.17524 R squared = 1 - Residual SS / Corrected 55 = .98648 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .265765467 .022804424 .217158987 .314371947 A. .021889159 .000655873 .020491200 .023287119 N 3.134065543 .258524024 2.583034630 3.685096456 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 .3154 .7042 .A .3154 1.0000 -.1477 N .7042 -.1477 1.0000 3.2. Cell 2 Nonlinear Regression Summary Statistics Dependent Variable PE_82 Source DF Sum of Squares Mean Square Regression 3 8.18419 2.72806 Residual 15 .01149 7.657700E-04 Uncorrected Total 18 8.19568 (Corrected Total) 17 1.15263 R squared = 1 - Residual SS / Corrected SS = .99003 149 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .265465274 .016785044 .229688800 .301241747 A .024563224 .000477704 .023545022 .025581426 N 3.943358281 .271261578 3.365177913 4.521538649 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 .3820 .6094 A .3820 1.0000 -.1490 N .6094 -.1490 1.0000 8.3. Cell 3 Nonlinear Regression Summary Statistics Dependent Variable PE_S3 Source DE Sum of Squares Mean Square Regression 3 8.46954 2.82318 Residual 15 .01051 7.008422E-04 Uncorrected Total 18 8.48006 (Corrected Total) 17 .94411 R squared = 1 - Residual SS / Corrected SS = .98887 Asymptotic 95 % Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .236929406 .025652718 .182251932 .291606881 A .023017836 .000700355 .021525065 .024510607 N 2.724943257 .165280030 2.372657212 3.077229303 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 .3217 .7502 A. .3217 1.0000 -.2121 N ' .7502 -.2121 1.0000 8.4. Cell 4 Nonlinear Regression Summary Statistics Dependent variable PE_S4 Source DF Sum of Squares Mean Square Regression 3 8.76843 2.92281 Residual 15 .02138 1.425210E-03 Uncorrected Total 18 8.78981 (Corrected Total) 17 .98161 R squared = 1 - Residual 55 / Corrected 33 = .97822 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .224802526 .025355605 .170758333 .278846720 A. .020010398 .000487883 .018970500 .021050296 N 3.908730532 .361630200 3.137934007 4.679527058 Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 .4754 .4657 A .4754 1.0000 -.2153 N .4657 -.2153 1.0000 150 8.5. Cell 5 Nonlinear Regression Summary Statistics Dependent Variable PE_SS Source DE Sum of Squares Mean Square Regression 3 8.37428 2.79143 Residual 15 .01095 7.299417E-04 Uncorrected Total 18 8.38523 (Corrected Total) 17 1.00163 R squared = 1 - Residual SS / Corrected SS = .98907 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .213951001 .023159461 .164587779 .263314223 A .019195945 .000491686 .018147941 .020243949 N 2.996441761 .184490669 2.603209209 3.389674313 ’ Asymptotic Correlation Matrix of the Parameter Estimates SR A N i SR 1.0000 .3411 .6886 I A .3411 1.0000 -.2401 L N .6888 -.2401 1.0000 8.6. Cell 6 Nonlinear Regression Summary Statistics Dependent variable PE_S6 Source DF Sum of Squares Mean Square Regression 3 8.31781 2.77260 Residual 15 .03086 2.0576088-03 Uncorrected Total 18 8.34867 (Corrected Total) 17 1.30710 R squared = 1 - Residual 85 / Corrected $5 = .97639 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .237359602 .024306009 .185552570 .289166635 .A .023323244 .000740512 .021744879 .024901608 N 3.721848523 .483569795 2.691143904 4.752553142 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 .2700 .5130 A. .2700 1.0000 -.3104 N .5130 -.3104 1.0000 8.7. Cell 7 Nonlinear Regression Summary Statistics Dependent Variable PE_S7 Source DF Sum of Squares Mean Square Regression 3 7.42435 2.47478 Residual 15 .05292 3.528294E-03 Uncorrected Total 18 7.47728 (Corrected Total) 17 1.12329 R Squared = 1 - Residual 58 / Corrected SS = .95288 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .221494731 .039091195 .138173821 .304815641 A- .020132438 .000656554 .018733025 .021531851 N 4.317164518 .764616195 2.687423678 5.946905358 151 .Asymptotic Correlation Matrix of the Parameter Estimates SR A SR 1.0000 .3580 A .3580 1.0000 N .6546 -.l672 8.8. P! All N .6546 -.1672 1.0000 Nonlinear Regression Summary Statistics DE Sum of Squares Source Regression 3 Residual 123 Uncorrected Total 126 (Corrected Total) 125 R squared = 1 - Residual SS / Corrected SS = Dependent Variable PE_S Mean Square 58.08456 19.36152 .33164 2.696268E-03 58.41620 7.73521 Asymptotic Std. Error .014684055 .000383119 .147475130 2.828453167 .95713 .Asymptotic 95 8 Confidence Interval Lower Upper .200429752 .258562125 .021088691 .022605412 3.412289123 .Asymptotic Correlation Matrix of the Parameter Estimates Parameter Estimate SR .229495938 A .021847051 N 3.120371145 SR SR 1.0000 A .2799 N .6713 C. Pooled SE and 88 data A. .2799 1.0000 -.2870 N .6713 -.2870 1.0000 Nonlinear Regression Summary Statistics Source Regression Residual Uncorrected Total (Corrected Total) R squared = l - Residual 83 / Corrected SS = DE Sum of Squares 3 168 171 170 Dependent variable S Mean Square 74.52789 24.84263 .40717 2.423657E-03 74.93507 12.31172 .96693 .Asymptotic 95 8 Confidence Interval Lower Upper 192365905 .235430670 021027223 .022352698 .110100837 2.803656840 3.238375722 Asymptotic Parameter Estimate Std. Error SR .213898288 .010906972 A. .021689960 .000335702 N 3.021016281 .Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 .1848 .6725 .A .1848 1.0000 -.3740 N .6725 -.3740 1.0000 152 8888 Result Air-Oil System A. Saturation Equilibrium A.1. Cell 1 Nonlinear Regression Summary Statistics Dependent Variable 551 Source DF Sum of Squares Mean Square Regression 3 3.14637 1.04879 Residual 6 1.654965E-03 2.758274E-04 Uncorrected Total 9 3.14803 (Corrected Total) 8 .73555 R squared = 1 - Residual 35 / Corrected $3 = .99775 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .216736677 .015467300 .178889556 .254583797 A. .062351208 .002559450 .056088460 .068613957 N 2.618260450 .140317205 2.274916619 2.961604281 .Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 -.O595 .7526 A. -.0595 1.0000 -.5531 N .7526 -.5531 1.0000 A12. Cell 2 Nonlinear Regression Summary Statistics Dependent Variable 552 Source DF Sum of Squares Mean Square Regression 3 3.04630 1.01543 Residual 6 1.396263E-03 2.327105E-04 Uncorrected Total 9 3.04769 (Corrected Total) 8 .61374 R squared = 1 - Residual SS / Corrected 53 = .99772 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .208232505 .019519642 .160469662 .255995348 .A .079441789 .004012099 .069624537 .089259041 N 2.134933754 .105663076 1.876385521 2.393481986 .Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 -.2380 .8413 A. -.2380 1.0000 -.6373 N .8413 -.6373 1.0000 Am3. Cell 3 Nonlinear Regression Summary Statistics Dependent variable SS3 Source DF Sum of Squares Mean Square Regression 3 3.29934 1.09978 Residual 6 2.0448078-03 3.4080128-04 Uncorrected Total 9 3.30139 (Corrected Total) 8 .69184 R squared = 1 - Residual 55 / Corrected SS = .99704 Asymptotic 95 8 Asymptotic Confidence Interval 153 Parameter Estimate Std. Error Lower Upper SR .237269389 .017900322 .193468879 .281069900 A .060570293 .002899210 .053476182 .067664404 N 2.574162336 .157311212 2.189235666 2.959089006 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 -.0662 .7579 A -.0662 1.0000 -.5611 N .7579 -.5611 1.0000 A.4. Cell 4 Nonlinear Regression Summary Statistics Dependent Variable SS4 Source DF Sum of Squares Mean Square Regression 3 3.31058 1.10353 Residual 6 2.165306E-03 3.608844E-04 Uncorrected Total 9 3.31275 (Corrected Total) 8 .67805 R squared = 1 - Residual SS / Corrected $8 = .99681 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .256779702 .016891314 .215448146 .298111257 A .060096911 .002837410 .053154018 .067039803 N 2.721194994 .176107274 2.290276017 3.152113971 Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 -.0316 .7354 .A -.0316 1.0000 -.5389 N .7354 -.5389 1.0000 A15. SI ALL Nonlinear Regression Summary Statistics Dependent variable SS_ALL Source DF Sum of Squares Mean Square Regression _ 3 12.79225 4.26408 Residual 33 .01760 5.334762E-04 Uncorrected Total 36 12.80986 (Corrected Total) 35 2.72315 R squared = 1 - Residual SS / Corrected SS = .99354 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .230852388 .011515977 .207422957 .254281820 A. .064636301 .001999256 .060568785 .068703818 N 2.495353457 .094577587 2.302933909 2.687773005 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 -.0941 .7733 A -.0941 1.0000 -.5712 N .7733 -.5712 1.0000 154 8. Presure'lquilibriumi 8.1. Cell 1 Nonlinear Regression Summary Statistics Dependent Variable PSl Source DF Sum of Squares Mean Square Regression 3 6.80869 2.26956 Residual 10 .02558 2.557503E-03 Uncorrected Total 13 6.83426 (Corrected Total) 12 .78736 R squared = 1 - Residual 55 / Corrected $5 = .96752 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .216224062 .064561620 .072371806 .360076315 A .067541880 .005856179 .054493500 .080590260 N 2.312055769 .301804838 1.639592684 2.984518853 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 .2226 .7835 A .2226 1.0000 -.2716 N .7835 -.2716 1.0000 8.2. Cell 2 Nonlinear Regression Summary Statistics Dependent variable P52 Source DF Sum of Squares Mean Square Regression 3 6.40744 2.13581 Residual 10 .01025 1.024532E-03 Uncorrected Total 13 6.41768 (Corrected Total) 12 .97341 R squared = 1 - Residual SS / Corrected SS = .98947 .Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .222122383 .022709674 .171522076 .272722689 A. .052948823 .001498861 .049609152 .056288493 N 4.068471978 .344772143 3.300271770 4.836672186 Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 .3767 .4367 A. .3767 1.0000 -.3795 N .4367 -.3795 1.0000 8.3. Cell 3 Nonlinear Regression Summary Statistics Dependent variable P33 Source DE Sum of Squares Mean Square Regression 3 6.54629 2.18276 Residual 10 5.794601E-03 5.794601E-04 Uncorrected Total 13 6.55409 (Corrected Total) 12 .85967 R squared = 1 - Residual 55 / Corrected 58 = .99326 155 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .240359759 .021526053 .192396725 .288322793 A .068707315 .001834661 .064619436 .072795193 N 3.177073165 .207895561 2.713852990 3.640293341 Asymptotic Correlation Matrix of the Parameter Estimates SR A N SR 1.0000 .3517 .6682 A .3517 1.0000 -.1998 N .6682 -.1998 1.0000 8.4. Cell 4 Nonlinear Regression Summary Statistics Dependent variable PS4 Source DE Sum of Squares Mean Square Regression 3 6.27705 2.09235 Residual 10 2.939531E-03 2.939531E-04 Uncorrected Total 13 6.27999 (Corrected Total) 12 .82741 R squared = 1 - Residual SS / Corrected 88 = .99645 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .270039190 .018093034 .229725398 .310352982 .A .058346729 .001846515 .054232436 .062461022 N 2.464993945 .122208989 2.192695348 2.737292541 Asymptotic Correlation Matrix of the Parameter Estimates SR .A N SR 1.0000 .0668 .7907 A. .0668 1.0000 -.4043 N .7907 -.4043 1.0000 8.5. P! All Nonlinear Regression Summary Statistics Dependent variable PséALL Source DF Sum of Squares Mean Square Regression 3 25.93174 8.64391 Residual 49 .15428 3.1486143-03 Uncorrected Total 52 26.08602 (Corrected Total) 51 3.45833 R squared = 1 - Residual SS / Corrected SS = .95539 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .256510471 .024897608 .206476855 .306544088 .A .063704804 .002514772 .058651180 .068758429 N 2.797045959 .206964647 2.381134929 3.212956988 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 .2279 .6752 .A .2279 1.0000 -.3471 N .6752 -.3471 1.0000 156 C. Pooled 83 and PB date Nonlinear Regression Summary Statistics Dependent Variable S_TOTAL Source DF Sum of Squares Mean Square Regression 3 38.72051 12.90684 Residual 85 .17537 2.063226E-03 Uncorrected Total 88 38.89588 (Corrected Total) 87 6.54283 R squared = 1 - Residual 55 / Corrected $5 = .97320 Asymptotic 95 8 Asymptotic Confidence Interval Parameter Estimate Std. Error Lower Upper SR .245726009 .014516142 .216864030 .274587988 A .064034385 .001806441 .060442697 .067626074 N 2.659540349 .123329373 2.414328515 2.904752183 Asymptotic Correlation Matrix of the Parameter Estimates SR A. N SR 1.0000 .0716 .7271 A .0716 1.0000 -.4429 N .7271 -.4429 1.0000 157 D. INTERFACIAL TENSION DATA D.1. D.2. D.3. D4. D5. D.6. Air-Water System - Metea Soil Oil-wter System - Metea Soil Air-Oil System - Metea Soil Air-Water System - Ottawa Sand Oil-Water System - Ottawa Sand Air-Oil System - Ottawa Sand 158 Air. (-) Water: air-water experiment (ME 1 14 and ME 1 15) me114#1la me 114#2/a Density= 0.974 0.9836 Mag Coef= 49.325 50 de 05 y de 08 7 109.3 59.8 71.85 111.6 61.8 71.34 109.1 60.2 70.02 111.4 61.8 70.75 109.1 60 70.63 1 1 1.5 62 70.45 110 60.6 71.48 111.1 61.7 70.17 109.9 60.7 70.87 110.9 62 68.72 109.9 60.5 71.49 110.9 61.6 69.89 110 61.5 68.79 111.3 61.9 70.16 110.2 61.3 69.95 111.1 62.1 69.00 110.2 61.6 69.07 111.1 61.9 69.59 110 61 70.26 110.7 61.2 70.50 110.1 61 70.56 110.7 61.1 70.80 110 61 .1 69.96 110.9 61.3 70.78 109.8 61 69.68 1 1 1 61.8 69.59 109.8 60.9 69.98 111 62 69.01 109.8 60.8 70.28 1 11 61.8 69.59 109.9 60.9 70.27 111 61.8 69.59 109.7 60.7 70.28 110.9 61.8 69.30 109.9 60.5 71.49 111 61.8 69.59 AVG 70.38 69.93 STDEV 0.802 0.707 me 114113/a me 114 #4Ia Density= 0.9818 0.9766 Mag Coeff- 50.2 49.975 ’ de 05 y de ds 1 110 61 68.38 111.1 62.5 67.44 109.9 61.4 66.94 111.1 62.2 68.29 109.9 61.2 67.52 111 62.1 68.30 110.2 61.4 67.79 111.1 62.6 67.16 110.1 61.4 67.51 111.1 62.7 66.88 110 61.1 68.09 110.9 62.6 66.61 110 61.4 67.22 110.9 61.9 68.59 109.5 61 66.96 110.8 61.7 68.88 110 61.3 67.51 111 61.8 69.16 110.2 61 68.95 111 62 68.58 110.2 61.5 67.50 110.9 61.8 68.88 110.1 61.5 67.22 110.9 61.9 68.59 109.8 61.1 67.52 110.3 62.5 65.23 109.8 61 67.81 110.3 62.4 65.51 109.5 61.4 65.83 110.4 62.4 65.78 110 61.8 66.10 110.8 62.2 67.45 109.9 61.8 65.82 110.9 62.2 67.73 110 61.8 66.10 110.9 62.3 67.45 AVG 67.26 67.58 STDEV 0.841 1.171 159 Air. (-) Water. air-water experiment (ME 114 and ME 115) me115#1la me115#2la Density= 0.9576 0.9738 Mag C0ef= 50 49.875 de ds 1 de ds 7 110.6 60.9 69.23 110 61 68.71 110.5 61.1 68.35 110 60.9 69.00 110.5 61.1 68. 35 110 61 68.71 110.9 61 69. 80 111.1 61.8 69.53 110.9 61.1 69.50 111.2 61.6 70.41 110.9 61.2 69.21 111.1 61.6 70.12 110.6 61 68.93 110.7 61.2 70.14 110.6 61.2 68.35 110.6 61.4 69.26 110.7 61.1 68.93 110.7 61.2 70.14 110.9 61.4 68.62 111 62 68. 66 111 61.6 68.32 110.9 62.2 67. 81 111 61.7 68.04 111 62 68.66 110.6 61.1 68.64 110.9 61.9 68.67 110.7 61 69. 22 110.9 61.8 68.96 110.6 61.2 68. 35 110.9 62 68. 38 110.7 61.1 68. 93 110.6 61.7 68.39 1 10.8 61.1 69.21 1 10.4 62 66.97 110.7 61.3 68.34 110.6 62 67.53 AVG 68.80 68.89 STDEV 0.478 0.912 AVG 68.65 me 11583/a me 115#4/a STDEV 1.51 Density= 0.9721 0.9804 N 144.00 Mag Coef= 46.66 50.325 t-95 3.411228 09 ds 7 06 ds 7 103.1 57.7 67.22 110.5 61.9 66.78 103.1 57.8 66. 92 110.6 62 66.78 103.1 58 66. 32 110.5 61.9 66.78 103.2 57.8 67. 22 112.1 62. 5 69.58 103.2 57.7 67.52 112.2 62. 5 69.87 103.3 57.8 67.52 112.2 62. 5 69.87 103.1 57.6 67.52 112.1 62. 5 69.58 103.1 57.8 66.92 112.3 62 .6 69.86 103.2 57.6 67.83 1 12.1 62 71.05 103.1 57.9 66.62 112.1 62. 7 69. 00 103.1 57.9 66.62 112.1 62.9 68.43 103.1 58 66.32 112.2 62.8 69. 00 103 57.9 66.32 1 12.4 62.5 70. 44 103 57.9 66.32 112.3 62.5 70.15 103 57.7 66.92 1 12.3 62.7 69.57 103.1 57.8 66.92 111.9 62.1 70.17 103 57.8 66.62 1 12 61.6 71.96 103 57.8 66. 62 111.9 61.8 71.06 AVG 66.90 69.44 STDEV 0.465 1.431 OiI: Air-oil experiment (ME 124) Water oil-water experiment (ME 125) Density= dag Coef= AVG STDEV Density= Jag Coef= #1 0.6028 43.5 de 93.5 93.5 93.5 91.5 91.5 91.6 96.9 96.9 96.9 #3 0.6088 42.817 de 93.5 93.5 93.5 90.5 90.7 90.5 92.7 92.6 92.5 AVG 53 53.1 53 52.5 52.5 52.2 54.5 55 54.7 52.1 52.3 52.4 51 .4 51 .2 51 .2 52.9 52.8 52.8 STDEV 38.15 37.98 38.15 35.40 35.40 36.11 41.82 40.84 41.42 38.38 2.365 41.45 41.04 40.83 38.95 37.71 37.33 38.29 38.67 38.10 38.93 1.819 160 0.5988 41 .75 de 91 91 .1 91 89.2 89.1 89.3 93.5 93.7 93.5 0.5918 42.95 93.4 93.2 93.3 94.1 93.1 93 93 95.1 95.1 95.1 95.7 95.9 95.8 96.1 96.2 98.2 52 51.7 51.9 51.2 51.5 51.2 52.6 52.2 52.9 ds AVG 54.1 53.9 52.8 52.7 52.8 54.3 54.1 54.5 54.5 54.5 53.9 55.1 54.8 55.1 STDEV AVG STDEV {-95 38.16 38.94 38.35 38.24 35.51 36.43 41.98 43.23 41.35 38.91 2.575 7 38.22 38.04 38.39 37.49 37 .49 37.87 38.03 38.03 37.84 38.99 39.38 38.81 39.75 40.14 41 .1 1 39.38 40.14 39.57 38.48 1 .454 38.82 1 .97815 45 4.590717 24.01 24. 1 8 24.85 25.25 25.25 25.25 25.25 25.73 25.59 25.02 0.571 25.53 24.85 25.53 25.17 25.51 25.15 25.56 24.84 25.03 25.22 Air. (-) Oil: air-oil experiment (ME 118) me 118 0a It 1 Density= 1.598 dag Coef= 46.824 de ds 60.1 39.8 60.1 39.7 60.1 39.4 60.5 39.5 60.5 39.5 60.5 39.5 60.5 39.5 60.4 39.1 60.5 39.3 AVG STDEV me 118 0a # 3 Density= 1.5746 wag Coef= 46.396 de ds 60.6 39.5 60.4 39.8 60.6 39.5 60.5 39.6 60.5 39.4 60.4 39.5 60.7 39.6 60.5 39.8 60.6 39.8 AVG STDEV 0.317 161 me 118 0a # 2 1.5982 46.432 de ds 60.9 39.9 60.9 39.9 60.9 39.7 60.7 40.1 60.9 39.7 60.9 39.9 60.9 40.0 60.9 40.0 60.9 39.9 me 116 0a # 4 1.5867 46.204 de 05 60.9 40.8 61.1 40.7 61.1 40.8 61.0 40.1 61.1 39.9 61.9 402 61.0 41.0 61.0 41.0 61.1 40.9 AVG STDEV N 1-95 25.78 25.78 28.13 25.06 28.13 25.78 25.82 25.82 25.78 25.74 0.297 24.38 24.92 24.78 25.71 28.25 27.33 24.28 24.28 24.60 25.18 0.999 25.29 0.89 27.00 1 .850298 Air : Water. Oil-Water experiment (OS 121 and OS 122) as 121 ow #1Ia Density= dag Coef= AVG STDEV Density= dag Coef= 46.275 AVG STDEV Density= dag Coef= AVG STDEV (-) 0.9844 45.975 de ds 102.8 56.1 102.6 56.6 102.7 56.4 102.6 56.9 102.6 56.5 102.6 56.9 103.0 57.0 102.9 57.2 102.9 56.9 as 121 ow #318 0.9769 de ds 103.3 57.3 103.5 56.8 103.5 57.0 103.7 57.3 103.5 57.5 103.6 57.1 103.4 57.1 103.3 57.7 103.3 57.1 as 122 CW #1Ia 0.9824 46.2 de 05 103.0 57.0 103.0 57.1 103.1 56.9 101.6 57.2 101.7 56.9 101.9 56.9 103.0 57.5 102.9 57.5 103.0 57.4 74.43 72.09 73.08 71 .10 72.42 71 .10 72.08 71 .08 72.06 72.18 1.028 70.58 72.84 72.17 71.83 70.55 72.18 71.53 89.30 71.21 71.35 1.017 71 .21 70.89 71 .88 88.25 87.47 88.09 89.81 89.30 89.93 89.40 1 .738 162 05 121 ow #2/a 0.9773 46 de ds 102.8 56.9 102.9 56.7 102.9 56.8 103.0 56.6 102.9 56.6 103.1 56.1 102.9 57.0 102.9 56.9 103.0 57.0 as 121 ow #4Ia 0.9824 46.325 de ds 103.8 57.2 103.9 57.1 103.9 56.9 103.4 57.5 103.5 57.5 103.2 57.4 103.8 57.9 103.9 57.9 104.0 57.6 as 122 0w # 21a 0.9832 46.25 de ds 103.1 56.7 103.1 57.0 103.1 56.9 103.1 56.7 103.0 56.6 103.0 57.0 103.0 57.0 103.1 56.8 103.1 56.9 71.15 72.12 71.79 72.78 72.45 74.81 71.14 71.46 71.46 72.13 1.089 72.73 73.38 74.06 70.48 70.79 70.17 70.48 70.77 72.06 71 .86 1 .358 72.42 71.43 71.78 72.42 72.43 71.1 1 71.11 72.09 71.78 71.84 0.507 Air: (-) 163 Water. Oil-Water experiment (OS 121 and OS 122) Density= Mag Coef= AVG STDEV es 122 ow # 3/a 0.979 47.175 de ds 106.0 58.5 106.0 58.8 106.0 58.7 106.0 59.0 106.0 58.9 105.9 58.9 106.0 59.0 106.0 58.7 105.9 59.0 7 72.60 71.64 71.96 71.01 71.32 71.01 71.01 71.96 70.70 71.47 0.580 as 122 0w #4/a 0.9816 47.675 de ds 105.9 58.0 105.9 57.9 105.8 58.0 105.9 57.5 105.7 57.9 105.9 57.8 105.9 58.0 105.9 57.8 105.9 . 58.1 AVG STDEV N t-95 7 72.57 72.89 72.25 74.22 72.28 73.22 72.57 73.22 72.24 72.83 0.81 3 72.14 1 .432 72 1 .881 245 164 Oil : Air-oil experiment (OS 123) Water: oil-water experiment (OS 12) -os1220w#1 os1220w#2 Density= 0.61 79 0.6048 .439 Coef= 40.825 44.1 de ds 1 de 05 y 83. 1 46.0 37.31 87.9 46.8 38.74 83.0 46.0 37.10 87.8 46.8 38.54 83-1 45.9 37.52 87.8 47.0 38.12 84.2 46.2 39.19 87.2 47.0 36.93 84.5 46.2 39.84 87.1 46.8 37.15 84.2 46.0 39.64 87.1 47.0 36.74 82.8 46.0 36.69 88.8 48.2 37.61 82.8 46.2 36.28 88.8 48.0 38.02 82.7 45.9 36.70 88.8 48.2 37.61 AVG 37.81 37.72 STDEV 1.292 0.658 es 122 ow #3 as 122 ow #4 Density= 0.6133 0.5827 hag Coef= 40.1 5 43.2 de 05 y de ds 7 88.1 50.9 38.49 91.8 50.1 39.79 88.1 51.6 37.14 91.9 50.1 39.99 88.0 51.0 38.09 92.0 50.1 40.19 86.1 50.1 36.08 89.6 48.9 37.90 86.1 50.0 36.27 89.3 48.7 37.72 86.0 50.0 36.08 89.4 49.0 37.32 87.9 51.0 37.89 91.9 50.1 39.99 87 .9 51.0 37.89 91.9 50.0 40.20 88.0 50.9 38.29 91.9 50.0 40.20 AVG 37.36 39.25 STDEV 0.928 1.152 AVG 38.03 STDEV 1.281291 N 36 t-95 3.00074 25.49 25.68 25.85 25.54 28.06 25.54 25.51 26.19 26.38 25.80 0.314 28.83 27.84 27. 88 26.09 28. 82 28.28 26.97 28. 80 28.97 27.14 Air : (-) Oil : air-oil experiment (OS 123) as 123 so #1 la Density= 1.6013 JIag Coef= 47.775 de ds 62.0 40. 5 62.1 40. 5 62.1 40. 4 622 40.7 62.3 40.5 62.2 40.7 62.1 40.6 62.1 402 62.2 40.2 AVG STDEV es 123 so # 3Ia Density= 1.5913 illag Coef= 46.3 de ds 62.9 40.5 62.6 40. 7 62.7 40. 8 61.9 40. 9 62.0 40. 7 62.0 40. 9 62.0 40 .5 62.0 40. 6 62.0 40. 5 AVG STDEV 0.825 165 05 123 a0 #2/a 1.5851 47.275 de ds 62.6 40 9 62.7 40. 9 62.7 40. 9 62.1 40.6 62.2 40. 4 62.1 40. 2 62.8 40.8 62.7 40.9 62.6 40.7 as 123 0a #4/a 1.5647 4625 de ds 61.9 40.1 61.8 40.0 61.8 40.3 61.5 40. 0 61.6 40.1 61.5 40.0 61.1 40.3 61.1 40.1 61.1 40.1 AVG STDEV N 1-95 26.26 28. 45 28.45 25.79 26.32 28.47 28.82 26.45 28.80 26. 40 0.262 27.07 27.05 28.52 26.45 28.47 26.45 25.18 25.50 25.50 26.24 0.853 28.40 0.75 38.00 1.7823 166 E. PRESSURE-SATURATION DATA OBTAINED USING PRESSURE EQUILIBRIUM METHOD E.1. Air-Water Data - Metea Soil. E2. Oil-Water Data - Metea Soil. E.3. Air-Oil Data - Metea Soil. B.4. Air-Water Data - Ottawa Sand. E.5. Oil-Water Data - Ottawa Sand E.6. Air-Oil Data - Ottawa Sand. 167 weenie.— hmedm mvoKN oNS Ax. mad. 3“de mmhd— v 3.:va whmwo :oKN mug» o\o VNA _vwdm mm. .3 cm 32.3; Soda mwcdm Nos. e\e A: .N- vNNdN mood— N 333.....— ooM.wo mmoém O05 .X. 8.2. wmmdm Maw: _ __ ASV Av 87 AV ADV A:v sz _— men-.2— :ow .583 633 3:23: .855 .5235 Oz 5:5 .5 ... 83 52:33 3.28: 8.2 88:5 855 _ «Sm—2 a 82. n Av a 8. a I As. a :88. » Av H :E: 8 .5.“ KO 262 Ede one: mod: :..m 36% wad? «Econ 36% 3.63 v 3.2. 3.3.. ood: nod muwmm «can... modem 8.8m 8.3m m 3.9 3.45 8:: and mod? madam 2.an 3.54 85.3 N h? 2.4.: and: «En ntmnm 34mm 2 don SE 3.43 _ on AV AV on Amv AV Amv Amv Anv AV 02:. a. 02:. a 5.5% 4 cups. 02:. ...8 £8 Es:— 022 or» .52.. ...: Sim RODS. Ess— Eo zs>o .5: 3.555 o: E REC—H3 5H0 <88; :8 as. 28% 168 c:.m.. 03.2 36.3 :36 «8.: nut: Ema- 3%.: ~32: emnN. 35.: 2nd— Nvoé mNnd— ohms owed not: Rod «3.7 and— a—md— wood. «mafi— Caed— hovd «8.: 52 .2 End *3}: $8.3 mmmfiu ciao— nmNd— nmhd 93.n— cond— NQCN #2 .N— «3.: vmvd we”: $3.2 966. 8a.: 29.2 owed 35.n— need— GNBN eta: 29.: 966 news: 35.: had. 33.n— mn—N— n34 Newsw— mm»;— cvud 03.: wood— wad 33.: Sad and. n2 .2 n2 .2 New; ”at: mend— nowd unad— gad— emod mafia 20.3 $8.». ”—5.2 8N: 8m; 85.2 shod— omoN 35.2 thS wood no.3 52.2 35.». mwNé— mend. ham; shofig nag: o:.m 392 08.2 :66— mn~.~_ mnud— chad. 48.2 unad— SQA ”—5.: .3.: E.— .m 89.2 find 2 _ .3 god— whod mad. wnwd_ mad n~n.— —m~ .: w: .2 8nd #2 .c— omN.n_ 2: .2 mad— eohd— 35.”. 52.2 ”Nuan— hmm; $5.3 anod— 03 .m smud— 02 .3 3.4.2 85.2 no.6.— Neadi anew— habé— 8m.— nmod— hund— enmd N: .3 mend— w— _ .c— cues: 8a.“: 35.». 3%.: 23.2 cmn._ ”med; owed— m-.m 3&2 enqu— nSd— coca: wad— soom. SYM— _nm.N~ Kr:— owed— owed— Nmmd anod— wu—A— NEd— wad— 30.: 296. :02 ZmA— One..— cwod— 85.: gram ”3.: was moms: 3.0.: hand— ne—d. :0: 3:: 8mg 85.: need— avd nmod ”wad Namd— homec— homd mead. mm— .2 «had cm..— nwedu hand Coed nond— ohmd- and .25 85.3 53.? ._o> 8.0.3 58.3 =e> 88.3 53.3 ._o> 96.3 .83.? 373—92 mare—_mz $33—82 gig—m2 169 .352 3.33. $.28. SS2... 5.2.8 .53.. 82.3. Sande. 2.. an an 2.... 2.... 2 ...... 2.... 2 3.. oh. 2 2.. 2.. n2 ...... 2... n2 2.... 8.... an 5.... 8.... n2 2.... 2.... 2.. ...... 8... 3. a... 3.... 2.. 3n- 3.... n9 :..”- o2- n3 8... 8... n8 n3- «.3. 2... 2.... Rd. «S 2.... a... 2.. 8d 8a 3.. 8.... 3.... 2.. on... S... 3. 2."- 8.”. no 8... ...... 2.. 8d. :..N. no 2.... S... «S 2.... 3.... 2.. :.e .3 2.. 2.... on... 2.. 2...- 35. n9 2.. 2.. RN ...... 2... n2 3.... ....e. an 8.... 2.... n2 3.." an.” an. «as a"... 2. 8.... 2.... mm. 3... an... 3. n..." S." 2. 2... :... nu R... :..... R a... a... up 2.. 3m 2 2.... 8.... 3 .....- .2.. 2 a... n5 2 .> 8......» .— AM...» 3...!» .— dfl$ :..; .— flvg :..; .— 2328). 91.53. «3.5:. ........m2 .833 170 cod and _mde 8.0 23.6 Sada 8.0 «and QNNdN 86 00nd wandu ~0.N00 ~86 Nnod 5N.05m god anmd 3.5m 306 93.0 5— .35 V006 35.— N— 60— 3h _ .c 3.06 50.3..” 2: .0 New ._ 2 .muw I: .o a: .N 00.5m _ awcd. 0m0.— 2nd: 3:... n05.~ 3...:— 526 wand 8.9.— Nm—d wand $62 n26 §.m 3.3 90. .c n00: ”~62 mm— .c 984 mm .3: c: .o 0&56 8.0m 52 .a *9..— nm.5w n5. .0 3.56 :.3 «3 .o ”on; 8.2— _2 .c m5~.~ ”~55 a5— .0 5:: .. 5nd” 0w— .c 93.0 Nn.n0 0: .0 2x: 05.8 o: .c cooN :60 *2 6 03.0 365 can... and 3.3 ncnd 80.0 860 53 .o and .3.—0 53.0 80.: «man Nefld «and wnfin 320 an: mu .3 nnmd 30A :60 5am... Nand 5&3. nnnd 03H.— 5— .2. EN... 9x: 20* c0~d 0:.— 2 .mn 3N6 52.0 3.? 05nd m8.— 0a.5n 8nd nmad 8.3. 8nd 3N; 05.3 3N6 m3.— madm «and wand mm .nm 3”: «2: 3.3. 3nd $0.2. 3.3. 3N6 02 .— 503. 3nd «.2 ._ 2.0..” 2.00 :00.— m— .3. w. md 85d 2 An and «and 00.3” and 30.0 50.5m mmnd 25.: ~53 and «5.0 2 An vamd a3.— ~5.5m vmmd o3 ._ N— .Nm 036 N50. 8.3 :66 034 an .5. $86 ace.— na.5N Nmmd ~86 519“ N02. Q5N._ wvdn 5mm... 25.— 8.3. and owad 05.0N 002. 8— A mndm 3nd 5!: «mm—N N56 «2 A 3.3. ”mad n8.— 5N.: man... 53.: 3.3 05m... 306 no.2 and ncad no.0,“ wand ”and 8.0. 55nd 8.0. 03.0 . 09°. «and acd. 00nd _1 9.— 9—2 30 9.— UE 3.. on U: E 9.— 02 30 —1 3T3 a: man: a: Ntui _mz :73 _NE 171 0m._ : 0N.¢0~ _560 N000 00.00 0~.N0 5— .00 00.50 0 0 .00 :.00 50.50 2 .3 50.0m 0 0 .mN 00.2 00.0 3.0:. SN: ,2... 0N.0 2.0 and 00.0 9.0 _0 .0 m0 .0 «0.0 00.0 «5.0 05.0 N00 00.0 00.0 00.0 00.0 00.— 5.2: w - $.02. 8.03 00 .mm— 00.0: 05.50 00.00 00.00 0m.0.v 8.: ~00.“ 0— .0... 5.00 00.0m 0000 00.: 00."— 00.0 00.0 00.0 05.0 05.0 ~00 90.0 50.0 3.0 00.0 00.0 00.— , ~33 00.50 0 3.00— 00.00 0N.55 :.00 00. _0 00.00 2 .m0 05.00 3.00 0— . _0 2.: 00.50 00.00 00.00 00.00 00.0 5— .0 00.0 :0 3.0 00.0 00 .0 00.0 00 .0 00.0 00.0 «5.0 55.0 00.0 00.0 00.0 «0.0 n— m a w : 3:082, -- 2.5% w, 172 Ex; .3.: 23h 3m .x. 25 ~83 892 MI... 83$; «:2. 22.: 20 $2.. SEN 3a.: .m 323; 82,... E .Ma :.0 .x. a..- and 95.: N 80:; 92.2. :02 8.0 .x. ”2 M52 35.2 H __ AEV ADV sz AV Adv A:v Aev £2.90 :00 .855 .333 3.3.—am .505» .5035 Oz v...... .9 ... EB .35.? 3:28: 2.2 38:8 .352: a a: u Av a 2 _ ,w u $0 a 98.2 n Av “ 5:: 8 .8 $5 302 and. 05.: S; Sam 0. :.fi :3: 8.8“ 23% n 8.? :.E 8.8.. «2.3 N Si 2 .2.. 8.9:. 8.23 m A0v A0v Amv ANV Adv oz:- Q Uzi A—Om GHMU<0 ESE >¢0 ZED Oz z: _ page? .58 H ‘ 4MP”: :00 3%. E033 0 : a: 173 'IV a .84- Sad— Emd— n33. _me— 3.9.2 :1? 35.2 SQZ chad- wand.— ~mn.3 52..— find— Sim— :wd ”no.2 3b.: head. 25.: 5%.: gram. «are: and— m—~.N 63.n— CNI Sc.» 2b.: m:.: 95.». some: 92:: find aqg cred nan.” EN: :9: Quad 39.2 shud— Eudu 2.12 8a.»: can; 9.5.3 25.: ”bed 20.: chad and 63.2 ahwd— hand. ”8.2 emad— onm._ 58.: «8.: :.wd mad— wofiv— ~36 $5.2 2%.: 3nd. eg.~_ SN: _mm._ 35.: p.96— mccd 80.2 2 “.2 ”cad 2%.: ~21: 39w. nnn.2 I _ .2 c5..— mnfa 3nd 8:5 and. 3&2 BAZZ $9.: 292 god. ”3.: v3.2 ~21 cs».«— :1: 2: .m and— ehod— 3GB 2.3“— bwaA— «nud- Vncd— :wN— nun; :1: Sma— wS .n 2.0.2 cad ~86— hwaA— $5.3 mund- Zwfl— 58.: «S.— 3nd— 93.: 2:5 Sad and 2 _ .3 36.2 Sad $3.»- 53.: nmcd_ 8w.— vno.: N2 .2 S _ .m N34: 29.2 _:.S 3nd— _on.2 53.». 36.3 _N~.n_ «we; mucd— _S .3 ”2 .m 20.3 Ebé. 2: .2 ENE 3b.: Maud. _N~.n_ an”: 20.— _~—.: wad— Xud 25.: 03.2 o— _ .c— 35.2 n93”— mmod. «3.: 3a.”; «No; ”3.2 3mN— :2 .m 30.".— wmflfl ”nud— mmin. 89$ Sad. :52 9nd— hmfi— de— 53.: 3nd 282 ~36— ?md— 39"— mvn.: :36. 3.92 36.: 63.— 53.: «find— 2n.m NZ .c— 036 ”3.3 were: _2 .3 Sad. 03.: _NNAZ chog vend- cad Zc.m nurse— cbnd- 3N.~ «:5 65.3 53.? as 95.3 .38.? =o> 20.3 933.3 «:5 660.3 53?? «tug—m2 Qua—m2 gun—5E4 $.2—m2 174 ,___ ._...____.— _._——— ~——_~—— 3.2» 2%.»? $.28? 83mg- ii; --,,§N.$» :22» as»: has»? - a.» a.» 2 2.2 2.2 a :a. 2.». n a»._ e: a N2 9: a» 2.» 2.» 2» 3.». 8.». a» :3. 3.». m2 _ 2.». n3. 3. z.» 8.» n» 3.... on... 2.. »».»- 2.»- 3. . fl». oz. 2» 8.» 8.» n» 8.». a». 3.. 2.». F». a» _ e3. 2.... a» 8.» 8.» a» as. as. n» on». 8.». a» :..»- 8.». n8 8... :... a» 8.». :.N- a... 2.9 8.». “S w 2...». 35. 2» z.» a.» 2» 2a... a». 3. 2.»- 8.»- a» ‘ 2.. 2» n2 2.» n.» a» 3.». 3.». a» :3. m»? an \ 3.» »».~ 2. a.» a.» n2 8.»- 2.». n2 3.: «3 n2 _ 2.» S.» 2. R.» 2..» a R». :..». fl 3.. »_._ n» n a.» a.» a 2.2 8.2 ...» :..»- .2. a 3.. »S a _ .... .> 3:; mm; 33$ .— £24» 3.3, .— ues 53$ .— - :2sz 2.252 :..»sz ;.»:m2 _ seesaw 175 #382. Need dddd dwmd dddd zed Snd owed awed _ n dwdd Qdmd 356% Nwdd mdmd 35.3.0 mmdd was; Nmnfion hddd QnNd mine“: 33d «2 A Sndnn owed add wowdnw 2: d 333 933mm 9: d 3N.— nhmfin— 9: d n2; wwdfimm ”add Sad Sedm— 2: d 55.4 dad: 2 ~d add $33 ”2d wood vafin— dd_d ”and 2.0.2: 53d 3*..— nn—NQ _3d cad «3.3 $2 d 3.5..— owfid: d3 d «9..— hmogd awn d N3.— Noméo mmNd dmmd «3.3. «m_ d nun; 53.8 _S d ddnd 8de ndmd nah; naméo emud «mm; 02.3 Nde and; dandn ”end 3 m; 35.; and a: A gain wand dmwd mafia w— Nd :_ ._ 805* and 3 M." wdhde nVNd c3”.— _nd.nn demd nae.— wnmdv nmmd new; n5”: _cmd n8.— o—«Nn momd nvwd nnmfic and cod; Sade 8nd 35d v3.9. find ”2 ._ ZN: had 3&4 mnnév edmd hand wnndv :.«d 3.: _mvdn gNd wdwd nhwfiv omNd 3n; 3%.: End Nan.— aoddv hand awn; vnmdm A..dmd 8m.— _n_ .3. M: md cw»: nnndfi ammd Qddd 3d,: bdm d «dad :.va wand 93d 53.9. den d god N593. «mm d m2 d 53.3 d~m d #3.— nhnfim _mmd 3m; _ _d. 3 nnmd and 3b.? enmd dam.— Nowde nmmd wnea deodm _nmd Nhhd dendm anmd d3.— beném 55nd «and dmcdm enmd 02: and— Nomd dm— ._ dbddm and Rad nnmfin nwmd n2”.— cdmd hum d and; wand find no»: N: .2 and «an; ched— eded v33 dhnd. mdvd wand- gmd ncdd. zed mend 03d on 02 3e 9.— 02 3% 9.— 0—2 39 0.— 0—2 36 STE—m: 372—82 NtfiZm—E 176 ....h... .r,’§¢.§.._.. 'I.‘ E «as... ..N... ~23 a... 2... .22 a... 222 8... 3.3 a... 2... 8.5 a... 3.5 a... 8:. ..N... a... 3.5 a... 8...... ...... 2.3. 2.... a... 2.2. a... 3.... 3... 3.3 a... 3.... 8.8 3.... 2...: 3.... 2.3. a... ...... was a... a: a... a...“ 8... x... :..: ..n... :.s a... 3.2 8... an d an. _ v bod mnfin 3d enfi¢ mod 3... 8.9. :... a... S... 2.3. a... S... 2...... a... :..: fl... 3.2. K... ...... a...” as... 2.2. R... 2.: ...... on... ......” 2.... 3.2. a... 3.3. 8... ...... 5... a... :..... 3... 8.2. 8... 2... 8.8 8... a...“ a... 2.: ...... and 2.2 mad bddn «dd OMHN cad 8.. ......n a... 3.2 a... $2 a... 8 m 8.6 8.— 8.° 8.“ ©N.O 8. _F $.52... n. m n. m n. . 8.2m: -. -- , «.3.sz w J— ; Fem: 177 a... .2... .3 . ASv- AV "Aflv mac; n EYE—u hnocwcau “302 voao: ~ 02..an unfin :.e .x. 8.7 2.5.6 Exxon 2.3 5v ammo: — Nada cnfiu «NS e\.. :.o wnAhn mm.:.n _Ndu in via; _ 2th :..cn Nvfi. .x. No; N¢.—hn 99.2% aw: cw ”an: — 29g 8.2 36 .x. a N- modem ~n.:.n hndn fl AEV Afiv ADV AEV An—V Av—V An_v ASV __ £32. .3 5.3 5.3 8.3.... «5...: .5... .... 5...: Oz ...-a ... .5»? .33.... 1333.. Is. ...: nosing-8 38.—5 :.. .2 8. .2 2 .2 2.... ”a? “a“ 8,: deS onfig Sum ”93% Nwdom :.wom mmghn hmdnn Sign No.39. a. coach and: Sue—N wed NOAMM undo“ mm. :m vadhn 590mm mofiom : .No». n M: .2 and: mate—N ocfi Sammm endow sign coax-m 2.60m mw.mom wade... N 3 db CON: 3 .EN Qua no.3” :..oom modem mwéhm :..mmn nwdom wads. H A. .v A... v on va Av Acv Amv Avv Amv A~v Av .224. a wad 3.5% a ...o ...o .2 a. 32:. 28 28 >22... Uzi o... ...: ...: ES. ...... :55. 5.5;» ...—2,2... E2... .5... 2.30 :2. 32.5... ......o .....mwaj Oz Ely ENE; 33 2.0.3 :53... ....> 2.93 can...» ....> 26.3 5...? ....> 2.0.3 5...? ......2 ......) . ~22.ms. $2.8... $2.2. .\ w s_—_ 179 2....ch 58...... 2.32.. 2.2.2.- h 22.2383... .252... $23... nun-"fl -- nil l . If}! I! illlll iii- :Iitlltllliiu- - 3 l2. Inlllil. Isl! IIIIlltllvllllvlliIuil! 2.... 2.... . ...... R... . . F..- .....a. .. a... 3.... .. ...... ...... .n ...... ....a R 8...- 2.... on ....a 2.. ..n R... 2.. ~... :... ...; c... :..... :..... ...: ...... 2.... 8. 2.. :.m ..2 ...... 8... 3. 2...- 2.... 2 Rs 2.... ..2 $.~ 2..~ ....~ 3... an... 8. 8.... .3. ...~ .....~ .....~ ....~ 2.. 2.. ..2 8... 2.... ...: 2...- 2.... 2 R... ...... ..2 ..n. a... 8. :... .2. n. :..... 3.... o... S... a... 8. ...... a... 2. 3... ...... ... :... :..... 2. ...... 2.... 2. ...... ...... ..n .3 8... ..n 8..- ...... .... ...... .2. ..m 2... ..~... e~ ...... 2.... n~ 25. .....a. ... 2... 2.... .... 2.... 3.... .. :..... :..... .. F..- .3. .. S... .....o .. .> 58$ .— 2..» 2...; .— .> 3......» .— u8.> 3...; .. 2.2.2). 2.2.8,. - $2.2... $2.2». :.....cwiu 180 89¢ vote Nnn. _ N «fin.— vcvd 2 N6" ancfi «omd gNAN 896 End :.ndN ..2.2~ .2.... «2.... .2....2 S... .2.... ..2..~..~ 8.... :~... .2...? .2...... «.2... Nch: 2: .c Kn; anus... _2 .c ”on; 25.2: 366 N36 3; ...n *N— .c 63 .m .3? «~. ... 2.«... «3.2. «:... .2.... 3...... 8.... n-.. .~...e~ 2.... 2...... «Guano o: .o ”we. ochdv an m .o mam; :.mdm n— u 6 8m.— ”:..QN a: .c thd 2....~n .2... ..««.. 2.2.. 2.... 2...~ .23.. 3.... «8.. ..2...~ 8.... Sn... 2 _.wm ow— .c 3w; "3.3 acmd one; Neda 93 .c Vncfi bafnm R— .c Kc.— cohd~ 2 N6 mafia m: _ .NN nnnd «a .N ahndn Q3 6 mam..— n—an n—Nd 856 592mm mend hm”; _Nad— moNd can; ”$de nuNd §.— ~36— mmNé 2.: .N _2 .NN NhNd 25.— Mde— wad :N._ 3.0mm ode n2 .N Sad— can... one.— «......2 8.... ~««.. 8...: 3m... .2.... «~68 .«~... .8... ~82 ..«~... 8.... nhN.w_ 2 m6 «Nod ”~92 3 m6 unnd SQON namd 364 8a.»: cam... mad RY: hand m8. vad— _mmd n2..— onndu 3 m6 «.8... «m9: ”and Sad __ .32 ~..«... «8... «2.: 2...... 8.... «2.2 «mm... ~8... ..«~... -n... 25... 3...: ..«~... 2.5... .292 8.... 3.... ~32 «em... .8... .S.« can... 28.. «2.2 .2.... ....~.. 3...... n2... 2.... ~«..... Sn... «...... ..32 can... :... :....2 «.«n... ..~«.. .2...... 2.2... m««.. ..m. ... Rm... «8... . . . .. .. .... ...... ...... ...... m... .... .....2 . ms. 22...: 3.2.9. 9.. u: 3.. 3.2 $2 _ _ _ W _ _ 181 r r. fiifim‘i’icr 2..2~ -... 3.3. .~... .....~..~ ..~.. 3.3 .~... .. 3.... .~... 3.3 3... «5... .~... :.3 3... m ~22. 3.. 3.3 3... a... 3... 3.3 2.... _ ....S 2.... 3.3 3... 3.3 3... ~3~ 3... ....3 3... 3.: 3.... «32. 3... 2....~ 3... ...«m ...... .«.3 3... 2.3 3... 3.- 3... 3.3 3... ~..- «3... «m.«~ 3... ......~ 3... . ..3 ...... 32 8... 3.3 3... ..«2 3... 2.- ...... ~«... ...... 3.3 3.... «..2 «.... ..2 .~... 3.: 3... 2...~ 3... ~22 3... ~.«. 3... 2...... 3... 3.3 2.... ....2 «S. «...: .«... 2.2 ~«... 32 .«... 3... ~«... 2.2 m«... .~.... ««... 32 n«... 3.. n«... ...: ««... 3.2 3... 3.2 3... .....« 3... _ .~.2 ...... :... 2.... «...: 3... 2.2 8... . ..~. 2... 2..« 2.... 2... 8... «2. 8... ,. .... ...... 3.. ...... 3.. ...... ...... ...... =Ha¥2 m n. w a m n. w .2). 3.2 .22 3.2 .ms. 1: 182 5o 9..-: . AVA? " Am—V «S. n :3... “2oz (42% <\Z¥ <\Z% <93 (\2% <93 <\Z% <\Z% v nahvvva a—vda 2.?" S.” e\e had endow 7.60m hmd— m 03.3. _ umndm .vnfiu QMN .x. mud 2 .won 2.50.». 8.: N 30va _ ahvda andN hmé .X. and 3&3 mason _Ndn — Aflv AIV Aflv ASV __ 831a «a... . a... 5.3 Oz I: '43 38:5 52.0— 526. Ed'— ood ”flab H u gig <\Z% <\7$ <\Z% <\Z% (9% <95 howcm wmdmn on. _ on 925* v 02.9. SN: bud—N and NmAmm mmfiom owdém woghm BEAR vcfiom :.Now m “3.9. me: 3mm: wage. NhAmm :.mem :.wom :.th modmn endow 26.3. N 2.2. _mN: OWEN hwd demm mofien 3.43m Nut—hm aadmn :..mcm NWNO». _ A. _v A: v on va Ahv on Anv Avv Amv ANV A_v cad cm...“ a .3.—2% a [a qfl. ..oh a 32:. ...o... 4.8 ELEM $22 a... ...: 2: ES. ...? 5:3 :..—.53 8655. Es... t... 2.55 :2. 37.35 ...—09.5555» Oz F: mun; .....6 5....-.» .. ....-.» 3...; .. we.-.» 2......» .. we-.. 5......» .. :.....ms. ...-:.mfl ...-2% ...-2.2.4 5.2:...5 $2.. $2.. $2.. 2.... 3...: 5.... ....~.... ~32 .~.... ..3... 3...... 3.2 $2.. $2.. $2.. ...... ~32 32.. 3......- -o... .3.... 2.... 3.2 33.2 $2.. $2.. $2.. 3.... 2~2 .3... .82. .32 ...~.~. ~.o... 2.3.. .....2 $2.. $2.. $2.. 3...... 3.... 2...: .... .2- ...~.~. .3.... 32. .....2 32.. $2.. $2.. $2.. 2...... .8... ..~2 .32- 3...... .32 .~... ~..~.... 2.3.. $2.. $2.. $2.. 32.. ..~.2 32.. .32- .32 3.,... 83.. 2.3.. .2.... $2.. $2.. $2.. .32.. .3.... 2...... 3.2- ~32 2.... .~... .....2 .~...... $2.. $2.. $2.. 2.2.. 2...... 2...... ......2- 2.... 3...: .3.. .2.... ...... $2.. $2.. $2.. 32.. 2...... ...~.~. .22. 3...... 3.2 .3.. :..... .....2 $2.. $2.. $2.. 32.. ...~.~. ~32 .32- 3..~. 23.... 3... :..... .3... r .32. .~.... ._o> vac-3 28.-3 20> .2...-3 28.-3 20> v5.3 28.-3 E ...... .ms. ...... 2.2 184 <9... «find. :5...— 3N6~ <9... <9... <.Z.. 23.3 8. .c 59.... an... 2 we. .c n3... 3.2.3. *8... can... J <9... <9... <9... 3...: ~26 no... 293 3.6 an. owodw 8.6 «:.. <.Z.. <9... <9... ~5an a... .c o. N; eondn em. 6 ha. A 3a.? «a. .c mum.— <9... <9... <9... 3...: 8. .c and SW? .2 6 2b.. acqmv N: .O can.” <9... <9... <9... N. m.mm N: .o 3;... 356m an. .O mm _ .N 02.6.” an. .o sad <9... <.Z.. <9... _ .m.n~ «a. .c ”and 33.3 haw... cg." panda «MN... Sufi <9... <9... <\Z.. $5.2 3m... «3.. 593 and mad 893 find 2%.. <9... <9... <9... _n. .3 oomd awed $53 and on...— omodu 8nd 3.: <9... <9... <9... .gdu wand ENN 3o.- ham... N. Na as." mmmd ooh.— <\Z.. <9... <9... .3.: .26 one. chad. own... an. .N 8nd. «mm... «3 .N <9... <9... «$23 nnmd n 3.2 chm... «and can... 9.— US. .5: um 02 3% em 0!. E 9.— US. .59 . ; $1.2m: - ,‘ ...»... .2: ...-.....mi. .2...... 185 2.32... 2...... 3m... ~12: .x. 3.~ ......” 2.3 ¢ 232... 2...... 3...... 3.,... ... .3...- .m... ....E ..m .33.... 3%.... .~.... .3.... .... 2.2 .~.... .~...m .~ 3.3.... .33 3...... 3.... .... ... ... .....3 3.3 .. A..v AV sz ADV A? A? $7 __ .5 35:5. 02 ...... .. .. ...: BEE... .5623. :..: 3.8.30 .2...... _— a 3....- " Av .m 3....- » AV ... ... . .. n Av . 5.5.8 .... .25 082 co.” m. .2 fl“ "n u 3.2. wniL om...»— mwd ...-.33 3...? 86% «5.8.. e 8.2. 3.3.. 2.3. a...” 3.33 ...-6...». 3.23 8.3.. m and. 9......— n. .3. SM mNdQn and? 2.8... 3.5.. N RN» «:..... 3.3: end #033 .2...-n 25.8.. 3.8.. . on va Ahv va Amv Avv Amv Amv Anv A.v uzi i 92.”: man—...... a 3:3. 62:. do. .28 Es... oz... 9.5 .52.. ...... -35.... :85...— E2... .55 286 ....2. 5.23:. oz 2... awn; ...-mu " E I" «N... 235 9.0.3 28.35 =o> .35-3 28?? 20> v5.3 938-3 33.2% Qua—ms. gazm—z :..-«:9». 187 «N6 34: end 3.: 8.3 nmd :.c $0.9m 36 3.6 3.3 3.9 3.0 end” end and and" cod 8.: 3M _ :..c mud mud «56 3.: N— .a had 36 8.: mad 8.: Ed 36 nod :..v «ad 86 84 86 8.— n. m n. m Q?— sz 3?— 5m; 8e... ”and 8: m Sad ”and 83m . 3.5 condom 2: .o 93d 2 n62 mad :3.— meadnn Sad «36 5:: e36 056 :53 um _ .o :66 «8.3. 2.6 «2..« ~nm.n2 waod as .m v8.3 22 .c 3.— .N «nude E: .o cg." 3a.? en _ .c than 2%.— m 3— .c «and 50.3“ $.— .o 3a.». wmode m2 .o 9: .m 58.3 mm. d entm 2.83 3: d 2 ~.n 03.3 «3 .c 3.— .m 3%»... 2.— .o find End— Eud So." 328 End an?” 893 Naud Zed mm _ ._ m 2 ad .3." EN: 3&6 www.~ 2.— .Nn amnd «nod «26" find cad $0.2 wmmd 3...... 3a.: ”86 03am 2 _ .2 send 3.— .m «3.2 uwmd a. _ .m 93.2 «and k: .n 3nd _ and and 39.2 and weed wvwd 2 n6 can...” 3&6— nmnd coca” on»? and mvvd 08.2 3nd maven 8— 6 3nd mam .N 32: «mad 3mg chad 3nd 304 gad— Smd «2 .N 3.:— 32. 8nd 8? 32. m3. :2 at”... 2 _ ._ 32: 28 n8.— noad 82. 9.2" $3 2.3 2.2 on: 83 2.”.— X: 83 . 9.3 :3 3.3 :..... 83. and 83.. 82. ”8.”- 82. 5.1 8.8 2 >6 o.— 02 3.. on US 3.. 9.— QE 3.. u.— o «33...: :1: :..-E L is; 32: :5 (<23fi (\Zn 412% <24.“ 4&2“ <95. <\Z* 1 33.8.— No.m: 8.3 no; .X. w_._- NNMN woNu m vammbw; 24.: MNQN mm: .x. 05.0 c~.- wfium N <\Zn «:2» 32% 4&2», 4.2% <\Zn $2» _ __ AEV Aflv sz AEV ASV A? $7 = Eng: :8 .533 .5235 35:2— .525» .555 he 2— 33 =0€33< 1323»— 852 Exec—BU 62:95 a a: n Av m a 98.” n 4v .5 3: n Amv ” 5:: 8 ..8 35 8oz Efifl <95 <\Z% <\Z% <\7$ <95 Q 3.2. mada— 2 NS :..ofin ooficv m 8N5 Neda. Sana— nQSn ctmow N <95 <93 <\Z¢ <\Z* «34* ~ Aav va Ahv Amv ANV AV a. d 9mm _q 62:. d8 58 :85: E2m an: 286 :2. 9. LE FIG—H3 AAHU <>»<.—.O :5 >3 5.5 2an 189 <93 <23 <93 <93 <93 <93 <93 <93 <93 <93 <93 <93 <\Z3 <23 <\Z3 <93 <93 <93 <93 <93 <93 <93 <\Z3 <93 <\Z3 <\Z3 <\Z3 3.3m INN— wNa.: 89.» ”3.2 353.2 <93 <93 <93 <93 <23 <93 och." «Na; _ ”2 .2 9&6 ”3.2 «and. <93 <93 <93 <93 <93 <93 coo.— chN— «8.: 93.6 20$ 2a.: <93 <93 <93 <93 <93 (<23 93.: «8.: «Nad— 934 2a.: SN: <93 <93 <93 <23 <23 <93 can .0 m mm ._ _ awed OX: 93.— _ n _ ed <93 <93 <93 <93 <23 <\Z3 and ”3.: 85.2 Sad 29»— cEB— <93 <93 <93 <93 <23 <93 can .c 35.2 van .2 356 85.2 «2..«— <93 <93 <93 <93 <93 <\Z3 can .c 3n .2 The.” cued «3.2 «mad <93 <93 <93 <93 <\Z3 <\Z3 ovmd had— oMo.n_ can .c n2 .2 :.Nd— <93 <\Z3 <93 <93 <93 <93 $36 98.2 62 .2 and :Nd— 3&2 <93 <93 <93 <93 <23 <\Z3 Sad 62 .m _ 33.: owvd 3&6— NQQZ <93 <93 <\Z3 <93 <\Z3 <93 8N6 03..— . 35.3 and Na..— _ SQ— _ <93 <93 <93 <93 <93 <93 :3 .o 35.2 39.3 and 93.— _ 3.3.6— <93 <93 <93 <93 <93 <\Z3 c3 .c «wed— ged cad 34.2 586 <93 <33 <93 <93 <93 <93 c3 .c awed 2 m.” c _ ~.o 586 can.» <93 <93 <93 <93 35.: >86 <93 .75 95.3 «33-3 35 9.0.3 53.3 =o> 93.3 9.835 20> 9.0.? 986.3 3.3 50 3.250 «3.350 3.250 190 ml antral-I LP, 4 . 8...... %_ 8... o 8.... c 8.. ...... c ._ _ ...... 8.. 8. 8..... 8.. 8. ...8. .... 8. 8.8. 8.. 8. A _ 8.8. 8.. 8. :.8. 8.. 8. 8.8. .... 8. 8.8. 8.. 8. _ . 8.8. 8.. 8.. 8.8.. 8... 88 8.8... 8.. 8.. 8.8 8... 8.. _ 8.8. 8.. 8. 8.8. 8.. 8. 8.8. 8.8 8. 8.8. 8.. 8. _fl _ 8.8. 8.. 8. ...8. 8.. 8. 8...... .... 8. 8.8. 8.. 8. m 8.8. 8.. 8. 8.8. 8.. 8. ...... 2.. 8. .2.... 8.. 8. m _ 8.8. 8.. 8. 8.8 8... 8. 8.8 .... 8. ...8 8... 8. _ _ 8.8 2.... 8 8.... R... 2. ...8 8.. 8 8.... ...... 2. _ u 8.8 .2. 8 8.8 8... 8 8.8 8.. 8 8.8 8... 8 m _ ...... 8... m. 8..... ...... n. S... 8... n. 8.8. ...... n. W _ :8 8... o. ...: ...... ... 8.. 8... ... 8.... ...... o. .. “ ...- 8... c .... 8... o 8... ...... a 8.. ...... c _ - .. 38$ .. Jr... 3...; .. ... .. 5......» .. L8H. 3......» .. _ 8.8.80 8.8.80 8.8.80 - 8.8.80 9888.10 191 <\Z3 mhudm 83.3 3nd 3— .NN <93 <23 <93 <\Z3 <93 <93 <93 <93 <93 <\Z3 <93 <93 <93 <93 <93 <93 32 Gun . quad can... 58.83 one... can... <93 <33 <93 <93 <23 <33 $5.6: ”Ned End. 3363 and 32 .c <93 <93 <93 <93 <23 <\Z3 _ ~36... _ 33.3 83.3 @862 «8.: 3 m6 <93 <\Z3 <93 <93 <23 <93 mnwdw 386 93.— 3565 n36 m8... <93 <93 <93 <93 <93 <93 3 _ .hm nncd Que.— _ 2.3. 33.3 8....— <93 <93 <93 <93 <93 <93 32 .n m 2.36 awn; _mvdm whcd «Q..— <93 <93 <93 <93 <23 <93 mm— ._ m 33.3 «2 .3 3 3.3 83 .3 «mod <93 <\Z3 <93 <93 <23 <23 N..— .eN mm. .3 mg. m 35.3 N3 .9 con.” <93 <93 <93 <93 <93 <93 nomdn «36 Sad 32 .2 82. «mud <93 <93 <93 <93 <23 <\Z3 on Now 936 Man; 3n .3 8N6 hncfi <93 <93 <93 <93 <93 <\Z3 ~N~.a_ wand 3o.— _nn .2 can... 80.— <93 <93 <93 <93 <93 <\Z3 unma— Nmmd hmhd 33.: m 3 m6 «and <93 <93 <93 <93 <93 <93 aflo— nm m6 need an.” man .c 30.9 <93 <93 <93 <93 <93 <\Z3 33.: m3“... and 3%.... _mm .o 33.6 <93 <93 <93 <93 <93 <93 unfin— _nm.c mad 33." men d 33.3 <93 <\Z3 <93 <93 <23 can .N 3nd Nuns. . onnd <93 <93 9.— m 3:. on m 35 oh m 39 u.— m 3‘ 33.350 3.350 [\l «3.350 3.250 192 2 .8m 590 23.83 .86 ohdhn wad Edna c— .c N352 2.6 8.32 26 3:3 «ad 35.? n— .c 3 Km 3— .c 2.3 2 .o 2 .nm :6 ”3.3 26 2.3 end 36” and w— .3 ~36 nhdn 335 OWNN and 2 .2 and nndm 33.6 3.3 who find— 36 3.2 3nd mNd— mad 8.: 86 3.2 36 a.” mad 2.33 336 3&6 nod 3.2 *6 and 36 8.9 cc.— ood 8.3 A m m m A 333:..w0 ".33—$0 «332.8 33278 193 ESL: 32.3: Run 2% °\.. 3d. 832 52 .mm 3 3505.3 bani. 33.3 o~.~ o\.. 3.3 2 3.2 33.3 m 83.5.. 32.3: ocmsm 3n.~ o3. 3.3- 33.2 30.3 N 33$.- 33 .3: 03.3 3d X n. .o w-._~ 33.: _ AEV ADV A3_v ASV ASV A:v AEV _— .855 3.3.—um 3232< .333: 3a: 183:8 32:95 a a: u Av 3 e8.» u Av a on. a n Av H 5:: 8 .8 as 082 H H" 3.2. we c2 «03% can an? 3.9% 36% Rd; 2.33 3 SE. 332 2.3— 2am 3.32 3.2% 3.3% and? 3.33 m Si. 363 8.8— 2:“ 3.an 8.9.... 8.3m 3:: 2.83 n «3.2. «333 3.32 and 93mm oflnvm $53 2.2... 3.83 _ on AV Ahv on Anv A3v Amv Amv ANV AV Uzi d 022 i ”.~...—24a a nut:— ozi £8 £8 Esm— 022 9.5 E: ....< 55% :85: Esfi 35 29:5 Ez— mflzgn Oz z: Emmy—H3 1315.5 n‘ H nz 95.3 .535 :9 95.3 58.3 8.880 , 8.880 «.1880 195 $2.2- 52.2 522- ~§.8_ 3.? £8.34”- ~83: N 82.8.». 332: 55 «3 c c2. 2:. a San :.m o 8.“- 85 c 3.8 m3 3. 2.8 a... on :.3 NS 3 8.2 S... on 8.8 as a: 3.3 8.. 9: 8.2 02 c2 «32 a: 8. «22 2a =2 8.8 8a 8N 8.2: a... 8a 5.8 ”3 SN ”2.2 :.n 8,“ 3.8” 2; 8.“ 3.8” 8.“ 8m 8.8m 8a 8m 3.5 3... 8.. 2.8.. 8... as. 2.8. S.“ 8.. 2:: 8+ 2: 2.8“ 8.“ 8w 3.8“ S.“ can ~38 n3 8n 2.8“ 3.“ 8“ 3.8 8.» as 3.3“ 8.“ o8 $.93 .3 as S8“ 86 c8 ... 3.?» .— ... 3...; .— lwa... ...-F» .— WE 33$ .— 3.280 432.8 ~32mo .228 8:550 196 ~o~.- e. . .- 3.6.2 6666 ..«~... -~.. ~ 63.2. 366 .86 6n~.~. ~ n~66 an. 6 53.8.. 686 3. 6 v~n6¢~ -66 .-6 a... .... 366 ~¢~6 mam. ~. 266 66~6 86.2% -66 2~6 .862 266 ~o~6 2.86.. ~86 ~36 08.. .. 62.6 2.. 6 -~.3~ -66 an. 6 62.. . n . 266 66. 6 6. ~66 ~86 we. 6 v3.2 ~36 ~v. 6 -6.2~ e~66 ..~6 h. ~.~6 ~36 n... 6 2...? 266 ....6 .2..“ 3.66 2. 6 .266. 266 8~6 3.3m 3.66 6... 6 ~o~6~ .86 8. 6 n3. . e 2.66 ~6~6 ~86... 636 n. «6 63.2. 9.66 o-6 S. .n~ ~86 66~6 ~$6~ .86 «.26 6. ~.-. ~36 «86 8...? 366 .26 26.2 $66 . .~6 63.2 A3...... 866 ~36: ~86 ~63 ~86~ ~86 mi... ~86~ .86 ~32 .~.... ~ $66 «.36 @862 ~86 ~n... ~36~ 2.66 636 .~n6~ .~. 6 .c~.~ -¢6~ 866 n3..~ 6~6.-. 686 §.~ w~c6~ ~86 3.... 63.... an. 6 6~6.~ ~26. e~. 6 $~.~ .38.. ~v. 6 .8. ~36. o. . 6 8~.~ 68.: 66. 6 266 «~02 60. 6 .56.. ~86... .2 6 2~.~ ~26~ 8. 6 .36 63.... ~.~6 36.. «26. nn. 6 $66 6366 3~6 3N. ~86. ~26 2... n.~.~. -~6 6666 326. 8.6 «26 8:6 -~6 ~36 366. a... 6 6~66 62.2 ~w~6 68.. 6-.n. ~.~6 3. .. ~.«.? n-6 686 news. ~6~6 6.36 n.~.~ . ~26 86.. 6.3.2 -~6 6.6.. 2...? 3&6 3. .. ~36. ~.~6 ~86 6....~. ~26 866 «2.... ~v~6 we. .. o~. 65 no~6 266 wens. -~6 .26 6... .~. cw~6 n8. n2... 8~6 ~66. 68.2 6-6 36.. 69.6. 3~6 no. .. v~. ... .26 2.66 .2.~. 226 26.. «~98 n6~6 .2... 2nd. ~c~6 366 ..~. ... e. ~6 2... ~o~... 66~6 .26 ~36... ..~6 36.. ”.~.... 526 ~86 2.66. -~6 2.8.. 8~.~. 3~6 F. .. 626.. e-6 68.. 8n.~. ~o~6 we... .~.... 626 «.36 .2 .2 ~26 ~36 v.~6~ nn~6 ~36 66m.~. 8~6 63.6 636. 626 Sn... .2.~. -~6 3w... 3.6.? .c~6 626 3...» n .~6 6~n6 «~66 ~o~6 866 .2.~. -~6 n36 2~6~ 62 6 2...... ~26 -~6 626 «86 ~o~6 266 ec~.~. ~¢~6 .26 8. .v. 226 ~36 ~26 6-6 2n... 2.6... 3.6 .86 3~.~. 6-6 236 02... n-6 636 .96 e-6 63.6 8..... .26 ~36 ~26 2~6 ¢-6 6-.~ ~26 2n... ~36 ~¢~6 e~n6 ~56 626 .-.v no~6 ~26. 626 n25 6n~6 9.. u: 3.. 9.. us. 3.. a. u: 3.. .... u: 3.. ._ 2328 2.8.8 $28 .328 197 n~.~—~ 3.2. 5.. .— 8.2. hndn c0; 2.6m 3.2 3.2 ~13 2.». 2.5. and. 26. and. 2.2 2.... and. 2.2 cm... 2.2 2.2 2.~. 2.2 2d. 2.2 nmd ...: ... c. .c ...c ~— .c ~. .c m. .c v. .c n. .c w. .c -.c 2.c 3c .nc vac and ~c.c sec .56 Sc and mac and ccc Ncc vac vac «cc 3.35 cad: $.ch vc.2~ «~66. cod... 2.3. 3c: ”can. ~c.-. was: ccdc. 3.: 2.3 .92 3.2 v. .05 cc. 2. ~98 chow 2.3. 2.2 692 2.3 cmé. was Nwd hcc ncc ccc cc.c c. .c N. .c e. .c .~.c “~.c n~.c 9c :6 ac and ~c.c Sc ch.c Xuc 2.6 mac and c6.c 36 36 36 ac 2.3.. 2.3 cwdn 3.2 cwdv 8.2 3.2 ~.ch on... 2.c~ end. nnd— n92 end. unfi— 36. no.2 ~03 .m.~. cud. 8.» 2.6 Fwd 3..» San n. m n. m m.828 . . .328, 198 E ....)a . ...-«5...... p. ...o of -.c . Afiv . AV r. Aflv NNQ. ..u :06...— “3oz Ncnhc. me. .Q. . cc.¢~ .m.n .x. and 592% cvdhn «0.3 .v bthhc. h. «.3 . cc: $.N .x. 8.... on .hhn 9...th "v.3 m . 3360.. 25.3. n. .QN and .x. and- cod? 3.an 8.3 N _ 2 «650.. can... . mvém «QM .X 3am. *6.an vmdhn ...NN . Ac.v A2V A2v Agv An.v Av.v A2v AN.V _ all. =8 ...-3 5.3 all... «583 3.2:! 33> oz _ ...-n ... 23.3 g ...-:8: ...: In autumn-lam. .5810 m Ear—$36 § 02:— £3 >30 ZH>O 448 ...—.53:— ...(m :55 Brad—u $3. $55 =5 Bé 52m :— mo 199 Fe... . .. _ .— cwnc cc..c. ccc.c cc..c 3nd con.» 2.6..» and 36.» 3nd 62.3 366 cc..c ccw.N. 66o... cmhc ncac. Nnhc. c6n.n 63.2 eth. cnhc ”3.2 8. .2 c6N.c 66c... 6cm. .. c666 Nnhc. anc. chNd 656.2 n66... ccnc 6c. .2 aNhN. cNN.c 6%... .oc. .. c3..c aNnc. 62 .c. cu. .m A~66... ccc.c. cnmc cNhN. 65... cm. .c .cc... gnc. cc6.c 62c. .86 c... .m ccc.c. moNc. chc 6...”... c6. ... c6. .c 6wnc. as. .c. cnmc .cwd .36 c2 .n m6N.c. cm”... cth c3... .~.....c. cm . .c w... .c. ”mud cth .6nd mm. .c c2 .n cmwd «de ch.c n.h.c. wand. c2 .c ”mud n66.» c6N.c mm. .c 26.x cN. .n wde mncd co. .c ”mac. 2nd cN. .c «2.2 $5.” c6N.c NcN.N. ch; cccd nnc... who.” cw. .c hNh.2 ccc.c cc. .c Kc... SN... 8. .c 2.66. c. c... cwcd N2 .c. Q36. ”2 .c 366. Mac... cc. .c chN. 52 6. c2 .c c5... 6am... chcfi «.86. 63.2 ”N. .c 36.5. N36. ccc.c 526. N2... cc..c 63...”. 3nd. c6c.n 63.2 .6... N26 N36. 62.... ccc.c Nn. ... 66N.c. ccc.c find. 66»... cmcd .6. ... 63.5 cN. .c 62.... 36c. .wcc 3N6. New.» ccc.c 66nd ac...» c6c.m 66nd «66.» cc. .c 366. nhNd 65cc 2%.... .co.« chcc 62.6 and cvcd .mbd wNm.n c. . .c Ncwd ”on... 65cc .c6.n con... c...c.c «3.6 Nnmd cvcd ”N00 .... m6 c. . .c ”cm“ 65.6 chcc ccmd mm. .n cbcc Nmmd 6. ..n cvcd ... n6 N86 c. ..c 66nd 62 .m chcc mm. .n th6 chcc 6. ..n 62.6 c6c.n chd 2.56 c. . .c 62 .n 8&6 chcc N36 58 .n nmcc «3.6 8nd cmcd 6666 6cm.n cc. .c cwwd Nhnd choc 56nd 66nd nncc mend cmmw cmcd 8nd mcmd cc. .c 66nd wmmd c6c.c Span 6N. .n mncc cnmfi n. . .n cmc.n ncm6 vhci. ccc.c wmmfi wccfi. c6c.c 6N. .n 5Nm6 mncc n. ..n 6cm... cmcd 686 m6...6 ccc.c wood 686 6nc.c 62am nnnfi. ..6c.c 356 8nd ch.n and mth gcc 6656. ..Nnfi cncc mung». nnmd .6c.c cond 3nd ch.n and ch.n ccc.c hNnd N. m6 35c «and w. . .n omcc 6cm.n cccd c.c.m naNfi. Rcd nwcc N. N... cncd 66cc a. . .n ch6 nmcc cocd ”656 ccc.n :.cf. Nnh6 ccc.c cncd $56 ccc.c ch.c c6c6 ch.c.. ..o> 9.0.3 €3.93 ..o> ...—9.3 986.3 =0> 98.3 :83.» ..o> 98.3 5.3 66.3.0 $.20 N*.Nm0 .....NEO 200 .88.. :88. :88 88.8 88.8... :.88 88...... 88.8 8... 8.... c 8... 8.... o 2... ....n a 8.. 8... c 8..: 8.. 8. 8.8. ...... 8. 8..: 8.8 8. :..: 8.. 8. 8.5 8... 8. 8.8 8... 8. 88. 8.9 8. 8.8 8.. 8. 8.8 8... 8 ...: 8... 8 8.8 :..n 8 :8 R... :. 8.3. 8... 8 .3.. 3.... 8 2:. 38 8 8..... 8... 8 8.8 8... 8 :..8 8... 8 8.8 8.8 8 3.8 8... 8 8.8 8... 8 8.8 8... 8 8.8 8.8 8 8.8 8... 8 8.8 2... 8 8.... :... 8 8.... :.n 8 .3.. 8... 8 8.... 8... o. 8... 8... ... 8.... 8.8 ... 8.... :... ... 8.8 a... a 8... 8... c 8.... ....n c 2.... 8... o 8..: 8.. 8. 8.8. 8.. 8. 8.8. 8.8 8. 8.8. 8.. 8. .. 33$ .— .«~.... 5...; .— me... 3......» H...— 3....> I... 8.8.0 8.8.0 8.88 8.8.8 1 .5889 201 ccc.c wemd $6.: ocod hemd EVAN ccc.c «and «3.: ccc.c 3nd 9: .2 cm. _ n 356 $36 3.: $86 mn— .c 3.3 356 «:6 2 do «86 and 3.2 Sad mend 3&0 Sod mNNd 3.3 «86 find 33m #36 «Rd mmdm ccc.c mend 3.3 3&6 and ~03 cncd has... «man «36 53.6 2.2 2.96 bond _a.nv 956 mmmd 3.3 39° 3.06 RMN n86 nnwd 2.: 89¢ 03.6 3.5” n36 8nd 3.3 :.cd mmvd mnNN tbd «Ned 3.3 Sod envd mafifl awed $26 3.2 Ebd 8nd v0.2 «wad 526 2.3 mocd mom; mafim v86 3nd 3.3 39: n2”.— 86— owed mend 3.2 N— _6 3nd Nada uhcd mafia an: 8— .o 593 8.2 2: .o 62 .v 3%: an _ d 62..— gd— #2 .c «chm 2k: 3: 6 Sn." 3.6 «2 .c mg.— vm.: ow— .5 9.— .m Nad— 536 I _ .n 36 3:6 mad an.” :2 6 35d 3.2 «26 nwoN Nad— m_~.¢ Sod aw.” and nwofl as 2N6 3.— .m 3 .2 03d mead Nag: hand 95.— mafi hand man; :..h ”had sac.— 36 Quad NOVA Nad— mhmd «3..— wad $56 «and mud 3N... «on; 3.» 2 m6 hmmd «ad and mead wad baud mnmd N56 Ind QMNd 2..« m— md N26 Nod «36 Bad wad 8nd :Nd N56 and *2 6 cm.” 2 m5 8 N6 mad 5:. w— «6 wad mcmd wand and cmmd wmmd em.” 22. _wmd wad 32. 3nd mod hand bend N56 and 2&6 QM.» vamd nmm. Ne.” and RN.— mod :2. «mm; and «36 ~36 3”.” ~35 mcwd Ne.» and wad 3.0 and 5nd mbd nmmd «36 XS and 9&6 mad :2. n— Nd mod «and :86 and 986 $26 ems ”end 536 -.w «and :md 86 wand 82. mmfi 3nd Sud 3.0 mnmd and Nag. ”and 32. boa 3nd 036 n— .n 3nd 9&6 3.0 on md cc~d 8H ~36 Rad bod mend «and :6 send and $6 anmd med mad 3:. 3nd 2 .9 send «3.: 3d _nmd Nemd «fin mend wand Ned ovmd «26 S .m and find end :2. and X“.— hcnd Naé nnmd 3.0. «and 2 .9 and u.— m .55 o.— m 2 oh m a u.— m 39 38.—SO 2:20 QANmO { 37:50 202 78.8 25 8.8 :5 55.8 85 :55 5.5 8.8 :5 8.8 :5 8.8 :5 8.8 25 _88 2 5 85.. 2 5 8.2 :5 8.8 2 5 _82 85 8.8 25 8.2 :5 8.8 25 55.: 85 8.8 25 55.2 85 8.8 85 8.2 85 8.8 25 55.2 85 8.2 85 _8.: 85 8.8 :5 8.: 85 55.: 85 "8.2 555 858 85 8.: 85 55.2 85 8.: 85 85. 85 8.2 85 $5 2.5 8.: 555 8.2 85 85 85 85 85 8.2 85 8.2 85 85 85 88 555 2.2 2.5 8.2 585 88 85 8.8 85 3.5 85 8.2 85 85 85 85 85 .85 85 85 85 85 85 85 85 _ 85 85 85 85 85 85 85 85 85 85 85 85 85 555 85 85 _ 85 85 85 85 85 85 85 85 _ .55 85 85 85 85 85 85 85 w .55 85 85 85 85 85 85 85 _ .55 85 85 85 85 85 85 85 .2 85 85 85 85 85 85 85 85 85 8.8 85 85 85 25 85 .85 85 88 85 S5 85 8... 85 _85 85 85 85 83 85 8... 85 r .85 85 85 85 :5 85 85 8.. 5 n. m m 8 a 8 m m 5 8.880 8.880 8.880 8.880 w _ . 203 £0 0.... -C . A~_v.A¢v u Av 35.. n :05...— ”302 02.3. and." tan .3 3.". 3.2.... 385. 2.5" c 82.0 _ N3 5.: 3.3 end o\e 3d. Pd? and? 5.3 m 32.3 _ «and: 35H 38 o\.. SN. 8.»? «~63. 3:: N 082.0 ~ «3...: n _ .3 8H.” .8 and- 2.0% 3.0% 3.3 _ Aav sz AV A2v 2v #7 AV AV __ ”man: :8 5.3 :35 8.3!. 8.3;. 3.14 5.3 02 313.3 3.3.: 1.3.3. :..: In 3:130: 38an «an 8 | 85:52., - Egg. - .. - HG - 0_ 8.2. 3.62 wadNN 8d ac —©m mafion madam and?“ 3d: no.2; Q 5_ vwdh hbdo— andmN awn 8A3” hwéan Egan sbfig nndun 2 .NOQ m “ ands. :..:— NQdNN 36 mmdwm endan seq—an ”mach _cdz —N.NOQ N H :.mh «n4: aaénn cad mmNcm ~N6an madam 3.59% m— .w—w Nuance — h A _ v A: v on va Ahv Acv Anv Avv Amv Auv A—v _ 05: d 02:. a. 5.828 a do .=0 .58 _- 95. 5.05 ...8 Eu uzi e: 50.. e: .52.. ...? maps; 5:3 :86: E25 _ EB zm>0 ...—z. 323:. ...—no 3.5.138 Oz W 2:. guy—Hg HARD 35m (348° 26m 3.0 .528 204 ”urn . hilfiflumlr’. 000 .0 0m00— 0000— 000.0 5050— N—00— 000.0 0000— 0N0.0— 000.0 000.2 0N0.N_ 00N.0 0000— 000.0 00N0 N00- 050.0 00N.0 0N0.0— 0—00 0— m0 0N0.N~ 5.0.: 0MNO 000.0 0 ~00 00— .0 050.0 00n.0 00N.0 00.0 500.0 00N0 53.: 20.0— 0NN.0 0— 0.0 N000 0— u .0 .00m .0 :00 0mN.0 50~ .0 000.0 0MN.0 n— 0.0— 000.0 00N0 0000— N0— 0— 00— .0 _0_ .m— 5N0.N~ 00— .0 500.0— 030— 0NN.0 0500“ m00.0_ 00— .0 N0_ 0— 000.0— 00— .0 5N0.N_ 00— .N— 00— .0 0~N.0— mm0.N— 00N.0 000.0— N000— 00— .0 000.0— 0N0.0 000.0 00— .N— 005.0 0N— .0 n00.N— 500.0 00- .0 N000— 005.0 5— _ .0 0N— .0— 050.: 000.0 000.0— _5— .0— 00— .0 000.0— 005.2 00- .0 NMO0N 0500— 0500 05m.0— 000.: N000 _5_ .0— N000— 000.0 0050— 0000— 0: .0 0500— 000N— 0000 000.: N000 000.0 N000— 00— .0 050.0 0000— 000.0 02 .0 000N— N500 0000 500 000.0 N000 000.0 000.0 500.0 0N0.0 ~000 0— _ .0 550 5—00 0000 00m0 000.0 000.0 0000 N00 000.0 800 000.0 00— .0 5— 0.0 000.0 500.0 500 N000 000.0 000.0 N—50 000.0 00.0 N000 00— .0 005.0 NN00 0000 N000 0500 000.0 N~50 ~000 N000 N000 0000 00— .0 NN00 _Nm0 M000 05m0 00— 0 000.0 _000 0MNO 000.0 0000 N: .0 000.0 :00 000.0 N000 00- 0 000.0 N000 0nN0 02.0 000.0 N5— .0 N000 N000 000.0 555.0 000.0 5050 :50 «00.0 000.0 _N00 000.0 0000 0N50 000.0 0050 0000 000.0 :50 NN00 N000 _N00 N050 000.0 0N50 m000 000.0 0m00 000.0 000.0 NN00 0000 000.0 N050 5000 000.0 0000 _000 000.0 000.0 _000 500.0 000 5000 000.0 5000 0000 000.0 _000 0500 000.0 00.0 0000 000.0 5000 _000 000.0 90.0 N000 000.0 0500 0500 000.0 0000 00m0 m000 _0m0 N5N0 N000 N000 50m0 000.0 05m0 00N0 000.0 000.0 _0— .0 000.0 N5N0 55— 0 000.0 5000 00N0 000.0 00N0 0~N0 050.0 —0— .0 0000 :00 55— .0 N000 0m00 00N0 0500 000.0 0—N0 0N~ .0 050.0 000.0 000.0 0N0.0 N000 050.0 NN0.0 050.0 _N00 000.0 0N— .0 0000 000.0 000.0 000.0 0N0.0 050.0 0N0.0 5—00 _N00 000.0 08.0 0000 0N0.0 000.0 000.0 0N5.0 0000.. 08.0 0N0.0 0000 .r as .393 5a-; as 29.3 as; j; a; as 293 5a-; 00.2.00 Qua—00 NaTNNaO :TNNBO 205 3.88: 2.2.2 .25? 2.2: 5E» é»... in? .3 85 c on.“ 2:. o 8.” 8.» o and 8... o _ :..2 2... cm 2.2. 25 on 2.2. 26 on 8.2. and 8 n 8.8. 8.. :2 3.8. :2 9: 852 8s 2: :62 8.. 8. m .3: :2 O: ”22 2.. an «5.2 as ca 5.2; 2.. .5 W 3.8. 8.. .2 3.2 8.. .2 $2 :3 82 :2 2: c2 _ 32. 26 a 3.2. MS a 2.8 8.2 a 2.: 8.: a w :..2 N3 2 2.2. 28 on 2.3 $6 an 8.2. a... R w 2.2. 25 2. 3.2 25 2. 3.2 2.... 8 8.2. 28 2. F 2.2 :.o 2 3.2 25 2 5.2 one 2 8.2 23 on . . 2.2 :... 2 :..2 c3 2.. 3.2 22 2 22 3o 2 “ ~ ..2 :.c 2 2.: 2.: 2 3.2 2.0 2 85 25 2 _ . .3 86 c on.” 85 a mi . 86 c :..? 25 o a ... 53$ .— LWE 36$ .— -.. 3......» .— Luoflm 8.....» m 3328 2.2.8 «328 .328 53:35.9 206 555.5 52.5 ~m~..~ .~5..~ 5%.: 25.2 ~.m.- 555.5 35 25.2 2.55 525 25.2 255 25.5 52.55 .555 25.5 ~5~..~ 255 92.5 ~22 :55 ~25 52.2 2.55 2~5 -55~ :55 555.. 52.2 ~555 225 52.5. 255 5555 ~555~ 95.5 255 25.2 555.5 5555 52.2 .85 52.5 .25 555.5 «25 ~22 255 2.... 2...: .85 5555 52.... .85 52.5 .25 555.5 525 555.2 555.5 ~25 25.... «~55 :55 5-.~. 5~55 555.5 .25 -55 .35 25.5. 5~55 .2.. 2...... 5555 .2.. -~.5. 25.5 555.5 52.5 3.55 .~...n 25.2 5555 58... 555.2 525 ~5~... ~35 ~25 2: ~25 2.5 52.5 25.5. 2.5 5555 555.2 55.5 23. 25.5 ~25 5.5... 52.5 2~5 52... 3.55 2~5 52... 5555 2~5 .2.~ ._ «~55 .25 2.3 52.5 $~5 ~5~.~ ~53 5-5 ~5..~ 52.5 55~5 5.5.~ 25... 525 2.5.5 525 52.5 -55 2.3 525 52.5 ~25 .25 25.5 __ 25... 52.5 255 552 .25 52.5 ~5~.~ 52.5 255 525 525 25.5 25.... ~25 «25 55.5 22.5 ~25 2.2 ~25 :3 52.5 52.5 52.5 -~.m 525 5.~5 55.5 ~25 .25 2...~ «25 83 52.“ 525 .25 .22 525 23 ~55... n25 -~5 ~25 52.5 23 52... ~25 -~5 52% ~25 2~5 ~52. 525 555.5 25.5 ~25 52.5 52.2 52.5 -~5 5-.m 52.5 5555 555... ~25 3.5 25.. ~25 .25 22.2 525 .25 5-.n ~25 255 ~55... 225 25.5 25.. 52.5 58.5 ~5~.~ 525 25.5 2 .... 52.5 3.55 ~53 “2.5 2.5.5 252 52.5 25.5 ~5~.~ ~25 555.5 fl 25.~ 525 .555 553 525 25.5 25.5 .25 25.5 ~5~.~ ..25 ~25 2~.~ .25 5555 35.2 525 .555 255 ~25 25.5 ~5~.~ 52.5 25.5 ~.n.~ ~25 25.5 552 52.5 55.5 555... 22.5 3.55 ~5~.~ 525 2.5 25.~ 2.5 355 $2 225 255 25... 32.5 .555 52.. ~25 ~85 2.... 525 2.5.5 8.2 52.5. 52.5 555... 52.5 25.5 52.. 525 52.5 5555 52.5 2.5 2.... 225 555.5 25... ~25 25.5 555.. 52.5 5555 5.5.5 525 ~25 25.. 525 2.5 25.2 525 .25 52.. 525 «2.5 .~...? 525 52.5 52.5 ~55.~ 525 2.5 525 5.. Us 3.. 5.. us. 3.. 5.. 02 3.. 5.. 02 3.. 5.5.2.8 2228 $250 .2280 207 no 9 and hm n end was. on... 86 3.6 No v and hm m on... 3.5 and 86 had No v a... N n S... as... 86 Sam and No m mad 5 n cad nmfi Na... and 3... NE.“ 36 h. n cod 2.5 86 ON.n Nad Nn m 36 .5 v a... nmé cad aNé 8... N.. a mad hm v Nod 36 mad am...” .36 NN n 3... he v cad nafi cad «Na 36 NN n had hr. ... a... now. had EN 86 N. n had hm v 3... mad had EN 36 Na N wad 5N e 3... mod had 3N had NEN mod n» m had 3.... mad EN 8... Nn N wad on...“ no... cod. 36 EN mad 3 N am... cN m mad ad wad 2... nod :.. 86 on N wad 02¢ cad ab.— 36 _c c a... on . 3... 9:. 92¢ a9. a... .m a 8.. co . 3... ea...“ 36 an; a... n. w n. w n. m n— w XVNN FwO nfiuNN FmO NtrNN wwO wt.“ ..wO 208 29:3 9.3: M33 me .x. 2: 3.: win w 28?: 8:: 3m 3m .x. as 8.2 3.2 F. SES 8?: 2:3 mm... .x. we.” 8.: no? N «233 83: 5.2 3..“ .x. 35 8.2 2.: _ AV G? A? A? AV 3!? __ .5 85:5 .5 ..c a: he :— 63 coau3n< 328m n32 copay—BU 12:55 nu. ma a a: u Av a a: u A3 a 2 _ a u Av ” 5:28 8m 85 802 flfl“ 2.? 8.8. -62 3d and? 8.30 and—n fiance w. 8.? 38— 26.2 cad Inn? 2.93 3.2m «33. m 2.2. 3%: 932 x; :..avn Swan 3; 3. $63 n 3.2. $62 2.2: and 86?. 2.3m 3.2m. 2.8» _ on va Abv on Anv Avv Amv Amv ANV A7 62:. d a a manilfi a. 5:5. Uzi 4.9. 58 E2» ozi QB .52.. 4,2 -35 :86: E2» b:— zm>c ...—z. @235 4;. Oz En. Emu—H3 .~th «SA BEE az 20.3 as.-.» ....> 9.0.3 5...-..» ....> 2.0.3 .53.. ....> 9.93 52...; 4.2.6.0 3&6... . .~....ng :..an . 210 AWN—Zn.— 35%.: wcfi # m~.N . . . cm.— mwd . . an. :2. and . . on; . end ...... . . ...... ...... 3...; 3......» 3...; :.....umfimo - l , - sahWMwwmoi- - ,- ...---rai -~.fimmmmp .- -i L - , I -.¥m~..mo 8.553. 211 «on... «8.3” c ”an... www.mm o 3nd .Nnén a own... 26.? 02.6.: wee... emu... 3o. 3. Rod .8... «no.3 .3... .3.: e562 «3... owed «3...: one... and 3.. .2. Rad n. _ .o 3&3 .36 3...: 52.5 «3... can... 3...? 8...: . .«.o 35.2. ”8... mm. .o «3...: ~35 ”um... and.” Rs... 3:. Swan 3o... w . n... 98. _ n a8... «85 No.93 one... :2 omm.m~ Sad and 2.3.5. «3... m2... .2.: a8... «8... ska: n8... one... 2 n.- «8... _Rd 3...? who... «m: 23H ccc.c mead 2o.“— Ngd «R.— 32: «he... 2 _ .— 3. .8 3...... 3a.” wad. need 2.... Sa... 086 m2 .n can... «8... «a. .m .2...... E _ .c ”an... 28.2 .2.... .2.... «3.: 2 _ .o .2.? 8N2 .2... can... .32... an. ... 2:. 98.. mm. .c «an... 3: .o en. .c 8m... find «2 .o 2.. ... Exam mo. ... ans... «.36 3. .o «2..«. cmms a. .o «$6 wond— 3. .c 25. Sad.” 2."... 3m... 93.... 82. man... ”and 22. ~33. 3nd 3N... Nun... .«.....NN EN... 95.. 3a.... .3... «hr... ”and n5... at: 02.5 mum... Rn.— ..vd. can... 3: cu. .m «mm... «5.— vuné am... «3.. 93.5 am... an.— uvnd. c. m... can ._ 03..» 2k... we»... and and 2.... «mud and «no.— Saé 3m... cm: on...“ ..«~... .~.... an.“ in... so“... «36 En... how.— 8. .2 c2... .3.. 8:. and 3m... and «mm... 2.2 53.... «mm... n3... .2...». 32. 03... 3.2“ 2 m6 3. .c 30.” 3nd 3N... 3m... 3m... 3.. .o at“... man... 3nd oafim n. m... 3. .c one." and ”mm... men... 3nd 86c find «2.: $2. sand 2 m... man... mag Nam... 2h... 2.”... «end mend «3d Rm... 2."... afiu a. m... «on... 36." 3m... «3... me”... can... 3nd RNN 8:. SN... and «am... .2... Rs.— onmd aha... new... $2. new... 8m... Sn... 3. .o 25.. wumd 3 ad «and own... «a. .c an?“ 9.2. mm. .c «on... mom... 3. .o 25.— oumd ban... was... Sn... 3.”... 52.6. Rm... .5... m9... 7 3m... and. an... R. .o. mom... 8:- awn... o.— 02 3.. on us. 3.. on US 3.. a.— US 3.. 1228 2.8.8 972.8 $328 212 no.8? n ...o ... ..o. ...o nndw v. ... undo. I ... no. ...... m ...o ... .wh .... 26.. v. ... wmfin n. ... 9&5 ofio :..... N. ... 2mg 9. ... nm...m c. ... Eon h v... 8.: a. ... aodn m. ... Nndn h. ... soda 2.0 and... N. ... and. a. ... NM.NN w. ... .nm ..Nd 2.3 u. ... .ad. cm... and. ..N... Ndm mud Na... and ma... «m... cm... mm... hvév No.0 ... .... nu... no... NM... wad. .m... adv nvd 3... mm... a. 6 9.... an... 3... 9mm 3... N”... on... «as mm... ~.«.... cm... had» cod .....o ...... 36 no... 36 8... aaNu NB... 8.6 ms... and 3... cm... 55... Stop Nu... m. .m an... and cm... enfi ...... m6. mod .....w an... an... 3... 0N6 3... ad. and ... .n ..w... an... a... nné am... for mud on... .a... an... ma... mad 3... 9m had en... ca... cod 8... a... 3... 0.1 had as." ea... 8d 3... a... a... .n cad «EN 3... 8N no... 36 ma... .0 cad NEN a... cod 3... en... a... NN 8.0 «ha ma... 2.. ”a... a... g... cm... wad 8.. 3... cm... 3... a... 8.. on... co... m... ...... cw... 3... as... ...... n— w n. a n. #0qu 0%.”me NXNwO _inanO