(1.西安理工大学土木博士后流动站；2.长安大学地测学院.成都理工学院工程地质研究所)

作者简介：范留明(1968-)，河南修武人，男，地质工程博士，研究方向为岩体稳定。

在水电工程中，侧裂面是控制拱坝坝肩抗滑稳定的一个重要力学边界。如果坝肩部位没有明显的陡倾断层通过，则可能构成侧裂面的结构面只能是某一走向的陡倾基体裂隙，作者将这组陡倾裂隙称之为侧裂结构面。侧裂结构面迹长是一个重要岩体结构参数，它直接反映了结构面的几何大小，所以侧裂结构面迹长是侧裂面边界条件分析中一项重要内容。但是，由于受野外露头的影响，可测量到的全迹长结构面样本数量往往很少，利用确定模型方法直接求柔构面平均迹长受到了很大的限制，因此基于概率模型的迹长估计方法便成为获柔构面平均迹长的一个主要途径。其基本原理是根据取样窗口中结构面可见迹长及其与窗口(或测线)的交切关系来推求平均迹长。尽管现有迹长估计方法较多，但是专门用于侧裂结构面迹长估计方法尚没见报道。在对西南某大型水电站工程侧裂面边界条件研究中发现，在某些情况下，各种迹长估计算法之间存在较大差别，迹长估计结果与实际情况出入很大。作者对这一现象进行了较深入的分析，并在此基础上，提出了广义H-H结构面迹长估计公式，理论和实践证明：该算法是正确可靠的。

1 结构面迹长估计常用方法

目前已有的迹长估计方法很多，常用的方法有：Kulatilake迹长估计法、Cruden迹长估计法、Laslett(1982)迹长估计法、H-H迹长估计法和倍半迹长方法(2L_0.5)等等。其中Kulatilake算法要求由倾向和倾角两个独立变量定义的结构面二维概率密度函数是已知的，所以实用性较差。Cruden公式假设迹长服从负指数分布，然而众多学者则认为迹长服从对数正态分布。鉴于上述理由，在实际计算中，主* 捎昧薍-H、Laslett和2L_0.5三种迹长估算方法。

Laslett基于极大似然原理，推导了由窗口数据估计结构面平均迹长的极大似然公式：

式中：X_i是第i条两端可见结构面可见迹长；Y_j是第j条一端可见结构面可见迹长；Z_k是第k条两端均不可见结构面可见迹长；n为两端可见结构面条数；m为一端可见结构面条数；p为两端均不可见结构面条数。

特别地，当两端可见结构面条数n和两端均不可见结构面条数p均为零时，上式变为简化为如下形式：

式中：n₀是两端均不可见的结构面条数；n₁是一端可见另一端不可见的结构面条数；n₂是两端可见的结构面条数；N是结构面总数，由下式确定：

2 侧裂结构面迹长估计

2.1 问题特点 拟建西南某大型水电站是一巨型规模的水力发电枢纽工程，大坝为双曲高拱坝，坝区侧裂结构面主要由随机分布的陡倾基体裂隙构成。按照工程设计和抗滑稳定性要求，确定左岸侧裂结构面走向变化范围为N5°W～N35°W，右岸为N75°E～N75°W，倾角均大于60°。侧裂结构面的野外现场调查主要在勘探平硐中进行，但是由于坝区绝大多数平硐沿垂直河谷或平行河谷方向施工，且硐高仅为2m，所以平硐布置方向和硐高皆不利于侧裂面迹长调查(图2).在坝肩抗力体局部地段，尤其是在距离坝肩抗力体上游约200m处的厂房区，布置了少数与侧裂面走向近平行的勘探平硐，在其硐顶可以获得一部分全迹长资料(图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方法，而2L_0.5方法的迹长估计值一般不超过2倍的硐高(4m).例如厂房PD13、PD45-4和PD45-5，H-H与Laslett迹长估计结果显然与实际情况不相符合。(3)H-H方法具有一定的“不稳定性”，主要表现在两个方面：一方面，当n₂接近N时，其迹长估计值明显偏小，例如PD69平硐，H-H公式计算迹长值尚不足2m；另一方面，当n₀较大时，其迹长估计值又明显增大，例如PD2、PD13、PD45-4和PD45-5等。

表1 各种算法估计迹长对比

3 广义H-H迹长估计法及其应用

3.1 广义H-H迹长估计方法 为了克服上述各种迹长估计方法的缺点，作者研究了一种迹长估计新方法。此方法能够正确解决Laslett对长迹长估计的局限性，有效消除H-H方法的不稳定性，理论上较其它迹长估计算法更完善、更成熟、更实用。其基本思想是：将结构面分成两组，其中一组由交切窗口的结构面组成，平均迹长lh按照推广的H-H迹长估计方法计算获得；另一组为两端可见结构面，其平均迹长l2等于全迹长均值，具体算法如下：

式中：l₂和l_h由式(6)和(7)确定：

式(7)是推广的H-H迹长估计方法，限于篇幅，推导过程在此从略。与H-H方法不同的是，对于给定一组统计样本，平均迹长不仅与窗口尺度有关，而且还与结构面在窗口内的视倾角φ有关，φ为结构面与窗口底边的夹角，详见图1所示。

式(5)的物理意义是：交切窗口结构面代表了长迹长结构面，迹长估计需要用概率方法解决；两端可见结构面代表了短迹长结构面，迹长估计问题已经成为确定性问题，可以利用统计方法直接求取。式(5)虽然与H-H算法密切相关，但是其实质上已经不是单纯的概率模型方法，而是概率模型方法与确定模型方法的有机结合，所以称式(5)～式(7)为广义H-H迹长估计公式。由于l2不可能为零，从而避免了当n2接近N时，H-H迹长估计值接近于零的情形。对于无限长窗口，由于忽略结构面与取样窗口两侧边的交切关系，所以式(5)可简化为如下形式：

3.2 广义H-H迹长估计方法检验 任何一种算法除了具有理论上的精确性、严密性和完备性外，更重要的是还必须接受实践的检验，因此如何验证广义H-H迹长估计算法的正确性是一个关键问题。按照圆盘理论，结构面走向方向和倾向方向迹长大致相当，那么侧裂结构面走向方向的全迹长平均值与其倾向方向平均迹长估计值应当相等，据此可以验证迹长估计的正确性。

表2 侧裂结构面估计迹长和全迹长对照

为此，采用式(8)对各勘探平硐侧裂结构面迹长进行计算，考虑到实际调查中侧裂结构面倾角大于60°，故取φ₁、φ₂分别为π/3和2π/3，计算结果如表2所列。由表2可见，除PD76差别较大(0.46m)外，多数不超过0.3m，估计迹长与全迹长总体上比较接近，由此可以证明，广义H-H迹长估计算法是正确的、可靠的和实用的。

3.3 广义H-H迹长估计方法应用 广义H-H迹长估计方法可以推广到坝肩抗力体部位，在实际计算时，仍然取φ₁、φ₂分别为π/3和2π/3，利用式(8)对坝肩抗力体各勘探平硐侧裂结构面迹长进行计算。为了便于各种迹长估计方法之间相互对比，将计算结果列入表1中。通过对比发现，其它迹长算法与广义H-H方法有如下关系：(1)当迹长短小时，广义H-H方法计算结果与Laslett接近，2L_0.5法和H-H方法迹长估计结果略微偏高，但无严格的对应关系。(2)当迹长较大时，Laslett和H-H方法计算结果显著偏大，而2L_0.5法略微偏小，例如PD45-4和PD45-5，只有广义H-H方法真实反映了侧裂结构面的实际迹长。(3)当n2接近N时，H-H方法完全失效(例如PD69)；多数情况下H-H方法迹长估计值较广义H-H方法偏大，主要原因是H-H迹长估计与结构面在窗口内的视倾角φ无关。(4)2L_0.5是Laslett(1982)公式在n1=N时的一种特殊情形，通常情况下，其估计迹长值与n₂和n₀的相对数量有关，且与n₂成正相关，即当n₂较n₀大，则估计迹长较广义H-H估计迹长偏高；反之，当n₂较n₀小时，估计迹长较广义H-H估计迹长偏低，这表明2L_0.5迹长估计属于有偏估计。

4 结 语

迹长是岩体结构面定量描述的重要指标之一，迹长估计也是岩体结构研究中的一项重要内容。但是由于现有的各种迹长估计算法不能够全面反映结构面(视)倾角与迹长的关系，所以无法正确估计侧裂结构面迹长。作者在对西南某大型水电站工程侧裂面边界条件研究中，提出了一种适用于侧裂结构面迹长估计的新方法，即广义H-H迹长估计算法，并用之于该水电工程的侧裂结构面迹长研究中，取得了很好的应用效果。

参 考 文 献：

[1] 陈剑平，肖树芳，王清。随机不连续面三维网络计算机模拟原理[M]。长春：东北师范大学出版社，1995.

[2] 范留明，柴贺军，黄润秋。岩体结构二维可视化模型及其应用研究[J]。地质灾害与环境保护，2000，11(4)：341～345.

[3] 黄国明。节理岩体及其工程应用[D]。成都：成都理工学院，1999.

[4] 黄润秋，胡卸文，陶连金，等。金沙江溪落渡水电工程岩体结构模型及其工程应用研究[C]。成都：成都理工学院，1997.