CN114252913A - 平面断层信息的识别方法及装置 - Google Patents
平面断层信息的识别方法及装置 Download PDFInfo
- Publication number
- CN114252913A CN114252913A CN202011022184.3A CN202011022184A CN114252913A CN 114252913 A CN114252913 A CN 114252913A CN 202011022184 A CN202011022184 A CN 202011022184A CN 114252913 A CN114252913 A CN 114252913A
- Authority
- CN
- China
- Prior art keywords
- image data
- difference
- slice
- initial
- processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012545 processing Methods 0.000 claims abstract description 89
- 238000001914 filtration Methods 0.000 claims abstract description 35
- 238000010606 normalization Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 238000010586 diagram Methods 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 230000001427 coherent effect Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 241000213006 Angelica dahurica Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003595 mist Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种平面断层信息的识别方法及装置,其中,该方法包括:获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,切片初始图像数据具有断层信息;对切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与切片初始图像数据之间的差异图像数据;对差异图像数据进行直方图均衡处理和阈值调整处理,以增强断层信息;根据断层信息增强后的图像数据识别切片初始图像数据中的断层信息。通过本发明,可以提高地震数据品质,增强局部断层区域与背景值的差异性,从而可以提高断层识别的精度。
Description
技术领域
本发明涉及地震数据处理领域,具体涉及一种平面断层信息的识别方法及装置。
背景技术
断层、裂缝、溶洞等是裂缝型油气藏如碳酸盐非均储层的重要组成部分,它们在油气成藏过程中发挥着重要作用。断层发育状况是控制储层油气“生、储、盖、圈、运、保”重要因素之一,因此世界上许多大型油气藏都与断层发育密切相关。同时,断层是重要的油气运移通道,它的存在既能够沟通分散的储集空间如溶洞体系形成规模化的储层,又能够为各个储集空间内部提供必要的渗透率。因此开展高精度的断层检测工作以支持后续的建模、油藏模拟等勘探开发研究,具有重要意义。
因此,断层识别是油气勘探中地震解释工作的重要一环。但是,常规叠后断层检测手段受制于地震数据品质的影响,长期以来存在两个方面的困扰。一方面,在地层破碎区域存在云雾状模糊信息对断层遮挡;另一方面,断层属性信息数值过于接近,无法区分相邻断层。
也就是说,由于地震数据品质较差,导致了无法有效地识别地震数据中的断层信息。
发明内容
有鉴于此,本发明提供一种平面断层信息的识别方法及装置,以解决上述提及的至少一个问题。
根据本发明的第一方面,提供一种平面断层信息的识别方法,所述方法包括:获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;对所述切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
根据本发明的第二方面,提供一种平面断层信息的识别装置,所述装置包括:数据获取单元,用于获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;中值滤波单元,用于对所述切片初始图像数据进行中值滤波处理;差异数据确定单元,用于确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;均衡处理单元,用于对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;识别单元,用于根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
根据本发明的第三方面,提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述方法的步骤。
根据本发明的第四方面,提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述方法的步骤。
由上述技术方案可知,通过对获取的目的层切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与切片初始图像数据之间的差异图像数据,之后对该差异图像数据进行直方图均衡处理和阈值调整处理,以增强断层信息,从而可以根据断层信息增强后的图像数据识别目的层切片初始图像数据中的断层信息,本技术方案通过提高地震数据品质,增强局部断层区域与背景值的差异性,从而可以提高断层识别的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本发明实施例的平面断层信息识别方法的流程图;
图2是根据本发明实施例的目的层切片初始图像;
图3是根据本发明实施例的对目的层切片初始图像进行增强处理后的图像;
图4是根据本发明实施例的平面断层信息识别装置的结构框图;
图5为本发明实施例的电子设备600的系统构成的示意框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在实现本发明的过程中,申请人发现如下相关技术:
断层检测技术的典型代表是相干系列技术。相干技术产生于上世纪90年代,目前已经发展三代算法。基于互相关的第一代算法(简称C1算法)由Bahorich和Frmer于1995年提出,利用多道相似性的第二代算法(简称C2算法)由Marfurt等于1998年提出,基于特征结构的第三代相干算法(简称C3算法)由Gersztenkorn和Marfurt提出。C1算法的应用前提条件较为苛刻,C2算法的缺陷是对波形敏感而且对横向振幅的变化不敏感,相比而言,C3算法弥补了上述两者的不足,它在偏移后的三维数据体中,对每一样点求得与周围数据的相干性,形成一个表征相干性的三维数据体,即计算时窗内的数据相干性。这样既可以压制连续性,突出不连续性,又可以定量地反映地震特征的横向变化,所得结果比地震水平切片的地质解释更直观。该方法主要应用于更客观、更细致的断层解释、河道、砂体及裂隙预测。但整体而言,相干属性技术的抗噪性能不佳,容易受到噪声污染,因此经常存在断层与周围背景对比度过低的情况,这导致结果图像对断层的刻画不够清晰,甚至影响解释人员对断层解释的判断。例如,在多期次构造运动形成的断层发育区域,偏移参数设置不当造成断层边界能量分散,导致断层信息呈现云雾状分布等等。
针对以上问题,本发明实施例提供了一种平面断层信息的识别方案,该方案以第三代相干属性二维切片(或者层位切片)数据为基础,通过增强局部断层区域与背景值的差异性,消除噪声污染,能够提高研究区的断层识别精度。以下结合附图来详细描述本发明实施例。
图1是根据本发明实施例的平面断层信息识别方法的流程图,如图1所示,该方法包括:
步骤101,获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息。
步骤102,对所述切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据。
在一个实施例中,可以先根据中值滤波处理后的切片图像数据与所述切片初始图像数据确定初始差异图像数据;之后,对所述初始差异图像数据进行归一化处理,得到灰度域的差异图像数据。
这里的对所述初始差异图像数据进行归一化处理,具体包括:根据所述初始差异图像数据确定差异最大值和差异最小值;之后根据所述差异最大值和差异最小值对所述初始差异图像数据进行归一化处理。
步骤103,对差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息。
具体地,可以对上述灰度域的差异图像数据进行直方图均衡处理;之后,根据预定规则对直方图均衡处理后的差异图像数据中的各像素进行阈值调整处理。
这里的预定规则可以是,根据断层与背景信息的差异来调整各像素的阈值。这里的背景信息是地层平缓、没有发生断裂现象时的相干属性数值,一般数值接近1。
步骤104,根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
通过对获取的目的层切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与切片初始图像数据之间的差异图像数据,之后对该差异图像数据进行直方图均衡处理和阈值调整处理,以增强断层信息,从而可以根据断层信息增强后的图像数据识别目的层切片初始图像数据中的断层信息,本发明实施例通过提高地震数据品质,增强局部断层区域与背景值的差异性,消除噪声污染,从而可以提高断层识别的精度。
为了更好的理解本发明,以下给出本发明实施例的平面断层信息识别的具体流程示例:
1)输入数据:该输入的数据为从三维地震数据体提取的第三代相干属性的二维时间切片(或者称为层位切片)图像U,U中各像素的值域浮动范围为[0,1]。
2)对相干属性的二维时间切片(或者层位切片)数据预处理:
a)对初始切片图像U进行中值滤波,得到中值滤波后的图像Um,滤波窗口大小可以根据断层尺度大小来设置,例如,可以设置为3*3、或者5*5等;
b)计算初始切片属性U与中值滤波后的图像Um中各像素的差值图像ΔU,其中,ΔU=U-Um;
c)将差值图像ΔU中各像素的值域归一化处理到灰度域的[0,255]范围内,得到归一化差异图像ΔUnormal。
具体地,计算公式如(1)所示,其中,round为四舍五入运算符,ΔUmin为ΔU的最小值,ΔUmax为ΔU的最大值,
ΔUnormal=round((ΔU-ΔUmin)/(ΔUmax-ΔUmin)*255) (1)
3)对归一化差异图像ΔUnormal进行常规的全局化直方图均衡化处理,得到图像Uenhanced。
通过直方图均衡处理,可以消除图像中云雾状的断层模糊区域,凸显该区域内部的线状断层结构。
4)根据断层与背景信息的差异,或者结合已知裂缝发育的井点处Uenhanced所展现的像素值,设置像素阈值,获得时间切片(或者层位切片)清晰的断层增强结果。
在实际操作中,在有井情况下,如果工区内有井,比如内A井,A井目的层的裂缝发育,假如A井处Uenhanced所展现的像素值为159,那么这里设置159为阈值,大于该数值的设置为255,小于等于该数值的都保留原数值。在无井情况下,根据地质经验,大型断层会在横向范围内有一定的影响范围,Uenhanced会在横向上远离断层的方向逐渐变大,此时选择合适横向距离处的Uenhanced数值作为阈值,之后做法参考有井的情况。
本发明实施例基于第三代相干属性结果的时间切片数据(或者层位切片数据),通过针对性的预处理方法,结合直方图均衡、阈值设置等手段进行增强处理,强化了断层的线状结构,提高了断层识别的精度。
为了进一步理解本发明实施例,以下给出一个实例。在该实例中,对某X工区三维叠后数据进行实验性断层增强分析,已知X工区内主断层分布趋势较为清晰,但是次生断层、溶洞等边界模糊不清。
根据本发明实施例,识别X工区中断层信息的流程具体包括:
1)提取X工区三维叠后地震数据体的第三代相干属性;
2)提取X工区第三代相干属性目的层切片图像数据,参见图2所示;
3)将提取的目的层切片属性执行本发明实施例提供的技术方案进行增强处理,增强处理后的图像可参见图3。
从图2和图3的对比分析可见,根据本发明实施例,研究区域断层识别精度获得了整体的提高,在原来相干结果存在大片模糊的区域,经过增强处理后,各级断层信息均得到了不同程度的提升,能够清晰展现断层的线性结构,从而可以提高断层识别的精度。
基于相似的发明构思,本发明实施例还提供一种平面断层信息的识别装置,优选地,该装置可用于实现上述方法实施例中的流程。
图4是该平面断层信息识别装置的结构框图,如图4所示,该装置包括:数据获取单元41、中值滤波单元42、差异数据确定单元43、均衡处理单元44和识别单元45,其中:
数据获取单元41,用于获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;
中值滤波单元42,用于对所述切片初始图像数据进行中值滤波处理;
差异数据确定单元43,用于确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;
均衡处理单元44,用于对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;
识别单元45,用于根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
通过中值滤波单元42对数据获取单元41获取的目的层切片初始图像数据进行中值滤波处理,差异数据确定单元43确定中值滤波处理后的切片图像数据与切片初始图像数据之间的差异图像数据,之后均衡处理单元44对该差异图像数据进行直方图均衡处理和阈值调整处理,以增强断层信息,从而识别单元45可以根据断层信息增强后的图像数据识别目的层切片初始图像数据中的断层信息,本发明实施例通过提高地震数据品质,增强局部断层区域与背景值的差异性,消除噪声污染,从而可以提高断层识别的精度。
具体地,上述差异数据确定单元43包括:初始差异确定模块和归一化模块,其中:
初始差异确定模块,用于根据中值滤波处理后的切片图像数据与所述切片初始图像数据确定初始差异图像数据;
归一化模块,用于对所述初始差异图像数据进行归一化处理,得到灰度域的差异图像数据。
在一个实施例中,上述归一化模块具体包括:差异值确定子模块和归一化子模块,其中,差异值确定子模块,用于根据所述初始差异图像数据确定差异最大值和差异最小值;归一化子模块,用于根据所述差异最大值和差异最小值对所述初始差异图像数据进行归一化处理。
上述均衡处理单元44具体包括:直方图均衡模块和阈值调整模块,其中:
直方图均衡模块,用于对所述灰度域的差异图像数据进行直方图均衡处理;
阈值调整模块,用于根据预定规则对直方图均衡处理后的差异图像数据中的各像素进行阈值调整处理。
上述各单元、各模块、各子模块的具体执行过程,可以参见上述方法实施例中的描述,此处不再赘述。
在实际操作中,上述各单元、各模块、各子模块可以组合设置、也可以单一设置,本发明不限于此。
本实施例还提供一种电子设备,该电子设备可以是台式计算机、平板电脑及移动终端等,本实施例不限于此。在本实施例中,该电子设备可以参照上述方法实施例进行实施及平面断层信息识别装置的实施例进行实施,其内容被合并于此,重复之处不再赘述。
图5为本发明实施例的电子设备600的系统构成的示意框图。如图5所示,该电子设备600可以包括中央处理器100和存储器140;存储器140耦合到中央处理器100。值得注意的是,该图是示例性的;还可以使用其他类型的结构,来补充或代替该结构,以实现电信功能或其他功能。
一实施例中,平面断层信息识别功能可以被集成到中央处理器100中。其中,中央处理器100可以被配置为进行如下控制:
获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;
对所述切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;
对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;
根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
从上述描述可知,本申请实施例提供的电子设备,通过对获取的目的层切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与切片初始图像数据之间的差异图像数据,之后对该差异图像数据进行直方图均衡处理和阈值调整处理,以增强断层信息,从而可以根据断层信息增强后的图像数据识别目的层切片初始图像数据中的断层信息,本发明实施例通过提高地震数据品质,增强局部断层区域与背景值的差异性,从而可以提高断层识别的精度。
在另一个实施方式中,平面断层信息识别装置可以与中央处理器100分开配置,例如可以将平面断层信息识别装置配置为与中央处理器100连接的芯片,通过中央处理器的控制来实现平面断层信息识别功能。
如图5所示,该电子设备600还可以包括:通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是,电子设备600也并不是必须要包括图5中所示的所有部件;此外,电子设备600还可以包括图5中没有示出的部件,可以参考现有技术。
如图5所示,中央处理器100有时也称为控制器或操作控件,可以包括微处理器或其他处理器装置和/或逻辑装置,该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中,存储器140,例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息,此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序,以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器,但并不限于此。
该存储器140可以是固态存储器,例如,只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器,其即使在断电时也保存信息,可被选择性地擦除且设有更多数据,该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142,该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143,该数据存储部143用于存储数据,例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100,以提供输入信号和接收输出信号,这可以和常规移动通信终端的情况相同。
基于不同的通信技术,在同一电子设备中,可以设置有多个通信模块110,如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132,以经由扬声器131提供音频输出,并接收来自麦克风132的音频输入,从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外,音频处理器130还耦合到中央处理器100,从而使得可以通过麦克风132能够在本机上录音,且使得可以通过扬声器131来播放本机上存储的声音。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时以实现上述平面断层信息识别方法的步骤。
综上所述,本发明实施例提供了一种平面断层信息识别方案,该方案抗噪性较强。通过对相干属性经过直方图均衡处理后,能够消除图像中云雾状的断层模糊区域,凸显该区域内部的线状断层结构,并且,计算过程简单,易于实现,能够实时提升二维平面内断层预测精度,且容易向三维扩展,支撑断层自动追踪等后续工作,具备挖掘部分隐含断层信息的能力。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的,因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种平面断层信息的识别方法,其特征在于,所述方法包括:
获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;
对所述切片初始图像数据进行中值滤波处理,并确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;
对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;
根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
2.根据权利要求1所述的方法,其特征在于,确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据包括:
根据中值滤波处理后的切片图像数据与所述切片初始图像数据确定初始差异图像数据;
对所述初始差异图像数据进行归一化处理,得到灰度域的差异图像数据。
3.根据权利要求2所述的方法,其特征在于,对所述初始差异图像数据进行归一化处理包括:
根据所述初始差异图像数据确定差异最大值和差异最小值;
根据所述差异最大值和差异最小值对所述初始差异图像数据进行归一化处理。
4.根据权利要求2所述的方法,其特征在于,对所述差异图像数据进行直方图均衡处理和阈值调整处理包括:
对所述灰度域的差异图像数据进行直方图均衡处理;
根据预定规则对直方图均衡处理后的差异图像数据中的各像素进行阈值调整处理。
5.一种平面断层信息的识别装置,其特征在于,所述装置包括:
数据获取单元,用于获取预定区域三维地震数据体的第三代相干属性目的层的切片初始图像数据,所述切片初始图像数据具有断层信息;
中值滤波单元,用于对所述切片初始图像数据进行中值滤波处理;
差异数据确定单元,用于确定中值滤波处理后的切片图像数据与所述切片初始图像数据之间的差异图像数据;
均衡处理单元,用于对所述差异图像数据进行直方图均衡处理和阈值调整处理,以增强所述断层信息;
识别单元,用于根据所述断层信息增强后的图像数据识别所述切片初始图像数据中的断层信息。
6.根据权利要求5所述的装置,其特征在于,所述差异数据确定单元包括:
初始差异确定模块,用于根据中值滤波处理后的切片图像数据与所述切片初始图像数据确定初始差异图像数据;
归一化模块,用于对所述初始差异图像数据进行归一化处理,得到灰度域的差异图像数据。
7.根据权利要求6所述的装置,其特征在于,所述归一化模块包括:
差异值确定子模块,用于根据所述初始差异图像数据确定差异最大值和差异最小值;
归一化子模块,用于根据所述差异最大值和差异最小值对所述初始差异图像数据进行归一化处理。
8.根据权利要求6所述的装置,其特征在于,所述均衡处理单元包括:
直方图均衡模块,用于对所述灰度域的差异图像数据进行直方图均衡处理;
阈值调整模块,用于根据预定规则对直方图均衡处理后的差异图像数据中的各像素进行阈值调整处理。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011022184.3A CN114252913A (zh) | 2020-09-25 | 2020-09-25 | 平面断层信息的识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011022184.3A CN114252913A (zh) | 2020-09-25 | 2020-09-25 | 平面断层信息的识别方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114252913A true CN114252913A (zh) | 2022-03-29 |
Family
ID=80790324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011022184.3A Pending CN114252913A (zh) | 2020-09-25 | 2020-09-25 | 平面断层信息的识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114252913A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130229891A1 (en) * | 2012-03-01 | 2013-09-05 | Transform Software & Services, Inc. | Method and system for image-guided fault extraction from a fault-enhanced seismic image |
CN103869358A (zh) * | 2014-02-18 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种基于直方图均衡的断层识别方法及设备 |
CN104459801A (zh) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | 用于识别断层的相干增强处理方法 |
CN107765301A (zh) * | 2017-10-13 | 2018-03-06 | 中国煤炭地质总局地球物理勘探研究院 | 煤层小断层的快速识别方法和装置 |
CN108957530A (zh) * | 2018-05-23 | 2018-12-07 | 电子科技大学 | 一种基于地震相干体切片的裂缝自动检测方法 |
CN109143348A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 三维地震数据断层强化处理方法 |
CN111582114A (zh) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种地震断层识别方法、装置、设备和存储介质 |
-
2020
- 2020-09-25 CN CN202011022184.3A patent/CN114252913A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130229891A1 (en) * | 2012-03-01 | 2013-09-05 | Transform Software & Services, Inc. | Method and system for image-guided fault extraction from a fault-enhanced seismic image |
CN103869358A (zh) * | 2014-02-18 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种基于直方图均衡的断层识别方法及设备 |
CN104459801A (zh) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | 用于识别断层的相干增强处理方法 |
CN109143348A (zh) * | 2017-06-28 | 2019-01-04 | 中国石油化工股份有限公司 | 三维地震数据断层强化处理方法 |
CN107765301A (zh) * | 2017-10-13 | 2018-03-06 | 中国煤炭地质总局地球物理勘探研究院 | 煤层小断层的快速识别方法和装置 |
CN108957530A (zh) * | 2018-05-23 | 2018-12-07 | 电子科技大学 | 一种基于地震相干体切片的裂缝自动检测方法 |
CN111582114A (zh) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | 一种地震断层识别方法、装置、设备和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106896406B (zh) | 基于阻抗域属性预测优质储层的方法和装置 | |
CN103792576A (zh) | 一种基于梯度结构张量的储层非均质检测方法及设备 | |
CN103869358A (zh) | 一种基于直方图均衡的断层识别方法及设备 | |
CN103278864A (zh) | 洞缝型储层的地质特征参数及分布的测定方法及装置 | |
CN112415591A (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
CN108629459B (zh) | 储层含烃孔隙的检测方法及装置 | |
CN108830140B (zh) | 一种基于电成像测井分形维数的火山岩岩性识别方法 | |
CN111582114B (zh) | 一种地震断层识别方法、装置、设备和存储介质 | |
CN108375785A (zh) | 裂缝带位置校正方法及装置 | |
CN114841927A (zh) | 一种页岩储层裂缝识别方法、装置、设备及存储介质 | |
US11467301B1 (en) | Method and system for identifying pore and fracture based on two-dimensional (2D) scan image of core | |
CN114252913A (zh) | 平面断层信息的识别方法及装置 | |
CN113077074B (zh) | 基于储层预测因子的储层预测方法及装置 | |
CN114255174A (zh) | 三维断层信息的识别方法及装置 | |
CN109948291A (zh) | 一种砂体不连续界线方向自适应识别方法 | |
CN114139229A (zh) | 断层增强方法和断层发育解释方法、存储介质及电子设备 | |
CN115184998A (zh) | 一种基于改进U-net神经网络的瑞利波频散曲线自动提取方法 | |
CN104314558B (zh) | 用斯通利波能量损失度判别储层有效性的方法 | |
CN112523740A (zh) | 用于超深海井下油气勘探的岩石类型识别方法及装置 | |
CN113586042A (zh) | 孔隙型碳酸盐岩动态渗透率确定方法及装置 | |
CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
CN114755740A (zh) | 岩石分布确定方法、装置、设备及介质 | |
CN108108512B (zh) | 一种储层岩相的表征方法及装置 | |
CN114578445A (zh) | 基于重力资料确定断裂位置的方法及装置 | |
CN113960657A (zh) | 地震数据特征获取方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |