Research on the Effect of Hollow Ditch in Homogeneous Soil on Vibration Insulation
-
摘要: 为分析均质地基中各参数变化对土体振动响应和空沟隔振效果的影响，通过有限元软件建立二维均质地基计算模型，分析空沟尺寸、激振频率等对土体的影响，研究了不同土体弹性模量、密度及泊松比下的振动加速度差异。结果表明，空沟深度为重要影响因素，空沟深度越深，隔振效果越好；空沟宽度对隔振效果几乎无影响；弹性模量对空沟隔振效果的影响较密度大，且弹性模量越大，隔振效果越差；密度对土体振动加速度和振幅衰减系数的影响较小；泊松比对空沟隔振效果的影响大于密度，但低于弹性模量，泊松比越大，振动加速度越小，隔振效果越差；激振频率越大，空沟隔振效果越好。Abstract: In order to analyze the specific effects of the changes of various parameters in the homogeneous foundation on the vibration isolation effect of the hollow trench, a two-dimensional homogeneous foundation calculation model was established through the finite element software. The model calculation not only analyzes the influence of the hollow trench size and the excitation frequency, but also The difference in vibration acceleration with different soil elastic modulus, density and Poisson’s ratio is discussed. The results show that the depth of the hollow trench is an important factor, and the deeper the hollow trench, the better the vibration isolation effect. The effect of the change on the vibration isolation effect is almost nothing; the effect of the elastic modulus on the vibration acceleration of the hollow trench is not strong, and it can be roughly expressed as the greater the elastic modulus, the greater the vibration acceleration. Density has little effect on vibration acceleration and amplitude attenuation coefficient of soil mass. Poisson's ratio has a greater influence on the isolation effect than the density, but it is lower than the elastic modulus. The larger the Poisson's ratio is, the smaller the vibration acceleration is, and the worse the isolation effect is. The larger the vibration frequency, the better the isolation effect.
-
引言
城市地球物理探测与传统的地球物理勘查差别很大，城市地质环境是被人类活动改造过的自然环境，实际情况会更加复杂。地下除地层、构造等天然地质体以外，还有地下管网系统、地下轨道交通、地下储存库等地下综合体。此外，城市地球物理场是动态演化的，在城市化前后甚至还会发生更加显著的变化（陈颙等，2003）。因此，城市地球物理探测具有探测深度浅、探测精度要求高、干扰因素多、场地受限等特点（李万伦等，2018）。目前国内外对目标层埋深仅有几十米（甚至十几米）的超浅层地震勘探经验不多，理论研究也不够深入。从地表至地下几十米的深度范围是各种人工地震勘探的盲区，但它对断层定位和活动性研究意义重大，近年来由于工程建设和防灾减灾的需要，已经开展了不少浅层地震勘探实践和技术研究（朱金芳等，2005）。在城市活断层探测中，用于第四纪覆盖区的活断层探测方法可分为地质学方法和地球物理方法，浅层或超浅层地震方法是地球物理方法调查城市活动断层研究中的重要手段之一（邓起东，2002；邓起东等，2003）。由于浅层地质条件以及新构造运动的特点，浅层或超浅层地震方法与常规地震勘探中采用的方法不完全相同。合理的数据采集方案是获得高分辨率地震剖面的重要基础；对浅层地震波传播速度的研究为数据处理中速度参考值的确定提供了理论依据。近年来，国内外浅层反射地震方法在城市活断层调查中的研究及其实例也相继发表，前人的研究成果反映了浅层地震方法的研究水平正在不断提高（Williams等，1995；潘纪顺等，2002；方盛明等，2002；刘保金等，2002；徐明才等，2005；段生全等，2005；杨晓平等，2007）。在利用浅层地震勘探方法确定隐伏活断层位置及其上断点埋深等工作中，一些地球物理工作者从地震波的激发、数据采集方法、室内资料处理与解释方法等方面进行了有益的探索，并取得了许多宝贵的经验，使得城市活断层浅层地震探测剖面质量和断层定位精度不断提高（何正勤等，2010；李大虎等，2010；杨歧焱等, 2015, 2016；许汉刚等，2016）。
青岛市地处胶东半岛西南部、黄海之滨，新构造运动以隆升为主，缺失第三系和第四系下部。第四纪晚期，受构造作用和海平面上升影响，接受十几米至二十余米厚以砂粘土为主的晚更新世—全新世河、海堆积物，直接覆盖于风化基岩（多为花岗岩）之上。东南部主要为花岗岩组成的低山丘陵，基岩在长期风化作用下形成了一定厚度的风化壳，并在山麓等和缓地带堆积有风化-坡积物，沿海港湾有第四纪晚期河海沉积；西北部主要是第四纪晚期的河海堆积平原或海湾，沉积有10—30m厚的第四系。因此，青岛市地震勘探属于超浅层地震勘探范畴。区域内发育若干北北东向断裂，活动性不甚强烈，地质地貌表现不很清楚，许多段落隐伏于第四系之下，断层的准确位置和活动性都不十分可靠（图 1）。与其他地区的有利探测环境相比，青岛市第四系很薄，横向物性变化大，垂向分层性差。对于不同的地震地质条件和探测目的，必须采用合理的野外采集方法和采集参数，因此，在一个新工区进行正式数据采集之前，应认真进行数据采集参数试验（何正勤等，2007）。本文通过在青岛王哥庄断裂上进行不同观测系统和采集参数的对比试验，获得了几米至二十几米深度范围内高分辨率的地震剖面，为隐伏活动断层位置及断层上断点的确定、钻孔剖面的布设和断层活动性的判定提供了可靠的科学依据。
图 1 青岛地区地质构造图Figure 1. Simplified geological map of the Qingdao areaF1：沧口断裂；F2：青岛山断裂；F3：劈石口断裂；F4：王哥庄断裂；F5：马山-王哥庄断裂1. 青岛地区地质构造概况
青岛及邻近地区在大地构造上属华北地台胶辽断隆，基底为元古代的变质岩系，经吕梁运动固结形成；古生代开始为地台发育时期，长期处于构造隆起状态，缺失古生代至三叠纪地层，致使结晶基底大面积出露。总体而言，青岛及邻近地区山势不高，以低山、丘陵和平原为主，海拔较高的山峰多为孤立分布，大致以沧口断裂为界划分为西北、东南2个不同的地貌分区。西北部地势较低，多为海湾、第四系冲洪积平原。由于新构造时期鲁东地区以隆升运动为主，青岛及邻近地区第四纪地层发育不良。许多地区基岩出露，部分地区沉积有晚更新世及全新世堆积物，局部地区发育有中更新世晚期的沉积物。市区东南部主要为花岗岩组成的低山丘陵，在山麓等和缓地带堆积有较厚的风化、坡积物，沿海港湾有第四纪晚期河海堆积；西北部是中生代地层组成的准平原和第四纪河海堆积平原（或海湾），沉积较厚的第四系，最大厚度近30米。
2. 超浅层地震勘探
2.1 青岛地区超浅层地震勘探的特点
（1）目的层比较浅，有效反射波大多在70ms以内，在反射波出现的时窗内，各种干扰波、特别是声波和面波很强，在接收反射波的后几道，出现折射波和反射波混合的情况，严重影响了有效反射波。
（2）在近震源处，震源干扰波条带内干扰严重，很难获取有效波，不同波阻抗界面形成的反射信号幅度远远小于干扰波振幅。
（3）青岛及附近地区沉积结构和地层厚度横向变化大，反射波速度和频率受此影响较大。
（4）青岛及附近地区覆盖层较薄，第四系底界面通常只有一个较强的风化界面，由于风化地层的风化程度是渐变的，风化基岩之间波阻抗差较小，第四系内部垂向分层差。基岩顶部常有数米厚的风化壳，有些基岩面甚至更浅，风化基岩与上覆地层之间的波阻抗差也不是很明显。
2.2 观测系统的设计和采集参数的选取
针对上述地震勘探特点，本文选取王哥庄断层（图 1中F4断裂）为实验对象，为了获得王哥庄断层从近地表至第四系底界面（深约30m）深度范围内的地震剖面，对同一深度范围按照探测精度由低到高的工作思路，在同一地震测线上采用不同的方法技术和探测参数，不同探测参数的地震剖面分别命名为TEST-a和TEST-b。地震数据采集使用德国DMT公司生产的SUMMIT数字地震仪，地震波激发采用了人工夯击和人工锤击，针对同一探测深度、不同精度的探测要求，采用的观测系统参数、地震波激发和数据采集参数见表 1。
表 1 观测系统参数和地震采集参数表Table 1. Survey parameters and seismic data acquisition details剖面编号 道间距/m 偏移距/m 接收道数 覆盖次数 激发震源 检波器固有频率/Hz 采样间隔/ms 记录长度/ms TEST-a 3 6 24 6 人工夯击 60 0.25 256 TEST-b 1 7 24 12 人工锤击 100 0.25 256 在不同道间距的超浅层地震记录上（图 2），图 2（a）、图 2（b）分别是TEST-a、TEST-b对应的超浅层地震共炮点记录。对比图 2（a）和图 2（b）可以看出，采用较大道间距、较小偏移距以及能量较大的震源可以获得一组很清楚的反射波组，但是近震源受干扰波影响严重，在超浅层不能获取有效反射波；采用较小道间距、较大偏移距以及能量较小震源，可以获取超浅层有效反射波，也可以减轻震源干扰波的影响，但是能量较弱，反射波信号相对变弱。
图 2 不同道间距的超浅层地震共炮点记录Figure 2. The common-shot gather of ultra-shallow seismic records in different trace intervals采用较小道间距和较弱震源，可以获取超浅层有效反射波，且垂向分层较好。由于采用小排列接收，道间距越小反射波能量衰减也越慢，因此受震源干扰波影响也变的越小，但是面波比较发育。由此可见，较小道间距和较小偏移距有利于保护反射波，延缓反射波的衰减，但是也带来了面波发育的问题；在保证探测深度要求的前提下，较弱的震源有利于避开震源干扰波的影响。
2.3 超浅层地震数据处理中的方法
采用Promax地震反射处理系统处理地震数据。在数据处理过程中，由于获得的原始记录的反射波近震源干扰严重，远震源反射波能量较强、干扰能量较弱等特点，因此选取记录信噪比高的远震反射波，可最大限度地提高地震资料的分辨率。其数据处理流程和方法主要包括：道编辑、静校正、真振幅恢复、地表一致性振幅处理、地表一致性反褶积、二维随机噪音衰减、f-k滤波、正常时差校正（NMO）、共中心点（CMP）叠加、叠后频率带通滤波和叠后剖面去噪等。
本文利用不同的滤波和处理方法来提高地震数据信噪比，图 3为同一记录数据在不同情况下进行f-k滤波获得的对比分析结果和频谱分析结果。其中，图 3（a）为原始数据和相对应的谱分析结果，可以看出在超浅层地震勘探中，地震反射波受浅层折射、面波和声波的影响，反射波信号相对较弱（图中红色圆圈），信噪比不高，通过频谱分析则显示地震数据主频有2个频段：50Hz和100Hz；图 3（b）为对原始数据应用f-k滤波后的结果及其频谱分析，由图可见应用f-k滤波后，数据中面波和声波得到有效减弱，反射波信号得到加强，信噪比也得到提高，频谱分析显示数据主频为100Hz，同时结果还显示面波频率范围与有效地震反射波接近，滤波困难，数据信噪比低；图 3（c）通过切除数据前几道信号，避免了近道面波对反射波的干扰，再进行f-k滤波，反射波信号明显增强，信噪比也得到进一步提高。对比结果表明，在浅层或超浅层地震勘探中，在取得信噪比较高的情况下，再对地震数据进行f-k滤波可以有效去除干扰波，提高地震反射波的信噪比。
图 3 原始数据及不同信噪比下的f-k滤波对比及其频谱分析图（一）Figure 3. Analysis between f-k filtering under different signal-to-noise ratio data and corresponding spectrum2.4 超浅层地震剖面特征
图 4、图 5分别是TEST-a、TEST-b测线经过数据处理后的地震反射时间及地质解释剖面（图中方向左西右东）。在图 4地震反射时间剖面中，丰富的地层界面反射波组出现在剖面100ms以上的时间范围内，存在1组起伏变化的反射波组，由2—3个强相位组成；在图 5地震反射时间剖面中，丰富的地层界面反射波组出现在剖面60ms以上的时间范围内，存在2组起伏变化的反射波组，分别由2—3个强相位组成。根据剖面上反射波同相轴的横向变化特征和上下反射波组关系，在该测线剖面上解释了1个断点，在图 4和图 5中标识为F。断层F位于TEST-a测线CDP423#和TEST-b测线CDP87#处，反射波同相轴出现了错动，断层为倾向西、高倾角逆断层，根据测区地质资料，判断该断层为王哥庄断裂。
图 4 TEST-a地震反射时间及地质解释剖面Figure 4. Stack time section of seismic reflection and the interpretation in TEST-a
图 5 TEST-b地震反射时间及地质解释剖面Figure 5. Stack time section of seismic reflection and the interpretation in TEST-b对比TEST-a、TEST-b测线的地震反射波叠加剖面发现，TEST-b测线在超浅层部分，特别是30ms以上部分有清晰的反射同相轴，说明小道距在超浅层地震勘探中，能揭示更多的断层结构细节。
3. 钻孔地质剖面及综合解
在王哥庄北TEST-a测线CDP423#和TEST-b测线CDP87#处的王哥庄断裂断点附近，布设了1组深19.1—23.5m、相邻孔距5—10m的钻孔（ZK1—ZK7），发现在ZK2与ZK6之间基岩顶面和第四系底部被断层错断4m左右，西升东降（图 6）。由图 6可见，除底部层位外，第四系各个层位均连续分布，未受断层构造变动，根据前期第四纪地层划分结果（包括测年结果），这套地层应为上更新统。因此，该处王哥庄断裂在晚更新世早期仍在活动。
图 6 TEST-a和TEST-b测线王哥庄断裂钻孔联合地质剖面Figure 6. The drilling geological section of WGZ fracture combined with the TEST-a and TEST-b survey lines4. 结论与讨论
（1）本文从城市活断层探测的角度出发，利用地质学、地球物理学等学科中的相关理论和方法，探讨超浅层地震在复杂地质构造背景下获得有效探测结果的前提条件。综合分析青岛地区地震勘探资料发现：①要依据扩展排列试验并参考观测系统参数选择原则，选好最小和最大偏移距，在超浅层探测中，近震源道记录中反射波与干扰波很难分辨，需要一定的偏移距才能将其分辨清楚；②最大偏移距不能过大，由于目标层埋深很浅，若排列长度过大，记录的多为折射波，所以需要采用短排列、小道间距进行数据采集。
（2）在对超浅层地震数据采集技术、资料处理与解释方法进行系统研究的基础上，以可清晰分辨出断层位置的典型剖面为例，证明在选择合适的观测系统和数据处理方法基础上，超浅层地震反射波法可以获取深度仅有十几米的地层反射信号，且大部分反射剖面都可较清楚地揭示出超浅部断层位置和断层特征。青岛王哥庄断裂地震勘探研究结果显示，超浅层地震反射方法在城市活断层探测中，特别是沿海覆盖层较浅的区域具有良好的应用前景。
（3）在对青岛王哥庄断裂超浅层地震反射时间剖面开展研究和解释的基础上，结合典型断点的钻孔联合剖面探测资料，研究了王哥庄断裂的准确位置、几何结构和活动性。研究表明，王哥庄断裂晚更新世以来可能也是活动的。
-
图 6 不同弹性模量下积分振幅衰减系数
Figure 6. Integrating amplitude attenuation coefficient under different elastic modulus
图 8 不同密度下积分振幅衰减系数
Figure 8. Integrating amplitude attenuation coefficient under different densities
图 10 不同泊松比下积分振幅衰减系数
Figure 10. Integral amplitude attenuation coefficient under different Poisson ratios
图 11 不同激振频率下积分振幅衰减系数
Figure 11. Integral amplitude attenuation coefficient under different frequencies
表 1 弹性模量计算工况
Table 1. Working condition of elastic modulus calculation
项目 工况1 工况2 工况3 弹性模量/MPa 60 141 25 压力波速/（m·s−1） 257.2 394.8 166.3 剪切波速/（m·s−1） 117.20 182.00 75.52 
下载: 导出CSV
表 2 土体密度计算工况
Table 2. Working conditions of soil density calculation
项目 工况1 工况2 工况3 密度/（kg·m−3） 1 600 1 800 2 000 压力波速/（m·s−1） 257.2 242.8 230.4 剪切波速/（m·s−1） 117.2 110.3 104.7
下载: 导出CSV
表 3 土体泊松比计算工况
Table 3. Calculation condition table of Poisson's ratio of soil
项目 工况1 工况2 工况3 土体泊松比 0.37 0.10 0.20 压力波速/（m·s−1） 257.2 195.8 204.2 剪切波速/（m·s−1） 117.2 130.6 125.0
下载: 导出CSV
表 4 激振频率计算工况
Table 4. Excitation frequency calculation condition table
项目 工况1 工况2 工况3 激振频率/Hz 40 70 100 压力波速/（m·s−1） 257.2 257.2 257.2 剪切波速/（m·s−1） 117.2 117.2 117.2
下载: 导出CSV
-
巴振宁, 王靖雅, 梁建文, 2016. 层状地基中隔振沟对移动列车荷载隔振研究-2.5维IBEM方法. 振动工程学报, 29（5）: 860—873 doi: 10.16385/j.cnki.issn.1004-4523.2016.05.013Ba Z. N. , Wang J. Y. , Liang J. W. , 2016. Reduction of train-induced vibrations by using a trench in a layered foundation. Journal of Vibration Engineering, 29(5): 860—873. (in Chinese) doi: 10.16385/j.cnki.issn.1004-4523.2016.05.013 陈昆, 贾霄, 刘彬等, 2014. 高速铁路空沟隔振措施隔振效果的有限元分析. 地震工程学报, 36（3）: 575—579 doi: 10.3969/j.issn.1000-0844.2014.03.0575Chen K. , Jia X. , Liu B. , et al. , 2014. Finite element analysis of reducing the high-speed railway vibration effect on the environment using open trenches. China Earthquake Engineering Journal, 36(3): 575—579. (in Chinese) doi: 10.3969/j.issn.1000-0844.2014.03.0575 陈曦仪, 2003. 含水膜包对振波阻隔效应之试验研究. 台湾: 国立成功大学. 邓亚虹, 夏唐代, 陈敬虞, 2006. 车辆动荷载作用下隔振沟响应增强区数值分析. 岩土工程学报, 28（12）: 2121—2127 doi: 10.3321/j.issn:1000-4548.2006.12.014Deng Y. H. , Xia T. D. , Chen J. Y. , 2006. Numerical analysis of response magnified area of vibration isolation trenches subjected to vehicle dynamic loads. Chinese Journal of Geotechnical Engineering, 28(12): 2121—2127. (in Chinese) doi: 10.3321/j.issn:1000-4548.2006.12.014 丁智, 张霄, 吴敏慧等, 2019. 双线地铁运营隔振沟屏障性能研究. 地震工程学报, 41（1）: 9—15Ding Z. , Zhang X. , Wu M. H. , et al. , 2019. Vibration isolation effect of isolation trenches during double-line subway operation. China Earthquake Engineering Journal, 41(1): 9—15. (in Chinese) 高广运, 张婧, 郭星宇, 2007. 土体参数对扩底墩竖向承载性状影响的数值分析. 岩土工程技术, 21（5）: 230—233, 238 doi: 10.3969/j.issn.1007-2993.2007.05.004Gao G. Y. , Zhang J. , Guo X. Y. , 2007. Numerical analysis of vertical bearing behavior of belled pier under different soil parameters. Geotechnical Engineering Technique, 21(5): 230—233, 238. (in Chinese) doi: 10.3969/j.issn.1007-2993.2007.05.004 高广运, 李宁, 何俊锋等, 2011. 列车移动荷载作用下饱和地基的地面振动特性分析. 振动与冲击, 30（6）: 86—92 doi: 10.3969/j.issn.1000-3835.2011.06.019Gao G. Y. , Li N. , He J. F. , et al. , 2011. Analysis of ground vibration generated by train moving loads on saturated soil. Journal of Vibration and Shock, 30(6): 86—92. (in Chinese) doi: 10.3969/j.issn.1000-3835.2011.06.019 关歆莹, 刘超, 2011. 地下铁道的振动及其控制措施的研究. 震灾防御技术, 6（1）: 77—84 doi: 10.3969/j.issn.1673-5722.2011.01.008Guan X. Y. , Liu C. , 2011. The study on the vibration of metro and its control measures. Technology for Earthquake Disaster Prevention, 6(1): 77—84. (in Chinese) doi: 10.3969/j.issn.1673-5722.2011.01.008 李建平, 张晓磊, 冯世进, 2020. 水位变化对多空沟屏障远场隔振效果的影响研究. 岩土力学, 41（9）: 3131—3138, 3147Li J. P. , Zhang X. L. , Feng S. J. , 2020. Research on the effect of the moving water table on passive isolation effectiveness using a multiple open trench barrier. Rock and Soil Mechanics, 41(9): 3131—3138, 3147. (in Chinese) 娄建武, 龙源, 周翔等, 2004. 爆炸波传播的沟槽效应及分析. 矿冶工程, 24（1）: 11—15 doi: 10.3969/j.issn.0253-6099.2004.01.004Lou J. W. , Long Y. , Zhou X. , et al. , 2004. Damping ditch effect in propagation of blasting seismic wave and its analysis. Mining and Metallurgical Engineering, 24(1): 11—15. (in Chinese) doi: 10.3969/j.issn.0253-6099.2004.01.004 时刚, 高广运, 2011. 饱和地基中二维空沟远场被动隔振研究. 振动与冲击, 30（9）: 157—162 doi: 10.3969/j.issn.1000-3835.2011.09.034Shi G. , Gao G. Y. , 2011. Two-dimensional analysis of open trench used as passive barriers in saturated soil. Journal of Vibration and Shock, 30(9): 157—162. (in Chinese) doi: 10.3969/j.issn.1000-3835.2011.09.034 王亚秀, 1989. 隔振沟对加速度衰减的影响. 内陆地震, 3（1）: 56-59. doi: 10.16256/j.issn.1001-8956.1989.01.009 肖世伟, 雷长顺, 郭超等, 2011. 高速铁路空沟隔振数值分析. 铁道工程学报, 28（8）: 23—29 doi: 10.3969/j.issn.1006-2106.2011.08.005Xiao S. W. , Lei C. S. , Guo C. , et al. , 2011. Numerical analysis of vibration reduction of open trench for high-speed railway. Journal of Railway Engineering Society, 28(8): 23—29. (in Chinese) doi: 10.3969/j.issn.1006-2106.2011.08.005 杨贝贝, 李伟华, 赵成刚等, 2017. 基于Comsol的饱和多孔介质动力方程的数值模拟及应用. 地震工程学报, 39（2）: 321—328Yang B. B. , Li W. H. , Zhao C. G. , et al. , 2017. Numerical simulation and application of dynamic equations of saturated porous media based on Comsol. China Earthquake Engineering Journal, 39(2): 321—328. (in Chinese) 姚锦宝, 夏禾, 胡敬梁, 2018. 空沟对列车运行产生的环境振动隔振效果研究. 中国铁道科学, 39（2）: 44—51Yao J. B. , Xia H. , Hu J. L. , 2018. Study on vibration isolation effect of open trench on environmental vibration induced by train operation. China Railway Science, 39(2): 44—51. (in Chinese) 叶茂, 曹保兴, 郑志华等, 2015. 列车动荷载下某古遗址隔振沟的减振效果研究. 应用力学学报, 32（4）: 682—688Ye M. , Cao B. X. , Zheng Z. H. , et al. , 2015. Effect of vibration isolation trench on an ancient ruin subjected to train dynamic loads. Chinese Journal of Applied Mechanics, 32(4): 682—688. (in Chinese) 郁雯, 刘杰, 刘航等, 2021. 连续型隔振屏障在高速铁路环境振动问题治理中的应用研究. 震灾防御技术, 16（1）: 157—164 doi: 10.11899/zzfy20210116Yu W. , Liu J. , Liu H. , et al. , 2021. Based on the infinite element boundary analysis, the study of the effect of continuous barrier on the environmental vibration of high-speed railway. Technology for Earthquake Disaster Prevention, 16(1): 157—164. (in Chinese) doi: 10.11899/zzfy20210116 周凤玺, 马强, 周志雄, 2020. 二维地基中空沟−波阻板联合隔振屏障分析. 岩土力学, 41（12）: 4087—4092, 4115Zhou F. X. , Ma Q. , Zhou Z. X. , 2020.2 D analysis of vibration-isolation efficiency of an open trench-wave impedence block barrier. Rock and Soil Mechanics, 41(12): 4087—4092, 4115. (in Chinese) -
下载:
点击查看大图
计量
- 文章访问数: 190
- HTML全文浏览量: 22
- PDF下载量: 35
- 被引次数: 0
