海洋学报

断层参数	参数值
震源位置	31.57°S, 71.67°W
震源深度	22.4 km
断层长度	212 km
断层宽度	79 km
倾角	19°
滑动角	83°
走向角	353°
平均滑移量	6.3 m

断层参数	参数值
震源位置	31.57°S, 71.67°W
震源深度	22.4 km
断层长度	212 km
断层宽度	79 km
倾角	19°
滑动角	83°
走向角	353°
平均滑移量	6.3 m

测站名称	位置	到达时间		先导波波幅
COQUIMBO CL	30.0°S,71.3°W	28	32	12.5%		104.0	96.9	6.8%
SAN ANTONIO CL	33.6°S,71.6°W	30	38	21.1%	73.4	65.3	11.0%
BUCALEMU CL	34.6°S,72.0°W	41	49	16.3%	48.5	46.9	3.3%
HUASCO CL	28.5°S,71.2°W	30	35	14.3%	46.9	46.2	1.5%
CONSTITUCION CL	35.4°S,72.5°W	53	63	15.9%	68.6	49.7	27.6%
CHANARAL CL	26.4°S,70.6°W	57	61	6.6%	41.8	42.5	1.6%
TALCAHUANO CL	36.7°S,73.1°W	104	108	3.7%	43.4	44.0	1.4%
TALTAL CL	25.4°S,70.5°W	67	67	0%	18.8	19.5	3.6%
JUAN FERNANDEZ	33.6°S,78.8°W	63	60	4.8%	45.4	40.9	9.9%
CORRAL CL	39.9°S,73.4°W	112	132	15.2%	26.1	24.4	6.5%
TOCOPILLA CL	22.1°S,70.2°W	93	92	1.1%	12.0	10.8	10.0%
SAN FELIX CL	26.3°S,80.1°W	80	81	1.2%	23.8	22.6	5.0%
MATARANI PE	17.0°S,72.1°W	141	140	0.7%	13.3	13.9	4.3%
CALLAO LA-PUNTA PE	12.1°S,77.2°W	206	218	5.5%	23.1	20.3	12.1%

测站名称	位置	到达时间		先导波波幅
COQUIMBO CL	30.0°S,71.3°W	28	32	12.5%		104.0	96.9	6.8%
SAN ANTONIO CL	33.6°S,71.6°W	30	38	21.1%	73.4	65.3	11.0%
BUCALEMU CL	34.6°S,72.0°W	41	49	16.3%	48.5	46.9	3.3%
HUASCO CL	28.5°S,71.2°W	30	35	14.3%	46.9	46.2	1.5%
CONSTITUCION CL	35.4°S,72.5°W	53	63	15.9%	68.6	49.7	27.6%
CHANARAL CL	26.4°S,70.6°W	57	61	6.6%	41.8	42.5	1.6%
TALCAHUANO CL	36.7°S,73.1°W	104	108	3.7%	43.4	44.0	1.4%
TALTAL CL	25.4°S,70.5°W	67	67	0%	18.8	19.5	3.6%
JUAN FERNANDEZ	33.6°S,78.8°W	63	60	4.8%	45.4	40.9	9.9%
CORRAL CL	39.9°S,73.4°W	112	132	15.2%	26.1	24.4	6.5%
TOCOPILLA CL	22.1°S,70.2°W	93	92	1.1%	12.0	10.8	10.0%
SAN FELIX CL	26.3°S,80.1°W	80	81	1.2%	23.8	22.6	5.0%
MATARANI PE	17.0°S,72.1°W	141	140	0.7%	13.3	13.9	4.3%
CALLAO LA-PUNTA PE	12.1°S,77.2°W	206	218	5.5%	23.1	20.3	12.1%

浮标号	位置	到达时间		先导波波幅
32402	26.7°S,74.0°W	46	39	15.2%		9.3	10.0	7.0%
32401	20.5°S,73.4°W	103	96	6.8%	5.5	5.1	7.3%
32412	18.0°S,86.4°W	174	171	1.7%	6.3	6.1	3.2%
32411	5.0°N,90.9°W	398	394	1.0%	2.5	1.9	24.0%
43413	11.0°N,100.1°W	476	473	0.6%	1.9	1.4	26.3%
43412	16.0°N,107.0°W	560	560	0%	2.5	1.9	24.0%
46411	39.3°N,127.1°W	837	839	0.2%	1.3	0.8	38.5%
46407	42.7°N,127.8°W	874	876	0.2%	1.1	1.3	15.4%
46404	45.8°N,128.8°W	910	910	0%	1.0	0.7	30.0%
46409	55.3°N,148.5°W	1 030	1 034	0.4%	1.2	1.8	33.3%
46403	52.7°N,157.0°W	1 032	1 038	0.6%	1.4	1.9	26.3%
46408	49.6°N,169.9°W	1 063	1 073	0.9%	1.6	2.2	27.3%
46413	48.0°N,174.2°W	1 074	1 082	0.7%	1.3	1.8	27.8%
51407	19.6°N,156.5°W	872	878	0.7%	1.8	1.9	5.3%
51425	9.5°S,176.3°W	879	891	1.3%	1.4	0.8	42.9%
52406	5.3°S,165.0°E	1 075	1 091	1.5%	0.3	0.4	25.0%
52401	19.3°N,155.8°E	1 220	1 238	1.5%	1.3	1.0	23.1%
52402	11.9°N,153.9°E	1 224	1 242	1.4%	1.1	0.7	36.4%
52403	4.0°N,145.5°E	1 300	1 318	1.4%	0.8	0.3	62.5%
21414	49.0°N,178.2°E	1 113	1 121	0.7%	1.1	1.4	21.4%

浮标号	位置	到达时间		先导波波幅
32402	26.7°S,74.0°W	46	39	15.2%		9.3	10.0	7.0%
32401	20.5°S,73.4°W	103	96	6.8%	5.5	5.1	7.3%
32412	18.0°S,86.4°W	174	171	1.7%	6.3	6.1	3.2%
32411	5.0°N,90.9°W	398	394	1.0%	2.5	1.9	24.0%
43413	11.0°N,100.1°W	476	473	0.6%	1.9	1.4	26.3%
43412	16.0°N,107.0°W	560	560	0%	2.5	1.9	24.0%
46411	39.3°N,127.1°W	837	839	0.2%	1.3	0.8	38.5%
46407	42.7°N,127.8°W	874	876	0.2%	1.1	1.3	15.4%
46404	45.8°N,128.8°W	910	910	0%	1.0	0.7	30.0%
46409	55.3°N,148.5°W	1 030	1 034	0.4%	1.2	1.8	33.3%
46403	52.7°N,157.0°W	1 032	1 038	0.6%	1.4	1.9	26.3%
46408	49.6°N,169.9°W	1 063	1 073	0.9%	1.6	2.2	27.3%
46413	48.0°N,174.2°W	1 074	1 082	0.7%	1.3	1.8	27.8%
51407	19.6°N,156.5°W	872	878	0.7%	1.8	1.9	5.3%
51425	9.5°S,176.3°W	879	891	1.3%	1.4	0.8	42.9%
52406	5.3°S,165.0°E	1 075	1 091	1.5%	0.3	0.4	25.0%
52401	19.3°N,155.8°E	1 220	1 238	1.5%	1.3	1.0	23.1%
52402	11.9°N,153.9°E	1 224	1 242	1.4%	1.1	0.7	36.4%
52403	4.0°N,145.5°E	1 300	1 318	1.4%	0.8	0.3	62.5%
21414	49.0°N,178.2°E	1 113	1 121	0.7%	1.1	1.4	21.4%

基于全和谐型浅水方程的有限体积海啸数值模型开发与应用

PDF下载

周文 ¹ , 王培涛 ² , 王岗 ¹^,^* , 于福江 ² , 郑金海 ¹ , 梁秋华 ³

海洋学报 | 论文 2021,43(5): 27-37

收起

海洋学报 | 论文 2021, 43(5): 27-37

基于全和谐型浅水方程的有限体积海啸数值模型开发与应用

全屏

周文¹, 王培涛², 王岗¹^,^*, 于福江², 郑金海¹, 梁秋华³

作者信息

¹海岸灾害及防护教育部重点实验室（河海大学），江苏南京 210098

²国家海洋环境预报中心，北京 100081

³拉夫堡大学建筑与土木工程学院，英国莱斯特郡 LE11 3TU

周文（1995-），男，江苏省靖江市人，主要从事海啸数值模型开发研究。E-mail：wenzhou_hhu@hhu.edu.cn

通讯作者:

王岗（1982-），男，河北省张家口市人，博士，副教授，主要从事港湾共振、水波模拟及海啸与洪水风险评估研究。E-mail：gangwang@hhu.edu.cn

Development and application of a finite volume tsunami numerical model based on the well-balanced shallow water equations

Wen Zhou¹, Peitao Wang², Gang Wang¹^,^*, Fujiang Yu², Jinhai Zheng¹, Qiuhua Liang³

Affiliations

¹Key Laboratory of Coastal Disaster and Protection (Hohai University), Ministry of Education, Nanjing 210098, China

²National Marine Environmental Forecasting Center, Beijing 100081, China

³School of Architecture, Building and Civil Engineering, Loughborough University, Leicestershire LE11 3TU, UK

出版时间: 2021-05-25 doi: 10.12284/hyxb2021095

文章导航

摘要

收起

数值模拟作为海啸预报的主要研究方法在海啸预警中起着关键作用。本文采用Godunov格式的有限体积方法，使用MUSCL-Hancock格式，并利用HLLC Riemann近似求解器计算单元界面上的流体通量，建立了球坐标系下二阶精度的海啸数值模型。模型所基于的全和谐型浅水方程保证了数值的稳定性，而地形重构方法实现了干湿边界的精准模拟。本文模拟了2015年9月16日智利M_w8.3级地震海啸，通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较，验证了模型对实际越洋海啸模拟预报的能力。

关键词

海啸数值模型 / 全和谐型浅水方程 / 有限体积法 / 干湿边界 / 智利海啸

Abstract

收起

Numerical simulation, as the major research method of tsunami, plays a key role in tsunami warning. The present paper develops a second-order accuracy numerical tsunami model in the spherical coordinate using the Godunov-type finite volume method and MUSCL-Hancock scheme. An HLLC approximate Riemann solver is employed to evaluate fluxes across cell interfaces. The well-balanced expression format of shallow water equations ensures the numerical stability, while the local topography reconstruction method is used to deal with the moving shoreline boundary. The model is used to investigate the propagation of the 16 September 2015 Chile tsunami. The capability of the model is verified by comparison with the observational data from 14 coastal tidal-gauge stations near Chile and 20 DART buoys covering the Pacific Ocean.

Key words

tsunami numerical model / well-balanced shallow water equations / finite volume method / moving shoreline boundary / Chile tsunami

引用本文

周文, 王培涛, 王岗, 于福江, 郑金海, 梁秋华. 基于全和谐型浅水方程的有限体积海啸数值模型开发与应用. 海洋学报, 2021 , 43 (5) : 27 -37 . DOI: 10.12284/hyxb2021095

Wen Zhou, Peitao Wang, Gang Wang, Fujiang Yu, Jinhai Zheng, Qiuhua Liang. Development and application of a finite volume tsunami numerical model based on the well-balanced shallow water equations[J]. Haiyang Xuebao, 2021 , 43 (5) : 27 -37 . DOI: 10.12284/hyxb2021095

正文

收起

1 引言

收起

海啸作为最具危险性和破坏力的海洋灾害之一，给沿海地区造成了巨大的经济损失和人员伤亡。如1960年5月22日智利瓦尔迪维亚（Valdivia）发生9.5级大地震，这也是人类历史上用仪器监测到的最强地震，随之引发的巨大海啸，在当地的海啸波波高达25 m，对附近一个居住着250万人口的城镇造成约1万人死亡、500亿～700亿美元的经济损失。此外，该海啸还影响了整个太平洋海域，导致了太平洋另一端的日本数百人丧生^[1]。

为了应对海啸灾害，人们在海啸监测、预报和致灾机理等方面做了大量研究。由于人类史上灾难性的海啸主要是由地震引起的，因此人们通过铺设海底地震仪、投放海啸浮标和设置沿岸验潮站等建立海啸监测网。地震仪和观测海底变形的大地传感器等通过海底线缆将数据传输给陆上监测中心。海啸浮标系统主要是美国国家海洋与大气管理局（NOAA）下属太平洋海洋环境实验室研制的海啸深海评估和报告浮标系统（Deep-ocean Assessment and Reporting of Tsunamis, DART），它是由海底压力记录仪、海表浮标和地球同步气象卫星构成的立体监测系统^[2]，目前已更新至第4代。目前全球海域设置了60多个DART和900余个验潮站^[3]，覆盖全面的大洋与近岸监测网为海啸预警以及海啸研究提供了大量可靠的原始数据。

随着计算机和数值计算技术的发展，数值模拟已成为海啸的主要研究方法并在海啸预警系统中起到了关键作用。目前常见的海啸模型可主要分为基于浅水方程的静压模型和考虑波浪动水压力的数值模型两类。第1类是基于浅水方程的数值模型，由于其计算效率高且在大多数模拟过程均能得到较好的结果，应用最为广泛。华盛顿大学LeVeque教授团队在1994年发布了用于求解双曲型方程的软件包Clawpack（Conservation Laws Package），且其版本一直在不断完善与更新^[4]。2010年GeoClaw^[5]作为Clawpack求解浅水方程的模块用于海啸模拟。GeoClaw包含笛卡尔坐标和球坐标两种形式，基于高分辨率激波捕捉的有限体积法进行数值求解，同时采用自适应网格嵌套技术并考虑了干湿边界的模拟。第2类常见的海啸模型主要是基于Boussinesq方程，考虑了波浪的频散性。美国特拉华大学建立了基于完全非线性Boussinesq型方程的波浪模型—FUNWAVE^[6]。该模型此后进一步更新为TVD^[7]及GPU^[8]版本。Yamazaki等^[9]建立了直接求解Navier-Stokes方程的非静压波浪模型模拟海啸的传播演化过程。Zhao等^[10]基于Green-Naghdi理论建立了考虑海底动态可变的非静压波浪模型模拟二维和三维的海啸传播过程。

虽然第2类海啸模型有更高的数值精度，但考虑到海啸预报的时效性，目前海啸模型中仍以浅水方程数值模型居多。早期的部分海啸模型如MOST^[11]等采用有限差分法（FDM）求解控制方程，其优点是简单成熟，方便构造高精度的计算格式，然而处理复杂网格时不够灵活^[12]。由于有限体积法（FVM）数学表达直观、物理涵义清晰、满足守恒规律且能方便地捕捉浅水方程中的间断解^[12]，已逐渐成为一种重要的数值模型方法。浅水方程的有限体积数值求解主要涉及源项与通量的处理以及动边界（干湿边界）的准确模拟。关于源项与通量的处理目前主要是通过不同的数值方法用以平衡有限体积法中的通量梯度和源项^[13-15]。而动边界的准确模拟目前主要基于Riemann近似解。HLL（Harten, Lax and van Leer）Riemann近似求解器^[16]在1983年提出，它根据有限体积单元界面处左右波速的预测值将界面通量值的计算划分为3个区域，能够很好地解决一维浅水流动问题。在此基础上，HLLC Riemann近似求解器^[17]考虑了二维问题中另一个方向的横波影响，将单元界面通量值的计算划分为4个区域，用以简便和精确地求解二维浅水波动问题。为了合理简洁地模拟动边界，Liang和Borthwick^[18]提出了全和谐型浅水方程，以自动满足控制方程通量与源项的守恒，并通过对地形重构抑制单元内的虚拟流动，开发了一种高性能强稳定性的浅水数值模型。该模型数值方法操作简便且自动满足单元体积各参量守恒，已被包括FUNWAVE等数值模型广泛采用^[7]。

由于灾害性海啸的影响通常波及整个大洋甚至全球范围，本文基于球坐标系下的全和谐型浅水方程建立了有限体积海啸数值模型。考虑到海啸不仅会对桥墩、防波堤等近海建筑物产生影响，也会造成近岸大面积的陆地淹水，因此海啸近岸爬高的精确模拟也是海啸数值模拟的重要技术要求。沿海地区往往是经济较为发达的区域，不仅存在着复杂的自然地形还涉及大量的人工建筑物。本文采用地形重构的方法，以模拟海啸在此复杂地形上的近岸爬高过程。本文所建立的海啸数值模型应用于2015年9月16日智利M_w8.3级地震海啸的模拟，以验证模型对越洋海啸模拟预报的能力。

2 控制方程

收起

越洋海啸的传播需要考虑地球的曲率和地转偏向力的影响，因此控制方程选取考虑科氏力的球坐标系下的非线性浅水方程，其矩阵形式为^[19]

(1) $\frac{{\partial {\boldsymbol{U}}}}{{\partial t}} + \frac{{\partial {\boldsymbol{F}}}}{{\partial \theta }} + \frac{{\partial {\boldsymbol{G}}}}{{\partial \varphi }} = {\boldsymbol{S}} ,$

式中，θ和φ分别为地球的经度和纬度；t为时间；向量U为变量；F和G分别是θ和φ方向的通量；S为源项；具体为

(2) ${\boldsymbol{U}} = \left[ {\begin{array}{*{20}{c}}\xi &{h{u_\theta }}&{h{u_\varphi }} \end{array}} \right]{\rm{ ,}}$

(3) ${\boldsymbol{F}} = \frac{1}{{r\cos \varphi }}\left[ {\begin{array}{*{20}{c}}{h{u_\theta }}&{hu_\theta ^2}&{h{u_\theta }{u_\varphi }} \end{array}} \right]{\rm{ ,}}$

(4) ${\boldsymbol{G}} = \frac{1}{r}\left[ {\begin{array}{*{20}{c}}{h{u_\varphi }}&{h{u_\varphi }{u_\theta }}&{hu_\varphi ^2} \end{array}} \right]{\rm{ ,}}$

(5) ${\boldsymbol{S}} = \left[ \begin{aligned} & \quad\quad\quad\quad\quad\quad\quad0 \\ & - \frac{{gh}}{{r\cos \varphi }}\frac{{\partial \xi }}{{\partial \theta }} + fh{u_\varphi } - {C_f}{u_\theta }\sqrt {{u_\theta }^2 + {u_\varphi }^2} \\ & - \frac{{gh}}{r}\frac{{\partial \xi }}{{\partial \varphi }} - fh{u_\theta } - {C_f}{u_\varphi }\sqrt {{u_\theta }^2 + {u_\varphi }^2} \\ \end{aligned} \right]{\rm{ ,}}$

式中，r为地球半径；ξ为相对于静止海平面的波面位移；h_s为静水深；h=ξ+h_s为总水深；u_θ和u_φ分别为θ和φ方向沿水深平均的速度。源项中第2、第3项分别对应科氏力和摩阻项，其中f=2ωsinφ为科氏力系数；ω为地球绕地轴自转的角速度；C_f = gn²/h^1/3为底摩阻系数，n是Manning系数。

为了平衡通量和源项，本文利用Liang和Borthwick^[18]提出的全和谐型浅水方程表达形式改写源项式（5）中波面的空间导数。以θ方向为例，有

(6) $h\frac{{\partial \xi }}{{\partial \theta }} = \frac{1}{2}\frac{{\partial ({\eta ^2} - 2\eta {z_{\rm{b}}})}}{{\partial \theta }} + \eta \frac{{\partial {z_{\rm{b}}}}}{{\partial \theta }}{\rm{ ,}}$

式中，z_b为基准线上的海底地形高程；η = z_b + h为波面高程。

在此不考虑z_b和h_s随时间t改变，从而有

(7) $\frac{{\partial \eta }}{{\partial t}} = \frac{{\partial h}}{{\partial t}} = \frac{{\partial \xi }}{{\partial t}}{\rm{ ,}}$

因此式（2）至式（5）的最终形式为

(8) ${\boldsymbol{U}} = \left[ {\eta \;\;\;h{u_\theta }\;\;\;h{u_\varphi }} \right]{\rm{ ,}}$

(9) ${\boldsymbol{F}} = \frac{1}{{r\cos \varphi }}\left[ {h{u_\theta }\;\;\;hu_\theta ^2 + \frac{1}{2}g({\eta ^2} - 2\eta {z_{\rm{b}}})\;\;\;h{u_\theta }{u_\varphi }} \right]{\rm{ ,}}$

(10) ${\boldsymbol{G}} = \frac{1}{r}\left[ {h{u_\varphi }\;\;\;h{u_\varphi }{u_\theta }\;\;\;hu_\varphi ^2 + \frac{1}{2}g({\eta ^2} - 2\eta {z_{\rm{b}}})} \right]{\rm{ ,}}$

(11) ${\boldsymbol{S}} = \left[ \begin{aligned} & \quad\quad\quad\quad\quad\quad\quad0 \\ & - \frac{{g\eta }}{{r\cos \varphi }}\frac{{\partial {z_{\rm{b}}}}}{{\partial \theta }} + fh{u_\varphi } - {C_f}{u_\theta }\sqrt {{u_\theta }^2 + {u_\varphi }^2} \\ & - \frac{{g\eta }}{r}\frac{{\partial {z_{\rm{b}}}}}{{\partial \varphi }} - fh{u_\theta } - {C_f}{u_\varphi }\sqrt {{u_\theta }^2 + {u_\varphi }^2} \\ \end{aligned} \right]{\rm{ ,}}$

式（9）至式（10）将地形梯度分离至方程右侧的源项中用以平衡方程左侧通量的变化，从数学表达上保证了方程的守恒形式。

此外，为方便模型用于小范围特别是物理实验的模拟，本文在建立球坐标模型的同时也开发了笛卡尔坐标海啸模型。由于两套方程的数值求解方法完全一致，后文将统一介绍。

3 数值方法

收起

本文基于均匀矩形空间单元，利用Godunov格式的有限体积法^[18]求解方程矩阵。基于MUSCL-Hancock格式，采用HLLC Riemann近似求解器^[17]计算单元界面上的流体通量，保证数值格式在时间和空间上均为二阶精度。

3.1 MUSCL-Hancock格式

3.1.1 预测步

首先在时间上求解半个时间步（Δt/2）对应的流体变量预测结果

(12) ${\boldsymbol{U}}_{i,j}^{n + 1/2} = {\boldsymbol{U}}_{i,j}^n - \frac{{\Delta t}}{{2\Delta \theta }}({{\boldsymbol{f}}_{\rm{e}}} - {{\boldsymbol{f}}_{\rm{w}}}) - \frac{{\Delta t}}{{2\Delta \varphi }}({{\boldsymbol{g}}_{\rm{n}}} - {{\boldsymbol{g}}_{\rm{s}}}) + \frac{{\Delta t}}{2}{S_{i,j}}{\rm{ ,}}$

式中，上标n代表第n时刻；i、j为有限体积单元标号；Δθ、Δφ表示球坐标系水平方向的单元大小；f_e、f_w，g_n、g_s分别代表单元东、西、北、南4个方向的流体界面通量。单元界面的流体变量利用MUSCL线性重构计算^[20]，以单元东西界面为例

(13) $ {\boldsymbol{U}}_{\rm{w}}^n = {\boldsymbol{U}}_{i,j}^n - \frac{1}{2}\varPsi (r)({\boldsymbol{U}}_{i,j}^n - {\boldsymbol{U}}_{i - 1,j}^n){\rm{ ,}} $

(14) $ {\boldsymbol{U}}_{\rm{e}}^n = {\boldsymbol{U}}_{i,j}^n + \frac{1}{2}\varPsi (r)({\boldsymbol{U}}_{i + 1,j}^n - {\boldsymbol{U}}_{i,j}^n){\rm{ }}{\rm{,}} $

式中，U_eⁿ、U_wⁿ为界面处的流体变量；Ψ(r)为用以抑制坡度较陡时可能产生数值振荡的坡度限制器^[12]，其中r为守恒变量差的比值

(15) $r = \frac{{{{\boldsymbol{U}}_{i + 1,j}} - {{\boldsymbol{U}}_{i,j}}}}{{{{\boldsymbol{U}}_{i,j}} - {{\boldsymbol{U}}_{i - 1,j}}}}{\rm{.}}$

此外，源项中海底地形底坡项直接利用中心差分法离散，科氏力项直接代入计算即可，而摩阻项后文详细介绍。

3.1.2 校正步

在校正步将变量递进至完整的时间Δt，具体为

(16) ${\boldsymbol{U}}_{i,j}^{n + 1} = {\boldsymbol{U}}_{i,j}^n - \frac{{\Delta t}}{{\Delta \theta }}({{\boldsymbol{F}}_{\rm{E}}} - {{\boldsymbol{F}}_{\rm{W}}}) - \frac{{\Delta t}}{{\Delta \varphi }}({{\boldsymbol{G}}_{\rm{N}}} - {{\boldsymbol{G}}_{\rm{S}}}) + \Delta t{S_{i,j}}{\rm{ ,}}$

式中，F_E、F_W、G_N、G_S分别为求解Riemann问题得到的单元四周的流体界面通量，单元界面处Riemann变量的求解依然利用MUSCL线性重构^[20]。

3.2 干湿边界

为了准确模拟波浪在陆地和海洋边界处的运动，本文采用Liang^[21]提出的地形重构方法。以单元东界面为例，为了确保计算过程中单元的水深恒为正，则地形重构为

(17) ${z_{b i + {1/ 2},j}} = \max (z_{b{\rm{ }}i + {1/ 2},j}^{\rm{L}},{\rm{ }}z_{{{\rm{b }}i + {1/ 2},j}}^{\rm{R}} ){\rm{ ,}}$

式中，上标L和R分别表示界面的左右侧。对应的水深为

(18) $ h_{i + {1/ 2},j}^{\rm{L}} = \max (0,{\rm{ }}\eta _{i + {1/ 2},j}^{\rm{L}} - {z_{{\rm{b }}i + {1/ 2},j}}){\rm{ ,}} $

(19) $h_{i + {1/ 2},j}^{\rm{R}} = \max (0,{\rm{ }}\eta _{i + {1/ 2},j}^{\rm{R}} - {z_{{\rm{ b}}i + {1/ 2},j}}){\rm{ }}{\rm{,}} $

综合式(17)至式(19)可求得Riemann变量η和q。

为了使干湿边界的处理不失一般性，考虑如图1干湿单元相邻的情形。单元(i, j)的东界面通量的计算仍以地形$z_{{\rm{b}}{\rm{ }}i + 1/2,j}^{\rm{R}}$为基线，而西界面的计算基线为水面线，假设水面静止，东侧地形高于西侧则水面不会发生流动，而数值计算时单元(i, j)内会引入从东至西的通量。为了避免虚假通量的引入，进一步修正z_b与η为

(20) $ {z_{{\rm{b}}{\rm{ }}i + {1/ 2},j}} = {z_{{\rm{b}}{\rm{ }}i + {1/ {2,j}}}} - \Delta {z_{i + {1/ 2},j}}{\rm{ ,}} $

(21) $\eta _{i + {1/ 2},j}^{\rm{L}} = \eta _{i + {1/ 2},j}^{\rm{L}} - \Delta {z_{i + {1/ {2,j}}}}{\rm{ ,}} $

(22) $ \eta _{i + {1/ 2},j}^{\rm{R}} = \eta _{i + {1/ 2},j}^{\rm{R}} - \Delta {z_{i + {1/ 2},j}}{\rm{ }}{\rm{,}} $

其中

(23) $\Delta {z_{i + {1/ {2,j}}}} = \max \left( {0,{\rm{ }} - {\eta ^L}_{i + {1/ 2},j} + {z_{{\rm{b}}{\rm{ }}i + {1/ 2},j}}} \right){\rm{ .}}$

经过上述修正，东西界面的计算基线得到统一，保证了通量计算的平衡性。

3.3 底摩阻

在利用MUSCL-Hancock格式数值求解后进一步利用分裂格式(splitting point-implicit scheme)对摩阻项离散^[22]，其等价于求解以下常微分方程

(24) $\frac{{{\rm{d}}{\boldsymbol{q}}}}{{{\rm{d}}t}} = {{\boldsymbol{S}}_f}{\rm{ ,}}$

式中，S_f为摩阻项，单宽流量q=[hu_θ hu_φ]^T，其随时间演进的计算为

(25) $ {{\boldsymbol{q}}^{n + 1}} = {{\boldsymbol{q}}^n} + \Delta t{\left(\frac{{{{\boldsymbol{S}}_f}}}{D}\right)^n}{\rm{ ,}} $

(26) $ D = 1 - \Delta t{\left(\frac{{\partial {{\boldsymbol{S}}_f}}}{{\partial {\boldsymbol{q}}}}\right)^n}{\rm{ }}{\rm{,}} $

式中，D为隐式系数，设F=(S_f/D)ⁿ。为了确保干湿边界附近数值计算的稳定性，对F设置阈值

(27) ${\boldsymbol{F}}\left\{ \begin{aligned} & \geqslant - {{\boldsymbol{q}}^n}/\Delta t{\rm{ }},\;\;{\rm{ }}{{\boldsymbol{q}}^n} \geqslant {\bf{0}}{\rm{ ,}} \\& \leqslant - {{\boldsymbol{q}}^n}/\Delta t{\rm{ , }}\;\;{{\boldsymbol{q}}^n} \leqslant {\bf{0}}{\rm{ }}{\rm{.}} \end{aligned} \right.$

其物理意义是摩阻作用至多会使流体停止运动而不产生逆向流动。

3.4 边界条件和稳定性准则

如图2以某个方向第一个单元为例，边界处虚拟单元的变量利用相邻单元的值近似处理，单宽流量q的取值根据边界类型而变化。对于开放边界，有

(28) ${\eta _0} = {\eta _1}, \;\; {z_{{\rm{b}}0}} = {z_{{\rm{b}}1}},\;\;{\rm{ }}{q_0} = {q_1}{\rm{ }},$

对于闭合边界、即全反射边界，有

(29) ${\eta _0} = {\eta _1},\;\;{\rm{ }}{z_{{\rm{b}}0}} = {z_{{\rm{b}}1}},\;\;{\rm{ }}{q_0} = - {q_1}{\rm{ }}{\rm{.}}$

为了保证数值计算的稳定，时间步Δt的确定遵循Courant-Friedrichs-Lewy (CFL)准则

(30) $\Delta t = C\min \left(\mathop {\min }\limits_{i,j} \frac{{\Delta {S_\theta }}}{{\left| {{u_\theta }} \right| + \sqrt {gh} }},\mathop {\min }\limits_{i,j} \frac{{\Delta {S_\varphi }}}{{\left| {{u_\varphi }} \right| + \sqrt {gh} }}\right){\rm{ }}{\rm{,}}$

式中，C为Courant系数；ΔS_θ和ΔS_φ分别为经度与纬度方向的单元空间长度。

3.5 海啸源模型

由于历史上灾害性海啸几乎都是由地震所引起的，所以海啸源模型主要考虑地震引起的海水波动过程，它是海啸数值计算的基础，其适用性直接影响海啸波的传播及海啸与海岸间相互作用。自Steketee^[23-24]把位错理论引入地震学后，许多地震学者针对不同的断层类型发展了不同的位错理论，使得计算地球同震形变的方法得到不断发展。Kajiura^[25]首次提出可以将静态的海底位移转化成海啸产生阶段自由表面的初始边界条件，使得地震海啸数值模拟过程进一步简化。为此，利用断层模型计算地震引发的海床位移量，从而估算出海底形变引起初始水面高度来初始化海啸数值模型已成为海啸模拟的惯用方法。目前，Mansinha和Smylie^[26]以及Okada^[27]基于弹性错移理论发展的两套断层模型被广泛应用，大量的研究和应用实例表明此类模型对逆冲断层地震海啸源计算具有较好的适用性^[28-30]。Okada模型^[27]计算海底同震位移需要输入7个震源参数，分别为震中位置、震源深度、断层走向角、倾角、滑动角、滑移量以及断层破裂的尺度（断层长度及宽度）。震源基本信息（震中位置、震源深度）可由地震速报系统测定，而震源机制角可由W-phase反演系统得到。震源的尺度和滑移量可参考震级与破裂尺度经验关系率估计^[31]。需要说明的是，上述方法在海底较为平缓且为典型的逆冲断层条件下较为适用。当地震引起陡坡大范围的水平运动或者地震引起断层水平运动比垂向运动大得多的走滑断层时应该考虑水平位移在海啸产生过程中的贡献。

4 2015年智利地震海啸数值模拟

收起

2015年9月16日22:54:32（UTC）位于纳斯卡板块和南美洲板块之间的南美洲俯冲带发生M_w8.3级地震。震中位于31.57°S, 71.67°W，距离智利伊亚佩尔市以西48 km，震源深度为22.4 km。该地震引发了太平洋范围的中等强度海啸，造成近场地区30多人伤亡，经济损失高达6亿美元^[32]。

4.1 海啸产生阶段模拟

考虑海啸预警业务化需求，通常假定海底地震具有一致震源机制特征，采用点源、单断层面和平均滑动量的假定以简化断层破裂的复杂性，忽略了断层破裂的局地特征以便最短时间内给出相对合理的震源破裂特征。自Kanamori和Rivera^[33]提出基于W震相测定地震矩张量的方法后，全球多个预警中心基于W震相研制出了中强震矩心矩张量反演系统。基本实现在震后8～15 min准确获取全球M_w≥6.5地震的震源机制解，从而有效提高了海啸定量预警及时性。表1为2015年智利M_w8.3级地震断层基本参数，其中地震断层破裂尺度和滑移量依据Wu等^[31]关系率估算所得，即地震断层长度约为212 km，宽度约为79 km，平均滑移量6.3 m。断层尺度量级与记录的深海海啸波长对比相当，估算的结果可信。选择USGS快速W-phase解作为该事件的震源机制解，从震源机制解可知该地震断层为典型的低倾角逆冲断层，因此可利用Okada弹性位错理论模型^[27]计算海啸传播模型所需的海表面初始形变（图3）。从海表面初始形变分布可知主要的形变呈轻微NNE向的扁状椭圆形分布，纬度方向跨度约2°，经度方向跨度约0.5°。

4.2 智利海域区域海啸模拟

智利海域区域海啸模拟范围为60°S～5°N、65°～110°W，地形数据来自ETOPO1 (https://www.ngdc.noaa.gov/mgg/global/global.html)，其中以平均海平面为基准，海域水深为负、陆地为正。模型采用分辨率为1′的均匀网格，时间步长Δt=2 s，整个模拟过程为8 h（图4）。Manning系数在水深小于100 m时取0.025，其余则取0。计算域四周设置为开放边界。本文模拟采用CPU为Intel(R) Core(TM) i7-7500U的笔记本电脑，整个模拟耗时50.5 h。

为了验证结果，本文选取了14个智利近岸测站的实测数据（各测站位置见图4）。模拟结果与测站实测比较见图5，各测站第一个测到的海啸先导波时间与波幅比较见表2。模拟与实测的海啸先导波到达时间及波幅均吻合良好。智利海域的海啸持续时间较长，且实测相较于模拟结果呈现出更多的微小波动现象，这可能是由于近岸剧烈变化的地形及岸线影响。此外，COQUIMBO测站所测最大波高达4.75 m，最大波不是先导波，且先导波后的几个波幅基本超过2.5 m，大于模拟结果。其可能的原因一方面是该测站位于震源附近，受复杂的地震影响，概化的海啸源模型没有充分反演海啸产生阶段一系列复杂的过程。另一方面是该测站所处位置地形变化剧烈，模型采用的网格未能精细刻画地形。

4.3 泛太平洋海啸模拟

泛太平洋海啸模拟区域为75°S～70°N，100°E～60°W，地形数据同样来自ETOPO1，采用空间分辨率为4′的均匀网格。Manning系数与4.2节区域海啸相同，计算域四周均为开放边界。模拟时间步长Δt=5 s，整个模拟过程为25 h，同样采用CPU为Intel（R）Core（TM）i7-7500U的笔记本电脑，用时28.0 h。

地震发生后不同时刻的瞬时波面过程见图6，第一个出现的海啸先导波到达时间见图7。海啸传播初期基本以震中为圆心向四周扩散（图6a）。当海啸经过海脊和海底山等地形时，由于其在水深较浅的地区传播速度变慢，导致海啸波峰线呈现V形，即在海脊和海底山附近的海啸波传播较慢而在深水处传播较快。由于太平洋NE侧水深较SW一侧深，所以海啸波峰线逐渐从NE−SW方向变化至SN向（图6d），其对应的先导波到达时间也呈现出NE侧快于SW侧的特点（图7）。此外，太平洋SW侧存在众多海脊、海底山及海岛，导致此处的海啸波较NE侧呈更加复杂的波面过程。

进一步选取20个分布于环太平洋沿岸及位于太平洋中部的DART浮标（见图4）验证本文所建立的数值模型，其中波面过程对比见图8、海啸先导波的到达时间和波幅对比见表3。由于从智利至太平洋东北侧的海域地形相对规则，海啸最大波通常为其先导波，且波形规则、持续时间较短。该海域的海啸模拟结果与实测吻合较好。相比于东北侧，太平洋西南侧的地形复杂，存在着海脊、海底山和海岛等地形。海啸在此受海底地形折射、绕射等影响，自由水面呈较为复杂的波动过程，海啸最大波不再是先导波。该海域本文模拟的前几个波动过程与实测基本吻合，但后续的波动过程无论在振幅还是相位上均存在一定差别。这可能是由于模型采用的网格未能精细刻画出该区域复杂的地形所致。

本文模拟的海啸先导波到达时间与实测结果随着传播距离增加其误差也逐渐增加，这可能是由于海啸波的频散效应所致。Glimsdal等^[34]提出一个与传播距离和时间呈正相关的“频散时间”以量化海啸的频散效应。这与本文模拟结果所显示的随着海啸传播距离和时间的增加先导波到达时间的误差逐渐增加是一致的。此外，实测资料显示在海啸较大波动到达之前会有一段幅度较小的减水过程，Watada等^[35]认为其与海底的弹性形变有关。由于本文的模型没有考虑该因素的影响，未能真实反演这一物理现象。

图9为智利海啸在太平洋的最大波幅分布。在智利附近海域，海啸最大波幅主要受断层走向影响，呈现出由智利沿岸向西扩散的趋势。在距离震源较远的区域，海啸最大波幅呈现由智利向日本延伸的趋势，且基本与海脊走向一致，这可能是海啸在海脊上形成所谓的海脊俘获波所致^[36]。

5 结论

收起

本文基于球坐标系下的全和谐型浅水方程，采用MUSCL-Hancock格式，并利用HLLC Riemann近似求解器计算单元界面上的流体通量，最终建立了二阶精度的Godunov格式有限体积海啸数值模型。由于近岸局部地区海啸淹没爬坡的精细化模拟和受灾评估中必将涉及到干湿边界的处理；同时，良好的干湿边界处理还有益于维持越洋海啸模拟的稳定性，因此，本文建立的模型包含了处理干湿边界的功能。在此基础上，本文利用弹性半空间位错理论得到的初始自由水面模拟了2015年9月16日智利M_w8.3级地震引发的区域海洋和泛太平洋海啸传播过程。分别通过与智利近岸14个测站和环太平洋20个DART浮标实测数据比较，进一步验证了模型模拟实际海啸过程的能力。

本文泛太平洋海啸的模拟采用空间分辨率为4′的均匀网格，在开阔海域和水深变化较缓的区域模拟结果与实测吻合良好，然而在一些岛屿及近岸受复杂地形及岸线影响的区域有一定误差，主要是由于模拟所采用的网格空间分辨率未能精细刻画这些因素的影响。为了兼顾实际模拟中的精度与效率，今后将在控制方程中考虑海水的压缩性、地震导致的海底弹性形变等因素的影响并结合多层嵌套与并行技术进一步提高模型的模拟能力。

基金

收起

国家重点研发计划（2017YFC1404205，2018YFC1407000）；国家自然科学基金面上项目（51579090）；中央高校基本科研业务费（2019B12214）

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Barrientos S E, Ward S N. The 1960 Chile earthquake: Inversion for slip distribution from surface deformation[J]. Geophysical Journal International, 1990, 103(3): 589−598.

Bernard E N, Meinig C. History and future of deep-ocean tsunami measurements[C]//Oceans'11 MTS/IEEE Kona. Waikoloa, HI, USA: IEEE, 2011: 1−7.

董杰, 田士政, 武文, 等. 全球海啸预警系统发展及其对我国的启示[J]. 海洋通报, 2019, 38(4): 368−378.

Dong Jie, Tian Shizheng, Wu Wen, et al. The development of global tsunami warning systems and their enlightenments to China[J]. Marine Science Bulletin, 2019, 38(4): 368−378.

Mandli K T, Ahmadia A J, Berger M, et al. Clawpack: Building an open source ecosystem for solving hyperbolic PDEs[J]. PeerJ Computer Science, 2016, 2: e68.

Berger M J, George D L, LeVeque R J, et al. The GeoClaw software for depth-averaged flows with adaptive refinement[J]. Advances in Water Resources, 2011, 34(9): 1195−1206.

Kirby J T, Wei G, Chen Q, et al. FUNWAVE 1.0: Fully nonlinear Boussinesq wave model-documentation and user's manual[R]. Network, USA: University of Delaware, 1998.

Shi Fengyan, Kirby J T, Harris J C, et al. A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundation[J]. Ocean Modelling, 2012, 43−44: 36−51.

Yuan Ye, Shi Fengyan, Kirby J T, et al. FUNWAVE-GPU: Multiple-GPU acceleration of a Boussinesq-type wave model[J]. Journal of Advances in Modeling Earth Systems, 2020, 12(5): e2019MS001957.

Yamazaki Y, Kowalik Z, Cheung K F. Depth-integrated, non-hydrostatic model for wave breaking and run-up[J]. International Journal for Numerical Methods in Fluids, 2009, 61(5): 473−497.

Zhao Binbin, Duan Wenyang, Webster W C. Tsunami simulation with Green-Naghdi theory[J]. Ocean Engineering, 2011, 38(2/3): 389−396.

Titov V V, Gonzalez F I. Implementation and testing of the Method of Splitting Tsunami (MOST) model[R]. Washington, USA: Pacific Marine Environmental Laboratory, NOAA, 1997.

余锡平. 近岸水波的数值方法[M]. 北京: 科学出版社, 2017.

Yu Xiping. Numerical Methods in Nearshore Wave[M]. Beijing: Science Press, 2017.

Toro E F. Shock-Capturing Methods for Free-Surface Shallow Flows[M]. Chichester, UK: John Wiley & Sons Inc, 2001.

Alias N A, Liang Qiuhua, Kesserwani G. A Godunov-type scheme for modelling 1D channel flow with varying width and topography[J]. Computers & Fluids, 2011, 46(1): 88−93.

Xia Xilin, Liang Qiuhua. A new efficient implicit scheme for discretising the stiff friction terms in the shallow water equations[J]. Advances in Water Resources, 2018, 117: 87−97.

Harten A, Lax P D, Van Leer B. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws[J]. Siam Review, 1983, 25(1): 35−61.

Toro E F, Spruce M, Speares W. Restoration of the contact surface in the HLL-Riemann solver[J]. Shock Waves, 1994, 4(1): 25−34.

Liang Qiuhua, Borthwick A G L. Adaptive quadtree simulation of shallow flows with wet-dry fronts over complex topography[J]. Computers & Fluids, 2009, 38(2): 221−234.

Baba T, Allgeyer S, Hossen J, et al. Accurate numerical simulation of the far-field tsunami caused by the 2011 Tohoku earthquake, including the effects of Boussinesq dispersion, seawater density stratification, elastic loading, and gravitational potential change[J]. Ocean Modelling, 2017, 111: 46−54.

Van Leer B. Towards the ultimate conservative difference scheme[J]. Journal of Computational Physics, 1997, 135(2): 229−248.

Liang Qiuhua. Flood simulation using a well-balanced shallow flow model[J]. Journal of Hydraulic Engineering, 2010, 136(9): 669−675.

Liang Qiuhua, Marche F. Numerical resolution of well-balanced shallow water equations with complex source terms[J]. Advances in Water Resources, 2009, 32(6): 873−884.

Steketee J A. On volterra’s dislocations in a semi-infinite elastic medium[J]. Canadian Journal of Physics, 1958, 36(2): 192−205.

Steketee J A. Some geophysical applications of the elasticity theory of dislocations[J]. Canadian Journal of Physics, 1958, 36(9): 1168−1198.

Kajiura K. The leading wave of a tsunami[J]. Bulletin of the Earthquake Research Institute, 1963, 41(3): 535−571.

Mansinha L, Smylie D E. The displacement field of inclined faults[J]. Bulletin of the Seismological Society of America, 1971, 61(5): 1433−1440.

Okada Y. Surface deformation due to shear and tensile faults in a half-space[J]. Bulletin of the Seismological Society of America, 1985, 75(4): 1135−1154.

于福江, 王培涛, 赵联大, 等. 2010年智利地震海啸数值模拟及其对我国沿海的影响分析[J]. 地球物理学报, 2011, 54(4): 918−925.

Yu Fujiang, Wang Peitao, Zhao Lianda, et al. Numerical simulation of 2010 Chile tsunami and its impact on Chinese coasts[J]. Chinese Journal of Geophysics, 2011, 54(4): 918−925.

王培涛, 于福江, 赵联大, 等. 2011年3月11日日本地震海啸越洋传播及对中国影响的数值分析[J]. 地球物理学报, 2012, 55(9): 3088−3096.

Wang Peitao, Yu Fujiang, Zhao Lianda, et al. Numerical analysis of tsunami propagating generated by the Japan M_w9.0 earthquake on Mar. 11 in 2011 and its impact on China coasts[J]. Chinese Journal of Geophysics, 2012, 55(9): 3088−3096.

闪迪, 王培涛, 任智源, 等. 有限断层模型在2015年9月16日智利M_w8.3级地震海啸数值模拟中的应用与评估[J]. 海洋学报, 2017, 39(11): 49−60.

Shan Di, Wang Peitao, Ren Zhiyuan, et al. Application and evaluation of the 16 September 2015 Illapel, Chile M_w8.3 earthquake finite fault rupture model from numerical simulation[J]. Haiyang Xuebao, 2017, 39(11): 49−60.

Wu Wennan, Zhao Li, Wu Y M. Empirical relationships between aftershock zone dimensions and moment magnitudes for plate boundary earthquakes in Taiwan[J]. Bulletin of the Seismological Society of America, 2013, 103(1): 424−436.

Hayes G P, Myers E K, Dewey J W, et al. Tectonic summaries of magnitude 7 and greater earthquakes from 2000 to 2015[R]. USA: USGS, 2016.

Kanamori H, Rivera L. Source inversion of W phase: Speeding up seismic tsunami warning[J]. Geophysical Journal International, 2008, 175(1): 222−238.

Glimsdal S, Pedersen G K, Harbitz C B, et al. Dispersion of tsunamis: Does it really matter?[J]. Natural Hazards and Earth System Sciences, 2013, 13(6): 1507−1526.

Watada S, Kusumoto S, Satake K. Traveltime delay and initial phase reversal of distant tsunamis coupled with the self-gravitating elastic earth[J]. Journal of Geophysical Research: Solid Earth, 2014, 119(5): 4287−4310.

王岗, 胡见, 王培涛, 等. 双曲余弦海脊上海啸俘获波的解析与数值研究[J]. 海洋学报, 2018, 40(5): 15−23.

Wang Gang, Hu Jian, Wang Peitao, et al. Analytical and numerical investigation of tsunami trapped waves over a hyperbolic-cosine squared ocean ridge[J]. Haiyang Xuebao, 2018, 40(5): 15−23.

2021年第43卷第5期

PDF下载

124

引用本文

BibTeX

文章信息

doi: 10.12284/hyxb2021095

接收时间：2020-05-24
首发时间：2026-02-26
出版时间：2021-05-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2020-05-24
修回日期：2020-10-19

基金

国家重点研发计划（2017YFC1404205，2018YFC1407000）；国家自然科学基金面上项目（51579090）；中央高校基本科研业务费（2019B12214）

作者信息

¹海岸灾害及防护教育部重点实验室（河海大学），江苏南京 210098

²国家海洋环境预报中心，北京 100081

³拉夫堡大学建筑与土木工程学院，英国莱斯特郡 LE11 3TU

通讯作者:

王岗（1982-），男，河北省张家口市人，博士，副教授，主要从事港湾共振、水波模拟及海啸与洪水风险评估研究。E-mail：gangwang@hhu.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/hyxb/CN/10.12284/hyxb2021095

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

断层参数	参数值
震源位置	31.57°S, 71.67°W
震源深度	22.4 km
断层长度	212 km
断层宽度	79 km
倾角	19°
滑动角	83°
走向角	353°
平均滑移量	6.3 m

断层参数

参数值

震源位置

31.57°S, 71.67°W

震源深度

22.4 km

断层长度

212 km

断层宽度

79 km

倾角

19°

滑动角

83°

走向角

353°

平均滑移量

6.3 m

测站名称	位置	到达时间		先导波波幅
COQUIMBO CL	30.0°S,71.3°W	28	32	12.5%		104.0	96.9	6.8%
SAN ANTONIO CL	33.6°S,71.6°W	30	38	21.1%	73.4	65.3	11.0%
BUCALEMU CL	34.6°S,72.0°W	41	49	16.3%	48.5	46.9	3.3%
HUASCO CL	28.5°S,71.2°W	30	35	14.3%	46.9	46.2	1.5%
CONSTITUCION CL	35.4°S,72.5°W	53	63	15.9%	68.6	49.7	27.6%
CHANARAL CL	26.4°S,70.6°W	57	61	6.6%	41.8	42.5	1.6%
TALCAHUANO CL	36.7°S,73.1°W	104	108	3.7%	43.4	44.0	1.4%
TALTAL CL	25.4°S,70.5°W	67	67	0%	18.8	19.5	3.6%
JUAN FERNANDEZ	33.6°S,78.8°W	63	60	4.8%	45.4	40.9	9.9%
CORRAL CL	39.9°S,73.4°W	112	132	15.2%	26.1	24.4	6.5%
TOCOPILLA CL	22.1°S,70.2°W	93	92	1.1%	12.0	10.8	10.0%
SAN FELIX CL	26.3°S,80.1°W	80	81	1.2%	23.8	22.6	5.0%
MATARANI PE	17.0°S,72.1°W	141	140	0.7%	13.3	13.9	4.3%
CALLAO LA-PUNTA PE	12.1°S,77.2°W	206	218	5.5%	23.1	20.3	12.1%

测站名称

位置

到达时间

先导波波幅

模拟/min

实测/min

相对误差

模拟/cm

实测/cm

相对误差

COQUIMBO CL

30.0°S,71.3°W

12.5%

104.0

96.9

6.8%

SAN ANTONIO CL

33.6°S,71.6°W

21.1%

73.4

65.3

11.0%

BUCALEMU CL

34.6°S,72.0°W

16.3%

48.5

46.9

3.3%

HUASCO CL

28.5°S,71.2°W

14.3%

46.9

46.2

1.5%

CONSTITUCION CL

35.4°S,72.5°W

15.9%

68.6

49.7

27.6%

CHANARAL CL

26.4°S,70.6°W

6.6%

41.8

42.5

1.6%

TALCAHUANO CL

36.7°S,73.1°W

104

108

3.7%

43.4

44.0

1.4%

TALTAL CL

25.4°S,70.5°W

18.8

19.5

3.6%

JUAN FERNANDEZ

33.6°S,78.8°W

4.8%

45.4

40.9

9.9%

CORRAL CL

39.9°S,73.4°W

112

132

15.2%

26.1

24.4

6.5%

TOCOPILLA CL

22.1°S,70.2°W

1.1%

12.0

10.8

10.0%

SAN FELIX CL

26.3°S,80.1°W

1.2%

23.8

22.6

5.0%

MATARANI PE

17.0°S,72.1°W

141

140

0.7%

13.3

13.9

4.3%

CALLAO LA-PUNTA PE

12.1°S,77.2°W

206

218

5.5%

23.1

20.3

12.1%

浮标号	位置	到达时间		先导波波幅
32402	26.7°S,74.0°W	46	39	15.2%		9.3	10.0	7.0%
32401	20.5°S,73.4°W	103	96	6.8%	5.5	5.1	7.3%
32412	18.0°S,86.4°W	174	171	1.7%	6.3	6.1	3.2%
32411	5.0°N,90.9°W	398	394	1.0%	2.5	1.9	24.0%
43413	11.0°N,100.1°W	476	473	0.6%	1.9	1.4	26.3%
43412	16.0°N,107.0°W	560	560	0%	2.5	1.9	24.0%
46411	39.3°N,127.1°W	837	839	0.2%	1.3	0.8	38.5%
46407	42.7°N,127.8°W	874	876	0.2%	1.1	1.3	15.4%
46404	45.8°N,128.8°W	910	910	0%	1.0	0.7	30.0%
46409	55.3°N,148.5°W	1 030	1 034	0.4%	1.2	1.8	33.3%
46403	52.7°N,157.0°W	1 032	1 038	0.6%	1.4	1.9	26.3%
46408	49.6°N,169.9°W	1 063	1 073	0.9%	1.6	2.2	27.3%
46413	48.0°N,174.2°W	1 074	1 082	0.7%	1.3	1.8	27.8%
51407	19.6°N,156.5°W	872	878	0.7%	1.8	1.9	5.3%
51425	9.5°S,176.3°W	879	891	1.3%	1.4	0.8	42.9%
52406	5.3°S,165.0°E	1 075	1 091	1.5%	0.3	0.4	25.0%
52401	19.3°N,155.8°E	1 220	1 238	1.5%	1.3	1.0	23.1%
52402	11.9°N,153.9°E	1 224	1 242	1.4%	1.1	0.7	36.4%
52403	4.0°N,145.5°E	1 300	1 318	1.4%	0.8	0.3	62.5%
21414	49.0°N,178.2°E	1 113	1 121	0.7%	1.1	1.4	21.4%

浮标号

位置

到达时间

先导波波幅

模拟/min

实测/min

相对误差

模拟/cm

实测/cm

相对误差

32402

26.7°S,74.0°W

15.2%

9.3

10.0

7.0%

32401

20.5°S,73.4°W

103

6.8%

5.5

5.1

7.3%

32412

18.0°S,86.4°W

174

171

1.7%

6.3

6.1

3.2%

32411

5.0°N,90.9°W

398

394

1.0%

2.5

1.9

24.0%

43413

11.0°N,100.1°W

476

473

0.6%

1.9

1.4

26.3%

43412

16.0°N,107.0°W

560

2.5

1.9

24.0%

46411

39.3°N,127.1°W

837

839

0.2%

1.3

0.8

38.5%

46407

42.7°N,127.8°W

874

876

0.2%

1.1

1.3

15.4%

46404

45.8°N,128.8°W

910

1.0

0.7

30.0%

46409

55.3°N,148.5°W

1 030

1 034

0.4%

1.2

1.8

33.3%

46403

52.7°N,157.0°W

1 032

1 038

0.6%

1.4

1.9

26.3%

46408

49.6°N,169.9°W

1 063

1 073

0.9%

1.6

2.2

27.3%

46413

48.0°N,174.2°W

1 074

1 082

0.7%

1.3

1.8

27.8%

51407

19.6°N,156.5°W

872

878

0.7%

1.8

1.9

5.3%

51425

9.5°S,176.3°W

879

891

1.3%

1.4

0.8

42.9%

52406

5.3°S,165.0°E

1 075

1 091

1.5%

0.3

0.4

25.0%

52401

19.3°N,155.8°E

1 220

1 238

1.5%

1.3

1.0

23.1%

52402

11.9°N,153.9°E

1 224

1 242

1.4%

1.1

0.7

36.4%

52403

4.0°N,145.5°E

1 300

1 318

1.4%

0.8

0.3

62.5%

21414

49.0°N,178.2°E

1 113

1 121

0.7%

1.1

1.4

21.4%