CN112966895A - Dynamic calibration method and device for crowd risk area based on internal energy and information entropy - Google Patents

Dynamic calibration method and device for crowd risk area based on internal energy and information entropy Download PDF

Info

Publication number
CN112966895A
CN112966895A CN202110142459.5A CN202110142459A CN112966895A CN 112966895 A CN112966895 A CN 112966895A CN 202110142459 A CN202110142459 A CN 202110142459A CN 112966895 A CN112966895 A CN 112966895A
Authority
CN
China
Prior art keywords
crowd
internal energy
grid
entropy
risk
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110142459.5A
Other languages
Chinese (zh)
Other versions
CN112966895B (en
Inventor
赵荣泳
王妍
贾萍
刘琼
李翠玲
张智舒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN202110142459.5A priority Critical patent/CN112966895B/en
Publication of CN112966895A publication Critical patent/CN112966895A/en
Application granted granted Critical
Publication of CN112966895B publication Critical patent/CN112966895B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Primary Health Care (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于内能和信息熵的人群风险区域动态标定方法及装置,所述方法包括以下步骤:构建人群运动流体动力学模型,基于该模型进行行人流动仿真;根据仿真结果,基于热焓原理计算人群内能;根据所述人群内能和仿真结果,基于热熵原理标定人群风险区域。与现有技术相比,本发明具有风险评估准确性高等优点。

Figure 202110142459

The invention relates to a method and device for dynamic calibration of crowd risk areas based on internal energy and information entropy. The method includes the following steps: constructing a crowd movement fluid dynamics model, and performing pedestrian flow simulation based on the model; The enthalpy principle is used to calculate the crowd internal energy; according to the crowd internal energy and the simulation results, the crowd risk area is demarcated based on the thermal entropy principle. Compared with the prior art, the present invention has the advantages of high risk assessment accuracy.

Figure 202110142459

Description

Dynamic calibration method and device for crowd risk area based on internal energy and information entropy
Technical Field
The invention relates to a crowd risk area judgment method, in particular to a dynamic crowd risk area calibration method and device based on internal energy and information entropy.
Background
Any disaster occurs because of the outbreak of potential risks in the system, and accurate risk assessment is carried out on the control places, so that the system has great significance for timely and effectively controlling the flow of large-scale pedestrians. Existing research methods are broadly divided into three categories: 1) statistical analysis based on experience. According to the method, various factors influencing the risk of the crowd are provided through a large amount of literature collection, and various factors are comprehensively considered, so that risk assessment is given. Such methods are basically a qualitative determination of risk. 2) State evaluation method based on experiments. The method mainly extracts and analyzes the flow data of the crowd on site, or simulates the flow of the crowd in a simulation mode, extracts the state characteristics of the speed, the density and the like of the crowd and carries out risk assessment and management and control guidance. 3) The quantitative determination method based on theoretical derivation quantifies the crowd movement indexes to form a risk criterion according to mature analysis theory.
To date, there are several deficiencies in determining the risk area of a population: 1) at present, most of the risk judgment and research of crowd evacuation stays on the theoretical basis, and the actual achievement falls to the ground. And the crowd gathering risk is defined from the theoretical derivation point of view, and a recognized and effective judgment method is not formed. 2) Overcrowding of people, i.e., high crowd density or fast movement of people, is generally considered to be a major factor in the development of pedaling events. However, it is not comprehensive to consider only the density and the speed, and it is difficult to accurately determine the risk region.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a dynamic crowd risk area calibration method and device based on internal energy and information entropy.
The purpose of the invention can be realized by the following technical scheme:
a dynamic calibration method for a crowd risk area based on internal energy and information entropy comprises the following steps:
constructing a crowd movement fluid dynamic model, and performing pedestrian flow simulation based on the model;
according to the simulation result, the internal energy of the crowd is calculated based on the enthalpy principle;
and calibrating the crowd risk area based on the thermal entropy principle according to the crowd internal energy and the simulation result.
Further, the crowd motion fluid dynamics model is represented as:
Figure BDA0002929375830000021
Figure BDA0002929375830000022
Figure BDA0002929375830000023
Figure BDA0002929375830000024
where ρ ist(x, y, t) represents pedestrian flow density,
Figure BDA0002929375830000025
indicates the pedestrian flow, Ue(x, y, ρ) represents a monotonically decreasing function of pedestrian stream density, c (x, y, t) is a cost function, and φ (x, y, t) is a cost bit function.
Further, the calculation of the crowd internal energy based on the enthalpy principle specifically includes:
dividing a research scene into networks, and calculating the energy of each grid by taking the grid as a unit, wherein the calculation formula is as follows:
Figure BDA0002929375830000026
wherein, R is the grid crowd internal energy, and rho and v respectively represent the crowd density and the crowd speed in the grid.
Further, the calibrating the crowd risk area based on the principle of thermal entropy specifically comprises:
calculating the motion entropy of the crowd in each grid, obtaining a crowd risk value by adopting a crowd aggregation risk evaluation model, and calibrating a crowd risk area, wherein the crowd aggregation risk evaluation model is expressed as follows:
Figure BDA0002929375830000027
wherein, High _ Ri,jFor the calculated crowd risk value, i and j are pedestrian coordinates corresponding to x and y axes, rho and v respectively represent the crowd density and the crowd speed in the grid, R is the internal energy of the crowd in the grid, and R is the internal energy of the crowd in the gridmaxIs the maximum of the internal energy, α + β ═ 1, p (d)1),p(d2),…,p(dn) For the crowd in the grid along 8 directions d1,d2,…,d8Probability of movement, Ei,jmax is the maximum value of the population motion entropy.
Further, the calculation formula of the crowd motion entropy is as follows.
Figure BDA0002929375830000028
The invention also provides a dynamic calibration device for the crowd risk area based on the internal energy and the information entropy, which comprises the following components:
the simulation module is used for constructing a crowd movement fluid dynamic model and carrying out pedestrian flow simulation based on the model;
the internal energy calculation module is used for calculating the internal energy of the crowd based on the enthalpy principle according to the simulation result;
and the calibration module is used for calibrating the crowd risk area based on the thermal entropy principle according to the crowd internal energy and the simulation result.
Further, in the simulation module, the constructed crowd movement fluid dynamic model is represented as:
Figure BDA0002929375830000031
u(x,y,t):=Ue(x,y,ρ)
Figure BDA0002929375830000032
Figure BDA0002929375830000033
where ρ ist(x, y, t) represents pedestrian flow density,
Figure BDA0002929375830000034
indicates the pedestrian flow, Ue(x, y, ρ) represents a monotonically decreasing function of pedestrian stream density, c (x, y, t) is a cost function, and φ (x, y, t) is a cost bit function.
Further, in the internal energy calculation module, the internal energy of the crowd calculated based on the enthalpy principle is specifically:
dividing a research scene into networks, and calculating the energy of each grid by taking the grid as a unit, wherein the calculation formula is as follows:
Figure BDA0002929375830000035
wherein, R is the grid crowd internal energy, and rho and v respectively represent the crowd density and the crowd speed in the grid.
Further, in the calibration module, calibrating the crowd risk region based on the principle of thermal entropy specifically includes:
calculating the motion entropy of the crowd in each grid, obtaining a crowd risk value by adopting a crowd aggregation risk evaluation model, and calibrating a crowd risk area, wherein the crowd aggregation risk evaluation model is expressed as follows:
Figure BDA0002929375830000036
wherein, High _ Ri,jFor the calculated crowd risk value, i and j are pedestrian coordinates corresponding to x and y axes, rho and v respectively represent the crowd density and the crowd speed in the grid, R is the internal energy of the crowd in the grid, and R is the internal energy of the crowd in the gridmaxIs the maximum value of internal energy, α + β ═1,p(d1),p(d2),…,p(dn) For the crowd in the grid along 8 directions d1,d2,…,d8Probability of movement, Ei,jmax is the maximum value of the population motion entropy.
Further, the calculation formula of the crowd motion entropy is as follows:
Figure BDA0002929375830000037
compared with the prior art, the invention has the following beneficial effects:
(1) an effective crowd gathering risk judgment method does not exist in the prior art, and factors considered by the existing risk evaluation model are not comprehensive enough. According to the invention, a more comprehensive risk dynamic evaluation model is established, so that the large-scale crowd risk area dynamic calibration is carried out, the reliability is high, and the application value is very high. The invention has certain research significance for further mining the internal mechanism of the movement of the crowd in the public place, can provide real-time guidance for crowd evacuation, and is more beneficial to timely managing and controlling the risk crowd.
(2) The crowd asphyxia degree is defined by applying the system internal energy principle, the three-dimensional characteristics of crowd density distribution, speed and movement direction are fully considered, the crowd enthalpy is provided to represent the crowd movement internal energy on the basis of the enthalpy, and the risk assessment accuracy is improved.
(3) High density populations can experience trample disasters and at relatively low densities, excessive chaotic population movements can also create dangerous events. The method is based on the information entropy principle, defines the motion entropy to express the chaos degree of the crowd motion, and improves the accuracy of risk assessment.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a schematic diagram of a grid for high risk area calibration of a population in accordance with the present invention;
fig. 3 is a schematic diagram of an evacuation layout of a waiting hall of a railway station according to an embodiment of the present invention;
FIG. 4 is a conventional population aggregation risk distribution based on population density distribution in an embodiment of the present invention;
FIG. 5 is a diagram illustrating a population group risk distribution based on a hybrid risk assessment model according to an embodiment of the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments. The present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation manner and a specific operation process are given, but the scope of the present invention is not limited to the following embodiments.
The invention provides a dynamic crowd risk area calibration method based on internal energy and information entropy, which takes a time axis as a driving shaft, reflects the change relation of plane distribution of crowd risks along with time, can set and calculate simulation step length according to specific precision requirements, and provides crowd enthalpy based on thermodynamic enthalpy to represent crowd density degree; motion entropy is provided based on thermodynamic entropy to represent the degree of crowd motion disorder, and the method comprises the following steps: constructing a crowd movement fluid dynamic model, and performing pedestrian flow simulation based on the model; according to the simulation result, the internal energy of the crowd is calculated based on the enthalpy principle; and calibrating the crowd risk area based on the thermal entropy principle according to the crowd internal energy and the simulation result.
The method comprises the following steps: crowd moving fluid dynamics modeling
In a macroscopic model of crowd evacuation, pedestrian flow is analogized into fluid, based on the assumption of continuous media, a mass conservation equation of density, speed and flow changing along with time and space is established, a partial differential equation continuously dependent on time and space is obtained, and a pedestrian flow model based on hydrodynamics is formed, wherein the model consists of two parts: continuity equations and pedestrian path selection conditions. The pedestrian path selection is to select the total instantaneous travel cost to the destination to be the minimum according to a reactive dynamic balance distribution principle, so that in a path selection balance state, a defined cost function meets an Eikonal equation, and pedestrians advance at the fastest speed reduced by a cost bit function along the direction of the negative gradient of the cost function. The constructed crowd movement fluid dynamic model is represented as:
Figure BDA0002929375830000051
u(x,y,t):=Ue(x,y,ρ) (2)
Figure BDA0002929375830000052
Figure BDA0002929375830000053
where ρ ist(x, y, t) represents pedestrian flow density,
Figure BDA0002929375830000054
indicates the pedestrian flow, Ue(x, y, ρ) represents a monotonically decreasing function of pedestrian flow density, c (x, y, t) is a cost function representing local travel cost per unit distance of travel, depending on the pedestrian operating conditions of the study scenario itself, and φ (x, y, t) is a cost bit function representing the total instantaneous travel cost to the destination.
Step two: crowd internal energy calculation based on enthalpy principle
Enthalpy is an important state parameter for characterizing the system energy of a substance in thermodynamics, and population enthalpy can be calculated based on enthalpy definition to indicate that population accumulation energy is high:
Figure BDA0002929375830000055
wherein E iskCan be used
Figure BDA0002929375830000056
The calculation is carried out, namely, the kinetic energy is used for replacing the internal energy,
Figure BDA0002929375830000057
instead of the internal pressure, F is the force and A is the area, then
Figure BDA0002929375830000058
In the above formula, the average weight m of the pedestrians in different countries is 65kg, and the maximum speed v of the crowd is 1.34m/s2In grid units. In the formula, ρ is the population density of the unit grid (1m × 1m), and V and a are the volume and area of the unit grid, respectively, and are both set to 1. Therefore, the formula for applying the enthalpy formula to the crowd energy calculation is as follows:
Figure BDA0002929375830000059
it represents the enthalpy of each grid, and when the energy in the crowd is too high, the system is easy to trample once disturbance occurs. According to the definition of the enthalpy, the internal energy distribution of the crowd is calculated, and the energy size of the internal energy distribution is related to the crowd density, the crowd speed and the crowd acceleration size in the grid.
Step three: crowd risk calibration based on thermal entropy principle
Entropy is derived from thermodynamics and is later introduced into the science of information, where entropy is often used to measure uncertainty and disorder in random variables. The thermal entropy theorem of the invention provides a linear combination of crowd energy and the chaos degree of crowd movement directions, and a crowd risk value is obtained by adopting a crowd aggregation risk evaluation model to calibrate a crowd risk area.
According to the crowd flow model, the direction of crowd movement depends on the decreasing direction of the gradient of the cost function. Fig. 1 is a schematic diagram of a grid calibrated in a high risk region of a population, wherein the gray value of the grid represents the cost function value of the grid, and the larger the gray value, the higher the cost function value. The movement direction of the crowd is divided into 8 directions d1,d2,…,d8The probability of occurrence is labeled as p (d)1),p(d2),…,p(dn) Therefore, the formula for calculating the motion entropy of the crowd is (8), and the cost function is d3The descending gradient of the direction is greatest. On the basis, the invention provides a more scientific and comprehensive crowdAnd (4) aggregating the risk assessment model as shown in formula (11).
Figure BDA0002929375830000061
Figure BDA0002929375830000062
Figure BDA0002929375830000063
Figure BDA0002929375830000064
Where α + β is 1, and α and β are empirical values selected according to different scene characteristics. High _ Ri,jThe calculated population risk value is in the range of (0, 1), and a larger value represents a larger risk of population aggregation, and a dangerous event such as treading is likely to occur. In this embodiment, the crowd risk value is higher than 0.8, namely, the crowd risk value is calibrated to be a high risk area, and important management and control and continuous observation are required.
The above functions, if implemented in the form of software functional units and sold or used as a separate product, may be stored in a computer-readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another embodiment, a dynamic calibration apparatus for a crowd risk area based on internal energy and information entropy is provided, which includes: the simulation module is used for constructing a crowd movement fluid dynamic model and carrying out pedestrian flow simulation based on the model; the internal energy calculation module is used for calculating the internal energy of the crowd based on the enthalpy principle according to the simulation result; and the calibration module is used for calibrating the crowd risk area based on the thermal entropy principle according to the crowd internal energy and the simulation result.
Examples
In the embodiment, the effectiveness of the method is verified by taking a train station waiting hall scene as a case. The area of the waiting hall of the railway station in the embodiment is about 11340 square meters, 1 ten thousand waiting people can be accommodated at the same time, and the arrangement of the waiting hall for evacuation is shown in fig. 3.
In the waiting hall scene, the initial number of people is 10000, the people are uniformly distributed in three waiting areas and are evacuated to eight emergency exits, and the simulation step length is 200. Fig. 4 is a result of a conventional crowd gathering risk model based on crowd density distribution, and fig. 5 is a result of a mixed risk assessment model based on internal energy and information entropy proposed by the present invention. By contrast, the population density near obstacles, intersections, and exits is not very high, but the population risk values in these areas are high due to the confusion of the population movement direction and speed.
The foregoing detailed description of the preferred embodiments of the invention has been presented. It should be understood that numerous modifications and variations could be devised by those skilled in the art in light of the present teachings without departing from the inventive concepts. Therefore, the technical solutions available to those skilled in the art through logic analysis, reasoning and limited experiments based on the prior art according to the concept of the present invention should be within the scope of protection defined by the claims.

Claims (10)

1.一种基于内能和信息熵的人群风险区域动态标定方法,其特征在于,包括以下步骤:1. a crowd risk area dynamic calibration method based on internal energy and information entropy, is characterized in that, comprises the following steps: 构建人群运动流体动力学模型,基于该模型进行行人流动仿真;Build a fluid dynamics model of crowd movement, and perform pedestrian flow simulation based on this model; 根据仿真结果,基于热焓原理计算人群内能;According to the simulation results, the internal energy of the crowd is calculated based on the enthalpy principle; 根据所述人群内能和仿真结果,基于热熵原理标定人群风险区域。According to the internal energy of the crowd and the simulation results, the crowd risk area is demarcated based on the principle of thermal entropy. 2.根据权利要求1所述的基于内能和信息熵的人群风险区域动态标定方法,其特征在于,所述人群运动流体动力学模型表示为:2. the crowd risk area dynamic calibration method based on internal energy and information entropy according to claim 1, is characterized in that, described crowd movement fluid dynamics model is expressed as:
Figure FDA0002929375820000011
Figure FDA0002929375820000011
u(x,y,t):=Ue(x,y,ρ)u(x, y, t):=U e (x, y, ρ)
Figure FDA0002929375820000012
Figure FDA0002929375820000012
Figure FDA0002929375820000013
Figure FDA0002929375820000013
其中,ρt(x,y,t)表示行人流密度,
Figure FDA0002929375820000014
表示行人流流量,Ue(x,y,ρ)表示行人流密度的单调递减函数,c(x,y,t)为费用函数,φ(x,y,t)为费用位函数。
Among them, ρ t (x, y, t) represents the pedestrian flow density,
Figure FDA0002929375820000014
Represents pedestrian flow, U e (x, y, ρ) represents the monotonically decreasing function of pedestrian flow density, c(x, y, t) is the cost function, and φ(x, y, t) is the cost function.
3.根据权利要求1所述的基于内能和信息熵的人群风险区域动态标定方法,其特征在于,所述基于热焓原理计算人群内能具体为:3. the crowd risk area dynamic calibration method based on internal energy and information entropy according to claim 1, is characterized in that, the described crowd internal energy calculation based on enthalpy principle is specifically: 对研究场景进行网络划分,以网格为单位,计算每个网格的能量大小,计算公式为:The research scene is divided into networks, and the grid is used as the unit to calculate the energy size of each grid. The calculation formula is:
Figure FDA0002929375820000015
Figure FDA0002929375820000015
其中,R为网格人群内能,ρ和v分别表示网格内人群密度和人群速度。Among them, R is the internal energy of the grid crowd, and ρ and v represent the crowd density and crowd velocity in the grid, respectively.
4.根据权利要求1所述的基于内能和信息熵的人群风险区域动态标定方法,其特征在于,所述基于热熵原理标定人群风险区域具体为:4. the crowd risk area dynamic calibration method based on internal energy and information entropy according to claim 1, is characterized in that, described crowd risk area calibration based on thermal entropy principle is specifically: 计算每个网格内的人群运动熵,采用人群聚集风险评估模型获得人群风险值,标定人群风险区域,所述人群聚集风险评估模型表示为:Calculate the crowd movement entropy in each grid, use the crowd aggregation risk assessment model to obtain the crowd risk value, and demarcate the crowd risk area. The crowd aggregation risk assessment model is expressed as:
Figure FDA0002929375820000016
Figure FDA0002929375820000016
其中,High_Ri,j为求得的人群风险值,i,j为x,y轴对应的行人坐标,ρ和v分别表示网格内人群密度和人群速度,R为网格人群内能,Rnax为内能最大值,α+β=1,p(d1),p(d2),…,p(dn)为网格内人群沿8个方向d1,d2,…,d8运动的概率,Ei,jmax为人群运动熵最大值。Among them, High_R i,j is the obtained crowd risk value, i,j is the pedestrian coordinate corresponding to the x,y axis, ρ and v represent the crowd density and crowd velocity in the grid, respectively, R is the internal energy of the grid crowd, R nax is the maximum value of internal energy, α+β=1, p(d 1 ), p(d 2 ),...,p(d n ) are the population in the grid along 8 directions d 1 ,d 2 ,...,d 8 The probability of movement, E i,j max is the maximum value of crowd movement entropy.
5.根据权利要求4所述的基于内能和信息熵的人群风险区域动态标定方法,其特征在于,所述人群运动熵的计算公式为:5. the crowd risk area dynamic calibration method based on internal energy and information entropy according to claim 4, is characterized in that, the calculation formula of described crowd movement entropy is:
Figure FDA0002929375820000021
Figure FDA0002929375820000021
6.一种基于内能和信息熵的人群风险区域动态标定装置,其特征在于,包括:6. A crowd risk area dynamic calibration device based on internal energy and information entropy, is characterized in that, comprising: 仿真模块,用于构建人群运动流体动力学模型,基于该模型进行行人流动仿真;The simulation module is used to construct a fluid dynamics model of crowd movement, and perform pedestrian flow simulation based on the model; 内能计算模块,用于根据仿真结果,基于热焓原理计算人群内能;The internal energy calculation module is used to calculate the internal energy of the crowd based on the enthalpy principle according to the simulation results; 标定模块,用于根据所述人群内能和仿真结果,基于热熵原理标定人群风险区域。The calibration module is used for calibrating the crowd risk area based on the thermal entropy principle according to the crowd internal energy and the simulation result. 7.根据权利要求6所述的基于内能和信息熵的人群风险区域动态标定装置,其特征在于,所述仿真模块中,构建的人群运动流体动力学模型表示为:7. the crowd risk area dynamic calibration device based on internal energy and information entropy according to claim 6, is characterized in that, in described simulation module, the crowd movement fluid dynamics model of construction is expressed as:
Figure FDA0002929375820000022
Figure FDA0002929375820000022
u(x,y,t):=Ue(x,y,ρ)u(x,y,t):=U e (x,y,ρ)
Figure FDA0002929375820000023
Figure FDA0002929375820000023
Figure FDA0002929375820000024
Figure FDA0002929375820000024
其中,ρt(x,y,t)表示行人流密度,
Figure FDA0002929375820000025
表示行人流流量,Ue(x,y,ρ)表示行人流密度的单调递减函数,c(x,y,t)为费用函数,φ(x,y,t)为费用位函数。
Among them, ρ t (x, y, t) represents the pedestrian flow density,
Figure FDA0002929375820000025
Represents pedestrian flow, U e (x, y, ρ) represents the monotonically decreasing function of pedestrian flow density, c(x, y, t) is the cost function, and φ(x, y, t) is the cost function.
8.根据权利要求6所述的基于内能和信息熵的人群风险区域动态标定装置,其特征在于,所述内能计算模块中,基于热焓原理计算人群内能具体为:8. The crowd risk area dynamic calibration device based on internal energy and information entropy according to claim 6, characterized in that, in the internal energy calculation module, calculating crowd internal energy based on the enthalpy principle is specifically: 对研究场景进行网络划分,以网格为单位,计算每个网格的能量大小,计算公式为:The research scene is divided into networks, and the grid is used as the unit to calculate the energy size of each grid. The calculation formula is:
Figure FDA0002929375820000026
Figure FDA0002929375820000026
其中,R为网格人群内能,ρ和v分别表示网格内人群密度和人群速度。Among them, R is the internal energy of the grid crowd, and ρ and v represent the crowd density and crowd velocity in the grid, respectively.
9.根据权利要求6所述的基于内能和信息熵的人群风险区域动态标定装置,其特征在于,所述标定模块中,基于热熵原理标定人群风险区域具体为:9. The crowd risk area dynamic calibration device based on internal energy and information entropy according to claim 6, is characterized in that, in the described calibration module, the crowd risk area calibration based on thermal entropy principle is specifically: 计算每个网格内的人群运动熵,采用人群聚集风险评估模型获得人群风险值,标定人群风险区域,所述人群聚集风险评估模型表示为:Calculate the crowd movement entropy in each grid, use the crowd aggregation risk assessment model to obtain the crowd risk value, and demarcate the crowd risk area. The crowd aggregation risk assessment model is expressed as:
Figure FDA0002929375820000027
Figure FDA0002929375820000027
其中,High_Ri,j为求得的人群风险值,i,j为x,y轴对应的行人坐标,ρ和v分别表示网格内人群密度和人群速度,R为网格人群内能,Rmax为内能最大值,α+β=1,p(d1),p(d2),…,p(dn)为网格内人群沿8个方向d1,d2,…,d8运动的概率,Ei,jmax为人群运动熵最大值。Among them, High_R i,j is the obtained crowd risk value, i,j is the pedestrian coordinate corresponding to the x,y axis, ρ and v represent the crowd density and crowd velocity in the grid, respectively, R is the internal energy of the grid crowd, R max is the maximum value of internal energy, α+β=1, p(d 1 ), p(d 2 ),...,p(d n ) are the population in the grid along 8 directions d 1 ,d 2 ,...,d 8 The probability of movement, E i,j max is the maximum value of crowd movement entropy.
10.根据权利要求9所述的基于内能和信息熵的人群风险区域动态标定装置,其特征在于,所述人群运动熵的计算公式为:10. The crowd risk area dynamic calibration device based on internal energy and information entropy according to claim 9, wherein the calculation formula of the crowd movement entropy is:
Figure FDA0002929375820000031
Figure FDA0002929375820000031
CN202110142459.5A 2021-02-02 2021-02-02 Dynamic calibration method and device for crowd risk area based on internal energy and information entropy Active CN112966895B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110142459.5A CN112966895B (en) 2021-02-02 2021-02-02 Dynamic calibration method and device for crowd risk area based on internal energy and information entropy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110142459.5A CN112966895B (en) 2021-02-02 2021-02-02 Dynamic calibration method and device for crowd risk area based on internal energy and information entropy

Publications (2)

Publication Number Publication Date
CN112966895A true CN112966895A (en) 2021-06-15
CN112966895B CN112966895B (en) 2022-09-16

Family

ID=76271813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110142459.5A Active CN112966895B (en) 2021-02-02 2021-02-02 Dynamic calibration method and device for crowd risk area based on internal energy and information entropy

Country Status (1)

Country Link
CN (1) CN112966895B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109241950A (en) * 2018-10-19 2019-01-18 杭州电子科技大学 A kind of crowd panic state identification method based on enthalpy Distribution Entropy
CN109408889A (en) * 2018-09-21 2019-03-01 同济大学 Macroscopical crowd panic measure and its application based on comentropy
CN109472049A (en) * 2018-09-29 2019-03-15 同济大学 Establishment method and application of dynamic model of macro-population panic propagation dynamic model

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109408889A (en) * 2018-09-21 2019-03-01 同济大学 Macroscopical crowd panic measure and its application based on comentropy
CN109472049A (en) * 2018-09-29 2019-03-15 同济大学 Establishment method and application of dynamic model of macro-population panic propagation dynamic model
CN109241950A (en) * 2018-10-19 2019-01-18 杭州电子科技大学 A kind of crowd panic state identification method based on enthalpy Distribution Entropy

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HAODONG YIN等: "An energy based method to measure the crowd safety", 《TRANSPORTATION RESEARCH PROCEDIA》 *
LING HUANG等: "Revisiting Hughes’ dynamic continuum model for pedestrian flow and the development of an efficient solution algorithm", 《TRANSPORTATION RESEARCH PART B》 *
李小红: "城市轨道交通车站客流群集规律及拥挤踩踏风险研究", 《中国博士学位论文全文数据库 工程科技II辑》 *

Also Published As

Publication number Publication date
CN112966895B (en) 2022-09-16

Similar Documents

Publication Publication Date Title
Yan et al. Distributionally consistent simulation of naturalistic driving environment for autonomous vehicle testing
Song et al. Optimization of wind turbine micro-siting for reducing the sensitivity of power generation to wind direction
CN101807224A (en) Mesoscopic-microcosmic integrated traffic simulation vehicle flow loading method
CN115512538B (en) Variable speed limiting control method and system based on carbon emission
CN107220447A (en) Crowd evacuation emulation method and system based on evacuation set of paths
Xie et al. Robust and efficient prediction of the collection efficiency in icing accretion simulation for 3D complex geometries using the Lagrangian approach I: An adaptive interpolation method based on the restricted radial basis functions
CN117725812B (en) A distributed robust optimization method for passenger evacuation routes in subway stations during flooding
CN117351521A (en) Digital twinning-based power transmission line bird detection method, system, medium and equipment
CN111291098A (en) Track mode mining method and device
CN115470731A (en) Method and system for predicting wind power of wind field based on microclimate
Sarmady et al. Simulation of pedestrian movements using fine grid cellular automata model
CN118334869A (en) Traffic state estimation method and equipment based on traffic flow theory and CTM
CN109190325A (en) Crowd evacuation Path Planning Simulation method based on the analysis of pedestrian's crowding
CN104573264B (en) The method that simulation airborne vehicle passes through low area
Xu et al. PIGAT: Physics-informed graph attention transformer for air traffic state prediction
CN105138773B (en) Pedestrian simulation method and device based on discrete element emulation platform
CN103226893A (en) Method and simulation system for describing traffic flow behavior based on two-dimensional macroscopic flow model
CN112966895A (en) Dynamic calibration method and device for crowd risk area based on internal energy and information entropy
CN112966363B (en) Random disturbance-oriented crowd flow simulation method and device in waiting hall of transportation hub
Xiang et al. Using SIR model to simulate emotion contagion in dynamic crowd aggregation process
CN109408889B (en) Macro-crowd panic measurement method based on information entropy and its application
Yang et al. The research of multi-resolution modeling and simulation of the emergency evacuation
CN117116096A (en) Airport delay prediction method and system based on multi-channel traffic images and deep CNN
CN106682390A (en) Lead-coating ice increase simulating computation method and system
Papathanasopoulou et al. Flexible car-following models incorporating information from adjacent lanes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant