1. A method for protecting and distributing explosives based on environmental gridding is characterized by comprising the following steps:

determining the environmental protection weight of each position point to be protected in the area to be protected;

constructing a protection capability model of the protection equipment, and determining a protection capability factor of each position point to be protected;

determining an outsourcing rectangle of the region to be protected, and carrying out grid division on the outsourcing rectangle by adopting unit grids; determining an effective factor corresponding to each unit grid according to whether the unit grid is in an effective area;

determining the normalized protection capability of the position point to be protected in the unit grid in which the position point to be protected is positioned according to the unit grid in which the position point to be protected is positioned and by combining the environment protection weight and the protection capability factor;

aiming at an outsourcing rectangle of an area to be protected, selecting unit grids at four vertex positions, and calculating the protection capability of any position point in the area to be protected by utilizing a bilinear interpolation mode according to the selected four unit grids; and laying protection equipment according to the calculated protection capability of any position point in the region to be protected.

2. The method according to claim 1, wherein the determining the environmental protection weight for each location point to be protected in the area to be protected specifically comprises:

w_f＝F_S/(F_A×Ar)

wherein w_fThe environmental protection weight of the position point to be protected is obtained; f_SThe number of times of the dangerous explosive materials appearing in unit time at the position point to be protected is shown; f_AThe number of times of the dangerous explosive materials in the whole area to be protected in unit time is counted; ar is the area of the position point to be protected in meters when the number of times of the dangerous explosive article in unit time of the position point to be protected is counted; wherein the unit time is selected to be one year or more.

3. The method according to claim 2, wherein the building of the protection capability model of the protection device and the determining of the protection capability factor of each location point to be protected are specifically:

wherein pij is a protective ability factor; pd represents the location of the flexible protective equipment setting; cij represents an index of the position point to be protected based on grid division; pd-cij represents the distance between the flexible protective equipment and the grid where the position point to be protected is located, and the unit is meter.

4. The method of claim 3, wherein the outsourcing rectangle is determined for the area to be protected and is gridded using the cell grid; determining an effective factor corresponding to each unit grid according to whether the unit grid is in the effective area, specifically:

wherein vij is an effective factor corresponding to the unit grid under the index of i and j; i is a grid mark of the unit grid in the horizontal direction; j is a grid mark of the unit grid in the vertical direction; u is an effective area;

if the 50% area and above of the unit grid are in the effective area, the effective factor of the unit grid is 1, otherwise, the effective factor is 0.

5. The method according to claim 4, wherein the determining, according to the cell grid where the position point to be protected is located, the normalized protection capability of the position point to be protected in the cell grid where the position point to be protected is located in combination with the environment protection weight and the protection capability factor, specifically comprises:

wherein F (i, j) represents the protection capability of the cell grid under the i, j index after normalization; f_maxCalculated vijxw within a grid representing the entire area to be protected_f/p_ijIs measured.

6. The method according to claim 5, wherein for the outsourcing rectangle of the area to be protected, a unit grid at four vertex positions is selected, and the protection capability of any position point in the area to be protected is calculated by using a bilinear interpolation mode according to the selected four unit grids, specifically:

the outsourcing rectangle for the area to be protected takes the lower right corner as an origin, the length of the outsourcing rectangle is an x axis, and the width of the outsourcing rectangle is a y axis, and a coordinate system is established;

selecting unit grids at four vertex positions, wherein the unit grids comprise a lower left corner grid c1, a lower right corner grid c2, an upper left corner grid c3 and an upper right corner grid c 4;

wherein the protective capacities of c 1-c 4 are known values, namely F (c1) to F (c 4);

the grid where the position point to be obtained is located is T, the grid projected by the T on the edge grid line where the c1 and the c2 are located together is c12, and the grid projected by the T on the edge grid line where the c3 and the c4 are located together is c 34; the protective capacities of c12 and c34 are F (c12) and F (c34), respectively

Wherein x10 is the x-axis coordinate corresponding to the center point of grid c 2; x0 is the x-axis coordinate corresponding to the center point of grid c 1; xi is the x-axis coordinate corresponding to the central point of the grid c 12;

formula for calculating protection capability f (T) of grid T:

wherein y10 is the y-axis coordinate corresponding to the center point of grid c 34; y0 is the y-axis coordinate corresponding to the center point of grid c 12; and yi is the coordinate of the central point of the grid T corresponding to the y axis.

7. The method of any of claims 1 to 6, further comprising the steps of:

evaluating the protection capability of the to-be-protected area with deployed protection equipment, wherein the evaluation includes a protection capability index and a protection equalization index;

the protective capacity index is prc:

wherein n is_iThe number of grids in the effective area, the protection capability of which is smaller than a set empirical value mu, is determined; n is_aRepresenting the number of all meshes within the active area;

the protection equalization index is eq:

the larger the equalization degree, the better

Wherein cij is all grids in the active area;is the average of the protective capacities of all grids in the effective areaA value; n is the number of all meshes in the active area.

Background

The current explosive protection equipment is usually placed in the environment according to the principle of personal experience and even where the explosive protection equipment is conveniently placed, and no scientific basis can be referred to. The evaluation of the hazard of the explosives usually occurs after the explosion, and the occurrence situation of the explosion is inverted according to the damage situation of the scene, which belongs to a post analysis summary behavior. There is currently a lack of pre-analysis for protection against explosives.

For the sudden protection work, the nearest protection equipment cannot be effectively organized to protect. On the basis of the networking of dangerous explosive equipment, the application and protection of the equipment can be more scientific, and the distribution of protective equipment is reasonably arranged according to the protection range of the equipment and the protection grade requirements of different positions of the airport environment.

The protection range of the explosive device is taken as an example of a G500 flexible explosion-proof barrel, and G500 is flexible protection equipment with the protection equivalent of 500G of TNT. G500 is about 60kg, and can be used independently or matched with a soft sanitary vehicle. The concept of protection range is basically not provided because the prior steel explosion-proof barrel/tank is too heavy to be moved. The protection scope is a new concept proposed after adopting a novel flexible protection device. The protection range refers to the range which can be reached by the protection equipment within 1 minute and is centered on the protection equipment, namely the range which can be protected within 1 minute.

For large public places, for example, airport environments include many places such as security checkpoints, counseling desks, toilets, shopping malls, restaurants, ticket offices, packing areas, customs, health stations, public rest areas, drinking water offices, public telephone offices, and the like. Through manual analysis, the difficulty is very high along with the improvement of the environmental complexity, and sometimes effective support is difficult to provide reasonably without depending on scientific means. With the development of protection technology in China, the protection equipment for the explosive substances is greatly improved, so that the effective utilization of the equipment has urgent needs and great significance. Meanwhile, with the development of network technology and the gradual application of the internet of things, the technology of the internet of things and big data is adopted in many industrial fields, and the layout method provides a feasible theoretical method for protecting networking and datamation.

Therefore, an efficient, scientific, reasonable and rapid scheme capable of optimizing the layout of protective equipment in a large public place is absent at present.

Disclosure of Invention

In view of the above, the invention provides an environmental gridding-based explosive protection layout method, which can efficiently and scientifically arrange protection equipment through advanced simulation calculation, so as to achieve the purposes of no weakness of full environmental protection and most reasonable and rapid disposal time.

In order to achieve the purpose, the technical scheme of the invention comprises the following steps:

and determining the environment protection weight of each position point to be protected in the area to be protected.

And constructing a protection capability model of the protection equipment, and determining a protection capability factor of each position point to be protected.

Determining an outer rectangular shape of the region to be protected, and performing grid division on the outer rectangular shape by adopting unit grids; and determining the effective factor corresponding to each unit grid according to whether the unit grid is in the effective area.

And determining the normalized protection capability of the position point to be protected in the unit grid according to the unit grid where the position point to be protected is located and by combining the environment protection weight and the protection capability factor.

Aiming at the outsourcing rectangle of the area to be protected, selecting unit grids at the positions of four vertexes, and calculating the protection capability of any position point in the area to be protected by utilizing a bilinear interpolation mode according to the selected four unit grids; and laying protection equipment according to the calculated protection capability of any position point in the region to be protected.

Further, determining an environmental protection weight for each to-be-protected location point in the to-be-protected area, specifically:

w_f＝F_S/(F_A×Ar)

wherein w_fThe environmental protection weight of the position point to be protected is obtained; f_SThe number of times of the dangerous explosive materials appearing in unit time at the position point to be protected is shown; f_AThe number of times of the dangerous explosive materials in the whole area to be protected in unit time is counted; ar is the area of the position point to be protected in meters when the number of times of the explosive substances in unit time of the position point to be protected is counted; wherein the unit time is selected to be one year or more.

Further, a protection capability model of the protection device is built, and a protection capability factor of each to-be-protected location point is determined, specifically:

wherein pij is a protective ability factor; pd represents the location of the flexible protective equipment setting; cij represents an index of the position point to be protected based on grid division; pd-cij represents the distance between the flexible protective equipment and the grid where the position point to be protected is located, and the unit is meter.

Further, determining an outer rectangular shape of the region to be protected, and performing grid division on the outer rectangular shape by adopting unit grids; determining an effective factor corresponding to each unit grid according to whether the unit grid is in the effective area, specifically:

wherein vij is an effective factor corresponding to the unit grid under the index of i and j; i is a grid mark of the unit grid in the horizontal direction; j is a grid mark of the unit grid in the vertical direction; u is the effective area.

If the 50% area and above of the unit grid are in the effective area, the effective factor of the unit grid is 1, otherwise, the effective factor is 0.

Further, according to the unit grid where the position point to be protected is located, the normalized protection capability of the position point to be protected in the unit grid where the position point to be protected is located is determined by combining the environmental protection weight and the protection capability factor, and specifically, the normalized protection capability is as follows:

wherein F (i, j) represents the protection capability of the cell grid under the i, j index after normalization; f_maxCalculated vijxw within a grid representing the entire area to be protected_f/p_ijIs measured.

Further, aiming at the outsourcing rectangle of the area to be protected, unit grids at four vertex positions are selected, and according to the selected four unit grids, the protection capability of any position point in the area to be protected is calculated by using a bilinear interpolation mode, specifically:

and aiming at the outsourcing rectangle of the area to be protected, a coordinate system is established by taking the lower right corner as an origin, and the length of the outsourcing rectangle is an x axis and the width of the outsourcing rectangle is a y axis.

Cell grids at four vertex positions are selected, including a lower left corner grid c1, a lower right corner grid c2, an upper left corner grid c3, and an upper right corner grid c 4.

Wherein the protective capacities of c 1-c 4 are known values, namely F (c1) to F (c 4).

The grid where the position point to be obtained is located is T, the grid projected by the T on the edge grid line where the c1 and the c2 are located together is c12, and the grid projected by the T on the edge grid line where the c3 and the c4 are located together is c 34; the protective capacities of c12 and c34 are F (c12) and F (c34), respectively

Wherein x10 is the x-axis coordinate corresponding to the center point of grid c 2; x0 is the x-axis coordinate corresponding to the center point of grid c 1; xi is the x-axis coordinate corresponding to the center point of grid c 12.

Formula for calculating protection capability f (T) of grid T:

wherein y10 is the y-axis coordinate corresponding to the center point of grid c 34; y0 is the y-axis coordinate corresponding to the center point of grid c 12; and yi is the coordinate of the central point of the grid T corresponding to the y axis.

Further, the method also comprises the following steps:

and evaluating the protection capability of the to-be-protected area with the deployed protection equipment, wherein the evaluation includes a protection capability index and a protection equalization index.

The protective capacity index is prc:

wherein n is_iThe number of grids in the effective area, the protection capability of which is smaller than a set empirical value mu, is determined; n is_aRepresenting the number of all meshes in the active area.

The protection equalization index is eq:

the larger the equalization degree, the better

Wherein cij is all grids in the active area;the average value of the protection capability of all grids in the effective area is obtained; n is the number of all meshes in the active area.

Has the advantages that:

according to the environment gridding based explosive protection layout method, the protection capability model of the protection equipment is constructed by setting the environment protection weight, providing the effective protection factor based on the protection capability and the protection distance of the protection equipment, and the protection parameters are quantized according to the environment by adopting a gridding layout distribution method according to the environment protection weight and the protection capability model, so that the scientific layout of the explosive protection equipment is further realized, and the optimization of protection and explosive disposal is achieved. The layout method is efficient and scientific in arranging the protective equipment, achieves the purposes of no weakness of full-environment protection and most reasonable and rapid disposal time, is suitable for simplified scenes, is also suitable for complex scenes, and has more obvious advantages compared with the traditional manual experience layout protective force along with the improvement of the complexity of the scenes.

Drawings

FIG. 1 is a plan view of an airport for layout in an embodiment of the present invention;

FIG. 2 is a diagram illustrating the distribution of guard weights in an embodiment of the present invention;

FIG. 3 is a schematic diagram of the protective capabilities of protective equipment in an embodiment of the invention;

FIG. 4 is a schematic diagram of a protective gridding partition method according to an embodiment of the present invention;

fig. 5 is a schematic diagram of layout grid refinement solution based on bilinear interpolation in the embodiment of the present invention.

Fig. 6 is a flowchart of a method for protecting a hazardous explosive material based on environmental gridding according to the present invention.

Detailed Description

The invention is described in detail below by way of example with reference to the accompanying drawings.

The embodiment of the invention takes an airport as an example, and a plan view of the airport is shown in figure 1, wherein the scene layout adopts a simplified description mode, and public facilities such as toilets, service desks, drinking water places and the like are hidden for convenient expression. The layout method is not only suitable for simplified scenes, but also suitable for complex scenes, and analysis shows that the layout method has more obvious advantages compared with the traditional method for laying out protection force through manual experience along with the improvement of the complexity of the scenes.

And S1, determining the environmental protection weight of each position point to be protected in the area to be protected.

An environmental protection weight is a probability that a threat explosive will appear in a certain place in the whole environment for a certain period of time.

As shown in fig. 2, a proposed algorithm for guard weights is proposed here:

w_f＝F_S/(F_A×Ar)

wherein w_fThe environmental protection weight of the position point to be protected is obtained; f_SThe number of times of the dangerous explosive materials appearing in unit time at the position point to be protected is shown; f_AThe number of times of the dangerous explosive materials in the whole area to be protected in unit time is counted; ar is the area of the position point to be protected in meters when the number of times of the explosive substances in unit time of the position point to be protected is counted; wherein the unit time is selected for one year or more, if the number of times is small, it is considered to increase the length of the unit time, for example, two years, five years, etc.

S2, constructing a protection capability model of the protection equipment, and determining a protection capability factor of each position point to be protected

The protection capability model of the protection device is described by taking the flexible explosion-proof barrel G500 as an example. G500 is the novel explosion-proof equipment that adopts flexible composite protection technology, possesses technical characterstic such as lightweight, flexible protection, fire prevention. The technical parameters are as follows:

height: 613mm

Outer diameter: 550mm

Total weight: 62kg of

Protection capability: 500g TNT +200g fragment

The method is more critical in the initial stage after finding out the explosive, and can effectively reduce the damage risk by protecting the explosive as soon as possible and reducing the exposure time of the explosive. 200 seconds are set as key time according to empirical values, the moving speed of the flexible anti-explosion barrel is 2 meters per second, namely the main action radius of the protection of the flexible anti-explosion barrel is 400 meters multiplied by the two.

A schematic diagram of the protection capabilities is shown in fig. 3.

And taking the maximum value of the protection capacity as the protection capacity of the point in the area with overlapped protection capacity (the area covered by the protection capacity of two or more explosion-proof barrels).

Wherein pij is a protective ability factor; pd represents the location of the flexible protective equipment setting; cij represents an index of the position point to be protected based on grid division; pd-cij represents the distance between the flexible protective equipment and the grid where the position point to be protected is located, and the unit is meter.

S3, determining an outer rectangular shape of the region to be protected, and performing grid division on the outer rectangular shape by adopting unit grids; and determining the effective factor corresponding to each unit grid according to whether the unit grid is in the effective area.

In the embodiment of the present invention, in order to analyze the overall protection layout in a quantitative manner, a method for distributing a protection layout based on grid division is provided, as shown in fig. 4. The grid requirement can cover the whole environment area by taking the lower left corner as the origin of coordinates, the horizontal axis as the X axis and the vertical axis as the Y axis. The interval value of the unit cell is determined according to the size of the whole area, the length of the example airport is 800 meters, the width of the airport is 240 meters, and the interval width of the unit cell is 20 meters.

Since the environment is generally an irregular rectangle, and therefore contains both valid regions (inside the airport) and invalid regions (outside the airport), the validity factor formula is as follows:

wherein vij is an effective factor corresponding to the unit grid under the index of i and j; i is a grid mark of the unit grid in the horizontal direction; j is a grid mark of the unit grid in the vertical direction; u is an effective area;

if the 50% area and above of the unit grid are in the effective area, the effective factor of the unit grid is 1, otherwise, the effective factor is 0.

The normalized protection formula in a grid can be expressed as:

wherein F (i, j) represents the protection capability of the cell grid under the i, j index after normalization; f_maxCalculated vijxw within a grid representing the entire area to be protected_f/p_ijIs measured.

And S4, aiming at the outsourcing rectangle of the area to be protected, selecting unit grids at the positions of four vertexes, and calculating the protection capability of any position point in the area to be protected by utilizing a bilinear interpolation mode according to the selected four unit grids.

In a case where the environment is relatively large, the number of the entire grid is relatively large, and it is difficult to perform fine-grained division for reducing the workload, for example, the width of the unit grid is 20 meters. Because the size of explosion-proof equipment is less, consider simultaneously that the location demand based on the thing networking, the positional information of the acquisition equipment that needs more meticulous to and the overall arrangement condition of protection strength. The current precision generally cannot meet the field requirement and needs to be refined.

Referring to the size of the explosion-proof equipment, the explosion-proof equipment is generally in the range of 1 x1 meter, the explosion-proof equipment is in the range of 2 x 2 meters, and the common minimum equipment of other airports is generally in the range of 1-2 meters, so the minimum grid unit is set to be 2 x 2 meters.

And (3) selecting a bilinear interpolation algorithm to finish the refinement of the grid according to the characteristics of the protection layout grid, wherein the specific implementation method of the bilinear interpolation is as follows:

aiming at an outsourcing rectangle of an area to be protected, a coordinate system is established by taking the lower right corner as an origin, and the length of the outsourcing rectangle is an x axis and the width of the outsourcing rectangle is a y axis;

selecting unit grids at four vertex positions, wherein the unit grids comprise a lower left corner grid c1, a lower right corner grid c2, an upper left corner grid c3 and an upper right corner grid c 4;

wherein the protective capacities of c 1-c 4 are known values, namely F (c1) to F (c 4);

the grid where the position point to be obtained is located is T, the grid projected by the T on the edge grid line where the c1 and the c2 are located together is c12, and the grid projected by the T on the edge grid line where the c3 and the c4 are located together is c 34; the protective capacities of c12 and c34 are F (c12) and F (c34), respectively

As shown in fig. 5, the diamond grid point c34 is obtained according to the triangles c3 and c4, the diamond grid point c12 is obtained according to the triangle grid points c1 and c2, and the five-pointed star grid point is obtained according to the diamond grid points c34 and c12, and the calculation formula is as follows:

wherein x10 is the x-axis coordinate corresponding to the center point of grid c 2; x0 is the x-axis coordinate corresponding to the center point of grid c 1; xi is the x-axis coordinate corresponding to the central point of the grid c 12;

formula for calculating protection capability f (T) of grid T:

wherein y10 is the y-axis coordinate corresponding to the center point of grid c 34; y0 is the y-axis coordinate corresponding to the center point of grid c 12; and yi is the coordinate of the central point of the grid T corresponding to the y axis.

For the special case of the edge point in fig. 5, only the single linear interpolation is adopted, and the principle is the same as above.

The embodiment of the invention also comprises the step of evaluating the protection capability of the to-be-protected area where the protection equipment is deployed, wherein the evaluation comprises a protection capability index and a protection equalization index.

The effective protection evaluation index represents whether the current protection layout can perform effective protection of the whole scene.

The protection equalization evaluation index represents whether the layout of protection force in the current protection layout is reasonable or not.

Protection capability characterization formula:

wherein: prc is effective protection evaluation index; ni: in the effective area, the number of grids with F (i, j) smaller than a certain empirical value mu, and the value of mu needs to be obtained according to the test. And na: the number of all meshes of the active area.

The value range of the effective protection evaluation index is 0-1, and the smaller the value is, the better the protection performance is.

Protection equalization index representation formula:

wherein: eq: protection equalization evaluation indexes; cij: all meshes within the active area;the average value of protection; n is the number of all grids in the effective area; the value range of the protection equalization index is 0-positive infinity, and the smaller the value is, the better the protection performance is.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.