CN109085647B - 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 - Google Patents
适用于断陷盆地陡坡带扇体顶面形变量的计算方法 Download PDFInfo
- Publication number
- CN109085647B CN109085647B CN201810855370.1A CN201810855370A CN109085647B CN 109085647 B CN109085647 B CN 109085647B CN 201810855370 A CN201810855370 A CN 201810855370A CN 109085647 B CN109085647 B CN 109085647B
- Authority
- CN
- China
- Prior art keywords
- fan
- fan body
- top surface
- calculating
- deformation quantity
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000008021 deposition Effects 0.000 claims abstract description 18
- 230000000694 effects Effects 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 230000001186 cumulative effect Effects 0.000 claims abstract description 8
- 239000004576 sand Substances 0.000 claims description 6
- 230000035772 mutation Effects 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000005056 compaction Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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
本发明提供一种适用于断陷盆地陡坡带扇体顶面形变量的计算方法,包括:利用三维地震资料、实钻井资料及速度资料,计算扇体沉积之后上覆地层累积落差H;基于地质年代表计算扇体沉积之后累积的地质时间ΔT;通过垂直于断层走向的二维剖面确定断层真倾角α;综合落差、沉积时间及断层倾角,计算扇体沉积之后累积活动速率S;基于实钻井资料,利用扇体砂岩厚度与地层厚度的比值计算扇体砂岩百分含量P;利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量,并依据形变量的大小进行扇体顶面产状演化类型的划分。该方法全面的考虑了扇体形变量的主要影响因素,并进行了量化表征,方法合理,具有较广的推广应用价值。
Description
技术领域
本发明涉及油田勘探开发技术领域,特别是涉及到一种适用于断陷盆地陡坡带扇体顶面形变量的计算方法。
背景技术
陡坡带的砂砾岩扇体是一种重要的勘探目标。受到边界断层切割及上覆地层压实的影响,砂砾岩扇体顶面往往会产生构造回倾等形变现象。这种构造回倾对于油气聚集成藏较为有利,而由于扇体顶面往往覆盖巨厚的高速砾岩层,目前的三维地震资料在成像上会有一定的误差,使得砂砾岩扇体产生上拉现象,这使得在砂砾岩体顶面构造描述中产生较大的误差。目前这种误差难以通过三维地震的方法进行彻底消除,需要借助其他的一些方法进行综合分析。目前砂砾岩扇体顶面形变量的计算方面还处于空白状态,一方面是由于在早期的勘探研究中扇体形变量的计算未引起足够的重视,另一方面则是由于扇体在沉积之后收到多起构造运动及差异压实等作用的影响,其顶面形变量的主控因素较为复杂,难以用较为简单的方法进行量化表征。为此我们发明了一种新的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种应用简便,既考虑了外因断层活动的影响,也同时考虑了扇体本身内在的影响,较为全面的考虑了扇体形变量的主要影响因素的适用于断陷盆地陡坡带扇体顶面形变量的计算方法。
本发明的目的可通过如下技术措施来实现:适用于断陷盆地陡坡带扇体顶面形变量的计算方法,该适用于断陷盆地陡坡带扇体顶面形变量的计算方法包括:步骤1,利用三维地震资料、实钻井资料及速度资料,计算扇体沉积之后上覆地层累积落差H;步骤2,基于地质年代表计算扇体沉积之后累积的地质时间ΔT;步骤3,通过垂直于断层走向的二维剖面确定断层真倾角α;步骤4,综合落差、沉积时间及断层倾角,计算扇体沉积之后累积活动速率S;步骤5,基于实钻井资料,利用扇体砂岩厚度与地层厚度的比值计算扇体砂岩百分含量P;步骤6,利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量,并依据形变量的大小进行扇体顶面产状演化类型的划分。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,综合分析研究区三维地震资料、速度资料、井资料,计算扇体沉积之后上覆地层的累积落差H。
在步骤4中,利用落差H与断面倾角α,累积沉积时间ΔT计算出扇体沉积之后累积活动速率S,其中S=H/sin(α)*ΔT。
在步骤5中,基于实钻井资料计算扇体砂岩百分含量P,其中砂岩百分含量P等于扇体内部每一个砂体的厚度之和与扇体厚度的比值。
在步骤6中,利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量M,其中M=S/P。
在步骤6中,基于扇体顶面形变量M的大小将扇体顶面产状演化类型划分为突变型,渐变型和不变型3种类型。
本发明中的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,在理清扇体顶面形变主控因素的基础上,选择了断层活动速率及扇体砂岩百分含量这两个参数对扇体顶面形变量进行量化表征,在最大程度简化的基础上,充分考虑了扇体顶面形变的主要控制因素。因此本发明对于早期砂砾岩扇体顶面形变量的量化表征具有重要的参考意义,同时也能结合三维地震最大程度的表征扇体顶面的构造形态,为砂砾岩体勘探中圈闭描述的有效性提供有利的支撑。
附图说明
图1为本发明的适用于断陷盆地陡坡带扇体顶面形变量的计算方法的一具体实施例的流程图;
图2为本发明的一具体实施例中扇体顶面落差及断面倾角的计算方法示意图;
图3为本发明的一具体实施例中扇体砂岩百分含量计算方法示意图;
图4为本发明的一具体实施例中扇体顶面形变量的表征图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的适用于断陷盆地陡坡带扇体顶面形变量的计算方法的流程图。
在步骤101,利用三维地震资料计算扇体沉积之后上覆地层累积落差H。综合分析研究区三维地震资料、速度资料、井资料等,计算扇体沉积之后上覆地层的累积落差H,如图2,流程进入到步骤102;
在步骤102,基于地质年代表计算扇体沉积之后累积的地质时间ΔT。根据地层年代表确定目的层顶底界面沉积的时间间隔ΔT。流程进入到步骤103;
在步骤103,通过垂直于断层的二维剖面确定断层倾角α,如图2。流程进入到步骤104;
在步骤104,综合落差、沉积时间及断层倾角,计算扇体沉积之后累积活动速率S。利用落差H与断面倾角α,累积沉积时间ΔT计算出扇体沉积之后累积活动速率S,其中S=H/sin(α)*ΔT,如图2。流程进入到步骤105;
在步骤105,基于实钻井资料计算扇体砂岩百分含量P,如图3。其中砂岩百分含量P等于扇体内部每一个砂体的厚度之和与扇体厚度的比值。流程进入步骤106;
在步骤106,最终利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量,并依据形变量的大小进行扇体顶面产状演化类型的划分。利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量M,其中M=S/P。基于扇体顶面形变量M的大小可将扇体顶面产状演化类型划分为突变型,渐变型和不变型3种类型。图4为本发明的一具体实施例中扇体顶面形变量的表征图。依据扇体形变量M的大小将扇体顶面产状的演化过程划分为突变型,渐变型和不变形等3种类型。该方法适用于断陷盆地陡坡带扇体顶面形变量的估算。流程结束。
本发明中的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,综合考虑了扇体沉积之后的累积落差及扇体本身的砂岩百分含量,提出了针对陡坡带砂砾岩扇体顶面形变量的表征参数M。这种方法在最大程度简化的基础上,充分考虑了扇体顶面形变的主要控制因素。本发明应用简便,既考虑了外因断层活动的影响,也同时考虑了扇体本身内在的影响,较为全面的考虑了扇体形变量的主要影响因素,并进行了量化表征,方法合理。在断陷盆地陡坡带砂砾岩体勘探中具有较广的推广应用价值。
Claims (5)
1.适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,该适用于断陷盆地陡坡带扇体顶面形变量的计算方法包括:
步骤1,利用三维地震资料、实钻井资料及速度资料,计算扇体沉积之后上覆地层累积落差H;
步骤2,基于地质年代表计算扇体沉积之后累积沉积时间ΔT;
步骤3,通过垂直于断层走向的二维剖面确定断层真倾角α;
步骤4,综合累积落差H、累积沉积时间ΔT及断层真倾角α计算扇体沉积之后的累积断层活动速率S,其中S=H/(sin(α)*ΔT);
步骤5,基于实钻井资料,利用扇体砂岩厚度与地层厚度的比值计算扇体砂岩百分含量P;
步骤6,利用累积断层活动速率S与扇体砂岩百分含量P的比值来量化表征扇体顶面形变量,并依据形变量的大小进行扇体顶面产状演化类型的划分;
2.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤1中,综合分析研究区三维地震资料、速度资料、井资料,计算扇体沉积之后上覆地层累积落差H。
3.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤5中,基于实钻井资料计算扇体砂岩百分含量P,其中扇体砂岩百分含量P等于扇体内部每一个砂体的厚度之和与扇体厚度的比值。
4.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤6中,利用累积断层活动速率S与扇体砂岩百分含量P的比值来量化表征扇体顶面形变量M,其中M=S/P。
5.根据权利要求4所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤6中,基于扇体顶面形变量M的大小将扇体顶面产状演化类型划分为突变型,渐变型和不变型3种类型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810855370.1A CN109085647B (zh) | 2018-07-30 | 2018-07-30 | 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810855370.1A CN109085647B (zh) | 2018-07-30 | 2018-07-30 | 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109085647A CN109085647A (zh) | 2018-12-25 |
CN109085647B true CN109085647B (zh) | 2021-05-04 |
Family
ID=64830997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810855370.1A Active CN109085647B (zh) | 2018-07-30 | 2018-07-30 | 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109085647B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376694B (zh) * | 2021-06-07 | 2024-01-23 | 中海石油(中国)有限公司深圳分公司 | 断陷盆地的活动性判断方法、装置、电子设备和存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105388526A (zh) * | 2015-10-28 | 2016-03-09 | 中国石油化工股份有限公司 | 一种陆相湖盆的古地质结构恢复方法 |
US9354338B1 (en) * | 2012-02-22 | 2016-05-31 | Westerngeco L.L.C. | Generating synthetic seismic traces |
CN106199709A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | 扭张断裂活动性的评价方法 |
CN106772602A (zh) * | 2016-12-22 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 生长断层活动强度定量表征方法 |
CN107015290A (zh) * | 2017-03-13 | 2017-08-04 | 西北大学 | 一种改造型断陷盆地原始沉积面貌恢复的方法 |
CN107742016A (zh) * | 2017-09-29 | 2018-02-27 | 中国石油化工股份有限公司 | 一种适用于上升盘地层缺失型断层的活动性计算方法 |
CN108073744A (zh) * | 2016-11-16 | 2018-05-25 | 中国石油化工股份有限公司 | 基于地质统计学的浊积扇推进距离预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO20121471A1 (no) * | 2012-12-06 | 2014-06-09 | Roxar Software Solutions As | Fremgangsmåte og system for presentasjon av seismisk informasjon |
-
2018
- 2018-07-30 CN CN201810855370.1A patent/CN109085647B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9354338B1 (en) * | 2012-02-22 | 2016-05-31 | Westerngeco L.L.C. | Generating synthetic seismic traces |
CN105388526A (zh) * | 2015-10-28 | 2016-03-09 | 中国石油化工股份有限公司 | 一种陆相湖盆的古地质结构恢复方法 |
CN106199709A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | 扭张断裂活动性的评价方法 |
CN108073744A (zh) * | 2016-11-16 | 2018-05-25 | 中国石油化工股份有限公司 | 基于地质统计学的浊积扇推进距离预测方法 |
CN106772602A (zh) * | 2016-12-22 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 生长断层活动强度定量表征方法 |
CN107015290A (zh) * | 2017-03-13 | 2017-08-04 | 西北大学 | 一种改造型断陷盆地原始沉积面貌恢复的方法 |
CN107742016A (zh) * | 2017-09-29 | 2018-02-27 | 中国石油化工股份有限公司 | 一种适用于上升盘地层缺失型断层的活动性计算方法 |
Non-Patent Citations (2)
Title |
---|
东营凹陷边界断层活动性与沉积演化耦合关系的量化表征;张伟忠 等;《中国石油大学学报》;20170831;第41卷(第4期);第18-26页 * |
孤北洼陷长堤断裂带扇三角洲砂砾岩体沉积特征及空间展布规律;邱隆伟 等;《石油与天然气地质》;20180228;第39卷(第1期);第175-182页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109085647A (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vollgger et al. | Regional dome evolution and its control on ore-grade distribution: Insights from 3D implicit modelling of the Navachab gold deposit, Namibia | |
CN105701319B (zh) | 一种水平井规则开发井网下的沉积微相建模方法 | |
CN105388526B (zh) | 一种陆相湖盆的古地质结构恢复方法 | |
CN109870719B (zh) | 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统 | |
CN105467464B (zh) | 近岸水下扇砂砾岩有效连通体划分和对比方法 | |
EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
CN109884707B (zh) | 近地表分层时深曲线静校正方法 | |
US20160298427A1 (en) | Continuum sedimentary basin modeling using particle dynamics simulations | |
CN110632657B (zh) | 泥岩涂抹型断层封闭性分析方法及装置 | |
CN107742016A (zh) | 一种适用于上升盘地层缺失型断层的活动性计算方法 | |
CN109085647B (zh) | 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 | |
CN112147684A (zh) | 表征同沉积断层活动强度的方法和装置 | |
CN114910953B (zh) | 一种复杂断裂带不规则断面的解释方法 | |
CN111638552A (zh) | 一种古地貌恢复方法 | |
CN107450096B (zh) | 气藏储层的确定方法和装置 | |
McKinnon | Landslide runout: statistical analysis of physical characteristics and model parameters | |
CN113050161B (zh) | 一种基于断层交面线的断裂走向断距图解计算方法 | |
CN109164486A (zh) | 一种油气有利目标区的识别评价方法及系统 | |
CN115685327A (zh) | 一种深水浊积砂体展布范围的预测方法 | |
Murray et al. | A geostatistical approach to analyzing gold distribution controlled by large-scale fault systems–an example from Côte d’Ivoire | |
Liu et al. | Study on multi-period palaeotectonic stress fields simulation and fractures distribution prediction in Lannigou gold mine, Guizhou | |
CN110568149B (zh) | 沉积盆地烃源岩生排烃史精细快速定量模拟方法 | |
CN114076982B (zh) | 一种基于波形特征差异的岩溶古地貌恢复方法及装置 | |
CN112213465B (zh) | 一种基于电成像测井恢复古熔岩流方位及火山机构原始形态的方法 | |
CN110443890A (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 |