Regional groundwater environmental risk assessment system and methods considering the vulnerability of risk receptors

Abstract

This invention discloses a regional groundwater environmental risk assessment system and method considering the vulnerability of risk receptors, relating to the field of hydrogeology. The system includes the following steps: S100, establishing a cross-timescale monitoring baseline, collecting groundwater recharge intensity sequences and hydrochemical sequences, extracting sedimentary pollution fingerprints through multidimensional correlation, and identifying germination areas with re-release characteristics; S200, under the constraint of the cross-timescale monitoring baseline, calculating the temporal coherence and phase residuals between the recharge disturbance sequence and the sedimentary pollution fingerprint, identifying key spatiotemporal locations triggering pollution release, and constructing a pollution release causal residual spectrum. This invention extracts pollution fingerprints through cross-timescale monitoring, identifies high-risk release areas, reconstructs pollution pathways by combining coherence analysis and phase inversion, sets threshold surfaces and precursor scores, generates pre-suppression parameters, and achieves pollution limiting and closed-loop control through frequency shifting and time inversion, thereby improving groundwater risk response and management capabilities.

Description

Technical Field

[0001] This invention relates to the field of hydrogeology, specifically to a regional groundwater environmental risk assessment system and method that considers the vulnerability of risk receptors. Background Technology

[0002] "Regional groundwater environmental risk assessment considering the vulnerability of risk receptors" refers to a risk assessment of potential pollution or damage to the groundwater environment in a given area. This assessment not only analyzes the intensity of pollution sources, groundwater transport characteristics, and environmental exposure pathways, but also comprehensively incorporates the vulnerability differences of affected entities (i.e., risk receptors, such as drinking water wells, agricultural irrigation water sources, and ecological wetland systems). Vulnerability here includes factors such as the receptor's sensitivity to pollutants, the availability of alternative water sources, and its capacity to withstand pollution. By combining risk receptor vulnerability with pollution risk sources and groundwater environmental characteristics, the spatial differences and severity of groundwater environmental risks can be more accurately reflected, thereby enabling regional-scale risk zoning, priority remediation zone designation, and scientific decision support. This approach emphasizes a holistic approach, starting from the "source-transport-receptor" chain, highlighting the impact of receptor characteristics on risk outcomes, and is more in line with the practical needs of environmental management and public health protection.

[0003] The existing technology has the following shortcomings: In existing technologies, some pollutants remain relatively dormant in groundwater sedimentary layers due to long-term adsorption, complexation, or depositional fixation, showing no significant migration or exposure risk in the short term. However, under the influence of a single, intense recharge event (such as rainstorm infiltration, sudden flood recharge, or unconventional artificial water replenishment), pollutants that were originally stable in the sedimentary layer may be rapidly leached and released into the groundwater flow field, forming a high-concentration transient pollution pulse. Such sudden releases often exceed the prediction range of conventional risk assessment models, causing vulnerable receptors to encounter pollution exposure exceeding their carrying capacity threshold in a very short time, resulting in delayed risk control measures and serious consequences.

[0004] The information disclosed in the background section is only intended to enhance the understanding of the background of this disclosure, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0005] The purpose of this invention is to provide a regional groundwater environmental risk assessment system and method that takes into account the vulnerability of risk receptors, so as to solve the problems in the background art mentioned above.

[0006] To achieve the above objectives, the present invention provides the following technical solution: a method for assessing the environmental risk of regional groundwater considering the vulnerability of risk receptors, comprising the following steps: S100 establishes cross-timescale monitoring baselines, collects groundwater recharge intensity sequences and hydrochemical sequences, extracts sedimentary pollution fingerprints through multidimensional correlation, and identifies germination areas with re-release characteristics; S200, under the constraint of cross-timescale monitoring baseline, calculates the temporal coherence and phase residual between the replenishment perturbation sequence and the deposited contamination fingerprint, identifies the key spatiotemporal locations that trigger contamination release, and constructs the contamination release causal residual spectrum; S300 constructs a phase conjugate inversion chain based on causal residual spectrum, combines historical data of supply disturbance for counterfactual playback, reconstructs pollution release path, and forms transient release threshold surface and precursor score; S400, based on transient release threshold surface and precursor score, injects feedforward constraint field into high-risk streamline region to regulate extraction intensity, replenishment rhythm and boundary water level, and generates pre-suppression control parameters. Under the influence of pre-suppression control parameters, the S500 initiates a time-inversion conformal mechanism, applies frequency misalignment and time-scale traction to adjust the groundwater flow line, and implements amplitude-limited rewriting of the release path to achieve pollution pulse reduction and risk closed-loop control.

[0007] Preferably, step S100 includes: In areas sensitive to groundwater pollution, a multi-layered and multi-depth observation well network is deployed to collect recharge intensity sequences and water chemistry sequences. The collected recharge intensity and hydrochemical sequences were normalized to analyze the synergistic relationship between recharge changes and hydrochemical responses, and recharge response features were extracted. Based on the replenishment response characteristics, the characteristics of contamination deposition are extracted to construct a fingerprint map of depositional contamination. By combining sedimentary pollution fingerprint information with groundwater flow direction, geological structure and recharge path for spatial mapping and cluster analysis, high-risk locations for pollutant re-release are identified, forming spatial triggering units for pollution release risk.

[0008] Preferably, step S200 includes: The recharge intensity sequence obtained from the observation wells is time-aligned with the pollution response sequence to establish a correspondence between recharge and response; Based on the correspondence between supply and response, time-frequency coherence features are extracted to analyze the linkage between supply disturbance and pollution response; Phase residual analysis is introduced based on time-frequency coherence characteristics to characterize the propagation delay behavior of disturbances and identify the response mechanisms of different pollutant releases; By combining coherence characteristics and phase residuals, a causal residual spectrum of pollution release is constructed, forming a dynamic causal structure between perturbation input and pollution response.

[0009] Preferably, when constructing the causal residual spectrum of pollution release, the replenishment disturbance and pollution response paths of different time windows in multiple observation wells are superimposed to form a spatiotemporal projection map of pollution release, which is used to identify the spatial propagation trajectory of key disturbance events and pollution responses.

[0010] Preferably, step S300 includes: The disturbance events and pollution response sequences are reverse-ordered to form the reverse propagation trajectory of the disturbance response; A phase conjugation processing mechanism is introduced based on the reverse propagation trajectory to construct a time-symmetric mapping trajectory for identifying sensitive areas of pollution release; Based on disturbance path, response path and conjugate trajectory, a three-dimensional groundwater release channel map is constructed to reconstruct the pollution release path; Based on the response characteristics of disturbance events, pollutant release threshold parameters are extracted and a precursor scoring factor system is established to form a quantitative basis for identifying pollution release risks.

[0011] Preferably, when constructing the precursor scoring factor system, the scoring factors include the frequency of water level fluctuations, the slope of pH changes, and the rate of increase in conductivity. By assigning different weights, a comprehensive score is formed to identify the trend of pollution release risk in advance.

[0012] Preferably, step S400 includes: By overlaying the scoring results with the groundwater flow direction map, high-risk groundwater flow line areas are identified and control units are delineated. In high-risk areas, analyze the driving relationship between extraction behavior and pollution response, adjust extraction intensity and cycle, and form extraction control parameters; A supply rhythm regulation model was constructed by combining the precursor score results, and supply intensity and frequency limits were set to form supply regulation parameters. A water level control strategy is constructed at the boundary of high-risk areas, linking extraction parameters, replenishment parameters, and boundary control results to generate a set of pre-inhibition control parameters for dynamic intervention.

[0013] Preferably, the linkage and integration of extraction control parameters, replenishment control parameters and boundary control results are based on the risk level classification standard set in the pollutant transient release threshold surface, and the extraction intensity adjustment range, replenishment rhythm intervention cycle and boundary water level variation range are set according to different risk levels.

[0014] Preferably, step S500 includes: A time-reversal driving field is constructed based on historical pollution response data to identify key transmission channels in the pollution release pathway; Based on the time-reversal driving field, a golden ratio frequency shifting strategy is implemented to regulate the extraction and replenishment rhythm in order to disturb the resonance of the pollution release cycle. A dual-mirror timescale traction mechanism is constructed by combining pollution response curves to guide streamline shift and extend pollution migration paths; Establish a control response limit table based on pollutant concentration fluctuations, fine-tune the control input parameters in real time, and complete the limit write-back and risk closed-loop control of the pollution release path.

[0015] A regional groundwater environmental risk assessment system considering the vulnerability of risk receptors includes a pollution fingerprint identification module, a disturbance response analysis module, a release path inversion module, a feedforward regulation generation module, and a streamline limiting closed-loop module. The pollution fingerprint identification module establishes a cross-timescale monitoring baseline, collects groundwater recharge intensity sequences and hydrochemical sequences, extracts sedimentary pollution fingerprints through multidimensional correlation, and identifies germination areas with re-release characteristics. The disturbance response analysis module calculates the temporal-frequency coherence and phase residual between the replenishment disturbance sequence and the deposited contamination fingerprint under the constraint of cross-timescale monitoring baseline, identifies the key spatiotemporal locations that trigger contamination release, and constructs the contamination release causal residual spectrum. The release path inversion module constructs a phase conjugate inversion chain based on the causal residual spectrum, combines historical data of replenishment disturbances to perform counterfactual playback, reconstructs the pollution release path, and forms a transient release threshold surface and precursor score. The feedforward control generation module, based on the transient release threshold surface and precursor score, injects a feedforward constraint field into the high-risk streamline region to regulate the extraction intensity, replenishment rhythm and boundary water level, and generate pre-suppression control parameters. The streamline limiting closed-loop module, under the action of pre-suppression control parameters, initiates the time-reversal conformal mechanism, applies frequency misalignment traction and time-scale traction to adjust the groundwater streamline, implements amplitude limiting backwriting of the release path, and realizes pollution pulse reduction and risk closed-loop control.

[0016] The technical effects and advantages provided by the present invention in the above technical solution are as follows: This invention extracts sedimentary pollution fingerprints through cross-timescale data fusion, accurately identifies easily released areas, and then uses time-frequency coherence and phase residuals to dynamically reveal the causal relationship between replenishment disturbances and pollution release. Based on this, it combines phase conjugate inversion and counterfactual playback mechanisms to reconstruct the release path and form a quantifiable threshold surface and precursor scoring system. Finally, through frequency misalignment and time inversion strategies, it achieves real-time amplitude-limited rewriting and closed-loop control of the pollution release process. Compared to traditional risk assessment methods based on steady-state simulation or concentration exceedance, this invention significantly improves the sensitivity and timeliness of response to sudden pollution behavior, and has significant application value and engineering prospects in precise groundwater pollution prevention and control, risk zoning and priority management, and sustainable water resource utilization. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this invention. For those skilled in the art, other drawings can be obtained based on these drawings.

[0018] Figure 1 The flowchart of the method for assessing regional groundwater environmental risk considering the vulnerability of risk receptors in this invention is shown.

[0019] Figure 2 This is a schematic diagram of the modules of the regional groundwater environmental risk assessment system that considers the vulnerability of risk receptors in this invention. Detailed Implementation

[0020] Exemplary embodiments will now be described more fully with reference to the accompanying drawings. However, these exemplary embodiments can be implemented in many forms and should not be construed as limited to the examples set forth herein; rather, they are provided so that the description of this disclosure will be more complete and fully convey the concept of the exemplary embodiments to those skilled in the art.

[0021] This invention provides, for example Figure 1 The method for assessing regional groundwater environmental risk considering the vulnerability of risk receptors, as shown, includes the following steps: S100, establish cross-timescale monitoring baseline, including collecting groundwater recharge intensity sequence and hydrochemical sequence, extracting pollutant deposition characteristics in groundwater sedimentary layers based on multidimensional correlation analysis, forming a sedimentary pollution fingerprint, and identifying germination areas with pollutant re-release characteristics; To identify the risk of re-release of contaminants from groundwater sedimentary layers, establishing accurate cross-timescale monitoring baselines is a crucial prerequisite. This process involves multiple stages, from data acquisition to contamination fingerprinting and risk zone delineation. The specific operational steps are detailed below, taking into account the actual implementation environment:

[0022] A network of observation wells will be deployed in the study area sensitive to groundwater pollution, employing a multi-layered and multi-depth distribution strategy to cover major aquifers and potential recharge pathways. Equipment will be deployed to continuously monitor groundwater recharge intensity and collect hydrochemical parameters in the target area. Recharge intensity monitoring mainly includes recording atmospheric precipitation, surface water infiltration rate, groundwater level dynamics, and artificial recharge activities, differentiating recharge processes based on daily and hourly scales. Hydrochemical parameter collection includes, but is not limited to, total dissolved solids, conductivity, redox potential, pH, concentrations of major anions and cations (such as nitrates, chloride ions, sodium ions, and calcium ions), concentrations of organic pollutants, and heavy metal content in groundwater. The sampling period will be consistent with the recharge intensity monitoring to ensure comparability across time scales. To ensure data representativeness and continuity, all sampling points should be equipped with automatic data recording devices and remote transmission units to ensure timeliness and completeness. The data collection period should be no less than three typical recharge cycles to fully cover seasonal recharge variations and the impact of extreme events.

[0023] After obtaining the recharge intensity sequence and hydrochemical parameter sequence, spatial normalization was performed on data from different monitoring wells to eliminate background biases caused by differences in well location, depth, and formation structure. Subsequently, the correlation between recharge intensity changes and hydrochemical parameter responses was analyzed point-by-point based on the synergistic trend of these indicators. Unlike existing technologies that only perform single-point time-series analysis of hydrochemical parameters, this invention introduces a cross-timescale analysis method to interpret hourly recharge fluctuations in conjunction with daily, weekly, and even seasonal hydrochemical responses. The analysis focuses on the following three types of correlation patterns among indicators: the first is the coupling change between redox potential and metal ion concentration, used to reveal the redox release process of deposited pollutants; the second is the response relationship between pH fluctuations and organic pollutant concentration, used to determine the release stability of organic pollutants under recharge disturbances; and the third is the synchronous response characteristics of conductivity and major inorganic ions, used to indicate the driving force of recharge on the overall ion migration of the sedimentary layer. Based on these analyses, a water chemical evolution trajectory map of each monitoring point under different replenishment events was constructed and superimposed with the corresponding replenishment intensity sequence. Its time response lag characteristics, response intensity peak, response duration and other features were extracted to form a traceable replenishment response feature database.

[0024] Based on a recharge response characteristic database, and focusing on the synergistic pattern between recharge intensity and hydrochemical response, pollutant deposition characteristics were extracted well by well to construct a depositional pollution fingerprint map. During the construction process, regions exhibiting significant nonlinear responses—time intervals where recharge intensity changes relatively little while hydrochemical parameters fluctuate dramatically—were prioritized for identification. These time intervals typically correspond to pollutants in a physically or chemically unstable state within groundwater sedimentary layers. A depositional pollution response weight matrix was formed by comprehensively evaluating the frequency of occurrence, the number of response parameter types, and peak amplitudes of these regions, and this matrix was used to classify the depositional modes of various pollutants. For example, the response characteristics of complexed heavy metals are characterized by a sharp decrease in redox potential and a simultaneous increase in metal ion concentration, while adsorbed organic pollutants exhibit a short-term release peak after a sudden pH change. These typical response characteristic parameters were combined and categorized to form a three-dimensional identifier map of pollutant type, binding mode, and release behavior, thus constructing a characteristic "fingerprint" of groundwater depositional pollution in the region.

[0025] After constructing the sedimentary pollution fingerprint, spatial inversion and risk localization were performed on the fingerprint information of the entire monitoring area to identify germination areas with the risk of pollutant re-release. This process combined multi-source information such as geological structure, hydrological boundaries, groundwater flow direction, and recharge pathways to spatially map the sedimentary pollution fingerprint and perform cluster analysis on high-intensity response areas. The clustering criteria included not only the frequency and intensity of fingerprint features but also the degree of consistency with changes in recharge intensity and the persistence of hydrochemical response peaks. Based on the clustering results, the distribution maps of main groundwater flow channels and concentrated recharge areas were overlaid to ultimately identify high-risk locations where multiple pollutants exhibit re-release behavior and delineate them as germination areas. These areas, as spatial triggering units for pollution release risk, play a core anchoring role in subsequent disturbance response analysis and pre-control strategy formulation.

[0026] S200, under the constraint of cross-timescale monitoring baseline, calculates the temporal and frequency coherence and phase residual between groundwater recharge disturbance sequence and sedimentary pollution fingerprint, analyzes the dynamic impact of groundwater recharge disturbance on sedimentary layer stability, identifies key spatiotemporal locations that trigger pollutant release, and constructs a causal residual spectrum of pollution release. To clarify the dynamic driving mechanism of groundwater recharge disturbance on the re-release of deposited pollutants, under the constraint of the established cross-timescale monitoring baseline, it is necessary to further quantify the dynamic relationship between recharge disturbance and deposited pollution, identify its key response mechanisms, and construct the spatiotemporal causal characteristics of pollution release accordingly. The specific implementation method of this process is described in detail below, combined with a real-world environment:

[0027] Based on the established sedimentary contamination fingerprint, the recharge intensity sequence and hydrochemical response sequence obtained from each observation well are synchronously paired to establish a time-to-time correspondence recharge-response data structure. The recharge intensity sequence includes changes in infiltration during heavy rainfall, the dynamic rise rate of groundwater level, and the frequency and duration of external water input; while the hydrochemical response sequence includes measured data reflecting changes in pollutant state, such as pollutant concentration, redox potential, conductivity, and total organic carbon content. Based on this, the recharge intensity sequence and pollutant response sequence undergo dual time-scale alignment. Specifically, through linear interpolation, sample interpolation, and time window sliding, data from different sampling frequencies are uniformly reconstructed to ensure that each time point contains accurate recharge input information and corresponding hydrochemical response information. The innovation of this step lies in the introduction of cross-scale integration processing, enabling daily, hourly, and minute-level perturbation signals to simultaneously participate in subsequent coherent analysis, overcoming the limitation of existing technologies that are limited to single-time-granularity modeling.

[0028] After aligning the supply-response data, a representative supply disturbance process is selected as the research window to conduct time-frequency linkage analysis of the disturbance-driven characteristics. This analysis focuses on identifying the lag response duration, peak response time, persistence, and periodic oscillation patterns of the pollution response sequence to the supply disturbance. Unlike existing technologies that only perform intensity analysis on the pollution response, this invention introduces a coherence analysis mechanism to quantify the linkage strength between the disturbance input and the pollution response in a high-temporal-resolution manner. Specifically, the supply disturbance sequence is used as the driving sequence, and the pollution response sequence as the driven sequence. Wavelet packet transform is used to expand them across multiple scales, extracting their common frequency components and coherence intensities in each frequency band. Simultaneously, by combining the fluctuations of the coherence spectrum in each time window, the periods when the pollution response intensity increases or decreases with the supply disturbance are identified, providing a basis for subsequent identification of key response intervals.

[0029] Based on the peak region of the coherent spectrum, a phase residual analysis mechanism is further introduced to characterize the propagation delay between the supply disturbance and the pollution response. The phase difference change between each pair of supply disturbance peaks and pollution response peaks is measured, and a continuous-time phase evolution sequence is established. Particular attention is paid to three scenarios after the supply disturbance: premature rise, synchronous enhancement, or delayed release of the pollution response, corresponding to different types of pollutant release mechanisms. For example, complexed heavy metals typically exhibit a delayed release response, while free organic matter may respond synchronously with the supply. Through the cumulative calculation of the phase residuals, media hindrance points or fast channels in the disturbance propagation path are further identified.

[0030] By combining the results of coherence spectroscopy and phase residuals, a causal residual spectrum of pollution release is constructed to comprehensively characterize the dynamic causal structure between disturbance input and pollution response. In implementation, firstly, the characteristics of replenishment disturbance fluctuations under different time windows are causally mapped to the corresponding hydrochemical responses, identifying the frequency, intensity, and spatial distribution patterns of pollution release events triggered by disturbances. Subsequently, by superimposing causal event paths from different observation points, a spatiotemporal projection map of pollution release is generated, clarifying the propagation trajectory of the pollution response triggered by each disturbance input along the time and spatial dimensions. The causal residual spectrum constructed based on this can be used to identify key disturbance events, anomalous response nodes, and potential risk triggering sources.

[0031] S300 constructs a phase conjugate inversion chain based on causal residual spectrum, combines historical groundwater recharge disturbance data for counterfactual playback, reconstructs the pollutant release process and release path, and forms a pollutant transient release threshold surface and release precursor scoring index system. Building upon the completed causal residual spectrum of pollution release, further inversion analysis is needed to reconstruct the entire pollution release process and establish a quantitative threshold and early warning indicator system for prediction and control, in order to achieve in-depth analysis and pathway tracing of pollutant release mechanisms. The following section details the implementation of this process, based on the results of the previous stage:

[0032] Based on the obtained causal residual spectrum of pollution release, key disturbance events and their corresponding pollution response sequences are extracted and reverse-ordered in the time dimension to construct a disturbance-response reverse propagation trajectory. To this end, each typical disturbance event (such as short-duration heavy rainfall, sudden groundwater level rise, artificial groundwater recharge, etc.) needs to be numbered and classified, and its pollution response path map formed in the causal residual spectrum needs to be extracted. Each path map contains multiple blocks, including the disturbance initiation point, response time delay segment, pollution response peak area, and recovery plateau segment. By structurally encoding the path maps corresponding to each disturbance event, a disturbance event database is generated, and the pollution response is traced back from the endpoint to the initiation point in reverse time order, initially forming the core framework of the inversion chain. During this process, the changing trend of pollution response intensity in the tracing path needs to be specifically marked, including the decrease in response velocity, the widening or convergence of the response peak, and the disappearance of the fluctuation cycle. This operation breaks through the path dependency of existing technologies based on forward evolution for process reconstruction, and has the independent reverse modeling capability to deduce driving factors from pollution results.

[0033] Based on the constructed disturbance-response reverse path, a phase conjugate processing mechanism is introduced to perform time-symmetric mapping on the pollution response sequence to simulate the ideal state of the groundwater system when no recharge disturbance occurs. In specific implementation, the time series of the pollution response is mirrored around the disturbance initiation point, generating a virtual response trajectory on the conjugate time axis. This trajectory represents the natural distribution state of pollutants in the groundwater system under undisturbed conditions, exhibiting strong stability, small fluctuation amplitude, and weak response delay. By comparing the conjugate trajectory with the actual disturbance response trajectory point by point, the abnormal growth range of the pollution response amplitude, response direction shift, and response rate caused by the disturbance event can be identified. Combining the overlapping areas of the reverse path and the conjugate trajectory, sensitive segments in the pollutant release process are further extracted, and key nodes with the smallest disturbance intensity but the largest pollution response are marked.

[0034] Based on a disturbance-response comparison dataset obtained through counterfactual playback, a refined reconstruction of pollution release pathways was conducted. This reconstruction, using time series as the main thread, combined with groundwater recharge disturbance characteristics, pollution response characteristics, and geological-hydrological boundary conditions, performed a spatial-temporal multidimensional mapping of the pollutant release process. In practice, the disturbance pathway, response pathway, and conjugate trajectory were projected onto the groundwater flow profile to construct a three-dimensional groundwater release channel map. This channel map reveals the path, time window, and force direction of pollutant migration from sedimentary layers to the water body, while also clarifying the displacement trend of groundwater streamlines under disturbance and the enrichment areas of pollutants in the aquifer. Through overlay analysis of release pathways from multiple disturbance events, common characteristic segments were extracted and categorized into main pollutant release channel areas, disturbance-induced areas, and sediment re-release sensitive areas.

[0035] After reconstructing the pollution release pathways, a threshold surface for transient pollutant release and a scoring index system for precursors were developed based on response characteristic data from various disturbance events. This step first establishes threshold classification standards according to different pollutant types, release mechanisms, and response intensities. For example, for heavy metal pollutants, the ratio of the amplitude of the redox potential jump to the increase in metal ion concentration is set as the response threshold; for organic pollutants, the rate of change of total organic carbon concentration within a specific time window is used as the judgment indicator. The data of these indicators under different disturbance scenarios are normalized and their spatial distribution is superimposed to form a threshold response surface. This threshold surface reveals the critical conditions required for pollutants to transition from a stable state to a migratory state under different recharge intensities, response lag times, and pollution concentrations. Based on this, a set of precursor scoring factors is constructed by statistically analyzing the micro-changes in hydrochemistry before multiple disturbance events, including the frequency of short-term water level fluctuations, the slope of slight pH shifts, and the rate of increase in conductivity fluctuation amplitude, and weights are assigned to each factor to form a comprehensive scoring system. This scoring system allows for the early identification of pollution release risk trends in actual monitoring, providing a predictive basis for subsequent proactive intervention and dynamic control.

[0036] S400, based on transient release threshold surface and precursor score results, sets feedforward constraint field for high-risk groundwater flow line areas, regulates groundwater extraction intensity, recharge rhythm and boundary water level, and generates pre-suppression control parameters to reduce the trend of pollutant re-release in advance. Building upon the completion of the construction of the transient release threshold surface for pollutants and the extraction of the pre-release precursor scoring index system, a feedforward control mechanism based on real-time risk status is needed to proactively intervene in and suppress the re-release trend of pollution in high-risk areas. The following section, building upon the results of the previous stage, details the specific implementation method for generating pre-suppression control parameters:

[0037] Based on the constructed transient release threshold surface and precursor scoring index system, the entire groundwater flow field is spatially projected. The risk score results of each monitoring point are superimposed on the regional groundwater flow direction map. Combined with geological profiles, hydraulic gradient distribution, and recharge boundary conditions, concentrated connected sections with risk scores significantly higher than the average threshold are identified. These sections are typically characterized by distribution along the main groundwater flow lines, being affected by multiple recharge disturbances, and having fault zones, fissure zones, or highly permeable sand layers in the geological structure. After identifying high-risk groundwater flow line areas, a feedforward control range is set with these areas as the core, and three-dimensional regional units are divided. Unlike existing technologies that typically use only pollution concentration distribution as the basis for control, this implementation method determines flow lines based on dynamic risk score results, which can more effectively identify potentially highly sensitive areas where pollutant release is imminent.

[0038] Within high-risk groundwater flow lines, a responsive simulation analysis of extraction behavior to pollutant release behavior is conducted. Specifically, based on historical groundwater level data, extraction well operation records, aquifer storage parameters, and flow boundary conditions, a groundwater flow disturbance-pollution response curve is constructed. The focus is on analyzing the triggering time and duration of the pollution response peak under different extraction intensities. By comparing changes in pollutant concentration, response lag period, and threshold breakthrough frequency in the release path before and after extraction, it is determined whether the current extraction behavior has a tendency to induce pollution release. If extraction-induced release behavior exists, the extraction control strategy needs to be restructured, including reducing extraction intensity, extending the extraction cycle, and adjusting the extraction time period. Furthermore, for areas with a composite distribution of deep and shallow aquifers, a vertical extraction regulation strategy should be implemented to avoid disturbance of the upper pollution migration zone caused by deep extraction. In this process, the hydraulic radius, recharge capacity, and spatial intersection with the pollution release path of each extraction well need to be clearly defined to form a precise and controllable extraction parameter database.

[0039] To explore the relationship between recharge rhythm and pollutant response, a recharge input regulation model based on precursor scores was constructed. The previously established precursor score index system identified short-term high-frequency recharge, sudden water replenishment events, or periodic high-fluctuation processes as having a significant inducing effect on pollutant release. Based on this, upper limits for recharge intensity adjustment, lower limits for time intervals, and cumulative fluctuation thresholds were set according to different precursor score intervals. For areas with a strong risk of recharge-induced pollution, the disturbance intensity of high-frequency recharge input on the depositional pollutant release pathway was reduced by limiting the duration of artificial recharge, dispersing the infiltration rate of natural precipitation, and establishing retention buffer zones. Furthermore, the hydraulic response characteristics of the recharge inlet boundary were comprehensively considered, and a water level response curve was constructed to analyze whether the groundwater level rise rate under a specific recharge input could exceed the critical point defined in the release threshold curve, thereby further guiding the fine-tuning design of the recharge rhythm.

[0040] After optimizing the extraction intensity and recharge rhythm, a boundary water level control strategy is constructed, centered on the aforementioned high-risk areas and considering the groundwater level boundary conditions, to further enhance the ability to suppress pollution release trends. Boundary water level, as the fundamental driving factor for groundwater flow direction and kinetic energy, often has a nonlinear impact on the release path due to even minor changes. In practice, the water level change trends at key locations along the regional boundary (such as upstream infiltration zones, downstream outflow points, and areas near regional geological barriers) should be monitored, and the direction of boundary conditions driving release behavior should be determined by combining the release threshold surface trend and streamline offset trend. Based on this, dynamic fine-tuning of the boundary water level should be implemented, for example, appropriately raising the downstream boundary water level during high-risk pollution release periods to suppress flow velocity, or lowering the upstream water level during the release response mitigation phase to reduce infiltration kinetic energy. Furthermore, for sites with lateral recharge characteristics, spatial shaping control of the boundary hydraulic gradient can be achieved through engineering measures such as constructing underground impermeable curtains and seepage regulation zones, thereby forming a three-dimensional, multi-boundary linkage control pattern. Finally, the data results from the three control dimensions of extraction intensity, replenishment rhythm, and boundary water level are linked and integrated to generate a set of pre-suppression control parameters for dynamic input, which can then be invoked and executed by the subsequent closed-loop regulation mechanism for pollutant release.

[0041] Under the action of pre-suppression control parameters, S500 initiates the time-reversal conformal regulation mechanism, applies the golden ratio frequency misalignment traction and dual mirror time scale traction, adjusts the groundwater flow line distribution in real time, implements amplitude-limited write-back control on the pollution release path, and completes the dynamic reduction of pollution pulses and risk closed-loop regulation. Based on the generation of pre-suppression control parameters, in order to achieve precise intervention and real-time response regulation of pollution release channels, it is necessary to further construct a groundwater flowline regulation mechanism based on a combination of time inversion and traction strategies. The following describes the implementation method of the entire process in conjunction with the aforementioned steps: With feedforward control parameters set and injected into the target groundwater flow field, a time-reversal driving field is constructed based on real-time monitoring data and historical pollution response data to reshape groundwater flow trends and deduce the actual migration paths of pollutants during previous disturbances. The key to this process lies in reversing the time series of recharge intensity, extraction rhythm, and boundary water level changes, and combining this with the dynamic fluctuation trends of pollutant concentration response curves to achieve the reverse evolution of the pollution release process. Unlike traditional forward modeling, this implementation does not focus on predicting future changes, but rather on the historical trajectory of pollution occurrences. Through time-reversal analysis, it identifies the dynamic input patterns leading to release behavior and their evolutionary mechanisms over time. In the inversion driving field, it is crucial to identify the "backward convergence" points of pollutant concentrations at different well locations—key nodes where multiple monitoring points show a consistent trend in pollutant concentration changes during time-reversal. These nodes typically represent high-flux transport channels or concentrated disturbance areas along the actual pollutant release path. Based on this inversion identification result, the dominant flow lines and convergence areas requiring regulatory intervention are further clarified.

[0042] Based on the construction of the time-reversal driving field, a golden ratio frequency misalignment strategy is introduced to dynamically and rhythmically disturb the groundwater flowline, thereby breaking the stable excitation mode formed by the existing disturbance on the pollutant release path. This strategy originates from the simulation of the asynchronous characteristics of disturbances under natural rhythms. That is, without disrupting the overall replenishment balance of the system, a slightly deviated difference in replenishment and extraction cycles is artificially created to achieve a "golden misalignment" of the disturbance frequency. In specific operations, the operating frequency of groundwater extraction and replenishment equipment is set so that the ratio of their operating cycles is close to the golden ratio (approximately 0.618:1), thereby preventing the continuous excitation of pollutant release behavior due to periodic resonance. In this process, the time nodes of the original release path also need to be mapped onto the frequency misalignment cycle to form a time axis offset table, which is used as the basis for adjusting the rhythmicity of the flowline. The core of this operation is not to directly intervene in the pollutants themselves, but to achieve soft interference with the conditions for pollutant migration by controlling the fluid dynamic conditions, thereby indirectly suppressing the further amplification of pollutant release flux.

[0043] While implementing a frequency-shifting traction strategy, a dual-mirror timescale traction mechanism is constructed to achieve directional reshaping of the pollution release path and dispersion of disturbance energy. The dual-mirror timescale refers to mirroring the pollution response curves before and after pollution release on the time axis, centered on the current observation time, thus forming two mimicking traction trajectories: time reversal and time extension. In specific implementation, the acceleration phase before the pollution release peak and the decay phase after the peak need to be mirrored and paired, and this should be used as the control basis to introduce directional offset design in the groundwater flowline distribution. For example, in the area where the pollution response acceleration is most intense in the early stage of the peak, a hydraulic slow-release section should be set up on the flowline to increase the length of the pollutant movement path or extend the migration time; while in the area where the pollution response rapidly subsides in the later stage of the peak, a "traction convergence" zone for pollutant residues should be formed by slightly increasing the recharge intensity or appropriately lowering the boundary water level, thereby hindering its further diffusion. By using dual-mirror timescales, we can not only actively extend and shift the flow path of pollutants, but also enhance the sensitivity of pollutant migration trends to regulatory behavior by leveraging the laws of time symmetry, thus significantly improving the timeliness and effectiveness of regulatory interventions.

[0044] After constructing and deploying the time-reversal driving field, frequency misalignment traction strategy, and dual-mirror timescale traction mechanism, the next stage is the limit-based write-back control of the pollution release path. This aims to achieve gradual reduction of pollution pulses and closed-loop risk stabilization through dynamic feedback of the control structure. In practice, dynamic monitoring of the concentration fluctuation range, release time window length, and hydraulic gradient changes before and after release along the existing release path is required. Based on this, a control response limit table is established. This table clarifies the acceptable range within which pollution release behavior must be suppressed under different control strategies. For example, in a high-risk zone, concentration fluctuations must not exceed the set maximum fluctuation rate, and the release cycle must not exceed the preset window value. Based on this, and according to real-time monitoring results, the control input parameters are written back in small increments. This involves fine-tuning the extraction rhythm, boundary water level, or recharge intensity based on the current pollution response results to form a feedback loop. All write-back actions must be performed without disrupting the overall stability of the flow field and to achieve maximum control effect with minimal energy consumption. Ultimately, through multiple rounds of limited write-back operations, the pollutant release pathways were stabilized and tended to be inactive, concentration fluctuations were controlled below the tolerable threshold, and the groundwater flow line was restored to a sustainable management level, thus completing the overall reduction and dynamic control loop of the pollution pulse.

[0045] This invention extracts sedimentary pollution fingerprints through cross-timescale data fusion, accurately identifies easily released areas, and then uses time-frequency coherence and phase residuals to dynamically reveal the causal relationship between replenishment disturbances and pollution release. Based on this, it combines phase conjugate inversion and counterfactual playback mechanisms to reconstruct the release path and form a quantifiable threshold surface and precursor scoring system. Finally, through frequency misalignment and time inversion strategies, it achieves real-time amplitude-limited rewriting and closed-loop control of the pollution release process. Compared to traditional risk assessment methods based on steady-state simulation or concentration exceedance, this invention significantly improves the sensitivity and timeliness of response to sudden pollution behavior, and has significant application value and engineering prospects in precise groundwater pollution prevention and control, risk zoning and priority management, and sustainable water resource utilization.

[0046] This invention provides, for example Figure 2 The regional groundwater environmental risk assessment system shown includes a pollution fingerprint identification module, a disturbance response analysis module, a release path inversion module, a feedforward regulation generation module, and a streamline limiting closed-loop module. The pollution fingerprint identification module establishes a cross-timescale monitoring baseline, collects groundwater recharge intensity sequences and hydrochemical sequences, extracts sedimentary pollution fingerprints through multidimensional correlation, and identifies germination areas with re-release characteristics. The disturbance response analysis module calculates the temporal-frequency coherence and phase residual between the replenishment disturbance sequence and the deposited contamination fingerprint under the constraint of cross-timescale monitoring baseline, identifies the key spatiotemporal locations that trigger contamination release, and constructs the contamination release causal residual spectrum. The release path inversion module constructs a phase conjugate inversion chain based on the causal residual spectrum, combines historical data of replenishment disturbances to perform counterfactual playback, reconstructs the pollution release path, and forms a transient release threshold surface and precursor score. The feedforward control generation module, based on the transient release threshold surface and precursor score, injects a feedforward constraint field into the high-risk streamline region to regulate the extraction intensity, replenishment rhythm and boundary water level, and generate pre-suppression control parameters. The streamline limiting closed-loop module, under the action of pre-suppression control parameters, initiates the time-reversal conformal mechanism, applies frequency misalignment traction and time-scale traction to adjust the groundwater streamline, implements amplitude limiting backwriting of the release path, and realizes pollution pulse reduction and risk closed-loop control.

[0047] The regional groundwater environmental risk assessment method considering risk receptor vulnerability provided in this embodiment of the invention is implemented through the aforementioned regional groundwater environmental risk assessment system considering risk receptor vulnerability. For details of the specific methods and procedures of the regional groundwater environmental risk assessment system considering risk receptor vulnerability, please refer to the embodiments of the aforementioned regional groundwater environmental risk assessment method considering risk receptor vulnerability, which will not be repeated here.

[0048] The foregoing has only described certain exemplary embodiments of the present invention by way of illustration. Undoubtedly, those skilled in the art can modify the described embodiments in various ways without departing from the spirit and scope of the present invention. Therefore, the foregoing drawings and descriptions are illustrative in nature and should not be construed as limiting the scope of protection of the claims of the present invention.

Claims

1. A method for assessing the environmental risk of groundwater in areas considering the vulnerability of risk receptors, characterized in that, Includes the following steps: S100 establishes cross-timescale monitoring baselines, collects groundwater recharge intensity sequences and hydrochemical sequences, extracts sedimentary pollution fingerprints through multidimensional correlation, and identifies germination areas with re-release characteristics; S200, under the constraint of cross-timescale monitoring baseline, calculates the temporal coherence and phase residual between the replenishment perturbation sequence and the deposited contamination fingerprint, identifies the key spatiotemporal locations that trigger contamination release, and constructs the contamination release causal residual spectrum; S300 constructs a phase conjugate inversion chain based on causal residual spectrum, combines historical data of supply disturbance for counterfactual playback, reconstructs pollution release path, and forms transient release threshold surface and precursor score; S400, based on transient release threshold surface and precursor score, injects feedforward constraint field into high-risk streamline region to regulate extraction intensity, replenishment rhythm and boundary water level, and generates pre-suppression control parameters. S500, under the action of pre-suppression control parameters, initiates the time inversion conformal mechanism, applies frequency misalignment traction and time scale traction, adjusts the groundwater flow line, and implements amplitude-limited rewriting of the release path.

2. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 1, characterized in that, Step S100 includes: In areas sensitive to groundwater pollution, a multi-layered and multi-depth observation well network is deployed to collect recharge intensity sequences and water chemistry sequences. The collected recharge intensity and hydrochemical sequences were normalized to analyze the synergistic relationship between recharge changes and hydrochemical responses, and recharge response features were extracted. Based on the replenishment response characteristics, the characteristics of contamination deposition are extracted to construct a fingerprint map of depositional contamination. By combining sedimentary pollution fingerprint information with groundwater flow direction, geological structure and recharge path for spatial mapping and cluster analysis, high-risk locations for pollutant re-release are identified, forming spatial triggering units for pollution release risk.

3. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 1, characterized in that, Step S200 includes: The recharge intensity sequence obtained from the observation wells is time-aligned with the pollution response sequence to establish a correspondence between recharge and response; Based on the correspondence between supply and response, time-frequency coherence features are extracted to analyze the linkage between supply disturbance and pollution response; Phase residual analysis is introduced based on time-frequency coherence characteristics to characterize the propagation delay behavior of disturbances and identify the response mechanisms of different pollutant releases; By combining coherence characteristics and phase residuals, a causal residual spectrum of pollution release is constructed, forming a dynamic causal structure between perturbation input and pollution response.

4. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 3, characterized in that, When constructing the causal residual spectrum of pollution release, the replenishment disturbance and pollution response paths of different time windows in multiple observation wells are superimposed to form a spatiotemporal projection map of pollution release, which is used to identify the spatial propagation trajectory of key disturbance events and pollution responses.

5. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 1, characterized in that, Step S300 includes: The disturbance events and pollution response sequences are reverse-ordered to form the reverse propagation trajectory of the disturbance response; A phase conjugation processing mechanism is introduced based on the reverse propagation trajectory to construct a time-symmetric mapping trajectory for identifying sensitive areas of pollution release; Based on disturbance path, response path and conjugate trajectory, a three-dimensional groundwater release channel map is constructed to reconstruct the pollution release path; Based on the response characteristics of disturbance events, pollutant release threshold parameters are extracted and a precursor scoring factor system is established to form a quantitative basis for identifying pollution release risks.

6. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 5, characterized in that, When constructing the precursor scoring factor system, the scoring factors include the frequency of water level fluctuations, the slope of pH changes, and the rate of increase in conductivity. By assigning different weights, a comprehensive score is formed to identify the trend of pollution release risk in advance.

7. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 1, characterized in that, Step S400 includes: By overlaying the scoring results with the groundwater flow direction map, high-risk groundwater flow line areas are identified and control units are delineated. In high-risk areas, analyze the driving relationship between extraction behavior and pollution response, adjust extraction intensity and cycle, and form extraction control parameters; A supply rhythm regulation model was constructed by combining the precursor score results, and supply intensity and frequency limits were set to form supply regulation parameters. A water level control strategy is constructed at the boundary of high-risk areas, linking extraction parameters, replenishment parameters, and boundary control results to generate a set of pre-inhibition control parameters for dynamic intervention.

8. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 7, characterized in that, The linkage and integration of extraction control parameters, replenishment control parameters and boundary control results are based on the risk level classification standard set in the pollutant transient release threshold surface, and the extraction intensity adjustment range, replenishment rhythm intervention cycle and boundary water level variation range are set according to different risk levels.

9. The method for regional groundwater environmental risk assessment considering the vulnerability of risk receptors according to claim 1, characterized in that, Step S500 includes: A time-reversal driving field is constructed based on historical pollution response data to identify key transmission channels in the pollution release pathway; Based on the time-reversal driving field, a golden ratio frequency shifting strategy is implemented to regulate the extraction and replenishment rhythm in order to disturb the resonance of the pollution release cycle. A dual-mirror timescale traction mechanism is constructed by combining pollution response curves to guide streamline shift and extend pollution migration paths; Establish a control response limit table based on pollutant concentration fluctuations, fine-tune the control input parameters in real time, and complete the limit write-back and risk closed-loop control of the pollution release path.

10. A regional groundwater environmental risk assessment system considering the vulnerability of risk receptors, used to implement the regional groundwater environmental risk assessment method considering the vulnerability of risk receptors as described in any one of claims 1-9, characterized in that, It includes a contamination fingerprint recognition module, a disturbance response analysis module, a release path inversion module, a feedforward control generation module, and a streamline limiting closed-loop module: The pollution fingerprint identification module establishes a cross-timescale monitoring baseline, collects groundwater recharge intensity sequences and hydrochemical sequences, extracts sedimentary pollution fingerprints through multidimensional correlation, and identifies germination areas with re-release characteristics. The disturbance response analysis module calculates the temporal-frequency coherence and phase residual between the replenishment disturbance sequence and the deposited contamination fingerprint under the constraint of cross-timescale monitoring baseline, identifies the key spatiotemporal locations that trigger contamination release, and constructs the contamination release causal residual spectrum. The release path inversion module constructs a phase conjugate inversion chain based on the causal residual spectrum, combines historical data of replenishment disturbances to perform counterfactual playback, reconstructs the pollution release path, and forms a transient release threshold surface and precursor score. The feedforward control generation module, based on the transient release threshold surface and precursor score, injects a feedforward constraint field into the high-risk streamline region to regulate the extraction intensity, replenishment rhythm and boundary water level, and generate pre-suppression control parameters. The streamline limiting closed-loop module, under the action of pre-suppression control parameters, initiates the time inversion conformal mechanism, applies frequency misalignment traction and time scale traction to adjust the groundwater streamline, and implements amplitude limiting backwriting of the release path.

Images