Nonlinear Seismic Response Analysis of Two-dimensional Depositional Basins
-
摘要: 沉积盆地对地震动具有明显的放大效应。利用ABAQUS软件及二次开发平台，编写了基于Davidenkov骨架曲线的VUMAT子程序，对饱和沉积盆地进行了非线性反应模拟。建立不同基岩面坡度的二维T型沉积盆地，输入4条不同频谱特性及幅值的近断层地震波，研究不同盆地基岩面坡度、饱和介质特性、地震波强度等因素对沉积盆地非线性反应规律的影响。研究结果表明，地震波作用可造成饱和沉积盆地发生局部液化，地震波强度越大，液化区分布范围越广，地表处出现了永久塑性变形；地震波强度的增大导致地表土层发生液化的时间提前，且下层土的液化会减弱地震波的传播速度，导致上层地表土发生液化的时间延后；盆地内地表最不利位置随着基岩面坡度的增大逐渐向中心处移动，不同工况下速度峰值放大系数为1.18～3.33。
-
关键词:
- Davidenkov本构模型 /
- VUMAT /
- 二维非线性分析 /
- 近断层沉积盆地 /
- 场地反应
Abstract: Near fault seismic motion can lead to significant seismic disasters in the near field area, particularly in regions with sedimentary terrain, exacerbating the potential disaster. In this paper, ABAQUS software and secondary development platform are employed to compile VUMAT subroutine based on Davidenkov skeleton curve, to conduct nonlinear simulation of saturated sedimentary basin, and to establish two-dimensional T-shaped Sedimentary basin with different bedrock surface slopes. Four near fault Seismic wave with different spectral characteristics are selected to input the site with different amplitudes, and the nonlinear response of the site under the action of different basin slopes and seismic wave intensities is studied. It is found that seismic wave can cause liquefaction in saturated sedimentary basin. The greater the strength of seismic wave is, the wider the distribution range of liquefaction area is, and permanent plastic deformation occurs on the surface; The increase of seismic wave strength can lead to the advance of liquefaction time of surface soil layer, and the liquefaction of subsoil will weaken the propagation strength of Seismic wave, resulting in the delay of liquefaction time of upper surface soil layer.As the bedrock surface slope increases, the most vulnerable location on the basin's surface shifts toward the center. Furthermore, the amplification coefficient of peak velocity ranges from 1.18 to 3.33 under various working conditions. -
引言
近年来，针对地震活动、断裂构造进行的地球化学探测逐渐兴起并发展（郑海刚等，2016；王喜龙等，2017），氢气作为断裂带主要气体成分之一，被认为与地震构造活动密切相关（杜乐天，2005；高小其等，2017；向阳等，2018）。众所周知，自然界氢气来源主要有以下几点：（1）断层错动：沿地壳下层塑性岩石（蛇纹岩）或硅酸岩，断层新鲜裂隙、断层错动过程中产生高浓度氢气（陈辉，1996；夏林圻等，2008）；（2）应力活动：主要来源为地壳深部岩石孔隙、裂隙中封存的氢气（孙小龙等，2017）；（3）生物化学作用也可产生氢气，但仅有微量，且存在于浅层地表（夏林圻等，2008；车用太等，2002）。地震发生前氢气浓度异常已被众多震例证实，且异常多出现在震前短临阶段（林元武等，1994；张彬等，2018）。国内外已有众多震例证明氢气在震前、震后均发生明显变化，地下氢气是灵敏反映地下应力应变的地球化学成分（高清武，1992；张培仁等，1993；车用太等，2002），断裂带氢气浓度与地震活动释放的能量呈正相关（Ryuichi，1984）。近年来，很多学者利用氢气研究地震、断层，取得了良好的效果（李源等，2018；康健等，2019）。氢气浓度观测是值得广泛推广的地震监测方法之一，也被认为是具有广阔前景的地震监测预报方法（车用太等，2002；刘耀炜等，2006；范雪芳等，2016）。
依兰-伊通断裂是东北地区规模最大的断裂，近年来在断裂通河段发现了1次古地震事件，探槽揭示距今1700年前发生过7级以上地震事件（闵伟等，2011），这一发现彻底颠覆了人们对该断裂“全新世以来不活动”的认识。本文利用地球化学探测方法研究依兰-伊通断裂（图 1），了解依兰-伊通断裂断层氢气时空变化特征，捕捉短期地震构造活动变化信息，建立氢气背景值及氢气与地震活动的耦合关系。
图 1 依兰-伊通断裂北段构造图Figure 1. Structural map of the northern section of Yilan-Yitong fault1. 构造背景
依兰-伊通断裂总体呈NNE-NE向展布，中国境内总长约1100km，是郯庐断裂带渤海以北区段的主要组成部分（徐嘉炜等，1992），也是划分中蒙、中朝和燕山3个活动地块的重要边界断裂（张培震等，2003）。该断裂基本由近于平行的东、西2条分支断裂（东支F4-2和西支F4-1）组成，构成宽10—20km的地堑，部分区段呈北东向沿松花江分布。
依兰-伊通断裂带在黑龙江境内一般分为北段（方正-萝北段）、中段（尚志-延寿-方正段）和南段（伊通-舒兰-尚志段）。依兰-伊通断裂北段地震活动明显强于中段和南段，沿着北段存在整个黑龙江省内最密集的小震活动，显示了较强烈的当前地壳活动特征。就北段而言，地震活动分布不均，萝北地区地震活动强于北段其他区域，该地区小震活跃，历史上发生过多次中强地震，如1963年6月22日M5.8地震。研究发现西支断裂较东支断裂活动性明显，西支通河段发现存在高2m左右的断层陡坎等地貌线性影像，延伸长约70km，逆走滑性质，并在祥顺乡南楼村探槽剖面发现全新世活动，并推测约在1700年前发生过强烈古地震事件（闵伟等，2011；疏鹏等，2014）。
2. 测线布设与测量方法
选择在依兰-伊通断裂北段西支断裂上进行研究，由方正至萝北选出8条剖面垂直断裂进行观测（图 1），由南至北依次为方正西测线、祥顺乡测线、松花江农场测线、香兰农场测线、汤原测线、鹤岗南测线、宝泉岭测线和延军农场4队测线。观测剖面垂直断裂，每条测线等距布设调整间距15—20调整间距个测点，点距一般为调整间距10—50m，具体间距以实际情况为准，可根据地形地貌和断层类型调整。如果未发现断层具体位置，则根据收集到的资料分析断层可能的位置，用于确定最佳观测位置。
观测仪器为便携式氢气分析仪（ATG-300H型），标定后空气检测值一般为0.5×10-6。测量时需在土壤中打观测孔，测孔直径约3cm，孔深约80cm。首先将取样器（麻花钻）吸气端插入孔中并封住气孔，然后用软管将取样器排气端与仪器相连，并启动仪器开始观测，读取最高值。为减少非构造因素的影响，选取土质条件相同的场地进行观测。
3. 观测结果
3.1 数据
经多期观测，得到断裂8条测线氢气浓度背景值。背景值采用各测点浓度值的几何平均值，这种方法能弱化极端数值对整体数据的影响，从而更真实地反映数据的整体水平，一般可体现研究区域客观氢气浓度背景。计算公式如下：
$$ G=\sqrt[n]{{x}_{1}·{x}_{2}·\cdots ·{x}_{n}}$$ (1) 式中，G为观测点所有样本几何平均值；${x_n}$为观测点第n个样本；n为观测点样本个数。
观测时间段为2017年6月至2018年9月，每条测线观测2—5期。其中，方正西测线、祥顺乡测线、汤原测线完成2期测量；香兰农场测线、宝泉岭测线、延军农场4队测线完成3期测量；鹤岗南测线完成4期测量；松花江农场测线完成5期测量。整体数据如表 1所示。
表 1 测线断层氢气统计结果Table 1. Hydrogen statistical results of line faults测线序号 测线名称 观测日期/年-月-日 测点个数 背景值/ppm 1 方正西 2017-09-12 15 62 2018-08-19 15 51 2 祥顺乡 2017-09-12 15 26 2018-08-18 15 260 3 松花江农场 2017-09-13 15 267 2017-11-02 15 106 2018-08-18 15 602 2018-09-13 15 245 2018-09-19 15 162 4 香兰农场 2018-01-19 30 47 2018-04-02 30 178 2018-08-17 15 359 5 汤原 2017-07-18 15 98 2018-08-15 15 168 6 鹤岗南 2017-07-19 15 744 2017-11-03 15 33 2018-08-12 15 214 2018-09-14 15 158 7 宝泉岭 2017-09-14 15 360 2018-08-13 15 569 2018-09-16 15 280 8 延军农场4队 2017-06-12 15 92 2017-07-20 15 168 2018-08-13 15 434 3.2 数据分析
将各条测线氢气浓度背景值按月统计，得到2017—2018年部分时段的背景值，如图 2所示。由图 2可知，香兰测线和延军农场4队测线氢气背景值由春季到夏季逐渐升高；松花江农场测线、鹤岗南测线和宝泉岭农场测线氢气背景值由夏季到秋季逐渐降低，9月氢气背景值较8月大幅下降。
由图 2可知，春季至夏季为氢气浓度上升段，夏季至秋季为氢气浓度下降段，每年各条测线最高值基本出现在8月前后。杜建国等（1998）发现土壤氡气浓度值存在冬夏高、春秋低的年变规律，并认为断层气在土壤中表现出的夏季高主要受土壤湿度的影响（王博等，2010；孙小龙等，2017），冬季高主要受冻土层的影响（康健等，2019）。土壤中的氡气与氢气源于地壳深部，运移路径与方式相似。图 2中数据虽不完整，但仍可显示氢气浓度具有夏季高、春秋低的特点。
将2018年观测数据与2017年同期数据进行比较，由于条件限制，使每年的观测时间难以保持一致，相对而言，2017、2018年7—9月（夏季）观测最集中，因此以7—9月为基础，最终得到夏季对比数据，如图 3所示。祥顺乡测线、鹤岗南测线、延军农场4队测线氢气浓度同期变化较显著，这可能与构造活动有关。其余5条测线氢气浓度同期变化不大，相关系数达0.87，可见背景值变化较稳定，如松花江农场测线2017年9月13日测得的背景值为267ppm，于2018年9月13日再次观测，测得背景值为245ppm，表明该地区2017、2018年氢气浓度非常稳定。同期变化较大的测线中，鹤岗南测线背景值2018年较2017年明显降低，而祥顺乡测线、延军农场4队测线背景值2018年较2017年明显升高。将依兰-伊通断裂北段观测到的每条剖面进行整体比较（图 3），可知鹤岗南测线2017年背景值最高，为742ppm，祥顺乡测线最低，为26ppm，整条断裂年平均值为279ppm；延军农场4队测线2018年背景值最高，为434ppm，汤原测线最低，为168ppm，整条断裂年平均值为268ppm。综上所述，2017、2018年平均值近乎持平，但从变化幅度上看，2018年较2017年更趋于稳定状态。
图 3 依兰-伊通断裂北段各条测线断层氢气统计结果Figure 3. Hydrogen statistical results of line faults of the northern section of Yilan-Yitong fault2018年8月在依兰-伊通断裂北段观测了完整的一期，所有测线均在同一时段内完成观测。图 4中红色数据为2018年8月各条测线氢气浓度分布特征，可知在本文研究的断裂区域内，松花江农场测线、香兰农场测线、宝泉岭测线和延军农场4队测线氢气浓度较高；蓝色数据为2017、2018年氢气浓度平均背景值，可知其变化趋势与氢气浓度变化趋势相似，说明2017、2018年氢气浓度变化较稳定，也可知夏季浓度值普遍高于年平均值。
图 4 2017—2018年依兰-伊通断裂北段氢气浓度特征划分图Figure 4. Characteristic partition of hydrogen concentration of the northern section of Yilan-Yitong fault from 2017 to 2018同一时刻的横向比较更能看出断裂氢气浓度分布特征，因此将整条断裂数据进行横向比较，根据2018年8月观测结果，按浓度值将断裂简单划分为方正-依兰段（平均氢气浓度＜250ppm）、依兰-汤原段（平均氢气浓度＞450ppm）、汤原-鹤岗段（平均氢气浓度＜250ppm）、鹤岗-萝北段（平均氢气浓度＞450ppm），如图 4所示，可见断层氢气在断裂不同地段释放量不同。
4. 机理分析
对断裂上测线进行环比分析，可知氢气浓度呈夏季高、春秋低的特点。土壤氢气浓度具有年变特征，造成这种情况的原因是温度、湿度变化（刘学领等，2011；康健等，2019）。研究表明微生物活动可产生氢气，而温度、湿度影响微生物活跃程度，温度低于10℃时大量微生物停止活动或死亡，温度为20—30℃时微生物最活跃（樊耀亭等，2002；Des Marais，2000；Kasting等，2006；梅冥相，2014；车用太等，2015）。湿度增加会减小土壤孔隙度，孔隙度变小会造成土壤中气体浓度增加（孙小龙等，2017）。因此，夏季温度、湿度最适合微生物生长，而黑龙江地区春秋季温度较低，且干燥少雨，微生物基本停止活动。因此土壤氢气浓度夏季高、春秋低的年变特征主要与微生物活动有关。
对断裂上测线进行同期分析，可知测线2017、2018年平均值近乎持平。特别是松花江农场测线2017、2018年观测日期一致，观测数值也近乎相等，说明该地区断层氢气在2017、2018年输出稳定。个别地区前后2年同期相差较大，分析原因如下：（1）2017、2018年降雨量不同，土壤湿度不同，造成土壤内孔隙比不同，影响氢气浓度变化，且湿度不同还影响微生物的活跃程度；（2）地下地震构造活动存在变化。
对断裂同一时段进行横向比较，可知每个区域氢气浓度变化很大，按浓度特征可简单划分为4个区域。将分区结果结合区内地震活动特征进行分析（图 5），可知地震活动有明显的分段性。图 5（a）显示氢气浓度高低与历史地震活动强弱的对应关系并不明显，图 5（b）显示氢气浓度高低与当前地震活动存在对应关系，说明氢气浓度高低与当前地震活动存在耦合关系（车用太等，2002；李晨桦等，2016）。特别是依兰-汤原段香兰农场测线附近，2018年1月连续发生3.2级和3.8级地震，与地震对应的分区氢气浓度偏高，说明断层气浓度反映当前地震构造活跃情况，与当前地震活动有直接关系，与历史地震甚至地震多发区无关。
图 5 依兰-伊通断裂北段地震分布图Figure 5. Earthquake distribution map of the northern section of Yilan-Yitong fault5. 结论与讨论
通过对依兰-伊通断裂北段氢气流动进行观测，得出以下结论：
（1）从时间纵向比较，依兰-伊通断裂北段2018年背景值与2017年较接近，受地震构造活动等因素影响局部升高，个别有下降；氢气浓度存在年变特征，主要受温度、湿度和微生物活动的影响，一年之中夏季氢气浓度最高，春秋季氢气浓度最低。
（2）从空间横向比较，依兰-伊通断裂北段呈现出分段性特征，氢气浓度在该段断裂上呈现出北高南低的特征，这与当前地震活动特征相符，高值地区对应地震活动较强的地区，低值地区对应地震活动较弱的地区，断层气浓度与当前地震构造活动有很好的耦合关系。
（3）本文没有连续观测数据和显著震例，无法捕捉断层氢气浓度在震前、震后的变化特征，但不难看出构造地球化学方法对当前地震构造活动映震效果明显。如果增加观测手段，多种断层气同步观测，相互印证，将极大提高观测结果的可信度。对重点地区进行连续观测可提高观测精度并及时掌握断层气动态变化，有利于监视地震构造活动。相信随着科技的发展，构造地球化学观测将在分析预报、短临跟踪、断层研究等领域发挥重要作用。
-
图 7 试验与仿真孔压时程曲线对比
Figure 7. Comparison of experimental and simulated pore pressure time-histories
图 10 基岩面坡度为30°时孔压比分布
Figure 10. Distribution of pore pressure ratio with the slope of the bedrock surface of 30°
图 11 基岩面坡度为45°时孔压比分布
Figure 11. Distribution map of pore pressure ratio with the slope of the bedrock surface of 45°
图 12 基岩面坡度为60°时孔压比分布
Figure 12. Distribution map of pore pressure ratio with the slope of the bedrock surface of 60°
图 13 测点水平向反应谱（集集WGK0.2 g，5%阻尼比）
Figure 13. Response spectrum of measurement points ①~⑦ in horizontal direction（Jiji WGK0.2 g, 5% damping ratio）
图 14 测点竖向加速度反应谱（集集WGK0.2 g，5%阻尼比）
Figure 14. Measurement points ①~⑦ vertical direction acceleration response spectrum (Jiji WGK0.2 g, 5% damping ratio)
表 1 场地模型参数
Table 1. Site model parameters
位置 层数 厚度/m S波波速/（m·s−1） 网格类型 密度/（kg·m−3） 弹性模量E/MPa 静泊松比$ \mathrm{\vartheta } $ 动泊松比$ {\mathrm{\vartheta }}_{\mathrm{d}} $ 盆地内 第1层 30 205.8 CPE4 R 2 120 231 0.444 0.49 第2层 40 300 1 800 453 0.4 — 第3层 30 400 1 900 851 0.4 — 盆地外 200 600 2 000 1 920 0.334 — 
下载: 导出CSV
表 2 模型参数
Table 2. Model parameter
直径D/m 高度H/m 相对密度Dr 网格类型 单元数/个 天然孔隙比 初始有效围压/kPa 动剪应力比 0.0391 0.08 0.5 C3 D8 R 96 0.885 100 1.05
下载: 导出CSV
表 3 子程序参数
Table 3. Subprogram parameters
$ \rho $/（kg·m−3） ${V}_{{\rm{s}}}$/（m·s−1） A/B $ {\mathrm{\vartheta }}_{\mathrm{d}} $/$ \mathrm{\vartheta } $ $ {\gamma }_{0} $ a C1/C2 m/n k/b 2 120 0.08 1.02/0.36 0.49/0.444 0.000 38 0.5 0.73/0.475 0.43/0.62 0.002 5/0.5
下载: 导出CSV
表 4 仿真与试验对比
Table 4. Comparison between simulation and experiment
对比结果 试验 仿真 轴向应变峰/% 3.2 2.9 破坏时间/s 70 71.4
下载: 导出CSV
-
董菲蕃,金丹丹,陈国兴,2013. 场地条件对地表地震动特征的影响. 防灾减灾工程学报,33(6):712−718.Dong F. F., Jin D. D., Chen G. X. 2013. Effect of site condition on ground motion characteristics. Journal of Disaster Prevention and Mitigation Engineering, 33(6): 712−718.(in Chinese) 董菲蕃,陈国兴,金丹丹,2014. 泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较. 防灾减灾工程学报,34(2):154−160.Dong F. F., Chen G. X., Jin D. D. 2014. Comparison between one-dimensional equivalent linear and two-dimensional nonlinear analysis on seismic effect of Quanzhou basin. Journal of Disaster Prevention and Mitigation Engineering, 34(2): 154−160.(in Chinese) 金丹丹,陈国兴,2012. 福州盆地地震效应特征的一、二维模型对比研究. 土木工程学报,45(S1):48−53.Jin D. D., Chen G. X. 2012. Large-scale two-dimensional nonlinear FE analysis vs. one-dimensional quivalent linearization analysis on seismic effect of Fuzhou Basin. China Civil Engineering Journal, 45(S1): 48−53.(in Chinese) 金丹丹,陈国兴,董菲蕃,2014. 多地貌单元复合场地非线性地震效应特征二维分析. 岩土力学,35(6):1818−1825.Jin D. D., Chen G. X., Dong F. F. 2014. 2D analysis of nonlinear seismic effect characteristics of multi-geomorphic composite site. Rock and Soil Mechanics, 35(6): 1818−1825.(in Chinese) 李雪强,2011. 沉积盆地地震效应研究. 哈尔滨:中国地震局工程力学研究所.Li X. Q.,2011. Study on seismic effect of sedimentary basin. Harbin:Institute of Engineering Mechanics,China Earthquake Administration. (in Chinese). 刘启方,袁一凡,金星等,2006. 近断层地震动的基本特征. 地震工程与工程振动,26(1):1−10.Liu Q. F., Yuan Y. F., Jin X., et al, 2006. Basic characteristics of near-fault ground motion. Earthquake Engineering and Engineering Vibration, 26(1): 1−10.(in Chinese) 刘启方,于彦彦,章旭斌,2013. 施甸盆地三维地震动研究. 地震工程与工程振动,33(4):54−60.Liu Q. F., Yu Y. Y., Zhang X. B. 2013. Three-dimensional ground motion simulation for Shidian Basin. Journal of Earthquake Engineering and Engineering Vibration, 33(4): 54−60.(in Chinese) 刘庆忠,金丹丹,陈国兴等,2015. 河漫滩-阶地复合地形场地的非线性地震反应特征. 应用基础与工程科学学报,23(1):136−144.Liu Q. Z., Jin D. D., Chen G. X., et al, 2015. Nolinear seismic site response characteristics of composite site of alluvial flat and terraces. Journal of Basic Science and Engineering, 23(1): 136−144.(in Chinese) 刘中宪,刘明珍,韩建斌,2017. 近断层沉积盆地强地震动谱元模拟. 世界地震工程,33(4):76−86.Liu Z. X., Liu M. Z., Han J. B. 2017. Spectral-element simulation of strong ground motion in the near-fault alluvial basin. World Earthquake Engineering, 33(4): 76−86.(in Chinese) 孙锐,袁晓铭,2004. 液化土层对地表加速度反应谱的影响. 世界地震工程,20(3):33−38. doi: 10.3969/j.issn.1007-6069.2004.03.006Sun R., Yuan X. M. 2004. Effect of soil liquefaction on response spectrum of surface acceleration. World Earthquake Engineering, 20(3): 33−38.(in Chinese) doi: 10.3969/j.issn.1007-6069.2004.03.006 王海云,谢礼立,2006. 近断层强地震动的特点. 哈尔滨工业大学学报,38(12):2070−2072,2076.Wang H. Y., Xie L. L. 2006. Characteristics of near-fault strong ground motions. Journal of Harbin Institute of Technology, 38(12): 2070−2072,2076.(in Chinese) 张宇翔,袁志祥,2010. 汶川8级地震陕西灾区震害特征分析. 地震研究,33(3):329−335.Zhang Y. X., Yuan Z. X. 2010. Analysis on the characteristics of earthquake damages of the M S8.0 Wenchuan earthquake in the affected area in Shaanxi. Journal of Seismological Research, 33(3): 329−335.(in Chinese) 赵丁凤,阮滨,陈国兴等,2017. 基于Davidenkov骨架曲线模型的修正不规则加卸载准则与等效剪应变算法及其验证. 岩土工程学报,39(5):888−895.Zhao D. F., Ruan B., Chen G. X., et al, 2017. Validation of modified irregular loading-unloading rules based on Davidenkov skeleton curve and its equivalent shear strain algorithm implemented in ABAQUS. Chinese Journal of Geotechnical Engineering, 39(5): 888−895.(in Chinese) Anderson J. G., Bodin P., Brune J. N., et al, 1986. Strong ground motion from the Michoacan, Mexico, earthquake. Science, 233(4786): 1043−1049. Liu Z. X., Qiao Y. F., Cheng X. L., et al, 2023. IBEM-FEM coupling method for full process nonlinear ground motion simulation of near-fault sedimentary basins. Soil Dynamics and Earthquake Engineering, 170: 107916. doi: 10.1016/j.soildyn.2023.107916 Lu R. Q., Jiang C. S., He D. F., et al, 2024. Seismogenic fault of the 2021 MS 6.0 Luxian induced earthquake in the Sichuan Basin, China constrained by high-resolution seismic reflection and dense seismic array. Journal of Structural Geology, 179: 105050. doi: 10.1016/j.jsg.2024.105050 Martin P. P., Seed H. B. 1982. One-dimensional dynamic ground response analyses. Journal of the Geotechnical Engineering Division, 108(7): 935−952. doi: 10.1061/AJGEB6.0001316 期刊类型引用(5)
1. 张婷. 基于微信小程序的江苏地震应急联络通讯录设计与实现. 地震科学进展. 2024(02): 159-163 . 
百度学术2. 符捷，陈梅，刘霞. 基于微信小程序的医院设施设备报修系统设计. 电子设计工程. 2024(19): 58-61+66 . 百度学术3. 程树岐，王西宝，刘海刚，王帅合，赵小贺. 地震中心站观测设备实时跟踪管理系统设计与开发. 防灾减灾学报. 2023(03): 52-57 . 百度学术4. 刘彦会. 基于自适应遗传算法的微信小程序界面信息布局设计方法. 湖北科技学院学报. 2023(06): 135-140 . 百度学术5. 王渊博，李晶，马世鹏. 基于微信小程序的评价工具设计. 湘潭大学学报(自然科学版). 2022(04): 110-117 . 百度学术其他类型引用(2)
-
下载:
点击查看大图
计量
- 文章访问数: 189
- HTML全文浏览量: 72
- PDF下载量: 42
- 被引次数: 7
