A cfd-based qra consequence simulation method

Abstract

The application discloses a kind of QRA consequence simulation methods based on CFD, it is related to environmental simulation technical field, for solving the problem of reducing the accuracy of evaluation, the problem of increasing environmental pollution risk, by calling historical database to obtain waste liquid discharge historical information and setting two types of propagation benchmark, setting detection time and initial condition of discharge and calculating hindering characteristics, in the detection time, carry out simulation in advance to generate detection result, then based on initial condition, carry out formal simulation, the detection result and hindering characteristics are predicted to each time point propagation characteristics, compare the prediction result with two types of propagation benchmark and input fuzzy rule output, based on the output result, evaluate the propagation consequence of waste liquid under initial condition and input into interactive interface, significantly improve the scientificity and reliability of industrial waste liquid discharge risk assessment, provide more reference basis for environmental safety decision and emergency management, reduce environmental pollution risk, improve risk prediction result accuracy.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of environmental simulation, more particularly, the present application relates to a CFD-based QRA consequence simulation method. BACKGROUND

[0002] With the continuous acceleration of industrialization, chemical, metallurgical, energy and other industries will produce different types of industrial waste liquid in the production and disposal process. Industrial waste liquid usually contains a large amount of harmful or high-risk chemical components. Once abnormal discharge, leakage or diffusion occurs, it will cause air, water and soil pollution. In recent years, computational fluid dynamics (CFD) has been gradually applied to environmental simulation. It can realistically reflect the diffusion process of pollutants in air or water through numerical calculation.

[0003] The prior art has the following disadvantages:

[0004] At present, the existing simulation research based on CFD mainly focuses on numerical simulation of single scene, lacks effective combination with historical database, lacks quantitative modeling of hindering factors under the initial condition of discharge, and it is difficult to provide prediction basis with reference value for environmental safety decision-making, resulting in insufficient reliability of simulation results in risk assessment, reduced evaluation accuracy and increased environmental pollution risk. Therefore, a CFD-based QRA consequence simulation method is proposed.

[0005] The above information disclosed in the background section is only used to enhance the understanding of the background of the present disclosure, and therefore it can include information that does not constitute prior art known to those of ordinary skill in the art. SUMMARY

[0006] In order to overcome the above-mentioned defects of the prior art, the embodiments of the present application provide a CFD-based QRA consequence simulation method, which uses a historical database-driven CFD multi-stage linkage simulation and a propagation benchmark comparison evaluation mechanism to solve the problems raised in the above background.

[0007] To achieve the above-mentioned purpose, the present application provides the following technical scheme, a CFD-based QRA consequence simulation method, comprising the following steps:

[0008] Step S1: Before simulating the discharge of industrial waste liquid, call the historical database to obtain waste liquid discharge historical information, and set two types of propagation benchmarks based on the waste liquid discharge historical information;

[0009] Step S2: Set the detection time and the initial condition of the industrial waste liquid discharge, calculate the discharge hindering characteristics based on the initial condition of the discharge, and perform a preliminary simulation according to the initial condition of the industrial waste liquid discharge within the detection time, detect the waste liquid discharge information and generate a detection result;

[0010] Step S3: After generating the detection result, formal simulation discharge is performed based on the initial condition of the industrial waste liquid discharge. During the formal simulation discharge, the detection result and the discharge hindering feature are comprehensively used to predict the propagation feature at each time point during the simulation discharge process;

[0011] Step S4: The comparison rule is set, and the prediction result is compared with the two types of propagation benchmarks in the comparison rule. Based on the comparison result, the propagation consequence of the industrial waste liquid under the initial condition of the discharge is evaluated and is input into the interactive interface.

[0012] In a preferred embodiment, in step S1, before the simulation discharge of the industrial waste liquid, a query period is preset and is divided into multiple query time points;

[0013] The historical database is called to obtain the waste liquid discharge history information of the industrial waste liquid in the last time. The waste liquid discharge history information includes the diffusion area and the water surface layer wind force at each query time point;

[0014] The diffusion areas at adjacent query time points are subtracted to obtain diffusion area changes, and the diffusion area changes are combined into a diffusion area set;

[0015] The water surface layer wind forces at adjacent query time points are subtracted to obtain water surface layer wind force changes, and the water surface layer wind force changes are combined into a water surface layer wind force set;

[0016] Based on the waste liquid discharge history information, two types of propagation benchmarks are set, including a diffusion benchmark and a power benchmark.

[0017] In a preferred embodiment, in step S1, the median of the diffusion area set is taken as an area change median, each diffusion area change is subtracted from the area change median and the absolute value is taken to obtain an area deviation value, and the median of the area deviation value is taken as an area deviation median;

[0018] After the area change median and the area deviation median are respectively summed and subtracted, a diffusion first threshold and a diffusion second threshold are obtained;

[0019] The average value of the diffusion area changes between the diffusion first threshold and the diffusion second threshold is taken as the diffusion benchmark.

[0020] In a preferred embodiment, in step S1, the median of the water surface layer wind force set is taken as a wind force median;

[0021] If the water surface layer wind force change is greater than or equal to the wind force median, the water surface layer wind force change is divided into an upper interval group;

[0022] Otherwise, the water surface layer wind force change is divided into a lower interval group;

[0023] The median of the upper layer interval group and the median of the lower layer interval group are taken as the upper wind change median and the lower wind change median, respectively.

[0024] The average of the water body surface layer wind change amount between the upper wind change median and the lower wind change median is taken as the assistance reference.

[0025] In a preferred embodiment, in step S2, the detection time is set and the time points are divided at preset time intervals, and the discharge initial conditions include the diffusion initial velocity and the propagation average resistance;

[0026] The diffusion initial velocity and the propagation average resistance are obtained from an industrial discharge database;

[0027] The propagation average resistance is standardized to obtain a resistance coefficient;

[0028] The product of the resistance coefficient and the preset first adjustment parameter and the preset second adjustment parameter is taken as the discharge resistance characteristics of the diffusion area and the discharge resistance characteristics of the water body surface layer wind.

[0029] In a preferred embodiment, in step S2, based on the diffusion initial velocity and the propagation average resistance, a preliminary simulation is performed within the detection time to obtain waste liquid discharge information at different time points within the detection time;

[0030] The diffusion area change amount and the water body surface layer wind change amount at each time point are calculated according to the waste liquid discharge information;

[0031] The diffusion area change amount and the water body surface layer wind change amount at each time point are taken as the detection results at the corresponding time points.

[0032] In a preferred embodiment, in step S3, after generating the detection results, based on the discharge initial conditions of the industrial waste liquid, the diffusion and propagation process of the industrial waste liquid in the water body is formally simulated and discharged;

[0033] The interval between adjacent time points of the formal simulation discharge is set to be consistent with the interval between adjacent time points of the preliminary simulation;

[0034] The discharge resistance characteristics are called, and at each time point of the formal simulation discharge, based on the detection results and the discharge resistance characteristics, comprehensive operation is performed to obtain the propagation characteristics at the time point, including the diffusion area change amount prediction value and the wind change amount prediction value;

[0035] The propagation characteristics of all time points in the formal simulation discharge process are taken as the prediction results.

[0036] In a preferred embodiment, in step S4, the comparison rule is set, and the prediction result is compared with the two types of propagation reference in the comparison rule, and the comparison rule is as follows:

[0037] The diffusion area change amount prediction value at each time point is compared with the diffusion reference, and the wind force change amount prediction value at each time point is compared with the assistance reference;

[0038] If the diffusion area change amount prediction value is greater than the diffusion reference and the wind force change amount prediction value is greater than the assistance reference, the time point is defined as an abnormal time point;

[0039] Otherwise, the time point is defined as a normal time point.

[0040] In a preferred embodiment, in step S4, the number of abnormal time points is counted, and the propagation consequences of the industrial waste liquid under the initial conditions of discharge are evaluated:

[0041] When the number of abnormal time points is less than the preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions is within a controllable range, and no abnormal propagation risk is triggered;

[0042] When the number of abnormal time points is greater than or equal to the preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions has a significant risk, and the propagation consequences exceed the safety range;

[0043] The evaluation result is transmitted to the interactive interface.

[0044] The technical effects and advantages of the present application are:

[0045] The present application obtains the waste liquid discharge historical information by calling the historical database, sets two types of propagation reference based on the waste liquid discharge historical information, sets the detection time and the initial conditions of the industrial waste liquid discharge and calculates the discharge hindering characteristics, performs a prior simulation according to the initial conditions of the industrial waste liquid discharge within the detection time, detects the waste liquid discharge information and generates a detection result, performs a formal simulation discharge based on the initial conditions of the industrial waste liquid discharge, and predicts the propagation characteristics of each time point in the simulation discharge process by comprehensively considering the detection result and the discharge hindering characteristics during the formal simulation discharge. The prediction result is compared with the two types of propagation reference, the comparison result is input into the fuzzy rule for output, the propagation consequences of the industrial waste liquid under the initial conditions of discharge are evaluated based on the output result and input into the interactive interface, the scientificity and reliability of the risk assessment of the industrial waste liquid discharge are significantly improved, a more reference basis is provided for environmental safety decision and emergency management, the environmental pollution risk is reduced, and the accuracy of the risk prediction result is improved. BRIEF DESCRIPTION OF DRAWINGS

[0046] Fig. 1 It is a flowchart for the implementation of the QRA consequence simulation method based on CFD.

[0047] Fig. 2 A flowchart of steps of a CFD-based QRA consequence simulation method of the present application. DETAILED DESCRIPTION

[0048] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative effort fall within the protection scope of the present application.

[0049] The present application obtains waste liquid discharge historical information by calling a historical database, sets two types of propagation benchmarks based on the waste liquid discharge historical information, sets a detection time and an initial condition of industrial waste liquid discharge and calculates discharge hindering characteristics, performs preliminary simulation within the detection time according to the initial condition of industrial waste liquid discharge, detects waste liquid discharge information and generates a detection result, performs formal simulation discharge based on the initial condition of industrial waste liquid discharge, predicts propagation characteristics at each time point in the simulation discharge process by comprehensively considering the detection result and the discharge hindering characteristics, compares the prediction result with the two types of propagation benchmarks, outputs the comparison result after passing through a fuzzy rule, evaluates the propagation consequence of industrial waste liquid under the initial condition of discharge based on the output result, and inputs the evaluation result into an interactive interface, which significantly improves the scientificity and reliability of industrial waste liquid discharge risk assessment and provides a more referable basis for environmental safety decision and emergency management.

[0050] Embodiment 1

[0051] A CFD-based QRA consequence simulation method, as shown in Figs. 1-2 includes the following steps:

[0052] Step S1: Before simulating discharge of industrial waste liquid, historical database is called to obtain waste liquid discharge historical information, and two types of propagation benchmarks are set based on the waste liquid discharge historical information;

[0053] Step S2: A detection time and an initial condition of industrial waste liquid discharge are set, discharge hindering characteristics are calculated based on the initial condition of discharge, and preliminary simulation is performed within the detection time according to the initial condition of industrial waste liquid discharge, and waste liquid discharge information is detected and a detection result is generated;

[0054] Step S3: After the detection result is generated, formal simulation discharge is performed based on the initial condition of industrial waste liquid discharge, and propagation characteristics at each time point in the simulation discharge process are predicted by comprehensively considering the detection result and the discharge hindering characteristics in the formal simulation discharge process;

[0055] Step S4: Set the comparison rule, compare the prediction result with the two types of propagation benchmark into the comparison rule, evaluate the propagation consequence of the industrial waste liquid under the initial condition of discharge based on the comparison result and transmit into the interactive interface.

[0056] The implementation is as follows:

[0057] In step S1, before the simulation discharge of the industrial waste liquid, the query period is preset and divided into multiple query times, the historical database is called to obtain the waste liquid discharge history information of the industrial waste liquid at the last time, and the waste liquid discharge history information includes the diffusion area and the water surface layer wind force at each query time;

[0058] The diffusion area refers to the waste liquid coverage area formed by the industrial waste liquid in the horizontal direction after being discharged into the water body; the water surface layer wind force refers to the wind force acting on the water surface layer;

[0059] The diffusion area of each adjacent query time is subtracted to obtain the diffusion area change amount, and the diffusion area change amount is combined into a diffusion area set;

[0060] The water surface layer wind force of each adjacent query time is subtracted to obtain the water surface layer wind force change amount, and the water surface layer wind force change amount is combined into a water surface layer wind force set.

[0061] It should be noted that the preset query period refers to a fixed time window set for retrieving the waste liquid discharge history information; the historical database refers to a data set for storing and managing the historical information related to the industrial waste liquid discharge, and in this embodiment, it is used to obtain the waste liquid discharge history information.

[0062] Based on the waste liquid discharge history information, two types of propagation benchmarks are set, including diffusion benchmark and power benchmark;

[0063] The median of the diffusion area set is taken as the area change median, each diffusion area change amount is subtracted from the area change median and the absolute value is taken to obtain the area deviation value, and the median of the area deviation value is taken as the area deviation median;

[0064] The area change median and the area deviation median are respectively summed and subtracted to obtain the diffusion first threshold and the diffusion second threshold;

[0065] The average value of the diffusion area change amount between the diffusion first threshold and the diffusion second threshold is taken as the diffusion benchmark, and the diffusion benchmark reflects the diffusion ability of the industrial waste liquid in the historical data.

[0066] The median of the water surface layer wind force set is taken as the wind force median;

[0067] If the water surface layer wind force change amount is greater than or equal to the wind force median, the water surface layer wind force change amount is divided into the upper interval group;

[0068] Conversely, the amount of wind change of the water surface layer is divided into the lower interval group;

[0069] The median of the upper interval group and the lower interval group is taken as the upper wind change median and the lower wind change median, respectively;

[0070] The average of the amount of wind change of the water surface layer between the upper wind change median and the lower wind change median is taken as the assistance reference, which reflects the average assistance ability of the water surface layer wind to the diffusion rate of industrial waste liquid.

[0071] In this step, by calling the industrial waste liquid discharge information in the historical database, the diffusion reference and the assistance reference constructed by the median and the deviation are further improved to enhance the reference reference and adaptability of the prediction of the propagation trend in the subsequent simulation process, and the propagation risk quantitative analysis is realized for the actual scene.

[0072] In step S2, the detection time is set and the time points are divided according to the preset time interval, and the initial condition of discharge refers to the initial diffusion speed and environmental resistance condition of the simulation of discharge behavior, including diffusion initial speed and propagation average resistance;

[0073] The propagation average resistance refers to the environmental resistance of the industrial waste liquid during the discharge process in the water body, and the diffusion initial speed refers to the speed of the industrial waste liquid when it is released from the discharge port. The diffusion initial speed gradually approaches zero in the process of the simulation, until the diffusion stops naturally;

[0074] The monitoring data in the process of industrial discharge of waste liquid is retrieved through the industrial discharge database, and the diffusion initial speed and the propagation average resistance are obtained according to the monitoring data;

[0075] It should be noted that the detection time is used to limit the time range of the simulation process, which can be determined according to the discharge time in the historical data; the preset time interval reference is used to divide the continuous detection time into multiple time points, which can be set according to the simulation accuracy; the industrial discharge database refers to a management platform for recording real-time monitoring data in the process of industrial waste liquid discharge, which is used to record and call the monitoring information of discharge behavior formed in the process of industrial waste liquid discharge, including diffusion initial speed and propagation average resistance.

[0076] The preset first adjustment parameter and the preset second adjustment parameter are used to standardize the propagation average resistance to obtain the resistance coefficient, and the product of the resistance coefficient and the preset first adjustment parameter and the preset second adjustment parameter is taken as the discharge resistance characteristics of the diffusion area and the discharge resistance characteristics of the water surface layer wind;

[0077] The emission hindering feature of the diffusion area refers to that, under the natural diffusion state of the industrial waste liquid, the diffusion area of the industrial waste liquid is limited by the propagation average resistance, and reflects the limiting effect of the propagation average resistance on the spatial diffusion process of the industrial waste liquid in the water body.

[0078] The emission hindering feature of the surface wind of the water body refers to that, under the natural diffusion state of the industrial waste liquid, the surface wind of the water body is inhibited by the propagation average resistance to the diffusion rate of the waste liquid.

[0079] It should be explained that the preset first adjustment parameter is a quantitative coefficient for reflecting the corresponding relationship between the propagation average resistance and the diffusion area change amount, and has the same dimension as the diffusion area change amount, which can be determined according to the corresponding relationship between the propagation average resistance and the diffusion area in the historical emission data. For example, the average ratio of multiple sets of propagation average resistance values and their corresponding diffusion area change amounts in the historical emission data can be set, and the specific setting is performed by professional personnel; the preset second adjustment parameter is a quantitative coefficient for reflecting the corresponding relationship between the propagation average resistance and the surface wind change amount of the water body, and has the same dimension as the surface wind change amount of the water body, which can be determined according to the monitoring ratio data analysis of the propagation average resistance and the surface wind of the water body under different wind field conditions. For example, the average ratio of multiple sets of propagation average resistance values and their corresponding surface wind change amounts of the water body in the historical emission data can be set; the standardization processing mode includes but is not limited to standard linear transformation based on interval scaling, Z-Score standardization method based on statistics or normalization method based on nonlinear mapping function, and the application method of standardization processing is not described here.

[0080] According to the initial condition of the industrial waste liquid emission, the simulation in advance is performed within the detection time, and the simulation in advance refers to the simulation of the emission process of the industrial waste liquid, which is used to obtain the waste liquid emission information at different time points within the detection time;

[0081] According to the waste liquid emission information, the diffusion area change amount and the surface wind change amount of the water body at each time point are calculated; the diffusion area change amount and the surface wind change amount of the water body at each time point are taken as the detection results at the corresponding time points.

[0082] In this step, the detection time and the initial condition of the emission are set, so that the retardation control condition and the propagation convergence trend of the industrial waste liquid diffusion process in the simulation stage are provided, and the detection results generated by the simulation in advance provide data support for the formal simulation, and improve the stability and accuracy of the subsequent prediction.

[0083] In step S3, after the detection results are generated, the diffusion and propagation process of the industrial waste liquid in the water body is formally simulated based on the initial condition of the industrial waste liquid emission.

[0084] It should be noted that the formal simulation discharge refers to the core simulation calculation process of the waste liquid diffusion process based on the initial conditions of the industrial waste liquid discharge after the preliminary simulation is completed and the test results are obtained.

[0085] In the formal simulation discharge process, the test results and the discharge hindering characteristics are comprehensively considered to predict the propagation characteristics at each time point. The time sequence of the formal simulation discharge is set to be consistent with the preliminary simulation, and the interval between adjacent time points is consistent with the interval between the preliminary simulation test time points. That is, if the preliminary simulation is performed once every second to collect data, the formal simulation also takes every second as the time step to ensure the one-to-one correspondence between the prediction results and the test results.

[0086] For each time point, first read the test results generated in the preliminary simulation:

[0087] ;

[0088] Among them, is the test result at the i-th time point, is the diffusion area change amount at the i-th time point, is the water surface wind change amount at the i-th time point.

[0089] Secondly, the discharge hindering characteristics are called. At each time point of the formal simulation discharge, the test results and the discharge hindering characteristics are comprehensively operated to obtain the propagation characteristics at the time point, including the diffusion area change amount prediction value and the wind change amount prediction value:

[0090] ;

[0091] Among them, is the diffusion area change amount prediction value at the i-th time point, is the discharge hindering characteristic of the diffusion area, is the wind change amount prediction value at the i-th time point, is the discharge hindering characteristic of the water surface wind.

[0092] For the diffusion area change amount prediction value, the covered area of the industrial waste liquid in the water body will only remain unchanged or expand when the industrial waste liquid diffuses in the water body. If , then directly assign to 0;

[0093] For the wind change amount prediction value, because the water surface wind change amount has natural fluctuations such as enhancement and weakening, if , the negative value result is retained to reflect the actual environmental state that the environmental resistance weakens the wind power boosting effect and the wind power shows a weakening trend;

[0094] The surface wind variation of the water body is an environmental condition and only affects the diffusion rate of the industrial waste liquid. When the surface wind variation of the water body at two continuous time points is zero, it indicates that the environment of the water body does not provide an external driving force for the diffusion process of the industrial waste liquid, and at this time, the industrial waste liquid is in a natural diffusion state.

[0095] In the above operation process, if the surface wind variation of the water body in the detection result satisfies , then the propagation feature is only determined by the detection result of the diffusion area minus the hindering feature, that is, it satisfies:

[0096] ;

[0097] When the diffusion area variation at two continuous time points is zero, it indicates that the diffusion range of the industrial waste liquid will no longer expand thereafter, that is, the diffusion process has reached a stationary state.

[0098] Further, if the diffusion area variation in the detection result satisfies , it is determined that the industrial waste liquid no longer diffuses at this time point, and further, the subsequent prediction operation is terminated.

[0099] The propagation features of all time points in the formal simulation discharge process are taken as the prediction results to participate in subsequent operations.

[0100] In step S4, a comparison rule is set, and the prediction results are compared with two types of propagation reference in the comparison rule. The comparison rule is as follows:

[0101] The propagation features of each time point in the prediction result are compared with the corresponding propagation reference, that is, the diffusion area variation prediction value of each time point is compared with the diffusion reference, and the wind variation prediction value of each time point is compared with the power reference.

[0102] If the diffusion area variation prediction value is greater than the diffusion reference and the wind variation prediction value is greater than the power reference, the time point is defined as an abnormal time point.

[0103] Otherwise, the time point is defined as a normal time point.

[0104] The number of abnormal time points is counted to evaluate the propagation consequences of the industrial waste liquid under the initial conditions of discharge:

[0105] When the number of abnormal time points is less than a preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions is within a controllable range and does not trigger an abnormal propagation risk.

[0106] When the number of abnormal time points is greater than or equal to the preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions has a significant risk, and the propagation consequences exceed the safety range.

[0107] It should be noted that the quantity threshold is dynamically set through historical operation data and optimization training, which is not described here.

[0108] The above evaluation results are transmitted to the interactive interface, so that the user can intuitively understand the propagation dynamics of the industrial waste liquid under the initial discharge condition, thereby making quick decisions, emergency disposal or adjusting the discharge scheme in actual application.

[0109] The interactive interface is a human-computer interaction platform for presenting and interacting with the prediction results of the propagation consequences of the industrial waste liquid, outputting the numerical results and state discrimination results generated in the simulation operation in a visual manner, and providing information query and feedback functions for the user.

[0110] It should be noted that the initial discharge conditions of the industrial waste liquid can be dynamically adjusted according to actual operation data to adapt to the simulation requirements in different scenarios.

[0111] By setting the comparison rule and the quantity threshold, the diffusion process of the industrial waste liquid under the initial discharge condition is quantitatively evaluated to realize accurate determination of abnormal diffusion risk; combined with the visual interactive interface, the user can intuitively understand the waste liquid propagation dynamics and make quick scientific decisions and emergency disposal, thereby significantly improving the reliability and response efficiency of environmental safety management.

[0112] Finally, it should be noted that in this document, the relationship terms such as first and second are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations.

[0113] Moreover, the terms "include", "contain" or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method, article or equipment including a series of elements not only includes those elements, but also includes other elements not explicitly listed or inherent to such process, method, article or equipment. Without more limitations, the element defined by the statement "includes a" does not exclude the presence of another identical element in the process, method, article or equipment including the element.

[0114] In this document, the singular forms "a", "an" and "the" can also include the plural forms unless the context clearly indicates otherwise. It should also be understood that the term "includes / contains" or "has" and the like specifies the presence of the stated features, integers, steps, operations, components, parts or combinations thereof, but does not exclude the presence or addition of one or more other features, integers, steps, operations, components, parts or combinations thereof. The term "and / or" used in this specification includes any and all combinations of the related listed items.

[0115] The various embodiments described in this specification are presented for the purpose of illustration and description. Each of the embodiments described in this specification can be implemented alone or in combination with others, and the same or similar parts are cross-referenced where applicable.

[0116] The above description of disclosed embodiments is intended to be illustrative and not restrictive. Many modifications of the embodiments described herein will obvious to those with skill and / or knowledge in the field, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the application. Therefore, the application is not to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A CFD-based QRA consequence simulation method, characterized in that: The method comprises the following steps: Step S1: Before simulating the industrial waste liquid, the historical database is called to obtain the waste liquid discharge history information, and two types of propagation benchmarks are set based on the waste liquid discharge history information; Step S2: The detection time and the initial condition of the industrial waste liquid discharge are set, the discharge hindering characteristics are calculated based on the initial condition, the preliminary simulation is performed within the detection time according to the initial condition of the industrial waste liquid discharge, the waste liquid discharge information is detected, and the detection result is generated; In step S2, the detection time is set and the time points are divided according to the preset time interval, the initial condition of the discharge includes the initial diffusion speed and the average propagation resistance; The initial diffusion speed and the average propagation resistance are obtained through the industrial discharge database; The average propagation resistance is standardized to obtain the resistance coefficient; The product of the resistance coefficient and the preset first adjustment parameter and the preset second adjustment parameter is taken as the discharge hindering characteristics of the diffusion area and the water surface wind force; Step S3: After the detection result is generated, the formal simulation discharge is performed based on the initial condition of the industrial waste liquid discharge, and the propagation characteristics of each time point in the simulation discharge process are predicted by comprehensively considering the detection result and the discharge hindering characteristics; Step S4: The comparison rule is set, the prediction result is compared with the two types of propagation benchmarks in the comparison rule, and the propagation consequences of the industrial waste liquid under the initial condition of the discharge are evaluated based on the comparison result and are input into the interactive interface; In step S4, the comparison rule is set, the prediction result is compared with the two types of propagation benchmarks in the comparison rule, and the comparison rule is as follows: The diffusion area change prediction value of each time point is compared with the diffusion benchmark, and the wind force change prediction value of each time point is compared with the assistance benchmark; If the diffusion area change prediction value is greater than the diffusion benchmark and the wind force change prediction value is greater than the assistance benchmark, the time point is defined as an abnormal time point; Otherwise, the time point is defined as a normal time point.

2. The CFD-based QRA consequence simulation method according to claim 1, wherein: In step S1, before simulating the industrial waste liquid, a query period is preset and divided into a plurality of query time points; The historical database is called to obtain the waste liquid discharge history information of the industrial waste liquid at the last time, and the waste liquid discharge history information includes the diffusion area and the water surface wind force at each query time point; The diffusion areas of adjacent query time points are subtracted to obtain diffusion area changes, and the diffusion area changes are combined into a diffusion area set; The water surface wind forces of adjacent query time points are subtracted to obtain water surface wind force changes, and the water surface wind force changes are combined into a water surface wind force set; Two types of propagation benchmarks are set based on the waste liquid discharge history information, and the propagation benchmarks include a diffusion benchmark and an assistance benchmark.

3. The CFD-based QRA consequence simulation method according to claim 2, wherein: In step S1, the median of the diffusion area set is taken as the area change median, the diffusion area changes are subtracted from the area change median to obtain area deviation values, and the median of the area deviation values is taken as the area deviation median. Summing and subtracting the median of the area change and the median of the area deviation respectively to obtain a diffusion first threshold and a diffusion second threshold; Taking the average of the diffusion area change between the diffusion first threshold and the diffusion second threshold as the diffusion reference.

4. The CFD-based QRA consequence simulation method according to claim 2, wherein: In step S1, the median of the wind force set of the water surface layer is taken as the wind force median; If the water surface layer wind force change is greater than or equal to the wind force median, the water surface layer wind force change is divided into the upper interval group; Otherwise, the water surface layer wind force change is divided into the lower interval group; The medians of the upper interval group and the lower interval group are taken as the upper wind force change median and the lower wind force change median respectively; The average of the water surface layer wind force change between the upper wind force change median and the lower wind force change median is taken as the assistance reference.

5. The CFD-based QRA consequence simulation method according to claim 1, wherein: In step S2, based on the diffusion initial velocity and the propagation average resistance, a preliminary simulation is performed within the detection time to obtain the waste liquid discharge information at different time points within the detection time; The diffusion area change and the water surface layer wind force change at each time point are calculated according to the waste liquid discharge information; The diffusion area change and the water surface layer wind force change at each time point are taken as the detection results at the corresponding time points.

6. The CFD-based QRA consequence simulation method according to claim 1, wherein: In step S3, after generating the detection results, based on the initial conditions of the industrial waste liquid discharge, a formal simulation discharge of the diffusion and propagation process of the industrial waste liquid in the water body is performed; The interval of adjacent time points in the formal simulation discharge is set to be consistent with the interval of adjacent time points in the preliminary simulation; The discharge hindering features are called, and at each time point of the formal simulation discharge, based on the detection results and the discharge hindering features, a comprehensive operation is performed to obtain the propagation features at the time point, including the diffusion area change prediction value and the wind force change prediction value; The propagation features at all time points in the formal simulation discharge process are taken as the prediction results.

7. The CFD-based QRA consequence simulation method according to claim 1, wherein: In step S4, the number of abnormal time points is counted to evaluate the propagation consequences of the industrial waste liquid under the initial conditions of the discharge: When the number of abnormal time points is less than a preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions is within a controllable range and does not trigger an abnormal propagation risk; When the number of abnormal time points is greater than or equal to the preset number threshold, it is determined that the diffusion process of the waste liquid under the set conditions has a significant risk and the propagation consequences exceed the safety range; The evaluation results are transmitted to the interactive interface.

Images