一种获取和计算工程岩体代表性RQD值的方法与流程

本发明涉及工程岩体质量评价与分级(分类)领域，具体涉及一种获取和计算工程岩体代表性rqd值的方法。

背景技术：

岩体工程涉及范围甚广，主要包括水利工程、采矿工程、土木工程、石油工程、交通工程、军事工程等。受长期地质作用，岩体本身具有复杂的几何特性和物理力学特性，加上赋存环境复杂，造成岩体工程评价和设计多依赖经验，难以把握。为此，依据规范的测试和试验，建立一种基础性岩体质量评价方法，实现工程岩体分级(分类)，并根据行业要求与经验提出支护措施，对提高岩体工程设计与施工的合理性具有重要的意义。现有的工程岩体质量评价与分级(分类)方法主要有rmr、mrmr、q、srm、bq、gsi、srmr、m-rmr、m-rmc、cci等，这些分级方法除gsi和bq方法外均涉及一个共同的输入参数rqd(rockqualitydesignation，岩石质量指标)。由此可见rqd对岩体质量评价与分级非常重要。

然而，很多工程师和研究人员意识到rqd有一大弊端——方向性，即不同方向上的rqd有不同的值。如图1所示，沿y轴方向上的rqd最小，值为85；沿x轴或z轴方向上的rqd最大，为100。当前若利用钻孔获取rqd时一般都取竖直方向，若利用测线法获得rqd值时一般都是取竖直方向或水平方向。这种做法致使针对同一工程岩体不同的人采用不同的方法得到的rqd值不一样，有较大的盲目性，获得rqd未必能反映出岩体的节理度，显然不合理。针对rqd的方向性这一问题，需要探寻出一种获取和计算工程岩体代表性rqd值的方法。

技术实现要素：

本发明要解决的技术问题是，克服现有技术中获取rqd值方法的不足，提供一种获取和计算工程岩体代表性rqd值的方法。

为解决技术问题，本发明的解决方案是：

提供一种获取和计算工程岩体代表性rqd值的方法，包括以下步骤：

(1)结构面产状测量

测量岩石结构面的产状，测量数据的数量应满足结构面分组的要求；

(2)结构面分组

对结构面产状测量数据进行结构面分组，其组数记为n；

(3)计算各组结构面的平均产状

对于每一个结构组，分别采用式(7)计算其相应的平均产状：

式中，i为结构面组号，δmi为第i组结构面的平均倾角，θmi为第i组结构面的平均倾向，(xri,yri,zri)是i组所有结构面的向上法向量的合向量，具体由式(8)确定：

其中，δij和θij是第i组的第j个结构面的倾角和倾向，ni是第i组的结构面数量；

(4)现场rqd的测量

在岩石的露头或开挖面上布设m条不同方向的测线，m≥n；测量相邻两个结构面与测线交点的长度lj，由式(9)计算每条测线的rqd值：

其中，l′j为长度大于0.1m的lj，l为测线的总长度，j为长度大于0.1m的lj的总数量；

同时，用罗盘测量每条测线的倾伏角与倾伏向；

(5)计算各条测线方向上的结构面线密度

将步骤(4)中计算得到的各条测线方向上的rqd值代入式(1)中，即可求出其相应的结构面线密度λ：

rqd＝100e^-0.1λ(0.1λ+1)(1)

(6)计算各组结构面平均法向量上的结构面线密度

将步骤(5)中得到的各条测线方向上的结构面线密度λ代入式(6)中，得到由m个式子构成的方程组：

其中，n为结构面组的个数；

若m＝n，则直接解出n组结构面组的平均法向量上的线密度；若m>n，则利用最小二次误差法解出n组结构面组的平均法向量上的线密度；

(7)求代表性rqd及方向

将得到的结构面线密度再代入式(1)中，则得到每条测线上的rqd值；找出其中最小的rqd值即为代表性rqd值，其对应(α，β)即为代表性rqd方向。

本发明中，所述步骤(4)中，在现场采集不同测线方向上的rqd值，且测线数量应大于或等结构面组数量。

本发明中，所述步骤(6)中，利用结构面的产状和测线方向上的线密度求得平均法向量上的结构面线密度。

本发明中，所述步骤(7)中，将最小rqd选为代表性rqd值，其方向选为代表性方向。

本发明中，所述步骤(1)中，测量结构面的产状是利用地质罗盘的接触手段，或摄影测量、三维激光扫描的非接触手段实现的。所述步骤(2)中，对结构面分组的方法是节理玫瑰花图分析法、极点等密度图分析法、聚类分析法或多参数聚类分析法。所述步骤(7)中，若将不同产状方向上的rqd绘制到赤平投影图中，再作等值线则得到典型的分析结果。

发明原理描述：

本发明提出的获取和计算工程岩体代表性rqd值的方法，其主要包括：

(1)如何选定代表性rqd的方向；

(2)根据已知几个方向上的rqd值，如何计算任意方向上的rqd值；

(3)如何确定代表性rqd的值。

代表性rqd的方向的选择原理：如图1中的岩体a在沿y轴方向上的rqd值最小(等于85)，沿x轴或z轴方向上的rqd值最大(等于100)；如图2中的岩体b在沿y轴方向上的rqd值最小(等于60)，沿x轴或z轴方向上的rqd值最大(等于100)。假如将x或z方向选为代表性rqd方向，则岩体a和b的rqd值相等，这样不能反映这两个岩体的节理度的差别；假如将y方向选为代表性rqd方向，则岩体a和b之间有最大的rqd差值(等于25)，这样能够最大程度地反映出这两个岩体的节理度的差别。因此，将rqd值达到最小值的方向选为rqd的代表性方向，在这个方向的rqd值(即最小值)即为代表性rqd值。

任意方向上的rqd值计算方法：rqd与结构面线密度λ之间的关系式如下：

rqd＝100e^-0.1λ(0.1λ+1)(1)

假设结构面为有限尺寸的圆盘且均匀分布在空间，若仅考虑一组(记为i)结构面，则结构面体密度λvi与沿测线l方向上的线密度λli有如下关系式(根据图3)：

其中，e(·)为函数的期望，l是沿测线l的向量，di是结构面的直径，ni是结构面法向量。若有两条测线l1和l2，由式(2)有：

进而有：

若将结构面组的平均法向量方向nmi赋予l2，某测线方向l赋予l1，则式(4)变为：

其中λmi是沿结构面组平均法向量方向上的结构面线密度。利用叠加原理，由式(5)可得考虑多组结构面时，沿某测线方向l上的结构面线密度与结构面组平均法向量上的结构面线密度之间的关系时如下：

其中n为结构面组的个数。若测得m(m≥n)个不同方向上的rqd值，由式(1)可求得m个不同方向上的结构面线密度，将m个线密度代入式(6)，则可得到m个方程。若m等于n，则可以直接解出n组结构面组的平均法向量上的线密度；若m>n，则可以利用最小二次误差法解出n组结构面组的平均法向量上的线密度。从而根据式(6)可求得任意方向上的线密度，进而利用式(1)可求得任意方向上的rqd。

确定代表性rqd的值：记测线的倾伏角为α，倾伏向为β，α的取值范围为[0，90°]，β的取值范围为[0，360°)。将[0，90°]平均分为r份，[0，360°)平均分为q份，从而可得到(r+1)q组(α，β)。对于每一组(α，β)分别求出相应的rqd值，找出其中最小的rqd值即为代表性rqd值，其对应(α，β)即为代表性rqd方向。

与现有技术相比，本发明的有益效果是：

1、本发明明确了最小rqd为代表性rqd，给出了根据已知几个方向上的rqd值计算任意方向上的rqd值的方法，从而消除了当前获取rqd方法的弊端；

2、相比现有rqd的获取方法，本发明只需额外测一些结构面的产状，现场方便测得。

附图说明

图1是介绍rqd方向性弊端的一个实例示意图和介绍rqd代表性方向选择原理的说明实例一示意图(单位：m)；

图2是介绍rqd代表性方向选择原理的说明实例二示意图(单位：m)；

图3是结构面体密度λvi与沿测线l方向上的线密度λli关系推导原理图；

图4是结构面涉及到的几个几何参数示意图(位置、倾向、倾角、尺寸)；

图5是本发明一种获取和计算工程岩体代表性rqd值的方法的执行流程图；

图6是露头或开挖面上测线布设示意图；

图7是是一个算例的rqd分析结果等值线图。

具体实施方式

下面结合附图，对本发明作进一步的详细说明。以下的具体实施步骤可以使本专业领域的技术人员更全面的了解本发明，但不以任何方式限制本发明。

如图5所示，本发明一种获取和计算工程岩体代表性rqd值的方法包括以下步骤：

1.结构面产状测量

利用接触(如地质罗盘)或非接触(如摄影测量、三维激光扫描等)手段测量岩石结构面的产状，对于数量没有特别的要求，满足一般的结构面分组要求即可。

2.结构面分组

本发明对结构面分组没有特别的要求，可以采用传统分组方法(如依据节理玫瑰花图或极点等密度图等进行分析)，也可采用近年来新发展的分组方法(如聚类分析法、多参数聚类分析法等)，最终获得的结构面组数记为n。

3.各组结构面平均产状的计算

对于每一个结构组，可分别采用以下公式计算其相应的平均产状：

其中i为结构面组号，δmi为第i组结构面的平均倾角，θmi为第i组结构面的平均倾向，(xri,yri,zri)是i组所有结构面的向上法向量的合向量，可由下式确定：

其中δij和θij是第i组的第j个结构面的倾角和倾向，ni是第i组的结构面数量。

4.现场rqd的测量

如图6所示，在露头或开挖面上布设m条不同方向的测线(m≥n)，通过测量相邻两个结构面与测线交点的长度lj(如图6中pj和pj+1两个交点长度)，由下式计算每条测线的rqd值：

其中，l′j为长度大于0.1m的lj，l为测线的总长度，j为长度大于0.1m的lj的总数量。同时用罗盘测量每条测线的倾伏角与倾伏向。

5.计算各条测线方向上的结构面线密度

将步骤4中计算得到的各条测线方向上的rqd值代入式(1)中，即可求出其相应的结构面线密度。

6.计算各组结构面平均法向量上的结构面线密度

将步骤5中得到的各条测线方向上的结构面线密度代入式(6)中可得到由m个式子构成的方程组，若m＝n，则可以直接解出n组结构面组的平均法向量上的线密度；若m>n，则可以利用最小二次误差法解出n组结构面组的平均法向量上的线密度。例如，若m＝6和n＝4，则4组结构面平均法向量上的结果面线密度λm1、λm2、λm3、λm4可由以下式子确定：

其中：

7.求代表性rqd及方向

将测线的倾伏角α的取值范围[0，90°]平均分为r份，将测线的倾伏向β的取值范围[0，360°)平均分为q份，从而可得到(r+1)q组(α，β)。例如，若r取18，q取72，则共得到(18+1)×72＝1368条测线，这些测线的产状分别为(0°,0°)、(0°,5°)、….、(0°,355°)、(5°,0°)、….、(85°,355°)、(90°,355°)。将这些值分别代入式(6)中可得到每条测线上的结构面线密度，再将得到的结构面线密度代入式(1)中则可得到每条测线上的rqd值，找出其中最小的rqd值即为代表性rqd值，其对应(α，β)即为代表性rqd方向。若将不同产状方向上的rqd绘制到赤平投影图中，再作等值线可得到典型的分析结果见图7。

注意：本发明的实际范围不仅包括上述所公开的具体实施例，还包括在权利要求书之下实施或者执行本发明的所有等效方案。