一种缝洞连通域检测与表征方法技术

本发明专利技术公开一种缝洞连通域检测与表征方法，该方法包括以下步骤：S1提取三维地震数据的相干并计算其第一特征值，构建第一特征值三维数据体；S2选取五邻域作为连通域分析点邻域；S3连通域标记，标记算法采用双队列递归搜索的一次扫描法，步骤如下：S31寻找所有非零值位置组成大队列；S32在大队列中建立五邻域、连通域小队列并标记；S33从第一个点开始判断大队列数据点是否与五邻域连通，若连通则标记；S34去除异常值输出标记切片，记录连通域编号；S4通过纵向扫描和横向扫描标记连通域边界，计算形心，以之为中心点计算其与边界点的夹角与距离，实现边界点逆时针排序；S5计算连通域面积。本发明专利技术可提高缝洞解释精度与效率，为储层评价提供基础支撑。

【技术实现步骤摘要】
一种缝洞连通域检测与表征方法
本专利技术涉及地震资料处理与储层描述领域，特别涉及一种基于相干体第一特征值、五邻域双队列递归搜索的缝洞连通域检测与表征方法。
技术介绍
缝洞型储层是一类很重要的油气储层，在国内外油气勘探开发中具有很重要的地位。它的特点是储层孔隙度、渗透率高，缝洞连通性好，油气储量和产量高。在我国西部油区，此类储层分布很多，塔里木、长庆、新疆多个高产油气田均属于此类储层。缝洞型储层剖面有比较典型的特征，就是串珠状地震反射特征，但这种反射特征较强时才便于人工解释和拾取；而在地震水平切片上，缝洞的解释也受到缝洞的大小、数量、复杂边界、背景干扰等限制。
技术实现思路
基于上述技术问题，本专利技术提出一种基于相干体第一特征值、五邻域双队列递归搜索的缝洞连通域检测与表征方法，该方法可提高储层预测的精度和效率，为缝洞型油气藏勘探开发提供依据。本专利技术所采用的技术解决方案是：一种缝洞连通域检测与表征方法，包括以下步骤：S1提取三维地震数据的相干并计算其第一特征值，构建第一特征值数据体；S2选取本文档来自技高网...

【技术保护点】
1.一种缝洞连通域检测与表征方法，其特征在于包括以下步骤：/nS1提取三维地震数据的相干并计算其第一特征值，构建第一特征值数据体；/nS2选取五邻域作为连通域分析点邻域；/nS3连通域标记，标记算法采用双队列递归搜索的一次扫描法，步骤如下：/nS31寻找所有非零值位置组成大队列；/nS32在大队列中建立五邻域、连通域小队列并标记；/nS33从第一个点开始判断大队列数据点是否与五邻域连通，若连通则标记，若不连通大队列长度减一，递归直至大队列长度为零；/nS34去除异常值，输出标记切片，记录连通域编号；/nS4通过纵向扫描和横向扫描标记连通域边界，计算形心，以之为中心点计算其与边界点的夹角与距离，...

【技术特征摘要】
1.一种缝洞连通域检测与表征方法，其特征在于包括以下步骤：
S1提取三维地震数据的相干并计算其第一特征值，构建第一特征值数据体；
S2选取五邻域作为连通域分析点邻域；
S3连通域标记，标记算法采用双队列递归搜索的一次扫描法，步骤如下：
S31寻找所有非零值位置组成大队列；
S32在大队列中建立五邻域、连通域小队列并标记；
S33从第一个点开始判断大队列数据点是否与五邻域连通，若连通则标记，若不连通大队列长度减一，递归直至大队列长度为零；
S34去除异常值，输出标记切片，记录连通域编号；
S4通过纵向扫描和横向扫描标记连通域边界，计算形心，以之为中心点计算其与边界点的夹角与距离，实现边界点逆时针排序；
S5计算连通域面积。


2.根据权利要求1所述的一种缝洞连通域检测与表征方法，其特征在于，步骤S1具体内容如下：
对于某地震三维数据体，取相邻J道N个样点组成一个N×J的地震子体，用矩阵DN×J来表达，由D构成矩阵CJ×J：



计算C的第一特征值λ1，构建第一特征值数据体。


3.根据权利要求1所述的一种缝洞连通域检测与表征方法，其特征在于，步骤S31具体内容如下：
按水平切片或地层切片提取第一特征值数据，设置阈值系数，进行阈值处理，将小于阈值的置零；
先按行i再按列j，找出所有非零值位置，组成大队列；线号定义为LINE(k)，道号定义为CDP(k)；k＝1,2,…,NP；k为搜索到的非零值数据的序列号，NP为总的非零值数据个数。


4.根据权利要求1所述的一种缝洞连通域检测与表征方法，其特征在于，步骤S32具体内容如下：
从大队列第一点，i＝LINE(1),j＝CDP(1)开始，判断与其五邻域点是否连通，建立对应的最长为5个点的小队列；设第一个连通域标记为F＝1，记第F个连通域所包含的点数初值为NL(F)＝0，邻域小队列队首位置L＝1，标记Y(i,j)＝F；
(1)判识左侧点；若点(i,j-1)的值X(i,j-1)大于零且未标记，即X(i,j-1)>0&Y(i,j-1)＝0，则置Y(i,j-1)＝F，并从队列里找出该点，置邻域小队列位置点L＝L+1，设该邻域点连通标志G1＝1，NL(F)＝NL(F)+1；否则，置该邻域点不连通标志G1＝0,执行下一步(2)；
(2)判识右侧点；若X(i,j+1)>0&Y(i,j+1)＝0,则置Y(i,j+1)＝F，置邻域小队列位置点L＝L+1，设该邻域点连通标志G2＝1，NL(F)＝NL(F)+1，因为此点紧跟点(i,j)之后，所以不必查找和置换；否则，置该邻域点不连通标志G2＝0,执行下一步(3)；
(3)判识左侧下点；若X(i+1,j-1)>0&Y(i+1,j-1)＝0,则置Y(i+1,j-1)＝F，并从队列里找出该点，换置到邻域小队列位置点L＝L+1，设该邻域点连通标志G3＝1，NL(F)＝NL(F)+1；否则置该邻域点不连通标志G3＝0，执行下一步(4)；
(4)判识下方点；若X(i+1,j)>0&Y(i+1,j)＝0,则置Y(i+1,j)＝F，并从队列里找出该点，置邻域小队列位置点L＝L+1，设该邻域点连通标志G4＝1，NL(F)＝NL(F)+1；否则置该邻域点不连通标志G4＝0,执行下一步(5)；
(5)判识右侧下点；若X(i+1,j+1)>0&Y(i+1,j+1)＝0,则置Y(i+1,j+1)＝F，并从队列里找出该点，置邻域小队列位置点L＝L+1，设该邻域点连通标...

【专利技术属性】
技术研发人员：张军华，王静，刘震，陈永芮，赵杰，胡陈康，李琴，王作乾，常健强，任陆庆，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37