CN109143354B - 一种地震波形特征分解的方法及装置 - Google Patents
一种地震波形特征分解的方法及装置 Download PDFInfo
- Publication number
- CN109143354B CN109143354B CN201810958255.7A CN201810958255A CN109143354B CN 109143354 B CN109143354 B CN 109143354B CN 201810958255 A CN201810958255 A CN 201810958255A CN 109143354 B CN109143354 B CN 109143354B
- Authority
- CN
- China
- Prior art keywords
- seismic
- data volume
- data
- radian
- horizon
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005070 sampling Methods 0.000 claims abstract description 54
- 230000008021 deposition Effects 0.000 claims abstract description 25
- 238000000354 decomposition reaction Methods 0.000 claims description 31
- 238000012545 processing Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient 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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
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
技术领域
本申请涉及地球物理勘探技术领域，特别涉及一种地震波形特征分解的方法及装置。
背景技术
随着勘探技术的发展，勘探的难度和深度也越来越大，对地震解释成果的精准度要求也越来越高。在地震数据采集时，由震源(炮点)激发，会产生一个地震子波，这个地震子波在地下介质中向下传播，当遇到一个反射界面时，地震子波会发生反射和透射，一部分地震子波的能量反射后改变方向，向上传播，到达地面接收点，并被布设在地面接收点的检波器接收。地下不同反射界面所发生的反射，依次会向地表传播，在地面的检波器会依次接收到这些地下反射的地震子波，构成了地震勘探仪器记录的地震数据。
目前，常见的地震波形特征分析技术可以分为两类：一是利用多道地震记录之间的波形特征变化，来识别储层、断层等；二是利用单道地震记录某一时间段内的波形，定性或定量地描述波形的特征，从而开展地震解释工作。地震的层位解释一般是根据同向轴，通过同一时期的波峰或者波谷来解释；沉积解释一般是根据同向轴和同向轴之间的关系，通过波形边界来解释。但是现有的方法都需要使用一段时间内的波形，波形中包含了峰值和边界在内，所以把波形的峰值和边界混在了一起，无法很好的区分开层位和沉积特征，从而导致地震解释工作具备较低的精度和效率。
发明内容
本申请实施例的目的是提供一种地震波形特征分解的方法及装置，以提高地震解释的精度和效率。
为达到上述目的，本申请实施例提供一种地震波形特征分解的方法，包括：
获取目标区域的地震数据，所述地震数据中包括至少一个采样点；
将所述地震数据转换为弧度域的地震数据体，所述弧度域的地震数据体中包括至少一个弧度值，所述弧度值由所述采样点的斜率确定；
将所述弧度域的地震数据体中的弧度值与指定阈值进行对比，并根据对比结果，将所述弧度域的地震数据体分解为层位数据体和沉积数据体；
基于所述层位数据体进行层位解释，并基于所述沉积数据体进行沉积解释。
本申请实施例的地震波形特征分解方法，将所述地震数据转换为弧度域的地震数据体，包括：
针对所述地震数据中除最后一个采样点之外的其他采样点，计算所述其他采样点的斜率，并对计算得到的所述斜率进行反正切运算，以得到所述其他采样点对应的弧度值。
本申请实施例的地震波形特征分解方法，所述方法还包括：
针对所述最后一个采样点进行反正切运算，以得到所述最后一个采样点对应的弧度值。
本申请实施例的地震波形特征分解方法，按照下面公式计算所述采样点的斜率L：
L＝(zi+1-zi)/(ti+1-ti)
其中zi表示第i个采样点的数据，ti表示第i个采样点的数据是在t时刻记录的。
本申请实施例的地震波形特征分解方法，将所述弧度域的地震数据体分解为层位数据体和沉积数据体，包括：
将绝对值小于所述指定阈值的弧度值确定为所述层位数据体，并将绝对值大于或者等于所述指定阈值的弧度值确定为所述沉积数据体。
本申请实施例的地震波形特征分解方法，所述方法还包括：基于所述层位数据体绘制地震层位数据体剖面图，基于所述沉积数据体绘制地震沉积数据体剖面图。
本申请实施例还提供一种地震波形特征分解的装置，该装置包括：
数据获取模块，用于获取目标区域的地震数据；
数据处理模块，用于计算采样点的斜率，再进行反正切运算，得到弧度域的地震数据体；
数据分解模块，用于将弧度域的地震数据体分解为层位数据体和沉积数据体；
地震解释模块，用于进行层位解释和沉积解释。
本申请实施例的地震波形特征分解的装置，所述数据处理模块，包括：
斜率计算单元，用于计算地震数据中采样点的斜率；
弧度计算单元，用于对采样点和采样点的斜率进行反正切运算。
本申请实施例的地震波形特征分解的装置，所述数据分解模块，包括：
筛选单元，用于将弧度值的绝对值与指定阈值进行比较；
图像处理单元，用于将层位数据体和沉积数据体分别绘制地震剖面图。
本申请实施例的地震波形特征分解的装置，所述地震解释模块，包括：
预处理单元，用于将地震子波和反射系数序列进行卷积，得到合成地震记录；
井震标定单元，用于将得到的合成地震记录和井旁地震进行对比；
解释单元，用于进行层位解释和沉积解释。
由上可见，本申请提供的技术方案进行了一种波形变化算法研究，将地震数据转化为对应的弧度值，通过设立阈值，将弧度域的地震数据体分解成两部分，分别为层位数据体和沉积数据体，从而把地震波形的峰值位置从地震波形中准确地拾取出来，从而避免地震波形中的峰值位置与边界位置相重叠。基于层位数据体和沉积数据体分别绘制的地震剖面图能够更好地反映出地层的层位特征和沉积特征。因此，本申请能够为地震解释工作提供更有针对性的数据体，提高了地震解释的精度和效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一种地震波形特征分解的流程图；
图2是本申请实施例中原始地震剖面的示意图；
图3是本申请实施例中计算采样点斜率的示意图；
图4是本申请实施例中地震波形分解的示意图；
图5是本申请实施例中波形分解后的地震层位数据体剖面的示意图；
图6是本申请实施例中波形分解后的地震沉积数据体剖面的示意图；
图7是本申请实施例中在波形分解后的地震层位数据体剖面上进行层位解释的示意图；
图8是本申请实施例中在波形分解后的地震沉积数据体剖面上进行沉积解释的示意图；
图9是本申请实施例中地震波形特征分解的装置的模块结构示意图。
具体实施方式
本申请实施例提供一种地震波形特征分解的方法及装置。
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本申请一部分实施方式，而不是全部的实施方式。基于本申请中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。
参考图1所示，本申请实施方式的地震波形特征分解的方法可以包括如下步骤：
S101：获取目标区域的地震数据，所述地震数据中包括至少一个采样点。
在本实施方式中，所述目标区域的地震数据是通过人工地面激发，产生一个地震子波，地震子波在地下传播，遇到介质性质不同的岩层分界面，地震子波发生反射与透射，在地面用检波器(地震波接受装置)进行地震信号的接收，并由地震勘探仪器采集和处理后的在时间域反映地下地层特征的数据。上述目标区域的地震数据具体可以是目标区域的三维地震数据体。上述三维地震数据体具体是一种通过三维地震勘探获取的地震数据。其中，上述三维地震勘探区别于二维地震勘探，通过三维地震勘探不仅能获得一张张地震剖面图，还能获得一个三维空间的数据体，上述三维地震数据体的信息点的密度通常可以达到12.5米×12.5米(即在12.5米×12.5米的面积内便采集一个数据)，相比于二维测线信息点的密度(通常最高为1千米×1千米)，三维地震勘探获得信息量丰富，地震剖面分辨率高。从上述三维地震数据体中提取关键数据，包括线数、道数，所述线数为地震数据体在X方向(东西方向或近东西方向)上记录的地震道个数，所述道数为地震数据体在Y方向(南北方向或近南北方向)上记录的地震道个数。当然，需要说明的是，上述所列举的三维地震资料只是为了更好地说明本申请实施方式。具体实施时，也可以根据目标区域的具体情况和精度要求，选择其他类型的地震数据作为目标区域的地震数据。对此，本申请不作限定。
地震数据有多种保存格式，比较常用的是SEG-Y格式，它是由SEG(Society ofExploration Geophysicists，勘探地球物理学家协会)提出的标准磁带数据格式之一，可以由三部分组成：400字节的二进制文件头、每个地震道的240字节道头和每个地震道的地震数据，它的数据都可以按地震道的顺序一个个存放。从每个地震道的240字节的道头中提取这个地震道的采样点数、采样间隔(以微秒表示)，并提取每一个地震道的数据，上述地震道的数据是对地震信号按采样间隔进行取样，再把这一系列采样点的离散振幅值记录下来。
参考图2所示，是通过在上述目标区域的三维地震数据体中选择其中一条测线，获得的地震剖面图。得到地震剖面图之前，先将上述测线中的每个地震道的数据在一个平面内以波形的形式表现出来，其中横坐标代表地震波振幅、纵坐标代表时间，地震数据的大小代表振幅的强弱，正负代表振动的方向，再将上面得到的地震道的波形按照检波点在地面的顺序依次排列。在一系列地震道中，同一反射波的相同相位在相邻地震道上的到达时间是相近的，每道记录的振幅也是相似的，因此相邻地震道之间波形相似，波形的峰值比较靠近，在地震剖面上相互叠套成串，一连串的峰值组成一条线，形成同相轴。通常为了更好的表现同相轴，地震剖面图采用变面积显示的方式，包括正极性填充、负极性填充、双极性填充等方式。本申请实施方式采用正极性填充的方式进行展示，具体的方法为将波形振幅为正(即正值采样点)的波形线与零值位置线围成的多个不规则多边形填充自定义的颜色(一般为黑色)，而在波形振幅为负(即负值采样点)只作曲线连接，将所有的波形都进行填充后得到图2所示的地震剖面图。当然，需要说明的是上述选择变面积显示只是为了更好地说明本申请的实施方式。具体实施时，也可以选择变密度显示等其他方式。对此，本申请不作限定。
S102：将所述地震数据转换为弧度域的地震数据体，所述弧度域的地震数据体中包括至少一个弧度值，所述弧度值由所述采样点的斜率确定。
参考图3所示是一个地震道的波形，设相邻的两个采样点为(ti,zi)、(ti+1,zi+1)，则第i个采样点的斜率L的计算公式为L＝(zi+1-zi)/(ti+1-ti)，其中，zi表示第i个采样点处的振幅，ti表示第i个采样点的数据是在t时刻被记录的。用此方法计算出地震道中第一个采样点到倒数第二个采样点的斜率，然后对采样点的斜率进行反正切运算求出弧度值。由于最后一个采样点的后面没有采样点，无法用上述方法计算斜率，所以对地震道的最后一个采样点直接求反正切运算求出弧度值。经过上述步骤，将这个地震道中所有采样点的振幅值都转换为了对应的弧度值。对三维地震数据体中其余的地震道都重复上述操作，从而得到了弧度域的地震数据体。
S103：将所述弧度域的地震数据体中的弧度值与指定阈值进行对比，并根据对比结果，将所述弧度域的地震数据体分解为层位数据体和沉积数据体。
在地震解释中，层位解释的基础就是寻找同一时期的同相轴，沉积解释的基础是同相轴之间的关系。参考图3所示，由于在原始的地震数据体中，每个地震道的波形都是连续的，所以较难准确地拾取同一时期的波峰或波谷。但是在实际应用中，位于波峰和波谷附近位置的采样点，通过上面介绍的计算采样点斜率的方法，它们的斜率值都在零附近，经过反正切运算后得到的弧度值也在零附近，所以通过选择一个合适的阈值，可以将上面得到的弧度域的地震数据体根据弧度值的绝对值进行分解，从而将峰值附近的采样点准确地挑选出来。所述阈值可以经过多次试验进行确定，也可以通过对历史地震数据进行分析后得到。例如，阈值的选择可以为绝对值0.1。当然，需要说明的是上述所列举的阈值的数值只是为了更好地说明本申请的实施方式。具体实施时，也可以根据具体的精度要求，选择其他数值作为上述阈值。对此，本申请不作限定。
确定了阈值后，将弧度域的地震数据体中绝对值小于阈值的弧度值提取出来，绝对值不小于阈值的弧度值舍弃，组成了弧度域的层位数据体；将绝对值不小于阈值的弧度值提取出来，绝对值小于阈值的弧度值舍弃，组成了弧度域的沉积数据体。从而将弧度域的地震数据体进行了分解。
参考图4所示，图中用虚线绘制的波形是基于原始的地震数据体，用实线绘制的波形是基于这个地震道的弧度域的地震数据体。上面介绍的利用阈值将弧度域的地震数据体进行分解的结果，表现在平面上的地震波形中，可以将图4中的实线波形进行分解，其中，在图4中用圆圈标记出来的部分是弧度域的层位数据体，实线波形上除此之外的部分是弧度域的沉积数据体，从而实现了将一个地震道的波形进行了分解。从图4中可以直观的看出，由于选择的阈值为0.1，所以用圆圈标记的部分位于实线波形与坐标轴的交点附近。从图4中还可以清楚的看出，本申请提出的地震波形特征分解方法并不是将原始的地震波形进行分段，而是在基于弧度域的地震数据体得到的波形上面进行采样点的选取，分别提取出在时间轴上对应原始波形不同位置的采样点。从图4可知，用圆圈标记的部分在时间轴的投影与原始波形中的峰值部分在时间轴的投影基本重合，准确地拾取了原始波形中的峰值位置，从而提取到了地震波的关键信息，便于更加准确快速地识别出同相轴。从图4还可以看出，由于弧度值有固定的值域范围，即-π/2～π/2，所以基于弧度域数据体绘制的波形，都分布在-π/2～π/2之间。
选择与图2相同的测线，将此测线中的地震道波形都按上述阈值和波形分解方法进行分解。选取弧度域的层位数据体进行剖面图的绘制，可以得到地震层位数据体剖面图，如图5所示；选取弧度域的沉积数据体进行剖面图的绘制，可以得到地震沉积数据体剖面图，如图6所示。
S104：基于所述层位数据体进行层位解释，并基于所述沉积数据体进行沉积解释。
本申请实施例中，为了进行地震的层位解释，首先需要进行层位标定。为此，还要获取目标区域的测井数据，所述测井数据应该至少包括声波测井和密度测井，以便于后面求取波阻抗。从声波测井曲线可获取速度，从密度测井曲线可获取密度，根据波阻抗为速度和密度的乘积，从而可得到波阻抗曲线。由于时间剖面的纵坐标是双程旅行时，而测井数据是以铅锤深度表示的，所以两者还需要进行时深转换，媒介就是地震波的传播速度，将波阻抗曲线进行时深转换，再由波阻抗曲线计算反射系数曲线。
一般认为，地震震源激发时所产生的地震波仅是一个延续时间极短的尖脉冲，随着尖脉冲在粘弹性介质中传播，尖脉冲的高频成分很快衰减，波形随之增长，便形成了地震子波，对上述反射系数曲线和地震子波进行卷积，得到合成地震记录。
根据所述合成地震记录和所述地震数据，还需要进行井震标定，目的是进一步明确哪个地震波形对应哪个地质层位，具体实施时，可以包括：将上述合成地震记录与井旁地震(即井旁地震数据)进行比较，并根据比较结果进行调整，以使得合成地震记录与井旁地震达到匹配。其中，上述合成地震记录与井旁地震达到匹配具体可以包括：合成地震记录的波峰位置与井旁地震波峰位置重合，合成地震记录的波谷位置与井旁地震的波谷位置重合。
经过上述步骤，完成了层位标定的过程，然后就可以进行层位解释，所述层位解释是对在地震记录上相位相同且为同一时期的地层进行追踪解释。经过波形分解后的地震层位数据体剖面，层位界限更加清晰，在剖面上表现为一条条细线，更有利于进行同相轴自动追踪和层位解释。并且由于常规数据体是连续变化的，在常规数据体上根据数据的波峰、波谷进行层位解释，往往不能准确拾取同一时期的波峰或波谷值，就会产生跳点，而波形分解后，层位数据体只有波峰或波谷处有值，其他位置可以认为是无值的，这样就能更准确的拾取同一时期的波峰或波谷值，而不会产生跳点。在图5中就已经知道了哪些线对应T1、T2，沿着这些线进行解释就得到T1、T2，参考图7所示。
层位解释完后，明确了地层的空间排布，再进行层与层之间的沉积特征解释，所述沉积解释是指对地层的接触关系和地震的反射特征进行解释，其中的接触关系是指新老地层或岩石在空间上的相互叠置状态，在图8中的最大湖泛面，前积扇，前积特征明显，更易解释。
图9是本申请地震特征波形分解装置的模块结构示意图。如图9所示，所述地震特征波形分解装置可以包括：数据获取模块100、数据处理模块200、数据分解模块300、地震解释模块400，下面对该装置进行具体说明。
数据获取模块，用于获取目标区域的地震数据；
数据处理模块，用于计算采样点的斜率，再进行反正切运算，得到弧度域的地震数据体；
数据分解模块，用于将弧度域的地震数据体分解为层位数据体和沉积数据体；
地震解释模块，用于进行层位解释和沉积解释。
本申请实施方式的地震波形特征分解的装置，所述数据处理模块，包括：
斜率计算单元，用于计算地震数据中采样点的斜率；
弧度计算单元，用于对采样点和采样点的斜率进行反正切运算。
本申请实施方式的地震波形特征分解的装置，所述数据分解模块，包括：
筛选单元，用于将弧度值的绝对值与指定阈值进行比较；
图像处理单元，用于将层位数据体和沉积数据体分别绘制地震剖面图。
本申请实施方式的地震波形特征分解的装置，所述地震解释模块，包括：
预处理单元，用于将地震子波和反射系数序列进行卷积，得到合成地震记录；
井震标定单元，用于将得到的合成地震记录和井旁地震进行对比；
解释单元，用于进行层位解释和沉积解释。
需要说明的是，上述的地震特征波形分解装置中的各个模块和单元所实现的具体功能，可参见上述的方法实施例，并可以实现上述方法实施例的技术效果，在此便不再赘述。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (9)
1.一种地震波形特征分解的方法，其特征在于，包括：
获取目标区域的地震数据，所述地震数据中包括至少一个采样点；
将所述地震数据转换为弧度域的地震数据体，所述弧度域的地震数据体中包括至少一个弧度值，所述弧度值由所述采样点的斜率确定；
将所述弧度域的地震数据体中的弧度值与指定阈值进行对比，并根据对比结果，将所述弧度域的地震数据体分解为层位数据体和沉积数据体；其中，将绝对值小于所述指定阈值的弧度值确定为所述层位数据体，并将绝对值大于或者等于所述指定阈值的弧度值确定为所述沉积数据体；
基于所述层位数据体进行层位解释，并基于所述沉积数据体进行沉积解释。
2.根据权利要求1所述的方法，其特征在于，将所述地震数据转换为弧度域的地震数据体，包括：
针对所述地震数据中除最后一个采样点之外的其他采样点，计算所述其他采样点的斜率，并对计算得到的所述斜率进行反正切运算，以得到所述其他采样点对应的弧度值。
3.根据权利要求2所述的方法，其特征在于，所述方法还包括：
针对所述最后一个采样点进行反正切运算，以得到所述最后一个采样点对应的弧度值。
4.根据权利要求1所述的方法，其特征在于，按照下面公式计算所述采样点的斜率L：
L＝(zi+1-zi)/(ti+1-ti)
其中zi表示第i个采样点的数据，ti表示第i个采样点的数据是在t时刻记录的。
5.根据权利要求1所述的方法，其特征在于，所述方法还包括：基于所述层位数据体绘制地震层位数据体剖面图，基于所述沉积数据体绘制地震沉积数据体剖面图。
6.一种地震波形特征分解的装置，其特征在于，该装置包括：
数据获取模块，用于获取目标区域的地震数据；
数据处理模块，用于计算采样点的斜率，再进行反正切运算，得到弧度域的地震数据体；
数据分解模块，用于将弧度域的地震数据体分解为层位数据体和沉积数据体；其中，将绝对值小于指定阈值的弧度值确定为所述层位数据体，并将绝对值大于或者等于所述指定阈值的弧度值确定为所述沉积数据体；
地震解释模块，用于进行层位解释和沉积解释。
7.根据权利要求6所述的装置，其特征在于，所述数据处理模块，包括：
斜率计算单元，用于计算地震数据中采样点的斜率；
弧度计算单元，用于对采样点和采样点的斜率进行反正切运算。
8.根据权利要求6所述的装置，其特征在于，所述数据分解模块，包括：
筛选单元，用于将弧度值的绝对值与指定阈值进行比较；
图像处理单元，用于将层位数据体和沉积数据体分别绘制地震剖面图。
9.根据权利要求6所述的装置，其特征在于，所述地震解释模块，包括：
预处理单元，用于将地震子波和反射系数序列进行卷积，得到合成地震记录；
井震标定单元，用于将得到的合成地震记录和井旁地震进行对比；
解释单元，用于进行层位解释和沉积解释。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810958255.7A CN109143354B (zh) | 2018-08-22 | 2018-08-22 | 一种地震波形特征分解的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810958255.7A CN109143354B (zh) | 2018-08-22 | 2018-08-22 | 一种地震波形特征分解的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109143354A CN109143354A (zh) | 2019-01-04 |
| CN109143354B true CN109143354B (zh) | 2020-03-10 |
Family
ID=64791055
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810958255.7A Active CN109143354B (zh) | 2018-08-22 | 2018-08-22 | 一种地震波形特征分解的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143354B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101676745A (zh) * | 2008-09-18 | 2010-03-24 | 核工业北京地质研究院 | 一种dct域内的重磁场分离方法 |
| CN103675897A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种地震绕射波分离成像方法 |
| CN106133556A (zh) * | 2014-01-10 | 2016-11-16 | 斯塔特伊石油公司 | 确定波场的分量 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9442209B2 (en) * | 2012-07-10 | 2016-09-13 | Pgs Geophysical As | Methods and systems for reconstruction of low frequency particle velocity wavefields and deghosting of seismic streamer data |
| US9594134B2 (en) * | 2014-01-22 | 2017-03-14 | Wisconsin Alumni Research Foundation | System and method for fully phase-encoded magnetic resonance imaging using multiband radio frequency excitation |
-
2018
- 2018-08-22 CN CN201810958255.7A patent/CN109143354B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101676745A (zh) * | 2008-09-18 | 2010-03-24 | 核工业北京地质研究院 | 一种dct域内的重磁场分离方法 |
| CN103675897A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种地震绕射波分离成像方法 |
| CN106133556A (zh) * | 2014-01-10 | 2016-11-16 | 斯塔特伊石油公司 | 确定波场的分量 |
Non-Patent Citations (1)
| Title |
|---|
| "地震数据的稀疏高斯束分解方法";刘鹏 等;《地球物理学报》;20131130;第56卷(第11期);第3887-3895页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109143354A (zh) | 2019-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA3076280C (en) | Method and apparatus for determining thin interlayers | |
| CN113759424B (zh) | 基于频谱分解和机器学习的岩溶储层充填分析方法和系统 | |
| US11402529B2 (en) | Identifying geologic features in a subterranean formation using seismic diffraction and refraction imaging | |
| US9075163B2 (en) | Interferometric seismic data processing | |
| CN102395902B (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
| CA2964893C (en) | Structure tensor constrained tomographic velocity analysis | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| CN108196304B (zh) | 一种多次波压制方法及装置 | |
| CN108363096B (zh) | 一种水陆检波器地震数据上下行波场分离方法及装置 | |
| CN106054251A (zh) | 一种初至波拾取方法及装置 | |
| CN107193045A (zh) | 一种地震数据处理方法及装置 | |
| CN103163567B (zh) | 基于能量的气枪子波气泡比分析方法及装置 | |
| US6985404B2 (en) | Method and apparatus of performing amplitude variation with offset analysis | |
| CN108897040B (zh) | 一种确定陆检与海底之间的耦合性的方法及装置 | |
| CN104199083A (zh) | 一种确定宽线地震观测系统的方法和装置 | |
| CN104502996A (zh) | 一种密度曲线校正方法及系统 | |
| CN104375178A (zh) | 一种碳酸盐岩缝洞储层预测方法和装置 | |
| EP3978961B1 (en) | System and method for quantitative seismic integration modeling workflow | |
| CN106772599B (zh) | 一种计算地层横波速度的方法及装置 | |
| CN109143354B (zh) | 一种地震波形特征分解的方法及装置 | |
| CN108035708A (zh) | 一种去除地层界面反射波的方法及装置 | |
| CN1308703C (zh) | 用于生成相似性图和计算反射体倾角的地震探测方法及设备 | |
| CN111965705B (zh) | 地震单炮记录地质层位的标定方法、装置、设备及介质 | |
| CN111025393B (zh) | 针对含薄煤层地层的储层预测方法、装置、设备及介质 | |
| CN114114421A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |