CN114185083B - 一种碎屑岩地层中断层封闭性的定量评价方法 - Google Patents

一种碎屑岩地层中断层封闭性的定量评价方法 Download PDF

Info

Publication number
CN114185083B
CN114185083B CN202111484074.3A CN202111484074A CN114185083B CN 114185083 B CN114185083 B CN 114185083B CN 202111484074 A CN202111484074 A CN 202111484074A CN 114185083 B CN114185083 B CN 114185083B
Authority
CN
China
Prior art keywords
fault
section
stratum
formula
stress
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
Application number
CN202111484074.3A
Other languages
English (en)
Other versions
CN114185083A (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.)
Chengdu North Petroleum Exploration And Development Technology Co ltd
Original Assignee
Chengdu North Petroleum Exploration And Development Technology Co ltd
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 Chengdu North Petroleum Exploration And Development Technology Co ltd filed Critical Chengdu North Petroleum Exploration And Development Technology Co ltd
Priority to CN202111484074.3A priority Critical patent/CN114185083B/zh
Publication of CN114185083A publication Critical patent/CN114185083A/zh
Application granted granted Critical
Publication of CN114185083B publication Critical patent/CN114185083B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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. analysis, for interpretation, for correction
    • 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. analysis, for interpretation, for correction
    • 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. analysis, for interpretation, for correction
    • 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. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

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)
Figure GDA0003793437640000021
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数。
针对不同时代地层，选择公式1或公式2进行计算。计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到。在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
(3)地应力场与地层流体压力模型：
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向及变化梯度、垂向应力的变化梯度与地层流体压力。
上述步骤(3)中，水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征，不同井之间显示出的水平最大主应力方向可能存在一定差异，通过取平均值的方法获得研究区水平最大主应力方向。水平最小主应力与水平最大主应力角度差为90°。
由于密度测井资料通常十分全面，几乎包含了全部井段。因此，垂向主应力大小通过对岩石密度的积分获得，计算公式为：
Figure GDA0003793437640000022
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m。
水平最大与最小主应力主要依靠偶极子声波测井值计算获得，计算公式为：
Figure GDA0003793437640000031
Figure GDA0003793437640000032
Figure GDA0003793437640000033
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，单位％；φc为临界孔隙度(一般取40)。
由于成本导致偶极子声波测井通常针对特定层位，研究层位可能无相关数据。对此，在有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深度水平最大、最小主应力值。(除此方法外，在样品充裕的情况下，也可通过岩石破裂法、声发射法、油田现场井壁崩落资料、储层压裂改造资料获得水平最大与最小主应力数据。)
地层流体压力计算主要针对地层中存在的超压，利用等效法进行计算，公式为：
Pz＝γwZe+γr(Z-Ze)＝ρrgZ-(ρr-ρw)gZe (公式8)
Figure GDA0003793437640000034
Figure GDA0003793437640000035
式中，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)。其中，断面泥质含量用来表征断层的岩性特征，滑动趋势和膨胀趋势用来表征断面的应力特征。断面泥质含量为沿断层滑过计算点的累计泥质含量，泥质含量越高，则认为断层岩泥质含量越高，封闭性越强，反之越弱。滑动趋势用来描述断层两盘沿断面滑动可能性的大小，滑动趋势越大，断层越趋于滑动，趋于不稳定，封闭能力越弱，反之，断层趋于稳定，封闭能力强。膨胀趋势用来表征断面的紧闭程度，膨胀趋势越大，断层趋于张开，封闭性越弱，反之，断层越紧闭，趋于封闭。
断面泥质含量计算公式为：
Figure GDA0003793437640000041
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断距，单位m。Vshi为步骤(2)中空间插值得到的紧邻断层地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得。
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-PZ) (公式12)
Figure GDA0003793437640000042
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
Figure GDA0003793437640000043
式中，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)
Figure GDA0003793437640000071
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数。
针对不同时代地层，选择公式1或公式2进行计算。计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到。在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
(3)地应力场与地层流体压力模型：
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向及变化梯度、垂向应力的变化梯度。
上述步骤(3)中，水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征，不同井之间显示出的水平最大主应力方向可能存在一定差异，通过取平均值的方法获得研究区水平最大主应力方向。水平最小主应力与水平最大主应力角度差为90°。
由于密度测井资料通常十分全面，几乎包含了全部井段。因此，垂向主应力大小通过对岩石密度的积分获得，计算公式为：
Figure GDA0003793437640000072
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m。
水平最大与最小主应力主要依靠偶极子声波测井值计算获得，计算公式为：
Figure GDA0003793437640000073
Figure GDA0003793437640000074
Figure GDA0003793437640000081
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，％；φc为临界孔隙度(一般取40％)。
由于成本导致偶极子声波测井通常针对特定层位，研究层位可能无相关数据。对此，在有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深度水平最大、最小主应力值，如图2a和图2b。除此方法外，在样品充裕的情况下，也可通过岩石破裂法、声发射法、油田现场井壁崩落资料、储层压裂改造资料获得水平最大与最小主应力数据。
地层流体压力计算主要针对地层中存在的超压，利用等效法进行计算，公式为：
Pz＝γwZe+γr(Z-Ze)＝ρrgZ-(ρr-ρw)gZe (公式8)
Figure GDA0003793437640000082
Figure GDA0003793437640000083
式中，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，计算公式为：
Figure GDA0003793437640000091
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断距，单位m。Vshi为步骤(2)中空间插值得到的紧邻断层地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得。
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-Pz) (公式12)
Figure GDA0003793437640000092
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
Figure GDA0003793437640000093
式中，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 (4)

1.一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，包括以下步骤：
(1)构建断面、层位模型
基于三维地震数据体，对目标断层及地层进行精细解释，刻画断层面在空间不同位置处的走向、倾向、倾角及断层断距特征；
(2)地层泥质含量模型
利用测井中反映泥质含量的GR曲线，对地层泥质含量进行计算，获得单井的泥质分布特征，通过插值的方法，对三维空间的泥质含量特征进行刻画；
(3)地应力场与地层流体压力模型
通过实际测试或测井资料计算地应力场特征，包括水平最大、最小主应力的大小和方向及变化梯度、垂向应力的变化梯度与地层流体压力；
水平最大主应力方向利用偶极子声波测井中的快横波方向进行表征，通过取平均值的方法获得研究区水平最大主应力方向；水平最小主应力与水平最大主应力角度差为90°；
垂向主应力大小通过对岩石密度的积分获得，计算公式为：
Figure FDA0003793437630000011
式中：P为目的深度的垂向应力，单位MPa；ρ(z)为岩层密度，由密度测井资料获得，单位kg/m3；g为重力加速度，9.8N/kg；h为目的层深度，单位m；
水平最大与最小主应力梯度依靠偶极子声波测井值对地应力的响应获得，计算公式为：
Figure FDA0003793437630000012
Figure FDA0003793437630000013
Figure FDA0003793437630000014
式中：σH为水平最大主应力，单位MPa；σh为水平最小主应力，单位MPa；βH和βh分别为水平最大应力和最小应力方向上的构造系数，无量纲；v和α分别为岩石的泊松比和Biot系数，无量纲；φ为岩石孔隙度，此处使用测井解释孔隙度，单位％；φc为临界孔隙度，单位％；
在有数据层段计算水平最大、最小主应力后，结合深度拟合垂向上变化梯度，从而获得不同深度水平最大、最小主应力值；
地层流体压力计算针对地层中存在的超压，利用等效法进行计算，公式为：
Figure FDA0003793437630000021
Figure FDA0003793437630000022
Figure FDA0003793437630000023
式中，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)断面参数刻画
计算的参数包括断面泥质含量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)
Figure FDA0003793437630000031
式中：GRlog为伽马测井读数，单位API；GRmax为纯泥岩时的伽马测井读数，单位API；GRmin为纯砂岩时的伽马测井读数，单位API；Igr为计算泥质含量，小数；Vsh为实际泥质含量，小数；
针对不同时代地层，选择公式1或公式2进行计算；计算时，纯泥岩和纯砂岩GR的获取，根据厚层纯泥岩、纯砂岩的平均值计算得到；在获得单井泥质含量的基础上，以步骤1中的地层解释结果为约束，进行三维空间插值，实现对地层中泥质含量的刻画。
4.根据权利要求1所述的一种碎屑岩地层中断层封闭性的定量评价方法，其特征在于，所述步骤(4)中断面泥质含量计算公式为：
Figure FDA0003793437630000041
式中，SGR为计算点的断面泥质含量，小数；n为划过目的点的砂泥岩层层数；ΔZi为滑过目的点的第i层岩层厚度，单位m；Vshi为滑过目的点的第i层岩层泥质含量，H为断层断距，单位m；Vshi为步骤(2)中空间插值得到的紧邻断层的地层中不同深度的泥质含量，断距H由步骤(1)的解释结果获得；
滑动趋势及膨胀趋势的计算公式为：
ST＝τ/(δn-PZ) (公式12)
Figure FDA0003793437630000042
σn＝σHcos2αn+σhcos2βn+σvcos2γn (公式14)
Figure FDA0003793437630000043
式中，ST为断面滑动趋势，小数；DT为膨胀趋势，小数；τ为作用于断面上的剪应力，单位MPa；δn为作用在断面上的正应力，单位MPa；PZ为孔隙流体压力，单位MPa，利用公式8计算；其中，δH为水平最大主应力，MPa；δh为水平最小主应力，MPa；σv为垂向主应力，MPa；αn为断面法线与σH的夹角，单位°；βn为断面法线与σh的夹角，单位°；γn为断面法线与σv的夹角，单位°。
CN202111484074.3A 2021-12-07 2021-12-07 一种碎屑岩地层中断层封闭性的定量评价方法 Active CN114185083B (zh)

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 CN114185083A (zh) 2022-03-15
CN114185083B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115630132B (zh) * 2022-12-13 2023-03-17 华夏天信物联科技有限公司 一种用于确定地层中空范围的数据处理方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103675906A (zh) * 2012-09-19 2014-03-26 中国石油化工股份有限公司 用于提高复杂碎屑岩储层预测精度的多控储层预测方法
CN107917865A (zh) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 一种致密砂岩储层多参数渗透率预测方法
CN108798655A (zh) * 2018-05-23 2018-11-13 成都理工大学 一种煤层气逸散的三气合采实验装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8861309B2 (en) * 2011-01-31 2014-10-14 Chevron U.S.A. Inc. Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context
CN105388530B (zh) * 2015-12-04 2020-05-19 中国石油化工股份有限公司 流线法油气运移模拟方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103675906A (zh) * 2012-09-19 2014-03-26 中国石油化工股份有限公司 用于提高复杂碎屑岩储层预测精度的多控储层预测方法
CN107917865A (zh) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 一种致密砂岩储层多参数渗透率预测方法
CN108798655A (zh) * 2018-05-23 2018-11-13 成都理工大学 一种煤层气逸散的三气合采实验装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
克拉苏构造带克深段构造特征及其石油地质意义;能源，等;《中国石油勘探》;20131231(第2期);1-6 *

Also Published As

Publication number Publication date
CN114185083A (zh) 2022-03-15

Similar Documents

Publication Publication Date Title
Satter et al. Reservoir engineering: the fundamentals, simulation, and management of conventional and unconventional recoveries
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
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
Alfataierge 3D modeling and characterization of hydraulic fracture efficiency integrated with 4D/9C time-lapse seismic interpretations in the Niobrara Formation, Wattenberg Field, Denver Basin
Deng et al. An integrated workflow for reservoir modeling and flow simulation of the nikanassin tight gas reservoir in the western Canada sedimentary basin
CN114185083B (zh) 一种碎屑岩地层中断层封闭性的定量评价方法
CN111206921A (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
EP2912582B1 (en) System, method and computer program product for evaluating and ranking geobodies using a euler characteristic
Rashid et al. The effect of fracturing on permeability in carbonate reservoir rocks
CN108195669B (zh) 油藏围压条件下岩石静态力学参数的矫正及预测方法
Sharaf et al. Reservoir description of Upper Cretaceous concretion-rich sandstone, Mississippi, USA: example from Tinsley Field
CN111608649A (zh) 外源补给型页岩气勘探有利区预测方法
Parsegov Physics-Based Forward Modeling of Multistage Hydraulic Fracturing in Unconventional Plays
Levon Workflow Development and Sensitivity Investigation of Offset Well-to-Well Interference through 3D Fracture Modeling and Reservoir Simulation in the Denver-Julesburg Basin
Sadeghvishkaei Modelling of geomechanics for informed hydraulic fracturing operations
Okwen et al. Storage Complex Modeling for CarbonSAFE Illinois-Macon County (Task 8)
Deng et al. Reservoir characterization and flow simulation of a low-permeability gas reservoir: an integrated approach for modelling the tommy lakes gas field
Branning Modeling Unconventional Water Production in the Northeast Elm Coulee Field, Bakken Formation
Arigbe Uncertainty reduction in reservoir parameters prediction from multiscale data using machine learning in deep offshore reservoirs.
Hamoud Reservoir Simulation Study of Libya’s O-Field NC-115 in Murzuq Basin

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