CN110531418B - Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging - Google Patents
Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging Download PDFInfo
- Publication number
- CN110531418B CN110531418B CN201910774268.3A CN201910774268A CN110531418B CN 110531418 B CN110531418 B CN 110531418B CN 201910774268 A CN201910774268 A CN 201910774268A CN 110531418 B CN110531418 B CN 110531418B
- Authority
- CN
- China
- Prior art keywords
- breakpoint
- coal seam
- group
- seismic
- rays
- 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.)
- Active
Links
- 230000010287 polarization Effects 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000003384 imaging method Methods 0.000 title claims abstract description 11
- 239000003245 coal Substances 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 abstract description 9
- 238000005259 measurement Methods 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000013277 forecasting method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
Images
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. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging, wherein a three-component geophone receives a fed-back seismic wave signal to obtain a break-off wing coal seam breakpoint diffraction wave group of each three-component geophone; respectively obtaining main polarization direction rays of each three-component detector by adopting a polarization analysis method based on a Hilbert transform complex covariance matrix for the breakpoint diffraction wave group; grouping every two rays, so that the intersection or the distance of the two rays in the main polarization directions of each group is the minimum; and obtaining the position coordinates of the intersection points or the position coordinates of the middle points of each group; finally, counting to obtain the three-dimensional space position of the breaking wing coal seam breakpoint where the most intersection points and the most middle points exist in each three-dimensional grid; the method disperses the error generated by each detector during detection in the grouping and counting modes, and finally determines the three-dimensional coordinates of the breaking point of the loss-wing coal seam; thereby improving the accuracy of detecting the breaking point of the coal seam with the breaking wing.
Description
Technical Field
The invention relates to a breakpoint fine detection method, in particular to a breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging.
Background
The fault structure is a hidden disaster-causing factor of coal mine water inrush, coal and gas outburst and other serious disasters. The fault not only controls the distribution and migration of deep water, but also influences the occurrence state of coal beds and gas. At present, the mine seismic wave method is most effective in the advanced detection of fault parameters (including inclination, dip angle, fall and the like). The main concern of coal mine actual production is fault fall parameters, because the fall determines fault scale and safety precaution measures.
The mine seismic method for fault advanced detection is classified into a reflection trough wave method, a reflector wave method and the like. The reflection trough wave method can accurately judge the position of the coal seam fault of the plate where the coal roadway is located; the reflector wave method can realize accurate prediction of fault plane inclination and dip angle; however, both methods cannot determine the position of the coal seam with the breaking wing, and further cannot invert fault fall parameters according to the position of the coal seam with the breaking wing where the coal roadway is located and the position of the coal seam with the breaking wing.
To solve this problem, the patent numbers are: 201310123575.8, the invention provides a method for advanced detection of a loss wing coal seam in a coal roadway based on a single shot-blast inspection observation system, wherein the position of the loss wing coal seam is inverted by utilizing the diffracted wave of the breakpoint of the coal seam, the method utilizes X, Z two-component receiving point seismic records to calculate the main polarization direction of the diffracted wave of the breakpoint, a ray is formed by the two-dimensional coordinate of the receiving point and the main polarization direction, an arc is formed by taking the receiving point as the center of a circle and the distance of the diffracted wave of the breakpoint as the radius, and the intersection of the ray and the arc is the position of the loss wing coal seam. Similar to the above patent, 201610316770.6, the chinese invention patent proposes a qualitative forecasting method for projecting coal seam roadway small fault advance, which determines the position of the break point by intersecting the rays of two receiving points of the top plate or the bottom plate. However, the above methods have great limitations, mainly represented by: 1. the method has extremely high requirement on the accuracy of determining the main polarization direction, and once the obtained direction has a slight error, the position deviation of the broken-loss wing coal seam obtained by intersection is very large, however, the solving error of the main polarization direction is very large due to the complex actual seismic record under the coal mine; 2. the method is a two-dimensional exploration method, and the breakpoint of the actual loss wing coal seam is in a three-dimensional space, but the three-dimensional coordinate position of the actual loss wing coal seam cannot be determined by the method. In summary, the existing methods can not meet the requirement of fine detection of the breaking point of the coal seam with the breaking wing.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging, which can determine the three-dimensional coordinates of the breakpoint of a loss-wing coal seam; and the influence of measurement errors on detection can be effectively reduced, so that the accuracy of detecting the breaking point of the loss-of-wing coal seam is improved.
In order to achieve the purpose, the invention adopts the technical scheme that: a breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging comprises the following specific steps:
A. arranging an excitation point at the position 2m below a bottom plate and at the head of a coal roadway, wherein the excitation point is used for exciting seismic waves; arranging an even number of three-component detectors on one lateral side of the roadway behind the excitation point, wherein the even number of three-component detectors are arranged at equal intervals and are positioned on the same straight line with the excitation point, and the straight line is parallel to the direction of the roadway;
B. exciting seismic waves to the front of the head at a shock point, and then receiving the fed back seismic waves by each three-component geophone in real time and transmitting the received seismic waves to a seismic recorder;
C. the seismic recorder draws a seismic waveform diagram according to seismic waves received by each three-component detector, and then obtains a breaking wing coal seam breakpoint diffraction wave group of each three-component detector according to the seismic waveform diagram;
D. c, processing the break point diffraction wave group of the break-loss wing coal seam obtained in the step C by a complex covariance matrix polarization analysis method based on Hilbert transform, and respectively obtaining main polarization direction rays of each three-component detector;
E. determining the position of the breaking point of the breaking wing coal seam according to Hilbert polarization imaging, wherein the specific process comprises the following steps:
firstly, establishing a space coordinate system, and setting the space coordinate of any detector R as (x)R,yR,zR) And the direction information of the azimuth angle and the inclination angle of the diffracted wave received by the detector R is (m, n, p), and a parameter equation is listed, wherein t represents time:
calculating the intersection point position of any two main polarization direction rays:
the above equation can be expressed in the form of the following matrix:
VT＝P (3)
wherein:
if two intersected main polarization direction rays can be obtained due to the influence of noise interference and polarization method errors on the seismic records obtained by the seismic recorder, recording the position coordinates of the intersection point; if two intersecting main polarization direction rays cannot be obtained, the objective function solution J of the formula (4) is minimized, wherein J is expressed as:
J＝||T-V-1P||2 (5)
solving using least squares, i.e.
T＝(VTV)-1VTP
Thereby obtaining the minimum distance of two non-intersected main polarization direction rays, and recording the position coordinates of the midpoint of the minimum distance;
the method is adopted to group the even number of main polarization direction rays pairwise, so that the intersection or the distance of the two main polarization direction rays in each group is minimum; and obtaining the position coordinates of the intersection points or the position coordinates of the middle points of each group;
reflecting the position coordinates of the intersection points or the position coordinates of the middle points obtained by each group in the established three-dimensional space coordinate system, and dividing the three-dimensional space coordinate system into three-dimensional space grids by taking 1m as a unit;
and thirdly, counting the number of the intersection points and the middle points falling into each three-dimensional space grid, and determining the three-dimensional space grid with the largest number of the intersection points and the middle points as the three-dimensional space position of the fault.
Compared with the prior art, the invention is provided with an even number of three-component detectors, the three-component detectors receive the fed-back seismic wave signals respectively, and then the breaking wing coal seam breakpoint diffraction wave group of each three-component detector is obtained; adopting a polarization analysis method based on a Hilbert transform complex covariance matrix to the breakpoint diffraction wave group to respectively obtain main polarization direction rays and polarization characteristic parameters of each three-component detector; determining the position of a breakpoint according to Hilbert polarization imaging; grouping every two rays, so that the intersection or the distance of the two rays in the main polarization directions of each group is the minimum; and obtaining the position coordinates of the intersection points or the position coordinates of the middle points of each group; finally, dividing three-dimensional space grids, and counting to obtain the grids with the most intersection points and the most middle points in each three-dimensional grid, namely the three-dimensional space positions of the breaking wing coal seam breakpoints; the method comprises the steps of setting a plurality of three-component detectors, grouping every two detectors, dividing grid statistics, dispersing errors generated during detection of each detector in the grouping and statistics mode, and finally determining three-dimensional coordinates of break-off wing coal seam breakpoints; and the influence of measurement errors on detection can be effectively reduced, so that the accuracy of detecting the breaking point of the loss-of-wing coal seam is improved.
Drawings
FIG. 1 is a diagram of the placement of three-component detectors and excitation points according to the present invention;
fig. 2 is a schematic projection diagram of each intersection point and midpoint in the Y direction in a three-dimensional coordinate system according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of fig. 2 after being divided into three-dimensional space grids.
In the figure: 1. three-component detector, 2, roadway, 3, head-on, 4, excitation point, 5 and breaking wing coal seam breakpoint.
Detailed Description
The present invention will be further explained below.
As shown in the figure, the right side of the head in fig. 1 is taken as the front for patent description, and the specific steps of the invention are as follows:
A. an excitation point 4 is arranged at the position 2m below a bottom plate of a head 3 of the coal roadway and is used for exciting seismic waves; arranging 90 three-component detectors 1 on one side of the roadway 2 behind the excitation point, wherein the 90 three-component detectors 1 are arranged at equal intervals and are positioned on the same straight line with the excitation point 4, and the straight line is parallel to the direction of the roadway;
B. exciting seismic waves to the head-on front at a shock excitation point 4, and then receiving the fed back seismic waves by each three-component geophone 1 in real time and transmitting the seismic waves to a seismic recorder;
C. the seismic recorder draws a seismic oscillogram according to seismic waves received by each three-component geophone 1, and then obtains a diffracted wave group of a breaking wing coal seam breakpoint 5 of each three-component geophone 1 according to the seismic oscillogram;
D. c, processing the diffracted wave group of the break-off wing coal seam breakpoint 5 obtained in the step C by using a complex covariance matrix polarization analysis method based on Hilbert transform, and respectively obtaining main polarization direction rays of each three-component detector 1;
E. determining the position of a breaking wing coal seam breakpoint 5 according to Hilbert polarization imaging, wherein the specific process comprises the following steps:
firstly, establishing a space coordinate system, and setting the space coordinate of any detector R as (x)R,yR,zR) And the direction information of the azimuth angle and the inclination angle of the diffracted wave received by the detector R is (m, n, p), and a parameter equation is listed, wherein t represents time:
calculating the intersection point position of any two main polarization direction rays:
the above equation can be expressed in the form of the following matrix:
VT＝P (3)
wherein:
if two intersected main polarization direction rays can be obtained due to the influence of noise interference and polarization method errors on the seismic records obtained by the seismic recorder, recording the position coordinates of the intersection point; if two intersecting main polarization direction rays cannot be obtained, the objective function solution J of the formula (4) is minimized, wherein J is expressed as:
J＝||T-V-1P||2 (5)
solving using least squares, i.e.
T＝(VTV)-1VTP
Thereby obtaining the minimum distance of two non-intersected main polarization direction rays, and recording the position coordinates of the midpoint of the minimum distance;
the method is adopted to group every two of the 90 main polarization direction rays, so that the intersection or the distance of the two main polarization direction rays in each group is minimum; and obtaining 45 position coordinates of the intersection points or the midpoint of each group;
reflecting the position coordinates of the 45 intersection points or the middle points obtained by each group in the established three-dimensional space coordinate system, and dividing the three-dimensional space coordinate system into three-dimensional space grids by taking 1m as a unit;
and thirdly, counting the number of the intersection points and the middle points falling into each three-dimensional space grid, and determining the three-dimensional space grid with the largest number of the intersection points and the middle points as the three-dimensional space position of the breaking wing coal seam breakpoint 5.
Claims (1)
1. A breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging is characterized by comprising the following specific steps:
A. arranging an excitation point at the position 2m below a bottom plate and at the head of a coal roadway, wherein the excitation point is used for exciting seismic waves; arranging an even number of three-component detectors on one lateral side of the roadway behind the excitation point, wherein the even number of three-component detectors are arranged at equal intervals and are positioned on the same straight line with the excitation point, and the straight line is parallel to the direction of the roadway;
B. exciting seismic waves to the head-on front at an excitation point, and then receiving the fed back seismic waves by each three-component geophone in real time and transmitting the received seismic waves to a seismic recorder;
C. the seismic recorder draws a seismic waveform diagram according to seismic waves received by the three-component detectors, and then obtains a breaking wing coal seam breakpoint diffraction wave group of each three-component detector according to the seismic waveform diagram;
D. c, processing the break point diffraction wave group of the break-loss wing coal seam obtained in the step C by adopting a Hilbert transform complex covariance matrix polarization analysis method to respectively obtain main polarization direction rays of each three-component detector;
E. determining the position of the breaking point of the breaking wing coal seam according to Hilbert polarization imaging, wherein the specific process comprises the following steps:
firstly, establishing a space coordinate system, and setting the space coordinate of any three-component detector R as (x)R,yR,zR) The azimuth and inclination direction information of the diffracted wave received by the three-component detector R is (m, n, p), and a parametric equation is listed, wherein t represents time:
calculating the intersection point position of any two main polarization direction rays:
the above equation can be expressed in the form of the following matrix:
VT＝P (3)
wherein:
if two intersected main polarization direction rays can be obtained due to the influence of noise interference and polarization method errors on the seismic records obtained by the seismic recorder, recording the position coordinates of the intersection points; if two intersecting main polarization direction rays cannot be obtained, the objective function solution J of the formula (4) is minimized, wherein J is expressed as:
J＝||T-V-1P||2 (5)
solving using least squares, i.e.
T＝(VTV)-1VTP
Thereby obtaining the minimum distance of two non-intersected main polarization direction rays, and recording the position coordinates of the midpoint of the minimum distance;
the method is adopted to group the even number of main polarization direction rays pairwise, so that the intersection or the distance of the two main polarization direction rays in each group is minimum; and obtaining the position coordinates of the intersection points or the position coordinates of the middle points of each group;
reflecting the position coordinates of the intersection points or the position coordinates of the middle points obtained by each group in the established three-dimensional space coordinate system, and dividing the three-dimensional space coordinate system into three-dimensional space grids by taking 1m as a unit;
and thirdly, counting the number of the intersection points and the middle points falling into each three-dimensional space grid, and determining the three-dimensional space grid with the largest number of the intersection points and the middle points as the three-dimensional space position of the breaking wing coal seam breakpoint.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910774268.3A CN110531418B (en) | 2019-08-21 | 2019-08-21 | Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910774268.3A CN110531418B (en) | 2019-08-21 | 2019-08-21 | Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110531418A CN110531418A (en) | 2019-12-03 |
| CN110531418B true CN110531418B (en) | 2020-11-20 |
Family
ID=68662503
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910774268.3A Active CN110531418B (en) | 2019-08-21 | 2019-08-21 | Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110531418B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113030691B (en) * | 2019-12-24 | 2022-07-19 | 芯恩（青岛）集成电路有限公司 | Chip electrical test probe head alignment method, system, storage medium and terminal |
| CN111965695B (en) * | 2020-08-17 | 2023-04-25 | 山西潞安环保能源开发股份有限公司五阳煤矿 | Small fault fall detection method based on reflection groove wave |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7463552B1 (en) * | 2003-11-25 | 2008-12-09 | Michael John Padgett | Method for deriving 3D output volumes using filters derived from flat spot direction vectors |
| CN101604029A (en) * | 2009-05-06 | 2009-12-16 | 孔庆丰 | Polarization treatment method of track type controllable earthquake source |
| CN101915939A (en) * | 2010-08-17 | 2010-12-15 | 中国科学院地质与地球物理研究所 | Surface wave pressing method |
| CN103217719A (en) * | 2013-04-11 | 2013-07-24 | 中国矿业大学 | Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system |
| CN105974472A (en) * | 2016-05-13 | 2016-09-28 | 中国矿业大学 | Tunnel advanced detection speed modeling method based on reflected signal |
| CN106019371A (en) * | 2016-05-13 | 2016-10-12 | 中国矿业大学 | Outburst coal seam roadway minor fault advanced qualitative forecast method |
| CN106154327A (en) * | 2016-08-17 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of method improving hidden fault recognizing precision |
| CN106443765A (en) * | 2016-08-30 | 2017-02-22 | 安徽惠洲地质安全研究院股份有限公司 | Urban engineering seismic detection comprehensive imaging method based on multi-component observation system |
-
2019
- 2019-08-21 CN CN201910774268.3A patent/CN110531418B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7463552B1 (en) * | 2003-11-25 | 2008-12-09 | Michael John Padgett | Method for deriving 3D output volumes using filters derived from flat spot direction vectors |
| CN101604029A (en) * | 2009-05-06 | 2009-12-16 | 孔庆丰 | Polarization treatment method of track type controllable earthquake source |
| CN101915939A (en) * | 2010-08-17 | 2010-12-15 | 中国科学院地质与地球物理研究所 | Surface wave pressing method |
| CN103217719A (en) * | 2013-04-11 | 2013-07-24 | 中国矿业大学 | Method of advanced detection of breaking-loss wing coal seam of coal road based on single offset pair observation system |
| CN105974472A (en) * | 2016-05-13 | 2016-09-28 | 中国矿业大学 | Tunnel advanced detection speed modeling method based on reflected signal |
| CN106019371A (en) * | 2016-05-13 | 2016-10-12 | 中国矿业大学 | Outburst coal seam roadway minor fault advanced qualitative forecast method |
| CN106154327A (en) * | 2016-08-17 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of method improving hidden fault recognizing precision |
| CN106443765A (en) * | 2016-08-30 | 2017-02-22 | 安徽惠洲地质安全研究院股份有限公司 | Urban engineering seismic detection comprehensive imaging method based on multi-component observation system |
Non-Patent Citations (1)
| Title |
|---|
| 反射槽波绕射偏移成像及应用;姬广忠 等;《煤田地质与勘探》;20170228;第45卷(第1期);第121-124页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110531418A (en) | 2019-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6553315B2 (en) | Seismic imaging using omni-azimuth seismic energy sources and directional sensing | |
| US3581273A (en) | Marine seismic exploration | |
| US6697737B2 (en) | Quality control cube for seismic data | |
| Liu et al. | A new 3D observation system designed for a seismic ahead prospecting method in tunneling | |
| CN110531418B (en) | Breakpoint three-dimensional fine positioning method based on Hilbert polarization imaging | |
| US6470276B1 (en) | Method of determining merit of offset and azimuth distributions in seismic data acquisition | |
| CN112578428B (en) | Scattering multi-wave advanced detection method based on roadway vertical virtual survey line | |
| CN109765606B (en) | Method for detecting nature of hidden fault of stope face based on reflected trough wave | |
| US5189643A (en) | Method of accurate fault location using common reflection point gathers | |
| CN106443765A (en) | Urban engineering seismic detection comprehensive imaging method based on multi-component observation system | |
| US6529832B1 (en) | Seismic imaging using omni-azimuth seismic energy sources and directional sensing | |
| Kurz et al. | Source localization | |
| US4813029A (en) | Geophone apparatus and a seismic exploration method | |
| CN110632667B (en) | Hidden collapse column advanced detection method based on shock wave shock condition | |
| US5307268A (en) | Method for using P-S converted waves as a subsalt imaging tool | |
| CN111103620A (en) | Three-dimensional offset imaging method for rock roadway advanced detection | |
| CN110531413B (en) | Advanced three-dimensional visual modeling method for small fault | |
| US6678617B2 (en) | Mapping subsurface open fractures using elastically nonlinear interaction of two surface-generated waves | |
| CN110531416B (en) | Fault determination method based on time-frequency domain polarization parameters of three-component reflection signals | |
| US7663971B2 (en) | Resonance scattering seismic method | |
| EP0297852A2 (en) | Method for real time display of marine seismic survey data coverage | |
| Alessandrini et al. | A simple P-wave polarization analysis: its application to earthquake location | |
| Dunkin et al. | Isochrons for a three-dimensional seismic system | |
| US20140297191A1 (en) | Predicting sensitivity to positioning for seismic surveys | |
| Spencer et al. | Vertical seismic profile polarization method for determining reflector orientation |
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 |