The concept of using dynamic excitation to enhance fluid flow in a porous medium began to arise in the mid-twentieth century. The initial spark of interest in the subject spurred numerous laboratory investigations throughout the latter half of the twentieth century to identify the mechanisms at work, and to develop field techniques for practical application of the technology. Several prominent laboratory and field studies have been published; however, there are some deficiencies that facilitate the need for further investigation. Groundwater flow and soil dynamics are two distinct areas of research. There is little in common between the two subjects and there is no consideration of soil dynamic properties in any of the reviewed papers. This study will attempt to bridge the gap between these two areas of research. The objective of this research is to attempt to determine how dynamic excitation of a soil matrix affects saturated single-phase fluid flow. This question is investigated through an extensive literature review of previous studies conducted on this topic, as well as through experimentation designed to replicate the mechanisms responsible for this phenomenon. Experimentation on coarse soil samples is conducted using a modified Stokoe-type resonant column device that allows a quantification of the effects of torsional and axial excitation, frequency of vibration, and strain level. This type of testing in the both the torsional and axial mode has never been conducted before using a resonant column; the Poisson ratios computed using the complimentary data has never been published in the literature.

The concept of using dynamic excitation to enhance fluid flow in a porous medium began to arise in the mid-twentieth century. The initial spark of interest in the subject spurred numerous laboratory investigations throughout the latter half of the twentieth century to identify the mechanisms at work, and to develop field techniques for practical application of the technology. Several prominent laboratory and field studies have been published; however, there are some deficiencies that facilitate the need for further investigation. Groundwater flow and soil dynamics are two distinct areas of research. There is little in common between the two subjects and there is no consideration of soil dynamic properties in any of the reviewed papers. This study will attempt to bridge the gap between these two areas of research. The objective of this research is to attempt to determine how dynamic excitation of a soil matrix affects saturated single-phase fluid flow. This question is investigated through an extensive literature review of previous studies conducted on this topic, as well as through experimentation designed to replicate the mechanisms responsible for this phenomenon. Experimentation on coarse soil samples is conducted using a modified Stokoe-type resonant column device that allows a quantification of the effects of torsional and axial excitation, frequency of vibration, and strain level. This type of testing in the both the torsional and axial mode has never been conducted before using a resonant column; the Poisson ratios computed using the complimentary data has never been published in the literature.

No description available.

No description available.

The flow of wood fiber suspensions plays an important role during the pulp and paper manufacture process. Considerable research has been carried out in the past 50 years to characterize the fiber suspension flow behavior and to monitor the fiber suspension flow during paper manufacture. However, the above research has been hampered by the lack of techniques to directly characterize fiber suspension flow fields because fibers and fiber flocs tend to interfere with instruments inserted into the flow. The fundamental studies in this thesis concentrated on three parts: (1) examine the feasibility of measuring wood fiber suspension flow by Pulsed Ultrasonic Doppler Velocimetry (PUDV), (2) apply PUDV to characterize fiber suspension flow behavior in a rectangular channel, (3) apply PUDV to measure the forming jet velocity profile along the jet thickness direction (ZD). In the first part, it is demonstrated that PUDV is an accurate technique for the velocity profile measurement of fiber suspension flow. The measurement has high repeatability and sensitivity. Suitable parameters should be selected in order to obtain the optimum measuring results.

This thesis presents a preconditioned solenoidal basis method to solve the algebraic system arising from the linearization and discretization of primitive variable formulations of Navier-Stokes equations for incompressible fluid flows. The system is restricted to a discrete divergence-free space which is constructed from the incompressibility constraint. This research work extends an earlier work on the solenoidal basis method for two-dimensional flows and three-dimensional flows that involved the construction of the solenoidal basis P using circulating flows or vortices on a uniform mesh. A localized algebraic scheme for constructing P is detailed using mixed finite elements on an unstructured mesh. A preconditioner which is motivated by the analysis of the reduced system is also presented. Benchmark simulations are conducted to analyze the performance of the proposed approach.

This research is focused on development and enhancement of the model of fluid flow in a formation with stress-dependent permeability. Several typical axi-symmetrical problems of fluid flow in a multi-layered reservoir with account for wellbore storage and skin have been solved numerically. The permeability was assumed to be a function of the vertical deformation of the reservoir. This deformation is the result of changing stress-strain state in the elastic system, comprised of the reservoir itself and the surrounding rock mass. The change in the stress-strain state of the system is induced by pressure change in the layers of the reservoir. Numerical results qualitatively agree with observed field behavior. Such behavior includes (1) deviation of an inflow performance curve from the straight-line relationship at pressures above bubble-point pressure, (2) time- and rate-dependence of well-testing derivative, (3) asymmetry of processes of production and of injection, and (4) inconsistent results between drawdown and buildup, or injection and falloff tests. Based on the results, a procedure to estimate the parameters of the suggested permeability model is proposed.

Mudstones compose nearly 70% of the volume of sedimentary basins, yet they are among the least studied of sedimentary rocks. Their low permeability and high compressibility contribute to overpressure around the world. Despite their fundamental importance in geologic processes and as seals for anthropogenic-related storage, a systematic, process-based understanding of the interactions between porosity, compressibility, permeability, and pore-size distribution in mudstones remains elusive. I use sediment mixtures composed of varying proportions of natural mudstone such as Boston Blue Clay or Nankai mudstone and silt-sized silica to study the effect of composition on permeability and compressibility during burial. First, to recreate natural conditions yet remove variability and soil disturbance, I resediment all mixtures in the laboratory to a total stress of 100 kPa. Second, in order to describe the systematic variation in permeability and compressibility with clay fraction, I uniaxially consolidate the resedimented samples to an effective stress equivalent to about 2 km of burial under hydrostatic conditions. Scanning electron microscope images provide insights on microstructure. My experiments illuminate the controls on mudstone permeability and compressibility. At a given porosity, vertical permeability increases by an order of magnitude for clay contents ranging from 59% to 34% by mass whereas compressibility reduces by half at a given vertical effective stress. I show that the pore structure can be described by a dual-porosity system, where one rock fraction is dominated by silt where large pores are present and the majority of flow occurs and the other fraction is dominated by clay where limited flow occurs. I use this concept to develop a coupled compressibility-permeability model in order to predict porosity, permeability, compressibility, and coefficient of consolidation. These results have fundamental implications for a range of problems in mudstones. They can be applied to carbon sequestration, hydrocarbon trapping, basin modeling, overpressure distribution and geometry as well as morphology of thrust belts, and an understanding of gas-shale behavior. / text

To understand the fluid temperature and pressure during the Ouachita orogeny, we used isotopic analysis of syntectonic veins and adjacent host material, quartz-calcite oxygen isotope thermometry and fluid inclusion analysis. The veins were at or near isotopic equilibrium with their host rocks; neither the host nor veins has been isotopically reset. The average isotopic variation in (delta18)O between vein and host is 2.4 plus/minus 1.7% and 0.7 plus/minus 1.7% for quartz and calcite, respectively. The temperature of vein formation from quartz-calcite oxygen isotope thermometry is about 210-430 degrees C. Although this is a large range, the temperature does not vary systematically in the exposed Ordovician through Mississippian rocks. The lack of isotopic difference between host and vein suggests that the host oxygen determined that of the veins. This in turn suggests that the fluid in the rocks did not change regionally. The vitrinite reflectance/temperature of the host rocks increases with restored stratigraphic depth more than that calculated with the quartz-calcite thermometer in veins. Fluid inclusion analysis in vein quartz constrains homogenization temperatures to be from 106-285 degrees C. Isochores from fluid inclusion analyses were constrained using quartz-calcite thermometry and vitrinite reflectance temperatures to calculate vein formation pressures of 0.3?4.7 kbars. These pressures correspond to vein formation depths up to 19 km, assuming an unduplicated stratigraphic section. Using burial curves and a reasonable range of geothermal gradients, vein formation ages are between 300 to 315 Ma, i.e., Early to Middle Pennsylvanian.

This thesis investigates normal faulting and its influence on fluid flow over a wide range of spatial and temporal scales using tunnel engineering geological logs, outcrop, surface fault traces, earthquakes, gravity, and volcanic ages. These data have been used to investigate the impact of faults on fluid flow (chapter 2), the geometry and kinematics of the Taupo Rift (chapter 3), the hydration and dehydration of the subducting Pacific plate and its influence on the Taupo Volcanic Zone (chapter 4), the migration of arc volcanism across the North Island over the 16 Myr and the associated changes in slab geometry (chapter 5) and the Pacific-Australia relative plate motion vectors since 38 Ma and their implications for arc volcanism and deformation along the Hikurangi margin (chapter 6). The results for each of these five chapters are presented in the five paragraphs below. Tunnels excavated along the margins of the southern Taupo Rift at depths < 500 m provide data on the spatial relationships between faulting and ground water flow. The geometry and hydraulic properties of fault-zones for Mesozoic basement and Miocene strata vary by several orders of magnitude approximating power-law distributions with the dimensions of these zones dependent on many factors including displacement, hostrock type and fault geometries. Despite fault-zones accounting for a small proportion of the total sample length (≤ 15%), localised flow of ground water into the tunnels occurs almost exclusively (≥ 91%) within, and immediately adjacent to, these zones. The spatial distribution and rate of flow from fault-zones are highly variable with typically ≤ 50% of fault-zones in any given orientation flowing. The entire basement dataset shows that 81% of the flow-rate occurs from fault-zones ≥ 10 m wide, with a third of the total flow-rate originating from a single fault-zone (i.e. the golden fracture). The higher flow rates for the largest faults are interpreted to arise because these structures are the most connected to other faults and to the ground surface. The structural geometry and kinematics of rifting is constrained by earthquake focal mechanisms and by geological slip and fault mapping. Comparison of present day geometry and kinematics of normal faulting in the Taupo Rift (α=76-84°) with intra-arc rifting in the Taranaki Basin and southern Havre Trough show, that for at least the last 4 Myr, the slab and the associated changes in its geometry have exerted a first-order control on the location, geometry, and extension direction of intra-arc rifting in the North Island. Second-order features of rifting in the central North Island include a clockwise ~20° northwards change in the strike of normal faults and trend of the extension direction. In the southern rift normal faults are parallel to, and potentially reactivate, Mesozoic basement fabric (e.g., faults and bedding). By contrast, in the northern rift faults diverge from basement fabric by up to 55° where focal mechanisms indicate that extension is achieved by oblique to right-lateral strike-slip along basement fabric and dip-slip on rift faults. Hydration and dehydration of the subducting Pacific plate is elucidated by earthquake densities and focal mechanisms within the slab. The hydration of the subducting plate varies spatially and is an important determinant for the location of arc volcanism in the overriding plate. The location and high volcanic productivity of the TVZ can be linked to the subduction water cycle, where hydration and subsequent dehydration of the subducting oceanic lithosphere is primarily accomplished by normal-faulting earthquakes. The anomalously high heat flow and volcanic productivity of the TVZ is spatially associated with high rates of seismicity in the underlying slab mantle at depths of 130-210 km which can be tracked back to high rates of deeply penetrating shallow intraplate seismicity at the trench in proximity to oceanic fluids. Dehydration of the slab mantle correlates with the location and productivity of active North Island volcanic centres, indicating this volcanism is controlled by fluids fluxing from the subducting plate. The ages and locations of arc volcanoes provide constraints on the migration of volcanism across the North Island over the last 20 Myr. Arc-front volcanoes have migrated southeast by 150 km in the last 8 Ma (185 km since 16 Ma) sub-parallel to the present active arc. Migration of the arc is interpreted to mainly reflect slab steepening and rollback. The strike of the Pacific plate beneath the North Island, imaged by Benioff zone seismicity (50-200 km) and positive mantle velocity anomalies (200-600 km) is parallel to the northeast trend of arc-front volcanism. Arc parallelism since 16 Ma is consistent with the view that the subducting plate beneath the North Island has not rotated clockwise about vertical axes which is in contrast to overriding plate vertical-axis rotations of ≥ 30º. Acceleration of arc-front migration rates (~4 mm/yr to ~18 mm/yr), eruption of high Mg# andesites, increasing eruption frequency and size, and uplift of the over-riding plate indicate an increase in the hydration, temperature, and size of the mantle wedge beneath the central North Island from ~7 Ma. Seafloor spreading data in conjunction with GPlates have been used to generate relative plate motion vectors across the Hikurangi margin since 38 Ma. Tracking the southern and down-dip limits of the seismically imaged Pacific slab beneath the New Zealand indicates arc volcanism in Northland from ~23 Ma and the Taranaki Basin between ~20 and 11 Ma requires Pacific plate subduction from at (or beyond) the northern North Island continental margin from at least 38 Ma to the present. Pacific plate motion in a west dipping subduction model shows a minimum horizontal transport distance of 285 km preceding the initiation of arc volcanism along the Northland-arc normal to the motion vector, a distance more than sufficient for self-sustaining subduction to occur. Arc-normal convergence rates along the Hikurangi margin doubled from 11 to 23 mm/yr between 20 and 16 Ma, increasing again by approximately a third between 8 and 6 Ma. This latest increase in arc-normal rates coincided with changes in relative plate motions along the entire SW Pacific plate boundary and steepening/rollback of the Pacific plate.