李 林

（中鐵十二局集團鐵路養(yǎng)護工程有限公司，拉薩 250014）

青藏高原海拔平均4 000 m 以上，大部分屬于高海拔區(qū)域，年平均氣溫低至-3 ℃～-7 ℃，導致了多年凍土的形成與保存。青藏鐵路全線穿過長達632 km 的多年凍土區(qū)，其中包括102 km 的融區(qū)［1］。鐵路路基左、右兩側(cè)變形的較大差異會導致路基傾斜、滑坡、坍塌和側(cè)向熱侵蝕等嚴重的工程問題，不僅影響乘客舒適度，嚴重時還會威脅行車安全，增加后期維修成本。因此，研究路基兩側(cè)變形的相關性以及主導變形的方位對鐵路養(yǎng)護部門進行路基維修與治理有極其重要的意義。

目前，國內(nèi)外已有多種指標可用來刻畫兩個隨機變量的相關程度，但都有一定局限性［2-3］。而Copula作為一種靈活且穩(wěn)健的統(tǒng)計工具，得到了廣泛的應用。唐家銀等［4］敘述到可通過Copula 產(chǎn)生多種相關性量化指標進行相關結(jié)構屬性的初判，且具有對分析對象的邊緣分布無要求、數(shù)學處理方便等許多獨特優(yōu)勢。Li 等［5］提出一種Copula 增強卷積技術，解決了交叉變量的相關性；XU 等［6］建立了以混合Copula 函數(shù)進行相關性分析的數(shù)學模型，解決風能與PV 功率輸出最佳容量分配方案。除了相關性上的應用，Chen等［7］以Copula 函數(shù)為基礎，構造了以殘差最小為目標的優(yōu)化模型，得到了退化相關可靠性模型；Zhao 等［8］利用Copula 函數(shù)將邊際故障分布綁定以求得子系統(tǒng)聯(lián)合分布；牛富俊等［9］基于實際觀測數(shù)據(jù)分析了過渡段沉降變形的相關因素并求出了相關系數(shù)；徐東升［10］以有限元軟件進行數(shù)值模擬，對凍土區(qū)路基的融沉進行了可靠性分析；袁莉芬等［11］基于Copula 函數(shù)針對光纖陀螺可靠性評估中的多性能退化量間的相關性，研究了光纖陀螺貯存可靠性評估方法。目前，國內(nèi)外針對多年凍土區(qū)鐵路觀測站點觀測的數(shù)據(jù)對路基兩側(cè)變形相關性的研究較少，基于Copula 函數(shù)對多年凍土區(qū)路基兩側(cè)變形進行相關性研究的文獻就更為鮮見。

本文以相關性數(shù)據(jù)分析的統(tǒng)計工具——Copula相關性理論為基礎，結(jié)合過去10年測得的凍土區(qū)路基特征的真實數(shù)據(jù)，擬合性選擇多年凍土區(qū)路基兩側(cè)變形相關性結(jié)構的最優(yōu)Copula，從而進行其相關性分析，并在線性回歸效應下實現(xiàn)路基兩側(cè)變形主導判定。研究成果可為保護凍土路基以及鐵路養(yǎng)護部門及時采取保護措施、維修與治理提供理論依據(jù)，還可用于路基變形的可靠性分析和預測。

1 路基兩側(cè)變形相關性分析的Copula理論基礎

定義1［12-13］：二維Copula 是一個函數(shù)C：I2=［0，1］2→I=［0，1］且 滿 足C（u，0）C（0，v）=0，C（u，1）=u，C（1,v）=v，當≤u1≤u2≤1、A≤v1≤v2≤1 的存在

定義1 實質(zhì)上是強調(diào)對任意u，v∈［0，1］，二維Copula 函數(shù)C（u，v）均保持非減的性質(zhì)。

定理1（Sklar 定理）：令（X，Y）是待考察隨機變量，其邊緣分布函數(shù)分別是F（x）、G（y），那么存在著唯一的二維Copula CXY（u，v），使得對任意（x，y）∈R2有:

定理1 使用二維Copula 函數(shù)考慮了隨機變量（x，y）的相關性，不僅克服傳統(tǒng)意義上各隨機變量相互獨立的限制條件，而且為二維隨機變量（X，Y）相關分析提供了框架基礎。

二維Copula 函數(shù)還可被看作是服從［0，1］上均勻分布的變量的聯(lián)合分布函數(shù)，當二維Copula 函數(shù)絕對連續(xù)時，存在密度函數(shù)C（u，v），由定理1 知（X，Y）的聯(lián)合密度函數(shù)為：

式中：f（x），g（y）——分別表示兩個隨機變量X，Y的邊緣密度函數(shù)，實質(zhì)上Copula 函數(shù)能夠把隨機變量邊緣分布進行組合研究，這種思路在青藏鐵路多年凍土區(qū)路基變形的探究中起到重要作用。

Copula 函數(shù)由生成函數(shù)決定，在研究中使用最多的 Copula 函數(shù)主要有橢圓 Copula 函數(shù)族和阿基米德Copula 函數(shù)族，幾種經(jīng)典Copula函數(shù)特征如表1所示。

表1 經(jīng)典Copula 函數(shù)特征表

2 路基兩側(cè)變形的相關性主導統(tǒng)計分析

2.1 基于Copula函數(shù)的路基兩側(cè)變形相關性統(tǒng)計分析處理算法

（1）選取各不同路基類型的觀測站點，對數(shù)據(jù)組利用箱線圖原理，即對N 組數(shù)據(jù)樣本進行異常值檢驗（xi，yi）過高或過低，利用異常值檢驗結(jié)果給出不同路基類型觀測站點左、右路基異常值點情況表，并用邊界值代替異常值。

（2）計算各不同觀測站點數(shù)據(jù)的Spearman 相關系數(shù)［15］，根據(jù)正負性初步判斷不同路基類型的路基兩側(cè)高度變形的相關性。

（3）根據(jù)邊緣樣本數(shù)據(jù)（xi，yi）擬合出邊緣分布的經(jīng)驗分布函數(shù)、，分別定義并做出這兩個變量經(jīng)驗分布函數(shù)的散點圖，初步判斷觀測站點左右兩側(cè)的樣本數(shù)據(jù)服從的分布函數(shù)并檢驗。

（4）由（1）中的變量X 和Y 的正（負）相關性，根據(jù)表2選擇正（負）屬性的備選Copula 族cb，b=1，2，…。

表2 相關屬性對應備選Copula 模型表

（5）采用極大似然估計法得到備選Copula 族中cb的參數(shù)向量的估計值［16］。

（6）基于各個Cb（u，v，）隨機生成理論的聯(lián)合估計樣本（ui，vi），i=1，2，…n，并做出散點圖。

（7）基于最小二乘誤差的原理，計算出各個Copula模型下的偏差距離，為：

選擇偏差距離最小的Copula 函數(shù)C0=｛b|min（db（ui，vi），b=1，2，…｝作為變量X 和Y 的相關結(jié)構。

（8）根據(jù)AIC 準則，即在熵的概念上提供一種權衡估計模型復雜度和擬合數(shù)據(jù)優(yōu)良性的標準，備選模型中AIC 最小的為最佳模型，對Copula C0做擬合優(yōu)度檢驗。

2.2 路基兩側(cè)的變形主導統(tǒng)計分析

2.2.1 相關性變形主導變量判斷理論基礎

借助偏度系數(shù)實現(xiàn)主導方向的判定［17］：對于線性回歸方程Y=β0+β1X+ε，通過數(shù)學期望及三階中心矩處理：

式中：ρXY——x與y 之間Pearson 相關系數(shù)，|ρXY|≤1時，則可推斷出≤|γX|。

2.2.2 線性回歸特征下主導變量判定算法

（1）建立兩個隨機樣本X 和Y 的線性回歸方程，存在如下兩種形式：

（2）對建立的線性回歸方程，基于被解釋變量Y使用最小二乘估計出回歸參數(shù)、。

（3）采用t 統(tǒng)計量的方法對方程的回歸系數(shù)進行顯著性檢驗。

（4）由樣本數(shù)據(jù)（xi，yi），i=1，2，…，n，估計出變量X 和Y 的樣本偏度系數(shù)為：

3 路基變形相關性分析及主導方向判定工程實踐

本文的現(xiàn)場數(shù)據(jù)檢測時間從2010年1月-2019年10月，跨度長達118 個月。選取路基類型為碎石護坡、片石護道的K 1 425+350 觀測點，路基類型為一般路基的K 1 447+750 觀測點，路基類型為一般路基、骨架護坡的K 1 449+750 觀測點和路基類型為一般路基、骨架護坡、碎石護坡、護道的K 1 568+800 觀測點，左右兩側(cè)路基高程變化數(shù)據(jù)，部分數(shù)據(jù)如表3所示。

表3 不同測點路基兩側(cè)高度變化（mm）

3.1 路基兩側(cè)變形相關性分析

（1）對各組數(shù)據(jù)進行異常值檢驗，檢驗結(jié)果為僅有K 1 449+750 處的樣本數(shù)據(jù)存在異常值，其余3 組樣本均無異常值。K 1 449+750 處異常值位置如圖1所示。

圖1 K1449+750 異常值檢驗圖

使用邊界值替代異常值后，分別計算4 組數(shù)據(jù)的Spearman 相關系數(shù)指標，結(jié)果如表4所示。

表4 不同路基觀測站點兩側(cè)高度變化Spearman 相關系數(shù)表

由表4可知，不同路基類型兩側(cè)高度的變形呈正相關且相關性極高，Spearman 相關系數(shù)在0.8 以上，極其顯著，說明青藏鐵路運行期間，多年凍土區(qū)路基兩側(cè)高度變形基本滿足一致性。

（2）分別做4 個觀測站點左右兩側(cè)的樣本數(shù)據(jù)概率密度函數(shù)估計圖，初步判斷兩側(cè)路基樣本數(shù)據(jù)服從的分布函數(shù)，結(jié)果如圖2所示：

圖2 4 組樣本數(shù)據(jù)的邊緣概率密度函數(shù)估計圖

由圖2可知，4 組樣本數(shù)據(jù)的邊緣概率密度函數(shù) 呈現(xiàn)集中性、對稱性和均勻變動性，初步判斷凍土各類型路基兩側(cè)高度的變形情況服從正態(tài)分布。利用K-S 檢驗法對其分布進行檢驗，結(jié)果如表5所示。

表5 K-S 正態(tài)性檢驗P 值結(jié)果表

由于K-S 檢驗的P 值均大于0.05，故可認為樣本邊緣分布為正態(tài)分布。

（3）4 組不同路基類型的左右兩側(cè)路基變形趨勢為正相關，由表2可知，即從正屬性的備選Copula 族中選擇左右兩側(cè)路基變形的最優(yōu)相關性結(jié)構。

（4）依次對不同路基類型的數(shù)據(jù)組進行2.1 節(jié)中的（4）～（6）步驟，擬合出最優(yōu)相關性結(jié)構，做4 組樣本的經(jīng)驗分布函數(shù)（，）其中的散點圖和基于最小偏差距離求得的Copula 函數(shù)生成的理論聯(lián)合估計隨機樣本（ui，vi）散點圖，如圖3所示。

圖3 根據(jù)最小偏差距離選出的最優(yōu)Copula 模型

（5）基于AIC 準則對每個正屬性備選Copula 族中的Copula 函數(shù)進行擬合優(yōu)度檢驗，結(jié)果如表6所示。

表6 最優(yōu)相關結(jié)構表

由表6可知，4 組樣本數(shù)據(jù)的相關結(jié)構中t-Copula 對應的AIC 值最小，即路基兩側(cè)變形的最優(yōu)相關結(jié)構為t-Copula 函數(shù)，這說明多年凍土區(qū)路基兩側(cè)變形具有對稱性且上、下尾部相關性較強的特點；不同觀測站點本身的地理因素和外界影響有所差異，t-Copula 的參數(shù)估計不同，導致尾部相關系數(shù)有差異；整體而言，兩側(cè)路基的變化量尾部顯著相關。

3.2 路基兩側(cè)變形主導性統(tǒng)計分析

根據(jù)2.2.2 節(jié)中所述在線性回歸特征下判斷主導變量的算法步驟，首先令多年凍土區(qū)左側(cè)路基的高度變化量為解釋變量，右側(cè)路基的高度變化量為響應變量，得到4 組樣本數(shù)據(jù)都滿足線性回歸方程，其回歸系數(shù)都顯著不為0，且F 統(tǒng)計量顯著?；谏鲜鼋Y(jié)論，做出線性回歸效應下路基兩側(cè)的主導方向相依的判斷，結(jié)果如表7所示。

從表7可以看出，多年凍土區(qū)路基兩側(cè)變形符合線性回歸方差的模型，由兩側(cè)樣本數(shù)據(jù)的偏度系數(shù)可看出，路基兩側(cè)的變形模式幾乎均等性相依，且不同的觀測站點由于內(nèi)部因素和外部影響的差異，變形主導的方向不同。整體而言，全線鐵路的變形主導方位占比接近1∶1。

表7 線性回歸效應下變形主導方位的判斷結(jié)果表

4 路基養(yǎng)護措施建議

基于本文對4 個不同觀測站點的路基左、右兩側(cè)變形及主導方向判斷的工程實踐結(jié)果，給出以下路基養(yǎng)護措施建議：

（1）加強對于多年凍土區(qū)的勘測和研究，全面掌握青藏高原多年凍土區(qū)和季節(jié)性凍土的分布情況，進行凍土現(xiàn)象的區(qū)域劃分工作，掌握可能出現(xiàn)地下冰、凍脹丘、熱融下沉等災害的區(qū)域。

（2）修建凍土區(qū)路基時，應該注意調(diào)整鐵路的弧度和路基的填充方式，避免路基沉降致因側(cè)承受過大的荷載，加速凍土的流變及融沉，導致路基下沉的問題出現(xiàn)。

（3）對于K 1 449+750 處觀測站點，從2010年1月-2019年10月左、右路基高度變化值中融沉凍脹異常點，即左側(cè)路基9 個、右側(cè)路基9 個異常值點處需要引起鐵路養(yǎng)護部門的重點監(jiān)測。

（4）青藏高原多年凍土區(qū)全線鐵路運行時，鐵路養(yǎng)護部門要實時監(jiān)測多年凍土區(qū)路基的高度變化，尤其是主導方位的變形情況，避免因不及時的整治導致另一側(cè)也出現(xiàn)嚴重的地質(zhì)災害。

5 結(jié)論

青藏高原獨特的地理環(huán)境導致多年凍土區(qū)沿線的路基都不同程度地出現(xiàn)了凍脹、融沉等地質(zhì)災害，研究多年凍土區(qū)路基兩側(cè)變形的相關性和對其變形的主導方位進行判斷十分有重要的意義。本文以青藏鐵路多年凍土區(qū)路基變形的數(shù)據(jù)為研究對象，利用箱線圖查找異常值的原理分別對其進行異常值檢驗，再基于Copula 理論給出了路基兩側(cè)變形相關性分析的建模步驟，確定路基兩側(cè)變形的相關性結(jié)構，并依據(jù)不同Copula 族特征對其進行相關性分析，實現(xiàn)多年凍土區(qū)路基兩側(cè)變形的主側(cè)方向判定，得出的主要結(jié)論有：

（1）基于Copula 理論確定凍土區(qū)路基兩側(cè)變形的最優(yōu)相關性結(jié)構為t-Copula，說明了多年凍土區(qū)路基兩側(cè)變形具有對稱性且上、下尾部相關性較強的特點，即若路基一側(cè)出現(xiàn)大幅沉降變形，則很容易導致另一側(cè)路基高程沉降幅度變大。

（2）通過相關性變形的變量主導方向判斷理論判別，實現(xiàn)了路基兩側(cè)主要沉降變形誘導側(cè)的確定，施工中應以此作為主要養(yǎng)護施工對象側(cè)。

（3）路基變形主要致因側(cè)應重點檢測力度，必要時采取加大熱棒、片石護道、片石路肩、土護道等防護措施，避免因不及時的整治而導致另一側(cè)也出現(xiàn)嚴重的地質(zhì)災害。本文相關性統(tǒng)計分析對保護凍土路基、設計路基形式以及鐵路勘察和養(yǎng)護、維修都有重要意義，但針對溫度和路基類型對多年凍土區(qū)路基兩側(cè)變形的相關性結(jié)構和影響的研究還需進一步擴展。