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.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

Intrathecal analgesia and restrictive perioperative fluid management within enhanced recovery pathway: hemodynamic implications.

BACKGROUND: Intrathecal analgesia and avoidance of perioperative fluid overload are key items within enhanced recovery pathways. Potential side effects include hypotension and renal dysfunction. STUDY DESIGN: From January 2010 until May 2010, all patients undergoing colorectal surgery within enhanced recovery pathways were included in this retrospective cohort study and were analyzed by intrathecal analgesia (IT) vs none (noIT). Primary outcomes measures were systolic and diastolic blood pressure, mean arterial pressure, and heart rate for 48 hours after surgery. Renal function was assessed by urine output and creatinine values. RESULTS: One hundred and sixty-three consecutive colorectal patients (127 IT and 36 noIT) were included in the analysis. Both patient groups showed low blood pressure values within the first 4 to 12 hours and a steady increase thereafter before return to baseline values after about 24 hours. Systolic and diastolic blood pressure and mean arterial pressure were significantly lower until 16 hours after surgery in patients having IT compared with the noIT group. Low urine output (<0.5 mL/kg/h) was reported in 11% vs 29% (IT vs noIT; p = 0.010) intraoperatively, 20% vs 11% (p = 0.387), 33% vs 22% (p = 0.304), and 31% vs 21% (p = 0.478) for postanesthesia care unit and postoperative days 1 and 2, respectively. Only 3 of 127 (2.4%) IT and 1 of 36 (2.8%) noIT patients had a transitory creatinine increase >50%; no patients required dialysis. CONCLUSIONS: Postoperative hypotension affects approximately 10% of patients within an enhanced recovery pathway and is slightly more pronounced in patients with IT. Hemodynamic depression persists for <20 hours after surgery; it has no measurable negative impact and therefore cannot justify detrimental postoperative fluid overload.

Prograde mica 40Ar/39Ar growth ages recorded in high pressure rocks (Syros, Cyclades, Greece)

Phengites from the eclogite and blueschist-facies sequences of the Cycladic island of Syros (Greece) have been dated by the in situ UV-laser ablation Ar-40/Ar-39 method. A massive, phengite-rich eclogite and an omphacite-rich metagabbro were investigated. The phengites are eubedral and coarse-grained (several 100 mum), strain-free and exhibit no evidence for late brittle deformation or recrystallization. Apparent ages in these samples range from 43 to 50 Ma for the phengite-rich eclogite and 42 to 52 Ma for the ompbacitic metagabbro. This large spread of ages is visible at all scales-within individual grains as well as in domains of several 100 mum and across the entire sample (ca. 2 cm). Such variations have been traditionally attributed to metamorphic cooling or the incorporation of excess argon. However, the textural equilibrium between the phengites and other high pressure phases and the subtle compositional variations within the phengites, especially the preservation of growth textures, alternatively suggest that the observed range in ages may reflect variations of radiogenic argon acquired during phengite formation and subsequent growth, thus dating a discrete event on the prograde path. This implies that the oldest phengite 40Ar/39Ar ages provide the best estimate of a minimum crystallization age, which is in agreement with recently reported U-Pb and Lu-Hf geochronological data. Our results are consistent with available stable isotope data and further suggest that, under fluid-restricted conditions, both stable and radiogenic isotopic systems can survive without significant isotopic exchange during subduction and exhumation from eclogite-facies P-T conditions. (C) 2004 Elsevier B.V. All rights reserved.

Pressure-dependent scaling scenarios in experiments of spontaneous imbibition

The scaling properties of the rough liquid-air interface formed in the spontaneous imbibition of a viscous liquid by a model porous medium are found to be very sensitive to the magnitude of the pressure difference applied at the liquid inlet. Interface fluctuations change from obeying intrinsic anomalous scaling at large negative pressure differences, to being super-rough with the same dynamic exponent z¿3 at less negative pressure differences, to finally obeying ordinary Family-Vicsek scaling with z¿2 at large positive pressure differences. This rich scenario reflects the relative importance on different length scales of capillary and permeability disorder, and the role of surface tension and viscous pressure in damping interface fluctuations.

Fluid dynamical prediction of changed v1 flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new "symmetrized" v1S(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Dissolution-reprecipitation of zircon at low-temperature, high-pressure conditions (Lanzo Massif, Italy)

An eclogite facies meta-plagiogranite from the Lanzo massif (western Alps, Italy) contains crystals of zircon intimately associated with allanite. Zircon displays different microtextures ranging from pristine, euhedral, and magmatic to fractured, porous varieties with mosaic zoning, and pervasive recrystallization into euhedral microcrystals. Fractures and voids in the recrystallized zircon microcrystals are mainly filled by high-pressure Na-rich pyroxene. Electron backscattered diffraction analysis revealed a similar crystallographic orientation for primary magmatic zircon crystals and microcrystals, with less than 2 degrees misorientation among neighboring microdomains. The textural change is coupled with chemical and isotopic modifications: recrystallized zircon domains contain significantly less Th and light- to mid-REE, but are richer in Sr than magmatic zircon crystals. Magmatic zircon preserves the protolith U-Pb age of 163.5 +/- 1.7 Ma, whereas zircon microcrystals have a mean age of 55 +/- 1 Ma. The coexisting allanite also contains inclusions of Na-rich pyroxene and has chemical features (elevated Sr and Ni contents and lack of Eu anomaly) indicating formation at high pressure. Despite being associated texturally with zircon, allanite yields a younger Th-Pb age of 46.5 +/- 3.0 Ma, suggesting that the Lanzo unit remained at relatively high pressure conditions for similar to 8 m.y. Zircon recrystallization proceeded with volume reduction and loss of material to an alkaline metamorphic fluid that acted as the agent for a coupled dissolution-reprecipitation process. Recrystallization occurred with minimum transport, in a low-strain environment, and was not significantly enhanced by metamictization. The source of the fluid for zircon recrystallization is most probably related to prograde devolatilization reactions in the surrounding serpentinite.

Interfacial instabilities of a fluid annulus in a rotating Hele-Shaw cell

We have studied the interfacial instabilities experienced by a liquid annulus as it moves radially in a circular Hele-Shaw cell rotating with angular velocity Omega. The instability of the leading interface (oil displacing air) is driven by the density difference in the presence of centrifugal forcing, while the instability of the trailing interface (air displacing oil) is driven by the large viscosity contrast. A linear stability analysis shows that the stability of the two interfaces is coupled through the pressure field already at a linear level. We have performed experiments in a dry cell and in a cell coated with a thin fluid layer on each plate, and found that the stability depends substantially on the wetting conditions at the leading interface. Our experimental results of the number of fingers resulting from the instability compare well with the predictions obtained through a numerical integration of the coupled equations derived from a linear stability analysis. Deep in the nonlinear regime we observe the emission of liquid droplets through the formation of thin filaments at the tip of outgrowing fingers.

Fluid evolution at the Variscan front in the vicinity of the Aachen thrust

Quartz-carbonate-chlorite veins were studied in borehole samples of the RWTH-1 well in Aachen. Veins formed in Devonian rocks in the footwall of the Aachen thrust during Variscan deformation and associated fluid flow. Primary fluid inclusions indicate subsolvus unmixing of a homogenous H(2)O-CO(2)-CH(4)-(N(2))-Na-(K)-Cl fluid into a H(2)O-Na-(K)-Cl solution and a vapour-rich CO(2)-(H(2)O, CH(4), N(2)) fluid. The aqueous end-member composition resembles that of metamorphic fluids of the Variscan front zone with salinities ranging from 4 to 7% NaCl equiv. and maximum homogenisation temperatures of close to 400A degrees C. Pressure estimates indicate a burial depth between 4,500 and 8,000 m at geothermal gradients between 50 and 75A degrees C/26 MPa, but pressure decrease to sublithostatic conditions is also indicated, probably as a consequence of fracture opening during episodic seismic activity. A second fluid system, mainly preserved in pseudo-secondary and secondary fluid inclusions, is characterised by fluid temperatures between 200 and 250A degrees C and salinities of < 5% NaCl equiv. Bulk stable isotope analyses of fluids released from vein quartz, calcite, and dolomite by decrepitation yielded delta D(H2O) values from -89 to -113 aEuro degrees, delta(13)C(CH4) from -26.9 to -28.9aEuro degrees (VPDB) and delta(13)C(CO2) from -12.8 to -23.3aEuro degrees (VPDB). The low delta D and delta(13)C range of the fluids is considered to be due to interaction with cracked hydrocarbons. The second fluid influx caused partial isotope exchange and disequilibrium. It is envisaged that an initial short lived flux of hot metamorphic fluids expelled from the epizonal metamorphic domains of the Stavelot-Venn massif. The metamorphic fluid was focused along major thrust faults of the Variscan front zone such as the Aachen thrust. A second fluid influx was introduced from formation waters in the footwall of the Aachen thrust as a consequence of progressive deformation. Mixing of the cooler and lower salinity formation water with the hot metamorphic fluid during episodic fluid trapping resulted in an evolving range of physicochemical fluid inclusion characteristics.

Transscleral Coulomb-controlled iontophoresis of methylprednisolone into the rabbit eye: influence of duration of treatment, current intensity and drug concentration on ocular tissue and fluid levels.

The major problems associated with the use of corticosteroids for the treatment of ocular diseases are their poor intraocular penetration to the posterior segment when administered locally and their secondary side effects when given systemically. To circumvent these problems more efficient methods and techniques of local delivery are being developed. The purposes of this study were: (1) to investigate the pharmacokinetics of intraocular penetration of hemisuccinate methyl prednisolone (HMP) after its delivery using the transscleral Coulomb controlled iontophoresis (CCI) system applied to the eye or after intravenous (i.v.) injection in the rabbit, (2) to test the safety of the CCI system for the treated eyes and (3) to compare the pharmacokinetic profiles of HMP intraocular distribution after CCI delivery to i.v. injection. For each parameter evaluated, six rabbit eyes were used. For the CCI system, two concentrations of HMP (62.5 and 150mg ml(-1)), various intensities of current and duration of treatment were analyzed. In rabbits serving as controls the HMP was infused in the CCI device but without applied electric current. For the i.v. delivery, HMP at 10mg kg(-1)as a 62.5mg ml(-1)solution was used. The rabbits were observed clinically for evidence of ocular toxicity. At various time points after the administration of drug, rabbits were killed and intraocular fluids and tissues were sampled for methylprednisolone (MP) concentrations by high pressure liquid chromatography (HPLC). Histology examinations were performed on six eyes of each group. Among groups that received CCI, the concentrations of MP increased in all ocular tissues and fluids in relation to the intensities of current used (0.4, 1.0 and 2.0mA/0.5cm(2)) and its duration (4 and 10min). Sustained and highest levels of MP were achieved in the choroid and the retina of rabbit eyes treated with the highest current and 10min duration of CCI. No clinical toxicity or histological lesions were observed following CCI. Negligible amounts of MP were found in ocular tissues in the CCI control group without application of current. Compared to i.v. administration, CCI achieved higher and more sustained tissue concentrations with negligible systemic absorption. These data demonstrate that high levels of MP can be safely achieved in intraocular tissues and fluids of the rabbit eye, using CCI. With this system, intraocular tissues levels of MP are higher than those achieved after i.v. injection. Furthermore, if needed, the drug levels achieved with CCI can be modulated as a function of current intensity and duration of treatment. CCI could therefore be used as an alternative method for the delivery of high levels of MP to the intraocular tissues of both the anterior and posterior segments.

A Variscan pressure-temperature-time path for the N-E Mont Blanc massif

The northeastern portion of the Mont Blanc massif in western Switzerland is predominantly comprised of the granitic rocks of the Mont Blanc intrusive suit, and the Mont Blanc basement gneisses. Within these metamorphic rocks are a variety of sub-economic Fe skarns. The mineral assemblages and fluid inclusions from these rocks have been used to derive age, pressure, temperature and fluid composition constraints for two Variscan events. Metamorphic hornblendes within the assemblages from the basement amphibolites and iron sk:lms have been dated using Ar-40/Ar-39, and indicate that these metamorphic events have a minimum age of approximately 334 Ma. Garnet-hornblende-plagioclase thermobarometry and stable isotope data obtained from the basement amphibolites are consistent with metamorphic temperatures in the range 515 to 580 degrees C, and pressures ranging from 5 to 8 kbar. Garnet-hornblende-magnetite thermobarometry and fluid inclusion studies indicate that the iron skarns formed at slightly lower temperatures, ranging from 400 to 500 degrees C in the presence of saline fluids at formational pressures similar to those experienced by the basement amphibolites. Late Paleozoic minimum uplift rates and geothermal gradients calculated using these data and the presence of Ladinien ichnofossils are on the order of 0.32 mm/year and 20 degrees C/km respectively. These uplift rates and geothermal gradients differ from those obtained from the neighbouring Aiguilles Rouges massif and indicate that these two massifs experienced different metamorphic conditions during the Carboniferous and Permian periods. During the early to late Carboniferous period the relative depths of the two massifs were reversed with the Aiguilles Rouges being initially unroofed at a much greater rate than the Mont Blanc, but experiencing relatively slower uplift rates near the termination of the Variscan orogeny.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Fluid dynamical prediction of changed v1-flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new"symmetrized" vS1(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Fluid flow and dispersion predictions in porous media by utilizing advanced numerical methods

Diplomityön tavoitteena oli tarkastella numeerisen virtauslaskennan avulla virtaukseen liittyviä ilmiöitä ja kaasun dispersiota. Diplomityön sisältö on jaettu viiteen osaan; johdantoon, teoriaan, katsaukseen virtauksen mallinnukseen huokoisessa materiaalissa liittyviin tutkimusselvityksiin, numeeriseen mallinnukseen sekä tulosten esittämiseen ja johtopäätöksiin. Diplomityön alussa kiinnitettiin huomiota erilaisiin kokeellisiin, numeerisiin ja teoreettisiin mallinnusmenetelmiin, joilla voidaan mallintaa virtausta huokoisessa materiaalissa. Kirjallisuusosassa tehtiin katsaus aikaisemmin julkaistuihin puoliempiirisiin ja empiirisiin tutkimusselvityksiin, jotka liittyvät huokoisen materiaalin aiheuttamaan painehäviöön. Numeerisessa virtauslaskenta osassa rakennettiin ja esitettiin huokoista materiaalia kuvaavat numeeriset mallit käyttäen kaupallista FLUENT -ohjelmistoa. Työn lopussa arvioitiin teorian, numeerisen virtauslaskennan ja kokeellisten tutkimusselvitysten tuloksia. Kolmiulotteisen huokoisen materiaalinnumeerisessa mallinnuksesta saadut tulokset vaikuttivat lupaavilta. Näiden tulosten perusteella tehtiin suosituksia ajatellen tulevaa virtauksen mallinnusta huokoisessa materiaalissa. Osa tässä diplomityössä esitetyistä tuloksista tullaan esittämään 55. Kanadan Kemiantekniikan konferenssissa Torontossa 1619 Lokakuussa 2005. ASME :n kansainvälisessä tekniikan alan julkaisussa. Työ on hyväksytty esitettäväksi esitettäväksi laskennallisen virtausmekaniikan (CFD) aihealueessa 'Peruskäsitteet'. Lisäksi työn yksityiskohtaiset tulokset tullaan lähettämään myös CES:n julkaisuun.

Primary fluid bolus therapy for infection-associated hypotension in the emergency department.

OBJECTIVES: The physiological changes associated with fluid bolus therapy (FBT) for patients with infection-associated hypotension in the emergency department (ED) are poorly understood. We describe the physiological outcomes of FBT in the first 6 hours (primary FBT) for patients presenting to the ED with infection-associated hypotension. METHODS: We studied 101 consecutive ED patients with infection and a systolic blood pressure (SBP)<100 mmHg who underwent FBT in the first 6 hours. RESULTS: We screened 1123 patients with infection and identified 101 eligible patients. The median primary FBT volume given was 1570 mL (interquartile range, 1000- 2490 mL). The average mean arterial pressure (MAP) did not change from admission to 6 hours in the whole cohort, or in patients who were hypotensive on arrival at the ED. However, the average MAP increased from its lowest value during the first 6 hours (66 mmHg [SD, 10 mmHg]) to its value at 6 hours (73 mmHg [SD, 12 mmHg]; P<0.001). The mean heart rate, body temperature, respiratory rate and plasma creatinine level decreased (P<0.05). In patients who were severely hypotensive (SBP<90 mmHg) on arrival at the ED, the MAP increased from 54 mmHg (SD, 8 mmHg) to 70 mmHg (SD, 14 mmHg) (P<0.001). At 6 hours, however, SBP was still <100 mmHg in 44 patients and <90 mmHg in 17 patients. When noradrenaline was used, in 10 patients, hypotension was corrected in all 10 and the MAP increased from 58 mmHg (SD, 9 mmHg) to 75 mmHg (SD, 13 mmHg). CONCLUSION: Among ED patients admitted to an Australian teaching hospital with infection, hypotension was uncommon. FBT for hypotension was limited in volumes given and failed to achieve a sustained SBP of >100 mmHg in 40% of cases. In contrast, noradrenaline therapy corrected hypotension in all patients who received it.