Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Blood flow, flow reserve, and glucose utilization in viable and nonviable myocardium in patients with ischemic cardiomyopathy.

PURPOSE: The aim of the study was to determine whether glucose uptake in viable myocardium of ischemic cardiomyopathy patients depends on rest myocardial blood flow (MBF) and the residual myocardial flow reserve (MFR). METHODS: Thirty-six patients with ischemic cardiomyopathy (left ventricular ejection fraction 25 ± 10 %) were studied with N-ammonia and F-fluorodeoxyglucose (FDG) positron emission tomography (PET). Twenty age-matched normals served as controls. Regional MBF was determined at rest and during dipyridamole hyperemia and regional FDG extraction was estimated from regional FDG to N-ammonia activity ratios. RESULTS: Rest MBF was reduced in viable (0.42 ± 0.18 ml/min per g) and nonviable regions (0.32 ± 0.09 ml/min per g) relative to remote regions (0.68 ± 0.23 ml/min per g, p < 0.001) and to normals (0.63 ± 0.13 ml/min per g). Dipyridamole raised MBFs in controls, remote, viable, and nonviable regions. MBFs at rest (p < 0.05) and stress (p < 0.05) in viable regions were significantly higher than that in nonviable regions, while MFRs did not differ significantly (p > 0.05). Compared to MFR in remote myocardium, MFRs in viable regions were similar (1.39 ± 0.56 vs 1.70 ± 0.45, p > 0.05) but were significantly lower in nonviable regions (1.23 ± 0.43, p < 0.001). Moreover, the FDG and thus glucose extraction was higher in viable than in remote (1.40 ± 0.14 vs 0.90 ± 0.20, p < 0.001) and in nonviable regions (1.13 ± 0.21, p < 0.001). The extraction of FDG in viable regions was independent of rest MBF but correlated inversely with MFRs (r =-0.424, p < 0.05). No correlation between the FDG extraction and MFR was observed in nonviable regions. CONCLUSION: As in the animal model, decreasing MFRs in viable myocardium are associated with increasing glucose extraction that likely reflects a metabolic adaptation of remodeling hibernating myocytes.

Plant DNA barcodes and the influence of gene flow.

Success of species assignment using DNA barcodes has been shown to vary among plant lineages because of a wide range of different factors. In this study, we confirm the theoretical prediction that gene flow influences species assignment with simulations and a literature survey. We show that the genome experiencing the highest gene flow is, in the majority of the cases, the best suited for species delimitation. Our results clearly suggest that, for most angiosperm groups, plastid markers will not be the most appropriate for use as DNA barcodes. We therefore advocate shifting the focus from plastid to nuclear markers to achieve an overall higher success using DNA barcodes.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Added prognostic value of myocardial blood flow quantitation in rubidium-82 positron emission tomography imaging.

AIMS: We studied the respective added value of the quantitative myocardial blood flow (MBF) and the myocardial flow reserve (MFR) as assessed with (82)Rb positron emission tomography (PET)/CT in predicting major adverse cardiovascular events (MACEs) in patients with suspected myocardial ischaemia. METHODS AND RESULTS: Myocardial perfusion images were analysed semi-quantitatively (SDS, summed difference score) and quantitatively (MBF, MFR) in 351 patients. Follow-up was completed in 335 patients and annualized MACE (cardiac death, myocardial infarction, revascularization, or hospitalization for congestive heart failure or de novo stable angor) rates were analysed with the Kaplan-Meier method in 318 patients after excluding 17 patients with early revascularizations (<60 days). Independent predictors of MACEs were identified by multivariate analysis. During a median follow-up of 624 days (inter-quartile range 540-697), 35 MACEs occurred. An annualized MACE rate was higher in patients with ischaemia (SDS >2) (n = 105) than those without [14% (95% CI = 9.1-22%) vs. 4.5% (2.7-7.4%), P < 0.0001]. The lowest MFR tertile group (MFR <1.8) had the highest MACE rate [16% (11-25%) vs. 2.9% (1.2-7.0%) and 4.3% (2.1-9.0%), P < 0.0001]. Similarly, the lowest stress MBF tertile group (MBF <1.8 mL/min/g) had the highest MACE rate [14% (9.2-22%) vs. 7.3% (4.2-13%) and 1.8% (0.6-5.5%), P = 0.0005]. Quantitation with stress MBF or MFR had a significant independent prognostic power in addition to semi-quantitative findings. The largest added value was conferred by combining stress MBF to SDS. This holds true even for patients without ischaemia. CONCLUSION: Perfusion findings in (82)Rb PET/CT are strong MACE outcome predictors. MBF quantification has an added value allowing further risk stratification in patients with normal and abnormal perfusion images.

Reduced atrial emptying after orthotopic heart transplantation masquerading as restrictive transmitral Doppler flow pattern?

BACKGROUND: An elevated early (E) to late (A) diastolic filling velocities ratio, typically seen in advanced diastolic dysfunction, has also been observed after cardioversion of atrial fibrillation as a consequence of the depressed left atrial (LA) contractility. We hypothesized that the impaired LA contractile function demonstrated after orthotopic cardiac transplantation (OCT) could also lead to this "pseudorestrictive" pattern. METHOD: E/A ratio related to the tissue Doppler early mitral annular velocity (Ea) as preload-independent index of LV relaxation was evaluated in all consecutive OCT patients between 2005 and 2007. RESULTS: The study population comprised 48 patients 97 ± 77 months after OCT. Thirty-two patients (67%) had an E/A ratio > 2. LV systolic function and myocardial relaxation assessed by the Ea velocity were similar compared to patients with normal ratio (61 ± 6% vs. 60 ± 12%, P = 0.854 and 15 ± 4 cm/s vs. 14 ± 3 cm/s, r = 0.15, P = 0.323, respectively). On the other hand, the proportion of the recipient and donor LA cuffs as estimated by the recipient/global LA area ratio and the LA emptying fraction significantly correlated with the E/A ratio (r = 0.40, P = 0.005 and r =-0.33, P = 0.022, respectively). CONCLUSION: Our study shows that there is a high prevalence of elevated E/A ratio after standard OCT which seems mainly related to reduced LA contractility. Recognition of this "pseudorestrictive" pattern may avoid misdiagnosis of diastolic dysfunction.

Calcium reabsorption in the distal tubule: regulation by sodium, pH, and flow.

We developed a mathematical model of Ca transport along the late distal convoluted tubule (DCT2) and the connecting tubule (CNT) to investigate the mechanisms that regulate Ca reabsorption in the DCT2-CNT. The model accounts for apical Ca influx across transient receptor potential vanilloid 5 (TRPV5) channels and basolateral Ca efflux via plasma membrane Ca-ATPase pumps and type 1 Na/Ca exchangers (NCX1). Model simulations reproduce experimentally observed variations in Ca uptake as a function of extracellular pH, Na, and Mg concentration. Our results indicate that amiloride enhances Ca reabsorption in the DCT2-CNT predominantly by increasing the driving force across NCX1, thereby stimulating Ca efflux. They also suggest that because aldosterone upregulates both apical and basolateral Na transport pathways, it has a lesser impact on Ca reabsorption than amiloride. Conversely, the model predicts that full NCX1 inhibition and parathyroidectomy each augment the Ca load delivered to the collecting duct severalfold. In addition, our results suggest that regulation of TRPV5 activity by luminal pH has a small impact, per se, on transepithelial Ca fluxes; the reduction in Ca reabsorption induced by metabolic acidosis likely stems from decreases in TRPV5 expression. In contrast, elevations in luminal Ca are predicted to significantly decrease TRPV5 activity via the Ca-sensing receptor. Nevertheless, following the administration of furosemide, the calcium-sensing receptor-mediated increase in Ca reabsorption in the DCT2-CNT is calculated to be insufficient to prevent hypercalciuria. Altogether, our model predicts complex interactions between calcium and sodium reabsorption in the DCT2-CNT.

Flow cytofluorometric analysis of the binding of Vicia villosa lectin to T lymphoblasts: lack of correlation with cytolytic function.

The relationship between the binding of Vicia villosa (VV) lectin and the expression of cytolytic function in T lymphoblasts has been investigated using flow cytofluorometric techniques. Spleen cells activated in vitro in 5-day mixed leukocyte cultures (MLC) were incubated sequentially with VV, rabbit anti-V antiserum, and fluoresceinated sheep anti-rabbit IgG. When these stained MLC cells were passed on a flow cytometer gated to exclude nonviable cells and small lymphocytes, a single heterogeneous peak of fluorescence was seen, as compared to control MLC cells that had not been incubated with VV. Fluorescence of lymphoblasts was dependent upon lectin dose and was eliminated when staining was performed in the presence of N-acetyl-D-galactosamine, the appropriate competitive sugar for VV. T cell blast populations activated against H-2, Mls, or parasite antigens all had comparable levels of fluorescence after staining with VV, although the cytolytic activity of these cells varied widely. Furthermore, when MLC lymphoblasts binding large or small amounts of VV were sorted on the basis of their relative fluorescence intensity and tested for cytolytic function, no appreciable difference in activity between the 2 populations was observed. These results are inconsistent with the hypothesis that VV binds selectively to cytolytic T lymphocytes.

Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.

Unilateral urinary flow impairment at the pelviureteral junction: outcome of renal function with respect to therapeutic strategy.

OBJECTIVES: To evaluate the renal function outcome in children with unilateral hydronephrosis and urinary flow impairment at the pelviureteral junction with respect to the therapeutic strategy. METHODS: We retrospectively selected 45 children with iodine-123-hippuran renography performed at diagnosis and after 3 or more years of follow-up. All children had bilateral nonobstructive pattern findings on diuretic renography at follow-up. Eleven children were treated conservatively, and 34 underwent unilateral pyeloplasty. Split and individual renal function, measured by an accumulation index, was computed from background-corrected renograms for the affected and contralateral kidneys at diagnosis and the follow-up examination. RESULTS: Of 11 children treated conservatively, 9 had normal bilateral function at diagnosis, all had reached normal function at follow-up. Of the 34 operated kidneys, 12 (38%) had initially normal function that remained normal at the follow-up examination, and 22 had impaired function that had normalized at the follow-up examination in 15 (68%). The function of the contralateral kidneys was increased in 5 of 8 children with persistently abnormal affected kidneys. Pyeloplasty was performed in 23 children (68%) and 11 children (32%) younger and older than 1 year, respectively. The function of the affected kidneys increased in both groups, but normalization occurred only in the younger children. CONCLUSIONS: Of the children selected for conservative treatment, 82% had normal bilateral renal function at diagnosis that was normal in all at the follow-up examination. Of the children treated surgically, 65% had initially impaired function of the affected kidney that improved in 87% after pyeloplasty. Normalization of function was observed only in children who were younger than 1 year old at surgery. Persistently low function of the affected kidney was compensated for by the contralateral one regardless of the age at surgery.

Topographic forcing of flow partition and flow structures at river bifurcations

This paper presents the predicted flow dynamics from the application of a Reynolds-averaged NavierStokes model to a series of bifurcation geometries with morphologies measured during previous flume experiments. The topography of the bifurcations consists of either plane or bedform-dominated beds which may or may not possess discordance between the two bifurcation distributaries. Numerical predictions are compared with experimental results to assess the ability of the numerical model to reproduce the division of flow into the bifurcation distributaries. The hydrodynamic model predicts: (1) diverting fluxes in the upstream channel which direct water into the distributaries; (2) super-elevation of the free surface induced at the bifurcation edge by pressure differences; and (3) counter-rotating secondary circulation cells which develop upstream of the apex of the bifurcation and move into the downstream channels, with water converging at the surface and diverging at the bed. When bedforms are not present, weak transversal fluxes characterize the upstream channel for almost its entire length, associated with clearly distinguishable secondary circulation cells, although these may be under-estimated by the turbulence model used in the solution. In the bedform dominated case, the same hydrodynamic conditions were not observed, with the bifurcation influence restricted and depth scale secondary circulation cells not forming. The results also demonstrate the dominant effect bed discordance has upon flow division between the two distributaries. Finally, results indicate that in bedform dominated rivers. Consequently, we suggest that sand-bed river bifurcations are more likely to have an influence that extends much further upstream and have a greater impact upon water distribution. This may contribute to observed morphological differences between sand-bedded and gravel-bedded braided river networks. Copyright (C) 2012 John Wiley & Sons, Ltd.

Seismic attenuation due to wave-induced fluid flow at microscopic and macroscopic scales

Wave-induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro-cracks or broken grain contacts, causes the stiffness of the so-called modified dry frame to be complex-valued and frequency-dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave-induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.