Impact of Different Source Incidence Angles on Seismic Response of Submarine Basin
-
摘要: 为研究地震作用对海底盆地的影响,在分析网格大小和PML边界厚度对计算精度影响的基础上,综合考虑了海水、盆地内不规则地形以及多层介质等因素,运用谱元法理论建立了礼乐盆地地震动分析二维模型,研究了盆地在不同震源入射角下的响应。结果表明:在考虑流固耦合的情况下,网格尺寸为四分之一波长时即可取得较高精度,PML边界厚度至少为2个波长时方可达到理想效果。震源的不同会影响放大系数的大小及频谱特征。当地震波从盆地中部入射时,随着震源入射角的增大,放大系数逐渐减小;在不同震源入射角下,盆地不同位置的加速度时程及频谱特征差异较大。当地震波从左侧入射时,放大系数的分布规律与地震波从盆地中部入射时有所不同;在不同的震源入射角下,盆地不同位置的加速度时程及频谱特征差异较小。在研究盆地地震相关问题时,海水、地形以及多层介质的影响是不可忽视的。Abstract: In order to study the influence of seismic action on the submarine basin, based on the analysis of the influence of grid size and PML boundary thickness on the calculation accuracy, a two-dimensional seismic analysis model of the Lile basin is established by using the spectral element method, taking into account the factors of seawater, irregular topography in the basin and multi-layer media, and the response of the basin at different incident angles of earthquake sources is studied. The results show that high accuracy can be achieved when the mesh size is a quarter wavelength and the PML boundary thickness is at least two wavelengths when considering the fluid-solid coupling. The magnification factor and spectrum characteristics are influenced by the source. When seismic waves are incident from the middle of the basin, the amplification coefficient decreases with the increase of the incident angle of the source. Acceleration time history and frequency spectrum characteristics differ greatly in different locations of the basin at different source incidence angles. When the seismic wave is incident from the left side, the distribution law of the amplification coefficient is different from that when the seismic wave is incident from the middle of the basin. Acceleration time history and spectrum characteristics differ slightly in different locations of the basin at different source incidence angles. The influence of seawater, topography and multi-layered media can not be ignored in the study of basin earthquake-related problems.
-
Key words:
- Basin /
- Spectral element method /
- Ground motion
-
引言
我国是地震灾害多发的国家 (高孟潭等,2008)。如何快速有效地获取和评估灾情信息已成为降低地震损失的重点研究内容 (崔丽萍等,2016)。目前,我国在地震后主要依赖于传统的搜救犬或生命探测仪等进行搜救工作。由于生命探测仪在工作时受距离、光线、废墟材质、温度和湿度等因素影响,在进行搜救时具有盲目性和不确定性,这在一定程度上限制了其应用。除了直接探测人体生命特征的探测技术,通过探测受灾人员身边的一些能发出信号的物品,以此来找到受灾人员的相关技术也逐渐成为生命探测领域的重要发展方向。
随着无线通信技术的发展,约80%的人持有带Wi-Fi的智能手机和平板电脑 (Chen等,2010a),绝大部分人都将手机随身携带。当地震引发建筑物坍塌时,手机和平板电脑成为散落在人们身边的仅有设备。通过对这些设备进行探测和定位,可以间接地找到被压埋的人员的大概位置,再配合生命探测设备,即可将人员及时救出。
目前手机检测方法分为被动式和主动式2种 (柴伦尚,2014)。由于被动式手机检测通过手机信号进行探测,要求手机处于联网状态,因此,该方法对处于通信网络中断状态的手机是无效的。地震发生后,房屋道路坍塌,各处线路中断,基站被破坏,加上暴雨等次生灾害的影响,往往导致通信链路及无线网络中断 (李世雄等,2016),手机等设备无法与基站进行通信,且被压埋者通常丧失行为能力,无法对手机进行操作或尝试连接网络,因此被动式手机检测方法不适合用于检测救援。主动式手机检测方法是指手机主动对外发射某种信号,外界通过对该信号进行采集和处理,再根据定位算法得到手机位置,代表技术主要有全球定位系统 (Global Positioning System,GPS)、辅助全球定位系统 (Assisted Global Positioning System,AGPS)、基站通讯系统 (Location Based Service,LBS) 及Wi-Fi定位技术 (Wi-Fi Positioning System)。GPS一般用于空旷的户外定位,进行复杂的室内定位时,存在较大误差;AGPS需要基站辅助GPS精确定位,而震后基站被毁,无法进行定位;LBS同样需要手机处于联网状态时才能用来对手机定位。Wi-Fi定位通过对手机发射的Wi-Fi信号进行采集与滤波处理,结合定位算法,得到手机位置。因构建Wi-Fi无线网络和使用矢量分析仪成本较低,且操作简易,本文采用基于主动手机探测Wi-Fi定位方法,对废墟信道进行测量分析。
研究基于Wi-Fi的手机定位救援方法,首先需要分析地震废墟环境下Wi-Fi无线信号的传播特征及相应的传播规律,为基于信号强度指示 (Received Signal Strength Indication,RSSI) 定位算法的运用提供理论支撑。
传统的无线信道研究主要包括对城市区域的室内、室外、室内到室外3种范围的传播模型研究 (Turin等,1972;Toeltsch等,2002),其结论表明,传统的无线信道传播模型并不适用于震后废墟环境。只有通过对震后废墟环境下的无线信道进行理论估计和实际测量,才能获得其传播特性,从而建立适用的无线信道模型。国外对灾害环境下的无线信道传播模型有一定研究,但都是针对移动通信标准的频段,如434MHz、868MHz (DiCarlofelice等,2015)、900MHz和1800MHz (Chai等,2013),而国内相关的研究较少。柴伦尚 (2014)提出“基于GSM和DCS体制的手机定位方法的无线信道研究”,也是针对移动通信频段电波信号传播的无线信道进行建模仿真,实现主动式手机探测定位方法。
本文对废墟下Wi-Fi 2.4G频段传输的无线信号进行测量,对数据处理后得到信道冲击响应,并得到功率的延迟分布,以分析信道大小尺度的表征参数。在实际废墟测量场景下,分析了信号的频率响应,最终建立了路径损耗模型。该模型分为自由空间损耗和穿透废墟的损耗两部分,然后对该模型进行了仿真验证。小尺度参数主要包括平均超量时延和均方根时延。
1. 废墟环境下无线电波传播特性
地震灾害发生后的废墟坍塌状况在受灾程度和建筑环境等因素的影响下各不相同,不同建筑材料对信号的损耗也各不相同,导致路径损耗衰减指数变化较大 (柴伦尚,2014;Jiao等,2014),因此无法用通用模型模拟不同的灾害环境,分析和提取有用的随机信道参数,并得到灾害现场信道的真实传播规律。
无线信号在传输过程中遇到障碍物时有可能发生直射、反射、绕射和散射现象 (郑金金,2013;Cheng等,2015)。基于这些现象的影响,信号在无线信道中的衰减一般分为大尺度衰减和小尺度衰落2种 (刘江涛,2011;Jayaram等,2014)。
1.1 大尺度衰减
大尺度衰减特性由无线电波传输路径损耗和阴影衰减组成。无线通信中,发射机和接收机的距离越远,接收到的信号强度越小,即信号强度的损耗随着距离的增大而增大,这种强度的变化是指数型的。根据路径损耗模型,如果已知发送端无线电波的强度和距离,则可以推算出接收端无线信号强度;如果已知收发两端的无线信号强度,则可以推算出发送机和接收机之间的距离,这种理论广泛用于定位算法中,例如基于RSSI的各种定位算法等。阴影衰减效应是指无线电波在传播过程中碰到较大的障碍物时 (如剧烈的地势变化以及高大建筑物阻挡),穿透障碍物导致其能量发生衰减。这两种衰减的产生都是由于信号传输距离较远,故为大尺度特性。
自由空间的路径损耗可表示为:
$$ {P_{\rm{r}}}(d) = \frac{{{P_{\rm{t}}}{G_{\rm{t}}}{G_{\rm{r}}}{\lambda ^2}}}{{{{(4{\rm{ \mathsf{ π} }})}^2}{d^2}L}} $$ (1) 式中,Pt表示发射端发射的电磁波的功率,Gt发射天线的增益,Gr接收天线的增益,λ表示电磁波的波长,d为接收天线和发射天线之间的距离,L表示系统损耗因子 (与传播无关)。接收和发送天线信号功率差值为损耗,记为PL-free ,换算成分贝为:
$$ {P_{L{\rm{ - free}}}} = - 10\lg ({P_{\rm{r}}}/{P_{\rm{t}}}) = - 10{\rm{lg}}\left[ {\frac{{{G_{\rm{t}}}{G_{\rm{r}}}{\lambda ^2}}}{{{{\left({4{\rm{ \mathsf{ π} }}} \right)}^2}{d^2}}}} \right]\left({{\rm{dB}}} \right) $$ (2) 发射天线和接收天线为单位增益Gt=Gr=1时可表示为:
$$ {P_{L{\rm{ - free}}}} = - 10\lg ({P_{\rm{r}}}/{P_{\rm{t}}}) = - 10\lg \left[ {\frac{{{\lambda ^2}}}{{{{\left({4{\rm{ \mathsf{ π} }}} \right)}^2}{d^2}}}} \right] $$ (3) 或表示为:
$$ {P_{L{\rm{ - free}}}} = - 10\lg \left[ {\frac{{{{(c/f)}^2}}}{{{{(4{\rm{ \mathsf{ π} }})}^2}{d^2}}}} \right] = 20\lg d + 20\lg f - 10\lg \left[ {\frac{{{c^2}}}{{{{(4{\rm{ \mathsf{ π} }})}^2}}}} \right] $$ (4) 式中c为光速,f表示频率。由式 (4) 可知路径损耗随收发天线间距离和信号频率的增大而增大。现有的研究中使用的信号最高频率为移动通信频段的1800MHz,因其频率高导致信号穿墙损耗严重,通信频段定位需要与基站通信,不适合地震等灾后环境。基于Wi-Fi的定位救援原理是在现场布排接入点 (Access Point,AP,即简单路由),激发手机与其通信,捕获信号并进行定位分析,相比基站定位,有较高时效性和可操作性。
1.2 小尺度衰减
相比于大尺度衰减,无线信道小尺度衰减特性指当待测设备移动较小的距离时,接收信号在短期内发生的非常小的变化,表现为接收信号在幅度、相位或时延上快速变化的现象。小尺度衰减主要是由于多径现象和多普勒现象造成的。无线电波在实际信道传输过程中,信号遇到障碍物时会发生散射、反射和绕射现象,产生多条路径,最终到达接收端的信号相当于多条路径上的信号叠加所得结果,即多径现象。多普勒现象是由于接收机和发射机之间相对运动造成的。
小尺度衰减有2个主要特点,分别为电磁波的时延扩展和无线信道的时变特点。多径现象使信号在时域表现为时延扩展,可用平均超量时延和均方根时延来表示;在频域表现为选择性衰减,可以用相关带宽来表示。在极小时间间隔内,废墟下无线信道可视为快速时不变系统 (Chai等,2013),所以平均超量时延及均方根时延都可以使用功率延迟分布 (PDP) 来计算。功率延迟分布为本地的瞬时功率分布的平均值,实际上常以信道的脉冲响应来计算,即:
$$ p(t) = E({\left| {h(t)} \right|^2}) $$ (5) 式中,p(t) 为功率延迟,h(t) 为信道脉冲响应。功率延迟分布的一阶矩为平均超量时延 (Mean Excess Delay),用τ表示:
$$ \overline \tau = \frac{{\sum {t \cdot p(t)} }}{{\sum {p(t)} }} = D $$ (6) 均方根 (Root Mean Square) 时延扩展τrms是平均超量时延的标准差:
$$ {\tau _{{\rm{rms}}}} = \sqrt {\overline {{\tau ^2}} - {{(\bar \tau )}^2}} = \sqrt {\frac{{\sum {{{\left( {t - D} \right)}^2} \cdot p(t)} }}{{\sum {p(t)} }}} $$ (7) 相关带宽用Bc表示,与选择的信道频率正相关,可以通过均方根时延τrms进行计算:
$$ {B_{\rm{c}}} = \frac{1}{{2{\rm{ \mathsf{ π} }} \cdot {\tau _{{\rm{rms}}}}}} $$ (8) 在无线信道的小尺度建模之中,莱斯分布和瑞利分布是使用最广泛的两种建模方式。
2. 废墟环境信道测量
当前主要的信道测量方法有2种 (孟娟,2009),即时域测量法和频域扫频测量法。前者通过发射端不断发送脉冲信号,用带通滤波器在接收端进行接收,再对数据进行一系列数学处理即可得到信道脉冲响应。后者在频域进行测量,在接收端捕捉穿过信道后的数据,再利用傅里叶反变换即可获得信道的复数响应。时域测量系统对发送信号的时间分辨率及接收端示波器采样率的要求较高,导致系统复杂度和成本较高,加上废墟环境本身存在诸多的影响因素,所以本文主要采用频域扫频测量方法。该方法操作较容易,设备需求少,需要的主要仪器为矢量网络分析仪 (Vector Network Analyzer,VNA),只要增加所测量的带宽即可提高系统分辨率。通过对不同类型的实际废墟场景下的无线信道进行频域测量,以及对测量数据进行处理,得到各类型废墟场景的无线信道参数表征和不同类型的信道转换函数 (Channel Transmission Function,CTF),再对其进行傅里叶反变换得到时域响应函数。
2.1 信道测量设计
由于所测环境为震后废墟场景,震后供电一般遭到破坏,再加上测量环境复杂,为了更快、更高效地确定该环境下信道模型,本文使用自带储电功能且方便携带的矢量分析仪,实验采用安捷伦手持式矢量网络分析仪N9918A。
本文运用频域扫频测量方法进行信道测量,测量系统框架如图 1所示。系统包括手持式矢量网络分析仪N9918A、功率放大器、一对收发天线。实验中所有的设备通过同轴电缆连接,所测数据则通过USB接口保存到U盘中。在无线信道测量实验中,天线参数的选择对结果影响较大。由于常用路由器天线在水平方向上是全向的,实验所用的路由器也为全向天线。Chai等 (2013)的研究表明,垂直极化的接收信号强度大于水平极化方式,故本实验选择使用垂直极化方式接收信号。
2.2 废墟环境无线信道测量
本文中进行信道测量的现场是“5·12”汶川地震的极重灾区——北川老县城的震后废墟堆。接收天线与矢量网络分析仪直接相连,发射天线固定在废墟内部高1m的支架上,测量时,接收天线固定不变,发射天线可随支架自由移动 (模拟救援定位过程)。废墟内部设置一个N×M的矩形方格阵,如图 2,间隔为1m。在测量时,方格阵尽可能覆盖整个废墟内部可能有人员生还的小空间。对方格阵每个节点上的值进行多次测量,取均值作为结果进行最终计算。
实际测量情况见图 3,测量场景为北川县曲山小学遗址。图中为小学3层教学楼,震后基本全部倒塌,选择废墟内部较大的一个空间进行测量,入口为图中所示C处,所测区域在图中C处的右侧,长7m,宽5m,高2m,设置4×6点阵,记为P1—P24。测量示意图如图 4,五角星表示接收天线与网格的相对位置,实际接收天线在测量空间入口处,即图 3中的C处。
3. 废墟信道参数分析
对矢量网络分析仪所采集的数据按式 (9) 进行计算,获得频域的信道传递函数:
$$ H(f) = \sum\limits_{k = 0}^{N - 1} {H(k)\delta (f - k\Delta f - {f_0})} $$ (9) 式中,Δf是频率步长,N是步数,f0是起始频率,H(k) 是测量时记录的数据。对于线性信道,离散脉冲响应h(n) 可使用离散傅里叶反变换IDFT由频率响应计算得出:
$$ h(n) = 1/N\sum\limits_{k = 0}^{N - 1} {H(k){{\rm{e}}^{\frac{{ - 2{\rm{ \mathsf{ π} }}jkn}}{N}}}} $$ (10) 式中,N代表样本数目。时域采样间隔T等于1/Δf,步长$\Delta t = T/N = 1\left({N \cdot \Delta f} \right)$。由公式 (11),运用信道脉冲响应h(n) 计算得到功率延迟分布,得到文中所示废墟环境下无线信道模型大小尺度参量模型。
3.1 信道大尺度参数分析
功率时延分布 (Power Delay Profile,PDP) 描述信道在时间上的色散,指发射功率经过多径信道的传输后在接收端的能量分布,假设为广义平稳非相关散射 (Wide Sense Stationary Uncorrelated Scattering,WSSUS) 信道,则可用下式计算:
$$ P(t) = \left\langle {\left| {h{{(t)}^2}} \right|} \right\rangle = {\left| {h(t)} \right|^2} = {\left| {\sum\limits_{n = 0}^{N - 1} {h(n)\delta (t - \Delta t)} } \right|^2} $$ (11) 式中,〈…〉表示取括号中值的平均值;N为采样数,本文中为201;Δt为时间步长,本文为10ns。
信道大尺度衰减用PL表示,包括无线信号在空气中的传播损耗和穿透废墟损耗,两者分别用Afree和Aruin表示。PL可表示为:
$$ {P_L} = {A_{{\rm{free}}}} + {A_{{\rm{ruin}}}} $$ (12) 通过对多次测量的数据进行滑动平均计算,得到测量场景下的频率响应,计算得到h(n),并对其运用滑动平均的方式进行处理。
P1和P24分别为测量中路径损耗最大和最小的2个点,其路径损耗情况如图 5所示。主要原因是P1距离天线最近,损耗较小,而P24离天线最远,损耗比较大。通过计算,路径损耗标准差在[1.38dB,4.14dB]区间变动。
在对场景进行信道特征实际测量时,记录了各个测量点收发端的相对位置,即能确定接收和发射天线之间距离,由式 (13) 可求得无线信号在传播过程中的损耗:
$$ {P_{L{\rm{ - free}}}} = - 10\lg \left[ {\frac{{{{\left({c/f} \right)}^2}}}{{{{\left({4{\rm{ \mathsf{ π} }}} \right)}^2}{d^2}}}} \right] = - 20\lg \left[ {\frac{c}{{4{\rm{ \mathsf{ π} }}fd}}} \right] $$ (13) 式中c是光速,f是无线信号频率,d是收发天线之间距离。结合式 (14) 可以计算得到信号穿透废墟损耗:
$$ {A_{{\rm{ruin}}}} = {P_L} - {A_{{\rm{free}}}} = {P_L} - {P_{L{\rm{ - free}}}} $$ (14) 在不同频段,不同障碍物对无线信号的阻碍效果是不同的,地震后,墙体坍塌后露出的钢筋支架、水泥、钢筋混凝土墙体、木头或合成材料造的家具等混合在一起,共同阻碍信号的传播。理论上,可以将混合材料分割后抽象为若干层障碍物叠加的模型,估计损耗可由式 (15) 表示:
$$ {A_{{\rm{ruin}}}}^{\rm{*}} = \sum\limits_{k = 1}^N {{a_k}{h_k}} $$ (15) 式中,ak表示第k层材料的损耗因子,hk表示第k层障碍物的厚度。信号经过不同材料时的穿透损耗是由材料的内部组成结构决定的,针对不同种类单位数量的材料对各频段信号的阻碍效果,一些文献 (Chai等,2013;Chen等,2010b) 已进行了详细的研究与说明。如图 6(柴伦尚,2014) 所示,可以推断出在2400MHz,钢筋混凝土层的穿透损耗为0.463—0.464dB/cm。
本文中测量的废墟场景下,收发天线间主要的阻挡物是钢筋混凝土材质的墙体,测量得到穿透损耗值为约0.92dB/cm。因此,在接下来的计算中,取频率在2450MHz时测得的穿透钢筋混凝土的损耗因子0.92dB/cm。由此,可用公式${A_{{\rm{ruin}}}} = {a_{{\rm{ferroconcrete}}}} \times {h_{{\rm{ferroconcrete}}}}$技术出穿透损耗预测值,式中aferroconcrete=0.92dB/cm,hferroconcrete=28cm为本文场景中的钢筋混凝土层的厚度,再根据公式 (14) 计算得到穿透损耗测量值。在2450MHz时的无线信号穿过钢筋混凝土层的预测值和测量值如图 7。
由图可知,实际测量值围绕预测值上下波动,但不是很剧烈。原因是实际测量环境下混凝土层不可能完全一样,不同点的厚度、密度及内部材料的构排都存在差值,属于不可避免的误差,加上震后混凝土层上碎土和碎石堆积以及裂缝等,造成电磁波穿过混凝土层时被吸收、散射、衍射和泄露,使各点的衰减值各不相同。在P1—P20,混凝土损耗测量值呈上升趋势,这是因为测量点的ID值越大越靠近横梁,信号穿透的横梁厚度有细微的变化;在P21—P24点位,测量值也呈起点不同的上升趋势,这是由于混凝土层并不是水平的,而是有一定的斜度,导致各点到所穿透点的距离不完全一样,存在细小差别。经过计算测量值与预测值误差约为3.81dB。
3.2 信道小尺度参数分析
本节主要研究小尺度特性的2个参数特征,分别为平均超量时延τ和均方根时延扩展τrms。平均超量时延是衡量信号多径分量的能量集中度,可以通过对多径分量幅值的平方进行加权平均计算得到:
$$ \bar \tau = \frac{{\sum {t \cdot P(t)} }}{{\sum {P(t)} }} = \frac{{\sum\limits_{n = 1}^N {{\tau _n}{{\left| {{\alpha _n}} \right|}^2}} }}{{\sum\limits_{n = 1}^N {\alpha _n^2} }} $$ (16) 式中,P(t) 代表测量数据中的功率延迟分布,已由式 (11) 求出,αn表示在测量中不同多径分量的幅度值,τn表示不同多径分量对应的各自的时延。
均方根时延用τrms来描述无线信道中多径效应引起的频率选择性衰减现象。数字通信系统中,如果2个脉冲的发送间隔Δt < τrms会导致相邻的脉冲发生串扰,使得接收端无法正确判别信号。为了避免码间串扰,要求Δt > τrms,而增大Δt会限制码速和系统传输速率。均方根时延可由式 (17) 来计算:
$$ {\tau _{{\rm{rms}}}} = \sqrt {\overline {{\tau ^2}} - {{\left({\overline \tau } \right)}^2}} = \sqrt {\frac{{\sum\limits_{n - 1}^N {\tau _n^2{{\left| {{\alpha _n}} \right|}^2}} }}{{\sum\limits_{n - 1}^N {\alpha _n^2} }} - {{\left({\frac{{\sum\limits_{n - 1}^N {{\tau _n}{{\left| {{\alpha _n}} \right|}^2}} }}{{\sum\limits_{n - 1}^N {\alpha _n^2} }}} \right)}^2}} $$ (17) 实际通信中,由于频带有限宽,信号码元都存在一定的扩展,接收到的每一个码元都会受到在该脉冲之前其他码元的扩展干扰,距离该脉冲越近,干扰越大。
4. 结论
本文通过对“5·12”汶川地震灾害区北川老县城遗址的实际测量,分析了地震废墟环境下无线信号传播模型的参数特性。主要从大尺度和小尺度衰减两方面进行分析。路径损耗主要通过反傅里叶变换进行计算和分析,并计算得到其标准差在[1.38dB,4.14dB]间波动。另外,研究了钢筋混凝土对Wi-Fi 2.4G无线信号的穿透损耗,通过分析其耗损因子为aferroconcrete=0.92dB/cm。小尺度衰减包括平均超量时延和均方根时延扩展的特性。分析表明,相比移动通信频段,Wi-Fi 2.4G信号频率较高,无论是在自由空间传播还是穿透废墟,都有比较快的衰减,限制了传播距离。但震后基站被毁,手机和基站通讯中断的情况下,无法采集手机GSM900和GSM1800信号进行定位。本文研究基于手机Wi-Fi的主动探测方法,可通过尽可能多的安置接入点来克服传播距离短的缺点,提高定位精度,这比其它被动手机检测更容易实现,同时具备便利和低成本的优势。本文的结果可为下一步进行震后手机定位探测研究提供理论基础。
-
表 1 模型介质参数
Table 1. Model media parameters
介质 $ \rho $/(kg·m−3) Vp/(m·s−1) Vs/(m·s−1) 海水 1000 1500 — 黏土 1650 1650 218 砂土 1800 1697 264 岩石 2100 2135 485 -
丁海平, 朱重洋, 于彦彦, 2017. P, SV波斜入射下凹陷地形地震动分布特征. 振动与冲击, 36(12): 88—92, 98 doi: 10.13465/j.cnki.jvs.2017.12.015Ding H. P. , Zhu C. Y. , Yu Y. Y. , 2017. Characteristic of ground motions of a canyon topography under inclined P and SV waves. Journal of Vibration and Shock, 36(12): 88—92, 98. (in Chinese) doi: 10.13465/j.cnki.jvs.2017.12.015 董阳, 2016. 谱元法在海洋声学中的应用. 哈尔滨: 哈尔滨工程大学.Dong Y. , 2016. The spectral element method applied to ocean acoustic. Harbin: Harbin Engineering University. (in Chinese) 孔曦骏, 李鸿晶, 2021. 基于谱元模型的带斜坡覆水场地地震动效应分析. 南京工业大学学报(自然科学版), 43(2): 264—272 doi: 10.3969/j.issn.1671-7627.2021.02.018Kong X. J. , Li H. J. , 2021. Spectral element analysis of seismic effects for overlying water fields with side slopes. Journal of Nanjing University of Technology (Natural Science Edition), 43(2): 264—272. (in Chinese) doi: 10.3969/j.issn.1671-7627.2021.02.018 李颂, 杨树春, 仝志刚等, 2012. 南海南部礼乐盆地深水区烃源岩生烃潜力研究. 天然气地球科学, 23(6): 1070—1076Li S. , Yang S. C. , Tong Z. G. , et al. , 2012. Study on the hydrocarbon generation potential of source rocks in Lile Basin, Southern South China Sea. Natural Gas Geoscience, 23(6): 1070—1076. (in Chinese) 李伟华, 赵成刚, 2003. 圆弧形凹陷饱和土场地对平面P波散射问题的解析解. 地球物理学报, 46(4): 539—546Li W. H. , Zhao C. G. , 2003. An analytical solution for the diffraction of plane P-waves by circular cylindrical canyons in a fluid-saturated porous media half space. Chinese Journal of Geophysics, 46(4): 539—546. (in Chinese) 刘少林, 杨顶辉, 徐锡伟等, 2021. 模拟地震波传播的三维逐元并行谱元法. 地球物理学报, 64(3): 993—1005Liu S. L. , Yang D. H. , Xu X. W. , et al. , 2021. Three-dimensional element-by-element parallel spectral-element method for seismic wave modeling. Chinese Journal of Geophysics, 64(3): 993—1005. (in Chinese) 彭浩天, 2018. 起伏地表下弹性波传播数值模拟方法对比研究. 成都: 西南石油大学. 苏波, 李怀良, 刘少林等, 2019. 修正辛格式有限元法的地震波场模拟. 地球物理学报, 62(4): 1440—1452Su B. , Li H. L. , Liu S. L. , et al. , 2019. Modified symplectic scheme with finite element method for seismic wavefield modeling. Chinese Journal of Geophysics, 62(4): 1440—1452. (in Chinese) 孙龙涛, 孙珍, 詹文欢等, 2010. 南沙海域礼乐盆地油气资源潜力. 地球科学——中国地质大学学报, 35(1): 137—145 doi: 10.3799/dqkx.2010.014Sun L. T. , Sun Z. , Zhan W. H. , et al. , 2010. Petroleum potential prediction of the Lile basin in Nansha. Earth Science—Journal of China University of Geosciences, 35(1): 137—145. (in Chinese) doi: 10.3799/dqkx.2010.014 谭文卓, 吴帮玉, 李博等, 2020. 梯形网格伪谱法地震波场模拟. 石油地球物理勘探, 55(6): 1282—1291 doi: 10.13810/j.cnki.issn.1000-7210.2020.06.014Tan W. Z. , Wu B. Y. , Li B. , et al. , 2020. Seismic wave simulation using a trapezoid grid pseudo-spectral method. Oil Geophysical Prospecting, 55(6): 1282—1291. (in Chinese) doi: 10.13810/j.cnki.issn.1000-7210.2020.06.014 万子轩, 2020. 基于谱元法的山体地形效应模拟研究. 成都: 成都理工大学.Wan Z. X. , 2020. Investigations of mountain topographic effects based on spectral element method. Chengdu: Chengdu University of Technology. (in Chinese) 魏成前, 2021. 平面波入射下二维成层盆地地震动放大特征研究. 苏州: 苏州科技大学.Wei C. Q. , 2021. Research on characteristics of ground motion amplification in two-dimensional layered basins under plane wave incidence. Suzhou: Suzhou University of Science and Technology. (in Chinese) 姚铭, 高刚, 周游等, 2017. 基于有限差分法的地震波数值模拟研究综述. 能源与环保, 39(10): 75—79, 85Yao M. , Gao G. , Zhou Y. , et al. , 2017. Summarization of numerical simulation of seismic wave based on minite difference method. China Energy and Environmental Protection, 39(10): 75—79, 85. (in Chinese) 禹乐, 于彦彦, 丁海平, 2020. 内源作用下盆地倾角对地表地震动放大特征的影响. 地震工程与工程振动, 40(5): 97—106 doi: 10.13197/j.eeev.2020.05.97.yul.010Yu L. , Yu Y. Y. , Ding H. P. , 2020. Effect of basin slope angle on ground motion amplification characteristics under internal point source. Earthquake Engineering and Engineering Vibration, 40(5): 97—106. (in Chinese) doi: 10.13197/j.eeev.2020.05.97.yul.010 张田升, 吴自银, 赵荻能等, 2019. 南海礼乐盆地海底麻坑地貌及成因分析. 海洋学报, 41(3): 106—120Zhang T. S. , Wu Z. Y. , Zhao D. N. , et al. , 2019. The morphologies and genesis of pockmarks in the Reed Basin, South China Sea. Haiyang Xuebao, 41(3): 106—120. (in Chinese) 赵成刚, 王磊, 李伟华, 2008. 具有饱和土沉积层的充水河谷对平面瑞雷波的散射. 地球物理学报, 51(5): 1567—1573 doi: 10.3321/j.issn:0001-5733.2008.05.032Zhao C. G. , Wang L. , Li W. H. , 2008. Scattering of plane Rayleigh waves by circular-arc alluvial valleys with saturated soil deposits and water layer. Chinese Journal of Geophysics, 51(5): 1567—1573. (in Chinese) doi: 10.3321/j.issn:0001-5733.2008.05.032 朱重洋, 丁海平, 于彦彦, 2016. SV波入射下有无覆盖层凹陷地形的地面运动比较. 苏州科技学院学报(工程技术版), 29(2): 33—37Zhu C. Y. , Ding H. P. , Yu Y. Y. , 2016. Comparison of canyon ground motion with or without covering soil layers under SV wave incidence. Journal of Suzhou University of Science and Technology (Engineering and Technology), 29(2): 33—37. (in Chinese) Antonietti P. F. , Ferroni A. , Mazzieri I. , et al. , 2018. Numerical modeling of seismic waves by discontinuous spectral element methods. ESAIM: Proceedings and Surveys, 61: 1—37. doi: 10.1051/proc/201861001 Chen B. K. , Du Y. J. , Shi Y. , et al. , 2021. Seismic analysis of isolated continuous bridge considering influence of seawater and site condition. Shock and Vibration, 2021: 7599715. Cristini P. , Komatitsch D. , 2012. Some illustrative examples of the use of a spectral-element method in ocean acoustics. The Journal of the Acoustical Society of America, 131(3): EL229—EL235. doi: 10.1121/1.3682459 Graves R. W. , Pitarka A. , Somerville P. G. , 1998. Ground-motion amplification in the Santa Monica area: effects of shallow basin-edge structure. Bulletin of the Seismological Society of America, 88(5): 1224—1242. doi: 10.1785/BSSA0880051224 Komatitsch D. , Barnes C. , Tromp J. , 2000. Wave propagation near a fluid-solid interface: a spectral-element approach. Geophysics, 65(2): 623—631. doi: 10.1190/1.1444758 Lee J., 2013. Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method. Physics and Chemistry of the Earth, Parts A/B/C, 63: 102—115. Patera A. T. , 1984. A spectral element method for fluid dynamics: laminar flow in a channel expansion. Journal of Computational Physics, 54(3): 468—488. doi: 10.1016/0021-9991(84)90128-1 期刊类型引用(1)
1. 赵煜,魏毛毛,孙艳萍,史一彤,陈文凯. Bootstrap抽样下基于贝叶斯估计的地震人员死亡评估研究. 地震研究. 2025(01): 132-140 . 百度学术
其他类型引用(4)
-