CN113777657B - 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程 - Google Patents

基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程 Download PDF

Info

Publication number
CN113777657B
CN113777657B CN202111047129.4A CN202111047129A CN113777657B CN 113777657 B CN113777657 B CN 113777657B CN 202111047129 A CN202111047129 A CN 202111047129A CN 113777657 B CN113777657 B CN 113777657B
Authority
CN
China
Prior art keywords
well
compaction
seismic
stratum
section
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
CN202111047129.4A
Other languages
English (en)
Other versions
CN113777657A (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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202111047129.4A priority Critical patent/CN113777657B/zh
Publication of CN113777657A publication Critical patent/CN113777657A/zh
Application granted granted Critical
Publication of CN113777657B publication Critical patent/CN113777657B/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/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • 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/301Analysis for determining seismic cross-sections or geostructures

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

本发明公开了基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程,具体方法流程是,合成记录标定确定时深关系,进行目的层顶底界面的地震解释;利用砂岩和泥岩的压实系数进行井点压实恢复,在单井压实厚度恢复的基础上,通过剖面多井建立地层压实综合恢复系数与深度关系,获得连井剖面上压实恢复的地层厚度;通过时深转换,得到时间域的压实恢复后地层厚度,利用其对地震道数据进行拉伸,进行插值、重采样处理,在处理后的地震剖面上进行地震相解释,得到差异压实恢复后的地震相。本发明为斜坡带差异压实背景下的地震相准确识别提供了一种简单、合理而有效的方法。

Description

基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程
技术领域
本发明涉及油气藏勘探开发技术领域,具体的说,是一种压实校正下的斜坡带连井剖面地震相表征的工业化流程。
背景技术
地震相是指在地震剖面上特定的地震反射及其包络,它是20世纪70年代建立起来的一套新的概念和地震资料地质解释方法,是地震地层学中的核心概念之一。地震相的提出以及常见的地震相类型的总结,建立了利用地震资料解释沉积体系的新方法和新思路,为后续发展起来的层序地层学提供了重要方法支撑,并且一直沿用至今。地震相分析将地震剖面上特定的规律性反射与沉积体系之间建立关系,为沉积体系解释提供了可靠的方法。
在实际应用中发现,地震相受到地震资料质量、地震分辨率、地层的岩石物理性质等多种因素的影响。然而,有一个重要因素的影响一直没有受到足够的关注,这个因素就是沉积期古地貌的影响,对于非剥蚀区来说,最大的影响来自于差异压实造成的局部地层倾角变化。
地震相分析中,将地震反射同相轴作为等时沉积界面,通过分析这些等时界面的延伸以及相互关系来识别和表征沉积体系及其内部结构。比如,对上超、前积等典型反射结构的识别帮助认识地层的沉积建造。在地震相分析过程中,等时沉积界面的倾角关系、相邻等时界面之间地层的厚度及其横向变化等因素都是识别沉积体系及其展布的重要依据。而上述重要因素除了与地层沉积本身相关外,还受到地层沉积之后差异压实等因素的影响。比如,一些大型的前积反射结构,受到后期的差异压实影响,倾角减小,反射变得平缓,已经难以识别出前积反射的地震相结构,造成前积反射不发育的地质假象。这种情况在斜坡区的影响尤其显著。真正能够指导沉积体系表征的,是沉积期的地层特征,因此,对沉积期古地貌进行恢复,在此基础上进行地震相分析,才能更好的反映真实的沉积体系及其内部结构。
差异压实校正在基于井资料的古地貌恢复中属于常规的流程性操作,但是在开展地震相研究中,压实差异性对地震相的影响没有引起关注,更没有相应的方法被引入地震相表征。
在井点少的情况下,如在海上油田或者陆上油田勘探及开发早期,大部分情况下,在地震相研究中并不考虑古地貌的影响,即使部分研究考虑了这一点,也都是进行简单的地震剖面层拉平,没有考虑压实差异对地震反射形态的变化引起的地震相解释不准确问题。本发明对斜坡区差异压实引起的古地貌变化进行恢复,然后进行地震相识别,形成一个新的基于差异压实恢复的地震相表征工业化流程。
发明内容
本发明提供了一种基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程,该方法具体的技术方案包括如下步骤:
步骤(1):在通过n口井的连井剖面上,n为自然数,逐口井制作合成地震记录进行井震标定,利用单井标定结果拟合研究区目的层段时深关系,将拟合的时深关系式记作F0
步骤(2):根据井震标定结果,在地震剖面上追踪解释目的层段的顶界面和底界面,将顶界面记为界面T,底界面记为界面B;
步骤(3):利用界面B和T,得到目的层时间域地层厚度,记为Ht
步骤(4):选取连井剖面上所有井的典型砂岩段和典型泥岩段,采用孔隙度变化拟合法,计算研究区的砂岩压实系数Cs和泥岩压实系数Cm;
步骤(5):选取连井剖面上的第i口井,i为不大于n的自然数,该井名称记作Wi,利用Wi井的录井岩性数据和步骤(5)得到的Cs及Cm,采用回剥法进行单井压实校正,得到Wi井处目的层的压实恢复沉积厚度Hwir
步骤(6):Wi井的目的层地层厚度记作Hwi,将Hwi除以步骤(5)求得的Hwir,得到Wi井处的地层压实综合恢复系数,记作Tci
步骤(7):按照步骤(5)和步骤(6)的方法,计算每口井处目的层段地层压实综合恢复系数;
步骤(8):利用步骤(7)得到的每口井处地层压实综合恢复系数与对应井点处的目的层顶面深度,建立深度与地层压实综合恢复系数的拟合关系式,将关系式记作F;
步骤(9):利用公式F与步骤(1)得到的时深关系F0,得到时间域的目的层顶面地震反射时间与地层压实综合恢复系数的关系式,将关系式记作F2
步骤(10):利用步骤(2)得到的目的层顶界面T和步骤(9)得到的关系式F2,获取连井剖面上沿着目的层顶面每个位置的地层压实综合恢复系数;
步骤(11):利用步骤(3)得到的时间域地层厚度Ht与步骤(10)得到的地层压实综合恢复系数,计算得到连井剖面上的时间域地层厚度压实恢复结果Tf
步骤(12):读取连井剖面上目的层段的地震道数据,按照步骤(11)得到的时间域地层厚度恢复结果Tf对剖面上每个地震道的长度进行拉伸,对拉伸后的目的层段地震道数据进行重新插值和采样;
步骤(13):对步骤(12)得到的处理后的目的层段地震数据进行剖面波形显示,在新的地震剖面上,利用传统的地震相模式,进行地震相识别与解释。
具体实施方式
按照上述发明内容中的步骤实施,具体包括以下步骤:
步骤(1):对连井剖面上的每一口井制作合成地震记录,进行井震标定,利用单井标定结果拟合研究区目的层段时深关系,将拟合的时深关系式记作F0
步骤(2):根据井震标定结果,在地震剖面上追踪解释目的层段的顶界面和底界面,将顶界面记为界面T,底界面记为界面B;
步骤(3):利用界面B和T,得到目的层时间域地层厚度,记为Ht
步骤(4):选取连井剖面上所有井的典型砂岩段和典型泥岩段,采用孔隙度变化拟合法,计算研究区的砂岩压实系数Cs和泥岩压实系数Cm;
步骤(5):选取连井剖面上的第i口井,i为不大于n的自然数,该井名称记作Wi,利用Wi井的录井岩性数据和步骤(5)得到的Cs及Cm,采用回剥法进行单井压实校正,得到Wi井处目的层的压实恢复沉积厚度Hwir
步骤(6):Wi井的目的层地层厚度记作Hwi,将Hwi除以步骤(5)求得的Hwir,得到Wi井处的地层压实综合恢复系数Tci
步骤(7):按照步骤(5)和步骤(6)的方法,计算每口井处目的层段地层压实综合恢复系数;
步骤(8):利用步骤(7)得到的每口井处地层压实综合恢复系数与对应井点处的目的层顶面深度,建立深度与地层压实综合恢复系数的拟合关系式,将关系式记作F;
步骤(9):利用公式F与步骤(1)得到的时深关系F0,得到时间域的目的层顶面地震反射时间与地层压实综合恢复系数的关系式,将关系式记作F2
步骤(10):利用步骤(2)得到的目的层顶界面T和步骤(9)得到的关系式F2,获取连井剖面上沿着目的层顶面每个位置的地层压实综合恢复系数;
步骤(11):利用步骤(3)得到的时间域地层厚度Ht与步骤(10)得到的地层压实综合恢复系数,计算得到连井剖面上的时间域地层厚度压实恢复结果Tf
步骤(12):读取连井剖面上目的层段的地震道数据,按照步骤(11)得到的时间域地层厚度恢复结果Tf对剖面上每个地震道的长度进行拉伸,对拉伸后的目的层段地震道数据进行重新插值和采样;
步骤(13):对步骤(12)得到的处理后的目的层段地震数据进行剖面波形显示,在新的地震剖面上,利用传统的地震相模式,进行地震相识别与解释。
实施例
本实施例的研究区位于辽河油田,位于辽河东部凹陷斜坡带,目的层段为古近系沙河街组三段(以下简称沙三段)地层,选取自西北向东南方向的连井剖面,该剖面穿过4口井,名称依次记作W1,W2,W3,W4
步骤(1):分别制作W1,W2,W3,W4四口井的合成地震记录,进行井震标定,利用单井标定结果拟合研究区目的层段时深关系,将拟合的时深关系式为Dep=1.429Tim-392.1,其中Dep为深度,单位m,Tim为地震反射时间,单位为ms;
步骤(2):根据井震标定结果,在地震剖面上追踪解释目的层段的顶界面和底界面,将顶界面记为界面T,底界面记为界面B;
步骤(3):将界面B减去T得到目的层时间域地层厚度,记为Ht
步骤(4):在四口井上选取6个典型砂岩段和4个典型泥岩段,采用孔隙度变化拟合法,计算研究区的砂岩压实系数为0.039和泥岩压实系数为0.056;
步骤(5):利用W1井的录井岩性数据、砂岩压实系数、泥岩压实系数,采用回剥法进行单井压实校正,得到W1井处目的层的压实恢复沉积厚度125.42m;
步骤(6):W1井目的层地层厚度为120m,将其除以步骤(5)求得的压实恢复厚度125.42m,得到W1井处的地层压实综合恢复系数0.957;
步骤(7):按照步骤(5)和步骤(6)的方法,计算每口井处目的层段地层压实综合恢复系数;
步骤(8):利用步骤(7)得到的每口井处地层压实综合恢复系数与对应井点处的目的层顶面深度,建立深度与地层压实综合恢复系数的拟合关系式,Rh=-2.67×10-5Dep+0.493,其中,Dep为深度,单位m,Rh为地层压实综合恢复系数;
步骤(9):利用步骤(8)得到的计算公式与步骤(1)得到的时深关系,得到时间域的目的层顶面地震反射时间与地层压实综合恢复系数的关系式,Rh=-3.815×10-5Tim+0.503,其中,Tim为地震反射时间,单位为ms;
步骤(10):利用步骤(2)得到的目的层顶界面T和步骤(9)得到的关系式,获取连井剖面上沿着目的层顶面每个位置的地层压实综合恢复系数;
步骤(11):利用步骤(3)得到的时间域地层厚度Ht与步骤(10)得到的地层压实综合恢复系数,计算得到连井剖面上的时间域地层厚度压实恢复结果Tf
步骤(12):读取连井剖面上目的层段的地震道数据,按照步骤(11)得到的时间域地层厚度恢复结果Tf对剖面上每个地震道的长度进行拉伸,对拉伸后的目的层段地震道数据进行重新插值和采样;
步骤(13):对步骤(12)得到的处理后的目的层段地震数据进行剖面波形显示,在新的地震剖面上,利用传统的地震相模式,进行地震相识别与解释,在剖面上识别出了中连续低角度地震前积相以及代表河道沉积的强振幅弱连续相等四种地震相类型。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (1)

1.基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程,其特征在于,包括以下步骤:
步骤(1):在通过n口井的连井剖面上,n为自然数,逐口井制作合成地震记录进行井震标定,利用单井标定结果拟合研究区目的层段时深关系,将拟合的时深关系式记作F0
步骤(2):根据井震标定结果,在地震剖面上追踪解释目的层段的顶界面和底界面,将顶界面记为界面T,底界面记为界面B;
步骤(3):利用界面B和T解释结果,得到目的层时间域地层厚度,记为Ht
步骤(4):选取连井剖面上所有井的典型砂岩段和典型泥岩段,采用孔隙度变化拟合法,计算研究区的砂岩压实系数Cs和泥岩压实系数Cm;
步骤(5):选取连井剖面上的第i口井,i为不大于n的自然数,该井名称记作Wi,利用Wi井的录井岩性数据和步骤(5)得到的Cs及Cm,采用回剥法进行单井压实校正,得到Wi井处目的层的压实恢复沉积厚度Hwir
步骤(6):Wi井的目的层地层厚度记作Hwi,将Hwi除以步骤(5)求得的Hwir,得到Wi井处的地层压实综合恢复系数,记作Tci
步骤(7):按照步骤(5)和步骤(6)的方法,计算每口井处目的层段地层压实综合恢复系数;
步骤(8):利用步骤(7)得到的每口井处地层压实综合恢复系数与对应井点处的目的层顶面深度,建立深度与地层压实综合恢复系数的拟合关系式,将关系式记作F;
步骤(9):利用步骤(8)得到的公式F与步骤(1)得到的时深关系F0,得到时间域的目的层顶面地震反射时间与地层压实综合恢复系数的关系式,将关系式记作F2
步骤(10):利用步骤(2)得到的目的层顶界面T和步骤(9)得到的关系式F2,获取连井剖面上沿着目的层顶面每个点处的地层压实综合恢复系数;
步骤(11):利用步骤(3)得到的时间域地层厚度Ht与步骤(10)得到的地层压实综合恢复系数,计算得到连井剖面上的时间域地层厚度压实恢复结果Tf
步骤(12):读取连井剖面上目的层段的地震道数据,按照步骤(11)得到的时间域地层厚度恢复结果Tf对剖面上每个地震道的长度进行拉伸,对拉伸后的目的层段地震道数据进行重新插值和采样;
步骤(13):对步骤(12)得到的处理后的目的层段地震数据进行剖面波形显示,在新的地震剖面上,利用传统的地震相模式,进行地震相识别与解释。
CN202111047129.4A 2021-09-07 2021-09-07 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程 Active CN113777657B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111047129.4A CN113777657B (zh) 2021-09-07 2021-09-07 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111047129.4A CN113777657B (zh) 2021-09-07 2021-09-07 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程

Publications (2)

Publication Number Publication Date
CN113777657A CN113777657A (zh) 2021-12-10
CN113777657B true CN113777657B (zh) 2023-10-17

Family

ID=78841773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111047129.4A Active CN113777657B (zh) 2021-09-07 2021-09-07 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程

Country Status (1)

Country Link
CN (1) CN113777657B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2170987A1 (en) * 1993-09-03 1995-03-09 Kenneth Rayvenor Lusty Armitage Evaluation of seismic sequence lithology
FR2893421A1 (fr) * 2005-11-14 2007-05-18 Inst Francais Du Petrole Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain.
CN105607145A (zh) * 2015-09-07 2016-05-25 中国石油天然气股份有限公司 碎屑岩地层的压实校正方法及装置
CN106094030A (zh) * 2016-08-24 2016-11-09 青岛海洋地质研究所 一种通过地震剖面定量恢复湖盆最大古水深的方法
CN106526679A (zh) * 2017-01-10 2017-03-22 中国石油大学(华东) 一种利用旋回分析法恢复不整合剥蚀厚度的方法
AU2020102025A4 (en) * 2020-07-02 2020-10-01 China University Of Geosciences, Beijing Method for Reconstructing Basin Paleogeomorphology
CN113189675A (zh) * 2021-05-13 2021-07-30 中国石油大学(华东) 一种砂岩压实作用数值模拟方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2170987A1 (en) * 1993-09-03 1995-03-09 Kenneth Rayvenor Lusty Armitage Evaluation of seismic sequence lithology
FR2893421A1 (fr) * 2005-11-14 2007-05-18 Inst Francais Du Petrole Methode d'evaluation quantitative des pressions de fluides et de detection des surpressions d'un milieu souterrain.
CN105607145A (zh) * 2015-09-07 2016-05-25 中国石油天然气股份有限公司 碎屑岩地层的压实校正方法及装置
CN106094030A (zh) * 2016-08-24 2016-11-09 青岛海洋地质研究所 一种通过地震剖面定量恢复湖盆最大古水深的方法
CN106526679A (zh) * 2017-01-10 2017-03-22 中国石油大学(华东) 一种利用旋回分析法恢复不整合剥蚀厚度的方法
AU2020102025A4 (en) * 2020-07-02 2020-10-01 China University Of Geosciences, Beijing Method for Reconstructing Basin Paleogeomorphology
CN113189675A (zh) * 2021-05-13 2021-07-30 中国石油大学(华东) 一种砂岩压实作用数值模拟方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Compaction trend versus seismic anisotropy in shaly formations;Pervukhina M 等;《Geophysical Prospecting》;第65卷(第5期);1351-1365 *
东营凹陷东南斜坡泥岩压实规律及超压分布;杨彬 等;《新疆石油地质》;第36卷(第6期);714-718+723 *
珠江口盆地番禺地区珠江组沉积前古地貌及其对沉积体系的控制;戴宗 等;《海相油气地质》(第03期);80-88 *

Also Published As

Publication number Publication date
CN113777657A (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
CN102053270B (zh) 一种基于沉积地层单元的地震相分析方法
US5226019A (en) Method of geophysical exploration
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
CN108802812A (zh) 一种井震融合的地层岩性反演方法
US9091786B2 (en) Image based effective medium modeling of the near surface earth formation
CN107121699A (zh) 一种地震相控制下的沉积微相识别方法
US20100118653A1 (en) Vertical seismic profiling velocity estimation method
CN107272062B (zh) 一种数据驱动的地下介质q场估计方法
US20130176822A1 (en) Method for estimating elastic parameters by inverting 4d seismic measurements
AU2007201922B2 (en) Method for sub-salt migration velocity analysis
CN111722284B (zh) 一种基于道集数据建立速度深度模型的方法
CN106526671B (zh) 一种碳酸盐岩储层地震波响应特征分析方法
CN112505754B (zh) 基于高精度层序格架模型的井震协同划分沉积微相的方法
CN113777657B (zh) 基于差异压实恢复的斜坡带连井剖面地震相表征工业化流程
CN109283577A (zh) 一种地震层位标定方法
CN109143335B (zh) 一种合成地震记录的制作方法、系统、介质及设备
CN108897039B (zh) 基于层位控制的dtw多波匹配方法
CN115598735A (zh) 一种基于速度重构的地震地貌恢复方法
CN107561581B (zh) 基于相关系数法建立多井模型的方法
CN113589365B (zh) 基于时频域信息的储层尖灭线描述方法
CN110673211B (zh) 一种基于测井与地震数据的品质因子建模方法
CN111060967A (zh) 一种井控速度场建模方法
CN107561580B (zh) 基于多层多井建立初始地质模型的方法
CN113589375A (zh) 基于倾斜层约束旅行时计算的vsp层速度反演方法
CN114706141B (zh) 一种基于沉积类比的地下曲流河带定量表征方法

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