CN114185083A - 一种碎屑岩地层中断层封闭性的定量评价方法 - Google Patents
一种碎屑岩地层中断层封闭性的定量评价方法 Download PDFInfo
- Publication number
- CN114185083A CN114185083A CN202111484074.3A CN202111484074A CN114185083A CN 114185083 A CN114185083 A CN 114185083A CN 202111484074 A CN202111484074 A CN 202111484074A CN 114185083 A CN114185083 A CN 114185083A
- Authority
- CN
- China
- Prior art keywords
- fault
- section
- stratum
- stress
- shale content
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000011435 rock Substances 0.000 title claims abstract description 45
- 238000007789 sealing Methods 0.000 title claims abstract description 36
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 36
- 239000012530 fluid Substances 0.000 claims abstract description 24
- 238000011161 development Methods 0.000 claims abstract description 13
- 230000000903 blocking effect Effects 0.000 claims abstract description 10
- 238000011156 evaluation Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 34
- 239000004215 Carbon black (E152) Substances 0.000 claims description 12
- 229930195733 hydrocarbon Natural products 0.000 claims description 12
- 150000002430 hydrocarbons Chemical class 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 10
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 6
- 238000012512 characterization method Methods 0.000 claims description 4
- 238000005056 compaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000007480 spreading Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000009933 burial Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 238000002224 dissection Methods 0.000 claims 1
- 239000010802 sludge Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 3
- 238000013210 evaluation model Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 238000011835 investigation Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 208000006440 Open Bite Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000002920 hazardous waste Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 210000001503 joint Anatomy 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000009919 sequestration Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000002699 waste material Substances 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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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
本发明提供一种碎屑岩地层中断层封闭性的定量评价方法，包括以下步骤：(1)构建断面、层位模型；(2)地层泥质含量模型；(3)地应力场与地层流体压力模型；(4)断面参数刻画；(5)评价参数提取；(6)断面封闭模型建立。本发明技术方案在操作过程中通过利用断面泥质含量、滑动趋势、膨胀趋势三个参数建立综合评价模型，综合考虑地层泥质含量、断层空间发育特征(走向、倾向、断距)、地应力以及地层流体压力等因素对断层封闭性的控制作用。对于全面、正确理解断层封闭主控因素、预测潜在断层封闭性具有重要帮助。
Description
技术领域
本发明属于油气资源地质勘探评价领域及生态环境领域，具体涉及一种碎屑岩地层中断 层封闭性的定量评价方法。
背景技术
断层是地应力作用下形成的分割地层的构造。在不同地质历史时期及不同地质条件下， 断层可具有封闭地层流体与输导地层流体两种不同的状态。具有较好封闭性的断层对于断层 油气藏的形成至关重要。此外，废弃油气田中具有强封闭能力的断层圈闭也是危险废弃物处 理以及二氧化碳捕获与封存的重要靶区。正确评价断层封闭性除具有重要经济价值外也具有 显著的社会价值。
随着人们对断层认识的加深，影响断层封闭能力的因素及评价方法也不断变化。早期认 为断层是一个剪切面，低渗透层与储层对接是导致断层封闭的主要原因，利用三角图解模拟 断面两侧地层接触特征是断层封闭性的主要研究、评价方法。随着研究的深入，人们意识到 断层不仅是一个简单平面，而是一个包含断层核和破碎带的二维结构，断面空间形态(倾向、 倾角、埋深)、规模、地层岩性、地层流体压力、地应力均对断层封闭性造成影响。评价对象 逐渐变为断裂带岩性、断面应力等要素。诸如SSF、SGR、断面正应力、断-储排替压力差等 参数逐渐被用来评价断层的封闭性。
然而，目前主要的评价模型大多将重点放在断层岩岩性及断面正应力特征(断面紧闭性)， 忽略沿断层方向的稳定性对封闭性的影响。普遍认为处于沿断面滑动状态的断层对流体起输 导作用(不具有封闭性)，趋于滑动(即不稳定)的断层封闭封闭性也弱。断层的稳定性缺乏 有效表征。另外，在多个参数均影响断层封闭性的背景下，明确各参数对断层封闭性的影响， 建立定量评价指标，对于准确认识断层封闭性具有重要意义。
发明内容
基于以上分析可以看出，现阶段缺乏基于断面岩性、稳定性及断面紧性综合评价断层封 闭性的方法，本发明公开一种碎屑岩地层中断层封闭性的定量评价方法。
具体实现技术方案如下：
(1)构建断面、层位模型：
主要基于三维地震数据体，对目标断层及地层进行精细解释，刻画断层面在空间不同位 置处的走向、倾向、倾角及断层断距特征；
上述步骤(1)中，利用地震解释软件，对三维地震数据体以3～5道为间隔进行解释，实 现对研究区断层、地层的高精度刻画，从而重点获得地层的空间分布，断层的走向、倾向、 倾角、断距。其中，地震数据解释间隔可根据三维数据体精度及断层发育规模进行调整，地 震采样间距大时减小解释间隔，断层规模小时也需减小地震数据解释间隔，反之，可增大解 释间隔。总体目标为以较高效率对目标断层、地层的展布特征进行刻画，构建三维地质模型。
(2)地层泥质含量模型：
利用测井中反映泥质含量的GR曲线，对地层泥质含量进行计算，获得单井的泥质分布 特征，通过插值的方法，对三维空间的泥质含量特征进行刻画；
上述步骤(2)中，首先利用测井曲线对单井泥质含量的分布进行计算，泥质含量计算主 要用公式：
Vsh＝0.083×(23.7×Igr-1)(第三系地层) (公式1)
Vsh＝0.33×(22.1×Igr-1)(前第三系地层) (公式2)
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数。
针对不同时代地层，选择公式1或公式2进行计算。计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到。在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
(3)地应力场与地层流体压力模型：
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向 及变化梯度、垂向应力的变化梯度与地层流体压力。
上述步骤(3)中，水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征， 不同井之间显示出的水平最大主应力方向可能存在一定差异，通过取平均值的方法获得研究 区水平最大主应力方向。水平最小主应力与水平最大主应力角度差为90°。
由于密度测井资料通常十分全面，几乎包含了全部井段。因此，垂向主应力大小通过对 岩石密度的积分获得，计算公式为：
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单 位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m。
水平最大与最小主应力主要依靠偶极子声波测井值计算获得，计算公式为：
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分 别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，单位％；φc为临界孔隙度(一般取40)。
由于成本导致偶极子声波测井通常针对特定层位，研究层位可能无相关数据。对此，在 有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深 度水平最大、最小主应力值。(除此方法外，在样品充裕的情况下，也可通过岩石破裂法、声 发射法、油田现场井壁崩落资料、储层压裂改造资料获得水平最大与最小主应力数据。)
地层流体压力计算主要针对地层中存在的超压，利用等效法进行计算，公式为：
Pz＝γwZe+γr(Z-Ze)＝ρrgZ-(ρτ-ρw)gZe (公式8)
式中，PZ为欠压实泥岩的地层流体压力，单位MPa；γw为等效深度的静水压力梯度，单 位MPa/m；Ze为欠压实泥岩所对应的等效深度，单位m；γr为埋藏深度与等效深度间的静岩 压力梯度，单位MPa/m；Z为欠压实泥岩的埋深，单位m；ρr为埋藏深度与等效深度间的岩层平均密度，单位kg/m3；g为重力加速度，取9.8N/kg；ρw为地层水密度，取1000kg/m3； Δt为欠压实泥岩的声波时差，单位μs/m；Δt0为地表的原始声波时差，单位μs/m；C为正 常压实泥岩的压实系数，单位m-1。
通过以上计算，获得水平最大、最小主应力大小及方向，垂向应力变化梯度以及地层流 体压力特征。
(4)断面参数刻画
上述步骤(4)中需要计算的参数包括断面泥质含量(SGR)、滑动趋势(ST)和膨胀趋势(DT)。其中，断面泥质含量用来表征断层的岩性特征，滑动趋势和膨胀趋势用来表征断面的应力特征。断面泥质含量为沿断层滑过计算点的累计泥质含量，泥质含量越高，则认为断层岩泥质含量越高，封闭性越强，反之越弱。滑动趋势用来描述断层两盘沿断面滑动可能性的大小，滑动趋势越大，断层越趋于滑动，趋于不稳定，封闭能力越弱，反之，断层趋于 稳定，封闭能力强。膨胀趋势用来表征断面的紧闭程度，膨胀趋势越大，断层趋于张开，封 闭性越弱，反之，断层越紧闭，趋于封闭。
断面泥质含量计算公式为：
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑 过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断 距，单位m。Vshi为步骤(2)中空间插值得到的紧邻断层地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得。
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-PZ) (公式12)
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
式中，ST为断面滑动趋势，小数；DT为膨胀趋势，小数；τ为作用于断面上的剪应力，单位MPa；δn为作用在断面上的正应力，单位MPa；PZ为孔隙流体压力，单位MPa，利用公 式8计算。其中，δH为水平最大主应力，MPa；δh为水平最小主应力，MPa；σv为垂向主应 力，MPa；αn为断面法线与σH的夹角，单位°；βn为断面法线与σh的夹角，单位°；γn为断面 法线与σv的夹角，单位°。
(5)评价参数提取
根据勘探实际情况，对典型断层型油气藏及不含油气的断层圈闭进行精细解剖，确定封 闭断层与不封闭断层发育范围，提取断面参数信息。其中，对于典型油气藏，基于测井解释、 试油分析明确油气分布层位，通过地震属性刻画油气空间范围，特别是油气藏与断层相接触 的范围。对比分析断面封闭油气藏范围及储层范围，确定封闭断面范围及不封闭断面范围， 进而提取封闭断层与不封闭断层断面参数。对于不含油气的断层圈闭，首先通过成藏条件分 析，明确储层发育特征、下部烃源岩的供烃能力以及顶部盖层封闭能力。选取发育储层、下 部烃源岩具有供烃能力、顶部盖层具有较强封闭能力，油气未成藏主要由于断层封闭性弱， 提取断面与储层相交范围内的断面参数作为不封闭断层参数特征。在参数提取过程中，以断 面和储层相交区域为参数提取范围，根据实际情况(工区大小及地震道间隔)每隔数个地震 道，提取断层的断面泥质含量、滑动趋势、膨胀趋势参数。
(6)断面封闭模型建立
在明确封闭断层和不封闭断层各参数的基础上，通过费希尔(Fisher)判别分析，建立基 于三参数的概率判别函数。函数建立过程中，随机抽取80％数据用于构建判别方程，剩余20％ 数据用于验证判别方程可靠性。
本发明技术方案在操作过程中通过利用断面泥质含量、滑动趋势、膨胀趋势三个参数建 立综合评价模型，综合考虑地层泥质含量、断层空间发育特征(走向、倾向、倾角、断距)、 地应力以及地层流体压力对断层封闭性的控制作用。对于全面、正确理解断层封闭主控因素、 预测潜在断层封闭性具有重要帮助。
附图说明
图1本发明方法的步骤框图；
图2a为本发明方法步骤的水平最大主应力(σH)拟合示意图；
图2b为本发明方法步骤的水平最小主应力(σh)拟合示意图；
图3为本发明方法步骤的断面上某点泥质含量计算示意图；
图4为本发明方法步骤的σn和τn计算时的各参数空间示意图；
图5a为本发明方法步骤的研究区(胜利油田)某断层的断面泥质含量三维展布；
图5b为本发明方法步骤的研究区(胜利油田)某断层的断面滑动趋势三维展布；
图5c为本发明方法步骤的研究区(胜利油田)某断层的断面膨胀趋势三维展布；
图6a为本发明方法步骤的研究区某断层附近含气储层通过测井解释示意图；
图6b为本发明方法步骤的研究区某断层附近过含气储层的地震剖面图；
图6c为本发明方法步骤的研究区某断层附近过含气储层的均方根振幅属性图(与6b为 同一剖面，提取均方根振幅属性后得到)；
图6d为本发明方法步骤的研究区某断层附近含气储层通过均方根振幅属性数据及体雕 刻刻画气藏三维空间分布特征示意图；
图7a为本发明方法步骤的研究区某断层上封闭断层范围示意图；
图7b为本发明方法步骤的研究区某断层下盘储层在断层上的投影示意图。
具体实施方式
下面结合具体实施例对本发明作进一步描述，在此发明的示意性实施例以及说明用来解 释本发明，但并不作为对本发明的限定。
下面结合附图与实例对本发明作进一步的详细说明；可以理解的是，此处所叙述的具体 实例仅仅用于解释本发明，而非对本发明的限定；另外还需要说明的是，为了便于描述，附 图中仅示出了与本发明相关的部分而非全部内容。
实施例选取济阳坳陷三合村油田内的断层进行说明。
如图1所示，碎屑岩地层中断层封闭性的定量评价方法，所述方法包括以下步骤：
(1)构建断面、层位模型：
主要基于三维地震数据体，对目标断层及地层进行精细解释，刻画断层在空间不同位置 处的走向、倾向、倾角及断层断距特征；
上述步骤(1)中，利用地震解释软件，对三维地震数据体以3～5道为间隔进行解释，实 现对研究区断层、地层的高精度刻画，从而重点获得地层的空间分布，断层的走向、倾向、 倾角、断距。其中，地震数据解释间隔可根据三维数据体精度及断层发育规模进行调整，地 震采样间距大时减小解释间隔，断层规模小时也需减小地震数据解释间隔，反之，可增大解 释间隔。总体目标为以较高效率完成对目标断层、地层的展布特征进行刻画，构建三维地质 模型。
(2)地层泥质含量模型：
利用测井中反映泥质含量的GR曲线，对地层泥质含量进行计算，获得单井的泥质分布 特征，通过插值的方法，对三维空间的泥质含量特征进行刻画，为后续断面泥质含量刻画做 准备；
上述步骤(2)中，首先利用测井曲线对单井泥质含量的分布进行计算，泥质含量计算主 要用公式：
Vsh＝0.083×(23.7×Igr-1)(第三系地层) (公式1)
Vsh＝0.33×(22.1×Igr-1)(前第三系地层) (公式2)
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数。
针对不同时代地层，选择公式1或公式2进行计算。计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到。在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
(3)地应力场与地层流体压力模型：
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向 及变化梯度、垂向应力的变化梯度。
上述步骤(3)中，水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征， 不同井之间显示出的水平最大主应力方向可能存在一定差异，通过取平均值的方法获得研究 区水平最大主应力方向。水平最小主应力与水平最大主应力角度差为90°。
由于密度测井资料通常十分全面，几乎包含了全部井段。因此，垂向主应力大小通过对 岩石密度的积分获得，计算公式为：
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单 位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m。
水平最大与最小主应力主要依靠偶极子声波测井值计算获得，计算公式为：
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分 别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，％；φc为临界孔隙度(一般取40％)。
由于成本导致偶极子声波测井通常针对特定层位，研究层位可能无相关数据。对此，在 有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深 度水平最大、最小主应力值，如图2a和图2b。除此方法外，在样品充裕的情况下，也可通 过岩石破裂法、声发射法、油田现场井壁崩落资料、储层压裂改造资料获得水平最大与最小 主应力数据。
地层流体压力计算主要针对地层中存在的超压，利用等效法进行计算，公式为：
Pz＝γwZe+γr(Z-Ze)＝ρrgZ-(ρr-ρw)gZe (公式8)
式中，PZ为欠压实泥岩的地层流体压力，单位MPa；γw为等效深度的静水压力梯度，单 位MPa/m；Ze为欠压实泥岩所对应的等效深度，单位m；γr为埋藏深度与等效深度间的静岩 压力梯度，单位MPa/m；Z为欠压实泥岩的埋深，单位m；ρr为埋藏深度与等效深度间的岩层平均密度，单位kg/m3；g为重力加速度，取9.8N/kg；ρw为地层水密度，取1000kg/m3； Δt为欠压实泥岩的声波时差，单位μs/m；Δt0为地表的原始声波时差，单位μs/m；C为正 常压实泥岩的压实系数，单位m-1。
通过以上计算，获得水平最大、最小主应力大小及方向，垂向应力变化梯度以及地层流 体压力特征。
(4)断面参数刻画
上述步骤(4)中需要计算的参数包括断面泥质含量(SGR)、滑动趋势(ST)和膨胀趋势(DT)。其中，断面泥质含量用来表征断层的岩性特征，滑动趋势和膨胀趋势用来表征断面的应力特征。断面泥质含量为沿断层滑过计算点的累计泥质含量，泥质含量越高，则认为断层岩泥质含量越高，封闭性越强，反之越弱。滑动趋势用来描述断层两盘沿断面滑动可能性的大小，滑动趋势越大，断层越趋于滑动，趋于不稳定，封闭能力越弱，反之，断层趋于 稳定，封闭能力强。膨胀趋势用来表征断面的紧闭程度，膨胀趋势越大，断层趋于张开，封 闭性越弱，反之，断层越紧闭，趋于封闭。
断面泥质含量原理如图3，计算公式为：
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑 过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断 距，单位m。Vshi为步骤(2)中空间插值得到的紧邻断层地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得。
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-PZ) (公式12)
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
式中，ST为断面滑动趋势，小数；DT为膨胀趋势，小数；τ为作用于断面上的剪应力，单位MPa；δn为作用在断面上的正应力，单位MPa；PZ为孔隙流体压力，单位MPa，利用公 式8计算。其中，δH为水平最大主应力，MPa；δh为水平最小主应力，MPa；σv为垂向主应 力，MPa；αn为断面法线与σH的夹角，单位°；βn为断面法线与σh的夹角，单位°；γn为断面 法线与σv的夹角，单位°。各参数空间意义见图4。以渤海湾盆地济阳坳陷三合村洼陷内的垦 116断层为例对各参数进行刻画，图5a到图5c为对该断面泥质含量、滑动趋势及膨胀趋势的 空间刻画结果。
(5)评价参数提取
根据勘探实际情况，对典型断层型油气藏及不含油气的断层圈闭进行精细解剖，确定封 闭断层与不封闭断层发育范围，提取断面参数信息。其中，对于典型油气藏，基于测井解释、 试油分析明确油气分布层位，通过地震属性刻画油气空间范围，特别是油气藏与断层相接触 的范围。对比分析断面封闭油气藏范围及储层范围，确定封闭断面范围及不封闭断面范围， 进而提取封闭断层与不封闭断层断面参数。对于不含油气的断层圈闭，首先通过成藏条件分 析，明确储层发育特征、下部烃源岩的供烃能力以及顶部盖层封闭能力。选取发育储层、下 部烃源岩具有供烃能力、顶部盖层具有较强封闭能力，油气未成藏主要由于断层封闭性弱的 断层，提取断面与储层相交范围内的断面参数作为不封闭断层参数特征。在参数提取过程中， 以断面和储层相交区域为范围，根据实际情况(工区大小及地震道间隔)每隔数个地震道， 提取断层的断面泥质含量、滑动趋势、膨胀趋势参数。
以渤海湾盆地济阳坳陷内的垦116断层及其附近气藏为例对提取方法进行说明。通过测 井解释与试油资料确定含气储层深度(图6a)。由地震剖面及地震均方根振幅属性(RMS) 剖面可知，气藏上倾方向被断层遮挡(图6b)。由地震RMS属性强度在断层外明显降低(不 含气)可知断层起封闭作用(图6c)。根据含气储层具有高RMS属性的特征应用体雕刻技术 对气藏形态进行刻画(图6d)获得气藏的空间展布特征。将气藏空间分布范围与断层空间分 布范围叠加可以明显看出具有封闭能力的断层的范围(图7a)。由断层下盘地层在断面投影 可以看出，气藏所在的砂体储层空间上连续稳定(图7b虚线范围内地层泥质含量均小于0.4)， 由此可确定断层上具有封闭能力和不具有封闭能力断层的范围。结合断面上泥质含量、滑动 趋势及膨胀趋势的三维分布特征(图5a到图5c)，可获得具有封闭能力和不具有封闭能力的 断层的断面参数特征。对于只发育断层圈闭而未发育油气藏的情况，同样选取砂体在断面上 的投影作为提取参数的范围，只是该范围内所有数据均为不封闭断层的参数。
(6)断面封闭模型建立
在明确封闭断层和不封闭断层各参数的基础上，通过费希尔(Fisher)判别分析，建立基 于三参数的概率判别函数。函数建立过程中，随机抽取80％数据用于构建判别方程，剩余20％ 数据用于验证判别方程可靠性。
以渤海湾盆地济阳坳陷三合村油田为例进行说明。通过方法(5)共获得封闭断层与不封 闭断层断面参数255组，其中封闭断层数据170组，不封闭断层85组。选取其中202组数据 (80％)构建分类函数，分类函数系数见表1。将构建的分类函数分析剩余20％数据并将判 别结果与实际结果进行对比可以看出，正确率达到了94％(表2)。由此可看出通过建立断面 泥质含量、膨胀趋势及滑动趋势三维模型，综合地质分析提取断面参数，建立封闭判别模型 能够有效对碎屑岩地层中断层的封闭性进行准确评价。
表1封闭断层与不封闭断层判别方程系数表
| 不封闭 | 封闭 | |
| SGR(断面泥质含量) | 73.089 | 93.566 |
| ST(滑动趋势) | 173.572 | 237.193 |
| DT(膨胀趋势) | 19.602 | 8.693 |
| 常数 | -38.614 | -55.901 |
表2断层封闭预测结果与实际结果对比
| 预测封闭 | 预测不封闭 | |
| 实际封闭 | 34 | 1 |
| 实际不封闭 | 2 | 15 |
上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公 开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述 研究区具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的 意义。
Claims (5)
1.一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，包括以下步骤：
(1)构建断面、层位模型
基于三维地震数据体，对目标断层及地层进行精细解释，刻画断层面在空间不同位置处的走向、倾向、倾角及断层断距特征；
(2)地层泥质含量模型
利用测井中反映泥质含量的GR曲线，对地层泥质含量进行计算，获得单井的泥质分布特征，通过插值的方法，对三维空间的泥质含量特征进行刻画；
(3)地应力场与地层流体压力模型
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向及变化梯度、垂向应力的变化梯度与地层流体压力；
(4)断面参数刻画
计算的参数包括断面泥质含量SGR、滑动趋势ST和膨胀趋势DT；其中，断面泥质含量用来表征断层的岩性特征，滑动趋势和膨胀趋势用来表征断面的应力特征；断面泥质含量为沿断层滑过计算点的累计泥质含量，泥质含量越高，则认为断层岩泥质含量越高，封闭性越强，反之越弱；滑动趋势用来描述断层两盘沿断面滑动可能性的大小，滑动趋势越大，断层越趋于滑动，趋于不稳定，封闭能力越弱，反之，断层趋于稳定，封闭能力强；膨胀趋势用来表征断面的紧闭程度，膨胀趋势越大，断层趋于张开，封闭性越弱，反之，断层越紧闭，趋于封闭；
(5)评价参数提取
根据勘探实际情况，对典型断层型油气藏及不含油气的断层圈闭进行精细解剖，确定封闭断层与不封闭断层发育范围，提取断面参数信息；其中，对于典型油气藏，基于测井解释、试油分析明确油气分布层位，通过地震属性刻画油气空间范围，特别是油气藏与断层相接触的范围；对比分析断面封闭油气藏范围及储层范围，确定封闭断面范围及不封闭断面范围，进而提取封闭断层与不封闭断层断面参数；对于断层圈闭，首先通过成藏条件分析，明确储层发育特征、下部烃源岩的供烃能力以及顶部盖层封闭能力；选取发育储层、下部烃源岩具有供烃能力、顶部盖层具有较强封闭能力，油气未成藏主要由于断层封闭性弱，提取断面与储层相交范围内的断面参数作为不封闭断层参数特征；在参数提取过程中，以断面和储层相交区域为范围，根据实际情况即工区大小及地震道间隔，每隔数个地震道，提取断层的断面泥质含量、滑动趋势、膨胀趋势参数；
(6)断面封闭模型建立
在明确封闭断层和不封闭断层各参数的基础上，通过费希尔Fisher判别分析，建立基于三参数的概率判别函数；函数建立过程中，随机抽取80％数据用于构建判别方程，剩余20％数据用于验证判别方程可靠性。
2.根据权要求1所述的一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，所述步骤(1)中，利用地震解释软件，对三维地震数据体以3～5道为间隔进行解释，实现对研究区断层、地层的高精度刻画，从而重点获得地层的空间分布，断层的走向、倾向、倾角、断距特征；其中，地震采集间隔根据三维数据体精度及断层发育规模进行调整，地震采样间距大时减小解释间隔，断层规模小时也需减小地震解释时的间隔，反之，增大解释间隔；总体目标为以较高效率对目标断层、地层的展布特征进行刻画，构建三维地质模型。
3.根据权要求1或2所述的一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，所述步骤(2)中，首先利用测井曲线对单井泥质含量的分布进行计算，泥质含量计算主要用公式：
第三系地层Vsh＝0.083×(23.7×Igr-1) (公式1)
前第三系地层Vsh＝0.33×(22.1×Igr-1) (公式2)
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数；
针对不同时代地层，选择公式1或公式2进行计算；计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到；在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
4.根据权要求3所述的一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，所述步骤(3)中，水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征，通过取平均值的方法获得研究区水平最大主应力方向；水平最小主应力与水平最大主应力角度差为90°；
垂向主应力大小通过对岩石密度的积分获得，计算公式为：
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m；
水平最大与最小主应力梯度主要依靠偶极子声波测井值对地应力的响应获得，计算公式为：
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，单位％；φc为临界孔隙度，单位％；
在有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深度水平最大、最小主应力值；
地层流体压力计算主要针对地层中存在的超压，利用等效法进行计算，公式为：
Pz＝γw+Ze+γr(Z-Ze)＝ρrgZ-(ρr-ρw)gZe (公式8)
式中，PZ为欠压实泥岩的地层流体压力，单位MPa；γw为等效深度的静水压力梯度，单位MPa/m；Ze为欠压实泥岩所对应的等效深度，单位m；γr为埋藏深度与等效深度间的静岩压力梯度，单位MPa/m；Z为欠压实泥岩的埋深，单位m；ρr为埋藏深度与等效深度间的岩层平均密度，单位kg/m3；g为重力加速度，取9.8N/kg；ρw为地层水密度，取1000kg/m3；△t为欠压实泥岩的声波时差，单位μs/m；△t0为地表的原始声波时差，单位μs/m；C为正常压实泥岩的压实系数，单位m-1；
通过以上计算，获得水平最大、最小主应力大小及方向，垂向应力变化梯度以及地层流体压力特征。
5.根据权要求4所述的一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，所述步骤(4)中断面泥质含量计算公式为：
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断距，单位m；Vshi为步骤(2)中空间插值得到的紧邻断层的地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得；
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-PZ) (公式12)
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
式中，ST为断面滑动趋势，小数；DT为膨胀趋势，小数；τ为作用于断面上的剪应力，单位MPa；δn为作用在断面上的正应力，单位MPa；PZ为孔隙流体压力，单位MPa，利用公式8计算；其中，δH为水平最大主应力，MPa；δh为水平最小主应力，MPa；σv为垂向主应力，MPa；αn为断面法线与σH的夹角，单位°；βn为断面法线与σh的夹角，单位°；γn为断面法线与σv的夹角，单位°。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111484074.3A CN114185083B (zh) | 2021-12-07 | 2021-12-07 | 一种碎屑岩地层中断层封闭性的定量评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111484074.3A CN114185083B (zh) | 2021-12-07 | 2021-12-07 | 一种碎屑岩地层中断层封闭性的定量评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114185083A true CN114185083A (zh) | 2022-03-15 |
| CN114185083B CN114185083B (zh) | 2022-11-22 |
Family
ID=80603596
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111484074.3A Active CN114185083B (zh) | 2021-12-07 | 2021-12-07 | 一种碎屑岩地层中断层封闭性的定量评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114185083B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115630132A (zh) * | 2022-12-13 | 2023-01-20 | 华夏天信物联科技有限公司 | 一种用于确定地层中空范围的数据处理方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
| CN103675906A (zh) * | 2012-09-19 | 2014-03-26 | 中国石油化工股份有限公司 | 用于提高复杂碎屑岩储层预测精度的多控储层预测方法 |
| CN105388530A (zh) * | 2015-12-04 | 2016-03-09 | 中国石油化工股份有限公司 | 流线法油气运移模拟方法 |
| CN107917865A (zh) * | 2016-10-11 | 2018-04-17 | 中国石油化工股份有限公司 | 一种致密砂岩储层多参数渗透率预测方法 |
| CN108798655A (zh) * | 2018-05-23 | 2018-11-13 | 成都理工大学 | 一种煤层气逸散的三气合采实验装置 |
-
2021
- 2021-12-07 CN CN202111484074.3A patent/CN114185083B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120195165A1 (en) * | 2011-01-31 | 2012-08-02 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
| CN103675906A (zh) * | 2012-09-19 | 2014-03-26 | 中国石油化工股份有限公司 | 用于提高复杂碎屑岩储层预测精度的多控储层预测方法 |
| CN105388530A (zh) * | 2015-12-04 | 2016-03-09 | 中国石油化工股份有限公司 | 流线法油气运移模拟方法 |
| CN107917865A (zh) * | 2016-10-11 | 2018-04-17 | 中国石油化工股份有限公司 | 一种致密砂岩储层多参数渗透率预测方法 |
| CN108798655A (zh) * | 2018-05-23 | 2018-11-13 | 成都理工大学 | 一种煤层气逸散的三气合采实验装置 |
Non-Patent Citations (1)
| Title |
|---|
| 能源，等: "克拉苏构造带克深段构造特征及其石油地质意义", 《中国石油勘探》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115630132A (zh) * | 2022-12-13 | 2023-01-20 | 华夏天信物联科技有限公司 | 一种用于确定地层中空范围的数据处理方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114185083B (zh) | 2022-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dell'Arciprete et al. | Comparison of three geostatistical methods for hydrofacies simulation: a test on alluvial sediments | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| US10495782B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| AlRassas et al. | CO2 storage capacity estimation under geological uncertainty using 3-D geological modeling of unconventional reservoir rocks in Shahejie Formation, block Nv32, China | |
| Greenberg et al. | Geologic carbon storage at a one million tonne demonstration project: Lessons learned from the Illinois Basin–Decatur Project | |
| Elahi et al. | Dynamic fracture characterization from tracer-test and flow-rate data with ensemble Kalman filter | |
| Wimmers et al. | Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir | |
| Cerveny et al. | Reducing uncertainty with fault-seal analysis | |
| Benetatos et al. | Coping with uncertainties through an automated workflow for 3D reservoir modelling of carbonate reservoirs | |
| Deng et al. | An integrated workflow for reservoir modeling and flow simulation of the nikanassin tight gas reservoir in the western Canada sedimentary basin | |
| Schwartz | Fracture pattern characterization of the Tensleep formation, Teapot Dome, Wyoming | |
| CN114185083B (zh) | 一种碎屑岩地层中断层封闭性的定量评价方法 | |
| Balch et al. | Recent developments in modeling: Farnsworth Texas, CO 2 EOR carbon sequestration project | |
| Krásny et al. | Introductory paper: Hydrogeology of fractured rocks from particular fractures to regional approaches: State-of-the-art and future challenges | |
| Rashid et al. | The effect of fracturing on permeability in carbonate reservoir rocks | |
| Jasim et al. | Specifying quality of a tight oil reservoir through 3-d reservoir modeling | |
| Jambayev | Discrete fracture network modeling for a carbonate reservoir | |
| Kalantari-Dahaghi | Systematic approach to numerical simulation and modelling of shale gas reservoirs | |
| Esmaeilpour et al. | Permeability and Water Saturation Characterization and Prediction in Wellington Oil Field Using Core Analysis and Seismic Inversion | |
| Wallace | Use of 3-dimensional dynamic modeling of CO₂ injection for comparison to regional static capacity assessments of Miocene sandstone reservoirs in the Texas State Waters, Gulf of Mexico | |
| Slatt | Geologic controls on reservoir quality | |
| EP2912582B1 (en) | System, method and computer program product for evaluating and ranking geobodies using a euler characteristic | |
| Sharaf et al. | Reservoir description of Upper Cretaceous concretion-rich sandstone, Mississippi, USA: example from Tinsley Field | |
| CN108195669B (zh) | 油藏围压条件下岩石静态力学参数的矫正及预测方法 | |
| Lønøy et al. | Paleokarst reservoir architecture derived from forward modelling of modern cave systems-a case study |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20221101 Address after: Floor 33, Unit 1, Building 1, Diamond Plaza, No. 10, Jianshe Road, Chenghua District, Chengdu, Sichuan 610051 Applicant after: CHENGDU NORTH PETROLEUM EXPLORATION AND DEVELOPMENT TECHNOLOGY Co.,Ltd. Address before: Three road 610000 Sichuan city of Chengdu province Chenghua District Erxian Qiaodong No. 1 Applicant before: Chengdu University of Technology |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |