Underground water stability judgment method and system

Abstract

The invention provides a groundwater stability judgment method and system. The method comprises the following steps: acquiring ground temperature data, seismic data and geological survey data of a geothermal area; judging whether the ground temperature data meets a preset ground temperature data condition or not, and if yes, judging that underground water in the geothermal area is unstable; if not, whether an earthquake exists in the geothermal area is judged, if the earthquake exists, whether the earthquake data meet preset earthquake data conditions is judged, and if yes, it is judged that underground water in the geothermal area is unstable; if not, whether the geological survey data meet the preset geological survey data condition or not is judged, and if yes, it is judged that the underground water in the geothermal area is unstable. According to the technical scheme provided by the invention, accurate identification of underground water stability risks is realized through triple layered judgment logic of ground temperature anomaly features, earthquake activity modes and geologic structure verification, the misjudgment rate is reduced, the judgment efficiency is improved, and the risk early warning capability is enhanced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of geological engineering, and particularly relates to a groundwater stability judgment method and system. BACKGROUND

[0002] Groundwater stability refers to the ability of a groundwater system to maintain a relatively balanced or predictable change state in terms of water quantity, water quality and dynamic characteristics (such as water level, flow rate, etc.) within a specific spatiotemporal range. Abnormal changes in groundwater level (such as rapid rise or fall) can significantly affect the stability of soil structures such as slopes and foundation pits, causing landslides, foundation pit failures, and other geological disasters. By monitoring groundwater stability, early warning of landslides, ground subsidence and other geological disasters can be provided to provide scientific basis for emergency response and engineering protection. In related technologies, only a single parameter such as ground temperature or seismic activity is used as the basis for judgment in judging groundwater stability, without considering the synergistic effect of ground temperature anomaly characteristics and seismic activity patterns, which can easily lead to misjudgment, and the accuracy of identifying and the efficiency of judging groundwater stability risk are low. Therefore, there is an urgent need to propose a scheme that can improve the accuracy of identifying and the efficiency of judging groundwater stability risk. SUMMARY

[0003] The present application provides a groundwater stability judgment method and system to at least solve the technical problem of low accuracy and efficiency in judging groundwater stability.

[0004] The first aspect of the present application provides a groundwater stability judgment method, which comprises:

[0005] Obtaining ground temperature data, seismic data and geological survey data of a geothermal area;

[0006] Judging whether the ground temperature data meets a preset ground temperature data condition, and if the ground temperature data meets the preset ground temperature data condition, determining that the groundwater in the geothermal area is unstable;

[0007] If the ground temperature data does not meet the preset ground temperature data condition, judging whether there is an earthquake in the geothermal area, and if there is an earthquake, judging whether the seismic data meets a preset seismic data condition, and if the seismic data meets the preset seismic data condition, determining that the groundwater in the geothermal area is unstable;

[0008] If the seismic data does not meet the preset seismic data condition, judging whether the geological survey data meets a preset geological survey data condition, and if the geological survey data meets the preset geological survey data condition, determining that the groundwater in the geothermal area is unstable.

[0009] Preferably, the preset geothermal data condition comprises: the geothermal temperature being higher than the average geothermal temperature of the background area of the geological survey by 60℃ or more, or the geothermal temperature being greater than the maximum value of the average geothermal temperature of all stations in the background area of the geological survey, and the number of stations with geothermal temperature greater than the average geothermal temperature being less than half of the total number of stations in the local area.

[0010] Preferably, the preset seismic data condition comprises: the number of earthquakes in the region within one month being more than 20, and the average magnitude of the earthquakes being less than or equal to 2.5, or the number of earthquakes in the region within one year being greater than or equal to 10, and the average magnitude of the earthquakes being less than or equal to 3.5.

[0011] Preferably, the preset geological survey data condition comprises: the existence of unstable strata in the stratum structure in the geological survey data.

[0012] Preferably, the geothermal region refers to a region within a range of 20℃ higher than the geothermal background value of the background area of the geological survey.

[0013] Preferably, the seismic data is seismic data obtained according to the number of earthquakes, the magnitude of the earthquakes, and the location of the earthquakes during a high-incidence period of earthquakes.

[0014] Preferably, the geological survey data refers to stratum structure data obtained during geological survey, and the geological survey comprises: taking a seismic peak line as an axis, taking a drill hole profile of a certain length on both sides of the seismic strike, and obtaining geological data at the profile.

[0015] The second aspect embodiment of the present application provides a groundwater stability judgment system, comprising:

[0016] A data acquisition module is configured to acquire geothermal data, seismic data, and geological survey data of a geothermal region.

[0017] A geothermal condition judgment module is configured to judge whether the geothermal data meets a preset geothermal data condition, and if the geothermal data meets the preset geothermal data condition, to determine that the groundwater in the geothermal region is unstable.

[0018] A seismic condition judgment module is configured to, if the geothermal data does not meet the preset geothermal data condition, judge whether there is an earthquake in the geothermal region, and if there is an earthquake, judge whether the seismic data meets a preset seismic data condition, and if the seismic data meets the preset seismic data condition, determine that the groundwater in the geothermal region is unstable.

[0019] A geological condition judgment module is configured to, if the seismic data does not meet the preset seismic data condition, judge whether the geological survey data meets a preset geological survey data condition, and if the geological survey data meets the preset geological survey data condition, determine that the groundwater in the geothermal region is unstable.

[0020] The third aspect of the present application provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method according to the first aspect.

[0021] The fourth aspect of the present application provides a computer readable storage medium, which stores a computer program, wherein the program is executed by a processor to implement the method according to the first aspect.

[0022] The embodiments of the present application provide at least the following beneficial effects:

[0023] The present application provides a groundwater stability judgment method and system, the method comprising: obtaining geothermal region ground temperature data, seismic data and geological survey data; judging whether the ground temperature data meets the preset ground temperature data condition, if the ground temperature data meets the preset ground temperature data condition, determining that the groundwater in the geothermal region is unstable; if the ground temperature data does not meet the preset ground temperature data condition, judging whether there is an earthquake in the geothermal region, if there is an earthquake, judging whether the seismic data meets the preset seismic data condition, if the seismic data meets the preset seismic data condition, determining that the groundwater in the geothermal region is unstable; if the seismic data does not meet the preset seismic data condition, judging whether the geological survey data meets the preset geological survey data condition, if the geological survey data meets the preset geological survey data condition, determining that the groundwater in the geothermal region is unstable. The technical scheme provided by the present application realizes accurate identification of groundwater stability risk by using three-layer judgment logic of ground temperature anomaly characteristics, seismic activity mode and geological structure verification, reduces the misjudgment rate, improves the judgment efficiency, and enhances the risk early warning capability.

[0024] Additional aspects and advantages of the present application will be made apparent by the following description and the accompanying drawings. BRIEF DESCRIPTION OF DRAWINGS

[0025] The above and / or additional aspects and advantages of the present application will become apparent and be readily understood from the following description, taken in conjunction with the accompanying drawings, in which:

[0026] Figure 1 A flowchart of a groundwater stability judgment method according to an embodiment of the present application is provided.

[0027] Figure 2 A structural block diagram of a groundwater stability judgment system according to an embodiment of the present application is provided. DETAILED DESCRIPTION

[0028] Embodiments of the present application are described below in detail with reference to the accompanying drawings, wherein the same or similar components or elements are denoted throughout the drawings by the same or similar reference numerals. The embodiments described below are exemplary and are intended to explain the present application, and should not be understood as limiting the present application.

[0029] A flowchart of a groundwater stability judgment method according to an embodiment of the present application is shown in FIG. 1, which includes the following steps. Figure 1 The method includes obtaining geothermal region geothermal data, seismic data and geological survey data; determining whether the geothermal data meets a preset geothermal data condition, and if the geothermal data meets the preset geothermal data condition, determining that the groundwater in the geothermal region is unstable; if the geothermal data does not meet the preset geothermal data condition, determining whether the geothermal region has an earthquake, and if there is an earthquake, determining whether the seismic data meets a preset seismic data condition, and if the seismic data meets the preset seismic data condition, determining that the groundwater in the geothermal region is unstable; if the seismic data does not meet the preset seismic data condition, determining whether the geological survey data meets a preset geological survey data condition, and if the geological survey data meets the preset geological survey data condition, determining that the groundwater in the geothermal region is unstable. The technical solution of the present application realizes accurate identification of groundwater stability risk by using three-layer judgment logic of geothermal anomaly characteristics, seismic activity mode and geological structure verification, reduces the misjudgment rate, improves the judgment efficiency and enhances the risk early warning capability.

[0030] A groundwater stability judgment method and system according to an embodiment of the present application are described below with reference to the accompanying drawings.

[0031] Embodiment One

[0032] Figure 1 A flowchart of a groundwater stability judgment method according to an embodiment of the present application is shown in FIG. 1, which includes the following steps. Figure 1 The method includes obtaining geothermal region geothermal data, seismic data and geological survey data; determining whether the geothermal data meets a preset geothermal data condition, and if the geothermal data meets the preset geothermal data condition, determining that the groundwater in the geothermal region is unstable; if the geothermal data does not meet the preset geothermal data condition, determining whether the geothermal region has an earthquake, and if there is an earthquake, determining whether the seismic data meets a preset seismic data condition, and if the seismic data meets the preset seismic data condition, determining that the groundwater in the geothermal region is unstable; if the seismic data does not meet the preset seismic data condition, determining whether the geological survey data meets a preset geological survey data condition, and if the geological survey data meets the preset geological survey data condition, determining that the groundwater in the geothermal region is unstable. The technical solution of the present application realizes accurate identification of groundwater stability risk by using three-layer judgment logic of geothermal anomaly characteristics, seismic activity mode and geological structure verification, reduces the misjudgment rate, improves the judgment efficiency and enhances the risk early warning capability.

[0033] Step 1: Obtain geothermal region geothermal data, seismic data and geological survey data.

[0034] Geothermal data is the core information for evaluating the thermal field distribution of geothermal areas, which can reveal the mechanism of groundwater flow and heat exchange. It can be divided into three categories: surface temperature, shallow geothermal temperature and deep well geothermal temperature. The specific collection methods are as follows: 1. Surface temperature data reflects the characteristics of surface thermal anomalies (such as temperature distribution in hot spring and geothermal display area). Collection method: select representative measuring points (cover different landforms and tectonic units) to ensure uniform distribution; use standardized thermometer or thermal infrared remote sensing equipment (such as thermal imager) to record surface temperature values; record measurement time, weather conditions (such as solar radiation, wind speed) and other auxiliary information simultaneously for subsequent correction; construct surface temperature field distribution map by interpolation method (such as Kriging interpolation) to identify thermal anomaly areas. 2. Shallow geothermal temperature data (<100 meters) represents the characteristics of near-surface heat conduction (such as soil heat flow and shallow thermal reservoir temperature). Collection method: place geothermometer or heat flow meter (such as thermistor probe) at 20-100 meters depth below the ground surface; ensure accurate calibration of the equipment and record temperature values at different depths; correct the surface thermal disturbance with meteorological data (such as annual mean temperature); use interpolation techniques to draw shallow geothermal isopleth map and calculate shallow geothermal gradient. 3. Deep well geothermal data (>100 meters) reflects the true temperature distribution of deep geothermal reservoir (such as bottom hole temperature and stratified temperature). Collection method: select representative well sites (covering target tectonic units) to ensure sufficient well depth (usually >500 meters); use deep well temperature measuring instrument or thermistor (such as high-precision digital temperature sensor) to measure temperature at different depths; record auxiliary information such as wellbore conditions (such as wellbore stability, fluid seepage state); construct deep geothermal field distribution map by interpolation or fitting method and calculate deep geothermal gradient.

[0035] Seismic data is used to infer the characteristics of underground structures (such as fault zones and thermal reservoir morphology) and thermal fluid migration. It is mainly obtained through active and passive seismic exploration: 1. Active seismic exploration, by artificially exciting seismic waves (such as explosives, controlled sources), recording the propagation characteristics of waves in underground rock layers (such as reflection coefficient, wave velocity). Collection method: place seismic source points (excitation interval covering main structures) and receivers (such as three-component geophones) in the target area; record the arrival time, amplitude and phase of reflected and refracted waves after exciting seismic waves; use inversion algorithms (such as full waveform inversion) to construct underground velocity structure model, identify thermal reservoir boundaries and fault zones. 2. Passive seismic exploration, monitoring the seismic waves excited by natural sources (such as microseisms, background noise), analyzing the dynamic characteristics of underground media (such as anisotropy, fluid permeability). Collection method: deploy high-sensitivity seismic network (covering main thermal display points in the region); record wave field data (such as frequency spectrum of background noise) excited by natural sources for a long time; infer the migration path of underground thermal fluid and the depth of thermal reservoir by wave velocity inversion or microseismic positioning technology.

[0036] Geological survey data reveals the tectonic background, lithological characteristics and thermal reservoir conditions of the geothermal region through field observation and data integration, including the following contents: 1. Tectonic and stratigraphic information, location, direction and scale of regional deep fault zone; lithology (such as sandstone, limestone, lava), thickness and permeability (such as sandstone permeability > 100 mD, volcanic clastic rock permeability < 10 mD). Collection method: Combine existing geological maps, profile maps and field geological observation (such as outcrop analysis) to verify the accuracy of stratigraphy and structure; supplement missing profile data to build 2D geological profile (such as covering the main thermal reservoir); input drilling data, stratigraphic information into modeling software (such as Petrel, GOCAD) to generate 3D geological model to visualize the structure of the thermal system. 2. Thermal display characteristics and mineral information, surface thermal display (such as hot springs, sinter, fumarole), fluid chemical characteristics (hot springs with temperature > 50℃, geothermal wells with temperature > 80℃); distribution of typical minerals (such as clay minerals, zeolite, calcareous silicate minerals) for inferring thermal reservoir evolution temperature (such as clay minerals indicating thermal reservoir temperature < 200℃, zeolite indicating > 200℃). Collection method: Systematic measurement of surface thermal display points (latitude, longitude, temperature, conductivity, pH, flow rate); collection of spring water samples (avoiding interference of mixing effect), testing of Na, K, Ca, Mg, SiO2, etc. components; analysis of mineral types (such as calcareous, siliceous) and alteration degree (such as chloritization, sericitization) to assist in thermal reservoir temperature estimation. 3. Geological hazards and heat source correlation, potential geological hazards (such as volcanic activity, landslide, flood risk) on geothermal system; heat source properties (such as magma chamber, radioactivity) and thermal reservoir supply conditions (such as groundwater circulation path). Collection method: Combine volcanic activity records (such as young igneous rock distribution since the Pliocene), slope stability analysis (such as ground penetrating radar detection); determine heat source type (magma or hydrothermal) through drilling data (such as bottom hole temperature, logging curve) and geochemical analysis (such as isotope ratio 18 O、 87 Sr / 86 Sr) .

[0037] Step 2: Determine whether the geothermal data meets the preset geothermal data condition. If the geothermal data meets the preset geothermal data condition, it is determined that the groundwater in the geothermal region is unstable.

[0038] In the embodiments of the present disclosure, the preset geothermal data condition includes: the geothermal temperature is higher than the average geothermal temperature of the background area of geological survey by more than 60℃, or the geothermal temperature is greater than the maximum value of the average geothermal temperature of all sites in the background area of geological survey, and the number of sites with geothermal temperature greater than the average geothermal temperature is less than half of the total number of local sites.

[0039] First, determine the geothermal temperature. If the geothermal temperature meets any of the following conditions, it is directly determined that the region is an unstable region, and no further seismic determination is needed:

[0040] Condition one: the ground temperature is higher than the average ground temperature of the geological survey background area by 60℃, that is:

[0041] T i > T avg + 60℃ (1)

[0042] wherein, T i : the ground temperature of the i-th site (target site);

[0043] T avg : the average ground temperature of the geological survey background area (average value of all background sites).

[0044] Condition two: the ground temperature is greater than the maximum value of the average ground temperature of the geological survey background area, and the number of sites with ground temperature greater than the average ground temperature is less than half of the total number of local sites, that is:

[0045] T i > max(T avg ) and M < N / 2 (2) max(T avg ): the maximum value of the average ground temperature of the geological survey background area (if there are multiple average value scenarios in the background area, such as regional statistics);

[0046] M: the number of sites in the local sites that satisfy T j > T avg (the number of sites with ground temperature higher than the background average value);

[0047] N: the total number of local sites.

[0048] In the embodiments of the present disclosure, the geothermal area refers to an area within a range of 20℃ higher than the ground temperature background value of the geological survey background area.

[0049] Step 3: If the ground temperature data does not meet the preset ground temperature data condition, it is determined whether the geothermal area has an earthquake, if there is an earthquake, it is determined whether the earthquake data meets the preset earthquake data condition, if the earthquake data meets the preset earthquake data condition, it is determined that the underground water of the geothermal area is unstable.

[0050] In the embodiments of the present disclosure, the preset earthquake data condition includes: the number of earthquakes in the area within one month is more than 20 times, and the average earthquake magnitude is less than or equal to 2.5, or the number of earthquakes in the area within one year is greater than or equal to 10 times, and the average earthquake magnitude is less than or equal to 3.5.

[0051] If the geothermal temperature does not meet any of the above conditions, it is determined whether there is an earthquake in the region, and if the earthquake meets any of the following conditions, it is directly determined that the region is an unstable region and does not need to be further determined by geological survey: Condition three: the number of earthquakes in the region within one month is more than 20 times, and the average magnitude of the earthquakes is less than or equal to 2.5; Condition four: the number of earthquakes in the region within one year is greater than or equal to 10 times, and the average magnitude of the earthquakes is less than or equal to 3.5.

[0052] That is:

[0053] The preset earthquake data condition can be expressed as the following logical formula:

[0054] (N month >20 and Mˉ month ≤2.5) or (N year ≥10 and Mˉ year ≤3.5) (3)

[0055] Wherein, N month : the number of earthquakes in the target region within one month (the statistical period is a natural month);

[0056] Mˉ month : the average magnitude of earthquakes in the target region within one month (the average magnitude of all earthquake events);

[0057] N year : the number of earthquakes in the target region within one year (the statistical period is a natural year);

[0058] Mˉ year : the average magnitude of earthquakes in the target region within one year (the average magnitude of all earthquake events).

[0059] Logical relationship explanation

[0060] Condition three: N month >20 and Mˉ month ≤2.5

[0061] The earthquake activity is frequent (more than 20 times) within one month, but the single earthquake magnitude is generally low (the average is less than or equal to 2.5).

[0062] Condition four: N year ≥10 and Mˉ year ≤3.5

[0063] The earthquake activity is relatively frequent (≥10 times) within one year, but the single earthquake magnitude is generally small (the average is less than or equal to 3.5).

[0064] Logical "or": meeting condition three or condition four is considered to meet the preset earthquake data condition.

[0065] For easy understanding, take an example: suppose the area data is:

[0066] The number of earthquakes in a month N month = 25, the average magnitude M month = 2.3;

[0067] The number of earthquakes in a year N year = 12, the average magnitude M year = 3.2.

[0068] Formula verification:

[0069] Condition three is established (25>20 and 2.3≤2.5);

[0070] Condition four is established (12≥10 and 3.2≤3.5);

[0071] Since condition three or condition four is established, the area meets the preset earthquake data condition.

[0072] In the embodiments of the present disclosure, the earthquake data is earthquake data in a high earthquake period obtained according to the number of earthquakes, the magnitude of earthquakes and the position of earthquakes.

[0073] The earthquake specifically refers to a measured microseismic event.

[0074] In the high earthquake period, the number of earthquakes within 5-30 kilometers, the magnitude of earthquakes and the position of earthquakes in a year are taken as the standards, and the data obtained by statistics is used to judge the earthquake situation.

[0075] The method for counting the number of earthquakes, the magnitude of earthquakes and the position of earthquakes is specifically: in the high earthquake period, within 50 kilometers from the earthquake source, the number of times of earthquakes with a magnitude less than or equal to 2.5 in a month is selected and recorded as the monthly earthquake number; the number of times of earthquakes with a magnitude less than or equal to 3.5 in a year is selected and recorded as the annual earthquake number; the average magnitude of earthquakes in a year is selected and recorded as the annual average earthquake magnitude.

[0076] The high earthquake period refers to a period of time in which the average magnitude of earthquakes in a high earthquake year is greater than or equal to 2.5, the number of earthquakes is greater than or equal to 7 times, the average magnitude of earthquakes is greater than the average magnitude of earthquakes in the previous year by more than 0.5, the average period of earthquakes is less than one year, and the distribution range of the epicenter of earthquakes is less than 50 kilometers, and this period of time is defined as the high earthquake period.

[0077] The specific steps of the method for determining the high earthquake period are as follows: 1) determining the baseline parameters of the high earthquake area in the year before the earthquake: the high earthquake period is 2 years, and the earthquakes in the high earthquake area are screened at two time points, i.e., one year before the high earthquake and one month before the start of the high earthquake, the average magnitude of the screened earthquakes one year before the high earthquake is recorded as the one-year baseline magnitude average, and the average period of the screened earthquakes one month before the start of the high earthquake is recorded as the one-month baseline period, the latitude, longitude and altitude of the earthquake in the high earthquake area are measured and recorded as the area latitude, area longitude and area altitude, the shallow source earthquake epicenter and the deepest source point in the high earthquake area are measured and recorded as the depth epicenter and the depth deepest source point; 2) determining the characteristic parameters of the high earthquake area in the high earthquake year: the high earthquake period is one year, and the earthquakes in the high earthquake area are screened at two time points, i.e., one month before the high earthquake year and half a year before the start of the earthquake, the average magnitude of the screened earthquakes one month before the high earthquake year is recorded as the half-year baseline magnitude average, and the average period of the screened earthquakes half a year before the high earthquake year is recorded as the half-year baseline period, the magnitude difference average and the period difference average of the high earthquake year compared with the previous year are calculated, the magnitude difference average = the one-year baseline magnitude average of the high earthquake year - the baseline magnitude average of the previous year, the period difference average = the one-month baseline period of the high earthquake year - the one-month baseline period of the previous year, the earthquake magnitude of the high earthquake year is recorded as the one-time earthquake average magnitude, the sum of the earthquake magnitudes in a seismic sequence is divided by the number of earthquakes in the seismic sequence, wherein a seismic sequence refers to all earthquakes occurring within 4-10 time ranges from the start of the earthquake to the calm period of the earthquake during the high earthquake period; the earthquake period of the high earthquake year is the sum of the time intervals between the end of the earthquake and the start of the next earthquake, and the sum of the earthquake periods in a seismic sequence is divided by the number of earthquakes in the seismic sequence; 3) determining whether the high earthquake area has occurred a high earthquake: the high earthquake area needs to meet the following conditions: condition one: the magnitude difference average is greater than or equal to 0.5 magnitude; condition two: the half-year baseline period is greater than the one-month baseline period.

[0078] At the start of the high earthquake period, the number of earthquakes, earthquake magnitude and earthquake location within one month are screened and recorded as the one-month earthquake number, the one-month earthquake magnitude average and the one-month earthquake epicenter location, respectively; the number of earthquakes, earthquake magnitude and earthquake location within half a year are screened and recorded as the half-year earthquake number, the half-year earthquake magnitude average and the half-year earthquake epicenter location, respectively; the number of earthquakes, earthquake magnitude and earthquake location within one year are screened and recorded as the one-year earthquake number, the one-year earthquake magnitude average and the one-year earthquake epicenter location, respectively.

[0079] During the high earthquake period, a one-time high earthquake event specifically refers to: the number of earthquakes in the high earthquake area within one month is greater than or equal to 20 times, the one-month earthquake magnitude average is less than or equal to 2.5 magnitude, and the distribution range of the earthquake epicenter in the one-time high earthquake event is within 10 KM.

[0080] In the earthquake high incidence period, the secondary earthquake high incidence event specifically refers to that the number of earthquakes in an earthquake high incidence area in a year is greater than or equal to 10, the average magnitude of earthquakes in a year is less than or equal to 3.5, and the distribution range of the earthquake epicenter in the secondary earthquake high incidence event is within 30 km.

[0081] The primary earthquake high incidence event and the secondary earthquake high incidence event both occur in the earthquake high incidence area.

[0082] The earthquake high incidence area is the next period of the earthquake high incidence period in the earthquake sequence, that is, the earthquake recovery period. Specifically, after the earthquake high incidence period, less than 1 earthquake occurs within 2 days after a primary earthquake, or the magnitude of the earthquake is greater than that of the previous earthquake by 1. The third earthquake high incidence event specifically refers to that, after a primary earthquake, 2 earthquakes or earthquakes with a magnitude difference of more than 3 occur within 50 km from the earthquake location.

[0083] If multiple earthquake high incidence events occur in the earthquake high incidence period, the number of earthquakes, the magnitude of earthquakes, and the location of earthquakes in the earthquake high incidence period are analyzed, and are recorded as the number of earthquakes in the earthquake high incidence period, the average magnitude of primary earthquakes in the earthquake high incidence period, and the epicenter location of primary earthquakes in the earthquake high incidence period, respectively.

[0084] Step 4: If the earthquake data does not meet the preset earthquake data condition, it is determined whether the geological survey data meets a preset geological survey data condition. If the geological survey data meets the preset geological survey data condition, it is determined that the underground water in the geothermal area is unstable.

[0085] In the embodiment of the present disclosure, the preset geological survey data condition includes that the rock structure in the geological survey data has an unstable stratum.

[0086] In the embodiment of the present disclosure, the geological survey data refers to rock structure data obtained when a geological survey is performed. The geological survey includes taking a drill hole profile of a specific length on both sides along the strike of an earthquake peak value line with the earthquake peak value line as an axis, and obtaining geological data at the profile.

[0087] If the geothermal temperature does not meet any of the above conditions, and the earthquake does not meet any of the above conditions, a geological survey is performed to determine the rock structure. According to the rock structure, it is determined whether the area has an unstable stratum. If it does, it is determined to be an unstable area.

[0088] When performing a geological survey to determine the rock structure, a drill hole profile of 1 km long is taken on both sides along the strike of the earthquake peak value line with the earthquake peak value line as an axis. If there is an unstable stratum, it is determined to be an unstable area.

[0089] The geological survey specifically refers to: taking a 1km long borehole profile directly above a seismic peak value line, and if there is an unstable stratum, determining it as an unstable area.

[0090] The method for determining the structure of the geological survey includes but is not limited to: borehole detection, geological sketch, borehole photography and measurement.

[0091] To sum up, the underground water stability judgment method provided in the embodiment realizes the accurate identification of the underground water stability risk, reduces the misjudgment rate, improves the judgment efficiency, and enhances the risk early warning capability.

[0092] Embodiment two

[0093] Figure 2 The structural diagram of an underground water stability judgment system provided according to an embodiment of the present application is shown in Figure 2 The system comprises:

[0094] The data acquisition module 100 is configured to acquire the geothermal region's geothermal data, seismic data and geological survey data.

[0095] The geothermal condition judgment module 200 is configured to judge whether the geothermal data meets the preset geothermal data condition, and if the geothermal data meets the preset geothermal data condition, determine that the underground water of the geothermal region is unstable.

[0096] The preset geothermal data condition includes: the geothermal temperature is higher than the average geothermal temperature of the background area by more than 60℃, or the geothermal temperature is greater than the maximum value of the average geothermal temperature of all stations in the background area, and the number of stations with geothermal temperature greater than the average geothermal temperature is less than half of the total number of stations in the local area.

[0097] The geothermal region refers to a region within a range of 20℃ higher than the background geothermal temperature of the background area.

[0098] The seismic condition judgment module 300 is configured to judge whether there is an earthquake in the geothermal region if the geothermal data does not meet the preset geothermal data condition, and if there is an earthquake, judge whether the seismic data meets the preset seismic data condition, and if the seismic data meets the preset seismic data condition, determine that the underground water of the geothermal region is unstable.

[0099] The preset seismic data condition includes: the number of earthquakes in the region within one month is more than 20 times, and the average magnitude of the earthquakes is less than or equal to 2.5, or the number of earthquakes in the region within one year is greater than or equal to 10 times, and the average magnitude of the earthquakes is less than or equal to 3.5.

[0100] The earthquake data is earthquake data during a high-incidence period of earthquakes, which is obtained according to the number of earthquakes, the magnitude of earthquakes and the location of earthquakes.

[0101] The geological condition judging module 400 is configured to judge whether the geological survey data meets a preset geological survey data condition if the earthquake data does not meet the preset earthquake data condition, and determine that the underground water in the geothermal region is unstable if the geological survey data meets the preset geological survey data condition.

[0102] The preset geological survey data condition includes that an unstable stratum exists in the rock stratum structure in the geological survey data.

[0103] The geological survey data refers to rock stratum structure data obtained during geological survey, and the geological survey includes: taking a seismic peak line as an axis, taking a drill hole profile of a specific length on both sides of a seismic strike, and obtaining geological data at the profile.

[0104] The third aspect of the present application provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method of the first aspect.

[0105] In summary, the underground water stability judging system provided by the present embodiment realizes accurate identification of the stability risk of underground water by using the triple-layered judging logic of the geothermal anomaly characteristics, the seismic activity mode and the geological structure verification, reduces the misjudgment rate, improves the judging efficiency, and enhances the risk early warning capability.

[0106] Embodiment three

[0107] To achieve the above-mentioned embodiments, the present disclosure further provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method of embodiment one.

[0108] Embodiment four

[0109] To achieve the above-mentioned embodiments, the present disclosure further provides a computer readable storage medium having a computer program stored thereon, wherein the program is executed by a processor to implement the method of embodiment one.

[0110] In the description of the application, the description of the terms "one embodiment", "some embodiments", "an example", "a specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the application. The illustrative description of the above terms in the specification does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any appropriate manner in any one or more embodiments or examples. In addition, different embodiments or examples described in the specification and the features of different embodiments or examples can be combined and combined by those skilled in the art without contradiction, if any.

[0111] Any process or method descriptions or descriptions of the flow diagrams in the flow charts described herein or otherwise described in this specification can be understood as representing the steps of a method or process, including one or more steps for implementing custom logic functions or processes, and the scope of the preferred embodiments of the present application includes additional implementation involving other steps, which can be performed at substantially the same time or in reverse order or in other order, according to the functions involved, which should be understood by those skilled in the art.

[0112] Although the embodiments of the application have been shown and described above, it should be understood that the above-described embodiments are exemplary and should not be construed as limiting the application, and those skilled in the art can make changes, modifications, replacements and variations to the above-described embodiments within the scope of the application.

Claims

1. A method for judging groundwater stability, characterized in that, The method includes: Acquire geothermal data, seismic data, and geological survey data for geothermal areas; Determine whether the ground temperature data meets the preset ground temperature data conditions. If the ground temperature data meets the preset ground temperature data conditions, then determine that the groundwater in the geothermal area is unstable. If the geothermal data does not meet the preset geothermal data conditions, it is determined whether there is an earthquake in the geothermal area. If there is an earthquake, it is determined whether the earthquake data meets the preset earthquake data conditions. If the earthquake data meets the preset earthquake data conditions, it is determined that the groundwater in the geothermal area is unstable. If the earthquake data does not meet the preset earthquake data conditions, then it is determined whether the geological survey data meets the preset geological survey data conditions. If the geological survey data meets the preset geological survey data conditions, then it is determined that the groundwater in the geothermal area is unstable.

2. The groundwater stability assessment method according to claim 1, characterized in that, The preset ground temperature data conditions include: the ground temperature is more than 60°C higher than the average ground temperature of the geological survey background area, or the ground temperature is greater than the maximum value of the average ground temperature of all stations in the geological survey background area, and the number of stations with ground temperatures greater than the average ground temperature is less than half of the total number of local stations.

3. The groundwater stability assessment method according to claim 1, characterized in that, The preset earthquake data conditions include: the number of earthquakes in the region within a month exceeds 20, and the average earthquake magnitude is less than or equal to 2.5; or the number of earthquakes in the region within a year is greater than or equal to 10, and the average earthquake magnitude is less than or equal to 3.

5.

4. The groundwater stability assessment method according to claim 1, characterized in that, The preset geological survey data conditions include: the geological survey data contains unstable strata.

5. The method for judging groundwater stability according to claim 1, characterized in that, The geothermal area refers to the region where the geothermal temperature is 20°C higher than the background value of the geological survey area.

6. The method for judging groundwater stability according to claim 1, characterized in that, The earthquake data mentioned are earthquake data from periods of high earthquake frequency, obtained by statistical analysis of earthquake frequency, magnitude, and location.

7. The method for judging groundwater stability according to claim 1, characterized in that, The geological survey data refers to the rock strata structure data obtained during the geological survey. The geological survey includes: taking borehole profiles of a specific length on both sides of the earthquake strike line as the axis, and obtaining geological data at the profiles.

8. A groundwater stability assessment system, characterized in that, The system includes: The data acquisition module is used to acquire geothermal data, seismic data, and geological survey data for the geothermal area. The geothermal condition judgment module is used to determine whether the geothermal data meets the preset geothermal data conditions. If the geothermal data meets the preset geothermal data conditions, the groundwater in the geothermal area is determined to be unstable. The earthquake condition judgment module is used to determine whether an earthquake exists in the geothermal area if the ground temperature data does not meet the preset ground temperature data conditions. If an earthquake exists, it determines whether the earthquake data meets the preset earthquake data conditions. If the earthquake data meets the preset earthquake data conditions, it determines that the groundwater in the geothermal area is unstable. The geological condition judgment module is used to determine whether the geological survey data meets the preset geological survey data conditions if the seismic data does not meet the preset seismic data conditions. If the geological survey data meets the preset geological survey data conditions, the module determines that the groundwater in the geothermal area is unstable.

9. An electronic device, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the method as described in any one of claims 1-7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the program is executed by the processor, it implements the method as described in any one of claims 1-7.

Images