CN104730578A - Automatic interpretation method and system for shallow earthquake single reflection data - Google Patents
Automatic interpretation method and system for shallow earthquake single reflection data Download PDFInfo
- Publication number
- CN104730578A CN104730578A CN201510106684.8A CN201510106684A CN104730578A CN 104730578 A CN104730578 A CN 104730578A CN 201510106684 A CN201510106684 A CN 201510106684A CN 104730578 A CN104730578 A CN 104730578A
- Authority
- CN
- China
- Prior art keywords
- data
- reflecting interface
- earthquake
- props
- curve
- 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 an automatic interpretation method and system for shallow earthquake single reflection data. The automatic interpretation method comprises the steps that 1 data pre-processing is conducted, and in particular relates to that a first processing is conducted on single reflection earthquake data; 2 travel times of reflected waves are collected, and particularly the travel times of reflected waves from reflecting interfaces are collected; 3 an X2-T2 curve is drew; 4 data interpretation is conducted, particularly the single earthquake reflection data is interpreted; 5 an achievement is obtained and outputing is conducted, and particularly the interpretation result is generated and the result is output. The automatic interpretation method and system for the shallow earthquake single reflection data are mainly applied to the field of water-powering engineering, water conservancy project, oil exploration, highway, engineering investigation in railway and the like. According to the automatic interpretation method and system for the shallow earthquake single reflection data, the working efficiency is improved greatly, and the precision of the interpretation result of shallow earthquake single reflection data is improved at the same time.
Description
Technical field
The invention belongs to engineering geophysics field, relate to engineering geophysics Shallow seismic reflection method data analysis field, particularly relate to a kind of shallow earthquake list and prop up reflectance data automatic interpretation method.Be mainly used in the fields such as hydroelectric project, hydraulic engineering, petroleum prospecting, highway, railway engineering prospecting.
Background technology
Shallow seismic reflection method is the method that engineering geophysics is commonly used, economical and efficient, it comprises single reflection, reflection multi-fold etc., wherein single reflection is one of effective physical prospecting means of all kinds of engineering simple place cover thickness and geological layering detection, is also the important means of test each reflecting interface overlying strata average velocity parameter.
In practical implementation, because the transversal inhomogeneity of shallow surficial geology character changes, under in shortage, boring provides speed parameter situation, usually need units differ use shallow earthquake list prop up reflection to carry out cover thickness and geological layering detection, handled by seismic reflection multi-fold number and decipher formation physical property (speed) layering mark and formation velocity parameter are provided.At present, shallow earthquake list props up that reflectance data decipher is many provides informative formation parameter by Manual analysis, and efficiency is extremely low, and the precision of data interpretation is not high.
Summary of the invention
The object of the present invention is to provide a kind of shallow earthquake list to prop up reflectance data automatic interpretation method, effectively increase work efficiency, improve the precision that shallow earthquake list props up reflectance data interpretation result simultaneously, thus solve the foregoing problems existed in prior art.
To achieve these goals, the technical solution used in the present invention is as follows:
A kind of shallow earthquake list props up reflectance data automatic interpretation method, comprises the following steps:
S1, data prediction, is specially and carries out rough handling to single seismic reflection data;
S2, picks up each Travel time, is specially the Travel time of pickup from each reflecting interface;
S3, draws X
2-T
2curve;
S4, data interpretation, is specially decipher list and props up earthquake reflective data;
S5, achievement and output, be specially and generate interpretation result and Output rusults.
Preferably, step S1 comprises the following steps:
S11, Data Format Transform, the original earthquake data by seismic instrument record is converted to required RPS formatted data;
S12, editing files head, namely edits the File header information in RPS formatted data after format conversion, and described File header information comprises: survey line number, excite pile No., receive pile No..
Preferably, step S2 is specially, and in oscillogram, according to the time-distance equation from each reflecting interface reflection wave, by the comparative analysis of ripple group, determines the reflection wave from different reflecting interface, and picks up Travel time.
Preferably, step S3 is specially, and first calculates each acceptance point respectively and is designated as X to the square value of shot point distance
2, each acceptance point is designated as T from the square value of the Travel time of different reflecting interface
2, then with X
2for horizontal ordinate, with T
2for ordinate is depicted as X
2-T
2curve.
Preferably, step S4 comprises the following steps:
S41, the average velocity of wave in decipher overlying strata, namely selects (the X of different reflecting interface
2, T
2) coordinate points asks for the average earthquake velocity of longitudinal wave of each reflecting interface overlying strata respectively;
S42, decipher reflecting interface buried depth, namely according to the buried depth of the different reflecting interface of each measuring point of the Geometrical propagation rule decipher of reflection wave.
Preferably, step S5 comprises the following steps:
S51, profiling, is about to the scale engineering sectional view of different reflecting interface buried depths by regulation of each measuring point obtained according to decipher;
S52, mark stratum average velocity, namely marks the information of seismic velocity fifty-fifty of various places layer on engineering sectional view;
S53, result output, namely selects, then exports in time section and engineering section.
Preferably, Travel time and X
2-T
2dynamically associating between curve dynamically updates X after making amendment Travel time in time
2-T
2curve.
Preferably, pickup or amendment Travel time can the T-X curves of this reflecting interface of dynamic drafting in time, the rationality of picking up when being reflected away by the hyperbolic curve form inspection by attributes of each reflecting interface T-X curve and correctness.
A kind of shallow earthquake list props up reflectance data automatic interpretation system, comprising:
Data preprocessing module, for carrying out rough handling to single seismic reflection data;
Each Travel time pickup model, for picking up the Travel time from each reflecting interface;
X
2-T
2drawing of Curve module, for drawing X
2-T
2curve;
Data interpretation module, props up earthquake reflective data for decipher list;
Achievement and output module, for generating interpretation result and Output rusults.
The invention has the beneficial effects as follows:
The present invention is mainly used in the fields such as hydroelectric project, hydraulic engineering, petroleum prospecting, highway, railway engineering prospecting.Prop up reflectance data automatic interpretation method and system by shallow earthquake list of the present invention, effectively increase work efficiency, improve the precision that shallow earthquake list props up reflectance data interpretation result simultaneously.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the embodiment of the present invention;
Fig. 2 picks up each reflecting interface Travel time and draws T-X curve figure in the embodiment of the present invention;
Fig. 3 is the X in the embodiment of the present invention
2-T
2curve, each reflecting interface overlying strata average velocity figure;
Fig. 4 is the decipher result map of the embodiment of the present invention.
Embodiment
In order to make object of the present invention, technical scheme and advantage clearly understand, below in conjunction with accompanying drawing, the present invention is further elaborated.Should be appreciated that embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
As shown in drawings, the invention discloses a kind of shallow earthquake list and prop up reflectance data automatic interpretation method and system, the method comprises the following steps:
1, data prediction A, carries out rough handling to single seismic reflection data;
2, pick up each interface Travel time B, pick up the Travel time from each reflecting interface;
3, X is drawn
2-T
2curve C, draws X
2-T
2curve;
4, data interpretation D, decipher list props up earthquake reflective data;
5, achievement and output E, generate interpretation result and Output rusults.
Data prediction A comprises Data Format Transform A1, editing files head A2 two subprocess.Wherein the original earthquake data of seismic instrument record is converted to " physical prospecting test data overall analysis system " specific RPS formatted data by Data Format Transform A1; Editing files head A2 is to the survey line in RPS formatted data after format conversion number, excite pile No., receive the File header information such as pile No. edits.
Pickup Travel time B is in oscillogram, according to the time-distance equation from each reflecting interface reflection wave, by the comparative analysis of ripple group, determines the reflection wave from different reflecting interface, and picks up Travel time.
Draw X
2-T
2curve C first calculates X respectively
2(each acceptance point is to the square value of shot point distance), T
2the value of (each acceptance point is from the square value of the Travel time on difference reflection circle), then thinks X
2horizontal ordinate, with T
2for ordinate is depicted as X
2-T
2curve.
Data interpretation D comprises the average velocity of wave D1 in decipher overlying strata and decipher reflecting interface buried depth D2 two subprocess, and wherein the average velocity of wave D1 in decipher overlying strata is the different reflecting interface (X of choose reasonable
2, T
2) coordinate points asks for the average earthquake velocity of longitudinal wave of each reflecting interface overlying strata respectively; Decipher reflecting interface buried depth D2 is the buried depth according to the different reflecting interface of each measuring point of the Geometrical propagation rule decipher of reflection wave.When earthquake list props up method for reflection for simple place exploration, the average velocity of wave D1 in decipher overlying strata is a critical process and the step of data interpretation D, and object is each reflection spot buried depth of decipher; When earthquake list props up method for reflection for providing the parameter of earthquake multi-fold data processing and decipher, the average velocity of wave D1 in decipher overlying strata is also one of object of data interpretation D.
Achievement and export E and comprise profiling E1, mark stratum average velocity E2 and result output E3 tri-subprocess, wherein profiling E1 is by the different reflecting interface buried depths of each measuring point that obtain according to the decipher scale engineering sectional view by regulation; Mark stratum average velocity E2 is the information of seismic velocity fifty-fifty marking various places layer on engineering sectional view; Result output E3 selects in time section and engineering section, then exports.
Travel time and X
2-T
2dynamically associating between curve dynamically updates X after making amendment Travel time in time
2-T
2curve; At X
2-T
2curve is selected the (X of different reflecting interface
2, T
2) solution can be followed the tracks of fast after point and translate this reflecting interface overlying strata average speed value, work efficiency is significantly higher than traditional analytical approach.
Pickup or amendment Travel time can the T-X curves of dynamic drafting (renewal) this reflecting interface in time, the rationality of picking up when being reflected away by the hyperbolic curve form inspection by attributes of each reflecting interface T-X curve and correctness.
Composition graphs 1 is described in detail to technical solution of the present invention by specific embodiment, as follows:
The present embodiment is that some hydropower station pre feasibility study stage W1 section 0-48 pile No. section shallow earthquake list props up reflectance data analysis.
The present embodiment comprises: data prediction A, pick up each interface Travel time B, draw X
2-T
2curve C, data interpretation D, achievement and output E five processes.
1, data prediction A:
The source data of the present embodiment is the SEG2 formatted data of U.S. R24 seismograph record, be " physical prospecting test data overall analysis system " specific RPS formatted data through format conversion, then to survey line number, excite pile No., receive the File header information such as pile No. and edit.
2, each interface Travel time B is picked up:
This example has picked up the Travel time (on accompanying drawing 2) of three reflecting interfaces.Present the different monolateral hyperbolic curve form of curvature from each reflecting interface reflection wave, by the comparative analysis of ripple group, determine the reflection wave from different reflecting interface, and pick up Travel time; While pickup Travel time, with shot point pile No. be initial point, reception pile No. is horizontal ordinate, Travel time draws the T-X curve (accompanying drawing 2 times) of corresponding reflection circle for ordinate.Dynamically associate between Travel time and T-X curve, namely in ripple figure, change Travel time, T-X curve then dynamically updates (redrawing).
3, X is drawn
2-T
2curve C:
This example depicts the X of three reflecting interfaces
2-T
2curve (under Fig. 3).The kinematical equation of reflection wave is hyperbolic function, is initial point with shot point pile No., receives the square value (X of pile No.
2) be the square value (T of horizontal ordinate, Travel time
2) be ordinate drafting X
2-T
2curve, object is that quadratic function has been converted linear function, is transformed to linear relationship by hyperbolic relation.Travel time and X
2-T
2dynamically associate between curve, namely in ripple figure, change Travel time, X
2-T
2curve then dynamically updates (redrawing).
4, data interpretation D:
The buried depth (accompanying drawing 4) of this example decipher overlying strata average velocity of three reflecting interfaces (accompanying drawing 3 times), two reflecting interfaces.At X
2-T
2curve is chosen the (X at same reflection interface
2, T
2) coordinate points, go out straight line with least square fitting, the subduplicate inverse of straight slope is exactly the average velocity of this reflecting interface overlying strata, and it is decipher reflecting interface buried depth and earthquake multi-fold data processing and the requisite parameter of decipher; According to reflection wave Geometrical propagation rule, solution the buried depth of each reflection spot can be translated.By each reflecting interface overlying strata average velocity of decipher and the buried depth of reflecting interface, can calculate each formation velocity according to Dix formula, formation velocity is the physical parameter of requirement of engineering.Reflection is positioned at the midpoint of shot point and acceptance point, and the reflecting interface length of therefore decipher is only the half of shot point to acceptance point distance farthest.
5, achievement and output E
This example is the engineering sectional view (Fig. 4) drawn according to interpretation result, and the speed parameter marked in figure is the average velocity on each stratum.In real work, design department is it is desirable that engineering section, and time section is also one of achievement submitted to, therefore this process free section and engineering section two options on result output.
By adopting technique scheme disclosed by the invention, obtain effect useful as follows:
The present invention is mainly used in the fields such as hydroelectric project, hydraulic engineering, petroleum prospecting, highway, railway engineering prospecting.Prop up reflectance data automatic interpretation method and system by shallow earthquake list of the present invention, effectively increase work efficiency, improve the precision that shallow earthquake list props up reflectance data interpretation result simultaneously.
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should look protection scope of the present invention.
Claims (9)
1. shallow earthquake list props up a reflectance data automatic interpretation method, it is characterized in that, comprises the following steps:
S1, data prediction, is specially and carries out rough handling to single seismic reflection data;
S2, picks up each Travel time, is specially the Travel time of pickup from each reflecting interface;
S3, draws X
2-T
2curve;
S4, data interpretation, is specially decipher list and props up earthquake reflective data;
S5, achievement and output, be specially and generate interpretation result and Output rusults.
2. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, and it is characterized in that, step S1 comprises the following steps:
S11, Data Format Transform, the original earthquake data by seismic instrument record is converted to required RPS formatted data;
S12, editing files head, namely edits the File header information in RPS formatted data after format conversion, and described File header information comprises: survey line number, excite pile No., receive pile No..
3. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, it is characterized in that, step S2 is specially, in oscillogram, according to the time-distance equation from each reflecting interface reflection wave, by the comparative analysis of ripple group, determine the reflection wave from different reflecting interface, and pick up Travel time.
4. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, and it is characterized in that, step S3 is specially, and first calculates each acceptance point respectively and is designated as X to the square value of shot point distance
2, each acceptance point is designated as T from the square value of the Travel time of different reflecting interface
2, then with X
2for horizontal ordinate, with T
2for ordinate is depicted as X
2-T
2curve.
5. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, and it is characterized in that, step S4 comprises the following steps:
S41, the average velocity of wave in decipher overlying strata, namely selects (the X of different reflecting interface
2, T
2) coordinate points asks for the average earthquake velocity of longitudinal wave of each reflecting interface overlying strata respectively;
S42, decipher reflecting interface buried depth, namely according to the buried depth of the different reflecting interface of each measuring point of the Geometrical propagation rule decipher of reflection wave.
6. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, and it is characterized in that, step S5 comprises the following steps:
S51, profiling, is about to the scale engineering sectional view of different reflecting interface buried depths by regulation of each measuring point obtained according to decipher;
S52, mark stratum average velocity, namely marks the information of seismic velocity fifty-fifty of various places layer on engineering sectional view;
S53, result output, namely selects, then exports in time section and engineering section.
7. shallow earthquake list according to claim 1 props up reflectance data automatic interpretation method, it is characterized in that, Travel time and X
2-T
2dynamically associating between curve dynamically updates X after making amendment Travel time in time
2-T
2curve.
8. shallow earthquake list according to claim 7 props up reflectance data automatic interpretation method, it is characterized in that, pickup or amendment Travel time can the T-X curves of this reflecting interface of dynamic drafting in time, the rationality of picking up when being reflected away by the hyperbolic curve form inspection by attributes of each reflecting interface T-X curve and correctness.
9. shallow earthquake list props up a reflectance data automatic interpretation system, it is characterized in that, comprising:
Data preprocessing module, for carrying out rough handling to single seismic reflection data;
Each Travel time pickup model, for picking up the Travel time from each reflecting interface;
X
2-T
2drawing of Curve module, for drawing X
2-T
2curve;
Data interpretation module, props up earthquake reflective data for decipher list;
Achievement and output module, for generating interpretation result and Output rusults.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510106684.8A CN104730578A (en) | 2015-03-11 | 2015-03-11 | Automatic interpretation method and system for shallow earthquake single reflection data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510106684.8A CN104730578A (en) | 2015-03-11 | 2015-03-11 | Automatic interpretation method and system for shallow earthquake single reflection data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104730578A true CN104730578A (en) | 2015-06-24 |
Family
ID=53454649
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510106684.8A Pending CN104730578A (en) | 2015-03-11 | 2015-03-11 | Automatic interpretation method and system for shallow earthquake single reflection data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104730578A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020116160A1 (en) * | 2000-12-18 | 2002-08-22 | Patrick Lailly | Method intended to obtain reflection travel times from an interpretation of migrated cylindrical wave seismic data |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN102565852A (en) * | 2010-12-08 | 2012-07-11 | 同济大学 | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir |
-
2015
- 2015-03-11 CN CN201510106684.8A patent/CN104730578A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020116160A1 (en) * | 2000-12-18 | 2002-08-22 | Patrick Lailly | Method intended to obtain reflection travel times from an interpretation of migrated cylindrical wave seismic data |
| US7069200B2 (en) * | 2000-12-18 | 2006-06-27 | Institut Francais Du Petrole | Method intended to obtain reflection travel times from an interpretation of migrated cylindrical wave seismic data |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN102565852A (en) * | 2010-12-08 | 2012-07-11 | 同济大学 | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir |
Non-Patent Citations (2)
| Title |
|---|
| 北京克浪石油技术有限公司: "《地震采集工程软件系统用户手册》", 15 May 2001 * |
| 施俊杰，等: "井间地震直达波和反射波走时联合层析反演", 《中国海洋大学学报》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102176054B (en) | Near-surface comprehensive information processing explanation method | |
| CN104181596B (en) | Geologic horizon automatic tracking method and device | |
| CN103424777B (en) | A kind of method that improves seismic imaging resolution ratio | |
| CN100549731C (en) | Demodulator probe secondary localization method | |
| CN101315427A (en) | Method and system for processing seismic exploration data of complex area | |
| CN102221709A (en) | Velocity analysis and dynamic correction method based on formation parameter information | |
| CN105259570A (en) | Seaborne time shifting earthquake receiving point displacement correction method | |
| CN102081167B (en) | Picking method of three-dimensional VSP data first-motion wave | |
| CN105259581B (en) | Seismic data time-depth conversion method | |
| CN107728227B (en) | A kind of well pattern maturation zone quick discrimination buried channel sand body method | |
| CN108121009A (en) | Development late stage complex fault block velocity-variant depth mapping method | |
| CN102043167A (en) | Curved earth surface seismic data velocity analysis method | |
| CN101452082A (en) | First arrival picking -up method for fractal seismic waves | |
| CN105866834B (en) | Single branch refraction data decomposition method and system based on automatic identification refracted wave starting point | |
| CN103513279B (en) | A kind of illumination analysis computing method based on seismic wave equation and calculation element | |
| CN104730578A (en) | Automatic interpretation method and system for shallow earthquake single reflection data | |
| CN103529473A (en) | Compensation method for differential compaction in seismic inversion process | |
| CN113420348B (en) | Method for rapidly drawing uranium mine exploration line profile map | |
| CN103984021B (en) | Method for analyzing migration velocity based on series of velocity error sections | |
| CN101609164B (en) | Method for conducting horizon calibration of ground converted wave data by using VSP converted wave data | |
| CN107356970B (en) | High-precision well seismic data matching method | |
| CN104280766B (en) | A kind of direct offset method of utilization local data lineups slope | |
| CN1523373A (en) | Static correction method for uplink penetration transformed wave computation transverse wave | |
| CN104199107A (en) | Vertical seismic multi-wave data based predication method and system for depth before drilling | |
| Cummings | Determination of depths to an irregular interface in shallow seismic refraction surveys using a pocket calculator |
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 |
Application publication date: 20150624 |
|
| RJ01 | Rejection of invention patent application after publication |