High-precision chaotic radial basis function neural network model:Data forecasting for the Earth electromagnetic signal before a strong earthquake

The Earth’s natural pulse electromagnetic field data consists typically of an underlying variation tendency of intensity and irregularities.The change tendency may be related to the occurrence of earthquake disasters.Forecasting of the underlying intensity trend plays an important role in the analysis of data and disaster monitoring.Combining chaos theory and the radial basis function neural network,this paper proposes a forecasting model of the chaotic radial basis function neural network to conduct underlying intensity trend forecasting by the Earth’s natural pulse electromagnetic field signal.The main strategy of this forecasting model is to obtain parameters as the basis for optimizing the radial basis function neural network and to forecast the reconstructed Earth’s natural pulse electromagnetic field data.In verification experiments,we employ the 3 and 6 days’data of two channels as training samples to forecast the 14 and 21-day Earth’s natural pulse electromagnetic field data respectively.According to the forecasting results and absolute error results,the chaotic radial basis function forecasting model can fit the fluctuation trend of the actual signal strength,effectively reduce the forecasting error compared with the traditional radial basis function model.Hence,this network may be useful for studying the characteristics of the Earth’s natural pulse electromagnetic field signal before a strong earthquake and we hope it can contribute to the electromagnetic anomaly monitoring before the earthquake.

Interval grey number sequence prediction by using non-homogenous exponential discrete grey forecasting model

This paper aims to study a new grey prediction approach and its solution for forecasting the main system variable whose accurate value could not be collected while the potential value set could be defined. Based on the traditional nonhomogenous discrete grey forecasting model（NDGM）, the interval grey number and its algebra operations are redefined and combined with the NDGM model to construct a new interval grey number sequence prediction approach. The solving principle of the model is analyzed, the new accuracy evaluation indices, i.e. mean absolute percentage error of mean value sequence（MAPEM） and mean percent of interval sequence simulating value set covered（MPSVSC）, are defined and, the procedure of the interval grey number sequence based the NDGM（IG-NDGM） is given out. Finally, a numerical case is used to test the modelling accuracy of the proposed model. Results show that the proposed approach could solve the interval grey number sequence prediction problem and it is much better than the traditional DGM（1,1） model and GM（1,1） model.

Finite Element Implementation of the Exponential Drucker-Prager Plasticity Model for Adhesive Joints

This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influence of hydrostatic pressure on adhesive strength was investigated by a modified Arcan fixture designed particularly to induce a different state of hydrostatic pressure within an adhesive layer.The developed user subroutine UMAT,which utilizes an associated plastic flow during a plastic deformation,can provide a good agreement between the simulations and the experimental data.Better numerical stability at highly positive hydrostatic pressure loads for a very high order of exponential function can also be achieved compared to when a non-associated flow is used.

Short-term forecast in the early stage of the COVID-19 outbreak in Italy. Application of a weighted and cumulative average daily growth rate to an exponential decay model

To estimate the size of the novel coronavirus(COVID-19)outbreak in the early stage in Italy,this paper introduces the cumulated and weighted average daily growth rate(WR)to evaluate an epidemic curve.On the basis of an exponential decay model(EDM),we provide estimations of the WR in four-time intervals from February 27 to April 07,2020.By calibrating the parameters of the EDM to the reported data in Hubei Province of China,we also attempt to forecast the evolution of the outbreak.We compare the EDM applied to WR and the Gompertz model,which is based on exponential decay and is often used to estimate cumulative events.Specifically,we assess the performance of each model to short-term forecast of the epidemic,and to predict the final epidemic size.Based on the official counts for confirmed cases,the model applied to data from February 27 until the 17th of March estimate that the cumulative number of infected in Italy could reach 131,280(with a credibility interval 71,415-263,501)by April 25(credibility interval April 12 to May 3).With the data available until the 24st of March the peak date should be reached on May 3(April 23 to May 23)with 197,179 cumulative infections expected(130,033e315,269);with data available until the 31st of March the peak should be reached on May 4(April 25 to May 18)with 202,210 cumulative infections expected(155.235 e270,737);with data available until the 07st of April the peak should be reached on May 3(April 26 toMay 11)with 191,586(160,861-232,023)cumulative infections expected.Based on the average mean absolute percentage error(MAPE),cumulated infections forecasts provided by the EDM applied to WR performed better across all scenarios than the Gompertz model.An exponential decay model applied to the cumulated and weighted average daily growth rate appears to be useful in estimating the number of cases and peak of the COVID-19 outbreak in Italy and the model was more reliable in the exponential growth phase.

A Medium/Long-Range Forecast of Pacific Subtropical High Based on Dynamic Statistic Model Reconstruction

Based on the 500-hPa geopotential height field series of T106 numerical forecast products, by empirical orthogonal function （EOF） time-space separation, and on the hypotheses of EOF space-models being stable, the EOF time coefficient series were taken as dynamical statistic model variables. The dynamic system reconstruction idea and genetic algorithm were introduced to make the dynamical model parameters optimized, and a nonlinear dynamic statistic model of EOF separating time coefficient series was established. By the model time integral and EOF time-space reconstruction, a medium/long-range forecast of subtropical high was carried out. The results show that the dynamical model forecast and T106 numerical forecast were approximately similar in the short-range forecast （≤5 days）, but in the medium/long-range forecast （≥5 days）, the forecast results of dynamical model was superior to that of T106 numerical products. A new method and idea were presented for diagnosing and forecasting complicated weathers such as subtropical high, and showed a better application outlook.

Modelling HIV/AIDS Cases in Zambia: A Comparative Study of the Impact of Mandatory HIV Testing

In this study, a time series modeling approach is used to determine an ARIMA model and advance counterfactual forecasting at a point of policy intervention. We consider monthly data of HIV/AIDS cases from the Ministry of Health (Copperbelt province) of Zambia, for the period 2010 to 2019 and have a total of 120 observations. Results indicate that ARIMA (1, 0, 0) is an adequate model which best fits the HIV/AIDS time series data and is, therefore, suitable for forecasting cases. The model predicts a reduction from an average of 3500 to 3177 representing 14.29% in HIV/AIDS cases from 2017 (year of policy activation) to 2019, but the actual recorded cases dropped from 3500 to 1514 accounting for 57.4% in the same time frame.

Optimal Quota-Share and Excess-of-Loss Reinsurance and Investment with Heston’s Stochastic Volatility Model

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.

AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.

Predator Population Dynamics Involving Exponential Integral Function When Prey Follows Gompertz Model

The current study investigates the predator-prey problem with assumptions that interaction of predation has a little or no effect on prey population growth and the prey’s grow rate is time dependent. The prey is assumed to follow the Gompertz growth model and the respective predator growth function is constructed by solving ordinary differential equations. The results show that the predator population model is found to be a function of the well known exponential integral function. The solution is also given in Taylor’s series. Simulation study shows that the predator population size eventually converges either to a finite positive limit or zero or diverges to positive infinity. Under certain conditions, the predator population converges to the asymptotic limit of the prey model. More results are included in the paper.

A Model-based Phenomenological Investigation of Char Combustion Kinetics through Thermogravimetry

Five coal char samples were burnt in thermobalance with ramp heating rate of 30 K/min. The pore structure of these char samples was studied through mercury intrusion method. Combined with the kinetic theory of gases, the data of surface area was used in fitting the results. As a result, the kinetic triplet was given. The analysis showed that five char samples share almost the same intrinsic activation energy of the overall reaction. The phenomenological implication of the derived combustion rate equation was given.

Improved grey prediction model based on exponential grey action quantity

With the passage of time, it has become important to investigate new methods for updating data to better fit the trends of the grey prediction model. The traditional GM（1,1） usually sets the grey action quantity as a constant. Therefore, it cannot effectively fit the dynamic characteristics of the sequence, which results in the grey model having a low precision. The linear grey action quantity model cannot represent the index change law. This paper presents a grey action quantity model, the exponential optimization grey model（EOGM（1,1））, based on the exponential type of grey action quantity; it is constructed based on the exponential characteristics of the grey prediction model. The model can fully reflect the exponential characteristics of the simulation series with time. The exponential sequence has a higher fitting accuracy. The optimized result is verified using a numerical example for the fluctuating sequence and a case study for the index of the tertiary industry＇s GDP. The results show that the model improves the precision of the grey forecasting model and reduces the prediction error.

Development of 3D top coal caving angle model for fully mechanized extra-thick coal seam mining

During high-intensity,fully mechanized mining of extra-thick coal seam,the top coal would cave to a certain 3D form.Based on the data collected during drilling,a 3D model of top coal caving surface space was established to determine the relationship between the location of the stope roof and the caving surface,enabling the mathematical computation of the top caving angle(φ).The drilling method was employed to measure the top caving angle on two extra-thick fully mechanized coal caving faces under the conditions of three geological structures,namely,no geological structure,igneous rock structure,and fault structure.The results show that the value of top caving angle could be accurately estimated on-site with the 9-parameter 3D top coal caving surface model built with the drilling method.This method is a novel on-site measurement that can be easily applied.Our findings reveal that the characteristics of the coal-rock in the two mining faces are different;yet their caving angles follow the ruleφ_(igneous rock structure)<φ_(no geological structure)<φ_(fault structure).Finally,through the data fitting with two indexes(the top coal uniaxial compressive strength and the top caving angle),it is found that the relationship between the two indexes satisfies an exponential decay function.

MODELLING OF TENSILE STRESS-STRAIN CURVE OF WOVEN FABRICS

A general shape of tensile stress-strain curves of woven fabrics is first recognised by puttingtested and predicted results together.An exponential function with two parameters is then selectedfor the prediction of tensile stress-strain relationship.The predicted results by using the proposedfunction show excellent agreement with experimental data.

Order Exponential Evaluation Model for Road Traffic Safety in City Clusters

This study presents an order exponential model for estimating road traffic safety in city clusters.The proposed model introduces the traffic flow intrinsic properties and uses the characteristics and regular patterns of traffic development to identify road traffic safety levels in city clusters.Additionally,an evaluation index system of city cluster road traffic safety was constructed based on the spatial and temporal distribution.Then Order Exponential Evaluation Model(OEEM),a comprehensive model using order exponent function for road traffic safety evaluation,was put forward,which considers the main characteristics and the generation process of traffic accidents.The model effectively controlled the unsafe behavior of the traffic system.It could define the levels of city cluster road traffic safety and dynamically detect road safety risk.The proposed model was verified with statistical data from three Chinese city clusters by comparing the common model for road traffic safety with an ideal model.The results indicate that the order exponent approach undertaken in this study can be extended and applied to other research topics and fields.

Dynamic Spatio-Temporal Modeling in Disease Mapping

Spatio-temporal models are valuable tools for disease mapping and understanding the geographical distribution of diseases and temporal dynamics. Spatio-temporal models have been proven empirically to be very complex and this complexity has led many to oversimply and model the spatial and temporal dependencies independently. Unlike common practice, this study formulated a new spatio-temporal model in a Bayesian hierarchical framework that accounts for spatial and temporal dependencies jointly. The spatial and temporal dependencies were dynamically modelled via the matern exponential covariance function. The temporal aspect was captured by the parameters of the exponential with a first-order autoregressive structure. Inferences about the parameters were obtained via Markov Chain Monte Carlo (MCMC) techniques and the spatio-temporal maps were obtained by mapping stable posterior means from the specific location and time from the best model that includes the significant risk factors. The model formulated was fitted to both simulation data and Kenya meningitis incidence data from 2013 to 2019 along with two covariates;Gross County Product (GCP) and average rainfall. The study found that both average rainfall and GCP had a significant positive association with meningitis occurrence. Also, regarding geographical distribution, the spatio-temporal maps showed that meningitis is not evenly distributed across the country as some counties reported a high number of cases compared with other counties.

高压管汇材料疲劳性能测试及P-S-N模型曲线的拟合

高压管汇作为压裂设备中的主要易损件之一,其失效危害较大。它的失效原因主要是疲劳、冲蚀、腐蚀或者材料缺陷引起的刺漏和爆裂,其中尤以疲劳失效最不可预估。目前,对于高压管汇材料的疲劳性能研究不够深入,为解决高压管汇材料疲劳寿命的准确描述问题,以某国产高压管汇材料为例,进行了一系列疲劳试验,并基于试验数据,采用多种分布模型和不同S-N模型进行拟合分析,得出综合评价拟合能力最强的P-S-N模型。结果表明,该材料在中长疲劳寿命区,Weibull三参数模型在7级应力水平下综合评价能力最好;在存活率分别为50%、90%、99%、99.9%时,指数S-N模型的拟合系数均大于0.98,拟合能力最好。得出的P-S-N模型曲线可以为高压管汇的疲劳寿命以及安全设计提供依据。

基于水电储能调节的风光水发电联合优化调度策略

为缓解新能源装机容量扩大引起的弃风弃光现象,在已有梯级水电上下电站之间加入储能泵站,提出风光水储短期优化调度策略。构建以风光水储系统负荷跟踪误差最小、梯级水电站发电量最大和梯级水电站发电耗水量最小的多目标优化调度模型;提出基于季节性自回归移动平均(seasonal auto-regressive lntegrated moving average, SARIMA)模型和Copula函数的风光出力预测模型作为优化调度模型的边界条件,通过SARIMA预测模型将风光出力历史数据分解为季节性分量、趋势分量以及随机噪声余项进行全天96个调度时段风光出力预测,并叠加上基于Copula函数生成风光出力预测误差,然后通过拉丁超立方采样以及K-means聚类进行场景生成和缩减得到5个风光出力场景。选取风光典型日出力数据为例进行算例分析,算例结果表明:所提预测模型较SARIMA模型可以显著提高预测准确度,模型预测风光出力均方根误差从33.34、229.49 MW分别下降至0.697、9.534 MW;所提优化调度策略可以在全年丰、平、枯水期有效减少弃风弃光现象,并可将过剩新能源中的50%转化为上级水库储存水能。

变形监测灰色预测模型对比及替代方法研究

将变形监测灰色预测模型分为传统GM(1,1)模型及其改进模型、非齐次灰色模型、GM(1,1)幂模型及其改进模型3种类型,以Origin拟合函数Exp2PModl、Exponential和SRichards2作为3类灰色预测模型的替代方法,基于理论研究和实例验证对比分析3类灰色预测模型及其替代方法。结果表明,3类灰色预测模型在拟合函数、有无极限值、适合等时距或非等时距建模和适用范围等方面存在显著差异,需要根据变形监测数据特征选择合适的灰色预测模型类别;与3类灰色预测模型相比,Origin拟合函数在参数求解和建模数据要求上更具优势,而且可以得到相当甚至更高的拟合或预测精度,除需要编程实现的特殊优化目标外,完全可以代替灰色预测模型用于变形监测。

基于改进扩散模型的温度预报

针对传统数值预报模式计算时间长和计算资源消耗大的问题,以及现有深度学习预报方法在温度预报结果上不精确,且预测结果模糊的问题,提出了一个新的温度预报模型。首先,设计了一个时空信息捕捉模块,将该模块捕获的长期依赖信息,作为扩散模型的生成条件,赋予扩散模型预报的能力;其次,设计了一个新的平衡损失函数,同时保护了扩散模型的生成能力和时空信息捕捉模块对时空信息的捕捉能力;最后,基于美国国家环境预报中心的再分析数据进行预报,与现有的深度学习方法相比,所提模型预报结果的质量在均方误差(mean square error,MSE)上降低了17.3%,在均方根误差(root mean square error,RMSE)上降低了9.14%,在峰值信噪比(peak signal to noise ratio,PSNR)上提升了5.1%。改进的扩散模型能有效地捕捉时空依赖的关系,有效地进行时空序列预测,效果优于其他对比方法。