CN109085647A - 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 - Google Patents
适用于断陷盆地陡坡带扇体顶面形变量的计算方法 Download PDFInfo
- Publication number
- CN109085647A CN109085647A CN201810855370.1A CN201810855370A CN109085647A CN 109085647 A CN109085647 A CN 109085647A CN 201810855370 A CN201810855370 A CN 201810855370A CN 109085647 A CN109085647 A CN 109085647A
- Authority
- CN
- China
- Prior art keywords
- fan body
- top surface
- deformation quantity
- body top
- fan
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 24
- 238000009825 accumulation Methods 0.000 claims abstract description 17
- 238000013139 quantization Methods 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 9
- 238000003325 tomography Methods 0.000 claims abstract description 5
- 238000004062 sedimentation Methods 0.000 claims abstract description 4
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 3
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 15
- 230000008569 process Effects 0.000 description 6
- 230000008021 deposition Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 240000007762 Ficus drupacea Species 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
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 (6)
1.适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,该适用于断陷盆地陡坡带扇体顶面形变量的计算方法包括:
步骤1,利用三维地震资料、实钻井资料及速度资料,计算扇体沉积之后上覆地层累积落差H;
步骤2,基于地质年代表计算扇体沉积之后累积的地质时间ΔT;
步骤3,通过垂直于断层走向的二维剖面确定断层真倾角α;
步骤4,综合落差、沉积时间及断层倾角,计算扇体沉积之后累积活动速率S;
步骤5,基于实钻井资料,利用扇体砂岩厚度与地层厚度的比值计算扇体砂岩百分含量P;
步骤6,利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量,并依据形变量的大小进行扇体顶面产状演化类型的划分。
2.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤1中,综合分析研究区三维地震资料、速度资料、井资料,计算扇体沉积之后上覆地层的累积落差H。
3.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤4中,利用落差H与断面倾角α,累积沉积时间ΔT计算出扇体沉积之后累积活动速率S,其中S=H/sin(α)*ΔT。
4.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤5中,基于实钻井资料计算扇体砂岩百分含量P,其中砂岩百分含量P等于扇体内部每一个砂体的厚度之和与扇体厚度的比值。
5.根据权利要求1所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤6中,利用累积断层活动速率S与扇体百分含量P的比值来量化表征扇体顶面形变量M,其中M=S/P。
6.根据权利要求5所述的适用于断陷盆地陡坡带扇体顶面形变量的计算方法,其特征在于,在步骤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 true CN109085647A (zh) | 2018-12-25 |
CN109085647B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376694A (zh) * | 2021-06-07 | 2021-09-10 | 中海石油(中国)有限公司深圳分公司 | 断陷盆地的活动性判断方法、装置、电子设备和存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140163890A1 (en) * | 2012-12-06 | 2014-06-12 | Roxar Software Solutions As | Method and system for presenting seismic information |
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 | 中国石油化工股份有限公司 | 基于地质统计学的浊积扇推进距离预测方法 |
-
2018
- 2018-07-30 CN CN201810855370.1A patent/CN109085647B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9354338B1 (en) * | 2012-02-22 | 2016-05-31 | Westerngeco L.L.C. | Generating synthetic seismic traces |
US20140163890A1 (en) * | 2012-12-06 | 2014-06-12 | Roxar Software Solutions As | Method and system for presenting seismic information |
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 |
---|
张伟忠 等: "东营凹陷边界断层活动性与沉积演化耦合关系的量化表征", 《中国石油大学学报》 * |
邱隆伟 等: "孤北洼陷长堤断裂带扇三角洲砂砾岩体沉积特征及空间展布规律", 《石油与天然气地质》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376694A (zh) * | 2021-06-07 | 2021-09-10 | 中海石油(中国)有限公司深圳分公司 | 断陷盆地的活动性判断方法、装置、电子设备和存储介质 |
CN113376694B (zh) * | 2021-06-07 | 2024-01-23 | 中海石油(中国)有限公司深圳分公司 | 断陷盆地的活动性判断方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109085647B (zh) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Macauley et al. | Slope channel sedimentary processes and stratigraphic stacking, Cretaceous Tres Pasos Formation slope system, Chilean Patagonia | |
CN105373648B (zh) | 一种砂泥岩互层型沉积体储层构型的建模方法 | |
CN108873066A (zh) | 弹性介质波动方程反射波旅行时反演方法 | |
CN105701319B (zh) | 一种水平井规则开发井网下的沉积微相建模方法 | |
CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN108508482A (zh) | 一种地下裂缝地震散射响应特征模拟方法 | |
CN102243678A (zh) | 一种基于沉积动力学反演的储集砂体分析方法 | |
CN105467464A (zh) | 近岸水下扇砂砾岩有效连通体划分和对比方法 | |
CN107870373A (zh) | 叠合盆地多期非连续性生烃潜力评价方法 | |
CN108376189B (zh) | 碎屑岩储层埋藏过程中成岩相演化的恢复方法 | |
CN107609254A (zh) | 细粒沉积物的沉积数值模拟方法 | |
CN105116449A (zh) | 一种弱反射储层的识别方法 | |
CN109085647A (zh) | 适用于断陷盆地陡坡带扇体顶面形变量的计算方法 | |
Tagliaferri et al. | Basin-scale, high-resolution three-dimensional facies modeling of tectonically confined turbidites: An example from the Firenzuola system (Marnoso-arenacea Formation, northern Apennines, Italy) | |
Yong et al. | Application of high frequency lake level change in the prediction of tight sandstone thin reservoir by sedimentary simulation | |
Xiaorong et al. | Dynamics of hydrocarbon accumulation in the west section of the northern margin of the Qaidam Basin, NW China | |
CN112329192B (zh) | 厚层砂岩油藏注聚均衡驱替设计方法 | |
CN100434934C (zh) | 重磁延拓回返垂直导数目标优化处理方法 | |
CN103675906B (zh) | 用于提高复杂碎屑岩储层预测精度的多控储层预测方法 | |
CN110688728B (zh) | 一种时间域内一维沉积物沉积特征及水体环境定量分析方法 | |
CN106842289B (zh) | 一种适用于测井约束反演的波阻抗曲线去压实处理方法 | |
CN105005078B (zh) | 一种基于速度模型的地震数据快速时深转换方法 | |
Molz et al. | A pseudo-fractal model for hydraulic property distributions in porous media | |
CN109931055A (zh) | 盆地深层复合成因的流体压力预测方法 | |
CN107290803A (zh) | 建立基于改进Pride-Lee模型的岩石物理新模型的方法 |
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 |