2024年3月9日发(作者：西门英悟)
流量和含沙量对水质参数影响分析
石 伟，王光谦
(清华大学 水利水电工程系)
摘要：为了研究近年来黄河下游水环境容量萎缩的内在机理，本文用流量级频率分析的方法，对1982～1984年黄河下游花园口、泺口水文站实测水质分析成果表中的耗氧量、溶解氧及其相应流量进行了分析，得到在平滩流量下黄河下游河道对污染物输运能力最强，而且此时其水环境容量也最大的结论。随着花园口、泺口水文站平滩流量由1958年的8000m/s和9200m/s降至1999年的4000m/s和3000m/s，黄河下游水环境容量萎缩是必然的。本文又通过统计分析得知含沙量对黄河下游水质监测数据中的耗氧量及溶解氧量影响不大。
关键词：水环境容量；平滩流量；黄河下游
基金项目：国家重点基础研究发展规划(973)项目(G19990436)；清华大学“985”校级重点科研项目软课题
作者简介：石伟(1967-)，女，内蒙古呼和浩特人，博士研究生，研究方向为水力学及河流动力学。
3333近年来，在人类大量引水、降水偏枯[1，2]等人为和自然因素作用下，黄河下游的来流明显减少，不仅表现出断流、洪水频繁、河道萎缩、主槽淤积、“小水大灾”[3，4]等水沙过程的变异，而且呈现出水环境容量萎缩、水体自净能力下降、水质恶化[3]，1996年山东黄河有1/4的时间是超标水，1997年上升为1/2，1998[5]年再度上升为近3/4，污染形势极为严峻。河流是污染物迁移载体，污染物是借助于河流中的水流和泥沙来输运的。黄河来流流量的大小影响着河流携带污染物的能力，影响着河流的水环境容量及其自净能力，夏军等人已证实增加来水量可以改善河段水环境容量[6]。进入黄河的泥沙，给黄河造成一定污染，同时由于其组成含有相当数量的粘土矿物、无机胶体和一定数量的有机胶体、有机无机复合胶体，对排入黄河的种类繁多的污染物具有显著的吸附作用，从而表现出净化水质的效应，高传德等人通过水质监测数据得到浑水高锰酸钾指数与其含沙量呈[7]正相关。本文首次采用流量级频率分析的方法研究流量对黄河下游耗氧量和溶解氧的影响，采用统计分析的方法分析含沙量对其耗氧量和溶解氧的影响，是对黄河下游水沙变化对其水质影响机理的初步探索。
1 水质实测数据的选取
耗氧量(COD)是指在一定氧化条件下，以一定的氧化剂氧化水或泥沙中还原性物质所消耗的氧化剂量。它可以作为水体有机物相对含量的指标，常用作反映水质污染程度的参数[7]。溶解氧是水中溶解态氧的浓度，是衡量水质好坏的一项重要指标。维持良好水质感官性状、水体的自净作用以及水生生物的生存需要有[8，9]氧。本文选取了1982～1984年花园口、泺口水文站实测水质分析成果表中的耗氧量、溶解氧及其相应的流量进行分析，含沙量在有相应时段洪水水文要素摘
录表时用相应时间、相应流量的含沙量，若无相应时段洪水文要素摘录表则用相应的日平均含沙量。
2 耗氧量和溶解氧量流量级频率分析
2.1 耗氧量和溶解氧量流量级频率分析方法 Wolman和Miller在1960年给出了河流有效流量的各级流量——频率分析的计算方法[10～13]：假设河流流量频率服从对数正态分布，输沙率作为流量功函数，则输沙效率即流量出现频率与输沙率乘积最大处的流量值为输沙的有效流量，就长时期内河流输沙量而言在该流量处最大[11～13]。本文借用此方法，假设河流流量频率服从对数正态分布，耗氧效率即流量出现频率与流量和耗氧量乘积最大处的流量值是长时期内河流耗氧量最大的流量，此分析方法作为耗氧量的流量级频率分析。如果将上述耗氧量换为溶解氧量，则是溶解氧的流量级频率分析。
2.2 耗氧量和溶解氧量流量级频率分析结果 图1为花园口、泺口水文站1982～1984年各级流量耗氧效率(即流量出现频率与流量和耗氧量乘积)、各级流量溶解氧效率(即流量出现频率与流量和溶解氧量乘积)，均以1000m3/s分级。从图中可以看到，花园口、泺口水文站在1982～1984年相应于耗氧效率(即流量出现频率与流量和耗氧量乘积)、溶解氧效率(即流量出现频率与流量和溶解氧量乘积)最大处的流量分别为：6348m3/s、5640～6348m3/s，二者基本相当；它们出现的天数均占总测量天数的7%.对照图2可以看到，花园口、泺口水文站在1982～1984年相应于耗氧效率和溶解氧效率最大处的流量6348m3/s和5640～6348m3/s与累计耗氧率及其累计溶解氧率曲线由陡变缓的转折点6460m3/s在数值上相近。也就是说，相应于耗氧效率和溶解氧效率最大处的流量也是黄河下游输送污染物最经济的流量。事实上，相应于耗氧效率和溶解氧效率最大处的流量6348m3/s和5640～6348m3/s与该时段花园口、泺口水文站的平滩流量相当。根据明渠流连续方程
Q=UBH (1)
式中：U为断面平均流速(m/s)，B为河面宽度(m)，H为水深(m).根据方程(1)，在平滩流量之下，当水深H和断面平均流速U随流量增大而增大时，河流宽度B保持不变[10]，因而B/H减小。当水流漫滩以后，由于滩地糙率比主槽糙率大好几倍，因而随着水位继续上涨，主槽和全断面的流速反而有所下降，当滩地水深约为主槽平均水深的35%时，断面流速达到最低值，自此以后，又随着水位的抬高而不断加大[10]。所以，相应的河道输运污染物能力在平滩流量之下随着流量的增加而增大，直至平滩流量时达到最大，在流量超过平滩流量以后，随着主槽和全断面流速的降低，输运污染物能力反而下降。
图1 花园口、泺口水文站各级流量耗氧效率及其溶解氧效率
图2 花园口、泺口水文站累计耗氧率及其
累计溶解氧率与流量关系
由于黄河下游河道在平滩流量下不仅对污染物输运能力最大，而且此时河道中的溶解氧也最多，因此在平滩流量下黄河下游河道的水环境容量为最大。在平滩流量时不仅输沙量最大[10]，而且输送污染物量也最大，为我们估算黄河下游生态需水量提供了依据，即根据平滩流量估算的黄河下游用水量不仅是输沙最经济的水量，也是输送污染物最经济的水量。
2.3 黄河下游水环境容量萎缩的原因 从表1可以看到，从1958年到1999年花园口和泺口水文站的平滩流量经历了由大到小到大又到小的过程[4]。1958年黄河下游为基本未受人类影响的自然状态，此时平滩流量较大，花园口和泺口33分别为8000m/s和9200m/s；1964年在三门峡水库的冲刷作用下，花园口平滩流量增大，为9000m3/s；1973年在三门峡水库滞洪排沙作用下，黄河下游平滩流量急剧减小，花园口和泺口的平滩流量降为3500m3/s和3100m3/s；1985年由于丰水少沙，黄河下游河道冲刷，平滩流量上升，花园口和泺口的平滩流量又上升为6900m3/s和6000m3/s；1985年以后，由于人为、自然等因素的作用，黄河下游来水来沙条件改变，河道严重萎缩，主槽淤积严重，平滩流量下降，花园口、泺口水文站平滩流量变为4000m3/s和3000m3/s.由于在平滩流量下黄河下游河道的水环境容量为最大，所以随着平滩流量的上述变化，水环境容量也会发生相应的变化。1985年以后，由于人类大量引水、降水偏枯等原因，黄河下游断流频繁，不仅使黄河下游的河床演变过程发生改变，河道萎缩、“小水大灾”、洪水频繁，而且也大大降低了其输运污染物的能力以及溶解氧的能力，表现为水环境容量萎缩、水质明显下降。与此同时，随着经济发展、人口增加，黄河流域废污水排放量与日俱增。据统计，20世纪70年代后期年排放废污水18.5亿t，80年代初为21.7亿t，90年代初已达32.6亿t，10年间增加50%以上[6]。与此相应，黄河下游必然污染情况越来越严重。
表1 典型年份花园口、泺口水文站平滩流量变化(m3/s)
时间
1958年汛后
1964年汛后
1973年汛前
花园口
8000
9000
3500
泺 口
9200
3500
3100
1985年汛前
1999年汛前
6900
4000
6000
3000
3 含沙量对耗氧量和溶解氧量的影响分析
考虑到耗氧量和溶解氧量与含沙量是同一个量纲，我们将花园口、泺口水文站1982～1984年实测水质分析成果表中的耗氧量(单位：mg/l)和溶解氧量(单位：mg/l)分别乘以1000后再分别除以含沙量(单位：kg/m3)得到两个无量纲的量，将这两个分别与耗氧量和溶解氧量有关的无量纲的量作为纵轴，均以含沙量为横轴作图，得到图3和图4.在图3中我们点绘了双曲线在第一象限的一支，其常数为1982～1984年花园口、泺口水文站所有实测耗氧量的算数平均值(约为1.46)乘以1000，横轴和纵轴同上所述；在图4中我们同样点绘了双曲线在第一象限的一支，其常数为相应溶解氧量的算数平均值(约为8.15)乘以1000，横轴和纵轴同上所述。从图3和图4中可以看出，实测点与双曲线吻合得很好，也就是说，耗氧量和溶解氧量可以看作是个常数，约等于其实测数据的算数平均值，而与含沙量无关。究其原因，与耗氧量的测定方法不无关系，因为多年来黄河上都是用去除泥沙后的清水测定耗氧量[7]，只代表了黄河水溶液中耗氧物质的含量。事实上，黄河水体中的泥沙吸附着大量的污染物，其吸附量必然与含沙量有关，也必然随着吸附量达到饱和而使水体受到二次污染。由于缺乏浑水的资料[7]，限制了我们对耗氧量和溶解氧量与含沙量关系的更深入的认识。
图4 花园口、泺口水文站溶解氧量与含沙量之图3 花园口、泺口水文站耗氧量与含沙量之比
比
(乘以1000)与含沙量的关系
(乘以1000)与含沙量的关系
4 结论
本文通过对黄河下游花园口、泺口水文站1982～1985年实测水质监测数据的各级流量频率分析，得到黄河下游河道在平滩流量下对污染物输运能力最大，此时河道中的溶解氧也最多，也就是说，在平滩流量下黄河下游河道的水环境容量最大。这一研究丰富和发展了前人的研究，即河流在平滩流量时不仅输沙能力[10]最大，而且输运污染物能力也最大。黄河下游水环境容量萎缩是平滩流量降低
的必然结果。同时，通过统计分析得知含沙量对黄河下游水质监测数据中的耗氧量及溶解氧量影响不大，即黄河下游水体水溶液的耗氧量和溶解氧基本不变。
参 考 文 献：
[1] 许炯心。人为季节性河流的初步研究[J]。地理研究，2000，19(3)：234-242.
[2] 刘昌明，成立。黄河干流下游断流的径流序列分析[J]。地理学报，2000，55(3)：257-265.
[3] 常炳炎，薛松贵，张会言，等。黄河流域水资源合理分配和优化调度[M]。郑州：黄河水利出版社，1998.
[4] 申冠卿，张晓华，李勇，等。1986年以来黄河下游水沙变化及河道演变分析[J]。人民黄河，2000，22(9)：10-11，16.
[5] 苏国良。黄河山东段水污染状况分析[J]。山东环境，1999(3)：44-45.
[6] 席家治，常炳炎，高传德。黄河水资源[M]。郑州：黄河水利出版社，1996.
[7] 高传德，崔树彬，李钰洪。水质参数与含沙量相关统计模式及其应用[J]。人民黄河，1986(2)：11-13.
[8] 夏军。区域水环境及生态环境质量评价——多级关联评估理论与应用[M]。武汉：武汉水利电力大学出版社，1999.
[9] 都昌杰，褚新生。环境监测水质分析基础[M]。哈尔滨：黑龙江科学技术出版社，1985.
[10] 钱宁，张仁，周志德。河床演变学[M]。北京：科学出版社，1987.
[11] Wolman M G,Miller J ude and frequency of forces in geomorphic processes[J]。Journal of Geology,1960,68:54-74.
[12] Nash D ive sediment-transporting discharge from magnitude-frequency
analysis［J］。Journal of Geology,1994,102:79-95.
[13] Orndorff R L,Whiting P ing effective discharge with S-PLUS[J]。Computers
& Geosciences,1999,25:559-565.
2024年3月9日发(作者：西门英悟)
流量和含沙量对水质参数影响分析
石 伟，王光谦
(清华大学 水利水电工程系)
摘要：为了研究近年来黄河下游水环境容量萎缩的内在机理，本文用流量级频率分析的方法，对1982～1984年黄河下游花园口、泺口水文站实测水质分析成果表中的耗氧量、溶解氧及其相应流量进行了分析，得到在平滩流量下黄河下游河道对污染物输运能力最强，而且此时其水环境容量也最大的结论。随着花园口、泺口水文站平滩流量由1958年的8000m/s和9200m/s降至1999年的4000m/s和3000m/s，黄河下游水环境容量萎缩是必然的。本文又通过统计分析得知含沙量对黄河下游水质监测数据中的耗氧量及溶解氧量影响不大。
关键词：水环境容量；平滩流量；黄河下游
基金项目：国家重点基础研究发展规划(973)项目(G19990436)；清华大学“985”校级重点科研项目软课题
作者简介：石伟(1967-)，女，内蒙古呼和浩特人，博士研究生，研究方向为水力学及河流动力学。
3333近年来，在人类大量引水、降水偏枯[1，2]等人为和自然因素作用下，黄河下游的来流明显减少，不仅表现出断流、洪水频繁、河道萎缩、主槽淤积、“小水大灾”[3，4]等水沙过程的变异，而且呈现出水环境容量萎缩、水体自净能力下降、水质恶化[3]，1996年山东黄河有1/4的时间是超标水，1997年上升为1/2，1998[5]年再度上升为近3/4，污染形势极为严峻。河流是污染物迁移载体，污染物是借助于河流中的水流和泥沙来输运的。黄河来流流量的大小影响着河流携带污染物的能力，影响着河流的水环境容量及其自净能力，夏军等人已证实增加来水量可以改善河段水环境容量[6]。进入黄河的泥沙，给黄河造成一定污染，同时由于其组成含有相当数量的粘土矿物、无机胶体和一定数量的有机胶体、有机无机复合胶体，对排入黄河的种类繁多的污染物具有显著的吸附作用，从而表现出净化水质的效应，高传德等人通过水质监测数据得到浑水高锰酸钾指数与其含沙量呈[7]正相关。本文首次采用流量级频率分析的方法研究流量对黄河下游耗氧量和溶解氧的影响，采用统计分析的方法分析含沙量对其耗氧量和溶解氧的影响，是对黄河下游水沙变化对其水质影响机理的初步探索。
1 水质实测数据的选取
耗氧量(COD)是指在一定氧化条件下，以一定的氧化剂氧化水或泥沙中还原性物质所消耗的氧化剂量。它可以作为水体有机物相对含量的指标，常用作反映水质污染程度的参数[7]。溶解氧是水中溶解态氧的浓度，是衡量水质好坏的一项重要指标。维持良好水质感官性状、水体的自净作用以及水生生物的生存需要有[8，9]氧。本文选取了1982～1984年花园口、泺口水文站实测水质分析成果表中的耗氧量、溶解氧及其相应的流量进行分析，含沙量在有相应时段洪水水文要素摘
录表时用相应时间、相应流量的含沙量，若无相应时段洪水文要素摘录表则用相应的日平均含沙量。
2 耗氧量和溶解氧量流量级频率分析
2.1 耗氧量和溶解氧量流量级频率分析方法 Wolman和Miller在1960年给出了河流有效流量的各级流量——频率分析的计算方法[10～13]：假设河流流量频率服从对数正态分布，输沙率作为流量功函数，则输沙效率即流量出现频率与输沙率乘积最大处的流量值为输沙的有效流量，就长时期内河流输沙量而言在该流量处最大[11～13]。本文借用此方法，假设河流流量频率服从对数正态分布，耗氧效率即流量出现频率与流量和耗氧量乘积最大处的流量值是长时期内河流耗氧量最大的流量，此分析方法作为耗氧量的流量级频率分析。如果将上述耗氧量换为溶解氧量，则是溶解氧的流量级频率分析。
2.2 耗氧量和溶解氧量流量级频率分析结果 图1为花园口、泺口水文站1982～1984年各级流量耗氧效率(即流量出现频率与流量和耗氧量乘积)、各级流量溶解氧效率(即流量出现频率与流量和溶解氧量乘积)，均以1000m3/s分级。从图中可以看到，花园口、泺口水文站在1982～1984年相应于耗氧效率(即流量出现频率与流量和耗氧量乘积)、溶解氧效率(即流量出现频率与流量和溶解氧量乘积)最大处的流量分别为：6348m3/s、5640～6348m3/s，二者基本相当；它们出现的天数均占总测量天数的7%.对照图2可以看到，花园口、泺口水文站在1982～1984年相应于耗氧效率和溶解氧效率最大处的流量6348m3/s和5640～6348m3/s与累计耗氧率及其累计溶解氧率曲线由陡变缓的转折点6460m3/s在数值上相近。也就是说，相应于耗氧效率和溶解氧效率最大处的流量也是黄河下游输送污染物最经济的流量。事实上，相应于耗氧效率和溶解氧效率最大处的流量6348m3/s和5640～6348m3/s与该时段花园口、泺口水文站的平滩流量相当。根据明渠流连续方程
Q=UBH (1)
式中：U为断面平均流速(m/s)，B为河面宽度(m)，H为水深(m).根据方程(1)，在平滩流量之下，当水深H和断面平均流速U随流量增大而增大时，河流宽度B保持不变[10]，因而B/H减小。当水流漫滩以后，由于滩地糙率比主槽糙率大好几倍，因而随着水位继续上涨，主槽和全断面的流速反而有所下降，当滩地水深约为主槽平均水深的35%时，断面流速达到最低值，自此以后，又随着水位的抬高而不断加大[10]。所以，相应的河道输运污染物能力在平滩流量之下随着流量的增加而增大，直至平滩流量时达到最大，在流量超过平滩流量以后，随着主槽和全断面流速的降低，输运污染物能力反而下降。
图1 花园口、泺口水文站各级流量耗氧效率及其溶解氧效率
图2 花园口、泺口水文站累计耗氧率及其
累计溶解氧率与流量关系
由于黄河下游河道在平滩流量下不仅对污染物输运能力最大，而且此时河道中的溶解氧也最多，因此在平滩流量下黄河下游河道的水环境容量为最大。在平滩流量时不仅输沙量最大[10]，而且输送污染物量也最大，为我们估算黄河下游生态需水量提供了依据，即根据平滩流量估算的黄河下游用水量不仅是输沙最经济的水量，也是输送污染物最经济的水量。
2.3 黄河下游水环境容量萎缩的原因 从表1可以看到，从1958年到1999年花园口和泺口水文站的平滩流量经历了由大到小到大又到小的过程[4]。1958年黄河下游为基本未受人类影响的自然状态，此时平滩流量较大，花园口和泺口33分别为8000m/s和9200m/s；1964年在三门峡水库的冲刷作用下，花园口平滩流量增大，为9000m3/s；1973年在三门峡水库滞洪排沙作用下，黄河下游平滩流量急剧减小，花园口和泺口的平滩流量降为3500m3/s和3100m3/s；1985年由于丰水少沙，黄河下游河道冲刷，平滩流量上升，花园口和泺口的平滩流量又上升为6900m3/s和6000m3/s；1985年以后，由于人为、自然等因素的作用，黄河下游来水来沙条件改变，河道严重萎缩，主槽淤积严重，平滩流量下降，花园口、泺口水文站平滩流量变为4000m3/s和3000m3/s.由于在平滩流量下黄河下游河道的水环境容量为最大，所以随着平滩流量的上述变化，水环境容量也会发生相应的变化。1985年以后，由于人类大量引水、降水偏枯等原因，黄河下游断流频繁，不仅使黄河下游的河床演变过程发生改变，河道萎缩、“小水大灾”、洪水频繁，而且也大大降低了其输运污染物的能力以及溶解氧的能力，表现为水环境容量萎缩、水质明显下降。与此同时，随着经济发展、人口增加，黄河流域废污水排放量与日俱增。据统计，20世纪70年代后期年排放废污水18.5亿t，80年代初为21.7亿t，90年代初已达32.6亿t，10年间增加50%以上[6]。与此相应，黄河下游必然污染情况越来越严重。
表1 典型年份花园口、泺口水文站平滩流量变化(m3/s)
时间
1958年汛后
1964年汛后
1973年汛前
花园口
8000
9000
3500
泺 口
9200
3500
3100
1985年汛前
1999年汛前
6900
4000
6000
3000
3 含沙量对耗氧量和溶解氧量的影响分析
考虑到耗氧量和溶解氧量与含沙量是同一个量纲，我们将花园口、泺口水文站1982～1984年实测水质分析成果表中的耗氧量(单位：mg/l)和溶解氧量(单位：mg/l)分别乘以1000后再分别除以含沙量(单位：kg/m3)得到两个无量纲的量，将这两个分别与耗氧量和溶解氧量有关的无量纲的量作为纵轴，均以含沙量为横轴作图，得到图3和图4.在图3中我们点绘了双曲线在第一象限的一支，其常数为1982～1984年花园口、泺口水文站所有实测耗氧量的算数平均值(约为1.46)乘以1000，横轴和纵轴同上所述；在图4中我们同样点绘了双曲线在第一象限的一支，其常数为相应溶解氧量的算数平均值(约为8.15)乘以1000，横轴和纵轴同上所述。从图3和图4中可以看出，实测点与双曲线吻合得很好，也就是说，耗氧量和溶解氧量可以看作是个常数，约等于其实测数据的算数平均值，而与含沙量无关。究其原因，与耗氧量的测定方法不无关系，因为多年来黄河上都是用去除泥沙后的清水测定耗氧量[7]，只代表了黄河水溶液中耗氧物质的含量。事实上，黄河水体中的泥沙吸附着大量的污染物，其吸附量必然与含沙量有关，也必然随着吸附量达到饱和而使水体受到二次污染。由于缺乏浑水的资料[7]，限制了我们对耗氧量和溶解氧量与含沙量关系的更深入的认识。
图4 花园口、泺口水文站溶解氧量与含沙量之图3 花园口、泺口水文站耗氧量与含沙量之比
比
(乘以1000)与含沙量的关系
(乘以1000)与含沙量的关系
4 结论
本文通过对黄河下游花园口、泺口水文站1982～1985年实测水质监测数据的各级流量频率分析，得到黄河下游河道在平滩流量下对污染物输运能力最大，此时河道中的溶解氧也最多，也就是说，在平滩流量下黄河下游河道的水环境容量最大。这一研究丰富和发展了前人的研究，即河流在平滩流量时不仅输沙能力[10]最大，而且输运污染物能力也最大。黄河下游水环境容量萎缩是平滩流量降低
的必然结果。同时，通过统计分析得知含沙量对黄河下游水质监测数据中的耗氧量及溶解氧量影响不大，即黄河下游水体水溶液的耗氧量和溶解氧基本不变。
参 考 文 献：
[1] 许炯心。人为季节性河流的初步研究[J]。地理研究，2000，19(3)：234-242.
[2] 刘昌明，成立。黄河干流下游断流的径流序列分析[J]。地理学报，2000，55(3)：257-265.
[3] 常炳炎，薛松贵，张会言，等。黄河流域水资源合理分配和优化调度[M]。郑州：黄河水利出版社，1998.
[4] 申冠卿，张晓华，李勇，等。1986年以来黄河下游水沙变化及河道演变分析[J]。人民黄河，2000，22(9)：10-11，16.
[5] 苏国良。黄河山东段水污染状况分析[J]。山东环境，1999(3)：44-45.
[6] 席家治，常炳炎，高传德。黄河水资源[M]。郑州：黄河水利出版社，1996.
[7] 高传德，崔树彬，李钰洪。水质参数与含沙量相关统计模式及其应用[J]。人民黄河，1986(2)：11-13.
[8] 夏军。区域水环境及生态环境质量评价——多级关联评估理论与应用[M]。武汉：武汉水利电力大学出版社，1999.
[9] 都昌杰，褚新生。环境监测水质分析基础[M]。哈尔滨：黑龙江科学技术出版社，1985.
[10] 钱宁，张仁，周志德。河床演变学[M]。北京：科学出版社，1987.
[11] Wolman M G,Miller J ude and frequency of forces in geomorphic processes[J]。Journal of Geology,1960,68:54-74.
[12] Nash D ive sediment-transporting discharge from magnitude-frequency
analysis［J］。Journal of Geology,1994,102:79-95.
[13] Orndorff R L,Whiting P ing effective discharge with S-PLUS[J]。Computers
& Geosciences,1999,25:559-565.