砂礫石料顆粒級(jí)配參數(shù)的多維威布爾分布模型

裴彥青，李建鋒，趙宇飛,劉 彪,李 銳

(1.新疆新華葉爾羌河流域水利水電開(kāi)發(fā)有限公司，新疆 喀什 844000;2.中國(guó)水利水電科學(xué)研究院，北京 100048;)

1 研究背景

砂礫石料作為一種無(wú)黏性的粗粒土材料，具有較優(yōu)的壓實(shí)性能、填筑密度大、承載力高、抗剪強(qiáng)度高、沉陷變形小等諸多優(yōu)良的工程特性[1- 4]。在中國(guó)地域遼闊的南疆地區(qū)，河床和戈壁沙漠中廣泛分布著儲(chǔ)量豐富的天然砂礫石，在當(dāng)?shù)厮婇_(kāi)發(fā)中，砂礫石壩也因享有就地取材、結(jié)構(gòu)簡(jiǎn)單、變形適應(yīng)性好、施工工序少等工程優(yōu)勢(shì)而成為當(dāng)?shù)刈畛Ｒ?jiàn)的壩型之一[5]。在大壩填筑碾壓過(guò)程中，掌握壩料的物理力學(xué)特性是有效控制土石壩填筑施工質(zhì)量的關(guān)鍵，對(duì)壩體質(zhì)量控制以及壩體安全穩(wěn)定運(yùn)行具有非常重要的意義。

在土石壩填筑施工過(guò)程中，不同粒徑壩料按照一定的比例混合起來(lái)，通過(guò)不同的碾壓工藝，使其達(dá)到符合要求的密實(shí)度，以滿足上部結(jié)構(gòu)對(duì)填方的要求，土石壩料的級(jí)配特征參數(shù)直接影響壩料填筑的壓實(shí)特性；而壩體材料的滲透特性主要是由孔隙結(jié)構(gòu)特征決定，級(jí)配是影響壩料孔隙結(jié)構(gòu)特征的主要因素。因此，研究砂礫石料在一定壓實(shí)度條件下的物理力學(xué)特性(壓實(shí)特性和滲透特性等)，首先應(yīng)對(duì)砂礫石料顆粒級(jí)配進(jìn)行深入研究[6]。從國(guó)內(nèi)外相關(guān)文獻(xiàn)來(lái)看，不少學(xué)者通過(guò)室內(nèi)試驗(yàn)或數(shù)值仿真等手段也開(kāi)展了相關(guān)研究，均指出顆粒級(jí)配對(duì)砂礫石料的壓實(shí)特性和滲透特性具有重要影響[7- 16]。如在級(jí)配表征壓實(shí)特性方面，李曉柱等[11]依托糯扎渡現(xiàn)場(chǎng)試驗(yàn)情況，引入了離散元數(shù)值方法來(lái)研究堆石體的碾壓特性，探討了粒徑分布對(duì)碾壓質(zhì)量的影響，并提出了最佳級(jí)配范圍；朱晟等[12]結(jié)合室內(nèi)試驗(yàn)和現(xiàn)場(chǎng)試驗(yàn)開(kāi)展了堆石壩料級(jí)配特征對(duì)壓實(shí)干密度的影響，并建立了具有分形分布的堆石料級(jí)配公式，證實(shí)了級(jí)配是影響堆石壩料壓實(shí)性能的主要因素。在級(jí)配表征滲透特性方面，王雙等[13]采用神經(jīng)網(wǎng)絡(luò)研究了碎石土級(jí)配對(duì)滲透系數(shù)的影響，研究發(fā)現(xiàn)P5含量(小于5mm的粒徑)[11]、不均勻系數(shù)或曲率系數(shù)與滲透系數(shù)之間存在較強(qiáng)的相關(guān)性；張國(guó)棟等[14]通過(guò)室內(nèi)試驗(yàn)研究了顆粒級(jí)配對(duì)滲透特性的影響，研究發(fā)現(xiàn)特征粒徑、不均勻系數(shù)和曲率系數(shù)與滲透系數(shù)之間的關(guān)系符合線性擬合關(guān)系。

然而現(xiàn)有的相關(guān)研究基本上是基于室內(nèi)試驗(yàn)或數(shù)值仿真等手段開(kāi)展的研究，通過(guò)設(shè)計(jì)模型試驗(yàn)獲取數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析的方法，由于試驗(yàn)中數(shù)據(jù)量較少且未能考慮施工中的不確定性因素的影響，試驗(yàn)數(shù)據(jù)往往較為理想，所得方法用于擬合實(shí)際工程中的問(wèn)題時(shí)通常難以較好地描述。此外相關(guān)研究?jī)H考慮了砂礫石料顆粒級(jí)配中的某幾個(gè)特征，未全面、細(xì)致考慮顆粒級(jí)配的變化，導(dǎo)致顆粒級(jí)配的代表性不足，雖然反映出顆粒級(jí)配對(duì)土的壓實(shí)性和滲透特性存在影響，但是獲得的相關(guān)影響規(guī)律尚不夠深入[17]。當(dāng)前試驗(yàn)方法仍然是確定其參數(shù)最直接、最可靠的方法，也是大多數(shù)理論和經(jīng)驗(yàn)估值的重要依據(jù)[15]。充分利用有限的試驗(yàn)數(shù)據(jù)，結(jié)合概率統(tǒng)計(jì)方法探索顆粒級(jí)配與力學(xué)參數(shù)間的關(guān)系，是目前較為實(shí)用的途徑。

根據(jù)已有的樣本數(shù)據(jù)去推斷總體分布時(shí)，一般需要將數(shù)據(jù)擬合成具體的概率分布模型，并估計(jì)其相應(yīng)的分布參數(shù)，在數(shù)理統(tǒng)計(jì)分析中威布爾分析模型對(duì)各種類型試驗(yàn)數(shù)據(jù)有極強(qiáng)的適應(yīng)能力，不僅能夠適用于大樣本數(shù)據(jù)，對(duì)小樣本數(shù)據(jù)也有較好的適用性[18]，因此本文選用威布爾分布作為壩料級(jí)配特征參數(shù)的總體分布模型。在本研究中，依托新疆阿爾塔什水利樞紐工程施工過(guò)程中獲取的大量檢測(cè)試驗(yàn)數(shù)據(jù)，在大量統(tǒng)計(jì)的基礎(chǔ)上，通過(guò)分析各參數(shù)與干密度和滲透系數(shù)之間的相關(guān)性，選取了與干密度和滲透系數(shù)相關(guān)性最密切的3個(gè)級(jí)配特征參數(shù)，并對(duì)其分別構(gòu)建了一維威布爾分布模型，求出對(duì)應(yīng)模型的3個(gè)控制參數(shù)：位置參數(shù)μ、尺度參數(shù)σ、形狀參數(shù)ξ，得到了在一定保證率條件下的砂礫石料級(jí)配參數(shù)估計(jì)值。在此基礎(chǔ)上，進(jìn)一步推導(dǎo)出了三維威布爾分布函數(shù)表達(dá)式，建立了能夠利用3個(gè)級(jí)配特征參數(shù)間相關(guān)關(guān)系反映砂礫石填料區(qū)干密度和滲透系數(shù)的三維威布爾分布模型，使干密度和滲透系數(shù)的表征方法更具綜合性，可用于粗粒料的質(zhì)量評(píng)估和預(yù)測(cè)。

2 多維威布爾分布模型的建立

2.1 一維威布爾分布模型的建立

在概率統(tǒng)計(jì)分析方面，威布爾分析模型能夠較好的擬合數(shù)據(jù)樣本點(diǎn)，適用于完全樣本、有偏樣本和截尾樣本的擬合，對(duì)各種類型試驗(yàn)數(shù)據(jù)有極強(qiáng)的適應(yīng)能力。三參數(shù)威布爾概率分布函數(shù)的表達(dá)式可表示為[19- 22]：

(1)

其密度函數(shù)為：

(2)

式中，F(xiàn)(x)—分布函數(shù)；f(x)—密度函數(shù)；x—隨機(jī)變量；ξ—形狀參數(shù)；σ—尺度參數(shù)；μ—位置參數(shù)。

威布爾分布參數(shù)求解常用的方法是最小二乘法[23- 25]，本文選用最小二乘法進(jìn)行參數(shù)計(jì)算，初始計(jì)算時(shí)μ設(shè)定為任意小于樣本數(shù)據(jù)最小值的數(shù)，進(jìn)而可求出形狀參數(shù)ξ和尺度參數(shù)σ：

(3)

此外在實(shí)際應(yīng)用時(shí)，往往需要給出在一定可靠度條件下的參數(shù)估計(jì)值，設(shè)參數(shù)估計(jì)的可靠度為R2，則一維威布爾分布在可靠度R2條件下的估計(jì)值為[27]：

2.2 三維威布爾分布模型的建立

實(shí)際工程中壩料級(jí)配的描述，是一系列顆粒粒徑及其含量的表達(dá)，這些參數(shù)在一定的尺度范圍內(nèi)存在著隨機(jī)性和符合統(tǒng)計(jì)規(guī)律的自相似性，且某些參數(shù)間存在著明顯的相關(guān)性[28]。因此，難以采用一維概率分布進(jìn)行全面描述，但目前多維概率分布中參數(shù)、相關(guān)系數(shù)確定以及假設(shè)檢驗(yàn)仍存在困難，限制了多維概率分布在壩料級(jí)配特性表達(dá)中的應(yīng)用。

近年來(lái)不少研究者逐步完善了多維復(fù)合極值分布理論[29- 31]，本研究受上述離散型多維復(fù)合極值分布模型建立思想的啟發(fā)[32- 35]，在一維威布爾分布的基礎(chǔ)上，構(gòu)建了基于威布爾分布的多維復(fù)合分布函數(shù)模型。對(duì)于三維威布爾分布而言，其概率分布函數(shù)應(yīng)該是每個(gè)變量的邊緣函數(shù)的組合，考慮三維威布爾概率分布中變量之間的相關(guān)系數(shù)，通過(guò)公式推導(dǎo)與數(shù)值擬合，得出的三維威布爾分布模型概率函數(shù)的數(shù)學(xué)表達(dá)式為：

(5)

式中，ξj、μj、σj—x，y，z的邊緣分布形狀參數(shù)、位置參數(shù)和尺度參數(shù)，j=1，2，3；0≤α，β≤1為由相關(guān)系數(shù)r得到的相關(guān)參數(shù)，α，β均為1時(shí)，x，y，z相互獨(dú)立。史道濟(jì)[36]通過(guò)矩估計(jì)法，給出了相關(guān)參數(shù)α，β的顯式表達(dá)式：

(6)

2.3 擬合優(yōu)度檢驗(yàn)

在評(píng)估上述所獲得的概率分布模型的擬合精度，本文采用了K-S檢驗(yàn)(Kolmogorov-Smirnov)來(lái)檢驗(yàn)指定的威布爾分布和實(shí)際分布是否相符，具體檢驗(yàn)步驟如下[37]：

(1)建立原假設(shè)H0:G(x)=F(x)，其中F(x)為假設(shè)的分布函數(shù)，G(x)為樣本觀測(cè)值得累積分布函數(shù)：G(xi)=i/n，i=1，2，…，n。

(2)計(jì)算檢驗(yàn)統(tǒng)計(jì)量：Dn=max{|Fn(x)-F(x)|}。

(3)在不同顯著性水平α下，比較Dn與K-S檢驗(yàn)臨界值表中的臨界值D(n，α)。若Dn

3 工程應(yīng)用

3.1 工程背景

阿爾塔什水利樞紐工程位于塔里木河源流之一的葉爾羌河干流山區(qū)下游河段的新疆維吾爾自治區(qū)克孜勒蘇柯?tīng)柨俗巫灾沃莅⒖颂湛h庫(kù)斯拉甫鄉(xiāng)境內(nèi)，是一座在保證向塔里木河干流生態(tài)供水目標(biāo)的前提下，承擔(dān)防洪、灌溉、發(fā)電等綜合利用任務(wù)的大型骨干水利樞紐工程。水庫(kù)工程正常蓄水位為1820.0m，水庫(kù)設(shè)計(jì)洪水位為1821.62m，校核洪水位為1823.69m，總庫(kù)容22.49億m3；電站裝機(jī)容量755MW。阿爾塔什水利樞紐工程為大(1)型一等工程。

阿爾塔什水庫(kù)混凝土面板砂礫石堆石壩壩體橫剖面圖如圖1所示，壩軸線全長(zhǎng)795.0m，壩頂高程1825.80m，壩頂寬度為12m，最大壩高164.8m，上游壩坡采用1∶1.7，下游壩坡坡度為1∶1.6。面板壩直接建造于河床深厚覆蓋層上，覆蓋層最大厚度94m。大壩抗震設(shè)計(jì)烈度為9度，100年超越概率2%的設(shè)計(jì)地震動(dòng)峰值加速度為320.6g。壩體填筑分區(qū)從上游至下游分別為上游蓋重區(qū)、上游鋪蓋區(qū)、混凝土面板、墊層料區(qū)、過(guò)渡料區(qū)、砂礫料區(qū)、利用料區(qū)、爆破料區(qū)。

圖1 阿爾塔什混凝土面板堆石壩壩體橫剖面圖

3.2 影響壩料物理力學(xué)特性的級(jí)配參數(shù)確定

筑壩砂礫填料的物理力學(xué)特性(壓實(shí)特性和滲透特性)具有級(jí)配相關(guān)性，但在砂礫石料施工質(zhì)量控制中，目前主要以含礫量表征級(jí)配對(duì)物理力學(xué)特性的影響，未考慮級(jí)配形狀、最大粒徑、不均勻系數(shù)、曲率系數(shù)的影響。已有研究表明，砂礫石料壓實(shí)特性和滲透特性受這些因素的影響較為顯著，而簡(jiǎn)化的室內(nèi)實(shí)驗(yàn)確定的級(jí)配參數(shù)難以準(zhǔn)確描述現(xiàn)場(chǎng)原級(jí)配料的干密度和滲透系數(shù)。在此基礎(chǔ)上推導(dǎo)出來(lái)的外延公式缺乏普遍適用性，應(yīng)用時(shí)存在較大的誤差。為了研究這些參數(shù)對(duì)其敏感性，本研究借助數(shù)理統(tǒng)計(jì)的方法，基于施工過(guò)程中大量的檢測(cè)試驗(yàn)數(shù)據(jù)，建立試驗(yàn)參數(shù)間的關(guān)系，選擇與干密度和滲透系數(shù)相關(guān)性最緊密的3個(gè)因素，并對(duì)這3個(gè)因素通過(guò)威布爾分布進(jìn)行擬合分析，推導(dǎo)出在一定保證率條件下的估計(jì)值。

在阿爾塔什水利樞紐工程堆石壩工程施工過(guò)程中，本文收集整理了堆石壩工程檢測(cè)資料中砂礫石料挖坑質(zhì)量檢測(cè)試驗(yàn)數(shù)據(jù)504組和滲透試驗(yàn)數(shù)據(jù)45組，每組試驗(yàn)中的數(shù)據(jù)包含：干密度(g/cm3)、滲透系數(shù)(cm/s)、礫石含量(%)、含泥量(<0.075mm/%)、曲率系數(shù)Cc、不均勻系數(shù)Cu和最大粒徑(mm)。由于大壩填筑實(shí)際施工過(guò)程中實(shí)時(shí)采集到的數(shù)據(jù)受衛(wèi)星點(diǎn)位、地形地貌等的影響，不可避免的會(huì)產(chǎn)生一定的異常值，需要將這些異常值從數(shù)據(jù)集中剔除。本研究采用箱形圖的四分位區(qū)間檢測(cè)異常值的方法，該方法從實(shí)際樣本數(shù)據(jù)出發(fā)，不限制數(shù)據(jù)服從某種特定的分布形式，可以比較簡(jiǎn)明直觀的反映數(shù)據(jù)真實(shí)的面貌，在異常值識(shí)別處理方面具有一定的優(yōu)越性。針對(duì)異常值數(shù)據(jù)本文采用剔除的方式處理，計(jì)算所得的干密度和滲透系數(shù)與其他級(jí)配特征參數(shù)的相關(guān)性r見(jiàn)表1，由表可知，與干密度相關(guān)性最密切的三個(gè)特征參數(shù)為：最大粒徑、礫石含量、曲率系數(shù)；與滲透系數(shù)相關(guān)性最密切的3個(gè)特征參數(shù)為：含泥量、礫石含量、曲率系數(shù)。這3個(gè)因素與干密度的相關(guān)性關(guān)系如圖2所示，限于篇幅，僅列出砂礫填筑料顆粒級(jí)配參數(shù)與干密度之間的關(guān)系。

圖2 砂礫石料顆粒級(jí)配特征參數(shù)相關(guān)散點(diǎn)圖

表1 砂礫石料物理力學(xué)特性相關(guān)性分析

3.3 基于一維威布爾分布的砂礫石級(jí)配參數(shù)擬合

根據(jù)第2.1節(jié)推導(dǎo)的威布爾分布擬合公式，編制相應(yīng)的程序?qū)εc干密度和滲透系數(shù)最為相關(guān)的級(jí)配特征參數(shù)逐個(gè)進(jìn)行擬合，進(jìn)而計(jì)算出威布爾分布的3個(gè)參數(shù)，構(gòu)建一維威布爾概率分布模型。分別對(duì)各參數(shù)的實(shí)際有效數(shù)據(jù)進(jìn)行一維威布爾分布擬合，限于篇幅，此處僅列出與干密度最為相關(guān)的級(jí)配特征參數(shù)的擬合結(jié)果如圖3—5所示，對(duì)應(yīng)的威布爾分布參數(shù)見(jiàn)表2—3，從擬合結(jié)果可以看出，威布爾分布擬合結(jié)果的相關(guān)系數(shù)r均超過(guò)0.98，擬合結(jié)果較為可觀。然后利用K-S檢驗(yàn)法分別對(duì)本文威布爾分布擬合得到的概率分布進(jìn)行檢驗(yàn)，K-S檢驗(yàn)法結(jié)果見(jiàn)表2—3所示，從表中可看出在在顯著性水平α=0.01條件下，威布爾分布擬合得到的各項(xiàng)級(jí)配參數(shù)的概率密度函數(shù)都通過(guò)了擬合良好性檢驗(yàn)，能夠真實(shí)反映隨機(jī)變量的概率分布情況。說(shuō)明本研究的威布爾分布模型能夠較好的適用于壩體填料試驗(yàn)數(shù)據(jù)的描述及參數(shù)估計(jì)，具有較好的應(yīng)用前景。

圖3 最大粒徑～維威布爾分布擬合

圖4 礫石含量～維威布爾分布擬合

圖5 曲率系數(shù)～維威布爾分布擬合

表2 表征壓實(shí)特性的級(jí)配威布爾分布模型參數(shù)

表3 表征滲透特性的級(jí)配威布爾分布模型參數(shù)

3.4 基于三維威布爾分布的砂礫石級(jí)配參數(shù)擬合

通過(guò)對(duì)干密度和滲透系數(shù)主要影響因素的相關(guān)性分析，得出了在砂礫石填料區(qū)與干密度和滲透系數(shù)各自相關(guān)性最密切的3個(gè)因素，但這3個(gè)參數(shù)間還同時(shí)存在著相關(guān)性。基于此，將這3個(gè)參數(shù)利用推導(dǎo)出的公式(5)建立三維威布爾分布模型。其中，此時(shí)的頻率分布不再是一維威布爾分布中的單個(gè)參數(shù)的頻率，而是同時(shí)滿足3個(gè)參數(shù)時(shí)的頻率。計(jì)算結(jié)果如圖5所示，圖中橫坐標(biāo)表示三維變量分組的個(gè)數(shù)，縱坐標(biāo)表示分布函數(shù)值，圖6(a)是砂礫石填料區(qū)基于最大粒徑、礫石含量、曲率系數(shù)間相關(guān)關(guān)系建立的三維威布爾分布模型，用于表征壩料的壓實(shí)特性；圖6(b)砂礫石填料區(qū)基于含泥量、礫石含量、曲率系數(shù)間相關(guān)關(guān)系建立的三維威布爾分布模型，用于表征壩料的滲透特性。由圖可知，三維威布爾分布較好地?cái)M合了試驗(yàn)結(jié)果，能夠同時(shí)考慮3個(gè)參數(shù)對(duì)壩料物理力學(xué)特性參數(shù)的影響，使壩料物理力學(xué)特性參數(shù)的表征方法更具綜合性。

圖6 三維威布爾分布擬合

4 結(jié)語(yǔ)

本文通過(guò)系統(tǒng)的概率理論推導(dǎo)研究，并結(jié)合阿爾塔什水利樞紐真實(shí)壩料篩分試驗(yàn)結(jié)果，進(jìn)行了砂礫石壩填料多維概率分布擬合研究與應(yīng)用探討。在構(gòu)建一維威布爾分布模型基礎(chǔ)上，進(jìn)一步推導(dǎo)出了三維威布爾分布函數(shù)表達(dá)式，建立了能夠利用顆粒最大粒徑、礫石含量、曲率系數(shù)間的相關(guān)關(guān)系反映砂礫石填料區(qū)壓實(shí)特性的三維威布爾分布模型，以及利用含泥量、礫石含量、曲率系數(shù)間相關(guān)關(guān)系反映砂礫石填料區(qū)滲透特性的三維威布爾分布模型，使壩料的干密度和滲透系數(shù)的表征方法更具綜合性，可用于砂礫填料的質(zhì)量和變形的評(píng)估和預(yù)測(cè)。