[发明专利]一种改进列联表法的岩体均质区划分方法

权利要求书

1.一种改进列联表法的岩体均质区划分方法，其特征在于，包括：

分别记录的两个区域的节理的产状、迹长数据，所述产状包括倾向和倾角；

改进施密特图建立新坐标系：所述施密特图为半球面在水平面的投影，投影包括中心圆和以所述中心圆为圆心的多个圆环；多个所述圆环和所述中心圆以半径由小至大的方式代表多个倾角带；每个倾角带内代表节理的倾角范围；每个所述倾角带内均等划分为多个倾向块；每个所述倾向块代表节理的倾向范围；在施密特图上增加纵向垂直于所述施密特图的迹长坐标轴，且迹长坐标轴用于显示迹长；所述迹长坐标轴按照迹长划分为n个迹长范围区间；所述倾向块、所述倾角带和所述迹长范围区间的交集构成胞元；

分别提取两个区域的节理的产状和迹长数据在所述新坐标系下投影得到两个投影图：相同倾角范围、相同倾向范围、相同迹长范围区间的节理位于其所对应的胞元内；

以两个所述投影图为基础建立列联表：以所述胞元为列、所述投影图数为行建立列联表，所述胞元中的数值为所述胞元中节理的投影频数；并根据所述列联表得到卡方检验的卡方值；

确定卡方检验的显著性水平和自由度；

根据卡方检验对两个区域的节理进行均质区划分。

2.根据权利要求1所述的一种改进列联表法的岩体均质区划分方法，其特征在于，所述卡方值χ²为：

式中：i为所述投影图，j为胞元，c为胞元总数，f为胞元中节理的观测频数；e为胞元中节理的期望频数，e＝(R_iC_j)/N，其中R为所述区域中的所有胞元中节理的总观测频数，C为两个所述投影图中同一胞元中的节理的总观测频数。

3.根据权利要求2所述的一种改进列联表法的岩体均质区划分方法，其特征在于，所述卡方检验的显著性水平的计算公式为：

f(n)＝ae^-a(n+b)+d

式中，系数a＝0.8078，b＝2.4926，d＝0.0000055。

4.根据权利要求1所述的一种改进列联表法的岩体均质区划分方法，其特征在于，所述自由度为kn-1，其中k为所述倾向快的数量与两个所述投影图在同一纵列中胞元均为空的纵列数的差值。

5.根据权利要求1所述的一种改进列联表法的岩体均质区划分方法，其特征在于，采用卡方检验的p值法进行两个区域的均质区判定。

6.根据权利要求1所述的一种改进列联表法的岩体均质区划分方法，其特征在于，分为三个倾角带，内圈的倾角带的倾角或等于0度，且小于28.1度；中间环的倾角带的倾角大于或等于28.1度，且小于59度；外圈环的倾角带的倾角大于或等于59度，且小于或等于90度；所述内圈的倾角带按照顺时针等分为四个倾向块；所述中间换的倾角带按照顺时针等分为12个倾向块；所述外圈的倾角带按照顺时针等分为18个倾向块。