CN102565852B - 针对储层含油气性检测的角度域叠前偏移数据处理方法 - Google Patents
针对储层含油气性检测的角度域叠前偏移数据处理方法 Download PDFInfo
- Publication number
- CN102565852B CN102565852B CN 201010579495 CN201010579495A CN102565852B CN 102565852 B CN102565852 B CN 102565852B CN 201010579495 CN201010579495 CN 201010579495 CN 201010579495 A CN201010579495 A CN 201010579495A CN 102565852 B CN102565852 B CN 102565852B
- Authority
- CN
- China
- Prior art keywords
- imaging
- angle
- angle domain
- cos
- alpha
- 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.)
- Expired - Fee Related
Links
Images
Abstract
本发明涉及一种针对储层含油气性检测的角度域叠前偏移数据处理方法,包括以下步骤:1)模拟地面地震,数据采集器采集地震波信号,按振幅相对保真处理流程进行预处理;2)输入设备根据储层预测需要选定局部目标成像范围,并对角度域成像数组进行初始化;3)数据采集器将步骤1)预处理后的地震波数据逐道输入到处理器,处理器根据中心点位置与偏移孔径限制,确定该道地震波对应的最大成像范围,并判断该道最大成像范围是否在选定的局部目标成像范围内,若为是,执行步骤4),若为否,输入下一道等。与现有技术相比,本发明具有根据地震波局部方向特征实现局部角度域叠前时间偏移成像,直接在偏移过程中一步获得高保真的入射角域共成像点道集等优点。
Description
技术领域
本发明涉及一种角度域叠前偏移数据处理方法,尤其是涉及一种针对储层含油气性检测的角度域叠前偏移数据处理方法。
背景技术
地震勘探是通过人工激发地震波,在地下传播的波场信号到达地表检波器或传感器被记录下来。这些波场记录中携带了地下地质结构、岩石弹性性质与孔隙流体性质有关的响应信息。根据地震波在地层岩石中的传播规律,利用合理的波场传播模型和地震数据处理、反演与解释方法,结合对研究区域的地质认识寻找石油、天然气等化石能源的地球物理探测方法。
地震勘探一方面利用观测波场的运动学信息获得地下宏观速度结构和构造图像,在此基础上利用波场的动力学信息估算地层弹性参数,预测储层及其含油气性质。因此,从地震资料中估计地层介质的弹性参数,并结合测井与岩石物理资料将弹性参数与岩性与流体性质联系起来,对油气储层预测与描述至关重要。
目前,油气预测技术仍然比较匮乏。从地震数据中恢复地层波阻抗界面对应反射振幅随偏移距变化(Amplitude versus Offset,简称AVO)或者反射振幅随入射角变化(Amplitude versus Incident-angle,简称AVA)分析或反演,以及物理基础相同的弹性阻抗(Elastic Impedance,简称EI)反演等叠前反演技术在储层与油气分布预测方面扮演重要角色。其基本原理是从地面观测地震记录通过恰当的处理恢复目标界面反射振幅随偏移距或入射角的变化特征,通过Zeoppritz方程或其近似形式将AVO/AVA响应与地层弹性参数与物性参数联系起来,通过分析或反演获得能够指示储层、区分流体类型与饱和情况的AVO属性或弹性参数等。该反演过程的数据基础和关键是获得相对保真的“角道集”或“共角度剖面”。常规做法是采用“两步法”,首先利用正常时差(NMO)校正或叠前时间偏移获得偏移距道集,再由偏移距-入射角映射获得“角道集”(Castagna and Backus,1993)。近来研究表明,第一步中使用的偏移距域成像会引起振幅在角度域的混叠(Zheng,2006),第二步中入射角的计算没有考虑下行与上行射线非对称情况,仅适用于水平层状介质(程玖兵等,2009)。这些限制降低了所获得“角道集”的AVA响应的保真度,会降低在此基础上的储层参数与流体参数反演精度。目前一些地质情况较复杂地区油气预测成功率不高与传统方法这些缺陷不无关系。AVO/AVA反演前处理常常采用的Kirchhoff叠前时间偏移算法只能得到偏移距域共成像点道集。为此,王楠等(2008)和程玖兵等(2009)提出了基于局部一维介质射线追踪走时表驱动叠前时间偏移角度域成像方法。该方法在横向非均匀性很弱的地区效率与精度都比较高。但随着横向非均匀性增强,需要根据速度模型实现多次局部一维射线追踪,然后在成像时插值获得实际位置对应射线路径的走时、角度与加权系数信息,这样做效率降低很多,缺少灵活性,而且精度与射线追踪的密度与速度变化特点有关。该类方法目前还没有扩展的各向异性介质情况。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种针对储层含油气性检测的角度域叠前偏移数据处理方法。
本发明的目的可以通过以下技术方案来实现:
一种针对储层含油气性检测的角度域叠前偏移数据处理方法,其特征在于,包括以下步骤:
1)模拟地面地震,数据采集器采集地震波信号,按振幅相对保真处理流程进行预处理;
2)输入设备根据储层预测需要选定局部目标成像范围,并对角度域成像数组进行初始化;
3)数据采集器将步骤1)预处理后的地震波数据逐道输入到处理器,处理器根据中心点位置与偏移孔径限制,确定该道地震波对应的最大成像范围,并判断该道最大成像范围是否在选定的局部目标成像范围内,若为是,执行步骤4),若为否,输入下一道;
4)处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理,判断该道地震波所有的成像时间的角度域成像处理是否完成,若为是,执行步骤5),若为否,继续执行步骤4);
5)处理器判断所有的地震波数据是否都理处完成,若为是,执行步骤6),若为否,返回步骤3);
6)得到局部成像空间的角度域共成像点道集和部分角度叠加剖面,基于角度域共成像点道集进行振幅随入射角变化AVA分析或反演,或基于部分角度叠加剖面进行弹性阻抗EI反演,获得反映储层及其含油气性的属性数据与图表,并通过显示器显示。
所述的步骤4)中的处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理过程如下:
1)从角度域成像空间开始时刻进行偏移脉冲响应计算:根据当前成像时间(τ)的偏移孔径范围以及中心点坐标,确定当前成像时间成像点横向分布区间;
2)在成像点横向分布区间内,针对每个成像点(x,y,τ),根据非双曲双平方根DSR时差方程计算炮点-成像点单程射线走时ts以及成像点-接收点单程射线走时tr,再利用公式:计算入射与散射射线参数ps与pr,再根据Snell定律计算射线的起飞角αs与αr,二维各向同性介质情况下按γ=(αs+αr)/2计算得到入射角,在三维各向同性介质情况下,通过
计算得到入射角,其中β是入射射线所在垂直平面与散射射线所在垂直平面的夹角。对于二维情况,通过如下公式计算几何扩散校正权系数:
对于三维情况,通过如下公式计算几何扩散校正权系数:
3)将输入地震道对应走时的振幅乘以几何扩散校正权系数后叠加到角度域成像空间(x,y,γ,τ),其中x,y,τ为时间域成像坐标,γ为入射角;完成当前成像时间孔径内所有点的角度域成像后,进入下一个成像时间,返回步骤1)。
所述的步骤1)中预处理包括叠前去噪、静校正、地表一致性校正。
所述的步骤2)中的局部目标成像范围包括成像空间起始iLine号与xLine号、起始自激自收时间范围。
与现有技术相比,本发明具有以下优点:
1)针对复杂地质情况,根据地震波局部方向特征实现局部角度域叠前时间偏移成像,直接在偏移过程中一步获得高保真的入射角域共成像点道集,可为储层及其含油气预测提供目的层“原位”的振幅随偏移距或入射角的变化资料,克服了传统“角道集”生成方法精度较低的问题。
2)考虑了大偏移距射线弯曲效应与速度各向异性,也符合目前AVO/AVA反演、EI反演主要基于叠前时间偏移成果数据的现状,拓展了它们的应用范围,能提高储层含油气性检测的准确度。技术流程实用性强,国内外尚无同类方法技术公开发表。
附图说明
图1为本发明的流程图;
图2为本发明的硬件结构示意图。
图3为本发明提取的小、中、大角度成像剖面图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
实施角度域叠前时间偏移处理之前,在常规处理流程中对叠前地震资料做到振幅相对保真的噪音压制、静校正以及地表一致性校正等准备工作。建议利用井资料基础上的AVO/AVA正演对角度域成像结果进行质量控制。
先由常规叠前时间偏移处理获得合理的速度模型供偏移使用,并转换生成相对准确的层速度供角度计算使用。然后按照本发明技术方案实现角度域叠前时间偏移成像,提高储层及其含油气性预测的精度,提高油气勘探开发的成功率。
如图1、图2所示,一种针对储层含油气性检测的角度域叠前偏移数据处理方法,包括以下步骤:
1)模拟地面地震,数据采集器1采集地震波信号,按振幅相对保真处理流程进行预处理;
2)输入设备4根据储层预测需要选定局部目标成像范围,并对角度域成像数组进行初始化;
3)数据采集器1将步骤1)预处理后的地震波数据逐道输入到处理器2,处理器2根据中心点位置与偏移孔径限制,确定该道地震波对应的最大成像范围,并判断该道最大成像范围是否在选定的局部目标成像范围内,若为是,执行步骤4),若为否,输入下一道;
4)处理器2对该道地震波进行每个成像时间孔径内所有点的角度域成像处理,判断该道地震波所有的成像时间的角度域成像处理是否完成,若为是,执行步骤5),若为否,继续执行步骤4);
5)处理器2判断所有的地震波数据是否都理处完成,若为是,执行步骤6),若为否,返回步骤3);
6)得到局部成像空间的角度域共成像点道集和部分角度叠加剖面,基于角度域共成像点道集进行振幅随入射角变化AVA分析或反演,或基于部分角度叠加剖面进行弹性阻抗EI反演,获得反映储层及其含油气性的属性数据与图表,并通过显示器3显示(如图3)。
所述的步骤4)中的处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理过程如下:
1)从角度域成像空间开始时刻进行偏移脉冲响应计算:根据当前成像时间(τ)的偏移孔径范围以及中心点坐标,确定当前成像时间成像点横向分布区间;
2)在成像点横向分布区间内,针对每个成像点(x,y,τ),根据非双曲双平方根DSR时差方程计算炮点-成像点单程射线走时ts以及成像点-接收点单程射线走时tr,再利用公式:计算入射与散射射线参数ps与pr,再根据Snell定律计算射线的起飞角αs与αr,二维各向同性介质情况下按γ=(αs+αr)/2计算得到入射角,在三维各向同性介质情况下,通过
计算得到入射角,其中β是入射射线所在垂直平面与散射射线所在垂直平面的夹角。对于二维情况,通过如下公式计算几何扩散校正权系数:
对于三维情况,通过如下公式计算几何扩散校正权系数:
3)将输入地震道对应走时的振幅乘以几何扩散校正权系数后叠加到角度域成像空间(x,y,γ,τ),其中x,y,τ为时间域成像坐标,γ为入射角;完成当前成像时间孔径内所有点的角度域成像后,进入下一个成像时间,返回步骤1)。
Claims (3)
1.一种针对储层含油气性检测的角度域叠前偏移数据处理方法,其特征在于,包括以下步骤:
1)模拟地面地震,数据采集器采集地震波信号,按振幅相对保真处理流程进行预处理;
2)输入设备根据储层预测需要选定局部目标成像范围,并对角度域成像数组进行初始化;
3)数据采集器将步骤1)预处理后的地震波数据逐道输入到处理器,处理器根据中心点位置与偏移孔径限制,确定该道地震波对应的最大成像范围,并判断该道最大成像范围是否在选定的局部目标成像范围内,若为是,执行步骤4),若为否,输入下一道;
4)处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理,判断该道地震波所有的成像时间的角度域成像处理是否完成,若为是,执行步骤5),若为否,继续执行步骤4);
5)处理器判断所有的地震波数据是否都理处完成,若为是,执行步骤6),若为否,返回步骤3);
6)得到局部成像空间的角度域共成像点道集和部分角度叠加剖面,基于角度域共成像点道集进行振幅随入射角变化AVA分析或反演,或基于部分角度叠加剖面进行弹性阻抗EI反演,获得反映储层及其含油气性的属性数据与图表,并通过显示器显示;
所述的步骤4)中的处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理过程如下:
(1)从角度域成像空间开始时刻进行偏移脉冲响应计算:根据当前成像时间τ的偏移孔径范围以及中心点坐标,确定当前成像时间成像点横向分布区间;
(2)在成像点横向分布区间内,针对每个成像点(x,y,τ),根据非双曲双平方根DSR时差方程计算炮点-成像点单程射线走时ts以及成像点-接收点单程射线走时tr,再利用公式:计算入射与散射射线参数ps与pr,再根据Snell定律计算射线的起飞角αs与αr,二维各向同性介质情况下按γ2D=(αs+αr)/2计算得到入射角;在三维各向同性介质情况下,通过
计算得到入射角,其中β是入射射线所在垂直平面与散射射线所在垂直平面的夹角;对于二维情况,通过如下公式计算几何扩散校正权系数:
对于三维情况,通过如下公式计算几何扩散校正权系数:
3)将输入地震道对应走时的振幅乘以几何扩散校正权系数后叠加到角度域成像空间(x,y,γ,τ),其中x,y,τ为时间域成像坐标,γ为入射角;完成当前成像时间孔径内所有点的角度域成像后,进入下一个成像时间,返回步骤1)。
2.根据权利要求1所述的一种针对储层含油气性检测的角度域叠前偏移数据处理方法,其特征在于,所述的步骤1)中预处理包括叠前去噪、静校正、地表一致性校正。
3.根据权利要求1所述的一种针对储层含油气性检测的角度域叠前偏移数据处理方法,其特征在于,所述的步骤2)中的局部目标成像范围包括成像空间起始iLine号与xLine号、起始自激自收时间范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010579495 CN102565852B (zh) | 2010-12-08 | 2010-12-08 | 针对储层含油气性检测的角度域叠前偏移数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010579495 CN102565852B (zh) | 2010-12-08 | 2010-12-08 | 针对储层含油气性检测的角度域叠前偏移数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102565852A CN102565852A (zh) | 2012-07-11 |
CN102565852B true CN102565852B (zh) | 2013-10-16 |
Family
ID=46411718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010579495 Expired - Fee Related CN102565852B (zh) | 2010-12-08 | 2010-12-08 | 针对储层含油气性检测的角度域叠前偏移数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102565852B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760603B (zh) * | 2014-01-28 | 2016-11-09 | 中国石油大学(北京) | 转换波地震数据的叠前时间偏移方法及装置 |
CN104730578A (zh) * | 2015-03-11 | 2015-06-24 | 中国电建集团北京勘测设计研究院有限公司 | 一种浅层地震单支反射数据自动解译方法和系统 |
CN106873031B (zh) * | 2017-02-15 | 2019-01-15 | 中国科学院地质与地球物理研究所 | 一种三维地震观测系统垂向分辨率定量分析评价方法 |
CN107561583A (zh) * | 2017-08-01 | 2018-01-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于高斯束叠前深度偏移的局部角度计算方法及成像方法 |
CN108415073B (zh) * | 2018-03-06 | 2019-09-13 | 中国科学院测量与地球物理研究所 | 角度域逆散射偏移成像方法及装置 |
US11644591B1 (en) * | 2021-09-28 | 2023-05-09 | China Petroleum & Chemical Corporation | Computer-implemented method and system for obtaining a subsurface stack image, subsurface angle gathers, and a subsurface velocity model, over an entire survey region having high velocity contrast geo-bodies |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO322089B1 (no) * | 2003-04-09 | 2006-08-14 | Norsar V Daglig Leder | Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder |
CN101281253B (zh) * | 2007-04-06 | 2010-06-30 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 |
US20100118652A1 (en) * | 2008-11-13 | 2010-05-13 | Jorg Friedrich Schneider | Determination of depth moveout and of residual radii of curvature in the common angle domain |
CN101900830B (zh) * | 2010-06-02 | 2012-01-11 | 西安石油大学 | 一种地震散射p-p波成像方法 |
-
2010
- 2010-12-08 CN CN 201010579495 patent/CN102565852B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102565852A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mi et al. | Horizontal resolution of multichannel analysis of surface waves | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
Tang et al. | Processing array acoustic-logging data to image near-borehole geologic structures | |
CN104216014B (zh) | 一种地震信号分频处理方法 | |
CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
CN101086535A (zh) | 表征油藏随时间演变的方法和程序 | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
CN102053263B (zh) | 调查表层结构的方法 | |
CN101738636B (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
Xia et al. | Resolution of high-frequency Rayleigh-wave data | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
CN107817526A (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
CN104360388A (zh) | 一种三维地震观测系统评价方法 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN104122582B (zh) | 利用叠加速度求取高精度地震波速度的方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN102866422B (zh) | 一种深度域地质实体模型生成方法 | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
Ning et al. | Application of multichannel analysis of passive surface waves method for fault investigation | |
Bakulin et al. | Seismic imaging of vertical array data acquired using smart DAS uphole acquisition system | |
CN102466818B (zh) | 一种利用井间地震数据对各向异性介质成像的方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131016 Termination date: 20161208 |