CN104181594A - 一种复杂沉积体边界的确定方法 - Google Patents

Abstract

本发明提供一种复杂沉积体边界的确定方法，以解决现有技术中对剥蚀和超覆识别精度低的问题；包括以下步骤：常规地震剖面追踪、地震属性分析和三角函数预测法。本发明应用资料更加广泛，能够识别出更小级别的剥蚀和超覆，减小油田开发部署面临的风险。且在平面上识别精度高，通过实际应用误差范围降低到100-300m。

Description

一种复杂沉积体边界的确定方法

技术领域

本发明涉及油气勘探开发储层保护技术领域，特别涉及一种复杂沉积体边界的确定方法。 

背景技术

由于构造的变化，地层在沉积过程中会经常出现剥蚀和超覆的现象，这对于油气分布与储量规模影响巨大。目前对于剥蚀和超覆一般是利用地震追踪的方法进行识别，但是精度较低，误差范围可以达到1-3Km。随着地球物理技术的发展，目前逐渐将地震属性信息加入，辅助地震追踪进行剥蚀和超覆的预测，精度虽然有所提高，但是提高幅度不大。 

发明内容

针对上述现有技术不足，本发明提供一种复杂沉积体边界的确定方法，以解决现有技术中对剥蚀和超覆识别精度低的问题。 

本发明的技术方案是这样实现的： 

一种复杂沉积体边界的确定方法，包括以下步骤： 

(1)高精度三维地震剖面追踪：首先对研究区所有钻井资料进行统层对比分析，在单井上识别出超覆和剥蚀的深度点；在完钻井层位标定基础上，通过高精度三维地震剖面同相轴追踪，确定地层超覆和剥蚀线的大体位置； 

(2)地震属性分析：通过对平均瞬时相位、平均振幅地震属性资料进行分析研究区储层的接触关系，进而判断地层超覆和剥蚀线的位置范围； 

(3)三角函数预测法：在地层构造解释的基础上，获得地层倾角，进而利用三角函数预测法进行地层超覆和剥蚀位置的判断。 

其中，优选地，所述三角函数预测法包括以下步骤： 

a.超覆线位置确定 

Ng底的倾角为α；1井与2井之间距离为△L，1井与2井Ng1到Ng底的距离分别为△H₁和△H₂； 

为了能够降低计算的误差，假设层位超覆位置距离1井与2井分别为x₁和x₂+△L，最终取x₁与x₂的平均值作为距离1井的层位超覆位置； 

ΔH₁＝x₁·tgα 

即 $x_1=\frac{\mathrm{\Δ}{H}_1}{\mathrm{tg\α}}$

ΔH_２＝(x_２＋ΔL)·tgα 

即 $x_2=\frac{\mathrm{\Δ}{H}_2}{\mathrm{tg\α}}-\mathrm{\ΔL}$

于是可得： 

$x=\frac{x_1+x_2}{2}=\frac{\mathrm{\Δ}{H}_1+\mathrm{\Δ}{H}_2}{2\mathrm{tg\α}}-\frac{\mathrm{\ΔL}}{2}$

X为超覆线的位置； 

b.剥蚀线位置确定 

Ng底的倾角为α′；ES3上的倾角为β；1井与2井之间距离为△L′，1井与2井Ng底到ES3上顶的距离分别为△H₁′和△H₂′； 

为了能够降低计算的误差，假设ES3上的剥蚀位置距离1井与2井分别为x₁′和x₂′+△L′，ES3上的剥蚀位置在1井与2井上的水平投影距离Ng底分别为y₁和y₂，最终取x₁′与x₂′的平均值作为距离1井的ES3上的剥蚀位置； 

y₁＝x₁′·tgα′             (1) 

$\mathrm{tg\β}=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{{x_1}^{\′}}---\left(2\right)$

即 ${x_1}^{\′}=\frac{y_1+\mathrm{\Δ}{H_1}^{\′}}{\mathrm{tg\β}}---\left(3\right)$

将式(3)代入式(1)可得： 

$y_1=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(4\right)$

于是可得： $y_1=\frac{\mathrm{\Δ}{H_1}^{\′}\·{\mathrm{tg\α}}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

进而可得： ${x_1}^{\′}=\frac{{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

y_２＝(x_２′+ΔL′)·tgα′      (5) 

$\mathrm{tg\β}=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{{x_2}^{\′}+{\mathrm{\ΔL}}^{\′}}---\left(6\right)$

即 ${x_2}^{\′}=\frac{{y_2}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}-{\mathrm{\ΔL}}^{\′}---\left(7\right)$

将式(7)代入式(5)可得： 

$y_2=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(8\right)$

于是可得： $y_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}\·{\mathrm{tg\α}}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

进而可得： $X_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}-{\mathrm{\Δl}}^{\′}$

于是可得： 

$x^{\′}=\frac{{x_1}^{\′}+{x_2}^{\′}}{2}=\frac{{{\mathrm{\ΔH}}_1}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{2(\mathrm{tg\β}-{\mathrm{tg\α}}^{\′})}-\frac{{\mathrm{\ΔL}}^{\′}}{2}$

X′为剥蚀线的位置。 

本发明的有益效果： 

1、本发明应用资料更加广泛，能够识别出更小级别的剥蚀和超覆，减小油田开发部署面临的风险。 

2、本发明在平面上识别精度高，通过实际应用误差范围降低到100-300m。 

说明书附图 

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。 

图1地层接触关系示意图； 

图2地层超覆示意图； 

图3地层剥蚀示意图。 

具体实施方式

下面对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。 

如图1、图2和图3所示，假如有一套沉积体如图1所示，上覆地层馆陶为超覆沉积，下伏地层沙一段及沙三段均有剥蚀。本实施例对于这种复杂沉积体的边界建立了一套确定方法。文字中Ng1、Ng2、Ng底以及Es1、Es2和Es3等是地层的名称，也可以统一使用1号地层、2号地层等等表示。 

根据实钻井情况，通过常规地震剖面追踪、地震属性分析，确定超覆及剥蚀位置的范围；结合钻井分层资料以及三角函数法，综合确定超覆及剥蚀位置点； 综合这两个步骤最终得到超覆线及剥蚀线的具体位置。 

一种复杂沉积体边界的确定方法，包括以下步骤： 

(1)高精度三维地震剖面追踪：首先对研究区所有钻井资料进行统层对比分析，在单井上识别出超覆和剥蚀的深度点；在完钻井层位标定基础上，通过高精度三维地震剖面同相轴追踪，确定地层超覆和剥蚀线的大体位置； 

(2)地震属性分析：通过对平均瞬时相位、平均振幅地震属性资料进行分析研究区储层的接触关系，进而判断地层超覆和剥蚀线的位置范围； 

(3)三角函数预测法：在地层构造解释的基础上，获得地层倾角，进而利用三角函数进行地层超覆和剥蚀位置的判断。 

其中，所述三角函数预测法包括以下步骤： 

a.超覆线位置确定 

Ng底的倾角为α；1井与2井之间距离为△L，1井与2井Ng1到Ng底的距离分别为△H₁和△H₂； 

为了能够降低计算的误差，假设层位超覆位置距离1井与2井分别为x₁和x₂+△L，最终取x₁与x₂的平均值作为距离1井的层位超覆位置； 

ΔH₁＝x₁·tgα 

即 $x_1=\frac{\mathrm{\Δ}{H}_1}{\mathrm{tg\α}}$

ΔH_２＝(x_２+ΔＬ)·tgα 

即 $x_2=\frac{\mathrm{\Δ}{H}_2}{\mathrm{tg\α}}-\mathrm{\ΔL}$

于是可得： 

$x=\frac{x_1+x_2}{2}=\frac{\mathrm{\Δ}{H}_1+\mathrm{\Δ}{H}_2}{2\mathrm{tg\α}}-\frac{\mathrm{\ΔL}}{2}$

X为超覆线的位置； 

b.剥蚀线位置确定 

Ng底的倾角为α′；ES3上的倾角为β；1井与2井之间距离为△L′，1井与2井Ng底到ES3上顶的距离分别为△H₁′和△H₂′； 

为了能够降低计算的误差，假设ES3上的剥蚀位置距离1井与2井分别为x₁′和x₂′+△L′，ES3上的剥蚀位置在1井与2井上的水平投影距离Ng底分别为y₁和y₂，最终取x₁′与x₂′的平均值作为距离1井的ES3上的剥蚀位置； 

y₁＝x₁′·tgα′         (1) 

$\mathrm{tg\β}=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{{x_1}^{\′}}---\left(2\right)$

即 ${x_1}^{\′}=\frac{y_1+\mathrm{\Δ}{H_1}^{\′}}{\mathrm{tg\β}}---\left(3\right)$

将式(3)代入式(1)可得： 

$y_1=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(4\right)$

于是可得： $y_1=\frac{{{\mathrm{\ΔH}}_1}^{\′}\·\mathrm{tg\α}\′}{\mathrm{tg\β}-\mathrm{tg\α}\′}$

进而可得： ${x_1}^{\′}=\frac{{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

y₂＝(x₂′+ΔL′)·tgα′       (5) 

$\mathrm{tg\β}=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{{x_2}^{\′}+{\mathrm{\ΔL}}^{\′}}---\left(6\right)$

即 ${x_2}^{\′}=\frac{{y_2}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}-{\mathrm{\ΔL}}^{\′}---\left(7\right)$

将式(7)代入式(5)可得： 

$y_2=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(8\right)$

于是可得： $y_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}\·{\mathrm{tg\α}}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

进而可得： $x_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}-{\mathrm{\ΔL}}^{\′}$

于是可得： 

$X^{\′}=\frac{{x_1}^{\′}+{x_2}^{\′}}{2}=\frac{{{\mathrm{\ΔH}}_1}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{2(\mathrm{tg\β}-{\mathrm{tg\α}}^{\′})}-\frac{{\mathrm{\ΔL}}^{\′}}{2}$

X′为剥蚀线的位置。 

将上述的方法应用于胜利油田某区块，效果如下： 

表1超覆位置确定误差表 

表2剥蚀位置确定误差表 

以上数据说明本发明的一种复杂沉积体边界的确定方法在实际应用中得到证实，对剥蚀和超覆识别精度较高。 

可知，生产情况说明，该层水洗程度并不强，驱油效率在33％左右是符合实际情况的，这也证明了校正方法是科学可靠的。 

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。 

Claims (2)

1.一种复杂沉积体边界的确定方法，其特征在于，包括以下步骤：

(1)高精度三维地震剖面追踪：首先对研究区所有钻井资料进行统层对比分析，在单井上识别出超覆和剥蚀的深度点；在完钻井层位标定基础上，通过高精度三维地震剖面同相轴追踪，确定地层超覆和剥蚀线的大体位置；

(2)地震属性分析：通过对平均瞬时相位、平均振幅地震属性资料进行分析研究区储层的接触关系，进而判断地层超覆和剥蚀线的位置范围；

(3)三角函数预测法：在地层构造解释的基础上，获得地层倾角，进而利用三角函数预测法进行地层超覆和剥蚀位置的判断。

2.根据权利要求1所述的一种复杂沉积体边界的确定方法，其特征在于，所述三角函数预测法包括以下步骤：

a.超覆线位置确定

Ng底的倾角为α；1井与2井之间距离为△L，1井与2井Ng1到Ng底的距离分别为△H₁和△H₂；

为了能够降低计算的误差，假设层位超覆位置距离1井与2井分别为x₁和x₂+△L，最终取x₁与x₂的平均值作为距离1井的层位超覆位置；

ΔH₁＝X₁·tgα

即 $x_1=\frac{\mathrm{\Δ}{H}_1}{\mathrm{tg\α}}$

ΔH₂＝(x₂+ΔL)·tgα

即 $x_2=\frac{\mathrm{\Δ}{H}_2}{\mathrm{tg\α}}-\mathrm{\ΔL}$

于是可得：

$x=\frac{x_1+x_2}{2}=\frac{\mathrm{\Δ}{H}_1+\mathrm{\Δ}{H}_2}{2\mathrm{tg\α}}-\frac{\mathrm{\ΔL}}{2}$

X为超覆线的位置；

b.剥蚀线位置确定

Ng底的倾角为α′；ES3上的倾角为β；1井与2井之间距离为△L′，1井与2井Ng底到ES3上顶的距离分别为△H₁′和△H₂′；

为了能够降低计算的误差，假设ES3上的剥蚀位置距离1井与2井分别为x₁′和x₂′+△L′，ES3上的剥蚀位置在1井与2井上的水平投影距离Ng底分别为y₁和y₂，最终取x₁′与x₂′的平均值作为距离1井的ES3上的剥蚀位置；

y₁＝x₁′·tgα′               (1)

$\mathrm{tg\β}=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{{x_1}^{\′}}---\left(2\right)$

即 ${x_1}^{\′}=\frac{y_1+\mathrm{\Δ}{H_1}^{\′}}{\mathrm{tg\β}}---\left(3\right)$

将式(3)代入式(1)可得：

$y_1=\frac{y_1+{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(4\right)$

于是可得： $y_1=\frac{\mathrm{\Δ}{H_1}^{\′}\·{\mathrm{tg\α}}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

进而可得： ${x_1}^{\′}=\frac{{{\mathrm{\ΔH}}_1}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}$

y_２＝(ｘ_２′＋ΔL′)·tgα′       (5)

$\mathrm{tg\β}=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{{x_2}^{\′}+{\mathrm{\ΔL}}^{\′}}---\left(6\right)$

即 ${x_2}^{\′}=\frac{{y_2}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}-{\mathrm{\ΔL}}^{\′}---\left(7\right)$

将式(7)代入式(5)可得：

$y_2=\frac{y_2+{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}}\·{\mathrm{tg\α}}^{\′}---\left(8\right)$

于是可得： $y_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}\·{\mathrm{tg\α}}^{\′}}{\mathrm{tg\β}-\mathrm{tg}{\mathrm{\α}}^{\′}}$

进而可得： $x_2=\frac{{{\mathrm{\ΔH}}_2}^{\′}}{\mathrm{tg\β}-{\mathrm{tg\α}}^{\′}}-{\mathrm{\ΔL}}^{\′}$

于是可得：

$x^{\′}=\frac{{x_1}^{\′}+{x_2}^{\′}}{2}=\frac{{{\mathrm{\ΔH}}_1}^{\′}+{{\mathrm{\ΔH}}_2}^{\′}}{2(\mathrm{tg\β}-{\mathrm{tg\α}}^{\′})}-\frac{{\mathrm{\ΔL}}^{\′}}{2}$

X′为剥蚀线的位置。