黑龙江省水库站径流还原计算方法及有关问题的探讨

1、黑龙江省水库站径流复原计算方法及有关问题的讨论摘要：针对黑龙江省人类活动对河川径流的影响，系统地阐述水库站月径流量的复原计算方法。并就稻田拦蓄水量、地下水开采对河流水文效应产生的影响及径流复原计算方法进展讨论。关键词：水库站径流量耗水量复原计算我省现有水库站12处。这些站的建立，为弥补根本站网定位观测的缺乏，扩大资料搜集范围，增强资料的完好性均发挥了重要作用。但由于这些站所测得径流资料受水利工程影响，在数量和时程分配上均不能代表河流的天然状况，因此在应用这些资料进展水资源评价和工程设计时，必须将测验断面的实测径流复原为天然河川径流量。本文针对我省水利工程和水库站的特点，对水库站月径流量的复原计

2、算方法进展全面阐述，并就径流量负值处理、稻田拦蓄水量问题、地下水开采对河川径流影响等进展了讨论。1水库站河川径流复原计算公式河川径流复原计算的方法很多。主要有分项复原法、模型法、经历公式法、径流双累积法和流域蒸发差值法。而进展月径流复原计算时，通常采用分项复原法。分项复原法是根据水量平衡原理建立程度衡公式，通过计算水利工程引起的增减水量推求测验断面的天然河川径流量。计算公式如下：天然=出库+农耗+工耗库蓄+渗漏+库蒸1式中：天然复原后的天然河川月径流量；出库出库实测月径流量；农耗农田灌溉月耗水量；工耗工业用水月耗水量；库蓄水库月蓄水变量；渗漏水库月渗漏损失量；库蒸水库增加的月水面蒸发损失量。2

3、分项水量复原计算2.1农田灌溉耗水量农田灌溉耗水量是指农田灌溉引水过程中，因蒸发消耗和渗漏损失掉而不能回归到河流的水量，为渠首引水量与回归入河水量之差。在计算农田灌溉耗水量之前，须首先弄清区域用水水源、用水区域和回归水之间的关系和相对位置，来判别应复原的水量。假如灌区引水口在测验断面上游，而灌区在测验断面下游时，那么灌溉耗水量即为渠首引水量。如灌区引水口和灌区均在测验断面上游时，那么灌溉耗水量应等于渠首引水量减去综合回归水量，即：2式中：引水渠首引水量；综回农田灌溉水综合回归水量，包括田渠下渗回归水量和田渠弃水量；净灌农田灌溉净用水量净灌=净定*F实灌；田回田间下渗回归水量；净定农田灌溉净定额

4、；F实灌实际灌溉面积；E渠蒸渠系引水、输水过程中增加的蒸发损失量。由于田回和E渠蒸可互相抵消一局部，因此灌溉耗水量计算可简化为下式：农耗=净定*F实灌32.1.1水稻净灌溉定额确实定我省农田灌溉以水田为主，水田灌溉耗水量占农田灌溉耗水量总量的80%以上，因此合理确定水田灌溉用水定额，是准确计算农田耗水量的关键。根据水量平衡原理：净定=泡田+E需水+渗漏-P有效4式中：泡田泡田期用水量；E需水水稻全生育期需水量；渗漏水稻全生育期田间渗漏量；P有效水稻生长期有效降水量P有效=P降雨*有效P降雨灌溉期降水量；有效降雨有效利用系数。上式中，泡田主要和土壤前期含水量和土地平整程度有关，渗漏主要取决于土壤

5、的浸透条件，E需水受当地气候情况影响，P有效主要受降水量的大小和时程分配制约。在地形地貌、气温、浸透条件相近的同一地区，影响泡田，渗漏，E需水大小的因素相对稳定，此时，同一地区的净定将主要取决于降水量的大小及其时程分配。根据上述分析，省水科院将全省划分为东、西两大平原和南、北两大山地四个农田灌溉分区，并根据全省28个水田灌溉实验站资料，给出了不同分区P=75%的水田净灌溉定额。本文在此根底上，通过计算机回归分析确定出如下计算水田净灌溉定额的经历公式：区：净定=286.6Pi2+302.1Pi+361.65区：净定=248.5Pi2+303.1Pi+314.66区及区：净定=243.3Pi2+2

6、91.9Pi+290.87式中：净定不同降水频率下的水田净灌溉定额；Pi年降水量频率。2.1.2水田灌溉月耗水量的计算首先求出水田所在农业分区的逐年面降水量，然后经频率分析确定复原年份的年降水量频率，代入式57得年耗水量，再乘以水田灌溉水量月分配系数即为月耗水量。水田灌溉水量月分配系数如表1。表1黑龙江省水稻灌溉水量月分配系数表月份分区567890.3610.1890.2090.1500.0920.3550.2050.2270.1510.0620.3530.1880.2080.1730.0780.4120.1610.1830.1630.0822.2工业用水耗水量2.2.1工业用水耗水量确实定方

7、法工业用水耗水量为工业用水取水量与工业废水入河排放量之差。其值包括两局部：一是用水户在消费过程中被产品带走、蒸发和渗漏掉的水量，称为用水消耗量；二是工业废水在排放过程中因渗漏和蒸发而耗损的水量，称为排水消耗量。以上各项水量的程度衡关系为：工耗8式中：工耗工业用水耗水量；取水工业用水取水量；河排工业废水入河排放量；用耗用水消耗量；排耗排水消耗量。我省工业废水一般无入河排放量观测资料，且工业废水、生活污水和雨洪径流往往在城市排污口中混杂排放，很难区分，因此工业废水入河排放量资料很难直接获得。为解决这一问题，笔者根据用水指标的定义和各项水量之间的程度衡关系，提出如下入河排放量的计算方法：排放=取水*

8、1-0*1-i/1-09河排=排放*1-10将式9、10代入式8，并经整理得：工耗=取水*1-1-0*1-i/1-0*1-11式中：排放用水户厂区工业废水排放量0调查起始年用水消耗率0=用耗/取水；0、i调查起始年、复原年的工业用水重复利用率；工业废水在入河过程中的渗漏、蒸发损失率=排耗/排放，明渠输水取0.200.30，管道输水取0.150.25。2.2.2工业用水月耗水量的计算将调查起始年的0、0及复原计算年份的i、带入式11，求出工业用水年耗水量再平均分配到年内各月，即可求出工业用水月耗水量。2.3水库渗漏损失量水库渗漏量分为坝基渗漏、库底渗漏和库岸侧渗三局部。在坝基渗漏较大的水库上，可

9、利用坝下反滤沟实测流量资料推算水库渗漏量。假如没有实测流量资料，可采用以下方法进展计算：渗漏=入库+降雨蓄变-出库-库蒸12式中：入库入库站实测月径流量；降雨库区水面月降雨量；蓄变水库月蓄水变量。在缺乏实测入库径流量的水库上，可选用封冻期入库流量为零或入库流量接近零的月份推算月渗漏水量。在河流封冻期降雨为零，由于枯水月份入库和库蒸均很小，且可互相抵消一局部，因此式12可进一步简化为：渗漏=蓄变-出库13应用式13计算水库月渗漏量，并与水库月平均水位点绘相关曲线，再以月平均水位在关系上查读月渗漏水量。2.4水库增加的水面蒸发损失量水库蓄水以后，库区由陆地蒸发转变为水面蒸发，其月净增的水面蒸发损失

10、量为：库蒸=1-E陆年/E601年*E601月*A月14式中：E陆年库区年陆地蒸发量；E601年库区年E601蒸发器蒸发量；E601月库区月E601蒸发器蒸发量；A月库区月平均水面面积。2.4.1E陆年的计算我省无陆面蒸发量实验资料。以往计算E陆年时，一般采用年降水量与年径流深之差P-R代替。其中年径流深R往往借用相邻地区资料，或以实测出库水量和水库蓄水变量代数和近似代替。由于以P-R代替E陆年是以忽略区域年蓄水量为前提的，方法本身不符合程度衡原理，加之R值多为近似值，因此用上述方法计算的E陆年会产生较大的误差。我省在第一次水资源评价时曾提出利用水热平衡法计算多年平均陆面蒸发量的方法，可应用这

11、一方法计算E陆年。1山丘区E陆年的计算E陆年=P降雨*10-a+0.1553/LP152平原区E陆年的计算E陆年=P降雨*-0.246+0.642/LP16表2山丘区a取值条件表流域平均高程()a2003000.3543004500.4044507500.4277500.475以上公式中：P降雨库区年降水量；太阳辐射平衡值；L蒸发潜热；a随流域平均高程变化的系数，取值条件如表2。2.4.2水库月蒸发损失量的计算利用式15、16求出E陆年后，再将其同E601年、E601月、A月一并代入式14即可求出水库月蒸发损失量。3径流复原计算中有关问题的讨论3.1稻田拦蓄水量问题3.1.1稻田拦蓄水量复原计

12、算的必要性稻田拦蓄水量是指稻田拦蓄利用产流降雨而引起河川径流减少的水量。目前，对于是否进展该水量复原计算，专家和学者一直持两种不同态度：一种观点认为，稻田拦蓄水量应归于下垫面条件变化对河川径流的影响，本着向后复原的原那么，不进展该水量的复原计算；另一种观点那么认为，稻田用水属于水利工程的延续，稻田拦蓄的水量就是产流降雨，稻田对降水的拦蓄利用，使河川径流在数量和时程上均发生显著变化，因此应进展该水量的复原计算。笔者倾向于后一种观点，主要理由有以下几条：1每次降雨过程中稻田周围的地面径流通过渠系直接聚集于稻田，并直接拦蓄了稻田面积上的降雨，因此稻田实际起到和水库一样的拦蓄作用，假如稻田不拦蓄，这局

13、部水量就会汇流于河道形成河川径流。2和水库不同的是，稻田拦蓄的水量除局部通过田间下渗回归河流外，其余大局部水量作为水资源被利用以作物蒸发和稞间蒸发的形式被就地消耗，使下垫面的产水量明显减少。3稻田拦蓄的水量数量可观。以我省东部的挠力河为例，农田区和一般平原区相比，同量级降水情况下，前者的径流深较后者减少2050，因此稻田拦蓄作用对河川径流的影响是显著和不容忽略的。3.1.2稻田拦蓄水量的计算方法稻田拦蓄作用下的产汇流机制相当复杂，计算拦蓄水量目前尚无成熟方法。本文以降水量为输入，产流量为输出，在忽略产汇流过程的前提下，用以下方法估算：拦蓄=P降雨*F实灌*径流1-有效*1-田间17式中：拦蓄因

14、稻田拦蓄利用而减少的径流量；P降雨灌溉期月降水量；田间灌溉水田间回归系数，取0.150.20；径流天然状态下的降雨径流系数，可由稻田建成前的月降雨径流关系曲线查算；表3降雨有效系数表降水量()有效50501501000.81502000.752000.7有效降雨有效利用系数，按表3取值。3.2地下水开采对河川径流的影响目前，在河川径流复原计算中，一般只考虑地表水开发利用对河川径流的影响，而不考虑地下水对河流产生的水文效应。是否进展复原计算，争论已久，至今尚无定论。地下水开采对河流的水文效应的影响，可概括为增水或减水两个方面。一方面，地下水开采产生地下水位下降后阶段性或持续性下降，地下水和地表水

15、补排关系发生变化，造成河流激化补给地下水，使河川径流量较天然状态有所减少；另一方面，经开采利用的地下水会以废污水的形式直接排入河道，使河川径流量增加。以福利屯站为例，安邦河定国山水库至测验断面的区间内，双鸭山市区每年有近2000104t的地下水废污水排入河流，其年排放量占80年代以来最枯年径流量的55%，月排放量占月最枯径流量的98%。地下水开采对河川径流的影响，可用下式估算：地下=B*K2H2T2-1H1T1-地排18式中：地下地下水开采使河川径流增加或减少的水量；B河道长度；K含水量浸透系数；H1、H2地下水开采前、开采后的地下水水力坡度；1、2地下水开采前、开采后的含水层影响带有效厚度；

16、T1、T2地下水开采前、开采后的河道渗漏补给时间；地排地下水废污水入河排放量。当地下较小，地下可忽略不计，不进展该水量的复原计算。当地下较大，且地下为正值时，将地下按月径流占补给期河流总径流量的百分比分配到各月，即可求出月复原水量；当地下较大且地下为负值时，将地下平均分配到年内各月，即为月复原水量。3.3径流量负值的处理在实际工作中，在封冻期的一些月份12月2月常常出现复原后的径流量为负值的不合理现象。对负值月份采用以下处理方法：1首先通过频率分析确定负值年份的径流量丰枯程度，然后选择与其接近的丰、平、枯年份作为典型年，直接移用典型年同一月份的径流量，将“负值处理为“正值。为保持年径流总量不变，需将“负值月份“借用的水量，在其后的45月份按比例扣除。2采用亏水比例法。选择降水产流量甚小的枯水时段，根据亏水比例公式t=0*K，先用不受工程影响的天然月径流分析K值。然后以受工程影响、径流复原计算后的汛末一个月的月径流量0为基数，再乘以K值，分别求出汛末到翌年2月份的各月径流量，用以代替“负值月份的径流量。4结语本文