CN107703541B - 地层倾角的确定方法和装置 - Google Patents
地层倾角的确定方法和装置 Download PDFInfo
- Publication number
- CN107703541B CN107703541B CN201710740537.5A CN201710740537A CN107703541B CN 107703541 B CN107703541 B CN 107703541B CN 201710740537 A CN201710740537 A CN 201710740537A CN 107703541 B CN107703541 B CN 107703541B
- Authority
- CN
- China
- Prior art keywords
- preset direction
- measured
- seismic data
- preset
- imaginary
- 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 53
- 238000010606 normalization Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 238000003325 tomography Methods 0.000 description 8
- 238000012512 characterization method Methods 0.000 description 5
- 238000011002 quantification Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000009466 transformation 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
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
技术领域
本发明涉及地质勘探技术领域，特别涉及一种地层倾角的确定方法和装置。
背景技术
地质勘探目标往往存在地层倾角，因而在以数据驱动为特点的地震资料处理、解释和反演的过程中，如果不考虑勘探目标实际存在的地层倾角，地质资料解释的精度可能会受到较大影响，甚至会导致错误的结果。因此，如何计算得到稳定而精确的地层倾角，进而沿地层倾角方向进行分析以得到更加准确的地质目标解释结果，具有重要的意义。
目前，可以通过局部倾斜叠加方法、平面波分解滤波器法、基于相干度量的扫描搜索法以及梯度结构张量法来确定地层倾角。然而，上述这些方法大多基于局部时窗分析来确定地层倾角，并没有将断层、不整合等地层边缘区域与地震同相轴连续性良好的区域区别对待，因而容易将上述地层边缘区域的倾角估计为断层两侧同相轴的平均倾角，最终导致计算精度较低。
发明内容
本发明提供了一种地层倾角的确定方法和装置，以解决现有技术中容易将地层边缘区域的倾角估计为断层两侧同相轴的平均倾角的问题。
本发明实施例提供了一种地层倾角的确定方法，可以包括：计算待测地震数据在多个预设方向中各个预设方向的总能量；确定多个总能量中的最大值，并确定所述最大值对应的方向；计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。
在一个实施例中，确定多个总能量中的最大值，可以包括：对所述多个预设方向中各个预设方向的总能量进行去噪处理；根据去噪处理后的各个预设方向的总能量，确定多个总能量中的最大值。
在一个实施例中，可以按照以下公式对所述多个预设方向中各个预设方向的总能量进行去噪处理：
其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，1≤io≤No，No表示所述预设方向的总数，n＝(n1,n2)，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量，T表示预设的去噪阈值。
在一个实施例中，所述预设的去噪阈值可以为TotalE(io,n)×10-6。
在一个实施例中，计算所述待测地震数据在多个预设方向中各个预设方向的总能量，可以包括：根据所述待测地震数据的曲波变换系数，计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量；计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和；根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。
在一个实施例中，根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量，可以包括：对所述实部和进行归一化处理，得到归一化后的实部和；对所述实部和进行归一化处理，得到归一化后的实部和；根据所述实部分量、所述虚部分量、所述归一化后的实部和以及所述归一化后的虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。
在一个实施例中，可以按照以下公式计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量：
其中，RCoef(is,io,k).r＝Coef(is,io,k).r，RCoef(is,io,k).i＝0，ICoef(is,io,k).r＝0，ICoef(is,io,k).i＝Coef(is,io,k).i，
其中，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，Seis2D(n)表示所述待测地震数据，n＝(n1,n2)，0≤n1≤Nt,0≤n2≤Nx表示所述待测地震数据的下标，Nt表示所述待测地震数据在时间方向的采样点数，Nx表示所述待测地震数据在空间方向的采样点数，Coef(is,io,k)表示所述曲波变换系数，is表示所述预设尺度，1≤is≤Ns，Ns表示所述预设尺度的最大值，io表示所述预设方向，1≤io≤No，No表示所述预设方向的总数，k＝(k1,k2)表示所述曲波变换系数的下标，表示的共轭复数，表示在预设尺度为is、预设方向为io、下标为k时的离散曲波，表示ψis,io,k(n1,n2)的共轭复数，ψis,io,k(n1,n2)表示在预设尺度为is、预设方向为io、下标为k时的离散逆曲波，RCoef(is,io,k)表示置复系数Coef(is,io,k)的虚部为零，ICoef(is,io,k)表示置复系数Coef(is,io,k)的实部为零。
在一个实施例中，可以按照以下公式计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量：
其中，TmpA(is,io,n)＝RS(is,io,n)·MeanSRS(io,n)+IS(is,io,n)·MeanSIS(io,n)，
TmpB(is,io,n)＝RS(is,io,n)·MeanSIS(io,n)-IS(is,io,n)·MeanSRS(io,n)；
其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，MeanSRS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的实部和，MeanSIS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的虚部和。
在一个实施例中，可以按照以下公式计算所述待测地震数据相应的地层倾角：
Dip(n)＝IO(n)·360No+180No
其中，Dip(n)表示所述待测地震数据在下标为n时相应的地层倾角，IO(n)表示预设方向为io、下标为n时的总能量TotalE(io,n)中的最大值对应的方向编号，No表示所述预设方向的总数。
在一个实施例中，还可以包括：在计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角之后，利用所述地层倾角进行地质勘探。
本发明实施例还提供了一种地层倾角的确定装置，可以包括：总能量计算模块，可以用于计算待测地震数据在多个预设方向中各个预设方向的总能量；方向确定模块，可以用于确定多个总能量中的最大值，并确定所述最大值对应的方向；地层倾角确定模块，可以用于计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。
在本发明实施例中，计算待测地震数据在多个预设方向中各个预设方向的总能量，并基于多个总能量中的最大值对应的方向，计算待测地震数据相应的地层倾角。由于在曲波域对待测地震数据进行处理时能够实现地震数据的稀疏表示，即，能够从多个尺度上对多个方向的地震数据反射特征的细节进行定量化表征，从而解决了现有技术中没有将断层、不整合等地层边缘区域与地震同相轴延展性良好的区域区别对待，所导致的计算精度下降的缺陷。根据上述定量化表征结果，确定待测地震数据中反射同相轴的地层倾角信息，从而在保证计算效率的同时保证较高准确度的局部倾角计算结果，从而提高了以数据驱动为特点地震资料处理、解释和反演方法的性能。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种地层倾角的确定方法流程图；
图2是本申请提供的待测叠后地震剖面的示意图；
图3是本申请提供的采用本申请所提出的地层倾角的确定方法确定某地区的地层倾角的示意图；
图4是本申请提供的利用确定的地层倾角进行平面波解构的结果示意图；
图5是本申请提供的一种地层倾角的确定装置的一种结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
考虑到现有技术中确定地层倾角时，由于没有将断层、不整合等地层边缘区域与地震同相轴延展性良好的区域区别对待，因而容易将上述地层边缘区域的倾角估计为断层两侧同相轴的平均倾角，最终导致计算精度较低的缺陷，发明人提出了将待测地震数据转换至能够从多个尺度上对多个方向的地震数据反射特征进行定量化表征的曲波域，通过选取待测地震数据在曲波域多个预设方向中各个预设方向的总能量中的最大值所对应的方向，计算得到地层倾角。基于此，提出了一种地层倾角的确定方法，如图1所示，可以包括以下步骤：
S101：计算待测地震数据在多个预设方向中各个预设方向的总能量。
在本实施例中，曲波域中的曲波可以由三个参量确定，分别为：尺度is、方向io以及空间位置k，k＝(k1,k2)。曲波具有可变的长度和宽度，并且具有较高的各向异性。可以利用曲波变换将待测地震数据变换至能够对待测地震数据进行稀疏表示的曲波域，并计算待测地震数据在多个预设方向中各个预设方向的总能量。具体的，可以包括以下步骤：
S1-1：根据待测地震数据的曲波变换系数，计算得到待测地震数据在多个预设尺度下多个预设方向下的实部分量以及虚部分量。
步骤1：可以根据曲波变换将待测地震数据变换到曲波域，得到曲波变换系数。其中，上述待测地震数据可以为叠后地震记录或者叠前地震记录。
在本申请的一个实施例中，可以按照以下公式计算曲波变换系数Coef：
其中，Seis2D(n)为输入的待测地震数据，n＝(n1,n2)，0≤n1≤Nt,0≤n2≤Nx为输入的待测地震数据的下标，Nt表示待测地震数据在时间方向的采样点数，Nx表示待测地震数据在空间方向的采样点数，Coef(is,io,k)表示所述待测地震数据在曲波域的预设尺度为is、预设方向为io、下标为k时的曲波变换系数，1≤is≤Ns，Ns表示预设尺度的最大值，io表示预设方向，k＝(k1,k2)为所述曲波变换系数的下标，为的共轭复数，其中为在预设尺度为is、预设方向为io、下标为k的离散曲波。
在本申请的一个实施例中，Ns可以为用户基于数据处理设备的内存等情况所设定的最大尺度数。例如，Ns可以为范围在3至6的整数值。当所设定的Ns越大，所确定的曲波变换精度越高，所需要的存储空间(即，内存)也就越大。
在本申请的一个实施例中，No可以为用户给定的最大方向数。例如，No可以为8、16或者32的整数值。当所设定的No越大，所确定的曲波变换精度越高，所需要的存储空间(即，内存)也就越大。
步骤2：根据上述曲波变换系数，通过曲波反变换获得待测地震数据各个预设尺度和多个预设方向下的实部分量以及虚部分量。
具体的，可以按照以下公式计算得到待测地震数据在多个预设尺度多个预设方向的实部分量以及虚部分量：
其中，RCoef(is,io,k).r＝Coef(is,io,k).r，RCoef(is,io,k).i＝0，ICoef(is,io,k).r＝0，ICoef(is,io,k).i＝Coef(is,io,k).i；
其中，RS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，表示ψis,io,k(n1,n2)的共轭复数，ψis,io,k(n1,n2)表示在预设尺度为is、预设方向为io、下标为k时的离散逆曲波，RCoef(is,io,k)表示置复系数Coef(is,io,k)的虚部为零，ICoef(is,io,k)表示置复系数Coef(is,io,k)的实部为零。其中，置复系数Coef(is,io,k)的虚部为零表示的是将复系数Coef(is,io,k)的虚部设置为0，置复系数Coef(is,io,k)的实部为零表示的是将复系数Coef(is,io,k)的实部设置为0。
S1-2：计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；
S1-3：计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和。
在本实施例中，可以根据待测地震数据在各个预设尺度和各个预设方向上的实部分量及虚部分量，计算各个预设方向上的实部和以及虚部和。
在本申请的一个实施例中，可以按照以下公式得到多个预设方向中各个预设方向的实部和：
其中，SRS(io,k)表示所述待测地震数据在预设方向为io、下标为n时的实部和。
在本申请的一个实施例中，可以按照以下公式得到多个预设方向中各个预设方向的虚部和：
其中，SIS(io,k)表示所述待测地震数据在预设方向为io、下标为n时的虚部和。
进一步地，在本申请的一个实施例中，在得到多个预设方向中各个预设方向的虚部和之后，还可以对实部和、虚部和进行归一化处理，得到归一化后的实部和、虚部和。
具体的，可以按照以下公式分别计算得到归一化后的实部和、虚部和：
其中，MeanSRS(io,n)表示归一化后的在预设方向为io、下标为n时的实部和，MeanSIS(io,n)表示归一化后的在预设方向为io、下标为n时的虚部和。
S1-4：根据实部分量、虚部分量、实部和以及虚部和，计算得到待测地震数据在多个预设方向中各个预设方向的总能量。
可以对实部和进行归一化处理，得到归一化后的实部和；对实部和进行归一化处理，得到归一化后的实部和；根据实部分量、虚部分量、归一化后的实部和以及归一化后的虚部和，计算得到待测地震数据在多个预设方向中各个预设方向的总能量。
其中，TmpA(is,io,n)＝RS(is,io,n)·MeanSRS(io,n)+IS(is,io,n)·MeanSIS(io,n)，
TmpB(is,io,n)＝RS(is,io,n)·MeanSIS(io,n)-IS(is,io,n)·MeanSRS(io,n)；
其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，RS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,k)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，MeanSRS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的实部和，MeanSIS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的虚部和。
在本申请的一个实施例中，在计算得到待测地震数据在多个预设方向中各个预设方向的总能量之后，可以对多个预设方向中各个预设方向的总能量进行去噪处理。
具体的，可以按照以下公式对多个预设方向中各个预设方向的总能量进行去噪处理：
其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，1≤io≤No，No表示所述预设方向的总数，n＝(n1,n2)，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量，T表示预设的去噪阈值。
在本实施例中，可以根据所述去噪前的总能量确定上述预设的去噪阈值。例如，在一个实施例中，上述去噪阈值T可以为TotalE(io,n)×10-6。
由于现有技术中，采用局部倾斜叠加方法、平面波分解滤波器法、基于相干度量的扫描搜索法以及梯度结构张量法来确定地层倾角时，对随机噪声没有有效的应对策略，当地震数据本身的信噪比较低时，往往难以正确估算出有效反射层的局部倾角信息。在本实施例中，通过上述阈值去噪处理，可以有效地降低随机噪声对待测地震数据处地层倾角计算准确性的影响。
在本实施例中，可以利用曲波变换将待测地震数据变换至能够对待测地震数据进行稀疏表示的曲波域，并在曲波域中进行具有良好的多尺度性、多方向性和各向异性的总能量计算，得到各个待测地震数据在各个预设方向以及各个预设尺度下的波前、同相轴等线状等细节特征，解决了现有技术中由于没有将断层、不整合等地层边缘区域与地震同相轴延展性良好的区域区别对待，所导致的将上述地层边缘区域的倾角估计为断层两侧同相轴的平均倾角的缺陷，能够更加有效地稀疏表示地震数据。
S102：确定多个总能量中的最大值，并确定最大值对应的方向。
S103：计算最大值对应方向的地层倾角，将确定的地层倾角作为待测地震数据相应的地层倾角。
在本实施例中，可以选取S101中所得到的各个方向的噪声压制后的总能量中的最大值所对应的方向，并基于该最大值所对应的方向，计算地层局部倾角。
在本实施例中，可以按照以下公式计算所述待测地震数据相应的地层倾角：
Dip(n)＝IO(n)·360No+180No
其中，Dip(n)表示待测地震数据在下标为n时相应的地层倾角，IO(n)表示预设方向为io、下标为n时的总能量TotalE(io,n)中的最大值对应的方向编号，No表示预设方向的总数。
在本申请的一个实施例中，No可以为用户给定的最大方向数。例如，No可以为8、16或者32的整数值。当所设定的No越大，所确定的曲波变换精度越高，所需要的存储空间(即，内存)也就越大。
在本实施例中，可以应用所得到的地层倾角进行地质勘探。即，可以将待测地震数据相应的地层倾角应用于波场分离与成像、构造导向滤波、反演参数建模、地层自动拾取与拉平等地质方面。
例如：在基于倾角控制的构造导向滤波中，将计算得到的原始地震资料的倾角信息作为约束条件，与各向异性扩散滤波技术相结合，来进行地震资料的保边处理，从而可以有效抑制随机噪声、增强同相轴的一致连续性、保持原始地震图像的结构信息、突出断裂的边界特征。
在本发明实施例中，计算待测地震数据在多个预设方向中各个预设方向的总能量，并基于多个总能量中的最大值对应的方向，计算待测地震数据相应的地层倾角。由于在曲波域对待测地震数据进行处理时能够实现地震数据的稀疏表示，即，能够从多个尺度上对多个方向的地震数据反射特征的细节进行定量化表征，从而解决了现有技术中没有将断层、不整合等地层边缘区域与地震同相轴延展性良好的区域区别对待，所导致的计算精度下降的缺陷。根据上述定量化表征结果，确定待测地震数据中反射同相轴的地层倾角信息，从而在保证计算效率的同时保证较高准确度的局部倾角计算结果，从而提高以数据驱动为特点地震资料处理、解释和反演方法的性能。
下面结合一个具体的实施例对上述地层倾角的确定方法进行具体说明，然而值得注意的是，该具体实施例仅是为了更好地说明本发明，并不构成对本发明的不当限定。
对某地区进行地层倾角的确定，如图2所示为输入的待测叠后地震剖面的示意图，计算图2所示的地震剖面的地层倾角，其中，用户所采用的最大预设尺度Ns为8，最大方向数No为32，预设的去噪阈值T为100。如图3所示为采用本发明实施例所提出的地层倾角的确定方法确定某地区的地层倾角示意图。如图4所示为利用确定的地层倾角进行平面波解构的结果示意图。从图3中颜色较深的起伏图形可知：图2所示的地震剖面相应的地层倾角空间连续性较好，抗噪能力强，特别是在断层及陡倾角构造位置的倾角刻画较为准确。利用图3所示的地层倾角对图2中的叠后地震剖面进行平面波解构，结果如图4，观察图4可知：叠后地震剖面中的有效反射信息基本上都被解构了，从而可以验证本发明所提出的地层倾角确定方法的准确性。
基于同一发明构思，本发明实施例中还提供了一种地层倾角的确定装置，如下面的实施例所述。由于地层倾角的确定装置解决问题的原理与地层倾角的确定方法相似，因此地层倾角的确定装置的实施可以参见地层倾角的确定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图5是本发明实施例的地层倾角的确定装置的一种结构框图，如图5所示，可以包括：总能量计算模块501、方向确定模块502、地层倾角确定模块503，下面对该结构进行说明。
总能量计算模块501，可以用于计算待测地震数据在多个预设方向中各个预设方向的总能量；
方向确定模块502，可以用于确定多个总能量中的最大值，并确定所述最大值对应的方向；
地层倾角确定模块503，可以用于计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角。
利用上述各实施例所提供的地层倾角的确定装置的实施方式，可以自动实施所述地层倾角的确定方法，对地层倾角进行预测，可以不需要实施人员的具体参与，直接输出地层倾角的预测结果，操作简单快捷，有效提高了用户体验。
需要说明的，上述地层倾角的确定装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本申请并不局限于必须是符合本申请实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取/计算/判断等获取的实施例，仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的单元、装置或模块等，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，移动终端，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述，各个实施例之间相同或相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种地层倾角的确定方法，其特征在于，包括：
计算待测地震数据在多个预设方向中各个预设方向的总能量；
确定多个总能量中的最大值，并确定所述最大值对应的方向；
计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角；
其中，计算所述待测地震数据在多个预设方向中各个预设方向的总能量，包括：
根据所述待测地震数据的曲波变换系数，计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量；
计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；
计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和；
根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。
2.如权利要求1所述的方法，其特征在于，确定多个总能量中的最大值，包括：
对所述多个预设方向中各个预设方向的总能量进行去噪处理；
根据去噪处理后的各个预设方向的总能量，确定多个总能量中的最大值。
3.如权利要求2所述的方法，其特征在于，按照以下公式对所述多个预设方向中各个预设方向的总能量进行去噪处理：
其中，ThreshTotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时进行去噪处理后的总能量，1≤io≤No，No表示所述预设方向的总数，n＝(n1,n2)，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量，T表示预设的去噪阈值。
4.如权利要求3所述的方法，其特征在于，所述预设的去噪阈值为TotalE(io,n)×10-6。
5.如权利要求1所述的方法，其特征在于，根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量，包括：
对所述实部和进行归一化处理，得到归一化后的实部和；
对所述虚部和进行归一化处理，得到归一化后的虚部和；
根据所述实部分量、所述虚部分量、所述归一化后的实部和以及所述归一化后的虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。
6.如权利要求1所述的方法，其特征在于，按照以下公式计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量：
其中，RS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，n＝(n1,n2)，0≤n1≤Nt,0≤n2≤Nx表示所述待测地震数据的下标，Nt表示所述待测地震数据在时间方向的采样点数，Nx表示所述待测地震数据在空间方向的采样点数，Coef(is,io,k)表示所述曲波变换系数，is表示所述预设尺度，1≤is≤Ns，io表示所述预设方向，1≤io≤No，k＝(k1,k2)表示所述曲波变换系数的下标，表示ψis,io,k(n1,n2)的共轭复数，ψis,io,k(n1,n2)表示在预设尺度为is、预设方向为io、下标为k时的离散逆曲波，RCoef(is,io,k)表示置Coef(is,io,k)的虚部为零，ICoef(is,io,k)表示置Coef(is,io,k)的实部为零。
7.如权利要求1所述的方法，其特征在于，按照以下公式计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量：
其中，TmpA(is,io,n)＝RS(is,io,n)·MeanSRS(io,n)+IS(is,io,n)·MeanSIS(io,n)，
TmpB(is,io,n)＝RS(is,io,n)·MeanSIS(io,n)-IS(is,io,n)·MeanSRS(io,n)；
其中，RS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的实部分量，IS(is,io,n)表示所述待测地震数据在预设尺度为is、预设方向为io、下标为n时的虚部分量，MeanSRS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的实部和，MeanSIS(io,n)表示进行标准化处理后的预设方向为io、下标为n时的虚部和，TotalE(io,n)表示所述待测地震数据在预设方向为io、下标为n时的总能量。
8.如权利要求1所述的方法，其特征在于，按照以下公式计算所述待测地震数据相应的地层倾角：
Dip(n)＝IO(n)·360/No+180/No
其中，Dip(n)表示所述待测地震数据在下标为n时相应的地层倾角，IO(n)表示预设方向为io、下标为n时的总能量TotalE(io,n)中的最大值对应的方向编号，No表示所述预设方向的总数。
9.如权利要求1所述的方法，其特征在于，在计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角之后，所述方法还包括：
利用所述地层倾角进行地质勘探。
10.一种地层倾角的确定装置，其特征在于，包括：
总能量计算模块，用于计算待测地震数据在多个预设方向中各个预设方向的总能量；
方向确定模块，用于确定多个总能量中的最大值，并确定所述最大值对应的方向；
地层倾角确定模块，用于计算所述最大值对应方向的地层倾角，将确定的地层倾角作为所述待测地震数据相应的地层倾角；
其中，所述总能量计算模块具体用于根据所述待测地震数据的曲波变换系数，计算得到所述待测地震数据在多个预设方向在多个预设尺度下的实部分量以及虚部分量；
计算所述多个预设方向中各个预设方向在多个预设尺度下的实部分量之和，得到多个预设方向中各个预设方向的实部和；
计算所述多个预设方向中各个预设方向在多个预设尺度下的虚部分量之和，得到多个预设方向中各个预设方向的虚部和；
根据所述实部分量、所述虚部分量、所述实部和以及所述虚部和，计算得到所述待测地震数据在多个预设方向中各个预设方向的总能量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710740537.5A CN107703541B (zh) | 2017-08-25 | 2017-08-25 | 地层倾角的确定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710740537.5A CN107703541B (zh) | 2017-08-25 | 2017-08-25 | 地层倾角的确定方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107703541A CN107703541A (zh) | 2018-02-16 |
| CN107703541B true CN107703541B (zh) | 2019-08-06 |
Family
ID=61169799
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710740537.5A Active CN107703541B (zh) | 2017-08-25 | 2017-08-25 | 地层倾角的确定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107703541B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111273351B (zh) * | 2019-11-21 | 2022-04-08 | 西安工业大学 | 用于地震资料去噪的结构导向方向广义全变差正则化方法 |
| CN113419282B (zh) * | 2021-07-19 | 2022-02-15 | 中国矿业大学(北京) | 垂直地震剖面波场分离方法和系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8463551B2 (en) * | 2009-11-17 | 2013-06-11 | Schlumberger Technology Corporation | Consistent dip estimation for seismic imaging |
| CN103257363B (zh) * | 2013-05-21 | 2015-09-30 | 西南石油大学 | 一种探测地下裂缝性储层中裂缝倾角的方法 |
| CN105259571B (zh) * | 2014-07-15 | 2018-03-09 | 中国石油化工股份有限公司 | 一种地层倾角检测方法 |
| CN106019368A (zh) * | 2016-06-17 | 2016-10-12 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于层位模型的倾角计算方法 |
-
2017
- 2017-08-25 CN CN201710740537.5A patent/CN107703541B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107703541A (zh) | 2018-02-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Calvet et al. | Lapse‐time dependence of coda Q: Anisotropic multiple‐scattering models and application to the Pyrenees | |
| Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
| Bignardi | The uncertainty of estimating the thickness of soft sediments with the HVSR method: A computational point of view on weak lateral variations | |
| Ripperger et al. | Variability of near-field ground motion from dynamic earthquake rupture simulations | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| CN105425292B (zh) | 一种油气预测方法及装置 | |
| CN106855640B (zh) | 一种海底电缆地震数据初至计算方法及装置 | |
| CN109100796A (zh) | 一种天然气水合物地震数据处理方法及装置 | |
| CN107607990B (zh) | 三分量检波器水平分量的方向检测方法和装置 | |
| CN108983285A (zh) | 一种基于矩张量的多种震源波场模拟方法及装置 | |
| CN110023790B (zh) | 地震采集几何全波形反演 | |
| CN109061734A (zh) | 一种地震数据的处理方法及装置 | |
| CN104656139A (zh) | 一种转换波的静校正方法及设备 | |
| CN107703541B (zh) | 地层倾角的确定方法和装置 | |
| US6141622A (en) | Seismic semblance/discontinuity method | |
| Zong et al. | Young’s modulus variation with azimuth for fracture-orientation estimation | |
| CN108181649A (zh) | 一种介质分界面振幅补偿方法及装置 | |
| Ktenidou et al. | Directional dependence of site effects observed near a basin edge at Aegion, Greece | |
| CN109100803B (zh) | 微断裂的确定方法和装置 | |
| CN107576985A (zh) | 一种地震反演的方法和装置 | |
| CN106353807B (zh) | 裂缝识别方法和装置 | |
| CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
| CN107870361B (zh) | 一种地震回折波层析成像方法、装置及终端设备 | |
| CN105301638B (zh) | 一种提取风化层底界面的方法和装置 | |
| CN105093300B (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 |