Selection of VS30 Empirical Estimation Model for Sites Consideration of Geological Condition Dependent−A Case of Xinjiang Area
-
摘要: 新疆地区地震频发且强度较大，震害较为严重，因此对该区进行地震危险性分析很有必要。场地效应是影响危险性分析结果的重要因素之一，不同场地条件对地震动的放大作用不同；常用的表征场地条件的重要参数有VS30，但在实际工程中通常由于经济等原因，钻孔深度无法达到30m，不能直接获取真实的VS30值，因此选用合理的VS30经验估计模型成为研究人员关注的热点问题。为了得到更准确的VS30估计值，本文共收集了新疆地区1307个工程钻孔数据，针对目前常用的5种经验模型，基于地质图信息将钻孔场地划分为6类年代单元，采用相关系数检验法分析回归模型的系数，采用对数正态分布检验方法分析钻孔的分布特征，选取各类单元最适用的模型。结果表明，相关系数检验法在最适用模型评价方面是失效的，而通过比较对数正态分布检验下的峰度系数和偏度系数可选取出各类地质单元最适用模型。最后根据模型选取结果，计算得到全部工程钻孔的VS30值，并利用克里金插值方法得出新疆地区VS30值的分布结果，研究结果可为新疆地区VS30的潜在使用者提供有价值的辅助信息，此分析方法也可为我国其他地区VS30经验估计模型的选用提供技术思路。Abstract: Earthquakes occurs in Xinjiang area frequently with high intensity and serious damage, that calls for seismic hazard analysis in this area certainly. The site effect is an important factor affecting the result of seismic hazard analysis. Different site conditions give different contributions on site amplification on earthquake ground motions. As we have known, VS30 is one of the important parameters used commonly to characterize site condition. In some cases, however, the drilling depth can not reach 30 m due to limited funds or other reasons so that the true value of VS30 cannot be obtained directly, which traises a popular research issue that how to choose a reasonable VS30 empirical estimation model. In this regard, a total of 1307 engineering boreholes in Xinjiang area have been collected. The boreholes are classified into six types based on their site surface geological age extracted from geological map. For five empirical models used commonly at present, both the correlation coefficient testing and logarithmic distribution testing methods are used to evaluate which model is the most suitable one for each type of geological unit. The results show that the correlation coefficient testing method is ineffective, but the kurtosis and skewness coefficients under the lognormal distribution testing method is capable of selecting the most suitable model for each geological unit. Finally, the VS30 values of all engineering boreholes are calculated with selected empirical models, and the spatial distribution of VS30 values in Xinjiang are obtained by Kriging interpolation, which can provide valuable site information for potential users of VS30 in this area. The method proposed in this study can also provide the technical base for selection of VS30 empirical estimation models in other regions of China.
-
Key words:
- Earthquake hazard /
- Geological age /
- VS30 /
- Correlation coefficient /
- Lognormal test
-
引言
地震动衰减关系是表征地震动参数随震级、距离和场地等因素变化规律的函数。地震动衰减关系作为工程地震学中最常用的估计或预测地震动的方法，被广泛用于工程地震学领域，包括地震危险性分析、地震区划、地震小区划和设计地震动确定等多方面(袁一凡等，2012)。地震动衰减关系通常采用相对固定的函数形式，大致由震源影响函数、距离影响函数、场地影响函数和修正项4个部分构成。衰减关系的系数值一般利用大量强震动记录资料统计确定。每个地区的地震构造环境、地壳介质以及地表和近地表局部场地条件等因素具有很大的差异性，因此，地震动衰减关系具有很强的地区性。欧洲(Akkar等，2014)、美国(Sadigh等，1997)、中国台湾(Lin等，2011)等国家和地区都获得基于强震动记录的适用于相应地区的地震动衰减关系。
何种地震动衰减关系更为适合1个特定的地区始终是研究的热点。Campbell(2016)利用NGA-West2中的数据，研究了3个针对欧洲和中东地区的地区性衰减关系和1个NGA衰减关系，分别比较了地震动衰减关系的预测中值和标准偏差，并进行了残差分析。认为Campbell-Bozorgnia衰减关系在更具有工程实用价值的大震和近场上优于另外3个衰减关系。Zafarani等(2017)为了找到更适合伊朗高原的小震和中震的地震动预测方程(GMPE)，使用3种不同的统计学方法，对8个GMPE进行了比较，其中只有Zafarani于2015年提出的衰减关系在全部的频率范围内均与数据较吻合。
虽然中国有很长的历史地震记载，但已获取强震动记录的地震事件不多，总体的强震动记录资料较少。因此，长久以来中国研究地震影响的分布往往采用烈度衰减关系。然而，地震烈度的空间分布虽然也能体现地震动强度的空间变化，反映其衰减特征，但它并非能直接反映地震动特性的相关物理量，甚至可能受到一些主观因素的干扰。经过多年来的发展，目前中国已经积累了一定量的强震动加速度记录，这或将使得我们有可能直接建立地震动加速度衰减关系(王继等，2008)。此外，也可以使用中国的强震动加速度记录，与目前国际上已有的衰减关系进行比较，找出其中较为适用的衰减关系，并进而研究对已有衰减关系进行修正处理的方法。司宏俊等(2014)比较了汶川地震和芦山地震中基岩和土层场地观测台站得到的加速度反应谱和基于日本强震动观测数据得出的基岩加速度反应谱衰减关系，发现两者有着相似的衰减特征，并通过拟合强震动记录加速度反应谱及其与衰减关系的残差，分析得到汶川地震和芦山地震的地震动记录台站的场地影响，与日本衰减关系中平均场地参数比较接近，VS30约为500m/s。Dangkua等(2018)考察了NGA-West2项目给出的地震动衰减关系对于中国地区的适用性，认为NGA-West2中的衰减关系对短周期地震动适用性较好，但在长周期范围结果偏高；此外，还利用对数最大似然法，为这些衰减关系分配了权重。
在地震动衰减关系中，峰值加速度(PGA)是主要被考虑的参数。地震动峰值加速度是指在一次特定地震事件中取得的1条加速度记录的绝对最大值，常被用于描述地震动的强度以及对结构反应的影响程度。虽然峰值加速度并非能完全控制结构反应的地震动参数，特别是对于较长自振周期的结构，但因为其是强震动记录中最容易取得的参数，且结构抗震设计规范中多采用标准形式的地震动反应谱与峰值加速度结合使用，因此，人们在进行衰减关系的研究时仍然十分关心且首先关注峰值加速度(Douglas，2003)。
1. 数据来源与分析
1.1 强震动记录
在地震动衰减关系研究中，为了关注有显著工程意义的地震动记录，同时也为了减小独立变量的误差，可以选取MW＞5地震的地震动记录进行研究(Douglas，2003)。因此，本文采用中国西部地区2012年3月9日—2017年9月16日间发生的42次MW＞5地震事件中共916条有效水平方向强震动记录(峰值加速度绝对值超过10gal)。震中位置和发震时间来自地震简报，震源机制和震级参考GCMT(Dziewonski等，1981；Ekström等，2012)，详细情况如表 1所示。研究中涉及的地震事件和地震台站的位置分别见图 2、图 3。
表 1 地震事件基本信息Table 1. Basic information of earthquake events事件编号 地震时间 纬度/°N 经度/°E 矩震级MW 1 2012年3月9日6:50 39.4 81.3 5.9 2 2012年9月7日11:19 27.5 104 5.6 3 2013年1月18日20:42 31 99.4 5.6 4 2013年3月3日13:41 25.9 99.7 5.4 5 2013年3月11日11:01 40.2 77.5 5.2 6 2013年3月29日13:01 43.4 86.8 5.4 7 2013年4月17日9:45 26.9 99.8 5.3 8 2013年4月20日8:02 30.3 103 6.6 9 2013年7月22日7:45 34.5 104.2 6.0 10 2013年8月30日13:27 43.8 87.6 5.3 11 2013年9月20日5:37 37.7 101.5 5.1 12 2014年5月24日4:49 25 97.8 5.7 13 2014年5月30日9:20 25 97.8 5.9 14 2014年7月9日5:52 39.3 78.3 5.0 15 2014年8月3日16:30 27.1 103.3 6.2 16 2014年8月17日6:07 28.1 103.5 5.1 17 2014年10月2日23:56 36.4 97.8 5.2 18 2014年10月7日21:49 23.4 100.5 6.1 19 2014年11月22日16:55 30.26 101.69 6.1 20 2014年11月25日23:19 30.18 101.73 5.7 21 2014年12月6日2:43 23.31 100.5 5.6 22 2014年12月6日18:20 23.31 100.5 5.5 23 2015年1月10日13:21 40.2 77.3 5.1 24 2015年3月1日18:24 23.5 98.9 5.3 25 2015年4月15日15:39 39.8 106.3 5.5 26 2015年7月3日9:07 37.6 78.2 6.4 27 2015年10月26日17:09 36.5 70.8 7.5 28 2015年11月23日5:02 38 100.4 5.2 29 2015年12月7日15:50 38.2 72.9 7.2 30 2016年1月14日5:18 42.19 84.12 5.1 31 2016年1月21日1:13 37.71 101.7 5.9 32 2016年5月18日0:48 26.1 99.53 5.0 33 2016年6月26日16:17 39.43 73.4 6.4 34 2016年9月23日1:23 30.1 99.61 5.2 35 2016年11月25日22:24 39.27 74.04 6.6 36 2016年12月8日13:15 43.83 86.35 6.0 37 2017年3月27日7:55 25.89 99.8 5.1 38 2017年6月3日18:11 37.99 103.56 5 39 2017年8月8日21:19 33.2 103.82 6.5 40 2017年8月9日7:27 44.27 82.89 6.5 41 2017年9月16日18:11 42.11 83.43 5.5 42 2017年9月30日14:14 32.27 105 5.1 1.2 观测台站场地VS30值的确定
VS30即场地地下30m深度内岩土层的平均剪切波速，是用于表征场地条件的重要物理参数。研究中涉及的强震动观测台站，其VS30值需要借助于已有的研究成果确定。因此，本文多数观测台站的VS30值采用了姜治军(2017)和俞畑等(2015)的研究成果，而对于已有研究成果中没有涉及的台站，其VS30值借助100km范围内最近台站的VS30或同类场地的VS30均值进行近似处理来确定。此外，本次研究还利用观测场地的VS30值对观测台站的场地类型进行划分。传统上将VS30=750m/s的场地作为基岩场地和土层场地的分界波速值，但因为VS30＞750m/s的台站十分缺乏，为此，本文以VS30=360m/s将场地划分两类，一类包括基岩和硬土场地，另一类包括中硬土和软土场地。这一场地分类处理与Sadigh等(1997)的做法一致。
文中选择的强震动记录大部分来自基岩和硬土场地，仅有少部分数据来自中硬土和软土场地(图 1)。
图 1 强震动记录的震级-震中距分布Figure 1. The magnitude v.s. epic-center distance of strong ground motion records used in this study2. 地震动峰值加速度衰减关系的比较
2.1 地震动衰减关系的选择
选择来自欧洲(Akkar等，2014)、美国(Sadigh等，1997)和中国台湾(Lin等，2011)的3个较有代表性的地区性地震动衰减关系。这些衰减关系的建立均未利用中国西部地区的强震动记录数据，在确定衰减关系系数时，是以当地的强震动记录数据为主，并加入其它地区发生的大地震的强震动记录进行补充。选择的衰减关系均采用浅层地壳地震的强震动记录进行推导，与本文采用的中国西部地区强震动记录的地震具有可比性。所选地震动衰减关系的具体情况见表 2。
表 2 所用地震动峰值加速度衰减关系及相关参数Table 2. Peak acceleration attenuation relations and their parameters adopted in this study衰减关系编号 地区 震级MW 震中距/km 参数 ASB14(Akkar等，2014) 泛欧和中东 4.0—7.6 ＜200 RJB，Repi，Rhyp，MW，VS30，FN，FR SCEMY97(Sadigh等，1997) 美国西部 3.8—7.4 ＜200 Rrup(Repi)，MW，FN，FS LLCS11(Lin等，2011) 中国台湾 3.5—7.6 ＜240 Rrup(Repi)，MW 表 2中，RJB为Joyner-Boore距离，Rhyp为震源距，Repi为震中距，Rrup为台站到破裂平面的距离(震级不太大时可近似为Repi)，MW为矩震级，VS30为30m处剪切波速，FN、FR和FS分别为正断层、逆冲断层和走滑断层的指标，一般为0；当且仅当震源机制为正断层时，FN=1；当且仅当震源机制为逆冲时，FR=1；当且仅当震源机制为走滑时，FS=1。
针对本文所用的地震动峰值加速度衰减关系的震级与震中距的适用范围，进一步将所采用的地震动记录限制在震级MW 5.0—7.5、震中距小于200km范围内，即只利用了震中距小于200km范围内的地震动记录，而不考虑震中距更大的那些地震动记录。
衰减关系LLCS11还考虑了断层上、下盘效应对地震动的影响，并分别为上、下盘给出了不同的系数。考虑到当地的地震构造背景，这种划分也许是有益的。但是，因为本文获得的地震动数据难以区分断层上、下盘，故本文将LLCS11中同种场地的上、下盘取平均使用。
本文采用的3个地震动衰减关系的峰值加速度衰减计算公式如下(峰值加速度单位为g)。
2.1.1 衰减关系ASB14
$$ \text{ln}(Y)=\text{ln}[{{Y}_{\text{REF}}}({{M}_{w}}, R, SoF)]+\text{ln}[S({{V}_{S30}}\text{, PG}{{\text{A}}_{\text{REF}}})]+\varepsilon \sigma $$ (1) $$ {\rm{ln}}({Y_{{\rm{REF}}}}) = \\ \left\{ \begin{array}{l} {a_1} + {a_2}({M_{\rm{W}}} - {c_1}) + {a_3}{(8.5 - {M_{\rm{W}}})^2} + [{a_4} + {a_5}({M_{\rm{W}}} - {c_1})]\ln (\sqrt {{R^2} + a_6^2}) + {a_8}{F_N} + {a_9}{F_R}, {M_{\rm{W}}} \le {c_{\rm{1}}}\\ {a_1} + {a_7}({M_{\rm{W}}} - {c_1}) + {a_3}{(8.5 - {M_{\rm{W}}})^2} + [{a_4} + {a_5}({M_{\rm{W}}} - {c_1})]\ln (\sqrt {{R^2} + a_6^2}) + {a_8}{F_N} + {a_9}{F_R}, {M_{\rm{W}}}{c_{\rm{1}}} \end{array} \right. $$ (2) $$ \ln \left( S \right) = \left\{ \begin{array}{l} {b_1}\ln \left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right) + {b_2}\ln \left[ {\frac{{{\rm{PG}}{{\rm{A}}_{{\rm{REF}}}} + c{{\left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right)}^n}}}{{\left( {{\rm{PG}}{{\rm{A}}_{{\rm{REF}}}} + c} \right){{\left( {{V_{{\rm{S}}30}}/{V_{{\rm{REF}}}}} \right)}^n}}}} \right],{V_{{\rm{S30}}}} \le {V_{{\rm{REF}}}}\\ {b_1}{\rm{ln}}\left[ {\frac{{{\rm{min}}\left( {{V_{{\rm{S}}30}},{V_{{\rm{CON}}}}} \right)}}{{{V_{{\rm{REF}}}}}}} \right],{V_{{\rm{S}}30}} > {V_{{\rm{REF}}}} \end{array} \right. $$ (3) 其中，PGAREF和VS30REF为参考场地的峰值地震动和VS30，参考场地的VS30为750m/s。YREF表示由衰减关系计算的参考场地的地震动变量，当计算参考场地的峰值地震动时，YREF=PGAREF。各系数取值为a1=2.52977，a2=0.0029，a3=-0.05496，a4=-1.31001，a5=0.2529，a6=7.5，a7=-0.5096，a8=-0.1091，a9=0.937，b1=-0.41997，b2=-0.28846，c=2.5，c1=6.75，n=3.2。
2.1.2 衰减关系SCEMY97
对于基岩场地：
$$ {\rm{ln}}(y) = {C_1} + {C_2}M + {C_3}{(8.5M)^{2.5}} + {C_4}\ln ({R_{rup}} + \exp ({C_5} + {C_6}M)) + {C_7}\ln ({R_{rup}} + 2) $$ (4) 其中，MW≤6.5：C1=-0.524，C2=1.0，C3=0.000，C4=-2.100，C5=1.29649，C6=0.250，C7=0.0；MW＞6.5：C1=-1.274，C2=1.1，C3=0.000，C4=-2.100，C5=-0.48451，C6=0.524，C7=0.0。
对于土层场地：
$$ {\rm{ln}}(y) = {C_1} + {C_2}M - {C_3}\ln ({R_{rup}} + {C_4}{e^{{C_5}M}}) + {C_6} + {C_7}{(8.5 - M)^{2.5}} $$ (5) 其中，对于走滑型地震，C1=-2.17；对于正断层地震和逆冲型地震，C1=-1.92；当MW≤6.5，C4=2.1861，C5=0.32；当MW＞6.5，C4=0.3825，C5=0.5882；C2=1.0，C3=1.70，C6=0.0，C7=0.0。
2.1.3 衰减关系LLCS11
对于基岩场地：
$$ \left\{ \begin{array}{l} {\rm{ln}}({\rm{PGA}}) = - 3.279 + 1.035M - 1.651\ln (R + 0.152{e^{0.623M}}), 上盘\\ {\rm{ln}}({\rm{PGA}}) = - 3.232 + 1.047M - 1.662\ln (R + 0.192{e^{0.630M}}), 下盘 \end{array} \right. $$ (6) 对于土层场地：
$$ \left\{ \begin{array}{l} {\rm{ln}}({\rm{PGA}}) = - 3.248 + 0.943M - 1.471\ln (R + 0.100{e^{0.648M}}), 上盘\\ {\rm{ln}}({\rm{PGA}}) = - 3.218 + 0.935M - 1.464\ln (R + 0.125{e^{0.650M}}), 下盘 \end{array} \right. $$ (7) 2.2 峰值加速度衰减关系计算值与实际记录比较
分别根据上述3个衰减关系公式计算峰值加速度，并与实际强震动记录中的峰值加速度进行比较，计算残差(R)。同时计算了残差在不同震中距和震级上的标准偏差(σ)，得到了残差随距离与震级的分布图像(图 4、图 5和图 6)。本文残差为实际记录中的峰值加速度减去衰减关系计算值，标准偏差为方差的平方根。
图 4 ASB14计算峰值加速度与实际强震动记录的残差随震级和距离的分布Figure 4. The distribution of residuals between peak ground acceleration calculated using ASB14 and recorded against magnitude and distance
图 5 SCEMY97计算峰值加速度与实际强震动记录的残差随距离和震级的分布Figure 5. The distribution of residuals between peak ground acceleration calculated using SCEMY97 and recorded against magnitude and distance
图 6 LLCS11计算峰值加速度与实际强震动记录的残差随距离和震级的分布Figure 6. The distribution of residuals between peak ground acceleration calculated using LLCS11 and recorded against magnitude and distance由图 4、图 5和图 6可以看出，大多数数据点的残差虽然不为0，但均位于[R-σ，R+σ]范围内，且残差的标准偏差与震级、距离关系不大，说明可能存在与距离和震级无关的误差。
SCEMY97和LLCS11整体上存在正的残差，即根据这2个衰减关系计算得到的峰值加速度小于中国西部地区的强震动记录值。这种趋势随距离和震级基本不变，可能意味着2个衰减关系需要引入1个常数项进行修正。ASB14虽然也有整体为正的残差，但残差随距离增大而增大，说明该衰减关系需要引入1个随距离增大而增大的项进行修正。总体来说，ASB14的残差小于SCEMY97和LLCS11。
此外，分别计算了ASB14、SCEMY97和LLCS11的残差与震级、距离的相关系数(表 3)，发现SCEMY97和LLCS11的残差与距离成中等相关，两者之间可能有中等程度的线性关系。这意味着若要对这些衰减关系进行修改，可以尝试引入线性项。考虑到相关系数的特点，这种线性关系可能需要更多的观测记录才能确定。
表 3 残差与震级、距离的相关系数Table 3. Residual-magnitude and residual-distance correlation coefficients衰减关系编号 残差与震级的相关系数 残差与距离的相关系数 ASB14 -0.121235 0.495025 SCEMY97 0.138204 0.64689 LLCS11 0.0942092 0.582699 3. 结论
观察上述3个地震动衰减(ASB14、SCEMY97和LLCS11)的残差与震级的关系，可以发现在地震动峰值加速度衰减关系中，考虑震级的平方项有重要的意义。衰减关系中震级的平方项与高频地震动的震级饱和现象有关，而这一现象在震级较大时表现更为明显。因为实际工程中更应关注的是震级较大的地震，所以在建立地震动峰值加速度衰减关系时，应该尽量采用具有震级的平方项的形式。
此外，上述3个地震动衰减关系的残差关于距离均存在一定程度的偏离。观察残差关于距离的分布，发现若对ASB14进行常数项的修正，对SCEMY97和LLCS11进行距离的线性项和常数项的修正，可能会使3个地震动衰减关系更适合于中国西部地区。
对于中国西部强震动记录数据，ASB14的表现优于SCEMY97和LLCS11，如果直接利用国外的地震动衰减关系，ASB14更适合于中国西部地区。然而总体来看，无论是随着震级还是距离的变化，采用3个衰减关系预测中国西部地区的地震动将会明显低估地震动峰值加速度。因此，发展强震动观测以积累更多的强震动记录，并基于中国不同区域的强震动记录资料直接统计确定相应的地震动衰减关系，仍然是需要努力的目标。
-
图 1 不同外推模型计算所得VS30估计值与实测值对比
Figure 1. Comparison of VS30 estimated values to measured values calculated by different extrapolation models
图 2 新疆地区数据所得拟合系数与原模型拟合系数对比
Figure 2. The comparison of fitting coefficients between the data in Xinjiang region to that of the original model
图 3 新疆地区地质单元及钻孔空间分布图
Figure 3. Geological unit map and location of the drilling sites in Xinjiang area
图 4 新疆地区工程钻孔在不同地质年代单元下钻孔数量统计结果
Figure 4. The statistical results of the number of drilling holes drilled under different geological age units in Xinjiang area
图 5 不同地质年代分类下钻孔VS30实测值与估计值之间的相关系数随深度变化曲线
Figure 5. Correlation coefficient with depth between the value of VS30 measured from drilled hole and the estimated VSZ calculated by using five extrapolated models under different geological age classification
图 6 新疆地区不同地质单元下工程钻孔实测VS30数据的对数正态检验峰度及偏度系数结果
Figure 6. Lognormal test kurt and skewness coefficient of engineering bore VS30 data measured under different geological units in Xinjiang
图 7 不同地质年代下各外推模型峰度及偏度系数检验结果图
Figure 7. Test results of the kurtosis and skewness coefficient of each extrapolated model under different geological ages
-
陈兰生, 冀昆, 温瑞智等, 2021. 考虑场地放大效应的全国地震危险分析. 地震工程与工程振动, 41（6）: 177—185Chen L. S. , Ji K. , Wen R. Z. , et al. , 2021. Seismic hazard analysis for China considering site amplification effect. Earthquake Engineering and Engineering Dynamics, 41(6): 177—185. (in Chinese) 陈应, 2020. 基于地震动参数及网格化工程场地的地震危险性分析——以成都市为例. 南昌: 东华理工大学.Chen Y. , 2020. Seismic hazard analysis based on ground motion parameters and gridding engineering site-A case study of Chengdu. Nanchang: East China University of Technology. (in Chinese) 华永超, 2020. 基于场地周期的土体场地放大效应研究. 哈尔滨: 中国地震局工程力学研究所.Hua Y. C. , 2020. Study on soil site amplification effect based on site period. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese) 李兵, 谢富仁, 黄金水等, 2022. 龙门山断裂带大邑地震空区地应力状态与地震危险性. 中国科学: 地球科学, 52（7）: 1409—1418.Li B. , Xie F. R. , Huang J. S. , et al. , 2022. In situ stress state and seismic hazard in the Dayi seismic gap of the Longmenshan thrust belt. Science China Earth Sciences, 65(7): 1388—1398. 米欣雪, 2022. 一种场地VS30经验估计最佳模型的选用方法——以新疆地区为例. 见: 第十一届全国地震工程学术会议.Mi X. X. , 2022. A method for selecting the optimal model of VS30 empirical estimation: case study in Xinjiang area. In: Proceedings of the 11 th National Conference on Earthquake Engineering. (in Chinese) 缪铨生, 1997. 概率与数理统计. 上海: 华东师范大学出版社. 田禹, 2012. 基于偏度和峰度的正态性检验. 上海: 上海交通大学.Tian Y. , 2012. Tests for normality based on skewness and kurtosis. Shanghai: Shanghai Jiao Tong University. (in Chinese) 尹迪, 董培育, 曹建玲等, 2022. 川滇地区地震危险性数值分析. 地球物理学报, 65（5）: 1612—1627 doi: 10.6038/cjg2022P0651Yin D. , Dong P. Y. , Cao J. L. , et al. , 2022. Numerical analysis of the seismic hazard in Sichuan-Yunnan region. Chinese Journal of Geophysics, 65(5): 1612—1627. (in Chinese) doi: 10.6038/cjg2022P0651 张雨婷, 2020. 基于我国地形与地质条件的场地分类方法研究. 哈尔滨: 中国地震局工程力学研究所.Zhang Y. T. , 2020. Site classification method based on topographical and geological characteristics in China. Harbin: Institute of Engineering Mechanics, China Earthquake Administration. (in Chinese) 翠川三郎, 野木淑裕, 2015. 深さ30 mまでの地盤の平均S波速度を深さの浅いデータから推定する方法について. 日本地震工学会論文集, 15（2）: 91—96.Midorikawa S., Nogi Y., 2015. Estimation of VS30 from shallow velocity profile. Journal of Japan Association for Earthquake Engineering, 15(2): 91—96. (in Japanese) Boore D. M. , 2004. Estimating V̄s(30) (or NEHRP site classes) from shallow velocity models (depths<30 m). Bulletin of the Seismological Society of America, 94(2): 591—597. doi: 10.1785/0120030105 Boore D. M. , Thompson E. M. , Cadet H. , 2011. Regional correlations of VS30 and velocities averaged over depths less than and greater than 30 meters. Bulletin of the Seismological Society of America, 101(6): 3046—3059. doi: 10.1785/0120110071 Dai Z. J. , Li X. J. , Hou C. L. , 2013. A shear-wave velocity model for VS30 estimation based on a conditional independence property. Bulletin of the Seismological Society of America, 103(6): 3354—3361. doi: 10.1785/0120130025 Kwak D. Y., Ancheta T. D., Mitra D., et al., 2017. Performance evaluation of VSZ-to-VS30 correlation methods using global VS profile database. In: Proceedings of the 3 rd International Conference on Performance-Based Design in Earthquake Geotechnical Engineering. Vancouver, 399. Thompson E. M., Wald D. J., 2012. Developing VS30 site-condition maps by combining observations with geologic and topographic constraints. In: Proceedings of the 15 th World Conference on Earthquake Engineering. Lisbon, Portugal: IEM. -
下载:
下载:
点击查看大图
计量
- 文章访问数: 91
- HTML全文浏览量: 62
- PDF下载量: 18
- 被引次数: 0
