-
海上风电场作为近年国家发展的重点，目前江苏、浙江、福建、广东、山东海域呈现大规模推进的形势。海洋环境下钻孔取样的扰动较大，室内试验获得的岩土参数可靠性不足。因此，海上工程地质勘察需采用对土体无扰动或轻微扰动的原位测试技术，而孔压静力触探试验(CPTU)因其适用于深水、能获得连续的曲线、精确性和可重复性良好等优点在海上工程勘察中得到广泛应用，在欧洲海上风电场工程中积累了较多的经验。据统计约95%的海洋土原位测试采用了CPT/CPTU(Lunne, 2001)[1]。
海洋静力触探分为海床式(Seabed CPT)、井下式(将探头从钻孔底部压入土中)，也可分为固定平台式和漂浮平台式。海床式静力触探(图1a)采用勘查船将静探仪沉至海床，依靠自身重力作反力贯入地层中，适用于深水，但贯入深度受到一定制约。井下静力触探系统(图1b)将钻探和CPT测试相结合，探头从钻孔底部贯入地层中，可在动态定位的大型综合勘探船上实施，也可在海上固定平台上实施，其贯入深度较少受限。海上固定式平台一般采用自升式平台，也可采用简易支架平台，其适用水深受到较大限制。
CPTU测试成果广泛用于海底地层划分、地层状态判别、土体参数分析、地基承载力估算和海底地质灾害评估。通过孔隙水压力对锥尖阻力进行修正，结合锥身侧阻力对地层进行划分并判断黏性土的稠度状态和无黏性土的密实状态，通过经验公式分析土体的重度、侧压力系数、渗透系数，原位状态下的压缩模量、杨氏模量和小应变剪切模量，黏性土的超固结比、不排水抗剪强度、固结系数，无黏性土的相对密实度和有效内摩擦角等参数。
HTML
-
CPTU探头在饱和土中贯入时，可测得土中的孔压，包括锥面(u1)、锥肩(u2)及摩擦筒尾部(u3)的孔压(图2)。既可测得贯入引起的超孔压，也可测得其随时间的消散过程。
Figure 2. Cone and pore pressure transducer of CPTU
采用锥肩的孔压u2对实测锥尖阻力qc进行修正，见下式(Lunne等，1997[2])。
((1)) ((2)) 式中：qt为修正后的锥尖阻力(kPa)；qc为实测的锥尖阻力(kPa)；u2为锥肩的孔隙水压力(图2，kPa)；a为有效面积比，大部分探头为0.55～0.9，常为0.8；Aa为顶柱的横截面积(图3，cm2)；Ac为锥底的横截面积(图3，cm2)。
Figure 3. Pore water pressure effects on measured parameters
也可采用锥尖的孔压u1，根据经验公式计算u2(Sandven，1990[3])，见下式：
((3)) 式中：u1为锥面的孔压(kPa)；u2为锥尖的孔压(kPa)；u0为静水压(kPa)；K为经验系数。
可采用下式对侧壁摩阻力进行修正(图3)：
((4)) 式中：ft为修正后的侧壁摩阻力(kPa)；u3为套筒尾部位置的孔隙水压力(kPa)；Ast为套筒顶部的横截面积(cm2)；Asb为套筒底部的横截面积(cm2)；As为侧壁摩擦筒表面积(cm2)。
一般情况下Ast＝ Asb时，ft基本与fs相等。
采用下列公式计算摩阻比Rf和孔压参数比Bq：
((5)) ((6)) 式中：σv0为土的竖向总应力(kPa)；u0为静水压(kPa)；Δu为超孔压(kPa)。
-
在不同土类中贯入的锥尖阻力qc、侧壁摩阻力fs、锥肩孔隙水压力u2的响应不同，可据此对地层进行辨别与划分。Robertson等(1986)[4]根据Douglas与Olsen(1981)[5]的土质分类，提出了土的分类图(图4a)，并提出了以qt与Bq的分类图(图4b)，但图4对深度大于30m的土层分类可能产生偏差。
Figure 4. Soil behavior type classification system from CPTU data (after Robertson et al. , 1986)
Robertson(1990)[6]进一步对锥尖阻力qt与摩阻比Rf进行了修正，考虑了上覆压力的影响，采用归一化锥尖阻力Qt、归一化摩阻比Ft、孔压参数比Bq对土进行分类(图5)。
Figure 5. Soil behavior type classification chart based on normalized CPTU data (Robertson, 1990)
Qt与Ft的计算如下：
((7)) ((8)) 式中：σ′v0为土的有效自重应力。
Jefferies与Davies(1993)[7]对Robertson提出的归一化土质分类图(图5)合并为同心圆形式，并采用土行为分类指数Ic对土质分类(表1)，Ic的表达式如下：
土类 Ic值范围 有机质黏土 Ic>3.22 黏土 2.82<Ic<3.22 粉质混合土 2.54<Ic<2.82 砂质混合土 1.90<Ic<2.54 砂土 1.25<Ic<1.90 砾砂 Ic<1.25 Table 1.
Soil behavior type classification based on SBT parameter Ic(Jefferies & Davies，1993) ((9)) Robertson(2009)[8]进一步修正了土行为分类指数Ic，同时Mayne(2014)[9]在Robertson(2009)的基础上细化了土的分类(表2)，修正的土行为分类指数见下式：
划分区域 土类 Ic值范围 1 灵敏性细粒土 N/A 2 有机质黏土 Ic>3.60 3 黏土 2.95<Ic<3.60 4 粉质混合物 2.60<Ic<2.95 5 砂质混合物 2.05<Ic<2.60 6 砂土 1.31<Ic<2.05 7 密实砂土至砾砂 Ic<1.25 8 非常密实的砂土至黏性砂土 N/A 9 非常坚硬的超固结黏土至粉土 N/A Table 2.
Soil behavior type classification based on the modified Ic(Robertson, 2010) ((10)) 表2中1、8、9区的计算方法可参照Mayne(2014)[9]的划分方法，灵敏性细粒土采用式(11)判别，密实砂性土和坚硬黏性土采用式(12)判别。当采用式(12)判定时，1.5%<Ft<4.5%时为区域8，Ft>4.5%时为区域9。
((11)) ((12)) 此外，还有概率分类法(Zhang与Tumay，2003[10])和聚类分析法(Hegazy，1998[11])，工程上应用较少。
-
Robertson(2009)[8]基于无量纲化的锥尖阻力qt/pa和摩阻比Rf提出了土的天然重度计算公式：
((13)) 式中：pa为参考大气压(100kPa)；γ为土的天然重度(kN/m3)；γw为水的重度(kN/m3)。
Mayne(2001)[12]提出土的天然重度与剪切波速、深度的经验公式：
((14)) 式中：Vs为剪切波速(m/s)；z为深度(m)。
-
砂性土的有效内摩擦角主要通过经验或半经验相关关系法(标定罐试验)分析，建立相对密实度与峰值摩擦角的关系。Schmertanan(1978)[13]提出了有效内摩擦角φ′与相对密实度Dr的区间图，如图6所示。
Figure 6. Relationship between φ′ and Dr(Schmertanan，1978)
Baldi等(1986)[14]根据Ticino砂土标定罐试验，建议采用下式分析相对密实度Dr：
((15)) 式中：C0、C1和C2为土的常数；σ′为有效应力，可采用平均有效应力或竖向有效应力。
Jamiolkowski等(2003)[15]在大量标定罐试验的基础上，采用直剪试验得到干砂和饱和湿砂相对密实度的关系：
((16)) ((17)) 式中：K0为侧向土压力系数；Dr(dry)为干砂的相对密实度；Dr(saturated)为饱和湿砂的相对密实度。
Mayne(2006)[16]提出了直接采用修正锥尖阻力qt估算有效内摩擦角的经验公式：
((18)) -
黏性土不排水抗剪强度su的计算可采用理论计算法和经验公式法。理论计算法主要有承载力理论法、孔穴扩张理论法、应力路径理论法等。经验公式法包括直接经验公式法与间接经验公式法。
直接经验公式法包括修正锥尖阻力法、有效锥尖阻力法和超静孔压法，计算公式分别为：
((19)) ((20)) ((21)) 式中：Nkt、Nke、NΔu均为经验系数。
对于修正锥尖阻力法，经验系数Nkt的取值具有一定的区域性，Robertson(2009)[8]建议Nkt取10～18。
对于有效锥尖阻力法，Senneset等(1982)[17]建议Nke＝9±3，但一般认为Nke取值较离散。
对于超静孔压法，从孔穴扩张理论(Vesic，1972[18]；Randolph与Wroth，1979[19]；Campanella等，1985[20])得到NΔu的理论值为2～20。Lunne等(1985)[21]发现NΔu与Bq有较好的相关性，根据室内CU试验反推的NΔu范围为4～10。
刘松玉等(2013)[22]对连云港海相黏性土进行了CPTU试验，认为修正锥尖阻力法对Nkt的拟合效果最好。部分国外学者认为超静孔压法适用于强度较低的黏性土，而软黏土的修正锥尖阻力会出现一定的不确定性。
间接经验公式法一般先估算黏性土的超固结比OCR。Mayne(1991)[23]采用孔穴扩张理论和临界状态理论提出了计算OCR的半经验公式：
((22)) 式中：M为临界状态线的斜率；σ′v0为土体的有效竖向应力。
Lunne等(1997)[24]提出采用归一化锥尖阻力来估算OCR：
((23)) 式中：σ′pc为先期固结压力。k一般为0.2～0.5，平均值为0.3。
Mayne等(2005)[25]给出了估算先期固结压力的简化公式：
((24)) ((25)) 式中：m为土类参数，随细颗粒含量增加而增大；σatm为大气压力(kPa)。
Robertson(2009)[8]在Lunne等(1997)[24]的基础上提出了OCR的估算公式：
((26)) 目前一般采用Lunne等(1997)[24]的k值法估算OCR，再采用SHANSEP法(Stress History and Normalised Soil Engineering Properties)估算不排水抗剪强度su，估算公式如下：
((27)) 式中：a为正常固结土不排水抗剪强度与有效竖向应力比值，一般可取0.25～0.33；b为反映超固结比OCR对黏土不排水强度影响的经验参数，一般可取0.65～1.0。
经验系数a与b在不同区域有较大差异，需事先采用室内模型试验或通过未扰动样的三轴不排水试验标定。
-
土的小应变剪切模量G0的估算可分别采用直接法和间接法，直接法通过CPTU的数据估算，间接法则先计算土的剪切波速Vs。
Pincus等(1993)[26]认为大多数黏性土的小应变剪切模量G0与孔隙比和锥尖阻力存在以下关系：
((28)) Jamiolkowski等(1988)[27]发现土的Dr和原位有效应力的变化控制着qc和G0，对未胶结和未老化的无黏性土，qc和G0呈现一定的相关性。基于标定罐试验与现场测试结果，Rix与Stokoe(1991)[28]提出了G0-qc-σ′v0的关系曲线(图7)，发现归一化
Figure 7. G0-qc-σ′v0 curve(Rix & Stokoe，1991)
采用带剪切波速测试的静力触探时，黏性土和砂性土均可下式估算小应变剪切模量G0：
((29)) 式中：ρ为密度。
Piratheepan(2002)[29]分别给出了砂性土和黏性土剪切波速Vs的估算公式：
((30)) ((31)) 式中：D为泥面以下土体的深度(m)。
Andrus等(2007)[30]对Piratheepan(2002)[29]的公式进行了合并，得到了可用于砂性土和黏性土的剪切波速Vs估算公式：
((32)) 式中：Sf为尺寸效应参数，可取0.92。
3.1 土的天然重度
3.2 砂性土的有效内摩擦角
3.3 黏性土的不排水抗剪强度
3.4 土的小应变刚度
-
某海上风电场位于海岸线约22 km处，基于CPTU测试的数据计算修正锥尖阻力qt与摩阻比Rf，得到不同孔位不同深度的土质分类情况如图8所示，其土类分层及描述如表3所示。图8基于ISO标准划分，而基于我国标准的分类图是今后研究的重点。
Figure 8. Soils classification of a offshore wind farm based on CPTU
土层 深度(海床以下) 土类 描述 顶部软土层 3～9 m 有机质黏土 非常软的黏土层，有机质含量5%～10% 顶部混合土层 10～36 m 粉质黏土 含有黏性土与粉性土、砂性土沉积，沉积物之间存在较多夹层和夹薄层现象，土的强度随着深度增加而增加 粉质黏土混粉砂 黏质粉土 砂质粉土 粉土 粉细砂 粉、黏性土层 37～48 m 淤泥质粉质黏土 含有黏性土与粉性土沉积，沉积物之间存在较多夹层和夹薄层现象，偶尔存在一定的砂土层 粉质黏土 粉质黏土混粉砂 黏质粉土 粉土 砂质粉土 底部软土层 48～56 m 淤泥质粉质黏土 强度与密度较低的粉质黏土 底部混合土层 56～77 m 粉质黏土 粉土与黏土过渡层，随着深度增加粉粒含量增加，逐渐从黏土过渡至粉土，粉质黏土较为多见，常见于砂土层下 粉质黏土混粉砂 粉土 细纱 砂质粉土 Table 3.
Overview of typical soil units and geotechnical characterization 基于CPTU对土体的物理力学参数进行分析，并与室内土工试验的结果对比，结合工程经验提出不同深度的土质参数取值，如表4所示。
序号 层名 层底埋深Z/ m 有效重度γ′/(kN·m-3) 有效内摩擦角φ′/(°) 不排水强度su/kPa 1 粉质粘土 4.7 7.0 — 12 2 粉土夹粉质粘土 5.7 7.4 27 17 3 粉土夹粉砂 6.7 8.3 28 35 4 粉质粘土夹粉土 9.1 8.0 — 30 5 粉砂 17.0 9.2 31 — 6 粉砂 19.6 9.3 32 — 7 粉质粘土 23.2 7.9 — 30 8 粉砂夹粉土 24.5 9.0 31 — 9 粉细砂 25.7 10.0 36 — 10 粉质粘土夹粉土 33.5 8.3 — 40～45 11 粉质粘土夹粉砂 37.8 8.4 — 45～50 12 粉质粘土 49.3 8.3 — 50～55 13 粉质粘土夹粉土 54.3 8.5 — 65 14 粉砂夹粉土 56.2 9.2 34 — 15 粉质粘土夹粉土 59.0 8.8 — 70 16 粉细砂 63.9 10.0 35 — 17 粉质粘土夹粉土 68.6 9.4 — 105 18 粉砂 70.1 9.8 34 — 19 粉质粘土夹粉土 75.0 9.3 — 95 Table 4.
Soils parameters of a offshore wind farm 根据CPTU测试的数据对不同深度土体的浮重度进行计算，如图9所示，可见对浅部松散的无黏性土层，基于CPTU数据的计算值小于室内试验值，而对深部密实的无黏性土层，基于CPTU数据的计算值大于室内试验值，反映了松砂取样振密而密砂取样振松的普遍规律。
Figure 9. Soil unit weight from CPTU data and laboratory test
对于黏性土的不排水抗剪强度，也可采用间接法对CPTU数据进行分析。采用经验参数估算的先期固结压力如图10所示。可见，估算的先期固结压力数据点较为离散，如根据建议的经验参数取值，可能会高估先期固结压力σ′pc，以此计算的不排水抗剪强度可能偏高，因此间接法的经验公式和经验参数取值，也需结合各海域地层的实际情况进行标定和论证。
Figure 10. Preconsolidation pressure based on CPTU data
基于CPTU数据分析的有效内摩擦角和不排水抗剪强度，与室内试验结果的对比如图10所示。对无黏性土，如图11a所示，基于CPTU分析的有效摩擦角与慢剪试验结果基本一致，采用CPTU的经验分析法是可行的。而对于黏性土，如图11b所示，采用修正锥尖阻力Nkt法得到的结果普遍比室内不固结不排水试验(UU)结果大，但小于室内固结快剪强度指标计算的强度值。因钻探取样卸除了原位有效应力，造成试样扰动，使试样进行室内三轴不固结不排水试验获得的不排水抗剪强度偏低，我国某些地基基础设计规范建议对试样进行预固结后再进行室内试验，预固结压力可选择有效自重应力或先期固结压力的一定比例。另一方面，基于CPTU数据估算黏性土的不排水抗剪强度，其经验系数不宜盲目套用国外的经验系数，否则可能高估黏性土的不排水抗剪强度，而应针对本地区的海相黏性土进行专门的标定。在室内采用模型试验，对不同含水率和不同不排水强度的黏性土进行CPTU测试，根据经验公式计算的不排水强度，与十字板、T形贯入仪或球形贯入仪测得的不排水强度进行对比，估算经验参数的范围，是可行的研究方法。
Figure 11. Comparison of effective internal friction angle and undrained shear strength from CPTU data and laboratory test
-
海洋工程勘察中CPTU原位测试技术具有独特的优点，已在我国东部、南部、北部海上风电场勘察中得到成功的应用。
国际上基于CPTU进行地层划分，判断无黏性土的密实状态和黏性土的稠度状态，估算土层的有效重度，分析无黏性土的有效摩擦角和黏性土的不排水抗剪强度等，均已有较为成熟的方法和经验。但由于各海域土层的区域性特征，不宜盲目套用国外的经验公式及其参数。
结合我国各海域的区域性地层特征，针对江苏海域典型的粉土、粉砂，浙江海域典型的软黏土，福建海域典型的花岗岩残积土等区域性土，通过模型试验和多种手段原位测试的对比，对经验公式进行验证或改进，并对经验参数进行标定是切实可行的。同时，基于我国标准的土分类图也需进一步研究。
DownLoad:
