一种基于7级以上地震的救灾物资储备库选址方法及系统

技术领域

本发明涉及救灾物资储备库选址技术领域，特别是涉及一种基于7级以上地震的救灾物资储备库选址方法及系统。

背景技术

在7级以上地震发生后，通往受灾区域的道路会不同程度的损毁，甚至市内的救灾物资储备库也在短时间内失效，这就导致了救灾物资无法通过常规方式送达受灾区域。在7级以上地震发生后，位于合适位置的救灾物资储备库显然能及时、足量为受灾区域的灾民提供有效的民生物资保障。因此，大型救灾物资储备库的选址对应对7级以上地震是十分关键的。

目前，救灾物资储备库的选址方法主要是：以最大区域覆盖法、最小时间法、P-中值等方法为基础，将救灾物资储备库建设成本、运输成本、采购成本等因素作为约束条件，利用遗传算法、改进型遗传算法进行选址方法模型的构建。但是上述救灾物资储备库的选址方法均应用于道路不被破坏、高度重视救灾成本的情景下，并不符合实际需求。

发明内容

本发明的目的是提供一种基于7级以上地震的救灾物资储备库选址方法及系统，以达到现实可行性的目的。

为实现上述目的，本发明提供了如下方案：

一种基于7级以上地震的救灾物资储备库选址方法，包括：

根据研究区域的地震文献资料，确定所述研究区域内地震震级为7级以上的地震风险区域；

根据所述地震风险区域和地震烈度衰减模型，确定所述研究区域内地震烈度为8级以上的烈度区域；

采用遥感识别建筑物技术识别所述烈度区域内的建筑物以生成建筑物分布图，然后采用八连通方式对所述建筑物分布图进行处理，以得到受灾区域；

确定直升机运送救灾物资时的最远安全飞行距离以及确定所述研究区域内救灾物资储备库的选址候选区域；

根据约束条件对所述选址候选区域进行筛选，以确定最佳选址区域；所述约束条件至少包括所述最远安全飞行距离和所述受灾区域之间的约束关系、以及所述研究区域的自然灾害、地形指数、交通条件和城市规划。

可选地，根据研究区域的地震文献资料，确定所述研究区域内地震震级为7级以上的地震风险区域，具体包括：

根据研究区域的地震科学研究的预判成果和统计规律，确定所述研究区域内的地震风险候选区域；

确定每个所述地震风险候选区域对应的最大地震震级；

将所述最大地震震级与设定阈值进行比较，并将大于或者等于所述设定阈值的所述最大地震震级对应的所述地震风险候选区域确定为地震震级为7级以上的地震风险区域；其中，所述设定阈值为7。

可选地，所述地震风险区域为椭圆形区域；根据所述地震风险区域和地震烈度衰减模型，确定所述研究区域内地震烈度为8级以上的烈度区域，具体包括：

根据所述地震风险区域和地震烈度衰减模型，确定地震烈度为8级时的椭圆等震线；

在所述地震风险区域上增加所述椭圆等震线后再进行椭圆形圈定，得到地震烈度为8级以上的烈度区域。

可选地，所述地震烈度衰减模型的表达式为：

I＝6.8053+1.2972Ms-4.7603log(R

I＝5.3315+1.2013Ms-4.1917log(R

其中，I为地震烈度，I＝8；Ms为地震震级；R

可选地，所述确定直升机运送救灾物资时的最远安全飞行距离，具体包括：

根据公式D

其中，直升机安全飞行的约束条件为：

v＝MCS；

q＝MP；

式中，D

可选地，确定所述研究区域内救灾物资储备库的选址候选区域，具体包括：

将所述研究区域内的第一区域和/或第二区域确定为救灾物资储备库的选址候选区域；所述第一区域为以高速公路入口为圆心，以5km为半径所确定的区域；所述第二区域为以火车货站为圆心，以5km为半径所确定的区域。

可选地，所述自然灾害至少包括洪灾；所述约束条件还包括人均GDP统计结果；根据约束条件对所述选址候选区域进行筛选，以确定最佳选址区域，具体包括：

根据所述最远安全飞行距离和所述受灾区域计算出多个第一目标候选区，并将符合所述选址候选区域的第一目标候选区域确定为第二目标候选区域；

根据所述研究区域内的洪灾区域和8级以上的烈度区域对所述第二目标候选区域进行筛选，并将符合所述研究区域内的非洪灾区域和非8级以上的烈度区域的第二目标候选区域确定为第三目标候选区域；

根据两年内所述研究区域的人均GDP统计结果对所述第三目标候选区域进行筛选，并将符合两年内所述研究区域的人均GDP统计结果的第三目标候选区域确定为第四目标候选区域；

根据所述研究区域的地形指数对所述第四目标候选区域进行筛选，并将符合所述研究区域的地形指数约束的第四目标候选区域确定为第五目标候选区域；

根据所述研究区域的交通条件对所述第五目标候选区域进行筛选，并将符合所述研究区域的交通条件约束的第五目标候选区域确定为第六目标候选区域；

根据所述研究区域的城市规划对所述第六目标候选区域进行筛选，并将符合所述研究区域的城市规划约束的第六目标候选区域确定为最佳选址区域。

一种基于7级以上地震的救灾物资储备库选址系统，包括：

地震风险区域确定模块，用于根据研究区域的地震文献资料，确定所述研究区域内地震震级为7级以上的地震风险区域；

烈度区域确定模块，用于根据所述地震风险区域和地震烈度衰减模型，确定所述研究区域内地震烈度为8级以上的烈度区域；

受灾区域生成模块，用于采用遥感识别建筑物技术识别所述烈度区域内的建筑物以生成建筑物分布图，然后采用八连通方式对所述建筑物分布图进行处理，以得到受灾区域；

最远安全飞行距离和选址候选区域确定模块，用于确定直升机运送救灾物资时的最远安全飞行距离以及确定所述研究区域内救灾物资储备库的选址候选区域；

最佳选址区域确定模块，用于根据约束条件对所述选址候选区域进行筛选，以确定最佳选址区域；所述约束条件至少包括所述最远安全飞行距离和所述受灾区域之间的约束关系、以及所述研究区域的自然灾害、地形指数、交通条件和城市规划。

可选地，所述地震风险区域确定模块，具体包括：

地震风险候选区域单元，用于根据研究区域的地震科学研究的预判成果和统计规律，确定所述研究区域内的地震风险候选区域；

最大地震震级确定单元，用于确定每个所述地震风险候选区域对应的最大地震震级；

地震风险区域确定单元，用于将所述最大地震震级与设定阈值进行比较，并将大于或者等于所述设定阈值的所述最大地震震级对应的所述地震风险候选区域确定为地震震级为7级以上的地震风险区域；其中，所述设定阈值为7。

可选地，所述自然灾害至少包括洪灾；所述约束条件还包括人均GDP统计结果；所述最佳选址区域确定模块，具体包括：

第二目标候选区域确定单元，用于根据所述最远安全飞行距离和所述受灾区域计算出多个第一目标候选区，并将符合所述选址候选区域的第一目标候选区域确定为第二目标候选区域；

第三目标候选区域确定单元，用于根据所述研究区域内的洪灾区域和8级以上的烈度区域对所述第二目标候选区域进行筛选，并将符合所述研究区域内的非洪灾区域和非8级以上的烈度区域的第二目标候选区域确定为第三目标候选区域；

第四目标候选区域确定单元，用于根据两年内所述研究区域的人均GDP统计结果对所述第三目标候选区域进行筛选，并将符合两年内所述研究区域的人均GDP统计结果的第三目标候选区域确定为第四目标候选区域；

第五目标候选区域确定单元，用于根据所述研究区域的地形指数对所述第四目标候选区域进行筛选，并将符合所述研究区域的地形指数约束的第四目标候选区域确定为第五目标候选区域；

第六目标候选区域确定单元，用于根据所述研究区域的交通条件对所述第五目标候选区域进行筛选，并将符合所述研究区域的交通条件约束的第五目标候选区域确定为第六目标候选区域；

最佳选址区域确定单元，用于根据所述研究区域的城市规划对所述第六目标候选区域进行筛选，并将符合所述研究区域的城市规划约束的第六目标候选区域确定为最佳选址区域。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

一是本发明通过研究区域的地震文献资料确定出地震烈度为8级以上的烈度区域，更有利于快速实现受灾区域选取。

二是本发明通过遥感识别建筑物技术对受灾区域内的建筑物识别，并通过八联通方式设置标签，有效保证了受灾区域选取的精准度。

三是通常情况下，是对救灾道路运输进行约束，而本发明是根据直升机机型对最远安全飞行距离进行求解，有效的提升了救灾物资运输的安全性和效率。

四是本发明通过对选址候选区的地形指数、灾害情况、市政规划等进行约束，有效保证了选址的现实可行性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种基于7级以上地震的救灾物资储备库选址方法的流程图；

图2为本发明实施例提供的8级以上的烈度区域选取结果图；

图3为本发明实施例提供的受灾区域标记结果图；

图4为本发明实施例提供的候选区排除自然灾害区域后的结果图；

图5为本发明实施例提供的根据城市规划指标选取参考区域后的结果图；

图6为本发明实施例提供的一种基于7级以上地震的救灾物资储备库选址系统的结构图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

目前，所采用的救灾物资储备库的选址方法，与地震科学研究的预判成果和统计规律结合较少，且也没有考虑地震后道路无法使用、临近救灾物资储备库失效等问题，鉴于此，本发明提供一种基于7级以上地震的救灾物资储备库选址方法及系统，首先通过研究区域的地震文献资料获得7级以上地震8级以上烈度的受灾区域；其次通过直升机设备确定安全最远飞行距离；然后将研究区内高速公路入口和火车货站为圆心，以5km为半径的区域作为选址候选区；最后根据选址候选区、受灾区域，并结合安全最远飞行距离、灾害、直升机起降、地形指数等条件进行选址，确保大型救灾物资储备库的选址能够符合实际需求，且区域内发生7级以上地震后，也能够使受灾区域在短时间获得民生物资的救助。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

实施例二

本实施例提供了一种基于7级以上地震的救灾物资储备库选址方法，如图1所示，包括如下步骤：

步骤100：根据研究区域的地震文献资料，确定所述研究区域内地震震级为7级以上的地震风险区域。

此步骤具体为：

首先根据研究区域的地震科学研究的预判成果和统计规律，确定所述研究区域内的地震风险候选区域。

其次确定每个所述地震风险候选区域对应的最大地震震级；

将所述最大地震震级与设定阈值进行比较，并将大于或者等于所述设定阈值的所述最大地震震级对应的所述地震风险候选区域确定为地震震级为7级以上的地震风险区域。

其中，所述设定阈值为7；所述地震风险区域为椭圆形区域。

步骤200：根据所述地震风险区域和地震烈度衰减模型，确定所述研究区域内地震烈度为8级以上的烈度区域。

此步骤具体为：

步骤A：根据所述地震风险区域和地震烈度衰减模型，确定地震烈度为8级时的椭圆等震线。

一个示例为：根据所述地震风险区域的边界位置，确定地震烈度为8级时的扩展边界。

根据地震烈度为8级时的扩展边界和地震烈度衰减模型，确定地震烈度为8级时的椭圆等震线。

所述地震烈度衰减模型的表达式为：

I＝6.8053+1.2972Ms-4.7603log(R

I＝5.3315+1.2013Ms-4.1917log(R

其中，I为地震烈度，I＝8；Ms为地震震级；R

步骤B：在所述地震风险区域上增加所述椭圆等震线后再进行椭圆形圈定，得到地震烈度为8级以上的烈度区域。

一个示例为：在所述地震风险区域上下左右增加为地震烈度为8级的椭圆等震线的长半轴、短半轴的长度，然后进行椭圆形圈定，得到如图2所示的地震烈度为8级以上的烈度区域。该地震烈度为8级以上的烈度区域为椭圆形区域。

步骤300：采用遥感识别建筑物技术识别所述烈度区域内的建筑物以生成建筑物分布图，然后采用八连通方式对所述建筑物分布图进行处理，以得到受灾区域。其中，所述受灾区域也为椭圆形区域。

一个示例为：选择IGBP全球植被分类体系，以城市和建筑区(分类代码为13)作为建筑区域提取的主要目标，填充像元，并在该区域周边设计缓冲区，接着通过目视译对城市周边的孤立像元进行判别，将村庄、工厂等地物类型纳入承灾体范围，最终生成建筑物分布图。

为消除孤立建筑的影响，对建筑物分布图进行腐蚀膨胀运算，并以八连通方法作为判断依据，将联通区域设置为承灾体，即得到如图3所示的受灾区域。

步骤400：确定直升机运送救灾物资时的最远安全飞行距离以及确定所述研究区域内救灾物资储备库的选址候选区域。

一个示例为：给定救灾物资运输直升机机型，考虑到直升机安全飞行要求，求出该运送救灾物资时的最远安全飞行距离。

所述最远安全飞行距离的计算公式如下所示；

D

直升机安全飞行的约束条件为：v＝MCS；

q＝MP；

式中，D

一个示例为：将所述研究区域内的第一区域和/或第二区域确定为救灾物资储备库的选址候选区域；所述第一区域为以高速公路入口为圆心，以5km为半径所确定的区域；所述第二区域为以火车货站为圆心，以5km为半径所确定的区域。

步骤500：根据约束条件对所述选址候选区域进行筛选，以确定最佳选址区域；所述约束条件至少包括所述最远安全飞行距离和所述受灾区域之间的约束关系、以及所述研究区域的自然灾害、地形指数、交通条件和城市规划。

其中，所述自然灾害至少包括洪灾；所述约束条件还包括人均GDP统计结果；此步骤具体包括：

步骤a：根据所述最远安全飞行距离和所述受灾区域计算出多个第一目标候选区，检验求出的第一目标候选区是否符合要求，即是否在选址候选区域内；如不符合要求，则放弃该第一目标候选区；如果符合要求，则保存该第一目标候选区，并将符合所述选址候选区域的第一目标候选区域确定为第二目标候选区域。

进一步地，根据所述最远安全飞行距离和所述受灾区域(该受灾区域也为椭圆形区域)计算第一目标候选区的公式如下所示：

a＝6378137m；

b＝6356752m；

y＝K(L-L

其中，L为经度，B为纬度，L

约束条件：

式中，n表示受灾区域的个数；(x

步骤b：在考虑研究区域内的洪灾区域和8级以上的烈度区域的情况下，统计出第二目标候选区域是否符合所述研究区域内的非洪灾区域和非8级以上的烈度区域；如不符合要求，则放弃该第二目标候选区域；如果符合要求，则保存该第二目标候选区域，并将符合所述研究区域内的非洪灾区域和非8级以上的烈度区域的第二目标候选区域确定为第三目标候选区域，具体如图4所示。

进一步地，步骤b中，所述洪灾区域为洪灾统计≥0.3次/a的区域。

步骤c：在考虑两年内研究区域的人均GDP统计结果的情况下，统计出第三目标候选区域是否符合要求；如不符合要求，则放弃该第三目标候选区域；如果符合要求，则保存该第三目标候选区域，并将符合两年内所述研究区域的人均GDP统计结果的第三目标候选区域确定为第四目标候选区域。

进一步地，步骤c中，此处的人均GDP统计结果为在市内排名非最后1名。

步骤d：在考虑研究区域的地形指数的情况下，检验求出第四目标候选区域是否符合要求；如不符合要求，则放弃该第四目标候选区域；如果符合要求，则保存该第四目标候选区域，并将符合所述研究区域的地形指数约束条件的第四目标候选区域确定为第五目标候选区域。

进一步地，步骤d中，所述地形指数求公式如下所示；

高程均值计算公式：

高程方差计算公式：

高程标准差计算公式：

其中，Xi为每个像元的高程，N为区域像元数量。

所述地形指数的约束方法实现过程为：

以5×5窗口为单元，候选区所含单元坡度需具有小于25％的单元，且坡度越小越优先；候选区所含单元高程均值位于研究区内所有单元的前50％；候选区所含单元高程标准差位于研究区内所有单元的后15％。

步骤e：根据所述研究区域的交通条件对所述第五目标候选区域进行筛选，并将符合所述研究区域的交通条件约束的第五目标候选区域确定为第六目标候选区域。

例如：在考虑交通便利因素的情况下，求出第五目标候选区域到达各选定区域坐标的道路距离，并根据各道路距离之和进行从小到大排序，选择前m个各道路距离之和所对应的第五目标候选区域，例如前10名，保存符合要求的第五目标候选区域，即得到第六目标候选区域。

进一步地，步骤d中，所述选定区域坐标的计算过程为：

将距离第五目标候选区域最近的道路运输交通枢纽、铁路运输的火车货站(或者高速公路路口)、机场作为选定区域。

在选定区域选取过程中，需排除上文涉及的自然灾害危险区域，交通枢纽点的选取需为两条高速的交叉点且至少有2个方向且判定为非地震危险区方向。

先以第五目标候选区域为高速公路路口为圆心的区域为例，计算第五目标候选区域到达各选定区域坐标的道路距离，具体如下所示：

a＝6378137m；

b＝6356752m；

y＝K(L-L

Z＝min∑

决策变量：

约束条件：

∑

∑

∑

符号用于研究的变量意义：

d为路网地理信息各点距离之和；

Q为路网地理信息数量之和；

q为路网信息点；

Z为目标函数，即最短距离之和；

j为第五目标候选区域，j∈J，J为第五目标候选区域集合；

h

f为交通枢纽，f∈F，F为交通枢纽集合

g为机场，g∈G，G为机场集合

p为火车货站，p∈P，P为火车货站集合

d

步骤f：根据所述研究区域的城市规划对所述第六目标候选区域进行筛选，并将符合所述研究区域的城市规划约束的第六目标候选区域确定为最佳选址区域，具体如图5所示。

例如：基于保存的候选区，根据城市规划指标，估划出参考候选区域；在考虑最大限度远离火源的情况下，得到理想候选区。

所述城市规划指标为：

以保山市为例，根据保山市的DEM坡度数据，结合《城市用地竖向规划规范》要求，将机场净空限制、水文条件、矿产资源分布、直升机起降等制约救灾物资储备库建设的因素纳入其中，得出区域建设用地情况。

需要说明一点：在最佳选址区域确定过程中，可能在某一步输出的结果为空集，进而造成无法筛选出最佳选址区域，在此情况下，首先要降低烈度区域的地震烈度，若仍然无法筛选出最佳选址区域则需要更换遥感识别建筑物技术，通过其他识别技术获取人类活动区域，即更新受灾区域，若还是无法筛选出最佳选址区域则需要更换研究区域。

本发明实施例提供的一种基于7级以上地震的救灾物资储备库选址方法，根据地震危险区域预判，确定地震烈度为8级以上的烈度区域；根据遥感识别建筑物技术识别8级以上烈度区域内的建筑物，并作为受灾区域；根据震后救灾物资运输的直升机设备，确定直升机运送救灾物资时的最远安全距离。将研究区域内以高速公路入口和火车货站为圆心，以5km为半径的区域设为选址候选区，在直升机运送救灾物资时的最远安全距离的情况下，求出符合到达所有受灾区域距离的候选区；基于研究区历史洪灾统计数据，划分出符合达到频次要求的候选区；基于候选区所在市各县人均GDP的统计数据，得出候选区所在县在市中的排名，得出符合达到名次要求的候选区；在考虑便于直升机起降的情况下，检验求出符合地形指数要求的候选区。如果候选区能够符合上述要求，则输出各候选区到达交通便利指标的陆运距离。在考虑城市规划指标的情况下，将前20％带入规划图进行判断。本发明实施例能够获得基于7级以上地震救灾物资储备库所在位置，保证选址工作能够应对区域内7级以上地震救灾物资运输和救援。

实施例二

如图6所示，本发明实施例提供的一种基于7级以上地震的救灾物资储备库选址系统，包括：

地震风险区域确定模块1，用于根据研究区域的地震文献资料，确定所述研究区域内地震震级为7级以上的地震风险区域。

烈度区域确定模块2，用于根据所述地震风险区域和地震烈度衰减模型，确定所述研究区域内地震烈度为8级以上的烈度区域。

受灾区域生成模块3，用于采用遥感识别建筑物技术识别所述烈度区域内的建筑物以生成建筑物分布图，然后采用八连通方式对所述建筑物分布图进行处理，以得到受灾区域。

最远安全飞行距离和选址候选区域确定模块4，用于确定直升机运送救灾物资时的最远安全飞行距离以及确定所述研究区域内救灾物资储备库的选址候选区域。

最佳选址区域确定模块5，用于根据约束条件对所述选址候选区域进行筛选，以确定最佳选址区域；所述约束条件至少包括所述最远安全飞行距离和所述受灾区域之间的约束关系、以及所述研究区域的自然灾害、地形指数、交通条件和城市规划。

其中，所述地震风险区域确定模块1，具体包括：

地震风险候选区域单元，用于根据研究区域的地震科学研究的预判成果和统计规律，确定所述研究区域内的地震风险候选区域。

最大地震震级确定单元，用于确定每个所述地震风险候选区域对应的最大地震震级。

地震风险区域确定单元，用于将所述最大地震震级与设定阈值进行比较，并将大于或者等于所述设定阈值的所述最大地震震级对应的所述地震风险候选区域确定为地震震级为7级以上的地震风险区域；其中，所述设定阈值为7。

所述自然灾害至少包括洪灾；所述约束条件还包括人均GDP统计结果；所述最佳选址区域确定模块5，具体包括：

第二目标候选区域确定单元，用于根据所述最远安全飞行距离和所述受灾区域计算出多个第一目标候选区，并将符合所述选址候选区域的第一目标候选区域确定为第二目标候选区域。

第三目标候选区域确定单元，用于根据所述研究区域内的洪灾区域和8级以上的烈度区域对所述第二目标候选区域进行筛选，并将符合所述研究区域内的非洪灾区域和非8级以上的烈度区域的第二目标候选区域确定为第三目标候选区域。

第四目标候选区域确定单元，用于根据两年内所述研究区域的人均GDP统计结果对所述第三目标候选区域进行筛选，并将符合两年内所述研究区域的人均GDP统计结果的第三目标候选区域确定为第四目标候选区域。

第五目标候选区域确定单元，用于根据所述研究区域的地形指数对所述第四目标候选区域进行筛选，并将符合所述研究区域的地形指数约束的第四目标候选区域确定为第五目标候选区域。

第六目标候选区域确定单元，用于根据所述研究区域的交通条件对所述第五目标候选区域进行筛选，并将符合所述研究区域的交通条件约束的第五目标候选区域确定为第六目标候选区域。

最佳选址区域确定单元，用于根据所述研究区域的城市规划对所述第六目标候选区域进行筛选，并将符合所述研究区域的城市规划约束的第六目标候选区域确定为最佳选址区域。

本发明根据文献资料确定地震风险区域，设定风险区最大发生地震震级和破坏性烈度等级(Ⅷ级)；根据判定地震风险区域的边界位置，确定烈度等级扩展边界；基于研究区的地震烈度衰减模型，求出地震Ⅷ级烈度区域轴距，所述烈度区域呈椭圆形；根据烈度区域，利用遥感识别技术求烈度区内建筑物；通过八联通方式给定受灾区域标签；给定救灾物资运输直升机机型；在考虑到直升机安全飞行要求，求出直升机最远安全飞行距离；给定研究区内高速公路入口和火车货站为圆心5km范围内是选址的候选区；基于受灾区域标签和最远安全飞行距离；在考虑研究区洪灾统计结果和Ⅷ级烈度区域的情况下，统计出候选区是否符合要求；在考虑两年内县人均GDP统计结果的情况下，统计出候选区是否符合要求；在考虑地形指数的情况下，检验求出候选区是否符合要求；在考虑交通便利因素的情况下，求出候选区到达各选定区域坐标的道路距离；根据各之和排序；基于保存的候选区，根据城市规划指标，估划出参考候选区域；在考虑最大限度远离火源的情况下，最终得到理想候选区。本发明基于受灾区域和直升机最远安全距离，将灾害情况、经济情况、地形指数作为约束条件，交通便利性作为优选条件，城市规划因素作为参考条件，保证救灾物资储备库选址与现实申报选址标准的有效性。

本发明弥补了基于7级以上地震大型救灾物资储备库选址方法的空缺，保障了救灾物资储备库选址的需求。本发明通过有效约束条件和优选条件求解，再通过参考条件进行分析就可以得到救灾物资储备库的理想选址区域，可以有效节省因现场勘测花费的人力、物力和时间，并为当前研究空缺的救灾物资储备库选址问题提供方法，有助于应急管理人员建库选址工作的快速推进。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理为对本发明的限制。