Correction to: Finite Element Solvers for Biot’s Poroelasticity Equations in Porous Media

Mathematical Geosciences - A correction to this paper has been published: https://doi.org/10.1007/s11004-021-09942-0     

Anisotropic Rock Poroelasticity Evolution in Ultra-low Permeability Sandstones under Pore Pressure, Confining Pressure, and Temperature: Experiments with Biot’s Coefficient

This study aimed to show anisotropic poroelasticity evolution in ultra-low permeability reservoirs under pore pressure, confining pressure, and temperature. Several groups of experiments examining Biot’s coefficient under different conditions were carried out. Results showed that Biot’s coefficient decreased with increased pore pressure, and the variation trend is linear, but the decreasing rate is variable between materials. Biot’s coefficient increased with increased confining pressure; the variation trend is linear, but the increasing rate varies by material as well. Generally, Biot’s coefficient remains stable with increased temperature. Lithology, clay mineral content, particle arrangement, and pore arrangement showed impacts on Biot’s coefficient. For strong hydrophilic clay minerals, expansion in water could result in a strong surface adsorption reaction, which could result in an increased fluid bulk modulus and higher Biot’s coefficient. For skeleton minerals with strong lipophilicity, such as quartz and feldspar, increased oil saturation will also result in an adsorption reaction, leading to increased fluid bulk modulus and a higher Biot’s coefficient. The study’s conclusions provide evidence of poroelasticity evolution of ultra-low permeability and help the enhancing oil recovery (EOR) process.     

Groundwater Behaviour in Layered Slope Rock Mass And Finite Element Simulation

The paper deals with the relatively weak permeable layer in a layered structure of a slope rock mass and itseffect on groundwater behaviour in the rock mass. When the whole groundwater table is declined in such typeof slope rock mass, the permeability of the relatively weak permeable layer in the middle will be several timeslower than that of the upper and lower permeable layers and consequently the water level in layers of differentpermeability will fall in different rates which will lead to a discontinuous phenomenon of the water body. Amethod of finite element simulation on a vertical two-dimensional profile is suggested for the study of the phe-nomenon and its application to practical situation is described. The calculation result shows that the existingweak permeable layer in a layered slope rock mass exerts unnegligible influence on the variation of water levelin the rock mass of a slope. The result also proves the applicability of the finite element simulation method.     

Analysis of the Wenchuan Ms8.0 Earthquake’s Co-seismic Stress and Displacement Change by Using the Finite Element Method

The variation of in situ stress before and after earthquakes is an issue studied by geologists. In this paper, on the basis of the fault slip dislocation model of Wenchuan Ms8.0 earthquake, the changes of co-seismic displacement and the distribution functions of stress tensor around the Longmen Shan fault zone are calculated. The results show that the co-seismic maximum surface displacement is 4.9 m in the horizontal direction and 6.5 m in the vertical direction, which is almost consistent with the on-site survey and GPS observations. The co-seismic maximum horizontal stress in the hanging wall and footwall decreased sharply as the distance from the Longmen Shan fault zone increased. However, the vertical stress and minimum horizontal stress increased in the footwall and in some areas of the hanging wall. The study of the co-seismic displacement and stress was mainly focused on the long and narrow region along the Longmen Shan fault zone, which coincides with the distribution of the earthquake aftershocks. Therefore, the co-seismic stress only affects the aftershocks, and does not affect distant faults and seismic activities. The results are almost consistent with in situ stress measurements at the two sites before and after Wenchuan Ms8.0 earthquake. Along the fault plane, the co-seismic shear stress in the dip direction is larger than that in the strike direction, which indicates that the faulting mechanism of the Longmen Shan fault zone is a dominant thrust with minor strike-slipping. The results can be used as a reference value for future studies of earthquake mechanisms.     

Multiphase Transport of Tritium in Unsaturated Porous Media—Bare and Vegetated Soils

Tritium is a short-lived radioactive isotope (T _1/2=12.33 yr) produced naturally in the atmosphere by cosmic radiation but also released into the atmosphere and hydrosphere by nuclear activities (nuclear power stations, radioactive waste disposal). Tritium of natural or anthropogenic origin may end up in soils through tritiated rain, and may eventually appear in groundwater. Tritium in groundwater can be re-emitted to the atmosphere through the vadose zone. The tritium concentration in soil varies sharply close to the ground surface and is very sensitive to many interrelated factors like rainfall amount, evapotranspiration rate, rooting depth and water table position, rendering the modeling a rather complex task. Among many existing codes, SOLVEG is a one-dimensional numerical model to simulate multiphase transport through the unsaturated zone. Processes include tritium diffusion in both, gas and liquid phase, advection and dispersion for tritium in liquid phase, radioactive decay and equilibrium partitioning between liquid and gas phase. For its application with bare or vegetated (perennial vegetation or crops) soil surfaces and shallow or deep groundwater levels (contaminated or non-contaminated aquifer) the model has been adapted in order to include ground cover, root growth and root water uptake. The current work describes the approach and results of the modeling of a tracer test with tritiated water (7.3×10⁸ Bq m⁻³) in a cultivated soil with an underlying 14 m deep unsaturated zone (non-contaminated). According to the simulation results, the soil’s natural attenuation process is governed by evapotranspiration and tritium re-emission. The latter process is due to a tritium concentration gradient between soil air and an atmospheric boundary layer at the soil surface. Re-emission generally occurs during night time, since at day time it is coupled with the evaporation process. Evapotranspiration and re-emission removed considerable quantities of tritium and limited penetration of surface-applied tritiated water in the vadose zone to no more than ∼1–2 m. After a period of 15 months tritium background concentration in soil was attained.     

Three dimensional Finite Element modeling of seismic soil–structure interaction in soft soil

Earthquakes in regions underlain by soft clay have amply demonstrated the detrimental effects of soil–structure interaction (SSI) in such settings. This paper describes a new three dimensional Finite Element model utilizing linear elastic single degree of freedom (SDOF) structure and a nonlinear elasto-plastic constitutive model for soil behavior in order to capture the nonlinear foundation–soil coupled response under seismic loadings. Results from an experimental SSI centrifuge test were used to verify the reliability of the numerical model followed by parametric studies to evaluate performance of linear elastic structures underlain by soft saturated clay. The results of parametric study demonstrate that rigid slender (tall) structures are highly susceptible to the SSI effects including alteration of natural frequency, foundation rocking and excessive base shear demand. Structure–foundation stiffness and aspect ratios were found to be crucial parameters controlling coupled foundation–structure performance in flexible-base structures. Furthermore, frequency content of input motion, site response and structure must be taken into account to avoid occurrence of resonance problem.     

Finite Element Analysis of Crossing Flow in an Intersected Fracture by Penalty Function Formulation

On the basis of Navier-Stokes equation, a digital analogue of crossing flow is carried out by means of finiteelement analysis performed through penalty function formulation. The results of the digital analogue are ingood agreement with that of the previous hydraulic physical simulation. They both demonstrate the effect of"deflection flow" occurring at the fracture intersection that the fracture-water crossing flow deflects △q intothe wide fracture.     

Dance and wellbeing in Vancouver’s ‘A Healthy City for All’

Through the lens of the dancing body, this paper examines practices of health and wellbeing produced in response to City of Vancouver urban governance policies. In particular, it calls attention to the legislative onslaught by city government in the years abutting the 2010 Winter Olympics to cultivate and manage healthy people, communities, and environments. In an effort to sell Vancouver’s ‘liveability’, I argue City of Vancouver endorsed a new legislative alliance that merged a conspicuously Anglo-American wellbeing lexicon, favouring individual responsibility and self-governance, with the performing arts industries. Drawing upon interviews and performance-based research, the paper illustrates how Karen Jamieson’s community dance project Connect, created for the In the Heart of the City festival, embodies Vancouver’s tri-level legislative ambitions to nurture A Healthy City For All. This materialised through the crafting of a dance-health body practice (healthy people), by choreographing a sense of belonging among ‘at risk’ communities (healthy communities), and in the uniting of the arts and health professions in the process of ‘cleaning up’ disenfranchised neighbourhoods (healthy environments). In bringing together scholarship on cultures of wellbeing and creative dance practice, the article contributes to understandings of how the health-seeking subject is embodied and performed. It also offers a productive critique of the exclusionary nature of urban health legislation, and of the contested role artists and arts festivals can play in nurturing urban wellbeing and normalising inequalities.     

Development of methods for constructing models for intra-year irregularity in the Earth’s rotation

The development of methods for the construction of stochastic, dynamical models for intra-year irregularity of the Earth’s rotation is considered. A correlational model based on harmonically additive and parametrically random, colored and broadband, gravitational-tidal perturbations from the Sun and Moon is developed. One-dimensional and multi-dimensional characteristic functions are found for the case of Gaussian and non-Gaussian colored and broadband fluctuations in the irregularity of the Earth’s rotation. Examples of computer modeling of the irregularity in the Earth’s rotation based on a priori and a posteriori IERS data are presented.     

China’s Progress in Prospecting for Concealed Deposits by Applying Geo-electrochemical Method

Please?refer?to?the?attachment(s)?for?more?details.     

Time-lag effects in Parker’s dynamo

The effects of a time lag in the magnetic quenching of the α effect is considered for an oscillating magnetic field in a Parker dynamo. The hypothesis of a parametric resonance in the system is justified, a modification of the solution is found, and the appearance of processes with periods much longer than the fundamental oscillation period is demonstrated.     

Breakthrough in the China’s Separation Technology for Extraction of Rare Earth Elements

Separation technology of rare earth elements （REEs）, as the critical step in the separation process, had long been fraught with technical difficulty. A research project conducted by Baotou Shibo Rare Earth Extraction ＆ Equipment Co. Ltd., Baotou REE Research Institute, and Baotou Steel ＆ REE Group Hi-Tech Co. Ltd （Inner Mongolia）, has successfully solved the problem using a centrifugal extractor and advanced techniques to achieve a key breakthrough.     

Development of tomographic technology for the Earth’s crust in the shelf regions

On the basis of experimental data obtained during a comprehensive experiment in Vityaz Bay of the Sea of Japan using onshore laser strainmeters and a low-frequency hydroacoustic emitter generating complex phase-manipulated signals with a central frequency of 33 Hz, we developed the basic principles of contactless tomography of the Earth’s crust in the shelf regions of various seas, including those covered by ice, making it possible to determine efficiently the structure and composition of the upper Earth’s crust under seas.

     

A survey-based analysis of the public’s willingness for disaster relief in China

Natural Hazards - Meteorological disasters frequently occur in China and around the world. These natural hazards can cause huge economic losses and threaten the personal safety of citizens. The...     

Accounting for the influence of the Earth’s sphericity in three-dimensional density modelling

A method for transformation of the three-dimensional regional “flat” density models of the Earth’s crust and upper mantle to the “spherical” models and vice versa is proposed. A computation algorithm and a method of meaningful comparison of the vertical component of the gravity field of both models are presented.     

What’s in a name? Epistemic perspectives and Payments for Ecosystem Services policies in Nicaragua

As Payments for Ecosystem Services (PES) continues to gain attention as a policy tool for securing efficient and effective environmental governance, a rising tide of criticism warns of the potentially detrimental social–ecological consequences of nature commodification and ‘green neoliberalism’. These concerns are also expressed at international policy fora, where the market rhetoric has met with political resistance from countries belonging to the ‘Bolivarian Alliance for the Peoples of Our America’ (ALBA). But despite this ideological opposition, some ALBA countries are increasingly integrating PES into their environmental policies. In this article we consider the reasons underlying this apparent contradiction and relate it to the notion of ‘epistemic circulation’. On the basis of a study on the evolution of PES-thinking in Nicaragua (an ALBA member) and a reassessment of the supposed ‘success’ of an influential pilot project, we shed light on the forces driving the adoption of particular PES modes and contextualise practical difficulties to endorsing more critical approaches to the tool. Instead of either ideologically rejecting PES as a neoliberal evil or embracing it uncritically as the new panacea, we argue that it is precisely through the socio-political processes surrounding environmental governance debates that the application of PES is shaped. In practice, it may either contribute to an imposed and dispossessing form of capitalism, or tend towards a more negotiated and socio-culturally embedded version of it. Only through its reconceptualisation based on political–cultural primacy rather than market-fetishism can PES achieve its true potential within a broader strategy towards improved environmental governance.     

Solution of Kepler’s equation for rectilinear motion

The properties of the solution of the kinematic equation (Kepler’s equation) presented as a series in powers of a function of time are considered for the case of unperturbed, rectilinear elliptical and rectilinear hyperbolic motion. Kepler’s equation for unperturbed, rectilinear elliptical motion has the form E − sinE = z ³/6, where E is the eccentric anomaly determining the position in the orbit, z ³/6 is the mean anomaly, which is proportional to the time measured from an encounter, and its solution can be represented as a series in powers of z. It is established that the coefficients of the series are positive. The asymptotic for the coefficients in the region of convergence of the series is found, which covers the entire orbit. The series continues to converge over the entire boundary of the circle of convergence. The kinematic equation for unperturbed, rectilinear hyperbolic motion has the form sinhH − H = ζ ³/6, where, as before, ζ ³ is proportional to time. The substitution E = iH, z = iζ reduces one equation to the other. The series for the solution in the hyperbolic case differs from the series for the elliptical solution only in its alternating-sign coefficients. However, the region of convergence covers only part of the orbit in the hyperbolic case.