CN109116424A - 地震波数据的低波数噪音分离方法、装置及存储介质 - Google Patents
地震波数据的低波数噪音分离方法、装置及存储介质 Download PDFInfo
- Publication number
- CN109116424A CN109116424A CN201811184185.0A CN201811184185A CN109116424A CN 109116424 A CN109116424 A CN 109116424A CN 201811184185 A CN201811184185 A CN 201811184185A CN 109116424 A CN109116424 A CN 109116424A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- depth domain
- data body
- domain seismic
- line number
- 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
- 238000000926 separation method Methods 0.000 title claims abstract description 39
- 238000012545 processing Methods 0.000 claims abstract description 105
- 238000012546 transfer Methods 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000009466 transformation Effects 0.000 claims abstract description 47
- 238000003384 imaging method Methods 0.000 claims abstract description 44
- 239000013589 supplement Substances 0.000 claims abstract description 22
- 238000005070 sampling Methods 0.000 claims description 144
- 230000008859 change Effects 0.000 claims description 67
- 230000006870 function Effects 0.000 claims description 21
- 230000000295 complement effect Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 239000007787 solid Substances 0.000 description 13
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000001914 filtration Methods 0.000 description 5
- 238000012804 iterative process Methods 0.000 description 5
- 101150044251 OGT gene Proteins 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 230000004927 fusion Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 108010001267 Protein Subunits Proteins 0.000 description 2
- 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
- 238000013461 design Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101100181929 Caenorhabditis elegans lin-3 gene Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 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
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
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条线数据的迭代结果：
dat(s+1)|i,j,k,1＝dat(s+1)|i,j,k,0+ai×{dat(s+1)|i+1,j,k,0+dat(s+1)|i-1,j,k,0-2×dat(s+1)|i,j,k,0}+aj×{dat(s)|i,j+1,k,0+dat(s+2)|i,j-1,k,0-2×dat(s+1)|i,j,k,0}+ak×{dat(s+1)|i,j,k+1,0+dat(s+1)|i,j,k-1,0-2×dat(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述dat(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，所述dat(s)|i,j+1,k,0为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据，所述dat(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，所述dat(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条线数据的迭代结果：
dat(s+1)|i,j,k,1＝dat(s+1)|i,j,k,0+ai×{dat(s+1)|i+1,j,k,0+dat(s+1)|i-1,j,k,0-2×dat(s+1)|i,j,k,0}+aj×{dat(s)|i,j+1,k,0+dat(s+2)|i,j-1,k,0-2×dat(s+1)|i,j,k,0}+ak×{dat(s+1)|i,j,k+1,0+dat(s+1)|i,j,k-1,0-2×dat(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述dat(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，所述dat(s)|i,j+1,k,0为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据，所述dat(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，所述dat(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条线数据的迭代结果：
dat(s+1)|i,j,k,1＝dat(s+1)|i,j,k,0+ai×{dat(s+1)|i+1,j,k,0+dat(s+1)|i-1,j,k,0-2×dat(s+1)|i,j,k,0}+aj×{dat(s)|i,j+1,k,0+dat(s+2)|i,j-1,k,0-2×dat(s+1)|i,j,k,0}+ak×{dat(s+1)|i,j,k+1,0+dat(s+1)|i,j,k-1,0-2×dat(s+1)|i,j,k,0}
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，dat(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，dat(s)|i,j+1,k,0为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据，dat(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k,1为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据的迭代结果，ai为i对应的坐标轴方向上的预设迭代系数，aj为j对应的坐标轴方向上的预设迭代系数，ak为k对应的坐标轴方向上的预设迭代系数。
上述所有可选技术方案，均可按照任意结合形成本发明的可选实施例，本发明实施例对此不再一一赘述。
图2是本发明实施例提供的另一种地震波数据的低波数噪音分离方法的流程示意图，该方法可以应用于终端中，该终端可以为手机、平板电脑或计算机等。参见图2，该方法包括如下步骤：
步骤201：获取目标工区的第一深度域地震波数据体。
需要说明的是，目标工区是指目标工作区域，具体可以为油藏上的任一待处理的工作区域。而且，目标工区包括多个样点，这多个样点是在目标工区上预先设置的样点，用于对目标工区的空间位置进行定位，且多个样点中的每个样点均可以存储目标工区在该样点位置处的深度域地震波数据。
需要说明的是，第一深度域地震波数据体是指对目标工区的原始深度域地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的深度域地震波数据体。在实际应用中，第一深度域地震波数据体可以由用户输入得到，也可以由其他设备发送得到，还可以先获取原始深度域地震波数据体，并通过逆时成像技术对原始深度域地震波数据体进行处理，得到带有低波数噪音的深度域地震波数据体，本发明实施例对此不做具体限定。另外，由于逆时成像技术是对叠前深度域地震波数据体进行处理，因此，得到的第一深度域地震波数据体也属于叠前深度域地震波数据体。
此外，由于逆时成像技术在偏移成像路径上是将不相关的炮、检点波场进行了互相关处理，因此，得到的第一深度域地震波数据体会引入低波数噪音，而在后续应用中由于第一深度域地震波数据体中存在低波数噪音，使第一深度域地震波数据体的处理结果很难实现保幅，从而影响了根据第一深度域地震波数据体确定目标工区的地质构造的结果。因此，可以对目标工区的第一深度域地震波数据体进行低波数噪音分离，以使第一深度域地震波数据体的处理结果得以保幅，提高根据第一深度域地震波数据体确定地震成像图像的精度。
步骤202：按照预设规则，对第一深度域地震波数据体进行数据补充，得到目标工区的第二深度域地震波数据体。
需要说明的是，在应用本发明实施例提供的地震波数据的低波数噪音分离方法时，为了便于在后续步骤中对地震波数据进行处理，可以先对第一深度域地震波数据体对应的目标工区中的多个样点进行补充，使补充后的多个样点在立体空间中排列时，组成的空间为规则的立方体，也即是该立体空间的边界为光滑的平面，并且没有分布在立体空间外的样点。例如，由于原始叠前地震波数据体对应的目标工区中的多个样点在立体空间中排列时，组成的空间可能是三棱柱或三棱锥等空间，因此，对原始叠前地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的深度域地震波数据体对应的目标工区中的多个样点在立体空间中排列时，组成的空间也是三棱柱或三棱锥等空间，这种情况下，由于本发明实施例的后续步骤需要在样点排列为规则的立方体空间中进行数据处理，便可对目标工区进行样点补充，使补充后的目标工区中的多个样点组成的立体空间为立方体，并对补充的多个样点进行地震波数据赋值，最后，将补充的立方体空间内多个样点中存储的地震波数据确定为第一深度域地震波数据。
也即是，通过对第一深度域地震波数据体进行数据补充，可以在第一深度域地震波数据体对应的目标工区中增加更多样点，使得增加样点后的目标工区中的多个样点在立体空间中排列时，组成的空间为立方体。并且，第二深度域地震波数据体中即包括第一深度域地震波数据体中的多个样点，也可包括增加的样点，且所增加的每个样点均具有对应的地震波数据。可选地，所增加的每个样点对应的地震波数据可以为零，当然，也可以为其他数据，本发明实施例对此不做限定。
具体地，对第一深度域地震波数据体进行数据补充，以得到第二深度域地震波数据体可以包括：确定多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；根据多个样点在预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，目标空间网格的外形为立方体，且目标空间网格的网格点包括多个样点；将目标网格点对应的深度域地震波数据赋值为零，目标网格点是指目标空间网格中除多个样点之外的网格点；将赋值后的目标空间网格中每个网格点对应的深度域地震波数据，确定为第二深度域地震波数据体。
需要说明的是，为了对第一深度域地震波数据体对应的目标工区中的多个样点在立体空间中进行准确定位，在对第一深度域地震波数据体进行数据补充前，还可以在第一深度域地震波数据体中设立以XYZ为坐标轴、以O为原点的预设空间直角坐标系。
具体地，可以在目标工区中，将立体空间中的一个边界线上的端点位置处的样点设置为预设空间直角坐标系的原点O，将立体空间中原点O所在的边界面上的互相垂直的两个方向分别设置为预设空间直角坐标系的X轴和Y轴，将垂直于立体空间中原点O所在的边界面的方向设置为预设空间直角坐标系的Z轴。例如，在第一深度域地震波数据体中，用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之间的样点对应的网格点确定为目标网格点，并将该目标网格点对应的深度域地震波数据赋值为零，进而将赋值后的目标空间网格中每个网格点对应的深度域地震波数据，确定为第二深度域地震波数据体。
步骤203：通过热传导变换方程，对第二深度域地震波数据体进行迭代处理，得到目标工区的低波数噪音深度域地震波数据体。
需要说明的是，热传导变换方程是对热传导方法中应用到描述区域内的温度随时间变化的偏微分方程进行变换，使该热传导变换方程可以用来描述目标网格中地震波数据随时间变化的偏微分方程。
具体的，热传导方程为由下述公式(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轴方向的迭代系数。
需要说明的是，在根据第一深度域地震波数据体确定第二深度域地震波数据体时，可以对第一深度域地震波数据体进行N次迭代，N为预设的迭代次数，通过预设N的数值，来控制迭代处理后的得到的第二深度域地震波数据体的准确性，预设迭代次数N的数值越高，根据第一深度域地震波数据体确定第二深度域地震波数据体时的迭代次数也就越多，得到的第二深度域地震波数据体更加准确，也更符合目标工区的真实地质构造。
另外，在迭代过程中，还可以设置判断步骤，判断迭代次数是否达到了预设迭代次数，如果迭代次数未达到预设迭代次数，则继续执行迭代步骤，如果迭代次数达到预设迭代次数，则停止迭代，并将迭代结果确定为第二深度域地震波数据体。
具体地，通过热传导变换方程，可以通过以下步骤对第二深度域地震波数据体进行迭代处理，得到目标工区的低波数噪音深度域地震波数据体。
步骤2031：令r＝1，将第一深度域地震波数据体确定为第r个地震波数据体，将第r个地震波数据体确定为目标地震波数据体。
需要说明的是，在对第一深度域地震波数据体进行N次迭代计算时，由于每次迭代过程均相同，因此，可以设置一个目标地震波数据体，将每次迭代开始前的深度域地震波数据体确定为目标地震波数据体，再对目标地震数据体进行迭代计算。在第一次迭代过程中，可以令r＝1，将第一深度域地震波数据体确定为第r个地震波数据体，再将第r个地震波数据体确定为目标地震波数据体，从而对目标地震波数据体执行迭代计算。
步骤2032：通过热传导变换方程，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体。
具体地，可以通过下述迭代步骤1)-3)，对目标地震波数据体进行迭代处理，得到第r+1个地震波数据体。
1)令s＝2，获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的第s-1、s、s+1条线数据，每条线数据包括在坐标轴方向上的坐标值相同的至少一个样点所对应的深度域地震波数据，利用热传导方程的离散函数，对第s-1、s、s+1条线数据进行迭代处理，得到第s条线数据的迭代结果。
需要说明的是，在通过热传导变换方程对目标地震波数据体进行迭代处理时，目标地震波数据体的边界位置处的线数据是无法进行迭代的，只能对除目标地震波数据体的边界位置处的线数据之外的线数据进行迭代处理。因此，在第一次获取目标地震波数据体中在预设空间直角坐标系的任一坐标轴方向上的三条线数据时，是从边界位置处开始按顺序进行获取三条线数据，而获取之后只能根据三条线数据和热传导变换方程可以对三条线数据中的第二条线数据进行一次迭代。例如，可以从X轴方向上坐标最小值对应的inline线数据开始，从小到大依次获取三条inline线数据，记为dat1、dat2和dat3，利用公式(2)示出的热传导变换方程对dat1、dat2和dat3进行迭代处理，便可得到dat2的迭代结果，并将dat2的迭代结果记为dat21。
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条线数据的迭代结果：
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，dat(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，dat(s)|i,j+1,k,1为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据的迭代结果，dat(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，dat(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)进行迭代处理。
例如，在对dat1、dat2和dat3通过公式(2)进行一次迭代处理后得到dat21，将dat1存入临时中间文件，用dat21替换dat1，用dat3替换dat2，读取dat4，并用dat4替换dat3，这样，对于线数据dat21、dat3和dat4，取迭代系数为0.1，则根据公式(3)可以对dat3在(2，2，2)坐标处的样点进行迭代处理：
dat(3)|2,2,2,1＝dat(3)|2,2,2,0+0.1×{dat(3)|3,2,2,0+dat(3)|1,2,2,0-2×dat(4)|2,2,2,0}+0.1×{dat(2)|2,3,2,1+dat(4)|2,1,2,0-2×dat(3)|2,2,2,0}+0.1×{dat(3)|2,2,3,0+dat(3)|2,2,1,0-2×dat(3)|2,2,2,0}
在对dat3在(2，2，2)坐标处的样点进行迭代处理后，还可以通过公式(3)对dat3在除(2，2，2)坐标之外的其他位置的样点所存储的地震波数据进行迭代处理，在完成对lin3上的除边界位置上的样点之外的其他所有样点的数据处理后，可以得到dat31。
需要说明的是，上述样点坐标、线数据条数、迭代系数均为示例性数据，实际应用中还可以为其他数据，本发明实施例对此不做具体限定。
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个地震波数据体。
步骤2033：判断r+1是否等于N，N为大于或等于3的整数。
步骤2034：若r+1不等于N，则令r＝r+1，并返回步骤2032。
也即是，令r＝r+1，并返回将第r个地震波数据体确定为目标地震波数据体的步骤，如此，可以继续对第r+1个地震波数据体进行迭代处理，得到第r+2个地震波数据体。
步骤2035：若r+1等于N，则将第r+1个地震波数据体确定为低波数噪音深度域地震波数据体。
需要说明的是，在根据第r个地震波数据体确定了第r+1个地震波数据体之后，为了进一步提高数据分离的准确性，使分离之后得到的低波数噪音深度域地震波数据体的精度更高，还可以对第r+1个地震波数据体按照步骤2032进行N次迭代处理。其中，N是预设的迭代次数，可以根据N的取值控制迭代后得到的不同精度的深度域反射波数据体。在实际应用中，N可以由用户输入得到，可以由其他设备发送得到，本发明实施例对此不做具体限定。例如，可以由用户在终端中输入N的数值，如50、80或100等，N的数值越大，则迭代处理过程所需要的时间越长，得到的低波数噪音深度域地震波数据体的精度也越高，本发明实施例对N的数值不做具体限定。
另外，为了将低波数噪音从第二深度域地震波数据体中的地震波数据中分离的更彻底，可以将N设置为大于或等于3的整数，也即是，对第二深度域地震波数据体按步骤2032所示出的方法，迭代至少两次，得到低波数噪音深度域地震波数据体。
此外，在判断r+1与N的关系时，当r+1不等于N时，则使r＝r+1，并重新根据将第r个地震波数据体确定目标地震波数据体的步骤确定新的目标地震波数据体，当r+1等于N时，则可以将第r+1个地震波数据体确定为低波数噪音深度域地震波数据体。
步骤204：将第二深度域地震波数据体和低波数噪音深度域地震波数据体之间的差值，确定为目标工区的去低波数噪音深度域地震波数据体。
需要说明的是，由于低波数噪音深度域地震波数据体中的每个样点所对应的地震波数据是对第二深度域地震波数据体中的每个样点所对应的地震波数据进行迭代处理后得到的，因此，低波数噪音深度域地震波数据体中的每个样点所对应的地震波数据和第二深度域地震波数据体中的每个样点所对应的地震波数据的数据类型是一致的，两者中的样点位置也是一致的。基于此，在确定了低波数噪音深度域地震波数据体之后，可以直接将第二深度域地震波数据体和低波数噪音深度域地震波数据体进行相减，并将第二深度域地震波数据体和低波数噪音深度域地震波数据体之间的差值，确定为目标工区的去低波数噪音深度域地震波数据体。
例如，在目标工区中，将第二深度域地震波数据体记为DATA1，第二深度域地震波数据体成像图如图3所示，利用本发明实施例提供的地震波数据的低波数噪音分离方法对第二深度域地震波数据体进行低波数噪音分离，得到低波数噪音深度域地震波数据体和去低波数噪音深度域地震波数据体，将低波数噪音深度域地震波数据体记为DATA2，低波数噪音深度域地震波数据体主要表现为对原始深度域地震波数据体进行逆时成像处理后得到的第一深度域地震波数据体中包含的低波数噪音，低波数噪音深度域地震波数据体成像图如图4所示，将第二深度域地震波数据体DATA1和低波数噪音深度域地震波数据体DATA2进行相减，得到去低波数噪音深度域地震波数据体将去低波数噪音深度域地震波数据体记为DATA3，则DATA3即为目标工区的去除去低波数噪音后的深度域地震波数据体，去低波数噪音地震波数据体成像图如图5所示，去低波数噪音成像剖面即为将低波数噪音分离后得到的深度域地震波数据体成像剖面。
在本发明实施例中，可以获取对目标工区的原始深度域地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的第一深度域地震波数据体，并按照预设规则，对第一深度域地震波数据体进行数据补充，得到目标工区的第二深度域地震波数据体，通过基于热传导方法推导得到，用于描述目标工区内的深度域地震波数据随时间的变化情况的热传导变换方程，对第二深度域地震波数据体进行迭代处理，得到目标工区的低波数噪音深度域地震波数据体，再用第二深度域地震波数据体减去低波数噪音深度域地震波数据体，便可得到目标工区的不含低波数噪音的去低波数噪音深度域地震波数据体，从而实现了将低波数噪音从采用逆时成像技术处理后得到的第一深度域地震波数据体中分离的目的。由于本发明实施例可以直接对采用逆时成像技术处理后得到的深度域地震波数据体进行低波数噪音分离，解决了相关技术中对深度域地震波数据体进行低波数噪音分离时，需要先将深度域地震波数据体转换为时间域地震波数据体，并对时间域地震波数据体进行低波数噪音分离，得到去低波数噪音时间域地震波数据体，最后将去低波数噪音时间域地震波数据体转换为去低波数噪音深度域地震波数据体，简化了深度域地震波数据体进行低波数噪音分离的步骤，提高了对地震波数据体中的地震波数据的低波数噪音分离的效率。并且，由于热传导变换方程是对地震波数据体中的地震波数据整体滤波，使得地震波数据体中的地震波数据得以保幅，提高了地震成像图像的精度。
图6是本发明实施例提供的一种地震波数据的低波数噪音分离装置的流程示意图。参见图6，该装置可以包括：
获取模块601，用于获取目标工区的第一深度域地震波数据体，第一深度域地震波数据体是指对目标工区的原始深度域地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的深度域地震波数据体。
补充模块602，用于按照预设规则，对第一深度域地震波数据体进行数据补充，得到目标工区的第二深度域地震波数据体。
处理模块603，用于通过热传导变换方程，对第二深度域地震波数据体进行迭代处理，得到目标工区的低波数噪音深度域地震波数据体，热传导变换方程是基于热传导方法推导得到，用于描述目标工区内的深度域地震波数据随时间的变化情况。
确定模块604，用于将第二深度域地震波数据体和低波数噪音深度域地震波数据体之间的差值，确定为目标工区的去低波数噪音深度域地震波数据体。
可选地，补充模块包括：
第一确定单元，用于确定多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
第二确定单元，用于根据多个样点在预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，目标空间网格的外形为立方体，且目标空间网格的网格点包括多个样点；
赋值单元，用于将目标网格点对应的深度域地震波数据赋值为零，目标网格点是指目标空间网格中除多个样点之外的网格点；
第三确定单元，用于将赋值后的目标空间网格中每个网格点对应的深度域地震波数据，确定为第二深度域地震波数据体。
可选地，处理模块包括：
第四确定单元，用于令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条线数据的迭代结果：
dat(s+1)|i,j,k,1＝dat(s+1)|i,j,k,0+ai×{dat(s+1)|i+1,j,k,0+dat(s+1)|i-1,j,k,0-2×dat(s+1)|i,j,k,0}+aj×{dat(s)|i,j+1,k,0+dat(s+2)|i,j-1,k,0-2×dat(s+1)|i,j,k,0}+ak×{dat(s+1)|i,j,k+1,0+dat(s+1)|i,j,k-1,0-2×dat(s+1)|i,j,k,0}
其中，i，j，k分别与预设空间直角坐标系的3个坐标轴方向对应，dat(s+1)|i,j,k,0为第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i+1,j,k,0为第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，dat(s+1)|i-1,j,k,0为第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，dat(s)|i,j+1,k,0为第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据，dat(s+2)|i,j-1,k,0为第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k+1,0为第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，dat(s+1)|i,j,k-1,0为第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，dat(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 (10)
1.一种地震波数据的低波数噪音分离方法，其特征在于，所述方法包括：
获取目标工区的第一深度域地震波数据体，所述第一深度域地震波数据体是指对所述目标工区的原始深度域地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的深度域地震波数据体；
按照预设规则，对所述第一深度域地震波数据体进行数据补充，得到所述目标工区的第二深度域地震波数据体；
通过热传导变换方程，对所述第二深度域地震波数据体进行迭代处理，得到所述目标工区的低波数噪音深度域地震波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域地震波数据随时间的变化情况；
将所述第二深度域地震波数据体和所述低波数噪音深度域地震波数据体之间的差值，确定为所述目标工区的去低波数噪音深度域地震波数据体。
2.根据权利要求1所述的方法，其特征在于，所述第一深度域地震波数据体包括所述目标工区中多个样点对应的深度域地震波数据，所述按照预设规则，对所述第一深度域地震波数据体进行数据补充，得到所述目标工区的第二深度域地震波数据体，包括：
确定所述多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
将目标网格点对应的深度域地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
将赋值后的所述目标空间网格中每个网格点对应的深度域地震波数据，确定为所述第二深度域地震波数据体。
3.根据权利要求1所述的方法，其特征在于，所述通过热传导变换方程，对所述第二深度域地震波数据体进行迭代处理，得到所述目标工区的低波数噪音深度域地震波数据体，包括：
令r＝1，将所述第一深度域地震波数据体确定为第r个地震波数据体；
将所述第r个地震波数据体确定为目标地震波数据体；
通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数；
若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述低波数噪音深度域地震波数据体。
4.根据权利要求3所述的方法，其特征在于，所述通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第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个地震波数据体。
5.根据权利要求4所述的方法，其特征在于，所述利用热传导方程的离散函数，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果，包括：
通过以下公式，对所述第s条线数据的迭代结果、所述第s+1条线数据、所述第s+2条线数据进行迭代处理，得到所述第s+1条线数据的迭代结果：
dat(s+1)|i,j,k,1＝dat(s+1)|i,j,k,0+ai×{dat(s+1)|i+1,j,k,0+dat(s+1)|i-1,j,k,0-2×dat(s+1)|i,j,k,0}
+aj×{dat(s)|i,j+1,k,0+dat(s+2)|i,j-1,k,0-2×dat(s+1)|i,j,k,0}
+ak×{dat(s+1)|i,j,k+1,0+dat(s+1)|i,j,k-1,0-2×dat(s+1)|i,j,k,0}
其中，i，j，k分别与所述预设空间直角坐标系的3个坐标轴方向对应，所述dat(s+1)|i,j,k,0为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i+1,j,k,0为所述第s+1条线数据中在坐标(i+1，j，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i-1,j,k,0为所述第s+1条线数据中在坐标(i-1，j，k)上的样点对应的深度域地震波数据，所述dat(s)|i,j+1,k,0为所述第s条线数据中在坐标(i，j+1，k)上的样点对应的深度域地震波数据，所述dat(s+2)|i,j-1,k,0为所述第s+2条线数据中在坐标(i，j-1，k)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k+1,0为所述第s+1条线数据中在坐标(i，j，k+1)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k-1,0为所述第s+1条线数据中在坐标(i，j，k-1)上的样点对应的深度域地震波数据，所述dat(s+1)|i,j,k,1为所述第s+1条线数据中在坐标(i，j，k)上的样点对应的深度域地震波数据的迭代结果，所述ai为所述i对应的坐标轴方向上的预设迭代系数，所述aj为所述j对应的坐标轴方向上的预设迭代系数，所述ak为所述k对应的坐标轴方向上的预设迭代系数。
6.一种地震波数据的低波数噪音分离装置，其特征在于，所述装置包括：
获取模块，用于获取目标工区的第一深度域地震波数据体，所述第一深度域地震波数据体是指对所述目标工区的原始深度域地震波数据体采用逆时成像技术处理后得到的带有低波数噪音的深度域地震波数据体；
补充模块，用于按照预设规则，对所述第一深度域地震波数据体进行数据补充，得到所述目标工区的第二深度域地震波数据体；
处理模块，用于通过热传导变换方程，对所述第二深度域地震波数据体进行迭代处理，得到所述目标工区的低波数噪音深度域地震波数据体，所述热传导变换方程是基于热传导方法推导得到，用于描述所述目标工区内的深度域地震波数据随时间的变化情况；
确定模块，用于将所述第二深度域地震波数据体和所述低波数噪音深度域地震波数据体之间的差值，确定为所述目标工区的去低波数噪音深度域地震波数据体。
7.根据权利要求6所述的方法，其特征在于，所述补充模块包括：
第一确定单元，用于确定所述多个样点在预设空间直角坐标系中的每个方向上的最大坐标值和最小坐标值；
第二确定单元，用于根据所述多个样点在所述预设空间直角坐标系中的每个坐标轴方向上的最大坐标值和最小坐标值，确定目标空间网格，所述目标空间网格的外形为立方体，且所述目标空间网格的网格点包括所述多个样点；
赋值单元，用于将目标网格点对应的深度域地震波数据赋值为零，所述目标网格点是指所述目标空间网格中除所述多个样点之外的网格点；
第三确定单元，用于将赋值后的所述目标空间网格中每个网格点对应的深度域地震波数据，确定为所述第二深度域地震波数据体。
8.根据权利要求6所述的方法，其特征在于，所述处理模块包括：
第四确定单元，用于令r＝1，将所述第一深度域地震波数据体确定为第r个地震波数据体；
第五确定单元，用于将所述第r个地震波数据体确定为目标地震波数据体；
处理单元，用于通过所述热传导变换方程，对所述目标地震波数据体进行迭代处理，得到第r+1个地震波数据体，并判断所述r+1是否等于N，所述N为大于或等于3的整数；
第一判断单元，用于若所述r+1不等于所述N，则令r＝r+1，并返回所述将所述第r个地震波数据体确定为目标地震波数据体的步骤；
第二判断单元，用于若所述r+1等于所述N，则将所述第r+1个地震波数据体确定为所述低波数噪音深度域地震波数据体。
9.一种地震波数据的低波数噪音分离装置，其特征在于，所述装置包括：
处理器及用于存储处理器可执行指令的存储器；
其中，所述处理器被配置为执行如权利要求1-5任一项所述的方法。
10.一种计算机可读存储介质，其特征在于，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-5任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811184185.0A CN109116424B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811184185.0A CN109116424B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109116424A true CN109116424A (zh) | 2019-01-01 |
| CN109116424B CN109116424B (zh) | 2020-06-09 |
Family
ID=64857854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811184185.0A Active CN109116424B (zh) | 2018-10-11 | 2018-10-11 | 地震波数据的低波数噪音分离方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109116424B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113031070A (zh) * | 2021-03-19 | 2021-06-25 | 大庆油田有限责任公司 | 一种深度域合成地震记录的制作方法 |
| CN113534258A (zh) * | 2020-04-21 | 2021-10-22 | 中国石油天然气股份有限公司 | 地震数据指纹特征处理方法及装置 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103076626A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种波场净化处理方法 |
| CN103489163A (zh) * | 2013-09-13 | 2014-01-01 | 电子科技大学 | 基于正则化混合范数滤波的地震图像结构导向降噪方法 |
| CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
| CN106597549A (zh) * | 2016-12-05 | 2017-04-26 | 中国石油大学(北京) | 一种反射波提取方法及装置、反射声波成像测井仪器 |
| WO2017149419A1 (en) * | 2016-03-04 | 2017-09-08 | Seismic Apparition Gmbh | Method for deghosting and redatuming operator estimation |
| CN107247290A (zh) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | 一种基于时空分数阶滤波的地震资料噪声压制方法 |
| CN107422381A (zh) * | 2017-09-18 | 2017-12-01 | 西南石油大学 | 一种基于eemd‑ica的地震低频信息流体预测方法 |
| CN108181657A (zh) * | 2017-12-28 | 2018-06-19 | 中国石油大学(北京) | 全波形反演梯度计算中分离偏移和层析成像模式的方法 |
| US20180196154A1 (en) * | 2017-01-12 | 2018-07-12 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
-
2018
- 2018-10-11 CN CN201811184185.0A patent/CN109116424B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103076626A (zh) * | 2011-10-25 | 2013-05-01 | 中国石油化工股份有限公司 | 一种波场净化处理方法 |
| CN103489163A (zh) * | 2013-09-13 | 2014-01-01 | 电子科技大学 | 基于正则化混合范数滤波的地震图像结构导向降噪方法 |
| CN105044771A (zh) * | 2015-08-05 | 2015-11-11 | 北京多分量地震技术研究院 | 基于有限差分法的三维tti双相介质地震波场数值模拟方法 |
| WO2017149419A1 (en) * | 2016-03-04 | 2017-09-08 | Seismic Apparition Gmbh | Method for deghosting and redatuming operator estimation |
| CN106597549A (zh) * | 2016-12-05 | 2017-04-26 | 中国石油大学(北京) | 一种反射波提取方法及装置、反射声波成像测井仪器 |
| US20180196154A1 (en) * | 2017-01-12 | 2018-07-12 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
| CN107247290A (zh) * | 2017-07-06 | 2017-10-13 | 西安交通大学 | 一种基于时空分数阶滤波的地震资料噪声压制方法 |
| CN107422381A (zh) * | 2017-09-18 | 2017-12-01 | 西南石油大学 | 一种基于eemd‑ica的地震低频信息流体预测方法 |
| CN108181657A (zh) * | 2017-12-28 | 2018-06-19 | 中国石油大学(北京) | 全波形反演梯度计算中分离偏移和层析成像模式的方法 |
Non-Patent Citations (2)
| Title |
|---|
| 陈可洋 等: "基于扩散滤波的多尺度分解和重构方法及应用初探", 《计算物理》 * |
| 陈可洋: "一种相对保幅的低频逆时噪声压制方法及其应用", 《油气藏评价与开发》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113534258A (zh) * | 2020-04-21 | 2021-10-22 | 中国石油天然气股份有限公司 | 地震数据指纹特征处理方法及装置 |
| CN113534258B (zh) * | 2020-04-21 | 2024-02-02 | 中国石油天然气股份有限公司 | 地震数据指纹特征处理方法及装置 |
| CN113031070A (zh) * | 2021-03-19 | 2021-06-25 | 大庆油田有限责任公司 | 一种深度域合成地震记录的制作方法 |
| CN113031070B (zh) * | 2021-03-19 | 2022-06-28 | 大庆油田有限责任公司 | 一种深度域合成地震记录的制作方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109116424B (zh) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110083791A (zh) | 目标群组检测方法、装置、计算机设备及存储介质 | |
| CN109308727A (zh) | 虚拟形象模型生成方法、装置及存储介质 | |
| CN109977775A (zh) | 关键点检测方法、装置、设备及可读存储介质 | |
| CN109887494A (zh) | 重构语音信号的方法和装置 | |
| CN109816042A (zh) | 数据分类模型训练的方法、装置、电子设备和存储介质 | |
| CN110064200A (zh) | 基于虚拟环境的物体构建方法、装置及可读存储介质 | |
| CN110047152A (zh) | 基于虚拟环境的物体构建方法、装置及可读存储介质 | |
| CN109558837A (zh) | 人脸关键点检测方法、装置及存储介质 | |
| CN109696706A (zh) | 盐构造检测方法、装置及存储介质 | |
| CN109886208A (zh) | 物体检测的方法、装置、计算机设备及存储介质 | |
| CN109166150A (zh) | 获取位姿的方法、装置存储介质 | |
| CN109192218A (zh) | 音频处理的方法和装置 | |
| CN110174699A (zh) | 确定断层封闭性的方法、装置及存储介质 | |
| CN110389380A (zh) | 地震剖面的同相轴自动拾取方法、装置及存储介质 | |
| CN109547843A (zh) | 对音视频进行处理的方法和装置 | |
| CN109116424A (zh) | 地震波数据的低波数噪音分离方法、装置及存储介质 | |
| CN109992685A (zh) | 一种检索图片的方法及装置 | |
| CN109218751A (zh) | 推荐音频的方法、装置及系统 | |
| CN109102811A (zh) | 音频指纹的生成方法、装置及存储介质 | |
| CN110097002A (zh) | 形近字确定方法、装置、计算机设备和存储介质 | |
| CN109189290A (zh) | 点击区域识别方法、装置和计算机可读存储介质 | |
| CN109660876A (zh) | 显示列表的方法和装置 | |
| CN110166275A (zh) | 信息处理方法、装置及存储介质 | |
| CN107884823B (zh) | 砂体连通图的绘制方法及装置 | |
| CN109833623A (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 |