CN105093330A - 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法 - Google Patents

Abstract

多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，属于利用地震技术勘探地下储集体技术领域。对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；应用动态波形匹配算法计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差，得到积累振幅差曲线；对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；对这些多分辨率分解曲线作时深转换，与测井资料划分的储层进行分析对比；研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系，进而确定地震横向追踪的属性。

Description

多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法

技术领域

本发明涉及利用地震技术勘探地下储集体技术领域，具体而言，涉及一种多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法。

背景技术

石油地震勘探的基本目标是根据人工震源激发的地震反射波了解地下的构造及相应地层属性。原始的地震资料经过一系列的滤波、叠加、偏移处理后形成相应的自激自收反射波记录数据体。在此基础上，结合区域地质背景在根据数据体绘制的剖面上进行层位追踪、断层解释及储层预测，以寻找地下有利油气储集的构造与圈闭。

在缝洞型储层预测方面，目前主要有两类方法：以单道为基础的波阻抗反演方法和多道相干或波形匹配方法。两者在寻找缝洞储集体方面均存在不足。前者存在的问题是，由于地层中缝洞发育，地层不再是层状的或准层状的，地震波反演的理论假设不再成立，导致应用效果不好。后者多道数据相干或动态波形匹配方法尽管考虑了地层在横向上的变化，但未考虑局部地层在纵向上的变化，地震资料提取属性不能与测井资料划分的储层段对比分析。

发明内容

本发明的主要目的在于提供一种多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，以解决现有技术中的缝洞型储层预测方法中不能同时获取缝洞型储层纵向与横向变化的问题。

为了克服现有技术的不足,本发明提供多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法。

多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法,含有以下步骤；

a.对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；

b.以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；

c.应用动态波形匹配算法计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差，得到积累振幅差曲线；

d.对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；

e.对这些多分辨率分解曲线作时深转换，与测井资料划分的储层进行分析对比；

f.研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系，进而确定地震横向追踪的属性；

g.根据确定的地震属性绘制连井剖面，研究井中储层所对应的横向分布；

h.根据不同井的产量特征及生产情况，推广应用到整个地震资料覆盖区，估算井旁缝洞储集体视体积及储量。

作为优选，所述叠加偏移处理是叠前深度偏移，或者是叠后时间或深度偏移。

作为优选，所述地震道波形匹配积累振幅差过程如下，设目的层段窗长为N个采样点，初始计算时从n＝1开始进行波形匹配计算振幅差，接着计算n＝2的波形匹配积累振幅差，一直计算到n＝N，最终得到一条波形匹配积累振幅差曲线f(n)(n＝1，2，…，N)。

作为优选，所述估算缝洞视体积的步骤如下：

(1)每一CDP点上的视厚度

$H_i=\underset{j}{\Σ}{H}_{ij}$

式中H_i是第i道的总厚度，H_ij是第i道j段的厚度。

(2)在地震属性数据体上对每一条线、每一CDP点变厚度求视体积V

$V=\underset{i}{\Σ}{S}_i{H}_i$

式中S_i是CDPi面元上的面积，H_i是第i道的总厚度。

一种多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，该方法包括：步骤a，.对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；步骤b，以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；步骤c，获取积累振幅差曲线；步骤d.，对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；步骤e，对上述多条曲线作时深转换，与测井资料划分的储层进行分析对比；步骤f，确定地震横向追踪的属性；步骤g，根据确定的地震属性绘制连井剖面，研究井中储层所对应的横向分布；步骤h，根据不同井的产量特征及生产情况，推广应用到整个地震资料覆盖区，估算井旁缝洞储集体视体积及储量。

进一步地，上述叠加偏移处理是叠前深度偏移，或者是叠后时间或深度偏移。

进一步地，上述步骤c包括：步骤c1，计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差；步骤c2，根据上述地震道波形匹配积累振幅差获取上述积累振幅差曲线。

进一步地，采用动态波形匹配算法计算上述步骤c1中的上述地震道波形匹配积累振幅差。

进一步地，上述地震道波形匹配积累振幅差过程如下，设目的层段窗长为N个采样点，初始计算时从n＝1开始进行波形匹配计算振幅差，接着计算n＝2的波形匹配积累振幅差，一直计算到n＝N，最终得到一条波形匹配积累振幅差曲线f(n)(n＝1，2，…，N)。

进一步地，采用二进小波分解对上述动态波形匹配积累振幅差进行上述多尺度分解。

进一步地，上述步骤f包括：步骤f1，研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系；步骤f2，根据上述步骤f1中的上述对应关系确定上述地震横向追踪的属性。

进一步地，上述估算缝洞视体积的步骤如下：步骤S1，采用公式计算每一CDP点上的视厚度，式中H_i是第i道的总厚度，H_ij是第i道j段的厚度；步骤S2，利用公式在地震属性数据体上对每一条线、每一CDP点变厚度求视体积V，式中S_i是CDPi面元上的面积，H_i是第i道的总厚度。

本发明针对已有方法的局限性，提出了一种同时表达缝洞型储层纵向与横向变化的方法。在面元地震道上计算动态波形匹配振幅差表达地层横向变化，逐步增加波形匹配的点数计算积累振幅差表达局部地层的纵向变化，进而把这种纵向变化通过二进小波分解谱表达出来，在已钻井的井点处与测井资料划分的储层段对比分析，能准确预测缝洞储集体在空间的分布。为石油勘探、开发中的评价井提供井位坐标、侧钻井侧钻方位、碳酸盐岩缝洞储层分布，也为碳酸盐岩缝洞型储层储量计算提供视体积等及相关图件。

附图说明

构成本申请的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1示出了根据本发明的本申请一种典型实施方式提供的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法的流程图；

图2示出了动态波形匹配振幅差表达缝洞纵、横向变化示意图，有为不同共深度点(CDP)地震道波形穿过缝洞体示意图、不同斜率的动态波形匹配积累振幅差示意图；

图3示出了TKXXX井测井资料解释的储层段与井点处地震资料波形匹配积累震幅差小波分解分量的对比图，其中有TKXXX井测井资料储层评价图、TKXXX井井点处地震资料波形匹配积累震幅差的多分辨率分解分量图；

图4a示出了TKXXX井东西方向积累震幅差一阶小波分量过井剖面图(东西方向)；

图4b为TKXXX井东西方向积累震幅差一阶小波分量过井剖面图(南北方向)；

图5示出了另一地区HaXXX-X井测井资料解释的储层段与井点处地震资料波形匹配积累震幅差小波分解分量的对比图，其中有测井资料储层评价结果、地震资料提取的积累震幅差的小波分解分量；

图6a示出了HAXXX-X井井点积累振幅差一阶小波分量剖面与均方根差平面分布图；

图6b示出了HAXXX-X井井点处积累振幅差一阶分量剖面与测井储层对比图。

具体实施方式

需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。

显然，本领域技术人员基于本发明的宗旨所做的许多修改和变化属于本发明的保护范围。

实施例1：如图1、图2、图3、图4a、图4b、图5、图6a、图6b所示，

多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，包括如下步骤：

a.对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；

原始地震资料的基本处理与解释采用已有的方法。其中偏移处理可以是叠前深度偏移或叠后时间或深度偏移。层位拾取在解释工作站上进行。这样就得到了叠加偏移后的三维数据体及解释的层位数据。

b.以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；

根据研究目的层的层位数据决定地震道数据的的起点，同时根据目的层的最大厚度及地震数据体的采样间隔定出波形匹配的采样点数N。两条相邻地震道之间的波形差异大小问题转化为动态规划中最优路径选择问题。

c.应用动态波形匹配算法计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差，得到积累振幅差曲线；

地震道波形匹配振幅差用来度量地层的横向变化的大小。若地层在横向上没有大的变化，则相邻地震道地震波形之间的差异不大，从而波形匹配振幅差就低。反之，若地层在横向上有较大的变化，如有断层、裂缝发育带或大溶洞附近等，则相邻地震道地震波形之间的差异变大，从而波形匹配振幅差就高。因此，动态波形匹配振幅差可以作为地层横向变化的一种度量。两者的关系为，地层横向变化小，则波形匹配振幅差低，地层横向变化大，则波形匹配振幅差高，如图3所示。由于缝洞储集体在纵向上(深度上)也是变化的，为了表达这种纵向上的变化，提出计算积累振幅差的方法来表达地震道波形振幅在深度上的差异。

设目的层段窗长为N个采样点。初始计算时从n＝1开始进行波形匹配计算振幅差，接着计算n＝2的波形匹配积累振幅差，一直计算到n＝N，这样就可以得到一条波形匹配积累振幅差曲线f(n)(n＝1，2，…，N)。如图2所示，表达了动态波形匹配振幅差表达缝洞纵、横向变化示意图；不同共深度点(CDP)地震道波形穿过缝洞体示意图；不同斜率的动态波形匹配积累振幅差示意图。

波形匹配积累振幅差曲线是随着采样点的增加而增加的增函数。在不同的深度段增加的快慢不同：若在纵向上(某个深度段)地层的横向变化激烈，则积累振幅差异变大，积累振幅差曲线的斜率变陡；若在某个深度段地层的横向变化平缓，则积累振幅差曲线斜率较小。由此可知，积累振幅差曲线同时表达了两个方面的信息：一方面是表达了地震道波形的横向变化大小；另一方面也表达了振幅差随深度变化的快慢，亦即表达了地层纵向上的非均质性。对于台地碳酸盐岩地层，地层横向上与纵向上的非均质性与缝洞储集体的发育相关。

d.对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；

地震道波形动态匹配积累振幅差是随采样点增加而增加的函数，其增加的快慢与缝洞体在纵向、横向上的发育相关。为了研究积累振幅差中的变化信息与井中划分的储层段的联系。对积累振幅差曲线进行二进小波分解，在不同小波分解分量上研究井中根据测井资料划分的储层段与不同阶小波分解谱的关系。

积累振幅差随深度(或时间)变化的函数f(x)在位置x对尺度为s的小波变换为

W_sf(x)＝f(x)*Ψ_s(x)

取尺度为s＝{2^j}_j∈Z，定义则函数f(x)对尺度为2j的小波变换[5]为

$W_{2^j}f\left(x\right)=f*{\mathrm{\Ψ}}_{2^j}\left(x\right)$

在频率域的形式为

${\hat{W}}_{2^j}f\left(w\right)=\hat{f}\left(w\right)\mathrm{\Ψ}\left(2^{j}w\right)$

设θ(x)为一平滑函数，令Ψ(x)为θ(x)的一阶导数，

$\mathrm{\Ψ}\left(x\right)=\frac{d\mathrm{\θ}\left(x\right)}{dx}$

记则对尺度为2j的小波变换为

$W_{2^j}f\left(x\right)=f\left(x\right)*{\mathrm{\Ψ}}_{2^j}\left(x\right)=f*\left(2^j\frac{{\mathrm{d\θ}}_{2^j}\left(x\right)}{dx}\right)=2^j\frac{d}{dx}(f*{\mathrm{\θ}}_{2^j})\left(x\right)$

小波变换的极大值提供了一种多尺度的检测函数奇性点位置的方法，且信号变化幅度越大，则相应小波变换的幅度值越大。对于实际问题，信号的可测分辨率是有限的。分辨率通常取有限值，即把j值限定在0与J值之间。20表示最高分辨率，2J表示最低分辨率。于是可得到信号小波变换的多分辨率表示。这种多分辨率表示对于检测信号的不同级次的变化来说有很好的性质。当j值较小时，用对函数f(x)光滑化的结果对f(x)大的突变部分的形态影响不大；而当j值较大时，则此光滑将会将f(x)的细小部分(如信号中的高频成分)消去而剩下尺度较大的突变信号(信号中大的变化部分)。因此，具体应用中可根据研究问题的需要选择某些j值而检测出某个级次的变化。在实现这一算法时采用具有一阶消失矩的小波函数ψ(x)，它是三次B-样条函数的一阶导数，其Fourier变换为

本发明中对地震道波形匹配积累振幅差选择分辨率高的j＝1表达地层纵、横向细微的变化。在积累振幅差曲线斜率变大的深度，其小波变换出现上升；在积累振幅差曲线斜率变小的深度，其小波变换出现下降。这样在缝洞发育的部位，积累振幅差曲线的小波变换出现一个单峰。参见图2。

e.对这些多分辨率分解曲线作时深转换，与测井资料划分的储层进行分析对比；

如前上述，在缝洞发育的部位，地震积累振幅差曲线的小波变换出现一个单峰。单峰的半宽度与储层的厚度相关。如果在此处钻探了一口井，测量了相应的测井资料，对测井资料进行处理、解释，则可以根据处理、解释结果划分储层段(厚度)。测井资料解释的储层段与井点处地震资料积累振幅差曲线的小波变换曲线对比，研究两者之间的对应关系。测井资料解释的大套缝洞储层段与地震积累振幅差曲线小波变换谱上的单峰相对应。通过井震对比追踪井旁缝洞储集体的发育范围，认识单井产量与井旁储集体大小的关系。

f.研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系，进而确定地震横向追踪的属性；

g.根据确定的地震属性绘制连井剖面，研究井中储层所对应的横向分布；

h.根据不同井的产量特征及生产情况，推广应用到整个地震资料覆盖区，估算井旁缝洞储集体视体积及储量。

由于缝洞储层纵、横向上的非均质性，平面上每一CDP点的视厚度Hi是不同的。估算缝洞视体积的的步骤如下：

(1)每一CDP点上的视厚度

$H_i=\underset{j}{\Σ}{H}_{ij}$

式中H_i是第i道的总厚度，H_ij是第i道j段的厚度。

(2)在地震属性数据体上对每一条线、每一CDP点变厚度求视体积V

$V=\underset{i}{\Σ}{S}_i{H}_i$

式中S_i是CDPi面元上的面积，H_i是第i道的总厚度。

下面为本发明方法的实际应用说明。

图3是TKXXX井测井资料解释的储层段与井点处地震资料波形匹配积累震幅差小波分解分量的对比图。图3中有TKXXX井测井资料储层评价图。图中，第一道为井径、自然伽马和去铀伽马曲线；第二道为深、浅侧向电阻率曲线；第三道为三孔隙度曲线，即密度测井曲线、中子测井曲线、声波测井曲线；第四道为谱曲线，即U、K、TH；第五道为深度索引；第六道为计算的孔隙成分曲线；第七道为相对连通孔隙度曲线，表达储层的有效性；第八道计算的孔隙结构指数。测井资料上显示有三大套储层，第一套储层为5554-5587米井段，厚度为33米；第二套储层为5604-5622米井段，厚度为18米；第三套储层为5650-5710米井段，厚度为60米。

图3中有TKXXX井井点处地震资料波形匹配积累震幅差的多分辨率分解分量图。图中的纵坐标为目的层顶点为起点的相对深度，图中的横坐标为小波分解分量的幅度。红色实线为该井点处地震资料动态波形匹配费用小波分解一阶分量，蓝色实线为二阶分量，鲜绿色实线为三阶分量，灰色实线为四阶分量变换结果。由地震资料积累费用小波变换原理可知，一阶小波分量大于0的部分是地震资料横向变化大的储层发育段(缝洞发育)。由图可见，地震资料提取的积累震幅差的小波分解属性也显示为三套储层，第一套层为风化壳下66-93米井段，厚度为27米；第二套层为为风化壳下126-147米井段，厚度为21米；第三套层为为风化壳下162-231米井段，厚度为69米。对比图也表明，井点处地震资料积累震幅差的一阶小波分量划分的厚度与测井资料划分的储层段有好的对应关系。

井旁储集体的发育范围影响着井中储层油气产量，地震波形匹配积累振幅差一阶小波分量可以追踪储集体横向上的发育范围。而井周储集体的发育范围与井中累计产量的有联系。图4a为TKXXX井东西方向积累震幅差一阶小波分量过井剖面图(东西方向)；图4b为TKXXX井东西方向积累震幅差一阶小波分量过井剖面图(南北方向)。由图可见，动态波形匹配积累振幅差方法处理的TKXXX井东西方向、南北方向储集体均较发育。5627-5687m井段,射孔酸压,初产，日产液171吨,日产油171吨,含水0％。五年，积累产油19.04万吨,累计产水万0.62万吨。5534.33-5576m井段,射孔酸压,初产，日产液97.4吨,日产油97.4吨,含水0％。一年，积累产油1.37万吨,累计产水1.14万吨。可见，井中储层段发育，根据地震资料累振幅差一阶小波分量表达的井旁对应储集体发育范围大，则油井的累计产量大。可见本发明对井周缝洞储层的横向追踪效果与井的实际产量相符，效果显著。

图5是另一地区HaXXX-X井测井资料解释的储层段与井点处地震资料波形匹配积累震幅差小波分解分量的对比图。有测井资料储层评价结果图。测井资料上显示有两套储层，第一套储层为6621-6650米井段，厚度为29米；第二套储层为6660-6685米井段，厚度为25米。有地震资料提取的积累震幅差的小波分解分量图。图中的纵坐标为地震资料深度偏移后的相对深度，图中的横坐标为小波分解分量的幅度。图中地震资料提取的积累震幅差的小波分解属性显示为两套储层。第一套6871-6903米，厚度约为32米；第二套6980-7008米，厚度约为28米。可见，井点处地震资料积累震幅差的一阶小波分量划分的厚度与测井资料划分的储层段有好的对应关系。

根据累振幅差一阶小波分量数据体估算井旁缝洞储集体的储量，如图6a和图6b所示，图6a是HAXXX-X井井点积累振幅差一阶小波分量剖面与均方根差平面分布图，图6b是HAXXX-X井井点处积累振幅差一阶分量剖面与测井储层对比图。计算地震属性体视体积的公式为

$V_1=\underset{i}{\Σ}{S}_i{H}_i---\left(1\right)$

式中，i为缝洞体平面属性面积求和.m2。由式(1)可以计算出储集空间体积。

V₂＝V₁×Φ(2)

式中，V1为地震属性体视体积，m3；Φ为储层平均孔隙度，小数。根据式(1)、式(2)即可求出油藏储量

N_o＝V₂×S_o×ρ_o/B_oi(3)

式中，V2为储集空间体积，m3；So为含油饱和度，小数；ρo为地面原油密度，t/m3；Boi为原油体积系数。

下面是HAXXX-X井区储量计算实例。

(I)地震属性视体积

HAXXX-X井缝洞体平面属性面积为0.100096×106m2，缝洞体纵向属性厚度约为300m，测井解释的储层段厚度为28.5m。根据(1)式HAXXX-X井区地震属性视体积为

V1＝523.1×104m3

(II)储集空间体积

HAXXX-X井储层平均孔隙度为0.021，储集空间体积为：

V2＝523.1×104×0.021＝11.1×104m³

(III)HAXXX-X井区油藏储量

HAXXX-X井储层含油饱和度取0.85，地面原油密度取0.886，原油体积系数取1.1478。则HAXXX-X井区的油藏储量为：

N0＝11.1×104×0.85×0.886/1.1478＝7.29×104t

该井累计产液量4.44×104t，累计时间587天(统计截止时间：2012-04-09)，产液速度75.66t/天。其中累产油3.914×104t；累产水0.080×104t；累产气0.0447×108m3。可见，该井的累计产液量与HAXXX-X井区油藏储量一致。

从以上的描述中，可以看出，本发明上述的实施例实现了如下技术效果：

本发明针对已有方法的局限性，提出了一种同时表达缝洞型储层纵向与横向变化的方法。在面元地震道上计算动态波形匹配振幅差表达地层横向变化，逐步增加波形匹配的点数计算积累振幅差表达局部地层的纵向变化，进而把这种纵向变化通过二进小波分解谱表达出来，在已钻井的井点处与测井资料划分的储层段对比分析，能准确预测缝洞储集体在空间的分布。为石油勘探、开发中的评价井提供井位坐标、侧钻井侧钻方位、碳酸盐岩缝洞储层分布，也为碳酸盐岩缝洞型储层储量计算提供视体积等及相关图件。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于包括如下步骤：

a.对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；

b.以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；

c.应用动态波形匹配算法计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差，得到积累振幅差曲线；

d.对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；

e.对这些多分辨率分解曲线作时深转换，与测井资料划分的储层进行分析对比；

f.研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系，进而确定地震横向追踪的属性；

g.根据确定的地震属性绘制连井剖面，研究井中储层所对应的横向分布；

h.根据不同井的产量特征及生产情况，推广应用到整个地震资料覆盖区，估算井旁缝洞储集体视体积及储量。

2.根据权利要求1所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于所述叠加偏移处理是叠前深度偏移，或者是叠后时间或深度偏移。

3.根据权利要求1所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于所述地震道波形匹配积累振幅差过程如下，设目的层段窗长为N个采样点，初始计算时从n＝1开始进行波形匹配计算振幅差，接着计算n＝2的波形匹配积累振幅差，一直计算到n＝N，最终得到一条波形匹配积累振幅差曲线f(n)(n＝1，2，…，N)。

4.根据权利要求1所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于所述估算缝洞视体积的步骤如下：

(1)每一CDP点上的视厚度：

$H_i=\underset{j}{\mathrm{\Σ}}{H}_{\mathrm{ij}}$

式中H_i是第i道的总厚度，H_ij是第i道j段的厚度。

(2)在地震属性数据体上对每一条线、每一CDP点变厚度求视体积V

$V=\underset{i}{\mathrm{\Σ}}{S}_i{H}_i$

式中S_i是CDPi面元上的面积，H_i是第i道的总厚度。

5.一种多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述方法包括如下步骤：

步骤a，.对原始地震资料进行资料处理和资料解释，得到叠加偏移处理后的三维数据体及解释拾取的层位数据；

步骤b，以地震解释拾取的层位数据为起点，得到井旁相邻地震道目的层数据；

步骤c，获取积累振幅差曲线；

步骤d，对动态波形匹配积累振幅差进行多尺度分解，将一条动态积累振幅差曲线变为多条曲线；

步骤e，对所述多条曲线作时深转换，与测井资料划分的储层进行分析对比；

步骤f，确定地震横向追踪的属性；

步骤g，根据确定的地震属性绘制连井剖面，研究井中储层所对应的横向分布；

步骤h，根据不同井的产量特征及生产情况，推广应用到整个地震资料覆盖区，估算井旁缝洞储集体视体积及储量。

6.根据权利要求5所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述叠加偏移处理是叠前深度偏移，或者是叠后时间或深度偏移。

7.根据权利要求5所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述步骤c包括：

步骤c1，计算以地震解释拾取的层位数据起点至终点面元相邻地震道波形匹配积累振幅差；

步骤c2，根据所述地震道波形匹配积累振幅差获取所述积累振幅差曲线。

8.根据权利要求7所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，采用动态波形匹配算法计算所述步骤c1中的所述地震道波形匹配积累振幅差。

9.根据权利要求8所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述地震道波形匹配积累振幅差过程如下，设目的层段窗长为N个采样点，初始计算时从n＝1开始进行波形匹配计算振幅差，接着计算n＝2的波形匹配积累振幅差，一直计算到n＝N，最终得到一条波形匹配积累振幅差曲线f(n)(n＝1，2，…，N)。

10.根据权利要求5所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，采用二进小波分解对所述动态波形匹配积累振幅差进行所述多尺度分解。

11.根据权利要求5所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述步骤f包括：

步骤f1，研究考察不同级次分辨率的地震属性与测井资料储层评价划分的厚度对应关系；

步骤f2，根据所述步骤f1中的所述对应关系确定所述地震横向追踪的属性。

12.根据权利要求5所述的多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法，其特征在于，所述估算缝洞视体积的步骤如下：

步骤S1，采用公式计算每一CDP点上的视厚度，式中H_i是第i道的总厚度，H_ij是第i道j段的厚度；

步骤S2，利用公式在地震属性数据体上对每一条线、每一CDP点变厚度求视体积V，式中S_i是CDPi面元上的面积，H_i是第i道的总厚度。