CN111650644B - 盐胶结砂岩定量预测方法 - Google Patents
盐胶结砂岩定量预测方法 Download PDFInfo
- Publication number
- CN111650644B CN111650644B CN202010650463.8A CN202010650463A CN111650644B CN 111650644 B CN111650644 B CN 111650644B CN 202010650463 A CN202010650463 A CN 202010650463A CN 111650644 B CN111650644 B CN 111650644B
- Authority
- CN
- China
- Prior art keywords
- salt
- sandstone
- wave impedance
- cemented
- amplitude ratio
- 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
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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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)盐胶结砂岩综合预测:结合振幅比和相对纵波阻抗两种方法的预测结果、盐胶结砂岩形成模式及地质认识对盐胶结进行综合预测,得到盐胶结砂岩发育范围综合预测图。本发明避免了传统预测方法的多解性强及难以预测盐含量的问题,预测结果更具有地质合理性。
Description
技术领域
本发明涉及油气地球物理勘探技术领域,具体是关于一种盐胶结砂岩定量预测方法。
背景技术
胶结作用是一种成岩作用,是使储层孔隙度降低的重要因素。根据胶结物的不同,常见的胶结类型主要有泥质胶结、钙质胶结、硅质胶结、铁质胶结。然而,在近些年西非、北海、西伯利亚东部等一些油气田发现中,盐胶结成为破坏砂岩储层使其致密的主要因素之一。盐胶结砂岩中石盐充填了孔隙,对储层破坏作用远大于其他胶结类型。因此,盐胶结分布直接影响着油气田储量分布及勘探和开发策略。
通常基于叠后资料开展储层预测方法有两种,一种是振幅属性的方法,一种是叠后纵波阻抗反演的方法。盐胶结发育区,储层段地层受盐含量、砂岩厚度、岩性组合等多种因素影响,地震响应特征复杂,无/低盐胶结砂岩和盐胶结致密砂岩的振幅响应阈值门限难以确定,因此振幅属性预测方法存在较大不确定性;而纵波阻抗受正常压实作用影响,随着埋深增加而增大,当目的层埋深相差较大时,从已钻井确定的纵波阻抗阈值门限难以区分存在埋深差异较大的盐胶结砂岩和无盐胶结砂岩及泥岩地层。两种现有储层预测技术用于识别盐胶结作用范围预测多解性强,且难以定量预测盐含量,亟待在潜在盐胶结发育区研发盐胶结砂岩定量预测方法,提高储层预测和钻探成功率。
发明内容
针对上述问题,本发明的目的是提供一种适用于叠后资料条件下的盐胶结砂岩定量预测方法,以避免传统预测方法的多解性强及难以预测盐含量的问题。
为实现上述目的,本发明采取以下技术方案:
本发明所述的盐胶结砂岩定量预测方法,包括以下步骤:
1)通过振幅比方法定量预测盐胶结砂岩的分布和盐含量;
2)通过相对纵波阻抗方法定量预测盐胶结砂岩的分布、厚度和盐含量;
3)盐胶结砂岩综合预测:结合振幅比和相对纵波阻抗两种方法的预测结果、盐胶结砂岩形成模式及地质认识对盐胶结进行综合预测,得到盐胶结砂岩发育范围综合预测图。
所述的盐胶结砂岩定量预测方法,优选地,所述步骤1)包括如下子步骤:
11)建立振幅比与砂岩盐含量之间的定量关系模版:
基于流体/固体替代技术,获取不同盐含量砂岩弹性参数,在此基础上通过对不同砂岩盐含量、厚度及岩性组合的地质模型开展地震正演,建立振幅比与砂岩盐含量之间的定量关系模版;
12)确定目的层顶面反射背景振幅值,计算目的层顶面振幅比属性:
统计资料可靠区目的层顶面最强反射振幅概率密度分布,将概率密度最大处的振幅作为该层背景无盐胶结储层顶面反射振幅,用目的层顶面反射最强振幅除以该层背景无盐胶结储层顶面反射振幅,得到目的层顶面振幅比属性;
13)定量确定阈值门限,求取基于振幅比的盐胶结砂岩分布图:
将盐含量与振幅比关系模版中盐含量为2%时对应的振幅比作为阈值门限,振幅比大于阈值门限的区域判别为盐胶结砂岩发育区,得到振幅比属性预测的盐胶结砂岩分布图;
14)根据盐含量与振幅比关系模版求取基于振幅比的盐含量分布图。
所述的盐胶结砂岩定量预测方法,优选地,所述振幅比为目的层顶面最强反射振幅与该层背景无盐胶结储层顶面反射振幅的比值。
所述的盐胶结砂岩定量预测方法,优选地,所述步骤2)包括如下子步骤:
21)对已钻井实测纵波阻抗进行高通滤波得到相对纵波阻抗,高通滤波参数根据地震频带确定,地震有效频带低值频率作为高通滤波起始频率;通过埋深与相对纵波阻抗交会分析,将盐胶结砂岩相对纵波阻抗属性的阈值门限标定为零,即相对纵波阻抗大于零则判别为盐胶结砂岩;
22)基于流体/固体替代技术,建立砂岩盐含量与相对纵波阻抗关系模版;
23)通过叠后反演求取研究区纵波阻抗体,对纵波阻抗体采用与井上分析相同参数进行高通滤波并通过已钻井标定得到相对纵波阻抗体。
24)提取目的层层间最大相对纵波阻抗属性,将最大相对纵波阻抗大于零的区域判别为盐胶结发育区,得到基于相对纵波阻抗的盐胶结砂岩分布图;
25)将层间相对纵波阻抗大于零的网格累加厚度,得到盐胶结砂岩厚度分布图;
26)根据盐含量与相对纵波阻抗关系模板,求取基于相对纵波阻抗的盐含量分布图。
所述的盐胶结砂岩定量预测方法,优选地,所述综合预测包括以下步骤:
①当振幅比预测结果与相对纵波阻抗预测结果相差小,且在边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,存在卤水灌入通道,符合盐胶结砂岩形成模式和地质认识时,则认为定量预测的盐胶结砂岩分布、厚度和盐含量可靠程度高;
②当振幅比预测结果与相对纵波阻抗预测结果相差大,认为定量预测存在多解性;
③当存在地震资料品质差,不满足定量预测资料基础的区域,如果紧邻的边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,则推测存在卤水灌入通道,该资料差区发育盐胶结;如果紧邻的边界断层处解释岩性对接为泥-砂对接,则推测不存在卤水灌入通道,该资料差区不发育盐胶结。
本发明由于采取以上技术方案,其具有以下优点:1、振幅比属性相对于振幅属性更加稳定,方法适用于不同批次处理的不同能量资料,避免了因资料能量差异导致盐胶结砂岩阈值门限难以定量选取的问题。2、相对纵波阻抗属性去除了正常压实影响,较纵波阻抗属性能更准确的识别盐胶结砂岩。3、基于振幅比-盐含量关系模版和相对纵波阻抗-盐含量关系模版,不仅可以定量确定盐胶结砂岩阈值门限,同时可以定量求取盐含量。4、本发明基于振幅比和相对纵波阻抗的盐胶结预测方法均为定量预测,可信度高。5、两种属性方法对比验证,并结合盐胶结形成机制、地质认识等因素综合预测盐胶结发育范围,使得预测结果更符合地质合理性。
附图说明
图1为西非某气田的盐胶结模式图;
图2为西非某气田区的璧心样本砂岩岩石薄片的示意图;其中图2(a)展示的是无盐胶结砂岩薄片;图2(b)展示的是盐胶结砂岩薄片。
图3为盐胶结砂岩定量预测方法流程图;
图4为不同盐含量砂岩地质模型和振幅响应特征图;其中图4(a)是不同盐含量砂岩地质模型图;图4(b)是图4(a)模型对应的地震反射记录图;图4(c)为提取的砂岩顶面最大反射振幅曲线图;图4(d)为振幅比曲线图。
图5为楔形砂岩模型及地震振幅响应特征图;其中,图5(a)是楔形砂岩地质模型图;图5(b)和图5(c)分别为楔形砂岩地质模型中砂岩为盐胶结砂岩和无盐胶结含气砂岩情况的地震记录图;图5(d)为对应的不同厚度砂岩的最大振幅曲线图;图5(e)为不同砂岩厚度情况下的振幅比图;
图6为厚层砂岩地质模型及地震反射记录图;其中,图6(a)为厚砂岩地质模型图,图6(b)和图6(c)分别为厚砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录图;
图7为薄互层砂岩地质模型及地震反射记录图;其中,图7(a)为薄互层砂岩地质模型图,图7(b)和图7(c)分别为薄互层砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录图;
图8为单薄层砂岩地质模型及地震反射记录图;其中,图8(a)为单薄层砂岩地质模型图,图8(b)和图8(c)分别为单薄层砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录图;
图9为厚砂岩、薄互层砂岩和单薄层砂岩三种岩性组合盐胶结砂岩和无盐胶结砂岩情况的反射振幅及振幅比图;
图10为盐含量与振幅比定量关系模版图;
图11为西非某气田目的层砂岩顶面最小反射振幅概率密度分布图(负极性资料);
图12为振幅比属性定量预测的盐胶结砂岩分布图;
图13为振幅比属性定量预测的盐含量图;
图14(a)为埋深和纵波阻抗交会图;图14(b)是埋深与相对纵波阻抗交会图;
图15为盐含量与相对纵波阻抗定量关系模板;
图16为叠后反演剖面图;其中,上、中和下剖面分别对应深度域地震剖面、叠后反演得到的纵波阻抗剖面和相对纵波阻抗剖面;
图17为相对纵波阻抗属性定量预测的盐胶结分布图;
图18为相对纵波阻抗属性定量预测的盐胶结砂岩厚度分布图;
图19为相对纵波阻抗属性定量预测的盐含量图;
图20为结合振幅比属性、相对纵波阻抗属性、盐胶结形成机制和地质认识综合预测的盐胶结砂岩发育范围图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明,以便更清楚理解本发明的目的、特点和优点。应理解的是,附图所示的实施例并不是对本发明范围的限制,而只是为了说明本发明技术方案的实质精神。
如图1所示,断层活动造成盐岩与盐下砂岩对接,卤水由盐岩向砂岩灌入是形成盐胶结的必要条件,而盐岩中存在泥岩并且与砂岩对接时阻止了卤水向砂岩灌入,盐浓度、流体势能、油气充注影响盐胶结平面范围。
如图2所示,其中,图2(a)为无盐胶结砂岩薄片图;图2(b)为盐胶结砂岩薄片图。从图2(a)和图2(b)的薄片图显示,无盐胶结砂岩发育大量粒间孔隙,孔隙连通,被流体充填,流体剪切模量为零;而在盐胶结砂岩中,固体石盐几乎充满了粒间孔隙,孔隙不再连通,孔隙充填物剪切模量不再为零,影响砂岩整体的弹性特征,进而势必影响砂岩地震振幅响应特征。
如图3所示,本发明盐胶结砂岩定量预测提出了两种方法。一种是振幅比方法,一种是相对纵波阻抗方法,最后结合两种方法预测结果进行综合预测。其中的方法步骤为:
步骤一、通过振幅比方法定量预测盐胶结砂岩的分布和盐含量;
其中步骤一包括如下步骤:
11)建立振幅比与砂岩盐含量之间的定量关系模版:
基于流体/固体替代技术,获取不同盐含量砂岩弹性参数,在此基础上通过对不同砂岩盐含量、厚度及岩性组合的地质模型开展地震正演,建立振幅比与砂岩盐含量之间的定量关系模版。
12)确定目的层顶面反射背景振幅值,计算目的层顶面振幅比属性:
统计资料可靠区目的层顶面最强反射振幅概率密度分布,将概率密度最大处的振幅作为背景振幅,用目的层顶面反射最强振幅除以背景振幅,得到目的层顶面振幅比属性。
13)定量确定阈值门限,求取基于振幅比的盐胶结砂岩分布图:
将盐含量与振幅比关系模版中盐含量为2%时对应的振幅比作为阈值门限,振幅比大于阈值门限的区域判别为盐胶结砂岩发育区,得到振幅比属性预测的盐胶结砂岩分布图;
14)根据盐含量与振幅比关系模版求取基于振幅比的盐含量分布图。
需要说明的是,本发明的振幅比为目的层顶面最强(波峰或波谷)反射振幅与该层无盐胶结砂岩储层顶面反射背景振幅的比值;定义盐含量大于2%时为盐胶结砂岩。
步骤二、通过相对纵波阻抗方法定量预测盐胶结砂岩的分布、厚度和盐含量。
其中步骤二包括如下步骤:
21)对已钻井实测纵波阻抗进行高通滤波得到相对纵波阻抗,高通滤波参数根据地震频带确定,地震有效频带低值频率作为高通滤波起始频率;通过埋深与相对纵波阻抗交会分析,将盐胶结砂岩相对纵波阻抗属性的阈值门限标定为零,即相对纵波阻抗大于零则判别为盐胶结砂岩;
22)基于流体/固体替代技术,建立砂岩盐含量与相对纵波阻抗关系模版;
23)通过叠后反演求取研究区纵波阻抗体,对纵波阻抗体采用与井上分析相同参数进行高通滤波并通过已钻井标定得到相对纵波阻抗体。
24)提取目的层层间最大相对纵波阻抗属性,将最大相对纵波阻抗大于零的区域判别为盐胶结发育区,得到基于相对纵波阻抗的盐胶结砂岩分布图;
25)将层间相对纵波阻抗大于零的网格累加厚度,得到盐胶结砂岩厚度分布图;
26)根据盐含量与相对纵波阻抗关系模板,求取基于相对纵波阻抗的盐含量分布图。
步骤三、盐胶结砂岩综合预测:
结合振幅比和相对纵波阻抗两种方法的预测结果、盐胶结砂岩形成模式及地质认识对盐胶结进行综合预测,得到盐胶结砂岩发育范围综合预测图。
上述实施例中,优选地,所述综合预测包括以下步骤
①当振幅比预测结果与相对纵波阻抗预测结果基本一致(相差小),且在边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,存在卤水灌入通道,符合盐胶结砂岩形成模式和地质认识时,则认为定量预测的盐胶结砂岩分布、厚度和盐含量可靠程度高;
②当振幅比预测结果与相对纵波阻抗预测结果相差大,认为定量预测存在多解性;
③当存在地震资料品质差,不满足定量预测资料基础的区域,如果紧邻的边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,则推测存在卤水灌入通道,该资料差区发育盐胶结;如果紧邻的边界断层处解释岩性对接为泥-砂对接,则推测不存在卤水灌入通道,该资料差区不发育盐胶结。最终得到盐胶结发育范围综合预测图。
如图4所示,其中,图4(a)为不同盐含量砂岩地质模型图;图4(b)为图4(a)模型对应的地震反射记录图;图4(c)为提取的砂岩顶面最大反射振幅曲线图;图4(d)为振幅比曲线图。
从图中可以反映出,该不同盐含量地质模型中砂岩厚度为50米,总孔隙度为16%,盐含量从左到右逐渐增加,盐含量从0%逐渐增加到16%;反射振幅和振幅比均随着盐含量增加而增大。盐含量为2%时对应振幅比为2.3,将其作为振幅比识别盐胶结砂岩的阈值门限。振幅比大于2.3时判别为盐胶结砂岩,反之为无盐或少盐砂岩。
如图5所示,其中,图5(a)为楔形砂岩地质模型图;图5(b)和图5(c)分别为楔形砂岩地质模型中砂岩为盐胶结砂岩和无盐胶结含气砂岩情况的地震记录;图5(d)为对应的不同厚度砂岩的最大振幅曲线;图5(e)为不同砂岩厚度情况下的振幅比;
楔形模型正演表明,盐胶结具有强振幅反射特征,而无盐胶结砂岩具有弱反射振幅特征。不同砂岩厚度对应的振幅比变化很小,振幅比方法稳定。
如图6-9所示,其中,图6(a)为厚砂岩地质模型,图6(b)和图6(c)分别为厚砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录;图7(a)为薄互层砂岩地质模型,图7(b)和图7(c)分别为薄互层砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录;图8(a)为单薄层砂岩地质模型,图8(b)和图8(c)分别为单薄层砂岩为盐胶结砂岩情况和无盐胶结含气砂岩情况的地震记录;图9为厚砂岩、薄互层砂岩和单薄层砂岩三种岩性组合盐胶结砂岩和无盐胶结砂岩情况的反射振幅及振幅比。
从图6-9中的厚砂层、薄互层和单薄层地质模型正演表明,盐胶结砂岩振幅远大于相同岩性组合的无盐胶结含气砂岩,不同岩性组合的振幅比相差很小,振幅比方法稳定。
如图10所示,建立了盐含量与振幅比的定量解释模板,该模板与图4(d)横、纵坐标相反,用于通过振幅比定量解释盐含量。盐胶结砂岩和无、低盐胶结砂岩盐含量划分门限为2%,对应振幅比阈值门限为2.3;图中标识了壁心实测的一个盐胶结样本点,该处盐含量9.6%,振幅比6.5。前面基于流体/固体替代的地震正演验证了振幅比主要与盐含量有关,受砂岩厚度和岩性组合的影响较小,因此该解释模板稳定。
下面通过实施例对本发明中振幅比预测盐胶结方法做进一步介绍。
实施例1:
如图11所示,地质认为西非某气田盐胶结为局部发育,即大部分区域为未发生盐胶结的含气砂岩,因此取目的层顶面反射振幅概率密度最大处的振幅作为背景振幅。该例中背景振幅为-1500(负极性资料)。提取目的层顶面最强反射振幅(该区储层顶面为正反射系数,资料为负极性资料,因此提取顶面反射波谷处振幅值,即最小振幅值),将顶面反射振幅除以背景振幅,得到振幅比属性,进而利用如图10所示的盐含量与振幅比的关系模板,求取盐胶结砂岩的盐含量。
如图12-13所示,将振幅比大于2.3的区域识别出来,判别为盐胶结砂岩发育区,得到基于振幅比的盐胶结砂岩分布图。进一步根据图10中所示模板,求取得到了基于振幅比的盐含量分布图。
下面结合附图和实施例对本发明相对纵波阻抗定量预测盐胶结进行详细的描述。
如图14所示,其中,图14(a)为埋深和纵波阻抗交会图;图14(b)是埋深与相对纵波阻抗交会图;图14(a)和图14(b)中黑框内为盐胶结砂岩。
图14(a)的埋深与纵波阻抗交会图表明,多井正常压实趋势一致,盐胶结砂岩具有相对高纵波阻抗特征。但因纵波阻抗受埋深影响大,纵波阻抗阈值难以准确预测盐胶结,本发明提出利用相对纵波阻抗区分。即对井上纵波阻抗进行高通滤波得到相对纵波阻抗,高通滤波参数根据地震频带确定,地震有效频带低值频率作为高通滤波起始频率。通过图14(b)的埋深与相对纵波阻抗交会分析,将盐胶结砂岩相对纵波阻抗属性的阈值门限标定为零,即相对纵波阻抗大于零则判别为盐胶结砂岩。
如图15所示,基于目的层段典型盐胶结砂岩壁心样本分析数据和流体/固体替代技术,建立了盐含量与相对纵波阻抗关系模板。该模板埋深采用样本点处埋深,离散样本点无法采用高通滤波方法,因此通过坐标转换方法等价地将纵波阻抗转换为相对纵波阻抗。基于该模板可通过相对纵波阻抗求取盐含量。
下面通过实施例对本发明中相对纵波阻抗定量预测盐胶结方法做进一步介绍。
实施例2:
如图16所示,通过叠后反演求取研究区纵波阻抗体,对纵波阻抗采用与井上分析相同参数进行高通滤波并通过已钻井标定得到相对纵波阻抗体。
如图17-19所示,提取目的层层间最大相对纵波阻抗属性,将最大相对纵波阻抗大于零的区域判别为盐胶结发育区,得到相对纵波阻抗预测的盐胶结分布图。将层间相对纵波阻抗大于零的网格累计叠加厚度,得到盐胶结砂岩厚度分布图;根据图15所示的盐含量与相对纵波阻抗定量关系模版,求取基于相对纵波阻抗预的盐含量分布图。
如图20所示,结合振幅比和相对纵波阻抗两种属性预测结果、盐胶结砂岩形成机制及地质认识对盐胶结发育范围进行综合预测。通过速度异常和地震相解释了低速泥岩在边界断层处的发育范围,在没有泥岩封隔的地方,砂岩与盐岩对接,卤水灌入发育盐胶结。利用振幅比和相对纵波阻抗两种方法预测的盐胶结发育范围一致,互相得到验证,降低了多解性。预测结果与已钻井一致,且分布范围符合盐胶结形成机制及地质认识,为储量评价及井位研究提供了重要参考。
综上所述,本发明与现有技术具有明显的区别,其主要为以下几点:1)本发明提出的振幅比属性方法相对于现有振幅属性方法更加稳定,可适用于不同批次处理的不同能量资料,避免了因资料能量差异导致盐胶结砂岩阈值门限难以定量选取的问题。2、相对纵波阻抗属性去除了正常压实影响,较纵波阻抗属性更能准确的识别盐胶结砂岩。3、本发明使用定量解释模版,不仅可以定量确定盐胶结砂岩阈值门限,同时可以定量求取盐含量,降低现有技术定性预测的不确定性;4、两种属性方法对比验证,并结合盐胶结形成机制、地质认识等因素综合预测盐胶结发育范围,使得预测结果更具有地质合理性。
另外,需要说明的是:本发明中的流体/固体替代技术是指:利用标定好的盐胶结砂岩岩石物理等效模型实现流体/固体替代,获取石盐替换为流体或流体替换为石盐情况的弹性参数。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (3)
1.一种盐胶结砂岩定量预测方法,其特征在于,包括以下步骤:
1)通过振幅比方法定量预测盐胶结砂岩的分布和盐含量;
2)通过相对纵波阻抗方法定量预测盐胶结砂岩的分布、厚度和盐含量;
3)盐胶结砂岩综合预测:结合振幅比和相对纵波阻抗两种方法的预测结果、盐胶结砂岩形成模式及地质认识对盐胶结进行综合预测,得到盐胶结砂岩发育范围综合预测图;
所述步骤1)包括如下子步骤:
11)建立振幅比与砂岩盐含量之间的定量关系模版:
基于流体/固体替代技术,获取不同盐含量砂岩弹性参数,在此基础上通过对不同砂岩盐含量、厚度及岩性组合的地质模型开展地震正演,建立振幅比与砂岩盐含量之间的定量关系模版;
12)确定目的层顶面反射背景振幅值,计算目的层顶面振幅比属性:
统计资料可靠区目的层顶面最强反射振幅概率密度分布,将概率密度最大处的振幅作为该层背景无盐胶结储层顶面反射振幅,用目的层顶面反射最强振幅除以该层背景无盐胶结储层顶面反射振幅,得到目的层顶面振幅比属性;
13)定量确定阈值门限,求取基于振幅比的盐胶结砂岩分布图:
将盐含量与振幅比关系模版中盐含量为2%时对应的振幅比作为阈值门限,振幅比大于阈值门限的区域判别为盐胶结砂岩发育区,得到振幅比属性预测的盐胶结砂岩分布图;
14)根据盐含量与振幅比关系模版求取基于振幅比的盐含量分布图。
2.根据权利要求1所述的盐胶结砂岩定量预测方法,其特征在于,所述步骤2)包括如下子步骤:
21)对已钻井实测纵波阻抗进行高通滤波得到相对纵波阻抗,高通滤波参数根据地震频带确定,地震有效频带低值频率作为高通滤波起始频率;通过埋深与相对纵波阻抗交会分析,将盐胶结砂岩相对纵波阻抗属性的阈值门限标定为零,即相对纵波阻抗大于零则判别为盐胶结砂岩;
22)基于流体/固体替代技术,建立砂岩盐含量与相对纵波阻抗关系模版;
23)通过叠后反演求取研究区纵波阻抗体,对纵波阻抗体采用与井上分析相同参数进行高通滤波并通过已钻井标定得到相对纵波阻抗体;
24)提取目的层层间最大相对纵波阻抗属性,将最大相对纵波阻抗大于零的区域判别为盐胶结发育区,得到基于相对纵波阻抗的盐胶结砂岩分布图;
25)将层间相对纵波阻抗大于零的网格累加厚度,得到盐胶结砂岩厚度分布图;
26)根据盐含量与相对纵波阻抗关系模板,求取基于相对纵波阻抗的盐含量分布图。
3.根据权利要求1所述的盐胶结砂岩定量预测方法,其特征在于,所述综合预测包括以下步骤:
①当振幅比预测结果与相对纵波阻抗预测结果相差小,且在边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,存在卤水灌入通道,符合盐胶结砂岩形成模式和地质认识时,则认为定量预测的盐胶结砂岩分布、厚度和盐含量可靠程度高;
②当振幅比预测结果与相对纵波阻抗预测结果相差大,认为定量预测存在多解性;
③当存在地震资料品质差,不满足定量预测资料基础的区域,如果紧邻的边界断层处根据地震速度和地震相解释岩性对接为盐-砂对接,则推测存在卤水灌入通道,该资料差区发育盐胶结;如果紧邻的边界断层处解释岩性对接为泥-砂对接,则推测不存在卤水灌入通道,该资料差区不发育盐胶结。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010650463.8A CN111650644B (zh) | 2020-07-08 | 2020-07-08 | 盐胶结砂岩定量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010650463.8A CN111650644B (zh) | 2020-07-08 | 2020-07-08 | 盐胶结砂岩定量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111650644A CN111650644A (zh) | 2020-09-11 |
CN111650644B true CN111650644B (zh) | 2022-03-22 |
Family
ID=72347844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010650463.8A Active CN111650644B (zh) | 2020-07-08 | 2020-07-08 | 盐胶结砂岩定量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111650644B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114624776B (zh) * | 2020-12-14 | 2024-03-19 | 中国石油化工股份有限公司 | 断层胶结带影响断层地震响应特征变化的定量分析方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884696A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于钙质胶结作用影响的砂岩横波速度计算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5081613A (en) * | 1988-09-27 | 1992-01-14 | Applied Geomechanics | Method of identification of well damage and downhole irregularities |
CN103527184B (zh) * | 2013-10-28 | 2016-04-20 | 北京大学 | 一种白云岩储层的预测方法和系统 |
CN103852787B (zh) * | 2014-02-24 | 2016-08-24 | 长江大学 | 一种砂岩储层成岩地震相表征方法 |
CN105005074A (zh) * | 2015-06-23 | 2015-10-28 | 成都理工大学 | 利用频变地震反射系数识别含气储层的方法 |
BR112021000838A2 (pt) * | 2018-07-27 | 2021-04-13 | Baker Hughes Holdings Llc | Avaliação de cimento através da tubulação com o uso de métodos sísmicos |
-
2020
- 2020-07-08 CN CN202010650463.8A patent/CN111650644B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884696A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于钙质胶结作用影响的砂岩横波速度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111650644A (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Verma et al. | Estimation of total organic carbon and brittleness volume | |
US9244182B2 (en) | Method of assessing hydrocarbon source rock candidate | |
NO20180189A1 (en) | Seismic constrained discrete fracture network | |
WO2021048597A1 (en) | An integrated geomechanics model for predicting hydrocarbon and migration pathways | |
CN103643949A (zh) | 一种储层含油气的定量预测方法及装置 | |
CN111271057A (zh) | 一种测井识别气层方法 | |
CN107966730A (zh) | 识别断陷湖盆陡坡带近岸水下扇沉积相带的方法 | |
Eshimokhai et al. | Reservoir characterization using seismic and well logs data (a case study of Niger Delta) | |
CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
CN107831540A (zh) | 储层物性参数直接提取新方法 | |
CN112946782B (zh) | 一种致密油气储渗体地震精细刻画方法 | |
CN111650644B (zh) | 盐胶结砂岩定量预测方法 | |
CN113093286B (zh) | 一种冷泉发育区储层非均质性的反演方法 | |
Davies | Permeability Modelling of a Sandstone Reservoir in Parts of the Niger Delta | |
CN106481337B (zh) | 超压顶界面的预测方法 | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
Masoud et al. | Reservoir Characterization and Geostatistical Model of the Cretaceous and Cambrian-Ordovician Reservoir Intervals, Meghil Field, Sirte Basin, Libya | |
Jageler et al. | Use of well logs and dipmeters in stratigraphic-trap exploration: geologic exploration methods | |
CN113495293B (zh) | 油藏流体预测方法及装置 | |
CN114509812B (zh) | 一种碎屑岩储层甜点区预测方法、装置、介质及电子设备 | |
CN110297280A (zh) | 一种碳酸盐岩超压空间分布特征的预测方法 | |
Crowley et al. | The Stag Oilfield | |
Deabl et al. | Evaluation of Petrophysical Properties Interpretations from Log Interpretation for Tertiary Reservoir/Ajeel Field | |
Adlakha | Fracture analysis and modelling of the South Arne field | |
Pradhanitasari | Petrophysical characterizations and shale distribution analysis of Kissenda reservoir sands, Gabon |
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 |