Method for identifying ecological effects of contaminated sites based on groundwater level and vegetation index

Abstract

The application discloses a pollution site ecological influence identification method based on underground water level and vegetation index, relates to the field of environmental science and technology, and comprises the following steps: using a water level analysis method and a vegetation analysis method to analyze an ecological influence analysis area, and obtaining water level coupling parameters and vegetation coupling parameters; establishing an ecological influence identification method; using the ecological influence identification method to identify ecological influence, and obtaining an ecological anomaly direction; the application is used for solving the problem that in the existing pollution site ecological influence identification method, in the aspect of pollution ecology driving mechanism analysis, there is lack of improvement in the coupling analysis of periodic fluctuation of underground water level and vegetation growth trend, so that the ecological influence on the site can only be judged by one-sided factors, and the problem of misjudgment or wrong judgment of ecological influence causes occurs.

Description

Technical Field

[0001] This invention relates to the field of environmental science and technology, specifically to a method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index. Background Technology

[0002] Ecological impact identification of contaminated sites is a process of systematically analyzing and judging the potential ecological impacts of contaminated sites and their surrounding ecosystems. It aims to identify the pollution source, the affected ecological objects, and the potential ecological effects, providing a basis for subsequent risk management and remediation. Contaminated sites typically include industrial sites, mining areas, areas of environmental pollution accidents, and electronic waste treatment plants.

[0003] Existing methods for identifying contaminated sites typically involve detecting and analyzing specific substances within the site. The resulting data is then filtered to obtain comprehensive evaluation indicators. These indicators are combined with expert scoring and weighted assessments to construct a comprehensive evaluation model. This model quantifies the impact of substances within the site, thereby identifying ecological impacts. While this approach can identify the ecological impact caused by a single substance, it lacks improvements in analyzing the driving mechanisms of pollution-related ecological impacts, particularly regarding the coupling analysis of groundwater level fluctuations and vegetation growth trends. This leads to judgments based on only one factor, failing to distinguish between natural climate fluctuations and vegetation degradation caused by pollution itself. This results in misjudgments or incorrect assessments of the causes of ecological impacts, as seen in methods published under CN. Patent application 116794185A discloses a method for identifying and evaluating odor substances in organically contaminated sites. This method accurately identifies odors in organically contaminated sites by precisely quantifying the odor impact in the soil. Other improvements for identifying contaminated sites typically focus on identifying heavy metal factors. However, in terms of analyzing the driving mechanisms of pollution ecology, there is a lack of improvement in the coupling analysis of periodic fluctuations in groundwater levels and vegetation growth trends. This results in the ability to judge the ecological impact of a site based on only one factor, failing to distinguish between natural climate fluctuations and vegetation degradation caused by pollution itself, leading to misjudgments or misdiagnoses of the causes of ecological impacts. Therefore, it is necessary to improve existing methods for identifying the ecological impacts of contaminated sites. Summary of the Invention

[0004] This invention aims to at least partially solve one of the technical problems in the prior art by proposing a method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index. This method addresses the shortcomings of existing methods for identifying the ecological impact of contaminated sites, which lack improvements in the coupling analysis of the periodic fluctuations of groundwater level and vegetation growth trends in identifying the driving mechanisms of pollution ecology. As a result, the ecological impact of a site can only be judged based on one-sided factors, and the method cannot distinguish between natural climate fluctuations and vegetation degradation caused by pollution itself, leading to misjudgments or misdiagnoses of the causes of ecological impact.

[0005] To achieve the above objectives, this application provides a method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index, comprising the following steps: Multiple regions under different ecological impact types were identified and designated as ecological impact analysis regions. Water level analysis and vegetation analysis methods were used to analyze each ecological impact analysis region, and the water level coupling parameters and vegetation coupling parameters were obtained from the analysis results. Based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas, an ecological impact identification method is established. Ecological impact identification methods are used to identify the ecological impact of contaminated sites and to obtain the direction of ecological anomalies based on the results of the ecological impact identification.

[0006] Furthermore, multiple areas under different ecological impact types were identified and denoted as ecological impact analysis areas, including: Multiple areas under different ecological impact types are identified and designated as ecological impact analysis areas; an area with a radius of Gkm centered on the contaminated site to be identified is designated as the identification association area. For any ecological impact analysis area: the ecological impact type of the ecological impact analysis area is recorded as the main type, and the area with the main type of ecological impact in the identified related areas is recorded as the related area of ​​the ecological impact analysis area; Obtain all related regions for the ecological impact analysis area; when no related regions exist for the ecological impact analysis area, record the ecological impact analysis area as an independent analysis area.

[0007] Furthermore, water level analysis methods include: For any ecological impact analysis area, k locations are randomly selected within the ecological impact analysis area and recorded as groundwater sampling points; groundwater levels are continuously acquired at each groundwater sampling point using a water level sensor equipped with a data acquisition instrument; For any groundwater sampling point: establish a plane rectangular coordinate system with the horizontal axis in meters and the vertical axis in h, and denote it as the first water level coordinate system; for any data collected by the water level sensor at the groundwater sampling point, record the groundwater level height detected in the collected data as the groundwater depth A, and record the time when the water level sensor uses the groundwater depth A as the detection result as the water level duration of the groundwater depth A.

[0008] Furthermore, water level analysis methods also include: In the first water level coordinate system, the groundwater depth A is used as the horizontal axis and the water level duration of groundwater depth A is used as the vertical axis to mark the point, and the resulting point is recorded as the water level duration point of groundwater depth A. Based on the continuous detection results obtained by the water level sensor, multiple groundwater depths are obtained at continuous water level points in the first water level coordinate system, and the continuous water level points for each groundwater depth are updated in real time. After the water level sensor is activated for t days, the curve obtained by fitting all water level duration points in the first water level coordinate system is recorded as the sampling point water level curve of the groundwater sampling point, and the water level duration points with the largest and smallest horizontal coordinates are recorded as the right boundary point and left boundary point of the sampling point water level curve.

[0009] Furthermore, water level analysis methods also include: Establish a Cartesian coordinate system with the unit of m for the horizontal axis and the unit of h for the vertical axis, and denote it as the second water level coordinate system. Place the sampling point water level curves of all groundwater sampling points in the second water level coordinate system. The horizontal coordinate of the right boundary point with the largest horizontal coordinate among all sampling point water level curves is marked as Xmax, and the horizontal coordinate of the left boundary point with the smallest horizontal coordinate among all sampling point water level curves is marked as Xmin. The highest and lowest points of the water level curves of all sampling points between [Xmin, Xmax] are recorded as the high persistence point and the low persistence point, respectively; straight lines parallel to the X-axis are drawn at the high persistence point and the low persistence point, respectively, and recorded as the high persistence boundary line and the low persistence boundary line. Let H be the total length of the water level curves of all sampling points between [Xmin, Xmax]. Using a straight line parallel to the X-axis, divide region α into k regions on average, and label them from bottom to top as water level analysis region SF1 to water level analysis region SFk, where region α is the region between the high continuity boundary and the low continuity boundary.

[0010] Furthermore, water level analysis methods also include: For any water level analysis area, the length of the water level curve of all sampling points in the water level analysis area is recorded as h, and h / H is recorded as the water level time ratio of the water level analysis area; the ordinate of the midpoint of the water level analysis area is marked as the water level identification time of the water level analysis area. Obtain the water level identification time and water level time ratio for all water level analysis areas; When the ecological impact analysis area is an independent analysis area, the water level identification time and water level time ratio of all water level analysis areas are recorded as the water level coupling parameters of the ecological impact analysis area. When there are related areas in the ecological impact analysis area, the water level coupling parameters of the ecological impact analysis area are set as follows: the water level identification time and water level time ratio of all water level analysis areas in the ecological impact analysis area, and the water level identification time and water level time ratio of all water level analysis areas in all related areas of the ecological impact analysis area.

[0011] Furthermore, vegetation analysis methods include: For any ecological impact analysis area, the vegetation index within the ecological impact analysis area is obtained using PIE remote sensing software; Based on satellite maps, the areas where all vegetation is located within the ecological impact analysis area are obtained and recorded as ecological vegetation areas; for any groundwater sampling point, the vegetation index at the groundwater sampling point is obtained using PIE remote sensing software and recorded as the sampling vegetation index. The closed interval formed by the maximum and minimum values ​​of the sampling vegetation index of all groundwater sampling points is denoted as the sampling index interval; The vegetation index and the sampling index interval of the ecological impact analysis area are denoted as the vegetation coupling parameters of the ecological impact analysis area.

[0012] Furthermore, methods for identifying ecological impacts include: When identifying the ecological impact of a contaminated site, groundwater sampling points are obtained within the contaminated site based on the water level analysis method. The corresponding water level analysis area of ​​the contaminated site, along with the water level identification time and water level time ratio of all water level analysis areas, are obtained within a second water level coordinate system. An array of all water level identification times from smallest to largest is recorded as the water level time array of the contaminated site, and an array of all water level time ratios is recorded as the water level ratio array of the contaminated site. For any water level time ratio β, the position of water level time ratio β within the water level ratio array is the same as the position of the water level identification time of water level time ratio β within the water level time array. Vegetation parameters and sampling index ranges for contaminated sites were obtained based on vegetation analysis methods.

[0013] Furthermore, ecological impact identification methods also include: For any ecological impact analysis area: When the ecological impact analysis area is an independent analysis area, based on the method of obtaining the water level time array and water level ratio array of the contaminated site, obtain the water level time array and water level ratio array of the ecological impact analysis area; Subtract the water level time array and water level ratio array of the contaminated site from the water level time array and water level ratio array of the ecological impact analysis area respectively, and the resulting arrays are denoted as ecological difference array A and ecological difference array B respectively; The absolute value of the sum of all values ​​in ecological difference array A is denoted as time difference value A, and the absolute value of the sum of all values ​​in ecological difference array B is denoted as water level difference value B; Time difference value A and water level difference value B are denoted as water level determination parameters of the ecological impact analysis area. When there are related areas in the ecological impact analysis area, the time difference value A and water level difference value B of the ecological impact analysis area, as well as the time difference value A and water level difference value B of the related areas of the ecological impact analysis area, are obtained respectively. The minimum value of all time difference values ​​A and the minimum value of all water level difference values ​​B are recorded as the water level determination parameters of the ecological impact analysis area.

[0014] Furthermore, ecological impact identification methods also include: Obtain water level determination parameters for all ecological impact areas, and record the ecological impact areas corresponding to the smallest time difference value A and the smallest water level difference value B among all water level determination parameters as the time-fitting area and the water level-fitting area, respectively. Among all ecological impact areas, the ecological impact area with the smallest difference between the vegetation index and the vegetation parameters of the contaminated site is denoted as the index-fitting area. Among all ecological impact areas, the ecological impact area with the most overlap between the sampling index interval and the sampling index interval of the contaminated site is recorded as the sampling fit area; When there is a sampling matching area among the time matching area, water level matching area, and index matching area, the ecological impact type of the sampling matching area is recorded as the site ecological type. When there is no sampling matching area in the time matching area, water level matching area, and index matching area, the ecological impact type corresponding to the area that appears most frequently in the time matching area, water level matching area, index matching area, and sampling matching area is recorded as the site ecological type. The site's ecological type is recorded as the ecological impact type of the contaminated site, and the direction of ecological anomalies within the site's ecological type is recorded as the ecological anomaly direction of the contaminated site.

[0015] The beneficial effects of this invention are as follows: This application first obtains multiple areas under different ecological impact types and labels them as ecological impact analysis areas; then, it analyzes each ecological impact analysis area using water level analysis and vegetation analysis methods, and obtains water level coupling parameters and vegetation coupling parameters from the analysis results. The advantage of this is that by using water level analysis and vegetation analysis methods to obtain the water level coupling parameters and vegetation coupling parameters corresponding to each ecological impact analysis area, it is possible to obtain the relevant characteristics of groundwater level and vegetation in polluted sites under different ecological impact types. This allows for the establishment of an ecological impact identification method in subsequent analyses, enabling accurate identification of the ecological impact type of the polluted site through the coupling analysis of the periodic fluctuations of groundwater level and vegetation growth trends within the polluted site. This avoids the problem of judging the ecological impact of the site based on only one factor. This application also establishes an ecological impact identification method based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas. Finally, the ecological impact identification method is used to identify the ecological impact of the contaminated site to be identified, and the direction of ecological anomalies is obtained based on the results of the ecological impact identification. The advantage of this is that by establishing an ecological impact identification method, it is possible to ensure that the results of the ecological impact identification of the contaminated site can be correlated with the groundwater level and vegetation growth status within the contaminated site. This allows for the effective identification of natural climate fluctuations and vegetation degradation caused by pollution itself, avoiding misjudgments or misidentifications of the causes of ecological impacts and improving the accuracy of ecological impact identification within the contaminated site. Attached Figure Description

[0016] Figure 1 This is a flowchart illustrating the steps of the method of the present invention; Figure 2 This is a schematic diagram of the second water level coordinate system of the present invention; Figure 3 This is a schematic diagram illustrating the water level-time ratio and the acquisition of water level indicator time according to the present invention. Figure 4 This is a schematic diagram of the electronic device of the present invention. Detailed Implementation

[0017] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0018] Example 1, please refer to Figure 1As shown, this application provides a method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index, including the following steps: Step S1: Obtain multiple regions under different ecological impact types and record them as ecological impact analysis regions; use water level analysis method and vegetation analysis method to analyze each ecological impact analysis region, and obtain water level coupling parameters and vegetation coupling parameters from the analysis results; Step S1 includes: Step S101, obtaining multiple areas under different ecological impact types and recording them as ecological impact analysis areas; and recording the area with a radius of Gkm centered on the polluted site to be identified as the identification association area; In the specific implementation process, the value of G can be determined according to the region where the contaminated site is located and the actual area of ​​the contaminated site. This ensures that the ecological impact analysis area obtained within the identified associated area is related to the ecological impact within the contaminated site, so as to assist in the judgment of the ecological impact type of the contaminated site in subsequent analysis based on the associated area of ​​the ecological impact analysis area. In the data analysis of this embodiment, the value of G is set to 10, that is, the area with a radius of 10km centered on the contaminated site is recorded as the identified associated area and used for the subsequent acquisition of associated areas. Step S102: For any ecological impact analysis area: Record the ecological impact type of the ecological impact analysis area as the main type, and record the area with the main ecological impact type in the identified related areas as the related areas of the ecological impact analysis area. Step S103: Obtain the associatable regions for all ecological impact analysis areas; when there are no associatable regions for an ecological impact analysis area, record the ecological impact analysis area as an independent analysis area.

[0019] Step S104, the water level analysis method includes: Step S1041, for any ecological impact analysis area, randomly select k locations within the ecological impact analysis area and record them as groundwater sampling points; use a water level sensor equipped with a data acquisition instrument to continuously acquire the groundwater level at each groundwater sampling point; In the specific implementation process, the value of k can be set according to the actual area of ​​the contaminated site, the contaminated area within the contaminated site, and the vegetation coverage area within the contaminated site. If the area of ​​the contaminated site is large, the contaminated area within the contaminated site is large, and the vegetation coverage area within the contaminated site is relatively scattered, the value of k can be increased to ensure more detailed data analysis of the groundwater level and vegetation growth status within the contaminated site in subsequent analysis. If the area of ​​the contaminated site is small, the contaminated area within the contaminated site is small, and the vegetation coverage area within the contaminated site is relatively concentrated, the value of k can be decreased to ensure that data related to the ecological impact within the contaminated site is fully acquired while improving the overall data analysis efficiency. In the data analysis of this embodiment, the value of k is set to 5. Step S1042: For any groundwater sampling point: Establish a plane rectangular coordinate system with the horizontal axis in m and the vertical axis in h, and denot it as the first water level coordinate system; For any data collected by the water level sensor at the groundwater sampling point, record the groundwater level height detected in the collected data as groundwater depth A, and record the time when the water level sensor uses groundwater depth A as the detection result as the water level duration of groundwater depth A.

[0020] The water level analysis method also includes: step S1043, in the first water level coordinate system, mark the groundwater depth A as the abscissa and the water level duration of groundwater depth A as the ordinate, and record the obtained mark as the water level duration point of groundwater depth A. Step S1044: Based on the detection results continuously obtained by the water level sensor, obtain multiple groundwater depth water level continuity points in the first water level coordinate system, and update the water level continuity points for each groundwater depth in real time. Step S1045: After the water level sensor is activated for t days, the curve obtained by fitting all water level duration points in the first water level coordinate system is recorded as the sampling point water level curve of the groundwater sampling point, and the water level duration points with the largest and smallest horizontal coordinates are recorded as the right boundary point and left boundary point of the sampling point water level curve. In the specific implementation process, the value of t can be determined according to the actual data analysis capability. In the data analysis of this embodiment, the value of t is set to 7, that is, after the water level sensor is started at the groundwater sampling point for 7 days, the curve obtained by fitting all water level duration points in the first water level coordinate system is recorded as the sampling point water level curve of the groundwater sampling point. For example, in the data analysis of this embodiment, a groundwater depth A is obtained as 6m, and after the water level sensor is started for 7 days, the water level duration of the groundwater depth of 6m is 300h in the detection results obtained by the water level sensor. Then, at this time, a water level duration point in the first water level coordinate system is (6, 300), which means that within 7 days, the groundwater depth of 6m is 300h in the area where the water level sensor is located.

[0021] The water level analysis method also includes: step S1046, establishing a plane rectangular coordinate system with the unit of the horizontal axis being m and the unit of the vertical axis being h, and denoting it as the second water level coordinate system; placing the sampling point water level curves of all groundwater sampling points in the second water level coordinate system, marking the horizontal coordinate of the right boundary point with the largest horizontal coordinate among all sampling point water level curves as Xmax, and marking the horizontal coordinate of the left boundary point with the smallest horizontal coordinate among all sampling point water level curves as Xmin; Step S1047: Record the highest and lowest points of the water level curves of all sampling points between [Xmin, Xmax] as high persistence points and low persistence points, respectively; draw straight lines parallel to the X-axis at the high persistence points and low persistence points, and record them as high persistence boundary lines and low persistence boundary lines, respectively. In the specific implementation process, for example in a data analysis, the second water level coordinate system obtained is as follows: Figure 2 As shown, where, Figure 2 All curves within are the water level curves for all groundwater sampling points; through analysis of... Figure 2 Analyzing the curves within the range, we can obtain Xmax and Xmin as 9.5m and 2m respectively. Therefore, the highest point Top and the lowest point Bot of the water level curves of all sampling points between [2m and 9.5m] can be recorded as the high continuity point and the low continuity point respectively. The straight lines drawn at Top and Bot parallel to the X-axis can be recorded as the high continuity boundary line and the low continuity boundary line respectively. Step S1048: The total length of the water level curves of all sampling points between [Xmin, Xmax] is denoted as H; using a straight line parallel to the X-axis, the region α is divided into k regions on average, and they are denoted from bottom to top as water level analysis region SF1 to water level analysis region SFk, where region α is the region between the high continuity boundary and the low continuity boundary.

[0022] The water level analysis method also includes: step S1049, for any water level analysis area, the length of the water level curve of all sampling points in the water level analysis area is recorded as h, and h / H is recorded as the water level time ratio of the water level analysis area; the ordinate of the midpoint of the water level analysis area is marked as the water level identification time of the water level analysis area; In the specific implementation process, such as in a data analysis, by Figure 2 The region α obtained by the second water level coordinate system in the middle is as follows: Figure 3 As shown, among which, Figure 3The region α is stripped to obtain a water level analysis region called region SF. The length of the water level curve of all sampling points in region SF is 3cm. Through data acquisition, the total length of the water level curve of all sampling points in region α between [2m, 9.5m] is 24cm. Therefore, it can be calculated that the water level time ratio corresponding to region SF is 0.125, and the water level marking time is 43.5. Step S1050: Obtain the water level identification time and water level time ratio for all water level analysis areas; Step S1051: When the ecological impact analysis area is an independent analysis area, the water level identification time and water level time ratio of all water level analysis areas are recorded as the water level coupling parameters of the ecological impact analysis area. Step S1052: When there are related areas in the ecological impact analysis area, the water level coupling parameters of the ecological impact analysis area are set as follows: the water level identification time and water level time ratio of all water level analysis areas in the ecological impact analysis area, and the water level identification time and water level time ratio of all water level analysis areas in all related areas of the ecological impact analysis area.

[0023] Step S106, the vegetation analysis method includes: Step S1061, for any ecological impact analysis area, use PIE remote sensing software to obtain the vegetation index within the ecological impact analysis area; In the specific implementation process, the software that can be used to obtain vegetation indices in actual data analysis can be used to obtain the vegetation indices in the vegetation analysis method. In addition to PIE remote sensing software, ArcGIS / QGIS, ENVI, ESA Sentinel Hub or various remote sensing cloud platforms can also be used to obtain the vegetation index of the area where the vegetation is located in the ecological impact analysis area within the vegetation analysis method. Step S1062: Based on satellite maps, obtain the areas where all vegetation is located within the ecological impact analysis area and record them as ecological vegetation areas; for any groundwater sampling point, use PIE remote sensing software to obtain the vegetation index at the groundwater sampling point and record it as the sampling vegetation index. Step S1063: The closed interval formed by the maximum and minimum values ​​of the sampling vegetation index of all groundwater sampling points is denoted as the sampling index interval. Step S1064: Record the vegetation index and sampling index interval of the ecological impact analysis area as the vegetation coupling parameter of the ecological impact analysis area. In the data analysis of this embodiment, for example, in a single data analysis, the vegetation index and sampling index range of the ecological impact analysis area obtained using PIE remote sensing software are 0.6 and [0.1, 0.9], respectively.

[0024] Step S2: Based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas, establish an ecological impact identification method; The ecological impact identification method includes: Step S201, when conducting ecological impact identification on the contaminated site to be identified, groundwater sampling points are obtained within the contaminated site based on the water level analysis method, and the water level analysis area corresponding to the contaminated site and the water level identification time and water level time ratio of all water level analysis areas are obtained in the second water level coordinate system; the array composed of all water level identification times from smallest to largest is recorded as the water level time array of the contaminated site, and the array composed of all water level time ratios is recorded as the water level ratio array of the contaminated site, wherein, for any water level time ratio β, the position of the water level time ratio β in the water level ratio array is the same as the position of the water level identification time of the water level time ratio β in the water level time array; In specific implementation, for example, in a data analysis, the groundwater level in the contaminated site is analyzed based on the water level analysis method. The resulting water level time array and water level ratio array corresponding to the contaminated site are [30, 90, 150, 210, 270] and [0.1, 0.4, 0.3, 0.1, 0.1], respectively. Step S202: Obtain vegetation parameters and sampling index ranges for the contaminated site based on vegetation analysis methods.

[0025] The ecological impact identification method further includes: step S203, for any ecological impact analysis area: when the ecological impact analysis area is an independent analysis area, based on the acquisition method of the water level time array and water level ratio array of the contaminated site, the water level time array and water level ratio array of the ecological impact analysis area are obtained; the water level time array and water level ratio array of the contaminated site are subtracted from the water level time array and water level ratio array of the ecological impact analysis area respectively, and the resulting arrays are denoted as ecological difference array A and ecological difference array B respectively; the absolute value of the sum of all values ​​in ecological difference array A is denoted as time difference value A, and the absolute value of the sum of all values ​​in ecological difference array B is denoted as water level difference value B; time difference value A and water level difference value B are denoted as water level determination parameters of the ecological impact analysis area. In the data analysis of this embodiment, for example, in a single data analysis, the water level time array and water level ratio array of an ecological impact analysis area with an ecological impact type of "waterlogging and pollution combined stress" are [35, 80, 140, 200, 280] and [0.1, 0.4, 0.3, 0.1, 0.1], respectively. Calculations show that after subtracting the water level time array and water level ratio array of the polluted site from the water level time array and water level ratio array of the ecological impact analysis area, the resulting ecological difference arrays A and B are [5...]. [,-10,-10,-10,-10] and [0,0,0,0,0]; through calculation, the time difference value A is 35 and the water level difference value B is 0; indicating that between the polluted site and the ecological impact analysis area with the ecological impact type of "waterlogging and pollution composite stress", the parameter corresponding to the duration of each level of groundwater level is 35, and the parameter corresponding to the proportion of each level of water level to the overall water level is 0; in the data analysis of this embodiment, the smaller the time difference value A and the water level difference value B, the smaller the difference in groundwater level between the polluted site and the ecological impact analysis area; Step S204: When there are related areas in the ecological impact analysis area, obtain the time difference value A and water level difference value B of the ecological impact analysis area, as well as the time difference value A and water level difference value B of the related areas in the ecological impact analysis area, and record the minimum value of all time difference values ​​A and the minimum value of all water level difference values ​​B as the water level determination parameter of the ecological impact analysis area. In the data analysis of this embodiment, for example, after a data analysis, the time difference value A and the water level difference value B of the ecological impact analysis area are 35 and 0, respectively. There is a correlated area in the ecological impact analysis area, and the time difference value A and the water level difference value B of the correlated area are 20 and 5, respectively. Through analysis, it can be seen that the time difference value A of 20 and the water level difference value B of 0 can be recorded as the water level determination parameters of the ecological impact analysis area.

[0026] The ecological impact identification method also includes: step S205, obtaining water level determination parameters for all ecological impact areas, and recording the ecological impact areas corresponding to the smallest time difference value A and the smallest water level difference value B among all water level determination parameters as time-matching areas and water level-matching areas, respectively. Step S206: Among all ecological impact areas, the ecological impact area with the smallest difference between the vegetation index and the vegetation parameter of the contaminated site is recorded as the index-matching area. Step S207: Among all ecological impact areas, the ecological impact area with the most overlap between the sampling index interval and the sampling index interval of the contaminated site is recorded as the sampling fit area. Step S208: When there is a sampling matching area in the time matching area, water level matching area and index matching area, the ecological impact type of the sampling matching area is recorded as the site ecological type. In the data analysis of this embodiment, for example, in a single data analysis, the ecological impact type of the sampling matching area is "waterlogging pollution complex stress". The ecological impact types corresponding to the time matching area, water level matching area and index matching area are "groundwater pollution diffusion stress", "waterlogging pollution complex stress" and "waterlogging pollution complex stress", respectively. Therefore, "waterlogging pollution complex stress" can be regarded as the ecological impact type of the contaminated site. Step S209: When there is no sampling matching area in the time matching area, water level matching area, and index matching area, the ecological impact type corresponding to the area that appears most frequently in the time matching area, water level matching area, index matching area, and sampling matching area is recorded as the site ecological type. In the data analysis of this embodiment, since the sampling fit area is the area that best matches the contaminated site by combining the vegetation distribution in various places within the contaminated site, the sampling fit area can be used as a decisive factor to screen the time fit area, water level fit area, and index fit area. If there is an area in the time fit area, water level fit area, and index fit area that is the same as the sampling fit area, the ecological impact type corresponding to that area can be recorded as the ecological impact type of the contaminated site. If there is no area in the time fit area, water level fit area, and index fit area that is the same as the sampling fit area, the ecological impact type of the contaminated site can be determined by combining the area that appears most frequently in the time fit area, water level fit area, index fit area, and sampling fit area. Alternatively, the time fit area, water level fit area, index fit area, and sampling fit area can be reported to assist staff in making more accurate screenings. Step S210: Record the site ecological type as the ecological impact type of the polluted site, and record the direction of ecological anomaly in the site ecological type as the ecological anomaly direction of the polluted site.

[0027] Step S3: Use the ecological impact identification method to identify the ecological impact of the contaminated site to be identified, and obtain the direction of ecological anomalies based on the results of the ecological impact identification.

[0028] Example 2, please refer to Figure 4 As shown, Figure 4A schematic diagram of an electronic device is provided, which may include a processor, a communication interface, a memory, and a communication bus. The processor, communication interface, and memory communicate with each other via the communication bus. The memory stores computer-readable instructions, which the processor can call. When the processor executes a computer-readable instruction, it performs steps such as those in a method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index, to achieve the following functions: First, multiple areas under different ecological impact types are identified and designated as ecological impact analysis areas; each ecological impact analysis area is analyzed using both groundwater level analysis and vegetation analysis methods, and the analysis results yield groundwater level coupling parameters and vegetation coupling parameters; then, based on the ecological impact types, groundwater level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas, an ecological impact identification method is established; finally, the ecological impact identification method is used to identify the ecological impact of the contaminated site to be identified, and the direction of ecological anomalies is obtained based on the results of the ecological impact identification.

[0029] Furthermore, when the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0030] Example 3: This application also provides a computer program product, which includes a computer program stored on a computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the above-described methods for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index. This method includes: first, acquiring multiple areas under different ecological impact types and labeling them as ecological impact analysis areas; analyzing each ecological impact analysis area using both water level analysis and vegetation analysis methods, and obtaining water level coupling parameters and vegetation coupling parameters from the analysis results; then, establishing an ecological impact identification method based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas; and finally, using the ecological impact identification method to identify the ecological impact of the contaminated site to be identified, and obtaining the direction of ecological anomalies based on the results of the ecological impact identification.

[0031] Example 4: This application also provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, it performs the steps of the above-mentioned method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index, to achieve the following functions: First, multiple areas under different ecological impact types are acquired and denoted as ecological impact analysis areas; each ecological impact analysis area is analyzed using water level analysis and vegetation analysis methods, and the water level coupling parameters and vegetation coupling parameters are obtained from the analysis results; then, based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas, an ecological impact identification method is established; finally, the ecological impact identification method is used to identify the ecological impact of the contaminated site to be identified, and the direction of ecological anomalies is obtained based on the results of the ecological impact identification.

[0032] Based on the above description of the embodiments, the embodiments of the present invention can be provided as methods, systems, or computer program products. Based on this understanding, the above technical solutions, in essence or in terms of their contribution to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or certain parts of the embodiments.

[0033] In the embodiments provided in this application, it should be understood that the disclosed system or method can be implemented in other ways. The embodiments described above are merely illustrative. For example, the division of modules or units is only a logical functional division, and there may be other division methods in actual implementation. Furthermore, multiple modules or units may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces. The indirect coupling or communication connection between systems, modules, and units may be electrical, mechanical, or other forms.

[0034] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.

Claims

1. A method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index, characterized in that, Includes the following steps: Multiple regions under different ecological impact types were identified and designated as ecological impact analysis regions. Water level analysis and vegetation analysis methods were used to analyze each ecological impact analysis region, and the water level coupling parameters and vegetation coupling parameters were obtained from the analysis results. Based on the ecological impact types, water level coupling parameters, and vegetation coupling parameters of all ecological impact analysis areas, an ecological impact identification method is established. Ecological impact identification methods are used to identify the ecological impact of contaminated sites and to obtain the direction of ecological anomalies based on the results of the ecological impact identification.

2. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 1, characterized in that, Multiple regions under different ecological impact types were identified and denoted as ecological impact analysis areas, including: Multiple areas under different ecological impact types are identified and designated as ecological impact analysis areas; an area with a radius of Gkm centered on the contaminated site to be identified is designated as the identification association area. For any ecological impact analysis area: the ecological impact type of the ecological impact analysis area is recorded as the main type, and the area with the main type of ecological impact in the identified related areas is recorded as the related area of ​​the ecological impact analysis area; Obtain all related regions for the ecological impact analysis area; when no related regions exist for the ecological impact analysis area, the ecological impact analysis area is recorded as an independent analysis area.

3. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 2, characterized in that, Water level analysis methods include: For any ecological impact analysis area, k locations are randomly selected within the ecological impact analysis area and recorded as groundwater sampling points; groundwater levels are continuously acquired at each groundwater sampling point using a water level sensor equipped with a data acquisition instrument; For any groundwater sampling point: establish a plane rectangular coordinate system with the horizontal axis in meters and the vertical axis in h, and denote it as the first water level coordinate system; for any data collected by the water level sensor at the groundwater sampling point, record the groundwater level height detected in the collected data as the groundwater depth A, and record the time when the water level sensor uses the groundwater depth A as the detection result as the water level duration of the groundwater depth A.

4. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 3, characterized in that, Water level analysis methods also include: In the first water level coordinate system, the groundwater depth A is used as the horizontal axis and the water level duration of groundwater depth A is used as the vertical axis to mark the point, and the resulting point is recorded as the water level duration point of groundwater depth A. Based on the continuous detection results obtained by the water level sensor, multiple groundwater depths are obtained at continuous water level points in the first water level coordinate system, and the continuous water level points for each groundwater depth are updated in real time. After the water level sensor is activated for t days, the curve obtained by fitting all water level duration points in the first water level coordinate system is recorded as the sampling point water level curve of the groundwater sampling point, and the water level duration points with the largest and smallest horizontal coordinates are recorded as the right boundary point and left boundary point of the sampling point water level curve.

5. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 4, characterized in that, Water level analysis methods also include: Establish a Cartesian coordinate system with the unit of m for the horizontal axis and the unit of h for the vertical axis, and denote it as the second water level coordinate system. Place the sampling point water level curves of all groundwater sampling points in the second water level coordinate system. The horizontal coordinate of the right boundary point with the largest horizontal coordinate among all sampling point water level curves is marked as Xmax, and the horizontal coordinate of the left boundary point with the smallest horizontal coordinate among all sampling point water level curves is marked as Xmin. The highest and lowest points of the water level curves of all sampling points between [Xmin, Xmax] are recorded as the high persistence point and the low persistence point, respectively; straight lines parallel to the X-axis are drawn at the high persistence point and the low persistence point, respectively, and recorded as the high persistence boundary line and the low persistence boundary line. Let H be the total length of the water level curves of all sampling points between [Xmin, Xmax]. Using a straight line parallel to the X-axis, divide region α into k regions on average, and label them from bottom to top as water level analysis region SF1 to water level analysis region SFk, where region α is the region between the high continuity boundary and the low continuity boundary.

6. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 5, characterized in that, Water level analysis methods also include: For any water level analysis area, the length of the water level curve of all sampling points in the water level analysis area is recorded as h, and h / H is recorded as the water level time ratio of the water level analysis area; the ordinate of the midpoint of the water level analysis area is marked as the water level identification time of the water level analysis area. Obtain the water level identification time and water level time ratio for all water level analysis areas; When the ecological impact analysis area is an independent analysis area, the water level identification time and water level time ratio of all water level analysis areas are recorded as the water level coupling parameters of the ecological impact analysis area. When there are related areas in the ecological impact analysis area, the water level coupling parameters of the ecological impact analysis area are set as follows: the water level identification time and water level time ratio of all water level analysis areas in the ecological impact analysis area, and the water level identification time and water level time ratio of all water level analysis areas in all related areas of the ecological impact analysis area.

7. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 6, characterized in that, Vegetation analysis methods include: For any ecological impact analysis area, the vegetation index within the ecological impact analysis area is obtained using PIE remote sensing software; Based on satellite maps, the areas where all vegetation is located within the ecological impact analysis area are obtained and recorded as ecological vegetation areas; for any groundwater sampling point, the vegetation index at the groundwater sampling point is obtained using PIE remote sensing software and recorded as the sampling vegetation index. The closed interval formed by the maximum and minimum values ​​of the sampling vegetation index of all groundwater sampling points is denoted as the sampling index interval; The vegetation index and the sampling index interval of the ecological impact analysis area are denoted as the vegetation coupling parameters of the ecological impact analysis area.

8. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 7, characterized in that, Ecological impact identification methods include: When identifying the ecological impact of a contaminated site, groundwater sampling points are obtained within the contaminated site based on the water level analysis method. The corresponding water level analysis area of ​​the contaminated site, along with the water level identification time and water level time ratio of all water level analysis areas, are obtained within a second water level coordinate system. An array of all water level identification times from smallest to largest is recorded as the water level time array of the contaminated site, and an array of all water level time ratios is recorded as the water level ratio array of the contaminated site. For any water level time ratio β, the position of water level time ratio β within the water level ratio array is the same as the position of the water level identification time of water level time ratio β within the water level time array. Vegetation parameters and sampling index ranges for contaminated sites were obtained based on vegetation analysis methods.

9. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 8, characterized in that, Ecological impact identification methods also include: For any ecological impact analysis area: When the ecological impact analysis area is an independent analysis area, based on the method of obtaining the water level time array and water level ratio array of the contaminated site, obtain the water level time array and water level ratio array of the ecological impact analysis area; Subtract the water level time array and water level ratio array of the contaminated site from the water level time array and water level ratio array of the ecological impact analysis area respectively, and the resulting arrays are denoted as ecological difference array A and ecological difference array B respectively; The absolute value of the sum of all values ​​in ecological difference array A is denoted as time difference value A, and the absolute value of the sum of all values ​​in ecological difference array B is denoted as water level difference value B; Time difference value A and water level difference value B are denoted as water level determination parameters of the ecological impact analysis area. When there are related areas in the ecological impact analysis area, the time difference value A and water level difference value B of the ecological impact analysis area, as well as the time difference value A and water level difference value B of the related areas of the ecological impact analysis area, are obtained respectively. The minimum value of all time difference values ​​A and the minimum value of all water level difference values ​​B are recorded as the water level determination parameters of the ecological impact analysis area.

10. The method for identifying the ecological impact of contaminated sites based on groundwater level and vegetation index according to claim 9, characterized in that, Ecological impact identification methods also include: Obtain water level determination parameters for all ecological impact areas, and record the ecological impact areas corresponding to the smallest time difference value A and the smallest water level difference value B among all water level determination parameters as the time-fitting area and the water level-fitting area, respectively. Among all ecological impact areas, the ecological impact area with the smallest difference between the vegetation index and the vegetation parameters of the contaminated site is denoted as the index-fitting area. Among all ecological impact areas, the ecological impact area with the most overlap between the sampling index interval and the sampling index interval of the contaminated site is recorded as the sampling fit area; When there is a sampling matching area among the time matching area, water level matching area, and index matching area, the ecological impact type of the sampling matching area is recorded as the site ecological type. When there is no sampling matching area in the time matching area, water level matching area, and index matching area, the ecological impact type corresponding to the area that appears most frequently in the time matching area, water level matching area, index matching area, and sampling matching area is recorded as the site ecological type. The site's ecological type is recorded as the ecological impact type of the contaminated site, and the direction of ecological anomalies within the site's ecological type is recorded as the ecological anomaly direction of the contaminated site.

Images