CN117371267A - Mining area water inrush quantitative prediction method based on ground transient electromagnetic - Google Patents
Mining area water inrush quantitative prediction method based on ground transient electromagnetic Download PDFInfo
- Publication number
- CN117371267A CN117371267A CN202311209718.7A CN202311209718A CN117371267A CN 117371267 A CN117371267 A CN 117371267A CN 202311209718 A CN202311209718 A CN 202311209718A CN 117371267 A CN117371267 A CN 117371267A
- Authority
- CN
- China
- Prior art keywords
- water
- resistivity
- transient electromagnetic
- core
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 198
- 230000001052 transient effect Effects 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005065 mining Methods 0.000 title claims abstract description 37
- 239000011435 rock Substances 0.000 claims abstract description 78
- 230000004044 response Effects 0.000 claims abstract description 23
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims description 18
- 238000004088 simulation Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- 238000007605 air drying Methods 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 230000005672 electromagnetic field Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000004445 quantitative analysis Methods 0.000 abstract 1
- 239000003245 coal Substances 0.000 description 16
- 238000005516 engineering process Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 8
- 238000005553 drilling Methods 0.000 description 8
- 238000011156 evaluation Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000011835 investigation Methods 0.000 description 6
- 238000011158 quantitative evaluation Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005086 pumping Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a mining area water inrush quantitative prediction method based on ground transient electromagnetic, which relates to the technical field of mining area water inrush quantitative prediction and comprises the following steps: collecting rock in a working surface area of a mining area, measuring the resistivity of the rock after the water-bearing interval is air-dried and the resistivity of the rock after the air-dried rock is injected with different water quantities by adopting a three-dimensional high-density resistivity method, further obtaining the resistivity of the rock under different water quantities, and establishing a quantitative relation between the resistivity of the rock and the water content based on a classical Arch formula; analyzing and calculating the established actual complex terrain model by adopting a three-dimensional vector finite element transient electromagnetic simulation algorithm to obtain the influence rule of complex terrain on the ground transient electromagnetic response; further, a mine water inrush quantitative prediction model based on ground transient electromagnetic is established to quantitatively predict water inrush of a water-bearing layer of a mining area; the method establishes a qualitative relation between the rock resistivity and the water content, and carries out accurate quantitative analysis and prediction on the water content of the rock stratum.
Description
Technical Field
The invention relates to the technical field of quantitative prediction of water inrush in mining areas, in particular to a mining area water inrush quantitative prediction method based on ground transient electromagnetic.
Background
The water hazard is a disaster next to gas explosion in serious and oversized accidents of the coal mine. The water damage mine disaster not only easily causes the major casualties of underground operators, but also is most prominent in the aspects of economic loss severity, accident rescue difficulty, time required for recovering mine production and the like. With the improvement and development of the water control technology level in recent years, the occurrence frequency of water burst accidents of coal mines is in a continuous trend under the double strong supervision of countries and places, but the water burst accidents are not effectively restrained at all.
The geophysical hydrogeological detection technology used as a 'head soldier' has a certain application effect in preventing and controlling water in coal mines, but still cannot meet the technical requirements of preventing and controlling water in coal mines, and the theory, method, technology, application and other aspects of the technology need further innovative development so as to realize accurate detection of water inrush dangerous sources, effectively solve the problem of water inrush disaster, improve the technical level of preventing and controlling water in coal mines, reduce the economic cost investment of preventing and controlling water in coal mines and improve the safety and economic benefit of coal mines.
At present, the electric method and the electromagnetic geophysical method mainly adopted in the hydrogeology investigation of the coal field are both used for deducing and explaining water-rich abnormality according to the relative size of the underground resistivity value, are both indirect detection of water-imparting property, and cannot be used for quantitatively evaluating the water content of the rock stratum, and the root cause is as follows: (1) the qualitative relation between the rock resistivity and the water content is not established; (2) the measuring signal is influenced by various factors such as topography, instrument stability, human interference and the like, and has larger deviation with the real response of underground rock; (3) polynomials of geophysical inversion.
Disclosure of Invention
Aiming at the defect that the prior art is indirect detection of water imparting capacity and cannot be used for quantitative evaluation of water content of a rock stratum, the invention provides a mining area water inrush quantitative prediction method based on ground transient electromagnetism, and the mining area water inrush quantitative prediction theory and method based on ground transient electromagnetism are established by acquiring the relation between rock resistivity and water content and the response influence rule of complex terrain conditions on ground transient electromagnetic signals, so that the quantitative evaluation of water filling capacity of a water-bearing stratum is realized, and the problem that the prior art cannot be used for quantitative evaluation of water content of the rock stratum is solved.
A mining area water inrush quantitative prediction method based on ground transient electromagnetic comprises the following steps:
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.
Further, the method further comprises the steps of measuring the resistivity of the water-containing layer sections of the core after air drying by adopting a three-dimensional high-density resistivity method and the resistivity of the core after air drying when different water amounts are injected into the core, testing the porosities of the cores of all the water-containing layer sections by adopting a core nuclear magnetic resonance imaging analyzer, and calculating the water content of the water-containing layer in the working surface area according to the porosities and the qualitative relation between the resistivity of the core and the water content.
Further, the resistivity of the core with different water contents is obtained by sequentially increasing the water injection quantity 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.
Further, a model of the actual complex terrain is built by using DEM elevation data in combination with COMSOL software.
Further, based on the classical Archie formula, according to the resistivity of the core under different water contents, a qualitative relation between the resistivity of the core and the water content is established, and the specific steps comprise:
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 0 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 0 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.
Further, performing wave impedance inversion according to the seismic exploration data of the working surface area 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.
Further, the geophysical model is taken as an initial model, inversion simulation is carried out on the corrected transient electromagnetic data, and the underground resistivity distribution is obtained, wherein 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 1 Representing the total number of model parameters, P represents the model parameters after the same attribute is combined, N 2 Representing the total number of model parameters after the same attribute is combined,is a partial derivative operator.
Further, in the process of performing inversion simulation on the corrected transient electromagnetic data by taking the geophysical model as an initial model, the same geological lithology unit is taken as the same inversion unit.
Further, a mining area water inrush quantitative prediction system based on ground transient electromagnetic comprises:
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.
The invention provides a mining area water inrush quantitative prediction method based on ground transient electromagnetic, which comprises the following steps of
The beneficial effects are that:
according to the method, based on a classical Archie formula, according to the resistivity of cores with different water contents, the relation between the rock resistivity and the water content is obtained, a three-dimensional vector finite element transient electromagnetic simulation algorithm is adopted to analyze and calculate an actual complex terrain model, the influence rule of complex terrain conditions on the response of ground transient electromagnetic signals is obtained, terrain correction is carried out on high-precision transient electromagnetic data of each measuring point in a working surface area according to the influence rule of complex terrain on the ground transient electromagnetic response, inversion simulation is carried out on the terrain corrected transient electromagnetic data, and further, a mine water inrush quantitative prediction theory and method based on the ground transient electromagnetic are established, so that the quantitative evaluation error of water filling quantity of a water-bearing layer is smaller than 30%.
Drawings
FIG. 1 is a flow chart of a mining area water inrush quantitative prediction method based on ground transient electromagnetic;
FIG. 2 is a flow chart of the calculation of the quantitative relationship between rock resistivity and water content in an embodiment of the invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments.
As shown in figure 1, the invention surrounds the goal of quantitative evaluation of water filling amount of an aquifer, adopts the means of rock physics, numerical simulation, field test and the like, and realizes the quantitative evaluation error of water filling amount of the aquifer less than 30% by researching the relation between rock resistivity and water content, the influence rule of complex terrain conditions on the response of ground transient electromagnetic signals and the high-precision joint inversion method of transient electromagnetic under constraint conditions, thereby establishing the mine water inrush quantitative prediction theory and method based on ground transient electromagnetic.
(1) The hydrogeological prospecting technology for coal mines comprises the following steps: the traditional hydrogeological exploration means of coal deposits are single, and the logging and water pumping tests are mainly used. The water filling conditions are researched from multiple aspects, multiple angles and multiple layers by comprehensively surveying by means of remote sensing hydrogeological survey, geophysical exploration, hydrogeological drilling, water pumping test, underground water flow field analysis by water chemistry and environmental isotope analysis, physical and mechanical analysis of rock and soil mass and the like. Firstly, remote sensing hydrogeology investigation and ground geophysical prospecting are carried out, the conditions of the underground (surface) water system for the path supplementing and arranging and the positions of the well fields in the system are known, and the water enrichment of the water-filled aquifer is primarily judged. And then carrying out series of works such as drilling and pumping tests. And drawing a water-resisting layer structure contrast diagram, a water-bearing layer thickness contour diagram and a brittle rock and plastic rock ratio contour diagram according to drilling and coring results and by combining coal field geological exploration drilling data, and clearly showing the lithology and lithology structures of the well field through the three diagrams. The larger the thickness of the aquifer and the larger the ratio of the brittle plastic rock, the stronger the water enrichment of the aquifer. The traditional water sample test and water level observation are only used for establishing a water chemistry background value and knowing the groundwater level, and are not used for carrying out analysis and research on water filling conditions. At present, on the basis of sound basic data, the information of a water chemistry field and an underground water flow field is fully analyzed and researched, and a series of water chemistry diagrams and flow field diagrams are drawn to analyze the condition of the underground water repair and drainage.
(2) Water burst evaluation: for the prediction of the water inflow of the mine, an analytic method and a comparison method are adopted in the past, the prediction result of the mine is relatively reliable for the coal seam with simple conditions or stone-carbon two-stack system, but the deviation of the calculation result is larger for the mine with complex conditions. At present, a three-dimensional numerical simulation method of underground water is mainly adopted to predict the water inflow and drainage of a mine under the exploitation condition. The three-dimensional visual model of the hydrogeologic concept is characterized by Modflow software, so that the hydrogeologic condition of the mining area can be intuitively and deeply known from a three-dimensional angle. The method can better solve the complex problems of heterogeneous anisotropy and the like. The development of various investigation means obtains a great deal of geographical information, and provides important basis for the evaluation and prediction of the water-rich nature and the water inrush danger of the water-filled aquifer. However, if the analysis of water damage is directly and sporadically applied, the information is locally scattered, and the whole effect and the value contained in the investigation data cannot be fully exerted. In recent years, an aquifer water enrichment evaluation method and a mine water inrush risk prediction method based on a multi-source geochemical information composite superposition technology are developed, a three-graph-double prediction method is adopted for roof water damage, and a vulnerability index method is adopted for floor water damage on the basis of the water inrush coefficient evaluation method. The principle of the advanced evaluation technology is that a mathematical model is utilized to systematically analyze, count, simulate and superimpose a large amount of geochemical information obtained by various exploration means, so that an integral evaluation result which accords with actual conditions is obtained. Summarizing the development rule of the water damage evaluation prediction technology: from qualitative to quantitative, from single factor to multifactor, from homogeneous, planar to heterogeneous three-dimensional visualization.
(3) Hydrogeological geophysical exploration techniques: the hydrogeological exploration firstly utilizes a geophysical prospecting means to carry out omnibearing detection, then adopts a drilling means to verify the geophysical prospecting indirect detection result, and finally carries out drainage treatment. Geophysical prospecting is a prior and macroscopic control means and plays a very important role. Through basic theoretical research, field test and analysis comparison for many years, a set of more mature and advanced detection technology is formed. In terms of space and time, a ground-well-roadway combined fine exploration technology is formed, an abnormal part is firstly defined by a ground geophysical prospecting, then underground geophysical prospecting is carried out to be close to a target for fine exploration, and finally, abnormal is verified by utilizing drilling and data. Transient electromagnetic has become the main technical means of hydrogeological investigation in coal field mining areas at present.
In summary, the current electric method and electromagnetic geophysical method mainly adopted in the hydrogeology investigation of the coal field are both used for deducing and explaining water-rich abnormality according to the relative magnitude of the underground resistivity value, are both indirect detection of water-imparting property, and cannot be used for quantitatively evaluating the water content of the rock stratum: (1) the qualitative relation between the rock resistivity and the water content is not established; (2) the measuring signal is influenced by various factors such as topography, instrument stability, human interference and the like, and has larger deviation with the real response of underground rock; (3) polynomials of geophysical inversion. Therefore, implementing geophysical-based water inrush quantitative prediction techniques must solve the 3 problems described above. For this purpose, the present invention is intended to carry out the following studies:
(1) Qualitative relationship between rock resistivity and rock water content: based on the Archie formula, the change of the rock resistivity value under different water content conditions is studied.
(2) Influence rule of complex terrain on transient electromagnetic signals: aiming at the current coal field transient electromagnetic hydrogeological exploration technical means, the transient electromagnetic signal response characteristics of the return line source and the electrical source under the complex qualitative condition are researched, and a foundation is provided for obtaining the transient electromagnetic response of the real rock stratum.
(3) The high-precision constraint inversion method of the ground transient electromagnetic comprises the following steps: by combining drilling data and seismic data, the ground transient electromagnetic multi-parameter high-precision constraint inversion method is developed, and the accurate interpretation of the longitudinal positioning, thickness and real resistivity of the electrical layer is solved.
The method specifically comprises the following steps:
(1) A typical coal mine is selected as a test mine area, and drilling and coring are performed on all stratum in a representative working surface area, and a typical core is provided for petrophysical experiments.
(2) Carrying out petrophysical test indoors, firstly adopting a three-dimensional high-density rock core resistivity measuring system to test the resistivity of all the obtained rock cores; then, testing the porosity and the water content of the rock cores of all water-bearing intervals by adopting a MesoMR rock core nuclear magnetic resonance imaging analyzer; and then air-drying the water-bearing stratum rock core to remove all water, and testing the resistivity of the air-dried water-free rock by adopting a three-dimensional high-density rock core resistivity measuring system again to establish a rock resistivity database of each stratum.
(3) The method comprises the steps of sequentially increasing water injection quantity and carrying out resistivity measurement on air-dried rock in a water injection displacement mode to obtain rock core resistivity under different water contents, and establishing a quantitative relation between the rock resistivity and the water content by using a polynomial fitting technology based on a classical Arch formula, wherein the quantitative relation is shown in fig. 2 and specifically comprises the following steps:
firstly, fitting rock electricity experimental data of a rock sample by using an Alqi formula (1) to obtain a and m
Wherein R is 0 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; a is lithology coefficient; m is the cementation index.
Then using Alqi formula (2) to fit the rock electricity experimental data of the rock sample to obtain b and n
Wherein I is a resistivity increase coefficient; r is R 0 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; b is lithology coefficient; n is the saturation index.
Multiplying equation (1) by equation (2):
since a and b are lithology coefficients, let ab=x, get:
in actual operation, experimental data can be used to fit x, m, n, thereby obtaining the quantitative relationship between resistivity and water.
(4) The method comprises the steps of using DEM elevation data and COMSOL software to establish an actual complex terrain model, calculating the ground transient electromagnetic field response of an electric source and a loop source under various fluctuating terrain conditions by adopting a three-dimensional vector finite element transient electromagnetic simulation algorithm, summarizing the influence rule of complex terrain on the ground transient electromagnetic response, and providing a foundation for terrain influence correction and underground stratum actual transient electromagnetic signal acquisition.
(5) And (3) taking the existing seismic data and logging data of the test area, carrying out wave impedance inversion on the processed seismic exploration data by using mature software Geoview three-dimensional inversion software, constructing a geological model below the test area, wherein the same attribute units in the inversion model are the same geological lithology, and endowing the same geological lithology with the same electrical property by combining the resistivity logging data, so as to construct the three-dimensional ground model.
(6) And (3) carrying out high-precision transient electromagnetic data acquisition in a measuring area, recording the elevation of each measuring point, and carrying out terrain correction on the data of each measuring point by combining the research result in the step (4).
(7) And (3) forward modeling is carried out on the transient electromagnetic data acquired in the step (6) and subjected to terrain correction by taking the geophysical model constructed in the step (5) as an initial model, and in the modeling process, the same geological lithology units are taken as the same inversion units, so that the number of inversion model units (as follows) is reduced, the inversion efficiency is improved, the electrical units below a measurement area are kept to have higher resolution, the influence of electromagnetic data volume effect is controlled to be reduced by an electrical boundary, and more accurate underground resistivity distribution is obtained.
(8) And (3) combining the relation among the porosity, the water content and the resistivity of the underground rock obtained in the step (2), adopting the resistivity obtained by inversion in the step (7), calculating the water content of the water-containing layer in the range of the area on the basis of delineating the water-rich area of the water-containing layer, and realizing the spatial distribution visualization of the water content in the range of the area according to a geological model constructed according to the seismic inversion result.
(9) A representative aquifer was selected, drilled to drain and the water inflow was recorded.
(10) And comparing the actual water inflow with the predicted water inflow, and verifying the prediction error.
(11) And establishing a mine water inrush quantitative prediction theory and technology based on ground transient electromagnetic, and realizing that the quantitative prediction error of the water quantity of the water-filled aquifer is less than 30%.
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art, who is within the scope of the present invention, should make equivalent substitutions or modifications according to the technical scheme of the present invention and the inventive concept thereof, and should be covered by the 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 0 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 0 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 1 Representing the total number of model parameters, P represents the model parameters after the same attribute is combined, N 2 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.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311209718.7A CN117371267B (en) | 2023-09-19 | 2023-09-19 | Mining area water inrush quantitative prediction method based on ground transient electromagnetic |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311209718.7A CN117371267B (en) | 2023-09-19 | 2023-09-19 | Mining area water inrush quantitative prediction method based on ground transient electromagnetic |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117371267A true CN117371267A (en) | 2024-01-09 |
| CN117371267B CN117371267B (en) | 2024-05-14 |
Family
ID=89390103
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311209718.7A Active CN117371267B (en) | 2023-09-19 | 2023-09-19 | Mining area water inrush quantitative prediction method based on ground transient electromagnetic |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117371267B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118392934A (en) * | 2024-04-25 | 2024-07-26 | 山东大学 | Quality monitoring method and system for mining reclamation soil |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105223612A (en) * | 2015-06-10 | 2016-01-06 | 中国矿业大学 | A kind of coal mine flood prediction and evaluation method based on earthquake information |
| CN111255446A (en) * | 2020-01-19 | 2020-06-09 | 长江大学 | Resistivity correction method based on stratum simulation |
| CN114839671A (en) * | 2022-03-28 | 2022-08-02 | 中国矿业大学 | Method for finely identifying coal measure stratum electrical interface by ground and ground roadway transient electromagnetic method |
| US20230031116A1 (en) * | 2022-09-27 | 2023-02-02 | Chengdu University Of Technology | Method of quantitative evaluation on structural disturbance characteristics of present in-situ geo-stress in deep shale gas reservoirs |
-
2023
- 2023-09-19 CN CN202311209718.7A patent/CN117371267B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105223612A (en) * | 2015-06-10 | 2016-01-06 | 中国矿业大学 | A kind of coal mine flood prediction and evaluation method based on earthquake information |
| CN111255446A (en) * | 2020-01-19 | 2020-06-09 | 长江大学 | Resistivity correction method based on stratum simulation |
| CN114839671A (en) * | 2022-03-28 | 2022-08-02 | 中国矿业大学 | Method for finely identifying coal measure stratum electrical interface by ground and ground roadway transient electromagnetic method |
| US20230031116A1 (en) * | 2022-09-27 | 2023-02-02 | Chengdu University Of Technology | Method of quantitative evaluation on structural disturbance characteristics of present in-situ geo-stress in deep shale gas reservoirs |
Non-Patent Citations (4)
| Title |
|---|
| 刘小平;刘森;钟秋;: "瞬变电磁法探测水文地质因素在复杂地表条件下的应用", 地下水, no. 01, 25 January 2014 (2014-01-25) * |
| 姜国庆;程久龙;王玉和;贾春梅;: "矿井瞬变电磁法探测工作面顶板岩层含水性", 煤矿安全, no. 10, 10 October 2010 (2010-10-10) * |
| 张军;: "矿井瞬变电磁法在顶板富水评价中的应用研究", 煤炭技术, no. 04, 10 April 2020 (2020-04-10) * |
| 田伟;邹冠贵;汤小明;曾葫;: "基于AVO技术与伪泊松比属性预测灰岩富水性的方法", 煤炭学报, no. 10, 15 October 2017 (2017-10-15) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118392934A (en) * | 2024-04-25 | 2024-07-26 | 山东大学 | Quality monitoring method and system for mining reclamation soil |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117371267B (en) | 2024-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10386531B2 (en) | Geological model analysis incorporating cross-well electromagnetic measurements | |
| CN104330838B (en) | Method applied to fine detection of sudden water burst channel in underground engineering and evaluation of grouting and plugging effects of sudden water burst channel | |
| Schmidt-Hattenberger et al. | Application of a Vertical Electrical Resistivity Array (VERA) for monitoring CO2 migration at the Ketzin site: First performance evaluation | |
| CN104331745B (en) | In oil-gas reservoir intrinsic fracture by stages, be divided into because of prediction and evaluation method | |
| CN117371267B (en) | Mining area water inrush quantitative prediction method based on ground transient electromagnetic | |
| Zappone et al. | Fault sealing and caprock integrity for CO 2 storage: an in-situ injection experiment | |
| CN101952744A (en) | Be used to explain the workflow based on model of deep investigation electromagnetic data | |
| CN104636980B (en) | Collect the geophysics characterizing method of condition for channel reservoir type oil gas | |
| US11789170B2 (en) | Induced seismicity | |
| CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
| CN105956271B (en) | A kind of tunneling shield construction causes the calculation method of formation displacement | |
| CN214366030U (en) | Horizontal well micro-seismic monitoring system based on distributed optical fiber sensing | |
| CN111352172A (en) | Method for acquiring spatial distribution position of uranium anomaly in sand body by well-seismic combination method | |
| Su et al. | Small fixed‐loop transient electromagnetic in tunnel forward geological prediction | |
| Al Hagrey | Electric study of fracture anisotropy at Falkenberg, Germany | |
| CN113484910A (en) | Tunnel advanced geological prediction method and system based on seismic interference method | |
| Jenkins et al. | Fit for purpose monitoring-a progress report on the CO2CRC Otway Stage 3 project | |
| CN113433595B (en) | Advanced prediction method based on natural electric field tunnel crack water | |
| Wang et al. | A new attempt for the detection of potential water-bearing structures in tunnels via hydraulic tomography: The first numerical investigation | |
| CN114839671A (en) | Method for finely identifying coal measure stratum electrical interface by ground and ground roadway transient electromagnetic method | |
| Wagner | New developments in electrical resistivity imaging with applications to geological CO₂ storage | |
| Li et al. | Experimental study on comprehensive geophysical advanced prediction of water-bearing structure of shaft | |
| Gale et al. | Progress in the hydrogeological characterization of the Stripa site | |
| Nakajima et al. | Geophysical monitoring at the Nagaoka pilot-scale CO2 injection site in Japan | |
| Land | Utilizing a digital model to determine the hydraulic properties of a layered aquifer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |