CN110940734B - 岩体内部异常结构监测与潜在地质灾害评估方法及*** - Google Patents
岩体内部异常结构监测与潜在地质灾害评估方法及*** Download PDFInfo
- Publication number
- CN110940734B CN110940734B CN201911226383.3A CN201911226383A CN110940734B CN 110940734 B CN110940734 B CN 110940734B CN 201911226383 A CN201911226383 A CN 201911226383A CN 110940734 B CN110940734 B CN 110940734B
- Authority
- CN
- China
- Prior art keywords
- rock mass
- ultrasonic
- geological disaster
- ultrasonic imaging
- potential geological
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/26—Arrangements for orientation or scanning by relative movement of the head and the sensor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种岩体内部异常结构监测与潜在地质灾害评估方法及***,解决了现有技术无法实现可视化监测和对潜在地质灾害准确评估的问题。一种岩体内部异常结构可视化监测与潜在地质灾害评估方法,获取目标区域内所有位置、角度的岩体结构超声成像剖面图,获得处理后的岩体结构超声成像剖面图;对所述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值L;根据所述S、N、L评估岩体内部异常结构潜在地质灾害程度。实现了岩体内部异常结构的可视化监测与潜在地质灾害的准确评估。
Description
技术领域
本发明涉及岩土工程中地质勘察与测量技术领域,尤其是涉及一种岩体内部异常结 构监测与潜在地质灾害评估方法及***。
背景技术
地质灾害往往始于局部岩土结构的侵蚀破坏与受力失衡,特别是水作用下的软弱结 构或土-岩交界区;因此,岩土体软弱结构区的地质勘察与监测诊断是防灾减灾的重要手段;传统的地质雷达或声波探测仪等设备,目前只能实现较大范围的地质勘查,针对 岩土工程中特定区域的软弱夹层或软弱结构面,传统方法在精度范围上难以达到,也无 法实现可视化监测和对潜在地质灾害准确的评估,严重影响了岩土工程或边坡工程的稳 定性和安全性,迫切需要一种新的勘察手段加以可视化监测与潜在地质灾害准确的评 估。
发明内容
本发明的目的在于至少克服上述一种技术不足,提出了一种岩体内部异常结构监测 与潜在地质灾害评估方法及***。
一方面,本发明提供了一种岩体内部异常结构监测与潜在地质灾害评估方法,括以 下步骤:
获取目标区域内所有位置、角度的岩体结构超声成像剖面图,对所述所有位置、角度的超声成像剖面图进行组合及拼接融合处理,获得处理后的岩体结构超声成像剖面 图;
对所述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区 域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值L;
根据所述S、N、L获取潜在地质灾害可能性参数,并根据所述地质灾害可能性参 数评估岩体内部异常结构潜在地质灾害程度。
进一步地,所述获取目标区域内多个位置角度的岩体结构超声成像剖面图,具体包 括,将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过旋转、滑动或者侧向倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描成像图,对所 述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角 度的岩体结构超声成像剖面图。
进一步地,根据所述S、N、L获取潜在地质灾害可能性参数,具体包括,利用公 式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的总区域,T为修正系数,N≥3。
进一步地,根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度,具体包括,地质灾害可能性参数越大岩体内部异常结构潜在地质灾害程度越大。
另一方面,本发明还提供了一种岩体内部异常结构监测与潜在地质灾害评估***, 包括超声成像剖面图获取设备、超声成像剖面图处理设备、图像特征进行统计模块以及 潜在地质灾害程度评估模块;
所述超声成像剖面图获取设备,用于获取目标区域内所有位置、角度的岩体结构超 声成像剖面图;
所述超声成像剖面图处理设备,用于对所述所有位置、角度的超声成像剖面图进行 组合及拼接融合处理,获得处理后的岩体结构超声成像剖面图;
所述图像特征进行统计模块,用于对所述处理后的岩体结构超声成像剖面图中二值 化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值 L;
所述潜在地质灾害程度评估模块,用于根据所述S、N、L获取潜在地质灾害可能 性参数,并根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度。
进一步地,所述超声成像剖面图获取设备包括超声扫描成像探头、推杆和图像处理 模块;超声成像剖面图获取设备获取目标区域内所有位置、角度的岩体结构超声成像剖面图,具体包括,
通过推杆将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过推杆旋转、 滑动或者侧向倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描 成像图,所述图像处理模块对所述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角度的岩体结构超声成像剖面图。
进一步地,所述潜在地质灾害程度评估模块根据所述S、N、L获取潜在地质灾害 可能性参数,具体包括,利用公式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的总区域,T为修正系数,N≥3。
进一步地,所述潜在地质灾害程度评估模块,根据所述地质灾害可能性参数评估岩 体内部异常结构潜在地质灾害程度,具体包括,所述地质灾害可能性参数越大,评定岩体内部异常结构潜在地质灾害程度越大。
与现有技术相比,本发明的有益效果包括:通过获取目标区域内所有位置、角度的岩体结构超声成像剖面图,对所述所有位置、角度的超声成像剖面图进行组合及拼接融 合处理,获得处理后的岩体结构超声成像剖面图;对所述处理后的岩体结构超声成像剖 面图中二值化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和 条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点 个数的比值L;根据所述S、N、L获取潜在地质灾害可能性参数,并根据所述地质灾害 可能性参数评估岩体内部异常结构潜在地质灾害程度;实现了岩体内部异常结构的可视 化监测与潜在地质灾害的准确评估。
附图说明
图1是本发明实施例1所述的岩体内部异常结构监测与潜在地质灾害评估方法的流 程示意图;
图2是本发明实施例1所述的岩体孔壁剖面的超声扫描成像图示意图;
图3是本发明实施例2所述的钻孔内孔壁岩体内部结构的可视化监测示意图;
图4是本发明实施例2所述的钻孔内孔壁岩体内部结构剖面成像效果图;
图5是本发明实施例3所述的超声扫描成像探头固定方式示意图;
图6是本发明实施例3所述的岩体内部结构的可视化监测示意图;
图7是本发明实施例3所述的三种超声扫描成像示意图;
图8是本发明实施例3所述的多个位置、角度获得的超声扫描图像示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1
本发明提供了一种岩体内部异常结构监测与潜在地质灾害评估方法,其流程示意图,如图1所示,所述岩体内部异常结构监测与潜在地质灾害评估方法,包括以下步骤:
获取目标区域内所有位置、角度的岩体结构超声成像剖面图,对所述所有位置、角度的超声成像剖面图进行组合及拼接融合处理,获得处理后的岩体结构超声成像剖面 图;
对所述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区 域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值L;
根据所述S、N、L获取潜在地质灾害可能性参数,并根据所述地质灾害可能性参 数评估岩体内部异常结构潜在地质灾害程度。
需要说明的是,所述岩体内部异常结构是指岩体结构面中除正常岩石之外的结构面、破碎带、孔洞、软弱夹层等结构,这些异常结构将在超声扫描图像中所呈现出来异 常斑点、条纹等成像特征;
优选的,所述获取目标区域内多个位置角度的岩体结构超声成像剖面图,具体包括,
将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过旋转、滑动或者侧向 倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描成像图,对所述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角 度的岩体结构超声成像剖面图。
具体实施时,针对钻孔内部的岩体结构,超声扫描成像仪的探头可以在钻井平台或 特定设备的帮助下,借助于推杆或滑竿和夹子等辅助装置,把超声扫描成像探头固定在支架上并放入钻孔内部岩体孔壁的目标区域中;超声扫描成像的探头可以实现上下左右移动和360度自由旋转;
由于钻孔内部一般都有水或者已注水,或者注入了特定耦合剂,基于超声相控阵成 像原理的超声扫描成像探头可以实现岩体孔壁内部的超声扫描成像;岩体孔壁剖面的超 声扫描成像图示意图,如图2所示;其中,超声波的频率和相关配置的调整,以超声波 可以穿透岩体内部结构1米左右并可以相对清晰的得到超声扫描成像图为准;
针对超声扫描成像仪实时获取的大量超声扫描成像图,可以利用图像处理方法自动 识别超声扫描成像图中的异常结构,并进行定性定量分析,具体的,首先利用超声扫描成像仪(超声扫描成像探头)直接获取超声扫描成像视频*.cin文件,然后,利用高斯滤 波方法和开闭运算方法去除超声扫描成像图中的高频和低频噪声信号,滤除噪声,保留 岩体结构中结构面、破碎带等异常结构的超声图像,得到异常结构的二值化图像,即为 岩体结构超声成像剖面图;接着统计分析二值化图像中的超声光斑点数和边界条纹个数 及其所占的总面积;最后,通过监测这些光斑和边界条纹的变化情况,从而实现岩体内 部异常结构的可视化监测。
优选的,根据所述S、N、L获取潜在地质灾害可能性参数,具体包括,利用公式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的总区域,T为修正系数,N≥3。
优选的,根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度, 具体包括,地质灾害可能性参数越大,评定岩体内部异常结构潜在地质灾害程度越大。
具体实施例时,针对所述处理后的岩体结构超声成像剖面图,统计分析每一个角度 位置岩体内部结构超声成像剖面图中二值化图像特征,把这些二值化特征统归为岩体内 部的异常结构,具体的,针对岩体内容结构超声成像剖面图中二值化图像特征,首先,统计非零区域的总像素点数,即等价于岩体内部剖面的异常结构总面积S;然后,统计 非零不连续光斑块和条纹带的总个数N,并统计所有光斑块和条纹带的边界长度与所有 光斑块和条纹带所占像素点个数的比值,简称之为占长比L;最后,监测二值化图像中 这些光斑块和条纹带的变化情况,可二值化统计特征和岩体内部异常结构中存在的对应 关系;
岩体内部异常结构与二值化统计特征之间的对应关系为:若统计的异常结构总面积 S数值大,则说明岩体中存在杂质碎块区域较多;若统计的光斑块和条纹带总个数N较多,则说明岩体中存在杂质碎块数量较多;若统计的光斑块和条纹带总占长比L较大, 则反映了岩体中存在较多的裂隙结构面或者较为规整的破碎区域;
因此,可以根据这些对应关系建立一个定性定量评估岩体结构潜在地质灾害的潜在 地质灾害可能性参数K,该参数K在一定条件下可以反映出所探测区域岩体结构潜在地质灾害程度;
潜在地质灾害可能性参数K为,
其中,T为修正系数,修正系数T值是根据理想条件下完整岩体结构中超声扫描成像图像处理结果来确定,具体是超声成像剖面图的图像对比度来确定,图像对比度越大 T值越小,同时,T还需满足理想条件下完整岩体结构的潜在地质灾害可能性参数K的 值为1,T取值范围为1~100;T值及超声扫描成像仪的标准参数率定校准过程需要在获 取超声成像剖面图之前进行;
当统计的光斑块和条纹带总个数N小于3时,潜在地质灾害可能性参数K的关系 式不成立,此时只能重新调整设备更换或者图像处理方法;
在N≥3的条件下,钻孔内部岩体异常结构潜在地质灾害与K存在一定的正比例关系,即K值越大,则潜在地质灾害就越大;
当K小于或等于1,则判定该区域存在地质灾害的可能性较小,岩体结构相对较为稳定;当K大于1并且小于10,则判定该区域存在一定的潜在地质灾害,岩体结构相 对较为破碎;若此时L值较大并且S/N较小,则该区域存在地质滑坡的风险;若此时L 值较小并且S/N较大,则该区域存在地质塌方的风险;该区域整体存在一定的潜在地质 灾害风险;当K大于10并且小于或等于100,则判定该区域存在的潜在地质灾害风险 较大,岩体结构相对不稳定;需要说明的是,当K大于100,则判定K值的计算过程存 在错误,修正系数T值的确定存在较大的误差;
实施例2
本发明提供了一种岩体内部异常结构监测与潜在地质灾害评估方法,采用Apogee1100全数字彩色多普勒超声诊断***为超声扫描成像仪,以5MHz的凸阵探头为最佳的 超声频率和产生扫描成像探头,在配置好该型号的超声扫描成像仪和探头后,把超声扫 描成像探头固定在推杆上,并放入在孔深为11.5米的钻孔下,此时,孔内充满了浑浊的 泥浆水,并且浓度较高,能够使超声扫描成像探头较好地接触钻孔孔壁,并在必要情况 下投入特制的泥浆粉充当超声探头与岩体孔壁的耦合剂,修正系数T值设置为8.27;所 述岩体内部异常结构监测与潜在地质灾害评估方法,包括以下步骤,
步骤S1、搭建钻孔内岩体孔壁超声扫描成像透试验***及其相关配件,比如现场测试的电源、线缆、支架、推杆等配件,并做好正式超声扫描成像前的准备工作,比如 调整好测试前超声扫描成像仪的最佳设置,使超声扫描成像仪能够立马进行最佳的工作 状态;
步骤S2、把超声扫描成像探头固定到第一根推杆上去,连续依次内管穿线接上五根推杆,并固定好第五根推杆的位置,读取第五推杆的读数和旋转角度分别为1.5米和 45度,则此时超声扫描成像探头在钻孔内的位置为9.5米和北偏东45度的位置方向;
步骤S3、连通电源,并打开超声扫描成像仪;此时可以在超声扫描成像仪的屏幕上看到钻孔内超声成像探头发出的超声波对孔壁内部结构的透射成像情况;比如可以看到钻孔内部的杂碎破损区域所形成的超声成像斑点和岩体结构面所形成的带状条纹结构;超声扫描成像***调整好后,若当前角度位置的超声扫描图像清晰可见之后,即可 保存该角度位置的超声视频图像,并通过数据线实时传输超声扫描视频图像*.cin文件;
步骤S4、针对传输回来的超声扫描成像视频图像*.cin文件,将超声扫描成像视频图像解析为一帧帧连续的单张超声扫描图像,针对每一帧超声扫描图像,分析图像中岩 体孔壁内部结构的成像特征及其对应关系;
首先,利用高斯滤波方法和开闭运算方法去除超声扫描图像中的高频和低频噪声信 号,滤除噪声,保留岩体结构中结构面、破碎带等异常结构的超声图像,得到异常结构的二值化图像;然后,统计岩体内部剖面的异常结构总面积;然后,统计非零不连续光 斑块和条纹带的总个数,并统计所有光斑块和条纹带的边界长度与所有光斑块和条纹带 所占像素点个数的比值,最后,通过监测这些光斑和边界条纹的变化情况,从而实现岩 体内部异常结构的可视化监测;
步骤S5、为了能够多方面多角度地监测钻孔内孔壁岩体内部的异常结构,故有次序缓慢地移动的推杆位置或旋转推杆,进而实现超声成像探头在钻孔内部孔壁的平移、 倾斜扫描成像和钻孔底部的旋转扫描成像,得到钻孔内不同位置不同角度下多个钻孔孔 壁岩体内部结构的超声扫描成像图像;依次保存不同位置不同角度下的多个超声扫描视 频图像*.cin文件和可视化监测结果;依照由深到浅、由大到小的顺序进行多个位置角度 岩体结构超声成像剖面图的排序组合,部分重叠区域进行必要的图像拼接和信息融合; 最后得到大范围的岩体内部结构信息和多维空间岩体剖面结构信息,并依次保存好图像 数据。
步骤6S、等到所有目标区域的测点全部完成后,抽拔出推杆及其固定部位的超声扫描成像探头,收拾好设备仪器,打包好所有的视频图像数据,以及深度位置和旋转角 度等记录文件信息。
步骤S7、根据现场打包回来视频图像数据和记录文件,分析统计每一个角度位置岩体结构超声成像剖面图中二值化图像特征,这些二值化特征即可统归为岩体内部的异常结构特征,并观测每个二值化图像特征与原始视频中对应的岩体结构之间的对应关 系。
步骤S8、视频图像数据和记录文件,以每个孔内角度位置为单位,依次统计每帧超声扫描图像的非零区域的总像素点数S1~Sn,非零不连续光斑块和条纹带的总个数 N1~Nn和光斑块和条纹带的的占长比L1~Ln,并计算它们的平均值得到S1~Sn的平均值为2568.31,N1~Nn的平均值为16.42,L1~Ln的平均值为738.29。钻孔内孔壁岩体内部结构 的可视化监测示意图,如图3所示。
步骤S9、根据潜在地质灾害可能性参数K的计算公式,超声扫描图像实现准备的标准值A=768*576,T=8.27,可得K的计算结果为2568.31÷(768×576)×738.29÷16.42×8.27≈2.16,此时S/N值为2568.31÷16.42=156.41,L值为738.29,即可判定钻孔内 孔壁上该区域的岩体结构相对较为破碎,该区域存在滑坡的风险,存在一定的潜在地质 灾害风险;钻孔内孔壁岩体内部结构剖面成像效果图,如图4所示。
实施例3
本发明实施例提供了一种岩体内部异常结构可视化监测与潜在地质灾害评估***, 包括超声成像剖面图获取设备、超声成像剖面图处理设备、图像特征进行统计模块以及 潜在地质灾害程度评估模块;
所述超声成像剖面图获取设备,用于获取目标区域内所有位置、角度的岩体结构超 声成像剖面图;
所述超声成像剖面图处理设备,用于对所述所有位置、角度的超声成像剖面图进行 组合及拼接融合处理,获得处理后的岩体结构超声成像剖面图;
所述图像特征进行统计模块,用于对所述处理后的岩体结构超声成像剖面图中二值 化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值 L;
所述潜在地质灾害程度评估模块,用于根据所述S、N、L获取潜在地质灾害可能 性参数,并根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度。
需要说明的是,所述超声成像剖面图获取设备具体包括B超仪、B超探头、推杆、 图像处理模块;
针对钻孔内部的岩体结构,超声扫描成像仪的探头可以在钻井平台或特定设备的帮 助下,借助于推杆或滑竿和夹子等辅助装置,把超声扫描成像探头固定在支架上并放入钻孔内部岩体孔壁的目标区域中,超声扫描成像探头可以实现上下左右移动和360度自 由旋转,其中,超声扫描成像探头固定方式示意图,如图5所示;
由于钻孔内部一般都有水或者已注水,或者注入了特定耦合剂,基于超声相控阵成 像原理的超声扫描成像探头可以实现岩体孔壁内部的超声扫描成像,其中,超声波的频率和相关配置的调整,以超声波可以穿透岩体内部结构1米左右并可以相对清晰的得到 超声扫描成像图为准;
针对超声扫描成像仪实时获取的大量超声扫描成像图,可以利用图像处理方法自动 识别超声扫描成像图中的异常结构,并进行定性定量分析,具体的,首先利用超声扫描成像仪(超声扫描成像探头)直接获取超声扫描成像视频*.cin文件,然后,利用高斯滤 波方法和开闭运算方法去除超声扫描成像图中的高频和低频噪声信号,滤除噪声,保留 岩体结构中结构面、破碎带等异常结构的超声图像,得到异常结构的二值化图像,即为 岩体结构超声成像剖面图;接着统计分析二值化图像中的超声光斑点数和边界条纹个数 及其所占的总面积;最后,通过监测这些光斑和边界条纹的变化情况,从而实现岩体内 部异常结构的可视化监测;
利用超声成像仪直接获取超声扫描成像视频*.cin文件;然后,利用高斯滤波方法和 开闭运算方法去除超声扫描图像中的高频和低频噪声信号,滤除噪声,保留岩体结构中结构面、破碎带等异常结构的超声图像,得到异常结构的二值化图像;接着统计分析二 值化图像中的超声光斑点数和边界条纹个数及其所占的总面积;最后,通过监测这些光 斑和边界条纹的变化情况,从而实现岩体内部异常结构的可视化监测;岩体内部结构的 可视化监测示意图,如图6所示;
优选的,所述超声成像剖面图获取设备包括超声扫描成像探头、推杆和图像处理模 块;超声成像剖面图获取设备获取目标区域内所有位置、角度的岩体结构超声成像剖面图,具体包括,
通过推杆将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过推杆旋转、 滑动或者侧向倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描 成像图,所述图像处理模块对所述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角度的岩体结构超声成像剖面图。
为了能够连续的并多方面的可视化监测岩体内部的异常结构,本发明实施例提供了 一种可使超声成像探头方便移动旋转的推杆。
该推杆具有以下特征:(1)每根推杆均为空心细不锈钢管,外径15mm,内径10mm,内部可窜带走线缆;(2)每根推杆均为2米,并从前端到后端依次标上深度刻度,并且 在后端标上360度方位角;(3)每根推杆后端均可与下一根推杆前端通过螺母连接;(4) 第一根推杆前端需连接超声成像探头,并可通过旋转螺母固定超声成像探头,实现探头 原地360度旋转,实现探头平移滑动和侧向倾斜。
通过推杆可实现超声成像探头在钻孔内部孔壁的平移、倾斜扫描成像和钻孔底部的 旋转扫描成像,三种超声扫描成像示意图,如图7(a)-(c)所示;
通过读取推杆的根数和最后一根推杆上刻度计算推杆的总长度,可知钻孔内超声成 像探头的准确位置;并读取最后一根推杆上旋转角度,可知钻孔内超声成像探头所对准的方向,据此可计算得出钻孔内超声成像探头的位置和方向等实时状态信息。通过探头 的位置和方向的实时信息来计算每个位置方向所获得的超声扫描图像对应于钻孔内壁 岩体内部的岩体结构信息。
根据推杆的总长度和旋转角度,得出孔内当前位置当前旋转角度的超声扫描图像, 进而得出该位置角度下的岩体结构超声成像剖面图,过一段时间后,旋转角度或移动到 下一个位置,根据该位置该旋转角度下的超声扫描图像,得出下一个位置角度的岩体结构超声成像剖面图,如此反复,依次得出目标区域内所有位置、角度下的岩体结构超声 成像剖面图;
通过旋转或平移超声扫描成像探头,获取多个位置、角度获得的超声扫描图像示意 图,如图8所示;针对所述处理后的岩体结构超声成像剖面图,统计分析每一个角度位置岩体内部结构超声成像剖面图中二值化图像特征;
优选的,所述潜在地质灾害程度评估模块根据所述S、N、L获取潜在地质灾害可 能性参数,具体包括,利用公式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的 总区域,T为修正系数,N≥3。
优选的,所述潜在地质灾害程度评估模块,根据所述地质灾害可能性参数评估岩体 内部异常结构潜在地质灾害程度,具体包括,所述地质灾害可能性参数越大,评定岩体内部异常结构潜在地质灾害程度越大。
需要说明的是,实施例1-3所述技术方案未重复描述之处可相互借鉴。
本发明公开了一种岩体内部异常结构监测与潜在地质灾害评估方法及***,通过获 取目标区域内所有位置、角度的岩体结构超声成像剖面图,对所述所有位置、角度的超声成像剖面图进行组合及拼接融合处理,获得处理后的岩体结构超声成像剖面图;对所 述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区域的总像 素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与 所有光斑块和条纹带所占像素点个数的比值L;根据所述S、N、L获取潜在地质灾害可 能性参数,并根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度; 实现了岩体内部异常结构的可视化监测与潜在地质灾害的准确评估。
本领域技术人员可以理解,实现上述实施例方法的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于计算机可读存储介质中。其 中,所述计算机可读存储介质为磁盘、光盘、只读存储记忆体或随机存储记忆体等。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本 发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保 护范围内。
Claims (6)
1.一种岩体内部异常结构监测与潜在地质灾害评估方法,其特征在于,包括以下步骤:
获取目标区域内所有位置、角度的岩体结构超声成像剖面图,对所述所有位置、角度的超声成像剖面图进行组合及拼接融合处理,获得处理后的岩体结构超声成像剖面图;
对所述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值L;
利用公式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的总区域,T为修正系数,T是超声成像剖面图的图像对比度来确定,图像对比度越大T值越小,N≥3,
根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度。
2.根据权利要求1所述的岩体内部异常结构监测与潜在地质灾害评估方法,其特征在于,所述获取目标区域内多个位置角度的岩体结构超声成像剖面图,具体包括,
将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过旋转、滑动或者侧向倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描成像图,对所述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角度的岩体结构超声成像剖面图。
3.根据权利要求2所述的岩体内部异常结构监测与潜在地质灾害评估方法,其特征在于,根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度,具体包括,所述地质灾害可能性参数越大,评定岩体内部异常结构潜在地质灾害程度越大。
4.一种岩体内部异常结构监测与潜在地质灾害评估***,其特征在于,包括超声成像剖面图获取设备、超声成像剖面图处理设备、图像特征进行统计模块以及潜在地质灾害程度评估模块;
所述超声成像剖面图获取设备,用于获取目标区域内所有位置、角度的岩体结构超声成像剖面图;
所述超声成像剖面图处理设备,用于对所述所有位置、角度的超声成像剖面图进行组合及拼接融合处理,获得处理后的岩体结构超声成像剖面图;
所述图像特征进行统计模块,用于对所述处理后的岩体结构超声成像剖面图中二值化图像特征进行统计,得到非零区域的总像素点数S、非零不连续光斑块和条纹带的总个数N、所有光斑块和条纹带的边界长度与所有光斑块和条纹带所占像素点个数的比值L;
所述潜在地质灾害程度评估模块,用于利用公式
获取潜在地质灾害可能性参数K,其中,A为超声扫描成像探头所探测覆盖的总区域,T为修正系数,T是超声成像剖面图的图像对比度来确定,图像对比度越大T值越小,N≥3,
根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度。
5.根据权利要求4所述的岩体内部异常结构监测与潜在地质灾害评估***,其特征在于,所述超声成像剖面图获取设备包括超声扫描成像探头、推杆和图像处理模块;超声成像剖面图获取设备获取目标区域内所有位置、角度的岩体结构超声成像剖面图,具体包括,
通过推杆将超声扫描成像探头入钻孔内部岩体孔壁的目标区域中,通过推杆旋转、滑动或者侧向倾斜所述超声扫描成像探头,获取目标区域内所有位置、角度的超声扫描成像图,所述图像处理模块对所述所有位置、角度的超声扫描成像图进行高频和低频噪声信号滤除,得到所有位置、角度的岩体结构超声成像剖面图。
6.根据权利要求5所述的岩体内部异常结构监测与潜在地质灾害评估***,其特征在于,所述潜在地质灾害程度评估模块,根据所述地质灾害可能性参数评估岩体内部异常结构潜在地质灾害程度,具体包括,所述地质灾害可能性参数越大,评定岩体内部异常结构潜在地质灾害程度越大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911226383.3A CN110940734B (zh) | 2019-12-04 | 2019-12-04 | 岩体内部异常结构监测与潜在地质灾害评估方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911226383.3A CN110940734B (zh) | 2019-12-04 | 2019-12-04 | 岩体内部异常结构监测与潜在地质灾害评估方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110940734A CN110940734A (zh) | 2020-03-31 |
CN110940734B true CN110940734B (zh) | 2021-03-30 |
Family
ID=69909837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911226383.3A Expired - Fee Related CN110940734B (zh) | 2019-12-04 | 2019-12-04 | 岩体内部异常结构监测与潜在地质灾害评估方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110940734B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101486033A (zh) * | 2009-03-09 | 2009-07-22 | 中国农业大学 | 一种果品分检装置及板栗分检方法 |
CN101555661A (zh) * | 2009-04-14 | 2009-10-14 | 山东农业大学 | 基于机器视觉技术的棉花异性纤维重量计量方法及*** |
CN102037492A (zh) * | 2008-05-23 | 2011-04-27 | 澳大利亚国立大学 | 图像数据处理 |
CN102636656A (zh) * | 2012-04-01 | 2012-08-15 | 长春迪瑞医疗科技股份有限公司 | 一种全自动尿有形成分分析仪的校准方法 |
CN102798583A (zh) * | 2012-07-13 | 2012-11-28 | 长安大学 | 一种基于改进的ferret的矿岩块度测量方法 |
CN103679675A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种面向水质定量遥感应用的遥感影像融合方法 |
CN103698397A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油化工股份有限公司 | 一种定量接触压力超声波检测***及其检测方法 |
CN105092589A (zh) * | 2015-07-07 | 2015-11-25 | 东北大学 | 一种胶囊头缺陷检测方法 |
CN107024416A (zh) * | 2017-04-07 | 2017-08-08 | 福建福光股份有限公司 | 结合相似性和不连续性的准圆形颗粒平均尺寸检测方法 |
CN206740654U (zh) * | 2017-05-15 | 2017-12-12 | 中国地质调查局油气资源调查中心 | 一种全直径页岩样品的快速综合评价装置 |
CN109614851A (zh) * | 2018-10-26 | 2019-04-12 | 上海市环境监测中心(上海长三角区域空气质量预测预报中心) | 大气污染颗粒物形态参数的图像检测方法 |
CN109931876A (zh) * | 2019-03-11 | 2019-06-25 | 三峡大学 | 一种通过摄影测量土工试样形变的方法 |
CN109946380A (zh) * | 2019-04-03 | 2019-06-28 | 石家庄铁道大学 | 无砟轨道砂浆层脱空检测方法、装置及终端设备 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO345502B1 (no) * | 2012-03-28 | 2021-03-08 | Logined Bv | Omforming av seismisk attributtfargemodell |
CN103325118A (zh) * | 2013-06-26 | 2013-09-25 | 中国石油大学(北京) | 一种获取碳酸盐岩岩心孔洞特征参数的方法及装置 |
CN106053475B (zh) * | 2016-05-24 | 2018-10-23 | 浙江工业大学 | 基于主动式全景视觉的隧道病害全断面动态快速检测装置 |
CN107064169B (zh) * | 2017-01-22 | 2021-11-09 | 太原科技大学 | 一种桥梁裂缝检测装置及检测方法 |
-
2019
- 2019-12-04 CN CN201911226383.3A patent/CN110940734B/zh not_active Expired - Fee Related
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102037492A (zh) * | 2008-05-23 | 2011-04-27 | 澳大利亚国立大学 | 图像数据处理 |
CN101486033A (zh) * | 2009-03-09 | 2009-07-22 | 中国农业大学 | 一种果品分检装置及板栗分检方法 |
CN101555661A (zh) * | 2009-04-14 | 2009-10-14 | 山东农业大学 | 基于机器视觉技术的棉花异性纤维重量计量方法及*** |
CN102636656A (zh) * | 2012-04-01 | 2012-08-15 | 长春迪瑞医疗科技股份有限公司 | 一种全自动尿有形成分分析仪的校准方法 |
CN102798583A (zh) * | 2012-07-13 | 2012-11-28 | 长安大学 | 一种基于改进的ferret的矿岩块度测量方法 |
CN103698397A (zh) * | 2012-09-27 | 2014-04-02 | 中国石油化工股份有限公司 | 一种定量接触压力超声波检测***及其检测方法 |
CN103679675A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种面向水质定量遥感应用的遥感影像融合方法 |
CN105092589A (zh) * | 2015-07-07 | 2015-11-25 | 东北大学 | 一种胶囊头缺陷检测方法 |
CN107024416A (zh) * | 2017-04-07 | 2017-08-08 | 福建福光股份有限公司 | 结合相似性和不连续性的准圆形颗粒平均尺寸检测方法 |
CN206740654U (zh) * | 2017-05-15 | 2017-12-12 | 中国地质调查局油气资源调查中心 | 一种全直径页岩样品的快速综合评价装置 |
CN109614851A (zh) * | 2018-10-26 | 2019-04-12 | 上海市环境监测中心(上海长三角区域空气质量预测预报中心) | 大气污染颗粒物形态参数的图像检测方法 |
CN109931876A (zh) * | 2019-03-11 | 2019-06-25 | 三峡大学 | 一种通过摄影测量土工试样形变的方法 |
CN109946380A (zh) * | 2019-04-03 | 2019-06-28 | 石家庄铁道大学 | 无砟轨道砂浆层脱空检测方法、装置及终端设备 |
Non-Patent Citations (2)
Title |
---|
基于成像测井技术的地下岩体裂隙分布规律预测研究;涂善波等;《河南水利与南水北调》;20190131;第83-85 * |
基于数字照相量测的岩体裂隙识别与表征初步研究;杨苏;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20180315;C038-1051 * |
Also Published As
Publication number | Publication date |
---|---|
CN110940734A (zh) | 2020-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2981205C (en) | Methods and means for creating three-dimensional borehole image data | |
CN102680575B (zh) | 一种复杂岩土介质的冲击映像方法及*** | |
BRPI0413387B1 (pt) | método para visualização e identificação de materiais em um conduto de transporte de fluidos | |
Leucci et al. | GPR and sonic tomography for structural restoration: the case of the cathedral of Tricarico | |
US7295928B2 (en) | Enhancing the quality and resolution of an image generated from single or multiple sources | |
CN109521479A (zh) | 一种地下连续墙渗漏水的非开挖检测方法 | |
JP2008500470A (ja) | アンカーボアホールを調査するための装置 | |
Sturzenegger et al. | Application and limitations of ground-based laser scanning in rock slope characterization | |
CN112965135B (zh) | 一种石窟崖体裂隙空间异质分布的无损探测综合方法 | |
JP4936387B2 (ja) | コンクリート構造物の劣化診断方法 | |
US20130083888A1 (en) | Apparatus for detecting volume of foreign substance existed in core of geological sample using computer tomography apparatus and method thereof | |
CN108776175A (zh) | 一种冻结壁平均抗压强度超声波检测方法 | |
CN108508093A (zh) | 一种工件缺陷高度的检测方法及*** | |
CN110940734B (zh) | 岩体内部异常结构监测与潜在地质灾害评估方法及*** | |
Simonin et al. | Case study of detection of artificial defects in an experimental pavement structure using 3D GPR systems | |
Wang et al. | Fine detection technology of rock mass structure based on borehole acousto-optic combined measurement | |
CN116519477A (zh) | 一种岩石失稳判别方法、装置及存储介质 | |
Van Knapen et al. | Identification and characterisation of rock mass discontinuity sets using 3D laser scanning | |
US3474879A (en) | Acoustic method for mapping the surface characteristics of a borehole | |
Broding | Volumetric scanning well logging | |
Zou et al. | The visual detection of inner abnormal structures and potential geological hazards assessment method for buried rock and soil mass | |
JP2020041405A (ja) | 凍結地盤の評価方法、及び凍結地盤の評価装置 | |
Karlöf et al. | Statistical techniques to select detection thresholds for peak signals in ice-core data | |
CN116183875B (zh) | 岩体波速的确定方法、装置和电子设备 | |
Song et al. | Morphological Reconstruction for the Internal Structure of In-Situ SoilRock Mixture Through Ultrasonic Transmission Imaging Experiments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210330 Termination date: 20211204 |