CN115902142A - Water body point source tracing method based on improved ensemble Kalman filtering algorithm - Google Patents

Water body point source tracing method based on improved ensemble Kalman filtering algorithm Download PDF

Info

Publication number
CN115902142A
CN115902142A CN202211630910.9A CN202211630910A CN115902142A CN 115902142 A CN115902142 A CN 115902142A CN 202211630910 A CN202211630910 A CN 202211630910A CN 115902142 A CN115902142 A CN 115902142A
Authority
CN
China
Prior art keywords
pollution source
model
parameters
assimilation
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211630910.9A
Other languages
Chinese (zh)
Inventor
孔俊
荆立
杨悦知
江朝华
洪淑娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN202211630910.9A priority Critical patent/CN115902142A/en
Publication of CN115902142A publication Critical patent/CN115902142A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a water point source tracing method based on an improved ensemble Kalman filtering algorithm, which comprises the following steps: investigating the background of pollutants in a research area where leakage accidents may occur; arranging pollutant concentration monitoring points according to the field conditions of the research area; according to the terrain and hydrological data of a research area, constructing a hydrodynamic model and a substance transport model of the area; improving an ensemble Kalman filtering algorithm, coupling and nesting the algorithm and a substance transportation model, and constructing a pollution source parameter assimilation model based on the coupling algorithm; according to the pollutant concentration data of the monitoring points transmitted back remotely in real time, carrying out assimilation calculation on pollution source parameters by adopting an assimilation model; and correcting the pollution source information according to the assimilation result, determining the assimilation value of the pollution source parameter when the parameter assimilation result tends to be stable, and predicting the pollutant propagation path in real time by using the assimilated pollution source parameter. The method can more quickly and accurately capture the position of the pollution source and identify the initial release time and release amount of the pollution source.

Description

Water body point source tracing method based on improved ensemble Kalman filtering algorithm
Technical Field
The invention relates to the technical field of surface water pollution source parameter optimization identification, in particular to a water point source tracing method based on an improved ensemble Kalman filtering algorithm.
Background
When an environmental pollution accident occurs suddenly, pollution source information such as the position, release time and release amount of the pollution source is rapidly and accurately obtained, which is very critical for the management, control, risk assessment and remediation of the accident. Because the timeliness and the urgency are very important to solve the problem of tracing the source of the pollutants, the efficiency and the precision of the method in identifying the parameters of the pollution source need to be concerned when the tracing method is applied. The traditional tracing method needs to distribute a large number of monitoring points on site and survey and research on the spot, has huge consumption on manpower and material resources, is lower in source searching efficiency, and cannot take remedial measures to the source in time. With the continuous upgrading of monitoring equipment, the remote transmission of monitoring data is popularized, and meanwhile, an optimization algorithm is matched, so that the pollution source can be quickly captured and identified under the condition of fewer monitoring points.
Disclosure of Invention
The invention aims to solve the technical problem of providing a water body point source tracing method based on an improved integrated Kalman filtering algorithm, which can more quickly and accurately capture the position of a pollution source and identify the initial release time and release amount of the pollution source, and provides technical support for the treatment work of the pollution source.
In order to solve the technical problem, the invention provides a water point source tracing method based on an improved ensemble Kalman filtering algorithm, which comprises the following steps:
step 1, investigating a pollutant background of a possible leakage accident in a research area;
step 2, reasonably arranging pollutant concentration monitoring points according to the field conditions of the research area;
step 3, constructing a hydrodynamic model and a substance transportation model of the region according to the terrain and hydrological data of the research region;
step 4, improving an ensemble Kalman filtering algorithm, coupling and nesting the algorithm and a substance transport model, and constructing a pollution source parameter assimilation model based on the coupling algorithm;
step 5, carrying out assimilation calculation on pollution source parameters by adopting an assimilation model according to the pollutant concentration data of the monitoring point transmitted remotely in real time;
and 6, dynamically correcting the pollution source information according to each synchronization result, determining the assimilation value of the pollution source parameter when the parameter assimilation result tends to be stable, and predicting the pollutant propagation path in real time by using the assimilated pollution source parameter.
Preferably, in step 1, the background of the pollutants in the investigation region, where a leakage accident may occur, is a potential pollutant type and risk hazard near the investigation region, so that a corresponding pollution monitoring device can be selected conveniently.
Preferably, in the step 2, pollutant concentration monitoring points including the arrangement positions and the number of the monitoring points are reasonably arranged according to the field conditions of the research area, if the flow field of the research area is unidirectional flow, only one or two monitoring points need to be arranged at the downstream, if the flow field of the research area is reciprocating flow, a plurality of monitoring points are arranged in the research area, and specifically, the monitoring points can be arranged near the potential pollution source area according to investigation or at equal intervals in the longitudinal direction or the transverse direction, so that the comprehensiveness of data acquisition is ensured.
Preferably, in step 3, the hydrodynamic model is constructed according to the terrain and water depth of the research area, the model is a planar two-dimensional model, the model calculation adopts an ELCIRC open source program based on Fortran language, the model grid adopts an unstructured grid, the grid dimension can be generally set to 10m-50m, and the local area can be encrypted according to the actual situation.
Preferably, in step 3, the material transport model is a planar two-dimensional model, and the control equation is represented by formula (1):
Figure BDA0004005725940000021
in the formula: c is the concentration of the substance; t is time; u, v are the flow velocity components in the x, y directions, respectively; h is water depth; k is an x, y spatial two-dimensional diffusion tensor whose value is represented by the following equation (2):
K xx =K L cos 2 (θ)+K T sin 2 (θ)(2a)
K xy =K xy =(K L -K T )cos(θ)sin(θ)(2b)
K yy =K L sin 2 9θ)+K T cos 2 (θ) (2 c) wherein: theta is an included angle between the flow direction and the x axis; k L 、K T Longitudinal and transverse diffusion coefficients, respectively, which are expressed in relation to the flow velocity by the following equation (3)
Figure BDA0004005725940000022
Figure BDA0004005725940000023
In the formula: c is the metabolization capacity coefficient; alpha and beta are respectively constant coefficients; g is the acceleration of gravity;
separating the convection term and the diffusion term in the formula (1) into two parts by an operator splitting method,
Figure BDA0004005725940000024
Figure BDA0004005725940000025
and (4) respectively carrying out differential solution on the formula by adopting a finite volume method, and finally obtaining the concentration of the released substances in the research area.
Preferably, in step 3, the mesh of the material transport model is the same as that of the hydrodynamic model, and the flow field of the material transport model is provided by calculation of the hydrodynamic model and used for driving material transport.
Preferably, in step 4, the assimilated parameters of the pollution source in the pollution source parameter assimilation model are respectively a pollution source position, an initial release time of the pollution source and a total release amount of the pollution source.
Preferably, in step 4, the improved ensemble kalman filtering algorithm is constructed by integrating correlation coefficients based on a standard ensemble kalman filtering algorithm, and the construction steps are as follows:
(1) Generating an initial set; respectively generating sets containing Nr members aiming at four pollution source parameters, randomly selecting initial values of the set members from independent uniform distribution, determining the upper limit and the lower limit of the uniform distribution by the acquired prior information, and representing a pollution source parameter vector by the following formula:
Figure BDA0004005725940000031
in the formula:
Figure BDA0004005725940000032
representing a pollution source parameter vector, wherein i is a number realized by a set, a superscript T represents a matrix transposition, and X, Y, T and M in the vector are respectively a horizontal coordinate, a vertical coordinate, initial release time and release amount of a pollution source space plane position;
(2) Predicting; putting the latest pollution source parameter value obtained at an observation time t-1 into a material transportation model, starting calculation from an initial time 0, stopping calculation until the observation time t, obtaining a simulation state of the observation time t at the moment, and updating the state from the time 0 to the time t-1, wherein the simulation state parameter refers to the pollutant concentration, and the prediction of the state variable is as follows:
Figure BDA0004005725940000033
in the formula: c 0 Is a state initial concentration profile; superscripts a and f represent predicted values and updated values, respectively; psi is a state transition operator, referred to herein as a material transport model;
for equation (6), a correlation coefficient R between the predicted concentration and the observed concentration in the observed time series is introduced as a state variable instead of the contaminant concentration, the correlation coefficient R being represented by the following equation:
Figure BDA0004005725940000034
in the formula: subscript j is a time series number; superscripts s and ob represent concentration analog and observed values respectively,
Figure BDA0004005725940000035
and &>
Figure BDA0004005725940000036
Are respectively C s And C ob Average value of (d);
therefore, the state prediction of equation (6) is represented by the following equation:
Figure BDA0004005725940000037
in the formula: f is a
Figure BDA0004005725940000038
Switch over to->
Figure BDA0004005725940000039
The conversion operator of (3);
due to the introduction of the correlation coefficient R, according to the migration characteristic of the pollutants, when the position and the release time of the pollution source are determined, if other situations remain unchanged, the pollutant release amount has no influence on the magnitude of the correlation coefficient, and the pollutant release amount M is separated from the pollution source parameter to be identified; pollution source parameter vector to be assimilated
Figure BDA0004005725940000041
Only two parameters of the contamination source position and the initial release time are included;
Figure BDA0004005725940000042
(3) Updating; before updating, normal score conversion is carried out on parameters, the parameters are selected from uniform distribution, the optimal parameter can be achieved by the Kalman filtering method when the parameters obey Gaussian distribution, the normal score conversion equation needs to be recalculated after each prediction step, and the conversion of system parameters is shown as the following formula:
Figure BDA0004005725940000043
in the formula: phi (phi) of t A conversion form of normal score at time t, the conversion being related to the cumulative distribution function; the upper line indicates the parameters after conversion of the normal score.
Then, updating the parameters; assuming that the number of observation points is Nob, each realization in the set is updated by the following equation;
Figure BDA0004005725940000044
in the formula:
Figure BDA0004005725940000045
and &>
Figure BDA0004005725940000046
Respectively updating parameters and predicting vectors after the conversion of the normal score at the moment t; r is t The autocorrelation coefficient of the observed pollutant concentration from time 0 to t; e.g. of the type t Is r t The error vector of (2); g t Is a Kalman gain matrix, and the expression is as follows:
Figure BDA0004005725940000047
in the formula: r t Is an error vector e t The covariance matrix of (a);
Figure BDA0004005725940000048
is contaminated byDye source parameter vector>
Figure BDA0004005725940000049
And predicted correlation coefficient>
Figure BDA00040057259400000410
The covariance of (a); />
Figure BDA00040057259400000411
Is the relevant coefficient->
Figure BDA00040057259400000412
An auto-covariance;
updated transformation source parameters
Figure BDA00040057259400000413
Requires a reverse conversion into ∑ in the prediction of step (2)>
Figure BDA00040057259400000414
By means of a conversion inverse function related to the cumulative distribution inverse function->
Figure BDA00040057259400000415
The assimilation vector at time t->
Figure BDA00040057259400000416
The parameters that contain the following final updates:
Figure BDA00040057259400000417
preferably, in the step 5, the assimilation calculation of the pollution source parameters comprises two parts, wherein one part is the assimilation calculation of the pollution source parameters including the position and the initial release time of the pollution source, and the other part is the parameter optimization calculation of the total release amount of the separated pollutants; assimilating calculation is carried out on the pollution source parameters including the position of the pollution source and the initial release time according to an expression (5) -an expression (13), and the optimization of the pollutant release total quantity parameters is calculated according to the following expression:
Figure BDA0004005725940000051
in the formula: RMSE is the average error between the process of simulating concentration and the process of observing concentration, i is the step length of observation time,
Figure BDA0004005725940000052
and &>
Figure BDA0004005725940000053
Respectively an observed concentration and a simulated concentration at an observation time step i, and n is a concentration time process data volume.
Preferably, in step 6, dynamically correcting the pollution source information is performed based on the data of each observation step, the corrected pollution source information is the position of the pollution source and the initial release time, the optimization calculation is performed according to the formula (14) when the assimilation results of the two parameters tend to be stable, and finally the optimized pollution source parameters are obtained to perform the real-time prediction of the pollutant propagation path.
The invention has the beneficial effects that: the invention provides a tracing method suitable for the point source of the surface water body on the basis of the transportation rule of pollutants released by the point source in the surface water body, by introducing relevant coefficient factors into an ensemble Kalman filtering method and coupling an improved ensemble Kalman filtering algorithm with a surface water substance transportation model.
Drawings
FIG. 1 is a schematic flow chart of the method of the present invention.
FIG. 2 is a schematic flow chart of a method of an improved ensemble Kalman filtering algorithm of the present invention.
Fig. 3 is a schematic diagram of a source tracing and assimilating result based on an improved ensemble kalman filtering algorithm.
Fig. 4 is a schematic diagram of a source tracing and assimilating result based on an improved ensemble kalman filtering algorithm.
Fig. 5 is a schematic view of a traceability homogenization result based on an improved ensemble kalman filter algorithm according to the present invention.
Fig. 6 is a schematic view of a traceability homogenization result based on an improved ensemble kalman filter algorithm according to the present invention.
Detailed Description
As shown in fig. 1, a water point source tracing method based on an improved ensemble kalman filtering algorithm includes the following steps:
step 1, investigating a pollutant background of a possible leakage accident in a research area;
step 2, reasonably arranging pollutant concentration monitoring points according to the field conditions of the research area;
step 3, constructing a regional hydrodynamic model and a substance transportation model according to the terrain and hydrological data of the research region;
step 4, improving an ensemble Kalman filtering algorithm, coupling and nesting the algorithm and a substance transport model, and constructing a pollution source parameter assimilation model based on the coupling algorithm;
step 5, carrying out assimilation calculation on pollution source parameters by adopting an assimilation model according to the pollutant concentration data of the monitoring point transmitted remotely in real time, wherein the assimilation calculation step is shown in figure 2;
and 6, dynamically correcting the pollution source information according to each synchronization result, determining the assimilation value of the pollution source parameter when the parameter assimilation result tends to be stable, and predicting the pollutant propagation path in real time by using the assimilated pollution source parameter.
Example (b):
the water point source tracing method based on the improved ensemble Kalman filtering algorithm comprises the following steps:
step 1: the background of pollutants in a possible leakage accident in a research area is investigated, the research area is a square area of 500m by 500m in the embodiment, the real position of a pollution source is (250 ), the initial release time is 100s, the total pollutant release amount is 10kg, and the pollutant index is COD.
And 2, step: pollutant concentration monitoring points are reasonably arranged according to the field conditions of the research area, in the embodiment, 3 monitoring points are arranged, the positions are (290, 260), (300, 255) and (330, 270), the pollutant concentration data of the observation points contain 20% of random errors, the observation time is 130 s-300 s, and the observation interval is 1s.
And step 3: and constructing a hydrodynamic model and a substance transportation model of the region according to the terrain and hydrological data of the research region. The hydrodynamic model provides a hydrodynamic field for the material transport model.
And 4, step 4: and improving an ensemble Kalman filtering algorithm, coupling and nesting the algorithm and a substance transport model, and constructing a pollution source parameter assimilation model based on the coupling algorithm.
And 5: and carrying out assimilation calculation on the pollution source parameters by adopting an assimilation model according to the pollutant concentration data of the monitoring point transmitted back in real time in a remote mode. As shown in the box diagrams of fig. 3 to 5, assimilation data of two parameters, i.e., the source position (abscissa and ordinate points) and the initial release time of the contamination source at the calculation time steps of 11 th, 20 th, 40 th, 60 th, 80 th, 100 th, 120 th, 140 th, and 160 th, respectively, are selected, and the central horizontal line of the box in the diagrams is the median of the set, it can be seen that the assimilation parameters tend to converge stably at the calculation time step of 100 th, and the assimilation results can maintain accuracy even when the contamination concentration data contains a random error of 20%.
And 6: and dynamically correcting the pollution source information according to each synchronization result, determining the assimilation value of the pollution source parameter when the parameter assimilation result tends to be stable, and performing real-time prediction on the pollutant propagation path by using the assimilated pollution source parameter. And selecting the position (horizontal and vertical coordinates) of the pollution source at the 100 th calculation time step and the assimilation set average value of the initial release time of the pollution source as an assimilation result, optimizing the parameter of the total pollutant release amount based on the result, wherein fig. 6 shows the optimization result of the total pollutant release amount, and the release value corresponding to the lowest RMSE value is selected as the final optimization value of the parameter. And performing forward simulation by using a material transport model according to the optimized pollution source parameters, and predicting a subsequent pollutant propagation path in real time.

Claims (10)

1. A water point source tracing method based on an improved ensemble Kalman filtering algorithm is characterized by comprising the following steps:
step 1, investigating a pollutant background of a possible leakage accident in a research area;
step 2, reasonably arranging pollutant concentration monitoring points according to the field conditions of the research area;
step 3, constructing a regional hydrodynamic model and a substance transportation model according to the terrain and hydrological data of the research region;
step 4, improving an ensemble Kalman filtering algorithm, coupling and nesting the algorithm and a substance transport model, and constructing a pollution source parameter assimilation model based on the coupling algorithm;
step 5, carrying out assimilation calculation on pollution source parameters by adopting an assimilation model according to the pollutant concentration data of the monitoring point transmitted back in real time in a remote mode;
and 6, dynamically correcting the pollution source information according to each synchronization result, determining the assimilation value of the pollution source parameter when the parameter assimilation result tends to be stable, and predicting the pollutant propagation path in real time by using the assimilated pollution source parameter.
2. The water body point source tracing method based on the improved ensemble Kalman filtering algorithm according to claim 1, wherein in step 1, the background of the pollutants possibly causing leakage accidents in the investigation research area is potential pollutant types and risk hazards near the investigation research area, and corresponding pollution monitoring equipment is convenient to select.
3. The water body point source tracing method based on the improved ensemble Kalman filtering algorithm according to claim 1, characterized in that in step 2, pollutant concentration monitoring points are reasonably arranged according to field conditions of a research area, including the arrangement positions and the number of the monitoring points, if the flow field of the research area is unidirectional flow, only one or two monitoring points are arranged at downstream, if the flow field of the research area is reciprocating flow, a plurality of monitoring points are arranged in the research area, and the comprehensiveness of acquired data is ensured according to the arrangement near a potential pollution source area or according to longitudinal or transverse equidistance.
4. The water body point source tracing method based on the improved ensemble Kalman filtering algorithm according to claim 1, characterized in that in step 3, the hydrodynamic model is constructed according to the terrain and water depth of the research area, and is a planar two-dimensional model, the model calculation adopts an ELCIRC source opening program based on Fortran language, the model mesh adopts an unstructured mesh, the mesh size is set to 10m-50m, and the local area is encrypted according to the actual situation.
5. The water body point source tracing method based on the improved ensemble Kalman filtering algorithm according to claim 1, wherein in step 3, the material transport model is a planar two-dimensional model, and the control equation is expressed by formula (1):
Figure FDA0004005725930000011
in the formula: c is the concentration of the substance; t is time; u and v are flow velocity components in the x and y directions respectively; h is water depth; k is an x, y spatial two-dimensional diffusion tensor whose value is represented by the following equation (2):
K xx =K L cos 2 (θ)+K T sin 2 (θ)(2a)K xy =K xy =(K L -K T )cos(θ)sin(θ)(2b)K yy =K L sin 2 (θ)+K T cos 2 (θ) (2 c) wherein: theta is an included angle between the flow direction and the x axis; k L 、K T Longitudinal and transverse diffusion coefficients, respectively, which are expressed in relation to the flow velocity by the following equation (3)
Figure FDA0004005725930000021
Figure FDA0004005725930000022
In the formula: c is the metabolization capacity coefficient; α, α) are constant coefficients, respectively;
separating the convection term and the diffusion term in the formula (1) into two parts by an operator splitting method,
Figure FDA0004005725930000023
/>
Figure FDA0004005725930000024
and (4) respectively carrying out differential solution on the formula by adopting a finite volume method, and finally obtaining the concentration of the released substances in the research area.
6. The improved ensemble kalman filter algorithm-based water point source tracing method according to claim 1, wherein in step 3, the mesh of the material transport model is the same as that of the hydrodynamic model, and the flow field of the material transport model is provided by calculation of the hydrodynamic model and used for driving material transport.
7. The method for tracing the source of the water body point based on the improved ensemble kalman filter algorithm according to claim 1, wherein in step 4, the assimilated parameters of the pollution source in the pollution source parameter assimilation model are the position of the pollution source, the initial release time of the pollution source and the total release amount of the pollution source, respectively.
8. The water body point source tracing method based on the improved Kalman filtering algorithm of claim 1, wherein in the step 4, the improved Kalman filtering algorithm is constructed by integrating correlation coefficients based on a standard Kalman filtering algorithm, and the construction steps are as follows:
(1) Generating an initial set; respectively generating sets containing Nr members aiming at three pollution source parameters, randomly selecting initial values of the set members from independent uniform distribution, determining the upper limit and the lower limit of the uniform distribution by the acquired prior information, and expressing a pollution source parameter vector by the following formula:
Figure FDA0004005725930000025
in the formula:
Figure FDA0004005725930000026
representing a pollution source parameter vector, wherein i is a number realized by a set, and superscript T represents matrix transposition;
(2) Predicting; putting the latest pollution source parameter value obtained at an observation time t-1 into a material transportation model, starting calculation from an initial time 0, stopping calculation until the observation time t, obtaining a simulation state of the observation time t at the moment, and updating the state from the time 0 to the time t-1, wherein the simulation state parameter refers to the pollutant concentration, and the prediction of the state variable is as follows:
Figure FDA0004005725930000031
in the formula: c 0 Is a state initial concentration profile; superscripts a and f represent predicted values and updated values, respectively; psi is a state transition operator, here a material transport model;
for equation (6), a correlation coefficient R between the predicted concentration and the observed concentration in the observed time series is introduced as a state variable instead of the contaminant concentration, the correlation coefficient R being represented by the following equation:
Figure FDA0004005725930000032
in the formula: subscript j is a time series number; superscripts s and ob represent the simulated and observed values respectively,
Figure FDA0004005725930000033
and &>
Figure FDA0004005725930000034
Are respectively C s And C ob Average value of (d);
therefore, the state prediction of equation (6) is represented by the following equation:
Figure FDA0004005725930000035
in the formula: f is a
Figure FDA0004005725930000036
Switch over to->
Figure FDA0004005725930000037
The conversion operator of (3);
due to the introduction of the correlation coefficient R, according to the migration characteristic of the pollutants, after the position and the release time of the pollution source are determined, if other situations remain unchanged, the pollutant release amount has no influence on the magnitude of the correlation coefficient, and the pollutant release amount M is separated from the parameters of the pollution source to be identified independently; pollution source parameter vector to be assimilated
Figure FDA0004005725930000038
Only two parameters of the position of the pollution source and the initial release time are included; />
Figure FDA0004005725930000039
(3) Updating; before updating, normal score conversion is carried out on parameters, the parameters are selected from uniform distribution, the optimal parameter can be achieved by the Kalman filtering method when the parameters obey Gaussian distribution, the normal score conversion equation needs to be recalculated after each prediction step, and the conversion of system parameters is shown as the following formula:
Figure FDA00040057259300000310
in the formula: phi t A conversion form of normal score at time t, the conversion being related to the cumulative distribution function;
then, updating the parameters; assuming that the number of observation points is Nob, each realization in the set is updated by the following equation;
Figure FDA0004005725930000041
in the formula:
Figure FDA0004005725930000042
and &>
Figure FDA0004005725930000043
Respectively updating parameters and predicting vectors after normal score conversion at the time t; r is t Is the observed contaminant concentration autocorrelation coefficient from time 0 to t; e.g. of the type t Is r t The error vector of (2); g t Is a Kalman gain matrix, and the expression is as follows:
Figure FDA0004005725930000044
in the formula: r t Is an error vector e t The covariance matrix of (a);
Figure FDA0004005725930000045
for the contamination source parameter vector->
Figure FDA0004005725930000046
And predicted correlation coefficient R f The covariance of (a); />
Figure FDA0004005725930000047
Is a correlation coefficient R f An auto-covariance;
updated translation source parameters
Figure FDA0004005725930000048
Requires a reverse conversion into &'s in step (2) prediction>
Figure FDA0004005725930000049
By means of a conversion inverse function related to the cumulative distribution inverse function->
Figure FDA00040057259300000410
Assimilation vector at time t>
Figure FDA00040057259300000411
The parameters that contain the following final updates:
Figure FDA00040057259300000412
9. the water body point source tracing method based on the improved ensemble Kalman filtering algorithm according to claim 1, characterized in that in step 5, the pollution source parameter assimilation calculation includes two parts, one part is the pollution source parameter assimilation calculation including the pollution source position and the initial release time, and the other part is the parameter optimization calculation of the total release amount of the separated pollutants; assimilating calculation is carried out on the pollution source parameters including the position of the pollution source and the initial release time according to an expression (5) -an expression (13), and the optimization of the pollutant release total quantity parameters is calculated according to the following expression:
Figure FDA00040057259300000413
in the formula: i is the observation time step size,
Figure FDA00040057259300000414
and &>
Figure FDA00040057259300000415
Respectively, an observed concentration and a simulated concentration at an observation time step i.
10. The water body point source tracing method based on the improved ensemble kalman filter algorithm as claimed in claim 1, wherein in step 6, the dynamic correction of the pollution source information is performed based on the data of each observation step, the corrected pollution source information is the position of the pollution source and the initial release time, the total amount of the pollution source release is optimized and calculated according to the formula (14) when the assimilation results of the two parameters tend to be stable, and finally the optimized pollution source parameters are obtained for the real-time prediction of the pollutant propagation path.
CN202211630910.9A 2022-12-19 2022-12-19 Water body point source tracing method based on improved ensemble Kalman filtering algorithm Pending CN115902142A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211630910.9A CN115902142A (en) 2022-12-19 2022-12-19 Water body point source tracing method based on improved ensemble Kalman filtering algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211630910.9A CN115902142A (en) 2022-12-19 2022-12-19 Water body point source tracing method based on improved ensemble Kalman filtering algorithm

Publications (1)

Publication Number Publication Date
CN115902142A true CN115902142A (en) 2023-04-04

Family

ID=86496118

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211630910.9A Pending CN115902142A (en) 2022-12-19 2022-12-19 Water body point source tracing method based on improved ensemble Kalman filtering algorithm

Country Status (1)

Country Link
CN (1) CN115902142A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116739386A (en) * 2023-08-10 2023-09-12 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) Multi-index fusion pollution tracing method, equipment and readable storage medium

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116739386A (en) * 2023-08-10 2023-09-12 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) Multi-index fusion pollution tracing method, equipment and readable storage medium
CN116739386B (en) * 2023-08-10 2024-03-08 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) Multi-index fusion pollution tracing method, equipment and readable storage medium

Similar Documents

Publication Publication Date Title
Pecchi et al. Species distribution modelling to support forest management. A literature review
Wang et al. Ecological drivers of spatial community dissimilarity, species replacement and species nestedness across temperate forests
CN104899298B (en) A kind of microblog emotional analysis method based on large-scale corpus feature learning
Sihag et al. Modeling unsaturated hydraulic conductivity by hybrid soft computing techniques
CN115902142A (en) Water body point source tracing method based on improved ensemble Kalman filtering algorithm
CN111178786B (en) Emission source position determining method and system for guaranteeing regional air quality
van Noordwijk et al. Concepts and methods for studying interactions of roots and soil structure
CN113267607A (en) Characteristic parameter identification system for field organic pollutant migration process
Viskari et al. Improving Yasso15 soil carbon model estimates with ensemble adjustment Kalman filter state data assimilation
CN108399313A (en) A kind of underground engineering injures Simulation & evaluation system and method
Safitri Comparison of students science process skills after using learning an experimental and virtual laboratory on Archimedes Laws
Albo et al. The Aerodyne Inverse Modeling System (AIMS): Source estimation applied to the FFT 07 experiment and to simulated mobile sensor data
Diop et al. The use of ALADIN model and MERRA-2 reanalysis to represent dust seasonal dry deposition from 2006 to 2010 in Senegal, West Africa
CN107016466A (en) The method and system that a kind of landscape pattern optimizing is built
Zeng et al. Optimal reduction of anthropogenic emissions for air pollution control and the retrieval of emission source from observed pollutants Ӏ. Application of incomplete adjoint operator
TAJIMA Root Phenotyping with Root Modeling: Towards Sustainable Rice Production
Viguria et al. Accuracy of vertical radial plume mapping technique in measuring lagoon gas emissions
Gaffar et al. Prediction of the Topographic Shape of the Ground Surface Using IDW Method through the Rectangular-Neighborhood Approach
Tarasov et al. Topsoil pollution forecasting using artificial neural networks on the example of the abnormally distributed heavy metal at Russian subarctic
Hammerling et al. Improving satellite monitoring of methane emissions: data science is fundamental to better emissions tracking
Viskari et al. Soil carbon estimates by Yasso15 model improved with state data assimilation
Ramos et al. Implementation of a Geographic Information System in the Chemistry Curriculum: An Exercise in Integrating Environmental Analysis and Assessment
Rodríguez Above ground biomass estimation in palm trees using terrestrial LiDAR and tree modelling
Wawrzynczak et al. Bayesian-based approach to application of the genetic algorithm to localize the abrupt atmospheric contamination source
Stettler et al. Terranimo®-a web-based tool for assessment of the risk of soil compaction due to agricultural field traffic.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination