Page 26 - 2023年第54卷第1期
P. 26

件 A的条件概率(或称为试验概率,对应于某一具体工程的试验数据);P(B?A)为事件 B(对应于统
                                                                                     i           i
              计值)的后验概率。如果已知大样本的先验分布 P(B)和试验概率 P(A?B),按式(8)和式(9)就可以
                                                                                 i
                                                              i
              对抗剪强度参数的统计结果做出新的估计(即图 14中的 “统计值”):
                                                         k
                                                     μ b ∑  P(B?A) μ i                                  (8)
                                                       =
                                                               i
                                                        i =1
                                                      k
                                                   2
                                                    =
                                                 σ b ∑  ( μ i μ b  2  i                                 (9)
                                                           - )P(B?A)
                                                     i =1
                                                               2
              式中:k为数据分段数;μ i          为数据分段中间值;μ b         、σ b 分别为 Bayes方法估算的均值和标准差。
                  下面通过一个实例来说明 Bayes方法的具体应用情况。某工程坝基辉绿岩弱风化带中有一组缓倾
              角夹泥裂隙,是控制坝基抗滑稳定的主要结构面。为此,进行了 9组现场裂隙面抗剪试验,得出黏聚
                                  5
              力 c的均值为 5.13 × 10 Pa,从经验判断此值明显偏大                 [4] 。现借助 Bayes公式重新进行估算,根据同类
              岩体裂隙面抗剪试验 26组的试验成果和该工程的 9组试验成果,计算得到的后验概率分布列于表 5
              中。黏聚力 c的大概率估计值为:c = 0.58 × l.5 + 0.18 × 4.1 + 0.13 × 6.7 + 0.09 × 9.3 + 0.02 × 11.9 + 0 × 14.5 = 3.577 ×
                5
              10 Pa。
                                           表 5 抗剪强度参数 c的 Bayes公式计算结果
                     c值分段           0.2~2.8    2.8~5.4    5.4~8      8~10.6   10.6~13.2  13.2~15.8   总和
                           5
                  分段中值× 10?Pa         1.5       4.1        6.7        9.3       11.9       14.5
                    后验概率             0.58       0.18       0.13       0.09      0.02        0
                          5
                   分段值× 10?Pa        0.867      0.73       0.89      0.828      0.262       0        3.577

                  (2) “标准值” 和 “设计值” 的确定方法。由 Bayes方法求得的统计值,虽然反映了抗剪强度参
              数大样本的统计特性,但并不能直接应用于工程设计中,譬如假定某一参数服从正态分布,其均值为
              0.56,当标准差为 0.08时,该参数小于均值的概率 P(x< μ )可达 50%,意味着采用该参数直接进入设
              计中会承担很大的风险。因此,“统计值” 还需通过取概率分位数予以折减后变为 “标准值”,再除
              以分项系数后才能作为设计采用值(“设计值”)。
                  考虑到岩体结构的复杂性和工程实际应用需求,抗剪强度参数的 “标准值” 可采用概率分布的
              0.2分位值(相应的概率保证率为 80%)。对应于不同的分布概型,“标准值” 可以按下面的公式计算。
                  当抗剪强度参数服从正态分布时:
                                                                                                       (10)
                                                        μ k μ b
                                                          = + K σ b
                  当抗剪强度参数服从对数正态分布时:
                                                     μ b
                                                 =                    2 ))                             (11)
                                                               槡
                                               μ k       exp(K ln(1 + δ b
                                                        2
                                                   槡
                                                    1 + δ b
                                   = ? ;K为与概率分位数有关的系数。实际工程中一般取概率分位数为 0.2,
              式中:μ k  为标准值;δ b σ b μ b
              此时 K =- 0 .84。
                  如前所述,标准值仍然不能直接用于实际工程设计中,还必须要采用适当的分项系数 γ进行折减
              后才能作为工程 “设计值”,具体计算公式如下:
                                                             μ k
                                                            =                                          (12)
                                                          μ d
                                                              γ
                  当抗剪强度参数服从正态分布时:
                                                               2
                                                           1 + Kδ b
                                                        γ =                                            (13)
                                                           1 + Kδ b
                                                               1
                  当抗剪强度参数服从对数正态分布时:
                                                             1
                                                    γ =                                                (14)
                                                                    ]
                                                       exp [(1 - K) δ b
                                                                2
                                                                                                 —  2 1 —
   21   22   23   24   25   26   27   28   29   30   31