A method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions

Abstract

This invention relates to the field of soil remediation technology, specifically to a method for dynamically formulating soil remediation targets based on protecting the groundwater conservation function. The method includes the following steps: acquiring remediation starting point and water level data; identifying lag periods; adjusting boundaries to generate dynamic maps; extracting water level gradients and seepage directions to prioritize interventions; determining fluctuations and water retention capacity to mark control zones; and overlaying concentration and water level trends to set a remediation timetable. In this invention, the lag relationship between interventions is clarified by extracting water level inflection points and infiltration differences. The boundary response direction is adjusted based on the lag period to enhance the adaptability of the control range. An intervention priority order is constructed by combining hydraulic gradients and seepage trends to highlight the remediation efficiency of continuous segments. Restricted areas are determined using water layer fluctuations and water retention capacity to accurately delineate concentration control locations. The remediation frequency and timetable are set by overlaying water level trends and pollution distribution, achieving dynamic coordination of intervention identification, boundary adjustment, and remediation rhythm.

Description

Technical Field

[0001] This invention relates to the field of soil remediation technology, and in particular to a method for dynamically formulating soil remediation targets based on protecting the function of groundwater conservation. Background Technology

[0002] Soil remediation technology involves the treatment and restoration of contaminated soil to reduce or eliminate the risks posed by harmful substances to the ecological environment and human health, and is a key component of environmental engineering. Core aspects of this technology include the identification and assessment of contaminated soil, the setting of remediation targets, the selection and implementation of remediation technologies, and post-remediation monitoring and management. It typically encompasses multiple methods such as physical remediation, chemical remediation, and bioremediation, and remediation plans are comprehensively formulated based on the type of pollutant, site characteristics, and environmental functional requirements. Particularly in the protection of regional water and soil ecosystems, the synergy between soil remediation and groundwater conservation functions is crucial and has become an important direction in soil environmental governance in recent years. Among them, the traditional method of dynamically setting soil remediation targets refers to setting remediation targets based on the type and degree of contaminated soil and dynamically adjusting the targets during implementation to adapt to changing environmental conditions and remediation progress. The technical issue it addresses is how to scientifically set and dynamically adjust soil remediation targets to effectively protect the groundwater conservation function. Traditionally, it usually adopts the method of regularly revising remediation targets based on pollutant concentration detection data, hydrogeological survey results and national or local soil environmental quality standards, combined with pollutant migration trend simulation results and risk assessment parameters. Common methods include site zoning strategies based on regional groundwater recharge intensity, pollution source response analysis under hydrological boundary conditions, and setting remediation paths and threshold standards by referring to empirical model parameters of areas with similar geological conditions.

[0003] Existing technologies rely on static data structures of phased detection results and geological survey parameters during the adjustment of remediation targets. They fail to build a continuous response discrimination system. In scenarios with frequent dynamic changes in groundwater or unclear timing of human intervention, it is difficult to determine in a timely manner whether the intervention has produced an effective response. This results in problems such as lag in target adjustment and solidification of control interfaces. For example, when hydrological boundaries change, the original remediation area lacks dynamic correction capabilities, causing misalignment between the timing of remediation intervention and the pollution diffusion path, which affects the accurate allocation of remediation resources and the coordinated maintenance of regional hydrological functions. Summary of the Invention

[0004] To address the technical problems existing in the prior art, embodiments of the present invention provide a method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions, comprising the following steps:

[0005] To achieve the above objectives, the present invention adopts the following technical solution: a method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions, comprising the following steps:

[0006] S1: Obtain the starting point and water level data of the groundwater conservation area, extract the three-day water level rise trend and locate the inflection point of change, record the interval with the restoration time point, detect the difference between evaporation and infiltration rates, and obtain the time lag period for the restoration intervention to take effect.

[0007] S2: Call the delayed period of the remediation intervention, select the boundary extension line, and determine whether the period is shorter than the response time limit. If it is shorter, move the boundary outward. If it is longer, freeze the boundary and mark the core area to obtain the dynamic boundary map of the pollution control area.

[0008] S3: Call the dynamic boundary map of the pollution control area, extract the consistency between water level gradient change and seepage direction, screen the continuous hydrodynamic area and sort the intervention order to obtain a priority sequence list of remediation intervention.

[0009] S4: Call the aforementioned priority sequence list for repair intervention, extract the difference between the upper and lower boundaries of the capillary water layer and calculate the fluctuation, detect the particle size and organic matter content, compare the fluctuation with the water retention capacity limit, and obtain the intervention concentration control section marking map;

[0010] S5: Call the intervention concentration control zone marker map, extract the pollutant distribution map and mark the high value segment, overlay the groundwater level rise trend map, determine the consistency of the water content trend, and obtain the dynamic formulation plan of soil remediation target.

[0011] As a further aspect of the present invention, the effective lag period of the remediation intervention includes the timing of the inflection point of water level change, the length of the lag interval, the magnitude of the rate difference, and the classification of lag types; the dynamic boundary map of the pollution control zone includes the boundary adjustment status, the location of the core retention area, the map correction range, and the response time comparison results; the priority sequence list of remediation intervention includes the frequency of gradient changes, the ranking of change rates, the seepage consistency ratio, the continuous hydrodynamic segment, and the intervention priority level; the intervention concentration control section marking map includes the fluctuation amplitude index, the estimated particle size ratio, the organic matter content level, the water retention capacity judgment result, and the restriction section identification; and the dynamic formulation scheme of soil remediation targets includes the distribution of high concentration segments, the direction of water content trend, the setting of remediation frequency, and the arrangement of time cycles.

[0012] As a further aspect of the present invention, the response time limit refers to the shortest time threshold required for hydrological changes after boundary expansion to generate a response to restoration intervention measures.

[0013] As a further aspect of the present invention, the step of screening continuous hydrodynamic regions and sorting the intervention order refers to identifying continuous hydrodynamic segments based on the consistency between the frequency of water level gradient changes and the seepage direction, and determining the intervention priority according to the level of water level abrupt change.

[0014] As a further aspect of the present invention, the specific steps of S1 are as follows:

[0015] S101: Obtain groundwater level monitoring data frames, extract three consecutive days of water level sequence, calculate daily increase, filter time periods of continuous water level rise, match the end date with the corresponding water level, and obtain a set of continuously rising water level markers.

[0016] S102: Call the continuously rising water level marker set, calculate the rate of increase, locate the time point corresponding to the maximum increase, determine the change in the direction of the increase before and after, filter the time when the trend turns, and obtain the time node of the water level change inflection point.

[0017] S103: Call the time node of the water level change inflection point and the timestamp of the repair start point, calculate the time difference between the two, remove the inflection point that is earlier than the repair start point, retain the time difference value and aggregate it to obtain the continuously rising water level marker set, and detect the difference between the evaporation rate and the infiltration rate to generate the repair intervention lag period.

[0018] As a further aspect of the present invention, the specific steps of S2 are as follows:

[0019] S201: Call the delayed period of the remediation intervention, select the original extension line position of the pollution control zone boundary, judge the difference between the delayed period and the boundary response time limit, filter the boundary segments with a delayed period shorter than the response time limit, and obtain the boundary extension adjustment segment set;

[0020] S202: Based on the set of boundary extension adjustment segments, perform an outward shift operation on the corresponding boundary, retain the boundary with a lag time longer than the response time limit and mark the core retention area, integrate the boundary changes, and generate a repair boundary adjustment layer;

[0021] S203: Call the repair boundary adjustment layer, overlay the original boundary map, identify the boundary change area, redraw the overall boundary line and mark the area attributes, extract the layer shape data, and generate the dynamic boundary map of the pollution control area.

[0022] As a further aspect of the present invention, the specific steps of S3 are as follows:

[0023] S301: Call the dynamic boundary map of the pollution control area, extract the water level data sequence between monitoring points, calculate the water level difference between adjacent monitoring points and count the frequency of changes in the positive and negative directions, and calculate the rate of change for the absolute value of continuous differences. Aggregate the frequency and rate data to obtain a set of water level gradient change indicators.

[0024] S302: Based on the water level gradient change index set, extract the seepage direction scalar value of adjacent monitoring segments, determine the consistency of direction between adjacent segments, aggregate consistent segments and arrange them in order, extract the set of segments with continuity, and generate a list of hydrodynamic continuity areas.

[0025] S303: Call the list of hydrodynamic continuity areas, extract the water level change difference within the area, classify and mark the difference value range, sort the area numbers from high to low according to the level, establish a mapping table between area and number and record the sorting position, and establish a priority sequence list for remediation intervention.

[0026] As a further aspect of the present invention, the specific steps of S4 are as follows:

[0027] S401: Call the marked segments in the repair intervention priority sequence list, extract the elevation values ​​of the upper and lower boundaries of the capillary water layer in the corresponding area, calculate the difference between the upper and lower boundaries, and calculate the difference between the maximum and minimum amplitudes of the continuous difference sequence to obtain the capillary water layer fluctuation amplitude value.

[0028] S402: Based on the capillary water layer fluctuation amplitude value, extract the soil particle size distribution and organic matter content data frames of the same segment, calculate the proportion of coarse particles and estimate the water storage capacity index in combination with the organic matter content, compare the index with the water storage capacity threshold, and generate a restricted area identification result set.

[0029] S403: Call the restricted area identification result set, locate the corresponding segment position, mark the area range in the pollution control area layer, and establish a mapping by combining the segment number and the control segment number, output the vector boundary and the marking attribute field, and generate the intervention concentration control segment marking map.

[0030] As a further aspect of the present invention, the specific steps of S5 are as follows:

[0031] S501: Call the area marked by the intervention concentration control section marker map, extract the pollutant concentration sequence in the soil depth direction of the corresponding area, identify the concentration peak position and make spatial annotation, combine the profile depth data to construct the map structure, and generate a vertical pollution concentration distribution map.

[0032] S502: Call the vertical pollution concentration distribution map, overlay the groundwater level rise trend change map, locate the concentration change section in the map and extract the corresponding spatial location, determine whether the groundwater water content trend direction at the location is consistent with the pollution trend direction, and obtain the trend matching relationship determination result.

[0033] S503: Based on the trend matching relationship determination results, set corresponding remediation cycles for trend-consistent and trend-inconsistent regions, establish a remediation frequency and time-history correspondence table based on spatial segment numbers, organize the intervention segment operation arrangements according to time sequence, and establish a dynamic formulation plan for soil remediation targets.

[0034] As a further aspect of the present invention, the water retention capacity threshold is the minimum water retention performance standard value used to determine whether a soil layer has the ability to retain water.

[0035] Compared with the prior art, the advantages and positive effects of the present invention are as follows:

[0036] In this invention, the intervention lag relationship is clarified by extracting the water level inflection point and the infiltration difference. The boundary response direction is adjusted according to the lag period to enhance the adaptability of the control range. The intervention priority order is constructed by combining the hydraulic gradient and the seepage trend to highlight the repair efficiency of continuous segments. The restricted area is determined by the water layer fluctuation and water retention capacity, and the concentration control position is accurately delineated. The repair frequency and time are set by superimposing the water level trend and the pollution distribution, so as to achieve dynamic coordination of intervention identification, boundary adjustment and repair rhythm. Attached Figure Description

[0037] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0038] Figure 1 This is a schematic diagram of the steps of the present invention;

[0039] Figure 2 This is a detailed schematic diagram of S1 of the present invention;

[0040] Figure 3 This is a detailed schematic diagram of S2 of the present invention;

[0041] Figure 4 This is a detailed schematic diagram of S3 of the present invention;

[0042] Figure 5 This is a detailed schematic diagram of S4 of the present invention;

[0043] Figure 6 This is a detailed schematic diagram of S5 of the present invention. Detailed Implementation

[0044] The technical solution of the present invention will now be described with reference to the accompanying drawings.

[0045] In embodiments of the present invention, words such as "exemplarily," "for example," etc., are used to indicate that something is an example, illustration, or description. Any embodiment or design described as "exemplary" in the present invention should not be construed as being more preferred or advantageous than other embodiments or designs. Specifically, the use of the word "exemplary" is intended to present the concept in a concrete manner. Furthermore, in embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one.

[0046] In the embodiments of this invention, the terms "image" and "picture" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning. Similarly, the terms "of," "corresponding (relevant)," and "corresponding" may sometimes be used interchangeably. It should be noted that, without emphasizing the distinction between them, they convey the same meaning.

[0047] In this embodiment of the invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1. When the difference is not emphasized, the meaning they express is the same.

[0048] To make the technical problems, technical solutions and advantages of the present invention clearer, a detailed description will be given below in conjunction with the accompanying drawings and specific embodiments.

[0049] Please see Figure 1 This invention provides a method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions, comprising the following steps:

[0050] S1: Obtain the starting point of the groundwater conservation area restoration and groundwater level monitoring data, extract the three-day water level rise trend and locate the inflection point of the maximum change, record the time interval between the inflection point and the restoration time point, detect the difference between the evaporation and vertical infiltration rate during the same period, compare the change range between the time interval and the rate difference, divide the lag relationship between intervention behavior and hydrological response, and obtain the lag period for the restoration intervention to take effect.

[0051] S2: Call the lag period of the remediation intervention, select the original extension line position of the pollution control zone boundary, and determine whether the lag period is shorter than the boundary response time limit. If it is shorter, move the remediation boundary outward. If it is longer, freeze the boundary and mark the core retention area. Correct the remediation treatment area corresponding to the boundary of the area to obtain the dynamic boundary map of the pollution control zone.

[0052] S3: Call the dynamic boundary map of the pollution control area, extract the frequency and rate of change of water level gradient between monitoring points, calculate the consistency ratio of seepage direction between adjacent segments and determine the trend direction, screen the hydrodynamic continuity area and sort the intervention order according to the water level change level, establish the order of intervention operations, and obtain the priority sequence list of remediation intervention.

[0053] S4: Call the marked segments in the priority sequence list of repair intervention, extract the difference between the upper and lower boundaries of the capillary water layer and calculate the fluctuation amplitude, detect the particle size ratio and organic matter content in the same segment to estimate the water retention capacity, compare the fluctuation amplitude with the retention capacity threshold and determine the restricted area, and obtain the intervention concentration regulation section marking map.

[0054] S5: Call the area marked by the intervention concentration control section marker map, extract the vertical distribution spectrum of pollutants in the soil layer and mark the high concentration range, overlay the groundwater level rise trend change map, determine whether the water content trend direction in the overlapping position is consistent, set the remediation frequency and time according to the trend matching situation, and obtain a dynamic soil remediation target formulation plan.

[0055] The lag period for remediation intervention includes the timing of water level change inflection points, lag interval length, rate difference magnitude, and lag type classification. The dynamic boundary map of the pollution control area includes the boundary adjustment status, core retention area location, map correction range, and response time comparison results. The priority sequence list for remediation intervention includes the frequency of gradient changes, change rate ranking, seepage consistency ratio, continuous hydrodynamic segments, and intervention priority level. The intervention concentration control zone marking map includes fluctuation amplitude indicators, particle size ratio estimation, organic matter content level, water retention capacity judgment results, and restricted zone identification. The dynamic formulation plan for soil remediation targets includes the distribution of high concentration segments, water content trend direction, remediation frequency setting, and time cycle arrangement.

[0056] Please see Figure 2 The specific steps of S1 are as follows:

[0057] S101: Obtain groundwater level monitoring data frames, extract three consecutive days of water level sequence, calculate daily increase, filter time periods of continuous water level rise, match the end date with the corresponding water level, and obtain a set of continuously rising water level markers.

[0058] Groundwater level monitoring data frames are typically acquired periodically by automatic water level monitoring equipment deployed in groundwater observation wells. Each record includes a timestamp, water level value, and station number. The data frame structure is arranged in ascending order by timestamp to ensure sequence continuity. Subsequently, time windows are divided into 3-day units, and a sliding window method with a step size of 1 day is used to extract water level sequences for consecutive three days. For each three-day sequence, the water level values ​​for the three days are recorded separately. The water level changes between the first and second days, and between the second and third days are calculated. The differences in water levels between the first and second days, and between the second and third days, are extracted, and the increase is calculated by direct subtraction. The calculation is performed group by group. For example, if the water level at a certain observation point is 3.12 meters on the first day, 3.24 meters on the second day, and 3.41 meters on the third day, the daily increases are 0.12 meters and 0.17 meters, respectively. When both increases are positive, it means that the water level is on an upward trend for two consecutive days. The timestamp and water level value corresponding to the third day are then grouped together as a record of the continuous rise endpoint in the current time window. After continuously traversing the entire set of monitoring data, multiple time nodes and water level values ​​of continuous water level rise can be generated. These paired values ​​constitute a continuous water level rise marker set, which is used to determine the inflection point time distribution of the water level change trend.

[0059] S102: Call the continuously rising water level marker set, calculate the rate of change of the increase, locate the time point corresponding to the maximum increase, judge the change of the direction of the increase before and after, filter the time when the trend turns, and obtain the time node of the water level change inflection point.

[0060] After arranging the continuously rising water level markers in chronological order, the water level values ​​and timestamps of two adjacent records are extracted. The increase in water level between the two records is calculated, and the current increase value is obtained. Then, the two adjacent increases are compared to determine whether the increase has an upward or downward trend. By continuously judging the changing trend of the difference between adjacent increases, it is possible to locate whether the increase shows a decreasing or increasing trend. For example, in the records of a certain observation station, the first continuously rising water level is 3.50 meters, the second is 3.68 meters, two days apart, with an increase of 0.18 meters. The next one is 3.79 meters, one day apart, with an increase of 0.11 meters. The increase has decreased from 0.18 to 0.11, which can be determined as a downward trend. By detecting this trend change process, the maximum increase can be identified. The time point where the amplitude occurs is compared with several increase records before and after it. If it changes from continuous increase to decrease, it indicates a turning point in the rising water level trend. The necessary condition for judging the turning point is that the direction of the increase changes, that is, from increase to decrease or reverses. The trend inflection point is identified by comparing the increase at fixed intervals. To avoid misjudgment, a minimum change threshold needs to be set. It is recommended to set it to a change amplitude of not less than 0.05 meters. Changes below this amplitude will be considered invalid fluctuations. In actual operation, this value can be set with reference to historical data statistics. For example, if the standard deviation of water level change throughout the year is 0.12 meters, then 0.05 meters can be used as the lower limit of the effective trend change. Finally, the time points that meet the above change characteristics are extracted to form the water level change inflection point time points.

[0061] S103: Call the time node of the water level change inflection point and the timestamp of the repair start point, calculate the time difference between the two, remove the inflection point that is earlier than the repair start point, retain the time difference value and aggregate it to obtain the set of continuously rising water level markers, and detect the difference between evaporation rate and infiltration rate to generate the lag period of the repair intervention.

[0062] Based on the set of inflection point time nodes for water level changes, each inflection point time is compared with the timestamp of the intervention start point to calculate the time interval between them. If the inflection point time is earlier than the intervention start point time, it is judged as a non-intervention-related inflection point and discarded. The remaining inflection point times are the potential response time nodes affected by the intervention. For each inflection point time that meets the conditions, the specific interval in days between it and the intervention start point is calculated and recorded as the basic data for the lagged response. For example, if the intervention start point time is May 1st and an inflection point is May 3rd, the corresponding lag days are 2 days. If the inflection point is April 29th, it is discarded because it is earlier than the intervention start point. To avoid calculation errors, it is necessary to... Dates are standardized to ensure that interval calculations are performed in a uniform time format. All eligible lag time points are then aggregated into a lag period set. To ensure the significance of identified lag responses, a lag threshold range can be set. For example, only lag inflection points within 1 to 10 days can be identified; those exceeding 10 days are considered indirect responses and excluded. This threshold can be set with reference to the typical effective time of interventions, combined with historical response records. If the historical average response time is 7 days and the maximum response time is 10 days, the lag analysis interval can be set between 1 and 10 days. Finally, the lag period for the intervention's effectiveness is output as the basic input data for subsequent analysis of the intervention's lag effect.

[0063] Please see Figure 3 The specific steps of S2 are as follows:

[0064] S201: Call the lag period of the remediation intervention, select the original extension line position of the pollution control zone boundary, judge the difference between the lag period and the boundary response time limit, filter the boundary segments with lag periods shorter than the response time limit, and obtain the set of boundary extension adjustment segments;

[0065] When calling the lag period of the repair intervention, it is necessary to traverse the lag days data, extract the spatial coordinate information corresponding to each lag period one by one, ensure that the period can be accurately mapped to the boundary area of the pollution control zone geographically, select the original set extension line position in the boundary layer. The original set extension line is a buffer segment reserved for the outer edge of the initial pollution control range. Each extension segment has a unique number and its corresponding spatial coordinate information, and the corresponding boundary response time limit is marked. The boundary response time limit refers to the theoretical minimum time for the pollutant impact to spread to this boundary segment set during the original assessment process. This response time limit can be determined comprehensively according to parameters such as the source release intensity, regional groundwater flow velocity, and formation permeability. It is recommended to set the interval as 6 to 10 days. When more than 80% of the data in the historical response time frequency distribution of a certain area is concentrated between 6 and 8 days, the upper limit can be taken as the response threshold. For example, if the boundary response time limit is set to 8 days, the difference between each lag period T_delay and the response time limit T_response of the boundary segment where it is located is judged. When T_delay is less than T_response, that is, T_delay < T_response, it means that the pollution water level reaction speed at this point is relatively fast and has not reached the response time that the boundary segment should have, which is regarded as the boundary segment not fully responding to the pollution diffusion and meets the conditions for extension adjustment. Otherwise, it is kept unchanged. For example, if the lag time of a certain point is 5 days and the response time limit of the boundary segment where it is located is 8 days, the difference is days, meeting the adjustment conditions; if the lag time is 9 days, it does not meet the screening conditions and is discarded. After repeating this process to traverse all lag points and boundary segments, finally, all boundary segments that satisfy T_delay less than T_response are extracted to generate a boundary extension adjustment section set.

[0066] S202: According to the boundary extension adjustment section set, perform an outward movement operation on the corresponding boundary, retain the boundaries with lag periods longer than the response time limit, mark the core retention area, integrate the boundary changes, and generate a repair boundary adjustment layer;

[0067] Based on the numbering and spatial location provided by the boundary extension adjustment segment set, each boundary segment is moved outward. The outward movement should be along the direction of pollution source diffusion, typically taking the normal direction of the boundary segment. The outward movement distance needs to comprehensively consider factors such as the actual migration capacity of pollutants, groundwater flow velocity, and terrain slope. It is recommended to set the offset range to 50 meters to 300 meters. In high-permeability areas, it is recommended to use the upper limit of the offset value, while in low-permeability clay areas, the lower limit can be used to ensure that the boundary extrapolation conforms to geological conditions. For example, if the outward movement distance of a high-permeability segment is set to 250 meters, then the original boundary coordinate point is translated 250 meters along the normal direction to generate a new boundary vertex, constructing a new boundary. After the boundary segment is processed, the next one is processed. During the processing, it is simultaneously determined whether the lag time of the boundary segment exceeds the response time limit. If the lag time is greater than or equal to the response time limit, no adjustment is made, and the boundary segment is directly marked as the core retention area. The identifier field "Retention Area = Yes" is set in the attribute table. In addition, to facilitate subsequent layer integration, the "Adjustment Area = Yes" field is added to all boundary segments that have undergone the outward movement operation. After all boundary segment operations are completed, the outward-moved boundary segments are merged with the retained boundary segments. Based on their connection relationship, a new continuous boundary line is reconstructed to achieve seamless splicing of layer boundary segments. The overall result is output as a repair boundary adjustment layer.

[0068] S203: Call the repair boundary adjustment layer, overlay the original boundary map, identify the boundary change area, redraw the overall boundary line and label the area attributes, extract the layer shape data, and generate the dynamic boundary map sheet of the pollution control area;

[0069] After calling the boundary adjustment layer, it needs to be spatially overlaid with the original boundary map. A joint spatial analysis is performed on the two layers using a GIS platform to identify boundary segments with non-overlapping coordinates as changed areas. These changed segments are then re-stitched with adjacent unchanged segments to construct a complete new boundary line. The redrawn boundary line needs to be displayed on the graphical interface with different colors or line types for manual identification. Simultaneously, all boundary segments are labeled with their attribute status, including fields such as "adjusted," "retained," "offset distance," and "adjustment direction." This attribute data is bound to the layer spatial data and stored uniformly in the boundary layer data table. When extracting layer shape data, the start and end point coordinates of each boundary segment must be extracted to ensure the original boundary segment geometry can be restored. Finally, all labeled boundary lines are merged and output as a dynamic boundary map of the pollution control area. This map should be in a standard GIS vector data format, such as shp or geojson, and include metadata information such as layer timestamps and version numbers. Each boundary segment in the map contains complete attribute fields and geometric information, which can be used for subsequent monitoring, analysis, and display.

[0070] Please see Figure 4 The specific steps of S3 are as follows:

[0071] S301: Call the dynamic boundary map of the pollution control area, extract the water level data sequence between monitoring points, calculate the water level difference between adjacent monitoring points and count the frequency of changes in the positive and negative directions, and calculate the rate of change for the absolute value of continuous differences. Aggregate the frequency and rate data to obtain a set of water level gradient change indicators.

[0072] First, the geographic coordinates and corresponding water level time series data of all monitoring points within the map coverage area need to be extracted. Each monitoring point should have a uniform time step, such as daily or hourly, to ensure comparability in the time dimension. The extracted data are used to construct a two-dimensional sequence matrix, with rows representing time and columns representing different monitoring points. Based on this, spatially adjacent monitoring point pairs are selected and paired according to the principle that the spatial distance does not exceed 200 meters or that there is an underground water channel between them. For example, monitoring points P1 and P2 are 120 meters apart, which constitutes a pair of effective comparison units. For each pair of monitoring points, the water level difference is calculated at each time step. The water level of the downstream point is subtracted from the water level of the upstream point to form a difference sequence. At the same time, the positive and negative directions of the difference are recorded. If the difference is positive, it indicates that... The water flow direction is from P1 to P2. If it is negative, the direction is reversed. Iterate through all time steps and record the frequency of positive and negative directions for each pair of monitoring points. If the positive direction occurs 8 times and the negative direction occurs 2 times, then the positive frequency of the monitoring point pair is recorded as 80% and the negative frequency as 20%. Then, the absolute value of the difference between each time step is subtracted from the adjacent time step to obtain the water level change rate. The water level change rate is the change amplitude of the absolute value of the water level difference between adjacent time steps. All rate values ​​are aggregated by taking the average, maximum value and standard deviation to form the change rate feature value of the monitoring point pair. Finally, the positive and negative frequency and change rate feature of all monitoring point pairs are integrated and summarized to output the water level gradient change index set.

[0073] S302: Based on the water level gradient change index set, extract the seepage direction scalar value of adjacent monitoring segments, determine the consistency of direction between adjacent segments, aggregate consistent segments and arrange them in order, extract the set of segments with continuity, and generate a list of hydrodynamic continuity areas.

[0074] Based on the water level difference and direction change information extracted from each pair of monitoring points using the water level gradient change index, the seepage direction scalar values ​​of adjacent point pairs constituting a monitoring segment are summarized. The seepage direction scalar values ​​are set according to the order of water level at the monitoring points: if the water level at the upstream point is higher than that at the downstream point, the direction scalar is set to positive; otherwise, it is set to negative. The direction scalar values ​​of all monitoring segments are compared. If the signs of the direction scalar values ​​between adjacent segments are consistent, their directions are considered consistent. If the scalar values ​​of three adjacent segments are all positive, the direction of that segment is considered continuous and consistent, and they are aggregated into a continuous segment. If the direction changes abruptly at some point, such as the scalar value changing from positive to negative, it is considered... If the continuity is interrupted, the segment will not be included in the current aggregation segment. A new continuous segment will be constructed from the beginning. For example, if the direction of segment AB is positive and the direction of BC is also positive, then ABC constitutes a continuous segment. If the direction of CD becomes negative, the new segment will start from D. The entire index set will be traversed to complete the aggregation of directional continuous segments. For each aggregated segment, the start point, end point and intermediate monitoring points will be recorded and arranged in spatial connection order to form a traceable path segment. All segments that meet the continuous direction condition will be integrated to generate a list of hydrodynamic continuous areas. This list should include the spatial location, length, start and end monitoring point numbers and direction scalar of each continuous segment.

[0075] S303: Call the list of hydrodynamic continuity areas, extract the water level change difference within the area, classify and mark the difference range, sort the area numbers from high to low according to the level, establish a mapping table between area and number and record the sorting position, and establish a priority sequence list for remediation intervention.

[0076] After retrieving the list of hydrodynamic continuity zones, a secondary extraction of water level time series data within each continuity zone is required. Particular attention should be paid to whether there are abrupt changes in water level differences within the continuity zone, i.e., the water level difference between two adjacent monitoring points at a certain time step is significantly higher than the regional average. To address this, the average and standard deviation of water level differences within the zone should be calculated first, and then abrupt water level changes should be screened. It is recommended that the threshold for determining abrupt changes be set to the regional average difference plus twice the standard deviation. For example, if the average water level difference in a continuity zone is 0.2 meters and the standard deviation is 0.05 meters, then the threshold should be set to 0.3 meters. Points exceeding this value are considered abrupt changes. After extracting all abrupt changes, their difference range should be analyzed. The process involves tiered processing, with recommended tiers divided into three groups: 0.3-0.5 meters, 0.5-0.8 meters, and above 0.8 meters. These are labeled as Tier 1, Tier 2, and Tier 3 mutation levels, with each level corresponding to a different intervention intensity. Regions are numbered and sorted from highest to lowest mutation level. The priority of these numbers can be set by multiplying the mutation level by a length-weighted method, meaning that regions with longer and more severe mutations are ranked higher. The ranking results are then linked to the corresponding regions to establish a one-to-one mapping table between regions and numbers. The position of each region within the entire ranking is recorded. The final output is a priority sequence list for repair interventions, which can be used for subsequent repair task scheduling.

[0077] Please see Figure 5 The specific steps of S4 are as follows:

[0078] S401: Call the marked segments in the priority sequence list of repair intervention, extract the elevation values ​​of the upper and lower boundaries of the capillary water layer in the corresponding area, calculate the difference between the upper and lower boundaries, and calculate the difference between the maximum and minimum amplitudes of the continuous difference sequence to obtain the capillary water layer fluctuation amplitude value.

[0079] When calling up the marked segments in the priority sequence list for repair intervention, the spatial number of the corresponding segment should first be retrieved one by one from the regional map sheet or database. Using this number as the search criterion, the elevation values ​​of the upper and lower boundaries of the capillary aquifer within the segment's coverage area should be extracted. This data comes from borehole stratification records or ground-penetrating radar profile data, and is generally recorded in segments at intervals of one meter to ensure that the units and datum planes of the upper and lower boundaries are consistent. Then, a one-to-one difference calculation is performed on the upper and lower boundary elevations at all sampling points to calculate the capillary aquifer thickness at each point. Finally, the thickness difference sequence is constructed by sorting the data by time or space. For this sequence, the maximum and minimum values ​​are extracted, representing the maximum and minimum thickness of the capillary water layer in the time or space scale, respectively. Then, a difference operation is performed on these two values ​​to obtain the capillary water layer fluctuation amplitude value. This value reflects the degree of vertical change of the capillary water layer within a certain period or section. For example, if the upper boundary elevation range of a certain segment is 63.0 to 64.2 meters and the lower boundary is 60.5 to 62.8 meters, then the minimum thickness is 1.2 meters, the maximum thickness is 3.0 meters, and the fluctuation amplitude value is 1.8 meters. This provides key basic parameters for subsequent identification of restricted areas and formulation of control strategies.

[0080] S402: Based on the capillary water layer fluctuation amplitude value, extract the soil particle size distribution and organic matter content data frames of the same segment, calculate the proportion of coarse particles and estimate the water storage capacity index in combination with the organic matter content, compare the index with the water storage capacity threshold, and generate a restricted area identification result set.

[0081] Extract soil particle size distribution data and organic matter content data frames within the same spatial range. Data sources include in-situ sampling particle size analysis reports or geotechnical test results. Particle size distribution data is typically expressed as a percentage by mass, covering particle sizes such as coarse sand, medium sand, fine sand, silt, and clay. The coarse-grained portions of these sizes are summarized by adding the proportions of coarse and medium sand to obtain the coarse-grained component percentage. This is then combined with the organic matter content data for the area, which is generally expressed as a percentage of total organic carbon. These two factors are combined to estimate the water storage capacity index. The estimation process involves multiplying the coarse-grained proportion by the organic matter content to obtain a dimensionless value. For example, if a segment has a coarse-grained proportion of 55% and an organic matter content of 2.5%, then the water storage capacity index is... The capacity index is 1.375. To make a judgment, the water storage capacity index needs to be compared with the water storage capacity threshold. The setting of this threshold value should be based on statistical analysis of the historical soil water storage performance of the pollution control area. It is recommended to select typical representative values ​​from historical data and set them according to percentiles. For example, if the 25th percentile water storage capacity index is 1.5, then 1.5 can be set as the threshold, indicating that segments with a value lower than 1.5 have weak water storage capacity under most conditions and need to be included in the restriction range. The index of all segments should be traversed and judged. The corresponding segment number and attribute information of the records with a value lower than the threshold should be extracted together to form a restricted area identification result set. The dataset should include fields such as segment number, index value, and whether it is lower than the threshold.

[0082] S403: Call the restricted area identification result set, locate the corresponding segment position, mark the area range in the pollution control area layer, and establish a mapping by combining the segment number and the control segment number, output the vector boundary and the marking attribute field, and generate the intervention concentration control segment marking map;

[0083] After retrieving the restricted area identification result set, the spatial extent of each segment should be located in the pollution control zone layer according to its segment number, and its spatial geometric boundary in the layer should be determined. The extracted spatial extent exists in the form of polygons or polylines. It is necessary to ensure that the layer coordinate system is consistent to guarantee accurate projection. At the same time, the restriction mark field value should be written into the layer attribute table. Setting the mark value to "1" indicates that the area is identified as a restricted segment. All marked areas should be sorted by segment number and mapped with the original boundary segment number of the control zone. The mapping method is to map each restricted segment to its corresponding pollution control zone segment. For example, if a restricted segment number is L-P009 and the corresponding control zone segment number is Z3-7, it will be recorded in the mapping table as... After establishing the mapping table, it is attached to the layer attribute table as an external association field, ultimately forming a layer file containing vector boundary graphics and attribute data. The layer attribute fields should at least include fields such as segment number, restriction marker, control section number, and water storage capacity index value. The output format is a standard GIS-supported vector data format such as SHP or GeoJSON, and can be used for spatial calling and regional identification in the downstream intervention area regulation strategy module, ultimately generating an intervention concentration regulation section marking map.

[0084] Please see Figure 6 The specific steps of S5 are as follows:

[0085] S501: Call the area marked by the intervention concentration control section marker map, extract the pollutant concentration sequence in the soil depth direction of the corresponding area, identify the concentration peak position and make spatial annotation, combine the profile depth data to construct the map structure, and generate a vertical pollution concentration distribution map.

[0086] When accessing the areas marked on the intervention concentration control zone map, the vector boundaries of each marked intervention area in geographic space should be located one by one. The pollution monitoring borehole or sampling point numbers within these boundaries should be extracted, and the corresponding borehole's soil depth-direction pollutant concentration sequence should be retrieved accordingly. This sequence records the pollutant content values ​​at different layers from the surface to the final borehole depth, typically in mg / kg or mg / L. Concentration values ​​are recorded at 0.5-meter or 1-meter intervals to construct a depth-concentration one-dimensional data curve. The concentration values ​​before and after each point on the curve are compared, and the point with the highest value is extracted as the peak concentration position for that borehole. The depth value and corresponding concentration at this position are recorded. The location is marked in space by combining the horizontal coordinates of the borehole with the peak depth to form a three-dimensional spatial coordinate system, generating a set of visualized points. After extracting the peak positions of all boreholes, geological profile data of each borehole is further superimposed to obtain information such as the burial depth range and lithological layering structure of each stratum in the profile. The relative position of the concentration peak in the profile is located. Then, combined with the three-dimensional distribution information of boreholes in the whole area, a spatial interpolation structure is constructed. The concentration of each point is interpolated to generate a planar map, which is drawn layer by layer according to the depth dimension to form a vertical pollution concentration distribution map. This map is output in 2D profile and 3D stereoscopic map formats, which is convenient for subsequent superposition and analysis with other environmental factors.

[0087] S502: Call the vertical pollution concentration distribution map, overlay the groundwater level rise trend change map, locate the concentration change section in the map and extract the corresponding spatial location, determine whether the groundwater water content trend direction at the location is consistent with the pollution trend direction, and obtain the trend matching relationship determination result.

[0088] After retrieving the vertical pollution concentration distribution map, it needs to be spatially overlaid with the groundwater level rise trend map. The groundwater level map records the groundwater level elevation at various points during different monitoring periods, and calculates the trend direction according to the time series. If the water level data for a certain area shows an upward trend for four consecutive periods, the trend direction is upward and marked as positive; conversely, if it is downward, it is marked as negative. Match the depth and horizontal position corresponding to the high concentration range in the map with the same position in the groundwater level map, extract the groundwater level trend indicator at that position, and determine the consistency of the direction with the pollution trend direction of that concentration range. The direction of pollution trend is determined based on the change of concentration value in the depth direction. For example, if the concentration value increases layer by layer from top to bottom in a certain concentration range, it is determined to be a downward migration trend. If the groundwater level is rising and the pollution is migrating downward, the trend direction is consistent; otherwise, it is inconsistent. The above steps are repeated to traverse all high concentration value ranges in the map to form a location-trend matching result set. In this result set, the spatial coordinates of each location, pollution migration trend, groundwater trend, and trend consistency judgment result are marked. If the trend is consistent, it is marked as "1" and if it is inconsistent, it is marked as "0". Finally, the trend matching relationship judgment result is generated.

[0089] S503: Based on the trend matching relationship determination results, set corresponding remediation cycles for trend-consistent and trend-inconsistent areas, establish a remediation frequency and time history correspondence table based on spatial segment number, organize the intervention segment operation arrangements according to time sequence, and establish a dynamic soil remediation target formulation plan.

[0090] Based on the trend matching results, an intervention frequency is set for the area where each marker point is located. If the trend is consistent (marked as "1"), it is designated as a synchronous response zone, and a higher frequency of repair cycle is assigned to this zone, such as once every 15 days. If the trend is inconsistent (marked as "0"), it is designated as a lagging response zone, and a delayed intervention cycle is set, such as once every 30 days. The frequency value is set with reference to the pollution response cycle of similar sections in historical repair projects. If the average repair cycle of the synchronous zone in historical data is 14 to 17 days, then the midpoint of 15 days can be taken as the set value. The lagging zone is extended by double the set value. The intervention period is 30 days. Subsequently, records of the corresponding remediation frequency and remediation cycle are generated for each segment number. The segment number is used as the primary key, and the remediation frequency (e.g., "15d") is used as the value item to establish a table corresponding to the remediation frequency and time schedule. This table is in a structured data format, and the fields include segment number, trend judgment, remediation cycle days, remediation frequency level, etc. Finally, the remediation arrangements of different segments are sorted by time axis, and segments with the same frequency are grouped into the same execution batch. The execution order is marked in the table, and the output is a complete operation arrangement table for the intervention segment. Based on this, a dynamic soil remediation target formulation plan is established. The plan should cover all identified segments and associate them with their intervention time schedule plans.

[0091] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions, characterized in that, Includes the following steps: S1: Obtain the starting point of the groundwater conservation area restoration and groundwater level monitoring data, extract the three-day continuous water level rise trend and locate the inflection point of the maximum change, record the time interval between the inflection point and the restoration time point, detect the difference between the evaporation and vertical infiltration rate during the same period, compare the change range between the time interval and the rate difference, divide the lag relationship between intervention behavior and hydrological response, and obtain the lag period of restoration intervention effectiveness. S2: Call the lag period of the remediation intervention, select the original extension line position of the pollution control zone boundary, and determine whether the lag period is less than the boundary response time limit. If it is less, move the boundary segment outward along the direction of pollution source diffusion. If it is greater than or equal to, freeze the boundary segment and mark it as the core retention area, correct the regional boundary, and obtain the dynamic boundary map of the pollution control zone. The original extension line is a buffer segment reserved at the outer edge of the initial pollution control range. Each extension segment has a unique number and its corresponding spatial coordinate information, and is marked with the corresponding boundary response time limit. The boundary response time limit refers to the theoretical shortest time threshold for the impact of pollutants to spread to the boundary segment, which is determined by comprehensively considering the pollution source release intensity, regional groundwater flow velocity and stratum permeability. S3: Call the dynamic boundary map of the pollution control area, extract the frequency and rate of change of water level gradient between monitoring points, calculate the consistency ratio of seepage direction between adjacent segments and determine the trend direction, screen the hydrodynamic continuity area and sort the intervention order according to the water level change level, establish the order of intervention operations, and obtain the priority sequence list of remediation intervention. S4: Call the marked segments in the priority sequence list of repair intervention, detect the particle size ratio and organic matter content in the same segment, estimate the water storage capacity index, compare the water storage capacity index with the water storage capacity threshold, when the water storage capacity index is lower than the water storage capacity threshold, identify the segment as a restricted area, generate a restricted area identification result set, and obtain the intervention concentration regulation section marking map; S5: Call the area marked by the intervention concentration control section marker map, extract the vertical distribution spectrum of pollutants in the soil layer, identify the location of the concentration peak and mark it spatially, overlay the groundwater level rise trend map, determine whether the groundwater moisture content trend direction in the overlapping position is consistent with the pollution trend direction, set the first remediation cycle for the area with consistent trend, set the second remediation cycle for the area with inconsistent trend, the first remediation cycle is shorter than the second remediation cycle, establish a remediation frequency and time history correspondence table according to the segment number, organize the operation arrangement of the intervention segment according to the time sequence, and obtain the dynamic formulation plan of soil remediation target.

2. The method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions according to claim 1, characterized in that, The effective lag period of the remediation intervention includes the timing of the inflection point of water level change, the length of the lag interval, and the rate difference; the dynamic boundary map of the pollution control zone includes the boundary adjustment status, the location of the core retention area, the map correction range, and the response time limit judgment result; the priority sequence list of remediation intervention includes the frequency of gradient changes, the ranking of change rates, the seepage consistency ratio, the continuous hydrodynamic segment, and the intervention priority level; the intervention concentration control section marking map includes the particle size ratio, organic matter content level, the result of the restricted area identification, and the restricted section marker; the dynamic formulation plan of soil remediation targets includes the distribution of concentration peak locations, the direction of water content trend, the setting of remediation frequency, and the arrangement of time cycles.

3. The method for dynamically formulating soil remediation targets based on the function of protecting and conserving groundwater, as described in claim 1, is characterized in that... The specific steps of S1 are as follows: S101: Obtain groundwater level monitoring data frames, extract three consecutive days of water level sequence, calculate daily increase, filter time periods of continuous water level rise, match the end date with the corresponding water level, and obtain a set of continuously rising water level markers. S102: Call the continuously rising water level marker set, calculate the rate of increase, locate the time point corresponding to the maximum increase, determine the change in the direction of the increase before and after, filter the time when the trend turns, and obtain the time node of the water level change inflection point. S103: Call the time node of the water level change inflection point and the timestamp of the repair start point, calculate the time difference between the two, remove the inflection point that is earlier than the repair start point, retain the time difference value and aggregate it, detect the difference between the evaporation rate and the infiltration rate, and generate the lag period of the repair intervention.

4. The method for dynamically formulating soil remediation targets based on protecting groundwater conservation functions according to claim 1, characterized in that, The specific steps of S2 are as follows: S201: Call the delayed period of the remediation intervention, select the original extension line position of the pollution control zone boundary, judge the difference between the delayed period and the boundary response time limit, filter the boundary segments with a delayed period less than the response time limit, and obtain the boundary extension adjustment segment set; S202: Based on the set of boundary extension adjustment segments, perform an outward shift operation on the corresponding boundary segments, retain the boundary segments with a lag time greater than or equal to the response time limit and mark the core retention area, merge the outward shifted boundary segments with the retained boundary segments, reconstruct a new continuous boundary line based on their connection relationship, and output the overall result as a repair boundary adjustment layer. S203: Call the repair boundary adjustment layer, overlay the original boundary map, identify the boundary change area, redraw the overall boundary line and mark the area attributes, extract the layer shape data, and generate the dynamic boundary map of the pollution control area.

5. The method for dynamically formulating soil remediation targets based on the function of protecting and conserving groundwater, as described in claim 1, is characterized in that... The specific steps for S3 are as follows: S301: Call the dynamic boundary map of the pollution control area, extract the water level data sequence between monitoring points, calculate the water level difference between adjacent monitoring points and count the frequency of changes in the positive and negative directions, and calculate the rate of change for the absolute value of continuous differences. Aggregate the frequency and rate data to obtain a set of water level gradient change indicators. S302: Based on the water level gradient change index set, extract the seepage direction scalar value of adjacent monitoring segments, determine the consistency of direction between adjacent segments, aggregate consistent segments and arrange them in order, extract the set of segments with continuity, and generate a list of hydrodynamic continuity areas. S303: Call the list of hydrodynamic continuity areas, extract the water level change difference within the area, classify and mark the difference value range, sort the area numbers from high to low according to the level, establish a mapping table between area and number and record the sorting position, and output the priority sequence list of repair intervention.

6. The method for dynamically formulating soil remediation targets based on the function of protecting and conserving groundwater sources according to claim 1, characterized in that, The specific steps of S4 are as follows: S401: Call the marked segment in the repair intervention priority sequence list, extract the elevation values ​​of the upper and lower boundaries of the capillary water layer in the corresponding area, calculate the difference between the upper and lower boundaries, and calculate the difference between the maximum and minimum amplitudes of the continuous difference sequence to obtain the capillary water layer fluctuation amplitude value. S402: Based on the capillary water layer fluctuation amplitude value, extract the soil particle size distribution and organic matter content data frames of the same segment, calculate the proportion of coarse particles and estimate the water storage capacity index in combination with the organic matter content; the water storage capacity threshold is the minimum water storage performance standard value used to determine whether the soil layer has the ability to retain water, compare the water storage capacity index with the water storage capacity threshold, when the water storage capacity index is lower than the water storage capacity threshold, identify the segment as a restricted area, and generate a restricted area identification result set; S403: Call the restricted area identification result set, locate the corresponding segment position, mark the area range in the pollution control area layer, and establish a mapping by combining the segment number and the control area number, output the vector boundary and the marking attribute field, and generate the intervention concentration control area marking map.

7. The method for dynamically formulating soil remediation targets based on the function of protecting and conserving groundwater sources according to claim 1, characterized in that, The specific steps of S5 are as follows: S501: Call the area marked by the intervention concentration control section marker map, extract the pollutant concentration sequence in the soil depth direction of the corresponding area, identify the concentration peak position and make spatial annotation, combine the profile depth data to construct the map structure, and generate a vertical pollution concentration distribution map. S502: Call the vertical pollution concentration distribution map, overlay the groundwater level rise trend change map, locate the concentration change section in the map and extract the corresponding spatial location, determine whether the groundwater water content trend direction at the location is consistent with the pollution trend direction, and obtain the trend matching relationship determination result. S503: Based on the trend matching relationship determination result, a first repair cycle is set for the trend-consistent area, and a second repair cycle is set for the trend-inconsistent area. The first repair cycle is shorter than the second repair cycle. A repair frequency and time history correspondence table is established according to the segment number. The repair arrangements of different segments are sorted according to the time axis. Segments with the same repair frequency are grouped into the same execution batch, and the execution order is marked in the table. An intervention segment operation arrangement table is output, and a dynamic soil remediation target formulation plan is established accordingly.

Images