发明名称 一种基于孔壁图像的相邻钻孔结构面连通性分析方法
摘要 本发明公开了一种基于孔壁图像的相邻钻孔结构面连通性分析方法,步骤是:获取结构面参数信息;局部坐标系下结构面特征点的坐标计算;全局坐标系下结构面特征点的坐标归一化处理;计算两两特征点之间的距离,初步分析结构面的连通性;连通性再判断。本发明对初步分析具备连通性的结构面进行特征匹配,进一步确认结构面的连通性。方法简单易行,稳定可靠,操作简便,工作效率高,能够快速准确地对钻孔间的结构面进行连通性分析。
申请公布号 CN105741329B 申请公布日期 2017.04.19
申请号 CN201610056069.5 申请日期 2016.01.27
申请人 中国科学院武汉岩土力学研究所 发明人 韩增强;王川婴;邹先坚
分类号 G06T7/40(2017.01)I 主分类号 G06T7/40(2017.01)I
代理机构 武汉宇晨专利事务所 42001 代理人 李鹏;王敏锋
主权项 一种基于孔壁图像的相邻钻孔结构面连通性分析方法,其特征在于,包括以下步骤:步骤1、根据第一钻孔的钻孔孔壁的平面展开图获得第一结构面倾向值α<sub>1</sub>、第一结构面倾角值β<sub>1</sub>、第一结构面位置深度d<sub>1</sub>;根据第二钻孔的钻孔孔壁的平面展开图获得第二结构面倾向值α<sub>2</sub>、第二结构面倾角值β<sub>2</sub>、第二结构面位置深度d<sub>2</sub>;步骤2、设定第一局部坐标系,第一局部坐标系的原点为第一钻孔孔口中心点,X轴正向指向正东方向,Y轴正向指向正北方向,Z轴正向垂直向上;设定第二局部坐标系,第二局部坐标系的原点为第二钻孔孔口中心点,X轴正向指向正东方向,Y轴正向指向正北方向,Z轴正向垂直向上,通过第一结构面倾向值α<sub>1</sub>、第一结构面倾角值β<sub>1</sub>、第一结构面位置深度d<sub>1</sub>获得获得第一局部坐标系下的第一结构面特征点F<sub>1L</sub>(X<sub>1L</sub>,Y<sub>1L</sub>,Z<sub>1L</sub>);通过第二结构面倾向值α<sub>2</sub>、第二结构面倾角值β<sub>2</sub>、第二结构面位置深度d<sub>2</sub>获得获得第二局部坐标系下的第二结构面特征点F<sub>2L</sub>(X<sub>2L</sub>,Y<sub>2L</sub>,Z<sub>2L</sub>);步骤3、选取一点作为全局坐标系原点,全局坐标系原点的X轴正向指向正东方向,Y轴正向指向正北方向,Z轴正向垂直向上,根据第一结构面特征点F<sub>1L</sub>(X<sub>1L</sub>,Y<sub>1L</sub>,Z<sub>1L</sub>)获得第一结构面所在平面在全局坐标系下的方程,进而获得第一结构面在全局坐标系下对应的特征点F<sub>1</sub>(X<sub>1</sub>,Y<sub>1</sub>,Z<sub>1</sub>);根据第二结构面特征点F<sub>2L</sub>(X<sub>2L</sub>,Y<sub>2L</sub>,Z<sub>2L</sub>)获得第二结构面所在平面在全局坐标系下的方程,进而获得第二结构面在全局坐标系下对应的特征点F<sub>2</sub>(X<sub>2</sub>,Y<sub>2</sub>,Z<sub>2</sub>);所述的第一结构面在全局坐标系下对应的特征点F<sub>1</sub>(X<sub>1</sub>,Y<sub>1</sub>,Z<sub>1</sub>)基于以下公式获得:<maths num="0001"><math><![CDATA[<mrow><msub><mi>X</mi><mn>1</mn></msub><mo>=</mo><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000021.GIF" wi="982" he="133" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>Y</mi><mn>1</mn></msub><mo>=</mo><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000022.GIF" wi="941" he="134" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>1</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>1</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000023.GIF" wi="966" he="135" /></maths>其中,X<sub>O1L</sub>,Y<sub>O1L</sub>,Z<sub>O1L</sub>为第一局部坐标系的原点在全局坐标系中的坐标;所述的第二结构面在全局坐标系下对应的特征点F<sub>2</sub>(X<sub>2</sub>,Y<sub>2</sub>,Z<sub>2</sub>)基于以下公式获得:<maths num="0004"><math><![CDATA[<mrow><msub><mi>X</mi><mn>2</mn></msub><mo>=</mo><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000024.GIF" wi="1022" he="127" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>Y</mi><mn>2</mn></msub><mo>=</mo><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000025.GIF" wi="985" he="135" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>Z</mi><mn>2</mn></msub><mo>=</mo><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>X</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Y</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mo>*</mo><msub><mi>Z</mi><mrow><mi>O</mi><mn>2</mn><mi>L</mi></mrow></msub></mrow><mrow><msup><msub><mi>X</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Y</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>Z</mi><mrow><mn>2</mn><mi>L</mi></mrow></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001215589470000026.GIF" wi="1004" he="135" /></maths>其中,X<sub>O2L</sub>,Y<sub>O2L</sub>,Z<sub>O2L</sub>为第二局部坐标系的原点在全局坐标系中的坐标,所述的第一局部坐标系下的第一结构面特征点F<sub>1L</sub>(X<sub>1L</sub>,Y<sub>1L</sub>,Z<sub>1L</sub>)基于以下公式获得:X<sub>1L</sub>=d<sub>1</sub>cosβ<sub>1</sub>sinβ<sub>1</sub>cosα<sub>1</sub>Y<sub>1L</sub>=d<sub>1</sub>cosβ<sub>1</sub>sinβ<sub>1</sub>sinα<sub>1</sub>Z<sub>1L</sub>=d<sub>1</sub>cosβ<sub>1</sub>cosβ<sub>1</sub>所述的第二局部坐标系下的第二结构面特征点F<sub>2L</sub>(X<sub>2L</sub>,Y<sub>2L</sub>,Z<sub>2L</sub>)基于以下公式获得:X<sub>2L</sub>=d<sub>2</sub>cosβ<sub>2</sub>sinβ<sub>2</sub>cosα<sub>2</sub>Y<sub>2L</sub>=d<sub>2</sub>cosβ<sub>2</sub>sinβ<sub>2</sub>sinα<sub>2</sub>Z<sub>2L</sub>=d<sub>2</sub>cosβ<sub>2</sub>cosβ<sub>2</sub>步骤4、设定特征点F<sub>1</sub>(X<sub>1</sub>,Y<sub>1</sub>,Z<sub>1</sub>)与特征点F<sub>2</sub>(X<sub>2</sub>,Y<sub>2</sub>,Z<sub>2</sub>)之间的距离为|F<sub>1</sub>F<sub>2</sub>|,若|F<sub>1</sub>F<sub>2</sub>|小于判别标准ε,则第一结构面与第二结构面初步分析具备连通性;若|F<sub>1</sub>F<sub>2</sub>|大于等于判别标准ε,则第一结构面与第二结构面不具备连通性。
地址 430071 湖北省武汉市武昌区水果湖街小洪山2号