基于HVSR谱比动态聚类的海域场地特性研究

刘献伟; 陈苏; 李小军; 傅磊; 胡进军; 孙浩

doi:10.11779/CJGE20211449

基于HVSR谱比动态聚类的海域场地特性研究 English Version

刘献伟^1,,
陈苏^2, ,,
李小军²,
傅磊¹,
胡进军³,
孙浩^{4, 5}

1.
中国地震局地球物理研究所, 北京 100081
2.
北京工业大学城市建设学部, 北京 100124
3.
中国地震局工程力学研究所, 黑龙江哈尔滨 150080
4.
中国人民大学高瓴人工智能学院, 北京 100872
5.
东北大学, 波士顿 02115

基金项目:

国家自然科学基金项目 51878626

国家自然科学基金项目 U1839202

中央级公益性科研院所指令性项目 DQJB20R13

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

详细信息

作者简介:
刘献伟(1997—)，男，内蒙古赤峰人，博士研究生，主要从事海域场地特性分析、场地分类及地震动衰减关系等研究工作。E-mail: liuxianwei97@163.com

通讯作者:
陈苏, E-mail: chensuchina@126.com

中图分类号: TU435
计量
- 文章访问数: 202
- HTML全文浏览量: 65
- PDF下载量: 72
出版历程
- 收稿日期: 2021-12-16
- 网络出版日期: 2023-02-03
- 刊出日期: 2022-12-31

Characteristics of marine site based on HVSR dynamic clustering method

LIU Xianwei^1,,
CHEN Su^2, ,,
LI Xiaojun²,
FU Lei¹,
HU Jinjun³,
SUN Hao^{4, 5}

1.
Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
2.
College of Civil Engineering, Beijing University of Technology, Beijing 100124, China
3.
Institute of Engineering Mechanics, China Earthquake Administration, Harbin 150080, China
4.
Renmin University of China, Beijing 100872, China
5.
Northeastern University, Boston 02115, America

摘要

摘要: 海域场地工程特性是海洋重大工程设计、建造、运维的重要基础，也是岩土地震工程的研究热点。收集了日本相模湾海域2000~2018年间6个海底台站的强地震动记录，筛选、前处理后采用Horizontal-to-vertical spectral ratio (HVSR)方法计算了各条记录的傅里叶振幅谱比值。通过动态聚类学习结合最优化分类，获取了相模湾海底台站的最优场地特性。结果表明，动态聚类学习方法针对时序型数据的聚类性能及鲁棒性强，可满足局部及全局分辨率；相模湾海区中台站场地环境可宏观分为海底峡谷区及海洋盆地区，其中峡谷区场地卓越周期显著小于盆地区，且高频衰减较缓；通过局部场地地形坡度，给出了海洋场地工程特性。研究成果可为海域地震动模拟、海洋工程建设等提供参考。
- 海域地震动 /
- HVSR方法 /
- 动态聚类 /
- 海域场地 /
- 场地分类
Abstract: The engineering characteristics of marine site are the important basis for the design, construction, operation and maintenance of major marine projects and the research focus of geotechnical earthquake engineering. The strong ground motion records of 6 seabed stations in Sagami Bay, Japan from 2000 to 2018 are collected, after screening and pre-processing, the horizontal-to-vertical spectral ratio (HVSR) method is used to calculate the Fourier amplitude spectrum ratio of each record. The optimal site characteristics of Sagami Bay are obtained through the dynamic clustering learning combined with the optimal classification. The results show that the dynamic clustering learning method has strong performance and robustness for time series data, and can meet the local and global resolutions. The site of Sagami Bay can be macroscopically divided into submarine canyon area and basin area. The dominant period of the canyon area is significantly smaller than that of the basin area, and the high frequency attenuation is relatively slow. The engineering characteristics of marine site are given by the topographic slope of the local site. The relevant research results may provide reference for simulation of sea ground motion and construction of marine projects.
- sea ground motion /
- horizontal-to-vertical spectral ratio method /
- dynamic clustering /
- marine site /
- site classification

HTML全文

0. 引言

近年来，海洋开发和利用不断加深，特别是中国海上风电、跨海大桥、海上石油平台、海底隧道、海底管道等重大海洋工程投入使用。与此同时，中国海域及邻区所处地震地质环境复杂，地震活动性较强，海域地震频发^[1]，海洋工程的抗震设防能力及重要性凸显，海域场地工程特性得到充分关注^[2]。针对中国海域及邻近区域的场地地震反应研究工作仍不足^[1]，对海域土层条件，场地特性及分类等研究还处在起步阶段^[3]，所以研究海域场地效应和海域场地分类可以为海洋工程抗震设防，海洋结构安全性评价提供帮助。

Nakamura^[4]提出基于地脉动的水平和竖向谱比（horizontal-to-vertical spectral ratio, HVSR）作为场地特征评估方法，国内外涌现很多基于HVSR的场地地震效应和场地分类的研究。其中Boore等^[5]认为水平向地震动与竖向地震动之比（H/V）可以消除震源、断层类型、传播路径等因素，用以分析场地效应；同时Boore等^[5]认为H/V谱比法可同时应用于陆域和海域场地的研究。Field等^[6]利用HVSR法研究了1994年北岭地震场地非线性放大效应。Bonilla等^[7]利用南加州地区的井下阵列钻孔记录对比了HVSR法和传统的场地特征研究方法，发现两种方法得到的场地卓越周期基本一致。Wen等^[8]利用HVSR法分析了集集地震的场地非线性反应。任叶飞等^[9]将广义线性反演法和HVSR法进行对比，得出HVSR方法可以应用于评估场地卓越周期。荣棉水等^[10]利用美国Garner Valley Downhole Array（GVDA）竖向观测台阵的记录，对比了传递函数法和HVSR法的差异，并指出在可忽略场地竖向放大的情况下，可以利用HVSR研究场地特性。Ren等^[11]分析了2008年汶川地震的场地非线性效应，并提出了5个关键参数来表征场地非线性。Dhakal等^[12]利用K-NET六座海底台站记录到的315次地震事件分析了海底台站的场地效应，当PGA > 50gal时，海底部分台站所在场地出现非线性。Kubo等^[13]基于日本DONET1观测网记录到的2016年日本Off-Mie地震，利用HVSR法分析了台站所在场地的场地效应，认为部分台站场地的非线性反应抑制了地震动的放大效应。Ji等^[14]将HVSR计算出的非线性响应程度（DNL）和非线性响应百分比（PNL）联合PGA和V_S30（或V_S20）作为变量，采用k-means方法识别非线性场地响应。Yaghmaei-Sabegh等^[15]利用2017年伊朗Ezgeleh地震所记录到的地震数据，提出了基于广义回归神经网络（generalized regression neural networks，GRNN）的HVSR分类方法，该方法仅利用地震记录就可进行场地类别的划分，并将其与Vs30场地分类方法进行对比，展现了GRNN方法的可靠性。周旭彤等^[16]利用HVSR方法对日本DONET1台网20座海底台站场地效应进行分析，研究结果表明，部分台站具有相似的场地特征，海底台站的布设方式和海底地形会对谱比结果产生影响。

本文基于日本K-NET台网的海底台站收集的地震数据，利用HVSR联合动态聚类的方法对海底台站的场地特征进行识别，给出了台站所在场地的分类结果，为研究海域场地效应提供参考。

1. 海域强地震动观测数据与谱比方法

1.1 海域强地震动观测现状及数据处理

海域强地震动观测与陆域相比，存在设备安置难度大、回收效率低及观测稳定性较差等技术难题，发展远落后于现有陆域观测水平及海洋工程建设需求。美国和日本较早开展海域地震动观测，中国台湾地区也布设了海底地震观测设备^[17]，美国在20世纪70年代（1977年）针对海底地震记录观测需求，提出了SEMS计划（seafloor earthquake monitor system，以下简称SEMS）^[3]，日本于90年代（1996年）在东京圈近海开展ETMC计划（earthquake and tsunami monitoring cable，以下简称ETMC）^[18]。美国SEMS计划开发了4代海底地震仪，安装在南加州海域，水深区间为数10~300 m。日本ETMC系统布置在相模湾，观测系统由海底强震仪及海啸压力传感器组构，水深区间在900~2000 m不等。本文收集、整理了日本K-net台网6个海底台站2000—2018年的海域地震动数据。海底台站信息见表 1，分布见图 1。

表 1 海底台站信息

Table 1. Information of seabed stations

台站编号	东经/(°)	北纬/(°)	海水深度/m
KNG201	139.9183	34.5956	2197
KNG202	139.8393	34.7396	2339
KNG203	139.6435	34.7983	902
KNG204	139.5711	34.8931	933
KNG205	139.4213	34.9413	1486
KNG206	139.3778	34.0966	1130

下载: 导出CSV

| 显示表格

图 1 海底台站与水深分布图

Figure 1. Distribution map of submarine stations and water depth

下载: 全尺寸图片幻灯片

为研究6个台站所在海域的场地效应，考虑噪声对地震记录的影响，本文对上述数据进行筛选，筛选原则如下：①震级大于4级且PGA < 20 cm/s²；②震源位置均在海域。据此选取338组地震记录，每组记录均包含三分量，震级震中距分布见图 2。由图 2可知，震级主要分布在4.0~5.5级，震中距主要分布在0~300 km。

图 2 震级震中距分布图

Figure 2. Distribution map of magnitude and distance

下载: 全尺寸图片幻灯片

针对原始地震记录可能出现的基线漂移、仪器噪声及环境噪声等信号偏差，对筛选后的海域地震动记录采用基线校正和滤波的方法进行处理，其中采用Boore提出的分段式基线对地震动记录进行校正^[19]，构建4阶Butterworth带通滤波器，截止频率为0.1~20 Hz，处理前后的地震记录对比见图 3，由图可知，未处理地震动将出现基线漂移，中、长周期信号失真。处理后的地震记录更符合实际物理过程。

图 3 地震记录处理

Figure 3. Processing of seismic records

下载: 全尺寸图片幻灯片

1.2 海域强地震动记录HVSR计算

在利用强震记录研究场地效应时，需尽量选取S波时段^[20]。姚鑫鑫等^[21]关于S波的截取对HVSR谱比结果影响的研究表明，截取S波可以更好地体现高频段场地放大效应。因此本文对筛选后的海域地震动记录截取S波到时，将Arias强度能量^[22]达到90%所对应的持时作为S波的结束时间。

对截取后的S波数据进行快速傅里叶变换（FFT），选取0.2~20 Hz频带范围，并利用Konno平滑方法^[23]对傅里叶振幅谱进行平滑，平滑系数b取20。图 4（a）表明，利用该平滑方法既不低估场地放大效应的幅值，又不影响卓越频率的识别，可以较好地体现各频段的信息。将平滑后的地震记录S波的EW、NS向利用几何平均合成水平向，并除以相应的UD向平滑后的傅里叶振幅谱，得到的HVSR结果见图 4（b）。

图 4 HVSR计算过程

Figure 4. Calculation process of HVSR

下载: 全尺寸图片幻灯片

2. 动态聚类学习原理与评价

2.1 动态k-means聚类学习原理

Zhang等^[24]完善了动态 $k$ -means聚类方法，针对反应谱数据开展了聚类分析，本文通过添加谱比序列标签，采用动态聚类算法实现HVSR谱比分类。动态聚类学习原理如下：针对 $m$ 个数据序列 $\boldsymbol{Z} =$ $\left\{\boldsymbol{z}_1, \boldsymbol{z}_2, \cdots, \boldsymbol{z}_m\right\}^{\mathrm{T}} \in \mathbb{R}^{m \times n}$ ，各序列由 $n$ 个样本点组成行向量，即 ${{\bf{z}}_i} = \{ {z_{i1}}, {z_{i2}}, \cdot \cdot \cdot , {z_{in}}\}$ ，其中 $i = 1, 2, \cdot \cdot \cdot , m$ 。通过聚类、迭代方法，将数据序列实现分类，即将计算序列 ${\bf{Z}}$ 和各类数据中心 $\boldsymbol{C}=\left\{\boldsymbol{c}_1, \boldsymbol{c}_2, \cdots, \boldsymbol{c}_k\right\}^{\mathrm{T}} \in \mathbb{R}^{k \times n}$ 的欧式距离聚类：

$P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})=\sum\limits_{p=1}^k \sum\limits_{i=1}^m \sum\limits_{j=1}^n \lambda_{i p} w_{p j}\left(z_{i j}-c_{p j}\right)^2 \text{，}$

(1)

式中： $\boldsymbol{c}_p=\left\{c_{p 1}, c_{p 2}, \cdots, c_{p n}\right\}$ 是 $p$ 个类中心，由n个样本点组成的行向量。 $\boldsymbol{U} \in {\mathbb{Z}^{z \times k}}$ 为二元矩阵，当 ${\lambda _{ip}} = 1$ 表示第 $i$ 个序列分配给 $p$ 类，并且每个序列仅可分配给同类，即 $\sum\nolimits_{p = 1}^k {{\lambda _{ip}} = 1}$ 。通过设置权重矩阵 $\boldsymbol{W}$ ，可更好识别局部数据跨度内时间序列的相关性，提高聚类性能，每个序列的权重之和为1，即 $\sum\nolimits_{j = 1}^n {{w_{pj}} = 1}$ 。为了平滑在相邻序列间的权重，防止陷入局部最小，设置约束正则项，定义如下：

$Q_w(\boldsymbol{W})=\frac{\alpha}{2} \sum\limits_{p=1}^k \sum\limits_{j=1}^{n-1}\left(w_{p j}-w_{p(j+1)}\right)^2 \text{，}$

(2)

式中， $\alpha$ 为平滑参数， $\alpha = \sum\nolimits_{i = 1}^m {\sum\nolimits_{j = 1}^n {({z_{ij}} - \sum\nolimits_{i = 1}^m {{z_{ij}}/} } }$ $m{)^2}$ 。

综上，动态聚类可以被定义为具有正则化目标函数的优化问题：

$\begin{aligned} & \{\hat{\boldsymbol{U}}, \hat{\boldsymbol{W}}, \hat{\boldsymbol{C}}\}=\underset{\{\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C}\}}{\arg \min } P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})+Q_{\mathrm{w}}(\boldsymbol{W}) \\ & \mathrm{s}.\text{t}.\left\{\begin{array}{l}{\displaystyle \sum _{i=1}^{m}{\lambda }_{ip}\ge {m}_{0}}且{\displaystyle \sum _{p=1}^{k}{\lambda }_{ip}=1};{\lambda }_{ip}\in \{0, 1\}\\ {\displaystyle \sum _{j=1}^{n}{w}_{pj}=1;}\qquad \qquad\qquad 0\le {w}_{pj}\le 1\end{array}\right\} \text{，} \end{aligned}$

(3)

式中， ${m_0}$ 为每类所存在的最小序列数。 $\sum\nolimits_{i = 1}^m {{\lambda _{ip}} \geqslant {m_0}}$ 代表 $p$ 类中的每一类的序列数至少是m₀。式（3）的第一项使得所有类内的加权总距离最小，第二项使相邻类间隔的权重更平滑，即类间的区别更明显。通过设置计算序列到各中心的欧式距离，当此距离越近表示序列越相似。

2.2 动态聚类学习流程

该方法旨在最小化目标函数式（3），式（3）的变量由二元矩阵 $\boldsymbol{U}$ ，权重矩阵 $\boldsymbol{W}$ 和聚类中心矩阵 $\boldsymbol{C}$ 构成。采用优化方法，即将其中两个变量作为已知量，采用信息更新方法优化第三变量，当迭代次数达到40次，目标函数（3）收敛条件为10^-6时停止更新。整个优化过程可以分为如下3步。

（1）给定假设的权重矩阵 $\boldsymbol{W}$ 和聚类中心 $\boldsymbol{C}$ ，当式（1）满足式（4），序列被分配到假设聚类中心的类中，各类中心与序列之间的权重距离被最小化。

${\displaystyle \sum _{i=1}^{m}{\lambda }_{ip}\ge {m}_{0}}且{\displaystyle \sum _{p=1}^{k}{\lambda }_{ip}=1};{\lambda }_{ip}\in \{0, 1\} 。$

(4)

（2）当确定权重矩阵 $\boldsymbol{W}$ 和二元矩阵 $\boldsymbol{U}$ 后，通过令式（1）的导数为0求解信息更新后的新聚类中心，具体过程如下：

$\frac{\partial P(\boldsymbol{U}, \boldsymbol{W}, \boldsymbol{C})}{\partial c_{p j}}=2 \sum\limits_{i=1}^n \lambda_{i p} w_{p j}\left(z_{i j}-c_{p j}\right)=0 \text{，}$

(5)

变换得

${c_{pj}} = {{\sum\limits_{i = 1}^m {{\lambda _{ip}}{c_{ij}}} } \mathord{\left/ {\vphantom {{\sum\limits_{i = 1}^m {{\lambda _{ip}}{c_{ij}}} } {\sum\limits_{i = 1}^m {{\lambda _{ip}}} }}} \right. } {\sum\limits_{i = 1}^m {{\lambda _{ip}}} }} \text{，}$

(6)

通过式（6）可得到更新后的聚类中心。

（3）通过第一、二步确定二元矩阵 $\boldsymbol{U}$ 和聚类中心 $\boldsymbol{C}$ ，此时当式（3）满足

$\sum\limits_{j = 1}^n {{w_{pj}} = 1} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} (0 \leqslant {w_{pj}} \leqslant 1) 。$

(7)

式（3）转化为计算权重矩阵 $\boldsymbol{W}$ 的二次规划问题，具体求解方法可以参考文献[25]。

2.3 动态k-means聚类方法验证

为验证该方法在HVSR分类中的性能，选取5%阻尼比的反应谱来进行方法验证。依据《建筑抗震设计规范：GB50011—2010》^[26]，在场地分类为III类条件下，生成两组设计反应谱，两类反应谱见图 5，图 5中第一类反应谱的特征周期为0.9 s，地震影响系数最大值为0.9；第二类反应谱特征周期为0.3 s，地震影响系数最大值为0.79。

图 5 设计反应谱

Figure 5. Design response spectra

下载: 全尺寸图片幻灯片

依据上述不同特性反应谱，从PEER强震动数据库中^[27]，分别针对各分类选取50条地震记录。图 6给出了数据的选取情况。

图 6 不同记录反应谱特性

Figure 6. Two clusters of response spectra

下载: 全尺寸图片幻灯片

将两组数据利用动态聚类学习方法进行分类测试，结果见图 7。由图 7可知，采用动态聚类学习方法划分的中心和反应谱均值无偏差，展现出良好的聚类性能及鲁棒性。

图 7 反应谱均值与中心匹配

Figure 7. Mean values of response spectra and cluster centroid

下载: 全尺寸图片幻灯片

3. 海域地震动谱比法分类

3.1 海域地震动谱比法分类依据

将HVSR结果利用上述动态聚类方法进行聚类。利用误差平方和准则评价聚类性能，表示聚类的畸变程度。给定数据序列 $\boldsymbol{Z}$ ，假设 $\boldsymbol{Z}$ 中的序列可分为 $k$ 类，即 ${\boldsymbol{Z}_1}, {\boldsymbol{Z}_2}, \cdot \cdot \cdot , {\boldsymbol{Z}_k}$ ；各类序列数分别为 ${m_1}, {m_2}, \cdot \cdot \cdot , {m_k}$ ；各类中心分别为 ${c_1}, {c_2}, \cdot \cdot \cdot , {c_k}$ ，则误差平方和准则函数为

$D = \sum\nolimits_{i = 1}^k {\sum\limits_{z \in {Z_i}} {{{(d(z, {c_i}))}^2}} } \text{，}$

(8)

式中， $D$ 为数据序列分为 $k$ 类时的每类数据序列与其中心序列的距离和。 $D$ 越小表示每类数据序列到对应中心的距离越小，越紧凑，分类性能越好。分别令 $k = 1, 2, \cdot \cdot \cdot , 8$ 进行聚类计算并计算 $D$ ，计算结果见。从中可以看出，随着分类数增大，误差平方和趋于稳定，结合肘部法则^[28]确定最优分类数为 $k = 3$ 。此外为保证每类结果具有统计意义，设置每类最少存在6组HVSR数据序列，即 ${m_0} = 6$ 。

图 8 误差平方和与分类数关系

Figure 8. Identification of optimal number of cluster at the elbow point

下载: 全尺寸图片幻灯片

3.2 海域场地特性分析

HVSR动态聚类结果见图 9，对应台站统计如图 10所示。由图 9，10可知：海域场地HVSR曲线无明显峰值或多峰值现象，结合肘部算法，相模湾海区最优化场地可分为三类，结合相模湾海区地貌^[29-30]，综合分为海底峡谷区（KNG203，KNG204）及海洋盆地区（KNG201，KNG202，KNG205，KNG206），两类地貌分类的HVSR谱比呈显著差异，海底峡谷区高频衰减慢且优势频段主要集中在中高频，台站KNG204代表的第一类主要集中在12 Hz，体现了海底峡谷中少量沉积的薄覆盖层环境；台站KNG203代表的第三类优势频带跨度较大（5~10 Hz），台站KNG201，KNG202，KNG205，KNG206则表现为中长周期发育，高频衰减快，呈现出典型的深厚覆盖层的海洋盆地特性。

图 9 HVSR聚类结果

Figure 9. Clustering results of HVSR

下载: 全尺寸图片幻灯片

图 10 各海域台站HVSR与类别对应

Figure 10. HVSR of each station corresponding to category

下载: 全尺寸图片幻灯片

其次，局部场地地形效应也显著影响海域地震动传播，以台站作为中心点，以类似地貌作为依据，通过采样点与标准点之间的距离（500 m）与高程计算局部场地地形坡度：

$i = h/l \times 100\% \text{，}$

(9)

式中，h为采样点与标准点的平均高程差，l为采样点与标准点距离。

各台站所在局部场地地形坡度见表 2。由表 2可知：海底峡谷区台站KNG203较KNG204场地局部地形平缓度高，存在局部沉积形成的浅覆盖层；海洋盆地区台站KNG201，KNG205，KNG206所在场地局部地形坡度基本一致，KNG202所在场地处于盆地边缘，存在基岩面拉伸导致的局部场地不均匀分布特性。

表 2 局部场地地形坡度

Table 2. Topographic slopes of local site

类别	台站	地貌	坡度	地形	场地特征
一	KNG204	海洋峡谷	4.62	陡峭	薄覆盖层
二	KNG201	海洋盆地	1.39	平缓	厚覆盖层
	KNG205		1.85
	KNG206		1.02
	KNG202		3.81	陡峭	薄覆盖层
三	KNG203	海洋峡谷	0.38	平缓	浅覆盖层

下载: 导出CSV

| 显示表格

4. 结论

本文采用动态聚类学习方法结合海域场地HVSR谱比，分析了动态聚类方法在地震动分类的适用性；给出了相模湾海区场地工程特性，结合海洋工程地质给出了海域地形坡度与场地分类关系。希望通过对有限的海域场地地震记录研究，为近海以及海洋工程结构和设施的设计提供合理、科学的工程场地参考。得到以下2点结论。

（1）动态聚类学习方法针对时序数据的聚类性能及鲁棒性强，可满足局部及全局分辨率，应用于场地工程特性分析中，可结合HVSR谱比，最优化、自动化的给出场地分类特性；

（2）HVSR动态聚类的海洋场地特性表明：台站KNG203和KNG204位于海底峡谷区，台站KNG201、KNG202、KNG205、KNG206位于海洋盆地区，两类地貌分类对海域地震动传播影响较大，海底峡谷区高频衰减慢且优势频段主要集中在中高频，通过局部场地地形坡度分析给出了各台站所在场地的工程特性。

图 1 海底台站与水深分布图

Figure 1. Distribution map of submarine stations and water depth

下载: 全尺寸图片幻灯片

图 2 震级震中距分布图

Figure 2. Distribution map of magnitude and distance

下载: 全尺寸图片幻灯片

图 3 地震记录处理

Figure 3. Processing of seismic records

下载: 全尺寸图片幻灯片

图 4 HVSR计算过程

Figure 4. Calculation process of HVSR

下载: 全尺寸图片幻灯片

图 5 设计反应谱

Figure 5. Design response spectra

下载: 全尺寸图片幻灯片

图 6 不同记录反应谱特性

Figure 6. Two clusters of response spectra

下载: 全尺寸图片幻灯片

图 7 反应谱均值与中心匹配

Figure 7. Mean values of response spectra and cluster centroid

下载: 全尺寸图片幻灯片

图 8 误差平方和与分类数关系

Figure 8. Identification of optimal number of cluster at the elbow point

下载: 全尺寸图片幻灯片

图 9 HVSR聚类结果

Figure 9. Clustering results of HVSR

下载: 全尺寸图片幻灯片

图 10 各海域台站HVSR与类别对应

Figure 10. HVSR of each station corresponding to category

下载: 全尺寸图片幻灯片

表 1 海底台站信息

Table 1 Information of seabed stations

台站编号	东经/(°)	北纬/(°)	海水深度/m
KNG201	139.9183	34.5956	2197
KNG202	139.8393	34.7396	2339
KNG203	139.6435	34.7983	902
KNG204	139.5711	34.8931	933
KNG205	139.4213	34.9413	1486
KNG206	139.3778	34.0966	1130

下载: 导出CSV

表 2 局部场地地形坡度

Table 2 Topographic slopes of local site

类别	台站	地貌	坡度	地形	场地特征
一	KNG204	海洋峡谷	4.62	陡峭	薄覆盖层
二	KNG201	海洋盆地	1.39	平缓	厚覆盖层
	KNG205		1.85
	KNG206		1.02
	KNG202		3.81	陡峭	薄覆盖层
三	KNG203	海洋峡谷	0.38	平缓	浅覆盖层

下载: 导出CSV

参考文献(30)

[1]	李小军, 李娜, 陈苏. 中国海域地震区划及关键问题研究[J]. 震灾防御技术, 2021, 16(1): 1-10. https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101001.htm LI Xiaojun, LI Na, CHEN Su. Study on seismic zoning in China Sea area and its key issues[J]. Technology for Earthquake Disaster Prevention, 2021, 16(1): 1-10. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101001.htm
[2]	陈苏, 周越, 李小军, 等. 近海域地震动的时频特征与工程特性[J]. 振动与冲击, 2018, 37(16): 227-233. https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201816033.htm CHEN Su, ZHOU Yue, LI Xiaojun, et al. Time-frequency and engineering characteristics on offshore ground motion[J]. Journal of Vibration and Shock, 2018, 37(16): 227-233. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZDCJ201816033.htm
[3]	胡进军, 刁红旗, 谢礼立. 海底强地震动观测及其特征的研究进展[J]. 地震工程与工程振动, 2013, 33(6): 1-8. https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC201306001.htm HU Jinjun, DIAO Hongqi, XIE Lili. Review of observation and characteristics of seafloor strong motion[J]. Journal of Earthquake Engineering and Engineering Vibration, 2013, 33(6): 1-8. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC201306001.htm
[4]	NAKAMURA Y. Method for dynamic characteristics estimation of subsurface using microtremor on the ground surface[J]. Quarterly Report of RTRI (Railway Technical Research Institute) (Japan), 1989, 30(1): 25-33.
[5]	BOORE D M, SMITH C E. Analysis of earthquake recordings obtained from the Seafloor Earthquake Measurement System (SEMS) instruments deployed off the coast of southern California[J]. Bulletin of the Seismological Society of America, 1999, 89(1): 260-274. doi: 10.1785/BSSA0890010260
[6]	FIELD E H, JOHNSON P A, BERESNEV I A, et al. Nonlinear ground-motion amplification by sediments during the 1994 Northridge earthquake[J]. Nature, 1997, 390(6660): 599-602. doi: 10.1038/37586
[7]	BONILLA L F. Borehole response studies at the garner valley downhole array, southern California[J]. Bulletin of the Seismological Society of America, 2002, 92(8): 3165-3179. doi: 10.1785/0120010235
[8]	WEN K L, CHANG T M, LIN C M, et al. Identification of nonlinear site response using the H/V spectral ratio method[J]. Terrestrial, Atmospheric and Oceanic Sciences, 2006, 17(3): 533. doi: 10.3319/TAO.2006.17.3.533(T)
[9]	任叶飞, 温瑞智, 山中浩明, 等. 运用广义反演法研究汶川地震场地效应[J]. 土木工程学报, 2013, 46(增刊2): 146-151. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2013S2025.htm REN Yefei, WEN Ruizhi, YAMANAKA H, et al. Research on site effect of Wenchuan Earthquake by using generalized inversion technique[J]. China Civil Engineering Journal, 2013, 46(S2): 146-151. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2013S2025.htm
[10]	荣棉水, 李小军, 王振明, 等. HVSR方法用于地震作用下场地效应分析的适用性研究[J]. 地球物理学报, 2016, 59(8): 2878-2891. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201608014.htm RONG Mianshui, LI Xiaojun, WANG Zhenming, et al. Applicability of HVSR in analysis of site-effects caused by earthquakes[J]. Chinese Journal of Geophysics, 2016, 59(8): 2878-2891. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201608014.htm
[11]	REN Y F, WEN R Z, YAO X X, et al. Five parameters for the evaluation of the soil nonlinearity during the Ms8.0 Wenchuan Earthquake using the HVSR method[J]. Earth, Planets and Space, 2017, 69(1): 1-17. doi: 10.1186/s40623-016-0587-x
[12]	DHAKAL Y P, AOI S, KUNUGI T, et al. Assessment of nonlinear site response at ocean bottom seismograph sites based on S-wave horizontal-to-vertical spectral ratios: a study at the Sagami Bay area K-NET sites in Japan[J]. Earth, Planets and Space, 2017, 69(1): 1-7. doi: 10.1186/s40623-016-0587-x
[13]	KUBO H, NAKAMURA T, SUZUKI W, et al. Ground‐motion characteristics and nonlinear soil response observed by DONET1 seafloor observation network during the 2016 southeast off‐Mie, Japan, earthquake[J]. Bulletin of the Seismological Society of America, 2019, 109(3): 976-986.
[14]	JI K, WEN R, REN Y, et al Nonlinear seismic site response classification using K-means clustering algorithm: case study of the September 6, 2018 Mw6.6 Hokkaido Iburi-Tobu earthquake, Japan[J]. Soil Dynamics and Earthquake Engineering, 2020, 128: 105907.
[15]	YAGHMAEI-SABEGH S, RUPAKHETY R. A new method of seismic site classification using HVSR curves: a case study of the 12 November 2017 Mw 7.3 Ezgeleh earthquake in Iran[J]. Engineering Geology, 2020, 270: 105574.
[16]	周旭彤, 胡进军, 谭景阳, 等. 基于HVSR的DONET1海底地震动场地效应研究[J]. 震灾防御技术, 2021, 16(1): 105-115. https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101011.htm ZHOU Xutong, HU Jinjun, TAN Jingyang, et al. The study of site effect of DONET1 offshore ground motions based on HVSR[J]. Technology for Earthquake Disaster Prevention, 2021, 16(1): 105-115. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101011.htm
[17]	赵纪东. 中国台湾将建设首座海底地震台站[J]. 地球科学进展, 2010, 25(2): 146. https://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ201002005.htm ZHAO Jidong. China will build the first submarine seismic station[J]. Advances in Earth Science, 2010, 25(2): 146. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DXJZ201002005.htm
[18]	EGUCHI T, FUJINAWA Y, FUJITA E, et al. A real-time observation network of ocean-bottom-seismometers deployed at the Sagami trough subduction zone, central Japan[J]. Marine Geophysical Researches, 1998, 20(2): 73-94.
[19]	BOORE D M, BOMMER J J. Processing of strong-motion accelerograms: needs, options and consequences[J]. Soil Dynamics and Earthquake Engineering, 2005, 25(2): 93-115.
[20]	李小军, 李娜, 王巨科, 等. 场地地震动水平/竖向谱比与地表/基底谱比差异及修正水平/竖向谱比法研究[J]. 震灾防御技术, 2021, 16(1): 81-90. https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101009.htm LI Xiaojun, LI Na, WANG Juke, et al. Difference between horizontal-to-vertical spectral ratio and surface-to-bedrock spectral ratio of strong-motion and modified horizontal-to-vertical spectral ratio method[J]. Technology for Earthquake Disaster Prevention, 2021, 16(1): 81-90. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY202101009.htm
[21]	姚鑫鑫, 任叶飞, 温瑞智, 等. 强震动记录H/V谱比法计算处理的若干关键环节[J]. 震灾防御技术, 2019, 14(4): 719-730. https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY201904003.htm YAO Xinxin, REN Yefei, WEN Ruizhi, et al. Some technical notes on the data processing of the spectral ratio based on the strong-motion records[J]. Technology for Earthquake Disaster Prevention, 2019, 14(4): 719-730. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZZFY201904003.htm
[22]	ARIAS A. Measure of Earthquake Intensity[R]. Massachusetts: Massachusetts Inst of Tech, 970.
[23]	KONNO K, OHMACHI T. Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor[J]. Bulletin of the Seismological Society of America, 1998, 88(1): 228-241.
[24]	ZHANG R Y, HAJJAR J, SUN H. Machine learning approach for sequence clustering with applications to ground-motion selection[J]. Journal of Engineering Mechanics, 2020, 146(6): 04020040.
[25]	FRANK M, WOLFE P. An algorithm for quadratic programming[J]. Naval Research Logistics Quarterly, 1956, 3(1/2): 95-110.
[26]	建筑抗震设计规范: GB50011—2010 [S]. 北京: 中国建筑工业出版社, 2010. Code for Seismic Design of Buildings: GB50011—2010[M]. Beijing: China Architecture & Building Press, 2010. (in Chinese).
[27]	CHIOU B, DARRAGH R, GREGOR N, et al. NGA project strong-motion database[J]. Earthquake Spectra, 2008, 24(1): 23-44.
[28]	KODINARIYA T, MAKWANA P R. Review on determining number of Cluster in K-Means Clustering[J]. International Journal, 2013, 1(6): 90-95.
[29]	EGUCHI T, FUJINAWA Y, FUJITA E, et al. A real-time observation network of ocean-bottom-seismometers deployed at the Sagami trough subduction zone, central Japan[J]. Marine Geophysical Researches, 1998, 20(2): 73-94.
[30]	UTASHIRO S, IWABUCHI Y. The submarine topography and geological structure of Sagami Bay, south coast of Honshu, Japan[J]. Journal of Geography (Chigaku Zasshi), 1971, 80(2): 77-88.

施引文献(15)

期刊类型引用(13)

1.	侯鹏，刘尧，黄毅. 地铁区间隧道防沉“微扰动”注浆加固技术研究. 城市建设理论研究(电子版). 2025(02): 173-175 . 百度学术
2.	梁发云，魏圣明，陈可，袁强，闫静雅，顾晓强. 基于UH模型的盾构隧道非线性地基反力确定方法. 岩土工程学报. 2025(02): 284-295 . 本站查看
3.	张玄，肖颖，梁寿，王禹玺，王睿，胡黎明. 地热源作用下盾构地铁隧道管片应力响应规律. 公路. 2025(02): 443-449 . 百度学术
4.	梁荣柱，张志伟，李忠超，符宇坤，曹世安，康成，肖铭钊，吴文兵. 考虑环间接头影响的盾构隧道纵向变形简化解. 岩土工程学报. 2024(02): 254-263 . 本站查看
5.	王朋路，王微笑，冯畅. 城市砂土隧道注浆效果评价研究. 武汉理工大学学报(交通科学与工程版). 2024(01): 180-184 . 百度学术
6.	申子龙. 地铁盾构隧道掘进过程注浆技术研究. 工程机械与维修. 2024(02): 162-164 . 百度学术
7.	王观璐，徐德鹏，张辉. 市政配套工程施工对既有地铁车站的影响分析. 现代城市轨道交通. 2024(04): 97-105 . 百度学术
8.	黄大维，罗仲睿，罗文俊，朱碧堂，刘家璇，赵梽錡. 地层注浆附加土压力形成及影响因素单元体试验研究. 岩石力学与工程学报. 2024(S1): 3520-3529 . 百度学术
9.	袁钊. 城市复杂环境条件下盾构隧道变形治理监控关键措施. 建设监理. 2024(07): 131-134 . 百度学术
10.	刘禹，刘性锋，李世堂，周静. 考虑多因素的深埋管线施工对临近地铁影响分析——以龙泉驿区玉竹路为例. 科技和产业. 2024(18): 327-332 . 百度学术
11.	王虎，叶胜，罗天霖，宋敦安，李栋. 深厚软土区地表附加荷载对顶管隧道变形的影响分析. 现代隧道技术. 2024(S1): 1011-1017 . 百度学术
12.	袁强，赵旭，陈诗艾，李建平，袁炳祥. 盾构隧道纵缝接头抗弯性能评价指标及极限承载分析. 现代隧道技术. 2024(S1): 432-444 . 百度学术
13.	涂玖林，王春森，何永辉. 市政道路斜交上跨既有地铁隧道影响研究. 交通节能与环保. 2023(05): 275-279+284 . 百度学术