Application of the UVA Oblique Photograghy Technique in the Field of Ludian Earthquake, Yunnan
-
摘要: 本文首先介绍了利用六旋翼无人机倾斜摄影系统和动力三角翼倾斜摄影系统在云南鲁甸地震现场开展倾斜摄影数据采集、三维建模的相关工作，利用建立的三维模型，分别对灾区房屋震害、滑坡、滚石、堰塞湖地震地质灾害进行了系统的分析研究。相关研究表明：倾斜摄影技术可以表现地震灾害场景和具体灾害特征，对遥感在地震灾情精细化分析及了解灾区建筑物的布局、破坏程度和恢复重建建筑现状方面具有重要的意义。该项工作是地震领域初次利用倾斜摄影技术开展的相关研究，对于进一步深入研究倾斜摄影新技术在地震灾害领域和其他相关领域的工作具有很好的借鉴作用。Abstract: In this paper we describe the use of six rotor UVA oblique photography system and power of delta wing oblique photography system, and carry out related work oblique photography data acquisition, 3D modeling of the Ludian earthquake in Yunnan. Based on the 3D model, earthquake disaster housing and seismic geological disasters of landslides, rockfall and earthquake-induced lake are analyzed. Our study has a good reference for further research of new technology in earthquake disaster tilt photography and other related works.
-
Key words:
- Yunnan Ludian earthquake /
- Oblique photography /
- Earthquake disaster
-
引言
据河北省地震台网测定，2018年2月12日18时31分在河北省廊坊市永清县(39.37°N，116.67°E)发生4.3级地震，震源深度20km。此次地震未造成人员伤亡，在震中区域造成个别房屋墙皮脱落，河北省廊坊市、固安市震感较强，北京市、天津市普遍有感。地震发生后，中国地震局工程力学研究所对永清4.3级地震强震动记录进行了远程和现场记录回收。强震动记录是地震能量在地表或建筑结构上的真实反映，是工程地震、地震学研究的基础数据资料(朱升初等，2018)，本文对此次地震构造背景及强震动记录进行了分析，对永清4.3级地震烈度采用建立虚拟台站方法进行了校正，为河北地区烈度调查、地震区划等提供基础资料。
1. 发震构造背景
永清4.3级地震区域构造主要以NE、NW向断裂为主，地震发生在廊固凹陷内(图 1)，廊固凹陷是冀中坳陷北部的一个古近系湖盆，整个凹陷大体呈NE走向，南北长约90km，东西宽20—40km，面积约2600km2。该凹陷是受区域伸展并伴有走滑作用形成的典型箕状凹陷，整个凹陷表现为北断南超、西断东超的构造格局。由于受NE向断裂构造控制，表现出整体不均匀下降的箕状断陷特征。该凹陷内曾于1621年发生廊坊永清$ 5\frac{1}{2} $级地震(国家地震局震害防御司，1995)。
图 1 永清4.3级地震构造背景Figure 1. The seismotectonic background map of the Yongqing MS4.3 earthquake(1)施庄断裂(2)南口山前断裂(3)南口-孙河断裂(4)小汤山-东北旺断裂(5)黄庄-高丽营断裂 (6)顺义-前门-良乡断裂(7)南苑-通县断裂(8)新夏垫断裂(9)夏垫断裂(10)大兴凸起东缘断裂 (11)蓟县山前断裂(12)宝坻断裂(13)桐柏断裂(14)河西务断裂(15)天津断裂(16)大城东断裂 (17)沧东断裂(18)海河断裂(19)牛驼镇凸起东缘断裂(20)徐水断裂(21)东垒子-涞水断裂 (22)永定河断裂(23)徐水南断裂本次地震发生在廊固凹陷东侧控制断裂河西务断裂带附近，该断裂北起河西务西北，向西南经大王务，延伸至永清东码头镇，全长近50km，走向NE，倾向SE，为上陡下缓的铲形正断层，断层产状与震源机制参数基本一致(王晓山等，2018)，所以永清4.3级地震发震构造可能为河西务断裂。周月玲等(2018)通过浅层地震勘探和钻孔联合剖面探测，结合年代样品测试，对河西务断裂活动性进行了综合研究，揭示断裂上断点埋深约150m或以浅，第四系底界面垂直断错为20—45m，断裂最新活动时代为晚更新世早期。
2. 观测台站及记录获取
2.1 强震动台站分布
本次地震强震动数据主要由中国地震局工程力学研究所提供，共获取3个分量自由场记录74组，常规处理后得到222条3分向加速度记录，获取记录的强震动台站均为土层台，记录台站分布如图 2所示，由图 2可知，触发台站主要分布在震中北部，南部无台站触发，这可能是因为首都圈台站较密集，而河北省境内强震动台站稀疏，或与仪器型号不同且很多台站老旧未改造有关。触发台站的震中距为29.4—185.5km，其中震中距＜50km的记录有5组，震中距为50—100km的有29组，震中距为100—150km的有37组，震中距为150—200km的有3组，距震中最近的台站为长子营台(11ZZY)，震中距为29.4km。虽然此次地震获取的强震动记录较全面，但由于震级较小，未造成大的灾害损失，且震中附近台站相对稀疏，能供工程使用的强震动记录较少。
图 2 永清4.3级地震触发台站分布图Figure 2. The distribution of seismic trigger station of Yongqing MS4.3 earthquake2.2 强震动记录获取
在强震动记录使用分析前，首先对波形数据进行消除基线漂移，即原始强震动记录通过减去事前记录平均值的方法进行零线调整(吕俊强等，2015)，通过分析计算得到震中距R＜50km的5个台站强震动记录基本信息及相关参数，见表 1。其中距震中最近的长子营台记录到的最大峰值加速度为7.9cm/s2，觅子店台记录到的东西向最大峰值加速度为36.1cm/s2、南北向最大峰值加速度为40.4cm/s2，红星中学台记录到的垂直方向最大峰值加速度为10.8cm/s2，5个台站加速度时程曲线如图 3。由表 1可知，离震中较近的长子营台(震中距29.4km)较离震中较远的觅子店台(震中距38.4km)和于家务台(震中距36.8km)各分量加速度小，影响地震动强度的因素除震中距外，还包括仪器型号、场地特征条件和地质构造、场地土层结构、覆盖层厚度、地形及断裂分布等(王海云等，2010；温瑞智等，2013；郭秋娜等，2018)。本文根据温瑞智等(2015)提出的适用于我国强震动台站场地划分经验曲线，采用H/V谱比法对本次地震触发台站进行分类，分类结果见图 4。由图 4可知，记录到的本次地震台站大部分为Ⅰ类场地(51.35%)，其次为Ⅱ类场地(35.14%)和Ⅲ类场地(13.51%)，其中，觅子店台和红星中学台为Ⅰ类场地，于家务台和张家湾中学台为Ⅱ类场地，长子营台为Ⅲ类场地。长子营台与觅子店台和于家务台记录差异较大，这可能与场地类型有关。
表 1 强震动记录(震中距＜50km)及相关参数Table 1. Strong motion records(epicenter distance < 50km) and related parameters台站名称 台站代码 仪器型号 震中距/km PGA/cm·s-2 EW SN UD 觅子店 11MZD REFTEK/SLJ-100 38.4 36.1 40.4 8.5 红星中学 11YHZ REFTEK/SLJ-100 43.4 -5.4 -4.6 10.8 于家务 11YJW ETNA/ES-T 36.8 -21.7 -18.5 -8.1 张家湾中学 11ZJW REFTEK/SLJ-100 47.9 -5.8 4.2 -4.4 长子营 11ZZY REFTEK/SLJ-100 29.4 5.0 -7.9 -3.5 3. 谱分析
加速度反应谱是目前国内外工程抗震设计的重要依据，反应谱形状、幅值及卓越周期的差异体现了场地条件、震级对地面运动反应谱的影响(温瑞智等，2013；徐钦等，2017)。为避免地震时建筑结构发生共振，确保结构安全，有必要对加速度反应谱成分进行分析。本文选取震中距＜50km的5个典型强震动台站进行频谱分析，阻尼比取0.05，得到本次地震EW、SN、UD向加速度反应谱(图 5)。
由图 5可知，觅子店台记录到了本次地震EW向和SN向加速度反应谱最大值，红星中学台记录到了本次地震SN向加速度反应谱最大值。觅子店台与长子营台在EW和SN向差别最大，觅子店台与张家湾中学台在SN向差别最大。5个典型台站记录到的加速度反应谱高频成分丰富，谱卓越周期主要为0.1—0.3s，0.3s后迅速下降。根据现场调查，震区主要房屋类型为土木结构、木房顶砌体结构、预制板砌体结构、钢结构简易房屋4种，本次地震除少数土木结构房屋原有裂缝扩大、墙皮开裂外，其他类型房屋未发现新开裂、屋顶掉瓦、墙皮掉灰现象，中小城市主要建筑物结构自振周期为0.3—1s(王文才等，2018)，因此本次地震对该震区自振周期范围内的建筑物破坏影响较小。通过对比不同震中距台站各方向反应谱特征周期发现，震中距对特征周期的影响不明显。
傅氏谱能反映地震波在频域上的幅值特性，地震动记录中包含多种振动频率，不同地区的主要振动频率不同，结构自身具有不同频率、幅频响应，因而在不同频率振动作用下，结构振动不同。为避免地震作用下结构发生共振现象，确保结构安全，有必要对地震频率成分进行分析(冉志杰等，2012)。
本文计算了永清4.3级地震震中距50km以内的5个典型台站强震动记录傅氏谱(图 6)，由图 6可知典型台站傅氏谱谱型主要为单峰值型，强震动记录EW分量傅氏谱最大峰值为3.5—8.5Hz，SN分量傅氏谱最大峰值为2.9—6.2Hz，UD分量傅氏谱最大峰值为2.1—9.8Hz。总体来看，UD分量频率范围分布较宽，EW分量傅氏谱最大峰值大于SN分量和UD分量傅氏谱峰值。一般而言，随着震中距的增大，地震动记录长周期(低频)越显著，但本次地震5个典型台站特征不明显，可能与仪器型号及场地特征条件有关。另外，EW分量傅氏谱最大峰值随着震中距的增大呈先减小后增大的趋势，而SN分量和UD分量无该特征。
本文以东西向为正方向，统计峰值加速度反应谱与永清4.3级地震触发台站同震中所处方位角和震中距的关系(图 7)。由图 7可知，峰值加速度反应谱在同一震中距下整体随着方位角的增大而增大，达70°左右后下降，这可能与附近河西务断裂构造走向有关；峰值加速度反应谱值随着震中距的增加而减小。
图 7 峰值加速度反应谱与震中距、方位角关系图Figure 7. The relation among the peak ground acceleration and epicentral distance、azimuth angles4. 地震烈度分析
地震烈度是衡量地震破坏和地震动强弱程度的重要指标。Wald等(1999)认为当烈度较小时，用峰值加速度计算烈度较合适，因为人对加速度较敏感。随着强震动观测技术的发展，由于以各观测点峰值加速度作为依据可真实反映地震发生时各点地面运动强度(田秀丰等，2013)，基于地震动记录快速确定地震仪器烈度，越来越受到研究者的重视(崔建文等，2016)。目前国内外基于强震动记录资料的地震烈度计算方法较多，不同方法利用不同的地震动参数，包括时间、频度和强度等，建立地震烈度与地震动参数之间的关系。本文采用美国ShakeMap烈度计算方法(李俊等，2010)，计算永清4.3级地震事件记录台站烈度，并采用克里格插值方法进行网格化插值，绘制等值线图(图 8)。Wald等(1999)的研究认为，当烈度低于Ⅴ度时峰值加速度PGA与烈度相关性明显高于峰值速度PGV，而当烈度为Ⅴ—Ⅶ度时，PGA和PGV与烈度相关程度相当，考虑本次地震烈度较低，故采用峰值加速度计算烈度。由图 9可知，本次地震仪器烈度等震线长轴沿NW走向，且极震区位于震中东北部，而不是震中附近，但由现场宏观烈度调查(图 8)与震源机制解(王晓山等，2018)可知该地震等震线长轴沿NE向展布，这种差异现象的出现可能与台站分布不均匀有关，震中南部无强震动记录，所以图 8中震中北部附近出现数据空白现象，不利于获取实际烈度长轴方向真实性。
图 8 永清4.3级地震校正前仪器烈度分布图Figure 8. The distribution of the un-corrected instrumental seismic intensity of Yongqing MS4.3 earthquake可靠快速生成地震烈度可为人员伤亡估计、应急救援决策和工程抢险修复决策提供重要依据，为解决由于缺少地震动观测值出现空白区域现象及烈度长轴方向和实际方向不一致的问题，本文首先利用永清4.3级地震中获取的所有台站加速度峰值，按地震动峰值衰减模型(王国新，2001)确定该地区地震动衰减关系，然后在整个地震区内建立空间随机假设台站进行补点插值，得到新的地震烈度分布图。地震动峰值衰减模型(王国新，2001)如下：
$$ \lg Y = {{\rm{C}}_1} + {{\rm{C}}_2}M + {{\rm{C}}_3}M\lg \left({R + {R_0}} \right) + {{\rm{C}}_4}\lg \left({R + {R_0}} \right) \pm \varepsilon $$ (1) 式中，Y为PGA(峰值加速度)，M为震级，R为震中距，C1、C2、C3、C4为待求常数，R0为近场饱和因子，ε为统计分析误差项。采用通用全局优化法拟合得到地震动峰值加速度随震级和震中距的变化规律，考虑本次地震震级较小，故选取震中距＜100km的Ⅰ类场地强震动数据进行拟合，其中水平向按EW、NS向加速度峰值矢量拟合回归分析，拟合曲线如图 9，计算得到的衰减关系式为：
$$ \lg Y = 1.6882 + 0.3926M - 0.1546M\lg \left({R + 5} \right) + 0.7000\lg \left({R + 5} \right) \pm 0.2287 $$ (2) 
图 9 永清4.3级地震现场调查烈度分布图Figure 9. The distribution of site investigation intensity of Yongqing MS4.3 earthquake峰值加速度随震中距的变化如图 10所示，由图 10可知峰值加速度随着震中距的增加而减小。同一震中距不同台站测得的峰值加速度离散性较大。通过拟合各台站实际PGA值，与俞言祥等(2013)地震动衰减关系进行对比，可知20km范围内2条曲线吻合较好，而20km外拟合曲线呈衰减相对缓慢状态，这可能与本次地震震级较小、记录台站均为土层台、远场有一定放大作用有关(王海云，2011)。由观测值得到的预测方程虽存在一定离散性，但整体趋势与实际结果保持一致。
图 10 峰值加速度随震中距的变化Figure 10. `The change of peak ground acceleration with epicentral distance在震区空白区域随机产生N个空间坐标的假设台站，按式(2)计算水平峰值加速度，按美国ShakeMap烈度计算方法转换成仪器烈度，绘制增加假设台站的烈度分布图(图 11)，由图 11可知，增加假设台站的仪器烈度分布图与现场烈度调查分布图极震区方向基本一致，进一步验证了河西务断裂为本次地震发震构造的可能，解决了台站分布不均匀导致地震烈度影响场计算出现缺值现象的问题，为以后缺少台站记录震区提供准确快速制作烈度分布图的思路，为震害调查和地震应急救援提供重要依据。
图 11 永清4.3级地震校正后烈度分布图Figure 11. The distribution of the corrected instrumental seismic intensity of Yongqing MS4.3 earthquake5. 结论
本文介绍了河北永清4.3级地震发震构造背景，并对国家强震动台网中心在本次地震中获取的加速度记录特征进行了分析，为河北省地震动特征和震害调查提供了一定参考，结论如下：
(1) 河北永清4.3级地震在廊固凹陷东侧控制断裂河西务断裂带附近，地震烈度极震区走向与该断裂展布方向一致，结合震源机制解初步判定本次地震发震构造为河西务断裂。
(2) 74个强震动台站获取到强震动加速度记录，震中距位于29.4—185.5km，其中觅子店台记录到的本次地震东西向最大峰值加速度为36.1cm/s2，垂直方向最大加速度为40.4cm/s2，于家务台记录到的南北向最大峰值加速度为18.5cm/s2。采用H/V谱比法对触发台站进行分类，可知Ⅰ类场地居多。
(3) 通过对5个典型台站反应谱分析发现永清4.3级地震加速度反应谱高频成分丰富，谱卓越周期主要为0.1—0.3s。傅氏谱谱型主要为单峰值型，UD分量频率范围分布较宽，EW分量傅氏谱最大峰值大于SN分量和UD分量。
(4) 通过分析峰值加速度反应谱与震中距、方位角的关系，发现峰值加速度反应谱整体随着方位角的增大而增大，达70°左右后下降，这可能与附近河西务断裂构造走向有关，有待进一步研究。其次，峰值加速度反应谱值整体随着震中距的增加而减小。
(5) 由本次地震仪器烈度分布图可知，台站分布不均匀性直接影响烈度分布的真实性，为避免此类问题，拟合出该地区地震动衰减关系，随机产生空间假设台站进行弥补，为以后缺少台站记录震区提供准确快速制作烈度分布图的思路，为震害调查和地震应急救援提供重要依据。由于收集到的强震动数据有限，文中仅采用本次地震获取的峰值加速度拟合当地地震动衰减关系，拟合结果可能存在误差，今后将收集更多的强震动数据进行拟合，以校正该地区地震动衰减关系。
-
图 1 中心投影与倾斜投影原理对比图
Figure 1. Comparison of principles between central projection and oblique projection
图 6 从另一角度展示滚石击中房屋，滑坡土方掩埋房屋
Figure 6. The building is destroyed by the landslide and rockfall from different view angle
图 7 滑坡掩埋房屋、滚石击中房屋（放大）
Figure 7. The building is destroyed by the landslide and rockfall(enlarged partial detail)
图 11 动力三角翼倾斜摄影的红石岩堰塞湖及其周边影像
Figure 11. The oblique photography of Hongshiyan earthquake-induced lake and its surrounding imagery by the power delta wing aircraft
表 1 无人机基本参数
Table 1. The basic parameters of UVA
飞行器性能 技术参数 总重量 ＜8kg（含两块电池） 飞行升限 1500m 载重与续航时间 载重2.5kg，续航时间≥30min 悬停精度 垂直方向±1m，水平方向±2m 最大倾斜角度 35° 最大升降速度、平飞速度 最大升降速度6m/s；最大平飞速度10m/s 巡航速度 ≥6m/s 平均维修时间 ≤30min 工作效率 0.5km2/架次 
下载: 导出CSV
表 2 云台和相机基本参数
Table 2. The basic parameters of the cloud platform and the camera
云台性能 技术参数 尺寸 210mm×265mm 重量 2.1kg CCD数量 5个 CCD尺寸 13.2mm×8.8mm 最小曝光间隔 2s 总像素 大于1亿 防水等级 IP67 曝光方式 定点曝光 POS 记录定点曝光POS信息 影像分辨率 2—10cm 定位精度 1m
下载: 导出CSV
表 3 动力三角翼基本参数
Table 3. The basic parameters of Delta-wing airplane
性能 指标 基本参数 巡航速度120km/h，最大速度150km/h
耗油12L/h，油箱容量50L
巡航航程702km
最大空机重量260kg，最大起飞重量450kg
翼展9.97m
起飞距离：小于60—100m适用环境 抗风性能：4级 保护装置 具有安全弹射伞防护功能 作业参数 作业高度：200—1500m
作业效率：4—6km2/h
成果精度：3—14cm
下载: 导出CSV
表 4 动力三角翼搭载云台和相机参数
Table 4. The Basic Parameters of the cloud platform and the camera on Delta-wing airplane
基本参数 总重≤25kg
平台尺寸：85mm×50mm×25mm
独立供电
统一的数据存储和传输接口曝光方式 相机定时曝光 摄影精度 0.03-0.20m 适用环境温度 -10℃—50℃ CCD大小 24mm×35mm 像素大小 4000万
下载: 导出CSV
-
期刊类型引用(1)
1. 李姜，李铁纯，张兴，郭彦宇. 农村人口流动对地震灾害人员伤亡评估的影响研究-以张家口为例. 震灾防御技术. 2022(03): 557-568 . 
本站查看其他类型引用(0)
-
下载:
点击查看大图
计量
- 文章访问数: 104
- HTML全文浏览量: 47
- PDF下载量: 6
- 被引次数: 1
