CN118131340A - Mine geophysical prospecting data analysis method based on electrical prospecting technology - Google Patents
Mine geophysical prospecting data analysis method based on electrical prospecting technology Download PDFInfo
- Publication number
- CN118131340A CN118131340A CN202410552301.9A CN202410552301A CN118131340A CN 118131340 A CN118131340 A CN 118131340A CN 202410552301 A CN202410552301 A CN 202410552301A CN 118131340 A CN118131340 A CN 118131340A
- Authority
- CN
- China
- Prior art keywords
- goaf
- region
- mine
- extremum
- resistivity
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000005516 engineering process Methods 0.000 title claims abstract description 28
- 238000007405 data analysis Methods 0.000 title claims abstract description 11
- 239000011435 rock Substances 0.000 claims description 29
- 230000002159 abnormal effect Effects 0.000 claims description 18
- 238000005065 mining Methods 0.000 claims description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 10
- 239000011707 mineral Substances 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004891 communication Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 3
- 239000000945 filler Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the field of electromagnetic exploration, in particular to a mine geophysical prospecting data analysis method based on an electrical prospecting technology, which comprises the following steps: acquiring mine geophysical prospecting data based on an electrical prospecting technology, and analyzing according to the geophysical prospecting data to obtain resistance data corresponding to each area of a mine; analyzing the resistance data, and extracting an extremum region; calculating the resistivity change difference of the extremum region, and determining the region property of the extremum region according to the resistivity change difference; if the regional property of the extremum region is a goaf, determining the position of the goaf according to geophysical prospecting data; and carrying out secondary survey on the goaf based on the electric exploration technology according to the position of the goaf, and determining the internal condition of the goaf. The invention improves the efficiency of locating and analyzing the goaf; the method solves the technical problems that the prior art needs a great deal of time and energy for positioning and analyzing the measurement result by technicians.
Description
Technical Field
The invention relates to the field of electromagnetic exploration, in particular to a mine geophysical prospecting data analysis method based on an electrical prospecting technology.
Background
The electrical prospecting is a geophysical prospecting method for prospecting geological structures and searching useful mineral products by observing and researching the space distribution rule and time characteristics of an electromagnetic field or an electrochemical field according to the difference of electromagnetic conductivity, magnetic permeability, dielectric property and electrochemical characteristics of various rocks or minerals in the crust; the electric prospecting technology is divided into a direct current method, namely a resistance method, a pulse transient field method and the like according to the time characteristics of an electromagnetic field; for exploring differences in electrical properties of subterranean formations.
In the exploitation and development process of metal mineral resources, the traditional goaf is mostly not effectively treated due to various reasons such as exploitation technology, geological environment and the like, and serious influences and even serious losses are caused on subsequent mine construction, safe exploitation, reserve calculation and the like. Therefore, the spatial position and occurrence condition of the goaf are determined to be important for comprehensive treatment and development of the mine.
Drawing a corresponding resistivity image according to different electrical characteristics of the underground rock stratum environment in a two-dimensional resistivity model image obtained based on an electrical exploration technology; goafs of different regions appear as high-resistance anomaly regions or low-resistance anomaly regions in the resistivity image; in the prior art, a professional related technician is required to analyze and identify a corresponding abnormal region in the abnormal resistivity profile, the abnormal region is easy to be interfered by typical properties of surrounding rock strata, and the manual analysis is easy to cause errors or misjudgment on the judgment of the region.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a mine geophysical prospecting data analysis method based on an electrical prospecting technology, and the adopted technical scheme is as follows:
a mine geophysical prospecting data analysis method based on an electrical prospecting technology comprises the following steps:
acquiring mine geophysical prospecting data based on an electrical prospecting technology, and analyzing according to the geophysical prospecting data to obtain resistance data corresponding to each area of a mine;
analyzing the resistance data, and extracting an extremum region;
calculating the resistivity change difference of the extremum region, and determining the region property of the extremum region according to the resistivity change difference;
If the regional property of the extremum region is a goaf, determining the position of the goaf according to the geophysical prospecting data;
and carrying out secondary survey on the goaf based on an electrical prospecting technology according to the position of the goaf, and determining the internal condition of the goaf.
Further, the acquiring mine geophysical prospecting data based on the electrical prospecting technology further comprises:
Measuring a set area above a mining area where mining exists by adopting a high-frequency electromagnetic sounding method based on an electrical prospecting technology;
and obtaining the mine geophysical prospecting data.
Further, the geophysical prospecting data at least includes: mine depth, mine width, and mine depth, and mine width.
Further, the analyzing the resistance data, and extracting the extremum region further includes:
Determining the total resistance range of the mine according to the resistance data, and selecting a preset extremum range as an extremum region; wherein the extremum region comprises: a maximum value region and a minimum value region;
And acquiring a connected domain formed by the corresponding extremum region, and constructing a corresponding connected domain set.
Further, the calculating the difference in resistivity change of the extremum region further includes:
according to the points in the connected domain set, points which are obtained from the edge of the outermost layer of the connected domain and the points which are opposite to the gradient change of the points are obtained outwards in the vertical direction of the edge of the outermost layer of the connected domain;
forming a plurality of line segments from points obtained from the edge of the outermost layer of the connected domain and points with opposite gradient changes;
and calculating the variation difference of the resistivity difference values of the unit distances on the line segments.
Further, the regional properties include at least: goaf, subsidence area and rock of the mineral seam of corresponding resistance.
Further, the determining the region property of the extremum region from the resistivity variation difference further comprises:
calculating the average resistivity change degree of the points inside the connected domain and the surrounding ore layer rock:
calculating the first possibility that the corresponding extremum region belongs to the mining area or the subsidence area according to the resistivity difference value variation difference and the average resistivity variation degree;
and if the first possibility is larger than a first preset threshold value, judging that the extremum area belongs to the goaf or the subsidence area.
Further, the determining the region property of the extremum region from the resistivity variation difference further comprises:
Defining a minimum circumscribed rectangle of the extremum region;
Acquiring the height and the width of the rectangle, and acquiring the mass center of the communication domain corresponding to the extremum region according to the difference of the resistivity change amplitude of the extremum region in the vertical upward direction;
calculating the maximum value of the resistivity change gradient according to the centroid of the connected domain;
Calculating a second possibility that the extremum area belongs to the goaf according to the height, the width and the maximum value of the resistivity change gradient of the rectangle;
And if the second possibility is larger than a second preset threshold value, determining the extremum area as a goaf.
Further, if the regional property of the extremum region is a goaf, determining the position of the goaf according to the geophysical prospecting data further includes:
And acquiring the corresponding connected domain of the goaf, obtaining the region range of the corresponding connected domain, and determining the position according to the depth and the width displayed in the resistivity profile.
Further, the secondary survey is performed on the goaf based on the electrical prospecting technology according to the position of the goaf, and determining the internal condition of the goaf further includes:
Measuring again by using the high-frequency electromagnetic sounding method with the position as a central point, and determining the specific position of the final goaf according to the two results;
when the goaf is expressed as a low-resistance abnormal region in the resistivity profile, the goaf is correspondingly filled with water, and the water contains more minerals;
when the corresponding goaf is represented as a high-resistance abnormal region in the resistivity profile, the interior of the corresponding goaf is preserved intact and is not filled with water.
The invention has the following beneficial effects: acquiring mine geophysical prospecting data based on an electrical prospecting technology, and analyzing according to the geophysical prospecting data to obtain resistance data corresponding to each area of a mine; analyzing the resistance data, and extracting an extremum region; calculating the resistivity change difference of the extremum region, and determining the region property of the extremum region according to the resistivity change difference; if the regional property of the extremum region is a goaf, determining the position of the goaf according to geophysical prospecting data; and carrying out secondary survey on the goaf based on the electric exploration technology according to the position of the goaf, and determining the internal condition of the goaf. According to the invention, the corresponding abnormal resistance area is obtained by carrying out image recognition on the related section two-dimensional resistivity model diagram obtained by electric exploration, the goaf position information and the internal structure information are obtained by combining the related goaf abnormal conditions, and the goaf positioning analysis efficiency is improved by carrying out positioning and component analysis on the goaf; the method solves the technical problems that the prior art needs a great deal of time and energy for positioning and analyzing the measurement result by technicians.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions and advantages of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for analyzing geophysical prospecting data in a mine based on an electrical prospecting technique according to an embodiment of the present invention;
FIG. 2 is a cross-sectional view of a two-dimensional resistivity model image in accordance with one embodiment of the invention;
FIG. 3 is a cross-sectional view of a two-dimensional resistivity model image, in accordance with one embodiment of the present invention;
FIG. 4 is a cross-sectional view of a two-dimensional resistivity model image in accordance with one embodiment of the invention;
FIG. 5 is a cross-sectional view of a two-dimensional resistivity model image, in accordance with one embodiment of the invention.
Detailed Description
In order to further describe the technical means and effects adopted by the invention to achieve the preset aim, the following is a detailed description of specific implementation, structure, characteristics and effects of the method for analyzing the geophysical prospecting data of the mine based on the electrical prospecting technology according to the invention, with reference to the accompanying drawings and the preferred embodiment. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The invention provides a specific scheme of a mine geophysical prospecting data analysis method based on an electrical prospecting technology, which is specifically described below with reference to the accompanying drawings.
In the prior art, mine geophysical prospecting data are acquired through an electrical prospecting technology, for example, a two-dimensional resistivity model image of a drawing section of bottom layer information of a relevant area is acquired through a high-frequency electromagnetic sounding method, a goaf can be presented in a high-resistance or low-resistance mode in the two-dimensional resistivity model image, and as a result, relevant technicians are required to analyze and position.
Specifically, a high-frequency electromagnetic sounding method is adopted to measure the upper part of a mined mining area, and the mining area is distributed in a multilayer structure due to complex geological environment, and the direction and the extension of the deployment are different; collapse occurs in some cases, water filling or other fillers exist in some cases, and the goaf resistance properties are greatly different due to different fillers; for example, when the goaf is left intact and substantially free of water, the goaf exhibits a typical high resistance compared to the surrounding rock; the electrical properties of goaf when it collapses or has filler are related to filler properties, and in the case of blocky rock it is highly resistive and in the case of crushed stone or more water it is a relatively low resistive body.
Referring to fig. 1, a flowchart of a method for analyzing geophysical prospecting data in a mine according to an embodiment of the present invention is shown, where the method includes:
and step S110, acquiring the geophysical prospecting data of the mine based on an electrical prospecting technology, and analyzing according to the geophysical prospecting data to obtain the resistance data corresponding to each area of the mine.
In an alternative embodiment, step S110 further includes: measuring a set area above a mining area where mining exists by adopting a high-frequency electromagnetic sounding method based on an electrical prospecting technology; and obtaining the mine geophysical prospecting data.
In an alternative embodiment, the geophysical prospecting data at least include: mine depth, mine width, and mine depth, and mine width.
In the step, a high-frequency electromagnetic sounding method is used for collecting two-dimensional resistivity model images of a profile drawn by bottom layer information of a relevant area, the high-frequency electromagnetic sounding method is used for measuring above a mining area where mining exists, the mining area is combined with local geological data and relevant known drilling data, a torsion fracture structure develops in the mining area, goafs are distributed in a multi-layer structure, and a resistivity abnormal area can be shown in a relevant resistivity profile obtained through detection. FIGS. 2-5 are cross-sectional views, one-four, of two-dimensional resistivity model images, as shown in FIGS. 2-5, with the horizontal axis being the mine width (distance/m) and the vertical axis being the mine depth (depth/m), with the shaded portions of the cross-sectional views showing the resistance (in ohm) corresponding to the depth and width.
And S120, analyzing the resistance data and extracting an extremum region.
In the step, a resistance abnormal region is positioned through a measured related resistivity profile, specifically, for the distribution of rock stratum of an underground mining area, the change among the resistivity contours in an unexplored region is small, the distribution is in continuous layered distribution, an underground goaf presents high resistance or low resistance abnormality in an image, the contours present irregular bulk closure conditions, the resistance change of the contour regions of a closure region and a peripheral region is large, and the distribution is irregular. Therefore, in order to locate the relevant goaf position information, the connected domain analysis is performed on the region composed of the high resistance and the low resistance, and the distribution and the resistance change of the surrounding region are judged. After the abnormal resistance area is determined, distinguishing the goaf and the collapse area according to the area condition of the connected area; the goaf and the surrounding resistance have larger change trend, the goaf area is smaller, the resistivity change degree of the collapse area in the vertical upward direction is smaller, and the resistivity change condition in other directions is larger.
In an alternative embodiment, step S120 further includes: determining a resistance range according to the resistance data, and selecting a preset extremum range as an extremum region; wherein the extremum region comprises: a maximum value region and a minimum value region; and acquiring a connected domain formed by the corresponding extremum region, and constructing a corresponding connected domain set.
Specifically, according to the obtained resistivity profile, the region where the resistance extremum is located is obtained, and according to the corresponding resistance range, the extremum range to be detected is selected, for example, in the profile, the resistance change range is between 0 and 2000, the extremum region is defined according to the two sides of the interval, for example, the preset extremum range is set to ten percent, and the region with ten percent on the two sides is defined as the extremum region; its minimum value region isThe maximum area is/>When the resistance change range is between 0 and 2000, 10% of the two sides is the minimum value area, namely 10% of 2000 is 0 to 200; the maximum value area takes 10% of 2000, namely 1800-2000; wherein, 0-200 and 1800-2000 are resistance values, the unit is omega, the connected domain formed by the corresponding extremum region is obtained, and the connected domain set formed by the extremum region is i=/>. Wherein,，/>Assigning a value to the resistor; /(I)Letters are assigned to connected domains.
And S130, calculating the resistivity change difference of the extremum region, and determining the region property of the extremum region according to the resistivity change difference.
In an alternative embodiment, step S130 further includes: according to the points in the connected domain set, the points which are obtained from the edge of the outermost layer of the connected domain and the points which are opposite to the gradient change of the points are obtained outwards in the vertical direction of the edge of the outermost layer of the connected domain; forming a plurality of line segments from points obtained from the edge of the outermost layer of the connected domain and points with opposite gradient changes; the change in the difference in resistivity per unit distance over a plurality of line segments is calculated.
In an alternative embodiment, the zone properties include at least: goaf, subsidence area and rock of the mineral seam of corresponding resistance.
That is, the regional properties of the extremum regions selected in the acquired resistivity profile may be goaf, subsidence area, and the corresponding resistivity seam rock. Wherein, the maximum value area and the minimum value area of the rock area of the ore layer are formed by the electrical characteristics of the corresponding rock stratum, and the rock area is continuously layered distributed on the ground morphology distribution, and the degree of change of the resistivity of the rock area and the surrounding rock stratum is small; the goaf and the subsidence area are greatly different from the surrounding rock stratum, and the goaf and the subsidence area have larger resistivity mutation degree in the resistivity profile. Therefore, firstly, the rock of the mineral seam and the goaf formed by the extremum resistance area are distinguished according to the corresponding resistivity abrupt change degree, and the subsidence area is formed.
In an alternative embodiment, step S130 further includes: calculating the average resistivity change degree of the points inside the connected domain and the surrounding ore layer rock: calculating the first possibility that the corresponding extremum region belongs to the mining area or the subsidence area according to the difference of the resistivity values and the average resistivity change degree; if the first probability is greater than a first preset threshold value, judging that the extremum area belongs to the goaf or the subsidence area.
Specifically, the rock area of the ore layer is continuously layered in the ground form distribution, the resistivity change difference of the rock area per unit distance is small, the goaf or the subsidence area is more disordered in the ground form distribution, and the resistivity change difference of the rock area per unit distance is large; therefore, the probability that the connected domain belongs to the goaf or the subsidence area is determined according to the resistance change difference of the unit distance; specifically, all the points are obtained on the edge of the outermost layer of the connected domain, the number of the obtained points is n, the points with opposite gradient changes are obtained outwards along the direction perpendicular to the edge of the outermost layer of the connected domain from the points obtained on the edge of the outermost layer of the connected domain, a line segment with the length of f meters is formed by the points with opposite gradient changes, and the formed n line segments have the change of the difference value of the resistance of f unit distances, wherein the formula is shown in the following formula (1):
；（1）
Wherein, Representing the change of the difference value of the resistivity according to the unit distance on the ith line segment, for determining the possibility that the connected domain belongs to the goaf or the subsidence area; /(I)To average the change of the resistivity difference value of the unit distance on a plurality of line segments; Representing the resistivity difference of the kth unit distance over the distance of f meters acquired at the nth point on the connected domain i; /(I) Is the average value of resistivity differences per unit distance on an outwards derived line segment corresponding to the nth point on the connected domain i,/>Is the standard deviation of resistivity differences per unit distance on the outwardly derived line segment corresponding to the nth point on connected domain i.
In particular, the method comprises the steps of,The kurtosis of the resistivity difference per unit distance of the nth point on the connected domain i is represented. When the connected domain belongs to a goaf or a collapse area, the obtained unit distance resistivity differences at n points are large in difference, and corresponding kurtosis is increased. When the connected domain belongs to the rock region of the ore layer, the resistivity difference per unit distance of n points obtained by continuously and lamellar distribution is uniform, and the kurtosis is low. So when/>The larger the probability that the connected domain belongs to the goaf or the subsidence area is, when/>The smaller the probability that the connected domain belongs to the goaf or the subsidence area.
Further, the maximum difference in resistivity of the line segment obtained at the nth point inside the ith connected domain in the cross-sectional view isThe average resistivity variation degree between the connected domain i and the surrounding rock layer is/>Calculate the average degree of resistivity change as/>Formula (2) below:
；（2）
Wherein, The average resistivity variation degree between the connected domain i and the surrounding rock stratum f meters; /(I)Representing the maximum difference of the resistivity of the obtained line segment in the section view at the nth point on the connected domain i; n is the number of all the points on the edge of the outermost layer of the connected domain.
When (when)When the resistivity of the connected domain and the surrounding rock stratum is larger, the probability of belonging to the goaf or the subsidence area is larger; when/>The smaller the resistivity of the connected domain and the surrounding rock layer, the smaller the degree of change of the resistivity, and the smaller the probability of belonging to a goaf or a subsidence area.
According to the obtained result, a first possibility that the corresponding extremum region belongs to the mining area or the subsidence area can be obtainedThe calculation formula is as follows (3):
；（3）
Wherein, Representing the possibility that the connected domain belongs to the goaf or the subsidence area according to the resistance change difference of the unit distance; /(I)Determining the possibility that the connected domain belongs to a goaf or a subsidence area for the average resistivity change degree of the connected domain i and the surrounding rock stratum f meter; if the first probability is greater than a first preset threshold value, judging that the extremum area belongs to the goaf or the subsidence area; the first preset threshold is set to z, then when/>When z, judging that the extremum area belongs to the goaf and/or the collapse area; the corresponding connected domain set is assigned as/>Wherein/>Assigning letters to the letters; the first preset threshold is set by a person skilled in the art according to the actual scenario.
In an alternative embodiment, step S130 further includes: defining a minimum circumscribed rectangle of the extremum region; acquiring the height and the width of a rectangle, and acquiring the mass center of a communication domain corresponding to the extremum region according to the difference of the resistivity change amplitude of the extremum region in the vertical upward direction; calculating the maximum value of the resistivity change gradient according to the mass center of the connected domain; calculating the second possibility that the extremum area belongs to the goaf according to the height, width and maximum value of the resistivity change gradient of the rectangle; and if the second possibility is larger than a second preset threshold value, determining the extremum area as the goaf.
The goaf and the subsidence area are further distinguished, according to the characteristic distinction of the goaf and the subsidence area, the subsidence area has larger resistivity variation amplitude in the horizontal direction and the vertical downward direction, the trend of large area and downward bulge is shown on the ground form, the resistivity variation amplitude of the goaf in each direction is larger, and the goaf is shown as smaller bulk distribution form on the ground form.
Therefore, the minimum circumscribed rectangle of the extremum region is defined, the height is h, the width is w, according to the circumscribed rectangle, the characteristic extraction technology in image processing is utilized, according to the difference of the resistivity change amplitude of the corresponding extremum region in the vertical upward direction, the centroid of the corresponding connected region is obtained, the point of which the gradient direction is opposite to the change is derived and obtained along the earth surface direction to form a line segment, and the maximum value of the corresponding resistivity change gradient is obtained; the formula for calculating the second possibility that the extremum region belongs to the goaf is as follows (4):
；（4）
Wherein, Representing a second possibility that a communicating region a formed by the corresponding goaf and the subsidence area belongs to the goaf; /(I)The maximum value of the resistivity change gradient of a line segment formed by points, the gradient direction of which is derived from the centroid of the connected domain a along the earth surface direction and is opposite to the change of the centroid; /(I)Representing the height of the smallest circumscribed rectangle of the corresponding connected domain a; /(I)The width of the smallest circumscribed rectangle corresponding to the connected domain a is indicated.
In the formula (4) of the present invention,Representing a second possibility that the connected domain belongs to the goaf according to the feature distinction of the goaf and the subsidence area in the ground and the corresponding resistivity in the vertical direction, in particular, when the gradient change is larger and the aspect ratio of the minimum circumscribed rectangle thereof is more approximate to 1, the second possibility that the connected domain belongs to the goaf is larger; setting the second preset threshold value as/>Then when/>Indicating that the connected domain belongs to the goaf; wherein the second preset threshold is set by a person skilled in the art according to the actual scenario.
And step 140, if the regional property of the extremum region is the goaf, determining the position of the goaf according to the geophysical prospecting data.
In an alternative embodiment, step S140 further includes: and acquiring the corresponding connected domain of the goaf, obtaining the region range of the corresponding connected domain, and determining the position according to the depth and the width displayed in the resistivity profile.
Specifically, if the extreme value region corresponding to the secondary goaf is a maximum value region, judging that the goaf is a high-resistance abnormal region; and if the corresponding connected region of the secondary goaf is a minimum value region, judging that the goaf is a low-resistance abnormal region.
And step S150, performing secondary survey on the goaf based on the electric prospecting technology according to the position of the goaf, and determining the internal condition of the goaf.
Specifically, according to the obtained corresponding connected domain of the goaf, obtaining the region range of the corresponding connected domain, and determining the position according to the depth and the width displayed in the resistivity profile; and then measuring by using the position as a central point by adopting a high-frequency electromagnetic sounding method, and finally determining the specific position of the goaf according to the two results.
In an alternative embodiment, step S150 further includes: measuring again by using a high-frequency electromagnetic sounding method with the position as a central point, and determining the accurate position of the final goaf according to the two results; when the goaf judged at the time shows a low-resistance abnormal region in the resistivity profile, the goaf is filled with water, and the water contains more minerals; when the corresponding goaf is represented as a high-resistance abnormal region in the resistivity profile, the interior of the corresponding goaf is preserved intact and is not filled with water.
In particular, if the secondary survey result and the primary survey result are too different, the accurate position of the goaf cannot be determined, and the secondary survey result is discarded. And repeating the steps after the current two-dimensional resistivity model image is acquired again for calculation. By adopting the method of the embodiment, the corresponding abnormal resistance area is obtained by carrying out image recognition on the related section two-dimensional resistivity model obtained by electric exploration, the goaf position information and the internal structure information are obtained by combining the related goaf abnormal conditions, and the goaf positioning analysis efficiency is improved by carrying out positioning and component analysis on the goaf; the method solves the technical problems that the prior art needs a great deal of time and energy for positioning and analyzing the measurement result by technicians.
It should be noted that: the sequence of the embodiments of the present invention is only for description, and does not represent the advantages and disadvantages of the embodiments. The processes depicted in the accompanying drawings do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
In this specification, each embodiment is described in a progressive manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments.
Claims (10)
1. The mine geophysical prospecting data analysis method based on the electrical prospecting technology is characterized by comprising the following steps of:
acquiring mine geophysical prospecting data based on an electrical prospecting technology, and analyzing according to the geophysical prospecting data to obtain resistance data corresponding to each area of a mine;
analyzing the resistance data, and extracting an extremum region;
calculating the resistivity change difference of the extremum region, and determining the region property of the extremum region according to the resistivity change difference;
If the regional property of the extremum region is a goaf, determining the position of the goaf according to the geophysical prospecting data;
and carrying out secondary survey on the goaf based on an electrical prospecting technology according to the position of the goaf, and determining the internal condition of the goaf.
2. The electrical prospecting technique based mine geophysical prospecting data analysis method according to claim 1, wherein the collecting mine geophysical prospecting data based on the electrical prospecting technique further comprises:
Measuring a set area above a mining area where mining exists by adopting a high-frequency electromagnetic sounding method based on an electrical prospecting technology;
and obtaining the mine geophysical prospecting data.
3. The method for analyzing geophysical prospecting data for a mine based on the electrical prospecting technique according to claim 2, wherein the geophysical prospecting data at least include: mine depth, mine width, and mine depth, and mine width.
4. The method for analyzing the geophysical prospecting data for a mine based on the electrical prospecting technique according to claim 3, wherein analyzing the resistance data, extracting the extremum region further comprises:
Determining the total resistance range of the mine according to the resistance data, and selecting a preset extremum range as an extremum region; wherein the extremum region comprises: a maximum value region and a minimum value region;
And acquiring a connected domain formed by the corresponding extremum region, and constructing a corresponding connected domain set.
5. The method of analyzing geophysical prospecting data for a mine based on the electrical prospecting technique according to claim 4, wherein said calculating the difference in resistivity change in the extremum region further comprises:
according to the points in the connected domain set, points which are obtained from the edge of the outermost layer of the connected domain and the points which are opposite to the gradient change of the points are obtained outwards in the vertical direction of the edge of the outermost layer of the connected domain;
forming a plurality of line segments from points obtained from the edge of the outermost layer of the connected domain and points with opposite gradient changes;
and calculating the variation difference of the resistivity difference values of the unit distances on the line segments.
6. The method for analysis of geophysical prospecting for a mine based on electrical prospecting techniques according to any one of claims 1 to 5, wherein the properties of the region comprise at least: goaf, subsidence area and rock of the mineral seam of corresponding resistance.
7. The electrical prospecting techniques based mine geophysical data analysis method according to claim 6, wherein determining the regional nature of the extremum region from the resistivity variation difference further comprises:
calculating the average resistivity change degree of the points inside the connected domain and the surrounding ore layer rock:
calculating the first possibility that the corresponding extremum region belongs to the mining area or the subsidence area according to the resistivity difference value variation difference and the average resistivity variation degree;
and if the first possibility is larger than a first preset threshold value, judging that the extremum area belongs to the goaf or the subsidence area.
8. The electrical prospecting techniques based mine geophysical data analysis method according to claim 7, wherein determining the regional nature of the extremum region from the resistivity variation difference further comprises:
Defining a minimum circumscribed rectangle of the extremum region;
Acquiring the height and the width of the rectangle, and acquiring the mass center of the communication domain corresponding to the extremum region according to the difference of the resistivity change amplitude of the extremum region in the vertical upward direction;
calculating the maximum value of the resistivity change gradient according to the centroid of the connected domain;
Calculating a second possibility that the extremum area belongs to the goaf according to the height, the width and the maximum value of the resistivity change gradient of the rectangle;
And if the second possibility is larger than a second preset threshold value, determining the extremum area as a goaf.
9. The method of claim 8, wherein if the regional property of the extremum region is a goaf, determining the location of the goaf based on the geophysical prospecting data further comprises:
And acquiring the corresponding connected domain of the goaf, obtaining the region range of the corresponding connected domain, and determining the position according to the depth and the width displayed in the resistivity profile.
10. The method for analyzing geophysical prospecting data for a mine based on the electrical prospecting technique according to claim 9, wherein said performing a secondary survey on the goaf based on the electrical prospecting technique according to the position of the goaf, determining the internal condition of the goaf further comprises:
Measuring again by using the high-frequency electromagnetic sounding method with the position as a central point, and determining the specific position of the final goaf according to the two results;
When the goaf is in a low-resistance abnormal region in the resistivity profile, water is filled in the goaf, and more mineral substances are contained in the water;
when the corresponding goaf is represented as a high-resistance abnormal region in the resistivity profile, the interior of the corresponding goaf is preserved intact and is not filled with water.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410552301.9A CN118131340B (en) | 2024-05-07 | 2024-05-07 | Mine geophysical prospecting data analysis method based on electrical prospecting technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410552301.9A CN118131340B (en) | 2024-05-07 | 2024-05-07 | Mine geophysical prospecting data analysis method based on electrical prospecting technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118131340A true CN118131340A (en) | 2024-06-04 |
| CN118131340B CN118131340B (en) | 2024-07-16 |
Family
ID=91244313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410552301.9A Active CN118131340B (en) | 2024-05-07 | 2024-05-07 | Mine geophysical prospecting data analysis method based on electrical prospecting technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118131340B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004301610A (en) * | 2003-03-31 | 2004-10-28 | Secom Co Ltd | Subterranean cavity detecting device |
| CN102520450A (en) * | 2011-11-16 | 2012-06-27 | 中国科学院地质与地球物理研究所 | Coal mine full-water goaf detection method |
| CN104965229A (en) * | 2015-06-10 | 2015-10-07 | 山东科技大学 | Method used for determining range of old mined-out area of coal mine |
| CN112462432A (en) * | 2020-11-16 | 2021-03-09 | 新疆维吾尔自治区煤田灭火工程局 | High-precision detection method and device for coal field fire area cavity |
| CN113550791A (en) * | 2021-07-15 | 2021-10-26 | 陕西麟北煤业开发有限责任公司 | Coal seam roof separation water detection method for coal mine control management |
| CN115980863A (en) * | 2022-12-16 | 2023-04-18 | 中国煤炭地质总局地球物理勘探研究院 | Exploration detection method and exploration device for geophysical engineering and application |
| CN116009111A (en) * | 2021-10-22 | 2023-04-25 | 中国石油化工股份有限公司 | High-density electrical method and micro-motion combined near-surface water-rich region exploration method |
| CN116299760A (en) * | 2023-02-20 | 2023-06-23 | 中国煤炭地质总局地球物理勘探研究院 | Fine processing method and system for detection data of coal mine goaf water accumulation area |
| JP2023154566A (en) * | 2022-04-07 | 2023-10-20 | 西日本技術開発株式会社 | Cavity determination device and cavity determination method and cavity determination program |
| CN117266933A (en) * | 2023-09-18 | 2023-12-22 | 国家能源投资集团有限责任公司 | Comprehensive monitoring system and method for coal mine goaf water migration storage simulation test |
| CN117974754A (en) * | 2022-10-19 | 2024-05-03 | 中国石油化工股份有限公司 | Method, electronic equipment and medium for quantitatively identifying reservoir karst cave based on image processing |
-
2024
- 2024-05-07 CN CN202410552301.9A patent/CN118131340B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004301610A (en) * | 2003-03-31 | 2004-10-28 | Secom Co Ltd | Subterranean cavity detecting device |
| CN102520450A (en) * | 2011-11-16 | 2012-06-27 | 中国科学院地质与地球物理研究所 | Coal mine full-water goaf detection method |
| CN104965229A (en) * | 2015-06-10 | 2015-10-07 | 山东科技大学 | Method used for determining range of old mined-out area of coal mine |
| CN112462432A (en) * | 2020-11-16 | 2021-03-09 | 新疆维吾尔自治区煤田灭火工程局 | High-precision detection method and device for coal field fire area cavity |
| CN113550791A (en) * | 2021-07-15 | 2021-10-26 | 陕西麟北煤业开发有限责任公司 | Coal seam roof separation water detection method for coal mine control management |
| CN116009111A (en) * | 2021-10-22 | 2023-04-25 | 中国石油化工股份有限公司 | High-density electrical method and micro-motion combined near-surface water-rich region exploration method |
| JP2023154566A (en) * | 2022-04-07 | 2023-10-20 | 西日本技術開発株式会社 | Cavity determination device and cavity determination method and cavity determination program |
| CN117974754A (en) * | 2022-10-19 | 2024-05-03 | 中国石油化工股份有限公司 | Method, electronic equipment and medium for quantitatively identifying reservoir karst cave based on image processing |
| CN115980863A (en) * | 2022-12-16 | 2023-04-18 | 中国煤炭地质总局地球物理勘探研究院 | Exploration detection method and exploration device for geophysical engineering and application |
| CN116299760A (en) * | 2023-02-20 | 2023-06-23 | 中国煤炭地质总局地球物理勘探研究院 | Fine processing method and system for detection data of coal mine goaf water accumulation area |
| CN117266933A (en) * | 2023-09-18 | 2023-12-22 | 国家能源投资集团有限责任公司 | Comprehensive monitoring system and method for coal mine goaf water migration storage simulation test |
Non-Patent Citations (2)
| Title |
|---|
| 严军;梁标;韩学长;刘亚南;: "瞬变电磁法在腾升煤矿探放水中的应用", 科技创新导报, no. 25, 1 September 2009 (2009-09-01) * |
| 马瑞花;蒋喜昆;: "综合电法勘探技术在探测奥陶岩溶水中的应用", 中国煤炭地质, no. 05, 25 May 2013 (2013-05-25) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118131340B (en) | 2024-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lu et al. | New insights into the carbonate karstic fault system and reservoir formation in the Southern Tahe area of the Tarim Basin | |
| Tian et al. | Multi-layered Ordovician paleokarst reservoir detection and spatial delineation: A case study in the Tahe Oilfield, Tarim Basin, Western China | |
| CN106526693B (en) | crack identification method and device | |
| US10310136B2 (en) | Lateral placement and completion design for improved well performance of unconventional reservoirs | |
| CN104597511B (en) | A kind of multilayer goaf ground tunnel transient electromagnetic detecting method | |
| CN110133737B (en) | Electromagnetism prediction method for blind mine | |
| CN112965141B (en) | Delineating method of ore-forming favorable section of uranium polymetallic ore | |
| CN113514886B (en) | Geological-seismic three-dimensional prediction method for beneficial part of sandstone-type uranium deposit mineralization | |
| Zhang et al. | Architecture characterization of Ordovician fault-controlled paleokarst carbonate reservoirs in Tuoputai, Tahe oilfield, Tarim Basin, NW China | |
| CN109425899B (en) | Method and device for predicting distribution of fractured zone of carbonate rock | |
| KR101495836B1 (en) | System for 3-dimensional electrical resistivity survey and surveying method for caving in section of tunnel using the same | |
| CN114386464B (en) | Deep learning extraction method for transient electromagnetic excitation information | |
| Liang et al. | Tectonic stress estimation from ultrasonic borehole image logs in a coal bed methane well, northeastern Qinshui Basin, China | |
| CN118131340B (en) | Mine geophysical prospecting data analysis method based on electrical prospecting technology | |
| Golebiowski et al. | Ambiguities in geophysical interpretation during fracture detection—case study from a limestone quarry (Lower Silesia Region, Poland) | |
| CN105700040A (en) | Bauxite positioning method and apparatus | |
| CN110118992B (en) | Method for exploring coal resources of fully-concealed deep coal field | |
| CN117492102A (en) | Detection method for ground water environment polluted range of in-situ uranium mine based on physical and electrical detection | |
| CN109991667B (en) | Method for quickly judging hydraulic connection between water-containing layers | |
| Marinenko et al. | Solving direct problems of electrical resistivity tomography for media with high-conductivity irregular-shaped heterogeneities by an example of a multiple well platform | |
| CN111208571B (en) | Well-ground combined detection method for multilayer mining water accumulation area | |
| CN104749645B (en) | A kind of method and apparatus of coal seam substrate interface undulation form detection in ash top difficult to understand | |
| Francke et al. | Preliminary evaluation of gpr for nickel laterite exploration | |
| Jageler et al. | Use of well logs and dipmeters in stratigraphic-trap exploration: geologic exploration methods | |
| Gönenç et al. | Geophysical Investigation of the Hematite Zones in Koçarlı-Demirtepe (Aydın/Turkey) |
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 |