波形分类的一些介绍.doc

用波形分类预测储层 在先前的练习，我们预测了地震属性对油气的影响。在下一节中，我们将比较属性分析的结果和用波形分类地震子波的特征更好的理解储层的特征。我们将在相同的测区用top_chalk_(output)层引导运行我们的分析。波形分类的输出被首先显示在平面图，然后我们探测结果，更多的细节在下章讨论RAVE时探索。介绍波形分类波形分类是地震波形分类分解使你自动的分组和绘制相似的波形。从任意选择的地震数据中决定烃源岩特性是困难的。在过去，人们可以通过测井曲线响应在地震剖面上的垂直时间确定烃源岩的位置，然后绘制通过测网地层。因为在垂向上由深度到时间的测井数据反演是不精确的，所以有关烃源岩的地震响应导致错误的解释。一个解释员怎样用远离井的地震道去推断岩性的变化？靠井旁地震道上下的烃源岩，你可以引导计算机去追踪。相似波组在平面图上显示子波的分类，在勘探区域识别岩性的变化是有可能的。追踪子波为了指导计算机追踪子波，你或者拾取一个恒定时窗或者地层，设置一个包含目标层时域的时窗长度。波形分类假设所有道的源子波是相似的，两个邻道由同一地层产生的子波是相似的。如果一个子波随距离的变化是可见的，那么一个地震同相轴已产生了。可视化探测子波数据是很容易的，但是数字化探测子波变化是不容易的。Manhattan距离Manhattan距离或单元块是一种求两点之间的相似度和非相似度的统计方法。这个对比认为从一个子波到另一个子波只可能沿着子波的相似度移动。对角线移动是不可能的。取两个等长的拥有N个采样点的子波作为所有N个采样相对应的采样的差的绝对值相加得到Manhattan距离。Manhattan距离被给出为：MA1-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应用你的选项。输出的波形分类层 输入一个输出波形分类的层名，也可在层列表中选一个输出波形分类层。你可以重写一个层，依靠点击清除存在的输出层。门限 用门限值限制一个组的子波数。门限值越小，在同一组的目标子波与参数子波越相似。门限的有效参数范围MinimumMaximumDefault0.011.00.4扩展长度（采样）在监督波形分类过程中，扩展长度起到理想的作用。用这个值补偿在用于引导追踪目的层的不一致性。扩展长度的有效参数范围MinimumMaximum Default 0 100 0练习：进行波形分类在下图中，曲线代表参考子波的位置，这个位置的选项是基于这个区域的数据质量，非常低的信噪比，并且因为我们想绘制在顶部地层结构与亮点相一致的近似值。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个部分：l 自动的子波分组l K方法分组调解l 重新分组(合并相似组)自动子波分组自动子波分组开始浏览输入的第一个完整子波。进程浏览下一个子波，然后用Manahattan距离去决定该子波是否与第一子波相似。如果M时在门限之内，子波被认为时第一组。如果M不是在最小的门限之内，子波被认为是一个新组，让每一个参考子波与新的子波做比较来重复该过程，或者认为子波是目前的组成或者认为子波是属于一个新组。依上可知，根据勘探区域的大小，将来有成百上千的组。某些组可能只有一个成员，而其它的将会有许多的成员。用每组子波数量的最小值参数来限制自动子波分组产生的组个数。任何一个组子波个数比最小组子波个数少的被认为是0组。该组被放在输出范围结束的最低位，依靠色标的选择。组根据相似程序有序的排列。当重新排序时，相似组被挨着按照数字的大小排列。K方法使用选项性在自动子波分组中K方法调解定义地震道应用输出分组。因为有许多组，所以来自自动子波分组的输出看起来自然，当在二维平面图上显示时，从自动子波分组的输出可以看到进程。除非输入的数据体是由无噪音的模型数据体。K方法调解的第一步是使每一组的子波数相同，这需要依靠增大所有子波的M值，并且分得的结果由成员的个数定。一旦每一组子波的平均值被计算出来，所有的原始子波与子波的平均值做比较，从而适应最合适的门限。如果原始子波没有与任何一个子波组相适应，则该组被认为是0组。通过K方法的不断迭代，比较原始子波与参考子波的过程被定义为1。用K方法迭代次数参数在控制这个过程。通常作为平均数2次迭代是充足的。重新分组的选项依靠以上参数的输入，层间子波相似的差异，数据的信噪比比，在K方法调解以后，可以有或多或少的组。如果有太多的组，你将从K方法调解或以自动子波重新分组中输出依靠合并相似组的方法重新分组，。当相似性比较参考子波时，合并相似组的关键时增加用在Manhattan参考的门限。在非监督（门限K方法）波形分类中，这个过程是最初的依靠从1组中提出的参考子波与该组的其它子波做比较从而得到在门限内的相似值。如果对于1组的参考子波设一个子波在门限内，则该组的其它子波也被认为是1组，而且标志也被认可。下面不被认可的组参考子波与仍然不认可的子波组作比较。相匹配的子波被认为是2组。这个过程重复为每一个参考子波定位，创建一个小的但拥有关系的组。为了近一步确定组是恰当的合并，K方法调解被应用。K方法调解将要求更多的重复，从而达到平衡，K方法调解是由K方法迭代次数参数控制的。断层的该参数是2。分组过程的最后一步是重复的用一个不断增大的门限值重新分组。门限被增加直到该组的数值比最大的组数值小或相等、参数选择过程波形分类非监督门限K方法，然后点击参数进入参数对话框。K方法调节 当激活了K方法调节，K方法调节提炼从自动子波分组步骤中输出的道分组。K方法迭代的次数比较目标子波和平均的K方法子波的过程定义为K方法的一次迭代。定义你希望运行的K方法迭代次数通常求平均值二次迭代就充分了。K方法迭代次数的有效参数范围MinimumMaximum Default 2 1000 2重新分组当重新分组选项被激活，重新分组合并从K方法调节或自动子波分组步骤中输出的相似组。重新分组阶段估算在当前门限限制和应用K方法调节的每一子波的相似度。这被定义为重新分组约定。控制K方法调节迭代数每一次重新分组约定有K方法重新分组迭代数参数。K方法重新分组迭代次数的有效参数的范围MinimumMaximum Default 1 100 2如果在第一次重新分组约定之后：分组数仍然超过组参数的最大值，由门限增量控制随门限限制增大，该过程被重复着。门限增量的有限参数范围MinimumMaximum Default 0.01 0.50 0.1这个方程式被继续直到组数字输出少于或等于在流程中定义的最大分组数。最大分组数的有效参数范围MinimumMaximum Default 2 5000 1练习非监督（K方法）波形分类在先前的练习中，我们应用波形分类与周围地震数据作比较，这些地震数据都是来自顶部标记层位亮点的单子波。在这项练习中，我们将会运行非监督（K方法）波形反演看它是否可以识别与生产油场相一致的地质特征。看该方法与监督反演做比较。另外，我们正在探索用非监督波形分类的显示问题。1、 在开始前，关闭所有的地震视图和平面视图窗口。2、 在叠后的主要对话框里，在监督波形分类后加一个空的过程。3、 取消监督波形分类过程，将非监督（门限K方法）波形分类放入空的过程。4、 打开参数对话框得到下面的窗口。集中注意力在对话框的下半部分。非监督波形分类对这些参数是非常敏感的。所以当选项参数值时应慎重选项。门限的微小变动和扩展长度参数的微小变化、将会导致子波分组的彻底变化。非监督波形分类反演运行时间比监督反演的长。在监督反演中参数测试应该发生或者在地震数据体子集前运行。非监督波形分类是在整个勘测过程中。这项练习将用断层设置。增加一个例子是为了与当前的目的做比较。如果时间充足，你可以用以前的参数进行实验。5、 打开层位监督追踪分析选项并打开参数对话框。6、 选项顶部一标记作为层位的名称，设定时间宽度为40毫秒。7、 选项一个输出层位名称。8、 运行。从非监督波形分类中显示输出。非监督（门限K方法）波形分类输出层位值作为组的值。这种方法中，在输出层位中发现的最大值范围是由最大组值在先前参数剖面中的参数限制的。注意：最大组值参数设置的组的上限值可能输出，但层位的结果将比这个限制要小一些。1、 打开平面视图显示在非监督波形分类运行中创建的输出层位。2、 设置一个光谱色标并使显示范围与之相适应而且从O开始，保留下来的分类是为那些不能分组的道，到10结束，这种运行方式很可能出现组的最大值。就像在先前页提到的，组参数的最大数值提供了一个输出的上限。从层位的下拉菜单中用计算极限的选项核查该层位的最大值和最小值。在这个例子中，波形分类已经将子波分为6组。重新设置显示范围并注意显示的影响。非监督（门限K方法）波形分类已经可识别生产油场子波与子波之间的相似性。那些与色斑区域相一致的区域的识别是依靠监督反演。非监督波形分类的另一个优点是它可以使同一测区中所有道进行比较并分组是基于每一个的相似性，而不是寻找与参考子波相匹配的子波分组方式。用监督波形分类的方法是一种绘制与已知生产区域参考图件的好方法，但是可能会出现没有发现的油区的参考组未被探测到。修订门限1、 返回到非监督波形分类流程并显示该项工作的参数。2、 改变门限值从0.4到0.1。这将会增加严格分组子波的灵敏性分析。3、 改变输出层位的名称并反馈参数的变化。当我们测试门限参数值时，用一个习惯的名称帮助保证不同层位的追踪。这项测试的参数应与下图相匹配。 请注意：在这种情况下，当你依靠减小门限增加波形分类的灵敏度时，你将增大了运行时间。当工作完成了，打开并显示创建门限为0.1的地震层位。因为该数据有一个好的信噪比，用一个门限为0.1的值可提供一个可接受的结果。（噪音数据将会导致更多的道被放在0组）分布在测区内的子波相似性被证实是依靠属于0组的子波最小质量。着色区域暗示是由箭头所指生产区相一致（对应）的区域。当我们修改子波分类的灵敏度时，注意这些区域的影响。1、 重复先前的练习，改变门限值为0.25并改变输出层位的名称。2、 在一个新的平面图中，显示结果。你的层位应该与下面的相类似注意一般的集中边界与先前的例子做比较。同时注意来自己改变区域代表道的组数。比较该层位与先前层位的最大值与最小值。组内的子波分类将会大范围变化当你修改门限值时。1、 作为最后一个例子，返回到叠后的主要对话框，改变门限为1。2、 重新命名输出地震层位并运行。再一次注意在子波分组时门限的影响。当门限是增加的，灵敏度将会降低。换句话说，当门限接近1，最终有更多的子波被分为同一组，测区的所有子波将会被分为同一组。注意：黑色区域是未拾取层位控制的区域。在波形分类中将不会在该层位空白区分组。结论在该例中箭头所指示的区域与生产区域相对应。虽然当我们修改我们分析的灵敏度时，事实上我们仍可以看到参数测试暗示了一个均质程度相当高的材料组成的这些区域。你将从这些结果中得到结论，生产区也包括相似度很高的物质。这个信息在决定生产策略和开采潜能估计中是非常有用的。非监督（门限K方法）波形分类的例子下面的例子是展示K方法调节的影响和从非监督波形分类的输出中重新分组选项的影响。产生的地震层位平面图件只用自动子波分组。61152输入道波形分类输出142个组类型加上0组。选项练习1：自动子波分组如果时间充足创建上面的显示，开始你的叠后过程如下：1、 设置输入到demo163dv。用整个测区。2、 选项非监督（门限K方法）波形分类作为过程。3、 用Top_chalk_output层位引导层位时间窗并选项输入门，时间宽度为40定义窗。4、 在断层设置中离开门限和每组子波最小值。5、 显示先前参数并重选K方法调节和重新分组选项。6、