Study on VS30 prediction of Datong Basin Based on the Analysis of Typical Geological Characteristics
-
摘要: 土体剪切波速是进行土层地震反应分析的动力学参数，对场地地震动参数确定具有重要意义。基于地质地貌分析，将大同盆地划分为5类典型地质单元。对盆地1429个钻孔剪切波速资料进行分析，探讨VS30与VS20的相关性，研究土体埋深、岩性、地质单元、标贯击数及密实度等地质特征对VS的影响，并基于地质单元、剪切波速比、密实度系数及第四系上部覆盖层厚度相关性分析给出土体VS30预测模型。研究结果表明，基于典型地质特征的VS30预测模型拟合优度R2>0.90，预测精度很高，对于离散性较大、直接拟合估算较差及无剪切波速场地来说，以区分地质单元及土体类型的方式进行VS30分解预测是良好的研究思路。首次在区分地质单元及土体类型的前提下提出剪切波速比及密实度系数，并将其与第四系上部覆盖层厚度综合应用于VS30预测研究。研究结果可为大同盆地城市防震减灾规划、震害预测、区域性地震安全评价提供重要技术支撑。Abstract: The shear wave velocity of soil is a dynamic parameter for seismic response analysis of soil layer,which is of great significance to determine the ground motion parameters of the site. The shear wave velocity data of 1429 boreholes in the basin are analyzed,the correlation between VS30 and VS20 is discussed,and the influence of geological characteristics such as soil burial depth,lithology,geological unit,SPT blow number and compactness on VS is studied. Based on the correlation analysis of geological unit,shear wave velocity ratio,compactness coefficient and thickness of Quaternary upper overburden,the soil VS30 prediction model is given. The results show that the goodness of fitting of the VS30 prediction model based on typical geological characteristics is greater than 0.90,and the prediction accuracy is very high. For the sites with large discreteness,poor direct fitting estimation and no shear wave velocity,it is a good research idea to distinguish geological units and soil types for VS30 decomposition prediction. For the first time,the shear wave velocity ratio and compactness coefficient are proposed on the premise of distinguishing geological units and soil types,and they are integrated with the thickness of Quaternary upper overburden for VS30 prediction. The research results can provide important technical support for urban earthquake prevention and mitigation planning,earthquake disaster prediction and regional earthquake safety evaluation in Datong Basin.
-
引言
土体剪切波速是重要的动力学参数，是局部场地条件的直接反映，能表征地层的软硬状态及相对密实程度，进行剪切波速研究对场地类别划分及地震动参数确定具有重要意义。20世纪90年代以来，随着城市发展进程加快，学者们开始采用地下30m的等效剪切波速VS30考虑场地效应，进而确定场地地震动影响(吕悦军等，2008；Abrahamson等，2008；黄雅虹等，2009；陈鲲等，2010；史大成等，2012)。对于波速测量深度小于30m的钻孔，部分学者尝试利用浅层波速外推深层波速，从而进行VS30预测(彭艳菊等，2009；喻畑等，2015；Xie等，2016；江志杰等，2018)。对于未进行剪切波速测试的地区，国外部分学者尝试建立土层年代、岩性、沉积环境等地质要素与VS30的相关性(Wald等，2006；Wills等，2015)，目前国内针对无剪切波速测试场地VS30的研究不多，因此如何确定土层地质年代、地层岩性、沉积环境、覆盖层厚度等地质特征与VS30的相关性，并建立适用于本地区的VS30经验预测模型是值得深入探讨的问题。
大同盆地近年来已进行了大量地质勘察、地震安全性评价及地震小区划工作，积累了丰富的钻孔测试数据。对大同盆地进行典型地质单元划分，通过分析岩性、地质单元、标贯击数、密实度、覆盖层厚度等典型地质特征，给出不同地质单元土体剪切波速分布特征及土体与埋深的统计关系，基于典型地质特征分析给出盆地VS30预测模型。研究结果可为大同盆地大震情景构建、地震危险性分析、震害预测、区域性地震安全评价提供基础数据，为震后快速评估提供技术支撑，并可作为大同地区建筑场地类别划分及城市防震减灾规划的重要参考依据。
1. 大同盆地新构造分区及第四系沉积特征
大同盆地位于鄂尔多斯地块东北缘，属汾渭断陷带次级构造单元(王乃梁等，1996；瞿伟等，2013)。盆地走向NE-NNE，长约225km，宽约60km，盆地内发育了大同-怀仁凹陷、马营庄凹陷等5个次级构造单元及3个沉降中心，其中马营庄凹陷新生界最大沉积厚度为3500m，第四系最大沉积厚度为900m，是盆地的沉降中心；朔州断阶新生界最大沉积厚度为350m，第四系最大沉积厚度为100m，是断陷最浅小的构造单元(王乃梁等，1996；张世民等，1997)，如图 1所示。
图 1 大同盆地典型地质单元分区及钻孔分布图Figure 1. Typical geological units and borehole distribution of Datong Basin大同盆地第四系沉积表现出明显的同时异相特征，从盆地周边的丘陵、洪积扇及倾斜平原向盆地中心过渡为冲湖积平原，由于桑干河及其支流的搬运，沉积物分布表现出显著的规律性，总体上盆地中心粉土、粉质粘土分布较多，由盆地中心向边缘过渡砾石粒径逐渐增大，砂砾石层逐渐增厚，直至山前洪积扇地带转变为卵砾石(李铁锋等，1993)。
2. 区域地质单元划分及钻孔波速资料分析
作者收集了近十年大同盆地内实施的地质钻孔，钻孔涉及地震小区划、地震安全性评价、地质勘察及其他相关项目，其中孔深30-50m的地质勘察钻孔542个，孔深50-80m的地震安全性评价及地震小区划钻孔872个，孔深80-100m的地震小区划钻孔2个，孔深100-150m的水文地质、煤田地质及地震小区划钻孔13个，共计1429个。由图 1可知，大多数钻孔位于大同市、朔州市及县城周边，大致以城区为中心呈不均匀发散状分布，其中大同-怀仁凹陷、朔州断阶构造单元钻孔密度最大，黄花梁凸起构造单元钻孔密度最小，钻孔分布特征如表 1所示。
表 1 大同盆地典型地质单元钻孔分布特征Table 1. Borehole distribution characteristics of typical geological units of Datong Basin地质单元类型 钻孔数量 VS30优势分布区间/(m·s-1) 钻孔分布特征 Qhal-l 292 224.5-272.6 分布于盆地中心Qh冲湖积平原及市区Ⅰ级阶地，钻孔分布面积面大，其中山阴、应县、怀仁及浑源县城钻孔密度较大 Qp3al-l 274 250.8-294.7 分布于市区河谷两岸，地貌上为Ⅱ级阶地，与Qp3al单元前缘相邻，除朔州市区东南部以外，盆地内该单元钻孔分布密度最大 Qhal+pl 221 268.4-306.9 分布于大同和朔州市区河床及山前倾斜平原前缘，地貌上与Qp3pl前缘及Qp3al后缘相邻，山阴、怀仁及浑源县城钻孔分布较密 Qp3al 423 274.6-332.4 分布于大同及朔州市区Ⅱ级以上高阶地，地貌上处于Qhpl及Qp3pl单元前缘，除市区钻孔分布稀疏外，其余地段钻孔密度均较大 Qp3pl 219 288.7-363.5 分布于山前倾斜洪积平原，地貌上与Qhpl及Qp3al单元后缘相邻，除市区钻孔分布较密外，其余地段钻孔密度较小 根据地质地貌及土体结构特征，遵循由盆地中心向边缘过渡的地质单元空间顺序变化规律，将大同盆地划分为5类典型地质单元，划分时主要遵循以下原则：盆地中心向盆地边缘过渡，地质单元由冲湖积平原向冲积平原及洪积倾斜平原过渡，上部覆盖层由全新统转变为上更新统，土体粒径具有逐渐增大的趋势，第四系上部土层结构表现为由软变硬的显著规律性。本文给出的分类方法覆盖了大同盆地除Qp3eol外的所有工程地质单元，由于基于地质地貌分区给出地质单元分类结果，基本反映了不同工程地质单元的土体差异性，即由盆地中心向盆地边缘过渡，岩土体物理性质呈现出明显不同的工程地质特征，因此本文给出的地质单元分类结果与盆地内工程地质单元类型具有良好的匹配性。
计算钻孔等效剪切波速VS20及VS30，如图 2所示，二者拟合优度R2=0.9289，表明VS20与VS30之间存在高度正相关性，对于已有20m测试数据的钻孔，其VS30可直接由VS20求得，且拟合精度较高。将VS30按5类地质单元统计后，给出VS30与地质单元关系图，由图 3可知，不同地质单元的VS30有不同的优势分布区间，其中Qhal-l单元VS30区间范围最小，其优势范围为224.5m/s≤VS30≤272.6m/s，反映了该单元土体地质年代较新、颗粒度较小、密实度较低，VS30变化幅度相对较小；而Qp3al单元VS30区间范围最大，其优势范围为288.7m/s≤VS30≤363.5m/s，表明该地质单元土体地质年代较老、颗粒度较大、实度较高，VS30变化幅度相对较大。5类单元VS30数值区间总体近似呈指数分布，R2=0.8232，表明VS30与地质单元具有较高的正相关性。
3. 土体剪切波速与典型地质特征的关系
剪切波速与土体类型、年代、埋深、密实度和颗粒度等地质特征有关(汪云龙等，2016；沈方铝等，2018；王琦等，2018)，本节重点分析土体剪切波速与埋深及密实度的关系。
3.1 大同盆地剪切波速与土体埋深的关系
基于盆地37295组剪切波速数据绘制了土体剪切波速随埋深的变化曲线，结果表明土体埋深为35-80m时，相比于二次多项式拟合结果，幂函数拟合预测结果总体比剪切波速实际值偏低，二次多项式拟合效果优于幂函数，综合对比后给出土体埋深与剪切波速之间的二次回归关系式，如图 4、表 2所示，本次回归分析利用了大同盆地大量基础测试资料，回归公式精度较高。分析可知卵砾石剪切波速最大，粉细砂剪切波速最小；粉质粘土与中粗砂剪切波速曲线相近，粉土与粉细砂剪切波速曲线相近；25-60m埋深区间内，中粗砂、粉质粘土、粉土及粉细砂4组剪切波速曲线呈逼近状态；卵砾石剪切波速随土体埋深离散性最大，总体拟合效果最差，而粉细砂剪切波速随土体埋深离散性最小，总体拟合效果最好。基于回归公式可实现不同埋深土体剪切波速的快速估计。
图 4 不同土体剪切波速-埋深关系图Figure 4. Relationship between shear wave velocity and buried depth of different soils表 2 不同土体剪切波速与埋深的相关性Table 2. Correlation between shear wave velocity and buried depth of different soils土体类型 VS=ah2+bh+c 拟合优度R2 a b c 卵砾石 -0.0299 6.7326 227.92 0.7591 中粗砂 -0.0275 6.4253 215.77 0.8755 粉质粘土 -0.0235 5.9766 216.12 0.8256 粉土 -0.0312 6.5989 197.99 0.8337 粉细砂 -0.0371 7.2017 178.39 0.9392 3.2 大同盆地剪切波速与土体密实度的关系
对盆地4310组标贯数据按照地质成因进行归类，给出了标贯击数、密实度与剪切波速关系图，见图 5-7。其中卵砾石、粉土、砂类土密实状态界限值参照《岩土工程勘察规范》(GB 50021-2001)(中华人民共和国建设部等，2009)，粉质粘土密实状态界限值参照《工程地质手册》(《工程地质手册》编委会，2018)。分析可知随标贯击数及密实度的增大，剪切波速呈递增趋势，洪积单元剪切波速增速最大，波速变化区间最大，冲湖积单元增速最小，波速变化区间最小；标贯击数较小时，剪切波速与标贯击数的关系表现为曲线型，随着标贯击数的增加或土体密实度的增大，剪切波速与标贯击数的关系近似为直线型，这种变化可用幂函数表示。
图 5 洪积单元剪切波速与标贯击数关系图Figure 5. Relationship between shear wave velocity and standard penetration number of diluvial units
图 6 冲积单元剪切波速与标贯击数关系图Figure 6. Relationship between shear wave velocity and standard penetration number of alluvial units
图 7 冲湖积单元剪切波速与标贯击数关系图Figure 7. Relationship between shear wave velocity and standard penetration number of alluvial-lacustrine units对剪切波速与标贯击数进行回归分析，如表 3所示，对幂函数模型进行拟合优度判定及F显著性检验，其中R2区间为[0.716，0.850]，表明拟合度一般-良，F检验Sig值均小于0.0001，表明模型差异性极显著。分析可知相同地质单元内不同土体密实状态一致时，其剪切波速近似呈一定比例，以粉土为基准，给出其他土体与粉土的平均剪切波速比。上述分析结果表明土体剪切波速与标贯击数及密实度存在良好的正相关性，可用式(1)近似表示三者之间的关系。
$$ {V_{\rm{s}}} = \left({\left\{ {\begin{array}{*{20}{c}} {1.09(卵砾石)}\\ {1.04(中粗砂)}\\ {1.02(粉质粘土)}\\ {1.0(粉土)}\\ {0.97(粉细砂)} \end{array}} \right. ~~~~+ \left\{ {\begin{array}{*{20}{c}} {0.01({\rm{pl}}, {N_{63.5}} \le 10;{\rm{al - l}}, {N_{63.5}} \ge 30)}\\ {0({\rm{pl, al - l}}, 10 < {N_{63.5}} < 30)}\\ { - 0.01({\rm{pl}}, {N_{63.5}} \ge 30;{\rm{al - l}}, {N_{63.5}} \le 10)} \end{array}} \right.} \right)\left({\begin{array}{*{20}{c}} {116.33{N_{63.5}}^{0.358}({\rm{pl}})}\\ {122.52{N_{63.5}}^{0.315}({\rm{al}})}\\ {108.86{N_{63.5}}^{0.321}({\rm{al - l}})} \end{array}} \right) $$ (1) 表 3 剪切波速回归分析及相同密实状态下的剪切波速比Table 3. Regression analysis of shear wave velocity and shear wave velocity ratio under the same compaction state地质单元 回归项 土体类型 卵砾石 中粗砂 粉细砂 粉质粘土 粉土 冲湖积单元al-l k1 110.42 109.40 105.23 104.02 108.86 k2 0.345 0.333 0.322 0.343 0.321 拟合优度R2 0.843 0.764 0.850 0.804 0.741 Sig(F显著性检验) ＜0.0001 波速比区间 1.048-1.106 1.021-1.049 0.968-0.970 0.986-1.035 1.000 波速比均值 1.087 1.039 0.969 1.018 1.000 冲积单元al k1 119.34 123.40 106.51 120.91 122.52 k2 0.346 0.320 0.344 0.326 0.315 拟合优度R2 0.815 0.780 0.825 0.742 0.716 Sig(F显著性检验) ＜0.0001 波速比区间 1.035-1.113 1.032-1.044 0.929-0.995 1.014-1.028 1.000 波速比均值 1.088 1.040 0.973 1.023 1.000 洪积单元pl k1 130.86 123.23 115.37 120.04 116.33 k2 0.348 0.352 0.353 0.356 0.358 拟合优度R2 0.770 0.755 0.762 0.764 0.772 Sig(F显著性检验) ＜0.0001 波速比区间 1.084-1.109 1.036-1.051 0.975-0.985 1.028-1.030 1.000 波速比均值 1.092 1.042 0.978 1.029 1.000 对土体剪切波速进行分段统计，发现平均剪切波速呈现一定比例关系，如表 4所示，除密实状态卵砾石、砂类土与粘性土统计结果相差较大外，松散-中密的卵砾石、砂类土及粘性土统计结果基本一致，这是由于卵砾石标贯击数界限值与砂类土及粘性土不同，卵砾石在标贯击数较小时优先进入密实状态，其剪切波速增速远大于其他土体。
表 4 不同密实状态土体剪切波速比Table 4. Soil shear wave velocity ratios in different compaction states土体类型 土体密实状态 松散(流-软塑) 稍密 可塑 中密(硬塑) 密实(坚硬) 卵砾石 区间 0.64-0.78 0.80-0.922 0.88-0.98 1.00 1.21-1.28 均值 0.70 0.86 0.94 1.00 1.26 其它土体 区间 0.63-0.78 0.79-0.92 0.90-0.99 1.00 1.16-1.18 均值 0.69 0.86 0.94 1.00 1.16 4. 土体等效剪切波速VS30预测分析
4.1 土体等效剪切波速VS30预测模型
对于大同盆地来说，孔深小于100m的钻孔达不到揭穿第四系覆盖层厚度的要求，但却能揭示第四系上部各统地层分界面，本文在收集此类钻孔的基础上，详细统计了第四系上部各统地层分界面埋深数据，基于地质单元类型、土体平均密实度、第四系上部覆盖层厚度及基岩覆盖层厚度给出土体30m等效剪切波速VS30预测模型，如下式：
$$ \begin{array}{l} \lg {V_{S30}} = \lg {\lambda _{(l)}} + \lg {D_{30(i, j)}} + \lg {\rm{(}}{S_{(H)}}{\rm{(}}{D_{30(i, j)}}\\ {\rm{ + 1) + }}{\sigma _{(i)}}({D_{30(i, j)}} - 1)) - 0.301 \end{array} $$ (2) 式中，VS30为30m等效剪切波速预测值；λ(l)为基岩覆盖层厚度影响系数，由表 5获得；D30(i, j)为30m厚度平均密实度系数，由式(4)获得；S(H)为30m等效剪切波速预估初始值，由式(5)获得；σ(i)为不同地质单元统计样本标准偏差，由表 6获得。
表 5 基岩覆盖层厚度影响系数Table 5. Influence coefficient of bedrock overburden thickness基岩覆盖层厚度l/m > 60 50-60 40-50 30-40 < 30 λ(l) 1.00 1.01 1.03 1.05 1.08 表 6 S(H)回归分析Table 6. Regression analysis of S(H)编号 地质单元类型 统计样本 回归系数 拟合优度 F显著性检验 样本量 标准差 a b c R2 F Sig 1 Qhal-l 292 7.48 -0.018 0.554 259.47 0.724 378.29 ＜0.0001 2 Qp3al-l 274 8.14 -0.006 -0.264 291.38 0.758 409.55 3 Qhal+pl 221 8.64 -0.002 -0.607 308.78 0.812 423.49 4 Qp3al 423 11.18 -0.006 -0.530 331.05 0.732 579.50 5 Qp3pl 219 16.56 -0.007 -0.969 365.98 0.769 352.33 (1) 基岩覆盖层厚度影响系数λ(l)
根据《建筑抗震设计规范》(GB 50011-2010)(中华人民共和国住房和城乡建设部，2016)条文，由基岩覆盖层厚度L、等效剪切波速VS20及特征周期Tg三者之间的关系图(见图 8)可知，随基岩覆盖层厚度的减小，场地类别由Ⅳ类变为Ⅰ类，特征周期由0.55s递增至0.25s，Tg减小意味着VS20增大，即土层VS20随着基岩覆盖层厚度的减小而增大，尤其在场地类别转变的临界点，其变化趋势最明显。本文基岩覆盖层厚度均指钻孔揭示的中等风化基岩面埋深，基岩覆盖层厚度影响系数如表 8所示。
(2) 土体平均密实度D30(i, j)
土体相对密实度系数及平均密实度系数见式(3)、(4)。
$$ \beta (i, j) = {X_i}{Y_j} = \left({\left\{ {\begin{array}{*{20}{c}} {1.09}\\ {1.04}\\ {1.02}\\ {1.0}\\ {0.97} \end{array}} \right. + \left\{ {\begin{array}{*{20}{c}} {0.01({\rm{pl}}, {N_{63.5}} \le 10;{\rm{al - l}}, {N_{63.5}} \ge 30)}\\ {0({\rm{pl}}, {\rm{al - l}}, 10 < {N_{63.5}} < 30)}\\ { - 0.01({\rm{pl}}, {N_{63.5}} \ge 30;{\rm{al - l}}, {N_{63.5}} \le 10)} \end{array}} \right.} \right)\left({\begin{array}{*{20}{c}} {{\rm{0}}{\rm{.65}}}(流塑粘土)\\ {{\rm{0}}{\rm{.75}}}(松散、软塑)\\ {{\rm{0}}{\rm{.86}}}(稍密)\\ {{\rm{0}}{\rm{.94}}}(可塑)\\ {{\rm{1}}{\rm{.00}}}(中密、硬塑)\\ {{\rm{1}}{\rm{.16}}}(密实、坚硬)\\ {{\rm{1}}{\rm{.26}}}(密实卵石) \end{array}} \right) $$ (3) $$ \lg {D_{30(i, j)}} = \lg \left({\sum\limits_{i = 1}^5 {\sum\limits_{j = 1}^6 {{\beta _{(i, j)}}{h_{(i, j)}} + 30} } } \right) - 1.778 $$ (4) 式中，D30(i, j)为30m土体平均密实度系数，h(i, j)为i类土体j状态下的土体厚度，β(i, j)为不同土体相对密实度系数。
(3) 土体30m等效剪切波速预估初始值S(H)
第四系上部覆盖层厚度指由钻孔揭示的第四系中更新统及以上分层厚度，作者首次引入第四系上部覆盖层厚度，并根据钻孔标志层埋深给出第四系上部覆盖层厚度统计结果，本次统计表明大同盆地第四系上部覆盖层厚度H取值区间为[10.8，60.1]。其VS30可由式(5)确定：
$$ {S_{\left(H \right)}} = a{H^2} + bH + c $$ (5) 式中，S(H)为30m等效剪切波速预估初始值；H为第四系上部覆盖层厚度；a、b、c为回归系数；R2为拟合优度。
图 9、表 6所示为5个地质单元的VS30统计结果，Sig＜0.0001，远小于显著性水平0.01，说明5个模型中的S(H)与H之间的差异性极显著，R2为(0.724，0.812)，说明两者中不被解释的变量较少，虽然在总体样本情况下VS30随H的变化呈现一定离散区间，二者拟合情况较差，但在区分地质单元类型的情况下二者呈明显负相关性，拟合程度较好。分析可知随着H的增大，不同单元土体VS30均呈减小趋势，其中Qp3pl单元VS30递减速率最大，Qhal-l单元VS30递减速率最小。
图 9 土体VS30与第四系上部覆盖层关系Figure 9. Relationship between VS30 and Quaternary upper overburden4.2 土体VS30模型预测结果误差分析
由图 9可知土体剪切波速具有明显的离散性，虽在横向上VS30与H具有一定的拟合效果，但纵向上看，即不考虑地质单元类型时，VS30与H拟合效果较差，相同的H对应于5类地质单元的VS30区间值，此区间范围＞100m/s，此时仅能根据VS30随H的变化趋势，给出二者之间上下限拟合曲线，上下限曲线之间整个区间的波速变化规律无法深入研究，若不考虑VS30随地质单元类型的空间变化，研究结果会造成巨大偏差。基于地质单元分解预测思路，给出每个单元及所有单元的分解预测结果。表 7、图 10为基于地质单元、岩性、H、N63.5及D30(i, j)按式(2)进行统计后的结果，其中岩性、H、λ(l)均由钻孔获得，D30(i, j)由标贯试验测得，886个钻孔直接基于钻孔揭示的地层及N63.5求得VS30，其余无N63.5数据的钻孔参照邻孔相应测试数据类比求得VS30。
表 7 预测值VS30与实测值VSE30拟合优度及残差描述统计Table 7. Fitting goodness and residual descriptive statistics of predicted value VS30 and measured value VSE30编号 地质单元类型 拟合优度R2 均方根误差
RMSE残差δ 平均值M 标准误差Se 标准差σ 置信度(95.0%) 1 Qhal-l 0.828 4.584 0.031 0.269 4.592 0.529 2 Qp3al-l 0.843 3.660 -0.035 0.219 3.486 0.432 3 Qhal+pl 0.904 3.523 -0.028 0.237 3.531 0.468 4 Qp3al 0.824 5.834 0.155 0.282 5.843 0.555 5 Qp3pl 0.865 7.100 -0.047 0.548 7.134 1.081 6 最终预测结果 0.973 5.376 0.010 0.142 5.379 0.278 
图 10 大同盆地不同地质单元土体VS30预测结果Figure 10. Prediction results of VS30 for different geological units of Datong Basin由表 7、图 10可知，5类地质单元VS30与VSE30的R2区间为[0.824，0.904]，R2＞0.80，分解预测最终R2为0.973，R2＞0.90，说明预测模型具有较高-极高的拟合优度。RMSE区间为[3.523，7.100]，表明二者偏差随VS30的增加而增大，但均小于8.0。预测模型最终RMSE为5.376，表明误差集中分布于5.376附近，模型预测效果较好。分析残差结果δ可知，M区间为[-0.035，0.155]，Se区间为[0.219，0.548]，σ区间为[3.486，7.134]，95%的置信度区间为[0.432，1.081]，表明Qp3pl单元残差误差最大，Qhal-l单元残差误差最小，Qp3al及Qp3pl单元残差离散程度及残差置信度(95.0%)较大，Qhal+pl及Qp3al-l单元残差离散程度及残差置信度(95.0%)较小。Qhal-l、Qp3al-l及Qhal+pl单元预测效果要优于Qp3al及Qp3pl单元。
本文是基于地质单元、钻孔岩性及原位标贯进行预测分析得到的研究结果，由于同时给出了土体剪切波速随埋深、标贯击数的变化曲线，因此可由埋深数据反推获得密实度，对于无标贯数据的钻孔场地，可根据本文给出的土体密实度曲线进行相关预测研究。
5. 结论
(1) 对大同盆地进行典型地质单元划分，对大同盆地VS30与VS20及地质单元类型的相关性进行讨论，分析由VS20估算VS30的可行性。对盆地VS与地质特征的关系进行讨论，表明在考虑土体及地质单元类型前提下，标贯击数、密实度、剪切波速三者具有良好的正相关性，给出了基于剪切波速比的VS预测模型，给出了剪切波速与埋深的二次回归关系式。
(2) 本文研究表明，在区分地质单元类型的前提下，VS30与H呈明显负相关性。首次引入剪切波速比、密实度系数及第四系上部覆盖层厚度，给出VS30分解预测模型，模型拟合优度R2＞0.90，具有较高-极高的拟合优度。对于离散性较大、直接拟合估算较差的场地，本文提供了一种无剪切波速测试场地估算预测VS30的研究思路。
(3) 本文基于200m/s≤VS30≤400m/s统计样本进行研究，缺少VS30＜200m/s及VS30＞400m/s的优势统计样本，一方面反映了大同盆地VS30的优势分布区间特征，另一方面也与钻孔分布密度不均有关，后续将进一步结合相关测试数据对研究模型进行改进。
(4) 首次提出不同地质单元土体剪切波速比及密实度系数概念，首次综合地表地质及孔内原位试验等典型地质特征给出VS30预测模型，对于无任何剪切波速测试数据的场地，本预测模型综合考虑地表地质及孔内原位地质等地质要素，具有良好的推广使用价值。
-
图 1 大同盆地典型地质单元分区及钻孔分布图
Figure 1. Typical geological units and borehole distribution of Datong Basin
图 4 不同土体剪切波速-埋深关系图
Figure 4. Relationship between shear wave velocity and buried depth of different soils
图 5 洪积单元剪切波速与标贯击数关系图
Figure 5. Relationship between shear wave velocity and standard penetration number of diluvial units
图 6 冲积单元剪切波速与标贯击数关系图
Figure 6. Relationship between shear wave velocity and standard penetration number of alluvial units
图 7 冲湖积单元剪切波速与标贯击数关系图
Figure 7. Relationship between shear wave velocity and standard penetration number of alluvial-lacustrine units
图 9 土体VS30与第四系上部覆盖层关系
Figure 9. Relationship between VS30 and Quaternary upper overburden
图 10 大同盆地不同地质单元土体VS30预测结果
Figure 10. Prediction results of VS30 for different geological units of Datong Basin
表 1 大同盆地典型地质单元钻孔分布特征
Table 1. Borehole distribution characteristics of typical geological units of Datong Basin
地质单元类型 钻孔数量 VS30优势分布区间/(m·s-1) 钻孔分布特征 Qhal-l 292 224.5-272.6 分布于盆地中心Qh冲湖积平原及市区Ⅰ级阶地，钻孔分布面积面大，其中山阴、应县、怀仁及浑源县城钻孔密度较大 Qp3al-l 274 250.8-294.7 分布于市区河谷两岸，地貌上为Ⅱ级阶地，与Qp3al单元前缘相邻，除朔州市区东南部以外，盆地内该单元钻孔分布密度最大 Qhal+pl 221 268.4-306.9 分布于大同和朔州市区河床及山前倾斜平原前缘，地貌上与Qp3pl前缘及Qp3al后缘相邻，山阴、怀仁及浑源县城钻孔分布较密 Qp3al 423 274.6-332.4 分布于大同及朔州市区Ⅱ级以上高阶地，地貌上处于Qhpl及Qp3pl单元前缘，除市区钻孔分布稀疏外，其余地段钻孔密度均较大 Qp3pl 219 288.7-363.5 分布于山前倾斜洪积平原，地貌上与Qhpl及Qp3al单元后缘相邻，除市区钻孔分布较密外，其余地段钻孔密度较小 
下载: 导出CSV
表 2 不同土体剪切波速与埋深的相关性
Table 2. Correlation between shear wave velocity and buried depth of different soils
土体类型 VS=ah2+bh+c 拟合优度R2 a b c 卵砾石 -0.0299 6.7326 227.92 0.7591 中粗砂 -0.0275 6.4253 215.77 0.8755 粉质粘土 -0.0235 5.9766 216.12 0.8256 粉土 -0.0312 6.5989 197.99 0.8337 粉细砂 -0.0371 7.2017 178.39 0.9392
下载: 导出CSV
表 3 剪切波速回归分析及相同密实状态下的剪切波速比
Table 3. Regression analysis of shear wave velocity and shear wave velocity ratio under the same compaction state
地质单元 回归项 土体类型 卵砾石 中粗砂 粉细砂 粉质粘土 粉土 冲湖积单元al-l k1 110.42 109.40 105.23 104.02 108.86 k2 0.345 0.333 0.322 0.343 0.321 拟合优度R2 0.843 0.764 0.850 0.804 0.741 Sig(F显著性检验) ＜0.0001 波速比区间 1.048-1.106 1.021-1.049 0.968-0.970 0.986-1.035 1.000 波速比均值 1.087 1.039 0.969 1.018 1.000 冲积单元al k1 119.34 123.40 106.51 120.91 122.52 k2 0.346 0.320 0.344 0.326 0.315 拟合优度R2 0.815 0.780 0.825 0.742 0.716 Sig(F显著性检验) ＜0.0001 波速比区间 1.035-1.113 1.032-1.044 0.929-0.995 1.014-1.028 1.000 波速比均值 1.088 1.040 0.973 1.023 1.000 洪积单元pl k1 130.86 123.23 115.37 120.04 116.33 k2 0.348 0.352 0.353 0.356 0.358 拟合优度R2 0.770 0.755 0.762 0.764 0.772 Sig(F显著性检验) ＜0.0001 波速比区间 1.084-1.109 1.036-1.051 0.975-0.985 1.028-1.030 1.000 波速比均值 1.092 1.042 0.978 1.029 1.000
下载: 导出CSV
表 4 不同密实状态土体剪切波速比
Table 4. Soil shear wave velocity ratios in different compaction states
土体类型 土体密实状态 松散(流-软塑) 稍密 可塑 中密(硬塑) 密实(坚硬) 卵砾石 区间 0.64-0.78 0.80-0.922 0.88-0.98 1.00 1.21-1.28 均值 0.70 0.86 0.94 1.00 1.26 其它土体 区间 0.63-0.78 0.79-0.92 0.90-0.99 1.00 1.16-1.18 均值 0.69 0.86 0.94 1.00 1.16
下载: 导出CSV
表 5 基岩覆盖层厚度影响系数
Table 5. Influence coefficient of bedrock overburden thickness
基岩覆盖层厚度l/m > 60 50-60 40-50 30-40 < 30 λ(l) 1.00 1.01 1.03 1.05 1.08
下载: 导出CSV
表 6 S(H)回归分析
Table 6. Regression analysis of S(H)
编号 地质单元类型 统计样本 回归系数 拟合优度 F显著性检验 样本量 标准差 a b c R2 F Sig 1 Qhal-l 292 7.48 -0.018 0.554 259.47 0.724 378.29 ＜0.0001 2 Qp3al-l 274 8.14 -0.006 -0.264 291.38 0.758 409.55 3 Qhal+pl 221 8.64 -0.002 -0.607 308.78 0.812 423.49 4 Qp3al 423 11.18 -0.006 -0.530 331.05 0.732 579.50 5 Qp3pl 219 16.56 -0.007 -0.969 365.98 0.769 352.33
下载: 导出CSV
表 7 预测值VS30与实测值VSE30拟合优度及残差描述统计
Table 7. Fitting goodness and residual descriptive statistics of predicted value VS30 and measured value VSE30
编号 地质单元类型 拟合优度R2 均方根误差
RMSE残差δ 平均值M 标准误差Se 标准差σ 置信度(95.0%) 1 Qhal-l 0.828 4.584 0.031 0.269 4.592 0.529 2 Qp3al-l 0.843 3.660 -0.035 0.219 3.486 0.432 3 Qhal+pl 0.904 3.523 -0.028 0.237 3.531 0.468 4 Qp3al 0.824 5.834 0.155 0.282 5.843 0.555 5 Qp3pl 0.865 7.100 -0.047 0.548 7.134 1.081 6 最终预测结果 0.973 5.376 0.010 0.142 5.379 0.278
下载: 导出CSV
-
陈鲲, 俞言祥, 高孟潭, 2010.考虑场地效应的ShakeMap系统研究.中国地震, 26(01):92-102. doi: 10.3969/j.issn.1001-4683.2010.01.009 《工程地质手册》编委会, 2018.工程地质手册(第五版).北京:中国建筑工业出版社. 黄雅虹, 吕悦军, 彭艳菊, 2009.国内外不同抗震设计规范中场地分类方法的内在关系研究.震灾防御技术, 4(1):80-90. doi: 10.3969/j.issn.1673-5722.2009.01.008 李铁锋, 任明达, 1993.大同盆地晚新生代环境演化特征.北京大学学报(自然科学版), 29(4):476-483. doi: 10.3321/j.issn:0479-8023.1993.04.001 江志杰, 彭艳菊, 方怡等, 2018.北京平原地区VS30估算模型适用性研究.震灾防御技术, 13(01):75-86. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180107&journal_id=zzfyjs 吕悦军, 彭艳菊, 兰景岩等, 2008.场地条件对地震动参数影响的关键问题.震灾防御技术, 3(02):126-135. doi: 10.3969/j.issn.1673-5722.2008.02.003 彭艳菊, 吕悦军, 黄雅虹等, 2009.工程地震中的场地分类方法及适用性评述.地震地质, 31(02):349-362. doi: 10.3969/j.issn.0253-4967.2009.02.016 沈方铝, 李培, 张颖等, 2018.福州市区土体剪切波速与土层深度的经验关系研究.地震工程学报, 40(S1):83-89. 史大成, 温瑞智, 杜春清, 2012.区域性场地VS30及峰值加速度放大系数估算方法.地震工程与工程振动, 32(04):40-46. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzgcygczd201204007 汪云龙, 袁晓铭, 轩浩等, 2016.级配砂石剪切波速与相对密度关系实验研究.自然灾害学报, 25(02):173-178. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zrzhxb201602021 王乃梁, 杨景春, 夏正锴等, 1996.山西地堑系新生代沉积与构造地貌.北京:科学出版社, 156-200. 王琦, 刘红帅, 郑桐等, 2018.天津地区覆盖土体剪切波速与埋深的相关性分析.地震工程与工程振动, 38(06):190-201. 喻畑, 李小军, 2015.四川、甘肃地区VS30经验估计研究.地震工程学报, 37(02):525-533. doi: 10.3969/j.issn.1000-0844.2015.02.0525 瞿伟, 王庆良, 张勤等, 2013.大同盆地现今地壳形变及应变分布特征.大地测量与地球动力学, 33(03):11-15. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201303003 张世民, 窦素芹, 杨景春, 1997.大同第四纪火山群的活动特点.地壳构造与地壳应力文集, 1997(00):103-111. http://www.cnki.com.cn/Article/CJFDTotal-SEIS199700013.htm 中华人民共和国建设部, 中华人民共和国国家质量监督检验检疫总局, 2009.GB 50021-2001岩土工程勘察规范(2009).北京: 中国建筑工业出版社. 中华人民共和国住房和城乡建设部, 2016.GB 50011-2010建筑抗震设计规范(2016).北京: 中国建筑工业出版社. Abrahamson N., Silva W., 2008. Summary of the Abrahamson & Silva NGA ground-motion relations. Earthquake Spectra, 24(1):67-98. doi: 10.1193/1.2924360 Xie J. J., Zimmaro P., Li X. J., et al., 2016. VS30 empirical prediction relationships based on a new soil-profile database for the Beijing plain area, China. Bulletin of the Seismological Society of America, 106(6):2843-2854. doi: 10.1785/0120160053 Wald D. J., Earle P. E., Lin K., et al., 2006. Challenges in rapid ground motion estimation for the prompt assessment of global urban earthquakes. Bull. Earthq. Res. Inst. Univ. Tokyo, 81(1):273-281. Wills C. J., Gutierrez C. I., Perez F. G., et al., 2015. A next generation VS30 map for California based on geology and topography.Bull. Seism. Soc. Am, 105(6):3083-3091. doi: 10.1785/0120150105 -
下载:
点击查看大图
计量
- 文章访问数: 204
- HTML全文浏览量: 78
- PDF下载量: 4
- 被引次数: 0
