一种结构面产状3D统计分布估计精度的经验预测方法与流程

本发明涉及工程地质学和岩石力学领域，具体涉及一种结构面产状3d统计分布估计精度的经验预测方法。

背景技术：

岩体中发育了相互连接的结构面(或称裂隙、节理)网络，揭示裂隙化岩体介质运动(如块体运动)、力学(如强度)、水力(如流体运移)行为规律往往依赖于分离结构面网络建模，而分离结构面网络建模的实现需要输入在3d空间内的准确产状统计分布。

产状3d统计分布不能在场地测量中直接获得，实际测量中往往获得的是在2d或1d尺度下的产状统计分布，如与露头面相交的结构面的产状统计分布(2d)，和与测线相交的结构面的产状统计分布(1d)。2d和1d的实测产状统计分布与产状3d统计分布存在差别，并不能反映产状3d统计分布。

fouché提出了一种由实测产状1d统计分布来估算产状3d统计分布的方法。由于逼真的分离结构面网络建模依赖于准确的产状3d统计分布估计，因此，是否能获得逼真的网络取决于产状3d统计分布的估算精度。目前，并不存在能够预测fouché方法精度的技术。因此，有必要开发一种预测fouché方法精度的技术。

技术实现要素：

有鉴于此，本发明提供了一种结构面产状3d统计分布估计精度的经验预测方法，该经验预测方法步骤简单实用、适用性强。

本发明提供一种结构面产状3d统计分布估计精度的经验预测方法，包括以下步骤：

s1，在岩体的新鲜露头面上布置一条或多条方向一致的测线，测量测线的倾伏向和倾伏角；

s2，捕获与所述测线相交的实测结构面，统计实测结构面的样本数量，测量实测结构面的产状，得到实测结构面的产状1d样本；

s3，利用所述实测结构面的样本数量和实测结构面的产状1d样本计算fisher常量；

s4，利用所述实测结构面的产状1d样本计算优势结构面的产状；

s5，利用所述测线的倾伏向、倾伏角及优势结构面的产状计算测线与优势结构面之间的夹角；

s6，利用所述fisher常量、测线与优势结构面之间的夹角及实测结构面的样本数量计算双尾渐近显著性概率，所述双尾渐近显著性概率的数值大小指示结构面产状3d统计分布的估计精度。

进一步地，步骤s3中，fisher常量的计算公式为：

式中，κ为fisher常量，n为实测结构面的样本数量，αi为第i个实测结构面的倾向，βi为第i个实测结构面的倾角。

进一步地，步骤s4中，优势结构面的产状包括平均倾向和平均倾角，平均倾向的计算公式为：

平均倾角的计算公式为：

式中，为平均倾向，为平均倾角，n为实测结构面的样本数量，αi为第i个实测结构面的倾向，βi为第i个实测结构面的倾角。

进一步地，步骤s5中，测线与优势结构面之间的夹角的计算公式为：

式中，θ为测线与优势结构面之间的夹角，为平均倾向，为平均倾角，ψ为测线的倾伏向，ζ为测线的倾伏角。

进一步地，步骤s6中，双尾渐近显著性概率的计算公式为：

式中，p为双尾渐近显著性概率，κ为fisher常量，θ为测线与优势结构面之间的夹角，n为实测结构面的样本数量。

进一步地，根据双尾渐近显著性概率的数值大小，按照自定义准则对估计精度进行分级，自定义准则为：p≤0.05说明估计精度为极低精度；0.05<p≤0.35说明估计精度为低精度；0.35<p≤0.65说明估计精度为中精度；0.65<p≤0.95说明估计精度为高精度；p>0.95说明估计精度为极高精度。

本发明提供的技术方案带来的有益效果是：

1.本发明提供的结构面产状3d统计分布估计精度的经验预测方法填补了相关技术的空白，使得fouché方法估计精度的预测成为可能，为相关研究提供路途；

2.对于大夹角或大样本(即样本数量较大)情况，本发明提供的经验预测方法可准确预测用于估算结构面产状3d统计分布的fouché方法的精度；

3.本发明提供的经验预测方法不需要严格的假设条件，所需参数较少，步骤简单实用，适用性强。

附图说明

图1是本发明一种结构面产状3d统计分布估计精度的经验预测方法的流程示意图。

图2是本发明一实施例中测线和结构面的示意图。

图3是本发明一实施例中测线的倾伏向和倾伏角的示意图。

图4是本发明一实施例中结构面的倾向和倾角的示意图。

图5是本发明一实施例中测线与优势结构面之间的夹角的示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。

请参考图1，本发明的实施例提供了一种结构面产状3d统计分布估计精度的经验预测方法，包括以下步骤：

步骤s1，在岩体的新鲜露头面上布置一条或多条方向一致的测线，然后采用机械式或电子式地质罗盘测量测线的倾伏向ψ和倾伏角ζ，倾伏向ψ的范围是0至360°，倾伏角ζ的范围是0至90°；

参考图2，一实施例中，楔形体表示岩体1，岩体1的上端面为植被覆盖面2，岩体1的右侧面为风化面3，岩体1的正面为新鲜露头面4，斜线表示测线5，新鲜露头面4内与测线5相交的线段为结构面迹线6；

岩体1存在植被覆盖面2、风化面3及新鲜露头面4，只选择新鲜露头面4，在新鲜露头面4上布置一条或多条方向一致的测线5，采用机械式或电子式地质罗盘测量测线5的倾伏向ψ和倾伏角ζ，见图3，第一向量51是测线5的指向大地的方向向量；第二向量52是第一向量51的水平分量；第三向量8为指北线，第一夹角53是测线5的倾伏角ζ，即第一向量51和第二向量52之间的夹角；第二夹角54是测线5的倾伏向ψ，即第三向量8按顺时针方向旋转至第二向量52所转过的角度。

步骤s2，捕获与测线相交的实测结构面，统计实测结构面的样本数量n，采用机械式或电子式地质罗盘测量实测结构面的产状，得到实测结构面的产状1d样本，实测结构面的产状1d样本表示为αi/βi，i＝1,2,3,l,n，αi为第i个实测结构面的倾向，βi为第i个实测结构面的倾角，倾向的范围是0至360°，倾角的范围是0至90°。

参考图4，一实施例中，圆形表示迹线6所在的实测结构面7，第四向量71是实测结构面7的指向天空的法向量；第五向量72是第四向量71的竖直分量；第六向量73是第四向量71的水平分量；第三夹角74是实测结构面7的倾角，即第四向量71和第五向量72之间的夹角；第四夹角75是实测结构面7的倾向，即第三向量8按顺时针方向旋转至第六向量73所转过的角度。

步骤s3，利用实测结构面的样本数量n和实测结构面的产状1d样本计算fisher常量κ，fisher常量κ的计算公式为：

fisher常量κ表征产状统计分布的离散度和实测结构面的平行度，fisher常量κ值的范围从0至无穷大，fisher常量κ值越大，产状分布越集中，实测结构面相互间越平行。

步骤s4，利用实测结构面的产状1d样本计算优势结构面的产状，优势结构面的产状包括平均倾向和平均倾角

平均倾向的计算公式为：

平均倾角的计算公式为：

步骤s5，利用测线的倾伏向、倾伏角及优势结构面的产状计算测线与优势结构面之间的夹角θ，夹角θ的计算公式为：

参考图5，一实施例中，圆形表示优势结构面9，四边形表示优势结构面9所在的平面91，直线92为测线5在平面91上的投影，第五夹角93为测线5与优势结构面9之间的夹角θ，即测线5与直线92之间的夹角。

步骤s6，利用fisher常量κ、测线与优势结构面之间的夹角θ、实测结构面的样本数量n计算双尾渐近显著性概率p，双尾渐近显著性概率p的计算公式为：

式中，双尾渐近显著性概率p的范围是0至1，双尾渐近显著性概率的数值大小指示结构面产状3d统计分布的估计精度，双尾渐近显著性概率p越大指示结构面产状3d统计分布的估计精度越高。

本发明按照双尾渐近显著性概率p的数值大小，根据自定义准则，对估计精度进行分级，自定义准则为：p≤0.05说明估计精度为极低精度；0.05<p≤0.35说明估计精度为低精度；0.35<p≤0.65说明估计精度为中精度；0.65<p≤0.95说明估计精度为高精度；p>0.95说明估计精度为极高精度。

上述自定义准则可以根据实际情况进行设定，并不限制于上述分级标准。

本发明提供的结构面产状3d统计分布估计精度的经验预测方法填补了相关技术的空白，使得fouché方法估计精度的预测成为可能，为相关研究提供路途；对于大夹角或大样本(即样本数量较大)情况，本发明提供的经验预测方法可准确预测用于估算结构面产状3d统计分布的fouché方法的精度；本发明提供的经验预测方法不需要严格的假设条件，所需参数较少，步骤简单实用，适用性强。

在不冲突的情况下，本文中上述实施例及实施例中的特征可以相互结合。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。