岔巴沟流域不同淤积状态下淤地坝淤积深度研究.docxVIP

? ? 岔巴沟流域不同淤积状态下淤地坝淤积深度研究 ? ? 马 力, 张秋芬, 倪用鑫, 李 勉, 吕锡芝 (1.黄河水利科学研究院, 郑州 450003; 2.河南省黄河流域生态环境保护与修复重点实验室, 郑州 450003) 目前，黄土高原地区水土保持作用有明显提高，但遭遇特大暴雨时仍会形成高含沙大洪水，引发重大灾情，严重威胁着下游人民的生活和生产安全[1]。淤地坝具备高效的保持水土、蓄水拦沙功能，在黄土高原地区分布广泛[2]，在暴雨洪水过程中，淤地坝具有明显的削峰和滞洪作用[3]。淤地坝在拦蓄泥沙的同时，也赋存了大量的侵蚀产沙信息，了解淤地坝泥沙沉积特征是坝地泥沙输移过程及规律研究的基础和前提[4-6]。许多学者利用坝地泥沙养分[5]、粒径[6]、干容重[7]、核素示踪等[8]指标方法分析淤地坝沉积泥沙信息，研究流域暴雨洪水侵蚀量、侵蚀产沙速率、沉积泥沙来源等[9-14]。淤地坝通常由坝体，排水设施或溢洪道构成，泥沙呈动态淤积状态，较难量测和计算，而闷葫芦坝由于无排水排沙出口，能够很好地揭示坝地的泥沙淤积深度和泥沙来源情况，是大部分学者关注的重点。相比之下，针对其他类型淤地坝泥沙淤积深度的研究则相对薄弱。然而通常情况下，流域内淤地坝类型分布多样，闷葫芦坝数量占比较少，单凭闷葫芦坝的研究结果难以准确反映小流域整体侵蚀产沙状态，分析不同淤积状态下的淤地坝淤积深度，有利于从多角度反映流域洪水过程泥沙输移状态，从而进一步揭示降雨对流域侵蚀的作用机理。因此，本研究以2017年“7·26”暴雨中心区域岔巴沟小流域为研究对象，通过对不同淤积状态下淤地坝内不同断面位置次洪水泥沙沉积层的量测，研究了次洪水条件下坝地泥沙沉积层淤积深度特征，为淤地坝泥沙沉积和输移过程的研究提供理论支撑。 1 研究区概况 岔巴沟(位于东经109°47′，北纬37°31′)是大理河流域的一级支流、无定河流域二级支流，属于黄土丘陵沟壑区，流域面积205 km2，沟道长26.5 km，流域平均宽度7.22 km，沟道密度1.05 km/km2，流域出口水文站为曹坪水文站，控制流域面积为187 km2。2017年7月25日20时—26日8时，陕北榆林地区无定河流域发生特大暴雨洪水，暴雨中心降雨量达252.3 mm。本研究调查了暴雨中心岔巴沟流域的大、中、小型淤地坝共148座(其中骨干坝31座，中型坝65座，小型坝52座)，选取有明显淤积的淤地坝38座，根据坝体排水排沙状态和是否具备放水工程等条件，将淤地坝分为4类：(1) 淤满坝(3座)。在暴雨洪水前已淤满，即超过拦泥库容的淤地坝；(2) 水毁坝(7座)。暴雨洪水过程中被侵蚀而产生垮塌或形成钻水洞的淤地坝；(3) 完好坝(19座)。工程完好，有放水设施(溢洪道、竖井等)且未被堵塞的淤地坝；(4) 堵塞坝(9座)。工程完好，无排水设施，或排水设施在洪水前被堵塞的淤地坝。淤地坝具体现状及分布情况见图1和图2。 图1 岔巴沟流域淤地坝类型分布 图2 不同淤积状态淤地坝示意图 2 研究方法与数据来源 2.1 泥沙淤积体积及坝地面积计算 采用野外调查定位与卫星影像观测相结合，对流域内淤地坝分布进行排查并编号(坝号)。根据淤地坝内地形及汇水情况，分别对坝地主沟、支沟设置若干淤积断面，采用激光测距仪测量并记录断面长度，分区块计算坝地淤积面积(区块合并后的面积称为坝地面积)。每个断面人工开挖3个探坑，将暴雨洪水形成的沉积旋回深度作为“7·26”洪水形成的泥沙沉积厚度。淤地坝的泥沙淤积体积V计算公式为： (1) (2) 式中:Si为所调查淤地坝坝地相邻两个断面之间的淤积面积(m2);hi为第i个断面的平均淤积厚度(m);n为断面个数;rj为第j个探坑的淤积厚度(m)。 坝地面积与流域坝控面积之比(以下简称B值)的计算公式为： B=b/S (3) 式中:b为淤地坝单坝的坝地面积;S为淤地坝单坝的坝控流域面积。 2.2 拦沙模数计算 将淤地坝坝址的GPS地理坐标转换到具有大地坐标系的地图上，利用ArcMap空间分析工具，将坝址作为倾泻点，计算出坝控流域面积。将每座淤地坝泥沙淤积量平摊到该淤地坝的坝控面积上，得到单坝拦沙模数[15]。计算公式为： MY=1.30Y/A (4) 式中:MY为单坝拦沙模数(t/km2);Y为淤积体积(m3);A为坝控面积(km2);1.30为土壤容重(g/cm3)。 2.3 雨强数据 降雨数据来源于《中华人民共和国水文年鉴》[16]，其中岔巴沟流域包含12个雨量站和1个水文站，最小降雨时段数据为2 h，坝控流域雨强根据雨量站2 h最大降雨量，利用GIS空间分析工具进行反距离加权插值计算获得。 3 结果与分析 3.1 坝地淤积深度对比关系 岔巴沟流域不同淤积状态下淤地坝坝地淤积深度差异较为明显(图3)。单坝