Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the sel...Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the selection and application of higher sensitivity seismic attribute parameters to improve the ability to discriminate fluid and lithology. To better integrate with seismic interpretation, we construct a template to highlight rock physics parameters in sensitivity space, providing guidance for the quantitative seismic interpretation of hydrocarbon-bearing reservoirs.展开更多
The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which compris...The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.展开更多
Population viability analysis(PVA) is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dy...Population viability analysis(PVA) is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.展开更多
Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about di...Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.展开更多
In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The for...In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate（ZooAE_N）,zooplankton basal metabolism rate（ZooBM）,maximum specific growth rate of zooplankton（μ_（20）） and maximum chlorophyll to carbon ratio（Chl2C_m） are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_（20＇）,half-saturation for phytoplankton NO_3 uptake（K_（NO_3）） and remineralization rate of small detritusN（SDeRRN） are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^（-1）.The K_（NO_3） and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5（mmol/m^3）^（-1）,and interval of SEdRRN is from 0.01 and 0.1 d^（-1）.The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.展开更多
In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through...In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.展开更多
Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to ...Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
Acoustic and electrical methods are commonly used to evaluate hydrate saturation based on P-wave velocity(Vp)and resistivity,respectively.We evaluate hydrate saturation using petrophysical parameters directly related ...Acoustic and electrical methods are commonly used to evaluate hydrate saturation based on P-wave velocity(Vp)and resistivity,respectively.We evaluate hydrate saturation using petrophysical parameters directly related to the presence of hydrates.Five petrophysical parameters sensitive to hydrate saturation were first analyzed using the equivalent medium rock physical model,logging intersection plots,and petrophysical parameter inversion.The simulated annealing global optimization method was then used to estimate the hydrate saturation profile in the Shenhu Area,China.The petrophysical parameters Vp,λρ,andλμ,which are associated with the rock elastic and shear moduli,are highly sensitive to hydrate saturation for an estimated saturation range of 0.1-0.44.This range is consistent with that obtained from the original well diameter curves.However,the parameters Vs andμρ,which are only related to the rock shear modulus,yield high hydrate saturation estimates of 0.22-0.43 and exhibit some deviations from the real-time data.Owing to its sensitivity,the Poisson’s ratio is least desired for hydrate evaluation among the studied parameters.The sensitivity of hydrate saturation depends on the petrophysical model used for studying hydrate physical properties and storage analysis.展开更多
Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dyn...Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.展开更多
It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continui...It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continuing variation on the result and the effects of exceeding probability and upper bound magnitude of potential source on the parameter sensitivity. Some useful cognitions were got. We have overall calculated the interactions among the parameters including spatial distribution function of earthquake and suggested the concept of relative interaction. The relative interaction is not great than 10% with annual exceeding probability 2×10 -3 . The relative interaction between upper bound magnitude and annual mean rate and the one between upper bound magnitude and space distribution function are 17% and 22% respectively for potantial sources higher upper bound magnitudes in annual exceeding probability 10 -4 . It is convenient to judge whether considering the interactions in actual application. The importance and necessity of calculating the result of the probabilistic distribution are indicated. The effect of the result distribution by different parameter level weights is analyzed. The whole picture of multi parameter sensitivity is revealed. It makes necessary preparation for the reasonable uncertainty correction.展开更多
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.展开更多
In this paper, we are interested to find the most sensitive parameter, local and global stability of ovarian tumor growth model. For sensitivity analysis, we use Latin Hypercube Sampling (LHS) method to generate sampl...In this paper, we are interested to find the most sensitive parameter, local and global stability of ovarian tumor growth model. For sensitivity analysis, we use Latin Hypercube Sampling (LHS) method to generate sample points and Partial Rank Correlation Coefficient (PRCC) method, uses those sample points to find out which parameters are important for the model. Based on our findings, we suggest some treatment strategies. We investigate the sensitivity of the parameters for tumor volume, <em>y</em>, cell nutrient density, <em>Q</em> and maximum tumor size, <em>ymax</em>. We also use Scatter Plot method using LHS samples to show the consistency of the results obtained by using PRCC. Moreover, we discuss the qualitative analysis of ovarian tumor growth model investigating the local and global stability.展开更多
This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenv...This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.展开更多
Under severe sea conditions, wave slamming on ships and marine engineering structures may lead to structural damage and casualties. Moreover, the strong nonlinearity inherent in the wave slamming process significantly...Under severe sea conditions, wave slamming on ships and marine engineering structures may lead to structural damage and casualties. Moreover, the strong nonlinearity inherent in the wave slamming process significantly limits the accuracy of numerical analyses and finite element simulations. Therefore, this paper takes a new type of floating wind turbine as an example and performs a physical model test on the wave slamming characteristics of this floating wind turbine.Based on a 1:50 model of the Pivot Buoy floating wind turbine, an experimental study is performed under the combined effects of wind-wave loads on the peak pressure, duration, and pressure distribution of slamming. First, two sets of mooring systems, the combined scheme and the full mooring chain scheme, are designed to conduct a series of experimental studies of model slamming under different wind and wave incidence angles, wave heights, and wave periods. By doing so, the slamming characteristics of the wind turbine can be obtained. Moreover, to solve the problem of the large pitch motion response of the prototype wind turbine, a set of vertically oscillating structures is designed,and the slamming pressure characteristics of the optimized model are also investigated through model tests.展开更多
Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （...Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （i） determine influential and non-influential parameters with respect to above ground biomass （AGB）, canopy cover （CC）, and grain yield of winter wheat in the Beijing area based on the AquaCrop model under different water treatments （rainfall, normal irrigation, and over-irrigation）; and （ii） generate an AquaCrop model that can be used in the Beijing area by setting non-influential parameters to fixed values and adjusting influential parameters according to the SA results. In this study, field experiments were conducted during the 2012-2013,2013-2014, and 2014-2015 winter wheat growing seasons at the National Precision Agriculture Demonstration Research Base in Beijing, China. The extended Fourier amplitude sensitivity test （EFAST） method was used to perform SA of the AquaCrop model using 42 crop parameters, in order to verify the SA results, data from the 2013-2014 growing season were used to calibrate the AquaCrop model, and data from 2012-2013 and 2014-2015 growing seasons were val- idated. For AGB and yield of winter wheat, the total order sensitivity analysis had more sensitive parameters than the first order sensitivity analysis. For the AGB time-series, parameter sensitivity was changed under different water treatments; in comparison with the non-stressful conditions （normal irrigation and over-irrigation）, there were more sensitive parameters under water stress （rainfall）, while root development parameters were more sensitive. For CC with time-series and yield, there were more sensitive parameters under water stress than under no water stress. Two parameters sets were selected to calibrate the AquaCrop model, one group of parameters were under water stress, and the others were under no water stress, there were two more sensitive parameters （growing degree-days （GDD） from sowing to the maximum rooting depth （root） and the maximum effective rooting depth （rtx）） under water stress than under no water stress. The results showed that there was higher accuracy under water stress than under no water stress. This study provides guidelines for AquaCrop model calibration and application in Beijing, China, as well providing guidance to simplify the AquaCrop model and improve its precision, especially when many parameters are used.展开更多
During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.Th...During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.The radical basis function（RBF）neural network model of a five-link biped robot is established,and two certain disturbances and a randomly uncertain disturbance are then mixed with the optimal torques in the network model to study the performance of the biped robot by several evaluation indices and a specific Poincar′e map.In contrast with the simulations,the response varies as desired under optimal inputting while the output is fluctuating in the situation of disturbance driving.Simulation results from noise inputting also show that the dynamics of the robot is less sensitive to the disturbance of knee joint input of the swing leg than those of the other three joints,the response errors of the biped will be increasing with higher disturbance levels,and especially there are larger output fluctuations in the knee and hip joints of the swing leg.展开更多
Traditionally,semi-rigid analysis of extended end plate connections is assumed to be deterministic parameters.Most of the moment-rotation models of the joint are determined by curve-fitting and statistically regressed...Traditionally,semi-rigid analysis of extended end plate connections is assumed to be deterministic parameters.Most of the moment-rotation models of the joint are determined by curve-fitting and statistically regressed with the test data.Parameters of these models have little physical meaning,and the effect of correlations in the parameters is neglected.This paper deals with a semi-rigid joint of extended end plate connections based on correlation and sensitivity analysis.A numerical study was carried out and the finite element analysis was validated compared with the experimental results.The relative moment-rotation model was concluded based on correlations in parameters and probabilistic sentivity analysis.The relationship between height and width of steel beam was partial correlation.The relationship between thickness of steel beam flange and thickness of steel beam web was also partial correlation.Height and width of steel beam,and distance of upper bolts to axle wire of steel beam section were important sensitive parameters for the extended end plate connections.展开更多
To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale co...To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.展开更多
The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slop...The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slope.The least squares support vector machine(LSSVM)evolved by a bacterial foraging optimization algorithm(BFOA)is used to establish a response surface model to express the mapping relationship between the intermittent joint parameters and the slope safety factor.The training samples are obtained from the numerical calculation based on the joint finite element method during this process.Considering the randomness of the intermittent joint parameters in the actual project,each parameter is evaluated at different locations on the site,and its distribution characteristics are counted.According to these statistical results,a large number of parameter combinations are obtained through Monte Carlo sampling.The trained machine learning mapping model is used to obtain the slope safety factor corresponding to each group,and these results are then used to obtain the slope reliability.When the research results were applied to slope disaster treatment along the Yalu River in China’s Jilin Province,it was found that the joint length and joint inclination angle both play key roles in rock slope stability,which should receive more attention in the slope treatment.In summary,this study establishes a method for evaluating the reliability of intermittent jointed rock slope based on an evolutionary SVM model,and its feasibility is verified by engineering application.展开更多
文摘Proper analysis of petrophysical parameter sensitivity from well log data can greatly improve the ability to discriminate hydrocarbon-bearing rocks. In this paper we discuss a petrophysical analysis method for the selection and application of higher sensitivity seismic attribute parameters to improve the ability to discriminate fluid and lithology. To better integrate with seismic interpretation, we construct a template to highlight rock physics parameters in sensitivity space, providing guidance for the quantitative seismic interpretation of hydrocarbon-bearing reservoirs.
基金Project(2015B37714)supported by the Fundamental Research Funds for the Central Universities of ChinaProject(51605005)supported by the National Natural Science Foundation of China+1 种基金Project(ZK16-03-03)supported by the Open Foundation of Jiangsu Wind Technology Center,ChinaProject([2013]56)supported by the First Group of 2011 Plan of Jiangsu Province,China
文摘The parameter sensitivities affecting the flutter speed of the NREL （National Renewable Energy Laboratory） 5-MW baseline HAWT （horizontal axis wind turbine） blades are analyzed. An aeroelastic model, which comprises an aerodynamic part to calculate the aerodynamic loads and a structural part to determine the structural dynamic responses, is established to describe the classical flutter of the blades. For the aerodynamic part, Theodorsen unsteady aerodynamics model is used. For the structural part, Lagrange’s equation is employed. The flutter speed is determined by introducing “V–g” method to the aeroelastic model, which converts the issue of classical flutter speed determination into an eigenvalue problem. Furthermore, the time domain aeroelastic response of the wind turbine blade section is obtained with employing Runge-Kutta method. The results show that four cases （i.e., reducing the blade torsional stiffness, moving the center of gravity or the elastic axis towards the trailing edge of the section, and placing the turbine in high air density area） will decrease the flutter speed. Therefore, the judicious selection of the four parameters （the torsional stiffness, the chordwise position of the center of gravity, the elastic axis position and air density） can increase the relative inflow speed at the blade section associated with the onset of flutter.
基金funded by the State Forestry Administration,China
文摘Population viability analysis(PVA) is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.
基金The first author would like to thank the China Scholarship Council(CSC No.201806430001)for the financial support for his study in Germany.The authors are grateful to the anonymous reviewer for the recommendations to improve the paper。
文摘Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.
基金The National Natural Science Foundation of China under contract Nos 41206023,41222038 and 41076011the National Basic Research Project(973 Program)of China under contract No.2011CB403606+2 种基金the China-Korea Joint Ocean Research Center"Cooperation on the Development of Basic Technologies for the Yellow Sea and East China Sea Operational Oceanographic System(YOOS)"the Public Science and Technology Research Funds Projects of Ocean under contrcat No.201205018the"Strategic Priority Research Program"of the Chinese Academy of Sciences,under contract No.XDA01020304
文摘In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas（CCS）.The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate（ZooAE_N）,zooplankton basal metabolism rate（ZooBM）,maximum specific growth rate of zooplankton（μ_（20）） and maximum chlorophyll to carbon ratio（Chl2C_m） are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_（20＇）,half-saturation for phytoplankton NO_3 uptake（K_（NO_3）） and remineralization rate of small detritusN（SDeRRN） are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^（-1）.The K_（NO_3） and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5（mmol/m^3）^（-1）,and interval of SEdRRN is from 0.01 and 0.1 d^（-1）.The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.
基金Supported by the National Natural Science Foundation of China(10272109)
文摘In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.
基金supported by the National Natural Science Foundation of China (11142013, 11172260 and 11072214)the Doctoral Fund of Ministry of Education of China (20110101110016)the Fundamental Research Funds for the Central Universities of China(2011QNA4001)
文摘Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
基金supported by the National Natural Science Foundation of China(Nos.91958206,41876053)the National Key Research and Development Plan(2017YFC0307401,2018YFC1405901)+1 种基金the Fund of Acoustics Science and Technology Laboratory(GK2050260214,GK2050260217,GK2050260218,KY10500180084,KY10500190031,6142108200202)Fundamental Research Funds for the Central Universities(HEUCFJ180503,201964016).
文摘Acoustic and electrical methods are commonly used to evaluate hydrate saturation based on P-wave velocity(Vp)and resistivity,respectively.We evaluate hydrate saturation using petrophysical parameters directly related to the presence of hydrates.Five petrophysical parameters sensitive to hydrate saturation were first analyzed using the equivalent medium rock physical model,logging intersection plots,and petrophysical parameter inversion.The simulated annealing global optimization method was then used to estimate the hydrate saturation profile in the Shenhu Area,China.The petrophysical parameters Vp,λρ,andλμ,which are associated with the rock elastic and shear moduli,are highly sensitive to hydrate saturation for an estimated saturation range of 0.1-0.44.This range is consistent with that obtained from the original well diameter curves.However,the parameters Vs andμρ,which are only related to the rock shear modulus,yield high hydrate saturation estimates of 0.22-0.43 and exhibit some deviations from the real-time data.Owing to its sensitivity,the Poisson’s ratio is least desired for hydrate evaluation among the studied parameters.The sensitivity of hydrate saturation depends on the petrophysical model used for studying hydrate physical properties and storage analysis.
文摘Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.
文摘It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continuing variation on the result and the effects of exceeding probability and upper bound magnitude of potential source on the parameter sensitivity. Some useful cognitions were got. We have overall calculated the interactions among the parameters including spatial distribution function of earthquake and suggested the concept of relative interaction. The relative interaction is not great than 10% with annual exceeding probability 2×10 -3 . The relative interaction between upper bound magnitude and annual mean rate and the one between upper bound magnitude and space distribution function are 17% and 22% respectively for potantial sources higher upper bound magnitudes in annual exceeding probability 10 -4 . It is convenient to judge whether considering the interactions in actual application. The importance and necessity of calculating the result of the probabilistic distribution are indicated. The effect of the result distribution by different parameter level weights is analyzed. The whole picture of multi parameter sensitivity is revealed. It makes necessary preparation for the reasonable uncertainty correction.
基金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.
文摘In this paper, we are interested to find the most sensitive parameter, local and global stability of ovarian tumor growth model. For sensitivity analysis, we use Latin Hypercube Sampling (LHS) method to generate sample points and Partial Rank Correlation Coefficient (PRCC) method, uses those sample points to find out which parameters are important for the model. Based on our findings, we suggest some treatment strategies. We investigate the sensitivity of the parameters for tumor volume, <em>y</em>, cell nutrient density, <em>Q</em> and maximum tumor size, <em>ymax</em>. We also use Scatter Plot method using LHS samples to show the consistency of the results obtained by using PRCC. Moreover, we discuss the qualitative analysis of ovarian tumor growth model investigating the local and global stability.
基金Project supported by the 985-Engineering Innovation of Graduate Students of Jilin Universitythe Science and Technology Development Foundation of Jilin Province(20070541)
文摘This paper presents methods for computing a second-order sensitivity matrix and the Hessian matrix of eigenvalues and eigenvectors of multiple parameter structures. Second-order perturbations of eigenvalues and eigenvectors are transformed into multiple parameter forms,and the second-order perturbation sensitivity matrices of eigenvalues and eigenvectors are developed.With these formulations,the efficient methods based on the second-order Taylor expansion and second-order perturbation are obtained to estimate changes of eigenvalues and eigenvectors when the design parameters are changed. The presented method avoids direct differential operation,and thus reduces difficulty for computing the second-order sensitivity matrices of eigenpairs.A numerical example is given to demonstrate application and accuracy of the proposed method.
基金financially supported by the National Natural Science Foundation of China (Grant No. 52071161)。
文摘Under severe sea conditions, wave slamming on ships and marine engineering structures may lead to structural damage and casualties. Moreover, the strong nonlinearity inherent in the wave slamming process significantly limits the accuracy of numerical analyses and finite element simulations. Therefore, this paper takes a new type of floating wind turbine as an example and performs a physical model test on the wave slamming characteristics of this floating wind turbine.Based on a 1:50 model of the Pivot Buoy floating wind turbine, an experimental study is performed under the combined effects of wind-wave loads on the peak pressure, duration, and pressure distribution of slamming. First, two sets of mooring systems, the combined scheme and the full mooring chain scheme, are designed to conduct a series of experimental studies of model slamming under different wind and wave incidence angles, wave heights, and wave periods. By doing so, the slamming characteristics of the wind turbine can be obtained. Moreover, to solve the problem of the large pitch motion response of the prototype wind turbine, a set of vertically oscillating structures is designed,and the slamming pressure characteristics of the optimized model are also investigated through model tests.
基金supported by the National Natural Science Foundation of China(41571416)the Natural Science Foundation of Beijing,China(4152019)the Beijing Academy of Agricultural and Forestry Sciences Innovation Capacity Construction Specific Projects,China(KJCX20150409)
文摘Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （i） determine influential and non-influential parameters with respect to above ground biomass （AGB）, canopy cover （CC）, and grain yield of winter wheat in the Beijing area based on the AquaCrop model under different water treatments （rainfall, normal irrigation, and over-irrigation）; and （ii） generate an AquaCrop model that can be used in the Beijing area by setting non-influential parameters to fixed values and adjusting influential parameters according to the SA results. In this study, field experiments were conducted during the 2012-2013,2013-2014, and 2014-2015 winter wheat growing seasons at the National Precision Agriculture Demonstration Research Base in Beijing, China. The extended Fourier amplitude sensitivity test （EFAST） method was used to perform SA of the AquaCrop model using 42 crop parameters, in order to verify the SA results, data from the 2013-2014 growing season were used to calibrate the AquaCrop model, and data from 2012-2013 and 2014-2015 growing seasons were val- idated. For AGB and yield of winter wheat, the total order sensitivity analysis had more sensitive parameters than the first order sensitivity analysis. For the AGB time-series, parameter sensitivity was changed under different water treatments; in comparison with the non-stressful conditions （normal irrigation and over-irrigation）, there were more sensitive parameters under water stress （rainfall）, while root development parameters were more sensitive. For CC with time-series and yield, there were more sensitive parameters under water stress than under no water stress. Two parameters sets were selected to calibrate the AquaCrop model, one group of parameters were under water stress, and the others were under no water stress, there were two more sensitive parameters （growing degree-days （GDD） from sowing to the maximum rooting depth （root） and the maximum effective rooting depth （rtx）） under water stress than under no water stress. The results showed that there was higher accuracy under water stress than under no water stress. This study provides guidelines for AquaCrop model calibration and application in Beijing, China, as well providing guidance to simplify the AquaCrop model and improve its precision, especially when many parameters are used.
基金supported by the Science Fund for Creative Research Groups of National Natural Science Foundation of China(51221004)the National Natural Science Foundation of China(11172260,11372270,and 51375434)+2 种基金the Higher School Specialized Research Fund for the Doctoral Program(20110101110016)the Science and technology project of Zhejiang Province(2013C31086)the Fundamental Research Funds forthe Central Universities of China(2013XZZX005)
文摘During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.The radical basis function（RBF）neural network model of a five-link biped robot is established,and two certain disturbances and a randomly uncertain disturbance are then mixed with the optimal torques in the network model to study the performance of the biped robot by several evaluation indices and a specific Poincar′e map.In contrast with the simulations,the response varies as desired under optimal inputting while the output is fluctuating in the situation of disturbance driving.Simulation results from noise inputting also show that the dynamics of the robot is less sensitive to the disturbance of knee joint input of the swing leg than those of the other three joints,the response errors of the biped will be increasing with higher disturbance levels,and especially there are larger output fluctuations in the knee and hip joints of the swing leg.
基金the National Natural Science Foundation of China(Nos.50978104,51008129 and 51178192)the China Postdoctoral Science Foundation (No.20110490895)+1 种基金the "Code for Design of Steel Structures" National Institute of Standards Management Group(No.GB500172010-04)the Fundamental Research Funds for the Central Universities of China (No.2012ZZ0099)
文摘Traditionally,semi-rigid analysis of extended end plate connections is assumed to be deterministic parameters.Most of the moment-rotation models of the joint are determined by curve-fitting and statistically regressed with the test data.Parameters of these models have little physical meaning,and the effect of correlations in the parameters is neglected.This paper deals with a semi-rigid joint of extended end plate connections based on correlation and sensitivity analysis.A numerical study was carried out and the finite element analysis was validated compared with the experimental results.The relative moment-rotation model was concluded based on correlations in parameters and probabilistic sentivity analysis.The relationship between height and width of steel beam was partial correlation.The relationship between thickness of steel beam flange and thickness of steel beam web was also partial correlation.Height and width of steel beam,and distance of upper bolts to axle wire of steel beam section were important sensitive parameters for the extended end plate connections.
基金supported by the National Key R&D Program of China(Grant No.2018YFC1405900)the Major Projects of National Science and Technology(Grant Nos.2016ZX05011-002,2016ZX05027-002-005)+3 种基金the National Natural Science Foundation of China(Grant No.41806073)the Natural Science Foundation of Shandong Province(Grant No.ZR2017BD014)Shandong Provincial Key Laboratory of Depositional Mineralization and Sedimentary Minerals,Shandong University of Science and Technology(Grant No.DMSM2017042)the Fundamental Research Funds for the Central Universities(Grant No.201964016)
文摘To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.
基金The authors sincerely appreciate the support from the National Natural Science Foundation of China[Grant Nos.51678101,52078093]Liaoning Revitalization Talents Program[Grant No.XLYC1905015]the Doctoral innovation Program of Dalian Maritime University[Grant No.BSCXXM016].
文摘The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slope.The least squares support vector machine(LSSVM)evolved by a bacterial foraging optimization algorithm(BFOA)is used to establish a response surface model to express the mapping relationship between the intermittent joint parameters and the slope safety factor.The training samples are obtained from the numerical calculation based on the joint finite element method during this process.Considering the randomness of the intermittent joint parameters in the actual project,each parameter is evaluated at different locations on the site,and its distribution characteristics are counted.According to these statistical results,a large number of parameter combinations are obtained through Monte Carlo sampling.The trained machine learning mapping model is used to obtain the slope safety factor corresponding to each group,and these results are then used to obtain the slope reliability.When the research results were applied to slope disaster treatment along the Yalu River in China’s Jilin Province,it was found that the joint length and joint inclination angle both play key roles in rock slope stability,which should receive more attention in the slope treatment.In summary,this study establishes a method for evaluating the reliability of intermittent jointed rock slope based on an evolutionary SVM model,and its feasibility is verified by engineering application.
