Coda wave Q values of shallow media in Binchuan basin using dense array data
-
摘要:
为得到更准确的浅部介质衰减参数，便于工程地震领域开展更精细的地震动模拟等研究，以气枪密集台阵资料为基础，联合天然地震资料，利用Sato单散射模型研究了云南宾川盆地浅层介质尾波衰减特征。结果显示，随着频率的增加，尾波Q值Qc整体呈增加趋势，符合Q值的频率依赖关系。在空间上，该地区Qc值分布具有明显的横向不均匀性，位于研究区中部宾川盆地的台站Qc值较低，而位于研究区南西和北东方向的山地丘陵地区的台站Qc值较高，与速度层析成像研究结果一致。研究区平均Qc频率依赖关系为$ { {Q}_{\mathrm{c}} （ f ） ＝28.04{f}^{1.07}}$，气枪震源密集台阵资料得到的Q0值较天然地震得到的值更低，证明使用密集台阵资料得到的结果反映了更浅部介质的衰减特征，而更高的频率依赖指数$\, \eta \,$值意味着浅层介质的非均匀性高于深部介质，符合实际情况。宾川盆地Q0值大于松辽盆地、华北盆地和中国大陆平均沉积层，而指数$\, \eta \,$小于松辽盆地和中国大陆，表明本研究所用资料反映的Q值信息介于近地表与深部介质之间。随着尾波流逝时间窗口的增大，Q0逐渐增大，而频率依赖性降低。因此，在研究区尺度较小且选择的天然地震事件震级不高的情况下，需要选取较小的流逝时间窗口以确保满足研究需要。此外，局部地形和地震地质构造变化可能导致衰减参数和相应的标准差相差较大。使用气枪震源和近震资料配合小尺度密集台阵，可以得到较准确的浅层介质衰减参数，为工程应用和浅部结构探测提供了新的思路。
Abstract:In order to obtain more accurate attenuation parameters of shallow medium and carry out more detailed ground motion simulation research in the field of engineering earthquakes, this paper uses the Sato single scattering model to study the coda attenuation characteristics of shallow medium in the Binchuan basin based on the dense array data and combined with natural seismic data. The results showed that with the increase of frequency, the coda Q value Qc showed an increasing trend, which was in line with the frequency dependence of the Q value. The spatial distribution of Qc values in this area has obvious lateral inhomogeneity. The stations located in the Binchuan basin in the central part of the study area have lower Qc values, while the stations located in the hilly areas of the study area in the southwest and northeast directions have higher Qc values, which is consistent with the velocity tomography results. The frequency dependence of the average Qc in the study area is $ {Q}_{{\rm{c}}} （ f ） ＝28.04{f}^{1.07} $. The Q0 value obtained by the dense array data is lower than that obtained by the natural earthquake, which proves that the results obtained by the dense array data reflects the attenuation characteristics of the shallower medium. And a higher $ \eta $ value means that the inhomogeneity of the shallow medium is higher than that of the deep medium, which is in line with the actual situation. The Q0 value of the Binchuan basin is larger than the Q0 value of the Songliang basin, the North China basin and the average sedimentary layer of the Chinese mainland, while the index $ \eta $ is smaller than the $ \eta $ value of the Songliao basin and the Chinese mainland, indicating that the information reflected by the dense array is between the near-surface and deep media. As the depth increases, the inhomogeneity of the medium gradually decreases, and the dependence of the Q value with frequency gradually weakens. Q0 increases with the lapse time window of coda, and the index $ \eta $ is the opposite. It is necessary to select a small lapse time window to ensure the accuracy of the results when the scale of study area is small and the magnitude of the selected natural earthquake events not large. In addition, local topographic and seismic-tectonic changes may lead to large differences in attenuation parameters and standard deviations. The useage of airgun source and near-field earthquake with small-scale dense seismic arrays can obtain more accurate attenuation parameters in shallow media, which provides a new idea for engineering application and shallow structure detection.
-
引言
圆柱型储罐因其具有容量大、用钢量少、占地面积小、结构和施工简单等优点已成为我国战略储备石油的主要设备。我国地震频发且烈度高，储罐在其作用下会随机运动，因而诱发流体剧烈的非线性晃荡（孙颖，2012），伴随的强烈抨击荷载和倾覆力矩影响着储罐的在位稳定性，储罐一旦破损，不仅会造成经济损失，也会带来严重的次生灾害，如环境污染和火灾事故（方浩等，2012）。因此，厘清地震载荷下储罐内流体的水动力特性对储罐的安全分析与结构优化设计具有重要的工程价值。
Virella等（2006）通过附加质量法分析了水平简谐载荷下圆柱型储罐内液体晃荡的脉冲模式，Chen等（2007）对圆柱型储罐在简谐载荷与地震载荷下的液体晃荡现象进行了数值研究和试验研究，其结果表明当外部载荷的主频接近于储液的固有频率时，外部载荷会诱发流体共振。杨杰等（2016）使用附加质量法对水平地震载荷下的大型渡槽中的流体非线性晃荡进行了研究。Caron等（2018）通过物理模型试验并结合数值模拟研究了圆柱型储罐在共振简谐载荷下自由液面的晃动情况。周利剑等（2019）通过物理模型试验与弹簧-质量法探究了地震频率对圆柱型储罐中流体运动的影响，其结果表明长周期地震引起的波高明显大于其它类型地震，且晃荡模态为多阶振型的叠加。Cheng等（2019）以及Hejazi和Mohammadi （2019）分别基于势流理论和物理模型实验分析了地震载荷下的液体晃荡现象，结果表明影响晃荡波高的主导因素是地震峰值加速度（peak ground acceleration，缩写为PGA）和主要频率。罗东雨等（2020）应用有限元软件ADINA研究16×104 m3液化天然气（liquefied natural gas，缩写为LNG）储罐在地震载荷下的响应，其结果显示地震频率是影响储液晃荡的主要因素。Hernandez-Hernandez等（2021）试验研究了PGA对储罐内流体晃荡的影响，其结果显示，最大环向应力与PGA线性相关。Jin等（2021）基于Navier-Stokes数值模型研究了天然地震对晃荡波高的影响规律，发现地震峰值速度（peak ground velocity，缩写为PGV）对峰值波高起着决定性作用。尽管地震是多向耦合的随机振动，但部分研究显示竖向地震对液体晃荡的影响并不显著（刘勇，2007；陈贵清等，2012；张如林等，2017）。
上述研究主要针对简谐共振或低频地震载荷下的液体晃荡现象，且多以附加质量法和弹簧-质量法来分析流固耦合问题，关于地震频率、PGA与PGV等因素对晃荡响应的影响规律还有待进一步揭示。为此，本研究拟采用计算流体力学软件Fluent对流体晃荡进行数值模拟，建立储罐与液体的耦合模型，选用包括低、中、高频在内的15种地震信号，探究地震信号的频率、PGA和PGV对储液的晃荡模式以及不同位置水动力特性的影响规律，以期为大型储罐的优化设计提供一定的参考。
1. 数值模型与验证
本研究采用Fluent进行数值模拟研究。流体在地震载荷下的晃荡通过多参考系法（multiple reference frame，缩写为MRF）实现。假设储液不可压，罐壁为刚性薄壁（Luo et al，2016），选用流体体积法（volume of fluid，缩写为VOF）追踪自由液面（Hirt，Nichols，1981）。
1.1 控制方程
不可压流动由连续性方程和Navier-Stokes方程描述，具体如下：
$$ \nabla \, {\bullet} \,\, {\boldsymbol{u}}＝0 \text{，} $$ (1) $$\frac{\partial \boldsymbol{u}}{\partial t} ＋ \nabla \, {\bullet} \,\, （ \boldsymbol{u} \boldsymbol{u} ） ＝－\frac{1}{\rho} \nabla p ＋ \nabla \, {\bullet} \,\, \left[\mu\left(\nabla \boldsymbol{u} ＋ （ \nabla \boldsymbol{u} ） ^{\mathrm{T}}\right) \right] ＋ f_i \text{，} $$ (2) 式中：ρ为流体密度；$\nabla $为哈密顿算子；速度矢量u＝$ ［ $u， v， w$ ］ $，其中u，v和w分别为液体在x，y和z方向上的速度分量；p为压强；μ为运动黏度；fi为单位体积流体所受到的体积力。选用多参考系法实现储罐内流体晃荡模块，控制方程可参照Shu等（2018）。
自由液面采用流体体积法进行追踪，通过求解以下通量方程来实现：
$$\begin{gathered}\frac{\partial \alpha_{\mathrm{w}}}{\partial t} ＋ \boldsymbol{U} \, {\bullet} \,\, \nabla \alpha_{\mathrm{w}}＝0 \text{，} \alpha_{\mathrm{a}}＋\alpha_{\mathrm{w}}＝1 \text{，} \end{gathered}$$ (3) 式中，αa和αw分别为控制体中空气和水的体积分数，αa＝1表示单元内为空气，αw＝1表示单元内为水，0＜αa＜1表示单元内为水和空气。
1.2 数值模型准确性验证
Chen等（2007）和Xue等（2019）研究了1995年神户M7地震载荷（图1a）下的流体晃荡现象，其中圆柱型储罐的尺寸为：半径R＝30 cm，静水水深H＝10 cm。本文对该工况进行数值模拟。图1b给出了罐壁右侧波高的试验结果（Chen et al，2007）和数值解（Xue et al，2019）与本文数值解的对比。尽管本文数值解与Xue等（2019）的数值解、Chen等（2007）的试验结果在峰值处有差异，但最大误差不超过12.4%，三者间几乎没有相位差，整体吻合得较好，说明本文所采用的数值方法具有较高的精度，可用于后续研究。
![]()
图 1 1995年神户M7地震的位移时程（a）以及神户地震载荷下罐壁右侧波高的时程曲线对比（b）Figure 1. Time history of ground displacement of Kobe seismic excitation （a） and comparison of time histories of wave heights at the right sidewall of the cylindrical tank under 1997 M7 Kobe seismic excitation （b）2. 地震的选取
本研究选取12个地震事件，涉及15种地震信号，相应的地震参数列于表1，相应的加速度时程如图2所示。PGA从0.081g增长到1.79g；PGV最低为0.22 m/s，最高可达1.76 m/s。地震的频率成分根据PGA与PGV的比值划分为低频、中频与高频等三大类（Tso et al，1992；Kianoush，Ghaemmaghami，2011）：PGA/PGV＞1.2为高频，0.8＜PGA/PGV＜1.2为中频，PGA/PGV＜0.8为低频。因此，低频、中频和高频地震分别有7、5和3个。
表 1 本研究涉及地震的相关参数Table 1. The characteristics of considered seismic excitations in this study地震信号 频率类别 主要频率/Hz PGA/g PGV/（m·s−1） PGA/PGV 芦山53Elk-EW 低频 0.425，0.625 0.081 0.377 0.215 芦山53Elk-SN 低频 0.325，0.4 0.224 0.866 0.259 芦山53Yjb-EW 低频 0.658，0.703 0.075 0.232 0.323 埃尔森特罗（El Centro）-EW 低频 0.449，0.842 0.214 0.48 0.446 神户（Kobe） 低频 0.269，0.537，0.879 0.611 1.27 0.48 埃尔津（Erzican） 低频 0.42，0.563，0.798 0.52 0.84 0.619 北岭-1（Northridge-1） 低频 0.353，0.873，1.08 0.607 0.8 0.754 唐山 中频 0.65，0.8，1.2 1.268 1.584 0.8 康奎特（Concrete） 中频 0.49，0.6，0.98，1.1 0.63 0.758 0.83 埃尔森特罗（El Centro）-SN 中频 0.561，0.839，1.16 0.349 0.38 0.91 塔夫脱（Taft） 中频 0.11，0.37，0.61，1.2 1.79 1.76 1.02 曼杜西诺角（Cape Mendocino） 中频 0.868，1.43，2.14 1.45 1.26 1.15 汶川卧龙-EW 高频 0.622，1.38，2.35 0.97 0.585 1.66 柯依那（Koyna）-SN 高频 0.87，1.55 0.585 0.31 1.89 芦山53Yjb-SN 高频 0.841，2.545 0.473 0.22 2.15 注：第一列中的53Elk，53Yjb和卧龙为站点名称，EW和SN为地震信号的东西分量和南北分量。 ![]()
图 2 本文选取地震的加速度时程曲线Figure 2. Time histories of considered seismic accelerations used in this study3. 地震载荷下的晃荡特性
本节对地震载荷下圆柱型储罐内的晃荡模式和水动力特性进行讨论。液体的一阶固有频率满足如下公式（Chen et al，2007）：
$$\begin{split} \omega_1＝\sqrt{\frac{2 \lambda_1 g}{d} \tanh \left(\lambda_1 \frac{h}{R}\right)} \text{，} \end{split}$$ (4) 式中：$\lambda_1 $为一阶贝塞尔函数的一阶导数，$\lambda_1 $＝1.841。本文中，储罐高HT＝15 m，直径d＝15 m，静水水深h＝10 m （Djermane et al，2014）。根据式（4）计算出ω1＝1.54 rad/s。储罐东侧布置9个压力测点，标记为P1，P2，···，P9 （分别距底部0 m，1 m，5 m，6 m，8 m，9.5 m，10 m，10.5 m，11 m）；四只数值波高仪排布在罐体的东、南、西、北四个方向，分别标记为We，Ws，Ww和Wn，如图3所示。
![]()
图 3 圆柱型储罐的整体视图（a）和顶视图（b）Figure 3. The overall-view （a） and top-view （b） geometries of the cylindrical tank3.1 网格无关性验证
基于网格法的数值仿真需要评估网格数对计算精度的影响。本文选用15万、21万、25万、28万和31万等五种网格数进行网格无关性验证。图4给出了芦山53Elk-SN载荷下，不同网格数下东侧（We）和北侧（Wn）处波高的对比。可见：当网格数由15万增至25万时，波高略微增大，进一步增至31万时，波高曲线几乎重合。综合考虑计算效率与精度，本文选取25万网格数开展后续研究。
![]()
图 4 芦山53Elk-SN载荷下储罐东侧We （a）和北侧Wn （b）处不同网格数下波高的对比Figure 4. Comparison of time histories of the wave heights at We （a） and Wn （b） from various grids under Lushan 53Elk SN excitation3.2 模拟结果与讨论
本节重点分析地震频率、PGV和PGA对晃荡波高、水动压的影响。
3.2.1 晃荡波高
不同地震载荷下晃荡波高的时程曲线如图5所示。主频接近储液固有频率的地震信号会引发共振晃荡（张如林等，2017；Jin et al，2021）。例如：芦山53Elk-SN、神户以及北岭-1地震信号的主频分别为0.325 Hz，0.269 Hz和0.353 Hz，接近于储液的一阶固有频率（0.245 Hz），相应的波高时程如图5b，e，g所示。从该图可见，三个工况均为共振波且远大于非共振工况下的波高。15种工况下地震的主频、PGV和东侧（We）处相应的峰值波高列于表2。一般情况下，中、高频地震载荷引起的晃荡波高小于低频地震引起的波高，但也存在一些特例。例如：低频的芦山53Yjb-EW诱发的峰值波高为0.17 m，而主频较高的低频的埃尔津地震和中频的曼杜西诺角地震却激发了1.06 m和1.35 m的波高，分别为前者的6.2倍和8倍。由表1可知，虽然芦山53Yjb-EW的主频低，但其PGV仅为0.232 m/s，而埃尔津和曼杜西诺角的PGV却达到0.84 m/s和1.26 m/s，分别为前者的3.6倍和5.4倍。由此可见，PGV对晃荡波高的影响不可忽略，有必要对地震主频和PGV对晃荡波高的影响进一步分析。
![]()
图 5 不同地震载荷下储罐东侧We和北侧Wn处的波高时程曲线（a） 芦山53Elk-EW；（b） 芦山53Elk-SN；（c） 芦山53Yjb-EW；（d） 埃尔森特罗-EW；（e） 神户地震；（f） 埃尔津；（g） 北岭-1；（h） 唐山；（i） 康奎特；（j） 埃尔森特罗-SN；（k） 塔夫脱；（l） 曼杜西诺角；（m） 汶川卧龙-EW；（n） 柯依娜-SN；（o） 芦山53Yjb-SNFigure 5. Time histories of wave heights under various seismic excitations（a） Lushan53Elk-EW；（b） Lushan53Elk-SN；（c） Lushan53Yjb-EW；（d） El Centro-EW；（e） Kobe；（f） Erzican；（g） Northridge-1；（h） Tangshan；（i） Concrete；（j） El Centro-SN；（k） Taft；（l） Cape Mendocino；（m） WenchuanWolong-EW；（n） Koyna-SN；（o） Lushan53Yjb-SN表 2 地震主频、PGV和东侧We处相应的峰值波高Table 2. The dominant frequency and PGV of seismic excitations and corresponding maximum wave heights at east side We of the tank地震信号 主频/Hz PGV/（m·s−1） 峰值波高/m 地震信号 主频/Hz PGV/（m·s−1） 峰值波高/m 塔夫脱 0.11 1.76 4.15 埃尔森特罗-SN 0.561 0.38 0.484 神户 0.269 1.27 2.02 汶川卧龙-EW 0.622 0.585 0.26 芦山53Elk-SN 0.325 0.866 3.574 唐山 0.65 1.584 2.3 北岭-1 0.353 0.8 1.48 芦山53Yjb-EW 0.658 0.232 0.17 埃尔津 0.42 0.84 1.06 芦山53Elk-SN 0.841 0.22 0.39 芦山53Elk-EW 0.425 0.377 0.713 曼杜西诺角 0.868 1.26 1.35 埃尔森特罗-EW 0.449 0.48 0.664 柯依那-SN 0.87 0.31 0.334 康奎特 0.49 0.758 0.675 峰值波高与地震主频的关系如图6所示。不难发现，峰值波高随地震主频的升高逐渐减小，但也有例外。例如：芦山53Elk-SN地震，其主频接近储液一阶固有频率，共振晃荡引起较大的波高，尽管其PGV仅为塔夫脱地震的一半，但其波高仍接近塔夫脱地震激发的波高；主频较低的埃尔森特罗-SN激发了0.484 m的波高，而主频较高的曼杜西诺角地震和唐山地震却诱发了1.35 m和2.3 m的波高，为前者的2.79倍和4.75倍，主要原因是：尽管曼杜西诺角地震和唐山地震的主频高于埃尔森特罗-SN，但二者的PGV较大，分别是埃尔森特罗-SN的3.32倍和4.17倍。因此，PGV也是影响晃荡的主要因素之一。
![]()
图 6 储罐东侧We处峰值波高与地震主频的关系Figure 6. The relationship between the maximum wave height at east side We of the tank and dominant frequency of earthquake储罐东侧We处的峰值波高与PGV的关系如图7所示。整体而言，峰值波高随PGV的增大而上升，二者呈正相关，如具有最大PGV的中频塔夫脱地震，激发的晃荡波最剧烈，波高可达4.15 m。当地震载荷的PGV相近时，主频低的地震诱发的波高较大，如：北岭-1地震与埃尔津地震的PGV相近，分别为0.8 m/s和0.84 m/s，由于二者的主频分别为0.353 Hz和0.563 Hz，差异明显，主频低的北岭-1激发的峰值波高（1.48 m）为埃尔津地震（1.06 m）的1.4倍。尽管神户地震的PGV （1.27 m/s）与曼杜西诺角的PGV （1.26 m/s）相近，二者引起的峰值波高分别为2.02 m和1.35 m，前者为后者的1.5倍，这主要源于其主频差异明显，分别为0.269 Hz和0.868 Hz。
![]()
图 7 储罐东侧We处峰值波高与PGV的关系Figure 7. The relationship between the maximum wave height at the east side We of the tank and PGV塔夫脱地震的主频（0.11 Hz）偏离储液的固有频率，塔夫脱地震载荷下部分时刻的自由液面如图8所示，可以看出：液面晃荡幅度大，但未出现波浪破碎、液体飞溅等强非线性行为，形成驻波的波面比较光滑，在40.0 s时刻，爬升至最高处。塔夫脱地震载荷下自由液面的波高相位如图9所示，可见：自由液面位移在载荷方向上变化明显，在非载荷方向上并不显著；在0 s至15 s时段（图9a），Ws处的自由液面位移与We处结果呈线性关系，即正向与负向的振幅几乎一致，储液沿外部载荷方向运动，表现为平面波（Faltinsen et al，2003）；在15 s至45 s时段（图9b，c），自由液面呈非线性关系，不对称性加剧，正向振幅接近负向振幅的两倍。
![]()
图 8 塔夫脱地震载荷下的自由液面Figure 8. Free surface snapshots under excitation of Taft earthquake![]()
图 9 塔夫脱地震载荷下自由液面的波高相位图Figure 9. Phase diagrams of wave heights under excitation of Taft earthquake（a） 0—15 s；（b） 15—30 s；（c） 20—45 s多数工况符合上述规律，但也存在特例，如芦山53Elk-SN。尽管其PGV在所有地震信号中处于中间值，却激发了较周围工况更为剧烈的晃荡波，主要原因是其主频与储液的一阶固有频率接近且具有较大PGV。部分时刻的自由液面如图10所示，可见：液体在储罐内形成了共振行进波，在第40 s和第43 s时刻形成了波斗较大的波面，最终发生了明显的水跃，造成波高局部增大；第44 s波浪破碎且局部飞溅；随后，爬升的液体落回自由液面，与向左运动的液体一起冲击罐壁，并在第45 s爬升至最高处，波高为3.58 m。该工况包含强非线性的物理现象，如液面翻转、破碎等，增加了流场的随机性（Wu，Chen，2009）。
![]()
图 10 芦山53Elk-SN地震载荷下的自由液面Figure 10. Free surface snapshots under Lushan 53Elk-SN excitation以上讨论表明：地震主频也是影响自由液面响应的主因之一，晃荡波高随地震主频的升高逐渐减小，当主频与储液的一阶固有频率相近时，易引起强烈的非线性共振晃荡；峰值波高与PGV呈较强的正相关；最剧烈的工况一般出现在低频且具有较大PGV的地震载荷，在设计储罐时应考虑减晃装置。
3.2.2 水动压
非破碎晃荡波对罐壁产生的水动压分为脉冲压力和对流压力两种模式（薛米安等，2019）。脉冲压力是由地震载荷直接导致的瞬时冲击压力，对流压力与自由液面运动密切相关。图11为峰值水动压沿罐壁面的分布。可见：除芦山53Elk-SN，其它工况基本一致，即峰值水动压从底部到自由液面逐步减小，且增长速率逐渐减缓，极大值总出现在底部。
![]()
图 11 峰值水动压沿罐壁的分布Figure 11. Distributions of peak hydrodynamic pressures along the tank sidewall图12为芦山53Elk-SN载荷下罐内水动压部分时刻的分布情况。可见：液体向右侧聚集后出现水跃，从43.2 s到43.4 s，液体裹挟空气抨击罐壁，自由液面附近出现峰值水动压（约为20.1 kPa）。其余工况以对流压力特性为主导。
![]()
图 12 芦山53Elk-SN地震载荷下罐内水动压分布Figure 12. Distributions of hydrodynamic pressures in the tank under Lushan 53Elk-SN excitation储罐上、中、下三个测点P6，P4和P2的水动压频谱特性如图13a和13b所示，可见：P6处出现多个谱峰，最大谱峰值出现在储液一阶固有频率，说明储罐上部水动压为低频的共振对流；P4处的最大谱峰值也出现在储液一阶固有频率，并且次要谱峰也出现在地震主要频率处，表明水动压为对流压力与脉冲压力的组合模式，但以对流压力为主导；P2处水动压的谱峰值与地震主要频率完全一致，表明其为脉冲冲击模式。
![]()
13 储罐内水动压的快速傅里叶分析结果13. Fast Fourier transform results of hydrodynamic pressures in the tank![]()
13 储罐内水动压的快速傅里叶分析结果13. Fast Fourier transform results of hydrodynamic pressures in the tankPGV，PGA与峰值水动压的关系如图14所示，可见：P4和P2处的峰值水动压与PGV的关系并不明显，但其与PGA呈明显的线性关系，且PGA较大的地震能激发较高的水动压；峰值水动压在下部（P2点）的增长速率明显大于中部（P4点），说明越靠近底部，水动压的脉冲特性越强，因此在储罐设计时，要充分考虑当地的场地条件（重点关注PGA较大的区域），并适当提高设计指标确保罐壁下部的强度；P6处水动压表现为共振对流模式，同属低频的埃尔津、芦山53Elk-SN及北岭-1地震具有相近的PGV （分别为0.84 m/s，0.866 m/s和0.8 m/s），但埃尔津的主频（0.563 Hz）远离储液一阶固有频率，而芦山53Elk-SN和北岭-1的主频（分别为0.325 Hz和0.353 Hz）更加接近储液一阶固有频率，因此，二者激发的共振晃荡产生的对流压力使得储罐上部（P6点）的水动压远大于非共振时的工况，其峰值分别为前者的2.51倍和3.35倍。上述讨论表明，地震主频是影响储罐上部水动压的主要因素，且其峰值随PGV的增大而上升。
![]()
图 14 各工况下PGV和PGA与峰值水动压的关系Figure 14. The relationships between PGV （or PGA） and peak hydrodynamic pressure for all cases4. 讨论与结论
本文基于计算流体力学软件Fluent，采用多参考系法对流体晃荡进行数值模拟研究，探究地震载荷下圆柱型储罐内储液的晃荡特性，得出以下结论：
1） 地震的频率成分可用于表征地震的主周期，与储罐模型共同决定晃荡波的共振模式，能反映出晃荡的共振状态，同时共振模式还受地震主频和PGV的影响。
2） 自由液面响应以对流为主，随着PGV的增加，波高整体增大，但存在特例，如：北岭-1与埃尔津地震（或神户与曼杜西诺角地震）的PGV相近，前者的主频低，诱发的波高是后者的1.4倍（或1.5倍）。最剧烈的工况为塔夫脱地震，其具有最大的PGV且主频最低。
3） 水动压表现为脉冲和对流两种模式。上部为对流模式，受地震主频和PGV影响；中部为对流和脉冲的组合模式；下部为脉冲模式，其峰值与PGA呈线性关系。该发现为储罐结构自上而下的整体振动响应设计提供了科学依据，有助于找到降低结构振动响应的工程方案。
4） 峰值水动压一般出现在罐壁底部，对于波面破碎的工况，峰值水动压则出现在自由液面附近。因此在储罐抗震设计时，应加强罐壁下部的强度，同时兼顾液面附近破碎冲击带来的强冲击。
本研究给出了极端工况的激发机制，即：地震属于低频且具有强PGV，以及水动力载荷峰值和频谱特性沿储罐壁面的演化规律。后续将针对此类工况下的结构振动分析和流体减晃研究，在液舱内布置刚性或柔性结构，以降低晃荡波高和水动力载荷。
-
![]()
图 1 密集台阵及震源位置分布图
Figure 1. Location of dense seismic array and distribution of hypocenters
![]()
图 2 研究区射线覆盖及信号叠加对比示意图
（a） 研究区内的射线分布；（b） 710号台站原始信号（上）和叠加后信号（下）
Figure 2. Distribution of rays in the study area and the comparison of signal before and after stacking
（a） Distribution of rays in the study area；（b） Original signal （upper） and stacked signal （lower） of the station No. 710
![]()
图 3 尾波Q值拟合示意图（以气枪信号为例）
（a） 11号台站叠加后信号，tc表示尾波起算时间，红色方框表示截取的尾波窗口；（b） 函数$\mathrm{l}\mathrm{n} ［ A （ f|r \text{，} t ） /K （ r \text{，} \alpha ） ］ $随时间的变化曲线，红色直线为最小二乘拟合线
Figure 3. Diagram of coda Q value fitting (take the airgun signal as an example)
（a） The stacked signal of the station No. 11，tc represents the the starting point of coda，the red box indicates the selected coda window；（b） The function $\mathrm{l}\mathrm{n} ［ A （ f|r \text{，} t ） /K （ r\text{，} \alpha ） ］ $ versus lapse time along with least squares fits of selected coda window，the red line represents the fitted line
![]()
图 4 不同中心频率下Qc值的空间分布及与速度层析成像结果对比
（a−e） 2—6 Hz频率尾波Qc值的空间分布，黑线表示宾川盆地范围；（f） 宾川地区0 km处速度层析成像结果（引自张云鹏等，2020）
Figure 4. Spatial distribution of Qc values at different frequencies and comparison with velocity tomography
（a−e） Spatial distribution of coda Qc values at 2−6 Hz frequency，the black line indicates the extent of Binchuan basin；（f） Velocity tomography results at 0 km depth in Binchuan area （from Zhang et al，2020）
![]()
图 5 六个台站位置（a）及其Qc值与频率的关系（b）
Figure 5. Location of the six stations （a） and the relationship between Qc and frequency for the stations （b）
![]()
图 6 30号台站Q0和频率依赖指数 $ \eta $ 随流逝时间的变化曲线
Figure 6. Values of Q0 and frequency-dependent coefficient $ \eta $ as function of lapse time for the station No. 30
![]()
图 7 全区域各台站尾波衰减参数及标准差
（a） 各台站Q0统计特征；（b） 各台站频率依赖指数 $ \eta $ 统计特征
Figure 7. Coda attenuation parameters and their standard deviations of each station in the study area
（a） The statistical characteristics of Q0 of each station；（b） The statistical characteristics of the frequency dependence coefficient $ \eta $ of each station
表 1 所用地震详细参数
Table 1 Parameters of the earthquakes using in the study
发震时刻 东经/° 北纬/° ML 震源
深度/km发震时刻 东经/° 北纬/° ML 震源
深度/km年-月-日 时:分:秒 年-月-日 时:分:秒 2017-03-26 16:24:03 100.563 25.852 0.4 17 2017-04-28 00:16:55 100.445 25.831 1.3 15 2017-03-27 03:28:14 100.394 25.928 0.8 8 2017-05-07 10:20:12 100.698 25.891 1.6 5 2017-03-31 04:54:56 100.594 25.829 0.4 18 2017-05-12 11:27:14 100.681 25.863 1.3 6 2017-04-03 17:48:41 100.409 25.815 1.0 14 2017-05-12 11:36:35 100.648 25.884 1.4 10 2017-04-12 21:05:22 100.673 25.659 0.5 9 2017-05-21 21:42:29 100.451 25.771 1.0 5 2017-04-13 10:41:05 100.613 25.908 1.7 5 2017-05-23 09:19:53 100.359 25.680 1.5 12 
下载: 导出CSV
表 2 30号台站不同流逝时间下Q0和频率依赖指数 $ \eta $ 的统计结果
Table 2 Q0 and frequency-dependent coefficient $ \mathrm{\eta } $ at different lapse times of the station No.30
流逝时间/s Q0 ${ \delta }{ Q_0}$ $ \eta $ ${\delta }\eta$ 5 26.00 15.55 1.057 0.342 10 38.11 18.59 0.988 0.242 15 50.60 24.27 0.937 0.227
下载: 导出CSV
表 3 Q0标准差较大的台站详细信息
Table 3 Information of the stations with larger standard deviation of Q0
台站编号 东经/° 北纬/° Q0 $\delta{ Q_0}$ $ \eta $ ${\delta }\eta$ 195 100.399 1 25.719 9 38 37 1.123 0.431 633 100.523 5 25.657 2 48 36 0.713 0.674 613 100.538 8 25.680 0 61 33 0.702 0.401 83 100.670 6 25.835 4 40 31 0.958 0.298 18 100.630 4 25.912 7 43 29 0.831 0.282
下载: 导出CSV
表 4 表3中相邻台站的Q0结果
Table 4 Results of Q0 for neighboring stations in Table 3
台站编号 东经/° 北纬/° Q0 $\delta{ Q_0}$ $ \eta $ ${\delta }\eta$ 196 100.417 5 25.722 4 25 5 1.148 0.198 632 100.501 0 25.660 6 34 5 0.912 0.218 752 100.530 1 25.680 8 18 6 1.195 0.202 81 100.639 7 25.840 2 31 16 0.857 0.347 17 100.620 1 25.901 3 23 11 1.107 0.335
下载: 导出CSV
-
陈婷,董建辉,王晓山. 2017. 唐山地区尾波 Q值时空分布特征[J]. 地震地磁观测与研究,38(2):51–56. Chen T,Dong J H,Wang X S. 2017. Temporal and spatial distribution characteristics of the coda wave Q value in Tangshan region[J]. Seismological and Geomagnetic Observation and Research,38(2):51–56 (in Chinese).
洪玉清,杨选. 2015. 利用Aki模型对新丰江水库地区尾波 Q值的研究[J]. 华南地震,35(3):66–71. doi: 10.13512/j.hndz.2015.03.010 Hong Y Q,Yang X. 2015. Q value research of coda wave in Xinfengjiang reservoir area based on Aki model[J]. South China Journal of Seismology,35(3):66–71 (in Chinese).
刘斌,Kern H,Popp T. 1998. 不同围压下孔隙度不同的干燥及水饱和岩样中的纵横波速度及衰减[J]. 地球物理学报,41(4):537–546. doi: 10.3321/j.issn:0001-5733.1998.04.012 Liu B,Kern H,Popp T. 1998. Velocities and attenuation of P- and S-waves in dry and wet rocks with different porosities under different confining pressures[J]. Acta Geophysica Sinica,41(4):537–546 (in Chinese).
罗睿洁,吴中海,黄小龙,黄小巾,周春景,田婷婷. 2015. 滇西北宾川地区主要活动断裂及其活动构造体系[J]. 地质通报,34(1):155–170. doi: 10.3969/j.issn.1671-2552.2015.01.013 Luo R J,Wu Z H,Huang X L,Huang X J,Zhou C J,Tian T T. 2015. The main active faults and the active tectonic system of Binchuan area,northwestern Yunnan[J]. Geological Bulletin of China,34(1):155–170 (in Chinese).
师海阔,曾宪伟,张立恒,贺永忠. 2016. 利用Aki模型对宁夏及邻区尾波 Q值分布特征的研究[J]. 地震工程学报,38(1):51–57. doi: 10.3969/j.issn.1000-0844.2016.01.0051 Shi H K,Zeng X W,Zhang L H,He Y Z. 2016. Distribution characteristics of Q values of seismic coda in Ningxia and neighboring area based on Aki model[J]. China Earthquake Engineering Journal,38(1):51–57 (in Chinese).
史水平,周斌,黄树生,阎春恒,郭培兰. 2020. 广西龙滩水库库区地震尾波衰减特征[J]. 地震学报,42(2):151–162. Shi S P,Zhou B,Huang S S,Yan C H,Guo P L. 2020. Characteristics of seismic coda attenuation in Longtan reservoir of Guangxi region[J]. Acta Seismologica Sinica,42(2):151–162 (in Chinese).
苏金波,王宝善,王海涛,王琼,冀战波. 2015. 利用大容量气枪震源资料研究北天山地区介质衰减特征[J]. 地震研究,38(4):598–605. Su J B,Wang B S,Wang H T,Wang Q,Ji Z B. 2015. Research on characteristic of seismic attenuation in the northern Tianshan area using seismic signal from airgun source[J]. Journal of Seismological Research,38(4):598–605 (in Chinese).
苏有锦. 2009. 云南地区地震波衰减(Q值)结构反演成像研究[D]. 合肥:中国科学技术大学:36−38. Su Y J. 2009. Inversion Tomography of the Seismic Wave Attenuation (Q Value) Structure in Yunnan Region[D]. Hefei:University of Science and Technology of China:36−38 (in Chinese).
孙业君,黄耘,王斌,王俊,李锋,江昊琳. 2014. 江苏地区尾波 Q值特征研究[J]. 地震,34(1):24–33. Sun Y J,Huang Y,Wang B,Wang J,Li F,Jiang H L. 2014. Characteristics of coda wave Q values in Jiangsu area[J]. Earthquake,34(1):24–33 (in Chinese).
王宝善,王伟涛,葛洪魁,徐平,王彬. 2011. 人工震源地下介质变化动态监测[J]. 地球科学进展,26(3):249–256. Wang B S,Wang W T,Ge H K,Xu P,Wang B. 2011. Monitoring subsurface changes with active sources[J]. Advances in Earth Science,26(3):249–256 (in Chinese).
王勤彩,陈章立,王中平,郑斯华. 2010. 云南地区散射衰减、吸收衰减及尾波衰减的综合研究[J]. 地球物理学进展,25(2):419–431. doi: 10.3969/j.issn.1004-2903.2010.02.007 Wang Q C,Chen Z L,Wang Z P,Zheng S H. 2010. Comprehensive study of scattering,intrinsic and coda attenuation in the Yunnan area[J]. Progress in Geophysics,25(2):419–431 (in Chinese).
王振宇,赵培培,薄景山. 2017. 地震动随机模拟方法主要影响参数分析[J]. 世界地震工程,33(3):34–41. Wang Z Y,Zhao P P,Bo J S. 2017. Analysis for the effects of main parameters on ground motions by stochastic simulation method[J]. World Earthquake Engineering,33(3):34–41 (in Chinese).
杨微,王宝善,葛洪魁,王伟涛,陈颙. 2013. 大容量气枪震源主动探测技术系统及试验研究[J]. 中国地震,29(4):399–410. Yang W,Wang B S,Ge H K,Wang W T,Chen Y. 2013. The active monitoring system with large volume airgun source and experiment[J]. Earthquake Research in China,29(4):399–410 (in Chinese).
张云鹏,王宝善,林国庆,王伟涛,杨微,吴中海. 2020. 利用密集台阵近震层析成像研究云南宾川上地壳速度结构[J]. 地球物理学报,63(9):3292–3306. Zhang Y P,Wang B S,Lin G Q,Wang W T,Yang W,Wu Z H. 2020. Upper crustal velocity structure of Binchuan,Yunnan revealed by dense array local seismic tomography[J]. Chinese Journal of Geophysics,63(9):3292–3306 (in Chinese).
赵秋芳,云美厚,朱丽波,李晓斌,李伟娜. 2019. 近地表 Q值测试方法研究进展与展望[J]. 石油地球物理勘探,54(6):1397–1418. doi: 10.13810/j.cnki.issn.1000-7210.2019.06.026 Zhao Q F,Yun M H,Zhu L B,Li X B,Li W N. 2019. Progress and outlook of near-surface quality factor Q measurement and inversion[J]. Oil Geophysical Prospecting,54(6):1397–1418 (in Chinese).
周连庆,赵翠萍,修济刚,陈章立,郑斯华. 2008. 利用天然地震研究地壳 Q值的方法和进展[J]. 国际地震动态,(2):1–11. Zhou L Q,Zhao C P,Xiu J G,Chen Z L,Zheng S H. 2008. Methods and developments of research on crustal Q value by using earthquakes[J]. Recent Developments in World Seismology,(2):1–11 (in Chinese).
周庆,虢顺民,向宏发. 2004. 滇西北地区潜在震源区的划分原则和方法[J]. 地震地质,26(4):761–771. doi: 10.3969/j.issn.0253-4967.2004.04.022 Zhou Q,Guo S M,Xiang H F. 2004. Principle and method of delineation of potential seismic sources in northeastern Yunnan Province[J]. Seismology and Geology,26(4):761–771 (in Chinese).
朱新运. 2011. 衰减、场地响应等地震波传播相关信息综合研究[D]. 北京:中国地震局地球物理研究所:66−79. Zhu X Y. 2011. Comprehensive Study on Seismic Wave Propagation Related Information Such as Attenuation and Site Response[D]. Beijing:Institute of Geophysics,China Earthquake Administration:66−79 (in Chinese).
朱新运,刘杰,张帆. 2006. 基于Aki模型的近震S波尾波 Q值求解及分析软件研制[J]. 地震研究,29(1):76–80. doi: 10.3969/j.issn.1000-0666.2006.01.015 Zhu X Y,Liu J,Zhang F. 2006. Development of Q-value calculating and processing software using S-wave coda in local earthquakes based on Aki model[J]. Journal of Seismological Research,29(1):76–80 (in Chinese).
Abdel-Fattah A K,Morsy M,El-Hady S,Kim K Y,Sami M. 2008. Intrinsic and scattering attenuation in the crust of the Abu Dabbab area in the eastern desert of Egypt[J]. Phys Earth Planet Inter,168(1/2):103–112.
Aki K. 1969. Analysis of the seismic coda of local earthquakes as scattered waves[J]. J Geophys Res,74(2):615–631. doi: 10.1029/JB074i002p00615
Aki K,Chouet B. 1975. Origin of coda waves:Source,attenuation,and scattering effects[J]. J Geophys Res,80(23):3322–3342. doi: 10.1029/JB080i023p03322
Aki K. 1980. Scattering and attenuation of shear waves in the lithosphere[J]. J Geophys Res: Solid Earth,86(B11):6496–6504.
Correig A M,Mitchell B J,Ortiz R. 1990. Seismicity and coda Q values in the eastern Pyrenees:First results from the La Cerdanya seismic network[J]. Pure Appl Geophys,132(1):311–329.
Das R,Mukhopadhyay S,Singh R K,Baidya P R. 2018. Lapse time and frequency-dependent coda wave attenuation for Delhi and its surrounding regions[J]. Tectonophysics, 738 − 739 :51−63.
Dinesh K,Sarkar I,Sriram V,Khattri K N. 2005. Estimation of the source parameters of the Himalaya earthquake of October 19,1991,average effective shear wave attenuation parameter and local site effects from accelerograms[J]. Tectonophysics,407:1–24
Domínguez T,Rebollar C J,Fabriol H. 1997. Attenuation of coda waves at the Cerro Prieto geothermal field,Baja California,Mexico[J]. Bull Seismol Soc Am,87(5):1368–1374. doi: 10.1785/BSSA0870051368
Herrmann R B. 1980. Q estimates using the coda of local earthquake[J]. Bull Seismol Soc Am,70(2):447–468. doi: 10.1785/BSSA0700020447
Jin A S,Aki K. 1988. Spatial and temporal correlation between coda Q and seismicity in China[J]. Bull Seismol Soc Am,78(2):741–769. doi: 10.1785/BSSA0780020741
Kosuga M. 1992. Dependence of coda Q on frequency and lapse time in the western Nagano region,central Japan[J]. J Phys Earth,40(2):421–445. doi: 10.4294/jpe1952.40.421
Ma’hood M,Hamzehloo H. 2009. Estimation of coda wave attenuation in east central Iran[J]. J Seismol,13(1):125–139. doi: 10.1007/s10950-008-9130-2
Padhy S,Subhadra N,Kayal J R. 2011. Frequency-dependent attenuation of body and coda waves in the Andaman Sea Basin[J]. Bull Seismol Soc Am,101(1):109–125. doi: 10.1785/0120100032
Pulli J J. 1984. Attenuation of coda waves in New England[J]. Bull Seismol Soc Am,74(4):1149–1166.
Rautian T G,Khalturin V I. 1978. The use of the coda for determination of the earthquake source spectrum[J]. Bull Seismol Soc Am,68(4):923–948. doi: 10.1785/BSSA0680040923
Roecker S W,Tucker B,King J,Hatzfeld D. 1982. Estimates of Q in central Asia as a function of frequency and depth using the coda of locally recorded earthquakes[J]. Bull Seismol Soc Am,72(1):129–149. doi: 10.1785/BSSA0720010129
Sato H. 1977. Energy propagation including scattering effects single isotropic scattering approximation[J]. J Phys Earth,25(1):27–41. doi: 10.4294/jpe1952.25.27
Su F,Aki K. 1995. Site amplification factors in central and southern California determined from coda waves[J]. Bull Seismol Soc Am,85(2):452–466.
Tian W,Li Z W. 2020. S wave anelastic attenuation of shallow sediments in mainland China[J]. Earth Space Sci,7(10):e2020EA001348. doi: 10.1029/2020EA001348
Toksöz M N,Johnston D H,Timur A. 1979. Attenuation of seismic waves in dry and saturated rocks:I. Laboratory measurements[J]. Geophysics,44(4):681–690. doi: 10.1190/1.1440969
Wang S,Li Z W. 2018. S-wave attenuation of the shallow sediments in the North China basin based on borehole seismograms of local earthquakes[J]. Geophys J Int, 214 :1391−1400.
Wang W,Shearer P M. 2019. An improved method to determine coda- Q,earthquake magnitude,and site amplification:Theory and application to southern California[J]. J Geophys Res: Solid Earth,124(1):578–598. doi: 10.1029/2018JB015961
Wong V,Rebollar C J,Munguía L. 2001. Attenuation of coda waves at the Tres Vírgenes Volcanic area,Baja California Sur,Mexico[J]. Bull Seismol Soc Am,91(4):683–693. doi: 10.1785/0120000025
计量
- 文章访问数: 133
- HTML全文浏览量: 52
- PDF下载量: 29
