The Preliminary Analysis of Airguns Seismic Data in the MS6.2 Taiwan Strait Earthquake Region
-
摘要: 对2019年在台湾海峡6.2级地震震区布设的N01、NE02测线地震测深剖面的Pg波走时进行层析成像反演，获得测线下方地壳上部二维速度结构。对N01测线单道反射地震测深剖面进行多次波衰减等处理，并与Pg波成像结果进行对比。研究结果表明，采用走时层析成像方法与单道反射地震测深剖面获得沿探测剖面沉积层上地壳基底形态特征等具有较好的对应关系。由于测线穿越多个地质构造单元和多条断裂带，走时层析成像和单道反射地震测深剖面结果综合显示研究区结晶基底面起伏较大，沉积层速度和厚度变化较剧烈，受台湾海峡西部新生代构造活动影响，显示了相应的断裂或不同地质构造单元在上部地壳内的结构特征。Abstract: The Pg wave travel time of N01 and NE02 deep seismic sounding profiles laid in the MS6.2 Taiwan Strait earthquake region in 2019 was inverted by tomography, and the two-dimensional velocity structure of the upper crust below the survey lines was obtained. The N01 single seismic section was processed by multiple attenuation, and compared with the Pg wave imaging results, there is a good correspondence between the morphological characteristics of the basement. As the survey lines pass through multiple geological tectonic units and multiple fault zones, the results of tomography and single seismic section show that the crystalline basement surface in the study area fluctuates greatly, and the velocity and thickness of the sedimentary layer changes violently. Affected by cenozoic tectonic activities in the western Taiwan Strait, it shows the structural characteristics of the corresponding faults or different geological tectonic units in the upper crust.
-
Key words:
- Taiwan strait /
- Deep seismic exploration /
- Basement structure /
- Seismic profile
-
引言
21世纪初以来，伴随着城市化进程的不断加快，城市中可以使用的空间越来越少，通过建设和开发地下空间工程来提高城市空间利用率已经成为世界性的发展趋势。地下工程包括地下综合管廊、地下停车场、轨道交通、地下污水处理厂、地下商业综合体等（钱七虎，2017），其对于解决城市化进程中出现的土地紧张、环境污染、交通拥堵和能源浪费等问题都发挥着积极作用（陈晓强等，2010）。然而，地下结构的地震安全性并没有引起人们的重视，以往学者们普遍认为地下结构完全被土体所包围，地震时地下结构比地面结构安全（陈国兴等，2016），因而地下结构的抗震设计并没有受到充分的重视，这也导致地下结构抗震研究一直停滞不前。在1995年日本阪神大地震中，地下结构发生了严重破坏（Iida等，1996），科研人员才开始真正重视地下结构的抗震问题。
一直以来，中国缺乏专门针对地下结构进行抗震设计的统一规范，各种规范中所考虑的方法与参数各不相同。早期，大部分地下结构抗震设计沿用地面结构的静力设计方法，比如《铁路工程抗震设计规范（GB 50111—2006）》（中华人民共和国建设部，2006），就是用地面结构抗震的思想进行地下结构设计，但地下结构在地震过程中受地基土约束，变形也受土体控制（林皋, 1990a, 1990b），在抗震原理上与地面结构大不相同（Hashash等，2001）。
美国和日本对地下结构抗震开展研究较早，基于结构变形受土体变形控制这一核心思想，率先提出了很多实用的关于地下结构抗震简化设计方法（权登州等，2015），包括自由场变形法（Wang，1993）、反应位移法（川岛一彦，1994）和柔度系数法（Penzien，2000）。其中以反应位移法推广度最高，但在反应位移法的使用过程中，地基弹簧系数直接决定了计算结果是否合理，弹簧系数的取值需要通过6次有限元计算才能获得，计算过程比较繁琐且不准确。在规避弹簧系数的取值问题上，国内外诸多学者提出了改进方法，如片山幾夫等（1985）提出了反应加速度法；Tateishi（2005）提出了强制反应位移法和改进方法；刘如山等（2007）提出了有限元反应应力法，即将一维土层剪应力沿竖向微分，再将其作为水平体荷载离散到有限元节点上；刘晶波等(2008, 2013)提出了地下结构Pushover分析法和整体式反应位移法，前者借鉴地面结构Pushover方法，施加倒三角形荷载同时推动带结构的土体和自由场土体，通过自由场土体变形控制位移，后者直接在挖除结构的土-结构边界上施加位移，再把反力施加到带有结构的土体。这些简化分析方法为地下结构抗震研究提供了思路，但其适用性和精度还有待进一步检验。
本文以有限元软件ABAQUS为计算平台，以1个双跨箱形地下结构为算例，通过反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度4种方法进行对比计算，再以动力时程分析的结果为参考，分析4种简化分析方法所得弯矩、剪力和相对位移的计算精度，验证4种简化分析方法的适用性与准确性。
1. 常用地下结构抗震分析方法
1.1 动力时程分析法
动力时程分析方法以有限元软件为计算平台，建立土-结构相互作用模型。其基本模型是将土分割为二维平面应变有限元，采用梁单元模拟地下结构并与土体连接，仅考虑由场地底部往上传播的剪切波，侧边界采用固定竖向、释放水平向的方式，底边界则固定竖向，再从模型的底部输入地震动，进而分析得到水平地震作用下土-地下结构的地震反应规律，计算模型如图 1所示。动力时程分析法是采用逐步积分的方式求解地下结构在地震荷载作用下每一时刻的动力响应，能够较好地考虑材料的非线性、地震动频谱、持时和幅值情况。动力时程分析方法能够较准确地反应地下结构的地震反应，该方法常被作为其它简化分析方法的校核标准，然而其对土体动力本构关系要求较高，计算非线性问题耗时较长，求解难度较大，因而对使用人员的专业知识要求较高，同时还需要具备较好的计算平台，通常只被应用于重大工程项目的抗震分析中。国内地下建筑抗震设计规范《地下铁道建筑结构抗震设计规范（DG/TJ 08-2064-2009）》（同济大学，2009）和《城市轨道交通结构抗震设计规范（GB 50909—2014）》（中华人民共和国住房和城乡建设部，2014）都引入了动力时程分析法。
1.2 反应位移法
反应位移法采用梁-弹簧单元建模（不包含土体），其中梁单元模拟地下结构，弹簧单元（正向和切向）模拟结构和地基土之间的相互作用。首先，计算出结构顶、底板位置处所对应土层最大位移差，再把此时刻对应的位移差施加在弹簧远离结构的一端，同时在自由场计算中提取出该时刻对应的结构惯性力和地层剪应力，并直接施加在梁单元模型上，采用静力的方法计算结构响应，其模型见图 2。反应位移法的计算结果依赖于地基弹簧刚度的取值，通过有限元计算弹簧刚度的方法也比较繁琐，且其取值是否合理也有待商榷。同样，《地下铁道建筑结构抗震设计规范（DG/TJ 08-2064-2009）》和《城市轨道交通结构抗震设计规范（GB 50909—2014）》都借鉴了反应位移法。
1.3 强制反应位移法
土体边界施加强制变形是反应位移法的1种简化形式，进行简化主要是为了规避弹簧刚度的取值问题，计算过程和反应位移法相同。首先，需要通过自由场计算得到地震荷载下结构顶板和底板所处土层的最大位移差，此时结构变形最大，即处于极限响应状态，然后将该时刻自由场计算模型中全部土层的相对位移提取出来，并以强制位移的模式施加在有限元模型侧边界上，将动力问题简化为静力问题，从而得到结构响应，其计算模型如图 3所示，其中W为地下结构宽度，D为地下结构到有限元模型侧边的距离。Tateishi（2005）指出，采用强制反应位移法处理动力作用虽然能够满足自由场应变在静力模型边界及附近准确传递，但远离模型边界处会因为土体阻尼的影响而发生衰减现象，尤以结构周围土层最为明显，应变值可衰减20%—50%。由于地下结构刚度通常大于土体刚度，因而结构变形比自由场变形小，如果选取的侧边距过小，土体阻尼作用将无法得到发挥，则相当于把自由场变形直接施加到了地下结构上，从而导致地下结构响应过大。由此可见，找到合适的侧边距是强制反应位移法准确求解结构响应的关键。
1.4 反应加速度法
反应位移法也可以通过反应加速度法进行简化，与强制反应位移法相同，均是为了规避弹簧刚度取值失真产生较大的计算误差。首先，通过自由场计算给出地震荷载下结构顶板和底板所处土层的最大位移差，然后把该时刻模型中全部土层的加速度取出，最后将该加速度结果通过水平体力的方式赋予整个模型土体和地下结构，通过静力的方法给出结构响应，其计算模型如图 4所示。近年来，反应加速度法在国内得以大力推广，《城市轨道交通结构抗震设计规范（GB 50909—2014）》首次引入了反应加速度法。
1.5 仅在土地施加加速度的简化方法
Tateishi（2005）提出了强制反应位移法的改进方法，对整个土层强制施加水平体力从而让土层变形与自由场变形一致，即仅把水平加速度按体力的形式施加到土体上。仅在土体施加水平加速度也是对反应加速度法的进一步简化，计算求解的过程和反应加速度法一致，唯一的区别在于反应加速度法把最大变形状态时的土层加速度以体力的形式施加到土体和结构上，而该方法则进一步简化了这一处理过程，其模型见图 5。由于地下结构一般内部挖空，仅有周围一圈衬砌，纵向单位长度的隧道质量远远小于原土体的质量，加速度在结构上产生的体力较小，因此理论上对结构的地震响应影响甚微（杨智勇等，2012）。
2. 算例分析
2.1 模型建立
以某单一土层中埋深7m的双跨箱形地下结构为算例，计算其在地震荷载下的响应，该地下结构横截面尺寸为8m×4m，墙厚0.8m，中柱厚0.8m。依次采用动力时程法、反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度这5种方法计算该地下结构的内力和中柱顶底相对变形。本例中结构的横截面抗震分析可按平面应变问题处理，土-结构接触面采用ABAQUS软件中自带的绑定接触，土体与结构参数见表 1。使用实体单元对土体进行建模，地下结构则使用梁单元进行建模并与土体连接，土体选取摩尔-库伦本构模型，结构则取为线弹性本构模型。
表 1 土体与结构的物理力学参数Table 1. Physical and mechanical parameters of soil and structure类型 密度/g·cm-3 弹性模量/MPa 泊松比 内摩擦角/° 粘聚力/kPa 土体 1.9 30 0.3 30 30 结构 2.3 30000 0.167 采用动力时程分析法时，分别将结构宽度和高度的4倍取为侧边距和底边距，计算结果用作校核，模型见图 6、图 7。反应位移法不建立土体，改变模型尺寸进行强制反应位移法、反应加速度法和仅在土体施加加速度方法的试算。
采用动力时程分析法时，侧边界固定竖向、水平向释放，底边界为粘性边界，从模型底部输入峰值加速度为0.2g的El-Centro波，如图 8所示。然后，再采用峰值加速度为0.2g的Kobe波（图 9）验证分析所得的规律。自由场模型与土-结构动力时程分析模型采用相同的土体本构模型和土体参数，0.2g峰值加速度下土层位移差最大时刻的相对位移分布和加速度分布见图 10、图 11。
图 10 自由场最不利水平相对位移分布Figure 10. Least favorable displacement distributions of free-field2.2 结构响应分析
2.2.1 El-Centro波下不同侧边距结果分析
首先，以计算出1的动力时程分析法结果作为后续简化分析方法的校核标准。根据《城市轨道交通结构抗震设计规范（GB 50909—2014）》要求，动力时程分析的侧边距应大于结构宽度的3倍及以上，本次计算模型侧边距取为结构宽度的4倍；其次，采用反应位移法单独建立梁-弹簧单元对结构响应进行分析；最后，改变模型土体尺寸，用3种需要建立土体的简化分析方法对结构响应进行计算。输入峰值加速度为0.2g的El-Centro波，则控制截面A、B、C、D的弯矩、剪力和顶、底相对位移计算结果见表 2、表 3和表 4，地下结构计算误差分析见图 12。
表 2 不同侧边距的弯矩计算结果Table 2. Bending moment of structure with different lateral marginsD/W 使用方法 弯矩/kN·m 截面A 截面B 截面C 截面D 1倍 强制反应位移法 292.9 -291.3 280.8 353.9 反应加速度法 319.4 -324.5 294.0 367.1 仅在土体施加加速度 296.5 -300.6 278.4 351.8 2倍 强制反应位移法 179.5 -181.3 200.7 282.1 反应加速度法 328.1 -33.3 301.0 372.2 仅在土体施加加速度 309.3 -313.6 287.9 360.0 3倍 强制反应位移法 100.4 -101.3 147.7 224.7 反应加速度法 331.7 -336.9 303.8 374.9 仅在土体施加加速度 314.5 -18.9 291.8 363.9 4倍 强制反应位移法 50.2 -50.4 113.8 187.4 反应加速度法 332.6 -37.8 304.3 375.5 仅在土体施加加速度 316.2 -20.6 292.8 365.0 反应位移法 286.5 -290.1 201.0 199.8 动力分析 282.5 -89.1 272.3 345.2 表 3 不同侧边距的剪力计算结果Table 3. Shear force of structure with different lateral marginsD/W 使用方法 剪力/kN 截面A 截面B 截面C 截面D 1倍 强制反应位移法 -182.6 -182.6 -213.5 -263.8 反应加速度法 -195.8 -204.0 -205.6 -271.7 仅在土体施加加速度 -186.6 -186.6 -196.8 -266.3 2倍 强制反应位移法 -112.8 -112.8 -169.0 -249.3 反应加速度法 -201.3 -209.4 -210.0 -274.3 仅在土体施加加速度 -194.7 -194.7 -202.0 -270.7 3倍 强制反应位移法 -63.0 -63.0 -144.4 -219.7 反应加速度法 -203.6 -211.7 -211.9 -276.2 仅在土体施加加速度 -198.0 -198.0 -204.2 -273.3 4倍 强制反应位移法 -31.4 -31.4 -128.2 -199.6 反应加速度法 -204.1 -212.3 -212.0 -276.4 仅在土体施加加速度 -199.0 -199.0 -204.5 -273.8 反应位移法 -174.8 -183.0 -106.2 -106.2 动力分析 -175.7 -180.6 -194.0 -266.7 表 4 不同侧边距的结构变形计算结果Table 4. Deformation of structure with different lateral marginsD/W 使用方法 中柱顶、底相对位移/mm 1倍 强制反应位移法 2.396 反应加速度法 5.011 仅在土体施加加速度 4.636 2倍 强制反应位移法 0.302 反应加速度法 5.340 仅在土体施加加速度 5.026 3倍 强制反应位移法 0.584 反应加速度法 5.354 仅在土体施加加速度 5.071 4倍 强制反应位移法 0.722 反应加速度法 5.346 仅在土体施加加速度 5.077 反应位移法 3.628 动力分析 5.596 
图 12 3种带土体简化分析方法的计算误差Figure 12. Calculation errors of three simplified methods with soil model在输入0.2g峰值加速度El-Centro波的情况下，反应位移法的计算精度在结构中柱部分较好。对比反应位移法和动力时程分析方法可知，弯矩误差在中柱上端为1.41%，下端最小仅为0.34%，剪力误差在中柱上端最小仅为0.51%，下端为1.32%。但是在结构的角部，计算误差较大，剪力和弯矩计算值都偏小，在C截面（右侧墙上端）弯矩误差值达到26.18%，剪力误差值达到45.25%，尤其D截面（左侧墙下端），弯矩误差值最大达到42.12%，剪力误差值最大达到60.17%。弯矩平均误差为13.8%，剪力平均误差23.3%。
强制反应位移法的计算精度依赖于静力模型土体侧边距的选取，随着土体侧边距取值的不断增加，地下结构的弯矩、剪力计算误差会相应的增大。当侧边距取为4倍结构宽度时，弯矩平均误差最大达到70.2%，剪力平均误差最大达到58.1%；然而，当侧边距取为1倍结构宽度时，弯矩平均误差仅为2.8%，剪力平均误差为2.5%。当侧边距为1倍结构宽度时，中柱顶、底相对变形误差为57.18%，侧边距为2倍结构宽度时误差最大为94.60%，但随着侧边距的增大，误差有所减小。
反应加速度法与仅在土体施加加速度的简化分析方法计算所得弯矩、剪力和中柱变形随着模型侧边距的增大，计算误差变化幅度很小。当侧边距取为1倍结构宽度时，反应加速度法的弯矩平均误差为10.1%，剪力平均误差为8.7%，中柱变形误差为10.5%，仅在土体施加加速度方法的弯矩平均误差为3.1%，剪力平均误差为2.4%，中柱变形误差为17.2%；当侧边距取为4倍结构宽度时，反应加速度法的弯矩平均误差为14.3%，剪力平均误差为12.7%，中柱变形误差为4.5%，仅在土体施加加速度方法的弯矩平均误差为9.2%，剪力平均误差为7.8%，中柱变形误差为9.3%。由此可见，以加速度模拟地震作用时土体侧边距对计算误差影响较小，2种方法变化趋势基本同步。由于仅在土体施加加速度的简化分析方法忽略了结构加速度，相同工况下该方法的计算数值略小于反应加速度法。施加加速度的2种方法对于结构变形的计算精度均较高，反应位移法次之，强制反应位移法的计算结果较不理想。
2.2.2 Kobe地震波验证分析
为了消除地震动不同带来的偶然性，改用0.2g峰值加速度的Kobe地震波验证上述分析规律。
从前文分析可知，反应位移法在应用过程中不需要建立土体；强制反应位移法在侧边距取为1倍结构宽度时计算误差较小；反应加速度法和仅在土体施加加速度的简化方法对于侧边界取值不敏感。因此，需要建立土体的3种静力模型统一采用1倍结构宽度作为侧边距。弯矩、剪力、中柱变形计算结果见表 5、表 6、表 7。
表 5 Kobe波下弯矩计算结果Table 5. Bending moment calculation results of structure under Kobe seismic wave使用方法 弯矩/kN·m 平均误差/% 截面A 截面B 截面C 截面D 动力分析 375.3 -383.7 332.1 410.6 反应位移法 400.2 -403.7 284.3 276.0 10.5 强制反应位移法 298.8 -310.4 275.9 364.3 18.0 反应加速度法 341.9 -349.2 308.4 384.7 8.0 仅在土体施加加速度 312.4 -317.7 288.0 364.5 15.0 表 6 Kobe波下剪力计算结果Table 6. Shear force calculation results of structure under Kobe seismic wave使用方法 剪力/kN 平均误差/% 截面A 截面B 截面C 截面D 动力分析 -230.8 -242.9 -218.3 -299.0 反应位移法 -242.6 -257.5 -159.9 -139.9 13.7 强制反应位移法 -190.4 -190.4 -180.2 -290.7 27.0 反应加速度法 -207.4 -222.2 -209.8 -281.2 4.4 仅在土体施加加速度 -196.9 -196.9 -198.1 -275.4 5.6 表 7 Kobe波下结构变形计算结果Table 7. Deformation calculation results of structure under Kobe seismic wave使用方法 中柱顶、底相对位移/mm 平均误差/% 动力分析 5.808 反应位移法 4.768 17.9 强制反应位移法 3.232 44.4 反应加速度法 5.371 7.5 仅在土体施加加速度 4.900 15.6 对于不同类型的地震动，4种简化分析方法的计算精度均比较稳定，弯矩误差在10%左右，反应加速度法平均弯矩误差最小，为8.0%，强制反应位移法平均弯矩误差最大，为18.0%；2种施加加速度简化分析方法的剪力误差极小，分别为4.4%和5.6%，强制反应位移法剪力误差最大，仅为27%；强制反应位移法的结构变形误差最大达到44.4%，其余3种方法精度较好，反应加速度法误差最小，仅为7.5%。
3. 4种简化分析方法的特点分析
4种简化分析方法的相同点在于脱离了动力分析，能够按静力的方式计算结构响应，使得计算效率大大提高，且在不同地震动输入条件下都能得到较好的精度；另外，4种简化分析方法均需要进行自由场动力时程分析，从而获取地层位移或者地层加速度作为静力计算的地震输入荷载，自由场变形是4种简化分析方法的基础，但是结构存在对自由场的影响却被忽略了。
4种简化分析方法的不同点在于，①采用反应位移法不需要建立模型土体，以弹簧单元模拟土-结构相互作用，但是地基弹簧刚度取值需要进行6次有限元计算，过程比较繁琐。并且，土体弹簧刚度值对结构内力计算有巨大影响，通过有限元计算弹簧刚度的方法无法反应土层自身的相互作用，计算过程中土体的等效弹簧之间是不相关的，造成的结果则是土层对结构特别是角部的约束作用减小，计算所得的角部内力值偏小。②强制反应位移法对静力模型侧边界宽度值的选取非常敏感，侧边距的选取对于结构内力计算精度有重大影响。从计算结果来看，采用1倍结构宽度作为侧边距时结构内力的计算误差最小。但是无论侧边距如何变化，对结构变形的计算误差均较大。在对结构内力进行校核时，采用该方法计算的结果较为可靠。③反应加速度法和仅在土层施加加速度的简化方法对于侧边界的选取不敏感，2种方法对于不同的侧边距得到的计算精度均比较稳定，计算精度随着侧边距的增加仅有轻微的变化，2种方法对结构内力、变形的计算误差均较小。
4. 结论
本文通过反应位移法、强制反应位移法、反应加速度法和仅在土体施加加速度4种简化分析方法，分别计算了1个双跨箱形地下结构在0.2g峰值加速度的EL-Centro波和Kobe波作用下的地震响应，再利用时程分析法进行校核，经过对比，得出以下结论：
（1）总体而言，4种简化分析方法计算得到的结构在地震作用下的内力、变形值接近，规律一致，不同类型的地震动作用下均能得到较好的计算精度，能合理地反映地下结构的动力特性。
（2）若要进一步推广使用反应位移法，则须使其基地弹簧系数的取值快捷、可靠。
（3）模拟地震等效荷载的自由场土层位移和土层加速度具有较好的计算精度，且便于施加，但是需要进一步考虑地下结构存在对自由场的影响。
（4）在使用强制反应位移法进行结构内力校核时，侧边距取为1倍结构宽度时效果最好，计算误差最小。而反应加速度法和仅在土体施加加速度的方法对侧边距的选取不敏感。
（5）强制反应位移法、反应加速度法和仅在土体施加加速度3种带土体模型的简化分析方法均省略了反应位移法中地基弹簧刚度系数的取值，同时更加真实地反应了土-结构间的协调作用，使用时比反应位移法更简便。
-
图 6 N01测线剖面基底二维速度结构
Figure 6. The basement two-dimensional velocity structure of N01 profile
图 7 NE02测线剖面基底二维速度结构
Figure 7. The basement two-dimensional velocity structure of NE02 profile
-
[1] 蔡锋, 曹超, 周兴华等, 2013. 中国近海海洋—海底地形地貌. 北京: 海洋出版社. [2] 陈园田, 谢志平, 1996. 台湾海峡的活动断裂与地震活动. 华南地震, 16(1): 57—62.Chen Y. T., Xie Z. P., 1996. Active faults and earthquakes in Taiwan Strait. South China Journal of Seismology, 16(1): 57—62. (in Chinese) [3] 邓起东, 闽伟, 晁洪太等, 2001. 渤海地区新生代构造与地震活动. 见: 卢演俦主编. 新构造与环境. 北京: 地震出版社, 218—233. [4] 丁祥焕, 1999. 福建东南沿海活动断裂与地震. 福州: 福建科学技术出版社. [5] 郭令智, 施央申, 马瑞士, 1983. 西太平洋中、新生代活动大陆边缘和岛弧构造的形成及演化. 地质学报, 57(1): 11—21.Guo L. Z., Shi Y. S., Ma R. S., 1983. On the formation and evolution of the Mesozoic-Cenozoic active continental margin and island arc tectonics of the western Pacific Ocean. Acta Geologica Sinica, 1983, 57(1): 11—21. (in Chinese) [6] 郭令智, 马瑞士, 施央申等, 1998. 论西太平洋活动大陆边缘中—新生代弧后盆地的分类和演化. 成都理工学院学报, 25(2): 134—144.Guo L. Z., Ma R. S., Shi Y. S., et al., 1998. On the generation and evolution of Mesozoic Cenozoic back arc basins of western Pacific active continental margin. Journal of Chengdu University of Technology, 25(2): 134—144. (in Chinese) [7] 黄昭, 王善雄, 2006. 台湾海峡滨海断裂带的构造特征与活动性. 大地测量与地球动力学, 26(3): 16—22.Huang Z., Wang S. X., 2006. Tectonic features and activity of Binhai Fault zone in Taiwan Strait. Journal of Geodesy and Geodynamics, 26(3): 16—22. (in Chinese) [8] 李赶先, 1994. 台湾海峡新生代断裂构造及其形成演化. 热带海洋, 13(2): 17—24.Li G. X., 1994. The formation and evolution of Cenozoic rifts Taiwan Strait. Tropic Oceanology, 13(2): 17—24. (in Chinese) [9] 戚筱俊, 1999. 台湾海峡地形、地质及地震概述. 西部探矿工程, 11(4): 4—6. [10] 单联瑜, 刘连升, 2007. 三维高阶快速步进波前重构旅行时计算精度分析. 石油地球物理勘探, 42(4): 413—417. doi: 10.3321/j.issn:1000-7210.2007.04.010Shan L. Y., Liu L. S., 2007. Analysis on computation precision of wavefront-reconstructing traveltime for 3-D high-order rapid stepping algorithm. Oil Geophysical Prospecting, 42(4): 413—417. (in Chinese) doi: 10.3321/j.issn:1000-7210.2007.04.010 [11] 王本善, 1995. 台湾海峡地区地质构造特征. 南方油气地质, 2(1): 57—64. [12] 王建, 赵明辉, 贺恩远等, 2014. 初至波层析成像的反演参数选取: 以南海中央次海盆三维地震探测数据为例. 热带海洋学报, 33(5): 74—83. doi: 10.3969/j.issn.1009-5470.2014.05.010Wang J., Zhao M. H., He E. Y., et al., 2014. The selection of optimal inversion parameters for first-arrival seismic tomography: an application to 3D seismic data from the central sub-basin of the South China Sea. Journal of Tropical Oceanography, 33(5): 74-83. (in Chinese) doi: 10.3969/j.issn.1009-5470.2014.05.010 [13] 王笋, 丘学林, 方伟华等, 2018. 台湾海峡西南部的海陆联合深地震探测资料特点与处理对策. 热带海洋学报, 37(2): 92—99.Wang S., Qiu X. L., Fang W. H., et al., 2018. Features of the onshore-offshore seismic data in Southwest Taiwan Strait and some countermeasures for data processing. Journal of Tropical Oceanography, 37(2): 92—99. (in Chinese) [14] 熊绍柏, 刘宏兵, 王有学等, 2002. 华南上地壳速度分布与基底、盖层构造研究. 地球物理学报, 45(6): 784—791. doi: 10.3321/j.issn:0001-5733.2002.06.005Xiong S. B., Liu H. B., Wang Y. X., et al., 2002. A study on velocity distribution in upper crust and tectonics of basement and cover in South China. Chinese Journal of Geophysics, 45(6): 784—791. (in Chinese) doi: 10.3321/j.issn:0001-5733.2002.06.005 [15] 闫培, 张艺峰, 李普春等, 2015. 台湾海峡深地震探测: HX9测线试验及初步成果. 中国地震, 31(1): 27—35. doi: 10.3969/j.issn.1001-4683.2015.01.003Yan P., Zhang Y. F., Li P. C., et al., 2015. Deep seismic exploration in the Taiwan Strait: the experiment of HX9 survey line and preliminary results. Earthquake Research in China, 31(1): 27—35. (in Chinese) doi: 10.3969/j.issn.1001-4683.2015.01.003 [16] 杨肖琪, 宋文隆, 陈承惠, 1996. 台湾海峡地质构造特征. 台湾海峡, 15(2): 127—136.Yang X. Q., Song W. L., Chen C. H., 1996. Characteristics of geological structure in Taiwan Strait. Journal of Oceanography in Taiwan Strait, 15(2): 127—136. (in Chinese) [17] 姚道平, 张艺峰, 闫培等, 2016. 台湾海峡大容量气枪震源海陆联测初探. 地震学报, 38(2): 167—178. doi: 10.11939/jass.2016.02.002Yao D. P., Zhang Y. F., Yan P., et al., 2016. A preliminary exploration into onshore-offshore seismic experiment by using large volume air-gun source in Taiwan Strait. Acta Seismologica Sinica, 38(2): 167—178. (in Chinese) doi: 10.11939/jass.2016.02.002 [18] 俞何兴, 陈汝勤, 1996. 台湾海域之沉积盆地. 台北: 台湾编译馆. [19] 俞何兴, 2006. 台湾海域海底地形及其地质意义. 台北: 台湾大学. [20] 詹文欢, 孙宗勋, 唐诚等, 2004. 华南滨海断裂带及其对台湾海峡地震活动的控制作用. 热带海洋学报, 23(4): 19—24. doi: 10.3969/j.issn.1009-5470.2004.04.003Zhan W. H., Sun Z. X., Tang C., et al., 2004. Littoral active fault belt of south China and its control on seismic activity in Taiwan Strait. Journal of Tropical Oceanography, 23(4): 19—24. (in Chinese) doi: 10.3969/j.issn.1009-5470.2004.04.003 [21] 张艺峰, 金星, 郭晓然等, 2017. 台湾海峡地区深部结构探测数据处理及初步结果. 见: 2017中国地球科学联合学术年会论文集. 北京: 中国地球物理学会. [22] Paige C. C., Saunders M. A., 1982. LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software, 8(1): 43—71. doi: 10.1145/355984.355989 [23] Rawlinson N., Sambridge M., 2004a. Wave front evolution in strongly heterogeneous layered media using the fast marching method. Geophysical Journal International, 156(3): 631—647. doi: 10.1111/j.1365-246X.2004.02153.x [24] Rawlinson N., Sambridge M., 2004b. Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 69(5): 1338—1350. doi: 10.1190/1.1801950 [25] Sethian J. A., 1996. Theory, algorithms, and applications of level set methods for propagating interfaces. Acta Numerica, 5: 309—395. doi: 10.1017/S0962492900002671 [26] Sethian J. A., Popovici A. M., 1999. 3-D traveltime computation using the fast marching method. Geophysics, 64(2): 516—523. doi: 10.1190/1.1444558 [27] Zelt C. A., Barton P. J., 1998. Three-dimensional seismic refraction tomography: a comparison of two methods applied to data from the Faeroe Basin. Journal of Geophysical Research: Solid Earth, 103(B4): 7187—7210. doi: 10.1029/97JB03536 期刊类型引用(7)
1. 张强，金鑫，刘英君. 浅埋地下结构反应位移法不同位置弹簧刚度系数的影响分析. 市政技术. 2023(02): 113-119 . 
百度学术2. 闫晓京，邢珏蕙，周忠发，朱忠义，孙枫然，周笋，段世昌，閤东东. 北京歌舞剧院大埋深地下结构设计. 建筑结构. 2023(18): 63-68+143 . 百度学术3. 程新俊，许翔，景立平，崔杰，梁海安. 一种基床系数计算方法及其在反应位移法中的应用. 岩土工程学报. 2023(12): 2604-2613 . 百度学术4. 刘春峰，屈永强，许紫刚. 用于地下结构地震反应分析的改进反应加速度法. 震灾防御技术. 2021(04): 710-716 . 本站查看5. 赵密，邵伟昂，黄景琦，杜修力，王媛. 圆形隧道横断面地震响应简化解析方法的数值验证. 工程力学. 2020(09): 84-93 . 百度学术6. 赵密，邵伟昂，黄景琦，杜修力，王媛. 衬砌厚度对隧道抗震简化解析方法精度的影响. 地震工程与工程振动. 2020(06): 47-60 . 百度学术7. 卜永红，权登州，祝河清，石晶. 砂土场地地铁车站抗震计算方法. 铁道建筑. 2019(11): 88-90+100 . 百度学术其他类型引用(8)
-
下载:
下载:
点击查看大图
计量
- 文章访问数: 169
- HTML全文浏览量: 70
- PDF下载量: 14
- 被引次数: 15
