Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

Molecular and isotopic characterization of biomarkers in the Frick Swiss Jura sediments: A palaeoenvironmental reconstruction on the northern Tethys margin

Molecular and stable carbon isotope compositions of source-specific hydrocarbons have been used to reconstruct palaeoenvironmental conditions during deposition of the Middle Hettangian to Upper Sinemurian sediments on the northern epicontinental Tethys margin, Frick Swiss Jura. Increasing algal, cyanobacterial and phytoplanktonic (i.e., dinoflagellate) contributions associated with the C-13-enrichment of cyanobacteria derivatives (i.e., hopanes and monomethylalkanes) suggest enhanced primary productivity upsection. This is related to the C-13-enrichment of dissolved CO2 in the upper layers and the progressive increase of depth and oxygenation of the water column. In the Middle Hettangian shallow-water environments (lagoon), the occurrence of green sulfur bacteria (Chlorobiaceae) derivatives indicates that the lower part of the water column was strictly anoxic and rich in H2S. Since these bacteria require very low light intensity to grow, these euxinic conditions may be extended up to the photic zone, allowing for anaerobic photosynthesis. Light penetration depth is most likely reduced by high productivity and/or turbidity in the photic zone. In these sediments, C-13-depleted hopanoids (-39.5 parts per thousand) are most likely associated with phototrophic purple sulfur bacteria utilizing isotopically light organic carbon at the base of the aerobic zone. These purple sulfur bacteria may have consumed the H2S used by Chlorobiaceae in the deeper layers and thus, sustained the algae and cyanobacteria productivity in the upper layers. The C-13-depleted carbonate (-13.3 parts per thousand) may be partially related to the anaerobic oxidation of the organic matter during bacterial sulfate-reduction. (c) 2006 Elsevier Ltd. All rights reserved.

Imaging molasse and quaternary sediments in Lake Geneva, Switzerland, with 3-D high-resolution seismic reflection methods: A case study

A high-resolution three-dimensional (3-D) seismic reflection survey was conducted in Lake Geneva, near the city of Lausanne, Switzerland, as part of a project for developing such seismic techniques. Using a single 48-channel streamer, the 3-D site with an area of 1200 m x 600 m was surveyed in 10 days. A variety of complex geologic structures (e.g. thrusts, folds, channel-fill) up to similar to150 m below the water bottom were obtained with a 15 in.(3) water gun. The 3-D data allowed the construction of an accurate velocity model and the distinction of five major seismic facies within the Lower Freshwater Molasse (Aquitanian) and the Quaternary sedimentary units. Additionally, the Plateau Molasse (PM) and Subalpine Molasse (SM) erosional surface, "La Paudeze" thrust fault (PM-SM boundary) and the thickness of Quaternary sediments were accurately delineated in 3-D.

Antibiotic resistant bacteria/genes dissemination in lacustrine sediments highly increased following cultural eutrophication of Lake Geneva (Switzerland).

This study investigates faecal indicator bacteria (FIB), multiple antibiotic resistant (MAR), and antibiotic resistance genes (ARGs), of sediment profiles from different parts of Lake Geneva (Switzerland) over the last decades. MARs consist to expose culturable Escherichia coli (EC) and Enterococcus (ENT) to mixed five antibiotics including Ampicillin, Tetracycline, Amoxicillin, Chloramphenicol and Erythromycin. Culture-independent is performed to assess the distribution of ARGs responsible for, β-lactams (blaTEM; Amoxicillin/Ampicillin), Streptomycin/Spectinomycin (aadA), Tetracycline (tet) Chloramphenicol (cmlA) and Vancomycin (van). Bacterial cultures reveal that in the sediments deposited following eutrophication of Lake Geneva in the 1970s, the percentage of MARs to five antibiotics varied from 0.12% to 4.6% and 0.016% to 11.6% of total culturable EC and ENT, respectively. In these organic-rich bacteria-contaminated sediments, the blaTEM resistant of FIB varied from 22% to 48% and 16% to 37% for EC and ENT respectively, whereas the positive PCR assays responsible for tested ARGs were observed for EC, ENT, and total DNA from all samples. The aadA resistance gene was amplified for all the sediment samples, including those not influenced by WWTP effluent water. Our results demonstrate that bacteria MARs and ARGs highly increased in the sediments contaminated with WWTP effluent following the cultural eutrophication of Lake Geneva. Hence, the human-induced changing limnological conditions highly enhanced the sediment microbial activity, and therein the spreading of antibiotic resistant bacteria and genes in this aquatic environment used to supply drinking water in a highly populated area. Furthermore, the presence of the antibiotic resistance gene aadA in all the studied samples points out a regional dissemination of this emerging contaminant in freshwater sediments since at least the late nineteenth century.

Mirdita Zone ophiolites and associated sediments in Albania reveal Neotethys Ocean origin

The Mirdita Ophiolite Zone in Albania is associated with widespread melanges containing components of up to nappe-size. We dated matrix and components of the melange by radiolarians, conodonts, and other taxa. The components consist of radiolarites, pelagic limestones and shallow-water limestones, all of Triassic age, as well as ophiolites. Triassic radiolarite as a primary cover of ophiolite material proves Middle Triassic onset of Mirdita ocean-floor formation. The melange contains a turbiditic radiolarite-rich matrix (''radiolaritic flysch''), dated as Late Bajocian to Early Oxfordian. It formed as a synorogenic sediment during west-directed thrusting of ophiolite and sediment-cover nappes representing ocean floor and underplated fragments of the western continental margin. The tectonic structures formed during these orogenic events (''Younger Kimmeridian or Eohellenic Orogeny'') are sealed by Late Jurassic platform carbonates. The geological history conforms with that of the Inner Dinarides and adjoining areas; we therefore correlate the Mirdita-Pindos Ophiolite Zone with the Vardar Zone and explain its present position by far-distance west-directed thrusting.

Characterization of fecal indicator bacteria in sediments cores from the largest freshwater lake of Western Europe (Lake Geneva, Switzerland)

This study characterized the fecal indicator bacteria (FIB), including Escherichia coli (E. coli) and Enteroccocus (ENT), disseminated over time in the Bay of Vidy, which is the most contaminated area of Lake Geneva. Sediments were collected from a site located at similar to 500 m from the present waste water treatment plant (WWTP) outlet pipe, in front of the former WWTP outlet pipe, which was located at only 300 m from the coastal recreational area (before 2001). E. coil and ENT were enumerated in sediment suspension using the membrane filter method. The FIB characterization was performed for human Enterococcus faecalis (E. faecalis) and Enterococcus faecium (E. faecium) and human specific bacteroides by PCR using specific primers and a matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). Bacterial cultures revealed that maximum values of 35.2 x 10(8) and 6.6 x 10(6) CFU g(-1) dry sediment for E. coil and ENT, respectively, were found in the sediments deposited following eutrophication of Lake Geneva in the 1970s. whereas the WWTP started operating in 1964. The same tendency was observed for the presence of human fecal pollution: the percentage of PCR amplification with primers ESP-1/ESP-2 for E. faecalis and E. faecium indicated that more than 90% of these bacteria were from human origin. Interestingly, the PCR assays for specific-human bacteroides HF183/HF134 were positive for DNA extracted from all isolated strains of sediment surrounding WWPT outlet pipe discharge. The MALDI-TOF MS confirmed the presence of general E. coli and predominance E. faecium in isolated strains. Our results demonstrated that human fecal bacteria highly increased in the sediments contaminated with WWTP effluent following the eutrophication of Lake Geneva. Additionally, other FIB cultivable strains from animals or adapted environmental strains were detected in the sediment of the bay. The approaches used in this research are valuable to assess the temporal distribution and the source of the human fecal pollution in aquatic environments. (C) 2011 Elsevier Inc. All rights reserved.

The role of soil in vegetated gravelly river braid plains: more than just a passive response?

This paper reviews the role of alluvial soils in vegetated gravelly river braid plains. When considering decadal time scales of river evolution, we argue that it becomes vital to consider soil development as an emergent property of the developing ecosystem. Soil processes have been relatively overlooked in accounts of the interactions between braided river processes and vegetation, although soils have been observed on vegetated fluvial landforms. We hypothesise that soil development plays a major role in the transition (speed and pathway) from a fresh sediment deposit to a vegetated soil-covered landform. Disturbance (erosion and/or deposition), vertical sediment structure (process history), vegetation succession, biological activity and water table fluctuation are seen as the main controls on early alluvial soil evolution. Erosion and deposition processes may not only act as soil disturbing agents, but also as suppliers of ecosystem resources, because of their role in delivering and changing access (e.g. through avulsion) to fluxes of water, fine sediments and organic matter. In turn, the associated initial ecosystem may influence further fluvial landform development, such as through the trapping of fine-grained sediments (e.g. sand) by the engineering action of vegetation and the deposit stabilisation by the developing above and belowground biomass. This may create a strong feedback between geomorphological processes, vegetation succession and soil evolution which we summarise in a conceptual model. We illustrate this model by an example from the Allondon River (CH) and identify the research questions that follow.

Geomodelling of a fluvial system with semi-supervised support vector regression

Fluvial deposits are a challenge for modelling flow in sub-surface reservoirs. Connectivity and continuity of permeable bodies have a major impact on fluid flow in porous media. Contemporary object-based and multipoint statistics methods face a problem of robust representation of connected structures. An alternative approach to model petrophysical properties is based on machine learning algorithm ? Support Vector Regression (SVR). Semi-supervised SVR is able to establish spatial connectivity taking into account the prior knowledge on natural similarities. SVR as a learning algorithm is robust to noise and captures dependencies from all available data. Semi-supervised SVR applied to a synthetic fluvial reservoir demonstrated robust results, which are well matched to the flow performance

Shallow-platform palaeoenvironmental conditions recorded in deep-shelf sediments: C and stable isotopes in Upper Jurassic sections of southern Germany (Oxfordian-Kimmeridgian)

A high-resolution carbon and oxygen isotope analysis of Late Oxfordian-Early Kimmeridgian deep-shelf sediments of southern Germany is combined with investigation of nannofossil assemblage composition and sedimentological interpretations in order to evaluate the impact of regional palaeoenvironmental conditions on isotopic composition of carbonates. This study suggests that carbonate mud was essentially derived from the Jura shallow platform environments and also that the isotopic signature of carbonates deposited in the Swabian Alb deep shelf indirectly expresses the palaeoenvironmental evolution of the platform. Short-term fluctuations in delta(13) C and delta(18)O are probably controlled by changes in salinity (fresh-water input versus evaporation) in platform environments. Long-term fluctuations in carbon and oxygen isotope record throughout the Late Oxfordian-Early Kimmeridgian result from the interplay of increasing temperature and decreasing humidity, which both control the trophic level. Changes from mesotrophic to oligotrophic conditions in platform environments and in the deep-shelf surface waters are inferred. During the Late Oxfordian (Bimammatum Subzone to Planula Zone), the delta(13)C curve displays a positive shift of about 1 parts per thousand, which is comparable in intensity to global perturbations of the carbon cycle. This evident isotopic shift has not been documented yet in other basinal settings. It can be reasonably explained by local palaeoenvironmental changes on the Jura platform (salinity, temperature, and nutrient availability) that controlled platform carbonate production, and the geochemistry of overlying waters. However, increasing carbonate production on the Jura platform and related positive delta(13)C shifts recorded in the Swabian Alb deep shelf are the regional signatures of climatic changes affecting other palaeogeographical domains of Europe in which the carbonate production increased throughout the Late Oxfordian. (C) 2002 Elsevier Science B.V. All rights reserved.

Stable isotopes in lake Geneva carbonate sediments and molluscs: Review and new data

New isotopic results on bulk carbonate and mollusc (gastropods and bivalves) samples from Lake Geneva (Switzerland), spanning the period from the Oldest Dryas to the present day, are compared with pre-existing stable isotope data. According to preliminary calibration of modern samples, Lake Geneva endogenic calcite precipitates at or near oxygen isotopic equilibrium with ambient water, confirming the potential of this large lake to record paleoenvironmental and paleoclimatic changes. The onset of endogenic calcite precipitation at the beginning of the Allerod biozone is clearly indicated by the oxygen isotopic signature of bulk carbonate. A large change in delta(13)C values occurs during the Preboreal. This carbon shift is likely to be due to a change in bioproductivity and/or to a `'catchment effect'', the contribution of biogenic CO2 from the catchment area to the dissolved inorganic carbon reservoir of the lake water becoming significant only during the Preboreal. Gastropods are confirmed as valuable for studies of changes in paleotemperature and in paleowater isotopic composition, despite the presence of a vital effect. Mineralogical evidence indicates an increased detrital influence upon sedimentation since the Subboreal time period. On the other hand, stable isotope measurements of Subatlantic carbonate sediments show values comparable to those of pure endogenic calcite and of gastropods (taking into account the vital effect). This apparent disagreement still remains difficult to explain.