CN110824555B - 地震能量均衡方法及装置、计算机可读存储介质 - Google Patents
地震能量均衡方法及装置、计算机可读存储介质 Download PDFInfo
- Publication number
- CN110824555B CN110824555B CN201810907206.0A CN201810907206A CN110824555B CN 110824555 B CN110824555 B CN 110824555B CN 201810907206 A CN201810907206 A CN 201810907206A CN 110824555 B CN110824555 B CN 110824555B
- Authority
- CN
- China
- Prior art keywords
- point data
- sampling point
- seismic
- energy value
- data
- 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 54
- 238000005070 sampling Methods 0.000 claims abstract description 235
- 239000000523 sample Substances 0.000 claims abstract description 71
- 238000004458 analytical method Methods 0.000 claims abstract description 56
- 238000001514 detection method Methods 0.000 claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000007781 pre-processing Methods 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000036961 partial effect Effects 0.000 description 14
- 239000011159 matrix material Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer Methods 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
技术领域
本发明涉及地震勘探技术领域，特别涉及一种地震能量均衡方法及装置、计算机可读存储介质。
背景技术
地震勘探技术是利用人工激发的地震波在地层内传播来勘探地下的地质情况的技术。在地面某处激发的地震波向地下传播时，遇到不同弹性的地层分界面就会产生反射波或折射波返回地面，用专门的仪器可记录这些波，以获取到地震数据，通过分析所记录的波的特点(如波的传播时间、振动形状等)，能较准确地测定这些界面的深度和形态，判断地层的岩性。
地震波的振幅随着传播距离(时间)的增大而不断衰减，且振幅的衰减还随频率的不同而变化，高频成份比低频成份衰减得更快，这使得地震数据中主频的频率降低，导致地震数据分辨率降低。为了提高地震数据的分辨率，需要对地震数据进行地震能量均衡处理。
目前在进行地震能量均衡时，通常是对地震数据中的多个地震道进行逐个处理，但是相邻地震道之间还是存在能量不均，同相轴不连续的问题，导致地震数据细节不清晰。
发明内容
本发明实施例提供了一种地震能量均衡方法及装置、计算机可读存储介质，能够解决现有地震能量均衡处理存在相邻地震道之间能量不均，同相轴不连续的问题。所述技术方案如下：
一方面，本发明实施例提供了一种地震能量均衡方法，所述方法包括：
获取地震数据，所述地震数据包括由空间上矩阵分布的多个检波点获取的多个地震道数据，所述多个检波点的采样起止时间和采样间隔均相同，每个所述地震道数据均包括对应的所述检波点沿时间顺序采集的多个采样点数据的集合；
计算当前中心采样点数据对应的均方根能量值，所述均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，所述多个检波点沿所述矩形窗口的长度方向和宽度方向排列，在所述矩形窗口内沿所述矩形窗口的长度和宽度方向均排列有奇数个检波点，每个所述检波点在所述预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，所述当前中心采样点数据为当前位于所述矩形窗口的中心的检波点在所述预设分析时窗的中间时刻获取的采样点数据；
计算所述当前中心采样点数据对应的均衡振幅能量值，所述均衡振幅能量值为所述当前中心采样点数据对应的所述振幅能量值和所述均方根能量值的商；
移动所述矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出所述地震数据中的所有中心采样点数据对应的均衡振幅能量值。
可选地，所述多个检波点共直线分布，所述矩形窗口内沿所述矩形窗口的长度方向分布有至少3个检波点，沿所述矩形窗口的宽度方向分布有1个检波点。
可选地，所述多个检波点呈网格分布，所述矩形窗口内沿所述矩形窗口的长度方向和宽度方向均至少分布有3个检波点。
可选地，所述方法还包括：
以边缘采样点数据对应的振幅能量值与参考能量值的商作为所述边缘采样点数据对应的均衡振幅能量值，所述边缘采样点数据为所有的所述采样点数据中除去所有中心采样点数据之外的采样点数据，所述参考能量值为与所述边缘采样点数据的采集时间相隔最近的中心采样点数据中，采集于与所述边缘采样点数据所对应的检波点相距最近的检波点的中心采样点数据的均方根能量值，或者与所述边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。
可选地，所述方法还包括：
在所述获取地震数据之后，且在所述计算当前中心采样点数据对应的均方根能量值之前，对所述地震数据进行预处理，所述预处理包括噪声压制、动补偿、静补偿中的至少一种。
另一方面，本发明实施例提供了一种地震能量均衡装置，所述地震能量均衡装置包括：
获取模块，用于获取地震数据，所述地震数据包括由空间上矩阵分布的多个检波点获取的多个地震道数据，所述多个检波点的采样起止时间和采样间隔均相同，每个所述地震道数据均包括对应的所述检波点沿时间顺序采集的多个采样点数据的集合；
处理模块，用于计算当前中心采样点数据对应的均方根能量值，计算所述当前中心采样点数据对应的均衡振幅能量值，且还用于移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出所述地震数据中的所有中心采样点数据对应的均衡振幅能量值，所述均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，所述多个检波点沿所述矩形窗口的长度方向和宽度方向排列，在所述矩形窗口内沿所述矩形窗口的长度和宽度方向均排列有奇数个检波点，每个所述检波点在所述预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，所述当前中心采样点数据为当前位于所述矩形窗口的中心的检波点在所述预设分析时窗的中间时刻获取的采样点数据，所述均衡振幅能量值为所述当前中心采样点数据对应的所述振幅能量值和所述均方根能量值的商。
可选地，所述多个检波点共直线分布，所述矩形窗口内沿所述矩形窗口的长度方向分布有至少3个检波点，沿所述矩形窗口的宽度方向分布有1个检波点。
可选地，所述多个检波点呈网格分布，所述矩形窗口内沿所述矩形窗口的长度方向和宽度方向均至少分布有3个检波点。
可选地，所述处理模块还用于以边缘采样点数据对应的振幅能量值与参考能量值的商作为所述边缘采样点数据对应的均衡振幅能量值，所述边缘采样点数据为所有的所述采样点数据中除去所有中心采样点数据之外的采样点数据，所述参考能量值为与所述边缘采样点数据的采集时间相隔最近的中心采样点数据中，采集于与所述边缘采样点数据所对应的检波点相距最近的检波点的中心采样点数据的均方根能量值，或者与所述边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。
又一方面，本发明实施例还提供了一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现如前所述的地震能量均衡方法。
本发明实施例提供的技术方案带来的有益效果至少包括：通过获取地震数据后，以当前中心采样点数据对应的振幅能量值与当前中心采样点数据对应的均方根能量值的商作为当前中心采样点数据对应的均衡振幅能量值，由于均方根能量值为矩形窗口内的所有检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，而矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，使得均方根能量值的计算不仅和与当前中心采样点数据属于同一地震道数据的采样点数据有关，也与相邻的地震道数据中的采样点数据有关，通过移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值，由于均方根振幅能量是采用相同大小的矩形窗口和相同大小的预设分析时窗计算的，其值是连续变化的，因此在进行地震能量均衡处理之后，可以使相邻地震道之间的能量更均匀，同向轴连续，使地震数据的细节显示的更加清晰。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种地震能量均衡方法流程图；
图2是本发明实施例提供的另一种地震能量均衡方法流程图；
图3是本发明实施例提供的地震数据的局部示意图；
图4是本发明实施例提供的另一种地震数据的局部示意图；
图5是图3所示的地震数据采用现有的地震能量均衡方法处理后的局部示意图；
图6是图3所示的地震数据采用图2所示的地震能量均衡方法处理后的局部示意图；
图7是采用现有的地震能量均衡方法处理后的地面地震数据的局部示意图；
图8是采用图2所示的地震能量均衡方法处理后的地面地震数据的局部示意图；
图9是本发明实施例提供的一种地震能量均衡装置的结构框图；
图10示出了本发明一个示例性实施例提供的计算机设备的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种地震能量均衡方法流程图。如图1所示，该方法包括：
S11：获取地震数据。
其中，地震数据包括由空间上矩阵分布的多个检波点获取的多个地震道数据，多个检波点的采样起止时间和采样间隔均相同，每个地震道数据均包括对应的检波点沿时间顺序采集的多个采样点数据的集合。
可选地，地震数据可以是道集数据，例如共炮点道集数据，也可以是叠加剖面数据，例如偏移叠加剖面数据、水平叠加剖面数据。
S12：计算当前中心采样点数据对应的均方根能量值。
其中，均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根。
多个检波点沿矩形窗口的长度方向和宽度方向排列，在矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，每个检波点在预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，当前中心采样点数据为当前位于矩形窗口的中心的检波点在预设分析时窗的中间时刻获取的采样点数据。
S13：计算当前中心采样点数据对应的均衡振幅能量值。
其中，均衡振幅能量值为当前中心采样点数据对应的振幅能量值和均方根能量值的商。
S14：移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值。
通过获取地震数据后，以当前中心采样点数据对应的振幅能量值与当前中心采样点数据对应的均方根能量值的商作为当前中心采样点数据对应的均衡振幅能量值，由于均方根能量值为矩形窗口内的所有检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，而矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，使得均方根能量值的计算不仅和与当前中心采样点数据属于同一地震道数据的采样点数据有关，也与相邻的地震道数据中的采样点数据有关，通过移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值，由于均方根振幅能量是采用相同大小的矩形窗口和相同大小的预设分析时窗计算的，其值是连续变化的，因此在进行地震能量均衡处理之后，可以使相邻地震道之间的能量更均匀，同向轴连续，使地震数据的细节显示的更加清晰。
图2是本发明实施例提供的另一种地震能量均衡方法流程图。如图2所示，该方法包括：
S21：获取地震数据，对地震数据进行预处理。
预处理可以包括噪声压制、动补偿、静补偿中的至少一种。通过对地震数据进行预处理，可以提高地震数据的分辨率，使地震能量均衡处理后的地震数据更加准确。
可选地，地震数据可以是道集数据，也可以是叠加剖面数据，例如偏移叠加剖面数据、水平叠加剖面数据，道集数据可以是共炮点道集数据、共检波点道集数据、共中心点道集数据、共成像点道集数据。以下以共炮点道集数据为例进行说明。
图3是本发明实施例提供的地震数据的局部示意图。图3中示出了3个共炮点道集数据的一部分，其中，横坐标t表示时间，纵坐标i表示空间距离，以其中的一个共炮点道集数据为例，该共炮点道集数据包括由空间上矩阵分布的多个检波点获取的多个地震道数据，例如，图3中显示了3个共炮点道集数据(在纵坐标方向上，每个1～40的区间表示一个共炮点道集数据，共3个共炮点道集数据)，每个共炮点道集数据包括40个地震道数据，每个地震道数据由一个检波点获取。40个检波点的采样起止时间和采样间隔均相同，每个地震道数据均包括对应的检波点沿时间顺序采集的多个采样点数据的集合。
该共炮点道集数据为二维地震数据，二维地震数据常用来确定油气储藏位置，确定油气钻探井位，是一种常见的地震数据形式，在获取二维地震数据时，多个检波点共直线分布。
通常地震数据包括的地震道数据总数I很多，每个地震道数据所包括的采样点数据T也很多，以图3中所示的一个共炮点道集数据为例，该共炮点道集数据包括的地震道数据总数I＝40，每个地震道数据包括的采样点数据总数T＝4000。
需要说明的是，在标识采样点数据时，并未严格对应每个采样点数据进行标识，例如图3，在纵坐标为1的地震道数据中，横坐标4和横坐标5对应的采样点数据表示检波点在相邻的两次采样中获取的数据，而在图像上，对应横坐标4～5之间实际可能包括有多个采样点数据。这是由于检波点实际的采样间隔很小，作图时难以严格按照采样顺序一一标识，因此为了便于说明，仅选取了其中的数个采样点数据进行标识，并假设所选取的采样点数据为相邻的数次采样获得，在实际处理中，应按照采样顺序对采样点数据逐个进行标识。
图4是本发明实施例提供的另一种地震数据的局部示意图。如图4所示，该地震数据为三维地震数据，在获取三维地震数据时，布置的多个检波点呈矩阵分布，且布置成矩形网格状。其中，t坐标表示时间，i坐标和j坐标分别表示在检波点分布的两个正交方向上的空间距离。
以图4中所示的地震数据为例，该地震数据在i坐标方向的地震道数据总数I＝80(图4中仅对i坐标方向的部分地震道数据进行了标识)，在j坐标方向的地震道数据总数J＝40(图4中仅对j坐标方向的部分地震道数据进行了标识)，每个地震道数据包括的采样点数据T＝4000(图4中仅对部分采样点数据进行了标识)。
S22：计算当前中心采样点数据对应的均方根能量值。
其中，均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根。
多个检波点沿矩形窗口的长度方向和宽度方向排列，在矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，每个检波点在预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，当前中心采样点数据为当前位于矩形窗口的中心的检波点在预设分析时窗的中间时刻获取的采样点数据。
在二维地震数据中，矩形窗口的长度方向沿检波点分布的方向(即图3中的纵坐标方向)。在矩形窗口内沿长度方向分布的检波点数量应不少于3个，这样在处理过程中才可以将相邻的检波点进行关联。
如图3所示，矩形I的边a即为矩形窗口在i-t坐标平面内的投影，矩形窗口内沿矩形窗口的长度方向(即图3中矩形Ⅰ的边a的延伸方向)分布有3个检波点，沿矩形窗口的宽度方向(即垂直于i-t坐标平面的方向)分布有1个检波点。在预设分析时窗范围内(即图3中矩形Ⅰ围成的区域内)，沿预设分析时窗方向(即图3中矩形Ⅰ的边b的延伸方向)分布有5个采样点数据，在图3所示的地震数据中，坐标为(4，2)的中心采样点数据的均方根能量值为矩形窗口b内的3个检波点在预设分析时窗内采集的共15个采样点数据对应的振幅能量值的均方根。
在三维地震数据中，矩形窗口内沿矩形窗口的长度方向和宽度方向均至少分布有3个检波点。如图4所示，矩形窗口(如图4中的矩形Ⅱ)的长度沿检波点分布的一个方向(如图4中的i方向)，矩形窗口的宽度沿检波点分布的另一个方向(如图4中的j方向)，预设分析时窗c沿时间方向(如图4中的t方向)。
在图4所示的地震数据中，坐标为(2，2，2)的中心采样点的均方根能量值为矩形窗口内的9个检波点在预设分析时窗内采集的共27个采样点数据对应的振幅能量值的均方根。
在图3和图4中所示的矩形窗口和预设分析时窗仅为示例，实际设置时可以通过增大矩形窗口的长度、宽度或预设分析时窗，改变计算一个中心采样点数据对应的均方根能量值所关联的采样点数据的个数。
在本实施例中，矩形窗口的长度、宽度指的是矩形窗口内沿长度、宽度方向排列的检波点的数量，例如在图3中，由于为二维地震数据，矩形窗口的宽度方向上仅排列有1个检波点，因此矩形窗口的宽度为1，该矩形窗口的长度方向排列有3个检波点，则矩形窗口的长度为3。在图4中，矩形窗口的长度和宽度方向上均排列有3个检波点，即长度和宽度均为3。
可选地，均方根能量值可以设置有下限值，例如1，对于计算出均方根能量值小于下限值的，可以将均方根能量值设定为下限值。在地震数据中，由于部分采样点数据对应的振幅能量值可能为0(如图3中横坐标为1、2或3的采样点数据)，因此可能会使得部分采样点的均方根能量值小于1，甚至等于0，为了能够准确计算均衡振幅能量值，可以对均方根能量值设定下限值，例如将下限值设置为1。
S23：计算当前中心采样点数据对应的均衡振幅能量值。
均衡振幅能量值为当前中心采样点数据对应的振幅能量值和均方根能量值的商。
以图3所示的地震数据为例，当前中心采样点数据(4，2)对应的振幅能量值与当前中心采样点数据(4，2)对应的均方根能量值的商即为当前中心采样点数据(4，2)对应的均衡振幅能量值。
S24：移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值。
在移动矩形窗口或预设分析时窗后，可以返回步骤S22，如此计算出下一中心采样点数据对应的均衡振幅能量值。
通过在空间方向上移动矩形窗口，在时间方向上移动预设分析时窗，采用步骤S22和步骤S23可以计算出地震数据中的所有采样点数据的均衡振幅能量值。矩形窗口和预设分析时窗可以分别移动，每次移动的步长可以为1。
对于给定的地震数据，中心采样点数据的数量与矩形窗口和预设分析时窗的尺寸有关，以图4为例，矩形窗口的尺寸为3×3，在矩形窗口移动的过程中，i坐标为1和i坐标为80的采样点数据、j坐标为1和j坐标为40的采样点数据都无法处于矩形窗口的中心，预设分析时窗的尺寸为3×3，在预设分析时窗移动的过程中，t坐标为1和4000的采样点数据都无法处于预设分析时窗的中心，因此中心采样点数据为i坐标为2～79，j坐标为2～39，t坐标为2～3999的采样点数据，若将矩形窗口的尺寸设置为5×3，预设分析时窗的尺寸保持不变，则中心采样点数据为i坐标为3～78，j坐标为2～39，t坐标为2～3999的采样点数据。
由于矩形窗口的尺寸(包括长度、宽度)相比于检波点的数量，预设分析时窗相比于每个地震道数据包括的采样点数据的数量小得多，因此中心采样点数据在所有的采样点数据中占比较大，通过中心采样点数据足已对地震数据进行其他处理，以进行地震数据的解释，因此可以只对所有中心采样点数据进行地震能量均衡。
S25：计算边缘采样点数据对应的均衡振幅能量值。
以边缘采样点数据对应的振幅能量值与参考能量值的商作为边缘采样点数据对应的均衡振幅能量值，边缘采样点数据为所有的采样点数据中除去所有中心采样点数据之外的采样点数据，参考能量值为与边缘采样点数据的采集时间相隔最近的中心采样点数据中，采集于与边缘采样点数据所对应的检波点相距最近的检波点的中心采样点数据的均方根能量值，或者与边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。以图3为例，在计算i坐标为2，t坐标为1的边缘采样点数据对应的均衡振幅能量值时，参考能量值可以是采集于同一检波点的i坐标为2，t坐标为3的中心采样点数据的均方根能量值；在计算i坐标为1，t坐标为4的边缘采样点数据对应的均衡振幅能量值时，参考能量值可以是i坐标为2，t坐标为4的中心采样点数据的均方根能量值。
通过采用参考能量值计算边缘采样点数据对应的均衡振幅能量值，可以对步骤S22～S24中无法计算的采样点数据也进行处理，这样可以使整个地震数据得到地震能量均衡处理。
步骤S25为可选步骤，地震数据通常包括较多数量的地震道数据，且每个地震道数据中的采样点数据数量也很多，在进行地震能量均衡时，当矩形窗口的尺寸远小于在矩形窗口的长度方向和宽度方向上分布的检波点的数量，且预设分析时窗内的采样点数据的个数n远小于每个地震道数据包括的采样点数据T时，可以忽略边缘采样点，而不执行步骤S25。远小于可以是指较小的量与较大的量的比值小于10％。
需要说明的是，图2所示的方法既适用于垂直地震剖面地震数据的地震能量均衡处理，也适用于地面地震数据的地震能量均衡处理。
图3所示的地震数据为垂直地震剖面地震数据，图5是图3所示的地震数据采用现有的地震能量均衡方法处理后的局部示意图，图6是图3所示的地震数据采用图2所示的地震能量均衡方法处理后的局部示意图，对比图5中的A处和图6中的B处，图6中的波形更加清晰规则，采用图2所示的方法处理后的地震数据具有更高的分辨率。
图7是采用现有的地震能量均衡方法处理后的地面地震数据的局部示意图，图8是采用图2所示的地震能量均衡方法处理后的地面地震数据的局部示意图，对比图7中的C处和图8中的D处，图8中沿D的宽度可以更加明显的区分每条波形，采用图2所示的方法处理后的地面地震数据也具有更高的分辨率。
图9是本发明实施例提供的一种地震能量均衡装置的结构框图，该地震能量均衡装置用于执行图1或图2所示的地震能量均衡方法。如图9所示，该地震能量均衡装置包括获取模块10和处理模块20。
其中，获取模块10用于获取地震数据，地震数据包括由空间上矩阵分布的多个检波点获取的多个地震道数据，多个检波点的采样起止时间和采样间隔均相同，每个地震道数据均包括对应的检波点沿时间顺序采集的多个采样点数据的集合。
处理模块20用于计算当前中心采样点数据对应的均方根能量值，计算当前中心采样点数据对应的均衡振幅能量值，且还用于移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值，均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，多个检波点沿矩形窗口的长度方向和宽度方向排列，在矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，每个检波点在预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，当前中心采样点数据为当前位于矩形窗口的中心的检波点在预设分析时窗的中间时刻获取的采样点数据，均衡振幅能量值为当前中心采样点数据对应的振幅能量值和均方根能量值的商。
可选地，地震数据可以为二维地震数据也可以为三维地震数据，在二维地震数据中，矩形窗口的长度方向沿检波点分布的方向。在矩形窗口内沿长度方向分布的检波点数量应不少于3个。在三维地震数据中，矩形窗口内沿矩形窗口的长度方向和宽度方向均至少分布有3个检波点。
获取模块10在获取地震数据时的方法可以参见前述方法实施例中的步骤S21。
可选地，获取模块10还可以用于对地震数据进行预处理，预处理包括噪声压制、动补偿、静补偿中的至少一种。通过对地震数据进行预处理，可以提高地震数据的分辨率，使地震能量均衡处理后的地震数据更加准确。
可选地，处理模块20还可以用于以边缘采样点数据对应的振幅能量值与参考能量值的商作为边缘采样点数据对应的均衡振幅能量值，边缘采样点数据为所有的采样点数据中除去所有中心采样点数据之外的采样点数据，参考能量值为与边缘采样点数据的采集时间相隔最近的中心采样点数据中，采集于与边缘采样点数据所对应的检波点相距最近的检波点的中心采样点数据的均方根能量值，或者与边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。通过采用参考能量值计算边缘采样点数据对应的均衡振幅能量值，可以对无法计算的采样点数据也进行处理，这样可以使整个地震数据得到地震能量均衡处理。
处理模块20在计算中心采样点数据对应的均衡振幅能量值和边缘采样点对应的均衡振幅能量值时的方法可以参见前述方法实施例中的步骤S22～S25。
通过获取地震数据后，以当前中心采样点数据对应的振幅能量值与当前中心采样点数据对应的均方根能量值的商作为当前中心采样点数据对应的均衡振幅能量值，由于均方根能量值为矩形窗口内的所有检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，而矩形窗口内沿矩形窗口的长度和宽度方向均排列有奇数个检波点，使得均方根能量值的计算不仅和与当前中心采样点数据属于同一地震道数据的采样点数据有关，也与相邻的地震道数据中的采样点数据有关，通过移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出地震数据中的所有中心采样点数据对应的均衡振幅能量值，由于均方根振幅能量是采用相同大小的矩形窗口和相同大小的预设分析时窗计算的，其值是连续变化的，因此在进行地震能量均衡处理之后，可以使相邻地震道之间的能量更均匀，同向轴连续，使地震数据的细节显示的更加清晰。
图10示出了本发明一个示例性实施例提供的计算机设备的结构框图。如图10所示，计算机设备300包括中央处理单元(CPU)301、包括随机存取存储器(RAM)302和只读存储器(ROM)303的系统存储器304，以及连接系统存储器304和中央处理单元301的系统总线305。计算机设备300还包括帮助计算机内的各个器件之间传输信息的基本输入/输出系统(I/O系统)306，和用于存储操作系统313、应用程序314和其他程序模块315的大容量存储设备307。
基本输入/输出系统306包括有用于显示信息的显示器308和用于用户输入信息的诸如鼠标、键盘之类的输入设备309。其中显示器308和输入设备309都通过连接到系统总线305的输入输出控制器310连接到中央处理单元301。基本输入/输出系统306还可以包括输入输出控制器310以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器310还提供输出到显示屏、打印机或其他类型的输出设备。
大容量存储设备307通过连接到系统总线305的大容量存储控制器(未示出)连接到中央处理单元301。大容量存储设备307及其相关联的计算机可读介质为计算机设备300提供非易失性存储。也就是说，大容量存储设备307可以包括诸如硬盘或者CD-ROM驱动器之类的计算机可读介质(未示出)。
不失一般性，计算机可读介质可以包括计算机存储介质和通信介质。计算机存储13介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括RAM、ROM、EPROM、EEPROM、闪存或其他固态存储其技术，CD-ROM、DVD或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知计算机存储介质不局限于上述几种。上述的系统存储器304和大容量存储设备307可以统称为存储器。
根据本发明的各种实施例，计算机设备300还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即计算机设备300可以通过连接在系统总线305上的网络接口单元311连接到网络312，或者说，也可以使用网络接口单元311来连接到其他类型的网络或远程计算机系统(未示出)。
上述存储器还包括一个或者一个以上的程序，一个或者一个以上程序存储于存储器中，被配置由CPU执行。所述一个或者一个以上程序包含用于进行本发明实施例提供的如前所述的地震能量均衡方法的指令。
本发明实施例还提供了一种非临时性计算机可读存储介质，当该存储介质中的指令由地震能量均衡装置的处理器执行时，使得地震能量均衡装置能够执行图1～2中任一项提供的地震能量均衡方法。
一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行图1～2其中任一项提供的地震能量均衡方法的指令。
需要说明的是：上述实施例提供的地震能量均衡装置在进行地震能量均衡时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的地震能量均衡装置与地震能量均衡方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种地震能量均衡方法，其特征在于，所述方法包括：
获取地震数据，所述地震数据包括由多个检波点获取的多个地震道数据，所述多个检波点共直线分布或者呈网格分布，所述多个检波点的采样起止时间和采样间隔均相同，每个所述地震道数据均包括对应的所述检波点沿时间顺序采集的多个采样点数据的集合；
计算当前中心采样点数据对应的均方根能量值，所述均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，所述多个检波点沿所述矩形窗口的长度方向和宽度方向排列，在所述矩形窗口内沿所述矩形窗口的长度和宽度方向均排列有奇数个检波点，每个所述检波点在所述预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，所述当前中心采样点数据为当前位于所述矩形窗口的中心的检波点在所述预设分析时窗的中间时刻获取的采样点数据；
计算所述当前中心采样点数据对应的均衡振幅能量值，所述均衡振幅能量值为所述当前中心采样点数据对应的所述振幅能量值和所述均方根能量值的商；
移动所述矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出所述地震数据中的所有中心采样点数据对应的均衡振幅能量值。
2.根据权利要求1所述的地震能量均衡方法，其特征在于，若所述多个检波点共直线分布，所述矩形窗口内沿所述矩形窗口的长度方向分布有至少3个检波点，沿所述矩形窗口的宽度方向分布有1个检波点。
3.根据权利要求1所述的地震能量均衡方法，其特征在于，若所述多个检波点呈网格分布，所述矩形窗口内沿所述矩形窗口的长度方向和宽度方向均至少分布有3个检波点。
4.根据权利要求1～3任一项所述的地震能量均衡方法，其特征在于，所述方法还包括：
以边缘采样点数据对应的振幅能量值与参考能量值的商作为所述边缘采样点数据对应的均衡振幅能量值，所述边缘采样点数据为所有的所述采样点数据中除去所有中心采样点数据之外的采样点数据，所述参考能量值为与所述边缘采样点数据的采集时间相隔最近的中心采样点数据中，由与所述边缘采样点数据所对应的检波点相距最近的检波点所采集的中心采样点数据的均方根能量值，或者与所述边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。
5.根据权利要求1～3任一项所述的地震能量均衡方法，其特征在于，所述方法还包括：
在所述获取地震数据之后，且在所述计算当前中心采样点数据对应的均方根能量值之前，对所述地震数据进行预处理，所述预处理包括噪声压制、动补偿、静补偿中的至少一种。
6.一种地震能量均衡装置，其特征在于，所述地震能量均衡装置包括：
获取模块，用于获取地震数据，所述地震数据包括由多个检波点获取的多个地震道数据，所述多个检波点共直线分布或者呈网格分布，所述多个检波点的采样起止时间和采样间隔均相同，每个所述地震道数据均包括对应的所述检波点沿时间顺序采集的多个采样点数据的集合；
处理模块，用于计算当前中心采样点数据对应的均方根能量值，计算所述当前中心采样点数据对应的均衡振幅能量值，且还用于移动矩形窗口或预设分析时窗，计算下一中心采样点数据对应的均衡振幅能量值，直至计算出所述地震数据中的所有中心采样点数据对应的均衡振幅能量值，所述均方根能量值为位于矩形窗口内的检波点在预设分析时窗内获取的多个采样点数据对应的振幅能量值的均方根，所述多个检波点沿所述矩形窗口的长度方向和宽度方向排列，在所述矩形窗口内沿所述矩形窗口的长度和宽度方向均排列有奇数个检波点，每个所述检波点在所述预设分析时窗内采集的采样点数据的个数为n，n为奇数且不小于3，所述当前中心采样点数据为当前位于所述矩形窗口的中心的检波点在所述预设分析时窗的中间时刻获取的采样点数据，所述均衡振幅能量值为所述当前中心采样点数据对应的所述振幅能量值和所述均方根能量值的商。
7.根据权利要求6所述的地震能量均衡装置，其特征在于，若所述多个检波点共直线分布，所述矩形窗口内沿所述矩形窗口的长度方向分布有至少3个检波点，沿所述矩形窗口的宽度方向分布有1个检波点。
8.根据权利要求6所述的地震能量均衡装置，其特征在于，若所述多个检波点呈网格分布，所述矩形窗口内沿所述矩形窗口的长度方向和宽度方向均至少分布有3个检波点。
9.根据权利要求7或8所述的地震能量均衡装置，其特征在于，所述处理模块还用于以边缘采样点数据对应的振幅能量值与参考能量值的商作为所述边缘采样点数据对应的均衡振幅能量值，所述边缘采样点数据为所有的所述采样点数据中除去所有中心采样点数据之外的采样点数据，所述参考能量值为与所述边缘采样点数据的采集时间相隔最近的中心采样点数据中，由与所述边缘采样点数据所对应的检波点相距最近的检波点所采集的中心采样点数据的均方根能量值，或者与所述边缘采样点数据采集于同一检波点的中心采样点数据的均方根能量值。
10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现如权利要求1至5任一项所述的地震能量均衡方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810907206.0A CN110824555B (zh) | 2018-08-09 | 2018-08-09 | 地震能量均衡方法及装置、计算机可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810907206.0A CN110824555B (zh) | 2018-08-09 | 2018-08-09 | 地震能量均衡方法及装置、计算机可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110824555A CN110824555A (zh) | 2020-02-21 |
| CN110824555B true CN110824555B (zh) | 2021-07-02 |
Family
ID=69541081
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810907206.0A Active CN110824555B (zh) | 2018-08-09 | 2018-08-09 | 地震能量均衡方法及装置、计算机可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110824555B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114114418B (zh) * | 2020-08-28 | 2024-06-18 | 中国石油化工股份有限公司 | 地震信号增益方法、装置、计算机设备和存储介质 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA1017844A (en) * | 1974-06-03 | 1977-09-20 | Clifford H. Ray | Method of seismic surveying by extracting and displaying seismic properties |
| CN103913767A (zh) * | 2013-08-13 | 2014-07-09 | 中国石油化工股份有限公司 | 多块三维地震数据叠前融合处理能量不均的解决方法 |
| CN103592682B (zh) * | 2013-10-15 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种地震波场能量均衡方法及处理装置 |
| CN105700012B (zh) * | 2016-01-28 | 2017-12-15 | 中国海洋石油总公司 | 一种海上非重复性时移地震分频互均衡处理方法 |
| US10338255B2 (en) * | 2016-04-15 | 2019-07-02 | Pgs Geophysical As | Noise attenuation |
| CN107843919B (zh) * | 2017-10-27 | 2020-02-07 | 中国石油天然气集团有限公司 | 一种微地震监测数据振幅均衡方法 |
-
2018
- 2018-08-09 CN CN201810907206.0A patent/CN110824555B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110824555A (zh) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11073631B2 (en) | Selection of pick values for automated near surface analysis | |
| US10705238B2 (en) | Method and apparatus for processing seismic data | |
| Grigoli et al. | Automated seismic event location by travel‐time stacking: An application to mining induced seismicity | |
| Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
| US10877175B2 (en) | Seismic acquisition geometry full-waveform inversion | |
| US10877174B2 (en) | Diving wave illumination using migration gathers | |
| MX2011003850A (es) | Estimado de señal de dominio de imagen a interferencia. | |
| Khoshnavaz et al. | Velocity-independent estimation of kinematic attributes in vertical transverse isotropy media using local slopes and predictive painting | |
| CN110824555B (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
| Zhang et al. | Automated microseismic event location by amplitude stacking and semblance | |
| CN111123361B (zh) | 垂直地震剖面地震数据规则化重建方法及装置、存储介质 | |
| Hardy | High resolution tomographic MVA with automation | |
| Ge et al. | Reverse travel time imaging of microseismic location | |
| EP2775321A2 (en) | Method and device for calculating time-shifts and time-strains in seismic data | |
| CN114114421A (zh) | 基于深度学习的导向自学习地震数据去噪方法及装置 | |
| CN111352157B (zh) | 一种横波静校正方法及系统 | |
| CN110858000B (zh) | 地震数据重构方法、装置、计算机设备及存储介质 | |
| US20090052279A1 (en) | Method and system for determining azimuth parameters for seismic data processing | |
| CN112649876A (zh) | 建立地震偏移速度模型的方法和装置 | |
| MX2011003852A (es) | Atributos de procesamiento de imagen con inversion de tiempo. | |
| CN114002741B (zh) | 叠前深度偏移方法及装置、计算机可读存储介质 | |
| WO2018093432A1 (en) | Automated mistie analysis and correction across two-dimensional ("2d") seismic surveys | |
| CN113433588B (zh) | 一种分偏移距扫描叠加的近地表速度分析方法 | |
| US20220397690A1 (en) | Detection and Removal of Delayed Seismic Travel Times Produced by Velocity Inversions | |
| CN116520428A (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 |