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.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Do seismic waves sense fracture connectivity?

A defining characteristic of fractured rocks is their very high level of seismic attenuation, which so far has been assumed to be mainly due to wave-induced fluid flow (WIFF) between the fractures and the pore space of the embedding matrix. Using oscillatory compressibility simulations based on the quasi-static poroelastic equations, we show that another important, and as of yet undocumented, manifestation of WIFF is at play in the presence of fracture connectivity. This additional energy loss is predominantly due to fluid flow within the connected fractures and is sensitive to their lengths, permeabilities, and intersection angles. Correspondingly, it contains key information on the governing hydraulic properties of fractured rock masses and hence should be accounted for whenever realistic seismic models of such media are needed.

New measurements of energy expenditure and physical activity in chronic kidney disease.

The accurate estimation of total daily energy expenditure (TEE) in chronic kidney patients is essential to allow the provision of nutritional requirements; however, it remains a challenge to collect actual physical activity and resting energy expenditure in maintenance dialysis patients. The direct measurement of TEE by direct calorimetry or doubly labeled water cannot be used easily so that, in clinical practice, TEE is usually estimated from resting energy expenditure and physical activity. Prediction equations may also be used to estimate resting energy expenditure; however, their use has been poorly documented in dialysis patients. Recently, a new system called SenseWear Armband (BodyMedia, Pittsburgh, PA) was developed to assess TEE, but so far no data have been published in chronic kidney disease patients. The aim of this review is to describe new measurements of energy expenditure and physical activity in chronic kidney disease patients.

Creation of an Age-Adjusted, Dual-Energy X-ray Absorptiometry-Derived Trabecular Bone Score Curve for the Lumbar Spine in Non-Hispanic US White Women.

The trabecular bone score (TBS, Med-Imaps, Pessac, France) is an index of bone microarchitecture texture extracted from anteroposterior dual-energy X-ray absorptiometry images of the spine. Previous studies have documented the ability of TBS of the spine to differentiate between women with and without fractures among age- and areal bone mineral density (aBMD)-matched controls, as well as to predict future fractures. In this cross-sectional analysis of data collected from 3 geographically dispersed facilities in the United States, we investigated age-related changes in the microarchitecture of lumbar vertebrae as assessed by TBS in a cohort of non-Hispanic US white American women. All subjects were 30 yr of age and older and had an L1-L4aBMDZ-score within ±2 SD of the population mean. Individuals were excluded if they had fractures, were on any osteoporosis treatment, or had any illness that would be expected to impact bone metabolism. All data were extracted from Prodigy dual-energy X-ray absorptiometry devices (GE-Lunar, Madison, WI). Cross-calibrations between the 3 participating centers were performed for TBS and aBMD. aBMD and TBS were evaluated for spine L1-L4 but also for all other possible vertebral combinations. To validate the cohort, a comparison between the aBMD normative data of our cohort and US non-Hispanic white Lunar data provided by the manufacturer was performed. A database of 619 non-Hispanic US white women, ages 30-90 yr, was created. aBMD normative data obtained from this cohort were not statistically different from the non-Hispanic US white Lunar normative data provided by the manufacturer (p = 0.30). This outcome thereby indirectly validates our cohort. TBS values at L1-L4 were weakly inversely correlated with body mass index (r = -0.17) and weight (r = -0.16) and not correlated with height. TBS values for all lumbar vertebral combinations decreased significantly with age. There was a linear decrease of 16.0% (-2.47 T-score) in TBS at L1-L4 between 45 and 90 yr of age (vs. -2.34 for aBMD). Microarchitectural loss rate increased after age 65 by 50% (-0.004 to -0.006). Similar results were obtained for other combinations of lumbar vertebra. TBS, an index of bone microarchitectural texture, decreases with advancing age in non-Hispanic US white women. Little change in TBS is observed between ages 30 and 45. Thereafter, a progressive decrease is observed with advancing age. The changes we observed in these American women are similar to that previously reported for a French population of white women (r(2) > 0.99). This reference database will facilitate the use of TBS to assess bone microarchitectural deterioration in clinical practice.

An interpretation of the deep seismic lines from the Penninic Alps of Valais (Switzerland)

The upper part of three deep seismic lines running across the Penninic Swiss Alps of Valais have been studied. Numerous reflectors illustrate the nappe structure of this internal part of the orogen. These reflectors, even at great depths (20-25 km), can be correlated with outcropping geological features and are most likely produced by lithological boundaries rather than by mylonites zones, which are hardly reflective in such an environment. Our interpretations, largely constrained by projections of the outcropping geology, have improved our knowledge of the deep structure of this segment of the Alpine belt, enhancing the importance of the backfolding and the crustal scale deformation phase which produced the Rawil-Valpelline depression and the Aar-Toce culmination. Furthermore we have here the possibility of correlating seismic patterns produced by ductile folds with the outcropping structures.

Enteral versus parenteral nutrition: comparison of energy metabolism in healthy subjects.

Continuous respiratory exchange measurements were performed on 10 healthy young women for 1 h before, 3 h during, and 3 h after either parenteral (iv) or intragastric (ig) administration of a nutrient mixture (52% glucose, 18% amino acid, and 30% lipid energy) infused at twice the postabsorptive resting energy expenditure (REE). REE rose from 0.98 +/- 0.02 (iv) and 0.99 +/- 0.02 kcal/min (ig) postabsorptively to 1.13 +/- 0.03 (iv) and 1.13 +/- 0.02 kcal/min (ig), resulting in nutrient-induced thermogenesis of 10 +/- 0.6 and 9.3 +/- 0.9%, respectively, when related to the metabolizable energy. The respiratory quotient rose from preinfusion values of 0.81 +/- 0.02 (iv) and 0.80 +/- 0.01 (ig) to 0.86 +/- 0.01 (iv) and 0.85 +/- 0.01 (ig). After nutrient administration the respiratory quotient fell significantly to below the preinfusion values. Plasma glucose and insulin concentrations rose during nutrient administration but were higher during the intravenous route. It is concluded that, although the response time to intragastric administration was delayed, the thermic effects and overall substrate oxidations were comparable during intravenous or intragastric administration, albeit, at lower plasma glucose and insulin concentrations via the intragastric route.

Macronutrients and energy balance in obesity.

Energy balance is the difference between metabolizable energy intake and total energy expenditure. Energy intake is difficult to measure accurately; changes in body weight, for example, are not a good measure of the adequacy of energy intake, because fluctuations in body weight are common even if the overall trend is toward weight loss. It is now customary to assess energy requirements indirectly from total energy expenditure. Total energy expenditure consists of basal metabolism, postprandial thermogenesis, and physical activity. Energy expenditure is related to both body weight and body composition. A reduction in total energy expenditure accompanies weight loss, because basal metabolic rate decreases with the loss of lean tissue mass. Similarly, with weight gain, there is an increase in basal metabolic rate, because lean tissue mass grows to support the increase in fat tissue mass. Excess energy intake over energy expenditure causes weight gain and an accompanying increase in total energy expenditure. Following a period of adaptation, total energy expenditure will match energy intake and body weight will stabilize at a higher level. This same relationship holds for weight loss. Respiratory quotient (measured in steady state) is an indication of the proportion of energy expenditure derived from fat and carbohydrate oxidation. Over long periods of time, fat balance is equivalent to energy balance, as an excess of fat intake over fat oxidation causes fat storage.

Social isolation causers mortality by disrupting energy homeostrasis in ants

Social deprivation can have negative effects on the lives of social animals, including humans, yet little is known about the mechanisms by which social withdrawal affects animal health. Here we show that in the carpenter ant Camponotus fellah, socially isolated workers have a greatly reduced life span relative to ants kept in groups of ten individuals. By using a new tracking system, we found that social isolation resulted in important behavioral changes and greatly increased locomotor activity. The higher activity of single ants and their increased propensity to leave the nest to move along the walls suggested that the increased mortality of isolated ants might stem from an imbalance of energy income and expenditure. This view was supported by the finding that while isolated ants ingested the same amount of food as grouped ants, they retained food in the crop, hence preventing its use as an energy source. Moreover, the difference in life span between single and grouped individuals vanished when ants were not fed. This study thus underlines the role of social interactions as key regulators of energy balance, which ultimately affects aging and health in a highly social organism.

Effects of a glucose meal on energy metabolism in patients with cirrhosis before and after liver transplantation.

HYPOTHESIS: Liver transplantation results in hepatic denervation. This may produce alterations of liver energy and substrate metabolism, which may contribute to weight gain after liver transplantation. DESIGN: Prospective clinical study. SETTING: Liver transplantation clinics in a university hospital. PATIENTS: Seven nondiabetic patients with cirrhosis were recruited while on a waiting list for liver transplantation. Seven healthy subjects were recruited as controls. INTERVENTION: Orthotopic liver transplantation. MAIN OUTCOME MEASURES: Evaluation of energy and substrate metabolism after ingestion of a glucose load with indirect calorimetry was performed before, 2 to 6 weeks after, and 5 to 19 months after transplantation. Whole-body glucose oxidation and storage and glucose-induced thermogenesis were calculated. RESULTS: Patients with cirrhosis had modestly elevated resting energy expenditure and normal glucose-induced thermogenesis and postprandial glucose oxidation and storage. These measures remained unchanged after liver transplantation despite a significant increase in postprandial glycemia. Patients, however, gained an average of 3 kg of body weight after 5 to 19 months compared with their weight before transplantation. CONCLUSION: Liver denervation secondary to transplantation does not lead to alterations of energy metabolism after ingestion of a glucose load.