CN110850472A - Variable offset distance advanced fault detection method based on shock wave excitation seismic source - Google Patents

Variable offset distance advanced fault detection method based on shock wave excitation seismic source Download PDF

Info

Publication number
CN110850472A
CN110850472A CN201910994193.XA CN201910994193A CN110850472A CN 110850472 A CN110850472 A CN 110850472A CN 201910994193 A CN201910994193 A CN 201910994193A CN 110850472 A CN110850472 A CN 110850472A
Authority
CN
China
Prior art keywords
seismic
seismic source
component
fault
shock wave
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
Application number
CN201910994193.XA
Other languages
Chinese (zh)
Other versions
CN110850472B (en
Inventor
王勃
孙华超
刘盛东
龚震
李兴兴
王志丹
邢世雨
丁昕
张恒
张振东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201910994193.XA priority Critical patent/CN110850472B/en
Publication of CN110850472A publication Critical patent/CN110850472A/en
Application granted granted Critical
Publication of CN110850472B publication Critical patent/CN110850472B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02—Generating seismic energy
    • G01V1/04—Details
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/20—Arrangements of receiving elements, e.g. geophone pattern

Abstract

The invention discloses a variable offset distance advanced fault detection method based on shock wave excitation seismic sources, which comprises the steps of exciting the seismic sources by adopting shock waves, placing the shock wave excitation seismic sources into a drill hole, utilizing the characteristic that the shock waves excite the seismic sources for multiple times of shock excitation, and respectively exciting by changing different hole depths so as to generate a plurality of seismic sources, wherein the offset distances between the plurality of seismic sources and a detector are different, a common shot point gather is obtained by a seismic observation system after each time of excitation, then each common shot point gather is processed, and the position of a maximum amplitude reflection point of a reflection channel wave generated by each seismic source on a fault plane of a local disk is obtained; finally, obtaining a regression linear equation through the position of each maximum amplitude reflection point, and displaying a straight line of the regression linear equation on an XOY surface of a coordinate system, wherein the position of the straight line is the position of a fault; therefore, the staff can visually acquire the fault strike parameters.

Description

Variable offset distance advanced fault detection method based on shock wave excitation seismic source
Technical Field
The invention relates to a mine geophysical detection method, in particular to a variable offset distance advanced fault detection method based on a shock wave excitation seismic source.
Background
Frequent coal roadway tunneling accidents always disturb the safety production of coal mines, wherein water inrush or coal and gas outburst disaster accidents caused by geological structures are mainly caused by the fact that the fault detection precision in front of the head is poor; the problem of 'the diameter of a hole' in underground drilling cannot fully understand the front geological structure, and the cost is high; the geophysical method can detect the front fault structure, wherein the seismic advanced detection method is the main method for coal mine geological exploration at present; at present, fault fall, inclination, dip angle and other parameter technologies are solved, but the fault strike parameters are difficult to calculate due to the limitation of roadway seismic advanced detection minimum offset aperture.
The seismic source in the seismic advance detection is an important link for determining whether the seismic advance detection is successful or not; the energy of a traditional percussion hammer seismic source is weak; although the explosive seismic source can meet the requirement of earthquake advanced detection, the explosive seismic source is limited to coal mine conditions and the approval process is complex, so that the development of a novel seismic source is imminent; the application number CN201811061977.9 patent provides a novel CO used for the exploration method of earthquake transmission groove waves of high gas mines2The cannon can meet the requirement of underground exploration of a coal mineThe seismic source is a clean seismic source, and the limitation of a high gas mine to the explosive seismic source is solved; but CO2The cannon excitation instrument is subjected to CO2The reservoir is limited and difficult to circulate through the borehole.
The method and the system for the advanced detection of the groove wave of the underground single-head roadway based on the hole-roadway combination, which are proposed by the patent No. CN108121010B, improve the accuracy of the groove wave reflected by the roadway by arranging the wave detectors in holes and arranging the wave detectors at left and right sides of the roadway in left and right coal seams, but the method needs 2 deep-hole wave detectors, has large workload, and needs to research a new method to reduce the number of drilled holes so as to improve the efficiency; in addition, the method also needs to manually explain the fault strike parameters according to the offset imaging section, the explanation difference of personnel with different technical levels is large, the explanation deviation is easy to cause, and a new method with intuitive detection results and easy identification is urgently needed.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a variable offset distance advanced fault detection method based on a shock wave excitation seismic source, which can display the fault position in a coordinate system by exciting the seismic source for multiple times, so that a worker can intuitively acquire fault trend parameters and the development of subsequent work is facilitated.
In order to achieve the purpose, the invention adopts the technical scheme that: a fault detection method based on shock wave excitation seismic source variable offset distance advance comprises the following specific steps:
the method comprises the following steps: constructing a drill hole at the position 2m away from the head of the coal roadway right side, wherein the length of the drill hole is N meters, and the drill hole is positioned in the middle of the coal seam and is parallel to a coal roadway top bottom plate;
step two: sending a shock wave excitation seismic source into the deepest part of the drilled hole;
step three: arranging m mounting holes on the right side of the coal roadway, wherein the distance between every two adjacent mounting holes is 3 meters, the positions of the mounting holes and the positions of drill holes are on the same horizontal plane, and a three-component detector is arranged in each mounting hole; each three-component detector is connected with an earthquake recorder, and a coal roadway tunneling advanced prediction observation system is established;
step four: shock wave excited vibrationThe source is excited for the first time from the deepest part of the borehole and then is excited every 2m from the inside to the outside for a total of n times to form a seismic source S1、S2、S3…SnEach three-component detector receives seismic waves generated by each seismic source in real time and transmits the seismic waves to a seismic recorder;
step five: establishing a three-dimensional coordinate system along a coal roadway, wherein the X direction is the direction towards the head, the Y direction is perpendicular to the coal wall and perpendicular to the direction towards the head, the Z direction is perpendicular to the top bottom plate of the roadway and perpendicular to the direction towards the head, the center of the roadway surface on the left boundary of the model is used as an origin to establish the three-dimensional coordinate system, and each seismic source position and each three-component detector position are included in the three-dimensional coordinate system;
step six: each seismic source generates a three-component common shot point gather, n three-component common shot point gathers are generated in total, each common shot point gather is analyzed, a reflection slot wave signal from the fault plane of the disc is selected from the common shot point gathers, and the maximum amplitude reflection point of each common shot point gather on the fault plane of the disc is obtained and is A1，A2，A3…An；
Step seven: all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinates of the reflection points are A1(x1，y1)、A2(x2，y2)，A3(x3，y3)…An(xn，yn) Solving a regression linear equation of the linear equation;
step eight: and forming a straight line on the XOY plane according to the regression linear equation obtained in the step seven, namely the fault position, so as to determine the fault trend.
Further, the specific process of the sixth step is as follows:
(1) extracting seismic source S1The generated common shot gather selects a time window t with three-component reflected channel wave signals1；
(2) From the time window t1Selecting three-component seismic waves of the duration t of the Ehry phase from the reflected channel wave signals to analyze;
(3) and (3) calculating the three-component seismic waves with the duration t in the step (2) to obtain the polarization parameters of the three-component seismic waves, and specifically comprising the following steps:
i, carrying out Hilbert transformation on X, Y, Z three-component seismic signals as follows:
Figure BDA0002239236520000031
Figure BDA0002239236520000032
in the formula: x (t), y (t), z (t) X, Y, Z components and symbols respectively
Figure BDA0002239236520000034
Which represents the Hilbert transform, is,
II, establishing a Hermitian matrix structure:
C(t)＝M*(t)·M(t)
wherein m (t) ([ hx (t)' hy (t) ()) hz (t) ], and the symbol denotes a complex conjugate transpose of the matrix;
III, solving to obtain the maximum characteristic value lambda1With main polarization direction centered on the eigenvalue λ1In the direction of the characteristic vector, the normalized characteristic vector is obtained and expressed as [ x ]1(t)，y1(t)，z1(t)]；
IV, the azimuth angle of the main polarization direction is expressed as:
Figure BDA0002239236520000036
v, carrying out amplitude calculation according to the parameters, wherein the formula is as follows:
where Ω represents all traces of n sources corresponding to m three-component detectors, A (t, f) represents a stacking function of scalar amplitudes at a certain frequencyiRepresents the amplitude energy, Ls+LRRepresenting the total path of propagation of the seismic waves, VCRepresenting the wave velocity of the trough wave airy phase, and x, y and z represent the positions of each seismic source and the three-component detector; weight function C of the main polarization directioniThe expression is as follows:
Figure BDA0002239236520000038
through CiLeft and right artifacts in amplitude calculation can be eliminated, so that whether the position of a reflection point is in the positive Y direction or the negative Y direction is directly determined;
(6) through the above (1) to (5), the position of the maximum amplitude reflection point of the reflection groove wave on the sectional plane of the disc can be determined and recorded as A1；
(7) Analyzing the common shot gather of other seismic sources, repeating the steps (1) to (6) to obtain other maximum amplitude reflection points A2，A3，A4…An。
Further, the specific process of the seventh step is as follows:
(1) all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinate is A1(x1，y1)、A2(x2，y2)、A3(x3，y3)…An(xn，yn)；
(2) Respectively obtaining:
Figure BDA0002239236520000041
Figure BDA0002239236520000042
Figure BDA0002239236520000043
Figure BDA0002239236520000045
Figure BDA0002239236520000046
(3) obtaining a regression linear equation through the parameters in the step (2):
Figure BDA0002239236520000047
thereby ensuring that as many reflection points as possible are on or near the straight line.
Compared with the prior art, the method adopts shock waves as seismic sources, the shock wave seismic sources are placed in the drill holes, the multiple excitation of the shock wave seismic sources is utilized, the multiple excitation is respectively carried out by changing different hole depths, so that a plurality of seismic sources are generated, the offset distances between the plurality of seismic sources and the detectors are different, after each excitation, the seismic observation system acquires a common shot point gather once, then each common shot point gather is processed, and the maximum amplitude reflection point position of reflection channel waves generated by each seismic source on a fault interface is acquired; finally, obtaining a regression linear equation through the position of each maximum amplitude reflection point, and displaying a straight line of the regression linear equation on an XOY surface of a coordinate system, wherein the position of the straight line is the position of a fault; therefore, the staff can visually know the detected fault position, and the subsequent work can be conveniently carried out. In addition, the seismic source adopted by the invention has fewer and safer inspection procedures than the traditional explosive seismic source, and is more CO-safe than the traditional explosive seismic source2The cannon seismic source can be recycled for multiple times; the method only needs to construct one drilling hole, and the efficiency is doubled compared with the prior art.
Drawings
FIG. 1 is a schematic view of the arrangement of the observation system in the present invention;
FIG. 2 is a scatter plot of a reflection point image obtained using the present invention.
In the figure: 1. the system comprises a coal roadway, a three-component detector, a detector 3, a head-on detector, a shock wave excitation seismic source 4 and a fault.
Detailed Description
The present invention will be further explained below.
As shown in the figure, the method comprises the following specific steps:
the method comprises the following steps: constructing a drill hole at the position of the right side of the coal roadway 1, which is 2m away from the head, wherein the length of the drill hole is N meters, and the drill hole is positioned in the middle of the coal seam and is parallel to the top floor of the coal roadway;
step two: sending the shock wave excitation seismic source 4 into the deepest part of the drilled hole;
step three: the method comprises the following steps that m mounting holes are distributed on the right side of a coal roadway 1, the distance between every two adjacent mounting holes is 3 m, the positions of all the mounting holes and a drill hole are located on the same horizontal plane, and a three-component detector 2 is uniformly distributed in each mounting hole; each three-component detector 2 is connected with an earthquake recorder, and a coal roadway tunneling advanced prediction observation system is established;
step four: the shock wave excitation seismic source 4 carries out first excitation from the deepest part of the drill hole, and then the shock wave excitation seismic source is excited once every 2m from inside to outside for n times in total to form a seismic source S1、S2、S3…SnEach three-component wave detector 2 receives seismic waves generated by each seismic source in real time and transmits the seismic waves to a seismic recorder;
step five: establishing a three-dimensional coordinate system along a coal roadway, wherein the X direction is the direction towards the head, the Y direction is perpendicular to the coal wall and perpendicular to the direction towards the head, the Z direction is perpendicular to the top bottom plate of the roadway and perpendicular to the direction towards the head, the center of the roadway surface on the left boundary of the model is used as an origin to establish the three-dimensional coordinate system, and each seismic source position and each three-component detector position are included in the three-dimensional coordinate system;
step six: each seismic source generates a three-component common shot point gather, n three-component common shot point gathers are generated in total, each common shot point gather is analyzed, a reflection slot wave signal from the fault plane of the disc is selected from the common shot point gathers, and each common shot point gather is solved in the local shot point gatherMaximum amplitude reflection points on the disc fault plane are respectively A1，A2，A3…An；
Step seven: all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinates of the reflection points are A1(x1，y1)、A2(x2，y2)，A3(x3，y3)…An(xn，yn) Solving a regression linear equation of the linear equation;
step eight: and forming a straight line on the XOY plane according to the regression linear equation obtained in the step seven, namely the position of the fault 5, so as to determine the trend of the fault 5.
Further, the specific process of the sixth step is as follows:
(1) extracting seismic source S1The generated common shot gather selects a time window t with three-component reflected channel wave signals1；
(2) From the time window t1Selecting three-component seismic waves of the duration t of the Ehry phase from the reflected channel wave signals to analyze;
(3) and (3) calculating the three-component seismic waves with the duration t in the step (2) to obtain the polarization parameters of the three-component seismic waves, and specifically comprising the following steps:
i, carrying out Hilbert transformation on X, Y, Z three-component seismic signals as follows:
Figure BDA0002239236520000061
Figure BDA0002239236520000062
Figure BDA0002239236520000063
in the formula: x (t), y (t), z (t) X, Y, Z components and symbols respectively
Figure BDA0002239236520000064
Which represents the Hilbert transform, is,
Figure BDA0002239236520000065
II, establishing a Hermitian matrix structure:
C(t)＝M*(t)·M(t)
wherein m (t) ([ hx (t)' hy (t) ()) hz (t) ], and the symbol denotes a complex conjugate transpose of the matrix;
III, solving to obtain the maximum characteristic value lambda1With main polarization direction centered on the eigenvalue λ1In the direction of the characteristic vector, the normalized characteristic vector is obtained and expressed as [ x ]1(t)，y1(t)，z1(t)]；
IV, the azimuth angle of the main polarization direction is expressed as:
Figure BDA0002239236520000066
v, carrying out amplitude calculation according to the parameters, wherein the formula is as follows:
Figure BDA0002239236520000067
where Ω represents all traces of n sources corresponding to m three-component detectors, A (t, f) represents a stacking function of scalar amplitudes at a certain frequencyiRepresents the amplitude energy, Ls+LRRepresenting the total path of propagation of the seismic waves, VCRepresenting the wave velocity of the trough wave airy phase, and x, y and z represent the positions of each seismic source and the three-component detector; weight function C of the main polarization directioniThe expression is as follows:
through CiLeft and right artifacts in amplitude calculation can be eliminated, so that whether the position of a reflection point is in the positive Y direction or the negative Y direction is directly determined;
(6) through the above (1) to (5), the position of the maximum amplitude reflection point of the reflection groove wave on the sectional plane of the disc can be determined and recorded as A1；
(7) Analyzing the common shot gather of other seismic sources, repeating the steps (1) to (6) to obtain other maximum amplitude reflection points A2，A3，A4…An。
Further, the specific process of the seventh step is as follows:
(1) all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinate is A1(x1，y1)、A2(x2，y2)、A3(x3，y3)…An(xn，yn)；
(2) Respectively obtaining:
Figure BDA0002239236520000072
Figure BDA0002239236520000073
Figure BDA0002239236520000075
Figure BDA0002239236520000076
Figure BDA0002239236520000077
(3) obtaining a regression linear equation through the parameters in the step (2):
Figure BDA0002239236520000078
thereby ensuring that as many reflection points as possible are on or near the straight line.

Claims (3)

1. A fault detection method based on shock wave excitation seismic source variable offset distance advance is characterized by comprising the following specific steps:
the method comprises the following steps: constructing a drill hole at the position 2m away from the head of the coal roadway right side, wherein the length of the drill hole is N meters, and the drill hole is positioned in the middle of the coal seam and is parallel to a coal roadway top bottom plate;
step two: sending a shock wave excitation seismic source into the deepest part of the drilled hole;
step three: arranging m mounting holes on the right side of the coal roadway, wherein the distance between every two adjacent mounting holes is 3 meters, the positions of the mounting holes and the positions of drill holes are on the same horizontal plane, and a three-component detector is arranged in each mounting hole; each three-component detector is connected with an earthquake recorder, and a coal roadway tunneling advanced prediction observation system is established;
step four: the shock wave excitation seismic source is excited for the first time from the deepest part of the drill hole, and then is excited every 2m from inside to outside for n times in total to form a seismic source S1、S2、S3...SnEach three-component detector receives seismic waves generated by each seismic source in real time and transmits the seismic waves to a seismic recorder;
step five: establishing a three-dimensional coordinate system along a coal roadway, wherein the X direction is the direction towards the head, the Y direction is perpendicular to the coal wall and perpendicular to the direction towards the head, the Z direction is perpendicular to the top bottom plate of the roadway and perpendicular to the direction towards the head, the center of the roadway surface on the left boundary of the model is used as an origin to establish the three-dimensional coordinate system, and each seismic source position and each three-component detector position are included in the three-dimensional coordinate system;
step six: each seismic source generates a three-component common shot point gather, n three-component common shot point gathers are generated in total, each common shot point gather is analyzed, a reflection slot wave signal from the fault plane of the disc is selected from the common shot point gathers, and the maximum amplitude reflection point of each common shot point gather on the fault plane of the disc is obtained and is A1，A2，A3...An；
Step seven: all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinates of the reflection points are A1(x1，y1)、A2(x2，y2)，A3(x3，y3)...An(xn，yn) Solving a regression linear equation of the linear equation;
step eight: and forming a straight line on the XOY plane according to the regression linear equation obtained in the step seven, namely the fault position, so as to determine the fault trend.
2. The method for variable offset advanced fault detection based on the shock wave excitation seismic source according to claim 1, wherein the specific process of the sixth step is as follows:
(1) extracting seismic source S1The generated common shot gather selects a time window t with three-component reflected channel wave signals1；
(2) From the time window t1Selecting three-component seismic waves of the duration t of the Ehry phase from the reflected channel wave signals to analyze;
(3) and (3) calculating the three-component seismic waves with the duration t in the step (2) to obtain the polarization parameters of the three-component seismic waves, and specifically comprising the following steps:
I. the Hilbert transform is performed on X, Y, Z three component seismic signals as follows:
Figure FDA0002239236510000022
Figure FDA0002239236510000023
in the formula: x (t), y (t), z (t) X, Y, Z components and symbols respectively
Figure FDA0002239236510000028
Which represents the Hilbert transform, is,
Figure FDA0002239236510000024
II. Establishing a Hermitian matrix structure:
C(t)＝M*(t)·M(t)
wherein m (t) ([ hx (t)' hy (t) ()) hz (t) ], and the symbol denotes a complex conjugate transpose of the matrix;
III, solving to obtain a maximum characteristic value lambda1With main polarization direction centered on the eigenvalue λ1In the direction of the characteristic vector, the normalized characteristic vector is obtained and expressed as [ x ]1(t)，y1(t)，z1(t)]；
IV, azimuth of the main polarization direction is expressed as:
Figure FDA0002239236510000025
v, carrying out amplitude calculation according to the parameters, wherein the formula is as follows:
where Ω represents all traces of n sources corresponding to m three-component detectors, A (t, f) represents a stacking function of scalar amplitudes at a certain frequencyiRepresents the amplitude energy, Ls+LRRepresenting the total path of propagation of the seismic waves, VCRepresenting the wave velocity of the trough wave airy phase, and x, y and z represent the positions of each seismic source and the three-component detector; weight function C of the main polarization directioniThe expression is as follows:
Figure FDA0002239236510000027
through CiLeft and right artifacts in amplitude calculation can be eliminated, so that whether the position of a reflection point is in the positive Y direction or the negative Y direction is directly determined;
(6) through the above (1) to (5), the position of the maximum amplitude reflection point of the reflection groove wave on the sectional plane of the disc can be determined and recorded as A1；
(7) The co-shot gathers of other seismic sources are analyzed,repeating the steps (1) to (6) to obtain other maximum amplitude reflection points A2，A3，A4…An。
3. The method for variable offset advanced fault detection based on shock wave excitation seismic source according to claim 1, wherein the concrete process of the seventh step is as follows:
(1) all maximum amplitude reflection points are classified into an XOY surface of a three-dimensional coordinate system, and the coordinate is A1(x1，y1)、A2(x2，y2)、A3(x3，y3)...An(xn，yn)；
(2) Respectively obtaining:
Figure FDA0002239236510000031
Figure FDA0002239236510000033
Figure FDA0002239236510000035
Figure FDA0002239236510000036
(3) obtaining a regression linear equation through the parameters in the step (2):
Figure FDA0002239236510000037
CN201910994193.XA 2019-10-18 2019-10-18 Variable offset distance advanced fault detection method based on shock wave excitation seismic source Active CN110850472B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910994193.XA CN110850472B (en) 2019-10-18 2019-10-18 Variable offset distance advanced fault detection method based on shock wave excitation seismic source

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910994193.XA CN110850472B (en) 2019-10-18 2019-10-18 Variable offset distance advanced fault detection method based on shock wave excitation seismic source

Publications (2)

Publication Number Publication Date
CN110850472A true CN110850472A (en) 2020-02-28
CN110850472B CN110850472B (en) 2021-07-02

Family

ID=69596549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910994193.XA Active CN110850472B (en) 2019-10-18 2019-10-18 Variable offset distance advanced fault detection method based on shock wave excitation seismic source

Country Status (1)

Country Link
CN (1) CN110850472B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363210A (en) * 2020-11-13 2021-02-12 福州华虹智能科技股份有限公司 Quantitative coal thickness prediction method based on joint inversion of wave velocity and attenuation coefficient of transmission channel waves
CN112904424A (en) * 2021-01-21 2021-06-04 中国科学院地理科学与资源研究所 Method and device for determining position of abnormal body and readable storage medium

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1992012444A1 (en) * 1991-01-02 1992-07-23 Western Atlas International, Inc. Method for determining geometry of subsurface features while drilling
KR100660563B1 (en) * 2005-09-30 2006-12-22 (주)지 텍크 Buoy-plate multi channel seismic survey system and method
CN102759746A (en) * 2011-04-28 2012-10-31 中国石油天然气集团公司 Method for inverting anisotropy parameters using variable offset vertical seismic profile data
CN103235334A (en) * 2013-04-11 2013-08-07 中国矿业大学 Low drainage roadway fault-throw advanced-detection method based on seismic signals
CN104166164A (en) * 2014-08-08 2014-11-26 山东科技大学 Three-component multiwave reflection three-dimensional earthquake advanced detection method for coal roadway tunneling geologic structure
JP2016133497A (en) * 2015-01-22 2016-07-25 株式会社セオコンプ System and method for underground structure survey
CN105974472A (en) * 2016-05-13 2016-09-28 中国矿业大学 Tunnel advanced detection speed modeling method based on reflected signal
CN108490488A (en) * 2018-03-19 2018-09-04 中煤科工集团西安研究院有限公司 It is a kind of based on when window energy ratio channel wave seismic scattering imaging method and device
CN109765606A (en) * 2019-01-28 2019-05-17 阳泉煤业（集团）有限责任公司 A kind of detection method of the hidden fault properties of stope based on reflection slot wave
CN110531417A (en) * 2019-08-21 2019-12-03 中国矿业大学 A kind of advanced multilayer speed fine modeling method based on polarization offset

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1992012444A1 (en) * 1991-01-02 1992-07-23 Western Atlas International, Inc. Method for determining geometry of subsurface features while drilling
KR100660563B1 (en) * 2005-09-30 2006-12-22 (주)지 텍크 Buoy-plate multi channel seismic survey system and method
CN102759746A (en) * 2011-04-28 2012-10-31 中国石油天然气集团公司 Method for inverting anisotropy parameters using variable offset vertical seismic profile data
CN103235334A (en) * 2013-04-11 2013-08-07 中国矿业大学 Low drainage roadway fault-throw advanced-detection method based on seismic signals
CN104166164A (en) * 2014-08-08 2014-11-26 山东科技大学 Three-component multiwave reflection three-dimensional earthquake advanced detection method for coal roadway tunneling geologic structure
JP2016133497A (en) * 2015-01-22 2016-07-25 株式会社セオコンプ System and method for underground structure survey
CN105974472A (en) * 2016-05-13 2016-09-28 中国矿业大学 Tunnel advanced detection speed modeling method based on reflected signal
CN108490488A (en) * 2018-03-19 2018-09-04 中煤科工集团西安研究院有限公司 It is a kind of based on when window energy ratio channel wave seismic scattering imaging method and device
CN109765606A (en) * 2019-01-28 2019-05-17 阳泉煤业（集团）有限责任公司 A kind of detection method of the hidden fault properties of stope based on reflection slot wave
CN110531417A (en) * 2019-08-21 2019-12-03 中国矿业大学 A kind of advanced multilayer speed fine modeling method based on polarization offset

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BO WANG，等: "Fine Imaging by Using Advanced Detection of Reflected Waves in Underground Coal Mine", 《EARTH SCIENCES RESEARCH JOURNAL》 *
胡泽安，等: "煤层槽波的极化特征及其滤波", 《合肥工业大学学报（自然科学版）》 *
金丹，等: "频率域极化滤波在槽波信号处理中的应用", 《煤炭学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363210A (en) * 2020-11-13 2021-02-12 福州华虹智能科技股份有限公司 Quantitative coal thickness prediction method based on joint inversion of wave velocity and attenuation coefficient of transmission channel waves
CN112363210B (en) * 2020-11-13 2023-10-13 福州华虹智能科技股份有限公司 Coal thickness quantitative prediction method based on transmission groove wave velocity and attenuation coefficient joint inversion
CN112904424A (en) * 2021-01-21 2021-06-04 中国科学院地理科学与资源研究所 Method and device for determining position of abnormal body and readable storage medium

Also Published As

Publication number Publication date
CN110850472B (en) 2021-07-02

Similar Documents

Publication Publication Date Title
Li et al. An overview of ahead geological prospecting in tunneling
CN106194159B (en) Mine inclination measurement while drilling exploration system and measuring method thereof
JP6982103B2 (en) Underground structure detection
CN101603423B (en) Direct current method for bed-parallel advanced detection of water bearing structure in mine tunnel
CN109613616B (en) Geological comprehensive three-dimensional detection method for shallow karst tunnel
CN109765606B (en) Method for detecting nature of hidden fault of stope face based on reflected trough wave
CN108957521B (en) Long-distance three-dimensional advanced geological prediction method for tunnel
CN102866417A (en) Device and method for seismic cross hole computed tomography (CT) detection and tomography of underground cave
CN105431612A (en) Drilling method and apparatus
CN110850472B (en) Variable offset distance advanced fault detection method based on shock wave excitation seismic source
CN105719433A (en) In-hole seismic wave based advanced prediction method
CN104880729A (en) Coal roadway anomalous structure advanced detection method based on continuous channel wave signal tracking
CN110632667B (en) Hidden collapse column advanced detection method based on shock wave shock condition
Luo et al. Tomographic imaging of rock conditions ahead of mining using the shearer as a seismic source—a feasibility study
CN208334667U (en) One kind is for soft rock tunnel Preceding geology predictor three-dimensional over long distances
CN102520441B (en) Method for acquiring earthquake exploration data through roadway oriented detection
CN104199110A (en) Method for channel wave earthquake stereoscopic advanced detection during underground coal mine support
CN117192615A (en) Method for detecting hidden geological structure in coal face based on transmission seismic wake wave
Zhou et al. Seeing coal‐seam top ahead of the drill bit through seismic‐while‐drilling
CN112130206A (en) Underground roadway full-space seismic detection system and detection method
CN110531413A (en) A kind of advanced Visualization Modeling method of craven fault
CN109991654A (en) A kind of Gas Outburst driving face gas pocket is with pick forward probe device and detection method
CN110531415A (en) A kind of three-dimensional craven fault forward probe method influenced using wall rock loosening ring
CN114839672A (en) Rapid advanced geological prediction method for tunnel capable of being measured while drilling in advanced horizontal drilling
Young et al. Seismic attenuation spectra in rock mass characterization; a case study in open-pit mining

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