CN107918151A - 基于敏感方位地震数据的断裂识别方法及装置 - Google Patents
基于敏感方位地震数据的断裂识别方法及装置 Download PDFInfo
- Publication number
- CN107918151A CN107918151A CN201710978827.3A CN201710978827A CN107918151A CN 107918151 A CN107918151 A CN 107918151A CN 201710978827 A CN201710978827 A CN 201710978827A CN 107918151 A CN107918151 A CN 107918151A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- fault
- target zone
- coherent
- sensitive
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000001427 coherent effect Effects 0.000 claims abstract description 95
- 230000007480 spreading Effects 0.000 claims abstract description 48
- 238000003892 spreading Methods 0.000 claims abstract description 48
- 238000011161 development Methods 0.000 claims abstract description 44
- 238000011160 research Methods 0.000 claims abstract description 29
- 230000004927 fusion Effects 0.000 claims abstract description 25
- 238000005520 cutting process Methods 0.000 claims abstract description 6
- 238000004590 computer program Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 238000010276 construction Methods 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 6
- 239000012634 fragment Substances 0.000 claims description 3
- 238000010030 laminating Methods 0.000 claims description 3
- 238000005336 cracking Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 14
- 238000003860 storage Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供了一种基于敏感方位地震数据的断裂识别方法及装置,该方法包括:确定研究区的断裂发育展布方向;将研究区的地震数据按照方位角划分成偶数个数据分组,并将数据分组沿相同走向方向两两叠加,获得多个新地震数据体;获取新地震数据体的目的层沿层相干切片;根据断裂发育展布方向及所述切片,确定对断裂成像敏感的新地震数据体;对对断裂成像敏感的新地震数据体作剖面断裂解释,并根据解释结果对相应的切片作平剖面吻合;融合吻合后的切片,使不同敏感方位识别出的断裂展现于同一张切片上,并根据其上断裂的接触关系作目的层断裂的平面组合,获得目的层断裂。本申请可提高复杂区块断裂识别的准确度。
Description
技术领域
本申请涉及地球物理勘探技术领域,尤其是涉及一种基于敏感方位地震数据的断裂识别方法及装置。
背景技术
在一些地下地质条件复杂区块(诸如尼日尔Bilma区块等),常有断块破碎、断裂的密集发育,特别是构造部位,往往多条断裂相互交错,断裂关系复杂,因而断裂的平面组合具有多解性,主断裂和次级断裂的关系也不十分明确,从而容易造成圈闭识别的不确定性。
由此可见,在复杂区块的圈闭识别中,断裂识别非常关键。目前常用的断裂识别方法有相干体、曲率体等。然而,对于地下地质条件复杂的区块,这些方法无法准确地识别出断裂,从而影响了圈闭识别的准确性。因此,目前亟需一种新的断裂识别方案,以实现对复杂区块断裂的准确识别。
发明内容
本申请实施例的目的在于提供一种基于敏感方位地震数据的断裂识别方法及装置,以实现对复杂区块断裂的准确识别。
为达到上述目的,一方面,本申请实施例提供了一种基于敏感方位地震数据的断裂识别方法,包括:
确定研究区的断裂发育展布方向;
将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
获取每个新地震数据体的目的层沿层相干切片;
根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
本申请实施例的基于敏感方位地震数据的断裂识别方法,将所述研究区的地震数据按照方位角划分成偶数个数据分组,包括:
以所述地震数据正北方向为零度基准,按照360度方位将所述地震数据按照方位角等分成偶数个数据分组。
本申请实施例的基于敏感方位地震数据的断裂识别方法,所述两两叠加至少包含沿垂直于所述断裂发育展布方向的叠加,以及沿平行于所述断裂发育展布方向的叠加。
本申请实施例的基于敏感方位地震数据的断裂识别方法,所述获取每个新地震数据体的目的层沿层相干切片,包括:
分别对所述多个新地震数据体进行倾角扫描,并利用倾角扫描结果进行构造导向滤波;
将构造导向滤波后的多个新地震数据体进行相干处理,对应获得多个相干数据体;
对每个相干数据体,沿目的层进行沿层相干提取,对应获得多个目的层沿层相干切片。
本申请实施例的基于敏感方位地震数据的断裂识别方法,所述根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体,包括:
确定其方位角与所述断裂发育展布方向垂直的目的层沿层相干切片,则该目的层沿层相干切片所对应新地震数据体为对断裂成像敏感的新地震数据体。
另一方面,本申请实施例还提供了一种基于敏感方位地震数据的断裂识别装置,包括:
断裂展布方向确定模块,用于确定研究区的断裂发育展布方向;
地震数据分组叠加模块,用于将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
目的层沿层相干切片获取模块,用于获取每个新地震数据体的目的层沿层相干切片;
敏感方位数据确定模块,用于根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
断裂解释及平剖面吻合模块,用于对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
切片融合及断裂获取模块,用于对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
本申请实施例的基于敏感方位地震数据的断裂识别装置,将所述研究区的地震数据按照方位角划分成偶数个数据分组,包括:
以所述地震数据正北方向为零度基准,按照360度方位将所述地震数据按照方位角等分成偶数个数据分组。
本申请实施例的基于敏感方位地震数据的断裂识别装置,所述两两叠加至少包含沿垂直于所述断裂发育展布方向的叠加,以及沿平行于所述断裂发育展布方向的叠加。
本申请实施例的基于敏感方位地震数据的断裂识别装置,所述获取每个新地震数据体的目的层沿层相干切片,包括:
分别对所述多个新地震数据体进行倾角扫描,并利用倾角扫描结果进行构造导向滤波;
将构造导向滤波后的多个新地震数据体进行相干处理,对应获得多个相干数据体;
对每个相干数据体,沿目的层进行沿层相干提取,对应获得多个目的层沿层相干切片。
本申请实施例的基于敏感方位地震数据的断裂识别装置,所述根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体,包括:
确定其方位角与所述断裂发育展布方向垂直的目的层沿层相干切片,则该目的层沿层相干切片所对应新地震数据体为对断裂成像敏感的新地震数据体。
再一方面,本申请实施例还提供了另一种基于敏感方位地震数据的断裂识别装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
确定研究区的断裂发育展布方向;
将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
获取每个新地震数据体的目的层沿层相干切片;
根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
由以上本申请实施例提供的技术方案可见,本申请实施例首先确定研究区的断裂发育展布方向;然后将研究区的地震数据按照方位角划分成偶数个数据分组,并将数据分组沿相同走向方向两两叠加,获得多个新地震数据体;其次获取新地震数据体的目的层沿层相干切片,并根据断裂发育展布方向及上述的切片,确定其中对断裂成像敏感的新地震数据体;然后对对断裂成像敏感的新地震数据体作剖面断裂解释,并根据解释结果对相应的切片作平剖面吻合;然后融合吻合后的切片,使不同敏感方位识别出的断裂展现于同一张切片上,并根据其上断裂的接触关系作目的层断裂的平面组合,获得目的层断裂。由于本申请实例是通过利用敏感方位角叠加的地震数据进行断裂识别,并通过沿层相干融合切片对识别出的断层进行精细的平面断裂组合,而得到目的层断裂的,因此,本申请实施例提高了复杂区块断裂识别的准确度,从而可为下一步圈闭识别奠定了良好的断裂解释基础,进而有利于提高圈闭的有效性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请一实施例的基于敏感方位地震数据的断裂识别方法的流程图;
图2为本申请一实施例中某研究区主要断裂走向示意图;
图3为本申请一实施例中分方位叠加与全叠加效果对比示意图;
图4为本申请一实施例中分方位叠加800ms相干切片对比示意图;
图5为本申请一实施例中分方位叠加E5目的层沿层相干切片对比示意图;
图6为本申请一实施例中垂直于断裂走向方位数据与全叠加剖面对比示意图;
图7为本申请一实施例中E5断裂特征值融合沿层相干与全叠加相干对比示意图;
图8为本申请一实施例中分方位叠加与全方位叠加数据断裂解释平面组合对比示意图;
图9为本申请一实施例的基于敏感方位地震数据的断裂识别装置的结构框图;
图10为本申请另一实施例的基于敏感方位地震数据的断裂识别装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
参考图1所示,本申请实施例的基于敏感方位地震数据的断裂识别方法可以包括以下步骤:
S101、确定研究区的断裂发育展布方向。
在本申请一些实施方式中,可首先对研究区以往的研究成果(例如区域地质概况、构造解释成果等)进行断裂研究,从而确定研究区断裂展布方向,并可以正北方向为0度基准确定研究区断裂发育展布方向的度数。以尼日尔某工区为例,通过前期研究发现本研究区区域断裂走向以北西--南东向展布为主,以正北方向为0度基准确定研究区断裂发育展布方向的度数大概在140°左右,垂直断裂方向在50°左右,如图2所示。
S102、将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体。
在本申请一些实施方式中,将不同方位的数据叠加时,首先要考虑信噪比,如果划分的数据分组过多,必然造成每份叠加资料的覆盖次数减少,从而容易造成信噪比降低;如果划分的数据分组过少,又会类似全方位数据叠加,从而容易降低不同方位的独特性和敏感性。因此,数据分组数量应综合考虑以上因素来确定。在本申请一示例性实施方式中,可以所述地震数据正北方向为零度基准,按照360度方位将所述地震数据按照方位角等分成n个(n为偶数,例如4,8,16等)数据分组。然后,将n份数据分组沿共走向方向两两相加可得到n/2个新地震数据体。
在本申请一示例性实施方式中,以尼日尔某工区为例,以正北方向为0度基准按360度方位将其宽方位角地震资料均分成8份,即345~30度、30~75度、75~120度、120~165度、165~210度、210~255度、255~300度、300~345度,将8份地震数据体沿共走向方向两两相加得到4个地震数据体,分别为数据体A、数据体B、数据体C、数据体D。所述数据体A为345~30度和165~210度叠加数据体,数据体B为30~75度和210~255度叠加数据体,数据体C为75~120度和255~300度叠加数据体,数据体D为120-165度和300-345度叠加数据体。例如图3所示,四个分方位角叠加剖面同全方位叠加剖面相比信噪比没有降低,从而保证了断层解释的精度,同时不同的方位又具有其独特性和敏感性。
其中,宽方位角地震资料具有更丰富的地质信息,可以提供方位角道集、偏移距道集等多维地震数据,因此可进行部分角度叠加。许多研究表明,对于断裂等具有各向异性强的地层,地震振幅和频率随方位角的变化呈规律性变化。如果岩石介质中的各向异性是由一组定向的断裂引起,那么,根据地震波传播理论,当P波在各向异性介质中平行或垂直于断裂方向传播时具有不同的旅行速度,从而导致P波振幅响应的变化和旅行时的差异。因此同一条断裂在不同方位角叠加数据体上的成像精度存在差异,对于断裂的成像最为清晰的方位角可称之为地震敏感方位角。
在本申请一些实施方式中,所述的两两叠加例如可以对角线方向上的两组数据分组间两两叠加,并且所述两两叠加至少包含沿垂直于所述断裂发育展布方向的叠加,以及沿平行于所述断裂发育展布方向的叠加。所述共同的断裂走向方向是指参与叠加的两个分组所共同具有的断裂走向方向,例如数据分组1和数据分组2两两分方位叠加时,对应的共同的断裂走向方向为数据分组1和数据分组2所共同具有的断裂走向方向。
S103、获取每个新地震数据体的目的层沿层相干切片。
在本申请一些实施方式中,所述获取每个新地震数据体的目的层沿层相干切片例如可以包括:分别对所述多个新地震数据体进行倾角扫描,并利用倾角扫描结果进行构造导向滤波;其次将构造导向滤波后的多个新地震数据体进行相干处理,对应获得多个相干数据体;然后对每个相干数据体,沿目的层进行沿层相干提取,对应获得多个目的层沿层相干切片。如此,n/2个新地震数据体便可以得到n/2个目的层沿层相干切片。
在本申请一示例性实施方式中,以尼日尔某工区为例,将数据体A、数据体B、数据体C、数据体D分别进行倾角扫描,并利用倾角扫描结果进行构造导向滤波,然后对构造导向滤波后的数据体进行相干处理,并分别对4个相干数据体沿目的层E5提取沿层相干,得到4个沿层相干切片E5-Coh-A、E5-Coh-B、E5-Coh-C、E5-Coh-D。例如图4所示,为某一时刻的相干切片,从30-75度与120-165度两个相干切片的椭圆和箭头处的对比可看出不同方位角叠加数据对不同展布方向的断层的敏感性不同。
S104、根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体。
研究表明,不同方位角数据体对断裂识别的敏感程度有所差异,主要表现是断点清晰程度不一样,当数据体方位角垂直于断裂展布方向时,相干切片上断点干脆清晰,当数据体方位角平行于断裂展布方向时,相干体切片上断点模糊不清,其它方向数据识别清晰程度介于以上两者之间。因此,在本申请一些实施方式中,所述根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体可以是:确定其方位角与所述断裂发育展布方向垂直的目的层沿层相干切片,则该目的层沿层相干切片所对应新地震数据体为对断裂成像敏感的新地震数据体。这样,就可以从n/2个新地震数据体中确定对断裂成像敏感的m个新地震数据体。
S105、对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合。
在本申请一些实施方式中,在确定对断裂成像敏感的m个新地震数据体后,可以以对应的目的层沿层相干切片为依据在剖面上进行精细断裂解释,并进行平面与剖面断裂解释的吻合。在本申请一示例性实施方式中,以尼日尔某工区为例,从步骤S101可以得知本研究区区域断裂走向以北西--南东向展布为主,即在分方位数据体B(30-75度和210-255度叠加数据体)上北西向展布的断裂成像更清晰,在数据体D(120-165度和300-345度叠加数据体)上出现一组成像清晰的北东向展布的小断裂,因此根据沿层相干的结果挑选对断裂敏感的分方位数据B和D进行断裂解释,即分别在30-75度和210-255度叠加数据体以及120-165度和300-345度叠加数据体上进行剖面断裂解释,并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合。图5为沿主要目的层E5提取的不同方位角叠加相干切片与全叠加相干切片对比,箭头处的北东向小断裂在全方位叠加相干切片上没有显示,在120-165度分方位角叠加相干切片上有清晰显示;图6为沿垂直于北东向小断层的任意剖面显示,在垂直于该小断层的120-165度方位角剖面上断裂成像清晰于全方位叠加剖面。
S106、对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
在本申请一示例性实施方式中,以尼日尔某工区为例,以沿层相干切片E5-Coh-B为基础图片,把沿层相干切片E5-Coh-D上成像清晰的北东向小断裂进行图像提取,把提取出的断裂图像与沿层相干切片E5-Coh-B进行切片融合形成新的沿层相干切片E5-Coh-BD,使北东向断裂和北西向断裂能更清晰的出现在同一张融合后的目的层沿层相干切片上,使断裂的接触关系更加清晰,并依据新的沿层相干切片E5-Coh-BD上断裂的接触关系进行E5目的层断裂的平面组合,即把上一步剖面上解释的断裂进行分组,把同一条断裂在平面上绘制出来,然后根据切片上不同断裂之间的接触关系进行断裂的平面组合,从而获得目的层断裂。进而为本区块下一步圈闭的识别奠定好的断裂解释基础。图7为30-75度和120-165度两个方位角相干融合切片与全方位相干切片对比,全叠加相干切片上北西向展布的区域大断裂成像清晰,北东向小断层几乎没有显示,而在两个角度融合切片上北西向区域大断裂以及北东向小断层成像均比较清晰,同时两组断裂之间的接触关系也比较清晰;图8为主要目的层E1和E5顶的断裂平面组合对比,通过分方位融合解释的断层同全方位叠加数据解释的断层平面组合相比,断裂解释更加精细,断裂之间的接触关系更加合理。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
参考图9所示,本申请实施例的一种基于敏感方位地震数据的断裂识别装置可以包括:
断裂展布方向确定模块91,可以用于确定研究区的断裂发育展布方向;
地震数据分组叠加模块92,可以用于将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
目的层沿层相干切片获取模块93,可以用于获取每个新地震数据体的目的层沿层相干切片;
敏感方位数据确定模块94,可以用于根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
断裂解释及平剖面吻合模块95,可以用于对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
切片融合及断裂获取模块96,可以用于对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
参考图10所示,本申请实施例的另一种基于敏感方位地震数据的断裂识别装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
确定研究区的断裂发育展布方向;
将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
获取每个新地震数据体的目的层沿层相干切片;
根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
同样,本申请实施例的装置与上述实施例的方法对应,因此,有关于本申请的装置细节,请参见上述实施例的方法,在此不再赘述。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、装置的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (11)
1.一种基于敏感方位地震数据的断裂识别方法,其特征在于,包括:
确定研究区的断裂发育展布方向;
将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
获取每个新地震数据体的目的层沿层相干切片;
根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
2.如权利要求1所述的基于敏感方位地震数据的断裂识别方法,其特征在于,将所述研究区的地震数据按照方位角划分成偶数个数据分组,包括:
以所述地震数据正北方向为零度基准,按照360度方位将所述地震数据按照方位角等分成偶数个数据分组。
3.如权利要求1所述的基于敏感方位地震数据的断裂识别方法,其特征在于,所述两两叠加至少包含沿垂直于所述断裂发育展布方向的叠加,以及沿平行于所述断裂发育展布方向的叠加。
4.如权利要求1所述的基于敏感方位地震数据的断裂识别方法,其特征在于,所述获取每个新地震数据体的目的层沿层相干切片,包括:
分别对所述多个新地震数据体进行倾角扫描,并利用倾角扫描结果进行构造导向滤波;
将构造导向滤波后的多个新地震数据体进行相干处理,对应获得多个相干数据体;
对每个相干数据体,沿目的层进行沿层相干提取,对应获得多个目的层沿层相干切片。
5.如权利要求1所述的基于敏感方位地震数据的断裂识别方法,其特征在于,所述根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体,包括:
确定其方位角与所述断裂发育展布方向垂直的目的层沿层相干切片,则该目的层沿层相干切片所对应新地震数据体为对断裂成像敏感的新地震数据体。
6.一种基于敏感方位地震数据的断裂识别装置,其特征在于,包括:
断裂展布方向确定模块,用于确定研究区的断裂发育展布方向;
地震数据分组叠加模块,用于将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
目的层沿层相干切片获取模块,用于获取每个新地震数据体的目的层沿层相干切片;
敏感方位数据确定模块,用于根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
断裂解释及平剖面吻合模块,用于对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
切片融合及断裂获取模块,用于对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
7.如权利要求6所述的基于敏感方位地震数据的断裂识别装置,其特征在于,将所述研究区的地震数据按照方位角划分成偶数个数据分组,包括:
以所述地震数据正北方向为零度基准,按照360度方位将所述地震数据按照方位角等分成偶数个数据分组。
8.如权利要求6所述的基于敏感方位地震数据的断裂识别装置,其特征在于,所述两两叠加至少包含沿垂直于所述断裂发育展布方向的叠加,以及沿平行于所述断裂发育展布方向的叠加。
9.如权利要求6所述的基于敏感方位地震数据的断裂识别装置,其特征在于,所述获取每个新地震数据体的目的层沿层相干切片,包括:
分别对所述多个新地震数据体进行倾角扫描,并利用倾角扫描结果进行构造导向滤波;
将构造导向滤波后的多个新地震数据体进行相干处理,对应获得多个相干数据体;
对每个相干数据体,沿目的层进行沿层相干提取,对应获得多个目的层沿层相干切片。
10.如权利要求6所述的基于敏感方位地震数据的断裂识别装置,其特征在于,所述根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体,包括:
确定其方位角与所述断裂发育展布方向垂直的目的层沿层相干切片,则该目的层沿层相干切片所对应新地震数据体为对断裂成像敏感的新地震数据体。
11.一种基于敏感方位地震数据的断裂识别装置,包括存储器、处理器、以及存储在所述存储器上的计算机程序,其特征在于,所述计算机程序被所述处理器运行时执行如下步骤:
确定研究区的断裂发育展布方向;
将所述研究区的地震数据按照方位角划分成偶数个数据分组,并将所述数据分组沿相同走向方向两两叠加,获得多个新地震数据体;
获取每个新地震数据体的目的层沿层相干切片;
根据所述断裂发育展布方向及所述目的层沿层相干切片,确定对断裂成像敏感的新地震数据体;
对所述对断裂成像敏感的新地震数据体进行剖面断裂解释,获得解释结果;并根据所述解释结果对相应的目的层沿层相干切片进行平面与剖面的断裂解释吻合;
对断裂解释吻合后的目的层沿层相干切片进行切片融合,使不同敏感方位识别出的断裂展现于同一张融合后的目的层沿层相干切片上,并根据所述融合后的目的层沿层相干切片上断裂的接触关系进行目的层断裂的平面组合,获得目的层断裂。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710978827.3A CN107918151B (zh) | 2017-10-19 | 2017-10-19 | 基于敏感方位地震数据的断裂识别方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710978827.3A CN107918151B (zh) | 2017-10-19 | 2017-10-19 | 基于敏感方位地震数据的断裂识别方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107918151A true CN107918151A (zh) | 2018-04-17 |
CN107918151B CN107918151B (zh) | 2019-10-11 |
Family
ID=61894830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710978827.3A Active CN107918151B (zh) | 2017-10-19 | 2017-10-19 | 基于敏感方位地震数据的断裂识别方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107918151B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于构造导向滤波的古生界潜山储层地震预测方法 |
CN111751871A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油天然气集团有限公司 | 三维地震沿层相干切片处理方法、装置及电子设备 |
CN111830564A (zh) * | 2019-04-18 | 2020-10-27 | 中国石油天然气股份有限公司 | 识别断层的方法和装置 |
CN112394407A (zh) * | 2020-10-28 | 2021-02-23 | 中国石油天然气集团有限公司 | 陡坡扇体的预测方法及装置 |
CN113126156A (zh) * | 2021-04-13 | 2021-07-16 | 中海石油(中国)有限公司 | 拉东域高角度断裂提取方法、装置、存储介质及设备 |
CN114252914A (zh) * | 2020-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | 一种确定断裂系统分布的方法和装置 |
CN115407401A (zh) * | 2021-05-28 | 2022-11-29 | 中国石油化工股份有限公司 | 断裂带结构类型的地震表征方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120212229A1 (en) * | 2008-11-10 | 2012-08-23 | Cbg Corporation | Azimuthally sensitive resistivity logging tool |
CN104267431A (zh) * | 2014-10-13 | 2015-01-07 | 成都北方石油勘探开发技术有限公司 | 一种地震构造解释的小断层获取方法 |
CN105403916A (zh) * | 2015-10-27 | 2016-03-16 | 中国石油天然气集团公司 | 储层预测方法和装置 |
CN106556861A (zh) * | 2015-09-24 | 2017-04-05 | 中国石油化工股份有限公司 | 一种基于全方位地震资料的方位avo反演方法 |
-
2017
- 2017-10-19 CN CN201710978827.3A patent/CN107918151B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120212229A1 (en) * | 2008-11-10 | 2012-08-23 | Cbg Corporation | Azimuthally sensitive resistivity logging tool |
CN104267431A (zh) * | 2014-10-13 | 2015-01-07 | 成都北方石油勘探开发技术有限公司 | 一种地震构造解释的小断层获取方法 |
CN106556861A (zh) * | 2015-09-24 | 2017-04-05 | 中国石油化工股份有限公司 | 一种基于全方位地震资料的方位avo反演方法 |
CN105403916A (zh) * | 2015-10-27 | 2016-03-16 | 中国石油天然气集团公司 | 储层预测方法和装置 |
Non-Patent Citations (1)
Title |
---|
陈志刚 等: ""两宽一高地震资料的敏感方位油气检测技术在乍得潜山油藏描述的应用"", 《地球物理学进展》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 基于构造导向滤波的古生界潜山储层地震预测方法 |
CN109884708B (zh) * | 2019-03-20 | 2021-09-07 | 中国石油化工股份有限公司 | 基于构造导向滤波的古生界潜山储层地震预测方法 |
CN111751871A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油天然气集团有限公司 | 三维地震沿层相干切片处理方法、装置及电子设备 |
CN111751871B (zh) * | 2019-03-29 | 2023-02-07 | 中国石油天然气集团有限公司 | 三维地震沿层相干切片处理方法、装置及电子设备 |
CN111830564A (zh) * | 2019-04-18 | 2020-10-27 | 中国石油天然气股份有限公司 | 识别断层的方法和装置 |
CN114252914A (zh) * | 2020-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | 一种确定断裂系统分布的方法和装置 |
CN112394407A (zh) * | 2020-10-28 | 2021-02-23 | 中国石油天然气集团有限公司 | 陡坡扇体的预测方法及装置 |
CN113126156A (zh) * | 2021-04-13 | 2021-07-16 | 中海石油(中国)有限公司 | 拉东域高角度断裂提取方法、装置、存储介质及设备 |
CN113126156B (zh) * | 2021-04-13 | 2023-02-24 | 中海石油(中国)有限公司 | 拉东域高角度断裂提取方法、装置、存储介质及设备 |
CN115407401A (zh) * | 2021-05-28 | 2022-11-29 | 中国石油化工股份有限公司 | 断裂带结构类型的地震表征方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107918151B (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107918151A (zh) | 基于敏感方位地震数据的断裂识别方法及装置 | |
Nasuti et al. | STDR: A novel approach for enhancing and edge detection of potential field data | |
Grigoli et al. | Pick-and waveform-based techniques for real-time detection of induced seismicity | |
Burjánek et al. | Instrumental evidence of normal mode rock slope vibration | |
Thurber et al. | Three-dimensional compressional wavespeed model, earthquake relocations, and focal mechanisms for the Parkfield, California, region | |
Kuo et al. | Strong ground motion and pulse‐like velocity observations in the near‐fault region of the 2018 Mw 6.4 Hualien, Taiwan, earthquake | |
Bowles et al. | Discovering marine terraces using airborne LiDAR along the Mendocino-Sonoma coast, northern California | |
CN105954797B (zh) | 地震数据的断裂识别方法和装置 | |
CN107918156A (zh) | 检测海底节点采集地震数据极性的方法及装置 | |
CN106990438A (zh) | 复杂地表区域微测井点位及井深确定方法及装置 | |
CN107144877B (zh) | 异常构造识别方法及装置 | |
CN105866832A (zh) | 深层致密砂岩气藏的低级序断层的识别方法和装置 | |
Stolte et al. | An experimental topographic amplification study at Los Alamos National Laboratory using ambient vibrations | |
Nicollin et al. | Anisotropy of electrical conductivity of the excavation damaged zone in the Mont Terri Underground Rock Laboratory | |
CN107656312A (zh) | 基于分方位角叠加的河道砂体预测方法及装置 | |
Saltogianni et al. | A two‐fault model of the 2003 Leucas (Aegean arc) earthquake based on topological inversion of GPS data | |
Midzi et al. | 1-D velocity model for use by the SANSN in earthquake location | |
Smith et al. | Microseismic evidence for bookshelf faulting in western Montana | |
CN106842318A (zh) | 微观溶蚀孔隙地球物理二维表征确定方法和装置 | |
Arosio et al. | Seismic noise measurements on unstable rock blocks: the case of Bismantova rock cliff | |
Sheng et al. | Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array | |
Ikawa et al. | A case study of seismic anisotropy analysis on multi-azimuth OBC seismic data in an offshore carbonate oil field, UAE | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
Mostaccio et al. | Improving seismic surveillance at Mt. Etna volcano by probabilistic earthquake location in a 3D model | |
Sicking et al. | Forecasting reservoir performance by mapping seismic emissions |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |