地震DNA地层超剥点线识别技术及应用

罗红梅

(1.中国石油化工股份有限公司胜利油田分公司勘探开发研究院,山东东营257015;2.中国石油化工股份有限公司胜利油田分公司博士后工作站,山东东营257015)

地震DNA地层超剥点线识别技术及应用

罗红梅1,2

(1.中国石油化工股份有限公司胜利油田分公司勘探开发研究院,山东东营257015;2.中国石油化工股份有限公司胜利油田分公司博士后工作站,山东东营257015)

摘要:地层油藏超剥点线的地球物理精细描述有利于提高含油气储层识别的可靠性和钻井成功率,但在地震描述中根据地震反射特征追踪的超剥点线与实际地层油藏的超剥点线往往存在较大的误差。介绍了地震DNA地层超剥点线识别技术的方法原理与实现步骤,通过将地震反射数据和蚂蚁体属性数据转换成相应的字符数据,采用地震DNA搜索方法将对应字符所表示的地震数据体中的超剥点识别出来。在此基础上,应用聚类算法将超剥点进行合理连线得到地层超剥线。模型合成地震记录的超剥点识别结果表明,地震DNA地层超剥点识别技术能够有效识别出超剥点,与传统瞬时相位识别的超剥点相比,该方法的识别精度更高。东营凹陷草桥地区的实际应用表明,该技术可以更加有效地描述超剥点线,实现了该区低序级不整合超剥边界的精细刻画,比基于瞬时相位属性描述的超剥线的识别精度明显提高。

关键词：地震DNA;字符数据;超剥点;聚类;超剥线

地层油藏是一种重要的隐蔽油气藏类型,在世界石油和天然气的产量和储量中,43%的石油储量和30%的天然气储量存在于地层圈闭中,但在地震描述中，根据地震反射特征追踪的地层尖灭线与实际地层油藏的尖灭线常常存在较大误差,地层油藏超剥线精细落实困难,导致钻探成功率低,仅在新疆的克拉玛依、准噶尔、南襄盆地的泌阳和渤海湾盆地等少数地区实现了成功钻探[1-4],如何准确落实地层油藏的超剥线位置是地层油藏勘探开发的关键。目前,在地层油藏勘探开发中,超剥点和超剥线描述技术发展较快,其中,地震反射夹角外推技术是利用不整合与地层的倾角差外推地层尖灭点的预测技术[5-6],该技术工作量大,局限性严重,难以得到广泛应用。利用地震属性进行超剥点和超剥线描述的技术[7-8],是沿剖面的不整合面或超覆层提取多种地震属性,通过对比分析确定用于地层圈闭识别的优势属性,如平均瞬时相位、波峰波谷振幅属性等,但是由于受地震资料分辨率限制,根据地震剖面和属性识别出来的地层超剥点与实际地层的超剥点之间存在着一定的误差。利用地质统计方法进行超剥点识别和描述的技术[9-10]，是对未剥蚀区的地层厚度数据拟合，然后外推确定超剥点位置，但该技术仅适用于被动大陆边缘、大型三角洲沉积区及坳陷型盆地。构造地质学中制作复原横剖面的地层厚度梯度法可以用于确定超剥点的位置[11],但其精度难以把握。周卫红等[12]、王章青等[13]、谢风猛等[14]、胡宇双等[15]应用精细地层对比、层序地层学方法确定地层超剥点位置,但其计算结果往往因人而异。波形分析技术[16]是通过计算“曼哈顿距离”来识别尖灭点的准确位置。所谓“曼哈顿距离”是指衡量地震波形与所选取参考子波非相似程度的参数，在尖灭点位置“曼哈顿距离”为最大值，其他位置则相对较小。地震DNA技术[17-19]能够有效提取数据体中的构造特点,为此,我们提出了地震DNA地层超剥点线识别技术,采用类似于生物遗传学基因组合识别的方式识别超剥点,并在此基础上,利用聚类算法将超剥线提取出来,最后用理论和实际数据验证本文提出的方法有效性。

1方法原理及实现

1.1地震DNA技术原理

在生物信息学中,经常用到DNA搜索技术去检测DNA分子中的碱基对序列,地震DNA技术就是受这种方法的启发而产生的一项新技术。该技术已经被成功应用于地震层位和层序的自动提取[18-19]。这种搜索技术的核心类似于文本搜索中的正则表达式。

地震DNA的实现过程主要有两步。第一步叫做“转化”,也就是将输入的数值转化为文本。通过指定每个字符为一个唯一范围的数据值将连续的数据转换为字符。例如:

(1)

式中:字符a,b,c分别对应于输入数据范围[min,-0.1),[-0.1,0.1),[0.1,max],min和max分别对应于输入数据的最小值和最大值。这样就完成了数据从数值到字符的转化过程,也就是将输入数据从模型域转化到DNA域(图1)。第二步就是设置搜索模式,这个搜索模式也叫做“基因”。众所周知,基因是由一系列“核苷酸”排序而成。目的是为了在庞大的DNA库中搜索出感兴趣的“基因”,实现方法是创建一个能够描述待寻找特征的正则表达式。图1中基因为a{nmin,nmax},b{nmin,nmax}和c{nmin,nmax},其中,字符a,b,c是字符转换中的字符,花括号中的数字表示对应字符重复次数的取值范围,花括号的左边表示最小重复次数,右边表示最大重复次数。例如一个简单的字符基因组合的正则表达式(DNA链):a{2,3}b{1,2}c{4}b{0,1}a{3}。如果使用在字符转换中描述的模式描述式(1),则地震数据范围在[min,-0.1)中的转换成字符a,范围在[-0.1,0.1)中的转换成字符b,范围在[0.1,max]中的字符转换成字符c。根据正则表达式(DNA链)的形式,字符a至少重复2次,至多重复3次,紧接着字符b至少重复1次,至多重复2次,然后字符c重复4次,接着字符b至少重复0次,至多重复1次,最后字符a至少重复3次,至多重复5次,即它们的组合有aabccccaaa,aaabccccaaa,aabbccccaaa,aaabbccccaaa,aabccccbaaa,aaabccccbaaa,aabbccccbaaa,aaabbccccbaaa这八种有效匹配方式。正则表达式搜索还允许字符为零的重复,这意味零重复也是一个有效的匹配,使用零重复这个特性,可以创建更复杂的正则表达式,让其中一些子特性在一个有效的匹配中出现或者不出现。

图1　地震DNA流程结构

地震DNA技术并不限定输入数据,其输入数据可以为地震振幅数据,也可以是任意一种地震属性数据。由于地震DNA可以在不同属性之中同时进行搜索,这种多属性的特征给此方法带来了很大的灵活性。

1.2地层超剥点识别方法

1.2.1方法的提出

传统的地震DNA搜索主要是应用于地震层位和层序的自动提取。因为地震层位是能产生连续反射的界面,在地震剖面中表现为相同或者相似的波形、波组或者波系结构。即传统的地震DNA搜索主要是搜索地震数据中相似的基因。而地层超剥点的识别与此特性相反。如图2所示,通过将地震反射振幅数据转换成字符数据,可以发现超剥点两侧明显都具有不同的地震DNA基因序列,或表现为基因类型突变,或表现为某个字符重复次数突变。基于此,我们创新性地将地震DNA技术应用于超剥点识别。

图2　东营凹陷超剥点处的DNA突变示意

1.2.2识别方法及其实现步骤

将地震反射振幅数据、相位数据或蚂蚁体属性数据转换成字符数据,发现字符序列即DNA在超剥点处表现为突变,即超剥点两侧DNA(基因类型或单个字符的重复次数)是不相同的,基于此特点,建立如图3所示的地震DNA超剥点识别流程,主要由下述8个步骤组成:

1) 输入三维地震数据;

2) 按照Inline或者Crossline方向读入二维地震剖面数据;

3) 输入不整合层位数据;

4) 地震数据预处理:根据地震资料的品质可选做,旨在提高信噪比和分辨率,如f-x域Cadzow滤波技术提高地震数据的信噪比[20];边界保持各向异性扩散滤波[21],使得超剥现象变得清楚;以及谱反演提高分辨率等[22];

5) 截取不整合面上下一定时窗内的地震数据或属性数据,并设定字符转换参数(字符转换个数与范围参数),将地震数据或属性数据转化为字符数据;

6) 设定搜索参数,包括相似度、容忍度和匹配度;

7) 搜索二维剖面的超剥点,记录超剥点信息;

8) 搜索三维数据的超剥点。

其中,字符转换个数决定图像的分辨率,而范围参数直接决定提取目标同相轴的完整性和清晰性。字符转换个数通常通过不同字符个数的对比试验来优选。对于范围参数,如果设定的值过小,追踪出来的同相轴中包含的干扰信号就会比较多,无法有效区分目标同相轴和干扰;如果设定的值过大,会导致目标同相轴不完整或不清晰,甚至没有目标同相轴,通常通过对比试验来优选范围参数。搜索参数中的相似度是对原始DNA算法的改进之一,该参数在原始DNA算法中是一种字符串的匹配,也就是给定一个字符串,它能将与该字符串相似的字符串搜索出来，即搜索的是相似字符串。而对于寻找的超剥点就不一样,超剥点两侧的字符串类型是不相似的,当搜索到突变点时,搜索停止,搜索目的就是寻找突变点。因此,在超剥点识别中该参数为需要检测的当前层位上下不同字符的个数,也就是表示需要检测的同相轴个数,可以通过在不整合层位之下字符串类型搭配的情况来寻找突变点。容忍度是对原始DNA算法的改进之二,表示横向上不满足相似条件的道数超过该值时被认为是超剥点,该参数的设定是为了减少误差,因为在搜索过程中难免会存在误差,假如字符串类型出现突变,而突变后的字符类型又与之前的一致,这时就得判断该突变点是超剥点还是因误差导致的。鉴于此,设定容忍度来有效控制误判的发生。匹配度是对原始DNA算法的改进之三,它表示的是相邻道的字符最大重复次数的差值即代表同相轴厚度的变化范围,一般在同相轴分叉的地方,同相轴的视周期较大,而分叉又是导致字符串突变的原因之一,在字符转换之后,字符的最大重复次数与同相轴的视周期息息相关,而在沿着不整合层位向下搜索同相轴的过程中,常常会出现字符串类型相同,但字符重复次数不同的情况,这就需要设定一个控制参数即匹配度来判断突变点。

图3　地震DNA超剥点识别流程

1.3地层超剥线识别方法

地层超剥线是由识别出来的超剥点连接而成,因此首先要对超剥点进行合理分类,将属于同一条超剥线的超剥点聚成一类,从而实现超剥线的提取。超剥点的聚类原则是将合理的超剥点分为一类,所谓合理,就是指超剥点之间的距离和方向满足一定的最优条件。由于超剥点之间的距离和方向对于聚类具有约束意义,因此可以根据这些约束条件对超剥点进行聚类。描述超剥点之间的接近程度用“距离”来度量。两个超剥点之间的距离越小,表示两者之间共同点越多,距离越大,共同点越小。在平面上,超剥点都是用线号和道号来表示,所以用欧式距离公式表示:

(2)

同时将两个类中最近的两个超剥点距离作为这两个类的距离,即最短距离原则[23]。

聚类流程如图4所示。包含5个步骤:

1) 输入超剥点数据;

2) 按照距离和方向的要求从第i个超剥点开始判断,i小于最大线号;

3) 给定一定条件进行聚类搜索,然后判断两超剥点间的距离是否满足条件,判断两超剥点间的倾角方向是否满足条件,否就回到步骤2),i=i+1;

4) 记录满足条件的超剥点信息;

5) 输出超剥线。

图4　超剥点聚类算法流程

2模型试算

图5为模拟实际地层沉积结构设计的地质模型与相应的地震响应,自下而上发育5套厚度逐渐减薄的砂体,砂岩速度为2600～2900m/s,泥岩速度为2400～2700m/s,频率为30Hz。不整合面角度和剥蚀地层的角度见图5a。图5b中还给出了5套砂体超剥层试算结果,很明显5套砂体实际的超剥点与地震响应记录上的超剥点存在着一定的误差,而且不难看出,随着砂体厚度的减薄,这种距离误差随之增大(表1)。

图5　地质模型(a)与相应的地震响应(b)

剥蚀地层厚度/m实际剥蚀点位置/m传统瞬时相位识别结果基于地震DNA识别结果剥蚀点位置/m误差/m剥蚀点位置/m误差/m5037045080409.539.530730830100777.047.020113012801501218.088.020137015201501470.0100.010155019003501774.5224.5

针对上述地质模型,将正演模拟得到的合成地震记录作为地震反射振幅数据输入,该数据未经过任何地震预处理,采用本文方法进行超剥点识别(图6)。图6中从上往下依次为地质模型的地震响应记录、地震DNA超剥点识别结果和传统瞬时相位属性的识别结果。对比发现,两种方法识别的结果都比较准确,每一个识别出来的点都有效对应于一个砂体的超剥点。为了更加清楚地说明地震DNA识别算法的精度,对识别结果进行了定量统计(表1),进一步分析了超剥点识别误差的大小和造成误差的原因。

不难发现,无论是传统瞬时相位还是基于地震DNA识别的结果都存在着一定的误差,这是由于在地震剖面上地层超剥点线附近砂层的厚度明显减薄,同时受地震资料分辨率的限制,地震反射往往提前变弱或者消失。从表1中也能看出随着剥蚀地层厚度的减小,传统瞬时相位与基于地震DNA识别结果的误差都是增大的。地层厚度越薄,地震资料的分辨率就会越低,识别误差就会越大。尽管如此,地震DNA地层超剥点识别技术比传统瞬时相位识别的精度要高。这是因为地震DNA方法通过字符的转换参数控制可以有效增加地层超剥点线附近的字符,从而间接增强了地震反射。对比图6中地震响应和搜索结果可以看出地层超剥点线附近地震反射明显增强。

图6　超剥点识别结果对比

为了进一步验证方法的有效性,建立了更为复杂的地质模型(图7)进行分析。该模型给出了两个不整合层位和不整合层位下超剥点的实际位置,超剥点的实际位置是以一个宽度范围给出的。其中,Hor1_1级不整合和Hor2_2级不整合面之下分别发育7套和5套单层厚度相当的砂体。对该模型的正演模拟地震记录进行地震DNA超剥点识别,图8和图9是分别对Hor1_1级不整合和Hor2_2级不整合层位下超剥点的识别结果。可以看出,对于1级不整合层位和2级不整合层位之下的每一个砂体,其对应的超剥点都被很好地识别出来,识别结果刻画模型砂体比较准确。

图7　复杂地质模型地震响应

图8　Hor1_1级不整合层位下超剥点的识别结果

图9　Hor2_2级不整合层位下超剥点的识别结果

3实例分析

3.1地震数据的预处理

选取东营凹陷草桥南部三维工区的实际资料进行了地层超剥点、线的识别与描述。在地震振幅数据体中,存在断层、超剥结构等等。断层处基因会发生突变,因此,在断裂系统比较发育的地区进行地震DNA超剥点线识别之前应首先进行断层的提取工作,从而消除断层的影响,目前断层的识别方法有很多[24],如相干体技术[25]、波形分析技术[26]等,而现今最常用的自动断层提取方法是蚂蚁追踪算法[27]。考虑到地震反射同相轴与断层数据体中断层的相似性,可以利用蚂蚁追踪算法直接对地震振幅数据体进行地震同相轴追踪,进而提高地震同相轴对尖灭点的分辨能力。

在进行超剥点线识别之前,首先对地震数据进行蚂蚁体追踪预处理。图10a为草桥地区东西向地震振幅数据剖面,图10b为蚂蚁追踪预处理的蚂蚁体属性剖面,对比可以看出,蚂蚁体属性剖面上同相轴连续性明显更高,这充分说明蚂蚁追踪算法对于提高地层连续性是有效的。以此数据体为基础,进行地震DNA分析,可以避免断层的干扰,更好地识别出真实的超剥点线。

图10　二维地震剖面(a)及蚂蚁体属性剖面(b)

3.2草桥地区超剥点线识别

首先利用常规瞬时相位方法对该地区地层超剥线进行识别,再应用本文提出的地震DNA地层超剥点线识别技术进行识别。比较两种方法的识别结果,以验证本文方法的正确性和有效性。

从草桥地区过4口已钻井的东西向地震剖面来看(图11),在沙二段底部发育了一个二级不整合T3(黄色标明的位置),为一个沉积不整合面,由于其角度较小,并且上下波阻抗差异小,地震上表现为连续性较差的中等振幅反射。从利用瞬时相位识别的二级不整合T3之下剥蚀地层边界结果来看(图12a),可以大致识别出几条属性相近可以连起来的条带,图中画出的趋势线就认为是地层超剥线的展布。但识别结果并不太清楚,无法确切、有效地刻画超剥线的具体位置。进一步地,先输入地震蚂蚁体数据,然后利用地震DNA地层超剥点线识别技术来识别二级不整合T3之下剥蚀地层的边界。首先识别出地层超剥点,进而利用超剥点聚类算法得到地层超剥线(图12b),图中展示的近东西向粗红线为断层,其两侧基因类型为ab-ba或ba-ab。其它颜色的细线为地层超剥线,刻画完整,展布特征明显。分析典型地震剖面的地震DNA识别结果(图13)可以看到,每个超剥点左右两侧的DNA链(基因类型和单个字符的重复次数)与断层两侧基因差异较大。因此,地震DNA地层超剥点线识别结果清晰,每一个识别出来的超剥点均有效对应一条超剥线,说明本文所提出方法应用于地层超剥点、超剥线的识别非常有效。

图11　过4口井近东西向地震剖面

图12　草桥地区地层超剥线瞬时相位识别(a)与地震DNA识别(b)结果对比

图13　草桥地区H线地震DNA识别结果

4结论

本文提出了一种地震DNA地层超剥点线识别与描述方法,该方法将振幅和蚂蚁体属性数据转换为字符数据,通过一个正则表达式的模式识别来进行字符数据的搜索以获得超剥点,并通过聚类提取超剥线。通过模型数据和实际资料的处理,以及与传统的瞬时相位识别方法对比可以看出,该方法可以提高地层超剥点和超剥线描述的精度,尤其对于剥蚀地层厚度小于20m的情况下,识别精度提高明显。

与传统的瞬时相位识别超剥点相比,地震DNA方法具有更强的灵活性。地震DNA方法不限定输入数据,也就是说可以引入任意一种对识别有帮助的属性进行处理,甚至可以多属性同时识别。在本文研究中,仅仅使用振幅数据和蚂蚁体数据进行处理,下一步可以考虑加入更多属性同步处理,比如加入相干体和曲率属性等。同时,地震DNA作为一种灵活性很强的方法技术,在地震资料解释工作中也具有广泛的应用前景,下一步可以考虑将其应用于某些特定地质体的自动识别。

参考文献

[1]周小进,陈跃,黄继文,等.库车坳陷南缘地层超覆圈闭特征及含油气前景[J].石油实验地质,2000,22(1):40-45

ZHOU X J,CHEN Y,HUANG J W,et al.Characteristics of stratigraphic overlap traps and petroleum prospect in the south margin of Kuqa depression[J].Experimental Petroleum Geology,2000,22(1):40-45

[2]陈宗清.川东石炭系地层—构造复合圈闭气藏的预测、分类和评价[J].石油学报,2000,21(3):25-30

CHEN Z Q.Forecast,classification and assesment for gas reservoirs of carboniferous stratigraphic-structural combination traps in eastern Sichuan[J].Acta Petrolei Sinica,2000,21(3):25-30

[3]李丕龙.陆相断陷盆地油气地质与勘探(卷五)[M].北京:石油工业出版社,2003:31-39

LI P L.Petroleum geology and exploration of continental fault basin (book five)[M].Beijing:Petroleum Industry Press,2003:31-39

[4]隋风贵,赵乐强,林会喜,等.济阳坳陷第三系地层油藏形成机制与组合模式[J].油气地质与采收率,2005,12(6):1-5

SUI F G,ZHAO L Q,LIN H X,et al.The forming mechanism and its combination mode of stratigraphic reservoir in Tertiary in Jiyang Depression[J].Petroleum Geology and Recovery Efficiency,2005,12(6):1-5

[5]苏朝光,闫昭岷,张营革,等.地层油藏超剥尖灭线夹角定量外推方法模型研究[J].地球物理学进展,2007,22(6):1841-1846

SU C G,YAN Z M,ZHANG Y G,et al.The study of quantitive extrapolation model for overlap-denudation pitching-out line stratigraphic reservoir[J].Progress in Geophysics,2007,22(6):1841-1846

[6]张福利,赵丽平.定量确定地层超覆线位置的地震反射夹角外推法[J].勘探地球物理进展,2008,31(4):278-280

ZHANG F L,ZHAO L P.Quantitative determination of overlap line position via extrapolation of angle of seismic reflection[J].Progress in Exploration Geophysics,2008,31(4):278-280

[7]管晓燕,穆星,刘军.相位分析技术在地层圈闭描述中的应用[J].油气地质与采收率,2007,14(6):56-58

GUAN X Y,MU X,LIU J.Application of the phase analysis technique for stratigraphic trap description[J].Petroleum Geology and Recovery Efficiency,2007,14(6):56-58

[8]鲁红利.哈山南缘浅层超剥带圈闭特征研究[J].内江科技,2012,33(12):150-146

LU H L.Research of pitching-out belt trap in the south shallow of Hasan[J].Neijiang Technology,2012,33(12):150-146

[9]钟广法,周祖翼,耿建华.有效厚度拟合定量恢复原始地层超覆点位置[J].沉积学报,2003,21(4):614-618

ZHONG G F,ZHOU Z Y,GENG J H.Extrapolating the original positions of removed onlap-offlap points using a piecewise-constructed model of effective stratal thickness[J].Acta Sedimentological Sinica,2003,21(4):614-618

[10]张营革,田建华.沾车地区馆陶组地层超覆油藏地震描述技术应用及效果分析[J].石油物探,2003,42(3):322-327

ZHANG Y G,TIAN J H.Seismic description technology for stratigraphic overlap reservoir of the Guantao Formation and its effect analysis in Zhanche area[J].Geophysical Prospecting for Petroleum,2003,42(3):322-327

[11]江涛,苗洪波,王芳,等.松辽盆地西部盆缘带地层超覆油藏形成条件[J].大庆石油地质与开发,2006,25(4):24-26

JIANG T,MIAO H B,WANG F,et al.Generation condition for overlap reservoir in border belt of western Songliao basin[J].Petroleum Geology and Oilfield Development in Daqing,2006,25(4):24-26

[12]周卫红,谢风猛,刘健.东营凹陷西北部地层超覆油藏勘探技术研究[J].承德石油高等专科学校学报,2005,7(3):23-26

ZHOU W H,XIE F M,LIU J.Exploration technology research on stratigraphic overlapped oil reservoirs in northwest of Dongying depression[J].Journal of Chengde Petroleum College,2005,7(3):23-26

[13]王章青,杨立战,林胜春.地震综合解释技术在地层超覆圈闭研究中的应用[J].江汉石油学院学报,2003,25(2):60-61

WANG Z Q,YANG L Z,LIN S C.Application of integrated seismic interpretation technique in studying stratigraphic overlap trap[J].Journal of Jianghan Petroleum Institute,2003,25(2):60-61

[14]谢风猛,李彦强.滨南油区地层超覆油藏的勘探[J].油气地质与采收率,2004,11(6):27-30

XIE F M,LI Y Q.Exploration of the stratigraphic overlap oil reservoirs in Binnan petroliferous province[J].Petroleum Geology and Recovery Efficiency,2004,11(6):27-30

[15]胡玉双,李黎明.层序地层学在松辽东部断陷有利勘探区带预测中的应用[J].大庆石油地质与开发,1999,18(2):10-12,26

HU Y S,LI L L.Application of sequence stratigraphy to the prediction of favorable exploration areas in the east rift of Songliao Basin[J].Petroleum Geology and Oilfield Development in Daqing,1999,18(2):10-12,26

[16]李海亮,陈启林,朱允辉,等.波形分析技术在识别小断层和剥蚀线中的应用[J].新疆石油地质,2007,28(2):232-234

LI H L,CHEN Q L,ZHU Y H,et al.Application of wave analysis to identifying minor fault and nip-out line[J].Xinjiang Petroleum Geology,2007,28(2):232-234

[17]BAKKE J,GRAMSTAD O,SØNNELAND L.Seismic DNA——a novel visually guided search method using non-local search and multi attribute data sets[J].Expanded Abstract of 82ndAnnual Internat SEG Mtg,2012:589-593

[18]GRAMSTAD O,BAKKE J,SØNNELAND L,et al.Simultaneous extraction of stratigraphic sequences using iterative seismic DNA detection[J].Interpretation,2014,2(4):T167-T176

[19]GRAMSTAD O,BAKKE J,SØNNELAND L.Seismic surface extraction using iterative seismic DNA detection[J].Expanded Abstract of 82ndAnnual Internat SEG Mtg,2012:600-604

[20]崔树果,朱凌燕,王建花.f-x域Cadzow技术分块压制随机噪声及其应用[J].石油物探,2012,51(1):43-50

CUI S G,ZHU L Y,WANG J H.Random noise attenuation with Cadzow technique inf-xdomain and its application[J].Geophysical Prospecting of Petroleum,2012,51(1):43-50

[21]严哲,顾汉明,蔡成国.基于各向异性扩散滤波的地震图像增强处理[J].石油地球物理勘探,2013,48

(3):390-394

YAN Z,GU H M,CAI C G.Seismic image enhancement based on anisotropic diffusion[J].Oil Geophysical Prospecting,2013,48(3):390-394

[22]柴新涛,李振春,韩文功,等.基于LSQR算法的谱反演方法研究[J].石油物探,2012,51(1):11 -18

CHAI X T,LI Z C,HAN W G,et al.Spectral inversion method analysis based on the LSQR algorithm[J].Geophysical Prospecting for Petroleum,2012,51(1):11-18

[23]向东进,李宏伟,刘小雅.实用多元统计分析[M].武汉:中国地质大学出版社,2005:1-184

XIANG D J,LI H W,LIU X Y.Applied multivariate statistical analysis[M].Wuhan:China University of Geosciences Press,2005:1-184

[24]王鹏,钟建华,张亚金,等.断层识别技术及其在肇州油田的应用[J].石油物探,2011,50(5):521-525

WANG P,ZHONG J H,ZHANG Y J,et al.Fault identification technique and its application in Zhaozhou oilfield[J].Geophysical Prospecting for Petroleum,2011,50(5):521-525

[25]孙夕平,杨国权.三维地震相干体技术在目标区沉积相研究中的应用[J].石油物探,2004,43(6):591-594

SUN X P,YANG G Q.Application of 3D seismic coherence cube technology to target area sedimentary facies research[J].Geophysical Prospecting for Petroleum,2004,43(6):591-594

[26]师永民,祁军,张成学,等.应用地震波形分析技术预测裂缝的方法探讨[J].石油物探,2005,44(2):128-130

SHI Y M,QI J,ZHANG C X,et al.Discussion of the method of applying seismic waveform analysis techniques to predict fracture[J].Geophysical Prospecting for Petroleum,2005,44(2):128-130

[27]严哲.三维地震断层自动识别与智能解释[D].武汉:中国地质大学,2010

YAN Z.Fault automatic recognition and intelligent interpretation in 3-D seismic data[D].Wuhan:China University of Geosciences,2010

(编辑：朱文杰)

Stratigraphic pitching-out point and line identification technology based on seismic DNA and its application

LUO Hongmei1,2

(1.ResearchInstituteofPetroleumExplorationandDevelopment,SinopecShengliOilfieldCompany,Dongying257015,China;2.PostdoctoralWorkstation,ShengliOilfield,SINOPEC,Dongying257015,China)

Abstract:The geophysical description of pitching-out point and line in stratigraphic reservoir is beneficial to improve the reliability of oil and gas reservoir recognition and success rate of well drilling.The pitching-out point and line is traced by seismic reflection characteristics.There is a large difference between the pitching-out point and line and the actual point and line.In this paper we introduce the method principle and implementation steps of stratigraphic pitching-out point and line identification technology based on seismic DNA.The seismic reflection data and ant-attribute data is converted into the corresponding character data,and then the pitching-out point from the character data indicated corresponding seismic data is identified using seismic DNA.On the basis,stratigraphic pitching-out points are reasonably linked to be the pitching-out line by using clustering.Model simulated records show that stratigraphic pitching-out point and line identification technology based on seismic DNA can effectively identify pitching-out point,and can improve the identification precision compared with the result of conventional instantaneous phase.The application result of actual data in Caoqiao area,Dongying depression shows that this technology can effectively describe pitching-out line and the fine description of the unconformity pitching-out boundary for thin layer is realized.

Keywords:seismic DNA,character data,pitching-out point,cluster,pitching-out line

收稿日期：2015-10-22;改回日期：2016-01-16。

作者简介：罗红梅(1973—),女,高级工程师,博士,主要从事地震地质数据综合处理与解释方面的研究。

基金项目：中国石油化工股份有限公司科技开发部“东营凹陷南部地层超剥带精细描述与目标评价”项目(P13039)资助。

中图分类号：P631

文献标识码：A

文章编号：1000-1441(2016)03-0414-11

DOI:10.3969/j.issn.1000-1441.2016.03.012

This project is financially supported by the Project of SINOPEC Development Planning Department (Grant No.P13039).