Fitness of Four-Parameter Beta Distribution Function for Forecasting Gold Reserve and Its Production Lifespan in Ghana

Ghana, renowned for its abundant gold reserves, plays a significant role in the global mining industry. Effective management and accurate forecasting of these reserves are vital for sustainable resource utilization and economic planning. Forecasting gold reserves and estimating their production lifespan are complex tasks that require robust statistical models capable of capturing the underlying dynamics of gold deposit accumulation and extraction. To this end, the four-parameter Beta distribution function emerges as a promising candidate due to its flexibility and ability to handle non-negative data. This research aims to investigate the fitness and applicability of the four-parameter Beta distribution function for forecasting Ghana’s gold reserves and estimating the production lifespan of this precious resource. The empirical paper relied mainly on quarterly secondary datasets on gold reserve between the years 2009 and 2022 secured from the Minerals Commission of Ghana, Accra. Several known statistical distributions including Beta, Weibull, Normal, Logistic and Gamma were explored with Maximum Likelihood Estimation (MLE) and evaluated using model selection criteria as AIC and BIC. Goodness of Fits were evaluated using Kolmogorov-Smirnov Test (K-S), Cramer-Von Mises Statistic and Anderson-Darling Statistic. Based on the analysis conducted, the four-parameter Beta distribution provided the best fit for gold reserve in Ghana. At a 99.9% confidence level and considering the current annual average gold production estimate of 3,700,031.248 to 4,302,647.888 ounces, the projected lifespan of gold production in Ghana extends to the year 1,953,765. This astounding estimate suggests that the country’s gold reserves are expected to sustain production for an extended period, providing a critical resource for economic development and supporting the mining industry well into the distant future.

Fitness of Four-Parameter Beta Distribution Function for Forecasting Gold Reserve and Its Production Lifespan in Ghana

Ghana, renowned for its abundant gold reserves, plays a significant role in the global mining industry. Effective management and accurate forecasting of these reserves are vital for sustainable resource utilization and economic planning. Forecasting gold reserves and estimating their production lifespan are complex tasks that require robust statistical models capable of capturing the underlying dynamics of gold deposit accumulation and extraction. To this end, the four-parameter Beta distribution function emerges as a promising candidate due to its flexibility and ability to handle non-negative data. This research aims to investigate the fitness and applicability of the four-parameter Beta distribution function for forecasting Ghana’s gold reserves and estimating the production lifespan of this precious resource. The empirical paper relied mainly on quarterly secondary datasets on gold reserve between the years 2009 and 2022 secured from the Minerals Commission of Ghana, Accra. Several known statistical distributions including Beta, Weibull, Normal, Logistic and Gamma were explored with Maximum Likelihood Estimation (MLE) and evaluated using model selection criteria as AIC and BIC. Goodness of Fits were evaluated using Kolmogorov-Smirnov Test (K-S), Cramer-Von Mises Statistic and Anderson-Darling Statistic. Based on the analysis conducted, the four-parameter Beta distribution provided the best fit for gold reserve in Ghana. At a 99.9% confidence level and considering the current annual average gold production estimate of 3,700,031.248 to 4,302,647.888 ounces, the projected lifespan of gold production in Ghana extends to the year 1,953,765. This astounding estimate suggests that the country’s gold reserves are expected to sustain production for an extended period, providing a critical resource for economic development and supporting the mining industry well into the distant future.

Influence of High-Density Bedding Plane Characteristics on Hydraulic Fracture Propagation in Shale Oil Reservoir

The existence of high-density bedding planes is a typical characteristic of shale oil reservoirs.Understanding the behavior of hydraulic fracturing in high-density laminated rocks is significant for promoting shale oil production.In this study,a hydraulic fracturing model considering tensile failure and frictional slip of the bedding planes is established within the framework of the unified pipe-interface element method(UP-IEM).The model developed for simulating the interaction between the hydraulic fracture and the bedding plane is validated by comparison with experimental results.The hydraulic fracturing patterns in sealed and unsealed bedding planes are compared.Additionally,the effects of differential stress,bedding plane permeability,spacing,and the friction coefficient of the bedding plane are investigated.The results showed that a single main fracture crossing the bedding planes is more likely to form in sealed bedding planes under high differential stress.The decrease in bedding plane permeability and the increase in the friction coefficient also promote the fracture propagating perpendicular to the bedding planes.Shale with high-density bedding planes has a poorer fracturing effect than that with low-density bedding planes,as the hydraulic fracture is prone to initiate and propagate along the bedding planes.Moreover,higher injection pressure is needed to maintain fracture propagation along the bedding.An increase in bedding density will lead to a smaller fracturing area.Fracturing fluid seepage into the bedding planes slows shale fracturing.It is recommended that increasing the injection flow rate,selecting alternative fracturing fluids,and employing multi-well/multi-cluster fracturing may be efficient methods to improve energy production in shale oil reservoirs.

Subplanes of PG( 2, q 3 ) and the Ruled Varieties V 2 5 of PG( 6,q )

In this note we study subplanes of order q of the projective plane Π=PG( 2, q 3 ) and the ruled varieties V 2 5 of Σ=PG( 6,q ) using the spatial representation of Π in Σ, by fixing a hyperplane Σ ′ with a regular spread of planes. First are shown some configurations of the affine q-subplanes. Then to prove that a variety V 2 5 of Σ represents a non-affine subplane of order q of Π, after having shown basic incidence properties of it, such a variety V 2 5 is constructed by choosing appropriately the two directrix curves in two complementary subspaces of Σ. The result can be translated into further incidence properties of the affine points of V 2 5 . Then a maximal bundle of varieties V 2 5 having in common one directrix cubic curve is constructed.

THE NONLINEAR STABILITY OF PLANE PARALLEL SHEAR FLOWS WITH RESPECT TO TILTED PERTURBATIONS

The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.

Nonlinear constitutive models of rock structural plane and their applications

Structural planes play an important role in controlling the stability of rock engineering,and the influence of structural planes should be considered in the design and construction process of rock engineering.In this paper,mechanical properties,constitutive theory,and numerical application of structural plane are studied by a combination method of laboratory tests,theoretical derivation,and program development.The test results reveal the change laws of various mechanical parameters under different roughness and normal stress.At the pre-peak stage,a non-stationary model of shear stiffness is established,and threedimensional empirical prediction models for initial shear stiffness and residual stage roughness are proposed.The nonlinear constitutive models are established based on elasto-plastic mechanics,and the algorithms of the models are developed based on the return mapping algorithm.According to a large number of statistical analysis results,empirical prediction models are proposed for model parameters expressed by structural plane characteristic parameters.Finally,the discrete element method(DEM)is chosen to embed the constitutive models for practical application.The running programs of the constitutive models have been compiled into the discrete element model library.The comparison results between the proposed model and the Mohr-Coulomb slip model show that the proposed model can better describe nonlinear changes at different stages,and the predicted shear strength,peak strain and shear stiffness are closer to the test results.The research results of the paper are conducive to the accurate evaluation of structural plane in rock engineering.

A novel method for simulating nuclear explosion with chemical explosion to form an approximate plane wave: Field test and numerical simulation

A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole explosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the underground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simulation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.

Ultrasound-guided serratus anterior plane block enhances postoperative analgesia and recovery in thoracoscopic surgery

BACKGROUND The serratus anterior muscle,located in the lateral aspect of the thorax,plays a crucial role in shoulder movement and stability.Thoracoscopic surgery,while minimally invasive,often results in significant postoperative pain,complicating patient recovery and potentially extending hospital stays.Traditional anesthesia methods may not adequately address this pain,leading to increased complications such as agitation due to inadequate pain management.AIM To evaluate the application value of ultrasound-guided serratus anterior plane block(SAPB)in patients undergoing thoracoscopic surgery,focusing on its effects on postoperative analgesia and rehabilitation.METHODS Eighty patients undergoing thoracoscopic surgery between August 2021 and December 2022 were randomly divided into two groups:An observation group receiving ultrasound-guided SAPB and a control group receiving standard care without SAPB.Both groups underwent general anesthesia and were monitored for blood pressure,heart rate(HR),oxygen saturation,and pulse.The primary outcomes measured included mean arterial pressure(MAP),HR,postoperative visual analogue scale(VAS)scores for pain,supplemental analgesic use,and incidence of agitation.RESULTS The observation group showed significantly lower cortisol and glucose concentrations at various time points post-operation compared to the control group,indicating reduced stress responses.Moreover,MAP and HR levels were lower in the observation group during and after surgery.VAS scores were significantly lower in the observation group at 1 h,4 h,6 h,and 12 h post-surgery,and the rates of analgesic supplementation and agitation were significantly reduced compared to the control group.CONCLUSION Ultrasound-guided SAPB significantly improves postoperative analgesia and reduces agitation in patients undergoing thoracoscopic surgery.This technique stabilizes perioperative vital signs,decreases the need for supplemental analgesics,and minimizes postoperative pain and stress responses,underscoring its high application value in enhancing patient recovery and rehabilitation post-thoracoscopy.

An 8-Node Plane Hybrid Element for StructuralMechanics Problems Based on the Hellinger-Reissner Variational Principle

The finite element method (FEM) plays a valuable role in computer modeling and is beneficial to the mechanicaldesign of various structural parts. However, the elements produced by conventional FEM are easily inaccurate andunstable when applied. Therefore, developing new elements within the framework of the generalized variationalprinciple is of great significance. In this paper, an 8-node plane hybrid finite element with 15 parameters (PHQ8-15β) is developed for structural mechanics problems based on the Hellinger-Reissner variational principle.According to the design principle of Pian, 15 unknown parameters are adopted in the selection of stress modes toavoid the zero energy modes.Meanwhile, the stress functions within each element satisfy both the equilibrium andthe compatibility relations of plane stress problems. Subsequently, numerical examples are presented to illustrate theeffectiveness and robustness of the proposed finite element. Numerical results show that various common lockingbehaviors of plane elements can be overcome. The PH-Q8-15β element has excellent performance in all benchmarkproblems, especially for structures with varying cross sections. Furthermore, in bending problems, the reasonablemesh shape of the new element for curved edge structures is analyzed in detail, which can be a useful means toimprove numerical accuracy.

Review on dynamic analysis of road pavements under moving vehicles and plane strain conditions

This paper reviews works on the dynamic analysis of flexible and rigid pavements under moving vehicles on the basis of continuum-based plane strain models and linear theories.The purpose of this review is to provide in-formation about the existing works on the subject,critically discuss them and make suggestions for further research.The reviewed papers are presented on the basis of the various models for pavement-vehicle systems and the various methods for dynamically analyzing these systems.Flexible pavements are modeled by a homogeneous or layered half-plane with isotropic or anisotropic and linear elastic,viscoelastic or poroelastic material behavior.Rigid pavements are modeled by a beam or plate on a homogeneous or layered half-plane with material properties like the ones for flexible pavements.The vehicles are modeled as concentrated or distributed over a finite area loads moving with constant or time dependent speed.The above pavement-vehicle models are dynamically analyzed by analytical,analytical/numerical or purely numerical methods working in the time or frequency domain.Representative examples are presented to illustrate the models and methods of analysis,demonstrate their merits and assess the effects of the various parameters on pavement response.The paper closes with con-clusions and suggestions for further research in the area.The significance of this research effort has to do with the presentation of the existing literature on the subject in a critical and easy to understand way with the aid of representative examples and the identification of new research areas.

Study on the Evolution of Plane Shape of Guangfu Ancestral Hall Buildings and the Construction of Type Genealogy in the Ming and Qing Dynasties

This paper studies the clan temples in four periods,namely,the development period(early Ming Dynasty and before),the shaping period(mid-Ming Dynasty to early Qing Dynasty),the maturity period(mid-Qing Dynasty),and the programmed period(late Qing Dynasty),analyzes and summarizes their prototypical planform features and attempts to summarize and construct a genealogy of the architectural characteristics of the plane shape of clan temples.During the Ming and Qing dynasties,Guangfu ancestral halls in general showed a journey of development from a centripetal-single-linear-multi-linear-grid pattern,forming the unique Guangsanlu classic pattern in the Guangfu area.

Erector Spinae Plane Block Combined with Serratus Anterior Plane Block Versus Thoracic Paravertebral Block for Postoperative Analgesia and Recovery After Thoracoscopic Surgery:A Randomized Controlled Non-inferiority Clinical Trial

Objective This study aimed to compare the postoperative analgesia and recovery of ultrasound-guided erector spinae plane block combined with serratus anterior plane block(ESPB combined with SAPB)versus thoracic paravertebral block(PVB)after thoracoscopic surgery.Methods Ninety-two patients who underwent video-assisted thoracoscopic surgery(VATS)were randomly divided into group S(n=46)and group P(n=46).After anesthesia induction,the same anesthesiologist performed ultrasound-guided ESPB at T5 and T7 levels combined with SAPB at the level of the fifth rib in the midaxillary line in group S and ultrasound-guided PVB at T5 and T7 levels in group P.Patients in both groups were given 40 mL of 0.4%ropivacaine.Eighty-six patients completed the study(group S,n=44;group P,n=42).The morphine consumption,visual analogue scale(VAS)scores at rest and coughing,and frequency of remedial analgesia were recorded at 1,2,4,8,and 24 h postoperatively.Pulmonary function parameters were recorded at 1,4,and 24 h postoperatively,and the quality of recovery(QoR)-15 score at 24 h postoperatively.The adverse effects,duration of chest tube drainage and length of stay were also recorded.Results The morphine consumption at postoperative 4 and 8 h and the incidence of ipsilateral shoulder pain(ISP)were significantly lower in group S than in group P.The QoR-15 questionnaire score at postoperative 24 h was significantly lower in group P than in group S(P<0.05).The morphine consumption was lower at 24 h postoperatively in group S than in group P,with no significant difference found yet.The morphine consumption at other observed times,VAS scores,pulmonary function parameters,frequency of remedial analgesia,duration of chest tube drainage,length of stay,and incidence of other adverse events were comparable between group S and group P.Conclusion Ultrasound-guided ESPB combined with SAPB is non-inferior to PVB in terms of morphine consumption at postoperative 24 h and postoperative recovery.But,this approach can significantly reduce morphine consumption in the early postoperative period(0–8 h)after thoracoscopy with lower incidence of ISP.It is a simpler and safer operation.

Static and free vibration analysis of four-parameter continuous grading elliptical sandwich plates

In the present study,the static and dynamic analyses of elliptical functionally graded sandwich(FGS)plates are investigated.The constituent materials of the sandwich plates are ceramic and metal so that the core is made of pure metal,while the face sheets consist of a combination of metal and ceramic according to a four-parameter power-law distribution.Different material profiles such as classic,symmetric,and asymmetric can be obtained using the applied generalized power-law distribution relation.The analysis is performed based on the classical laminated plate theory(CLPT)and the Ritz method.The effects of four parameters in the material distribution relation as well as different geometric parameters on the deflection and natural frequencies of elliptical FGS plates are studied.The results of this study show that with a proper distribution of materials,the optimal static and dynamic behavior can be achieved.The results also indicate that the generalized power-law distribution has significant effects on the natural frequencies of elliptical FGS plates.For example,although the frequency parameter of a plate with ceramic face sheets is more than the one with metal face sheets,the use of larger amounts of ceramic does not necessarily increase the natural frequency of the structure.

Crystal plane induced in-situ electrochemical activation of manganese-based cathode enable long-term aqueous zinc-ion batteries

Rapid capacity decay and sluggish reaction kinetics are major barriers hindering the applications of manganese-based cathode materials for aqueous zinc-ion batteries.Herein,the effects of crystal plane on the in-situ transformation behavior and electrochemical performance of manganese-based cathode is discussed.A comprehensive discussion manifests that the exposed(100)crystal plane is beneficial to the phase transformation from tunnel-structured MnO_(2) to layer-structured ZnMn_(3)O_(7)·3H_(2)O,which plays a critical role for the high reactivity,high capacity,fast diffusion kinetics and long cycling stability.Additionally,a two-stage zinc storage mechanism can be demonstrated,involving continuous activation reaction and phase transition reaction.As expected,it exhibits a high capacity of 275 mAh g^(-1)at 100 mA g^(-1),a superior durability over 1000 cycles and good rate capability.This study may open new windows toward developing advanced cathodes for ZIBs,and facilitate the applications of ZIBs in large-scale energy storage system.

Anti-Corrosion and Reconstruction of Surface Crystal Plane for Zn Anodes by an Advanced Metal Passivation Technique

For the aqueous Zn-ion battery,dendrite formation,corrosion,and interfacial parasitic reactions are major issues,which greatly inhibits their practical application.How to develop a method of Zn construction or treatment to solve these issues for Zn anodes are still great challenges.Herein,a simple and cheap metal passivation technique is proposed for Zn anodes from a corrosion science perspective.Similar to the metal anticorrosion engineering,the formed interfacial protective layer in a chemical way can sufficiently solve the corrosion issues.Furthermore,the proposed passivity approach can reconstruct Zn surface-preferred crystal planes,exposing more(002)planes and improving surface hydrophilicity,which inhibits the formation of Zn dendrites and hydrogen evolution effectively.As expected,the passivated Zn achieves outstanding cycling life(1914 h)with low voltage polarization(<40 mV).Even at 6 mA cm^(−2) and 3 mA h cm^(−2),it can achieve stable Zn deposition over 460 h.The treated Zn anode coupled with MnO_(2) cathode shows prominently reinforced full batteries service life,making it a potential Zn anode candidate for excellent performance aqueous Zn-ion batteries.The proposed passivation approach provides a guideline for other metal electrodes preparation in various batteries and establishes the connections between corrosion science and batteries.

Event plane determination from the zero degree calorimeter at the cooling storage ring external-target experiment

The Cooling Storage Ring external-target experiment(CEE)spectrometer is used to study the nuclear matter created in heavy-ion collisions at√sNN=2.1-2.4 GeV with the aim to reveal the quantum chromodynamics phase structure in the high-baryon-density region.Collective flow is considered an effective probe for evaluating the properties of media during high-energy nuclear collisions.One of the main functions of the zero-degree calorimeter(ZDC),a subdetector system in the CEE,is to determine the reaction plane in heavy-ion collisions.This step is crucial for measuring the collective flow and other reaction-plane-related analyses.In this paper,we illustrate the procedures for event plane determination using the ZDC.Finally,isospin-dependent quantum molecular dynamics model-based predictions of the rapidity dependence of the directed and elliptical flows for p,d,t,3He,and 4He,produced in 2.1 GeV U+U collisions,are presented.

Dynamic damage evolution of bank slopes with serrated structural planes considering the deteriorated rock mass and frequent reservoirinduced earthquakes

To investigate the dynamic damage evolution characteristics of bank slopes with serrated structural planes,the shaking table model test and the numerical simulation were utilized.The main findings indicate that under continuous seismic loads,the deformation of the bank slope increased,particularly around the hydro-fluctuation belt,accompanying by the pore water pressure rising.The soil pressure increased and then decreased showed dynamic variation characteristics.As the undulation angle of the serrated structural planes increased(30°, 45°, and 60°),the failure modes were climbing,climbinggnawing,and gnawing respectively.The first-order natural frequency was used to calculate the damage degree(Dd)of the bank slope.During microseisms and small earthquakes,it was discovered that the evolution of Dd followed the“S”shape,which was fitted by a logic function.Additionally,the quadratic function was used to fit the Dd during moderately strong earthquakes.Through the numerical simulation,the variation characteristics of safety factors(Sf)for slopes with serrated structural planes and slopes with straight structural planes were compared.Under continuous seismic loads,the Sf of slopes with straight structural planes reduce stalely,whereas the Sf for slopes with serrated structural planes was greater than the former and the reduction rate was increasing.

Position Optimization for Mobile Relay Communications with Constrained Energy over a Whole Horizontal Plane

This paper concerns the position optimization problem of a mobile relay over a whole horizontal plane.This problem is important because the position of a mobile relay directly affects the end-to-end performance,e.g.,reliability,connectivity,and data rate.In this paper,we propose a new position optimization scheme of a mobile relay over a whole horizontal plane based on the one-bit feedback information from the destination node,which improves the performance over the prior scheme whose position of the mobile relay is optimized over a fixed orbit.In the proposed scheme,the mobile relay is equipped with merely one single onboard antenna.Moreover,no prior information about the positions of both the source node and the destination node is required.Thus,the proposed scheme can work at low network resources scenario,which is particularly suitable for mobile relay communication with constrained energy,e.g.,the communications in a disaster area where the infrastructure is heavily damaged,volcano monitoring,and wireless powered communication networking.According to the characteristics of the proposed scheme,we further design two heuristic implementations to calculate the optimal position of a mobile relay over a whole horizontal plane.The first implementation has better steady performance whereas the second implementation has better convergence speed.We implement the proposed scheme and conduct an extensive performance comparison between the proposed scheme and prior schemes to verify the advantages of the proposed scheme.

Study of beam divergence and thrust vector eccentricity characteristics of the Hall thruster based on dual Faraday probe array planes and its applications

The accurate knowledge of the thrust vector eccentricity and beam divergence characteristics of Hall thrusters are of significant engineering value for the beneficial integration and successful application of Hall thrusters on spacecraft.For the characteristics of the plume bipolar diffusion due to the annular discharge channel of the Hall thruster,a Gaussian-fitted method for thrust vector deviation angle and beam divergence of Hall thrusters based on dual Faraday probe array planes was proposed in respect of the Hall thruster beam characteristics.The results show that the ratios of the deviation between the maximum and minimum values of the beam divergence angle and the thrust vector eccentricity angle using a Gaussian fit to the optimized Faraday probe dual plane to the mean value are 1.4%and 11.5%,respectively.The optimized thrust vector eccentricity angle obtained has been substantially improved,by approximately 20%.The beam divergence angle calculated using a Gaussian fitting to the optimized Faraday probe dual plane is approximately identical to the non-optimized one.The beam divergence and thrust vector eccentricity angles for different anode mass flow rates were obtained by averaging the beam divergence and thrust vector eccentricity angles calculated by the dual-plane,Gaussian-fitted ion current density method for different cross-sections.The study not only allows for an immediate and effective tool for determining the design of thrust vector adjustment mechanisms of spacecraft with different power Hall thrusters but also for characterizing the 3D spatial distribution of the Hall thruster plume.

Two-dimensional plane strain consolidation of unsaturated soils considering the depth-dependent stress

In practical engineering,the total vertical stress in the soil layer is not constant due to stress diffusion,and varies with time and depth.Therefore,the purpose of this paper is to investigate the effect of stress diffusion on the two-dimensional(2D)plane strain consolidation properties of unsaturated soils when the stress varies with time and depth.A series of semi-analytical solutions in terms of excess pore air and water pressures and settlement for 2D plane strain consolidation of unsaturated soils can be derived with the joint use of Laplace transform and Fourier sine series expansion.Then,the inverse Laplace transform of the semi-analytical solution is given in the time domain using a self-programmed code based on Crump’s method.The reliability of the obtained solutions is proved by the degeneration.Finally,the 2D plots of excess pore pressures and the curves of settlement varying with time,considering different physical parameters of unsaturated soil stratum and depth-dependent stress,are depicted and analyzed to study the 2D plane strain consolidation properties of unsaturated soils subjected to the depthdependent stress.