Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

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.

Long-term citalopram administration reduces responsiveness of HPA axis in patients with major depression: relationship with S-citalopram concentrations in plasma and cerebrospinal fluid (CSF) and clinical response.

RATIONALE: A dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis is a well-documented neurobiological finding in major depression. Moreover, clinically effective therapy with antidepressant drugs may normalize the HPA axis activity. OBJECTIVE: The aim of this study was to test whether citalopram (R/S-CIT) affects the function of the HPA axis in patients with major depression (DSM IV). METHODS: Twenty depressed patients (11 women and 9 men) were challenged with a combined dexamethasone (DEX) suppression and corticotropin-releasing hormone (CRH) stimulation test (DEX/CRH test) following a placebo week and after 2, 4, and 16 weeks of 40 mg/day R/S-CIT treatment. RESULTS: The results show a time-dependent reduction of adrenocorticotrophic hormone (ACTH) and cortisol response during the DEX/CRH test both in treatment responders and nonresponders within 16 weeks. There was a significant relationship between post-DEX baseline cortisol levels (measured before administration of CRH) and severity of depression at pretreatment baseline. Multiple linear regression analyses were performed to identify the impact of psychopathology and hormonal stress responsiveness and R/S-CIT concentrations in plasma and cerebrospinal fluid (CSF). The magnitude of decrease in cortisol responsivity from pretreatment baseline to week 4 on drug [delta-area under the curve (AUC) cortisol] was a significant predictor (p<0.0001) of the degree of symptom improvement following 16 weeks on drug (i.e., decrease in HAM-D21 total score). The model demonstrated that the interaction of CSF S-CIT concentrations and clinical improvement was the most powerful predictor of AUC cortisol responsiveness. CONCLUSION: The present study shows that decreased AUC cortisol was highly associated with S-CIT concentrations in plasma and CSF. Therefore, our data suggest that the CSF or plasma S-CIT concentrations rather than the R/S-CIT dose should be considered as an indicator of the selective serotonergic reuptake inhibitors (SSRIs) effect on HPA axis responsiveness as measured by AUC cortisol response.

Ultra high performance supercritical fluid chromatography coupled with tandem mass spectrometry for screening of doping agents. II: Analysis of biological samples.

The potential and applicability of UHPSFC-MS/MS for anti-doping screening in urine samples were tested for the first time. For this purpose, a group of 110 doping agents with diverse physicochemical properties was analyzed using two separation techniques, namely UHPLC-MS/MS and UHPSFC-MS/MS in both ESI+ and ESI- modes. The two approaches were compared in terms of selectivity, sensitivity, linearity and matrix effects. As expected, very diverse retentions and selectivities were obtained in UHPLC and UHPSFC, proving a good complementarity of these analytical strategies. In both conditions, acceptable peak shapes and MS detection capabilities were obtained within 7min analysis time, enabling the application of these two methods for screening purposes. Method sensitivity was found comparable for 46% of tested compounds, while higher sensitivity was observed for 21% of tested compounds in UHPLC-MS/MS and for 32% in UHPSFC-MS/MS. The latter demonstrated a lower susceptibility to matrix effects, which were mostly observed as signal suppression. In the case of UHPLC-MS/MS, more serious matrix effects were observed, leading typically to signal enhancement and the matrix effect was also concentration dependent, i.e., more significant matrix effects occurred at the lowest concentrations.

Bradykinin in blood and cerebrospinal fluid after acute cerebral lesions: correlations with cerebral edema and intracranial pressure.

Abstract Bradykinin (BK) was shown to stimulate the production of physiologically active metabolites, blood-brain barrier disruption, and brain edema. The aim of this prospective study was to measure BK concentrations in blood and cerebrospinal fluid (CSF) of patients with traumatic brain injury (TBI), subarachnoid hemorrhage (SAH), intracerebral hemorrhage (ICH), and ischemic stroke and to correlate BK levels with the extent of cerebral edema and intracranial pressure (ICP). Blood and CSF samples of 29 patients suffering from acute cerebral lesions (TBI, 7; SAH,: 10; ICH, 8; ischemic stroke, 4) were collected for up to 8 days after insult. Seven patients with lumbar drainage were used as controls. Edema (5-point scale), ICP, and the GCS (Glasgow Coma Score) at the time of sample withdrawal were correlated with BK concentrations. Though all plasma-BK samples were not significantly elevated, CSF-BK levels of all patients were significantly elevated in overall (n=73) and early (≤72 h) measurements (n=55; 4.3±6.9 and 5.6±8.9 fmol/mL), compared to 1.2±0.7 fmol/mL of controls (p=0.05 and 0.006). Within 72 h after ictus, patients suffering from TBI (p=0.01), ICH (p=0.001), and ischemic stroke (p=0.02) showed significant increases. CSF-BK concentrations correlated with extent of edema formation (r=0.53; p<0.001) and with ICP (r=0.49; p<0.001). Our results demonstrate that acute cerebral lesions are associated with increased CSF-BK levels. Especially after TBI, subarachnoid and intracerebral hemorrhage CSF-BK levels correlate with extent of edema evolution and ICP. BK-blocking agents may turn out to be effective remedies in brain injuries.

Low-Dose Vascular Photodynamic Therapy Decreases Tumor Interstitial Fluid Pressure, which Promotes Liposomal Doxorubicin Distribution in a Murine Sarcoma Metastasis Model.

INTRODUCTION: Solid tumors are known to have an abnormal vasculature that limits the distribution of chemotherapy. We have recently shown that tumor vessel modulation by low-dose photodynamic therapy (L-PDT) could improve the uptake of macromolecular chemotherapeutic agents such as liposomal doxorubicin (Liporubicin) administered subsequently. However, how this occurs is unknown. Convection, the main mechanism for drug transport between the intravascular and extravascular spaces, is mostly related to interstitial fluid pressure (IFP) and tumor blood flow (TBF). Here, we determined the changes of tumor and surrounding lung IFP and TBF before, during, and after vascular L-PDT. We also evaluated the effect of these changes on the distribution of Liporubicin administered intravenously (IV) in a lung sarcoma metastasis model. MATERIALS AND METHODS: A syngeneic methylcholanthrene-induced sarcoma cell line was implanted subpleurally in the lung of Fischer rats. Tumor/surrounding lung IFP and TBF changes induced by L-PDT were determined using the wick-in-needle technique and laser Doppler flowmetry, respectively. The spatial distribution of Liporubicin in tumor and lung tissues following IV drug administration was then assessed in L-PDT-pretreated animals and controls (no L-PDT) by epifluorescence microscopy. RESULTS: L-PDT significantly decreased tumor but not lung IFP compared to controls (no L-PDT) without affecting TBF. These conditions were associated with a significant improvement in Liporubicin distribution in tumor tissues compared to controls (P < .05). DISCUSSION: L-PDT specifically enhanced convection in blood vessels of tumor but not of normal lung tissue, which was associated with a significant improvement of Liporubicin distribution in tumors compared to controls.

Association of N-acetylaspartate and cerebrospinal fluid Aβ42 in dementia.

The interplay of amyloid and mitochondrial function is considered crucial in the pathophysiology of Alzheimer's disease (AD). We tested the association of the putative marker of mitochondrial function N-acetylaspartate (NAA) as measured by proton magnetic resonance spectroscopy within the medial temporal lobe and cerebrospinal fluid amyoid-β42 (Aβ42), total Tau and pTau181. 109 patients were recruited in a multicenter study (40 mild AD patients, 14 non-AD dementia patients, 29 mild cognitive impairment (MCI) AD-type patients, 26 MCI of non-AD type patients). NAA correlated with Aβ42 within the AD group. Since the NAA concentration is coupled to neuronal mitochondrial function, the correlation between NAA and Aβ42 may reflect the interaction between disrupted mitochondrial pathways and amyloid production.

Optimized venous return with a self-expanding cannula: from computational fluid dynamics to clinical application.

The Smart canula concept allows for collapsed cannula insertion, and self-expansion within a vein of the body. (A) Computational fluid dynamics, and (B) bovine experiments (76+/-3.8 kg) were performed for comparative analyses, prior to (C) the first clinical application. For an 18F access, a given flow of 4 l/min (A) resulted in a pressure drop of 49 mmHg for smart cannula versus 140 mmHg for control. The corresponding Reynolds numbers are 680 versus 1170, respectively. (B) For an access of 28F, the maximal flow for smart cannula was 5.8+/-0.5 l/min versus 4.0+/-0.1 l/min for standard (P<0.0001), for 24F 5.5+/-0.6 l/min versus 3.2+/-0.4 l/min (P<0.0001), and for 20F 4.1+/-0.3 l/min versus 1.6+/-0.3 l/min (P<0.0001). The flow obtained with the smart cannula was 270+/-45% (20F), 172+/-26% (24F), and 134+/-13% (28F) of standard (one-way ANOVA, P=0.014). (C) First clinical application (1.42 m2) with a smart cannula showed 3.55 l/min (100% predicted) without additional fluids. All three assessment steps confirm the superior performance of the smart cannula design.

Immuno-reactive proteins from Mycobacterium immunogenum useful for serodiagnosis of metalworking fluid hypersensitivity pneumonitis.

Metalworking fluid-associated hypersensitivity pneumonitis (MWF-HP) is a pulmonary disease caused by inhaling microorganisms present in the metalworking fluids used in the industrial sector. Mycobacterium immunogenum is the main etiological agent. Among the clinical, radiological and biological tools used for diagnosis, serological tests are important. The aim of this study was to identify immunogenic proteins in M. immunogenum and to use recombinant antigens for serological diagnosis of MWF-HP. Immunogenic proteins were detected by two-dimensional Western blot and candidate proteins were identified by mass spectrometry. Recombinant antigens were expressed in Escherichia coli and tested by enzyme-linked immunosorbent assay (ELISA) with the sera of 14 subjects with MWF-HP and 12 asymptomatic controls exposed to M. immunogenum. From the 350 spots visualized by two-dimensional gel electrophoresis with M. immunogenum extract, 6 immunogenic proteins were selected to be expressed as recombinant antigens. Acyl-CoA dehydrogenase antigen allowed for the best discrimination of MWF-HP cases against controls with an area under the receiver operating characteristics (ROC) curve of 0.930 (95% CI=0.820-1), a sensitivity of 100% and a specificity of 83% for the optimum threshold. Other recombinant antigens correspond to acyl-CoA dehydrogenase FadE, cytosol aminopeptidase, dihydrolipoyl dehydrogenase, serine hydroxymethyltransferase and superoxide dismutase. This is the first time that recombinant antigens have been used for the serodiagnosis of hypersensitivity pneumonitis. The availability of recombinant antigens makes it possible to develop standardized serological tests which in turn could simplify diagnosis, thus making it less invasive.

Characteristics of infections associated with external ventricular drains of cerebrospinal fluid.

OBJECTIVES: Manifestations of external ventricular drain (EVD) - associated infections overlap with those of the underlying neurosurgical conditions. We analyzed characteristics of EVD-associated infections. METHODS: We included patients aged ≥18 years with EVD-associated infections from 1997 to 2008, using modified CDC criteria for nosocomial infections. Hospital charts were reviewed retrospectively and the in-hospital outcome was evaluated. RESULTS: Forty-eight patients with EVD-associated infections were included (median age, 52 years, range 20-74 years). The median EVD-indwelling time was 7 days (range, 1-39 days) and EVD-associated infection occurred 6 days after insertion (range, 1-17 days). In 23% of patients, meningitis occurred 1-10 days after EVD removal. Fever >38 °C was present in 79% of patients, but Glasgow Coma Scale (GCS) scores were reduced in only 29%, and headache, vomiting and/or neck stiffness were present in only 31%. The median cerebrospinal fluid (CSF) leukocyte count was higher at onset of EVD-associated infection than at EVD insertion (175 × 10(6)/l versus 46 × 10(6)/l, p = 0.021), but other CSF parameters did not differ significantly. The most commonly implicated organisms were coagulase-negative staphylococci (63%) and Propionibacterium acnes (15%). CONCLUSIONS: Fever and increased CSF leukocytes should raise the suspicion of EVD-associated infection, which may occur up to 10 days after removal of EVD.

Cerebrospinal fluid antimicroglial antibodies in Alzheimer disease: a putative marker of an ongoing inflammatory process.

Immunocompetent microglia play an important role in the pathogenesis of Alzheimer's disease (AD). Antimicroglial antibodies in the cerebrospinal fluid (CSF) in clinically diagnosed AD patients have been previously recorded. Here, we report the results of the analysis of the CSF from 38 autopsy cases: 7 with definite AD; 14 with mild and 10 with moderate Alzheimer's type pathology; and 7 controls. Antimicroglial antibodies were identified in 70% of patients with definite AD, in 80% of patients with moderate and in 28% of patients with mild Alzheimer's type pathology. CSF antimicroglial antibodies were not observed in any of the control cases. The results show that CSF antimicroglial antibodies are present in the majority of patients with definite AD and also in cases with moderate Alzheimer's type changes. They may also indicate dysregulation of microglial function. Together with previous observations, these findings indicate that compromised immune defense mechanisms play an important role in the pathogenesis of AD.