基于GPS数据计算岩层矢量真厚度的方法与流程

本发明涉及地质矿藏储量技术领域，尤其涉及一种基于gps数据计算岩层矢量真厚度的方法。

背景技术：

在大比例尺地质剖面测量和固体矿床储量计算中经常遇到岩层厚度计算的问题。可见，岩层厚度的计算十分重要。岩层厚度计算可分为标量法和矢量法。我国早期的岩层厚度计算法以列昂夫斯基公式为主，这些方法的计算结果只有大小，没有方向，而且计算麻烦、速度慢、易出错、效率低。

因此需要对岩层厚度计算公式进行改进，提供一种新的于gps数据计算岩层矢量真厚度的方法。

技术实现要素：

本发明为了解决现有技术问题，设计一种基于gps数据计算岩层矢量真厚度的方法，可使岩层厚度测量完全摆脱了传统用的测绳及对通视条件要求。

本发明提供一种基于gps数据计算岩层矢量真厚度的方法，包括：

步骤1，在产状或岩性变化处定位并进行地质记录，该记录包括在岩石断层或褶皱处从一岩层层面到另一岩层层面做垂线计算出各岩层矢量真厚度，岩层矢量真厚度采用(x)式计算公式计算：

其中，某次测量(两点之间)所控制的某一岩层矢量真厚度，x：两gps经度差对应的x轴上的距离，y：两gps纬度差对应的y轴上的距离，z：两gps点高程差，λ：岩层倾向方位角，β：岩层真倾角；

步骤2，将样品采集计算结果进行编号；

步骤3，将编号资料输入计算机，并制作出地质剖面图或柱状图。

优选的，在步骤1中，地质记录包括岩性描述和产状量取。

优选的，在步骤1中，所述定位采用罗盘和gps定位仪。

与相关技术相比，本发明提供的基于gps数据计算岩层矢量真厚度的方法具有以下有益效果：

一、具有十分广泛的应用前景，迎合了地质测量电子信息化大革命的趋势；避免了导线长度、坡度量取记录中出现差错和造假；可以兼容全站仪等各种光电测距仪器所获得的数据；gps数据不存在多点累计误差，其精度相对传统法更具优势；

二、gps原始数据即为gps电子信息数据，更加适用无纸化办公的趋势；

三、采用由gps数据定点填制的地质图，在远离剖面线区域内，合理挑选一些有岩层产状的岩层分界点，则可在图上截取制作虚拟地质剖面图，计算出真实可靠的岩层真厚度，简单快捷，效率高。

附图说明

图1为根据本发明一实施例的导线与坐标关系平面图；

图2为根据本发明一实施例的导线立体空间关系图；

图3为根据本发明一实施例的当ε＞0时，导线在oce上的投影示意图；

图4为根据本发明一实施例的当ε＜时，导线在oce上的投影示意图；

图5根据本发明一实施例的当ε≥0时，导线在oce面上的投影示意图；

图6为根据本发明一实施例的导线在垂直于岩层走向铅垂面上的投影线与平面夹角大于岩层倾角时剖面示意图；

图7为根据本发明一实施例的导线在垂直于岩层走向铅垂面上的投影线与平面夹角小于岩层倾角时剖面示意图；

图8为根据本发明一实施例的x-y与x’-y’坐标系及导线关系平面示意图；

图9为根据本发明一实施例的导线在垂直于岩层走向的铅垂面上的投影线与水平面夹角大于岩层倾角时的截面示意图；

图10为根据本发明一实施例的导线在垂直于岩层走向的铅垂面上的投影线与水平面夹角小于岩层倾角时的截面示意图。

具体实施方式

以下将参考附图并结合实施例来详细说明本发明。

实施例

一种基于gps数据计算岩层矢量真厚度的方法，包括：

步骤1，在产状或岩性变化处定位并进行地质记录，定位采用gps定位仪，地质记录包括岩性描述和使用罗盘量取产状。岩层矢量真厚度采用(x)式计算公式计算：

其中，某次测量(两点之间)所控制的某一岩层矢量真厚度，x：两gps经度差对应的x轴上的距离，y：两gps纬度差对应的y轴上的距离，z：两gps点高程差，λ：岩层倾向方位角，β：岩层真倾角。

步骤2，将样品采集计算结果进行编号；

步骤3，将编号资料输入计算机，并制作出地质剖面图或柱状图。步骤1中(x)式推导过程：

首先，将(a)式简化为(b)式：

现有技术岩层矢量真厚度的采用(a)式计算：

其中：

某次测量(两点之间)所控制的某一岩层矢量真厚度

l：该次测量的导线长度(两点之间的距离)

β：所测岩层的真倾角

ε：导线前进方向上的地面坡度角(后一点对于前一点的坡角)

λ：所测岩层的倾向方位角

γ：该次测量中导线前进方向方位角(后一点对于前一点的方位角)

请同时参阅图1和图2，其中图1为根据本发明一实施例的导线与坐标关系平面图，图2为根据本发明一实施例的导线立体空间关系图。o为起点，a为终点，o-a相当于传统方式中的导线号n-n+1中的n号点和n+1号点。也即，gps中的n点和n+1点。

ob为o点所在水平面和o点所在岩层层面之相交线(岩层走向线)；同样，ca为a点所在水平面和a点所在岩层层面之相交线。

过o点，做垂直于ob的平面，交ca线于c、过c点做铅垂线交o点所在水平面于e；过a点，做垂直于ac的平面，交ob线于b、过a点做铅垂线交o点所在水平面于d。

现以o为坐标原点，oe为y′轴，以岩层倾向方向为正，也即导线前进方向与岩层倾向之间的夹角为锐角时为y′值为正，钝角时为负；ob为x′轴，oe在水平面上按顺时针旋转90度方向为正，过o点且垂直于o所在的水平面的直线为z′轴(上正下负)，如此建立临时过渡立体坐标系。对于o-a导线来说，y′＝oe、x′＝ob、z′＝of＝ec；oe方向代表所测岩层的倾向、od方向代表导线前进方向、角eod代表所测岩层的倾向方位角与导线前进方向方位角之间的夹角、角aod代表导线前进方向上的地面坡角。

由此可知，o点所在岩层层面到a点(也即c点)所在岩层层面之间的距离，就等于两岩层层面在oce(bad)面上的交线之间的距离。

如果按两层面之间的上下关系约定并赋予正负内涵，该距离就是两层面之间的矢量真厚度。

将(a)式坐标数据化：

(b)式是建立在导线长度、地面坡度、导线前进方向等众多数据基础上的一个简单的过渡公式。x′轴上的信息和岩层倾向等都无法识别和读取，具有一定的局限性。但这种局限可通过其扩展(x)式得到补充和完善。

由此完成了将(a)式简化为(b)式。

(b)式证明过程如下：

水平岩层

因β＝0，(b)式变为l恒为正。当ε为零时，sinε＝0.导线沿层面前进，岩层厚度为零；当ε为正时，l·sinε必为正值，表示自下部岩层向上部岩层所测的厚度。反之亦然。

直立岩层

因β＝90，(b)式变为由图2可见，obde为一水平面，aced铅垂面为a点所在的岩层层面，oe⊥aced，点o到aced面的距离等于oe，oe等于y′，其正负由人为规定的岩层倾向确定。

倾斜岩层

1)导线前进方向与岩层走向相同：请同时参阅图2、图3、图4，其中，图3为根据本发明一实施例的当ε＞0时，导线在oce上的投影示意图；图4为根据本发明一实施例的当ε＜时，导线在oce上的投影示意图。

此时ed与ob重合，ad＝co，y′＝0，

(1)ε＝0，sinε＝0，则在走向方向上水平前进，of＝oc＝0，两层面重叠，两点之间所控制的厚度为零。

(2)ε＞0，如图3所示，sinε为正，岩层厚度

of＝oc·cosβ＝l·sinε·cosβ亦为正；当ε＝90时，l＝oc，岩层厚度of＝l·cosβ。of为自下部岩层向上部岩层之厚度。

(3)ε＜0，如图4所示，sinε为负，岩层厚度

of＝oc·cosβ＝l·sinε·cosβ亦为负；当ε＝-90时，l＝oc，岩层厚度of＝-l·cosβ。of为自上部岩层向下部岩层之厚度。

2)导线前进方向与岩层倾向之间的夹角为锐角

请参阅图5根据本发明一实施例的当ε≥0时，导线在oce面上的投影示意图。

(1)ε≥0，如图5所示，oe为y′，恒为正；ec为z′。岩层厚度因y′、z′、sinβ、cosβ均为正值，故为正。oc为自下部岩层向上部岩层所测的岩层真厚度。ε＝0时，z′＝0，y′＝og。当ε＝90时，y′＝0，z′＝oj。

(2)ε＜0，jc为y′；oj为z′为负值。岩层厚度即为：

当导线在垂直于岩层走向线的铅垂面上的投影线与水平面的夹角ocj大于岩层倾角icj时，如图6所示，图6为根据本发明一实施例的导线在垂直于岩层走向铅垂面上的投影线与平面夹角大于岩层倾角时剖面示意图。ji始终包含于mj之中，其绝对值小于mj。为负值。oc为自上部岩层向下部岩层所测的岩层真厚度；当导线在垂直于岩层走向线的铅垂面上的投影线与水平面的夹角ocj小于岩层倾角icj时，如图7所示，图7为根据本发明一实施例的导线在垂直于岩层走向铅垂面上的投影线与平面夹角小于岩层倾角时剖面示意图。mj始终包含于ji之中，其绝对值小于ji。为正值。oc为自下部岩层向上部岩层所测的岩层真厚度；当导线在垂直于走向线的铅垂面上的交线与水平面的夹角ocj等于岩层倾角时，ji＝mj＝0，岩层厚度为零。

3)导线前进方向与岩层倾向之间的夹角为钝角

(1)通过颠倒起点-终点的方式将导线前进方向与岩层倾向之间的夹角由锐角变为钝角。借助锐角的证明来证明(b)式同样适用于钝角。

起点-终点的颠倒，必然伴随着y′值和z′值(地面坡角不为零时)正负的改变。岩层厚度的的绝对值虽然不变，但矢量也随之颠倒。

①地面坡角不为零

②地面坡角为零

当地面坡角为零，则z′必然为零

同理，

(2)通过对岩层倾向做180度旋转的方式将导线前进方向与岩层倾向之间的夹角由锐角变为钝角。

此种情况为(b)式的一个盲点。但从图8、图9、图10中可以看出，其中，图8为根据本发明一实施例的x-y与x’-y’坐标系及导线关系平面示意图；图9为根据本发明一实施例的导线在垂直于岩层走向的铅垂面上的投影线与水平面夹角大于岩层倾角时的截面示意图；图10为根据本发明一实施例的导线在垂直于岩层走向的铅垂面上的投影线与水平面夹角小于岩层倾角时的截面示意图。

如果把y′转换到x-y坐标系上，问题就会得到解决：

y′＝hf+fa′＝x·sinλ+y·cosλ

因为sin(π+λ)＝-sinλ；cos(π+λ)＝-cosλ

所以岩层倾向做180度旋转后，

当起点岩层为下，终点为上部岩层时，gl为正，反之为负。

再由(b)式扩展为(x)式：

gps坐标的y轴与地理坐标的0度(360度)重合，(b)式中的y′轴指向地理坐标的岩层倾向方位角，二者之间的夹角正好等于岩层倾向方位角。

因为z坐标是通用的，z′＝z，所以根据平面坐标转换公式：

x′＝x·cosλ-y·sinλ

y′＝y·cosλ+x·sinλ

(b)式：

最后得到(x)式：

因为(b)式和(c)式包括了所有可能出现的各种岩层与两测点之间的关系。利用y·cosλ+x·sinλ＝y′，即可把(x)式变为(b)式，(b)式已经被证明；当λ＞180时，(x)式就变为(c)式。所以(x)式可行。

对比例

用(a)式、(x)式分别对同一岩层的岩层厚度进行计算，其结果显示(x)式更为精确，从另一角度也证明了(x)式完全适用于岩层真厚度的计算。计算结果对比分析表见表1、表2。

表1：(a)式与(x)式实测数据计算结果对比表

表2：(a)式与(x)式虚拟数据计算结果对比表

从以上可以看出，gps法是建立在岩层厚度的矢量法计算通式的基础，本文通过直接或间接方法均能建立该法，由此说明本文建立的借助于gps数据的新的岩层真厚度矢量计算公式具有普适性，且计算更为简单、精准。

以上提及的术语解释如下：

岩层真厚度：指自一岩层层面上的点a向另一层面做垂线交于b点，ab之间的距离即为该两层面所包含的岩层厚度。

水平岩层是指：两岩层层面在铅垂面上的交线之间的距离。

直立岩层：指两岩层层面在水平面上的交线之间的距离。

倾斜岩层：指两岩层层面在垂直于走向线的铅垂面上的交线之间的距离。

岩层矢量真厚度：指规定由下部岩层到上部岩层所测得的岩层真厚度为正厚度，由上而下则为负厚度；反之亦然。

与相关技术相比，本发明提供的基于gps数据计算岩层矢量真厚度的方法具有以下有益效果：

一、具有十分广泛的应用前景，迎合了地质测量电子信息化大革命的趋势；避免了导线长度、坡度量取记录中出现差错和造假；可以兼容全站仪等各种光电测距仪器所获得的数据；gps数据不存在多点累计误差，其精度相对传统法更具优势；

二、gps原始数据即为gps电子信息数据，更加适用无纸化办公的趋势；

三、采用由gps数据定点填制的地质图，在远离剖面线区域内，合理挑选一些有岩层产状的岩层分界点，则可在图上截取制作虚拟地质剖面图，计算出真实可靠的岩层真厚度，简单快捷，效果高。

以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利保护范围内。