CN109991659A - A kind of post-stack seismic data high frequency sequence method for automatic tracking - Google Patents
A kind of post-stack seismic data high frequency sequence method for automatic tracking Download PDFInfo
- Publication number
- CN109991659A CN109991659A CN201910270711.3A CN201910270711A CN109991659A CN 109991659 A CN109991659 A CN 109991659A CN 201910270711 A CN201910270711 A CN 201910270711A CN 109991659 A CN109991659 A CN 109991659A
- Authority
- CN
- China
- Prior art keywords
- high frequency
- post
- interface
- seismic data
- phase
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 239000002344 surface layer Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004062 sedimentation 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
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 discloses a kind of post-stack seismic data high frequency sequence method for automatic tracking comprising following steps: Step 1: denoising to post-stack seismic data body and calculating instantaneous phase body；Step 2: tracking third-band sequence interface or system tract interface；Step 3: using the third-band sequence interface or system tract interface tracked in step 2 as screen work, segmentation calculate in step 1 calculated instantaneous phase body phase unwrapping body, be then combined with；Step 4: contour surface is extracted on phase unwrapping body after merging, as high frequency sequence.Using method of the invention, can automatic tracing high frequency sequence, it is high-efficient, and there is accurate isochronism.
Description
Technical field
The present invention relates to petroleum exploration technical field more particularly to a kind of post-stack seismic data high frequency sequence automatic tracing sides
Method.
Background technique
In seismic data explanation, high frequency sequence refers to four, Pyatyi sequence.High frequency sequence whens accurately equal is for utilizing
When slice carry out sedimentation analysis and portray lithosomic body and be of great significance.
High frequency sequence explains to be mainly manual interpretation at present.The high frequency sequence in work area is established according to well-log information in work area,
Or high frequency sequence is determined by defining simple parallel relation.This completes high frequency sequence explanations.
The method that conventional highfrequency sequence is explained has subjectivity if then depending on the explanation experience of people by well-log information
Property, and take time and effort；If by define parallel relation determine, it is not necessarily consistent with seismic event therefore high-efficient but
Accuracy is low.
Summary of the invention
Place in view of the above-mentioned deficiencies in the prior art, the present invention provide a kind of post-stack seismic data high frequency sequence from
Dynamic method for tracing.The post-stack seismic data high frequency sequence method for automatic tracking the following steps are included:
Step 1: being denoised to post-stack seismic data body and calculating instantaneous phase body；
Step 2: tracking third-band sequence interface or system tract interface；
Step 3: segmentation calculates step using the third-band sequence interface or system tract interface tracked in step 2 as screen work
In one calculated instantaneous phase body phase unwrapping body, be then combined with；
Step 4: contour surface is extracted on phase unwrapping body after merging, as high frequency sequence.
Further, the step 2 specifically:
By the seed point of pickup, with linear interpolation method or method for automatic tracking, track out third-band sequence interface or
System tract interface.
Preferably, the step 3 specifically:
Objective function
In formula (1), φijIt is unknown quantity for the element of true phase Φ the i-th row jth column,WithRespectively indicate wink
The difference of Shi Xiangwei, i.e., are as follows:
Wherein, instantaneous phase ψ, size are M × N, and subscript x and y respectively indicate the difference on horizontal and vertical direction, and W is
Operator is wrapped up, is the modulo operation of 2 π；
φ is unknown quantity in the objective function, available by minimizing objective function；
φ is the optimal solution in global scope, meets lower relation of plane:
(φi+1,j-2φi,j+φi-1,j)+(φi,j-1-2φi,j+φi,j-1)=ρi,j (4)
Wherein:
Formula (5) is Poisson equation
Discrete form；The discrete Poisson equation shaped like formula (6) is solved, Neumann boundary condition is defined:
Under the boundary condition of above formula defined, instantaneous phase is carried out ρ is calculated according to formula (7)i,j, reciprocity later
Formula (4) both sides carry out discrete cosine transform (DCT) simultaneously, and the boundary condition of convolution (7), arrangement obtain:
In formula,WithRespectively φi,jAnd ρi,jDiscrete cosine transform, it is finally rightCarry out inverse discrete cosine transform
φ can be acquiredi,j, i.e., the phase be finally unfolded.
Further, in the step 3, phase unwrapping body is calculated using segmented mode, is tracked in step 2
Third-band sequence interface or system tract interface constrain phase unwrapping, and phase unwrapping body is made to have accurate isochronism.
Further, in the step 4, phase unwrapping body contour surface corresponds to seismic event isochronous surface, therefore can be with
As sequence whens waiting.
Post-stack seismic data high frequency sequence method for automatic tracking of the invention has the following beneficial effects:
Post-stack seismic data high frequency sequence method for automatic tracking provided by the present invention can be efficient compared to manual interpretation
Explain high frequency sequence in ground；Due to being based on earthquake instantaneous phase, the equivalent surface layer position of the phase body of expansion is extracted, is had same with earthquake
The consistent isochronism of phase axis.Therefore post-stack seismic data high frequency sequence method for automatic tracking of the invention is than traditional based on simple
The method that parallel relation is assumed is more accurate.
Detailed description of the invention
Fig. 1 is the flow diagram of post-stack seismic data high frequency sequence method for automatic tracking provided in an embodiment of the present invention；
Fig. 2 is the graphic software platform form at post-stack seismic data body (one section of display) and the third-band sequence interface explained
Figure；
Fig. 3 is that instantaneous phase body is segmented in Sequence Stratigraphic Formation, and each section is sought phase unwrapping body respectively, then will knot
Fruit merges the graphic software platform form figure of (one section of display)；
Fig. 4 is the contour surface extracted on phase unwrapping body (one section of display), i.e. the graphic software platform shape of high frequency sequence
Formula figure.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiments of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
Member's every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
The present embodiment provides a kind of post-stack seismic data high frequency sequence method for automatic tracking referring to FIG. 1 to FIG. 4, this is folded
Afterwards seismic data high frequency sequence method for automatic tracking as shown in Figure 1 the following steps are included:
S101, post-stack seismic data body is denoised and calculates instantaneous phase body；
S102, tracking third-band sequence interface or system tract interface；
S103, using the third-band sequence interface tracked in S102 or system tract interface as screen work, segmentation calculate S101 in institute
The phase unwrapping body of calculated instantaneous phase body, is then combined with；
Contour surface is extracted on S104, phase unwrapping body after merging, as high frequency sequence.
Further, S102 specifically:
By the seed point of pickup, with linear interpolation method or method for automatic tracking, track out third-band sequence interface or
System tract interface.
Preferably, S103 specifically:
Objective function
In formula (1), φijIt is unknown quantity for the element of true phase Φ the i-th row jth column,WithRespectively indicate wink
The difference of Shi Xiangwei, i.e., are as follows:
Wherein, instantaneous phase ψ, size are M × N, and subscript x and y respectively indicate the difference on horizontal and vertical direction, and W is
Operator is wrapped up, is the modulo operation of 2 π；
φ is unknown quantity in the objective function, available by minimizing objective function；
φ is the optimal solution in global scope, meets lower relation of plane:
(φi+1,j-2φi,j+φi-1,j)+(φi,j-1-2φi,j+φi,j-1)=ρi,j (4)
Wherein:
Formula (5) is Poisson equation
Discrete form；The discrete Poisson equation shaped like formula (6) is solved, Neumann boundary condition is defined:
Under the boundary condition of above formula defined, instantaneous phase is carried out ρ is calculated according to formula (7)i,j, reciprocity later
Formula (4) both sides carry out discrete cosine transform (DCT) simultaneously, and the boundary condition of convolution (7), arrangement obtain:
In formula,WithRespectively φi,jAnd ρi,jDiscrete cosine transform, it is finally rightCarry out inverse discrete cosine transform
φ can be acquiredi,j, i.e., the phase be finally unfolded.
Further, in the S103, phase unwrapping body is calculated using segmented mode, is with the three-level tracked in S102
Sequence interface or system tract interface constrain phase unwrapping, and phase unwrapping body is made to have accurate isochronism.In S104, phase unwrapping body
Contour surface corresponds to seismic event isochronous surface, therefore can be used as sequence whens waiting.
Post-stack seismic data high frequency sequence method for automatic tracking provided by the present embodiment can be high compared to manual interpretation
Explain high frequency sequence in effect ground；Due to being based on earthquake instantaneous phase, the equivalent surface layer position of the phase body of expansion is extracted, is had and earthquake
The consistent isochronism of lineups.Therefore post-stack seismic data high frequency sequence method for automatic tracking of the invention is than traditional based on letter
The method that single parallel relation is assumed is more accurate.
In addition, it should be noted that, the terms "include", "comprise" or its any other variant are intended in the present embodiment
Cover non-exclusive inclusion, so that including that process, method, article or the terminal device of a series of elements not only include
Those elements, but also including other elements that are not explicitly listed, or further include for this process, method, article or
The intrinsic element of person's terminal device.In the absence of more restrictions, the element limited by sentence "including a ...",
Be not precluded is including that there is also other identical elements in the process, method of the element, article or terminal device.
The above is only the preferred embodiment of the present invention, are not intended to restrict the invention, and those skilled in the art are come
It says, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any modification, equivalent
Replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (5)
1. a kind of post-stack seismic data high frequency sequence method for automatic tracking, which is characterized in that the described method comprises the following steps:
Step 1: being denoised to post-stack seismic data body and calculating instantaneous phase body；
Step 2: tracking third-band sequence interface or system tract interface；
Step 3: segmentation calculates in step 1 using the third-band sequence interface or system tract interface tracked in step 2 as screen work
Calculated instantaneous phase body phase unwrapping body, be then combined with；
Step 4: contour surface is extracted on phase unwrapping body after merging, as high frequency sequence.
2. post-stack seismic data high frequency sequence method for automatic tracking as described in claim 1, which is characterized in that the step 2
Specifically:
Third-band sequence interface or system are tracked out with linear interpolation method or method for automatic tracking by the seed point of pickup
Domain interface.
3. post-stack seismic data high frequency sequence method for automatic tracking as described in claim 1, which is characterized in that the step 3
Specifically:
Objective function
In formula (1), φijIt is unknown quantity for the element of true phase Φ the i-th row jth column,WithRespectively indicate instantaneous phase
The difference of position, i.e., are as follows:
Wherein, instantaneous phase ψ, size are M × N, and subscript x and y respectively indicate the difference on horizontal and vertical direction, and W is package
Operator is the modulo operation of 2 π；
φ is unknown quantity in the objective function, available by minimizing objective function；
φ is the optimal solution in global scope, meets lower relation of plane:
(φi+1,j-2φi,j+φi-1,j)+(φi,j-1-2φi,j+φi,j-1)=ρi,j (4)
Wherein:
Formula (5) is Poisson equation
Discrete form；The discrete Poisson equation shaped like formula (6) is solved, Neumann boundary condition is defined:
Under the boundary condition of above formula defined, instantaneous phase is carried out ρ is calculated according to formula (7)i,j, peer-to-peer (4) later
Both sides carry out discrete cosine transform (DCT) simultaneously, and the boundary condition of convolution (7), arrangement obtain:
In formula,WithRespectively φi,jAnd ρi,jDiscrete cosine transform, it is finally rightCarry out inverse discrete cosine transform
Acquire φi,j, i.e., the phase be finally unfolded.
4. post-stack seismic data high frequency sequence method for automatic tracking as claimed in claim 3, which is characterized in that the step 3
In, using segmented mode calculate phase unwrapping body, be with the third-band sequence interface or system tract interface tracked in step 2 about
Beam phase unwrapping makes phase unwrapping body have accurate isochronism.
5. post-stack seismic data high frequency sequence method for automatic tracking as described in claim 1, which is characterized in that the step 4
In, phase unwrapping body contour surface corresponds to seismic event isochronous surface, therefore can be used as sequence whens waiting.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910270711.3A CN109991659A (en) | 2019-04-04 | 2019-04-04 | A kind of post-stack seismic data high frequency sequence method for automatic tracking |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910270711.3A CN109991659A (en) | 2019-04-04 | 2019-04-04 | A kind of post-stack seismic data high frequency sequence method for automatic tracking |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109991659A true CN109991659A (en) | 2019-07-09 |
Family
ID=67132397
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910270711.3A Pending CN109991659A (en) | 2019-04-04 | 2019-04-04 | A kind of post-stack seismic data high frequency sequence method for automatic tracking |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109991659A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113945973A (en) * | 2020-07-17 | 2022-01-18 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6594585B1 (en) * | 1999-06-17 | 2003-07-15 | Bp Corporation North America, Inc. | Method of frequency domain seismic attribute generation |
| CN104597494A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method and device for analyzing seismic formation body |
| US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
-
2019
- 2019-04-04 CN CN201910270711.3A patent/CN109991659A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6594585B1 (en) * | 1999-06-17 | 2003-07-15 | Bp Corporation North America, Inc. | Method of frequency domain seismic attribute generation |
| CN104597494A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method and device for analyzing seismic formation body |
| US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
Non-Patent Citations (5)
| Title |
|---|
| MARCÍLIO CASTRO DE MATOS，等: "Stratigraphic discontinuities mapped through joint time-frequency seismic phase unwrapping", 《SEG》 * |
| 单蕊，等: "地震沉积学在川东北YB地区沉积相分析中的应用", 《煤田地质与勘探》 * |
| 徐长贵: "渤海海域低勘探程度区古近系岩性圈闭预测", 《中国博士学位论文全文数据库》 * |
| 胡春光: "利用相移显微干涉术和频闪成像技术研究MEMS离面运动", 《中国优秀硕士学位论文全文数据库》 * |
| 陈茂山: "地震地层体及其分析方法", 《石油地球物理勘探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113945973A (en) * | 2020-07-17 | 2022-01-18 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
| CN113945973B (en) * | 2020-07-17 | 2024-04-09 | 中国石油化工股份有限公司 | Reservoir characteristic analysis method, storage medium and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hale | A method for estimating apparent displacement vectors from time-lapse seismic images | |
| FR2974636B1 (en) | DEVICE AND METHOD FOR ADAPTIVE ANGULAR EXTRAPOLATION FOR THREE-DIMENSIONAL REVERSE MIGRATION ANGULAR GAPS IN TIME | |
| CN107817526B (en) | Prestack seismic gather segmented amplitude energy compensation method and system | |
| US10234581B2 (en) | System and method for high resolution seismic imaging | |
| WO2015160652A1 (en) | Generating subterranean imaging data based on vertical seismic profile data | |
| US9581709B2 (en) | Suppressing 4D-noise by weighted stacking of simultaneously acquired wave-fields | |
| EP2773984A1 (en) | Methods and devices for transformation of collected data for improved visualization capability | |
| CN107193045B (en) | A kind of seismic data processing technique and device | |
| SONG et al. | Research status and progress in multiple elimination | |
| KR20180067650A (en) | FWI model domain angular stacks with amplitude preservation | |
| WO2013152221A1 (en) | Converting a first acquired data subset to a second acquired data subset | |
| CN105093318B (en) | A kind of adaptive wave equation wave field extrapolation static correcting method | |
| CN109991659A (en) | A kind of post-stack seismic data high frequency sequence method for automatic tracking | |
| CN111694053A (en) | First arrival picking method and device | |
| Shragge et al. | Wave-equation migration from topography | |
| US9442207B1 (en) | System and method for computing residual moveout from seismic images | |
| CN105445795B (en) | The velocity analysis method and its system of tilted interface | |
| CN108802822B (en) | The direct prestack time migration method of guarantor's width and device in direction anisotropy medium | |
| Xiong et al. | Extracting a high-quality data space for stereo-tomography based on a 3D structure tensor algorithm and kinematic de-migration | |
| LIU et al. | Investigation on reverse-time migration of VSP data | |
| CN103852786B (en) | Reverse time migration imaging method and system applied to land seismic data | |
| LI et al. | Seismic blind deconvolution method based on generalized Gaussian distribution | |
| CN117908121A (en) | Seismic profile dip angle extraction method, system and electronic equipment | |
| Zinck et al. | Local dip transformation for fast seismic horizon reconstruction | |
| CN107238863B (en) | Method and device for judging authenticity of anticline structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190709 |
|
| RJ01 | Rejection of invention patent application after publication |