我国平原地区幅员辽阔、城市林立,人口密集、经济发达,防震减灾工作尤为重要。高邮市位于长江下游的苏北平原,以高邮市城区为例,研究这类场地的地震动效应特征,对影响地震动参数小区划成果的几个主要因素进行总结与探讨。考虑多方因素进行钻孔布设;增加地震动概率水平;利用统计分析理论,补充数据,建立土层力学模型;利用地表峰值加速度放大倍数等,研究如何减轻等效线性化方法低估高频地震动的缺陷。 China’s flat areas are vast, and as a largely urbanized, densely populated, economically developed country, the job of earthquake prevention and disaster reduction is especially important. Gaoyou lies in north Jiangsu plain which located in the lower reaches of Yangtze river. Taking Gaoyou city as an example, this paper studies the characteristics of ground motion, and summarizes and discusses several main factors influencing ground motion parameter microzonation. Borehole layout after considering many factors, increasing the level of ground motion probability, using statistical analysis theory, supplementing data, building soil mechanics model, using surface peak acceleration to magnify times, to study on how to alleviate the defect of equivalence method on underestimating high frequency ground motion.
—以高邮市城区为例
陶小三,杨伟林,李细兵,邵斌
江苏省地震工程研究院,江苏 南京
收稿日期:2017年8月21日;录用日期:2017年9月3日;发布日期:2017年9月11日
我国平原地区幅员辽阔、城市林立,人口密集、经济发达,防震减灾工作尤为重要。高邮市位于长江下游的苏北平原,以高邮市城区为例,研究这类场地的地震动效应特征,对影响地震动参数小区划成果的几个主要因素进行总结与探讨。考虑多方因素进行钻孔布设;增加地震动概率水平;利用统计分析理论,补充数据,建立土层力学模型;利用地表峰值加速度放大倍数等,研究如何减轻等效线性化方法低估高频地震动的缺陷。
关键词 :平原地区,城市,高邮,地震动参数小区划
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
我国平原面积约为112万平方公里,占全国面积的12%。这些平原主要是由江、河、湖、海的泥沙冲积而成,地势坦荡、水网稠密、土壤肥沃。平原地区城市林立、人口密集、经济发达,开展相关城市的地震小区划工作,有利于政府实施土地利用规划,有利于提高城市防震减灾能力,减轻地震灾害、经济损失、人员伤亡。
高邮市位于长江下游的苏北平原,地势低洼,场地上覆土层较厚。近年来,尽管在高邮市对部分重大工程开展过场地地震安全性评价工作,但这些工程数量上相对较少,仅限局部地区,对整个高邮市而言控制作用有限。高邮市临高邮湖,河网密布,场地不均,地质条件复杂,研究程度不足等原因,土地开发利用和工程设施建设面临一些风险和不确定性,需开展地震小区划等基础性工作,为经济建设和区域开发提供技术依据。本文以高邮市城区为例,研究这类场地的地震动效应特征,对影响地震动参数小区划成果的几个主要因素(钻孔布设、土层力学模型的建立、土层地震反应分析)进行总结与探讨。
高邮市地处长江下游——南黄海地震带的苏北地震构造区内,高邮市及其邻近地区历史上也曾发生过多次破坏性地震,如1624年2月10日江苏扬州6级地震,中、远场强震对高邮市城区最大影响烈度为VII度。2012年7月20日在江苏高邮、宝应交界还发生了M4.9级地震,造成了人员伤亡和财产损失。近场区现代小震频发,1990年在高邮南有小震群活动。高邮市区及其邻近地区地质构造比较复杂,分布有多条第四纪活动的断裂,具备发生中强破坏性地震的构造条件。
高邮市位于苏北平原西南缘,为第四系覆盖区,地形平坦,第四系沉积厚度约180 m,上新世以来形成的沉积物(N + Q)厚度约900 m,基岩为古近系粉砂层(E)。
岩土工程勘察主要是钻探勘察,用于探测场地工程地质条件,划分工程地质单元,建立土层地震反应分析模型,是地震小区划工作中的一项重要内容,直接影响到工作量、经费投入、工作质量。平原城市经济发达,基础设施较为完善,重点或重大工程多,故岩土勘察资料非常丰富,应充分收集勘察资料。《工程场地地震安全性评价宣贯教材》 [
在收集、整理、分析高邮市城区的工程地质、水文地质、地形地貌、抗震防灾规划、城市总体规划等资料的基础上,初步划分高邮市城区为2个工程地质单元分区。《工程场地地震安全性评价宣贯教材》 [
《中国地震动参数区划图》 [
地震动参数小区划,是由地震危险性分析确定各计算点的基岩地震动后,开展场地地震动效应分析,得到地表地震动参数的分布情况,再结合地震小区划范围内的地形、地貌和岩土性质的特点,综合分析后,勾划出具有一定概率水平的地震动参数小区划图并确定各分区的设计地震动参数。
城市是一个系统工程,工程种类很多,抗震设计时有需要考虑50年概率水平的,如建构物等;有需要考虑100年概率水平的,如大型桥梁、地铁工程等。《中国地震动参数区划图》 [
地震小区划项目投资较大(几百万元),工期较长(几个月,甚至1年多),如果多提供几种实用概率水平地震动参数,其野外工作的投入是一样的,仅需增加室内工作,从而提高资金使用率、成果使用范围。高邮市城区地震小区划工作中,项目组提供了7种概率水平的地震动参数(表1),具体工程可根据其重要性及工程所允许承担的风险,选用相应的地震动参数进行抗震设计。
地质体的沉积环境复杂多变,故土层性质在纵横向均有变化。横向均匀场地一般简化为成层土层场地,采用一维波动模型进行场地地震反应分析;横向非均匀场地,应简化为二维或三维模型,采用有限单元法(FEM )、有限差分法(FDM)及边界单元法(BEM)等进行场地地震反应分析。
平原地区城市范围各工程地质单元内地层分布基本稳定,土层水平成层较好,一般可简化为一维波动模型,常采用一维等效线性波动法进行土层地震反应分析。如局部范围内岩土层水平界面起伏较大,宜增加二维或三维场地模型。一维等效线性波动法中,建立土层力学模型所需场地土的静、动力性能参数有:土层剪切波速、土的密度、动剪切模量及阻尼比与剪应变关系曲线等。
概率水平 | 中国地震动参数区划图 | 地震小区划(规范) | 高邮市城区地震小区划 |
---|---|---|---|
50年 | 63%、10%、2% | 63%、10%、2% | 63%、10%、2% |
100年 | 1% | 63%、10%、3%、1% |
表1. 地震区划图的概率水平表
土层剪切波速的确定,一般采用原位测试结果。但有时由于受时间、费用、设备和技术条件等方面因素的制约,可能会没有条件测试或测试数据不全。部分平原地区沉积时间较短,控制性钻孔深度102 m终孔时,其土层剪切波速依然小于500 m/s;有时为提高工作质量、工作精度,需收集、补充部分非控制性钻孔进行土层地震反应分析,非控制性钻孔的波速数据往往不全。《工程场地地震安全性评价》 [
模型1:
模型2:
模型3:
式中Vs以m/s计,H以m计,a、b、c等为待定系数。
研究表明,Vs~H的经验关系有很强的地区性,同一地区同一土类剪切波速与土层深度的相关性较好,有足够的工程精度。在土层剪切波速资料不足时,应收集相关资料或钻探深孔进行原位测试,再采用最小二乘法进行某土层剪切波速Vs与土层深度H的回归分析,用相关系数R进行显著性检验,并根据剩余标准差SD值选择较优的Vs~H经验关系,进而推出各土层的剪切波速。
动剪切模量比、阻尼比与动剪应变关系曲线由现场取得典型原状土样,再室内试验(共振柱或动三轴试验)确定。很多学者针对某一地区土的动剪切模量比和阻尼比进行了大量试验研究,并取得一些较为有价值的研究成果,认为同一地区同类土的动剪切模量比、阻尼比与动剪应变关系曲线变化不大 [
其中:A、B、γ。、n、λ。是与土性有关的拟合参数;λmin是与初始动剪切模量Gmax相对应的土的最小阻尼比,与土的性质、固结状态等因素有关。
随着防震减灾工作的推进,高邮市城区积累了大量各类土动剪切模量比和阻尼比与剪应变关系的数据。为减小数据离散性,提高土层地震反应计算效率、计算精度,项目组采用式(4)和式(5)分别对G/Gmax~γ和λ~γ关系试验数据进行回归拟合。统计样本位于同一城区、样本数量较多、全是共振柱试验数据,故统计分析结果(表2)离散性很小,其R-Square均达到0.98以上,因而给出的统计结果具有较高的可靠性、适用性,可用于相关工作的借鉴参考。
场地土层是一个物理滤波器,对不同方向传来的入射波有滤波性能。场地软弱土层(剪切波速较低)对高频信号有滤波性能,对低频信号则有放大作用;与此相反,场地坚硬土层(剪切波速较高)对低频信号有滤波性能,但对高频信号则有放大作用。场地土层愈厚,滤波作用愈明显。厚层软弱场地土的地震动以长周期为主,柔性结构物(如高耸结构物)的震害比较严重;薄层坚硬场地土的地震动以短周期为主,刚性结构物的震害比较严重。
平原地区主要分布在大河两岸和濒临海洋的地区,水网密布,沉积土层很深,有的甚至达到数千米。浅部土层基本为软土,建筑场地类别一般为III类和IV类,地震动以长周期为主,IV类场地的地表Amax(加速度峰值)或反应谱平台较III类场地低,但特征周期更长(图1),显示出不同的场地地震动效应。
研究土层地震反应的理论计算方法中,影响较大的有等效线性化方法和非线性方法两类。前者如中国地震局系统的ESE软件包、XQH软件包等;后者包括有效应力和总应力方法。许多研究者注意到,当基岩输入地震动较强时,等效线性法高频段算得的频响放大倍率比实测结果偏低。在传统的等效线性化波动法中,等效剪应变通常固定取最大剪应变的0.65倍,对于等效线性化方法的改进,集中于对等效剪
土类名称 | 参数 | γ | |||||||
---|---|---|---|---|---|---|---|---|---|
5e-6 | 1e-5 | 5e-5 | 1e-4 | 5e-4 | 1e-3 | 5e-3 | 1e-2 | ||
1.2-2粉质粘土:软塑 | G/Gmax | 0.9906 | 0.9809 | 0.9065 | 0.8271 | 0.4891 | 0.3267 | 0.0926 | 0.0497 |
λ | 0.0195 | 0.0230 | 0.0428 | 0.0594 | 0.1155 | 0.1387 | 0.1698 | 0.1753 | |
2.2-3淤泥质粉质粘土:流塑 | G/Gmax | 0.9921 | 0.9842 | 0.9246 | 0.8591 | 0.5473 | 0.3762 | 0.1073 | 0.0566 |
λ | 0.0187 | 0.0246 | 0.0497 | 0.0675 | 0.1226 | 0.1451 | 0.1753 | 0.1805 | |
3.2-4粉土:中密 | G/Gmax | 0.9892 | 0.9788 | 0.9041 | 0.8267 | 0.4992 | 0.3381 | 0.0984 | 0.0532 |
λ | 0.0196 | 0.0248 | 0.0483 | 0.0654 | 0.1170 | 0.1372 | 0.1640 | 0.1687 | |
4.3-1粉质粘土:可塑 | G/Gmax | 0.9910 | 0.9830 | 0.9230 | 0.8600 | 0.5600 | 0.3900 | 0.1130 | 0.0590 |
λ | 0.0184 | 0.0227 | 0.0422 | 0.0570 | 0.1069 | 0.1291 | 0.1609 | 0.1666 | |
5.4-1粉质粘土:可塑 | G/Gmax | 0.9899 | 0.9806 | 0.9152 | 0.8471 | 0.5378 | 0.3715 | 0.1085 | 0.0579 |
λ | 0.0189 | 0.0237 | 0.0451 | 0.0608 | 0.1109 | 0.1321 | 0.1615 | 0.1668 | |
6.4-1b粉土:中密 | G/Gmax | 0.9910 | 0.9830 | 0.9273 | 0.8687 | 0.5846 | 0.4167 | 0.1268 | 0.0678 |
λ | 0.0176 | 0.0220 | 0.0412 | 0.0554 | 0.1030 | 0.1246 | 0.1568 | 0.1628 | |
7.5-1粉质粘土:可塑 | G/Gmax | 0.9917 | 0.9839 | 0.9281 | 0.8685 | 0.5823 | 0.4165 | 0.1311 | 0.0717 |
λ | 0.0171 | 0.0214 | 0.0405 | 0.0548 | 0.1022 | 0.1235 | 0.1552 | 0.1612 | |
8.6-1粉质粘土:硬塑 | G/Gmax | 0.9934 | 0.9873 | 0.9439 | 0.8966 | 0.6473 | 0.4826 | 0.1603 | 0.0876 |
λ | 0.0151 | 0.0183 | 0.0331 | 0.0448 | 0.0880 | 0.1102 | 0.1473 | 0.1550 |
表2. 高邮市城区场地土的G/Gmax~γ和λ~γ典型值
图1. 不同建筑场地地表加速度反应谱曲线图
应变的确定。已有的研究工作大体有两种途径 [
式中,X为折减系数,M为震级,R为震中距。
本文以高邮市城区地震小区划中某计算控制点为例,来研究如何减轻等效线性化方法低估高频地震动的缺陷。进行土层地震反应分析时,为减小不同输入地震波的离散性,规范要求每一概率水平至少输入3条地震波,现对此计算控制点力学模型输入6条超越概率100年1%的基岩地震波。6条地震波的地表加速度峰值分别为0.287 g、0.292 g、0.295 g、0.302 g、0.309 g、0.316 g,最大值与最小值相差约10%,平均值为0.300 g。式(6)和(7)两类改善等效线性化方法存在复杂性和不确定性,《工程场地地震安全性评价宣贯教材》 [
图2. 地表加速度反应谱曲线图
根据地震危险性分析结果,该计算控制点超越概率100年1%等效震级M为6级,R为18 Km,基岩加速度峰值为0.217 g。由沈建文等建议的EQLM(M,R)方案 [
不同超越概率水平 | |||
---|---|---|---|
50年2% | 100年3% | 100年1% | |
基岩地震动Amax(g) | 0.148 | 0.161 | 0.217 |
地表地震动Amax(g) | 0.218 | 0.232 | 0.310 |
放大倍数 | 1.47 | 1.44 | 1.42 |
表3. 地表地震动峰值加速度的放大倍数
采用上述方法,从表3可知,100年1%地表Amax取0.310 g比较合适,取0.314 g就不对了。最终确定此计算控制点的地震动参数为(图2(b)):Amax = 0.310 g,βmax = 2.5,αmax = 0.775,Tg = 0.90 s。相对于常规方法,Amax或αmax提高了3%,有效地减缓了加速度反应谱规准的带帽现象,减轻了传统等效线性化方法的低估缺陷。
地震动参数小区划的结果可以是分区图或者等值线图,二者各有优缺点。分区图法编图简单,与工程地质单元结合紧密,便于工程应用,此种方法用得比较多。鉴于《中国地震动参数区划图》 [
本文对影响高邮市城区地震动参数小区划结果的几个主要因素,进行了总结和探讨,成果可供同类平原地区城市借鉴。
1) 本文分析平原地区工程地质条件和地震动特性的关系,不同场地条件体现出不同的地震动效应,二者有密切的依存关系。为了解和探测场地工程地质条件,必须开展钻探勘察。钻孔布设的合理与否,直接影响到工作量、经费投入、工作质量。钻孔布设应结合工程地质单元划分、城市规划重点工程分布(如商务区、医院、大型基础设施等),并注意适当均匀分布(钻孔位置间隔1 Km内),这样最终成果既符合规范要求,又更有针对性、实用性。
2) 近年来,我国城市发展极为迅速,各类新结构、新类型的建设工程层出不穷,对地震小区划工作提出更高的要求。如多概率水平的地震动参数、超长周期反应谱、不同阻尼比的反应谱等等,需要认真加以研究。为提高资金利用率、成果使用范围,本文建议地震动参数小区划提供7个概率水平地震动参数。
3) 土层剪切波速、土层动力非线性关系是土层力学模型的重要基础资料。土层剪切波速的确定,一般采用原位测试结果。当数据不全时,可根据土的常规物理力学性能等,采用经验关系确定相应土层的剪切波速。为减小数据离散性,提高计算效率、计算精度,本文分别对高邮市城区各层土动剪切模量比和阻尼比与剪应变关系的试验数据进行统计回归分析,得到高邮市城区典型场地土的动力非线性参数,统计结果离散性很小,具有较高的可靠性、适用性。
4) 平原地区土层地震反应分析中,等效线性化方法运用较多,但其存在当基岩输入地震动较强时,高频段算得的频响放大倍率比实际场地的实测结果偏低。研究人员提出EQLM(M)和EQLM(M,R)两种理论改善方案,有其合理性,但实际应用时由于参数的复杂性可能产生较大的不确定性。本文提出利用地表峰值加速度放大倍数来综合控制,方法简单,有其实用性。
陶小三,杨伟林,李细兵,邵 斌. 平原地区城市地震动参数小区划研究—以高邮市城区为例Study on the Ground Motion Parameter Microzonation for Cities in the Flat—Taking Gaoyou City as an Example[J]. 应用物理, 2017, 07(09): 261-269. http://dx.doi.org/10.12677/APP.2017.79032