大尺度范围的场地分类处理方法、装置及电子设备与流程

1.本技术属于地质评估技术领域，具体涉及一种大尺度范围的场地分类处理方法、装置及电子设备。


背景技术：

2.场地条件对工程震害和地震动特性有重要影响。世界各国一般性建、构筑物的抗震设计规范中，均采用场地分类的方法来考虑工程场地对设计地震动参数的影响。我国采用覆盖土层的等效剪切波速v
se
和覆盖土层厚度的双指标来进行场地分类。实际中，一般的可在目标施工场地钻测试孔，基于剪切波波速测试仪及相关设备进行实际测试，进而基于测试数据来得到等效剪切波速v
se
和覆盖土层厚度，从而基于预设国家区域场地分类的分类标准对目标施工场地进行场地分类。
3.而在一些场景中，需要对大尺度范围的场地进行场地分类，以得到对应的场地分类结果数据，如在我国快速发展的地震巨灾保险行业中，需要相关数据进行地震风险量化，以确定保险费率等。类似场景中，基于地域广大、标的目标众多等原因，对大尺度范围的场地分类，显然不宜适用现有的实际钻孔测试来进行场地分类的方法。


技术实现要素：

4.为至少在一定程度上克服相关技术中存在的问题，本技术提供一种大尺度范围的场地分类处理方法、装置及电子设备，以解决如何有效实现大尺度范围的场地分类的技术问题。
5.为实现以上目的，本技术采用如下技术方案：
6.第一方面，
7.本技术提供一种大尺度范围的场地分类处理方法，所述场地分类处理方法包括：
8.获取目标场地的地形地貌特征数据；
9.根据所述地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；
10.在预设的场地分类映射规则约束下，根据所述30m剪切波速值确定所述目标场地的场地分类；
11.其中，所述场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，所述场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。
12.可选地，所述对实测钻孔波速数据的统计分析，包括：
13.获取样本区域的实测钻孔波速剖面数据，根据所述实测钻孔波速剖面数据确定区域内各样本钻孔的覆盖土层厚度数据，并计算得到样本钻孔的覆盖土层等效剪切波速v
se
数据；
14.基于预设国家区域场地分类标准，根据所述覆盖土层厚度数据和所述覆盖土层等效剪切波速v
se
数据，确定区域内各样本钻孔所在场地的分类类别；
15.计算各样本钻孔的30m剪切波速值，并根据所述分类类别对样本钻孔进行分组，基
于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定所述场地分类映射规则。
16.可选地，所述基于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定所述场地分类映射规则，包括：
17.根据各组中样本钻孔的30m剪切波速值，分别对应计算30m剪切波速的均值和标准差，基于正态分布假设，进而对应计算确定满足预设置信度的置信区间，将各组所对应的置信区间作为统计均值建议结果；
18.基于所述统计均值建议结果，对于相邻组的置信区间，针对两置信区间之间的数据重叠部分，采用二分类法进行正确率考察，将正确率最高的样本值确定为分类界限；
19.根据各所述分类界限对统计均值建议结果中置信区间进行界限确定，以得到所述场地分类映射规则。
20.可选地，所述预设置信度为95.4％。
21.可选地，所述地形地貌特征数据包括卫星获取的数字高程dem数据。
22.可选地，所述调用v
s30
预测模型确定目标场地的30m剪切波速值的过程，包括：
23.根据目标场地的数字高程dem数据，采用gmt软件提供的grd—gradient算法工具进行计算处理，得到目标场地的场址地形坡度数据；
24.根据所述场址地形坡度数据，调用wald&allen模型确定目标场地的30m剪切波速值。
25.可选地，所述预设国家区域场地分类类别为中国场地分类类别。
26.第二方面，
27.本技术提供一种大尺度范围的场地分类处理装置，所述场地分类处理装置包括：
28.获取模块，用于获取目标场地的地形地貌特征数据；
29.第一处理模块，用于根据所述地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；
30.第二处理模块，用于在预设的场地分类映射规则约束下，根据所述30m剪切波速值确定所述目标场地的场地分类；
31.其中，所述场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，所述场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。
32.第三方面，
33.本技术提供一种电子设备，包括：
34.存储器，其上存储有可执行程序；
35.处理器，用于执行所述存储器中的所述可执行程序，以实现上述所述方法的步骤。
36.本技术采用以上技术方案，至少具备以下有益效果：
37.本技术的大尺度范围的场地分类处理方法，通过获取目标场地的地形地貌特征数据；根据地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；在预设的场地分类映射规则约束下，根据30m剪切波速值确定目标场地的场地分类；方案中，场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别(如中国场地分类类别)的映射关系，场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。本技术的技术方案中，在对大尺度范围的目标场地进行场地分类的场景中，无需钻孔实测，直接基于获取的目标场地的地形地貌特征数据确定30m剪切波速，进而可根据用于表征30m剪切波速与预设
国家区域场地分类类别的映射关系，确定得到场地分类结果，从而有效实现了大尺度范围的场地分类。
38.本发明的其他优点、目标，和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。
附图说明
39.附图用来提供对本技术的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本技术实施例的附图与本技术的实施例一起用于解释本技术的技术方案，但并不构成对本技术技术方案的限制。
40.图1为本技术一个实施例提供的大尺度范围的场地分类处理方法的流程示意图；
41.图2为本技术一实施例中进行统计分析所采用的钻孔资料剖面深度分布示意图；
42.图3为本技术一实施例中中国现行场地分类标准与美国nehrp分类标准的对比示意图；
43.图4为本技术一实施例中分类界限确定过程中正确率考察的示意说明图；
44.图5为本技术一个实施例提供的大尺度范围的场地分类处理装置的结构示意图；
45.图6为本技术一个实施例提供的电子设备的结构示意图。
具体实施方式
46.为使本技术的目的、技术方案和优点更加清楚，下面将对本技术的技术方案进行详细的描述。显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式，都属于本技术所保护的范围。
47.如背景技术中所述，在一些场景中，需要对大尺度范围的场地进行场地分类，以得到对应的场地分类结果数据，如在我国快速发展的地震巨灾保险行业中，需要相关数据进行地震风险量化，以确定保险费率等。类似场景中，基于地域广大、标的目标众多等原因，对大尺度范围的场地分类，显然不宜适用现有的实际钻孔测试来进行场地分类的方法。
48.针对于此，本技术提出一种大尺度范围的场地分类处理方法，以解决如何有效实现大尺度范围的场地分类的技术问题。
49.在一实施例中，如图1所示，本技术提出的大尺度范围的场地分类处理方法，包括：
50.步骤s110，获取目标场地的地形地貌特征数据，举例而言，这里地形地貌特征数据为卫星获取的地形地貌资料，如卫星获取的数字高程dem数据；
51.之后进行步骤s120，根据地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；
52.这里需要说明的是，在世界范围内，v
s30
(也称为30m剪切波速)被广泛应用于地分类和场地影响模型中，且基于场地地形地貌数据可以利用经验模型快速给出场地v
s30
估计值；
53.举例而言，可根据获得的目标场地的数字高程dem数据，采用gmt软件提供的grd—gradient算法工具进行计算处理，得到目标场地的场址地形坡度数据；进而根据该场址地
形坡度数据，调用wald&allen模型(如wald&allen(2007)模型)确定目标场地的30m剪切波速，在实际实施中，可采用python编程进行该过程的实现。
54.在步骤s120中的得到目标场地的30m剪切波速值后，就可进行步骤s130，在预设的场地分类映射规则约束下，根据30m剪切波速值确定目标场地的场地分类；
55.这里的场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，该场地分类映射规则基于对实测钻孔波速数据的统计分析而确定，举例而言，该实施例中，预设国家区域场地分类类别为中国场地分类类别，即该实施例的技术方案用于大尺度范围的中国场地分类。
56.本技术的技术方案中，在对大尺度范围的目标场地进行场地分类的场景中，无需钻孔实测，直接基于获取的目标场地的地形地貌特征数据确定30m剪切波速，进而可根据用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，确定得到场地分类结果，从而有效实现了大尺度范围的场地分类。
57.为便于理解本技术的技术方案，下面再对场地分类映射规则的确定过程进行一下介绍说明。
58.在一实施例中，预设国家区域场地分类类别为中国场地分类类别，为得到用于表征30m剪切波速与预设国家区域场地分类类别的映射关系的场地分类映射规则，需对实测钻孔波速数据的进行统计分析，该统计分析过程包括：
59.首先，获取样本区域的实测钻孔波速剖面数据，根据实测钻孔波速剖面数据确定区域内各样本钻孔的覆盖土层厚度数据，并计算得到样本钻孔的覆盖土层等效剪切波速v
se
数据；
60.举例而言，该实施例中，利用收集的中国、日本和美国加州地区实测钻孔波速剖面数据(如图2所示给出了采用的中国、日本和美国加州地区共2652个实测钻孔波速数据的剖面深度分布情况)，获取场地覆盖土层厚度d，并计算场地覆盖土层等效剪切波速v
se
(相关计算方式可现于现有公开技术资料，本技术对此不进行详述)。
61.进而可基于预设国家区域场地分类标准(如该实施例中，为下表1所示的中国场地分类标准)，根据覆盖土层厚度数据和覆盖土层等效剪切波速v
se
数据，确定区域内各样本钻孔所在场地的分类类别。
62.表1我国现行建筑抗震规范gb50011-2010(2016版)的场地分类标准
[0063][0064]
之后，计算各样本钻孔的30m剪切波速值(如可基于实测钻孔波速剖面数据采用现有相关方法计算得到各样本钻孔的v
s30
)，并根据分类类别对样本钻孔进行分组，如基于表1进行分类得到i0、i1、ii、iii、iv五个分类类别，则相应的将样本钻孔对应分为五个组；
[0065]
进而基于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定场地分类映射规则。
[0066]
具体的，该实施例中，基于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定场地分类映射规则，包括：
[0067]
(1)、根据各组中样本钻孔的30m剪切波速值，分别对应计算30m剪切波速的均值和标准差，基于正态分布假设，进而对应计算确定满足预设置信度的置信区间，将各组所对应的置信区间作为统计均值建议结果，具体如下表2所示：
[0068]
表2中国规范各类场地对应的v
s30
统计分布结果
[0069][0070]
容易理解的是，该实施例中预设置信度为95.4％，基于正态分布特性，置信区间为表2中所示的(μ-2σ,μ+2σ)，而附图3给出了i0、i1、ii、iii与iv类场地(各组)的v
s30
统计分布以及与美国场地分类标准的对比结果。
[0071]
(2)、根据表2和附图3中各类场地的v
s30
分布结果，考虑到相邻的场地类别之间v
s30
分布范围存在重叠，如i1类((431,1278))和ii类((202,637))场地之间，需要界定相邻的场地类别之间的v
s30
分类界限。
[0072]
作为一种具体的分类界限确定方式，该实施例中，基于统计均值建议结果，对于相邻组的置信区间，针对两置信区间之间的数据重叠部分，采用二分类法进行正确率考察，将正确率最高的样本值确定为分类界限；
[0073]
下面对二分类法在这里的具体应用方式进行一下举例说明：
[0074]
如i0类与i1类，取v
s30
在900至1200m/s之间的不同值，作为i0类与i1类的分类界限进行场地分类，将以上不同分类方式的分类结果与真实场地分类结果(指基于表2，对2652个钻孔中相应钻孔所对应的分类结果)进行比较，考察分类结果相对于真实场地分类的正确率，选取分类正确率最优的结果确定为推荐分类界限；
[0075]
作为典型的二分类问题，这里以i0类为正类，i1类为负类，采用精确率(分类正确的正类个数/判定为正类的样本个数)与召回率(分类正确的正类个数/全部实际的正类个数)作为评价指标，并采用两者的调和平均值作为反映正确率的指标，在900至1200m/s区间内以10m/s为步长计算每次取值的正确率，最终将其峰值对应的最优v
s30
值1140m/s作为i0与i1的分类界限；
[0076]
采用类似方式对各相邻类的分类界限进行确定，如附图4所示为分类界限确定过程中正确率考察的示意说明图，图4中显示了针对相应四个重叠区间的情景，场地分类的正确率随相应选取的v
s30
界限值的变化曲线。
[0077]
(3)、最后在(1)、(2)的基础上，根据各分类界限对统计均值建议结果中置信区间进行界限确定，以得到场地分类映射规则，如在该实施例中，得到的场地分类映射规则如下表3所示：
[0078]
表3推荐的基于v
s30
的中国场地分类方案
[0079][0080]
至此，场地分类映射规则的确定过程完成。
[0081]
本技术的技术方案，是一种基于30m剪切波速v
s30
的预设国家区域场地分类方法，可利用卫星获取的地形地貌数据快速给出基于预设国家区域现行分类标准的宏观场地分类，为强震后震动图绘制、应急救援和灾害评估提供场地分类依据。该方法通过大量实测钻孔波速资料的统计分析，建立现行的预设国家区域(如中国)场地分类与v
s30
的统计关系，减少了模型结果的不确定性，提高了宏观场地分类的准确度。
[0082]
应用本技术的技术方案，能够在缺少实际测量波速资料时，基于卫星获取的场地地形地貌特征快速给出针对预设国家区域现行场地分类标准的宏观场地分类结果。
[0083]
图5为一实施例提供的大尺度范围的场地分类处理装置的结构示意图，如图5所示，该场地分类处理装置300包括：
[0084]
获取模块301，用于获取目标场地的地形地貌特征数据；
[0085]
第一处理模块302，用于根据地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；
[0086]
第二处理模块303，用于在预设的场地分类映射规则约束下，根据30m剪切波速值确定目标场地的场地分类；
[0087]
其中，场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。
[0088]
关于上述相关实施例中的场地分类处理装置300，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0089]
图6为本技术一个实施例提供的电子设备的结构示意图，如图6所示，该电子设备400包括：
[0090]
存储器401，其上存储有可执行程序；
[0091]
处理器402，用于执行存储器401中的可执行程序，以实现上述方法的步骤。
[0092]
关于上述实施例中的电子设备400，其处理器402执行存储器401中的程序的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
[0093]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉该技术的人员在本发明所揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

技术特征：
1.一种大尺度范围的场地分类处理方法，其特征在于，包括：获取目标场地的地形地貌特征数据；根据所述地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；在预设的场地分类映射规则约束下，根据所述30m剪切波速值确定所述目标场地的场地分类；其中，所述场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，所述场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。2.根据权利要求1所述的场地分类处理方法，其特征在于，所述对实测钻孔波速数据的统计分析，包括：获取样本区域的实测钻孔波速剖面数据，根据所述实测钻孔波速剖面数据确定区域内各样本钻孔的覆盖土层厚度数据，并计算得到样本钻孔的覆盖土层等效剪切波速v
se
数据；基于预设国家区域场地分类标准，根据所述覆盖土层厚度数据和所述覆盖土层等效剪切波速v
se
数据，确定区域内各样本钻孔所在场地的分类类别；计算各样本钻孔的30m剪切波速值，并根据所述分类类别对样本钻孔进行分组，基于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定所述场地分类映射规则。3.根据权利要求2所述的场地分类处理方法，其特征在于，所述基于各组所对应的样本钻孔的30m剪切波速的统计分布数据，确定所述场地分类映射规则，包括：根据各组中样本钻孔的30m剪切波速值，分别对应计算30m剪切波速的均值和标准差，基于正态分布假设，进而对应计算确定满足预设置信度的置信区间，将各组所对应的置信区间作为统计均值建议结果；基于所述统计均值建议结果，对于相邻组的置信区间，针对两置信区间之间的数据重叠部分，采用二分类法进行正确率考察，将正确率最高的样本值确定为分类界限；根据各所述分类界限对统计均值建议结果中置信区间进行界限确定，以得到所述场地分类映射规则。4.根据权利要求3所述的场地分类处理方法，其特征在于，所述预设置信度为95.4％。5.根据权利要求1所述的场地分类处理方法，其特征在于，所述地形地貌特征数据包括卫星获取的数字高程dem数据。6.根据权利要求5所述的场地分类处理方法，其特征在于，所述调用v
s30
预测模型确定目标场地的30m剪切波速值的过程，包括：根据目标场地的数字高程dem数据，采用gmt软件提供的grd—gradient算法工具进行计算处理，得到目标场地的场址地形坡度数据；根据所述场址地形坡度数据，调用wald&allen模型确定目标场地的30m剪切波速值。7.根据权利要求1至6中任一项所述的场地分类处理方法，其特征在于，所述预设国家区域场地分类类别为中国场地分类类别。8.一种大尺度范围的场地分类处理装置，其特征在于，所述场地分类处理装置包括：获取模块，用于获取目标场地的地形地貌特征数据；第一处理模块，用于根据所述地形地貌特征数据，调用v
s30
预测模型确定目标场地的30m剪切波速值；第二处理模块，用于在预设的场地分类映射规则约束下，根据所述30m剪切波速值确定
所述目标场地的场地分类；其中，所述场地分类映射规则用于表征30m剪切波速与预设国家区域场地分类类别的映射关系，所述场地分类映射规则基于对实测钻孔波速数据的统计分析而确定。9.一种电子设备，其特征在于，包括：存储器，其上存储有可执行程序；处理器，用于执行所述存储器中的所述可执行程序，以实现权利要求1-7中任一项所述方法的步骤。

技术总结
本申请涉及大尺度范围的场地分类处理方法、装置及电子设备，属于地质评估技术领域，本申请的方法包括：获取目标场地的地形地貌特征数据；根据所述地形地貌特征数据，调用V


技术研发人员：谢俊举 李柯苇 李小军 温增平 安昭
受保护的技术使用者：中国地震局地球物理研究所
技术研发日：2022.09.01
技术公布日：2023/8/24