基于边缘阈值分割的钻孔图像RQD自动分析方法研究

李清波; 杜朋召

doi:10.11779/CJGE202011022

基于边缘阈值分割的钻孔图像RQD自动分析方法研究 English Version

黄河勘测规划设计研究院有限公司,河南　郑州　450003

基金项目:

国家重点研发计划 2018YFC1508703

详细信息

作者简介:
李清波(1964—),男,教授级高级工程师,主要从事工程地质及水文地质勘察与研究工作。E-mail：liqb_yrec@163.com

中图分类号: TU45
计量
- 文章访问数: 370
- HTML全文浏览量: 33
- PDF下载量: 179
出版历程
- 收稿日期: 2020-06-17
- 网络出版日期: 2022-12-05
- 刊出日期: 2020-10-31

Automatic RQD analysis method based on information recognition of borehole images

Yellow River Engineering Consulting Co., Ltd., Zhengzhou 450003, China

摘要

摘要: RQD是评价岩体完整性的重要指标,传统获取方法受钻进工艺、取芯质量和回次长度影响较大,不能客观反映岩体质量。针对该问题,提出了一种基于边缘阈值分割的钻孔图像RQD自动分析方法。该方法首先对钻孔图像进行预处理和边缘阈值分割,实现图像中结构面和破碎带目标筛分；然后,对筛分出的目标区域进行编码,并进行连通区合并分析,确定影响RQD分析的有效目标；最后,提取有效目标中结构面和破碎带的位置、宽度信息,计算钻孔岩体RQD值。以泾河东庄水利枢纽工程坝址区钻孔为例,进行图像RQD自动分析,结果与钻孔图像、波速测试结果一致性较好,提高了钻孔RQD统计的准确性。该方法丰富了RQD获取途径,为评价岩体完整性提供了一种快速有效的方法。
- 钻孔图像 /
- RQD /
- 边缘信息 /
- 阈值分割 /
- 自动分析
Abstract: RQD is an important index for evaluating the integrity of rock masses. The traditional methods are greatly affected by drilling process, core quality and run length, and cannot objectively reflect the quality of rock masses. To solve this problem, An automatic RQD analysis method based on the edge threshold segmentation of borehole images is proposed. This method first performs preprocessing and edge threshold segmentation on the borehole images to achieve target screening. Then, encode the screened target area is encoded, and the connected areas are merged to determine the effective targets that affect the RQD analysis. Finally, the location and width of the structural plane and fracture zone are extracted, and the RQD of boreholes is calculated. For a case study of the Dongzhuang Dam at Jinghe River, the automatic RQD analysis is performed, and the achieved results are in good agreement with the borehole images and the coefficient of rock mass integrity. The proposed method improves the accuracy of RQD statistics, enriches the acquisition of RQD, and provides a fast and effective method for evaluating the rock mass integrity.
- borehole image /
- RQD /
- edge information /
- threshold segmentation /
- automatic analysis

HTML全文

0. 引言

在建筑工程领域中,对于水平荷载作用下桩基的受力分析,常见有极限地基反力法、弹性地基反力法、复合地基反力法、弹性理论法和p-y曲线法等^[1-2]。弹性地基反力法又包括地基系数常数法、k法、c法、m法以及吴恒立^[3]的双参数法。张有龄给出了地基系数为常数时的桩身响应解析解,N.B.ypdh与众多学者给出了桩身内力与变形的幂级数解,更有采用^[4]纽玛克法、有限差分法与有限元法来求解桩身内力与变形。

上述常用方法中对于多层地基情况的处理略显粗糙,如目前建筑桩基^[6]与公路桥涵桩基领域^[7]最常采用各层地基按其地基系数以权重进行折算,得到一个地基系数的等效值。近年来,Pise^[5]对双层地基水平受荷桩进行了数值求解,赵明华等^[8-10]对成层地基中桩的受力与试验做了大量工作,并尝试用无网格法分析计算,戴自航等^[11]采用有限元与有限差分进行数值计算,竺明星等^[12]利用矩阵传递法依次求解多层地基中的桩身各点内力,詹红志等^[13]也采用类似矩阵传递方法对抗滑桩嵌固段多层岩层进行了计算。

本文不同于先前学者从桩身形函数利用幂级数角度出发,引入张氏法的解析解函数形式,利用节点内力变形连续条件,建立全桩全节点统一矩阵线性方程,引入边界条件后一次性求解所有节点的变位与内力,并将该计算方法应用于多层地基桩基的水平响应计算。

1. 线性方程解力学模型

1.1 力学模型建立

不同于竺明星等^{[12, 14]}建立三参数地基系数模型并利用Laplace变换求解桩身响应的方法,本文在理论推导过程中不特别假定地基系数的分布模式,但考虑到设计人员使用上的便利性,以单层地基m值与多层地基m值分别演示计算过程。

根据Winkler理论,假定地基是服从胡克定律的弹性体,且每层地基厚度为h_j,如图1所示。

图 1 线性方程力学模型

Figure 1. Linear equation mechanics model

下载: 全尺寸图片幻灯片

将桩身沿深度方向分成 $n$ 段,桩单元依次编号为1,2,···, $n$ ,桩结点编号为0,1,2,···, $n$ ,结点对应各自坐标值。记桩身水平位移为 $y (z)$ ,桩身转角为 $φ (z)$ ,桩身弯矩为 $M (z)$ ,桩身剪力为 $F (z)$ , $M_{0}, F_{0}$ 表示桩顶作用的弯矩与水平力, $M_{i, z_{j}}, F_{i, z_{j}}$ 表示第i桩单元的 $z_{j}$ 节点处的弯矩与剪力,对于第 $i$ 桩单元,第 $i$ 段内地基系数 $k_{i}$ 以该段内的积分中值定理为原则,即

$k_{i} = \frac{\int_{z_{i - 1}}^{z_{i}} k (z) d z}{z_{i} - z_{i - 1}}$ 。

(1)

约定弯矩以桩左侧受拉为正,剪力以使桩顺时针转动方向为正,水平位移以坐标正向为正,而截面转角以逆时针转动为正。

1.2 微分方程及解答

对于第 $i$ 段,满足如下微分方程：

$E I \frac{d^{4} y_{i}}{d z^{4}} + k_{i} b y_{i} = 0$ 。

(2)

令 $A_{i} = \sqrt[4]{\frac{k_{i} b}{4 E I}}$ ,则可得到第 $i$ 段 $z \in [z_{i}, z_{i - 1}]$ 的挠曲线方程 $y_{i} (z)$ 解析解与挠曲线各阶导数：

$\begin{matrix} y_{i} (z) = C_{i}_{1} e^{- A_{i} z} \sin (A_{i} z) + C_{i}_{2} e^{- A_{i} z} \cos (A_{i} z) + \\ C_{i}_{3} e^{A_{i} z} \sin (A_{i} z) + C_{i 4} e^{A_{i} z} \cos (A_{i} z) \end{matrix}$ 。

(3)

将以上各函数表达式整理成矩阵形式,如式（4）：

$\begin{matrix} y_{i} = [\begin{matrix} f_{i 1} (z) & f_{i 2} (z) & f_{i 3} (z) & f_{i 4} (z) \end{matrix}] \cdot \\ {[\begin{matrix} C_{i 1} & C_{i 2} & C_{i 3} & C_{i 4} \end{matrix}]}^{T}, \\ y^{(1)}_{i} = [\begin{matrix} g_{i 1} (z) & g_{i 2} (z) & g_{i 3} (z) & g_{i 4} (z) \end{matrix}] \cdot \\ {[\begin{matrix} C_{i 1} & C_{i 2} & C_{i 3} & C_{i 4} \end{matrix}]}^{T}, \\ y^{(2)}_{i} = [\begin{matrix} p_{i 1} (z) & p_{i 2} (z) & p_{i 3} (z) & p_{i 4} (z) \end{matrix}] \cdot \\ {[\begin{matrix} C_{i 1} & C_{i 2} & C_{i 3} & C_{i 4} \end{matrix}]}^{T}, \\ y^{(3)}_{i} = [\begin{matrix} q_{i 1} (z) & q_{i 2} (z) & q_{i 3} (z) & q_{i 4} (z) \end{matrix}] \cdot \\ {[\begin{matrix} C_{i 1} & C_{i 2} & C_{i 3} & C_{i 4} \end{matrix}]}^{T} 。 \end{matrix}}$

(4)

令各系数矩阵表达式如下：

$\begin{array}{l} [P^{*}_{i, z_{i - 1}}] = [\begin{matrix} E I p_{i 1, z_{i - 1}} & E I p_{i 2, z_{i - 1}} & E I p_{i 3, z_{i - 1}} & E I p_{i 4, z_{i - 1}} \end{matrix}], \\ [Q^{*}_{i, z_{i - 1}}] = [\begin{matrix} E I q_{i 1, z_{i - 1}} & E I q_{i 2, z_{i - 1}} & E I q_{i 3, z_{i - 1}} & E I q_{i 4, z_{i - 1}} \end{matrix}], \\ [f^{*}_{i, z_{i - 1}}] = [\begin{matrix} f_{i 1, z_{i - 1}} & f_{i 2, z_{i - 1}} & f_{i 3, z_{i - 1}} & f_{i 4, z_{i - 1}} \end{matrix}], \\ [g^{*}_{i, z_{i - 1}}] = [\begin{matrix} g_{i 1, z_{i - 1}} & g_{i 2, z_{i - 1}} & g_{i 3, z_{i - 1}} & g_{i 4, z_{i - 1}} \end{matrix}] 。 \end{array}}$

(5)

则桩身 $n$ 段各个节点处的弯矩、剪力、挠度与转角记成： $[B] = [V^{*}] \cdot [C^{*}]$ ,其中

$[B] = [\begin{matrix} M_{1, 0} & ... & M_{n, n - 1} \\ M_{1, 1} & ... & M_{n, n} \\ F_{1, 0} & ... & F_{n, n - 1} \\ F_{1, 1} & ... & F_{n, n} \\ y_{1, 0} & ... & y_{n, n - 1} \\ y_{1, 1} & ... & y_{n, n} \\ φ_{1, 0} & ... & φ_{n, n - 1} \\ φ_{1, 1} & ... & φ_{n, n} \end{matrix}]$ ,

(6)

(7)

$[C^{*}] = {[\begin{matrix} [C_{1, j}] & [0] & ... & [0] \\ [0] & [C_{2, j}] & ... & [0] \\ [0] & [0] & ... & [0] \\ [0] & [0] & ... & [C_{n, j}] \end{matrix}]}_{j = 1, 2, 3, 4}$ 。

(8)

1.3 全桩线性方程组的构建与解答

保证桩身每一结点处内力与位移是连续的,以此思路建立桩身全结点的线性方程组如下：

$\begin{array}{l} M_{i, z_{i}} = M_{i + 1, z_{i}}, \\ F_{i, z_{i}} = F_{i + 1, z_{i}}, \\ y_{i, z_{i}} = y_{i + 1, z_{i}}, \\ φ_{i, z_{i}} = φ_{i + 1, z_{i}} 。 \end{array}}$

(9)

最终记成如下矩阵形式： $[ξ^{*}] \cdot [C] = [H]$ ,其中

$[C] = {[\begin{array}{l} \begin{matrix} [C_{1, j}] \\ [C_{2, j}] \\ [C_{3, j}] \\ ... \\ [C_{n - 1}_{, j}] \end{matrix} \\ [C_{n}_{, j}] \end{array}]}_{j = 1, 2, 3, 4}, [H] = [\begin{matrix} M_{0} \\ F_{0} \\ 0 \\ ... \\ 0 \\ H_{n - 1, z_{n}} \\ H_{n, z_{n}} \end{matrix}],$

(10)

$[ξ^{*}] = [\begin{matrix} [P^{*}_{1, z_{0}}] & [0] & [0] \\ [Q^{*}_{1, z_{0}}] & [0] & [0] \\ [P^{*}_{1, z_{1}}] & - [P^{*}_{2, z_{1}}] & [0] \\ [Q^{*}_{1, z_{1}}] & - [Q^{*}_{2, z_{1}}] & [0] \\ [f^{*}_{1, z_{1}}] & - [f^{*}_{2, z_{1}}] & [0] \\ [g^{*}_{1, z_{1}}] & - [g^{*}_{2, z_{1}}] & [0] \\ [0] & ... & ... \\ [0] & [0] & [ξ^{*}_{4 n - 1, n}] \\ [0] & [0] & [ξ^{*}_{4 n, n}] \end{matrix}]$ 。

(11)

上式矩阵运算表示了全桩全结点内力与位移值需要满足该线性方程组,引入桩顶与桩端的边界条件后,等式右侧矩阵也为常数阵,这样可通过Gauss消元等多种方法求解线性方程组,解得桩身每一段的四组参数 $C_{i 1}, C_{i 2}, C_{i 3}, C_{i 4}$ ,再将系数C矩阵回代式（8）即可。

2. 基于m法的设计计算与结果验证

2.1 单层地基算例验证

某建筑物^[2]采用桩基基础,直径d=1.5 m,埋入并支持在非岩石类土中,入土深度h=15 m,桩头在地面处自由,作用有水平荷载H₀=60 kN和M₀=700 kN·m,C25级混凝土的弹性模量E_c=2.8×10⁴ MPa= 2.8×10⁷ kN/m²,地基的反力系数的比例系数m=9400 kN/m⁴,土的内摩擦角 $φ$ =22°,黏聚力 $c = 15 {kN/m}^{2}$ ,重度 $γ = 20 k N / m^{3}$ 。

$b = K_{φ} K_{0} d = 0.9 (1.5 + 1) = 2.25 m$ ,

$E I = \frac{0.85 \times 2.8 \times 10^{7} \times π \times {1.5}^{4}}{64} = 59.3 \times 10^{5} kN \cdot m^{2}$ 。

分别将n取5,10,15,20进行了桩身内力与位移计算,本文法计算结果与传统m法计算的桩身弯矩绘制成曲线图,如图2所示。按规范法计算桩身最大弯矩为766.9 kN·m,将桩等分20段后桩身弯矩最大值为762.8 kN·m,相比规范法误差0.53%。从上图2看出,桩身弯矩随深度增加总体呈现先上升后下降的过程,弯矩极值出现在距离桩顶(1～3)d范围之间（d为桩身直径）,弯矩零点位于距离桩顶6d位置左右。

图 2 不同分段数桩身弯矩

Figure 2. Moments of pile body with different numbers of sections

下载: 全尺寸图片幻灯片

图3,4中看出3种端部约束下的弯矩、位移曲线重合度较高,仅在桩端附近处弯矩曲线出现了分叉发展的趋势。位移零点出现在距离桩顶4d位置处,相比桩身弯矩的6d变化范围缩小了33.3%。本算例所得的桩身弯矩极值、弯矩与位移零点所在的桩身位置符合目前国内外学者的研究结果,如赵明华等^[15-16]曾建议桩影响范围取3～5d,冯忠居等^[17]建议取（2～8）d等。

图 3 不同桩端约束的桩身弯矩

Figure 3. Moments of pile body restrained by different pile ends

下载: 全尺寸图片幻灯片

图 4 不同桩端约束的桩身位移

Figure 4. Displacements of pile body with different pile end constraints

下载: 全尺寸图片幻灯片

2.2 多层地基算例验证

某圆形^{[12, 18]}截面灌注桩^[10]桩径d=1.0 m,地面处桩顶剪力Q =150 kN,弯矩M =0,桩的弹性模量E =2.1675×10 kN/m²。桩侧有两层地基土体：第一层为流塑状回填土,层厚为2.0 m,相应的地基反力系数m为3000 kN/m⁴；第二层为硬塑状黏性土,桩身在该层土体中的长度为10.0 m,相应的地基反力系数为20000 kN/m⁴。

表1为不同计算方法的计算结果,从中可知本文线性方程解与精确解之间的桩顶位移误差为2.3%,最大弯矩误差为0.285%。

表 1 不同计算方法结果对比

Table 1. Results of different calculation methods

计算方法	桩顶位移/mm	最大弯矩/(kN·m)	最大弯矩位置/m
精确解^[18]	4.3735	336.81	3
规范解^[6]	3.0156	238.95	3
挠度加权换算^[9]	4.2093	343.83	3
有限差分法^[11]	4.2989	337.60	3
杆系有限元法^[11]	4.2520	336.07	3
矩阵传递法^[12]	4.2990	337.60	3
本文法	4.2703	335.85	3

下载: 导出CSV

| 显示表格

笔者在本算例基础上,改变地层情况再次进行桩身弯矩与位移计算,分别将地层视为全为上层土的单一地层与全为下层土的单一地层（简称“上层土地基”与“下层土地基”）,将计算结果分别绘制成曲线图5,6用以对比分析。由图5可以看出,桩身弯矩极值出现在距离桩顶（1～5）d之间,本例中双层地基情况与“下层土地基”情况都在8d位置处达到了弯矩零点；“下层土地基”与“两层土地基”均在深度5d处为位移零点。

图 5 不同地层情况桩身弯矩

Figure 5. Moments of pile body in different strata

下载: 全尺寸图片幻灯片

图 6 不同地层情况桩身位移

Figure 6. Pile displacements in different strata

下载: 全尺寸图片幻灯片

3. 结论

假设地基为弹性材料,分段建立梁挠曲线微分方程,通过结点内力与位移的连续条件一次性建立桩身全结点的线性方程组,求解得到各点内力与位移。以两个算例验证了线性方程解法在单层地基与多层地基中桩身响应计算的正确性,并对桩底不同边界条件、桩身周围不同地层进行了计算与讨论,得出如下结论：

（1）在桩顶水平荷载的作用下,桩身弯矩最大值出现在距离桩顶（1～5）d范围内,桩顶附近土层抗力越差,最大弯矩所出现的位置将越深。

（2）在距离桩顶（6～8）d位置附近将出现弯矩函数零点,且下降段所处区间受桩中部土层的抗力大小控制。

（3）桩身位移最大值出现在桩顶,距桩顶5d位置处出现位移零点。桩端不同的边界条件对桩身的位移影响较小,而桩周土层的抗力大小对桩身的位移起到控制作用。

（4）同一种情况下的桩身水平响应,其弯矩零点所出现的位置将比位移零点所出现的位置滞后（2～3）d。

图 1 钻孔二维平面图像和三维柱状图像

Figure 1. 2D planar and 3D columnar images of borehole

下载: 全尺寸图片幻灯片

图 2 钻孔结构面图像和破碎岩体图像

Figure 2. Structural plane and broken rock mass images of borehole

下载: 全尺寸图片幻灯片

图 3 分析流程图

Figure 3. Flow chart of proposed analysis method in this study

下载: 全尺寸图片幻灯片

图 4 图像梯度直方图

Figure 4. Histogram of image gradient

下载: 全尺寸图片幻灯片

图 5 图像分割结果

Figure 5. Results of image segmentation

下载: 全尺寸图片幻灯片

图 6 连通区合并结果

Figure 6. Results of merging connected regions

下载: 全尺寸图片幻灯片

图 7 ZK408钻孔原始图像

Figure 7. Original image in ZK408

下载: 全尺寸图片幻灯片

图 8 钻孔RQD实测值与计算值对比图

Figure 8. Comparison between measured and calculated values of RQD

下载: 全尺寸图片幻灯片

图 9 钻孔16.5~17.3 m段孔壁图像与岩芯图像

Figure 9. Image of borehole rock mass and core in 16.5 ~ 17.3 m

下载: 全尺寸图片幻灯片

图 10 钻孔RQD与完整性系数对比图

Figure 10. Comparison between RQD and integrity coefficient

下载: 全尺寸图片幻灯片

表 1 钻孔图像分割参数表

Table 1 Parameters for image segmentation

图像深度/m	梯度均值	梯度标准差	梯度点	灰度均值
8~10	0.1045	0.0309	0.1554	142.77
30~32	0.1285	0.0911	0.2196	124.52

下载: 导出CSV

表 2 岩体图像信息识别表

Table 2 Recognition of image information of rock mass

图像深度/m	连通区编号	像素数目	最小行号	最大行号	所占列数
8~10	1	338	45	117	70
	2	117	411	429	45
	3	280	545	576	67
30~32	1	3455	335	524	68
30~32	2	134	588	600	42

下载: 导出CSV

表 3 钻孔RQD分析结果

Table 3 Results of RQD analysis

图像编号	编号	起深度/m	止深度/m	中心线深度/m	宽度/m	RQD
图7（a）	结构面1	8.15	8.39	8.27	0.017	96.5
	结构面2	9.37	9.43	9.40	0.022
	结构面3	9.82	9.92	9.87	0.031
图7（b）	破碎带1	31.12	31.75	31.43	0.630	66.5
图7（b）	结构面4	31.96	32.00	31.98	0.016	66.5

下载: 导出CSV

参考文献(18)

[1]	DEERE D U. Technical description of rock cores for engineering purposes[J]. Rock Mechanics and Engineering Geology, 1963, 1(1): 18-22.
[2]	PALMSTROM A. Measurements of and correlations between blocksize and rock quality designation (RQD)[J]. Tunn Undergr SpaceTechnol, 2005, 20: 362-377. doi: 10.1016/j.tust.2005.01.005
[3]	OZTURK C A, NASUF E. Geostatistical assessment of rock zones for tunneling[J]. Tunneling and Underground Space Technology, 2002(17): 275-285.
[4]	王国欣, 肖树芳, 陈剑平, 等. 不连续面三维网络在RQD中的应用研究[J]. 岩石力学与工程学报, 2002, 21(12): 1761-1764. doi: 10.3321/j.issn:1000-6915.2002.12.003 WANG Guo-xin, XIAO Shu-fang, CHEN Jian-ping, et al. Study on application of 3D network of discontinuities to RQD[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(12): 1761-1764. (in Chinese) doi: 10.3321/j.issn:1000-6915.2002.12.003
[5]	陈剑平, 王清, 赵红亮, 等. 窗口测线法获取岩体RQD[J]. 岩石力学与工程学报, 2004, 23(9): 1491-1495. doi: 10.3321/j.issn:1000-6915.2004.09.016 CHEN Jian-ping, WANG Qing, ZHAO Hong-liang, et al. Obtaining RQD of rock mass by sampling window method[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(9): 1491-1495. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.09.016
[6]	徐伟, 胡新丽, 黄磊, 等. 结构面三维网络模拟计算RQD及精度对比研究[J]. 岩石力学与工程学报, 2012, 31(4): 822-833. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201204024.htm XU Wei, HU Xin-li, HUANG Lei, et al. Research on RQD of rock mass calculated by three-dimensional discontinuity network simulation method and its accuracy comparison[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(4): 822-833. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201204024.htm
[7]	崔小平, 秦趣, 马士彬. 三维激光扫描技术获取岩体RQD值中的应用[J]. 六盘水师范高等专科学校学报, 2011, 23(3): 29-32. doi: 10.3969/j.issn.1671-055X.2011.03.007 CUI Xiao-ping, QIN Qu, MA Shi-bin. The application of 3d laser scanning technology in rock RQD data[J]. Journal of Liupanshui Teachers College, 2011, 23(3): 29-32. (in Chinese) doi: 10.3969/j.issn.1671-055X.2011.03.007
[8]	KAYABASI A, YESILOGLU-GULTEKIN GOKCEOGLU C. Use of non-linear prediction tools to assess rock masspermeability using various discontinuity parameters[J]. Engineering Geology, 2015, 185(5): 1-9.
[9]	Azimia , Abdolazim . A New method for improving the RQD determination of rock core in borehole[J]. Rock Mechanics & Rock Engineering, 2016, 49(4): 1559-1566.
[10]	胡瀚. 基于数字近景摄影测量的岩体RQD获取[D]. 长春: 吉林大学, 2018. HU Han. Rock Mass RQD Acquisition Based on Digital Close-range Photogrammetry[D]. Changchun: Jinlin University, 2018. (in Chinese)
[11]	THAPA B B, HUGHETT P, KARASAKI K. Semiautomatic analysis of rock fracture orientations from borehole wall images[J]. Geophysics, 1997, 62(1): 129-137. doi: 10.1190/1.1444112
[12]	吴剑, 冯少孔, 李宏阶. 钻孔成像中结构面自动判读技术研究[J]. 岩土力学, 2011, 32(3): 951-957. doi: 10.3969/j.issn.1000-7598.2011.03.051 WU Jian, FENG Shao-kong, LI Hong-jie. Study of automatically extracting structural plane parameters from borehole images[J]. Rock and Soil Mechanics, 2011, 32(3): 951-957. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.03.051
[13]	曹洋兵, 晏鄂川, 胡德新, 等. 岩体结构面产状测量的钻孔摄像技术及其可靠性[J]. 地球科学(中国地质大学学报), 2014, 39(4): 473-480. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201404010.htm CAO Yang-bing, YAN E-chuan, HU De-xin, et al. Calculation methods of rock mass discontinuity orientation measured by borehole camera technology and technology reliability[J]. Journal of China University of Geosciences: Earth Science, 2014, 39(4): 473-480. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201404010.htm
[14]	王川婴, 邹先坚, 韩增强, 等. 基于特征函数的钻孔图像结构面识别方法[J]. 岩土力学, 2017, 38(10): 3062-3066. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710038.htm WANG Chuan-ying, ZOU Xian-jian, HAN Zeng-qiang, et al. Recognition of structural planes in borehole image based on characteristic function[J]. Rock and Soil Mechanics, 2017, 38(10): 3062-3066. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710038.htm
[15]	葛云峰, 钟鹏, 唐辉明, 等. 基于钻孔图像的岩体结构面几何信息智能测量[J]. 岩土力学, 2019, 40(11): 4467-4476. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201911038.htm Intelligent measurement on geometric information of rock discontinuities based on borehole image[J]. Rock and Soil Mechanics, 2019, 40(11): 4467-4476. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201911038.htm
[16]	王川婴, Tim LAW K. 钻孔摄像技术的发展与现状[J]. 岩石力学与工程学报, 2005, 24(19): 3440-3448. doi: 10.3321/j.issn:1000-6915.2005.19.006 WANG Chuan-ying, Tim LAW K. Review of borehole camera technology[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(19): 3440-3448. (in Chinese) doi: 10.3321/j.issn:1000-6915.2005.19.006
[17]	汪进超, 王川婴, 胡胜, 等. 孔壁钻孔图像的结构面参数提取方法研究[J]. 岩土力学, 2017, 38(10): 3074-3080. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710040.htm WANG Jin-chao, WANG Chuan-ying, HU Sheng. A new method for extraction of parameters of structural surface in borehole images[J]. Rock and Soil Mechanics, 2017, 38(10): 3074-3080. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710040.htm
[18]	赵芳, 栾晓明, 孙越. 数字图像几种边缘检测算子检测比较分析[J]. 计算机应用, 2004, 24(9): 28-30. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDHJ200903024.htm ZHAO Fang, LUAN Xiao-ming, SUN Yue. Edge detection operators in digital image processing[J]. Computer Application, 2004, 24(9): 28-30. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDHJ200903024.htm