Liquefaction And Pore Water Pressure Generation In Sand - A Cyclic Strain Approach

This paper presents the results of laboratory investigation carried out on Ahmedabad sand on the liquefaction and pore water pressure generation during strain controled cyclic loading. Laboratory experiments were carried out on representative natural sand samples (base sand) collected from earthquake-affected area of Ahmedabad City of Gujarat State in India. A series of strain controled cyclic triaxial tests were carried out on isotropically compressed samples to study the influence of different parameters such as shear strain amplitude, initial effective confining pressure, relative density and percentage of non-plastic fines on the behavior of liquefaction and pore water pressure generation. It has been observed from the laboratory investigation that the potential for liquefaction of the sandy soils depends on the shear strain amplitude, initial relative density, initial effective confining pressure and non-plastic fines. In addition, an empirical relationship between pore pressure ratio and cycle ratio independent of the number of cycles of loading, relative density, confining pressure, amplitude of shear strain and non-plastic fines has been proposed.

Laboratory determination of coefficient of consolidation from pore water pressure measurement

In the recent past, many studies have been carried out on the determination of coefficient of consolidation (c(v)) from the time (t)-deformation (d) data obtained from conventional consolidation tests. Several researchers have also proposed different curve fitting procedures for determining cv from the t-d data. It is anticipated that the cv values obtained from the t-d data may be influenced by initial and secondary compressions. Nevertheless, the pore water pressure data measured during the consolidation process will be independent of initial and secondary compressions. In this study, the conventional Asaoka (1978) method is extended to evaluate cv and end-of-primary (EOP) consolidation from the pore water pressure data measured from laboratory experiments. Laboratory experiments were carried out on the modified one-dimensional consolidation apparatus on different remoulded clay samples measuring pore water pressure during the consolidation process. The cv and EOP computed from the proposed approach have been compared with the results of the t-d data and found to be in good agreement.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

The Liquefaction and Displacement of Highly Saturated Sand under Water Pressure Oscillation

In order to investigate the characteristics of water wave induced liquefaction in highly saturated sand in vertical direction, a one-dimensional model of highly saturated sand to water pressure oscillation is presented based oil the two-phase continuous media theory. The development of the effective stresses and the liquefaction thickness are analyzed. It is shown that water pressure oscillating loading affects liquefaction severely and the developing rate of liquefaction increases with the decreasing of the sand strength or the increasing of the loading strength. It is shown also that there is obvious phase lag in the sand Column. If the sand permeability is non-uniform, the pore pressure and the strain rise sharply at which the smallest permeability occurs. This solution may explain Why the fracture occurs in the sand column in some conditions.

Strength reduction of till under dynamic pore pressure condition

Cuttings in heavily overconsolidated clays are known to be susceptible to progressive deformation caused by creep and fatigue that usually begins at the toe of the slope. The progressive deformation leads to strength reduction with time at constant stress (or called softening) and could be accelerated by fluctuation of groundwater level associated with more extreme rainfall events predicted through climate change. The purpose of this paper is to assess the mechanism of progressive deformation due to creep and fatigue using element testing on samples of till. The samples were subjected to fully drained loading and the deviator stresses were held constant at various percentages of peak failure stress, while the pore water pressure was kept static or dynamic (fluctuating ±5 kPa) over a period of time. The results have shown that the samples experienced significant deformation even at a higher factor of safety (i.e. the failure deviator stress/deviator stress at which the pore water pressure was fluctuated) under pore water pressure dynamics.

Squeezing of Pore-water from Core Samples of DGR Boreholes: Feasibility Study

Three archived core samples from boreholes DGR-4, DGR-5 and DGR-6 from the Salina F Unit, Queenston Formation and the Georgian Bay Formation were subjected to squeezing tests at pressures of up to 500 MPa. Two samples did not yield any water, while a total of 0.88 g pore water was obtained from a clay-rich sample from the Blue Mountain Formation (water content = 2.8 wt.%, porosity = 8 %). This water mass was sufficient for a full chemical and water-isotope analysis – the first direct determination of pore-water composition in rocks from the DGR boreholes. The results are generally in reasonable agreement with those of independent methods, or the observed differences can be explained. Ancillary investigations included the determination of water content, densities and mineralogy, aqueous extraction of squeezed cores, and SEM investigations to characterise the microtexture of unsqueezed and squeezed rock materials. It is concluded that squeezing is a promising method of pore-water extraction and characterisation and is recommended as an alternative method for future studies. Selection criteria for potentially squeezable samples include high clay-mineral content (correlating in a high water content) and low carbonate content (low stiffness, limited cementation). Potential artefacts of the method, such as ion filtration or pressure solution, should be explored and quantified in future efforts.

Pore-water squeezing from indurated shales

High-pressure mechanical squeezing was applied to sample pore waters from a sequence of highly indurated and overconsolidated sedimentary rocks in a drillcore from a deep borehole in NE Switzerland. The rocks are generally rich in clay minerals (28–71 wt.%), with low water contents of 3.5–5.6 wt.%, resulting in extremely low hydraulic conductivities of 10− 14–10− 13 m/s. First pore-water samples could generally be taken at 200 MPa, and further aliquots were obtained at 300, 400 and 500 MPa. Chemical and isotopic compositions of squeezed waters evolve with increasing pressure. Decreasing concentrations of Cl−, Br−, Na+ and K+ are explained by ion filtration due to the collapse of the pore space during squeezing. Increasing concentrations of Ca2 + and Mg2 + are considered to be a consequence of pressure-dependent solubilities of carbonate minerals in combination with sorption/desorption reactions. The pressure dependence was studied by model calculations considering equilibrium with carbonate minerals and the exchanger population on clay surfaces, and the trends observed in the experiments could be confirmed. The compositions of the squeezed waters were compared with results of independent methods, such as aqueous extraction and in-situ sampling of ground and pore waters. On this basis, it is concluded that the chemical and isotopic composition of pore water squeezed at the lowest pressure of 200 MPa closely represents that of the in-situ pore water. The feasibility of sampling pore waters with water contents down to 3.5 wt.% and possibly less opens new perspectives for studies targeted at palaeo-hydrogeological investigations using pore-water compositions in aquitards as geochemical archives.

Pore fluid pressure at DSDP Leg 84 Holes

Evidence of considerable overpressuring of pore fluids in the sediment drilled during Leg 84 was obtained from direct measurement of pressure by two methods. The first involved measurement of back pressure when the annulus of the drill hole became constricted with unremoved drill cuttings or constriction was caused by plastic inflow of the drill hole walls. The second involved measurement of pressure ahead of the bit in conjunction with in situ water samples and heat flow. All measurements indicated abnormally high pore pressure even in slope deposits of the Middle America Trench off Guatemala.

Pore water acetate and hydrogen, helium and methane gas concentrations in ODP Leg 204 sediments

Acetate and hydrogen concentrations in pore fluids were measured in samples taken at seven sites from southern Hydrate Ridge (SHR) offshore Oregon, USA. Acetate concentrations ranged from 3.17 to 2515 µM. The maximum acetate concentrations occurred at Site 1251, which was drilled on a slope basin to the east of SHR at depths just above the bottom-simulating reflector (BSR) that marks the boundary of gas hydrate stability. Acetate maxima and localized high acetate concentrations occurred at the BSR at all sites and frequently corresponded with areas of gas hydrate accumulation, suggesting an empirical relationship. Acetate concentrations were typically at a minimum near the seafloor and above the sulfate/methane interface, where sulfate-reducing bacteria may consume acetate. Hydrogen concentrations in pressure core samples ranged from 16.45 to 1036 parts per million by volume (ppmv). In some cases, hydrogen and acetate concentrations were elevated concurrently, suggesting a positive correlation. However, sampling of hydrogen was limited in comparison to acetate, so any relationships between the two analytes, if present, were difficult to discern.

Pore-water chemistry of ODP Leg 115 samples

Analyses of the Sr2+ concentrations of interstitial fluids obtained from sediments squeezed during Leg 115 were used to estimate the rates and total amount of recrystallization of biogenic carbonates. The total amount of recrystallization calculated using this method varies from less than 1 % in sediments at Site 706 to more than 40% at Site 709 in sediments of 47 Ma. Five of the sites drilled during Leg 115 (Sites 707 through 711) were drilled in a depth transect within a restricted geographic area so that theoretically they received similar amounts of sediment input. Of these, the maximum rate of recrystallization occurred in the upper 50 m of Site 710 (3812 m). The amount of recrystallization decreased with increasing water depth at Sites 708 (4096 m) and 711 (4428 m), presumably as a result of the fact that most of the reactive calcium carbonate was dissolved before burial. We also observed significant alkalinity deficits at many of these sites, a condition which most likely resulted from the precipitation of calcium carbonate either in the sedimentary column, or during retrieval of the core. Precipitation of CaCO3 as a result of pressure changes during core retrieval was confirmed by the comparison of Ca2+ and alkalinity from water samples obtained using the in-situ sampler and squeezed from the sediments. At Sites 707 and 716, the shallowest sites, no calcium or alkalinity deficits were present. In spite of our estimations of as much as 45% recrystallization at Site 709, all the carbonate sites exhibited what would be previously considered conservative Ca2+/Mg2+ profiles, which varied from -1 to -0.5. By virtue of the position of these sites relative to known basaltic basement or through the actual penetration of basalt (i.e., Sites 706, 707 and 712), these sites are all known to be underlain by basalt. Our results suggest, therefore, that more positive Ca2 + /Mg2+ gradients cannot necessarily be used as indicators of the nature of basement material.

Permeable coral reef sediment dissolution driven by elevated pCO2 and pore water advection

Ocean acidification (OA) is expected to drive the transition of coral reef ecosystems from net calcium carbonate (CaCO3) precipitating to net dissolving within the next century. Although permeable sediments represent the largest reservoir of CaCO3 in coral reefs, the dissolution of shallow CaCO3 sands under future pCO2 levels has not been measured under natural conditions. In situ, advective chamber incubations under elevated pCO2 (~800 µatm) shifted the sediments from net precipitating to net dissolving. Pore water advection more than doubled dissolution rates (1.10 g CaCO3/m**2/day) when compared to diffusive conditions (0.42 g CaCO3/m**2 /day). Sediment dissolution could reduce net ecosystem calcification rates of the Heron Island lagoon by 8% within the next century, which is equivalent to a 25% reduction in the global average calcification rate of coral lagoons. The dissolution of CaCO3 sediments needs to be taken into account in order to address how OA will impact the net accretion of coral reefs under future predicted increases in CO2.