[发明专利]隧道岩层产状三维测量方法

摘要

本发明公开了一种隧道岩层产状三维测量方法，该方法以Google Earth提供的海量DEM、影像等地理信息数据为数据源，利用Google Earth COM API实现从Google Earth上获取岩层分界点数据，进行岩层面拟合，计算岩层产状，并通过标记语言实现岩层面以及岩层产状要素的可视化显示。基于Google earth平台隧道岩层产状的三维测量方法相比用地质罗盘进行野外勘测，不仅省时省力，而且扩展了勘测范围，对自然条件恶劣、岩层跨度范围大的区域，亦可方便测出其岩层产状要素，是对传统测量方法的一个很好补充。

主权项

一种隧道岩层产状三维测量方法，其特征在于，它包括如下步骤：步骤1：在Google Earth软件数据库中人工选取多个待测地表岩层分界点的大地坐标，选取的待测地表岩层分界点的个数不少于3个；步骤2：采用高斯投影将上述待测地表岩层分界点大地坐标转换为平面坐标；步骤3：在步骤2的平面坐标系中任意选取三个点确定基准平面方程：A0X+B0Y+C0Z+1=0                  （1）其中，A0、B0、C0是基准平面方程系数，x，y，z是三维空间直角坐标；步骤4：根据上述基准平面拟合出一个近似平面，即根据上述基准平面方程确定拟合平面方程，任意选取多个点所拟合的平面方程为AX+BY+CZ+1=0              （2）则拟合平面方程和基准平面方程之间满足如下关系：A=A0+A＇，B=B0+B＇，C=C0+C＇，其中A'、B'、C'为误差；步骤5：确定待测地表岩层分界点(Xi,Yi,Zi)的拟合平面方程，将A=A0+A＇，B=B0+B＇，C=C0+C＇代入步骤4中确定的拟合平面方程AX+BY+CZ+1=0，得到待测地表岩层分界点(Xi,Yi,Zi)的拟合平面方程Vi=A＇Xi+B＇Yi+C＇Zi+A0Xi+B0Yi+C0Zi+1=A＇Xi+B＇Yi+C＇Zi+Mi  （3）其中，Mi是将测量点i的坐标代入基准平面方程所得的误差，为常数；步骤6：根据最小二乘法，用待测地表岩层分界点(Xi,Yi,Zi)拟合计算岩层面方程，应使： S = Σ i = 1 n ( A ′ X i + B ′ Y i + C ′ Z i + M i ) 2 - - - ( 4 ) 最小，为求得最小值，S分别对A＇,B＇,C＇求导，并令其等于0，得： d ( Σ i = 1 n V i 2 ) d A ′ = 0 , d ( Σ i = 1 n V i 2 ) d B ′ = 0 , d ( Σ i = 1 n V i 2 ) d C ′ = 0 . - - - ( 5 ) 将式（3）代入式（5）展开，得到线性方程组： A ′ Σ i = 1 n X i 2 + B ′ Σ i = 1 n X i Y + C ′ Σ i = 1 n X i Z i + Σ i = 1 n M i X i = 0 , A ′ Σ i = 1 n X i Y + B ′ Σ i = 1 n Y i 2 + C ′ Σ i = 1 n Y i Z i + Σ i = 1 n M i Y i = 0 , A ′ Σ i = 1 n X i Z i + B ′ Σ i = 1 n Y i Z i + C ′ Σ i = 1 n Z i 2 + Σ i = 1 n M i Z i = 0 . - - - ( 6 ) 解线性方程组（6），可求得A＇,B＇,C＇，将A＇,B＇,C＇代入A=A0+A＇，B=B0+B＇，C=C0+C＇，拟合出岩层面方程为：(A0+A＇)X+(B0+B＇)Y+(C0+C＇)Z+1=0   （7）其中，Z为岩层边界点的高程；步骤7：得到待测岩层在空间分布状态的产状要素描述信息，上述要素包括岩层的走向、倾向和倾角，其中，待测岩层的走向为岩层面与水平面的交线；待测岩层的倾向为岩层面上最大倾斜线在水平面上投影所指的方向，计算方法是在上述岩层走向的基础上加或减90°；待测岩层的倾角为岩层面与水平面的夹角；步骤8：获取待测隧道中线上任意一点的大地坐标和方向，步骤9：根据待测地表岩层面向待测隧道中线延伸相交的位置，确定待测隧道与待测地表岩层面的分界线，其中待测隧道中线起点 和终点的三维空间坐标分别为(X1,Y1,Z1)和(X2,Y2,Z2)，则待测隧道中线方程为： X - X 1 X 2 - X 1 = Y - Y Y 2 - Y 1 = Z - Z 1 Z 2 - Z 1 - - - ( 10 ) 求解方程组： X - X 1 X 2 - X 1 = Y - Y Y 2 - Y 1 = Z - Z 1 Z 2 - Z 1 ( A 0 + A ′ ) X + ( B 0 + B ′ ) Y + ( C 0 + C ′ ) Z + 1 = 0 - - - ( 11 ) 方程组的解即为待测地表岩层面与待测隧道的中线在三维空间中的交点，依次计算出各交点到线路中线上的投影里程，即得到待测隧道岩层分界点里程；步骤10：在空间坐标系下对式（7）的岩层面方程进行可视化操作，在该可视化操作中须为其拟定一个边界，取步骤1中待测岩层分界点平面坐标（x、y）的最大值和最小值点作为边界点，记为（MaxX、MaxY）、（MaxX、MinY）、（MinX、MaxY）和（MinX、MinY）,连接上述四个点（MaxX、MaxY）、（MaxX、MinY）、（MinX、MaxY）和（MinX、MinY）所组成的平面，即为待测岩层面边界；将上述四个边界点数据代入式（7）的岩层面方程，得到边界点的高程Z；考虑到Google Earth软件采用的是大地坐标系，因此还须进行高斯投影坐标反算，得到待测岩层边界点的大地坐标；步骤11：编辑Google Earth软件中的标记语言文件，在标记语言文件的“多边形”标签中添加步骤10中获得的待测岩层边界点的大地坐标数据，并在标记语言文件的“描述”标签中添加步骤7中得到的待测岩层的走向、倾向和倾角信息，同时，进行标记语言文件的实时生成与加载，完成计算结果的实时可视化显示，这样就实现了隧道岩层产状三维测量。