Control Factors of Solutional Voids in Feldspars and Favorable Zone Forecast of Chang 2 Oil Reservoir Group in the Midwest Ordos Basin

The development of pores in a clastic reservoir is one of the most important research subjects in oil-gas exploration and development, whereas the many reasons for the formation of secondary porosity have increased the degree of difficulty in such research. Thus the research aims are to discover the controlling factors of solutional voids in feldspars and to predict favorable regions for these voids. Macroscopic and systematic researches into the relationship between the kaolinite content in the feldspar solutional void developed area of the Chang 2 reservoir group of the Triassic Yanchang Formation in the Midwest Ordos Basin and the solutional void in feldspar have been made, and from this it can be determined that the kaolinite content has an indicative function to the distribution of the solutional void in feldspar. Solutional void in feldspar is relatively well developed at the area where kaolinite content is high. Although the factors affecting kaolinite content are complicated, yet that of the research area is mainly affected by the impact of the leaching atmospheric water acting on the palaeogeomorphology. Three favorable zone belts for the development of solutional voids in feldspars are forecasted on the basis of restoration of palaeogeomorphology.

GLOBAL SOLUTIONS TO 1D COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM WITH DENSITY-DEPENDENT VISCOSITY AND FREE-BOUNDARY

This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

Smart prediction of liquefaction-induced lateral spreading

The prediction of liquefaction-induced lateral spreading/displacement(Dh)is a challenging task for civil/geotechnical engineers.In this study,a new approach is proposed to predict Dh using gene expression programming(GEP).Based on statistical reasoning,individual models were developed for two topographies:free-face and gently sloping ground.Along with a comparison with conventional approaches for predicting the Dh,four additional regression-based soft computing models,i.e.Gaussian process regression(GPR),relevance vector machine(RVM),sequential minimal optimization regression(SMOR),and M5-tree,were developed and compared with the GEP model.The results indicate that the GEP models predict Dh with less bias,as evidenced by the root mean square error(RMSE)and mean absolute error(MAE)for training(i.e.1.092 and 0.815;and 0.643 and 0.526)and for testing(i.e.0.89 and 0.705;and 0.773 and 0.573)in free-face and gently sloping ground topographies,respectively.The overall performance for the free-face topology was ranked as follows:GEP>RVM>M5-tree>GPR>SMOR,with a total score of 40,32,24,15,and 10,respectively.For the gently sloping condition,the performance was ranked as follows:GEP>RVM>GPR>M5-tree>SMOR with a total score of 40,32,21,19,and 8,respectively.Finally,the results of the sensitivity analysis showed that for both free-face and gently sloping ground,the liquefiable layer thickness(T_(15))was the major parameter with percentage deterioration(%D)value of 99.15 and 90.72,respectively.

Construction of Dynamic Alloy Interfaces for Uniform Li Deposition in Li-Metal Batteries

It is well accepted that a lithiophilic interface can effectively regulate Li deposition behaviors,but the influence of the lithiophilic interface is gradually diminished upon continuous Li deposition that completely isolates Li from the lithiophilic metals.Herein,we perform in-depth studies on the creation of dynamic alloy interfaces upon Li deposition,arising from the exceptionally high diffusion coefficient of Hg in the amalgam solid solution.As a comparison,other metals such as Au,Ag,and Zn have typical diffusion coefficients of 10-20 orders of magnitude lower than that of Hg in the similar solid solution phases.This difference induces compact Li deposition pattern with an amalgam substrate even with a high areal capacity of 55 mAh cm^(-2).This finding provides new insight into the rational design of Li anode substrate for the stable cycling of Li metal batteries.

Wave interaction for a generalized higher-dimensional Boussinesq equation describing the nonlinear Rossby waves

Based on an algebraically Rossby solitary waves evolution model,namely an extended(2+1)-dimensional Boussinesq equation,we firstly introduced a special transformation and utilized the Hirota method,which enable us to obtain multi-complexiton solutions and explore the interaction among the solutions.These wave functions are then employed to infer the influence of background flow on the propagation of Rossby waves,as well as the characteristics of propagation in multi-wave running processes.Additionally,we generated stereogram drawings and projection figures to visually represent these solutions.The dynamical behavior of these solutions is thoroughly examined through analytical and graphical analyses.Furthermore,we investigated the influence of the generalized beta effect and the Coriolis parameter on the evolution of Rossby waves.

Localized wave solutions and interactions of the (2+1)-dimensional Hirota-Satsuma-Ito equation

This paper studies the(2+1)-dimensional Hirota-Satsuma-Ito equation.Based on an associated Hirota bilinear form,lump-type solution,two types of interaction solutions,and breather wave solution of the(2+1)-dimensional Hirota-Satsuma-Ito equation are obtained,which are all related to the seed solution of the equation.It is interesting that the rogue wave is aroused by the interaction between one-lump soliton and a pair of resonance stripe solitons,and the fusion and fission phenomena are also found in the interaction between lump solitons and one-stripe soliton.Furthermore,the breather wave solution is also obtained by reducing the two-soliton solutions.The trajectory and period of the one-order breather wave are analyzed.The corresponding dynamical characteristics are demonstrated by the graphs.

PDE Standardization Analysis and Solution of TypicalMechanics Problems

A numerical approach is an effective means of solving boundary value problems(BVPs).This study focuses on physical problems with general partial differential equations(PDEs).It investigates the solution approach through the standard forms of the PDE module in COMSOL.Two typical mechanics problems are exemplified:The deflection of a thin plate,which can be addressed with the dedicated finite element module,and the stress of a pure bending beamthat cannot be tackled.The procedure for the two problems regarding the three standard forms required by the PDE module is detailed.The results were in good agreement with the literature,indicating that the PDE module provides a promising means to solve complex PDEs,especially for those a dedicated finite element module has yet to be developed.

Dynamics of Nonlinear Waves in(2+1)-Dimensional Extended Boiti-Leon-Manna-Pempinelli Equation

Based on the Hirota bilinear method,this study derived N-soliton solutions,breather solutions,lump solutions and interaction solutions for the(2+1)-dimensional extended Boiti-Leon-Manna-Pempinelli equation.The dynamical characteristics of these solutions were displayed through graphical,particularly revealing fusion and ssion phenomena in the interaction of lump and the one-stripe soliton.

Solution approximations for a mathematical model of relativistic electrons with beta derivative

The main aim of this paper is to obtain the exact and semi-analytical solutions of the nonlinear Klein-Fock-Gordon(KFG)equation which is a model of relativistic electrons arising in the laser thermonuclear fusion with beta derivative.For this purpose,both the modified extended tanh-function(mETF)method and the homotopy analysis method(HAM)are used.While applying the mETF the chain rule for beta derivative and complex wave transform are used for obtaining the exact solution.The advantage of this procedure is that discretization or normalization is not required.By applying the mETF,the exact solutions are obtained.Also,by applying the HAM semi-analytical results for the considered equation are acquired.In HAM?curve gives us a chance to find the suitable value of the for the convergence of the solution series.Also,comparative graphical representations are given to show the effectiveness,reliability of the methods.The results show that the m ETF and HAM are reliable and applicable tools for obtaining the solutions of non-linear fractional partial differential equations that involve beta derivative.This study can bring a new perspective for studies on fractional differential equations.On the other hand,it can be said that scientists can apply the considered methods for different mathematical models arising in physics,chemistry,engineering,social sciences and etc.which involves fractional differentiation.Briefly the results may cause a new insight who studies on relativistic electron modelling.

Malargüe Site Remediation:A Successful Solution for Uranium Mill Tailings

Comisión Nacional de Energía Atómica (CNEA) has the responsibility for restoring uranium mining facilities once the operations have finished.CNEA,within its Environmental Program and in compliance with its legal responsibilities,decides to implement a restoration project for all sites related to the mining and processing of uranium ores.The Malargüe Site is located within the Province of Mendoza in the city of Malargüe.It is the first site to successfully complete its remediation.The activities consist of relocation of tailings to an engineering repository.The tailings management(encapsulation) and rehabilitation of the area was finished in June 2017.The remediation alternative for the ore tailings was selected after conducting comparative studies and submitted the project to the society for consideration.The objective of the encapsulation of the mineral tails is to isolate them from the environment,also proceeding with the decontamination and rehabilitation of the area (landscaping,post-closure monitoring and 20 years monitoring period).Encapsulation consisted of the construction of a containment cell for the mine tailings,to isolate them and prevent pollutants from entering the environment through the transfer routes.To clean the impacted areas,the soil was removed,it was incorporated into the encapsulation,and the filling was carried out with natural soils from the area.Remediation prevents radon transfer to the environment,as ^(222)Ra is an alpha emitter with a half-life of four days,which produces its own radioactive progeny.Radon progeny are solids,and when a ^(222)Ra nucleus emits an alpha particle into the air,the resulting ^(218)Po nucleus,momentarily electrically charged,adheres to any dust particle.Remediation prevents the discharge into the air containing radon and also containing dust particles charged with intensely radioactive radon progeny.The tasks mentioned make it possible to decrease radon emanation,reduce radiological risks to the public and prevent the entry of rainwater into the system.In addition,the containment system prevents the discharge of contaminated liquids into the environment,avoiding contamination of the groundwater.All these activities are according to the concepts of sustainability.

Hamiltonian system for the inhomogeneous plane elasticity of dodecagonal quasicrystal plates and its analytical solutions

A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions of the phonon and phason displacements are obtained further for the quasicrystal plates. In addition, the effectiveness of the approach is verified by comparison with the data of the finite integral transformation method.

THE LIMITING PROFILE OF SOLUTIONS FOR SEMILINEAR ELLIPTIC SYSTEMS WITH A SHRINKING SELF-FOCUSING CORE

In this paper,we consider the semilinear elliptic equation systems{△u+u=αQ_(n)(x)|u|^(α-2)|v|^(β)u in R^(N),-△v+v=βQ(x)|u|^(α)|v|^(β-2)v in R^(N),where N≥3,α,β>1,α+β<2^(*),2^(*)=2N/N-2 and Q_(n) are bounded given functions whose self-focusing cores{x∈R^(N)|Q_(n)(x)>0} shrink to a set with finitely many points as n→∞.Motivated by the work of Fang and Wang[13],we use variational methods to study the limiting profile of ground state solutions which are concentrated at one point of the set with finitely many points,and we build the localized concentrated bound state solutions for the above equation systems.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.

Effect of inorganic salt impurities on seeded precipitation of silica hydrate from sodium silicate solution

To clarify the precipitation of silica hydrate from the real desilication solutions of aluminosilicate solid wastes by adding seeds and improve integrated waste utilization,the seeded precipitation was studied using synthesized sodium silicate solution containing different inorganic salt impurities.The results show that sodium chloride,sodium sulfate,sodium carbonate,or calcium chloride can change the siloxy group structure.The number of high-polymeric siloxy groups decreases with increasing sodium chloride or sodium sulfate concentration,which is detrimental to seeded precipitation.Calcium chloride favors the polymerization of silicate ions,and even the chain groups precipitate with the precipitation of high-polymeric sheet and cage-like siloxy groups.The introduced sodium cations in sodium carbonate render a more open network structure of high-polymeric siloxy groups,although the carbonate ions favor the polymerization of siloxy groups.No matter how the four impurities affect the siloxy group structure,the precipitates are always amorphous opal-A silica hydrate.

GLOBAL WEAK SOLUTIONS FOR AN ATTRACTION-REPULSION CHEMOTAXIS SYSTEM WITH p-LAPLACIAN DIFFUSION AND LOGISTIC SOURCE

This paper is concerned with the following attraction-repulsion chemotaxis system with p-Laplacian diffusion and logistic source:■The system here is under a homogenous Neumann boundary condition in a bounded domainΩ ■ R^(n)(n≥2),with χ,ξ,α,β,γ,δ,k_(1),k_(2)> 0,p> 2.In addition,the function f is smooth and satisfies that f(s)≤κ-μs~l for all s≥0,with κ ∈ R,μ> 0,l> 1.It is shown that(ⅰ)if l> max{2k_(1),(2k_(1)n)/(2+n)+1/(p-1)},then system possesses a global bounded weak solution and(ⅱ)if k_(2)> max{2k_(1)-1,(2k_(1)n)/(2+n)+(2-p)/(p-1)} with l> 2,then system possesses a global bounded weak solution.

A FLEXIBLE OBJECTIVE-CONSTRAINT APPROACH AND A NEW ALGORITHM FOR CONSTRUCTING THE PARETO FRONT OF MULTIOBJECTIVE OPTIMIZATION PROBLEMS

In this article, a novel scalarization technique, called the improved objective-constraint approach, is introduced to find efficient solutions of a given multiobjective programming problem. The presented scalarized problem extends the objective-constraint problem. It is demonstrated that how adding variables to the scalarized problem, can lead to find conditions for (weakly, properly) Pareto optimal solutions. Applying the obtained necessary and sufficient conditions, two algorithms for generating the Pareto front approximation of bi-objective and three-objective programming problems are designed. These algorithms are easy to implement and can achieve an even approximation of (weakly, properly) Pareto optimal solutions. These algorithms can be generalized for optimization problems with more than three criterion functions, too. The effectiveness and capability of the algorithms are demonstrated in test problems.

Solute transport and geochemical modeling of the coastal quaternary aquifer, Delta Dahab Basin, South Sinai, Egypt

The wadi dahab delta is in a dry, arid coastal zone within Egypt’s south Sinai Peninsula’s eastern portion. The primary water source is the Quaternary coastal alluvial aquifer. The groundwater salinity varies from 890to 8213 mg/L, with a mean value of 3417 mg/L. The dissolved major ions have been used to calculate the seawater mixing index(SWMI) using a linear equation that discriminates the groundwater mostly affected by water–rock interaction(SWMI 1>) and other samples mixed with Seawater(SWMI < 1). The isotopic composition of groundwater for specifically chosen groundwater samples ranges from-0.645‰ to +5.212‰ for δ^(18)O and from-9.582‰ to + 22.778‰ for δ^(2)H, where the seawater represented by a Red Sea water sample(δ^(18)O + 1.64‰-δ^(2)H + 9.80‰) and reject brine water are considerably enriched the isotopic groundwater values. The geochemical NETPATH model constrained by the dissolved significant ions, isotopes, and the rock aquifer forming minerals as phases indicate the mixing percent with the seawater ranges from 9% to 97% of seawater from 91% to 3% of original recharge water. According to the SEAWAT 3-D flow models, seawater has penetrated the Northeastern Dahab delta aquifer, with the intrusion zone extending1500 m inland. The salt dissolution, upwelling of saline water, recharge from the upstream mountain block, and seawater encroachment are the primary aspects contributing to the deterioration of groundwater quality. These findings may have significance for effective groundwater withdrawal management in arid locations worldwide with similar hydrogeological systems.

Porous-DeepONet:Learning the Solution Operators of Parametric Reactive Transport Equations in Porous Media

Reactive transport equations in porous media are critical in various scientific and engineering disciplines,but solving these equations can be computationally expensive when exploring different scenarios,such as varying porous structures and initial or boundary conditions.The deep operator network(DeepONet)has emerged as a popular deep learning framework for solving parametric partial differential equations.However,applying the DeepONet to porous media presents significant challenges due to its limited capability to extract representative features from intricate structures.To address this issue,we propose the Porous-DeepONet,a simple yet highly effective extension of the DeepONet framework that leverages convolutional neural networks(CNNs)to learn the solution operators of parametric reactive transport equations in porous media.By incorporating CNNs,we can effectively capture the intricate features of porous media,enabling accurate and efficient learning of the solution operators.We demonstrate the effectiveness of the Porous-DeepONet in accurately and rapidly learning the solution operators of parametric reactive transport equations with various boundary conditions,multiple phases,and multiphysical fields through five examples.This approach offers significant computational savings,potentially reducing the computation time by 50–1000 times compared with the finite-element method.Our work may provide a robust alternative for solving parametric reactive transport equations in porous media,paving the way for exploring complex phenomena in porous media.

An approximate analytical model for unconventional reservoir considering variable matrix blocks and simultaneous matrix depletion

In regard to unconventional oil reservoirs,the transient dual-porosity and triple-porosity models have been adopted to describe the fluid flow in the complex fracture network.It has been proven to cause inaccurate production evaluations because of the absence of matrix-macrofracture communication.In addition,most of the existing models are solved analytically based on Laplace transform and numerical inversion.Hence,an approximate analytical solution is derived directly in real-time space considering variable matrix blocks and simultaneous matrix depletion.To simplify the derivation,the simultaneous matrix depletion is divided into two parts:one part feeding the macrofractures and the other part feeding the microfractures.Then,a series of partial differential equations(PDEs)describing the transient flow and boundary conditions are constructed and solved analytically by integration.Finally,a relationship between oil rate and production time in real-time space is obtained.The new model is verified against classical analytical models.When the microfracture system and matrix-macrofracture communication is neglected,the result of the new model agrees with those obtained with the dual-porosity and triple-porosity model,respectively.Certainly,the new model also has an excellent agreement with the numerical model.The model is then applied to two actual tight oil wells completed in western Canada sedimentary basin.After identifying the flow regime,the solution suitably matches the field production data,and the model parameters are determined.Through these output parameters,we can accurately forecast the production and even estimate the petrophysical properties.