CN113093277B - 基于时效的隐蔽油气藏目标增强解释性处理方法 - Google Patents

基于时效的隐蔽油气藏目标增强解释性处理方法 Download PDF

Info

Publication number
CN113093277B
CN113093277B CN202110377582.5A CN202110377582A CN113093277B CN 113093277 B CN113093277 B CN 113093277B CN 202110377582 A CN202110377582 A CN 202110377582A CN 113093277 B CN113093277 B CN 113093277B
Authority
CN
China
Prior art keywords
strong reflection
target
data
seismic
filtering
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.)
Active
Application number
CN202110377582.5A
Other languages
English (en)
Other versions
CN113093277A (zh
Inventor
丁建强
张志林
魏冰
徐雷良
吕公河
李彩文
张卫华
王梅兰
刘斌
刘学勇
霍晗勇
邸志欣
朱静
马双斌
段伟伟
杨明
闫玉莎
申龙斌
曲志鹏
裴海
王月梅
何京国
徐钰
周雄波
刘明波
韩生全
和芬芬
谢金萍
张媛
王鑫
郝晓敏
赖进明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sinopec Oilfield Service Corp
Sinopec Petroleum Engineering Geophysics Co Ltd
Original Assignee
Sinopec Oilfield Service Corp
Sinopec Petroleum Engineering Geophysics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinopec Oilfield Service Corp, Sinopec Petroleum Engineering Geophysics Co Ltd filed Critical Sinopec Oilfield Service Corp
Priority to CN202110377582.5A priority Critical patent/CN113093277B/zh
Publication of CN113093277A publication Critical patent/CN113093277A/zh
Application granted granted Critical
Publication of CN113093277B publication Critical patent/CN113093277B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • G01V2210/324Filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase

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)采用倾角导向滤波、反Q滤波、中值滤波以及不同频段带通滤波相结合,突出目标体特征;(4)利用小波变换、调谐体有效地识别隐藏在强反射层中的目标体。本发明能大幅提高储层预测吻合率、断溶体单井能力以及断溶体钻遇率,为勘探目标的优选,井位优化部署,提供强有力的技术支撑。

Description

基于时效的隐蔽油气藏目标增强解释性处理方法
技术领域
本发明涉及一种基于时效的隐蔽油气藏目标增强解释性处理方法,属于地震勘探领域。
背景技术
隐蔽性油气藏在我国陆地上油气储量和产量中都占有非常重要的地位,已经成为中国油气储量增长的重要方向,碳酸盐岩油气藏是西部主要油藏类型之一,部分断控缝洞型油气藏地质目标体识别难度大,一方面是碳酸盐岩储层上奥陶统与中奥陶统波组抗差异大,对下伏目标层有很大的屏蔽作用;另一方面断裂成像不清晰,常规方法难以识别,导致断裂性质、规模以及期次演化认识不清,规模储层地震识别模式还不完善,制约了该区进一步的发展。因此,对隐蔽地质目标体进行针对性目标增强解释性处理,突显地质目标尤为重要,直接影响着层位、断裂解释甚至是后期储层预测的效果。而在现实工作中很容易忽略了这个重要环节,一方面是缺少有效的突出地质目标体的解释性处理方法,另一方面,由于工作时间节点的限制,对解释性处理重视性不够,没有给予这一环节足够的时间,因此,根据技术和工作需求进行有效地解释性处理尤为重要,为此,我们发明了一种基于时效的隐蔽油气藏目标增强解释性处理方法,解决了以上技术问题,能够根据工作安排进行有效地解释性处理,前景广阔,具有非常重要的意义。
发明内容
本发明的目的在于提供一种基于时效的隐蔽油气藏目标增强解释性处理方法,以克服以前采用单一的滤波方法,无法满足碳酸盐岩油气藏断裂、裂缝和洞穴不同储层目标的技术要求,也无法适用工作节点有长有短弹性的需求,针对性和实用性都不够,不能在技术人员之间进行大范围内推广应用。
为实现上述目的,本发明采取的技术方案是:一种基于时效的隐蔽油气藏目标增强解释性处理方法,其包括如下步骤:
(1)通过对多期处理成果进行对比分析确定最终成果数据;
(2)利用匹配追踪算法削弱成果数据中的强反射层,以输入层位为基准,根据最佳匹配准则匹配局部最优子波,对原始地震记录进行分解,提取与最优子波相一致的强反射特征,然后将其与地震信号强反射层进行消除,通过多次迭代,使残差最小,然后用原始地震记录减去匹配出的强反射记录,达到最终分离强反射层的目的;
(3)在对成果数据进行强反射分离后,采用倾角导向滤波、反Q滤波、中值滤波以及不同频段带通滤波相结合,抑制各种干扰,提高目标层信噪比和成像效果,同时消除断裂、缝洞体假象,突出目标体特征;
(4)利用小波变换、调谐体有效地识别隐藏在强反射层中的目标体。
进一步的,所述步骤(1)中选择数据信噪高、目标层信息丰富的数据作为最终成果数据,以用于解释性处理。
进一步的,所述步骤(2)的具体步骤如下:
(2.1)对强反射同相轴进行波形、能量、频率以及相位分析,求取强反射位置的地震子波参数,从而确定强反射同相轴的特征;
(2.2)建立一个超级的子波库,取一地震道数据利用复数道分析求其频率和相位作为初始值,然后在子波库中反复搜索,求最优子波,并进行数据重构,在重构地震道数据与原始数据误差最小的条件下,将最优子波的频谱作为地震道的频谱;
(2.3)利用最终优选的地震子波参数进行全区分离,然后再结合井资料以及实际勘探的静、动态油气藏参数,对分离后的数据进行后续验证分析,预测研究区有利目标。
本发明的有益效果是:本发明能够根据不同工作需求优选有效地解释性处理方法,可操作空间大,灵活适用,推广性强,前景广阔,具有非常重要的意义。该方法的完善和推广能够大幅度提高储层预测吻合率、断溶体单井能力以及断溶体钻遇率,为勘探目标的优选,井位优化部署,提供强有力的技术支撑。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
图1是研究区强反射分离前后效果图;
图2是paradigm软件中dip and azimuth参数测试对比图;
图3是paradigm软件中DSE参数测试对比图;
图4是News软件中针对断裂参数测试对比图;
图5是目标体典型剖面及频谱分析图;
图6是地震张量调谐体和小波分频记录;
图7是DSE处理前后地震剖面及频谱对比图;
图8是解释性处理前后AFE属性对比图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本发明提供的一种实施例:一种基于时效的隐蔽油气藏目标增强解释性处理方法,其包括如下步骤:
步骤1,通过对多期处理成果进行对比分析确定最终成果数据;其中选择数据信噪高、目标层信息丰富的数据作为最终成果数据,以用于解释性处理。
步骤2,利用匹配追踪算法削弱成果数据中的强反射层,以输入层位为基准,根据最佳匹配准则匹配局部最优子波,对原始地震记录进行分解,提取与最优子波相一致的强反射特征,然后将其与地震信号强反射层进行消除,通过多次迭代,使残差最小,然后用原始地震记录减去匹配出的强反射记录,达到最终分离强反射层的目的(图1至图3)。
具体做法如下:
(1)对强反射同相轴进行波形、能量、频率以及相位分析,求取强反射位置的地震子波参数,从而确定强反射同相轴的特征。
(2)建立一个超级的子波库,取一地震道数据利用复数道分析求其频率和相位作为初始值,然后在子波库中反复搜索,求最优子波,并进行数据重构,在重构地震道数据与原始数据误差最小的条件下,
将最优子波的频谱作为地震道的频谱。
(3)利用最终优选的地震子波参数进行全区分离,然后再结合井资料以及实际勘探的静、动态油气藏参数,对分离后的数据进行后续验证分析,预测研究区有利目标。
步骤3,在对成果数据进行强反射分离后,采用倾角导向滤波、反Q滤波、中值滤波以及不同频段带通滤波相结合,可以抑制各种干扰,提高目标层信噪比和成像效果,同时可以消除断裂、缝洞体假象,突出目标体特征。
(1)倾角导向滤波利用倾角及方位角的变化计算相邻道的相似性,提高地震横向信噪比,提高目标层连续性,由于其采用计算倾角方位角的方式,其结果对断层的刻画能力明显增强,有利于对断层及裂缝的后续研究,但计算周期较长,一般10G数据一个星期左右。反Q滤波能够提高缝洞体成像效果,大约需要2个小时,带通滤波可以消除高频和低频干扰,提高目标层信噪比,中值滤波可用于衰减地震数据中的相干波场,抑制线性干扰,butterworth滤波可以较好地削弱低频和高频干扰,突出目标体的发射特征,这三种方法计算速度快,一般5-10分钟。
(2)根据工作的安排,如果有足够的时间,采用倾角导向滤波、反Q滤波以及不同频段带通滤波相结合进行目标增强性处理可以有效地提高不同级别断裂和目的层缝洞的地震成像效果,如果时间紧张,需要短时间完成任务,采用中值滤波、带通滤波以及butterworth滤波突出缝洞体特征(图2至图4)。
步骤4,利用小波变换、调谐体有效地识别隐藏在强反射层中的目标体。
(1)对目标层进行频谱分析,确定地质目标有效频宽,比如奥陶系有效频宽为5-50,进行5Hz、10Hz、15Hz、20Hz、25Hz、30Hz、40Hz、50Hz小波分频分析。
(2)选取目标层时窗进行调谐体分析,并结合地震张量属性进行5Hz、10Hz、15Hz、20Hz、25Hz、30Hz、40Hz、50Hz小波分频分析(图5和图6)。
通过本发明对新疆YL地区解释性处理前后的对比,DSE处理前后地震剖面及频谱对比图,从图7上可以看出,处理前后地震剖面的特征有了明显的改善,噪音得到压制,剖面上的串珠的强振幅特征保持较好,断裂的断点干脆清晰。从频谱处理前后对比(红线为处理前结果,蓝线为处理后结果)来看,数据的振幅谱基本上没有变化,说明处理结果对于振幅谱在保幅方面没有变化,处理结果不会改变地震数据的振幅谱特征。从解释性处理前后提取的AFE属性对比分析也可以看出,解释性处理后断层和裂缝的清晰程度上有了明显的提高(图8)。
以上显示和描述了本发明的基本原理、主要特征和优点。本领域的普通技术人员应该了解,上述实施例不以任何形式限制本发明的保护范围,凡采用等同替换等方式所获得的技术方案,均落于本发明的保护范围内。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。

Claims (3)

1.一种基于时效的隐蔽油气藏目标增强解释性处理方法,其特征在于包括如下步骤:
(1)通过对多期处理成果进行对比分析确定最终成果数据;
(2)利用匹配追踪算法削弱成果数据中的强反射层,以输入层位为基准,根据最佳匹配准则匹配局部最优子波,对原始地震记录进行分解,提取与最优子波相一致的强反射特征,然后将其与地震信号强反射层进行消除,通过多次迭代,使残差最小,然后用原始地震记录减去匹配出的强反射记录,达到最终分离强反射层的目的;
(3)在对成果数据进行强反射分离后,采用倾角导向滤波、反Q滤波、中值滤波以及不同频段带通滤波相结合,抑制各种干扰,提高目标层信噪比和成像效果,同时消除断裂、缝洞体假象,突出目标体特征;
(4)利用小波变换、调谐体有效地识别隐藏在强反射层中的目标体。
2.根据权利要求1所述的一种基于时效的隐蔽油气藏目标增强解释性处理方法,其特征在于,所述步骤(1)中选择数据信噪高、目标层信息丰富的处理成果数据作为最终成果数据,以用于解释性处理。
3.根据权利要求1所述的一种基于时效的隐蔽油气藏目标增强解释性处理方法,其特征在于,所述步骤(2)的具体步骤如下:
(2.1)对强反射同相轴进行波形、能量、频率以及相位分析,求取强反射位置的地震子波参数,从而确定强反射同相轴的特征;
(2.2)建立一个超级的子波库,取一地震道数据利用复数道分析求其频率和相位作为初始值,然后在子波库中反复搜索,求最优子波,并进行数据重构,在重构地震道数据与原始数据误差最小的条件下,将最优子波的频谱作为地震道的频谱;
(2.3)利用最终优选的地震子波参数进行全区分离,然后再结合井资料以及实际勘探的静、动态油气藏参数,对分离后的数据进行后续验证分析,预测研究区有利目标。
CN202110377582.5A 2021-04-08 2021-04-08 基于时效的隐蔽油气藏目标增强解释性处理方法 Active CN113093277B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110377582.5A CN113093277B (zh) 2021-04-08 2021-04-08 基于时效的隐蔽油气藏目标增强解释性处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110377582.5A CN113093277B (zh) 2021-04-08 2021-04-08 基于时效的隐蔽油气藏目标增强解释性处理方法

Publications (2)

Publication Number Publication Date
CN113093277A CN113093277A (zh) 2021-07-09
CN113093277B true CN113093277B (zh) 2022-06-24

Family

ID=76675092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110377582.5A Active CN113093277B (zh) 2021-04-08 2021-04-08 基于时效的隐蔽油气藏目标增强解释性处理方法

Country Status (1)

Country Link
CN (1) CN113093277B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115877460B (zh) * 2023-02-28 2023-05-05 福瑞升(成都)科技有限公司 一种碳酸盐岩岩溶缝洞型储层的增强方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6789018B1 (en) * 2003-08-29 2004-09-07 Nonlinear Seismic Imaging, Inc. Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave
CN103527184A (zh) * 2013-10-28 2014-01-22 北京大学 一种白云岩储层的预测方法和系统
CN109507732A (zh) * 2018-11-07 2019-03-22 中石化石油工程技术服务有限公司 一种基于成像道集的绕射波分离成像方法
CN110794461A (zh) * 2018-08-01 2020-02-14 中国石油化工股份有限公司 一种强地震反射界面下的小尺度缝洞信息凸显方法
CN112180447A (zh) * 2019-07-04 2021-01-05 中国石油天然气集团有限公司 一种消除储层强反射屏蔽的方法及系统
CN112346121A (zh) * 2020-11-04 2021-02-09 东北石油大学 一种基于全波形的储集体地层分离处理方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6789018B1 (en) * 2003-08-29 2004-09-07 Nonlinear Seismic Imaging, Inc. Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave
CN103527184A (zh) * 2013-10-28 2014-01-22 北京大学 一种白云岩储层的预测方法和系统
CN110794461A (zh) * 2018-08-01 2020-02-14 中国石油化工股份有限公司 一种强地震反射界面下的小尺度缝洞信息凸显方法
CN109507732A (zh) * 2018-11-07 2019-03-22 中石化石油工程技术服务有限公司 一种基于成像道集的绕射波分离成像方法
CN112180447A (zh) * 2019-07-04 2021-01-05 中国石油天然气集团有限公司 一种消除储层强反射屏蔽的方法及系统
CN112346121A (zh) * 2020-11-04 2021-02-09 东北石油大学 一种基于全波形的储集体地层分离处理方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
卡尔曼最优估计下独立变量分析地震反射系数提取方法;刘菲菲等;《地球物理学进展》;20161215(第06期);全文 *
地震属性分析技术在储层预测中的应用;刘晶等;《长江大学学报(自科版)》;20150115(第02期);全文 *
海相碳酸盐岩储层地震预测技术进展及应用实效;常少英等;《海相油气地质》;20200325(第01期);全文 *

Also Published As

Publication number Publication date
CN113093277A (zh) 2021-07-09

Similar Documents

Publication Publication Date Title
CN107991706B (zh) 基于小波包多重阈值和改进经验模态分解的煤层水力压裂微震信号联合降噪方法
Wu et al. Noise attenuation for 2-D seismic data by radial-trace time-frequency peak filtering
CN111158044A (zh) 潜山裂缝储集体油藏预测方法
CN113093277B (zh) 基于时效的隐蔽油气藏目标增强解释性处理方法
CN111007569A (zh) 一种集成经验模态分解的样本熵阈值微地震信号降噪方法
CN113887398A (zh) 一种基于变分模态分解和奇异谱分析的gpr信号去噪方法
Li et al. Experimental analysis and application of sparsity constrained deconvolution
CN110333533B (zh) 基于匹配追踪和rgb色彩融合的亮点油气检测方法及系统
CN111257933B (zh) 基于低频阴影现象的油气藏预测新方法
Xudong et al. Pre-stack gather optimization technology based on an improved bidimensional empirical mode decomposition method
Yin et al. Development and application of a novel hybrid CEEMDAN-LWT denoising approach in microseismic monitoring
Jiang et al. Hydrocarbon detection based on empirical mode decomposition, teager-kaiser energy, and the cepstrum
CN110749923A (zh) 一种基于范数方程提高分辨率的反褶积方法
CN112130209B (zh) 岩溶储层预测方法及装置
CN112346121B (zh) 一种基于全波形的储集体地层分离处理方法
Wang et al. Application of multichannel Wiener filters to the suppression of ambient seismic noise in passive seismic arrays
CN113702821A (zh) 一种提取gis局部放电信号的方法及系统
CN114428280A (zh) 一种地震数据低频信息补偿方法及其应用
CN114966828A (zh) 一种断溶体储层地震反射特征增强方法
CN111290020A (zh) 基于构造滤波处理、分频属性融合断裂检测方法及装置
CN112200069A (zh) 时频域谱减法和经验模态分解联合的隧道滤波方法及系统
CN112444866B (zh) 提高超深层地震数据分辨率的方法及存储介质
CN117724153B (zh) 一种基于多窗口级联交互的岩性识别方法
Peng et al. A nonparametric method for automatic denoising of microseismic data
Ma et al. A Global and Multi-Scale Denoising Method Based on Generative Adversarial Network for DAS VSP Data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant