In the assessment of car insurance claims,the claim rate for car insurance presents a highly skewed probability distribution,which is typically modeled using Tweedie distribution.The traditional approach to obtaining ...In the assessment of car insurance claims,the claim rate for car insurance presents a highly skewed probability distribution,which is typically modeled using Tweedie distribution.The traditional approach to obtaining the Tweedie regression model involves training on a centralized dataset,when the data is provided by multiple parties,training a privacy-preserving Tweedie regression model without exchanging raw data becomes a challenge.To address this issue,this study introduces a novel vertical federated learning-based Tweedie regression algorithm for multi-party auto insurance rate setting in data silos.The algorithm can keep sensitive data locally and uses privacy-preserving techniques to achieve intersection operations between the two parties holding the data.After determining which entities are shared,the participants train the model locally using the shared entity data to obtain the local generalized linear model intermediate parameters.The homomorphic encryption algorithms are introduced to interact with and update the model intermediate parameters to collaboratively complete the joint training of the car insurance rate-setting model.Performance tests on two publicly available datasets show that the proposed federated Tweedie regression algorithm can effectively generate Tweedie regression models that leverage the value of data fromboth partieswithout exchanging data.The assessment results of the scheme approach those of the Tweedie regressionmodel learned fromcentralized data,and outperformthe Tweedie regressionmodel learned independently by a single party.展开更多
Expanding photovoltaic(PV)resources in rural-grid areas is an essential means to augment the share of solar energy in the energy landscape,aligning with the“carbon peaking and carbon neutrality”objectives.However,ru...Expanding photovoltaic(PV)resources in rural-grid areas is an essential means to augment the share of solar energy in the energy landscape,aligning with the“carbon peaking and carbon neutrality”objectives.However,rural power grids often lack digitalization;thus,the load distribution within these areas is not fully known.This hinders the calculation of the available PV capacity and deduction of node voltages.This study proposes a load-distribution modeling approach based on remote-sensing image recognition in pursuit of a scientific framework for developing distributed PV resources in rural grid areas.First,houses in remote-sensing images are accurately recognized using deep-learning techniques based on the YOLOv5 model.The distribution of the houses is then used to estimate the load distribution in the grid area.Next,equally spaced and clustered distribution models are used to adaptively determine the location of the nodes and load power in the distribution lines.Finally,by calculating the connectivity matrix of the nodes,a minimum spanning tree is extracted,the topology of the network is constructed,and the node parameters of the load-distribution model are calculated.The proposed scheme is implemented in a software package and its efficacy is demonstrated by analyzing typical remote-sensing images of rural grid areas.The results underscore the ability of the proposed approach to effectively discern the distribution-line structure and compute the node parameters,thereby offering vital support for determining PV access capability.展开更多
In this paper, platoons of autonomous vehicles operating in urban road networks are considered. From a methodological point of view, the problem of interest consists of formally characterizing vehicle state trajectory...In this paper, platoons of autonomous vehicles operating in urban road networks are considered. From a methodological point of view, the problem of interest consists of formally characterizing vehicle state trajectory tubes by means of routing decisions complying with traffic congestion criteria. To this end, a novel distributed control architecture is conceived by taking advantage of two methodologies: deep reinforcement learning and model predictive control. On one hand, the routing decisions are obtained by using a distributed reinforcement learning algorithm that exploits available traffic data at each road junction. On the other hand, a bank of model predictive controllers is in charge of computing the more adequate control action for each involved vehicle. Such tasks are here combined into a single framework:the deep reinforcement learning output(action) is translated into a set-point to be tracked by the model predictive controller;conversely, the current vehicle position, resulting from the application of the control move, is exploited by the deep reinforcement learning unit for improving its reliability. The main novelty of the proposed solution lies in its hybrid nature: on one hand it fully exploits deep reinforcement learning capabilities for decisionmaking purposes;on the other hand, time-varying hard constraints are always satisfied during the dynamical platoon evolution imposed by the computed routing decisions. To efficiently evaluate the performance of the proposed control architecture, a co-design procedure, involving the SUMO and MATLAB platforms, is implemented so that complex operating environments can be used, and the information coming from road maps(links,junctions, obstacles, semaphores, etc.) and vehicle state trajectories can be shared and exchanged. Finally by considering as operating scenario a real entire city block and a platoon of eleven vehicles described by double-integrator models, several simulations have been performed with the aim to put in light the main f eatures of the proposed approach. Moreover, it is important to underline that in different operating scenarios the proposed reinforcement learning scheme is capable of significantly reducing traffic congestion phenomena when compared with well-reputed competitors.展开更多
The generation and propagation mechanism of strong nonlinear waves in the South China Sea is an essential research area. In this study, the third-generation wave model WAVEWATCH Ⅲ is employed to simulate wave fields ...The generation and propagation mechanism of strong nonlinear waves in the South China Sea is an essential research area. In this study, the third-generation wave model WAVEWATCH Ⅲ is employed to simulate wave fields under extreme sea states. The model, integrating the ST6 source term, is validated against observed data, demonstrating its credibility. The spatial distribution of the occurrence probability of strong nonlinear waves during typhoons is shown, and the waves in the straits and the northeastern part of the South China Sea show strong nonlinear characteristics. The high-order spectral model HOS-ocean is employed to simulate the random wave surface series beneath five different platform areas. The waves during the typhoon exhibit strong nonlinear characteristics, and freak waves exist. The space-varying probability model is established to describe the short-term probability distribution of nonlinear wave series. The exceedance probability distributions of the wave surface beneath different platform areas are compared and analyzed. The results show that with an increase in the platform area, the probability of a strong nonlinear wave beneath the platform increases.展开更多
Crude oil scheduling optimization is an effective method to enhance the economic benefits of oil refining.But uncertainties,including uncertain demands of crude distillation units(CDUs),might make the production plans...Crude oil scheduling optimization is an effective method to enhance the economic benefits of oil refining.But uncertainties,including uncertain demands of crude distillation units(CDUs),might make the production plans made by the traditional deterministic optimization models infeasible.A data-driven Wasserstein distributionally robust chance-constrained(WDRCC)optimization approach is proposed in this paper to deal with demand uncertainty in crude oil scheduling.First,a new deterministic crude oil scheduling optimization model is developed as the basis of this approach.The Wasserstein distance is then used to build ambiguity sets from historical data to describe the possible realizations of probability distributions of uncertain demands.A cross-validation method is advanced to choose suitable radii for these ambiguity sets.The deterministic model is reformulated as a WDRCC optimization model for crude oil scheduling to guarantee the demand constraints hold with a desired high probability even in the worst situation in ambiguity sets.The proposed WDRCC model is transferred into an equivalent conditional value-at-risk representation and further derived as a mixed-integer nonlinear programming counterpart.Industrial case studies from a real-world refinery are conducted to show the effectiveness of the proposed method.Out-of-sample tests demonstrate that the solution of the WDRCC model is more robust than those of the deterministic model and the chance-constrained model.展开更多
This paper considers setting different dips for different sub-faults to fit the actual rupture situation based on the fault rupture of the 2013 Lushan M_(S)7.0 earthquake.Meanwhile,combined with the coseismic GNSS dat...This paper considers setting different dips for different sub-faults to fit the actual rupture situation based on the fault rupture of the 2013 Lushan M_(S)7.0 earthquake.Meanwhile,combined with the coseismic GNSS data of the Lushan earthquake,the source parameters and sliding distribution of the Lushan earthquake fault are inversed.Firstly,we use the gradient based optimizer(GBO)in nonlinear inversion to obtain the source parameters of this seismic fault.The inversion results indicate that the strike of the fault is 206.52°,the dip is 44.10°,the length is 21.92 km,and the depth is 12.79 km.To refine the sliding distribution of the seismic fault,the seismic fault is divided into 3×3 sub-faults.Then,we fix the central sub-fault dip of 44.10°;the dip of other sub-faults is obtained by iteration.After that,the model is further divided into a fault layer model composed of 23×19 sub fault slices,and using the Matlab fitting function is used to fit the dip of the 23×19 sub faults.Finally,the Lushan seismic fault plane is established as a shovel structure with steep upper and gentle lower,steep south and gentle north.The slip distribution inversion results indicate that the depth of the slip peak is 13 km,the corresponding maximum slip momentum is 0.67 m,the seismic moment is 1.10×10^(19)N·m and the corresponding moment magnitude is MW6.66.The results above are consistent with the research results of seismology.展开更多
Neutron-skin thickness is a key parameter for a neutron-rich nucleus;however,it is difficult to determine.In the framework of the Lanzhou Quantum Molecular Dynamics(LQMD)model,a possible probe for the neutron-skin thi...Neutron-skin thickness is a key parameter for a neutron-rich nucleus;however,it is difficult to determine.In the framework of the Lanzhou Quantum Molecular Dynamics(LQMD)model,a possible probe for the neutron-skin thickness(δ_(np))of neutron-rich ^(48)Ca was studied in the 140A MeV ^(48)Ca+^(9)Be projectile fragmentation reaction based on the parallel momentum distribution(p∥)of the residual fragments.A Fermi-type density distribution was employed to initiate the neutron density distributions in the LQMD simulations.A combined Gaussian function with different width parameters for the left side(Γ_(L))and the right side(Γ_(R))in the distribution was used to describe the p∥of the residual fragments.Taking neutron-rich sulfur isotopes as examples,Γ_(L) shows a sensitive correlation withδ_(np) of ^(48)Ca,and is proposed as a probe for determining the neutron skin thickness of the projectile nucleus.展开更多
Snow interacting with a high-speed train can cause the formation of ice in the train bogie region and affect its safety.In this study,a wind-snow multiphase numerical approach is introduced for high-speed train bogies...Snow interacting with a high-speed train can cause the formation of ice in the train bogie region and affect its safety.In this study,a wind-snow multiphase numerical approach is introduced for high-speed train bogies on the basis of the Euler-Lagrange discrete phase model.A particle-wall impact criterion is implemented to account for the presence of snow particles on the surface.Subsequently,numerical simulations are conducted,considering various snow particle diameter distributions and densities.The research results indicate that when the particle diameter is relatively small,the distribution of snow particles in the bogie cavity is relatively uniform.However,as the particle diameter increases,the snow particles in the bogie cavity are mainly located in the rear wheel pairs of the bogie.When the more realistic Rosin-Rammler diameter distribution is applied to snow particles,the positions of snow particles with different diameters vary in the bogie cavity.More precisely,smaller diameter particles are primarily located in the front and upper parts of the bogie cavity,while larger diameter snow particles accumulate at the rear and in the lower parts of the bogie cavity.展开更多
In recent years,Meloidogyne enterolobii has emerged as a major parasitic nematode infesting many plants in tropical or subtropical areas.However,the regions of potential distribution and the main contributing environm...In recent years,Meloidogyne enterolobii has emerged as a major parasitic nematode infesting many plants in tropical or subtropical areas.However,the regions of potential distribution and the main contributing environmental variables for this nematode are unclear.Under the current climate scenario,we predicted the potential geographic distributions of M.enterolobii worldwide and in China using a Maximum Entropy(MaxEnt)model with the occurrence data of this species.Furthermore,the potential distributions of M.enterolobii were projected under three future climate scenarios(BCC-CSM2-MR,CanESM5 and CNRM-CM6-1)for the periods 2050s and 2090s.Changes in the potential distribution were also predicted under different climate conditions.The results showed that highly suitable regions for M.enterolobii were concentrated in Africa,South America,Asia,and North America between latitudes 30°S to 30°N.Bio16(precipitation of the wettest quarter),bio10(mean temperature of the warmest quarter),and bio11(mean temperature of the coldest quarter)were the variables contributing most in predicting potential distributions of M.enterolobii.In addition,the potential suitable areas for M.enterolobii will shift toward higher latitudes under future climate scenarios.This study provides a theoretical basis for controlling and managing this nematode.展开更多
Saharan dust represents more than 50%of the total desert dust emitted around the globe and its radiative effect significantly affects the atmospheric circulation at a continental scale.Previous studies on dust vertica...Saharan dust represents more than 50%of the total desert dust emitted around the globe and its radiative effect significantly affects the atmospheric circulation at a continental scale.Previous studies on dust vertical distribution and the Saharan Air Layer(SAL)showed some shortcomings that could be attributed to imperfect representation of the effects of deep convection and scavenging.The authors investigate here the role of deep convective transport and scavenging on the vertical distribution of mineral dust over Western Africa.Using multi-year(2006-2010)simulations performed with the variable-resolution(zoomed)version of the LMDZ climate model.Simulations are compared with aerosol amounts recorded by the Aerosol Robotic Network(AERONET)and with vertical profiles of the Cloud-Aerosol Lidar with Orthogonal Polarization(CALIOP)measurements.LMDZ allows a thorough examination of the respective roles of deep convective transport,convective and stratiform scavenging,boundary layer transport,and advection processes on the vertical mineral dust distribution over Western Africa.The comparison of simulated dust Aerosol Optical Depth(AOD)and distribution with measurements suggest that scavenging in deep convection and subsequent re-evaporation of dusty rainfall in the lower troposphere are critical processes for explaining the vertical distribution of desert dust.These processes play a key role in maintaining a well-defined dust layer with a sharp transition at the top of the SAL and in establishing the seasonal cycle of dust distribution.This vertical distribution is further reshaped offshore in the Inter-Tropical Convergence Zone(ITCZ)over the Atlantic Ocean by marine boundary layer turbulent and convective transport and wet deposition at the surface.展开更多
This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interf...This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (T<sub>react</sub>). In addition, to validate the proposed heat transfer model, T<sub>react</sub> obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O<sub>2</sub> and H<sub>2</sub>O from the inlet to the outlet is more even with the increase in T<sub>ini</sub> due to the lower performance of O<sub>2</sub> reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in T<sub>ini</sub> and the value of current density is smaller with the increase in T<sub>ini </sub>due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in T<sub>react</sub> from the inlet to the outlet is more even with the increase in T<sub>ini</sub> irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in T<sub>ini </sub>due to the lower performance of O<sub>2</sub> reduction reaction.展开更多
Abiotic factors play an important role in species localisation,but biotic and anthropogenic predictors must also be considered in distribution modelling for models to be biologically meaningful.In this study,we formal...Abiotic factors play an important role in species localisation,but biotic and anthropogenic predictors must also be considered in distribution modelling for models to be biologically meaningful.In this study,we formalised the biotic predictors of nesting sites for four threatened Caucasian vultures by including species distribution models(wild ungulates,nesting tree species)as biotic layers in the vulture Maxent models.Maxent was applied in the R dismo package and the best set of the model parameters were defined in the R ENMeval package.Performance metrics were continuous Boyce index,Akaike's information criterion,the area under receiver operating curve and true skill statistics.We also calculated and evaluated the null models.Kernel density estimation method was applied to assess the overlap of vulture ecological niches in the environmental space.The accessibility of anthropogenic food resources was estimated using the Path Distance measure that considers elevation gradient.The availability of pine forests(Scots Pine)and wild ungulates(Alpine Chamois and Caucasian Goat)contributed the most(29.6%and 34.3%)to Cinereous Vulture(Aegypius monachus)nesting site model.Wild ungulate distribution also contributed significantly(about 46%)to the Bearded Vulture(Gypaetus barbatus)model.This scavenger nests in the highlands of the Caucasus at a minimum distance of 5–10 km from anthropogenic facilities.In contrast,livestock as a food source was most important in colony distribution of Griffon Vulture(Gyps fulvus).The contribution of distances to settlements and agricultural facilities to the model was 45%.The optimal distance from Egyptian Vulture(Neophron percnopterus)nesting sites to settlements was only 3–10 km,to livestock facilities no more than 15 km with the factor contribution of about 57%.Excluding the wild ungulate availability,the ecological niches of studied vultures overlapped significantly.Despite similar foraging and nesting requirements,Caucasian vultures are not pronounced nesting and trophic competitors due to the abundance of nesting sites,anthropogenic food sources and successful niche sharing.展开更多
Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to pr...Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to predict the stiffness and shear distribution of glulam beams in elastic-plastic stage,and consequently,the failure of such glulam beams cannot be predicted effectively.To address these issues,an analytical method considering material nonlinearity was proposed for glulam beams,and the calculating equations of deflection and shear stress distribution for different failure modes were established.The proposed method was verified by experiments and numerical models under the corresponding conditions.Results showed that the theoretical calculations were in good agreement with experimental and numerical results,indicating that the equations proposed in this paper were reliable and accurate for such glulam beams with wood material in the elastic-plastic stage ignoring the influence of mechanic properties in radial and tangential directions of wood.Furthermore,the experimental results reported by the previous studies indicated that the method was applicable and could be used as a theoretical reference for predicting the failure of glulam beams.展开更多
Micro sliding phenomenon widely exists in the operation process of mechanical systems,and the micro sliding friction mechanism is always a research hotspot.In this work,based on the total reflection method,a measuring...Micro sliding phenomenon widely exists in the operation process of mechanical systems,and the micro sliding friction mechanism is always a research hotspot.In this work,based on the total reflection method,a measuring device for interface contact behavior under two-dimensional(2D)vibration is built.The stress distribution is characterized by the light intensity distribution of the contact image,and the interface contact behavior in the 2D vibration process is studied.It is found that the vibration angle of the normal direction of the contact surface and its fluctuation affect the interface friction coefficient,the tangential stiffness,and the fluctuation amplitude of the stress distribution.Then they will affect the change of friction state and energy dissipation in the process of micro sliding.Further,an improved micro sliding friction model is proposed based on the experimental analysis,with the nonlinear change of contact parameters caused by the normal contact stress distribution fluctuation taken into account.This model considers the interface tangential stiffness fluctuation,friction coefficient hysteresis,and stress distribution fluctuation,whose simulation results are consistent well with the experimental results.It is found that considering the nonlinear effect of a certain contact parameter alone may bring a greater error to the prediction of friction behavior.Only by integrating multiple contact parameters can the accuracy of friction prediction is improved.展开更多
A physics-based analytical expression that predicts the charge,electrical field and potential distributions along the gated region of the GaN HEMT channel has been developed.Unlike the gradual channel approximation(GC...A physics-based analytical expression that predicts the charge,electrical field and potential distributions along the gated region of the GaN HEMT channel has been developed.Unlike the gradual channel approximation(GCA),the proposed model considers the non-uniform variation of the concentration under the gated region as a function of terminal applied volt-ages.In addition,the model can capture the influence of mobility and channel temperature on the charge distribution trend.The comparison with the hydrodynamic(HD)numerical simulation showed a high agreement of the proposed model with numerical data for different bias conditions considering the self-heating and quantization of the electron concentration.The ana-lytical nature of the model allows us to reduce the computational and time cost of the simulation.Also,it can be used as a core expression to develop a complete physics-based transistorⅣmodel without GCA limitation.展开更多
Species distribution models have been widely used to explore suitable habitats of species,the impact of climate change on the distribution of suitable habitats of species,and the construction of ecological reserves.Th...Species distribution models have been widely used to explore suitable habitats of species,the impact of climate change on the distribution of suitable habitats of species,and the construction of ecological reserves.This paper introduced species distribution models commonly used in biodiversity analysis,as well as model performance evaluation indexes,challenges in the application of species distribution models,and finally prospected the development trend of research on species distribution models.展开更多
For studying the carbon thermal reduction rules of titanium in hot metal and providing a theoretical basis for the blast furnace(BF) hearth protection, the distribution behavior of titanium between low-titanium slag...For studying the carbon thermal reduction rules of titanium in hot metal and providing a theoretical basis for the blast furnace(BF) hearth protection, the distribution behavior of titanium between low-titanium slag system of CaO-SiO2-MgO-Al2O3-TiO2 and hot metal was studied using analytical reagents in a temperature range from 1350 °C to 1600 °C. Through high temperature melting, rapid quenching, chemical analysis and thermodynamic model calculating, the results showed that the increase of reaction temperature, which improved the titanium distribution L(Ti) and lowered the system activity coefficient γsys, leads to the rise of equilibrium constant. Combined with Wagner and congregated electron phase models, the data obtained in distribution experiments were used to fit out the Gibbs free energy formula of titanium carbothermic reduction. Finally, the relations between the contents of Si and Ti in hot metal and the titanium load to reach the minimum w(Ti) for the formation of Ti C were given.展开更多
For performance optimization such as placement,interconnect synthesis,and routing, an efficient and accurate interconnect delay metric is critical,even in design tools development like design for yield (DFY) and des...For performance optimization such as placement,interconnect synthesis,and routing, an efficient and accurate interconnect delay metric is critical,even in design tools development like design for yield (DFY) and design for manufacture (DFM). In the nanometer regime, the recently proposed delay models for RLC interconnects based on statistical probability density function (PDF)interpretation such as PRIMO,H-gamma,WED and RLD bridge the gap between accuracy and efficiency. However, these models always require table look-up when operating. In this paper, a novel delay model based on the Birnbaum-Saunders distribution (BSD) is presented. BSD can accomplish interconnect delay estimation fast and accurately without table look-up operations. Furthermore, it only needs the first two moments to match. Experimental results in 90nm technology show that BSD is robust, easy to implement,efficient,and accurate.展开更多
The classical source-to-trap petroleum system concept only considers the migration and accumulation of conventional oil and gas in traps driven dominantly by buoyance in a basin,although revised and improved,even some...The classical source-to-trap petroleum system concept only considers the migration and accumulation of conventional oil and gas in traps driven dominantly by buoyance in a basin,although revised and improved,even some new concepts as composite petroleum system,total petroleum system,total composite petroleum system,were proposed,but they do not account for the vast unconventional oil and gas reservoirs within the system,which is not formed and distributed in traps dominantly by buoyancedriven.Therefore,the petroleum system concept is no longer adequate in dealing with all the oil and gas accumulations in a basin where significant amount of the unconventional oil and gas resources are present in addition to the conventional oil and gas accumulations.This paper looked into and analyzed the distribution characteristics of conventional and unconventional oil/gas reservoirs and their differences and correlations in petroliferous basins in China and North America,and then proposed whole petroleum system(WPS)concept,the WPS is defined as a natural system that encompasses all the conventional and unconventional oil and gas,reservoirs and resources originated from organic matter in source rocks,the geological elements and processes involving the formation,evolution,and distribution of these oil and gas,reservoirs and resources.It is found in the WPS that there are three kinds of hydrocarbons dynamic fields,three kinds of original hydrocarbons,three kinds of reservoir rocks,and the coupling of these three essential elements lead to the basic ordered distribution model of shale oil/gas reservoirs contacting or interbeded with tight oil/gas reservoirs and separated conventional oil/gas reservoirs from source rocks upward,which is expressed as“S\T-C”.Abnormal conditions lead to other three special ordered distribution models:The first is that with shale oil/gas reservoirs separated from tight oil/gas reservoirs.The second is that with two direction ordered distributions from source upward and downward.The third is with lateral distribution from source outside.展开更多
Due to the fact that the emergency medicine distribution is vital to the quick response to urgent demand when an epidemic occurs, the optimal vaccine distribution approach is explored according to the epidemic diffusi...Due to the fact that the emergency medicine distribution is vital to the quick response to urgent demand when an epidemic occurs, the optimal vaccine distribution approach is explored according to the epidemic diffusion rule and different urgency degrees of affected areas with the background of the epidemic outbreak in a given region. First, the SIQR (susceptible, infected, quarantined,recovered) epidemic model with pulse vaccination is introduced to describe the epidemic diffusion rule and obtain the demanded vaccine in each pulse. Based on the SIQR model, the affected areas are clustered by using the self-organizing map (SOM) neutral network to qualify the results. Then, a dynamic vaccine distribution model is formulated, incorporating the results of clustering the affected areas with the goals of both reducing the transportation cost and decreasing the unsatisfied demand for the emergency logistics network. Numerical study with twenty affected areas and four distribution centers is carried out. The corresponding numerical results indicate that the proposed approach can make an outstanding contribution to controlling the affected areas with a relatively high degree of urgency, and the comparison results prove that the performance of the clustering method is superior to that of the non-clustering method on controlling epidemic diffusion.展开更多
基金This research was funded by the National Natural Science Foundation of China(No.62272124)the National Key Research and Development Program of China(No.2022YFB2701401)+3 种基金Guizhou Province Science and Technology Plan Project(Grant Nos.Qiankehe Paltform Talent[2020]5017)The Research Project of Guizhou University for Talent Introduction(No.[2020]61)the Cultivation Project of Guizhou University(No.[2019]56)the Open Fund of Key Laboratory of Advanced Manufacturing Technology,Ministry of Education(GZUAMT2021KF[01]).
文摘In the assessment of car insurance claims,the claim rate for car insurance presents a highly skewed probability distribution,which is typically modeled using Tweedie distribution.The traditional approach to obtaining the Tweedie regression model involves training on a centralized dataset,when the data is provided by multiple parties,training a privacy-preserving Tweedie regression model without exchanging raw data becomes a challenge.To address this issue,this study introduces a novel vertical federated learning-based Tweedie regression algorithm for multi-party auto insurance rate setting in data silos.The algorithm can keep sensitive data locally and uses privacy-preserving techniques to achieve intersection operations between the two parties holding the data.After determining which entities are shared,the participants train the model locally using the shared entity data to obtain the local generalized linear model intermediate parameters.The homomorphic encryption algorithms are introduced to interact with and update the model intermediate parameters to collaboratively complete the joint training of the car insurance rate-setting model.Performance tests on two publicly available datasets show that the proposed federated Tweedie regression algorithm can effectively generate Tweedie regression models that leverage the value of data fromboth partieswithout exchanging data.The assessment results of the scheme approach those of the Tweedie regressionmodel learned fromcentralized data,and outperformthe Tweedie regressionmodel learned independently by a single party.
基金supported by the State Grid Science&Technology Project of China(5400-202224153A-1-1-ZN).
文摘Expanding photovoltaic(PV)resources in rural-grid areas is an essential means to augment the share of solar energy in the energy landscape,aligning with the“carbon peaking and carbon neutrality”objectives.However,rural power grids often lack digitalization;thus,the load distribution within these areas is not fully known.This hinders the calculation of the available PV capacity and deduction of node voltages.This study proposes a load-distribution modeling approach based on remote-sensing image recognition in pursuit of a scientific framework for developing distributed PV resources in rural grid areas.First,houses in remote-sensing images are accurately recognized using deep-learning techniques based on the YOLOv5 model.The distribution of the houses is then used to estimate the load distribution in the grid area.Next,equally spaced and clustered distribution models are used to adaptively determine the location of the nodes and load power in the distribution lines.Finally,by calculating the connectivity matrix of the nodes,a minimum spanning tree is extracted,the topology of the network is constructed,and the node parameters of the load-distribution model are calculated.The proposed scheme is implemented in a software package and its efficacy is demonstrated by analyzing typical remote-sensing images of rural grid areas.The results underscore the ability of the proposed approach to effectively discern the distribution-line structure and compute the node parameters,thereby offering vital support for determining PV access capability.
文摘In this paper, platoons of autonomous vehicles operating in urban road networks are considered. From a methodological point of view, the problem of interest consists of formally characterizing vehicle state trajectory tubes by means of routing decisions complying with traffic congestion criteria. To this end, a novel distributed control architecture is conceived by taking advantage of two methodologies: deep reinforcement learning and model predictive control. On one hand, the routing decisions are obtained by using a distributed reinforcement learning algorithm that exploits available traffic data at each road junction. On the other hand, a bank of model predictive controllers is in charge of computing the more adequate control action for each involved vehicle. Such tasks are here combined into a single framework:the deep reinforcement learning output(action) is translated into a set-point to be tracked by the model predictive controller;conversely, the current vehicle position, resulting from the application of the control move, is exploited by the deep reinforcement learning unit for improving its reliability. The main novelty of the proposed solution lies in its hybrid nature: on one hand it fully exploits deep reinforcement learning capabilities for decisionmaking purposes;on the other hand, time-varying hard constraints are always satisfied during the dynamical platoon evolution imposed by the computed routing decisions. To efficiently evaluate the performance of the proposed control architecture, a co-design procedure, involving the SUMO and MATLAB platforms, is implemented so that complex operating environments can be used, and the information coming from road maps(links,junctions, obstacles, semaphores, etc.) and vehicle state trajectories can be shared and exchanged. Finally by considering as operating scenario a real entire city block and a platoon of eleven vehicles described by double-integrator models, several simulations have been performed with the aim to put in light the main f eatures of the proposed approach. Moreover, it is important to underline that in different operating scenarios the proposed reinforcement learning scheme is capable of significantly reducing traffic congestion phenomena when compared with well-reputed competitors.
基金financially supported by the National Key R&D Program of China(No.2022YFC3104205)the National Natural Science Foundation of China(No.42377457).
文摘The generation and propagation mechanism of strong nonlinear waves in the South China Sea is an essential research area. In this study, the third-generation wave model WAVEWATCH Ⅲ is employed to simulate wave fields under extreme sea states. The model, integrating the ST6 source term, is validated against observed data, demonstrating its credibility. The spatial distribution of the occurrence probability of strong nonlinear waves during typhoons is shown, and the waves in the straits and the northeastern part of the South China Sea show strong nonlinear characteristics. The high-order spectral model HOS-ocean is employed to simulate the random wave surface series beneath five different platform areas. The waves during the typhoon exhibit strong nonlinear characteristics, and freak waves exist. The space-varying probability model is established to describe the short-term probability distribution of nonlinear wave series. The exceedance probability distributions of the wave surface beneath different platform areas are compared and analyzed. The results show that with an increase in the platform area, the probability of a strong nonlinear wave beneath the platform increases.
基金the supports from National Natural Science Foundation of China(61988101,62073142,22178103)National Natural Science Fund for Distinguished Young Scholars(61925305)International(Regional)Cooperation and Exchange Project(61720106008)。
文摘Crude oil scheduling optimization is an effective method to enhance the economic benefits of oil refining.But uncertainties,including uncertain demands of crude distillation units(CDUs),might make the production plans made by the traditional deterministic optimization models infeasible.A data-driven Wasserstein distributionally robust chance-constrained(WDRCC)optimization approach is proposed in this paper to deal with demand uncertainty in crude oil scheduling.First,a new deterministic crude oil scheduling optimization model is developed as the basis of this approach.The Wasserstein distance is then used to build ambiguity sets from historical data to describe the possible realizations of probability distributions of uncertain demands.A cross-validation method is advanced to choose suitable radii for these ambiguity sets.The deterministic model is reformulated as a WDRCC optimization model for crude oil scheduling to guarantee the demand constraints hold with a desired high probability even in the worst situation in ambiguity sets.The proposed WDRCC model is transferred into an equivalent conditional value-at-risk representation and further derived as a mixed-integer nonlinear programming counterpart.Industrial case studies from a real-world refinery are conducted to show the effectiveness of the proposed method.Out-of-sample tests demonstrate that the solution of the WDRCC model is more robust than those of the deterministic model and the chance-constrained model.
基金funded by the National Natural Science Foundation of China(42174011)。
文摘This paper considers setting different dips for different sub-faults to fit the actual rupture situation based on the fault rupture of the 2013 Lushan M_(S)7.0 earthquake.Meanwhile,combined with the coseismic GNSS data of the Lushan earthquake,the source parameters and sliding distribution of the Lushan earthquake fault are inversed.Firstly,we use the gradient based optimizer(GBO)in nonlinear inversion to obtain the source parameters of this seismic fault.The inversion results indicate that the strike of the fault is 206.52°,the dip is 44.10°,the length is 21.92 km,and the depth is 12.79 km.To refine the sliding distribution of the seismic fault,the seismic fault is divided into 3×3 sub-faults.Then,we fix the central sub-fault dip of 44.10°;the dip of other sub-faults is obtained by iteration.After that,the model is further divided into a fault layer model composed of 23×19 sub fault slices,and using the Matlab fitting function is used to fit the dip of the 23×19 sub faults.Finally,the Lushan seismic fault plane is established as a shovel structure with steep upper and gentle lower,steep south and gentle north.The slip distribution inversion results indicate that the depth of the slip peak is 13 km,the corresponding maximum slip momentum is 0.67 m,the seismic moment is 1.10×10^(19)N·m and the corresponding moment magnitude is MW6.66.The results above are consistent with the research results of seismology.
基金the National Natural Science Foundation of China(Nos.12375123,11975091,and 12305130)the Natural Science Foundation of Henan Province(No.242300421048)+1 种基金China Postdoctoral Science Foundation(No.2023M731016)Henan Postdoctoral Foundation(No.HN2022164).
文摘Neutron-skin thickness is a key parameter for a neutron-rich nucleus;however,it is difficult to determine.In the framework of the Lanzhou Quantum Molecular Dynamics(LQMD)model,a possible probe for the neutron-skin thickness(δ_(np))of neutron-rich ^(48)Ca was studied in the 140A MeV ^(48)Ca+^(9)Be projectile fragmentation reaction based on the parallel momentum distribution(p∥)of the residual fragments.A Fermi-type density distribution was employed to initiate the neutron density distributions in the LQMD simulations.A combined Gaussian function with different width parameters for the left side(Γ_(L))and the right side(Γ_(R))in the distribution was used to describe the p∥of the residual fragments.Taking neutron-rich sulfur isotopes as examples,Γ_(L) shows a sensitive correlation withδ_(np) of ^(48)Ca,and is proposed as a probe for determining the neutron skin thickness of the projectile nucleus.
基金Natural Science Foundation of Shandong Province(Grant No.ZR2022ME180),the National Natural Science Foundation of China(Grant No.51705267).
文摘Snow interacting with a high-speed train can cause the formation of ice in the train bogie region and affect its safety.In this study,a wind-snow multiphase numerical approach is introduced for high-speed train bogies on the basis of the Euler-Lagrange discrete phase model.A particle-wall impact criterion is implemented to account for the presence of snow particles on the surface.Subsequently,numerical simulations are conducted,considering various snow particle diameter distributions and densities.The research results indicate that when the particle diameter is relatively small,the distribution of snow particles in the bogie cavity is relatively uniform.However,as the particle diameter increases,the snow particles in the bogie cavity are mainly located in the rear wheel pairs of the bogie.When the more realistic Rosin-Rammler diameter distribution is applied to snow particles,the positions of snow particles with different diameters vary in the bogie cavity.More precisely,smaller diameter particles are primarily located in the front and upper parts of the bogie cavity,while larger diameter snow particles accumulate at the rear and in the lower parts of the bogie cavity.
基金supported by the Key R&D Project of Shaanxi Province,China(2020ZDLNY07-06)the Science and Technology Program of Shaanxi Academy of Sciences(2022k-11).
文摘In recent years,Meloidogyne enterolobii has emerged as a major parasitic nematode infesting many plants in tropical or subtropical areas.However,the regions of potential distribution and the main contributing environmental variables for this nematode are unclear.Under the current climate scenario,we predicted the potential geographic distributions of M.enterolobii worldwide and in China using a Maximum Entropy(MaxEnt)model with the occurrence data of this species.Furthermore,the potential distributions of M.enterolobii were projected under three future climate scenarios(BCC-CSM2-MR,CanESM5 and CNRM-CM6-1)for the periods 2050s and 2090s.Changes in the potential distribution were also predicted under different climate conditions.The results showed that highly suitable regions for M.enterolobii were concentrated in Africa,South America,Asia,and North America between latitudes 30°S to 30°N.Bio16(precipitation of the wettest quarter),bio10(mean temperature of the warmest quarter),and bio11(mean temperature of the coldest quarter)were the variables contributing most in predicting potential distributions of M.enterolobii.In addition,the potential suitable areas for M.enterolobii will shift toward higher latitudes under future climate scenarios.This study provides a theoretical basis for controlling and managing this nematode.
基金The authors wish to thank the Ecosystem Approach to the management of fisheries and the marine environment in the West African Waters(AWA)project.They also acknowledge support from the international joint laboratory ECLAIRS.The Laboratoire de Météorologie Dynamique(LMD)and the Global Challenges Research Fund(GCRF)African Science for Weather Information and Techniques(SWIFT)Programme.NASA,CNES,and ICARE are acknowledged for providing access to CALIOP and Sun photometer AERONET data.
文摘Saharan dust represents more than 50%of the total desert dust emitted around the globe and its radiative effect significantly affects the atmospheric circulation at a continental scale.Previous studies on dust vertical distribution and the Saharan Air Layer(SAL)showed some shortcomings that could be attributed to imperfect representation of the effects of deep convection and scavenging.The authors investigate here the role of deep convective transport and scavenging on the vertical distribution of mineral dust over Western Africa.Using multi-year(2006-2010)simulations performed with the variable-resolution(zoomed)version of the LMDZ climate model.Simulations are compared with aerosol amounts recorded by the Aerosol Robotic Network(AERONET)and with vertical profiles of the Cloud-Aerosol Lidar with Orthogonal Polarization(CALIOP)measurements.LMDZ allows a thorough examination of the respective roles of deep convective transport,convective and stratiform scavenging,boundary layer transport,and advection processes on the vertical mineral dust distribution over Western Africa.The comparison of simulated dust Aerosol Optical Depth(AOD)and distribution with measurements suggest that scavenging in deep convection and subsequent re-evaporation of dusty rainfall in the lower troposphere are critical processes for explaining the vertical distribution of desert dust.These processes play a key role in maintaining a well-defined dust layer with a sharp transition at the top of the SAL and in establishing the seasonal cycle of dust distribution.This vertical distribution is further reshaped offshore in the Inter-Tropical Convergence Zone(ITCZ)over the Atlantic Ocean by marine boundary layer turbulent and convective transport and wet deposition at the surface.
文摘This study is to understand the impact of operating conditions, especially initial operation temperature (T<sub>ini</sub>) which is set in a high temperature range, on the temperature profile of the interface between the polymer electrolyte membrane (PEM) and the catalyst layer at the cathode (i.e., the reaction surface) in a single cell of polymer electrolyte fuel cell (PEFC). A 1D multi-plate heat transfer model based on the temperature data of the separator measured using the thermograph in a power generation experiment was developed to evaluate the reaction surface temperature (T<sub>react</sub>). In addition, to validate the proposed heat transfer model, T<sub>react</sub> obtained from the model was compared with that from the 3D numerical simulation using CFD software COMSOL Multiphysics which solves the continuity equation, Brinkman equation, Maxwell-Stefan equation, Butler-Volmer equation as well as heat transfer equation. As a result, the temperature gap between the results obtained by 1D heat transfer model and those obtained by 3D numerical simulation is below approximately 0.5 K. The simulation results show the change in the molar concentration of O<sub>2</sub> and H<sub>2</sub>O from the inlet to the outlet is more even with the increase in T<sub>ini</sub> due to the lower performance of O<sub>2</sub> reduction reaction. The change in the current density from the inlet to the outlet is more even with the increase in T<sub>ini</sub> and the value of current density is smaller with the increase in T<sub>ini </sub>due to the increase in ohmic over-potential and concentration over-potential. It is revealed that the change in T<sub>react</sub> from the inlet to the outlet is more even with the increase in T<sub>ini</sub> irrespective of heat transfer model. This is because the generated heat from the power generation is lower with the increase in T<sub>ini </sub>due to the lower performance of O<sub>2</sub> reduction reaction.
基金the State Assignment,project 075-00347-19-00(Patterns of the spatiotemporal dynamics of meadow and forest ecosystems in mountainous areas(Russian Western and Central Caucasus)WWF's‘Save the Forest-Home of Raptors’project(2020-2022).
文摘Abiotic factors play an important role in species localisation,but biotic and anthropogenic predictors must also be considered in distribution modelling for models to be biologically meaningful.In this study,we formalised the biotic predictors of nesting sites for four threatened Caucasian vultures by including species distribution models(wild ungulates,nesting tree species)as biotic layers in the vulture Maxent models.Maxent was applied in the R dismo package and the best set of the model parameters were defined in the R ENMeval package.Performance metrics were continuous Boyce index,Akaike's information criterion,the area under receiver operating curve and true skill statistics.We also calculated and evaluated the null models.Kernel density estimation method was applied to assess the overlap of vulture ecological niches in the environmental space.The accessibility of anthropogenic food resources was estimated using the Path Distance measure that considers elevation gradient.The availability of pine forests(Scots Pine)and wild ungulates(Alpine Chamois and Caucasian Goat)contributed the most(29.6%and 34.3%)to Cinereous Vulture(Aegypius monachus)nesting site model.Wild ungulate distribution also contributed significantly(about 46%)to the Bearded Vulture(Gypaetus barbatus)model.This scavenger nests in the highlands of the Caucasus at a minimum distance of 5–10 km from anthropogenic facilities.In contrast,livestock as a food source was most important in colony distribution of Griffon Vulture(Gyps fulvus).The contribution of distances to settlements and agricultural facilities to the model was 45%.The optimal distance from Egyptian Vulture(Neophron percnopterus)nesting sites to settlements was only 3–10 km,to livestock facilities no more than 15 km with the factor contribution of about 57%.Excluding the wild ungulate availability,the ecological niches of studied vultures overlapped significantly.Despite similar foraging and nesting requirements,Caucasian vultures are not pronounced nesting and trophic competitors due to the abundance of nesting sites,anthropogenic food sources and successful niche sharing.
基金support from High-Level Natural ScienceFoundation of Hainan Province of China (Grant No. 2019RC055)National Natural Science Foundation ofChina (Grant No. 51808176) and the Project Funded by the National First-Class Disciplines (PNFD).
文摘Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to predict the stiffness and shear distribution of glulam beams in elastic-plastic stage,and consequently,the failure of such glulam beams cannot be predicted effectively.To address these issues,an analytical method considering material nonlinearity was proposed for glulam beams,and the calculating equations of deflection and shear stress distribution for different failure modes were established.The proposed method was verified by experiments and numerical models under the corresponding conditions.Results showed that the theoretical calculations were in good agreement with experimental and numerical results,indicating that the equations proposed in this paper were reliable and accurate for such glulam beams with wood material in the elastic-plastic stage ignoring the influence of mechanic properties in radial and tangential directions of wood.Furthermore,the experimental results reported by the previous studies indicated that the method was applicable and could be used as a theoretical reference for predicting the failure of glulam beams.
基金Project supported by the National Natural Science Foundation of China(Grant No.11872033)the Beijing Natural Science Foundation,China(Grant No.3172017)。
文摘Micro sliding phenomenon widely exists in the operation process of mechanical systems,and the micro sliding friction mechanism is always a research hotspot.In this work,based on the total reflection method,a measuring device for interface contact behavior under two-dimensional(2D)vibration is built.The stress distribution is characterized by the light intensity distribution of the contact image,and the interface contact behavior in the 2D vibration process is studied.It is found that the vibration angle of the normal direction of the contact surface and its fluctuation affect the interface friction coefficient,the tangential stiffness,and the fluctuation amplitude of the stress distribution.Then they will affect the change of friction state and energy dissipation in the process of micro sliding.Further,an improved micro sliding friction model is proposed based on the experimental analysis,with the nonlinear change of contact parameters caused by the normal contact stress distribution fluctuation taken into account.This model considers the interface tangential stiffness fluctuation,friction coefficient hysteresis,and stress distribution fluctuation,whose simulation results are consistent well with the experimental results.It is found that considering the nonlinear effect of a certain contact parameter alone may bring a greater error to the prediction of friction behavior.Only by integrating multiple contact parameters can the accuracy of friction prediction is improved.
基金This work was supported by the National Natural Science Foundation of China(NSFC)under Grant 61774141.
文摘A physics-based analytical expression that predicts the charge,electrical field and potential distributions along the gated region of the GaN HEMT channel has been developed.Unlike the gradual channel approximation(GCA),the proposed model considers the non-uniform variation of the concentration under the gated region as a function of terminal applied volt-ages.In addition,the model can capture the influence of mobility and channel temperature on the charge distribution trend.The comparison with the hydrodynamic(HD)numerical simulation showed a high agreement of the proposed model with numerical data for different bias conditions considering the self-heating and quantization of the electron concentration.The ana-lytical nature of the model allows us to reduce the computational and time cost of the simulation.Also,it can be used as a core expression to develop a complete physics-based transistorⅣmodel without GCA limitation.
基金Supported by Natural Science Foundation of Hunan Province (2021JJ30375)Natural Science Foundation of Hunan Provincial Department of Education (20A275)Science and Technology Innovation Team Project of Hunan Province (201937924).
文摘Species distribution models have been widely used to explore suitable habitats of species,the impact of climate change on the distribution of suitable habitats of species,and the construction of ecological reserves.This paper introduced species distribution models commonly used in biodiversity analysis,as well as model performance evaluation indexes,challenges in the application of species distribution models,and finally prospected the development trend of research on species distribution models.
基金Project(2012CB720401)supported by the National Basic Research Program of ChinaProject(2011BAC01B02)supported by the National Key Technology R&D Program of China
文摘For studying the carbon thermal reduction rules of titanium in hot metal and providing a theoretical basis for the blast furnace(BF) hearth protection, the distribution behavior of titanium between low-titanium slag system of CaO-SiO2-MgO-Al2O3-TiO2 and hot metal was studied using analytical reagents in a temperature range from 1350 °C to 1600 °C. Through high temperature melting, rapid quenching, chemical analysis and thermodynamic model calculating, the results showed that the increase of reaction temperature, which improved the titanium distribution L(Ti) and lowered the system activity coefficient γsys, leads to the rise of equilibrium constant. Combined with Wagner and congregated electron phase models, the data obtained in distribution experiments were used to fit out the Gibbs free energy formula of titanium carbothermic reduction. Finally, the relations between the contents of Si and Ti in hot metal and the titanium load to reach the minimum w(Ti) for the formation of Ti C were given.
文摘For performance optimization such as placement,interconnect synthesis,and routing, an efficient and accurate interconnect delay metric is critical,even in design tools development like design for yield (DFY) and design for manufacture (DFM). In the nanometer regime, the recently proposed delay models for RLC interconnects based on statistical probability density function (PDF)interpretation such as PRIMO,H-gamma,WED and RLD bridge the gap between accuracy and efficiency. However, these models always require table look-up when operating. In this paper, a novel delay model based on the Birnbaum-Saunders distribution (BSD) is presented. BSD can accomplish interconnect delay estimation fast and accurately without table look-up operations. Furthermore, it only needs the first two moments to match. Experimental results in 90nm technology show that BSD is robust, easy to implement,efficient,and accurate.
基金This work was supported by the major science and technology projects of CNPC during the“14th five-year plan”(Grant number 2021DJ0101)。
文摘The classical source-to-trap petroleum system concept only considers the migration and accumulation of conventional oil and gas in traps driven dominantly by buoyance in a basin,although revised and improved,even some new concepts as composite petroleum system,total petroleum system,total composite petroleum system,were proposed,but they do not account for the vast unconventional oil and gas reservoirs within the system,which is not formed and distributed in traps dominantly by buoyancedriven.Therefore,the petroleum system concept is no longer adequate in dealing with all the oil and gas accumulations in a basin where significant amount of the unconventional oil and gas resources are present in addition to the conventional oil and gas accumulations.This paper looked into and analyzed the distribution characteristics of conventional and unconventional oil/gas reservoirs and their differences and correlations in petroliferous basins in China and North America,and then proposed whole petroleum system(WPS)concept,the WPS is defined as a natural system that encompasses all the conventional and unconventional oil and gas,reservoirs and resources originated from organic matter in source rocks,the geological elements and processes involving the formation,evolution,and distribution of these oil and gas,reservoirs and resources.It is found in the WPS that there are three kinds of hydrocarbons dynamic fields,three kinds of original hydrocarbons,three kinds of reservoir rocks,and the coupling of these three essential elements lead to the basic ordered distribution model of shale oil/gas reservoirs contacting or interbeded with tight oil/gas reservoirs and separated conventional oil/gas reservoirs from source rocks upward,which is expressed as“S\T-C”.Abnormal conditions lead to other three special ordered distribution models:The first is that with shale oil/gas reservoirs separated from tight oil/gas reservoirs.The second is that with two direction ordered distributions from source upward and downward.The third is with lateral distribution from source outside.
基金The National Natural Science Foundation of China (No.70671021)
文摘Due to the fact that the emergency medicine distribution is vital to the quick response to urgent demand when an epidemic occurs, the optimal vaccine distribution approach is explored according to the epidemic diffusion rule and different urgency degrees of affected areas with the background of the epidemic outbreak in a given region. First, the SIQR (susceptible, infected, quarantined,recovered) epidemic model with pulse vaccination is introduced to describe the epidemic diffusion rule and obtain the demanded vaccine in each pulse. Based on the SIQR model, the affected areas are clustered by using the self-organizing map (SOM) neutral network to qualify the results. Then, a dynamic vaccine distribution model is formulated, incorporating the results of clustering the affected areas with the goals of both reducing the transportation cost and decreasing the unsatisfied demand for the emergency logistics network. Numerical study with twenty affected areas and four distribution centers is carried out. The corresponding numerical results indicate that the proposed approach can make an outstanding contribution to controlling the affected areas with a relatively high degree of urgency, and the comparison results prove that the performance of the clustering method is superior to that of the non-clustering method on controlling epidemic diffusion.