Environmental risk assessment method for dangerous goods container yard in ports

Abstract

The present invention relates to the field of risk assessment technology, and specifically to a method for assessing the environmental risk of a port dangerous goods container yard, including a dangerous goods container storage statistics module, a risk source positioning module, an environmental analysis module, and a risk comprehensive assessment module. The present invention can calculate a risk assessment value through multiple analysis methods, and the risk assessment value is affected by multiple factors such as the degree of dangerous and harmful characteristics of dangerous goods in the dangerous goods container stack, the storage order, and the distribution of personnel and key equipment and facilities in the surrounding area of ​​the yard. The calculated risk assessment value can be combined with multiple factors to comprehensively assess the degree of danger and harm after the dangerous goods have an accident to generate flammable, toxic and harmful gases. The assessment result is more efficient and more comprehensive than the assessment method in the prior art, which is helpful to predict and judge the safety impact of the on-site operating environment caused by the accident after the dangerous situation occurs.

Description

Technical Field

[0001] The invention relates to the technical field of risk assessment, and in particular to an environmental risk assessment method for a port dangerous goods container yard. Background Art

[0002] A port dangerous goods container yard is a special place for storing dangerous goods containers in the open air within the port area. It is mainly engaged in the lifting, stacking and storage of dangerous goods containers. The scale of a port dangerous goods container yard (container capacity, number of ground containers, yard area, etc.) is usually determined based on the port's dangerous goods container throughput, category, average storage period, etc. Dangerous goods containers are usually stored in different storage areas of the yard according to the nature and category of the dangerous goods.

[0003] Dangerous goods container refers to a container loaded with explosive, flammable, toxic and other dangerous goods. Dangerous goods containers loaded with different types of dangerous goods need to be stacked separately. In the prior art, Chinese Patent Publication No. CN117610939B discloses a method and system for assessing the environmental risk of dangerous goods container yards in ports. By establishing a new pollution factor diffusion model, the concentration of pollution factors at different locations is calculated, thereby assessing the risk of groundwater environment being affected when dangerous goods leak and optimizing the assessment efficiency. However, it does not assess the safety impact of toxic gases after the leakage of dangerous goods containers on the surrounding personnel. If the toxic gases generated by the leakage of dangerous goods are not assessed in time and emergency disposal measures are not taken, they may further cause poisoning accidents and cause casualties. Therefore, after the dangerous goods container leaks, the consequences of the accident are not reduced or mitigated, and the safety risks of major and serious accidents in port operations are prevented and resolved. It is necessary to conduct a risk assessment on the safety impact of toxic gases on the surrounding personnel. Summary of the invention

[0004] In view of the above-mentioned shortcomings of the prior art, the present invention provides a port dangerous goods container yard environmental risk assessment method, which can effectively solve the problem in the prior art that it is difficult to conduct risk prediction and assessment of the possible poisoning accidents caused by dangerous goods leakage in dangerous goods containers.

[0005] To achieve the above objectives, the present invention is implemented through the following technical solutions:

[0006] The present invention provides a method for assessing the environmental risk of a port dangerous goods container yard, which at least comprises: a dangerous goods container storage statistics module, a risk source positioning module, an environmental analysis module, and a risk comprehensive assessment module, wherein the yard range is recorded as a storage area, and different types of dangerous goods containers are stacked separately to form a dangerous goods container pile;

[0007] Dangerous goods container storage statistics module, obtain the type of each dangerous goods container stack in the storage area and record it as the dangerous type, obtain the dangerous type of each dangerous goods container stack and the loading volume of each container and calculate its dangerous unit value DE, obtain the number of containers in each dangerous goods container stack and record it as the stack number n, substitute it into the formula R danger =n*DE to calculate the hazard diffusion parameter R of each dangerous goods container pile danger ;

[0008] A pile of dangerous goods containers that generates toxic gases due to an accident shall be recorded as an accident pile;

[0009] The risk source location module obtains the specific location of the accident pile in the storage area, records the specific location as the danger occurrence point, and records the time when the danger occurrence point is determined as the occurrence time;

[0010] The environmental analysis module includes a ring analysis unit and a sector analysis unit. The ring analysis unit divides a certain range around the storage area into multiple ring areas Q i , and analyze each annular area, and calculate its life impact index δ based on the resident data in the annular area i The fan-shaped analysis unit monitors the wind direction angle in the storage area to construct multiple initial wind directions and corresponding fan-shaped areas, records the wind speed regularly, and calculates the life impact index δ corresponding to each fan-shaped area based on the human distribution characteristics and wind speed in each fan-shaped area. p ′ and record the duration T of each initial wind direction p ;

[0011] The comprehensive risk assessment module obtains the hazard diffusion parameter R of dangerous goods container stacks danger , each annular area Q i Life impact index δ i And the life impact index δ of different fan-shaped areas p ′, and the risk assessment value DS is obtained through comprehensive analysis.

[0012] Furthermore, the risk unit value calculation process is as follows:

[0013] Get the gas type generated by the chemical reaction corresponding to the hazardous type, set a toxicity value D, assign the toxicity value according to the acute toxicity category to which the gas type belongs. When the category to which the gas type belongs is N, set its toxicity value to 6-N. When there are multiple gas types with acute toxicity, select the gas type with the largest corresponding toxicity value as the main toxicity source, and use the toxicity value of this gas type as the toxicity value of the hazardous type. Obtain the loading volume G per box and the toxicity value D and substitute them into the formula DE=k*G*10 D Calculation is performed to obtain the danger unit value DE, where k is a preset proportional coefficient.

[0014] Furthermore, the process of dividing the annular area by the annular analysis unit is as follows:

[0015] The hazard diffusion parameter of the dangerous goods container stack corresponding to the hazard occurrence point is obtained and recorded as R′ danger , obtain the toxicity value of the dangerous goods container pile corresponding to the danger point and record it as R′ danger Substitute into the formula The maximum impact range R1 is obtained by calculation, where h is the preset height value and r is the preset proportional coefficient. A circular area is drawn with the danger point in the storage area as the center and the maximum impact range R1 as the radius. Draw a circular boundary line to divide the circular area into 5 annular areas.

[0016] Furthermore, the process of constructing the fan-shaped area by the fan-shaped analysis unit is as follows:

[0017] The wind direction and wind speed of the storage area are obtained regularly, and the wind direction at the time of occurrence is recorded as the initial wind direction. The vertex of the danger occurrence point and the initial wind direction are used as the angle bisector to construct an area angle of P°, where P is the preset value. danger , Substitute into the formula Calculation is performed to obtain the diffusion influence distance R2, and a circular outline is drawn with the vertex of the regional angle as the center and the diffusion distance R1 as the radius, and the circular outline intersects with the regional angle to form a fan-shaped area;

[0018] The angle difference between the current wind direction and the initial wind direction is calculated regularly. When the angle difference is less than or equal to the preset wind direction change threshold, the fan-shaped area corresponding to the initial wind direction is still used as the range area for calculating the life impact index. When the angle difference is greater than the preset wind direction change threshold, the current wind direction is used as the initial wind direction and the fan-shaped area is reconstructed, and the fan-shaped area is used as the range area for calculating the life impact index.

[0019] Furthermore, the life impact index calculation process is as follows:

[0020] The area to be calculated is recorded as the target area. When the target area is a plurality of annular areas divided by the annular analysis unit, they are recorded as Q1, Q2, Q3, Q4, and Q5 in descending order according to the outer diameter. The number of people NUM1 in each annular area is counted, and the total number of commercial buildings, office buildings, industrial buildings, and educational buildings in each annular area is counted as the number of buildings NUM2. The number of people NUM1 and the number of buildings NUM2 are normalized and substituted into the formula δ i =e 6-i*(η1*NUM1+η2*NUM2) to obtain each annular area Q i Life impact index δ i , where i = 1, 2, 3, 4, 5, η1, η2 are preset weight coefficients, and e is a natural constant;

[0021] When the target area is a fan-shaped area constructed by the fan-shaped analysis unit, the number of initial wind directions q within the time range from the occurrence time to the current time is obtained, and the duration T of each initial wind direction is calculated. p , P represents the sequence number of the initial wind direction, P = 1, 2, 3, ..., q, and each initial wind direction is analyzed separately:

[0022] Plot duration T p A line graph of wind speed changes over time within the range is drawn, and the area of ​​the graph enclosed by the line graph and the horizontal axis of the coordinate is calculated and recorded as the wind force impact value F. The number of people NUM1′ and the number of buildings NUM2′ in the sector area are obtained, and the formula δ′=F*(η1*NUM1′+η2*NUM2′) is substituted for calculation to obtain the life impact index δ of the sector area. p ′.

[0023] Furthermore, the comprehensive risk assessment module obtains the risk diffusion parameter R of the dangerous goods container stack corresponding to the risk occurrence point. danger , each annular area Q i Life impact index δ i And the life impact index δ of different fan-shaped areas p ′, substitute it into the formula to calculate The risk assessment value DS is obtained, where μ1, μ2, and μ3 are all preset weight coefficients.

[0024] Compared with the known prior art, the technical solution provided by the present invention has the following beneficial effects:

[0025] 1. The dangerous goods container stockpile statistics module of the present invention regards the dangerous goods container stack formed by stacking dangerous goods containers as a separate analysis object, and predicts the types of flammable gases and toxic gases generated when accidents such as leakage, fire, explosion, etc. occur according to the types of goods loaded therein, and then predicts the degree of harm that may be caused after the accident based on the degree of dangerous and harmful characteristics of the dangerous goods (toxicity intensity) and the order of magnitude of the cargo stockpile. This prediction method is based on the principles of chemical reactions and physical reactions, and can correctly and appropriately calculate and analyze the risks that the dangerous goods container stack can cause in the worst case, thereby facilitating the staff to understand the risk level of the dangerous goods container stack.

[0026] 2. The environmental analysis module of the present invention is based on the diffusion law of gas molecules in natural windless conditions and the diffusion law under the influence of wind. By calculating the life impact index in the target area, it can help staff determine the degree of impact and the danger level of flammable gases and toxic gases generated by the leakage of dangerous goods in dangerous goods container piles on human life activities in the surrounding area under the conditions without wind influence and with wind influence. It predicts the impact of dangerous goods on the air environment through mathematical means and combines the distribution of humans to further analyze the diffusion impact of toxic gases in two situations, thereby improving the efficiency of emergency response analysis after an accident occurs and facilitating staff to make reasonable judgments and decisions.

[0027] 3. The risk assessment value of the present invention is affected by multiple factors such as the degree of dangerous and harmful characteristics of dangerous goods in the dangerous goods container pile (toxicity intensity), the storage quantity level, and the personnel in the surrounding area of ​​the yard. The risk assessment value obtained by calculation can combine multiple factors to comprehensively evaluate the degree of harm after the dangerous goods have an accident and generate toxic gases. The evaluation result is more efficient and comprehensive than the evaluation method in the prior art, which helps to predict and judge the impact on the air environment caused by the dangerous situation after it occurs. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention, and for ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0029] Figure 1 It is the overall module block diagram of the present invention. DETAILED DESCRIPTION

[0030] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0031] The present invention will be further described below in conjunction with the embodiments.

[0032] Risk refers to the probability or degree of damage or destruction to a system or activity caused by a possible outcome under certain conditions. It is usually used to describe situations that may lead to loss or injury. Dangerous goods in dangerous goods containers will generate toxic and harmful gases that harm the surrounding environment when they leak or react chemically. Dumping, damage, combustion and explosion and other reasons may cause leakage or reaction of dangerous goods. Therefore, after the above events occur, timely assessment of environmental hazard risks will help managers understand the level of accidents and the extent of their impact.

[0033] See also Figure 1 The environmental risk assessment method for dangerous goods container yards in ports at least includes: dangerous goods container storage statistics module, risk source location module, environmental analysis module, and risk comprehensive assessment module.

[0034] The yard area is recorded as the storage area, and the containers are stacked in order in the storage area. Different types of dangerous goods containers are stacked separately to form different dangerous goods container piles, that is, dangerous goods containers of the same type are stacked together, while dangerous goods containers of different types are stacked separately. This practice can ensure the safety of port operations. Because different types of dangerous goods may have different physical and chemical properties, and may even react with each other to cause danger, storing them separately can avoid potential fires, explosions and other accidents, and ensure the safety of workers and the integrity of goods.

[0035] The dangerous goods container stockpile statistics module obtains the type of each dangerous goods container stack in the stockpile area as the dangerous type, and calculates its dangerous unit value according to the dangerous type of each dangerous goods container stack and the loading amount of each container (the average weight of dangerous goods loaded in each 20-foot standard container). The calculation process of the dangerous unit value is as follows:

[0036] Get the gas type generated by the chemical reaction corresponding to the hazardous type, set a toxicity value D, assign the toxicity value according to the acute toxicity category to which the gas type belongs. When the category to which the gas type belongs is N, set its toxicity value to 6-N. That is, when the category to which the gas type belongs is 5, its toxicity value is set to 1. When the category to which the gas type belongs is 4, its toxicity value is set to 2. And so on. When there are multiple gas types with acute toxicity (there may be multiple toxicity values ​​of different sizes at this time), select the gas type with the largest corresponding toxicity value as the main toxicity source, and use the toxicity value of the gas type as the toxicity value of the hazardous type. Obtain the loading volume G and toxicity value D per box and substitute them into the formula DE=k*G*10 DCalculation is performed to obtain the danger unit value DE, where k is the preset proportional coefficient. The danger unit value is affected by the weight of each box and the type of chemical corresponding to the danger type. When the gas generated by the corresponding chemical after physical or chemical reaction is more toxic, or the weight of each box is higher, the corresponding danger unit value is larger, and the ability to cause air environment risk is stronger.

[0037] It should be noted that gas types include not only gas types generated by chemical reactions with standard air (i.e. air with a content of nitrogen: 78.09%, oxygen: 20.95%, rare gases: 0.932%, carbon dioxide: 0.034%, water vapor and impurities: 0.02% and other dust), but also gas types generated by physical reactions of chemicals corresponding to the hazardous type (such as sublimation, gasification and other physical phenomena), as well as gas types generated by natural decomposition. The classification standard for acute toxicity categories adopts the Chinese national standard (GB 20592-2006) for classification. (GB 20592-2006) provides a detailed classification method, which divides the acute toxicity of gases into five categories based on toxicity data through the oral, dermal or inhalation routes, and each category has corresponding warning label elements and graphic symbols.

[0038] Obtain the number of containers corresponding to each hazard type in the dangerous goods container stack as the stack number n, and substitute the stack number n and the hazard unit value DE of the corresponding hazard type into the formula R danger =n*DE, and the hazard diffusion parameter R corresponding to each hazard type of dangerous goods container stack is obtained. danger .

[0039] It should be noted that the hazard diffusion parameters reflect the types of dangerous goods in the dangerous goods container pile, the degree of dangerous and harmful characteristics of the goods (toxicity intensity), and the storage order limit. The more toxic the gas produced by the dangerous goods or the larger the order of magnitude of the dangerous goods, the greater the potential risk, and therefore the larger the hazard diffusion parameters. The hazard diffusion parameters help staff to judge the magnitude and degree of the hazard of the dangerous goods container pile.

[0040] The dangerous goods container pile that generates toxic gas in an accident is recorded as an accident pile. The risk source positioning module obtains the specific location of the accident pile in the storage area, records the specific location as the danger occurrence point, and records the time when the danger occurrence point is determined as the occurrence time.

[0041] The environmental analysis module includes a ring analysis unit and a sector analysis unit. The ring analysis unit divides a certain range around the storage area into multiple ring areas, analyzes each ring area, and calculates its life impact index based on the resident data in the ring area. The sector analysis unit constructs a sector area according to the wind direction in the storage area and calculates the life impact index in the sector area.

[0042] The process of dividing the annular area by the annular analysis unit is as follows:

[0043] The hazard diffusion parameter of the dangerous goods container stack corresponding to the hazard occurrence point is obtained and recorded as R′ danger , obtain the toxicity value of the dangerous goods container pile corresponding to the danger point and record it as R′ danger Substitute into the formula , and the maximum impact range R1 is obtained, where h is the preset height value and r is the preset proportional coefficient, which represents the general height range reached by human life activities (that is, most of human life is below this height range, usually 100 meters). A circular area is drawn with the danger point in the storage area as the center and the maximum impact range R1 as the radius. Draw a circular boundary line and divide the circular area into 5 annular areas (the circular area to which the danger point belongs is regarded as a special ring with an inner diameter of 0).

[0044] It should be noted that, without considering the influence of wind, toxic gases will naturally diffuse slowly around, and the principle of gas diffusion is based on the thermal motion and mutual collision of gas molecules, as well as the existence of concentration gradient. Therefore, toxic gases will affect the environment in the surrounding circular area, and the concentration gradient of toxic gases will decrease with the increase of the distance from the danger point. After it is reduced to a certain range, it will no longer cause any impact. The maximum range refers to the farthest affected range of the toxic gas in the dangerous goods container pile. By demarcating different annular areas, the gradient range of the degree of impact can be better divided.

[0045] The process of constructing the fan-shaped area by the fan-shaped analysis unit is as follows:

[0046] The wind direction and wind speed of the storage area are obtained regularly, and the wind direction at the time of occurrence is recorded as the initial wind direction. The vertex of the danger occurrence point and the initial wind direction are used as the angle bisector to construct an area angle of P°, where P is a preset value (90 in a specific embodiment), and R′ danger , Substitute into the formula Calculation is performed to obtain the diffusion influence distance R2, and a circular outline is drawn with the vertex of the regional angle as the center and the diffusion distance R1 as the radius, and the circular outline intersects with the regional angle to form a fan-shaped area;

[0047] The angle difference between the current wind direction and the initial wind direction is calculated regularly. When the angle difference is less than or equal to the preset wind direction change threshold, the fan-shaped area corresponding to the initial wind direction is still used as the range area for calculating the life impact index. When the angle difference is greater than the preset wind direction change threshold, the current wind direction is used as the initial wind direction and the fan-shaped area is reconstructed, and the fan-shaped area is used as the range area for calculating the life impact index. In a specific embodiment, the wind direction change threshold is 30°.

[0048] It should be noted that under the action of wind, toxic gases will slowly diffuse in a fan shape, and the direction of diffusion is basically consistent with the wind direction (but it does not rule out that some gases will diffuse to the surroundings under the influence of other factors). Constructing fan-shaped areas based on wind direction can help staff better judge the air range affected by toxic gases based on wind factors, thereby judging the size and degree of the hazard.

[0049] The life impact index calculation process is as follows:

[0050] The area to be calculated is recorded as the target area. When the target area is a plurality of annular areas divided by the annular analysis unit, the plurality of annular areas are recorded as Q1, Q2, Q3, Q4, and Q5 in descending order according to the outer diameter. The number of personnel NUM1 in each annular area is counted, and the total number of commercial buildings, office buildings, industrial buildings, and educational buildings in each annular area is counted as the number of buildings NUM2. The number of personnel NUM1 and the number of buildings NUM2 are normalized and substituted into the formula δ i =e 6-i *(η1*NUM1+η2*NUM2) to obtain each annular area Q i Life impact index δ i , where i=1, 2, 3, 4, 5 correspond to regions Q1, Q2, Q3, Q4, Q5 respectively, η1, η2 are preset weight coefficients, and e is a natural constant;

[0051] It should be noted that the impact of the annular area is negatively correlated with the distance from the danger point, and the higher the distance, the higher the level of hazard.

[0052] When the target area is a fan-shaped area constructed by the fan-shaped analysis unit (there may be only one or more fan-shaped areas), the number of initial wind directions q within the time range from the occurrence time to the current time is obtained, and the duration T of each initial wind direction is calculated. p (that is, the time difference between the initial wind direction and the time when the next initial wind direction appears or the current time), P represents the sequence number of the initial wind direction, P = 1, 2, 3, ..., q, and each initial wind direction is analyzed separately: the duration T is plotted pA line graph of wind speed changes over time within the range is drawn, and the area of ​​the graph enclosed by the line graph and the horizontal axis of the coordinate is calculated and recorded as the wind force impact value F. The number of people NUM1′ and the number of buildings NUM2′ in the sector area are obtained, and the formula δ′=F*(η1*NUM1′+η2*NUM2′) is substituted for calculation to obtain the life impact index δ′ of the sector area. p .

[0053] It should be noted that the degree of impact in the fan-shaped area is positively correlated with the wind speed. The stronger the wind speed and the longer it lasts, the faster it spreads through the air, thus causing more serious and rapid air pollution, resulting in a higher level of hazard.

[0054] By calculating the life impact index in the target area, it can help workers determine the degree of impact and danger level of toxic gases generated by the leakage of dangerous goods in dangerous goods container piles on human life activities in the surrounding area under conditions without wind influence and with wind influence. By mathematically predicting the impact of dangerous goods on the air environment and combining it with the distribution of humans, the diffusion impact of toxic gases in two situations is further analyzed, which improves the efficiency of emergency response analysis after an accident occurs and facilitates workers to make reasonable judgments and decisions.

[0055] Comprehensive risk assessment module, to obtain the risk diffusion parameter R of the dangerous goods container stack corresponding to the risk occurrence point danger , each annular area Q i Life impact index δ i And the life impact index δ of different fan-shaped areas p ′, substitute it into the formula to calculate The risk assessment value DS is obtained, where μ1, μ2, and μ3 are all preset weight coefficients. The larger the risk assessment value, the greater the impact on the air environment and the greater the harm to human activities.

[0056] It should be noted that the risk assessment value is affected by many factors such as the degree of dangerous and harmful characteristics of dangerous goods in the dangerous goods container pile (toxicity intensity), the storage quantity, the personnel in the surrounding area of ​​the yard, etc. The risk assessment value obtained by calculation can combine various factors to comprehensively evaluate the degree of harm caused by the generation of flammable gases and toxic gases after an accident involving dangerous goods. The assessment result is more efficient, accurate, and comprehensive than the assessment method in the prior art, which helps to predict and judge the possible poisoning accidents after the dangerous goods container leakage accident occurs.

[0057] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements will not cause the essence of the corresponding technical solutions to deviate from the protection scope of the technical solutions of the embodiments of the present invention.

Claims

1. The environmental risk assessment method for dangerous goods container yards in ports is characterized by: include: The dangerous goods container storage statistics module, risk source location module, environmental analysis module, and comprehensive risk assessment module record the yard range as the storage area, and different types of dangerous goods containers are stacked separately to form dangerous goods container piles; Dangerous goods container storage statistics module, obtain the type of each dangerous goods container stack in the storage area and record it as the dangerous type, obtain the dangerous type of each dangerous goods container stack and the loading volume of each container and calculate its dangerous unit value DE, obtain the number of containers in each dangerous goods container stack and record it as the stack number n, substitute it into the formula R danger =n*DE to calculate the hazard diffusion parameter R of each dangerous goods container pile danger ; A pile of dangerous goods containers that generates toxic gases due to an accident shall be recorded as an accident pile; The risk source location module obtains the specific location of the accident pile in the storage area, records the specific location as the danger occurrence point, and records the time when the danger occurrence point is determined as the occurrence time; The environmental analysis module includes a ring analysis unit and a sector analysis unit. The ring analysis unit divides a certain range around the storage area into multiple ring areas Q i , and analyze each annular area, and calculate its life impact index δ based on the resident data in the annular area i The fan-shaped analysis unit monitors the wind direction angle in the storage area to construct multiple initial wind directions and corresponding fan-shaped areas, records the wind speed regularly, and calculates the life impact index δ corresponding to each fan-shaped area based on the human distribution characteristics and wind speed in each fan-shaped area. p ′ and record the duration T of each initial wind direction p ; The life impact index calculation process is as follows: The area to be calculated is recorded as the target area. When the target area is a plurality of annular areas divided by the annular analysis unit, they are recorded as Q1, Q2, Q3, Q4, and Q5 in descending order according to the outer diameter. The number of people NUM1 in each annular area is counted, and the total number of commercial buildings, office buildings, industrial buildings, and educational buildings in each annular area is counted as the number of buildings NUM2. The number of people NUM1 and the number of buildings NUM2 are normalized and substituted into the formula δ i =e 6-i *(η1*NUM1+η2*NUM2) to obtain each annular area Q i Life impact index δ i , where i = 1, 2, 3, 4, 5, η1, η2 are preset weight coefficients, and e is a natural constant; When the target area is a fan-shaped area constructed by the fan-shaped analysis unit, the number of initial wind directions q within the time range from the occurrence time to the current time is obtained, and the duration T of each initial wind direction is calculated. p , P represents the sequence number of the initial wind direction, P = 1, 2, 3, ..., q, and each initial wind direction is analyzed separately: Plot duration T p A line graph of wind speed changes over time within the range is drawn, and the area of ​​the graph enclosed by the line graph and the horizontal axis of the coordinate is calculated and recorded as the wind force impact value F. The number of people NUM1 and the number of buildings NUM2′ in the sector area are obtained, and the formula δ′=F*(η1*NUM1′+η2*NUM2′) is substituted for calculation to obtain the life impact index δ of the sector area. p ′; The comprehensive risk assessment module obtains the hazard diffusion parameter R of dangerous goods container stacks danger , each annular area Q i Life impact index δ i And the life impact index δ of different fan-shaped areas p ′, and the risk assessment value DS is obtained through comprehensive analysis.

2. The method for assessing the environmental risk of a port dangerous goods container yard according to claim 1 is characterized in that: The calculation process of the risk unit value is as follows: Get the gas type generated by the dangerous chemical reaction corresponding to the dangerous type, set a toxicity value D, assign the toxicity value according to the acute toxicity category to which the gas type belongs. When the category to which the gas type belongs is N, set its toxicity value to 6-N. When there are multiple gas types with acute toxicity, select the gas type with the largest corresponding toxicity value as the main toxicity source, and use the toxicity value of this gas type as the toxicity value of the dangerous type. Obtain the loading volume G per box and the toxicity value D and substitute them into the formula DE=k*G*10 D Calculation is performed to obtain the danger unit value DE, where k is a preset proportional coefficient.

3. The method for assessing the environmental risk of a port dangerous goods container yard according to claim 1 is characterized in that: The process of dividing the annular area by the annular analysis unit is as follows: The hazard diffusion parameter of the dangerous goods container stack corresponding to the hazard occurrence point is obtained and recorded as R′ danger , obtain the toxicity value of the dangerous goods container pile corresponding to the danger point and record it as R′ danger Substitute into the formula The maximum impact range R1 is obtained by calculation, where h is the preset height value and r is the preset proportional coefficient. A circular area is drawn with the danger point in the storage area as the center and the maximum impact range R1 as the radius. Draw a circular boundary line to divide the circular area into 5 annular areas.

4. The method for assessing the environmental risk of a port dangerous goods container yard according to claim 3 is characterized in that: The process of constructing the fan-shaped area by the fan-shaped analysis unit is as follows: The wind direction and wind speed of the storage area are obtained regularly, and the wind direction at the time of occurrence is recorded as the initial wind direction. The vertex of the danger occurrence point and the initial wind direction are used as the angle bisector to construct an area angle of P°, where P is the preset value. danger , Substitute into the formula Calculation is performed to obtain the diffusion influence distance R2, and a circular outline is drawn with the vertex of the regional angle as the center and the diffusion distance R1 as the radius, and the circular outline intersects with the regional angle to form a fan-shaped area; The angle difference between the current wind direction and the initial wind direction is calculated regularly. When the angle difference is less than or equal to the preset wind direction change threshold, the fan-shaped area corresponding to the initial wind direction is still used as the range area for calculating the life impact index. When the angle difference is greater than the preset wind direction change threshold, the current wind direction is used as the initial wind direction and the fan-shaped area is reconstructed, and the fan-shaped area is used as the range area for calculating the life impact index.

5. The method for assessing the environmental risk of a port dangerous goods container yard according to claim 1 is characterized in that: The risk comprehensive assessment module obtains the dangerous and harmful gas diffusion parameter R of the dangerous goods container stack corresponding to the dangerous occurrence point danger , each annular area Q i Life impact index δ i And the life impact index δ of different fan-shaped areas p ′, substitute it into the formula to calculate The risk assessment value DS is obtained, where μ1, μ2, and μ3 are all preset weight coefficients.

Images