CN113077907B - Novel infectious disease road infection probability calculation model - Google Patents
Novel infectious disease road infection probability calculation model Download PDFInfo
- Publication number
- CN113077907B CN113077907B CN202110317372.7A CN202110317372A CN113077907B CN 113077907 B CN113077907 B CN 113077907B CN 202110317372 A CN202110317372 A CN 202110317372A CN 113077907 B CN113077907 B CN 113077907B
- Authority
- CN
- China
- Prior art keywords
- model
- infection
- people
- street
- time
- 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.)
- Active
Links
- 208000015181 infectious disease Diseases 0.000 title claims abstract description 129
- 208000035473 Communicable disease Diseases 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 238000000034 method Methods 0.000 claims abstract description 35
- 241000711573 Coronaviridae Species 0.000 claims abstract description 29
- 241000700605 Viruses Species 0.000 claims abstract description 19
- 238000009792 diffusion process Methods 0.000 claims abstract description 10
- 230000002458 infectious effect Effects 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000011835 investigation Methods 0.000 claims description 9
- 230000002829 reductive effect Effects 0.000 claims description 7
- 201000010099 disease Diseases 0.000 claims description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 241001678559 COVID-19 virus Species 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- OWFJMIVZYSDULZ-PXOLEDIWSA-N (4s,4ar,5s,5ar,6s,12ar)-4-(dimethylamino)-1,5,6,10,11,12a-hexahydroxy-6-methyl-3,12-dioxo-4,4a,5,5a-tetrahydrotetracene-2-carboxamide Chemical compound C1=CC=C2[C@](O)(C)[C@H]3[C@H](O)[C@H]4[C@H](N(C)C)C(=O)C(C(N)=O)=C(O)[C@@]4(O)C(=O)C3=C(O)C2=C1O OWFJMIVZYSDULZ-PXOLEDIWSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000005541 medical transmission Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 230000009385 viral infection Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/80—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for detecting, monitoring or modelling epidemics or pandemics, e.g. flu
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Public Health (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
The invention provides a travel infection probability calculation method suitable for novel coronaviruses, which comprises the following steps: establishing a road infection probability database of the selected area; constructing a novel infectious disease model SEIRK based on the SEIR model; building a street pedestrian flow model; according to the road infection probability database, a novel infectious disease model SEIRK and a street pedestrian flow model are synthesized, a novel infectious disease road infection probability calculation model is constructed, and further the travel infection probability of the novel coronavirus is predicted and calculated. Compared with the traditional virus diffusion model for calculating the infection probability of the infectious viruses, the method comprehensively considers the influence of the density of pedestrians on the road, the wearing of the mask by the pedestrians, the distance from the living place of the ill person and the like on the novel infectious viruses, can reduce the risk of the pedestrians going out to be infected to a greater extent, and has important significance in the aspects of inhibiting epidemic spread and the like.
Description
Technical Field
The invention belongs to the field of infectious disease model construction, and particularly relates to a travel infection probability calculation method suitable for novel coronaviruses.
Background
There are many practical virus diffusion models, however, most of these models are not very specific to viruses, and the main problems are as follows:
1. the diffusion model considers that the content is too single, focuses on the virus itself and ignores the influence of the external environment on the virus, so that the simulated virus diffusion is biased;
2. the diffusion model concentrates on the change of the number of people transmitted along with time, ignores the influence caused by free population flow and national preventive measures, and cannot react according to the data of virus real-time transformation;
3. the diffusion model focuses on macroscopic infection of viruses, and based on evolution of specific data at specific time, the transformation probability of the virus cannot be specified to a certain time;
disclosure of Invention
In order to solve the technical problems, the invention provides a method for calculating the infection probability of a novel coronavirus on a road, which is applicable to a plurality of similar infectious diseases including 2019-nCoV (2019 novel coronavirus), aims to solve the problem that the existing infectious disease spread prediction model can not truly predict the infection probability of the infectious disease on the road, and provides reference for preventing the infectious disease.
The technical scheme for solving the problems is as follows: the travel infection probability calculation method suitable for the novel coronavirus comprises the following steps of:
1) Establishing a street infection probability database of the selected area;
2) Constructing a novel infectious disease model SEIRK based on the SEIR model according to the acquired data to simulate novel infectious disease infection;
3) Constructing a street pedestrian flow model according to the acquired data;
4) And constructing a novel infectious disease traveling infection probability calculation model according to the traveling infection probability database and the comprehensive infection model SEIRK and the street pedestrian flow model.
1. The method for calculating the traveling infection probability applicable to the novel coronavirus comprises the following specific steps of:
1-1) survey the number of patients in the region to be surveyed and the geographical position of the residence thereof according to the time line;
1-2) drawing a distance-infection rate simulation curve according to the data obtained by investigation;
1-2-1) acquisition of patient affected areas
Taking the living place of the patient in the data as the circle center, and taking the infection inner diameter r as the radius to form a circular area: s=pi r 2 This area is the patient affected area.
1-2-2) construction of street infection Rate model
Investigation of the patient affected area as a subject has found that the risk of infection with new coronavirus in the area within 1 meter of the place where the infected person resides is reduced by a factor of 80%; the risk of infection with new coronavirus in the area 1 m or more from the place where the infected person resides is reduced by 50% by index, so that a calculation formula of a street infection rate model (x is distance (unit: m), and Y is infection rate) is obtained as follows:
1-2-3) distance-infection rate simulation curves were plotted
Distance-infection rate simulated curves were plotted with distance on the x-axis and infection rate on the y-axis.
1-3) establishing a road infection probability database according to the data obtained by the distance-infection rate simulation curve, wherein if the influence areas of the patients coincide, the influence values are accumulated.
2. The method for calculating the traveling infection probability applicable to the novel coronavirus comprises the following specific steps of:
2-1) determining model considerations:
2-1-1) S-a susceptible person, i.e., a population susceptible to infection but not yet suffering from disease;
2-1-2) E-viral latency carrier, the latency being the group of people with disease latency who are ill but asymptomatic, the novel coronavirus also being infectious during latency;
2-1-3) I-the infected person, the infected person is the infected person;
2-1-4) R is a rehabilitating person, the rehabilitating person heals the patient, and the model considers that the rehabilitating person cannot be infected again;
2-1-5) K, which is a resistent, namely the number of people wearing the mask in the crowd, and the model considers that the wearing of the mask in epidemic situation can inhibit the spread of viruses to a certain extent, thereby effectively reducing the risk of infecting new coronaviruses;
2-2) according to the consideration, the virus diffusion simulation is carried out by taking time t as a variable, and the formula is calculated as follows:
2-2-1) counting the number of susceptible people on day t
2-2-2) counting the number of people with latent form on day t
2-2-3) counting the number of infected persons on day t
2-2-4) calculating the number of rehabilitative persons on day t
2-2-5) calculating the number of preventable persons on day t
Wherein S (t) is the number of susceptible persons on the t th day, R (t) is the number of recovered persons on the t th day, I (t) is the number of infected persons on the t th day, E (t) is the number of latent persons on the t th day, R is the number of persons contacting the patient when no strict measures are taken in the initial stage, R2 is the number of persons contacting the patient when strict measures are taken after epidemic situation occurs, b is the infection probability when no strict measures are taken in the initial stage, b2 is the infection probability when strict measures are taken after epidemic situation occurs, N is the current time pedestrian flow obtained according to a street pedestrian flow model, y is the recovery probability, and c is the probability of having resistant persons in the crowd (namely the probability of wearing a mask in the crowd).
3. The method for calculating the traveling infection probability applicable to the novel coronavirus comprises the following specific steps of:
3-1) obtaining street traffic
3-1-1) acquiring the number of user locations
First, dividing a street into n regionsA field recording the range O of each area i (i=1, 2, …, n), after the user allows to obtain his longitude and latitude coordinates, looking at which area the user coordinates are in at this time, so as to obtain the people flow data of each area of the street;
3-1-2) directly obtaining the flow of street people
And transmitting the shot video pictures to a server by using the existing cameras of each traffic road section, obtaining the number of pedestrians through a set of embedded intelligent video screen analysis algorithm, and carrying out density conversion according to the number of pedestrians obtained through recognition.
The number of people conversion formula is as follows:
N=m*c
m=p/m 0
wherein N is the total population number obtained by conversion, m is the number of people in a unit length (pedestrian density), p is the average number of people which can be identified by all cameras entering the current street, and m 0 The street length taken by the camera is calculated, and c is the actual street length, so that the people flow number is finally calculated.
3-2) correcting the people flow data with a linear least squares method
3-2-1) software and hardware data combining
The people flow data obtained by software and the people flow data obtained by hardware are combined into a group of two-dimensional data (x t ,y t ), t=1,2,…,l,x t The time t is represented, the people flow data counted by the user positioning is obtained by software, and x is t Different from each other, y t The time t is represented by the people flow data acquired by the hardware equipment, and the time l is represented by the time of acquiring the last group of data; the least square method is applied to obtain a function y=f (x) so that the sum of the square difference of the least square method result and the actual hardware obtained flow of people in the recorded momentMinimum;
3-2-2) the people flow data x obtained by the software at the time t+1 t+1 Substituting f (x), y t+1 Acquiring people flow data for the corrected software;
3-3) time series prediction
The people flow of each street is time, the people flow is predicted by adopting a 'simple sequence time average method' in a time sequence prediction method, and the calculation steps are as follows:
3-3-1) collecting time series sample data a at l' times daily over m days i′j′ Where i 'denotes the date number, i' =1, 2, …, m, j 'denotes the time number, j' =1, 2, …, l, and a is recovered after correction i′j′ ;
3-3-2) calculating the arithmetic mean of the m days l' time instants
3-3-3) predictive computation
When the time series is listed by time, the people flow predictive value y at the j-th time of the predictive date m+1 day m+1,j The method comprises the following steps:
3-4) using the obtained people flow data after the hardware acquired people flow and the software system acquired people flow are processed by the least square method as the result of the people flow model.
4. The method for calculating the traveling infection probability of the novel coronavirus comprises the following specific steps of:
4-1) determining model parameters:
w: acquiring the safety distance between people according to the investigation;
x: pedestrian position is based on the nearest patient residence distance;
t: predicting the occurrence time of the infectious disease on the same day;
pp: acquiring road infection probability data from the travel infection probability database;
i (t): the number of people infected on the t-th day obtained from the new infectious disease model SEIRK;
k (t): the number of preventable persons on the t-th day obtained from the new infectious disease model SEIRK;
e (t): the number of people with latent people on the t-th day obtained from a new infectious disease model SEIRK;
n: and obtaining the pedestrian traffic at the current moment according to the street pedestrian traffic model.
4-2) dynamic and static combination calculation of the probability of travel infection of novel coronaviruses
If x > w, the calculation formula is as follows:
P=(0.5*(2-w)+1)*P′
otherwise, the calculation formula is as follows:
the invention has the beneficial effects that:
1. the method comprehensively considers the influence of multiple factors such as environment, people flow, distance, wearing condition of the mask of pedestrians and the like on the spread of viruses of the infectious diseases, and improves the accuracy of calculating the probability of virus infection.
2. The model is based on the SEIR model, and provides a SEIRK model, and new parameters (resistant) are added to perform diffusion simulation, so that the predicted result is more approximate to the actual value.
3. The pedestrian identifier is introduced to detect pedestrians, the images shot by the cameras are processed mainly by the cameras at the entrances and exits of each street, and the number of passenger flows is calculated by means of density conversion of pedestrians.
4. The street traffic is calculated and corrected by adopting a mode of combining hardware and software, so that the credibility of error increasing data can be reduced.
5. The method can be applied to construction and simulation of a complex virus propagation prediction model, can reduce the risk of pedestrian outing infection, and has important significance in the aspects of inhibiting epidemic spread and the like.
Drawings
FIG. 1 is an overall flow chart of the present invention;
FIG. 2 is a flowchart of a street infection probability database for creating a selected area;
FIG. 3 is a graph of distance-infection rate modeling of the present invention;
FIG. 4 is a flow chart for constructing a novel infection model SEIRK;
FIG. 5 is a simulated graph of an infectious disease transmission curve in an embodiment of the present invention;
FIG. 6 is a flow chart for building a street pedestrian traffic model;
FIG. 7 is a flow chart of a method model for constructing a traveling infection probability calculation for a novel coronavirus.
Detailed description of the preferred embodiments
The invention is further described below with reference to the drawings and examples.
As shown in fig. 1, the present invention provides a method for constructing a method for calculating a probability of traveling infection applicable to a novel coronavirus, comprising the steps of:
1) Building a street infection probability database of the selected area, as shown in fig. 2:
1-1) survey the number of patients in the region to be surveyed and the geographical position of the residence thereof according to the time line;
1-2) drawing a distance-infection rate simulation curve according to the data obtained by investigation;
1-2-1) acquisition of patient affected areas
Taking the living place of the patient in the data as the circle center, and taking the infection inner diameter r as the radius to form a circular area: s=pi r 2 This area is the patient affected area.
1-2-2) construction of street infection Rate model
Investigation of the patient affected area as a subject has found that the risk of infection with new coronavirus in the area within 1 meter of the place where the infected person resides is reduced by a factor of 80%; the risk of infecting a new coronavirus in an area 1 meter or more from the place where the infected person resides is reduced by an index of 50%, so that a street infection rate model calculation formula (x is distance (unit: m), Y is infection rate) is obtained as follows:
1-2-3) distance-infection rate simulation curves were plotted
As shown in fig. 3, a distance-infection rate simulation curve is drawn with the distance as the x-axis and the infection rate as the y-axis.
1-3) establishing a road infection probability database according to the data obtained by the distance-infection rate simulation curve, wherein if the influence areas of the patients coincide, the influence values are accumulated.
2) And constructing a novel infection model SEIRK based on the SEIR model according to the acquired data, wherein the specific flow is shown in figure 4.
2-1) determining model considerations
2-1-1) S-a susceptible person, i.e., a population susceptible to infection but not yet suffering from disease;
2-1-2) E-viral latency carrier, the latency being the group of people with disease latency who are ill but asymptomatic, the novel coronavirus also being infectious during latency;
2-1-3) I-the infected person, the infected person is the infected person;
2-1-4) R is a rehabilitating person, the rehabilitating person heals the patient, and the model considers that the rehabilitating person cannot be infected again;
2-1-5) K, which is a resistent, the resistent is the number of people wearing the mask in the crowd, and the model considers that the wearing of the mask in epidemic situation can inhibit the spread of virus to a certain extent, thereby effectively reducing the risk of infecting new coronavirus.
2-2) according to the consideration, the virus diffusion simulation is carried out by taking time t as a variable, and the formula is calculated as follows:
2-2-1) counting the number of susceptible people on day t
2-2-2) counting the number of people with latent form on day t
2-2-3) counting the number of infected persons on day t
2-2-4) calculating the number of rehabilitative persons on day t
2-2-5) calculating the number of preventable persons on day t
Wherein S (t) is the number of susceptible persons on the t th day, R (t) is the number of recovered persons on the t th day, I (t) is the number of infected persons on the t th day, E (t) is the number of latent persons on the t th day, R is the number of persons contacting the patient when no strict measures are taken in the initial stage, R2 is the number of persons contacting the patient when strict measures are taken after epidemic situation occurs, b is the infection probability when no strict measures are taken in the initial stage, b2 is the infection probability when strict measures are taken after epidemic situation occurs, N is the current time pedestrian flow obtained according to a street pedestrian flow model, y is the recovery probability, c is the probability of having resistant persons in the crowd (i.e. the probability of wearing a mask in the crowd), and a specific simulation curve is shown in FIG. 5.
3) The street pedestrian traffic model is constructed according to the acquired data, as shown in fig. 6, and the specific steps are as follows:
3-1) obtaining street traffic
3-1-1) acquiring the number of user locations
Firstly dividing a street into n areas, and recording the range O of each area i (i=1, 2, …, n), after the user allows to obtain his longitude and latitude coordinates, looking at which area the user coordinates are in at this time, so as to obtain the people flow data of each area of the street;
3-1-2) directly obtaining the street traffic:
and transmitting the shot video pictures to a server by using the existing cameras of each traffic road section, obtaining the number of pedestrians through a set of embedded intelligent video screen analysis algorithm, and carrying out density conversion according to the number of pedestrians obtained through recognition.
The number of people conversion formula is as follows:
N=m*c
m=p/m 0
wherein N is the total population number obtained by conversion, m is the number of people in a unit length (pedestrian density), p is the average number of people which can be identified by all cameras entering the current street, and m 0 The street length taken by the camera is c is the real street length, and finally the people flow number is calculated;
3-2) correcting the people flow data with a linear least squares method
3-2-1) software and hardware data combining
The people flow data obtained by software and the people flow data obtained by hardware are combined into a group of two-dimensional data (x t ,y t ), t=1,2,…,l,x t The time t is represented, the people flow data counted by the user positioning is obtained by software, and x is t Different from each other, y t The time t is represented by the people flow data acquired by the hardware equipment, and the time l is represented by the time of acquiring the last group of data; the least square method is applied to obtain a function y=f (x) so that the sum of the square difference of the least square method result and the actual hardware obtained flow of people in the recorded momentMinimum;
3-2-2) the people flow data x obtained by the software at the time t+1 t+1 Substituting f (x), y t+1 Acquiring people flow data for the corrected software;
3-3) time series prediction
The people flow of each street is time, the people flow is predicted by adopting a 'simple sequence time average method' in a time sequence prediction method, and the calculation steps are as follows:
3-3-1) collecting time series sample data a at l' times daily over m days i′j′ I' represents a date number,i ' =1, 2, & m, j ' represents the time sequence number, j ' =1, 2, & l, and a is retrieved after correction i′j′ ;
3-3-2) calculating the arithmetic mean of the m days l' time instants
3-3-3) predictive computation
When the time series is listed by time, the people flow predictive value y at the j-th time of the predictive date m+1 day m+1 The method comprises the following steps:
3-4) using the obtained people flow data after the hardware acquired people flow and the software system acquired people flow are processed by the least square method as the result of the people flow model.
4) Constructing a novel infectious disease travel transmission according to a road infectious disease probability database and a comprehensive infectious disease model SEIRK and a street pedestrian flow model
As shown in fig. 7, the probability calculation model specifically includes the following steps:
4-1) determination of model parameters
W: acquiring the safety distance between people according to the investigation;
x: pedestrian position is based on the nearest patient residence distance;
t: predicting the occurrence time of the infectious disease on the same day;
pp: acquiring road infection probability data from the travel infection probability database;
i (t): the number of people infected on the t-th day obtained from the new infectious disease model SEIRK;
k (t): preventable persons on day t obtained from the new infectious disease model SEIRK;
e (t): the number of people with latent people on the t-th day obtained from a new infectious disease model SEIRK;
n: and obtaining the pedestrian traffic at the current moment according to the street pedestrian traffic model.
4-2) dynamic and static combination calculation of the probability of travel infection of novel coronaviruses
If x > w, the calculation formula is as follows:
P=(0.5*(2-w)+1)*P′
otherwise, the calculation formula is as follows:
preferably, in step three, r=20, b=0.03, a=0.1, r2=30, b2=0.03, y=0.1.
Preferably, at step tetran=20, the probability of infection with infectious disease is highest.
Preferably, step five sets w to 2 meters.
The previous description of the embodiments is provided to facilitate a person of ordinary skill in the art in order to make and use the present invention. It will be apparent to those skilled in the art that various modifications can be readily made to these embodiments and the generic principles described herein may be applied to other embodiments without the use of the inventive faculty. Therefore, the present invention is not limited to the above-described embodiments, and those skilled in the art, based on the present disclosure, should make improvements and modifications without departing from the scope of the present invention.
Claims (3)
1. The travel infection probability calculation method suitable for the novel coronavirus comprises the following steps of:
1) Establishing a street infection probability database of the selected area;
2) Constructing a novel infectious disease model SEIRK based on the SEIR model according to the acquired data to simulate the infectious process of the novel infectious disease;
the step 2) comprises the following specific steps:
2-1) determining model considerations:
2-1-1) S-a susceptible person, i.e., a population susceptible to infection but not yet suffering from disease;
2-1-2) E-viral latency carrier, the latency being the group of people with disease latency who are ill but asymptomatic, the novel coronavirus also being infectious during latency;
2-1-3) I-the infected person, the infected person is the infected person;
2-1-4) R is a rehabilitating person, the rehabilitating person heals the patient, and the model considers that the rehabilitating person cannot be infected again;
2-1-5) K, which is a resistent, namely the number of people wearing the mask in the crowd, and the model considers that the wearing of the mask in epidemic situation can inhibit the spread of viruses to a certain extent, thereby effectively reducing the risk of infecting new coronaviruses;
2-2) according to the consideration, the virus diffusion simulation is carried out by taking time t as a variable, and the formula is calculated as follows:
2-2-1) counting the number of susceptible people on day t
2-2-2) counting the number of people with latent form on day t
2-2-3) counting the number of infected persons on day t
2-2-4) calculating the number of rehabilitative persons on day t
2-2-5) calculating the number of preventable persons on day t
Wherein S (t) is the number of susceptible persons on the t th day, R (t) is the number of recovered persons on the t th day, I (t) is the number of infected persons on the t th day, E (t) is the number of latent persons on the t th day, R is the number of persons contacting the patient when no strict measures are taken in the initial stage, R2 is the number of persons contacting the patient when strict measures are taken after epidemic situation occurs, b is the infection probability when no strict measures are taken in the initial stage, b2 is the infection probability when strict measures are taken after epidemic situation occurs, N is the current time pedestrian flow obtained according to a street pedestrian flow model, y is the recovery probability, and c is the probability of persons with resistance in the crowd, namely the probability of persons wearing masks in the crowd;
3) Constructing a street pedestrian flow model according to the acquired data;
4) According to the travel infection probability database, a comprehensive infection model SEIRK and a street pedestrian flow model are combined to construct a novel infectious disease travel infection probability calculation model;
the specific steps of the step 4) are as follows:
4-1) determination of model parameters
W: acquiring the safety distance between people according to the investigation;
x: pedestrian position is based on the nearest patient residence distance;
t: predicting the occurrence time of the infectious disease on the same day;
pp: acquiring road infection probability data from the travel infection probability database;
i (t): the number of people infected on the t-th day obtained from the new infectious disease model SEIRK;
k (t): the number of preventable persons on the t-th day obtained from the new infectious disease model SEIRK;
e (t): the number of people with latent people on the t-th day obtained from a new infectious disease model SEIRK;
n: acquiring pedestrian traffic at the current moment according to the street pedestrian traffic model;
4-2) dynamic and static combination calculation of the probability of travel infection of novel coronaviruses
If x > w, the calculation formula is as follows:
P=(0.5*(2-w)+1)*P′
otherwise, the calculation formula is as follows:
2. the method for calculating the traveling infection probability applicable to the novel coronavirus according to claim 1, wherein the specific steps of the step 1) are as follows:
1-1) survey the number of patients in the region to be surveyed and the geographical position of the residence thereof according to the time line;
1-2) drawing a distance-infection rate simulation curve according to the data obtained by investigation;
1-2-1) acquisition of patient affected areas
Taking the living place of the patient in the data as the circle center, and taking the infection inner diameter r as the radius to form a circular area: s=pi r 2 The region is the affected region of the patient;
1-2-2) construction of street infection Rate model
Investigation of the patient affected area as a subject has found that the risk of infection with new coronavirus in the area within 1 meter of the place where the infected person resides is reduced by a factor of 80%; the risk of infection with new coronavirus in the area 1 m or more from the place where the infected person resides is reduced by 50% by index, so that a calculation formula of a street infection rate model (x is distance (unit: m), and Y is infection rate) is obtained as follows:
1-2-3) distance-infection rate simulation curves were plotted
Drawing a distance-infection rate simulation curve by taking the distance as an x axis and the infection rate as a y axis;
1-3) establishing a road infection probability database according to the data obtained by the distance-infection rate simulation curve, wherein if the influence areas of the patients coincide, the influence values are accumulated.
3. The method for calculating the traveling infection probability applicable to the novel coronavirus according to claim 1, wherein the step 3) specifically comprises the following steps:
3-1) obtaining street traffic
3-1-1) acquiring the number of user locations
Firstly dividing a street into n areas, and recording the range O of each area i I=1, 2, …, n, after the user allows to obtain the longitude and latitude coordinates of the user, checking which area the user coordinates are in at the moment, and obtaining the people flow data of each area of the street;
3-1-2) hardware to obtain street traffic
Transmitting the shot video pictures to a server by using the existing cameras of each traffic road section, obtaining the number of pedestrians through a set of embedded intelligent video screen analysis algorithm, and carrying out density conversion according to the number of pedestrians obtained through recognition;
the density conversion formula is as follows:
N=m*c
m=p/m 0
wherein N is the total population number obtained by conversion, m is the number of pedestrians in unit length, namely the density of pedestrians, p is the average number of identifiable persons entering all cameras on the current street, and m 0 The street length taken by the camera is c is the real street length, and finally the people flow number is calculated;
3-2) correcting the people flow data with a linear least squares method
3-2-1) software and hardware data combining
The people flow data obtained by software and the people flow data obtained by hardware are combined into a group of two-dimensional data (x t ,y t ),t=1,2,…,l,x i The time t is represented, the people flow data counted by the user positioning is obtained by software, and x is t Different from each other, y t The time t is represented by the people flow data acquired by the hardware equipment, and the time l is represented by the time of acquiring the last group of data; using least square method energyObtaining a function y=f (x) to enable the sum of the least square method result and the square difference of the flow of people obtained by the actual hardware in the recorded momentMinimum;
3-2-2) the people flow data x obtained by the software at the time t+1 t+1 Substituting f (x), y t+1 Acquiring people flow data for the corrected software;
3-3) time series prediction
The people flow of each street is time, the people flow is predicted by adopting a 'simple sequence time average method' in a time sequence prediction method, and the calculation steps are as follows:
3-3-1) collecting time series sample data a at l' times daily over m days i′j′ Where i 'denotes the date number, i' =1, 2, …, m, j 'denotes the time number, j' =1, 2, …, l, and a is recovered after correction i′j′ ;
3-3-2) calculating the arithmetic mean of the m days l' time instants
3-3-3) predictive computation
When the time series is listed by time, the people flow predictive value y at the j-th time of the predictive date m+1 day m+1,j The method comprises the following steps:
3-4) using the obtained people flow data after the hardware acquired people flow and the software system acquired people flow are processed by the least square method as the result of the people flow model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110317372.7A CN113077907B (en) | 2021-03-25 | 2021-03-25 | Novel infectious disease road infection probability calculation model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110317372.7A CN113077907B (en) | 2021-03-25 | 2021-03-25 | Novel infectious disease road infection probability calculation model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113077907A CN113077907A (en) | 2021-07-06 |
CN113077907B true CN113077907B (en) | 2024-03-19 |
Family
ID=76610737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110317372.7A Active CN113077907B (en) | 2021-03-25 | 2021-03-25 | Novel infectious disease road infection probability calculation model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113077907B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113693568A (en) * | 2021-10-26 | 2021-11-26 | 广州市广美电子科技有限公司 | Infectious disease prevention challenge device |
CN114512245B (en) * | 2022-02-17 | 2022-11-15 | 四川省疾病预防控制中心 | New infectious disease model establishing method and computer device |
CN115223728A (en) * | 2022-09-21 | 2022-10-21 | 南京汉卫公共卫生研究院有限公司 | Infectious disease prediction method and system based on big data |
CN117290200B (en) * | 2023-11-23 | 2024-01-23 | 奇点数联(北京)科技有限公司 | Abnormality degree analysis system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101960504B1 (en) * | 2017-12-18 | 2019-07-15 | 세종대학교산학협력단 | Apparatus and method for epidemic spread prediction modeling |
CN111430041A (en) * | 2020-03-26 | 2020-07-17 | 北京懿医云科技有限公司 | Infectious disease epidemic situation prediction method and device, storage medium and electronic equipment |
CN112164470A (en) * | 2020-08-28 | 2021-01-01 | 南京星耀智能科技有限公司 | Event-driven infectious disease prediction method |
CN112382405A (en) * | 2020-11-02 | 2021-02-19 | 北京机电工程研究所 | Method, device, equipment and medium for predicting and evaluating structured new coronary pneumonia epidemic situation |
CN112529329A (en) * | 2020-12-21 | 2021-03-19 | 广东珠江智联信息科技股份有限公司 | Infectious disease prediction method based on BP algorithm and SEIR model |
CN112542250A (en) * | 2020-11-04 | 2021-03-23 | 温州大学 | Global new coronavirus transmission prediction method based on optimized SEIRD model |
-
2021
- 2021-03-25 CN CN202110317372.7A patent/CN113077907B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101960504B1 (en) * | 2017-12-18 | 2019-07-15 | 세종대학교산학협력단 | Apparatus and method for epidemic spread prediction modeling |
CN111430041A (en) * | 2020-03-26 | 2020-07-17 | 北京懿医云科技有限公司 | Infectious disease epidemic situation prediction method and device, storage medium and electronic equipment |
CN112164470A (en) * | 2020-08-28 | 2021-01-01 | 南京星耀智能科技有限公司 | Event-driven infectious disease prediction method |
CN112382405A (en) * | 2020-11-02 | 2021-02-19 | 北京机电工程研究所 | Method, device, equipment and medium for predicting and evaluating structured new coronary pneumonia epidemic situation |
CN112542250A (en) * | 2020-11-04 | 2021-03-23 | 温州大学 | Global new coronavirus transmission prediction method based on optimized SEIRD model |
CN112529329A (en) * | 2020-12-21 | 2021-03-19 | 广东珠江智联信息科技股份有限公司 | Infectious disease prediction method based on BP algorithm and SEIR model |
Non-Patent Citations (2)
Title |
---|
基于SEIR模型的新型冠状病毒疫情实证分析;陈柯宇;;统计与管理(06);全文 * |
基于SEIR模型视角的重大公共卫生事件中伪科学网络谣言的传播治理:以新冠肺炎疫情为例;陈安;王子君;陈樱花;;科技导报(04);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113077907A (en) | 2021-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113077907B (en) | Novel infectious disease road infection probability calculation model | |
CN111027525B (en) | Method, device and system for tracking potential infected persons in public places during epidemic situation | |
Li et al. | Multilevel modelling of built environment characteristics related to neighbourhood walking activity in older adults | |
CN113192645B (en) | Discrete grid-based infectious disease space-time diffusion evolution and artificial prevention and control simulation method | |
Karaman et al. | Development of smart camera systems based on artificial intelligence network for social distance detection to fight against COVID-19 | |
CN113784292B (en) | Urban epidemic situation risk assessment method based on mobile phone signaling data | |
Ohkusa et al. | Real-time estimation and prediction for pandemic A/H1N1 (2009) in Japan | |
CN109668563A (en) | Processing method and processing device based on indoor track | |
CN115394455A (en) | Infectious disease time-space diffusion prediction method and device based on spatial clustering discrete grid | |
Doorley et al. | Mobility and COVID-19 in Andorra: Country-scale analysis of high-resolution mobility patterns and infection spread | |
Harnik et al. | Parks: How far is too far | |
CN113096785B (en) | Epidemic situation period medical resource allocation method | |
US20230170096A1 (en) | Systems, devices, and methods for managing contact instances of persons of interest | |
de Montigny et al. | A spatial analysis of the physical and social environmental correlates of discarded needles | |
Berke et al. | Using mobile phone data to estimate dynamic population changes and improve the understanding of a pandemic: A case study in Andorra | |
CN116453696B (en) | Respiratory tract disease infection risk prediction method based on personnel space-time distribution model | |
CN110555923A (en) | Personnel registration method and system | |
CN112259248A (en) | COVID-19 overseas input risk assessment and prediction method | |
Choi et al. | A face mask wearing detection technique to protect against coronavirus in metro stations | |
CN108596473B (en) | Drainage basin pollution management system and new accounting method | |
CN114708985B (en) | Respiratory infectious disease early warning method and system | |
Kaur et al. | Quantifying personal exposure to spatio-temporally distributed air pollutants using mobile sensors | |
CN112216144B (en) | Big data line-moving planning analysis method based on hospital parking lot | |
Guerrieri et al. | Social distance evaluation in transportation systems and other public spaces using deep learning | |
Goto et al. | Human vulnerability index for evaluating tsunami evacuation capability of communities |
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 |