1994年发生在台湾海峡的一次地震海啸的数值模拟

建立了一个地震海啸数值模式,模式包含越洋海啸传播部分和近岸海啸变形部分,在越洋海啸传播部分中采用线性浅水方程,使用蛙跃格式求解,并且选择合适的空间步长与时间步长,使差分格式中产生的数值频散与包辛尼斯克方程中的物理频散一致,这样在不影响海啸数值计算精度的前提下,节省了计算机的机时与内存.在近岸海啸变形部分的计算中,考虑了非线性对流项与海底摩擦项.同时该模式采用了多重网格嵌套技术,提高了所关心地区的计算精度.利用这个地震海啸模式模拟了1994年发生在台湾海峡的一次地震海啸,结果与观测记录较吻合.这个模型已用于我国沿海核电站可能最大地震海啸的数值计算.     

一次北印度洋海啸的数值模拟研究

对2012年4月11日(世界时)发生在北印度洋的一次地震海啸进行了模拟研究,给出了海啸模拟所需地震参数的获取途径,以及海啸波初始场的建立方法,分析了该海啸过程最大波幅分布和海啸波传播波时分布,并与测站实测值进行了对比。结果表明,海啸波的能量传播方向存在不均匀性,能量传播集中方向的波幅较大,具有更强的危害性。     

一个基于GPU并行加速的海啸数值模型

地震海啸通常发生在大洋板块向陆地板块俯冲的区域,距离震源最近的国家和地区往往在震后5～20 min之内就会遭受到海啸袭击。因此,及时的海啸预警和准确的海啸预报结果对于民众和决策者都至关重要。为了提升海啸预警效率,缩短海啸预报时间,本研究对COMCOT海啸数值模型进行了基于图形计算单元GPU的二次并行开发。将原模型中海啸传播计算模块通过CUDA＿C语言编写内核函数整体移植到GPU上并行加速,CPU负责模型其他代码的执行。为了减少CPU和GPU之间的数据通信,将吸收边界和变量更新函数一并改写。仅在需要输出的时间节点,GPU向CPU传递结果,其他时间步长,CPU和GPU之间只有指令和少量参数传输,基本可视为零耗时。基于GPU并行加速的COMCOT较串行版本效率提升超过67倍,加速性能显著优于基于CPU共享内存的OpenMP并行版本。交叉使用常水深和真实地形,采用均匀滑移海啸源和有限元海啸源对模型的计算结果进行了较为全面的分析检验,相对误差最大不超过1%,为大范围的越洋海啸实时计算提供了有力工具。     

台湾海峡一次海啸的初步数值模拟

应用一个三维斜压陆架海模式--HAMSOM模式,首先对台湾海峡内的背景潮汐场进行了数值模拟,随后将一个参数化公式作为强迫条件,对1994年发生在海峡内的一次海啸进行模拟,结果与实测数据比较吻合.还分析了海啸波在海峡沿岸的分布情况及其对沿岸区域的影响状况,结果表明该次海啸波动对海峡沿岸区域的影响不大.     

基于数值模拟的渤海海域地震海啸危险性定量化研究

根据地震海啸产生的条件,结合渤海海域的地形特征、地质构造、地震学特征和历史地震及海啸记录对渤海海域潜在的地震海啸进行了数值模拟研究。分析了渤海可能引发地震海啸的震源区域,讨论了渤海发生海啸灾害的可能性。文中通过数值模拟再现了渤海历史上几次规模较大的地震事件可能引发的海啸情景,研究分析了可能的地震海啸在渤海及周边海域的传播过程及波动特征.地震海啸传播模型采用基于四叉树原理的自适应网格加密技术,有效解决了局部分辨率与计算效率之间的矛盾。数值计算包括地震海啸产生及传播过程。利用该模型对渤海潜在的地震海啸进行了数值计算,基于数值计算结果定量阐述了渤海海域潜在地震海啸对渤海局部岸段及北黄海沿岸的影响,给出了渤海可能地震海啸危险性划分;研究结果将为我国海啸危险性分析和海啸预警技术研究工作提供技术支持。     

河北省海啸灾害风险评估和区划研究

基于康奈尔多电网耦合海啸模型（COMCOT）数值模式建立研究区海啸数值模型,利用模型对6个地震海啸源情景场进行模拟计算,通过对各海啸源下数值模拟结果的综合分析,确定研究区海啸灾害危险性分布;根据土地利用类型数据资料确定研究区承灾体脆弱性分布;通过危险性、脆弱性评估成果,综合确定研究区海啸灾害风险性评估和区划成果。结果表明：河北省遭受海啸灾害风险最高等级为Ⅱ级,主要分布在唐山市和秦皇岛市沿海区域;沧州市沿海区域海啸灾害风险等级较小,为Ⅲ级和Ⅳ级。     

深圳海域潮汐海啸波耦合数值研究

以COMCOT海啸模式和TPXO7.1全球潮汐模式为基础,采用三层嵌套网格,建立了南海海啸与潮汐耦合计算模型,分析深圳海域海啸和潮汐相互作用。潮汐计算结果与实测数据吻合较好,高、低潮位平均误差小于15 cm,20 cm;在潮汐验证的基础上,以马尼拉海沟潜在地震海啸源为案例,进行8.0,9.0级地震海啸与潮汐耦合情景模拟计算,计算结果表明,9级地震海啸在深圳海域外海波高为140~150 cm,如先行波为正波发生在高潮时将产生异常高潮位,负波发生在低潮时将产生异常低潮位,线性叠加计算结果偏大,在25.0 cm之内,到达时间差异小于6 min。     

2017年9月8日墨西哥沿岸Mw8.2级地震海啸观测数据分析与模拟

2017年9月8日4时49分（UTC）,墨西哥瓦哈卡州沿岸海域（15.21°N,93.64°W）发生M_w8.2级地震,震源深度30 km。强震在该海域引发海啸,海啸对震源附近数百千米范围内造成了严重影响。位于太平洋上的多个海啸监测网络捕捉到了海啸信号并详细记录了此次海啸的传播过程。本文选用了近场2个DART浮标和6个验潮站的水位数据,通过潮汐调和分析和滤波分离出海啸信号,对近场海啸特征值进行了统计分析,并采用小波变换分析方法进一步分析了海啸的波频特征。基于Okada弹性位错理论断层模型计算得到了强震引发的海底形变分布,并采用MOST海啸模式对本次海啸事件近场传播特征进行了模拟,模拟结果与观测吻合较好。最后,基于实测和模拟结果,详细分析了此次地震海啸的近场分布特征,发现除受海啸源的强度和几何分布特征影响外,近岸海啸波还主要受地形特征控制,在与特定地形相互作用后波幅产生放大效应,会进一步加剧海啸造成的灾害。     

有限断层模型在2015年9月16日智利Mw8.3级地震海啸数值模拟中的应用与评估

2015年9月16日22时54分（当地时间）智利中部近岸发生M_w8.3级地震,震源深度25 km。同时,强震的破裂区长200 km,宽100 km,随之产生了中等强度的越洋海啸。海啸影响了智利沿岸近700 km的区域,局部地区监测到近5 m的海啸波幅和超过13 m的海啸爬坡高度。太平洋区域的40多个海啸浮标及200多个近岸潮位观测站详细记录了此次海啸的越洋传播过程,为详细研究此次海啸近场及远场传播及演化规律提供了珍贵的数据。本文选择有限断层模型和自适应网格海啸数值模型建立了既可以兼顾越洋海啸的计算效率又可以实现近场海啸精细化模拟的高分辨率海啸模型。模拟对比分析了海啸的越洋传播特征,结果表明采用所建立的模型可以较好地再现远场及近场海啸特征,特别是对近场海啸的模拟结果非常理想。表明有限断层可以较好地约束近场、特别是局部区域的破裂特征,可为海啸预警提供更加精确的震源信息,结合高分辨率的海啸数值预报模式实现海啸传播特征的精细化预报。本文结合观测数据与数值模拟结果初步分析了海啸波的频散特征及其对模型结果的影响。同时对观测中典型的海啸波特征进行的简要的总结。谱分析结果表明海啸波的能量主要分布在10~50 min周期域内。这些波特征提取是现行海啸预警信息中未涉及,但又十分重要的预警参数。进一步对这些波动特征的详细研究将为海啸预警信息及预警产品的完善提供技术支撑。     

精细化地震海啸波流实时预警系统研究与应用

基于Okada模型和非线性浅水波模型,结合高精度多层嵌套网格针对我国浙江沿海的温州和台州地区建立了越洋–近海–局部的精细化地震海啸波流实时预警系统,近岸的分辨率为900 m。该预警系统包括了并行化的数值计算模块,基于Python 2D绘图库的计算结果可视化处理模块,以及通过Python语言将所有经过数值计算的图形与动画产品集成在一个网页上的产品集成模块。一旦地震发生,该系统可根据地震的震源参数信息在10 min内完成数值计算、可视化处理,以及产品集成。选取2011年日本东北9.0级地震海啸结合实测数值对该系统进行模拟验证,进一步应用该系统模拟计算了日本南海海槽和琉球海沟潜在极端海啸的影响规律。结果表明,该预警系统可有效地提高地震海啸实时预警的时效性和准确度,为海啸的预警、减灾,以及辅助决策提供科学依据。     

越洋海啸波传播数学模型的研究

应用内外解匹配的方法和驻相法推导了柱坐标系下地震引起的水面波动方程的解,即推导了地震海啸波生成与传播的理论方程,并对越洋地震海啸理论方程进行了求解。基于Boussinesq方程出发,建立了二维越洋海啸波传播数学模型,并对越洋海啸进行了数值计算,计算方法采用有限差分方法,差分格式采用交替方向隐格式(Alternating Direction Implicit即ADI方法)。利用越洋海啸计算模式对发生在大洋中的地震海啸进行了模拟,将数值模拟结果与地震海啸波理论方程的计算结果进行了比较,两种计算结果吻合较好。     

温州瓯江口浅滩地区越洋海啸影响评估计算

日本“3·11”地震海啸事件发生后,为了避免灾难重演,各滨海国家在加强海啸基础理论研究、改进海啸预警系统的同时,还应对沿海重大工程及重点保障目标进行地震海啸灾害风险排查及再评估工作;对在建的重大基础设施和社会经济功能区划应进行全面的地震海啸安全论证.在此背景下,该文首先概括总结了我国东南沿海的地震海啸风险及历史海啸事件时空分布,简要介绍了越洋海啸传播特征.海啸源选取基于潜在可能最大海啸,选取环太平洋地震带上的潜在地震海啸源,进行温州瓯江口地区越洋海啸影响评估计算.海啸数值计算模型采用美国康奈尔大学的COMCOT模型,利用该模型对2010年智利海啸、2011年日本海啸进行了近场、远场模拟验证,计算结果与观测数据吻合良好,模型可信.应用联合国教科文组织政府间海洋学委员会(UNESCO/IOC)太平洋海啸预警系统的海啸危险性等级标准,结合评估计算结果,对瓯江口浅滩地区海啸危险性进行等级划分,获得了该地区的海啸危险性初步评估结果.结果表明:在所选的10个潜在或历史海啸源产生的越洋海啸对研究区域的影响均小于100 cm,此规模的海啸不易对该地区造成灾害性影响.研究结果对于指导该地区的海啸灾害风险评估及风险排查具有一定的参考价值.     

局地海啸波传播数学模型的研究

推导了地震海啸波传播的理论方程,并时局地地震海啸情况下的理论方程进行了求解.基于Boussinesq方程出发,建立了二维局地海啸渡传播数学模型;对局地海啸进行了数值计算,计算方法采用有限差分方法.差分格式采用交替方向隐格式,即ADI方法.利用该模型对发生在台湾海峡的一次地震海啸进行了模拟计算.将计算结果与理论方程的计算结果以及实际情况进行了对比.计算结果较好.为局地海啸波传播提供了一种模拟方法.     

日本南海海槽地震海啸对我国东海沿岸影响

本文利用数值模拟技术重现了1707年宝永地震海啸的传播过程,定量分析了我国东海沿岸海啸时空分布特征。计算结果表明,地震发生2.5小时后海啸波传至东海陆架,震后6小时浙江沿海地区遭到海啸的袭击,沿岸最大海啸波高为0.8米。通过海啸波在东海大陆架传播时海底地形与波幅的关系,研究分析了东海陆架缓变地形下海啸放大效应,为及时判断沿海可能的海啸强度和受灾程度提供了便捷的估算方法。此外,本文还评估了南海海槽发生极端地震时,中国东海沿岸的海啸危险性,为东海区域针对日本南海海槽进行海啸预警和减灾评估提供定量科学的参考。