CN110261908A - Effective converted wave stacking image method - Google Patents
Effective converted wave stacking image method Download PDFInfo
- Publication number
- CN110261908A CN110261908A CN201910572701.5A CN201910572701A CN110261908A CN 110261908 A CN110261908 A CN 110261908A CN 201910572701 A CN201910572701 A CN 201910572701A CN 110261908 A CN110261908 A CN 110261908A
- Authority
- CN
- China
- Prior art keywords
- wave
- velocity ratio
- wave velocity
- converted
- converted wave
- 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 15
- 238000003384 imaging method Methods 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 claims 1
- 238000000605 extraction Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 13
- 230000008707 rearrangement Effects 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000010586 diagram 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
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
Description
技术领域technical field
本发明涉及油气勘探地震资料处理领域,特别是涉及到一种有效的转换波叠加成像方法。The invention relates to the field of oil and gas exploration seismic data processing, in particular to an effective converted wave superposition imaging method.
背景技术Background technique
利用纵横波速度比抽取共渐近点道集进行转换波的叠加是一种有效的转换波成像方法。一般的利用转换波共渐近点道集进行叠加成像是通过测试选取叠加效果最好的一个速度比进行转换波的叠加成像。It is an effective converted-wave imaging method to extract co-asymptotic point gathers by using the velocity ratio of P-to-S waves to superimpose converted waves. The common asymptotic point gathers of converted waves for stacking imaging is to select a velocity ratio with the best stacking effect through testing to perform stacking imaging of converted waves.
一般的做法往往是对于同一个纵横波速度比转换波地震数据中有的叠加段剖面效果好,有的转换波叠加段叠加效果不好,不能保证转换波的地震所有叠加段的叠加效果都好。The general practice is that for the same P-wave velocity ratio converted wave seismic data, some stacked sections have better stacking effects, and some converted-wave stacked sections have poor stacking effects, and it cannot be guaranteed that the stacking effects of all stacked sections of converted wave seismic data are good. .
利用某种方法能够使转换波所有地震叠加段的叠加效果都好是转换波叠加成像的必然需求。为此我们发明了一种有效的转换波叠加成像方法,解决了以上技术问题。It is an inevitable requirement for converted wave stacking imaging to use a certain method to make the stacking effect of all seismic stacking segments of converted wave good. For this reason, we invented an effective converted wave superposition imaging method, which solved the above technical problems.
发明内容Contents of the invention
本发明的目的是提供利用一种能够使转换波所有地震叠加段的叠加效果都好,能够实现转换波地震数据叠加成像的有效的转换波叠加成像方法。The purpose of the present invention is to provide an effective converted wave stacking imaging method which can make the stacking effect of all the seismic stacking sections of the converted wave good and realize the stacking imaging of converted wave seismic data.
本发明的目的可通过如下技术措施来实现:一种有效的转换波叠加成像方法,该一种有效的转换波叠加成像方法包括:步骤1:输入转换波地震数据;步骤2:确定纵横波速度比的范围;步骤3:确定纵横波速度比的间隔大小;步骤4:利用纵横波速度比的最大值,最小值和间隔计算纵横波速度比;步骤5:抽取转换波的共渐近点道集;步骤6:利用不同速度比的共渐近点道集重排得到新的共渐近点道集;步骤7:利用新的共渐近点道集进行叠加成像得到转换波的叠加剖面。The purpose of the present invention can be achieved by the following technical measures: an effective converted wave superposition imaging method, which comprises: step 1: input converted wave seismic data; step 2: determine the longitudinal and transverse wave velocity The range of the ratio; Step 3: Determine the interval size of the P-to-S wave velocity ratio; Step 4: Use the maximum value, minimum value and interval of the P-to-S wave velocity ratio to calculate the P-to-S wave velocity ratio; Step 5: Extract the co-asymptotic point channel of the converted wave Step 6: Use the co-asymptotic point gathers rearrangement with different velocity ratios to obtain a new co-asymptotic point gather; Step 7: Use the new co-asymptotic point gathers to perform stack imaging to obtain the stacked profile of converted waves.
本发明的目的还可通过如下技术措施来实现:The purpose of the present invention can also be achieved through the following technical measures:
在步骤1中,输入的原始转换波地震数据为有一定形态规律特征的地震信号。In step 1, the input original converted wave seismic data is a seismic signal with certain morphological characteristics.
在步骤2中,利用转换波地震数据测试纵横波速度比,确定纵横波速度比的范围。In step 2, the P-to-S wave velocity ratio is tested using the converted wave seismic data to determine the range of the P-to-S wave velocity ratio.
在步骤3中,根据确定的纵横波速度比范围,确定纵横波速度比的间隔。In step 3, the interval of the ratio of compressional and shear waves is determined according to the determined range of the ratio of compressional and shear waves.
在步骤4中,利用确定的纵横波速度比范围和纵横波速度比间隔计算不同的纵横波速度比是通过以下公式完成:In step 4, using the determined P-to-S wave velocity ratio range and P-to-S wave velocity ratio interval to calculate different P-to-S wave velocity ratios is completed by the following formula:
γ=γmin+Δγ和γ≤γmax计算不同的纵横波速度比,其中γ是需要计算的纵横波速度比,γmin是确定的最小纵横波速度比,Δγ是确定的纵横波速度比间隔,γmax是确定的最大纵横波速度比。γ=γ min +Δγ and γ≤γ max calculate different P-to-S wave velocity ratios, where γ is the P-to-S wave velocity ratio to be calculated, γ min is the determined minimum P-to-S wave velocity ratio, and Δγ is the determined P-to-S wave velocity ratio interval , γ max is the determined maximum P-to-S wave velocity ratio.
在步骤5中,利用计算的纵横波速度比抽取共渐近点道集是通过以下公式计算的:In step 5, the co-asymptotic point gather is extracted using the calculated P-to-S wave velocity ratio by the following formula:
其中xacp是共渐近点的坐标,γ是纵横波速度比,xoffset是炮检距。where x acp is the coordinate of the co-asymptotic point, γ is the ratio of P-to-S wave velocity, and x offset is the offset.
该有效的转换波叠加成像方法还包括,在步骤7之后,输出最终的转换波叠加地震数据。The effective converted wave stacking imaging method further includes, after step 7, outputting the final converted wave stacking seismic data.
本发明中的有效的转换波叠加成像方法,克服了现有转换波叠加方法对于同一个速度比可以使某一段转换波叠加剖面效果好,某一段转换波波叠加剖面效果不好,不能使所有的转换波叠加剖面效果都好的缺陷,利用不同速度比抽取的转换波共渐近点道集确定新的转换波共渐近点道集,实现转换波的叠加成像。The effective converted wave stacking imaging method in the present invention overcomes that the existing converted wave stacking method can make a certain section of converted wave stacking section effect good for the same speed ratio, but a certain section of converted wave stacking section has a bad effect and cannot make all The converted wave superimposition profile effect is good, and the converted wave common asymptotic point gathers extracted by different velocity ratios are used to determine a new converted wave common asymptotic point gather, so as to realize the superposition imaging of converted waves.
附图说明Description of drawings
图1为本发明的有效的转换波叠加成像方法的一具体实施例的流程图;Fig. 1 is the flowchart of a specific embodiment of the effective converted wave superposition imaging method of the present invention;
图2为本发明的一具体实施例中利用一个速度比抽取的转换波共渐近点道集的叠加剖面示意图;Fig. 2 is a superimposed cross-sectional schematic diagram of converted wave common asymptotic point gathers extracted by using a velocity ratio in a specific embodiment of the present invention;
图3为本发明的一具体实施例中利用不同速度比抽取转换波共渐近道集重排得到新的转换波共渐近点道集的叠加剖面示意图。Fig. 3 is a schematic cross-sectional view of a new converted-wave co-asymptotic point gather extracted by using different speed ratios to extract and rearrange converted-wave co-asymptotic point gathers in a specific embodiment of the present invention.
具体实施方式Detailed ways
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。In order to make the above and other objects, features and advantages of the present invention more comprehensible, the preferred embodiments are listed below and shown in the accompanying drawings in detail as follows.
本发明的一种有效的转换波叠加成像方法,能够使所有的转换波叠加剖面效果都好,实现转换波地震数据的叠加成像。其包括的主要步骤是:An effective converted wave stacking imaging method of the present invention can make all converted wave stacking sections have good effects, and realize stacking imaging of converted wave seismic data. The main steps it includes are:
步骤1,输入的原始转换波地震数据为有一定形态规律特征的地震信号;Step 1, the input original converted wave seismic data is a seismic signal with certain morphological characteristics;
步骤2,利用转换波地震数据测试纵横波速度比,确定纵横波速度比的范围;Step 2, using the converted wave seismic data to test the P-to-S wave velocity ratio, and determine the range of the P-to-S wave velocity ratio;
步骤3,根据确定的纵横波速度比范围,确定纵横波速度比的间隔;Step 3, according to the determined P-to-S wave velocity ratio range, determine the interval of P-to-S wave velocity ratio;
步骤4,利用确定的纵横波速度比范围和纵横波速度比间隔计算不同的纵横波速度比是通过以下公式完成:Step 4, using the determined P-to-S wave velocity ratio range and P-to-S wave velocity ratio interval to calculate different P-to-S wave velocity ratios is completed by the following formula:
γ=γmin+Δγ和γ≤γmax计算不同的纵横波速度比,其中γ是需要计算的纵横波速度比,γmin是确定的最小纵横波速度比,Δγ是确定的纵横波速度比间隔,γmax是确定的最大纵横波速度比;γ=γ min +Δγ and γ≤γ max calculate different P-to-S wave velocity ratios, where γ is the P-to-S wave velocity ratio to be calculated, γ min is the determined minimum P-to-S wave velocity ratio, and Δγ is the determined P-to-S wave velocity ratio interval , γ max is the determined maximum P-to-S wave velocity ratio;
步骤5,利用步骤4中计算的纵横波速度比抽取共渐近点道集是通过以下公式计算的:Step 5, using the P-to-S wave velocity ratio calculated in Step 4 to extract the common asymptotic point gather is calculated by the following formula:
其中xacp是共渐近点的坐标,γ是纵横波速度比,xoffset是炮检距;where x acp is the coordinate of the co-asymptotic point, γ is the ratio of P-to-S wave velocity, and x offset is the offset;
步骤6,利用不同速度比的共渐近点道集重排得到新的共渐近点道集;Step 6, using the co-asymptotic point gather rearrangement of different speed ratios to obtain a new co-asymptotic point gather;
步骤7,利用新的共渐近点道集进行叠加成像得到转换波的叠加剖面;Step 7, use the new co-asymptotic point gathers to perform stacking imaging to obtain the stacking section of the converted wave;
步骤8,输出最终的转换波叠加地震数据。Step 8, outputting the final converted wave stacked seismic data.
如图1所示,图1为本发明的一种有效的转换波叠加成像方法的流程图。As shown in FIG. 1 , FIG. 1 is a flow chart of an effective converted wave superposition imaging method of the present invention.
步骤101,输入需要处理的转换波地震数据;Step 101, input converted wave seismic data to be processed;
步骤102,利用转换波地震数据测试纵横波速度比,确定该转换波数据的纵横波速度比的范围;Step 102, using the converted wave seismic data to test the P-to-S wave velocity ratio, and determine the range of the P-to-S wave velocity ratio of the converted wave data;
步骤103,根据确定的纵横波速度比范围,确定纵横波速度比的间隔;Step 103, according to the determined range of the P-to-S wave velocity ratio, determine the interval of the P-to-S wave velocity ratio;
步骤104,利用确定的纵横波速度比范围和纵横波速度比间隔计算不同的纵横波速度比是通过以下公式完成:Step 104, using the determined P-to-S wave speed ratio range and P-to-S wave speed ratio interval to calculate different P-to-S wave speed ratios is completed by the following formula:
γ=γmin+Δγ和γ≤γmax计算不同的纵横波速度比,其中γ是需要计算的纵横波速度比,γmin是确定的最小纵横波速度比,Δγ是确定的纵横波速度比间隔,γmax是确定的最大纵横波速度比;γ=γ min +Δγ and γ≤γ max calculate different P-to-S wave velocity ratios, where γ is the P-to-S wave velocity ratio to be calculated, γ min is the determined minimum P-to-S wave velocity ratio, and Δγ is the determined P-to-S wave velocity ratio interval , γ max is the determined maximum P-to-S wave velocity ratio;
步骤105,利用步骤104中计算的速度比抽取公渐近点道集是通过以下公式计算的:Step 105, using the speed ratio calculated in step 104 to extract the common asymptotic point gather is calculated by the following formula:
其中xacp是共渐近点的坐标,γ是纵横波速度比,xoffset是炮检距;where x acp is the coordinate of the co-asymptotic point, γ is the ratio of P-to-S wave velocity, and x offset is the offset;
步骤106,利用不同速度比的共渐近点道集重排得到新的共渐近点道集;Step 106, using the rearrangement of co-asymptotic point gathers with different speed ratios to obtain a new co-asymptotic point gather;
步骤107,利用新的共渐近点道集进行叠加成像得到转换波的叠加剖面;Step 107, using the new co-asymptotic point gathers to perform superimposed imaging to obtain the superimposed profile of converted waves;
步骤108,输出最终的转换波叠加地震数据如图3所示。利用一个速度比抽取转换波共渐近点道集叠加成像的剖面如图2所示。Step 108, outputting the final converted wave stacked seismic data as shown in FIG. 3 . Figure 2 shows the section of superimposed imaging of converted wave co-asymptotic point gathers extracted using a velocity ratio.
本发明的有效的转换波叠加成像方法是正确的,一个速度比抽取的转换波共渐近点道集的叠加剖面没有由不同速度比抽取的转换波共渐近道集重排得到新的转换波共渐近点道集叠加剖面的效果好。本发明的实际资料转换波叠加剖面效果好。这是由于本发明利用不同速度比抽取的转换波共渐近点道集重排得到的新的转换波共渐近点道集能够使转换波所有叠加段叠加剖面叠加效果好,能够实现转换波地震数据叠加成像,在转换波叠加成像上优于常规转换波共渐近点道集叠加的方法。The effective converted-wave stacking imaging method of the present invention is correct, and the stacking profile of the converted-wave co-asymptotic point gathers extracted by a velocity ratio does not obtain a new conversion by rearranging the converted-wave co-asymptotic point gathers extracted by different velocity ratios The effect of wave co-asymptotic point gather stacking section is good. The actual data converted wave stacking profile of the present invention has good effect. This is because the new converted wave common asymptotic point gathers obtained by rearranging the converted wave common asymptotic point gathers extracted by different speed ratios can make the superimposition effect of all superposition sections of converted waves good, and can realize the conversion wave Seismic data stacking imaging is superior to conventional converted wave co-asymptotic point gather stacking method in converted wave stacking imaging.
Claims (7)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572701.5A CN110261908A (en) | 2019-06-27 | 2019-06-27 | Effective converted wave stacking image method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572701.5A CN110261908A (en) | 2019-06-27 | 2019-06-27 | Effective converted wave stacking image method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110261908A true CN110261908A (en) | 2019-09-20 |
Family
ID=67922721
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910572701.5A Pending CN110261908A (en) | 2019-06-27 | 2019-06-27 | Effective converted wave stacking image method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110261908A (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6128580A (en) * | 1998-04-17 | 2000-10-03 | Bp Amoco Corporation | Converted-wave processing in many-layered anisotropic media |
| CN101923177A (en) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | Velocity analysis method and device for multi-wave seismic data |
| CN103091709A (en) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | Method and device for acquiring longitudinal wave and converted wave seismic data time matching relation |
| US20130176823A1 (en) * | 2012-01-05 | 2013-07-11 | Cggveritas Services Sa | Simultaneous joint estimation of the p-p and p-s residual statics |
-
2019
- 2019-06-27 CN CN201910572701.5A patent/CN110261908A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6128580A (en) * | 1998-04-17 | 2000-10-03 | Bp Amoco Corporation | Converted-wave processing in many-layered anisotropic media |
| CN101923177A (en) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | Velocity analysis method and device for multi-wave seismic data |
| US20130176823A1 (en) * | 2012-01-05 | 2013-07-11 | Cggveritas Services Sa | Simultaneous joint estimation of the p-p and p-s residual statics |
| CN103091709A (en) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | Method and device for acquiring longitudinal wave and converted wave seismic data time matching relation |
Non-Patent Citations (2)
| Title |
|---|
| 李国发,等: "转换波地震资料处理的关键问题与解决方法", 《中国矿业大学学报》 * |
| 王高成: "P-SV转换波提取与成像", 《中国优秀硕士学位论文全文数据库》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107807390B (en) | Seismic data processing method and system | |
| CN106290002B (en) | The whole process detection method of rock type I crack propagation based on three-point bending test | |
| CN104932010B (en) | A kind of diffracted wave separation method based on the sparse Radon transformation of shortcut fringing | |
| CN105044769B (en) | The method for improving the resolution ratio of seismic signal | |
| CN103217710B (en) | Method and system for obtaining domain gathers with common offset and common azimuth | |
| CN104391327B (en) | A kind of marine inclined shaft crosshole seismic prestack reverse-time depth migration imaging method | |
| CN101776768A (en) | Anisotropic speed analysis and dynamic correction method | |
| CN102721977A (en) | High-precision converted wave anisotropy stacking velocity analysis method | |
| CN111239832B (en) | 3D3C-VSP imaging processing method, device and equipment | |
| CN103744116A (en) | Prestack gather all-time-domain phase-equalization correction method | |
| CN110261908A (en) | Effective converted wave stacking image method | |
| JP2017209288A5 (en) | ||
| CN106908836A (en) | Method and system for collecting and suppressing footprints | |
| CN104914483B (en) | A kind of South Pole ground wire quality evaluating method that data are observed based on multi-source | |
| CN103984016A (en) | Converted wave anisotropy amplitude various angle gather extracting method | |
| CN106680875A (en) | Analysis method and apparatus of wide-azimuth seismic migration speed | |
| CN104112273A (en) | Image quality evaluation method based on log Gabor filter phase similarity | |
| Zhang et al. | Automatic salt dome detection using U-net | |
| CN105223630B (en) | Omnibearing observation systematic parameter Demonstration Method based on geological model | |
| CN106569278B (en) | A kind of similar coherent velocity of multiple tracks composes computational methods | |
| CN106772596A (en) | A kind of method and device for determining pre-stack time migration velocity field | |
| CN112925020B (en) | A Parametric Array Shallow Profile Intelligent Energy Compensation Method | |
| CN105487111B (en) | A kind of automatic position angle gathers partial stack method | |
| CN105259573A (en) | A method and device for acquiring seismic sediment body information | |
| CN104483701B (en) | Anisotropy characterization method and anisotropy characterization system based on all-direction common imaging point gather |
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: 20190920 |
|
| RJ01 | Rejection of invention patent application after publication |