Coupling multipoint flux mixed finite element methodswith continuous Galerkin methods for poroelasticity

We study the numerical approximation on irregular domains with general grids of the system of poroelasticity, which describes fluid flow in deformable porous media. The flow equation is discretized by a multipoint flux mixed finite element method and the displacements are approximated by a continuous Galerkin finite element method. First-order convergence in space and time is established in appropriate norms for the pressure, velocity, and displacement. Numerical results are presented that illustrate the behavior of the method.     

Geochemical features and origin of natural gas in heavy oil area of the Western Slope,Songliao Basin,China

The Western Slope of the Songliao Basin is rich in heavy oil resources (>70 × 10⁸ bbl), around which there are shallow gas reservoirs (∼1.0 × 10¹² m³). The gas is dominated by methane with a dryness over 0.99, and the non-hydrocarbon component being overwelmingly nitrogen. Carbon isotope composition of methane and its homologs is depleted in ¹³C, with δ¹³C₁ values being in the range of −55‰ to −75‰, δ¹³C₂ being in the range of −40‰ to −53‰ and δ¹³C₃ being in the range of −30‰ to −42‰, respectively. These values differ significantly from those solution gases source in the Daqing oilfield. This study concludes that heavy oils along the Western Slope were derived from mature source rocks in the Qijia-Gulong Depression, that were biodegraded. The low reservoir temperature (30–50 °C) and low salinity of formation water with neutral to alkaline pH (NaHCO₃) appeared ideal for microbial activity and thus biodegradation. Natural gas along the Western Slope appears mainly to have originated from biodegradation and the formation of heavy oil. This origin is suggested by the heavy δ¹³C of CO₂ (−18.78‰ to 0.95‰) which suggests that the methane was produced via fermentation as the terminal decomposition stage of the oil.     

Laser Ablation (U-Th-Sm)/He Dating of Zircons: Analytical Age Bias and its Consequence for Studying Detrital Zircon Populations

The in situ (U-Th-Sm)/He and U/Pb laser-ablation double-dating procedure is a valuable method that can provide a large dataset relatively efficiently in contrast with conventional bulk helium thermochronometry. In this study, we evaluate the potential age error associated with the double ablation procedure and report the in situ (U-Th-Sm)/He double-ablation dating of 249 zircons from the Fish Canyon Tuff locality. With LA-ICP-MS pseudo-depth profiling and 3D numerical modelling, we show that the concentric double-ablation procedure in minerals with U-Th-Sm zoning can generate a significant (U-Th-Sm)/He age error (positive or negative), resulting in over-scattering and/or an offset of the mean age. Pseudo-depth profiling is insufficient to predict the individual age error, partly because of the superimposed ablations. To evaluate the consequence of this inherent bias, we confront a synthetic age distribution to the error expected for U-Th-Sm zoned zircons analysed with double-ablation (U-Th-Sm)/He thermochronometry. As expected, a strong age bias causes the spreading of peak ages, downgrading the original signal. Yet, the throughput of the ablation-based method can allow intra- and inter-sample peak age identification and comparison, and the coupling of (U-Th-Sm)/He and U/Pb ages extends our ability to deconvolute a multimodal age spectrum.     

Performance studies of the fixed stress split algorithm for immiscible two-phase flow coupled with linear poromechanics

In this work, we measure the performance of the fixed stress split algorithm for the immiscible water-oil flow coupled with linear poromechanics. The two-phase flow equations are solved on general hexahedral elements using the multipoint flux mixed finite element method whereas the poromechanics equations are discretized using the conforming Galerkin method. We introduce a rigorous calculation of the update in poroelastic properties during the iterative solution of the coupled system equations. The effects of the coupling parameter on the performance of the fixed stress algorithm is demonstrated in two field studies: the Frio oil reservoir and the Cranfield injection site.

     

Probabilistic estimates of future changes in California temperature and precipitation using statistical and dynamical downscaling

Sixteen global general circulation models were used to develop probabilistic projections of temperature (T) and precipitation (P) changes over California by the 2060s. The global models were downscaled with two statistical techniques and three nested dynamical regional climate models, although not all global models were downscaled with all techniques. Both monthly and daily timescale changes in T and P are addressed, the latter being important for a range of applications in energy use, water management, and agriculture. The T changes tend to agree more across downscaling techniques than the P changes. Year-to-year natural internal climate variability is roughly of similar magnitude to the projected T changes. In the monthly average, July temperatures shift enough that that the hottest July found in any simulation over the historical period becomes a modestly cool July in the future period. Januarys as cold as any found in the historical period are still found in the 2060s, but the median and maximum monthly average temperatures increase notably. Annual and seasonal P changes are small compared to interannual or intermodel variability. However, the annual change is composed of seasonally varying changes that are themselves much larger, but tend to cancel in the annual mean. Winters show modestly wetter conditions in the North of the state, while spring and autumn show less precipitation. The dynamical downscaling techniques project increasing precipitation in the Southeastern part of the state, which is influenced by the North American monsoon, a feature that is not captured by the statistical downscaling.