CN105740510A - Simulation system and method of evacuation crowd behavior based on grid-density-relation - Google Patents

Simulation system and method of evacuation crowd behavior based on grid-density-relation Download PDF

Info

Publication number
CN105740510A
CN105740510A CN201610044959.4A CN201610044959A CN105740510A CN 105740510 A CN105740510 A CN 105740510A CN 201610044959 A CN201610044959 A CN 201610044959A CN 105740510 A CN105740510 A CN 105740510A
Authority
CN
China
Prior art keywords
grid
crowd
scene
evacuation
relation
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
CN201610044959.4A
Other languages
Chinese (zh)
Other versions
CN105740510B (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.)
Shandong Normal University
Original Assignee
Shandong Normal 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 Shandong Normal University filed Critical Shandong Normal University
Priority to CN201610044959.4A priority Critical patent/CN105740510B/en
Publication of CN105740510A publication Critical patent/CN105740510A/en
Application granted granted Critical
Publication of CN105740510B publication Critical patent/CN105740510B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • 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
    • G06Q50/265Personal security, identity or safety

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Development Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Security & Cryptography (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Alarm Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a simulation system and method of an evacuation crowd behavior based on grid-density-relation. The system comprises an evacuation scene model building unit, an evacuation scene global path planning unit, a crowd activity generation unit and a crowd simulation unit, wherein the evacuation scene model building unit is used for extracting characteristics of an evacuation scene according to structural information of the evacuation scene to obtain a three-dimensional model of the evacuation scene; the evacuation scene global path planning unit used for performing global path planning on the evacuation scene according to the three-dimensional model of the evacuation scene, and calculating global paths of all entrances in the evacuation scene; the crowd activity generation unit used for performing crowd grouping according to the quantity of all entrances, a relationship between individuals in the crowd and individual location information in the evacuation scene; and the crowd simulation unit used for calculating the real-time speed of individuals in a simulation process according to a social force model, and in the simulation process, correcting the moving speed of each individual in the same group in real time, so as to keep the advancing consistency of each group, and realize behavior simulation of the evacuation crowd.

Description

Evacuation behavior simulation system and method based on grid-density-relation
Technical field
The present invention relates to evacuation emulation field, particularly relate to a kind of Evacuation behavior simulation system based on grid-density-relation And method.
Background technology
Large-scale gathering of people is modern city, or even a key character of modern society.The place of densely populated place is hidden The biggest potential safety hazard, is susceptible to accident, once causes the fear of crowd, and the accident of crowded trample occurs, and injures and deaths are past Toward the most heavy.Accelerating development especially with Urbanization in China, the activity of large-scale personnel's high aggregation increases, people The crowded trample accident of group is easier to occur.
Crowd evacuation under accident is a complicated system engineering, generally uses the mode evacuating experiment in actual applications Obtain and suitably evacuate prediction scheme.Which has the features such as with strong points, informative.But, owing to there is personal security Cannot ensure, test the input inevitably problem such as big, Computer Simulation becomes having most of crowd evacuation under research accident The method of effect.
Computer Simulation crowd evacuation problem is utilized to need in view of numerous factors affecting group movement, according to social psychology Viewpoint, colony can be divided into the small group with different characteristic.And inside a small group, there is some common between individuality Feature formed a group.Such as, modal group includes family relation group, classmate's friend group, colleague group etc.. These groups are that tie connects by one or more social relations.
Observable it turns out that, in the crowd evacuation under critical situation, the behavior of same group personnel always keeps consistent.He Reciprocally negotiate the direction of motion, and colony may be left the most together.So, each individuality in same group Can attempt to adjust the movement velocity of oneself, the concordance advanced with holding with whole group.
The group of relation between the most non-consideration crowd and the formation impact on motion in existing crowd evacuation emulation method.Although, Helbing and Molnar proposes social force model and simulates crowd behaviour according to Newtonian mechanics formula, it is possible to reappear some well existing As, such as " more think the slowest " and the arch phenomenon in exit.But the behavior of pedestrian of its over-simplification, does not examine Relation and the group of the formation impact on motion between worry crowd.
Dirk Helbing, according to crowd behaviour feature, establishes social force model based on Newtonian mechanics.Social force refers to one The power applied it by local environment (including the people in environment and thing) when people moves, not acts directly on him The power of physical concept.According to the different motivation of pedestrian and influence in the environment in social force model, have three kinds of active forces Impact: active force between driving force, people and people and people and the intermolecular forces of barrier.The force action of these power is in row People, produces an acceleration.In whole individual walking process, and between individuality and individuality, there is the work of certain power all the time With.Such as, driving force can guide individual towards target direction advance;Before individual Body contact, the active force between people and people Make the individuality in crowd avoid mutually to collide;Active force between human and environment makes the individuality in crowd avoid and barrier collision. This stage can explain by classical Newton's second law.
Therefore, in order to solve crowd's safety and evacuating personnel and carry out the problem that crowd evacuation is drilled safely, need one to consider The crowd behaviour analogue system and method thereof that group affects moving, reaches, according to detection building evacuation performance, to safeguard colony's row For order, optimize evacuation process and improve the purpose of evacuation efficiency.
Summary of the invention
The invention provides a kind of Evacuation behavior simulation system and method based on grid-density-relation.The present invention is examining Consider the characteristic of human relation, by the relation value of weighting and the density value of weighting, crowd is carried out packet and simulate crowd evacuation appearance Naturally packet and the arch phenomenon in crowd exit, for the safety instructing crowd evacuation to drill.
For achieving the above object, the present invention is by the following technical solutions:
A kind of Evacuation behavior simulation system based on grid-density-relation, including:
Evacuating model of place and set up unit, it for carrying out extracting the feature of evacuation scene according to the structural information evacuating scene, To the threedimensional model evacuating scene;
Evacuating scene global path planning unit, it is for carrying out global path according to the threedimensional model evacuating scene to evacuating scene Planning, calculates and evacuates the global path of all gateways in scene;
Crowd activity's signal generating unit, it is for according to evacuating in scene in all gateways quantity, crowd between individuality relation and individual Body position information carries out crowd's packet;
Crowd simulation unit, it is for calculating real-time speed individual in simulation process according to social force model, in the mistake of emulation Cheng Zhong, revises the movement velocity with individuality each in group, in real time with the concordance keeping each group to advance, it is achieved Evacuation Behavior simulation.
Described crowd activity's signal generating unit, including:
Crowd's individuality positional information generation module, it is crowd activity in the threedimensional model setting up demographic data collection and evacuate scene Mapping relations between region, generate the individual position in crowd activity region of crowd according to mapping relations;
Evacuate scene grid and divide module, determine crowd's number of packet, to evacuation according to evacuating all gateways quantity in scene In the threedimensional model of scene, crowd activity region carries out grid division;
Lattice relationship-density value computing module, it is for according to the individual amount in each grid cell and the relation between individuality Value, calculates mesh-density value and the lattice relationship value of each grid cell, and is weighted mesh-density value and lattice relationship value Superposition, it is thus achieved that the lattice relationship-density value of each grid cell;
Crowd's grouping module, its for by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual Body and the physical distance of group switching centre, be assigned to individual place grid in the group at each boundary fitting place, complete the packet of crowd.
A kind of emulation mode of Evacuation behavior simulation system based on grid-density-relation, including:
Step (1): carry out extracting the feature of evacuation scene according to the structural information evacuating scene, obtains evacuating the three-dimensional mould of scene Type;
Step (2): the threedimensional model according to evacuating scene carries out global path planning to evacuating scene, calculates and evacuates institute in scene The global path of there is any discrepancy mouth;
Step (3): determine according to density in evacuating scene of relation between individuality, crowd in crowd and individual positional information Crowd's grouping information;
Step (4): calculate real-time speed individual in simulation process according to social force model, during emulation, in real time Revise the movement velocity with individuality each in group, with the concordance keeping whole group to advance, it is achieved the behavior simulation of Evacuation.
Described step (3) determines the process of crowd's grouping information, including:
Step (3.1): set up the mapping relations between crowd activity region in the threedimensional model of demographic data collection and evacuation scene, The individual position in crowd activity region of crowd is generated according to mapping relations;
Step (3.2): determine crowd's number of packet, to the three-dimensional evacuating scene according to evacuating all gateways quantity in scene In model, crowd activity region carries out grid division;
Step (3.3): according to the individual amount in each grid cell and the relation value between individuality, calculates each grid cell Mesh-density value and lattice relationship value, and mesh-density value and lattice relationship value are weighted superposition, it is thus achieved that each grid list Lattice relationship-the density value of unit;
Step (3.4): by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual with in group The physical distance of the heart, is assigned to individual place grid in the group at each boundary fitting place, completes the packet of crowd.
In described step (3.4), individual place grid is assigned to the process in the group at each boundary fitting place, including:
Step (3.4.1): according to the size of the lattice relationship-density value of each grid cell, from big to small grid cell is carried out Sequence, is formed and treats packet queue;
Step (3.4.2): choose and treat that in packet queue, first grid is as boundary fitting, starts from second grid until treating point Last grid in forming a team to arrange, determines whether the adjacent grid of boundary fitting one by one, the most then this grid is incorporated to core Grid, and treating that packet queue leaves out this grid;Will merge after boundary fitting from treating that packet queue moves on to packet queue, and Mark group number;
Step (3.4.3): repeating step (3.4.2), until treating that in packet queue, meshes number is 0, packet terminates;Otherwise enter Enter next step;
Step (3.4.4): treat that in packet queue, meshes number is non-zero, and packet queue meshes number is equal to fixed packet Number, then according to treating each individuality and the physical distance of group switching centre in packet queue in the grid of packet queue, divide corresponding grid It is fitted in the group at each boundary fitting place, until treating that in packet queue, meshes number is 0.
In described step (3.3), the weight being weighted mesh-density value and lattice relationship value is respectively 0.5 and 0.5.
Described step (3.4.1) treats that in packet queue, the density value of grid is nonzero value.
In described step (3.4), packet count is equal to 3 times that export number of plane domain to be evacuated.
In described step (3.4), the unit grid number dividing crowd activity region in the threedimensional model of evacuation scene is equal to 10 times of packet count.
In described step (3.3), the lattice relationship value of each grid cell is equal to the relation between all individualities in this grid cell Value sum is divided by 2.
The invention have the benefit that
(1) the Evacuation behavior simulation system based on grid-density-relation of the present invention consider the relation between crowd and The group the formed impact on motion, carries out packet by the density value of relation value and the weighting of weighting and simulates crowd and dredge crowd Shed existing packet and the arch phenomenon in crowd exit naturally, crowd evacuation is drilled safely to the foundation providing important, also may be used To detect building evacuation performance, optimize actual evacuation process and improve evacuation efficiency;
(2) region such as the crowd evacuation emulation method of the present invention is applicable to have door plane room and large-scale plane venue, and wouldn't Consider the impact of barrier, the plane domain such as square of multiple outlet can be expanded to.
Accompanying drawing explanation
Fig. 1 is the structural representation of the Evacuation behavior simulation system based on grid-density-relation of the present invention;
Fig. 2 is the schematic flow sheet of the Evacuation behavior simulation method based on grid-density-relation of the present invention;
Fig. 3 is crowd activity's administrative division map of the evacuation scene of four doors in the embodiment of the present invention;
Fig. 4 is that the crowd activity region evacuating scene in the embodiment of the present invention is divided into 120 grids;
Fig. 5 is mesh-density schematic diagram in the embodiment of the present invention;
Fig. 6 is lattice relationship value schematic diagram in the embodiment of the present invention;
Fig. 7 is the whole adjacent grid schematic diagram of boundary fitting in the embodiment of the present invention;
Fig. 8 is that in the embodiment of the present invention, the face of boundary fitting adjoins grid schematic diagram;
Fig. 9 is that in the embodiment of the present invention, the point of boundary fitting adjoins grid schematic diagram;
Figure 10 is the crowd's schematic diagram after initializing in the embodiment of the present invention;
Figure 11 is crowd's grouping process schematic diagram in crowd evacuation emulation in the embodiment of the present invention;
Figure 12 is to be grouped backward outlet in the embodiment of the present invention to move, close to the state diagram of outlet.
Detailed description of the invention
The present invention will be further described with embodiment below in conjunction with the accompanying drawings:
In the Evacuation behavior simulation system based on grid-density-relation of the present invention:
(1) demographic data collection P={pi, i=1,2 ..., N}, pi=(pi1,pi2,…,pim) represent that in data set, i-th is individual, pij Represent the jth attribute that i-th is individual.
(2) relation between crowd's individuality:
For two individual plAnd pm, relation value represents the tightness degree between two individualities, and value is 0 to 1, with rel (pi, pm) Represent.
If two individualities are irrelevant, then relation value is 0, if mothers and sons or lover's relation, then relation value levels off to 1.In order to During reflection emergency evacuation, search guide or the leader of a group, everyone and the guide in tourist party is set up pass It is, and relation is 1, same reason, between individuality and the leader of each collective, the most all opening relationships, and relation value It is 1.Treat that the relation span between the individuality of Evacuation data set is:
(3) the individual mapping relations with the crowd activity region evacuating scene of crowd:
The present invention evacuates the crowd activity region of scene and is defined as S=L × W, is an a length of L, the rectangle plane district of a width of W Territory, point can (x y) represents above each with the coordinate being made up of two numbers.
Wherein, P → S be mapped as M: treat the Evacuation data set P mapping M={p to plane domaini(x,y)};
piIt is to treat that the i-th in Evacuation data set is individual, 1≤i≤N, (x, y) is the coordinate on two dimensional surface, 1≤x≤W, 1≤y≤L, pi(x, y) represent the individual position on plane domain of i-th be (x, y).
(4) crowd's packet count and the determination of number of meshes:
According to the outlet number E in the crowd activity region evacuating scene, determine packet count T and number of meshes G to be divided: as Shown in Fig. 3, evacuating scene plane domain has 4 outlets, packet count T=E × 3, mainly considers a corresponding outlet, can There to be the group in three directions, left, middle and right to evacuate to outlet.
Grid number G=T × 10, this value is to be drawn by experience, it is also possible to dynamically adjust.
Wherein, the process of stress and strain model is:
Plane domain S is evenly dividing for h × k rectangular mesh (being divided into square), h × k=G as far as possible.
Division methods is:
S is an a length of L, and the rectangle plane region of a width of W is divided into G grid cell, the face of the most each grid cell Long-pending a=(L × W)/G, the length of side isH=L/e, k=W/e.
As shown in Figure 4, such as, one long by 300, and the rectangle plane region of wide 250 is divided into 120 grids, each net The area a=(300 × 250)/120, a=625 of lattice, the length of sideE=25, region is laterally being divided into H=300/25, h=12, longitudinally divided k=10.
(5) grid cell:
Without loss of generality, if the coordinate in the plane domain lower left corner is (0,0), each rectangular mesh ci,jSit with the lower left corner and the upper right corner Mark is expressed asxi=xi-1+ e, yj=yj-1+ e, 1≤i≤h, 1≤j≤k, the most each rectangle net Sound of laughing be a grid cell.
(6) mesh-density:
Mesh-density is the individual sum in this grid cell.
For given grid cell ci,j, 1≤i≤h, 1≤j≤k, statistics has how many each and every one body drops in this grid cell. Individual pi(x y) falls at grid ci,jIn, the coordinate range of x, y existJudge in interval, i.e. xi-1< x≤xi, yi-1< y≤yi
As it is shown in figure 5, mesh-density den (ci,j)=count (pi(x,y)),xi-1< x≤xi, yi-1< y≤yi, count is meter Number function.
(7) relation value of grid:
The relation value of grid: for given grid cell ci,j, 1≤i≤h, grid cell ci,jRelation value between interior individuality it With divided by 2.
As shown in Figure 6, lattice relationship value r e l ( c i , j ) = Σ ( x 1 , y 1 ) ∈ c i , j ( x 2 , y 2 ) ∈ c i , j r e l ( p l ( x 1 , y 1 ) , p m ( x 2 , y 2 ) ) / 2 , xi-1< x1≤xi, yi-1< y1≤yi, xi-1< x2≤xi, yi-1< y2≤yi
Grid cell ci,jRelation value sum between interior individuality is divided by 2, if being because pl(x1,y1) and pm(x2,y2) relevant System, then pm(x2,y2) and pl(x1,y1) also there is relation, and be symmetrical relations, therefore, during summation, pl(x1,y1) with pm(x2,y2Relation between) calculates twice.
(8) lattice relationship-density value:
Lattice relationship-density value is the relation value density value sum with weighting of weighting,
I.e. RelDen (ci,j)=w1×rel(ci,j)+w2×den(ci,j)。
Wherein, w1,w2It is the weight of relation and density respectively, typically takes 0.5,0.5, if higher, then to packet rate request w2Weight increase, w1Weight reduce.
According to the rate request of the crowd's packet to crowd evacuation emulation model, the weighting weight of lattice relationship value is more than 0.5, grid The weighting weight of density value is less than 0.5.
(9) relationship gap value:
Two individual pl(x1,y1), pm(x2,y2Relationship gap value between) is the relation value physical distance value with weighting of weighting Sum, it may be assumed that
RelDist(pl(x1,y1),pm(x2,y2))=w1×rel(pl(x1,y1),pm(x2,y2))+w2×dist(pl(x1,y1),pm(x2,y2)) Wherein: d i s t ( p l ( x 1 , y 1 ) , p m ( x 2 , y 2 ) ) = ( x 2 - x 1 ) 2 + ( y 2 - y 1 ) 2 .
As it is shown in figure 1, the Evacuation behavior simulation system based on grid-density-relation of the present invention, including:
Evacuating model of place and set up unit, it for carrying out extracting the feature of evacuation scene according to the structural information evacuating scene, To the threedimensional model evacuating scene;
Evacuating scene global path planning unit, it is for carrying out global path according to the threedimensional model evacuating scene to evacuating scene Planning, calculates and evacuates the global path of all gateways in scene;
Crowd activity's signal generating unit, it is for according to evacuating in scene in all gateways quantity, crowd between individuality relation and individual Body position information carries out crowd's packet;
Crowd simulation unit, it is for calculating real-time speed individual in simulation process according to social force model, in the mistake of emulation Cheng Zhong, revises the movement velocity with individuality each in group, in real time with the concordance keeping each group to advance, it is achieved Evacuation Behavior simulation.
Further, crowd activity's signal generating unit, including:
Crowd's individuality positional information generation module, it is crowd activity in the threedimensional model setting up demographic data collection and evacuate scene Mapping relations between region, generate the individual position in crowd activity region of crowd according to mapping relations;
Evacuate scene grid and divide module, determine crowd's number of packet, to evacuation according to evacuating all gateways quantity in scene In the threedimensional model of scene, crowd activity region carries out grid division;
Lattice relationship-density value computing module, it is for according to the individual amount in each grid cell and the relation between individuality Value, calculates mesh-density value and the lattice relationship value of each grid cell, and is weighted mesh-density value and lattice relationship value Superposition, it is thus achieved that the lattice relationship-density value of each grid cell;
Crowd's grouping module, its for by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual Body and the physical distance of group switching centre, be assigned to individual place grid in the group at each boundary fitting place, complete the packet of crowd.
As in figure 2 it is shown, the emulation mode of the Evacuation behavior simulation system based on grid-density-relation of the present invention, including:
Step (1): carry out extracting the feature of evacuation scene according to the structural information evacuating scene, obtains evacuating the three-dimensional mould of scene Type;
Step (2): the threedimensional model according to evacuating scene carries out global path planning to evacuating scene, calculates and evacuates institute in scene The global path of there is any discrepancy mouth;
Step (3): determine according to density in evacuating scene of relation between individuality, crowd in crowd and individual positional information Crowd's grouping information;
Step (4): calculate real-time speed individual in simulation process according to social force model, during emulation, in real time Revise the movement velocity with individuality each in group, with the concordance keeping whole group to advance, it is achieved the behavior simulation of Evacuation.
Further, step (3) determines the process of crowd's grouping information, including:
Step (3.1): set up the mapping relations between crowd activity region in the threedimensional model of demographic data collection and evacuation scene, The individual position in crowd activity region of crowd is generated according to mapping relations;
Step (3.2): determine crowd's number of packet, to the three-dimensional evacuating scene according to evacuating all gateways quantity in scene In model, crowd activity region carries out grid division;
Step (3.3): according to the individual amount in each grid cell and the relation value between individuality, calculates each grid cell Mesh-density value and lattice relationship value, and mesh-density value and lattice relationship value are weighted superposition, it is thus achieved that each grid list Lattice relationship-the density value of unit;
Step (3.4): by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual with in group The physical distance of the heart, is assigned to individual place grid in the group at each boundary fitting place, completes the packet of crowd.
The grid higher than all of its neighbor lattice relationship-density value is referred to as boundary fitting, with core (ci,j) represent.Wherein, adjacent grid It is divided into the most adjacent grid, face to adjoin grid and the adjacent grid of point, the most as shown in Figure 7, Figure 8 and Figure 9.
Treat each element wq in packet queue WaitQueueiA non-NULL grid cell, i.e. den (wqi)≠0。
Element fq in packet queue FinishQueueiIt is it has been determined that the grid cell of affiliated group.
In step (3.4), individual place grid is assigned to the process in the group at each boundary fitting place, including:
The first step: sequence: according to the lattice relationship-density value of weighting, toward low, grid is ranked up from high, forms queue WaitQueue, has u element in row.Density value is the grid of zero, does not i.e. have individuality in this grid, does not enter queue.
Second step: boundary fitting core (ci,j) choose and merge: choose and treat first grid conduct in packet queue WaitQueue Boundary fitting, from the beginning of second grid, proceeds to the u-1 grid, compares with boundary fitting one by one, if certain Grid is the adjacent grid of boundary fitting, then this grid is incorporated to boundary fitting, u=u-1;After having circulated, the core after merging Grid is from treating that packet queue WaitQueue moves on to packet queue FinishQueue, and marks group number, u=u-1;
3rd step: repeat second step, until treating in packet queue WaitQueue that meshes number is 0 or packet queue FinishQueue meshes number is T;
4th step: if treating that packet queue WaitQueue meshes number is 0, the grid in packet queue FinishQueue the most It is exactly to individual packet, turns to the 5th step;Otherwise, if meshes number is T in FinishQueue, and in WaitQueue Meshes number is not 0, then, from WaitQueue, take out grid one by one, and to individuality each in grid, calculate one by one with Central point personal relationship-the distance value of boundary fitting in FinishQueue, is incorporated to the core minimum with central point personal relationship-distance value Heart grid, until WaitQueue element number is 0.
5th step: end of packet process.
Provide below an emulation embodiment:
By carrying out compuman group's evacuation emulation on the plane domain of 300 people 300 × 250, as shown in Figure 10~Figure 12.At figure In 10, the state after initialization, the individual expression of same grayscale: there is relation between individuality.In fig. 11, divide from crowd After group, each group of state to outlet movement is it can be seen that the individuality of same grayscale divides in identical group.Can from Figure 12 Going out, in outlet moving process, related individuality substantially remains in same group, and can try with each individuality in group Figure adjusts the movement velocity of oneself, the concordance advanced with whole group with holding.
Although the detailed description of the invention of the present invention is described by the above-mentioned accompanying drawing that combines, but not limit to scope System, one of ordinary skill in the art should be understood that on the basis of technical scheme, and those skilled in the art need not pay Go out various amendments or deformation that creative work can make still within protection scope of the present invention.

Claims (10)

1. an Evacuation behavior simulation system based on grid-density-relation, it is characterised in that including:
Evacuating model of place and set up unit, it for carrying out extracting the feature of evacuation scene according to the structural information evacuating scene, To the threedimensional model evacuating scene;
Evacuating scene global path planning unit, it is for carrying out global path according to the threedimensional model evacuating scene to evacuating scene Planning, calculates and evacuates the global path of all gateways in scene;
Crowd activity's signal generating unit, it is for according to evacuating in scene in all gateways quantity, crowd between individuality relation and individual Body position information carries out crowd's packet;
Crowd simulation unit, it is for calculating real-time speed individual in simulation process according to social force model, in the mistake of emulation Cheng Zhong, revises the movement velocity with individuality each in group, in real time with the concordance keeping each group to advance, it is achieved Evacuation Behavior simulation.
A kind of Evacuation behavior simulation system based on grid-density-relation, it is characterised in that Described crowd activity's signal generating unit, including:
Crowd's individuality positional information generation module, it is crowd activity in the threedimensional model setting up demographic data collection and evacuate scene Mapping relations between region, generate the individual position in crowd activity region of crowd according to mapping relations;
Evacuate scene grid and divide module, determine crowd's number of packet, to evacuation according to evacuating all gateways quantity in scene In the threedimensional model of scene, crowd activity region carries out grid division;
Lattice relationship-density value computing module, it is for according to the individual amount in each grid cell and the relation between individuality Value, calculates mesh-density value and the lattice relationship value of each grid cell, and is weighted mesh-density value and lattice relationship value Superposition, it is thus achieved that the lattice relationship-density value of each grid cell;
Crowd's grouping module, its for by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual Body and the physical distance of group switching centre, be assigned to individual place grid in the group at each boundary fitting place, complete the packet of crowd.
3. the emulation of an Evacuation behavior simulation system based on the grid-density-relation as described in claim 1-2 is arbitrary Method, it is characterised in that including:
Step (1): carry out extracting the feature of evacuation scene according to the structural information evacuating scene, obtains evacuating the three-dimensional mould of scene Type;
Step (2): the threedimensional model according to evacuating scene carries out global path planning to evacuating scene, calculates and evacuates institute in scene The global path of there is any discrepancy mouth;
Step (3): determine according to density in evacuating scene of relation between individuality, crowd in crowd and individual positional information Crowd's grouping information;
Step (4): calculate real-time speed individual in simulation process according to social force model, during emulation, in real time Revise the movement velocity with individuality each in group, with the concordance keeping whole group to advance, it is achieved the behavior simulation of Evacuation.
4. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 3, it is special Levy and be, described step (3) determines the process of crowd's grouping information, including:
Step (3.1): set up the mapping relations between crowd activity region in the threedimensional model of demographic data collection and evacuation scene, The individual position in crowd activity region of crowd is generated according to mapping relations;
Step (3.2): determine crowd's number of packet, to the three-dimensional evacuating scene according to evacuating all gateways quantity in scene In model, crowd activity region carries out grid division;
Step (3.3): according to the individual amount in each grid cell and the relation value between individuality, calculates each grid cell Mesh-density value and lattice relationship value, and mesh-density value and lattice relationship value are weighted superposition, it is thus achieved that each grid list Lattice relationship-the density value of unit;
Step (3.4): by the grid that lattice relationship-density value is maximum be boundary fitting to build group switching centre, according to individual with in group The physical distance of the heart, is assigned to individual place grid in the group at each boundary fitting place, completes the packet of crowd.
5. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 4, it is special Levy and be, in described step (3.4), individual place grid is assigned to the process in the group at each boundary fitting place, including:
Step (3.4.1): according to the size of the lattice relationship-density value of each grid cell, from big to small grid cell is carried out Sequence, is formed and treats packet queue;
Step (3.4.2): choose and treat that in packet queue, first grid is as boundary fitting, starts from second grid until treating point Last grid in forming a team to arrange, determines whether the adjacent grid of boundary fitting one by one, the most then this grid is incorporated to core Grid, and treating that packet queue leaves out this grid;Will merge after boundary fitting from treating that packet queue moves on to packet queue, and Mark group number;
Step (3.4.3): repeating step (3.4.2), until treating that in packet queue, meshes number is 0, packet terminates;Otherwise enter Enter next step;
Step (3.4.4): treat that in packet queue, meshes number is non-zero, and packet queue meshes number is equal to fixed packet Number, then according to treating each individuality and the physical distance of group switching centre in packet queue in the grid of packet queue, divide corresponding grid It is fitted in the group at each boundary fitting place, until treating that in packet queue, meshes number is 0.
6. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 4, it is special Levying and be, in described step (3.3), the weight being weighted mesh-density value and lattice relationship value is respectively 0.5 and 0.5.
7. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 5, it is special Levy and be, described step (3.4.1) treats that in packet queue, the density value of grid is nonzero value.
8. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 4, it is special Levying and be, in described step (3.4), packet count is equal to 3 times that export number of plane domain to be evacuated.
9. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 4, it is special Levy and be, in described step (3.4), the unit grid number that crowd activity region in the threedimensional model of evacuation scene is divided Equal to packet count 10 times.
10. the emulation mode of Evacuation behavior simulation system based on grid-density-relation as claimed in claim 4, its It is characterised by, in described step (3.3), between the lattice relationship value of each grid cell is equal to all individualities in this grid cell Relation value sum divided by 2.
CN201610044959.4A 2016-01-22 2016-01-22 Evacuation crowd behaviour analogue system based on grid-density-relationship and its method Expired - Fee Related CN105740510B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610044959.4A CN105740510B (en) 2016-01-22 2016-01-22 Evacuation crowd behaviour analogue system based on grid-density-relationship and its method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610044959.4A CN105740510B (en) 2016-01-22 2016-01-22 Evacuation crowd behaviour analogue system based on grid-density-relationship and its method

Publications (2)

Publication Number Publication Date
CN105740510A true CN105740510A (en) 2016-07-06
CN105740510B CN105740510B (en) 2018-08-31

Family

ID=56247559

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610044959.4A Expired - Fee Related CN105740510B (en) 2016-01-22 2016-01-22 Evacuation crowd behaviour analogue system based on grid-density-relationship and its method

Country Status (1)

Country Link
CN (1) CN105740510B (en)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202751A (en) * 2016-07-15 2016-12-07 山东师范大学 Crowd evacuation emulation method that Cultural Algorithm is combined with particle swarm optimization and system
CN106815422A (en) * 2017-01-11 2017-06-09 山东师范大学 Consider the emulation mode and system of alleviating excessive compression in dense population of arching
CN107256307A (en) * 2017-06-09 2017-10-17 山东师范大学 The crowd evacuation emulation method and system of knowledge based navigation
CN107423480A (en) * 2017-05-19 2017-12-01 山东师范大学 A kind of crowd movement's emulation mode and system
CN107463751A (en) * 2017-08-10 2017-12-12 山东师范大学 A kind of crowd based on DBSCAN clustering algorithms by half is grouped evacuation emulation method and system
CN108153966A (en) * 2017-12-22 2018-06-12 中山大学 The real-time evacuation emulation method of multiple exit potential energy field model crowd based on GPU
CN108416381A (en) * 2018-03-01 2018-08-17 国家海洋局第海洋研究所 A kind of multi-density clustering method towards three-dimensional point set
CN109063268A (en) * 2018-07-05 2018-12-21 同济大学 The iterative emulation mode of crowd evacuation macromodel
CN109101694A (en) * 2018-07-16 2018-12-28 山东师范大学 A kind of the crowd behaviour emulation mode and system of the guidance of safe escape mark
CN110688749A (en) * 2019-09-19 2020-01-14 山东师范大学 Crowd evacuation simulation method and system
CN111339668A (en) * 2020-02-28 2020-06-26 西南交通大学 Crowd evacuation visualization method based on emotion cognition
CN111400914A (en) * 2020-03-17 2020-07-10 山东师范大学 Vision-field-based crowd grouping and crowd evacuation simulation system and method
CN113221379A (en) * 2021-05-28 2021-08-06 视伴科技(北京)有限公司 Method and device for planning number of viewing seats of event activity
CN113902168A (en) * 2021-09-08 2022-01-07 长江信达软件技术(武汉)有限责任公司 Flood emergency evacuation method based on real-time crowd data fusion
CN114545932A (en) * 2022-01-28 2022-05-27 北京交通大学 Path correction method, device and equipment and readable storage medium
CN115017742A (en) * 2022-08-08 2022-09-06 西安深信科创信息技术有限公司 Automatic driving test scene generation method, device, equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090306946A1 (en) * 2008-04-08 2009-12-10 Norman I Badler Methods and systems for simulation and representation of agents in a high-density autonomous crowd
CN103679264A (en) * 2013-12-23 2014-03-26 山东师范大学 Crowd evacuation path planning method based on artificial fish swarm algorithm
US20140372348A1 (en) * 2011-12-15 2014-12-18 Northeastern University Real-time anomaly detection of crowd behavior using multi-sensor information
CN104239636A (en) * 2014-09-16 2014-12-24 北京航空航天大学 Fire emergency evacuation simulation method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090306946A1 (en) * 2008-04-08 2009-12-10 Norman I Badler Methods and systems for simulation and representation of agents in a high-density autonomous crowd
US20140372348A1 (en) * 2011-12-15 2014-12-18 Northeastern University Real-time anomaly detection of crowd behavior using multi-sensor information
CN103679264A (en) * 2013-12-23 2014-03-26 山东师范大学 Crowd evacuation path planning method based on artificial fish swarm algorithm
CN104239636A (en) * 2014-09-16 2014-12-24 北京航空航天大学 Fire emergency evacuation simulation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JOAO E.ALMEIDA等: "Crowd Simulation Modeling Applied to Emergency and Evacuation Simulations using Multi-Agent Systems", 《COMPUTER SCIENCE》 *
胡春雨: "基于群体智能的人群疏散建模与仿真方法研究", 《中国优秀硕士学位论文全文数据库(信息科技辑)》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202751B (en) * 2016-07-15 2019-03-19 山东师范大学 Crowd evacuation emulation method and system of the Cultural Algorithm in conjunction with particle swarm optimization
CN106202751A (en) * 2016-07-15 2016-12-07 山东师范大学 Crowd evacuation emulation method that Cultural Algorithm is combined with particle swarm optimization and system
CN106815422A (en) * 2017-01-11 2017-06-09 山东师范大学 Consider the emulation mode and system of alleviating excessive compression in dense population of arching
CN107423480A (en) * 2017-05-19 2017-12-01 山东师范大学 A kind of crowd movement's emulation mode and system
CN107256307A (en) * 2017-06-09 2017-10-17 山东师范大学 The crowd evacuation emulation method and system of knowledge based navigation
CN107463751A (en) * 2017-08-10 2017-12-12 山东师范大学 A kind of crowd based on DBSCAN clustering algorithms by half is grouped evacuation emulation method and system
CN108153966B (en) * 2017-12-22 2021-06-15 中山大学 GPU-based multi-outlet potential energy field model crowd real-time evacuation simulation method
CN108153966A (en) * 2017-12-22 2018-06-12 中山大学 The real-time evacuation emulation method of multiple exit potential energy field model crowd based on GPU
CN108416381B (en) * 2018-03-01 2021-07-16 国家海洋局第一海洋研究所 Multi-density clustering method for three-dimensional point set
CN108416381A (en) * 2018-03-01 2018-08-17 国家海洋局第海洋研究所 A kind of multi-density clustering method towards three-dimensional point set
CN109063268B (en) * 2018-07-05 2021-03-26 同济大学 Iterative simulation method of crowd evacuation macro model
CN109063268A (en) * 2018-07-05 2018-12-21 同济大学 The iterative emulation mode of crowd evacuation macromodel
CN109101694B (en) * 2018-07-16 2019-05-28 山东师范大学 A kind of the crowd behaviour emulation mode and system of the guidance of safe escape mark
CN109101694A (en) * 2018-07-16 2018-12-28 山东师范大学 A kind of the crowd behaviour emulation mode and system of the guidance of safe escape mark
CN110688749A (en) * 2019-09-19 2020-01-14 山东师范大学 Crowd evacuation simulation method and system
CN110688749B (en) * 2019-09-19 2023-04-11 山东师范大学 Crowd evacuation simulation method and system
CN111339668A (en) * 2020-02-28 2020-06-26 西南交通大学 Crowd evacuation visualization method based on emotion cognition
CN111400914A (en) * 2020-03-17 2020-07-10 山东师范大学 Vision-field-based crowd grouping and crowd evacuation simulation system and method
CN113221379A (en) * 2021-05-28 2021-08-06 视伴科技(北京)有限公司 Method and device for planning number of viewing seats of event activity
CN113902168A (en) * 2021-09-08 2022-01-07 长江信达软件技术(武汉)有限责任公司 Flood emergency evacuation method based on real-time crowd data fusion
CN114545932A (en) * 2022-01-28 2022-05-27 北京交通大学 Path correction method, device and equipment and readable storage medium
CN115017742A (en) * 2022-08-08 2022-09-06 西安深信科创信息技术有限公司 Automatic driving test scene generation method, device, equipment and storage medium

Also Published As

Publication number Publication date
CN105740510B (en) 2018-08-31

Similar Documents

Publication Publication Date Title
CN105740510A (en) Simulation system and method of evacuation crowd behavior based on grid-density-relation
CN105740514B (en) A kind of large-scale crowd evacuation computer simulation system and its method
CN106227958B (en) Group evacuation simulation system and method combining artificial bee colony with social force model
Li et al. A grouping method based on grid density and relationship for crowd evacuation simulation
CN107464021B (en) A kind of crowd evacuation emulation method based on intensified learning, device
Djehiche et al. A mean-field game of evacuation in multilevel building
Li et al. Relationship between behavior aggressiveness and pedestrian dynamics using behavior-based cellular automata model
CN108491598A (en) A kind of crowd evacuation emulation method and system based on path planning
CN105468801A (en) Simulation method and system for crowd evacuation in public place
CN109543285B (en) Crowd evacuation simulation method and system integrating data driving and reinforcement learning
CN105808852A (en) Indoor pedestrian microscopic simulation method based on cellular automaton
CN105550484A (en) Dynamic data driven modeling and simulation method for crowd evacuation and transportation behaviors in sudden-onset disaster
CN107423480A (en) A kind of crowd movement's emulation mode and system
CN105701314A (en) Complex crowd evacuation behavior simulation method based on self-adaption intelligent agent model
Zhang et al. Multi-agent based modeling and simulating for evacuation process in stadium
CN108153966A (en) The real-time evacuation emulation method of multiple exit potential energy field model crowd based on GPU
CN108446469A (en) A kind of group's behavior evacuation emulation method and device based on video drive
CN101567017A (en) Urban-evacuation simulation method based on multi-resolution images
CN111339668B (en) Crowd evacuation visualization method based on emotion cognition
Zhang et al. Knowledge-based crowd motion for the unfamiliar environment
Li et al. Lane-design for mixed pedestrian flow in T-shaped passage
Liao et al. A kinetic modeling of crowd evacuation with several groups in complex venues
CN112257258B (en) Crowd motion simulation method based on Boids cluster motion model
CN113299068B (en) Traffic network congestion state prediction method and system
Yamamoto et al. New approach for pedestrian dynamics by real-coded cellular automata (RCA)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180831

Termination date: 20220122

CF01 Termination of patent right due to non-payment of annual fee