基于数字图像技术的土石混合体表征单元体尺度的确定方法技术

本发明专利技术公开了一种基于数字图像技术的土石混合体表征单元体尺度的确定方法，首先对获取到的彩色图像进行去噪和二值化，得到表示土体和石块的二值图像；然后根据原始图像尺寸确定尺寸等级，获取不同尺寸的二值图像；再提取每个图像中的图形边界并进行光滑处理，得到用于有限元网格划分的文件，将得到的文件转换为数值计算模型，进行数值试验，得到对应的等效宏观力学参数；最后，对试验结果进行整理，根据精度要求确定土石混合体的REV尺度。本发明专利技术克服了传统基于元胞自动模拟土石混合体方法的无法表现块石颗粒边界形状的缺点，和传统随机模型无法体现现场实际土石分布情况的缺点，得到的土石混合体的REV尺度充分考虑了现场土石混合体的特征，并且具有效率高、精度高、成本低等优点。

【技术实现步骤摘要】

本专利技术涉及一种确定土石混合体表征单元体尺度的方法，具体涉及一种基于数字图像技术的土石混合体表征单元体尺度的确定方法。
技术介绍
土石混合体是一种土体与碎石的混合物。天然形成的土石混合体多为残坡积物、崩坡积物和冲洪积物等，在三峡库区，川西即青藏高原，东南沿海福建、香港等地区广泛分布。在工程建设过程中的土石坝、矿山开采产生的尾矿坝、排土场、碎石土地基、路基等也存在大量的土石混合体。这些土石混合体是岩土工程领域经常碰到而又必须妥善处理的地质载体。土石混合体作为土体和碎石的混合体，是一种介于均质土体和碎裂岩体之间的特殊工程地质材料。土石混合体的宏观力学特性不仅取决于土、石本身的力学性质，也与土、石的相对比例，碎石形态等密切相关。因此合理确定土石混合体的等效力学参数，对岩土工程的设计、施工及变形预测和稳定性评价具有十分重要的意义。土石混合体等效参数具有明显的尺寸效应。即土石混合体的等效参数计算值或实验值随土石混合体模型尺度的变化而变化，但当尺度增大到某一临界值后，等效参数将趋近于常数，这一临界尺度就是该土石混合体的表征体元(Representativeelementvolume,REV)尺度。从理论上将，只有岩土体模型达到其REV尺度，其相关数值分析或实验所得到的参数才与宏观土石混合体参数具有可等效性，也才能反映实际土石混合体的力学性质，REV尺度是保持岩体的力学参数基本稳定时的岩体最小尺度。目前学者关于土石混合体的REV尺度的研究，多采用基于元胞自动模拟的土石混合体数值模型的方法，也有依据现场岩土体成分的统计结果，并构建随机数值模型的方法。数字图像技术(Digitalimageprocessing,DIP)是一种将数字图像转换为数字形式，并且通过各种数学算法提取出其重要信息的技术。这种技术已经被广泛地应用于工程、计算机科学、信息科学、统计学、物理学以及其他学科，并且还有很广阔的发展前景。岩土材料是一种典型的具有复杂细观结构的不均匀材料，应用数字图像技术，岩土材料内部的多相性或微观信息可以被大量地提取。在近二十年里，文献调查显示融合数字图像技术的计算方法对于这种不均匀材料的分析非常有效。
技术实现思路
专利技术目的：针对如何确定土石混合体表征体元尺度的问题，本专利技术提出一种基于数字图像技术的土石混合体表征体元尺度的确定方法。技术方案：本专利技术提出一种基于数字图像技术的土石混合体表征体元尺度的确定方法，包括如下步骤：步骤1：获取实际尺寸足够大的土石混合体的彩色数字图像图形；所述步骤1中土石混合体的彩色数字图像文件要求实际尺寸足够大(图形实际尺寸为长×宽＝Lx×Ly)。通过数码相机等获取的数字图像文件格式主要有JEPG、PNG和BMP。彩色图像在RGB色彩系统中有三个整数分别代表红、绿和蓝三种颜色。如果一个长×宽＝Lx×Ly的图像有M×N个像素点，那么x方向单位像素点所对应的实际长度为Sx＝Lx/N，同理y方向单位像素点所对应的实际长度为Sy＝Ly/M。步骤2：对获取的彩色数字图像进行去噪和二值化处理，得到表示土体和石块的二值图像，具体包括如下步骤：(2.1)对RGB彩色图像进行中值滤波去噪，具体为将R、G、B通道对应的图像分别进行中值滤波，再将处理后的结果融合；(2.2)将彩色图像从RGB空间转换到HIS空间，尽管RGB色彩系统在日常生活中应用非常普遍，但却不利于计算机程序进行图像识别，本专利技术采用便于计算机进行色彩辨识与操作的HIS色彩系统。由于数字图像在色相(H)与亮度(I)方面区分土体和石块不明显，而采用饱和度(S)可以明显区分二者，因此选择在饱和度(S)下的图像；(2.3)根据图像在HIS空间中S分量的频率分布情况设定二值化阈值，将图像二值化，具体为根据饱和度(S)值在HIS色彩空间中的频率分布直方图，设定其出现最大频率所对应饱和度(S)值作为数字图像二值化的阈值。(2.4)在土石区分不明显的部分(比如有时岩石表面会被土壤覆盖)，通过观察原始彩色图像中的岩土分布情况对二值图像中的土石边界进行手动优化，得到更符合实际情况的土石分布的二值图像，并利用画图软件辅助修改。步骤3：根据原始图像尺寸确定m个尺寸等级，并在每个尺寸等级下选取与原二值图像含石率相同且边长为Li的n个图像，得到m×n个不同尺寸的二值图像，具体包括如下步骤：(3.1)选定尺寸等级Li,1≤i≤m，其中0＜L1＜L2＜…＜Lm＜L，L＝min(Lx,Ly)，Lx是原始图像的长，Ly是原始图像的款，L1尽量小，Lm尽量接近L；(3.2)在每一个尺寸等级Li下，选取与原二值图像含石率相同且边长为Li的n个正方形图形。步骤4：对步骤3中得到的所有二值图像提取图形边界，并将边界图像转化为矢量图形，具体包括如下步骤：(4.1)利用4连通域标记方法提取二值图像中表示岩石或土体其中一种的像素点集合；(4.2)对属于该集合的每一个点根据位置拆分为四条线，记录四条线的坐标；(4.3)对所有线的集合进行运算，删除重合的线，得到该像素点集合的外边框；(4.4)将外边框进行排序处理，形成首尾相连的锯齿状边界；(4.5)对锯齿状边界进行平滑处理；(4.6)根据一个像素点的实际尺寸，将边界图像文件转换为矢量图形文件。步骤5：将步骤4中得到的矢量图形划分有限元网格，并将网格图像转换为数值计算模型进行数值试验，得到对应的等效宏观力学参数，具体包括如下步骤：(5.1)利用有限元网格划分软件(如Ansys，HyperWorks和Gmsh等)进行网格划分；(5.2)针对每一个图形，分别通过数值方法模拟岩土力学实验，如单轴压缩试验、三轴压缩试验和直剪试验等(其中三轴压缩试验在二维的情况下描述为双轴压缩试验)，得到每一个图形所对应的的等效宏观力学参数，如等效粘结力、等效内摩擦角、等效切线弹性模量、等效平均弹性模量等。(6)对步骤5中得到的数值试验结果进行整理，根据精度要求确定土石混合体的REV尺度，具体包括如下步骤：(6.1)计算每个尺寸等级下等效宏观力学参数的均值和变异系数(CV)值，以等效弹性模量作为待求等效宏观参数为例，尺寸等级Li下均值的计算公式为：Ei=Σj=1nEijn,(i=1,2,...,m)---(1)]]>尺寸等级Li下变异系数(CV)值的计算公式为：CVi=SDiEi×100%,(i=1,2,...,m)---(2)]]>其中，本文档来自技高网...

【技术保护点】
一种基于数字图像技术的土石混合体表征单元体尺度的确定方法，其特征在于包括如下步骤：(1)获取实际尺寸足够大的土石混合体的彩色数字图像；(2)对获取的彩色数字图像进行去噪和二值化处理，得到表示土体和石块的二值图像；(3)根据原始图像尺寸确定m个尺寸等级，并在每个尺寸等级下选取与原二值图像含石率相同且边长为Li的n个图像，得到m×n个不同尺寸的二值图像；(4)对步骤3中得到的所有二值图像提取图形边界，并将边界图像转化为矢量图形；(5)将步骤4中得到的矢量图形划分有限元网格，并将网格图像转换为数值计算模型进行数值试验，得到对应的等效宏观力学参数；(6)对步骤5中得到的数值试验结果进行整理，根据精度要求确定土石混合体的REV尺度。

【技术特征摘要】
1.一种基于数字图像技术的土石混合体表征单元体尺度的确定方法，其特
征在于包括如下步骤：
(1)获取实际尺寸足够大的土石混合体的彩色数字图像；
(2)对获取的彩色数字图像进行去噪和二值化处理，得到表示土体和石块
的二值图像；
(3)根据原始图像尺寸确定m个尺寸等级，并在每个尺寸等级下选取与原
二值图像含石率相同且边长为Li的n个图像，得到m×n个不同尺寸的二值图像；
(4)对步骤3中得到的所有二值图像提取图形边界，并将边界图像转化为
矢量图形；
(5)将步骤4中得到的矢量图形划分有限元网格，并将网格图像转换为数
值计算模型进行数值试验，得到对应的等效宏观力学参数；
(6)对步骤5中得到的数值试验结果进行整理，根据精度要求确定土石混
合体的REV尺度。
2.根据权利要求1所述的基于数字图像技术的土石混合体表征单元体尺度
的确定方法，其特征在于，所述步骤2具体包括如下步骤：
(2.1)对RGB彩色图像进行中值滤波去噪；
(2.2)将彩色图像从RGB空间转换到HIS空间；
(2.3)根据图像在HIS空间中S分量的频率分布情况设定二值化阈值，将
图像二值化；
(2.4)在土石区分不明显的部分，通过观察原始彩色图像中的岩土分布情
况对二值图像中的土石边界进行手动优化，得到更符合实际情况的土石分布的二
值图像，并利用画图软件辅助修改。
3.根据权利要求2所述的基于数字图像技术的土石混合体表征单元体尺度
的确定方法，其特征在于，所述步骤2.3中确定图像二值化阈值的方法具体为：
根据饱和度值在HIS色彩空间中的频率分布直方图，设定其出现最大频率所对
应饱和度值作为图像二值化的阈值。
4.根据权利要求1所述的基于数字图像技术的土石混合体表征单元体尺度
的确定方法，其特征在于，所述步骤3具体包括如下步骤：
(3.1)选定尺寸等级Li,1≤i≤m，其中0＜L1＜L2＜…＜Lm＜L，

\tL＝min(Lx,Ly)，Lx是原始图像的长，Ly是原始图像的款，L1尽量小，Lm尽量接
近L；
(3.2)在每一个尺寸等级Li下，选取与原二值图像含石率相同且边长为Li的
n个正方形图形。
5.根据权利要求1所述的基于数字图像技术...

【专利技术属性】
技术研发人员：杨兰兰，徐卫亚，孟庆祥，闫龙，张强，张久长，俞隽，王如宾，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32