Regional Geomorphic Characteristics and its factors Controlling of the Maerkang M6.0 Earthquake Swarm in Sichuan Province, Eastern Tibetan Plateau
-
摘要: 2022年6月10日四川马尔康6.0级震群地震发生在青藏高原东缘，位于强烈隆升的龙门山及其西缘高原腹地过渡地带，断裂构造复杂，水系发育，小震与中强地震活动频发。基于90 m×90 m空间分辨率的SRTM3数据，利用数字地貌分析方法提取了马尔康6.0级震群所在的脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数等地貌参数，结合降水和岩性数据进行综合分析。研究结果表明，脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数整体表现出西北低、东南高的特征。构造活动是导致区域地貌参数呈现差异性的主要原因，地貌参数高值主要分布于龙日坝断裂两侧，沿松岗断裂两侧展布，龙日坝断裂东侧是构造变形与差异性隆升强烈的龙门山地区；龙日坝断裂西侧是多条断裂交汇的构造复杂区，中强地震活跃，地貌破坏与重塑作用相对强烈，这可能是地貌参数高值的主要原因。Abstract: On June 10, 2022, the Maerkang M6.0 earthquake swarm occurred inside the Baryan Har Block, eastern Tibetan Plateau. The earthquake swarm located in the transition zone between the strongly uplifted Longmen Mountains and the plateau hinterland, which is characterized by complex structures, well-developed water systems, frequent occurrence of moderate and small earthquakes. Based on SRTM3 data with spatial resolution of 90 m × 90 m, the slope, relief amplitude (RA), hypsometric integral (HI) and channel steepness index (ksn) of Jiaomuzu River Basin in the Maerkang M6.0 earthquake swarm areas are extracted by digital geomorphological analysis method. Combined with precipitation and lithology data, the results show that the slope, RA, HI and ksn values of Jiaomuzu River Basin showed the characteristics of low northwest and high southeast. Tectonic activity is the main reason for the difference of regional geomorphic parameters. The high values of the geomorphological parameters are mainly distributed on both sides of the Longriba fault, and along both sides of Songgang fault. The eastern side of the Longriba fault is the Longmen Mountains area that is characterized by the strong tectonic deformation and significant differential uplift. The western side of the Longriba fault is a complex structural area with a characteristics of multiple faults intersect and frequent occurrence of moderate earthquakes, where the geomorphological damage and remodeling are relatively strong. These may be the main reasons for the high values of geomorphic parameters in our study region.
-
引言
软土在我国沿海一带分布很广，如长江三角洲、珠江三角洲、渤海湾以及浙江、福建沿海地区等都存在海相或湖相沉积的软土。软土具有松软、孔隙比大、天然含水量高、压缩性高、强度低、渗透性小和结构性灵敏等特点，在受扰动或震动后极易丧失强度 (吴跃东，2008)。修建在软土上的建筑物在地震荷载作用下，瞬间内出现突发性超量沉陷，以及不均匀沉降，即震陷，从而导致建筑物显著倾斜。有时将这种沉陷称为地震引起的永久变形或残余变形。1976年唐山大地震、1940年埃尔森特罗地震、1960年智利地震、1957年和1985年的墨西哥大地震以及1957年日本宫城冲地震均出现了严重的软土震陷现象，并造成了严重的震害 (中国科学研究院工程力学研究所，1979；刘恢先，1986)。2008年汶川特大地震后在附近地区发现多处由软土地基引起的震害：汉旺镇某商品房建在软土地基上，房屋为砖柱结构，地面发生不均匀沉降，从而导致上部结构破坏；都江堰市电信大楼台阶沉降，最大沉降约8cm。另外，此次地震中发生多处软土震陷现象，其中以库坝的坝体震陷较为多见 (江席苗，2009)。震害现场显示，软土震陷是导致震害的主要原因之一，这引起了国内外学者的广泛关注。他们采用动三轴试验、离心振动台试验等方法对软土震陷问题开展了一系列的研究，取得了一定的成果 (张建毅等，2012；陈国兴等，2004；李冬等，2011；陈建峰等，2008；Narikawa等，1999)。但现有的软土震陷量计算方法在工程应用方面还存在一定的不足之处，我国相关规范中也未对软土震陷计算做出明确的规定。《岩土工程勘察规范》(GB50021—2001) 和《建筑抗震设计规范》(GB50011—2010) 只是在条文中说明可以参考已有的研究成果进行计算估计。最新的《软土地区工程地质勘察规范》(JGJ83—2011) 对震陷的判别和震陷量的估算做出了相应的规定，即场区的基本烈度为Ⅶ度或大于Ⅶ度时，对于采用天然地基的建筑物，甲级类建筑物和对沉降有严格要求的乙级类建筑物应进行专门的震陷分析计算；对沉降无特殊要求的乙级建筑物和对沉降敏感的丙级建筑物参照规范给出的参考值或根据地区经验确定，但是关于震陷量计算方法也没有做出具体规定。
综上所述，软土震陷问题早已被许多地震实例所证实。传统的结构性模型中所建立的土的结构性模型一般都较为复杂，参数较多，并且参数需通过大量的试验获得，对大型工程项目进行软土震陷量计算是可行的，但是对于中小型项目来说仍存在很大的困难，因此在实际应用中会受到一定的限制。本文在前人工作基础上，从理论分析、动三轴试验和震害验证等方面分析软弱土动力特性，提出一种基于动三轴试验结果的、计算准确、参数少并易测得的、工程技术人员易操作的软土震陷简化计算方法。
1. 软土震陷的简化计算方法
目前，分析土体震陷的方法主要分为2大类：第一类是基于Newmark提出的屈服加速度概念的有限滑动位移法；第二类是整体变形分析法，该方法主要涉及等价节点力模型和模量软化模型。国内外很多学者大多数都基于上述2类研究方法开展研究 (Monismith等，1975；Lee，1974；谢君斐等，1981；1988；郁寿松等，1989；周健等，2000；王建华等，1993；孟上九等，2002；Seed等，1966)，发展及改进了相关软化模型、残余应变模型，但很多方法都需要采用有限元方法进行计算。这些方法在工程中对较大型的项目是可行的，但计算和试验成本较高，使其在中小型项目中的应用存在困难。杨石红等 (1997)提出了软弱地基土层的震陷简化计算方法，该方法相对二维有限元计算方法简化了许多，但计算步骤繁杂，推导过程存在缺少试验数据支撑的缺陷，且土的应变势是引用其他研究成果，因此有必要在此基础上结合动三轴试验结果对该方法进行改进。本文提出基于“软化模型”概念的软土震陷简化计算方法，简化思路是通过只计算土层在静力和地震荷载下的压缩量，然后附加修正系数来估算土层的震陷量。本文提出的基本计算步骤简述如下：
(1) 根据实际需要确定震陷量计算所对应的地震动水平；
(2) 确定各层土的应力状态：计算土的自重应力、地震剪应力；
(3) 静三轴固结排水剪切试验确定土的邓肯参数；
(4) 确定各层土的变形模量；
(5) 计算土层沉降量；
(6) 根据地震动水平选取合适的震陷量修正系数；
(7) 确定场地的震陷量。
2. 震陷量修正系数
目前，震陷计算需要确定土层在地震作用下的残余应变值，但其方法较为复杂。本计算方法采用带有区域经验性质的震陷量修正系数进行计算，以避免对每个工程都要进行大量的动三轴试验确定残余应变值，从而达到简化计算的目的，使一般工程技术人员易于掌握使用。本文的震陷量修正系数定义为：在考虑残余应变情况下，计算土层静力和地震荷载下土层的压缩量与只计算土层在静力作用下的压缩量的比值。本文的研究思路为：首先选择典型软土场地，其次进行场地勘察取样，并将原状土样进行土物理力学参数、静动三轴试验，最后根据试验数据分别计算上述2种情况下的压缩量，进行对比后得出震陷量修正系数。
根据以往震害调查结果，本研究将研究场地选定在天津滨海新区南海路。该场地临近1976年唐山7.8级地震中震陷较严重的塘沽望海楼及建岗村，因此既具有现实研究意义，又方便对计算方法所得结果作对比检验。该场地为华北平原东部滨海平原地貌，属海洋与陆相交互沉积地层，地形平坦，受人类活动影响，表层土为人工填土。本研究的钻探取样和运输均按照《岩土勘察规范》的要求进行，采用薄壁取土器以保证取土质量，并进行了钻孔波速测试等相关工作。钻孔的等效剪切波速VS=143m/s，覆盖层厚度超过80m，根据《建筑抗震设计规范》(GB50011－2010)，判定该场地土类型属软弱土，场地类别为Ⅳ类。共钻取了6m、6.5m、7m、7.5m、9.5m、10.5m、11.5m、12.5m、13.5m共9个深度处原状土样。
2.1 静力学试验
为了掌握土样的静力学特性，对所采原状土样进行了基本物理力学参数的测定，密度测试采用环刀法，含水率测试采用烘干法，液塑限测试采用联合测定法。从试验结果可知，所取土样的孔隙比都大于1，并且含水率都大于液限，根据《岩土勘察规范》相关规定试验土样为软弱土，因此所取土样可以用于软土震陷动三轴试验研究。为了得到软土的邓肯参数，通过土的三轴固结排水剪切试验量测天津滨海地区软土在3种围压下的应力差、轴向变形和体变，以此来确定土的邓肯参数。试验结果如表 1所示。
表 1 天津软土的邓肯参数Table 1. Duncan parameters of soft soil in Tianjin土类 Ks/kPa Ns 摩擦角/° 粘聚力/kPa 破坏比 淤泥质土 1237 0.465 12 20 0.478 淤泥质粉质粘土 2930 0.39 12 44 0.687 2.2 动三轴试验
动三轴试验采用美国GCTS公司所生产的STX-200型电液伺服控制双向动三轴测试系统，配备CATS试验控制软件。仪器最大轴向加载力为100kN，最大围压2MPa，加载频率0.01—20Hz。试验用土为钻探中所取得的原状软土，现场取土后用样盒封装后运至实验室保存，试验时采用削土器将原状土削成直径为38mm、高为76mm的土柱。装样后分级施加固结压力对试样进行固结，固结完成后在不排水的条件下施加轴向动荷载激振。在激振过程中，观察试样状态并利用控制软件记录振次、轴压、应变等数据。
试验主要研究的是天津软土的残余应变随振次的发展规律，共对6.0m、6.5m、7.0m、7.5m、9.5m、11.5m、12.5m、13.5m共8个不同深度的16组土样进行动三轴试验。采用100kPa的固结压力对试样进行等压固结，动荷载采用频率为1Hz的正弦波，测定试样在不同动应力加载下的残余应变值。为获取更多的数据，当试样动应变达到15%或振次达到200周时结束试验。
2.2.1 原始激振现象
由于原状试样的土质不均匀，虽然按规范划分均属于淤泥质土，但在制样过程中发现土样存在一定差别，可划分为以下3类，如图 1所示：
(1) 第一类：试样土质较均匀，整体均为淤泥质土。灰色、质软，密度为1.7—1.8g/cm3。
(2) 第二类：试样上下土质不均，夹少量的粉土薄层，密度为1.8—1.9g/cm3。
(3) 第三类：试样土质不均，有互层，夹大量的贝屑。密度为1.8—1.9g/cm3。
其中，埋深为9.5m、13.5m的土样属于第一类土；埋深为6.0m、6.5m、7.0m、7.5m、11.5m的土样属于第二类土；埋深为12.5m的土样属于第三类土。因此，在激振过程中试样变形也表现为3类相似但有一定区别的情况：
(1) 第一类土在激振过程中变形较为均匀，在动应力的作用下呈现一种整体变形的状态，强度较低。
(2) 第二类土在激振过程中的变形主要集中在淤泥段，粉土夹层的变形很小。
(3) 第三类土在激振过程中变形情况与第一类土类似，也比较均匀。但因为土中夹杂着大量的贝屑，故其强度较高，激振结束后变形小于前两类土。
2.2.2 残余应变与振次的关系
残余应变的定义是指动应力卸除后土体的残余变形，试验中为动应力作用后试样的变形量与作用前试样高度的比值 (郁寿松等，1989)。通过动三轴试验，测试了天津滨海软土的残余应变与振次的关系，并利用Monismith等 (1975)提出的指数模型来拟合试样残余应变与振次的关系为：
$$ {\varepsilon _{\rm{P}}} = A{N^b} $$ (1) 式中，εp为残余应变，N为循环荷载振次，A、b为试验参数。
图 2和图 3给出了不同埋深土样在50kPa和60kPa动应力作用下残余应变随振次的变化关系、拟合结果以及95%的置信区间。总体上，可以看出动应力对于残余应变增长的影响较大，随着动应力的增大，残余应变也有着明显的增大。不同埋深的土残余应变随着振次的变化曲线不相同，即使是同一土类反应也不尽相同。在动应力的作用下，除了受到动应力大小的影响之外，试样的残余应变随振次的发展受土质的影响较大，当土质类似时，试样的残余应变大多处于同一水平。
图 2 50kPa动应力下试样残余应变与振次关系曲线Figure 2. Relationship of residual strain and cyclic number under 50kPa
图 3 60kPa动应力下试样残余应变与振次关系曲线Figure 3. Relationship of residual strain and cyclic number under 60kPa由表 2、表 3可以看出，绝大多数拟合结果的相关系数R2都在0.99以上，因此指数模型可以比较好地反映出等压固结情况下试样残余应变与振次之间的关系，用这种简单的模型来进行残余应变与振次关系试验结果的整理分析是可行的。由于软土具有较强的结构性，在振次较少时，试样以弹性变形为主，在达到一定振次之后会出现残余应变的突变，该模型对于突变点的反应较不明显。随着振次的逐渐增多，残余应变的发展十分符合指数增长的规律。
表 2 50kPa动应力下各组拟合成果Table 2. Fitting result of residual strain under 50kPa试样深度/m 拟合参数 相关系数 A b 6.5 0.0014 1.72 0.9941 7.0 0.00018 1.95 0.987 7.5 0.000035 2.47 0.997 9.5 0.00010 2.53 0.9987 11.5 0.00010 2.31 0.995 表 3 60kPa动应力下各组拟合成果Table 3. Fitting result of residual strain under 60kPa试样深度/m 拟合参数 相关系数 A b 6.0 0.0040 2.28 0.9975 7.0 0.00040 1.96 0.9949 7.5 0.0040 1.77 0.9985 11.5 0.00010 2.52 0.995 12.5 0.00060 2.142 0.9988 13.5 0.0020 2.14 0.997 2.3 软土震陷的计算
2.3.1 地震动水平确定
地震荷载是一种随机荷载，具有时域和频域上的非平稳性等特征。在计算中，通常将地震荷载按照一定的规则等效为简谐振动。对于天然土层在地震作用下因土体自重引起的剪应力，采用Seed等 (1966)提出的等效方法：
$$ {\tau _{{\rm{tn}}}} = 0.65{\gamma _{\rm{d}}}{{{a_{\max }}} \over {\rm{g}}}\sum\limits_{i = 1}^n {{\gamma _i}\Delta {h_i}} $$ (2) (2)其中γd=1-0.0133hZ；hz为土层埋深；amax为地面最大加速度，对烈度Ⅶ、Ⅷ、Ⅸ度区可分别取0.125、0.25、0.4g；g为重力加速度；γi为土的天然重度；Δhi为各土层的厚度。则σd土层的动应力可为：
$$ {\sigma _{\rm{d}}} = 2{\tau _{{\rm{tn}}}} $$ (3) 根据上述方法确定勘察场地钻孔土层的竖向应力。动应力简谐振动的振动次数采用Seed等 (1966)提出的与震级相对应的等效破坏振次来确定。历史上在天津地区引发较大震陷的地震为1976年7月28日发生的唐山大地震，震级为7.8级。根据表 4选取20次和30次为震陷计算的等效振次。
表 4 不同震级下等效振动次数Table 4. Equivalent vibration frequency under different magnitude震级 6.5 7 7.5 8 振次 8 12 20 30 2.3.2 残余应变确定
根据对天津滨海地区软土的动三轴试验结果，得到其在7.5级和8级地震作用下的残余应变，如表 5所示。
表 5 天津地区软土残余应变Table 5. Residual strain of soft soil in Tianjin序号 深度 层厚 土类 7.5级地震影响下残余应变 8级地震影响下残余应变 1 4.7 4.7 填土 - - 2 5.5 0.8 填土 - - 3 8.3 2.8 淤粉 0.01 0.2 4 9 0.7 粉粘 0.1 0.3 5 12 3 淤粘 0.1 0.3 6 15.8 3.8 淤粘 0.2 0.35 7 17.5 1.7 粉粘 - - 8 20 2.5 粉粘 - - 2.3.3 土层震前模量
鉴于淤泥质软土的变形表现为明显的非线性，故震前土层的模量采用邓肯模型来表示。根据土的三轴固结排水剪切试验及土层埋藏情况，钻孔土层的震前模量如表 6所示。
表 6 土层震前模量Table 6. Soil modulus before earthquake序号 深度/m 层厚/m 土类 大主应力/kPa 小主应力/kPa 震前模量/kPa 1 4.7 4.7 杂填土 42 28 9571 2 5.5 0.8 素填土 50 34 18103 3 8.3 2.8 淤粉粘 76 50 14632 4 9 0.7 粉粘 81 54 13723 5 12 3 淤粘 103 69 14428 6 15.8 3.8 淤粘 131 88 15186 7 17.5 1.7 粉粘 147 98 16207 8 20 2.5 粉粘 169 113 14232 2.3.4 土层震后模量
震后模量的计算基于模量软化模型 (郁寿松等，1989)，该模型中土的刚度由Ki和Kp两部分组成：
$$ {K_{{\rm{ip}}}} = {1 \over {{1 \over {{K_{\rm{i}}}}} + {1 \over {{K_{\rm{p}}}}}}} $$ (6) 在计算时，用土的初始模量Ei和拟割线模量Ep代替刚度Ki和Kp，得到振动后土的软化模量：
$$ {E_{{\rm{ip}}}} = {1 \over {{1 \over {{E_{\rm{i}}}}} + {1 \over {{E_{\rm{p}}}}}}} $$ (7) $$ {E_{\rm{p}}} = {{{\sigma _{{\rm{as}}}}} \over {{\varepsilon _{\rm{p}}}}} $$ (8) 式中，σas为地震前土体的静偏应力；εp为地震引起的土单元残余应变势，可由动三轴试验获取。根据本方法计算出7.5级和8级地震影响下土层的软化模量如表 7所示。
表 7 土层震后模量Table 7. Soil modulus after earthquake序号 深度/m 层厚/m 土类 震前模量/kPa 7.5级地震影响下割线模量/kPa 8级地震影响下割线模量/kPa 7.5级地震影响下软化模量/kPa 8级地震影响下软化模量/kPa 1 4.7 4.7 杂填土 9571 - 9571 9571 9571 2 5.5 0.8 素填土 18103 - 18103 18103 18103 3 8.3 2.8 淤粉粘 14632 251667 14632 13828 6765 4 9 0.7 粉粘 13723 27033 13723 9103 5439 5 12 3 淤粘 14428 34433 14428 10168 6392 6 15.8 3.8 淤粘 15186 21903 15186 8968 6861 7 17.5 1.7 粉粘 16207 - 16207 16207 16207 8 20 2.5 粉粘 14232 - 14232 14232 14232 需要说明的是，为考虑淤泥质软土的非线性变形特性，本文中震前土层的模量采用邓肯模型来表示，若工程中不具备试验条件，可用压缩模量代替。
2.4 震陷量修正系数
根据上述方法，计算出该钻孔在7.5级和8级地震影响下地面峰值加速度为0.4g时自由地面的沉降量分别为4.9cm和12.9cm，具体计算结果如表 8所示。
表 8 地震影响下自由地表震陷Table 8. Seismic settlement of free ground under earthquake motion序号 深度/m 层厚/m 土类 7.5级地震影响下各层震陷量/m 7.5级地震影响下总沉震陷量/cm 8级地震影响下各层震陷量/m 8级地震影响下总震陷量/m 1 4.7 4.7 填土 - 4.9 - 12.9 2 5.5 0.8 填土 - - 3 8.3 2.8 淤粉 0.001 0.025 4 9 0.7 粉粘 0.003 0.009 5 12 3 淤粘 0.013 0.039 6 15.8 3.8 淤粘 0.032 0.056 7 17.5 1.7 粉粘 - - 8 20 2.5 粉粘 - - 对于不考虑模量软化的自由场地沉降值，可得到该钻孔处自由场地的沉降量20cm，具体计算结果如表 9所示。
表 9 不考虑模量软化的自由场地沉降量Table 9. Free ground settlement without consideration of soften model序号 深度/m 层厚/m 土类 竖向应力/kPa 动应力/kPa 震前模量/kPa 应变值/% 沉降值/cm 1 4.7 4.7 杂填土 42.3 20.6 9571 0.007 20 2 5.5 0.8 素填土 50.3 24.2 18103 0.004 3 8.3 2.8 淤粉粘 75.5 34.9 14632 0.008 4 9 0.7 粉粘 81.1 37.1 13723 0.009 5 12 3 淤粘 103.3 45.1 14428 0.010 6 15.8 3.8 淤粘 131.42 54.0 15186 0.012 7 17.5 1.7 粉粘 146.72 58.5 16207 0.013 8 20 2.5 粉粘 169.22 64.6 14232 0.016 根据震陷量修正系数定义可得出，在7.5级地震作用下的震陷量大约是不考虑模量软化时的四分之一，在8级地震作用下的震陷量大约是不考虑模量软化时的二分之一。采用相同方法计算得到了地面峰值加速度为0.25g时自由场地震陷量修正系数，经估算天津滨海地区自由场地软土震陷量修正系数如表 10所示。
表 10 天津滨海地区自由场地软土震陷量修正系数Table 10. The settlement correction coefficient of free soft soil site in Tianjin地面峰值加速度/g 震级/M 修正系数C 0.25 7.5 0.2 8 0.6 0.4 7.5 0.25 8 0.5 需要说明的是，由于在求残余应变势时，动应力加载的振次采用Seed等 (1966)提出的与震级相对应的等效破坏振次来确定，因此上述表格给出的是震陷修正系数与震级的关系。我国抗震设防是以烈度形式给出，因此本表要结合该地区地震安全性评价中潜在地震危险震级和设防烈度联合使用，通过潜在地震危险震级分析获得该地区可能遭受的震级，也可通过历史地震的影响确定震级，设防烈度为Ⅷ度的地区选用0.25g，Ⅸ度区选用0.4g，从而确定震陷修正系数。
3. 震害实例对比
1976年唐山大地震中，天津滨海地区的地基普遍产生了不同程度的震陷，其中天津新港地区望海楼住宅区和建港村住宅区的震陷现象尤为明显。这些震陷事例对于回答软土遭受振动能否构成震陷，以及确定软粘土地基形成的条件和检验估算地基震陷的方法，提供了丰富而宝贵的资料。塘沽望海楼建筑群主要建筑为3、4层砖混结构住宅楼，基础形式为片筏，埋深0.6m，基底附加压力44kPa，观测沉降量17—22cm (刘恢先，1986)。1978年唐山7.8级地震中，天津塘沽地区烈度达到Ⅷ度，地面加速度峰值0.25g。
根据不考虑模量软化的简化算法，计算出的震陷量15.2cm。根据上文提出的经验系数，该地区修正系数取0.6，则该地区自由地表震陷量约为9.1cm。由于本计算只考虑自由场地的震陷量，故对望海楼地区地表震陷量的估算小于唐山地震中的实际震陷量。通过对比结果分析如下：
(1) 根据震害观测资料，塘沽新港地区的自然地面的震陷值约为8cm，北塘食品库震后桩头出露地表约12cm，与本经验系数估算结果吻合较好。
(2) 根据震害调查和观测结果可知实际地震作用下地面的震陷是不均匀的，本方法计算所得的是一个钻孔震陷结果，因此本估算方法只能代表在一定情况下的平均值。
(3) 本计算方法避免了在残余应变确定过程中所需进行的大量试验以及复杂的有限元计算，只需简单计算就可以对场地的震陷量进行初步的估算，方便快捷。
(4) 本方法存在比较多的假设条件，所给出经验系数仅是在特定情况下的自由地面震陷估计值。今后应对不同的场地类型、建筑物结构和基础形式进行分类，系统地给出震陷量估算的参考公式，形成科学的震陷量估算评价体系。
4. 结论及讨论
本文在前人工作的基础上，通过对天津滨海新区典型软土的静、动力学试验，分析了软土残余应变的变化规律，提出了结合地区经验系数的软土震陷估算方法，并利用1976年唐山7.8级地震中天津塘沽地区软土震陷观测资料对计算方法进行了验证。结果表明，本方法对于自由地表震陷量的估算是可行的，对于中小工程中软土震陷量的估算具有一定的参考价值。本文提出的震陷量修正系数省去了动三轴试验繁琐的操作，减少了相关的工作量，简化了软土震陷震陷量的计算方法，易被一般工程技术人员掌握。
本文基于软土的静、动力学试验，对现有软土计算方法中存在的难点进行了初步的研究。由于作者水平和时间的限制，还有大量问题有待于进一步的研究。本文的震陷试验是在等压固结的条件下进行的，对于残余应变的变化规律仅做了定性的研究，因此有必要再通过大量试验，确定出残余应变随动力条件变化的经验公式以及其参数与土样物理条件之间的关系；试验中动荷载是依据Seed等 (1966)提出的等效方法简化而来，然而研究表明动荷载的加载方式对土的震陷量有较大的影响，有必要提出一种更加适用的等效方法；动三轴试验确定的仅仅是土单元的相关参数，而尺度效应对于试验结果有较大的影响，因此应进行一些模型试验，如离心机振动台试验等，来更好的模拟土的原位应力状态；本文仅是用经验系数的方法对天津软土震陷计算进行了初步研究，未能真实反映出其物理机制，下一步研究应从其震陷机理入手，提出更能反映真实情况的震陷计算方法；本文给出的震陷修正系数表是基于震陷修正系数与震级的关系，不能很好地与我国的抗震设防烈度相对应，应用时需要确定影响震级，使用时仍有不便，应通过进一步研究给出与我国抗震设防烈度相对应的震陷修正系数表。
-
图 1 研究区位置及周边主要活动断裂分布
Figure 1. The distribution map of the main active faults, earthquakes and the location of the study area
图 2 研究区地形地貌
Figure 2. The topography map showing the rives, active faults and earthquake in the study area
-
常直杨, 王建, 白世彪等, 2015. 面积高程积分值计算方法的比较. 干旱区资源与环境, 29（3）: 171—175Chang Z. Y. , Wang J. , Bai S. B. , et al. , 2015. Comparison of hypsometric integral methods. Journal of Arid Land Resources and Environment, 29(3): 171—175. (in Chinese) 陈长云, 任金卫, 孟国杰等, 2013. 巴颜喀拉块体东部活动块体的划分、形变特征及构造意义. 地球物理学报, 56（12）: 4125—4141Chen C. Y. , Ren J. W. , Meng G. J. , et al. , 2013. Division, deformation and tectonic implication of active blocks in the eastern segment of Bayan Har block. Chinese Journal of Geophysics, 56(12): 4125—4141. (in Chinese) 陈彦杰, 郑光佑, 宋国城, 2005. 面积尺度与空间分布对流域面积高度积分及其地质意义的影响. 地理学报, （39）: 53—69Chen Y. J. , Zheng G. Y. , Song G. C. , 2005. Influence of area and space dependence for hypsometric integral and its geological implications. Journal of Geographical Science, (39): 53—69. (in Chinese) 高青, 李奋生, 邵崇建, 2018. 汶川地震驱动的剥蚀作用及对龙门山地貌演化的影响. 成都理工大学学报（自然科学版）, 45（3）: 386—392Gao Q. , Li F. S. , Shao C. J. , 2018. The erosion driven by Wenchuan earthquake and its impact on the geomorphic evolution of Longmenshan, Sichuan, China. Journal of Chengdu University of Technology (Science & Technology Edition), 45(3): 386—392. (in Chinese) 韩海辉, 高婷, 易欢等, 2012. 基于变点分析法提取地势起伏度——以青藏高原为例. 地理科学, 32（1）: 101—104Han H. H. , Gao T. , Yi H. , et al. , 2012. Extraction of relief amplitude based on Change Point Method——A case study on the Tibetan Plateau. Scientia Geographica Sinica, 32(1): 101—104. (in Chinese) 何建军, 任俊杰, 丁锐等, 2016. 青藏高原东缘龙日坝断裂带南段晚第四纪活动及其构造意义. 震灾防御技术, 11（4）: 707—721He J. J. , Ren J. J. , Ding R. , et al. , 2016. Late quaternary activity of the Southern Segment of Longriba fault zone in Eastern Tibet and its tectonic implications. Technology for Earthquake Disaster Prevention, 11(4): 707—721. (in Chinese) 洪艳, 赵银兵, 王运生等, 2019. 利用面积高程积分方法研究龙门山断裂带地貌隆升特征. 科学技术与工程, 19（30）: 43—51Hong Y. , Zhao Y. B. , Wang Y. S. , et al. , 2019. Study on the geomorphic uplift characteristics of the Longmenshan fault zone using hypsometric integral method. Science Technology and Engineering, 19(30): 43—51. (in Chinese) 贾营营, 付碧宏, 王岩等, 2010. 青藏高原东缘龙门山断裂带晚新生代构造地貌生长及水系响应. 第四纪研究, 30（4）: 825—836Jia Y. Y. , Fu B. H. , Wang Y. , et al. , 2010. Late Cenozoic tectono-geomorphic growth and drainage response in the Longmen Shan fault zone, east margin of Tibet. Quaternary Sciences, 30(4): 825—836. (in Chinese) 李奋生, 赵国华, 李勇等, 2015. 龙门山地区水系发育特征及其对青藏高原东缘隆升的指示. 地质论评, 61（2）: 345—355Li F. S. , Zhao G. H. , Li Y. , et al. , 2015. The characteristics of drainage development in Longmen Mountains Area and its indication to the uplift of eastern margin of Qinghai-Xizang (Tibet) Plateau. Geological Review, 61(2): 345—355. (in Chinese) 李明, 张世民, 姜大伟等, 2019. 龙门山中北段流域地貌特征及其构造意义. 震灾防御技术, 14（3）: 640—651Li M. , Zhang S. M. , Jiang D. W. , et al. , 2019. Geomorphic features and tectonic significance of the middle and North Section of Longmenshan. Technology for Earthquake Disaster Prevention, 14(3): 640—651. (in Chinese) 李伟, 谢超, 程宏宾等, 2021. 利用水系方位角和GPS数据研究龙门山后山断裂运动模式. 四川地震, （2）: 7—10Li W. , Xie C. , Cheng H. B. , et al. , 2021. Study on the fault movement model of the backmountain of Longmen Mountain based on the azimuth angle and GPS data. Earthquake Research in Sichuan, (2): 7—10. (in Chinese) 李勇, 周荣军, 赵国华等, 2015. 汶川地震驱动的隆升、剥蚀作用与龙门山地貌生长——以映秀红椿沟为例. 成都理工大学学报（自然科学版）, 42（1）: 5—17Li Y. , Zhou R. J. , Zhao G. H. , et al. , 2015. Uplift and erosion driven by Wenchuan earthquake and their effects on geomorphic growth of Longmen Mountains——A case study of Hongchun gully in Yingxiu, China. Journal of Chengdu University of Technology (Science & Technology Edition), 42(1): 5—17. (in Chinese). 梁明剑, 周荣军, 闫亮等, 2014. 青海达日断裂中段构造活动与地貌发育的响应关系探讨. 地震地质, 36（1）: 28—38Liang M. J. , Zhou R. J. , Yan L. , et al. , 2014. The relationships between neotectonic activity of the middle segment of Dari Fault and its geomorphological response, Qinghai Province, China. Seismology and Geology, 36(1): 28—38. (in Chinese) 梁明剑, 杨耀, 杜方等, 2020. 青海达日断裂中段晚第四纪活动性与1947年M7¾地震地表破裂带再研究. 地震地质, 42（3）: 703—714Liang M. J. , Yang Y. , Du F. , et al. , 2020. Late Quaternary activity of the central segment of the Dari Fault and restudy of the surface rupture zone of the 1974 M7¾ Dari earthquake, Qinghai Province. Seismology and Geology, 42(3): 703—714. (in Chinese) 梁欧博, 任俊杰, 吕延武, 2018. 涪江流域河流地貌特征对虎牙断裂带活动性的响应. 地震地质, 40（1）: 42—56Liang O. B. , Ren J. J. , Lü Y. W. , 2018. The response of fluvial geomorphologic characteristics of the Fujiang drainge basin to activity of the Huya fault zone. Seismology and Geology, 40(1): 42—56. (in Chinese) 刘静, 曾令森, 丁林等, 2009. 青藏高原东南缘构造地貌、活动构造和下地壳流动假说. 地质科学, 44（4）: 1227—1255Liu J. , Zeng L. S. , Ding L. , et al. , 2009. Tectonic geomorphology, active tectonics and lower crustal channel flow hypothesis of the southeastern Tibetan Plateau. Chinese Journal of Geology, 44(4): 1227—1255. (in Chinese) 苏琦, 2015. 青藏高原东北缘典型流域地貌参数分析与构造变形探讨. 兰州: 中国地震局兰州地震研究所.Su Q. , 2015. Study on typical drainage basins along northeastern Tibetan Plateau and its tectonic deformation. Lanzhou: China Earthquake Administration Lanzhou Institute of Seismology. (in Chinese) 孙东, 王道永, 吴德超等, 2010. 马尔康巴拉水电站近场区主要断裂活动性及对工程的影响. 工程地质学报, 18（6）: 940—949Sun D. , Wang D. Y. , Wu D. C. , et al. , 2010. Activity and effect of main faults in near field of Bala Hydropower station in Maerkang. Journal of Engineering Geology, 18(6): 940—949. (in Chinese) 熊仁伟, 任金卫, 张军龙等, 2010. 玛多—甘德断裂甘德段晚第四纪活动特征. 地震, 30（4）: 65—73Xiong R. W. , Ren J. W. , Zhang J. L. , et al. , 2010. Late Quaternary active characteristics of the Gande Segment in the Maduo-Gande Fault zone. Earthquake, 30(4): 65—73. (in Chinese) 徐锡伟, 闻学泽, 陈桂华等, 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义. 中国科学 D辑: 地球科学, 38（5）: 529—542.Xu X. W. , Wen X. Z. , Chen G. H. , et al. , 2008. Discovery of the Longriba fault zone in eastern Bayan Har block, China and its tectonic implication. Science in China Series D: Earth Sciences, 51(9): 1209—1223. 闫亮, 李勇, 赵国华等, 2018. 青藏高原东缘龙门山构造带晚第四纪构造隆升作用的河流地貌响应. 第四纪研究, 38（1）: 232—246Yan L. , Li Y. , Zhao G. H. , et al. , 2018. The uplift process and the geomorphological features of the rivers system in the Longmen Shan since Late Quaternary. Quaternary Sciences, 38(1): 232—246. (in Chinese) 詹艳, 梁明剑, 孙翔宇等, 2021.2021年5月22日青海玛多MS7.4地震深部环境及发震构造模式. 地球物理学报, 64（7）: 2232—2252Zhan Y. , Liang M. J. , Sun X. Y. , et al. , 2021. Deep structure and seismogenic pattern of the 2021.5. 22 Madoi (Qinghai) MS7.4 earthquake.Chinese Journal of Geophysics,64(7):2232—2252.(in Chinese) 张会平, 2006. 青藏高原东缘、东北缘典型地区晚新生代地貌过程研究[D]. 北京: 中国地质大学.Zhang H. P. , 2006. Study on late Cenozoic geomorphic processes of typical regions along the eastern and northeastern Tibetan margins. Beijing: China University of Geosciences. (in Chinese) 张培震, 闻学泽, 徐锡伟等, 2009.2008年汶川8.0级特大地震孕育和发生的多单元组合模式. 科学通报, 54（7）: 944—953 doi: 10.1360/csb2009-54-7-944Zhang P. Z. , Wen X. Z. , Xu X. W. , et al. , 2009. Tectonic model of the great Wenchuan earthquake of May 12, 2008, Sichuan, China. Chinese Science Bulletin, 54(7): 944—953. (in Chinese) doi: 10.1360/csb2009-54-7-944 张岳桥, 杨农, 施炜等, 2008. 青藏高原东缘新构造及其对汶川地震的控制作用. 地质学报, 82（12）: 1668—1678Zhang Y. Q. , Yang N. , Shi W. , et al. , 2008. Neotectonics of eastern Tibet and its control on the Wenchuan Earthquake. Acta Geologica sinica, 82(12): 1668—1678. (in Chinese) 郑立龙, 2015. 汶川地震后绵远河流域地表隆升与侵蚀过程研究[D]. 成都: 成都理工大学.Zheng L. L. , 2015. The research of surface uplift and erosion process of Mianyuan river drainage basin after the Wenchuan Earthquake. Chengdu: Chengdu University of Technology. (in Chinese) 周文英, 梁明剑, 左洪等, 2022. 濑溪河流域地貌对泸县MS6.0地震发震构造变形的响应及其指示作用. 地震学报, 44（2）: 250—259Zhou W. Y. , Liang M. J. , Zuo H. , et al. , 2022. Response of the topography of Laixihe drainage to the structural deformation induced by MS6.0 Luxian earthquake and its indication. Acta Seismologica Sinica, 44(2): 250—259. (in Chinese) Chen Y. C. , Sung Q. C. , Cheng K. Y. , 2003. Along-strike variations of morphotectonic features in the Western Foothills of Taiwan: tectonic implications based on stream-gradient and hypsometric analysis. Geomorphology, 56(1—2): 109—137. doi: 10.1016/S0169-555X(03)00059-X Clark M. K. , Royden L. H. , 2000. Topographic ooze: building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703—706. doi: 10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2 Fu B. H. , Shi P. L. , Guo H. D. , et al. , 2011. Surface deformation related to the 2008 Wenchuan earthquake, and mountain building of the Longmen Shan, eastern Tibetan Plateau. Journal of Asian Earth Sciences, 40(4): 805—824. doi: 10.1016/j.jseaes.2010.11.011 Howard A. D. , Kerby G. , 1983. Channel changes in Badlands. GSA Bulletin, 94(6): 739—752. doi: 10.1130/0016-7606(1983)94<739:CCIB>2.0.CO;2 Hubbard J. , Shaw J. H. , 2009. Uplift of the Longmen Shan and Tibetan Plateau, and the 2008 Wenchuan (M=7.9) earthquake. Nature, 458(7235): 194—197. doi: 10.1038/nature07837 Kirby E. , Whipple K. X. , Tang W. Q. , et al. , 2003. Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from bedrock channel longitudinal profiles. Journal of Geophysical Research: Solid Earth, 108(B4): 2217. Kirby E., Ouimet W., 2011. Tectonic geomorphology along the eastern margin of Tibet: Insights into the pattern and processes of active deformation adjacent to the Sichuan Basin. Geological Society, London, Special Publications, 353(1): 165—188. Liu-Zeng J., Wen L., Oskin M., et al., 2011. Focused modern denudation of the Longmen Shan margin, eastern Tibetan Plateau. Geochemistry, Geophysics, Geosystems, 12(11): Q11007. Meng Q. R. , Hu J. M. , Wang E. , et al. , 2006. Late Cenozoic denudation by large-magnitude landslides in the eastern edge of Tibetan Plateau. Earth and Planetary Science Letters, 243(1—2): 252—267. doi: 10.1016/j.jpgl.2005.12.008 Parker R. N. , Densmore A. L. , Rosser N. J. , et al. , 2011. Mass wasting triggered by the 2008 Wenchuan earthquake is greater than orogenic growth. Nature Geoscience, 4(7): 449—452. doi: 10.1038/ngeo1154 Pike R. J. , Wilson S. E. , 1971. Elevation-relief ratio, hypsometric integral, and geomorphic area-altitude analysis. GSA Bulletin, 82(4): 1079—1084. doi: 10.1130/0016-7606(1971)82[1079:ERHIAG]2.0.CO;2 Ren J. J. , Xu X. W. , Yeats R. S. , et al. , 2013. Latest Quaternary paleoseismology and slip rates of the Longriba fault zone, eastern Tibet: Implications for fault behavior and strain partitioning. Tectonics, 32(2): 216—238. doi: 10.1002/tect.20029 Royden L. H. , Burchfiel B. C. , King R. W. , et al. , 1997. Surface deformation and lower crustal flow in Eastern Tibet. Science, 276(5313): 788—790. doi: 10.1126/science.276.5313.788 Schwanghart W. , Scherler D. , 2014. TopoToolbox 2 – MATLAB-based software for topographic analysis and modeling in Earth surface sciences. Earth Surface Dynamics, 2(1): 1—7. doi: 10.5194/esurf-2-1-2014 Snyder N. P. , Whipple K. X. , Tucker G. E. , et al. , 2000. Landscape response to tectonic forcing: Digital elevation model analysis of stream profiles in the Mendocino triple junction region, northern California. Geological Society of America Bulletin, 112(8): 1250—1263. doi: 10.1130/0016-7606(2000)112<1250:LRTTFD>2.0.CO;2 Strahler A. N. , 1952. Hypsometric (area-altitude) analysis of erosional topography. GSA Bulletin, 63(11): 1117—1142. doi: 10.1130/0016-7606(1952)63[1117:HAAOET]2.0.CO;2 Wang Y. Z. , Zhang H. P. , Zheng D. W. , et al. , 2017. Coupling slope–area analysis, integral approach and statistic tests to steady-state bedrock river profile analysis. Earth Surface Dynamics, 5(1): 145—160 doi: 10.5194/esurf-5-145-2017 Whipple K., Wobus C. W., Crosby B. T., et al., 2007. New tools for quantitative geomorphology: extraction and interpretation of stream profiles from digital topographic data [EB/OL]. (2007-01-01)[2023-11-10].https://www.researchgate.net/profile/Cameron-Wobus/publication/268056101_New_Tools_for_Quantitative_Geomorphology_Extraction_and_Interpretation_of_Stream_Profiles_from_Digital_Topographic_Data/links/54f9bf970cf29a9fbd7c51d4/New-Tools-for-Quantitative-Geomorphology-Extraction-and-Interpretation-of-Stream-Profiles-from-Digital-Topographic-Data.pdf. -
下载:
下载:
点击查看大图
计量
- 文章访问数: 142
- HTML全文浏览量: 43
- PDF下载量: 18
- 被引次数: 0
