基于建造大数据的盾构隧道施工标准工时测算及评分方法与流程

1.本发明涉及工时测算及评价指标建立领域，具体涉及一种基于建造大数据的盾构隧道施工标准工时测算及评分方法。


背景技术：

2.在地铁的隧道盾构过程中，施工部门会面临着多样复杂的工作环境，盾构施工的进度受到多个因素的影响，如地质条件、设备状态、政策变动等。目前企业大多数是根据工作经验或行业内的专家评分进行评价指标测算，其评价结果受到建造主观因素的影响，且对于评分的工作量也要求较高。
3.现有的地质分类方式主要是根据地质纵面图和地层分类标准，对地铁隧道中每种地质进行归类。该方法划分出的地层类别繁杂，各类别具有不同的盾构环境和施工条件，不利于制定工作标准和规划工程进度。
4.此外，工时测算一般依据企业相关人员的已有经验和专业知识，对不同地质条件下的标准工时进行粗略地制定，但难以从数据的角度完全验证及支撑最后选定的标准工时，且该方法难以考虑到施工现场的实际客观条件。
5.目前，研究的评价指标建立方法主要有模糊综合评价法、fhw法等。模糊综合评价法需制定出每个指标的评价标准和加分权重，计算复杂，对指标权重矢量的确定主观性较强；fhw法通过组织评分专家进行独立评价和交流后评价，评价过程繁锁且需要多人参与，建立评价指标的成本较高。


技术实现要素：

6.本发明的目的是提供一种基于建造大数据的盾构隧道施工标准工时测算及评分方法，以建立客观的评价体系，给盾构施工提供科学依据，可用于对不同施工队进行科学合理的考核和选拔。
7.为了实现上述任务，本发明采用以下技术方案：
8.一种基于建造大数据的盾构隧道施工标准工时测算及评分方法，包括以下步骤：
9.步骤1，对盾构掘进参数产生影响的不同外界因素分别与掘进参数进行相关性检验，确定对掘进参数的状态产生主要影响的外界因素为地质因素；
10.步骤2，对应用相同盾构法的隧道区间对应的原始掘进参数进行标准化，在此基础上对其中的异常值进行界定及剔除；对剔除异常值后的掘进参数对应的地质类别符号进行统一；
11.步骤3，抽取掘进参数中的单环推进时间、停机时间以及单位时间内的推进环数作为聚类特征，利用聚类算法，得到各类掘进参数及其对应的地质类别最终的距离、分类内容以及分类数目；同时通过分析其显著性可以判断该聚类结果是否合理，若结果不显著，则需进一步调整聚类数目，再次聚类；
12.步骤4，首先根据聚类结果，将得到的每一种地质分类中的所有单环总时间数据进
行描述，得到每一分类中单环总时间的各项特征参数，包括均值、标准差、中位数、间距；其次，运用直方图对每一分类中的单环总时间进行统计并得到统计分布；通过分析和比较各项分布特征，选择单环总时间的ω百分位数作为该分类的单环标准工时；
13.以制定的单环标准工时为分界，分别在每一类地层中划分出总工时的及格区间和不及格区间，其中及格/不及格样本为每一分类及格/不及格区间中的各环掘进参数及其对应的地质类别；
14.步骤5，建立逻辑回归模型得到每个变量的影响系数，进而利用评分卡转换得到最终的评分公式；运用评分公式计算出每环盾构掘进的实际工时对应的分数，通过分数反映不同施工队的管理水平并将其作为考核选拔的参考依据；其中所述变量为单环推进时间、单环停机时间、单环总时间、环数和地质类别。
15.进一步地，所述确定对掘进参数的状态产生主要影响的外界因素为地质因素，包括：
16.先设定β＝{β1，β2，β3，...，β
n1
}、γ＝{γ1，γ2，γ3，...，γ
n2
}两总体，其中β为地质因素中的岩性因素对应的掘进数据总体；γ为地质因素中土的软硬度因素对应的掘进数据总体；掘进数据总体包括单环推进时间、单环停机时间、单环总时间、环数；在两数据总体中分别抽取容量为n1，n2的样本，且设两样本独立，并假定n1≠n2，现将这n1 n2的观察值放在一起，按自小到大的次序排列，每个实验值在序列中的次序叫做该值的秩，然后将属于第β组数据的秩相加，其和记为r1，同理可求得第γ组数据的r2，则r1和r2是离散型随机变量，且有是离散型随机变量，且有
17.当两个样本容量都大于10时，秩和t的分布接近于正态分布，采用z检验，其基本公式为：
18.t为容量较小样本的秩和；
19.最后对于给定的显著性水平，确定z值的临界值，比较临界值和z值的大小，判断两总体是否有显著性差别；若结果显著，则说明该地质因素的变化会影响掘进参数水平的状态，将其作为评分模型的变量。
20.进一步地，对原始掘进参数进行标准化采用z-score方法，对异常值进行界定及剔除采用三倍标准差法；在地质类别符号进行统一的过程中，将异常值进行界定及剔除后剩余的掘进参数对应的地质类别符号逐一修改，统一同类地质的标注形式，确定所有地质类别不存在同类而不同标注符号的情况。
21.进一步地，所述聚类算法为：
22.输入：n个环样本的特征向量{x
(1)
，x
(2)
，
…
，x
(n)
}，x(i)为第i环样本，环样本为步骤2处理后得到的掘进参数及其对应的地质类别；
23.输出：分类个数k和分类结果；
24.步骤a，指定需要划分的类的个数k；
25.步骤b，随机初始化k个聚类中心μj，j＝1，2，...，k；
26.步骤c，将类划分为t＝1，2，
…
，k，c
t
为同一类特征向量的集合；
27.步骤d，对于i＝1，2，
…
m，计算样本x(i)和各个质心向量μj(j＝1，2，
…
，k)的距离：将x(i)标记最小的为d
ij
所对应的类别λi，此时更新，此时更新
28.步骤e，对于i＝1，2，
…
k，对cj中所有的样本点重新计算新的质心中所有的样本点重新计算新的质心
29.步骤f，如果所有的k个质心向量都没有发生变化，则转到步骤c，有变化则转到步骤g；
30.步骤g，输出聚类结果，包括各分类的单环推进时间、单环停机时间、环数及单环总时间这些掘进参数以及其对应的地质类别。
31.进一步地，每一类地层中及格样本与不及格样本的比例为ω：1-ω；若i为及格样本，则及定义为yi＝1；若i为不及格样本，则定义为yi＝0。
32.进一步地，所述逻辑回归模型表示为：
33.设条件概率p(y＝1|v)＝p为根据观测量相对于某事件v发生的概率；那么逻辑回归模型可以表示为：
[0034][0035]
其中，v＝(v1，v2，
…
，v5)，向量vi为掘进参数：单环推进时间、单环停机时间、单环总时间、环数及地质类别，i＝1，2，3，4，5；g(v)＝w0 w1v1
…
w5v5，w0，w1，w2，w3，w4，w5为系数；那么在v条件下y不发生的概率为：
[0036][0037]
比值称为事件的发生比，记为odds；对其取对数得到：
[0038][0039]
进一步地，评分公式为：score＝a-b*ln(odds)，其中a与b是常数；评分卡的评分范围是0～1000；
[0040]
通过以下两个式子计算得出a、b的值，得到最终的评分公式。
[0041]
1000＝a-b*ln(odds
第一名
)
[0042]
0＝a-b*ln(odds
最后一名
)
[0043]
其中，odds
第一名
、odds
最后一名
为计算出的排名第一、最后的odds。
[0044]
与现有技术相比，本发明具有以下技术特点：
[0045]
1.本发明面向地质特征，合理制定盾构掘进的标准工时。地铁盾构隧道所经过的地层类别较为复杂，具有不同程度的施工难度，直接影响到工人的施工进度和效率。利用k-means等聚类方法，对不同岩性和软硬度的地层进行划分，得到清晰明确的地层聚类结果。在此基础上，可进一步分析各个地层类别的掘进数据分布，制定出对应的掘进标准时长，为盾构建设的进度规划及施工绩效考核提供考察依据。
[0046]
2.本发明面对多种复杂环境，建立科学合理的盾构施工评价指标。在明确地质、设
备等客观影响因素的基础上，建立施工绩效评价指标，引入sigmoid函数和修改过的信用评分卡，根据不同地质特征得到相对应的评价公式，更加有针对性地考核施工队的盾构工作进度和绩效，帮助企业对多个施工队进行精准科学地考核。
附图说明
[0047]
图1为本发明方法的流程示意图；
[0048]
图2为实施例中聚类算法流程图；
[0049]
图3为实施例中单环总时间直方图；
[0050]
图4为实施例中及格样本分布图；
[0051]
图5为实施例中不及格样本分布图。
具体实施方式
[0052]
参见附图，本发明提出了一种基于建造大数据的盾构隧道施工标准工时测算及评分方法，包括以下步骤：
[0053]
步骤1，确定对掘进参数的状态产生主要影响的客观因素
[0054]
对掘进参数与外界影响因素进行相关性检验，确定对掘进参数的状态产生主要影响的客观因素；本实施例中该因素为地质因素，不同地质会对掘进施工产生不同的难度，导致掘进速率水平发生变化。
[0055]
在地铁建设工程中，对盾构掘进参数产生影响作用的外界因素有多种。首先是盾构施工环境，包括地质因素(地质类别，例如地质的岩性和土的软硬度)、天气因素、经过的外界建筑物等客观因素；第二是施工队的管理水平，包括工程进度、设备材料供给等主观因素；最后是政策变动影响。这些因素的改变都可能对盾构施工的进度产生不同程度的影响。从大数据中得到的规律可发现，相关影响因素难以保证服从正态分布。因此，需将以上的几种外界影响因素与掘进参数做相关性分析和秩检验，分析数据之间的相关性关系，进而确定对掘进参数水平产生主要客观影响的外界影响因素。
[0056]
本发明中采用的相关性检验方法如下：
[0057]
先设定β＝{β1，β2，β3，...，β
n1
}、γ＝{γ1，γ2，γ3，...，γ
n2
}两总体，其中β为地质影响因素中的岩性因素对应的掘进数据总体；γ为相同影响因素中土的软硬度因素对应的掘进数据总体。掘进数据总体包括单环推进时间、单环停机时间、单环总时间、环数。在两总体中分别抽取容量为n1，n2的样本，且设两样本独立，并假定n1≠n2，现将这n1 n2的观察值放在一起，按自小到大的次序排列，每个实验值在序列中的次序叫做该值的秩，然后将属于第β组数据的秩相加，其和记为r1，同理可求得第γ组数据的r2。
[0058]
显然，r1和r2是离散型随机变量，且确
[0059]
当两个样本容量都大于10时，秩和t的分布接近于正态分布，因此可以用z检验，其基本公式为：
[0060]
t为容量较小样本的秩和。
[0061]
最后对于给定的显著性水平，确定z值的临界值，比较两者大小，判断两总体是否
有显著性差别。若结果显著，则说明该外界因素的变化会影响掘进参数水平的状态，可将其作为评分模型的变量；若不显著，则忽略该因素对掘进参数水平的影响，重新选择影响因素。
[0062]
由于天气、建筑物及政策因素影响作用不稳定，地质因素对盾构的影响较为明显且稳定；通过上述相关性检验方法，确定地质因素作为掘进速度的客观影响因素。
[0063]
步骤2，对掘进参数进行预处理
[0064]
本步骤先对应用相同盾构法的隧道区间对应的原始掘进参数进行标准化，在此基础上采用三倍标准差法对其中的异常值进行界定及剔除；另外，由于步骤1确定掘进参数的外界影响因素为地质因素，因此需对剔除异常值后的掘进参数对应的地质类别符号进行统一，避免由于标注符号的不一致导致多个相同地质类别在聚类时被识别为不同类别，出现错误的聚类结果。
[0065]
2.1数据标准化处理
[0066]
z-score也叫标准分数，它是以标准差为尺子去度量某一原始分数偏离平均数的距离，从而确定这一数据在全体数据中的位置。在本方案中，为了准确界定掘进参数中的离群值所在范围，利用z-score法将各环掘进参数进行标准化，其转化的公式为：
[0067][0068]
式中，x为原始掘进数据，包括单环推进时间、单环停机时间及单环总时间，为平均数，σ为标准差。
[0069]
标准化后，原始数据改为均值为0，方差为1的分布。
[0070]
2.2异常值的界定及剔除
[0071]
针对步骤2.1处理得到的掘进数据中，存在部分异常值，包括过大值、过小值、零值及空值。这些异常值的存在会改变接下来的相关数据处理结果，影响对数据处理结果的判断和分析，因此需要将其剔除。本方案采用的方法是工业工程中的三倍标准差法。筛选标准化后的两组z数据，将不在-3σ到3σ范围内的盾构掘进数据界定为异常值，即剔除在-3至3范围之外的异常值。经计算，在各个区间中，被剔除的异常值所占原有总数据量的比例均在2.5％左右，说明异常值的界定范围选取恰当。
[0072]
2.3地质标注符号的统一化
[0073]
盾构施工每一环的工时与地质状况相关性较大。由于施工过程中的记录不当，收集到的地质类别数据中可能存在以下几种错误：同类地层所在不同环号的地质符号标注不统一；同类地层所在不同环号的多种地质符号标注顺序不同；同类地层所在同一环号的同类地质符号重复标注；同类地层所在不同环号的地质符号记录形式不一致等。
[0074]
以上出现的掘进参数记录错误现象，导致同一类地层由于标注符号不统一而被划分到多个类别中，出现虚假的地质类别数。因此，需要将步骤2.3处理后剩余的掘进参数所对应的地质类别符号逐一修改，统一同类地质的标注形式，确定所有地质类别不存在同类而不同标注符号的情况，避免在步骤3地质聚类过程中出现错误的聚类结果。
[0075]
步骤3，对预处理后的掘进参数进行聚类分析
[0076]
根据地质领域知识对完成预处理的地质类别及其对应的掘进参数进行分类，在此基础上逐一对每个分类采用k-means算法，聚合掘进参数水平相同的所有地质类别及掘进
参数，得到最终的掘进参数聚类结果。
[0077]
3.1确定聚类特征
[0078]
抽取掘进参数中的单环推进时间、停机时间以及单位时间内的推进环数作为聚类特征，能够体现盾构掘进速度的状态。在本发明中，已将其与地质因素、政策变动等几个外界因素进行相关性检验，检验结果发现三个聚类特征均与这几个影响因素存在显著性差别，说明盾构掘进速度的状态会随着地质因素变化、施工队的管理水平不同而发生改变，其中地质因素为客观影响因素，管理水平为主观影响因素。本发明是在相同客观影响因素的作用下，对施工水平进行评价，因此可将单环的推进时间、停机时间和单位时间内的推进环数作为聚类特征，对地质类别进行聚类，在分析聚类结果的过程中逐步去除偶然性因素，从而得出具有真实性、客观性的最终结果。
[0079]
3.2聚类分析
[0080]
盾构掘进速度分类是盾构施工标准工时制定研究的基本工作，聚类方法的合理性会影响最终结果的准确性。目前对盾构施工的研究，都是依据盾构掘进时间的特征对掘进速度进行分类，这种分类方法无法很好地保证后续盾构标准工时制定的正确性，因此要得到合理的掘进速度分类，则需要研究出一套面向地质特征的盾构掘进速度分类方法。本发明在初步确定适合实际情况的聚类数目，利用改进的k-means聚类算法，得到各类掘进参数及其对应的地质类别最终的距离，分类内容以及分类数目k。同时通过分析其显著性可以判断该聚类结果是否合理，若结果不显著，则需进一步调整聚类数目，再次聚类直至结果具有较高的显著性。
[0081]
(1)k-means聚类
[0082]
k-means聚类又称快速聚类，属于覆盖数据划分聚类算法。在k-means聚类结果中，每个样本点都属于唯一的类别。k-means聚类算法的变量是数值变量，其采用划分原则进行聚类。划分原则，就是将样本空间随机划分为多个区域，根据定义的中心距离，将全部样本点分配到“封闭”域中，形成初始聚类。通过反复计算各类别之间的中心距离，对样本点进行分类，得到最终的结果。
[0083]
(2)改进的k-means聚类的算法为：
[0084]
输入：n个环样本的特征向量{x
(1)
，x
(2)
，
…
，x
(n)
}，x(i)为第i环样本，环样本为步骤2预处理后得到的掘进参数及其对应的地质类别。
[0085]
输出：分类个数k和分类结果
[0086]
步骤a：指定需要划分的类的个数k
[0087]
步骤b：随机初始化k个聚类中心μj，j＝1，2，...，k
[0088]
步骤c：将类划分为t＝1，2，
…
，k，c
t
为同一类特征向量的集合。
[0089]
步骤d：对于i＝1，2，
…
m，计算样本x(i)和各个质心向量μj(j＝1，2，
…
，k)的距离：将x(i)标记最小的为d
ij
所对应的类别λi，此时更新，此时更新
[0090]
步骤e：对于i＝1，2，
…
k，对cj中所有的样本点重新计算新的质心中所有的样本点重新计算新的质心
[0091]
步骤f：如果所有的k个质心向量都没有发生变化，则转到步骤c，有变化则转到步骤g；
[0092]
步骤g：输出聚类结果，包括各分类的单环推进时间、单环停机时间、环数及单环总时间等掘进参数以及其对应的地质类别。
[0093]
步骤4，制定标准工时和及格区间
[0094]
4.1制定标准工时
[0095]
关于单环标准工时的测算，首先根据聚类结果，将得到的每一种地质分类中的所有单环总时间数据进行描述，得到每一分类中单环总时间的各项特征参数，包括均值、标准差、中位数、间距等参数。其次，运用直方图对每一分类中的单环总时间进行统计并得到统计分布。通过分析和比较各项分布特征，选择单环总时间的ω百分位数作为该分类的单环标准工时；ω为预设参数。
[0096]
4.2定义及格区间
[0097]
以制定的单环标准工时为分界，分别在每一类地层中划分出总工时的及格区间和不及格区间，其中及格(不及格)样本为每一分类及格(不及格)区间中的各环掘进参数及其对应的地质类别。每一类地层中及格样本与不及格样本的比例为ω：1-ω。若i为及格样本，则定义为yi＝1；若i为不及格样本，则定义为yi＝0。
[0098]
步骤5，建立逻辑回归模型并进行评分
[0099]
本发明评价的对象为施工单位的施工速度效率，反映施工效率的直接指标是盾构的单环总时间，而主要影响施工快慢的因素为地质情况、设备因素和施工队的管理水平，地质情况和设备因素是客观方面的影响因素，因此实际上需要评价的对象其实是施工队的管理水平。
[0100]
步骤4制定了盾构隧道的单环标准工时，利用标准工时划分每一类掘进参数总体数据的及格区间和不及格区间，并定义及格样本不及格样本；在该步骤中，建立逻辑回归模型得到每个变量xi的影响系数b，进而利用评分卡转换得到最终的评分公式。运用该评分公式可计算出每环盾构掘进的实际工时对应的分数，通过分数反映不同施工队的管理水平并将其作为考核选拔的参考依据，具体如下：
[0101]
目前，在设备条件相同的前提下，利用经过步骤4处理后得到的及格和不及格样本建立逻辑回归模型。其中，为了保证评分公式的合理性和公平性，经过步骤1分析确定，除了施工进度等主观影响因素及部分不稳定的客观影响因素，将单环推进时间、单环停机时间、单环总时间、环数等掘进参数和地质类别共5个变量纳入模型中。
[0102]
5.1建立逻辑回归模型
[0103]
逻辑回归模型其本身是一个非线性回归模型，经过logit转换(连接函数)将相应变量y和线性自变量相联系，可以得到一个线性的形式，使用线性回归模型对参数进行估计，所以说逻辑回归模型是一个广义线性回归模型。
[0104]
经过步骤4得到及格样本为yi＝1；不及格样本为yi＝0。
[0105]
对于具有5个独立变量和地质类别的向量v＝(v1，v2，
…
，v5)，向量vi为掘进参数(单环推进时间、单环停机时间、单环总时间、环数)及地质类别。设条件概率p(y＝1|v)＝p为根据观测量相对于某事件v发生的概率。那么逻辑回归模型可以表示为：
[0106]
[0107]
这里称为逻辑函数。其中g(v)＝w0 w1v1
…
wnvn，那么在v条件下y不发生的概率为：
[0108][0109]
所以这个比值称为事件的发生比，简记为odds。对odds取对数得到：
[0110][0111]
5.2评分卡转换
[0112]
评分卡的背后逻辑是odds的变动与评分变动的映射(把odds映射为评分)。转换得到的评分公式为：score＝a-b*ln(odds)，其中a与b是常数。本发明采用的评分卡的评分范围是0～1000。
[0113]
通过以下两个式子计算得出a、b的值，得到最终的评分公式。
[0114]
1000＝a-b*ln(odds
第一名
)
[0115]
0＝a-b*ln(odds
最后一名
)
[0116]
5.3合理性检验
[0117]
本发明通过评分分布图和ks检验两个方面对模型进行合理性的检验。
[0118]
(1)评分分布图
[0119]
通过评分分布图可直观分析及格样本和不及格样本的分数分布情况。若两类样本的图象众数比较靠近，则说明该评分模型对两类样本的区分度高；若两者的图象众数之间存在较大的距离，则该模型的样本区分度不足。
[0120]
(2)ks检验
[0121]
本发明采用ks检验对该评分模型进行验证，ks检验是比较f
y＝1
(z)与f
y＝0
(z)两个样本分布的检验方法。
[0122]
若检验结果显著，则说明两类样本服从不同分布，该评分模型对这两类样本具有较高的区分度，可反映实际施工水平；若不显著，则需对评分模型的变量进行调整，重新选择更能准确反映掘进参数实际水平的变量及对掘进参数产生客观的影响作用的因素。
[0123]
通过该评分模型得到每环的实际施工分数，可反映不同施工队的管理水平，可将其作为工程队考核选拔的参考依据，进一步建立施工绩效评估系统；另外对施工水平的准确把握有利于工程进度的有效规划，降低项目成本。
[0124]
实施例：
[0125]
1.掘进速度影响因素分析
[0126]
在地铁建设工程中，对盾构掘进参数产生影响作用的外界因素有多种，包括地质环境、设备状态及政策变动等。地层类别的不同会影响盾构掘进的速度；设备状态是否正常会影响施工进程；政策变动会影响原定的工程规划和进度，导致不同情况下的掘进工时不一。
[0127]
采用秩检验的方法，将以上几个方面的影响因素与掘进数据，包括推进时间、停机时间等进行相关性检验，得到的结果如下表所示(以岩性的检验结果为例)，说明四个影响
因素与掘进数据之间均具有显著的相关性。
[0128]
表1显著性分析结果(岩性-推进时间)
[0129][0130]
表2显著性分析结果(岩性-停机时间)
[0131][0132]
2.掘进数据预处理
[0133]
2.1数据标准化处理
[0134]
为了准确界定掘进数据中的离群值所在范围，本发明采用z-score法对盾构掘进数据进行标准化处理。得到相对应的两组z数据，包括“z推进时间”、“z停机时间”。标准化后的掘进数据变为均值为0，方差为1的分布。
[0135]
2.2异常值的界定及剔除
[0136]
在采集到所有掘进数据中，存在部分异常值，包括过大值、过小值、零值及空值。本发明采用“三倍标准差法”，对超出-3σ至3σ范围的数据界定为异常值，即在两组z数据中，剔除数值大小超出-3至3范围的数据。经计算，最终被剔除的异常值所占原有总数据量的比例均在2.5％左右，说明异常值的界定范围选取恰当。
[0137]
2.3地质标注符号的统一化
[0138]
由于施工过程中的记录不当，收集到的地质类别数据中存在多种记录错误，导致同一类地层由于标注符号不统一而重复出现，被划分为多个类别，出现虚假的地质类别数。因此，需要将其记录形式统一化，减少重复的地质类别，最终确定所有区间的基本地质类别数为315种。
[0139]
3.地质聚类分析
[0140]
本方案进行聚类研究的数据是抽取相同盾构法对应区间的掘进参数。
[0141]
根据领域知识将全部地层划分为6类，标记为一、二、
…
、六类。
[0142]
按照上述聚类算法过程，对同一盾构法对应的每一类地层进行聚类分析，得到以
下聚类结果：
[0143]
表3最终聚类结果
[0144][0145]
由上表可知，在已有的6个地层分类维度的基础上分别对其地质类别进行逐一聚类，得到的最终聚类数共有14类。经分析，每一分类中各自包含的地质类别较为接近，且各个分类之间具有显著性差别，说明该聚类结果符合客观性和合理性等方面的要求。
[0146]
4.标准工时的测算
[0147]
4.1制定标准推进时间
[0148]
推进时间描述：将第2类地层的推进时间数据进行描述统计，得到如下结果：
[0149]
(1)描述统计
[0150]
表4推进时间描述
[0151][0152]
(2)百分位数表
[0153]
表5推进时间百分位数表
[0154][0155]
4.2统计推进时间分布
[0156]
统计推进时间数据得到的统计分布如表5所示。
[0157]
第2类地层的推进数据共有6128条，其推进时间的直方图如图3所示。
[0158]
4.3标准推进时间的确定
[0159]
通过分析工时的分布特征以及结合实际的盾构施工水平，最终选取推进时间的九十分位数作为标准推进工时，那么如上述分析可知，第2类地层的标准推进时间为49分钟/环。
[0160]
4.4定义及格样本和不及格样本
[0161]
确定第2类地层的标准工时后，则可划分出及格区间和不及格区间，两者比例为9∶1。由此可知在14类地层中，每类地层的及格样本y＝1所占比例为90％，不及格样本y＝0所占比例为10％。
[0162]
5.评价指标的制定
[0163]
5.1建立逻辑回归模型
[0164]
在现有地层聚类结果的基础上，分别对各类地层划分及格样本与不及格样本，其
中及格样本y＝1；不及格样本y＝0。将y作为因变量，单环推进时间、单环停机时间、单环总时间、地质类别及环数作为自变量，建立逻辑回归模型，结果如下表所示：
[0165]
表6逻辑回归输出表
[0166][0167][0168]
设推进时间为x1，停机时间为x2，总时间为x3，环数为x5[0169]
地质类别：
[0170]
由下列公式：
[0171]
ln(odds)＝g(x)＝w0 w1x1 w2x2 w3x3 w4x4 w5x5[0172]
可计算出每环odds的对数，方便接下来评分卡的转换，结果下表所示：
[0173]
表7评分表
[0174][0175]
5.2评分卡转换
[0176]
通过以下两个式子：
[0177]
1000＝a-b*ln(odds
第一名
)
[0178]
0＝a-b*ln(odds
最后一名
)
[0179]
计算得出a、b的值，其中a＝398.770，b＝-2.796，得到最终的评分公式为：
[0180]
score＝398.770-2.796
[0181]
*(6.903 0.001x
1-0.009x
2-0.044x3 340.607x
401
[0182]
83.597x
402
158.693x
403
254.954x
404
3.639x
405
[0183]
268.159x
406
0.750x
407
0.049x
408
0*x
409
237.798x
410
[0184]
237.757x
411
13.688x
412
12.351x
413
267.890x
414
[0185]
1.285x5)
[0186]
其中，x1为单环推进时间，x2为单环停机时间，x3为单环总时间，x
401
为第1类地层，x
402
为第2类地层，x
403
为第3类地层，x
404
为第4类地层，x
405
为第5类地层，x
406
为第6类地层，x
407
为第7类地层，x
408
为第8类地层，x
409
为第9类地层，x
410
为第10类地层，x
411
为第11类地层，x
412
为第12类地层，x
413
为第13类地层，x
414
为第14类地层，x5为单位时间内的环数。
[0187]
对现有评价对象评分后的结果如下表所示：
[0188]
表8评分结果
[0189][0190]
5.3模型的合理性检验
[0191]
(1)评分分布图
[0192]
不及格样本y＝0和及格样本y＝1所得分数的分布图如图4、图5所示。
[0193]
通过对及格样本和不及格样本分数分布图的分析，可知两者分布图山峰距离较远，说明本模型具有较好的区分度。
[0194]
(2)ks检验
[0195]
对不及格样本y＝0和及格样本y＝1所得分数进行ks检验分析，得到的结果如下表所示：
[0196]
表9ks检验结果
[0197][0198]
对ks检验结果进行分析后，可知显著性为0，不及格样本和及格样本服从不同的分布，具有显著性差异，说明本评分模型具有显著的区分度。
[0199]
附录：
[0200]
表10变量定义表
[0201][0202]