水利工程论文-侧裂结构面迹长估计方法研究.doc

水利工程论文-侧裂结构面迹长估计方法研究摘要：侧裂面是控制拱坝抗滑稳定的重要力学边界之一，可能构成侧裂面的陡倾结构面称为侧裂结构面。如果侧裂结构面为随机基体裂缝，那么其迹长估计就是侧裂面边界条件分析中的一项重要内容。作者在分析H-H、Laslett和2L0.5三种迹长估计方法的基础上，综合确定模型方法和概率模型方法，提出了一种适用于侧裂结构面迹长估计的新方法广义H-H迹长估计算法。并利用此方法对西南某大型水电站工程侧裂结构面迹长进行了研究，证明了该方法的正确性。关键词：侧裂面侧裂结构面迹长估计广义H-H算法在水电工程中，侧裂面是控制拱坝坝肩抗滑稳定的一个重要力学边界。如果坝肩部位没有明显的陡倾断层通过，则可能构成侧裂面的结构面只能是某一走向的陡倾基体裂隙，作者将这组陡倾裂隙称之为侧裂结构面。侧裂结构面迹长是一个重要岩体结构参数，它直接反映了结构面的几何大小，所以侧裂结构面迹长是侧裂面边界条件分析中一项重要内容。但是，由于受野外露头的影响，可测量到的全迹长结构面样本数量往往很少，利用确定模型方法直接求柔构面平均迹长受到了很大的限制，因此基于概率模型的迹长估计方法便成为获柔构面平均迹长的一个主要途径。其基本原理是根据取样窗口中结构面可见迹长及其与窗口(或测线)的交切关系来推求平均迹长。尽管现有迹长估计方法较多，但是专门用于侧裂结构面迹长估计方法尚没见报道。在对西南某大型水电站工程侧裂面边界条件研究中发现，在某些情况下，各种迹长估计算法之间存在较大差别，迹长估计结果与实际情况出入很大。作者对这一现象进行了较深入的分析，并在此基础上，提出了广义H-H结构面迹长估计公式，理论和实践证明：该算法是正确可靠的。1结构面迹长估计常用方法目前已有的迹长估计方法很多，常用的方法有：Kulatilake迹长估计法、Cruden迹长估计法、Laslett(1982)迹长估计法、H-H迹长估计法和倍半迹长方法(2L0.5)等等。其中Kulatilake算法要求由倾向和倾角两个独立变量定义的结构面二维概率密度函数是已知的，所以实用性较差。Cruden公式假设迹长服从负指数分布，然而众多学者则认为迹长服从对数正态分布。鉴于上述理由，在实际计算中，主*捎昧薍-H、Laslett和2L0.5三种迹长估算方法。Laslett基于极大似然原理，推导了由窗口数据估计结构面平均迹长的极大似然公式：(1)式中：Xi是第i条两端可见结构面可见迹长；Yj是第j条一端可见结构面可见迹长；Zk是第k条两端均不可见结构面可见迹长；n为两端可见结构面条数；m为一端可见结构面条数；p为两端均不可见结构面条数。特别地，当两端可见结构面条数n和两端均不可见结构面条数p均为零时，上式变为简化为如下形式：(2)上式的物理含义是结构面平均迹长等于可见迹长数学期望的两倍，故称之为倍半迹长估计公式(简记为2L0.5).Laslett极大似然迹长估计没有考虑到测量窗口尺度对迹长估计的影响，1997年黄润秋、黄国明提出的平均迹长计算公式则弥补了这一缺陷。在宽度为w、高度为h的取样窗内，结构面与窗口的几何交切形式共有3种：(a)两端可见；(b)一端可见；(c)两端均不可见(图1).根据结构面交切窗口概率与测量窗口尺度之间的关系，黄润秋、黄国明推导了裂隙中点服从均匀分布条件下的迹长估算公式(简称H-H迹长估算公式)：图1结构面与窗口交切关系l=n1+2n0/2Nwh/w+h(3)式中：n0是两端均不可见的结构面条数；n1是一端可见另一端不可见的结构面条数；n2是两端可见的结构面条数；N是结构面总数，由下式确定：N=n0+n1+n2(4)2侧裂结构面迹长估计2.1问题特点拟建西南某大型水电站是一巨型规模的水力发电枢纽工程，大坝为双曲高拱坝，坝区侧裂结构面主要由随机分布的陡倾基体裂隙构成。按照工程设计和抗滑稳定性要求，确定左岸侧裂结构面走向变化范围为N5WN35W，右岸为N75EN75W，倾角均大于60。侧裂结构面的野外现场调查主要在勘探平硐中进行，但是由于坝区绝大多数平硐沿垂直河谷或平行河谷方向施工，且硐高仅为2m，所以平硐布置方向和硐高皆不利于侧裂面迹长调查(图2).在坝肩抗力体局部地段，尤其是在距离坝肩抗力体上游约200m处的厂房区，布置了少数与侧裂面走向近平行的勘探平硐，在其硐顶可以获得一部分全迹长资料(图3).图2坝肩抗力体低高程平硐布置及其侧裂结构面分布图3厂房区勘探平硐布置示意坝区工程勘探布置状况以及侧裂结构面发育特点决定了迹长研究应该采用以概率模型方法为主，确定性统计方法为辅的技术路线，具体方法如下：对于那些可以测量到全迹长侧裂结构面的平硐，其平均迹长可以直接利用统计方法计算获得，并将此作为迹长估计的真实值。在这些平硐附近局部范围内，往往还沿着其它方向布置一些勘探平硐，在其某一侧壁上，可获取各种与侧壁顶、底相互交切的侧裂结构面样本，侧裂结构面迹长可以利用各种概率模型方法计算获得。通过估计值与真实值之间的相互对比，判断迹长估计方法的正确性和可靠性。如果迹长估计方法是正确的，则可将这一方法推广至整个坝区侧裂结构面的迹长估计。2.2侧裂结构面迹长估计由于侧裂结构面现场调查主要是在勘探平硐中进行的，所以以平硐为单位分别计算侧裂结构面迹长，具体算法采用式(1)式(3).因篇幅所限，表1列出了坝区部分平硐侧裂结构面迹长计算结果。由表1可知，各种方法迹长估计值有如下几个特点：(1)当迹长短小时，在多数情况下，H-H方法与Laslett方法计算结果比较接近，两者之差一般不超过0.5m，2L0.5比两者均略大，它基本上反映了其平均迹长的上限。例如PD17、PD18-1、PD32、PD52以及PD45-2等。(2)当迹长较大时，H-H方法与Laslett方法迹长估计值均显著偏大。由于Laslett算法不受窗口尺度限制，所以其增加幅度明显高于H-H方法，而2L0.5方法的迹长估计值一般不超过2倍的硐高(4m).例如厂房PD13、PD45-4和PD45-5，H-H与Laslett迹长估计结果显然与实际情况不相符合。(3)H-H方法具有一定的“不稳定性”，主要表现在两个方面：一方面，当n2接近N时，其迹长估计值明显偏小，例如PD69平硐，H-H公式计算迹长值尚不足2m；另一方面，当n0较大时，其迹长估计值又明显增大，例如PD2、PD13、PD45-4和PD45-5等。表1各种算法估计迹长对比硐号Nn2n1n0H-H法迹长/mLaslett法迹长/m2L0.5法迹长/m广义H-H法迹长/mPD144152902.061.772.631.70PD2111913.102.332.172.18PD132011274.014.382.462.75PD16532021122.662.352.512.10PD1754262531.791.702.761.66PD18-1157711.861.782.311.81PD3229141231.921.522.171.69PD45-224111301.661.372.271.52PD522151512.492.092.501.95PD691716100.180.903.530.89PD69-129101182.882.602.372.34PD45-4*68025435.0910.082.703.42PD45-5*2709185.0911.182.273.49*厂房区勘探平硐。3广义H-H迹长估计法及其应用3.1广义H-H迹长估计方法为了克服上述各种迹长估计方法的缺点，作者研究了一种迹长估计新方法。此方法能够正确解决Laslett对长迹长估计的局限性，有效消除H-H方法的不稳定性，理论上较其它迹长估计算法更完善、更成熟、更实用。其基本思想是：将结构面分成两组，其中一组由交切窗口的结构面组成，平均迹长lh按照推广的H-H迹长估计方法计算获得；另一组为两端可见结构面，其平均迹长l2等于全迹长均值，具体算法如下：l=n2/Nl2+(1-n2/N)lh(5)式中：l2和lh由式(6)和(7)确定：(6)(7)式(7)是推广的H-H迹长估计方法，限于篇幅，推导过程在此从略。与H-H方法不同的是，对于给定一组统计样本，平均迹长不仅与窗口尺度有关，而且还与结构面在窗口内的视倾角有关，为结构面与窗口底边的夹角，详见图1所示。式(5)的物理意义是：交切窗口结构面代表了长迹长结构面，迹长估计需要用概率方法解决；两端可见结构面代表了短迹长结构面，迹长估计问题已经成为确定性问题，可以利用统计方法直接求取。式(5)虽然与H-H算法密切相关，但是其实质上已经不是单纯的概率模型方法，而是概率模型方法与确定模型方法的有机结合，所以称式(5)式(7)为广义H-H迹长估计公式。由于l2不可能为零，从而避免了当n2接近N时，H-H迹长估计值接近于零的情形。对于无限长窗口，由于忽略结构面与取样窗口两侧边的交切关系，所以式(5)可简化为如下形式：l=1/N2-1/cos1-cos2(n1+2n0)h+(8)3.2广义H-H迹长估计方法检验任何一种算法除了具有理论上的精确性、严密性和完备性外，更重要的是还必须接受实践的检验，因此如何验证广义H-H迹长估计算法的正确性是一个关键问题。按照圆盘理论，结构面走向方向和倾向方向迹长大致相当，那么侧裂结构面走向方向的全迹长平均值与其倾向方向平均迹长估计值应当相等，据此可以验证迹长估计的正确性。表2侧裂结构面估计迹长和全迹长对照基于概率模型的估计迹长基于确定模型的平均迹长硐号岩性岩性段长/m平均迹长/m硐号岩性岩性段长/m平均迹长/mPD45-4角砾熔岩2633.42PD45-3角砾熔岩1453.40PD45-5角砾熔岩693.49PD75角砾熔岩1811.85PD75-3角砾熔岩241.49PD75-1角砾熔岩531.82PD87交斑状玄武岩3301.81PD87-1斑状玄武岩1201.70PD87-3斑状玄武岩1201.24PD87-2斑状玄武岩1401.35PD18-5玄武岩922.68PD18-3玄武岩5042.84PD18-6玄武岩1402.76PD56玄武岩1081.65PD56-1玄武岩211.85PD76玄武岩1531.98PD76交玄武岩1481.52为此，采用式(8)对各勘探平硐侧裂结构面迹长进行计算，考虑到实际调查中侧裂结构面倾角大于60，故取1、2分别为/3和2/3，计算结果如表2所列。由表2可见，除PD76差别较大(0.46m)外，多数不超过0.3m，估计迹长与全迹长总体上比较接近，由此可以证明，广义H-H迹长估计算法是正确的、可靠的和实用的。3.3广义H-H迹长估计方法应用广义H-H迹长估计方法可以推广到坝肩抗力体部位，在实际计算时，仍然取1、2分别为/3和2/3，利用式(8)对坝肩抗力体各勘探平硐侧裂结构面迹长进行计算。为了便于各种迹长估计方法之间相互对比，将计算结果列入表1中。通过对比发现，其它迹长算法与广义H-H方法有如下关系：(1)当迹长短小时，广义H-H方法计算结果与Laslett接近，2L0.5法和H-H方法迹长估计结果略微偏高，但无严格的对应关系。(2)当迹长较大时，Laslett和H-H方法计算结果显著偏大，而2L0.5法略微偏小，例如PD45-4和PD45-5，只有广义H-H方法真实反映了侧裂结构面的实际迹长。(3)当n2接近N时，H-H方法完全失效(例如PD69)；多数情况下H-H方法迹长估计值较广义H-H方法偏大，主要原因是H-H迹长估计与结构面在窗口内的视倾角无关。(4)2L