CN104965221A - Mine advanced detection earthquake data superposing processing method - Google Patents
Mine advanced detection earthquake data superposing processing method Download PDFInfo
- Publication number
- CN104965221A CN104965221A CN201510124046.9A CN201510124046A CN104965221A CN 104965221 A CN104965221 A CN 104965221A CN 201510124046 A CN201510124046 A CN 201510124046A CN 104965221 A CN104965221 A CN 104965221A
- Authority
- CN
- China
- Prior art keywords
- data
- wave field
- apparent velocity
- negative apparent
- earthquake
- 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
Abstract
The invention discloses a mine advanced detection earthquake data superposing processing method, for solving the problems of improved signal-to-noise ratio of reflection waves and low imaging quality of multi-purpose layers by use of a negative apparent velocity method. The method is characterized by comprising the following steps: TSP superposing processing of a sampling rate, i.e., performing time-invariable processing transformation on earthquake reflection data (digital earthquake tape records) acquired by a gallery or a tunnel to obtain a time sampling rate as 1ms or 2ms; separation of an earthquake data positive-negative apparent velocity wave field, i.e., performing separation processing on the wave field comprising positive and negative apparent velocities to obtain a seismic wave field with a negative apparent velocity; and dynamic correction and superposition, i.e., performing dynamic correction on the reflection profile of the negative apparent velocity, performing gather drawing on dynamically corrected data and then carrying out the superposition. The method is wide in application scope, is not restricted to complex tomography and large inclination angle and tendency changes and can improve the accuracy in evaluating construction exploration in front of the gallery.
Description
Technical field
The present invention relates to the digital processing method of the forward probe seismic data of artificial digging laneway in a kind of mine.The method is to TSP seismic data through high precision wave field separation, and superposition is converted into superposition time data body A
stk(x, t), to A
stk(x, t) makes an explanation, and realizes the method for mine forward probe seismic data overlap-add procedure, belongs to two-dimensional seismic survey data processing method.
Background technology
For a long time, in the explanation of TSP seismic data, use mutual bracketing method and the reflection wave migration method of simple primary wave and reflection line-ups always, these methods are and directly identify reflection wave, the signal to noise ratio (S/N ratio) gathering reflection wave due to negative apparent velocity method is not high, and wave field is complicated; Reflection wave migration method is the achievement of a single track substantially.Explanation Accuracy when explanation personnel utilize these methods to work is not high, often easily produces larger error or mistake.Particularly for multilayer zone of interest inclination angle and tendency change large time, exist and carry out the work efficiency of data interpretation and all low problem of Explanation Accuracy.
Summary of the invention
The invention provides a kind of method of mine forward probe seismic data overlap-add procedure, solve the technical matters that reflection wave signal to noise ratio (S/N ratio) improves and the image quality of many zone of interest is low of negative apparent velocity method.
The present invention solves above technical matters by the following technical programs:
The ultimate principle that the technical solution used in the present invention utilizes TSP seismic data to carry out shunting collection stacking image is: the TSP earthquake reflective data (digital seismics magnetic recording) actual acquisition arrived, at the enterprising row relax of seismic processing system, obtains TSP earthquake superposition time section.Make an explanation processing the TSP stacked section obtained; Obtain the position of front, tunnel reflecting interface.The setting of TSP seismic data sampling rate: TSP time section of earthquake figure is the foundation drawing of process time-stacked section figure, and this figure is the original graph of actual acquisition; Generally process based on the original graph of actual acquisition, current process software be all for sampling rate be 1ms or 2ms establishment software; Mine gathers sampling rate general less 25us, 50us, 62.5us or 100us etc. of data, when process for sampling rate adopts, constant criterion is carried out, that is: S=V × T, in formula, S is earthquake wave trajectory, unit is rice, V is earthquake velocity of wave propagation, and T is the hourage of seismic event.The sampling rate of 62.5us is changed into 1ms to be equivalent to be multiplied by 16 on the right of above formula, the left side of above formula is multiplied by 16, distance adds 16 times in other words simultaneously, and frequency reduces 16 times.Can analogize for other sampling rates.The wave field separation of TSP seismic data: utilize the feature of TSP data up-downgoing ripple to be separated, obtain upward traveling wave.Upward traveling wave is carried out bin planning according to the CDP size of setting, carries out velocity sweeping superposition.
The technical scheme that the present invention takes for its technical matters of solution is: (1), TSP overlap-add procedure sampling rate: it is 1ms or 2ms that earthquake reflective data (digital seismics magnetic recording) the elapsed time unchangeability process conversion collected in tunnel or tunnel obtains Temporal sampling; (2), seismic data positive and negative apparent velocity wave field separation: utilize two-dimentional wave field separation software will to be separated process to the wave field containing positive and negative apparent velocity, obtain the seismic wave field of negative apparent velocity; (3), normal moveout correction with superpose: normal moveout correction is carried out to the reflection profile of negative apparent velocity, the data after normal moveout correction is taken out collection according to common reflection point, then superposes.
Further describe preferred version of the present invention below, be described with four parts of carrying out this method in Geovecteur Plus disposal system: the Temporal sampling of (1), existing disposal system is 1ms or 2ms mostly, 1ms or 2ms is defined as to earthquake data sampling rate and namely obtains earthquake data set A (x, t).The sampling rate of each seismic trace reset, specific practice is: TSP time section of earthquake figure is the foundation drawing of process time-stacked section figure, and this figure is the original graph of actual acquisition; Generally process based on the original graph of actual acquisition, current process software be all for sampling rate be 1ms or 2ms establishment software; Mine gathers sampling rate general less 25us, 50us, 62.5us or 100us etc. of data, when process for sampling rate adopts, constant criterion is carried out, that is: S=V × T, in formula, S is earthquake wave trajectory, unit is rice, V is earthquake velocity of wave propagation, and T is the hourage of seismic event.The sampling rate of 62.5us is changed into 1ms to be equivalent to be multiplied by 16 on the right of above formula, the left side of above formula is multiplied by 16, distance adds 16 times in other words simultaneously, and frequency reduces 16 times.Can analogize for other sampling rates.By reset sampling rate single-channel seismic data acquisition form two-dimension earthquake time data body T (x, y, t) together; (2), trace header T (x, y, t) is put to data volume and set up space attribute, by earthquake data set T (x, y, t), arrange according to road sequence, wave field separation is carried out to the data set arranging sequence, obtain reflected waveform data collection Ts (x, y, t), to Ts (x, y, t) reflected waveform data collection takes out collection Tsn (x, y, t) according to CRP gather; (3), normal moveout correction is carried out to reflection spot road collection Ts (x, y, t), then after many normal moveout correction, data set superposition Ts1 (x, y, t) forms stacked section, implement filtering output after superposition again and obtain achievement section Tstack (x, y, t); (4), the making of negative apparent velocity method interpretation result map: on GeoFrame seismic interpretation system, the comparative interpretation of layer position, tracking and time and depth transfer are carried out to two-dimentional poststack achievement section Tstack (x, y, t), explain and namely complete.
The invention has the beneficial effects as follows: adopt the present invention to utilize 2-d seismic data to carry out overlap-add procedure, 2-d seismic data can be completed and be treated to stacked profile map, and complete two-dimension earthquake and bear apparent velocity method interpretation result map; Reduce the labour intensity of explanation personnel, increase work efficiency and mapping precision, solve the problem of front, tunnel many zone of interest imaging simultaneously.Overcome complicated fault, change of pitch angle greatly on the impact of seismic data interpretation.Solve inclination angle between many zone of interest, tendency change large condition under high precision explain problem; Restriction applied widely, not complicated by tomography, inclination angle tendency change is large; The accuracy evaluating the structure exploration of front, tunnel can be improved.
Accompanying drawing explanation
Fig. 1 is the process flow diagram utilizing mine forward probe seismic data to complete the method for overlap-add procedure formation method;
Fig. 2 gathers master map;
Fig. 3 is X, the Y-component of acquired original data;
Fig. 4 is the radial component after polarization filtering;
Fig. 5 is the reflection wave after radial component wave field separation;
Fig. 6 is through the stacked profile map that mine forward probe seismic data completes overlap-add procedure formation method, and in the sectional view after overlap-add procedure, 64ms and 102ms place finds reflection wave;
In FIG: A, mine earthquake source book load; B, the polarization filtering of mine seismic data; C, sampling rate reset; D, deconvolution; The planning of E, road collection and overlap-add procedure; F, output time section.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in detail.
Embodiment 1:
Certain mine forward probe seismic data overlap-add procedure, this mine work area complex structure, in order to solve high precision one-tenth figure, verify the problem of subterranean coal distribution and the invention that produces, be an example of practical probation, treatment step is as follows:
(1), the Temporal sampling of existing disposal system is 1ms or 2ms mostly, is defined as 1ms or 2ms namely obtains earthquake data set A (x, t) to earthquake data sampling rate;
(2), by earthquake data set A (x, t), arrange according to road sequence, wave field separation is carried out to the data set arranging sequence, obtains reflected waveform data collection, according to CRP gather, collection is taken out to reflected waveform data collection;
(3), to reflection spot road collection carry out normal moveout correction, then after many normal moveout correction, data set superposition forms stacked section, implements filtering again after superposition;
(4) be, finally by processing profiles sampling rate according to the return to primitive age sampling rate discussed above, namely process completes.
Claims (1)
1. a mine forward probe seismic data superposition processing method, is characterized in that following steps:
(1), TSP overlap-add procedure sampling rate: it is 1ms or 2ms that earthquake reflective data (digital seismics magnetic recording) the elapsed time unchangeability process conversion collected in tunnel or tunnel obtains Temporal sampling;
(2), seismic data positive and negative apparent velocity wave field separation: utilize two-dimentional wave field separation software will to be separated process to the wave field containing positive and negative apparent velocity, obtain the seismic wave field of negative apparent velocity;
(3), normal moveout correction with superpose: normal moveout correction is carried out to the reflection profile of negative apparent velocity, the data after normal moveout correction is taken out collection according to common reflection point, then superposes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510124046.9A CN104965221A (en) | 2015-03-22 | 2015-03-22 | Mine advanced detection earthquake data superposing processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510124046.9A CN104965221A (en) | 2015-03-22 | 2015-03-22 | Mine advanced detection earthquake data superposing processing method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104965221A true CN104965221A (en) | 2015-10-07 |
Family
ID=54219261
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510124046.9A Pending CN104965221A (en) | 2015-03-22 | 2015-03-22 | Mine advanced detection earthquake data superposing processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104965221A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040088385A1 (en) * | 2002-11-01 | 2004-05-06 | Hexago Inc. | Method and apparatus for connecting IPV4 devices through an IPV6 network using a tunnel setup protocol |
| WO2009091979A2 (en) * | 2008-01-18 | 2009-07-23 | Baker Hughes Incorporated | Method for 3-c 3-d wavefield decomposition |
| CN101798923A (en) * | 2010-03-30 | 2010-08-11 | 煤炭科学研究总院重庆研究院 | System and method for remote control coal mine evacuation working face advance detection and forecasting |
| CN102495434A (en) * | 2011-11-25 | 2012-06-13 | 成都畅达通地下工程科技发展有限公司 | Advanced geological prediction method for underground engineering |
-
2015
- 2015-03-22 CN CN201510124046.9A patent/CN104965221A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040088385A1 (en) * | 2002-11-01 | 2004-05-06 | Hexago Inc. | Method and apparatus for connecting IPV4 devices through an IPV6 network using a tunnel setup protocol |
| WO2009091979A2 (en) * | 2008-01-18 | 2009-07-23 | Baker Hughes Incorporated | Method for 3-c 3-d wavefield decomposition |
| CN101798923A (en) * | 2010-03-30 | 2010-08-11 | 煤炭科学研究总院重庆研究院 | System and method for remote control coal mine evacuation working face advance detection and forecasting |
| CN102495434A (en) * | 2011-11-25 | 2012-06-13 | 成都畅达通地下工程科技发展有限公司 | Advanced geological prediction method for underground engineering |
Non-Patent Citations (1)
| Title |
|---|
| 杜立志: "隧道施工地质地震波法超前探测技术研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111239802B (en) | Deep learning speed modeling method based on seismic reflection waveform and velocity spectrum | |
| CN102540250B (en) | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method | |
| CN102395902B (en) | Seismic imaging systems and methods employing a fast target-oriented illumination calculation | |
| CN102937721B (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
| CN106094032B (en) | A kind of method for building formation velocity model | |
| WO2017024523A1 (en) | Inversion method for ray elastic parameter | |
| CN101915938B (en) | Offset imaging method and device for converted waves | |
| CN104570125A (en) | Method for improving imaging velocity model precision by utilizing well data | |
| NO20111499A1 (en) | Expanding coverage of VSP / CDP imaging using first-order descending multiples | |
| CN102841375A (en) | Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition | |
| CN103336297B (en) | Micro rupture vector scan method | |
| CN106443766A (en) | 3 dimensional tunnel earthquake advance seismic method | |
| CN105372705B (en) | A kind of strata slicing method based on many ripple data | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| WO2012064605A1 (en) | Microseismic source location estimation method with high resolution using green's functions | |
| CN105093281A (en) | Earthquake multi-wave modeling method under inverse framework | |
| CN107831542A (en) | DDW high accuracy depths domain well shake matching process | |
| CN103149588B (en) | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration | |
| CN102901984B (en) | Method for constructing true earth surface dip angle trace gathers of seismic data | |
| CN104142518A (en) | Method for analyzing pre-stack time migration response of seismological observation system | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN103645497A (en) | Emergence angle based multi-component wave field separation method | |
| CN104932021A (en) | Constrained tomography speed modeling method based on reverse ray tracing | |
| CN106199704A (en) | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method | |
| CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151007 |
|
| WD01 | Invention patent application deemed withdrawn after publication |