Establishing probabilistic transformation models for geotechnical design parameters using multivariate Gaussian Copula
-
摘要: 岩土工程中经常需要基于直接测量参数去预测设计参数,概率转换模型是确定设计参数真实分布范围的有效工具。常用的基于多维正态分布的概率转换模型构建方法容易引起误差且受边缘分布类型限制。为此,提出了基于多维Gaussian Copula的岩土体设计参数概率转换模型构建方法。以全球黏土数据库CLAY/6/535为例验证了所提方法的有效性,分别推导了不排水抗剪强度和超固结比与孔压静力触探试验指标之间的概率转换模型。结果表明:基于多维Gaussian Copula的概率转换模型构建方法可以将相关结构与边缘分布分开构建,该方法不受岩土体参数的边缘分布类型限制,有效地避免了误差由边缘分布向相关结构的传播。在构建的概率转换模型中,岩土体设计参数的不确定性、相关性与直接测量参数的数目以及直接测量参数与设计参数之间的相关性成反比。
-
关键词:
- 直接测量参数 /
- 设计参数 /
- 概率转换模型 /
- 多维Gaussian Copula /
- 联合分布模型
Abstract: In geotechnical practice, it is common to transform the measured parameters to the design ones. It is also known that the probabilistic transformation models provide an effective tool for predicting the actual range of the design parameters. The commonly-used methods for establishing probabilistic transformation models based on multivariate normal distribution may induce large errors and have a limitation of incorporating various types of marginal distribution for soil parameters. In this study, a new method for establishing the probabilistic transformation models for geotechnical design parameters based on the multivariate Gaussian Copula is proposed. The global clay database CLAY/6/535 compiled in the literature is employed to verify the effectiveness of the proposed method. The probabilistic transformation models from CPTU indices to undrained shear strength and OCR are then derived. The results indicate that by modeling the marginal distribution and dependence structure individually, the proposed method removes the limitation of incorporating various types of marginal distribution and avoids error propagation from marginal distribution to dependence structure. For the proposed probabilistic transformation models, the uncertainty and correlation of design parameters are inversely proportional to the number of the measured parameters and the correlation between the measured and design parameters. -
0. 引言
岩土工程勘察中通常采用多种试验方法测定不同的岩土体参数[1],如天然含水率(含水率试验)、液性指数(界限含水率试验)、标准贯入锤击数(标准贯入试验)、比贯入阻力(单桥静力触探试验)、锥尖阻力与侧壁摩阻力(双桥静力触探试验)等[2]。在众多岩土体参数中,部分参数是岩土工程设计与安全评价的直接输入参数,如水力参数(饱和渗透系数和非饱和土水特征曲线参数)、变形参数(弹性模量和泊松比)、强度参数(黏聚力、内摩擦角和不排水抗剪强度),本文称之为设计参数。一般来说,受技术、经济、场地等条件限制,设计参数的试验数据在岩土工程勘察中不易获得,具有小样本特性[3-5]。研究表明,设计参数的小样本特性会给岩土工程设计与安全评价带来较大的不确定性[6-7],这对确定合理的设计方案、准确评估工程安全性带来了严峻的挑战。尽管如此,对于同一个场地或同一类型岩土体,不同岩土体参数之间存在一定的相关性,因此可以利用其它易获得的岩土体参数预测或估计设计参数,从而以最小的成本获得大量的设计参数试验数据,实现岩土工程的安全设计。为此,本文将这些易获得的用来预测设计参数的岩土体参数称为直接测量参数,如天然含水率、液性指数、标准贯入锤击数、比贯入阻力等。
在岩土工程实践中,为了利用直接测量参数预测设计参数,通常的做法是基于区域或全球数据库建立两者的经验回归方程[8],该经验回归方程也称为设计参数与直接测量参数之间的转换模型。目前,国内外许多学者构建了大量的设计参数与单个或多个直接测量参数之间的转换模型,如Bjerrum[9]基于挪威海洋黏土土样得出了归一化的不排水抗剪强度与塑性指数之间的转换模型。Kulhawy等[10]系统地总结了土体强度参数与常用的室内试验、原位试验参数之间的转换模型。Singh等[11]基于韩国釜山市2个地区的黏土土样试验数据,建立了超固结比与原位直剪试验不排水抗剪强度、孔压静力触探试验锥尖阻力、扁铲侧胀试验指标之间的转换模型。在国内,马海鹏等[12]基于上海地区19个工程201组试验数据,分别推求了黏土、粉质黏土、粉砂3类土中原位直剪试验抗剪强度与单桥静力触探试验比贯入阻力之间的转换模型。屈若枫等[13]基于武汉地区长江一级阶地400余组软土样本,分别研究了淤泥质粉质黏土层黏聚力和内摩擦角与液限、含水率、天然孔隙比之间的转换模型。宋许根等[14]基于珠海西部中心城区大面积深厚软土室内试验1100多组数据、原位试验2600多组数据,研究了含水率与塑性指数、液性指数、压缩模量、无侧限抗压强度、先期固结压力等物理力学参数之间的转换模型。郭林坪等[15]基于天津港地区40个钻孔试验数据,分别建立了天津滨海新区黏性土压缩模量和抗剪强度参数与孔隙比、含水率、密度之间的转换模型。上述研究建立的转换模型为确定性转换模型,即认为设计参数与直接测量参数之间存在确定性函数关系(经验回归方程)。然而受函数形式、试验误差、岩土材料的不均匀性等因素的影响,确定性转换模型计算的设计参数预测值与实测值之间往往存在偏差,且差值随着直接测量参数数值的变化呈现随机波动的特征。设计参数预测值的不确定性已经引起部分学者的注意[16],而考虑这种不确定性影响的有效方法就是建立设计参数与直接测量参数之间的概率转换模型,并推导设计参数均值和标准差与直接测量参数之间的函数关系。
在构建设计参数与直接测量参数之间概率转换模型方面,不少学者也开展了大量的探索研究。Phoon等[17-18]将设计参数的不确定性来源归纳为3类:岩土体参数的固有不确定性、试验误差和转换模型不确定性,在此基础上用随机变量表征上述3类不确定性,通过二阶泰勒展开建立了设计参数不确定性与上述3类不确定性之间的计算公式。基于这种方法,Kulhawy等[10]给出了常用确定性转换模型中设计参数不确定性的表达式,这便构成了设计参数概率转换模型的雏形。Ching等[19-21]在上述研究的基础上对岩土体参数之间相关机制进行了假设,通过已有确定性转换模型推导出设计参数与各个直接测量参数之间的协方差矩阵,采用多维正态分布得出了设计参数概率转换模型。Ching等[22-23]、D’Ignazio等[24]、Liu等[25]、Zou等[26]进一步基于设计参数与直接测量参数的多维试验数据,采用多维正态分布建立设计参数与直接测量参数之间的多维联合分布模型,在此基础上推导出给定直接测量参数条件下设计参数的条件分布及其均值和标准差,将设计参数的条件分布及其均值和标准差作为设计参数的概率转换模型。这种概率转换模型构建方法的优点在于不需要对岩土体参数之间相关机制进行假设,它是目前设计参数概率转换模型构建中较为常用的方法。然而,该方法缺点在于它不能有效地将相关结构的构建与边缘分布的构建分离开来,因而会将边缘分布的识别误差引入相关结构的构造过程中,从而导致设计参数与直接测量参数之间的多维联合分布模型和设计参数概率转换模型的构建存在误差。此外,为了能够得到设计参数条件分布及其均值和标准差的解析表达式,这种方法要求边缘分布和正态分布有直接转换关系,然而岩土体参数的常用边缘分布类型中除了对数正态分布、Johnson分布族等分布与正态分布有直接转换关系外,其余的如极值I型分布、威布尔分布和Beta分布等与正态分布没有直接的转换关系[27-29]。因此,岩土体参数边缘分布的限制使得该方法在岩土工程中的应用受到明显的限制。
Copula理论是构建联合分布模型的重要工具,近年来在金融、环境、水文等领域有着广泛的应用。Copula理论将联合分布分解为边缘分布和相关结构分别进行构造,它不受边缘分布类型的限制,且避免了直接构造联合分布时引入的误差。这些优点使得Copula理论被越来越多地应用于研究二维和多维岩土体参数之间的相关性[30-31]。针对前述基于多维正态分布的设计参数概率转换模型构建方法的缺点,唐小松[32]提出了基于Copula理论的设计参数概率转换模型构建方法。该方法首先利用Copula理论建立设计参数与直接测量参数之间的多维联合分布模型,然后基于Joe[33]提出的条件概率公式,推导了给定直接测量参数条件下设计参数的条件概率密度函数的一般公式。虽然该方法不受边缘分布类型限制,但是该方法得到的表征设计参数概率转换模型的条件概率密度函数没有解析表达式,对于实际工程中关心的设计参数均值和标准差等统计特征只能通过积分求解,难以给出设计参数的均值和标准差为直接测量参数函数的解析表达式,因此无法满足实际工程中根据直接测量参数推求设计参数概率分布的需求。
本文针对上述问题,研究了基于多维Gaussian Copula的一维和多维设计参数概率转换模型构建方法。在此基础上,验证了基于Gaussian Copula的设计参数概率转换模型构建方法的有效性。基于现有黏土实测数据库,分别建立了黏土不排水抗剪强度和超固结比与孔压静力触探试验指标之间的概率转换模型。
1. 多维GaussianCopula
1.1 多维Gaussian Copula的定义
为了便于读者理解本文内容,这里简要介绍多维Gaussian Copula的基本定义。根据Sklar定理[34],在给定n维变量Y1, Y2, …, Yn的边缘分布函数Fi(yi)和联合分布函数F1, 2, …, n(y1, y2, …, yn)的情况下,有且仅有一个Copula函数满足[34]:
F1, 2, ⋯, n(y1, y2, ⋯, yn)=C1, 2,⋯, n(u1, u2, ⋯, un)。 (1) 式中 ui = Fi(yi)为定义在[0, 1]区间内的均匀分布变量;C1, 2, …, n(u1, u2, …, un)为n维Copula函数。联合概率密度函数f1, 2, …, n(y1, y2, …, yn)与Copula的概率密度函数D1, 2, …, n(u1, u2, …, un)有如下关系[34]:
f1, 2, ⋯, n(y1, y2, ⋯, yn)=D1, 2, ⋯, n(u1, u2, ⋯, un)n∏i=1fi(yi), (2) 式中,fi(yi)为Y1, Y2, …, Yn的边缘概率密度函数。由式(1)和(2)可知,多变量联合分布函数与其对应Copula函数相等,而联合概率密度函数可表示为各变量边缘概率密度函数与对应Copula的概率密度函数乘积。Sklar定理将联合分布模型的构造分为两部分:各变量边缘分布的构造和变量间Copula函数的构造。上述两部分的构造相互独立、可分开进行,从而摆脱了传统多维分布模型将联合分布函数直接整体构造的局限,为多变量联合分布模型的构造提供了有效的工具。
对于多维Gaussian Copula而言,令均值为0、协方差矩阵为θ的n维标准正态分布的累积概率分布函数为Φ1, 2, …, n(x1, x2,…, xn; θ),其中协方差矩阵θ对角线上数值为1、且为半正定矩阵,那么n维Gaussian Copula函数的表达式为[35]
C1, 2, ⋯, n(u1, u2, ⋯, un;θ)=Φ1, 2, ⋯, n(Φ−1(u1), Φ−1(u2), ⋯, Φ−1(un); θ), (3) 式中,
ui=Φ(xi) (1 ≤ i ≤ n)。式(3)即是由多维标准正态分布的累积概率分布函数推导的多维Gaussian Copula函数的表达式。由式(3)可知,多维Gaussian Copula的相关参数矩阵θ 等于所对应的多维标准正态分布的协方差矩阵。根据n维Gaussian Copula函数的表达式,可得n维Gaussian Copula的概率密度函数为[35]
D1, 2, ⋯n(u1, u2, ⋯, un; θ) =ϕ1, 2, ⋯, n(x1, x2, ⋯, xn; θ)n∏i=1ϕ(xi) =|θ|−12exp(−12XT(θ−1−I)X), (4) 式中,X = (x1, x2, …, xn)为多维Gaussian Copula的自变量(u1, u2, …, un)通过一维标准正态分布函数的逆函数转换得到的多维标准正态分布变量,其中xi=Φ-1(ui)(i = 1, …, n);I为单位矩阵。
1.2 多维Gaussian Copula相关参数的估计
多维Gaussian Copula相关参数估计有多种方法,其中Pearson线性相关系数方法与Spearman秩相关系数方法是常用的方法,这里主要介绍上述两种方法。
岩土体参数Yi,Yj的Pearson相关系数定义为
ρij=Cov(YiYj)σiσj, (5) 式中,Cov(Yi, Yj)为岩土体参数Yi,Yj之间的协方差,
σi ,σj 分别为岩土体参数Yi,Yj的标准差。岩土体参数Yi,Yj的Pearson相关系数
ρij 与相关参数θ 之间的关系则可表示为[36]ρij=∫∞−∞∫∞−∞(yi−μiσi)(yj−μjσj)× D(Fi(yi), Fj(yj); θij)fi(yi)fj(yj)dyidyj, (6) 式中
μi ,μj 分别表示岩土体参数Yi,Yj的均值;D(Fi(yi),Fj(yj))表示岩土体参数Yi,Yj的二维Gaussian Copula概率密度函数;fi(yi),fj(yj)分别表示岩土体参数Yi,Yj的边缘概率密度函数。在岩土体参数实测数据已知的情况下,式(5)可改写为[36]
ρij=N∑k=1(yik−¯yi)(yjk−¯yj)√N∑k=1(yik−¯yi)2√N∑k=1(yjk−¯yj)2, (7) 式中,{yik, yjk, k = 1, 2, …, N}表示岩土体参数Yi,Yj的一组样本容量为N的实测数据。因此,已知岩土体参数的一组实测数据,根据式(6),(7)便可计算出多维岩土体参数两两间的所有相关参数,如此便估计出了多维Gaussian Copula函数的相关参数矩阵θ。
Spearman秩相关系数是利用两个变量秩的一致性来定义变量之间的相关性,其定义为[36]
rij=3Pr[(Yi−Y′i)(Yj−Y″j)>0] −Pr[(Yi−Y′i)(Yj−Y″j )<0],. (8) 式中,(Yi, Yj)、(Yi′, Yj′)、(Yi′′, Yj′′)为岩土体参数Yi,Yj的3个独立同分布向量。Spearman秩相关系数代表的是岩土体参数的一个向量(Yi, Yj),在Yi,Yj两个维度上分别与任意一个独立同分布向量在这个方向上的一致性概率的3倍与不一致概率之差。Spearman秩相关系数
rij 与Copula函数相关参数θij之间转换关系为[36]rij=12∫10∫10C(ui, uj; θij)duiduj−3。 (9) 当Copula函数类型为Gaussian Copula时,式(9)进一步简化为[36]
θij=2sin(π6×rij)。 (10) 同理,在岩土体参数实测数据已知的情况下,式(8)可改写为[37]
rij=12N∑k=1(rank(yik)rank(yjk))N(N−1)(N+1)−3N+1N−1, (11) 式中,rank(yik)为岩土体参数实测值yik在所有样本{yik, k = 1, 2, …, N}中的秩次[38]。
2. 设计参数概率转换模型构建方法
若n维岩土体参数向量Y = (Y1, Y2, …, Yn)包含设计参数向量Yu = (Yu1, Yu2, …, Yul)与直接测量参数向量Yk = (Yk1, Yk2, …, Ykm),在已知直接测量参数的实测值yk = (yk1, yk2, …, ykm)时,设计参数满足Yui ≤ yui(i = 1, 2, …, l)的概率可表示为
F(yu|yk)=Pr(Yu1≤yu1, ⋯, Yul≤yulYk1=yk1, ⋯, Ykm=ykm)。 (12) 鉴于Φ-1(Fi(∙))为单调递增函数,当式(12)右侧所有的岩土体参数及其观测值进行单调递增变换时,其概率应当保持不变,即
F(yu|yk)=Pr(Φ−1(Fu1(Yu1))≤Φ−1(Fu1(yu1)), ⋯, Φ−1(Ful(Yul))≤Φ−1(Ful(yul))Φ−1(Fk1(Yk1))=Φ−1(Fk1(yk1)), ⋯, Φ−1(Fkm(Ykm))=Φ−1(Fkm(ykm)))。 (13) 令xi = Φ-1(Fi(yi)),即向量X = (X1, X2, …, Xn)与岩土体参数向量Y = (Y1, Y2, …, Yn)满足等概率变换关系Φ(xi) = Fi(yi),则式(13)可写为
F(yu|yk)=Pr(Xu1≤xu1, ⋯, Xul≤xulXk1≤xk1, ⋯, Xkm≤xkm)=FN(xu|xk), (14) 式中 Xu = (Xu1, Xu2, …, Xul)由设计参数向量等概率变换而来,本文称之为设计变量;Xk = (Xk1, Xk2, …, Xkm)由直接测量参数向量等概率变换而来,本文称之为直接测量变量。
由上述公式可知,已知直接测量参数时设计参数的累积概率在进行等概率变换时不会发生改变。特别地,当向量Y所服从的Copula函数类型为Gaussian Copula时,向量X服从均值为0、相关系数矩阵为Gaussian Copula函数相关参数矩阵θ的多维标准正态分布,即X~N(0, θ)。因此,设计参数的条件概率,等于正态分布空间内设计变量在直接测量变量已知条件下的条件概率。基于正态空间内的设计参数条件分布在正态空间内产生模拟样本,并将它由正态空间等概率变换至原始空间,当样本量足够大时样本的累积概率分布将逼近于设计参数的累积条件概率分布,且样本的统计量为设计参数统计量的有效估计,这就是本文所提方法的核心思路。
对于多维正态分布向量X = (X1, X2, …, Xn)而言,在直接测量变量Xk = (Xk1, Xk2, …, Xkm)已知条件下,设计变量Xu = (Xu1, Xu2, …, Xul)的条件分布服从l维正态分布[39]:
Xu|Xk~Nl{(μuT+ΣukΣ−1k(XkT−μkT))T,Σu−ΣukΣ−1kΣukT}。 (15) 式中 上标T代表矩阵的转置;μu和Σu分别代表设计变量Xu的均值向量和协方差矩阵;
μk 和Σk分别代表直接测量变量Xk的均值向量和协方差矩阵;Σuk代表设计变量Xu与直接测量变量Xk间的协方差矩阵。由上文可知,X服从多维标准正态分布,故Xu ~ N(0, θu)、Xk ~ N(0, θk)、Σuk = θuk,因此式(15)可以写为Xu|Xk~ Nl((θukθ−1kXTk)T,θu−θukθ−1kθTuk)。 (16) 将θuk在整体相关参数矩阵θ中的位置标明如下:
θ=[θuθukθkuθk]=[θuθukθTukθk]。 (17) 因此,设计变量Xu在直接测量变量Xk已知的情况下条件分布的均值向量和协方差矩阵为
μ′u=(θukθ−1kXkT)T ,Σ′u=θu−θukθ−1kθTuk , (18) 式中
μ′u 和Σ′u 分别代表在Xk已知条件下设计变量Xu的均值向量和协方差矩阵,从而便于与其原始分布的均值向量μu 和协方差矩阵Σu相区分。在已知l维正态分布的均值向量和协方差矩阵的条件下,相应的l维正态分布函数就唯一确定了。直接测量变量Xk可由直接测量参数Yk通过等概率变换xki = Φ-1(Fki(yki))得到,进而根据已知的直接测量变量和已构建的联合概率分布即可由式(16)确定设计变量多维正态空间内的条件分布函数。基于多维正态空间内的条件分布函数产生模拟样本,再由等概率变换转换到原始空间内,据此可确定设计参数在原始空间内的条件分布及其统计量。多维正态分布除了采用均值向量μ′u 和协方差矩阵Σ′u 表示外,还可以采用均值向量μ′u 、标准差向量σ′u 和相关系数矩阵θ′u 表示,其中标准差向量元素σ′ui=√Σ′uii ,Σ′uii 为协方差矩阵Σ′u 第i行第i列的元素,相关系数矩阵元素θ′uij 由下式给出θ′uij=Σ′uij√Σ′uiiΣ′ujj (1≤i, j≤l)。 (19) 综上可得设计参数概率转换模型构建步骤如下:
(1)构造多维岩土体参数的联合概率分布,确定参数的边缘分布函数Fi(xi)及多维Gaussian Copula的相关参数矩阵
θ 。(2)将直接测量参数向量Yk = (Yk1, Yk2, …, Ykm)通过等概率变换xki = Φ-1(Fki(yki))转换为直接测量变量Xk = (Xk1, Xk2, …, Xkm)。
(3)根据式(16)确定正态空间内设计参数的条件分布,并产生模拟样本Ẋu = (Ẋu1, Ẋu2, …, Ẋul)。
(4)将模拟样本通过等概率变换ẏui =Fui-1(Φ(ẋui))(i = 1, 2, …, l)转换至原始分布空间,可得原始空间内设计参数的模拟样本Ẏu = (Ẏu1, Ẏu2, …, Ẏul),从而根据设计参数的模拟样本估计其条件分布及统计量。
3. 算例
3.1 全球黏土数据库CLAY/6/535
本文以文献Ching等[23]搜集的全球黏土数据库CLAY/6/535为例验证所提方法的有效性。CLAY/6/535数据库包含归一化不排水抗剪强度
su/σ′v 、超固结比OCR、归一化锥尖阻力(qt−σv)/σ′v 、归一化有效锥尖阻力(qt−u2)/σ′v 、归一化超孔隙水压力(u2−u0)/σ′v 和孔压比Bq共6个黏土参数,其中前2个参数为黏土的应力和强度参数,后4个参数则是孔压静力触探试验的试验指标。CLAY/6/535数据库涵盖了中国香港、新加坡、马来西亚、美国、加拿大、委内瑞拉、巴西、瑞士、英国、意大利、挪威等25个地区40个场地共535组实测数据。文献[23]采用Johnson分布族对CLAY/6/535数据库中黏土参数的边缘分布进行拟合,相应的边缘分布及统计特征如表1所示。表中SU和SB都是Johnson分布族成员,它们包括4个参数ax,ay,bx和by,相应的定义域分别为[-∞, ∞]和[by, ay + by]。表 1 CLAY/6/535数据库土体参数边缘分布及统计特征Table 1. Statistics and marginal distributions for five soil parameters in CLAY/6/535 database土体参数 均值 变异系数 边缘分布类型 边缘分布参数 ax ay bx by 归一化不排水抗剪强度Y1 = su/σ′v 0.641 0.596 Johnson SU 1.222 -1.742 0.141 0.250 超固结比Y2 = OCR 2.353 0.657 Johnson SB 0.709 1.887 12.724 0.954 归一化锥尖阻力Y3 = (qt−σv)/σ′v 9.350 0.678 Johnson SU 1.033 -1.438 1.723 4.157 归一化有效锥尖阻力Y4 = (qt−u2)/σ′v 5.280 0.885 Johnson SU 0.989 -1.593 0.868 1.638 归一化超孔隙水压力Y5 = (u2−u0)/σ′v 4.709 0.574 Johnson SU 0.971 -0.762 1.116 3.123 孔压比Y6 = Bq 0.556 0.338 Johnson SU 2.961 0.049 0.544 0.570 3.2 岩土体参数联合分布模型的构建
由Sklar定理[34]可知,岩土体参数多维联合分布模型的构建分为2步:第1步,构建各岩土体参数的边缘分布;第2步,构建岩土体参数的Gaussian Copula函数。因此,本节将首先构建CLAY/6/535数据库六维岩土体参数的边缘分布。由于文献[23]在收集数据库的同时还构建了岩土体参数的边缘分布,这里沿用文献[23]的边缘分布模型,如表1所示。可以看出,除孔压比Bq的变异系数较小以外,其余岩土体参数的变异系数均在0.5以上,这说明全球黏土数据库的数据来源较广、变化范围较大的特点。
为研究岩土体参数间相关关系,首先将CLAY/6/ 535数据库中六维岩土体参数的实测数据样本点{(y1i, y2i, y3i, y4i, y5i, y6i), i = 1, 2, …, 535}转换为均匀分布空间内的样本点{(u1i, u2i, u3i, u4i, u5i, u6i), i = 1, 2, …, 535},进而可得各参数在均匀分布空间内两两间的实测数据散点图,如图1所示。可以看出,CLAY/6/535数据库中岩土体参数间存在不同程度的相关性,除{Y1,Y6},{Y2,Y6},{Y5,Y6}这3组岩土体参数间相关性略低以外,其余岩土体参数间相关性均较强,这也充分说明研究岩土体参数联合分布模型的必要性。
采用式(7)和(11)可计算出CLAY/6/535数据库中各岩土体参数间Pearson相关系数与Spearman相关系数,进一步通过式(6)和(10)可计算出Gaussian Copula的相关参数矩阵,如表2所示。可以看出,{Y1,Y6},{Y2,Y6},{Y5,Y6}这3组岩土体参数间的相关系数较低,而其余岩土体参数间相关系数绝对值均在0.4以上,这些结果与图1结果相互印证,进一步表明数据库中岩土体参数间具有较强的相关性。表2还给出了两种相关系数方法所构造Gaussian Copula的AIC指标和BIC指标,可以看出Spearman秩相关系数方法的AIC指标与BIC指标明显小于Pearson相关系数方法,表明秩相关系数方法所构造的相关结构的拟合效果优于Pearson相关系数方法。
表 2 Pearson相关系数方法与Spearman相关系数方法的比较Table 2. Comparison between Pearson correlation coefficient method and Spearman rank correlation coefficient methodPearson相关系数方法 Spearman秩相关系数方法 相关系数矩阵 [10.6310.6390.5550.527−0.2760.63110.5920.5030.562−0.1500.6390.59210.9120.739−0.4210.5550.5030.91210.439−0.6240.5270.5620.7390.43910.165−0.276−0.150−0.421−0.6240.1651] [10.6520.6790.6420.466−0.3040.65210.5820.5090.538−0.1150.6790.58210.8480.682−0.4760.6420.5090.84810.346−0.7450.4660.5380.6820.34610.221−0.304−0.115−0.476−0.7450.2211] Gaussian Copula相关参数矩阵 [10.5900.6180.5900.514−0.2870.59010.5600.5270.521−0.1490.6180.56010.8450.673−0.4500.5900.5270.84510.474−0.7300.5140.5210.6730.47410.167−0.287−0.149−0.450−0.7300.1671] [10.6690.7000.6600.483−0.3170.66910.6000.5260.556−0.1200.7000.60010.8600.699−0.4930.6600.5260.86010.360−0.7600.4830.5560.6990.36010.231−0.317−0.120−0.493−0.7600.2311] AIC -2894.1 -3171.4 BIC -2829.8 -3107.2 3.3 岩土体设计参数概率转换模型构建方法的有效性验证
基于上述分析,这里选用Spearman秩相关系数方法构造CLAY/6/535数据库联合分布模型和设计参数概率转换模型。考虑到不排水抗剪强度和超固结比分别是边坡稳定性分析和工程中考虑应力历史进行沉降计算的重要参数,分别将归一化不排水抗剪强度
su/σ′v 和超固结比OCR作为设计参数,而将其余4个孔压静力触探试验指标作为直接测量参数。在给出设计参数概率转换模型前,首先采用Leave-one-out交叉验证方法验证设计参数概率转换模型构建方法的有效性,结果如图2所示。Leave-one-out交叉验证方法[40]是统计学中经典的穷举交叉验证方法,其验证过程为:在样本集中抽取一个样本作为测试集,其余的样本作为训练集,计算建立在训练集上的模型在测试集上的精度,重复上述过程,直至遍历所有样本为止。该方法只抽取一个样本作为测试集,不会由于测试集的高度重合而产生过拟合现象,在样本较多时测试精度较高。如前所述,本文算例所采用的全球黏土数据库CLAY/6/535由全球25个地区40个场地的六维岩土体参数实测数据所组成。为此,在Leave-one-out交叉验证过程中,每次选取一个场地的岩土体参数试验数据作为测试集,剩余的39个场地的岩土体参数试验数据作为训练集建立岩土体参数的联合分布模型和设计参数概率转换模型,然后利用测试集中的试验数据验证训练集所建设计参数概率转换模型的有效性。重复上述过程40次,即可遍历所有场地试验数据。可以看出,代表设计参数实测值的红色圆点大部分落入概率转换模型95%置信区间内,且部分实测值与概率转换模型的中位值较为接近,显示了良好的拟合效果。图2(a)中设计参数概率转换模型的95%置信区间下界部分值小于零,这是由于设计参数归一化不排水抗剪强度在归一化之后数值较小,Johnson分布族将其识别为定义域为[-∞, ∞]的Johnson SU分布,故而在条件概率较小时出现负值,表现在图2(a)中则为代表95%置信区间下界的黑色实线出现间断。在图2(b)中,95%置信区间下界则规律性地在1附近波动,这是由于CLAY/6/535数据库中的黏土绝大多数为正常固结土和超固结土,OCR通常大于等于1,因此Johnson分布族将其识别为有下界的Johnson SB分布,表现在图2(b)中则为95%置信区间的下界总是接近或大于1。总体来说,所提方法构建的概率转换模型与设计参数实测值较为接近,从而验证了设计参数概率转换模型构建方法的有效性。
3.4 基于孔压静力触探试验指标的不排水抗剪强度和超固结比概率转换模型
除了CLAY/6/535数据库揭示的孔压静力触探试验指标与土的不排水抗剪强度和超固结比具有较强的相关性外,中国现行规范也为基于孔压静力触探试验指标预测土的不排水抗剪强度和超固结比提供了参考依据。如《岩土工程勘察规范》[1]在条文说明中明确指出:“利用静探资料可估算土的强度参数、浅基或桩基的承载力、砂土或粉土的液化。只要经验关系经过检验已证实是可靠的,利用静探资料可以提供有关设计参数”、“利用孔压静探资料有可能评定土的应力历史,这方面还有待于积累经验”。为此,本文采用所提设计参数概率转换模型构建方法建立了已知不同孔压静力触探试验指标情况下土的不排水抗剪强度和超固结比的概率转换模型,结果如表3所示。由式(16)可知,正态空间内不排水抗剪强度和超固结比在给定孔压静力触探试验指标条件下服从二维正态分布,因此表3中给出了式(18)计算的不排水抗剪强度和超固结比的均值、方差以及两者的相关系数。可以看出,不排水抗剪强度和超固结比的均值都是给定孔压静力触探试验指标的函数,而方差和相关系数则是常数(也受给定孔压静力触探试验指标的影响)。因此,给定孔压静力触探试验指标的具体数值,利用表3中的概率转换模型,即可同时得到设计参数的条件边缘分布和条件联合分布。此外,设计参数的不确定性(方差)和相关性(相关系数)随工程勘察信息(直接测量参数的数目)的增加而减小,且减小的幅度与直接测量参数和设计参数之间的相关性成正比。在直接测量参数数目不变的情况下,直接测量参数与设计参数之间相关性越强,设计参数的不确定性和相关性就越小,这充分说明了引入与设计参数相关性较高的直接测量参数对减小设计参数不确定性、提高概率转换模型精度的重要性。
表 3 基于孔压静力触探试验指标的不排水抗剪强度和超固结比概率转换模型Table 3. Probabilistic transformation models for undrained shear strength and overconsolidation ratio using CPTU test data直接测量参数 设计参数 su/σ′v 设计参数OCR su/σ′v 与OCR相关系数θ′均值 μ′1 方差 σ′21 均值 μ′2 方差 σ′22 X3 0.697X3 0.515 0.600X3 0.640 0.438 X4 0.660X4 0.564 0.526X4 0.723 0.504 X5 0.483X5 0.767 0.556X5 0.691 0.551 X6 -0.317X6 0.900 -0.120X6 0.986 0.671 X3, X4 0.495X3 + 0.235X4 0.500 0.567X3 + 0.0392X4 0.639 0.440 X3, X6 0.714X3 + 0.035X6 0.514 0.715X3 + 0.233X6 0.599 0.442 X4, X6 0.993X4 + 0.438X6 0.483 1.03X4 + 0.664X6 0.537 0.391 X5, X6 0.587X5-0.452X6 0.573 0.616X5-0.262X6 0.626 0.482 X3, X4, X5 0.223X3 + 0.403X4 + 0.182X5 0.491 0.505X3-0.187X4 + 0.505X5 0.566 0.422 X3, X5, X6 1.07X3-0.332X5 + 0.288X6 0.506 0.682X3 + 0.0314X5 + 0.209X6 0.598 0.447 X3, X4, X5, X6 0.764X3 + 0.717X4-0.474X5 + 0.715X6 0.450 0.336X3 + 0.809X4-0.129X5 + 0.691X6 0.527 0.375 4. 结论
本文提出了基于多维Gaussian Copula的岩土体设计参数概率转换模型构建方法。基于全球黏土数据库CLAY/6/535,应用所提方法建立了可广泛适用于正常固结至超固结状态下粉质黏土或黏土的强度应力参数概率转换模型。主要结论如下:
(1)与传统的基于多维正态分布的概率转换模型构建方法相比,基于多维Gaussian Copula的概率转换模型构建方法可以将相关结构与边缘分布分开构建,该方法不受岩土体参数边缘分布类型限制,有效地避免了误差由边缘分布向相关结构的传播,从而为岩土体设计参数概率转换模型的构建提供了有效途径。
(2)常用的Pearson线性相关系数方法在岩土体参数多维联合分布模型中适用性有限,以Spearman相关系数方法为代表的秩相关系数方法简单实用,避免了复杂边缘分布条件下的积分运算,同时所构造的联合分布模型拟合效果好,建议优先采用。
(3)采用所提方法构建的土体强度应力参数概率转换模型中,岩土体设计参数的不确定性、互相关性与已知直接测量参数的数目成反比,即直接测量参数越多,设计参数的不确定性越小、互相关性越弱。此外,设计参数不确定性、互相关性也与直接测量参数、设计参数间的相关性成反比,即直接测量参数与设计参数间的相关性越强,设计参数的不确定性越小、互相关性越弱。
-
表 1 CLAY/6/535数据库土体参数边缘分布及统计特征
Table 1 Statistics and marginal distributions for five soil parameters in CLAY/6/535 database
土体参数 均值 变异系数 边缘分布类型 边缘分布参数 ax ay bx by 归一化不排水抗剪强度Y1 = 0.641 0.596 Johnson SU 1.222 -1.742 0.141 0.250 超固结比Y2 = OCR 2.353 0.657 Johnson SB 0.709 1.887 12.724 0.954 归一化锥尖阻力Y3 = 9.350 0.678 Johnson SU 1.033 -1.438 1.723 4.157 归一化有效锥尖阻力Y4 = 5.280 0.885 Johnson SU 0.989 -1.593 0.868 1.638 归一化超孔隙水压力Y5 = 4.709 0.574 Johnson SU 0.971 -0.762 1.116 3.123 孔压比Y6 = Bq 0.556 0.338 Johnson SU 2.961 0.049 0.544 0.570 表 2 Pearson相关系数方法与Spearman相关系数方法的比较
Table 2 Comparison between Pearson correlation coefficient method and Spearman rank correlation coefficient method
Pearson相关系数方法 Spearman秩相关系数方法 相关系数矩阵 Gaussian Copula相关参数矩阵 AIC -2894.1 -3171.4 BIC -2829.8 -3107.2 表 3 基于孔压静力触探试验指标的不排水抗剪强度和超固结比概率转换模型
Table 3 Probabilistic transformation models for undrained shear strength and overconsolidation ratio using CPTU test data
直接测量参数 设计参数 设计参数OCR 与OCR相关系数θ′ 均值 方差 均值 方差 X3 0.697X3 0.515 0.600X3 0.640 0.438 X4 0.660X4 0.564 0.526X4 0.723 0.504 X5 0.483X5 0.767 0.556X5 0.691 0.551 X6 -0.317X6 0.900 -0.120X6 0.986 0.671 X3, X4 0.495X3 + 0.235X4 0.500 0.567X3 + 0.0392X4 0.639 0.440 X3, X6 0.714X3 + 0.035X6 0.514 0.715X3 + 0.233X6 0.599 0.442 X4, X6 0.993X4 + 0.438X6 0.483 1.03X4 + 0.664X6 0.537 0.391 X5, X6 0.587X5-0.452X6 0.573 0.616X5-0.262X6 0.626 0.482 X3, X4, X5 0.223X3 + 0.403X4 + 0.182X5 0.491 0.505X3-0.187X4 + 0.505X5 0.566 0.422 X3, X5, X6 1.07X3-0.332X5 + 0.288X6 0.506 0.682X3 + 0.0314X5 + 0.209X6 0.598 0.447 X3, X4, X5, X6 0.764X3 + 0.717X4-0.474X5 + 0.715X6 0.450 0.336X3 + 0.809X4-0.129X5 + 0.691X6 0.527 0.375 -
[1] 岩土工程勘察规范:GB 50021—2001[S]. 2009. National Code for Investigation of Geotechnical Engineering: GB 50021—2001[S]. 2009. (in Chinese)
[2] 土工试验方法标准:GB/T 50123—2019[S]. 2019. Standard for Geotechnical Testing Method: GB/T 50123—2019[S]. 2019. (in Chinese)
[3] 张博庭. 用有限比较法进行拟合优度检验[J]. 岩土工程学报, 1991, 13(6): 84-91. doi: 10.3321/j.issn:1000-4548.1991.06.011 ZHANG Bo-ting. Test of goodness of fit using finite comparison method[J]. Chinese Journal of Geotechnical Engineering, 1991, 13(6): 84-91. (in Chinese) doi: 10.3321/j.issn:1000-4548.1991.06.011
[4] 张广文, 刘令瑶. 确定随机变量概率分布参数的推广Bayes法[J]. 岩土工程学报, 1995, 17(3): 91-94. doi: 10.3321/j.issn:1000-4548.1995.03.015 ZHANG Guang-wen, LIU Lin-yao. Extend Bayes method of probability distribution parameters in conforming random variable[J]. Chinese Journal of Geotechnical Engineering, 1995, 17(3): 91-94. (in Chinese) doi: 10.3321/j.issn:1000-4548.1995.03.015
[5] 唐小松, 李典庆, 周创兵, 等. 不完备概率信息条件下边坡可靠度分析方法[J]. 岩土工程学报, 2013, 35(6): 1027-1034. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201306007.htm TANG Xiao-song, LI Dian-qing, ZHOU Chuang-bing, et al. Reliability analysis of slopes with incomplete probability information[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(6): 1027-1034. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201306007.htm
[6] TANG X S, LI D Q, RONG G, et al. Impact of copula selection on geotechnical reliability under incomplete probability information[J]. Computers and Geotechnics, 2013, 49: 264-278. doi: 10.1016/j.compgeo.2012.12.002
[7] LI D Q, TANG X S, ZHOU C B, et al. Characterization of uncertainty in probabilistic model using bootstrap method and its application to reliability of piles[J]. Applied Mathematical Modelling, 2015, 39(17): 5310-5326. doi: 10.1016/j.apm.2015.03.027
[8] General Principles on Reliability for Structures: ISO2394: 2015[S]. 2015.
[9] BJERRUM L. Geotechnical properties of Norwegian marine clays[J]. Géotechnique, 1954, 4(2): 49-69. doi: 10.1680/geot.1954.4.2.49
[10] KULHAWY F H, MAYNE P W. Manual on Estimating Soil Properties for Foundation Design[R]. Ithaca: Cornell University, 1990.
[11] SINGH V K, CHUNG S G. Evaluation of overconsolidation ratios from laboratory and in situ tests on Busan clay[J]. Engineering Geology, 2015, 199: 38-47. doi: 10.1016/j.enggeo.2015.10.006
[12] 马海鹏, 陈祖煜, 于沭. 上海地区土体抗剪强度与静力触探比贯入阻力相关关系研究[J]. 岩土力学, 2014, 35(2): 536-542. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201402034.htm MA Hai-peng, CHEN Zu-yu, YU Shu. Correlations of soil shear strength with specific penetration resistance of CPT in Shanghai area[J]. Rock and Soil Mechanics, 2014, 35(2): 536-542. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201402034.htm
[13] 屈若枫, 徐光黎, 王金峰, 等. 武汉地区典型软土物理力学指标间的相关性研究[J]. 岩土工程学报, 2014, 36(增刊2): 113-119. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2014S2021.htm QU Ruo-feng, XU Guang-li, WANG Jin-feng, et al. Correlations of physical and mechanical properties of typical soft soils in Wuhan[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(S2): 113-119. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2014S2021.htm
[14] 宋许根, 王志勇, 柏威伟, 等. 珠海西部中心城区大面积深厚软土工程特性研究[J]. 岩石力学与工程学报, 2019, 38(7): 1434-1451. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201907013.htm SONG Xu-gen, WANG Zhi-yong, BAI Wei-wei, et al. Study on engineering characteristics of large-scale deep soft soil in the central area of western Zhuhai[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(7): 1434-1451. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201907013.htm
[15] 郭林坪, 穆坤, 杨爱武, 等. 天津滨海新区黏性土物理力学性质指标相关性[J]. 地下空间与工程学报, 2019, 15(2): 452-457. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201902020.htm GUO Lin-ping, MU Kun, YANG Ai-wu, et al. Correlation analysis on physic-mechanical properties of clayey soil in Binhai new area in Tianjin[J]. Chinese Journal of Underground Space and Engineering, 2019, 15(2): 452-457. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201902020.htm
[16] MAYNE P W, PEUCHEN J, BALTOUKAS D. Piezocone evaluation of undrained strength in soft to firm offshore clays[C]//Frontiers in Offshore Geotechnics III: Proceedings of the 3rd International Symposium on Frontiers in Offshore Geotechnics, 2015, Oslo.
[17] PHOON K K. Reliability-based design of foundations for transmission line structures[D]. Ithaca: Cornell University, 1995.
[18] PHOON K K, KULHAWY F H. Evaluation of geotechnical property variability[J]. Canadian Geotechnical Journal, 1999, 36(4): 625-639. doi: 10.1139/t99-039
[19] CHING J, PHOON K K, CHEN Y C. Reducing shear strength uncertainties in clays by multivariate correlations[J]. Canadian Geotechnical Journal, 2010, 47(1): 16-33. doi: 10.1139/T09-074
[20] CHING J, CHEN J R, YEH J Y, et al. Updating uncertainties in friction angles of clean sands[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2012, 138(2): 217-229. doi: 10.1061/(ASCE)GT.1943-5606.0000573
[21] CHING J, PHOON K K. Establishment of generic transformations for geotechnical design parameters[J]. Structural Safety, 2012, 35: 52-62. doi: 10.1016/j.strusafe.2011.12.003
[22] CHING J, PHOON K K. Modeling parameters of structured clays as a multivariate normal distribution[J]. Canadian Geotechnical Journal, 2012, 49(5): 522-545. doi: 10.1139/t2012-015
[23] CHING J, PHOON K K, CHEN C H. Modeling piezocone cone penetration (CPTU) parameters of clays as a multivariate normal distribution[J]. Canadian Geotechnical Journal, 2014, 51(1): 77-91.
[24] D’IGNAZIO M, PHOON K K, TAN S A, et al. Correlations for undrained shear strength of fnnish soft clays[J]. Canadian Geotechnical Journal, 2016, 53(10): 1628-1645.
[25] LIU S, ZOU H, CAI G, et al. Multivariate correlation among resilient modulus and cone penetration test parameters of cohesive subgrade soils[J]. Engineering Geology, 2016, 209: 128-142.
[26] ZOU H, LIU S, CAI G, et al. Multivariate correlation analysis of seismic piezocone penetration (SCPTU) parameters and design properties of Jiangsu quaternary cohesive soils[J]. Engineering Geology, 2017, 228: 11-38.
[27] LUMB P. Safety factors and the probability distribution of soil strength[J]. Canadian Geotechnical Journal, 1970, 7(3): 225-242.
[28] 苏永华, 何满潮, 孙晓明. 大子样岩土随机参数统计方法[J]. 岩土工程学报, 2001, 23(1): 117-119. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200101026.htm SU Yong-hua, HE Man-chao, SUN Xiao-ming. Approach on asymptotic approximations of polynomials for probability density function of geotechnics random parameters[J]. Chinese Journal of Geotechnical Engineering, 2001, 23(1): 117-119. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200101026.htm
[29] 陈立宏, 陈祖煜, 刘金梅. 土体抗剪强度指标的概率分布类型研究[J]. 岩土力学, 2005, 26(1): 37-40. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200501008.htm CHEN Li-hong, CHEN Zu-yu, LIU Jin-mei. Probability distribution of soil strength[J]. Rock and Soil Mechanics, 2005, 26(1): 37-40. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200501008.htm
[30] 张蕾, 唐小松, 李典庆, 等. 基于Copula函数的岩土结构物系统可靠度分析[J]. 岩土力学, 2016, 37(1): 193-202. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201601024.htm ZHANG Lei, TANG Xiao-song, LI Dian-qing, et al. System reliability analysis of geotechnical structures based on the Copula function[J]. Rock and Soil Mechanics, 2016, 37(1): 193-202. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201601024.htm
[31] WANG F, LI H. Distribution modeling for reliability analysis: Impact of multiple dependences and probability model selection[J]. Applied Mathematical Modelling, 2018, 59: 483-499.
[32] 唐小松. 基于Copula理论的岩土体参数不确定性建模与可靠度分析[D]. 武汉: 武汉大学, 2014. TANG Xiao-song. Uncertainty Modeling of Correlated Geotechnical Parameters and Reliability Analysis Using Copulas[D]. Wuhan: Wuhan University, 2014. (in Chinese)
[33] JOE H. Families of m-variate distributions with given margins and m(m-1)/2 bivariate dependence parameters[J]. Lecture Notes-Monograph Series, 1996: 120-141.
[34] SKLAR A. N-dimensional distribution functions and their margins [J]. Publications of the Institute of Statistics of the University of Paris, 1959, 8: 229-231. (in French)
[35] CHERUBINI U, LUCIANO E, VECCHIATO W. Copula Methods in Finance[M]. New York: John Wiley & Sons, 2004.
[36] MARI D D, KOTZ S. Correlation and Dependence[M]. Singapore: World Scientific, 2001.
[37] NELSEN R B. An introduction to Copulas[M]. 2nd ed. New York: Springer, 2006.
[38] IMAN R L, CONOVER W J. A distribution-free approach to inducing rank correlation among input variables[J]. Communications in Statistics-Simulation and Computation, 1982, 11(3): 311-334.
[39] DING P. On the conditional distribution of the multivariate t distribution[J]. The American Statistician, 2016, 70(3): 293-295.
[40] KEARNS M, RON D. Algorithmic stability and sanity-check bounds for leave-one-out cross-validation[J]. Neural Computation, 1999, 11(6): 1427-1453.
-
期刊类型引用(2)
1. 王军,刘志明,蔡国军,叶飞龙,宋小进. 基于砂土界面剪切试验的自传感压电土工电缆监测效果评价. 岩土工程学报. 2023(10): 2023-2031 . 本站查看
2. 崔新壮,包振昊,郝建文,杜业峰,张圣琦,李向阳,赵延涛. 重载公路路基动力响应现场测试与三维空间分布规律. 中国公路学报. 2023(10): 75-83 . 百度学术
其他类型引用(2)