用波形分类预测储层在先前的练习,我们预测了地震属性对油气的影响。
在下一节中,我们将比较属性分析的结果和用波形分类地震子波的特征更好的理解储层的特征。
我们将在相同的测区用top_chalk_(output)层引导运行我们的分析。
波形分类的输出被首先显示在平面图,然后我们探测结果,更多的细节在下章讨论RAVE时探索。
介绍波形分类波形分类是地震波形分类分解使你自动的分组和绘制相似的波形。
从任意选择的地震数据中决定烃源岩特性是困难的。
在过去,人们可以通过测井曲线响应在地震剖面上的垂直时间确定烃源岩的位置,然后绘制通过测网地层。
因为在垂向上由深度到时间的测井数据反演是不精确的,所以有关烃源岩的地震响应导致错误的解释。
一个解释员怎样用远离井的地震道去推断岩性的变化?靠井旁地震道上下的烃源岩,你可以引导计算机去追踪。
相似波组在平面图上显示子波的分类,在勘探区域识别岩性的变化是有可能的。
追踪子波为了指导计算机追踪子波,你或者拾取一个恒定时窗或者地层,设置一个包含目标层时域的时窗长度。
波形分类假设所有道的源子波是相似的,两个邻道由同一地层产生的子波是相似的。
如果一个子波随距离的变化是可见的,那么一个地震同相轴已产生了。
可视化探测子波数据是很容易的,但是数字化探测子波变化是不容易的。
Manhattan距离Manhattan距离或单元块是一种求两点之间的相似度和非相似度的统计方法。
这个对比认为从一个子波到另一个子波只可能沿着子波的相似度移动。
对角线移动是不可能的。
取两个等长的拥有N个采样点的子波作为所有N个采样相对应的采样的差的绝对值相加得到Manhattan距离。
Manhattan距离被给出为:M=∣A1-B1∣+∣A2-B2∣+…+∣A11-B11∣这儿M是Manhattan距离A是与井曲线闭合的参考子波B是层引导的目标子波N是每一个子波的采样点数,上述采样点数为11两个完全相同的子波将使Manhattan距离值为0,而那些完全不同的子波将导致一个大的Manhattan距离值。
已知这个信息,你可以在层引导的目标子波位置上绘制Manhattan平面图。
在灰度显示中相似度可以被清晰的看到:白区表示M=0,黑区表示M为最大的计算值。
白区域体现层引导的目标子波与A很相似而黑区域体现子波与子波A极不相似。
现在这种情况是多口井存在。
Manhattan距离方法必须被修改,从而可进行子波分组会分类。
为了得到或好或坏的近似值每一目标子波依靠参考子波。
应用Manhattan距离比较每一参考子波与目标子波,然后将Manhattan测量值输出,Manhattan距离是子波之间的采样差异值的总和。
这个总值被分给两个子波的各自采样值。
因此,Manhattan距离形式如下:M=总值/标准测量分母总值=∣A1-B1∣+∣A2-B2∣+…+∣A11-B11∣分母=(︱A1∣+∣A2∣)+…+(︱A11︱+︱B1︱+︱B2︱+︱B11︱…)两个完全相同的子波将导致Manhattan距离为零。
不完全相同的子波将产生正的小于或等于Manhattan距离。
因此大部分的不连续数据,例如断层,是非常接近+1的。
如果参考子波和目标子波的M值是0.25或更小,那么那些子波被分在同一个组或类。
作为一组的最大M值门限依靠数据的质量(信噪比),如果数据包括很多的随机噪声,则该组的门限值应该加+0.3。
作为一个没有噪音的模型数据体,则该组子波的门限值应减小0.05。
拥有这些信息,你就可以用参考子波确定分类,然后将周围的目标子波分到不同的参考子波组。
如果一个目标子波属于更多的参考子波组,将目标子波分配给M值最小的一组。
因为组可以重叠,所以确定一个较小的门限值并重新运行子波波形分类。
把与参考子波没有高程度近似的子波放在一个特别的组。
注意:在非指导标准K方法选择算法中不用门限值去确定子波的分组。
运行选择为识别相似子波和给相似子波分类波形分类提供下面的选择。
监督选项给你提供目标区域或已知的生产参考子波的波形分类,修改和绘制相似子波的位置。
非监督选项用第一决定参考子波自动将子波分类,然后识别相似子波给相似子波分类。
非监督选项运行K方法用门限K方法或标准K方法算法重新定义子波分组。
非监督选项比监督选项用大量的算法。
在运行非监督选项之前用监督选项帮助参数估计。
监督子波追踪选项要求确定参考子波的位置。
每一个参考位置由纵测线和横测线上的数字来定义。
由子波数字来识别,用常时地层或拾取地震地层及时的定位。
参考子波应该在数据方面与井位和目标体区域相对应。
一个参考子波数据与发现的任何目标子波时相似的。
方法学参考子波被输入需要依靠纵测线、非纵测线和子波标号,输入格式是:#非纵测线#纵测线:子波井#一旦你已选项了参考子波的位置,详细说明一个分析窗。
你有选项提供一个常时窗或预先拾取地层的选项权。
你将为常时窗选项开始和结束的时间输入一个参考的时间和输入一个时间窗。
注意:为了计算一个Manhattan距离值,子波分类要求每一道的时间值。
任何一个地震地层的之间将产生输出子波地层。
你由权选项或创建一个输出子波类型地层。
为了创建一个新的输出波形类型地层,在输出波形类型地层区给一个新地层定一个名,敲击列表选项存在的地层。
波形分类提供了两种输出格式:子波个数和近似值。
子波个数选项将输出一个为每一个参考子波加拒绝类的子波数字。
色标现实范围应该时从0到1,1表示非常的相似,0则表示完全的不相似。
地震图件的色标显示范围应该被适当的设置。
门限值是M的等同物,变化范围从0到1。
用门限值去限制进入一个特定子波分类的地震道数字。
门限值越小,为了质量分类的目标子波越相似。
一般情况,用作分组的最大门限值M依靠数据的质量。
如果数据中包括许多的随机噪音,这组子波的门限值应该增大0.3。
对于一个无噪音的理想数据体,门限值应该减小0.05。
一个值为0.2是断层。
在经常监督追踪的地层中扩展长度对不一致做补偿是非常有用的。
输入一个目标子波扩展在开始采样之上,结束采样之下的采样数字。
因为过程时间要求波组随扩展长度呈线性上升。
所以不要用过大的扩展值。
参数当你选项程序→波形分类→监督,然后点参数,将出现一个相应对话框分析窗点开常时窗创建波形分类地层,要求该地层是通过地震数据体沿常时切片的。
点参数输出参考时间和中心时窗宽度。
点开层时窗,沿一个存在的地震层创建一个波形分类地层。
当你点出层时窗时,将出现下面的层时窗的对话框。
点击列表在层选项框中选一个层。
靠近断层,时间门限时是中心对称的。
用一个记号或定一个值去设置时窗宽度。
为特定一个层的上下时间代替中心时窗门限,将中心门限关闭,用特定的标记或特定的值确定层的上下时间,点OK应用你的选项。
输出的波形分类层输入一个输出波形分类的层名,也可在层列表中选一个输出波形分类层。
你可以重写一个层,依靠点击清除存在的输出层。
门限用门限值限制一个组的子波数。
门限值越小,在同一组的目标子波与参数子波越相似。
扩展长度(采样)在监督波形分类过程中,扩展长度起到理想的作用。
用这个值补偿在用于引导追踪目的层的不一致性。
练习:进行波形分类在下图中,曲线代表参考子波的位置,这个位置的选项是基于这个区域的数据质量,非常低的信噪比,并且因为我们想绘制在顶部地层结构与亮点相一致的近似值。
1、返回到叠后剖面估计属性的复制流程。
2、证实输出文件被设置为demo16(3dv)区域扩展被定为整个勘探区。
3、波形分类输出地层地震图件,关闭所有的输出文件选项。
4、在过程领域,用波形分类→监督代替属性提取过程,显示参数对话框。
5、定义参考子波位置。
(非纵测线:纵测线:子波数字/)。
该例将用非纵测线为120,纵测线号为195。
如果时间允许,你可以选项额外的或交互的参考位置,看输出地层之间怎样相互影响。
6、打开地层时间窗选项并打开分析窗参数和对话框。
7、选项top_chalk_(output)作为地层名,设置时间宽度为10。
8、在输出地层领域,输入一个名字作为输出地层。
9、选项近似值作为输出的类型,为扩展长度用断层值。
10、当你满足以上设置时,运行工作。
从波型分类中显示输出在运行监督波形分类中,我们选项输出近似值。
波形分类将接受数值,0根本无相似;1:理想的相似,对于数据体中的每一道,为了看到输出的地层,在地震视图中色标的显示范围对话框必须与该范围想匹配。
1、当叠后工作已经完成打开地层平面视图并用波形分类创建的显示地层。
2、设置色标适应显示范围。
3、打开一个新平面视图在相反的显示器,并显示Avg_Inst _phase地层,设置色标显示范围与创建属性地层的值相互匹配。
比较波形分类和平均瞬时相位。
从比较两层显示中可以发现,在波形分类的显示中有另外的一些有用信息,而在属性显示或地震显示中是不存在的。
三个广泛区域指示与参数子波位置非常相似的区域。
并且与当前的天然气生产区域相一致。
非监督(门限K方法)非监督方法在没有井控时被应用。
因为它将会为了相似子波花费很多时间去测算每一条测线和每一道,在计算相似子波之前波形分类将自动把参考子波分组。
方法学非监督(门限K方法)波形分类追踪并将子波分为了3个部分:●自动的子波分组●K方法分组调解●重新分组(合并相似组)自动子波分组自动子波分组开始浏览输入的第一个完整子波。
进程浏览下一个子波,然后用Manahattan距离去决定该子波是否与第一子波相似。
如果M时在门限之内,子波被认为时第一组。
如果M不是在最小的门限之内,子波被认为是一个新组,让每一个参考子波与新的子波做比较来重复该过程,或者认为子波是目前的组成或者认为子波是属于一个新组。
依上可知,根据勘探区域的大小,将来有成百上千的组。
某些组可能只有一个成员,而其它的将会有许多的成员。
用每组子波数量的最小值参数来限制自动子波分组产生的组个数。
任何一个组子波个数比最小组子波个数少的被认为是0组。
该组被放在输出范围结束的最低位,依靠色标的选择。
组根据相似程序有序的排列。
当重新排序时,相似组被挨着按照数字的大小排列。
K方法使用选项性在自动子波分组中K方法调解定义地震道应用输出分组。
因为有许多组,所以来自自动子波分组的输出看起来自然,当在二维平面图上显示时,从自动子波分组的输出可以看到进程。
除非输入的数据体是由无噪音的模型数据体。
K方法调解的第一步是使每一组的子波数相同,这需要依靠增大所有子波的M值,并且分得的结果由成员的个数定。
一旦每一组子波的平均值被计算出来,所有的原始子波与子波的平均值做比较,从而适应最合适的门限。
如果原始子波没有与任何一个子波组相适应,则该组被认为是0组。
通过K方法的不断迭代,比较原始子波与参考子波的过程被定义为1。
用K 方法迭代次数参数在控制这个过程。
通常作为平均数2次迭代是充足的。
重新分组的选项依靠以上参数的输入,层间子波相似的差异,数据的信噪比比,在K方法调解以后,可以有或多或少的组。