节理产状的动态聚类分析以及最大距离法统计方法的研究

成绩：题目节理产状的动态聚类分析以及最大距离法统计方法的研究专业：构造地质学姓名：学号：班级：节理产状的动态聚类分析以及最大距离法统计方法的研究摘要：裂隙岩体内的节理往往是变化的，有时甚至变化很大。测得大量节理产状以后，对其分组是一项基础性的工作，但是如何分组目前还没有很好的理论和方法。采用动态聚类分析的方法，将节理产状的样本数据划分为不同的簇，通过不同簇的概率模型计算及其簇心的分布特性比较说明了簇数分类的适度性。形成了解决这类问题的系统实用的方法。又提出了一种合理、科学和简便的节理统计方法——最大距离法。其主要思路是对野外采集的节理产状数据做最优分组,使组间距离最大,而组内距离最小。该方法逻辑简单、明了,操作过程简便,而且能有效地剔除原始数据中误差较大的数据,结果精确。通过应用,表明该统计方法正确、可靠。关键词：节理产状；聚类分析；最大距离方法；地质统计。1引言岩体中的节理裂隙系统往往由几个产状不同的节理组合成的，在野外测得的节理的倾向和倾角一般而言都是变化的，有的甚至变化很大，在测得大量节理的产状后，对大量的数据进行分组或分类是一项基础性的工作。在实际的工程中，有关学者用统计学的方法对岩石裂隙岩体围岩分类、岩体裂化程度、岩体裂隙分形技术及岩体裂隙的几何参数概率模型进行了研究，本文不考虑地理空间的相关性，将样本数据看作是相互独立的，用统计学中k均值聚类分析的方法对节理分组。聚类分析是统计学上研究分类问题的一种方法，它的任务是把所有的样本数据分配到若干的簇，使得同一个簇的样本数据聚集在簇中心的周围，它们之间距离比较近，而不同簇样本数据之间的距离比较远。但如何评判节理的分组是否合理？评判的依据又是如何？目前还没有很好的理论和方法。本文的目的就是结合工程意义，用k均值聚类分析的方法，将节理产状的样本数据划分为不同的组。用系统聚类法聚类，样品一旦划到某个类以后就不变了，这要求分类的方法比较准确，另一方面系统聚类法要存贮距离矩阵，尤其当研究的样品较多时，需占用很多的存贮单元，计算方法中的迭代法德思想给我们以启发，能否先给一个粗糙的初始分类，然后用某种原则进行修改，直至分类比较合理为止。采用这种思想产生的聚类法叫做动态聚类法。图1.1动态聚类分析流程图动态聚类分析的计算步骤（1）首先随机指定k个簇中心，用欧氏距离计算每个样本数据距簇中心的距离；（2）将每个样本数据分配到距它最近的簇中心，得到k个簇；（3）分别计算各簇中所有样本数据的均值，把它们作为各簇新的簇中心；（4）重复计算步骤（2）和步骤（3）直到k个簇中心的位置都固定，簇的分配也固定，簇中心是它的均值。2SAS软件编程本文选取了，在野外选取的20组节理的产状作为实验的数据，运用sas软件，做了动态聚类分析，得出了一下的结果。2.1Sas软件程序title;goptionsftext="宋体";datad20;inputgroup$x1-x2;cards;1127922401235094306105192116149972371381912933081016121119910129971326410142322615151216897173338181761019344720997;procprintdata=b20;run;procstandardmean=0std=1data=d20out=sta20;varx1-x2;run;procfastclusdata=sta20out=out20maxc=4list;varx1-x2;idgroup;run;procsortdata=out20;bycluster;run;procprintdata=out20;vargroupclusterDISTANCE;run;proccandiscdata=out20out=can20;varx1-x2;classcluster;run;procplotdata=can20;plotcan2*can1=cluster;run;quit;2.2sas程序得出的部分图表见附图2.3分析结果通过动态聚类分析，将20组数据分别分为2类，3类，4类数据，得出分为4类是较合理的分类。其中第14组数据单独为一类的分类较不合理，可能是野外数据采集时引起的观测数据的不合理，并不符合节理的发育特征受到区域性的应力场影响的原则。所以将数据剔除。最终得到的结果是分为三类。进行了与前人所做的该地区的节理面产状等密度图分析。分为三类符合当地的实际情况。由于该地区早期受到了北-西向的应力场的作用，发育了较多北西向或者正北倾向的节理，即分类中有9个节理分到了一类当中。图2.1节理面产状等密度图3最大距离法统计方法最大距离法的基本原则是使分组后的节理组间距离最大，而组内距离最小。现以节理倾向为例，介绍最大距离法的数学模型:设有n个节理，按其倾向从0°～360°排列。以X表示节理倾向，Xi表示第i个节理的倾向。用{i，…j}表示由第i个节理至第j个节理终止的节理组，其中1≤i≤j≤n。组间距离以Se=(Sj+1-Xj)表示，即组间距离为相邻两组中前一组最小值与后一组最大值之差。组内距离用组内极差d(i，j)=Xj-Xi表示。n个节理倾向总极差S总=Xn-X1，为一常数。设n个节理被分为k组，以Xk，0和Xk，1分别表第k组节理中第1个和最后一个节理的倾向。则内距离差和为S组内=(Xn-Xk，0)+(Xk-1，1-Xk-1，0)+…+(X1，1-X1)组间距离和为S组间=(Xk，0-Xk-1，1)+(Xk-1，0-Xk-2，1)+(X2，0-X1，1)S组内+S组间=(Xk，0-Xk-1，1)+(Xk-1，0-Xk-2，+…+(X2，0-X1，1)+(Xn-Xk，0)(Xk-1，1-Xk-1，0)+…+(X1，1-X1)=S总因此当S组间最大时，S组内最小。由于组间距离为相邻节理倾向之差，如果把n1个(Xi+1-Xi)(1≤in≤-1)按从大到小排列，把个节理分成k组只需找前k-1个(Xi+1-Xi)所对的Xi作为分组点即可。用最大距离法进行节理统计，整个过程可由电子表格MicrosoftExcel，MicrsoftAccess完成，操作非常简便。也可选取不同的统计量，做成柱状图、饼状图等，使结果更加清晰。现以某地一个测量点上得到的节理产状资料[3]为例(表1)，用最大距离法对该节理根据其倾向分组。操作过程如下(1)对节理按倾向由小到大排序(2)计算排序后相邻节理倾向之差(3)对上一步计算出来的差按从小到大排序(4)剔除误差较大的数据，并找出最大距离，对节理进行分组(5)选择一定的统计量，计算并作图。图3.1节理数据图3.2第一分组点左右的节理产状柱状图图3.3第二分组点左右的节理产状柱状图图3.3第三分组点左右节理产状柱状图可得2个最大距离(42°，39°)，将该区节理根据其倾向分为3组。取?值为25°，由于已分组的3组组内极差均在25°以内，故本区节理仅有3组。计算每组节理倾向的平均方位，分别为64°，188°，300°。其对应的频度百分比为30%，45%，25%。比较符合前面所做的动态聚类的分析结果。4结论（1）本文聚类分成的四簇时，通过sas程序所做的检验，是比较好的分类方法。在结果上结合实际，剔除了一组数据。得到最终的分类。（2）簇分得太多会使样本数据在簇内的分布不连续（3）聚类的簇数据较少时，一般需要50个左右的数据，动态聚类分析得出的结果并不是十分的准确，用最大距离法得出的结果更贴合实际。但是当样本数据较多时，动态聚类分类方案更为准确。附图：参考文献：[1]章荣岫.节理等角度统计法[J].桂林冶金地质学院报,1989,9(2):219-222.[2]王青.节理裂隙统计方法研究[J].勘察科学技术,1992,(2):27-31.[3]孙宪春等.节理产状分组的k值聚类分析及其分组结果的费歇尔分布验证法[J].岩土力学，2008，29，533-536.[4]范雷等.节理岩体结构面产状的动态聚类分析[J].岩土力学，2007,28（11）：2406-2410.[5]徐云峰等.一种简便的节理统计方法[J].地质找矿丛论，2003（12）：18（24）：262-265.80196单片机IP研究与实现,TN914.42AT89S52单片机实验系统的开发与应用,TG155.1F406基于单片机的LED三维动态信息显示系统,O536TG174.444基于单片机的IGBT光伏充电控制器的研究,TV732.1TV312基于89C52单片机的印刷品色彩质量检测系统的研究,TP391.41基于单片机+CPLD体系结构的信标机设计,TU858.3TN915.62基于单片机SPCE061A的汽车空调控制系统,TM774TM621.3带有IEEE488接口的通用单片机系统方案设计与研究,TN015基于VC的单片机软件式开发平台,TG155.1F406基于VB的单片机虚拟实验软件的研究与开发,TG155.1F406采用单片机的电阻点焊智能控制器开发,TG155.1F406基于51系列单片机的PROFIBUS-DP智能从站研究,TG155.1F406八位单片机以太网接入研究与实现,TG155.1F406基于单片机与Internet的数控机床远程监控系统的研发,R319TP319基于单片机和DSP控制的医用输液泵的研究,U467.11基于单片机控制新型逆变稳压电源的设计与仿真,F426.22TP311.52基于8位单片机的摩托车发动机电控单元软硬件的开发,TB61基于430单片机的变压器监控终端的研究,TG155.1F406逆变点焊单片机控制系统研究,TG131TG113.14单片机控制数字变量柱塞泵的研究,F426.22TP311.52基于单片机控制的高通量药物筛选及检测系统开发,R730.55R734.2MCS8051以及DS80C320单片机软核的设计,TP391基于AVR单片机的应用设计实践,TN015LPC2210单片机的KGW脉冲固体激光掩膜加工控制系统研究,TG131TG113.14基于单片机控制的交流伺服系统的多梳栉经编机的研究,TN916TP31780C196单片机在铁路客车发电机控制系统中的应用研究,TP368.1TP393基于单片机的工程车辆3参数自动换档技术研究,F426.22TP311.52削方制材机摇尺机构单片机控制装置的研制,TH213.68XC196单片机集成开发环境的研制,F426.22TP311.52基于单片机与PC的光电靶测试系统研究,O536TG174.444手机和单片机控制系统的理论与应用研究,TG155.1F406基于单片机数控实验教学绘图仪研究,TN916TP317基于单片机控制的脉冲电化学齿轮修形研究,R319TP319基于AT89S52单片机的三相电度表研究,TP274.2基于MSP430单片机的嵌入式网络终端,TH812TP368.1基于MSP430单片机地下车库通风控制系统设计,TP273.2TG333.26基于PIC单片机的血压测量计设计,TP311.5TM910.6基于单片机的标记打印机的研究与开发,TP277TG156.82基于C8051F单片机和CANbus的航空三轴伺服转台控制系统的设计与研究,TP368.1V217.2HYPERL