地下水资源特点及分类计算地下水允许开采量主要方.ppt

这个公式就是用来描述泉域含水体这个“转换装置”的数学模型，可利用它来预测泉的流量和评价大中型供水井区的地下水资源，以及预测矿井涌水量。 在实际运用时，首先根据输入Ut（大气降水）与输出Qt（泉流量或开采量或矿井涌水量）的实际观测序列，识别系统即确定权函数或权序列W?。确定的方法就是根据最小平方准则，即要求计算与实测的误差总和为最小，用电子计算机解权序列的线性方程组便可得到。有了权函数，便可用该系统来进行预报。用一个例子来说明其具体方法。 §2 计算地下水允许开采量的主要方法 2．实例 河北省邯郸地区南单元的岩溶水从黑龙洞泉群排泄，泉域如图10一10所示。西安煤研所和武汉地质学院应用系统理论法对黑龙洞泉群进行了具体计算与分析。其步骤如下： §2 计算地下水允许开采量的主要方法 第一步：进行条件分析： 泉流量和降水量的多年观测结果表明，泉群流量的变化规律与降水量的变化规律有明显的一致性。但泉流量峰值出现的时间滞后于雨季1一2个月。这说明岩溶蓄水盆地对降水补给的地下水有调节作用。丰水年补给充沛，使泉流量增大，其影响可达三年之久。通过水文动态分析和泉水流量自然消耗计算，一次补给的地下水，需要两年半至三年的时间才能基本泄尽，即一次降水对泉流量的影响可达30—36个月。 §2 计算地下水允许开采量的主要方法 第二步：对降水量、泉流量实测资料的处理： 在本泉域范围内有十多个水文气象站，取其中有代表性的五个站资料的算术平均值。据当地情况，把三日连续降水量在20mm和两日10mm以上者定为“有效降水量”（即有能力补给灰岩地下水的降水量），从而得出各站的逐月降水量，再求平均值用于计算。非有效降水量作零处理。 本区地下水埋深一般在30m以上，故可不考虑地下水的蒸发消耗。但人工开采量应考虑，计算时给于适当处理。 §2 计算地下水允许开采量的主要方法 第三步：确定权序列：用最小平方估计准则确定，即： 式中： N––––观测数据的个数，且N≥n＞k≥0； Qt计––––计算出的t时刻泉流量； Qt测––––实际测得的t时刻泉流量。 于是得到权序列W的线性方程组：  §2 计算地下水允许开采量的主要方法 式中： §2 计算地下水允许开采量的主要方法 解此线性方程组，就可以得到权序列Wk，Wk＋1，……，Wn。这是最小平方估计意义下的最优权序列。利用它即可用任意给定的权序列长度m（m＝n一k十1）个月的有效降水量资料，按前述流量预测公式： 就可算出泉群的预测流量来。 §2 计算地下水允许开采量的主要方法 第四步：权序列长度的确定： 权序列的长度m是指某月降水量对泉流量起影响作用的时间长短，现采用月为时间单位。如果某月降水量自k月开始影响泉流量至第n月影响消逝，那么，权序列长度m＝n一k十1。因此，要确定m必须确定k和n。 k值，可根据动态曲线分析来确定。从泉流量与有效降水量动态曲线图上可见，降水的当月对泉流量就有影响，但不甚明显，一般是当月降水对下一个月的泉流量有明显的影响，所以取k＝1。n的大小确定比k困难一些，可以通过水文地质条件的分析，结合动态资料确定。经初步确定后，再用试算法来确定，其方法如下。 §2 计算地下水允许开采量的主要方法 先将求流量的公式改写成； 式中，右边第一个和式表示降水还没有影响到泉的流量，故应为零；第三个和式表示第n十1个月以前的降水对t月泉流量的影响。很明显，如果n选择得当，则第三个和式也应该接近于零。 §2 计算地下水允许开采量的主要方法 若令 则 由于Qt′是一个趋于零的微小量，对Qt的影响很小，所以可假定为一个很小的常量，设为a，则Qt′＝|a|，于是有 §2 计算地下水允许开采量的主要方法 对上式用最小平方估计的方法进行处理后，就得到一个求权函数的正则方程，即： 式中： 用此式通过试算的办法选择合适的n，使a的绝对值达到预先指定的一个微小量（本例取0.5rn3/s，与观测平均误差相当），则相应的一列权函数即为所求。 §2 计算地下水允许开采量的主要方法 第五步：对不合理权函数的处理： 从水文地质意义而言，权函数均应大于零，但在试算中可能出现一些负值或零值。分析其原因，可能是因降水量与泉流量之间的关系不匹配，或由于泉流量的实测值有较大的测量误差，或因权序列过长，或列入一些非有效降雨量等。从以上几方面原因对资料进行分析校核，取好 k值和n值的范围，再重新试算。 在本实例计算中，权序列长度最小取到12个月，最长取到37个月，此时末端出现一个负值。按误差平方和最小的原则最后确定了权序列的长度为33个月，全为正值。 §2