一种采煤沉陷区蓄滞洪效果的定量计算方法
【专利摘要】本发明公开了一种采煤沉陷区蓄滞洪效果的定量计算方法，获取与所述采煤沉陷区相连的河道天然来水日平均流量，根据河道天然来水日平均流量计算河道日天然来水量；获取所述采煤沉陷区设计洪水位对应的库容，获取所述采煤沉陷区死库容，获取与所述采煤沉陷区相连的河道汛限水位对应的洪量；根据所述河道天然来水日平均流量、设计洪水位对应的库容、死库容、河道汛限水位对应的洪量通过水量平衡关系定量计算采煤沉陷区蓄滞洪效果。本发明免去长时间大量开展各种观测实验才能推出采煤沉陷区的蓄滞洪效果，节省了工作成本，降低了研究难度，更加直观地展现了给定条件下采煤沉陷区蓄洪削峰的效果，对缓解当地的洪涝灾情有积极作用。
【专利说明】
-种采煤沉陷区蓄滞洪效果的定量计算方法
技术领域
[0001] 本发明设及水利技术领域，尤其设及一种采煤沉陷区蓄滞洪效果的定量计算方 法。
【背景技术】
[0002] 采煤沉陷区是由于煤层开采形成的采空区破坏了周围岩体原始应力平衡，造成地 表沉陷，使平原矿区地表大面积积水，不仅减少了可用±地资源，对地面设施、生态环境和 当地居民生活也有不同程度的影响。
[0003] 部分采煤沉陷区周边气候复杂多变，年际降水量变化大且年内分配极不均匀，形 成了降水少易旱、有降水易溃、强降水易洪的典型区域旱溃特征;利用采煤沉陷区进行蓄滞 洪，对建立科学防洪调度体系，解决沉陷区周边洪溃问题，缓解区域防讯减灾的压力有都具 有积极的作用。
[0004] 现阶段，国外对采煤沉陷区的研究主要集中于对±地复星及综合治理，主要是将 沉陷区法地进行回填、整平、恢复±壤结构等工作后，把沉陷区开发为林地、草地、农地、生 物栖息地等有益自然环境的场所；国内部分采煤沉陷区多用来作为储水水库，综合治理方 法、地质环境影响、沉陷区水资源综合利用是学者们研究的重点，针对沉陷区水资源配置的 研究较多，采煤沉陷区蓄滞洪效果方面的研究较少，特别是采煤沉陷区对洪水的削峰效果 定量计算化研究并不多见。
【发明内容】
[0005] 本发明目的就是为了弥补已有技术的缺陷，提供一种采煤沉陷区蓄滞洪效果的定 量计算方法，便于评估不同沉陷情景下沉陷区对各级别洪水的蓄滞效果。
[0006] 本发明是通过W下技术方案实现的：
[0007] -种采煤沉陷区蓄滞洪效果的定量计算方法，
[000引获取采煤沉陷区相连的河道天然来水日平均流量:所述河道天然来水日平均流量 为24小时内通过河流某一断面的水量的算术平均值；
[0009] 获取采煤沉陷区设计洪水位对应的库容:所述设计洪水位对应的库容为当遇到设 计标准洪水时，采煤沉陷区经调洪后达到的最高水位W下的蓄水容积；
[0010] 获取采煤沉陷区死库容:所述死库容为采煤沉陷区在正常运用情况下，允许消落 到的最低水位W下的蓄水容积；
[0011] 获取采煤沉陷区相连的河道讯限水位对应的洪量:所述讯限水位对应的洪量为河 流按照W往讯期的最大水流量和河岸的最大承受水量设定的水位对应流量与时间的乘积， 时间按24小时计算；
[0012] 根据所述的河道天然来水日平均流量、设计洪水位对应的库容、死库容、河道讯限 水位对应的洪量通过水量平衡关系定量计算采煤沉陷区蓄滞洪效果。
[0013] 所述的采煤沉陷区还包括:河道区、沉陷区、供水区；
[0014] 所述的采煤沉陷区蓄滞洪效果包括:沉陷区蓄水效果、沉陷区理论下泄量、供水区 取用水量和沉陷区库容调节变化量。
[0015] 所述的水量平衡关系为：
[0016] Vt-Vy = Qr-Qu-Qre
[0017]其中：Vt为当日沉陷区蓄水量，Vy为上一日沉陷区蓄水量，Qr为河道日天然来水量， Qu为沉陷区当日用水量，Qre为河道下泄量。
[001 引 Qr二化? A t
[0019] 其中qr为河道天然来水日平均流量，At为时间步长，按24小时计算；
[0020] 所述的沉陷区蓄水效果为：
[0021] 当Qr《Qf时，所述沉陷区不起蓄所述河道区来水，沉陷区理论下泄量不为零；
[0022] 当Qr〉Qf时，所述沉陷区起蓄所述河道区来水，沉陷区理论下泄量为零；
[0023]
[0024]
[0025]
[0026]
[0027]
[002引
[0029]
[0030] 其中：qf为河道区防洪讯限水位对应天然来水日平均流量，Qf为河道区防洪讯限水 位对应的日天然来水量，Vmax为沉陷区最大库容量，Vmin为沉陷区死库容量，如为沉陷区下泄 量，Qmu为供水区当日最大取用水量，Q为供水区理论供水量，V为沉陷区即时库容。
[0031] 通过对所述采煤沉陷区上游的水文站对洪水的规模进行预报，当讯期来水即将超 过所述河道区讯限水位对应的洪量时，打开所述采煤沉陷区与河道区上游连接的闽一，将 讯期洪水蓄入所述采煤沉陷区内，与此同时关闭用于连接所述供水区的闽=，待所述采煤 沉陷区内的库容即将达到设计洪水位对应的库容时，打开与所述采煤沉陷区W及河道区下 游连接的闽二。当洪峰过境时，关闭闽一和闽二，打开闽=，将所述采煤沉陷区内的蓄水有 序供给生活生产所用，在下一个洪峰到来之前科学调度调节库容，W备拦蓄洪峰。
[0032] 本发明的优点是:本发明免去长时间大量开展各种观测实验才能推出采煤沉陷区 的蓄滞洪效果，节省了工作成本，降低了研究难度，更加直观地展现了给定条件下采煤沉陷 区蓄洪削峰的效果，对缓解当地的洪溃灾情有积极作用，同时可方便地方矿业集团和水务 部口制定利用采煤沉陷区蓄滞洪进行防讯调度的方案，减轻洪溃灾害对当地的破坏。
【附图说明】
[0033] 图1为本发明实施例平面示意图。图2为本发明实施例调洪过程图。
[0034] 图3为本发明实施例2015年沉陷情景下某沉陷区5年一遇调洪过程图。
[0035] 图4为本发明实施例2020年沉陷情景下某沉陷区5年一遇调洪过程图。
[0036] 图5为本发明实施例2030年沉陷情景下某沉陷区5年一遇调洪过程图。
[0037] 图6为本发明实施例2015年沉陷情景下某沉陷区10年一遇调洪过程图。
[0038] 图7为本发明实施例2020年沉陷情景下某沉陷区10年一遇调洪过程图。
[0039] 图8为本发明实施例2030年沉陷情景下某沉陷区10年一遇调洪过程图。
[0040] 图9为本发明实施例2015年沉陷情景下某沉陷区20年一遇调洪过程图。
[0041] 图10为本发明实施例2020年沉陷情景下某沉陷区20年一遇调洪过程图。
[0042] 图11为本发明实施例2030年沉陷情景下某沉陷区20年一遇调洪过程图。
[0043] 图12为本发明实施例2015年沉陷情景下某沉陷区50年一遇调洪过程图。
[0044] 图13为本发明实施例2020年沉陷情景下某沉陷区50年一遇调洪过程图。
[0045] 图14为本发明实施例2030年沉陷情景下某沉陷区50年一遇调洪过程图。
[0046] 图15为本发明实施例调蓄前后蓄滞洪量对比示意图。
【具体实施方式】
[0047] -种采煤沉陷区蓄滞洪效果的定量计算方法，
[0048] 获取采煤沉陷区相连的河道天然来水日平均流量:所述河道天然来水日平均流量 为24小时内通过河流某一断面的水量的算术平均值；
[0049] 获取采煤沉陷区设计洪水位对应的库容:所述设计洪水位对应的库容为当遇到设 计标准洪水时，采煤沉陷区经调洪后达到的最高水位W下的蓄水容积；
[0050] 获取采煤沉陷区死库容:所述死库容为采煤沉陷区在正常运用情况下，允许消落 到的最低水位W下的蓄水容积；
[0051] 获取采煤沉陷区相连的河道讯限水位对应的洪量:所述讯限水位对应的洪量为河 流按照W往讯期的最大水流量和河岸的最大承受水量设定的水位对应流量与时间的乘积， 时间按24小时计算；
[0052] 根据所述的河道天然来水日平均流量、设计洪水位对应的库容、死库容、河道讯限 水位对应的洪量通过水量平衡关系定量计算采煤沉陷区蓄滞洪效果。
[0053] 所述的采煤沉陷区还包括:河道区、沉陷区、供水区；
[0054] 所述的采煤沉陷区蓄滞洪效果包括:沉陷区蓄水效果、沉陷区理论下泄量、供水区 取用水量和沉陷区库容调节变化量。
[0055]所述的水量平衡关系为：
[0化6] Vt-Vy = Qr-Qu-Qre
[0057]其中：Vt为当日沉陷区蓄水量，Vy为上一日沉陷区蓄水量，Qr为河道日天然来水量， Qu为沉陷区当日用水量，Qre为河道下泄量。
[0化引 Qr二化? At
[0059] 其中qr为河道天然来水日平均流量，A t为时间步长，按24小时计算；
[0060] 所述的沉陷区蓄水效果为：
[0061] 当Qr《Qf时，所述沉陷区不起蓄所述河道区来水，沉陷区理论下泄量不为零；
[0062] 当Qr〉Qf时，所述沉陷区起蓄所述河道区来水，沉陷区理论下泄量为零；
[0063]
[0064]
[00 化]
[0066]
[0067]
[006引
[0069]
[0070] 其中：Qf为河道区防洪巧限水位对应天然来水日平均流量，Qf为河道区防洪巧限水 位对应的日天然来水量，Vmax为沉陷区最大库容量，Vmin为沉陷区死库容量，如为沉陷区下泄 量，Qmu为供水区当日最大取用水量，Q为供水区理论供水量，V为沉陷区即时库容。
[0071] 图1为本发明采煤沉陷区蓄滞洪效果定量计算实施例平面示意图，如图1所示，所 述采煤沉陷区包括河道区1、沉陷区2、供水区3 =大区域，本实施例的方法可W包括：
[0072] 通过河道区上游水文站预报W及流量计4测定获取与所述采煤沉陷区相连的河道 天然来水日平均流量qr见附表1，获取所述采煤沉陷区设计洪水位对应的库容Vmax,获取所 述采煤沉陷区死库容Vmin,获取与所述采煤沉陷区相连的河道讯限水位对应的洪量Qf;
[0073]
[
[0075] 表1实施例一河道天然来水日平均流量表
[0076] 根据所述河道天然来水日平均流量qr计算为河道日天然来水量Qr。
[0077] 当河道区河道日天然来水量Qr即将超过防洪讯限水位对应的日天然来水量Qf时打 开所述沉陷区与河道区上游连接的闽一5,将讯期洪水蓄入所述沉陷区内，与此同时关闭用 于连接所述供水区的闽=7，此时河道区下泄量为Qre ;
[0078] 待所述沉陷区内的库容V即将达到设计洪水位对应的库容Vmax时，打开与所述沉陷 区W及河道区下游连接的闽二6进行泄水，所述沉陷区泄水量为Q。，；
[0079] 当洪峰过境，即河道区河道日天然来水量Qr低于防洪讯限水位对应的日天然来水 量Qf时，关闭闽一5和闽二6,打开闽=7,将所述沉陷区内的蓄水有序排至所述供水区，供生 产生活所用，相应的供水区取用水量为Q。
[0080] 待所述沉陷区内的库容V即将达到死库容Vmin时，关闭闽S7,相应的沉陷区当日用 水量为亂，在下一个洪峰到来之前科学调度调节库容，W备拦蓄洪峰，相应的采煤沉陷区蓄 滞洪效果计算见附表2。
[0081 ]附表2实施例一采煤沉陷区蓄滞洪效果计算表 [0082]
[0084]
[0085] 注：采煤沉陷区巧计洪水位对应的库容Vmax = 2432万，所述采煤沉陷区死库容 Vmin = 317万m3,与所述采煤沉陷区相连的河道讯限水位对应的洪量Qf = 547万m3,供水区当 日最大取用水重Qmu = 20万m3。
[0086] 图2为本发明采煤沉陷区蓄滞洪效果定量计算实施例蓄洪削峰效果图，如图2所 示，可W直观看出实施例采煤沉陷区对洪水的调蓄效果。
[0087] 图3-图14为本发明采煤沉陷区蓄滞洪效果定量计算方法实施例通过ENVI遥感技 术，从面积、深度、库容对某采煤沉陷区蓄滞空间历史演变规律进行研究得到分别在2015、 2020、2030年巧巾不同沉陷情景下，经过该采煤沉陷区分别拦蓄5年一遇、10年一遇、20年一 遇、50年一遇的洪水相应的调洪过程图。
[0088] 如图3-图14所示，实施例某采煤沉陷区在2015、2020、2030S种不同蓄滞空间下对 5、10、20、50年一遇的洪水均有一定的蓄洪削峰作用。由附表3可^看出实施例二某采煤沉 陷区在=种不同沉陷情景下对不同频率洪水的具体削减效果。
[0089] 附表3实施例二采煤沉陷区不同沉陷情景下对不同频率洪水削效果表
[0non1
[0091] 图15为本发明采煤沉陷区蓄滞洪效果定量计算方法实施例调蓄前后蓄滞洪量对 比示意图，如图15所示，在不同沉陷情景下经实施例某采煤沉陷区调蓄后，各级别洪量均相 应减少，调蓄前后洪量对比见附表4。
[0092] 附表4实施例二某采煤沉陷区调蓄前后蓄洪量对比表单位:亿m3 [mwl
【主权项】
1. 一种采煤沉陷区蓄滞洪效果的定量计算方法，其特征在于： 获取采煤沉陷区相连的河道天然来水日平均流量:所述河道天然来水日平均流量为24 小时内通过河流某一断面的水量的算术平均值； 获取采煤沉陷区设计洪水位对应的库容:所述设计洪水位对应的库容为当遇到设计标 准洪水时，采煤沉陷区经调洪后达到的最高水位以下的蓄水容积； 获取采煤沉陷区死库容:所述死库容为采煤沉陷区在正常运用情况下，允许消落到的 最低水位以下的蓄水容积； 获取采煤沉陷区相连的河道汛限水位对应的洪量:所述汛限水位对应的洪量为河流按 照以往汛期的最大水流量和河岸的最大承受水量设定的水位对应流量与时间的乘积，时间 按24小时计算； 根据所述的河道天然来水日平均流量、设计洪水位对应的库容、死库容、河道汛限水位 对应的洪量通过水量平衡关系定量计算采煤沉陷区蓄滞洪效果。2. 根据权利要求1所述的一种采煤沉陷区蓄滞洪效果的定量计算方法，其特征在于:所 述的采煤沉陷区还包括:河道区、沉陷区、供水区； 所述的采煤沉陷区蓄滞洪效果包括:沉陷区蓄水效果、沉陷区理论下泄量、供水区取用 水量和沉陷区库容调节变化量。3. 根据权利要求2所述的一种采煤沉陷区蓄滞洪效果的定量计算方法，其特征在于:所 述的水量平衡关系为： Vt-Vy= Qr-Qu-Qre 其中：Vt为当日沉陷区蓄水量，Vy为上一日沉陷区蓄水量，Qr为河道日天然来水量，Q uS 沉陷区当日用水量，Qre为河道下泄量。4. 根据权利要求3所述的一种采煤沉陷区蓄滞洪效果的定量计算方法，其特征在于： Qr = qr · At 其中qr为河道天然来水日平均流量，△ t为时间步长，按24小时计算； 所述的沉陷区蓄水效果为： 当QrSQf时，所述沉陷区不起蓄所述河道区来水，沉陷区理论下泄量不为零； 当Qr>Qf时，所述沉陷区起蓄所述河道区来水，沉陷区理论下泄量为零；所沭的沉陷冈理论下泄量为：所述的供水区取用水量为：其中：qf为河道区防洪汛限水位对应天然来水日平均流量，Qf为河道区防洪汛限水位对 应的日天然来水量，Vmax为沉陷区最大库容量，Vmin为沉陷区死库容量，Q。为沉陷区下泄量， Qmu为供水区当日最大取用水量，Q为供水区理论供水量，V为沉陷区即时库容。
【文档编号】G06F19/00GK106021947SQ201610370923
【公开日】2016年10月12日
【申请日】2016年5月25日
【发明人】朱梅, 程铁杰, 郑佳重, 李晓乐, 李雪凌
【申请人】安徽农业大学