Magnitude Simulation of Aftershocks Based on G-R Relationfor MAT Earthquake in Earthquake Emergency Response Exercise
-
摘要: 在现有的日常地震演练过程中，与应急救援密切相关的地震余震信息产品较为缺乏，直接影响发震构造的判断以及影响场修正等关键环节。本文从计算机系统提供的均匀分布随机数出发，运用反函数法模拟生成余震序列，并进行系统检验，证实该方法产生的余震序列满足G-R频次关系。模拟生成的余震震级数据既有助于增强地震应急救援演练的现实性，也有助于丰富地震应急宣传产品，提升地震部门的履职能力。Abstract: Insufficient aftershock information, which is one of the key factors to emergency rescue, directly affects the essential steps such as seismogenic structure determination and influence field correction in the daily exercise against earthquake. In this paper we simulate and generate aftershock sequence by using inverse function method on uniformly distributed random numbers generated by computer system, Then we carry out system examination with the simulation results, and finally prove that the aftershock sequence generated by this method satisfies the magnitude distribution from G-R relation. The simulated aftershock sequence not only helps to enhance the reality of emergency rescue exercise against, but also contributes to enrich promotional products of earthquake emergency response and improves the executing capacity of earthquake departments.
-
引言
随着城市地下空间的不断开发和利用，基坑工程数量日益增多，为分析基坑引起的变形(李平等，2010)，数值模拟方法得到广泛应用。目前，针对基坑开挖分析的岩土工程软件较多，如ABAQUS、ANSYS、FELAC、ADINA等，如果对软件适用性缺乏了解、应用不当，将造成分析结果不可靠，甚至产生错误结论。Hashash等(2010)、Grande(1998)、Potts等(2001)及徐中华等(2010)根据不同本构模型下基坑开挖变形结果与实测资料的对比分析，得出具有实用价值的研究成果，但这些成果多适用于特定土体条件下的工程，普遍使用性较差。Potts等(2002)指出，用于基坑工程数值模拟的土体本构模型既应能刻画分析问题特点，又应简单易实现。蒋明镜等(2012)基于PFC2D和FLAC3D软件对Mohr-Coulomb模型、Drucker-Prager模型在基坑分析中的适用性进行探讨。秦会来等(2012)对比分析基于ABAQUS软件修正剑桥模型和Mohr-Coulomb模型的基坑开挖二维数值模拟结果，认为Mohr-Coulomb模型不能反映加卸载模量差异和土体压硬性，故不适用于基坑开挖引起的变形模拟。但目前鲜有开展ABAQUS软件内置本构模型(Mohr-Coulomb模型、Drucker-Prager模型、修正剑桥模型)的基坑开挖变形三维数值模拟可靠性分析。笔者以ABAQUS软件为平台，分析内置本构模型对基坑开挖变形数值分析结果的影响，讨论内置本构模型数值模拟结果的合理性，在此基础上分析基于修正剑桥模型的开挖方式和二维简化分析模型的影响，并对比三维模型和二维模型模拟结果，论证基坑稳定性分析二维简化分析模型的适用条件，本文相关分析结论对基于ABAQUS软件的基坑开挖变形数值分析具有一定参考价值。
1. 模型建立
为简化分析，在弹性均匀半空间内，建立长120m、宽120m、高20m的分析模型，顶面为自由边界，侧面与底面均为约束边界(图 1)。考虑边界效应对基坑开挖变形静力分析的影响，基坑开挖于分析模型地表中央，且基坑平面长20m、宽20m，设计开挖深度5m，可确保消除边界效应(管俊峰等，2017)。分析模型介质为黏性土，其中土体弹性模量取25.2MPa，粘聚力30kPa，密度1.8g/cm3，对数体积模量k取0.04，泊松比v取0.35，λ取0.2，应力比M取1.2，初始屈服面а0取0，β默认为1，流应力比K取1。初始孔隙比取1。(李广信，2004)。
分析模型采用八节点六面体(费康等，2013)单元(C3D8、全局尺寸5m×5m)进行离散，由于基坑处及附近土体变形较大，采用较小尺寸单元(C3D8、尺寸1m×1m×1m)进行离散，离散过程中采用单精度偏移措施将小尺寸单元逐渐过渡到大尺寸单元，分析模型有限元离散大样如图 2。
2. 数值模拟结果分析
根据开挖基坑实际自应力状态，采用以下分析步骤模拟整个施工过程：①首先对模型施加重力荷载，使模型在自重应力作用下稳定建立自重应力场，并清除历史上自重作用造成的位移；②采用ABAQUS软件中的单元生死功能(Model change)开挖相应土体，并进行基坑开挖变形分析。变形分析结果分别为基坑开挖引起的地表沉降变形、基坑侧壁水平位移及坑底隆起变形，地表沉降变形、基坑侧壁水平位移及坑底隆起变形沿地表路径DE、坑底路径BC、基坑侧壁路径AB给出(图 3)。
2.1 本构模型对基坑变形模拟结果的影响
不同本构模型下基坑开挖产生的地表沉降变形、基坑侧壁水平位移及坑底隆起变形分别如图 4-图 6所示。
图 4 地表路径DE沉降模拟结果(不同本构模型)Figure 4. DE settlement simulation results of ground pit surface path (Different constitutive models)
图 5 基坑侧壁路径AB水平位移模拟结果(不同本构模型)Figure 5. Simulation results of horizontal displacement of foundation pit sidewall path AB (Different constructive models)
图 6 坑底路径BC隆起模拟结果(不同本构模型)Figure 6. Simulation results of BC uplift of the bottom pit path (Different constructive models)由图 4-图 6可知，在基坑开挖未支护的情形下，Mohr-Coulomb本构模型与Drucker- Prager本构模型模拟结果基本一致，基坑侧壁地表边缘在开挖后出现隆起变形，侧壁水平变形形态为内胀腹形，坑底隆起形态呈拱形，底角处隆起量最小，中部隆起量最大，这主要因为开挖基坑周边的土体未出现塑性变形区，Mohr-Coulomb模型与Drucker-Prager模型模拟基坑开挖在弹性变形范围内的变形量相等；修正剑桥模型模拟得到的地表变形在基坑侧壁地表边缘也为隆起变形，但侧壁顶隆起量小于Mohr-Coulomb本构模型与Drucker-Prager本构模型模拟结果，最大值出现在边缘附近，基坑侧壁水平位移变化表现为向基坑倾覆，坑底隆起形态呈拱形，底角处隆起量最小，中部隆起量最大，这是因为以Mohr-Coulomb为屈服或破坏准则的理想弹塑性本构模型不能反映加卸载模量的差异及土体压硬性，由此可知，修正剑桥模型沉降区域大于Mohr-Coulomb本构模型与Drucker-Prager本构模型模拟结果，坑底隆起差异不明显，且基坑侧壁水平位移变形特征较合理，表明修正剑桥模型较Mohr-Coulomb本构模型与Drucker-Prager本构模型更适用于模拟基坑开挖变形分析。
2.2 三维模型与二维简化分析模型对比
为简便计算，实际分析中常对大型基坑问题进行简化处理，将实际的三维问题简化为二维问题。为检验此处理方法的可靠性，笔者建立长300m、宽120m、高20m的分析模型，应用ABAQUS软件进行基坑开挖变形分析时，设计的开挖宽度、深度、荷载边界条件、有限单元网格划分与前文模型一致，本构模型均采用修正剑桥模型模拟，改变基坑模型长度，使长宽比分别为1:1、3:1、5:1、7:1及10:1。变形分析结果分别为基坑开挖引起的地表沉降变形、基坑侧壁水平位移及坑底隆起变形，观测路径与前文路径相同(图 3)。长宽比为10:1的分析模型如图 7。
为简化分析，常取基坑长边中轴线截面(图 7(a)中A′B′C′)为二维简化分析模型。分析模型采用四节点四边形单元(CPE4、近似全局布种取5)进行离散，由于基坑处及附近土体变形较大，采用较小尺寸单元(CPE4、近似单元布种取1)进行离散，离散过程中采用单精度偏移措施将小尺寸单元逐渐过渡到大尺寸单元。
采用ABAQUS软件对不同长宽比的三维基坑模型及二维简化分析模型进行基坑开挖变形分析，结果如图 8-图 10。
图 8 地表路径DE沉降模拟结果(不同分析模型)Figure 8. DE settlement simulation results of ground pit surface path (Different constitutive models)
图 9 基坑侧壁路径AB水平位移模拟结果(不同分析模型)Figure 9. Simulation results of horizontal displacement of foundation pit sidewall path AB (Different constitutive models)
图 10 坑底路径BC隆起模拟结果(不同分析模型)Figure 10. Simulation results of BC uplift of the bottom pit path由图 8-图 10可知，当基坑长宽比小于5:1时，二维简化分析模型基坑变形模拟结果明显区别于三维模型模拟结果；长宽比大于等于5:1时，二维简化分析模型基坑变形模拟结果与三维模型模拟结果相近；基坑长宽比达10:1时，二维简化分析模型基坑变形模拟结果与三维模型模拟结果一致。由此可认为，当基坑长宽比不小于5:1时，可采用二维简化分析模型模拟基坑变形。
3. 结论
本文通过分析基坑开挖引起的地表沉降变形、基坑侧壁水平位移及坑底隆起变形，研究ABAQUS软件内置本构模型及二维简化分析模型对基坑开挖数值模拟结果的影响，并给出基坑开挖三维模型简化成二维模型的适用条件，主要研究结论如下：
(1) 使用Mohr-Coulomb模型和Drucker-Prager模型进行三维基坑模拟时，如果土体未发生屈服变形，这两种本构模型模拟结果一致。
(2) 修正剑桥模型沉降区域大于Mohr-Coulomb本构模型与Drucker-Prager本构模型模拟结果，坑底隆起差异不明显，且基坑侧壁水平位移变形特征较合理，表明修正剑桥模型较Mohr-Coulomb本构模型与Drucker-Prager本构模型更适用于模拟基坑开挖变形分析。
(3) 当长宽比不小于5:1时，可采用二维简化分析模型进行基坑开挖变形分析。
-
表 1 模拟的余震频次及拟合频次对比表
Table 1. Comparison of simulated with fitted frequency of aftershocks
震级M 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0 4.4 4.8 5.2 5.6 模拟频次 48964 25291 12710 6405 3303 1615 859 428 229 105 58 33 拟合频次 48261 24662 12602 6440 3291 1682 859 439 224 115 59 30 相对误差 1.4% 2.5% 0.8% 0.5% 0.4% 4.1% 0 2.6% 2.2% 9.5% 1.7% 9.1% 
下载: 导出CSV
表 2 不同b值的拟合效果
Table 2. The fitting results of different b value
初始b值 0.5555 0.6180 0.7123 0.8234 拟合b值 0.5560 0.6209 0.7084 0.8174 b值相对误差 0.09% 0.47% 0.55% 0.73% 相关系数 0.9999 0.9999 0.9994 0.9991
下载: 导出CSV
表 3 不同随机数种子的拟合效果
Table 3. The fitting results of different seeds
随机数种子 12322 350003 1234567 19491001 20080808 拟合b值 0.7214 0.7162 0.7364 0.7379 0.7278 b值相对误差 1.18% 1.89% 0.88% 1.08% 0.30% 相关系数 0.9991 0.9990 0.9994 0.9995 0.9999
下载: 导出CSV
-
傅征祥, 吕晓健, 邵辉成等, 2008.中国大陆及其分区余震序列b值的统计特征分析.地震, 28(3):1-7. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizn200803002&dbname=CJFD&dbcode=CJFQ 蒋海昆, 曲延军, 李永莉等, 2006.中国大陆中强地震余震序列的部分统计特征.地球物理学报, 49(4):1110-1117. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx200604023&dbname=CJFD&dbcode=CJFQ 刘善琪, 李永兵, 田会全等, 2013.影响b值计算误差的Monte Carlo实验研究.地震, 33(4):135-144. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizn201304016&dbname=CJFD&dbcode=CJFQ 马宏生, 刘杰, 张国民等, 2005.基于G-R关系的应变积累释放模型研究中国大陆强震的分区活动.地震学报, 27(4):355-366. http://d.wanfangdata.com.cn/Periodical_dizhen200504001.aspx 任雪梅, 高孟潭, 刘爱文等, 2009.1900年以来我国西南地区强余震统计特征.震灾防御技术, 4(2):200-208. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20090209&journal_id=zzfyjs 王菊, 王和明, 徐海龙等, 2015.任意分布随机数的FPGA实现.火力与指挥控制, 40(4):173-175, 180. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hlyz201504041&dbname=CJFD&dbcode=CJFQ 王伟锞, 李志强, 李晓丽, 2011.利用余震法快速判定宏观震中的研究.震灾防御技术, 6(1):36-48. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20110104&journal_id=zzfyjs 吴开统, 焦远碧, 吕培苓等, 1990.地震序列概论.北京:北京大学出版社. 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会, 2005.GB/T 18207.2-2005防震减灾术语第2部分:专业术语.北京:中国标准出版社. 周蕙兰, 房桂荣, 章爱娣等, 1982.余震序列的持续时间.地震学报, 4(1):45-54. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxb198201004&dbname=CJFD&dbcode=CJFQ Bath M., 1965. Lateral in homogeneities of the upper mantle. Tectonophysics, 2(6):483-514. doi: 10.1016/0040-1951(65)90003-X Eob B. C., Soo S. J., 1992. Monte-Carlo simulation of earthquake sequence in the time and magnitude space. The Journal of Engineering Geology, 2(2):147-154. http://www.dbpia.co.kr/Article/NODE00766673 Gutenberg B., Richter C. F., 1944. Frequency of earthquakes in California. Bulletin of the Seismological Society of America, 34(4):185-188. 期刊类型引用(5)
1. 孙庆，王志云，邢朋飞，冯文强，杨兴东，丁勇. 盾构机刀盘泥饼性能及机械臂攀附试验研究. 低温建筑技术. 2024(03): 121-124 . 
百度学术2. 蔡烙. 基于可靠度指标约束的斜拉桥索力多目标优化研究. 湖南交通科技. 2024(04): 99-103 . 百度学术3. 杨忠平，刘浩宇，周小涵，张益铭，常嘉卓. 深厚淤泥土深长基坑开挖对邻近建筑的影响. 地下空间与工程学报. 2022(03): 1015-1024 . 百度学术4. 马玉岩，闫飞跃，程丽娟. 不同本构模型对近接施工力学响应的影响研究. 水电站设计. 2022(04): 23-27+33 . 百度学术5. 朱晓良. 道路下管线沟槽回填加固技术模拟研究. 河南科技. 2021(04): 124-127 . 百度学术其他类型引用(6)
-
下载:
点击查看大图
计量
- 文章访问数: 90
- HTML全文浏览量: 77
- PDF下载量: 1
- 被引次数: 11
