基于数字图像处理的土石坝坝料合格性智能检测方法

赵宇飞，刘 彪，王 毅，孟 亮，刘必旺

(1.中国水利水电科学研究院 流域水循环模拟与调控国家重点实验室，北京 100038；2.中国水利水电第八工程局有限公司，湖南 长沙 410007；3.中国水利水电第六工程局有限公司，辽宁 沈阳 110179)

1 研究背景

按照常规的土石坝工程建设施工管理模式，土石坝坝料合格性检测通常是在大坝单元工程填筑碾压结束后通过挖坑获得碾压之后的坝料，然后进行筛分试验得到级配曲线，通过比较级配曲线是否在设计单位提供的上下包络线范围之内，进而判断坝料是否满足设计要求[1]。这样的检测和合格性评价方法存在采样率低，操作过程繁琐，检测过程受人为因素影响大等缺点以致检测结果代表性差[2-4]，更为严重的是这样的检验属于典型的事后检验，一旦坝料级配不合格，需要将碾压完成的坝料挖除，然后换合格坝料重新进行填筑施工，浪费大量的人力物力。因此，在实际土石坝施工中如何快速获取坝料的级配特性，并与设计级配上下包络线相关级配特性指标进行比对，进而判断坝料是否合格，对于大坝填筑施工来说具有重要的意义，也是目前土石坝填筑施工精细化智能控制的重要研究方面。

随着计算机和人工智能技术的发展，利用数字图像处理技术实现颗粒级配的智能检测已经在选矿、道路填料生产以及混凝土骨料检验等工程领域中得到了广泛的应用。发展至今，逐渐形成了基于颗粒三维形态捕捉的级配检测方法和基于二维图像形态提取转换的级配检测方法。其中基于颗粒三维形态捕捉的级配检测方法常采用三维激光扫描技术量化颗粒的真实形态[5-7]，算法具有较高的鲁棒性，已有研究表明该方法比二维形态量化方法更加的准确有效[8]，但是该方法需要精密的设备，成本较高，操作复杂，检测效率低，多适用于试验室环境检测，不适用于环境复杂的现场施工检测[9]。相比而言基于二维图像形态提取转换的级配检测方法检测效率高，操作简单，该方法由颗粒二维图像轮廓提取和基于二维形态特征的级配转换两部分组成[10]，适合于施工现场级配检测，受到了众多学者的研究。

一些学者提出采用基于经典图像分割技术的二维图像形态级配检测方法，如张国英[11]、于沭[12]等通过采用图像分割技术实现了颗粒的二维轮廓提取，而后采用椭球拟合获取颗粒的三维体积模型，建立了颗粒的粒度检测与分析系统；蔡改贫等[13]采用基于分水岭分割算法对在线采集到的矿石颗粒进行形态提取，然后采用击中与击不中算法统计不同粒径的矿石颗粒数量，实现了矿石粒度的在线实时检测；宫兴等[14]采用基于阈值化分割的方法实现了对沥青混合料的平面形态提取，然后提出了将沥青混合料的平面级配作为输入，设计级配作为模型的输出，建立了BP神经网络来对沥青混合料的设计级配进行检测，为级配检测提供了新的思路。以上基于经典图像分割技术的级配检测方法虽然识别速度快，但是识别结果高度依赖于关键参数的手动选择，对参数的选取较为敏感，若参数选取不当会导致欠分割或过分割的现象，此外针对复杂成像条件以及粘连颗粒图像的分割精度仍需进一步的提高[15]。近年来，随着深度学习的异军突起，基于深度卷积神经网络(Deep Convolution Neural Network，DCNN)的图像分割方法在计算机视觉和医学图像分析领域占据主导地位，也为土石坝料二维图像轮廓提取带来了新的机遇。该方法通过建立深度学习神经网络模型可以自动提取隐含的抽象特征，获取泛化能力强的形态特征，从而自主解决经典图像处理方法的限制[16]，与经典的图像分割算法相比，分割精度有了极大的提升。不少学者也逐渐提出了很多复杂的以及不同方法组合的算法来提高分割的精度[16-19]，虽然分割精度高，但是该些方法的效率较低，对计算机的硬件性能要求较高，分割的准确率高度依赖于模型结构设计和样本标记，而土石料图像中存在的颗粒数量较多，人工标记颗粒轮廓的过程费时费力，模型训练耗时时间长。相比深度学习的方法而言，目前实际工程中常采用基于经典图像分割技术的二维图像形态级配检测方法。

综上所述虽然图像识别技术在诸多工程领域得到了广泛的应用，但是该技术在土石坝坝料级配智能检测方面却鲜有相关研究，主要原因如下：(1)坝料颗粒粒径大小差异大，形状也各不相同；(2)所拍摄的图像含有杂质，会有重叠、噪声和伪影，且图像往往受到光照不一致的影响，导致图像的亮度不均衡；(3)图像在拍摄过程中需要固定焦距或者一定的比例尺，没有明显的边界，将使颗粒尺寸的分析变得异常困难。这些外界的干扰降低了目标与背景的对比度，严重影响了土石坝坝料图像的有效分析。因此，坝料级配智能检测既是一项富有挑战性的任务，又是目前亟待解决的工程问题。考虑到实际施工现场的环境条件以及检测效率，本文采用基于二维图像形态提取转换的级配检测方法来实现土石坝坝料的级配智能检测，本文方法整体的逻辑框架如图1所示。首先提出了采用融合空间信息的直觉模糊C均值(SIFCM)算法实现了土石坝坝料的图像分割，有效改善了标准模糊C均值算法对噪声敏感的缺陷，然后采用了椭球体积来重构土石坝料的三维体积形态，最后建立了BP神经网络模型修正图像识别获取的级配曲线以获得真实条件下坝料的级配特征参数，从而获得4个坝料合格性评价指标：最大粒径、P5含量、曲率系数Cc和不均匀系数Cu，为坝料合格性快速判别和建立坝料碾压施工过程实时压实质量多层次综合评价模型提供了重要的输入信息。

图1 本文方法逻辑框架

2 二维图像形态提取

2.1 图像获取的条件在水利水电工程中，大坝填筑料的粒径组成是影响工程施工质量的重要因素，工程设计中对坝料的级配有明确的要求。采用图像识别技术对常规的土石坝坝料级配分析时，需要对坝料的数字图像进行图像分割获取颗粒轮廓，进而确定坝料级配特征。因此数字图像的获取对级配特征的识别分析至关重要。

本研究中，在获取数字图像时，需要在施工现场摊铺结束后的坝料上放置1 m×1 m的铁框作为比例标尺框，拍摄图像时要有较为充足的光照条件，保证拍摄图像的位置处于比例标尺框正上方，尽量避免比例标尺框内存在阴影，尽可能的提高数字图像的整体质量，最大程度上为图像识别减少外界干扰的影响。如图2所示，现场拍摄的图像经过编制的程序处理可实现比例标尺框内图像的快速提取。通过将比例标尺框真实尺寸和图片像素对比，可建立真实尺寸与图像像素间的转换关系，实现图像中坝料颗粒像素与实际尺寸之间的换算。此外碾压结束之后进行碾压质量坑检时，尽量选择在坝料数字图像获取的位置处进行，这样有利于进行图像识别得到的级配曲线与现场挖坑试验获得的级配曲线(以下简称为试验级配曲线)间的对比研究，并且可以将对比结果应用于其他级配检测结果修正，提高级配检测的准确度。

图2 提取标志物框内图像

2.2 图像预处理在实际施工现场通常是通过相机来获取颗粒图像，成像过程会受到外界复杂环境的影响，如光照不均匀、拍摄器件不稳定以及其他颗粒的影响等，导致图像中存在噪声信息，严重影响识别的精度。为了提高土石坝坝料数字图像分析的精度和稳定性，需要对采集到的坝料原始图像进行一定的图像预处理来消除噪声、增强目标信息，主要包括空间滤波和频率域滤波。本文在空间滤波方法中选用中值滤波消除椒盐噪声，双边滤波消除高斯噪声。在频率域滤波方法中选用同态滤波(Homomorphic Filtering，HF)方法[20-23]修正图像中存在的亮度不均匀、目标物体阴影以及边缘细节信息雾化等问题。

如图3(a)是经过空间滤波处理后的土石坝坝料数字图像，从图中可以看出图片由于光照不均匀等外界干扰影响而存在部分区域有伪阴影和细节特征模糊的现象，直接对其进行二值化处理后，图3(d)中区域A和C内的信息丢失，严重影响图像识别的精度。如果想要修正图像中偏暗的区域，必须设计合适的同态滤波器使得暗区的亮度饱和，提高整体的对比度。本文通过调整滤波器参数，最终设计了融合两种同态滤波器HF1和HF2的方法来修正图像降质的问题，两种滤波器的参数见表1。经过HF1滤波后虽然区域A和C丢失的信息未补全，但是图像的整体亮度有所提高，中大型颗粒轮廓分明，类似B部分的前景区域边缘细节特征得到加强，更有利于图像轮廓的提取，减少粘连颗粒现象的发生；经过HF2滤波后，图像整体的对比度有了很大幅度的提高，区域A和C丢失的前景信息得到修正补全，且小颗粒坝料轮廓分明，但是由于亮度提高过大，会存在临近的大颗粒坝料出现粘连的现象，使其边缘信息模糊化。因此本文设计了两种滤波器，采用HF1滤波器识别中大型颗粒和边缘细节特征信息，HF2滤波器修正补全图像中存在的伪阴影区域信息并用于识别零散的小颗粒物体。综合使用两种滤波器处理后的图像亮度和对比度均得到了有效的提升，图像质量得到改善，且目标的边缘没有受到损失，效果较佳。

图3 同态滤波器设计

表1 同态滤波器参数

2.3 融合空间信息的直觉模糊C均值聚类算法的坝料图像二值化在土石坝坝料图像处理过程中，我们只对图像中坝料颗粒所处的区域和形态感兴趣，这部分我们希望获取的目标区域常称为前景区域，而其他无关的区域称为背景区域。为了将前景和背景区域分开，需要将图像预处理后的灰度图像进行二值化处理从而生成二值图。在坝料图像中，通常我们所关注的坝料颗粒的亮度要高于颗粒间孔隙的亮度，且孔隙在灰度图像中偏暗色，因此采用基于灰度准则的图像分割方法可以较好的实现颗粒形态的提取。

标准模糊C均值(FCM)算法由于本身只考虑了图像灰度的相似性，而没有考虑图像像素邻域信息，导致图像分割的结果容易受到噪声以及伪阴影的影响，不能得到很好的分割效果[24]。为了改善FCM算法分割存在的缺陷，本文提出一种融合空间信息的直觉模糊C均值聚类(SIFCM)算法进行土石坝料的图像分割。该算法首先按照标准FCM算法的逻辑迭代更新隶属度矩阵和聚类中心，然后考虑空间信息对聚类结果的影响，将空间信息以空间函数的形式纳入隶属度矩阵中，其中空间函数为所考虑的像素邻域内所有隶属函数的总和，最后引入直觉模糊集中的直觉指数和非隶属度函数，更好地反映聚类的模糊问题。

2.3.1 基于FCM算法的图像二值化 FCM算法最早由Dunn提出，并由Bezdek将之推广[25]，被广泛地应用于模式识别、数据挖掘、图像分割以及计算机视觉等众多领域。该方法是一种基于划分的聚类算法，允许一个数据属于两个或多个类别，其基本思想是使被划分到同一簇对象之间相似度最大，而不同簇之间的相似度最小。算法定义如下：假设X={x1，x2，…，xN}为待聚类数字图像所有像素点灰度值的数据集，FCM算法将每个对象xj(1≤j≤N)分配给C个类别，使用C×N的隶属矩阵U={uij}表示对象与类别之间的隶属关系，其中uij表示第j个对象属于第i个类的隶属程度。针对X的模糊划分空间定义为[26]：

(1)

FCM算法的目标是计算给定数据集X的隶属度矩阵和聚类中心V={v1，v2，…，vC}，目标函数将欧氏距离作为需要最小化的代价函数，表示为：

(2)

式中：m为控制模糊程度的加权指数，当m趋于1时，分区接近硬聚类，m趋于无穷大时，分区将产生分离得很差的集群，建议值在1.5和2.5之间[27]，本研究使用默认值2；|xj-vi|为第j个像素点灰度值与第i个聚类中心之间的欧几里得距离。模糊划分是通过对目标函数的迭代优化来实现的，当目标函数对uij和vi的一阶导数为零时，Jm达到最小化，保证了组内相似度最高，组间相似度最低的聚类原则，此时隶属矩阵和聚类中心的更新方程为[28]：

(3)

(4)

算法迭代优化的终止条件为：

(5)

式中：α为一个很小的常数表示误差阈值，本文取0.05；t为迭代次数，迭代优化过程最终收敛于Jm的局部极小值点或鞍点[29]。

2.3.2 空间函数 图像的重要特征之一是相邻像素具有高度相关性，即这些相邻像素具有相似的特征值，属于同一个簇的概率很大。这种空间关系在聚类中很重要，但在标准FCM算法中没有充分利用这种空间信息。为了充分考虑空间信息的影响，将其定义为空间函数[24]：

(6)

式中：NB(xj)为空间域中以像素xj为中心的正方形窗口，本文采用5×5的像素领域；空间函数hij与隶属度函数一样，为像素xj属于簇i的概率。如果一个像素周围的大部分邻域属于同一个类别，那么它的空间函数就很大。然后将空间函数纳入隶属度函数，如下所示：

(7)

式中p和q为控制两个函数相对重要性的参数。在同质区域，由于没有噪声的影响，空间函数对原始的隶属度仅起到强化的作用，不会对聚类结果造成改变。但是若有噪声像素的存在，式(7)综合考虑了噪声点邻域内的像素来降低噪声对聚类的影响权重。考虑空间信息后，改进算法可以实现对噪声区域或伪阴影中容易错误分类的像素纠错处理[30]。参数p和q的空间FCM表示为SFCMp，q，当参数p=1，q=0时SFCM1，0与传统的FCM算法相同。针对本研究中土石坝坝料数字图像分割而言，经过多次不同参数的结果相比，最终确定p=1，q=3的空间FCM算法分割效果最佳。

2.3.3 直觉模糊集理论 Atanassov[31]提出的直觉模糊集(Intuitionistic Fuzzy Set，IFS)理论考虑了数字图像所有像素点灰度值集合的隶属度u和非隶属度v，设论域X上的直觉模糊集A可以用数学方法表示为：

A={(x，uA(x)，vA(x)，πA(x))|x∈X}

(8)

式中：uA(x)∈[0，1]为x属于A的隶属度；vA(x)∈[0，1]为x属于A的非隶属度；πA(x)∈[0，1]为直觉指数，表征由于缺乏知识和信息引起的犹豫程度；对于x∈X，均满足uA(x)+vA(x)+πA(x)=1。在直觉模糊集构造时，通常采用Sugeno[32]生成函数来生成非隶属度函数，表达式如下：

(9)

式中λ为模糊补算子。则对应的直觉指数为：

(10)

将式(10)加入到模糊隶属度矩阵中，得到

(11)

基于FCM算法的土石坝料图像分割在无监督的情况下，可以自适应迭代获得最终图像的二值化图，结果如图4(a)所示，采用本文提出的基于SIFCM算法的土石坝料图像二值化结果如图4(b)所示。与标准FCM算法处理得到的坝料图像二值化结果相比，由于本文算法综合考虑了图像的空间信息以及直觉模糊集理论，算法对噪声的敏感程度显著降低，提高了聚类的效果和精度，且尽可能的保留了图像的细节信息，性能要优于标准FCM算法，通过对石块的识别结果对比验证了本文算法的准确性和优越性。

图4 二值化图对比

2.4 形态学操作由于土石坝料形状随机、纹理多变、成像过程中存在复杂的施工环境以及设备不稳定等诸多外界因素的干扰，不可避免的会在图像前景和背景区域产生噪声，还会存在颗粒粘连的现象。为了对含噪声以及粘连颗粒的二值化图进行优化，本文设计了如图5所示的形态学滤波器。首先进行孔洞填充以及一次开闭运算初步消除前景和背景区域中的噪声，然后选用不同结构元素进行多轮形态学操作并做合并处理，最后综合考虑两种同态滤波器的工作特性，分别设置不同的面积阈值提取坝料颗粒轮廓。由于HF1滤波器处理后的图像对于大颗粒坝料识别效果较好，对于小粒径坝料以及伪阴影部分颗粒会存在部分信息丢失的现象，因此在提取经HF1同态滤波器处理后的坝料二值化图时设置最小面积阈值Amin作为下限，此处Amin取所有粒径最佳椭圆拟合面积的均值减去一个标准差，提取大于Amin的所有颗粒(图6(a))，可以看出上述方法可以有效分割图像中粘连颗粒(类似区域B)的部分；HF2滤波器提高了图像整体的对比度和亮度，不仅修正补全了伪阴影面积而且对图像中小颗粒坝料前景区域的识别有了整体提高，不过也会造成较大颗粒的粘连，因此在提取二值化图时设置最大面积阈值Amax作为上限，此处Amax取所有粒径最佳椭圆拟合面积的均值加上一个标准差，提取小于Amax的所有颗粒(图6(b))，可以看出该方法可以有效识别小颗粒和伪阴影区域的颗粒。最后对基于HF1和HF2处理后的图像中提取到的颗粒进行融合得到整体坝料二值图(图6(c))，可以看出经过本文设计的形态学滤波器处理后坝料图像中颗粒轮廓识别效果较好，有效解决了伪阴影和粘连颗粒识别的难题。

图5 形态学滤波器

图6 形态学处理后的二值化图

2.5 边缘检测边缘检测的方法是利用前景与背景区域内灰度高度相似，而边缘处灰度不连续的理论将图像边界提取出来。Canny边缘检测是目前非常流行的算法之一，主要由图像降噪、梯度求解、非极大值抑制和阈值筛选四阶段构成[33]。采用Canny边缘检测对经过HF1滤波器和HF2滤波器处理后的土石坝料图像分别进行Canny边缘检测，最后将识别出的轮廓取并集得到检测结果见图7所示。本章采用基于SIFCM算法的图像分割方法实现了土石坝料的图像分割，最终的分割结果如图8所示，可以看出分割精度较高，也为准确进行坝料形态特征级配转换从而获得坝料级配曲线提供了较高的技术支撑。

图7 边缘轮廓提取

图8 最终图像分割结果

3 土石坝料级配转换与修正模型

3.1 三维体积重构的级配转换从采集到的图像中可以看出土石坝料颗粒二维形态整体呈现出类椭圆状，在以往的土石坝料、沥青集料以及采矿工程中颗粒块体的研究中，也常用椭圆形态拟合颗粒的真实形态。由于数字图像所反映的只是土石坝料摊铺后表面所出露的情况，而坝料的级配曲线是依据各粒径范围内的土石坝料颗粒含量的质量占比计算得到的，涉及到颗粒的三维形态特征，因此本文采用椭球体积来重构土石坝料的三维形态特征。

Kemeny等[34-36]在矿石粒径检测研究中发现，通过筛网的颗粒等效粒径可以表示为：

(12)

(13)

(14)

式中：a和b分别为颗粒轮廓最佳拟合外接椭圆的长轴和短轴；C和A分别为最佳拟合外接椭圆的周长和面积。

张国英等[11]在此基础上建立了矿石颗粒的三维等效椭球体积公式，即：

(15)

假设本研究获取到的坝料颗粒的密度均相同，则各粒径范围内土石坝料颗粒的级配质量百分比GQP为：

(16)

式中：ρ为坝料颗粒的密度；V为坝料颗粒的体积；A为最佳拟合外接椭圆面积。

3.2 基于BP神经网络的级配修正模型本文所采用的基于二维图像形态提取转换的级配检测方法中，二维图像的分割结果不可避免的会与原始图像间存在一定的差异，此外土石坝料的三维体积重构也是基于一定的假设基础上建立的，因此通过数字图像识别出的级配曲线与施工现场挖坑检测得到的试验级配曲线之间存在一定的误差。为了修正两者之间的误差，提高坝料级配智能检测结果的准确性，本文通过对碾压结束后挖坑检测位置处拍摄的坝料图像识别得到的级配曲线与实际挖坑筛分得到的级配曲线进行对比分析，建立了基于BP神经网络的级配修正模型。

本文所获取的坝料图像和挖坑检测级配数据来源于辽宁清原抽水蓄能电站下水库主、次堆区，共获取到24组挖坑检测的试验级配曲线(图9)和对应位置处的坝料图像。由于本文选取的挖坑检测资料和图像均来自主次堆石区，图像中最大粒径可识别到200 mm，图像中以较大的石块居多，因此在土石料级配曲线求解时，选取的最小过筛尺寸为5 mm，取d=5、10、20、40、60、80、100和200 mm作为过筛粒径，然后采用本文提出的基于二维图像形态提取转换的级配检测方法可求得挖坑检测位置处过筛粒径对应的级配质量百分比GQP′，将过筛粒径d和图像识别求解的GQP′作为模型的输入，现场挖坑检测获得的GQP作为模型的输出，采用BP神经网络建立图像识别级配曲线与试验级配曲线之间的修正模型。

图9 现场挖坑检测试验级配曲线

图10 预测值和实测值的对比

级配修正模型采用4层BP神经网络，输入层有2个节点，分别代表着d和GQP′输入参数，两个隐藏层的节点数均为10，输出层有1个节点，学习率η取为0.001，期望误差goal取0.02，最大训练次数epoch取为200次。为了消除输入参数中不同数据量纲之间的差异对预测模型的影响，对数据进行归一化处理。将24×8组数据中的90%作为训练集，剩余10%数据作为测试集，利用训练集对级配修正模型进行训练，从而得到最优超参数；再将测试集中d和GQP′当作自变量输入到训练好的模型，输出的因变量即为预测GQP。图10为级配修正模型测试组预测值和真实值的比较，两者的决定系数R2为0.956，表明采用人工神经网络方法建立的级配修正模型具有的精度满足实际工程需要。

为了进一步验证模型的修正效果，将现场挖坑检测得到的试验级配曲线作为标准，取测试集中一张图片进行级配智能检测，坝料图像的分割结果如图11所示，获得坝料二维形态特征后进行颗粒三维体积重构求得图像的级配曲线，然后将其代入训练好的BP级配修正模型，修正后的级配曲线与试验级配曲线对比如图12(a)所示。本文图像识别修正级配与试验级配间误差如表2所示，相对误差小于5%，说明本文建立的级配修正模型具有较好的可行性。

图11 测试集中某原始图像和分割结果

表2 图像识别修正级配与试验级配间误差

上述通过本文提出的SIFCM_BP算法的级配智能识别修正模型可以在大坝坝料进场之前快速识别坝料5～200 mm的级配曲线，为了进一步补全级配曲线中丢失的部分，本文将图像识别得到的曲线与施工现场挖坑检测获得的试验级配数据库中的曲线依次进行比对，以两条曲线各粒径绝对误差的累计和作为目标函数f(x)，可以表示为：

(17)

式中：d1=5 mm，d2=10 mm，…，dn=200 mm；GQP″i为粒径di对应的图像识别修正后获得的级配质量百分比；GQPi为粒径di对应的现场试验获得的级配质量百分比。选取目标函数f(x)达到最小值时对应的试验级配曲线来补全本文图像识别获得的级配曲线中丢失的粒径小于5 mm以及大于200 mm的部分级配数据，补全后的级配曲线如图12(b)所示。通过对级配曲线进行特征分析，可以得到用于坝料合格性检测的4个指标：最大粒径、P5含量、曲率系数Cc和不均匀系数Cu。

图12 图像识别获取的级配曲线

3.3 相关问题讨论本文建立了基于SIFCM_BP算法的级配智能识别修正模型，实现了坝料级配特征参数的快速可靠获取，并在此基础上完成了坝料的合格性检测。虽然取得了一定的成果，但是本文提出的基于数字图像的坝料合格性快速检测方法仍处于初步研究阶段。由于本文选取的挖坑检测资料和图像均来自主次堆石区，图像中以较大的石块居多，因此在土石料级配曲线求解时，选取的最小过筛尺寸为5 mm，通过本文方法可以修正补全图像识别获得的级配曲线中丢失的粒径小于5 mm以及大于200 mm的部分级配数据。但是对于粉土、黏土来说，其粒径尺寸太小，难以利用常规堆石坝主堆区爆破料、砂砾石料数字图像获取技术进行级配特征识别的图像分析。本文在后续的研究过程中，将在已有成果基础上，尝试利用电镜扫描等技术开展粉土、黏土颗粒级配特征获取的数字图像分析方法研究。

此外由于实际工程中，一个单元工程碾压结束后只进行一个挖坑质量检测，因此本文收集到的现场挖坑检测获得的试验级配数据和对应位置处的数字图像数据量较少，修正模型的预测精度还有待进一步的提升。随着单元工程不断填筑碾压完成，图像识别获得的级配数据以及挖坑检测获得的试验级配数据会不断积累，用于级配修正模型训练的数据库会变得逐渐丰富和完善，这样通过本文方法获得的图像识别结果更加接近于真实条件下的坝料全级配特征曲线，对于本文图像识别未能得到的坝料缺失部分级配的补全精度也会逐步的提高，因此本文方法具有不断进化学习的工程优势。本文图像识别方法是基于经典图像分割技术基础上改进的二维图像形态级配检测方法，虽然识别速度快，识别精度有所提高，但是一些关键参数的选择仍然依赖于人工选取，若想进一步提升识别精度，可尝试采用基于深度卷积神经网络的图像分割方法进行工程实例应用，并在人工标记颗粒轮廓方面寻找新的快速标记方法，以解决该方法耗时长的弊端，更有利于在实际工程中的应用。

另外，需要说明的是，在实际工程中，堆石坝坝料的开采是需要进行爆破试验完成的，试验前会综合考虑岩体强度和结构面发育程度等因素对爆破参数(如炮孔间距、装药量、引爆延时等)进行调整，保证开采得到的堆石料能够最大程度上满足坝料级配的要求。然而本文坝料合格性检测是在坝料均匀摊铺后通过对采集到的数字图像进行智能识别得到的级配参数，因此，暂未考虑岩石分类及岩性的影响。在后续的研究中，将收集堆石料料场结构面发育及岩石强度特性的信息，对料场开采区域的岩石分类及岩性进行确定，然后在级配的学习进化方面考虑堆石料颗粒特性(形状特性)参数与级配特征参数的自适应与自生长的优化能力。

4 结论

在土石坝填筑工程中，坝料的级配特性一般都是在碾压完成之后进行挖坑检测得到的，级配不合格的坝料难以在出料场以及摊铺碾压过程中及时清理出场，造成施工机械和人工的浪费。针对现有级配智能感知技术较差的问题，本文提出了一种SIFCM_BP算法的级配智能识别修正模型实现了坝料级配特征参数的快速可靠获取，为合理有效地检测坝料的合格性提供了新的技术手段。主要结论如下：

(1)采用了空间滤波和同态滤波相融合的技术为土石坝坝料的快速准确分割提供了技术支撑；经设计的HF1同态滤波器处理后的图像对于大颗粒坝料以及粘连颗粒的识别效果较好；HF2同态滤波器处理后的图像不仅修正补全了伪阴影面积而且对图像中小颗粒坝料前景区域的识别有了整体提高，综合使用两种同态滤波器可以提高图像识别的整体精度。

(2)针对FCM算法进行图像分割时结果易受到噪声以及伪阴影影响的缺陷，本文在标准FCM算法的基础上考虑了空间信息对聚类结果的影响，将空间信息以空间函数的形式纳入隶属度矩阵中，然后引入直觉模糊集中的直觉指数和非隶属度函数，更好的反映聚类的模糊问题，较大幅度的提高了算法对噪声的鲁棒性。与标准的FCM算法相比，图像的分割精度有了明显的提升。

(3)在获得土石坝料图像的二维形态特征后，通过建立等效椭球体积的方法实现了土石坝料的三维体积重构，获得了基于图像识别技术的坝料级配曲线，然后建立了基于BP神经网络的级配修正模型，提高了坝料级配智能检测结果的准确性。对本文SIFCM_BP模型获得的级配曲线进行特征分析，最终获得评价坝料合格性4个指标：最大粒径、P5含量、曲率系数Cc和不均匀系数Cu，为坝料合格性快速判别和本文后续建立坝料压实质量多层次综合评价模型提供了基础信息，有利于施工前的坝料质量控制与施工过程的质量评价。