CN106405647B - 一种关于沉积地层厚度的调谐反演方法 - Google Patents

一种关于沉积地层厚度的调谐反演方法 Download PDF

Info

Publication number
CN106405647B
CN106405647B CN201611195830.XA CN201611195830A CN106405647B CN 106405647 B CN106405647 B CN 106405647B CN 201611195830 A CN201611195830 A CN 201611195830A CN 106405647 B CN106405647 B CN 106405647B
Authority
CN
China
Prior art keywords
thickness
frequency
frequency dividing
data
seismic
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
CN201611195830.XA
Other languages
English (en)
Other versions
CN106405647A (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.)
Chengdu Jingshi Petroleum Science & Technology Co Ltd
Original Assignee
Chengdu Jingshi Petroleum Science & Technology 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 Chengdu Jingshi Petroleum Science & Technology Co Ltd filed Critical Chengdu Jingshi Petroleum Science & Technology Co Ltd
Priority to CN201611195830.XA priority Critical patent/CN106405647B/zh
Publication of CN106405647A publication Critical patent/CN106405647A/zh
Application granted granted Critical
Publication of CN106405647B publication Critical patent/CN106405647B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

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:分析地震数据的主频,并根据地震频谱中提取子波,计算对应主频下的波长,λ=1/f,其中λ为波长,f为频率;步骤2:建立一个地层厚度由0到λ的地层模型,将提取子波进行负90度相位旋转后,对地层模型进行褶积运算,得到正演的不同厚度对应的薄层响应原子;步骤3:在地震有效频带范围内,利用分频工具将正演的不同厚度薄层响应原子进行分频获得N个分频体,在N个体内分别计算不同厚度下的最大振幅值,构建不同频率下最大振幅与地层厚度的一个响应曲线,分别对其做归一化处理,可以得到分频振幅与厚度的模板曲线,N为大于1的自然数。

Description

一种关于沉积地层厚度的调谐反演方法
技术领域
本发明涉及石油勘探领域,具体涉及一种关于沉积地层厚度的调谐反演方法。
背景技术
从寻找石油到利用石油,大致要经过四个主要环节,即寻找、开采、输送和加工,这四个环节一般又分别称为“石油勘探”、“油田开发”、“油气集输”和“石油炼制”。“石油勘探”有许多方法,但地下是否有油,最终要靠钻井来证实。一个国家在钻井技术上的进步程度,往往反映了这个国家石油工业的发展状况,因此,有的国家竞相宣布本国钻了世界上第一口油井,以表示他们在石油工业发展上迈出了最早的一步。“油田开发”指的是用钻井的办法证实了油气的分布范围,并且油井可以投入生产而形成一定生产规模。从这个意义上说,1821年四川富顺县自流井气田的开发是世界上最早的天然气田。“油气集输”技术也随着油气的开发应运而生,公元1875年左右,自流井气田采用当地盛产的竹子为原料,去节打通,外用麻布缠绕涂以桐油,连接成我们现在称呼的“输气管道”,总长二、三百里,在当时的自流井地区,绵延交织的管线翻越丘陵,穿过沟涧,形成输气网络,使天然气的应用从井的附近延伸到远距离的盐灶,推动了气田的开发,使当时的天然气达到年产7000多万立方米。
开采石油的第一关是勘探油田。今天的石油地质学家使用重力仪、磁力仪等仪器来寻找新的石油储藏。储层厚度是油气资源评价的重要参数。目前仅通过反演的方法提取储层厚度,但受限于初始模型及,时空变的压实关系,很难得到准确的厚度预测结果。
发明内容
本发明所要解决的技术问题是明确对厚度敏感的振幅优势区间,在不同的区间内分别预测厚度,目的在于提供一种关于沉积地层厚度的调谐反演方法,解决明确对厚度敏感的振幅优势区间,在不同的区间内分别预测厚度的问题。
本发明通过下述技术方案实现:
一种关于沉积地层厚度的调谐反演方法,包括以下步骤,步骤1:分析地震数据的主频,并根据地震频谱中提取子波,计算对应主频下的波长,λ=1/f,其中λ为波长,f为频率;步骤2:建立一个地层厚度由0到λ的地层模型,将提取子波进行负90度相位旋转后,对地层模型进行褶积运算,得到正演的不同厚度对应的薄层响应原子;步骤3:在地震有效频带范围内,利用分频工具将正演的不同厚度薄层响应原子进行分频获得N个分频体,在N个体内分别计算不同厚度下的最大振幅值,构建不同频率下最大振幅与地层厚度的一个响应曲线,分别对其做归一化处理,可以得到分频振幅与厚度的模板曲线,N为大于1的自然数;步骤4:在步骤3不同频率模板曲线的调谐振幅以内,利用多项式拟合的方式,可获得不同频率下振幅与厚度的确定性关系,因为不同频率模板曲线,都存在一个调谐振幅,归一化后为1,在调谐振幅以内,地层厚度与振幅值成单调正比关系,并且获得不同频率下振幅与厚度的确定性关系后,高频曲线可预测厚度范围小,集中在薄层;低频曲线可预测范围大,集中在厚层;步骤5:将地震数据进行负90度相位转换,使得波峰、波谷信息代表一个地层信息,利用分频工具对地震数据进行分频,获得N个分频数据体;并统计某一个层段内最大振幅,是根据层位提取某一段地震数据来统计,利用分频最大振幅对数据体进行归一化计算;步骤6:将N个分频归一化数据体,利用模板曲线的多项式关系将振幅转化为厚度,得到分频厚度数据体;步骤7:每个分频体的最大振幅位置,即这个分频体计算的厚度最为可信位置,所以将步骤6中得到的N个分频厚度数据体,利用反距离加权方法,计算出加权系数,以地震体中的一道为计算单元,将分频厚度数据体进行加权平均,获得地震厚度体。地层是由一系列薄层反射系数对组成,地震数据体相当于是一簇薄层反射系数与子波褶积的结果。这种效应相当于子波对不同砂层厚度进行了一次滤波。当子波主频与薄层厚度匹配时,地震波形振幅会被加强,发生调谐,显示最大振幅效应,对应子波的主频称为调谐频率,最大振幅称为调谐振幅。在调谐振幅以内,振幅随厚度是成正比关系的,因此我们可以利用这种关系进行厚度的预测。通常情况下,地震主频为地震的一个调谐频率,常规厚度预测以主频的计算的1/4波长为地震极限分辨率。利用一个固定子波对不同厚度地层进行滤波形成薄层响应原子序列,经过分频之后的薄层响应原子序列也存在调谐振幅与厚度的关系,且不同分频下振幅可预测范围不一样。低频调谐振幅可预测更厚地层,同时高频振幅曲线在薄层区梯度更大,对薄地层厚度更为敏感。因为我们可以利用多个分频体来综合预测地层厚度。本发明以理论模型正演为指导,构造薄层响应模型,制作相应分频振幅与与厚度关系的响应模板,地震进行小波分频并将其归一化后后,利用模板关系,将分频地震振幅转化为厚度分析不同频率下振幅的变化,不同频率下振幅对厚度敏感区间不同,将多个频率预测的厚度体进行加权叠加,提高厚度预测范围和预测精度;本发明预测厚度完全依赖于地震数据的分频振幅,是一种厚度属性数据体,不同于用测井数据、初始模型进行波阻抗反演,然后卡门槛的方式进行厚度预测,是一种简便有效,同时也是相对精确的厚度预测方法,解决了无井地区无法进行厚度预测的难题。
所述步骤3中的分频工具包括小波分频、S变换或匹配追踪。三种分频方式都可以在地震有效频带范围内,将正演的不同厚度薄层响应原子进行分频。
所述分频工具采用小波分频。小波分频工具是时频分析的一种方式,S变换及匹配追踪也可用进行尝试,但是前者输出时间域分频,后者是输出频率域分频,时间域分频方式稳定性和抗噪性更好。
所述步骤5中利用分频工具对地震数据进行分频,获得3个分频数据体。单一的频率所预测的厚度范围有限,高频或低频的数据体一方面扩大可预测厚度范围,同时增加了预测厚度的稳定性,数据体为低频、中频、高频数据体,本发明的数据体个数是不限制的,一般采用低频、中频、高频三个数据体,这选取中间值、最高值和最低值,是比较合理的范围
本发明与现有技术相比,具有如下的优点和有益效果:
1、本发明一种关于沉积地层厚度的调谐反演方法,步骤1、2、3、4是利用模型正演制作模板,步骤5、6、7将地震数据转化为地层厚度,完全不用测井数据、初始模型就可以达到厚度预测目的;
2、本发明一种关于沉积地层厚度的调谐反演方法,最终成果就是厚度剖面,不需要通过常规有井反演的方式,通过岩石物理门槛值求取和目的层时窗选取等繁琐工作来获得厚度数据体;
3、本发明一种关于沉积地层厚度的调谐反演方法,获得厚度数据体,相当于一种厚度属性剖面,有利于直接反映沉积厚度变化情况,更有利于地震沉积学切片分析技术在石油勘探领域的推广。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明所述步骤4中调谐振幅归一化前示意图;
图2为本发明所述步骤4中调谐振幅归一后前示意图;
图3为本发明所述步骤5中厚度预测优势区间示意图;
图4为本发明分频算厚度示意图;
图5为本发明模型厚度对比图示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例
如图1至图5所示,本发明一种关于沉积地层厚度的调谐反演方法,以模型为算例,描述厚度预测的过程,包括以下步骤:
步骤1:给定一个30Hz主频的雷克子波,计算主频对应波长λ;
步骤2:建立一个地层厚度由0到λ的地层模型,将雷克子波进行负90度相位旋转,利对地层模型进行褶积运算,正演的不同厚度对应的薄层响应原子;
步骤3:在地震有效频带范围内,利用小波分频工具将正演的不同厚度薄层响应原子进行分频获得3个分频体10Hz、20Hz、40Hz,在3个体内分别计算不同厚度下的最大振幅值如图1,构建不同频率下最大振幅与地层厚度的一个响应曲线,分别对其做归一化处理,不考虑分频之后振幅的差异,构成分频振幅与厚度的模板曲线;
步骤4:不同频率模板曲线,都存在一个调谐振幅,归一化后为1,,在调谐振幅以内,地层厚度与振幅值成单调正比关系,利用多项式拟合的方式,可获得不同频率下振幅与厚度的确定性关系,40Hz高频曲线可预测厚度范围小,集中在薄层,10Hz低频曲线可预测范围大,集中在厚层;
步骤5:将地震数据进行-90度相位转换,使得波峰、波谷信息代表一个地层信息,利用小波分频工具对地震数据进行分频,获得10Hz、20Hz、40Hz分频数据体;统计某一个层段内最大振幅,利用分频最大振幅对数据体进行归一化计算;
步骤6:将3个分频归一化数据体,利用模板曲线的多项式关系将振幅转化为厚度,形成分频厚度数据体;
步骤7:利用反距离加权方法,将每个分频厚度体的加权系数计算出来,哪个分频体最大振幅即这个分频体计算的厚度最为可信,计算权重最大。以地震体中的一道为计算单元,将分频厚度数据体进行加权平均,获得最终的地震厚度体。10Hz、20Hz和40Hz的分频体能准确预测19ms以内地层厚度,预测厚度与实际模型厚度计算误差不大。10Hz分频体预测厚度可预测分辨率最大能达到19ms,因此要分辨大于19ms的厚层,需更低频率的分频体来计算。采用本发明所述的一种调谐厚度反演方法特点:不用测井数据、初始模型就可以达到厚度预测目的;本发明最终成果在剖面上直接显示厚度,不需要进行岩石物理门槛值求取和目的层时窗选取等繁琐工作;本发明最终成果能在剖面上直接反映沉积厚度变化情况,更有利于地震沉积学切片分析技术在石油勘探领域的推广。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种关于沉积地层厚度的调谐反演方法,其特征在于:包括以下步骤,
步骤1:分析地震数据的主频,并根据地震频谱中提取子波,计算对应主频下的波长,λ=1/f,其中λ为波长,f为频率;
步骤2:建立一个地层厚度由0到λ的地层模型,将提取子波进行负90度相位旋转后,对地层模型进行褶积运算,得到正演的不同厚度对应的薄层响应原子;
步骤3:在地震有效频带范围内,利用分频工具将正演的不同厚度薄层响应原子进行分频获得N个分频体,在N个体内分别计算不同厚度下的最大振幅值,构建不同频率下最大振幅与地层厚度的一个响应曲线,分别对其做归一化处理,得到分频振幅与厚度的模板曲线,N为大于1的自然数;
步骤4:在步骤3不同频率模板曲线的调谐振幅以内,利用多项式拟合的方式,可获得不同频率下振幅与厚度的确定性关系;
步骤5:将地震数据进行负90度相位转换,使得波峰、波谷信息代表一个地层信息,利用分频工具对地震数据进行分频,获得N个分频数据体;并统计某一个层段内最大振幅,利用分频最大振幅对数据体进行归一化计算;
步骤6:将N个分频归一化数据体,利用模板曲线的多项式关系将振幅转化为厚度,得到分频厚度数据体;
步骤7:将步骤6中得到的N个分频厚度数据体,利用反距离加权方法,计算出加权系数,以地震体中的一道为计算单元,将分频厚度数据体进行加权平均,获得地震厚度体。
2.根据权利要求1所述的一种关于沉积地层厚度的调谐反演方法,其特征在于:所述步骤3中的分频工具包括小波分频、S变换或匹配追踪。
3.根据权利要求2所述的一种关于沉积地层厚度的调谐反演方法,其特征在于:所述分频工具采用小波分频。
4.根据权利要求1所述的一种关于沉积地层厚度的调谐反演方法,其特征在于:所述步骤5中利用分频工具对地震数据进行分频,获得3个分频数据体。
CN201611195830.XA 2016-12-22 2016-12-22 一种关于沉积地层厚度的调谐反演方法 Active CN106405647B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611195830.XA CN106405647B (zh) 2016-12-22 2016-12-22 一种关于沉积地层厚度的调谐反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611195830.XA CN106405647B (zh) 2016-12-22 2016-12-22 一种关于沉积地层厚度的调谐反演方法

Publications (2)

Publication Number Publication Date
CN106405647A CN106405647A (zh) 2017-02-15
CN106405647B true CN106405647B (zh) 2018-04-13

Family

ID=58087917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611195830.XA Active CN106405647B (zh) 2016-12-22 2016-12-22 一种关于沉积地层厚度的调谐反演方法

Country Status (1)

Country Link
CN (1) CN106405647B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109655904A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 一种多尺度溶洞识别方法及系统
CN110907994B (zh) * 2018-09-18 2021-09-17 中国石油化工股份有限公司 一种基于分频的波形反演储层预测方法
CN111221038B (zh) * 2018-11-26 2022-02-01 中国石油天然气股份有限公司 薄储层厚度定量预测的方法和装置
CN110568490B (zh) * 2019-08-21 2021-09-17 中国石油化工股份有限公司 一种高速层顶薄储层的识别方法
CN113419283B (zh) * 2020-05-28 2023-04-07 中国海洋石油集团有限公司 一种基于分频目标迭代体的调谐振幅反演方法
CN114019566B (zh) * 2021-11-04 2024-06-14 北京诺克斯达石油科技有限公司 连续面块地层切片薄互层预测方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1210591A (zh) * 1996-12-06 1999-03-10 阿莫科公司 用于地震分析的频谱解译
US6131071A (en) * 1996-12-06 2000-10-10 Bp Amoco Corporation Spectral decomposition for seismic interpretation
CN103412332A (zh) * 2013-01-22 2013-11-27 中国地质大学(北京) 一种定量计算薄储层厚度的方法
CN103558635A (zh) * 2013-10-30 2014-02-05 北京诺克斯达石油科技有限公司 基于偶函数地震响应以估算薄层厚度的方法及装置
CN104297791A (zh) * 2014-09-25 2015-01-21 中国石油天然气股份有限公司 一种基于地震优势频率的反演方法及系统
CN105093306A (zh) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 一种地球物理勘探中储层自动解释与厚度求取方法
CN105092343A (zh) * 2014-05-04 2015-11-25 中国石油化工股份有限公司 一种基于叠前道集去除薄层调谐效应的方法
WO2016168373A1 (en) * 2015-04-17 2016-10-20 Schlumberger Technology Corporation Thin bed tuning frequency and thickness estimation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1210591A (zh) * 1996-12-06 1999-03-10 阿莫科公司 用于地震分析的频谱解译
US6131071A (en) * 1996-12-06 2000-10-10 Bp Amoco Corporation Spectral decomposition for seismic interpretation
CN103412332A (zh) * 2013-01-22 2013-11-27 中国地质大学(北京) 一种定量计算薄储层厚度的方法
CN103558635A (zh) * 2013-10-30 2014-02-05 北京诺克斯达石油科技有限公司 基于偶函数地震响应以估算薄层厚度的方法及装置
CN105092343A (zh) * 2014-05-04 2015-11-25 中国石油化工股份有限公司 一种基于叠前道集去除薄层调谐效应的方法
CN105093306A (zh) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 一种地球物理勘探中储层自动解释与厚度求取方法
CN104297791A (zh) * 2014-09-25 2015-01-21 中国石油天然气股份有限公司 一种基于地震优势频率的反演方法及系统
WO2016168373A1 (en) * 2015-04-17 2016-10-20 Schlumberger Technology Corporation Thin bed tuning frequency and thickness estimation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
分频反演方法及应用;于建国 等;《石油地球物理勘探》;20060430;第41卷(第2期);第193-197页 *
基于调谐频率与分频处理的高分辨率反演技术;杨贵详;《石油物探》;20060531;第45卷(第3期);第242-244页 *

Also Published As

Publication number Publication date
CN106405647A (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN106405647B (zh) 一种关于沉积地层厚度的调谐反演方法
US10942293B2 (en) Rock physics based method of integrated subsurface reservoir characterization for use in optimized stimulation design of horizontal wells
CN106154323B (zh) 基于地震拓频处理的相控随机反演薄储层预测方法
Bosch et al. Petrophysical seismic inversion conditioned to well-log data: Methods and application to a gas reservoir
CN104502997B (zh) 一种利用裂缝密度曲线预测裂缝密度体的方法
US11181653B2 (en) Reservoir characterization utilizing ReSampled seismic data
CN106291677B (zh) 一种基于匹配追踪方法的叠后声波阻抗反演方法
CN104122581B (zh) 一种叠后声波阻抗反演方法
CN103439740B (zh) 基于偶极地震子波多重积分的相对阻抗预测的方法及装置
CN105738952B (zh) 一种水平井区储层岩石相建模方法
CN107219553A (zh) 基于gr分频反演的暗河充填预测方法
CN106443770A (zh) 一种页岩气地质甜点的预测方法
CN108508489B (zh) 一种基于波形微变化匹配的地震反演方法
Artun et al. Intelligent seismic inversion workflow for high-resolution reservoir characterization
Qifeng et al. Logging-while-drilling formation dip interpretation based on long short-term memory
Chitralekha et al. Application of the ensemble Kalman filter for characterization and history matching of unconventional oil reservoirs
Dawuda et al. Geologic modeling and ensemble-based history matching for evaluating co2 sequestration potential in point bar reservoirs
CN105785443B (zh) 高精度道积分计算相对波阻抗的方法
Semenov et al. Application of group method of data handling for geological modeling of vankor field
US11320565B2 (en) Petrophysical field evaluation using self-organized map
CN114075973B (zh) 地层元素测井曲线重构方法及装置
CN107703544B (zh) 基于地质统计学的叠前振幅随偏移距变化油气预测方法
Worthington Maximizing the effectiveness of integrated reservoir studies: practical approaches to improving the process and results
Huang et al. New seismic attribute: Fractal scaling exponent based on gray detrended fluctuation analysis
CN113640880B (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
GR01 Patent grant
GR01 Patent grant