Deep coal geologic structure detection method based on multi-source electrical method coupling

Abstract

The invention belongs to the technical field of coal geological exploration, and discloses a deep coal geological structure detection method based on multi-source electrical method coupling, which comprises the following steps: dividing a to-be-detected area into n sub-detection areas, and continuously detecting each sub-detection area for a period of time, then, collecting time-varying data of resistivity and time-varying data of polarizability of each detection area, and respectively calculating a resistivity difference coefficient and a polarizability difference coefficient of each area in combination with environment temperature and humidity of soil; the geological structure difference of each area is preliminarily judged according to the resistivity difference coefficient and the polarizability difference coefficient, then the geological structure difference of each area is further judged in combination with the resistivity difference coefficient and the polarizability difference coefficient, and the geological structure difference of the detection area is provided for exploration personnel. Explorers can make more sufficient preparations during exploration conveniently, and exploration efficiency and accuracy are improved.

Description

Technical Field

[0001] The invention belongs to the technical field of coal geological exploration, and in particular relates to a deep coal geological structure detection method based on multi-source electrical coupling. Background Art

[0002] As coal mining extends to deeper areas, accurately understanding the deep coal geological structure is crucial for the safe and efficient mining of coal resources.

[0003] The formation and distribution of coal resources are closely related to geological conditions. Geological differences may lead to changes in the thickness, stability and continuity of coal seams, thereby affecting the reserve assessment of coal resources. Therefore, the influence of geological conditions needs to be fully considered in the process of reserve assessment to ensure the accuracy of the assessment results. Most of the existing deep coal geological structure detection methods based on multi-source electrical coupling monitor the geological structure in real time through monitoring equipment, and judge the coal geological structure through the real-time value of the parameters. This method of monitoring has poor sensitivity and cannot achieve comprehensive fusion and efficient use of information, provide more on-site information to explorers, and help explorers to improve exploration technology and improve the accuracy of reserve prediction. Therefore, the present invention proposes a deep coal geological structure detection method based on multi-source electrical coupling, which aims to improve the efficiency and accuracy of coal geological structure detection. Summary of the invention

[0004] The purpose of the present invention is to provide a deep coal geological structure detection method based on multi-source electrical coupling to solve the problems faced in the above-mentioned background technology.

[0005] The purpose of the present invention can be achieved through the following technical solutions:

[0006] A method for detecting deep coal geological structures based on multi-source electrical coupling, the method comprising the following steps:

[0007] Step S1, dividing the area to be detected into n sub-detection areas, using geological radar to detect shallow geological structures in each detection area, and preliminarily delineating the approximate range of the deep target layer;

[0008] Step S2, arranging multi-source electrical detection electrodes in n sub-detection areas respectively;

[0009] Step S3, start the multi-source electrical detection equipment, continue the detection time for Δt hours, and collect the underground medium resistivity distribution data and polarizability data of each detection area respectively;

[0010] Step S4: real-time monitoring of temperature and humidity data of the area to be detected;

[0011] Step S5: Perform fusion analysis on the resistivity distribution data and the polarizability data, and judge the geological structure based on the analysis results.

[0012] As a further description of the technical solution of the present invention, the specific working process of step S3 includes:

[0013] Obtaining the resistivity variation data of the i-th detection area over time and the polarizability variation data of the i-th detection area over time during detection;

[0014] According to the resistivity variation data of the i-th detection area and the polarizability variation data of the i-th detection area, a resistivity variation function of the i-th detection area and a polarizability variation function of the i-th detection area are respectively constructed;

[0015] Among them, i belongs to n.

[0016] As a further description of the technical solution of the present invention, the specific working process of step S5 includes:

[0017] Construct a calculation model for the resistivity difference coefficient of the i-th detection area, and the expression is:

[0018]

[0019] Construct a calculation model for the polarizability difference coefficient of the i-th detection area, and the expression is:

[0020]

[0021] In the formula, ρ i is the resistivity difference coefficient of the i-th detection area, σ i is the polarizability difference coefficient of the i-th detection area, R i (t) is the resistivity function of the i-th detection area over time, W i (t) is the time-varying function of the polarizability of the ith detection area, t1 is the initial acquisition time, t2 is the end acquisition time, t2-t1=Δt, T i is the soil temperature of the ith detection area, H i is the soil moisture of the ith detection area, T0 is the standard soil temperature preset by the system, H0 is the standard soil moisture preset by the system, and α and β are the weight coefficients corresponding to temperature and humidity respectively.

[0022] As a further description of the technical solution of the present invention, the resistivity difference coefficient ρ of the i-th detection area is respectively i and the polarizability difference coefficient σ of the i-th detection area iCompared with the set corresponding threshold interval, if any difference coefficient does not meet the corresponding threshold interval, it means that the geological structure of the i-th detection area is quite different from that of other areas. If there is no difference coefficient that does not meet the corresponding threshold interval, the soil difference is further analyzed.

[0023] As a further description of the technical solution of the present invention, the specific working process of further analyzing the soil differences includes:

[0024] Construct the calculation model of the overall difference coefficient of the i-th detection area, and the expression is:

[0025]

[0026] In the formula, is the overall difference coefficient of the i-th detection area, [ρ1, ρ2] is the threshold interval corresponding to the resistivity difference coefficient preset by the system, [σ1, σ2] is the threshold interval corresponding to the polarizability difference coefficient preset by the system, δ and θ are the weight coefficients corresponding to the resistivity and polarizability, respectively.

[0027] As a further description of the technical solution of the present invention, the overall difference coefficient of the i-th detection area is Compared with the overall difference coefficient threshold interval set by the system, if the overall difference coefficient of the i-th detection area If the comparison does not fall within the threshold interval of the overall difference coefficient set by the system, it means that the i-th detection area is quite different from other areas.

[0028] As a further description of the technical solution of the present invention, the specific working process of step S5 also includes:

[0029] The i-th detection area is divided into m sub-areas. According to the resistivity variation function of the i-th detection area and the polarizability variation function of the i-th detection area, the resistivity R of the j-th sub-area is obtained respectively. ij and the polarizability W of the jth subregion ij ;

[0030] The resistivity R of the jth sub-region is ij and the polarizability W of the jth subregion ij Compared with the corresponding standard values ​​R0 and W0 set by the system, the geological structure of the i-th detection area is analyzed according to the comparison results, where j belongs to m.

[0031] As a further description of the technical solution of the present invention, the specific process of analyzing the geological structure of the i-th detection area according to the comparison result includes:

[0032] If the resistivity of the jth sub-region R ijIf the resistivity of the jth sub-region R ij The polarizability W of the jth sub-region is greater than or equal to R0. ij If it is greater than W0, it is judged that the coal seam may contain a large amount of metal sulfide electronic conductor minerals;

[0033] If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is less than W0, it is judged that the area may be a common rock formation;

[0034] If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is greater than W0, it is judged that there may be underground water flow channels or special hydrogeological conditions in the area.

[0035] As a further description of the scheme of the present invention, the specific working process of step S5 also includes: converting the resistivity variation function of the i-th detection area with time and the polarizability variation function of the i-th detection area with time into the resistivity variation function of the i-th detection area with sub-region and the polarizability variation function of the i-th detection area, and then constructing a rectangular coordinate system xoy, and respectively producing the resistivity variation function curve of the i-th detection area with sub-region and the polarizability variation function curve of the i-th detection area in the coordinate system, and then further judging the geological structure of the i-th detection area according to the resistivity variation function curve of the i-th detection area with sub-region and the polarizability variation function curve of the i-th detection area;

[0036] As a further description of the solution of the present invention, the process of further determining the geological structure of the i-th detection area includes:

[0037] When the resistivity of a certain area suddenly decreases and the polarizability increases, it is likely that this is a fault fracture zone. This is because the presence of groundwater reduces the resistivity, while the filling of minerals increases the polarizability. On the contrary, if the resistivity suddenly changes but the polarizability does not change significantly, the fault may be mainly filled with water, with less mineral content.

[0038] Beneficial effects of the present invention:

[0039] 1. The present invention divides the area to be detected into n sub-detection areas, detects each sub-detection area for a continuous detection time of Δt hours, and then collects the resistivity variation data and polarizability variation data of each detection area over time, and then calculates the resistivity difference coefficient and polarizability difference coefficient of each area according to the resistivity variation data and polarizability variation data of each detection area and the environmental temperature and humidity of the soil, and makes a preliminary judgment on the geological structure difference of each area according to the resistivity difference coefficient and polarizability difference coefficient, and then makes a further judgment on the geological structure difference of each area in combination with the resistivity difference coefficient and polarizability difference coefficient, and provides the geological structure difference of the detection area to the exploration personnel, so that the exploration personnel can make more adequate preparations during exploration and improve the exploration efficiency and accuracy. ;

[0040] 2. The present invention uses a variety of electrical techniques to overcome the limitations of a single electrical method in deep detection. Through data fusion and collaborative analysis, it significantly improves the resolution of deep coal geological structure detection and can accurately identify small-scale structures and stratigraphic anomalies.

[0041] Of course, any product implementing the present invention does not necessarily need to achieve all of the advantages described above at the same time. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for describing the embodiments are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0043] Figure 1 It is a partial flow chart of the deep coal geological structure detection method based on multi-source electrical coupling of the present invention. DETAILED DESCRIPTION

[0044] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0045] See also Figure 1 As shown, a method for detecting deep coal geological structures based on multi-source electrical coupling is disclosed, and the method comprises the following steps:

[0046] Step S1, dividing the area to be detected into n sub-detection areas, using geological radar to detect shallow geological structures in each detection area, and preliminarily delineating the approximate range of the deep target layer;

[0047] Step S2, arranging multi-source electrical detection electrodes in n sub-detection areas respectively;

[0048] Step S3, start the multi-source electrical detection equipment, continue the detection time for Δt hours, and collect the underground medium resistivity distribution data and polarizability data of each detection area respectively;

[0049] Step S4: real-time monitoring of temperature and humidity data of the area to be detected;

[0050] Step S5: Perform fusion analysis on the resistivity distribution data and the polarizability data, and judge the geological structure based on the analysis results.

[0051] Through the above technical scheme, the present invention divides the area to be detected into n sub-detection areas, and continuously detects each sub-detection area for Δt hours, then collects the resistivity change data and polarizability change data of each detection area over time, and then calculates the resistivity difference coefficient and polarizability difference coefficient of each area according to the resistivity change data and polarizability change data of each detection area and the ambient temperature and humidity of the soil, and makes a preliminary judgment on the geological structure difference of each area according to the resistivity difference coefficient and the polarizability difference coefficient, and then makes a further judgment on the geological structure difference of each area in combination with the resistivity difference coefficient and the polarizability difference coefficient, and provides the geological structure difference of the detection area to the exploration personnel, so that the exploration personnel can make more adequate preparations during exploration and improve the exploration efficiency and accuracy.

[0052] The specific working process of step S3 includes:

[0053] Obtaining the resistivity variation data of the i-th detection area over time and the polarizability variation data of the i-th detection area over time during detection;

[0054] According to the resistivity variation data of the i-th detection area and the polarizability variation data of the i-th detection area, a resistivity variation function of the i-th detection area and a polarizability variation function of the i-th detection area are respectively constructed;

[0055] Among them, i belongs to n.

[0056] The specific working process of step S5 includes:

[0057] Construct a calculation model for the resistivity difference coefficient of the i-th detection area, and the expression is:

[0058]

[0059] Construct a calculation model for the polarizability difference coefficient of the i-th detection area, and the expression is:

[0060]

[0061] In the formula, ρ i is the resistivity difference coefficient of the i-th detection area, σ i is the polarizability difference coefficient of the i-th detection area, R i (t) is the resistivity function of the i-th detection area over time, W i (t) is the time-varying function of the polarizability of the ith detection area, t1 is the initial acquisition time, t2 is the end acquisition time, t2-t1=Δt, T i is the soil temperature of the ith detection area, H i is the soil moisture of the ith detection area, T0 is the standard soil temperature preset by the system, H0 is the standard soil moisture preset by the system, and α and β are the weight coefficients corresponding to temperature and humidity respectively.

[0062] The resistivity difference coefficient ρ of the i-th detection area is i and the polarizability difference coefficient σ of the i-th detection area i Compared with the set corresponding threshold interval, if any difference coefficient does not meet the corresponding threshold interval, it means that the geological structure of the i-th detection area is quite different from that of other areas. If there is no difference coefficient that does not meet the corresponding threshold interval, the soil difference is further analyzed.

[0063] Through the above technical solution, this embodiment provides a method for judging the geological differences of the detection area. First, the resistivity change function of the i-th detection area and the polarizability change function of the i-th detection area are obtained respectively, and then the temperature and humidity of the soil in the detection area are obtained in real time to eliminate the influence of environmental factors on the detection signal. Then, the formulas are used to respectively Calculate the resistivity difference coefficient and formula of the i-th detection area Calculate the polarizability difference coefficient of the i-th detection area. Finally, the resistivity difference coefficient ρ of the i-th detection area is i and the polarizability difference coefficient σ of the i-th detection area i Compared with the set corresponding threshold interval, if any difference coefficient does not satisfy the corresponding threshold interval, it means that the geological structure of the i-th detection area is quite different from that of other areas.

[0064] The specific working process of further analyzing soil differences includes:

[0065] Construct the calculation model of the overall difference coefficient of the i-th detection area, and the expression is:

[0066]

[0067] In the formula, is the overall difference coefficient of the i-th detection area, [ρ1, ρ2] is the threshold interval corresponding to the resistivity difference coefficient preset by the system, [σ1, σ2] is the threshold interval corresponding to the polarizability difference coefficient preset by the system, δ and θ are the weight coefficients corresponding to the resistivity and polarizability, respectively.

[0068] The overall difference coefficient of the i-th detection area Compared with the overall difference coefficient threshold interval set by the system, if the overall difference coefficient of the i-th detection area If the comparison does not fall within the threshold interval of the overall difference coefficient set by the system, it means that the i-th detection area is quite different from other areas.

[0069] Through the above technical solution, this embodiment combines the resistivity difference coefficient and the polarizability difference coefficient to judge the overall geological difference of the detection area, and uses the formula Calculate the overall difference coefficient of the ith detection area and convert the overall difference coefficient of the ith detection area Compared with the overall difference coefficient threshold interval set by the system, if the overall difference coefficient of the i-th detection area If the comparison does not fall within the threshold interval of the overall difference coefficient set by the system, it means that the i-th detection area is quite different from other areas.

[0070] The specific working process of step S5 also includes:

[0071] The i-th detection area is divided into m sub-areas. According to the resistivity variation function of the i-th detection area and the polarizability variation function of the i-th detection area, the resistivity R of the j-th sub-area is obtained respectively. ij and the polarizability W of the jth subregion ij ;

[0072] The resistivity R of the jth sub-region is ij and the polarizability W of the jth subregion ij Compared with the corresponding standard values ​​R0 and W0 set by the system, the geological structure of the i-th detection area is analyzed according to the comparison results, where j belongs to m.

[0073] The specific process of analyzing the geological structure of the i-th detection area according to the comparison results includes:

[0074] If the resistivity of the jth sub-region R ij If the resistivity of the jth sub-region R ij The polarizability W of the jth sub-region is greater than or equal to R0. ij If it is greater than W0, it is judged that the coal seam may contain a large amount of metal sulfide electronic conductor minerals;

[0075] If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is less than W0, it is judged that the area may be a common rock formation;

[0076] If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is greater than W0, it is judged that there may be underground water flow channels or special hydrogeological conditions in the area;

[0077] The present invention can also convert the resistivity variation function of the i-th detection area with time and the polarizability variation function of the i-th detection area with time into the resistivity variation function of the i-th detection area with sub-area and the polarizability variation function of the i-th detection area, and then construct a rectangular coordinate system xoy, and respectively generate the resistivity variation function curve of the i-th detection area with sub-area and the polarizability variation function curve of the i-th detection area in the coordinate system, and then further judge the geological structure of the i-th detection area according to the resistivity variation function curve of the i-th detection area with sub-area and the polarizability variation function curve of the i-th detection area;

[0078] When the resistivity of a certain area suddenly decreases and the polarizability increases, it is likely that this is a fault fracture zone. This is because the presence of groundwater reduces the resistivity, while the filling of minerals increases the polarizability. On the contrary, if the resistivity suddenly changes but the polarizability does not change significantly, the fault may be mainly filled with water, with less mineral content.

[0079] The above contents are merely examples and explanations of the concept of the present invention. The technicians in this technical field may make various modifications or additions to the specific embodiments described or replace them in a similar manner. As long as they do not deviate from the concept of the invention or exceed the scope defined by the claims, they should all fall within the protection scope of the present invention.

Claims

1. A method for detecting deep coal geological structures based on multi-source electrical coupling, characterized in that: The method comprises the following steps: Step S1, dividing the area to be detected into n sub-detection areas, using geological radar to detect shallow geological structures in each detection area, and preliminarily delineating the approximate range of the deep target layer; Step S2, arranging multi-source electrical detection electrodes in n sub-detection areas respectively; Step S3, start the multi-source electrical detection equipment, continue the detection time for Δt hours, and collect the resistivity distribution data and polarizability data of the underground medium in each detection area; Step S4: real-time monitoring of temperature and humidity data of the area to be detected; Step S5: Perform fusion analysis on the resistivity distribution data and the polarizability data, and judge the geological structure based on the analysis results.

2. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 1 is characterized in that: The specific working process of step S3 includes: Obtaining the resistivity variation data of the i-th detection area over time and the polarizability variation data of the i-th detection area over time during detection; According to the resistivity variation data of the i-th detection area and the polarizability variation data of the i-th detection area, a resistivity variation function of the i-th detection area and a polarizability variation function of the i-th detection area are respectively constructed; Among them, i belongs to n.

3. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 1 is characterized in that: The specific working process of step S5 includes: Construct a calculation model for the resistivity difference coefficient of the i-th detection area, and the expression is: Construct a calculation model for the polarizability difference coefficient of the i-th detection area, and the expression is: In the formula, ρ i is the resistivity difference coefficient of the i-th detection area, σ i is the polarizability difference coefficient of the i-th detection area, R i (t) is the resistivity function of the i-th detection area over time, W i (t) is the time-varying function of the polarizability of the ith detection area, t1 is the initial acquisition time, t2 is the end acquisition time, t2-t1=Δt, T i is the soil temperature of the ith detection area, H i is the soil moisture of the ith detection area, T0 is the standard soil temperature preset by the system, H0 is the standard soil moisture preset by the system, and α and β are the weight coefficients corresponding to temperature and humidity respectively.

4. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 3 is characterized in that: The resistivity difference coefficient ρ of the i-th detection area is i and the polarizability difference coefficient σ of the i-th detection area i Compared with the set corresponding threshold interval, if any difference coefficient does not meet the corresponding threshold interval, it means that the geological structure of the i-th detection area is quite different from that of other areas. If there is no difference coefficient that does not meet the corresponding threshold interval, the soil difference is further analyzed.

5. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 4 is characterized in that: The specific working process of further analyzing soil differences includes: Construct the calculation model of the overall difference coefficient of the i-th detection area, and the expression is: In the formula, is the overall difference coefficient of the i-th detection area, [ρ1, ρ2] is the threshold interval corresponding to the resistivity difference coefficient preset by the system, [σ1, σ2] is the threshold interval corresponding to the polarizability difference coefficient preset by the system, δ and θ are the weight coefficients corresponding to the resistivity and polarizability, respectively.

6. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 4 is characterized in that: The overall difference coefficient of the i-th detection area Compared with the overall difference coefficient threshold interval set by the system, if the overall difference coefficient of the i-th detection area If the comparison does not fall within the threshold interval of the overall difference coefficient set by the system, it means that the i-th detection area is quite different from other areas.

7. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 1, characterized in that: The specific working process of step S5 also includes: The i-th detection area is divided into m sub-areas. According to the resistivity variation function of the i-th detection area and the polarizability variation function of the i-th detection area, the resistivity R of the j-th sub-area is obtained respectively. ij and the polarizability W of the jth subregion ij ; The resistivity R of the jth sub-region is ij and the polarizability W of the jth subregion ij Compared with the corresponding standard values ​​R0 and W0 set by the system, the geological structure of the i-th detection area is analyzed according to the comparison results, where j belongs to m.

8. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 7, characterized in that: The specific process of analyzing the geological structure of the i-th detection area according to the comparison results includes: The specific process of analyzing the geological structure of the i-th detection area according to the comparison results includes: If the resistivity of the jth sub-region R ij If the resistivity of the jth sub-region R ij The polarizability W of the jth sub-region is greater than or equal to R0. ij If it is greater than W0, it is judged that the coal seam may contain a large amount of metal sulfide electronic conductor minerals; If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is less than W0, it is judged that the area may be a common rock formation; If the resistivity of the jth sub-region R ij Less than R0 and the polarizability W of the jth sub-region ij If it is greater than W0, it is judged that there may be underground water flow channels or special hydrogeological conditions in the area.

9. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 8, characterized in that: The specific working process of step S5 also includes: converting the resistivity variation function of the i-th detection area with time and the polarizability variation function of the i-th detection area with time into the resistivity variation function of the i-th detection area with sub-area and the polarizability variation function of the i-th detection area, and then constructing a rectangular coordinate system xoy, and generating the resistivity variation function curve of the i-th detection area with sub-area and the polarizability variation function curve of the i-th detection area in the coordinate system respectively, and then further judging the geological structure of the i-th detection area according to the resistivity variation function curve of the i-th detection area with sub-area and the polarizability variation function curve of the i-th detection area.

10. The method for detecting deep coal geological structures based on multi-source electrical coupling according to claim 9, characterized in that: The process of further judging the geological structure of the i-th detection area includes: When the resistivity of a certain area suddenly decreases and the polarizability increases, it is likely that this is a fault fracture zone. This is because the presence of groundwater reduces the resistivity, while the filling of minerals increases the polarizability. On the contrary, if the resistivity suddenly changes but the polarizability does not change significantly, the fault may be mainly filled with water, with less mineral content.

Images