摘要:建立简化的岛礁-潟湖-海水系统地震反应分析模型,重点分析考虑潟湖影响下岛礁场地的地震反应规律。研究结果表明,潟湖对岛礁临湖侧水平向地震反应的放大作用较为显著,考虑潟湖影响时临湖侧角点的峰值加速度是不考虑潟湖影响时该点峰值加速度的138.6%;随着潟湖水深的增加,临湖侧的地震反应进一步增大;入射波经岛礁-潟湖-海水系统传至礁坪时,中周期的频率成分被过滤,反应谱峰值有向短周期方向移动的趋势。本文研究结果可为进一步确定岛礁场地地震参数提供参考。
加入收藏
南海岛礁位于环太平洋地震带西侧,岛上建筑面临严峻的地震灾害风险。作为岛礁工程抗震设计的基础,岛礁场地的地震反应问题逐渐受到重视。胡进军等探讨南海岛礁场地地震安全性研究涉及的关键问题[1],并分别利用一维土层反应分析模型[2]和浅部岛礁工程地质二维剖面模型[3],进行地震波垂直入射下岛礁地震反应研究;陈国兴等[4]考虑珊瑚砂的动力非线性特性和近场截断处的人工边界条件,利用珊瑚岛礁浅部二维地震反应分析模型,研究岛礁场地地震峰值加速度放大规律和地表加速度反应谱等。刘晶波等开发一种可模拟无限海水层波动辐射效应的流体介质人工边界[5],并提出适用于岛礁-海水系统地震波输入的人工边界子结构法[6],建立岛礁-海水系统地震反应分析模型[7],在验证其准确性的基础上研究P波和SV波垂直入射下典型岛礁场地的地震反应规律[8,9]。但是,以往研究中均未考虑潟湖对岛礁场地地震反应的影响。根据文献[10],多数南海岛礁顶部均存在开放或封闭式的潟湖,潟湖水深一般为5~30m。在地震作用下,岛礁地震反应可能受潟湖地形和湖水动力耦合作用的影响,有必要对其进行系统研究。
鉴于此,在前人研究的基础上,建立简化的岛礁-潟湖-海水系统地震反应分析模型,分析其在SV波垂直入射下的地震反应,并评价潟湖对岛礁场地地震反应的影响。
1、模型与方法
1.1无限域流体和固体介质模拟
结合文献[10,11,12,13],建立局部成层的岛礁-潟湖-海水系统近场有限元模型,见图1。岛礁礁体介质的材料和几何参数见表1,潟湖水深为15m,宽度为2400m。固体与流体介质分别采用通用有限元软件ANSYS中的SOLID45和FLUID80单元模拟,通过耦合流-固界面的法向自由度模拟岛礁-海水和岛礁-潟湖之间的相互作用[14]。此外,作为一种外源激励下的开放系统,除流-固耦合效应外,岛礁场地的地震反应模拟还涉及2个关键环节:一是无限域流体和固体介质波动辐射效应的模拟,二是地震波输入。刘晶波等在前期研究中提出一种空间解耦的流体介质人工边界条件[5],并验证其在岛礁-海水系统地震反应分析中的适用性和计算精度[7]。本文采用该流体介质人工边界模拟无限海水层的波动辐射效应,利用文献[5]推荐的基于ANSYS的COMBIN40单元建立模型。同时,采用一致黏弹性人工边界单元模拟无限海床的波动辐射效应,流、固介质人工边界参数的确定原则与文献[8,9]一致。
图1岛礁-潟湖-海水系统近场有限元模型
1.2地震波输入方法
地震波是一种外源载荷,将其合理有效地输入计算系统是决定岛礁场地地震反应计算精度的关键环节。目前,地震波输入一般是将入射地震波转化为等效输入地震载荷[15,16,17],从而避免地震波输入过程影响人工边界条件对外行散射波的吸收作用。本文采用文献[6]提出的人工边界子结构法进行岛礁-潟湖-海水系统的地震波输入,实施过程如下。
(1)参照文献[9],在岛礁-海水系统有限元模型中截取相对应的人工边界子结构模型,见图2。
(2)基于行波理论求解与入射波对应的自由波场。对人工边界子结构模型的节点分别施加相应的自由波场位移时程,并对该子结构模型进行动力分析,求得人工边界节点的反力,即为实现地震波输入的等效输入地震载荷。
(3)对岛礁-潟湖-海水系统有限元模型的人工边界节点施加由步骤(2)获得的等效输入地震载荷,完成岛礁场地的地震波输入,通过动力计算获得岛礁-潟湖-海水系统的地震反应。
图2地震波输入的人工边界子结构模型
选用Kobe波、Koceali波、LomaPrieta波和Northridge波(见图3)作为输入地震SV波,并将以上地震的峰值加速度调幅至1m/s2,此时与其相对应的平坦场地自由地表峰值加速度为2m/s2。采用人工边界子结构法将上述地震SV波垂直输入至岛礁-潟湖-海水系统有限元模型。
图3输入地震波的加速度时程和傅里叶幅值谱
2、模拟结果与分析
2.1加速度峰值放大效应
研究考虑潟湖影响的岛礁场地地震加速度峰值放大效应。建立几何和材料参数与原模型相同但不包含潟湖的岛礁-海水系统有限元模型作为对照,比较不同地震SV波垂直入射下岛礁顶部的峰值加速度amax。计算模型和入射地震波均空间对称,因此仅需考察岛礁模型右端由临湖侧至临海侧地震波峰值加速度的空间分布情况,岛礁场地加速度的空间分布见图4。在地震SV波垂直入射下,岛礁顶部水平方向地震反应由临湖侧至临海侧呈先减小后增大的趋势,最大值出现在临湖侧或临海侧角点上,此时峰值加速度最大可达7.6m/s2;岛礁顶部竖向地震反应相对较小,分布也较为平缓,加速度峰值位于1.0m/s2附近。潟湖对岛礁临湖侧水平向地震反应的放大作用较为显著,考虑潟湖影响时临湖侧角点的峰值加速度是不考虑潟湖影响时该点峰值加速度的138.6%。
图4岛礁场地峰值加速度放大的空间分布
2.2地表加速度反应谱分析
取测点A、B和C分别为岛礁顶部右侧礁坪的临湖角点、中点和临海角点。计算测点A、B和C的地表加速度反应谱,与相同地震波入射下平坦场地自由地表的加速度反应谱进行对比,结果见图5,其中β为动力因数。由此可见,岛礁顶部加速度反应谱均表现出单峰特征,由临湖侧至临海侧反应谱峰值呈增大趋势,且与平坦场地自由地表相比,反应谱峰值有向短周期方向移动的趋势,震动周期为0.2~0.3s时的谱加速度反应尤为显著,而当周期为0.4~1.7s时岛礁场地的谱加速度小于平坦场地自由地表,说明入射波经岛礁-海水系统传至礁坪时,中周期的频率成分有所减少。
图5岛礁场地的地表加速度反应谱
2.3潟湖深度对岛礁场地地震反应的影响
鉴于南海岛礁的潟湖水深存在较大差异[10],进一步研究潟湖水深对岛礁场地地震反应的影响。保持模型中其他参数不变,改变潟湖水深,令其分别为10、15、20和25m,比较不同模型计算得到的礁顶峰值加速度的空间分布,结果见图6。由此可见,潟湖水深对岛礁场地地震反应的影响仅限于临湖侧角点附近,影响区域的宽度一般不超过100m。随着潟湖水深的增加,临湖侧的地震反应增大,LomaPrieta波入射下该趋势尤为明显;当潟湖水深达25m时,临湖侧角点的加速度峰值最大可达9.0m/s2。
图6潟湖水深度对加速度的影响
3、结论
基于前期的研究基础,建立简化的岛礁-潟湖-海水系统地震反应分析模型,重点研究考虑潟湖影响下岛礁场地的地震加速度放大效应和地表加速度反应谱特性,得到以下结论。
(1)在地震SV波垂直入射下,岛礁顶部水平向地震反应由临湖侧至临海侧呈先减小后增大的趋势,潟湖对岛礁临湖侧水平向地震反应的放大作用较为显著,考虑潟湖影响时的临湖侧角点处的峰值加速度是不考虑潟湖影响时的138.6%。
(2)入射波经岛礁-海水系统传递至礁坪时,中周期的频率成分有所减少,礁坪处的反应谱峰值有向短周期方向移动的趋势。震动周期为0.2~0.3s时,谱加速度反应尤为显著。
(3)潟湖深度对岛礁场地地震反应的影响仅限于临湖侧角点附近,对于本文模型,影响区域的宽度一般不超过100m。随着潟湖水深的增加,临湖侧的地震反应增大。
参考文献:
[1]胡进军,徐长琦,谢礼立,等.南海岛礁场地地震稳定性研究中的关键问题探讨[J].地震工程学报,2018,40(2):92-100.
[2]胡进军,李天男,谢礼立,等.脉冲型地震动作用下典型珊瑚岛礁的场地放大研究[J].世界地震工程,2017,33(4):1-10.
[3]胡进军,徐长琦,李琼林,等.典型岛礁场地的地震效应初探[J].地震工程与工程振动,2018,36(6):18-25.
[4]陈国兴,朱翔,赵丁凤,等.珊瑚岛礁场地非线性地震反应特征分析[J].岩土工程学报,2019,41(3):405-413.
[5]刘晶波,宝鑫,谭辉,等.波动问题中流体介质的动力人工边界[J].力学学报,2017,49(6):1418-1427.
[6]刘晶波,谭辉,宝鑫,等.土-结构动力相互作用分析中基于人工边界子结构的地震波动输入方法[J].力学学报,2018,50(1):32-43.
[7]宝鑫,刘晶波,韩鹏飞,等.流体人工边界在海域场地地震分析中的应用[J].地震工程与工程震动,2018,38(S2):95-105..
[8]宝鑫,刘晶波,王东洋,等.P波垂直入射下海域岛礁场地动力反应分析[J].工程力学,2019,36(S1):1-7.
[10]王新志.南沙群岛珊瑚礁工程地质特性及大型工程建设可行性研究[D].武汉:中国科学院研究生院(武汉岩土力学研究所),2008.
[11]詹文欢,刘以宣,钟建强,等.南海南部活动断裂与灾害性地质初步研究[J].海洋地质与第四纪地质,1995,15(3):1-9.
[12]孙宗勋,卢博.南沙群岛珊瑚礁灰岩弹性波性质的研究[J].工程地质学报,1999,7(2):175-180.
[13]孙宗勋,詹文欢,朱俊江.南沙群岛永暑礁岩体稳定性评价[J].海洋地质与第四纪地质,2003,23(4):9-14.
[14]宝鑫,刘晶波.考虑流-固耦合效应的含液容器动力响应有限元分析方法[J].核动力工程,2017,38(2):111-114.
[15]刘晶波,吕彦东.结构-地基动力相互作用问题分析的一种直接方法[J].土木工程学报,1998,31(3):55-64.
刘晶波,宝鑫,李述涛,赵启明,王菲,王东洋.考虑潟湖影响的岛礁场地地震反应分析[J].计算机辅助工程,2020,29(02):46-50+69.
基金:国家自然科学基金(51878384);国家重点研发计划(2018YFC1504305).
分享:
随着社会经济不断发展,社会化进程不断加快,国际油田勘探工程也得到了较快的发展,全球工业生产,社会经济的发展对于石油的需求量越来越大,石油产业发展的相关技术研究也得到了业界的关注。目前,国际油田开采存在受到地质环境限制,石油勘探及开采项目开展会受到地质环境的限制,增加了石油勘探的难度。
2024-08-31目前我国地下金属矿山的开采正处于“由浅至深、由易至难、由富至贫”的关键转型期,在一个矿区内采用单一采矿方法进行回采的传统模式已难以实现矿产资源安全、经济、高效、充分回收的目标。近年来,协同开采理念的提出为化解这一难题提供了一条重要途径该理念的核心在于将不同采矿方法在同一矿区内进行最大程度整合。
2024-06-25随着全球化进程的加快和社会的不断发展,工业、军事等领域对矿产资源需求不断增大,但我国主要资源开采带的矿产储量逐渐减少,亟需对我国复杂地形的矿产资源进行勘探开采,因此需要更高性能和更加灵活的地层勘探仪器来探明地层结构,从而保障资源开采过程中人员和设备的安全。
2024-01-03随着全球化进程的加快和社会的不断发展,工业、军事等领域对矿产资源需求不断增大,但我国主要资源开采带的矿产储量逐渐减少,亟需对我国复杂地形的矿产资源进行勘探开采,因此需要更高性能和更加灵活的地层勘探仪器来探明地层结构,从而保障资源开采过程中人员和设备的安全。
2024-01-03油气管道是目前最重要的能源输送方式之一,我国陆上油气管道总里程长、分布广阔、沿线地质环境复杂多变。由于各方面的考虑,管道可能敷设在滑坡地带,也有可能在管道敷设后由于工程活动或者地址变动等种种原因导致沿线发生滑坡。而滑坡灾害对管道危害巨大,能够挤压管道导致管道变形甚至破坏,进而导致管道输送的油气泄露,
2023-12-29由于我国采矿技术的提升,以及煤矿开采范围的扩大,部分矿业公司在经营管理与生产活动中将面临着很多的安全难题,比如对安全技术的了解不足、对安全方面的资金支持欠缺等,使采矿过程中地质灾害问题频频出现。不论是矿业企业的管理者还是其他管理人员,都应当提升采矿过程中的安全生产能力,积极采取正确可行的安全措施,减少矿井地质灾害的风险,维护矿井职工的安全,确保我国矿业的可持续、安全、平稳的发展。
2023-12-06据中国地震台网正式测定,当地时间2023-02-06,土耳其发生4次地震,其中2次震级高达7.8。研究表明,地震前后大气中的水汽含量会出现异常变化[1,2,3,4,5]。大气可降水总量(total precipitable water, TPW)是表征大气中水汽含量的重要指标,而卫星遥感是反演TPW的重要手段之一[6,7]。
2023-11-29天水市是典型的山间河谷盆地型城市。区内黄土分布广,厚度大,垂直节理和裂隙发育,降雨集中分布,致使沟壑纵横,地貌地形支离破碎,地质环境十分脆弱[1,2,3],在自然因素以及人类工程活动影响下,滑坡、泥石流、崩塌、不稳定斜坡失稳等地质灾害现象频繁发生[4]。
2023-11-08在工程建设和岩土施工过程中,地质灾害尤其是滑坡问题时有发生,这主要是由于边坡的稳定性受到了岩土层变化的影响,因此,提高边坡的稳定性对预防滑坡问题至关重要。随着科学技术的不断发展,我国岩土施工技术也在不断地改进中,岩土施工技术的合理运用,对有效治理地质灾害起决定性作用,要求相关技术人员必须认识到科学合理应用其施工技术的重要意义。
2023-11-08乐矿能源(原乐平矿务局,简称乐矿)地处赣东北景德镇市境内,其下属的煤矿均分布在萍乐坳陷东北端,多在20世纪50、60年代建设。乐平矿区所属煤矿煤炭资源分布不均,开采条件相对较差,矿井规模都比较小,而且煤层瓦斯含量高、灾害严重,矿区矿井水文地质条件也较复杂。
2023-11-08我要评论
期刊名称:地学前缘
期刊人气:4058
主管单位:中华人民共和国教育部
主办单位:中国地质大学(北京),北京大学
出版地方:北京
专业分类:地质
国际刊号:1005-2321
国内刊号:11-3370/P
创刊时间:1994年
发行周期:双月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:2.546
影响因子:0.659
影响因子:0.149
影响因子:1.209
影响因子:2.866
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!