CN103217712A - Earthquake velocity analysis method along geologic structure trend - Google Patents
Earthquake velocity analysis method along geologic structure trend Download PDFInfo
- Publication number
- CN103217712A CN103217712A CN2012100171129A CN201210017112A CN103217712A CN 103217712 A CN103217712 A CN 103217712A CN 2012100171129 A CN2012100171129 A CN 2012100171129A CN 201210017112 A CN201210017112 A CN 201210017112A CN 103217712 A CN103217712 A CN 103217712A
- Authority
- CN
- China
- Prior art keywords
- velocity
- velocity analysis
- tectonic structure
- along
- elliptical area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an earthquake velocity analysis method along a geologic structure trend. The earthquake velocity analysis method along the geologic structure trend includes the steps of determining the geologic structure trend in a velocity analysis area, enabling velocity analysis surface element distribution to fit in with the geologic structure trend; centering on a single velocity analysis point for the single velocity analysis point, increasing the trace gathering number along the geologic structure trend, selecting an optimum elliptical face element combination short-shaft length; and increasing the trace gathering number along the geologic structure trend, and selecting the optimum elliptical face element combination short-shaft length. The earthquake velocity analysis method solves the problem in the prior art that the influence of the underground geologic structure on the velocity analysis is ignored, and has the advantages of strengthening focus capacity of a velocity spectrum of a place with a complex geologic structure, being capable of improving accuracy of velocity pickup, and thus improving stacking effect and migration imaging effect of the complex geologic structure.
Description
Technical field
The present invention relates to the seismic data process field, particularly relate to a kind of seismic velocity analysis methods.
Background technology
Seismic imaging is more and more higher to the analysis precision requirement of speed, and the combination by optimal speed analysis bin can improve degree of covering, improves the quality of velocity spectrum.Rectangular surfaces unit combination mode is generally chosen in conventional speeds analysis at present, and the choosing just simply along ordinate or crossline direction of the long axis direction of bin combination, has ignored the influence of subsurface geological structure to the velocity analysis bin.For this reason, we have invented and a kind ofly new have moved towards seismic velocity analysis methods along tectonic structure, have solved above problem.
Summary of the invention
The purpose of this invention is to provide a kind of improve complex geological structure velocity analysis precision move towards seismic velocity analysis methods along tectonic structure.
Purpose of the present invention can realize by following technical measures: move towards seismic velocity analysis methods along tectonic structure, this moves towards seismic velocity analysis methods along tectonic structure and comprises: determine the tectonic structure trend in the velocity analysis district, it is consistent with this tectonic structure trend that the velocity analysis bin is distributed; Is the center to single velocity analysis point with this velocity analysis point, increases road collection number along the tectonic structure tendency, chooses the first combination of best elliptical area minor axis length; And move towards to increase road collection number along this tectonic structure, choose best elliptical area unit combination long axis length.
Purpose of the present invention also can realize by following technical measures:
In making the velocity analysis bin distribution step consistent with this tectonic structure trend, this single velocity analysis point adopts the combination of elliptical area unit, each long axis of ellipse direction is consistent with this tectonic structure trend in this elliptical area unit combination, and short-axis direction is vertical with this tectonic structure trend.
In the step of choosing the first combination of best elliptical area minor axis length, fixing and moving towards direction road collection number along this tectonic structure is 1, increase road collection number along tectonic structure tendency direction again, generate a plurality of velocity spectrums, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these a plurality of velocity spectrums, the road collection number that produces this velocity spectrum is the minor axis length of this best elliptical area unit combination.
In the step of choosing the first combination of best elliptical area long axis length, fixing is this best elliptical area unit combination minor axis length along this tectonic structure tendency direction road collection number, move towards to increase road collection number along this tectonic structure again, generate a plurality of velocity spectrums, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these a plurality of velocity spectrums, the road collection number that produces this velocity spectrum is the long axis length of this best elliptical area unit combination.
Moving towards seismic velocity analysis methods along tectonic structure and analyzed of the influence of the trend of subsurface geological structure among the present invention to velocity analysis bin direction, improve the velocity analysis precision at complex geological structure place, improved the velocity spectrum resolution at complex geological structure place.
Description of drawings
Fig. 1 is for moving towards the process flow diagram of seismic velocity analysis methods along tectonic structure in one embodiment of the invention;
Fig. 2 left side is for adopting the velocity analysis bin distribution schematic diagram of one embodiment of the invention, and Fig. 2 right side is that conventional speeds is analyzed the bin distribution schematic diagram;
The synoptic diagram that Fig. 3 chooses for the road collection number of the major axis of single velocity analysis null ellipse bin combination in one embodiment of the invention and short-axis direction;
Fig. 4 is the velocity spectrum comparison diagram that conventional speeds analysis and the velocity analysis of adopting one embodiment of the invention produce in one embodiment of the invention.
Embodiment
For above-mentioned and other purposes, feature and advantage of the present invention can be become apparent, cited below particularlyly go out preferred embodiment, and cooperate appended graphicly, be described in detail below.
Fig. 1 is for moving towards the process flow diagram of seismic velocity analysis methods along tectonic structure in one embodiment of the invention.In step 101, according to the tectonic structure trend in definite analysis area such as figure such as T0 such as grade of work area zone of interest.Flow process enters into step 102.
In step 102, single velocity analysis point adopts the combination of elliptical area unit, and it is consistent with the tectonic structure trend that fixed speed is analyzed elliptical area unit combination long axis direction, and the ellipse short shaft direction is vertical with the tectonic structure trend.Flow process enters into step 103.
In step 103, it is 1 that fixed speed is analyzed elliptical area unit combination long axis direction road collection number.Flow process enters into step 104.
In step 104, constantly increase elliptical area unit combination short-axis direction road collection number, and generate corresponding
A series of velocity spectrums.Flow process enters into step 105.
In step 105, from a series of velocity spectrums that generate, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution, the road collection number that produces this velocity spectrum promptly is defined as the best minor axis length of elliptical area unit combination, and fixedly the road collection number of the short-axis direction of elliptical area unit combination is best minor axis length.Flow process enters into step 106.
In step 106, constantly increase elliptical area unit combination long axis direction road collection number, and generate corresponding a series of velocity spectrum.Flow process enters into step 106.
In step 107, from a series of velocity spectrums that generate, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution, the road collection number that produces this velocity spectrum promptly is defined as the best long axis length of elliptical area unit combination, and fixedly the road collection number of the short-axis direction of elliptical area unit combination is best long axis length.Can carry out along the SEISMIC VELOCTTY ANALYSIS AND of tectonic structure trend with the road collection formation speed spectrum that step 105 and step 107 are determined.Flow process finishes.
As shown in Figure 2, Fig. 2 left side is for adopting the velocity analysis bin distribution schematic diagram of one embodiment of the invention, and Fig. 2 right side is that conventional speeds is analyzed the bin distribution schematic diagram.Wherein, curve is the tectonic structure trend among the figure.Among the present invention to move towards seismic velocity analysis methods along tectonic structure be the tectonic structure trend of determining earlier in the velocity analysis district, the velocity analysis bin is distributed consistent with the tectonic structure trend (long axis direction and the tectonic structure of the elliptical area unit combination of velocity analysis move towards consistent), determine the minor axis of elliptical area unit combination and the road collection number of long axis direction again.
The synoptic diagram that Fig. 3 chooses for the road collection number of the minor axis of single velocity analysis null ellipse bin combination in one embodiment of the invention and long axis direction.Wherein, middle black circle is the velocity analysis point, and lattice is a common midpoint.Left side figure is that first fixed length shaft length is 1 to choose minor axis, and right side figure is that fixedly minor axis length is 3 to choose major axis.Among the present invention to move towards seismic velocity analysis methods along tectonic structure be the center with the velocity analysis point, fixing earlier to move towards direction road collection number along tectonic structure be 1, increases road collection number (seeing Fig. 3 left side) along tectonic structure tendency direction again, generates a series of velocity spectrums.Choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these velocity spectrums, the road collection number that produces this velocity spectrum is the minor axis length of elliptical area unit combination.Fixing is minor axis length along tectonic structure tendency direction road collection number, moves towards to increase road collection number (seeing Fig. 3 right side) along tectonic structure again, generates a series of velocity spectrums.Choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these velocity spectrums, the road collection number that produces this degree spectrum of speed is the long axis length of elliptical area unit combination.
Fig. 4 is conventional speeds analysis and the speed that adopts one embodiment of the invention in one embodiment of the invention
The velocity spectrum comparison diagram that degree analysis produces.As can be seen from Figure 4, move towards the focusing that seismic velocity analysis methods has strengthened the velocity spectrum at complex geological structure place along tectonic structure among the present invention, can improve the precision of velocity pick, and then improve the stack and the migration imaging effect of complex geological structure.
Claims (4)
1. move towards seismic velocity analysis methods along tectonic structure, it is characterized in that, this moves towards seismic velocity analysis methods along tectonic structure and comprises:
Determine the tectonic structure trend in the velocity analysis district, it is consistent with this tectonic structure trend that the velocity analysis bin is distributed;
Is the center to single velocity analysis point with this velocity analysis point, increases road collection number along the tectonic structure tendency, chooses the first combination of best elliptical area minor axis length; And
Move towards to increase road collection number along this tectonic structure, choose the first combination of best elliptical area long axis length.
2. according to claim 1ly move towards seismic velocity analysis methods along tectonic structure, it is characterized in that, in making the velocity analysis bin distribution step consistent with this tectonic structure trend, this single velocity analysis point adopts the combination of elliptical area unit, each long axis of ellipse direction is consistent with this tectonic structure trend in this elliptical area unit combination, and short-axis direction is vertical with this tectonic structure trend.
3. according to claim 1ly move towards seismic velocity analysis methods along tectonic structure, it is characterized in that, in the step of choosing the first combination of best elliptical area minor axis length, fixing and moving towards direction road collection number along this tectonic structure is 1, increase road collection number along tectonic structure tendency direction again, generate a plurality of velocity spectrums, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these a plurality of velocity spectrums, the road collection number that produces this velocity spectrum is the minor axis length of this best elliptical area unit combination.
4. according to claim 1ly move towards seismic velocity analysis methods along tectonic structure, it is characterized in that, in the step of choosing the first combination of best elliptical area long axis length, fixing is this best elliptical area unit combination minor axis length along this tectonic structure tendency direction road collection number, move towards to increase road collection number along this tectonic structure again, generate a plurality of velocity spectrums, choose signal to noise ratio (S/N ratio) and the highest velocity spectrum of resolution from these a plurality of velocity spectrums, the road collection number that produces this velocity spectrum is the long axis length of this best elliptical area unit combination.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210017112.9A CN103217712B (en) | 2012-01-19 | 2012-01-19 | Earthquake velocity analysis method along geologic structure trend |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210017112.9A CN103217712B (en) | 2012-01-19 | 2012-01-19 | Earthquake velocity analysis method along geologic structure trend |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103217712A true CN103217712A (en) | 2013-07-24 |
| CN103217712B CN103217712B (en) | 2015-06-17 |
Family
ID=48815664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210017112.9A Active CN103217712B (en) | 2012-01-19 | 2012-01-19 | Earthquake velocity analysis method along geologic structure trend |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103217712B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884693A (en) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | Adaptively move towards normal-moveout spectrum acquiring method and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2148838C1 (en) * | 1999-06-16 | 2000-05-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Seismic data processing technique |
| CN101290355A (en) * | 2007-04-17 | 2008-10-22 | 中国石油集团东方地球物理勘探有限责任公司 | Method for determining carbonate rock hydrocarbon reservoir structure |
-
2012
- 2012-01-19 CN CN201210017112.9A patent/CN103217712B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2148838C1 (en) * | 1999-06-16 | 2000-05-10 | Нижне-Волжский научно-исследовательский институт геологии и геофизики | Seismic data processing technique |
| CN101290355A (en) * | 2007-04-17 | 2008-10-22 | 中国石油集团东方地球物理勘探有限责任公司 | Method for determining carbonate rock hydrocarbon reservoir structure |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884693A (en) * | 2019-01-18 | 2019-06-14 | 潜能恒信能源技术股份有限公司 | Adaptively move towards normal-moveout spectrum acquiring method and system |
| CN109884693B (en) * | 2019-01-18 | 2021-10-15 | 潜能恒信能源技术股份有限公司 | Self-adaptive trend velocity spectrum solving method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103217712B (en) | 2015-06-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102221709B (en) | Velocity analysis and dynamic correction method based on formation parameter information | |
| CN102455439B (en) | Diffracted wave field separation method based on Kirchhoff integral method | |
| CN102053261B (en) | Method for processing seismic data | |
| CN102928873B (en) | Method for positioning ground micro-seismic based on four-dimensional energy focusing | |
| CN1873441A (en) | Method for suppressing background noise not related to excitation source in seismic data enhancement | |
| CN101598808B (en) | Method for improving image quality of seismic data | |
| CN101598803B (en) | Method for directly obtaining stacked section of converted wave | |
| CN105607124A (en) | Seismic-wave near-surface stratum quality factor compensation method and device | |
| CN103645497B (en) | Emergence angle based multi-component wave field separation method | |
| CN102565857A (en) | Automatic remaining dynamic correction method | |
| CN104280777A (en) | Method for suppressing interference of seismic data multiples on land | |
| Chambers | Using DAS to investigate traffic patterns at Brady Hot Springs, Nevada, USA | |
| CN102053260B (en) | Method for acquiring azimuth velocity of primary wave and method for processing earthquake data | |
| CN104133244B (en) | A kind of near surface static corrections processing method of focus shear wave exploration | |
| CN103576197A (en) | Method for extracting converted wave angle channel set | |
| CN103217712A (en) | Earthquake velocity analysis method along geologic structure trend | |
| CN103645498B (en) | A kind of method of direct generation common-conversion point gather | |
| CN104765066B (en) | The scanning focused imaging method of earthquake three-dimensional velocity of wave | |
| CN110082820B (en) | Explosive source mixed distributed broadband excitation method | |
| CN104536048B (en) | A kind of automatic excision seismic wave ripple method extremely | |
| CN106338762B (en) | Anti- focal imaging method is rebuild based on the model trace that data field is oriented to | |
| CN107797140B (en) | Beading seismic response features focus optimization method and system | |
| Wu et al. | Offset VSP P wave and converted S wave salt proximity migration study | |
| CN104166157A (en) | Transverse wave processing method based on explosive source excitation three-component receiving on special condition | |
| CN105093316B (en) | Declined geological data static correcting method along the surface array of well track search |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |