Intelligent deduction method for spatial occurrence state of karst aquifer facing uncertainty propagation

Abstract

The application discloses a karst aquifer space occurrence state intelligent deduction method for uncertainty propagation, relates to the technical field of aquifer space occurrence state deduction, and comprises the following steps: collecting three-dimensional tomographic scanning response information corresponding to different buried depth positions during continuous rainfall, recording wave velocity offset amplitude, water pressure fluctuation state and seepage transfer track corresponding to the time advancing direction in the fissure area, and generating a time-varying distribution map corresponding to the underground water migration state; extracting an abnormal drift area corresponding to adjacent moments in the time-varying distribution map, and tracking the wave velocity jump position along the underground water migration direction. The application continuously tracks the underground abnormal propagation state, dynamically correlates the wave velocity offset, water pressure change and seepage track, realizes continuous identification of the underground abnormal propagation path and the hidden water-conducting migration state, and thus improves the dynamic deduction accuracy of the underground water-rich boundary and reduces the boundary distortion influence.

Description

Technical Field

[0001] This application relates to the field of spatial occurrence state inference technology of aquifers, specifically to an intelligent inference method for spatial occurrence state of karst aquifers oriented towards uncertainty propagation. Background Technology

[0002] The spatial occurrence state of karst aquifers refers to the distribution, connectivity, water abundance, extension direction of dissolution channels, and coupling occurrence characteristics between fissures, caves, and underground rivers in the three-dimensional geological environment of the aquifer space in the underground karst medium. It not only reflects the geometric structure of the aquifer medium in the spatial scale, but also reflects the dynamic retention and migration laws of groundwater under different strata, different burial depths, and different hydraulic conditions.

[0003] The intelligent simulation process for the spatial occurrence state of karst aquifers in response to the propagation of uncertainty refers to the continuous multi-angle and multi-scale detection of the underground medium using three-dimensional tomography technology. This process acquires wave velocity response, electrical differences, density changes, and water-bearing anomalies at different depths, forming an iteratively updatable underground structural tomographic map in three-dimensional space. Subsequently, combined with borehole exposure information, hydrological observation data, and historical water inrush evolution trajectories, spatial correlation propagation analysis is performed on the anomalous areas in the tomographic map. During the simulation process, the diffusion paths of uncertainties caused by detection errors, medium heterogeneity, and local data gaps are simultaneously tracked, allowing various uncertain information to be continuously transmitted along the underground karst network. Based on the response correlation strength between different regions, the spatial connectivity between hidden caves, fissure channels, and water-rich areas is dynamically reconstructed, thereby gradually forming a three-dimensional occurrence result of karst aquifers with probability distribution characteristics. This enables intelligent prediction of underground water-rich structures, water-conducting channels, and high-risk water inrush areas.

[0004] The existing technology has the following shortcomings: In the intelligent simulation of the spatial occurrence state of karst aquifers facing uncertainty propagation, when the water flow inside karst fissures is continuously replenished by rainfall, a local pulse-like enhancement of flow velocity will occur. The direction of groundwater migration will rapidly deflect in a short period of time, and the original hydraulic propagation rhythm inside the fissures will be suddenly disrupted, resulting in a dynamic shift in the abnormal response trajectory obtained by three-dimensional tomography. As the abnormal response continues to propagate, the uncertainty propagation path will undergo abrupt changes in direction during the simulation process, which will rapidly activate areas that were originally in a low-risk state and gradually form hidden water-conducting channels. Ultimately, this will lead to a significant distortion in the determination of groundwater-rich boundaries, seriously affecting the accuracy of the simulation of water inrush hazard areas.

[0005] 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

[0006] The purpose of this application is to provide an intelligent method for extrapolating the spatial occurrence state of karst aquifers in the face of uncertainty propagation, so as to solve the problems mentioned in the background art.

[0007] To achieve the above objectives, this application provides the following technical solution: an intelligent method for extrapolating the spatial occurrence state of karst aquifers oriented towards uncertainty propagation, comprising the following steps: Three-dimensional tomographic scan response information at different burial depths during continuous rainfall was collected. Wave velocity shift amplitude, water pressure fluctuation state and seepage transfer trajectory were recorded along the time progression direction within the fracture area, and a time-varying distribution map of the corresponding groundwater migration pattern was generated. Extract the abnormal drift regions corresponding to adjacent moments in the time-varying distribution map, track the wave velocity jump positions along the groundwater migration direction, and spatially arrange the continuously expanding offset regions according to the order of advancement to form a dynamic evolution region corresponding to the pulse seepage changes; A progressive mapping is performed on the hydraulic disturbance range corresponding to each moment in the dynamic evolution area. The deflection positions that repeatedly appear in a short period of time are continuously projected to the deep region, and the diffusion amplitude of the abnormal drift direction at different time positions is recorded to generate an underground anomaly migration direction map. Read the deflection areas that are continuously expanding within the underground anomaly migration direction map, perform synchronous compression on the diffusion trajectories corresponding to multiple time locations, and retain the migration direction information corresponding to the areas with a sudden increase in diffusion speed during the compression process to obtain the cluster areas corresponding to the hidden water-rich migration state. Extract the continuously expanding migration boundary within the aggregation area, perform gradual traction on the boundary change trend along the time progression, and adjust the spatial orientation of the underground water-rich area during the traction process to complete the dynamic simulation of the spatial occurrence state of the karst aquifer.

[0008] Preferably, the time-varying distribution map of the corresponding groundwater migration pattern includes: Cyclic detection was performed on the scanning areas corresponding to different burial depths, and the spatial response status at each time location was recorded to construct a set of underground dynamic response information. Extract the wave velocity offset amplitude, water pressure fluctuation state and seepage transfer trajectory corresponding to the underground dynamic response information set, and spatially connect the wave velocity change location, water pressure propagation path and seepage movement trajectory to obtain dynamic response correlation information; The wave velocity offset region, water pressure propagation path and seepage expansion path corresponding to the dynamic response association information are continuously expanded to obtain the continuous evolution state of underground dynamic changes. The results of underground response changes corresponding to the continuous evolution of underground dynamic changes are uniformly integrated to form a time-varying distribution map of the corresponding groundwater migration pattern.

[0009] Preferably, spatially connecting the locations of wave velocity changes, water pressure propagation paths, and seepage movement trajectories includes: longitudinally connecting the extension paths formed at different burial depths in the same water pressure propagation direction, and laterally connecting the water pressure change areas corresponding to multiple time locations.

[0010] Preferably, the dynamic evolution region forming the corresponding pulse seepage change includes: Extract underground anomaly regions that shift in position at adjacent times within the time-varying distribution map, and synchronously record the spatial shift direction corresponding to the anomaly drift regions to obtain a set of anomaly drift regions; Based on the set of abnormal drift regions, the wave speed jump positions formed by adjacent time positions are continuously extracted, and the spatial movement direction corresponding to the wave speed jump position is expanded segment by segment to form a continuous movement path. Based on the continuous movement path, the continuously expanding offset region inside the abnormal drift region is unfolded time by time, and the positions where the offset region expands and changes are spatially connected to construct a continuous advancement zone. By unifying and integrating the spatial migration paths, wave velocity jump paths, and migration extension regions formed at multiple time locations around the continuous advancement zone, a dynamic evolution region corresponding to the pulse seepage changes is formed.

[0011] Preferably, continuously extracting wave velocity jump positions formed by adjacent time positions includes: recording positions where the wave velocity change amplitude exceeds the continuous change range of the surrounding area as wave velocity jump positions, and continuously connecting the spatial expansion range formed by the same wave velocity jump area at multiple time positions.

[0012] Preferably, generating a subsurface anomaly migration directional map includes: The wave velocity change region, water pressure propagation region, and seepage expansion region corresponding to different time locations within the dynamic evolution region are obtained, and the underground anomaly change boundary is continuously extended and spatially connected to obtain the progressive mapping region. Under the progressive mapping region, the positions where the direction changes repeatedly occur within the same spatial range are recorded, and the deflection path is extended and connected vertically layer by layer to construct a deep continuous projection path; Based on the deep continuous projection path, the diffusion range corresponding to the anomaly drift direction is continuously extended, and the diffusion amplitude formed at multiple time locations is uniformly arranged to generate an underground anomaly migration direction map.

[0013] Preferably, recording locations where directional changes occur repeatedly within the same spatial range includes: comparing underground anomaly areas corresponding to multiple adjacent time locations segment by segment, recording locations where directional changes occur repeatedly within the same spatial range as repeated deflection locations, and spatially connecting the deflection direction corresponding to the previous time location with the deflection direction corresponding to the next time location.

[0014] Preferably, the clustered areas corresponding to the latent water-rich migration state include: Read the deflection area that continuously expands along the direction of groundwater migration in the map, and expand and spatially connect the underground anomaly expansion boundary corresponding to different time locations to output the continuously expanding deflection area. For the continuously expanding deflection area, extract the underground anomaly diffusion trajectory corresponding to multiple time locations, and perform gradual convergence and edge contraction on the underground anomaly diffusion trajectories with similar propagation directions to obtain the synchronously compressed trajectory area. The system records the advance distance corresponding to the underground anomaly expansion boundary in the synchronous compression trajectory region, extracts the locations where the advance distance suddenly increases, and spatially connects the underground anomaly propagation directions corresponding to the regions where the diffusion velocity suddenly increases to construct the migration pointing region. The underground anomaly propagation directions corresponding to the migration direction area and the synchronous compression trajectory area are collected, and the underground anomaly accumulation locations are spatially connected and uniformly arranged to obtain the accumulation area corresponding to the hidden water-rich migration state.

[0015] Preferably, performing gradual convergence and edge contraction on underground anomalous diffusion trajectories with similar propagation directions includes: spatially corresponding the underground anomalous diffusion trajectory corresponding to the previous time position with the underground anomalous diffusion trajectory corresponding to the next time position, and performing edge contraction on underground anomalous diffusion trajectories with continuously increasing spatial deviation direction.

[0016] Preferably, the dynamic simulation of the spatial occurrence state of karst aquifers includes: Expand the underground anomaly cluster edges corresponding to different time locations within the cluster area, and perform spatial location correspondence and segment-by-segment connection on the underground anomaly cluster edges to construct a continuous migration boundary; Connect the boundary advancement directions corresponding to the continuous migration boundaries, and perform segmented traction and spatial connection on the boundary change trend to obtain the progressive traction path; Map the underground anomaly accumulation direction corresponding to the progressive traction path, and perform intersection extraction and segmented connection on the underground anomaly accumulation direction to output the spatial pointing position corresponding to the underground water-rich area; By integrating the spatial orientation of the underground water-rich area with the underground anomaly expansion direction corresponding to the progressive traction path, and spatially connecting the locations that form a continuous expansion relationship within the underground water-rich area, the dynamic deduction of the spatial occurrence state of the karst aquifer is completed.

[0017] The technical effects and advantages provided by this application in the above technical solution are as follows: This application continuously tracks the propagation state of underground anomalies during continuous rainfall, dynamically correlates wave velocity shift amplitude, water pressure fluctuation state, and seepage transfer trajectory during the same time process, so that the migration direction change process of groundwater within the fracture area can be continuously expressed, thereby improving the ability to identify underground anomaly propagation paths, avoiding the problem of delayed judgment of water-rich areas caused by rapid deflection of groundwater migration direction, and improving the accuracy of dynamic prediction of the water inrush hazard range in complex karst areas.

[0018] This application performs progressive mapping and synchronous compression processing on the migration direction of underground anomalies to continuously converge the expression of the diffusion trajectory formed during the propagation of underground anomalies, and retains the corresponding migration direction information when the diffusion speed changes rapidly. This enables the hidden water-conducting migration state formed in the deep underground region to be identified in advance, thereby improving the dynamic adjustment capability of the underground water-rich boundary and reducing the impact of boundary distortion on the spatial endowment state inference results during the propagation of underground anomalies. Attached Figure Description

[0019] 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 application. For those skilled in the art, other drawings can be obtained based on these drawings.

[0020] Figure 1 This is the overall flowchart of this application. Detailed Implementation

[0021] 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.

[0022] This application provides, as follows: Figure 1 The intelligent simulation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation, as shown, includes the following steps: Step 1: Collect three-dimensional tomographic scan response information corresponding to different burial depths during continuous rainfall, record the wave velocity shift amplitude, water pressure fluctuation state and seepage transfer trajectory in the fracture area along the time progression direction, and generate a time-varying distribution map of the corresponding groundwater migration trend. The specific steps for generating a time-varying distribution map of the corresponding groundwater migration pattern are as follows: After continuous rainfall entered the underground karst area, three-dimensional tomographic scanning response information was continuously collected at different burial depths. When groundwater began to seep into the fracture area, the scanning area corresponding to each burial depth was cyclically probed in the order of time progression, and the spatial response state corresponding to the current time and location was recorded after each probe was completed.

[0023] As continuous rainfall progresses, the scanning areas corresponding to different burial depths are divided into multiple continuously distributed underground response intervals. The wave velocity change location, electrical response location, and water-bearing anomaly location within each underground response interval are recorded synchronously. At the same time, the spatial coordinates of each underground response interval at the current time location are uniformly calibrated, so that the three-dimensional tomographic scanning response information corresponding to different burial depth areas at the same time location can form a continuous correlation.

[0024] As continuous rainfall continues to enter the underground fissure area, the same area is scanned and collected repeatedly between adjacent time points. The response areas that change between the previous time point and the next time point are continuously preserved, and the underground response status corresponding to each time point is arranged in sequence according to the time progression direction, thereby forming a set of underground dynamic response information corresponding to the continuous rainfall process.

[0025] After completing the continuous recording of the underground dynamic response information set, the wave velocity shift amplitude corresponding to different time locations is extracted region by region. During the migration of groundwater along the fracture area to the depth, the locations where wave velocity changes occur between adjacent time locations are continuously tracked, and the movement direction of the wave velocity change location within the spatial range is recorded.

[0026] As continuous rainfall enters different stages, the areas of groundwater pressure change are tracked synchronously. The water pressure fluctuations at each time point are continuously arranged according to the direction of groundwater migration, and the propagation paths formed by the water pressure fluctuations at different burial depths are spatially connected. During the continuous flow of groundwater, the seepage areas that undergo positional changes within the fractured areas are continuously marked. The movement trajectories of the same seepage area at multiple time points are pieced together segment by segment, so that the flow direction of groundwater within the fractured areas can be continuously unfolded along the time progression direction.

[0027] Subsequently, the spatial locations corresponding to wave velocity shift amplitude, water pressure fluctuation state, and seepage transfer trajectory are uniformly correlated to form a synchronous association of the three types of underground response information at the same time location. The association results corresponding to different time locations are then continuously combined according to the time progression direction to obtain the dynamic response association information of the corresponding groundwater migration process.

[0028] It should be noted that: Water pressure fluctuation refers to the dynamic changes in groundwater pressure within a fractured region at different times and spatial locations during continuous rainfall replenishment and groundwater migration. This includes the rise, fall, fluctuation, local enhancement, continuous diffusion, and propagation and transfer of groundwater pressure. It is mainly used to characterize the pressure propagation state formed during the flow of groundwater within the fractured region, as well as the hydraulic disturbance characteristics corresponding to the expansion of groundwater into deeper areas along different migration directions. Thus, it reflects the dynamic changes in the flow activity, propagation direction, and range of underground anomalies at different time stages.

[0029] After forming dynamic response correlation information of the groundwater migration process, the changes in groundwater response at different time locations are continuously unfolded. The wave velocity shift area corresponding to the early time location is extended layer by layer to the subsequent time location, and the newly emerging wave velocity change area in the subsequent time location is spatially connected, so that the wave velocity changes generated during the groundwater migration process can form a continuous propagation zone between multiple time locations. As the groundwater flows deeper along the fracture area, the positional changes of water pressure fluctuations in the spatial range are continuously mapped. The extension paths formed by the same water pressure propagation direction at different burial depths are vertically connected, and the water pressure change areas corresponding to multiple time locations are horizontally connected, so that the pressure propagation process formed by the groundwater during continuous flow can be continuously presented in the spatial range.

[0030] As groundwater continues to migrate within the fractured area, the movement paths formed by the seepage transfer trajectory at different time locations are continuously superimposed. The seepage movement direction corresponding to the previous time location is connected segment by segment with the seepage expansion direction corresponding to the next time location, so that the migration trajectory formed by groundwater at different time stages can form a continuous expansion path in the spatial range.

[0031] Subsequently, the wave velocity shift region, water pressure propagation path, and seepage expansion path are uniformly mapped to the same underground space and unfolded sequentially according to the time progression direction, so that the dynamic changes formed during the groundwater migration process can form a continuous evolution state in both space and time.

[0032] After the underground dynamic changes form a continuous evolution state, the migration direction of groundwater at different time locations is unfolded as a whole, the migration path formed by groundwater at the early time location is continuously connected with the migration path formed at the later time location, and the locations where the migration direction of groundwater changes are recorded within the underground space.

[0033] During the continuous rainfall process, the changes in wave velocity shift, water pressure propagation, and seepage transfer of groundwater within the fractured area are synchronously correlated. The underground response areas corresponding to multiple time locations are continuously arranged according to the direction of groundwater migration, so that the spatial migration state of groundwater at different time stages can be continuously expressed within the same spatial range. As groundwater continues to expand into deeper areas, the migration paths of positional changes within the underground response areas are continuously extended, and the migration directions of groundwater at different burial depths are connected layer by layer, so that the spatial flow process of groundwater within the fractured area can continue to unfold along the time progression direction. Subsequently, the underground response change results corresponding to multiple time locations are unified and integrated to form a time-varying distribution map of the corresponding groundwater migration pattern within the underground space.

[0034] Step 2: Extract the abnormal drift regions corresponding to adjacent time points in the time-varying distribution map, track the wave velocity jump positions along the groundwater migration direction, and spatially arrange the continuously expanding offset regions according to the order of advancement to form a dynamic evolution region corresponding to the pulse seepage changes. The dynamic evolution region corresponding to the pulse seepage changes is formed through the following steps: After the time-varying distribution map of the corresponding groundwater migration pattern is formed, the underground response areas corresponding to adjacent moments in the time-varying distribution map are unfolded moment by moment. During the continuous rainfall entering the underground fissure area, the response changes formed at different time locations of the same spatial location are continuously compared, and the underground anomaly areas that have moved in location are extracted separately from the time-varying distribution map.

[0035] As groundwater flows along the fractured region into deeper areas, the location of the anomalous response at one moment is spatially correlated with the location of the anomalous response at the next moment. The area where displacement changes between the two time positions is recorded as the anomalous drift area, and the spatial movement direction corresponding to the anomalous drift area is recorded simultaneously. As continuous rainfall continues, underground anomalous areas that repeatedly shift spatially between multiple time positions are continuously preserved. The positional change process of the same anomalous area at different time positions is unfolded sequentially according to the time progression direction, so that the anomalous propagation state formed by groundwater within the fractured region can be continuously presented along the time progression direction.

[0036] During the continuous migration of groundwater to deeper areas, the abnormal drift areas formed at different burial depths are longitudinally connected, and the positional offset states formed in the shallow areas and the deep areas are spatially connected. This allows the abnormal propagation process of groundwater at different depths to extend continuously within the same underground space, thereby forming a set of abnormal drift areas corresponding to the groundwater migration process.

[0037] After forming a set of anomalous drift regions, the wave velocity response regions with positional changes within the anomalous drift regions are tracked layer by layer. During the continuous flow of groundwater, the locations where wave velocity abrupt changes occur between adjacent time positions are continuously extracted, and the locations where the wave velocity change amplitude exceeds the continuous change range of the surrounding area are recorded as wave velocity jump locations. During the migration of groundwater along the fracture region, the spatial movement direction corresponding to the wave velocity jump location is unfolded segment by segment, and the wave velocity jump region formed at the previous time position is spatially connected with the wave velocity jump region formed at the next time position, so that the wave velocity changes formed by groundwater at different time stages can form a continuous movement path within the underground space.

[0038] As groundwater continues to diffuse within the fractured region, wave velocity jump locations formed at multiple time points are continuously arranged according to the groundwater migration direction. Wave velocity jump areas formed at the forefront of the groundwater flow direction are arranged in the forward spatial position, and wave velocity jump areas formed subsequently are arranged in the subsequent extension direction. This allows the wave velocity propagation process formed by groundwater at different time stages to continuously advance along the underground space. During the migration of groundwater to different burial depth areas, the areas where wave velocity jump locations expand and change are recorded synchronously. The spatial expansion range formed by the same wave velocity jump area at multiple time points is continuously connected, allowing the wave velocity change trajectory formed by groundwater within the fractured region to form a continuous expansion state within the underground space.

[0039] After a continuous movement path is formed at the location of the wave velocity jump, the continuously expanding offset area within the abnormal drift area is unfolded time-by-time. As the groundwater continues to advance into deeper areas, the changes in the edge of the offset area formed at previous and subsequent time positions are continuously recorded, and the newly expanded locations within the offset area are synchronously preserved. As the direction of groundwater migration changes, the locations of rapid advancement within the offset area are continuously marked, and the rapid advancement area is spatially connected with the surrounding offset areas that are expanding and changing. This allows the spatial expansion process of groundwater within the fractured area to form a continuous advancement zone within the underground space.

[0040] As groundwater continues to penetrate deeper regions, the expansion direction of the same displacement area at multiple time points is unfolded segment by segment. The displacement areas formed at earlier time points are arranged at the forefront of the groundwater migration direction, and the displacement areas formed at later time points are arranged at the subsequent extension positions of the groundwater migration direction. This allows the spatial displacement process of groundwater at different time stages to continuously advance along the underground space. During the continuous migration of groundwater, the displacement areas corresponding to multiple time points are spatially connected, and the displacement expansion areas formed at adjacent time points are continuously linked. This allows the spatial expansion state of groundwater within the fracture area to form a continuous propagation process within the underground space.

[0041] After the continuously expanding offset region forms a continuous propagation process, the spatial offset states corresponding to multiple time positions are unfolded as a whole. As groundwater continues to flow along the fracture region, the direction of advancement of positional changes within the offset region is extended segment by segment, and the spatial offset paths formed at previous and subsequent time positions are continuously connected, so that the abnormal propagation state formed by groundwater at different time stages can form a continuous diffusion trajectory within the underground space. As groundwater continues to enter different burial depth areas, the movement direction formed by the wave velocity jump position, the expansion direction formed by the offset region, and the advancement direction formed by groundwater migration are synchronously corresponded. The spatial offset regions formed at multiple time positions are uniformly arranged according to the order of advancement, so that the abnormal drift state formed by groundwater during continuous rainfall can be continuously presented within the same underground space.

[0042] During the continuous diffusion of groundwater into deeper areas, the locations where advancement changes occur within the spatial range are continuously preserved, and the underground anomaly propagation paths formed between different burial depths are longitudinally connected, so that the pulse seepage changes formed by groundwater within the fracture area can continue to expand along the time progression direction; subsequently, the spatial offset paths, wave velocity jump paths, and offset expansion areas formed at multiple time locations are uniformly integrated to form a dynamic evolution area corresponding to the pulse seepage changes within the underground space.

[0043] Step 3: Perform progressive mapping on the hydraulic disturbance range corresponding to each moment in the dynamic evolution area, continuously project the deflection positions that repeatedly appear in a short period of time to the deep area, and record the diffusion amplitude of the abnormal drift direction at different time positions to generate an underground anomaly migration direction map. The specific steps for generating a subsurface anomaly migration map are as follows: After forming a dynamic evolution region corresponding to the pulse seepage changes, the range of underground anomalies at different time locations within the dynamic evolution region is continuously expanded. During the continuous migration of groundwater along the fracture region, the wave velocity change region, water pressure propagation region, and seepage expansion region formed at the same time location are simultaneously extracted, and the positional correspondence of multiple underground anomaly regions within the spatial range is uniformly recorded. As groundwater continues to advance into deeper regions, the boundaries of underground anomalies formed at different time locations are expanded layer by layer, continuously extending the boundaries of underground anomalies formed at the previous time location to the next time location, and synchronously connecting the newly formed underground anomaly expansion regions at the next time location, so that the anomaly propagation process formed by groundwater at different time stages can form a continuous advancing state within the underground space. As the direction of groundwater migration continuously changes between different burial depths, the range of underground anomaly expansion formed at different time locations is progressively mapped, continuously expanding the underground anomaly expansion state formed in the shallow region along the direction of groundwater migration to the deeper region, while spatially connecting the newly formed expansion regions in the deep region, so that the anomaly propagation process formed by groundwater within the fracture region can form a continuous expansion path within different depth ranges.

[0044] During the continuous flow of groundwater, the range of underground anomalies formed at each time point is arranged continuously in chronological order. The underground anomaly regions formed at earlier time points are arranged at the front of the groundwater migration direction, and the underground anomaly regions formed at later time points are arranged at the extension of the groundwater migration direction. This allows the range of hydraulic disturbances formed by groundwater during continuous rainfall to form a continuous progressive state along the time progression direction, thereby forming a progressive mapping region corresponding to the underground anomaly propagation process.

[0045] After the progressive mapping region of the underground anomaly propagation process is formed, the deflection positions that repeatedly appear in a short period of time are continuously extracted. As the groundwater continues to advance along the fracture region, the underground anomaly regions with positional changes formed at multiple adjacent time positions are compared segment by segment, and the positions with repeated directional changes within the same spatial range are recorded as repeated deflection positions. As the groundwater migration direction changes, the underground anomaly propagation direction corresponding to the repeated deflection position is continuously unfolded, and the deflection direction formed at the previous time position is spatially connected with the deflection direction formed at the next time position, so that the changes in the propagation direction formed by the groundwater at different time stages can form a continuous deflection trajectory within the underground space.

[0046] As groundwater continues to expand into deeper regions, the spatial movement direction formed by repeated deflection points is extended layer by layer. The deflection paths formed in the shallow region are continuously projected into the deep region, and the propagation directions of newly formed underground anomalies in the deep region are synchronously continued. This allows the propagation changes of groundwater within the fractured area to continue to advance along the underground space. During the groundwater migration process, the deflection paths formed at multiple time points are continuously arranged according to the groundwater migration direction. The deflection paths formed at the forefront of the groundwater flow are arranged at the front position of the underground space, and the deflection paths formed subsequently are arranged in the subsequent extension directions. This allows the deflection changes of groundwater at different time stages to form a continuous projection state within the underground space. Subsequently, the deflection paths formed at multiple time points are vertically connected, allowing the propagation directions of underground anomalies formed by groundwater at different burial depths to continue to expand along the deep region, thereby forming a deep continuous projection path corresponding to the underground anomaly propagation process.

[0047] After the formation of a continuous projection path in the deep subsurface, the diffusion amplitude of the anomaly drift direction at different time locations is recorded hour by hour. During the continuous migration of groundwater along the fracture area, the underground anomaly expansion boundary formed at previous and subsequent time locations is continuously extracted, and the positional change distance of the underground anomaly expansion boundary within the spatial range is recorded. During the continuous change of the groundwater propagation direction, the anomaly drift direction formed at different time locations is simultaneously unfolded, and the expansion state formed at multiple time locations by the same anomaly drift direction is connected segment by segment, so that the underground anomaly propagation process formed by groundwater at different time stages can form a continuous diffusion trajectory within the underground space.

[0048] As groundwater continues to penetrate deeper regions, the diffusion range corresponding to the anomaly drift direction is continuously extended, spatially connecting the diffusion areas formed at earlier and later time points. This allows the underground anomaly propagation state formed within the fractured area to continuously expand along the time progression. During the continuous flow of groundwater, the diffusion amplitudes formed at different time points are uniformly arranged according to the groundwater migration direction. Locations with faster underground anomaly expansion rates are synchronously correlated with locations where the underground anomaly expansion rates continuously change, enabling the underground anomaly propagation direction formed at different time stages to form a continuously changing state within the underground space. Subsequently, the anomaly drift directions, deep continuous projection paths, and diffusion amplitude change areas formed at multiple time points are uniformly integrated to generate an underground anomaly migration direction map within the underground space.

[0049] Step 4: Read the deflection areas that are continuously expanding in the underground anomaly migration direction map, perform synchronous compression on the diffusion trajectories corresponding to multiple time locations, and retain the migration direction information corresponding to the areas with a sudden increase in diffusion speed during the compression process to obtain the clustering areas of the corresponding hidden water-rich migration state. The specific steps to obtain the clustering areas corresponding to the hidden water-rich migration state are as follows: After the formation of the underground anomaly migration pattern, the deflection areas that continuously expand along the direction of groundwater migration within the pattern are continuously read. During the continuous rainfall entering the fracture area, the underground anomaly expansion boundaries formed at different time points are unfolded moment by moment, and the underground anomaly expansion range formed at each time point is projected onto the same underground space area. As groundwater flows along the fracture area to deeper areas, the spatial positions of the deflection areas formed at the previous time point and the deflection areas formed at the next time point are spatially correlated. The underground anomaly areas that form a continuous expansion relationship between the two time points are connected segment by segment, and the advancing direction and distance of the underground anomaly expansion boundary in the underground space are recorded. As the groundwater continues to expand into deeper areas, the deflection areas formed at multiple time points are extended longitudinally, and the underground anomaly expansion path formed in the shallow area is continuously advanced into the deep area. The propagation direction of newly formed underground anomalies in the deep area is synchronously continued, so that the underground anomaly propagation state formed by groundwater within the fracture area can form a continuous expansion path along the underground space.

[0050] During the continuous migration of groundwater, the spatial boundaries of underground anomalies formed at different time locations are connected. The underground anomaly expansion areas formed at adjacent time locations are arranged in chronological order. The deflection areas formed at the forefront of groundwater migration are arranged at the front of the underground space, and the subsequently formed deflection areas are arranged sequentially at the extension positions of the groundwater migration direction. This allows the underground anomaly propagation direction to form a continuous expansion state within the underground space during continuous rainfall. Subsequently, the underground anomaly expansion directions formed at multiple time locations are uniformly unfolded, and the underground anomaly propagation paths formed at different burial depths are continuously connected, thereby forming the continuously expanding deflection areas corresponding to the underground anomaly propagation process.

[0051] After the formation of the continuously expanding deflection zone, the underground anomalous diffusion trajectories corresponding to multiple time locations are continuously extracted. During the continuous migration of groundwater along the fracture area, the underground anomalous diffusion path formed at each time location is unfolded segment by segment, and the positional offset direction of the underground anomalous diffusion path within the underground space is recorded synchronously. During the continuous change of the groundwater propagation direction, the underground anomalous diffusion trajectory formed at the previous time location is spatially correlated with the underground anomalous diffusion trajectory formed at the next time location. The underground anomalous diffusion trajectories with similar propagation directions are gradually converged into the same underground space, while the underground anomalous diffusion trajectories with continuously increasing spatial deviation directions are edge-contracted, so that the underground anomalous propagation process formed by groundwater at different time stages can form a continuous compression path along the underground space.

[0052] As groundwater continues to advance into deeper areas, it synchronously compresses the underground anomaly diffusion trajectories formed at multiple time points, gradually bringing the underground anomaly diffusion paths formed at earlier time points closer to those formed at later time points, and continuously pulling the new diffusion directions formed at later time points into the compressed area, so that the underground anomaly propagation state formed by groundwater in the fracture area can form a continuous convergence process along the underground space.

[0053] During the continuous flow of groundwater, the underground anomaly diffusion trajectories after synchronous compression are longitudinally connected, spatially connecting the compression trajectories formed in shallow and deep regions. At the same time, the changes in the direction of underground anomaly diffusion formed at different time locations during the compression process are recorded, so that the direction of underground anomaly propagation formed by groundwater at different burial depths can form a continuous convergence state along the underground space. Subsequently, the synchronous compression trajectories formed at multiple time locations are uniformly arranged, so that the underground anomaly diffusion path forms a continuous contraction area within the underground space, thereby forming the synchronous compression trajectory area corresponding to the underground anomaly propagation process.

[0054] After the synchronous compression trajectory area is formed, the underground anomaly expansion boundary that changes position during the compression process is continuously unfolded. As the groundwater continues to advance into the deeper area, the underground anomaly expansion distance formed at different time positions is recorded moment by moment, and the underground anomaly expansion distance corresponding to each time position is spatially correlated with the underground anomaly expansion distance formed at the previous time position.

[0055] During rapid changes in the direction of groundwater migration, underground anomaly expansion areas that rapidly advance within a short period are continuously extracted. Locations where the advance distance suddenly increases between adjacent time points are recorded as areas of rapid increase in diffusion velocity. As groundwater continues to penetrate deeper areas, the underground anomaly propagation direction corresponding to the areas of rapid increase in diffusion velocity is unfolded layer by layer. The rapid advance direction formed at the previous time point is spatially connected with the rapid advance direction formed at the next time point, enabling the rapid propagation state of underground anomalies formed by groundwater at different time stages to form a continuous advancement path along the underground space. During continuous groundwater flow, the positional change direction of the areas of rapid increase in diffusion velocity is synchronously preserved. The rapid advance directions formed at multiple time points are continuously arranged according to the groundwater migration direction. Simultaneously, the positional change range and extension length of the rapid advance direction within the underground space are recorded, enabling the rapid propagation state of underground anomalies formed by groundwater within the fractured area to form a continuous extension along the underground space.

[0056] Subsequently, the regions with a sudden increase in diffusion rate formed at multiple time locations were vertically connected, and the rapid propagation direction formed in the shallow region was spatially connected with the rapid propagation direction formed in the deep region, thereby forming the migration direction region corresponding to the underground abnormal rapid propagation process.

[0057] After the migration-direction area is formed, the synchronous compression trajectory area and the migration-direction area are continuously merged. During the continuous flow of groundwater along the fracture area, the diffusion directions of underground anomalies formed at multiple time points are uniformly unfolded, and the locations where the diffusion directions of underground anomalies intersect within the underground space are synchronously extracted. During the continuous expansion of groundwater into deeper areas, the locations where the underground anomaly propagation directions continuously converge are connected segment by segment. The underground anomaly aggregation directions formed at the previous time point are spatially connected with the underground anomaly aggregation directions formed at the next time point, so that the underground anomaly propagation states formed by groundwater at different time stages can form a continuous aggregation path along the underground space.

[0058] As groundwater continues to migrate, the underground anomaly propagation areas that repeatedly intersect at multiple time locations are continuously preserved. The locations where the underground anomaly propagation directions repeatedly overlap are recorded as underground anomaly accumulation locations. The expansion areas formed by the underground anomaly accumulation locations within the underground space are simultaneously expanded, so that the underground anomaly propagation state formed by groundwater within the fracture area can form a continuous convergence process along the underground space.

[0059] During the continuous flow of groundwater, the migration direction information corresponding to the area with a sudden increase in diffusion rate is spatially fused with the synchronously compressed underground anomaly diffusion trajectory. At the same time, the location change direction of the underground anomaly accumulation location is continuously corresponded, so that the rapid propagation direction of underground anomalies formed by groundwater at different time stages and the convergence direction of underground anomalies can form a continuous accumulation state within the same underground space.

[0060] Subsequently, the underground anomalous clusters formed at multiple time points were uniformly arranged, and the locations within the underground anomalous clusters that formed continuous expansion relationships were spatially connected to obtain the clusters corresponding to the hidden water-rich migration states.

[0061] Step 5: Extract the continuously expanding migration boundary within the aggregation area, perform gradual traction on the boundary change trend along the time progression, and adjust the spatial orientation position corresponding to the underground water-rich area during the traction process to complete the dynamic deduction of the spatial occurrence state of the karst aquifer. During the traction and propulsion process, the spatial orientation of the underground water-rich area is adjusted to complete the dynamic simulation of the spatial occurrence state of the karst aquifer. The specific steps are as follows: After the formation of the accumulation area corresponding to the hidden water-rich migration state, the migration boundary of the accumulation area, which continuously expands along the direction of groundwater migration, is unfolded moment by moment. During the continuous rainfall entering the fissure area, the range of underground anomaly accumulation formed at different time positions is continuously read, and the edge of underground anomaly accumulation formed at each time position is projected into a unified underground space. During the flow of groundwater along the fissure area to the deeper area, the spatial position of the underground anomaly accumulation edge formed at the previous time position is matched with the edge of underground anomaly accumulation formed at the next time position. The underground anomaly accumulation areas that form a continuous expansion relationship between the two time positions are connected segment by segment. At the same time, the advancing direction, advancing distance, and edge expansion range of the underground anomaly accumulation edge in the underground space are recorded. As the groundwater continues to spread to the deeper area, the underground anomaly accumulation edges formed at multiple time positions are extended longitudinally. The underground anomaly accumulation expansion path formed in the shallow area is continuously advanced to the deeper area, and the new accumulation direction formed in the deeper area is synchronously continued, so that the underground anomaly convergence state formed by groundwater in the fissure area can form a continuous expansion path along the underground space.

[0062] During the continuous migration of groundwater, the edges of underground anomalies formed at different time locations are uniformly arranged in chronological order. The edges formed at the forefront of groundwater migration are arranged at the front of the underground space, and the edges formed subsequently are arranged sequentially at the extension positions of the groundwater migration direction. This ensures that the underground anomaly convergence direction continuously expands within the underground space during continuous rainfall. Subsequently, the underground anomaly convergence edges formed at multiple time locations are uniformly unfolded, and the underground anomaly convergence paths formed at different burial depths are continuously connected. At the same time, the spatial offset direction and edge extension length of the underground anomaly convergence edges between different time locations are recorded, thereby forming the continuous migration boundary of the corresponding underground anomaly convergence process.

[0063] After the continuous migration boundary is formed, the direction of boundary change at different time locations is unfolded segment by segment. As groundwater continues to flow along the fracture area, the direction of boundary advancement formed at the previous time location is continuously correlated with the direction of boundary advancement formed at the next time location, and the direction of position change of the underground anomaly accumulation edge within the underground space is recorded synchronously. As the direction of groundwater migration continues to change, the trend of boundary change formed at multiple time locations is spatially arranged according to the time progression. The direction of boundary advancement formed at the earlier time location is gradually pulled towards the direction of boundary advancement formed at the later time location. At the same time, the direction of expansion of the newly formed boundary is continuously converged into the interior of the existing advancement path, so that the underground anomaly accumulation state formed by groundwater at different time stages can form a continuous and gradual change path along the underground space.

[0064] As groundwater continues to advance into deeper areas, a gradual traction process is implemented on the boundary change areas formed at multiple time points. The location change paths formed by the underground anomaly accumulation edge at different time stages are extended layer by layer, and the locations where the direction deflection occurs during the boundary advancement are synchronously connected. This allows the underground anomaly expansion state formed by groundwater within the fracture area to form a continuous traction process along the underground space.

[0065] During the continuous flow of groundwater, the boundary advancement paths after gradual traction are longitudinally connected, spatially connecting the boundary traction directions formed in the shallow and deep regions. Simultaneously, the expansion location, advancement length, and directional change range of the boundary change area within the underground space are recorded, enabling the underground anomaly accumulation state formed by groundwater at different burial depths to form a continuous extension along the underground space. Subsequently, the boundary traction directions formed at multiple time points are uniformly arranged, synchronously preserving the locations where repeated expansion relationships occur during the boundary advancement process. The boundary change paths that form continuous advancement relationships between multiple time points are spatially connected, enabling the underground anomaly accumulation boundary to form a continuous advancement area within the underground space, thus forming a gradual traction path corresponding to the underground anomaly expansion process.

[0066] After the gradual traction path is formed, the spatial orientation of the underground water-rich area is continuously adjusted. As the groundwater continues to migrate along the fracture area, the underground anomaly accumulation direction formed at different time positions is uniformly unfolded, and the intersection of the underground anomaly accumulation direction within the underground space is extracted synchronously. As the groundwater continues to expand into deeper areas, the positions where the underground anomaly propagation direction continues to converge are connected segment by segment. The underground anomaly accumulation direction formed at the previous time position is spatially connected with the underground anomaly accumulation direction formed at the next time position, so that the underground anomaly convergence state formed by the groundwater at different time stages can form a continuous accumulation path along the underground space.

[0067] As groundwater continues to migrate, underground anomaly accumulation areas that repeatedly intersect at multiple time locations are continuously preserved. The locations where the directions of underground anomaly accumulation repeatedly overlap are recorded as the spatial pointing locations corresponding to underground water-rich areas. The spatial expansion range corresponding to underground water-rich areas is simultaneously expanded, so that the underground anomaly propagation state formed by groundwater in the fracture area can form a continuous convergence process along the underground space.

[0068] During the continuous flow of groundwater, the progressive traction path is spatially integrated with the spatial orientation of the groundwater-rich area. At the same time, the direction of positional change of the groundwater-rich area is continuously correlated, so that the direction of underground anomaly expansion and convergence formed by groundwater at different time stages can form a continuous migration state within the same underground space. Subsequently, the spatial orientation of the groundwater-rich area formed at multiple time points is uniformly arranged, and the positions of the continuous expansion relationship within the groundwater-rich area are spatially connected. At the same time, the expansion boundary, migration direction and spatial extension range formed by the groundwater-rich area at different time points are recorded, thereby completing the dynamic inference of the spatial occurrence state of the karst aquifer.

[0069] The foregoing has only described certain exemplary embodiments of this application 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 this application. 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 this application.

Claims

1. An intelligent method for extrapolating the spatial occurrence state of karst aquifers in the face of uncertainty propagation, characterized in that, The steps include the following: Three-dimensional tomographic scan response information at different burial depths during continuous rainfall was collected. Wave velocity shift amplitude, water pressure fluctuation state and seepage transfer trajectory were recorded along the time progression direction within the fracture area, and a time-varying distribution map of the corresponding groundwater migration pattern was generated. Extract the abnormal drift regions corresponding to adjacent moments in the time-varying distribution map, track the wave velocity jump positions along the groundwater migration direction, and spatially arrange the continuously expanding offset regions according to the order of advancement to form a dynamic evolution region corresponding to the pulse seepage changes; A progressive mapping is performed on the hydraulic disturbance range corresponding to each moment in the dynamic evolution area. The deflection positions that repeatedly appear in a short period of time are continuously projected to the deep region, and the diffusion amplitude of the abnormal drift direction at different time positions is recorded to generate an underground anomaly migration direction map. Read the deflection areas that are continuously expanding within the underground anomaly migration direction map, perform synchronous compression on the diffusion trajectories corresponding to multiple time locations, and retain the migration direction information corresponding to the areas with a sudden increase in diffusion speed during the compression process to obtain the cluster areas corresponding to the hidden water-rich migration state. Extract the continuously expanding migration boundary within the aggregation area, perform gradual traction on the boundary change trend along the time progression, and adjust the spatial orientation of the underground water-rich area during the traction process to complete the dynamic deduction of the spatial occurrence state of the karst aquifer. An abnormal drift area is an underground anomaly area that undergoes displacement between two time locations after spatially corresponding the location where an abnormal response occurs at one time and the location where an abnormal response occurs at the next time. The location of a sudden jump in wave velocity is: a position where the amplitude of the wave velocity change exceeds the range of continuous change in the surrounding area; The specific process of progressive mapping is as follows: the underground anomaly expansion state formed in the shallow region is continuously unfolded towards the deep region along the direction of groundwater migration, while the newly expanded region formed in the deep region is spatially connected, and the range of underground anomaly changes formed at each time position is continuously arranged in the order of time progression. The specific process of continuous projection is as follows: the spatial movement direction formed by repeated deflection positions is extended layer by layer, the deflection path formed in the shallow area is continuously projected to the deep area, and the propagation direction of newly formed underground anomalies in the deep area is synchronously continued. During the groundwater migration process, the deflection paths formed at multiple time points are arranged continuously according to the groundwater migration direction. The deflection paths formed at the front of the groundwater flow are arranged at the front of the underground space, and the deflection paths formed later are arranged in the subsequent extension direction. Then, the deflection paths formed at multiple time points are connected longitudinally to form a deep continuous projection path corresponding to the underground anomaly propagation process. The specific process of synchronous compression is as follows: spatially correspond the underground anomaly diffusion trajectory formed at the previous time position with the underground anomaly diffusion trajectory formed at the next time position, gradually converge the underground anomaly diffusion trajectories with similar propagation directions to the same underground space range, and at the same time, perform edge contraction on the underground anomaly diffusion trajectory with continuously increasing spatial deviation direction. The specific process of gradual traction is as follows: the boundary change trends formed at multiple time positions are spatially arranged according to the time progression sequence, the boundary advancement direction formed at the early time positions is gradually pulled towards the boundary advancement direction formed at subsequent time positions, and at the same time, the expansion direction of the newly formed boundary is continuously converged into the interior of the existing advancement path.

2. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 1, characterized in that, The specific steps for generating a time-varying distribution map of the corresponding groundwater migration pattern are as follows: Cyclic detection was performed on the scanning areas corresponding to different burial depths, and the spatial response status at each time location was recorded to construct a set of underground dynamic response information. Extract the wave velocity offset amplitude, water pressure fluctuation state and seepage transfer trajectory corresponding to the underground dynamic response information set, and spatially connect the wave velocity change location, water pressure propagation path and seepage movement trajectory to obtain dynamic response correlation information; The wave velocity offset region, water pressure propagation path and seepage expansion path corresponding to the dynamic response association information are continuously expanded to obtain the continuous evolution state of underground dynamic changes. The results of underground response changes corresponding to the continuous evolution of underground dynamic changes are uniformly integrated to form a time-varying distribution map of the corresponding groundwater migration pattern.

3. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 2, characterized in that, Spatial connection of wave velocity change location, water pressure propagation path and seepage movement trajectory, specifically including: longitudinal connection of extension paths formed by different burial depths in the same water pressure propagation direction, and lateral connection of water pressure change areas corresponding to multiple time locations.

4. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 2, characterized in that, The dynamic evolution region corresponding to the pulse seepage changes is formed through the following steps: Extract underground anomaly regions that shift in position at adjacent times within the time-varying distribution map, and synchronously record the spatial shift direction corresponding to the anomaly drift regions to obtain a set of anomaly drift regions; Based on the set of abnormal drift regions, the wave speed jump positions formed by adjacent time positions are continuously extracted, and the spatial movement direction corresponding to the wave speed jump position is expanded segment by segment to form a continuous movement path. Based on the continuous movement path, the continuously expanding offset region inside the abnormal drift region is unfolded time by time, and the positions where the offset region expands and changes are spatially connected to construct a continuous advancement zone. By unifying and integrating the spatial migration paths, wave velocity jump paths, and migration extension regions formed at multiple time locations around the continuous advancement zone, a dynamic evolution region corresponding to the pulse seepage changes is formed.

5. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 4, characterized in that, The wave velocity jump positions formed by adjacent time positions are continuously extracted, including: recording the position where the wave velocity change amplitude exceeds the continuous change range of the surrounding area as the wave velocity jump position, and continuously connecting the spatial expansion range formed by the same wave velocity jump area at multiple time positions.

6. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 4, characterized in that, The specific steps for generating a subsurface anomaly migration map are as follows: The wave velocity change region, water pressure propagation region, and seepage expansion region corresponding to different time locations within the dynamic evolution region are obtained, and the underground anomaly change boundary is continuously extended and spatially connected to obtain the progressive mapping region. Under the progressive mapping region, the positions where the direction changes repeatedly occur within the same spatial range are recorded, and the deflection path is extended and connected vertically layer by layer to construct a deep continuous projection path; Based on the deep continuous projection path, the diffusion range corresponding to the anomaly drift direction is continuously extended, and the diffusion amplitude formed at multiple time locations is uniformly arranged to generate an underground anomaly migration direction map.

7. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 6, characterized in that, Record the locations where directional changes occur repeatedly within the same spatial range, including: comparing underground anomaly areas corresponding to multiple adjacent time locations segment by segment, recording the locations where directional changes occur repeatedly within the same spatial range as repeated deflection locations, and spatially connecting the deflection direction corresponding to the previous time location with the deflection direction corresponding to the next time location.

8. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 6, characterized in that, The specific steps to obtain the clustering area corresponding to the hidden water-rich migration state are as follows: Read the deflection area that continuously expands along the direction of groundwater migration in the map, and expand and spatially connect the underground anomaly expansion boundary corresponding to different time locations to output the continuously expanding deflection area. For the continuously expanding deflection area, extract the underground anomaly diffusion trajectory corresponding to multiple time locations, and perform gradual convergence and edge contraction on the underground anomaly diffusion trajectories with similar propagation directions to obtain the synchronously compressed trajectory area. The system records the advance distance corresponding to the underground anomaly expansion boundary in the synchronous compression trajectory region, extracts the locations where the advance distance suddenly increases, and spatially connects the underground anomaly propagation directions corresponding to the regions where the diffusion velocity suddenly increases to construct the migration pointing region. The underground anomaly propagation directions corresponding to the migration direction area and the synchronous compression trajectory area are collected, and the underground anomaly accumulation locations are spatially connected and uniformly arranged to obtain the accumulation area corresponding to the hidden water-rich migration state.

9. The intelligent extrapolation method for the spatial occurrence state of karst aquifers oriented towards uncertainty propagation as described in claim 8, characterized in that, The specific steps for dynamically extrapolating the spatial occurrence state of karst aquifers are as follows: Expand the underground anomaly cluster edges corresponding to different time locations within the cluster area, and perform spatial location correspondence and segment-by-segment connection on the underground anomaly cluster edges to construct a continuous migration boundary; Connect the boundary advancement directions corresponding to the continuous migration boundaries, and perform segmented traction and spatial connection on the boundary change trend to obtain the progressive traction path; Map the underground anomaly accumulation direction corresponding to the progressive traction path, and perform intersection extraction and segmented connection on the underground anomaly accumulation direction to output the spatial pointing position corresponding to the underground water-rich area; By integrating the spatial orientation of the underground water-rich area with the underground anomaly expansion direction corresponding to the progressive traction path, and spatially connecting the locations that form a continuous expansion relationship within the underground water-rich area, the dynamic deduction of the spatial occurrence state of the karst aquifer is completed.

Images