CN103809215B - 一种厚冲积层矿区土体沉陷响应的获取方法 - Google Patents

Abstract

本发明提供一种厚冲积层矿区土体沉陷响应的获取方法，包括：获取厚冲积层矿区的地表沉陷参数；获取与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数；根据所述厚冲积层矿区的地表沉陷参数和所述土体下边界的参考参数计算得到所述厚冲积层矿区的土体内部沉陷特征参数，以对所述厚冲积层矿区的土体沉陷响应进行分析。本发明提供的厚冲积层矿区土体沉陷响应的获取方法，用于获取厚冲积层矿区土体内部的沉陷响应，以对厚冲积层矿区地表沉陷的机理进行研究。

Description

一种厚冲积层矿区土体沉陷响应的获取方法

技术领域

本发明涉及矿山开采技术，尤其涉及一种厚冲积层矿区土体沉陷响应的获取方法。

背景技术

在我国的煤矿开采区中，矿区地质从地表往下依次为冲积层和基岩，其中，基岩包括上覆岩层、煤层和底板岩层。按照矿区的地质特征将矿区分为不同的种类，以冲积层的厚度特征为例，通常将冲积层厚度小于50m的矿区称为薄冲积层矿区，将冲积层厚度大于100m的矿区称为厚冲积层矿区，将冲积层厚度大于300米的矿区称为巨厚冲积层矿区。在对煤层进行开采的过程中以及开采结束后，少了煤层的支撑作用，煤层上方的冲积层、上覆岩层会发生沉陷现象，沉陷现象对地表的建筑设施的稳定性和安全性有一定的影响，因此，需要对矿区的沉陷特征参数进行获取，以对沉陷特征进行分析，进而采取适当的手段以避免沉陷现象对建筑设施造成严重的影响。

目前所采用的沉陷特征参数的获取方法主要是针对薄冲积层矿区或无冲积层矿区来设定的，而对于厚冲积层矿区，由于其冲积层较厚，对薄冲积层矿区或无冲积层矿区所采用的参数获取方法并不能适用。目前，对厚冲积层矿区土体内部沉陷特征参数的获取方法主要采用室内实验和数值模拟等方式，这种方式仅能对矿区表面现象进行分析，例如获取到地表下沉系数以对地表沉陷的范围进行分析等，而不能对厚冲积层矿区中复杂的岩土体的变化特征进行监测。因此，获取厚冲积层矿区土体内部沉陷响应的特征参数，以对其地表沉陷响应进行分析成为本领域中亟待解决的一大难题。

发明内容

本发明提供一种厚冲积层矿区土体沉陷响应的获取方法，用于对厚冲积层矿区土体内部的沉陷特征参数进行获取，以对厚冲积层矿区的土体沉陷响应进行分析。

本发明实施例提供一种厚冲积层矿区土体沉陷响应的获取方法，包括：

获取厚冲积层矿区的地表沉陷参数；

获取与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数；

根据所述厚冲积层矿区的地表沉陷参数和所述土体下边界的参考参数计算得到所述厚冲积层矿区的土体内部沉陷特征参数，以对所述厚冲积层矿区的土体沉陷响应进行分析。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，获取与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数，包括：

确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区；

获取所述薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区，包括：

采集厚冲积层矿区基岩的钻孔柱状标本作为标准标本；

采集待识别薄冲积层矿区基岩的钻孔柱状标本作为待识别标本；

分别提取所述标准标本和待识别标本的岩性特征；

根据所述标准标本和待识别标本的岩性特征确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，提取所述标准标本的岩性特征，包括：

将所述标准标本分为n段岩层，沿着从基岩至冲积层的方向上，每段岩层之间的距离S_i=i（i=1,2，…，n）；

根据每一段岩层对应的岩性评价系数按照如下公式计算得到所述标准标本的岩性评价系数P_i，作为所述标准标本的岩性特征：

$P_i=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{M}_i{Q}_i}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{M}_i},$

其中，M_i表示n段岩层中第i段岩层的法线厚度，Q_i表示n段岩层中第i段岩层对应的岩性评价系数。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，提取所述待识别标本的岩性特征，包括：

将所述待识别标本分为n段岩层，沿着从基岩至冲积层的方向上，每段岩层之间的距离S_j=j（j=1,2，…，n）；

根据每一段岩层对应的岩性评价系数按照如下公式计算得到所述待识别标本的岩性评价系数P_j，作为所述待识别标本的岩性特征：

$P_j=\frac{\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{M}_j{Q}_j}{\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{M}_j},$

其中，M_j表示n段岩层中第j段岩层的法线厚度，Q_j表示n段岩层中第j段岩层对应的岩性评价系数。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，根据所述标准标本和待识别标本的岩性特征确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区，包括：

根据所述标准标本和待识别标本的岩性评价系数计算所述标准标本和待识别标本的相似度；

将所述相似度与标准值进行比较，将所述相似度与标准值的差值最小的待识别标本对应的薄冲积层矿区作为与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，根据所述标准标本和待识别标本的岩性评价系数计算所述标准标本和待识别标本的相似度，包括：

根据所述标准标本和待识别标本的岩性评价系数计算得到所述标准标本和待识别标本之间的隶属度；

根据所述隶属度计算得到所述标准标本和待识别标本的相似度。

如上所述的厚冲积层矿区土体沉陷响应的获取方法，根据所述标准标本和待识别标本的岩性评价系数计算得到所述标准标本和待识别标本之间的隶属度,包括：

根据所述标准标本和待识别标本的岩性评价系数按照如下公式计算得到所述标准标本和待识别标本之间第i段（i=1,2，…，n）岩层的隶属度u_k：

$u_i=1.0-\frac{|P_i-P_j|}{P_i}.$

如上所述的厚冲积层矿区土体沉陷响应的获取方法，根据所述隶属度计算得到所述标准标本和待识别标本的相似度，包括：

根据所述隶属度按照如下公式计算所述标准标本与m个待识别标本的相似度C_k（k=1,2，…，m）：

$C_k=\frac{\mathrm{\Σ}{u}_i}{n}.$

本发明实施例提供的技术方案，通过选取一个与厚冲积层矿区的基岩最为相似的薄冲积层矿区，分别获取厚冲积层矿区和薄冲积层矿区的地表沉陷参数，并将薄冲积层矿区的地表沉陷参数作为厚冲积层矿区冲积层的土体下边界条件，将厚冲积层矿区的地表沉陷参数作为冲积层的土体上边界条件，根据冲积层土体上边界条件和下边界条件对厚冲积层矿区的土体沉陷响应特征及地表沉陷特殊性进行计算和分析，得到厚冲积层矿区地表沉陷响应特征参数，进而分析厚冲积层矿区的地表沉陷机理，以对厚冲积层矿区地表沉陷的特殊性进行全面的解释，解决了现有技术中的一大难题。

附图说明

图1为本发明实施例提供的厚冲积层矿区土体沉陷响应的获取方法的流程图；

图2为矿区地层分布示意图；

图3为上覆岩层中各岩层的分布示意图；

图4为厚冲积层地表下沉的结构示意图。

具体实施方式

图1为本发明实施例提供的厚冲积层矿区土体沉陷响应的获取方法的流程图。本实施例提供一种厚冲积层矿区土体沉陷响应的获取方法，能够对厚冲积层矿区土体的沉陷响应及地表沉陷特殊性机理进行***研究，以对厚冲积层矿区地表沉陷的特殊性进行全面的解释，如图1所示，该方法可以包括：

步骤10、获取厚冲积层矿区的地表沉陷参数。

步骤20、获取与厚冲积层矿区的基岩具有相似特征的薄冲积层矿区的地表沉陷参数，作为厚冲积层矿区土体下边界的参考参数。

步骤30、根据厚冲积层矿区的地表沉陷参数和土体下边界的参考参数计算得到厚冲积层矿区的土体内部沉陷特征参数，以对厚冲积层矿区的土体沉陷响应进行分析。

其中，获取厚冲积层矿区的地表沉陷参数具体可在厚冲积层矿区的地表建立观测站，对厚冲积层矿区的地表进行测量，得到地表沉陷参数，具体可按照现有技术中常用的技术方案来实现，例如在开采影响范围内的地表，设定一系列相互联系的观测点，在采矿工作过程中，定期测量各观测点的空间坐标以及相对位置的变化量，通过对实测数据的反演分析得到实测的地表沉陷参数，作为厚冲积层土体上边界条件。

然后确定一个与厚冲积层矿区的基岩特征较相似的薄冲积层矿区，本实施例中，可确定一个薄冲积层矿区，也可以确定一个无冲积层矿区，该薄冲积层矿区指的是冲积层的厚度小于50米的矿区，无冲积层矿区指的是没有冲积层的矿区。本实施例仅以薄冲积层矿区为例对技术方案进行详细说明，本领域技术人员可以根据本实施例提供的技术方案也可以选取无冲积层矿区来实现。无论是厚冲积层矿区还是薄冲积层矿区的冲积层均为土体结构。

在确定出一个与厚冲积层矿区的基岩特征较相似的薄冲积层矿区之后，获取该薄冲积层矿区的地表沉陷参数，作为厚冲积层矿区冲积层土体下边界（即：冲积层与基岩的交界面）的参考参数。具体可在薄冲积层矿区的地表建立观测站，对薄冲积层矿区的地表进行测量，得到地表沉陷参数，具体实现方式可与在厚冲积层的地表建立观测站的方式相同，根据薄冲积层矿区表面各观测点的空间坐标以及相对位置的变化量，通过反演分析出实测地表沉陷参数，即为薄冲积层矿区的地表沉陷参数，以此作为厚冲积层土体的下边界条件。

图2为矿区地层分布示意图，如图2所示，矿区地层从地表向下依次分为冲积层土体和基岩，其中，基岩分为上覆岩层、煤层和底板岩层。由于薄冲积层矿区的冲积层相对很薄，若其基岩与厚冲积层矿区的基岩特征相类似，则薄冲积层矿区的地表变化特性与厚冲积层矿区基岩上表面（即冲积层下表面）的变化特性相类似，则可以将薄冲积层矿区的地表沉陷参数作为厚冲积层矿区冲积层下表面的沉陷参数（即冲积层土体的下边界条件），且与厚冲积层矿区的地表沉陷参数（即冲积层土体的上边界条件）相结合进行分析，以得到厚冲积层矿区的土体内部沉陷响应特征，以实现对厚冲积层矿区的土体沉陷响应进行分析。

上述确定与厚冲积层矿区的基岩特征较相似的薄冲积层矿区，采用的方式有很多中，本实施例提供一种可以实现的方式：

首先，采集厚冲积层矿区基岩的钻孔柱状标本作为标准标本，采集待识别薄冲积层矿区基岩的钻孔柱状标本作为待识别标本，然后分别提取标准标本和待识别标本的岩性特征，根据标准标本和待识别标本的岩性特征确定与厚冲积层矿区的基岩具有相似特征的薄冲积层矿区。

本实施例提供的技术方案是对基岩特性的比较，其中，采集厚冲积层矿区的钻孔柱状标本的方式可以在厚冲积层矿区的地表向下进行钻孔，然后去除采集得到的柱状标本中的冲积层，余下的即为基岩，作为标准标本。

对于薄冲积层矿区的确定过程，可选取多个薄冲积层矿区进行标本采集，分别与厚冲积层矿区的标准标本进行比较，确定一个最相似的。具体可在各薄冲积层矿区的地表向下进行钻孔，然后去除采集得到的柱状标本中的冲积层，余下的即为基岩，作为待识别标本。

分别提取标准标本和待识别标本的岩性特征，具体可以采用如下技术方案来实现：

（1）提取标准标本的岩性特征，可将标准标本分为n段岩层，沿着从基岩至冲积层的方向上，以煤层的上表面为原点，向上每段岩层之间的距离S_i=i（i=1,2，…，n）。然后根据每一段岩层对应的岩性评价系数按照如下公式计算得到标准标本的岩性评价系数P_i，作为标准标本的岩性特征：

$P_i=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{M}_i{Q}_i}{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{M}_i},$

其中，M_i表示n段岩层中第i段岩层的法线厚度，Q_i表示n段岩层中第i段岩层对应的岩性评价系数，可参照表一来对Q_i进行设定。例如当第2段岩层为硬粘土质片岩时，则Q₂可以取值为0.4。

表一岩性评价系数参考值

岩石名称	岩性评价系数
		硬的石灰岩、硬砂岩、硬大理石、不硬的花岗岩	0.0
较硬的石灰岩、砂岩、大理石	0.05
		普通砂岩、铁矿石	0.1
砂质页岩、片状砂岩	0.2
		硬粘土质片岩、不硬的砂岩和石灰岩、软砾岩	0.4
各种页岩（不坚硬的）、致密泥灰岩	0.6
		软页岩、很软石灰岩、无烟煤、普通泥灰岩	0.8
破碎页岩、烟煤、粘土（致密的）	0.9
		软砂质粘土、黄土、软砾石、松散砂层	1.0

（2）提取待识别标本的岩性特征，可以将待识别标本分为n段岩层，沿着从基岩至冲积层的方向上，以煤层的上表面为原点，向上每段岩层之间的距离S_j=j（j=1,2，…，n）。根据每一段岩层对应的岩性评价系数按照如下公式计算得到待识别标本的岩性评价系数P_j，作为待识别标本的岩性特征：

$P_j=\frac{\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{M}_j{Q}_j}{\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{M}_j},$

其中，M_j表示n段岩层中第j段岩层的法线厚度，Q_j表示n段岩层中第j段岩层对应的岩性评价系数，也可以参照表一来设定。

上述对标准标本和待识别标本按照等差级数进行分段，是考虑到上覆岩层中不同岩层的下沉作用不同，如图3所示，图3为上覆岩层中各岩层的分布示意图。上覆岩层从上至下依次可以分为松散冲积层带、弯曲带、裂缝带和冒落带，采集到的标本中，距离煤层越近的岩层受开采影响越大，对裂缝带内岩体的影响越大。

在提取到标准标本和待识别标本的岩性特征之后，根据标准标本和待识别标本的岩性特征确定与厚冲积层矿区的基岩具有相似特征的薄冲积层矿区，具体可以根据标准标本的岩性评价系数P_i和待识别标本的岩性评价系数P_j计算标准标本和待识别标本的相似度。

计算相似度的方式有很多种，例如采用如下方式：根据标准标本的岩性评价系数P_i和待识别标本的岩性评价系数P_j计算得到标准标本和待识别标本之间的隶属度，如：按照如下公式计算得到标准标本和待识别标本之间第i段（i=1,2，…，n）岩层的隶属度u_i：

$u_i=1.0-\frac{|P_i-P_j|}{P_i}.$

其中，i=j。然后根据隶属度u_i计算得到标准标本和各待识别标本之间的相似度，具体可按照如下公式计算标准标本与m个待识别标本的相似度C_k（k=1,2，…，m）：

$C_k=\frac{\mathrm{\Σ}{u}_i}{n}.$

将相似度C_k与标准值进行比较，该标准值可以为1.00。将相似度C_k与1.00的差值最小的待识别标本对应的薄冲积层矿区视为与厚冲积层矿区的基岩做接近，作为与厚冲积层矿区的基岩具有相似特征的薄冲积层矿区。

根据上述内容提供的技术方案，在确定最接近厚冲积层矿区基岩的薄冲积层矿区之后，在薄冲积层矿区的地表建立观测站，获取地表沉陷参数。然后将薄冲积层矿区的地表沉陷参数作为厚冲积层矿区中冲积层土体下边界的沉陷参数，与厚冲积层矿区的地表沉陷参数相结合进行分析，以得到厚冲积层矿区的土体沉陷响应特征。

具体可参照图4，图4为厚冲积层地表沉陷的结构示意图，图4中将冲积层和基岩分开，仅用于对本实施例的技术方案进行说明。在实际矿区中，冲积层的下表面即为基岩的上表面（也即冲积层土体与基岩的交界面）。在煤层被采出后形成采空区，少了煤层的支撑作用，上覆岩层发生冒落、离层及弯曲，然后逐渐向上传递至冲积层土体，使得冲积层土体下沉。冲积层土体发生下沉会导致土体跟随下沉、失水固结下沉以及土体压密下沉三方面的响应。且厚冲积层土体以荷载的形式向下传递直至采空区，由于作用在基岩及采空区的应力增加，受到冲积层荷载影响的采空区冒落带内破碎岩体更加压实，岩石碎胀性减小，裂缝带内沿层裂隙及离层趋于闭合，置换出来的空间传递到地表，引起岩土体协同作用下沉响应，增加地表的下沉量。

根据冲积层上表面和下表面的沉陷参数，可以按照土体的物理力学性质，用力学的方法计算得到各土层的下沉情况。厚冲积层矿区的土体内部沉陷特征参数可以为多种参数，具体可参照现有技术中对薄冲积层矿区或无冲积层矿区进行分析所采用的各项参数，其具体的分析方法也可参照现有技术中常用的方法来实现。例如可采用现有技术中常用的概率积分法根据薄冲积层矿区的实测沉陷参数计算厚冲积层矿区的地表下沉值（相当于薄冲积层矿区的采深变大），得到厚冲积层矿区的土体跟随下沉；根据厚冲积层矿区地表实测下沉和所计算的跟随下沉，求其差值得到相应的冲积层土体沉陷响应曲线。然后根据冲积层下表面的沉陷参数和土体跟随下沉参数之间的差值，得到土体附加下沉参数以及相应的曲线，根据该附加下沉参数以及曲线可以对失水固结下沉、土体压密下沉以及岩土协同作用下沉的情况进行分析。

本实施例提供的技术方案，通过选取一个与厚冲积层矿区的基岩最为相似的薄冲积层矿区，分别获取厚冲积层矿区和薄冲积层矿区的地表沉陷参数，并将薄冲积层矿区的地表沉陷参数作为厚冲积层矿区冲积层的土体下边界条件，将厚冲积层矿区的地表沉陷参数作为冲积层的土体上边界条件，根据冲积层土体上边界条件和下边界条件对厚冲积层矿区的土体沉陷响应特征及地表沉陷特殊性进行计算和分析，得到厚冲积层矿区土体内部沉陷响应特征参数，进而分析厚冲积层矿区的地表沉陷机理，以对厚冲积层矿区地表沉陷的特殊性进行全面的解释，解决了现有技术中的一大难题。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (6)

1.一种厚冲积层矿区土体沉陷响应的获取方法，其特征在于，包括：

获取厚冲积层矿区的地表沉陷参数；

获取与所述厚冲积层矿区的基岩具有相似岩性特征的薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数；

根据所述厚冲积层矿区的地表沉陷参数和所述土体下边界的参考参数计算得到所述厚冲积层矿区的土体内部沉陷特征参数，以对所述厚冲积层矿区的土体沉陷响应进行分析；

其中，获取与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数，包括：

确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区；

获取所述薄冲积层矿区的地表沉陷参数，作为所述厚冲积层矿区土体下边界的参考参数；

其中，确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区，包括：

采集厚冲积层矿区基岩的钻孔柱状标本作为标准标本；

采集待识别薄冲积层矿区基岩的钻孔柱状标本作为待识别标本；

分别提取所述标准标本和待识别标本的岩性特征；

根据所述标准标本和待识别标本的岩性特征确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区；

其中，提取所述标准标本的岩性特征，包括：

将所述标准标本分为n段岩层，沿着从基岩至冲积层的方向上，每段岩层之间的距离S_i＝i，其中i＝1,2，…，n；

根据每一段岩层对应的岩性评价系数按照如下公式计算得到所述标准标本的岩性评价系数P_i，作为所述标准标本的岩性特征：

$P_i=\frac{\underset{i=1}{\overset{n}{\Σ}}{M}_i{Q}_i}{\underset{i=1}{\overset{n}{\Σ}}{M}_i},$

其中，M_i表示n段岩层中第i段岩层的法线厚度，M_i＝S_i·cosα,α表示第i段岩层的岩层倾角，Q_i表示n段岩层中第i段岩层对应的岩性评价系数。

2.根据权利要求1所述的厚冲积层矿区土体沉陷响应的获取方法，其特征在于，提取所述待识别标本的岩性特征，包括：

将所述待识别标本分为n段岩层，沿着从基岩至冲积层的方向上，每段岩层之间的距离S_j＝j，其中j＝1,2，…，n；

根据每一段岩层对应的岩性评价系数按照如下公式计算得到所述待识别标本的岩性评价系数P_j，作为所述待识别标本的岩性特征：

$P_j=\frac{\underset{j=1}{\overset{n}{\Σ}}{M}_j{Q}_j}{\underset{j=1}{\overset{n}{\Σ}}{M}_j},$

其中，M_j表示n段岩层中第j段岩层的法线厚度，M_j＝S_j·cosα,α表示第j段岩层的岩层倾角，Q_j表示n段岩层中第j段岩层对应的岩性评价系数。

3.根据权利要求2所述的厚冲积层矿区土体沉陷响应的获取方法，其特征在于，根据所述标准标本和待识别标本的岩性特征确定与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区，包括：

根据所述标准标本和待识别标本的岩性评价系数计算所述标准标本和待识别标本的相似度；

将所述相似度与标准值进行比较，将所述相似度与标准值的差值最小的待识别标本对应的薄冲积层矿区作为与所述厚冲积层矿区的基岩具有相似特征的薄冲积层矿区。

4.根据权利要求3所述的厚冲积层矿区土体沉陷响应的获取方法，其特征在于，根据所述标准标本和待识别标本的岩性评价系数计算所述标准标本和待识别标本的相似度，包括：

根据所述标准标本和待识别标本的岩性评价系数计算得到所述标准标本和待识别标本之间的隶属度；

根据所述隶属度计算得到所述标准标本和待识别标本的相似度。

5.根据权利要求4所述的厚冲积层矿区土体沉陷响应的获取方法，其特征在于，根据所述标准标本和待识别标本的岩性评价系数计算得到所述标准标本和待识别标本之间的隶属度,包括：

根据所述标准标本和待识别标本的岩性评价系数按照如下公式计算得到所述标准标本和待识别标本之间第i段岩层的隶属度u_i，其中i＝1,2，…，n：

$u_i=1.0-\frac{\left|P_i-P_j\right|}{P_i}.$

6.根据权利要求5所述的厚冲积层矿区土体沉陷响应的获取方法，其特征在于，根据所述隶属度计算得到所述标准标本和待识别标本的相似度，包括：

根据所述隶属度按照如下公式计算所述标准标本与m个待识别标本的相似度C_k，其中k＝1,2，…，m：

$C_k=\frac{\Σu_i}{n}.$