Relationship between Shear Wave Velocity and Depth for Conventional Soils in the Bohai Sea Area
-
摘要: 土层的剪切波速是岩土地震工程中重要的物理量,本文利用多年来在渤海海域地震安全性评价中积累的资料,研究了渤海常见土类剪切波速和埋深的关系。利用非线性最小二乘法,采用指数函数、一次函数、二次函数、幂函数、"幂函数+常数函数"、"幂函数+一次函数" 6种回归模型对各类土的剪切波速和埋深的关系进行了回归分析,以拟合优度以及最小二乘拟合的误差平方和为评价指标对比了各种模型拟合效果的优劣。结果表明,"幂函数+一次函数"回归模型的拟合效果最好。此外,本文给出了该海域常见的7类土在此回归模型下的拟合公式的系数,以供工程中参考。Abstract: Shear wave velocity of soil layers is an important physical index in geotechnical earthquake engineering. Based on the experimental data accumulated over the years in the seismic safety evaluation research projects of the Bohai sea area, relationship between shear wave velocity and depth of conventional soils has been studied. In order to achieve the best fitting of experimental data, 6 regression models including exponential function, linear function, quadratic function, power function, power function combined with constant function, power function combined with linear function, have been used for regression analysis through nonlinear least square fitting method. The goodness-of-fit has been used as an index to evaluate the matching effect. The results showed that, the proposed form of function which combining power function with linear function has the best matching effect for experimental velocity data. Furthermore, the fitting coefficients of proposed function for 7 kinds of conventional soils in the Bohai Sea area have been presented for reference in engineering.
-
Key words:
- The Bohai Sea /
- Shear wave velocity /
- Depth /
- Regression analysis /
- Fitting
-
引言
台站管理是区域测震台网中心的一项基础性工作,涉及台站规划、建设和运行维护的全过程,其中,运维阶段的管理是整个管理工作的中心环节。台站运维阶段管理的目标是通过一系列管理手段的实施,提高台网的实时运行率和数据可用性,保证台网业务系统的稳定运行。台站运维工作的主要内容是台站基础设施、专业设备的维修维护和升级改造,在工作过程中会积累大量的工作日志和技术档案,这些记录是反映台站运行情况的基础资料,具有重要的应用价值。比如,通过定期对台站故障频次进行分类统计和分析,可以获悉导致台站断记的主要因素,从而为台站的升级改造提供决策依据。测震台站通常是无人值守台,台站维护需要借助计算机网络和应用软件完成。近年来,随着中国地震背景场探测项目、地震烈度速报和预警项目等一系列项目的实施,测震台网的规模日趋庞大,台站维护任务日益繁重,传统的记录和管理方式则较为简单、粗放,不利于资料的重复使用,研制专用软件,为台站运维阶段管理过程提供服务,提高管理效能,已成为当务之急。为了提高台站管理工作的效率、提高运维档案的利用率,我们研发了一套B/S结构的测震台站运维管理系统,支持用户在PC端利用通用的Web浏览器管理台站数据,在移动终端上远程检查台站运行状态。与传统C/S结构的应用程序相比,这种瘦客户端的应用架构具有更好的可维护性、可扩展性、可移植性和较低的使用门槛。
WebGIS是GIS技术与Web技术的有机结合,是在网络环境下传输、存储、处理、分析、显示与应用地理空间信息的分布式系统(孟令奎等,2005)。开源WebGIS是源代码开放的WebGIS平台(魏波等,2009),与商业平台相比,具有低成本、可定制、跨平台的特点。WebGIS被广泛应用于社会各个领域(刘吉夫等,2003)。近年来,在地震行业,也涌现出一系列基于WebGIS的应用,涉及监测预报、应急救援、震害防御3大体系(赵曦等,2014;胡斌等,2014;屈佳等,2014;班亚等,2014)。在系统中,我们利用GeoServer搭建了开源的台站电子地图服务器,利用OpenLayers实现了台站电子地图的远程接入。
1. 需求分析
1.1 测震台站的维护工作
在区域测震台网中心,台站运维的日常工作通常包括检查状态、发现故障、远程维修、监控服务器/工作站、通信链路等。台站发生故障时,首先判断故障类型(表 1),然后登记台站维护日志。台站故障总体可分为断记和病态运行2种情况。断记是指实时数据流中断,导致断记的主要因素有外部供电停止、通信链路中断、数据采集器无输出等。病态运行是指数据流未中断,但仪器的健康状态指标超限,如地震计零点漂移、数据采集器的钟差过大等。仪器的病态运行会降低观测数据的质量,影响台网业务系统的资料产出,因此对仪器的健康状态进行检查是十分必要的。通过检查,不仅可以纠正病态运行,还可以预防断记。例如,可以对供电电压进行监视,当电压低于11V时,通常意味着外部供电停止,这时可以通过及时更换蓄电池来保证运行。有些故障可以通过远程维修解决,例如地震计零点漂移、数据采集器无输出等。无法远程修复的故障必须进行现场维修,涉及到设备更换的,应当记录仪器管理日志。远程或现场维修以后,应当及时更新台站维护日志。在维护过程中,要记录4类信息:① 台站基本参数,主要内容包括台站代码、地理坐标、场地类型、仪器类型、启用时间等;② 台站维护日志,主要内容包括台站的故障类型、处置类型、故障持续时间等。③ 仪器管理日志,主要内容包括仪器的标识、序列号、台站、启用时间等;④ 台站仪器运行状态,包括电压、温度、钟差、存储空间、机械零位等11个参数(董一兵等,2013)。
表 1 台站故障类型编码Table 1. List of codes of fault type序号 代码 描述 1 1000 正常 2 1101 外部供电停止 3 1102 供电设备故障 4 1199 其他供电影响 5 1201 地震计无输出 6 1202 地震计零点漂移 7 1209 地震计不明原因故障 8 1211 数据采集器无输出 9 1212 数据采集器时钟错误 10 1213 数据采集器参数设置错误 11 1219 数据采集器不明原因故障 12 1231 脉冲标定无输出 13 1232 脉冲标定异常 14 1299 未知原因的专业设备故障 15 1301 通信链路中断 16 1302 通信设备故障 17 1701 基础设施改造 18 1702 盗窃 19 1703 观测环境变化 20 1799 其它人为干扰 21 1801 雷电 22 1899 其它天灾影响 23 1998 原因不明的数据中断 24 1999 其他因素 在以往的台站维护工作中,由于缺少专用的软件平台,各环节通常需要借助不同的工具完成,环节之间缺乏有效衔接和集成,导致工作程序繁琐、工作效率不高。例如,检查仪器状态要使用厂家提供的客户端软件,登记故障要使用JOPENS系统的台站日志,现场维修、仪器更新要借助纸介质或电子表格进行记录,这些过程会逐渐积累大量独立存储的档案资料,这些资料具有重要的参考价值,但难于保存、利用困难,尤其是纸介质记录,更容易散失。采用集成化的台站运维管理系统之后,可以对各类台站、仪器以及运维档案进行集中存储和管理,提高运维工作效率和技术档案的利用率。
1.2 功能和性能要求
系统应能够对台站运维过程中形成的各类信息进行管理,并提供人机交互式的图形用户界面,主要功能包括:
(1)测震台站运维信息管理:支持对测震台站基本参数、台站维护日志和仪器管理日志的增、删、改、查操作。
(2)测震台站仪器状态监视:对测震台站健康状态的信息进行定时采集和发布,对异常及时报警。
(3)测震台站故障统计:对测震台站的故障类型、处置类型按频次进行统计,输出统计图表。
(4)测震台站WebGIS:支持缩放操作,提供地图配置选项。
系统的性能应满足下列要求:
(1)健壮性:系统运行稳定;具有可靠的容错机制。
(2)安全性:系统应与测震台网目前的技术系统良好兼容;提供认证机制,限制对重要数据的访问。
(3)并发性:支持多用户并发操作。
(4)易用性:提供友好的图形用户界面。
2. 系统设计
2.1 网络结构
系统的网络拓扑结构如图 1所示,服务器端采用分布式部署方式。采用分布式部署时,需要分别搭建几种服务器环境,包括Web、WebGIS、消息和数据库服务器。其中,Web服务器使用Tomcat 6.0,它是Web应用程序的容器,打包后的Web源程序要部署到它的webapps目录下才能执行(Chetty,2009);WebGIS服务器使用GeoServer 2.5,用于提供电子地图;消息服务器使用ActiveMQ 5.10,用于管理台站仪器状态消息队列——服务器端程序从数据采集器接收状态数据,然后传递给ActiveMQ,客户端程序从ActiveMQ订阅状态消息;数据库服务器使用PostgreSQL 8.4+PostGIS,PostgreSQL存储属性数据,PostGIS是它的空间扩展,用于存储空间数据。
2.2 系统架构
系统采用分布式的多层架构(图 2)。在服务器端,仪器状态数据采集服务从数据采集器查询得到仪器的状态数据,并传递给ActiveMQ;ActiveMQ管理名为SOH.TOPIC的主题队列,缓存状态数据并持久化写入PostgreSQL数据库;PostgreSQL存储属性数据,包括台站基本参数、台站维护日志、仪器管理日志和台站仪器运行状态;PostGIS存储空间数据,例如行政区边界和台站分布等;测震台站管理业务逻辑,包括台站运维信息管理、仪器状态监视和故障统计等功能模块,都被部署到Tomcat容器中,并发布为Web服务,在前端使用JSP调用;GeoServer将PostGIS数据源发布为WebGIS服务(Iacovella等,2013),在前端使用Openlayers调用(Gratier等,2015)。在客户端,授权用户使用通用的Web浏览器即可访问系统。
2.3 仪器状态监视模块的设计
目前,在测震台网业务系统中,已有对台站实时地震波形数据进行监视、报警,并对台站实时运行率和数据完整率进行统计的应用软件,但对台站健康状态的数据仍缺乏高效的利用。近年来,各区域测震台网对通信链路、避雷系统和地震仪器相继进行了升级改造和更新换代,台网的实时运行率和数据完整率得到了可靠保证。现阶段,着力提高数据质量应成为台网运维工作的重点之一,数据质量直接影响台网的监测效能。仪器的健康状态是影响数据质量的重要因素,例如,数据采集器的供电电压、机箱温度直接影响其正常运行,钟差影响授时精度,地震计的零位飘移降低系统的动态范围等等,因此,监视仪器的健康状态十分必要。JOPENS系统在数据库中定义了健康状态数据表(Soh.table)的结构,其中,能够直接从数据采集器获取的数据及说明如表 2所示。表 2中定义的正常值区间是软件的缺省配置,表示监测值落入区间内为正常,否则为异常。
表 2 台站仪器状态参数表Table 2. SOH parameters序号 名称 含义 正常值区间 单位 1 Pwr_input 供电电压 [-9,18] V 2 Pwr_output 输出电压 [-9,18] V 3 Pwr_battery 电池电压 [0,9] V 4 Dsp_temp 主板温度 [-20,60] ℃ 5 Clk_diff 钟差 [-100,100] μs 6 Frq_diff 频差 [-10,10] PPB 7 Drive_free_space 存储器可用空间 [20,100] % 8 Z_mp_volts Z通道零位 [-1,1] V 9 N_mp_volts N通道零位 [-1,1] V 10 E_mp_volts E通道零位 [-1,1] V 本模块由3个服务构成:仪器状态采集服务、消息服务和Web服务。其中,仪器状态采集服务与台站的数据采集器进行通信,收集状态数据,并交给消息服务(董一兵等,2015);消息服务管理状态消息的队列,支持数据持久化选项;Web服务提供用户访问状态数据的界面。
2.4 用户操作流程
根据测震台站日常维护工作的一般程序设计了系统用户的操作流程(图 3)。用户登录系统后,可以并行地执行台站信息管理、仪器状态检查和故障统计等3类操作,当发现台站故障时,首先需要记录台站维护日志;然后,根据故障类型进行维修,若更换了仪器,需要记录仪器管理日志;最后,当台站恢复正常后,更新台站日志,将处置类型、故障结束时间、日志内容和备注等信息补充完整。
3. 系统实现与应用
3.1 系统主页
在客户端,用户在浏览器中键入服务器地址即可进入系统主页(图 4)。主页顶部是系统LOGO和快捷按钮区域,中部是测震台站专题电子地图。地图左侧提供了“鱼骨”控件,可以对地图进行缩放;地图右侧是WebGIS配置区域,用户可选择关心的图层予以显示。无论是局域网用户还是移动用户,只要与服务器连通,即可访问系统,便于现场工作人员远程检查台站运行状态。
3.2 测震台站运维信息管理
单击主页顶部的“测震台站管理”快捷按钮,可从地图左侧弹出功能区域(图 5),自上而下依次是台站参数管理、仪器管理和台站维护日志管理。在“台站参数管理”区域,缺省显示测震台站列表,用户可以通过单击字段名称对列表记录进行排序,也可以修改单页显示的记录条数。工具栏上提供了基于“台站代码”关键字的模糊检索功能。单击某1条台站记录可以使WebGIS定位到该台站。缺省状态下,用户以访客身份进入系统,仅能浏览部分数据,没有修改数据的权限。要使用系统的全部功能必须以管理员的身份登录后才可通过单击“新增”、“编辑”和“删除”按钮执行相应操作。仪器管理和台站维护日志管理的界面与图 4类似,不再赘述。利用这几项功能,当台站参数发生变更、仪器进行了置换或者更新、台站基础设施更新时,台站维护人员可以将这些变更记录到数据库中。当需要了解某台站的资料时,可以利用模糊检索查询该台站的历史记录。
3.3 对测震台站仪器状态的监视
在台站列表双击某台站即可打开该台站的仪器状态监视页面(图 6)。页面上提供了2种视图,即上部的数据列表和下部的仪表盘,2种方式同步刷新。根据系统配置的阈值范围(表 2),仪表盘上的绿色区域表示正常,红色区域表示异常,指针随区域改变颜色,用户可以通过观察仪表的读数和指针的偏移了解台站仪器的工作状态,也可以自定义正常值的区间。比如,在实际工作中,测震仪器供电电压的正常值在12V左右,当遇到外部供电停止时,供电电压会逐渐下降,降至10V以下时,数采将不能正常工作,因此可以将供电电压的阈值区间设置为[10, 13],这样一来,当某台站的供电电压低于10V时,系统将发出报警提示。系统钟差和地震计机械零位也是对数据质量影响很大的因素,用户应当重点关注。
3.4 测震台站故障统计
从“测震台站管理”功能区可打开台站故障统计页面。用户设置“台网”、“台站”、“故障类型/处置类型”、“开始时间”和“结束时间”之后单击“统计”按钮,系统将按照用户要求对台站维护记录进行统计并输出统计图表,用户可以获得各类故障的分布情况,并从中找出主要故障。图 7是河北测震台网2014年故障统计饼图,从图中可以看出,外部供电停止在所有故障类型中占比重最高,达到48%,这提示用户下一步应重视台站的市电供应情况,及时更新UPS和蓄电池。通信链路中断的情况占比达到13%,也属于比较多发的故障类型,这提示用户应当加强与网络运营商的沟通和协调,督促运营商提高通信链路质量,加快故障处理速率。
3.5 系统应用
河北是一个震害严重的区域(王想,2016),提高河北测震台网运维工作的效率有利于保证台网的监测能力,使之更好服务于防震减灾事业。2015年以来,系统在河北测震台网进行了应用。我们将2014年的台站故障记录导入系统,并进行了统计,结果表明,外部供电停止和通信链路中断是导致台站断记的最主要因素,占到故障总数的60%。2015年,台网有针对性地对台站进行了升级改造,更新了台站的供电设备并改进了通信设备的供电方式,使得因这2种因素导致的断记得到有效控制,台网的实时运行率和数据完整率均较上年有所提高。系统的应用,切实提高了台站维护工作的效率,为保障台网运行率和数据质量发挥了积极作用。
4. 结论
本文介绍了综合利用Java Web和GeoServer开发技术,设计、开发的一套测震台站运维管理系统。系统采用B/S结构,免去了客户端程序的安装和维护,增强了系统的灵活性和可维护性,降低了使用门槛和维护成本。系统完全采用成熟的开源技术开发,在保证性能的同时控制了研发成本;通过引进开源的WebGIS开发技术GeoServer+PostGIS+Openlayers,实现了功能丰富、性能稳定、界面友好的台站专题电子地图;测震台站运维信息管理实现了对台站参数、仪器参数和台站维护日志的管理;测震台站仪器状态监视实现了对台站状态参数的采集、存储、发布与报警,为台站数据质量的实时监控提供了工具;测震台站故障统计实现了对台站故障的分类统计,可为台站建设、更新提供参考依据。系统在河北测震台网的应用切实提高了台网的实时运行率和数据完整率,具有一定的推广应用价值。
-
表 1 各类土在各回归模型下的拟合优度
Table 1. Goodness-of-fits of 7 kinds of soils under different regression models
拟合函数 土类 粉质粘土 砂质粉砂 粉砂 粉质细砂 细砂 砂质粉土 粉土 R2 RSS /105 R2 RSS /104 R2 RSS /104 R2 RSS /105 R2 RSS /104 R2 RSS /103 R2 RSS /104 (1)v=aebh 0.861 3.46 0.865 4.81 0.883 3.20 0.873 9.83 0.900 3.26 0.925 9.36 0.865 2.45 (2)v=ah+b 0.915 2.10 0.901 3.51 0.913 2.37 0.909 7.04 0.916 2.73 0.962 4.79 0.928 1.31 (3)v=ah2+bh+c 0.915 2.10 0.901 3.51 0.913 2.37 0.909 7.04 0.916 2.73 0.962 4.79 0.928 1.31 (4)v=ahb 0.935 1.61 0.926 2.64 0.901 2.71 0.946 4.19 0.922 2.55 0.992 1.03 0.979 0.38 (5)v=ahb+c 0.943 1.42 0.932 2.42 0.919 2.22 0.949 3.91 0.928 2.33 0.992 0.94 0.981 0.34 (6)v=ahb+ch+d 0.944 1.40 0.934 2.33 0.920 2.20 0.950 3.82 0.931 2.25 0.993 0.91 0.981 0.34 注:R2表示拟合优度,RSS表示误差平方和,其单位为m2/s2。 表 2 各土类采用v=ahb+ch+d形式拟合的系数
Table 2. Fitting coefficients for 7 kinds of soils while using the formula v=ahb+ch+d
v=ahb+ch+d拟合系数 土类 粉质粘土 砂质粉砂 粉砂 粉质细砂 细砂 砂质粉土 粉土 a 76.4017 123.9201 141.3937 120.0582 145.2684 101.4440 49.1360 b 0.2760 0.1907 0.1190 0.2241 0.1630 0.2769 0.4209 c 1.4209 1.1374 1.7891 0.7949 1.1223 0.4958 1.6224e-6 d 20.2951 0.0013 0.0049 3.0568e-5 5.2199e-5 1.0829e-5 62.3314 -
陈义兰, 吴永亭, 刘晓瑜等, 2013.渤海海底地形特征.海洋科学进展, 31(1):75-82. http://www.cnki.com.cn/Article/CJFDTOTAL-HBHH201301010.htm 程祖锋, 李萍, 李燕等, 1997.深圳地区部分岩土类型剪切波速与深度的关系分析.工程地质学报, 5(2):163-168. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ702.010.htm 狄圣杰, 汪明元, 张昆等, 2014.江苏近海地层原位剪切波速相关特性及预测方法研究.海洋学报, 36(3):127-133. http://www.cnki.com.cn/Article/CJFDTOTAL-SEAC201403014.htm 高玉峰, 刘汉龙, 2003.合肥膨胀土剪切波速的特征分析.岩土工程学报, 25(3):371-373. http://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200303026.htm 胡钧, 杜坚, 1997.上海场地土的脉动特性及剪切波速特征.物探与化探, 21(4):305-307, 299. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH199704010.htm 刘红帅, 郑桐, 齐文浩等, 2010.常规土类剪切波速与埋深的关系分析.岩土工程学报, 32(7):1142-1149. http://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201007029.htm 齐鑫, 丁浩, 2012.下辽河平原区剪切波速与土层埋深关系分析.世界地震工程, 28(3):151-156. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201203028.htm 邱志刚, 薄景山, 罗奇峰, 2011.土壤剪切波速与埋深关系的统计分析.世界地震工程, 27(3):81-88. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC201103013.htm 苏旭, 杨丽萍, 邹海宁, 2009.西宁市区土层剪切波速随深度变化特征.高原地震, 21(1):12-16. http://www.cnki.com.cn/Article/CJFDTOTAL-GYDZ200901001.htm 徐秦, 薛茜, 徐睿, 1992.浅论曲线拟合中的相关指数R2.中国卫生统计, 9(6):44-45. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=zgwt199206022&dbname=CJFD&dbcode=CJFQ 战吉艳, 陈国兴, 刘建达, 2009.苏州城区深软场地土剪切波速与土层深度的经验关系.世界地震工程, 25(2):11-17. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC200902002.htm 张世强, 2002.曲线回归的拟合优度指标的探讨.中国卫生统计, 19(1):9-11. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGWT200201003.htm 赵淑芳, 杨宏亮, 2012.场地剪切波速的特征分析.云南大学学报(自然科学版), 34(S2):267-271. http://www.cnki.com.cn/Article/CJFDTOTAL-YNDZ2012S2045.htm 中国科学院海洋研究所海洋地质研究室, 1985.渤海地质.北京:科学出版社, 1-11. 周锡元, 王广军, 苏经宇, 1990.场地·地基·设计地震.北京:地震出版社, 28-172. Hardin B. O., Drnevich V. P., 1972. Shear modulus and damping in soils:design equations and curves. Journal of Soil Mechanics & Foundations Division, 98(7):667-692. Lew M., Chieruzzi R., Campbell K. W., et al, 1981. Correlations of seismic velocity with depth. In:Proceedings of the 1st International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics(ICRA GEESD) U.S.A. St. Louis, Missouri, 1981. University of Missouri-Rolla. -