[发明专利]一种基于分区异步长解的地震波动分析方法

摘要

本发明公开了一种基于分区异步长解的地震波动分析方法，包括S1、按照网格尺寸将计算区域分为I区域和II区域；S2、II区域按照计算时步取Δt＝Δts进行求解，当计算时步达到Δt＝Δtn时，对I、II区域进行求解，确定波动传播的初始位移、初始速度和初始加速度；S3、根据初始时刻t的加速度值，计算得到t+Δt时刻波动传播的位移与速度的预测值；S4、根据所述波动传播的位移和速度预测值，采用动力平衡方程求解t+Δt时刻的加速度，并确定该时刻波动传播的位移和速度；S5、重复步骤S2‑S4，直至全部时间结束。

主权项

1.一种基于分区异步长解的地震波动分析方法，其特征在于：包括S1、按照网格尺寸将计算区域分为I区域和II区域；S2、II区域按照计算时步取Δt＝Δt_s进行求解，当计算时步达到Δt＝Δt_n时，对I、II区域进行求解，确定波动传播的初始位移、初始速度和初始加速度；其中，Δt为时间步长，Δt_s为网格尺寸较小区域的计算时步，Δt_n为网格尺度较大区域的计算步长，Δt_n＝nΔt_s；S3、根据初始时刻t的加速度值，计算得到t+Δt时刻波动传播的位移与速度的预测值；S4、根据所述波动传播的位移和速度预测值，采用动力平衡方程求解t+Δt时刻的加速度，并确定该时刻波动传播的位移和速度；S5、重复步骤S2‑S4，直至全部时间结束。