CN102890289A - 一种高精度地震波层速度反演方法 - Google Patents
一种高精度地震波层速度反演方法 Download PDFInfo
- Publication number
- CN102890289A CN102890289A CN2011102035891A CN201110203589A CN102890289A CN 102890289 A CN102890289 A CN 102890289A CN 2011102035891 A CN2011102035891 A CN 2011102035891A CN 201110203589 A CN201110203589 A CN 201110203589A CN 102890289 A CN102890289 A CN 102890289A
- Authority
- CN
- China
- Prior art keywords
- fourier transform
- spectrum
- arrival
- phase
- frequency
- 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
本发明是勘探数据处理的高精度的地震波层速度反演方法,井下接收数据变换得到振幅谱,将高频振幅谱清零再反傅里叶变换,保留低频信号记录,拾取每道初至,在深度方向进行傅里叶变换得到F-K谱,将上行波对应的F-K谱清零,在波数方向和频率方向反傅里叶变换得到下行波,从每道的下行波中沿初至往后开一时窗,傅里叶变换得到每个频率的相位谱,求取相邻道中相同频率的相位谱之差,进行最小二乘拟合得到相邻道初至差,相邻道的深度差除以初至差得到层速度。本发明具有很强抗随即干扰能力,对激发子波和底层的吸收衰减情况适应能力强,不仅算法简单反演精度高。
Description
技术领域
本发明涉及地震勘探数据处理技术,本发明目的是提供一种利用地震数据的相位谱属性,且具有很强抗噪性的高精度的地震波层速度反演方法。
背景技术
在地震勘探中,零偏垂直地震剖面(VSP)由于炮点离井口很近,可认为不同深度接收到的下行直达波具有相同的传播路径,能直接利用不同深度地震记录的初至时间反演地震波的层速度,为地面地震处理以及解释提供有力手段。而目前的VSP资料处理中,都是先对地震记录在时间域进行重采样,提高采样率,再进行初至拾取。其具体做法是利用该道地震记录中相邻几个样点的值内插出更多的样点值,比如将1毫秒采样率的信号内插成0.0001毫秒采样率的信号,然后再进行初至拾取,利用获得的初至,计算相邻道之间的时差,再用相邻道的深度差除以时间差得到层速度。由于这种提高采样率的方法仅仅只是一种数学算法,其内插出的信号并不是真实信号,不能提高原始数据的分辨率,同时内插结果受内插方法和资料信噪比的影响很大,无法真正起稿初至拾取的精度,导致反演的层速度不够精确,制约了VSP资料的应用价值。
发明内容
本发明目的是提供一种利用地震数据的相位谱属性,且具有很强抗噪性的高精度地震波层速度反演方法。
本发明具体步骤包括:
1)地面震源激发,井下检波器接收得到VSP数据;
2)对地震记录进行傅里叶变换得到振幅谱,将高频部分的振幅谱清零,然后反傅里叶变换,得到只保留低频信号的地震记录,然后拾取每道地震记录的初至;
3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱,然后在深度方向再进行一次傅里叶变换得到F-K谱,并将上行波对应的F-K谱清零,然后分别在波数方向和频率方向进行反傅里叶变换得到下行波;
4)从每道的下行波中沿初至往后开一时窗,对该时窗内的信号进行傅里叶变换,得到每个频率的相位谱;
5)重复步骤4),得到所有道下行波每个频率的相位谱;
6)求取相邻道中相同频率的相位谱之差;
7)对步骤6)中的相位谱之差进行最小二乘拟合,得到相邻道初至差;
8)相邻道的深度差除以相邻道初至差得到高精度的层速度。
步骤6)所述的相邻道相位谱之差计算如下式:
则有:θi+1(ω)=2πfΔti+1
式中,si(t-Δti+1)为深度为hi的检波器得到的下行波为si(t),深度为hi+1的检波器得到的下行波;
Si+1(ω)为si+1(t)的傅里叶变换;
θi+1(ω)为si+1(t)的傅里叶变换的相位谱
Δti+1为第i+1道和第i道初至差;
j为复数单位。
本发明具有很强抗随即干扰能力,对激发子波和底层的吸收衰减情况适应能力强,不仅算法简单还能大大节约工作量,且反演的层速度精度很高。
附图说明
图1为本发明截取的时窗内信号对应的频谱;
图2为相邻道中相同频率的相位谱之差;
图3为反演层速度对比。
具体实施方式
以下结合附图详细说明本发明。
本发明是一种利用地震数据的相位谱属性,且具有很强抗噪性的高精度层速度反演方法。
本发明具体实现步骤如下:
1)地面震源激发,井下检波器接收得到VSP数据;
2)对地震记录进行傅里叶变换得到振幅谱,将高频部分的振幅谱清零,然后反傅里叶变换,得到只保留低频信号的地震记录,然后拾取每道地震记录的初至;
3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱,然后在深度方向再进行一次傅里叶变换得到F-K谱,并将上行波对应的F-K谱清零,然后分别在波数方向和频率方向进行反傅里叶变换得到下行波;
4)从每道的下行波中沿初至往后开一时窗,对该时窗内的信号进行傅里叶变换,得到每个频率的相位谱;图1截取的该时窗内信号对应的频谱;
5)重复步骤4),得到所有道下行波每个频率的相位谱;
6)求取相邻道中相同频率的相位谱之差;图2为相邻道中相同频率的相位谱之差;
设深度为hi的检波器得到的下行波为si(t),深度为hi+1的检波器得到的下行波为si(t-Δti+1),根据傅里叶变换的性质:信号时间域的时移对应其傅里叶变换频率域的相移;有:
则有:θi+1(ω)=2πfΔti+1
式中,Si+1(ω)为si+1(t)的傅里叶变换;
θi+1(ω)为si+1(t)的傅里叶变换的相位谱
Δti+1为第i+1道和第i道初至差;
j为复数单位。
(红色为加入部分)
7)对步骤6)中的相位谱之差进行最小二乘拟合,得到相邻道初至差;
8)相邻道的深度差除以相邻道初至差得到高精度的层速度。图3是反演层速度对比;蓝线是常规的时间域初至拾取反演得到的层速度,粉红线为本发明反演的层速度,其明显消除了人为拾取造成的异常点。
Claims (2)
1.一种高精度地震波层速度反演方法,特点是采用以下步骤:
1)地面震源激发,井下检波器接收得到VSP数据;
2)对地震记录进行傅里叶变换得到振幅谱,将高频部分的振幅谱清零,然后反傅里叶变换,得到只保留低频信号的地震记录,然后拾取每道地震记录的初至;
3)对步骤2)中的地震记录进行傅里叶变换得到振幅谱,然后在深度方向再进行一次傅里叶变换得到F-K谱,并将上行波对应的F-K谱清零,然后分别在波数方向和频率方向进行反傅里叶变换得到下行波;
4)从每道的下行波中沿初至往后开一时窗,对该时窗内的信号进行傅里叶变换,得到每个频率的相位谱;
5)重复步骤4),得到所有道下行波每个频率的相位谱;
6)求取相邻道中相同频率的相位谱之差;
7)对步骤6)中的相位谱之差进行最小二乘拟合,得到相邻道初至差;
8)相邻道的深度差除以相邻道初至差得到高精度的层速度。
2.根据权利要求1所述的方法,特点是步骤6)所述的相邻道相位谱之差计算如下式:
则有:θi+1(ω)=2πfΔti+1
式中,si(t-Δti+1)为深度为hi的检波器得到的下行波为si(t),深度为hi+1的检波器得到的下行波;
Si+1(ω)为si+1(t)的傅里叶变换;
θi+1(ω)为si+1(t)的傅里叶变换的相位谱
Δti+1为第i+1道和第i道初至差;
j为复数单位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110203589.1A CN102890289B (zh) | 2011-07-20 | 2011-07-20 | 一种高精度地震波层速度反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110203589.1A CN102890289B (zh) | 2011-07-20 | 2011-07-20 | 一种高精度地震波层速度反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102890289A true CN102890289A (zh) | 2013-01-23 |
CN102890289B CN102890289B (zh) | 2015-05-27 |
Family
ID=47533844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110203589.1A Active CN102890289B (zh) | 2011-07-20 | 2011-07-20 | 一种高精度地震波层速度反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102890289B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080225641A1 (en) * | 2003-08-15 | 2008-09-18 | Westerngeco, L.L.C. | Multi-Component Seismic Data Processing |
CN101363916A (zh) * | 2007-08-10 | 2009-02-11 | 中国石油集团东方地球物理勘探有限责任公司 | 一种高精度折射静校正数据反演方法 |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
-
2011
- 2011-07-20 CN CN201110203589.1A patent/CN102890289B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080225641A1 (en) * | 2003-08-15 | 2008-09-18 | Westerngeco, L.L.C. | Multi-Component Seismic Data Processing |
CN101363916A (zh) * | 2007-08-10 | 2009-02-11 | 中国石油集团东方地球物理勘探有限责任公司 | 一种高精度折射静校正数据反演方法 |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102890289B (zh) | 2015-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103376464B (zh) | 一种地层品质因子反演方法 | |
Hayashi et al. | CMP cross-correlation analysis of multi-channel surface-wave data | |
Le Feuvre et al. | Passive multi-channel analysis of surface waves with cross-correlations and beamforming. Application to a sea dike | |
Roberts et al. | Resolving a velocity inversion at the geotechnical scale using the microtremor (passive seismic) survey method | |
CN101598809A (zh) | 一种自适应消除线性规则噪声以及多次波干扰的方法 | |
WO2016187252A1 (en) | Surface wave tomography using sparse data acquisition | |
CN104533396A (zh) | 一种远探测声波的处理方法 | |
Hayashi et al. | CMP spatial autocorrelation analysis of multichannel passive surface-wave data | |
CN106526678A (zh) | 一种反射声波测井的波场分离方法及装置 | |
EA004486B1 (ru) | Способ ослабления поверхностной волны | |
Zhang et al. | Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise | |
Cheng et al. | Testing of a permanent orbital surface source and distributed acoustic sensing for monitoring of unconventional reservoirs: Preliminary results from the Eagle Ford Shale | |
CN102230973A (zh) | 一种三维分步傅立叶粘滞声波深度偏移方法 | |
CN109581481B (zh) | 一种便携式高频可控震源地震信号谐波干扰消除方法 | |
Baradello et al. | Fast method to transform chirp envelope data into pseudo-seismic data | |
CN102323618A (zh) | 基于分数阶傅里叶变换的相干噪声抑制方法 | |
CN108121010A (zh) | 基于孔巷联合的井下独头巷道槽波超前探测方法及系统 | |
Liberty | Hammer seismic reflection imaging in an urban environment | |
JP5517258B2 (ja) | 鉛直アレイ地震計を利用したq値測定方法 | |
Zhang et al. | Interval Q inversion based on zero-offset VSP data and applications | |
CN102890288B (zh) | 一种地震波层速度反演方法 | |
CN114415234B (zh) | 基于主动源面波频散和h/v确定浅地表横波速度的方法 | |
CN102890289B (zh) | 一种高精度地震波层速度反演方法 | |
US6263285B1 (en) | Amplitude spectra estimation | |
Galiana-Merino et al. | Seismic noise array measurements using broadband stations and vertical geophones: preliminary outcomes for the suitability on fk analysis |
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 |