• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

基于计算机视觉的岩石裂隙识别表征与软件研制

李元海, 徐晓华, 朱鸿鹄, 杨硕, 唐晓杰, 赵万勇

李元海, 徐晓华, 朱鸿鹄, 杨硕, 唐晓杰, 赵万勇. 基于计算机视觉的岩石裂隙识别表征与软件研制[J]. 岩土工程学报, 2024, 46(3): 459-469. DOI: 10.11779/CJGE20221239
引用本文: 李元海, 徐晓华, 朱鸿鹄, 杨硕, 唐晓杰, 赵万勇. 基于计算机视觉的岩石裂隙识别表征与软件研制[J]. 岩土工程学报, 2024, 46(3): 459-469. DOI: 10.11779/CJGE20221239
LI Yuanhai, XU Xiaohua, ZHU Honghu, YANG Shuo, TANG Xiaojie, ZHAO Wanyong. Identification and characterization of rock fractures based on computer vision and software development[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(3): 459-469. DOI: 10.11779/CJGE20221239
Citation: LI Yuanhai, XU Xiaohua, ZHU Honghu, YANG Shuo, TANG Xiaojie, ZHAO Wanyong. Identification and characterization of rock fractures based on computer vision and software development[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(3): 459-469. DOI: 10.11779/CJGE20221239

基于计算机视觉的岩石裂隙识别表征与软件研制  English Version

基金项目: 

国家自然科学基金项目 52274141

国家重点研发项目 2022YFC3003304

徐州市重点研发项目 KC21310

详细信息
    作者简介:

    李元海(1969—),男,博士,教授,博士生导师,主要从事岩土工程计算机视觉量测和隧道工程围岩稳定性方面的研究。E-mail: Lyh@cumt.edu.cn

    通讯作者:

    徐晓华, E-mail: 1750904260@qq.com

  • 中图分类号: TU45

Identification and characterization of rock fractures based on computer vision and software development

  • 摘要: 岩石裂隙特征是评判岩体结构及其完整性的核心指标,也是评估岩石工程安全稳定性的重要因素。针对岩石裂隙识别,采用深度学习方法,通过引入混合注意力机制对Unet模型进行了改进,有效提高了岩石裂隙识别的精度。针对交叉岩石裂隙的分离与特征提取,提出了一种基于迹线方向判定的裂隙分离与表征算法,依据裂隙分离的结果形式,采用重合追踪法或断裂追踪法分离交叉裂隙骨架,继而使用微分累加法、方框法、线性回归法求得裂隙的长度、宽度及倾角等几何特征指标。基于提出的算法,研制了一套具有图形用户界面的岩石裂隙图像智能识别与表征软件系统,实现了从深度学习模型参数选择、模型训练、裂隙识别、量化分析到结果可视化的完整功能。最后对岩石裂隙识别与分离表征算法的性能进行了评判,结果表明,改进Unet模型对复杂分布的裂隙识别效果最好,其总体识别性能要优于其他网络;骨架分离算法对常见类型交叉裂隙能够取得预期结果,表征算法对分离裂隙与交叉裂隙的表征精度高,对实际岩石裂隙图像的应用效果较好。研究成果可为基于计算机视觉的岩石工程试验与岩体结构检测技术研发提供参考依据。
    Abstract: The characteristics of rock fractures are the core indices for evaluating the structure and integrity of rock masses, and are also the important factors for evaluating the safety and stability of rock engineering. For the rock fracture recognition based on deep learning, the Unet model is improved by introducing the convolutional block attention module, effectively improving the accuracy of rock fracture recognition. For separation and feature extraction of intersection fractures, a fracture separation and characterization algorithm based on the trace direction judgment is proposed. According to the characteristics of fracture separation, the overlapping tracing method or fracture tracing method is used to separate the intersection fracture skeleton. Then, the differential accumulation method, box method and linear regression method are employed to calculate the geometric characteristic indices such as the length, width and dip angle of the fractures. Based on the proposed algorithm, a set of software system for the intelligent identification and characterization of rock fracture images with a graphical user interface is developed, achieving complete functions from deep learning model parameter selection, model training, fracture recognition, quantitative analysis to result visualization. Finally, the performance of rock fracture recognition and separation characterization algorithms is evaluated. The results show that the improved Unet model has the best recognition performance for complex distribution fractures, and its overall recognition performance is superior to that of other networks. The skeleton separation algorithm achieves the expected results for the separation of common types of intersection fractures, and the characterization algorithm has high accuracy in characterizing the separation and intersection fractures. It has good application effects on actual rock fracture images. The research can provide reference for rock engineering tests and detection of rock mass structure by using the computer vision.
  • 高水平放射性废物(以下简称高放废物)的安全处置,是能源、生态可持续发展不可忽视的重要研究内容,受到了核工业国家的高度重视。目前高放废物地质处置库多采用包含处置罐(贮存高放废物)、缓冲材料以及围岩(即包容地质处置设施的稳定地质体或岩层)的“多重屏障系统”进行设计[1],其中围岩对工程的长期稳定和安全具有关键的控制作用。由于地质构造运动的影响,围岩中不可避免地会存在各种随机分布的天然裂隙。另外处置工程在建设过程中也会导致围岩中产生裂隙或原生裂隙的扩展,形成开挖损伤区,对长久安全地封存放射性废物产生不利影响。因此,含裂隙围岩的力学行为特性研究,是高放废物地质处置工程稳定性和安全性评价的重要课题。

    工程岩体中所存在的大量裂隙会很大程度地降低岩体的强度,而地下工程的失稳破坏也与裂隙的扩展、贯通密切相关[2-3]。针对含节理灰岩试样,王在泉等[4]通过加轴压、卸围压的应力控制试验和三轴压缩试验,研究了节理岩体的强度、变形特征和破坏特征。结合声发射监测技术,王强等[5]对含弱面的辉绿岩开展了单轴压缩试验研究,分析了含不同角度弱面试件的变形、强度特性和声发射特征。针对天然裂隙试样取样困难、试验可重复性差等缺点,国内外学者通过采用激光切割、水射流等方法在真实的完整岩石的基础上预制裂隙,制成裂隙岩体,并开展了大量的研究工作[6-9]。罗丹旎等[10]通过真三轴岩爆试验,研究了不同产状预制单裂隙花岗岩的破坏模式、强度变形和声发射特征,得出了裂隙花岗岩的剪切破裂型和张拉板裂型的两种岩爆机制。Wong等[11]通过含三维表面裂隙的大理岩单轴压缩试验,研究了裂隙倾角、深度等参数对裂纹扩展的影响规律。谢其泰等[12]针对预制单裂隙砂岩试样开展了单轴压缩试验研究,分析了裂隙倾角对裂纹扩展机制的影响规律。综合而言,机械切割法或高压水射流法是制备预制裂隙的常用方法,但采用这两种方法制备的非贯通裂隙往往不可避免地存在一定的开度。由于原岩材料的刚度较大、具有较强的抗变形能力,导致外力即使达到了抗压强度,其预制裂隙也可能无法发生闭合。再者,通过切割机或高压水刀制成的裂隙切口较为平整,难以反映裂隙表面粗糙度对岩体力学强度的影响。

    本文以中国高放废物地质处置首选候选围岩甘肃北山花岗岩为研究对象,通过对压裂法获得的含不同倾角贯通单裂隙的花岗岩试样开展单轴压缩及三轴压缩试验研究,结合空间三维点云扫描信息,揭示了围压和裂隙倾角对花岗岩力学行为的影响规律,提出了沿裂隙面滑移破坏试样抗压强度与滑移前后裂隙面分形维数差值之间的经验关系式。

    本文针对完整花岗岩和含不同倾角预制单裂隙的花岗岩试样,开展了单轴压缩试验以及围压为2,5,10 MPa条件下的三轴压缩试验。试验设备为法国TOP公司所生产的岩石高温三轴流变试验系统,该系统额定荷载1500 kN;最大围压和渗透压均为60 MPa;试验最高温度达200℃。试验所用花岗岩取自甘肃北山新场地表,为中细粒花岗闪长岩。根据相关试验规程及试验方法标准,首先将岩样加工成Φ50 mm×H100 mm的标准试样,作为本文中的完整试样。随后,选取部分完整试样放置于自行研制的岩石试样通透裂隙制造模具中,通过压裂的方法制备含倾角30°,45°,60°贯通单裂隙的花岗岩试样若干,如图 1所示。通过对完整岩石压裂的方法制备的贯通裂隙面具有一定的粗糙度,并且裂隙面上下两个部分能够自然耦合接触,即使是大倾角也不会发生滑落。

    图  1  含不同倾角预制单裂隙的花岗岩试样
    Figure  1.  Granite samples with pre-existing single fractures at different inclination angles

    为降低预制裂隙面初始闭合可能引起较大变形的影响,试验开始前首先采用20 MPa左右的偏应力将试样预压紧。但考虑到倾角60°裂隙试样的单轴抗压强度较低,因此对该预制裂隙面的预压紧是在装样时人工完成的。试验过程中,首先采用5 MPa/min的加载速率将围压加载至预定值并保持不变,随后采用2 MPa/min的速率施加轴向荷载,直至试样残余变形阶段。试样的轴向变形采用2个线性位移传感器进行测量,横向变形采用环向位移引伸计测量。此外,试验前后采用Handy SCAN 3D 700便携式高速激光扫描设备对试样表面及预制裂隙面进行扫描建模,以获得不同试验条件下的裂纹分布特征。

    图 2~5为完整花岗岩及不同倾角裂隙花岗岩在单轴和三轴压缩试验后的试样照片,以及通过三维激光扫描获得的试样表面的裂纹分布特征。可以看出,完整花岗岩在单轴压缩条件下主要发生的是脆性的张拉劈裂破坏,试样表面分布诸多的竖向裂纹。2 MPa围压条件下的破坏是劈裂破坏和剪切破坏共存,但仍以劈裂破坏占主导地位,与单轴条件下不同的是,试样表面的裂纹数量显著减少。当围压超过5 MPa时,试样的破坏形式由劈裂破坏过渡到剪切破坏,且剪切破裂面和最大主应力之间的夹角随围压的增大而有所增加。

    图  2  不同围压条件下完整试样及表面裂纹
    Figure  2.  Intact specimens and cracks on surface under different confining pressures

    倾角30°裂隙花岗岩单轴和三轴条件下的破坏形式与完整岩石的破坏形式有诸多相似之处(图 3)。新生裂纹穿过预制裂隙面沿加载方向向试样两端扩展,最终在试样表面形成两条较为显著的主裂纹。在单轴条件下主要为竖向的劈裂破坏;在5 MPa和10 MPa围压条件下为典型的剪切破坏;而在2 MPa的低围压条件下,是以劈裂破坏为主,劈裂破坏和剪切破坏并存的一种破坏方式。不同的是,当微裂纹扩展至预制裂隙面时,预制裂隙面的存在导致试样在受压过程中能量释放面的增大,因此裂隙花岗岩试样表面的次生裂纹数量明显减少。

    图  3  不同围压条件下倾角30°试样及表面裂纹
    Figure  3.  Granite with fractures inclined at 30° and cracks on surface under different confining pressures

    对于倾角45°的裂隙花岗岩(图 4),单轴条件下的新生裂纹少量穿过预制裂隙面,但未贯通至试样上下两个端面。新生主裂纹主要分布于预制裂隙面靠近试样端面的区域,两侧的主裂纹未直接连通,而是由预制裂隙面隔开呈现倾斜的“Z”字型。在围压2 MPa和5 MPa条件下,新生主裂纹穿过预制裂隙面贯通至试样两端,在试样表面形成较为规则的剪切面。而当围压升高至10 MPa条件时,破坏形式发生转变,呈现出与单轴条件下相似的倾斜“Z”字型的剪切-张拉破坏,破坏形式为穿裂隙面破坏和沿裂隙面滑移并存的复合破坏模式。

    图  4  不同围压条件下倾角45°试样及表面裂纹
    Figure  4.  Granite with fractures inclined at 45° and cracks on surface under different confining pressures

    对于倾角60°的裂隙花岗岩(图 5),在单轴和三轴条件下的破坏形式较为统一,均是沿预制裂隙面发生滑移破坏。在单轴及低围压(2 MPa)条件下,仅在预制裂隙面靠近试样端部的区域观测到少许裂纹,这是由于在预制裂隙面摩擦滑移的过程中,预制裂隙长轴端部位置可能形成薄弱的凌空面,从而发生岩片的脱落。而在高围压条件下,升高的围压会对试样薄弱面的张拉剥落产生制约作用,因此破坏后试样表面除预制裂隙外几乎观测不到新生裂纹。

    图  5  不同围压条件下倾角60°试样及表面裂纹
    Figure  5.  Granite with fractures inclined at 60° and cracks on surface under different confining pressures

    图 6为不同围压条件下完整及裂隙花岗岩变形破坏过程中的应力-应变曲线。单轴条件下,倾角30°裂隙花岗岩发生穿裂隙面破坏,应力-应变曲线的弹性阶段缩小,较早地进入裂纹稳定扩展阶段。当轴向应力约为试件峰值应力的81.67%时,体积应变开始呈现出由压缩向膨胀的过渡,试样进入裂纹不稳定发展阶段。峰后阶段应力迅速回落,表明30°倾角裂隙的存在没有改变试样脆性劈裂破坏的属性,试样的变形主要为岩石材料的压缩变形。倾角60°裂隙花岗岩峰前的变形呈现出明显的屈服阶段,损伤应力非常接近峰值应力,平均约为峰值应力的88.43%。但在由该应力水平发展至峰值应力的时间段内,轴向及环向变形迅速增加,滑移特征明显。峰值后应力水平未出现立即跌落,倾角60°的预制裂隙面使试样在峰值滑移破坏后仍是两个单独承载的块体,这两个块体继续滑移致使轴向、横向变形的迅速增长和轴向荷载的缓慢下降,试样的变形主要为结构面的变形。

    图  6  不同围压条件下裂隙花岗岩应力-应变曲线
    Figure  6.  Stress-strain curves of fractured granite under different confining pressures

    对于倾角45°裂隙花岗岩,其破坏形式较为复杂,呈现出穿裂隙面破坏和沿裂隙面滑移并存的复合破坏模式。对应在应力-应变曲线上,倾角45°裂隙花岗岩的变形也介于30°和60°裂隙花岗岩的变形之间,其变形特征更是揉合了之前二者的变形特征,试样的变形是包含了岩石的压缩变形和结构面变形的耦合变形。首先,倾角45°裂隙花岗岩的峰前变形特征更接近于倾角30°裂隙花岗岩的变形,只是塑性变形阶段更为明显。但倾角45°裂隙花岗岩的损伤应力也非常接近峰值应力,平均值约为峰值应力均值的88.67%,这一点与倾角60°裂隙花岗岩的变形特征更为相似,不同的是倾角60°裂隙花岗岩的损伤应力后试件的轴向应变和横向应变在很短的时间发生了较大的变形。其次,倾角45°裂隙花岗岩的峰后的变形特征与倾角60°的更为接近,峰后未出现应力水平的立刻跌落。

    在三轴应力条件下,低倾角(30°,45°)裂隙花岗岩最终的破坏特征与完整岩石的破坏模式相似,其变形发展曲线与完整花岗岩的变形发展特征也相近,说明围压的存在限制了预制裂隙面的劣化影响,试样变形的主控因素为岩石材料的压缩变形。随着围压的增加,完整及低倾角裂隙花岗岩的变形逐渐从脆性向塑性转化。对于倾角60°裂隙花岗岩,三轴压缩条件下试样变形的主控因素为结构面的变形,峰值附近的非线性变形阶段随围压的升高显著减弱。与单轴条件下相比,2 MPa和5 MPa围压条件下损伤应力占峰值应力的比值显著降低,分别为73.91%和63.54%,当围压升高至10 MPa时,该比值又上升为83.32%。

    基于试样的应力-应变曲线,获得不同围压条件下试样的峰值应力对应的轴向应变,如图 7所示。

    图  7  峰值破坏时的轴向应变
    Figure  7.  Axial strains at peak failure

    结果表明,裂隙花岗岩峰值破坏时的轴向应变随围压的升高呈现出线性增加的趋势。低倾角裂隙花岗岩峰值破坏时的轴向应变与完整岩石的相差不大,在围压0,2,5,10 MPa条件下该值约为0.38%,0.49%,0.59%和0.75%。但当裂隙倾角增大至60°时,该值在前述围压条件下分别为0.39%,0.38%,0.49%,0.54%,三轴条件下峰值破坏时的轴向应变与完整及低倾角裂隙花岗岩相比存在一个明显的下降,且围压越高,下降越显著。

    不同围压条件下裂隙花岗岩的抗压强度如图 8所示。单轴压缩条件下(图 8(a)),完整花岗岩及不同倾角的裂隙花岗岩各开展了4组试样的测试。测试所得完整花岗岩的平均单轴抗压强度为155.12 MPa,倾角30°,45°,60°裂隙花岗岩的平均单轴抗压强度分别为142.76,109.86,43.93 MPa,裂隙倾角越大,抗压强度越低。结果还发现,对于完整花岗岩和低倾角(30°和45°)裂隙花岗岩,不同试样的抗压强度的差异性并不显著。但对于倾角60°的裂隙花岗岩,不同试样的抗压强度的差异性较大。鉴于此,在三轴压缩时,完整花岗岩和低倾角裂隙花岗岩在各个围压条件下仅开展了1组测试,但倾角60°裂隙花岗岩分别开展了3组测试。

    图  8  不同围压条件下裂隙花岗岩抗压强度
    Figure  8.  Compressive strengths of fractured granite under different confining pressures

    结果表明(图 8(b)),完整花岗岩及低倾角裂隙花岗岩的三轴抗压强度随围压升高呈现出线性增长的趋势,且围压的存在限制了裂隙对岩石强度的劣化作用,导致三轴压缩条件下低倾角裂隙花岗岩的抗压强度略低于完整岩石的抗压强度但整体相差不大,试样三轴抗压强度的平均值在2,5,10 MPa围压条件下分别为189.84,225.20,267.28 MPa。此时试样的破坏模式与完整岩石的三轴破坏模式相似,新生裂纹大多是穿过裂隙面贯通至试样的上下两个端面。当裂隙倾角增大至60°时,在前述围压条件下的三轴抗压强度分别为60.54,82.94,195.78 MPa。试样的抗压强度随围压升高呈现出非线性增长的趋势,且各围压条件下的抗压强度均显著低于完整花岗岩及低倾角裂隙花岗岩的抗压强度,相对应地试样发生的是沿裂隙面滑移的破坏模式。此外,由于围压的限制作用,不同试样抗压强度的差异性有所降低。

    图 9~11为裂隙花岗岩单轴及单轴压缩试验后预制裂隙面的状态。除试样表面呈现出分布特征不同的裂隙以外,试样内部的预制裂隙面也呈现出显著的变化。为分析预制裂隙面粗糙度的变化特征,在试验前后分别采用三维激光扫描仪对预制裂隙面进行了扫描,获得了预制裂隙面的三维点云数据。

    图  9  试验后倾角30°预制裂隙面
    Figure  9.  Fracture surfaces inclined at 30° after tests
    图  10  试验后倾角45°预制裂隙面
    Figure  10.  Fracture surfaces inclined at 45° after tests
    图  11  试验后倾角60°预制裂隙面
    Figure  11.  Fracture surfaces inclined at 60° after tests

    由于低倾角裂隙花岗岩试样表面的新生裂纹或穿过预制裂隙面贯通至试样上下两个端面,或由预制裂隙面隔开形成倾斜的“Z”字型,对应的预制裂隙面上的裂纹为张拉型或剪切型。因此预制裂隙面在单轴压缩条件下通常较为破碎,而在三轴条件下裂隙面又较宽,在拆样的过程中不可避免地会造成预制裂隙面上细碎岩块的脱落,从而对三维激光扫描结果造成误差,直接影响了试验后分形维数计算的准确性,因此本研究中未给出倾角30°和45°裂隙面的分形维数。而对于倾角60°的裂隙花岗岩,其破坏模式为沿裂隙面滑移破坏,试验后预制裂隙面较完整,有利于定量衡量试验前后粗糙度的变化。另外,在上一节的分析中提到,不同倾角60°裂隙花岗岩试样的抗压强度的差异性较大,因此本节基于三维激光扫描结果,重点对倾角60°预制裂隙面粗糙度的变化进行分析。

    关于裂隙面几何特征的表征方法,较为经典的是Barton提出的经典剖面法,包括用于确定结构面粗糙度系数JRC的10条标准剖面[13]。但由于采用该方法来比照衡量被测表面的几何形貌时,可能会因为较强的主观性而导致判断结果产生比较大的偏差。于是,相关领域的学者在此基础上提出了不同的JRC描述方法[14-16]。在定量化描述裂隙岩体结构面几何特征方面,20世纪70年代创立的分形几何,又为裂隙面的衡量提供了一个启示性的理论方法和手段[17]。根据以往的研究结果,分形维数可以作为描述岩体中裂隙面几何形貌特征复杂程度的定量指标[18-21]。与经典几何相比,分形维数具有无标度的特点,其描述的裂隙面的几何特征不受测量尺度和精度的限制。基于上述分析,本研究采用Zhou等[22]提出的立方体覆盖法计算试验前后预制裂隙面的分形维数,探讨裂隙面抗压强度与其破坏特征之间的关系。如图 12所示,通过改变单元格边长δ的大小,得到不同覆盖立方体的个数N(δ)。如果裂隙面的粗糙度表现出分形特征,则边长δ与个数N(δ)之间存在如下关系:

    N(δ)δDf
    (1)
    图  12  立方体覆盖法计算分形维数示意图
    Figure  12.  Schematic diagram of calculating fractal dimension by cube covering method

    式中:Df为表征裂隙面粗糙度的分形维数。

    单轴压缩或三轴压缩试验后,倾角60°预制裂隙面上均可观测到大片的摩擦滑移区,对应于三维点云图中(图 13),可明显地看出试验后预制裂隙面变得更加平缓。采用图 12的立方体覆盖法计算试验前后预制裂隙面的分形维数,结果如图 14所示。研究表明,不同围压条件下,预制裂隙面发生摩擦滑移破坏后其分形维数与试验前相比均有所下降;且围压越高,试验前后预制裂隙面的分形维数差值越大,说明摩擦滑移破坏后预制裂隙面粗糙度的变化越显著。图 14(b)为单轴及三轴压缩条件下试样抗压强度与试验前后裂隙面分形维数差的关系。研究发现,试样的抗压强度与试验前后的分形维数差值符合幂函数y=axb+c的增长关系,式中变量yx分别代表抗压强度和试验前后裂隙面的分形维数差值,参数abc为拟合参数。分形维数相差越大,试样的抗压强度越高,这也从侧面解释了相同试验条件下不同倾角60°裂隙花岗岩试样抗压强度存在差异性的原因。

    图  13  倾角60°预制裂隙面三维点云建模(σ3= 10 MPa)
    Figure  13.  Three-dimensional point cloud modeling of fracture surface inclined at 60° (σ3= 10 MPa)
    图  14  倾角60°预制裂隙面分形维数
    Figure  14.  Fractal dimensions of fracture surface inclined at 60°

    本文通过完整花岗岩和含不同倾角预制裂隙的花岗岩试样的单轴压缩和三轴压缩试验研究,揭示了围压及裂隙倾角对花岗岩力学行为的影响机制,得到以下3点结论。

    (1)花岗岩的破坏形式受围压与裂隙倾角双重因素的共同影响。完整花岗岩的破坏形式随围压的升高由劈裂破坏向剪切破坏过渡,且围压越高剪切破裂面和最大主应力之间的夹角越大。倾角30°裂隙花岗岩的新生裂纹穿过预制裂隙面向试样两端扩展形成穿裂隙面的破坏模式,与完整花岗岩的破坏形式相近但裂隙数量明显减少。倾角60°裂隙花岗岩在0~10 MPa围压条件下均发生沿预制裂隙面的摩擦滑移破坏,试样表面除预制裂隙外几乎观测不到新生裂纹。倾角45°裂隙花岗岩在2 MPa和5 MPa围压条件下,新生主裂纹穿过预制裂隙面贯通至试样两端,在试样表面形成较为规则的剪切面;但在单轴及10 MPa围压条件下,新生主裂纹沿预制裂隙面呈现出倾斜“Z”字型,破坏形式为穿裂隙面破坏和沿裂隙面滑移并存的复合破坏模式。

    (2)单轴压缩条件下,裂隙倾角越大,试样峰值破坏前的非线性屈服特征越明显。围压的存在限制了预制裂隙面的劣化作用。在三轴应力条件下,低倾角裂隙花岗岩的变形发展曲线与完整花岗岩的变形发展特征相近,试样峰值破坏时的轴向应变在围压0,2,5,10 MPa条件下分别约为0.38%,0.49%,0.59%,0.75%。对于倾角60°裂隙花岗岩,随围压的升高,峰值附近的非线性变形阶段显著减弱,且峰值破坏时对应的轴向应变与完整及低倾角裂隙花岗岩相比存在一个明显的下降,且围压越高,下降越显著。

    (3)倾角越大,裂隙花岗岩的单轴抗压强度越低。三轴压缩条件下低倾角裂隙花岗岩的抗压强度略低于完整岩石的抗压强度,但整体相差不大。倾角60°裂隙花岗岩的抗压强度随围压升高呈现出非线性增长的趋势,且各围压条件下的抗压强度均显著低于完整及低倾角裂隙花岗岩的抗压强度。倾角60°预制裂隙面发生摩擦滑移破坏后分形维数有所下降;且围压越高,试验前后预制裂隙面的分形维数差值越大。试样的抗压强度与试验前后的分形维数差值符合幂函数增长关系,分形维数相差越大,试样的抗压强度越高。

  • 图  1   数据集内部分原图与对应的标签图

    Figure  1.   Partial original and label images in dataset

    图  2   MRC-Unet模型结构示意图

    Figure  2.   Structural diagram of MRC-Unet model

    图  3   CBAM结构示意图

    Figure  3.   Schematic diagram of CBAM structure

    图  4   Loss-Epoch关系曲线

    Figure  4.   Loss-Epoch relation curves

    图  5   裂隙骨架特征点的判别原理

    Figure  5.   Discriminant principle of characteristic points of fracture skeleton

    图  6   不同类型裂隙交叉点的判别

    Figure  6.   Distinguishing of different types of intersection points of fractures

    图  7   交叉裂隙两种分离结果示意图

    Figure  7.   Schematic diagram of separation results of intersection fractures

    图  8   交叉点处裂隙追踪方向判定

    Figure  8.   Distinguishing of fracture tracking direction at intersection point

    图  9   主干裂隙判别原理

    Figure  9.   Principle of discrimination of main fracture

    图  10   裂隙长度计算示意图

    Figure  10.   Schematic diagram of calculation of crack length

    图  11   裂隙宽度计算示意图

    Figure  11.   Schematic diagram of calculation of crack width

    图  12   软件功能构架图

    Figure  12.   Function framework of software

    图  13   PhotoDetector程序的模型选择与训练界面

    Figure  13.   Model selection and training interface of PhotoDetector program

    图  14   PhotoDetector程序的裂隙识别界面

    Figure  14.   Fracture identification interface of PhotoDetector program

    图  15   PhotoObjects的裂隙表征界面

    Figure  15.   Fracture characterization interface of PhotoObjects program

    图  16   模拟裂隙长度与宽度的计算结果

    Figure  16.   Calculated results of simulated fracture length and width

    图  17   实际裂隙表征结果示意图

    Figure  17.   Schematic diagram of characterization results of actual fractures

    表  1   岩石裂隙基本类型划分

    Table  1   Classification of basic types of rock fractures

    裂隙类别 裂隙图像
    单条 多条
    分离型
    Y型 X型
    单点交叉型
    多点交叉型
    下载: 导出CSV

    表  2   各类算法裂隙识别效果比较

    Table  2   Comparison of fracture recognition effects of various algorithms

    类别 原图 基准图 VGG16 Unet 改进Unet
    类别1
    类别2
    类别3
    下载: 导出CSV

    表  3   不同算法的裂隙识别效果定量评价

    Table  3   Quantitative evaluation of fracture recognition effects under different algorithms

    评价指标 算法模型
    VGG16 Unet 改进Unet
    F1分数 0.924 0.933 0.943
    Iou 0.859 0.874 0.893
    下载: 导出CSV

    表  4   常见交叉裂隙骨架分离结果

    Table  4   Skeleton separation results of common intersection fractures

    图像类别 裂隙骨架原图 重合追踪分析结果 断裂追踪分析结果
    图像A(Y型)
    图像B(X型)
    图像C(双Y型)
    图像D(三角交叉型)
    下载: 导出CSV
  • [1] 张紫杉, 王述红, 王鹏宇, 等. 岩坡坡面裂隙网络智能识别与参数提取[J]. 岩土工程学报, 2021, 43(12): 2240-2248. doi: 10.11779/CJGE202112010

    ZHANG Zishan, WANG Shuhong, WANG Pengyu, et al. Intelligent identification and extraction of geometric parameters for surface fracture networks of rocky slopes[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(12): 2240-2248. (in Chinese) doi: 10.11779/CJGE202112010

    [2] 李术才, 刘洪亮, 李利平, 等. 基于数码图像的掌子面岩体结构量化表征方法及工程应用[J]. 岩石力学与工程学报, 2017, 36(1): 1-9. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201701001.htm

    LIU Shucai, LIU Hongliang, LI Liping, et al. A quantitative method for rock structure at working faces of tunnels based on digital images and its application[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(1): 1-9. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201701001.htm

    [3]

    LI Y H, TANG X J, ZHU H H. Optimization of the digital image correlation method for deformation measurement of geomaterials[J]. Acta Geotechnica, 2022, 17(12): 5721-5737. doi: 10.1007/s11440-022-01646-x

    [4] 许文涛, 李晓昭, 章杨松, 等. 基于摄影测量系统的岩体结构面精细识别表征及应用[J]. 测绘学报, 2022, 51(10): 2093-2106. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202210009.htm

    XU Wentao, LI Xiaozhao, ZHANG Yangsong, et al. Fine identification and characterization of rock mass discontinuities and its application using a digital photogrammetry system[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(10): 2093-2106. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB202210009.htm

    [5]

    LECUN Y, BENGIO Y, HINTON G. Deep learning[J]. Nature, 2015, 521(7553): 436-444. doi: 10.1038/nature14539

    [6]

    LONG J, SHELHAMER E, DARRELL T. Fully convolutional networks for semantic segmentation[C]//Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Boston, 2015: 3431-3440.

    [7] 黄宏伟, 李庆桐. 基于深度学习的盾构隧道渗漏水病害图像识别[J]. 岩石力学与工程学报, 2017, 36(12): 2861-2871. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201712001.htm

    HUANG Hongwei, LI Qintong. Image recognition for water leakage in shield tunnel based on deep learning[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(12): 2861-2871. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201712001.htm

    [8]

    ZHU H, AZARAFZA M, AKGUN H. Deep learning-based key-block classification framework for discontinuous rock slopes[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2022, 14: 1131-1139. doi: 10.1016/j.jrmge.2022.06.007

    [9] 薛东杰, 唐麒淳, 王傲, 等. 基于FCN的岩石混凝土裂隙几何智能识别[J]. 岩石力学与工程学报, 2019, 38(增刊2): 3393-3403. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2019S2014.htm

    XUE Dongjie, TANG Qichun, WANG Ao, et al. FCN-based intelligent identification of crack geometry in rock or concrete[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(S2): 3393-3403. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2019S2014.htm

    [10] 刘金杉, 李元海, 卢昱杰, 等. 基于深度学习的隧道衬砌表观病害模拟检测系统[J]. 哈尔滨工业大学学报, 2022, 54(5): 24-33.

    LIU Jinshan, LI Yuanhai, LU Yujie, et al. Tunnel lining surface defect simulation and detection system based on deep learning[J]. Journal of Harbin Institute of Technology, 2022, 54(5): 24-33. (in Chinese)

    [11]

    SONG Q, WU Y, XIN X, et al. Real-time tunnel crack analysis system via deep learning[J]. IEEE Access, 2019, 7: 64186-64197. doi: 10.1109/ACCESS.2019.2916330

    [12]

    XIANG X, ZHANG Y, EL SADDIK A. Pavement crack detection network based on pyramid structure and attention mechanism[J]. IET Image Processing, 2020, 14(8): 1580-1586. doi: 10.1049/iet-ipr.2019.0973

    [13] 刘春, 王宝军, 施斌, 等. 基于数字图像识别的岩土体裂隙形态参数分析方法[J]. 岩土工程学报, 2008, 30(9): 1383-1388. doi: 10.3321/j.issn:1000-4548.2008.09.021

    LIU Chun, WANG Bao-jun, SHI Bin, et al. Analytic method of morphological parameters of cracks for rock and soil based on image processing and recognition[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(9): 1383-1388. (in Chinese) doi: 10.3321/j.issn:1000-4548.2008.09.021

    [14] 王军祥, 曾相森, 徐晨晖, 等. 基于图像处理技术的岩体裂隙定量识别方法研究[J]. 地下空间与工程学报, 2022, 18(2): 446-457. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202202012.htm

    WANG Junxiang, ZENG Xiangsen, XU Chenhui, et al. Study on quantitative identification method of rock fracture based on image processing technology[J]. Chinese Journal of Underground Space and Engineering, 2022, 18(2): 446-457. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202202012.htm

  • 期刊类型引用(2)

    1. 黄勇钧,李世文,刘先林,邵羽,翁剑,任中俊,吕玺琳. 裂隙花岗岩力学特性试验及破坏条件理论模型. 力学季刊. 2025(01): 162-172 . 百度学术
    2. 戎虎仁,曹海云,董浩,王占盛,张佳瑶. 双孔洞倾角对红砂岩强度及其破裂分形特征影响. 公路交通科技. 2024(10): 173-183 . 百度学术

    其他类型引用(1)

图(17)  /  表(4)
计量
  • 文章访问数:  1101
  • HTML全文浏览量:  133
  • PDF下载量:  388
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-10-09
  • 网络出版日期:  2024-03-14
  • 刊出日期:  2024-02-29

目录

/

返回文章
返回