Dissertations / Theses on the topic 'Subsurface stress'

Aujourd’hui, l’un des principaux défis dans l’industrie pétrolière, et particulièrement dans le domaine de l’exploration, est l’exploitation des nouvelles ressources dans des zones structuralement complexes.Nous savons que la géométrie et le glissement le long des failles actives modifie la distribution locale des contraintes. La connaissance du champ de contrainte perturbé actuel est importante pour l’étude des tremblements de Terre, pour la planification de forages ainsi que pour la prédiction de la fracturation induite par l’hydro-fracturation et la prédiction de la réactivation des fractures. Les contraintes perturbées passées, quant à elles sont responsables du développement des fractures naturelles (préexistantes). La détection et la modélisation de celles-ci sont essentielles tant dans le domaine pétrolier (migration et piégeage des fluides) pour une exploitation efficace et à moindre coût des réserves naturelles.Comprendre et quantifier le développement spatial et temporel de ces contraintes a un impact économique non négligeable. L'analyse des paléo-contraintes a été introduite dans un premier temps par Anderson (1905 & 1942) de manière intuitive, puis dans le milieu du siècle dernier, Wallace (1951) et Bott ( 1959) proposèrent les simples postulats que le champ de contrainte est homogène et que la direction du rejet est parallèle à la traction projetée sur le plan de faille. Beaucoup de méthodes d’inversion de contraintes reposent aujourd’hui sur ce principe.Afin d’étudier la validité de l’hypothèse Wallace et Bott, une comparaison avec les vecteurs de glissement générés à partir de modèles numériques (BEM) est effectuée. En testant l’influence de multiples paramètres (géométrie, contraintes au limites du modèle, friction, coefficient de poisson, demi-espace, pression de fluide dans la faille), il est montré que les failles à géométries complexes soumises à certaines contraintes aux limites peuvent engendrer des vecteurs glissements présentant des écarts important avec les la plus grande contraintes cisaillantes résolue sur le plan de faille. A l’inverse, la présence d’une forte friction de glissement permet, dans certaines conditions, de valider l’hypothèse de Wallace et Bott. On s’attache ensuite à comparer les résultats des inversions de contraintes basées sur l’hypothèse de Wallace et Bott (appelé méthode d’inversion classique de contraintes) avec une méthode géomécanique. Pour cela, une faille à géométrie complexe est utilisée dans une étude de sensibilité (contraintes aux limites, friction, échantillonnage) permettant d’analyser l’incertitude des résultats des deux méthodes d’inversion. Cette analyse est ensuite confrontée à l’étude d’un cas de terrain, montrant les avantages et inconvénients des méthodes d’inversions classiques de contraintes.Un des principaux défis de l’industrie pétrolière est l’exploitation des ressources des zones structuralement complexes telles que les réservoirs naturellement fracturés. Connaitre l’état de contraintes hétérogène passé permet d’optimiser la modélisation de ces fractures naturelles. Etant donné que les glissements sur les failles sont difficiles à observer dans les réservoirs pétroliers, les données de d’orientation de fractures (joints, failles, stylolites) sont naturellement prises en compte lors de l’inversion des contraintes. On montre, en utilisant divers exemples de terrain et d’industrie, que dans de tels cas, l’utilisation d’inversions basée sur la mécanique est beaucoup plus appropriée. Cependant, il est parfois difficile de déterminer le type cinématique de fracture observée le long des puits, et très souvent, les zones étudiées ont subi une tectonique polyphasée. La dernière partie vise donc à résoudre le problème des données de types cinématiques non identifiables (joints, failles, stylolites…) et étend parallèlement l’inversion mécanique des contraintes à la séparation de phases tectoniques
Today, one of the main challenges in the oil industry, especially during the exploration phase, is the exploitation of new resources in structurally complex areas such as naturally fractured reservoirs, salt diapirs, mountain ranges, and unconventional reservoirs.We know that the geometry and sliding along active faults modifies the local stress distribution. Knowing the present day perturbed stress field is important for the study of earthquakes, for the planning of the borehole drilling and stability as well as for the prediction of fractures induced by hydro-fracturing and reactivation of natural fractures. In the other side, perturbed paleostress are responsible for the development of (pre-existing) natural fractures. The detection and modeling of the latter, are essential both in the oil industry (migration and trapping of fluids) for a cost efficient recovery of natural reserves.Understanding and quantifying the spatial and temporal development of the stress distribution has a significant economic and environmental impact. The analysis of paleo-constraints was intuitively introduced first by Anderson (1905 & 1942), then in the middle of the last century, Wallace (1951) and Bott (1959) proposed the simple hypothesis that (i) The stress field is homogeneous in space and constant in time, and that (ii) the slip direction is parallel to the traction projected on the fault plane which gives the direction of the shear stress. Many stress inversion methods are based on this hypothesis while recent studies raise doubts as to their compatibility with rock mechanics.In order to investigate the validity of the Wallace and Bott hypothesis, a comparison with vectors of slip generated with numerical models (BEM) is performed. By testing the influence of multiple parameters (geometry, boundary conditions, friction, Poisson’s coefficient , half-space, fault fluid pressure), it is shown that the complex geometry faults subject to specific boundary conditions can yield slip vectors with significant discrepancies with the maximum shear stress resolved on the fault plane. Conversely, the presence of a high sliding friction, allows under certain conditions, to validate the hypothesis of Wallace and Bott.We then focus on the task to compare the results of stress inversions based on the assumption of Wallace and Bott (called classical stress inversion methods) to a geomechanical method. For this, a complex fault geometry is used in a sensitivity analysis (boundary conditions, friction, sampling) to evaluate the uncertainty of the results of the two inversion methods. This analysis is then compared to a case study, Chimney Rock (Utah, USA), showing the advantages and disadvantages of the classical stress inversion methods.One of the main challenges of the oil industry is the exploitation of resource in structurally complex oil fields such as naturally fractured reservoirs. Knowing the heterogeneous paleostress allows to optimize the modeling of these natural fractures. Since slip on faults is hardly observed in petroleum reservoirs, fracture orientation data (joints, faults, stylolites) are naturally taken into account during the inversion of stresses. It is shown, using various field and industry examples, that in such cases the use of mechanical stress inversions is much more appropriate.However, it is sometimes difficult to determine the fracture kinematics observed along wellbores, and very often the studied regions underwent multiple tectonic phases. The final section aims to address the problem of data with unknown kinematic (joints, faults, stylolites ...) and expends the mechanical stress inversion to the separation of tectonic phases

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences, 2013.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 259-271).
We study the effect of fractures on reservoir characterization and subsurface rock property measurements using seismic data. Based on the scale of a fracture relative to seismic wavelength, we divide the dissertation into two parts: larger scale fractures and microcracks. In the first part, we study the sensitivity of seismic waves and their time-lapse changes in hydraulic fracturing to the geometrical and mechanical properties of fractures that have dimensions comparable to the seismic wavelength. Through our analysis, we give the general seismic response of a fracture with a linear slip boundary and introduce the fracture sensitivity wave equation for optimal time-lapse survey design. Based on the characteristics of scattering from fractures, we develop an approach to determine the fracture properties using scattered seismic waves. The applicability and accuracy of our method is validated through both numerical simulations and laboratory experiments. Application of our approach to the Emilio Field shows that two orthogonal fracture systems exist and the field data results are consistent with well data. In the second part, we study the effects of microcracks and in situ stress on the formation properties measured from borehole sonic logging. Formation property measurements in a borehole could be biased by the borehole stress concentration, which alters the near wellbore formation properties from their original state. To study this problem, we first develop an iterative approach, which combines a rock physics model and a finite-element method, to calculate the stress-dependent elastic properties of the rock around a borehole when it is subjected to an anisotropic stress loading. The validity of this approach is demonstrated through a laboratory experiment on a Berea sandstone sample. We then use the model obtained from the first step and a finite-difference method to simulate the acoustic response in a borehole. We compare our numerical results with published laboratory acoustic wave measurements of the azimuthal velocity variations along a borehole under uniaxial loading and find very good agreement. Our results show that the variation of P-wave velocity versus azimuth is different from the presumed cosine behavior due to the preference of the wavefield to propagate through a higher velocity region.
by Xinding Fang.
Ph.D.

Biologically-mediated subsurface remediation by biofilm communities is a poorly understood process that is spatially and temporally dynamic. Two microbial responses, catabolism and the stress response glutathione-gated potassium efflux (GGKE), to benzene, pentachlorophenol (PCP), or Cd exposure were studied in up-flow sand columns to examine the contribution of each response to the overall functional response of a subsurface biofilm. Benzene was catabolized in the aerobic zone, and did not activate the GGKE response, and exhibited the highest biomass concentrations of all columns. PCP was not catabolized during this study, but was found to elicit two responses, oxidative phosphorylation uncoupling and GGKE, that appeared to be concentration dependent. Oxidative uncoupling was the controlling metabolic response up to 10 mg/L PCP, while the GGKE stress response was activated near 20 mg/L PCP. PCP column biomass did not show long-term biomass detachment, although immediate detachment occurred during initial GGKE activation. Cd column biomass activated the GGKE response as perturbing Cd concentrations increased. Extracellular polymeric substance (EPS)-Cd complexation was a possible detoxification mechanism, as biomass concentrations did not decrease with increasing Cd concentration, and increased as Cd concentrations decreased. Results of this study suggested that the increased exposure of electrophilic contaminants to sand column biomass did not cause biomass detachment.
Master of Science

Subsurface volume loss takes place in many geotechnical situations, and it is inherently accompanied by complex stress and displacement fields that may influence the performance of engineered geosystems. This research is a deformation-centered analysis, it depends on soil compressibility and it is implemented using finite elements. Soil stiffness plays a central role in predicting ground deformation. First, an enhanced Terzaghi’s soil compressibility model is proposed to satisfy asymptotic conditions at low and high stress levels with a small number of physically meaningful parameters. Then, the difference between small and large strain stiffness is explored using published small and large-strain stress-strain data. Typically, emphasis is placed on the laboratory-measured stiffness or compressibility; however, there are pronounced differences between laboratory measurements and field values, in part due to seating effects that prevail in small-thickness oedometer specimens. Many geosystems are subjected to repetitive loads; volumetric strains induced by drained repetitive ko-loads are experimentally investigated to identify shakedown and associated terminal density. The finite element numerical simulation environment is used to explore the effect of localized subsurface mass loss on free-surface deformation and shallow foundations settlement and bearing capacity. A stress relaxation module is developed to reproduce the change in stress associated to dissolution features and soft zone formation. The comprehensive parametric study is summarized in terms of dimensionless ratios that can be readily used for engineering applications. Field settlement data gathered at the Savannah River Site SRS are back-analyzed to compare measured values with predictions based on in situ shear wave velocity and strain-dependent stiffness reduction. The calibrated model is used to estimate additional settlements due to the pre-existing cavities, new cavities, and potential seismic events during the design life of the facility.

Objective: The goal of the present study was to investigate the effect of chairside surface treatments on biaxial flexural strength and subsurface damage of monolithic zirconia ceramics. Methods: Specimens (15x15x1.2 mm3) were prepared by sectioning from commercially available zirconia blanks (BruxZirTM) and sintering according to manufacturer's recommendations. Fully dense specimens were randomly divided into five groups (n=30) and treated as follows; 1) as-sintered (AS) 2) air abraded with 50 μm alumina fine particles (AAF), 3) air abraded with 250 μm coarse alumina particles (AAC), 4) ground (G), and 5) ground and polished (GP) to mimic chairside and dental laboratory treatments. Microstructural changes were thoroughly characterized by optical and scanning electron microscopy, surface profilometry and atomic force microscopy. Crystalline phases and their depth profile were investigated by x-ray diffraction (XRD) and grazing incidence x-ray diffraction (GIXRD). Results were analyzed by Kruskal-Wallis test and Tukey's adjustment for multiple comparisons. A 0.05 level of significance was used. Reliability was evaluated by Weibull analysis. Results: All treatment groups exhibited a significant difference in mean surface roughness (Rq) compared to the as-sintered group (p<0.05). The AAC group showed the highest surface roughness at 1.08 ± 0.17 μm, followed by the G, AAF and AS groups. The GP group exhibited the lowest surface roughness. The group air abraded with fine particles showed the highest mean biaxial flexural strength (1662.62 ± 202.58 MPa), but was not different from the ground and polished group (1567.19 ± 209.76 MPa). The groups air abraded with coarse particles or ground with diamond bur exhibited comparable mean biaxial flexural strength at 1371.37 ±147.62 MPa and 1356.98 ±196.77 MPa, respectively. The as-sintered group had the lowest mean biaxial flexural strength at 1202.29 ±141.92 MPa. The depth of compressive stress layer, measured by GIXRD was approximately 50 μm in the AAF group, followed by the AAC group with ~35 μm, ~10 μm for the ground group and ~5 μm for the ground and polished group. Deep subsurface cracks were observed in the AAC group (~80 μm in depth) and G group (~25 μm in depth), whereas shallower flaws were present in the AAF and GP groups at 10 and 3 μm, respectively. Weibull analysis represented a greater reliability in zirconia specimens treated with air abrasion groups. Conclusions: Surface treatments induced the t-m transformation in 3Y-TZP and associated development of compressive stresses to a depth that varied with the severity of the treatment performed. GIXRD revealed that AAF led to the thickest compressive stress layer, followed by AAC, G and GP. SEM revealed that subsurface damage was most severe with AAC, followed by G, AAF and GP. We propose that the flexural strength results can be explained by the difference between the depth of the compressive stress layer induced by the transformation and the depth of the subsurface flaws.

Doctor of Philosophy
Department of Biological & Agricultural Engineering
Aleksey Y. Sheshukov
As the world faces an increasing demand for food due to the growing global population and the pernicious effects of land degradation, there is a need to overcome this challenge by using sustainable management practices for agricultural productions. One of the problems, which sustainable agriculture seeks to address, is the loss of topsoil due to soil erosion. Changing weather patterns also contribute to the average annual rainfall across the globe with an excess precipitation, which creates runoff and causes soil erosion. One of the significant yet less studied types of soil erosion is ephemeral gully erosion. Formed by the concentrated overland flow during intensive rainfall events, ephemeral gullies are channels on agricultural fields that can be removed by tillage operations but appear at the same location every year. Even though simplified ephemeral gully models estimate soil losses, they do not account for complicated hydrological and soil erosion processes of channel formations. The purpose of this research work is to investigate sediment sources and develop tools that can predict ephemeral gully erosion more efficiently. To achieve this goal, an experimental study was conducted on an agricultural field in central Kansas by tracking channel development, monitoring soil moisture content, and recording the amount of rainfall. Runoff and sediment loads from contributing catchment and critical and actual shear stresses were estimated by the computer model, and conclusions were made on the effect of saturation dynamics on the erosion processes. Furthermore, a two-dimensional subsurface water flow and soil erosion model was developed with the variable soil erodibility parameters which account for the subsurface fluxes and the effects on the soil detachment process. The model was applied to study the impacts of variable soil erodibility parameters on the erosion process for different soils and various antecedent soil moisture conditions. Also developed to estimate the soil losses at the field scale was an integrated spatially-distributed ephemeral gully model with dynamic time-dependent channel development. The model showed good fit by matching the experimental data. The results from this work can be used to advance the research of soil erosion prediction from concentrated flow channels and ephemeral gullies formed on agricultural fields.

Dynamically recrystallized quartz microstructures preserved in contact aureoles allow for stress and strain rate estimates associated with penetrative deformation of rocks surrounding pluton margins. Microstructural analysis of the Harkless quartzites surrounding the western margin of Papoose Flat pluton indicates that recrystallization occurred by grain boundary migration with mean recrystallized grain size ranging from 86-225 µm. The application of three calibrated piezometers results in differential stress estimates between ~11 and ~29 MPa. Published wet-quartzite dislocation creep flow laws combined with deformation temperature, water fugacity, and differential stress estimates infer strain rates that range from 1.2 x 10⁻¹⁴ s⁻¹ to 2.3 x 10⁻¹² s⁻¹. In order to analyze 3D subsurface fault geometry along map-pattern curves (salients and recesses), a structure contour map of the Moine thrust, extending from the North Coast southwards to the Dundonnel area, was constructed from 1:50,000 scale British Geological Survey (BGS) maps by correlating between elevation control points constrained by the intersection of the fault trace with topographic contours. The structure contour map indicates significant lateral variation in fault geometry along the Moine thrust, with recesses associated with antiformal corrugations in the subsurface and salients characterized by planar geometries or broad synformal corrugations. Additionally, structure contour maps constructed on the Glencoul thrust, as depicted by original BGS maps confirms that the thrust segments to the NE and SW of Loch Glencoul are part of the same structure, rather than different structures separated by a lateral ramp as shown on more recent BGS maps.
Master of Science

The primary objectives of this study were to address how stream temperature is influenced by (1) spatial variability in energy exchanges, (2) reach-scale stream-subsurface water interactions and (3) the net radiation dynamics associated with standing dead riparian vegetation. Stream temperature, riparian microclimate, and hydrology were characterized for a 1.5 km reach of Fishtrap Creek, located north of Kamloops, British Columbia. Within-reach air temperature and humidity variability was small, while wind speed, net radiation and surface-subsurface interactions exhibited considerable spatially variability. The field data were used to drive a deterministic energy budget model to predict stream temperature. The model was evaluated against measured stream temperature and performed well. The model indicated that the spatially complex hydrology was a significant control on the observed stream temperature patterns. A modelling exercise using three canopy cover scenarios revealed that post-disturbance standing dead trees reduce daytime net radiation reaching the stream surface by one third compared to complete vegetation removal. However, standing dead trees doubled daytime net radiation reaching the stream compared to pre-wildfire conditions. The results of this study have highlighted the need to account for the spatial variability of energy exchange processes, specifically net radiation and surface-subsurface water interactions, when understanding and predicting stream thermal regimes.

The gravel environment of 30 brown trout (Salmo trutta) redds and adjacent nonredd sites in two western North Carolina streams were studied during the incubation period in 1979-1980 and 1980-1981. Intragravel water temperature, dissolved oxygen concentration, and percent oxygen saturation were highly correlated with surface water measurements, indicating that intragravel water is of surface origin. Permeability ranged from 250 to 149,350 cm/hr and averaged 6,150 cm/hr. Apparent velocity varied from 0 to 1,000 cm/hr and averaged 30 cm/hr. Permeability in redds was significantly greater than at nonredd sites. No significant differences in apparent velocity were found between redd and nonredd sites. No consistent differences in permeability or apparent velocity were found between streams or over time. Permeability and apparent velocity decreased significantly with depth. Freeze cores were collected from redd and nonredd sites and divided into three 10-cm layers for analysis. Geometric mean diameter, sorting coefficient, fredle index, percent fines <2.00 mm, and percent porosity were highly variable and averaged 11.8 mm, 2.8, 4.2, 17.0 percent, and 19.0 percent, respectively. No significant differences were found among factors tested. Correlations between these gravel indices and permeability and apparent velocity were low. The gravel and intragravel environments appeared to be adequate for larval survival. Measurements did not reveal any clear trends during the incubation period. Brown trout did not by choice or redd construction appear to select or create (by redd construction) a subsurface environment different from the surrounding stream bed.
M.S.

Preferential flow paths are areas of substantially higher permeability than surrounding media. Macropores and soil pipes are a type of preferential flow path where conduit-like voids in the subsurface are typically greater than three millimeters in diameter. They are known to occur in agricultural and forest soils, often as a result of biological and physical processes. Macropores also exist in stream banks and have the potential to enhance the exchange of water and solutes between the channel and riparian groundwater, yet the geographic distribution of bank macropores is unknown. Here we determined the abundance, distribution, and geometry of naturally occurring surface-connected macropores in the banks of 20 streams across five physiographic provinces in the Eastern United States. We identified a total of 1,748 macropores, which were present in all 20 streams, with 3.8 cm average width, 3.3 cm average height, 11.5 cm average depth, and 27.9 cm average height above water surface elevation. Macropore abundance, distribution and geometry were statistically different between physiographic provinces, stream orders, and soil textures, with the latter being the most important. Macropores tended to be larger and more abundant in soils with a high cohesiveness and a low hydraulic conductivity compared to soils with a low cohesiveness and high hydraulic conductivity. As a result, streams with greater longitudinal heterogeneity of soil texture also had greater heterogeneity of macropore density. However, macropore size and height above baseflow water surface elevation also increased with stream order and therefore stream size. This work represents the first attempt to characterize macropores across a variety of riverine systems and presents evidence that macropores may play an important role in hyporheic exchange within stream banks. These results may have water quality implications, where macropores may enhance hyporheic exchange yet reduce the filtering capacity of riparian buffer zones.
MS

Water quality in streams around the world continues to be degraded by a series of human activities that feed pollutants into the vulnerable stream ecosystem via surface and subsurface runoff. This continues to accelerate global biodiversity and habitat losses within the stream environments and across entire watersheds with net adverse effects on public health and the ability of communities and ecosystems to adapt or become resilient to the prevalent impacts of climate change. One commonly used approach for protecting stream water from pollution is the use of vegetated riparian buffer zones to mitigate pollutants in surface and subsurface runoff prior to runoff entry into the stream channel. The optimal success of this approach requires land and water resource managers to understand the mechanisms by which riparian buffer zones function and the full range of factors that influence the effectiveness of riparian buffer vegetation in abating stream water pollution. Despite this need, resource managers in different geographical locations around the world still struggle to understand the linkages between riparian vegetation and stream chemical quality. This literature review therefore sought to synthesize findings from various scientific articles on the ways in which the major attributes of riparian vegetation [type, age, width, restoration and shading effect] influence the effectiveness of riparian vegetation in protecting the chemical quality of water in streams. This was aimed at generating conclusions and perspectives that could improve academic knowledge and natural resource managers’ understanding of the intricate linkage between riparian vegetation and changes in water chemistry. The study finds that the factors of riparian vegetation type, age, width, restoration and shading effects require due consideration in the development of riparian buffer zone and stream water chemical quality management interventions. I find that these factors require a high degree of integration, triangulation and context-specificity to achieve the objectives of riparian management intervention. I further find that stream water quality decision-making processes need to combine riparian vegetation-based approaches with other measures for mitigating and containing the spillage of pollutants at the source.

Presentation was conducted via Skype

Accurate models to reliably predict sediment and chemical transport in watershed water systems enhance the ability of environmental scientists, engineers and decision makers to analyze the impact of contamination problems and to evaluate the efficacy of alternative remediation techniques and management strategies prior to incurring expense in the field. This dissertation presents the conceptual and mathematical development of a general numerical model simulating (1) sediment and reactive chemical transport in river/stream networks of watershed systems; (2) sediment and reactive chemical transport in overland shallow water of watershed systems; and (3) reactive chemical transport in three-dimensional subsurface systems. Through the decomposition of the system of species transport equations via Gauss-Jordan column reduction of the reaction network, fast reactions and slow reactions are decoupled, which enables robust numerical integrations. Species reactive transport equations are transformed into two sets: nonlinear algebraic equations representing equilibrium reactions and transport equations of kinetic-variables in terms of kinetically controlled reaction rates. As a result, the model uses kinetic-variables instead of biogeochemical species as primary dependent variables, which reduces the number of transport equations and simplifies reaction terms in these equations. For each time step, we first solve the advective-dispersive transport of kinetic-variables. We then solve the reactive chemical system node by node to yield concentrations of all species. In order to obtain accurate, efficient and robust computations, five numerical options are provided to solve the advective-dispersive transport equations; and three coupling strategies are given to deal with the reactive chemistry. Verification examples are compared with analytical solutions to demonstrate the numerical accuracy of the code and to emphasize the need of implementing various numerical options and coupling strategies to deal with different types of problems for different application circumstances. Validation examples are presented to evaluate the ability of the model to replicate behavior observed in real systems. Hypothetical examples with complex reaction networks are employed to demonstrate the design capability of the model to handle field-scale problems involving both kinetic and equilibrium reactions. The deficiency of current practices in the water quality modeling is discussed and potential improvements over current practices using this model are addressed.
Ph.D.
Department of Civil and Environmental Engineering
Engineering and Computer Science
Civil Engineering

Compte tenu de l’évolution démographique et climatique planétaire, la gestion de la ressource en eau constitue un défi majeur auquel la communauté internationale devra faire face au cours du XXIème siècle. A cet effet, l'identification de la continuité hydrique entre les eaux de surface et les eaux souterraines permet l'introduction de la notion de gestion intégrée de la ressource. L'application de ce principe de gestion au bassin de la Seine, à travers l'estimation des échanges nappe-rivière, est rendue possible par la mise en pratique du concept d'interface nappe-rivière emboitées au sein de travaux de modélisation. Pour cela une procédure de modélisation multi-échelles peut être mise en place. Elle vise à intégrer des informations locales au sein de modélisation à une échelle supérieure. Dans ce mémoire, une procédure de modélisation multi-échelles est mise en œuvre. Ce protocole de modélisation est initié par une estimation des flux d'eau régionaux au sein de l'hydrosystème Seine. La cohérence globale de ces flux est garantie par le développement d'une méthodologie de calibration de modèles couplés en deux étapes. Ensuite les informations locales, que sont les hétérogénéités de la plaine alluviale de la Bassée et de la représentation des interfaces nappe-rivière du réseau secondaire, sont intégrées au modèle régional par une procédure de modélisation emboitée et de changement d'échelle des paramètres hydrauliques. La mise en place de cette procédure a finalement permis l'estimation fine des échanges nappe-rivière sur la quasi-totalité (83%) du réseau hydrographique naturel du bassin de la Seine et ainsi de répondre aux recommandations de gestion intégrée de la ressource faites par la directive cadre sur l'eau
Given the current climate and anthropogenic evolutions, water management becomes one of the greatest challenges of the 21st Century. For that purpose, by identifying hydraulic continuity between surface and subsurface water, the concept of integrated water management can be introduced. In this work this management concept is applied on the Seine basin by quantizing hydrological processes occuring at the nested stream-aquifer interface. The implementatin of the nested interface concept can bedone through multi-scale modeling. This modelling procedure, aimed at embody the local characteristics of the interfaces (such as structural or hydrodynamic heterogeneities) in large scale models. A multi-scale modelling procedures is applied to the regional Seine basin model (70000 km²) in order, to study the hydrodynamic behaviour of the Bassée alluvial plain, and to quantify the stream-aquifer exchanged fluxes at the basin scale. The modelling protocol is initiated with regionals fluxes estimation over Seine hydrosystem. Regional fluxes consistency are assured by a two-step calibration procedure of fully coupled models. Then, the local characteristics of the Bassée alluvial plain, are implemented in the regional model by nested modelling methodology associated with upscaling procedure of hydraulics properties. Finally, the multi-scale modelling procedure lead to quantify distributed stream-aquifer exchanged water fluxes over 83% of the natural river network of the Seine basin, and thus, achieve to answer the integrated water resources management recommandations of the water framework directive

碩士
國立中央大學
地球物理研究所
100
This study uses subsurface geological data from wells drilled in the past 40 years to estimate the in-situ stress from Hsinchu to Taichung area. The gradient of vertical stress (Sv) estimated from the formation density log data in Tiehchenshan (TCS) and Chuhuangkeng (CHK) region is about 23 MPa/km. From the repeat formation tests (RFT) collected in Chingtsaohu (CTH), Paoshan (PS), Chinshui (CS), Chuhuangkeng and Tiehchenshan regions, pore pressure (Pp) in sandstone of individual structure can be obtained. Additional background pore pressure values in shale can be estimated from wireline Sonic logs. The top of over-pressure as well as amount of erosional uplift can be inferred based in the plot of pore pressures vs. depth. The Chuhuankeng anticline is highest structure in the Hsinchu and Miaoli area in terms of elevater of the over-pressure zone. Hydraulic fracturing data including leak-off tests (LOTs) and fluid injection in cased and cemented wellbores (mini-fracs) in study area provide a good estimate the minimum horizontal stress (Shmin), which has a gradient about 17~19 MPa/km. The magnitude of maximum horizontal stress (SHmax) was constrained by frictional limits from Anderson’s theory of faulting and by presence of drilling induced tensile fractures. The gradient of the SHmax is about 26~35 MPa/km assuming frictional coefficient (μ) of 0.6. The measured In-situ stress magnitudes indicate a predominantly strike-slip stress regime in Hsinchu and Miaoli area. Applying the methodology of judging borehole breakouts in World Stress Map (WMS), orientations of maximum horizontal stresses were calculated from four-arm caliper tools surveyed in 8 wells. The maximum horizontal stress is predominantly orientated in NW-SE. Detailed SHmax orientation will be affected by surface topography and subsurface structures that wells encountered. The magnitude or the orientation of the in-situ stresses vary at shallow depths. The stress state flips from strike-slip and normal fault, and the SHmax (or Shmin) orientations can also change accordingly.