[发明专利]一种岩体结构面极点图和产状玫瑰花图智能绘制系统

权利要求书

1.一种岩体结构面极点图和产状玫瑰花图智能绘制系统，其特征在于，所述方法包括以下步骤：

1：数据智能导入模块

用于将野外测量的结构面数据，智能导入到软件系统中；

2：模糊等价聚类算法智能计算模块

基于模糊等价聚类算法，智能实现结构面产状的模糊等价聚类，过程如下：

2.1：设结构面的实测样本数为N，第i个样本表示为(x_i1，x_i2)，x_i1为结构面倾向，x_i2为结构面倾角，其模糊关系矩阵R为：

矩阵中的元素r_ij为第i个样本与第j个样本间的相似系数，表征其两者间的相似程度；r_ij越大表示样本i与样本j越相似；

2.2：计算相似系数r_ij：

式中，i＝1，2...N；j＝1，2...N；c为计算参数(0≤c≤1)，适当选取c值使得r_ij在[0，1]中分散开来；

2.3：求解闭包t(R)：

2.4：进行结构面分组判断：模糊矩阵乘法运算步骤与普通矩阵乘法相似，不同的是，并非先两项相乘后再相加，而是先取小后取大，若C＝AB，则C中的元素n级模糊关系矩阵R，即为n个R连乘；即：

当Rⁿ＝Rⁿ⁺¹＝Rⁿ⁺²＝…时 (5)

有，模糊等价矩阵t(R)＝Rⁿ (6)

取定截集水平λ∈[0，1]，若t(R)中r_ij≥λ，则结构面i和j属于同一类；即

r_ij≥λ (7)

3：极点图智能绘制模块

根据聚类结果和结构面分组结果，智能绘制出结构面产状极点图，过程如下；

3.1：采用下半球等角度投影方法绘制极点图；

3.2：将以倾向α_d和倾角β_d表示的节理产状数据转换为以节理单位法向量表示的结构面产状数据，设α_n和β_n分别为结构面单位法向量的倾伏向和倾伏角，对于任意结构面的单位法向量表示为X＝(x₁，x₂，x₃)，此时半球面上每个点都对应一个节理产状，公式为：

X＝(x₁，x₂，x₃) (8)

α_d∈(0，360)，β_d∈(0，90) (11)

3.3：得到以单位法向量表示的结构面产状数据；

3.4：基于结构面法向产状数据，根据结构面空间赤平投影图的纵剖面原理，A’点为该平面法线的赤面投影，结合赤平投影原理，计算出A’在赤平投影图上的坐标x_n和y_n，公式如下：

3.5：求解出所有结构面法线的赤平投影坐标点(x_n，y_n)；

3.6：绘制一条直径为单位长度的基圆，绘制出铅直和水平两条直径，并标出E、S、W、N；

3.7：将所有结构面的赤平投影坐标(x_n，y_n)，绘制在基圆图上；

3.8：智能实现结构面产状极点图的绘制；

4：结构面智能统计分析模块

用于对聚类后的结构面智能统计分析，获取每组结构面的倾向、倾角、迹长、间距、断距的均值和方差，过程如下：

4.1：确定样本分区区间m；

4.2：求解样本极差

4.3：计算每个分区区间M_m：

4.4：确定样本落在每个分区区间里的概率，先利用计算机循环语言统计落在每一个区间的样本个数N_m，结合样本总数N，计算样本数概率P_m：

4.5：求解样本均值

4.6：求解样本方差S²，其中S为标准差：

4.7：根据概率P_m值，自动绘制出每组结构面的倾向、倾角、迹长、间距和断距的概率分布形态；

5：数据智能输出模块

智能输出结构面产状的分类信息；

6：走向玫瑰花图智能绘制模块

根据走向玫瑰花图绘制方法，智能绘制出结构面走向玫瑰花图，过程如下；

6.1：将节理走向数据，换算成北东和北西方向，按节理走向方位角大小依次排序，每隔α＝10°角度进行分组，每组命名T_i；

T_i＝{α，α+9°} (19)

α＝10(i-1) (20)

i∈(1，10)∪(27，36) (21)

6.2：统计每组节理的数目和每组节理的平均走向

i∈(1，10)∪(27，36) (23)

6.3：根据作图的大小和各组节理数目，选取一定长度的线段代表一组节理，确定线段的比例尺L_T；

6.4：以等于按比例尺L_T表示的、数目最多的一组节理的线段长度为半径，作半圆，过圆心作南北线及东西线，在圆周上标明方位角；

6.5：对每组节理T_i，按平均走向为方位角，在半圆上做出记号，自圆心向圆周记号点半径方向，按组内数目和比例尺L_T定出一点该点即表示该组节理平均走向和节理数目；

6.6：顺次连接和如其中某组节理数为零，则连线回到圆心，再由圆心引出与下一组相连；

6.7：智能绘制出节理走向玫瑰花图；

7：倾向玫瑰花图智能绘制模块

根据倾向玫瑰花图绘制方法，智能绘制出结构面倾向玫瑰花图，过程如下；

7.1：将节理倾向数据，按节理倾向方位角大小依次排序，每隔θ＝10°角度进行分组，每组命名D_j；

D_j＝{θ，θ+9°} (26)

θ＝10(j-1) (27)

j∈(1，36) (28)

7.2：统计每组节理的数目和每组节理的平均倾向

j∈(1，36) (30)

7.3：根据作图的大小和各组节理数目，选取一定长度的线段代表一组节理，确定线段的比例尺L_D；

7.4：以等于按比例尺L_D表示的、数目最多的一组节理的线段长度为半径作圆，过圆心作南北线及东西线，在圆周上标明方位角；

7.5：对每组节理D_j，按平均倾向为方位角，在半圆上做出记号，自圆心向圆周记号点半径方向，按组内数目和比例尺L_D定出一点该点即表示该组节理平均倾向和节理数目；

7.6：顺次连接和如其中某组节理数为零，则连线回到圆心，再由圆心引出与下一组相连；

7.7：智能绘制出结构面倾向玫瑰花图；

8：倾角玫瑰花图智能绘制模块

根据倾角玫瑰花图绘制方法，智能绘制出结构面倾角玫瑰花图，过程如下；

8.1：将节理数据，先按节理倾向方位角大小依次排序，每隔θ＝10°角度进行分组，每组命名Q_j；

Q_j＝{θ，θ+9°} (33)

θ＝10(j-1) (34)

j∈(1，36) (35)

8.2：统计每组节理的数目每组节理的平均倾向和平均倾角

j∈(1，36) (37)

8.3：根据作图的大小和各组节理数目，选取一定的倾角角度代表一组节理，确定线段的比例尺L_Q；

8.4：以等于按比例尺L_Q表示的、数目最多的一组节理的线段倾角为半径作圆，过圆心作南北线及东西线，在圆周上标明方位角；

8.5：对每组节理D_j，按平均倾向为方位角，在半圆上做出记号，自圆心向圆周记号点半径方向，按组内数目和比例尺L_Q定出一点该点即表示该组节理平均倾角和节理数目；

8.6：顺次连接种如其中某组节理数为零，则连线回到圆心，再由圆心引出与下一组相连；

8.7：智能绘制出节理倾角玫瑰花图。