基于多维Gaussian Copula的岩土体设计参数概率转换模型构建方法

李典庆; 吕天健; 唐小松

doi:10.11779/CJGE202109003

基于多维Gaussian Copula的岩土体设计参数概率转换模型构建方法 English Version

李典庆^{1, 2,},
吕天健^{1, 2},
唐小松^{1, 2}

1.
武汉大学水资源与水电工程科学国家重点实验室,湖北　武汉　430072
2.
武汉大学工程风险与防灾研究所,湖北　武汉　430072

基金项目:

国家重点研发计划项目 2017YFC1501301

国家自然科学基金项目 51879204

国家自然科学基金项目 41861144022

详细信息

作者简介:
李典庆(1975— ),男,博士,教授,主要从事水工岩土工程可靠性与风险控制方面的研究。E-mail：dianqing@whu.edu.cn

中图分类号: TU431
计量
- 文章访问数: 311
- HTML全文浏览量: 16
- PDF下载量: 255
出版历程
- 收稿日期: 2020-05-17
- 网络出版日期: 2022-12-02
- 刊出日期: 2021-08-31

Establishing probabilistic transformation models for geotechnical design parameters using multivariate Gaussian Copula

LI Dian-qing^{1, 2,},
LÜ Tian-jian^{1, 2},
TANG Xiao-song^{1, 2}

1.
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
2.
Institute of Engineering Risk and Disaster Prevention, Wuhan University, Wuhan 430072, China

摘要

摘要: 岩土工程中经常需要基于直接测量参数去预测设计参数,概率转换模型是确定设计参数真实分布范围的有效工具。常用的基于多维正态分布的概率转换模型构建方法容易引起误差且受边缘分布类型限制。为此,提出了基于多维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.
- measured parameter /
- design parameter /
- probabilistic transformation model /
- multivariate Gaussian Copula /
- joint distribution model

HTML全文

0. 引言

成都市主城八区均下伏中等风化泥岩,埋深一般为-0.5～-35m,层厚十几米到几十米不等^[1-2],工程性质差,承载力估算困难。随着成渝地区双城经济圈建设的深度推进,成都地区已开建或筹建的高层、超高层建筑的基础持力层均涉及中等风化泥岩,研究该类岩体的承载能力具有一定的工程意义。

江苏^[3]、广东^[4]、重庆^[5]等地在实践中逐渐修正了规范^[6]对软岩地基承载力的取值建议,提出了不同岩石单轴抗压强度分区所对应的折减系数及承载力；长沙、兰州等地以岩基载荷试验为主,辅以点荷载试验、旁压试验等,提出了泥岩承载力的界定标准^[7-9]；并通过不同类型试验结果的对比,探讨了影响地基承载力取值的因素^[10-11],指出确定地基承载力时应考虑岩体的物理力学性质、成因、年代、地下水及其变化情况、建筑物对沉降及差异沉降的要求等^[12]。成都地标DB51/T5026—2001^[13]建议中等风化泥岩地基承载力特征值为500～1000 kPa,也可根据室内天然单轴抗压强度折减法确定地基承载力,工程实践表明规范条文存在两点不足。一方面是建议的承载力取值偏低,成都地区中风化泥岩极限承载力值试验值普遍分布在2000～6000 kPa,且主要集中在4000～6000 kPa^[14],按照安全系数3考虑,其承载力特征值亦大于1000 kPa；另一方面,折减法确定承载力合宜性差,折减系数ψ的取值无既定条款可循,往往带有过多的人为因素,即使对同一类型、同一条件、同一环境的地基岩体,不同的设计人员也会取不同的ψ值,有时甚至出入较大,在没有原位试验的情况下,所提出的承载力缺乏依据,带有一定的随意性,对工程的安全性没有切实保障。在成都地区更为合理的取值研究至今鲜有突破,受到不少诟病。

综上,本文通过对24个场地约100余处原位平板载荷试验与同条件岩石单轴抗压强度试验工程案例进行对比分析,讨论了成都地区中等风化泥岩计算公式中ψ值取值及其与承载力特征值的相关关系。进一步,结合成都高新区在建的某超高层建筑项目,采用原位平板载荷试验（9组）联合钻孔旁压试验法（17孔）针对性的探讨了该类岩体承载力特性,并分析了承压板尺寸等因素对承载能力发挥程度的影响。以上研究旨在探寻单轴抗压强度或旁压试验结果作为承载力设计依据的可能。

1. 成都中等风化泥岩特性

成都地区红层属于盆西小区^[15],是最具代表性的近水平红层分布区。区内中等风化泥岩以J_3p侏罗系蓬莱镇组上统和K_2g白垩系中统灌口组为主（见图1）,其中J_3p泥岩大范围的分布在成都天府新区高楼山—秦皇寺中央商务区—倒石桥一带,其他区域均为K_2g泥岩,从城区东南方至西北方埋深为10～30 m,深度逐渐加深。

图 1 成都地区中等风化泥岩分布略图

Figure 1. Distribution of moderately weathered mudstone in Chengdu

下载: 全尺寸图片幻灯片

区内中等风化泥岩呈紫红色,层理清晰,风化裂隙较发育,巨厚层构造,整体结构。锤击易碎、声哑,用镐难挖掘,岩芯钻方可钻进。岩芯采取率达90%以上（图2）。通过对岩芯的X衍射和矿物成份分析,可知区内泥岩及其风化物含有超50%的SiO₂、Al₂O₃、Fe₂O₃,其不同的组合构成了以蒙脱石、伊利石为主的片状黏土矿物。表1为岩体物理力学参数统计表。

图 2 成都地区典型中等风化泥岩岩芯

Figure 2. Typical moderately weathered mudstone core in Chengdu

下载: 全尺寸图片幻灯片

表 1 中等风化泥岩物理力学性质统计表

Table 1. Physical and mechanical properties of rock

统计指标	天然密度/(g·cm^-3)	单轴抗压强度/MPa		黏聚力/MPa		内摩擦角/(°)
统计指标	天然密度/(g·cm^-3)	天然状态	饱和状态	天然状态	饱和状态	天然状态	饱和状态
最大值	2.62	9.67	8.20	1.2	0.85	40.1	38.6
最小值	2.30	2.71	1.56	0.4	0.23	33.7	31.4
平均值	2.48	6.00	3.88	0.8	0.49	37.9	35.5
样本容量	141	51	39	28	17	28	17

下载: 导出CSV

| 显示表格

2. 原位平板载荷试验与岩石单轴抗压强度试验对比

目前成都地区确定岩石地基承载力特征值主要为岩石单轴抗压强度折减法和原位平板岩基载荷试验法：①第一种方法简便但计算结果与实际地基承载力有较大出入；②第二种计算方法得到的结果普遍被认为是真实的,但费时费力。

笔者整理了24个工程实例中约100个中等风化泥岩的原位平板载荷试验点的相关测试数据,开展了2种方法确定地基承载力的对比分析,结果见表2,试验点覆盖了成都主城八区,具有较好的代表性。需要说明的是,文中所示意的原位平板载荷试验点绝大多数未达到明显的极限荷载,因此岩基的承载力特征值确定标准：①极限荷载f_u/3与比例界限荷载的较小值；②未达到极限状态取比例界限荷载。

表 2 岩基载荷试验与单轴抗压强度试验对比

Table 2. Comparison between load tests and uniaxial compressive strength tests

工程名称	工程地点	f_rk/MPa	f_a/kPa	ψ	停载位移/mm	点数
隆鑫	成华区	4.00	2044	0.51	5.420	3
润都科技	成华区	4.00	1981	0.49	6.760	3
蓉华·上林	成华区	3.05	1841	0.60	3.990	8
			1773	0.58	9.730
			1861	0.61	8.680
科创中心	成华区	5.00	2263	0.45	18.370	8
东方荟	锦江区	5.00	2157	0.43	3.640	4
东大街项目	锦江区	3.07	1900	0.62	3.520	3
仁恒置地广场	锦江区	3.43	1990	0.58	—	2
航天科技大厦	锦江区	6.00	2880	0.48	—	3
西部国金	锦江区	1.48	843	0.57	10.410	3
融御大楼	锦江区	7.30	3010	0.41	0.562	3
海桐三期	锦江区	7.90	2959	0.37	8.300	4
滨江30号地块	高新区	6.80	2740	0.40	10.700	3
半岛城邦	高新区	4.50	2021	0.45	6.600	9
龙湖世纪城	高新区	2.79	1800	0.65	33.190	8
川大科技园	双流区	3.92	1894	0.48	5.290	3
水电·云立方	龙泉驿	5.98	2400	0.40	5.500	3
视高天府公园	天府新区	6.48	2400	0.37	7.320	3
学府城	天府新区	4.99	2290	0.46	1.420	4
学府城	天府新区	4.99	2266	0.45	2.720	2
仁中风华一期	天府新区	5.64	2400	0.43	5.980	6
驿都水岸丽园	天府新区	5.58	2524	0.45	17.500	1
驿都水岸丽园	天府新区	5.00	2162	0.43	18.630	2
东方希望天祥	武侯区	4.39	2000	0.50	1.160	3
海洋乐园	武侯区	6.55	2500	0.38	4.260	6
中邮金融中心	武侯区	5.74	2400	0.42	2.200	3
田园晶品	新都区	3.1	1697	0.55	3.720	6
注：f_rk为岩石天然单轴抗压强度标准值；f_a为承载力特征值；ψ为折减系数。

下载: 导出CSV

| 显示表格

表2中,成都主城八区中等风化泥岩地基承载力特征值试验值70%以上集中在1500～2500 kPa,较规范DB51T5026建议最大值普遍大1.5倍以上,统计原位平板载荷试验与同条件岩石单轴抗压强度试验结果得出折减系数分布在0.37～0.70,总体值大于0.5（表3）；综合可见,折减系数呈现出随单轴抗压强度增加而递减的趋势。

表 3 折减系数和承载力特征值建议值

Table 3. Recommended values of reduction coefficient and characteristic value of bearing capacity

f_rk/MPa	ψ	f_a/kPa	说明
<2	0.7	500~1400	符号意义同表2,岩石单轴抗压强度为天然状态
2~4	0.5~0.7	1400~2000
4~6	0.4~0.5	2000~2400
6~8	0.37~0.4	2400~3000
>8	0.37	>3000

下载: 导出CSV

| 显示表格

表3所列的统计结果对于岩石单轴抗压强度小于2 MPa和大于8 MPa的岩体,由于没有试验样本,可分别取统计数据中折减系数最高区段和最低区段的ψ值,即取0.7,0.37,承载力特征值最小值按现行规范^[12]最小值500 kPa考虑。同时,上述统计结果多数是在未加载到岩基破坏的试验条件下得到的,实际的折减系数比该统计值还要大。但是,仍要注意的是,表3中岩石单轴抗压强度为天然状态,若当基岩浸水软化时,需要适当考虑折减系数减小。

3. 承载力的影响因素分析

进一步依托成都高新区在建的某超高层建筑项目,探讨不同承压板直径对研究区中等风化泥岩地基承载力取值的影响。

3.1 场地条件

该超高层建筑位于成都市南部天府新区秦皇寺中央商务区。场地地貌单元属宽缓浅丘,为剥蚀型浅丘陵地貌,建筑基础底板持力层为侏罗系蓬莱镇组（J_3p）中等风化泥岩。泥岩天然密度2.48 g/cm³、自然状态下单轴抗压强度平均值6.00 MPa、饱和状态单轴抗压强度平均值3.88 MPa。三轴天然抗剪强度参数平均值内摩擦角为37.9°,黏聚力为0.8 MPa。同时,根据钻孔波速测得岩体波速为2200～3200 m/s；室内岩块波速平均值为2370 m/s。

3.2 试验方案

承压板载荷试验在场地开挖的3口深井平硐中进行（编号分别为SJ01,SJ02,SJ03）。试验点点位见图3所示。

图 3 试验点位示意图

Figure 3. Schematic diagram of test sites

下载: 全尺寸图片幻灯片

图3中,SJ01深30 m（高程459.06 m）、SJ02深36 m（高程451.15 m）、SJ03深39 m（高程447.57 m）,井底平硐长×宽×高为2 m×2 m×8 m。每一平硐内分别开展承压板直径R为300,500,800 mm的载荷试验各1组,共9组。载荷试验操作过程严格参考规范^[16],开挖至目标层后对试验点进行粗砂找平,放置承压板、千斤顶、百分表等试验设备。检查无误后采用慢速维持加载法加载,加载分8～12级,每级加载后连续三次读数之差均不大于0.01 mm可施加下一级荷载；当某级荷载下24 h变形速率不能达到相对稳定标准或本级荷载的变形量大于前级荷载变形量的5倍或承压板周边的岩土出现明显隆起或径向裂缝持续发展等现象终止加载。试验现场照片见图4。

图 4 岩基平板载荷试验过程照片

Figure 4. Photos of bearing plate tests

下载: 全尺寸图片幻灯片

试验数据绘制p–s曲线,曲线起始直线段终点的前一级荷载为比例界限荷载；终载前一级荷载为极限荷载。将极限荷载除以3的安全系数,所得值与对应于比例界限的荷载相比较,取小值定为承载力特征值。

3.3 试验结果分析

典型p–s曲线见图5所示。

图 5 p–s曲线（以02深井300 mm承压板为例）

Figure 5. p–s curves (deep well No. 2, plate diameter 300 mm)

下载: 全尺寸图片幻灯片

（1）曲线整体上表现出变形随着荷载的增大而逐渐增大。

（2）以SJ02中直径300 mm载荷试验结果为例,试验荷载分12级加载,试验结束荷载为7800 kPa。加载结束时,位移持续增加且不能稳定,累计变形量为4.182 mm。从曲线上均可直接读出比例界限荷载和极限荷载。曲线第一个拐点（2400 kPa）对应的前一级荷载（2100 kPa）为比例界限荷载；极限荷载为曲线第二个拐点（7800 kPa）对应的前一级荷载（6300 kPa）。极限荷载的1/3与比例界限荷载对比取小值为承载力特征值,即为2100 kPa。

承载板直径与承载力特征值之间的关系见图6所示。其他试验点试验结果统计见表4所示。

图 6 承载板直径与承载力特征值的关系

Figure 6. Plate diameter vs. bearing capacity

下载: 全尺寸图片幻灯片

表 4 研究区中等风化泥岩承载力特征值试验值统计表（承压板载荷试验）

Table 4. In-suit test values of bearing capacity of moderately weathered mudstone in study area

试验点号竖井编号	承压板直径D/mm	原位试验				极限值/3对应荷载f_a/kPa	岩体完整性指数K_v
试验点号竖井编号	承压板直径D/mm	最大加载量/kPa	比例界限对应荷载/kPa	终止条件	裂缝情况	极限值/3对应荷载f_a/kPa	岩体完整性指数K_v
SJ01 （459.06 m）	300	9000	3000	趋近试验装置系统极限	无	2400	0.46
	500	9600	2400	位移持续增大,不能稳定	无	2400	0.46
	800	6300	2700	沉降增大,大于前一级沉降量的5倍	有	2300	0.43
SJ02 （451.15 m）	300	7800	2400	变形不能保持稳定	有	2100	0.39
	500	10800	2400	变形不能保持稳定	无	2000	0.39
	800	11700	2700	沉降增大,大于前一级沉降量的2倍	有	2100	0.39
SJ03 （447.54 m）	300	11800	3600	趋近试验装置系统极限	有	2800	0.49
	500	9600	3000	趋近试验装置系统极限	有	3000	0.50
	800	10800	3600	支硐上壁无法提供足够反力	无	2700	0.43

下载: 导出CSV

| 显示表格

对于同一平硐,中等风化泥岩地基承载力特征值差异小于10%。具体来说,SJ01深井泥岩承载力特征值为2300～2400 kPa；SJ02深井泥岩承载力特征值为2000～2100 kPa；SJ03深井泥岩承载力特征值为2700～3000 kPa。说明,相似性状的泥岩,载荷板尺寸对于其地基承载力性能影响程度不大（见图7所示）,可推测基础宽度对于同一性状的中等风化泥岩地基承载力的影响不明显。

图 7 旁压试验过程照片

Figure 7. Photos of pressuremeter tests

下载: 全尺寸图片幻灯片

不同平硐内试验结果略有差异,与岩体个体差异尤其是岩体完整性和矿物组分不无关系,试验表明：三口平硐内泥岩岩体完整性指数分别为0.43～0.46,0.39,0.43～0.50；黏土矿物含量分别为29%～30%,36%～47%,17%～25%。上述差异导致不同平硐内中等风化泥岩地基承载性能出现差异。另外,岩石地基承载性能差异亦是试验环境条件、操作人为误差等因素综合影响所造成,需进一步详细甄别和探讨。

4. 承载力确定方法讨论

4.1 承载力特征值修正折减法

上述超高层建筑项目中,3口深井平硐中进行承压板载荷试验同时亦在300 mm压板试验点位附近各取3个试样进行了室内单轴抗压强度试验（试验参照规范^[16]）,对修正折减系数法进行验证验证。根据（DB51/T5026—2001）^[13]和（GB50007—2011）^[6]中确定承载力特征值 $f_{a}$ 公式为

$f_{a} = ψ f_{rk}$ ,

(1)

式中, $f_{a}$ 为岩石地基承载力特征值（kPa）, $f_{rk}$ 为岩石饱和单轴抗压强度标准值（kPa）,软岩取为天然单轴抗压强度标准值）, $ψ$ 为折减系数,由当地经验确定。

结合前述研究,对式（1）参数进行修正,得

$f_{a} = ψ^{'} f_{rk}$ ,

(2)

式中, $f_{a}$ 为岩石地基承载力特征值（kPa）, $f_{rk}$ 为岩石饱和单轴抗压强度标准值（kPa,软岩取为天然单轴抗压强度标准值）, $ψ^{'}$ 为修正折减系数（按表3选取）。

试验及验证结果见表5所示。

表 5 修正折减系数法验证结果

Table 5. Test results of modified reduction coefficient method

竖井编号	编号	f_a/kPa	f_rk/MPa	平均f_rk/MPa	ψ	计算值	差异/%
SJ01	1	2400	4.9	5.96	0.40	2384	0.6
	2		5.8
	3		7.19
SJ02	1	2100	3.47	4	0.50	2000	4.7
	2		3.78
	3		4.39
SJ03	1	2800	6.89	7.1	0.39	2769	1.1
	2		8.31
	3		8.41

下载: 导出CSV

| 显示表格

从表5中计算结果可见,本文所提出的折减系数修正法所估算的岩基承载力特征值与原位载荷试验结果相近,差异在5%以内,效果较好,可以作为承载力计算的依据。根据单轴抗压强度分区不同其修正系数可按表3选取,取值区间可采用线性内插来确定。

4.2 旁压试验适宜性探讨

旁压试验其实质是横向载荷试验。软岩中旁压临塑荷载既与岩石单轴抗压强度有关又与岩体的裂隙发育程度有关,利用旁压试验能更好地反映工程岩体的实际强度^[10-11]。近年来,在长沙等地得到广泛应用,然而其在成都地区中等风化泥岩中的适宜性仍需进一步探讨。在上述项目场地中14口勘察钻孔+3口对比钻孔中旁压试验（见图5）,测试深度15～30 m不等。其中,3口对比钻孔TL05,JK14,JK12紧邻SJ01,SJ02,SJ03使试验结果具有可对比性。

采用PM-2B型预钻式旁压仪进行试验（见图7）,主机尺寸23 cm×36 cm×90 cm,最大压力6 MPa,试验操作严格参考规范[17,18],根据试验结果绘制P–V曲线,进一步确定初始压力P₀和临塑压力P_f。两个压力差值为承载力特征值f_ak。

旁压试验结果见表6所示。

表 6 旁压试验结果统计表

Table 6. Results of pressuremeter tests

钻孔编号	测试深度/m	P_f/kPa	P₀/kPa	f_a/kPa	钻孔波速/(m·s^-1)	岩体完整性指数K_v
JK04	28.5	3002.4	625.2	2377.2	2465	0.380
JK03	21.0	1399.5	350.0	1049.5	1873	0.231
JK10	20.0	1294.2	310.5	983.7	2023	0.256
JK12	13.0	1495.5	311.2	1184.3	2160	0.292
JK12	20.0	2670.0	698.5	1971.5	2204	0.304
JK14	28.0	3753.0	850.6	2902.4	2912	0.507
TL01	26.5	3790.2	849.6	2940.6	2445	0.414
TL03	18.5	4085.1	1120.3	2964.8	2754	0.530
TL05	28.0	2828.4	475.6	2352.8	2542	0.447
TL11	20.5	3715.9	850.6	2865.3	3013	0.469
TL13	16.5	3342.3	942.3	2400.0	2555	0.452
TL13	23.0	2209.1	487.2	1721.9	2241	0.348
TL14	17.5	2300.5	422.1	1878.4	2613	0.473
TL16	14.5	2489.0	475.6	2013.4	2178	0.296
TL16	20.5	2737.9	475.6	2262.3	2824	0.498
TL18	23.0	2590.6	432.2	2158.4	2523	0.441
TL22	23.0	2856.7	511.2	2345.5	2930	0.415
TL25	24.5	2707.0	530.1	2176.9	2347	0.381
TL25	26.0	4773.5	681.2	4092.3	3094	0.663
TL35	23.0	3032.5	625.2	2407.3	2659	0.490
TL37	18.0	2494.4	550.0	1944.4	2690	0.501

下载: 导出CSV

| 显示表格

从表6中可见,中等风化泥岩地基承载力特征值测试结果最小值为983.7 kPa,最大值为4092.3 kPa,范围较大,这与每一测试孔测试时周围环境、成孔质量、岩体质量、测试深度等有一定关系。因测试时环境和成孔质量属于随机无量化指标的影响因素。进一步对比表4,6中深井及其对比钻孔试验结果,以探究旁压试验适宜性。以深井平硐内的承压板试验结果为基准。对比结果见图8所示。

图 8 试验结果对比图

Figure 8. Comparison between test results

下载: 全尺寸图片幻灯片

对比结果显示,各测点的泥岩均为中等风化泥岩,测试深度、岩体完整性均相似,两种方法测试结果极为相近。尤其是深井SJ01和钻孔TL05,测试深度分别为30 m和28 m,岩体完整性指数均为0.46左右,得到的地基承载力特征值2300 kPa范围,可以认为采用承压板载荷试验、钻孔旁压试验确定的泥岩地基承载力特征值f_a可互为印证。

5. 结论

本文通过对24个场地约100个原位平板载荷试验与同条件岩石单轴抗压强度试验工程案例进行对比分析,讨论了成都地区中等风化泥岩计算公式中折减系数取值。进一步,依托成都高新区在建的某超高层建筑项目,探讨了修正折减系数法和旁压试验结果作为设计依据的建议。

（1）成都主城八区中等风化泥岩地基承载力特征值主要集中在1500～2500 kPa。

（2）统计原位平板载荷试验与同条件岩石单轴抗压强度试验结果得出折减系数分布在0.37～0.70,总体值大于0.5。

（3）本文所提出的折减系数修正法所估算的岩基承载力特征值在实际工程中实用性较强,其估算结果与原位载荷试验结果差异在5%以内,效果较好,可以作为承载力计算的依据。

（4）对于相似性状的泥岩,泥岩地基承载力受承压板直径影响不显著；且承压板载荷试验、钻孔旁压试验确定的地基承载力特征值可互为印证；

（5）上述成果,可用于该地区同类岩石地基参考,但是岩土体物理力学参数间的相关性取决于样本容量,对上述相关内容的细化有待更多的试验对比和检验及完善。

图 1 实测数据在标准均匀分布空间内的散点图

Figure 1. Measured data points in standard uniform space for soil parameters in CLAY/6/535 database

下载: 全尺寸图片幻灯片

图 2 设计参数概率转换模型构建方法有效性验证

Figure 2. Validation of proposed method for establishing probabilistic transformation models for design parameters

下载: 全尺寸图片幻灯片

表 1 CLAY/6/535数据库土体参数边缘分布及统计特征

Table 1 Statistics and marginal distributions for five soil parameters in CLAY/6/535 database

土体参数	均值	变异系数	边缘分布类型	边缘分布参数
土体参数	均值	变异系数	边缘分布类型	a_x	a_y	b_x	b_y
归一化不排水抗剪强度Y₁ = $s_{u} / {σ^{'}}_{v}$	0.641	0.596	Johnson SU	1.222	-1.742	0.141	0.250
超固结比Y₂ = OCR	2.353	0.657	Johnson SB	0.709	1.887	12.724	0.954
归一化锥尖阻力Y₃ = $(q_{t} - σ_{v}) / {σ^{'}}_{v}$	9.350	0.678	Johnson SU	1.033	-1.438	1.723	4.157
归一化有效锥尖阻力Y₄ = $(q_{t} - u_{2}) / {σ^{'}}_{v}$	5.280	0.885	Johnson SU	0.989	-1.593	0.868	1.638
归一化超孔隙水压力Y₅ = $(u_{2} - u_{0}) / {σ^{'}}_{v}$	4.709	0.574	Johnson SU	0.971	-0.762	1.116	3.123
孔压比Y₆ = B_q	0.556	0.338	Johnson SU	2.961	0.049	0.544	0.570

下载: 导出CSV

表 2 Pearson相关系数方法与Spearman相关系数方法的比较

Table 2 Comparison between Pearson correlation coefficient method and Spearman rank correlation coefficient method

	Pearson相关系数方法	Spearman秩相关系数方法
相关系数矩阵	$[\begin{matrix} \begin{matrix} 1 \\ 0.631 \\ 0.639 \\ 0.555 \\ 0.527 \\ - 0.276 \end{matrix} & \begin{matrix} 0.631 \\ 1 \\ 0.592 \\ 0.503 \\ 0.562 \\ - 0.150 \end{matrix} & \begin{matrix} 0.639 \\ 0.592 \\ 1 \\ 0.912 \\ 0.739 \\ - 0.421 \end{matrix} & \begin{matrix} \begin{matrix} 0.555 \\ 0.503 \\ 0.912 \\ 1 \\ 0.439 \\ - 0.624 \end{matrix} & \begin{matrix} 0.527 \\ 0.562 \\ 0.739 \\ 0.439 \\ 1 \\ 0.165 \end{matrix} & \begin{matrix} - 0.276 \\ - 0.150 \\ - 0.421 \\ - 0.624 \\ 0.165 \\ 1 \end{matrix} \end{matrix} \end{matrix}]$	$[\begin{matrix} \begin{matrix} 1 \\ 0.652 \\ 0.679 \\ 0.642 \\ 0.466 \\ - 0.304 \end{matrix} & \begin{matrix} 0.652 \\ 1 \\ 0.582 \\ 0.509 \\ 0.538 \\ - 0.115 \end{matrix} & \begin{matrix} 0.679 \\ 0.582 \\ 1 \\ 0.848 \\ 0.682 \\ - 0.476 \end{matrix} & \begin{matrix} \begin{matrix} 0.642 \\ 0.509 \\ 0.848 \\ 1 \\ 0.346 \\ - 0.745 \end{matrix} & \begin{matrix} 0.466 \\ 0.538 \\ 0.682 \\ 0.346 \\ 1 \\ 0.221 \end{matrix} & \begin{matrix} - 0.304 \\ - 0.115 \\ - 0.476 \\ - 0.745 \\ 0.221 \\ 1 \end{matrix} \end{matrix} \end{matrix}]$
Gaussian Copula相关参数矩阵	$[\begin{matrix} \begin{matrix} 1 \\ 0.590 \\ 0.618 \\ 0.590 \\ 0.514 \\ - 0.287 \end{matrix} & \begin{matrix} 0.590 \\ 1 \\ 0.560 \\ 0.527 \\ 0.521 \\ - 0.149 \end{matrix} & \begin{matrix} 0.618 \\ 0.560 \\ 1 \\ 0.845 \\ 0.673 \\ - 0.450 \end{matrix} & \begin{matrix} \begin{matrix} 0.590 \\ 0.527 \\ 0.845 \\ 1 \\ 0.474 \\ - 0.730 \end{matrix} & \begin{matrix} 0.514 \\ 0.521 \\ 0.673 \\ 0.474 \\ 1 \\ 0.167 \end{matrix} & \begin{matrix} - 0.287 \\ - 0.149 \\ - 0.450 \\ - 0.730 \\ 0.167 \\ 1 \end{matrix} \end{matrix} \end{matrix}]$	$[\begin{matrix} \begin{matrix} 1 \\ 0.669 \\ 0.700 \\ 0.660 \\ 0.483 \\ - 0.317 \end{matrix} & \begin{matrix} 0.669 \\ 1 \\ 0.600 \\ 0.526 \\ 0.556 \\ - 0.120 \end{matrix} & \begin{matrix} 0.700 \\ 0.600 \\ 1 \\ 0.860 \\ 0.699 \\ - 0.493 \end{matrix} & \begin{matrix} \begin{matrix} 0.660 \\ 0.526 \\ 0.860 \\ 1 \\ 0.360 \\ - 0.760 \end{matrix} & \begin{matrix} 0.483 \\ 0.556 \\ 0.699 \\ 0.360 \\ 1 \\ 0.231 \end{matrix} & \begin{matrix} - 0.317 \\ - 0.120 \\ - 0.493 \\ - 0.760 \\ 0.231 \\ 1 \end{matrix} \end{matrix} \end{matrix}]$
AIC	-2894.1	-3171.4
BIC	-2829.8	-3107.2

下载: 导出CSV

表 3 基于孔压静力触探试验指标的不排水抗剪强度和超固结比概率转换模型

Table 3 Probabilistic transformation models for undrained shear strength and overconsolidation ratio using CPTU test data

直接测量参数	设计参数 $s_{u} / {σ^{'}}_{v}$		设计参数OCR		$s_{u} / {σ^{'}}_{v}$ 与OCR相关系数θ′
直接测量参数	均值 ${μ^{'}}_{1}$	方差 ${σ^{'}}_{1}^{2}$	均值 ${μ^{'}}_{2}$	方差 ${σ^{'}}_{2}^{2}$	$s_{u} / {σ^{'}}_{v}$ 与OCR相关系数θ′
X₃	0.697X₃	0.515	0.600X₃	0.640	0.438
X₄	0.660X₄	0.564	0.526X₄	0.723	0.504
X₅	0.483X₅	0.767	0.556X₅	0.691	0.551
X₆	-0.317X₆	0.900	-0.120X₆	0.986	0.671
X₃, X₄	0.495X₃ + 0.235X₄	0.500	0.567X₃ + 0.0392X₄	0.639	0.440
X₃, X₆	0.714X₃ + 0.035X₆	0.514	0.715X₃ + 0.233X₆	0.599	0.442
X₄, X₆	0.993X₄ + 0.438X₆	0.483	1.03X₄ + 0.664X₆	0.537	0.391
X₅, X₆	0.587X₅-0.452X₆	0.573	0.616X₅-0.262X₆	0.626	0.482
X₃, X₄, X₅	0.223X₃ + 0.403X₄ + 0.182X₅	0.491	0.505X₃-0.187X₄ + 0.505X₅	0.566	0.422
X₃, X₅, X₆	1.07X₃-0.332X₅ + 0.288X₆	0.506	0.682X₃ + 0.0314X₅ + 0.209X₆	0.598	0.447
X₃, X₄, X₅, X₆	0.764X₃ + 0.717X₄-0.474X₅ + 0.715X₆	0.450	0.336X₃ + 0.809X₄-0.129X₅ + 0.691X₆	0.527	0.375

下载: 导出CSV

参考文献(40)

[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.