Formation of Water-Bearing Defects in Olivine in the Presence of Water–Hydrocarbon Fluid at 6.3 GPa and 1200°C

The main trends of water dissolution in Fe-bearing olivine have been investigated in the olivine–H₂O–hydrocarbon fluid system in experiments at a pressure of 6.3 GPa, a temperature of 1200°C, and hydrogen fugacity ( fH₂) buffered by the Mo–MoO₂ equilibrium. The content and contribution of ОH defects of different types in Fe-bearing olivines depend on the composition of reduced fluids in the system. As the fraction of hydrocarbons in the fluid increases, the H₂O content in olivine crystals decreases from 900 to 160–180 ppm, while the ОН absorption peaks become lower at high frequencies and occupy a larger part of the infrared spectrum in the low-frequency region. According to the experimental results, even the deepest seated mantle olivines with OH defects were not equilibrated with a fluid rich in light alkanes or oxygenated hydrocarbons.     

Foidite and Meimechite Lavas of Polar Siberia (Some Questions of Petrogenesis)

Doklady Earth Sciences - For the Permian–Triassic foidite and meimechite lavas of Polar Siberia, both the whole-rock petrochemistry and geochemistry and that of melt inclusions in olivine...     

The Stress–Strain State of Recent Structures in the Northeastern Sector of the Russian Arctic Region

Doklady Earth Sciences - Complex research to determine the stress–strain state of the Earth’s crust and the types of seismotectonic destruction for the northeastern sector of the...     

The Fragmentation in a Protoplanetary Disk and the Properties of the Resultant Self-Gravitating Clumps

We study the fragmentation properties in the protoplanetary disk and properties of the resultant self-gravitating clumps using our newly constructed disk model. Our disk model includes the mass inflall term from a molecular cloud core and the photoevaporation winds effect. We adopt the conventional fragmentation criterion to judge whether a protoplanetary disk can fragment. In this work, we follow our previous work to investigate the properties of the resultant self-gravitating clumps. In our calculation, the initial masses of the resultant self-gravitating clumps lie in the range of tens of M_J to more than one hundred of M_J, where M_J is the Jupiter mass. These initial masses can seemingly account for the masses of extrasolar planets in magnitude. We also calculate the subsequent gas accretion of clumps in 1.27 × 10⁴ yr after the formation of self-gravitating clumps. We find that the subsequent gas accretion of self-gravitating clumps is very efficient, and the clump masses grow to hundreds of M_J and the physical radii R_c of clumps increase to about 10 AU. Additionally, we also calculate the orbital migration of clumps. We find that most clumps have short migration timescale to be accreted onto the protostar, and only a small fraction of clumps have long migration timescale (>10⁶ yr) to successfully become gas giant planets. These results are consistent with previous studies.     

Spectropolarimetric Observations of Magnetic White Dwarfs with the SAO 6-m Telescope

The results of spectropolarimetric observations of a number of magnetic white dwarfs obtained on the 6-m optical telescope of the Special Astrophysical Observatory are presented. The observations were carried out using the SCORPIO focal aperture-ratio reducer in a spectropolarimetric regime. Two characteristic dependences of the degree of polarization on the wavelength are observed. For one group of objects, the degree of linear polarization grows with wavelength, suggesting that the alignment of atoms and molecules in Rydberg states in the atmosphere of the white dwarf due to the action of its magnetic field influences scattering processes. The second group of objects displays an increase in the degrees of both linear and circular polarization with wavelength, providing evidence for the presence of protoplanetary disks around these magnetic white dwarfs, in which the alignment of circumstellar grains leads to the observed behavior.     

Paleoearthquakes in the Uimon basin (Gorny Altai)

Paleoseismological studies confirm that the Uimon basin is thrust by its northern mountain border along the active South Terekta fault. The latest motion along the fault in the 7-8th centuries AD induced an earthquake with a magnitude of M_w= 7.4-7.7 and a shaking intensity of I = 9-11 on the MSK-64 scale. The same fault generated another event (M > 7, I = 9-10), possibly, about 16 kyr ago, which triggered gravity sliding. The rockslide dammed the Uimon valley and produced a lake, where lacustrine deposition began about 14 ± 1 kyr ago, and a later M > 7 (I = 9-10) earthquake at ~ 6 ka caused the dam collapse and the lake drainage. Traces of much older earthquakes that occurred within the Uimon basin are detectable from secondary deformation structures (seismites) in soft sediments deposited during the drainage of a Late Pleistocene ice-dammed lake between 100 and 90 ka and in ~ 77 ka alluvium. The magnitude and intensity of these paleoearthquakes were at least M > 5.0-5.5 and I > 6-7.     

Adjoint analysis of Buckley-Leverett and two-phase flow equations

This paper analyzes the adjoint equations and boundary conditions for porous media flow models, specifically the Buckley-Leverett equation, and the compressible two-phase flow equations in mass conservation form. An adjoint analysis of a general scalar hyperbolic conservation law whose primal solutions include a shock jump is initially presented, and the results are later specialized to the Buckley-Leverett equation. The non-convexity of the Buckley-Leverett flux function results in adjoint characteristics that are parallel to the shock front upstream of the shock and emerge from the shock front downstream of the shock. Thus, in contrast to the behavior of Burgers’ equation where the adjoint is continuous at a shock, the Buckley-Leverett adjoint, in general, contains a discontinuous jump across the shock. Discrete adjoint solutions from space-time discontinuous Galerkin finite element approximations of the Buckley-Leverett equation are shown to be consistent with the derived closed-form analytical solutions. Furthermore, a general result relating the adjoint equations for different (though equivalent) primal equations is used to relate the two-phase flow adjoints to the Buckley-Leverett adjoint. Adjoint solutions from space-time discontinuous Galerkin finite element approximations of the two-phase flow equations are observed to obey this relationship.     

High-dimensional geostatistical history matching

Hydrocarbon reservoir modelling and characterisation is a challenging subject within the oil and gas industry due to the lack of well data and the natural heterogeneities of the Earth’s subsurface. Integrating historical production data into the geo-modelling workflow, commonly designated by history matching, allows better reservoir characterisation and the possibility of predicting the reservoir behaviour. We present herein a geostatistical-based multi-objective history matching methodology. It starts with the generation of an initial ensemble of the subsurface petrophysical property of interest through stochastic sequential simulation. Each model is then ranked according the match between its dynamic response, after fluid flow simulation, and the observed available historical production data. This enables building regionalised Pareto fronts and the definition of a large ensemble of optimal subsurface Earth models that fit all the observed production data without compromising the exploration of the uncertainty space. The proposed geostatistical multi-objective history matching technique is successfully implemented in a benchmark synthetic reservoir dataset, the PUNQ-S3, where 12 objectives are targeted.     

Constraint energy minimizing generalized multiscale finite element method in the mixed formulation

This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast.     

Fahlore and Sphalerite from the Darasun Gold Deposit in the Eastern Transbaikal Region,Russia: I. Mineral Assemblages and Intergrowths,Chemical Composition,and Its Evolution

The mineral assemblages, mode of occurrence, and chemical compositions of coexisting fahlore and sphalerite from the Darasun gold deposit have been described. Three generations of fahlore and three generations of sphalerite have been recognized. The FeS content in sphalerite coexisting with fahlore ranges from 0.8 to 9.4 mol %. The complete solid solution series Fe-tetrahedrite–Zn-tetrahedrite–Fe-tennantite–Zn-tennantite reflected in Sb/(Sb + As) and Fe/(Fe + Zn) ratios ranging from 0 to 0.97 and from 0.07 to 1.00, respectively, with a predominant negative relationship between these ratios has been identified for the first time at the deposit. Stepped, oscillatory, and combined stepped-oscillatory growth zonings within fahlore grains and heterogeneous aggregates of fahlore have been found. Fahlore is enriched in As with respect to Sb, and Zn-tetrahedrite is followed by Fe- and Zn-tennantite from early to late generation; Zn-tetrahedrite is followed by Fe-tennantite in zoned grains and overgrown rims; sphalerite crystallized at decreased temperature and sulfur fugacity. The evolution of the chemical composition of fahlores was caused by the evolving temperature, fluid salinity, and conditions of metal migration.