Gas diffusion translation calculation method and device, computer equipment and storage medium

Abstract

The embodiment of the invention discloses a gas diffusion translation calculation method and device, computer equipment and a storage medium. The method comprises the following steps: setting gas diffusion related parameters of each standard puff; initializing an observation space, an observation duration and an observation time point, and calculating a downwind direction and a crosswind direction; calculating a concentration matrix sequence of the standard puff of each ship at each time point within the observation duration to obtain a concentration matrix corresponding to each ship; and calculating a puff group concentration matrix of a plurality of ships at a specified observation point in the observation space according to the concentration matrix by adopting a standard puff multi-translation multiplexing mode. By implementing the method provided by the embodiment of the invention, the algorithm complexity during multi-ship continuous emission calculation can be greatly simplified, the repeated calculation amount is reduced, the memory space occupation is reduced, and the calculation time consumption is reduced.

Description

Technical Field

[0001] The present invention relates to a gas monitoring method, and more specifically to a gas diffusion translation calculation method, device, computer equipment and storage medium. Background Art

[0002] According to the puff emission diffusion theory, a puff emitted by a ship at time t0 will diffuse in the air over time t. i The gradual spread of As shown. Among them, represents the emission intensity at the spatial point (x0, y0, z0) at time t0, u is the wind speed, are the diffusion coefficients of the puff in the x, y, and z directions, respectively, and H e is the height of the emission source.

[0003] For a ship in continuous motion, there are countless smoke groups with extremely short time intervals along its trajectory, and their diffusion overlaps with each other. Therefore, the ship k is n ) The emission concentration of the smoke group continuously emitted during the time period C kt That is, for all observation points t i concentration The sum of Similarly, the continuous emission diffusion concentration C of multiple ships in the same navigation area t is the concentration C of the smoke group emissions from multiple emission sources kt The sum of C t =∑ k C kt .

[0004] It can be seen from this that when calculating continuous emissions from multiple ships, it is necessary to repeatedly calculate and superimpose the instantaneous emissions of a single ship or a single point, which makes the entire algorithm highly complex, requires a large amount of repeated calculations, occupies a lot of memory space, and takes a long time to calculate, which poses a major technical obstacle to the research and real-time monitoring of the diffusion of pollution emissions from multiple ships in a region.

[0005] Therefore, it is necessary to design a new method to simplify the algorithm complexity when calculating continuous emissions of multiple ships, reduce the amount of repeated calculations, reduce memory space usage, and reduce calculation time. Summary of the invention

[0006] The purpose of the present invention is to overcome the defects of the prior art and provide a gas diffusion translation calculation method, device, computer equipment and storage medium.

[0007] To achieve the above object, the present invention adopts the following technical solution: a gas diffusion translation calculation method, comprising:

[0008] Set gas diffusion related parameters for each standard puff;

[0009] Initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction;

[0010] Calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation period to obtain the concentration matrix corresponding to each ship;

[0011] The smoke puff group concentration matrix of several ships at the designated observation point in the observation space is calculated according to the concentration matrix by adopting the method of multiple translation and reuse of standard smoke puffs.

[0012] A further technical solution is: the gas diffusion related parameters include wind direction, wind speed, emission amount, emission point and diffusion factor.

[0013] A further technical solution is: the calculation of the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship includes:

[0014] Based on the concentration resolution, the standard smoke puffs of each ship at each time point within the observation time are calculated and filled into the initialized concentration matrix in sequence with longitude, latitude, and concentration vector to obtain the concentration matrix corresponding to each ship.

[0015] Its further technical solution is: the concentration matrix is C i represents the concentration of m spatial points observed at observation time i, where is the spatial coordinate location of observation point j corresponding to time point i, is the coordinate point The corresponding smoke puff diffusion concentration, j∈(1,2,…,m).

[0016] A further technical solution is: the method of using multiple translation and reuse of standard smoke puffs is used to calculate the smoke puff group concentration matrix of several ships at a specified observation point in the observation space according to the concentration matrix, including:

[0017] Determine the concentration matrix corresponding to the ship at a certain observation time point, and calculate the diffusion concentration matrix by multiple translation and reuse of standard smoke puffs;

[0018] All diffusion concentration matrices are superimposed to obtain the smoke group concentration matrix of several ships at a specified observation point in the observation space.

[0019] Its further technical solution is: the diffusion concentration matrix is in, is the coordinate of the instantaneous emission point of the standard puff, are the coordinates of the instantaneous emission point of the target puff, and (Δp1, Δp2) are the translation of the target puff relative to the reference puff.

[0020] A further technical solution is: all diffusion concentration matrices are superimposed to obtain the smoke group concentration matrix of several ships at a specified observation point in the observation space, including:

[0021] All diffusion concentration matrices are connected in series to form an intermediate matrix;

[0022] Set the latitude and longitude resolution;

[0023] Unifying all longitude and latitude coordinates in the intermediate matrix to the longitude and latitude resolution;

[0024] The rows with the same longitude and latitude values ​​in the unified intermediate matrix are merged into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0025] The present invention also provides a gas diffusion translation calculation device, comprising:

[0026] A setting unit, used to set gas diffusion related parameters of each standard puff;

[0027] Initialization unit, used to initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction;

[0028] The first calculation unit is used to calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship;

[0029] The second calculation unit is used to calculate the smoke group concentration matrix of several ships at the designated observation point in the observation space according to the concentration matrix by using the standard smoke puff multiple translation multiplexing method.

[0030] The present invention further provides a computer device, which includes a memory and a processor. The memory stores a computer program, and the processor implements the above method when executing the computer program.

[0031] The present invention also provides a storage medium, wherein the storage medium stores a computer program, and the computer program implements the above method when executed by a processor.

[0032] Compared with the prior art, the beneficial effect of the present invention is as follows: the present invention sets the diffusion parameters and observation conditions of the smoke puff, then calculates the emission concentration matrix of the ship within the specified observation time, and finally uses the multiple translation multiplexing method to calculate the smoke puff concentration matrix of multiple ships in the observation space, thereby simplifying the algorithm complexity when performing continuous emission calculations for multiple ships, reducing the amount of repeated calculations, reducing the memory space occupied, and reducing the calculation time.

[0033] The present invention is further described below in conjunction with the accompanying drawings and specific embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying any creative work.

[0035] Figure 1 A schematic diagram of an application scenario of the gas diffusion translation calculation method provided in an embodiment of the present invention;

[0036] Figure 2 A schematic diagram of a flow chart of a gas diffusion translation calculation method provided in an embodiment of the present invention;

[0037] Figure 3 A schematic diagram of a sub-process of a gas diffusion translation calculation method provided in an embodiment of the present invention;

[0038] Figure 4 A schematic diagram of a sub-process of a gas diffusion translation calculation method provided in an embodiment of the present invention;

[0039] Figure 5 A computer schematic diagram of a gas diffusion translation calculation method provided in an embodiment of the present invention;

[0040] Figure 6 A schematic diagram of the average time consumption ratio between the traditional algorithm and the gas diffusion translation calculation method provided in an embodiment of the present invention;

[0041] Figure 7 A schematic block diagram of a gas diffusion translation calculation device provided by an embodiment of the present invention;

[0042] Figure 8 A schematic block diagram of a computer device provided for an embodiment of the present invention;

[0043] Fig. 9 A schematic diagram of a wind direction projection diagram provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0044] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0045] It should be understood that when used in this specification and the appended claims, the terms "include" and "comprises" indicate the presence of described features, integers, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or combinations thereof.

[0046] It should also be understood that the terms used in this specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include plural forms.

[0047] It should be further understood that the term "and / or" used in the present description and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes these combinations.

[0048] See also Figure 1 and Figure 2 , Figure 1 A schematic diagram of an application scenario of the gas diffusion translation calculation method provided in an embodiment of the present invention. Figure 2 A schematic flow chart of a gas diffusion translation calculation method provided in an embodiment of the present invention. The gas diffusion translation calculation method is applied in a server. The server interacts with the terminal for data, and separates the calculation of the standard puff, and its time complexity is O(T). For the entire diffusion process of a group of smoke puffs continuously emitted by k emission sources, there are a total of k*T smoke puffs, and the translation calculation is performed on each puff, so the complexity of the entire algorithm is O(T+KT), which can be simplified to O(KT). Using the traditional calculation formula requires constantly calling the concentration matrix corresponding to the standard puff to calculate the diffusion concentration of the target area, so the complexity of the entire algorithm is O(KT 2 ). This has a great advantage in terms of time complexity and can achieve rapid calculation of the diffusion process of smoke clusters.

[0049] Figure 2 FIG. 1 is a flow chart of a gas diffusion translation calculation method provided by an embodiment of the present invention. Figure 2 As shown, the method includes the following steps S110 to S140.

[0050] S110. Setting gas diffusion related parameters of each standard puff.

[0051] In this embodiment, the gas diffusion related parameters include wind direction, wind speed, emission amount, emission point and diffusion factor.

[0052] Assume that the emission and diffusion process of each puff is independent. The diffusion process of the puff after emission is only related to the initial emission (emission point, emission amount), wind (wind force, wind direction), atmospheric diffusion parameters (diffusion coefficients of x, y, z axes) and time, and has nothing to do with the subsequent movement state of the emission source. For similar exhaust gas emissions, their meteorological conditions and atmospheric diffusion parameters are relatively stable within a relatively small observation range. It can be considered that the diffusion of the puff is only related to the initial emission source point and time parameters.

[0053] The standard puff can be any puff selected as a diffusion reference.

[0054] S120, initializing the observation space, observation duration, and observation time point, and calculating the downwind direction and the crosswind direction.

[0055] In this embodiment, the spatial range and time range of the observation area and the time interval of the observation time points need to be initialized. Calculating the downwind direction and the crosswind direction helps to determine the propagation path of the pollutants under different wind direction conditions.

[0056] Specifically, downwind direction usually refers to the direction of pollutant propagation along the wind direction at a given wind speed and wind direction; crosswind direction refers to the diffusion of pollutants in the direction perpendicular to the main wind direction. Smoke puffs always diffuse along the downwind direction. Assuming the wind direction θ, the downwind direction is β = θ-180°. The projection calculation of the downwind direction on the x and y axes is as follows Fig. 9 shown.

[0057] According to the angle relationship, when 0°≤θ-180°<180°, u x =-u*sinθ, when 180°≤θ-180°<360°, u x =u*sinθ; when 270°≤θ-180° or θ-180°<90°, u x =-u*cosθ; when 90°≤θ-180°<270°, u x =u*cosθ.

[0058] S130, calculating the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time, so as to obtain the concentration matrix corresponding to each ship.

[0059] In this embodiment, the concentration matrix corresponding to each ship refers to the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time.

[0060] Specifically, based on the concentration resolution, the standard smoke puffs of each ship at each time point within the observation period are calculated and filled into the initialized concentration matrix in sequence with longitude, latitude, and concentration vector to obtain the concentration matrix corresponding to each ship.

[0061] Wherein, the concentration matrix is C i represents the concentration of m spatial points observed at observation time i, where is the spatial coordinate location of observation point j corresponding to time point i, is the coordinate point The corresponding smoke puff diffusion concentration, j∈(1,2,…,m). Based on the concentration resolution, the smoke puffs in the target area are represented in the form of a concentration matrix consisting of m rows and 3 columns of vectors (longitude, latitude, concentration) to form a concentration matrix.

[0062] S130, using a method of multiple translation and reuse of standard smoke puffs, and calculating the smoke puff group concentration matrix of several ships at designated observation points in the observation space according to the concentration matrix.

[0063] In this embodiment, the smoke group concentration matrix of a plurality of ships at a designated observation point in the observation space refers to the entire concentration matrix of the plurality of ships at the designated observation point in the observation space.

[0064] See also Figure 5 , the two instantaneous emission puffs during the operation of the emission source, under the same atmospheric diffusion parameters and wind influence, follow the same diffusion function, the difference is the location of the puff and the duration of the diffusion process, the instantaneous emission of any point at a specific time is equivalent to the result of coordinate translation of the standard puff at the corresponding diffusion time. The diffusion concentration matrix can be calculated by translation calculation.

[0065] In one embodiment, see Figure 3 , the above-mentioned step S130 may include steps S131 to S132.

[0066] S131. Determine the concentration matrix corresponding to the ship at a certain observation time point, and calculate the diffusion concentration matrix by using a method of multiple translation and reuse of standard smoke puffs.

[0067] In this embodiment, the diffusion concentration matrix refers to the concentration matrix after the smoke puff diffuses.

[0068] Specifically, the diffusion concentration matrix is in, is the coordinate of the instantaneous emission point of the standard puff, is the instantaneous emission point coordinate of the target puff, (Δp1, Δp2) is the translation of the target puff relative to the reference puff; the puff translation calculation can be realized directly through matrix linear calculation.

[0069] S132. All diffusion concentration matrices are superimposed to obtain smoke group concentration matrices of several ships at designated observation points in the observation space.

[0070] In one embodiment, see Figure 4 , the above-mentioned step S132 may include steps S1321 to S1324.

[0071] S1321. Connect all diffusion concentration matrices in series to form an intermediate matrix.

[0072] In this embodiment, the intermediate matrix refers to a matrix formed by connecting and expanding n puff concentration matrices in series.

[0073] Although the emission and diffusion of smoke puffs are independent of each other, they will cross and overlap in the same observation area, and concentration superposition processing is required. That is, the pollution concentration corresponding to a point in the air is the sum of the concentrations of each smoke puff at that point. The calculated concentration matrices have different numbers of rows but the same number of columns. First, simply concatenate n puff concentration matrices together, where: is the diffusion concentration matrix.

[0074] S1322. Set the longitude and latitude resolution.

[0075] In this embodiment, as the number of smoke puffs increases, the matrix size will also increase dramatically, which will create a huge burden on storage. At the same time, although too fine a longitude and latitude resolution will promote display accuracy, it will have no effect on concentration calculation. Therefore, a certain longitude and latitude resolution is set, such as e=1e-5.

[0076] S1323, unifying all longitude and latitude coordinates in the intermediate matrix to longitude and latitude resolution;

[0077] S1324. The rows with the same longitude and latitude values ​​in the unified intermediate matrix are merged into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0078] In this embodiment, all longitude and latitude coordinates of the concentration matrix are unified to the corresponding resolution according to the rounding principle, and then the rows of the concentration matrix with the same longitude and latitude values ​​are merged into one row, and the corresponding concentration value is the sum of the concentration values ​​of the merged rows, thereby obtaining a greatly simplified diffusion concentration matrix of all ship emission smoke groups in the observation area at the set resolution.

[0079] In specific applications, set wind direction u, wind speed udir, emission volume Q, emission point [x0, y0, H0], diffusion factor σx , σ y , σ z ; Observation duration T, observation time point i; Initialize a 3D zero matrix C (T * N * 3); Calculate the downwind direction d wind and the crosswind direction c wind , Set a counter t = 1; When t is less than T, calculate and t = t + 1, and re - judge whether t is less than T; When t is not less than T, the concentration matrix sequence of the whole process of the standard emission plume diffusion can be obtained.

[0080] After that, repeatedly refer to the standard plume concentration matrix sequence and use the fast translation method to calculate the diffusion concentration of the plumes emitted by different ships at different time points. Calculate the diffusion duration t (observation time point minus emission time point, t < T) after a certain plume emission of ship k at observation time point i, then the corresponding concentration matrix sequence C of the plume at observation time point T i subset C t after position translation calculation. Round all longitude and latitude coordinates of the concentration matrix to the standard resolution according to the set observation accuracy, and then merge the rows with the same longitude and latitude values of the concentration matrix into 1 row, and the concentration value is the sum of the concentration values of the merged rows, so as to simplify the concentration matrix without affecting the accuracy.

[0081] Loop through N ships. When n < N, judge whether the plume s of ship n is less than the number of concentration matrixes of its corresponding standard concentration subset C , that is, the number of plumes. When s is less than C the number of, intercept the subset, perform translation calculation and accuracy processing according to the above method, and superimpose the diffusion concentrations. Set s equal to s + 1, and return to judge whether s is less than C the number of. When s is not less than C the number of, all the concentration calculations of this ship are completed, set n = n + 1, calculate the concentration matrix of the next ship, and execute the judgment that n is less than N. When n is not less than N, it indicates that all the concentration translations and superposition calculations of all ships and all plumes are completed, then output the superimposed plume result.

[0082] According to the result of comparing the number of emission sources and the observation window radius angle, the method of this embodiment is significantly superior to the traditional algorithm in terms of time consumption. In the multi - source emission diffusion algorithm, the number of emission sources directly affects the time efficiency of the algorithm, so it is an important index to measure the algorithm efficiency. By comparing the average time - consumption ratio of the traditional algorithm and the new algorithm, the improvement degree of the new algorithm in time efficiency is clearly shown.

[0083] To explore the relationship between the efficiency improvement and the number of emission sources, a trend graph of the time - consumption ratio (R) changing with the number of emission sources (N) is constructed, as Figure 6 As shown. The curve shows the change in the average time consumption ratio in the two stages. In the first stage, when the number of emission sources increased from 1 to 5, the average time consumption ratio of the traditional algorithm and the new algorithm decreased from about 200 to about 145. Due to the small number of emission sources, the advantage of the new algorithm in dealing with multi-source emissions is not obvious. In the second stage, when the number of emission sources increased from 5 to 20, the average time consumption ratio increased significantly, showing the advantage of the algorithm in large-scale emission processing. Therefore, judging from the trend of the curve, the time efficiency of the new algorithm is linearly improved compared with the traditional algorithm when dealing with more emission sources.

[0084] Sulfide, nitrogen oxide, carbon oxide, and suspended particulate matter emitted by ships and vehicles are the main components of air pollution and health impacts. In particular, carbon emissions are directly related to the carbon content of fuel and combustion power. The method of this embodiment can be used in the fields of science and technology, equipment, policy research, and supervision of pollution emission diffusion, and has broad application prospects.

[0085] The above-mentioned gas diffusion translation calculation method sets the diffusion parameters and observation conditions of the smoke puff, then calculates the emission concentration matrix of the ship within the specified observation time, and finally uses the multiple translation multiplexing method to calculate the smoke puff concentration matrix of multiple ships in the observation space, thereby simplifying the algorithm complexity when performing continuous emission calculations for multiple ships, reducing the amount of repeated calculations, reducing memory space occupancy, and reducing calculation time.

[0086] Figure 7 is a schematic block diagram of a gas diffusion translation calculation device 300 provided in an embodiment of the present invention. Figure 7 As shown, corresponding to the above gas diffusion translation calculation method, the present invention also provides a gas diffusion translation calculation device 300. The gas diffusion translation calculation device 300 includes a unit for executing the above gas diffusion translation calculation method, and the device can be configured in a server. Figure 7 The gas diffusion translation calculation device 300 includes a setting unit 301, an initialization unit 302, a first calculation unit 303 and a second calculation unit 304.

[0087] The setting unit 301 is used to set the gas diffusion related parameters of each standard smoke puff; the initialization unit 302 is used to initialize the observation space, observation time, observation time point, and calculate the downwind direction and the crosswind direction; the first calculation unit 303 is used to calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship; the second calculation unit 304 is used to use the standard smoke puffs to multiplex and calculate the smoke puff group concentration matrix of several ships at the specified observation point in the observation space according to the concentration matrix.

[0088] In one embodiment, the first calculation unit 303 is used to calculate the standard smoke puff of each ship at each time point within the observation period based on the concentration resolution, and fill the initialized concentration matrix with longitude, latitude, and concentration vector in sequence to obtain the concentration matrix corresponding to each ship.

[0089] In one embodiment, the second calculation unit 304 includes a diffusion concentration matrix calculation subunit and a superposition subunit.

[0090] The diffusion concentration matrix calculation subunit is used to determine the concentration matrix corresponding to the ship at a certain observation point, and calculate the diffusion concentration matrix by multiple translation and reuse of standard smoke puffs; the superposition subunit is used to superimpose all diffusion concentration matrices to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0091] In one embodiment, the superposition subunit includes a series connection module, a resolution setting module, a unification module and a merging module.

[0092] A series connection module is used to connect all diffusion concentration matrices in series to form an intermediate matrix; a resolution setting module is used to set the longitude and latitude resolution; a unification module is used to unify all longitude and latitude coordinates in the intermediate matrix to the longitude and latitude resolution; a merging module is used to merge the rows with the same longitude and latitude values ​​in the unified intermediate matrix into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0093] It should be noted that those skilled in the art can clearly understand that the specific implementation process of the above-mentioned gas diffusion translation calculation device 300 and each unit can refer to the corresponding description in the aforementioned method embodiment, and for the convenience and brevity of description, it will not be repeated here.

[0094] The gas diffusion translation calculation device 300 can be implemented in the form of a computer program. The computer program can be used in Figure 8 Runs on the computer device shown.

[0095] See also Figure 8 , Figure 8 5 is a schematic block diagram of a computer device provided in an embodiment of the present application. The computer device 500 may be a server, wherein the server may be an independent server or a server cluster composed of multiple servers.

[0096] See also Figure 8The computer device 500 includes a processor 502 , a memory and a network interface 505 connected via a system bus 501 , wherein the memory may include a non-volatile storage medium 503 and an internal memory 504 .

[0097] The non-volatile storage medium 503 can store an operating system 5031 and a computer program 5032. The computer program 5032 includes program instructions, and when the program instructions are executed, the processor 502 can execute a gas diffusion translation calculation method.

[0098] The processor 502 is used to provide computing and control capabilities to support the operation of the entire computer device 500 .

[0099] The internal memory 504 provides an environment for the operation of the computer program 5032 in the non-volatile storage medium 503. When the computer program 5032 is executed by the processor 502, the processor 502 can execute a gas diffusion translation calculation method.

[0100] The network interface 505 is used to communicate with other devices over the network. Figure 8 The structure shown in the figure is merely a block diagram of a portion of the structure related to the solution of the present application, and does not constitute a limitation on the computer device 500 to which the solution of the present application is applied. The specific computer device 500 may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0101] The processor 502 is used to run the computer program 5032 stored in the memory to implement the following steps:

[0102] Set the gas diffusion related parameters of each standard smoke puff; initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction; calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation duration to obtain the concentration matrix corresponding to each ship; use the standard smoke puffs to shift and reuse multiple times, and calculate the smoke puff group concentration matrix of several ships at the specified observation point in the observation space according to the concentration matrix.

[0103] The gas diffusion related parameters include wind direction, wind speed, emission amount, emission point and diffusion factor.

[0104] In one embodiment, the processor 502 specifically implements the following steps when implementing the step of calculating the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship:

[0105] Based on the concentration resolution, the standard smoke puffs of each ship at each time point within the observation time are calculated and filled into the initialized concentration matrix in sequence with longitude, latitude, and concentration vector to obtain the concentration matrix corresponding to each ship.

[0106] Wherein, the concentration matrix is C i represents the concentration of m spatial points observed at observation time i, where is the spatial coordinate location of observation point j corresponding to time point i, is the coordinate point The corresponding smoke puff diffusion concentration, j∈(1,2,…,m).

[0107] In one embodiment, when the processor 502 implements the step of calculating the smoke puff group concentration matrix of a plurality of ships at a specified observation point in the observation space according to the concentration matrix by using the standard smoke puff multiple translation multiplexing method, the processor 502 specifically implements the following steps:

[0108] Determine the concentration matrix corresponding to the ship at a certain observation point, and calculate the diffusion concentration matrix by multiple translation and reuse of standard smoke puffs; superimpose all diffusion concentration matrices to obtain the smoke puff group concentration matrix of several ships at the specified observation point in the observation space.

[0109] Among them, the diffusion concentration matrix is in, is the coordinate of the instantaneous emission point of the standard puff, are the coordinates of the instantaneous emission point of the target puff, and (Δp1, Δp2) are the translation of the target puff relative to the reference puff.

[0110] In one embodiment, when the processor 502 implements the step of superimposing all diffusion concentration matrices to obtain the smoke group concentration matrices of a plurality of ships at a specified observation point in the observation space, the processor 502 specifically implements the following steps:

[0111] All diffusion concentration matrices are connected in series to form an intermediate matrix; the longitude and latitude resolution is set; all longitude and latitude coordinates in the intermediate matrix are unified to the longitude and latitude resolution; the rows with the same longitude and latitude values ​​in the unified intermediate matrix are merged into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0112] It should be understood that in the embodiment of the present application, the processor 502 may be a central processing unit (CPU), and the processor 502 may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0113] It can be understood by those skilled in the art that all or part of the processes in the method for implementing the above embodiment can be completed by instructing the relevant hardware through a computer program. The computer program includes program instructions, and the computer program can be stored in a storage medium, which is a computer-readable storage medium. The program instructions are executed by at least one processor in the computer system to implement the process steps of the embodiment of the above method.

[0114] Therefore, the present invention also provides a storage medium. The storage medium may be a computer-readable storage medium. The storage medium stores a computer program, wherein when the computer program is executed by a processor, the processor executes the following steps:

[0115] Set the gas diffusion related parameters of each standard smoke puff; initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction; calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation duration to obtain the concentration matrix corresponding to each ship; use the standard smoke puffs to shift and reuse multiple times, and calculate the smoke puff group concentration matrix of several ships at the specified observation point in the observation space according to the concentration matrix.

[0116] The gas diffusion related parameters include wind direction, wind speed, emission amount, emission point and diffusion factor.

[0117] In one embodiment, when the processor executes the computer program to implement the step of calculating the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship, the processor specifically implements the following steps:

[0118] Based on the concentration resolution, the standard smoke puffs of each ship at each time point within the observation time are calculated and filled into the initialized concentration matrix in sequence with longitude, latitude, and concentration vector to obtain the concentration matrix corresponding to each ship.

[0119] Wherein, the concentration matrix is C i represents the concentration of m spatial points observed at observation time i, where is the spatial coordinate location of observation point j corresponding to time point i, is the coordinate point The corresponding smoke puff diffusion concentration, j∈(1,2,…,m).

[0120] In one embodiment, when the processor executes the computer program to implement the step of calculating the smoke group concentration matrix of several ships at a specified observation point in the observation space according to the concentration matrix by using the standard smoke puff multiple translation multiplexing method, the processor specifically implements the following steps:

[0121] Determine the concentration matrix corresponding to the ship at a certain observation point, and calculate the diffusion concentration matrix by multiple translation and reuse of standard smoke puffs; superimpose all diffusion concentration matrices to obtain the smoke puff group concentration matrix of several ships at the specified observation point in the observation space.

[0122] Among them, the diffusion concentration matrix is in, is the coordinate of the instantaneous emission point of the standard puff, are the coordinates of the instantaneous emission point of the target puff, and (Δp1, Δp2) are the translation of the target puff relative to the reference puff.

[0123] In one embodiment, when the processor executes the computer program to implement the step of superimposing all diffusion concentration matrices to obtain smoke group concentration matrices of a number of ships at a specified observation point in the observation space, the processor specifically implements the following steps:

[0124] All diffusion concentration matrices are connected in series to form an intermediate matrix; the longitude and latitude resolution is set; all longitude and latitude coordinates in the intermediate matrix are unified to the longitude and latitude resolution; the rows with the same longitude and latitude values ​​in the unified intermediate matrix are merged into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

[0125] The storage medium may be a USB flash drive, a mobile hard disk, a read-only memory (ROM), a magnetic disk, or an optical disk, etc., which are computer-readable storage media that can store program codes.

[0126] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0127] In the 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 only schematic. For example, the division of each unit is only a logical function division, and there may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed.

[0128] The steps in the method of the embodiment of the present invention can be adjusted in order, combined and deleted according to actual needs. The units in the device of the embodiment of the present invention can be combined, divided and deleted according to actual needs. In addition, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0129] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, terminal, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention.

[0130] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed by the present invention, and these modifications or replacements should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be based on the protection scope of the claims.

Claims

1. Gas diffusion translation calculation method, characterized in that: include: Set gas diffusion related parameters for each standard puff; Initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction; Calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation period to obtain the concentration matrix corresponding to each ship; The smoke puff group concentration matrix of several ships at the designated observation point in the observation space is calculated according to the concentration matrix by adopting the method of multiple translation and reuse of standard smoke puffs.

2. The gas diffusion translation calculation method according to claim 1, characterized in that: The gas diffusion related parameters include wind direction, wind speed, emission amount, emission point and diffusion factor.

3. The gas diffusion translation calculation method according to claim 2, characterized in that: The calculation of the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship includes: Based on the concentration resolution, the standard smoke puffs of each ship at each time point within the observation time are calculated and filled into the initialized concentration matrix in sequence with longitude, latitude, and concentration vector to obtain the concentration matrix corresponding to each ship.

4. The gas diffusion translation calculation method according to claim 3, characterized in that: The concentration matrix is C i represents the concentration of m spatial points observed at observation time i, where is the spatial coordinate location of observation point j corresponding to time point i, is the coordinate point The corresponding smoke puff diffusion concentration, j∈(1,2,…,m).

5. The gas diffusion translation calculation method according to claim 4, characterized in that: The method of using multiple translation and reuse of standard smoke puffs to calculate the smoke puff group concentration matrix of several ships at a designated observation point in the observation space according to the concentration matrix includes: Determine the concentration matrix corresponding to the ship at a certain observation time point, and calculate the diffusion concentration matrix by multiple translation and reuse of standard smoke puffs; All diffusion concentration matrices are superimposed to obtain the smoke group concentration matrix of several ships at a specified observation point in the observation space.

6. The gas diffusion translation calculation method according to claim 5, characterized in that: The diffusion concentration matrix is in, is the coordinate of the instantaneous emission point of the standard puff, are the coordinates of the instantaneous emission point of the target puff, and (Δp1, Δp2) are the translation of the target puff relative to the reference puff.

7. The gas diffusion translation calculation method according to claim 6, characterized in that: All diffusion concentration matrices are superimposed to obtain the smoke group concentration matrix of several ships at a specified observation point in the observation space, including: All diffusion concentration matrices are connected in series to form an intermediate matrix; Set the latitude and longitude resolution; Unifying all longitude and latitude coordinates in the intermediate matrix to the longitude and latitude resolution; The rows with the same longitude and latitude values ​​in the unified intermediate matrix are merged into one row, and the concentration value in the merged row is set to the sum of the concentration values ​​in the rows with the same longitude and latitude values ​​in the unified intermediate matrix, so as to obtain the smoke group concentration matrix of several ships at the specified observation point in the observation space.

8. A gas diffusion translation calculation device, characterized in that: include: A setting unit, used to set gas diffusion related parameters of each standard puff; Initialization unit, used to initialize the observation space, observation duration, observation time point, and calculate the downwind direction and crosswind direction; The first calculation unit is used to calculate the concentration matrix sequence of the standard smoke puffs of each ship at each time point within the observation time to obtain the concentration matrix corresponding to each ship; The second calculation unit is used to calculate the smoke group concentration matrix of several ships at the designated observation point in the observation space according to the concentration matrix by using the standard smoke puff multiple translation multiplexing method.

9. A computer device, characterized in that: The computer device comprises a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images