Seismic response analysis method of offshore site based on the seismic wave input technique of hybrid wave field
-
摘要: 基于人工边界子结构模型，提出一种利用混合波场实现近海场地中地震P波和SV波垂直输入的方法。该方法中用于波动输入的混合波场由计算模型两侧截断边界的自由波场和底面边界的入射波场构成，避免了不规则近海场地的自由波场求解。同时采用基于声流体单元的流固耦合算法模拟场地-海水动力相互作用，利用流体介质人工边界和黏弹性人工边界单元模拟无限流、固介质的辐射阻尼，建立了近海场地动力相互作用分析模型。将该模型与基于混合波场的地震波动输入方法相结合，提出了近海场地地震反应时域整体分析方法。数值算例表明，本文的计算模型和分析方法可以较为准确地实现地震波在近海场地中的输入问题，对于近海场地及海洋工程的地震反应分析具有良好的适用性。Abstract: The seismic response analysis of offshore site is the foundation for the seismic safety evaluation of ocean engineering. Seismic wave input is a key part of site seismic response analysis. However, the slope of the offshore site inclines to the ocean, which interacts with the seawater directly, leading to natural difficulties in solving the free wave field and thus provides a challenge to the existing seismic wave input methods based on the free wave fields. In this study, we adopted and improved the seismic wave input method based on the substructure of artificial boundaries, and constructed a hybrid input wave field with the free wave field on two lateral boundaries and the incident wave field on the bottom boundary of the calculation model. On that basis, an offshore site seismic wave input method for vertically incident P and SV waves is proposed. The proposed method avoids calculating the free wave field of irregular offshore sites, therefore it can simplify the seismic response analysis of offshore engineering considering characteristics of offshore site and the soil-structure interaction. Meanwhile, the fluid-structure coupling algorithm based on the acoustic fluid element is used to simulate the dynamic interaction between the offshore site and the seawater, and the fluid dynamic artificial boundary and the uniform viscoelastic artificial boundary element are used to simulate the radiation damping of the infinite fluid and solid medium. On that basis, the dynamic interaction model of offshore site is established. Combine this model with the proposed seismic wave input method based on the hybrid wave field, thereby an overall time-domain analysis method for the seismic response of offshore sites is established. Numerical examples show that the analysis method and numerical model proposed in this study can accurately realize the input of seismic waves in offshore sites, and have good applicability for seismic response analysis of offshore engineering.
-
引言
水氡观测是最为广泛的前兆测项之一（刘耀炜等，2015）。连续观测地下氡含量变化，可捕捉到地震孕育与发生的信息（姚玉霞等，2014；张磊等，2016）。氡的常规测量通常采用静电计电离法、硫化锌闪烁室法等，我国地震监测台网常用的测氡仪器有：FD-125型室内氡钍分析仪、FD-105K型静电计测氡仪器，JDZ-1型自动测氡仪、SD-1型双道自动测氡仪等。国内自主生产的闪烁室法测氡仪在地震系统中长期使用，为地震监测预报提供了有用信息（周红艳等，2015），并积累了大量的观测资料。然而，闪烁室法测氡仪自身存在难以克服的缺点，如：观测环节繁琐、氡长期积累导致降本底耗时、固体源活度失准、仪器校准不准确、计算过程复杂等，这些缺陷已经严重影响到观测资料的内在质量（任宏微等，2016）。
近年来，随着电子技术和计算机的飞速发展，用于粒子探测的核仪器性能也得到迅速提升，国内外成立了多家专业仪器生产厂（杨明太，2011）。诸多学者对新型仪器进行了测试研究，期望通过引入新仪器来代替传统测氡仪，以提供高质量的地震观测数据。DDL-1气氡仪是郑州晶微科技有限公司生产的自动化、网络化、数字化的测氡仪，其原理是利用电离法对氡气含量进行检测。该仪器性能稳定、一致性好，能够正确反映地下水中氡含量的真实变化（起卫罗等，2017，2019）。
从1984年至今，嘉峪关台一直使用FD-125闪烁室法测氡仪观测断层气氡，曾成功预报了2002年12月14日玉门MS5.9地震，2003年10月25日山丹、民乐MS6.1和MS5.8地震，为地震监测提供了宝贵资料。但随着仪器老化，观测数据波动增大，其映震效能逐渐下降。鉴于嘉峪关断层气氡观测点的实际需求，本文引进DDL-1气氡仪，对仪器性能、仪器校准、观测数据稳定性及连续观测等方面进行试验研究，对其替代传统测氡仪的可行性进行讨论。
1. 测氡仪工作原理
1） FD-125氡钍分析仪基于闪烁法原理进行测量。当氡气进入闪烁室，其衰变过程中释放的α粒子与ZnS （Ag）晶体发生碰撞，激发ZnS （Ag）原子释放出光子。这些光子被光电倍增管捕获，进而产生光电子，完成光电转换过程。闪烁室内α粒子的数量与氡气浓度之间存在正比关系，即氡气浓度与闪光频率成正比。因此，通过记录光电倍增管输出的脉冲频率，可以计算出闪烁室内的氡浓度（姚玉霞等，2017）。
2） DDL-1测氡仪基于电离法原理进行测量。利用氡气的放射性能电离周围介质的特性。当氡气进入气氡仪的电离室，由氡衰变产生的α粒子和RaA气子体引起电流的累积。氡气持续从进气口流入电离室，在特定正电场的作用下，产生的杂散离子形成向中心流动的电离离子流。这些定向的离子流在电离室中心的接收极上累积，其数量与电离室内氡气的浓度成正比。浓度越高，接收极上累积的离子流就越多，而电离室内的氡气浓度与观测井或断层中的氡气浓度直接相关。累积的离子流经过放大和模数转换后，输出至主机进行测值显示、存储和传输，实现氡气含量检测的目的。
测氡仪的主要技术指标和参数能反映该仪器的性能质量，FD-125氡钍分析仪和DDL-1气氡仪的主要技术参数对比列于表1，两套仪器出厂时的各项技术指标均执行中国地震局地震水文地球化学观测技术规范要求（中国地震局，2014）。由表可以看出，相较于FD-125氡钍分析仪，DDL-1气氡仪的采样率和灵敏度均更高。
表 1 FD-125氡钍分析仪和DDL-1气氡仪主要技术指标及参数Table 1. Main technical specifications and parameters of FD-125 radon-thoron analyzer and DDL-1 gas radon meter仪器型号 检测对象 灵敏度 稳定度 采样率 环境温度/℃ 环境湿度 电源 FD-125氡钍
分析仪气体样品 ≥2.0 cpm/pCi/L ≥90% 1次/10分钟 0—45 ≤80% AC 200 V—240 V （50 Hz） DDL-1气氡仪 气体样品 0.1 Bq/L ≥90% 1次/分钟 0—40 ＜80% AC 200 V—240 V （50 Hz）；
DC 10 V—15 V2. 监测点概况
监测点位于嘉峪关市西北方向6 km处的嘉峪关断层中段。嘉峪关断层总体走向为北偏西35°—40°，倾向西南，倾角介于73°—85°之间，属右旋走滑逆冲型断层。该监测点是嘉峪关断层气固定监测点，自1987年8月起被选定并投入运行，每日在采样点进行样本采集，并将样本带回观测室测量氡气浓度。采样井的建造方法如下：先向下挖掘深1.2 m，直径10 cm的孔洞，随后使用炸药爆破形成一个深约1.8 m、直径30 cm的柱形坑作为气体收集坑。将直径为30 cm的玻璃漏斗倒扣在集气口上，漏斗颈部布设两个通气导管，其中一根穿过漏斗颈延伸至坑底，用于深孔观测，另一根导管则刚好穿过漏斗颈，用于浅孔观测。结构示意如图1所示。此外，导气管外部包裹铜质套管以提供保护，并进行密封处理。导气管直径为2—3 mm，通过转换接头与乳胶管相连后引出地面。整个系统与大气相通，深孔观测深度为3 m，浅孔观测深度为1.2 m。
自1987年至今，嘉峪关监测点使用FD-125氡钍分析仪观测的气氡浓度年变动态表现为典型的夏高冬低型。通常情况下，1月或2月氡浓度位于谷值，少数情况下谷值会出现在3月，峰值则通常出现在7月或8月。氡浓度曲线的趋势性转折和突跳异常对祁连山地区M＞6.0地震有较高的预报效能，可作为衡量该地区地震活动性的重要指标。
3. 实验过程与分析
3.1 标定校准
根据氡观测技术规范要求（中国地震局，2014），观测仪器在使用中必须定期校准。本实验按照规范要求在同一时间用FD-
3024 固体源对DDL-1和FD-125进行校准，得到FD-125氡钍分析仪的K值为0.008 43 Bq/（脉冲·min−1），DDL-1气氡仪的K值为0.445 3 Bq/L，校准结果列于表2。对比整个校准过程，两套仪器均需要连续标定3次，取平均值得到仪器K值。但FD-125标定过程较繁琐，闪烁室降本底耗时，需要三天才能完成标定。而DDL-1过程简单且费时较少，可连续标定，自动计算平均K值，只需5小时就可以完成标定。本次校准结果，满足本次实验的要求。表 2 FD-125和DDL-1的 校准结果Table 2. Calibration results of the FD-125 and DDL-1仪器名称 标定次数 温度
/℃湿度 气压
/hPa本底
/（Bq·L−1）测值
/（Bq·L−1）各次校准
K值新K值 各次校准
相对误差新-旧K值
相对误差FD-125
氡钍分析仪第一次 19.0 37% 829.7 3.804 9 21 952 0.008 44 0.008 43 0.1% 1.3% 第二次 18.0 38% 829.1 4.227 7 22 353 0.008 30 −1.5% 第三次 18.0 38% 829.1 5.073 2 21 683 0.008 56 1.5% DDL-1
气氡仪第一次 19.5 37% 826.8 2.437 7 40.527 2 0.455 5 0.445 3 2.3% 1.5% 第二次 19.5 37% 826.8 2.526 7 40.743 4 0.453 1 2.3% 第三次 19.5 37% 826.8 3.945 9 43.185 4 0.427 5 4.0% 注：表中 FD-125 氡钍分析仪的K值单位为Bq/（脉冲·min−1），DDL-1 气氡仪的K值单位为Bq/L. 3.2 降本底对比
由于FD-125氡钍分析仪已使用多年，闪烁室内氡长期积累，降本底较困难。FD-125氡钍分析仪2021年全年的本底曲线如图2a所示，由图可以看出，本底曲线呈冬高夏低的形态，夏季气温较高时降本底费时较少，冬季气温低时降本底耗时多，且很难降到4 Bq/L以下。DDL-1气氡仪本底曲线如图2b所示，由图可以看出，DDL-1气氡仪在气温较低时降本底容易，甚至可以降到0 Bq/L，但在夏季温度较高时降本底困难。
![]()
图 2 FD-125氡钍分析仪（a）和DDL-1气氡仪 （b）的本底曲线图Figure 2. Background value curves of FD-125 radon-thoron analyzer （a） and DDL-1 gas radon meter （b）3.3 数据可靠性实验
嘉峪关断层气氡日常观测，是在取样点用扩散瓶取样后带回观测室，按照定时、定点、定量的方式进行观测。2021年5月11日至9月30日期间，对DDL-1气氡仪和FD-125氡钍分析仪采用上述相同进样方式，静置一小时读取测值，共积累143天的数据。
两套仪器的观测值如图3所示，由图可见两套仪器测值同步性差，呈现出负相关，相关系数为−0.507 4。其原因可能是DDL-1气氡仪的传感器体积较大，为0.7 L，进气、排气管道较长，而气样较少，气样在传感器及管道内分布不均，以致测值误差偏大。
![]()
图 3 采用相同进样方式时DDL-1气氡仪与FD-125氡钍分析仪的观测结果对比图Figure 3. Comparison of the observation results of DDL-1 gas radon meter and FD-125 radon-thoron analyzer with the same sample injection method为减小误差，从10月1日开始改变DDL-1气氡仪进样方式，从取样点取回两扩散瓶气样，一瓶用于FD-125观测，另一瓶通过DDL-1配套的循环泵进行循环进样观测，用止血钳夹住传感器的进气口，循环泵抽气使传感器和橡胶管道呈负压状态，然后把装有气样的扩散瓶接入循环装置，取开止血钳，形成闭合的循环回路，连接方式如图4所示。连续观测一小时后读取测值，采用循环进样方法后，共累积了71天数据。将DDL-1气氡仪与FD-125氡钍分析仪测值进行对比，结果如图5所示。由图可见，改为循环进样后，DDL-1气氡仪的观测数据质量有明显提高，两套仪器观测值同步性较好，相关系数为0.868，测值变化较平缓。但FD-125的部分转折变化尚未同步，表明该部分变化可能为人工进样引起的误差。通过对比表明，DDL-1气氡仪的观测资料更为可靠，且仪器稳定性高。
![]()
图 4 循环进样观测装置示意图Figure 4. Schematic diagram of the cyclic injection observation device![]()
图 5 采用循环进样方法DDL-1气氡仪与FD-125氡钍分析仪观测结果对比图Figure 5. Comparison of the observation results of DDL-1 gas radon meter using cyclic sample injection method and FD-125 radon-thoron analyzer3.4 连续观测试验
将DDL-1气氡仪架设到取样点，依次以取气口—干燥管—传感器—循环泵进行连接，循环泵排气口接通外界空气。2021年12月11日9时开始观测，初期测值较平稳约22 Bq/L，19时30分左右，测值持续下降至1 Bq/L左右，经检查地下井内导气管（直径2—3 mm）内结冰，堵塞导气管。13日取样后重新接通观测，晚间再次结冰。由于当地冬季气温最低可达−20 ℃，观测点在野外，无法采用保温措施使温度达0 ℃以上。为此，我们等待气温回升后再次进行试验。于2022年3月30日再次在取样点进行，试验结果如图6所示，可以看出仪器架设初期测值较高，随后开始下降，至28 Bq/L后趋于平稳，但数小时后测值持续下降。同时FD-125的测值由50 Bq/L下降到10 Bq/L左右，浅孔测值同步降低，严重影响到正常观测，试验停止。考虑到在降低采样率的情况下，仪器配套的真空泵仍然持续抽气，无法控制定时开关，因此没有再降低采样率进行试验。
![]()
图 6 3月30日—4月1日嘉峪关监测点利用DDL-1气氡仪进行连续观测曲线Figure 6. Continuous observation curve at Jiayuguan monitoring point using DDL-1 gas radon meter （March 30th to April 1st）试验结果表明：使用循环泵连续抽气观测，导致地下井内气体密度减小，并从周围断层裂隙中抽取气体，在地下结构无变化的情况下，从裂隙中排出的气量较小，导致测值持续降低；该测点深孔和浅孔观测在同一个井中，由导气管的深度不同区分深孔和浅孔，深孔在连续抽气时会导致测值降低，同时致使浅孔测值也同步降低。
4. 讨论与结论
本文对 DDL-1 气氡仪和 FD-125 测氡仪在仪器工作原理、主要技术参数、校准流程及观测数据稳定性等方面深入对比，探讨其性能表现及适用性。
从仪器设计来看，DDL-1气氡仪各项技术指标和参数满足地下流体氡测量的基本要求，其优势明显，采样率较高，可实现连续观测，且支持远程传输与参数控制，为数字化观测提供便利。在仪器校准方面，操作简便，能迅速完成校准并自动计算K值，极大地提高了工作效率；稳定性方面，该仪器稳定性高，能够体现数据的真实变化。但在降本底方面存在严重缺陷，当气温升高时，本底降低的难度增大，随着采样率的提升，仪器的灵敏度下降，影响测量的精准度。
连续测量时，测值持续降低，分析原因，可能是DDL-1 气氡仪连续抽气会改变地下井内的气体密度，从周围断层裂隙抽取气体，在地下结构稳定的情况下，断层气排出量逐渐减少，致使测量值持续降低。鉴于此，建议厂家对 DDL-1 气氡仪进行升级优化，为其配备具有定时切换进气源功能的智能设备，在低采样率测量时，通过抽取空气冲洗传感器来降低本底，以保障数据的可靠性与稳定性，从而更好地满足不同观测场景的需求。
综上，认为DDL-1 气氡仪设计具有先进性，在校准与稳定性方面表现良好，但也存在降本底困难、高采样率下灵敏度会降低，并且在特定环境下连续观测受限等问题（冬季导气管易结冰，无法进行连续观测）。因此，该仪器还不适用于嘉峪关断层气氡的连续观测。
-
![]()
图 2 流体介质动力人工边界示意图
Figure 2. Schematic diagram of dynamic artificial boundaries of fluid medium
![]()
图 3 基于混合波场的近海场地地震动输入方法
Figure 3. Seismic wave input method of offshore site based on hybrid wave field
![]()
图 4 基于混合波场的近海场地地震波动输入方法实现步骤
（a） 利用子结构模型求解等效地震荷载；（b） 在近海场地模型中施加等效地震荷载
Figure 4. Implementation steps of seismic wave input method of offshore site based on hybrid wave field
（a） Solving equivalent seismic loads through the substructure model；（b） Applying equivalent seismic loads on the offshore site model
![]()
图 7 SV波入射下近海场地的水平（a）和垂直（b）位移波形
Figure 7. Displacement waveforms in horizontal （a） and vertical （b） directions on offshore site under incident SV wave
![]()
图 8 SV波入射下近海场地位移峰值放大系数的空间分布
Figure 8. Spatial distributions of peak displacement amplification coefficient on offshore site under incident SV wave
![]()
图 9 P波入射下近海场地的水平（a）和垂直（b）位移波形
Figure 9. Displacement waveforms in horizontal （a） and vertical （b） directions on offshore site under incident P wave
![]()
图 10 P波入射下近海场地位移峰值放大系数的空间分布
Figure 10. Spatial distributions of peak displacement amplification coefficient on offshore site under incident P wave
表 1 介质材料参数
Table 1 Material parameters of media
流体介质 固体介质 密度ρF
/（kg·m−3）声速cF
/（m·s−1）密度ρ
/（kg·m−3）剪切波波速cS
/（m·s−1）泊松比γ 1000 1435.27 2 000 200 0.25 
下载: 导出CSV
-
宝鑫,刘晶波. 2017. 考虑流-固耦合效应的含液容器动力响应有限元分析方法[J]. 核动力工程,38(2):111–114. Bao X,Liu J B. 2017. Dynamic finite element analysis methods for liquid container considering fluid-structure interaction[J]. Nuclear Power Engineering,38(2):111–114 (in Chinese).
陈苏,周越,李小军,傅磊. 2018. 近海域地震动的时频特征与工程特性[J]. 振动与冲击,37(16):227–233. Chen S,Zhou Y,Li X J,Fu L. 2018. Time-frequency and engineering characteristics on offshore ground motion[J]. Journal of Vibration and Shock,37(16):227–233 (in Chinese).
陈少林,柯小飞,张洪翔. 2019. 海洋地震工程流固耦合问题统一计算框架[J]. 力学学报,51(2):594–606. doi: 10.6052/0459-1879-18-333 Chen S L,Ke X F,Zhang H X. 2019. A unified computational framework for fluid-solid coupling in marine earthquake engineering[J]. Chinese Journal of Theoretical and Applied Mechanics,51(2):594–606 (in Chinese).
杜修力. 2009. 工程波动理论与方法[M]. 北京: 科学出版社: 215–216. Du X L. 2009. Theories and Methods of Wave Motion for Engineering[M]. Beijing: Science Press: 215–216 （in Chinese）.
冯启民. 1990. 海洋工程场址地震动的分析方法[J]. 地震工程与工程振动,10(1):81–88. Feng Q M. 1990. Analytical method of seismic motion at offshore engineering site[J]. Earthquake Engineering and Engineering Vibration,10(1):81–88 (in Chinese).
胡进军,郑旭,郝彦春,谢礼立. 2017. 俯冲带地震动特征及其衰减规律探讨[J]. 地球物理学报,60(5):1773–1787. doi: 10.6038/cjg20170514 Hu J J,Zheng X,Hao Y C,Xie L L. 2017. Characterization of strong motion of subduction earthquakes and its attenuation relationship[J]. Chinese Journal of Geophysics,60(5):1773–1787 (in Chinese).
李小军. 2006. 海域工程场地地震安全性评价的特殊问题[J]. 震灾防御技术,1(2):97–104. doi: 10.3969/j.issn.1673-5722.2006.02.002 Li X J. 2006. Special problems on evaluation of seismic safety for offshore engineering site[J]. Technology for Earthquake Disaster Prevention,1(2):97–104 (in Chinese).
刘晶波,吕彦东. 1998. 结构-地基动力相互作用问题分析的一种直接方法[J]. 土木工程学报,31(3):55–64. Liu J B,Lü Y D. 1998. A direct method for analysis of dynamic soil-structure interaction[J]. China Civil Engineering Journal,31(3):55–64 (in Chinese).
刘晶波,谷音,杜义欣. 2006. 一致粘弹性人工边界及粘弹性边界单元[J]. 岩土工程学报,28(9):1070–1075. doi: 10.3321/j.issn:1000-4548.2006.09.004 Liu J B,Gu Y,Du Y X. 2006. Consistent viscous-spring artificial boundaries and viscous-spring boundary elements[J]. Chinese Journal of Geotechnical Engineering,28(9):1070–1075 (in Chinese).
刘晶波,宝鑫,谭辉,王建平,郭东. 2017. 波动问题中流体介质的动力人工边界[J]. 力学学报,49(6):1418–1427. doi: 10.6052/0459-1879-17-199 Liu J B,Bao X,Tan H,Wang J P,Guo D. 2017. Dynamical artificial boundary for fluid medium in wave motion problems[J]. Chinese Journal of Theoretical and Applied Mechanics,49(6):1418–1427 (in Chinese).
刘晶波,谭辉,宝鑫,王东洋,李述涛. 2018. 土-结构动力相互作用分析中基于人工边界子结构的地震波动输入方法[J]. 力学学报,50(1):32–43. doi: 10.6052/0459-1879-17-336 Liu J B,Tan H,Bao X,Wang D Y,Li S T. 2018. The seismic wave input method for soil-structure dynamic interaction analysis based on the substructure of artificial boundaries[J]. Chinese Journal of Theoretical and Applied Mechanics,50(1):32–43 (in Chinese).
荣棉水,李小军,卢滔,黄雅虹,吕悦军. 2013. 对含厚软表层海域工程场地设计地震动参数确定的一点建议[J]. 地震学报,35(2):262–271. doi: 10.3969/j.issn.0253-3782.2013.02.012 Rong M S,Li X J,Lu T,Huang Y H,Lü Y J. 2013. Suggestion on determination of design ground motion parameters for offshore engineering sites with deep soft surface layers[J]. Acta Seismologica Sinica,35(2):262–271 (in Chinese).
杨铭,胡进军,谭景阳,公茂盛. 2020. 日本DONET1海域地震动数据及其特征初步分析[J]. 地震工程与工程振动,40(3):139–147. Yang M,Hu J J,Tan J Y,Gong M S. 2020. Offshore ground motion data in DONET1 of Japan and preliminary analysis on its characteristics[J]. Earthquake Engineering and Engineering Vibration,40(3):139–147 (in Chinese).
郑天愉,姚振兴,谢礼立. 1985. 海底强地面运动计算[J]. 地震工程与工程振动,5(3):13–22. Zheng T Y,Yao Z X,Xie L L. 1985. Strong motions of ocean bottom[J]. Earthquake Engineering and Engineering Vibration,5(3):13–22 (in Chinese).
朱镜清. 1988. 地震作用下海水与海床土的耦合运动[J]. 地震工程与工程振动,8(2):37–43. Zhu J Q. 1988. Coupled motion between sea water and sea bed-soil under earthquake action[J]. Earthquake Engineering and Engineering Vibration,8(2):37–43 (in Chinese).
Bao X,Liu J B,Li S T,Wang F,Wang P G. 2020. Seismic response analysis of the reef-seawater system under obliquely incident P and SV waves[J]. Ocean Eng,200:107021. doi: 10.1016/j.oceaneng.2020.107021
Bielak J,Loukakis K,Hisada Y,Yoshimura C. 2003. Domain reduction method for three-dimensional earthquake modeling in localized regions,part I:Theory[J]. Bull Seismol Soc Am,93(2):817–824. doi: 10.1785/0120010251
Brekhovskikh L M. 1980. Waves in Layered Media[M]. 2nd ed. New York: Academic Press.
Chen B K,Wang D S,Li H N,Sun Z G,Li C. 2017. Vertical-to-horizontal response spectral ratio for offshore ground motions:Analysis and simplified design equation[J]. J Cent South Univ,24(1):203–216. doi: 10.1007/s11771-017-3421-0
Li C,Hao H,Li H N,Bi K M,Chen B K. 2017. Modeling and simulation of spatially correlated ground motions at multiple onshore and offshore sites[J]. J Earthq Eng,21(3):359–383. doi: 10.1080/13632469.2016.1172375
Lindsay R B. 1939. Filtration of oblique elastic waves in stratified media[J]. J Acoust Soc Am,11(2):178–183. doi: 10.1121/1.1916021
Link G,Kaltenbacher M,Breuer M,Döllinger M. 2009. A 2D finite-element scheme for fluid-solid-acoustic interactions and its application to human phonation[J]. Comput Methods Appl Mech Eng,198(41/42/43/44):3321–3334.
Okamoto T,Takenaka H,Nakamura T,Hara T. 2017. FDM simulation of earthquakes off western Kyushu,Japan,using a land-ocean unified 3D structure model[J]. Earth Planets Space,69:88. doi: 10.1186/s40623-017-0672-9
Yoshimura C,Bielak J,Hisada Y,Fernández A. 2003. Domain reduction method for three-dimensional earthquake modeling in localized regions,part II:Verification and applications[J]. Bull Seismol Soc Am,93(2):825–841. doi: 10.1785/0120010252
-
期刊类型引用(8)
1. 叶青，王晓，杜学彬，解滔，范晔，周振贵，刘高川. 中国地震井下地电阻率研究进展. 吉林大学学报(地球科学版). 2022(03): 669-683 . 
百度学术
2. Guoze ZHAO，Xuemin ZHANG，Juntao CAI，Yan ZHAN，Qinzhong MA，Ji TANG，Xuebin DU，Bing HAN，Lifeng WANG，Xiaobin CHEN，Qibin XIAO，Xiangyu SUN，Zeyi DONG，Jijun WANG，Jihong ZHANG，Ye FAN，Tao YE. A review of seismo-electromagnetic research in China. Science China(Earth Sciences). 2022(07): 1229-1246 . 必应学术
3. 赵国泽，张学民，蔡军涛，詹艳，马钦忠，汤吉，杜学彬，韩冰，王立凤，陈小斌，肖骑彬，孙翔宇，董泽义，王继军，张继红，范晔，叶涛. 中国地震电磁研究现状和发展趋势. 中国科学:地球科学. 2022(08): 1499-1515 . 百度学术
4. 鲍海英，夏忠，毕雪梅，李鸿宇，王昕. 江苏海安顾庄观测站井下地电阻率观测影响系数分析. 地震. 2022(03): 192-203 . 百度学术
5. 黄明威，张琪，李梦莹，杨牧萍. 新城子井下地电阻率观测影响系数分析. 防灾减灾学报. 2020(01): 49-56 . 百度学术
6. 高曙德，罗维斌，张博，武善艺，李娜，曹玲玲，窦喜英. 编码源地电观测对干扰信号的抑制试验——以通渭台观测为例. 地震工程学报. 2020(05): 1096-1103 . 百度学术
7. 毛先进，段炜，庄儒新，杨玲英，赵晋民. 井下地电阻率观测中布极参数的确定方法. 地震研究. 2019(01): 96-101 . 百度学术
8. 樊晓春，解滔，吴帆，袁慎杰. 井下地电阻率观测影响系数分析——以江宁地震台为例. 中国地震. 2019(02): 347-358 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 347
- HTML全文浏览量: 135
- PDF下载量: 81
- 被引次数: 9
