渭河咸阳段非饱和层状沉积物中水分分布特征

常琛朝，程东会,2†，钱康

(1.长安大学，环境科学与工程学院，710054，西安；2.长安大学 旱区地下水文与生态效应教育部重点实验室，710054，西安)

渭河咸阳段非饱和层状沉积物中水分分布特征

常琛朝1，程东会1,2†，钱康1

(1.长安大学，环境科学与工程学院，710054，西安；2.长安大学 旱区地下水文与生态效应教育部重点实验室，710054，西安)

层状结构是河流相松散沉积物中最常见的沉积结构，通常表现出一些特有的水力学性质，对水分运动过程及沉积物的持水性能产生影响。本文通过对渭河陕西咸阳段河漫滩一个层状沉积物剖面的含水率分布特征进行原位观测，讨论层状沉积物中含水率的分布特征，及其与粒度和毛细壁垒效应的关系。结果表明：在以水平层理和交错层理为主的粗细相间的层状河流沉积物中，天然含水率也呈干润相间的层状分布特征，各层含水率的大小主要与粒度成负相关；沉积物粗细相间的层状分布为毛细壁垒的形成提供条件，在降雨入渗和土壤水蒸发过程中，毛细壁垒效应对水分运动的控制会导致细粒层中含水率高于粗粒层。层状松散沉积物中水分的这种分布规律，为理解土壤持水性能提供了新的思路，对研究土壤含水率与降雨入渗、蒸发等水文过程的关系具有重要意义。

河流沉积物； 层状结构； 含水率； 毛细壁垒

层状结构是河流相松散沉积物中最常见的沉积结构，也广泛分布于冰川、风和湖泊等地质营力作用下的沉积物中，是沉积物非均质的一种重要形式。这种类型的沉积物表现出一些特有的水力学性质，如渗透系数各向异性，介面效应和毛细壁垒效应等。这些特性对水分的入渗和蒸发，污染物的迁移和转化以及水土保持等方面均有重要影响[1]。

层状介质的入渗与均质介质入渗有很大差异，持水能力也显著不同，因此，层状介质中复杂的水分入渗运动受到水文地质及环境研究者的广泛关注。层状介质中水分运动存在一种被称为毛细壁垒效应的现象[2]，即在非饱和层状沉积物中，当细粒介质层上覆于相对较粗的介质层时，粗粒层阻碍水分向下运动，而使细粒层持有更多的水分[2-3]。只有当水分在细粒层积累使界面处吸力达到某一特定值，通常认为达到下伏粗粒层的进水值，水分才会开始进入下部粗粒层[4]。由于毛细壁垒效应阻碍水分下渗，因此会增加上覆细粒层的含水率[5-6]。毛细壁垒效应的强弱以上覆细粒层持水能力增量大小为标志，在相同水分入渗条件下，相邻2层介质粒径差异越明显，则上覆细粒层持水能力增量越大，毛细壁垒效应越强[7-8]。Stormont等[9]的试验也得出同样结论，并用定量方法计算出毛细壁垒效应使细粒介质含水率的增加量；但是，目前毛细壁垒效应的大部分研究是室内简单层状介质入渗实验取得的，对于野外天然层状沉积物的研究很少。

笔者通过一个天然河流相非饱和松散层状沉积物中含水率、介质粒度、密度等性质的原位测定，研究天然状态下层状松散沉积物中水分分布特征及其与各层沉积物粒度的关系，在此基础上分析天然层状沉积物中毛细壁垒效应对水分分布的控制作用。

1 研究区概况

试验场地位于渭河陕西咸阳段的一处河漫滩上，地理位置为 E 108°49′47″，N 34°22′39″。该地区年平均气温13.2 ℃，年均降雨量540 mm左右，60%左右的降雨发生在7、8和9月[10]，而年均潜在蒸发量500 mm左右，属于典型的半干旱地区[11]。渭河在该地发育宽广的河漫滩，均具有明显的层状沉积特征，滩面较平坦，表层为粉砂,下为含砾石的细、中砂。河漫滩开垦指数高，以建造堤坝和河滩公园为主[12]。河漫滩地下水位埋深<4 m[13]。

2 材料与方法

实测剖面是在该河漫滩一个人工挖掘的陡坎基础上改造而成的，剖面保留完整且层理特征明显，具有典型的河流沉积特征。沿陡坎坡面向前挖掘4.0 m左右形成一个高度为1.42 m的垂直沉积物剖面，以避免陡坎坡面的侧壁效应对剖面测量产生影响。剖面上水平层理和交错层理构成典型的河漫滩二元沉积结构。根据肉眼可辨识的层面，将剖面上的沉积物分为11个层系，除上部扰动层外，从上到下按英文字母顺序命名，依次从A到K(图1)。同一层中不同取样点用该层字母加数字下标进行标记。最底部K层是以粗砂为主的水平层理，其上部J层和I层中有明显的“上截下斜”的交错层理特征，层系界面为平面且彼此不平行。G～D层具有近似水平的层理，层系界面呈曲面且纹层与之一致。剖面上部的A、B和C层为水平层理，层中的层理产状与层面产状基本相同。有些层会形成尖灭，如G层左侧和J层。

黑圆点为含水率测量和取样的位置，在三角形标记的位置除含水率测量和取样外还进行了吸力测量。图中的大写字母和数字下标为测量和取样位置编号。图中铁锹长度为1 m。Black dots denote positions of water content measurement and sampling,and the suction was measured additionally on the position denoted by triangle; the capital letters and numbers in subscripts denote the numbering of measurement and sampling locations.The length of shovel is 1 m. 图1 研究剖面的沉积结构Fig.1 Layer sedimentary structure on the profile

2.1 试验设计

在所选剖面上测量A～K层各取样点天然含水率、粒度分布和密度。另外，J层与上部 I层及下部K层的粒度差异明显，具备形成毛细壁垒的条件，且J层厚度较大，便于测量探头的安装；因此在I层(I1和I2)、K层(K3)和J层(J1、J2、J3、J4、J5和J6)进行沉积物的吸力测量。

2.2 试验方法及数据来源

2015年12月4日在渭河滩进行试验，此前至少20 d左右没有降雨或降雪过程。含水率测量采用土壤含水率传感器(EC-5,Decagon Devinces Inc,USA)，测量精度为0.001 cm3/cm3，测量误差为±0.01 cm3/cm3。吸力测量采用土壤水势传感器(MPS-6,Decagon Devinces Inc,USA)，测量精度为0.1 kPa，测量误差为±2 kPa。含水率和吸力测量数据均采用在线采集和储存系统(EM50,Eecagon Devinces Inc,USA)。在野外测量前72 h，先在室内将水势传感器预埋在含水率0.05 cm3/cm3左右(与野外介质的含水率相当)的细砂中并达到平衡，这样大大减少野外测量时水势传感器与介质之间的平衡时间。在野外原位测量吸力时，13 h即达到平衡。为方便讨论，笔者将吸力单位转化为cmH2O(1 kPa=10.2 cmH2O)。野外测量时，设置的吸力测量频率为1次/20 min，在野外连续测量15 h，取平衡后吸力实测值的平均值作为该测量点的吸力值。在吸力测量结束时，进行含水率原位测量，测量频率为每2 min1次，每个测量点连续测量10 min，取平衡后含水率随时间无变化的值作为该取样点实测含水率。

密度采用重量法测定。在测量位置上用体积为98.18 cm3的环刀取原状样品，每个位置取2个共40个样品，装袋密封，在实验室烘干后称量求其密度。另外，利用筛分法，将不同筛孔直径的标准筛固定于振筛机上，筛分完成后通过称量得到样品的粒度分布。

3 结果与分析

3.1 粒度和密度分布

研究剖面各层沉积物中值粒径(d50)在0.16～0.39 mm之间，总体上各层介质呈现粗细相间的特征。粒度相差较大的相邻层分别为D层(粗)和E层(细)、G层(细)和H层(粗)、J层(细)和K层(粗)(图2)。剖面中最粗层为H层，最细层位于J层底部。

观测剖面上密度变化范围介于1.32～1.53 g/cm3之间。与相邻层中值粒径大小相间相对应，各层密度也存在大小相间分布的特征；但趋势与粒径变化趋势相反，粗粒层密度较小，而与其相邻的细粒层密度相对较大(图2)。

3.2 吸力特征

所测9个采样点(图1)的吸力范围为98～122 cm H2O，相邻层的吸力值相差并不大。9个吸力测量点实测吸力与含水率如表1。介质的含水率与吸力的关系可用水土特征曲线描述，不同介质的水土特征曲线在同一吸力下对应不同的含水率，如采样点J1和J6，吸力相等，含水率却不相同，分别为0.012 cm3/cm3和0.110 cm3/cm3。

3.3 含水率分布

刚挖掘出的新鲜剖面表观呈现湿润特征，含水率分层现象不明显。经过不到30 min的风干后，剖面表层可以观察到明显的颜色差异，含水率测定结果表明深色层(A、C、E、G和I)的含水率均高于相邻的浅色层(B、D、F、H和K)(图2)。剖面图像整体呈现干润相间的特征(图1)。在剖面20个取样点中，I层左侧取样点I1含水率最高，为0.148 cm3/cm3，而最底部K层3个取样点含水率最低，均为0.01 cm3/cm3。相邻2层含水率差别较大的为取样点I1和J1，含水率相差0.136 cm3/cm3，J4和K2相差0.117 cm3/cm3，J5和J6相差0.1 cm3/cm3。

A-K are corresponding to the points in Fig.1.图2 沉积物剖面各层含水率、中值粒径和密度Fig.2 Water content,median grain size and bulk density in each layer of the sediments profile

表1 剖面中I、J和K层实测吸力值与实测含水率Tab.1 Measured suction and water content for the samples in the layer I,J and K

注：以上3组取样点分别对应图1中3组垂向吸力实测点。Note:3 sampling points belong to 1 group,and each of them corresponds to actual measuring point denoted by triangle in Fig.1.

图2中含水率与中值粒径(d50)随土层深度变化趋势相反，为更加定量的描述二者的关系，对剖面20个取样点(图1)所测得的含水率和d50数据进行统计分析，可得含水率大小与d50成负相关关系(图3)，密度该规律本质上与土壤非饱和理论一致，反映含水率和介质粒度的相互关系。在天然层状沉积物中，假设沉积物各层水分运动已达平衡，即总水头(位置水头与吸力水头之和)相等，相邻2层界面附近位置水头忽略不计，即认为相邻层吸力水头相等。由水土特征曲线可得，同一吸力下，细介质含水率必然高于粗介质含水率，也即处于平衡状态时的层状沉积物中，相邻2层细粒层含水率大，粗粒层含水率小，这就是相邻层含水率差别较大的原因之一。

图3 样品中值粒径与含水率的关系Fig.3 Relationship between median grain size and corresponding water content for all the samples

即使在同一层内，由于沉积物粒径的变化也会引起含水率的变化。取样点J1、J2和J4位于同一层面，从图4样品的粒度分布曲线可以看出，3个样品由细到粗为J4、J2、J1，含水率依次为0.127、0.061、0.012 cm3/cm3，介质越细，含水率越高。

图4 J1、J2和J4取样点处样品的粒度分布曲线Fig.4 Cumulative grainsize distribution curve of three sediment samples (J1,J2 and J4)

4 讨论

4.1 “毛细壁垒”效应对水分分布的控制作用

层状沉积物中各层含水率受介质粒度的影响外，层状结构本身也是控制水分分布的另一重要因素。天然状态下，上细下粗的层状结构为毛细壁垒的形成提供了条件。在水分入渗时，水分不易穿越下伏粗粒层下渗，而是在上覆细层积累[4]。这时候的下伏粗粒层充当“毛细壁垒”，阻碍水分的流动[14]。Iqbal[7]的研究指出，相邻2层粒度差异越明显，毛细壁垒效应越强，上覆细粒层持水量越多。

毛细壁垒效应对天然沉积物水分分布的控制可通过相同粒度层含水率不同来论述。取样点中的J2和G具有相同的粒度分布，但含水率却不同，分别为0.06和0.08 cm3/cm3。它们在各自的毛细壁垒结构中是位于上部的细层，下伏粗层分别是K1和H，在入渗情况下，水分在细层J2中积累，其吸力随含水率的增加而减小，达到K1的进水值时，水分才会开始进入K1层[4]。由中值粒径可以判断K1比H更细(表2)，因此K1的进水值更高。与G相比J2可以更早的达到相应下层的进水值，壁垒效果更弱，水分在J2层的积累量更少，含水率更低。一般情况下，在上细下粗结构中，如果上覆细粒层为同种介质，对应于下伏粗粒层介质越细的含水率越低，这与Stormont等[8]通过土柱试验研究毛细壁垒效应的结果一致。

表2 取样点J2、K1、G、 H处的含水率与d50Tab.2 Water content and median grain size of four sediment samples (J2,K1,G and H)

J1和J3均位于毛细壁垒的粗层，与之相对的上覆细粒层分别为I1和I2。取样点J1和J3粒径分布相同，含水率大小却不相等(表3)，在入渗情况下，水分会先在上覆细层I1和I2中积累。由于I1比I2更细，因此水分在I1层积累达到比I2层更高的含水率，才能使吸力等于下伏粗层的进水值(下伏粗层J1和J3的进水值相等)，水分才会开始下渗。对于J1和I1这一毛细壁垒结构，水分更难进入下伏粗层J1，因此与J3相比，J1含水率更低。一般情况下，在上细下粗结构中，如果下伏粗层为同种介质，上覆细粒层越细，所对应的粗粒层含水率越低。

表3 取样点I1、 J1、I2、J3处的含水率与d50Tab.3 Water content and median grain size of four sediment samples (I1,J1,I2 and J3)

用粒度与毛细壁垒效应共同解释天然层状沉积物中水分分布，与实际测量结果吻合。在天然层状沉积物中，各层沉积物粒度大小和毛细壁垒效应共同控制含水率分布。

4.2 毛细壁垒结构中粗粒层和细粒层的相对性

同一种介质，即可充当毛细壁垒结构中的粗粒层，也可在另一毛细壁垒中充当细粒层。如图5所示，取样点F和I3的概率累积曲线基本重合，可视为相同粒径分布的介质，F层在E层和F层构成的毛细壁垒中为下伏粗层，I3层在I3层和J5层构成的毛细壁垒结构中为上覆细层；因此，粒径分布相同的2层，由于在毛细壁垒结构中充当的角色不同，含水率也会不同。

图5 E、I3、F和 J5取样点样品的粒度分布曲线Fig.5 Cumulative grainsize distribution curve of sediment samples at E,I3,F and J5

4.3 不同水文过程的毛细壁垒效应分析

野外层状沉积物中水分运动主要与降雨入渗和土壤水蒸发这2个水文过程有关。天然水分分布是在不同水文过程中，在沉积物粒度及层状结构的影响下，水分运动达到平衡后的结果。

已有室内实验通过模拟降雨来观察简单层状沉积物中降雨强度和时间对毛细壁垒效应的影响[15]。基于此，我们进一步分析降雨时野外复杂层状沉积物中毛细壁垒效应对水分运动的影响。降雨入渗时，由于沉积物各层粒度粗细相间的特征，一定存在相邻2层上细下粗的结构。当水分从上覆细粒层入渗时，下伏粗粒层形成水分运动的壁垒，使水分在细层积累。界面处吸力水头达到粗粒层的进水值时，水分才会开始进入粗层。当水流至下一个上细下粗结构，也会出现相同的水分入渗特点。这与Zornberg等[6]对降雨时毛细壁垒效应的描述一致。虽然他用细粒层和粗粒层渗透系数的不同来解释下伏粗粒层对水分的阻碍作用，但本质上与进水值的原理一致。野外层状沉积物中水分天然分布正是入渗平衡后的结果，水力学性质不同的两层形成的毛细壁垒效应，阻碍水分的重分布过程，使得细粒层出现较高的含水率[16]。

沉积物蒸发过程为水分向上运动的过程，此时上粗下细的层状结构是毛细壁垒形成的基础。已有室内研究发现，上砂下黏型层状介质的土壤累积蒸发量小于上黏下砂型层状土壤[17]，也有研究者通过土柱实验发现土柱中粗介质夹层存在可显著减少蒸发量，主要原因为上粗下细层状结构中的毛细壁垒效应抑制了水分的蒸发作用[18-19]。在野外层状沉积物的水分蒸发过程中，水分从下伏细粒层开始向上运动，上覆粗粒层作为水分运动的壁垒，抑制了水分的蒸发，且蒸发过程散失的水分主要来源于粗粒介质[17]。在一定程度上，野外层状沉积物粗细相间的结构能够减少水分的蒸发，提高沉积物的持水性。

5 结论

渭河咸阳段河流松散沉积物具有明显的粗细相间的层状特征，水分分布也呈现干层与湿层相间的层状分布特征，各层含水率的大小与粒度成负相关。另外，沉积物粗细相间的层状特征为毛细壁垒的形成提供条件，使水分分布也受到毛细壁垒效应的控制。在降雨入渗和土壤水蒸发过程中，毛细壁垒对水分运动的控制会导致细粒层持有更多的水分。降雨入渗时，毛细壁垒会阻碍水分入渗，而蒸发过程中，毛细壁垒也会阻碍蒸发。总体上层状松散沉积物中水分的分布，受粒度和沉积结构二者共同控制。层状松散沉积物中这种水分分布规律，为垃圾填埋场防渗系统，土壤再利用及水土保持工程的设计提供了新的思路。

[1] 程东会，钱康，常琛朝.水分入渗条件下层状非饱和松散沉积物中的毛细壁垒效应研究[J].水文地质工程地质,2016,43 (4):20.CHENG Donghui,QIAN Kang,CHANG Chenchao.Capillary barrier effects during the infiltration process in layered unsaturated unconsolidated sediments [J].Hydrogeology & Engineering Geology,2016,43 (4):20.

[2] STORMONT C J.The effect of constant anisotropy on capillary barrier performance [J].Water Resources Research,1995,31(3):783.

[3] NAHLAWI H,BOUAZZA A,KODIKARA J.Surface water infiltration in a 1-dimensional soil-geotextile column [C]∥Proceedings of the 10th Australia New Zealand Conference on Geomechanics,Brisbane,Australia.2007:21.

[4] HILLEL D,BAKER R S.A descriptive theory of fingering during infiltration into layered soils[J].Soil Science,1988,146(146):207.

[5] MCCARTNEY J S,ZORNBERG J G.Effects of infiltration and evaporation on geosynthetic capillary barrier performance [J].Canadian Geotechnical Journal,2010,47(11):1201.

[6] ZORNBERG J G,BOUAZZA A,MCCARTNEY J S.Geosynthetic capillary barriers:Current state of knowledge [J].Geosynthetics International,2010,17(5):273.

[7] IQBAL M Z.Effects of layered heterogeneity in subsurface geologic materials on solute transport under field conditions:A case study from northeastern Iowa,USA [J].Hydrogeology Journal,2000,8(3):257.

[8] STORMONT J C,ANDERSON C E.Capillary barrier effect from underlying coarser soil layer [J].Journal of Geotechnical & Geoenvironmental Engineering,1999,125(8):641.

[9] STORMONT J C,MORRIS C E.Method to estimate water storage capacity of capillary barriers [J].Journal of Geotechnical & Geoenvironmental Engineering,1999,125(10):918.

[10] 周晓红，赵景波，王国祥.渭河高陵段河漫滩沉积物指示的洪水演变[J].干旱地区农业研究,2007,25(5):226.ZHOU Xiaohong,ZHAO Jingbo,WANG Guoxiang.The sediment of floodplain instruction flood change in Gaoling part of the Weihe River [J].Agricultural Research in the Arid Areas,2007,25(5):226.

[11] 刘闻.基于模型的水文模拟及径流响应分析：以渭河流域关中段为例[D].西安：西北大学,2014：10.LIU Wen.Hydrological simulation and runoff response analysis based on SWAT Model:a case study of Weihe River in Guanzhong area [D].Xi′an:Northwest University,2014:10.

[12] 牛俊杰，赵景波，马莉，等.西安北郊渭河河漫滩沉积与洪水事件[J].地理研究,2010,29(8):1484.NIU Junjie,ZHAO Jingbo,MA Li,et al.The research on the sediment of floodplain and flood events of Weihe River in the northern suburbs of Xi′an [J].Geographical Research,2010,29(8):1484.

[13] 董艳慧，周维博，卜卿，等.RBF网络在西安渭滨地下水位埋深预测中的应用[J].节水灌溉,2012(12):66.DONG Yanhui,ZHOU Weibo,BU Qing,et al.Aplication of RBF neural network to forecast depth of groundwater level in Weibin River wellhead [J].Water Saving Irrigation,2012(12):66.

[14] HANKS R J,BOWERS S A.Numerical solution of the moisture flow equation for infiltration into layered soils[J].Soil Science Society of America Journal,1962,26(6):530.

[15] MANCARELLA D,SIMEONE V.Capillary barrier effects in unsaturated layered soils,with special reference to the pyroclastic veneer of the Pizzo d′Alvano,Campania,Italy[J].Bulletin of Engineering Geology and the Environment,2012,71(4):791.

[16] ZETTL J D,BARBOUR S L,HUANG Mingbin,et al.Influence of textural layering on field capacity of coarse soils[J].Canadian Journal of Soil Science,2011,91(2):133.

[17] HUANG Mingbin,BRUCH P G,BARBOUR S L.Evaporation and water redistribution in layered unsaturated soil profiles[J].Vadose Zone Journal,2013,12(1).doi:10.2136/vzj2012.0108.

[18] SHOKRI N,LEHMANN P,OR D.Evaporation from layered porous media[J].Journal of Geophysical Research Atmospheres,2010,115(B6):258.

[19] 任利东,黄明斌.砂性层状土柱蒸发过程实验与数值模拟[J].土壤学报,2014(6):1282.REN Lidong,HUANG Mingbin.Experiment and numerical simulation of soil evaporation from layered sandy soil columns[J].Acta Pedologica Sinica,2014(6):1282.

Water content distribution of unsaturated layered sediments of Weihe River in Xianyang section

CHANG Chenchao1,CHENG Donghui1,2,QIAN Kang1

(1.School of Environment Science and Engineering,Chang′an University,710054,Xi′an,China; 2.Key Laboratory of Subsurface Hydrology and Ecological Effects in Arid Region (Chang′an University),Ministry of Education,710054,Xi′an,China)

[Background] Layer structure is the most common sedimentary structure in alluvial unconsolidated sediments,which usually exhibits special hydraulic properties to influence water flow and water retention ability.In this paper,a situ investigation on the water content distribution,grain-size and bulk density for a layered sediment profile was conducted at the floodplain of Weihe River in Xianyang,in the center of Shaanxi,then the effects of the grain-size and the capillary barriers on natural water content were discussed.[Methods] The experiment site outcrop was mapped into 11 macro layers having different texture and structures.With each layer,well developed laminations were observed.We investigated the characteristics of water content and its corresponding influence factors including soil bulk density,grain-size and matric suction in each layer,based on sampling analysis and the measurement of EC-5 Small Soil Moisture Sensor and MPS-6 Calibrated Water Potential Sensor.Firstly,we analyzed the relationship between grain-size and water content in each layer,based on this,the effect of capillary barriers on water content was discussed.Then,we analyzed the influence of capillary barriers in the process of precipitation infiltration and soil water evaporation.[Results] The results illustrated that in the layered alluvial sediments with alternating layers of coarse and fine sands,natural water content exhibited in a layered distribution on a form of interbedding high and low water content.The magnitude of the water content was negatively correlated with the grain-size,and was positively correlated with the bulk density.The water content of site outcrop was between 0.010 cm3/cm3-0.148 cm3/cm3,the median size ranged from 0.16 mm to 0.39 mm and the bulk density ranged from 1.32 g/cm3to 1.53 g/cm3.The results also demonstrated that the layered sediments with alternating layers of coarse and fine sands formed the capillary barriers in the process of water flow,therefore by which the water content distribution was controlled.In both process of precipitation infiltration and soil water evaporation,the capillary barriers effects will lead to the water content in fine layer obviously higher than that in coarse layer.[Conculsions] In conclusion,capillary barriers become a significant factor to control water flow in unsaturated soils.At the process of precipitation infiltration,capillary barriers have the effect on permeability reduction,which can be used in designing landfill anti-seepage system,optimized drainage layers and so on.Additionally,insights are potentially useful for designing mulching strategies and capillary barriers aimed at reducing evaporative losses,increasing water retention ability of soil,which is significantly important in arid and semi-arid regions.

alluvial sediments; layer structure; water content; capillary barriers

2016-09-06

2017-06-17

项目名称：中央高校基本科研业务费专项资金“点状入渗条件下层状沉积物中水分运动的毛细壁垒和穿透效应”(310829162015)

常琛朝(1991—)，女，硕士研究生。主要研究方向：地下水环境。E-mail:452293490@qq.com

†通信作者简介：程东会(1969—)，男，教授。主要研究方向：地下水资源与环境。E-mail:chdhbsh@chd.edu.cn

S152.7

A

2096-2673(2017)04-0104-07

10.16843/j.sswc.2017.04.013