An Improved Approach for Near-fault Pulse-like Ground Motion Simulation
-
摘要: 本文基于小波包技术的随机地震动模拟方法，提出一种改进的参数化随机近断层脉冲型地震动模拟方法。然后，通过识别和提取近断层脉冲型地震动数据库中脉冲型地震动的特征参数，建立了基于震源、传播路径和场地特征等参数的脉冲模型参数预测方程。最后，通过模拟实际记录和误差分析检验了改进的模拟方法的有效性。结果表明：应用改进的模拟方法得到的地震动时程无论在波形、频率特性还是峰值上均与实际记录具有较好的一致性。改进的模拟方法在保留地震动时频非平稳性的基础上，能够有效地提高近断层脉冲型地震动的模拟效果，并且能够很好地体现脉冲型地震动的主要特征。Abstract: In this paper, an improved parameterized stochastic simulation method of near-fault pulse-like ground motion is proposed based on the stochastic ground motion simulation method using the wavelet packet method. Then, by identifying and extracting the characteristic parameters of pulse-like ground motions in a set near-fault ground motions database, the prediction equation of pulse model parameters in terms of the earthquake source, path and site characteristic parameters is established. Finally, the improved simulation method is verified by simulating recorded pulse-like ground motions and error analysis. Our results illustrate that simulated time histories of pulse-like ground motions by improved simulation method is well consistent with the recorded one in terms of waveform, frequency characteristic and peak value. The improved simulation method can effectively improve the simulation results of near-fault ground motion under retain time-frequency nonstationary of ground motion, and can also well reflect the main characteristics of pulse-like ground motion, which has very important engineering application value.
-
引言
镇江地区历史上发生过多次中强震，地震预警系统的建设有一定意义，2015年12月镇江市地震预警系统建成并投入使用，该系统分为2个主要部分，即地震波接收台点（下文称烈度台点）、预警信号发出台点（下文称预警终端）。但该系统建成后，镇江地区未发生过中强震，有感地震震级也较小，仅在很小的范围产生影响，预警系统并没有产生一定作用，因此该系统使用方（目前主要是学校师生）对预警系统的实际应用效果并不了解，地震部门对预警终端的增补也缺少相应的依据。为了解决这些问题，本文根据镇江地区预警台点分布情况及地震动峰值加速度和地震烈度衰减关系，从理论上对该系统的预警效能进行分析，确定目前本市地震预警系统的布点合理性及不同地点布设的预警终端对本市地震的理论预警时间。
1. 地震预警系统
地震预警指地震发生后，对震中附近一定距离区域尚未感知到但即将有可能感知到地震的民众发布的警报。地震波从地震发生地（震源）传播到周边一定范围，根据距离远近需要几秒到几十秒，如果在地震发生地布设地震探测仪器，地震发生时仪器最先接收纵波，经过短暂时间，经计算机自动处理，估算出地震三要素，通过电磁波信号以比横波更快的速度（3×105km/s）发送给离震中一定范围外的民众，达到地震预警目的。
地震预警系统虽能发挥一些社会效用，但仍存在一些亟待解决的问题，主要包含技术层面与非技术层面的问题。技术层面主要为预警信号发出区离震中太近存在盲区，无法预警，而离太远发出的预警信号无实际意义。本文主要分析技术层面的问题，同时结合镇江地区地震背景和各台点建设情况，分析各预警终端有效预警时间。
地震预警时间计算方法可用下式表示（郭凯等，2016）：
$$T_{\mathrm{war}}=\frac{S_{\mathrm{hyp}}}{v_{\mathrm{s}}}-\frac{S_{\mathrm{sta}}}{v_{\mathrm{p}}}-T_{\mathrm{dat}}-T_{\mathrm{cen}}-T_{\mathrm{iss}} $$ (1) 式中，Twar为地震预警时间，Shyp为震源与目标区域的距离，Ssta为震源到台站的距离（假设3个台触发阈值才发出预警信息），vs为S波波速，vP为P波波速，Tdat为预警系统计算所需P波数据时长，Tcen为数据传输到数据处理中心的时延和系统计算时间之和；Tiss为预警信息发布时间，震源视为点源。
纵波速度快，离台站最近，最先到达烈度台点，而破坏性较大的横波速度慢，离预警地区较远，较迟到达预警地区（图 1）。
当地震预警时间Twar＞0时，预警有效；当Twar≤0时，预警无效，根据式（1）可以计算出预警盲区，在震中附近一定区域无法预警。
但如果距震中过远，地震波到达时造成的地震烈度等级较低，有的甚至无法为人所感知，这种预警无实际意义。
2. 镇江市及周边地区地震预警系统建设情况
目前，镇江地区共有烈度台点17个，预警终端8个。烈度台点装有烈度仪，仪器感知地震动时向数据处理中心发出地震动信号，但该台站并不对外发出预警信号，预警信号主要由预警终端发出，该终端配备2-3个60w的喇叭，当接收到预警信号时，通过喇叭向使用者发出报警。为便于维护和管理，烈度仪安装在各辖市区原有地震台、强震台、宏观观测点及各乡镇镇政府内。预警终端在各辖市区选择1-2所学校布设，测试使用。
3. 预警盲区与台点间距
假设需3个台点触发阈值，计算机自动处理并开始计算地震三要素，最好的情况是地震发生地点在3个台构成的三角形中心下方。因研究区域较小，将地表近似视为平面，假设台点间距相同，3个台点构成等边三角形，则：
$$S_{\mathrm{sta}}^{2}=h^{2}+\frac{1}{3} a_{0}^{2} $$ (2) 式中，a0为台点间距，h为震源深度。
根据式（1），当Twar≤0时无法预警，设Twar=0，则式（1）变为：
$$S_{\mathrm{hyp}}=v_{\mathrm{s}}\left(\frac{S_{\mathrm{sta}}}{v_{\mathrm{p}}}+T_{\mathrm{dat}}+T_{\mathrm{cen}}+T_{\mathrm{iss}}\right) $$ (3) 各参数可赋值为：Tdat=3s，Tcen=2s，Tiss=1s，vP=5.8km/s，vs=3.6km/s。
设震中距为Δepi，则：
$$\mathit{\Delta }_{\text {epi }}=\sqrt{S_{\text {hyp }}^{2}-h^{2}} $$ (4) 可计算出不同台网密度（台点间距）对应的预警盲区（震中距），计算结果如图 2所示。
图 2 台点间距与预警盲区关系曲线Figure 2. Relation curves between station density and warning blind area计算得出台点间距达1km时，预警盲区为25.96km，随台点间距增大，盲区扩大。
根据镇江市烈度台点分布，台点间距近似用下式计算：
$$a_{0}=\frac{a_{12}+a_{22}+, \ldots, +a_{172}}{34} $$ (5) 其中，a12+a+2+, …, +a172为各台到最近2个台点的距离和。
台点间距离计算方法如下：
设两点经纬度分别为A（WA，JA）、B（WB，JB），根据大地坐标系与直角坐标系换算关系（孔祥元等，2001），A点直角坐标为：
$$\left\{\begin{array}{l} X_{\mathrm{A}}=\frac{a \cos W_{\mathrm{A}}}{W} \cos J_{\mathrm{A}} \\ Y_{\mathrm{A}}=\frac{a \cos W_{\mathrm{A}}}{W} \sin J_{\mathrm{A}} \\ Z_{\mathrm{A}}=\frac{b \sin W_{\mathrm{A}}}{V} \end{array}\right. $$ (6) B点直角坐标为：
$$\left\{\begin{array}{l} X_{\mathrm{B}}=\frac{a \cos W_{\mathrm{B}}}{W} \cos J_{\mathrm{B}} \\ Y_{\mathrm{B}}=\frac{a \cos W_{\mathrm{B}}}{W} \sin J_{\mathrm{B}} \\ Z_{\mathrm{B}}=\frac{b \sin W_{\mathrm{B}}}{V} \end{array}\right. $$ (7) 式中，a为椭圆长半轴，b为椭圆短半轴，W、V为含椭圆偏心率的辅助函数。
为简化计算，将地球考虑为正球体而非椭球体，W、V取值为1，式（1）、式（2）中a、b则为地球半径，因a、b仅参与弧度计算，可将a、b假设为1。
AB间弧长即地面AB两点间的球面距离△AB：
$${\mathit{\Delta }_{{\rm{AB}}}} = {\rm{R}}C = {\rm{R}} \times {\mathop{\rm acos}\nolimits} \left({\frac{{a_1^2 + b_1^2 - c_1^2}}{{2{a_1}{b_1}}}} \right) $$ (8) 式中，C为OA与OB直线间夹角（弧度），R为地球半径，a1为A点到球心的距离，b1为B点到球心的距离，c1为A、B点直线距离。
$$c_{1}^{2}=\left(X_{\mathrm{A}}-X_{\mathrm{B}}\right)^{2}+\left(Y_{\mathrm{A}}-Y_{\mathrm{B}}\right)^{2}+\left(Z_{\mathrm{A}}-Z_{\mathrm{B}}\right)^{2} $$ (9) 地球半径R取6371.004km。计算得出镇江市烈度台点平均台点间距a0为16.3km，预警盲区为28.4km，即预警系统对震中28.4km范围内的地区无法预警。如果台点间距缩至10km，盲区范围为震中27km范围内，盲区面积缩小10%，缩小范围较小。总体说来，镇江市烈度速报台站台点间距较合理。
4. 预警有效范围
当预警地区地震动峰值加速度或地震烈度达到一定程度时，认为地震预警有效果。根据《中国地震烈度表》（GB/T 17742-2008）（中华人民共和国国家质量监督检验检疫总局等，2009），地震动峰值加速度达22cm/s2或地震烈度达Ⅴ度时，室内绝大多数人有感，房屋颤动，不稳定器物摇动或翻倒，因东部沿海地区民众对地震关注度不高，容忍程度较低，地震烈度为Ⅴ度时预警已有一定意义。根据东部地区地震动峰值加速度（PGA）衰减关系（王国新等，2014）和中强地震活动区地震烈度衰减关系（卢建旗等，2009），可计算不同震级时地震动峰值加速度和地震烈度与震中距的关系如下：
$$\ln \left(Y_{1}\right)=5.304+1.7196 M_{S}-2.5903 \ln \left(L+2.789 \mathrm{e}^{0.451 M_{\mathrm{S}}}\right) $$ (10) $$\ln \left(Y_{2}\right)=2.726+1.3467 M_{S}-1.7636 \ln \left(L+1.046 \mathrm{e}^{0.451 M_{\mathrm{S}}}\right) $$ (11) $$I_{\mathrm{a}}=5.841+1.071 M_{S}-3.657 \log (L+15) $$ (12) $$I_{\mathrm{b}}=3.944+1.071 M_{S}-2.845 \log (L+7) $$ (13) 其中，Y1、Y2分别为长、短轴地震动峰值加速度，Ia、Ib分别为长、短轴地震烈度，MS为震级，L为震中距（km）。
发生MS5.0以下地震时，在前文所述的28.4km盲区外，地震烈度达不到Ⅴ度，地震动峰值加速度达不到22cm/s2，对预警来说意义不大。因此根据长轴公式（10）、（12）绘制MS5.0、MS5.5、MS6.0地震时地震动峰值加速度、地震烈度和震中距关系曲线图（图 3）。由图 3可知，震中距越大，地震动峰值加速度越小，地震烈度越小；当发生MS5.0地震时，震中距38km内地震动峰值加速度可达22cm/s2，震中距34km内地震烈度达Ⅴ度；当发生MS5.5地震时，震中距57km内地震动峰值加速度可达22cm/s2，震中距54km内地震烈度达Ⅴ度；当发生MS6.0地震时，震中距84km内地震动峰值加速度可达22cm/s2，震中距82km内地震烈度达Ⅴ度。因此，本预警系统对于MS5.0地震预警有效范围为34-38km，对于MS5.5地震预警有效范围为54-57km，对于MS6.0地震预警有效范围为82-84km。
图 3 地震动峰值加速度和地震烈度衰减关系曲线Figure 3. Curves of seismic peak acceleration and seismic intensity attenuation laws with various magnitudes5. 镇江市地震预警系统预警有效性分析
镇江市处于华北地震区长江下游-黄海地震带，属中强地震活动带。根据《中国历史强震目录（公元前23世纪-公元1911年）》（国家地震局震害防御司，1995）、《中国近代地震目录（公元1912年-1990年，MS≥4.7）》（中国地震局震害防御司，1999），镇江行政区域内共发生4次破坏性地震（MS≥4.7)，分别为1630年2月4日江苏镇江句容$4 \frac{3}{4} $级地震、1872年7月24日江苏镇江西$4 \frac{3}{4}$级地震、1913年4月3日江苏镇江MS5.5地震、1930年1月3日江苏镇江MS5.5地震。因此，镇江地区存在发生MS5.0以上地震的背景。
根据1970年以来的中国台网地震目录，1970年至今，镇江市虽未发生过MS5.0以上地震，但在距镇江行政边界90km范围内，发生过2次MS5.0以上地震（1974年4月22日溧阳MS5.5地震、1979年7月9日溧阳MS6.0地震）、2次接近MS5.0地震（1979年7月11日溧阳MS4.7地震、2012年7月20日高邮MS4.9地震）。因此，在镇江及周边地区建设地震预警系统具有一定意义。
计算有意义的最大预警时间时，假设地震发生在距各预警终端最远的3个烈度台点中心点，计算步骤如下：
（1）根据各烈度台点经纬度得出每3个烈度台点中心点经纬度，再根据各预警终端经纬度求出各预警终端与各中心点的距离，可假设离各预警终端最远的中心点为震中，计算出震中与预警终端的震中距Δepi，根据震源深度h（本文取10km）计算出震源距Shyp；
（2）假设地震P波到达离其最近的第3个台站时开始进行数据处理，则根据震中经纬度（3个烈度台点中心点经纬度）和3个烈度台点经纬度，求得震中与烈度台点的距离，取其最大值，即为震中距Δsta，再根据震源深度h（本文取10km）计算出震源距Ssta；
（3）根据公式（1）计算出最大预警时间Tmax。
（4）引入前述预警有效范围，当地震发生时，如果Δepi超过该震级对应的最大预警有效范围ΔMS，则Δepi取值为该震级对应的最大预警有效范围ΔMS，Ssta可根据公式（2）取本市平均台点间距近似计算，再根据公式（1）计算出有意义的最大预警时间Tmea-MS。如果Δepi＞ΔMS，则有意义的最大预警时间Tmea-MS=Tmax。
计算结果见表 1。
表 1 各地震预警终端有效预警时间Table 1. Effective earthquake early warning time of the earthquake early warning terminals预警终端 Δsta3/km Δepi/km Tmea-Ms5.0/s Tmea-Ms5.5/s Tmea-Ms6.0/s 镇江市地震台 17.3 49.2 2.5 4.5 4.5 镇江实验学校万科魅力之城分校 17.3 53.5 2.5 5.7 5.7 江滨实验小学 17.3 57.8 2.5 7.7 6.9 辛丰小学 17.3 55.8 2.5 6.3 6.3 丹阳华南实验学校 17.3 46.8 2.5 3.9 3.9 扬中市外国语小学 17.3 80.8 2.5 7.7 13.1 句容实验小学 16.0 61.8 2.5 7.7 8.1 镇江新区实验小学 17.3 68.7 2.5 7.7 9.8 根据计算结果，当镇江市发生MS5.0地震时，因各地震预警终端最远震中距在该震级预警有效范围外，因此计算得出的有效预警时间相同，即可提前0-2.5s预警本地地震烈度在Ⅴ度或地震动峰值加速度达22cm/s2的地震；当镇江市发生MS5.5地震时，4个预警终端最远震中距大于震级预警有效范围，可提前0-7.7s预警本地地震烈度在Ⅴ度或地震动峰值加速度达22cm/s2的地震，另外4个预警终端最多可提前3.9-6.3s预警本地地震烈度在Ⅴ度或地震动峰值加速度达22cm/s2的地震；当镇江市发生MS6.0地震时，各预警终端最远震中距均在该震级预警有效范围内，最多可提前3.9-13.1s预警本地烈度在Ⅴ度或地震动峰值加速度达22cm/s2的地震。
6. 结论
（1）根据三台触发最优结果为地震发生在三台中心点的原则，台点间距为1km时，预警盲区为25.96km，且预警盲区随台点间距的增大而增大。根据镇江市各烈度台点至最近的2个烈度台点的距离，计算得出镇江市地震烈度台点平均台点间距为16.3km，预警盲区为28.4km。
（2）根据地震动峰值加速度和地震烈度随震中距的衰减关系，计算得出发生不同震级地震时地震动峰值加速度和地震烈度与震中距的关系。为使当地地震预警具有一定意义，本文假设当地地震动峰值加速度达22cm/s2或地震烈度达Ⅴ度时预警有效，计算知：发生MS5.0以下地震时，在28.4km盲区外，地震烈度达不到Ⅴ度，地震动峰值加速度达不到22cm/s2，对预警来说意义不大；当发生MS5.0地震时，震中距38km内地震动峰值加速度能达22cm/s2，震中距34km地震内烈度达Ⅴ度；当发生MS5.5地震时，震中距57km内地震动峰值加速度能达22cm/s2，震中距54km内地震烈度达Ⅴ度；当发生MS6.0地震时，震中距84km内地震动峰值加速度能达22cm/s2，震中距82km内地震烈度达Ⅴ度。
（3）根据各烈度台点和预警终端经纬度，再引入预警有效范围，可计算得到发生不同震级地震时，不同地点预警终端最多可提前2.5-13.1s预警本地地震烈度为Ⅴ度或地震动峰值加速度达22cm/s2的地震。
-
图 1 基于小波包的地震动模拟方法流程
Figure 1. Procedures for simulating ground motions using wavelet packet method
图 3 1979年Imperial Valley地震NGA#180水平向脉冲地震动模拟过程
Figure 3. Simulation process of NGA#180 horizontal pulse-like ground motions from Imperial Valley earthquake in 1979
图 4 NGA#182和NGA#1084的实际观测和模拟地震动时程
Figure 4. Acceleration, velocity and displacement time histories of observed and simulated ground motions
图 5 NGA#182和NGA#1084的实际观测和模拟地震动的加速度反应谱（5%阻尼比）
Figure 5. Acceleration response spectra at 5% damping ratio of observed and simulated ground motions
图 6 脉冲地震数据库中实际观测记录与模拟时程的加速度反应谱的残差均值及±1倍标准差
Figure 6. Mean residual bias and the±one standard deviation of near-fault stations observed and simulated in the pulse-like ground motions database
表 1 拟合参数γ和υ的分布系数和限定条件
Table 1. Boundary condition and coefficient of marginal distributions fitted to parameters γ and υ
参数 单位 符合分布 选定下限 选定上限 合适的分布参数 均值 标准差 γ burr 2.0 3.0 α=2.49，c=21.80，k=4.20 2.3 0.14 υ π uniform 0 2 1.0 0.58 
下载: 导出CSV
表 2 MP03脉冲模型参数的回归系数和标准差
Table 2. Regression coefficients and standard deviations of the MP03 model parameters
参数 α0 α1 α2 β1 β2 β3 β4 β5 ϕ τ σ Vp 1.6784 1.6152 1.6643 0.2201 -0.3060 -0.3230 -0.1231 0.0054 0.132 0.108 0.330 Tp -0.2312 -0.1828 0.1677 0.1762 -0.1295 0.0043 0.261 0.136 0.294 t0 -0.8739 -0.7214 -0.4478 0.2708 0.0551 -0.0455 0.0028 0.285 0.089 0.308
下载: 导出CSV
-
樊剑, 涂家祥, 吕超等, 2008.采用时频滤波技术的近断层脉冲地震人工模拟.华中科技大学学报(自然科学版), 36(11):116-119. http://www.cnki.com.cn/Article/CJFDTotal-HZLG200811044.htm 李亚楠, 2016.工程用地震动模拟随机性方法研究.大连: 大连理工大学. http://cdmd.cnki.com.cn/Article/CDMD-10141-1016204382.htm 罗全波, 陈学良, 高孟潭等, 2018.近断层速度脉冲与震源机制的关系浅析.震灾防御技术, 13(3):646-661. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180316&journal_id=zzfyjs 田玉基, 杨庆山, 卢明奇, 2007.近断层脉冲型地震动的模拟方法.地震学报, 29(1):77-84. http://d.old.wanfangdata.com.cn/Periodical/dizhen200701009 Atkinson G. M., Assatourians K., Boore D. M., et al., 2009. A guide to differences between stochastic point-source and stochastic finite-fault simulations. Bulletin of the Seismological Society of America, 99(6):3192-3201. doi: 10.1785/0120090058 Baker J. W., 2007. Quantitative classification of near-fault ground motions using wavelet analysis. Bulletin of the Seismological Society of America, 97(5):1486-1501. doi: 10.1785/0120060255 Beresnev I. A., Atkinson G. M., 1997. Modeling finite-fault radiation from the ωn spectrum. Bulletin of the Seismological Society of America, 87(1):67-84. Boore D. M., 1983. Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra. Bulletin of the Seismological Society of America, 73(6A):1865-1894. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=0aecda90-ca1d-4056-aa8d-57610c62fdb5 Bray J. D., Rodriguez-Marek A., 2004. Characterization of forward-directivity ground motions in the near-fault region. Soil Dynamics and Earthquake Engineering, 24(11):815-828. doi: 10.1016/j.soildyn.2004.05.001 Campbell K. W., 2008. Hybrid empirical ground motion model for PGA and 5% damped linear elastic response spectra from shallow crustal earthquakes in stable continental regions: Example for eastern North America. In: Proceeding of the 14th World Conference Earthquake Engineering. Beijing, China: WCEE. Dabaghi M. N., 2014. Stochastic modeling and simulation of near-fault ground motions for performance-based earthquake engineering. Berkeley: University of California. https://escholarship.org/uc/item/8db18756 Dabaghi M., Der Kiureghian A., 2017. Stochastic model for simulation of near-fault ground motions. Earthquake Engineering & Structural Dynamics, 46(6):963-984. http://cn.bing.com/academic/profile?id=080c7eb3415d4053e615dc23bd911f6a&encoded=0&v=paper_preview&mkt=zh-cn Joyner W. B., Boore D. M., 1993. Methods for regression analysis of strong-motion data. Bulletin of the Seismological Society of America, 83(2):469-487. http://cn.bing.com/academic/profile?id=862575101d927d3d6b3b88e6526957a0&encoded=0&v=paper_preview&mkt=zh-cn Li Y. N., Wang G. X., 2016. An improved approach for nonstationary strong ground motion simulation. Pure and Applied Geophysics, 173(5):1607-1626. doi: 10.1007/s00024-015-1189-4 Mavroeidis G. P., Papageorgiou A. S., 2003. A mathematical representation of near-fault ground motions. Bulletin of the Seismological Society of America, 93(3):1099-1131. doi: 10.1785/0120020100 Motazedian D., Atkinson G. M., 2005. Stochastic finite-fault modeling based on a dynamic corner frequency. Bulletin of the Seismological Society of America, 95(3):995-1010. doi: 10.1785/0120030207 Somerville P. G., Smith N. F., Graves R. W., et al., 1997. Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture directivity. Seismological Research Letters, 68(1):199-222. doi: 10.1785/gssrl.68.1.199 Yamamoto Y., Baker J. W., 2013. Stochastic model for earthquake ground motion using wavelet packets. Bulletin of the Seismological Society of America, 103(6):3044-3056. doi: 10.1785/0120120312 Yang D. X., Zhou J. L., 2015. A stochastic model and synthesis for near-fault impulsive ground motions. Earthquake Engineering & Structural Dynamics, 44(2):243-264. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=720a58f1cf843d0534c204535de9022c 期刊类型引用(4)
1. 闫新义，陈勇. 2022年7月20日昌吉M_S4.8地震预警处理结果分析. 内陆地震. 2024(03): 223-230 . 
百度学术2. 闫新义，李艳永，赵石柱，陈勇. 新疆天山中段地震预警结果分析. 内陆地震. 2023(01): 44-51 . 百度学术3. 王亚玲，毛国良，蒋宏毅，王宁，蔡玲玲. 河北预警网地震预警处理系统运行现状. 华北地震科学. 2023(03): 43-51 . 百度学术4. 闫新义，李艳永. 新疆地震预警系统预警处理能力分析——以拜城M_S5.4地震为例. 内陆地震. 2021(02): 143-149 . 百度学术其他类型引用(1)
-
下载:
点击查看大图
计量
- 文章访问数: 184
- HTML全文浏览量: 26
- PDF下载量: 13
- 被引次数: 5
