CN101046516A - 利用垂直地震剖面和微测井进行地震信号补偿方法 - Google Patents
利用垂直地震剖面和微测井进行地震信号补偿方法 Download PDFInfo
- Publication number
- CN101046516A CN101046516A CN 200710027732 CN200710027732A CN101046516A CN 101046516 A CN101046516 A CN 101046516A CN 200710027732 CN200710027732 CN 200710027732 CN 200710027732 A CN200710027732 A CN 200710027732A CN 101046516 A CN101046516 A CN 101046516A
- Authority
- CN
- China
- Prior art keywords
- seismic
- signal
- seismic signal
- vsp
- micro logging
- 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
本发明提供了一种利用垂直地震剖面(VSP)和微测井对低、降速带近地表地震信号进行恢复的方法,该方法把近地表低、降速带看成一个信号传输***,通过双井微测井信号可获得地震信号在不同深度或者不同时间在低、降速带内进行传输时的传输函数,从而获得对地表接收到的衰减了的地震信号进行恢复的补偿滤波器函数;同时,通过不同地段接收的不同深度的VSP信号可获得并推广得到低、降速带不同深度或者不同时间的传输函数和补偿滤波器函数。由此,恢复出低、降速带不同地段接收到的地震信号,达到了提高地震信号分辨率目的,为岩性地震勘探提供了有利的工具。该法对于所有地表为低、降速带的地区均适用。
Description
技术领域:
本发明涉及一种利用垂直地震剖面(VSP)和微测井进行地震信号补偿方法,特别是利用VSP和微测井等信息分析近地表低、降速带地震信号的衰减规律,并对该区信号进行补偿以提高地震信号分辨率的有效方法。
背景技术:
石油地震勘探利用人工方法产生地震波,通过研究地震波在地下岩层中的传播情况,勘察地下地质构造和岩层变化,寻找油气田。
随着地震勘探不断发展,地震资料采集逐步向地表、地下地质情况异常复杂的地区延伸,由于地震地质条件的影响以及勘探技术水平的限制,像黄土塬、沙漠等地区表层为巨厚的低、降速带所覆盖,地震波在低、降速带内传播时会发生吸收、散射等衰减,地震子波不断变形增长,尤其是地震子波在近地表传播时地震子波衰减及散射非常大,近地表每米的衰减量是高速层每米的500~1000多倍。近地表高程不同、地层岩性不同,地震波的衰减规律则不相同。
由于近地表衰减严重,导致接收到的地震资料干扰严重、信噪比低、资料品质较差,用该资料难以获得精确的地下构造形态。为此,须对地震信号进行补偿以恢复地震信号。
目前对地震信号补偿的技术较多。文章“董世学等,检波器—地表耦合***对地震记录的影响,地球物理学报,2001,44(增刊):161~169”中利用特殊耦合检波器接收地震信号,拓宽了地震信号频带,提高了信噪比和分辨率。文章“姚振兴等,用于深度域地震剖面衰减与频散补偿的反Q滤波方法,地球物理学报,2003,46(2):229~233。”研究基于Futterman模型的反Q滤波方法,补偿了地震信号的频散与振幅。文章“Satinder C et al.High-frequencyrestoration of surface seismic data,The Leading Edge,2003,22(8):730~738.”利用VSP数据求取不同深度地层的反滤波器补偿高频成分。微VSP采用“井口激发、井中接收”的方式,初至波信噪比相对较高。国家发明专利(申请号:200510049341.9)“利用特殊耦合检波器进行地震数据接收的最佳耦合与匹配方法”和文章“石战结等,沙漠地区地震检波器耦合的高频信号匹配滤波技术”提出了一种利用特殊耦合检波实现与地表最佳耦合匹配的方法,该法使在沙漠地区的地震数据无论是在叠前还是在叠后其能量和频率都得到加强,提高了地震信号的分辨率。国家发明专利(专利公开号:CN1673775A)“利用地震微测井进行地震信号高频补偿方法”和文章“田钢等,利用微测井资料补偿地震数据的高频成分,石油地球物理勘探,2005,40(5):546~549。”提出利用微测井直达波信号求出不同厚度沙漠地区沙层的吸收补偿反滤波器来补偿地震数据的高频成分;这种微测井采用“井中激发、地表接收”,在地表接收时,可采用普通检波器,只要将检波器尾锥***硬土即可保证检波器与介质的紧密耦合,极大地减小了检波器接收时的能量损耗,但是,由于井中介质为空气,在黄土塬地区井壁内介质的岩性一般为疏松的黄土,故空气与黄土的波阻抗相差很大,有效波入射到井壁时反射波极强,投射波极弱,透射波从井壁内侧到达地表检波点处还会有一定的衰减,这样到达检波点处的有效波就更弱了。一般而言,微测井与VSP相比较而言,微测井的初至波比VSP的初至波干扰更严重。
发明内容:
本发明的目的是:针对现有技术的不足,提供一种利用VSP和微测井进行地震信号补偿的方法。
为实现上述目的,本发明的技术方案是:由于VSP微测井/双井微测井横穿遭受风化作用的低、降速带,在VSP微测井信号/双井微测井信号中比地表接收到的地震信号将包含更多的高频信号。通过分析在双井微测井中不同深度的地震信号衰减情况,可以以此为依据对地表接收的地震信号进行恢复。
下面介绍利用微测井信号对地震信号进行补偿的原理:
步骤1:在钻VSP微测井单深井(可用生产井代替)和双井微测井(在构造复杂地带),通过地震震源在地表激发,在VSP微测井或双井微测井深度分别为d1米、d2米、.......、dN米(N为整数,N≥3)处获得地震接收信号;对地震信号按照下行波、上行波、PS波等进行划分,假设xi(t)表示下行波在深度为di米时所接收到的地震信号的振幅,其中i=1、2、....、N。此时,xi(t)即为该深度时所接收到的直达波。设xN(t)表示下行波在低、降速带以下,深度为dN米时所接收到的高速层地震信号的振幅。
步骤2:用垂直地震剖面方法分离掉高速层中xi(t)的虚反射波,用高速层的速度将...、xn-2(t)、xn-1(t)、xn+1(t)、xn+2(t)、...静校正到与xn(t)地震记录相同,然后将...、xn-2(t)、xn-1(t)、xn(t)、xn+1(t)、xn+2(t)、...叠加提高直达波xn(t)的信噪比;此处对井中地震记录进行静校正是指将x0(t)与xN(t)之间的初至时间差去掉。
步骤3:建立井中微测井地震在近地表传播的数学物理方程,计算出补偿滤波器函数。
把近地表低、降速带看成一个信号传输***,信号从深度di米处传输到深度dj处,其***函数记为hj(t)(j=1、2、....、N),则***函数hj(t)也表示地震信号从深度di米到dj米处的衰减程度;那么,对于深度为dj米处所接收到的信号xj(t)可表示为:
xj(t)=xi(t)*hj(t) i=1,.......,N j=1,.......,N i≠j
(1)
对(1)式进行傅氏变换有:
Xj(f)=Xi(f)×Hj(f) (2)
对于***函数的逆算子
称为补偿滤波器函数,该参数可用于补偿信号的衰减。有:
xi(t)=xj(t)*gj(t) (3)
对(3)式进行傅氏变换有:
Xi(f)=Xj(f)Gj(f) (4)
因此,可得:
步骤4:用补偿滤波器函数和微测井周围的地面地震记录作叠积运算补偿地震信号。
yj(t)=wj(t)*gj(t) (7)
Yj(f)=Wj(f)Gj(f) (8)
其中,wj(t)、Wj(f)分别是地面地震记录的时域和频域信号,yj(f)、Yj(f)为近地表衰减的地震波经过补偿后获得的时域和频域信号。注意,对于补偿滤波器函数的选择一定要根据该地面地震记录观测点低、降速带的厚度dj确定。微测井周围是指微测井附近几百米到几十公里的范围内,一般是指与微测井处地面高程和岩性变化都不大的区域。
步骤5:对求得的地震记录进行信号处理及解释。
本发明就是利用地震波在近地表中传播衰减的规律,对地面接收到的衰减的地震波进行补偿,恢复到相当于高速层中小药量激发和低速层之下的高速层中接收的地震波,从而实现提高地震分辨率的目的。
具体实施方式:
把低、降速带看成一个信号传输***,该***的传输函数反映了不同深度、不同岩性地震信号的衰减情况。本发明对不同深度得到的VSP或者双井微测井记录进行保留初至波的编辑处理,得到不同深度地震信号频率和幅度的变化情况;其初至波的幅度和波长的变化可以用来估计在该低、降速带中在该深度时频率成分的变化情况,获得***传输函数。由此,对于不同深度初至波的滤波器响应,可根据时域信号或者频域信号计算出补偿滤波器函数,由此可补偿地面接收到的地震信号的衰减,实现对低降速带地震信号的恢复。
附图说明:
图1是说明地震微VSP测井示意图。VSP微测井方法原理与常规井中激发微测井原理一样,都是利用直达波在地层中的传播规律,通过改变检波器观测深度,获取不同的初至时间,根据深度换算成垂直时间,拟合时深曲线得到分层的速度和厚度。VSP微测井工艺,可以在地震生产激发井中进行而不影响生产井使用,可以大面积应用和推广。
图2是双井微测井原理图。双井微测井就是在一定的距离内钻2口打穿低降速带的井,2井间的距离一般为4m,如果表层低降速带内含有流砂层,则应适当地加大2井之间的距离,以保证2井间的原始地层不会被破坏。一口井内布设激发点,另一口井布设检波器。另外,在激发井井口附近布设地面检波器,地面和地下接收用的检波器型号相同,用10Hz或40Hz的单只检波器均可以。地面接收一般采用6~23道,每道一个检波器,检波器呈扇形或直角形,井检距一般采用1~4m。激发***由电缆、***及引线构成,把***固定在专用微测井电缆上的指定位置,一次性下入井中,山深至浅逐点激发。井中激发点的间距要考虑到介质的速度、厚度、检波器的灵敏度、仪器的采样间隔和动态范围等,激发点的间距过大会造成控制点数减少甚至丢层,反之,冗余的采样会造成浪费。接收井井下检波器深度和最深的激发点从相同。地面检波器采集的信号主要是地震波的初至,挑选初至无干扰且波形起跳干脆的道用于资料处理。
图3为本发明专利方法产生特定地段低、降速带地震信号衰减补偿滤波器的框图。
图4为本发明专利方法获得一般地段低、降速带地震信号衰减补偿的框图。通过双井微测井获得特定地段低、降速带不同深度地震信号衰减规律,利用生产炮井做VSP微测井获得整个工区不同低、降速带地震信号的衰减规律,由于生产炮井一般比较浅,对超过激发井深处的地震信号衰减规律则通过双井微测井获得的规律与之比较获得,因此,将双井微测井信号和VSP信号联合使用可获得整个低、降速带不同地段、不同深度的地震信号衰减规律和地震信号补偿滤波器,最终可以恢复出没有经过低、降速带衰减的地震信号,该地震信号无疑有更高的频率和信噪比,以此恢复后的地震信号进行地震信号处理可提高地震信号的分辨率,为岩性地震勘探油气提供有利工具。
Claims (5)
1、本发明涉及一种利用垂直地震剖面和微测井进行地震信号补偿方法,其特征在于包括如下步骤:
步骤1:利用双井微测井曲线获得特定地段的地震信号下行波;
步骤2:用垂直地震剖面方法分离掉高速层中xi(t)的虚反射波,用高速层的速度将…、xn-2(t)、xn-1(t)、xn+1(t)、xn+2(t)、…静校正到与xn(t)地震记录相同,然后将…、xn-2(t)、xn-1(t)、xn(t)、xn+1(t)、xn+2(t)、…叠加提高直达波xn(t)的信噪比;此处对井中地震记录进行静校正是指将x0(t)与xN(t)之间的初至时间差去掉。
步骤3:计算特定地段的地震信号初至波在不同深度时幅度和频率的变化情况,建立井中微测井地震在近地表传播的数学物理方程,计算出补偿滤波器函数。
步骤4:利用VSP信号获得不同地段地震信号的初至波;
步骤5:由VSP信号的初至波确定不同地段的地震信号补偿滤波器函数;
步骤6:将不同地段的地震信号衰减恢复滤波器函数与特定地段地震信号衰减恢复滤波器函数进行比较,获得整个观察区域的地震信号补偿滤波器函数;
步骤7:用补偿滤波器函数和所有观察区的地面地震记录作叠积运算补偿地震信号,消除低、降速带地震信号衰减。
2、如权利要求1所说,一种利用垂直地震剖面(VSP)和微测井进行地震信号补偿方法,其特征在于对于地质情况比较单一的地区直接采用步骤4、5、7就可以对低、降速带地震信号进行补偿。
3、如权利要求1和权利要求2所说,一种利用垂直地震剖面(VSP)和微测井对地震信号进行补偿的方法,其特征在于:该方法可用于近地表为低、降速带的所有地区。
4、如权利要求1和权利要求2所说,一种利用垂直地震剖面(VSP)和微测井对地震信号进行补偿的方法,其特征在于:该方法中步骤1和步骤4所采用的地震信号可以是单井微测井中获得的近地表垂直地震记录,也可以是双井微测井获得的近地表垂直地震记录。
5、如权利要求1和权利要求2所说,一种利用垂直地震剖面(VSP)和微测井对地震信号进行补偿的方法,其特征在于:该方法中步骤7采用的地面地震记录可以是地面地震检波点地震记录,也可以是检波点地震记录静校正结果,还可以是地震水平叠加后的地震剖面记录,仍可以是地震偏移后的地震剖面记录。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200710027732XA CN100552472C (zh) | 2007-04-22 | 2007-04-22 | 利用垂直地震剖面和微测井进行地震信号补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200710027732XA CN100552472C (zh) | 2007-04-22 | 2007-04-22 | 利用垂直地震剖面和微测井进行地震信号补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101046516A true CN101046516A (zh) | 2007-10-03 |
CN100552472C CN100552472C (zh) | 2009-10-21 |
Family
ID=38771277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200710027732XA Expired - Fee Related CN100552472C (zh) | 2007-04-22 | 2007-04-22 | 利用垂直地震剖面和微测井进行地震信号补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100552472C (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102109612B (zh) * | 2009-12-23 | 2012-10-17 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
CN103210323A (zh) * | 2011-03-31 | 2013-07-17 | 雪佛龙美国公司 | 处理地震数据的***和方法 |
WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN112147684A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
-
2007
- 2007-04-22 CN CNB200710027732XA patent/CN100552472C/zh not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102109612B (zh) * | 2009-12-23 | 2012-10-17 | 中国石油天然气集团公司 | 一种地震波吸收衰减补偿方法 |
CN103210323A (zh) * | 2011-03-31 | 2013-07-17 | 雪佛龙美国公司 | 处理地震数据的***和方法 |
CN103210323B (zh) * | 2011-03-31 | 2016-08-10 | 雪佛龙美国公司 | 处理地震数据的***和方法 |
WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
CN112147684A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
CN112147684B (zh) * | 2019-06-28 | 2023-12-26 | 中国石油天然气股份有限公司 | 表征同沉积断层活动强度的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN100552472C (zh) | 2009-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
Mi et al. | Estimating near-surface shear-wave-velocity structures via multichannel analysis of Rayleigh and love waves: An experiment at the Boise hydrogeophysical research site | |
CN100552471C (zh) | 利用垂直地震剖面和双井微测井提高地震分辨率方法 | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
CN112746837A (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集***及方法 | |
CN111239798A (zh) | 光纤声波传感井地地震数据联采***和井驱数据处理方法 | |
CN112130195A (zh) | 基于分布式光纤声波传感的时移vsp数据采集***及方法 | |
Ghose et al. | Integrating S-wave seismic-reflection data and cone penetration test data using a multiangle multiscale approach | |
Zhao et al. | Key seismic survey technologies for deep complex geological structures: A case study of the northern section of the Longmenshan Fault Fold Belt in the Sichuan Basin | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
Sloan et al. | Void detection using near-surface seismic methods | |
CN214576965U (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集*** | |
CN213813970U (zh) | 基于分布式光纤声波传感的时移vsp数据采集*** | |
Khalil et al. | Engineering applications of seismic refraction method: A field example at Wadi Wardan, Northeast Gulf of Suez, Sinai, Egypt | |
Mari | Well seismic surveying and acoustic logging | |
CN109188528B (zh) | 井间弹性波层析成像***和方法 | |
CN100552472C (zh) | 利用垂直地震剖面和微测井进行地震信号补偿方法 | |
Greenwood et al. | Characterization and imaging of a hydrothermally active near-vertical fault zone in crystalline rocks based on hydrophone VSP data | |
Okpoli et al. | Integration of Well logs and seismic attribute analysis in reservoir identification on PGS field onshore Niger Delta, Nigeria | |
Young et al. | Near-surface, SH-wave surveys in unconsolidated, alluvial sediments | |
Juhlin et al. | Imaging of groundwater resources in glacial deposits using high-resolution reflection seismics, Sweden | |
CN104597497A (zh) | 一种基于叠前瞬时频率属性分析的储层烃类预测方法 | |
Ugbor et al. | Hydrocarbon prospectivity and risk assessment of “Bob” field central swamp depobelt, onshore Niger Delta Basin, Nigeria | |
Schwenk | Constrained parameterization of the multichannel analysis of surface waves approach with application at Yuma Proving Ground, Arizona | |
Oldenborger et al. | Radar determination of the spatial structure of hydraulic conductivity |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091021 Termination date: 20110422 |