地震波信息确定方法及装置、计算机可读存储介质

1.本技术涉及地震波技术领域，尤其涉及一种地震波信息确定方法及装置、计算机可读存储介质。


背景技术：

2.在油汽勘探等场景，勘探转头所在的深度几百至几千米的转井中，沿转井的垂直深层方向从上至下方向以预设间距依次排列有多个地震检波器。通过从地面的可控震源产生到达转井口的人工地震波，地震波通过转井口到达转井中设置的各个三分量检波器，由各检波器采集对应的地震波信号的x、y、z的分量信号。通过对检波器采集的三分量地震波信号数据进行分析，对地震波场进行恢复，从而可了解转井所在地下层构造情况。
3.然而，初至地震波信号弱容易受到噪音掩盖，对应深处设置的检波器，根据初至波信号难以准确确定检波器的方位角，导致不能够有效恢复地震波场信息，由此无法准确了解地下层构造。
4.如何准确确定地震检波器方位的方位角，实现地震波信息的有效恢复，是本技术所要解决的技术问题。


技术实现要素：

5.本技术实施例的目的是提供一种车地震波信息确定方法及装置、计算机可读存储介质，用以解决检波器方位角不能准确获取导致地震波场信息恢复困难的问题。
6.为了解决上述技术问题，本说明书是这样实现的：
7.第一方面，提供了一种地震波信息确定方法，包括：
8.确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器；
9.根据各浅层检波器在对应初至时窗内采集的目标地震波的水平分量信号和预设函数，确定对应初至时窗内采集的水平分量信号的偏振方向，以得到各浅层检波器的方位角；
10.根据各深层检波器在预设采集时窗内采集的所述目标地震波的水平分量信号的标量信号同相轴倾角，和基于同相轴倾角的各深层检波器与前向相邻的检波器的水平分量信号相关性，确定各深层检波器的方位角，其中所述预设采集时窗包含所述初至时窗；
11.根据各检波器的方位角和在所述预设采集时窗内对应采集的所述目标地震波的水平分量信号，确定所述目标地震波的径向地震波分量和切向地震波分量。
12.可选的，确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器，包括：
13.根据各检波器在对应初至时窗内采集的所述目标地震波的水平分量信号，确定各检波器的椭圆极化率；
14.根据椭圆极化率，确定各检波器中的浅层检波器和深层检波器。
15.可选的，所述目标地震波的水平分量信号包括横波信号和纵波信号，
16.根据各检波器在对应初至时窗内采集的所述目标地震波的水平分量信号，确定各检波器的椭圆极化率，包括：
17.分别计算目标检波器在对应的目标初至时窗内采集的各横波信号的平均值和各纵波信号的平均值；
18.根据各横波信号、各纵波信号、所述横波信号平均值和所述纵波信号平均值，确定所述目标检波器采集的所述目标地震波的水平分量信号对应的协方差矩阵；
19.根据所述协方差矩阵的最大特征值和最小特征值的比值，确定所述目标检波器的椭圆极化率。
20.可选的，确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器，包括：
21.将沿所述垂直深度方向排列在最上方的检波器确定为浅层检波器；
22.将沿所述垂直深度方向排列在所述最上方的检波器下方的各检波器确定为深层检波器。
23.可选的，根据各浅层检波器在对应初至时窗内采集的目标地震波的水平分量信号和预设函数，确定对应初至时窗内采集的水平分量信号的偏振方向，包括：
24.分别计算目标浅层检波器在对应的初至时窗内采集的各横波信号的平均值和各纵波信号的平均值；
25.根据各横波信号、各纵波信号、所述横波信号平均值和所述纵波信号平均值，确定所述目标浅层检波器采集的所述目标地震波的水平分量信号对应的协方差矩阵；
26.根据所述协方差矩阵的最大特征值对应的特征向量，确定所述目标浅层检波器采集的所述目标地震波的水平分量信号的偏振方向。
27.可选的，确定各深层检波器的方位角，包括：
28.根据目标深层检波器在所述预设采集时窗内目标时刻采集的所述目标地震波的水平分量信号，确定目标标量信号；
29.确定所述目标标量信号的同相轴倾角；
30.分别确定不同方位角下，基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性；
31.根据相关性程度最大时对应的方位角，确定所述目标深层检波器的方位角。
32.可选的，分别确定不同方位角下，基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性，包括：
33.根据所述目标深层检波器在所述预设采集时窗内各时刻采集的所述目标地震波的水平分量信号，确定各时刻对应的标量信号以及不同方位角下所述目标地震波的径向地震波分量和切向地震波分量；
34.确定各时刻的标量信号对应的同相轴倾角；
35.根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性；
36.根据所述径向地震波分量相关性与所述切向地震波分量相关性之和，确定所述水平分量信号相关性。
37.可选的，根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性，包括：
38.根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，确定在所述预设采集时窗内各时刻基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性；
39.确定各时刻中对应水平分量信号相关性最大的目标时刻；
40.根据所述目标时刻确定约束时窗，其中所述约束时窗的时长小于所述预设采集时窗；
41.根据所述目标深层检波器在所述约束时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性。
42.第二方面，提供了一种地震波信息确定装置，包括存储器和与所述存储器电连接的处理器，所述存储器存储有可在所述处理器运行的计算机程序，该计算机程序被该处理器执行时实现如第一方面所述的方法的步骤。
43.第三方面，提供了一种计算机可读存储介质，该计算机可读存储介质上存储计算机程序，该计算机程序被处理器执行时实现如第一方面所述的方法的步骤。
44.在本技术实施例中，通过确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器；根据各浅层检波器在对应初至时窗内采集的目标地震波的水平分量信号和预设函数，确定对应初至时窗内采集的水平分量信号的偏振方向，以得到各浅层检波器的方位角；根据各深层检波器在预设采集时窗内采集的所述目标地震波的水平分量信号的标量信号同相轴倾角，和基于同相轴倾角的各深层检波器与前向相邻的检波器的水平分量信号相关性，确定各深层检波器的方位角；根据各检波器的方位角和在所述预设采集时窗内对应采集的所述目标地震波的水平分量信号，确定所述目标地震波的径向地震波分量和切向地震波分量，由此，可以避免因噪声的掩盖使得初至地震波难以准确地获得检波器方位的问题，提高得到的各个检波器方位角的精度，从而能够准确有效地恢复目标地震波的波场信息。
附图说明
45.此处所说明的附图用来提供对本技术的进一步理解，构成本技术的一部分，本技术的示意性实施例及其说明用于解释本技术，并不构成对本技术的不当限定。在附图中：
46.图1是本技术实施例的地震波信息确定方法的流程示意图。
47.图2是本技术一个示例的地震波信息确定方法的流程示意图。
48.图3是本技术实施例的地震波信息确定装置的结构方框图。
具体实施方式
49.下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。本技术中附图编号仅用于区分方案中的各个步骤，不用于限定各个步骤的执行顺序，具体执行顺序以说明书中描述为准。
50.为了解决现有技术中存在的问题，本技术实施例提供一种地震波信息确定方法，图1是本技术实施例的地震波信息确定方法的流程示意图。
51.如图1所示，包括以下步骤：
52.步骤102，确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器；
53.步骤104，根据各浅层检波器在对应初至时窗内采集的目标地震波的水平分量信号和预设函数，确定对应初至时窗内采集的水平分量信号的偏振方向，以得到各浅层检波器的方位角；
54.步骤106，根据各深层检波器在预设采集时窗内采集的所述目标地震波的水平分量信号的标量信号同相轴倾角，和基于同相轴倾角的各深层检波器与前向相邻的检波器的水平分量信号相关性，确定各深层检波器的方位角，其中所述预设采集时窗包含所述初至时窗；
55.步骤108，根据各检波器的方位角和在所述预设采集时窗内对应采集的所述目标地震波的水平分量信号，确定所述目标地震波的径向地震波分量和切向地震波分量。
56.转井中设置有多个用于采集地震波信号的检波器，这里检波器为三分量检波器，用于采集地震波沿x、y、z三个分量的信号。x和y分量信号统称为地震波的水平分量信号，z称为地震波的垂直分量信号。这些检波器沿设置在转井所在的纵向深度空间中，沿转井的垂直深度方向从上至下按预设间距依次排列。
57.初至时窗是根据目标地震波最开始到达每个检波器的初至时刻确定的，每个检波器具有固定的初至时刻。根据每个检波器的初至时刻定义其对于的一个初至时窗。如果对沿垂直深度方向从上至下排列的多个检波器进行编号，则最上方的编号为第1个检波器，往下依次为第2、3、4......、n个检波器。
58.例如第10个检波器的初至时刻为第100毫秒，设定10毫秒的初至时窗长度，则可以第100毫秒为中心点，得到对应第95毫秒到第105毫秒的初至时间窗。
59.预设采集时窗是各检波器采集地震波信号的时间，例如5秒，在这5秒内，所有检波器均需要从第0时刻开始采集对应的地震波信号，直至5秒的采集时间结束。在地震波未到达对应的检波器时，对应采集时间点的地震波信号为0，直至地震波信号开始到达后。也即，预设采集时窗包含各检波器的初至时窗，不同检波器的初至时窗对应位于该全量的采集时窗的不同时间点。
60.在本技术实施例中，首先将多个检波器划分为浅层和深层部分。顾名思义，浅层检波器为设置在井中上方距离地面井口较近或较浅位置的检波器，深层检波器为设置在井中下方距离地面井口较远或较深位置的检波器。
61.基于上述实施例提供的方案，可选的，在一个实施例中，上述步骤102确定沿垂直
深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器，包括：根据各检波器在对应初至时窗内采集的所述目标地震波的水平分量信号，确定各检波器的椭圆极化率；根据椭圆极化率，确定各检波器中的浅层检波器和深层检波器。
62.在该实施例中，通过计算各检波器的椭圆极化率，来确定对应检波器位于转井的深层还是浅层，椭圆极化率表征了地震波信号的线性偏振强度。
63.椭圆极化率越小，表示地震波信号的线性偏振越强，也即检波器位于浅层，地震波信号较强不容易受到噪音掩盖；反之，椭圆极化率越大，表示地震波信号的线性偏振越弱，也即检波器位于深层，地震波信号较弱容易受到噪音掩盖。
64.如上文所述，地震波包括x和y分量的水平分量信号、和z分量的垂直分量信号，x表示横波，y表示纵波，即水平分量信号包括x横波分量信号和y纵波分量信号。水平分量的标量信号不受检波器方位角的影响，同时也可以体现地震波水平分量偏振波的同相轴倾角和连续性，因此，本技术实施例利用地震波水平分量进行对应检波器的方位角求取。
65.根据各检波器在对应初至时窗内采集的所述目标地震波的水平分量信号，确定各检波器的椭圆极化率，包括：分别计算目标检波器在对应的目标初至时窗内采集的各横波信号的平均值和各纵波信号的平均值；根据各横波信号、各纵波信号、所述横波信号平均值和所述纵波信号平均值，确定所述目标检波器采集的所述目标地震波的水平分量信号对应的协方差矩阵；根据所述协方差矩阵的最大特征值和最小特征值的比值，确定所述目标检波器的椭圆极化率。
66.对于每个检波器，分别计算其对应初至时窗内各采集时间点采集的x分量信号和y分量信号，然后分别计算x横波分量信号的平均值和y纵波分量信号的平均值，如以下公式(1)所示：
[0067][0068]
其中，i表示目标检波器的编号，例如第i个检波器，t表示第i个检波器在其对应初至时窗的采集时间点，xi(t)和yi(t)分别表示第i个检波器在t时刻采集的x横波分量信号和y纵波分量信号，t
fi
表示第i个检波器的初至时刻，w表示初至时窗对应时长的一半，n表示初至时窗内对应的采集时间点数量。
[0069]
在得到各x横波分量信号、各y纵波分量信号分别对应的平均值后，根据各x横波分量信号、各y纵波分量信号、及其平均值，构建协方差矩阵，如下面公式(2)所示：
[0070][0071]
根据构建的协方差矩阵所对应的特征向量，可以计算得到该矩阵的最大特征值和最小特征值，其中一个检波器在其对应的初至时窗内具有一个最大特征值和一个最小特征值。
value decomposition，svd)算法。
[0085]
对于根据椭圆极化率划分浅层和深层检波器的实施例，浅层检波器可能包括多个，则对于每个浅层检波器，分别基于对应深度位置的检波器在对应初至时窗采集的x横波分量和y纵波分量，计算该检波器对应的方位角。
[0086]
对于根据深度位置划分最上方的检波器为浅层检波器的实施例，浅层检波器仅包括一个，则根据该检波器在对应初至时窗采集的x横波分量和y纵波分量，计算该检波器对应的方位角。
[0087]
如上文所述，对于深层检波器，基于预设采集时窗内采集的全量地震波的水平分量信号进行方位角确定。
[0088]
可选的，确定各深层检波器的方位角，包括：根据目标深层检波器在所述预设采集时窗内目标时刻采集的所述目标地震波的水平分量信号，确定目标标量信号；确定所述目标标量信号的同相轴倾角；分别确定不同方位角下，基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性；根据相关性程度最大时对应的方位角，确定所述目标深层检波器的方位角。
[0089]
前向相邻的检波器表示与当前检波器相邻但更靠近进口的检波器，这里相邻的检波器可以包括多个，例如当前检波器为第10个检波器，则在前向相邻的检波器包括3个时，为第7、8、9个检波器；在前向相邻的检波器包括5个时，则为第5、6、7、8、9个检波器。
[0090]
在该实施例中，根据深层检波器在预设采集时窗内采集的地震波的水平分量信号的标量信号、标量信号对应的同相轴倾角、以及目标深层检波器与前向相邻的检波器基于同相轴倾角计算的水平分量信号相关性，来确定目标深层检波器的方位角。深层检波器的各检波器非各自独立计算方位角，需要依赖前向相邻的检波器来计算。
[0091]
首先，利用检波器采集的目标地震波的水平分量信号，即x横波分量信号和y纵波分量信号，计算水平分量信号的标量信号，例如以下公式(3)所示：
[0092][0093]
在一个实施例中，水平分量信号的标量大小也可以用xi(t)和yi(t)的平方和表示，公式(3)表示目标检波器所在的深度域和采集的水平分量信号的时间域的标量数值。
[0094]
地震波的水平分量的标量信号的同相轴是连续的，并且其同相轴的倾角表明了地震波的视速度。因此，本技术实施例中通过计算同相轴倾角来作为确定检波器方位角的约束条件。
[0095]
对于对应某个采集时间点的目标标量信号，可以通过下面的公式(4)计算其同相轴倾角：
[0096][0097]
其中，fft表示快速傅立叶变换(fast fourier transform)函数，h
ht
表示希尔伯特变换函数。
[0098]
通过地震波水平分量信号的标量信号(si(t))深度域方向i(即，对应第i个检波器的深度)和时间域方向t(即，对应第i个检波器采集的水平分量数据的时刻)进行希尔伯特变换，可获得较为精确的标量信号同相轴倾角。
[0099]
上述公式(4)表示一个深度位置的检波器在预设采集时窗内的不同采集时间点采集的水平分量信号，对应计算出一个同相轴倾角。
[0100]
然后，基于目标深层检波器在采集时窗内对应计算的同相轴倾角，采取不同的方位角，分别计算目标深层检波器的水平分量信号与其前向相邻的检波器的水平分量信号的相关性。
[0101]
可选的，分别确定不同方位角下，基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性，包括：根据所述目标深层检波器在所述预设采集时窗内各时刻采集的所述目标地震波的水平分量信号，确定各时刻对应的标量信号以及不同方位角下所述目标地震波的径向地震波分量和切向地震波分量；确定各时刻的标量信号对应的同相轴倾角；根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性；根据所述径向地震波分量相关性与所述切向地震波分量相关性之和，确定所述水平分量信号相关性。
[0102]
假设目标深度检波器i(即，第i个检波器)的方位角为αd，则根据目标深度检波器i在目标采集时间点t采集的水平分量信号，恢复采集的水平分量信号对应的径向地震波分量ri(t)和切向地震波分量ti(t)，例如以下公式(5)所示：
[0103][0104]
目标深度检波器i在目标采集时间点t采集的水平分量信号对应的标量信号同相轴可以根据上述公式(4)计算出，目标深层检波器i与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性可以通过以下公式(6)计算：
[0105][0106]
其中，m表示与目标深层检波器i前向相邻的检波器的数量，例如取前向相邻的3个检波器，或者5个检波器等；τ表示采集时窗，w表示目标采集时间点之前的w个样点，-w表示目标采集时间点之后的w个样点。δx为目标深度检波器i与前向相邻的某检波器之间的的间距大小，t
imax
表示预设采集时窗中目标深度检波器i与前向相邻的检波器的水平分量信号相关性最强的采集时间点，h0表示区分浅层检波器和深层检波器的分界点检波器编号。
[0107]
若以全时窗的预设采集时窗中的所有采集时间点位置进行不同角度的相关性计算，则需要大量的运算量。本技术实施例中，为提高运算的速度，减少大量不必要的运算，选择目标深度检波器i与前向相邻的检波器的水平分量信号相关性最强的采集时间点，再根据该相关性最强的采集时间点，确定对应的时窗，并进行目标深层检波器i与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性计算。
[0108]
可选的，根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波
分量相关性，包括：根据所述目标深层检波器在所述预设采集时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，确定在所述预设采集时窗内各时刻基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的水平分量信号相关性；确定各时刻中对应水平分量信号相关性最大的目标时刻；根据所述目标时刻确定约束时窗，其中所述约束时窗的时长小于所述预设采集时窗；根据所述目标深层检波器在所述约束时窗内各时刻的同相轴倾角、所述目标深层检波器与前向相邻的检波器沿所述垂直深度方向的间距，分别确定基于同相轴倾角的所述目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性。
[0109]
目标深层检波器i与前向相邻的检波器的水平分量信号相关性可以通过以下公式(7)确定：
[0110][0111]
第i个深度检波器的水平分量信号相关性最大值为：
[0112]
c(t
imax
)＝max{ci(t)}
ꢀꢀ
公式(8)。
[0113]
该相关性最大值对应的采集时窗内的时间点，即为目标深度检波器i与前向相邻的检波器的水平分量信号相关性最强的采集时间点t
imax
。
[0114]
以该时间点t
imax
为中心点，可以得到预设时窗的约束时窗，例如前后包括w个采集时间点，该约束时窗远小于全量的采集时窗。通过在约束时窗内进行目标深层检波器与前向相邻的检波器的径向地震波分量相关性、切向地震波分量相关性，可以在保证相关性准确度的前提下显著减少运算量，提高运算效率。
[0115]
返回公式(6)，在取不同的方位角的情况下，根据目标深层检波器i在该时窗内各时刻的同相轴倾角、目标深层检波器i与前向相邻的m个检波器沿垂直深度方向的间距δx，可以分别计算出目标深层检波器i与前向相邻的m个检波器的径向地震波分量相关性和切向地震波分量相关性。
[0116]
不同方位角可以在0-360度的范围内取值，例如按照5度的角度间隔依次选取并进行相关性计算。由此，可以得到对应方位角下，目标深层检波器i与前向相邻的m个检波器的径向地震波分量相关性和切向地震波分量相关性及其相关性之和。
[0117]
如果转井下排列的一列检波器方位角是相同的，则接收的目标地震波信号相似，目标地震波信号按同一个波长往前传播。当径向地震波分量相关性和切向地震波分量相关性的和取得最大值时，此时相关性和最大值对应的方位角αd即为目标深层检波器i的方位角。
[0118]
由此，浅层检波器的方位角和深层检波器的方位角都已确定，则根据公式(5)，各个检波器的方位角分别同时与各自采集的目标地震波的水平分量信号进行旋转计算，可以将采集的原始水平分量信号恢复为原始的目标地震波的径向地震波分量和切向地震波分量。也即，获得目标地震波r分量的sv波和t分量的sh波。
[0119]
现在参考图2，图2是本技术一个示例的地震波信息确定方法的流程示意图。
[0120]
在该实施例中，确定各检波器中的浅层检波器和深层检波器是根据椭圆极化率实
现的。
[0121]
如图2所示，包括以下步骤：
[0122]
步骤202，采集初至时窗各检波器的目标地震波水平分量信号；
[0123]
步骤204，根据各检波器的椭圆极化率确定浅层检波器和深层检波器；
[0124]
步骤206，对于浅层检波器，根据矩阵奇异值分解算法计算检波器方位角；
[0125]
步骤208，对于深层检波器，计算水平分量信号的标量大小；
[0126]
步骤210，计算水平分量信号的标量信号的同相轴倾角；
[0127]
步骤212，根据最强相关性的时间点选择约束时窗；
[0128]
步骤214，扫描检波器方位角，在约束时窗内计算各检波器对应的径向地震波分量相关性与所述切向地震波分量相关性之和；
[0129]
步骤216，根据最大相关性之和输出各检波器方位角；
[0130]
步骤218，根据各检波器方位角，计算目标地震波r分量和地震波t分量。
[0131]
在本技术实施例中，通过确定沿垂直深度方向从上至下按预设间距依次排列的多个检波器中的浅层检波器和深层检波器；根据各浅层检波器在对应初至时窗内采集的目标地震波的水平分量信号和预设函数，确定对应初至时窗内采集的水平分量信号的偏振方向，以得到各浅层检波器的方位角；根据各深层检波器在预设采集时窗内采集的所述目标地震波的水平分量信号的标量信号同相轴倾角，和基于同相轴倾角的各深层检波器与前向相邻的检波器的水平分量信号相关性，确定各深层检波器的方位角；根据各检波器的方位角和在所述预设采集时窗内对应采集的所述目标地震波的水平分量信号，确定所述目标地震波的径向地震波分量和切向地震波分量，由此，可以避免因噪声的掩盖使得初至地震波难以准确地获得检波器方位的问题，提高得到的各个检波器方位角的精度，从而能够准确有效地恢复目标地震波的波场信息。
[0132]
可选的，本技术实施例还提供一种地震波信息确定装置，如图2所示，地震波信息确定装置2000包括存储器2200和与所述存储器2200电连接的处理器2400，所述存储器2200存储有可在所述处理器2400运行的计算机程序，所述计算机程序被所述处理器执行时实现上述任意一种地震波信息确定方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。
[0133]
本技术实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述任意一种地震波信息确定方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(read-only memory，简称rom)、随机存取存储器(random access memory，简称ram)、磁碟或者光盘等。
[0134]
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
[0135]
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下
前者是更佳的实施方式。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本技术各个实施例所述的方法。
[0136]
上面结合附图对本技术的实施例进行了描述，但是本技术并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本技术的启示下，在不脱离本技术宗旨和权利要求所保护的范围情况下，还可做出很多形式，均属于本技术的保护之内。