花崗巖孤石統(tǒng)計方法及其在樁基工程中的應用

張 曉，曾 強，李 勇，張永定

(1.廣東省電力設計研究院，廣東 廣州 510663；2.中山大學，廣東 廣州 510275)

1 概述

在花崗巖地區(qū)進行嵌巖樁施工時，往往難以區(qū)分孤石和基巖面，因此花崗巖孤石是影響成樁質量的重要因素。現(xiàn)行規(guī)程規(guī)范并無花崗巖孤石的評價方法，因此研究花崗巖孤石的評價方法具有必要性和緊迫性。

與孤石類似的隱伏地質體為巖溶，《火力發(fā)電廠巖土工程勘測技術規(guī)程》(DL/T 5074)以巖溶點密度、鉆孔遇洞率、鉆孔線溶率評價巖溶發(fā)育程度。借鑒巖溶發(fā)育程度的評價方法，類似地提出“遇孤石率”和“線孤石率”，這兩個統(tǒng)計量，以孤石在水平向和垂向的分布特性來評價孤石發(fā)育程度，針對超前鉆所關心的入巖深度問題，本文對孤石厚度進行了統(tǒng)計。

由于孤石之間是相互獨立的，其水平向、豎向分布特性具有一定的統(tǒng)計意義。本文根據(jù)某電廠一期工程勘測成果進行反分析，通過概率統(tǒng)計原理，分析水平向遇孤石率和豎向孤石厚度、線孤石率的分布類型，估算孤石厚度的置信區(qū)間，并將統(tǒng)計成果與超前鉆資料進行比對以檢驗評價方法的適宜性。

2 水平向分布特性——遇孤石率的統(tǒng)計方法

遇孤石率是指孤石在水平方向的分布特性，為考察樣本數(shù)對統(tǒng)計精度的影響，以下對工程場地進行了分區(qū)統(tǒng)計和合并統(tǒng)計。

2.1 分區(qū)統(tǒng)計

以隨機變量X表示鉆孔是否遇孤石，遇孤石為1，未遇孤石為0，則單個鉆孔可視為貝努利試驗，遇孤石的事件服從概率為p的兩點分布?？赏ㄟ^區(qū)內(nèi)各鉆孔的Xi值估計參數(shù)

區(qū)內(nèi)的鉆孔相當于n重貝努利試驗，因此遇孤石的孔數(shù)x服從概率為p的二項分布如式1：

2.2 全部統(tǒng)計區(qū)

2.3 統(tǒng)計方法的應用和驗證

將煤場區(qū)劃分為若干個邊長300m的正方形區(qū)域，各分區(qū)遇孤石率的統(tǒng)計以及檢驗結果見表1。對整個煤場區(qū)遇孤石率的統(tǒng)計以及檢驗結果見表2。

表1 煤場區(qū)分區(qū)統(tǒng)計結果

表2 煤場區(qū)、主廠區(qū)統(tǒng)計結果

由表1，各分區(qū)k值與預測值的偏離較大；即使在鉆孔數(shù)量最多分區(qū)，其施工圖階段的估計量與超前鉆階段算得的遇孤石率k/n仍存在較大誤差；以顯著水平α=0.05檢驗假設p=，多個分區(qū)檢驗結果落在拒絕域。雖然第一個分區(qū)在大樣本(n＞30)的情況下，估計量與超前鉆k/n值仍存在一定的偏差，說明區(qū)n＜30時，估計量與p真值的偏差較大，不宜用其數(shù)學期望預測遇孤石孔數(shù)。

3 豎向分布特性——孤石揭露厚度、線孤石率的統(tǒng)計方法

3.1 統(tǒng)計量的選定及其概率分布

(1)孤石厚度

孤石厚度在本文有兩個統(tǒng)計量，首先是單孔揭露孤石的總厚度H，H反映單個鉆孔揭露多個孤石的情況；其次是單個孤石的揭露厚度H'。

孤石厚度H密度函數(shù)的形式未知，可通過核密度估計獲得任意h的總體密度函式中n為樣本數(shù)；K為某個分布的密度函數(shù)，不同的K對核密度估計影響不大，因此采用標準正態(tài)分布；L為區(qū)間長度，稱為帶寬，其取值影響密度函數(shù)曲線的光滑程度，K為標準正態(tài)分布時最佳帶寬為L=1.0Sn-0.2(S為樣本標準差)。

根據(jù)該電廠一期工程超前鉆階段的H、H'統(tǒng)計成果，作密度函數(shù)、分布函數(shù)的核密度估計曲線，與理論正態(tài)分布進行對比見圖1～圖2。

圖1 單孔孤石總厚度H密度函數(shù)f(H)和分布函數(shù)F(H)曲線

圖2 單個孤石厚度H'密度函數(shù)f(H')和分布函數(shù)F(H')曲線

H、H'均受到節(jié)理裂隙、風化程度、地形地貌以及鉆孔與孤石相對位置等多種因素影響，而每一個影響因素都不起決定作用，且這些影響因素是可以疊加的，根據(jù)概率論原理，可認為H、H'近似服從正態(tài)分布。圖1～圖2顯示，H、H'核密度估計曲線適宜用嚴格非負的對數(shù)正態(tài)分布來描述[3]，對數(shù)正態(tài)分布密度函數(shù)見式3：

為了簡化分析，擬采用正態(tài)分布近似地描述H、H'的分布規(guī)律。

下面對H、H'進行χ2擬合優(yōu)度檢驗，作假設檢驗1：H服從正態(tài)分布N(, S2H)，假設檢驗2：H'服從正態(tài)分布N(, S2H')，顯著水平α均為0.05，結果見表3、表4。

表3 H分布函數(shù)的假設檢驗

表4 H'分布函數(shù)的假設檢驗

χ2擬合優(yōu)度檢驗結果顯示H、H'均服從正態(tài)分布。

通過上述分析，得到以下認識：

①H、H'服從對數(shù)正態(tài)分布，近似服從正態(tài)分布；

②分布函數(shù)曲線圖中正態(tài)分布曲線、經(jīng)驗分布曲線、核密度估計的分布曲線交點在F軸0.25、0.88～0.90點，即置信度90%的正態(tài)分布單側界限值與實測值較接近，在F∈(0.25, 0.90)的區(qū)間內(nèi)，正態(tài)分布的估值大于實測值，超出這個區(qū)間，正態(tài)分布的估值就小于實測值。

(2)孤石所占的線性比例

以下簡稱“線孤石率”或y，把y定義為孤石總厚度H與基巖面以上風化帶厚度T(含孤石)的比值，y的χ2及Kolmogorov-Smirnov擬合優(yōu)度檢驗顯示y不服從正態(tài)分布，根據(jù)超前鉆階段的y統(tǒng)計成果，作密度函數(shù)、分布函數(shù)的核密度估計曲線，與理論正態(tài)分布、指數(shù)正態(tài)分布進行對比見圖3。

圖3 線孤石率y密度函數(shù)f(y)和分布函數(shù)F(y)曲線

由圖3，y核密度估計曲線與對數(shù)正態(tài)分布曲線擬合得非常好，可知y服從指數(shù)正態(tài)分布，其分布函數(shù)如式3，通過y數(shù)據(jù)樣本進行參數(shù)估計，式3中μ=2.5827，σ=0.9319。

據(jù)相關的研究，正態(tài)、對數(shù)正態(tài)分布特性的區(qū)別唯一受控于變異系數(shù)，變異系數(shù)＜0.3時兩者差異不大。由于y的變異系數(shù)達到了0.95，因此正態(tài)、指數(shù)正態(tài)分布曲線存在一定的誤差。

由圖3，y正態(tài)分布曲線、經(jīng)驗分布曲線、核密度估計的分布曲線交點在F(y)軸0.25、0.90兩點，即置信度90%的正態(tài)分布單側界限值與實測值較接近；在F(y)∈(0.25, 0.90)的區(qū)間內(nèi)，正態(tài)分布的y估值大于實測值；超出這個區(qū)間，正態(tài)分布的y估值就小于實測值，如F(y)=0.95時正態(tài)分布y估值與實際相差達0.075。

綜上，在F(y)∈(0.25, 0.90)的區(qū)間內(nèi)以正態(tài)分布近似地描述線孤石率y是較保守的、可行的，為簡化分析，仍以正態(tài)分布近似地描述y的分布規(guī)律。

3.2 統(tǒng)計方法及實例

超前鉆入巖深度應達到多少，才能以較大的把握確定孔底是落在基巖，而不是落在孤石上呢？這個問題可以通過H'的均值μ、樣本標準差S，以及μ的置信區(qū)間來解答。

若H'服從理論正態(tài)分布，則H'應＜μ+σ(單側置信度88%)，由3.1節(jié)的討論可知，以施工圖資料估計的存在一定誤差，不宜直接使用，但H'仍可以其均值的置信區(qū)間來描述。

由于總體方差σ2未知，求H '均值的置信區(qū)間，首先作標準化變換其中則統(tǒng)計量T服從t(n-1)分布；根據(jù)3.1節(jié)的討論，取單側顯著水平α=0.1，單側置信區(qū)間界限值為標準值，表示H'的均值有90%的可能落在單側置信區(qū)間標準值是H'的可靠性估值，上述方法與《巖土工程勘察規(guī)范》中標準值的計算方法相近，不同的是規(guī)范中顯著水平取0.05，學生氏函數(shù)界限值采用擬合公式計算，而本文顯著水平取0.1，學生氏函數(shù)界限值t1-α是直接計算而得。

同理可得H、y的統(tǒng)計方法。

對煤場區(qū)的H'、H、y分階段統(tǒng)計見表5。

表5 煤場區(qū)H'、H、y分階段統(tǒng)計

4 結論

(1)單個鉆孔是否遇孤石的事件可視為努利試驗，而多個鉆孔遇孤石孔數(shù)x服從遇孤石概率為p的二項分布，其數(shù)學期望np可用于預測超前鉆/樁基施工遇孤石的孔/樁數(shù)。

(2)孤石厚度、線孤石率服從對數(shù)正態(tài)分布，為簡化分析，可近似地以正態(tài)分布進行統(tǒng)計，通過其平均值，樣本標準差，以及可靠度估值，可為超前鉆及樁基施工提供依據(jù)。

[1]梁飛豹，等.應用統(tǒng)計方法[M].北京：北京大學出版社，2010.

[2]張繼周，繆林昌.巖土參數(shù)概率分布類型及其選擇標準[J].巖石力學與工程學報，2009，28(S2).

[3]復旦大學.概率論第二冊[M].北京：高等教育出版社，1979.