The S shaped log log drawdown time curve typical of pumping tests in unconfined aquifers is reinvestigated via numerical experiments. We examine the temporal and spatial evolution of the rate of change in storage in an unconfined aquifer during pumping. This evolution is related to the transition of water release mechanisms from the expansion of water and compaction of the porous medium to the drainage of water from the unsaturated zone above the initial water table and initially saturated pores as the water table falls during the pumping of the aquifer. We conclude that the transition of the water release mechanisms and vertical flow are the cause of the S shaped drawdown time. Cross-correlation analysis is then employed to examine the relationship between hydraulic properties of an unconfined aquifer and pressure observations. The analysis reveals that head observed in the saturated zone at late times along a streamline is positively correlated with the conductivity (K(s)) of the region upstream of the observation location, and negatively correlated with the K(s) of the region downstream of the observation location along the same streamline. Besides, head observations in the saturated zone at the early time are positively correlated with specific storage (S(s)) in a narrow region between the observation and pumping locations. At intermediate and late times, the head positively correlates with the heterogeneity of α (pore-size distribution parameter) in a thin disk-shaped unsaturated region above the pumping and observation locations. Saturated water content θ(s) in the vadose zone directly above the pumping and monitoring locations is found positively correlated with the head observations during the intermediate times and late times.In the end, a stochastic inverse estimation is conducted to jointly interpret a sequential pumping test in a three dimensional unconfined aquifer. K(s), S(s), θ(s) and α are estimated at the same time. The estimated results capture the pattern of the heterogeneous parameters as well as the details with a smooth distribution. The estimated heterogeneous parameter fields produce better head predictions than the traditional homogeneous method.

Taludes formados por solos na condição não saturada são muito comuns no Brasil. Nestes casos a poro-pressão de água contribui para a estabilidade desses taludes. Este trabalho pretende avaliar o efeito do tipo de revestimento superficial na alteração do perfil de poro-pressão de água provocado por fluxos gerados a partir de precipitações, evaporação e evapo-transpiração, e conseqüentemente na estabilidade de talude experimental. São apresentadas aqui as características geológico-geotécnica, geomorfológica, climática e pluviométrica de campo experimental implantado na zona leste de São Paulo/SP, onde foram instalados e monitorados tensiômetros, pluviômetro e medidor de nível dágua por um período oito meses. Resultados de ensaios de caracterização, resistência e obtenção das propriedades hidráulicas dos solos que compõem o maciço estudado são analisados, assim como ensaios realizados com a argamassa de revestimento do talude, ensaios para a análise do funcionamento dos tensímetros e tubos tensiométricos utilizados no monitoramento, e outros ensaios realizados in situ. O monitoramento realizado indicou coerência entre precipitações e variação da poropressão de água, e variação do nível dágua do terreno, em ambos os tipos de revestimento superficial, vegetação natural e argamassa, nos dois tipos de solo presentes no maciço. Análises de estabilidade do talude indicaram o efeito do tipo de revestimento superficial e das variações sazonais de poro-pressão de água no fator de segurança ao escorregamento do talude. / Slopes formed by unsaturated soils are very common in Brazil. In these cases porewater pressure contribute for the slope stability behavior. This work aims at evaluating the effects of types of superficial coating on the profile of pore-water pressure changes, due to flow generated by precipitation, evaporation and evapotranspiration, and consequently on the stability of a experimental slope. It is presented here the geological-geotechnical, geomorfological, climatic and pluvial characterizations of the experimental field located in the east zone of São Paulo/SP, where have been installed and monitored tensiometers, pluviometer and groundwater level measurer during eight months. Results of tests of characterization, shear strength, and for obtaining the hydraulic properties of the in situ soils are analyzed, as well as the tests carried out with de mortar applied over the slope, as well as the functionality tests for the tensimeters and tensiometric tubes used on the monitoring, as well as other in situ tests. Monitoring indicated coherence between precipitation and change in pore-water pressure profile, and groundwater level, in both types of superficial coatings, natural vegetation and mortar, and in both types of in situ soils. Stability analyzes indicated the effect of the type of superficial coating and of the variation of pore-water pressure profile, during the experiment, on the factor of safety of the experimental slope.

Taludes formados por solos na condição não saturada são muito comuns no Brasil. Nestes casos a poro-pressão de água contribui para a estabilidade desses taludes. Este trabalho pretende avaliar o efeito do tipo de revestimento superficial na alteração do perfil de poro-pressão de água provocado por fluxos gerados a partir de precipitações, evaporação e evapo-transpiração, e conseqüentemente na estabilidade de talude experimental. São apresentadas aqui as características geológico-geotécnica, geomorfológica, climática e pluviométrica de campo experimental implantado na zona leste de São Paulo/SP, onde foram instalados e monitorados tensiômetros, pluviômetro e medidor de nível dágua por um período oito meses. Resultados de ensaios de caracterização, resistência e obtenção das propriedades hidráulicas dos solos que compõem o maciço estudado são analisados, assim como ensaios realizados com a argamassa de revestimento do talude, ensaios para a análise do funcionamento dos tensímetros e tubos tensiométricos utilizados no monitoramento, e outros ensaios realizados in situ. O monitoramento realizado indicou coerência entre precipitações e variação da poropressão de água, e variação do nível dágua do terreno, em ambos os tipos de revestimento superficial, vegetação natural e argamassa, nos dois tipos de solo presentes no maciço. Análises de estabilidade do talude indicaram o efeito do tipo de revestimento superficial e das variações sazonais de poro-pressão de água no fator de segurança ao escorregamento do talude. / Slopes formed by unsaturated soils are very common in Brazil. In these cases porewater pressure contribute for the slope stability behavior. This work aims at evaluating the effects of types of superficial coating on the profile of pore-water pressure changes, due to flow generated by precipitation, evaporation and evapotranspiration, and consequently on the stability of a experimental slope. It is presented here the geological-geotechnical, geomorfological, climatic and pluvial characterizations of the experimental field located in the east zone of São Paulo/SP, where have been installed and monitored tensiometers, pluviometer and groundwater level measurer during eight months. Results of tests of characterization, shear strength, and for obtaining the hydraulic properties of the in situ soils are analyzed, as well as the tests carried out with de mortar applied over the slope, as well as the functionality tests for the tensimeters and tensiometric tubes used on the monitoring, as well as other in situ tests. Monitoring indicated coherence between precipitation and change in pore-water pressure profile, and groundwater level, in both types of superficial coatings, natural vegetation and mortar, and in both types of in situ soils. Stability analyzes indicated the effect of the type of superficial coating and of the variation of pore-water pressure profile, during the experiment, on the factor of safety of the experimental slope.

Cette thèse vise à étudier expérimentalement, analytiquement et numériquement, les conséquences de variations et d'oscillations hydrodynamiques à forte variabilité temporelle en milieux poreux partiellement saturés. Les problèmes que nous étudions comportent des surfaces libres tant à l'extérieur qu'à l'intérieur des milieux poreux, celles-ci étant définies comme des isosurfaces de pression d'eau égale à la pression atmosphérique (Pwater = Patm). Les différentes études expérimentales réalisées en laboratoire sont, respectivement : une expérience d'imbibition dans une boite à sable avec effets capillaires importants; la transmission d'oscillations de la surface libre à travers un massif sableux intercalaire dans un petit canal à houle (IMFT, Toulouse); l'étude de la dynamique et de la propagation des oscillations des niveaux d'eau dans un grand canal à houle (HYDRALAB, Barcelone), partiellement recouvert d'un fond sableux incliné, avec mesures de niveaux d'eau en pleine eau et sous le sable, et mesures du fond sableux (érosion/dépôts). Pour les études théoriques, nous avons développés des solutions analytiques linéarisées. Un exemple de problème traité analytiquement est: l'équation linéarisée de Dupuit-Boussinesq (D-B) transitoire à surface libre, en hypothèse d'écoulements plans et vidange/remplissage instantané : oscillations forcées, transmission et dissipation d'ondes à travers une boite à sable rectangulaire. Nous avons aussi développé une solution de l'équation faiblement non linéaire de Dupuit- Boussinesq (D-B) pour étudier le problème d'imbibition avec variation abrupte du niveau d'eau amont (suivi temporel du front de saturation). Nous avons pu étudier les différents types de problèmes transitoires liés aux expériences citées plus haut par simulation numérique. En particulier, nous avons simulé des écoulements partiellement saturés et insaturés, en coupe verticale, à l'aide d'un code de calcul (BIGFLOW 3D) qui résoud l'équation de Richards généralisée en régime transitoire. Nous avons ainsi étudié numériquement en régime non saturé, l'expérience d'imbibition dans un sable initialement sec à frontières verticales (IMFT sandbox), puis l'expérience de propagation d'ondes dans le grand canal à houle de Barcelone (laboratoire HYDRALAB) comportant une plage de sable inclinée, avec un couplage complètement intégré entre les zones micro-poreuse (sable) et “macro-poreuse” (pleine eau). Pour analyser les résultats de cette dernière expérience et les comparer aux simulations, nous avons utilisé plusieurs méthodes de traitement et d'analyse des signaux : analyse de Fourier (spectres de fréquences) ; ondelettes discrètes multi-résolution (Daubechies) ; analyses corrélatoires simple et croisée. Ces méthodes sont combinées avec des méthodes de préfiltrage pour estimer dérives et résidus (moyennes mobiles ; ondelettes multi-résolution). Cette analyse des signaux a permis de comprendre et quantifier la propagation à travers une plage de sable. Au total, les différentes approches de modélisation mis en oeuvre, associé à des procédures de calage en situation de couplage transitoire non linéaire ont permis de reproduire globalement les phénomènes de propagation de teneur en eau et de niveau d'eau dans les différentes configurations étudiées. / This thesis aims at investigating experimentally, analytically and numerically, the consequences of hydrodynamic variations and oscillations with high temporal variability in partially saturated porous media. The problems investigated in this work involve “free surfaces” both outside and inside the porous media, the free surface being defined as the “atmospheric” water pressure isosurface (Pwater = Patm). The laboratory experiments studied in this work are, respectively: Lateral imbibition in a dry sand box with significant capillary effects; Transmission of oscillations of the free surface through a vertical sand box placed in a small wave canal (IMFT, Toulouse); Dynamics of free surface oscillations and wave propagation in a large wave canal (HYDRALAB, Barcelona), partially covered with sand, with measurements of both open water and groundwater levels, and of sand topography (erosion / deposition). For theoretical studies, we have developed linearized analytical solutions. Here is a sample problem that was treated analytically in this work: The linearized equation of Dupuit-Boussinesq (DB) for transient free surface flow, assuming horizontal flow and instantaneous wetting/drainage of the unsaturated zone: forced oscillations, wave transmission and dissipation through a rectangular sandbox. We also developed a weakly nonlinear solution of the Dupuit-Boussinesq equation to study the sudden imbibition (temporal monitoring of the wetting front). We have studied the different types of transient flow problems related to the experiments cited above by numerical simulation. In particular, we have simulated unsaturated or partially saturated transient flows in vertical cross-section, using a computer code (BIGFLOW 3D) which solves a generalized version of Richards’ equation. Thus, using the Richards / BIGFLOW 3D model, we have studied numerically the experiment of unsaturated imbibition in a dry sand (IMFT sandbox), and then, with the same model, we have also studied the partially saturated wave propagation experiment in the large Barcelona wave canal (HYDRALAB laboratory), focusing on the sloping sandy beach, with coupling between the micro-porous zone (sand) and the “macro-porous” zone (open water). To interpret the results of the latter experiment and compare them to simulations, we use several methods of signal analyzis and signal processing, such as: Fourier analysis, discrete multi-resolution wavelets (Daubechies), auto and cross-correlation functions. These methods are combined with pre-filtering methods to estimate trends and residuals (moving averages; discrete wavelet analyses). This signal analyzis has allowed us to interpret and quantify water propagation phenomena through a sandy beach. To sum up, different modeling approaches, combined with model calibration procedures, were applied to transient nonlinear coupled flow problems. These approaches have allowed us to reproduce globally the water content distributions and water level propagation in the different configurations studied in this work.