Method and system for analyzing and diagnosing pollution intensity of riverway and lake in flood season

Abstract

The invention relates to the technical field of computers, and discloses a river and lake flood season pollution intensity analysis and diagnosis method and system, and the method comprises the steps: carrying out the gridding of sampling points, and obtaining the pollution data parameters of each sampling point; acquiring a corresponding pollution intensity estimation value according to the pollution data parameter; judging whether the pollution intensity estimation value is abnormal or not, and if yes, obtaining a pollution source estimation range according to the pollution intensity estimation value and the position corresponding to the sampling point; acquiring coastal surface soil data in the pollution source estimation range; and determining an accurate pollution source according to the matching degree of the coastal surface soil data and the pollution data parameters, thereby greatly improving the efficiency and accuracy of determining the pollution source and reducing the workload of personnel.

Description

Technical Field

[0001] The present invention relates to the field of computer technology, and in particular to a method and system for analyzing and diagnosing pollution intensity in rivers and lakes during flood season. Background Art

[0002] Pollution intensity during flood season refers to the ratio of the concentration of primary pollutants in a section during flood season to the target concentration limit of assessment. The flood season covers the period of the year when floods occur regularly in rivers and the period when the river water rises and falls to a certain water level.

[0003] The pollution intensity during the flood season mainly reflects the gap between the pollution level during the flood season in the monitoring section and the water quality target. It is conducive to accurately identifying the situation where the water environment quality is good in normal times and the pollutant concentration increases sharply during the flood season, and to clarify the responsibility for non-point source pollution prevention and control in relevant administrative regions.

[0004] For the cross section, when there is no external pollution, the water quality data is basically stable, but the surface runoff formed by precipitation will introduce pollution sources on land into the river, which will cause the water quality index data to increase. Therefore, how to accurately detect pollution, accurately locate the source of pollution, and deal with it in time is the main problem we are facing. Summary of the invention

[0005] The purpose of the present invention is to provide a method and system for analyzing and diagnosing pollution intensity in rivers and lakes during flood season, so as to solve the following technical problems:

[0006] How to provide a diagnostic method and system that can accurately detect pollution and determine the location of pollution sources through flood season pollution intensity analysis.

[0007] The purpose of the present invention can be achieved through the following technical solutions:

[0008] The analysis and diagnosis methods of pollution intensity of rivers and lakes during flood season include:

[0009] Grid the sampling points and obtain the pollution data parameters of each sampling point;

[0010] Obtaining a corresponding pollution intensity estimation value according to the pollution data parameter;

[0011] Determine whether the pollution intensity estimate is abnormal, and if so, obtain the estimated range of the pollution source according to the pollution intensity estimate and the position of the corresponding sampling point;

[0012] Obtaining coastal surface soil data within the estimated range of the pollution source;

[0013] The exact pollution source is confirmed based on the degree of matching between the coastal surface soil data and the pollution data parameters.

[0014] As a further solution of the present invention: the method of obtaining the pollution data parameters of each sampling point by gridding the sampling points includes:

[0015] Evenly arrange multiple sampling modules in the target waters;

[0016] Wherein, the sampling module comprises a floating sampling unit and a plurality of underwater sampling units connected in sequence, and the floating sampling unit is connected to one of the underwater sampling units at the end;

[0017] Fixing the floating sampling unit at the edge position;

[0018] For one of the sampling modules, the flow velocity and flow direction of different water faults are obtained through the relative position change data between the floating sampling unit and the corresponding underwater sampling unit and the relative position change data between adjacent underwater sampling units. Each of the sampling modules is provided with one underwater sampling unit in a water fault.

[0019] The underwater sampling unit is also used to obtain the pollution data parameters.

[0020] As a further solution of the present invention: the method for obtaining the corresponding pollution intensity estimation value according to the pollution data parameter includes:

[0021] According to the pollution data parameters collected by the underwater sampling unit, the corresponding pollution index data W is obtained in real time i (t);

[0022]

[0023] in, is the preset weight coefficient corresponding to the concentration of the qth water quality index, WQPI q (t) represents the pollution index corresponding to the concentration of the qth water quality indicator detected by the underwater sampling unit at the i coordinate at time t;

[0024] According to the pollution index data W i (t) obtaining the pollution intensity estimation value F(t);

[0025]

[0026] Among them, n i is the serial number of the underwater sampling unit at the i coordinate, k is the total number of the underwater sampling units, is the position weight corresponding to the underwater sampling unit at the i coordinate.

[0027] As a further solution of the present invention: the method for determining whether the pollution intensity estimate is abnormal includes:

[0028] If F(t)≥F sh (t) or H ≥ δk, it is considered abnormal, otherwise it is normal;

[0029] Among them, F sh (t) is the preset strength threshold, H is the value of W i (t)≥W ith (t) is the number of underwater sampling units, δ is a preset coefficient, 0.3<δ≤0.6.

[0030] As a further solution of the present invention: the method for obtaining the estimated range of the pollution source according to the estimated value of the pollution intensity and the position of the corresponding sampling point includes:

[0031] Obtain the first position change i1(t) of all the underwater sampling units in the designated water fault where the average pollution index data changes the most within the period Δt;

[0032] Get the number of W that occurs within the Δt period i (t)≥W ith (t) a second position change i2(t) of the underwater sampling unit;

[0033] Connect the locations of all the underwater sampling units in the designated water fault to form the first pollution response area, and W i (t)≥W ith (t) The positions of the outermost underwater sampling units are connected to form a second pollution response area;

[0034] Obtaining a first comprehensive trend direction according to the first position change i1(t), extending the first pollution response area in the first comprehensive trend direction, and taking the intersection of the extended area and the water bank as the first pollution source area;

[0035] Obtaining a second comprehensive trend direction according to the second position change i2(t), extending the second pollution response area in the second comprehensive trend direction, and taking the intersection of the extension and the water bank as the second pollution source area;

[0036] The combined area of ​​the first pollution source area and the second pollution source area is used as the pollution source estimation range.

[0037] As a further solution of the present invention: the method for obtaining the first comprehensive trend direction includes:

[0038] Draw a pollution intensity gradient map according to the pollution intensity estimation value of each underwater sampling unit in the designated water fault;

[0039] For the underwater sampling unit within the highest pollution intensity range in the pollution intensity gradient map, obtaining the water flow direction when the positive change rate of the pollution intensity estimation value corresponding to the underwater sampling unit is the largest;

[0040] The water flow direction corresponding to each underwater sampling unit in the designated water fault is superimposed with weight vectors, and the corresponding weight is positively correlated with the magnitude of the corresponding positive change rate to obtain the first comprehensive trend direction.

[0041] The pollution intensity analysis and diagnosis system for rivers and lakes during flood season includes:

[0042] The sampling module is arranged in a grid pattern and is used to obtain pollution data parameters at each sampling point;

[0043] A pollution intensity assessment module, used for obtaining a corresponding pollution intensity estimation value according to the pollution data parameter;

[0044] A pollution range determination module, used to determine whether the pollution intensity estimation value is abnormal, and if so, to obtain the pollution source estimation range according to the pollution intensity estimation value and the position of the corresponding sampling point;

[0045] A field data acquisition module, used to acquire coastal surface soil data within the estimated range of the pollution source;

[0046] The association management module is used to confirm the accurate pollution source according to the matching degree between the coastal surface soil data and the pollution data parameters.

[0047] The beneficial effects of the present invention are as follows: the present invention can comprehensively obtain pollution data parameters at different horizontal positions and different vertical positions in the water area through gridded sampling points, and then obtain corresponding pollution intensity estimation values ​​to determine whether there are abnormalities; once an abnormality is found, the estimated range of the pollution source can be obtained according to the pre-gridded sampling points to determine the approximate location and area of ​​the pollution source, and then the coastal surface soil data within the estimated range of the pollution source can be obtained online to narrow the evaluation range and determine the specific location of the pollution source, thereby greatly improving the efficiency and accuracy of pollution source determination and reducing the workload of personnel. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The present invention will be further described below in conjunction with the accompanying drawings.

[0049] Figure 1 It is a schematic diagram of the method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season in the present invention. DETAILED DESCRIPTION

[0050] 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 only 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.

[0051] See also Figure 1 As shown, the present invention is a method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season, comprising:

[0052] Grid the sampling points and obtain the pollution data parameters of each sampling point;

[0053] Obtaining a corresponding pollution intensity estimation value according to the pollution data parameter;

[0054] Determine whether the pollution intensity estimate is abnormal, and if so, obtain the estimated range of the pollution source according to the pollution intensity estimate and the position of the corresponding sampling point;

[0055] Obtaining coastal surface soil data within the estimated range of the pollution source;

[0056] The exact pollution source is confirmed based on the degree of matching between the coastal surface soil data and the pollution data parameters.

[0057] In this embodiment of the present invention, the pollution data parameters at different horizontal positions and different vertical positions in the water area can be comprehensively obtained through gridded sampling points, and then the corresponding pollution intensity estimation value can be obtained to determine whether there is an abnormality; once an abnormality is found, the estimated range of the pollution source can be obtained according to the pre-gridded sampling points to determine the approximate location and area of ​​the pollution source, and then the coastal surface soil data within the estimated range of the pollution source can be obtained online to narrow the evaluation range and determine the specific location of the pollution source, thereby greatly improving the efficiency and accuracy of pollution source determination and reducing personnel workload.

[0058] As a further solution of the present invention: the method of obtaining the pollution data parameters of each sampling point by gridding the sampling points includes:

[0059] Evenly arrange multiple sampling modules in the target waters;

[0060] Wherein, the sampling module comprises a floating sampling unit and a plurality of underwater sampling units connected in sequence, and the floating sampling unit is connected to one of the underwater sampling units at the end;

[0061] Fixing the floating sampling unit at the edge position;

[0062] For one of the sampling modules, the flow velocity and flow direction of different water faults are obtained through the relative position change data between the floating sampling unit and the corresponding underwater sampling unit and the relative position change data between adjacent underwater sampling units. One underwater sampling unit is correspondingly arranged in one water fault for each sampling module; the underwater sampling unit is also used to obtain the pollution data parameters.

[0063] The floating sampling unit and the corresponding underwater sampling unit are equipped with a GPS positioning unit. When performing the initial grid layout, an initial position is set for each floating sampling unit and the corresponding underwater sampling unit, so the above-mentioned relative position change data can be obtained based on the real-time position information.

[0064] As a further solution of the present invention: the method for obtaining the corresponding pollution intensity estimation value according to the pollution data parameter includes:

[0065] According to the pollution data parameters collected by the underwater sampling unit, the corresponding pollution index data W is obtained in real time i (t);

[0066]

[0067] in, is the preset weight coefficient corresponding to the concentration of the qth water quality index, WQPI q (t) represents the pollution index corresponding to the qth water quality index concentration detected by the underwater sampling unit at the i coordinate at time t; the type item corresponding to the water quality index concentration varies according to the actual situation and can be reasonably set according to the specific situation. In the present invention, at least m water quality index concentrations can be set, such as ammonia nitrogen concentration, total phosphorus concentration, COD concentration, and permanganate index concentration.

[0068] According to the pollution index data W i (t) obtaining the pollution intensity estimation value F(t);

[0069]

[0070] Among them, n i is the serial number of the underwater sampling unit at the i coordinate, k is the total number of the underwater sampling units, is the position weight corresponding to the underwater sampling unit at the i coordinate.

[0071] As a further solution of the present invention: the method for determining whether the pollution intensity estimate is abnormal includes:

[0072] If F(t)≥F sh (t) or H ≥ δk, it is considered abnormal, otherwise it is normal;

[0073] Among them, F sh (t) is the preset strength threshold, H is the value of W i (t)≥W ith (t) is the number of underwater sampling units, δ is a preset coefficient, 0.3<δ≤0.6.

[0074] As a further solution of the present invention: the method for obtaining the estimated range of the pollution source according to the estimated value of the pollution intensity and the position of the corresponding sampling point includes:

[0075] Obtain the first position change i1(t) of all the underwater sampling units in the designated water fault where the average pollution index data changes the most within the period Δt;

[0076] Get the number of W that occurs within the Δt period i (t)≥W ith (t) a second position change i2(t) of the underwater sampling unit;

[0077] Connect the locations of all the underwater sampling units in the designated water fault to form the first pollution response area, and W i (t)≥W ith (t) The positions of the outermost underwater sampling units are connected to form a second pollution response area;

[0078] Obtaining a first comprehensive trend direction according to the first position change i1(t), extending the first pollution response area in the first comprehensive trend direction, and taking the intersection of the extended area and the water bank as the first pollution source area;

[0079] Obtaining a second comprehensive trend direction according to the second position change i2(t), extending the second pollution response area in the second comprehensive trend direction, and taking the intersection of the extension and the water bank as the second pollution source area;

[0080] The combined area of ​​the first pollution source area and the second pollution source area is used as the pollution source estimation range.

[0081] As a further solution of the present invention: the method for obtaining the first comprehensive trend direction includes:

[0082] Draw a pollution intensity gradient map according to the pollution intensity estimation value of each underwater sampling unit in the designated water fault;

[0083] For the underwater sampling unit within the highest pollution intensity range in the pollution intensity gradient map, obtaining the water flow direction when the positive change rate of the pollution intensity estimation value corresponding to the underwater sampling unit is the largest;

[0084] The water flow direction corresponding to each underwater sampling unit in the designated water fault is superimposed with weight vectors, and the corresponding weight is positively correlated with the magnitude of the corresponding positive change rate to obtain the first comprehensive trend direction.

[0085] The pollution intensity analysis and diagnosis system for rivers and lakes during flood season includes:

[0086] The sampling module is arranged in a grid pattern and is used to obtain pollution data parameters at each sampling point;

[0087] A pollution intensity assessment module, used for obtaining a corresponding pollution intensity estimation value according to the pollution data parameter;

[0088] A pollution range determination module, used to determine whether the pollution intensity estimation value is abnormal, and if so, to obtain the pollution source estimation range according to the pollution intensity estimation value and the position of the corresponding sampling point;

[0089] A field data acquisition module, used to acquire coastal surface soil data within the estimated range of the pollution source;

[0090] The association management module is used to confirm the accurate pollution source according to the matching degree between the coastal surface soil data and the pollution data parameters.

[0091] The above is a detailed description of an embodiment of the present invention, but the content is only a preferred embodiment of the present invention and cannot be considered to limit the scope of implementation of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the patent coverage of the present invention.

Claims

1. A method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season, characterized in that: include: Grid the sampling points and obtain the pollution data parameters of each sampling point; Obtaining a corresponding pollution intensity estimation value according to the pollution data parameter; Determine whether the pollution intensity estimate is abnormal, and if so, obtain the estimated range of the pollution source according to the pollution intensity estimate and the position of the corresponding sampling point; Obtaining coastal surface soil data within the estimated range of the pollution source; The exact pollution source is confirmed based on the degree of matching between the coastal surface soil data and the pollution data parameters.

2. The method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season according to claim 1 is characterized in that: The method of gridding sampling points and obtaining pollution data parameters of each sampling point includes: Arrange multiple sampling modules evenly in the target waters; Wherein, the sampling module comprises a floating sampling unit and a plurality of underwater sampling units connected in sequence, and the floating sampling unit is connected to one of the underwater sampling units at the end; Fixing the floating sampling unit at the edge position; For one of the sampling modules, the flow velocity and flow direction of different water faults are obtained through the relative position change data between the floating sampling unit and the corresponding underwater sampling unit and the relative position change data between adjacent underwater sampling units. Each of the sampling modules is provided with one underwater sampling unit in a water fault. The underwater sampling unit is also used to obtain the pollution data parameters.

3. The method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season according to claim 2 is characterized in that: The method for obtaining the corresponding pollution intensity estimation value according to the pollution data parameter includes: According to the pollution data parameters collected by the underwater sampling unit, the corresponding pollution index data W is obtained in real time i (t); in, is the preset weight coefficient corresponding to the concentration of the qth water quality index, WQPI q (t) represents the pollution index corresponding to the concentration of the qth water quality indicator detected by the underwater sampling unit at the i coordinate at time t; According to the pollution index data W i (t) obtaining the pollution intensity estimation value F(t); Among them, n i is the serial number of the underwater sampling unit at the i coordinate, k is the total number of the underwater sampling units, is the position weight corresponding to the underwater sampling unit at the i coordinate.

4. The method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season according to claim 3 is characterized in that: The method for determining whether the pollution intensity estimation value is abnormal includes: If F(t)≥F sh (t) or H ≥ δk, it is considered abnormal, otherwise it is normal; Among them, F sh (t) is the preset strength threshold, H is the value of W i (t)≥W ith (t) is the number of underwater sampling units, δ is a preset coefficient, 0.3<δ≤0.

6.

5. The method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season according to claim 4 is characterized in that: The method for obtaining the estimated range of the pollution source according to the estimated value of the pollution intensity and the position of the corresponding sampling point includes: Obtain the first position change i1(t) of all the underwater sampling units in the designated water fault where the average pollution index data changes the most within the period Δt; Get the number of W that occurs within the Δt period i (t)≥W ith (t) a second position change i2(t) of the underwater sampling unit; Connect the locations of all the underwater sampling units in the designated water fault to form the first pollution response area, and W i (t)≥W ith (t) The positions of the outermost underwater sampling units are connected to form a second pollution response area; Obtaining a first comprehensive trend direction according to the first position change i1(t), extending the first pollution response area in the first comprehensive trend direction, and taking the intersection of the extended area and the water bank as the first pollution source area; Obtaining a second comprehensive trend direction according to the second position change i2(t), extending the second pollution response area in the second comprehensive trend direction, and taking the intersection of the extension and the water bank as the second pollution source area; The combined area of ​​the first pollution source area and the second pollution source area is used as the pollution source estimation range.

6. The method for analyzing and diagnosing pollution intensity of rivers and lakes during flood season according to claim 5 is characterized in that: The method for obtaining the first comprehensive trend direction includes: Draw a pollution intensity gradient map according to the pollution intensity estimation value of each underwater sampling unit in the designated water fault; For the underwater sampling unit within the highest pollution intensity range in the pollution intensity gradient map, obtaining the water flow direction when the positive change rate of the pollution intensity estimation value corresponding to the underwater sampling unit is the largest; The water flow direction corresponding to each underwater sampling unit in the designated water fault is superimposed with weight vectors, and the corresponding weight is positively correlated with the magnitude of the corresponding positive change rate to obtain the first comprehensive trend direction.

7. A system for analyzing and diagnosing pollution intensity during flood season in rivers and lakes, which adopts a method for analyzing and diagnosing pollution intensity during flood season in rivers and lakes as claimed in any one of claims 1 to 6, characterized in that: include: The sampling module is arranged in a grid pattern and is used to obtain pollution data parameters at each sampling point; A pollution intensity assessment module, used for obtaining a corresponding pollution intensity estimation value according to the pollution data parameter; A pollution range determination module, used to determine whether the pollution intensity estimation value is abnormal, and if so, to obtain the pollution source estimation range according to the pollution intensity estimation value and the position of the corresponding sampling point; A field data acquisition module, used to acquire coastal surface soil data within the estimated range of the pollution source; The association management module is used to confirm the accurate pollution source according to the matching degree between the coastal surface soil data and the pollution data parameters.

Images