矿区开采沉陷农用土地质量空间变化研究.pdf

第3 3 卷第2 期 中国矿业大学学报 V 0 1 ．3 3N o ．2 2 0 0 4 年3 月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 y M a r ．2 0 0 4 文章螭号1 0 0 0 1 9 6 4 2 0 0 4 0 20 2 1 3 0 6 矿区开采沉陷农用土地质量空间变化研究 卞正富 中国矿业大学环境与测绘学院，江苏徐州2 2 1 0 0 8 摘要提出利用土地环境务件、土壤物理特性、土壤养分因子反映矿区土地质量．通过徐州矿区 的试验研究，分析了开采沉陷引起的土地环境条件、土壤物理特性、养分因子的变化规律，得出 开采沉陷后土地产生附加坡度，排灌标准明显下降，农田林网防护范围缩小，地下水位埋深相对 变浅I 土壤物理特性中下坡位置土壤含水率明显升高，土壤客重几乎没有变化，总孔隙度在下坡 处明显变小；而土壤养分因子上坡各项指标普遍低于对照田，下坡位置要高于对照田，上、下坡的 养分因子平均值与开采沉陷前没有显著差异．利用熵流指数模型分析知开采沉陷后排灌条件 差、防护林防护范围小、地下水位埋深不理想是导致开采沉陷区土地生产力下降的主要原因． 关键词开采沉陷，土地利用} 土壤环境；土地复垦 中圈分类号；T D3 2 7文献标识码A C h a n g eo fA g r i c u l t u r a lL a n dQ u a l i t yd u et oM i n i n gS u b s i d e n c e B I A NZ h e n g f u S c h o o lo fE n v i r o n m e n t a la n dS p a t m lI n f o r m a t i c sE n g i n e e r i n g ．C U M T ，X u z h o u ，J i a n g s u2 2 1 0 0 8 C h i n a A l №t r a e tlI ti 8p r o p o s e dt h a tm i n e dl a n dq u a h t yc a nb ee s t i m a t e db ye n v i r o n m e n t a lc o n d i t i o n s ，s o i l c o n d f f i o n s ．s u c ha sn u t r i e n tc o n t e n t s 。p h y s i c a lf a c t o r s ．‘T h el a w so ft h ec h a n g eo fe n v i r o n m e n t a l c o n d i t i o n sa n ds o i ic o n d i t i o n si ns u b s i d e da r e aw e r ep r o p o s e db ym e a n so ff i e l ds u r v e y ．T h er e s u l t s s h o wt h a t1 a d d i t m n a ls l o p ei sf o r m e d ，t h el e v e lo fd r m n a g ea n di r r i g a t i o nd e c r e a s e ss t r i k i n g l y ，t h e d e f e n s ea r e ab yn e to ff a r m l a n df o r e s t r yr e d u c e s ．a n dt h el e v e lo fp h r e a t i cw a t e rb e c o m e ss h a l l o w a f t e rs u b s i d e d l2 s o i l w a t e rc o n t e n t i n l o Ws l o p e i n c r e a s e s ，b u l kd e n s i t y n e a r l y k e e p sc o n s t a n t ，a n d t o t a lp o r o s i t yd i m i n i s h e sa td o w ns l o p e ；3 a 1 1 i n d i r c so fs o i ln u t r i e n ti nt h ec o n t r a s tf a r m l a n da r e l o w e rt h a nt h o s ei nd o w ns l o p ea n dh i g h e rt h a nt h o s ei nu ps l o p e ，b u tt h ea v e r a g eo fs o i ln u t r i e n t i n d i c e si nd o w ns l o p ea n du ps l o p ei sn e a r l ye q u a lt oc o n t r a s tf a r m l a n d ．T h ep o o rc o n d k i o n so f d r a i n a g ea n di r r i g a t i o n ．t h es m a l lp o r t i o no fd e f e n s ea r e ab yn e to ff a r m l a n df o r e s t r y ，a n dt h e s h a d o wl e v e lo fp h r e a t i rw a t e ri ns u b s i d e da r e aa r et h ep r i m a r yf a c t o r st oc a u s el a n dp r o d u c t i v i t y r e d u c t i o nb ym ∞t n so fe n t r o p ye x c h a n g em o d e l ． K e yw o r d s m i m n gs u b s i d e n c e ll a n du t i l i z a t i o n } s o f te n v i r o n m e n t } l a n dr e c l a m a t i o n 地下开采产生下沉盆地．徐州等矿区由于潜水 位埋藏浅，下沉盆地通常积水，这是煤炭地下开采 导致土地减产、甚至绝收的主要现象．但是在不积 水的坡地也存在以下现象1 土地减产减收边界 大于塌陷边界，如图1 a 所示；2 土地减产减收边 界小于塌陷边界，如图1 b 所示} 3 在塌陷边界附 近出现增产增收现象，如图1 c 所示． 图1 中的塌陷边界是指下沉值等于1 0m m 的 等值线围成的边界，土地减产减收边界通常是矿山 企业征地或赔偿边界． 出现上述现象的原因是开采沉陷导致土地质 量因子变化，从而导致土地生产力变化．但土地质 收藕日翔l2 0 0 3 1 0 一1 0 薹盒珥目t 全国优秀博士论文基金项目 2 0 0 0 4 4 作着简介。卞正直 1 9 6 5 ， ．男．江苏省建瑚县人，中国矿业太学教授t 工学博士，从事土地赍源、土地复星与生态重建方面的研氪 万方数据 中国矿业大学学报第3 3 卷 ＼竺 ／一 飞＼竺 ／胪 积h 竺 彳盛 图1 土地产量增减边界与塌陷边界的关系 F i g ．1T h eb o u n d a r i e so fs u b s i d e n c ea n dd e c r e a s e o ri n c r e a s e o fc r o pp r o d u c t i o n ．塌陷边界I ◆．减产减收边界 量因子变化规律如何 沉陷区土地生产力变化规律 又如何 文献E l i 曾将出现图1 a 现象的原因假设为 土地沉陷后，地下水位的变化、土壤肥力的下降、耕 作管理措施跟不上等，因而在塌陷边界以外土地生 产力也下降；而出现图1 c 现象的原因可假设为土 地沉陷对土壤环境有一定的改善作用，即开采沉陷 在边界附近为正影响，因而使塌陷边界附近土地生 产力有所提高．RTH e t z l e r 和RGD a r m o d y 曾对 消除开采沉陷影响的措施对土壤及作物产量的影 响进行了研究口] ，国内外还有许多关于研究开采沉 陷对潜水位影响、对地表裂隙的影响、对土壤物理 特性、化学特性影响的研究报道[ 3 - 7 ] ，但均没有建立 起这些影响与沉陷区土地生产力的关系．本文通过 试验与理论分析探讨了开采沉陷对土地质量因子 及其对土地生产力的影响规律． 1原理、研究场地与试验方法 1 ．1 采用的原理 文献[ 8 ] 提出影响土地生产力大小的因素分为 三类；投入、环境与土壤．三要素对土地生产力的影 响机制如图2 所示，即；投入是必要的输入，土壤是 营养库和投入要素转化为生产力的场所．环境则起 激励或抑制作用，作用于土壤要素，三要素共同作 用的结果是土地具备生产功能，三要素缺一不可， 且存在一最佳匹配状态使得土地获得最大生产力． 这里土地的投人要素是指肥料等生产资料、人工及 臂稚库 图2 影响土地生产力的三要素及其关系 F i g ．2T h r e ee l e m e n t sa n dt h e i rr e l a t i o n s h i p s a f f e c t i n gh n dp r o d u c t i v i t y 管理的投入等，环境要素则是指土地的地形条件、 排灌条件、道路与防护林网条件、地下水位埋深等， 土壤要素是指土壤养分、土壤物理特性等． 在投入条件相同的情况下，环境与土壤要素就 成为决定土地生产力大小的两个重要因素． 1 ．2 试验地点与场地条件 徐州矿区姹城矿魏庄沉陷区 简称S ， 和大黄 山矿西翼沉陷区 简称s 。 均属于高潜水位类型， 试验地点均未积水，且能耕种，但开采沉陷后有一 定的坡度，S 、坡度在1 ／1 0 0 左右，S 坡度在1 ／2 0 0 左右．s ，，s 两场地均为稻麦两熟地，取样日期为6 月初麦收后．s 。场地小麦产量42 0 0k g ／h m 2 左右， S z 场地小麦产量15 0 0k g ／h m 2 左右．S 。对照田产 量约48 0 0k g ／h m 2 ，S 2 对照田产量在60 0 0 k g ／h m 2 以上． 1 ．3 试验方法 试验过程中，重点对环境条件进行了野外调查 与测定，调查与测定的内容包括地形坡度、排灌设 施及农田道路林网的完好度及地下水位埋深．土壤 条件测定了物理特性和养分状况，具体指标包括 土壤含水量、容重、速效氮、速效磷、速效钾、有机质 含量． 土地环境条件的调查主要通过1 1 00 0 0 地 形图和土地利用现状图获取信息，潜水位调查通过 野外挖坑，2 d 后观察坑中水位的方法．土壤含水量 测定采用小铝盒法，容重测定采用环刀法，比重测 定采用比重瓶法，养分测定利用常规方法测定．取 样地点分别为上坡和下坡及对照田，如图3 所示． h N 束沉陷对j f l f 田 S 卜器一 图3 场地取样点位置示意 F i g ．3 L o c a t i o no fs a m p l i n g 土壤含水率Ⅳ，土壤容重d ，土壤密度p ，土壤 总孔隙度P ，分别用式 1 ～ 4 计算 W [ 铝盒 湿土重 一 铝盒 干土重 3 ／ [ 铝盒 干土重 一铝盒重] 1 0 0 ％ ， 1 d 一 M x l 0 0 ／[ 1 0 0 W V ] g ／e r a 3 ， 2 P 一 Ⅳ2 一W 1 I [ W 4 W 2 一W 1 Ⅳ3 ] t C 时水 的密度， 3 尝一 卜 万方数据 第2 期卞正富矿区开采沉陷农用土地质量空间变化研究 2 1 5 P Pd ／P 1 0 0 ％ ， 4 式中M 为环刀中湿土重，g ；V 为环刀体积，c m 3 ， 试验用环刀体积V 一1 0 0c m 3 ；Ⅳ。为比重瓶重，g ； Ⅳ为比重瓶内放人烘干土 过1m m 土壤筛 后的 重量，g ；Ⅳ。为称重Ⅳ后蒸溜水在砂浴上煮沸去 除土和水中的空气，再加满蒸馏水，称取的重量，g ； Ⅳ。为比重瓶中水与土倒出洗净，用同样室温的蒸 溜水装满，称取的重量，g ． 2 结果分析 2 ．1 开采沉陷前后土地环境条件对比 两个试验场地开采沉陷前后土地环境条件均 有较大的变化，变化的共同点是土地产生附加坡 度，生产道路、沟渠、农田林网均遭受破坏，地下水 位埋深相对变浅．两个场地的变化对比如表1 ．S 。 场地开采沉陷尚未稳定，在其南侧的塌陷地基本稳 定，且已经过复垦改造，改造时开挖了魏庄大沟，魏 庄大沟不仅可作为两侧农田的排水通道，还可兼作 蓄水灌溉水源．S 场地开采沉陷已经稳定，塌陷深 度相对较大，沉陷后土地未采取任何治理措施，但 尚能利用，只是产量水平较低． 表1 试验研究场地土地环境条件对比 T a b l e1 C o m p a r j 8 0 no fl a n dc o n d i t i o n s a te x p e r i m e n t a ls i t e 试验场地坡度坡向絮掌誊勰鼗夸髅牌 2 ．2 开采沉陷前后土壤条件对比 开采沉陷前后由于开采沉陷的影响、地形的变 化，利用措施的差异等，土壤条件也发生了微小的 变化，具体指标如表2 ， 表2 试验研究场地壤条件对比 T a b l e2 C o m p a r i s o no fs o i lf a c t o r sa te x p e r i m e n t a ls i t e s 由表2 看出开采沉陷前后土壤物理特性中土 壤含水率变化较大，尤其是下坡位置土壤含水率明 显升高，土壤容重几乎没有变化，只有S ，场地下坡 土壤容重变大，可能与魏庄大沟施工时运土机械碾 压或生土混入有关，总孔隙度只在下坡处明显变 小，这与地表径流或灌溉退水将上坡土壤细微颗粒 带到下坡有关．土壤养分因子开采沉陷前后的变化 规律较为明显，即上坡各项指标普遍低于对照田、 而下坡位置要高于对照田，这一现象说明开采沉陷 形成的地面坡度变化对土壤养分的空间分布存在 显著影响． 2 ．3 土地生产力与土地环境条件、土壤条件的相 关关系 作物产量是反映土地生产力水平的综合指标， 即反映了土地环境条件、土壤条件、投入水平 生产 资料的投入、土地经营管理水平等 三方面的综合 影响．从实际效果看s ，场地的产量水平较对照田 差别不是很太，而S 场地的产量水平与对照田差 别则较大．对照表1 ，表2 两个试验研究场地开采 沉蹈前后土地环境条件与土壤条件的变化与产量 水平变化可以发现试验场地土地生产力水平的下 降是土地生产的环境条件和土壤条件变化共同作 用的后果，尽管S 。场地开采沉陷后土壤条件除了 速效磷与土壤含水率较差外，其他指标均好于s 。 场地，且两个场地下沉值均在2m 左右，但是其生 产力水平较s 。场地下降明显得多，其原因主要是 土地利用的环境条件的变化对土地生产力影响较 大所致．S 。场地的潜水位埋深只有1 ．1m ，更为重 要的是排灌设施已经失修多年，从而导致土壤含水 率偏高，不利于小麦生长．由此我们还可以看出土 地利用的环境条件与土壤条件对土地生产力同等 重要，S ，场地南侧经过土地复垦、修建的魏庄太 沟对改善S 。场地的土地利用环境具有明显的效 果． 2 ．4 土地生产力与其影响因子之间关系的定性分 析 如前所述，影响土地生产力的因素包括土地利 用环境因素、土壤因素与投入因素，一方面三类因 万方数据 中国矿业大学学报第3 3 卷 素中的单个因子应该符合作物生长条件，另一方 面，三类因素存在一相互匹配的关系，才能使得土 地具有较高的生产力水平．对开采沉陷土地而言， 通常土地的环境条件与土壤物理特性变化较大，只 要变化后仍然符合作物生长的有利条件，通常对土 地生产力影响就不大，相反则影响显著．因此土地 生产力水平可用下面的函数表示 P 西L ／』 z ，J ] ， 5 式中毒为土地生产力函数} ，，为环境因子、土壤因 子、投入因子在o ，∥ 处的取值． 对于开采沉陷前后土地生产力的变化可用图 4 所示的概念模型表示．图4 中，，m o ，y 为 z ，y 点处J 因子的初始值； h ，y 为 z ，_ 点处j 因子 沉陷后的值；凡为o ，y 点影响土地生产力的初始 因子值组成的集合；F 为 z ，， 点沉陷后与F 。相对 应的因子值组成的集合；P 。为土地初始状态的生 产力} P 则为沉陷后土地的生产力． 秘∥ 班≯ 一一F 巡厂 图4 土地沉陷对土地生产力影响的概念模型 F i g ．4C o n c e p t u a lm o d e lf o rt h ee f f e c to fl a n d s u b s i d e n c e0 1 3 ．1 a n dp r o d u c t i v i t y 只要F 集合整体好于或接近F 0 ，沉陷后土地 生产力变化就可能提高或下降不大，相反土地生产 力将明显下降．这就是对出现图1 所示现象的定性 解释． 2 ．5 土地生产力与其影响因子之间关系的定量分 析 反映土地生产力与其影响因子之间定量关系 的方法较多，如模糊数学方法、逐级修正模型以及 卡尔曼滤波方法等o ] ，但能够较好反映土地生产过 程中环境因子、投入因子与土壤因子对土地生产力 影响的模型有熵流模型o ] ．下面利用熵流模型对两 个试验场地的土地生产力进行评估．熵流模型用下 面的指数表示 眦 一竞撕踞 一扩世掣时 - 1 1 皿高剑时帮‰ 6 式中 一1 8 为符号函数，E t ，S f 好于初始状 态，n 取为l ，否则取为2 ；w ，～∞。为不同因子的影 响程度；E o ，s o ，， o 分别为环境、土壤、投入 因子的初始状态，可取对照农田或沉陷前土地的状 态值} E f ，s } ，， f 则为不同状态 如沉陷后、复 垦后 环境、土壤、投人因子的取值． 不难得出系统初始熵丘 o 一0 ，且熵流各分 量的符号符合熵流的物理意义，即能量流人为正， 流出为负，物质流出为正，物质流人为负．如果 E A t 为负值，说明土地生产力水平高于初始状态 或对照状态 ，反之较初始状态差． 这里假设两个场地的投入水平相同，可只比较 土壤因子与环境因子的影响．于是熵流指数为 E e o 一妻w , S E , 一 一妒坦掣蚶 一1 ⋯7 一”拯等i 掣蚧 7 利用式 7 和表3 数据，将小麦理想生长条件 看成是初始状态，首先将表3 中的1 1 个因子作归 一化处理，坡度、水位埋深、土壤容重、速效氮、磷、 钾及有机质含量归一化处理分别利用式 8 ～ 1 4 的适应性函数计算[ 1 0 3 ．坡向朝南为1 ，朝北为0 ；排 涝标准大于等于1 0a 一遇为1 ，5a 一遇为0 ．5 ，小 于5a 一遇为o ～o ．5 ；灌溉保证率为1 0 0o A 时取值 为l ，5 0 ％时取值为0 ．5 ；农田林网防护范围是指在 树高1 5m 的情况下，防护范围可达2 0 0m ，取值为 1 ，如果树高为5m ，则防护范围达到6 7m ，取值为 0 ．3 3 ．由于土壤含水率、总孔隙度与地下水位埋深、 土壤容重存在一定的相关关系．所以这里未考虑这 两个指标． 。J 1 s ≤1 ／4 0 0 ， ⋯ 。“一1 1 ／ 4 0 0 S S 1 ／4 0 0 ． ’ f 是H I ．5 ， s 。№一 2 ， 龟8 。8 三- 0 。6 。8 罴d 鲨 i { ．f 1 a N ≥1 0 0 ， 5 “一I o ．0 1 口N 州∈[ o ，1 0 0 3 ， f 1a P ≥1 5 ， P 一{ 0 ．0 5 a P O ．2 5a P E [ 5 ，1 5 ， 【0 ．1 a P a P E [ o ，5 ， f 1a K ≥2 0 0 ， s 堪一 3 1 0 - 3 4 K 0 ．4a K 6 [ 1 0 0 ，2 0 0 ， 1 0 ．0 7 a K a K E [ o 。1 0 0 ， f 1O M ≥2 ， ．J o ．5 8 8 2 0 M O M E o ．3 ，2 ] ， 6 “一1 o ．1 7 6 5 ”7 6 2 L 0 O M ≤0 。3 ， 9 1 0 式中．s 蜘，s 帆№，＆，S ．N ， P ，S , x ，s 洲分别为坡度、 ㈣ m Ⅲ Ⅲ 万方数据 第2 期卞正富；矿区开采沉陷农用土地质量空间变化研究2 1 7 水位埋深、土壤容重、速效氮、磷、钾及有机质含量表示坡度、水位埋深、土壤容重、速效氮、磷、钾及有 归一化后的取值，s ，H ，d ，a N ，a P ，a K ，O M 则分别机质含量的实际取值． 裹3 计算螭流指数的初始因子值 T a b l e3I n i t i a lv a l u eo ff a c t o r sf o rc a l c u l a t i n gi n d e xo fe n t r o p yf l o w 将土地环境条件与土壤条件对土地生产力的 影响程度看成是相同的，并对环境条件、土壤条件 中具体指标的影响程度分解如表4 中的w 值，计 算两个场地开采沉陷前后的熵流指数结果见表4 ． 裹4 计算熵流指数的因子归一化后的取值 T a b l e4T h ev a 岫o fn o r m a l i z e df a c t o r s 对比熵流指数与小麦的实际产量水平发现两 个场地开采沉陷前后的土地利用环境条件与土壤 条件均未达到理想状态，所以熵流指数均为正值， 且开采沉陷后均较开采沉陷前太，说明开采沉陷后 土地生产力水平均下降，但是s 。场地开采沉陷后 的熵流指数与开采沉陷前的熵流指数相差不大，如 果开采沉陷后的坡度因子和农田防护林网防护范 围因子达到开采沉陷前的水平，熵流指数将与开采 沉陷前持平并略好于开采沉陷前．s ，场地之所以 开采沉陷后的产量下降不大，与其灌排设施完善和 地下水位埋深适宜是密切相关的．对照s 。场地开 采沉陷前后的变化，发现排灌条件、防护林防护范 围、地下水位埋深均较开采沉陷前差，导致土地生 产力下降明显． 由本试验研究还可推断如果开采沉陷后 通 过土地整治等措施，如s 。场地开挖魏庄大沟 使得 影响土地生产力的各项指标较开采沉陷前变得更 加理想，土地生产力水平将高于开采沉陷前；本试 验研究中尚未考虑投入因子的影响，参照式 6 发 现如果通过适当的投入，降低熵流指数值，实现沉 陷后土地生产力水平高于开采沉陷前也是完全可 能的． 万方数据 中国矿业大学学报第3 3 卷 3 结论 1 徐州矿区开采沉陷前后土地环境条件均有 较大的变化，土地环境条件的变化对土地生产力水 平的影响是显著的．这些变化包括土地产生附加坡 度，排灌标准明显下降，农田林网防护范围缩小，地 下水位埋深相对变浅． 2 通过两个试验场地的调查与取样分析发 现，徐州矿区开采沉陷前后由于地形的变化、土地 利用措施与经营管理的差异等，土壤物理特性中土 壤含水率变化较大，尤其是下坡位置土壤含水率明 显升高，土壤容重几乎没有变化，总孔隙度在下坡 处明显变小．土壤养分因子开采沉陷前后的变化规 律较为明显，即上坡各项指标普遍低于对照田、而 下坡位置要高于对照田，而上下坡的养分因子平均 值与开采沉陷前没有显著差异． 3 作物产量是反映土地生产力水平的综合指 标，即反映了土地环境条件、土壤条件、投入水平 生产资料的投入、土地经营管理水平等 三方面的 综合影响．从试验研究场地看S ，场地的产量水平 较对照田差别不是很大，而S 场地的产量水平与 对照田差别则较大，试验场地土地生产力水平的下 降是土地生产的环境条件和土壤条件变化共同作 用的后果． 4 土地生产力及其影响因素之间的关系可用 熵流指数模型表示．该模型较好地反映了开采沉陷 前后土地生产力的差异．利用熵流指数模型分析 知S 。场地开采沉陷后排灌条件、防护林防护范 围、地下水位埋深均较开采沉陷前差，是导致其土 地生产力下降的主要原因．同时通过本试验研究可 以推断如果通过适当的投入，降低熵流指数值，实 现沉陷后土地生产力水平高于开采沉陷前也是完 全可能的． 参考文献 [ 1 3 卞正富．矿区土地复垦界面要素的演替规律及其调 控研究E D ] ．徐州中国矿业大学采矿系，1 9 9 7 ． [ 2 ] H e t z l e rRT ，D a r m o d yRG ．C o a lm i n es u b s i d e n c e m i t i g a t i o n E f f e c t so n s o i la n dC r o pY i e l d s [ A ] ． R o b e r tED ，B a r n h i s e lRI ，D a r m o d yRG ． P r o c e e d i n go f 1 9 9 2 N a t i o n a ls y m p o s i u mo nP r i m e F a r m l a n dR e c l a m a t i o n [ c ] ．U r b a n a U n i v e r s i t yo f 1 U | n o i sP r e s s ，1 9 9 2 ．1 2 9 1 3 6 ． E 3 ] S e l l sDE ，D a r m o d yRG ，S i m m o n sFW ．T h ee f f e c t s o fc o a lm i n es u b s i d e n c eo ns o i lm a c r o p o r s i t ya n dw a t e r f l o w [ A ] ．R o b e r tED 。B a r n h i s e lRI ．D a r m o d yRG ． P r o c e e d i n go f 1 9 9 2 N a t i o n a lS y m p o s i u mo nP r i m e F a r m l a n dR e c l a m a t i o n [ c ] ．U r b a n a ；U n i v e r s i t yo f l m n o i sP r e s s 。1 9 9 2 ．1 3 7 1 4 6 6 ． [ 4 ] R o o s e n d a a lV a nDJ ．L o n g w a l lm i n es u b s i d e n c eo f f a r m l a n di nS o u t b e r nI l l i n o i s N e a rs u r f a c ef r a c t u r i n g a n da s s o c i a t e dh y d r o g e o l o g i c a le f f e c t s [ A ] ．R o b e r tE D ，B a r n h i s e lRI ，D a r m o d yRG ．P r o c e d h go f1 9 9 2 N a t i o n a ls y m p o s i u mo nP r i m eF a r m l a n dR e o l a m a r i o n [ c ] ．U r b a n a U n i v e r s i t yo fI l l i n o i sP r e s s ，1 9 9 2 ．1 4 7 一 l5 8 ． E S ] 陈霓乾，邓喀中，赵志海，等．开采沉陷对耕地土壤物 理特性影响的空间变化规律口] ．煤炭学报，1 9 9 9 ，2 4 6 5 8 6 5 9 0 ． [ 6 ] 卞正富．开采沉陷对潜水环境的影响及其控制口] ． 有色金属，1 9 9 9 ，5 1 1 ；4 - 7 ． [ 7 ] D a t t aKK ，J o n gC ．A d v e r s ee f f e c to fw a t e rl o g g i n g a n ds o i le , a b n i t yo nc r o pa n dl a n dp r o d u c t i v i t yi n n o r t h w e s tr e g i o no fH a r y a n a ，I n d i a [ J ] ．A g r i c u l t u r a l W a t e r M a n a g e m e n t ，2 0 0 2 5 7 2 2 3 2 3 8 ． [ 8 ] 卞正富张国良．矿山复垦土地生产力的耗散结构模 型口] ．煤炭学报，1 9 9 8 ，2 3 6 6 6 3 6 6 8 ． [ 9 ] C h e nM ．D e v e l o p m e n ta n dt h eE n v i r o n m e n t [ D ] ． B e r k l e y U n l v e r s i t yo fC a l i f o r n i a ，2 0 0 1 ． [ 1 0 ] 卞正富，张国良．矿山复垦土壤生产力指数的修正 模型口] ．土壤学报，2 0 0 0 ，3 7 1 1 2 4 1 3 0 ． 责任编辑陈其泰 万方数据