CN109116415B - 地震波数据分离方法、装置及存储介质 - Google Patents
地震波数据分离方法、装置及存储介质 Download PDFInfo
- Publication number
- CN109116415B CN109116415B CN201811184216.2A CN201811184216A CN109116415B CN 109116415 B CN109116415 B CN 109116415B CN 201811184216 A CN201811184216 A CN 201811184216A CN 109116415 B CN109116415 B CN 109116415B
- Authority
- CN
- China
- Prior art keywords
- wave data
- seismic wave
- line
- data
- depth domain
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 56
- 238000012545 processing Methods 0.000 claims abstract description 108
- 230000009466 transformation Effects 0.000 claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 44
- 230000008859 change Effects 0.000 claims abstract description 13
- 238000009795 derivation Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 156
- 230000006870 function Effects 0.000 claims description 24
- 238000004590 computer program Methods 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 description 15
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 101100181929 Caenorhabditis elegans lin-3 gene Proteins 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 230000001131 transforming effect Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 230000009469 supplementation Effects 0.000 description 4
- 101100074846 Caenorhabditis elegans lin-2 gene Proteins 0.000 description 3
- 101100497386 Mus musculus Cask gene Proteins 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 108091053735 lin-4 stem-loop Proteins 0.000 description 1
- 108091032363 lin-4-1 stem-loop Proteins 0.000 description 1
- 108091028008 lin-4-2 stem-loop Proteins 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000000126 substance Substances 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
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
技术领域
本发明涉及地震勘探技术领域，特别涉及地震波数据分离方法、装置及存储介质。
背景技术
地震波是指从震源产生并向四周辐射的弹性波。在结构简单的层状地层界面处地震波会发生波场规则反射，而在构造复杂的块状地层界面处地震波则发生波场杂乱散射，因此，采集到的地震波数据中通常包含反射波数据和散射波数据。实际应用中，可以对地震波数据进行成像处理，得到地震波数据成像，然后根据地震波数据成像对地质构造进行研究。但是，在对包含反射波数据和散射波数据的地震波数据进行成像处理时，根据反射波数据形成的反射波成像和根据散射波数据形成的散射波成像会同时生成在地震波数据成像中，而且，由于散射波数据表征的地震波能量较低，散射波成像还会被反射波成像遮挡，从而影响对地质构造中块状地层的研究准确度。基于此，需要事先对地震波数据中的反射波数据和散射波数据进行分离，以便能够对分离后的反射波数据和散射波数据分别进行处理，提高地质构造研究的准确度。
目前，对原始叠前地震波数据中的叠前反射波数据和叠前散射波数据进行分离是分离地震波数据的主要手段。具体地，可以先对待研究的目标地层进行放炮作业，以收集放炮作用产生的原始叠前地震波数据体，然后，对原始叠前地震波数据体进行时间偏移，得到时间域叠前地震波数据体，之后，分别利用反射波波动方程和散射波波动方程对时间域叠前地震波数据体进行迭代计算，得到时间域叠前反射波数据体和时间域叠前散射波数据体，并将时间域叠前反射波数据体和时间域叠前散射波数据体，作为从地震波数据中分离的反射波数据体和散射波数据体。之后，分别对时间域叠前反射波数据体和时间域叠前散射波数据体进行空间偏移，得到深度域叠前反射波数据体和深度域叠前散射波数据体，再分别对深度域叠前反射波数据体和深度域叠前散射波数据体进行叠加处理，得到深度域叠后反射波数据体和深度域叠后散射波数据体。其中，深度域叠后反射波数据体和深度域叠后散射波数据体是真正能够用于对目标地层的地质构造进行研究的数据。
在实现本发明的过程中，本发明人发现现有技术中至少存在以下问题：
相关技术中是对时间域叠前地震波数据进行分离，分离后的数据还需要分别进行叠加处理，才能得到分离后的数据各自对应的叠后地震波数据。而且，在对时间域叠前地震波数据进行分离时，需要分别利用反射波波动方程和散射波波动方程对时间域叠前地震波数据各进行迭代计算，计算过程复杂，分离效率较低。
发明内容
本发明实施例提供了一种地震波数据分离方法、装置及存储介质，可以用来解决相关技术中对地震波数据体中的反射波数据和散射波数据分离时计算过程复杂、分离效率较低的问题。所述技术方案如下：
第一方面，提供了一种地震波数据分离方法，所述方法包括：
获取目标工区的第一深度域叠后地震波数据体，所述第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，所述第二深度域叠后地震波数据体是根据所述目标工区的原始叠前地震波数据体确定得到；
通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况；
将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体。
可选地，所述获取目标工区的第一深度域叠后地震波数据体，包括：
获取所述第二深度域叠后地震波数据体，所述第二深度域叠后地震波数据体包括所述多个样点对应的深度域叠后地震波数据，所述多个样点位于所述目标工区中；
确定所述多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
将目标网格点对应的深度域叠后地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
将赋值后的所述目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为所述第一深度域叠后地震波数据体。
可选地，所述通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据，包括：
令r＝1，将所述第一深度域叠后地震波数据体确定为第r个地震波数据体；
将所述第r个地震波数据体确定为目标地震波数据体；
通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数；
若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述深度域叠后反射波数据体。
可选地，所述通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，包括：
令s＝2，获取所述目标地震波数据体中在所述预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在所述坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
利用热传导方程的离散函数，对所述第s-1、s、s+1条线数据进行迭代处理，得到所述第s条线数据的迭代结果；
获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据；
利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果；
判断所述s+2是否等于M，所述M为所述目标地震波数据体在所述坐标轴方向上的线数据的条数；
若所述s+2不等于所述M，则令所述s＝s+1，返回获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据的步骤；
若所述s+2等于所述M，则将依次排列的第1条线数据、所述第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及所述第M条线数据，确定为所述第r+1个地震波数据体。
可选地，所述利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果，包括：
通过以下公式，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,0+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述lin(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s)|i,j+1,k,1为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述lin(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k,1为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述ai为所述i对应的坐标轴方向上的预设迭代系数，所述aj为所述j对应的坐标轴方向上的预设迭代系数，所述ak为所述k对应的坐标轴方向上的预设迭代系数。
第二方面，提供了一种地震波数据分离装置，所述装置包括：
获取模块，用于获取目标工区的第一深度域叠后地震波数据体，所述第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，所述第二深度域叠后地震波数据体是根据所述目标工区的原始叠前地震波数据体确定得到；
处理模块，用于通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况；
确定模块，用于将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体。
可选地，所述获取模块包括：
第一获取单元，用于获取所述第二深度域叠后地震波数据体，所述第二深度域叠后地震波数据体包括多个样点对应的深度域叠后地震波数据，所述多个样点位于所述目标工区中；
第一确定单元，用于确定所述多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
第二确定单元，用于根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
赋值单元，用于将目标网格点对应的深度域叠后地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
第三确定单元，用于将赋值后的所述目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为所述第一深度域叠后地震波数据体。
可选地，所述处理模块包括：
第四确定单元，用于令r＝1，将所述第一深度域叠后地震波数据体确定为第r个地震波数据体；
第五确定单元，用于将所述第r个地震波数据体确定为目标地震波数据体；
第一处理单元，用于通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数；
第一判断单元，用于若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
第二判断单元，用于若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述深度域叠后反射波数据体。
可选地，所述第一处理单元包括：
第一获取子单元，用于令s＝2，获取所述目标地震波数据体中在所述预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在所述坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
第一处理子单元，用于利用热传导方程的离散函数，对所述第s-1、s、s+1条线数据进行迭代处理，得到所述第s条线数据的迭代结果；
第二获取子单元，用于获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据；
第二处理子单元，用于利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果；
判断子单元，用于判断所述s+2是否等于M，所述M为所述目标地震波数据体在所述坐标轴方向上的线数据的条数；
若所述s+2不等于所述M，则令所述s＝s+1，返回获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据的步骤；
若所述s+2等于所述M，则将依次排列的第1条线数据、所述第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及所述第M条线数据，确定为所述第r+1个地震波数据体。
可选地，所述第二处理子单元具体用于：
通过以下公式，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,1+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述lin(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s)|i,j+1,k,1为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述lin(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k,1为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述ai为所述i对应的坐标轴方向上的预设迭代系数，所述aj为所述j对应的坐标轴方向上的预设迭代系数，所述ak为所述k对应的坐标轴方向上的预设迭代系数。
第三方面，提供了一种地震波数据分离装置，所述装置包括：
处理器及用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为执行上述第一方面所述的任一项方法。
第四方面，提供了一种计算机可读存储介质，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现上述第一方面提供的任一所述的方法。
本发明实施例提供的技术方案至少带来以下有益效果：在本发明实施例中，可以获取经过对目标工区的原始叠前地震波数据体进行变换后得到的第一深度域叠后地震波数据体，先通过基于热传导方法推导得到的用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况的热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到深度域叠后反射波数据体，再将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体，从而实现了对深度域叠后地震波数据体的直接分离。由于本发明实施例对第一深度域叠后地震波数据进行多次迭代计算后，便可从地震波数据中分离出叠后反射波数据，因此分离过程简单，分离效率高。并且，通过热传导变换方程对第一深度域叠后地震波数据进行数据分离，还可以在减小数据分离运算量的基础上，提高数据分离结果的准确度，使分离后的反射波数据和散射波数据更符合目标工区的真实地质构造。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种地震波数据分离方法的流程示意图；
图2是本发明实施例提供的另一种地震波数据分离方法的流程示意图；
图3是目标工区的第一深度域叠后地震波数据体成像图；
图4是利用本发明实施例分离得到的目标工区的深度域叠后反射波数据体成像图；
图5是利用本发明实施例分离得到的目标工区的深度域叠后散射波数据体成像图；
图6是本发明实施例提供的一种地震波数据分离装置的结构示意图；
图7是本发明实施例提供的一种终端700的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
在对本发明实施例进行详细的解释说明之前，先对本发明实施例中涉及到的名词、应用场景及系统架构分别进行解释说明。
首先，对本发明实施例中涉及到的名词进行介绍。
深度域叠后地震波数据体
深度域叠后地震波数据体是指通过对收集到的地震波在地层中的传播速度的数据进行静校正、动校正、叠加和偏移等处理后得到的，可以反应地震波振幅与地层深度之间的关系的地震波数据体。
热传导方法
热传导方法是指用描述区域内的温度随时间变化的偏微分方程来确定该区域内不同分布位置的温度的方法。
图1是本发明实施例提供的一种地震波数据分离方法的流程示意图。参见图1，该方法包括如下步骤：
步骤101：获取目标工区的第一深度域叠后地震波数据体，第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，第二深度域叠后地震波数据体是根据目标工区的原始叠前地震波数据体确定得到。
步骤102：通过热传导变换方程，对第一深度域叠后地震波数据体进行迭代处理，得到目标工区的深度域叠后反射波数据体，热传导变换方程是基于热传导方法推导得到，用于描述目标工区内的深度域叠后地震波数据随时间的变化情况。
步骤103：将第一深度域叠后地震波数据体和深度域叠后反射波数据体之间的差值，确定为目标工区的深度域叠后散射波数据体。
在本发明实施例中，可以获取经过对目标工区的原始叠前地震波数据体进行变换后得到的第一深度域叠后地震波数据体，先通过基于热传导方法推导得到的用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况的热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到深度域叠后反射波数据体，再将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体，从而实现了对深度域叠后地震波数据体的直接分离。由于本发明实施例对第一深度域叠后地震波数据进行多次迭代计算后，便可从地震波数据中分离出叠后反射波数据，因此分离过程简单，分离效率高。并且，通过热传导变换方程对第一深度域叠后地震波数据进行数据分离，还可以在减小数据分离运算量的基础上，提高数据分离结果的准确度，使分离后的反射波数据和散射波数据更符合目标工区的真实地质构造。
可选地，获取目标工区的第一深度域叠后地震波数据体，包括：
获取第二深度域叠后地震波数据体，第二深度域叠后地震波数据体包括多个样点对应的深度域叠后地震波数据，所述多个样点位于所述目标工区中；
确定多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
根据多个样点在预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，目标空间网格的外形为立方体，且目标空间网格的网格点包括多个样点；
将目标网格点对应的深度域叠后地震波数据赋值为零，目标网格点是指目标空间网格中除多个样点之外的网格点；
将赋值后的目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为第一深度域叠后地震波数据体。
可选地，通过热传导变换方程，对第一深度域叠后地震波数据体进行迭代处理，得到目标工区的深度域叠后反射波数据，包括：
令r＝1，将第一深度域叠后地震波数据体确定为第r个地震波数据体；
将第r个地震波数据体确定为目标地震波数据体；
通过热传导变换方程，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断r+1是否等于N，N为大于或等于3的整数；
若r+1不等于N，则令r＝r+1，并返回将第r个地震波数据体确定为目标地震波数据体的步骤；
若r+1等于N，则将第r+1个地震波数据体确定为深度域叠后反射波数据体。
可选地，通过热传导变换方程，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，包括：
令s＝2，获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
利用热传导方程的离散函数，对第s-1、s、s+1条线数据进行迭代处理，得到第s条线数据的迭代结果；
获取目标地震波数据体中在坐标轴方向上的第s+2条线数据；
利用热传导方程的离散函数，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果；
判断s+2是否等于M，M为目标地震波数据体在坐标轴方向上的线数据的条数；
若s+2不等于M，则令s＝s+1，返回获取目标地震波数据体中在坐标轴方向上的第s+2条线数据的步骤；
若s+2等于M，则将依次排列的第1条线数据、第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及第M条线数据，确定为第r+1个地震波数据体。
可选地，利用热传导方程的离散函数，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果，包括：
通过以下公式，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,0+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，lin(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，lin(s)|i,j+1,k,1为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，lin(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k,1为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，ai为i对应的坐标轴方向上的预设迭代系数，aj为j对应的坐标轴方向上的预设迭代系数，ak为k对应的坐标轴方向上的预设迭代系数。
上述所有可选技术方案，均可按照任意结合形成本发明的可选实施例，本发明实施例对此不再一一赘述。
图2是本发明实施例提供的另一种地震波数据分离方法的流程示意图，该方法可以应用于终端中，该终端可以为手机、平板电脑或计算机等。参见图2，该方法包括如下步骤：
步骤201：获取目标工区的第一深度域叠后地震波数据体。
需要说明的是，目标工区是指目标工作区域，具体可以为油藏上的任一待处理的工作区域。而且，目标工区包括多个样点，这多个样点是在目标工区上预先设置的样点，用于对目标工区的空间位置进行定位，且多个样点中的每个样点均可以存储目标工区在该样点位置处的深度域叠后地震波数据。
再者，第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，比如对第二深度域叠后地震波数据体进行数据补充得到。第二深度域叠后地震波数据体是根据目标工区的原始叠前地震波数据体确定得到，比如对目标工区的原始叠前地震波数据体进行过预处理得到。而且，第二深度域叠后地震波数据体包括该多个样点对应的地震波数据。
实际应用中，目标工区的第二深度域叠后地震波数据体可以由用户输入得到，也可以由其他设备发送得到，也可以由终端对目标工区的原始叠前地震波数据体进行预处理得到。本发明的实施例对第二深度域叠后地震波数据体的获取方法不做具体限定。
可选地，在终端获取到目标工区的原始叠前地震波数据体之后，可以将目标工区的原始叠前地震波数据体转换为深度域叠前地震波数据体，再对深度域叠前地震波数据体进行静校正、动校正、叠加和偏移等处理后，便可得到第二深度域叠后地震波数据体。
其中，可以使用时间域叠前地震波数据体向深度域叠前地震波数据体转换的转换公式，对目标工区的原始叠前地震波数据体进行转换，该转换公式可以为landmark、Geoframe等商业性软件内建立的转换公式。由于原始叠前地震波数据体向第二深度域叠后地震波数据体转换过程为本领域技术人员的常用技术手段，本发明实施例对此不做具体解释，也不对原始叠前地震波数据体向第二深度域叠后地震波数据体转换的转换方式做具体限定。
需要说明的是，在应用本发明实施例提供的地震波数据分离方法时，为了便于在后续步骤中对地震波数据进行处理，可以先对第二深度域叠后地震波数据体对应的目标工区中的多个样点进行补充，使补充后的多个样点在立体空间中排布时，组成的空间为规则的立方体，也即是该立体空间的边界为光滑的平面，并且没有分布在立体空间外的样点。例如，由于原始叠前地震波数据体对应的目标工区中的多个样点在立体空间中排列时，组成的空间可能是三棱柱或三棱锥等空间，因此，根据原始叠前地震波数据体确定的第二深度域叠后地震波数据体对应的目标工区中的多个样点在立体空间中排列时，组成的空间也是三棱柱或三棱锥等空间，这种情况下，由于本发明实施例的后续步骤需要在样点排列为规则的立方体空间中进行数据处理，便可对目标工区进行样点补充，使补充后的目标工区中的多个样点组成的立体空间为立方体，并对补充的多个样点进行地震波数据赋值，最后，将补充的立方体空间内多个样点中存储的地震波数据确定为第一深度域叠后地震波数据。
也即是，通过对第二深度域叠后地震波数据体进行数据补充，可以在第二深度域叠后地震波数据体对应的目标工区中增加更多样点，使得增加样点后的目标工区中的多个样点在立体空间中排列时，组成的空间为立方体。并且，第一深度域叠后地震波数据体中即包括第二深度域叠后地震波数据体中的多个样点，也可包括增加的样点，且所增加的每个样点均具有对应的地震波数据。可选地，所增加的每个样点对应的地震波数据可以为零，当然，也可以为其他数据，本发明实施例对此不做限定。
具体地，对第二深度域叠后地震波数据体进行数据补充，得到第一深度域叠后地震波数据体可以包括：获取第二深度域叠后地震波数据体，第二深度域叠后地震波数据体包括目标工区中的多个样点对应的深度域叠后地震波数据；获取多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；根据多个样点在预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，目标空间网格的外形为立方体，且目标空间网格的网格点包括多个样点；将目标空间网格中的目标网格点对应的深度域叠后地震波数据赋值为零，将赋值后的目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为第一深度域叠后地震波数据体，目标网格点是指目标空间网格中除多个样点之外的网格点。
需要说明的是，为了对第二深度域叠后地震波数据体对应的目标工区中的多个样点在立体空间中进行准确定位，在对第二深度域叠后地震波数据体进行数据补充前，还可以在第二深度域叠后地震波数据体中设立以XYZ为坐标轴、以O为原点的预设空间直角坐标系。具体地，可以在目标工区中，将立体空间中的一个边界线上的端点位置处的样点设置为预设空间直角坐标系的原点O，将立体空间中原点O所在的边界面上的互相垂直的两个方向分别设置为预设空间直角坐标系的X轴和Y轴，将垂直于立体空间中原点O所在的边界面的方向设置为预设空间直角坐标系的Z轴。实际应用中，Z轴可以用于指示深度，X轴和Y轴所组成的平面与Z轴垂直，即X轴、Y轴用于指示与深度垂直的平面上相互垂直的任意两个方向。
在一个实施例中，X轴可以用inline号指示，Y轴用xline号指示，Z轴用K号指示。例如，在第二深度域叠后地震波数据体中，用inline号和xline号记录地震波数据在空间中的排列，inline号和xline号在平面上互相垂直，且inline号和xline号包括有记录深度的K号，也即是，第二深度域叠后地震波数据体对应的目标工区中的多个样点是以inline号、xline号和K号为标记点记录该样点在立体空间中的位置，因此，可以在由inline号、xline号和K号组成的目标工区的立体空间中，将inline号、xline号和K号数值最小的样点设置为空间直角坐标系的原点O，将inline号从小到大的方向设置为X轴方向，将xline从小到大的方向设置为Y轴方向，并将K号从小到大的方向设置为Z轴方向。
另外，在建立空间直角坐标系后，第二深度域叠后地震波数据体对应的目标工区中的每个样点在该空间直角坐标系中均对应有一个坐标，根据第二深度域叠后地震波数据体对应的目标工区中的多个样点在该空间直角坐标系中对应的坐标，可以确定该空间直角坐标系中的三个坐标轴中的每个坐标轴方向上的最大坐标值和最小坐标值。进而可以根据多个样点在每个坐标轴方向上的坐标值范围，确定出一个外形为立方体的目标空间网格。
需要说明的是，对于由多个样点在各个坐标轴方向上的坐标值范围组成的目标空间网格而言，目标空间网格中包括的网格点的数量可能会大于或等于该多个样点的数量。当目标空间网格中包括的网格点的数量等于该多个样点的数量时，可以直接将第二深度域叠后地震波数据体确定为第一深度域叠后地震波数据体；当目标空间网格中包括的网格点的数量大于该多个样点的数量时，再将目标空间网格中未与第二深度域叠后地震波数据体中的样点对应的网格点确定为目标网格点，并将目标网格点对应的深度域叠后地震波数据赋值为零，最后，将赋值后的目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为第一深度域叠后地震波数据体。
例如，在一个目标工区中，将空间直角坐标系的原点O建立在inline号为0、xline号为0、K号为0的样点上，X轴为平行于inline号的方向，Y轴为平行于xline号的方向，Z轴为平行于K号的方向，在空间直角坐标系中，inline号在X轴方向的最大坐标值为1480，最小坐标值为0，xline号在Y轴方向的最大坐标值为590，最小坐标值为0，K号在Z轴方向的最大坐标值为1601，最小坐标值为0。此时，根据多个样点在每个坐标轴方向上的坐标值范围，确定出的目标空间网格中包括的网格点的数量等于第二深度域叠后地震波数据体中的样点的数量，因此，可以直接将该第二深度域叠后地震波数据体确定为第一深度域叠后地震波数据体。
又如，在另一个目标工区中，同样将空间直角坐标系的原点O建立在inline号为0、xline号为0、K号为0的样点上，X轴为平行于inline号的方向，Y轴为平行于xline号的方向，Z轴为平行于K号的方向，在该空间直角坐标系中，inline号在X轴方向的最大坐标值为1480，最小坐标值为0，xline号在Y轴方向的最大坐标值为590，最小坐标值为0，并且，当inline号在X轴方向的坐标值为0-1180时，K号在Z轴方向的最大坐标值为1601，最小坐标值为0当inline号在X轴方向的坐标值为1181-1480时，K号在Z轴方向的最大坐标值为1550，最小坐标值为0。此时，根据多个样点在每个坐标轴方向上的坐标值范围，确定出的目标空间网格中包括的网格点的数量大于第二深度域叠后地震波数据体中的样点的数量，因此，将inline号在X轴方向的坐标值为2000-2300，Y轴方向的坐标值为960-1550，Z轴方向的坐标值为1550-1601之间的样点对应的网格点确定为目标网格点，并将该目标网格点对应的深度域叠后地震波数据赋值为零，进而将赋值后的目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为第一深度域叠后地震波数据体。
步骤202：通过热传导变换方程，对第一深度域叠后地震波数据体进行迭代处理，得到目标工区的深度域叠后反射波数据体。
需要说明的是，热传导变换方程是对热传导方法中应用到描述区域内的温度随时间变化的偏微分方程进行变换，使该热传导变换方程可以用来描述目标网格中地震波数据随时间变化的偏微分方程。
具体的，热传导方程为由下述公式(1)示出：
其中，
是空间中一点的温度对时间的变化率，q是热扩散率，
是温度对三个空间坐标轴的二次导数。
通过将上述热传导方程中的温度变换为深度域叠后地震波数据体中的地震波数据，也即是地震波振幅数据，从而可以得到热传导变换方程，热传导变换方程由下述公式(2)示出。
gi,j,k＝L|i,j,k+ax×Lxx|i,j,k+ay×Lyy|i,j,k+az×Lzz|i,j,k) (2)
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，gi,j,k为坐标(i，j，k)上的样点对应的经过热传导变换方程计算后的深度域叠后地震波数据，L|i,j,k为坐标(i，j，k)上的样点对应的原始深度域叠后地震波数据，Lxx|i,j,k为坐标(i，j，k)上的样点沿空间直接坐标系中的X轴方向的二阶导数，Lyy|i,j,k为坐标(i，j，k)上的样点沿空间直接坐标系中的Y轴方向的二阶导数，Lzz|i,j,k为坐标(i，j，k)上的样点沿空间直接坐标系中的Z轴方向的二阶导数，ax为沿空间直接坐标系中的X轴方向的迭代系数，ay为沿空间直接坐标系中的Y轴方向的迭代系数，az为沿空间直接坐标系中的Z轴方向的迭代系数。
具体地，在确定了热传导变换方程后，可以通过以下步骤对第一深度域叠后地震波数据体进行迭代处理，得到目标工区的深度域叠后反射波数据：
步骤2021：令r＝1，将第一深度域叠后地震波数据体确定为第r个地震波数据体。
也即是，先将r赋值为1，并将第一深度域叠后地震波数据体作为第1个地震波数据体进行处理。
步骤2022：将第r个地震波数据体确定为目标地震波数据体，通过热传导变换方程，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体。
需要说明的是，在根据第一深度域叠后地震波数据体确定目标工区的深度域叠后反射波数据时，需要对第一深度域叠后地震波数据体进行N次迭代，以使迭代处理后的得到的目标工区的深度域叠后反射波数据更加准确，也更符合目标工区的真实地质构造。在对第一深度域叠后地震波数据体进行迭代处理时，可以通过预设迭代次数N的数值，并判断迭代次数是否达到了预设迭代次数，如果迭代次数未达到预设迭代次数，则继续执行迭代步骤，如果迭代次数达到预设迭代次数，则停止迭代，并将迭代结果确定为目标工区的深度域叠后反射波数据。
具体地，在进行迭代处理时，可以通过下述迭代步骤1)-3)，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体。
1)令s＝2，获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据，利用热传导方程的离散函数，对第s-1、s、s+1条线数据进行迭代处理，得到第s条线数据的迭代结果。
需要说明的是，在通过热传导变换方程对目标地震波数据体进行迭代处理时，目标地震波数据体的边界位置处的线数据是无法进行迭代的，只能对除目标地震波数据体的边界位置处的线数据之外的线数据进行迭代处理。因此，在第一次获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的三条线数据时，是从边界位置处开始按顺序进行获取三条线数据，而获取之后根据热传导变换方程可以对三条线数据中的第二条线数据进行一次迭代。例如，可以从X轴方向上坐标最小值对应的inline线数据开始，从小到大依次获取三条inline线数据，记为lin1、lin2和lin3，利用公式(2)示出的热传导变换方程对lin1、lin2和lin3进行迭代处理，便可得到lin2的迭代结果，并将lin2的迭代结果记为lin21。
2)获取目标地震波数据体中在坐标轴方向上的第s+2条线数据，利用热传导方程的离散函数，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果。
在对第s-1、s、s+1条线数据完成一次迭代处理后，可以得到第s条线数据的迭代结果，并将第s-1条线数据和第s条线数据的迭代结果存入用于保存线数据的临时中间文件。之后，可以对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果，并继续将第s+1条线数据的迭代结果存入临时中间文件。
可选地，可以通过以下公式(3)，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,1+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0} (3)
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，lin(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，lin(s)|i,j+1,k,1为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，lin(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k,1为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，ai为i对应的坐标轴方向上的预设迭代系数，aj为j对应的坐标轴方向上的预设迭代系数，ak为k对应的坐标轴方向上的预设迭代系数。
需要说明的是，该公式为一个3×3×3的空间网格处理模型，在利用该公式对多条线数据进行迭代处理时，可以依次读取3条线数据，在对该3条线数据中的第2条线数据上除边界位置上的样点之外的每个样点中存储的地震波数据完成一次迭代处理后，便可以得到第2条线数据的迭代结果，之后，将3条线数据的位置进行变换，将第1条线数据存入临时中间文件，用第2条线数据的迭代结果替换第1条线数据，用第3条线数据替换第2条线数据的迭代结果，重新读取第4条线数据，并用第4条线数据替换第3条线数据，从而可以继续对第2条线数据的迭代结果、第3条线数据和第4条线数据利用公式(3)进行迭代处理。
例如，在对lin1、lin2和lin3通过公式(2)进行一次迭代处理后得到lin21，将lin1存入临时中间文件，用lin21替换lin1，用lin3替换lin2，读取lin4，并用lin4替换lin3，这样，对于线数据lin21、lin3和lin4，取迭代系数为0.1，则根据公式(3)可以对lin3在(2，2，2)坐标处的样点进行迭代处理，得到如下结果：
lin(3)|2,2,2,1＝lin(3)|2,2,2,0+0.1×{lin(3)|3,2,2,0+lin(3)|1,2,2,0-2×lin|2,2,2,0}+0.1×{lin(2)|2,3,2,1+lin(4)|2,1,2,0-2×lin(3)|2,2,2,0}+0.1×{lin(3)|2,2,3,0+lin(3)|2,2,1,0-2×lin(3)|2,2,2,0} (4)
在对lin3在(2，2，2)坐标处的样点进行迭代处理后，还可以通过公式(3)对lin3在除(2，2，2)坐标之外的其他位置的样点所存储的地震波数据进行迭代处理，在完成对lin3上的除边界位置上的样点之外的其他所有样点的数据处理后，可以得到lin31。
需要说明的是，上述样点坐标、线数据条数、迭代系数均为示例性数据，实际应用中还可以为其他数据，本发明实施例对此不做具体限定。
(3)判断s+2是否等于M，M为目标地震波数据体在坐标轴方向上的线数据的条数，若s+2不等于M，则令s＝s+1，返回获取目标地震波数据体中在坐标轴方向上的第s+2条线数据的步骤，若s+2等于M，则将依次排列的第1条线数据、第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及第M条线数据，确定为第r+1个地震波数据体。
需要说明的是，因为在第r个地震波数据体向第r+1个地震波数据体迭代过程中，两个边界位置处的线数据不进行迭代处理，只对除两个边界位置处的线数据之外的其他线数据进行一次迭代处理，因此，当第r个地震波数据体中有M条线数据时，最多可以迭代M-2次，便可得到第r+1个地震波数据体。因此，可以为迭代过程设置一个判断条件，当s+2等于M时，则说明已经将第r个地震波数据体中的所有线数据迭代完成一次；而当当s+2小于M时，则说明第r个地震波数据体中还有线数据未经过迭代，还需要根据步骤2)继续对第r个地震波数据体进行迭代处理，直至当s+2等于M时，完成迭代全部迭代处理，此时，可以将保存在临时中间文件中的数据作为数据体确定为第r+1个地震波数据体。
步骤2023：判断r+1是否等于N，N为大于或等于3的整数。
步骤2024：若r+1不等于N，则令r＝r+1，并返回步骤2022。
也即是，令r＝r+1，并返回将第r个地震波数据体确定为目标地震波数据体的步骤，如此，可以继续对第r+1个地震波数据体进行迭代处理，得到第r+2个地震波数据体。
步骤2025：若r+1等于N，则将第r+1个地震波数据体确定为深度域叠后反射波数据体。
需要说明的是，在根据第r个地震波数据体确定了第r+1个地震波数据体之后，为了进一步提高数据分离的准确性，使分离之后得到的深度域叠后反射波数据体更符合目标工区的真实地质构造，还可以对第r+1个地震波数据体按照步骤2021进行N次迭代处理。其中，N是预设的迭代次数，可以根据N的取值控制迭代后得到的不同精度的深度域叠后反射波数据体。在实际应用中，N可以由用户输入得到，可以由其他设备发送得到，本发明实施例对此不做具体限定。例如，可以由用户在终端中输入N的数值，如10、20或30等，N的数值越大，则迭代处理过程所需要的时间越长，得到的深度域叠后反射波数据体的精度也越高，本发明实施例对N的数值不做具体限定。
另外，为了使深度域叠后反射波数据体可以准确反映目标工区的真实地质构造，可以将N设置为大于或等于3的整数，也即是，对第一深度域叠后地震波数据体按步骤2021所示出的方法，迭代至少两次，得到深度域叠后反射波数据体。
具体地，在判断r+1与N的关系时，当r+1不等于N时，则使r＝r+1，并重新根据将第r个地震波数据体确定目标地震波数据体的步骤确定新的目标地震波数据体，当r+1等于N时，则可以将第r+1个地震波数据体确定为深度域叠后反射波数据体。
步骤203：将第一深度域叠后地震波数据体和深度域叠后反射波数据体之间的差值，确定为目标工区的深度域叠后散射波数据体。
需要说明的是，由于深度域叠后反射波数据体中的每个样点所对应的地震波数据是对第一深度域叠后地震波数据体中的每个样点所对应的地震波数据进行迭代处理后得到的，因此，深度域叠后反射波数据体中的每个样点所对应的地震波数据和深度域叠后反射波数据体中的每个样点所对应的地震波数据的数据类型是一致的，两者中的样点位置也是一致的。并且，在假设第一深度域叠后地震波数据体中不含有噪音或噪音能量很弱的基础上，可以认为第一深度域叠后地震波数据体中只含有深度域叠后反射波数据体和深度域叠后散射波数据体。基于此，在确定了深度域叠后反射波数据体之后，可以直接将第一深度域叠后地震波数据体和深度域叠后反射波数据体进行相减，并将第一深度域叠后地震波数据体和深度域叠后反射波数据体之间的差值，确定为目标工区的深度域叠后散射波数据体。
例如，在目标工区中，第一深度域叠后地震波数据体成像图如图3所示，利用本发明实施例提供的地震波数据分离方法对第一深度域叠后地震波数据体进行地震波数据分离，得到深度域叠后反射波数据体和深度域叠后散射波数据体，将深度域叠后反射波数据体记为DATA2，深度域叠后反射波数据体主要表现为反射波能量数据体，代表第一深度域叠后地震波数据体DATA1的骨架信息，具有更好的信噪比和连续性，其深度域成像剖面如图4所示，将第一深度域叠后地震波数据体DATA1和深度域叠后反射波数据体深度域叠后反射波数据体DATA2进行相减，将减掉部分数据体记为DATA3，则DATA3即为目标工区的深度域叠后散射波数据体，主要表现为散射波能量数据体，代表第一深度域叠后地震波数据体DATA1的细节信息，其深度域成像剖面如图5所示，最终实现深度域叠后反射波数据体和深度域叠后散射波数据体从第一深度域叠后地震波数据体中分离。
在本发明实施例中，可以获取经过对目标工区的原始叠前地震波数据体进行变换后得到的第一深度域叠后地震波数据体，先通过基于热传导方法推导得到的用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况的热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到深度域叠后反射波数据体，再将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体，从而实现了对深度域叠后地震波数据体的直接分离。由于本发明实施例对第一深度域叠后地震波数据进行多次迭代计算后，便可从地震波数据中分离出叠后反射波数据，因此分离过程简单，分离效率高。并且，通过热传导变换方程对第一深度域叠后地震波数据进行数据分离，还可以在减小数据分离运算量的基础上，提高数据分离结果的准确度，使分离后的反射波数据和散射波数据更符合目标工区的真实地质构造。
图6是本发明实施例提供的一种地震波数据分离装置的流程示意图。参见图6，该装置可以包括：
获取模块601，用于获取目标工区的第一深度域叠后地震波数据体，第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，第二深度域叠后地震波数据体是根据目标工区的原始叠前地震波数据体确定得到；
处理模块602，用于通过热传导变换方程，对第一深度域叠后地震波数据体进行迭代处理，得到目标工区的深度域叠后反射波数据体，热传导变换方程是基于热传导方法推导得到，用于描述目标工区内的深度域叠后地震波数据随时间的变化情况；
确定模块603，用于将第一深度域叠后地震波数据体和深度域叠后反射波数据体之间的差值，确定为目标工区的深度域叠后散射波数据体。
可选地，获取模块包括：
第一获取单元，用于获取第二深度域叠后地震波数据体，第二深度域叠后地震波数据体包括多个样点对应的深度域叠后地震波数据，多个样点位于目标工区中；
第二获取单元，用于确定多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
第一确定单元，用于根据多个样点在预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，目标空间网格的外形为立方体，且目标空间网格的网格点包括多个样点；
赋值单元，用于将目标网格点对应的深度域叠后地震波数据赋值为零，目标网格点是指目标空间网格中除多个样点之外的网格点；
第二确定单元，用于将赋值后的目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为第一深度域叠后地震波数据体。
可选地，处理模块包括：
第三确定单元，用于令r＝1，将第一深度域叠后地震波数据体确定为第r个地震波数据体；
第四确定单元，用于将第r个地震波数据体确定为目标地震波数据体；
第一处理单元，用于通过热传导变换方程，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断r+1是否等于N，N为大于或等于3的整数；
第一判断单元，用于若r+1不等于N，则令r＝r+1，并返回将第r个地震波数据体确定为目标地震波数据体的步骤；
第二判断单元，用于若r+1等于N，则将第r+1个地震波数据体确定为深度域叠后反射波数据体。
可选地，第一处理单元包括：
第一获取子单元，用于令s＝2，获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
第一处理子单元，用于利用热传导方程的离散函数，对第s-1、s、s+1条线数据进行迭代处理，得到第s条线数据的迭代结果；
第二获取子单元，用于获取目标地震波数据体中在坐标轴方向上的第s+2条线数据；
第二处理子单元，用于利用热传导方程的离散函数，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果；
判断子单元，用于判断s+2是否等于M，M为目标地震波数据体在坐标轴方向上的线数据的条数；
若s+2不等于M，则令s＝s+1，返回获取目标地震波数据体中在坐标轴方向上的第s+2条线数据的步骤；
若s+2等于M，则将依次排列的第1条线数据、第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及第M条线数据，确定为第r+1个地震波数据体。
可选地，第二处理子单元具体用于：
通过以下公式，对第s条线数据的迭代结果、第s+1条线数据、第s+2条线数据进行迭代处理，得到第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,1+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，lin(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，lin(s)|i,j+1,k,1为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，lin(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，lin(s+1)|i,j,k,1为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，ai为i对应的坐标轴方向上的预设迭代系数，aj为j对应的坐标轴方向上的预设迭代系数，ak为k对应的坐标轴方向上的预设迭代系数。
在本发明实施例中，可以获取经过对目标工区的原始叠前地震波数据体进行变换后得到的第一深度域叠后地震波数据体，先通过基于热传导方法推导得到的用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况的热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到深度域叠后反射波数据体，再将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体，从而实现了对深度域叠后地震波数据体的直接分离。由于本发明实施例对第一深度域叠后地震波数据进行多次迭代计算后，便可从地震波数据中分离出叠后反射波数据，因此分离过程简单，分离效率高。并且，通过热传导变换方程对第一深度域叠后地震波数据进行数据分离，还可以在减小数据分离运算量的基础上，提高数据分离结果的准确度，使分离后的反射波数据和散射波数据更符合目标工区的真实地质构造。
需要说明的是：上述实施例提供的地震波数据分离装置在分离地震波数据时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的地震波数据装置与确定地震波数据方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
图7是本发明实施例提供的一种终端700的结构示意图。该终端700可以是：智能手机、平板电脑、MP3播放器(Moving Picture Experts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts Group Audio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端700还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
通常，终端700包括有：处理器701和存储器702。
处理器701可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器701可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器701也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器701可以在集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器701还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器702可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器702还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器702中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器701所执行以实现本申请中方法实施例提供的地震波数据分离方法。
在一些实施例中，终端700还可选包括有：外围设备接口703和至少一个外围设备。处理器701、存储器702和外围设备接口703之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口703相连。具体地，外围设备包括：射频电路704、触摸显示屏704、摄像头706、音频电路707、定位组件708和电源709中的至少一种。
外围设备接口703可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器701和存储器702。在一些实施例中，处理器701、存储器702和外围设备接口703被集成在同一芯片或电路板上；在一些其他实施例中，处理器701、存储器702和外围设备接口703中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路704用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路704通过电磁信号与通信网络以及其他通信设备进行通信。射频电路704将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路704包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路704可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及4G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路704还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏704用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏705是触摸显示屏时，显示屏705还具有采集在显示屏705的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器701进行处理。此时，显示屏705还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏705可以为一个，设置终端700的前面板；在另一些实施例中，显示屏705可以为至少两个，分别设置在终端700的不同表面或呈折叠设计；在再一些实施例中，显示屏705可以是柔性显示屏，设置在终端700的弯曲表面上或折叠面上。甚至，显示屏705还可以设置成非矩形的不规则图形，也即异形屏。显示屏705可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
摄像头组件706用于采集图像或视频。可选地，摄像头组件706包括前置摄像头和后置摄像头。通常，前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件706还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
音频电路707可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器701进行处理，或者输入至射频电路704以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端700的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器701或射频电路704的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路707还可以包括耳机插孔。
定位组件708用于定位终端700的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件708可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。
电源709用于为终端700中的各个组件进行供电。电源709可以是交流电、直流电、一次性电池或可充电电池。当电源709包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
在一些实施例中，终端700还包括有一个或多个传感器710。该一个或多个传感器710包括但不限于：加速度传感器711、陀螺仪传感器712、压力传感器713、指纹传感器714、光学传感器715以及接近传感器716。
加速度传感器711可以检测以终端700建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器711可以用于检测重力加速度在三个坐标轴上的分量。处理器701可以根据加速度传感器711采集的重力加速度信号，控制触摸显示屏705以横向视图或纵向视图进行用户界面的显示。加速度传感器711还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器712可以检测终端700的机体方向及转动角度，陀螺仪传感器712可以与加速度传感器711协同采集用户对终端700的3D动作。处理器701根据陀螺仪传感器712采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器713可以设置在终端700的侧边框和/或触摸显示屏705的下层。当压力传感器713设置在终端700的侧边框时，可以检测用户对终端700的握持信号，由处理器701根据压力传感器713采集的握持信号进行左右手识别或快捷操作。当压力传感器713设置在触摸显示屏705的下层时，由处理器701根据用户对触摸显示屏705的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器714用于采集用户的指纹，由处理器701根据指纹传感器714采集到的指纹识别用户的身份，或者，由指纹传感器714根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器701授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器714可以被设置终端700的正面、背面或侧面。当终端700上设置有物理按键或厂商Logo时，指纹传感器714可以与物理按键或厂商Logo集成在一起。
光学传感器715用于采集环境光强度。在一个实施例中，处理器701可以根据光学传感器715采集的环境光强度，控制触摸显示屏705的显示亮度。具体地，当环境光强度较高时，调高触摸显示屏705的显示亮度；当环境光强度较低时，调低触摸显示屏705的显示亮度。在另一个实施例中，处理器701还可以根据光学传感器715采集的环境光强度，动态调整摄像头组件706的拍摄参数。
接近传感器716，也称距离传感器，通常设置在终端700的前面板。接近传感器716用于采集用户与终端700的正面之间的距离。在一个实施例中，当接近传感器716检测到用户与终端700的正面之间的距离逐渐变小时，由处理器701控制触摸显示屏705从亮屏状态切换为息屏状态；当接近传感器716检测到用户与终端700的正面之间的距离逐渐变大时，由处理器701控制触摸显示屏705从息屏状态切换为亮屏状态。
也即是，本发明实施例不仅提供了一种终端，包括处理器和用于存储处理器可执行指令的存储器，其中，处理器被配置为执行图1或图2所示的实施例中的方法，而且，本发明实施例还提供了一种计算机可读存储介质，该存储介质内存储有计算机程序，该计算机程序被处理器执行时可以实现图1或图2所示的实施例中的地震波数据分离方法。
本领域技术人员可以理解，图7中示出的结构并不构成对终端700的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种地震波数据分离方法，其特征在于，所述方法包括：
获取目标工区的第一深度域叠后地震波数据体，所述第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，所述第二深度域叠后地震波数据体是根据所述目标工区的原始叠前地震波数据体确定得到；
通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况；
将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体；
其中，所述通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据，包括：
令r＝1，将所述第一深度域叠后地震波数据体确定为第r个地震波数据体；
将所述第r个地震波数据体确定为目标地震波数据体；
通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数，所述N为预设迭代次数；
若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述深度域叠后反射波数据体；
其中，所述通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，包括：
令s＝2，获取所述目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在所述坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
利用热传导方程的离散函数，对所述第s-1、s、s+1条线数据进行迭代处理，得到所述第s条线数据的迭代结果；
获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据；
利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果；
判断所述s+2是否等于M，所述M为所述目标地震波数据体在所述坐标轴方向上的线数据的条数；
若所述s+2不等于所述M，则令所述s＝s+1，返回获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据的步骤；
若所述s+2等于所述M，则将依次排列的第1条线数据、所述第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及所述第M条线数据，确定为所述第r+1个地震波数据体；
其中，所述利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果，包括：
通过以下公式，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,1+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述lin(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s)|i,j+1,k,1为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述lin(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k,1为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述ai为所述i对应的坐标轴方向上的预设迭代系数，所述aj为所述j对应的坐标轴方向上的预设迭代系数，所述ak为所述k对应的坐标轴方向上的预设迭代系数。
2.根据权利要求1所述的方法，其特征在于，所述获取目标工区的第一深度域叠后地震波数据体，包括：
获取所述第二深度域叠后地震波数据体，所述第二深度域叠后地震波数据体包括多个样点对应的深度域叠后地震波数据，所述多个样点位于所述目标工区中；
确定所述多个样点在所述预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
将目标网格点对应的深度域叠后地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
将赋值后的所述目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为所述第一深度域叠后地震波数据体。
3.一种地震波数据分离装置，其特征在于，所述装置包括：
获取模块，用于获取目标工区的第一深度域叠后地震波数据体，所述第一深度域叠后地震波数据体是根据第二深度域叠后地震波数据体确定得到，所述第二深度域叠后地震波数据体是根据所述目标工区的原始叠前地震波数据体确定得到；
处理模块，用于通过热传导变换方程，对所述第一深度域叠后地震波数据体进行迭代处理，得到所述目标工区的深度域叠后反射波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域叠后地震波数据随时间的变化情况；
确定模块，用于将所述第一深度域叠后地震波数据体和所述深度域叠后反射波数据体之间的差值，确定为所述目标工区的深度域叠后散射波数据体；
其中，所述处理模块包括：
第一确定单元，用于令r＝1，将所述第一深度域叠后地震波数据体确定为第r个地震波数据体；
第二确定单元，用于将所述第r个地震波数据体确定为目标地震波数据体；
第一处理单元，用于通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数，所述N为预设迭代次数；
第一判断单元，用于若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
第二判断单元，用于若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述深度域叠后反射波数据体；
其中，所述通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，包括：
令s＝2，获取所述目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在所述坐标轴方向上的坐标值相同的至少一个样点所对应的深度域叠后地震波数据；
利用热传导方程的离散函数，对所述第s-1、s、s+1条线数据进行迭代处理，得到所述第s条线数据的迭代结果；
获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据；
利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果；
判断所述s+2是否等于M，所述M为所述目标地震波数据体在所述坐标轴方向上的线数据的条数；
若所述s+2不等于所述M，则令所述s＝s+1，返回获取所述目标地震波数据体中在所述坐标轴方向上的第s+2条线数据的步骤；
若所述s+2等于所述M，则将依次排列的第1条线数据、所述第1条线数据和第M条线数据之间的每条线数据的迭代结果，以及所述第M条线数据，确定为所述第r+1个地震波数据体；
其中，所述利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果，包括：
通过以下公式，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果：
lin(s+1)|i,j,k,1＝lin(s+1)|i,j,k,0+ai×{lin(s+1)|i+1,j,k,0+lin(s+1)|i-1,j,k,0-2×lin(s+1)|i,j,k,0}+aj×{lin(s)|i,j+1,k,1+lin(s+2)|i,j-1,k,0-2×lin(s+1)|i,j,k,0}+ak×{lin(s+1)|i,j,k+1,0+lin(s+1)|i,j,k-1,0-2×lin(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述lin(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域叠后地震波数据，所述lin(s)|i,j+1,k,1为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述lin(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域叠后地震波数据，所述lin(s+1)|i,j,k,1为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域叠后地震波数据的迭代结果，所述ai为所述i对应的坐标轴方向上的预设迭代系数，所述aj为所述j对应的坐标轴方向上的预设迭代系数，所述ak为所述k对应的坐标轴方向上的预设迭代系数。
4.根据权利要求3所述的装置，其特征在于，所述获取模块包括：
第一获取单元，用于获取所述第二深度域叠后地震波数据体，所述第二深度域叠后地震波数据体包括多个样点对应的深度域叠后地震波数据，所述多个样点位于所述目标工区中；
第三确定单元，用于确定所述多个样点在所述预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
第四确定单元，用于根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
赋值单元，用于将目标网格点对应的深度域叠后地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
第五确定单元，用于将赋值后的所述目标空间网格中每个网格点对应的深度域叠后地震波数据，确定为所述第一深度域叠后地震波数据体。
5.一种地震波数据分离装置，其特征在于，所述装置包括：
处理器及用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为执行如权利要求1-2任一项所述的方法。
6.一种计算机可读存储介质，其特征在于，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-2任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811184216.2A CN109116415B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据分离方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811184216.2A CN109116415B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据分离方法、装置及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109116415A CN109116415A (zh) | 2019-01-01 |
| CN109116415B true CN109116415B (zh) | 2020-06-09 |
Family
ID=64857923
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811184216.2A Active CN109116415B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据分离方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109116415B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104090300A (zh) * | 2014-07-18 | 2014-10-08 | 成都晶石石油科技有限公司 | 置信叠前扩散滤波方法 |
| CN104730571A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种利用绕射再聚焦识别小尺度地质体的方法与装置 |
| CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103076626B (zh) * | 2011-10-25 | 2016-02-03 | 中国石油化工股份有限公司 | 一种波场净化处理方法 |
| CN103489163B (zh) * | 2013-09-13 | 2016-01-20 | 电子科技大学 | 基于正则化混合范数滤波的地震图像结构导向降噪方法 |
| GB2547942A (en) * | 2016-03-04 | 2017-09-06 | Robertsson Johan | Method for deghosting and redatuming operator estimation |
| CN106597549B (zh) * | 2016-12-05 | 2017-11-14 | 山东天元信息技术股份有限公司 | 一种反射波提取方法及装置、反射声波成像测井仪器 |
| US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
-
2018
- 2018-10-11 CN CN201811184216.2A patent/CN109116415B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104090300A (zh) * | 2014-07-18 | 2014-10-08 | 成都晶石石油科技有限公司 | 置信叠前扩散滤波方法 |
| CN104730571A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种利用绕射再聚焦识别小尺度地质体的方法与装置 |
| CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
Non-Patent Citations (1)
| Title |
|---|
| 地震反射波与散射波波场分离方法初探;陈可洋 等;《岩性油气藏》;20130430;第25卷(第2期);第76-81页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109116415A (zh) | 2019-01-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200334854A1 (en) | Position and attitude determining method and apparatus, smart device, and storage medium | |
| CN108682038B (zh) | 位姿确定方法、装置及存储介质 | |
| CN111768454B (zh) | 位姿确定方法、装置、设备及存储介质 | |
| CN110986930B (zh) | 设备定位方法、装置、电子设备及存储介质 | |
| CN109166150B (zh) | 获取位姿的方法、装置存储介质 | |
| CN109886208B (zh) | 物体检测的方法、装置、计算机设备及存储介质 | |
| CN109696706B (zh) | 盐构造检测方法、装置及存储介质 | |
| CN110174699B (zh) | 确定断层封闭性的方法、装置及存储介质 | |
| CN111897429A (zh) | 图像显示方法、装置、计算机设备及存储介质 | |
| CN110389380B (zh) | 地震剖面的同相轴自动拾取方法、装置及存储介质 | |
| CN111005715A (zh) | 确定气井产量的方法、装置及存储介质 | |
| CN112922589A (zh) | 尖灭线确定方法、装置、终端及存储介质 | |
| CN109116424B (zh) | 地震波数据的低波数噪音分离方法、装置及存储介质 | |
| CN111127541A (zh) | 车辆尺寸的确定方法、装置及存储介质 | |
| CN111928861B (zh) | 地图构建方法及装置 | |
| CN113160031A (zh) | 图像处理方法、装置、电子设备及存储介质 | |
| CN110673214B (zh) | 预测水平井的入口靶点和末端点的深度的方法和装置 | |
| CN109116415B (zh) | 地震波数据分离方法、装置及存储介质 | |
| CN109388732B (zh) | 音乐地图的生成和显示方法、装置及存储介质 | |
| CN110443841B (zh) | 地面深度的测量方法、装置及系统 | |
| CN113009573A (zh) | 地震数据采集方法、装置、计算机设备及存储介质 | |
| CN110058306B (zh) | 确定三维速度体的方法、装置及计算机存储介质 | |
| CN111650637B (zh) | 地震层位解释方法及装置 | |
| CN112182986B (zh) | 确定面积波及信息的方法、装置、设备及存储介质 | |
| CN113689484B (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 |