Method and device for tracing substances in sea area, computer equipment and storage medium

Abstract

The invention provides a method and device for tracing substances in a sea area, computer equipment and a storage medium, and belongs to the field of sea area pollution tracing. The method comprises the following steps: pre-constructing a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field according to position information of a sea area to be verified; setting at least two observation stations in a sea area to be verified, and obtaining a monitoring substance concentration time sequence of each observation station; determining starting and ending moments of a reverse flow field according to the monitoring substance concentration time sequence, and solving an adjoint equation according to the reverse flow field, the starting and ending moments and the position data of the observation station to obtain an adjoint substance concentration space distribution state and an adjoint substance concentration time sequence; and according to the monitoring substance concentration time sequence, the accompanying substance concentration space distribution state and the accompanying substance concentration time sequence, determining the position information of the target pollution source, the release intensity of pollutants released by the target pollution source and the change trend of the release intensity. The method can achieve the effect of improving the accuracy and efficiency of marine pollutant traceability.

Description

Technical Field

[0001] The present application relates to the field of tracing the source of sea area pollution. Specifically, it relates to a method, a device, a computer device, and a storage medium for tracing substances in a sea area. Background Art

[0002] As one of the environmental protection topics that have attracted much attention nowadays, the water safety problem in coastal areas is often affected by various human activities, resulting in a series of water environment problems. Human activities such as sewage discharge and ship oil spill will seriously threaten the balance of the marine ecological environment. In addition, the changes within the ocean cannot be ignored. By determining the location where abnormal signals occur in the sea area, the changes in ocean dynamics and ecosystems can be better understood. However, the ocean area is vast, and monitoring methods such as cruising or establishing observation stations can only cover a limited ocean area. In addition, the marine environment is changeable, and traditional monitoring methods are difficult to capture the migration paths of substances or signals in the ocean in real time and accurately, which makes it difficult to directly determine the location where abnormal signals occur based on monitoring data.

[0003] Currently, the Geographic Information System (GIS) technology is often used to obtain high-resolution image data through satellite remote sensing, UAV aerial photography, etc. to locate the occurrence of abnormal signals in the sea area, or through the fingerprint analysis method, that is, using three-dimensional fluorescence spectroscopy technology to analyze dissolved organic matter or isotopes in different waters, and realizing pollutant source tracing by calculating the quantitative ratio.

[0004] However, when tracing the source of sea area pollutants based on related technologies, for the analysis of abnormal events in non-surface waters, the image data has certain limitations, which leads to poor accuracy in tracing the source of sea area pollutants; the fingerprint analysis method requires analyzing a large amount of marine sample data, which leads to low efficiency in tracing the source of marine pollutants. Therefore, the solutions of related technologies have problems of poor accuracy and low efficiency in tracing the source of marine pollutants. Summary of the Invention

[0005] The purpose of the present application is to provide a method, a device, a computer device, and a storage medium for tracing substances in a sea area, which can achieve the effect of improving the accuracy and efficiency of tracing the source of marine pollutants.

[0006] The embodiments of the present application are implemented as follows:

[0007] In the first aspect of the embodiments of the present application, a method for tracing substances in a sea area is provided, and the method includes:

[0008] According to the location information of the sea area to be verified, a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field are pre-constructed;

[0009] At least two observation stations are pre - established in the sea area to be verified, and the time series of the monitored substance concentration at each observation station is obtained. The time series of the monitored substance concentration includes the monitored substance concentration changing with time.

[0010] According to the time series of the monitored substance concentration at each observation station, determine the start and end times of the reverse flow field corresponding to each observation station. And according to the reverse flow field, the start and end times, and the position data of each observation station, solve the adjoint equation corresponding to the convection - diffusion equation for each observation station to obtain the spatial distribution state of the adjoint substance concentration and the time series of the adjoint substance concentration at each observation station. The time series of the adjoint substance concentration includes the adjoint substance concentration changing with time.

[0011] According to the time series of the monitored substance concentration at each observation station, the spatial distribution state of the adjoint substance concentration at each observation station, and the time series of the adjoint substance concentration, determine the position information of the target pollution source, the release intensity of the pollutants released by the target pollution source, and the change trend of the release intensity.

[0012] As a possible implementation method, determining the start and end times of the reverse flow field corresponding to each observation station according to the time series of the monitored substance concentration at each observation station includes:

[0013] According to the last moment of the time series of the monitored substance concentration at each observation station, determine the start time of releasing the adjoint substance at each observation station;

[0014] According to the initial moment of the time series of the monitored substance concentration at each observation station, determine the end time of releasing the adjoint substance at each observation station;

[0015] According to the start time and the end time, determine the start and end times of the reverse flow field corresponding to the release of the adjoint substance at each observation station.

[0016] As a possible implementation method, the above - mentioned adjoint equation is as follows:

[0017]

[0018] where c i * is the substance concentration of the adjoint substance released by the observation station i at each grid point in the three - dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three - dimensional orthogonal grid, κ is the mixing coefficient. At the position of the observation station i, let at other positions outside the position of the observation station i i is a positive integer.

[0019] As a possible implementation method, determining the position information of the target pollution source and the release intensity of the pollutants released by the target pollution source according to the time series of the monitored substance concentration at each observation station, the spatial distribution state of the adjoint substance concentration at each observation station, and the time series of the adjoint substance concentration includes:

[0020]

[0021] where c i is the mass concentration of the pollutant released by the target pollution source monitored by the observation station i in the sea area to be verified, and c i * is the mass concentration of the adjoint substance released by the observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q is the release intensity of the pollutant released by the target pollution source at each grid point in the three-dimensional orthogonal grid, (x0, y0) is the location information of the target pollution source, t0 is the initial moment of the monitoring mass concentration time series of the observation station i, T is the last moment of the monitoring mass concentration time series of the observation station i, and at the location of the observation station i, let other locations outside the location of the observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

[0022] As a possible implementation manner, before determining the location information of the target pollution source and the release intensity of the pollutant released by the target pollution source according to the monitoring mass concentration time series of each observation station, the spatial distribution state of the adjoint substance concentration of each observation station, and the adjoint substance concentration time series, it further includes:

[0023]

[0024] where C.V is the release intensity difference coefficient of the corresponding pollutant detected by I observation stations at each grid point in the three-dimensional orthogonal grid, and Q i is the release intensity of the pollutant released by the target pollution source detected by the observation station i at each grid point in the three-dimensional orthogonal grid, is the average release intensity of the pollutant released by the target pollution source detected by I observation stations at each grid point in the three-dimensional orthogonal grid, and I is the set of observation stations.

[0025] As a possible implementation manner, determining the change trend of the release intensity of the pollutant released by the target pollution source according to the monitoring mass concentration time series of each observation station, the spatial distribution state of the adjoint substance concentration of each observation station, and the adjoint substance concentration time series includes:

[0026] Performing time slicing processing on the monitoring mass concentration time series of the observation station to obtain multiple time periods;

[0027] According to the time series of the monitored substance concentrations, the time series of the accompanying substance concentrations, and the location information of the target pollution source at each observation station during each time period, determine the segment release intensity of the pollutants released by the target pollution source during each time period, and determine the change trend of the release intensity of the pollutants released by the target pollution source according to the segment release intensity of the pollutants released by the target pollution source during each time period.

[0028] As a possible implementation, according to the time series of the monitored substance concentrations, the time series of the accompanying substance concentrations, and the location information of the target pollution source at each observation station during each time period, determine the segment release intensity of the pollutants released by the target pollution source during each time period, and determine the change trend of the release intensity of the pollutants released by the target pollution source according to the segment release intensity of the pollutants released by the target pollution source during each time period, including:

[0029]

[0030] where c is the substance concentration of each grid point of the pollutants released by the target pollution source monitored by observation station i in the three-dimensional orthogonal grid, and c * i is the substance concentration of the accompanying substance released by observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q1 to Q n is the segment release intensity of the pollutants released by the target pollution source, (x0, y0) is the location information of the target pollution source, and at the location of observation station i, let other locations outside the location of observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

[0031] In the second aspect of the embodiments of the present application, a device for tracing substances in the sea area is provided. The device includes:

[0032] A construction module for pre-constructing a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field according to the location information of the sea area to be verified;

[0033] An acquisition module for setting up at least one observation station in the pre-constructed three-dimensional orthogonal grid and acquiring the time series of the monitored substance concentrations of each observation station. The time series of the monitored substance concentrations includes multiple monitored substance concentrations that change with time;

[0034] A determination module, configured to determine the time-reversed flow field data of each observation station according to the time series of the monitored substance concentrations of each observation station, and solve the adjoint equation corresponding to the convection-diffusion equation for each observation station according to the time-reversed flow field data and the position data of each observation station, so as to obtain the spatial distribution state of the adjoint substance concentration and the time series of the adjoint substance concentration of each observation station, and the time series of the adjoint substance concentration includes a plurality of adjoint substance concentrations that vary with time;

[0035] The determination module is further configured to determine the position information of the target pollution source, the release intensity of the pollutants released by the target pollution source, and the change trend of the release intensity according to the time series of the monitored substance concentrations of each observation station, the spatial distribution state of the adjoint substance concentration of each observation station, and the time series of the adjoint substance concentration.

[0036] In a third aspect of the embodiments of the present application, a computer device is provided. The computer device includes: a memory and a processor. A computer program that can run on the processor is stored in the memory. When the processor executes the computer program, the steps of the method for tracing substances in the sea area described in the first aspect are implemented.

[0037] In a fourth aspect of the embodiments of the present application, a computer-readable storage medium is provided. A computer program is stored on the computer-readable storage medium. When the computer program is executed by the processor, the steps of the method for tracing substances in the sea area described in the first aspect are implemented.

[0038] The beneficial effects of the embodiments of the present application include:

[0039] A method for tracing substances in a sea area provided by an embodiment of the present application, based on the position information of the sea area to be verified, pre-constructs a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field for the sea area to be verified, so as to obtain a plurality of regular grid areas and water flow laws in the sea area to be verified; sets at least two observation stations in the sea area to be verified, and obtains the time series of the monitored substance concentration monitored by each observation station; according to the time series of the monitored substance concentration monitored by each observation station, reverses the flow field data in the corresponding time period in reverse order and in the reverse direction, and obtains the start and end times of the reversed flow field, and according to the reversed flow field, start and end times and position data of each observation station, solves the adjoint equation corresponding to the convection-diffusion equation for each observation station, and determines the spatial distribution state and time series of the adjoint substance concentration of the adjoint substance released by each observation station based on the solution process of the adjoint equation; according to the time series of the monitored substance concentration monitored by each observation station and the spatial distribution state and time series of the adjoint substance concentration of the adjoint substance released by each observation station, solves the release intensity of the pollutants released by the target pollution source and the position of the target pollution source in the sea area to be verified, and further determines the change trend of the release intensity of the pollutants released by the target pollution source. Among them, by setting a limited number of observation stations in the sea area to be verified to sample the time series of the substance concentration of the pollutants in the sea area to be verified, and by numerically simulating the experiment corresponding to the number of observation stations to solve the adjoint equation in reverse to obtain the adjoint substance concentration, the repeated calculation process in the process of tracing the target pollution source is reduced, thereby improving the tracing efficiency of the target pollution source. In this way, the accuracy and efficiency of tracing marine pollutants can be improved. Description of the Drawings

[0040] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the accompanying drawings required for the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0041] Figure 1 It is a flowchart of the first method for tracing substances in a sea area provided by an embodiment of the present application;

[0042] Figure 2 It is a schematic diagram of a target pollution source releasing pollutants provided by an embodiment of the present application;

[0043] Figure 3 It is a schematic diagram of an observation station releasing an adjoint substance provided by an embodiment of the present application;

[0044] Figure 4 It is a schematic diagram of the position distribution of a target pollution source and an observation station in a sea area provided by an embodiment of the present application;

[0045] Figure 5 This is a flowchart of the second method for tracing substances in the sea area provided by the embodiments of the present application;

[0046] Figure 6 This is a schematic diagram of the release intensity of pollutants released by the target pollution source detected by the first observation station provided by the embodiments of the present application;

[0047] Figure 7 This is a schematic diagram of the release intensity of pollutants released by the target pollution source detected by the second observation station provided by the embodiments of the present application;

[0048] Figure 8 This is a schematic diagram of the release intensity of pollutants released by the target pollution source detected by the third observation station provided by the embodiments of the present application;

[0049] Figure 9 This is a schematic diagram for determining the location of the target pollution source provided by the embodiments of the present application;

[0050] Figure 10 This is a flowchart of the third method for tracing substances in the sea area provided by the embodiments of the present application;

[0051] Figure 11 This is a graph of the change trend of the release intensity of pollutants released by the target pollution source provided by the embodiments of the present application;

[0052] Figure 12 This is a schematic structural diagram of a device for tracing substances in the sea area provided by the embodiments of the present application;

[0053] Figure 13 This is a schematic structural diagram of a computer device provided by the embodiments of the present application. Detailed implementation manners

[0054] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Apparently, the described embodiments are some but not all of the embodiments of the present application. Usually, the components of the embodiments of the present application described and illustrated herein can be arranged and designed in various different configurations.

[0055] Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the present application claimed, but merely represents selected embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present application without creative efforts shall fall within the scope of protection of the present application.

[0056] At present, geographic information system technology is often used for tracing the sources of pollutants in sea areas. High-resolution image data is obtained through satellite remote sensing, UAV aerial photography, etc. to locate the occurrence of abnormal signals in the sea area; or the dissolved organic matter or isotopes in different waters are analyzed by fingerprint analysis methods to calculate the quantitative proportion for tracing the sources of pollutants in the sea area. However, for the analysis of abnormal events in non-surface waters, there are certain limitations in the image data, which will lead to poor accuracy in tracing the sources of pollutants in the sea area. In addition, the fingerprint analysis method requires analyzing a large amount of marine sample data, which results in the problem of low efficiency in tracing the sources of pollutants in the sea area.

[0057] To this end, the embodiments of the present application provide a method for tracing the sources of substances in a sea area. A three-dimensional orthogonal grid and a three-dimensional hydrodynamic flow field are pre-constructed for the sea area to be verified, and at least one observation station is set up in the three-dimensional orthogonal grid; the time series of monitored substances at each observation station are obtained, and the time-reversed flow field data of each observation station are obtained by reverse calculation according to the time series of monitored substances at each monitoring station; according to the time-reversed flow field data of each observation station and the position data of each observation station, the adjoint equation corresponding to the convection-diffusion equation is solved for each observation station to obtain the spatial distribution state of the adjoint substance concentration and the time series of the adjoint substance concentration at each observation station; according to the time series of the monitored substance concentration, the spatial distribution state of the adjoint substance concentration, and the time series of the adjoint substance concentration at each observation station, the position information of the target pollution source, the release intensity of the pollutants released by the target pollution source, and the change trend of the release intensity are determined. In this way, the accuracy and efficiency of tracing the sources of marine pollutants can be improved.

[0058] The method for tracing the sources of substances in the sea area provided by the embodiments of the present application will be explained in detail below.

[0059] Figure 1 The flowchart of a method for tracing the sources of substances in a sea area provided by the present application, which can be applied to computer devices, terminal devices, and third-party servers. Refer to Figure 1 The embodiments of the present application provide a method for tracing the sources of substances in a sea area, and the method includes:

[0060] S101. According to the position information of the sea area to be verified, a three-dimensional orthogonal grid and a three-dimensional hydrodynamic flow field are pre-constructed.

[0061] Optionally, the position information of the sea area to be verified refers to the longitude and latitude information of the sea area to be verified, and the position information of the sea area to be verified includes geographical information such as the water depth, shoreline, and terrain of the sea area to be verified. Among them, the sea area to be verified can be any area in any ocean in the world, and the sea area to be verified can be arbitrarily selected by the user, and the present application does not make specific delineations in this regard.

[0062] Optionally, a three-dimensional orthogonal grid is constructed for the sea area to be verified according to the position information of the sea area to be verified. The irregular sea area to be verified is divided into multiple regular rectangular grid areas through the three-dimensional orthogonal grid, and each grid area has a uniquely determined three-dimensional coordinate. The position of each grid area in the sea area to be verified can be determined through the coordinates. Among them, the grid size of the three-dimensional orthogonal grid can be 100 meters, 300 meters, 500 meters, etc. The smaller the grid size of the three-dimensional orthogonal grid, the more three-dimensional orthogonal grids are obtained by dividing the sea area to be verified, and the more accurate the position of the pollutants in the sea area is indexed. The grid size of the three-dimensional orthogonal grid can be flexibly changed according to the actual situation of the sea area to be verified, and the present application does not make specific limitations on this.

[0063] Optionally, the three-dimensional orthogonal grid of the sea area to be verified can provide a computer grid for the hydrodynamic simulation and source tracing analysis of the sea area to be verified. Through hydrodynamic simulation operations and observational data, the three-dimensional hydrodynamic flow field data of the sea area to be verified are obtained. The three-dimensional hydrodynamic flow field data include the three-dimensional flow velocity, mixing coefficient, etc. of the sea area to be verified. The flow movement laws of each grid area are quantified through the three-dimensional hydrodynamic flow field, and based on the three-dimensional hydrodynamic flow field, the spatio-temporal distributions of parameters such as the velocity, pressure, temperature, and turbulence intensity of the water flow in the horizontal and vertical directions in each grid area can be determined.

[0064] S102. At least two observation stations are pre-established in the sea area to be verified, and the time series of the monitored substance concentration of each observation station is obtained. The time series of the monitored substance concentration includes the monitored substance concentration that changes with time.

[0065] Optionally, the observation station is a monitoring station for tracing the pollutants released by the pollution source in the sea area to be verified. The observation station releases a pollutant's accompanying tracer and then, based on the concentration distribution state of the accompanying tracer in the three-dimensional orthogonal grid, reversely deduces the spatial distribution state of the pollutant.

[0066] Optionally, at least two observation stations are deployed in the sea area to be verified in advance. By monitoring the pollutants in the sea area to be verified through multiple observation stations, the position of the pollution source in the sea area to be verified can be indexed more accurately, and the contingency of monitoring pollutants by a single observation station can be effectively avoided. Among them, the position of the observation station in the sea area to be verified can be set according to the actual situation of the sea area to be verified. For example, if the sea area to be verified is a sea area near a population, the observation station can be deployed within the range of human activities, etc. The present application does not make specific limitations on this.

[0067] It should be noted that when the position of the observation station is set on the pollutant transmission path, the observation station can provide higher-quality source tracing data. Therefore, the position of the observation station is mostly set in the downstream area of the sea area to be verified.

[0068] Optionally, the monitored substance concentration refers to the substance concentration of pollutants passing through the position of the observation station in the sea area to be verified, which is monitored in real time by the observation station. The monitored substance concentration time series refers to the substance concentration of pollutants monitored by the observation station that changes over time. Among them, the observation station monitors the pollutants in the sea area to be verified in real time, and stores and remembers the substance concentration of the monitored pollutants according to time to obtain the monitored substance concentration time series. It should be noted that the monitored substance concentration can be stored separately according to the position of the grid corresponding to the observation station; storing and remembering the substance concentration of pollutants based on time can analyze the trend of pollutant release from the pollution source over time. It should be noted that continuous monitoring equipment is installed in the observation station to automatically sample and measure the substance concentration of pollutants in the sea area to be verified at continuous time or at intervals.

[0069] Optionally, the number of observation stations can be 2, 3, 5, etc. In this embodiment of the application, the number of observation stations is taken as 3 as an example, which does not mean that only three observation stations can be set in the sea area to be verified. The application does not make specific limitations on this. Among them, the target pollution source can be conservative substances such as petroleum pollutants, plastic pollutants, heavy metal pollutants, and industrial chemical pollutants. See Figure 2 , the observation stations are set on the downstream transmission path of the pollutants released by the target pollution source, and the observation stations monitor the polluting substances released by the target pollution source in real time, Figure 2 which is used to characterize that the target pollution source is releasing pollutants in the sea area to be verified.

[0070] It should be noted that the monitored substance concentration time series can characterize how the substance concentration of the monitored substance changes over time at each observation station, can reflect the propagation direction of the monitored substance in the sea area to be verified and the influence of the environmental factors in the sea area to be verified on the propagation of the monitored substance, and can provide a reliable basis for solving the subsequent time-reversed flow field data.

[0071] S103. According to the monitored substance concentration time series of each observation station, determine the start and end times of the reverse flow field corresponding to each observation station, and according to the reverse flow field, the start and end times, and the position data of each observation station, solve the adjoint equation corresponding to the convection-diffusion equation of each observation station to obtain the spatial distribution state of the adjoint substance concentration and the adjoint substance concentration time series of each observation station. The adjoint substance concentration time series includes the adjoint substance concentration that changes over time.

[0072] Optionally, the reverse flow field is a spatio-temporal data set obtained by reversing the time sequence and the direction on the basis of the pre-obtained three-dimensional hydrodynamic flow field data according to the start time of the monitored substance concentration time series. For example, if the monitored substance concentration time series records the substance concentration time series of pollutants from the first day to the fifteenth day, the reverse time sequence is the backward time sequence from the fifteenth day to the first day; the direction reversal means that the flow velocity magnitude remains unchanged while the flow direction is reversed by 180° (for example: the westward flow becomes eastward after reversal).

[0073] Optionally, the position data of the observation station refers to the three-dimensional spatial position of the observation station in the sea area to be verified. According to the reverse flow field corresponding to each observation station, the start and end times of the reverse flow field, and the position data of each observation station, the adjoint equation corresponding to the convection-diffusion equation of each observation station is solved. Among them, the convection-diffusion equation is a partial differential equation that describes the combined action of convection and diffusion of pollutants or adjoint substances in the sea area to be verified. Convection means that pollutants or adjoint substances migrate with seawater, and diffusion means that pollutants or adjoint substances diffuse with seawater movement. It should be noted that the position data of the observation station is also the release position of the adjoint substance.

[0074] Optionally, after multiplying the convection-diffusion equation by the adjoint substance concentration and performing integration in time and space, the relationship equation between the convection-diffusion equation and the adjoint equation is obtained, and then the adjoint equation of the convection-diffusion equation is obtained.

[0075] Optionally, the adjoint substance refers to the substance released by the observation station to trace pollutants in the sea area to be verified. The spatial distribution state of the adjoint substance concentration refers to the substance concentration distribution state of the adjoint substance released by the observation station in each grid area of the three-dimensional orthogonal grid. The time series of the adjoint substance concentration refers to the sequence of the adjoint substance concentration changing with time in each grid area of the sea area to be verified. Among them, the adjoint substance concentration refers to the substance concentration of the adjoint substance released by the observation station to trace the pollution source in each grid area of the three-dimensional orthogonal grid.

[0076] Optionally, the adjoint equation is solved backward based on the time reverse flow field data to obtain the substance concentration distribution of the adjoint substance released by the observation station in the sea area to be verified. See Figure 3 , and the observation stations all release adjoint substances in the sea area to be verified, Figure 3 which is used to characterize the propagation process of the adjoint substance released by the observation station backward to the target pollution source. It should be noted that the observation stations in the sea area to be verified do not release adjoint substances to the target pollution source at the same time. The adjoint substances released by each observation station to the target pollution source do not interfere with each other. Based on the number of observation stations set in the sea area to be verified, the adjoint substances are released to the target pollution source the corresponding number of times to obtain the corresponding number of spatial distribution states of the adjoint substance concentration and the time series of the adjoint substance concentration.

[0077] S104. Determine the position information of the target pollution source, the release intensity of the pollutants released by the target pollution source, and the change trend of the release intensity according to the time series of the monitored substance concentration of each observation station, the spatial distribution state of the adjoint substance concentration of each observation station, and the time series of the adjoint substance concentration.

[0078] Optionally, based on the time series of the monitored substance concentration at each observation station, the spatial distribution state of the adjoint substance concentration at each observation station, and the time series of the adjoint substance concentration, the release intensity of the pollutant released in each grid area of the three-dimensional orthogonal grid in the sea area to be verified can be determined. Then, according to the difference coefficient of the release intensity of the pollutant released obtained from l observation stations in each grid area of the three-dimensional orthogonal grid in the sea area to be verified, the position of the target pollution source in the three-dimensional orthogonal grid can be traced. Among them, the smaller the difference coefficient of the release intensity, the closer the grid area is to the position where the target pollution source is located.

[0079] Optionally, the position information of the target pollution source refers to the coordinates of the target pollution source traced by the observation station in the sea area to be verified; the release intensity refers to the intensity of the pollutant released in each grid area of the three-dimensional orthogonal grid in the sea area to be verified, and its meaning is that the grid point releases pollutants corresponding to Q i releases polluting substances, and the total amount of pollutants monitored by observation station i is consistent with the total amount of pollutants monitored by observation station i when the target pollution source releases polluting substances with the obtained release intensity Q; the release intensity change trend refers to the trend of the release intensity of the pollutant released by the target pollution source changing with time.

[0080] In the embodiment of the present application, based on the position information of the sea area to be verified, a three-dimensional orthogonal grid and a three-dimensional hydrodynamic flow field are pre-constructed for the sea area to be verified to obtain a plurality of regular grid areas and water flow laws in the sea area to be verified; at least two observation stations are set in the sea area to be verified, and the time series of the monitored substance concentration monitored by each observation station is obtained; according to the time series of the monitored substance concentration monitored by each observation station, the flow field data in the corresponding time period are reversed in time and reversed in direction, and the start and end times of the reversed flow field are obtained. Then, according to the reversed flow field, the start and end times, and the position data of each observation station, the adjoint equation corresponding to the convection-diffusion equation of each observation station is solved, and based on the solution process of the adjoint equation, the spatial distribution state of the adjoint substance concentration and the time series of the adjoint substance concentration of the adjoint substance released by each observation station are determined; according to the time series of the monitored substance concentration monitored by each observation station, the spatial distribution state of the adjoint substance concentration of the adjoint substance released by each observation station, and the time series of the adjoint substance concentration, the release intensity of the pollutant released by the target pollution source and the position of the target pollution source in the sea area to be verified are solved, and further the release intensity change trend of the pollutant released by the target pollution source is determined. Among them, by setting a limited number of observation stations in the sea area to be verified to sample the time series of the substance concentration of the pollutant in the sea area to be verified, and by numerically simulating the experiment corresponding to the number of observation stations to solve the adjoint equation in reverse to obtain the adjoint substance concentration, the repeated calculation process in the process of tracing the target pollution source is reduced, thereby improving the tracing efficiency of the target pollution source. In this way, the effect of improving the accuracy and efficiency of tracing marine pollutants can be achieved.

[0081] Figure 4A schematic diagram of the position distribution of the target pollution source and the observation stations in the sea area provided by this application is shown in Figure 4 , the observation stations are distributed in the downstream area of the target pollution source, and the distances between the observation stations and the target pollution source are different. The setting positions of the observation stations in this application are not limited to this.

[0082] In an alternative implementation, as shown in Figure 5 , the operation of "determining the start and end times of the reverse flow field corresponding to each observation station according to the time series of the monitored substance concentrations of each observation station" in step S103 can specifically be:

[0083] S501. Determine the start time of releasing the accompanying substance for each observation station according to the last moment of the time series of the monitored substance concentrations of each observation station.

[0084] Optionally, take the last moment of the time series of the monitored substance concentrations of each observation station as the start time of releasing the accompanying substance for each observation station. For example, if the last moment of the time series of the monitored substance concentrations monitored by observation station P1 is the 15th day, then the start time of releasing the accompanying substance for observation station P1 is the 15th day. Among them, the last moment is used to indicate the last time point when the target pollutant in the sea area to be verified releases pollutants, and the start time is used to indicate the start time point of the reverse release of the accompanying substance by the observation station.

[0085] S502. Determine the end time of releasing the accompanying substance for each observation station according to the initial moment of the time series of the monitored substance concentrations of each observation station.

[0086] Optionally, take the initial moment of the time series of the monitored substance concentrations of each observation station as the end time of releasing the accompanying substance for each observation station. For example, if the initial moment of the time series of the monitored substance concentrations monitored by observation station P1 is the 1st day, then the end time of releasing the accompanying substance for observation station P1 is the 1st day. Among them, the initial moment is used to indicate the first time point when the target pollutant in the sea area to be verified releases pollutants, and the end time is used to indicate the last time point of the reverse release of the accompanying substance by the observation station.

[0087] It should be noted that the time series of releasing the accompanying substance by observation station P1 is the reverse time series of the time series of the pollutants monitored by observation station P1.

[0088] S503. Determine the start and end times of the reverse flow field corresponding to the release of the accompanying substance by each observation station according to the start time and the end time.

[0089] In an alternative implementation, the above-mentioned adjoint equation is as shown in the following formula (1):

[0090]

[0091] Among them, c i * is the substance concentration of the adjoint substance released by observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient. At the position of observation station i, let Positions other than the position of observation station i i is a positive integer.

[0092] Optionally, observation station i refers to any one of the observation stations set in the three-dimensional orthogonal grid of the sea area to be verified.

[0093] Optionally, the substance concentration of the adjoint substance released by each observation station at each grid point in the three-dimensional orthogonal grid of the sea area to be verified can be obtained through the above formula (1).

[0094] In an optional implementation manner, the position information of the target pollution source and the release intensity of the pollutants released by the target pollution source in the above step S104 are obtained through the following formula (2), and formula (2) is as follows:

[0095]

[0096] Among them, c i is the substance concentration of the pollutants released by the target pollution source monitored by observation station i in the sea area to be verified, c i * is the substance concentration of the adjoint substance released by observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q is the release intensity of the pollutants released by the target pollution source at each grid point in the three-dimensional orthogonal grid, (x0, y0) is the position information of the target pollution source, t0 is the initial moment of the monitoring substance concentration time series of observation station i, T is the last moment of the monitoring substance concentration time series of observation station i. At the position of observation station i, let Positions other than the position of observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

[0097] Optionally, based on the monitoring substance concentration time series monitored by each observation station and the adjoint substance concentration time series of the adjoint substances released by each observation station, a numerical value of the release intensity is obtained at each grid point, and the position of the target pollution source in the sea area to be verified is indexed according to the difference coefficient of the release intensity Qi determined by each observation station at the same grid point.

[0098] Optionally, substituting the time series of the monitored substance concentration monitored by the observation station, the horizontal coordinates of each grid point in the three-dimensional orthogonal grid, and the time series of the adjoint substance concentration of the adjoint substance released by the observation station into the above formula (2), the pollutant release intensity of each grid area in the three-dimensional orthogonal grid of the sea area to be verified can be obtained.

[0099] In an alternative embodiment, the release intensity values of the pollutants detected by each observation station in the grid area where the target pollution source is located should theoretically satisfy formula (3), and formula (3) is as follows:

[0100] Q1(x0,y0)=Q2(x0,y0)=…=Q n (x0,y0)(3)

[0101] Wherein, Q i (x0,y0) is used to represent the release intensity of the pollutants released by the target pollution source detected by the observation station i at the grid point (x0,y0) where the target pollution source is located.

[0102] An alternative embodiment Figure 6 is the distribution map of the possible pollutant release intensity in the three-dimensional orthogonal grid of the sea area to be verified obtained by the observation station P1 based on the above formula (2), Figure 7 is the distribution map of the possible pollutant release intensity in the three-dimensional orthogonal grid of the sea area to be verified obtained by the observation station P2 based on the above formula (2), Figure 8 is the distribution map of the possible pollutant release intensity in the three-dimensional orthogonal grid of the sea area to be verified obtained by the observation station P3 based on the above formula (2). See Figure 6 、 Figure 7 and Figure 8 , and the release intensity of the pollutants released by the target pollution source detected by each observation station in the grid area where the target pollution source is located is approximately equal.

[0103] In an alternative embodiment, before determining the location information of the target pollution source and the release intensity of the pollutants released by the target pollution source through the above formula (2), it is also necessary to determine the coefficient of variation of the release intensity of the pollutants released by each observation station detected by the target pollution source through the following formula (4), and formula (4) is as follows:

[0104]

[0105] Wherein, C.V is the coefficient of variation of the release intensity of the corresponding pollutants detected by I observation stations at each grid point in the three-dimensional orthogonal grid. By finding the small value area of C.V, the location of the target pollution source is determined; Q i is the release intensity of the pollutants released by the target pollution source detected by the observation station i at each grid point in the three-dimensional orthogonal grid, and its meaning is that the release intensity corresponding to this grid point is Q iRelease pollutants, and the total amount of pollutants monitored by the observation station i is consistent with the total amount of pollutants monitored by the observation station i when the target pollution source releases pollutants with the required release intensity Q. is the average release intensity of the target pollution source releasing pollutants detected at each grid point of the I observation stations in the three-dimensional orthogonal grid. After determining the location of the target pollution source, the corresponding is the release intensity of the required pollutants; I is the set of observation stations.

[0106] Optionally, the release intensity difference coefficient is used to describe the difference in the release intensity of the target pollution source releasing pollutants detected by different observation stations in the same grid area. The smaller the value of the release intensity difference coefficient, the smaller the difference in the release intensity detected between the observation stations. By limiting the range of the release intensity variation coefficient, grid areas with relatively small differences in the inverted release intensity are screened out, that is, the location of the grid area where the target pollution source is located is determined, and the average value of the release intensities detected by each observation in the grid area where the target pollution source is located is used as the release intensity of the target pollution source releasing pollutants.

[0107] Optionally, see Figure 9 , within the area where the release intensity coefficient C.V is set to be less than 15%, the target pollution source is located at Figure 9 the position shown. The release intensity coefficient C.V can also be set to 5%, 10%, 15%, etc. This application does not make specific limitations on this.

[0108] In a possible implementation manner, see Figure 10 , the operation of step S104 can specifically be:

[0109] S1001. Perform time slicing processing on the time series of the monitored substance concentration of the observation station to obtain multiple time periods.

[0110] Optionally, perform time slicing on the time series of the monitored substance concentration monitored by the observation station, for example, slice it into n time periods such as T0 - T1, T1 - T2,..., T n-1 -T n and so on.

[0111] S1002. Determine the segment release intensity of the target pollution source releasing pollutants in each time period according to the time series of the monitored substance concentration, the time series of the accompanying substance concentration, and the location information of the target pollution source at each observation station in each time period, and determine the change trend of the release intensity of the target pollution source releasing pollutants according to the segment release intensity of the target pollution source releasing pollutants in each time period.

[0112] Optionally, the segment release intensity refers to the release intensity of the target pollution source releasing pollutants in each time period after slicing. According to the continuous segment release intensity, the change trend of the release intensity of the target pollution source releasing pollutants can be determined.

[0113] In an alternative embodiment, the variation trend of the release intensity of the pollutants released by the target pollution source in the above step S1002 can be achieved by formula (5), and formula (5) is as follows:

[0114]

[0115] where c i is the mass concentration of the pollutants released by the target pollution source monitored by the observation station i in the sea area to be verified, c i * is the mass concentration of the accompanying substance released by the observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q1 to Q n are the segment release intensities of the pollutants released by the target pollution source in different time periods, (x0, y0) is the location information of the target pollution source, and at the location of the observation station i, let other positions outside the location of the observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

[0116] In an alternative embodiment, if the target pollution source releases pollutants with a release intensity of Q1 in the time period T0 - T1, the target pollution source releases pollutants with a release intensity of Q2 in the time period T1 - T2, and the target pollution source releases pollutants with a release intensity of Q n-1 -T n with a release intensity of Q n where the release intensity Q1 of the first time period can be obtained from the above formula (5), and the relational equation satisfies the following formula (6) in the time period T0 - T2:

[0117]

[0118] Optionally, the release intensity Q1 obtained from the above formula (5) is substituted into formula (6) to obtain the release intensity Q2 of the second time period T1 - T2. It should be noted that c in the equation of formula (6) i * needs to be solved again through numerical simulation experiments. Taking the moment T2 as the first moment of the flow field, the simulation duration is from T2 to T0. And so on, the segment release intensity of each time period is solved to obtain the release intensities of multiple segments.

[0119] Figure 11 is a graph of the variation trend of the release intensity of the pollutants released by the target pollution source provided by this application. See Figure 11, the solid line represents the change in the release intensity of pollutants released from a preset target pollution source, and the dashed line represents the change in the release source intensity of pollutants released from the target pollution source detected by the observation station through an inversion experiment.

[0120] The following describes the devices, equipment, computer-readable storage media, etc. for implementing the method for tracing substances in the sea area provided by the present application. For the specific implementation process and technical effects, please refer to the above, and will not be repeated below.

[0121] Figure 12 It is a schematic structural diagram of a device for tracing substances in a sea area provided by an embodiment of the present application. Refer to Figure 12 , the device includes:

[0122] A construction module 1201, configured to pre-construct a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field according to the position information of the sea area to be verified;

[0123] An acquisition module 1202, configured to pre-establish at least two observation stations in the sea area to be verified and acquire the monitoring substance concentration time series of each observation station. The monitoring substance concentration time series includes a plurality of monitoring substance concentrations that change with time;

[0124] A determination module 1203, configured to determine the start and end times of the reverse flow field corresponding to each observation station according to the monitoring substance concentration time series of each observation station, and solve the adjoint equation corresponding to the convection-diffusion equation for each observation station according to the reverse flow field, the start and end times, and the position data of each observation station, so as to obtain the spatial distribution state of the adjoint substance concentration and the adjoint substance concentration time series of each observation station. The adjoint substance concentration time series includes the adjoint substance concentration that changes with time;

[0125] The determination module 1203 is further configured to determine the position information of the target pollution source, the release intensity of pollutants released from the target pollution source, and the release intensity change trend according to the monitoring substance concentration time series of each observation station, the spatial distribution state of the adjoint substance concentration of each observation station, and the adjoint substance concentration time series.

[0126] In an optional implementation manner, the above determination module 1203 may specifically be configured to:

[0127] Determine the start time of releasing the adjoint substance for each observation station according to the last moment of the monitoring substance concentration time series of each observation station;

[0128] Determine the end time of releasing the adjoint substance for each observation station according to the initial moment of the monitoring substance concentration time series of each observation station;

[0129] Determine the start and end times of the reverse flow field corresponding to the release of the adjoint substance for each observation station according to the start time and the end time.

[0130] In an alternative embodiment, the determining module 1203 is specifically further configured to:

[0131] Perform time slicing processing on the time series of the monitored substance concentration of the observation station to obtain multiple time periods;

[0132] According to the time series of the monitored substance concentration, the time series of the accompanying substance concentration, and the location information of the target pollution source in each time period of each observation station, determine the segment release intensity of the pollutants released by the target pollution source in each time period, and according to the segment release intensity of the pollutants released by the target pollution source in each time period, determine the release intensity change trend of the pollutants released by the target pollution source.

[0133] The above device is used to execute the method provided in the foregoing embodiment, and its implementation principle and technical effect are similar, and will not be elaborated here.

[0134] The above modules may be one or more integrated circuits configured to implement the above method. For example: one or more application specific integrated circuits (ASICs), or, one or more microprocessors, or, one or more field programmable gate arrays (FPGAs), etc. Again, when the above certain module is implemented in the form of a processing element scheduling program code, the processing element may be a general-purpose processor, such as a central processing unit (CPU) or other processors that can call program code. Again, these modules may be integrated together and implemented in the form of a system-on-a-chip (SOC).

[0135] Figure 13 It is a schematic structural diagram of a computer device provided by an embodiment of the present application. Refer to Figure 13 , the computer device includes: a memory 1301 and a processor 1302. A computer program that can run on the processor 1302 is stored in the memory 1301. When the processor 1302 executes the computer program, the steps in any of the above method embodiments are implemented.

[0136] An embodiment of the present application further provides a computer-readable storage medium. The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments can be implemented.

[0137] Optionally, the present application further provides a program product, such as a computer-readable storage medium, including a program, which is used to execute the method embodiment of tracing the source of substances in any of the above sea areas when executed by a processor.

[0138] In several embodiments provided by the present invention, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections between each other can be through some interfaces. The indirect couplings or communication connections of devices or units can be in electrical, mechanical or other forms.

[0139] The units described as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place, or they can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0140] In addition, the functional units in each embodiment of the present invention can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above-mentioned integrated units can be implemented in the form of hardware, or in the form of a combination of hardware and software functional units.

[0141] The above-mentioned integrated units implemented in the form of software functional units can be stored in a computer-readable storage medium. The above-mentioned software functional units stored in a storage medium include several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor (English: processor) to execute some steps of the methods in each embodiment of the present invention. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (English: Read-Only Memory, abbreviated as: ROM), random access memories (English: Random Access Memory, abbreviated as: RAM), magnetic disks or optical discs that can store program codes.

[0142] The above is only the specific implementation manner of this application, but the protection scope of this application is not limited thereto. Any person skilled in the art can easily think of changes or substitutions within the technical scope disclosed in this application, and all should be covered by the protection scope of this application. Therefore, the protection scope of this application should be subject to the protection scope of the claims.

[0143] The above are only the preferred embodiments of the present application and are not intended to limit the present application. For those skilled in the art, various modifications and variations can be made to the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included within the protection scope of the present application.

Claims

1. A method for tracing substances in the sea area, characterized in that, The method includes: Pre - construct a three - dimensional orthogonal grid and a hydrodynamic three - dimensional flow field according to the position information of the sea area to be verified; Pre - set at least two observation stations in the sea area to be verified, and obtain the time series of the monitored substance concentration at each of the observation stations, where the time series of the monitored substance concentration includes the monitored substance concentration varying with time; According to the time series of the monitored substance concentration at each of the observation stations, determine the start and end times of the reverse flow field corresponding to each of the observation stations, and solve the adjoint equation corresponding to the convection - diffusion equation for each of the observation stations based on the reverse flow field, the start and end times, and the position data of each of the observation stations, so as to obtain the spatial distribution state of the adjoint substance concentration and the time series of the adjoint substance concentration at each of the observation stations, where the time series of the adjoint substance concentration includes the adjoint substance concentration varying with time; According to the time series of the monitored substance concentration at each of the observation stations, the spatial distribution state of the adjoint substance concentration at each of the observation stations, and the time series of the adjoint substance concentration, determine the position information of the target pollution source, the release intensity of the pollutants released by the target pollution source, and the change trend of the release intensity.

2. The method for tracing substances in the sea area according to claim 1, characterized in that The step of determining the start and end times of the reverse flow field corresponding to each of the observation stations according to the time series of the monitored substance concentration at each of the observation stations includes: Determine the start time of releasing the adjoint substance at each of the observation stations according to the last moment of the time series of the monitored substance concentration at each of the observation stations; Determine the end time of releasing the adjoint substance at each of the observation stations according to the initial moment of the time series of the monitored substance concentration at each of the observation stations; Determine the start and end times of the reverse flow field corresponding to the release of the adjoint substance at each of the observation stations according to the start time and the end time.

3. The method for tracing substances in the sea area according to claim 1, characterized in that, The adjoint equation is as follows: Among them, c i * is the substance concentration at each grid point of the adjoint substance released by the observation station i in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, and at the position of the observation station i, let at other positions outside the position of the observation station i i is a positive integer.

4. The method for tracing substances in the sea area according to claim 1, characterized in that The step of determining the position information of the target pollution source and the release intensity of the pollutants released by the target pollution source according to the time series of the monitored substance concentration at each of the observation stations, the spatial distribution state of the adjoint substance concentration at each of the observation stations, and the time series of the adjoint substance concentration includes: where c i is the substance concentration of the pollutant released by the target pollution source monitored by the observation station i in the sea area to be verified, c i * is the substance concentration of the accompanying substance released by the observation station i at each grid point in the three-dimensional orthogonal grid, t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q is the release intensity of the pollutant released by the target pollution source at each grid point in the three-dimensional orthogonal grid, (x0, y0) is the location information of the target pollution source, t0 is the initial moment of the monitoring substance concentration time series of the observation station i, T is the last moment of the monitoring substance concentration time series of the observation station i, at the location of the observation station i, let other locations outside the location of the observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

5. The method for tracing substances in the sea area according to claim 4, characterized in that, Before the step of determining the position information of the target pollution source and the release intensity of the pollutants released by the target pollution source according to the time series of the monitored substance concentration at each of the observation stations, the spatial distribution state of the adjoint substance concentration at each of the observation stations, and the time series of the adjoint substance concentration, it further includes: Among them, C.V is the difference coefficient of the release intensity of the corresponding pollutants detected at each grid point in the three-dimensional orthogonal grid by I observation stations, and Q i is the release intensity of the pollutants detected at each grid point in the three-dimensional orthogonal grid by the observation station i, is the average release intensity detected at each grid point in the three-dimensional orthogonal grid by I observation stations, and I is the set of observation stations.

6. The method for tracing substances in the sea area according to claim 1, characterized in that The step of determining the change trend of the release intensity of the pollutants released by the target pollution source according to the time series of the monitored substance concentration at each of the observation stations, the spatial distribution state of the adjoint substance concentration at each of the observation stations, and the time series of the adjoint substance concentration includes: Perform time slicing processing on the time series of the monitored substance concentration at the observation stations to obtain multiple time periods; According to the time series of the monitored substance concentration, the time series of the adjoint substance concentration at each observation station within each time period, and the position information of the target pollution source, determine the segment release intensity of the pollutants released by the target pollution source within each time period, and determine the change trend of the release intensity of the pollutants released by the target pollution source according to the segment release intensity of the pollutants released by the target pollution source within each time period.

7. The method for tracing substances in the sea area according to claim 6, wherein Determining the segment release intensity of pollutants released by the target pollution source in each time period according to the monitoring substance concentration time series, the accompanying substance concentration time series of each observation station in each time period, and the location information of the target pollution source, and determining the release intensity change trend of pollutants released by the target pollution source according to the segment release intensity of pollutants released by the target pollution source in each time period, including: Among them, c i is the mass concentration of the pollutant released by the target pollution source monitored by the observation station i in the sea area to be verified. c i * is the mass concentration of the accompanying substance released by the observation station i at each grid point in the three-dimensional orthogonal grid. t is the time, u is the seawater flow velocity at each grid point in the three-dimensional orthogonal grid, κ is the mixing coefficient, Q1 to Q n is the sectional release intensity of the pollutant released by the target pollution source. (x0, y0) is the location information of the target pollution source. At the location of the observation station i, let other positions outside the location of the observation station i i is a positive integer, and (x, y) are the coordinates of each grid point in the three-dimensional orthogonal grid.

8. An apparatus for tracing substances in a sea area, characterized in that, The device includes: A construction module for pre-constructing a three-dimensional orthogonal grid and a hydrodynamic three-dimensional flow field according to the location information of the sea area to be verified; An acquisition module for pre-establishing at least two observation stations in the sea area to be verified and acquiring the monitoring substance concentration time series of each of the observation stations, where the monitoring substance concentration time series includes a plurality of monitoring substance concentrations that change over time; A determination module for determining the start and end times of the reverse flow field corresponding to each of the observation stations according to the monitoring substance concentration time series of each of the observation stations, and solving the adjoint equation corresponding to the convection-diffusion equation for each of the observation stations according to the reverse flow field, the start and end times, and the location data of each of the observation stations, so as to obtain the spatial distribution state of the adjoint substance concentration and the adjoint substance concentration time series of each of the observation stations, where the adjoint substance concentration time series includes the adjoint substance concentration that changes over time; The determination module is further configured to determine the location information of the target pollution source, the release intensity of pollutants released by the target pollution source, and the release intensity change trend according to the monitoring substance concentration time series of each of the observation stations, the spatial distribution state of the adjoint substance concentration of each of the observation stations, and the adjoint substance concentration time series; 9. A computer device, characterized in that, Including: A memory and a processor, where the memory stores a computer program that can run on the processor, and when the processor executes the computer program, the steps of the method described in any one of claims 1 to 7 above are implemented.

10. A computer-readable storage medium, characterized in that, A computer program is stored on the computer-readable storage medium, and when the computer program is executed by the processor, the steps of the method described in any one of claims 1 to 7 are implemented.

Images