一种岩溶流域耦合水文模型预测方法

技术领域

本发明属于水文学流域水文模型水量预测预报技术领域，具体涉及一种岩溶流域耦合水文模型预测方法。

背景技术

岩溶流域由于其特殊的含水系统结构和含水介质的极不均一性，其水文运动规律与一般流域不同。经研究发现西南岩溶地区的地下河系统与地表河流具有很多相似特点；主要表现在一方面构成地下河系统的岩溶管道结构(网)与地表水系网络十分类似，另一方面地下河系统对降雨过程和大小响应迅速，也具有地表河流流量大、流速快、比降大等特点，且岩溶泉水流量衰减曲线呈现多段式曲线，对应多重岩溶含水介质特征，地下河与地表水之间有时候往往形为一体、相互转换，共同构成当地水系网络。

在资料缺乏且不便于采集资料数据的岩溶地区，无法为岩溶水系统涌水量提供有效的监测或预测数据，不便于在岩溶地区施工。

发明内容

针对上述技术问题，本发明的实施例提供了一种在资料缺乏地区可监测计预测岩溶地区其流域耦合水文模型的建立方法及其应用。

为解决上述技术问题，本发明的实施例提供了一种岩溶流域耦合水文模型预测方法，包括以下步骤：

(1)根据遥感、地理信息系统、全球定位系统在岩溶区采集的信息，下载矢量等高线制作岩溶区流域的数字高程模型，并对得到的模型进行处理得到岩溶地区流域的原始基础数据；

(2)将得到的原始基础数据建立Topmodel模型，并输入地形指数、降雨数据、蒸发数据计算地表水进入地下水的总量；

(3)采用高分辨率流量监测装置监测岩溶地区流域的流量，识别流域的岩溶含水介质特征，并建立多重岩溶含水介质水箱模型，通过所述水箱模型计算地下流域裂隙/管道的出流系数值；

(4)根据得到的地表水补给、地下流域裂隙/管道的出流系数值，利用上述Topmodel模型、Tank模型耦合模拟得到初步的耦合水文模型，经对参数率定后得到校正后的耦合水文模型；

(5)将得到的耦合水文模型应用于模拟岩溶区流域的降雨洪水过程曲线，预测岩溶地区暗河洪水量。

与相关技术比较，本发明的实施例提供的技术方案带来的有益效果是：本发明的岩溶流域耦合水文模型的建立方法及其应用，建立耦合水文模型方法简单，建立的耦合水文模型，既可模拟地表洪水过程曲线，又可以通过数据模拟地下水系统的多重含水介质洪水退水衰减过程，实现了岩溶地区地表水与地下水量耦合计算，可很好地模拟岩溶地区孔隙介质、裂隙介质、小管道、大管道、暗河等多重岩溶含水介质的分段式衰退规律，实现准确地预报岩溶地区大暴雨时暗河洪水衰退，具有更大的洪水预报应用价值。

附图说明

图1是本发明实施例的耦合水文模型建立方法流程图；

图2是本发明实施例的流量衰减曲线及分段示意图；

图3是本发明实施例的复式矩形堰流量监测站结构示意图；

图4是本发明实施例建立的耦合水文模型的物理概念示意图；

图5是本发明实施例孔家湾子流域数字高程模型图；

图6是本发明实施例子流域地形指数分布图；

图7是本发明实施例孔家湾子流域洪水量过程模拟曲线图；

其中：复式矩形流量堰21、监测装置22、钻孔水位监测站3、左侧堰体4、右侧堰体5、基底6、复式矩形堰口7、第一层薄壁堰口71、第二层砌砖堰口72、不锈钢板73、水位监测装置孔8、不锈钢管81、小孔82。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。

实施例一

请参考图1，本发明的实施例提供了一种岩溶流域耦合水文模型预测方法，包括以下步骤：

(1)根据遥感、地理信息系统、全球定位系统在岩溶区采集的地理地貌信息及图像，根据采集的信息及图像处理得到相应的矢量等高线，根据矢量等高线制作岩溶区流域的数字高程模型，并对得到的数字高程模型进行处理得到岩溶地区流域的原始基础数据；在岩溶区，人工不便于采集地理形势信息，利用遥感、地理信息系统、全球定位系统采集的地理地貌信息及图像，提高了采集信息的准确性，降低了劳动强度；

所述原始基础数据包括地形指数数据、地形指数概率分布曲线；

(2)将得到的原始基础数据建立Topmodel模型，并输入地形指数、降雨数据、蒸发数据计算地表水进入地下水的总量Q_q；其中降雨数据、蒸发数据根据气象监测的数据得到；

(3)采用高分辨率流量监测装置监测岩溶地区流域的流量，识别流域的岩溶含水介质特征，并建立多重岩溶含水介质水箱模型，通过所述水箱模型计算地下流域裂隙/管道的出流系数；

所述岩溶含水介质特征通过暗河流量衰减方程(1)中的衰减系数α识别；

其中：衰减期任意时刻—t；衰减开始时刻—t₀；t时刻相应的流量—Q_t；t₀时刻相应的流量—Q₀；衰减系数—α；

得到衰减系数方程(2)为：

其中，α的范围为n×10^-1～n×10^-4；

参照附图2，由于岩溶含水介质的高度非均质性，根据衰减系数值将岩溶水的衰减动态分解为若干个衰减段，判断岩溶水水流状态；

AB段：曲线较陡，α值较大，在n×10^-1～n×10^-2之间，表明在流量衰减初期各种泄水通道的总和，但水量主要来自大型的岩溶管道和地下暗河或洞穴的快速排泄，地下水流速大，流量衰减快，持续时间较短，仅几天至十几天，水流常呈现紊流状态；

BC段：曲线坡度较AB段有所变缓，α值也相应减小，一般在n×10^-2～n×10^-3之间；表明相应的反映来自大型岩溶管道和洞穴的水量有限，主要排泄的是岩溶化的大裂隙及其他溶洞裂隙系统的水，由于流量衰减速度减小，使得该段可以保持较长期的衰减趋势；

CD段：坡度变得更缓，其α值也更小，大多在n×10^-3～n×10^-4之间，表明地下水水力坡度大大减缓，以层流为主，主要是排泄细小裂隙、层间裂隙、节理中的储水；由于排泄速度进一步减慢，因而延长期比前两亚动态更长；

DE段：曲线趋于水平，α值最小，一般为n×10^-4数量级甚至更小，相当于较稳定的排泄充满于细小裂隙系统及洞穴充填物孔隙中的水；

参照附图3，其中，流量Q通过岩溶区的复式矩形流量监测站监测数据得到；具体地，所述复式矩形堰流量监测站包括复式矩形流量堰21及监测装置22，所述复式矩形流量堰21包括左侧堰体4、右侧堰体5和基底6，所述左侧堰体4与右侧堰体5为对称结构，所述左侧堰体4包括紧邻河沟岸的矩形结构及呈阶梯型的矩形结构；所述左侧堰体4、右侧堰体5和基底6形成复式矩形堰口7，所述监测装置22位于所述复式矩形堰口7的上游；所述监测装置22的监测步长为5min，水位监测的精度为1mm，通过复式矩形流量堰21的监测的水位变化得到各级流域系统的水流量、降雨量数据信息；

所述复式矩形堰口7的上游挖设有水位监测装置孔8，所述水位监测装置孔8中埋设不锈钢管81，所述不锈钢管81的管身上均匀设置小孔82，所述不锈钢管81中安装所述监测装置22；所述不锈钢管81的外径为5cm，壁厚4mm；

所述复式矩形堰口7包括第一层薄壁堰口71、第二层砌砖堰口72，所述第一层薄壁堰口71采用不锈钢板73制作成薄壁堰口；所述不锈钢板的厚度为3mm～5mm；所述第二层砌砖堰口72的厚度为9cm～11cm；

流量根据所述复式矩形流量堰21的堰前水位与最大堰顶高度之间的关系采用不同的计算方式；当所述堰前水位小于最大堰顶高度采用公式(3)计算流量，所述堰前水位大于最大堰顶高度采用公式(4)计算流量；

其中，Q为流量，单位为每秒每立方米；m为流量系数；B₁为所述第一层薄壁堰口71的宽度，单位为米；B₂为所述第二层砌砖堰口72的宽度，单位为米；g为重力加速度；H为堰前水位，单位为米；P₁为所述第一层薄壁堰口71对应的小堰上游坎高，单位为米；h₁为最大堰顶高度，单位为米；P₂为所述第二层砌砖堰口72对应的大堰上游坎高，单位为米；

通过所述复式矩形流量堰监测得到的流量值通过公式(3)或公式(4)计算得到暗河流量衰减系数，即地下流域的管道/裂隙出流系数值；

(4)统计气象数据，根据得到的地表水进入地下水总量Q_q、地下流域的管道/裂隙出流系数值及统计的气象数据，利用上述Topmodel模型、水箱模型耦合模拟得到初步的耦合水文模型，经对参数率定后得到矫正后的耦合水文模型；

具体地，统计岩溶区流域的降雨量、降雨强度数据，得到各流域的汇流滞后时间，对Topmodel模型进行参数率定，根据降雨量、降雨强度得到洪水期的流量衰减规律、岩溶水系统含水介质特征及各储水空间所占比例，对所述水箱模型的出流系数进行率定，并将率定后的参数输入所述耦合水文模型中得到校正后的耦合水文模型；

利用地下暗河或流量衰减曲线分析含水层储水空间的性质及其各占总储水量的比例；

由dV＝-Q_tdt(5)

当t＝0时V＝0(6)

故

如果衰减曲线是由几个亚动态迭加，则应为其积分之和，各亚动态贮水量(V_i)对总贮水量(V₀)的百分比为：

得到所述耦合水文模型为：

Q_q＝S*ξ*z(13)

其中，Q_总为地下暗河流量模拟值，Q_q为地表水进入地下水总量；S为流域面积系数；z为水位参数；ξ为流域尺度参数；α₁为微小裂隙/孔隙出流系数；z₁为微小裂隙/孔隙含水介质空间体积数量参数；α₂为裂隙/小管道出流系数；z₂为裂隙/小管道含水介质空间体积数量参数；α₃为大管道/暗河出流系数；z₃为大管道/暗河含水介质空间体积数量参数。所述流域面积系数S、流域尺度参数ξ、微小裂隙/孔隙含水介质空间体积数量参数z₁、裂隙/小管道含水介质空间体积数量参数z₂、大管道/暗河含水介质空间体积数量参数z₃均通过监测统计数据获得。

本发明实施例建立的耦合水文模型参数较少易于计算，与相关技术中的水文模型或传统的数值模拟方法比较，克服了传统数值模拟方法在空间结构上高度的非均质性、空间上的不均匀性和分散性的岩溶地区运用时的不足与适应性差的问题，更加适合用于水文监测资料缺乏的研究区；在洪水预测预报过程中极大地降低了劳动强度和经济成本；将建立的耦合水文模型应用于岩溶地区，模拟岩溶地区暗河洪水流量，很好地展现了孔隙介质、裂隙介质、小管道、大管道、暗河等多重岩溶含水介质的分段式衰退规律，使得在预测预报岩溶地区大暴雨时暗河洪水衰退规律更加准确，具有更大的洪水预报应用价值。

(5)将得到的耦合水文模型应用于模拟岩溶区流域的降雨洪水过程曲线，预测岩溶地区暗河洪水量。

实施例二

参照附图5～7，根据本发明实施例耦合水文模型预测方法，对孔家湾地区的地形地貌采集信息及图像，获得数字高程模型图及地形指数分布曲线图；根据计算得到的相关数据、统计数据和监测数据，采用上述耦合水文模型得到孔家湾子流域洪水量过程模拟曲线图预测孔家湾子流域的洪水流量。

在本文中，所涉及的前、后、上、下等方位词是以附图中零部件位于图中以及零部件相互之间的位置来定义的，只是为了表达技术方案的清楚及方便。应当理解，所述方位词的使用不应限制本申请请求保护的范围。

在不冲突的情况下，本文中上述实施例及实施例中的特征可以相互结合。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。