CN105738942A - Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process - Google Patents
Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process Download PDFInfo
- Publication number
- CN105738942A CN105738942A CN201610134734.8A CN201610134734A CN105738942A CN 105738942 A CN105738942 A CN 105738942A CN 201610134734 A CN201610134734 A CN 201610134734A CN 105738942 A CN105738942 A CN 105738942A
- Authority
- CN
- China
- Prior art keywords
- whilst
- tour
- slowness
- mine rock
- data
- 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.)
- Pending
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
Abstract
The invention discloses a mine rock structure evolution three-dimension continuous exploration system and a method for exploitation process. The method comprises steps of arranging a detector and an impact seismic source on two sides of the rock, transmitting the detection data of the detector and the impact seismic source to a data collection unit, and transmitting the data to the data center and the processing system for processing. The mine rock structure evolution three-dimension continuous exploration system and method for exploitation process can perform automatic detection in the three-dimension monitoring space, which is continuous in time, time saving and labor saving, and small in error. The mine rock structure evolution three-dimension continuous exploration system and method for exploitation process solvethe problem in the prior art the current detection method is incontinuity in the time, and is time-labor wasting by the staff operation.
Description
Technical field
The field of the invention is ore and rock stability monitoring and the analysis field of metal mine.
Background technology
After metal mine enters deep mining, ore deposit rock stress raises, the bad stability of engineering, it is easy to the accident such as wall caving, rock burst occurs.Rock stability is by rock and structures shape, and along with the structure carrying out rock mass of exploitation is also ceaselessly changing, therefore understands fully the change optimization important role at the advanced prediction of accident, mining Design of structure.Have studied a lot of detection methods and means for these people, as borehole camera, ultrasonic acoustic wave detect, but these method of testings are discontinuous in time, by human users, waste time and energy.
Summary of the invention
In order to solve above-mentioned technical problem, the invention provides a kind of three-dimensional continuous probe System and method for of recovery process Structure of mine rock evolution, cymoscope and impact origin it is provided with in rock mass both sides, the detection data of cymoscope and impact origin are transferred in data acquisition unit, and be sent in data to data center and process system and process, solve the method for testing existed in prior art discontinuous in time, by human users, the technical problem wasted time and energy.
To achieve these goals, the technical solution used in the present invention is: the three-dimensional continuous probe method of recovery process Structure of mine rock evolution, the steps include:
1), the three-dimensional continuous probe system of recovery process Structure of mine rock evolution is installed on rock mass；
2), send time signal by time service unit timing, make cymoscope keep Tong Bu with the clock of impact origin；
3), by the data message of impact origin and the Real-time Collection waveform transfer of cymoscope in data acquisition unit, and upload to data be sent to data center with in process system；
4), data process:
4a), determine each impact origin of system and the position coordinates of each cymoscope, select grid cell, search coverage is carried out stress and strain model, forms equivalently-sized grid cell；
4b), ray path is processed as straight line, total duration of ray and the ratio calculated of distance as initial slowness, initial slowness is assigned to each grid cell, as initial model；By calculating whilst on tour, carry out ray tracing, try to achieve ray length in each grid；
4c), initial slowness in the ray length obtained and grid is adopted to represent the whilst on tour of known point, by actual observation to the whilst on tour of whilst on tour and calculating do difference, draw error whilst on tour, thus being finally inversed by slowness error, thus obtaining the correction value of slowness in grid；
4d), complete first time iteration, then by based on the correction value of slowness, iteration, when calculate the error of whilst on tour and observation whilst on tour within the scope of restriction time, export the slowness of this iteration；
5), impact origin regularly produce the seismic wave of different frequency, and impact origin moves back and forth within the specific limits；Repeat step 3) to step 4), it is achieved the detection automatically continuously of rock mass structure, and obtain the changing value of Structure of mine rock in detection range by the slowness of iteration, utilize interpolation algorithm drawing three-dimensional cloud atlas, the evolutionary process of Structure of mine rock in acquisition region.
Step 1) described in recovery process Structure of mine rock develop three-dimensional continuous probe system particularly as follows: be provided with the impact origin be more than or equal on the rock mass side crag needing detection, and accessed time service unit and data acquisition unit by coaxial cable or optical fiber；The rock mass opposite side crag needing detection is installed the cymoscope corresponding with each impact origin, and is accessed time service unit and data acquisition unit by coaxial cable or optical fiber；Time service unit and data acquisition unit connect data center and the system of process.
Described impact origin is from portable impact origin.
Described step 4a) described in select NE particularly as follows: medium to be divided into length and width and the high square grid cell being vibroseis generator average distance 1.5 times.
Described step 4b) described in calculate whilst on tour particularly as follows: ripple from focus in media as well along specifying paths to arrive the time used by specified point.
Having the beneficial effects that of the invention: the invention provides a kind of three-dimensional continuous probe System and method for of recovery process Structure of mine rock evolution, cymoscope and impact origin it is provided with in rock mass both sides, the detection data of cymoscope and impact origin are transferred in data acquisition unit, and are sent in data to data center and process system and process.The present invention passes through said system and method, continuously automatically detection in three-dimension monitor space, Time Continuous, time saving and energy saving and Effect on Detecting good, and error is little.
Accompanying drawing explanation
Fig. 1: create system for the present invention.
Fig. 2: for flow chart of data processing figure of the present invention.
Fig. 3: for ray tracing in the inventive method to pretreatment process.
Fig. 4: for linear interpolation geometrical relationship schematic diagram in the inventive method.
Detailed description of the invention
The three-dimensional continuous probe method of recovery process Structure of mine rock evolution, the steps include:
1), the three-dimensional continuous probe system 1 of recovery process Structure of mine rock evolution is installed on rock mass, described recovery process Structure of mine rock develops three-dimensional continuous probe system 1 particularly as follows: be provided with the impact origin 4 be more than or equal on the rock mass side crag needing detection, and is accessed time service unit 2 and data acquisition unit 5 by coaxial cable or optical fiber；The rock mass opposite side crag needing detection is installed the cymoscope 3 corresponding with each impact origin 4, and is accessed time service unit 2 and data acquisition unit 5 by coaxial cable or optical fiber；Time service unit 2 and data acquisition unit 5 connect data center and process system 6.Impact origin 4 is from portable impact origin.
2), send time signal by time service unit 2 timing, make cymoscope 3 keep Tong Bu with the clock of impact origin 4；
3), by the data message of impact origin 4 and the Real-time Collection waveform transfer of cymoscope 3 in data acquisition unit 5, and upload to data be sent to data center with in process system 6；
4), data process:
4a), each impact origin 4 of system and the position coordinates of each cymoscope 3 are determined, select grid cell (medium is divided into length and width and the high square grid cell being vibroseis generator average distance 1.5 times), search coverage is carried out stress and strain model, forms equivalently-sized grid cell；
4b), ray path is processed as straight line, total duration of ray and the ratio calculated of distance as initial slowness, initial slowness is assigned to each grid cell, as initial model；By calculating whilst on tour, whilst on tour is that ripple arrives the time used by specified point along appointment paths in media as well from focus, carries out ray tracing, tries to achieve ray length in each grid；
4c), initial slowness in the ray length obtained and grid is adopted to represent the whilst on tour of known point, by actual observation to the whilst on tour of whilst on tour and calculating do difference, draw error whilst on tour, thus being finally inversed by slowness error, thus obtaining the correction value of slowness in grid；
4d), complete first time iteration, then by based on the correction value of slowness, iteration, when calculate the error of whilst on tour and observation whilst on tour within the scope of restriction time, export the slowness of this iteration；
5, impact origin 4 regularly produces the seismic wave of different frequency, and impact origin 4 moves back and forth within the specific limits；Repeat step 3 and arrive step 4, it is achieved the detection automatically continuously of rock mass structure, and obtain the changing value of Structure of mine rock in detection range by the slowness of iteration, utilize interpolation algorithm drawing three-dimensional cloud atlas, obtain the evolutionary process of Structure of mine rock in region.
Step during concrete application is:
1, a number of certainly portable impact origin 4 is installed on the rock mass side crag needing detection, and is accessed time service unit 2 and data acquisition unit 5 by coaxial cable or optical fiber.The rock mass opposite side crag needing detection is installed a number of earthquake wave detector 3, and is accessed time service unit 2 and data acquisition unit 5 by coaxial cable or optical fiber.Time service unit 2 and data acquisition unit 5 are connected data center and process system 6 by network system.
2, time service unit 2 timing sends time signal, makes cymoscope 3 keep Tong Bu with the clock from portable impact origin 4.
3 impact origins produce seismic wave, and its frequency and position, temporal information are transferred in data acquisition unit 5.Cymoscope 3 Real-time Collection Wave data, is sent in data acquisition unit 5.5 data collected of data acquisition unit are sent in data center and process system 6.
4 data handling procedures and algorithm: determine each focus of observation system, each detector position coordinate；Select suitable grid cell (suitable grid refers to medium is divided into length and width and the high square grid cell being vibroseis generator average distance 1.5 times), search coverage is carried out stress and strain model, forms the grid cell that yardstick one by one is identical.
Ray path is processed into straight line, and total duration of ray and the ratio of total distance are as initial slowness.This initial slowness is assigned to each grid cell, as initial model.Then pass through calculating whilst on tour, carry out ray tracing, try to achieve the length of ray in each grid.The whilst on tour of known point is represented by slowness initial in the length obtained and grid, the whilst on tour of the whilst on tour of actual observation Yu calculating is done difference, draw error whilst on tour, then being finally inversed by slowness error, thus obtaining the correction value of slowness in grid, now first time iteration completes, again based on the correction value of slowness, iteration, when the error calculating whilst on tour and observe whilst on tour is when being satisfied with in scope, exports the slowness of this iteration.Wherein,
4-1 three-dimensional ray tracing (as shown in Figure 3) step is:
4-1-1 is to pre-treatment
(1) each whilst on tour calculating point on the elementary boundary of shot point (focus) place is calculated, in Fig. 3 shown in a；
(2), on coordinate W direction, each whilst on tour calculating point on shot point unit column each unit border is calculated, in Fig. 3 shown in b；
(3), on coordinate U direction, each whilst on tour calculating point on each row each unit border, shot point unit rear is calculated from front to back line by line；
(4), on coordinate U direction, each whilst on tour calculating point on each row each unit border, shot point unit front is recalculated from back to front line by line；
(5) on coordinate W direction, by the left-to-right whilst on tour recalculating by column and respectively calculating point on the right side of shot point unit on each row each unit border, and its minima is taken with the contrast of former result of calculation；
(6) on coordinate W direction, the right side a to left side recalculate by column and on the left of shot point unit, rank each whilst on tour calculating point on each unit border, and take its minima with the contrast of former result of calculation；
(5), on coordinate V direction, the whilst on tour respectively calculating point above the layer of shot point unit place on each layer each unit border is successively calculated from down to up；
(6), on coordinate V direction, the whilst on tour respectively calculating point below the layer of shot point unit place on each layer each unit border is successively calculated from top to bottom；
(7) on coordinate W direction, by the left-to-right whilst on tour successively calculating and respectively calculating point on the right side of the layer of shot point unit place on each row each unit border, and its minima is taken with the contrast of former result of calculation；
(8) on coordinate W direction, the right side a to left side successively calculate the whilst on tour respectively calculating point on the left of the layer of shot point unit place on each row each unit border, and take its minima with the contrast of former result of calculation；
(7), on coordinate U direction, successively calculate the whilst on tour respectively calculating point on rear side of the layer of shot point unit place on each row each unit border from front to back, and take its minima with the contrast of former result of calculation；
(8), on coordinate U direction, successively calculate the whilst on tour respectively calculating point on front side of the layer of shot point unit place on each row each unit border from back to front, and take its minima with the contrast of former result of calculation；
4-1-2 is to post processing
According to whilst on tour result of calculation calculated on point each on each unit border, choosing the ray path with minimum traveltimes, processing procedure is as follows:
(1) intersecting point coordinate of minimum traveltimes and ray and the grid surface receiving point is calculated；
(2) intersecting point coordinate tried to achieve with step (1) is for new reception point, forwards step (1) to, till shot point unit；
(3) ray intersection tried to achieve is coupled together with shot point, reception point, namely complete the tracking of a ray；
(4) receive dot information according to other, obtain other ray whilst on tours by (1), (2), (3) step.
Calculating minimum traveltimes in step 4-1-1 and step 4-1-2 (1), as shown in Figure 4, algorithm is its linear interpolation geometrical relationship:
Three-dimensional speed media model is separated into uniform cube grid cell, in each grid cell, is set to constant velocity.Linear distribution when this algorithm is assumed to walk on grid cell boundary face, in boundary face, linear interpolation when being walked by adjacent 4 grid nodes (summit) when walking of arbitrfary point is tried to achieve.Ray path within each cubic units is straight line, is illustrated in figure 4 a cubic units, A, B, C, D be known walk time node, be intended to ask arrive through A, B, C, D face F point minimum walk time and the position of intersection point E.
Rectangular parallelepiped grid the unit step-length in u, v, w direction respectively du, dv, dw, it is grid cell initial point with A (uA, vA, wA) point, the coordinate of some F is (uF, vF, wF), the slowness of grid cell is s (inverse of speed), difference TA, TB, TC, TD point E relative coordinate in grid when walking of known A, B, C, D is (u, v, 0), obtain:
When F node is walked it is:
Available steepest descent method solves above formula, obtains the value of u, v, TF.
4-2 inversion algorithm process:
Bidiagonalization matrix A, namely makes with U and V orthogonal matrix:
A=UBVT
Wherein
(1) initialize
β1u1=b, α1v1=ATu1,w1=v1,x0=0, d1=0
u1, b is m dimensional vector, d1,w1,v1,x0It is n-dimensional vector, β1,α1,γ1,ρ1It it is real constant.
(2) QR decomposes
vi+1=λ di+1-βi+1vi,αi+1vi+1=ATui+1-vi+1
(3) amendment parameter
si=βi+1/ρi,θi+1=siαi+1
(4) iterative
wi+1=vi+1-(θi+1/ρi)wi
If cond (A) >=CONLIM, then terminate iteration.
5 impact origins regularly produce the seismic wave of different frequency, and impact origin moves back and forth within the specific limits.Repeat step 3 and step 4, the detection automatically continuously of rock mass structure can be realized, calculate the changing value of Structure of mine rock in monitoring range, utilize interpolation algorithm drawing three-dimensional cloud atlas, it is thus achieved that the evolutionary process of region surrounding rock structure.
Claims (5)
1. the three-dimensional continuous probe method of recovery process Structure of mine rock evolution, the steps include:
1), the three-dimensional continuous probe system (1) of recovery process Structure of mine rock evolution is installed on rock mass；
2), send time signal by time service unit (2) timing, make cymoscope (3) keep Tong Bu with the clock of impact origin (4)；
3), by the data message of impact origin (4) and the Real-time Collection waveform transfer of cymoscope (3) in data acquisition unit (5), and upload to data be sent to data center with in process system (6)；
4), data process:
4a), determine each impact origin (4) of system and the position coordinates of each cymoscope (3), select grid cell, search coverage is carried out stress and strain model, forms equivalently-sized grid cell；
4b), ray path is processed as straight line, total duration of ray and the ratio calculated of distance as initial slowness, initial slowness is assigned to each grid cell, as initial model；By calculating whilst on tour, carry out ray tracing, try to achieve ray length in each grid；
4c), adopt initial slowness in the ray length that obtains and grid to represent the whilst on tour of known point, by actual observation to the whilst on tour of whilst on tour and calculating do difference, draw error whilst on tour, thus being finally inversed by slowness error, thus obtaining the correction value of slowness in grid；
4d), complete first time iteration, then by based on the correction value of slowness, iteration, when calculate the error of whilst on tour and observation whilst on tour within the scope of restriction time, export the slowness of this iteration；
5), impact origin (4) regularly produce the seismic wave of different frequency, and impact origin (4) moves back and forth within the specific limits；Repeat step 3) to step 4), it is achieved the detection automatically continuously of rock mass structure, and obtain the changing value of Structure of mine rock in detection range by the slowness of iteration, utilize interpolation algorithm drawing three-dimensional cloud atlas, the evolutionary process of Structure of mine rock in acquisition region.
2. the three-dimensional continuous probe method of recovery process Structure of mine rock according to claim 1 evolution, it is characterized in that: the recovery process Structure of mine rock described in step 1) develops three-dimensional continuous probe system (1) particularly as follows: be provided with the impact origin (4) be more than or equal on the rock mass side crag needing detection, and is accessed time service unit (2) and data acquisition unit (5) by coaxial cable or optical fiber；The rock mass opposite side crag needing detection is installed the cymoscope (3) corresponding with each impact origin (4), and is accessed time service unit (2) and data acquisition unit (5) by coaxial cable or optical fiber；Time service unit (2) and data acquisition unit (5) connect data center and the system of process (6).
3. the three-dimensional continuous probe method of recovery process Structure of mine rock according to claim 2 evolution, it is characterised in that: described impact origin (4) is from portable impact origin.
4. the three-dimensional continuous probe method of recovery process Structure of mine rock according to claim 1 evolution, it is characterised in that: the selection NE described in described step 4a) particularly as follows: be divided into length and width and the high square grid cell being vibroseis generator average distance 1.5 times by medium.
5. the three-dimensional continuous probe method of recovery process Structure of mine rock according to claim 1 evolution, it is characterised in that: the calculating whilst on tour described in described step 4b) is particularly as follows: ripple specifies the time used by paths arrival specified point from focus edge in media as well.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610134734.8A CN105738942A (en) | 2016-03-09 | 2016-03-09 | Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610134734.8A CN105738942A (en) | 2016-03-09 | 2016-03-09 | Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105738942A true CN105738942A (en) | 2016-07-06 |
Family
ID=56250219
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610134734.8A Pending CN105738942A (en) | 2016-03-09 | 2016-03-09 | Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105738942A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267813A (en) * | 2021-06-18 | 2021-08-17 | 重庆大学 | Velocity structure and seismic source position joint inversion method, system, terminal and readable storage medium |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5978313A (en) * | 1997-09-30 | 1999-11-02 | Trimble Navigaiton Limited | Time synchronization for seismic exploration system |
| CN101706585A (en) * | 2009-10-28 | 2010-05-12 | 中南大学 | Method for electrically forecasting danger in underground tunneling engineering |
| CN201508417U (en) * | 2009-09-22 | 2010-06-16 | 山东大学 | Tunnel or pit advance geological prediction induced polarization time-domain receiver equipment |
| CN102520441A (en) * | 2011-11-30 | 2012-06-27 | 朱德兵 | Method for acquiring earthquake exploration data through roadway oriented detection |
| CN102879468A (en) * | 2012-09-18 | 2013-01-16 | 东北大学 | Double-bending-element ultrasonic sensing test device and method for evaluating rock damage |
| CN103122764A (en) * | 2013-02-05 | 2013-05-29 | 中国矿业大学 | Automatic mobile type transient electromagnetic field detection method and device |
| CN103389523A (en) * | 2013-08-01 | 2013-11-13 | 中国建筑第四工程局有限公司 | Method and system for forecasting tunnel geology |
| CN103698810A (en) * | 2013-12-10 | 2014-04-02 | 山东科技大学 | Hybrid network minimum travel time ray tracing tomography method |
| CN104678428A (en) * | 2015-03-11 | 2015-06-03 | 山东大学 | Tunnel boring machine rock breaking seismic source and active source three-dimensional seismic combined advanced detection system |
-
2016
- 2016-03-09 CN CN201610134734.8A patent/CN105738942A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5978313A (en) * | 1997-09-30 | 1999-11-02 | Trimble Navigaiton Limited | Time synchronization for seismic exploration system |
| CN201508417U (en) * | 2009-09-22 | 2010-06-16 | 山东大学 | Tunnel or pit advance geological prediction induced polarization time-domain receiver equipment |
| CN101706585A (en) * | 2009-10-28 | 2010-05-12 | 中南大学 | Method for electrically forecasting danger in underground tunneling engineering |
| CN102520441A (en) * | 2011-11-30 | 2012-06-27 | 朱德兵 | Method for acquiring earthquake exploration data through roadway oriented detection |
| CN102879468A (en) * | 2012-09-18 | 2013-01-16 | 东北大学 | Double-bending-element ultrasonic sensing test device and method for evaluating rock damage |
| CN103122764A (en) * | 2013-02-05 | 2013-05-29 | 中国矿业大学 | Automatic mobile type transient electromagnetic field detection method and device |
| CN103389523A (en) * | 2013-08-01 | 2013-11-13 | 中国建筑第四工程局有限公司 | Method and system for forecasting tunnel geology |
| CN103698810A (en) * | 2013-12-10 | 2014-04-02 | 山东科技大学 | Hybrid network minimum travel time ray tracing tomography method |
| CN104678428A (en) * | 2015-03-11 | 2015-06-03 | 山东大学 | Tunnel boring machine rock breaking seismic source and active source three-dimensional seismic combined advanced detection system |
Non-Patent Citations (2)
| Title |
|---|
| 段成龙等: "跨孔地震CT技术在地铁工程施工溶洞探测方面的应用", 《地质论评》 * |
| 物探监督编写组: "《物探监督》", 31 October 2012 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267813A (en) * | 2021-06-18 | 2021-08-17 | 重庆大学 | Velocity structure and seismic source position joint inversion method, system, terminal and readable storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104133245B (en) | The static correcting method and system of a kind of seismic data | |
| CN108064348A (en) | Chromatography conversion method when a kind of earthquake based on two spots ray tracing is walked | |
| CN103105624B (en) | Longitudinal and transversal wave time difference positioning method based on base data technology | |
| US11042676B2 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN105954796B (en) | A kind of method and apparatus for the hypocentral location for determining microseism | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| CN104166161A (en) | Method and device for predicating fractures based on elliptical velocity inversion of anisotropism | |
| CN106814391A (en) | Ground micro-seismic state event location method based on Fresnel zone tomographic inversion | |
| CN109459787B (en) | coal mine underground structure imaging method and system based on seismic channel wave full-waveform inversion | |
| CN105445789A (en) | Three-dimensional Fresnel volume travel-time tomographic method based on multiple reflected refraction wave constraint | |
| CN104181593B (en) | Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device | |
| CN113484910B (en) | Tunnel advanced geological prediction method and system based on seismic interferometry | |
| CN106199704B (en) | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method | |
| CN107817516A (en) | Near surface modeling method and system based on preliminary wave information | |
| CN105607119B (en) | Near-surface model construction method and static correction value acquiring method | |
| CN103364823B (en) | Vibroseis real-time positioning and the system of analysis | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| CN110515122B (en) | Forward grid search positioning and micro-seismic signal identification method and device | |
| CN109655890A (en) | A kind of shallow mid-deep strata joint chromatography inversion speed modeling method of Depth Domain and system | |
| CN103698810A (en) | Hybrid network minimum travel time ray tracing tomography method | |
| CN105487106B (en) | A kind of benefit big gun method based on the illumination of Gaussian ray bundle target zone energy | |
| CN107870361B (en) | A kind of earthquake diving Wave chromatography imaging method, device and terminal device | |
| CN105738942A (en) | Mine rock structure evolution three-dimension continuous exploration system and method for exploitation process | |
| CN107817524A (en) | The method and apparatus of three-dimensional seismic tomography |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160706 |