一种地震震中位置确定方法、装置、设备及可读存储介质

1.本发明涉及地震预警领域，具体而言，涉及一种地震震中位置确定方法、装置、设备及可读存储介质。


背景技术：

2.地震是一种随机的、破坏性极强的自然灾害。地震预警系统的主要是通过p波和s波的速度差来实现的。震中定位是地震预警系统中的关键参数，一方面是影响了震级估算的精度，另一方面也影响了地震预警处置范围的精度。目前震中定位方法中通过经验参数得到的震中距且精度较差，现有波速走时法的假定波速为相同，该方法假定地层的p波速度都是恒定的，因此需要一种考虑不同土层、地壳层的波速差异提出的空间异化波速走时快速震中定位方法，提高震中定位的精度。


技术实现要素：

3.本发明的目的在于提供一种地震震中位置确定方法、装置、设备及可读存储介质，以改善上述问题。为了实现上述目的，本发明采取的技术方案如下：
4.一方面，本技术提供了一种地震震中位置确定方法，所述方法包括：获取第一信息和第二信息，所述第一信息为历史采集到的地勘资料信息，所述第二信息包括至少一个地震台站采集到的地震数据信息和至少一个地震台站的位置信息；基于所述第一信息建立三维模型，得到三维地质模型；将所述第二信息发送至方位角确定模块，得到第三信息，所述第三信息为每个地震台站与所述地震震中的方位角信息；将所述第二信息和所述第三信息发送至训练好的距离预测模型进行处理，得到第四信息，所述第四信息为地震震中的地面位置信息，所述距离预测模型为预测每个所述地震台站与地震震中的地面位置的距离的预测模型；基于所述第一信息、所述第二信息和所述三维地质模型进行地震震中位置三维坐标计算，得到在所述三维地震模型中地震震中位置的三维坐标。
5.可选地，基于所述第一信息建立三维模型，得到三维地质模型，包括：将第一信息中所有的地形等高线图进行格式转化，将转化格式后的地形等高线图进行xoy平面框选处理，得到框选后的地形等高线图；将第一信息中的离散高程数据和所述框选后的地形登高线图发送至midas软件中进行处理，得到地表三维模型；将测量第一信息的测量点作为坐标原点，将第一信息中地层分布图发送至所述地表三维模型进行耦合处理，得到至少一个三维地质单元模型，所述三维地质单元模型为包含了每个测量点的地表上方的地层信息的三维模型；将至少一个所述三维地质单元模型内的区域信息进行离散化处理，得到至少一个三维地质柱状单元，并基于线性插值法对每个所述三维地质柱状单元进行线性插值处理，得到三维地质模型，所述三维地质模型为包含了地表上方和地表下方的地层信息的三维模型。
6.可选地，将所述第二信息发送至方位角确定模块，得到第三信息，包括：调用第二信息内至少一个地震台站的位置信息和每个地震台站采集到的加速度数据；将每个所述地
震台站采集到的加速度数据进行基线校正处理，得到每个校正后的加速度数据；判断每个所述校正后的加速度数据是否大于0，并根据判断结果和每个地震台站的位置确定第一震中方向，所述第一震中方向为所述地震震中位置分别位于每个地震台站的方向；将每个校正后地震台站采集到的加速度数据转化为每个地震台站采集到的加速度数据的向量值，并基于所述第一震中方向和所述向量值进行计算，得到每个地震台站与所述地震震中的方位角。
7.可选地，将所述第二信息和所述第三信息发送至训练好的距离预测模型进行处理，得到第四信息，包括：基于高斯核函数建立svm模型，并使用粒子群优化算法对所述svm模型进行优化，得到优化后的svm模型；将预设的历史地震数据发送至所述优化后的svm模型进行训练，得到训练好的svm模型；将所述第二信息内至少一个地震台站采集到的地震波波形数据发送至训练好的svm模型进行处理，得到每个地震台站与地震震中位置的距离信息，基于第三信息和每个地震台站与地震震中位置的距离信息，确定地震震中的地面位置信息。
8.第二方面，本技术实施例提供了一种地震震中位置确定装置，所述装置包括：第一获取单元，用于获取第一信息和第二信息，所述第一信息为历史采集到的地勘资料信息，所述第二信息包括至少一个地震台站采集到的地震数据信息和至少一个地震台站的位置信息；第一处理单元，用于基于所述第一信息建立三维模型，得到三维地质模型；第一发送单元，用于将所述第二信息发送至方位角确定模块，得到第三信息，所述第三信息为每个地震台站与所述地震震中的方位角信息；第二处理单元，用于将所述第二信息和所述第三信息发送至训练好的距离预测模型进行处理，得到第四信息，所述第四信息为地震震中的地面位置信息，所述距离预测模型为预测每个所述地震台站与地震震中的地面位置的距离的预测模型；第三处理单元，用于基于所述第一信息、所述第二信息和所述三维地质模型进行地震震中位置三维坐标计算，得到在所述三维地震模型中地震震中位置的三维坐标。
9.可选地，所述装置包括：
10.第一处理子单元，用于将第一信息中所有的地形等高线图进行格式转化，将转化格式后的地形等高线图进行xoy平面框选处理，得到框选后的地形等高线图；第二处理子单元，用于将第一信息中的离散高程数据和所述框选后的地形登高线图发送至midas软件中进行处理，得到地表三维模型；第三处理子单元，用于将测量第一信息的测量点作为坐标原点，将第一信息中地层分布图发送至所述地表三维模型进行耦合处理，得到至少一个三维地质单元模型，所述三维地质单元模型为包含了每个测量点的地表上方的地层信息的三维模型；第四处理子单元，用于将至少一个所述三维地质单元模型内的区域信息进行离散化处理，得到至少一个三维地质柱状单元，并基于线性插值法对每个所述三维地质柱状单元进行线性插值处理，得到三维地质模型，所述三维地质模型为包含了地表上方和地表下方的地层信息的三维模型。
11.可选地，所述装置包括：
12.第一调用子单元，用于调用第二信息内至少一个地震台站的位置信息和每个地震台站采集到的加速度数据；第五处理子单元，用于将每个所述地震台站采集到的加速度数据进行基线校正处理，得到每个校正后的加速度数据；第一判断子单元，用于判断每个所述校正后的加速度数据是否大于0，并根据判断结果和每个地震台站的位置确定第一震中方
向，所述第一震中方向为所述地震震中位置分别位于每个地震台站的方向；第一转化子单元，用于将每个校正后地震台站采集到的加速度数据转化为每个地震台站采集到的加速度数据的向量值，并基于所述第一震中方向和所述向量值进行计算，得到每个地震台站与所述地震震中的方位角。
13.可选地，所述装置包括：
14.第一优化子单元，用于基于高斯核函数建立svm模型，并使用粒子群优化算法对所述svm模型进行优化，得到优化后的svm模型；第一训练子单元，用于将预设的历史地震数据发送至所述优化后的svm模型进行训练，得到训练好的svm模型；第六处理子单元，用于将所述第二信息内至少一个地震台站采集到的地震波波形数据发送至训练好的svm模型进行处理，得到每个地震台站与地震震中位置的距离信息，第七处理子单元，用于基于第三信息和每个地震台站与地震震中位置的距离信息，确定地震震中的地面位置信息。
15.第三方面，本技术实施例提供了一种地震震中位置确定设备，所述设备包括存储器和处理器。存储器用于存储计算机程序；处理器用于执行所述计算机程序时实现上述地震震中位置确定方法的步骤。
16.第四方面，本技术实施例提供了一种可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现上述地震震中位置确定方法的步骤。
17.本发明的有益效果为：
18.本发明通过提出一种根据地勘资料将地下土层转化为三维地质模型的方法，并可以根据地勘资料得到地震p波在不同土层的传播速度及入射角、折射角，可选的，可以实现三维地质系统的快速化、信息化、可视化和数值化处理。
19.本发明通过提出一种单台站震中定位方法，包括单台站震中方位角计算，采用向量方法将台站假象为单一质点，则加速度数据可转化为质点运动轨迹，通过对质点运动轨迹进行分析得到震中方位角，并通过转化向量和的方式提出了震中方位角的计算。
20.本发明通过提出一种考虑空间异化波速的震中定位方法，主要包括从地勘资料获得不同土层的拉梅常数与密度，并建立地震p波走时方程，得到震中位置到地震台站的计算公式，通过建立多元非线性方程，求解多元非线性方程，得到震中位置坐标，并进行震中位置坐标转化，得到震中的精确位置。
21.本发明通过上述方法快速的的到震中的精确位置，减少人力物力的同时，在地震来临时增加政府机构的反应时间，减少地震造成的损失，可以更快速的保障人民群众的生命安全。
22.本发明的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
23.为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
24.图1为本发明实施例中所述的一种地震震中位置确定方法流程示意图；
25.图2为本发明实施例中所述的一种地震震中位置确定装置结构示意图；
26.图3是本发明实施例中所述的一种地震震中位置确定设备结构示意图。
具体实施方式
27.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
28.应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
29.实施例1
30.如图1所示，本实施例提供了一种地震震中位置确定方法，其所述方法包括步骤s1、步骤s2、步骤s3、步骤s4和步骤s5。
31.步骤s1、获取第一信息和第二信息，所述第一信息为历史采集到的地勘资料信息，所述第二信息包括至少一个地震台站采集到的地震数据信息和至少一个地震台站的位置信息；
32.步骤s2、基于所述第一信息建立三维模型，得到三维地质模型；
33.步骤s3、将所述第二信息发送至方位角确定模块，得到第三信息，所述第三信息为每个地震台站与所述地震震中的方位角信息；
34.步骤s4、将所述第二信息和所述第三信息发送至训练好的距离预测模型进行处理，得到第四信息，所述第四信息为地震震中的地面位置信息，所述距离预测模型为预测每个所述地震台站与地震震中的地面位置的距离的预测模型；
35.步骤s5、基于所述第一信息、所述第二信息和所述三维地质模型进行地震震中位置三维坐标计算，得到在所述三维地震模型中地震震中位置的三维坐标。
36.可以理解的是本发明通过提出一种根据地勘资料将地下土层转化为三维地质模型的方法，并可以根据地勘资料得到地震p波在不同土层的传播速度及入射角、折射角，可选的，可以实现三维地质系统的快速化、信息化、可视化和数值化处理。
37.可以理解的是本发明通过提出一种单台站震中定位方法，包括单台站震中方位角计算，采用向量方法将台站假象为单一质点，则加速度数据可转化为质点运动轨迹，通过对质点运动轨迹进行分析得到震中方位角，并通过转化向量和的方式提出了震中方位角的计算。
38.可以理解的是本发明通过提出一种考虑空间异化波速的震中定位方法，主要包括从地勘资料获得不同土层的拉梅常数与密度，并建立地震p波走时方程，得到震中位置到地震台站的计算公式，通过建立多元非线性方程，求解多元非线性方程，得到震中位置坐标，
并进行震中位置坐标转化，得到震中的精确位置。
39.在本公开的一种具体实施方式中，所述步骤s2包括步骤s21、步骤s22、步骤s23和步骤s24。
40.步骤s21、将第一信息中所有的地形等高线图进行格式转化，将转化格式后的地形等高线图进行xoy平面框选处理，得到框选后的地形等高线图；
41.步骤s22、将第一信息中的离散高程数据和所述框选后的地形登高线图发送至midas软件中进行处理，得到地表三维模型；
42.步骤s23、将测量第一信息的测量点作为坐标原点，将第一信息中地层分布图发送至所述地表三维模型进行耦合处理，得到至少一个三维地质单元模型，所述三维地质单元模型为包含了每个测量点的地表上方的地层信息的三维模型；
43.步骤s24、将至少一个所述三维地质单元模型内的区域信息进行离散化处理，得到至少一个三维地质柱状单元，并基于线性插值法对每个所述三维地质柱状单元进行线性插值处理，得到三维地质模型，所述三维地质模型为包含了地表上方和地表下方的地层信息的三维模型。
44.可以理解的是本发明通过地勘资料中的不同采集点的数据进行三维建模，得到多个三维建模图形，并通过线性插值法对多个采集点的三维建模图形进行耦合，组合成一个整体的三维地质模型。
45.在本公开的一种具体实施方式中，所述步骤s3包括步骤s31、步骤s32、步骤s33和步骤s34。
46.步骤s31、调用第二信息内至少一个地震台站的位置信息和每个地震台站采集到的加速度数据；
47.步骤s32、将每个所述地震台站采集到的加速度数据进行基线校正处理，得到每个校正后的加速度数据；
48.步骤s33、判断每个所述校正后的加速度数据是否大于0，并根据判断结果和每个地震台站的位置确定第一震中方向，所述第一震中方向为所述地震震中位置分别位于每个地震台站的方向；
49.步骤s34、将每个校正后地震台站采集到的加速度数据转化为每个地震台站采集到的加速度数据的向量值，并基于所述第一震中方向和所述向量值进行计算，得到每个地震台站与所述地震震中的方位角。
50.可以理解的是本发明通过调用多个地震台站的位置和所述地震台站采集到的加速度数据，然后实用校正公式进行校正，然后通过对将加速度数据转化为正东方向、正南方向、正西方向和正北方向的向量值，判断所述加速度数据的方向，进而反向推理的到地震震中方向。
51.可以理解的是校正公式为:
52.n＝n
1-n253.n1为任一采样点采集到的加速度数据，n2为第一个采集到数据的采样点采集到的加速度数据，n为任一采样点校正后的加速度数据。
54.在本公开的一种具体实施方式中，所述步骤s4包括步骤s41、步骤s42、步骤s43和步骤s44。
55.步骤s41、基于高斯核函数建立svm模型，并使用粒子群优化算法对所述svm模型进行优化，得到优化后的svm模型；
56.步骤s42、将预设的历史地震数据发送至所述优化后的svm模型进行训练，得到训练好的svm模型；
57.步骤s43、将所述第二信息内至少一个地震台站采集到的地震波波形数据发送至训练好的svm模型进行处理，得到每个地震台站与地震震中位置的距离信息，
58.步骤s44、基于第三信息和每个地震台站与地震震中位置的距离信息，确定地震震中的地面位置信息。
59.可以理解的是本发明基于例子群优化算法对svm模型进行优化，进而就计算出每个地震台站与地震震中的距离，并通过第三信息和每个地震台站与地震震中的距离推算出地震震中在地面的位置。
60.在本公开的一种具体实施方式中，所述步骤s41包括步骤s411、步骤s412、步骤s413和步骤s414。
61.步骤s411、将所述svm模型内的预设的惩罚参数和预设的高斯核函数参数作为粒子初始位置的两个坐标进行初始化处理，得到至少一个初始化的粒子初始位置；
62.步骤s412、基于预设的适应性函数计算每个粒子的适应度值，并将预设初始速度和粒子初始位置进行更新；
63.步骤s413、判断所述更新次数是否大于预设的最大迭代次数，若所述更新次数小于所述预设的最大迭代次数，则重复上一步骤进行迭代运算，若所述更新次数大于或等于所述预设的最大迭代次数，则将所述每个粒子的适应度值进行对比，将所述的适应度值最大的粒子初始位置作为svm模型内优化后的惩罚参数和优化后的高斯核函数参数；
64.步骤s414、将所述优化后的惩罚参数和所述优化后的高斯核函数参数带入所述svm模型，得到优化后的svm模型。
65.可以理解的是本发明通过对svm模型内的预设的惩罚参数和预设的高斯核函数参数进行迭代处理，得到最优的惩罚参数和高斯核函数参数，进而将所述svm模型进行优化，得到优化后的svm模型。
66.在本公开的一种具体实施方式中，所述步骤s5包括步骤s51、步骤s42、步骤s43和步骤s44。
67.步骤s51、基于所述第一信息和所述三维地质模型，确定所述三维地质模型中至少一个土层的拉美常数和至少一个不同土层的介质密度；
68.步骤s52、基于至少一个所述土层的拉美常数和至少一个所述土层的介质密度计算地震p波在至少一个所述土层的传播速度，并计算出p波在至少一个所述土层间传播时产生的折射角；
69.步骤s53、调用所述第二信息内至少一个地震台站接收到地震p波的时间；
70.步骤s54、基于所述第二信息内至少一个地震台站接收到地震p波的时间和所述每个所述地震p波在每个土层间传播时产生的折射角建立地震震中位置预测方程组；
71.步骤s55、将所述地震震中的第一位置作为初始解，带入所述地震震中位置预测方程组进行迭代计算，得到地震震中位置坐标；
72.步骤s56、将所述地震震中位置坐标转化为经纬度坐标，并发送至所述三维地震模
型内进行标记，得到所述三维地震模型中地震震中的坐标位置。
73.可以理解的是本发明通过计算地震p波在不同土层件的传播折射角，进而精确计算地震p波的传播路线，确定地震震中的深度，并根据所述地震中的深度确定地震震中的坐标位置。
74.可以理解的是首先要根据地震震中到达每个台站的地震p波确定的走时方程，并根据所述走时方程确定地震震中位置预测方程组，所述走时方程为
75.ti＝∑v
ihi
cotai76.其中，ti为所述地震p波到达第i个台站的时间，vi为地震p波在到达第i个台站时通过所有土层的速度，hi为地震p波在到达第i个台站时通过所有土层的高度，ai地震p波在到达第i个台站时通过的所有土层的折射角度。
77.可以理解的是位置预测方程组为
[0078][0079]
其中f1表示地震p波冲第一个台站到第二个台站的走时差方程、f2表示地震p波冲第一个台站到第三个台站的走时差方程、f3表示地震p波冲第一个台站到第四个台站的走时差方程；t1为所述地震p波到达第1个台站的时间、t2为所述地震p波到达第2个台站的时间、t3为所述地震p波到达第3个台站的时间、t4为所述地震p波到达第4个台站的时间；v
1i
为地震p波在到达第1个台站时通过的所有土层的速度、v
2i
为地震p波在到达第2个台站时通过的所有土层的速度、v
3i
为地震p波在到达第3个台站时通过的所有土层的速度、v
4i
为地震p波在到达第4个台站时通过的所有土层的速度；h
1i
为地震p波在到达第1个台站时通过所有土层的高度、h
2i
为地震p波在到达第2个台站时通过所有土层的高度、h
3i
为地震p波在到达第3个台站时通过所有土层的高度、h
4i
为地震p波在到达第4个台站时通过所有土层的高度；α
1i
为地震p波在到达第1个台站时通过的所有土层的折射角度、α
2i
为地震p波在到达第2个台站时通过的所有土层的折射角度、α
3i
为地震p波在到达第3个台站时通过的所有土层的折射角度、α
4i
为地震p波在到达第4个台站时通过的所有土层的折射角度。
[0080]
其中，位置预测方程组中的f(x)＝0，进而构建常微分方程自治系统：
[0081][0082]
可以理解的是求解位置预测方程组的方法在于对f(x)＝0进行求解，其中p(x)为对f(x)进行对x求导并转置，其中x0为震源初始位置坐标(x,y,z)，其中x,y(该位置会对应一个初始位置，按照该位置的土层分层及地震p波的传播路径进行计算)为地震震中的地面位置信息中的震中位置，z为假定深度10km(该10km为大部分地震发生在地下10km)，确定位置后，然后采用迭代方法进行求解，当其x满足f(x)＝0时，则得到方程的解，即为震中位置的精确位置。
[0083]
实施例2
[0084]
如图2所示，本实施例提供了一种地震震中位置确定装置，所述装置包括第一获取
单元701、第一处理单元702、第一发送单元703、第二处理单元704和第三处理单元705。
[0085]
第一获取单元701，用于获取第一信息和第二信息，所述第一信息为历史采集到的地勘资料信息，所述第二信息包括至少一个地震台站采集到的地震数据信息和至少一个地震台站的位置信息；
[0086]
第一处理单元702，用于基于所述第一信息建立三维模型，得到三维地质模型；
[0087]
第一发送单元703，用于将所述第二信息发送至方位角确定模块，得到第三信息，所述第三信息为每个地震台站与所述地震震中的方位角信息；
[0088]
第二处理单元704，用于将所述第二信息和所述第三信息发送至训练好的距离预测模型进行处理，得到第四信息，所述第四信息为地震震中的地面位置信息，所述距离预测模型为预测每个所述地震台站与地震震中的地面位置的距离的预测模型；
[0089]
第三处理单元705，用于基于所述第一信息、所述第二信息和所述三维地质模型进行地震震中位置三维坐标计算，得到在所述三维地震模型中地震震中位置的三维坐标。
[0090]
在本公开的一种具体实施方式中，所述第一处理单元702包括第一处理子单元7021、第二处理子单元7022、第三处理子单元7023和第四处理子单元7024。
[0091]
第一处理子单元7021，用于将第一信息中所有的地形等高线图进行格式转化，将转化格式后的地形等高线图进行xoy平面框选处理，得到框选后的地形等高线图；
[0092]
第二处理子单元7022，用于将第一信息中的离散高程数据和所述框选后的地形登高线图发送至midas软件中进行处理，得到地表三维模型；
[0093]
第三处理子单元7023，用于将测量第一信息的测量点作为坐标原点，将第一信息中地层分布图发送至所述地表三维模型进行耦合处理，得到至少一个三维地质单元模型，所述三维地质单元模型为包含了每个测量点的地表上方的地层信息的三维模型；
[0094]
第四处理子单元7024，用于将至少一个所述三维地质单元模型内的区域信息进行离散化处理，得到至少一个三维地质柱状单元，并基于线性插值法对每个所述三维地质柱状单元进行线性插值处理，得到三维地质模型，所述三维地质模型为包含了地表上方和地表下方的地层信息的三维模型。
[0095]
在本公开的一种具体实施方式中，所述第一发送单元703包括第一调用子单元7031、第五处理子单元7032、第一判断子单元7033和第一转化子单元7034。
[0096]
第一调用子单元7031，用于调用第二信息内至少一个地震台站的位置信息和每个地震台站采集到的加速度数据；
[0097]
第五处理子单元7032，用于将每个所述地震台站采集到的加速度数据进行基线校正处理，得到每个校正后的加速度数据；
[0098]
第一判断子单元7033，用于判断每个所述校正后的加速度数据是否大于0，并根据判断结果和每个地震台站的位置确定第一震中方向，所述第一震中方向为所述地震震中位置分别位于每个地震台站的方向；
[0099]
第一转化子单元7034，用于将每个校正后地震台站采集到的加速度数据转化为每个地震台站采集到的加速度数据的向量值，并基于所述第一震中方向和所述向量值进行计算，得到每个地震台站与所述地震震中的方位角。
[0100]
在本公开的一种具体实施方式中，所述第二处理单元704包括第一优化子单元7041、第一训练子单元7042、第六处理子单元7043和第七处理子单元7044。
[0101]
第一优化子单元7041，用于基于高斯核函数建立svm模型，并使用粒子群优化算法对所述svm模型进行优化，得到优化后的svm模型；
[0102]
第一训练子单元7042，用于将预设的历史地震数据发送至所述优化后的svm模型进行训练，得到训练好的svm模型；
[0103]
第六处理子单元7043，用于将所述第二信息内至少一个地震台站采集到的地震波波形数据发送至训练好的svm模型进行处理，得到每个地震台站与地震震中位置的距离信息，
[0104]
第七处理子单元7044，用于基于第三信息和每个地震台站与地震震中位置的距离信息，确定地震震中的地面位置信息。
[0105]
在本公开的一种具体实施方式中，所述第一优化子单元7041包括第八处理子单元70411、第一计算子单元70412、第二判断子单元70413和第二优化子单元70414。
[0106]
第八处理子单元70411，用于将所述svm模型内的预设的惩罚参数和预设的高斯核函数参数作为粒子初始位置的两个坐标进行初始化处理，得到至少一个初始化的粒子初始位置；
[0107]
第一计算子单元70412，用于基于预设的适应性函数计算每个粒子的适应度值，并将预设初始速度和粒子初始位置进行更新；
[0108]
第二判断子单元70413，用于判断所述更新次数是否大于预设的最大迭代次数，若所述更新次数小于所述预设的最大迭代次数，则重复上一步骤进行迭代运算，若所述更新次数大于或等于所述预设的最大迭代次数，则将所述每个粒子的适应度值进行对比，将所述的适应度值最大的粒子初始位置作为svm模型内优化后的惩罚参数和优化后的高斯核函数参数；
[0109]
第二优化子单元70414，用于将所述优化后的惩罚参数和所述优化后的高斯核函数参数带入所述svm模型，得到优化后的svm模型。
[0110]
在本公开的一种具体实施方式中，所述第三处理单元705包括第九处理子单元7051、第二计算子单元7052、第二调用子单元7053、第十处理子单元7054、第三计算子单元7055和第二转化子单元7056
[0111]
第九处理子单元7051，用于基于所述第一信息和所述三维地质模型，确定所述三维地质模型中至少一个土层的拉美常数和至少一个不同土层的介质密度；
[0112]
第二计算子单元7052，用于基于至少一个所述土层的拉美常数和至少一个所述土层的介质密度计算地震p波在至少一个所述土层的传播速度，并计算出p波在至少一个所述土层间传播时产生的折射角；
[0113]
第二调用子单元7053，用于调用所述第二信息内至少一个地震台站接收到地震p波的时间；
[0114]
第十处理子单元7054，用于基于所述第二信息内至少一个地震台站接收到地震p波的时间和所述每个所述地震p波在每个土层间传播时产生的折射角建立地震震中位置预测方程组；
[0115]
第三计算子单元7055，用于将所述地震震中的第一位置作为初始解，带入所述地震震中位置预测方程组进行迭代计算，得到地震震中位置坐标；
[0116]
第二转化子单元7056，用于将所述地震震中位置坐标转化为经纬度坐标，并发送
至所述三维地震模型内进行标记，得到所述三维地震模型中地震震中的坐标位置。
[0117]
需要说明的是，关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0118]
实施例3
[0119]
相应于上面的方法实施例，本公开实施例还提供了一种地震震中位置确定设备，下文描述的一种地震震中位置确定设备与上文描述的一种地震震中位置确定方法可相互对应参照。
[0120]
图3是根据一示例性实施例示出的一种地震震中位置确定设备800的框图。如图3所示，该地震震中位置确定设备800可以包括：处理器801，存储器802。该地震震中位置确定设备800还可以包括多媒体组件803，输入/输出(i/o)接口804，以及通信组件805中的一者或多者。
[0121]
其中，处理器801用于控制该地震震中位置确定设备800的整体操作，以完成上述的地震震中位置确定方法中的全部或部分步骤。存储器802用于存储各种类型的数据以支持在该地震震中位置确定设备800的操作，这些数据例如可以包括用于在该地震震中位置确定设备800上操作的任何应用程序或方法的指令，以及应用程序相关的数据，例如联系人数据、收发的消息、图片、音频、视频等等。该存储器802可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。多媒体组件803可以包括屏幕和音频组件。其中屏幕例如可以是触摸屏，音频组件用于输出和/或输入音频信号。例如，音频组件可以包括一个麦克风，麦克风用于接收外部音频信号。所接收的音频信号可以被进一步存储在存储器802或通过通信组件805发送。音频组件还包括至少一个扬声器，用于输出音频信号。i/o接口804为处理器801和其他接口模块之间提供接口，上述其他接口模块可以是键盘，鼠标，按钮等。这些按钮可以是虚拟按钮或者实体按钮。通信组件805用于该地震震中位置确定设备800与其他设备之间进行有线或无线通信。无线通信，例如wi-fi，蓝牙，近场通信(near fieldcommunication，简称nfc)，2g、3g或4g，或它们中的一种或几种的组合，因此相应的该通信组件805可以包括：wi-fi模块，蓝牙模块，nfc模块。
[0122]
在一示例性实施例中，地震震中位置确定设备800可以被一个或多个应用专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digitalsignal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述的一种地震震中位置确定方法。
[0123]
在另一示例性实施例中，还提供了一种包括程序指令的计算机可读存储介质，该程序指令被处理器执行时实现上述的地震震中位置确定方法的步骤。例如，该计算机可读存储介质可以为上述包括程序指令的存储器802，上述程序指令可由地震震中位置确定设
备800的处理器801执行以完成上述的地震震中位置确定方法。
[0124]
实施例4
[0125]
相应于上面的方法实施例，本公开实施例还提供了一种可读存储介质，下文描述的一种可读存储介质与上文描述的一种地震震中位置确定方法可相互对应参照。
[0126]
一种可读存储介质，可读存储介质上存储有计算机程序，计算机程序被处理器执行时实现上述方法实施例的地震震中位置确定方法的步骤。
[0127]
该可读存储介质具体可以为u盘、移动硬盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、磁碟或者光盘等各种可存储程序代码的可读存储介质。
[0128]
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
[0129]
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。