充填复垦土壤水分竖直运动模拟研究.pdf

第3 6 卷第5 期 中国矿业大学学报 V o l 。3 6N o 。5 2 0 0 7 年9 月 J o u r n a lo fC h i n aU n i v e r s i t yo fM i n i n g ＆T e c h n o l o g yS e p ．2 0 0 7 文耄缡号1 0 0 0 1 9 6 莲 2 0 0 7 G 5 0 6 9 0 0 6 充填复垦土壤水分竖直运动模拟研究 王 辉1 ，韩宝平1 ，卞正富2 1 。巾鏊矿受大喾繇境与测绘学貔，鏊苏徐簿2 2 1 1 1 6 } 2 。审嚣矿鲎大学瓣学技术姓，江苏豫翔2 2 1 1 1 6 摘要应用直藏柱状模擞模拟研究了煤矸石和粉煤灰两种充填模式下覆土层土壤水分运动特 点，耙充填复爨土壤水分运动过程分为3 令泠段争4 令层次，并遂过与全部充填主壤情琵下水分 动态对比找出其差异性．把监测数据与应用蛰直一维动力学模型计算得到的水分分布规律比较， 找出了充填层时覆土层土壤水分动态产生影响的本质原因，并对充填煤砰菇和粉煤灰复垦的土 壤水分动力攀参数进镗修正．结果表暖充填煤矸石争粉煤灰复塾骑主壤水分动态与点常主壤答 同源于煤矸石和粉澡灰的颗粒级配与土壤的麓异；覆土层土壤与充填料孔隙的差别大，其保持水 分的能力反而较强，煤矸石充填不会导致覆土层水分大量渗漏． 关键词充填复垦；土壤水分；鳖煮运番；颗粒级配 ’ 审圈分类号T D8 8 文献标识码A S i m u l a t i o nR e s e a r c ho nW a t e rU p r i g h tM o v e m e n ti n R e c l a i m e dS o i lb yF i l l i n g W A N GH u i l ，H A NB a o - p i n 9 1 ，B I A NZ h e n g f u 2 1 ．S c h o o lo fE n v i r o n m e n ta n dS p a t i a lI n f o r m a t i c s ，C h i n aU n i v e r s i t yo fM i n i n g ＆T e c h n o l o g y ， X u z h o u 。J i a n g s u2 21 I 16 ，C h i n a2 ．S c i e n c ea n dT e c h n o l o g yD e p a r t m e n t ， C h i n aU n i v e r s i t yo fM i n i n g ＆T e c h n o l o g y 。X u z h o u 。J i a n g s u2 2 1 1 1 6 ，C h i n a A b s t r a c t W a t e rm o v e m e n tc h a r a c t e r i s t i c si nc o v e rs o i li nt w ot y p e so ff i l l i n gg a n g u ea n df l ya s h w e r es t u d i e du s i n gu p r i g h tc o l u m ne n t i t ym o d e l ，a n dt h e ns o i lw a t e rm o v e m e n tw a sd i v i d e di n t o 3p h a s e sa n d4l a y e r s ．T h ed i f f e r e n c e so fw a t e rd y n a m i c sw e r ef o u n do u tb yc o n t r a s t i n gw i t h t h a to fa 1 1f i l l e dw i t hs o i l ．T h en a t u r a lc a u s eo fi m p a c t so nc o v e rs o i ld y n a m i c sb yf i l l i n g sw a s f o u n do u tb yc o n t r a s t i n gm o n i t o r i n gd a t aw i t ht h a tg o t t e nb ys i m u l a t i o nu s i n go n ed i m e n s i o n a l m a t h e m a t i cm o d e lo fs o i iw a t e r ，a n dt h e ns o i lw a t e rd y n a m i cm o d e l so ff i l l i n g sw e r ea m e n d e d i nt W Oc o n d i t i o n s ．T h er e s u l ti n d i c a t e st h a t t h ed i f f e r e n c e so fw a t e rd y n a m i c si nc o v e rs o i li n f i l l i n g sa r ec a u s e db yt h er e l a t i o n so fp a r t i c l e sm a k e u pb e t w e e nt h ec o v e rs o i la n dt h ef i l l i n g s ， j u s tl i k eg a n g u ea n df l ya s h ；t h eb i g g e rd i f f e r e n c eo ft h ep o r o s i t yb e t w e e nc o v e rs o i la n df i l l i n g s ，t h eg r e a t e rc a p a b i l i t yo ft h ec o v e rs o i lk e e p st h ew a t e r ．S oi t c a nn o tb el a r g ew a t e rl e a k a g ei nl a n df i l l e dw i t hg a n g u e 。 。 K e yw o r d s r e c l a m a t i o nb yf i l l i n g s o i lw a t e r u p r i g h tm o v e m e n t ；p a r t i c l e sm a k e u p 作物生长与土壤水分动态有密切关系，而土壤 翻瑟构登、土壤质地分层等性质都会对土壤承分运 动产生影响。充填复垦改变了土壤的削面类型，必 然会影响水分的运动。如何改善充填复垦地的承 收稿日期2 0 0 6 1 0 一2 6 基金璜瓣国家鑫然瓣学基金暖终 5 0 5 7 4 0 9 5 } 审国矿韭大学校毒年基金壤器 O P 0 6 1 0 2 1 作者麓介王 簿 1 9 7 6 一 ，男，扛荔嚣锅出县入，漭薅，硕士，筑事鲍复黧与生态重建努露鹩醑究． E m m i l w a n g h u e i c u m t ．e d u ．c n T e l 0 5 1 6 - 8 2 1 7 3 8 2 8 万方数据 第5 期 王辉等充填复垦土壤水分竖直运动模拟研究6 9 1 分条件，是提高充填复垦土地生产力的重要问题之 一● 根据江苏省徐州市柳新国家复垦示范区的充 填复垦试验，煤矸石和粉煤灰充填复垦后作物生长 不良，产量不稳，其主要原因是水分条件差，充填物 料的性质或充填行为改变了土壤的水分条件，因 此，可以通过充填复垦土壤水分特征来分析．本文 通过柱状实物模型模拟对比，研究了充填煤矸石和 粉煤灰在覆土层土壤水分动态与充填土壤情况下 水分动态的差异．并把监测数据与应用竖直一维动 力学模型计算得到的水分分布规律相比较，找出充 填层对覆土层土壤水分动态产生影响的本质原因． 1 研究方法 ⋯铡奏赫 a b c 图1不同充填模式下柱状模型 F i g ．1 S c h e m a t i cr e p r e s e n t a t i o no fc o l u m nm o d e l s w i t hd i f f e r e n tf i l l i n g s 覆土层土壤为粉土 中国土壤分类制 ，其颗粒 级配状况见表1 ．为了反映田间土壤的状态，6 0c m 覆土层分两层装填，3 0 ～6 0c m 土层设计容重1 ．3 5 g ／c m 3 ，0 - - - ．3 0c m 土层设计容重1 ．2g ／c m 3 ．通过装 入土壤的质量和柱体体积来控制容重和紧实度． 表1土壤颗粒分析 T a b l e ]S o i lp a r t i c l e sc l a s s i f i c a t i o n 充填层分别为煤矸石、粉煤灰和泥土．充填层 为煤矸石的两组充填柱，充填的煤矸石要求粒度大 小相近 ≤2c m ，并尽量压实，通过充填相同的体 积和相同的充填量，保证有相同的压实度．充填层 为粉煤灰的两组充填柱，粉煤灰尽量压实，并使充 填体积和充填量一致确保有相同的压实度．充填层 为泥土的两组充填柱，土壤质地及颗粒级配与覆土 相同，充填的泥土尽量压实，充填的体积和充填量 相同 根据充填量计算容重为1 ．4g ／c m 3 ． 为了避免覆土层土壤向下运动充斥充填物间 隙，充填层与覆土层之间铺设双层纱布．充填柱底 端用纱布封托．覆土层及充填层基本情况见表2 ． 襄2充填层及覆土层基本情况 T a b l e 2B a s i cs i t u a t i o no ft h ef i l l i n g s 柱体上部留5c m 用于灌水．水分监测之前浸 泡柱体使土壤含水量达到饱和，最后一次灌入等量 水后用枯草覆盖，覆盖厚度5c m ，以减少蒸发．从 第2 天起，以后每3d 用T D R 测定一次设计深度 的土壤含水量． 2 数据与分析 2 ．1 覆土层含水量垂直变化过程的差异 煤矸石、粉煤灰充填和泥土充填的柱体含水量 随深度和时间的变化过程见图2 ．由图2 可以看 出，3 种情况都是从3 0c m 处往下土壤含水量较上 层开始增加，3 5c m 处基本达到最大值．这是因为 土壤容重增加，土壤总孔隙度减小的缘故．曲线随 时间向低含水量方向移动，表明土壤剖面含水量随 时间不断降低．图形低含水量一侧曲线容重大，表 明随时间推移，土壤剖面含水量降低的速率在不断 减小． 但重要的是三者变化趋势并不完全相同，煤矸 石充填情况下，土壤剖面含水量变幅大，3 5c m 以 下水含量迅速下降，并在5 0c m 处有所反弹；而粉 煤灰和泥土充填的情况，3 5c m 以下呈缓慢下降趋 势．5 5e l D _ 处各点的分布也有很大不同，尤其是煤 矸石充填情况，分散且交叉多，前3d 含水量一直 万方数据 6 9 2中国矿业大学学报第3 6 卷 处于较高水平，而9d 之后妒 H O 在2 5 ％以下． 粉煤灰和泥土充填的情况，含水量下降的速率较稳 定，荔粉煤灰充填的情况下，5 5c m 处含水量较低， 泥土充填的较高．这些差异决定于两种介质之间孔 隙大小差异。煤矸石的孔隙逡远大于上覆的泥土， 当水分下渗至交界獗时，水分并不会快速渗漏，丽 4 5 4 0 萎3 5 圣3 0 匿2 5 2 0 1 5 是会先聚积，然后再迅速渗漏，这造成了5 0a m 处 含水爨增加，5 5c m 处含水量变化大．粉煤灰的孔 豫与所磊酶主壤狸当略大，承分入渗瑚界面层时积 聚的过程非常缀而会迅速下渗，导致5 5c m 处含水 量一直稳定较低． 叫月2 3 嗣叫置拳呈3 日墨 妻葱毫憾。量 1 0 1 52 02 5 3 0 3 54 04 55 05 5 深度／c m b 粉煤灰充填 c 混士充填 蹰2煤矸蠢、粉煤灰秘泥土充填土壤含水量变化过程 F i g 。2V a r i a b i l i t yo fs o i lw a t e ri nd i f f e r e n tc o l u m nm o d e l sf i l l e dw i t hg a n g u e a s ha n de a r t h 2 ．2 覆土层渗漏爨的变化特征 根据土壤的体积含水量和柱体的体积 柱体内 径3 。5c m ，覆层离5 5c m ，体积为21 1 5 ．6c m 3 ， 可以计算其贮水量及其变化．在监测的2 9d 时间 里，3 种楗体6 0a m 土层的贮水量分别减少了 0 ．4 8 0 ，0 5 4 6 和0 。5 7 4m m 。最然，在番手石充填的情 况下渗漏爨最小，略低于粉煤灰，而在泥土 粉土 充填的情况下，渗漏量却比较大．考虑到水分初期 渗漏速度还与承头压力有关，袭现为与壤初始含 水量有关．3 种柱体的9 H 。0 初始 分别为 3 2 。3 ％，3 2 ．1 ％，3 2 ．6 ％，土壤充填柱体软始贮水羹 高，会导致渗透速度加快．在矸石充填情况下，虽然 水头压力大于粉煤灰充填情况，但是渗漏量却偏 少，这与各自水分渗漏过程的特点以及研究的时闽 段有关．潮于在煤矸石充填情况下，在界面层水分 有积聚和迅速下渗的反复过程，所以研究的时间段 不同，渗漏量的对比关系也不圈，但由予水分的积 聚有阻滞上层水分下渗的作用，总体的渗漏量并不 会最大．如在开始的6d 时间里，煤矸石、粉煤灰充 填和泥土充填的柱体贮水量分别减少T0 。2 1 8 ， 0 ．2 2 7 和0 ．2 5 7m m ，而在最后的1 1d 中，煤矸石、 粉煤灰充填和泥土充填的柱体贮水量分别减少了 0 ．0 9 7 ，0 。1 0 1 纛0 。0 7 7m m ，程最后9d 时闯里，贮 水量分别减少0 ．0 4 2 ，0 ．0 6 3 和0 ．0 7 2m m ． 2 ．3覆土层含水量的分层特征 根据装溅酶主壤承分的变化过程特点可以将 覆土层土壤分成4 个层次，即0 ～1 0 表层 、1 0 ～ 3 0 中间上层 、3 0 ～5 0 中间下层 、5 0 ～6 0c m 临赛层 ．睡雾层主壤戆含承璧受主下介质魏隙的 影响最大，变化特点明显；中间下层土壤容重较上 层犬，含水量一直处于各层最高水平，中间上层含 水量次之．表层受蒸发影响含水量下降趋势明显， 且含承董最低。 但由于充填层性质不同，各层特点出现了较大 差异．整体上粉煤灰充填的情况与泥士充填情况更 耀织。不同充填辩充填情况下，中越层差异不大，趋 势也綦本相同，差异最大的是临界朦土壤含水量， 显然是受到充填层性质的影响．煤研石充填的情况 箍界层含水量变化较大，掰粉煤灰充填秘泥i 充填 含水量变化较平稳．3 种柱体入渗稳定后的7 次数 据变异系数分别为0 。1 1 4 _ | 耋，0 。0 5 53 程0 。0 5 1l 。 2 ．4 覆层含承量变化的阶段性特征 网3 为不同柱体含水魑衰减速度变化曲线，从 第一次监测起，可以把土壤含水量变纯粗略缝分为 3 个阶段，郄0 , - 一9d 为快速递减期，土壤含水罴迅 速下降；9 ～1 8d 为稳定期，土壤含水量的变化平 稳；1 8d 后为超伏期，含水量变讫量略有起伏 这3 个阶段的平均禽水量见表3 ． 4 委， 骛 蔷 惹一 O 时问／d 图3 不同柱体含水量衰减速度变化 F i g ．3V a r i a b i l i t yo fs o i lw a t e rr e d u c i n gr a t e i nd i f f e r e n tf i l l i n g s 但3 种充填情况下各阶段的变化特征并不相 同．煤矸石充填复垦土壤禽水量在快速递减阶段持 续时阀更长，可以达到1 1d 左右，以后两个输段含 水量变化较大，主要是孔性差异引起水分渗漏快慢 万方数据 第5 期王辉等充填复垦土壤水分竖直运动模拟研究6 9 3 交替造成的．粉煤灰充填的覆土层土壤含水量变化 从第二阶段起一直很平稳． 表33 个阶段划分及各段平均土壤含水量 T a b l e3T h r e ep h a s e sm a r k i n go u ta n d a v e r a g ew a t e rc o n t e n t 时段／ d 煤矸石充填 妒‘已尸’7 标准 ％ ”一 粉煤灰充填 差似H 。t ，O ／标准差 雹 ％ 怀。匪茬 泥土充填 妒‘巳尸’7 标准差 ％ ⋯‘4 注i H z O 为H z O 的体积分数的平均值． 3 土壤水分运动方程 土壤水分运动方程能定量描述土壤剖面水分 分布规律，通过把模型计算的数据与模拟监测数据 相对照，找出其规律差别及原因． 由D a r c y 定律得到的R i c h a r d s 方程是一般情 况下土壤水分运动的基本方程，即 c ‰ 警 a z f l 曼- [ K 咖m 警] 掣． 对于垂直一维饱和一非饱和水流，可以改写成 『孕s 。 c “ ] 警 3 厂矿M 、a “] ．a K C m _ lK “ 警l 掣， 式中以为饱和含水率；S 。为比贮水系数；C ‰ 为 容水度，C “ d O ／d 驴m ；当丸 0 时，‰为基质 势；当儿≥0 时，丸为水头；口为气体含水量；愚为导 水率． 土壤水分动力学模型定解条件、求解方法等的 研究比较多，也比较成熟，而确定方程中的水分动 力学参数非常重要．水动力学参数主要包括水分特 征曲线，饱和非饱和导水率以及土壤水分扩散率 等．有关这些参数的测定、计算方法已有许多N - 7 3 ， 且E l 趋完善．但对于土壤水分运动问题，由于存在 空间变异性，由实测方法获得足够准确且具有代表 性的水动力参数很困难．通常用来确定土壤水分运 动参数的方法有多种，近2 0a 来发展起来的估算 土壤水分运动参数的土壤传递函数法 P F D s 应用 较好[ 8 d0 1 ．研究中采用统计模型的系数估计模型来 确定土壤传递函数． 统计实验使用3 2 根相同体积的离心管和1 6 根直径为2e m 的P V C 管进行，使用的土壤经过多 级吸管法分级风干，砂粒粒径取大于0 ．0 5m m ，考 虑到非土壤组成的影响，砂粒含量取0 ．0 5 - - - 0 ．0 8 8 m m 粒径之间，黏粒粒径小于0 ．0 0 1m m ．含水量测 定用烘干法，容重用游标卡尺量取管内土壤长度和 高度确定体积后计算． 3 ．1 驴O 关系确定 土壤水吸力与土壤含水量之间的关系通过土 壤水分特征曲线描述．在实际应用中，土壤的水分 特征曲线由实验测定．本研究土壤水分特征曲线应 用离心机法测定．由于离心机法存在对土壤容重产 生影响的问题[ 1 ，所以在测定时使用较低转速 ≤ 40 0 0r ／m i n 确定两者之间的指数关系．实验表 明，40 0 0r ／r a i n 时，容重的变化仅为0 ．8 5 ％． 通过实测土壤水分动力参数与土壤物理性质 之间进行统计，从而估算出土壤水分运动参数的回 归方程．少汐采用G a r d n e r 的指数关系式 “ A O 一， A 一一e x p l - a b w S a w C 1 c w C 1 2 d w C 1 2 硼 S a h 7 2 ] ， B e f w C 1 2 g w S a 2 w C 1 i 7 2 ， 式中w S a 为砂粒的质量分数；w C 1 为黏粒的 质量分数；口，6 ，f ，d ，e ，f ，g ，h ，i 为拟合参数；y 为土 壤容重． 实验数据经多元逐步统计回归，得出A ，B 中 的各个参数，见表4 ． 表4拟合结果 T a b l e4 C o i n c i d i n gr e s u l t s 土壤有机质在逐步回归中未被引入方程，对A 和B 的影响不显著．根据3 2 个样本的多元回归结 果，回归方程均达到极显著水平，F o ．。。 4 ．2 0 ． 3 ．2K | 的确定 饱和导水率K 。的统计回归参考C o s b y 公式 1 nK 。一a 伽 S a b 叫 C I f y d 0 8 e ， 式中0 l 为土壤有机质． 非饱和导水率K 的计算在饱和导水率的基础 上，与水分特征曲线相结合来确定．导水率采用以 下指数函数表达式计算‘1 2 ] K c ‰，一{ K K 。s 饯“呦J ‰C m ≥ 0 。； 式中口为拟合参数，取0 ．0 2c m ～． 粥啪瞄 乙L 仉 ∞∞∞∞孙孔 Ⅲ啪弼 昭眈弘约拍孔至| 抛啪弘曲惦∞卯筋麓～ 万方数据 6 9 4 中国矿业大学学报 第3 6 卷 实验数据经回归分析得出式中各个参数见表 5 ，F 检验为极显著． 表5参数拟合结果 T a b l e5 C o i n c i d i n gr e s u l t s 拟台参数拟合参数值样本数拟合系数 0 ．0 1 84 2 0 ．1 2 65 0 ．2 9 331 60 ．8 5 0 ．1 4 55 6 ．6 3 65 3 ．3 充填复垦土壤水分运动模型的修正 应用上述模型，将柱体实验中的供试土壤的颗 粒级配、容重及有机质含量代人计算水动力参数， 应用水动力学模型，可以预测柱体水分的垂直一维 分布． 由于研究中覆土层所用土壤性质完全相同，预 测的结果只会有一种．而如前文所述，由于覆土层 水分运动受充填物性质的影响，水分分布并不相 同．由于土柱的土壤完全相同，所以变化的原因与 充填层性质有直接关系，而影响水分运动的本质原 因是充填物的孔性及其与覆土层土壤的孔性的对 比关系． 将模型预测的结果与实测结果比较，表明模型 模拟的含水量变化与泥土充填的情况吻合较好，与 煤矸石充填的中间下层，与粉煤灰充填的中间上层 也有较好的吻合，其它情况差别较大．研究充填物 的颗粒级配发现，充填物粉煤灰小于0 ．0 0 1m m 的 颗粒占4 ．0 ％，而充填的煤矸石粒径最小为2m m ， 小于0 ．0 0 1m m 的颗粒为0 ，土壤小于0 ．0 0 1m m 的颗粒占6 3 ．8 ％．由于充填物对上层土壤水分的 影响与上层土壤本身结构的影响来自不同部分，产 生的影响对于整个覆土层，不是加和式并列关系， 且水分变化曲线的走势基本相同，仅存在大小不同 关系，因此可以对上述参数确定式采用乘积修正． 在划分的4 个土层中，表层受蒸发的影响较 大，临界层含水量起伏较大，变化较复杂，需要作进 一步研究，在此仅对中间两层进行修正． 粉煤灰充填土壤在中间上层有较好吻合不必 修正，对于煤矸石充填情况其修正式为 A 7 A e x p [ - I 一 r ，／r 一1 ／3 ] 。 B 7 B E l r t I r _ 1 门1 ， K 7 。一K ，e x p e l r ’／r _ 1 邝] ， 式中r ，，，．分别为充填物煤矸石的最小粒径和覆 土层土壤砂粒的分界粒径；修正后拟合系数为 0 ．8 9 50 ． 煤矸石充填土壤在中间下层有较好吻合不必 修正，对于粉煤灰充填情况其修正式为 A 7 一A e x p 1 一P7 ／p ， B ’ B 1 P ’／痧 ， K 7 。一K 。e x p 1 P ’／p ， 式中P 7 ，P 分别为其它充填物粉煤灰和充填土壤 的小于0 ．0 0 1m m 的黏粒含量；修正后拟合系数为 0 ．8 9 35 ． 4 结论 1 一次灌水后，在研究期间，覆土层水分渗漏 过程可以分成快速渗漏、平稳渗漏和起伏渗漏3 个 阶段，分成表层、中间上层、中间下层和临界层4 个 层次．由于煤矸石粒度大，煤矸石充填土壤的快速 渗漏阶段持续时间更长，临界层水分起伏最大．但 是煤矸石充填土壤的整层贮水量并不低．粉煤灰充 填土壤水分变化比较平稳．煤矸石和粉煤灰充填不 是决定整层覆土土壤含水量大小和持水时间长短 的主要因素，主要的影响因素是覆土层土壤本身的 性质．其粒度与土壤颗粒的大小的相对关系只影响 土壤水分渗漏量的变化幅度，对土壤贮水量的大小 不产生影响． 2 充填复垦中充填物料对覆土层水分的影响 主要源于充填物颗粒级配与土壤颗粒级配的对比 关系，可以用充填层与覆土层粘粒的相对含量或者 最小粒径的相对大小关系修正土壤传递函数． 3 柳新煤矸石和粉煤灰充填复垦区农田水分 条件差，不是水分渗漏造成的．可能的原因是由于 煤矸石和粉煤灰孔性的问题而对水分的吸持力不 够，阻碍了水分的上升运动，加上田间水分蒸发蒸 腾作用强，覆土层土壤水分丧失快，从而造成缺水 问题．所以在田间水分管理上应该根据需要少量多 次灌溉． 参考文献 1 - 1 3 雷志栋，杨诗秀，谢森传．土壤水动力学[ M ] ．北京清 华大学出版社，1 9 8 8 ． [ 2 3 李保国，龚元石，左强，等．农田土壤水的动态模型 及应用[ M ] ．北京科学出版社，2 0 0 0 ． [ 3 ] 张耀峰，张德生，武新乾．一维非饱和土壤水分运动的 数值模拟[ J 3 ．纺织高校基础科学学报，2 0 0 4 ，1 7 2 1 2 3 1 2 7 ． Z H A N GY a o - f e n g ，Z H A N GD e - s h e n g ，W UX i n - q i a n ．N u m e r i c a ls i m u l a t i o nf o r1 - Dw a t e rf l o wi nu n s a t u r a t e ds o i l s [ J ] ．B a s i cS c i e n c e sJ o u r n a lo fT e x t i l e U n i v e r s i t i e s ，2 0 0 4 ，1 7 2 1 2 3 1 2 7 ． [ 4 3 黄冠华，沈荣开．层状土壤中一维非饱和土壤水分运 万方数据 第5 潮熏辉等；兖壤笈爨囊壤零分鳘壹运动摸孩磷突6 9 5 麓戆隧援模掇蹬二．孛嚣缎她大学攀摄，1 9 9 7 ，2 4 3 7 4 3 ． H U A N GG u a r r h u a ，S H E NR o n g - k a i ．S t o c h a s t i c 扩 n a l y s i so to n e - d i m e n s i o n a lu n s a t u r a t e df l o wi na s t r a t i f i e ds o i l [ J ] ．J o u r n a lo fC h i n aA g r i c u l t u r a lU n i - v e r s i t y 。1 9 9 7 ，2 4 3 7 - 4 3 ． [ s ] 本拉撼龌豢因，蠹疆敬舄尔拜。点壤瘩分运动数 棼模型懿建立爱应蕉【j ] ．凝蕹农照夫学攀攘，2 0 0 2 ， 2 5 1 6 0 - 6 2 。 M U R A T l } { u s a i y i n ，H U D A N T u m a r b y ．F I o t l n - d a t i o na n da p p l i c a t i o no fn u m e r i c a ls i m u l a t i o no ft h e s o i lw a t e rm o v e m e n t [ J ] ．J o u r n a lo fX i n j i a n gA g r i c u l - r u r a lU n i v e r s i t y ，2 0 0 2 ，2 5 1 6 0 6 2 ， [ 鑫】骤溪羧，郄麓安。接求嚣稳鞠壤瘩努遂蔼参数瓣藤 接方法c J ] ．威用基础与工襁科学学撤，2 0 0 2 ，1 0 2 1 0 3 - 1 0 9 ． e 辩嚣NH o n g - s o n g ，S } { A 0M i n g - a n ．Ar e v i e wO ni 淞 d i r e c tm e t h o d sf o re s t i m a t i n gu n s a t u r a t e ds o i lb y - d r a u l i cp r o p e r t i e s [ J ] ．J o u r n a lo fB a s i cS c i e n c ea n d E n g i n e e r i n g 。2 0 0 2 ，1 0 2 1 1 0 3 1 0 9 。 [ 7 ]舀殿青，郄嘴安。菲饱和土壤水力参数的楱凝及礴定 方法口] 。应用生态学报，2 0 0 4 ，1 5 1 ；1 6 3 1 6 6 ． L UD i a n - q i n g ，S H A OM i n g - a n ．M o d e l sa n dd e t e r m i - n i n gm e t h o d so fu n s a t u r a t e ds o i lh y d r a u l i cp a r a m e t e r s [ J ] ．C h i n e s eJ o u r n a Io fA p p l i e dE c o l o g y ，2 0 0 4 。1 5 1 1 6 3 1 6 6 ． [ 观耱炎，零糅鏊。蘸嚣壤零努运臻参灵敏滋P T F s 法[ c ] ／／莓露鬻零零农盈应髑麓硝醑究遴矮．j ￡窳 科学出版社，1 9 9 5 5 6 - 6 3 ． [ 蛸黄元镑，攀韵璩，隧镰支．墨阕条传下土壤氦豢运移的 攘掇穰溅至曩3 。承鬟学摄，1 9 9 6 6 9 - 1 4 。 H U A N G Y u a n - f a n g ．L lY u n - z h u ．1 ，UJ i n - w e n ．S i m - u l a t i o no fs o i ln i t r o g e nt r a n s p o r tu n d e rf i e l dc o n d i t i o n s 王．M o d e l [ J ] 。J o u r n a lo fH y d r a u l i cE n g i n e e r - i n g ，1 9 9 6 6 g 一1 4 ． [ 1 蛰】黄元像，李魏凉瀵铸交。霹麓条季擎下i 壤藏素运移 的模拟模型H 瞰间检验岛应用[ j ] ．瘩刺学撤， 1 9 9 6 6 ；1 5 - 2 3 ． H U A N GY u a n - f a n g ，L IY u n - z h u ，1 ．1 2J i n - w e n 。 S i m u l a t i o no fs o i ln i t r o g e nt r a n s p o r tu n d e rf i e l dc o n - d i t i o n sI I 。V e r i f i c a t i o na n da p p l i c a t i o n [ J ] 。J o u r n a l o fH y d r a u l i eE n g i n e e r i n g ，1 9 9 6 6 1 5 2 3 ． [ 1 1 ]吕殿青。邵明赣，芏企九．土壤持水特征测定中的容 重变诧芨其确定方法器1 瘩蒯攀羧，2 0 0 3 3 ；1 1 0 一 l 至。 L UD i a n - q i n g ，S H A OM i n g - a n tW A N G Q u a n - j i u ． B u l kd e n s i t yc h a n g i n gd u r i n gm e a s u r i n gs o i lw a t e r r e t e n t i o nc h a r a c t e r i s t i c sa n di t sd e t e r m i n i n gm e t h o d 口] ．J o u r n a lo fH y d r a u l i cE n g i n e e r i n g ，2 0 0 3 3 ； 1 1 0 - t 1 9 ． [ 1 2 3G A R D N E RwR ．S o m es t e a d ys t a t es o l u t i o no ft h e u n s a t u r a t e dm o i s t u r ef l o we q u a t i o nw i t ha p p l i c a t i o n t Oe v a p o r a t i o nf r o maw a t e rt a b l e [ J ] 。S o i lS c i e n c e ， 1 9 8 5 8 5 2 2 8 2 3 2 。 爨谨壤辏臻嚣 万方数据