CN107831538A - 调节构造活动强度的定量描述方法 - Google Patents

调节构造活动强度的定量描述方法 Download PDF

Info

Publication number
CN107831538A
CN107831538A CN201710099351.6A CN201710099351A CN107831538A CN 107831538 A CN107831538 A CN 107831538A CN 201710099351 A CN201710099351 A CN 201710099351A CN 107831538 A CN107831538 A CN 107831538A
Authority
CN
China
Prior art keywords
transfer tectonics
activity
tectonics
transfer
activity intensity
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
CN201710099351.6A
Other languages
English (en)
Other versions
CN107831538B (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 Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
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 Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201710099351.6A priority Critical patent/CN107831538B/zh
Publication of CN107831538A publication Critical patent/CN107831538A/zh
Application granted granted Critical
Publication of CN107831538B publication Critical patent/CN107831538B/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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes

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)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种调节构造活动强度的定量描述方法,包括:步骤1,根据三维地震资料、测井资料、录井资料和分析化验资料确定与调节构造带相关的主断层或主构造的分布特征;步骤2,确定与调节构造带相关的主断层或主构造的活动强度;步骤3,根据三维地震资料、测井资料、录井资料和分析化验资料确定各关键时期调节构造带在平面上的几何特征;步骤4,结合步骤3得到的结果确定调节构造带参与形变的地质体体积;步骤5,根据步骤2和步骤4得到的结果,计算得出表征调节构造带活动强度的参数,即调解率。该方法反映了调节构造的活动性,表征其在不同地质时期,不同构造位置的活动强度特征,对指导类似构造位置的油气勘探具有深远的意义。

Description

调节构造活动强度的定量描述方法
技术领域
本发明涉及油田勘探开发技术领域,特别是涉及到一种调节构造带活动强度的定量化描述方法。
背景技术
调节构造带地质结构复杂,相关的研究较少,目前有横向伸展量和纵向伸展率这两个概念,是从调节构造的变形程度上对其进行定量化的研究,并未与相关的地质体活动程度相关联,不能清晰的描述调节构造带地质体的活动强度。在实际工作中除了需要表征调节构造的变形程度,还需要与构造活动强度相结合,为此发明了一种调节构造活动强度的定量化描述方法,即调节率,解决了以上技术问题。
发明内容
本发明的目的是提供一种调节构造活动强度的定量化描述方法,达到量化标准,建立适用于各类调节构造的活动强度定量化描述方法。
本发明的目的可通过如下技术措施来实现:调节构造活动强度的定量描述方法,该调节构造活动强度的定量描述方法包括:步骤1,根据三维地震资料、测井资料、录井资料和分析化验资料确定与调节构造带相关的主断层或主构造的分布特征;步骤2,确定与调节构造带相关的主断层或主构造的活动强度;步骤3,根据三维地震资料、测井资料、录井资料和分析化验资料确定各关键时期调节构造带在平面上的几何特征;步骤4,结合步骤3得到的结果确定调节构造带参与形变的地质体体积;步骤5,根据步骤2和步骤4得到的结果,计算得出表征调节构造带活动强度的参数,即调解率。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,根据三维地震资料、测井资料、录井资料和分析化验资料,对盆地断层的构造样式和组合特征进行分析,并建立盆地静态模型,针对调节构造带所涉及的各构造带和断层进行分析,得到各构造带和断层在各个关键时期的平面分布特征,同时重点对盆地静态模型进行分析,得到各构造带和断层的动态演化分析结果。
在步骤2中,通过综合分析,确定主断层与主构造的分布特征,明确构造与断层的形成机制;通过平衡剖面分析、盆地静态模型分析和构造应力场模拟,对表征关键时期内的断层活动性的参数进行求取,求取的参数包括表征垂向活动性的断层平均古落差,表征平面活动性的伸展率以及表征走滑盆地内走滑活动性的走滑量。
在步骤3中,进行调节构造带几何学特征分析,包括:对各反射层调节构造带进行分级、对调节构造涉及的断层进行倾向和走向的统计、对调节构造带的平面-剖面构造样式进行分析、对调节构造的类型进行划分、成图。
在步骤4中,得到调节构造带参与形变的地质体在平面与纵向上的分布特征,通过数学手段对该地质体进行合理的简化,从而得到易于计算的参与形变的地质体体积。
在步骤4中,将参与形变的地质体简化为横截面为椭圆形的柱状体,其长轴长度为断层端点之间的距离L,半短轴长度随应力背景而变化,在同一地质历史时期内可视为是常值a,圆周率值为π,因此其截面积S=π·a·L/2,该地质柱状体的高度为研究需要的地层单元厚度h,从而得到参与形变的地质体体积V,
V=S·h=h·π·a·L/2 (1)
其中,V—参与形变的地质体体积;S—椭圆截面的截面积;
h—研究需要的地层单元厚度h;π—圆周率;
a—椭圆截面的半短轴长度;L—椭圆截面的长轴长度。
在步骤5中,将主断层或主构造的活动性参数与地质体体积的比值作为调节率,来表征调节构造活动强度,主断层或主构造的活动性参数用R表示,则调节率Tr的计算公式可表示为:
将公式(2)中与应力背景相关的系数为简化为A,其中π指圆周率,a为椭圆截面的半短轴长度,则公式(2)表示为:
在步骤5中,计算得到调解率用来表征各种类型的调节构造活动强度,使用不同类型的活动性参数R将得到不同含义的调节率。
本发明中的调节构造活动强度的定量化描述方法,重点对研究区平面及剖面构造样式、断裂与构造的时空演化特征进行全面的分析,得出调节构造在各个关键时期的平面分布特征,以及相关的主构造和主断层的活动性参数,在一系列工作的基础上开展研究区内构造演化的动态分析,通过建立关键时期调节构造的平面分布图,最终结合主断层和主构造的活动性参数得出调节率来表征调节构造的活动强度。该方法有助于真实地反映调节构造的活动性,表征其在不同地质时期,不同构造位置的活动强度特征,从而为下一步分析调节构造的演化规律以及对油气成藏的控制作用,提供了研究基础。对指导类似构造位置的油气勘探具有深远的意义。
附图说明
图1为本发明的调节构造活动强度的定量描述方法的一具体实施例的流程图;
图2为本发明的一具体实施例中参与形变的地质体面积计算示意图;
图3为本发明的一具体实施例中参与形变的地质体体积计算示意图;
图4为本发明的一具体实施例中不同类型调节构造的各种调节率分类柱状图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的调节构造活动强度的定量描述方法的流程图。
在步骤101,根据三维地震资料、测井资料、录井资料和分析化验资料确定与该调节构造带相关的主断层或主构造的分布特征。具体说来,根据三维地震资料、测井资料、录井资料和分析化验资料,对盆地断层的构造样式和组合特征进行分析,并建立盆地静态模型,针对调节构造带所涉及的各构造带和断层进行分析,得到各构造带和断层在各个关键时期的平面分布特征,同时重点对盆地静态模型进行分析,得到各构造带和断层的动态演化分析结果。流程进入到步骤102。
在步骤102,运用多种方法确定与该调节构造带相关的主断层或主构造的活动强度。通过综合分析,确定主断层与主构造的分布特征,明确构造与断层的形成机制;通过平衡剖面分析、盆地静态模型分析和构造应力场模拟,对表征关键时期内的断层活动性的参数进行求取,求取的参数包括表征垂向活动性的断层平均古落差,表征平面活动性的伸展率以及表征走滑盆地内走滑活动性的走滑量,流程进入到步骤103。
在步骤103,根据三维地震资料、测井资料、录井资料和分析化验资料确定各关键时期调节构造带在平面上的几何特征。
通过综合分析,确定各关键时期调节构造带在平面上的几何特征。进行的调节构造带几何学特征分析包括:对各反射层调节构造带进行分级、对调节构造涉及的断层进行倾向和走向的统计、对调节构造带的平面-剖面构造样式进行分析、对调节构造的类型进行划分、成图。在一实施例中,针对惠民凹陷构造进行静态研究工作:各关键时期断裂分级、断层倾向、断裂走向统计、平面-剖面构造样式分析、断裂复杂性统计、基底断裂分析及成图。
在步骤104,结合步骤103得到的结果确定调节构造带参与形变的地质体体积。得到调节构造带参与形变的地质体在平面与纵向上的分布特征,通过数学手段对该地质体进行合理的简化,从而得到易于计算的参与形变的地质体体积。即得到调节构造带参与形变的地质体在平面与纵向上的分布特征,将参与形变的地质体简化为横截面为椭圆形的柱状体,其长轴长度为断层端点之间的距离L,半短轴长度随应力背景而变化,在同一地质历史时期内可视为是常值a,因此其表面积S=π·a·L/2,如图2。
随着地质体形变的过程进行,参与形变的地质体范围也会发生变化,两个断层最初没有叠置的范围,其椭圆形的截面积如图2-a所示,随着形变加剧,两条断层端点逐渐靠近,其椭圆形的截面积如图2-b所示,随着形变进一步加剧,断层之间开始叠置,其椭圆形的截面积如图2-c所示。该地质柱状体的高度为研究需要的地层单元厚度h,从而得到参与形变的地质体体积V=S·h=h·π·a·L/2,如图3。流程进入到步骤105。
在步骤105,根据步骤102和步骤104得到的结果,通过公式计算得出表征该调节构造带活动强度的参数,即调解率。将主断层或主构造的活动性参数与地质体体积的比值作为调节率,来表征调节构造活动强度。主断层或主构造的活动性参数用R表示,则调节率Tr的计算公式可表示为:
将其中与应力背景相关的简化为A,其中π指圆周率,a为椭圆截面的半短轴长度,则公式表示为:
该方法适用于各种类型的调节构造。使用不同类型的活动性参数R将得到表征不同含义的调节率,如图4。流程结束。
本发明中的调节构造活动强度的定量化描述方法主要是针对盆地中调节构造的,结合主断层和主构造的活动性参数得出调节率来表征调节构造的活动强度。该方法有助于真实地反映调节构造的活动性,表征各种类型的调节构造在不同地质时期,不同构造位置的活动强度特征,从而为下一步分析调节构造的演化规律以及对油气成藏的控制作用,提供了研究基础。对指导类似构造位置的油气勘探具有深远的意义。该项研究成果独具创新性,可填补国内外研究空白。

Claims (8)

1.调节构造活动强度的定量描述方法,其特征在于,该调节构造活动强度的定量描述方法包括:
步骤1,根据三维地震资料、测井资料、录井资料和分析化验资料确定与调节构造带相关的主断层或主构造的分布特征;
步骤2,确定与调节构造带相关的主断层或主构造的活动强度;
步骤3,根据三维地震资料、测井资料、录井资料和分析化验资料确定各关键时期调节构造带在平面上的几何特征;
步骤4,结合步骤3得到的结果确定调节构造带参与形变的地质体体积;
步骤5,根据步骤2和步骤4得到的结果,计算得出表征调节构造带活动强度的参数,即调解率。
2.根据权利要求1所述的调节构造活动强度的定量描述方法,其特征在于,在步骤1中,根据三维地震资料、测井资料、录井资料和分析化验资料,对盆地断层的构造样式和组合特征进行分析,并建立盆地静态模型,针对调节构造带所涉及的各构造带和断层进行分析,得到各构造带和断层在各个关键时期的平面分布特征,同时重点对盆地静态模型进行分析,得到各构造带和断层的动态演化分析结果。
3.根据权利要求1所述的调节构造活动强度的定量描述方法,其特征在于,在步骤2中,通过综合分析,确定主断层与主构造的分布特征,明确构造与断层的形成机制;通过平衡剖面分析、盆地静态模型分析和构造应力场模拟,对表征关键时期内的断层活动性的参数进行求取,求取的参数包括表征垂向活动性的断层平均古落差,表征平面活动性的伸展率以及表征走滑盆地内走滑活动性的走滑量。
4.根据权利要求1所述的调节构造活动强度的定量描述方法,其特征在于,在步骤3中,进行调节构造带几何学特征分析,包括:对各反射层调节构造带进行分级、对调节构造涉及的断层进行倾向和走向的统计、对调节构造带的平面-剖面构造样式进行分析、对调节构造的类型进行划分、成图。
5.根据权利要求1所述的调节构造活动强度的定量描述方法,其特征在于,在步骤4中,得到调节构造带参与形变的地质体在平面与纵向上的分布特征,通过数学手段对该地质体进行合理的简化,从而得到易于计算的参与形变的地质体体积。
6.根据权利要求5所述的调节构造活动强度的定量描述方法,其特征在于,在步骤4中,将参与形变的地质体简化为横截面为椭圆形的柱状体,其长轴长度为断层端点之间的距离L,半短轴长度随应力背景而变化,在同一地质历史时期内可视为是常值a,圆周率值为π,因此其截面积S=π·a·L/2,该地质柱状体的高度为研究需要的地层单元厚度h,从而得到参与形变的地质体体积V,
V=S·h=h·π·a·L/2 (1)
其中,V—参与形变的地质体体积;S—椭圆截面的截面积;
h—研究需要的地层单元厚度h;π—圆周率;
a—椭圆截面的半短轴长度;L—椭圆截面的长轴长度。
7.根据权利要求6所述的调节构造活动强度的定量描述方法,其特征在于,在步骤5中,将主断层或主构造的活动性参数与地质体体积的比值作为调节率,来表征调节构造活动强度,主断层或主构造的活动性参数用R表示,则调节率Tr的计算公式可表示为:
<mrow> <mi>T</mi> <mi>r</mi> <mo>=</mo> <mfrac> <mi>R</mi> <mrow> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <mi>a</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mi>L</mi> <mn>2</mn> </mfrac> <mo>&amp;CenterDot;</mo> <mi>h</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
将公式(2)中与应力背景相关的系数为简化为A,其中π指圆周率,a为椭圆截面的半短轴长度,则公式(2)表示为:
<mrow> <mi>T</mi> <mi>r</mi> <mo>=</mo> <mi>A</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mi>R</mi> <mrow> <mi>L</mi> <mo>&amp;CenterDot;</mo> <mi>h</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
8.根据权利要求7所述的调节构造活动强度的定量描述方法,其特征在于,在步骤5中,计算得到调解率用来表征各种类型的调节构造活动强度,使用不同类型的活动性参数R将得到不同含义的调节率。
CN201710099351.6A 2017-02-23 2017-02-23 调节构造活动强度的定量描述方法 Active CN107831538B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710099351.6A CN107831538B (zh) 2017-02-23 2017-02-23 调节构造活动强度的定量描述方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710099351.6A CN107831538B (zh) 2017-02-23 2017-02-23 调节构造活动强度的定量描述方法

Publications (2)

Publication Number Publication Date
CN107831538A true CN107831538A (zh) 2018-03-23
CN107831538B CN107831538B (zh) 2019-07-12

Family

ID=61643088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710099351.6A Active CN107831538B (zh) 2017-02-23 2017-02-23 调节构造活动强度的定量描述方法

Country Status (1)

Country Link
CN (1) CN107831538B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110850468A (zh) * 2019-10-09 2020-02-28 中国石油大学(华东) 一种定量表征弯曲走滑断层增压、释压强度的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130238297A1 (en) * 2012-03-09 2013-09-12 Schlumberger Technology Corporation Faulted geological structures containing unconformities
CN106324701A (zh) * 2016-10-10 2017-01-11 中国海洋石油总公司 基于断‑隆联动含油气盆地分析方法
CN106443783A (zh) * 2016-10-31 2017-02-22 中国地质大学(北京) 一种基于断层活动性的多期次裂缝定量预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130238297A1 (en) * 2012-03-09 2013-09-12 Schlumberger Technology Corporation Faulted geological structures containing unconformities
CN106324701A (zh) * 2016-10-10 2017-01-11 中国海洋石油总公司 基于断‑隆联动含油气盆地分析方法
CN106443783A (zh) * 2016-10-31 2017-02-22 中国地质大学(北京) 一种基于断层活动性的多期次裂缝定量预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邹娟 等: ""构造活动强度划分断陷盆地构造区划—以饶阳凹陷为例"", 《石油学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110850468A (zh) * 2019-10-09 2020-02-28 中国石油大学(华东) 一种定量表征弯曲走滑断层增压、释压强度的方法

Also Published As

Publication number Publication date
CN107831538B (zh) 2019-07-12

Similar Documents

Publication Publication Date Title
CN107290782B (zh) 储层孔隙度、含水饱和度和泥质含量参数同时反演新方法
CN108072902B (zh) 一种碳酸盐岩岩石物理模型的建模方法
CN102426390B (zh) 一种非均质泥砂岩储层储量确定方法
CN109447532A (zh) 一种基于数据驱动的油藏井间连通性确定方法
CN109061728B (zh) 一种滩坝砂体精细预测方法
CN106556866B (zh) 一种利用岩石物理模板对储层定量解释的方法
CN104239666A (zh) 一种基于层次分析法的煤层气综合评价方法
CN104750884A (zh) 基于多因素非线性回归的页岩油气富集指数定量评价方法
Shi et al. Brittleness index prediction in shale gas reservoirs based on efficient network models
CN103728659A (zh) 一种提高地下岩溶探测精度的方法
CN109184677A (zh) 用于非均质互层砂体的储层评价方法
CN109490965A (zh) 一种定量评价地层非均匀性的方法及装置
CN105445441A (zh) 储层含气饱和度确定方法及装置
CN104122581B (zh) 一种叠后声波阻抗反演方法
CN105652316A (zh) 一种基于裂缝模型的智能优化地震多属性融合方法
CN102455438A (zh) 碳酸盐岩缝洞型储层体积预测方法
CN104570065B (zh) 一种利用地震波阻抗定量反演孔隙度的方法
CN107780923A (zh) 一种基于泥质校正的含水饱和度模型的建立、仿真方法
CN106443800A (zh) 一种中基性火山岩的核磁孔隙度校正方法
CN103376468A (zh) 基于神经网络函数逼近算法的储层参数定量表征方法
CN107728205A (zh) 一种地层压力预测方法
CN107831538B (zh) 调节构造活动强度的定量描述方法
CN115877447A (zh) 一种直平联合井网条件下地震约束三维地质建模的储层预测方法
CN104991280A (zh) 一种基于微观孔喉结构指数的射孔层段选择方法及装置
CN106908838A (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