CN106019367A - Method for improving coherent energy spectrum resolution ratio and modeling method for increasing prestack depth migration speed - Google Patents
Method for improving coherent energy spectrum resolution ratio and modeling method for increasing prestack depth migration speed Download PDFInfo
- Publication number
- CN106019367A CN106019367A CN201610333110.9A CN201610333110A CN106019367A CN 106019367 A CN106019367 A CN 106019367A CN 201610333110 A CN201610333110 A CN 201610333110A CN 106019367 A CN106019367 A CN 106019367A
- Authority
- CN
- China
- Prior art keywords
- depth
- coherent energy
- intervals
- energy spectrum
- modal 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.)
- Granted
Links
- 230000001427 coherent effect Effects 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000001228 spectrum Methods 0.000 title claims abstract description 42
- 238000013508 migration Methods 0.000 title claims abstract description 36
- 230000005012 migration Effects 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for improving coherent energy spectrum resolution ratio and a modeling method for increasing prestack depth migration speed. The method for improving the coherent energy spectrum resolution ratio includes scanning a prestack depth migration imaging channel set and determining coherent energy spectrum data of a depth domain channel set; dividing the depth domain channel set into a plurality of depth sections; performing n power and absolute value operation on coherent energy spectrum data of one specific depth section and all depth sections with depths greater than that of the specific depth section, and obtaining new coherent energy spectrum data of the specific depth section and all the depth sections with depths greater than that of the specific depth section, wherein n is a positive number greater than or equal to 1 and the n corresponding to the coherent energy spectrum of the depth sections subjected to N power operation increases along with the increase of the depths of the depth sections. By adopting the method provided by the invention, the depth domain channel set coherent energy spectrum with the good resolution ratio can be obtained and the formation of a more reasonable prestack depth migration speed model is facilitated.
Description
Technical field
The present invention relates to technical field of geophysical exploration, specifically, relate to a kind of raising seismic data
The method of pre-stack depth migration imaging road collection coherent energy spectral resolution and a kind of pre-stack depth migration speed
Degree modeling method, is particularly advantageous in improving the resolution of different depth coherent energy spectrum.
Background technology
In prior art, in process of seismic data processing, need to carry out pre-stack depth migration velocity modeling.
During pre-stack depth migration velocity modeling, by seismic data pre-stack depth migration imaging road collection is entered
Row scanning, forms coherent energy spectrum.
But, in actual process, owing to the coherent energy spectrum of Depth Domain road collection is along with the change of the degree of depth
Changing and change, when the degree of depth increases, the resolution of coherent energy spectrum is more and more lower.This is the energy that is concerned with at present
The greatest difficulty that amount spectrum scanning is faced.
Summary of the invention
It is an object of the invention to solve at least one in deficiencies of the prior art.
The problems referred to above existed for prior art, inventor, through research, finds: due to prestack depth
The feature of migration imaging road collection self, along with the increase of the degree of depth, the frequency of road collection reduces, thus causes sweeping
Retouch the coherent energy spectrum valid value range obtained wider, it is difficult to distinguish and pick up.
In view of this, the present invention proposes a kind of seismic data pre-stack depth migration imaging road collection that improves and is concerned with
The method of energy spectrum resolution.Described method includes: to seismic data pre-stack depth migration gained imaging road
Collection is scanned, and obtains determining the coherent energy modal data of Depth Domain road collection；The side from shallow to deep along stratum
It is divided into multiple depth intervals to by described Depth Domain road collection；To in the plurality of depth intervals a certain deeply
The coherent energy modal data in degree interval and all depth intervals more deeper than described a certain depth intervals
(that is, more deeper all other depth intervals than described a certain depth intervals in the plurality of depth intervals)
Coherent energy modal data carry out n power computing and signed magnitude arithmetic(al) respectively, obtain described a certain degree of depth district
Between new coherent energy modal data and all depth intervals more deeper than described a certain depth intervals
New coherent energy modal data, wherein, all degree n are the positive number not less than 1, and carry out n time
The degree n corresponding to coherent energy modal data in the different depth interval of side's computing is along with the depth intervals degree of depth
Intensification and increase.
Additionally, present invention also offers a kind of pre-stack depth migration velocity modeling method.Described modeling method
Including seismic data carrying out a series of process to form the degree of depth-rate pattern, and described modeling method is adopted
The different depth improved obtained by the collection scanning of pre-stack depth migration imaging road by method as above is concerned with
The resolution of energy spectrum.
Compared with prior art, the beneficial effect comprise that and be obtained in that there is good discrimination power
The coherent energy spectrum of Depth Domain road collection, is particularly advantageous in improving different depth and the deeper Depth Domain road of the degree of depth
The resolution of the coherent energy spectrum of collection；Advantageously form more reasonably pre-stack depth migration rate pattern.
Detailed description of the invention
Hereinafter, exemplary embodiment will be combined to describe the raising seismic data prestack of the present invention in detail
The method of depth migration imaging road collection coherent energy spectral resolution and pre-stack depth migration velocity modeling method.
In one exemplary embodiment of the present invention, improve seismic data pre-stack depth migration imaging road collection
The method of coherent energy spectral resolution comprises the steps that
(1) the imaging road collection formed according to seismic data pre-stack depth migration is scanned, obtains really
Determine the coherent energy modal data of Depth Domain road collection.
(2) along stratum direction (that is, the direction in earth's surface sensing the earth's core) from shallow to deep by described Depth Domain
Road collection is in turn divided into multiple depth intervals adjacent one another are.Such as, described Depth Domain road collection can be divided
For two or three above depth intervals, it is preferable that be divided into more than four depth intervals.Such as, edge
Direction from shallow to deep, stratum, multiple depth intervals can be designated as the 1st depth intervals ..., i-th degree of depth
Interval ..., m depth intervals, m is positive integer, and m is not less than 3,1≤i < m.
Depth intervals can be divided into two and above depth intervals according to actual needs.Divide
During depth intervals, can be to start to certain degree of depth from earth's surface, as first depth intervals, so
The rear depth intervals that can be divided into one or more below.Can also certain below earth's surface deep
Degree starts, and this degree of depth depth below scope is marked off two or more depth intervals.It addition,
Demarcation interval also can according to according to the feature of depth migration road collection before prestack seismic data so that in a degree of depth district
The energy level of interior road collection coherent energy spectrum is close.Additionally, geological data prestack depth migration road collection
Feature can be different with the difference of geologic province, but for certain specific geologic province (such as, Sichuan
Chuan Zhong somewhere, basin, somewhere, basin, Xinjiang) for, its geological data prestack depth migration road collection
Feature determines substantially.
(3) to the coherent energy modal data of a certain depth intervals in multiple depth intervals and divide
To the plurality of depth intervals in deeper own than this depth intervals (the most aforementioned a certain depth intervals)
The coherent energy modal data of other depth intervals carries out n power computing and signed magnitude arithmetic(al) (n power respectively
Computing and signed magnitude arithmetic(al) without priority point), obtain the new coherent energy spectrum number of this depth intervals described
According to and the new coherent energy modal data of all depth intervals more deeper than this depth intervals, wherein,
All degree n are the positive number not less than 1, and interval being concerned with of different depth carrying out n power computing
Degree n corresponding to energy spectra data increases along with the intensification of the depth intervals degree of depth.Such as, to i-th
The coherent energy modal data of depth intervals carries out niPower computing and signed magnitude arithmetic(al), obtain the i-th degree of depth district
Between new coherent energy modal data；Same mode, the coherent energy spectrum number to i+1 depth intervals
According to carrying out ni+1Power computing and signed magnitude arithmetic(al), obtain the new coherent energy spectrum number of i+1 depth intervals
According to ..., the coherent energy modal data of m depth intervals is carried out nmPower computing and the fortune that takes absolute value
Calculate, obtain the corresponding coherent energy modal data of m depth intervals, wherein, ni、ni+1、……、nm
It is the positive number not less than 1, and ni<ni+1<……<nm.Thus, obtain the new of m depth intervals
Coherent energy modal data.
Additionally, in another exemplary embodiment of the present invention, improve seismic data pre-stack depth migration
The method of imaging road collection coherent energy spectral resolution, on the basis of above-mentioned exemplary embodiment, also includes
Include the relevant energy new at least one or two the above imaging road collection in all new coherent energy modal data
Amount modal data is normalized.It addition, the method for the present invention also can be to all new coherent energy spectrums
Each coherent energy modal data in data and do not carry out the coherent energy of depth intervals of n power computing
Modal data is normalized respectively.After being normalized, it is possible to improve seismic data further
The resolution of pre-stack depth migration imaging road collection coherent energy spectrum.
Additionally, another exemplary embodiment of the present invention is a kind of pre-stack depth migration velocity modeling method.
Described modeling method includes carrying out seismic data a series of conventional treatment to form the degree of depth-rate pattern, and
And described modeling method employing method as above improves seismic data pre-stack depth migration imaging road collection
The resolution of the different depth coherent energy spectrum obtained by scanning, thus form better the folding of reasonability
Front depth migration rate pattern.
The exemplary embodiment of the present invention is described in detail below in conjunction with concrete example.
By conventional method (such as, the mode of fitting of a polynomial), pre-stack depth migration imaging road collection is entered
Row scanning, obtains the coherent energy spectrum of different depth.
For different depth, use different n power computings, and the relevant energy that n power computing is obtained
Amount spectrum takes absolute value computing so that the resolution of coherent energy spectrum is improved.
Such as: the degree of depth of selected Depth Domain road collection is divided into 5 intervals (that is, m=5), Mei Gequ
Between the energy level of coherent energy spectrum close.The direction from shallow to deep along stratum, the degree of depth that division is obtained
Interval is designated as depth intervals one, depth intervals two, depth intervals three, depth intervals four and degree of depth district successively
Between five.
The computing mode that depth intervals one uses is:I.e. be equivalent to | x |, be same energy level with x.
Now n1=1.
The computing mode that depth intervals two uses is:Now n2=3/2.
The computing mode that depth intervals three uses is: x2.Now n3=2.
The computing mode that depth intervals four uses is:Now n4=5/2.
The computing mode that depth intervals five uses is: | x3|.Now n5=3.
Here enumerate one group of depth intervals energy spectra data simplified, the effect of this example is described.Such as,
The coherent energy modal data of depth intervals one is { 1,5,9,5,1}；The coherent energy spectrum of depth intervals three
Data are { 1,1.5,2,2.1,2.5,2.1,1.5,1}.
In depth intervals one, it is easy to pickup to ceiling capacity 9.And in depth intervals three, do not allow
It is easy to identify to ceiling capacity 2.5, because its neighbouring energy is close.After the method using this example processes,
That is, this coherent energy modal data is carried out quadratic power process, coherent energy modal data become 1,2.25,4,
4.41,6.25,4.41,2.25,1}, then be just more prone to identifies 6.25, its correspondence is to need
Road set attribute parameter to be picked up.
The coherent energy spectrum through processing, output generates new coherent energy spectrum file.
It addition, also can be in the depth intervals one, two, three, four or five of n power computing
Or coherent energy modal data multiple, even all interval is normalized.
In sum, the method for the present invention efficiently solves the coherent energy spectrum of Depth Domain road collection along with the degree of depth
Change and change, and when the degree of depth increases, the problem that the resolution of coherent energy spectrum is the lowest,
The coherent energy spectrum of the Depth Domain road collection with good discrimination power obtained, is particularly advantageous in improving difference deep
The resolution of the coherent energy spectrum of degree and the deeper Depth Domain road collection of the degree of depth.It addition, the method institute of the present invention
It is more reasonable that the coherent energy spectrum of the Depth Domain road collection with fine resolution formed can advantageously form
The degree of depth-rate pattern.
The method of the present invention has the effect that can significantly improve deep layer coherent energy spectral resolution, it is simple to hands
Move or automatic Picking.Additionally, the method for the present invention can be at pre-stack depth migration velocity modeling along layer phase
Dry energy spectrum obtains good application effect in calculating.
Although describe the present invention already in connection with exemplary embodiment above, but ordinary skill people
Member, can be to above-described embodiment it should be clear that in the case of without departing from spirit and scope by the claims
Carry out various amendment.
Claims (4)
1. the method improving seismic data pre-stack depth migration imaging road collection coherent energy spectral resolution,
It is characterized in that, said method comprising the steps of:
Seismic data pre-stack depth migration gained imaging road collection is scanned, obtains determining Depth Domain road collection
Coherent energy modal data；
Along stratum, described Depth Domain road collection is divided into multiple depth intervals by direction from shallow to deep；
Coherent energy modal data and ratio to a certain depth intervals in the plurality of depth intervals are described
The coherent energy modal data of the deeper all depth intervals of a certain depth intervals carry out respectively n power computing and
Signed magnitude arithmetic(al), obtain described a certain depth intervals new coherent energy modal data and than described certain
The new coherent energy modal data of the one deeper all depth intervals of depth intervals, wherein, all degree n
It is the positive number not less than 1, and carries out the coherent energy modal data in the different depth interval of n power computing
Corresponding degree n increases along with the intensification of the depth intervals degree of depth.
Raising seismic data pre-stack depth migration imaging road the most according to claim 1 collection coherent energy
The method of spectral resolution, described method also includes at least one in all new coherent energy modal data
Or the new coherent energy modal data of two or more is normalized.
Raising seismic data pre-stack depth migration imaging road the most according to claim 1 collection coherent energy
The method of spectral resolution, wherein, is divided into two or three above depth intervals by described Depth Domain road collection.
4. a pre-stack depth migration velocity modeling method, described modeling method includes carrying out seismic data
A series of process are to form the degree of depth-rate pattern, it is characterised in that described modeling method uses right such as to want
The method described in any one in 1 to 3 is asked to sweep to improve seismic data pre-stack depth migration imaging road collection
Retouch the resolution of obtained different depth coherent energy spectrum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610333110.9A CN106019367B (en) | 2016-05-19 | 2016-05-19 | Improve coherent energy spectrum resolution method and pre-stack depth migration velocity modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610333110.9A CN106019367B (en) | 2016-05-19 | 2016-05-19 | Improve coherent energy spectrum resolution method and pre-stack depth migration velocity modeling method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106019367A true CN106019367A (en) | 2016-10-12 |
| CN106019367B CN106019367B (en) | 2018-01-16 |
Family
ID=57098554
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610333110.9A Active CN106019367B (en) | 2016-05-19 | 2016-05-19 | Improve coherent energy spectrum resolution method and pre-stack depth migration velocity modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106019367B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976715A (en) * | 2017-11-16 | 2018-05-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Mountain complex structure residual static correction method and reflection in-phase axis dislocation processing method |
| CN113050101A (en) * | 2019-12-26 | 2021-06-29 | 宁波飞芯电子科技有限公司 | Coherent signal receiving device, method and coherent detection system |
| CN113156495A (en) * | 2020-01-07 | 2021-07-23 | 中国石油天然气集团有限公司 | Grid chromatography inversion reflection point determination method and device |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090319190A1 (en) * | 2008-06-19 | 2009-12-24 | Chevron U.S.A. Inc. | System and method for seismic trace analysis |
| CN102129083A (en) * | 2010-12-15 | 2011-07-20 | 中国石油集团川庆钻探工程有限公司 | Coherent energy spectrum acquisition method |
| CN103364828A (en) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | Third-generation coherence cube algorithm fast realization method based on recursion and power method |
-
2016
- 2016-05-19 CN CN201610333110.9A patent/CN106019367B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090319190A1 (en) * | 2008-06-19 | 2009-12-24 | Chevron U.S.A. Inc. | System and method for seismic trace analysis |
| CN102129083A (en) * | 2010-12-15 | 2011-07-20 | 中国石油集团川庆钻探工程有限公司 | Coherent energy spectrum acquisition method |
| CN103364828A (en) * | 2013-07-01 | 2013-10-23 | 西安交通大学 | Third-generation coherence cube algorithm fast realization method based on recursion and power method |
Non-Patent Citations (3)
| Title |
|---|
| MIRKO VAN DER BAAN,ET AL.: "深层地震反射数据相干能量识别和重建", 《石油物探译丛》 * |
| 杜世通: "地震属性分析", 《油气地球物理》 * |
| 许自强，等: "叠前深度偏移在深水崎岖海底地震资料的应用", 《工程地球物理学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976715A (en) * | 2017-11-16 | 2018-05-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Mountain complex structure residual static correction method and reflection in-phase axis dislocation processing method |
| CN113050101A (en) * | 2019-12-26 | 2021-06-29 | 宁波飞芯电子科技有限公司 | Coherent signal receiving device, method and coherent detection system |
| CN113156495A (en) * | 2020-01-07 | 2021-07-23 | 中国石油天然气集团有限公司 | Grid chromatography inversion reflection point determination method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106019367B (en) | 2018-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105974479B (en) | The chromatography 2D/3D anisotropy Depth Domain velocity modeling method of GPU space lattice body | |
| CN102901985B (en) | A kind of Depth Domain interval velocity modification method being applicable to relief surface | |
| CN103969683B (en) | Method for picking position faces in batched mode based on constraint in three-dimensional seismic interpretation | |
| CN103424777B (en) | A kind of method that improves seismic imaging resolution ratio | |
| CN104620132B (en) | For inverse time migration shoots the signal enhancing (DeSSeRT) that piles up using diversity | |
| CN106019367A (en) | Method for improving coherent energy spectrum resolution ratio and modeling method for increasing prestack depth migration speed | |
| CN106886047A (en) | A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio | |
| CN107219554A (en) | The automatic obtaining method of the Value of residual static correction of land seismic data | |
| CN109444959A (en) | Full range high-precision interval velocity field method for building up | |
| CN105093281A (en) | Earthquake multi-wave modeling method under inverse framework | |
| CN106405653A (en) | Unconformity stratum denudation amount recovery method | |
| CN106646613A (en) | Depth-domain multi-scale well control modeling and imaging combination processing method | |
| CN104155694B (en) | A kind of residual static corrections reflecting converted shear wave common geophone stack section | |
| CN104090299B (en) | Surface seismic data amplitude compensation method based on VSP primary waves | |
| CN107390266A (en) | Speed update method and pre-stack depth migration velocity modeling method based on angle gathers | |
| CN106443777B (en) | The acquisition methods and converted-wave velocity analysis method of converted wave anisotropic parameters spectrum | |
| CN102338885B (en) | Three-component VSP data first arrival time automatic pick method | |
| CN103984016B (en) | Converted wave anisotropy Amplitudeversusangle road collection abstracting method | |
| CN103728666A (en) | Method for determining replacement speed of static correction middle reference plane | |
| CN105259578B (en) | Formation velocity based on vertical seismic profile data determines method | |
| CN106125136B (en) | A kind of common point main road set creation method of relief surface | |
| CN106990434A (en) | Ellipse expansion converted wave imaging method and system | |
| CN106443829A (en) | Method and apparatus for constructing near-surface model | |
| CN111142159B (en) | Seismic velocity modeling method and system for middle-deep layer abnormal body | |
| CN106908838A (en) | The method for building target area stratum inelastic attenuation quality factor three-dimensional model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder |
Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Patentee after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Address before: No. 216, No. 216, Huayang Avenue, Huayang Town, Shuangliu County, Shuangliu County, Sichuan Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180404 Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201110 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Patentee before: BGP Inc., China National Petroleum Corp. |
|
| TR01 | Transfer of patent right |