CN101419292A - 采用纵波源多分量地震数据生成横波地震剖面的方法 - Google Patents
采用纵波源多分量地震数据生成横波地震剖面的方法 Download PDFInfo
- Publication number
- CN101419292A CN101419292A CNA2007101763255A CN200710176325A CN101419292A CN 101419292 A CN101419292 A CN 101419292A CN A2007101763255 A CNA2007101763255 A CN A2007101763255A CN 200710176325 A CN200710176325 A CN 200710176325A CN 101419292 A CN101419292 A CN 101419292A
- Authority
- CN
- China
- Prior art keywords
- wave
- shear wave
- component
- earthquake
- stack
- 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
Abstract
本发明是石油地震勘探纵波源多分量地震数据生成横波地震剖面的方法。步骤是:采用多次覆盖采集三分量地震数据;激发横波井中接收;得到纵波、横波野外静校正量和剩余静校正量进行横波速度分析;根据实际采集记录的波场特征确定切除参数、反褶积参数、滤波参数,完成常规的叠前处理;根据构造及速度场复杂度选择叠加+叠后偏移或叠前偏移,生成剖面。本发明直接利用纵波源激发三分量检波器接收得到的数据提取全程反射横波成像,不需要进行专门的横波采集,利用横波速度动校正后求取炮点横波剩余静校正量,并与速度分析过程一起迭代,得到高静校正精度的横波成像剖面及与之对应的速度场等参数。
Description
技术领域
本发明涉及石油地震勘探数据处理,具体是采用纵波源多分量地震数据生成横波地震剖面的方法。
背景技术
油田开发在很大程度上依赖于钻探的成功率,而钻探成功率的提高与地震勘探成果的精度和可靠性密切相关。目前的勘探技术主要利用纵波勘探来获得地下构造信息,再综合利用其它已有资料确定钻探位置。对于岩性油气藏或裂缝油气藏,利用横波或转换波资料可以获得更多地下物性和岩性信息,从而减少勘探的多解性和开发风险,提高钻探成功率,更有效地促进油田的开发效率。
目前在进行多波多分量勘探时,主要采用纵波激发(通过爆炸源或地面可控震源激发),三分量检波器接收的方式。这种观测方式接收的有效地震波场主要为PP和PS全程反射波。实际上,多分量单点接收时,记录波场非常复杂,各种类型的反射波都可能被观察到,包括在近地表强横波阻抗差界面产生的转换波和由于激发点非均质性产生的横波及与之对应的SS型反射横波和SP型转换波,这些波场一般能量弱,且不稳定。
用三分量检波器进行接收得到矢量地震波场信息,可以通过转换处理的方法获取包含更多的地质构造、岩性变化和油气分布信息的地震剖面,也是油气田开发的有效地球物理勘探方法。
转换波地震资料处理的主要任务是利用野外三分量检波器获得的三分量资料,进行野外采集的转换波地震资料进行波场分离、静校正、非双曲速度分析、抽CCP道集、叠加和叠后偏移等处理,或利用转换波叠前偏移来获得能反映地下地层结构的多波(横波)地震剖面图,反映地层岩性变化的地震波振幅、频率、传播速度等信息,供寻找有利油气圈闭、确定井位使用。
横波地震资料的获取一般采用横波可控震源激发,也可以采用特殊爆炸方式激发。由于横波不受近地表潜水面起伏的影响,纵波和横波的静校正量之间几乎没有关系,用纵波静校正量乘以某一系数来近似计算横波静校正量通常是不正确的。因此,从纵波资料中反演得到横波剖面的方法不直接,反演会存在多解性和不确定性。
发明内容
本发明目的在于提供一种提高地震横波资料可利用性和可靠性,利用纵波源激发三分量检波器接收地震数据进行处理进而得到横波剖面的采用纵波源多分量地震数据生成横波地震剖面的方法。
本发明具体步骤包括:
1)采用爆炸源激发,三分量检波器接收,多次覆盖采集三分量地震数据;
步骤1所述的采集需要通过已知的地质资料获得井深及激发岩性,使激发的横波能量达到较大,同时兼顾转换波反射能量。
2)采用重锤测向敲击方式激发横波,井中三分量检波器接收,通过拾取初至波得到表层纵波及横波速度,建立工区表层纵波、横波速度模型;
步骤2所述的激发横波采用横波可控震源激发。
步骤2所述的接收为地面三分量检波器接收。
3)根据表层纵波、横波速度模型及已知基准面高程,得到纵波、横波野外静校正量;
4)对三分量地震数据进行前期处理,利用野外纵波、横波静校正后转换波数据,采用共检波点叠加反射波平滑法求取横波检波点剩余静校正量;
步骤4所述的前期处理包括观测系统加载,置道头,水平X分量旋转到炮检方向,Y分量旋转到垂直于炮检方向。
5)将横波野外静校正量和横波检波点剩余静校正量应用于水平分量,用通常的速度分析方法进行横波速度分析;
6)利用处理后的横波速度得到的叠加速度对静校正后的水平分量进行动校正,然后利用共炮点叠加反射波平滑静校正方法求取炮点横波剩余静校正;
7)应用炮点横波剩余静校正后重新进行横波速度分析,再重新求取炮点、检波点横波剩余静校正;
步骤7)所述的横波速度分析和剩余静校正要进行迭代处理,直到剩余量小于目标值或叠加效果不再有明显改善。
步骤7)所述的剩余量小于4毫秒。
8)根据实际采集记录的波场特征确定切除参数、反褶积参数、滤波参数,完成常规的叠前处理;
步骤8所述的确定切除参数时将面波区之外的能量切除掉。
步骤8所述的确定滤波时要注意保持低频,滤掉信噪比过低的低频成分。
9)根据构造及速度场复杂度选择叠加+叠后偏移或叠前偏移,生成剖面。
步骤9所述的偏移是对构造平缓的地区,先叠加后再进行叠后偏移。
步骤9所述的偏移是对构造复杂区则要进行叠前时间偏移。
步骤9所述的偏移是对构造和速度场均复杂的区域则要进行叠前深度偏移。
本发明通过对水平分量的处理,得到可用于地质研究的横波成像剖面及与之对应的速度场等参数。具有如下效果:
直接利用纵波源激发三分量检波器接收得到的数据提取全程反射横波,并通过适当的处理流程使之成像,不需要进行专门的横波采集。以往采集的转换波资料也可以用该方法提取横波资料。
以横波速度表层调查为基础,水平分量处理中先借助转换波求取检波点横波剩余静校正量,再利用横波速度动校正后求取炮点横波剩余静校正量,并与速度分析过程一起迭代,提高静校正精度。
附图说明
图1为本发明纵波源激发X分量单炮记录;
图2A为横波三分量微测井的垂直分量,2B为水平径向分量记录图;
图3为纵波表层平均速度模型图;
图4为横波表层平均速度模型图;
图5为横波叠加剖面图;
图6A为本发明纵波、B转换波与C横波叠加剖面对比图。
具体实施方式
本发明具体实施步骤为:
1、野外采用爆炸源激发,三分量检波器接收,通过试验选择合适的井深及激发岩性,使激发的横波能量达到较大,同时兼顾转换波反射能量,完成三分量地震数据采集(图1);
2、采用重锤对受压枕木测向敲击方式或横波可控震源激发横波,井中或地面三分量检波器按一定间隔接收(图2),通过拾取初至波完成表层纵波及横波速度解释,在施工区域通过多点调查建立工区表层纵波、横波速度模型(图3、图4);
3、根据表层纵波、横波速度模型及给定的基准面高程,计算得到纵波、横波野外静校正量;
4、对三分量地震数据进行预处理(包括观测系统加载,置道头,水平X分量旋转到炮检方向,Y分量旋转到垂直于炮检方向),应用野外静校正后在转换波数据上利用共检波点叠加反射波平滑法求取横波检波点剩余静校正量;
5、将横波野外静校正量和横波检波点剩余静校正量应用于水平分量,进行横波速度分析;
6、利用横波速度分析得到的叠加速度对静校正后的水平分量进行动校正,然后利用共炮点叠加反射波平滑静校正方法求取炮点横波剩余静校正;
7、应用炮点横波剩余静校正后重新进行速度分析,再重新求取炮点、检波点横波剩余静校正。速度分析和剩余静校正要进行迭代处理,直到剩余量小于某一给定值或叠加效果不再有明显改善为止;
8、选择合适的切除参数、反褶积参数、滤波参数完成叠前处理。由于横波速度低,反射能量主要集中在面波区之内,因此应将面波区之外的能量切除掉,另外,横波衰减较快,高频能量弱,滤波时要注意保持低频,滤掉信噪比过低的低频成分;
9、根据构造及速度场复杂度选择叠加+叠后偏移或叠前偏移,输出成果剖面(图5)。对于构造平缓的地区,先叠加后再进行叠后偏移就可满足勘探要求,对于构造复杂区则要进行叠前时间偏移,对构造和速度场均复杂的区域则要进行叠前深度偏移。果剖面。经过本发明处理得到信噪比和分辨率均比较理想的横波剖面,与纵波和转换波剖面相比(图6),横波剖面没有受到该区含气地层造成的衰减吸收等因素的影响,准确揭示了精细的构造特征。
Claims (12)
1、一种采用纵波源多分量地震数据生成横波地震剖面的方法,具体步骤包括:
1)采用爆炸源激发,三分量检波器接收,多次覆盖采集三分量地震数据;
2)采用重锤测向敲击方式激发横波,井中三分量检波器接收,通过拾取初至波得到表层纵波及横波速度,建立工区表层纵波、横波速度模型;
3)根据表层纵波、横波速度模型及已知基准面高程,得到纵波、横波野外静校正量;
4)对三分量地震数据进行前期处理,利用野外纵波、横波静校正后转换波数据,采用共检波点叠加反射波平滑法求取横波检波点剩余静校正量;
5)将横波野外静校正量和横波检波点剩余静校正量应用于水平分量,用通常的速度分析方法进行横波速度分析;
6)利用处理后的横波速度得到的叠加速度对静校正后的水平分量进行动校正,然后利用共炮点叠加反射波平滑静校正方法求取炮点横波剩余静校正;
7)应用炮点横波剩余静校正后重新进行横波速度分析,再重新求取炮点、检波点横波剩余静校正;
8)根据实际采集记录的波场特征确定切除参数、反褶积参数、滤波参数,完成常规的叠前处理;
9)根据构造及速度场复杂度选择叠加+叠后偏移或叠前偏移,生成剖面。
2、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤1所述的采集需要通过已知的地质资料获得井深及激发岩性,使激发的横波能量达到较大,同时兼顾转换波反射能量。
3、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤2所述的激发横波采用横波可控震源激发。
4、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤2所述的接收为地面三分量检波器接收。
5、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤4所述的前期处理包括观测系统加载,置道头,水平X分量旋转到炮检方向,Y分量旋转到垂直于炮检方向。
6、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤7)所述的横波速度分析和剩余静校正要进行迭代处理,直到剩余量小于目标值或叠加效果不再有明显改善。
7、根据权利要求1或6所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤7)所述的剩余量小于4毫秒。
8、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤8所述的确定切除参数时将面波区之外的能量切除掉。
9、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤8所述的确定滤波时要注意保持低频,滤掉信噪比过低的低频成分。
10、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤9所述的偏移是对构造平缓的地区,先叠加后再进行叠后偏移。
11、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤9所述的偏移是对构造复杂区则要进行叠前时间偏移。
12、根据权利要求1所述的采用纵波源多分量地震数据生成横波地震剖面的方法,其特征在于:步骤9所述的偏移是对构造和速度场均复杂的区域则要进行叠前深度偏移。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101763255A CN101419292B (zh) | 2007-10-25 | 2007-10-25 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101763255A CN101419292B (zh) | 2007-10-25 | 2007-10-25 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101419292A true CN101419292A (zh) | 2009-04-29 |
CN101419292B CN101419292B (zh) | 2011-05-25 |
Family
ID=40630176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101763255A Active CN101419292B (zh) | 2007-10-25 | 2007-10-25 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101419292B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102221709A (zh) * | 2011-06-01 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于地层参数信息的速度分析与动校正方法 |
CN101852866B (zh) * | 2010-05-13 | 2012-01-04 | 中国石油天然气股份有限公司 | 一种叠后地震数据滤波方法 |
CN103439744A (zh) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 三分量表层数据解释方法 |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN103675916A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种三分量检波器埋置方向高精度校正的方法 |
CN104166157A (zh) * | 2013-05-16 | 2014-11-26 | 中国石油化工股份有限公司 | 基于特殊条件下炸药震源激发三分量接收的横波处理方法 |
CN105093324A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种剩余静校正方法 |
WO2016041186A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种确定均匀介质中最佳激发井深的方法 |
CN110295892A (zh) * | 2019-06-21 | 2019-10-01 | 中国石油大学(北京) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 |
CN110612462A (zh) * | 2017-05-24 | 2019-12-24 | 通用电气公司 | 用于从井孔进行地层评估的系统和方法 |
CN111751876A (zh) * | 2020-07-07 | 2020-10-09 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置 |
CN112363223A (zh) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | 横波速度的获取方法及装置 |
CN112379443A (zh) * | 2020-11-03 | 2021-02-19 | 中国石油天然气集团有限公司 | 纵横波可控震源微测井表层分析系统及方法 |
CN112462426A (zh) * | 2020-11-02 | 2021-03-09 | 中国石油天然气集团有限公司 | 横波矢量静校正方法及装置 |
CN113031073A (zh) * | 2019-12-24 | 2021-06-25 | 中国石油天然气股份有限公司 | 叠前深度偏移速度确定方法、装置及设备 |
CN114415234A (zh) * | 2022-01-24 | 2022-04-29 | 西南交通大学 | 基于主动源面波频散和h/v确定浅地表横波速度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4736349A (en) * | 1987-04-24 | 1988-04-05 | Mobil Oil Corporation | Method for estimating shear wave reflection data from acquired compressional wave reflection data |
-
2007
- 2007-10-25 CN CN2007101763255A patent/CN101419292B/zh active Active
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101852866B (zh) * | 2010-05-13 | 2012-01-04 | 中国石油天然气股份有限公司 | 一种叠后地震数据滤波方法 |
CN102221709B (zh) * | 2011-06-01 | 2013-05-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于地层参数信息的速度分析与动校正方法 |
CN102221709A (zh) * | 2011-06-01 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于地层参数信息的速度分析与动校正方法 |
CN103675916B (zh) * | 2012-09-04 | 2016-12-21 | 中国石油天然气集团公司 | 一种三分量检波器埋置方向高精度校正的方法 |
CN103675916A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种三分量检波器埋置方向高精度校正的方法 |
CN104166157A (zh) * | 2013-05-16 | 2014-11-26 | 中国石油化工股份有限公司 | 基于特殊条件下炸药震源激发三分量接收的横波处理方法 |
CN103439744A (zh) * | 2013-08-27 | 2013-12-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 三分量表层数据解释方法 |
CN103487834B (zh) * | 2013-09-09 | 2016-07-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN105093324A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种剩余静校正方法 |
CN105093324B (zh) * | 2014-05-21 | 2017-06-16 | 中国石油化工股份有限公司 | 一种剩余静校正方法 |
WO2016041186A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种确定均匀介质中最佳激发井深的方法 |
CN110612462A (zh) * | 2017-05-24 | 2019-12-24 | 通用电气公司 | 用于从井孔进行地层评估的系统和方法 |
CN110295892A (zh) * | 2019-06-21 | 2019-10-01 | 中国石油大学(北京) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 |
CN110295892B (zh) * | 2019-06-21 | 2020-07-10 | 中国石油大学(北京) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 |
CN113031073A (zh) * | 2019-12-24 | 2021-06-25 | 中国石油天然气股份有限公司 | 叠前深度偏移速度确定方法、装置及设备 |
CN113031073B (zh) * | 2019-12-24 | 2024-03-01 | 中国石油天然气股份有限公司 | 叠前深度偏移速度确定方法、装置及设备 |
CN111751876A (zh) * | 2020-07-07 | 2020-10-09 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置 |
CN111751876B (zh) * | 2020-07-07 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp转换横波单程波叠前深度偏移方法和装置 |
CN112363223A (zh) * | 2020-10-29 | 2021-02-12 | 中国石油天然气集团有限公司 | 横波速度的获取方法及装置 |
CN112462426A (zh) * | 2020-11-02 | 2021-03-09 | 中国石油天然气集团有限公司 | 横波矢量静校正方法及装置 |
CN112379443A (zh) * | 2020-11-03 | 2021-02-19 | 中国石油天然气集团有限公司 | 纵横波可控震源微测井表层分析系统及方法 |
CN114415234A (zh) * | 2022-01-24 | 2022-04-29 | 西南交通大学 | 基于主动源面波频散和h/v确定浅地表横波速度的方法 |
CN114415234B (zh) * | 2022-01-24 | 2023-05-12 | 西南交通大学 | 基于主动源面波频散和h/v确定浅地表横波速度的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101419292B (zh) | 2011-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101419292B (zh) | 采用纵波源多分量地震数据生成横波地震剖面的方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN101598805B (zh) | 一种多分量地震资料层位对比和标定方法 | |
CN101281253B (zh) | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 | |
O'Neill et al. | Full-waveform P-SV reflectivity inversion of surface waves for shallow engineering applications | |
CN102176054B (zh) | 近地表综合信息处理解释方法 | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
US6055482A (en) | Method of seismic signal processing | |
CN103116184B (zh) | 一种表层调查中计算横波静校正量的方法 | |
CN101329405A (zh) | 一种简单的多参数地震反演方法 | |
CN101551466A (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
Witten et al. | Microseismic image-domain velocity inversion: Marcellus Shale case study | |
Onyebueke et al. | High-resolution shallow reflection seismic integrated with other geophysical methods for hydrogeological prospecting in the Nylsvley Nature Reserve, South Africa | |
Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
Kamei et al. | Application of waveform tomography to a crooked-line 2D land seismic data set | |
US20140257707A1 (en) | Progressive 3D Vertical Seismic Profiling Method | |
Ding et al. | Reverse time migration (RTM) imaging of iron-oxide deposits in the Ludvika mining area, Sweden | |
US20120269035A1 (en) | Evaluating Prospects from P-Wave Seismic Data Using S-Wave Vertical Shear Profile Data | |
Tselentis et al. | Application of a high-resolution seismic investigation in a Greek coal mine | |
Francese et al. | High-resolution images of shallow aquifers—A challenge in near-surface seismology | |
CN109839659B (zh) | 一种对叠前深度偏移剖面进行迭代优化的方法 | |
CN101290356A (zh) | 一种炸药震源勘探的辅助方法 | |
Zong et al. | Fractured basement imaging using random-space-shift reverse time migration: A vertical seismic profile survey in the Bohai Bay Basin, China | |
Mari et al. | Well seismic surveying |
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 |