渤海常见土类剪切波速与埋深关系分析1

荣棉水 吕悦军 蒋其峰 修立伟



渤海常见土类剪切波速与埋深关系分析1

荣棉水1）吕悦军1）蒋其峰2）修立伟1）

1）中国地震局地壳应力研究所，北京 100085 2）山东省地震局，济南 250014

土层的剪切波速是岩土地震工程中重要的物理量，本文利用多年来在渤海海域地震安全性评价中积累的资料，研究了渤海常见土类剪切波速和埋深的关系。利用非线性最小二乘法，采用指数函数、一次函数、二次函数、幂函数、“幂函数+常数函数”、“幂函数+一次函数”6种回归模型对各类土的剪切波速和埋深的关系进行了回归分析，以拟合优度以及最小二乘拟合的误差平方和为评价指标对比了各种模型拟合效果的优劣。结果表明，“幂函数+一次函数”回归模型的拟合效果最好。此外，本文给出了该海域常见的7类土在此回归模型下的拟合公式的系数，以供工程中参考。

渤海 剪切波速 埋深 回归分析 拟合

引言

土层的剪切波速在岩土地震工程中有广泛的应用（刘红帅等，2010）。地震工程中，土层剪切波速是确定场地类别、场地地震反应分析、计算动力学参数的重要参数。在工程实践中，土层剪切波速缺失的情况时有发生，这就需要依据经验公式对剪切波速进行科学合理的计算，如何确定合理的经验公式就成为重要的课题。土层的剪切波速与很多参数有关，包括土的成分、结构、类别、密度、埋深等（Hardin等，1972），而研究表明，在某个区域范围内，土层的剪切波速与埋深有很强的相关性，这就为计算剪切波速提供了经验性的手段。曾有不少学者研究过不同地区土体的剪切波速与埋深的经验关系，如周锡元等（1990）对沿海城市地区、程祖锋等（1997）对深圳地区、高玉峰等（2003）对合肥地区、战吉艳等（2009）对苏州城区、苏旭等（2009）对西宁市、邱志刚等（2011）对黑龙江地区、齐鑫等（2012）对下辽河平原、赵淑芳等（2012）对昆明呈贡新区均开展了相关的研究。这些研究表明，不同地区由于土体沉积历史等因素的不同，导致拟合经验公式系数存在很大差异。此外，在剪切波速与埋深的函数关系方面，统计分析过的学者存在不同的认识。如Lew等（1981）曾指出剪切波速与埋深之间呈幂函数关系；周锡元等（1990）给出的沿海城市场地土剪切波速与埋深呈幂函数形式；苏旭等（2009）拟合西宁市剪切波速与埋深的关系时也采用了幂函数形式；战吉艳等（2009）对苏州城区深软场地土的统计采用了线性函数、幂函数分段形式；刘红帅等（2010）曾研究了淤泥质土、碎石土、砂土、粉土、黏性土剪切波速和埋深的关系，认为除淤泥质土外，其他常规土类剪切波速与埋深间的关系用一元二次多项式拟合精度最高；邱志刚等（2011）、齐鑫等（2012）的工作中则采用了包括线性模型、幂函数模型、一元二次多项式模型在内的多种模型；狄圣杰等（2014）则曾用广义回归神经网络预测的方法对剪切波速与埋深的关系进行了探讨。综合以上研究可知，区域土体剪切波速与埋深的经验函数关系并没有一个普遍公认的函数形式，何种函数更为接近实测结果仍值得研究，其结果可能也与所研究的区域有关。渤海蕴藏着丰富的油气资源，随着对资源进行开发利用，渤海中建设了一系列的海洋工程项目，按照《建筑工程抗震设防分类标准》，这些工程建筑需要重点设防。在过去的10年中，渤海海域开展了大量的工程建设，如跨海桥梁、隧道、海洋平台等蓬勃开展，积累了丰富的土层剪切波速资料，具备了研究剪切波速与埋深关系的基础。本文在前人研究成果的基础上将尝试构造新的回归模型，对比筛选出最优的回归模型并给出相应的系数。

1　渤海地质地貌概况

渤海是一个近封闭的内海，为陆架浅海盆地，由辽东湾、渤海湾、莱州湾、渤海中央盆地以及渤海海峡5部分组成。渤海周围有黄河、海河、辽河、滦河等含沙量很大的河流流入，使得渤海水深较浅、地形平缓，海底从莱州湾、渤海湾和辽东湾3个海湾向渤中盆地及渤海海峡倾斜，坡度平缓，平均坡度0.13‰，渤海水深较浅，平均水深18m（陈义兰等，2013）。渤海位于新华夏构造第二沉降带内，是营口隆起、华北凹陷、鲁西隆起与郯庐断裂带的交汇区。渤海分布有厚层第四纪沉积物，在渤中盆地中，第四纪的沉积厚度约600m以上，在渤海湾西岸第四纪的厚度约300m。在更新世期间，渤海曾发生过若干次海进和海退，所以渤海地区更新世地层中海相地层和陆相地层交替出现。在10万年以来的地层中，共发现有3次海侵事件，从下而上分别称为沧州海侵、献县海侵、黄骅海侵。全新世沉积伴随着玉木冰期的最盛时期结束后的全新世海侵形成，在渤中盆地一般厚约8—13m，在辽东湾西侧，其厚度较浅，约2—4m，在渤海西岸一般厚12—14m，表现为灰色—暗灰色沉积层，在黄河和海河口附近海域，受近岸物质的影响，一般表现为黄色沉积层（中国科学院海洋研究所海洋地质研究室，1985）。

2　数据来源

本文所用土层剪切波速和埋深数据来自于中国地震局在渤海海域工程地震安全性评价中近十余年来积累的资料，这些资料涵盖了渤海海域锦州、渤中、南堡、锦西、辽东、秦皇岛、垦利、辽东、岐口等油田33处场地的钻孔剪切波速资料，取样钻孔的分布如图1。这些海域钻孔一般深约120m，钻孔土类以细粒土类为主，包括粉质粘土、粉质细砂、粉砂、粉土、砂质粉土、细砂，其中粉质粘土和粉质细砂分布最广。剪切波速均为获取土样后在实验室测得，埋深为取样位置的深度，资料的可靠性很强。对资料进行统计后，按照土类的不同，对数据进行了分类汇总。根据离散点的趋势对个别异常数据进行了剔除，数据仍余500多组。所统计的土类为渤海常见土类，包括粉质粘土、砂质粉砂、粉砂、粉质细砂、细砂、砂质粉土、粉土共7类。

3　剪切波速测量方法

剪切波速测量装置由三轴压力室、剪切波发射与接收传感器以及DB4型超声测量仪组成，如图2。具体测量过程是：①制备直径3.91cm、高8.0cm的原状土样，套上橡皮膜后安装进三轴压力室内；②根据土样所在土层的上覆有效压力，给土样施加等向固结压力（即图2中围压3与轴向应力1相等），使土样排水固结；③待其固结后，对土样输入剪切波，利用仪器测量出剪切波通过土样的旅行时；④土样的高度除以旅行时即得到土样的剪切波速。

4　数据处理

剪切波速有随深度的增加而增大的规律，不同的学者对不同地区的研究不尽相同，现在已有的经验公式的形式有指数函数、一次函数、二次函数、幂函数、“幂函数+常数函数”（胡钧等，1997）。其回归模型具体如下：

（1）指数函数模型：=e；

（2）一次函数模型：=+；

（3）二次函数模型：=2++；

（4）幂函数模型：=ah；

（5）“幂函数+常数函数”模型：=ah+；

其中，、、为拟合系数，为土体埋深。作者首先按土类分类统计了渤海常见土类的剪切波速和相应的埋深，然后分别用上述5种回归模型对剪切波速和埋深间的关系进行回归分析，并计算相应的拟合优度和最小二乘法拟合的误差平方和，其结果示于表1。拟合优度是用来度量变量在多大程度上能够被所估计的回归方程所解释的指标。对于一组样本观测值（X，Y），=1，2，…，，若得到以下样本回归线，以一次函数举例：

则的第个观测值与样本均值的离差，可分解为2个部分：

形式，则ESS、RSS和交叉项的值会发生变化。已有文献研究（徐秦等，1992）证实，若拟合曲线为线性函数形式，如式（1）或回归模型（2）等，交叉项为零。若拟合曲线为非线性函数，则交叉项为正、负或零均有可能。一般情况下，若实际观测点离样本回归线越近，则ESS在TSS中所占比重越大，RSS在TSS中占比越小。本文将拟合优度2定义如下：

2=1-RSS/TSS （4）

这一拟合优度定义主要出于以下几点考虑：①线性与非线性曲线通用的拟合优度指标虽引起了一些学者（徐秦等，1992；张世强，2002）的注意，但简单、合理适用的指标仍在讨论过程中；②解决实际问题时，所拟合的回归曲线优劣标准为，在相同参数个数条件下，拟 合值越接近观测值，则认为拟合得越好。若把{Y}与{}视为维空间中的点，若{Y}与{}距离越小，则可认为用{}估计{Y}越好。式（4）中RSS则恰好表示维空间中{Y}与{} 的距离，具有清晰的几何意义。RSS/TSS则体现了相对误差的概念。文献（张世强，2002） 中将RNL=1-做为线性与非线性曲线通用的拟合优度指标，但作者认为，这一定 义与式（4）并无本质区别；③在线性回归中，式（4）即为统计学中常用的可决系数；④RSS主要表示回归线不能解释的部分，ESS与交叉项均与回归线能解释的部分有关，则式（4）可以看做总体平方和之中，与回归线能解释的部分有关的成分所占的比例。

在式（3）中交叉项大于或等于零的情况下，2（此时其范围为0—1）越接近1，说明实际观测点离样本越近，拟合优度越高。此外，本文回归到用最小二乘法确定拟合函数中各个系数时采用观测值与拟合值误差平方和最小的原则上来，也计算了每一个2相应的RSS值，作为评价拟合效果的另一个指标。回归模型（1）—（5）的拟合结果分别示于图3—图9中，表示拟合效果的数值列于表1。表1中的数值表明，在拟合参数个数相同的情况下，幂函数拟合效果最佳。此外，从图中不难发现，当埋深较浅时，如<20m时，剪切波速随埋深的变化呈现幂函数变化特征；当埋深较大时，如>40m时，剪切波速随埋深变化更接近一次函数特征。因此，在前人研究的基础上，作者综合幂函数与一次函数的特点，构造了新的“幂函数+一次函数”的表达式形式，以期获得更好的拟合效果。函数模型如下：

（6）“幂函数+一次函数”模型：=ah++其中，、、、为拟合系数，为土体埋深。本文选用Matlab绘制统计数据的散点图，并利用fit函数的非线性最小二乘法拟合功能对数据进行拟合，计算出回归曲线的系数、拟合优度2、误差平方和，最后，对比各回归模型的拟合优度，并给出最优的回归模型的系数，以供工程中参考。

5　数据处理结果与分析

经过程序处理和计算，得到了7类土分别在6种回归模型下的剪切波速-埋深的表达式，以及相应的拟合优度、误差平方和（表1）。从表中看出，从模型（1）至模型（6），拟合优度越来越接近1，同时误差平方和逐渐减小，表明从模型（1）至模型（6）拟合效果越来越好，而模型（6）的拟合效果最好。笔者也曾尝试4参数的三次函数或更高次函数形式的拟合，难以获得较好拟合效果。综合以上拟合结果可知本文构造的“幂函数+一次函数”的回归模型具有很好的适用性。为直观地对比回归结果的优劣，本文给出粉质粘土、粉质细砂等7类常见土类在回归模型（1）—（6）下的拟合结果（图3—图9）。可以看出模型（1）至模型（3）对埋深较浅的土样的剪切波速拟合效果不理想，模型（4）对埋深较浅的土样的剪切波速的拟合效果很好，而埋深较大时拟合效果不佳，在这一方面，模型（5）、模型（6）依次有了更进一步的改进，模型（6）对试验数据具有最佳的拟合效果。渤海海域各类土在模型（6）下的各项系数见表2。

表1　各类土在各回归模型下的拟合优度

续表

注：2表示拟合优度，RSS表示误差平方和，其单位为m2/s2。

表2　各土类采用v=ahb+ch+d形式拟合的系数

6　结论

本文统计了渤海海域常见土类的剪切波速和埋深数据，研究了剪切波速与埋深之间的经验关系，构造了“幂函数+一次函数”形式的回归模型。通过与其他形式的回归模型拟合优度及误差平方和的对比，证实该回归模型具有很好的适用性。此外，本文给出了渤海各常见土类在该回归模型下的系数。本文的工作可为海域工程中缺失剪切波速数据的场地提供了重要的参考依据，也为钻至一定深度后终孔并推测下方土层剪切波速提供了有效的推测手段。

致谢：感谢匿名审稿专家提出的中肯的修改意见和建议。本文的基础数据多引用自笔者参与和收集的所在单位的多项海域建设项目工程场地地震安全性评价报告，作者向参与这些海域项目的技术人员表示感谢，如有引用不当之处还请谅解。

陈义兰，吴永亭，刘晓瑜等，2013．渤海海底地形特征．海洋科学进展，31（1）：75—82．

程祖锋，李萍，李燕等，1997．深圳地区部分岩土类型剪切波速与深度的关系分析．工程地质学报，5（2）：163—168．

狄圣杰，汪明元，张昆等，2014．江苏近海地层原位剪切波速相关特性及预测方法研究．海洋学报，36（3）：127—133．

高玉峰，刘汉龙，2003．合肥膨胀土剪切波速的特征分析．岩土工程学报，25（3）：371—373．

胡钧，杜坚，1997．上海场地土的脉动特性及剪切波速特征．物探与化探，21（4）：305—307，299．

刘红帅，郑桐，齐文浩等，2010．常规土类剪切波速与埋深的关系分析．岩土工程学报，32（7）：1142—1149．

齐鑫，丁浩，2012．下辽河平原区剪切波速与土层埋深关系分析．世界地震工程，28（3）：151—156．

邱志刚，薄景山，罗奇峰，2011．土壤剪切波速与埋深关系的统计分析．世界地震工程，27（3）：81—88．

苏旭，杨丽萍，邹海宁，2009．西宁市区土层剪切波速随深度变化特征．高原地震，21（1）：12—16．

徐秦，薛茜，徐睿，1992．浅论曲线拟合中的相关指数2．中国卫生统计，9（6）：44—45．

战吉艳，陈国兴，刘建达，2009．苏州城区深软场地土剪切波速与土层深度的经验关系．世界地震工程，25（2）：11—17．

张世强，2002．曲线回归的拟合优度指标的探讨．中国卫生统计，19（1）：9—11．

赵淑芳，杨宏亮，2012．场地剪切波速的特征分析．云南大学学报（自然科学版），34（S2）：267—271．

中国科学院海洋研究所海洋地质研究室，1985．渤海地质．北京：科学出版社，1—11．

周锡元，王广军，苏经宇，1990．场地·地基·设计地震．北京：地震出版社，28—172．

Hardin B. O., Drnevich V. P., 1972. Shear modulus and damping in soils: design equations and curves. Journal of Soil Mechanics & Foundations Division, 98(7): 667—692.

Lew M., Chieruzzi R., Campbell K. W., et al, 1981. Correlations of seismic velocity with depth. In: Proceedings of the 1st International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics (ICRA GEESD) U.S.A. St. Louis, Missouri, 1981. University of Missouri—Rolla.

Relationship between Shear Wave Velocity and Depth for Conventional Soils in the Bohai Sea Area

Rong Mianshui1), Lv Yuejun1), Jiang Qifeng2)and Xiu Liwei1)

1) Institute of Crustal Dynamics, China Earthquake Administration, Beijing 100085, China 2) Earthquake Administration of Shandong Province, Ji'nan 250014, China

Shear wave velocity of soil layers is an important physical index in geotechnical earthquake engineering. Based on the experimental data accumulated over the years in the seismic safety evaluation research projects of the Bohai sea area, relationship between shear wave velocity and depth of conventional soils has been studied. In order to achieve the best fitting of experimental data, 6 regression models including exponential function, linear function, quadratic function, power function, power function combined with constant function, power function combined with linear function, have been used for regression analysis through nonlinear least square fitting method. The goodness-of-fit has been used as an index to evaluate the matching effect. The results showed that, the proposed form of function which combining power function with linear function has the best matching effect for experimental velocity data. Furthermore, the fitting coefficients of proposed function for 7 kinds of conventional soils in the Bohai Sea area have been presented for reference in engineering.

The Bohai Sea; Shear wave velocity; Depth; Regression analysis; Fitting

1基金项目 中央公益型院所基本业务专项项目资助（ZDJ2014-07、ZDJ2015-02）

2015-12-19

荣棉水，男，生于1982年。博士，高工，硕士生导师。主要从事场地地震反应分析、强震动观测等方面的研究。E-mail：waltrong@126.com

荣棉水，吕悦军，蒋其峰，修立伟，2017．渤海常见土类剪切波速与埋深关系分析．震灾防御技术，12（2）：288—297． doi：10.11899/zzfy20170205