CN106094028A - 基于位移场的深度域地震属性提取方法 - Google Patents
基于位移场的深度域地震属性提取方法 Download PDFInfo
- Publication number
- CN106094028A CN106094028A CN201610458965.4A CN201610458965A CN106094028A CN 106094028 A CN106094028 A CN 106094028A CN 201610458965 A CN201610458965 A CN 201610458965A CN 106094028 A CN106094028 A CN 106094028A
- Authority
- CN
- China
- Prior art keywords
- displacement field
- depth domain
- attribute
- seismic
- method based
- 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
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 31
- 210000002706 plastid Anatomy 0.000 claims abstract description 24
- 239000000284 extract Substances 0.000 claims abstract description 15
- 239000010410 layer Substances 0.000 claims description 21
- 230000008859 change Effects 0.000 claims description 15
- 239000011229 interlayer Substances 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000013075 data extraction Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 10
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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
本发明提供一种基于位移场的深度域地震属性提取方法,该基于位移场的深度域地震属性提取方法包括:步骤1,在深度域地震数据上进行地质层位解释和断层解释;步骤2,在深度域中,分析地质体的位移场信息;步骤3,进行基于位移场的深度域地震属性提取;步骤4,通过深度域地震属性结果分析地质体的分布规律。该基于位移场的深度域地震属性提取方法提供一个对应解释层位的包含地质体信息的深度位移场,以此为提取属性时时窗定义的基础。通过位移场提取属性,从而得到地质体的分布特征。结果更加客观、准确。
Description
技术领域
本发明涉及油田勘探开发技术领域,特别是涉及到一种基于位移场的深度域地震属性提取方法。
背景技术
叠前深度偏移技术的研究是近10多年来全球油气地球物理勘探领域的热点,是目前地震资料处理技术中归位相对准确的一种偏移方法,偏移的作用是使地下界面的地震反射波归位到正确的空间位置,最终得到真实反映地下界面形态的地震图像。特别是在速度复杂或者构造复杂地区具有较明显的优势,主要用来解决盐丘、复杂基底、复杂断块、生物礁、泥底辟、流体底辟等复杂地质体的波场成像问题。随着偏移算法的成熟、计算机能力的进步,规模化的叠前深度偏移成为可能。深度域地震剖面与时间域地震剖面相比,具有地下构造真实、直观、便于解释等特点,并且无需进行时深转换就能得到初始深度构造。目前,叠前深度偏移已经是工业界的标准处理流程。
深度域叠前地震属性技术近几年较引人关注,其原因从主观来讲,是人们对储层描述精度要求的提高;从客观来看,则主要是由于硬件条件的改善,如计算机技术的发展,并由此带动地震成像技术,叠前地震属性提取和分析等相关技术的发展。在这样一个大的环境下,开展深度域叠前地震属性技术的探索与应用十分必要。基于波场外推理论的波动方程偏移可以使复杂构造精确成像,通常所用的成像条件是在每一深度的零偏移距抽取零时波场。Mosher和Foster(2000)提出了叠前深度偏移共角度成像条件,在每一深度对于一常射线参数抽取零时波场,一系列的零时偏移距射线参数地震道便生成共角度(偏移距射线参数)成像道集,该成像道集与偏移后的道集是等价的,可以用来做深度域属性分析。Biondi(1996)提出了共方位地震数据在叠前偏移中的作用,Prucha(1999)等提出了基于波动方程偏移形成共角度成像道集,Mosher(1996)等分析了偏移成像对AVO的影响。
胡中平等(2004)结合国内外研究动态,分析了叠前深度域地震属性提取的可行性、必要性以及相对于CMP道集技术的优越性。在复杂地质条件下用常规CMP道集进行地震属性分析会存在精度等方面的问题。用叠前深度域地震属性技术可以完成复杂构造条件下高精度储层预测。关于叠前深度域地震属性研究的实例不多,但此方法已经引起国内外地球物理学家的关注。初步分析研究表明,叠前深度域地震属性技术可行,具有广阔的应用前景。在深度域频谱分析研究方面,陈茂山等(1999)在傅里叶变换基础上研究测井资料的频谱信息特征及其与高频层序的关系,以用于层序地层学研究中的高频层序划分。提出了两种基于测井资料的深度域频谱分析方法—深度域频率扫描和滑动窗频谱分析,在阿南凹陷层序地层学研究中取得了良好的效果。蔡瑞等(2008)通过深度域谱分解技术的应用研究,分析孔隙变化在敏感波数振幅上的响应,揭示各目的层内部岩性横向变化特征。在深度域,同一套地层速度相同,反射波组有相似性,波数一振幅特征也应相似,当波数振幅特征发生改变时,说明地层岩性横向上已经发生变化。由于地震剖面从浅至深平均速度由小变大,大地吸收衰减由弱变强,频率特征由强变弱,因此,纵向上不同时代地层振幅有波数变化,不同波数剖面可反映各目的层岩性横向变化特征。
为此我们发明了一种新的基于位移场的深度域地震属性提取方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种通过包含地质体的位移场的地震属性提取,从而得到地质体的分布特征的基于位移场的深度域地震属性提取方法。
本发明的目的可通过如下技术措施来实现:基于位移场的深度域地震属性提取方法,该基于位移场的深度域地震属性提取方法包括:步骤1,在深度域地震数据上进行地质层位解释和断层解释;步骤2,在深度域中,分析地质体的位移场信息;步骤3,进行基于位移场的深度域地震属性提取;步骤4,通过深度域地震属性结果分析地质体的分布规律。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,在深度域地震数据上,通过精细标定把砂层组的地震反射同相轴从地震剖面上识别出来,以区别其它反射轴,通过精细标定建立不同岩性、电性的地层与地震反射特征的对应关系,最终结合地质认识进行地质层位解释和断层解释。
在步骤2中,深度域中同一套地层速度相同,反射波组有相似性,波数一振幅特征也应相似,当波数振幅特征发生改变时,说明地层岩性横向上已经发生变化,由于地震剖面从浅至深平均速度由小变大,大地吸收衰减由弱变强,频率特征由强变弱,因此,纵向上不同时代地层振幅有波数变化,不同波数剖面可反映各目的层岩性横向变化特征;分析地质体结构特征,确定地质体的位移场特性。
在步骤3中,当通过步骤2确定地质体的位移场后,进行深度域地震属性提取,其中需要计算的地震属性包括体属性、沿层属性、沿层体属性、层间属性。
在步骤3中,在计算沿层属性时,根据地质解释层位或移动解释层位进行地震数据提取;在计算层间体属性时,根据地质体位移场的信息,包含地质体位移场的提取层间地震属性,通过地质体的结构特征,确定地质体的位移场特性,通过位移场特性来反应波数和幅度确定沿层时窗或层间深度大小,提取深度域地震的幅度,波数和相位及众多的派生属性。
在步骤4中,通过在深度域中,提取层间属性或是提取基于位移场的地震属性的结果,分析地质体的平面展布规律。
本发明中的基于位移场的深度域地震属性提取方法,该方法利用地质体信息的位移场提取地震属性,偏移的作用是使地下界面的地震反射波归位到正确的空间位置,最终得到真实反映地下界面形态的地震图像。时间域地震属性主要包含了振幅、频率、相位、吸收衰减及波形等。而在深度域子波发生变化,系统本身会改变信号形状,许多时间域适用的属性分析方法并不能直接应用于深度域,特别是频谱类属性。深度域地震反射的畸变使得属性提取时相应的深度窗口的定义也发生变化,本发明提供一个对应解释层位的包含地质体信息的深度位移场,以此为提取属性时时窗定义的基础,综合地质认识,通过位移场进行属性提取,从而得到地质体的分布特征。结果更加客观、准确。
附图说明
图1为本发明的基于位移场的深度域地震属性提取方法的一具体实施例的流程图;
图2为本发明的一具体实施例中深度域地震解释层位剖面示意图;
图3为本发明的一具体实施例中基于位移场的深度域地震瞬变波数剖面示意图;
图4为本发明的一具体实施例中基于位移场的深度域地震属性平面示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的基于位移场的深度域地震属性提取方法的流程图。
在步骤101,在深度域地震数据上进行地质层位解释和断层解释。
在深度域地震数据上,通过精细标定把砂层组的地震反射同相轴从地震剖面上识别出来,以区别其它反射轴,通过精细标定建立不同岩性、电性的地层与地震反射特征的对应关系,最终结合地质认识进行地质层位解释和断层解释(如图2)。流程进入到步骤102。
在步骤102,深度域中,分析地质体的位移场信息。
在深度域,同一套地层速度相同,反射波组有相似性,波数一振幅特征也应相似,当波数振幅特征发生改变时,说明地层岩性横向上已经发生变化。由于地震剖面从浅至深平均速度由小变大,大地吸收衰减由弱变强,频率特征由强变弱,因此,纵向上不同时代地层振幅有波数变化,不同波数剖面(如图3)可反映各目的层岩性横向变化特征。分析地质体结构特征,确定地质体的位移场特性。流程进入到步骤103。
在步骤103,通过基于位移场的深度域地震属性提取。
通过步骤102中,确定地质体的位移场后,进行深度域地震属性提取,其中需要计算的地震属性包括体属性、沿层属性、沿层体属性、在上下两层之间,并且有明显的特征,在深度域中提取层间的地震体属性。在计算沿层属性时,根据地质解释层位或移动解释层位进行地震数据提取;在计算层间体属性时,根据地质体位移场的信息,包含地质体位移场的提取层间地震属性。通过地质体的结构特征,确定地质体的位移场特性,通过位移场特性来反应波数和幅度确定沿层时窗或层间深度大小,提取深度域地震的幅度,波数和相位及众多的派生属性。流程进入到步骤104。
在步骤104,通过深度域地震属性结果分析地质体的分布规律。
通过在深度域中,提取层间属性或是提取基于位移场的地震属性(如图4)的结果。通过平面的属性,分析目的层的沉积体系,从而分析地质体的平面展布规律,为研究人员提供直观、客观、准确的结果。
本发明中的基于位移场的深度域地震属性提取方法,通过在深度域中,解释地质层位和断层的基础上,对研究区中地质体的位移场信息分析,确定地质体是在层间还是在上下两层之间的,如果是地质体位于层间,就提取层间地震属性即可,如果地质体位于解释层位中上下位置,就通过地质体的结构特征,确定地质体的位移场特性,通过位移场特性来反应波数和幅度确定沿层深度时窗大小,提取深度域地震的幅度,波数和相位及众多的派生属性。最终对深度域的地震属性结果进行分析,得出地质体的平面展布规律,为研究人员提供直观、客观、准确的研究结果。
Claims (6)
1.基于位移场的深度域地震属性提取方法,其特征在于,该基于位移场的深度域地震属性提取方法包括:
步骤1,在深度域地震数据上进行地质层位解释和断层解释;
步骤2,在深度域中,分析地质体的位移场信息;
步骤3,进行基于位移场的深度域地震属性提取;
步骤4,通过深度域地震属性结果分析地质体的分布规律。
2.根据权利要求1所述的基于位移场的深度域地震属性提取方法,其特征在于,在步骤1中,在深度域地震数据上,通过精细标定把砂层组的地震反射同相轴从地震剖面上识别出来,以区别其它反射轴,通过精细标定建立不同岩性、电性的地层与地震反射特征的对应关系,最终结合地质认识进行地质层位解释和断层解释。
3.根据权利要求1所述的基于位移场的深度域地震属性提取方法,在步骤2中,深度域中同一套地层速度相同,反射波组有相似性,波数一振幅特征也应相似,当波数振幅特征发生改变时,说明地层岩性横向上已经发生变化,由于地震剖面从浅至深平均速度由小变大,大地吸收衰减由弱变强,频率特征由强变弱,因此,纵向上不同时代地层振幅有波数变化,不同波数剖面可反映各目的层岩性横向变化特征;分析地质体结构特征,确定地质体的位移场特性。
4.根据权利要求1所述的基于位移场的深度域地震属性提取方法,其特征在于,在步骤3中,当通过步骤2确定地质体的位移场后,进行深度域地震属性提取,其中需要计算的地震属性包括体属性、沿层属性、沿层体属性、层间属性。
5.根据权利要求4所述的基于位移场的深度域地震属性提取方法,其特征在于,在步骤3中,在计算沿层属性时,根据地质解释层位或移动解释层位进行地震数据提取;在计算层间体属性时,根据地质体位移场的信息,包含地质体位移场的提取层间地震属性,通过地质体的结构特征,确定地质体的位移场特性,通过位移场特性来反应波数和幅度确定沿层时窗或层间深度大小,提取深度域地震的幅度,波数和相位及众多的派生属性。
6.根据权利要求1所述的基于位移场的深度域地震属性提取方法,其特征在于,在步骤4中,通过在深度域中,提取层间属性或是提取基于位移场的地震属性的结果,分析地质体的平面展布规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610458965.4A CN106094028A (zh) | 2016-06-22 | 2016-06-22 | 基于位移场的深度域地震属性提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610458965.4A CN106094028A (zh) | 2016-06-22 | 2016-06-22 | 基于位移场的深度域地震属性提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106094028A true CN106094028A (zh) | 2016-11-09 |
Family
ID=57239016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610458965.4A Pending CN106094028A (zh) | 2016-06-22 | 2016-06-22 | 基于位移场的深度域地震属性提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106094028A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107831539A (zh) * | 2017-08-16 | 2018-03-23 | 中国石油化工股份有限公司 | 层间积分超剥线识别及提取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070247973A1 (en) * | 2006-04-21 | 2007-10-25 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
CN104884974A (zh) * | 2012-12-05 | 2015-09-02 | 兰德马克绘图国际公司 | 利用人工神经网络进行3d地震数据深度转换的系统和方法 |
CN105353412A (zh) * | 2015-12-14 | 2016-02-24 | 中国石油大学(华东) | 一种井震联合平均速度场的计算方法及系统 |
-
2016
- 2016-06-22 CN CN201610458965.4A patent/CN106094028A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070247973A1 (en) * | 2006-04-21 | 2007-10-25 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
CN104884974A (zh) * | 2012-12-05 | 2015-09-02 | 兰德马克绘图国际公司 | 利用人工神经网络进行3d地震数据深度转换的系统和方法 |
MX2015005627A (es) * | 2012-12-05 | 2016-02-03 | Landmark Graphics Corp | Sistemas y metodos para conversion de profundidad de datos sismicos 3d utlizando redes neuronales artificiales. |
CN105353412A (zh) * | 2015-12-14 | 2016-02-24 | 中国石油大学(华东) | 一种井震联合平均速度场的计算方法及系统 |
Non-Patent Citations (3)
Title |
---|
胡中平等: "叠前深度域地震属性技术研究及展望", 《石油物探》 * |
蔡瑞 等: "深度域波数谱剖面特征分析", 《中国地球物理学会第二十四届年会论文集》 * |
隋淑玲: "井间地震资料的属性分析及应用", 《石油勘探与开发》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107831539A (zh) * | 2017-08-16 | 2018-03-23 | 中国石油化工股份有限公司 | 层间积分超剥线识别及提取方法 |
CN107831539B (zh) * | 2017-08-16 | 2019-04-09 | 中国石油化工股份有限公司 | 层间积分超剥线识别及提取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3400463B1 (en) | Amplitude inversion on partitioned depth image gathers using point spread function | |
Chopra et al. | Emerging and future trends in seismic attributes | |
West et al. | Interactive seismic facies classification using textural attributes and neural networks | |
Al-Dossary et al. | 3D volumetric multispectral estimates of reflector curvature and rotation | |
EP3488267B1 (en) | Seismic spectral balancing | |
Zhang et al. | Passive seismic tomography using induced seismicity at a petroleum field in Oman | |
Burnett et al. | Seismic diffraction interpretation at Piceance Creek | |
CN109541685B (zh) | 一种河道砂体识别方法 | |
CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN105425299A (zh) | 确定地层裂缝分布的方法和装置 | |
US11561312B2 (en) | Mapping near-surface heterogeneities in a subterranean formation | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
Barthwal et al. | Passive seismic tomography using recorded microseismicity: Application to mining-induced seismicity | |
CN115877449A (zh) | 用于在勘测区域内获得地下堆叠图像的计算机实现方法 | |
Martorana et al. | Assessment of a geological model by surface wave analyses | |
CN103076630B (zh) | 一种基于弹性阻抗梯度的油气检测方法 | |
CN105093266A (zh) | 地震地质综合描述地层超覆线方法 | |
Elebiju et al. | Investigation of links between Precambrian basement structure and Paleozoic strata in the Fort Worth Basin, Texas, USA, using high-resolution aeromagnetic (HRAM) data and seismic attributes | |
CN105954800B (zh) | 利用地震纹理特性识别小断层的方法 | |
Liu et al. | Seismic characterization of fault and fractures in deep buried carbonate reservoirs using CNN-LSTM based deep neural networks | |
Zhou et al. | Seismic diffraction imaging for improved coal structure detection by a structure-oriented moving average error filter | |
CN106094028A (zh) | 基于位移场的深度域地震属性提取方法 | |
CN104820237A (zh) | 计算地层的方差体的方法 | |
US11531129B2 (en) | Picking seismic stacking velocity based on structures in a subterranean formation | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161109 |
|
RJ01 | Rejection of invention patent application after publication |