Study on Azimuth-related Site Amplification Effect of Seismic Ground Motion
-
摘要: 现有场地对地震动的放大作用分析中,常采用观测得到正交方向的HVSR和SBSR的矢量合成结果表征场地放大作用,而不考虑放大作用的方向相关性。为研究不同地震作用下场地放大作用的方向相关性,采用多向HVSR和多向SBSR及谱比峰值方向相关分布概率,以11次地震事件获得的地表、井下基岩记录作为数据,开展响嘡场地效应台阵场地对地震动放大作用的方向相关性分析。研究结果表明,常用的矢量合成法分析结果与多向谱比上包络曲线法分析结果存在差异,前者存在忽略部分峰值的可能;多向谱比的极坐标等高线图可较直观地体现场地对地震动放大作用的方向相关性特征,但不同多向谱比表现存在差异;多向谱比峰值方向分布概率可较好地定量分析放大作用峰值在不同方向上的差异,但不同地震事件存在差异,不同的多向谱比存在差异。可采用多向谱比上包络曲线法分析场地放大作用,除可避免矢量合成法引起的误估外,也可更好地分析实际场地频率-方向相关的放大作用。
-
关键词:
- 多向HVSR /
- 多向SBSR /
- 谱比包络 /
- 场地效应 /
- 峰值方向相关分布概率
Abstract: In recent analysis on site amplification effect of seismic ground motion, the square root of HVSR and SBSR obtained from seismic observation in orthogonal direction was often used to characterize the site amplification effect, however, the azimuth-related site amplification effect of seismic ground motion is not considered. To discuss this effect, the azimuth-related HVSRs, SBSRs, and distribution probability of peak values on the spectral ratios were selected, and the surface and borehole seismic records on Xiangtang Array site in 11 earthquake events were collected, finally, the azimuth-related site amplification effect was analyzed. The results show: (1) the results from the square root algorithm have the risk of ignoring some frequencies which are different with the proposed upper envelope algorithm in this work, (2) the contour of azimuth-related spectral ratios in polar coordinates could directly represent the directionality characteristics of site amplification effect, but the results from different multi-directional spectral ratios are different, (3) the azimuth-related of peak values could characterize in quantity the difference in azimuth of site amplification effect, but the probabilities were various with the different multi-directional spectral ratio and the data from different event. According to the results and discussion in this work, the azimuth-related spectral ratios could be used in the analysis of site amplification effect, for avoiding the under-estimate of amplification effect by the square root algorithm, and the method could be used to understand the frequency-direction correlation of site amplification effect better. -
引言
建筑类型分类是抗震设计的基础,也是地震风险研究的依据。对不同建筑物进行分类研究,编制建筑分类清单并进行易损性分析,根据所得破坏概率矩阵,可对建筑群及城市风险进行预测和评估。面向地震保险的房屋建筑分类需满足区域化、高效化及信息化等要求,与传统的建筑分类方式存在较大区别。在地震保险应用中,城市建筑物的整体情况是至关重要的数据资料,相比通过人力调查,利用遥感影像获取相关信息具有省时省力的优势。随着遥感技术的成熟,遥感影像的分辨率已经非常高,其图像获取及应用对地震保险有重要意义。遥感数据的不断丰富为地震保险领域的研究提供了数据支持,但从海量数据中提取有用信息也是一项挑战。面对海量图像数据,依靠人工对图像进行分类及标注的管理方式需耗费大量人力资源,采用计算机自动识别的方法显得尤为重要。利用传统图像处理方法提取的图像底层特征所建模型的类别语义与用户高层语义之间的语义鸿沟,使得计算机图像的自动解译变得异常困难,如何跨越语义鸿沟,使计算机不仅能独立完成像素(区域)级的分类,而且能学习图像的高层语义,自动完成图像级的分类和识别,成为该学科领域的一个热点和难点问题。主流遥感图像的分辨率已达分米级,以往行之有效的一系列针对中低分辨率的模型及方法在高分辨率下效果并不理想。对高分辨率遥感影像中的目标进行匹配、分类、识别等操作,需要更多地考虑诸如纹理、统计量分布、几何形状、空间上下文等计算机视觉特征。同时,由于遥感图像分辨率的提高,图像细节不断增多,结构复杂性也随之增加,可提取的结构特性也将更加丰富。
深度学习的概念源于人工神经网络研究,含多隐层的多层感知器(MLP)就是一种深度学习结构。深度学习通过组合低层特征形成更加抽象的高层表示(属性类别或特征),以发现数据的分布式特征表示。近几年,深度学习在机器学习领域取得重大突破,使得图像识别、语音识别精度大幅度提高(Mikolov等,2011;Ciodaro等,2012;Hinton等,2012;Krizhevsky等,2012;Farabe等,2013;Helmstaedter等,2013;Sainath等,2013;Tompson等,2014;Szegedy等,2014;Ma等,2015)。这方面的发展主要基于算法的逐渐成熟和数据的不断丰富。深度结构(涉及多个非线性处理单元层)非凸目标代价函数中普遍存在的局部最小是训练困难的主要来源。基于深信度网(DBN)提出非监督贪心逐层训练算法(Bengio等,1994, 2005),此外提出的卷积神经网络(CNNS)是第一个真正多层结构学习算法,它利用空间相对关系减少参数数目以提高BP训练性能(Le Cun等,1990)。利用深度学习技术从遥感影像中提取结构特性的相关研究较少,现有研究主要集中在利用固定的图像处理算法从遥感影像中提取建筑物相关信息,如利用阴影几何模型和阴影特征实现阴影检测,基于阴影的强度和几何形状实现影像的识别和分类(Jiang等,1994)。夏浩铭等通过提取纹理特征结合BP神经网络实现对遥感地物进行判定分类(夏浩铭等,2012);利用支持向量机的方法,通过训练分类器实现阴影区域的检测(Choi等,2010)。现有建筑物结构分类方法主要有现场调查、遥感影像人工解译、影像光谱特征分析、影像多特征融合等(李强等,2016;杜浩国等,2018;李金香等,2019),与机器学习方法相比人工干预较为显著。采用较为先进的图像分类及处理技术,可以提高分类精度及效率。
本文基于深度学习方法构建建筑高度、建筑功能和建筑结构类型识别模型,完成基于机器视觉的建筑高度模型;运用训练好的基于Xception神经网络深度学习模型和基于机器视觉的模型,对北京地区的建筑高度进行模型测试,此方法可为地震保险分析提供重要的基础数据支持。
1. 基于深度学习的建筑信息识别模型
1.1 数据来源及深度学习网络识别模型
本文所采用的街景图像主要来自于北京市主城区,获取内容如下:
(1)目标城市主城区道路分布图(shp格式),主要包括空间位置、道路长度、道路级别等信息,坐标系为WGS84。主城区建筑分布图(shp格式),包括空间位置信息(经纬度、轮廓范围,坐标系为WGS84)和房屋属性信息(编号、存在状态、地址、房屋层数、结构类型、建筑年代、面积、周长等)。
(2)选取所需关注因素:房屋编号、经纬度、房屋层数、结构类型、建筑年代、面积、周长等,对相应数据进行筛选和清洗。
(3)利用百度地图JavaScript API,将地理信息数据由WGS84坐标转化为百度地图BD09中可以使用的地理坐标。
(4)考虑建筑物中心点坐标以及全景的位置,选择视角、视野等参数,利用百度地图中的全景静态图服务批量获取街景图。建筑数据来源及入库流程如图 1所示。
本项目对建筑高度的识别主要通过提取建筑物的层高特征进行判断,将建筑物按照层高分为1层,2-3层,4-6层以及7层以上。其中训练集为15000张,验证集为1000张,测试集为895张,通过将训练集输入至深度学习模型进行计算,最后模型在测试集的准确率为52.4%,识别受限于树木、道路遮挡等状况。本文通过Xception神经网络中的卷积神经网络进行建筑特征图像识别,包括Entry flow、Middle flow、Exit flow三部分,其中Entry flow包含8个conv,Middle flow包含3×8=24个conv,Exit flow包含4个conv,共36层,层与层之间采用残差形式连接。Xception神经网络模型识别建筑信息处理流程如图 2所示。
1.2 模型参数对结果的影响
本文构建的机器学习模型共有5个卷积层,图 3(a)、(b)、(c)分别为冻结全部卷积层参数、冻结前四个卷积层参数、冻结前三个卷积层参数后的计算结果。表 1是冻结不同卷积层参数的结果对比。由图可知,当冻结全部卷积层参数时,只训练全连接层数,网络是发散的,因此训练获得参数无法使用;当冻结前四个卷积层参数或冻结前三个卷积层参数时,均可取得较好的效果,测试精度和测试误差相似,但冻结前四个卷积层参数相比冻结前三个卷积层参数,其收敛更快,能更快达到期望的效果。
表 1 冻结不同卷积层参数计算结果Table 1. The experimental results of freezing different layer parameters冻结层数 每次迭代时间/s 测试集误差 测试准确度 测试误差变化趋势 冻结全部卷积层参数 37 网络发散 — — 冻结前四层卷积层参数 36 0.4400 0.8342 在40次迭代后趋于稳定 冻结前三层卷积层参数 37 0.4189 0.8459 在60次迭代后趋于稳定 2. 基于机器视觉的建筑信息识别模型
2.1 图像来源及处理方法
首先从OSM网站2下载北京地区的地理信息,并将下载的OSM格式数据转化为json格式,再从已有信息中提取所有建筑物坐标和建筑物类型,共计57302栋建筑物。对相应数据进行处理,可将建筑物信息保存为包括建筑物角点坐标、建筑物中心位置坐标的一系列文件,建筑物中心位置坐标通过建筑物角点坐标计算得到;其次,通过百度地图开放网站获取建筑物对应的相机经纬度坐标、偏航角和俯仰角等,通过API接口批量获取对应建筑物的街景图。经统计,存在相机位置的建筑物数量为24437栋(包含损坏图片)。按照相机投影原理,将OSM数据中建筑物的轮廓坐标投影到街景图中,根据两者之间的匹配关系确定建筑物高度,从OSM数据可获得正阳门轮廓中四个角点的经纬度坐标分别为(x1,y1)、(x2,y2)、(x3,y3)、(x4,y4),及对应的相机位置为(xc,yc),根据相机位置和建筑物位置,选取四个角中的两个角点作为投影坐标,用P表示:
2 https://www.openstreetmap.org
$$\boldsymbol{P}{\rm{ = }}\left[ \begin{gathered} {x_{\rm{1}}}, {y_{\rm{1}}} \\ {x_{\rm{2}}}, {y_{\rm{2}}} \\ \end{gathered} \right]$$ (1) 获取的街景图宽1024像素,长512像素,偏航角α由建筑物中心位置和相机的相对位置确定,俯仰角β统一设置为0°,翻滚角γ为0°,确定相机外参矩阵及内参矩阵,得到相机投影坐标。外参矩阵和内参矩阵如式(2)—(6)所示:
$$ \boldsymbol{R}=\boldsymbol{R}_{y} \ · \ \boldsymbol{R}_{x} \ · \ \boldsymbol{R}_{z}$$ (2) $${\boldsymbol{R}_z} = \left[ {\begin{array}{*{20}{c}} {\cos (\alpha)}&{ - \sin (\alpha)}&0 \\ {\sin (\alpha)}&{\cos (\alpha)}&0 \\ {0}&{ 0}&1 \end{array} } \right]$$ (3) $${\boldsymbol{R}_x} = \left[ {\begin{array}{*{20}{c}} 1&0&0 \\ 0&{\cos (\beta)}&{ - \sin (\beta)} \\ 0&{\sin ((\beta))}&{\cos (\beta)} \end{array}} \right]$$ (4) $${\boldsymbol{R}_y} = \left[ {\begin{array}{*{20}{c}} {\cos (\gamma)}&0&{\sin (\gamma)} \\ {0}&1&{0} \\ { - \sin (\gamma)}&0&{\cos (\gamma)} \end{array} } \right]$$ (5) $$\boldsymbol{K}=\left[ \begin{matrix} {}^{f}\!\!\diagup\!\!{}_{{{d}_{x}}}\; & \ \ \ 0 & {{u}_{0}} \\ 0 & {}^{f}\!\!\diagup\!\!{}_{{{d}_{y}}}\; & {{v}_{0}} \\ 0 & \ \ \ 0 & 1 \\ \end{matrix} \right]$$ (6) 式中,f为焦距,dx和dy为尺度因子。
世界坐标到相机坐标的映射关系如式7所示:
$$\lambda p = {\rm{[}}K{\rm{|}}{\boldsymbol{O}_{\rm{3}}}{\rm{]}}\left[ {\begin{array}{*{20}{c}} \boldsymbol{R}&{ - \boldsymbol{R}C} \\ {\boldsymbol{O}_{\rm{3}}^{\rm{T}}}&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} P \\ 1 \end{array} } \right]$$ (7) 式中,$\lambda p$为投影到街景图上的坐标,P为投影坐标(辅以高度信息),K为相机内参数,C为相机中心坐标,O3表示零矩阵。
2.2 模型识别效果
为测试并验证本文模型利用街景图识别建筑高度的准确性和适用性,从OSM数据库下载北京市的建筑轮廓数据,并通过上述模型对采集的百度街景图进行识别计算,共获得14056栋建筑的高度。为对比模型计算结果与第三方数据公司提供的近似真实建筑高度的差别,本文按照楼层数量对建筑进行分类分析,具体分为1-9层(高度<30米)和10层及以上(高度≥30米)两类。模型测算高度与真实高度的相关性和误差如图 4-6所示,其中图 4为所有建筑数据集合,图 5、图 6分别为不同高度数据集合。由图可知,模型测算高度与误差之间有较好的相关性。在建筑高度为1-9层(小于30米)时,机器学习获取的建筑高度与误差之间的线性拟合相关性R2达到0.9337,随着建筑高度增加,本模型测算出的建筑高度误差增大,即本模型在识别建筑高度时,对于中低层建筑识别精度较高,建筑过高会增加识别的误差。
3. 结论
本文采用深度学习方法从影像数据中提取面向地震保险需求的建筑特性数据,构建基于深度学习方法的建筑高度识别模型,运用训练好的基于Xception神经网络深度学习模型对北京地区的建筑高度进行测试。结果表明,从技术方法角度,识别是可行的,但从数据源角度,街景图像的背景复杂,车辆、行人和树木可能会干扰提取图像特征,如何进行图像去噪将是后期研究的重点;采用机器视觉技术进行建筑高度识别,在建筑高度小于30米时,模型测算高度与误差之间的线性拟合相关性较高,其线性拟合给出的是误差的平均值,误差的平均值不代表识别精度的高低,主要反应的是整体估计的偏差。楼层越低,模型总体估计越大,楼层越高,模型总体估计越小。新方法和传统方法的结合使用可为地震保险分析提供重要的基础数据支持。
-
表 1 本文选用记录对应的地震事件相关参数
Table 1. Parameter of the earthquake events selected in this work
发震时间/(年-月-日) 震中位置 震级ML 震源深度/km 震源关系 PGA/Gal 纬度/° 经度/° 地点 震中距离/km 方位角/° EW向 NS向 2001-07-11 北纬39.75 东经118.28 陡河 3.5 5 36.4 278.9 6.22 8.30 2001-10-22 北纬39.78 东经118.70 滦县 3.1 — 8.9 0 18.41 11.58 2001-12-28 北纬39.67 东经118.65 滦县 4.2 12 5.4 232.2 97.96 161.71 2002-01-30 北纬39.68 东经118.25 唐山 3.3 — 38.7 266.9 3.10 4.75 2002-02-09 北纬39.77 东经118.37 陡河 4.0 — 29.3 285.5 13.61 18.51 2002-05-06 北纬39.80 东经118.77 滦县 3.8 — 12.6 28.4 48.80 60.00 2002-05-19 北纬39.37 东经117.98 丰南 4.7 — 71.9 239.6 5.84 6.83 2002-06-06 北纬39.70 东经118.75 滦县 2.7 6 4.3 90.0 5.80 6.90 2002-06-20 北纬39.88 东经118.78 卢龙 2.8 9 21.1 18.9 8.95 7.32 2002-07-20 北纬39.85 东经118.75 沙河驿 3.3 10 17.2 14.4 6.30 12.40 2002-11-13 北纬39.77 东经118.63 滦县 3.0 15 9.8 322.4 8.65 7.55 -
李红光, 冷崴, 2019. H/V谱比方法在强震记录应用中几个影响因素的初步研究. 地震工程学报, 41(1): 169—176 doi: 10.3969/j.issn.1000-0844.2019.01.169Li H. G. , Leng W. , 2019. Preliminary study on the influencing factors of the horizontal-to-vertical spectral ratio applied in strong motion. China Earthquake Engineering Journal, 41(1): 169—176. (in Chinese) doi: 10.3969/j.issn.1000-0844.2019.01.169 卢滔, 周正华, 周雍年等, 2006. 关于Nakamura方法有效性的讨论. 地震工程与工程振动, 26(1): 43—48 doi: 10.3969/j.issn.1000-1301.2006.01.007Lu T. , Zhou Z. H. , Zhou Y. N. , et al. , 2006. Discussion on validation of Nakamura's technique. Earthquake Engineering and Engineering Vibration, 26(1): 43—48. (in Chinese) doi: 10.3969/j.issn.1000-1301.2006.01.007 荣棉水, 李小军, 王振明等, 2016. HVSR方法用于地震作用下场地效应分析的适用性研究. 地球物理学报, 59(8): 2878—2891 doi: 10.6038/cjg20160814Rong M. S. , Li X. J. , Wang Z. M. , et al. , 2016. Applicability of HVSR in analysis of site-effects caused by earthquakes. Chinese Journal of Geophysics, 59(8): 2878—2891. (in Chinese) doi: 10.6038/cjg20160814 姚鑫鑫, 任叶飞, 温瑞智等, 2019. 强震动记录H/V谱比法计算处理的若干关键环节. 震灾防御技术, 14(4): 719—730 doi: 10.11899/zzfy20190403Yao X. X. , Ren Y. F. , Wen R. Z. , et al. , 2019. Some technical notes on the data processing of the spectral ratio based on the strong-motion records. Technology for Earthquake Disaster Prevention, 14(4): 719—730. (in Chinese) doi: 10.11899/zzfy20190403 张照鹏, 荣棉水, 卢滔等, 2019. HVSR谱比法应用于强震数据时基本假定合理性讨论. 地震研究, 42(4): 538—545 doi: 10.3969/j.issn.1000-0666.2019.04.012Zhang Z. P. , Rong M. S. , Lu T. , et al. , 2019. Discussion on rationality of basic assumptions when applying HVSR method to strong earthquake data. Journal of Seismological Research, 42(4): 538—545. (in Chinese) doi: 10.3969/j.issn.1000-0666.2019.04.012 D'Amico V. , Picozzi M. , Baliva F. , et al. , 2008. Ambient noise measurements for preliminary site-effects characterization in the urban area of Florence, Italy. Bulletin of the Seismological Society of America, 98(3): 1373—1388. doi: 10.1785/0120070231 Del Gaudioa V. , Zhao B. , Luo Y. H. , et al. , 2019. Seismic response of steep slopes inferred from ambient noise and accelerometer recordings: the case of Dadu River valley, China. Engineering Geology, 259: 105197. doi: 10.1016/j.enggeo.2019.105197 Konno K. , Ohmachi T. , 1998. Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bulletin of the Seismological Society of America, 88(1): 228—241. doi: 10.1785/BSSA0880010228 Lee C. T. , Cheng C. T. , Liao C. W. , et al. , 2001. Site classification of Taiwan free-field strong-motion stations. Bulletin of the Seismological Society of America, 91(5): 1283—1297. Li X. J. , Li N. , Rong M. S. , et al. , 2022. Novel evaluation method for site effect on earthquake ground motion based on modified horizontal to vertical spectral ratio. Frontiers in Earth Science, 10: 938514. doi: 10.3389/feart.2022.938514 Nakamura Y. , 1989. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Quarterly Reports of RTRI, 30(1): 25—33. Rong M. S. , Wang Z. M. , Woolery E. W. , et al. , 2016. Nonlinear site response from the strong ground-motion recordings in western China. Soil Dynamics and Earthquake Engineering, 82: 99—110. doi: 10.1016/j.soildyn.2015.12.001 Wang Z. M. , Carpenter N. S. , Woolery E. W. , 2019. Horizontal-to-vertical spectral ratio of S-waves and SH-wave transfer functions at the vertical seismic and strong-motion arrays in the Central United States. Journal of Applied Geophysics, 162: 64—71. doi: 10.1016/j.jappgeo.2018.10.017 Wen R. Z. , Ren Y. F. , Shi D. C. , 2011. Improved HVSR site classification method for free-field strong motion stations validated with Wenchuan aftershock recordings. Earthquake Engineering and Engineering Vibration, 10(3): 325—337. doi: 10.1007/s11803-011-0069-x 期刊类型引用(0)
其他类型引用(2)
-