Clastic rock reservoir is the main reservoir type in the oil and gas field.Archie formula or various conductive models developed on the basis of Archie’s formula are usually used to interpret this kind of reservoir,a...Clastic rock reservoir is the main reservoir type in the oil and gas field.Archie formula or various conductive models developed on the basis of Archie’s formula are usually used to interpret this kind of reservoir,and the three-water model is widely used as well.However,there are many parameters in the threewater model,and some of them are difficult to determine.Most of the determination methods are based on the statistics of large amount of experimental data.In this study,the authors determine the value of the parameters of the new three-water model based on the nuclear magnetic data and the genetic optimization algorithm.The relative error between the resistivity calculated based on these parameters and the resistivity measured experimentally at 100%water content is 0.9024.The method studied in this paper can be easily applied without much experimental data.It can provide reference for other regions to determine the parameters of the new three-water model.展开更多
The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The bas...The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.展开更多
Typically,a computer has infectivity as soon as it is infected.It is a reality that no antivirus programming can identify and eliminate all kinds of viruses,suggesting that infections would persevere on the Internet.T...Typically,a computer has infectivity as soon as it is infected.It is a reality that no antivirus programming can identify and eliminate all kinds of viruses,suggesting that infections would persevere on the Internet.To understand the dynamics of the virus propagation in a better way,a computer virus spread model with fuzzy parameters is presented in this work.It is assumed that all infected computers do not have the same contribution to the virus transmission process and each computer has a different degree of infectivity,which depends on the quantity of virus.Considering this,the parametersβandγbeing functions of the computer virus load,are considered fuzzy numbers.Using fuzzy theory helps us understand the spread of computer viruses more realistically as these parameters have fixed values in classical models.The essential features of the model,like reproduction number and equilibrium analysis,are discussed in fuzzy senses.Moreover,with fuzziness,two numerical methods,the forward Euler technique,and a nonstandard finite difference(NSFD)scheme,respectively,are developed and analyzed.In the evidence of the numerical simulations,the proposed NSFD method preserves the main features of the dynamic system.It can be considered a reliable tool to predict such types of solutions.展开更多
Model parameters estimation is a pivotal issue for runoff modeling in ungauged catchments.The nonlinear relationship between model parameters and catchment descriptors is a major obstacle for parameter regionalization...Model parameters estimation is a pivotal issue for runoff modeling in ungauged catchments.The nonlinear relationship between model parameters and catchment descriptors is a major obstacle for parameter regionalization,which is the most widely used approach.Runoff modeling was studied in 38 catchments located in the Yellow–Huai–Hai River Basin(YHHRB).The values of the Nash–Sutcliffe efficiency coefficient(NSE),coefficient of determination(R2),and percent bias(PBIAS)indicated the acceptable performance of the soil and water assessment tool(SWAT)model in the YHHRB.Nine descriptors belonging to the categories of climate,soil,vegetation,and topography were used to express the catchment characteristics related to the hydrological processes.The quantitative relationships between the parameters of the SWAT model and the catchment descriptors were analyzed by six regression-based models,including linear regression(LR)equations,support vector regression(SVR),random forest(RF),k-nearest neighbor(kNN),decision tree(DT),and radial basis function(RBF).Each of the 38 catchments was assumed to be an ungauged catchment in turn.Then,the parameters in each target catchment were estimated by the constructed regression models based on the remaining 37 donor catchments.Furthermore,the similaritybased regionalization scheme was used for comparison with the regression-based approach.The results indicated that the runoff with the highest accuracy was modeled by the SVR-based scheme in ungauged catchments.Compared with the traditional LR-based approach,the accuracy of the runoff modeling in ungauged catchments was improved by the machine learning algorithms because of the outstanding capability to deal with nonlinear relationships.The performances of different approaches were similar in humid regions,while the advantages of the machine learning techniques were more evident in arid regions.When the study area contained nested catchments,the best result was calculated with the similarity-based parameter regionalization scheme because of the high catchment density and short spatial distance.The new findings could improve flood forecasting and water resources planning in regions that lack observed data.展开更多
We discuss formulas and techniques for finding maximum-likelihood estimators of parameters of autoregressive (with particular emphasis on Markov and Yule) models, computing their asymptotic variance-covariance matrix ...We discuss formulas and techniques for finding maximum-likelihood estimators of parameters of autoregressive (with particular emphasis on Markov and Yule) models, computing their asymptotic variance-covariance matrix and displaying the resulting confidence regions;Monte Carlo simulation is then used to establish the accuracy of the corresponding level of confidence. The results indicate that a direct application of the Central Limit Theorem yields errors too large to be acceptable;instead, we recommend using a technique based directly on the natural logarithm of the likelihood function, verifying its substantially higher accuracy. Our study is then extended to the case of estimating only a subset of a model’s parameters, when the remaining ones (called nuisance) are of no interest to us.展开更多
The effects of the milling parameters involving shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM modeling in the planetary ball mill were investigated. BPDR was varied from 1 to 10....The effects of the milling parameters involving shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM modeling in the planetary ball mill were investigated. BPDR was varied from 1 to 10. The results revealed that the size and shape of the powder particles do not give a significant change in simulation results when BPDR attains maximum value of 10. The increasing of BPDR leads to the increase of simulation time and size. Hence, the effect of change of the powder particle shape on the calculated data size is not significant. The results also revealed that the increasing rotation speed increases impact energy between powder particles.展开更多
To estimate the parameters of the mixed additive and multiplicative(MAM)random error model using the weighted least squares iterative algorithm that requires derivation of the complex weight array,we introduce a deriv...To estimate the parameters of the mixed additive and multiplicative(MAM)random error model using the weighted least squares iterative algorithm that requires derivation of the complex weight array,we introduce a derivative-free cat swarm optimization for parameter estimation.We embed the Powell method,which uses conjugate direction acceleration and does not need to derive the objective function,into the original cat swarm optimization to accelerate its convergence speed and search accuracy.We use the ordinary least squares,weighted least squares,original cat swarm optimization,particle swarm algorithm and improved cat swarm optimization to estimate the parameters of the straight-line fitting MAM model with lower nonlinearity and the DEM MAM model with higher nonlinearity,respectively.The experimental results show that the improved cat swarm optimization has faster convergence speed,higher search accuracy,and better stability than the original cat swarm optimization and the particle swarm algorithm.At the same time,the improved cat swarm optimization can obtain results consistent with the weighted least squares method based on the objective function only while avoiding multiple complex weight array derivations.The method in this paper provides a new idea for theoretical research on parameter estimation of MAM error models.展开更多
The small-scale drilling technique can be a fast and reliable method to estimate rock strength parameters. It needs to link the operational drilling parameters and strength properties of rock. The parameters such as b...The small-scale drilling technique can be a fast and reliable method to estimate rock strength parameters. It needs to link the operational drilling parameters and strength properties of rock. The parameters such as bit geometry, bit movement, contact frictions and crushed zone affect the estimated parameters.An analytical model considering operational drilling data and effective parameters can be used for these purposes. In this research, an analytical model was developed based on limit equilibrium of forces in a Tshaped drag bit considering the effective parameters such as bit geometry, crushed zone and contact frictions in drilling process. Based on the model, a method was used to estimate rock strength parameters such as cohesion, internal friction angle and uniaxial compressive strength of different rock types from operational drilling data. Some drilling tests were conducted by a portable and powerful drilling machine which was developed for this work. The obtained results for strength properties of different rock types from the drilling experiments based on the proposed model are in good agreement with the results of standard tests. Experimental results show that the contact friction between the cutting face and rock is close to that between bit end wearing face and rock due to the same bit material. In this case,the strength parameters, especially internal friction angle and cohesion, are estimated only by using a blunt bit drilling data and the bit bluntness does not affect the estimated results.展开更多
Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes ...Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.展开更多
By using the widely used JOHNSON-COOK model and the gradient-dependent plasticity to consider microstructural effect beyond the occurrence of shear strain localization,the distributions of local plastic shear strain a...By using the widely used JOHNSON-COOK model and the gradient-dependent plasticity to consider microstructural effect beyond the occurrence of shear strain localization,the distributions of local plastic shear strain and deformation in adiabatic shear band(ASB)were analyzed.The peak local plastic shear strain is proportional to the average plastic shear strain,while it is inversely proportional to the critical plastic shear strain corresponding to the peak flow shear stress.The relative plastic shear deformation between the top and base of ASB depends on the thickness of ASB and the average plastic shear strain.A parametric study was carried out to study the influence of constitutive parameters on shear strain localization.Higher values of static shear strength and work to heat conversion factor lead to lower critical plastic shear strain so that the shear localization is more apparent at the same average plastic shear strain.Higher values of strain-hardening exponent,strain rate sensitive coefficient,melting point,thermal capacity and mass density result in higher critical plastic shear strain,leading to less apparent shear localization at the same average plastic shear strain.The strain rate sensitive coefficient has a minor influence on the critical plastic shear strain,the distributions of local plastic shear strain and deformation in ASB.The effect of strain-hardening modulus on the critical plastic shear strain is not monotonous.When the maximum critical plastic shear strain is reached,the least apparent shear localization occurs.展开更多
The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows ...The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.展开更多
Transmission line(TL)Parameter Identification(PI)method plays an essential role in the transmission system.The existing PI methods usually have two limitations:(1)These methods only model for single TL,and can not con...Transmission line(TL)Parameter Identification(PI)method plays an essential role in the transmission system.The existing PI methods usually have two limitations:(1)These methods only model for single TL,and can not consider the topology connection of multiple branches for simultaneous identification.(2)Transient bad data is ignored by methods,and the random selection of terminal section data may cause the distortion of PI and have serious consequences.Therefore,a multi-task PI model considering multiple TLs’spatial constraints and massive electrical section data is proposed in this paper.The Graph Attention Network module is used to draw a single TL into a node and calculate its influence coefficient in the transmission network.Multi-Task strategy of Hard Parameter Sharing is used to identify the conductance ofmultiple branches simultaneously.Experiments show that themethod has good accuracy and robustness.Due to the consideration of spatial constraints,the method can also obtain more accurate conductance values under different training and testing conditions.展开更多
Herein,iron oxide/hydroxides deposits(gossans)were utilized,for the first time,in the fabrication of magnetite nanoparticles(MNPs)to load modified coal(MC).The as-synthesized MNPs@MC composite was characterized via di...Herein,iron oxide/hydroxides deposits(gossans)were utilized,for the first time,in the fabrication of magnetite nanoparticles(MNPs)to load modified coal(MC).The as-synthesized MNPs@MC composite was characterized via different techniques and utilized for the Cr(Ⅵ)remediation.Experimental studies supported by theoretical treatment were applied to offer a new overview of the Cr(Ⅵ)adsorption geometry and mechanism at 25-45℃.Experimental results suggested that the Cr(Ⅵ)uptake was mainly governed by adsorption-reduction coupled mechanism.The Langmuir model fitted well the Cr(Ⅵ)adsorption data with maximum adsorption capacities extended from 115.24 to 129.63 mg·g^(-1).Theoretical calculations indicated that Cr(Ⅵ)ions were adsorbed on the MNPs@MC following the theory of the advanced monolayer statistical model.The number of ions removed per site ranged from 1.88 to1.23 suggesting the involvement of vertical geometry and multi-ionic mechanism at all temperatures.The increment of the active sites density and the adsorption capacity at saturation with improving temperature reflected an endothermic process.Energetically,the Cr(Ⅵ)adsorption was controlled by physical forces as the adsorption energies were less than 40 kJ·mol^(-1).The calculated free enthalpy,entropy.and internal energy explained the spontaneous nature and the viability of Cr(Ⅵ)adsorption on the MNPs@MC adsorbent.These results offer a new approach in utilizing the iron-rich deposits as gossans in the preparation of magnetic and low-cost adsorbents for wastewater remediation.展开更多
Minimizing parameter uncertainty is crucial in the application of hydrologic models.Isotopic information in various hydrologic components of the water cycle can expand our knowledge of the dynamics of water flow in th...Minimizing parameter uncertainty is crucial in the application of hydrologic models.Isotopic information in various hydrologic components of the water cycle can expand our knowledge of the dynamics of water flow in the system,provide additional information for parameter estimation,and improve parameter identifiability.This study combined the Philip infiltration model with an isotopic mixing model using an isotopic mass balance approach for estimating parameters in the Philip infiltration model.Two approaches to parameter estimation were compared:(a) using isotopic information to determine the soil water transmission and then hydrologic information to estimate the soil sorptivity,and(b) using hydrologic information to determine the soil water transmission and the soil sorptivity.Results of parameter estimation were verified through a rainfall infiltration experiment in a laboratory under rainfall with constant isotopic compositions and uniform initial soil water content conditions.Experimental results showed that approach(a),using isotopic and hydrologic information,estimated the soil water transmission in the Philip infiltration model in a manner that matched measured values well.The results of parameter estimation of approach(a) were better than those of approach(b).It was also found that the analytical precision of hydrogen and oxygen stable isotopes had a significant effect on parameter estimation using isotopic information.展开更多
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 calculation results of marine environmental design parameters obtained from different data sampling methods,model distributions,and parameter estimation methods often vary greatly.To better analyze the uncertainti...The calculation results of marine environmental design parameters obtained from different data sampling methods,model distributions,and parameter estimation methods often vary greatly.To better analyze the uncertainties in the calculation of marine environmental design parameters,a general model uncertainty assessment method is necessary.We proposed a new multivariate model uncertainty assessment method for the calculation of marine environmental design parameters.The method divides the overall model uncertainty into two categories:aleatory uncertainty and epistemic uncertainty.The aleatory uncertainty of the model is obtained by analyzing the influence of the number and the dispersion degree of samples on the information entropy of the model.The epistemic uncertainty of the model is calculated using the information entropy of the model itself and the prediction error.The advantages of this method are that it does not require many-year-observation data for the marine environmental elements,and the method can be used to analyze any specific factors that cause model uncertainty.Results show that by applying the method to the South China Sea,the aleatory uncertainty of the model increases with the number of samples and then stabilizes.A positive correlation was revealed between the dispersion of the samples and the aleatory uncertainty of the model.Both the distribution of the model and the parameter estimation results of the model have significant effects on the epistemic uncertainty of the model.When the goodness-of-fit of the model is relatively close,the best model can be selected according to the criterion of the lowest overall uncertainty of the models,which can both ensure a better model fit and avoid too much uncertainty in the model calculation results.The presented multivariate model uncertainty assessment method provides a criterion to measure the advantages and disadvantages of the marine environmental design parameter calculation model from the aspect of uncertainty,which is of great significance to analyze the uncertainties in the calculation of marine environmental design parameters and improve the accuracy of the calculation results.展开更多
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.展开更多
Considering the fact that the original two-parameter LCM model can only be used to investigate rainfall losses during the runoff period because the initial abstraction is not included, the LCM model was redefined as a...Considering the fact that the original two-parameter LCM model can only be used to investigate rainfall losses during the runoff period because the initial abstraction is not included, the LCM model was redefined as a three-parameter model, including the initial abstraction coefficient λ, the initial abstraction I_a, and the rainfall loss coefficient R. The improved LCM model is superior to the original two-parameter model, which only includes r and R, where r is the initial rainfall loss index and can be calculated with λ using the Soil Conservation Service curve number(SCS-CN) method, with r=(1+λ). The trial method was used to determine the parameter values of the improved LCM model at the watershed scale for 15 flood events in the Hongde Basin in China. The results show that larger r values are associated with smaller R values, and the parameter R ranges widely from 0.5 to 2.0. In order to improve the practicability of the LCM model, r = 0.833 with λ= 0.2 is reasonable for simplifying calculation. When the LCM model is applied to arid and semi-arid regions, rainfall without yielding runoff should be deducted from the total rainfall for more accurate estimation of rainfall-runoff.展开更多
A memory-type control chart utilizes previous information for chart construction.An example of a memory-type chart is an exponentially-weighted moving average(EWMA)control chart.The EWMA control chart is well-known an...A memory-type control chart utilizes previous information for chart construction.An example of a memory-type chart is an exponentially-weighted moving average(EWMA)control chart.The EWMA control chart is well-known and widely employed by practitioners for monitoring small and moderate process mean shifts.Meanwhile,the EWMA median chart is robust against outliers.In light of this,the economic model of the EWMA and EWMA median control charts are commonly considered.This study aims to investigate the effect of cost parameters on the out-of-control average run lengthðARL_(1)Þin implementing EWMA and EWMA median control charts.The economic model was used to compute the ARL_(1) parameter.The 14 input parameters were identified and the analysis was carried out based on the one-parameter-at-a-time basis.When the input parameters change based on a predetermined percentage,the ARL_(1) is affected.According to the results of the EWMA chart,nine input parameters had an effect andfive input parameters had no effect on the ARL_(1) parameter.Further,only seven of the 14 input parameters had an effect on the ARL_(1) of the EWMA median chart.However,the effect of each input parameter on the ARL_(1) was different.Moreover,the ARL_(1) for the EWMA median chart was smaller than the EWMA chart.This analysis is crucial to observe and determine the input parameters that have a significant impact on the ARL_(1) of the EMWA and EWMA median control charts.Hence,practitioners can obtain an overview of the influence of the input parameters on the ARL_(1) when implementing the EWMA and EWMA median control charts.展开更多
This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been pu...This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.展开更多
文摘Clastic rock reservoir is the main reservoir type in the oil and gas field.Archie formula or various conductive models developed on the basis of Archie’s formula are usually used to interpret this kind of reservoir,and the three-water model is widely used as well.However,there are many parameters in the threewater model,and some of them are difficult to determine.Most of the determination methods are based on the statistics of large amount of experimental data.In this study,the authors determine the value of the parameters of the new three-water model based on the nuclear magnetic data and the genetic optimization algorithm.The relative error between the resistivity calculated based on these parameters and the resistivity measured experimentally at 100%water content is 0.9024.The method studied in this paper can be easily applied without much experimental data.It can provide reference for other regions to determine the parameters of the new three-water model.
基金funded by the National Natural Science Foundation of China(Grant No.12272247)National Key Project(Grant No.GJXM92579)Major Research and Development Project of Metallurgical Corporation of China Ltd.in the Non-Steel Field(Grant No.2021-5).
文摘The tensile-shear interactive damage(TSID)model is a novel and powerful constitutive model for rock-like materials.This study proposes a methodology to calibrate the TSID model parameters to simulate sandstone.The basic parameters of sandstone are determined through a series of static and dynamic tests,including uniaxial compression,Brazilian disc,triaxial compression under varying confining pressures,hydrostatic compression,and dynamic compression and tensile tests with a split Hopkinson pressure bar.Based on the sandstone test results from this study and previous research,a step-by-step procedure for parameter calibration is outlined,which accounts for the categories of the strength surface,equation of state(EOS),strain rate effect,and damage.The calibrated parameters are verified through numerical tests that correspond to the experimental loading conditions.Consistency between numerical results and experimental data indicates the precision and reliability of the calibrated parameters.The methodology presented in this study is scientifically sound,straightforward,and essential for improving the TSID model.Furthermore,it has the potential to contribute to other rock constitutive models,particularly new user-defined models.
文摘Typically,a computer has infectivity as soon as it is infected.It is a reality that no antivirus programming can identify and eliminate all kinds of viruses,suggesting that infections would persevere on the Internet.To understand the dynamics of the virus propagation in a better way,a computer virus spread model with fuzzy parameters is presented in this work.It is assumed that all infected computers do not have the same contribution to the virus transmission process and each computer has a different degree of infectivity,which depends on the quantity of virus.Considering this,the parametersβandγbeing functions of the computer virus load,are considered fuzzy numbers.Using fuzzy theory helps us understand the spread of computer viruses more realistically as these parameters have fixed values in classical models.The essential features of the model,like reproduction number and equilibrium analysis,are discussed in fuzzy senses.Moreover,with fuzziness,two numerical methods,the forward Euler technique,and a nonstandard finite difference(NSFD)scheme,respectively,are developed and analyzed.In the evidence of the numerical simulations,the proposed NSFD method preserves the main features of the dynamic system.It can be considered a reliable tool to predict such types of solutions.
基金funded by the National Key Research and Development Program of China(2017YFA0605002,2017YFA0605004,and 2016YFA0601501)the National Natural Science Foundation of China(41961124007,51779145,and 41830863)“Six top talents”in Jiangsu Province(RJFW-031)。
文摘Model parameters estimation is a pivotal issue for runoff modeling in ungauged catchments.The nonlinear relationship between model parameters and catchment descriptors is a major obstacle for parameter regionalization,which is the most widely used approach.Runoff modeling was studied in 38 catchments located in the Yellow–Huai–Hai River Basin(YHHRB).The values of the Nash–Sutcliffe efficiency coefficient(NSE),coefficient of determination(R2),and percent bias(PBIAS)indicated the acceptable performance of the soil and water assessment tool(SWAT)model in the YHHRB.Nine descriptors belonging to the categories of climate,soil,vegetation,and topography were used to express the catchment characteristics related to the hydrological processes.The quantitative relationships between the parameters of the SWAT model and the catchment descriptors were analyzed by six regression-based models,including linear regression(LR)equations,support vector regression(SVR),random forest(RF),k-nearest neighbor(kNN),decision tree(DT),and radial basis function(RBF).Each of the 38 catchments was assumed to be an ungauged catchment in turn.Then,the parameters in each target catchment were estimated by the constructed regression models based on the remaining 37 donor catchments.Furthermore,the similaritybased regionalization scheme was used for comparison with the regression-based approach.The results indicated that the runoff with the highest accuracy was modeled by the SVR-based scheme in ungauged catchments.Compared with the traditional LR-based approach,the accuracy of the runoff modeling in ungauged catchments was improved by the machine learning algorithms because of the outstanding capability to deal with nonlinear relationships.The performances of different approaches were similar in humid regions,while the advantages of the machine learning techniques were more evident in arid regions.When the study area contained nested catchments,the best result was calculated with the similarity-based parameter regionalization scheme because of the high catchment density and short spatial distance.The new findings could improve flood forecasting and water resources planning in regions that lack observed data.
文摘We discuss formulas and techniques for finding maximum-likelihood estimators of parameters of autoregressive (with particular emphasis on Markov and Yule) models, computing their asymptotic variance-covariance matrix and displaying the resulting confidence regions;Monte Carlo simulation is then used to establish the accuracy of the corresponding level of confidence. The results indicate that a direct application of the Central Limit Theorem yields errors too large to be acceptable;instead, we recommend using a technique based directly on the natural logarithm of the likelihood function, verifying its substantially higher accuracy. Our study is then extended to the case of estimating only a subset of a model’s parameters, when the remaining ones (called nuisance) are of no interest to us.
文摘The effects of the milling parameters involving shape of powder particles, rotation speed, and ball-to-powder diameter (BPDR) on DEM modeling in the planetary ball mill were investigated. BPDR was varied from 1 to 10. The results revealed that the size and shape of the powder particles do not give a significant change in simulation results when BPDR attains maximum value of 10. The increasing of BPDR leads to the increase of simulation time and size. Hence, the effect of change of the powder particle shape on the calculated data size is not significant. The results also revealed that the increasing rotation speed increases impact energy between powder particles.
基金supported by the National Natural Science Foundation of China(No.42174011 and No.41874001).
文摘To estimate the parameters of the mixed additive and multiplicative(MAM)random error model using the weighted least squares iterative algorithm that requires derivation of the complex weight array,we introduce a derivative-free cat swarm optimization for parameter estimation.We embed the Powell method,which uses conjugate direction acceleration and does not need to derive the objective function,into the original cat swarm optimization to accelerate its convergence speed and search accuracy.We use the ordinary least squares,weighted least squares,original cat swarm optimization,particle swarm algorithm and improved cat swarm optimization to estimate the parameters of the straight-line fitting MAM model with lower nonlinearity and the DEM MAM model with higher nonlinearity,respectively.The experimental results show that the improved cat swarm optimization has faster convergence speed,higher search accuracy,and better stability than the original cat swarm optimization and the particle swarm algorithm.At the same time,the improved cat swarm optimization can obtain results consistent with the weighted least squares method based on the objective function only while avoiding multiple complex weight array derivations.The method in this paper provides a new idea for theoretical research on parameter estimation of MAM error models.
文摘The small-scale drilling technique can be a fast and reliable method to estimate rock strength parameters. It needs to link the operational drilling parameters and strength properties of rock. The parameters such as bit geometry, bit movement, contact frictions and crushed zone affect the estimated parameters.An analytical model considering operational drilling data and effective parameters can be used for these purposes. In this research, an analytical model was developed based on limit equilibrium of forces in a Tshaped drag bit considering the effective parameters such as bit geometry, crushed zone and contact frictions in drilling process. Based on the model, a method was used to estimate rock strength parameters such as cohesion, internal friction angle and uniaxial compressive strength of different rock types from operational drilling data. Some drilling tests were conducted by a portable and powerful drilling machine which was developed for this work. The obtained results for strength properties of different rock types from the drilling experiments based on the proposed model are in good agreement with the results of standard tests. Experimental results show that the contact friction between the cutting face and rock is close to that between bit end wearing face and rock due to the same bit material. In this case,the strength parameters, especially internal friction angle and cohesion, are estimated only by using a blunt bit drilling data and the bit bluntness does not affect the estimated results.
基金funded by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC)the authors would like to acknowledge the Niobec mine (Saint-Honoré, QuébecQuébec)
文摘Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.
基金Project(2004F052) supported by the Education Department of Liaoning Province,China
文摘By using the widely used JOHNSON-COOK model and the gradient-dependent plasticity to consider microstructural effect beyond the occurrence of shear strain localization,the distributions of local plastic shear strain and deformation in adiabatic shear band(ASB)were analyzed.The peak local plastic shear strain is proportional to the average plastic shear strain,while it is inversely proportional to the critical plastic shear strain corresponding to the peak flow shear stress.The relative plastic shear deformation between the top and base of ASB depends on the thickness of ASB and the average plastic shear strain.A parametric study was carried out to study the influence of constitutive parameters on shear strain localization.Higher values of static shear strength and work to heat conversion factor lead to lower critical plastic shear strain so that the shear localization is more apparent at the same average plastic shear strain.Higher values of strain-hardening exponent,strain rate sensitive coefficient,melting point,thermal capacity and mass density result in higher critical plastic shear strain,leading to less apparent shear localization at the same average plastic shear strain.The strain rate sensitive coefficient has a minor influence on the critical plastic shear strain,the distributions of local plastic shear strain and deformation in ASB.The effect of strain-hardening modulus on the critical plastic shear strain is not monotonous.When the maximum critical plastic shear strain is reached,the least apparent shear localization occurs.
文摘The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.
基金supported by the National Natural Science Foundation of PR China(42075130)the Postgraduate Research and Innovation Project of Jiangsu Province(1534052101133).
文摘Transmission line(TL)Parameter Identification(PI)method plays an essential role in the transmission system.The existing PI methods usually have two limitations:(1)These methods only model for single TL,and can not consider the topology connection of multiple branches for simultaneous identification.(2)Transient bad data is ignored by methods,and the random selection of terminal section data may cause the distortion of PI and have serious consequences.Therefore,a multi-task PI model considering multiple TLs’spatial constraints and massive electrical section data is proposed in this paper.The Graph Attention Network module is used to draw a single TL into a node and calculate its influence coefficient in the transmission network.Multi-Task strategy of Hard Parameter Sharing is used to identify the conductance ofmultiple branches simultaneously.Experiments show that themethod has good accuracy and robustness.Due to the consideration of spatial constraints,the method can also obtain more accurate conductance values under different training and testing conditions.
基金supported by Researchers Supporting Project number(RSP2023R455),King Saud University,Riyadh,Saudi Arabia。
文摘Herein,iron oxide/hydroxides deposits(gossans)were utilized,for the first time,in the fabrication of magnetite nanoparticles(MNPs)to load modified coal(MC).The as-synthesized MNPs@MC composite was characterized via different techniques and utilized for the Cr(Ⅵ)remediation.Experimental studies supported by theoretical treatment were applied to offer a new overview of the Cr(Ⅵ)adsorption geometry and mechanism at 25-45℃.Experimental results suggested that the Cr(Ⅵ)uptake was mainly governed by adsorption-reduction coupled mechanism.The Langmuir model fitted well the Cr(Ⅵ)adsorption data with maximum adsorption capacities extended from 115.24 to 129.63 mg·g^(-1).Theoretical calculations indicated that Cr(Ⅵ)ions were adsorbed on the MNPs@MC following the theory of the advanced monolayer statistical model.The number of ions removed per site ranged from 1.88 to1.23 suggesting the involvement of vertical geometry and multi-ionic mechanism at all temperatures.The increment of the active sites density and the adsorption capacity at saturation with improving temperature reflected an endothermic process.Energetically,the Cr(Ⅵ)adsorption was controlled by physical forces as the adsorption energies were less than 40 kJ·mol^(-1).The calculated free enthalpy,entropy.and internal energy explained the spontaneous nature and the viability of Cr(Ⅵ)adsorption on the MNPs@MC adsorbent.These results offer a new approach in utilizing the iron-rich deposits as gossans in the preparation of magnetic and low-cost adsorbents for wastewater remediation.
基金supported by the National Natural Science Foundation of China(Grant No.51279057)
文摘Minimizing parameter uncertainty is crucial in the application of hydrologic models.Isotopic information in various hydrologic components of the water cycle can expand our knowledge of the dynamics of water flow in the system,provide additional information for parameter estimation,and improve parameter identifiability.This study combined the Philip infiltration model with an isotopic mixing model using an isotopic mass balance approach for estimating parameters in the Philip infiltration model.Two approaches to parameter estimation were compared:(a) using isotopic information to determine the soil water transmission and then hydrologic information to estimate the soil sorptivity,and(b) using hydrologic information to determine the soil water transmission and the soil sorptivity.Results of parameter estimation were verified through a rainfall infiltration experiment in a laboratory under rainfall with constant isotopic compositions and uniform initial soil water content conditions.Experimental results showed that approach(a),using isotopic and hydrologic information,estimated the soil water transmission in the Philip infiltration model in a manner that matched measured values well.The results of parameter estimation of approach(a) were better than those of approach(b).It was also found that the analytical precision of hydrogen and oxygen stable isotopes had a significant effect on parameter estimation using isotopic information.
基金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.
基金Supported by the National Natural Science Foundation of China(No.52071306)the Natural Science Foundation of Shandong Province(No.ZR2019MEE050)。
文摘The calculation results of marine environmental design parameters obtained from different data sampling methods,model distributions,and parameter estimation methods often vary greatly.To better analyze the uncertainties in the calculation of marine environmental design parameters,a general model uncertainty assessment method is necessary.We proposed a new multivariate model uncertainty assessment method for the calculation of marine environmental design parameters.The method divides the overall model uncertainty into two categories:aleatory uncertainty and epistemic uncertainty.The aleatory uncertainty of the model is obtained by analyzing the influence of the number and the dispersion degree of samples on the information entropy of the model.The epistemic uncertainty of the model is calculated using the information entropy of the model itself and the prediction error.The advantages of this method are that it does not require many-year-observation data for the marine environmental elements,and the method can be used to analyze any specific factors that cause model uncertainty.Results show that by applying the method to the South China Sea,the aleatory uncertainty of the model increases with the number of samples and then stabilizes.A positive correlation was revealed between the dispersion of the samples and the aleatory uncertainty of the model.Both the distribution of the model and the parameter estimation results of the model have significant effects on the epistemic uncertainty of the model.When the goodness-of-fit of the model is relatively close,the best model can be selected according to the criterion of the lowest overall uncertainty of the models,which can both ensure a better model fit and avoid too much uncertainty in the model calculation results.The presented multivariate model uncertainty assessment method provides a criterion to measure the advantages and disadvantages of the marine environmental design parameter calculation model from the aspect of uncertainty,which is of great significance to analyze the uncertainties in the calculation of marine environmental design parameters and improve the accuracy of the calculation results.
基金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(Grants No.41271048 and 41330529)
文摘Considering the fact that the original two-parameter LCM model can only be used to investigate rainfall losses during the runoff period because the initial abstraction is not included, the LCM model was redefined as a three-parameter model, including the initial abstraction coefficient λ, the initial abstraction I_a, and the rainfall loss coefficient R. The improved LCM model is superior to the original two-parameter model, which only includes r and R, where r is the initial rainfall loss index and can be calculated with λ using the Soil Conservation Service curve number(SCS-CN) method, with r=(1+λ). The trial method was used to determine the parameter values of the improved LCM model at the watershed scale for 15 flood events in the Hongde Basin in China. The results show that larger r values are associated with smaller R values, and the parameter R ranges widely from 0.5 to 2.0. In order to improve the practicability of the LCM model, r = 0.833 with λ= 0.2 is reasonable for simplifying calculation. When the LCM model is applied to arid and semi-arid regions, rainfall without yielding runoff should be deducted from the total rainfall for more accurate estimation of rainfall-runoff.
基金funded by the Universiti Kebangsaan Malaysia,Geran Galakan Penyelidikan,GGP-2020-040.
文摘A memory-type control chart utilizes previous information for chart construction.An example of a memory-type chart is an exponentially-weighted moving average(EWMA)control chart.The EWMA control chart is well-known and widely employed by practitioners for monitoring small and moderate process mean shifts.Meanwhile,the EWMA median chart is robust against outliers.In light of this,the economic model of the EWMA and EWMA median control charts are commonly considered.This study aims to investigate the effect of cost parameters on the out-of-control average run lengthðARL_(1)Þin implementing EWMA and EWMA median control charts.The economic model was used to compute the ARL_(1) parameter.The 14 input parameters were identified and the analysis was carried out based on the one-parameter-at-a-time basis.When the input parameters change based on a predetermined percentage,the ARL_(1) is affected.According to the results of the EWMA chart,nine input parameters had an effect andfive input parameters had no effect on the ARL_(1) parameter.Further,only seven of the 14 input parameters had an effect on the ARL_(1) of the EWMA median chart.However,the effect of each input parameter on the ARL_(1) was different.Moreover,the ARL_(1) for the EWMA median chart was smaller than the EWMA chart.This analysis is crucial to observe and determine the input parameters that have a significant impact on the ARL_(1) of the EMWA and EWMA median control charts.Hence,practitioners can obtain an overview of the influence of the input parameters on the ARL_(1) when implementing the EWMA and EWMA median control charts.
文摘This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.
