CN102419457B - 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法 - Google Patents

一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法 Download PDF

Info

Publication number
CN102419457B
CN102419457B CN 201110216593 CN201110216593A CN102419457B CN 102419457 B CN102419457 B CN 102419457B CN 201110216593 CN201110216593 CN 201110216593 CN 201110216593 A CN201110216593 A CN 201110216593A CN 102419457 B CN102419457 B CN 102419457B
Authority
CN
China
Prior art keywords
image
structural surface
drilling hole
television
structural plane
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
Application number
CN 201110216593
Other languages
English (en)
Other versions
CN102419457A (zh
Inventor
黄达
张永兴
金华辉
杨海清
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN 201110216593 priority Critical patent/CN102419457B/zh
Publication of CN102419457A publication Critical patent/CN102419457A/zh
Application granted granted Critical
Publication of CN102419457B publication Critical patent/CN102419457B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法,包括以下顺序步骤:一、在钻孔口上标出地理方位点,并进行钻孔电视及图像提取;二、在电视平展图像上标出相应的地理方位点,并建立二维坐标系;三、将需要确定产状的岩体结构面在图像上标出,选取结构面上不共线三点并且量测其平面坐标值;四、将图像按钻孔形状卷成圆筒,竖直放置,建立三维直角坐标系;五、运用空间向量几何理论,计算结构面法线的向量,计算该结构面的产状。本发明克服了三钻孔确定深部结构面费用昂贵及同一结构面辨识误差和其它物探方法探测相对小规模结构面精度很差等不足,充分利用了单个钻孔的结构面信息,适应于各类结构面产状的统计分析。

Description

一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法
技术领域
本发明涉及一种适用于深部岩体结构面产状测量及统计的方法。
背景技术
本发明中所提到的“深部”是指地表以下或开挖面以内一定深度,即采用常规的探勘等(如探洞和探槽等)措施不便施工或无法到达的深部岩体,但钻探技术完全可以到达的地下岩体。目前,对于深部岩体结构面产状的测量,主要是根据相邻三个钻孔揭露的同一结构面空间位置,通过空间几何理论计算得出。这一方法存在三个方面的缺点:不同钻孔揭露的同一结构面的识别易于误判;三个钻孔并不一定能同时揭露出一些规模不大的结构面;深部钻探费用昂贵,三个钻孔成本过高。现也有采用高密度电法和EH4电磁探测成像系统等高科技手段测量深部相对较大规模的矿脉或断层等产状,但这类方法仅仅适应于大规模的地质结构面,对于小规模的裂隙性结构面并不适宜。
发明内容
本发明的目的是提供一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法。
为了实现上述目的,采用以下技术方案:一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法,其特征在于包括以下顺序步骤:
一、在钻孔口上标出地理方位点,并进行钻孔电视及图像提取;
二、在电视平展图像上标出相应的地理方位点,并建立二维坐标系;
三、将需要确定产状的岩体结构面在图像上标出,选取结构面上不共线三点并且量测其平面坐标值;
四、将图像按钻孔形状卷成圆筒,竖直放置,建立三维直角坐标系;
五、运用空间向量几何理论,计算结构面法线的向量n(n1,n2,n3),计算该结构面的产状;
倾角:
Figure BDA0000079698470000021
倾向:
Figure BDA0000079698470000022
式中:m表示点N(n1,n2)所在象限数(1、2、3、4),当在第一象限时,m值为1,在第二象限时,m值为2,在第三象限时,m值为3,在第四象限时,m值为4。
本发明克服了三钻孔确定深部结构面费用昂贵及同一结构面辨识误差和其它物探方法探测相对小规模结构面精度很差等不足,本发明的方法不仅经济地充分利用了单个钻孔的结构面信息,而且适应于各类结构面(特别是小规模裂隙等)产状的统计分析。由于电视图像为钻孔表壁成像,不存在钻孔扰动影响,非常客观地反映了结构面的局部产状特征。
附图说明
图1含有结构面的钻孔电视图像平展示意图及二维坐标系;
图2含有结构面的钻孔电视图像圆柱图及三维空间直角坐标系;
图3空间向量示意图;
图4结构面法向量投影图;
图5实际钻孔电视图像a;
图6实际钻孔电视图像b。
图中标记为:件1为结构面;件2为钻孔表壁;件3为三维直角坐标系(x,y,z);件4为二维直角坐标系(d,h);A、B、C为待测结构面钻孔电视图像上的点;n为待测结构面的法向矢量。
具体实施方式
下面结合附图及具体实施例对本发明做进一步说明。
具体的步骤如下:
(1)首先对钻孔冲洗干净,并最好用鼓风机吹干孔壁。采用罗盘在孔口标出地理方位(NESW)点位置。按钻孔电视摄像设备操作规程完成成像工作。
(2)在钻孔电视平展图像上标注钻孔孔口相应标注的地理方位点位置,在平展图上建立相应的图像深度与孔周位置二维直角坐标系(d,h),坐标原点为平展图中地理方位N点。如图1所示,其中H为钻孔的深度。
(3)将需要确定产状的结构面在电视图像用笔标出,选取该结构面上任意三个点A、B、C,并且在如图1中量测和读其相应的位置坐标:A(A1,A2),B(B1,B2),C(C1,C2),量测三点在钻孔方向的深度。
(4)将已标注方位点的电视图像按钻孔形状卷成圆筒,在钻孔口部圆心建立三维空间直角坐标系,如图2所示。以圆筒筒口圆心为坐标原点、正北地理方位为坐标轴y方向、坐标轴z竖直向上。
(5)运用数学几何换算,将平展图像的二维坐标转化为相应的空间三维直角坐标。如二维坐标A(A1,A2),B(B1,B2),C(C1,C2)转化成相应的空间三维直角坐标为:
A ( sin ( A 1 / r ) · r , cos ( A 1 / r ) · r , - A 2 ) B ( sin ( B 1 / r ) · r , cos ( B 1 / r ) · r , - B 2 ) C ( sin ( C 1 / r ) · r , cos ( C 1 / r ) · r , - C 2 ) , 记为 A ( a 1 , a 2 , a 3 ) B ( b 1 , b 2 , b 3 ) C ( c 1 , c 2 , c 3 ) - - - ( 1 )
其中:r为钻孔的内径;
(6)运用空间向量几何理论,计算结构面法线的向量n(x,y,z),如图3所示。通过结构面在空间直角坐标系中三个点的坐标:A(a1,a2,a3),B(b1,b2,b3),C(c1,c2,c3),求取向量
Figure BDA0000079698470000043
Figure BDA0000079698470000044
向量积即为法向量n(x,y,z):
x = ( a 2 - b 2 ) ( a 3 - c 3 ) - ( a 3 - b 3 ) ( a 2 - c 2 ) y = ( a 3 - b 3 ) ( a 1 - c 1 ) - ( a 1 - b 1 ) ( a 3 - c 3 ) z = ( a 1 - b 1 ) ( a 2 - c 2 ) - ( a 2 - b 2 ) ( a 1 - c 1 ) - - - ( 2 )
为便于记录,将向量n(x,y,z)记录为n(n1,n2,n3)。
(7)运用数学向量关系,计算该结构面的产状。
确定倾角:倾角即为该结构面与水平面xy平面的夹角α(锐角),求得夹角(二面角):
Figure BDA0000079698470000047
确定倾向:将结构面的法向量n(n1,n2,n3)(由于结构面倾向的唯一性,若n3为负则取n(-n1,-n2,-n3))投影到水平面xy平面,求得xy平面向量N(n1,n2),向量N(n1,n2)与y轴正方向(正北方向)夹角即为倾向方位角θ(如图4所示),可由下式计算获得θ:
Figure BDA0000079698470000051
式中:m表示点N(n1,n2)所在象限数(1、2、3、4),如图4所示,
如第1象限m=1,第四象限m=4。
以上二维和三维坐标原点可以任意选取,其确定产状方法及原理一样。上面步骤中分别选取地理方位N点及圆筒顶部中心为坐标原点,是为了计算简单和规范化。
实例1:
根据现场钻孔实际情况;钻孔深度为12m,孔径为130mm钻孔电视图像进行处理(如图5所示),读取图像中结构面中不共线三个点的位置坐标(单位为m):A1(0,-4.25)、B1(0.21,-4.35)、C1(0.3266,-4.3);进行坐标转化代人公式(1)得到空间直角坐标:A1(0,0.065,-4.25)、B1(-0.0058,-0.0647,-4.35)、C1(-0.0619,0.02,-4.3);代人法向量公式(2)求得结构面法向量为n(0.001985,0.00648,-0.0083);代人公式(3)求取倾角为39.3°;截取平面投影图求得N(n1,n2)为:N(-19.85,-64.8),位于第3象限;代人公式(4)求得其倾向为197°;则确定结构面产状为197°∠39.3°。
实例2:
根据现场钻孔实际情况,提取竖直钻孔深度为12m,孔径为130mm钻孔电视图像进行处理(如图6所示),读取图像中结构面不共线三个点的位置坐标(单位为m):A2(0,-10.7)、B2(0.16,-10.55)、C2(0.3266,-10.65);进行坐标转化代人公式(1)得到空间直角坐标:A2(0,0.065,-10.7)、B2(0.0408,-0.0505,-10.55)、C2(-0.0619,0.02,-10.65);代人法向量公式(2)求得结构面法向量为;n(0.000975,-0.011325,-0.008985);代人公式(3)求取倾角为51.68°;截取平面投影图求得N(n1,n2)为:N(-9.75,113.25),位于第4象限;代人公式(4)求得其倾向为315°;则确定两结构面产状为315°∠51.68°。
岩体结构面并不一定呈现平面展布,可能存在起伏或弯折等特征。由于钻孔孔径大小的局限性,因此根据单一竖向钻孔图像确定的结构面产状为局部产状特征,但是如果根据邻近多个钻孔图像的结构面局部产状,就可以非常准确地确定结构面的空间展布规律及整体产状特征。
本发明所述的方法使用操作方便,通过简单的测量和程序计算就能准确得出深部岩体结构面的产状,而且可以对钻孔揭露的全孔岩体结构面产状进行测量和统计分析。钻孔电视图像一般是物探领域用来分析岩体的破碎程度的勘探方法,本发明是对已有资料附加应用,其不但易操作和较准确,且充分利用了钻孔电视图像,并可编写成简单的程序软件,适合推广应用。

Claims (1)

1.一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法,其特征在于包括以下顺序步骤:
一、在钻孔口上标出地理方位点,并进行钻孔电视及图像提取;
二、在电视平展图像上标出相应的地理方位点,并建立二维坐标系;
三、将需要确定产状的岩体结构面在图像上标出,选取结构面上不共线三点并且量测其平面坐标值;
四、将图像按钻孔形状卷成圆筒,竖直放置,建立三维直角坐标系,运用数学几何换算,将平展图像的二维坐标转化为相应的空间三维直角坐标;
五、运用空间向量几何理论,计算结构面法线的向量n(n1,n2,n3),计算该结构面的产状;
倾角: α = π 2 - arcsin ( ( 0,0,1 ) ( n 1 , n 2 , n 3 ) | 0,0,1 | | n 1 , n 2 , n 3 | ) arccos ( ( 0,0,1 ) ( n 1 , n 2 , n 3 ) | 0,0,1 | | n 1 , n 2 , n 3 | )
倾向: 
Figure FDA0000252876693
式中:m表示点(n1,n2)所在象限数,可取的值为1、2、3或4。
CN 201110216593 2011-07-29 2011-07-29 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法 Expired - Fee Related CN102419457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110216593 CN102419457B (zh) 2011-07-29 2011-07-29 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110216593 CN102419457B (zh) 2011-07-29 2011-07-29 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法

Publications (2)

Publication Number Publication Date
CN102419457A CN102419457A (zh) 2012-04-18
CN102419457B true CN102419457B (zh) 2013-04-03

Family

ID=45943950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110216593 Expired - Fee Related CN102419457B (zh) 2011-07-29 2011-07-29 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法

Country Status (1)

Country Link
CN (1) CN102419457B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116190B (zh) * 2013-01-26 2015-11-25 袁雄飞 井中三分量磁测3d成图方法
CN103983228A (zh) * 2014-06-04 2014-08-13 长江工程地球物理勘测武汉有限公司 利用电子岩芯测量钻孔裂隙产状的方法
CN106053256A (zh) * 2016-08-09 2016-10-26 南华大学 一种岩体结构面抗剪强度指标计算的方法
CN109236275A (zh) * 2018-09-26 2019-01-18 中国电建集团成都勘测设计研究院有限公司 基于钻孔电视摄影结构面产状的解译方法
CN109611082B (zh) * 2018-12-29 2022-05-24 武汉固德超前高新科技研发有限公司 钻孔电视的双孔测量方法
CN113236345B (zh) * 2021-06-17 2022-09-02 西安科技大学 一种钻孔裂隙可视化系统的设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4357660A (en) * 1973-05-01 1982-11-02 Schlumberger Technology Corporation Formation dip and azimuth processing technique
US5675488A (en) * 1994-05-12 1997-10-07 Halliburton Energy Services, Inc. Location determination using vector measurements
CN101882171A (zh) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 一种交互式隧道及围岩体三维模型快速构建方法
CN102042822A (zh) * 2009-10-22 2011-05-04 宝钢集团上海梅山有限公司 一种巷道断面测量方法
CN102096944A (zh) * 2010-11-25 2011-06-15 中国水电顾问集团西北勘测设计研究院 地质体结构面三维激光扫描点云识别方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4357660A (en) * 1973-05-01 1982-11-02 Schlumberger Technology Corporation Formation dip and azimuth processing technique
US5675488A (en) * 1994-05-12 1997-10-07 Halliburton Energy Services, Inc. Location determination using vector measurements
CN102042822A (zh) * 2009-10-22 2011-05-04 宝钢集团上海梅山有限公司 一种巷道断面测量方法
CN101882171A (zh) * 2010-05-21 2010-11-10 中交第二公路勘察设计研究院有限公司 一种交互式隧道及围岩体三维模型快速构建方法
CN102096944A (zh) * 2010-11-25 2011-06-15 中国水电顾问集团西北勘测设计研究院 地质体结构面三维激光扫描点云识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
三峡地下厂房岩体结构及围岩稳定性控制;胡宇庭等;《地下空间与工程学报》;20110630(第03期);526-533页 *
胡宇庭等.三峡地下厂房岩体结构及围岩稳定性控制.《地下空间与工程学报》.2011,(第03期),526-533页.
锦屏Ⅰ级水电站地下厂房施工期围岩变形开裂特征及地质力学机制研究;黄润秋等;《岩石力学与工程学报》;20110131(第01期);23-35页 *
黄润秋等.锦屏Ⅰ级水电站地下厂房施工期围岩变形开裂特征及地质力学机制研究.《岩石力学与工程学报》.2011,(第01期),23-35页.

Also Published As

Publication number Publication date
CN102419457A (zh) 2012-04-18

Similar Documents

Publication Publication Date Title
CN102419457B (zh) 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法
CN105550448B (zh) 基于钻井轨迹设计参数的预钻井三维井眼建模方法及装置
CN106643792B (zh) 惯性测量单元和地磁传感器整体标定装置及标定方法
CN103697854A (zh) 一种非接触式结构面产状测量方法
CN109883450B (zh) 埋地钢质管道内检测器磁标的定位方法
CN104060982B (zh) 测距式井下钻孔开孔方位角测量方法
CN106595567A (zh) 地质结构面产状测量方法
CN202420501U (zh) 一种gps rtk中测量隐蔽点点位坐标的辅助测量装置
CN102108856B (zh) 小角度井斜姿态测量方法及装置
CN106988726A (zh) 高精度的井眼轨迹监测方法
CN103353612B (zh) 一种地下目标物体的测量定位设备及测量定位方法
CN105652303B (zh) 一种基于Pad等移动电子设备的山区野外采样快速不接触定位方法
CN105066949A (zh) 巷道空间形变测量装置及基于该装置的测量方法
Zhao et al. A detection system for pipeline direction based on shielded geomagnetic field
CN105424013A (zh) 一种地面坐标的测量方法
CN105890579A (zh) 一种地下管线要素的快捷测量仪器与技术方法
CN104406566B (zh) 一种矿井智能综合地质测量仪及其测量方法
CN107504959B (zh) 利用倾斜航空影像测量房屋墙基轮廓的方法
CN104978476A (zh) 利用智能手机进行室内地图现场补测的方法
CN104111063B (zh) 一种基于磁场的无线三维倾角传感器及其检测方法
CN110471122A (zh) 一种基于高密度电法温纳装置的三维电阻率测深应用方法
CN105842426A (zh) 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法
CN107192326B (zh) 基于gps数据计算岩层矢量真厚度的方法
CN102621583B (zh) 井中磁场三分量梯度测量对井外磁性体定向、定位的方法
CN203375934U (zh) 基坑围护结构分层水平位移的测量装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130403

Termination date: 20150729

EXPY Termination of patent right or utility model