CN117371267A - Mining area water inrush quantitative prediction method based on ground transient electromagnetic - Google Patents

Abstract

The invention discloses a quantitative prediction method for water inrush in mining areas based on ground transient electromagnetism, and relates to the technical field of quantitative prediction of water inrush in mining areas. The resistivity after air-drying and the resistivity when different amounts of water are injected into the air-dried rock, and then the resistivity of the rock under different water contents is obtained. Based on the classic Archie formula, a quantitative relationship between rock resistivity and water content is established; a three-dimensional vector is used The finite element transient electromagnetic simulation algorithm analyzes and calculates the actual complex terrain model established to obtain the influence of complex terrain on the transient electromagnetic response of the ground; and then establishes a quantitative prediction model of mine water inrush based on ground transient electromagnetic to predict the water content of the mining area. Quantitative prediction of layer water inrush; this method establishes a qualitative relationship between rock resistivity and water content, and conducts accurate quantitative analysis and prediction of rock layer water content.

Description

A quantitative prediction method for water inrush in mining areas based on ground transient electromagnetism

Technical field

The invention relates to the technical field of quantitative prediction of water inrush in mining areas, and specifically relates to a method for quantitative prediction of water inrush in mining areas based on ground transient electromagnetism.

Background technique

Water damage is second only to gas explosion in serious and serious accidents in coal mines. Water-damaged mine disasters are not only likely to cause heavy casualties to underground workers, but are also the most prominent among coal mine disasters in terms of the severity of economic losses, difficulty in accident rescue and rescue, and the time required to restore mine production. In recent years, with the improvement and development of water prevention and control technology, and under the strong supervision of the state and local governments, although the number of water inrush accidents in coal mines has shown a downward trend year after year, it has not yet been fundamentally and effectively controlled.

As a "vanguard", geophysical hydrogeological detection technology has achieved certain application results in preventing and controlling water in coal mines, but it still cannot meet the technical requirements of preventing and controlling water in coal mines. Its theory, methods, technology and applications need further innovation and development to ensure that Achieve accurate detection of water inrush hazard sources, effectively solve the problem of water inrush disasters, improve the technical level of water inrush control in coal mines, reduce the economic cost of water inrush control in coal mines, and improve coal mine safety and economic benefits.

At present, the electrical and electromagnetic geophysical methods mainly used in coalfield hydrogeological exploration are all based on the relative magnitude of underground resistivity values to infer and explain water-rich anomalies. They are indirect detection of water-bearing properties and cannot quantitatively evaluate the water content of rock formations. , the fundamental reasons are: ① The qualitative relationship between rock resistivity and water content has not been established; ② The measurement signal is affected by many factors such as terrain, instrument stability, human interference, etc., and there is a large deviation from the true response of underground rocks; ③ Geophysics Multiple interpretations of inversion.

Contents of the invention

In view of the shortcomings of the existing technologies that are all indirect detection of water-bearing properties and cannot quantitatively evaluate the water content of rock formations, the present invention provides a quantitative prediction method for water inrush in mining areas based on ground transient electromagnetism, by obtaining rock resistivity and moisture content. relationship and the influence of complex terrain conditions on the response of ground transient electromagnetic signals. Establish a quantitative prediction theory and method for mine water inrush based on ground transient electromagnetics to achieve a quantitative evaluation of the water filling capacity of aquifers, thereby solving the problem of existing technology that cannot carry out rock formation analysis. Problems with quantitative evaluation of moisture content.

A quantitative prediction method for water inrush in mining areas based on ground transient electromagnetics, including the following steps:

Collect cores from the working face area of the mining area;

The three-dimensional high-density resistivity method is used to measure the resistivity of the water-bearing section of the rock core after it is air-dried and the resistivity of the air-dried core when different amounts of water are injected into it, and then the resistivity of the core under different water contents is obtained;

Based on the classic Archie formula and based on the resistivity of rock cores under different water contents, a qualitative relationship between the resistivity and water content of the rock core is established;

Establish an actual complex terrain model, use a three-dimensional vector finite element transient electromagnetic simulation algorithm to analyze and calculate the actual complex terrain model, and obtain the ground transient electromagnetic field response of electrical sources and loop sources under various undulating terrain conditions, thus obtaining Find out the influence of complex terrain on the transient electromagnetic response of the ground;

Based on the seismic exploration data in the working surface area, Geoview three-dimensional inversion software is used to conduct wave impedance inversion, and then construct a geophysical model;

According to the influence of complex terrain on the transient electromagnetic response of the ground, terrain correction is performed on the transient electromagnetic data of each measuring point in the working surface area;

Using the geophysical model as the initial model, the terrain-corrected transient electromagnetic data are inverted and simulated to obtain the underground resistivity distribution;

According to the qualitative relationship between the resistivity and water content of the rock core, the underground resistivity distribution obtained by inversion is used to calculate the water content of the aquifer in the mining area to be measured, thereby quantitatively predicting the water inrush of the aquifer in the mining area.

Further, it also includes using a core nuclear magnetic resonance imaging analyzer after using the three-dimensional high-density resistivity method to measure the resistivity of the water-bearing section of the rock core after air-drying and the resistivity of the air-dried core when different amounts of water are injected. Test the porosity of the cores in all aquifer sections, and calculate the water content of the aquifer in the working surface area based on the qualitative relationship between porosity and the resistivity and water content of the cores.

Further, the resistivity of the core under different water contents was obtained by using water injection displacement method to sequentially increase the amount of water injected into the air-dried core, and the three-dimensional high-density resistivity method was used to measure the resistance of the core under different water contents. Rate.

Furthermore, the actual complex terrain model was established by using DEM elevation data in combination with COMSOL software.

Furthermore, based on the classic Archie formula, according to the resistivity of the core under different water contents, a qualitative relationship between the resistivity of the core and the water content was established. The specific steps include:

Archie's formula is used to fit the rock electrical test data of rock samples, and the lithology coefficient a and cementation index m are obtained:

In the formula, F is the formation factor, R ₀ is the resistivity of the rock sample containing 100% water, Ω·m; R _w is the resistivity of the water contained in the rock sample, Ω·m; φ is the porosity;

Archie's formula is then used to fit the rock electrical test data of the rock sample to obtain the lithology coefficient b and saturation index n:

In the formula, I is the resistivity increase coefficient; R ₀ is the resistivity of the 100% water-containing rock sample, Ω·m; R _t is the resistivity of the rock sample when it is partially hydrated, Ω·m; S _w is the water saturation of the rock. ;

And then get:

Since a and b are both lithology coefficients, let ab=x, we get:

Use the test data to fit x, m, n to obtain the quantitative relationship between resistivity and water.

Further, the geoview three-dimensional inversion software is used to perform wave impedance inversion based on the seismic exploration data in the working surface area, thereby constructing a geophysical model; which includes the following steps:

Geoview three-dimensional inversion software is used to conduct wave impedance inversion and construct a geological model beneath the survey area;

The same attribute units in the geological model are the same geological lithology. Combined with the resistivity logging data, the same geological lithology is given the same electrical properties to construct a geophysical model.

Furthermore, using the geophysical model as the initial model, the corrected transient electromagnetic data are inverted and simulated to obtain the underground resistivity distribution. The calculation formula is:

Among them, d _c is the theoretical value calculated by the current model, the subscript i is the i-th data, M represents the total number of observation data, m represents the model parameters, j in the table below represents the j-th model parameter, and N ₁ represents the model parameter. The total number, P represents the model parameters after merging the same attributes, N ₂ represents the total number of model parameters after merging the same attributes, is the partial derivative operator.

Furthermore, in the process of inverting and simulating the corrected transient electromagnetic data using the geophysical model as the initial model, the same geological and lithological units are used as the same inversion unit.

Furthermore, a quantitative prediction system for water inrush in mining areas based on ground transient electromagnetics includes:

The collection module is used to collect cores from the working face area of the mining area;

The test module is used to use the three-dimensional high-density resistivity method to measure the resistivity of the water-bearing section of the rock core after it is air-dried and the resistivity of the air-dried core when different amounts of water are injected into it, so as to obtain the resistivity of the rock core under different water contents;

The module for establishing the relationship between core resistivity and water content is used to establish a qualitative relationship between the core resistivity and water content based on the classic Archie formula and based on the resistivity of the core under different water contents;

The analysis module is used to establish actual complex terrain models. It uses a three-dimensional vector finite element transient electromagnetic simulation algorithm to analyze and calculate the actual complex terrain models, and obtains ground transient electromagnetic signals of electrical sources and loop sources under various undulating terrain conditions. Magnetic field response, thereby obtaining the influence of complex terrain on the transient electromagnetic response of the ground;

The geophysical model building module is used to perform wave impedance inversion using Geoview three-dimensional inversion software based on seismic exploration data in the working surface area, and then construct a geophysical model;

The correction module is used to perform terrain correction on the transient electromagnetic data of each measuring point in the working surface area based on the influence of complex terrain on the transient electromagnetic response of the ground;

The resistivity distribution acquisition module is used to perform inversion simulation on the terrain-corrected transient electromagnetic data using the geophysical model as the initial model to obtain the underground resistivity distribution;

The prediction module is used to calculate the aquifer water content of the mining area to be measured based on the qualitative relationship between the resistivity and water content of the core and the underground resistivity distribution obtained by inversion, thereby quantitatively predicting the water inrush of the aquifer in the mining area.

The present invention provides a quantitative prediction method for water inrush in mining areas based on ground transient electromagnetism, which has the following features:

Beneficial effects:

This invention is based on the classic Archie formula, and based on the resistivity of rock cores under different water contents, the relationship between rock resistivity and water content is obtained, and the actual complex terrain model is analyzed and calculated by using a three-dimensional vector finite element transient electromagnetic simulation algorithm. The influence rules of complex terrain conditions on the ground transient electromagnetic signal response are obtained. According to the influence rules of complex terrain conditions on the ground transient electromagnetic response, terrain correction is performed on the high-precision transient electromagnetic data of each measuring point in the working surface area. The terrain-corrected transient electromagnetic data are Conduct inversion simulation of variable electromagnetic data, and then establish a theory and method for quantitative prediction of water inrush in mines based on ground transient electromagnetics, achieving a quantitative evaluation error of less than 30% for aquifer water filling.

Description of the drawings

Figure 1 is a flow chart of the quantitative prediction method for water inrush in mining areas based on ground transient electromagnetics according to the present invention;

Figure 2 is a flow chart for calculating the quantitative relationship between rock resistivity and water content in the embodiment of the present invention.

Detailed ways

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some of the embodiments of the present invention, rather than all the embodiments.

As shown in Figure 1, this invention focuses on the goal of "quantitative evaluation of aquifer water filling" and uses rock physics, numerical simulation, field testing and other methods to study the relationship between rock resistivity and water content, complex terrain conditions The high-precision joint inversion method of transient electromagnetic signals under the influencing rules and constrained conditions of ground transient electromagnetic signals is then established to establish a theory and method for quantitative prediction of water inrush in mines based on ground transient electromagnetics, achieving a quantitative evaluation error of less than 30% for aquifer water filling. %.

(1) Coal mine hydrogeological exploration technology: In the past, coal mine hydrogeological exploration methods were relatively simple, mainly based on well logging and pumping tests. At present, remote sensing hydrogeological survey, geophysical prospecting, hydrogeological drilling, pumping test, water chemistry and environmental isotope analysis, groundwater flow field analysis, rock and soil physical and mechanical analysis and other means are used to carry out comprehensive exploration from multiple aspects, angles and levels. Study water-filled conditions. Firstly, remote sensing hydrogeological survey and ground geophysical prospecting were carried out to understand the replenishment and drainage conditions of the underground (surface) water system and the location of the well field in the system, and to initially determine the water richness of the water-filled aquifer. Then a series of work such as drilling and pumping tests will be carried out. Based on the drilling coring results and combined with the coalfield geological exploration drilling data, we draw a comparison map of the aquifer structure, aquifer thickness contour map, and a contour map of the ratio of brittle rock to plastic rock. Through these three maps, the well field The lithology and lithofacies structure are clearly displayed. The greater the thickness of the aquifer and the greater the ratio of brittle-to-plastic rock, the stronger the water-richness of the aquifer. In the past, water sample testing and water level observation were only used to establish water chemistry background values and understand groundwater levels, but were not used to conduct analysis and research on water filling conditions. At present, on the basis of sound basic data, the hydrochemical field and groundwater flow field information are fully analyzed and studied, and a series of hydrochemical diagrams and flow field diagrams are drawn to analyze groundwater replenishment and drainage conditions.

(2) Water inrush evaluation: For the prediction of mine water inflow, analytical and analogous methods have been used in the past. The prediction results are relatively reliable for mines with simple conditions or Carboniferous-Permian coal seams. However, for mines with complex conditions, the calculation results deviate greatly. At present, the three-dimensional numerical simulation method of groundwater is mainly used to predict the amount of water inflow and drainage volume in mines under mining conditions. By using Modflow software to depict a three-dimensional visualization model of hydrogeological concepts, you can intuitively have a deeper understanding of the hydrogeological conditions of the mining area from a three-dimensional perspective. This method can better solve complex problems such as heterogeneous anisotropy. A large amount of geoscientific information has been obtained through various exploration methods, which provides an important basis for the evaluation and prediction of the water-richness and water inrush risk of water-filled aquifers. However, this information is partially and fragmented. If it is applied directly and fragmented to analyze water damage, it will be very superficial and cannot give full play to the overall effect and the value contained in the exploration data. In recent years, aquifer water richness evaluation methods and mine water inrush risk prediction methods based on multi-source geoscientific information composite overlay technology have been developed. The roof water damage adopts the "three pictures-double prediction method", and the floor water damage is based on the "water inrush coefficient" evaluation method. Based on the "vulnerability index method". The principle of these advanced evaluation technologies is to use mathematical models to systematically analyze, count, simulate and superimpose a large amount of geoscientific information obtained by various exploration methods, so as to obtain an overall evaluation result that meets actual conditions. Summarize the development rules of water damage assessment and prediction technology: from qualitative to quantitative, from single factor to multiple factors, from homogeneous, planar to heterogeneous three-dimensional visualization.

(3) Hydrogeological and geophysical detection technology: Hydrogeological exploration first uses geophysical means to conduct all-round detection, then uses drilling means to verify the indirect detection results of geophysical prospecting, and finally conducts water drainage and treatment. Geophysical prospecting is an advanced and macro-control method and assumes very important responsibilities. After years of basic theoretical research, field testing and analysis and comparison, a relatively mature and advanced detection technology has been formed. From the perspective of time and space, a ground-well-lane joint fine exploration technology has been formed. First, ground geophysical exploration is performed to delineate abnormal locations, and then underground geophysical exploration is performed to get close to the target for fine detection. Finally, drilling and data are used to verify anomalies. At present, transient electromagnetic field has become the main technical means for hydrogeological exploration in coalfield mining areas.

To sum up, the current electromagnetic and electromagnetic geophysical methods mainly used in coalfield hydrogeological exploration are all based on the relative size of underground resistivity values to infer and explain water-rich anomalies. They are all indirect detection of water-bearing properties and cannot detect rock formations. The fundamental reasons for quantitative evaluation of water content: ① The qualitative relationship between rock resistivity and water content has not been established; ② The measurement signal is affected by many factors such as terrain, instrument stability, human interference, etc., and there is a large deviation from the true response of underground rocks; ③ Multiple solutions in geophysical inversion. Therefore, to implement quantitative water inrush prediction technology based on geophysics, the above three problems must be solved. To this end, the present invention intends to carry out the following research:

(1) Qualitative relationship between rock resistivity and rock water content: Based on Archie's formula, the changes in rock resistivity values under different water content conditions are studied.

(2) The influence of complex terrain on transient electromagnetic signals: Based on the current coalfield transient electromagnetic hydrogeological exploration technology, the response characteristics of transient electromagnetic signals of loop sources and electrical sources under complex qualitative conditions are studied in order to obtain the real rock formation instantaneous signals. Provides the basis for variable electromagnetic response.

(3) Ground transient electromagnetic high-precision constrained inversion method: Combined with drilling data and seismic data, the ground transient electromagnetic multi-parameter high-precision constrained inversion method was developed to solve the longitudinal positioning, thickness and true resistance of the electrical layer. Accurate explanation.

Specifically, it includes the following steps:

(1) Select a typical coal mine as the test mining area, drill and core all strata in the representative working face area, and catalog them to provide typical cores for petrophysical experiments.

(2) Perform petrophysical testing indoors. First, use a three-dimensional high-density core resistivity measurement system to test the resistivity of all cores obtained; then use the MesoMR core nuclear magnetic resonance imaging analyzer to test the pores of all water-bearing section cores. degree and moisture content; then air-dry the cores of the water-bearing sections to eliminate all water, and then use a three-dimensional high-density core resistivity measurement system to test the resistivity of the air-dried rocks without water, and establish a rock resistivity database for each layer.

(3) Use the water injection method to displace the air-dried rocks, gradually increase the amount of water injected, and conduct resistivity measurements to obtain the core resistivity under different water contents. Based on the classic Archie formula, polynomial fitting technology is used to establish the rock resistivity. The quantitative relationship with water content is shown in Figure 2, specifically including:

First, Archie’s formula (1) is used to fit the rock electrical experimental data of rock samples to obtain a and m.

In the formula, R ₀ is the resistivity of the rock sample containing 100% water, Ω·m; R _w is the resistivity of the water contained in the rock sample, Ω·m; φ is the porosity; a is the lithology coefficient; m is the cementation index.

Archie's formula (2) is then used to fit the rock electrical experimental data of the rock sample to obtain b and n.

In the formula, I is the resistivity increase coefficient; R ₀ is the resistivity of the 100% water-containing rock sample, Ω·m; R _t is the resistivity of the rock sample when it is partially hydrated, Ω·m; S _w is the water saturation of the rock. ; b is the lithology coefficient; n is the saturation index.

Multiply equation (1) by equation (2):

Since a and b are both lithology coefficients, let ab=x, we get:

In actual operation, experimental data can be used to fit x, m, n to obtain the quantitative relationship between resistivity and water.

(4) Use DEM elevation data combined with COMSOL software to establish an actual complex terrain model, and use a three-dimensional vector finite element transient electromagnetic simulation algorithm to calculate the ground transient electromagnetic field response of electrical sources and loop sources under various undulating terrain conditions. , summarizes the influence of complex terrain on the transient electromagnetic response of the ground, and provides a basis for the correction of terrain effects and the acquisition of real transient electromagnetic signals from underground rock formations.

(5) Retrieve the existing seismic data and well logging data in the test area, use the mature software Geoview three-dimensional inversion software to perform wave impedance inversion on the processed seismic exploration data, build a geological model below the test area, and invert the The same attribute units have the same geological lithology. Combined with the resistivity logging data, the same geological lithology is given the same electrical properties to construct a three-dimensional geoelectric model.

(6) Carry out high-precision transient electromagnetic data collection in the measurement area, record the elevation of each measurement point, and perform terrain correction on each measurement point data based on the research results in (4).

(7) Taking the geophysical model constructed in (5) as the initial model, perform a forward simulation on the transient electromagnetic data collected in (6) and corrected by terrain. During the simulation process, the same geological lithology unit is used as the initial model. The same inversion unit is used to reduce the number of inversion model units (as shown below), improve the inversion efficiency, and maintain a higher resolution of the electrical units below the measurement area to control the electrical boundaries and reduce the impact of the electromagnetic data volume effect. Obtain a more accurate underground resistivity distribution.

(8) Combined with the relationship between underground rock porosity, water content and resistivity obtained in (2), using the resistivity obtained by the inversion in (7), and on the basis of delineating the water-rich area of the aquifer, the calculation is The water volume of the aquifer within the measurement area is measured, and the geological model constructed based on the seismic inversion results enables visualization of the spatial distribution of water volume within the measurement area.

(9) Select representative aquifers, drill to drain them and record the water inflow.

(10) Compare the actual water inflow volume with the predicted water volume to verify the prediction error.

(11) Establish theory and technology for quantitative prediction of water inrush in mines based on ground transient electromagnetism, and achieve a quantitative prediction error of less than 30% for water content in water-filled aquifers.

The above are only preferred specific embodiments of the present invention, but the protection scope of the present invention is not limited thereto. Any person familiar with the technical field can, within the technical scope disclosed in the present invention, implement the technical solutions of the present invention. Equivalent substitutions or changes of the inventive concept thereof shall be included in the protection scope of the present invention.

Claims (9)

1. The mining area water inrush quantitative prediction method based on the ground transient electromagnetic is characterized by comprising the following steps of:

collecting cores of working surface areas of mining areas;

measuring the resistivity of the core after the water-containing layer section of the core is air-dried by adopting a three-dimensional high-density resistivity method, and injecting different water quantities into the air-dried core, so as to obtain the resistivity of the core under different water contents;

based on a classical Archie formula, establishing a qualitative relation between the resistivity and the water content of the core according to the resistivity of the core under different water contents;

establishing an actual complex terrain model, and analyzing and calculating the actual complex terrain model by adopting a three-dimensional vector finite element transient electromagnetic simulation algorithm to obtain the ground transient electromagnetic field responses of an electric source and a loop source under various undulating terrain conditions, thereby obtaining the influence rule of the complex terrain on the ground transient electromagnetic response;

performing wave impedance inversion by using Geoview three-dimensional inversion software according to the seismic exploration data of the working surface area, and further constructing a geophysical model;

carrying out terrain correction on transient electromagnetic data of each measuring point in the working surface area according to the influence rule of complex terrain on the transient electromagnetic response of the ground;

inversion simulation is carried out on the transient electromagnetic data corrected by the terrain by taking the geophysical model as an initial model, so that the underground resistivity distribution is obtained;

according to the qualitative relation between the resistivity and the water content of the rock core, the water-bearing layer water quantity of the mining area to be detected is calculated by adopting the underground resistivity distribution obtained by inversion, so that quantitative prediction is carried out on the water-bearing layer water burst of the mining area.

2. The mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 1, further comprising measuring the resistivity of the core water-containing interval after air-drying and the resistivity of the core after air-drying when injecting different water amounts into the core after air-drying by using a three-dimensional high-density resistivity method, measuring the porosities of the cores of all water-containing intervals by using a core nuclear magnetic resonance imaging analyzer, and calculating the water amount of the water-containing layer in the working surface area according to the porosities and the qualitative relation between the resistivity of the cores and the water content.

3. The mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 1, wherein the resistivity of the core with different water contents is obtained by sequentially increasing the water injection amount in the air-dried core by adopting a water injection displacement mode, and the resistivity of the core with different water contents is measured by a three-dimensional high-density resistivity method.

4. A mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 1, characterized in that the actual complex terrain model is built by using DEM elevation data in combination with COMSOL software.

5. The mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 1, wherein the establishing a qualitative relation between the resistivity of the core and the water content based on the classical alchi formula according to the resistivity of the core with different water contents comprises the following specific steps:

fitting the rock electrical measurement data of the rock sample by using an Archie formula, and solving a lithology coefficient a and a cementation index m:

wherein F is a stratum factor, R ₀ Is the resistivity of a rock sample containing 100% water, omega.m; r is R _w Is the resistivity of the water contained in the rock sample, Ω·m; phi is the porosity;

then, fitting the rock electrical measurement data of the rock sample by using an Archie formula to obtain a lithology coefficient b and a saturation index n:

wherein I is a resistivity increase coefficient; r is R ₀ Is the resistivity of a rock sample containing 100% water, omega.m; r is R _t Is the resistivity of the rock sample part when water is contained, and is omega.m; s is S _w Is the water saturation of the rock;

and then obtain:

since a and b are lithology coefficients, let ab=x, get:

the quantitative relationship between resistivity and water is obtained by fitting x, m, n using the test data.

6. The mining area water burst quantitative prediction method based on the ground transient electromagnetic according to claim 1, wherein the seismic exploration data of the working surface area is subjected to wave impedance inversion by adopting Geoview three-dimensional inversion software so as to construct a geophysical model; which comprises the following steps:

carrying out wave impedance inversion by adopting Geoview three-dimensional inversion software, and constructing a geological model below the region;

the same attribute units in the geologic model are the same geologic lithology, and the same electrical property is given to the same geologic lithology by combining resistivity logging data, so that the geophysical model is constructed.

7. The mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 1, wherein the geophysical model is used as an initial model, and inversion simulation is performed on corrected transient electromagnetic data to obtain underground resistivity distribution, and the calculation formula is as follows:

wherein d _c The theoretical value calculated for the current model is characterized in that the subscript i is the ith data, M represents the total number of observed data, M represents model parameters, the following table j represents the jth model parameter, and N ₁ Representing the total number of model parameters, P represents the model parameters after the same attribute is combined, N ₂ Representing the total number of model parameters after the same attribute is combined,is a partial derivative operator.

8. The mining area water inrush quantitative prediction method based on ground transient electromagnetic according to claim 7, wherein in the process of performing inversion simulation on corrected transient electromagnetic data by using a geophysical model as an initial model, the same geological lithology unit is used as the same inversion unit.

9. A mining area water inrush quantitative prediction system based on ground transient electromagnetic, comprising:

the acquisition module is used for acquiring the core of the working surface area of the mining area;

the testing module is used for measuring the resistivity of the core after the water-containing interval is air-dried and the resistivity of the core after the air-dried is injected with different water quantities by adopting a three-dimensional high-density resistivity method, so as to obtain the resistivity of the core under different water contents;

the rock core resistivity and water content relation building module is used for building a qualitative relation between the resistivity of the rock core and the water content according to the resistivity of the rock core under different water contents based on a classical Archie formula;

the analysis module is used for establishing an actual complex terrain model, analyzing and calculating the actual complex terrain model by adopting a three-dimensional vector finite element transient electromagnetic simulation algorithm to obtain the ground transient electromagnetic field responses of the electric source and the loop source under various fluctuating terrain conditions, so as to obtain the influence rule of the complex terrain on the ground transient electromagnetic response;

the geophysical model construction module is used for carrying out wave impedance inversion by adopting Geoview three-dimensional inversion software according to the seismic exploration data of the working surface area so as to construct a geophysical model;

the correction module is used for carrying out terrain correction on transient electromagnetic data of each measuring point in the working surface area according to the influence rule of complex terrain on the transient electromagnetic response of the ground;

the resistivity distribution acquisition module is used for carrying out inversion simulation on the transient electromagnetic data corrected by the terrain by taking the geophysical model as an initial model to obtain underground resistivity distribution;

and the prediction module is used for calculating the water content of the aquifer of the mining area to be detected by adopting inversion to obtain underground resistivity distribution according to a qualitative relation between the resistivity and the water content of the rock core, so as to quantitatively predict the water burst of the aquifer of the mining area.