CN111610561A - 定量建立海平面变化及海进-海退曲线计算方法 - Google Patents

定量建立海平面变化及海进-海退曲线计算方法 Download PDF

Info

Publication number
CN111610561A
CN111610561A CN202010504873.1A CN202010504873A CN111610561A CN 111610561 A CN111610561 A CN 111610561A CN 202010504873 A CN202010504873 A CN 202010504873A CN 111610561 A CN111610561 A CN 111610561A
Authority
CN
China
Prior art keywords
sea
advance
deposition
sea level
retreat
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.)
Granted
Application number
CN202010504873.1A
Other languages
English (en)
Other versions
CN111610561B (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 Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
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 Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202010504873.1A priority Critical patent/CN111610561B/zh
Priority to AU2020101482A priority patent/AU2020101482A4/en
Publication of CN111610561A publication Critical patent/CN111610561A/zh
Application granted granted Critical
Publication of CN111610561B publication Critical patent/CN111610561B/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. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • 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/308Time lapse or 4D effects, e.g. production related effects to the formation
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Fluid Mechanics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了定量建立海平面变化及海进‑海退曲线计算方法,包括以下步骤:利用高分辨率的跨盆地的连续地震资料,建立目的层的高精度层序地层格架;通过对地震剖面和钻井、测井资料的综合分析,揭示区内沉积体系及沉积相的发育分布;选择目的层顶部或其上覆地层的一个古沉积水平面进行拉平;利用区内的时‑深数据,将时间剖面转换为深度剖面。本发明更精确地建立海平面变化曲线和海进‑海退曲线,从而探讨相对海平面变化和沉积物供给、构造沉降等对大陆边缘沉积层序结构和陆架边缘生长的控制作用;揭示大陆架的生长过程,有利于预测三角洲及斜坡扇砂体的分布,对于深水油气的勘探预测有重要的经济价值。

Description

定量建立海平面变化及海进-海退曲线计算方法
技术领域
本发明涉及海平面变化及海进-海退曲线计算方法技术领域,尤其涉及定量建立海平面变化及海进-海退曲线计算方法。
背景技术
近年来,有关滨岸/陆架边缘迁移轨迹的层序分析方法得到了广泛的关注。滨岸或陆架边缘轨迹(shoreline and shelf-edge trajectories)是指滨岸或陆架边缘坡折点在沉积倾向上的迁移轨迹,轨迹分析可客观地描述沉积体系随时间的迁移变化,并可反映出陆架边缘的生长过程。事实上,陆架边缘的生长就是在三角洲体系不断地进积-退积过程中实现的。
滨岸迁移轨迹具有垂向和侧向两个分量,垂向上的变化反映了海平面的上升、下降,而侧向上的变化反映了海进-海退的变化。因此,依据对滨岸迁移轨迹的垂向和侧向位移量的计算,可定量重建研究区海平面变化及海进-海退曲线。
目前,国内外研究海平面变化的方法主要有上超点法、地球化学法、古地磁法、古生物法等。依据滨岸迁移轨迹定量建立海平面曲线的研究与应用仍在起步阶段。
因此急需开发出定量建立海平面变化及海进-海退曲线计算方法,来解决上述问题。
发明内容
本发明要解决的技术问题是克服现有的缺陷,提供定量建立海平面变化及海进-海退曲线计算方法,更精确地建立海平面变化曲线和海进-海退曲线,从而探讨相对海平面变化和沉积物供给、构造沉降等对大陆边缘沉积层序结构和陆架边缘生长的控制作用;通过对滨岸迁移轨迹的建立,揭示陆架的生长过程,有利于预测三角洲及斜坡扇砂体的分布,对于深水油气的勘探预测有重要的经济价值,可以有效解决背景技术中的问题。
为了解决上述技术问题,本发明提供了如下的技术方案:
本发明提供定量建立海平面变化及海进-海退曲线计算方法,包括以下步骤:
S01:利用高分辨率的跨盆地的连续地震资料,建立目的层的高精度层序地层格架;
S02:通过对地震剖面和钻井、测井资料的综合分析,揭示区内沉积体系及沉积相的发育分布;
S03:选择目的层顶部或其上覆地层的一个古沉积水平面进行拉平;
S04:利用区内的时-深数据,将时间剖面转换为深度剖面;
S05:识别滨岸坡折点,在地震上通常可依据前积倾斜体的顶积层与前积层间的拐点确定;
S06:依据识别出的坡折点,绘制滨岸迁移轨迹;
S07:在等时地层单元内,分别计算滨岸迁移轨迹曲线的垂向位移量y和侧向位移量h;
S08:分别统计各坡折点的垂向深度yn和侧向坐标hn,则两个坡折点间的垂向位移量y=yn-yn-1,侧向位移量h=hn-hn-1,因此海平面上升时y记为正值,海平面下降时y为记负值,海进时h记为正值,海退时h记为负值;
S09:对沉积层厚度进行去压实校正,其中,校正后的垂向位移量v可通过下式计算(Allen P.A.and Allen J.R.,1992):
Figure BDA0002526163810000031
S10:确定三级层序的时限t后,假设沉积速率不变,以各沉积体的沉积厚度di估算得到沉积时间ti
ti=di/(d1+d2+d3+…+di+…+dn)×t (2)
计算各沉积体的沉积年龄Ai,可依据公式3进行计算:
Ai=T1+(t1+t2+…+ti) (3);
S11:分别建立海平面升降曲线和海进-海退曲线。
作为优选方案,步骤S10中,所述Ai为目的沉积体的年龄,T1为目的沉积体所在的三级层序底界面的年龄,t1、t2…ti分别为该层序内前i个沉积体各自的沉积时间。
作为优选方案,步骤S03中,所述古沉积水平面进行拉平包括滨岸平原或三角洲平原沉积及其延伸面,以恢复沉积期古斜坡,包括大陆斜坡或滨岸斜坡或三角洲斜坡的地貌形态。
作为优选方案,还包括横坐标分别为海平面垂直升降高度v(m)和海进/海退距离h(km),分别依据步骤S09和步骤S08得出;纵坐标为沉积年龄A(Ma),依据步骤S10得出。
作为优选方案,步骤S05还包括,或是通过对近滨和前滨沉积相的识别,将二者的分界作为滨岸坡折的位置。
作为优选方案,步骤S09中,
Figure BDA0002526163810000041
为表面孔隙度,c为压实系数,沉积层顶、底深度分别为yn和yn-1,去压实校正后沉积层的顶、底深度分别为vn和vn-1
本发明中提供的一个或多个技术方案,至少具有如下技术效果或者优点:
更精确地建立海平面变化曲线和海进-海退曲线,从而探讨相对海平面变化和沉积物供给、构造沉降等对大陆边缘沉积层序结构和陆架边缘生长的控制作用;通过对滨岸迁移轨迹的建立,揭示陆架的生长过程,有利于预测三角洲及斜坡扇砂体的分布,对于深水油气的勘探预测有重要的经济价值。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
在附图中:
图1是本发明实施例中定量建立海平面变化及海进-海退曲线计算方法中滨岸迁移轨迹垂向、侧向分量位移量的定量计算示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
为了更好地理解上述技术方案,下面将结合说明书附图以及具体实施方式对上述技术方案进行详细的说明。
实施例:
请参阅图1,本实施例提供定量建立海平面变化及海进-海退曲线计算方法,包括以下步骤:
步骤1:利用高分辨率的跨盆地的连续地震资料,建立目的层的高精度层序地层格架;
步骤2:通过对地震剖面和钻井、测井资料的综合分析,揭示区内沉积体系及沉积相的发育分布;
步骤3:选择目的层顶部或其上覆地层的一个古沉积水平面进行拉平,如滨岸平原或三角洲平原沉积及其延伸面等,以恢复沉积期古斜坡(大陆斜坡或滨岸斜坡或三角洲斜坡)的地貌形态;
步骤4:利用区内的时-深数据,将时间剖面转换为深度剖面;
步骤5:识别滨岸坡折点,在地震上通常可依据前积倾斜体的顶积层与前积层间的拐点确定,或是通过对近滨和前滨沉积相的识别,将二者的分界作为滨岸坡折的位置;
步骤6:依据识别出的坡折点,绘制滨岸迁移轨迹;
步骤7:在等时地层单元内,分别计算滨岸迁移轨迹曲线的垂向位移量y和侧向位移量h;
步骤8:分别统计各坡折点的垂向深度yn和侧向坐标hn,则两个坡折点间的垂向位移量y=yn-yn-1,侧向位移量h=hn-hn-1,因此海平面上升时y记为正值,海平面下降时y记为负值,海进时h记为正值,海退时h记为负值;
步骤9:对沉积层厚度进行去压实校正,其中,
Figure BDA0002526163810000051
为表面孔隙度,c为压实系数,主要与岩性有关;沉积层顶、底深度分别为yn和yn-1(在步骤8中已得到统计),去压实校正后沉积层的顶、底深度分别为vn和vn-1,则校正后的垂向位移量v可通过公式1计算(Allen P.A.andAllen J.R.,1992):
Figure BDA0002526163810000061
步骤10:确定三级层序的时限t后,假设沉积速率不变,以各沉积体的沉积厚度di估算得到沉积时间ti
ti=di/(d1+d2+d3+…+di+…+dn)×t (2)
计算各沉积体的沉积年龄Ai(Ma),可依据公式3进行计算:
Ai=T1+(t1+t2+…+ti) (3)
步骤11:Ai为目的沉积体的年龄,T1为目的沉积体所在的三级层序底界面的年龄,t1、t2…ti分别为该层序内前i个沉积体各自的沉积时间,依据步骤10计算得出;
步骤12:分别建立海平面升降曲线和海进-海退曲线,横坐标分别为海平面垂直升降高度v(m)和海进/海退距离h(km),分别依据步骤9和步骤8得出;纵坐标为沉积年龄A(Ma),依据步骤10得出。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.定量建立海平面变化及海进-海退曲线计算方法,其特征在于:包括以下步骤:
S01:利用高分辨率的跨盆地的连续地震资料,建立目的层的高精度层序地层格架;
S02:通过对地震剖面和钻井、测井资料的综合分析,揭示区内沉积体系及沉积相的发育分布;
S03:选择目的层顶部或其上覆地层的一个古沉积水平面进行拉平;
S04:利用区内的时-深数据,将时间剖面转换为深度剖面;
S05:识别滨岸坡折点,在地震上通常可依据前积倾斜体的顶积层与前积层间的拐点确定;
S06:依据识别出的坡折点,绘制滨岸迁移轨迹;
S07:在等时地层单元内,分别计算滨岸迁移轨迹曲线的垂向位移量y和侧向位移量h;
S08:分别统计各坡折点的垂向深度yn和侧向坐标hn,则两个坡折点间的垂向位移量y=yn-yn-1,侧向位移量h=hn-hn-1,因此海平面上升时y记为正值,海平面下降时y记为负值,海进时h记为正值,海退时h记为负值;
S09:对沉积层厚度进行去压实校正,其中,校正后的垂向位移量v可通过下式计算:
Figure FDA0002526163800000011
S10:确定三级层序的时限t后,假设沉积速率不变,以各沉积体的沉积厚度di估算得到沉积时间ti
ti=di/(d1+d2+d3+…+di+…+dn)×t (2)
计算各沉积体的沉积年龄Ai,可依据公式3进行计算:
Ai=T1+(t1+t2+…+ti) (3);
S11:分别建立海平面升降曲线和海进-海退曲线。
2.根据权利要求1所述的定量建立海平面变化及海进-海退曲线计算方法,其特征在于:步骤S10中,所述Ai为目的沉积体的年龄,T1为目的沉积体所在的三级层序底界面的年龄,t1、t2…ti分别为该层序内前i个沉积体各自的沉积时间。
3.根据权利要求2所述的定量建立海平面变化及海进-海退曲线计算方法,其特征在于:步骤S03中,所述古沉积水平面进行拉平包括滨岸平原或三角洲平原沉积及其延伸面,以恢复沉积期古斜坡,包括大陆斜坡或滨岸斜坡或三角洲斜坡的地貌形态。
4.根据权利要求1所述的定量建立海平面变化及海进-海退曲线计算方法,其特征在于:还包括横坐标分别为海平面垂直升降高度v(m)和海进/海退距离h(km),分别依据步骤S09和步骤S08得出;纵坐标为沉积年龄A(Ma),依据步骤S10得出。
5.根据权利要求1所述的定量建立海平面变化及海进-海退曲线计算方法,其特征在于:步骤S05还包括,或是通过对前滨和近滨沉积相的识别,将二者的分界作为滨岸坡折的位置。
6.根据权利要求1所述的定量建立海平面变化及海进-海退曲线计算方法,其特征在于:步骤S09中,φ0为表面孔隙度,c为压实系数,沉积层顶、底深度分别为yn和yn-1,去压实校正后沉积层的顶、底深度分别为vn和vn-1
CN202010504873.1A 2020-06-05 2020-06-05 定量建立海平面变化及海进-海退曲线计算方法 Active CN111610561B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010504873.1A CN111610561B (zh) 2020-06-05 2020-06-05 定量建立海平面变化及海进-海退曲线计算方法
AU2020101482A AU2020101482A4 (en) 2020-06-05 2020-07-24 Calculation method for quantitatively establishing sea level change and transgression-regression curves

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010504873.1A CN111610561B (zh) 2020-06-05 2020-06-05 定量建立海平面变化及海进-海退曲线计算方法

Publications (2)

Publication Number Publication Date
CN111610561A true CN111610561A (zh) 2020-09-01
CN111610561B CN111610561B (zh) 2021-05-07

Family

ID=72200312

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010504873.1A Active CN111610561B (zh) 2020-06-05 2020-06-05 定量建立海平面变化及海进-海退曲线计算方法

Country Status (2)

Country Link
CN (1) CN111610561B (zh)
AU (1) AU2020101482A4 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112861322A (zh) * 2021-01-15 2021-05-28 哈尔滨工程大学 一种海底阶梯式地貌演化定量分析方法及系统
CN113343161A (zh) * 2021-05-31 2021-09-03 西南石油大学 一种低位域陆架边缘三角洲的定量表征方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114594526B (zh) * 2020-12-03 2023-07-11 中国石油化工股份有限公司 河流相砂体井间平面连通性定量评价方法
CN114114461A (zh) * 2021-12-16 2022-03-01 中国石油大学(北京) 基于陆缘分类的海底扇富砂性风险评价方法及装置
CN114910953B (zh) * 2022-04-15 2023-03-17 中国石油化工股份有限公司 一种复杂断裂带不规则断面的解释方法
CN115373024B (zh) * 2022-08-09 2023-04-18 中国科学院南海海洋研究所 基于地层记录沉降反演被动陆缘地壳结构的方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005019868A1 (en) * 2003-08-23 2005-03-03 Westerngeco Seismic Holdings Ltd. Multiple attenuation method
US20050207844A1 (en) * 2002-06-28 2005-09-22 Paolo Boccotti Oscillating water column wave energy converter incorporated into caisson breakwater
CN106526689A (zh) * 2016-11-25 2017-03-22 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种定量恢复湖盆古水深的方法
CN107131858A (zh) * 2017-05-09 2017-09-05 中国地质大学(北京) 一种计算湖岸带碎浪坝沉积厚度的方法
CN109116421A (zh) * 2018-10-23 2019-01-01 中海石油(中国)有限公司 一种辫状河三角洲储层统计学反演参数确定方法
US10213145B1 (en) * 2015-10-01 2019-02-26 Cerner Innovation, Inc. Context-aware post traumatic stress disorder monitoring and intervention

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050207844A1 (en) * 2002-06-28 2005-09-22 Paolo Boccotti Oscillating water column wave energy converter incorporated into caisson breakwater
WO2005019868A1 (en) * 2003-08-23 2005-03-03 Westerngeco Seismic Holdings Ltd. Multiple attenuation method
US10213145B1 (en) * 2015-10-01 2019-02-26 Cerner Innovation, Inc. Context-aware post traumatic stress disorder monitoring and intervention
CN106526689A (zh) * 2016-11-25 2017-03-22 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种定量恢复湖盆古水深的方法
CN107131858A (zh) * 2017-05-09 2017-09-05 中国地质大学(北京) 一种计算湖岸带碎浪坝沉积厚度的方法
CN109116421A (zh) * 2018-10-23 2019-01-01 中海石油(中国)有限公司 一种辫状河三角洲储层统计学反演参数确定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
林畅松 等: "塔里木盆地奥陶纪碳酸盐岩台地的层序结构演化及其对盆地过程的响应", 《沉积学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112861322A (zh) * 2021-01-15 2021-05-28 哈尔滨工程大学 一种海底阶梯式地貌演化定量分析方法及系统
CN113343161A (zh) * 2021-05-31 2021-09-03 西南石油大学 一种低位域陆架边缘三角洲的定量表征方法

Also Published As

Publication number Publication date
AU2020101482A4 (en) 2020-09-03
CN111610561B (zh) 2021-05-07

Similar Documents

Publication Publication Date Title
CN111610561B (zh) 定量建立海平面变化及海进-海退曲线计算方法
CN111766630B (zh) 一种盆地古地貌的恢复方法
McCaffrey et al. Process controls on the development of stratigraphic trap potential on the margins of confined turbidite systems and aids to reservoir evaluation
NL1029338C2 (nl) Werkwijze voor het simuleren van de sedimentaire afzetting in een bekken, waarbij de dikten van de afzettingsreeksen wordt gerespecteerd.
CN104765065B (zh) 一种确定沉积地层古地貌的方法及装置
CN107991714B (zh) 基于湖盆古地形恢复的定量方法
Fernandez-Blanco et al. Cenozoic tectonics of the Tuz Gölü Basin (Central Anatolian Plateau, Turkey)
Eide et al. Sedimentology and reservoir properties of tabular and erosive offshore transition deposits in wave-dominated, shallow-marine strata: Book Cliffs, USA
CN115267909B (zh) 一种地表构造沉降计算方法及装置
SCHEIBER-ENSLIN et al. New depth maps of the Main Karoo Basin, used to explore the Cape isostatic anomaly, South Africa
CN106886047A (zh) 一种接收函数和重力联合反演地壳厚度和波速比的方法
Sarhan et al. Late Miocene normal faulting beneath the northern Nile Delta: NNW propagation of the Gulf of Suez Rift
CN101609167B (zh) 基于起伏地表的井间地震波动方程叠前深度偏移成像方法
CN113093286B (zh) 一种冷泉发育区储层非均质性的反演方法
CN108267797B (zh) 同生正断层发育区古地貌的恢复方法
CN107247289B (zh) 一种河流相地层的中期沉积旋回末期等时界面恢复方法
Onyenanu et al. Geometry, distribution and fill of erosional scours in a heterolithic, distal lower shoreface sandstone reservoir analogue: Grassy Member, Blackhawk Formation, Book Cliffs, Utah, USA
Rowan et al. Three-dimensional geometry and evolution of a composite, multilevel salt system, western Eugene Island, offshore Louisiana
CN112180463A (zh) 一种过渡层段砂岩展布的预测方法
CN116243380A (zh) 一种基于地层倾角的古地貌恢复方法
CN116738121B (zh) 基于古水深和剥蚀恢复计算盆地构造沉降的方法和装置
CN116859466B (zh) 井震时深批量标定方法、装置、电子设备及介质
CN116931089B (zh) 基于厚度分区的时深标定方法、装置、电子设备及介质
CN110967745A (zh) 一种火成岩的深度域速度建模方法
Biswal et al. Integrated Facies Modeling of fluvio-deltaic environment using seismic attributes and analogue training image

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