Analysis of Lateral Fluid Gradients From DFA Measurements and Simulation of Reservoir Fluid Mixing Processes Over Geologic Time

Author(s):  
Qing Chen ◽  
◽  
Morten Kristensen ◽  
Yngve Bolstad Johansen ◽  
Vladislav Achourov ◽  
...  

Downhole fluid analysis (DFA) is one pillar of reservoir fluid geodynamics (RFG). DFA measurements provide both vertical and lateral fluid gradient data. These gradients, especially the asphaltene gradient derived from accurate optical density (OD) measurements, are critical in thermodynamic analysis to assess equilibration level and identify RFG processes. Recently, an RFG study was conducted using DFA and laboratory data from an oil field in the Norwegian North Sea. Fluid OD gradients show equilibrated asphaltenes in most of the reservoir, with a lateral variation of 20%. This indicates connectivity, which is confirmed by three years of production data. Two outliers are off the asphaltene equilibrium curve implying isolated sections, one each on the extreme east and west flank. Their asphaltene fraction varies by a factor of six. Such a difference reveals that different charge fluids entered the reservoir, and the equilibrated asphaltenes are the result of an after-charge mixing process. Meanwhile, different gas-oil contacts (GOCs) exist in the reservoir, indicating a lateral solution-gas gradient. Geochemistry analysis shows the same level of mild biodegradation in all the fluid samples, including those from two isolated sections. This means that biodegraded oil spills into the whole reservoir with little or no in-reservoir biodegradation. Furthermore, lateral asphaltene gradients at different times after charge have been preserved; it was a factor of six in asphaltenes content initially and is now 20% in the present day. This unique data set provides a valuable constraint to simulate reservoir fluid after-charge mixing processes to present day, aiming to investigate the factors impacting the evolution of lateral composition gradients in geologic time in a connected reservoir. Numerical simulations were performed over geologic time in reservoirs filled by oil with a lateral density gradient, which imitates the lateral compositional gradients in the gas-oil ratio (GOR) and asphaltenes measured in the above oil field. Simulations show that this lateral gradient creates lateral differential pressures and causes a countercurrent fluid flow forming a convection cell. In reservoirs with realistic vertical-to-horizontal aspect ratios, such fluid flows are not rapid, and lateral gradients can be partially retained in moderate geologic times. Additionally, diffusion was included in the simulation. The reservoir model was initialized with two GOCs producing subtle lateral GOR and density gradients. The simulated mixing process transports gas from higher GOR regions to lower GOR regions and reduces the GOC difference. However, the flux of solution gas transport is small. Consequently, we conclude that lateral GOR and asphaltene gradients can persist for moderate geologic time, which is consistent with observation from the field.

2020 ◽  
Author(s):  
Qing Chen ◽  
◽  
Morten Kristensen ◽  
Yngve Bolstad Johansen ◽  
Vladislav Achourov ◽  
...  

2015 ◽  
Vol 15 (1) ◽  
pp. 253-272 ◽  
Author(s):  
M. R. Canagaratna ◽  
J. L. Jimenez ◽  
J. H. Kroll ◽  
Q. Chen ◽  
S. H. Kessler ◽  
...  

Abstract. Elemental compositions of organic aerosol (OA) particles provide useful constraints on OA sources, chemical evolution, and effects. The Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) is widely used to measure OA elemental composition. This study evaluates AMS measurements of atomic oxygen-to-carbon (O : C), hydrogen-to-carbon (H : C), and organic mass-to-organic carbon (OM : OC) ratios, and of carbon oxidation state (OS C) for a vastly expanded laboratory data set of multifunctional oxidized OA standards. For the expanded standard data set, the method introduced by Aiken et al. (2008), which uses experimentally measured ion intensities at all ions to determine elemental ratios (referred to here as "Aiken-Explicit"), reproduces known O : C and H : C ratio values within 20% (average absolute value of relative errors) and 12%, respectively. The more commonly used method, which uses empirically estimated H2O+ and CO+ ion intensities to avoid gas phase air interferences at these ions (referred to here as "Aiken-Ambient"), reproduces O : C and H : C of multifunctional oxidized species within 28 and 14% of known values. The values from the latter method are systematically biased low, however, with larger biases observed for alcohols and simple diacids. A detailed examination of the H2O+, CO+, and CO2+ fragments in the high-resolution mass spectra of the standard compounds indicates that the Aiken-Ambient method underestimates the CO+ and especially H2O+ produced from many oxidized species. Combined AMS–vacuum ultraviolet (VUV) ionization measurements indicate that these ions are produced by dehydration and decarboxylation on the AMS vaporizer (usually operated at 600 °C). Thermal decomposition is observed to be efficient at vaporizer temperatures down to 200 °C. These results are used together to develop an "Improved-Ambient" elemental analysis method for AMS spectra measured in air. The Improved-Ambient method uses specific ion fragments as markers to correct for molecular functionality-dependent systematic biases and reproduces known O : C (H : C) ratios of individual oxidized standards within 28% (13%) of the known molecular values. The error in Improved-Ambient O : C (H : C) values is smaller for theoretical standard mixtures of the oxidized organic standards, which are more representative of the complex mix of species present in ambient OA. For ambient OA, the Improved-Ambient method produces O : C (H : C) values that are 27% (11%) larger than previously published Aiken-Ambient values; a corresponding increase of 9% is observed for OM : OC values. These results imply that ambient OA has a higher relative oxygen content than previously estimated. The OS C values calculated for ambient OA by the two methods agree well, however (average relative difference of 0.06 OS C units). This indicates that OS C is a more robust metric of oxidation than O : C, likely since OS C is not affected by hydration or dehydration, either in the atmosphere or during analysis.


Universe ◽  
2021 ◽  
Vol 7 (2) ◽  
pp. 25
Author(s):  
Umberto Battino ◽  
Claudia Lederer-Woods ◽  
Borbála Cseh ◽  
Pavel Denissenkov ◽  
Falk Herwig

The slow neutron-capture process (s-process) efficiency in low-mass AGB stars (1.5 < M/M⊙ < 3) critically depends on how mixing processes in stellar interiors are handled, which is still affected by considerable uncertainties. In this work, we compute the evolution and nucleosynthesis of low-mass AGB stars at low metallicities using the MESA stellar evolution code. The combined data set includes models with initial masses Mini/M⊙=2 and 3 for initial metallicities Z=0.001 and 0.002. The nucleosynthesis was calculated for all relevant isotopes by post-processing with the NuGrid mppnp code. Using these models, we show the impact of the uncertainties affecting the main mixing processes on heavy element nucleosynthesis, such as convection and mixing at convective boundaries. We finally compare our theoretical predictions with observed surface abundances on low-metallicity stars. We find that mixing at the interface between the He-intershell and the CO-core has a critical impact on the s-process at low metallicities, and its importance is comparable to convective boundary mixing processes under the convective envelope, which determine the formation and size of the 13C-pocket. Additionally, our results indicate that models with very low to no mixing below the He-intershell during thermal pulses, and with a 13C-pocket size of at least ∼3 × 10−4 M⊙, are strongly favored in reproducing observations. Online access to complete yield data tables is also provided.


2021 ◽  
Author(s):  
Elizabeth Ruiz ◽  
Brandon Thibodeaux ◽  
Christopher Dorion ◽  
Herman Mukisa ◽  
Majid Faskhoodi ◽  
...  

Abstract Optimized geomodeling and history matching of production data is presented by utilizing an integrated rock and fluid workflow. Facies identification is performed by use of image logs and other geological information. In addition, image logs are used to help define structural geodynamic processes that occurred in the reservoir. Methods of reservoir fluid geodynamics are used to assess the extent of fluid compositional equilibrium, especially the asphaltenes, and thereby the extent of connectivity in these facies. Geochemical determinations are shown to be consistent with measurements of compositional thermodynamic equilibrium. The ability to develop the geo-scenario of the reservoir, the coherent evolution of rock and contained fluids in the reservoir over geologic time, improves the robustness of the geomodel. In particular, the sequence of oil charge, compositional equilibrium, fault block throw, and primary biogenic gas charge are established in this middle Pliocene reservoir with implications for production, field extension,and local basin exploration. History matching of production data prove the accuracy of the geomodel; nevertheless, refinements to the geomodel and improved history matching were obtained by expanded deterministic property estimation from wireline log and other data. Theearly connection of fluid data, both thermodynamic and geochemical, with relevant facies andtheir properties determination enables a more facile method to incorporate this data into the geomodel. Logging data from future wells in the field can be imported into the geomodel allowingdeterministic optimization of this model long after production has commenced. While each reservoir is unique with its own idiosyncrasies, the workflow presented here is generally applicable to all reservoirs and always improves reservoir understanding.


Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede

Field ‘K’ is a mature field in the coastal swamp onshore Niger delta, which has been producing since 1960. As a huge producing field with some potential for further sustainable production, field monitoring is therefore important in the identification of areas of unproduced hydrocarbon. This can be achieved by comparing production data with the corresponding changes in acoustic impedance observed in the maps generated from base survey (initial 3D seismic) and monitor seismic survey (4D seismic) across the field. This will enable the 4D seismic data set to be used for mapping reservoir details such as advancing water front and un-swept zones. The availability of good quality onshore time-lapse seismic data for Field ‘K’ acquired in 1987 and 2002 provided the opportunity to evaluate the effect of changes in reservoir fluid saturations on time-lapse amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic impedance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir fluid saturations were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where reservoir porosity is 0.25%. In this paper, we demonstrated the use of repeat Seismic to delineate swept zones and areas hit with water override in a producing onshore reservoir.


2011 ◽  
Vol 134 (2) ◽  
Author(s):  
V. V. Dharaiya ◽  
S. G. Kandlikar

Study of fluid flow characteristics at microscale is gaining importance with shrinking device sizes. Better understanding of fluid flow and heat transfer in microchannels will have important implications in electronic chip cooling, heat exchangers, MEMS, and microfluidic devices. Due to short lengths employed in microchannels, entrance header effects can be significant and need to be investigated. In this work, three dimensional model of microchannels, with aspect ratios (α = a/b) ranging from 0.1 to 10, are numerically simulated using CFD software tool fluent. Heat transfer effects in the entrance region of microchannel are presented by plotting average Nusselt number as a function of nondimensional axial length x*. The numerical simulations with both circumferential and axial uniform heat flux (H2) boundary conditions are validated for existing data set for four wall heat flux case. Large numerical data sets are generated in this work for rectangular cross-sectional microchannels with heating on three walls, two opposing walls, one wall, and two adjacent walls under H2 boundary condition. This information can provide better understanding and insight into the transport processes in the microchannels. Although the results are seen as relevant in microscale applications, they are applicable to any sized channels. Based on the numerical results obtained for the whole range, generalized correlations for Nusselt numbers as a function of channel aspect ratio are presented for all the cases. The predicted correlations for Nusselt numbers can be very useful resource for the design and optimization of microchannel heat sinks and other microfluidic devices.


2021 ◽  
Author(s):  
Daiki Nomura ◽  
Alison Web ◽  
Yuhong Li ◽  
Manuel Dall’osto ◽  
Katrin Schmidt ◽  
...  

&lt;p&gt;We undertook a lead survey during the international drift campaign MOSAiC, Leg 5 (from 22 August to 17 September 2020) to understand the effects of lead width variation, re-freezing, and mixing events on lead water vertical structure. At the beginning of the survey period, the freshwater layer was occupied for the top 1 m depth and there were strong vertical gradients in temperature, salinity, and dissolved oxygen (DO) within 1 m depth: from 0.0&amp;#176;C to &amp;#8211;1.6&amp;#176;C for temperature, from 0.0 to 31.4 psu for salinity, and 10.5 to 13.5 mg L&lt;sup&gt;&amp;#8211;1&lt;/sup&gt; for DO. A strong DO minimum layer corresponded with a salinity of 25 psu, and usually occurred at the freshwater&amp;#8211;seawater interface at approx. 1 m depth, most likely as a result of an accumulation of organic matter and ongoing degradation/respiration processes at this interface. However, during the survey period, these strong gradients weakened and reduced the freshwater layer thickness (FLT). In the first half of the sampling period (until 4 September), FLT changed due to variations in lead width: as lead width increased, FLT decreased due to a stretching of the freshwater layer. In the second half of the sampling period, FLT was controlled by the surface ice formation (re-freezing) and mixing processes along the lower boundary of the freshwater layer. Surface ice formation removed freshwater and the formation of surface ice (about 0.2 m thick) explains 20% of the reduction of FLT. The remaining 80% of the reduction of FLT was due to the mixing process within the water column that was initiated by cooling and re-freezing. This mixing process diluted the salinity from 31.6 to 29.3 psu in the water below freshwater layer towards the end of the survey period. Our results indicate that lead water structure can change rapidly and dynamically and that this has significant effects on the biogeochemical exchange between lead systems and the atmosphere.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document