Algorithm of Accurate Location of Buried Personnel in Special Environment after Earthquake
-
摘要: 针对目前对震后压埋人员定位精度较低、探测设备成本高且易受环境影响等不足,提出适用于压埋环境特性的压埋人员手机WiFi定位方法,通过衰减因子模型对WiFi探针获取的RSSI数据进行距离解算,结合简化压埋环境内部信号传输方式,采用高斯-卡尔曼滤波对获取的RSSI数据进行处理,通过模型测定的距离,利用改进附有参数的加权最小二乘平差方法,结合粒子群优化算法,最终得到压埋人员手机平面坐标位置。研究结果表明,该方法具有较高精度,在10 m×10 m范围内其平面坐标定位误差在0.3 m左右,可为震后压埋人员应急救援提供辅助决策。Abstract: In view of the current shortcomings of the positioning accuracy of the buried people after the earthquake, the high cost of the detection equipment and the vulnerability to environmental impact, in this paper we propose a WiFi Positioning Method to search for the buried people's mobile phones, which is suitable for the characteristics of the signals from buried environment. The RSSI data obtained by the WiFi probe is calculated by the attenuation factor model. Combined with the simplified internal signal transmission mode of the buried environment, the obtained RSSI data is processed by Gauss Kalman filter. Through the distance measured by the model, using the improved weighted least square adjustment method with parameters, combined with particle swarm optimization algorithm, it is capab le to obtain the plane coordinate position of the buried personnel's mobile phone.The experimental results show that this method has high accuracy, with the average error of plane coordinate positioning about 0.3 m within the range of 10 m × 10 m, which can provide auxiliary decision-making for emergency rescue of buried personnel after earthquake.
-
引言
大量震例表明,活断层不仅是产生地震的根源,而且当地震发生时断层沿线的建筑物破坏、人员伤亡最为严重,其损失明显大于断层两侧的其它区域(邓起东等,2007;刘保金等,2008;酆少英等,2010)。因此,查明地下断层的性质、活动性及其空间展布特征,对城市的防震减灾工作,在规划建设中有效避让地震活断层,减轻活断层的破坏作用等方面具有重要的意义。城市周边的活断层一般被第四纪松散沉积物所覆盖,落差较小、埋深较浅,仅依靠地表地质地貌工作已不能对断层定位和活动性进行全面研究。目前,浅层地震反射波法是活断层探测的主要方法之一,通过反射地震剖面可以直观地判定断层的存在和形态,该方法在大多数城市活断层探测中均取得了良好的效果(刘保金等,2007;赵成彬等,2011)。
汤东断裂是汤阴地堑的东界断裂,对汤阴地堑的生成、演化和发展具有重要的控制作用。很多学者对汤阴地堑和汤东断裂作了大量的研究工作,韩慕康等(1980)通过对汤阴地堑及其周边进行地震地质调查,发现汤阴地堑只有南部是地堑,而北部是由一条东侧断裂(汤东断裂)构成的半地堑;杨承先(1984)利用石油、煤田地震反射剖面并结合部分钻井资料,提出汤阴地堑形成于新近纪,其北部主要受东缘的汤东断裂控制,南部受东西边缘的汤东、汤西断裂控制,且汤东、汤西断裂新构造活动明显;任青芳等(1998)利用深地震测深资料获得了汤阴地堑及其邻区的二维速度结构剖面,认为汤阴地堑的基底断裂发育,壳幔结构复杂,深部界面与速度等值线起伏变化大,且壳内有低速层,汤阴地堑及邻区存在发生中强地震的深部构造背景;刘保金等(2012)采用深、浅地震反射探测方法,获得了汤阴地堑及其邻区的地壳精细结构和断裂的深、浅构造关系。这些研究成果对了解汤阴地堑的形成和演化、壳幔结构特征以及断裂的深、浅构造关系提供了重要的基础资料。本文在分析该区已有深、浅地震探测剖面资料的基础上,跨汤东断裂开展了高分辨率的浅层地震剖面探测,取得了高信噪比的浅层地震反射叠加剖面图像,分析讨论了汤东断裂的近地表特征、性质及其活动性。研究结果为汤阴地堑的活动构造研究以及研究区的地震危险性评价提供了地震学依据。
1. 地质构造概况和地震剖面位置
汤阴地堑位于太行山与华北平原的过渡带上,是太行山隆起和内黄隆起间的一个北东向凹陷,东西分别受汤东断裂和汤西断裂的控制,南北分别以新乡-封丘断裂和安阳断裂为界(刘保金等,2012)(图 1)。汤阴地堑整体呈北北东走向,南北长100km,东西宽15km,面积1500km2,属新生代的地堑,在新生代以强烈下沉为主要构造活动,整个地堑内被巨厚的新生代地层所覆盖,其东西两侧汤西断裂和汤东断裂的发育状况不尽相同,形成南部呈断陷型地堑、北部呈半断陷型地堑的不规则地堑(梁生正等,2003;高战武等,2014;郑建彬,2015)。
汤东断裂是汤阴地堑东界的控制断裂,断裂全长约90km,走向北北东,倾向西,为上陡下缓的正断层。断裂以东的内黄隆起为太古界至古生界组成的向东倾伏的隆起,地表被新近系和第四系所覆盖,局部地段古生界基岩残山直接出露;断裂西盘的汤阴地堑则沉积了巨厚的古近系、新近系、第四系。断裂北段控制了汤阴次级凹陷的发育,断层附近的古近系、新近系和第四系最大厚度分别为3100m、1400m和40m左右,而在断裂下盘的内黄隆起上,第四系很薄,有的地方新近系出露。断裂发育在古生代及其以前的地层中,向上切割了古生代—新生代地层,第四纪以来仍然制约着断裂两侧的地层沉积,该断裂是第四纪活动断层。
在汤阴地堑已开展过一条长约45km的深地震反射探测剖面(图 1),剖面揭示了汤阴地堑及其周边的地壳深浅构造特征(刘保金等,2012)。根据深地震反射探测剖面,把该区域地壳自上而下分为上地壳、中地壳、下地壳和壳幔过渡带4个部分。上地壳又可以细分为沉积盖层和结晶基底,沉积盖层主要是新生代、中生代和古生代沉积岩系,在地震剖面上具有良好的反射性质;结晶基底可能由元古代、太古代的变质岩组成,在地震剖面探测中表现为一些反射能量较弱、且无规律可循的短小反射。深地震反射资料揭示汤东断裂是由2条具有不同构造特征的断层组成,2条断层倾向为西,东侧断层倾角较陡,向上错断了新近纪底界面后归并到西侧断层上,向下延伸至上地壳底部;西侧断层在剖面上具有铲型正断层的特征,断层向上错断了新近纪内部地层,向下错断了所有的沉积地层,延伸至7—8km的结晶基底附近。深地震反射资料较好地揭示了汤东断裂的中深部构造特征,但不能很好地反映断裂在第四纪内部延展特征,为了获得汤东断裂在第四纪内部的构造特征,在垂直断裂的走向进行了高分辨浅层地震反射探测,得到了汤东断裂清晰的近地表构造特征。
已有资料表明汤东断裂有着明显的第四纪活动特征,为了查明其位置、性质及浅部构造特征,在汤阴县东部跨汤东断裂进行了高分辨的浅层地震探测。沿着小屯村北的汤屯线自西向东布设小屯村测线,起点位于距第四街工业园70m处,终点位于距汤屯线5km里程碑东约500m处,测线全长2438m。古贤镇测线沿着省道S302自西向东布设,起点位于古贤镇北的加油站东侧,终点位于省道100.5km里程碑附近,测线全长3278m。
2. 地震数据采集和资料处理
为了压制在浅层地震数据采集时不可避免的车辆、行人及邻近居民点的干扰,地震波激发采用美国产M612-18型可控震源,数据接收采用德国SUMMIT遥测数字地震仪,采样间隔0.5ms、记录长度2s。根据现场试验确定的观测系统参数为:2m道间距、最小偏移距20m、200道接收、20次覆盖。
数据处理采用资料Focus反射地震数据处理软件包,根据本次地震勘探的工作特点和资料特征,采用了几何扩散补偿、叠前去噪、静校正、动校正、剩余静校正、CMP叠加、叠后去噪等处理模块,建立了图 2所示的主要数据处理流程。
该工作区地震资料上面波、声波以及各种随机干扰波较为严重,因此压制地震记录中各种高、低频干扰,提高地震剖面的信噪比和分辨率是本次资料处理的重点。通常根据噪音在地震记录上的不同特征,需采用不同的滤波方法来压制噪音提高资料的信噪比(张军华等,2006;何银娟等,2015)。本次数据处理主要采用了如下去噪方法:
(1)一维带通滤波和二维f-k滤波。根据原始单炮记录上面波、声波等干扰波能量较强的特点,资料处理时采用带通滤波方法对面波及其他随机低频干扰波进行了压制,采用二维f-k滤波对声波等线性干扰波进行滤除。对比图 3(a)、(b)2张单炮纪录可以看到,经带通和f-k滤波后,面波、声波、随机干扰波明显变弱或者消失,反射波能量变强,反射波同相轴清晰可见,单炮纪录的信噪比得到了较大的提高。
(2)矢量分解压噪。对于干扰信号和有信号频率相近的单炮记录,根据干扰信号偏离有效信号所形成的夹角采用曲边时窗矢量分解压噪的方法,可达到既去除干扰信号又很好地保护有效信号的目的(何银娟等,2015)。从图 3(b)、(c)2张单炮纪录可以看出,经过矢量分解压噪后,面波、声波基本消除,随机干扰波得到了进一步的压制,有效波同相轴变得更加连续(图 3中蓝圈处),资料的信噪比进一步提高。
(3)倾角滤波和f-x域噪声衰减。虽然叠前采用了各种压制干扰波的方法,但经过CMP叠加后剖面上仍然存在一些干扰波。为了消除叠后的线性干扰波和随机干扰波,本次叠后数据处理主要采用倾角滤波和f-x域随机噪声衰减技术,达到消除干扰信号、增强有效信号的目的。
3. 地震反射剖面特征
3.1 小屯村测线浅层地震反射剖面特征
图 4所示的小屯村测线地震反射叠加时间剖面具有较高的信噪比和分辨率,剖面揭示的多组浅部地层反射波均能被连续可靠追踪,且界面起伏变化形态也非常清楚。根据剖面的反射波组特征,共识别出11组能量较强、横向可连续追踪的反射震相(TQ、T1—T9、TN)。依据剖面反射波震相特征和本区域地质资料,把剖面上反射波TQ解释为第四纪地层的底界面反射,T1—T9解释为新近纪内部地层的界面反射,TN解释为新近纪地层的底界面反射。从剖面横向上看,反射波TQ自东向西略微向下倾伏,在剖面上表现出强反射震相特征。新近纪地层界面反射波以桩号1540m为界,东西两侧反射震相特征差异明显,桩号1540m以东反射震相较少,新近系沉积较薄,新近系底界面埋深较浅,可能是由于新近纪时期内黄隆起区强烈的构造活动,导致了地层的抬升与剥蚀;而桩号1540m以西反射震相丰富,可识别出多组新近纪地层界面反射波(T1—T9、TN),新近纪地层分层性好、反射波能量较强、同相轴连续性好,说明新近纪时期汤阴地堑在该区域的沉积环境较为稳定。
剖面所揭示的断裂结构和构造特征非常清楚。在剖面桩号1540m附近,从剖面上可以看到一个由浅到深的反射特征分界线,其东西两侧地层产状和反射波组特征明显不同,根据这些特征不难判定是断裂在地震剖面上的反映。在断裂下降盘一侧,地层界面反射波丰富,第四系底界面反射波TQ能量较强,同相轴分段连续性好,新近纪地层自东向西倾伏,其底界面TN埋深较深,约800m;断层上升盘一侧,地层基本呈水平展布,第四系底界面反射波TQ能量较弱,新近纪地层界面反射波明显少于下降盘一侧,其底界埋深较浅,约480m。根据区域地质资料可知该断裂为汤东断裂,断裂向上错断了第四系底界面TQ,上断点埋深约50m,是第四纪隐伏活动断裂。在汤东断裂FP3西侧还可以看到与其倾向相反的正断层FP1、FP2,断层FP1、FP2是汤东断裂的分支断层,分别在TWT 600ms和TWT 400ms归并到主断层FP3上。
3.2 古贤镇测线浅层地震反射剖面特征
图 5给出了古贤镇测线的浅层地震反射波叠加时间剖面,剖面具有较高的信噪比和分辨率,断裂构造特征也非常清晰。对比图 4、图 5可以看出,2条地震反射叠加时间剖面上反射波震相特征和汤东断裂的结构和构造特征均表现出了较大的相似性,剖面西段(汤东断裂FP7以西)反射震相丰富、反射波同相轴连续性好、反射波能量较强,表现出沉积盆地连续沉积的特征;而在剖面东段(汤东断裂FP7以东)新近纪地层沉积较薄,其底界面埋深明显浅于西段,新近纪地层之下为能量较弱、不连续的短小反射,可能是元古代、太古代变质岩系。图 5中汤东断裂主断层FP7为西倾铲型正断层,断层错断了剖面上的所有沉积层,上断点埋深较浅,约为30m;在其西侧有3条反向次级正断层,与主断层呈“Y”型组合特征。在主断层的东侧还有1条次级断层FP8,其倾向为西,倾角较陡,向上错断了新近纪内部地层T2,上覆盖的第四纪地层界面反射波未发现扭曲、错断。整体上看,2条浅层地震反射叠加时间剖面上揭示的反射波震相东西段差异特征及汤东断裂形态与浅部构造特征均表现出较大的一致性。
图 5 古贤镇测线反射叠加时间剖面和深度解释剖面图(刘保金等,2012)Figure 5. Stacked section and depth profile of the shallow seismic reflection in the Guxian town line(Liu et al., 2012)4. 结论和讨论
本项研究通过采用高分辨率的浅层地震勘探方法和选择恰当的数据处理技术,获得了太行山山前汤东断裂的位置、性质及其活动性。汤东断裂为走向NNE、倾向NWW的正断层,是汤阴地堑的主控边界断裂。深地震反射资料显示汤东断裂是由2条具有不同构造特征的断层组成,2条断层倾向均为西,东侧断层倾角较陡,向上错断了新近纪底界面后归并到西侧断层上,向下延伸至上地壳底部;西侧断层在剖面上具有铲型正断层的特征,断层向上错断了新近纪内部地层,向下错断了所有的沉积地层,延伸至7—8km的结晶基底附近,说明汤阴地堑新生代沉积主要受该断层的控制。浅层地震反射剖面揭示汤东断裂西侧断层浅部是由2—3条断层组成的Y字形构造,断层向上错断了第四系的底界TQ,其上断点埋深约30—50m,结合该区域钻孔资料可知汤东断裂错断了更新世地层,是更新世以来的隐伏活动断裂。
根据浅层地震剖面特征并结合区域地震及地质资料可知,汤东断裂属于汤阴地堑与内黄隆起的边界断裂,断裂两侧地层沉积特征差异明显,断裂下降盘新近纪地层界面反射丰富,地层横向上表现为两端向中部倾伏的凹陷型构造特征;断裂上升盘地层基本呈水平展布,新近纪内部地层缺失T6—T9,其底界埋深较浅,约480m,揭示了新近纪时期断裂构造活动剧烈,内黄隆起抬升,地层受到了明显的风化剥蚀,下降盘汤阴地堑沉降显著,新近纪地层沉积巨厚,其底界面埋深约800m。
浅层地震勘探方法是探测地层浅部结构、隐伏断裂浅部构造特征和对隐伏断裂进行准确定位的有效技术手段,但仅根据浅层地震剖面无法获知断裂在深部的结构情况及与其他断裂的交切关系,也无法确定地震剖面上反射层的地层属性和地层年代。因此,在浅层地震勘探解释时,还需结合深地震反射资料、区域地质、钻孔以及其他地震地质资料,才能得到断裂的深、浅构造特征,更好地理解断裂的形成和演化过程,并且结合地质、钻孔资料有助于进一步提高地震剖面解释的准确度。
致谢: 本文写作过程中得到了刘保金研究员的指导以及秦晶晶、邓小娟、刘增祺等同事的帮助,在此表示衷心的感谢。 -
表 1 模拟压埋环境下平面坐标定位结果(单位:米)
Table 1. Results of plane positioning in simulated buried environment(Unit:m)
点号 方案1 方案2 定位结果 测距中误差 定位误差 定位结果 测距中误差 定位误差 真实位置 P1 (4.86,1.79) 1.63 3.21 (5.18,4.86) 0.25 0.23 (5,5) P2 (5.21,0.58) 2.3 3.51 (6.14,4.21) 0.32 0.25 (6,4) P3 (5.42,4.37) 2.8 1.73 (6.28,5.84) 0.54 0.32 (6,6) P4 (3.80,3.42) 2.5 2.59 (4.51,5.86) 0.43 0.53 (4,6) -
柴伦尚,2014. 基于灾害现场定位搜救技术的无线信道研究. 杭州:浙江大学.Chai L. S., 2014. Research of radio channel based on a localization rescue system in disasters. Hangzhou:Zhejiang University. (in Chinese) 陈国良,张言哲,汪云甲等. 2015. WiFi-PDR室内组合定位的无迹卡尔曼滤波算法. 测绘学报, 44 (12):1314−1321Chen G. L., Zhang Y. Z., Wang Y. J., et al. 2015. Unscented Kalman filter algorithm for WiFi-PDR integrated indoor positioning. Acta Geodaetica et Cartographica Sinica, 44 (12): 1314−1321. (in Chinese) 陈群,高星伟,郭杭等. 2012. ZigBee网络中RSSI的测距模型分析. 测绘科学, 37(5):38−39, 75.Chen Q., Gao X. W., Guo H., et al. 2012. Analysis of RSSI ranging model at ZigBee network. Science of Surveying and Mapping, 37(5): 38−39, 75. (in Chinese) 陈雪纯,张入千,王立新. 2018. 基于Wi-Fi信号强度衰减模型的地震埋压人员定位方法研究. 华南地震, 38(4):10−16.Chen X. C., Zhang R. Q., Wang L. X. 2018. Research on positioning method of earthquake buried personnel based on signal strength attenuation model of wireless network. South China Journal of Seismology, 38(4): 10−16. (in Chinese) 何成文,袁运斌,潭冰峰. 2021. 一种基于UWB TDOA定位模式的迭代最小二乘算法. 大地测量与地球动力学, 41(8):806−809.He C. W., Yuan Y. B., Tan B. F. 2021. An iterative least squares algorithm based on UWB TDOA positioning model. Journal of Geodesy and Geodynamics, 41(8): 806−809. (in Chinese) 胡洋,田忠,伍习光等. 2018. 基于RSSI的无源信号定位算法研究. 自动化仪表,39(6):31−34.Hu Y., Tian Z., Wu X. G., et al. 2018. Study on the passive signal positioning algorithm based on RSSI. Process Automation Instrumentation, 39(6): 31−34. (in Chinese) 景裕,曹育森,朱明明等. 2021. 非接触式生命探测技术研究现状与发展. 中国医疗设备,36(1):1−4.Jing Y., Cao Y. S., Zhu M. M., et al. 2021. Research status and development of non-contact life detection technology. China Medical Devices, 36(1): 1−4. (in Chinese) 李惠芬,蒋向前,李柱. 2004. 高斯滤波稳健性能的研究与改进. 仪器仪表学报,25(5):633−637.Li H. F., Jiang X. Q., Li Z. 2004. Study and improvement of robust performance of Gaussian filtering. Chinese Journal of Scientific Instrument, 25(5): 633−637. (in Chinese) 李亦纲,张媛,李志伟. 2010. 地震现场倒塌建筑物的搜救策略研究. 震灾防御技术,5(4):477−483.Li Y. G., Zhang Y., Li Z. W. 2010. Research of search and rescue plan for collapsed buildings on earthquake disaster site. Technology for Earthquake Disaster Prevention, 5(4): 477−483. (in Chinese) 李桢,黄劲松,2016. 基于RSSI抗差滤波的WiFi定位. 武汉大学学报·信息科学版,41(3):361−366.Li Z., Huang J. S. 2016. WiFi positioning using robust filtering with RSSI. Geomatics and Information Science of Wuhan University, 41(3): 361−366. (in Chinese) 刘晶晶,宁宝坤,吕瑞瑞等,2017. 震后典型建筑物倒塌分类及救援特点分析. 震灾防御技术,12(1):220−229.Liu J. J., Ning B. K., Lv R. R., et al. 2017. Classification of typical building collapse and analysis of rescue characteristics after the earthquake. Technology for Earthquake Disaster Prevention, 12(1): 220−229. (in Chinese) 罗尚,肖东升,2020. 一种面向震后压埋人员的RSSI快速定位算法. 测绘科学,45(6):142−149.Luo S., Xiao D. S. 2020. An RSSI fast localization algorithm for buried personnel after earthquake. Science of Surveying and Mapping, 45(6): 142−149. (in Chinese) 裴凌,刘东辉,钱久超,2017. 室内定位技术与应用综述. 导航定位与授时,4(3):1−10.Pei L., Liu D. H., Qian J. C. 2017. A survey of indoor positioning technology and application. Navigation Positioning and Timing, 4(3): 1−10. (in Chinese) 武有文,洪利,姚振静等,2017. 地震废墟环境下Wi-Fi无线信号传播模型参数测量. 震灾防御技术,12(1):210−219.Wu Y. W., Hong L., Yao Z. J., et al. 2017. Parameter measurement of Wi-Fi radio signal propagating in constrction ruins by earthquake. Technology for Earthquake Disaster Prevention, 12(1): 210−219. (in Chinese) 肖东升,钟祖峰,2021. 基于人类动力学的地震压埋人员先验在室率研究. 灾害学,36(1):54−59.Xiao D. S., Zhong Z. F. 2021. Research on indoor rate of the earthquake buried personnel based on human dynamics. Journal of Catastrophology, 36(1): 54−59. (in Chinese) 赵大龙,白凤山,董思宇等,2015. 一种基于卡尔曼和线性插值滤波的改进三角质心定位算法. 传感技术学报,28(7):1086−1090.Zhao D. L., Bai F. S., Dong S. Y., et al. 2015. An improved triangle centroid location algorithm based on Kalman filtering and linear interpolation. Chinese Journal of Sensors and Actuators, 28(7): 1086−1090. (in Chinese) Ahn S., Lee G., Han D., 2018. A location estimating method of buried victims in collapsing area using Wi-Fi signals. In: Proceedings of the 2nd International Conference on Vision, Image and Signal Processing. Las Vegas: ACM, 49. Bouet M., Pujolle G. 2009. L-VIRT: range-free 3-D localization of RFID tags based on topological constraints. Computer Communications, 32(13−14): 1485−1494. Burian A., Kuosmanen P. 2002. Tuning the smoothness of the recursive median filter. IEEE Transactions on Signal Processing, 50(7): 1631−1639. Choi J. S., Lee H., Elmasri R., et al., 2009. Localization systems using passive UHF RFID. In: 2009 Fifth International Joint Conference on INC, IMS and IDC. Seoul: IEEE, 1727−1732. Hata M. 1980. Empirical formula for propagation loss in land mobile radio services. IEEE Transactions on Vehicular Technology, 29(3): 317−325. doi: 10.1109/T-VT.1980.23859 Joho D., Plagemann C., Burgard W., 2009. Modeling RFID signal strength and tag detection for localization and mapping. In: 2009 IEEE International Conference on Robotics and Automation. Kobe: IEEE, 3160−3165. Lima A. G. M., Menezes L. F., 2005. Motley-Keenan model adjusted to the thickness of the wall. In: SBMO/IEEE MTT-S International Conference on Microwave and Optoelectronics. Brasilia: IEEE, 180−182. Moon H., Kim C., Lee W. 2016. A UAV based 3-D positioning framework for detecting locations of buried persons in collapsed disaster area. The International Archives of the Photogrammetry Remote Sensing and Spatial Information Sciences, XLI-B8: 121−124. Pan J. J., Tang Y. Y., Pan B. C., 2007. The algorithm of fast mean filtering. In: 2007 International Conference on Wavelet Analysis and Pattern Recognition. Beijing: IEEE, 244−248. Tang S. Y., Shu X. M., Shen S. F., et al. 2014. Study of personnel positioning in large area based on pseudo base station. Procedia Engineering, 71: 481−485. -