Pollution influence range defining method and device, storage medium and electronic equipment

Abstract

The invention provides a pollution influence range defining method and device, a storage medium and electronic equipment. The method comprises the steps that the particle mass of each pollutant particle in M pollutant particles is determined; determining the diffusion simulation result of each pollutant particle in the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, and recording the position information of each pollutant particle under the current drift diffusion simulation time step; on the basis of the diffusion simulation result of each pollutant particle, the grid area mass concentration of each concentration judgment grid area in the target area is updated; and when the grid area mass concentration of each concentration judgment grid area does not reach a diffusion simulation stop condition, continuing drifting diffusion simulation to obtain at least one piece of position information of each initial pollutant particle in the initial pollutant particle set so as to determine a pollution influence range. The embodiment of the invention can improve the accuracy of the pollution influence range.

Description

Technical Field

[0001] The present invention relates to the field of computer technology, and in particular to a method, device, storage medium and electronic equipment for defining a pollution impact range. Background Art

[0002] At present, the definition of pollution impact scope (also known as pollution accident impact scope) has been widely used, especially in water ecological loss assessment (i.e. water pollution accident loss assessment) for water pollution accidents, so as to use the pollution impact scope for water pollution accident loss assessment, etc. However, related technologies usually use manual experience to define the pollution impact scope, resulting in low accuracy of the pollution impact scope. Based on this, there is currently no good solution for how to define the pollution impact scope to improve the accuracy of the pollution impact scope. Summary of the invention

[0003] In view of this, the embodiments of the present invention provide a method, device, storage medium and electronic device for defining the scope of pollution impact, so as to solve the problem of low accuracy of the pollution impact scope defined by the relevant technology; that is, the embodiments of the present invention can realize the definition of the pollution impact scope by simulating the drift and diffusion of pollutant particles, which can effectively improve the accuracy of the pollution impact scope, that is, the embodiments of the present invention can realize high-precision definition of the impact scope of sudden water pollution accidents. Then, when the pollution impact scope defined by the embodiments of the present invention is used for water pollution accident loss assessment, it can provide technical support for the accurate quantitative assessment of water ecological losses, that is, the pollution impact scope with higher accuracy can be used for water pollution accident loss assessment, thereby improving the accuracy of water pollution accident loss assessment.

[0004] According to one aspect of an embodiment of the present invention, a method for defining a pollution impact range is provided, the method comprising: Acquire target flow field data of a target area within a target time range, and acquire target meteorological field data of the target area within the target time range; Determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutant particles, and record the position information of each pollutant particle among the P pollutant particles at the current drift diffusion simulation time step; wherein one diffusion simulation result includes the particle mass and position information of one pollutant particle, and M and P are both positive integers; Based on the diffusion simulation result of each of the P pollutant particles, respectively updating the grid area mass concentration of each concentration determination grid area in the target area; When the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; Based on the at least one position information of each initial pollutant particle, a pollution impact range is determined.

[0005] According to another aspect of an embodiment of the present invention, a device for defining a pollution impact range is provided, the device comprising: An acquisition unit, used for acquiring target flow field data of a target area within a target time range, and acquiring target meteorological field data of the target area within the target time range; A processing unit, used to determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutant particles, and record the position information of each pollutant particle among the P pollutant particles at the current drift diffusion simulation time step; wherein one diffusion simulation result includes the particle mass and position information of one pollutant particle, and M and P are both positive integers; The processing unit is further used to update the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each pollutant particle in the P pollutant particles; The processing unit is further configured to continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; The processing unit is further used to determine the pollution impact range based on at least one position information of each initial pollutant particle.

[0006] According to another aspect of an embodiment of the present invention, an electronic device is provided, comprising a processor and a memory storing a program, wherein the program comprises instructions, which, when executed by the processor, enable the processor to perform the above-mentioned method.

[0007] According to another aspect of an embodiment of the present invention, a non-transitory computer-readable storage medium storing computer instructions is provided, wherein the computer instructions are used to enable a computer to execute the above-mentioned method.

[0008] The embodiment of the present invention can determine the particle mass of each pollutant particle in the M pollutant particles after obtaining the target flow field data of the target area within the target time range and the target meteorological field data of the target area within the target time range; and determine the diffusion simulation result of each pollutant particle in the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, and record the position information of each pollutant particle in the P pollutant particles at the current drift diffusion simulation time step; wherein, a diffusion simulation result includes the particle mass and position information of a pollutant particle, and M and P are both positive integers. Further, based on the diffusion simulation result of each pollutant particle in the P pollutant particles, the grid area mass concentration of each concentration determination grid area in the target area can be updated respectively; when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, the drift diffusion simulation is continued based on the diffusion simulation result of each pollutant particle in the P pollutant particles, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set. Based on this, the pollution impact range can be determined based on at least one position information of each initial pollutant particle. It can be seen that the embodiment of the present invention can achieve the definition of the pollution impact range by simulating the drift and diffusion of pollutant particles, which can effectively improve the accuracy of the pollution impact range. In other words, the embodiment of the present invention can achieve high-precision definition of the impact range of sudden water pollution accidents; based on this, when the pollution impact range defined by the embodiment of the present invention is used for water pollution accident loss assessment, it can provide technical support for the accurate quantitative assessment of water ecological losses, that is, the pollution impact range with higher accuracy can be used for water pollution accident loss assessment, thereby improving the accuracy of water pollution accident loss assessment. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] Further details, features and advantages of the invention are disclosed in the following description of exemplary embodiments in conjunction with the accompanying drawings, in which: Figure 1 A schematic flow chart of a method for defining a pollution impact range according to an exemplary embodiment of the present invention is shown; Figure 2 A schematic flow chart of another method for defining a pollution impact range according to an exemplary embodiment of the present invention is shown; Figure 3A schematic block diagram of a device for defining a pollution impact range according to an exemplary embodiment of the present invention is shown; Figure 4 A block diagram of an exemplary electronic device that can be used to implement an embodiment of the present invention is shown. DETAILED DESCRIPTION

[0010] Embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although certain embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms and should not be construed as being limited to the embodiments described herein, which are instead provided for a more thorough and complete understanding of the present invention. It should be understood that the drawings and embodiments of the present invention are only for exemplary purposes and are not intended to limit the scope of protection of the present invention.

[0011] It should be understood that the various steps described in the method embodiments of the present invention may be performed in different orders and / or in parallel. In addition, the method embodiments may include additional steps and / or omit the steps shown. The scope of the present invention is not limited in this respect.

[0012] The term "including" and its variations used in this document are open inclusions, that is, "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one other embodiment"; the term "some embodiments" means "at least some embodiments". Relevant definitions of other terms will be given in the following description. It should be noted that the concepts of "first", "second", etc. mentioned in the present invention are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0013] It should be noted that the modifications of "one" and "plurality" mentioned in the present invention are illustrative rather than restrictive, and those skilled in the art should understand that, unless otherwise clearly indicated in the context, it should be understood as "one or more".

[0014] The names of the messages or information exchanged between multiple devices in the embodiments of the present invention are only used for illustrative purposes, and are not used to limit the scope of these messages or information.

[0015] It should be noted that the execution subject of the pollution impact scope definition method provided in the embodiment of the present invention may be one or more electronic devices, and the present invention does not limit this; wherein, the electronic device may be a terminal (i.e., a client) or a server, then when the execution subject includes multiple electronic devices, and the multiple electronic devices include at least one terminal and at least one server, the pollution impact scope definition method provided in the embodiment of the present invention may be jointly executed by the terminal and the server. Accordingly, the terminals mentioned here may include but are not limited to: smart phones, tablet computers, laptops, desktop computers, intelligent voice interaction devices, smart home appliances, and the like. The server mentioned here may be an independent physical server, or a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides cloud services, cloud databases, cloud computing (cloud computing), cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN (Content Delivery Network), and basic cloud computing services such as big data and artificial intelligence platforms, and the like.

[0016] Based on the above description, an embodiment of the present invention proposes a method for defining a pollution impact range, which can be executed by the electronic device (terminal or server) mentioned above; or, the method for defining a pollution impact range can be executed by the terminal and the server together. For the sake of convenience, the following description will be given by taking the electronic device executing the method for defining a pollution impact range as an example; Figure 1 As shown, the pollution impact range definition method may include the following steps S101-S105: S101, obtaining target flow field data of a target area within a target time range, and obtaining target meteorological field data of a target area within a target time range.

[0017] Among them, the target area can be any jurisdiction area, that is, it can be any water area, and the embodiment of the present invention does not limit this. Optionally, the target time range can be any time range, and the embodiment of the present invention does not limit this; it should be understood that the target time range is greater than or equal to the time range when the grid area mass concentration of each concentration determination grid area below reaches the diffusion simulation stop condition, such as the target time range may include the time range from the monitoring of the water pollution accident (i.e., water pollution incident) to the time when the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition. Optionally, the target time range may include one or more time sub-ranges, that is, the target time range may be composed of one or more time sub-ranges, then the flow field data of the target area in each time sub-range of the one or more time sub-ranges may be obtained respectively to obtain the target flow field data of the target area within the target time range, and the meteorological field data of the target area in each time sub-range may be obtained respectively to obtain the target meteorological field data of the target area within the target time range. For example, the meteorological field data and flow field data may be obtained on a daily basis, at which time the meteorological field data and flow field data may be continuously obtained during the drift diffusion simulation process to obtain the target meteorological field data and target flow field data; alternatively, the target time range may be set to a longer time range, so as to directly simulate the target meteorological field data and target flow field data of the target area within the target time range, and so on.

[0018] In an embodiment of the present invention, the target flow field data may include sub-flow field data of each flow field grid area in the target area within the target time range. One sub-flow field data may include flow field information of the corresponding flow field grid area at each flow field calculation time (i.e., flow field time step) within the target time range. One flow field information may include at least one of the following: the water level of the corresponding flow field grid area at the corresponding flow field calculation time, the flow velocity component (also called flow velocity) in each direction in at least one direction (such as the east direction and the north direction, etc.), salinity, water temperature and suspended sediment concentration, etc.; the embodiment of the present invention is not limited to this.

[0019] Correspondingly, the target meteorological field data may include sub-meteorological field data of each meteorological field grid area in the target area within the target time range, and one sub-meteorological field data may include meteorological field information of each meteorological field calculation time (i.e., meteorological field time step) of the corresponding meteorological field grid area within the target time range, and one meteorological field information may include at least one of the following: at least one wind speed, atmospheric pressure, and temperature of the corresponding meteorological field grid area at the corresponding meteorological field calculation time, etc.; the embodiment of the present invention does not limit this. Among them, one wind speed corresponds to one direction, that is, one meteorological field information may include the wind speed (i.e., wind speed component) under each wind direction in at least one wind direction.

[0020] Optionally, each flow field grid area and each meteorological field grid area in the target area may be the same or different, and the embodiment of the present invention does not limit this; that is, the size of a flow grid area and the size of a meteorological field grid area may be the same or different. Correspondingly, each flow field calculation time and each meteorological field calculation time in the target time range may be the same or different, and the embodiment of the present invention does not limit this; that is, the time interval between any two flow field calculation times may be the same as or different from the time interval between any two meteorological field calculation times, that is, the flow field time step and the meteorological field time step may be the same or different, and the embodiment of the present invention does not limit this. In other embodiments, a flow field information may also be the flow field information of a flow field grid area within a flow field sub-time range, and a meteorological field information may also be the meteorological field information of a meteorological field grid area within a meteorological field sub-time range, and the embodiment of the present invention does not limit this; for the convenience of explanation, the following description is taken as an example that a flow field information is the flow field information of a flow field grid area at a flow field calculation time, and a meteorological field information is the meteorological field information of a meteorological field grid area at a meteorological field calculation time.

[0021] It should be noted that each flow field grid area in the target area can be determined based on the shoreline direction and / or water terrain distribution of the target area, and the embodiment of the present invention does not limit the specific division method of the flow field grid area. Correspondingly, the embodiment of the present invention does not limit the specific division method of the meteorological field grid area.

[0022] Optionally, the target meteorological field data may be obtained in the following ways, including but not limited to: The first acquisition method: when acquiring the target meteorological field data of the target area within the target time range, the electronic device can download the meteorological field data file when detecting the start-up requirement of the target task. The meteorological field data file may include the forecast meteorological field data source of the target area within the target time range; and convert the format of the meteorological field data file to obtain the target meteorological field data. The target meteorological field data supports the network common data format (NetCDF format). It can be seen that the embodiment of the present invention can convert the acquired meteorological field data file into the target meteorological field data readable by the numerical model (such as the water quality simulation model) through the developed file conversion program, and the target meteorological field data can be in the network common data format. Optionally, the electronic device can obtain the meteorological field data download link and use the meteorological field data download link to download the meteorological field data file. Optionally, the meteorological field data file may include one or more meteorological field data sub-files, that is, one or more meteorological field data sub-files may be downloaded separately, and a meteorological field data sub-file may include the forecast meteorological field data source of the target area within a time sub-range, so that each meteorological field data sub-file in the one or more meteorological field data sub-files is format converted separately to obtain the meteorological field data of the target area within each time sub-range, so as to obtain the target meteorological field data.

[0023] Optionally, the target task may refer to a task for defining the scope of pollution impact for the target area. Optionally, the electronic device may include a system task planning tool, and trigger the start-up requirements of the target task through the system task planning tool, such as automatically triggering the start-up requirements of the target task at a specified time every day, thereby triggering the meteorological field data acquisition script to obtain the target meteorological field data, etc., or automatically triggering the start-up requirements of the target task at a preset interval, or triggering the start-up requirements of the target task when a sudden water pollution accident is detected (such as triggering the start-up requirements of the target task when the concentration monitoring data of the target pollutant is detected to be greater than the preset concentration monitoring threshold), etc.; illustratively, the task can be planned based on the system task planning tool to perform daily timed triggering, daily timed download of the next 10 days of meteorological field data, and put it in a specified folder. Optionally, the specified folder, specified time, preset interval duration, and preset concentration monitoring threshold can all be set according to experience, or can be set according to actual needs, and the embodiment of the present invention does not limit this. Optionally, the target pollutant can be any pollutant, and the embodiment of the present invention does not limit this.

[0024] Optionally, the target meteorological field data can be obtained based on a forecast meteorological field data source of a global atmospheric model (Global Forecast System, GFS, also known as the global forecast system); wherein, when the meteorological field data only includes wind field data, it can also be obtained based on a forecast wind field data source.

[0025] The second acquisition method: The electronic device can obtain the meteorological data of the target area, and determine the target meteorological field data based on the meteorological data through the Weather Research and Forecasting Model (WRF), etc.

[0026] Optionally, the target flow field data may be obtained in the following ways, including but not limited to: The first acquisition method: when acquiring the target flow field data of the target area within the target time range, the electronic device can acquire the initial boundary data source, which may include at least one historical data to be processed; and preprocess the data of each historical data to be processed in the at least one historical data to be processed respectively, to obtain the preprocessed boundary data corresponding to each historical data to be processed, and one preprocessed boundary data includes a boundary data of the target area within the target time range; further, based on the preprocessed boundary data corresponding to each historical data to be processed, the target boundary data source can be determined, and the target hydrodynamic model can be called to determine the target flow field data based on the target boundary data source; wherein the target flow field data may include the sub-flow field data of each flow field grid area in the target area within the target time range. Specifically, when determining the target boundary data source based on the preprocessed boundary data corresponding to each historical data to be processed, the preprocessed boundary data corresponding to each historical data to be processed can be added to the target boundary data source, so that the target boundary data source includes the preprocessed boundary data corresponding to each historical data to be processed.

[0027] Among them, at least one historical data to be processed may include historical data to be processed at the upstream boundary and / or historical data to be processed at the downstream boundary; optionally, one historical data to be processed may include hydrological data of each flow field grid area in any boundary (such as the upstream boundary or the downstream boundary) of the target area at each monitoring time among multiple monitoring times (that is, one historical data to be processed may correspond to a flow field grid area in any boundary), and one hydrological data may include at least one hydrological sub-information, and at least one hydrological sub-information may include but is not limited to: water level and flow velocity in each direction in at least one direction. Optionally, when performing data preprocessing on each of at least one of the historical data to be processed to obtain preprocessing boundary data corresponding to each of the historical data to be processed, for any of the at least one of the historical data to be processed, any of the flow field grid areas corresponding to any of the historical data to be processed, and any of the hydrological sub-information in at least one of the hydrological sub-information, the electronic device can determine any hydrological sub-information of any flow field grid area at each monitoring time from any of the historical data to be processed, and perform regression analysis on any of the hydrological sub-information of any of the flow field grid areas at each monitoring time, so as to obtain any of the hydrological sub-information of any of the flow field grid areas at each flow field calculation time, so as to realize data preprocessing of any of the hydrological sub-information of any of the flow field grid areas at each monitoring time, thereby realizing data preprocessing of any of the historical data to be processed; in this case, the preprocessing boundary data corresponding to any of the historical data to be processed may include any of the hydrological sub-information of any of the flow field grid areas at each flow field calculation time. Based on this, the preprocessing boundary data corresponding to any historical data to be processed may include any hydrological sub-information of any flow field grid area at each flow field calculation time in any time sub-range, that is, the preprocessing boundary data corresponding to any historical data to be processed may include the preprocessing boundary sub-data within any time sub-range (the preprocessing boundary sub-data within any time range may include any hydrological sub-information of any flow field grid area at each flow field calculation time in any time sub-range), so that based on the preprocessing boundary sub-data within any time sub-range in each historical data to be processed, the boundary data source corresponding to any time sub-range can be determined (the target boundary data source may include the boundary data source corresponding to each time sub-range), so as to call the target hydrodynamic model, and determine the flow field data of the target area within any time sub-range based on the boundary data source corresponding to any time sub-range, so as to obtain the target flow field data.

[0028] Optionally, when the flow field data of the upstream boundary and the downstream boundary of the target area are similar, the number of the historical data to be processed in at least one historical data to be processed may be one (such as only including the historical data to be processed of the upstream boundary or the historical data to be processed of the downstream boundary); illustratively, at least one historical data to be processed may include the historical data to be processed of the upstream boundary, in which case the pre-processed boundary data corresponding to the historical data to be processed of the downstream boundary may be the same as the pre-processed boundary data corresponding to the historical data to be processed of the upstream boundary. Optionally, when the downstream boundary is a tidal water area, at least one historical data to be processed may include the historical data to be processed of the upstream boundary; in this case, the initial boundary data source may also include the hydrological data in each flow field grid area in the downstream boundary of the target area, and the hydrological data in each flow field grid area in the downstream boundary of the target area may be added to the target boundary data source, so that the target boundary data source includes the hydrological data in each flow field grid area in the downstream boundary of the target area. It can be seen that the embodiment of the present invention can design the upstream and downstream boundaries of the hydrodynamic model based on the water characteristics of the target area, that is, design the flow field data of the upstream and downstream boundaries.

[0029] Optionally, the electronic device may include a data preprocessing module. In this case, the electronic device may perform data preprocessing through the data preprocessing module.

[0030] Optionally, when calling the target hydrodynamic model and determining the target flow field data based on the target boundary data source, the electronic device may call the target hydrodynamic model and determine the target flow field data based only on the target boundary data source, or may determine the target flow field data based on the target boundary data source and the target meteorological field data, and so on; the embodiment of the present invention is not limited to this. Optionally, when determining the target flow field data based on the target boundary data source, the intermediate flow field data may be determined based on the target boundary data source, and the intermediate flow field data may be converted into the target flow field data in the hydrodynamic flow field data format (such as the network general data format), such as the target flow field data may be supported as the network general data format, so as to facilitate the rapid calling of the water quality simulation model (also known as the pollutant tracing model).

[0031] In an embodiment of the present invention, the electronic device may also obtain an initial hydrodynamic model corresponding to the target area, and call the initial hydrodynamic model to predict the calibrated flow field data of the target area within the historical time range, that is, call the initial hydrodynamic model, and predict the calibrated flow field data of the target area within the historical time range based on the boundary data and / or meteorological field data of the target area within the historical time range; optionally, the historical time range may be any time range, which is not limited in the embodiment of the present invention. Then, the initial hydrodynamic model may be calibrated based on the calibrated flow field data and the benchmark data of the target area within the historical time range to obtain a calibrated hydrodynamic model; wherein the benchmark data may include the monitoring hydrological data of each monitoring point in the target area at each historical monitoring moment in the historical time range, in which case, the predicted hydrological data of each monitoring point (such as one monitoring point may correspond to one flow field grid area) at each historical monitoring moment may be determined from the calibrated flow field data, and based on the difference between the predicted hydrological data of each monitoring point at each historical monitoring moment and the monitored hydrological data of each monitoring point at each historical monitoring moment, the initial hydrodynamic model may be calibrated to obtain a calibrated hydrodynamic model. Optionally, the initial hydrodynamic model may be constructed based on an Euler grid field algorithm (ie, Euler algorithm).

[0032] For example, when the sum of the predicted flow velocities of each monitoring point at each historical monitoring moment is greater than the sum of the monitored flow velocities of each monitoring point at each historical monitoring moment, the bottom friction coefficient in the initial hydrodynamic model can be increased; or, when the sum of the predicted flow velocities of each monitoring point at each historical monitoring moment is less than the sum of the monitored flow velocities of each monitoring point at each historical monitoring moment, the above-mentioned bottom friction coefficient can be reduced, and so on; the embodiments of the present invention are not limited to this.

[0033] Furthermore, the target hydrodynamic model can be determined based on the calibrated hydrodynamic model; specifically, the electronic device can verify the calibrated hydrodynamic model. If the verification result meets the model determination conditions (such as the sum of the predicted flow rates of each monitoring point at each verification historical monitoring time, and the difference between the sum of the monitored flow rates of each monitoring point at each verification historical monitoring time is less than the preset verification difference threshold, etc.), the calibrated hydrodynamic model can be used as the target hydrodynamic model; if the verification result does not meet the model determination conditions, the calibrated hydrodynamic model can continue to be calibrated until the verification result meets the model determination conditions to obtain the target hydrodynamic model. Optionally, the preset verification difference threshold can be set according to experience or according to actual needs, and the embodiment of the present invention is not limited to this.

[0034] Optionally, a hydrodynamic model may be any one of hydrodynamic models such as POM, ECOMSED, ECOM-si, FVCOM, ROMS and DELFT3D, which is not limited in the embodiment of the present invention.

[0035] Optionally, the electronic device may also include an operation status judgment module. In this case, the operation status of the target hydrodynamic model can be monitored through the operation status judgment module; when an operation error or other operation status is detected, the target hydrodynamic model can be re-run to obtain target flow field data, and so on.

[0036] The second acquisition method: the electronic device can also obtain a flow field data download link, and use the flow field data downloaded based on the flow field data download link as the target flow field data, and so on.

[0037] S102, determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutants, and record the position information of each pollutant particle among the P pollutants in the current drift diffusion simulation time step.

[0038] Wherein, M and P are both positive integers; optionally, a diffusion simulation result may include particle mass and position information of a pollutant particle, etc., which is not limited in the embodiment of the present invention.

[0039] Optionally, the current drift diffusion simulation time step may be the fth drift diffusion simulation time step, f∈[1, F], F may be the total number of drift diffusion simulation time steps (i.e., the total number of drift diffusion simulation processes), and F is a positive integer. Then when the value of f is 1, the M pollutant particles may be the initial pollutant particle set, i.e., at this time, the M pollutant particles may include all pollutant particles in the initial pollutant particle set, i.e., each initial pollutant particle may be regarded as a pollutant particle in the M pollutant particles. Optionally, the electronic device may initialize each initial pollutant particle, such as assigning an initial value to the particle mass, area, volume, proportion, and mobile weight coefficient of each initial pollutant particle; that is, the initial pollutant particle set may be initialized, and the initial pollutant particle set may be determined based on the concentration monitoring data of the target pollutant, such as determining the number and particle mass of the initial pollutant particles based on the concentration monitoring data, and the embodiment of the present invention does not limit this. Based on this, the target pollutant entering the surface water may be regarded as a set of a certain number of pollutant particles, and each pollutant particle may have a certain particle mass.

[0040] In an embodiment of the present invention, the above-mentioned determination of the particle mass of each pollutant particle among the M pollutant particles may refer to updating the particle mass of each pollutant particle among the M pollutant particles, and may also be referred to as performing a particle mass update operation on the particle mass of each pollutant particle; that is, the electronic device may perform a particle mass update operation on the particle mass of each pollutant particle.

[0041] Based on this, when the value of f is 1, before updating the particle mass of each pollutant particle, the particle mass of each pollutant particle can be the initialized particle mass, such as the particle mass of each pollutant particle can be determined based on the concentration monitoring data of the target pollutant, and at this time, each pollutant particle can be the initial pollutant particle. Then correspondingly, when the value of f is greater than 1, before executing the particle mass update operation under the f-th drift diffusion simulation time step, the particle mass of each pollutant particle can be the updated particle mass of each pollutant particle under the f-1th drift diffusion simulation time step, that is, it can be the result after executing the particle mass update operation under the f-1th drift diffusion simulation time step. It can be seen that the M pollutant particles are all pollutant particles in the initial pollutant particle set.

[0042] Then, correspondingly, the particle mass of the pollutant particles can be updated at each drift-diffusion simulation time step; optionally, the simulation time step corresponding to each drift-diffusion simulation time step can be set according to experience or according to actual needs, which is not limited in the embodiment of the present invention. Among them, one drift-diffusion simulation time step corresponds to one drift-diffusion simulation process, that is, the particle mass of the pollutant particles can be updated at each drift-diffusion simulation process.

[0043] S103 , based on the diffusion simulation result of each of the P pollutant particles, respectively updating the grid area mass concentration of each concentration determination grid area in the target area.

[0044] Optionally, each concentration determination grid area in the target area may be set based on experience or actual needs, and the embodiment of the present invention does not limit this; that is, the embodiment of the present invention does not limit the specific division method of the concentration determination grid area in the target area.

[0045] Optionally, before starting the first drift-diffusion simulation (i.e., before the first drift-diffusion simulation time step), the grid area mass concentration of each concentration determination grid area can be initialized to a preset area mass concentration, such as initialized to 0; optionally, the preset area mass concentration can be set according to experience or according to actual needs, and the embodiment of the present invention is not limited to this.

[0046] S104, when the grid area mass concentration of each concentration determination grid area has not reached the diffusion simulation stopping condition, continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stopping condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set.

[0047] At each drift-diffusion simulation time step, a drift-diffusion simulation based on the corresponding diffusion simulation results may include a particle mass update operation and a drift-diffusion simulation of the corresponding pollutant particles, that is, a drift-diffusion simulation time step may correspond to a drift-diffusion simulation process, and a drift-diffusion simulation process may include a particle mass update operation and a drift-diffusion simulation of the corresponding pollutant particles, so as to realize a particle tracing simulation (that is, a drift-diffusion simulation process) at each drift-diffusion simulation time step.

[0048] Optionally, while continuing to perform drift diffusion simulation based on the diffusion simulation results of each of the P pollutant particles, particle mass update operation and drift diffusion simulation can continue to be performed on each of the P pollutant particles until the drift diffusion simulation process at the last drift diffusion simulation time step is completed (that is, the drift diffusion simulation process when the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition). A drift diffusion simulation process may include particle mass update operation and drift diffusion simulation; it can be seen that the embodiment of the present invention can realize drift diffusion simulation of pollutant particles through a pollutant tracing model based on the Lagrangian particle algorithm.

[0049] S105: Determine the pollution impact range based on at least one position information of each initial pollutant particle.

[0050] Among them, the pollution impact range can support the loss assessment of water pollution accidents, that is, it can be used to assess the losses of water pollution accidents. Optionally, after determining the loss assessment data (also called damage assessment data, which may include the amount of loss, etc.) through the pollution impact range, the amount of engineering and funds required for the restoration of water pollution accidents can also be assessed from aspects such as the restoration of the water ecological environment baseline and the restoration of the biomass of aquatic species based on the loss assessment data, that is, the above-mentioned pollution impact range can also be used to assess the restoration data of water pollution accidents (such as the amount of engineering and funds required, etc.), etc.; the embodiment of the present invention does not limit this.

[0051] Optionally, when applying the pollution impact range to the loss assessment of water pollution accidents, the damage (i.e., loss assessment data) caused by the water pollution accident to the biological species in the target area (i.e., the main aquatic organisms) can be calculated based on the aquatic species in the target area (i.e., the biological species in the target area) and biomass, the pollution impact range and the hazard level, etc.; optionally, the hazard identification of characteristic pollutants in sudden water pollution accidents can be carried out based on the classification standards for hazardous substances to the water environment in the "Globally Harmonized System of Classification and Labelling of Chemicals" to determine the hazard level; or, the median lethal concentration or median effect concentration of pollutants for typical representative species in the accident water area (i.e., the target area) can be determined by toxicity database retrieval or toxicology testing, thereby determining the hazard level of the target pollutant to the water environment, etc.; this is not limited to the embodiments of the present invention.

[0052] Accordingly, representative ecosystem service value loss indicators can be selected from the three major service categories of supply, regulation, and culture according to the ecosystem service characteristics of biological species in the target area, so as to adopt the market value method to evaluate the value loss of supply indicators through the loss of ecological products; adopt the replacement cost method and shadow engineering method to evaluate the value loss of regulation indicators through the cost required to artificially replace natural regulation; adopt the analogy analysis method to evaluate the value loss of cultural and leisure and entertainment indicators through changes in consumer willingness to pay, etc.; based on this, all value losses caused by the current water pollution accident can be determined based on the value loss of each indicator. It can be seen that the embodiment of the present invention can apply the pollution impact range with higher accuracy to the loss assessment of water pollution accidents to improve the accuracy of loss assessment data and value loss data, thereby effectively improving the accuracy of water pollution accident loss assessment.

[0053] The embodiment of the present invention can determine the particle mass of each pollutant particle in the M pollutant particles after obtaining the target flow field data of the target area within the target time range and the target meteorological field data of the target area within the target time range; and determine the diffusion simulation result of each pollutant particle in the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, and record the position information of each pollutant particle in the P pollutant particles at the current drift diffusion simulation time step; wherein, a diffusion simulation result includes the particle mass and position information of a pollutant particle, and M and P are both positive integers. Further, based on the diffusion simulation result of each pollutant particle in the P pollutant particles, the grid area mass concentration of each concentration determination grid area in the target area can be updated respectively; when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, the drift diffusion simulation is continued based on the diffusion simulation result of each pollutant particle in the P pollutant particles, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set. Based on this, the pollution impact range can be determined based on at least one position information of each initial pollutant particle. It can be seen that the embodiment of the present invention can achieve the definition of the pollution impact range by simulating the drift and diffusion of pollutant particles, which can effectively improve the accuracy of the pollution impact range. In other words, the embodiment of the present invention can achieve high-precision definition of the impact range of sudden water pollution accidents; based on this, when the pollution impact range defined by the embodiment of the present invention is used for water pollution accident loss assessment, it can provide technical support for the accurate quantitative assessment of water ecological losses, that is, the pollution impact range with higher accuracy can be used for water pollution accident loss assessment, thereby improving the accuracy of water pollution accident loss assessment.

[0054] Based on the above description, the embodiment of the present invention also proposes a more specific method for defining the pollution impact range. Accordingly, the pollution impact range definition method can be executed by the electronic device (terminal or server) mentioned above; or, the pollution impact range definition method can be executed by the terminal and the server together. For the sake of convenience, the following description will be based on the example of an electronic device executing the pollution impact range definition method; please refer to Figure 2 The pollution impact range definition method may include the following steps S201-S208: S201, obtaining target flow field data of a target area within a target time range, and obtaining target meteorological field data of a target area within a target time range.

[0055] S202, determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutants, and record the position information of each pollutant particle among the P pollutants in the current drift diffusion simulation time step.

[0056] Optionally, when determining the particle mass of each pollutant particle among M pollutant particles, the electronic device may determine the particle impact calculation method of each pollutant particle among the M pollutant particles, and determine the particle impact data of the corresponding pollutant particles according to the particle impact calculation method of each pollutant particle, wherein one particle impact data includes the impact sub-data under each particle impact calculation method in the particle impact calculation method of the corresponding pollutant particle; based on this, the particle mass of each pollutant particle may be updated based on the particle impact data of each pollutant particle to determine the particle mass of each pollutant particle among the M pollutant particles.

[0057] Among them, the particle impact calculation method of a pollutant particle may include at least one of the following: a volatilization calculation method, a degradation calculation method, and an evaporation calculation method, etc., which is not limited in the embodiment of the present invention. In the embodiment of the present invention, for any pollutant particle among the M pollutant particles, the particle impact calculation method of any pollutant particle can be determined based on the position information of any pollutant particle; specifically, a position information may include the depth (i.e., the vertical position) and horizontal position of the corresponding pollutant particle, etc. In this case, the particle impact calculation method of any pollutant particle can be determined based on the depth of any pollutant particle.

[0058] Optionally, the electronic device may determine that any pollutant particle is located on the water surface or in the water body based on the depth of any pollutant particle; if any pollutant particle is located on the water surface, the particle impact calculation method for determining any pollutant particle may include but is not limited to a volatilization calculation method, an evaporation calculation method, and a degradation calculation method; if any pollutant particle is located in the water body, the particle impact calculation method for determining any pollutant particle may include but is not limited to a degradation calculation method. Specifically, the electronic device may determine that any pollutant particle is located on the water surface or in the water body based on the depth of any pollutant particle and the water surface depth threshold; if the depth of any pollutant particle is less than or equal to the water surface depth threshold, it may be determined that any pollutant particle is located on the water surface; if the depth of any pollutant particle is greater than the water surface depth threshold, it may be determined that any pollutant particle is located in the water body. Optionally, the water surface depth threshold may be determined based on the simulation time step and the wave height, or may be set according to experience, or may be set according to actual needs, and the embodiment of the present invention does not limit this; optionally, the wave height may be set according to experience, or may be set according to actual needs, and the embodiment of the present invention does not limit this. Optionally, if the water meter depth threshold is determined based on the simulation time step and wave height, the electronic device may use formula 1.1 to calculate the water meter depth threshold: Formula 1.1 Among them, D z can refer to the wave height and ∆t can refer to the simulation time step.

[0059] It should be noted that the particle impact data of a pollutant particle may be the particle impact data of the corresponding pollutant particle at the current drift diffusion simulation time step (i.e., the fth drift diffusion simulation time step), that is, the electronic device may respectively determine the particle impact data of each pollutant particle at the current drift diffusion simulation time step. Based on this, the particle impact data of a pollutant particle may include but is not limited to at least one of the following: the volatilization mass, degradation mass, and evaporation mass of the corresponding pollutant particle.

[0060] Specifically, when determining the particle impact data of the corresponding pollutant particles according to the particle impact calculation method of each pollutant particle, for any pollutant particle among the M pollutant particles, and any particle impact calculation method in the particle impact calculation method of any pollutant particle, if any particle impact calculation method is a volatilization calculation method, the impact sub-data of any pollutant particle under the volatilization calculation method (i.e., the volatilization mass of any pollutant particle) can be determined according to the volatilization calculation method, and the particle impact data of any pollutant particle includes the impact sub-data of any pollutant particle under any particle impact calculation method. Optionally, the electronic device can use formula 1.2 to calculate the volatilization mass of any pollutant particle: Formula 1.2 Among them, d m is the volatile mass of any pollutant particle in the current drift diffusion simulation time step, M0 is the particle mass of any pollutant particle, d is the above-mentioned water surface depth threshold, and K7 can be the mass transfer coefficient; optionally, one pollutant can correspond to one mass transfer coefficient, and the mass transfer coefficients corresponding to different pollutants can be set according to experience or according to actual needs; based on this, K7 can be the mass transfer coefficient corresponding to the target pollutant.

[0061] Correspondingly, if the impact calculation method of any particle is a degradation calculation method, then the electronic device can determine the impact sub-data of any pollutant particle under the degradation calculation method according to the degradation calculation method, that is, the degradation mass of any pollutant particle; optionally, the electronic device can use formula 1.3 to calculate the degradation result of any pollutant particle on the tth day: Formula 1.3 Among them, M t It can be the degradation result of any pollutant particle on the tth day, k can be the decay rate of the substance, and t can be the number of days in the current drift diffusion simulation time step (for example, if the current drift diffusion simulation time step is a drift diffusion simulation time step within the first day of the target time range, then the value of t can be 1). Optionally, the degradation mass of any pollutant particle (i.e., the degradation mass of any pollutant particle under the current drift diffusion simulation time step) can be (the particle mass of any pollutant particle - the degradation result of any pollutant particle on the tth day) / the total number of drift diffusion simulation time steps in one day, where the particle mass of any pollutant particle refers to the updated particle mass of any pollutant particle under the f-1th drift diffusion simulation time step when the current drift diffusion simulation time step is the fth drift diffusion simulation time step; or, it can be (the updated particle mass or the initialization particle mass of any pollutant particle under the last drift diffusion simulation time step on the t-1th day - the degradation result of any pollutant particle on the tth day) / the total number of drift diffusion simulation time steps in one day, etc., which is not limited in the embodiments of the present invention. Optionally, one pollutant may correspond to one decay rate, and the decay rates of different pollutants may be set based on experience or actual needs, which is not limited in the embodiment of the present invention.

[0062] Then, accordingly, if the impact calculation method of any particle is the evaporation calculation method, then the electronic device can determine the impact sub-data of any pollutant particle under the evaporation calculation method according to the evaporation calculation method, that is, the evaporation mass of any pollutant particle; optionally, the electronic device can calculate the evaporation rate of any pollutant particle, and then use the evaporation rate of any pollutant particle to calculate the evaporation mass of any pollutant particle, that is, the multiplication result between the evaporation rate of any pollutant particle and the simulation time step can be used as the evaporation mass of any pollutant particle. Specifically, the electronic device can use formula 1.4 to calculate the evaporation rate of any pollutant particle: Formula 1.4 Where E is the evaporation rate of any pollutant particle (g / m 2 h, i.e. grams per square meter per hour), K e M can be the evaporation rate (m / h), W It can be expressed as molar mass (g / mole), P VP K may be the saturated vapor pressure corresponding to the target pollutant, R may be the gas constant, and T may be the water temperature. Optionally, the electronic device may use formula 1.5 to calculate K e : Formula 1.5 Among them, U can be the wind speed at 10m (m / h, i.e. meters per hour), that is, the wind speed 10 meters above the water surface; L can be the diameter of the cylinder or sphere corresponding to any pollutant particle in any mass state (m); S C It can be the Schmidt number (i.e. critical point value) of the target pollutant in the air; optionally, one pollutant can correspond to one S C , S of different pollutants C It can be set according to experience or according to actual needs, and the embodiment of the present invention does not limit this.

[0063] Further, for any of the M pollutant particles, if the particle impact data of any pollutant particle includes the volatilization mass of any pollutant particle, the volatilization mass of any pollutant particle may be subtracted from the particle mass of any pollutant particle to update the particle mass of any pollutant particle. If the particle impact data of any pollutant particle includes the evaporation mass of any pollutant particle, the evaporation mass of any pollutant particle may be subtracted from the particle mass of any pollutant particle to update the particle mass of any pollutant particle. If the particle impact data of any pollutant particle includes the degradation mass of any pollutant particle, the degradation mass of any pollutant particle may be subtracted from the particle mass of any pollutant particle to update the particle mass of any pollutant particle. In other words, the particle mass of any pollutant particle may be subtracted from each influencing sub-data in the particle impact data of any pollutant particle to update the particle mass of each pollutant particle based on the particle impact data of each pollutant particle.

[0064] In other embodiments, the electronic device may also update the particle mass of any pollutant particle immediately according to the influence sub-data of any pollutant particle under any particle impact calculation mode after obtaining the influence sub-data of any pollutant particle under any particle impact calculation mode (such as the volatile mass of any pollutant particle), thereby continuing to calculate the influence sub-data of any pollutant particle under the next particle impact calculation mode according to the updated particle mass of any pollutant particle, thereby further updating the particle mass of any pollutant particle, and so on; the present invention is not limited to this.

[0065] Optionally, when determining the diffusion simulation result of each of the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data, and the target meteorological field data, an active mass threshold can be determined, and based on the particle mass of each pollutant particle, P pollutant particles with a particle mass greater than the active mass threshold can be determined from the M pollutant particles; based on this, drift diffusion simulation can be performed on each of the P pollutant particles based on the target flow field data and the target meteorological field data, respectively, to obtain the diffusion simulation result of each of the P pollutant particles. Optionally, the active mass threshold can be set according to experience or according to actual needs, and the embodiment of the present invention is not limited to this.

[0066] Optionally, when drift diffusion simulation is performed on each of the P pollutant particles based on the target flow field data and the target meteorological field data to obtain the diffusion simulation result of each of the P pollutant particles, for any of the P pollutant particles, the electronic device can capture the surrounding flow velocity data and wind speed data of any pollutant particle in the current drift diffusion simulation time step in real time based on the target flow field data and the target meteorological field data, so as to obtain the flow field information of each of the at least one adjacent flow field grid area corresponding to any pollutant particle, and obtain the meteorological field information (such as wind field information) of each of the at least one adjacent meteorological field grid area corresponding to any pollutant particle, so as to perform weighted summation on the flow field information of each adjacent flow field grid area to obtain the flow field information of any pollutant particle (such as the flow velocity in each direction in at least one direction), and perform weighted summation on the meteorological field information of each adjacent meteorological field grid area to obtain the meteorological field information of any pollutant particle (such as the wind speed in each direction in at least one direction). Based on this, drift diffusion simulation can be performed on any pollutant particle based on the flow field information and meteorological field information of any pollutant particle.

[0067] Optionally, an adjacent flow field grid area of ​​any pollutant particle may be: a flow field grid area whose distance to any pollutant particle is less than a preset flow field distance threshold, or the first Y flow field grid areas closest to any pollutant particle, Y is a positive integer, and so on; this is not limited in the embodiment of the present invention. Correspondingly, an adjacent meteorological field grid area of ​​any pollutant particle may be: a meteorological field grid area whose distance to any pollutant particle is less than a preset meteorological field distance threshold, or the first V meteorological field grid areas closest to any pollutant particle, V is a positive integer, and so on; this is not limited in the embodiment of the present invention. Optionally, the preset flow field distance threshold and the preset meteorological field distance threshold may be set according to experience or according to actual needs, and this is not limited in the embodiment of the present invention.

[0068] Optionally, the flow field information of any adjacent flow field grid area can be determined based on the flow field information of any adjacent flow field grid area at the first flow field time step (i.e., the previous flow field time step) and the second flow field time step (i.e., the next flow field time step) corresponding to the current drift-diffusion simulation time step; specifically, the flow field information of any adjacent flow field grid area at the first flow field time step and the flow field information at the second flow field time step can be weighted summed according to the spatial distance to obtain the flow field information of any adjacent flow field grid area. Optionally, the weight of the flow field information of any adjacent flow field grid area at the first flow field time step can be: the ratio between the first distance and the time step distance, the first distance can be the distance between the current drift-diffusion simulation time step and the first flow field time step, and the time step distance can be the distance between the first flow field time step and the second flow field time step; correspondingly, the weight of the flow field information of any adjacent flow field grid area at the second flow field time step can be: the ratio between the second distance and the time step distance, and the second distance can be the distance between the current drift-diffusion simulation time step and the second flow field time step. Similarly, the meteorological field information of any adjacent meteorological field grid area can be determined based on the meteorological field information of any adjacent meteorological field grid area at the first meteorological field time step and the second meteorological field time step corresponding to the current drift diffusion simulation time step. The method for determining the meteorological field information of any adjacent meteorological field grid area can be the same as the method for determining the flow field information of any adjacent flow field grid area, and the embodiments of the present invention will not be repeated here.

[0069] Optionally, when weighted summing is performed on the flow field information of each adjacent flow field grid area to obtain the flow field information of any pollutant particle, the weights of the flow field information of each adjacent flow field grid area may be the same (such as mean operation) or different, and the embodiment of the present invention does not limit this. Exemplarily, the electronic device may also determine the weights of the flow field information of each adjacent flow field grid area based on the distance between each adjacent flow field grid area and any pollutant particle; taking any adjacent flow field grid area in at least one adjacent flow field grid area as an example, the weight of the flow field information of any adjacent flow field grid area may be the ratio between the third distance and the flow field grid distance, the third distance may be the distance between any adjacent flow field grid area and any pollutant particle, and the flow field grid distance may be the sum of the distances between each adjacent flow field grid area and any pollutant particle. Correspondingly, the weights of the meteorological field information of each adjacent meteorological field grid area may be the same or different, and the embodiment of the present invention does not limit this; wherein, the specific determination method of the weights of the meteorological field information of each adjacent meteorological field grid area may be the same as the specific determination method of the weights of the flow field information of each adjacent flow field grid area, and the embodiment of the present invention will not be repeated here.

[0070] Optionally, the electronic device may also be assisted by a fixed shoreline to limit the drift and diffusion simulation of pollutant particles to always be within the target area, so as to realize the drift and diffusion simulation process of pollutant particles.

[0071] S203 , based on the diffusion simulation result of each of the P pollutant particles, respectively updating the grid area mass concentration of each concentration determination grid area in the target area.

[0072] Specifically, for any concentration determination grid area in the target area, the electronic device can determine each pollutant particle in any concentration determination grid area from the P pollutant particles based on the position information of each pollutant particle among the P pollutant particles; and update the grid area mass concentration of any concentration determination grid area based on the particle mass of each pollutant particle in any concentration determination grid area.

[0073] Optionally, the particle mass of each pollutant particle in any concentration determination grid area can be summed to obtain the total mass of particles in any concentration determination grid area, and the total mass of particles in any concentration determination grid area and the area volume of any concentration determination grid area can be divided to obtain the division result of any concentration determination grid area, so as to update the grid area mass concentration of any concentration determination grid area using the division result of any concentration determination grid area; or, when the area volumes of each concentration determination grid area are unit volumes, the total mass of particles in any concentration determination grid area can be used to update the grid area mass concentration of any concentration determination grid area, and so on; this is not limited to the embodiments of the present invention.

[0074] It can be seen that the embodiment of the present invention can use the diffusion simulation results of each pollutant particle (such as each of the P pollutant particles mentioned above) that survives the current drift-diffusion simulation time step to update the grid area mass concentration of each concentration determination grid area respectively; that is, the embodiment of the present invention can use the diffusion simulation results of each pollutant particle that survives each drift-diffusion simulation time step to update the grid area mass concentration of each concentration determination grid area respectively after each drift-diffusion simulation. Among them, the pollutant particles that survive a drift-diffusion simulation time step may include pollutant particles whose particle mass is greater than the active mass threshold value in the corresponding drift-diffusion simulation time step.

[0075] S204, determining the maximum tolerance concentration of biological species in the target area for the target pollutant, and determining the mass concentration of the grid area based on the grid area mass concentration of each concentration to determine the mass concentration of the area to be determined.

[0076] Optionally, the above maximum tolerance concentration may be set according to experience or according to actual needs, which is not limited in the embodiment of the present invention.

[0077] In one implementation, the electronic device may use the maximum value among the mass concentrations of the grid areas of each concentration determination grid area as the mass concentration of the area to be determined.

[0078] In another embodiment, the electronic device may screen out Q concentration determination grid areas from the target area based on the grid area mass concentrations of each concentration determination grid area, and use the average of the grid area mass concentrations of each concentration determination grid area in the Q concentration determination grid areas as the mass concentration of the area to be determined; wherein the grid area mass concentration of any concentration determination grid area among the Q concentration determination grid areas is greater than the grid area mass concentration of any concentration determination grid area in the target area except the Q concentration determination grid areas, and Q is a positive integer.

[0079] Optionally, the Q concentration determination grid areas may include all concentration determination grid areas in the target area whose grid area mass concentration is greater than a preset grid area mass concentration threshold, that is, all concentration determination grid areas in the target area whose grid area mass concentration is greater than a preset grid area mass concentration threshold may be used as concentration determination grid areas in the Q concentration determination grid areas; or, the concentration determination grid areas may be sorted in descending order according to the grid area mass concentration of each concentration determination grid area to obtain a concentration determination grid area sorting result, and the first Q concentration determination grid areas may be selected from the concentration determination grid area sorting result to select Q concentration determination grid areas from the target area, and so on; the embodiment of the present invention is not limited to this. Optionally, the preset grid area mass concentration threshold may be set according to experience or according to actual needs, and the embodiment of the present invention is not limited to this.

[0080] S205, if the mass concentration of the area to be determined is less than or equal to the maximum tolerance concentration, determining whether the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition.

[0081] S206: If the mass concentration of the area to be determined is greater than the maximum tolerance concentration, it is determined that the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition.

[0082] S207, when the grid area mass concentration of each concentration determination grid area has not reached the diffusion simulation stop condition, continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set.

[0083] Among them, the position information of an initial pollutant particle can be the position information of the corresponding initial pollutant particle in a drift diffusion simulation time step. It should be understood that the pollutant particles that survive in any drift diffusion simulation time step (such as M pollutant particles and P pollutant particles, etc.) are all pollutant particles in the initial pollutant particle set, that is, they are all a set of pollutant particles screened out from the initial pollutant particle set in each drift diffusion simulation time step, that is, any pollutant particle is an initial pollutant particle, that is, a pollutant particle in M ​​pollutant particles and a pollutant particle in P pollutant particles can be the same initial pollutant particle, such as the same initial pollutant particle with different diffusion simulation results; optionally, the diffusion simulation results of the same initial pollutant particle in different drift diffusion simulation time steps can be different.

[0084] In an embodiment of the present invention, the electronic device may record the position information of each initial pollutant particle at each drift-diffusion simulation time step, that is, the position information of pollutant particles that survive at each drift-diffusion simulation time step; optionally, when the particle mass of any pollutant particle (i.e., any initial pollutant particle) is less than or equal to the active mass threshold, the position information of any pollutant particle at subsequent drift-diffusion simulation time steps may remain unchanged or be empty, etc., and the embodiment of the present invention does not limit this. Based on this, at least one position information of any initial pollutant particle may include the position information of any initial pollutant particle at each drift-diffusion simulation time step in at least one drift-diffusion simulation time step; optionally, at least one drift-diffusion simulation time step may include all drift-diffusion simulation time steps when the particle mass of any initial pollutant particle is greater than the active mass threshold, and may also include all drift-diffusion simulation time steps in the entire simulation process, etc., and the embodiment of the present invention does not limit this.

[0085] S208: Determine the pollution impact range based on at least one position information of each initial pollutant particle.

[0086] In one embodiment, the electronic device may determine the pollution impact envelope based on at least one position information of each initial pollutant particle; and determine the pollution impact range based on the pollution impact envelope. Optionally, the electronic device may solve the envelope of all initial pollutant particles through a convhull function (a convex hull function) to achieve the determination of the pollution impact envelope. Accordingly, the polygon inside the pollution impact envelope may be divided into a number of trapezoids for integral solution, thereby achieving the determination of the pollution impact range based on the pollution impact envelope.

[0087] In another embodiment, the electronic device can count the number of times particles pass through each concentration determination grid area based on at least one position information of each initial pollutant particle. The number of times particles pass through a concentration determination grid area can be equal to the number of particles located in the corresponding concentration determination grid area in at least one position information of each initial pollutant particle, that is, the sum of the number of times each initial pollutant particle passes through the corresponding concentration determination grid area; in this case, all concentration determination grid areas with a particle passing number greater than a preset particle passing number threshold can be selected from the target area, and the sum of the area volumes between all selected concentration determination grid areas is used as the pollution impact range. Optionally, the preset particle passing number threshold can be set according to experience or according to actual needs, and the embodiment of the present invention does not limit this.

[0088] In another embodiment, the electronic device may record the grid area mass concentration of each concentration determination grid area at the current drift diffusion simulation time step, that is, after each update of the grid area mass concentration of each concentration determination grid area, the grid area mass concentration of each concentration determination grid area is recorded, thereby obtaining the grid area mass concentration recording data of each concentration determination grid area, and the grid area mass concentration recording data of a concentration determination grid area may include the grid area mass concentration of the corresponding concentration determination grid area at each drift diffusion simulation time step; in this case, based on the grid area mass concentration recording data of each concentration determination grid area, all concentration determination grid areas in the target area where the grid area mass concentration is greater than the preset concentration selection threshold in the grid area mass concentration recording data can be determined, and the sum of the area volumes between all the determined concentration determination grid areas is used as the pollution impact range, etc. Optionally, the preset concentration selection threshold can be set according to experience or according to actual needs, and the embodiment of the present invention is not limited to this.

[0089] The embodiment of the present invention can determine the particle mass of each pollutant particle in the M pollutant particles after obtaining the target flow field data of the target area within the target time range and the target meteorological field data of the target area within the target time range; and determine the diffusion simulation result of each pollutant particle in the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, and record the position information of each pollutant particle in the P pollutant particles at the current drift diffusion simulation time step; and can update the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each pollutant particle in the P pollutant particles. Based on this, the maximum tolerance concentration of biological species in the target area for the target pollutant can be determined, and the mass concentration of the area to be determined can be determined based on the grid area mass concentration of each concentration determination grid area; if the mass concentration of the area to be determined is less than or equal to the maximum tolerance concentration, it is determined that the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition; if the mass concentration of the area to be determined is greater than the maximum tolerance concentration, it is determined that the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition. Then correspondingly, when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, the drift diffusion simulation is continued based on the diffusion simulation result of each pollutant particle in the P pollutant particles, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; based on this, the pollution impact range can be determined based on at least one position information of each initial pollutant particle. It can be seen that the embodiment of the present invention can achieve high-precision definition of the impact range of sudden water pollution accidents through the target hydrodynamic model and the pollutant tracing model based on the Lagrangian particle algorithm; that is, the embodiment of the present invention can achieve the definition of the pollution impact range through the drift diffusion simulation of pollutant particles, which can effectively improve the accuracy of the pollution impact range.

[0090] Based on the description of the relevant embodiments of the above-mentioned pollution impact range definition method, the embodiment of the present invention further proposes a pollution impact range definition device, which can be a computer program (including program code) running in an electronic device; Figure 3 As shown, the pollution impact range defining device may include an acquisition unit 301 and a processing unit 302. The pollution impact range defining device may execute Figure 1 or Figure 2 The pollution impact scope definition method shown, that is, the pollution impact scope definition device can run the above units: An acquisition unit 301 is used to acquire target flow field data of a target area within a target time range, and to acquire target meteorological field data of the target area within the target time range; The processing unit 302 is used to determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutant particles, and record the position information of each pollutant particle among the P pollutant particles at the current drift diffusion simulation time step; wherein one diffusion simulation result includes the particle mass and position information of one pollutant particle, and M and P are both positive integers; The processing unit 302 is further configured to update the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each of the P pollutant particles; The processing unit 302 is further configured to continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; The processing unit 302 is further configured to determine a pollution impact range based on at least one position information of each initial pollutant particle.

[0091] In one embodiment, when the processing unit 302 determines the diffusion simulation result of each of the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data, and the target meteorological field data, it can be specifically used to: Determining an active mass threshold, and based on the particle mass of each of the pollutant particles, determining P pollutant particles having a particle mass greater than the active mass threshold from the M pollutant particles; Based on the target flow field data and the target meteorological field data, drift diffusion simulation is performed on each of the P pollutant particles to obtain a diffusion simulation result of each of the P pollutant particles.

[0092] In another embodiment, when the processing unit 302 updates the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each of the P pollutant particles, it can be specifically used to: For any concentration determination grid area in the target area, based on the position information of each pollutant particle among the P pollutant particles, determine each pollutant particle in the any concentration determination grid area from the P pollutant particles; Based on the particle mass of each pollutant particle in any concentration determination grid area, the grid area mass concentration of any concentration determination grid area is updated.

[0093] In another implementation, the processing unit 302 may also be used to: Determine the maximum tolerance concentration of biological species in the target area for the target pollutant, and determine the mass concentration of the area to be determined based on the grid area mass concentration of each concentration determination grid area; If the mass concentration of the area to be determined is less than or equal to the maximum tolerance concentration, determining that the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition; If the mass concentration of the to-be-determined area is greater than the maximum tolerated concentration, it is determined that the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition.

[0094] In another embodiment, when the processing unit 302 determines the mass concentration of the area to be determined based on the grid area mass concentrations of each concentration determination grid area, it can be specifically used to: The maximum value among the mass concentrations of the grid areas of each concentration determination grid area is used as the mass concentration of the area to be determined; or Based on the grid area mass concentrations of the various concentration determination grid areas, Q concentration determination grid areas are screened out from the target area, and the average of the grid area mass concentrations of the various concentration determination grid areas in the Q concentration determination grid areas is used as the mass concentration of the area to be determined; wherein, the grid area mass concentration of any concentration determination grid area among the Q concentration determination grid areas is greater than the grid area mass concentration of any concentration determination grid area in the target area except the Q concentration determination grid areas, and Q is a positive integer.

[0095] In another embodiment, when determining the pollution impact range based on the at least one position information of each initial pollutant particle, the processing unit 302 may be specifically configured to: Determining a pollution impact envelope based on at least one position information of each initial pollutant particle; Based on the pollution impact envelope, a pollution impact range is determined.

[0096] In another embodiment, when determining the particle mass of each of the M pollutant particles, the processing unit 302 may be specifically configured to: Determine a particle impact calculation method for each of the M pollutant particles, and determine particle impact data of the corresponding pollutant particles according to the particle impact calculation method of each of the pollutant particles, wherein one particle impact data includes impact sub-data under each particle impact calculation method in the particle impact calculation method of the corresponding pollutant particle; Based on the particle impact data of each pollutant particle, the particle mass of each pollutant particle is updated to achieve the determination of the particle mass of each pollutant particle in the M pollutant particles.

[0097] According to one embodiment of the present invention, Figure 1 or Figure 2 Each step involved in the method shown can be Figure 3 The pollution impact scoping shown is performed by each unit in the device. For example, Figure 1 The step S101 shown in FIG. 1 can be performed by Figure 3 The acquisition unit 301 shown in FIG. 1 is executed, and steps S102-S105 can be performed by Figure 3 The processing unit 302 shown in is executed. For example, Figure 2 The step S201 shown in FIG. 1 can be performed by Figure 3 The acquisition unit 301 shown in FIG. 1 is executed, and steps S202-S208 can be performed by Figure 3 The processing unit 302 shown in executes, and so on.

[0098] According to another embodiment of the present invention, Figure 3 Each unit in the pollution impact range definition device shown can be separately or completely combined into one or several other units to form a unit, or one (or some) of the units can be further divided into multiple functionally smaller units to form a unit, which can achieve the same operation without affecting the realization of the technical effects of the embodiments of the present invention. The above-mentioned units are divided based on logical functions. In practical applications, the function of one unit can also be implemented by multiple units, or the functions of multiple units can be implemented by one unit. In other embodiments of the present invention, any pollution impact range definition device may also include other units. In practical applications, these functions can also be implemented with the assistance of other units, and can be implemented by the collaboration of multiple units.

[0099] According to another embodiment of the present invention, the program can be executed by running a program on a general electronic device such as a computer including a central processing unit (CPU), a random access memory medium (RAM), a read-only memory medium (ROM), and other processing elements and storage elements. Figure 1 or Figure 2 A computer program (including program code) for each step involved in the corresponding method shown in Figure 3The pollution impact scope defining device shown in the embodiment of the present invention is used to implement the pollution impact scope defining method of the embodiment of the present invention. The computer program can be recorded on a computer storage medium, for example, and loaded into the above-mentioned electronic device through the computer storage medium and run therein.

[0100] Based on the description of the above method embodiment and device embodiment, the exemplary embodiment of the present invention further provides an electronic device, including: at least one processor; and a memory connected to the at least one processor in communication. The memory stores a computer program that can be executed by the at least one processor, and the computer program is used to enable the electronic device to perform the method according to the embodiment of the present invention when executed by the at least one processor.

[0101] Exemplary embodiments of the present invention also provide a non-transitory computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor of a computer, is used to cause the computer to perform a method according to an embodiment of the present invention.

[0102] An exemplary embodiment of the present invention further provides a computer program product, comprising a computer program, wherein when the computer program is executed by a processor of a computer, the computer is used to enable the computer to perform a method according to an embodiment of the present invention.

[0103] refer to Figure 4 , a block diagram of an electronic device 400 that can be used as a server or client of the present invention will now be described, which is an example of a hardware device that can be applied to various aspects of the present invention. The electronic device is intended to represent various forms of digital electronic computer devices, such as laptop computers, desktop computers, workbenches, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present invention described and / or required herein.

[0104] like Figure 4 As shown, the electronic device 400 includes a computing unit 401, which can perform various appropriate actions and processes according to a computer program stored in a read-only memory (ROM) 402 or a computer program loaded from a storage unit 408 into a random access memory (RAM) 403. In the RAM 403, various programs and data required for the operation of the electronic device 400 can also be stored. The computing unit 401, the ROM 402, and the RAM 403 are connected to each other via a bus 404. An input / output (I / O) interface 405 is also connected to the bus 404.

[0105] A plurality of components in the electronic device 400 are connected to the I / O interface 405, including: an input unit 406, an output unit 407, a storage unit 408, and a communication unit 409. The input unit 406 may be any type of device capable of inputting information to the electronic device 400, and the input unit 406 may receive input digital or character information, and generate key signal inputs related to user settings and / or function control of the electronic device. The output unit 407 may be any type of device capable of presenting information, and may include, but is not limited to, a display, a speaker, a video / audio output terminal, a vibrator, and / or a printer. The storage unit 408 may include, but is not limited to, a disk, an optical disk. The communication unit 409 allows the electronic device 400 to exchange information / data with other devices via a computer network such as the Internet and / or various telecommunication networks, and may include, but is not limited to, a modem, a network card, an infrared communication device, a wireless communication transceiver, and / or a chipset, such as a Bluetooth™ device, a WiFi device, a WiMax device, a cellular communication device, and / or the like.

[0106] The computing unit 401 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the computing unit 401 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various computing units running machine learning model algorithms, digital signal processors (DSPs), and any appropriate processors, controllers, microcontrollers, etc. The computing unit 401 performs the various methods and processes described above. For example, in some embodiments, the pollution impact scope definition method may be implemented as a computer software program, which is tangibly contained in a machine-readable medium, such as a storage unit 408. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 400 via the ROM 402 and / or the communication unit 409. In some embodiments, the computing unit 401 may be configured to perform the pollution impact scope definition method in any other appropriate manner (e.g., by means of firmware).

[0107] The program code for implementing the method of the present invention can be written in any combination of one or more programming languages. These program codes can be provided to a processor or controller of a general-purpose computer, a special-purpose computer or other programmable data processing device, so that the program code, when executed by the processor or controller, enables the functions / operations specified in the flow chart and / or block diagram to be implemented. The program code can be executed entirely on the machine, partially on the machine, partially on the machine as a stand-alone software package and partially on a remote machine, or entirely on a remote machine or server.

[0108] In the context of the present invention, a machine-readable medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, device, or equipment. A machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, device, or equipment, or any suitable combination of the foregoing. A more specific example of a machine-readable storage medium may include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0109] As used herein, the terms "machine-readable medium" and "computer-readable medium" refer to any computer program product, apparatus, and / or device (e.g., disk, optical disk, memory, programmable logic device (PLD)) for providing machine instructions and / or data to a programmable processor, including a machine-readable medium that receives machine instructions as a machine-readable signal. The term "machine-readable signal" refers to any signal for providing machine instructions and / or data to a programmable processor.

[0110] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the computer. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0111] The systems and techniques described herein may be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer with a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system may be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include: a local area network (LAN), a wide area network (WAN), and the Internet.

[0112] A computer system may include clients and servers. Clients and servers are generally remote from each other and usually interact through a communication network. The relationship of client and server is generated by computer programs running on respective computers and having a client-server relationship to each other.

[0113] Furthermore, it should be understood that what is disclosed above is only a preferred embodiment of the present invention, and certainly cannot be used to limit the scope of rights of the present invention. Therefore, equivalent changes made according to the claims of the present invention are still within the scope of the present invention.

Claims

1. A method for defining the scope of pollution impact, characterized in that: include: Acquire target flow field data of a target area within a target time range, and acquire target meteorological field data of the target area within the target time range; Determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutant particles, and record the position information of each pollutant particle among the P pollutant particles at the current drift diffusion simulation time step; wherein one diffusion simulation result includes the particle mass and position information of one pollutant particle, and M and P are both positive integers; Based on the diffusion simulation result of each of the P pollutant particles, respectively updating the grid area mass concentration of each concentration determination grid area in the target area; When the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; Based on the at least one position information of each initial pollutant particle, a pollution impact range is determined.

2. The method according to claim 1, characterized in that The step of determining the diffusion simulation result of each of the P pollutant particles based on the particle mass of each pollutant particle, the target flow field data, and the target meteorological field data includes: Determining an active mass threshold, and based on the particle mass of each of the pollutant particles, determining P pollutant particles having a particle mass greater than the active mass threshold from the M pollutant particles; Based on the target flow field data and the target meteorological field data, drift diffusion simulation is performed on each of the P pollutant particles to obtain a diffusion simulation result of each of the P pollutant particles.

3. The method according to claim 1 or 2, characterized in that: The updating of the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each of the P pollutant particles comprises: For any concentration determination grid area in the target area, based on the position information of each pollutant particle among the P pollutant particles, determine each pollutant particle in the any concentration determination grid area from the P pollutant particles; Based on the particle mass of each pollutant particle in any concentration determination grid area, the grid area mass concentration of any concentration determination grid area is updated.

4. The method according to claim 1 or 2, characterized in that: The method further comprises: Determine the maximum tolerance concentration of biological species in the target area for the target pollutant, and determine the mass concentration of the area to be determined based on the grid area mass concentration of each concentration determination grid area; If the mass concentration of the area to be determined is less than or equal to the maximum tolerance concentration, determining that the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition; If the mass concentration of the to-be-determined area is greater than the maximum tolerated concentration, it is determined that the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition.

5. The method according to claim 4, characterized in that The determining the mass concentration of the area to be determined based on the mass concentration of the grid area of ​​each concentration determination grid area comprises: The maximum value among the mass concentrations of the grid areas of each concentration determination grid area is used as the mass concentration of the area to be determined; or Based on the grid area mass concentrations of the various concentration determination grid areas, Q concentration determination grid areas are screened out from the target area, and the average of the grid area mass concentrations of the various concentration determination grid areas in the Q concentration determination grid areas is used as the mass concentration of the area to be determined; wherein, the grid area mass concentration of any concentration determination grid area among the Q concentration determination grid areas is greater than the grid area mass concentration of any concentration determination grid area in the target area except the Q concentration determination grid areas, and Q is a positive integer.

6. The method according to claim 1 or 2, characterized in that: The determining of the pollution impact range based on at least one position information of each initial pollutant particle includes: Determining a pollution impact envelope based on at least one position information of each initial pollutant particle; Based on the pollution impact envelope, a pollution impact range is determined.

7. The method according to claim 1 or 2, characterized in that: The determining of the particle mass of each of the M pollutant particles comprises: Determine a particle impact calculation method for each of the M pollutant particles, and determine particle impact data of the corresponding pollutant particles according to the particle impact calculation method of each of the pollutant particles, wherein one particle impact data includes impact sub-data under each particle impact calculation method in the particle impact calculation method of the corresponding pollutant particle; Based on the particle impact data of each pollutant particle, the particle mass of each pollutant particle is updated to achieve the determination of the particle mass of each pollutant particle in the M pollutant particles.

8. A pollution impact range defining device, characterized in that: The device comprises: An acquisition unit, used for acquiring target flow field data of a target area within a target time range, and acquiring target meteorological field data of the target area within the target time range; A processing unit, used to determine the particle mass of each pollutant particle among the M pollutant particles; and based on the particle mass of each pollutant particle, the target flow field data and the target meteorological field data, determine the diffusion simulation result of each pollutant particle among the P pollutant particles, and record the position information of each pollutant particle among the P pollutant particles at the current drift diffusion simulation time step; wherein one diffusion simulation result includes the particle mass and position information of one pollutant particle, and M and P are both positive integers; The processing unit is further used to update the grid area mass concentration of each concentration determination grid area in the target area based on the diffusion simulation result of each pollutant particle in the P pollutant particles; The processing unit is further configured to continue to perform drift diffusion simulation based on the diffusion simulation result of each pollutant particle in the P pollutant particles when the grid area mass concentration of each concentration determination grid area does not reach the diffusion simulation stop condition, until the grid area mass concentration of each concentration determination grid area reaches the diffusion simulation stop condition, thereby obtaining at least one position information of each initial pollutant particle in the initial pollutant particle set, and the M pollutant particles are all pollutant particles in the initial pollutant particle set; The processing unit is further used to determine the pollution impact range based on at least one position information of each initial pollutant particle.

9. An electronic device, characterized in that: include: processor; as well as Memory for storing programs, The program includes instructions, which, when executed by the processor, cause the processor to perform the method according to any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium storing computer instructions, wherein: The computer instructions are used to make a computer execute the method according to any one of claims 1-7.

Images