主权项 |
一种基于孔壁图像的相邻钻孔结构面连通性分析方法,其特征在于,包括以下步骤:步骤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>|大于等于判别标准ε,则第一结构面与第二结构面不具备连通性。 |