THE CALCULATION METHOD FOR OSMOTIC SUCTION OF SALINE SOLUTION
-
摘要: 膨润土垫层在高放废料处置库环境中与围岩接触，力学性能会受到围岩裂隙中所含盐溶液的影响。盐溶液的渗透吸力会在膨润土上产生类似于竖向荷载作用的附加应力，量化盐溶液对膨润土力学性能的影响对评估地下处置库的安全性具有重要意义。溶液渗透吸力系数作为计算渗透吸力的关键，目前需要通过较为复杂的实验测得，对工程实际应用造成了阻碍。通过引入Debye-Hückel公式，提出含单价离子电解质、2-2型电解质及混合电解质溶液的渗透吸力系数及渗透吸力的计算方法。基于Debye-Hückel公式，分析溶剂种类和温度对渗透吸力系数的影响，结果表明：溶剂极性越大，渗透吸力系数越大；温度越高，渗透吸力系数越小。
-
关键词:
- 膨润土 /
- 盐溶液 /
- 渗透吸力 /
- Debye-Hückel公式
Abstract: When bentonite cushion contacts with surrounding rocks in a disposal repository of highly radioactive waste, its mechanical properties will be affected by the saline solution contained in surrounding rock cracks, and the osmotic suction of salt solution will add additional stress similar to vertical load on bentonite. Therefore, quantifying the effect of saline solution on the mechanical properties of bentonite is of great significance for evaluating the safety of underground repository. Osmotic coefficients as the key to calculate the osmotic suction are usually obtained by complex tests, which obstructs the practical engineering applications. The Debye-Hückel formula are introduced to calculate the osmotic coefficients and osmotic suction for monovalent ion electrolyte, 2-2 type electrolyte and mixed electrolyte solutions. The impact of various factors on solution suction coefficients, such as solvent type, solute type, solute concentration and temperature are analyzed. Then the calculation method of osmotic suctions of different solutions are obtained. Based on the Debye-Hückel formula, the influence of solvent type and temperature on the osmotic suction coefficient is analyzed. The results show that the larger the solvent polarity, the larger the osmotic suction coefficient; the higher the temperature, the smaller the osmotic suction coefficient.-
Key words:
- bentonite /
- saline solution /
- osmotic suction /
- Debye-Hücke formula
-
东濮凹陷位于渤海湾盆地西南临清坳陷的南部, 凹陷内发育规模和活动周期不尽相同的断裂, 这些断裂的长期继承性发育不仅控制了东濮凹陷的形成和发育, 也控制了盆地内油气分布与富集[1~7]。濮卫环洼带隶属于东濮凹陷中央隆起带, 为北部卫东-文东断裂系和濮城-陈营断裂系相向下落所形成的“断洼型”负向构造单元, 地质结构复杂, 储层空间展布变化大。目前, 构造高部位的勘探程度较高, 而洼陷带勘探程度还较低。因此, 深入研究该区断裂特征及其对油气运聚的作用, 对滚动勘探开发和隐蔽油气藏勘探具有重要意义。渤海湾及其周缘乃至华北地区油气成藏条件与油气运聚过程研究日亦深入, 已取得许多新成果[8~19]可以借鉴。
1. 断裂发育特征
1.1 卫东-文东断裂系
卫东-文东断裂系由复杂的卫东断裂破碎带和文东断层、文51断层等组成(图 1、表 1)。
图 1 濮卫环洼带主要断裂分布图Figure 1. Distribution chart of main faults in the Pu-Wei ring depression表 1 濮卫环洼带主要断裂要素表Table 1. Main fault elements in Pu-Wei ring depression
(1) 卫东断层
由3条呈雁行式排列的断层组成, 是一组东倾的滑脱性大断层, 北起古云集地区, 南至卫70井附近消失。断层在沙二下亚段沉积晚期产生, 至东营组沉积末期停止活动。断层走向NNE, 倾向SEE, 倾角15~50°, 落差为50~1130m。剖面上, 在沙三段盐层顺层滑动, 并在盐层内消失, 向下不切割基底。卫东断层下降盘伴生一系列反向“Y”字型次级断层, 从深层到浅层多期发育, 明显控制沙二段、沙一段及东营组等地层的沉积。
(2) 文东断层
文东断层是环洼带南部的西界断层, 向南断层落差逐渐增大。在研究区内, 断层走向NE、倾向SE, 延伸长度5km, 剖面上表现为上缓下陡, 落差为30~250m, 沙三下亚段沉积末期开始活动, 沙三中沉积期强烈活动, 至东营组沉积末期停止活动。
(3) 文51断层
属于卫东断层的伴生断层, 与卫东断层斜交, 呈反“ Y”形, 断层走向NEE, 倾向NWW, 倾角10~40°, 落差50~260m, 南大北小, 延伸长度约2.5km, 于沙二上至沙一期活动, 沙二上沉积期对沉积的控制作用明显。
1.2 濮城-陈营断裂系
濮城-陈营断裂系展布于濮城和陈营地区, 发育濮城南断层、濮67断层、濮31断层、濮138断层和陈营断层(见图 1、表 1)。
(1) 濮城南断层
断层走向NE, 倾向NW, 倾角25~60°, 落差50~560m, 延伸长度大于8.5km。断层活动时间为沙三中亚段早期至东营组沉积期。
(2) 濮67断层
濮67断层位于濮城-陈营断裂系的中段, 与濮城南断层和陈营断层首尾相接, 走向NNE, 倾向NWW, 延伸长度约10km。剖面上具有上陡下缓的特点, 倾角20~60°, 断层落差为50~360m, 北大南小。产生于沙三中亚段沉积期, 结束于沙二上亚段沉积末期。
(3) 濮31断层
濮31断层断面下陡上缓, 为濮67断层的伴生断层, 从沙三中后期一直活动至东营组沉积期。断层走向NNE、倾向NWW, 倾角40~60°, 落差50~200m, 延伸长度约10km。
(4) 濮138断层
位于濮31断层的西侧, 断层走向NNE向, 倾向NWW, 倾角40~60°, 落差为50~350m, 延伸长度约8km。剖面上, 向深部交接到濮67断层上。断层活动时间为沙三中-上亚段沉积期。
(5) 陈营断层
陈营断层是一条切割基底的区域性控制沉积断层, SN走向、倾向正西, 断层落差为50 ~1000m, 南端小、北端大, 延伸长度约10km。断裂活动开始于沙三下亚段沉积早期, 结束于沙二上亚段沉积晚期, 沙三下-中期剧烈活动。陈营断层与濮67断层呈雁行式排列, 向北斜交于古云集的云3断层。
2. 断层活动性分析
目前主要采用断层生长指数、断层落差、断层活动速率等参数来定量表示断层的活动性[20~22]。濮卫环洼带主要断层生长指数和断层落差统计表明, 同一条断层的不同部位、甚至同一部位的不同时期, 其活动性都存在着较大的差异。
卫东断层是卫东-文东断裂系的主体, 沙二上亚段至东营组沉积期间, 持续较强的断裂活动, 是一条长期活动的控洼断层。卫东三号断层活动范围最广泛, 活动强度大, 断层最大落差高达1130m。东营组沉积期, 卫东断层的活动性最强, 断层中段生长指数高达3.5, 在601测线上断层落差达518m。沿走向分段性明显, 断层中段活动强度最大, 向两端活动性减弱(图 2a、b)。
图 2 濮卫环带主要断层生长指数与落差直方图Figure 2. The growth index and paleo-fall bar chart of the main faults in the Pu-Wei ring depression濮城南断层位于濮城-陈营断裂系南端, 活动时间较长, 但整体活动性较弱, 最大落差为560m。沙三中-上亚段沉积期是该断裂主要活动期, 活动范围最广泛, 最大生长指数为1.8, 901测线上的断层落差达460m。沙二下开始活动强度逐渐减弱, 活动范围也缩小至中南部, 至东营期仅在1001测线附近活动, 落差为60m (图 2c、d)。
濮67断层活动时间相对较短, 沙三中-上亚段沉积时期断层活动最为强烈, 沙二下次之, 沙二上活动强度最小。沙三中-上亚段沉积时期, 在701~801测线附近活动性最强, 断层生长指数近1.8, 最大落差达360m, 为濮卫环洼带早期的控洼断裂之一(图 2e、f)。
濮31断层活动时间长, 不同时期、不同部位的断层活动强度变化较大。从断层生长指数来看, 沙三中-上沉积期断层活动性最强, 沙二下开始断层活动相对减弱, 沙一和东营期仅在741测线附近活动, 断层活动微弱。(图 2g)。
濮138断层活动时间较短, 沙三中亚段沉积期断层活动强度弱, 至沙三上亚段沉积期断层活动较为强烈, 且向西南方向活动性增强(图 2h)。
3. 断裂与油气运聚的关系
3.1 断层与储集体的配置关系
断层和连通砂体是濮卫环洼带最主要的运移通道, 二者共同构成了断层型、断-砂组合型油气输导体系, 它们在平面分布以及成藏模式上均具有较大的差异。
(1) 断层输导体系
该输导体系类型分布较为局限, 仅分布于文明寨的主体高部位、濮城断阶带和卫城卫东断层附近, 其中以文明寨主体部位最具代表性。文明寨主体构造位置较高, 自身生烃能力差, 区内断裂发育, 由于卫东断层沟通了烃源岩和储集层, 油气沿断层面向上垂向运移, 在适当部位聚集成藏, 油气的垂向运移距离远远大于侧向运移距离。
据地层水矿化度对比分析, 濮卫环洼带沙二段地层水具有异常高矿化度, 其值与下部沙三段地层水矿化度较接近。平面上, 异常高矿化度的CaCl2型地层水主要沿濮31、濮城南、卫东断层分布(图 3), 表明这些断层是流体运移的主要通道, 它们沟通了深部的地层, 使高矿化度地层水沿断裂带向上运移。
图 3 濮卫环洼带沙二段异常高矿化度地层水平面分布图Figure 3. Distribution of exceedingly high-mineralization water in Sha-2 Member in the Pu-Wei ring depression然而, 断层输导体系若输送大量的油气, 往往需要断层沟通有效烃源岩, 且在烃源岩层系内存在大量的砂体, 这些砂体可以汇聚生烃层系的油气, 然后借助于断层的沟通作用向上部垂向运移, 运移至地层倾角不大的地层进入断层面两侧砂体中形成油气藏。
(2) 断-砂复合输导体系
断层与砂体相互配置构成了濮卫环洼带运载油气的主要输导体系, 分布范围较广, 在卫城和濮城的主体部位最具代表性。目前已经发现的濮城油田、濮城南地区的文213含油气区等, 均具有此类运聚特征。
濮卫环洼带砂岩厚度较大, 断-砂匹配关系良好, 构成了良好的输导体系。该区域的多套砂岩储层呈层状分布于洼陷内及斜坡带, 纵向上每套输导层之上均有泥岩和盐岩盖层封隔, 以侧向输导为主(图 4)。此时, 砂体主要起运移通道的作用, 而断层主要起阻挡油气, 并调节油气运移层系和方向的作用。因此, 随着油气运移层系的不断变浅, 运移动力具有明显的衰减, 在断层活动停止、封堵条件较好的地区聚集成藏, 形成主断层附近的构造-地层油气藏和远离主断层的边部岩性油气藏。
3.2 断层活动时间与油气运聚时间的匹配关系
断层活动时间与大规模油气运聚时间的匹配关系, 决定了断层在油气运移成藏中的贡献, 同时也决定了油气藏的分布特征[23]。濮卫环洼带烃源岩生排烃史、油气藏饱和压力及流体包裹体分析表明, 东营组沉积期沙三段烃源岩已经进入生烃门限, 是油气生成和成藏高峰期。东营组地层的剥蚀, 使烃源岩热演化生烃受到抑制甚至停止。新近纪以来研究区再次下沉接受沉积, 至新近纪晚期发生明显的二次生烃作用, 但此时构造活动弱, 生排烃量相对较小①。故该地区成藏关键时期是在东营组沉积时期[24]。
① 高平.濮卫洼陷岩性油气藏成藏规律与精细勘探[R].河南濮阳:中原油田分公司勘探开发科学研究院, 2008.
濮卫环洼带主要断层活动时间与油气运移期对比分析表明(表 2), 卫东断层活动时间长, 油气大规模运移时期强烈活动, 与成藏时间匹配好, 是油气从深层向浅层运移的重要通道。文东、濮31及濮城南等断层活动强度微弱, 活动范围局限, 但与油气大规模运移时间相匹配, 可充当油气运移的有效通道。文51、濮67、濮138和陈营等断层活动时间较早, 东营期已经停止活动或基本停止活动, 主要对油气的聚集成藏起圈闭遮挡的作用。
表 2 濮卫环洼带主要断层活动时间与油气运移期对比表Table 2. Comparison of growth index to main fault and migration phase of the hydrocarbon
此外, 断层本身活动的不均衡性也使油气运聚成藏情况复杂化。同一条断层(或断层的某一段)在不同时期或不同部位, 由于活动强度的不一致, 对油气成藏所起的作用也不尽相同, 某一时期作为油气运移的通道, 输导油气, 另一时期则成为阻止油气运移的屏障, 使油气在断层附近聚集成藏。
4. 结论
(1) 濮卫环洼带主要发育卫东-文东断裂系和濮城-陈营断裂系, 均为较大规模的二级断裂。卫东-文东断裂系由复杂的卫东断裂破碎带和文东断层、文51断层等组成, 濮城-陈营断裂系自东而西发育陈营断层、濮67断层、濮城南断层、濮31断层和濮138断层。
(2) 卫东、文东、濮城南及濮31等断层活动时间较长, 活动强度不一致, 与成藏时间匹配好, 是油气从深层向浅层运移的重要通道, 文东、濮31及濮城南等断层可充当油气运移的有效通道。
(3) 断层和连通砂体构成了断层型、断-砂组合型油气输导体系。断层型输导体系沟通了烃源岩和储集层, 增强了油气垂向运移能力, 断-砂组合型油气输导体系以侧向输导为主, 断层主要起阻挡油气, 并调节油气运移层系和方向的作用。这些结论对本区油气田滚动勘探开发有重要意义。
-
图 1 钠盐在水溶液中渗透吸力系数与溶质浓度关系
Figure 1. Relationship between osmotic suction coefficient and solute concentration for sodium salt in agueous solution
图 2 NaCl和CaCl2的渗透吸力与溶质浓度的关系
Figure 2. Relationship between osmotic suction and solute concentration of NaCl and CaCl2
图 3 三种硫酸盐水溶液中的渗透吸力系数与溶质浓度的关系
Figure 3. Relationship between osmotic suction coefficient and solute concentration for three kinds of sulfate solutions
图 4 NaCl-KCl混合溶液渗透吸力系数与浓度的关系
Figure 4. Relationship between osmotic suction coefficient and concentration for NaCl-KCl mixed solutions
图 5 渗透吸力系数随温度变化关系
Figure 5. Relationship between osmotic suction coefficient and temperature
表 1 不同种类含单价离子电解质的渗透吸力系数
Table 1. Osmotic suction coefficients of different monovalent ionic electrolytes
c(mol/kg) 1-1型 2-1型 3-1型 1-2型 1-3型 1-4型 1-5型 KCl Ca(NO3)2 CaCl2 AlCl3 Na2SO4 K3PO4 Na4P2O7 Na5P3O4 0.01 0.967 0.904 0.907 0.842 0.898 0.819 0.739 0.667 0.05 0.940 0.850 0.864 0.810 0.828 0.738 0.644 0.576 0.1 0.926 0.831 0.856 0.819 0.793 0.708 0.589 0.510 0.2 0.913 0.820 0.862 0.850 0.754 0.680 0.520 0.432 0.3 0.906 0.818 0.876 0.891 0.728 0.667 0.489 0.397 0.4 0.902 0.819 0.895 0.942 0.707 0.662 0.483 0.376 0.5 0.900 0.823 0.916 1.002 0.691 0.661 0.496 - 0.7 0.898 0.836 0.964 1.143 0.665 0.666 - - 1 0.898 0.861 1.048 1.390 0.641 - - - 1.6 0.905 0.917 1.243 1.953 0.622 - - - 2 0.913 0.951 1.386 - 0.624 - - - 2.5 0.924 0.904 1.572 - 0.638 - - - 3 0.937 0.850 - - 0.662 - - - 4 0.965 0.831 - - 0.734 - - - 4.8 0.989 0.820 - - - - - - 
下载: 导出CSV
表 2 不同溶质的渗透吸力系数
Table 2. Osmotic suction coefficients of different solutes
c(mol/kg) MgSO4 ZnSO4 CuSO4 CaSO4 0.001 0.894 0.895 0.888 0.910 0.005 0.791 0.793 0.778 0.792 0.01 0.741 0.741 0.725 - 0.05 0.635 0.623 0.611 - 0.1 0.595 0.573 0.560 - 0.5 0.525 0.482 0.467 - 1 0.527 0.477 0.463 - 1.4 0.562 0.510 0.488 2 0.660 0.606 - - 3 0.913 0.860 - - 3.5 - 1.022 -
下载: 导出CSV
表 3 不同溶质的Debye-Hückel系数
Table 3. Debye-Hückel coefficients of different solvents
溶剂 εr ρsl(kg/m3) A 水 78.4 1000 0.393 甲醇 32.63 786.6 1.297 乙醇 24.35 785 2.009 丙酮 20.7 784.35 2.562 乙酸乙酯 6.02 894.5 17.447 苯 2.28 873.6 73.974
下载: 导出CSV
表 4 水在不同温度中的D-H系数
Table 4. Debye-Hückel coefficients of water under different temperatures
温度/K 相对介电常数εr 参数A 273 87.74 0.3781 278 85.76 0.3808 283 83.83 0.3836 288 81.95 0.3866 293 80.1 0.3899 294 79.73 0.3906 296 79.02 0.3918 298 78.3 0.3933 300 77.6 0.3946 303 76.55 0.3968 308 74.83 0.4006 313 73.15 0.4046 318 71.15 0.4119 323 69.91 0.4131
下载: 导出CSV
-
[1] 王驹, 陈伟明, 苏锐, 等.高放废物地质处置及其若干关键科学问题[J].岩石力学与工程学报, 2006, 4(25):801~812. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb200604015WANG Ju, CHEN Weiming, SU Rui, et al. Geological disposal of high-level radioactive waste and its key scientific issues[J]. Chinese Journal of Rock Mechancis and Engineering, 2006, 4(25):801~812. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb200604015 [2] 孙德安, 张龙.盐溶液饱和高庙子膨润土膨胀特性及预测[J].岩土力学, 2013, 34(10):2790~2795. http://d.old.wanfangdata.com.cn/Periodical/ytlx201310007SUN Dean, ZHANG Long. Swelling characteristics of Gaomiaozi bentonite saturated by salt solution and their prediction[J]. Rock and Soil Mechanics, 2013, 34(10):2790~2795. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/ytlx201310007 [3] Sun Z M, Yu J, Zheng S L, et al. Effect of salt in aqueous solution on the swelling and water-retention capacity of bentonite[J]. Advanced Materials Research, 2011, 194~196:2039~2045. doi: 10.4028/www.scientific.net/AMR.194-196 [4] Alawaji H A. Swell and compressibility characteristics of sand-bentonite mixtures inundated with liquids[J]. Applied Clay Science, 1999, 15(3/4):411~430. http://cn.bing.com/academic/profile?id=11109b69d94d9eb124098a127c33ed9b&encoded=0&v=paper_preview&mkt=zh-cn [5] Zhang L, Sun D A, Jia D. Shear strength of GMZ07 bentonite and its mixture with sand saturated with saline solution[J]. Applied Clay Science, 2016, 132~133:24~32. doi: 10.1016/j.clay.2016.08.004 [6] Tabiatnejad B, Siddiqua S, Siemens G. Impact of pore fluid salinity on the mechanical behavior of unsaturated bentonite-sand mixture[J]. Environmental Earth Sciences, 2016, 75(22):1434. doi: 10.1007/s12665-016-6246-5 [7] Maio D C. Shear strength of clays and clayey soils: the influence of pore fluid composition[A]. Loret B, Huyghe J M. Chemo-Mechanical Couplings in Porous Media Geomechanics and Biomechanics[M]. Vienna: Springer, 2004, 45~55. [8] Rao S M, Thyagaraj T. Swell-compression behaviour of compacted clays under chemical gradients[J]. Canadian Geotechnical Journal, 2007, 44(5):520~532. doi: 10.1139/t07-002 [9] Barbour S L, Fredlund D G. Mechanisms of osmotic flow and volume change in clay soils[J]. Canadian Geotechnial Journal, 1989, 26(4):551~562. doi: 10.1139/t89-068 [10] Xu Y F, Xiang G S, Jiang H, et al. Role of osmotic suction in volume change of clays in salt solution[J]. Applied Clay Science, 2014, 101:354~361. doi: 10.1016/j.clay.2014.09.006 [11] Miller D J, Nelson J D. Osmotic suction in unsaturated soil mechanics[A]. Proceedings of the Fourth International Conference on Unsaturated Soils[C]. Arizona:American Society of Civil Engineers, 2006, 1382~1393. [12] Apelblat A, Dov M, Wisniak J, et al. The vapour pressure of water over saturated aqueous solutions of malic, tartaric, and citric acids, at temperatures from 288 K to 323 K[J]. The Journal of Chemical Thermodynamics, 1995, 27(1):35~41. doi: 10.1006/jcht.1995.0004 [13] Wijmans J G, Baker R W. The solution-diffusion model:a review[J]. Journal of Membrane Science, 1995, 107(1~2):1~21. doi: 10.1016/0376-7388(95)00102-I [14] Fernández D P, Goodwin A R H, Lemmon E W, et al. A formulation for the static permittivity of water and steam at temperatures from 238 K to 873 K at pressures up to 1200 MPa, including derivatives and Debye-Hückel coefficients[J]. Journal of Physical and Chemical Reference Data, 1997, 26(4):1125~1166. doi: 10.1063/1.555997 [15] Scatchard G. Excess free energy and related properties of solutions containing electrolytes[J]. Journal of the American Chemical Society, 1968, 90(12):3124~3127. doi: 10.1021/ja01014a027 [16] Pitzer K S, Mayorga G. Thermodynamics of electrolytes. Ⅱ. activity and osmotic coefficients for strong electrolytes with One or both ions univalent[J]. Pitzer, 1973, 77(19):2300~2308. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ026170697/ [17] Pitzer K S, Mayorga G. Thermodynamics of electrolytes. Ⅲ. activity and osmotic coefficients for 2~2 electrolytes[J]. Journal of Solution Chemistry, 1974, 3(7):539~546. doi: 10.1007/BF00648138 [18] Hamer W J, Wu Y C. Osmotic coefficients and mean activity coefficients of uni-univalent electrolytes in water at 25℃[J]. Journal of Physical and Chemical Reference Data, 1972, 1(4):1047~1099. doi: 10.1063/1.3253108 [19] Clarke E C W, Glew D N. Evaluation of the thermodynamic functions for aqueous sodium chloride from equilibrium and calorimetric measurements below 154℃[J]. Journal of Physical and Chemical Reference Data, 1985, 14(2):489~610. doi: 10.1063/1.555730 [20] Guendouzi M E L, Mounir A, Dinane A. Water activity, osmotic and activity coefficients of aqueous solutions of Li2SO4, Na2SO4, K2SO4, (NH4)2SO4, MgSO4, MnSO4, NiSO4, CuSO4, and ZnSO4 at T=298.15 K[J]. The Journal of Chemical Thermodynamics, 2003, 35(2):209~220. doi: 10.1016/S0021-9614(02)00315-4 [21] Castellanos E, Villar M V, Romero E, et al. Chemical impact on the hydro-mechanical behaviour of high-density FEBEX bentonite[J]. Physics and Chemistry of the Earth, 2008, 33(S1):S516-S526. https://www.researchgate.net/publication/223791226_Chemical_impact_on_the_hydro-mechanical_behavior_of_high-density_FEBEX_bentonite [22] Pitzer K S, Kim J J. Thermodynamics of electrolytes. Ⅳ. Activity and osmotic coefficients for mixed electrolytes[J]. Journal of the American Chemical Society, 1974, 96(18):5701~5707. doi: 10.1021/ja00825a004 [23] Dinane A, El Guendouzi M, Mounir A. Hygrometric determination of water activities, osmotic and activity coefficients of (NaCl+KCl)(aq) at T=298.15 K[J]. The Journal of Chemical Thermodynamics, 2002, 34(4):423~441. doi: 10.1006/jcht.2001.0845 -
下载:
下载:
计量
- 文章访问数: 479
- HTML全文浏览量: 141
- PDF下载量: 26
- 被引次数: 0
