Uncertainty Quantification Of A Semi-Supervised Support Vector Regression Reservoir Model

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. Recent advances in machine learning offer a novel approach to model spatial distribution of petrophysical properties in complex reservoirs alternative to geostatistics. The approach is based of semisupervised learning, which handles both ?labelled? observed data and ?unlabelled? data, which have no measured value but describe prior knowledge and other relevant data in forms of manifolds in the input space where the modelled property is continuous. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic geological features and describe stochastic variability and non-uniqueness of spatial properties. On the other hand, it is able to capture and preserve key spatial dependencies such as connectivity of high permeability geo-bodies, which is often difficult in contemporary petroleum reservoir studies. Semi-supervised SVR as a data driven algorithm is designed to integrate various kind of conditioning information and learn dependences from it. The semi-supervised SVR model is able to balance signal/noise levels and control the prior belief in available data. In this work, stochastic semi-supervised SVR geomodel is integrated into Bayesian framework to quantify uncertainty of reservoir production with multiple models fitted to past dynamic observations (production history). Multiple history matched models are obtained using stochastic sampling and/or MCMC-based inference algorithms, which evaluate posterior probability distribution. Uncertainty of the model is described by posterior probability of the model parameters that represent key geological properties: spatial correlation size, continuity strength, smoothness/variability of spatial property distribution. The developed approach is illustrated with a fluvial reservoir case. The resulting probabilistic production forecasts are described by uncertainty envelopes. The paper compares the performance of the models with different combinations of unknown parameters and discusses sensitivity issues.

Deer as a potential wildlife reservoir for Parachlamydia species.

Wildlife populations represent an important reservoir for emerging pathogens and trans-boundary livestock diseases. However, detailed information relating to the occurrence of endemic pathogens such as those of the order Chlamydiales in such populations is lacking. During the hunting season of 2008, 863 samples (including blood, conjunctival swabs, internal organs and faeces) were collected in the Eastern Swiss Alps from 99 free-living red deer (Cervus elaphus) and 64 free-living roe deer (Capreolus capreolus) and tested using ELISA, PCR and immunohistochemistry for members of the family Chlamydiaceae and the genus Parachlamydia. Parachlamydia spp. were detected in the conjunctival swabs, faeces and internal organs of both species of deer (2.4% positive, with a further 29.5% inconclusive). The very low occurrence of Chlamydiaceae (2.5%) was in line with serological data (0.7% seroprevalence for Chlamydia abortus). Further investigations are required to elucidate the zoonotic potential, pathogenicity, and distribution of Parachlamydia spp. in wild ruminants.

Insights into low fish mercury bioaccumulation in a mercury-contaminated reservoir, Guizhou, China.

We examined Hg biogeochemistry in Baihua Reservoir, a system affected by industrial wastewater containing mercury (Hg). As expected, we found high levels of total Hg (THg, 664-7421 ng g(-1)) and monomethylmercury (MMHg, 3-21 ng g(-1)) in the surface sediments (0-10 cm). In the water column, both THg and MMHg showed strong vertical variations with higher concentrations in the anoxic layer (>4m) than in the oxic layer (0-4 m), which was most pronounced for the dissolved MMHg (p < 0.001). However, mercury levels in biota samples (mostly cyprinid fish) were one order of magnitude lower than common regulatory values (i.e. 0.3-0.5 mg kg(-1)) for human consumption. We identified three main reasons to explain the low fish Hg bioaccumulation: disconnection of the aquatic food web from the high MMHg zone, simple food web structures, and biodilution effect at the base of the food chain in this eutrophic reservoir.

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.

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.

Reservoir compensation releases : impact on the macroinvertebrate community of the Derwent River, Northumberland, UK-A longitudinal

River regulation for the purposes of public water supply causes the flow regime downstream of a dam to change. Traditionally, in the UK, such regulation was accompanied by requirements for reservoir releases to compensate downstream water users (e.g. industry) for the loss of natural flow (compensation flows). In this article, we compare a unique pre-impoundment macroinvertebrate data set for a regulated upland river with survey data post-impoundment. This allows a longitudinal assessment of the response of the system to regulation. The Derwent River, Northumberland, was impounded in 1966. Impacts on the hydrological regime were quantified by comparing long-term hydrographs, flow duration curves, flow ranges and flashiness indices for the pre-impoundment and post-impoundment periods. The comparison of changes in macroinvertebrate richness and diversity post-impoundment showed that the change in flow regime has had limited effect on the ecological community structure. The flow regime of the Derwent River has become less flashy with fewer extreme events, and the richness and the diversity of macroinvertebrates have, in some cases, increased and at worst have not deteriorated. We suggest that this reflects the strict compensation regime, which has guaranteed minimum flows at all times. Copyright (c) 2012 John Wiley & Sons, Ltd.

Long-term antiretroviral treatment initiated at primary HIV-1 infection affects the size, composition, and decay kinetics of the reservoir of HIV-1-infected CD4 T cells.

Initiation of antiretroviral therapy during the earliest stages of HIV-1 infection may limit the seeding of a long-lasting viral reservoir, but long-term effects of early antiretroviral treatment initiation remain unknown. Here, we analyzed immunological and virological characteristics of nine patients who started antiretroviral therapy at primary HIV-1 infection and remained on suppressive treatment for >10 years; patients with similar treatment duration but initiation of suppressive therapy during chronic HIV-1 infection served as controls. We observed that independently of the timing of treatment initiation, HIV-1 DNA in CD4 T cells decayed primarily during the initial 3 to 4 years of treatment. However, in patients who started antiretroviral therapy in early infection, this decay occurred faster and was more pronounced, leading to substantially lower levels of cell-associated HIV-1 DNA after long-term treatment. Despite this smaller size, the viral CD4 T cell reservoir in persons with early treatment initiation consisted more dominantly of the long-lasting central-memory and T memory stem cells. HIV-1-specific T cell responses remained continuously detectable during antiretroviral therapy, independently of the timing of treatment initiation. Together, these data suggest that early HIV-1 treatment initiation, even when continued for >10 years, is unlikely to lead to viral eradication, but the presence of low viral reservoirs and durable HIV-1 T cell responses may make such patients good candidates for future interventional studies aiming at HIV-1 eradication and cure. IMPORTANCE: Antiretroviral therapy can effectively suppress HIV-1 replication to undetectable levels; however, HIV-1 can persist despite treatment, and viral replication rapidly rebounds when treatment is discontinued. This is mainly due to the presence of latently infected CD4 T cells, which are not susceptible to antiretroviral drugs. Starting treatment in the earliest stages of HIV-1 infection can limit the number of these latently infected cells, raising the possibility that these viral reservoirs are naturally eliminated if suppressive antiretroviral treatment is continued for extremely long periods of time. Here, we analyzed nine patients who started on antiretroviral therapy within the earliest weeks of the disease and continued treatment for more than 10 years. Our data show that early treatment accelerated the decay of infected CD4 T cells and led to very low residual levels of detectable HIV-1 after long-term therapy, levels that were otherwise detectable in patients who are able to maintain a spontaneous, drug-free control of HIV-1 replication. Thus, long-term antiretroviral treatment started during early infection cannot eliminate HIV-1, but the reduced reservoirs of HIV-1 infected cells in such patients may increase their chances to respond to clinical interventions aiming at inducing a drug-free remission of HIV-1 infection.

Deposits of the sandy braided South Saskatchewan River: Implications for the use of modern analogs in reconstructing channel dimensions in reservoir characterization

Estimation of the dimensions of fluvial geobodies from core data is a notoriously difficult problem in reservoir modeling. To try and improve such estimates and, hence, reduce uncertainty in geomodels, data on dunes, unit bars, cross-bar channels, and compound bars and their associated deposits are presented herein from the sand-bed braided South Saskatchewan River, Canada. These data are used to test models that relate the scale of the formative bed forms to the dimensions of the preserved deposits and, therefore, provide an insight as to how such deposits may be preserved over geologic time. The preservation of bed-form geometry is quantified by comparing the Alluvial architecture above and below the maximum erosion depth of the modem channel deposits. This comparison shows that there is no significant difference in the mean set thickness of dune cross-strata above and below the basal erosion surface of the contemporary channel, thus suggesting that dimensional relationships between dune deposits and the formative bed-form dimensions are likely to be valid from both recent and older deposits. The data show that estimates of mean bankfull flow depth derived from dune, unit bar, and cross-bar channel deposits are all very similar. Thus, the use of all these metrics together can provide a useful check that all components and scales of the alluvial architecture have been identified correctly when building reservoir models. The data also highlight several practical issues with identifying and applying data relating to cross-strata. For example, the deposits of unit bars were found to be severely truncated in length and width, with only approximately 10% of the mean bar-form length remaining, and thus making identification in section difficult. For similar reasons, the deposits of compound bars were found to be especially difficult to recognize, and hence, estimates of channel depth based on this method may be problematic. Where only core data are available (i.e., no outcrop data exist), formative flow depths are suggested to be best reconstructed using cross-strata formed by dunes. However, theoretical relationships between the distribution of set thicknesses and formative dune height are found to result in slight overestimates of the latter and, hence, mean bankfull flow depths derived from these measurements. This article illustrates that the preservation of fluvial cross-strata and, thus, the paleohydraulic inferences that can be drawn from them, are a function of the ratio of the size and migration rate of bed forms and the time scale of aggradation and channel migration. These factors must thus be considered when deciding on appropriate length:thickness ratios for the purposes of object-based modeling in reservoir characterization.

Large slope deformations detection and monitoring along shores of the Potrerillos dam reservoir, Argentina, based on a small-baseline InSAR approach

The Argentina National Road 7 that crosses the Andes Cordillera within the Mendoza province to connect Santiago de Chile and Buenos Aires is particularly affected by natural hazards requiring risk management. Integrated in a research plan that intends to produce landslide susceptibility maps, we aimed in this study to detect large slope movements by applying a satellite radar interferometric analysis using Envisat data, acquired between 2005 and 2010. We were finally able to identify two large slope deformations in sandstone and clay deposits along gentle shores of the Potrerillos dam reservoir, with cumulated displacements higher than 25mm in 5years and towards the reservoir. There is also a body of evidences that these large slope deformations are actually influenced by the seasonal reservoir level variations. This study shows that very detailed information, such as surface displacements and above all water level variation, can be extracted from spaceborne remote sensing techniques; nevertheless, the limitations of InSAR for the present dataset are discussed here. Such analysis can then lead to further field investigations to understand more precisely the destabilising processes acting on these slope deformations.