专利名称：基于地震纹分析的储层检测方法
技术领域：
本发明涉及石油地球物理勘探领域，更具体地涉及一种基于地震纹分析的储层检测方法。
背景技术：
利用地震数据进行岩性或烃类识别是石油物探研究人员的追求目标，同时也是难点问题。虽然现在已经发展了很多新方法、新技术，如亮点、平点、暗点、多属性融合、AV0、 AVA等方法技术，也取得了成功的应用实例。但这些方法的成功应用都必须具备较多的前提条件，首要是要有高保真的数据(如果能有高保真的横波数据更好)，其次是要求储层界面平缓并有一定的规模。随着勘探地区越来越复杂和勘探精度要求越来越高，对具体的某一勘探地区而言，这些条件不一定都满足，因而无法满足勘探需求。因此，需要一种适合不同地区条件的更简单高效的储层检测技术。
发明内容
针对上述现有技术的问题，本发明提供一种基于地震纹分析的储层检测方法，能够更简单高效地确定储层特征。本发明提供了一种基于地震纹分析的储层检测方法，包括以下步骤标定目标层位；针对目标层位确定基于声学系数的地震纹参数，所述声学系数包括倒谱系数、线性预测倒谱系数、梅尔倒谱系数中的至少一种；通过所述地震纹参数确定所述目标层位的储层特征。优选地，在本发明的各实施例中，所述确定地震纹参数包括通过针对所述目标层位的时窗计算地震参数；根据所述地震参数的一阶或多阶声学系数的值的分布，优选地进一步根据所述声学系数的交会图，从所述地震参数中确定出所述地震纹参数。优选地，在本发明的各实施例中，通过针对所述目标层位的井信息和/或非井信息，从所述地震参数中确定出所述地震纹参数。优选地，在本发明的各实施例中，所述储层特征包括以下中的至少一种岩性，烃类特征，裂缝，断层，河道。优选地，在本发明的各实施例中，所述基于地震纹的储层检测方法进一步包括通过以下中的至少一种信息标定所述目标层位地质信息，井信息，非井信息，实际地震记录信息，合成地震记录信息，层位记录信息。优选地，在本发明的各实施例中，所述基于地震纹的储层检测方法进一步包括对于标定的所述目标层位，确定非井区域的所述地震纹参数。优选地，在本发明的各实施例中，所述确定地震纹参数包括对地震记录信息作傅里叶变换或Z变换到频率域中并在所述频率域计算地震能量，对所述地震能量的绝对值取对数并对所述对数逆傅里叶变换或逆Z变换到时间域，得到不同阶的所述地震纹参数。优选地，在本发明的各实施例中，通过所述地震纹参数确定所述目标层位的储层特征包括通过对不同阶的所述地震纹参数的声学系数的交会图分析，确定所述目标层位中对应于不同地质结构的所述地震纹参数的差异，由此确定不同所述目标层位中的地质结构的性质。优选地，在本发明的各实施例中，所述地震纹参数包括线性预测系数，并采用自相关法计算所述线性预测系数。优选地，在本发明的各实施例中，所述目标层位包括油气矿层或金属矿层或无机非金属矿层。通过本发明提供的基于地震纹的储层检测方法，借助源自地震数据的地震纹参数，能够由此更简单高效地确定储层特征，例如岩性或烃类性质。
为了更清楚地说明本发明实施例或现有技术中的技术方案，以下将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，以下描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图所示实施例得到其它的实施例及其附图。图1为正演地震记录的示意图。图2为利用本发明的实施例的方案得到的1阶倒谱系数。图3为利用本发明的实施例的方案得到的2阶倒谱系数。图4为图2的1阶倒谱系数与图3的2阶倒谱系数的交会图。图5为根据本发明的实施例的基于地震纹的储层检测方法的流程示意图。
具体实施例方式众所周知，不同的储层特征(例如岩性或烃类性质)能够引起地震信号的不同变化，但是从常规的地震信号或现有的地震属性图上很难分辨出这种微小的变化。因此，利用标识性强，容错性好、能够反应地震信号微弱变化的地震特征参数确定储层特征(例如岩性或烃类性质)，正是本发明的关键之处。地震信号与声音信号在本质上是一致的，因而也可基于声波方程进行表征。声纹识别技术的核心是利用标识性强、容错性好的声纹参数，目前已成功进入实用阶段，因此， 可以借鉴声纹识别技术来分析地震信号，从而达到检测/确定储层特征(例如岩性或烃类性质)的目的。在本发明的实施例中，利用地震纹(seismic print)参数作为地震特征参数反映地震信号的微弱变化/差异，以确定导致这种微弱变化/差异的相应的储层特征(例如岩性或烃类性质)。“地震纹”对应于地震记录(例如，实际地震记录或合成地震记录)上的起伏波纹，称能表征所述储层的传播介质的地质或地质地球物理属性的地震纹或地震纹组合可作为地震纹特征，而称能表征地震纹特征的几何学、动力学、信号学参数可作为地震纹参数，由此称能表征所述地震纹参数的设计、获取、筛选、标定，单参数及组合参数的图形、图像化表示等就是地震纹分析。地震纹分析的一般步骤是确定目标层位的地震信号，确定合适的时窗，进行多种地震参数的获取和计算，并从中筛选出标识性强、容错性好的地震纹参数，以参数曲线或T/ D(时间/深度)_参数/数据融合参数图像的形式表示地震纹参数，并对其相应的地质、地球物理含义进行判别。以下将结合附图对本发明各实施例的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例， 本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施例，都属于本发明所保护的范围。本发明提供了一种基于地震纹的储层检测方法，包括以下步骤(参见图5中所示步骤 110、120、130)标定目标层位；针对目标层位确定基于声学系数的地震纹参数，所述声学系数包括倒谱系数、线性预测倒谱系数、梅尔倒谱系数中的至少一种；通过所述地震纹参数确定所述目标层位的储层特征。由此，针对需要勘探的目标层位(其中可包括含矿层、含水层或含气层等)，借助适合的地震数据来获取地震纹参数，进而利用地震纹参数确定目标层位的储层特征，例如岩性或烃类性质，如果目标层位的储层特征符合特定要求，则可确认为适合的储层。在本发明的各实施例中，优选地，所述确定地震纹参数包括通过针对所述目标层位的时窗计算地震参数；根据所述地震参数的一阶或多阶声学系数的值的分布，优选地进一步根据所述声学系数的交会图，从所述地震参数中确定出所述地震纹参数。例如，在如图2和3中所示的1阶和2阶倒谱系数的值的分布中，可看出地震信号的微弱变化/差异被有效放大而变得更明显，有利于分辨出目标层位的地质性质；而且进一步优选地，还可在此基础上构建声学系数交会图，例如在图4中所示的倒谱系数交会图， 从中可更清楚地看出在分布上的差异。优选地，在本发明的各实施例中，通过针对所述目标层位的井信息和/或非井信息，从所述地震参数中确定出所述地震纹参数。通过获知井信息和/或非井信息，可更有利于分辨出目标层位的地质性质。在确定地震纹参数时，可以根据实际需要计算目标层间的地震纹参数，也可以计算整条剖面的地震纹参数剖面，甚至可以计算3D地震数据沿层切片的地震纹参数。由此，可从计算出的多个初始地震纹中筛选出标识性强、容错性好的地震纹参数， 用于确定所述目标层位的储层特征。在本发明的各实施例中，优选地，所述储层特征包括以下中的至少一种岩性，烃类特征，裂缝，断层，河道。这样，通过本发明的实施例，能够区分出地震信号的微弱变化，由此不仅能够用来识别地质状况(例如岩性或烃类性质)，而且也可以检测裂缝、河道等微小地质尺度引起的地震信号变化，从而有助于确定所述目标层位的储层特征。在本发明的各实施例中，优选地，所述基于地震纹的储层检测方法进一步包括通过以下中的至少一种信息标定所述目标层位地质信息，井信息，非井信息，地震记录信息(例如实际地震记录信息或合成地震记录信息)，层位记录信息。在本发明的各实施例中，优选地，所述基于地震纹的储层检测方法进一步包括对于标定的所述目标层位，确定非井区域的所述地震纹参数。通过获知非井区域的所述地震纹参数，可更有利于分辨出目标层位的地质性质。在本发明的各实施例中，优选地，所述确定地震纹参数包括对地震记录信息作傅里叶变换或Z变换到频率域中并在所述频率域计算地震能量，对所述地震能量的绝对值取对数并对所述对数逆傅里叶变换或逆Z变换到时间域，得到不同阶的所述地震纹参数。由此，通过对地震记录信息(例如实际地震信息或合成地震信息)进行计算处理，可得到多个地震信息，然后从这些地震信息中可以筛选出标识性强、容错性好的地震信息作为适于进行储层检测的地震纹参数。在本发明的各实施例中，优选地，所述通过所述地震纹参数确定所述目标层位的储层特征包括通过对不同阶的所述地震纹参数的声学系数的交会图分析，确定所述目标层位中对应于不同地质结构的所述地震纹参数的差异，由此确定不同所述目标层位中的地质结构的性质。例如，可借助1、2阶倒谱系数交会图的分布，确定由不同岩性、流体等引起的地震声纹参数区别，从而进行岩性、烃类检测。在本发明的各实施例中，优选地，所述目标层位包括油气矿层或金属矿层或无机非金属矿层。 在本发明的各实施例中，可采用倒谱系数，线性预测倒谱系数，梅尔倒谱系数中的至少一种作为地震纹参数。具体如下1、倒谱系数将地震信号作为时间序列的信号x(n)而进行Z变换，变换后的值的绝对值取对数，再做逆Z变换，得χ (η)的复倒谱，即χ (η)的复倒谱定义为χ{η) = Z1 [In | Z(x(n)) |]( 1 )显然，一个时间序列的复倒谱仍然是一个时间序列，因此倒谱分析是一种同态分析。在实际应用中，优选地，常用快速傅里叶变换(FFT)代替Z变换。以上得到就是倒谱系数。2、线性预测倒谱系数线性预测倒谱系数在某种程度上就是对线性预测系数作倒谱分析。在地震信号中，相邻采样值之间存在很大的相关性，某时刻的信号在很大程度上可以利用此前的观测样值预测获得，即某一时刻的采样值可以利用此时刻之前数个时刻的采样值的线性组合来逼近。而线性预测倒谱系数兼顾了线性预测系数和倒谱系数的优点，在一定情况下更能够反应出地震信号的微弱差异。其计算的核心问题是求解线性预测系数的系数线性方程组。 常规的求解方法有自相关法、协方差法和格型法等。优选地，采用自相关法，能够保证系统的稳定性，计算精度也较高，并且具有高效的递推算法。采用自相关法计算线性预测系数时，可采用Durbin-To印Iz算法，该算法的计算步骤如下①计算自相关系数Rn(j)，j = 0，1，· · ·，ρ ；(2)E(0) = Rn(O)；③i = 1，开始按如下公式进行递推运算
权利要求
1.一种基于地震纹分析的储层检测方法，其特征在于，包括以下步骤标定目标层位；针对目标层位确定基于声学系数的地震纹参数，所述声学系数包括倒谱系数、线性预测倒谱系数、梅尔倒谱系数中的至少一种；通过所述地震纹参数确定所述目标层位的储层特征。
2.如权利要求1所述的基于地震纹分析的储层检测方法，其特征在于，所述确定地震纹参数包括通过针对所述目标层位的时窗计算地震参数；根据所述地震参数的一阶或多阶声学系数的值的分布，优选地进一步根据所述声学系数的交会图，从所述地震参数中确定出所述地震纹参数。
3.如权利要求1或2所述的基于地震纹分析的储层检测方法，其特征在于，通过针对所述目标层位的井信息和/或非井信息，从所述地震参数中确定出所述地震纹参数。
4.如权利要求1至3中任一项所述的基于地震纹分析的储层检测方法，其特征在于，所述储层特征包括以下中的至少一种岩性，烃类特征，裂缝，断层，河道。
5.如权利要求1至4中任一项所述的基于地震纹分析的储层检测方法，其特征在于，进一步包括通过以下中的至少一种信息标定所述目标层位地质信息，井信息，非井信息，实际地震记录信息，合成地震记录信息，层位记录信息。
6.如权利要求1至5中任一项所述的基于地震纹分析的储层检测方法，其特征在于，进一步包括对于标定的所述目标层位，确定非井区域的所述地震纹参数。
7.如权利要求1至6中任一项所述的基于地震纹分析的储层检测方法，其特征在于，所述确定地震纹参数包括对地震记录信息作傅里叶变换或Z变换到频率域中并在所述频率域计算地震能量，对所述地震能量的绝对值取对数并对所述对数逆傅里叶变换或逆Z变换到时间域，得到不同阶的所述地震纹参数。
8.如权利要求1至7中任一项所述的基于地震纹分析的储层检测方法，其特征在于，通过所述地震纹参数确定所述目标层位的储层特征包括通过对不同阶的所述地震纹参数的声学系数的交会图分析，确定所述目标层位中对应于不同地质结构的所述地震纹参数的差异，由此确定不同所述目标层位中的地质结构的性质。
9.如权利要求1至8中任一项所述的基于地震纹分析的储层检测方法，其特征在于，所述地震纹参数包括线性预测系数，并采用自相关法计算所述线性预测系数。
10.如权利要求1至9中任一项所述的基于地震纹分析的储层检测方法，其特征在于， 所述目标层位包括油气矿层或金属矿层或无机非金属矿层。
全文摘要
本发明涉及一种基于地震纹分析的储层检测方法，包括以下步骤针对目标层位确定基于声学系数的地震纹参数，所述声学系数包括倒谱系数、梅尔倒谱系数、线性预测倒谱系数中的至少一种；通过所述地震纹参数确定所述目标层位的储层特征。通过本发明提供的储层检测方法，借助源自地震数据的地震纹参数，能够由此更简单高效地确定储层特征，例如岩性或烃类性质。
文档编号G01V1/30GK102253414SQ20111016561
公开日2011年11月23日 申请日期2011年6月20日 优先权日2011年6月20日
发明者曹俊兴, 田仁飞 申请人:成都理工大学