We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived ...We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.展开更多
Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying i...Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying in sizes and lifespans,significantly influence the distribution of fluid velocities within the flow.Subsequently,the rapid velocity fluctuations in highly turbulent flows lead to elevated shear and normal stress levels.For this reason,to meticulously study these dynamics,more often than not,physical modeling is employed for studying the impact of turbulent flows on the stability and longevity of nearby structures.Despite the effectiveness of physical modeling,various monitoring challenges arise,including flow disruption,the necessity for concurrent gauging at multiple locations,and the duration of measurements.Addressing these challenges,image velocimetry emerges as an ideal method in fluid mechanics,particularly for studying turbulent flows.To account for measurement duration,a probabilistic approach utilizing a probability density function(PDF)is suggested to mitigate uncertainty in estimated average and maximum values.However,it becomes evident that deriving the PDF is not straightforward for all turbulence-induced stresses.In response,this study proposes a novel approach by combining image velocimetry with a stochastic model to provide a generic yet accurate description of flow dynamics in such applications.This integration enables an approach based on the probability of failure,facilitating a more comprehensive analysis of turbulent flows.Such an approach is essential for estimating both short-and long-term stresses on hydraulic constructions under assessment.展开更多
Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currentl...Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currently,least squares(LS)+auto-regressive(AR)hybrid method is one of the main techniques of PM prediction.Besides,the weighted LS+AR hybrid method performs well for PM short-term prediction.However,the corresponding covariance information of LS fitting residuals deserves further exploration in the AR model.In this study,we have derived a modified stochastic model for the LS+AR hybrid method,namely the weighted LS+weighted AR hybrid method.By using the PM data products of IERS EOP 14 C04,the numerical results indicate that for PM short-term forecasting,the proposed weighted LS+weighted AR hybrid method shows an advantage over both the LS+AR hybrid method and the weighted LS+AR hybrid method.Compared to the mean absolute errors(MAEs)of PMX/PMY sho rt-term prediction of the LS+AR hybrid method and the weighted LS+AR hybrid method,the weighted LS+weighted AR hybrid method shows average improvements of 6.61%/12.08%and 0.24%/11.65%,respectively.Besides,for the slopes of the linear regression lines fitted to the errors of each method,the growth of the prediction error of the proposed method is slower than that of the other two methods.展开更多
The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whiteno...The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whitenoise and Gaussian colored noise are introduced into a tumor growth model under immune surveillance. Asfollows, the long-time evolution of the tumor characterized by the Stationary Probability Density (SPD) and MFPTis obtained in theory on the basis of the Approximated Fokker-Planck Equation (AFPE). Herein the recurrenceof the tumor from the extinction state to the tumor-present state is more concerned in this paper. A moreefficient algorithmof Back-Propagation Neural Network (BPNN) is utilized in order to testify the correction of thetheoretical SPDandMFPT.With the existence of aweak signal, the functional relationship between Signal-to-NoiseRatio (SNR), noise intensities and correlation time is also studied. Numerical results show that both multiplicativeGaussian colored noise and additive Gaussian white noise can promote the extinction of the tumors, and themultiplicative Gaussian colored noise can lead to the resonance-like peak on MFPT curves, while the increasingintensity of the additiveGaussian white noise results in theminimum of MFPT. In addition, the correlation timesare negatively correlated with MFPT. As for the SNR, we find the intensities of both the Gaussian white noise andthe Gaussian colored noise, as well as their correlation intensity can induce SR. Especially, SNR is monotonouslyincreased in the case ofGaussian white noisewith the change of the correlation time.At last, the optimal parametersin BPNN structure are analyzed for MFPT from three aspects: the penalty factors, the number of neural networklayers and the number of nodes in each layer.展开更多
A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epi...A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epidemiological impact of vaccine booster doses on the co-dynamics of viral hepatitis B and COVID-19.The model is fitted to real COVID-19 data from Pakistan.The proposed model incorporates logistic growth and saturated incidence functions.Rigorous analyses using the tools of stochastic calculus,are performed to study appropriate conditions for the existence of unique global solutions,stationary distribution in the sense of ergodicity and disease extinction.The stochastic threshold estimated from the data fitting is given by:R_(0)^(S)=3.0651.Numerical assessments are implemented to illustrate the impact of double-dose vaccination and saturated incidence functions on the dynamics of both diseases.The effects of stochastic white noise intensities are also highlighted.展开更多
The multiple-input multiple-output(MIMO)-enabled beamforming technology offers great data rate and channel quality for next-generation communication.In this paper,we propose a beam channel model and enable it with tim...The multiple-input multiple-output(MIMO)-enabled beamforming technology offers great data rate and channel quality for next-generation communication.In this paper,we propose a beam channel model and enable it with time-varying simulation capability by adopting the stochastic geometry theory.First,clusters are generated located within transceivers'beam ranges based on the Mate?rn hardcore Poisson cluster process.The line-of-sight,singlebounce,and double-bounce components are calculated when generating the complex channel impulse response.Furthermore,we elaborate on the expressions of channel links based on the propagation-graph theory.A birth-death process consisting of the effects of beams and cluster velocities is also formulated.Numerical simulation results prove that the proposed model can capture the channel non-stationarity.Besides,the non-reciprocal beam patterns yield severe channel dispersion compared to the reciprocal patterns.展开更多
Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challeng...Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.展开更多
Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate b...Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.展开更多
In this paper we study optimal advertising problems that model the introduction of a new product into the market in the presence of carryover effects of the advertisement and with memory effects in the level of goodwi...In this paper we study optimal advertising problems that model the introduction of a new product into the market in the presence of carryover effects of the advertisement and with memory effects in the level of goodwill. In particular, we let the dynamics of the product goodwill to depend on the past, and also on past advertising efforts. We treat the problem by means of the stochastic Pontryagin maximum principle, that here is considered for a class of problems where in the state equation either the state or the control depend on the past. Moreover the control acts on the martingale term and the space of controls U can be chosen to be non-convex but now the space of controls U can be chosen to be non-convex. The maximum principle is thus formulated using a first-order adjoint Backward Stochastic Differential Equations (BSDEs), which can be explicitly computed due to the specific characteristics of the model, and a second-order adjoint relation.展开更多
In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochast...In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochastic system by computing the Lyapunov exponent of the linearized system. Further, the global stability of the stochastic model is analyzed based on the singular boundary theory. Moreover, we prove that the model undergoes a Hopf bifurcation and a pitchfork bifurcation. Finally, several numerical examples are provided to illustrate the theoretical results. .展开更多
Ross’ epidemic model describing the transmission of malaria uses two classes of infection, one for humans and one for mosquitoes. This paper presents a stochastic extension of a deterministic vector-borne epidemic mo...Ross’ epidemic model describing the transmission of malaria uses two classes of infection, one for humans and one for mosquitoes. This paper presents a stochastic extension of a deterministic vector-borne epidemic model based only on the class of human infectious. The consistency of the model is established by proving that the stochastic delay differential equation describing the model has a unique positive global solution. The extinction of the disease is studied through the analysis of the stability of the disease-free equilibrium state and the persistence of the model. Finally, we introduce some numerical simulations to illustrate the obtained results.展开更多
Aiming at the problem that the intermediate potential part of the traditional bistable stochastic resonance model cannot be adjusted independently, a new composite stochastic resonance(NCSR) model is proposed by combi...Aiming at the problem that the intermediate potential part of the traditional bistable stochastic resonance model cannot be adjusted independently, a new composite stochastic resonance(NCSR) model is proposed by combining the Woods–Saxon(WS) model and the improved piecewise bistable model. The model retains the characteristics of the independent parameters of WS model and the improved piecewise model has no output saturation, all the parameters in the new model have no coupling characteristics. Under α stable noise environment, the new model is used to detect periodic signal and aperiodic signal, the detection results indicate that the new model has higher noise utilization and better detection effect.Finally, the new model is applied to image denoising, the results showed that under the same conditions, the output peak signal-to-noise ratio(PSNR) and the correlation number of NCSR method is higher than that of other commonly used linear denoising methods and improved piecewise SR methods, the effectiveness of the new model is verified.展开更多
In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean ...In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean to the west and the highly saline waters of the Saloum River to the north. In this estuarine and deltaic zone with its very low relief, the hydraulic loads in the water tables are generally close to zero or even negative, creating a reversal of the natural flow and encouraging saline intrusion into this system, which makes it very vulnerable. This study concerns the implementation of a numerical model of saline intrusion to provide a better understanding of the vulnerability of the water table by analyzing the variability of the freshwater/saltwater interface. The Modflow-2005 code is used to simulate saline intrusion using the SWI2 module, coupled with the GRASS (Geographic Resources Analysis Support System) software under the Linux operating system with the steep interface approach. The probable expansion of the wedge is studied in three scenarios, taking into account its position relative to the bedrock at 1 m, 5 m and 10 m. Simulations carried out under imposed potential and river conditions, based on variations in groundwater reserves using two effective porosity values, 10−1 and 10−2, show that the water table is highly vulnerable in the northwest sector. The probable expansion of the wedge increases as the storage coefficient decreases and is more marked with river conditions in the areas surrounding the Saloum River, reaching 6 km with a probability of 1. The probability of the wedge reaching a certain degree of expansion decreases from 1 to 0.5, and then cancels out as it moves inland. The probable position of the wedge is limited to 500 m or even 1 km depending on the corner around the coast to the southwest and in the southern zone. This modelling, carried out under natural conditions, will be developed further, taking into account climatic parameters and pumping from wells and boreholes.展开更多
An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is...An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is governed by Heston's stochastic volatility(SV)model.With the objective of maximizing the expected index utility of the terminal wealth of the insurance company,by using the classical tools of stochastic optimal control,the explicit expressions for optimal strategies and optimal value functions are derived.An interesting conclusion is found that it is better to buy one reinsurance than two under the assumption of this paper.Moreover,some numerical simulations and sensitivity analysis are provided.展开更多
Markov modeling of HIV/AIDS progression was done under the assumption that the state holding time (waiting time) had a constant hazard. This paper discusses the properties of the hazard function of the Exponential dis...Markov modeling of HIV/AIDS progression was done under the assumption that the state holding time (waiting time) had a constant hazard. This paper discusses the properties of the hazard function of the Exponential distributions and its modifications namely;Parameter proportion hazard (PH) and Accelerated failure time models (AFT) and their effectiveness in modeling the state holding time in Markov modeling of HIV/AIDS progression with and without risk factors. Patients were categorized by gender and age with female gender being the baseline. Data simulated using R software was fitted to each model, and the model parameters were estimated. The estimated P and Z values were then used to test the null hypothesis that the state waiting time data followed an Exponential distribution. Model identification criteria;Akaike information criteria (AIC), Bayesian information criteria (BIC), log-likelihood (LL), and R2 were used to evaluate the performance of the models. For the Survival Regression model, P and Z values supported the non-rejection of the null hypothesis for mixed gender without interaction and supported the rejection of the same for mixed gender with interaction term and males aged 50 - 60 years. Both Parameters supported the non-rejection of the null hypothesis in the rest of the age groups. For Gender male with interaction both P and Z values supported rejection in all the age groups except the age group 20 - 30 years. For Cox Proportional hazard and AFT models, both P and Z values supported the non-rejection of the null hypothesis across all age groups. The P-values for the three models supported different decisions for and against the Null hypothesis with AFT and Cox values supporting similar decisions in most of the age groups. Among the models considered, the regression assumption provided a superior fit based on (AIC), (BIC), (LL), and R2 Model identification criteria. This was particularly evident in age and gender subgroups where the data exhibited non-proportional hazards and violated the assumptions required for the Cox Proportional Hazard model. Moreover, the simplicity of the regression model, along with its ability to capture essential state transitions without over fitting, made it a more appropriate choice.展开更多
The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic ...The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic flow analysis,simulation,autonomous vehicle development,etc.Two-dimensional(2D)vehicular movement is basically stochastic and is the result of interactions between a driver's behavior and a vehicle's characteristics.Current microscopic models either neglect 2D noise,or overlook vehicle dynamics.The modeling capabilities,thus,are limited,so that stochastic lateral movement cannot be reproduced.The present research extends an intelligent driver model(IDM)by explicitly considering both vehicle dynamics and 2D noises to formulate a stochastic 2D IDM model,with vehicle dynamics based on the stochastic differential equation(SDE)theory.Control inputs from the vehicle include the steer rate and longitudinal acceleration,both of which are developed based on an idea from a traditional intelligent driver model.The stochastic stability condition is analyzed on the basis of Lyapunov theory.Numerical analysis is used to assess the two cases:(i)when a vehicle accelerates from a standstill and(ii)when a platoon of vehicles follow a leader with a stop-and-go speed profile,the formation of congestion and subsequent dispersion are simulated.The results show that the model can reproduce the stochastic 2D trajectories of the vehicle and the marginal distribution of lateral movement.The proposed model can be used in both a simulation platform and a behavioral analysis of a human driver in traffic flow.展开更多
Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the...Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the evolution.However,although there has been a lot of work on stochastic epidemic models,most of them focus mainly on qualitative properties,which makes us somewhat ignore the original meaning of the parameter value.In this paper we extend the classic susceptible-infectious-removed(SIR)epidemic model by adding a white noise excitation and then we utilize the large deviation theory to quantitatively study the long-term coexistence exit problem with epidemic.Finally,in order to extend the meaning of parameters in the corresponding deterministic system,we tentatively introduce two new thresholds which then prove rational.展开更多
Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a n...Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a new non-linear generalized model to describe Cyber-Physical Systems.This model includes unknown multivariable discrete and continuous-time functions and different multiplicative noises to represent the evolution of physical processes and randomeffects in the physical and computationalworlds.Besides,the digitalization stage in hardware devices is represented too.Attackers and most critical sparse sensor attacks are described through a stochastic process.The reconstruction and protectionmechanisms are based on aweighted stochasticmodel.Error probability in data samples is estimated through different indicators commonly employed in non-linear dynamics(such as the Fourier transform,first-return maps,or the probability density function).A decision algorithm calculates the final reconstructed value considering the previous error probability.An experimental validation based on simulation tools and real deployments is also carried out.Both,the new technology performance and scalability are studied.Results prove that the proposed solution protects Cyber-Physical Systems against up to 92%of attacks and perturbations,with a computational delay below 2.5 s.The proposed model shows a linear complexity,as recursive or iterative structures are not employed,just algebraic and probabilistic functions.In conclusion,the new model and reconstructionmechanism can protect successfully Cyber-Physical Systems against sparse sensor attacks,even in dense or pervasive deployments and scenarios.展开更多
The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system ...The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.展开更多
This paper presents a study on a new rumor propagation model with nonlinear propagation rate and secondary propagation rate. We divide the total population into three groups, the ignorant, the spreader and the aware. ...This paper presents a study on a new rumor propagation model with nonlinear propagation rate and secondary propagation rate. We divide the total population into three groups, the ignorant, the spreader and the aware. The nonlinear incidence rate describes the psychological impact of certain serious rumors on social groups when the number of individuals spreading rumors becomes larger. The main contributions of this work are the development of a new rumor propagation model and some results of deterministic and stochastic analysis of the rumor propagation model. The results show the influence of nonlinear propagation rate and stochastic fluctuation on the dynamic behavior of the rumor propagation model by using Lyapunov function method and stochastic related knowledge. Numerical examples and simulation results are given to illustrate the results obtained.展开更多
基金Equinor for financing the R&D projectthe Institute of Science and Technology of Petroleum Geophysics of Brazil for supporting this research。
文摘We apply stochastic seismic inversion and Bayesian facies classification for porosity modeling and igneous rock identification in the presalt interval of the Santos Basin. This integration of seismic and well-derived information enhances reservoir characterization. Stochastic inversion and Bayesian classification are powerful tools because they permit addressing the uncertainties in the model. We used the ES-MDA algorithm to achieve the realizations equivalent to the percentiles P10, P50, and P90 of acoustic impedance, a novel method for acoustic inversion in presalt. The facies were divided into five: reservoir 1,reservoir 2, tight carbonates, clayey rocks, and igneous rocks. To deal with the overlaps in acoustic impedance values of facies, we included geological information using a priori probability, indicating that structural highs are reservoir-dominated. To illustrate our approach, we conducted porosity modeling using facies-related rock-physics models for rock-physics inversion in an area with a well drilled in a coquina bank and evaluated the thickness and extension of an igneous intrusion near the carbonate-salt interface. The modeled porosity and the classified seismic facies are in good agreement with the ones observed in the wells. Notably, the coquinas bank presents an improvement in the porosity towards the top. The a priori probability model was crucial for limiting the clayey rocks to the structural lows. In Well B, the hit rate of the igneous rock in the three scenarios is higher than 60%, showing an excellent thickness-prediction capability.
文摘Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying in sizes and lifespans,significantly influence the distribution of fluid velocities within the flow.Subsequently,the rapid velocity fluctuations in highly turbulent flows lead to elevated shear and normal stress levels.For this reason,to meticulously study these dynamics,more often than not,physical modeling is employed for studying the impact of turbulent flows on the stability and longevity of nearby structures.Despite the effectiveness of physical modeling,various monitoring challenges arise,including flow disruption,the necessity for concurrent gauging at multiple locations,and the duration of measurements.Addressing these challenges,image velocimetry emerges as an ideal method in fluid mechanics,particularly for studying turbulent flows.To account for measurement duration,a probabilistic approach utilizing a probability density function(PDF)is suggested to mitigate uncertainty in estimated average and maximum values.However,it becomes evident that deriving the PDF is not straightforward for all turbulence-induced stresses.In response,this study proposes a novel approach by combining image velocimetry with a stochastic model to provide a generic yet accurate description of flow dynamics in such applications.This integration enables an approach based on the probability of failure,facilitating a more comprehensive analysis of turbulent flows.Such an approach is essential for estimating both short-and long-term stresses on hydraulic constructions under assessment.
基金supported by National Natural Science Foundation of China,China(No.42004016)HuBei Natural Science Fund,China(No.2020CFB329)+1 种基金HuNan Natural Science Fund,China(No.2023JJ60559,2023JJ60560)the State Key Laboratory of Geodesy and Earth’s Dynamics self-deployment project,China(No.S21L6101)。
文摘Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currently,least squares(LS)+auto-regressive(AR)hybrid method is one of the main techniques of PM prediction.Besides,the weighted LS+AR hybrid method performs well for PM short-term prediction.However,the corresponding covariance information of LS fitting residuals deserves further exploration in the AR model.In this study,we have derived a modified stochastic model for the LS+AR hybrid method,namely the weighted LS+weighted AR hybrid method.By using the PM data products of IERS EOP 14 C04,the numerical results indicate that for PM short-term forecasting,the proposed weighted LS+weighted AR hybrid method shows an advantage over both the LS+AR hybrid method and the weighted LS+AR hybrid method.Compared to the mean absolute errors(MAEs)of PMX/PMY sho rt-term prediction of the LS+AR hybrid method and the weighted LS+AR hybrid method,the weighted LS+weighted AR hybrid method shows average improvements of 6.61%/12.08%and 0.24%/11.65%,respectively.Besides,for the slopes of the linear regression lines fitted to the errors of each method,the growth of the prediction error of the proposed method is slower than that of the other two methods.
基金National Natural Science Foundation of China(Nos.12272283,12172266).
文摘The Mean First-Passage Time (MFPT) and Stochastic Resonance (SR) of a stochastic tumor-immune model withnoise perturbation are discussed in this paper. Firstly, considering environmental perturbation, Gaussian whitenoise and Gaussian colored noise are introduced into a tumor growth model under immune surveillance. Asfollows, the long-time evolution of the tumor characterized by the Stationary Probability Density (SPD) and MFPTis obtained in theory on the basis of the Approximated Fokker-Planck Equation (AFPE). Herein the recurrenceof the tumor from the extinction state to the tumor-present state is more concerned in this paper. A moreefficient algorithmof Back-Propagation Neural Network (BPNN) is utilized in order to testify the correction of thetheoretical SPDandMFPT.With the existence of aweak signal, the functional relationship between Signal-to-NoiseRatio (SNR), noise intensities and correlation time is also studied. Numerical results show that both multiplicativeGaussian colored noise and additive Gaussian white noise can promote the extinction of the tumors, and themultiplicative Gaussian colored noise can lead to the resonance-like peak on MFPT curves, while the increasingintensity of the additiveGaussian white noise results in theminimum of MFPT. In addition, the correlation timesare negatively correlated with MFPT. As for the SNR, we find the intensities of both the Gaussian white noise andthe Gaussian colored noise, as well as their correlation intensity can induce SR. Especially, SNR is monotonouslyincreased in the case ofGaussian white noisewith the change of the correlation time.At last, the optimal parametersin BPNN structure are analyzed for MFPT from three aspects: the penalty factors, the number of neural networklayers and the number of nodes in each layer.
文摘A patient co-infected with COVID-19 and viral hepatitis B can be atmore risk of severe complications than the one infected with a single infection.This study develops a comprehensive stochastic model to assess the epidemiological impact of vaccine booster doses on the co-dynamics of viral hepatitis B and COVID-19.The model is fitted to real COVID-19 data from Pakistan.The proposed model incorporates logistic growth and saturated incidence functions.Rigorous analyses using the tools of stochastic calculus,are performed to study appropriate conditions for the existence of unique global solutions,stationary distribution in the sense of ergodicity and disease extinction.The stochastic threshold estimated from the data fitting is given by:R_(0)^(S)=3.0651.Numerical assessments are implemented to illustrate the impact of double-dose vaccination and saturated incidence functions on the dynamics of both diseases.The effects of stochastic white noise intensities are also highlighted.
基金supported by the National Key R&D Program of China under grant 2020YFB1804901the National Natural Science Foundation of China under grant 62341102。
文摘The multiple-input multiple-output(MIMO)-enabled beamforming technology offers great data rate and channel quality for next-generation communication.In this paper,we propose a beam channel model and enable it with time-varying simulation capability by adopting the stochastic geometry theory.First,clusters are generated located within transceivers'beam ranges based on the Mate?rn hardcore Poisson cluster process.The line-of-sight,singlebounce,and double-bounce components are calculated when generating the complex channel impulse response.Furthermore,we elaborate on the expressions of channel links based on the propagation-graph theory.A birth-death process consisting of the effects of beams and cluster velocities is also formulated.Numerical simulation results prove that the proposed model can capture the channel non-stationarity.Besides,the non-reciprocal beam patterns yield severe channel dispersion compared to the reciprocal patterns.
基金Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Project under Grant Number RGP2/302/45supported by the Deanship of Scientific Research,Vice Presidency forGraduate Studies and Scientific Research,King Faisal University,Saudi Arabia(Grant Number A426).
文摘Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.
基金supported by the National Natural Science Foundation of China(Grant No.52109010)the Postdoctoral Science Foundation of China(Grant No.2021M701047)the China National Postdoctoral Program for Innovative Talents(Grant No.BX20200113).
文摘Copula functions have been widely used in stochastic simulation and prediction of streamflow.However,existing models are usually limited to single two-dimensional or three-dimensional copulas with the same bivariate block for all months.To address this limitation,this study developed a mixed D-vine copula-based conditional quantile model that can capture temporal correlations.This model can generate streamflow by selecting different historical streamflow variables as the conditions for different months and by exploiting the conditional quantile functions of streamflows in different months with mixed D-vine copulas.The up-to-down sequential method,which couples the maximum weight approach with the Akaike information criteria and the maximum likelihood approach,was used to determine the structures of multivariate Dvine copulas.The developed model was used in a case study to synthesize the monthly streamflow at the Tangnaihai hydrological station,the inflow control station of the Longyangxia Reservoir in the Yellow River Basin.The results showed that the developed model outperformed the commonly used bivariate copula model in terms of the performance in simulating the seasonality and interannual variability of streamflow.This model provides useful information for water-related natural hazard risk assessment and integrated water resources management and utilization.
文摘In this paper we study optimal advertising problems that model the introduction of a new product into the market in the presence of carryover effects of the advertisement and with memory effects in the level of goodwill. In particular, we let the dynamics of the product goodwill to depend on the past, and also on past advertising efforts. We treat the problem by means of the stochastic Pontryagin maximum principle, that here is considered for a class of problems where in the state equation either the state or the control depend on the past. Moreover the control acts on the martingale term and the space of controls U can be chosen to be non-convex but now the space of controls U can be chosen to be non-convex. The maximum principle is thus formulated using a first-order adjoint Backward Stochastic Differential Equations (BSDEs), which can be explicitly computed due to the specific characteristics of the model, and a second-order adjoint relation.
文摘In this paper, we investigate an SIS model with treatment and immigration. Firstly, the two-dimensional model is simplified by using the stochastic averaging method. Then, we derive the local stability of the stochastic system by computing the Lyapunov exponent of the linearized system. Further, the global stability of the stochastic model is analyzed based on the singular boundary theory. Moreover, we prove that the model undergoes a Hopf bifurcation and a pitchfork bifurcation. Finally, several numerical examples are provided to illustrate the theoretical results. .
文摘Ross’ epidemic model describing the transmission of malaria uses two classes of infection, one for humans and one for mosquitoes. This paper presents a stochastic extension of a deterministic vector-borne epidemic model based only on the class of human infectious. The consistency of the model is established by proving that the stochastic delay differential equation describing the model has a unique positive global solution. The extinction of the disease is studied through the analysis of the stability of the disease-free equilibrium state and the persistence of the model. Finally, we introduce some numerical simulations to illustrate the obtained results.
基金Project supported by the National Natural Science Foundation of China (Grant No.62371388)the Key Research and Development Projects in Shaanxi Province,China (Grant No.2023-YBGY-044)。
文摘Aiming at the problem that the intermediate potential part of the traditional bistable stochastic resonance model cannot be adjusted independently, a new composite stochastic resonance(NCSR) model is proposed by combining the Woods–Saxon(WS) model and the improved piecewise bistable model. The model retains the characteristics of the independent parameters of WS model and the improved piecewise model has no output saturation, all the parameters in the new model have no coupling characteristics. Under α stable noise environment, the new model is used to detect periodic signal and aperiodic signal, the detection results indicate that the new model has higher noise utilization and better detection effect.Finally, the new model is applied to image denoising, the results showed that under the same conditions, the output peak signal-to-noise ratio(PSNR) and the correlation number of NCSR method is higher than that of other commonly used linear denoising methods and improved piecewise SR methods, the effectiveness of the new model is verified.
文摘In the Saloum region of central-western Senegal, water needs are essentially met by tapping an underground aquifer associated with the sandy-clay formations of the Continental Terminal, in contact with both the ocean to the west and the highly saline waters of the Saloum River to the north. In this estuarine and deltaic zone with its very low relief, the hydraulic loads in the water tables are generally close to zero or even negative, creating a reversal of the natural flow and encouraging saline intrusion into this system, which makes it very vulnerable. This study concerns the implementation of a numerical model of saline intrusion to provide a better understanding of the vulnerability of the water table by analyzing the variability of the freshwater/saltwater interface. The Modflow-2005 code is used to simulate saline intrusion using the SWI2 module, coupled with the GRASS (Geographic Resources Analysis Support System) software under the Linux operating system with the steep interface approach. The probable expansion of the wedge is studied in three scenarios, taking into account its position relative to the bedrock at 1 m, 5 m and 10 m. Simulations carried out under imposed potential and river conditions, based on variations in groundwater reserves using two effective porosity values, 10−1 and 10−2, show that the water table is highly vulnerable in the northwest sector. The probable expansion of the wedge increases as the storage coefficient decreases and is more marked with river conditions in the areas surrounding the Saloum River, reaching 6 km with a probability of 1. The probability of the wedge reaching a certain degree of expansion decreases from 1 to 0.5, and then cancels out as it moves inland. The probable position of the wedge is limited to 500 m or even 1 km depending on the corner around the coast to the southwest and in the southern zone. This modelling, carried out under natural conditions, will be developed further, taking into account climatic parameters and pumping from wells and boreholes.
基金National Natural Science Foundation of China(No.62073071)Fundamental Research Funds for the Central Universities and Graduate Student Innovation Fund of Donghua University,China(No.CUSF-DH-D-2021045)。
文摘An optimal quota-share and excess-of-loss reinsurance and investment problem is studied for an insurer who is allowed to invest in a risk-free asset and a risky asset.Especially the price process of the risky asset is governed by Heston's stochastic volatility(SV)model.With the objective of maximizing the expected index utility of the terminal wealth of the insurance company,by using the classical tools of stochastic optimal control,the explicit expressions for optimal strategies and optimal value functions are derived.An interesting conclusion is found that it is better to buy one reinsurance than two under the assumption of this paper.Moreover,some numerical simulations and sensitivity analysis are provided.
文摘Markov modeling of HIV/AIDS progression was done under the assumption that the state holding time (waiting time) had a constant hazard. This paper discusses the properties of the hazard function of the Exponential distributions and its modifications namely;Parameter proportion hazard (PH) and Accelerated failure time models (AFT) and their effectiveness in modeling the state holding time in Markov modeling of HIV/AIDS progression with and without risk factors. Patients were categorized by gender and age with female gender being the baseline. Data simulated using R software was fitted to each model, and the model parameters were estimated. The estimated P and Z values were then used to test the null hypothesis that the state waiting time data followed an Exponential distribution. Model identification criteria;Akaike information criteria (AIC), Bayesian information criteria (BIC), log-likelihood (LL), and R2 were used to evaluate the performance of the models. For the Survival Regression model, P and Z values supported the non-rejection of the null hypothesis for mixed gender without interaction and supported the rejection of the same for mixed gender with interaction term and males aged 50 - 60 years. Both Parameters supported the non-rejection of the null hypothesis in the rest of the age groups. For Gender male with interaction both P and Z values supported rejection in all the age groups except the age group 20 - 30 years. For Cox Proportional hazard and AFT models, both P and Z values supported the non-rejection of the null hypothesis across all age groups. The P-values for the three models supported different decisions for and against the Null hypothesis with AFT and Cox values supporting similar decisions in most of the age groups. Among the models considered, the regression assumption provided a superior fit based on (AIC), (BIC), (LL), and R2 Model identification criteria. This was particularly evident in age and gender subgroups where the data exhibited non-proportional hazards and violated the assumptions required for the Cox Proportional Hazard model. Moreover, the simplicity of the regression model, along with its ability to capture essential state transitions without over fitting, made it a more appropriate choice.
基金Project supported by the National Key Research and Development Program of China(Grant No.2021YFE0194400)the National Natural Science Foundation of China(Grant Nos.52272314 and 52131202)+1 种基金the Fund for Humanities and Social Science from the Ministry of Education of China(Grant No.21YJCZH116)the Public Welfare Scientific Research Project(Grant No.LGF22E080007)。
文摘The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic flow analysis,simulation,autonomous vehicle development,etc.Two-dimensional(2D)vehicular movement is basically stochastic and is the result of interactions between a driver's behavior and a vehicle's characteristics.Current microscopic models either neglect 2D noise,or overlook vehicle dynamics.The modeling capabilities,thus,are limited,so that stochastic lateral movement cannot be reproduced.The present research extends an intelligent driver model(IDM)by explicitly considering both vehicle dynamics and 2D noises to formulate a stochastic 2D IDM model,with vehicle dynamics based on the stochastic differential equation(SDE)theory.Control inputs from the vehicle include the steer rate and longitudinal acceleration,both of which are developed based on an idea from a traditional intelligent driver model.The stochastic stability condition is analyzed on the basis of Lyapunov theory.Numerical analysis is used to assess the two cases:(i)when a vehicle accelerates from a standstill and(ii)when a platoon of vehicles follow a leader with a stop-and-go speed profile,the formation of congestion and subsequent dispersion are simulated.The results show that the model can reproduce the stochastic 2D trajectories of the vehicle and the marginal distribution of lateral movement.The proposed model can be used in both a simulation platform and a behavioral analysis of a human driver in traffic flow.
基金supported by the National Natural Science Foundation of China(No.12172167)。
文摘Nonlinearity and randomness are both the essential attributes for the real world,and the case is the same for the models of infectious diseases,for which the deterministic models can not give a complete picture of the evolution.However,although there has been a lot of work on stochastic epidemic models,most of them focus mainly on qualitative properties,which makes us somewhat ignore the original meaning of the parameter value.In this paper we extend the classic susceptible-infectious-removed(SIR)epidemic model by adding a white noise excitation and then we utilize the large deviation theory to quantitatively study the long-term coexistence exit problem with epidemic.Finally,in order to extend the meaning of parameters in the corresponding deterministic system,we tentatively introduce two new thresholds which then prove rational.
基金supported by Comunidad de Madrid within the framework of the Multiannual Agreement with Universidad Politécnica de Madrid to encourage research by young doctors(PRINCE).
文摘Cyber-Physical Systems are very vulnerable to sparse sensor attacks.But current protection mechanisms employ linear and deterministic models which cannot detect attacks precisely.Therefore,in this paper,we propose a new non-linear generalized model to describe Cyber-Physical Systems.This model includes unknown multivariable discrete and continuous-time functions and different multiplicative noises to represent the evolution of physical processes and randomeffects in the physical and computationalworlds.Besides,the digitalization stage in hardware devices is represented too.Attackers and most critical sparse sensor attacks are described through a stochastic process.The reconstruction and protectionmechanisms are based on aweighted stochasticmodel.Error probability in data samples is estimated through different indicators commonly employed in non-linear dynamics(such as the Fourier transform,first-return maps,or the probability density function).A decision algorithm calculates the final reconstructed value considering the previous error probability.An experimental validation based on simulation tools and real deployments is also carried out.Both,the new technology performance and scalability are studied.Results prove that the proposed solution protects Cyber-Physical Systems against up to 92%of attacks and perturbations,with a computational delay below 2.5 s.The proposed model shows a linear complexity,as recursive or iterative structures are not employed,just algebraic and probabilistic functions.In conclusion,the new model and reconstructionmechanism can protect successfully Cyber-Physical Systems against sparse sensor attacks,even in dense or pervasive deployments and scenarios.
文摘The asymptotic stability of two species stochastic Lotka-Volterra model is explored in this paper. Firstly, the Lotka-Volterra model with random parameter is built and reduced into the equivalent deterministic system by orthogonal polynomial approximation. Then, the linear stability theory and Routh-Hurwitz criterion for nonlinear deterministic systems are applied to the equivalent one. At last, at the aid of Lyapunov second method, we obtain that as the random intensity or statistical parameter of random variable is changed, the stability about stochastic Lotka-Volterra model is different from the deterministic system.
文摘This paper presents a study on a new rumor propagation model with nonlinear propagation rate and secondary propagation rate. We divide the total population into three groups, the ignorant, the spreader and the aware. The nonlinear incidence rate describes the psychological impact of certain serious rumors on social groups when the number of individuals spreading rumors becomes larger. The main contributions of this work are the development of a new rumor propagation model and some results of deterministic and stochastic analysis of the rumor propagation model. The results show the influence of nonlinear propagation rate and stochastic fluctuation on the dynamic behavior of the rumor propagation model by using Lyapunov function method and stochastic related knowledge. Numerical examples and simulation results are given to illustrate the results obtained.