The phenomenology involved in severe accidents in nuclear reactors is highly complex.Currently,integrated analysis programs used for severe accident analysis heavily rely on custom empirical parameters,which introduce...The phenomenology involved in severe accidents in nuclear reactors is highly complex.Currently,integrated analysis programs used for severe accident analysis heavily rely on custom empirical parameters,which introduce considerable uncertainty.Therefore,in recent years,the field of severe accidents has shifted its focus toward applying uncertainty analysis methods to quantify uncertainty in safety assessment programs,known as“best estimate plus uncertainty(BEPU).”This approach aids in enhancing our comprehension of these programs and their further development and improvement.This study concentrates on a third-generation pressurized water reactor equipped with advanced active and passive mitigation strategies.Through an Integrated Severe Accident Analysis Program(ISAA),numerical modeling and uncertainty analysis were conducted on severe accidents resulting from large break loss of coolant accidents.Seventeen uncertainty parameters of the ISAA program were meticulously screened.Using Wilks'formula,the developed uncertainty program code,SAUP,was employed to carry out Latin hypercube sampling,while ISAA was employed to execute batch calculations.Statistical analysis was then conducted on two figures of merit,namely hydrogen generation and the release of fission products within the pressure vessel.Uncertainty calculations revealed that hydrogen production and the fraction of fission product released exhibited a normal distribution,ranging from 182.784 to 330.664 kg and from 15.6 to 84.3%,respectively.The ratio of hydrogen production to reactor thermal power fell within the range of 0.0578–0.105.A sensitivity analysis was performed for uncertain input parameters,revealing significant correlations between the failure temperature of the cladding oxide layer,maximum melt flow rate,size of the particulate debris,and porosity of the debris with both hydrogen generation and the release of fission products.展开更多
Within this work,we perform a sensitivity analysis to determine the influence of the material input parameters on the pressure in an isotropic porous solid cylinder.We provide a step-by-step guide to obtain the analyt...Within this work,we perform a sensitivity analysis to determine the influence of the material input parameters on the pressure in an isotropic porous solid cylinder.We provide a step-by-step guide to obtain the analytical solution for a porous isotropic elastic cylinder in terms of the pressure,stresses,and elastic displacement.We obtain the solution by performing a Laplace transform on the governing equations,which are those of Biot's poroelasticity in cylindrical polar coordinates.We enforce radial boundary conditions and obtain the solution in the Laplace transformed domain before reverting back to the time domain.The sensitivity analysis is then carried out,considering only the derived pressure solution.This analysis finds that the time t,Biot's modulus M,and Poisson's ratio ν have the highest influence on the pressure whereas the initial value of pressure P_(0) plays a very little role.展开更多
Traditional global sensitivity analysis(GSA)neglects the epistemic uncertainties associated with the probabilistic characteristics(i.e.type of distribution type and its parameters)of input rock properties emanating du...Traditional global sensitivity analysis(GSA)neglects the epistemic uncertainties associated with the probabilistic characteristics(i.e.type of distribution type and its parameters)of input rock properties emanating due to the small size of datasets while mapping the relative importance of properties to the model response.This paper proposes an augmented Bayesian multi-model inference(BMMI)coupled with GSA methodology(BMMI-GSA)to address this issue by estimating the imprecision in the momentindependent sensitivity indices of rock structures arising from the small size of input data.The methodology employs BMMI to quantify the epistemic uncertainties associated with model type and parameters of input properties.The estimated uncertainties are propagated in estimating imprecision in moment-independent Borgonovo’s indices by employing a reweighting approach on candidate probabilistic models.The proposed methodology is showcased for a rock slope prone to stress-controlled failure in the Himalayan region of India.The proposed methodology was superior to the conventional GSA(neglects all epistemic uncertainties)and Bayesian coupled GSA(B-GSA)(neglects model uncertainty)due to its capability to incorporate the uncertainties in both model type and parameters of properties.Imprecise Borgonovo’s indices estimated via proposed methodology provide the confidence intervals of the sensitivity indices instead of their fixed-point estimates,which makes the user more informed in the data collection efforts.Analyses performed with the varying sample sizes suggested that the uncertainties in sensitivity indices reduce significantly with the increasing sample sizes.The accurate importance ranking of properties was only possible via samples of large sizes.Further,the impact of the prior knowledge in terms of prior ranges and distributions was significant;hence,any related assumption should be made carefully.展开更多
Battery production is crucial for determining the quality of electrode,which in turn affects the manufactured battery performance.As battery production is complicated with strongly coupled intermediate and control par...Battery production is crucial for determining the quality of electrode,which in turn affects the manufactured battery performance.As battery production is complicated with strongly coupled intermediate and control parameters,an efficient solution that can perform a reliable sensitivity analysis of the production terms of interest and forecast key battery properties in the early production phase is urgently required.This paper performs detailed sensitivity analysis of key production terms on determining the properties of manufactured battery electrode via advanced data-driven modelling.To be specific,an explainable neural network named generalized additive model with structured interaction(GAM-SI)is designed to predict two key battery properties,including electrode mass loading and porosity,while the effects of four early production terms on manufactured batteries are explained and analysed.The experimental results reveal that the proposed method is able to accurately predict battery electrode properties in the mixing and coating stages.In addition,the importance ratio ranking,global interpretation and local interpretation of both the main effects and pairwise interactions can be effectively visualized by the designed neural network.Due to the merits of interpretability,the proposed GAM-SI can help engineers gain important insights for understanding complicated production behavior,further benefitting smart battery production.展开更多
This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are...This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are elucidated geometrically from the perspective of expanding ellipsoids.Based on this geometric interpretation,the QFOSM is further extended to estimate sensitivity indices and assess the significance of various uncertain parameters involved in the slope system.The proposed method has the advantage of computational simplicity,akin to the conventional first-order second-moment method(FOSM),while providing estimation accuracy close to that of the first-order reliability method(FORM).Its performance is demonstrated with a numerical example and three slope examples.The results show that the proposed method can efficiently estimate the slope reliability and simultaneously evaluate the sensitivity of the uncertain parameters.The proposed method does not involve complex optimization or iteration required by the FORM.It can provide a valuable complement to the existing approximate reliability analysis methods,offering rapid sensitivity evaluation and slope reliability analysis.展开更多
This paper presents an engineering system approach using a 2D model of conservation of mass to study the dynamics of ozone and concerned chemical species in the stratosphere.By considering all fourteen photolysis,ozon...This paper presents an engineering system approach using a 2D model of conservation of mass to study the dynamics of ozone and concerned chemical species in the stratosphere.By considering all fourteen photolysis,ozone-generating,and-depleting chemical reactions,the model calculated the transient,spatial changes of ozone under different physical-chemical-radiative conditions.Validation against the measured data demonstrated good accuracy,close match of our model with the observed ozone concentrations at both 20°S and 90°N locations.The deviation in the average concentration was less than 1% and in ozone profiles less than 17%.The impacts of various chlorine-(Cl),nitrogen oxides-(NO_(x)),and bromine-(Br)depleting cycles on ozone concentrations and distribution were investigated.The chlorine catalytic depleting cycle was found to exhibit the most significant impact on ozone dynamics,confirming the key role of chlorine in the problem of ozone depletion.Sensitivity analysis was conducted with levels of 25%,50%,100%,200%,and 400% of the baseline value.The combined cycles(Cl+NO_(x)+Br)showed the most significant influence on ozone behavior.The total ozone abundance above the South Pole could decrease by a small 3%,from 281 DU(Dubson Units)to 273 DU for the 25% level,or by a huge thinning of 60%to 114 DU for the 400% concentration level.When the level of chlorine gases increased beyond 200%,it would cause ozone depletion to a level of ozone hole(below 220 DU).The 2D Ozone Model presented in this paper demonstrates robustness,convenience,efficiency,and executability for analyzing complex ozone phenomena in the stratosphere.展开更多
Case history investigations have shown that pile foundations are more critically damaged in liquefiable soils than non-liquefiable soils.This study examines the differences in seismic response of pile foundations in l...Case history investigations have shown that pile foundations are more critically damaged in liquefiable soils than non-liquefiable soils.This study examines the differences in seismic response of pile foundations in liquefiable and non-liquefiable soils and their sensitivity to numerical model parameters.A two-dimensional finite element(FE)model is developed to simulate the experiment of a single pile foundation centrifuge in liquefiable soil subjected to earthquake motions and is validated against real-world test results.The differences in soil-pile seismic response of liquefiable and non-liquefiable soils are explored.Specifically,the first-order second-moment method(FOSM)is used for sensitivity analysis of the seismic response.The results show significant differences in seismic response for a soil-pile system between liquefiable and non-liquefiable soil.The seismic responses are found to be significantly larger in liquefiable soil than in non-liquefiable soil.Moreover,the pile bending moment was mainly affected by the kinematic effect in liquefiable soil,while the inertial effect was more significant in non-liquefiable soil.The controlling parameters of seismic response were PGA,soil density,and friction angle in liquefiable soil,while the pile bending moment was mainly controlled by PGA,the friction angle of soil,and shear modulus of loose sand in non-liquefiable soil.展开更多
Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately ...Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately predict the productivity of hydrate-bearing sediments(HBS).The multi-phase seepage parameters of HBS include permeability,porosity,which is closely related to permeability,and hydrate saturation,which has a direct impact on hydrate content.Existing research has shown that these multi-phase seepage parameters have a great impact on HBS productivity.Permeability directly affects the transmission of pressure-drop and discharge of methane gas,porosity and initial hydrate saturation affect the amount of hydrate decomposition and transmission process of pressure-drop,and also indirectly affect temperature variation of the reservoir.Considering the spatial heterogeneity of multi-phase seepage parameters,a depressurization production model with layered heterogeneity is established based on the clayey silt hydrate reservoir at W11 station in the Shenhu Sea area of the South China Sea.Tough+Hydrate software was used to calculate the production model;the process of gas production and seepage parameter evolution under different multi-phase seepage conditions were obtained.A sensitivity analysis of the parameters affecting the reservoir productivity was conducted so that:(a)a HBS model with layered heterogeneity can better describe the transmission process of pressure and thermal compensation mechanism of hydrate reservoir;(b)considering the multi-phase seepage parameter heterogeneity,the influence degrees of the parameters on HBS productivity were permeability,porosity and initial hydrate saturation,in order from large to small,and the influence of permeability was significantly greater than that of other parameters;(c)the production potential of the clayey silt reservoir should not only be determined by hydrate content or seepage capacity,but also by the comprehensive effect of the two;and(d)time scales need to be considered when studying the effects of changes in multi-phase seepage parameters on HBS productivity.展开更多
Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to ...Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.展开更多
Liquefaction is one of the most destructive phenomena caused by earthquakes,which has been studied in the issues of potential,triggering and hazard analysis.The strain energy approach is a common method to investigate...Liquefaction is one of the most destructive phenomena caused by earthquakes,which has been studied in the issues of potential,triggering and hazard analysis.The strain energy approach is a common method to investigate liquefaction potential.In this study,two Artificial Neural Network(ANN)models were developed to estimate the liquefaction resistance of sandy soil based on the capacity strain energy concept(W)by using laboratory test data.A large database was collected from the literature.One group of the dataset was utilized for validating the process in order to prevent overtraining the presented model.To investigate the complex influence of fine content(FC)on liquefaction resistance,according to previous studies,the second database was arranged by samples with FC of less than 28%and was used to train the second ANN model.Then,two presented ANN models in this study,in addition to four extra available models,were applied to an additional 20 new samples for comparing their results to show the capability and accuracy of the presented models herein.Furthermore,a parametric sensitivity analysis was performed through Monte Carlo Simulation(MCS)to evaluate the effects of parameters and their uncertainties on the liquefaction resistance of soils.According to the results,the developed models provide a higher accuracy prediction performance than the previously publishedmodels.The sensitivity analysis illustrated that the uncertainties of grading parameters significantly affect the liquefaction resistance of soils.展开更多
For laser cladding a large temperature gradient easily weakened the surface quality by generating cracks and irregular coating surfaces,which in turn affected the bearing capacity and corrosion resistance of coatings ...For laser cladding a large temperature gradient easily weakened the surface quality by generating cracks and irregular coating surfaces,which in turn affected the bearing capacity and corrosion resistance of coatings in the rapid heating and cooling process.The response surface methodology(RSM)was used to predict coating cracks by changing the powder ratio,energy density,and preheating temperature,which obtained the relevant mathematical model.After that,the sensitivity of the crack length to process parameters was analyzed based on the sensitivity analysis method.The effect of Ni60/WC composite powder process parameters on the surface quality was revealed in laser cladding.The crack length first decreased and then increased,and the Smooth decreased with the increased powder ratio.The crack length and Smooth increased with the increased energy density.The crack length decreased and Smooth increased with the increased preheating temperature.Sensitivity analysis showed that the crack length and Smooth were the most sensitive to the powder ratio.Therefore,the process parameters were reasonably selected to control the surface quality.The mathematical model and sensitivity analysis method in the work could improve the surface quality,which provided a theoretical basis for the prediction and control of laser cladding cracks.展开更多
Simulations and predictions using numerical models show considerable uncertainties,and parameter uncertainty is one of the most important sources.It is impractical to improve the simulation and prediction abilities by...Simulations and predictions using numerical models show considerable uncertainties,and parameter uncertainty is one of the most important sources.It is impractical to improve the simulation and prediction abilities by reducing the uncertainties of all parameters.Therefore,identifying the sensitive parameters or parameter combinations is crucial.This study proposes a novel approach:conditional nonlinear optimal perturbations sensitivity analysis(CNOPSA)method.The CNOPSA method fully considers the nonlinear synergistic effects of parameters in the whole parameter space and quantitatively estimates the maximum effects of parameter uncertainties,prone to extreme events.Results of the analytical g-function test indicate that the CNOPSA method can effectively identify the sensitivity of variables.Numerical results of the theoretical five-variable grassland ecosystem model show that the maximum influence of the simulated wilted biomass caused by parameter uncertainty can be estimated and computed by employing the CNOPSA method.The identified sensitive parameters can easily change the simulation or prediction of the wilted biomass,which affects the transformation of the grassland state in the grassland ecosystem.The variance-based approach may underestimate the parameter sensitivity because it only considers the influence of limited parameter samples from a statistical view.This study verifies that the CNOPSA method is effective and feasible for exploring the important and sensitive physical parameters or parameter combinations in numerical models.展开更多
The structural response of a single-layer reticulated dome to external explosions is shaped by many variables,and the associated uncertainties imply non-deterministic results.Existing deterministic methods for predict...The structural response of a single-layer reticulated dome to external explosions is shaped by many variables,and the associated uncertainties imply non-deterministic results.Existing deterministic methods for predicting the consequences of specific explosions do not account for these uncertainties.Therefore,the impact of the uncertainties associated with these input variables on the structures’response needs to be studied and quantified.In this study,a parametric uncertainty analysis was conducted first.Then,local and global sensitivity analyses were carried out to identify the drivers of the structural dynamic response.A probabilistic structural response model was established based on sensitive variables and a reasonable sample size.Furthermore,some deterministic empirical methods for explosion-resistance design,including the plane blast load model of CONWEP,the curved blast load model under the 50%assurance level,and the 20%mass-increased method,were used for evaluating their reliability.The results of the analyses revealed that the structural response of a single-layer reticulated dome to an external blast loading is lognormally distributed.Evidently,the MB0.5 method based on the curved reflector load model yielded results with a relatively stable assurance rate and reliability,but CONWEP did not;thus,the 1.2MB0.5 method can be used for making high-confidence simple predictions.In addition,the results indicated that the structural response is very sensitive to the explosion parameters.Based on these results,it is suggested that for explosion proofing,setting up a defensive barrier is more effective than structural strengthening.展开更多
The anti-sliding stability of a gravity dam along its foundation surface is a key problem in the design of gravity dams.In this study,a sensitivity analysis framework was proposed for investigating the factors affecti...The anti-sliding stability of a gravity dam along its foundation surface is a key problem in the design of gravity dams.In this study,a sensitivity analysis framework was proposed for investigating the factors affecting gravity dam anti-sliding stability along the foundation surface.According to the design specifications,the loads and factors affecting the stability of a gravity dam were comprehensively selected.Afterwards,the sensitivity of the factors was preliminarily analyzed using the Sobol method with Latin hypercube sampling.Then,the results of the sensitivity analysis were verified with those obtained using the Garson method.Finally,the effects of different sampling methods,probability distribution types of factor samples,and ranges of factor values on the analysis results were evaluated.A case study of a typical gravity dam in Yunnan Province of China showed that the dominant factors affecting the gravity dam anti-sliding stability were the anti-shear cohesion,upstream and downstream water levels,anti-shear friction coefficient,uplift pressure reduction coefficient,concrete density,and silt height.Choice of sampling methods showed no significant effect,but the probability distribution type and the range of factor values greatly affected the analysis results.Therefore,these two elements should be sufficiently considered to improve the reliability of the dam anti-sliding stability analysis.展开更多
Aiming at optimizing the energy consumption of HVAC,an energy conservation optimization method was proposed for HVAC systems based on the sensitivity analysis(SA),named the sensitivity analysis combination method(SAC)...Aiming at optimizing the energy consumption of HVAC,an energy conservation optimization method was proposed for HVAC systems based on the sensitivity analysis(SA),named the sensitivity analysis combination method(SAC).Based on the SA,neural network and the related settings about energy conservation of HVAC systems,such as cooling water temperature,chilled water temperature and supply air temperature,were optimized.Moreover,based on the data of the existing HVAC system,various optimal control methods ofHVAC systems were tested and evaluated by a simulated HVAC system in TRNSYS.The results show that the proposed SA combination method can reduce significant computational load while maintaining an equivalent energy performance compared with traditional optimal control methods.展开更多
A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed...A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed by stoma in leaves. In the photosynthesis module of these LSMs, variations of parameters arising from diversity in plant functional types(PFTs) and climate remain unclear. Identifying sensitive parameters among all photosynthetic parameters before parameter estimation can not only reduce operation cost, but also improve the usability of photosynthesis models worldwide. Here, we analyzed 13 parameters of a biochemically-based photosynthesis model(FvCB), implemented in many LSMs, using two sensitivity analysis(SA) methods(i.e., the Sobol’ method and the Morris method) for setting up the parameter ensemble. Three different model performance metrics, i.e.,Root Mean Squared Error(RMSE), Nash Sutcliffe efficiency(NSE), and Standard Deviation(STDEV) were introduced for model assessment and sensitive parameters identification. The results showed that among all photosynthetic parameters only a small portion of parameters were sensitive, and the sensitive parameters were different across plant functional types: maximum rate of Rubisco activity(Vcmax25), maximum electron transport rate(Jmax25), triose phosphate use rate(TPU) and dark respiration in light(Rd) were sensitive in broad leafevergreen trees(BET), broad leaf-deciduous trees(BDT) and needle leaf-evergreen trees(NET), while only Vcmax25and TPU are sensitive in short vegetation(SV), dwarf trees and shrubs(DTS), and agriculture and grassland(AG). The two sensitivity analysis methods suggested a strong SA coherence;in contrast, different model performance metrics led to different SA results. This misfit suggests that more accurate values of sensitive parameters, specifically, species specific and seasonal variable parameters, are required to improve the performance of the FvCB model.展开更多
The prediction of processor performance has important referencesignificance for future processors. Both the accuracy and rationality of theprediction results are required. The hierarchical belief rule base (HBRB)can i...The prediction of processor performance has important referencesignificance for future processors. Both the accuracy and rationality of theprediction results are required. The hierarchical belief rule base (HBRB)can initially provide a solution to low prediction accuracy. However, theinterpretability of the model and the traceability of the results still warrantfurther investigation. Therefore, a processor performance prediction methodbased on interpretable hierarchical belief rule base (HBRB-I) and globalsensitivity analysis (GSA) is proposed. The method can yield more reliableprediction results. Evidence reasoning (ER) is firstly used to evaluate thehistorical data of the processor, followed by a performance prediction modelwith interpretability constraints that is constructed based on HBRB-I. Then,the whale optimization algorithm (WOA) is used to optimize the parameters.Furthermore, to test the interpretability of the performance predictionprocess, GSA is used to analyze the relationship between the input and thepredicted output indicators. Finally, based on the UCI database processordataset, the effectiveness and superiority of the method are verified. Accordingto our experiments, our prediction method generates more reliable andaccurate estimations than traditional models.展开更多
The absorber is the key unit in the post-combustion monoethanolamine(MEA)-based carbon dioxide(CO_(2))capture process.A rate-based dynamic model for the absorber is developed and validated using steady-state experimen...The absorber is the key unit in the post-combustion monoethanolamine(MEA)-based carbon dioxide(CO_(2))capture process.A rate-based dynamic model for the absorber is developed and validated using steady-state experimental data reported in open literature.Sensitivity analysis is performed with respect to important model parameters associated with the reaction,mass transport and phy-sical property relationships.Then,a singular value decomposition(SVD)-based subspace parameter estimation method is proposed to improve the model accu-racy.Finally,dynamic simulations are carried out to investigate the effects of the feed rate of lean MEA solution and the flue inlet conditions.Simulation results indicate that the established dynamic model can reasonably reflect the physical behavior of the absorber.Some new insights are gained from the simulation results.展开更多
Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensi...Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.展开更多
Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable s...Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable stress sensitivity characterization models is still limited.In this study,three commonly used stress sensitivity models for shale oil reservoirs were considered,and experiments on representative core samples were conducted.By fitting and comparing the data,the“exponential model”was identified as a characterization model that accurately represents stress sensitivity in continental shale oil reservoirs.To validate the accuracy of the model,a two-phase seepage mathematical model for shale oil reservoirs coupled with the exponential model was introduced.The model was discretely solved using the finite volume method,and its accuracy was verified through the commercial simulator CMG.The study evaluated the productivity of a typical horizontal well under different engineering,geological,and fracture conditions.The results indicate that considering stress sensitivity leads to a 13.57%reduction in production for the same matrix permeability.Additionally,as the fracture half-length and the number of fractures increase,and the bottomhole flowing pressure decreases,the reservoir stress sensitivity becomes higher.展开更多
基金This work was supported financially by the National Natural Science Foundation of China(No.12375176).
文摘The phenomenology involved in severe accidents in nuclear reactors is highly complex.Currently,integrated analysis programs used for severe accident analysis heavily rely on custom empirical parameters,which introduce considerable uncertainty.Therefore,in recent years,the field of severe accidents has shifted its focus toward applying uncertainty analysis methods to quantify uncertainty in safety assessment programs,known as“best estimate plus uncertainty(BEPU).”This approach aids in enhancing our comprehension of these programs and their further development and improvement.This study concentrates on a third-generation pressurized water reactor equipped with advanced active and passive mitigation strategies.Through an Integrated Severe Accident Analysis Program(ISAA),numerical modeling and uncertainty analysis were conducted on severe accidents resulting from large break loss of coolant accidents.Seventeen uncertainty parameters of the ISAA program were meticulously screened.Using Wilks'formula,the developed uncertainty program code,SAUP,was employed to carry out Latin hypercube sampling,while ISAA was employed to execute batch calculations.Statistical analysis was then conducted on two figures of merit,namely hydrogen generation and the release of fission products within the pressure vessel.Uncertainty calculations revealed that hydrogen production and the fraction of fission product released exhibited a normal distribution,ranging from 182.784 to 330.664 kg and from 15.6 to 84.3%,respectively.The ratio of hydrogen production to reactor thermal power fell within the range of 0.0578–0.105.A sensitivity analysis was performed for uncertain input parameters,revealing significant correlations between the failure temperature of the cladding oxide layer,maximum melt flow rate,size of the particulate debris,and porosity of the debris with both hydrogen generation and the release of fission products.
基金Project supported by the Engineering and Physical Sciences Research Council of U. K.(Nos. EP/S030875/1, EP/T017899/1, and EP/T517896/1)。
文摘Within this work,we perform a sensitivity analysis to determine the influence of the material input parameters on the pressure in an isotropic porous solid cylinder.We provide a step-by-step guide to obtain the analytical solution for a porous isotropic elastic cylinder in terms of the pressure,stresses,and elastic displacement.We obtain the solution by performing a Laplace transform on the governing equations,which are those of Biot's poroelasticity in cylindrical polar coordinates.We enforce radial boundary conditions and obtain the solution in the Laplace transformed domain before reverting back to the time domain.The sensitivity analysis is then carried out,considering only the derived pressure solution.This analysis finds that the time t,Biot's modulus M,and Poisson's ratio ν have the highest influence on the pressure whereas the initial value of pressure P_(0) plays a very little role.
文摘Traditional global sensitivity analysis(GSA)neglects the epistemic uncertainties associated with the probabilistic characteristics(i.e.type of distribution type and its parameters)of input rock properties emanating due to the small size of datasets while mapping the relative importance of properties to the model response.This paper proposes an augmented Bayesian multi-model inference(BMMI)coupled with GSA methodology(BMMI-GSA)to address this issue by estimating the imprecision in the momentindependent sensitivity indices of rock structures arising from the small size of input data.The methodology employs BMMI to quantify the epistemic uncertainties associated with model type and parameters of input properties.The estimated uncertainties are propagated in estimating imprecision in moment-independent Borgonovo’s indices by employing a reweighting approach on candidate probabilistic models.The proposed methodology is showcased for a rock slope prone to stress-controlled failure in the Himalayan region of India.The proposed methodology was superior to the conventional GSA(neglects all epistemic uncertainties)and Bayesian coupled GSA(B-GSA)(neglects model uncertainty)due to its capability to incorporate the uncertainties in both model type and parameters of properties.Imprecise Borgonovo’s indices estimated via proposed methodology provide the confidence intervals of the sensitivity indices instead of their fixed-point estimates,which makes the user more informed in the data collection efforts.Analyses performed with the varying sample sizes suggested that the uncertainties in sensitivity indices reduce significantly with the increasing sample sizes.The accurate importance ranking of properties was only possible via samples of large sizes.Further,the impact of the prior knowledge in terms of prior ranges and distributions was significant;hence,any related assumption should be made carefully.
基金supported by the National Natural Science Foundation of China (62373224,62333013,U23A20327)。
文摘Battery production is crucial for determining the quality of electrode,which in turn affects the manufactured battery performance.As battery production is complicated with strongly coupled intermediate and control parameters,an efficient solution that can perform a reliable sensitivity analysis of the production terms of interest and forecast key battery properties in the early production phase is urgently required.This paper performs detailed sensitivity analysis of key production terms on determining the properties of manufactured battery electrode via advanced data-driven modelling.To be specific,an explainable neural network named generalized additive model with structured interaction(GAM-SI)is designed to predict two key battery properties,including electrode mass loading and porosity,while the effects of four early production terms on manufactured batteries are explained and analysed.The experimental results reveal that the proposed method is able to accurately predict battery electrode properties in the mixing and coating stages.In addition,the importance ratio ranking,global interpretation and local interpretation of both the main effects and pairwise interactions can be effectively visualized by the designed neural network.Due to the merits of interpretability,the proposed GAM-SI can help engineers gain important insights for understanding complicated production behavior,further benefitting smart battery production.
基金supported by the National Natural Science Foundation of China(Grant Nos.52109144,52025094 and 52222905).
文摘This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are elucidated geometrically from the perspective of expanding ellipsoids.Based on this geometric interpretation,the QFOSM is further extended to estimate sensitivity indices and assess the significance of various uncertain parameters involved in the slope system.The proposed method has the advantage of computational simplicity,akin to the conventional first-order second-moment method(FOSM),while providing estimation accuracy close to that of the first-order reliability method(FORM).Its performance is demonstrated with a numerical example and three slope examples.The results show that the proposed method can efficiently estimate the slope reliability and simultaneously evaluate the sensitivity of the uncertain parameters.The proposed method does not involve complex optimization or iteration required by the FORM.It can provide a valuable complement to the existing approximate reliability analysis methods,offering rapid sensitivity evaluation and slope reliability analysis.
文摘This paper presents an engineering system approach using a 2D model of conservation of mass to study the dynamics of ozone and concerned chemical species in the stratosphere.By considering all fourteen photolysis,ozone-generating,and-depleting chemical reactions,the model calculated the transient,spatial changes of ozone under different physical-chemical-radiative conditions.Validation against the measured data demonstrated good accuracy,close match of our model with the observed ozone concentrations at both 20°S and 90°N locations.The deviation in the average concentration was less than 1% and in ozone profiles less than 17%.The impacts of various chlorine-(Cl),nitrogen oxides-(NO_(x)),and bromine-(Br)depleting cycles on ozone concentrations and distribution were investigated.The chlorine catalytic depleting cycle was found to exhibit the most significant impact on ozone dynamics,confirming the key role of chlorine in the problem of ozone depletion.Sensitivity analysis was conducted with levels of 25%,50%,100%,200%,and 400% of the baseline value.The combined cycles(Cl+NO_(x)+Br)showed the most significant influence on ozone behavior.The total ozone abundance above the South Pole could decrease by a small 3%,from 281 DU(Dubson Units)to 273 DU for the 25% level,or by a huge thinning of 60%to 114 DU for the 400% concentration level.When the level of chlorine gases increased beyond 200%,it would cause ozone depletion to a level of ozone hole(below 220 DU).The 2D Ozone Model presented in this paper demonstrates robustness,convenience,efficiency,and executability for analyzing complex ozone phenomena in the stratosphere.
基金National Science Foundation for Excellent Young Scholars of China under Grant No.51722801National Natural Science Foundation of China under Grant Nos.51808006 and 52078016。
文摘Case history investigations have shown that pile foundations are more critically damaged in liquefiable soils than non-liquefiable soils.This study examines the differences in seismic response of pile foundations in liquefiable and non-liquefiable soils and their sensitivity to numerical model parameters.A two-dimensional finite element(FE)model is developed to simulate the experiment of a single pile foundation centrifuge in liquefiable soil subjected to earthquake motions and is validated against real-world test results.The differences in soil-pile seismic response of liquefiable and non-liquefiable soils are explored.Specifically,the first-order second-moment method(FOSM)is used for sensitivity analysis of the seismic response.The results show significant differences in seismic response for a soil-pile system between liquefiable and non-liquefiable soil.The seismic responses are found to be significantly larger in liquefiable soil than in non-liquefiable soil.Moreover,the pile bending moment was mainly affected by the kinematic effect in liquefiable soil,while the inertial effect was more significant in non-liquefiable soil.The controlling parameters of seismic response were PGA,soil density,and friction angle in liquefiable soil,while the pile bending moment was mainly controlled by PGA,the friction angle of soil,and shear modulus of loose sand in non-liquefiable soil.
基金supported by the National Natural Science Foundation of China(Grant Nos.42276224,and 42206230)the Jilin Scientific and Technological Development Program(Grant No.20190303083SF)+2 种基金the International Cooperation Key Laboratory of Underground Energy Development and Geological Restoration(Grant No.YDZJ202102CXJD014)the Interdisciplinary Integration and Innovation Project of JLU(Grant No.JLUXKJC2021ZZ18)the Graduate Innovation Fund of Jilin University(Grant No.2023CX100)。
文摘Natural gas hydrate(NGH)is an important future resource for the 21st century and a strategic resource with potential for commercial development in the third energy transition.It is of great significance to accurately predict the productivity of hydrate-bearing sediments(HBS).The multi-phase seepage parameters of HBS include permeability,porosity,which is closely related to permeability,and hydrate saturation,which has a direct impact on hydrate content.Existing research has shown that these multi-phase seepage parameters have a great impact on HBS productivity.Permeability directly affects the transmission of pressure-drop and discharge of methane gas,porosity and initial hydrate saturation affect the amount of hydrate decomposition and transmission process of pressure-drop,and also indirectly affect temperature variation of the reservoir.Considering the spatial heterogeneity of multi-phase seepage parameters,a depressurization production model with layered heterogeneity is established based on the clayey silt hydrate reservoir at W11 station in the Shenhu Sea area of the South China Sea.Tough+Hydrate software was used to calculate the production model;the process of gas production and seepage parameter evolution under different multi-phase seepage conditions were obtained.A sensitivity analysis of the parameters affecting the reservoir productivity was conducted so that:(a)a HBS model with layered heterogeneity can better describe the transmission process of pressure and thermal compensation mechanism of hydrate reservoir;(b)considering the multi-phase seepage parameter heterogeneity,the influence degrees of the parameters on HBS productivity were permeability,porosity and initial hydrate saturation,in order from large to small,and the influence of permeability was significantly greater than that of other parameters;(c)the production potential of the clayey silt reservoir should not only be determined by hydrate content or seepage capacity,but also by the comprehensive effect of the two;and(d)time scales need to be considered when studying the effects of changes in multi-phase seepage parameters on HBS productivity.
文摘Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.
基金supported by the Scientific Innovation Group for Youths of Sichuan Province under Grant No.2019JDTD0017。
文摘Liquefaction is one of the most destructive phenomena caused by earthquakes,which has been studied in the issues of potential,triggering and hazard analysis.The strain energy approach is a common method to investigate liquefaction potential.In this study,two Artificial Neural Network(ANN)models were developed to estimate the liquefaction resistance of sandy soil based on the capacity strain energy concept(W)by using laboratory test data.A large database was collected from the literature.One group of the dataset was utilized for validating the process in order to prevent overtraining the presented model.To investigate the complex influence of fine content(FC)on liquefaction resistance,according to previous studies,the second database was arranged by samples with FC of less than 28%and was used to train the second ANN model.Then,two presented ANN models in this study,in addition to four extra available models,were applied to an additional 20 new samples for comparing their results to show the capability and accuracy of the presented models herein.Furthermore,a parametric sensitivity analysis was performed through Monte Carlo Simulation(MCS)to evaluate the effects of parameters and their uncertainties on the liquefaction resistance of soils.According to the results,the developed models provide a higher accuracy prediction performance than the previously publishedmodels.The sensitivity analysis illustrated that the uncertainties of grading parameters significantly affect the liquefaction resistance of soils.
基金supported by Science and Technology Major Project of Fujian Province(Grant No.2020HZ03018).
文摘For laser cladding a large temperature gradient easily weakened the surface quality by generating cracks and irregular coating surfaces,which in turn affected the bearing capacity and corrosion resistance of coatings in the rapid heating and cooling process.The response surface methodology(RSM)was used to predict coating cracks by changing the powder ratio,energy density,and preheating temperature,which obtained the relevant mathematical model.After that,the sensitivity of the crack length to process parameters was analyzed based on the sensitivity analysis method.The effect of Ni60/WC composite powder process parameters on the surface quality was revealed in laser cladding.The crack length first decreased and then increased,and the Smooth decreased with the increased powder ratio.The crack length and Smooth increased with the increased energy density.The crack length decreased and Smooth increased with the increased preheating temperature.Sensitivity analysis showed that the crack length and Smooth were the most sensitive to the powder ratio.Therefore,the process parameters were reasonably selected to control the surface quality.The mathematical model and sensitivity analysis method in the work could improve the surface quality,which provided a theoretical basis for the prediction and control of laser cladding cracks.
基金supported by the National Nature Science Foundation of China(41975132)the Guangdong Major Project of Basic and Applied Basic Research(Grant No.2020B0301030004).
文摘Simulations and predictions using numerical models show considerable uncertainties,and parameter uncertainty is one of the most important sources.It is impractical to improve the simulation and prediction abilities by reducing the uncertainties of all parameters.Therefore,identifying the sensitive parameters or parameter combinations is crucial.This study proposes a novel approach:conditional nonlinear optimal perturbations sensitivity analysis(CNOPSA)method.The CNOPSA method fully considers the nonlinear synergistic effects of parameters in the whole parameter space and quantitatively estimates the maximum effects of parameter uncertainties,prone to extreme events.Results of the analytical g-function test indicate that the CNOPSA method can effectively identify the sensitivity of variables.Numerical results of the theoretical five-variable grassland ecosystem model show that the maximum influence of the simulated wilted biomass caused by parameter uncertainty can be estimated and computed by employing the CNOPSA method.The identified sensitive parameters can easily change the simulation or prediction of the wilted biomass,which affects the transformation of the grassland state in the grassland ecosystem.The variance-based approach may underestimate the parameter sensitivity because it only considers the influence of limited parameter samples from a statistical view.This study verifies that the CNOPSA method is effective and feasible for exploring the important and sensitive physical parameters or parameter combinations in numerical models.
基金the financial support from the China Postdoctora Science Foundation (project No. 2021M690406)the financial supports from the National Natural Science Foundation of China (project Nos. 51708521, 51778183)
文摘The structural response of a single-layer reticulated dome to external explosions is shaped by many variables,and the associated uncertainties imply non-deterministic results.Existing deterministic methods for predicting the consequences of specific explosions do not account for these uncertainties.Therefore,the impact of the uncertainties associated with these input variables on the structures’response needs to be studied and quantified.In this study,a parametric uncertainty analysis was conducted first.Then,local and global sensitivity analyses were carried out to identify the drivers of the structural dynamic response.A probabilistic structural response model was established based on sensitive variables and a reasonable sample size.Furthermore,some deterministic empirical methods for explosion-resistance design,including the plane blast load model of CONWEP,the curved blast load model under the 50%assurance level,and the 20%mass-increased method,were used for evaluating their reliability.The results of the analyses revealed that the structural response of a single-layer reticulated dome to an external blast loading is lognormally distributed.Evidently,the MB0.5 method based on the curved reflector load model yielded results with a relatively stable assurance rate and reliability,but CONWEP did not;thus,the 1.2MB0.5 method can be used for making high-confidence simple predictions.In addition,the results indicated that the structural response is very sensitive to the explosion parameters.Based on these results,it is suggested that for explosion proofing,setting up a defensive barrier is more effective than structural strengthening.
基金supported by the National Natural Science Foundation of China(Grant No.52079120).
文摘The anti-sliding stability of a gravity dam along its foundation surface is a key problem in the design of gravity dams.In this study,a sensitivity analysis framework was proposed for investigating the factors affecting gravity dam anti-sliding stability along the foundation surface.According to the design specifications,the loads and factors affecting the stability of a gravity dam were comprehensively selected.Afterwards,the sensitivity of the factors was preliminarily analyzed using the Sobol method with Latin hypercube sampling.Then,the results of the sensitivity analysis were verified with those obtained using the Garson method.Finally,the effects of different sampling methods,probability distribution types of factor samples,and ranges of factor values on the analysis results were evaluated.A case study of a typical gravity dam in Yunnan Province of China showed that the dominant factors affecting the gravity dam anti-sliding stability were the anti-shear cohesion,upstream and downstream water levels,anti-shear friction coefficient,uplift pressure reduction coefficient,concrete density,and silt height.Choice of sampling methods showed no significant effect,but the probability distribution type and the range of factor values greatly affected the analysis results.Therefore,these two elements should be sufficiently considered to improve the reliability of the dam anti-sliding stability analysis.
基金supported by National Key R&D Program of China(No.2020YFC2006602)National Natural Science Foundation of China(Nos.62072324,61876217,61876121,61772357)+1 种基金University Natural Science Foundation of Jiangsu Province(No.21KJA520005)Primary Research and Development Plan of Jiangsu Province(No.BE2020026).
文摘Aiming at optimizing the energy consumption of HVAC,an energy conservation optimization method was proposed for HVAC systems based on the sensitivity analysis(SA),named the sensitivity analysis combination method(SAC).Based on the SA,neural network and the related settings about energy conservation of HVAC systems,such as cooling water temperature,chilled water temperature and supply air temperature,were optimized.Moreover,based on the data of the existing HVAC system,various optimal control methods ofHVAC systems were tested and evaluated by a simulated HVAC system in TRNSYS.The results show that the proposed SA combination method can reduce significant computational load while maintaining an equivalent energy performance compared with traditional optimal control methods.
基金supported by the CAS"Light of West China"Program (No.[2020]82)Key technology projects of Inner Mongolia Autonomous Region (Grant No.2020GG0306)+1 种基金Science and Technology Plan Projects of Alxa League (Grant No.AMY2020-18)Natural Science Foundation of Gansu Province (No.21JR7RA038).
文摘A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed by stoma in leaves. In the photosynthesis module of these LSMs, variations of parameters arising from diversity in plant functional types(PFTs) and climate remain unclear. Identifying sensitive parameters among all photosynthetic parameters before parameter estimation can not only reduce operation cost, but also improve the usability of photosynthesis models worldwide. Here, we analyzed 13 parameters of a biochemically-based photosynthesis model(FvCB), implemented in many LSMs, using two sensitivity analysis(SA) methods(i.e., the Sobol’ method and the Morris method) for setting up the parameter ensemble. Three different model performance metrics, i.e.,Root Mean Squared Error(RMSE), Nash Sutcliffe efficiency(NSE), and Standard Deviation(STDEV) were introduced for model assessment and sensitive parameters identification. The results showed that among all photosynthetic parameters only a small portion of parameters were sensitive, and the sensitive parameters were different across plant functional types: maximum rate of Rubisco activity(Vcmax25), maximum electron transport rate(Jmax25), triose phosphate use rate(TPU) and dark respiration in light(Rd) were sensitive in broad leafevergreen trees(BET), broad leaf-deciduous trees(BDT) and needle leaf-evergreen trees(NET), while only Vcmax25and TPU are sensitive in short vegetation(SV), dwarf trees and shrubs(DTS), and agriculture and grassland(AG). The two sensitivity analysis methods suggested a strong SA coherence;in contrast, different model performance metrics led to different SA results. This misfit suggests that more accurate values of sensitive parameters, specifically, species specific and seasonal variable parameters, are required to improve the performance of the FvCB model.
基金This work is supported in part by the Postdoctoral Science Foundation of China under Grant No.2020M683736in part by the Teaching reform project of higher education in Heilongjiang Province under Grant No.SJGY20210456in part by the Natural Science Foundation of Heilongjiang Province of China under Grant No.LH2021F038.
文摘The prediction of processor performance has important referencesignificance for future processors. Both the accuracy and rationality of theprediction results are required. The hierarchical belief rule base (HBRB)can initially provide a solution to low prediction accuracy. However, theinterpretability of the model and the traceability of the results still warrantfurther investigation. Therefore, a processor performance prediction methodbased on interpretable hierarchical belief rule base (HBRB-I) and globalsensitivity analysis (GSA) is proposed. The method can yield more reliableprediction results. Evidence reasoning (ER) is firstly used to evaluate thehistorical data of the processor, followed by a performance prediction modelwith interpretability constraints that is constructed based on HBRB-I. Then,the whale optimization algorithm (WOA) is used to optimize the parameters.Furthermore, to test the interpretability of the performance predictionprocess, GSA is used to analyze the relationship between the input and thepredicted output indicators. Finally, based on the UCI database processordataset, the effectiveness and superiority of the method are verified. Accordingto our experiments, our prediction method generates more reliable andaccurate estimations than traditional models.
基金The work was financially supported by Basic Public Welfare research Plan of Zhejiang Province(LGG19F030006)Key Laboratory of Intelligent Manufacturing Quality Big Data Tracing and Analysis of Zhejiang Province,China Jiliang University(Grant No.ZNZZSZ–CJLU2022–04)the Key Research and Development Program of Ningbo(2022Z165).
文摘The absorber is the key unit in the post-combustion monoethanolamine(MEA)-based carbon dioxide(CO_(2))capture process.A rate-based dynamic model for the absorber is developed and validated using steady-state experimental data reported in open literature.Sensitivity analysis is performed with respect to important model parameters associated with the reaction,mass transport and phy-sical property relationships.Then,a singular value decomposition(SVD)-based subspace parameter estimation method is proposed to improve the model accu-racy.Finally,dynamic simulations are carried out to investigate the effects of the feed rate of lean MEA solution and the flue inlet conditions.Simulation results indicate that the established dynamic model can reasonably reflect the physical behavior of the absorber.Some new insights are gained from the simulation results.
文摘Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.
基金supported by the China Postdoctoral Science Foundation(2021M702304)Natural Science Foundation of Shandong Province(ZR2021QE260).
文摘Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable stress sensitivity characterization models is still limited.In this study,three commonly used stress sensitivity models for shale oil reservoirs were considered,and experiments on representative core samples were conducted.By fitting and comparing the data,the“exponential model”was identified as a characterization model that accurately represents stress sensitivity in continental shale oil reservoirs.To validate the accuracy of the model,a two-phase seepage mathematical model for shale oil reservoirs coupled with the exponential model was introduced.The model was discretely solved using the finite volume method,and its accuracy was verified through the commercial simulator CMG.The study evaluated the productivity of a typical horizontal well under different engineering,geological,and fracture conditions.The results indicate that considering stress sensitivity leads to a 13.57%reduction in production for the same matrix permeability.Additionally,as the fracture half-length and the number of fractures increase,and the bottomhole flowing pressure decreases,the reservoir stress sensitivity becomes higher.