测线取样法引起的岩体结构面几何偏差纠正

【摘要】重要工程构筑物的建造和运行例如公路、隧道、坝基、港口地基、桥基、高楼大多涉及到岩体。理想人造材料是连续和均匀的介质,而岩体是一个复杂的、经历了长期有些是上千亿年的物理地质和工程地质作用的地质体,因此岩体形成各种不均匀和不连续的结构,岩体结构主要由完整岩块和结构面两个部分组成。结构面定义为包含各种地质成因的不连续面,它包括层面、断层、裂隙、裂纹、节理和其它相对于周围岩石具有一些共同的特征如较低的剪切强度、微小的抗拉强度和较好的流体渗透性的力学缺陷。大量的工程实例证明结构面的几何特征很大程度上控制了工程岩体的力学性能如变形模式和稳定性。结构面的几何特征通常采用窗取样法和线取样法两种取样技术获得。本文只考虑线取样法,它包括钻孔法和测线法,集中研究测线法。测线法是统计与测线相交的结构面,而不相交的结构面则不统计,该方法被广泛使用在露头面上,尽量选择一个清晰、平坦的、相对于结构面尺寸和间距更大的岩石面,出露点一般分布于海滩陡岸、峡谷、路堑、矿场和露天矿山等,在选取测量面时应保证测量面岩体和结构面能够代表该场地的整个特征。很多自然的和开挖的岩体表面由节理面和断层面组成,当这些面被选为测量面时,此时有必要布置多条不同方位的测线进行测量,以便于全面统计三维空间内结构面网络。测线法不可避免的存在产状取样偏差,即测线优先交切与它交角比较大的结构面。Terzaghi在1965年首次对此偏差进行纠正,她用产状观测频率除以交角的正弦来获取纠正后的频率。这种方法自从提出,已经被广泛运用。Terzaghi方法应用在盲区内(与测线交角小于30°)时是无效的。1.前人提出了大量避免盲区的方法,然而,盲区以外(与测线交角大于30°)的纠偏效果却很少被注意。目前为止,盲区以外的低效来源仍然是未知,同时提高盲区以外效果的方法也未见。2.网格化投影图是Terzaghi方法的内部步骤之一,它一直被前人所使用,以致于人们习惯认为网格化是Terzaghi方法的必要步骤之一。其实,目前为止,网格化的必要性还有待检验,网格化是否会给纠编结果带来误差还有待研究。3.Tang在2013年最新的研究指出,即使采用优化措施,Terzaghi方法纠偏结果仍然具有很高的误差。目前缺乏高效的纠偏方法。本文研究主要包含三个部分：1.寻找Terzaghi方法应用于盲区外时纠偏低效来源,给出提高纠偏效果的建议。2.对于Terzaghi的方法,评估网格化的必要性。如果评估结果显示网格化是不必要的,测试非网格化代替网格化的可行性,并对这两种情况下得到的纠偏结果进行效果比较。3.基于产状在三维空间内的概率分布解,提出一种比Terzaghi方法更有效的纠偏方法。这三个方面的研究方法分别是：1.我们怀疑低效来自于Terzaghi公式。因此为探索来源,我们需要Terzaghi公式的整个推导过程。鉴于Terzaghi在1965年给出的推导较为简单,我们利用解析几何学、概率论和积分学的方法作出了更为详细和合理的推导过程。我们试图在推导中找到可能带来误差的步骤。为提高效果,我们给出两个参数(网格大小和样本密度)的建议值。网格大小是Terzaghi方法中的一个设定参数,采样密度是样本数除以倾向的上下限之差、再除以倾角的上下限之差后的一个参数。(1)通过试验测得这两个参数一系列组合的效果大小。(2)将具有最高效果的一对参数选为建议值。(3)采用四川省汶川地区实测的长石砂岩结构面产状数据验证网格大小建议值和样本密度建议值的合理性。第一,用测线法实测获取倾向和倾角数据；第二,用Terzaghi方法纠偏实测数据；第三,采用纠偏后数据,建立结构面的随机三维网络模型。在该模型中,利用与野外实测测线产状一致的虚拟测线获取与之相交结构面的倾向和倾角。为区分野外实测产状,我们将从模型中获取的产状命名为“模型产状”；第四,利用Kolmogorov-Smimov双样本检验法,测试不同网格大小下实测产状与模型产状的概率分布差异。2.评价无网格化时的误差量级。(1)有必要测试无网格化的可行性。(2)基于前述的实测产状数据比较网格化与不网格化的误差。利用纠偏结果构建结构面三维网络模型,从模型中得到55个测量产状数据,用Kolmogorov-Smimov双样本检验法测试了原始产状概率分布和模拟产状概率分布之间的误差。另一方面,利用结构面随机网络模拟开展精度测试,此模拟与实测数据相比拥有更多测线产状和样本数的组合。第一,建立结构面、露头和测线模型,并且获取与测线相交的结构面产状。为区分这些相交的产状,岩体三维空间内的产状被称为“原始产状”。然后,利用Terzaghi方法,分别在网格化和非网格化的两种情况下对获取的产状数据进行频率纠偏。紧接着,用Kolmogorov-Smimov双样本检验法测试原始产状和纠偏后产状之间的概率分布差异。(3)比较网格化与非网格化两种情况下概率分布差异大小,如果非网格化差异较小,则讨论非网格化代替网格化的可行性。3.提出一种更更有效的纠偏方法。(1)用积分方法推导产状的概率密度函数(为数值解)。(2)然后在此函数基础上提出了包括2个假设和5个步骤的数值解法。两个假设是：(a)假设每一分组内的产状概率分布与直径概率分布是相互独立的；(b)假设在岩体三维空间中,每一分组内的两个产状元素(倾向和倾角)概率分布是相互独立的。同样的假设也应适用于与测线或钻孔相交的结构面产状。五个步骤包括：(a)分组。可利用Dips软件实现。需注意的是每个分组都应各自进行分析和纠偏。(b)检验倾向与倾角之间的独立性。由前面假设可知倾向与倾角之间须保持独立,因此有必要检验实测产状的倾向与倾角是否满足这个假设。独立性测试方法有许多种,Pearson开方(x2)检验法就是其中之一。(c)根据推导的数值解纠正倾向与倾角的概率分布。(d)判断纠偏后的产状概率分布类型。(e)计算修正后的产状分布参数。此处参数选择需先考虑概率分布类型,已经由步骤4确定。例如,对于正态分布,参数为均值和标准差(或方差)；对于对数正态分布,参数为均值和对数标准差(或方差)；对于均匀分布,参数为下限和上限；对于指数分布,参数为均值。确定好参数形式后,可以利用Matlab或SPSS通过最小二乘法得到参数值。对于步骤4,在一般的统计分析中,有两种用于判断概率分布类型的方式：从概率密度曲线形状判断或从累积分布曲线形状判断。为了在两种方式之间选取较合适者,在4种产状分布类型下比较两种方式有效性。比较的程序如下：(a)不妨假定产状和其他为建模所必需的参数(如大小、密度、张开度)。其中假设的产状分布包括4种类型,分别为正态分布,对数正态分布,均匀分布和指数分布；为每种分布类型设定7种样本数,分别为50,100,150,200,300,500和1000。(b)输入这些参数,采用随机模拟方法建立结构面的三维网络模型。根据4种不同的概率分布,建立4个不同的结构面网络模型。(c)获取模型中与测线相交的结构面的产状数据。针对每种模型,采用7种不同的样本数进行测量。因此,对于4中模型,共得到28个系列的产状数据。(d)由于提出的方法是基于实测倾向和倾角之间独立性假设,因此,在纠偏前就得采用Pearson开方(x2)检验法测试此独立性。(e)为得到概率密度曲线,纠偏28个产状系列；为得到累积分布曲线,也纠偏这些系列。然后对比两者的有效性。(3)利用前述实测产状数据比较提出方法与Terzaghi方法的纠偏效果。研究得到如下结论：1.更详细、合理地推导了Terzaghi的纠偏公式。推导过程中发现,在盲区以外,Terzaghi方法由于公式中的近似替代而产生理论误差。因此,为了提高纠偏效果,通过试验和实例验证给出了网格大小和样本密度的建议值,即网格大小为2°×2°以及样本密度为0.05°-2时误差最小。2.作为Terzaghi的方法一个内部步骤,网格化投影图往往会给纠偏结果引入一个相当大的误差。试验证明非网格化是可行的。最后,网格化与非网格化的纠偏结果准确性比较表明,纠偏结果在无网格化的情况下比在网格化时更精确。3.在推导产状在三维空间的概率分布函数的过程中,发现很难求取其解析解。所以推导出了数值近似解。基于数值近似解,提出了包含两个假设和五个步骤的产状纠偏的数值解方法。4.对于提出的方法,概率密度曲线和累积分布曲线的判断分布类型效果对比显示,根据概率密度曲线很难清楚地判断产状分布类型,相反地,可以很容易地从累积分布曲线中看出。这表明,在分布类型的判断上,累积分布曲线比概率密度曲线更有效。所以相比于概率密度曲线和它的函数(概率密度函数),累积分布曲线和它的函数(累积分布函数)是最佳分布形式判断方式。5.对于提出的方法,在各种样本数下的纠偏效果对比试验表明：(1)当倾向或倾角服从指数分布时,该方法是低效的。这说明该方法不适合指数分布的产状。(2)当倾向或倾角服从均匀分布时该方法是相当有效的。且增加样本数量几乎不影响其效果。表明该方法适合均匀分布的产状。(3)当产状服从正态分布且样本数量超过500,增加样品数很难提高效果。所以应该限制样本数量使其低于500。(4)当倾向或倾角服从对数正态分布且样本数量超过150时,增加样品数不能显著提高效果。所以150为对数正态分布最佳的样本数量。由于在测量和纠偏前未知分布类型,最佳样本数量的判断可能是上述四个情况之一。因此,结合上述四种结果后,为保证提出的方法高效性,最佳样本数量应约为150。6.与Terzaghi方法相比,提出的方法更有效。原因可能是Terzaghi方法是假设在每个计数圆里的所有结构面是平行的。然而事实上,并不是所有的结构面都是平行的。这种假设与事实的不符造成了Terzaghi方法低效。而提出的方法不基于这个假设,所以其纠偏效果会更高。
【作者】黄磊；
【导师】唐辉明；
【作者基本信息】中国地质大学，地质工程，2014，博士
【关键词】岩体；结构面；测线法；几何偏差；纠正；

【参考文献】
[1]戈余丽.320T液压平板车结构有限元分析及其优化研究[D].江苏科技大学,机械制造及其自动化,2013,硕士.
[2]吴镇宇.高强度结构钢低温韧性研究[D].西安建筑科技大学,钢铁冶金,2013,硕士.
[3]常玥.严歌苓笔下的男性形象研究[D].吉林大学,中国现当代文学,2013,硕士.
[4]李爱民.三维弹性接触问题的数值模拟技术及其应用研究[D].西北工业大学,机械制造及其自动化,2004,硕士.
[5]罗培骁.鲤春病毒血症病毒G蛋白单抗和多抗的制备及特性鉴定[D].华中农业大学,水产养殖,2014,硕士.
[6]王玲.国际视点[J].微纳电子技术,2011,02:138-140.
[7]杨丹.平面Bonnesen-型等径不等式与Bonnesen-型的Ros等周不等式[D].西南大学,基础数学,2013,硕士.
[8]李龙飞.已烯雌酚磁性分子印迹聚合物的制备及其在牛奶检测中的应用[D].齐鲁工业大学,食品工程,2014,硕士.
[9]王敏.机动车预登记系统的设计与实现[D].燕山大学,计算机技术,2014,硕士.
[10]夏新海.物流配送车辆调度优化研究[D].武汉理工大学,交通运输规划与管理,2004,硕士.
[11]王娟.二十世纪二十年代历史小说的现代性[D].南京师范大学,中国现当代文学,2011,硕士.
[12]王栋.葡萄糖辅助的不对称合成反应及其在有机合成中应用的研究[D].杭州师范大学,有机化学,2007,硕士.
[13]李琪.石墨烯/四氧化三钴纳米电极材料的制备及其电容性质研究[D].陕西师范大学,2014.
[14]朱芳来,韩正之.非线性系统指数型降维观测器设计[J].自动化学报,2002,06:1015-1018.
[15]李伏超.SNP对共刺激分子CD274、CD80在胃癌中表达调控研究[D].苏州大学,消化病学,2014,硕士.
[16]唐宇.无外接电容型LDO环路稳定性的研究[D].西南交通大学,电路与系统,2014,硕士.
[17]王赛.青霉素发酵过程软测量建模方法的研究[D].东北大学,控制理论与控制工程,2010,硕士.
[18]钟绍华,杨通胜.野外地震记录干扰的监视与排除[J].石油地球物理勘探,1995,S1:164-169+186.
[19]陈歆.试述海关卡口监管现代化建设必由之路[J].上海海关高等专科学校学报,2003,02:43-45+55.
[20]牛成诚.光刻机双驱动微动台温度场研究[D].哈尔滨工业大学,机械制造及其自动化,2013,硕士.
[21]高玮.工业遗产改造中的文化景观整合与表达[D].合肥工业大学,设计艺术学,2011,硕士.
[22]张满彩.病理语音的特征提取与分类研究[D].哈尔滨工业大学,计算机技术,2014,硕士.
[23]沈伟杰.气膜冷却孔型作用的实验研究[D].浙江理工大学,车辆工程,2014,硕士.
[24]周晓霞.宏观审慎监管对我国商业银行经营的影响[D].浙江工商大学,金融,2013,硕士.
[25]张颖.现代农业园景观与CIS可识别性研究[D].南京农业大学,风景园林,2013,硕士.
[26]宫姗姗.济南市都市农业发展研究[D].华中师范大学,农业推广,2014,硕士.
[27]王晓强.VIENNA整流器的研究[D].北京交通大学,2014.
[28]史晓琳.单板式配电自动化通讯前置机的研究[D].华北电力大学（北京）,电力系统及其自动化,2003,硕士.
[29]刘文盼.基于Buffer插入的互连线性能优化[D].哈尔滨工业大学,微电子学与固体电子学,2013,硕士.
[30]李静海.防空导弹聚焦破片式战斗部对空中目标的毁伤[J].飞航导弹.2006(01)
[31]何明,傅向华,马兆丰.基于不完备信息系统的RoughSet决策规则提取方法[J].计算机应用,2003,11:6-8.
[32]张建.海洋环境下CFRP及CFRP-混凝土界面性能的耐久性研究[D].大连理工大学,结构工程,2013,硕士.
[33]马越.刑讯逼供行为的法律心理学分析[D].江西农业大学,思想政治教育,2012,硕士.
[34]陈先昌.基于卷积神经网络的深度学习算法与应用研究[D].浙江工商大学,电子与通信工程（专业学位）,2014,硕士.
[35]邵海洋.上海市典型介质中邻苯二甲酸酯类环境激素的研究[D].上海大学,环境工程,2013,硕士.
[36]张晓航,郭天松.数据挖掘过程中的五大错误[J].通信企业管理,2005,10:73-74.
[37]何建.金融市场已实现波动率预测研究[D].浙江工商大学,数量经济学,2013,硕士.
[38]方桦.论数码图像艺术基本特性及其思维方法[D].湖南师范大学,美术课程教学论,2004,硕士.
[39]刘虎兴.面向服务的铁路调度系统集成及复杂业务处理研究[D].北京交通大学,2010.
[40]VivienLoua（卫恩）.FCA-based E-Government System Modeling and Prototyping System Implementation[D].大连海事大学,计算机科学与技术,2013,硕士.
[41]程晓俊,朱卫平.浅谈超链挖掘技术在搜索引擎中的应用[J].福建电脑,2005,08:11-12.
[42]张嗣瀛,杨光红.一类复杂控制系统的全息控制[J].控制与决策,1996,04:501-505.
[43]李雨蕭.《常山貞石志》跋尾研究[D].西南大学,漢語言文字學,2013,硕士.
[44]喻莹.阮曲《幽远的歌声》解析[D].中国音乐学院,音乐学,2013,硕士.
[45]赵巍.基于EtherCAT总线的运动控制器通信系统设计与实现[D].哈尔滨工业大学,机械电子工程,2013,硕士.
[46]刘岩.齐市地区农村土地承包纠纷的调查研究[D].黑龙江大学,法律,2012,硕士.
[47]郑新春.金藤清痹颗粒治疗类风湿关节炎（毒热内蕴、湿瘀阻络证）临床研究[D].湖北中医学院,中医内科学,2004,硕士.
[48]朱振宇.工程车辆自动变速智能控制系统试验研究[D].吉林大学,2004.
[49]高博.大学生抑郁情绪的动物辅助干预研究[D].大连医科大学,应用心理学,2012,硕士.
[50]周峰.聚（己二酸—四甘醇酯）增塑剂的合成与应用[D].武汉纺织大学,纺织材料与纺织品设计,2013,硕士.