期刊文献+
共找到18篇文章
< 1 >
每页显示 20 50 100
一、二维耦合模型在茅洲河水环境治理中的应用 被引量:9
1
作者 郭聪 施家月 +1 位作者 张亚力 郁关明 《环境影响评价》 2019年第4期59-62,共4页
为了给茅洲河水环境综合整治工程方案优化提供依据,定量评估工程实施对水质改善效果,采用DHI MIKE Flood耦合MIKE 11与MIKE 21,建立了茅洲河一维模型与口外海域二维模型相耦合的水动力与水质数学模型。利用水文与水质同步监测资料对模... 为了给茅洲河水环境综合整治工程方案优化提供依据,定量评估工程实施对水质改善效果,采用DHI MIKE Flood耦合MIKE 11与MIKE 21,建立了茅洲河一维模型与口外海域二维模型相耦合的水动力与水质数学模型。利用水文与水质同步监测资料对模型进行了率定验证,验证结果表明本模型较好地反映了茅洲河水动力与氨氮浓度变化过程,可用于茅洲河水环境治理工程的水质模拟预测。结合水环境治理工程各阶段入河污染负荷,利用所建模型对茅洲河及主要一级、二级支流的氨氮浓度进行了模拟。结果表明,在水环境综合治理工程实施后,茅洲河氨氮浓度明显降低,水质可得到显著改善,至2017年水体黑臭时间明显下降,至2020年可全面消除黑臭,并可达到V类水。 展开更多
关键词 一、二维耦合模型 水环境治理 茅洲河
在线阅读 下载PDF
一、二维耦合模型在梯级水电站大坝溃坝分析中的应用 被引量:2
2
作者 张凤山 杜德进 杨彦龙 《大坝与安全》 2023年第1期7-12,共6页
合理论证溃坝洪水灾害影响对大坝风险评价和应急管理具有重要意义。针对梯级水电站大坝连溃风险,建立一、二维耦合模型,对溃坝洪水演进规律及大坝下游平原洪泛区淹没影响开展分析。以W江梯级水电站为例,考虑汛期与非汛期水库不同运行情... 合理论证溃坝洪水灾害影响对大坝风险评价和应急管理具有重要意义。针对梯级水电站大坝连溃风险,建立一、二维耦合模型,对溃坝洪水演进规律及大坝下游平原洪泛区淹没影响开展分析。以W江梯级水电站为例,考虑汛期与非汛期水库不同运行情况下两种工况的溃坝影响,从典型断面洪峰流量、最高水位以及峰现时间等角度分析溃坝洪水演进规律。在此基础上,以洪泛区最大淹没水深为代表性指标,对溃坝洪水风险等级进行了划分,并绘制洪水风险图,工况1、工况2高风险区面积分别达到431.31 km^(2)和575.76 km^(2),表明溃坝洪水对下游地区影响严重。相关成果可为应急预案制定和防灾减灾提供重要参考。 展开更多
关键词 溃坝洪水 梯级大坝 一、二维耦合模型
在线阅读 下载PDF
一、二维耦合数学模型在感潮河网洪水风险图编制中的应用 被引量:8
3
作者 杨莉玲 宋利祥 +3 位作者 邓军涛 徐爽 胡晓张 孙倩雯 《长江科学院院报》 CSCD 北大核心 2017年第9期36-40,共5页
为了解决感潮河网径流、潮汐交汇,动力复杂,溃决洪水难以用经验公式准确概化的问题,建立直接以溃口为耦合断面的河网一维、保护区二维侧向耦合模型,将感潮河网与保护区一体化,避免环境因素及经验参数的不确定性带来的溃口流量估算误差... 为了解决感潮河网径流、潮汐交汇,动力复杂,溃决洪水难以用经验公式准确概化的问题,建立直接以溃口为耦合断面的河网一维、保护区二维侧向耦合模型,将感潮河网与保护区一体化,避免环境因素及经验参数的不确定性带来的溃口流量估算误差。典型算例和中顺大围溃决洪水情景模拟表明:洪水自溃口集中喷射出后分散流向围内,流态受围内下垫面影响显著,溃口水位、流量随外江潮位涨落而起伏变化,由溃口流量过程线计算所得溃口水量与根据淹没区各单元的面积和水深计算的围内总水量一致。模拟成果直接反映下垫面、水头差、潮位涨落,及溃口流态对溃决洪水的综合影响,反映出本耦合模型计算溃口流量及对溃决洪水模拟的合理性、可信性,具有较好的应用前景。 展开更多
关键词 溃决洪水 一、二耦合数学模型 感潮河网 洪水风险图 侧向联解
在线阅读 下载PDF
考虑复杂河工建筑物的漫堤洪水一、二维水动力耦合模型 被引量:6
4
作者 张慧萍 田福昌 +1 位作者 耿庆柱 王建家 《水电能源科学》 北大核心 2018年第9期81-85,共5页
针对景观坝和过水路面等河工建筑物严重影响河道排洪泄洪问题,基于MIKE FLOOD平台构建漫堤洪水一、二维水动力耦合模型,模拟分析复杂河工建筑物壅水与漫溢洪水淹没演进过程。采用Abbott六点隐式格式和显式有限体积法分别离散求解一、二... 针对景观坝和过水路面等河工建筑物严重影响河道排洪泄洪问题,基于MIKE FLOOD平台构建漫堤洪水一、二维水动力耦合模型,模拟分析复杂河工建筑物壅水与漫溢洪水淹没演进过程。采用Abbott六点隐式格式和显式有限体积法分别离散求解一、二维水流运动方程,利用干湿水深理论优化模型运算效率,将河道内过水路面与景观坝作为宽顶堰概化处理,同步计算桥涵淹没出流。将模型应用于宁夏茹河(店洼水库至小河河口)河道跨河工程阻水效果及漫堤洪水的数值模拟。结果表明,模拟结果较真实地反映了河工建筑物对水流的阻碍效应及桥涵动态过流效果。研究成果可为宁夏茹河及其他类似河道的整治工程规划设计及风险评价提供参考。 展开更多
关键词 一、二水动力耦合模型 河工建筑物 干湿水深 漫堤洪水 有限体积法
在线阅读 下载PDF
防洪保护区水动力一二维精细化模拟模型及应用 被引量:3
5
作者 张庆梓 刘小龙 +1 位作者 陈俊鸿 彭思韦 《治淮》 2016年第9期19-21,共3页
为准确合理地进行洪水演进模拟,完成防洪保护区的洪水风险图编制工作,本文建立了能够对防洪保护区溃堤洪水进行模拟的一、二维耦合水动力学模型。一维水动力学模型采用有限差分法求解,二维水动力学模型采用二维有限体积法求解浅水流方... 为准确合理地进行洪水演进模拟,完成防洪保护区的洪水风险图编制工作,本文建立了能够对防洪保护区溃堤洪水进行模拟的一、二维耦合水动力学模型。一维水动力学模型采用有限差分法求解,二维水动力学模型采用二维有限体积法求解浅水流方程组。模拟结果表明,与传统方法相比该模型能有效提高复杂区域内模型计算效率,所得信息可为防洪决策部门提供参考依据。 展开更多
关键词 水动力学 一、二维耦合模型 洪水风险分析 药湖联圩
在线阅读 下载PDF
MIKE耦合模型在内河航道水面线计算中的应用
6
作者 张会杰 周昂 《河南科技》 2024年第17期37-40,共4页
【目的】在航道水面线计算过程中,常常会有水闸等水工建筑物对水产生不可忽略的阻碍作用。为研究含有水工建筑物的耦合模型在航道水面线计算中的适用性及准确性,以汝河为例,建立数学模型对河道典型工况进行研究。【方法】通过实测河道... 【目的】在航道水面线计算过程中,常常会有水闸等水工建筑物对水产生不可忽略的阻碍作用。为研究含有水工建筑物的耦合模型在航道水面线计算中的适用性及准确性,以汝河为例,建立数学模型对河道典型工况进行研究。【方法】通过实测河道横断面,创新性地建立闸位处地形的一、二维河道耦合模型,并在模型率定后进行水面线计算分析。【结果】分析闸位处的流场分布情况可知,模拟结果能够体现水闸结构对过闸水流的阻碍作用,从而影响闸前闸后的水位变化。【结论】一、二维河道耦合模型结果率定和适用性较好,该方法的应用为含有水工建筑物的水面线计算提供了参考。 展开更多
关键词 水面线 水工建筑物 一、二维耦合模型
在线阅读 下载PDF
陆面-河网耦合模型在极端风暴潮淹没危险性分析中的应用 被引量:1
7
作者 滕飞 李路 +3 位作者 贝竹园 王军 邱桔斐 章晓洁 《华东师范大学学报(自然科学版)》 CAS CSCD 北大核心 2021年第3期138-150,共13页
在风暴潮灾害中,台风、天文高潮位、区间暴雨和上游洪水等因子存在同时出现的可能,对沿海地区的防汛安全形成了较为严重的威胁.双因子、多因子影响下的复合风暴潮灾害成为沿海地区灾害风险评估的主要研究对象.风暴潮水通过漫堤和溃口进... 在风暴潮灾害中,台风、天文高潮位、区间暴雨和上游洪水等因子存在同时出现的可能,对沿海地区的防汛安全形成了较为严重的威胁.双因子、多因子影响下的复合风暴潮灾害成为沿海地区灾害风险评估的主要研究对象.风暴潮水通过漫堤和溃口进入陆域后,河网对潮水的调蓄和运输作用直接影响对受灾区域灾害危险性的评估结果.本研究选取上海市金山区作为典型区域,基于MIKE11和MIKE21模型分别建立金山区一维河网模型和二维陆面模型,并进行耦合计算,模拟分析在风暴潮、台风、区间暴雨和流域洪水的综合影响下,金山区地面积水变化过程,为区域的灾害危险性分析提供依据.模型计算结果表明,在考虑陆域河网调蓄作用后,风暴潮引起的金山区地面积水大幅度减小,区域内整体淹没分布出现一定程度的改变.综合考虑风暴潮及区间暴雨和流域洪水与仅考虑风暴潮条件相比,金山区中部和北部大部分区域危险性等级降低,西北角危险性等级升高. 展开更多
关键词 风暴潮 河网 一、二模型耦合
在线阅读 下载PDF
溃堤洪水的一维、二维耦合水动力模型及应用 被引量:38
8
作者 张大伟 李丹勋 +1 位作者 陈稚聪 王兴奎 《水力发电学报》 EI CSCD 北大核心 2010年第2期149-154,共6页
为了减小因溃堤洪水而造成的人员生命和财产损失,本文建立了能够对天然河道溃堤水流进行模拟的一维、二维耦合数学模型。一维水流模型采用Preissmann格式离散求解,二维水流模型采用基于非结构网格的Roe格式离散求解,在一、二维模型的链... 为了减小因溃堤洪水而造成的人员生命和财产损失,本文建立了能够对天然河道溃堤水流进行模拟的一维、二维耦合数学模型。一维水流模型采用Preissmann格式离散求解,二维水流模型采用基于非结构网格的Roe格式离散求解,在一、二维模型的链接处采用堰流公式实现河槽内外水流的交互。通过典型测试算例验证了本耦合模型的可靠性和稳定性。检验后的耦合模型为大江大河的溃堤洪水计算提供了基础,将该模型成功应用于松花江哈尔滨段松北分洪区溃堤洪水的模拟,所得信息可以为防洪决策部门提供参考依据。 展开更多
关键词 水动力学 溃堤洪水 一、二维耦合模型 松北分洪区
原文传递
基于GIS与水动力模型融合技术的山洪沟溃堤风险分析 被引量:2
9
作者 田俊玲 田福昌 +2 位作者 李天浩 孙倩 李雪飞 《水电能源科学》 北大核心 2018年第2期64-67,120,共5页
为加强洪灾预防与应对措施建设,融合GIS和一、二维水动力数值计算方法,建立了具有实际地形的山洪沟溃堤风险分析精细化仿真模型,并将其应用于宁夏贺兰山东麓小峰沟溃堤洪水演进模拟。采用非结构化网格与内边界处理技术合理优化了复杂边... 为加强洪灾预防与应对措施建设,融合GIS和一、二维水动力数值计算方法,建立了具有实际地形的山洪沟溃堤风险分析精细化仿真模型,并将其应用于宁夏贺兰山东麓小峰沟溃堤洪水演进模拟。采用非结构化网格与内边界处理技术合理优化了复杂边界与线状地物阻挡水效果,精确模拟了小峰沟遭遇100年一遇洪水导洪堤溃决造成的淹没风险,基于GIS分析了洪水演进不同时刻的风险分布特征与影响情况,实现了GIS与水动力数值模型的有机融合。计算结果基本合理可靠,可为该地区沟道治理、风险管理与灾防决策提供参考。 展开更多
关键词 GIS 一、二维耦合模型 精细化仿真 风险分析 山洪沟 溃堤洪水
在线阅读 下载PDF
宁晋泊和大陆泽蓄滞洪区洪水淹没历时及洪水风险分析 被引量:17
10
作者 王晓磊 韩会玲 李洪晶 《水电能源科学》 北大核心 2013年第8期59-62,共4页
为加强宁晋泊和大陆泽蓄滞洪区的安全建设规划,利用MIKE FLOOD软件包建立宁晋泊和大陆泽蓄滞洪区的一、二维耦合模型,采用"96.8"历史洪水淹没资料进行模型验证,对大陆泽任县环水村不同时刻的模拟水位和实测水位进行了对比,模... 为加强宁晋泊和大陆泽蓄滞洪区的安全建设规划,利用MIKE FLOOD软件包建立宁晋泊和大陆泽蓄滞洪区的一、二维耦合模型,采用"96.8"历史洪水淹没资料进行模型验证,对大陆泽任县环水村不同时刻的模拟水位和实测水位进行了对比,模拟结果与实际基本吻合。利用该模型模拟了蓄滞洪区50年一遇的洪水淹没过程,得到了洪水淹没范围、水深、历时等洪水淹没特征信息,并以最大淹没水深单指标来进行洪水风险分区,绘制洪水风险图,为蓄滞洪区的安全建设规划提供了重要依据。 展开更多
关键词 蓄滞洪区 数值模拟 一、二维耦合模型 洪水淹没风险 MIKE FLOOD
在线阅读 下载PDF
西三洼洪水演进数值模拟及洪水风险分析 被引量:3
11
作者 李大鸣 罗珊 +2 位作者 范丽虹 张弘强 熊明明 《水利水电技术》 CSCD 北大核心 2018年第8期78-86,共9页
为模拟西三洼蓄滞洪区洪水演进,对其洪水风险进行分析,以一、二维非恒定流控制方程为基本理论,采用有限体积法,结合西三洼地形资料及大清河流域水文资料,建立一、二维耦合的洪水演进数学模型。采用63·8洪水实测资料对模型进行验证... 为模拟西三洼蓄滞洪区洪水演进,对其洪水风险进行分析,以一、二维非恒定流控制方程为基本理论,采用有限体积法,结合西三洼地形资料及大清河流域水文资料,建立一、二维耦合的洪水演进数学模型。采用63·8洪水实测资料对模型进行验证,结果基本吻合。对不同工况及重现期下的设计洪水进行洪水演进模拟计算并绘制洪水风险图。研究表明:该模型可模拟西三洼蓄滞洪区各分洪口门开扒时间、各时刻淹没面积、淹没水深及滞洪量等水情信息及其变化过程;规划情况(增加安全区)的滞洪总量大于现状情况滞洪总量;发生百年一遇洪水时,西三洼蓄滞洪区77%地区都遭受洪水侵害,其中57%区域为危险区,11%区域为重灾区,5%区域为中灾区。研究成果可为防汛部门制定科学合理的防洪减灾规划提供参考依据。 展开更多
关键词 西三洼蓄滞洪区 洪水演进 一、二维耦合模型 有限体积法
在线阅读 下载PDF
埒子口滩涂围垦对沂北地区防洪除涝的影响
12
作者 杨晨霞 曹命凯 +1 位作者 袁尧 洪继龙 《中国水运(下半月)》 2015年第11期219-221,共3页
随着沿海经济的发展,围垦已成为解决土地资源稀缺的主要方法。埒子口是沂北地区防洪排涝的主要入海口门之一,埒子口周边分布有善后河闸、烧香河闸、车轴河闸、五图闸、图西闸五座排涝闸,共有烧香支河、古泊善后河、车轴河等在此汇流入... 随着沿海经济的发展,围垦已成为解决土地资源稀缺的主要方法。埒子口是沂北地区防洪排涝的主要入海口门之一,埒子口周边分布有善后河闸、烧香河闸、车轴河闸、五图闸、图西闸五座排涝闸,共有烧香支河、古泊善后河、车轴河等在此汇流入海。应用MIKE11、MIKE21建立沂北地区及埒子口近海域的一维、二维耦合模型,利用模型对围垦前后河道潮流场进行模拟,分析围垦前后河道内水动力条件的变化以及对上游诸闸过流能力的影响。相较于围垦前,围垦临近区域的涨落潮平均流速基本均有所减小,而河道的水动力条件变化则很小,围垦对上游诸闸过流基本无影响。 展开更多
关键词 埒子口 滩涂围垦 一、二维耦合模型
在线阅读 下载PDF
塔基建设对蓄滞洪区洪水演进影响评价——以蒙西—天津南输变电工程为例 被引量:4
13
作者 李大鸣 田顺发 +2 位作者 张弘强 李彦卿 陈靖 《水利水电技术》 CSCD 北大核心 2018年第10期28-36,共9页
为研究特高压输电线路的塔基建设对蓄滞洪区洪水演进的影响,以一、二维非恒定流方程为基础,建立了一维河道嵌套于二维大清河系蓄滞区的洪水演进数学模型。应用该模型模拟了工程建设前后的洪水淹没过程,同时对塔基所在网格按照工程建设... 为研究特高压输电线路的塔基建设对蓄滞洪区洪水演进的影响,以一、二维非恒定流方程为基础,建立了一维河道嵌套于二维大清河系蓄滞区的洪水演进数学模型。应用该模型模拟了工程建设前后的洪水淹没过程,同时对塔基所在网格按照工程建设前后的12个方案进行局部加密计算。结果表明:计算结果与实测洪水资料的洪水到达时间、峰现时间相差在1%以内,表明模型可靠;工程建设前后的洪水演进过程、淹没范围、最大水深分布、流场分布基本相同;塔基建设对水位影响不大,产生的最大壅水为0. 022 m,但是对局部流速影响较大,最大减少0. 545 6 m/s,占建设前流速的56%,造成冲刷深度达1. 3 m;塔基建设对200年一遇洪水的流速减小量是100年一遇洪水的2. 5倍,圆-方柱塔基对局部水流形态的影响是圆柱式的2倍。研究成果可为评价塔基建设对蓄滞洪区的防洪影响提供参考。 展开更多
关键词 塔基建设 蓄滞洪区 洪水演进 一、二耦合数学模型 有限体积法
在线阅读 下载PDF
黄河卫宁段河南溃堤洪水风险及沟渠倒灌影响 被引量:3
14
作者 岳志春 周跃华 +1 位作者 曹鲁赣 RIZWAN Qadir 《人民黄河》 CAS 北大核心 2021年第4期67-71,共5页
黄河卫宁段河南溃堤洪水灾害风险较大,引排水沟渠倒灌危害较为严重,威胁沿岸交通、通信、水利工程等基础设施以及人民群众生命财产安全,为加强洪灾预防与应对措施建设,基于一、二维水动力学数值模拟方法及沟渠倒灌影响评估方法,采用非... 黄河卫宁段河南溃堤洪水灾害风险较大,引排水沟渠倒灌危害较为严重,威胁沿岸交通、通信、水利工程等基础设施以及人民群众生命财产安全,为加强洪灾预防与应对措施建设,基于一、二维水动力学数值模拟方法及沟渠倒灌影响评估方法,采用非结构化网格与内边界处理技术合理优化了复杂边界与线状地物阻挡水效果,建立了考虑排水沟渠倒灌影响的溃堤风险分析精细化仿真模型,并将其应用于黄河卫宁段河南保护区溃堤洪水演进模拟,并评估了引排水沟渠倒灌影响。结果表明:所建模型能较精确地模拟卫宁段河南保护区遭遇100 a一遇洪水堤防溃决造成的淹没风险,当卫宁段河南保护区遭遇100 a一遇洪水在泉眼山处溃堤时,跃进渠、七星渠和张裕沟将不发生倒灌退水现象,南河子沟将发生倒灌退水现象,南河子沟倒灌长度为3910 m。 展开更多
关键词 一、二维耦合模型 精细化仿真 倒灌 引排水沟渠 溃堤洪水 黄河卫宁段
在线阅读 下载PDF
肇庆市大旺第一水厂一级抽水泵站防洪方案比较
15
作者 黄兆玮 陈刚 +2 位作者 刘霞 梁海涛 徐辉荣 《电网与清洁能源》 2006年第S1期93-97,共5页
肇庆市大旺第一水厂一级抽水泵站工程的泵房建于大旺围外河道滩地内,为满足堤围防洪需要,提出把泵房向河道方向移出的修改方案,方案修改后阻流面积增大,但上游壅水水位的基本不变,泵站处跌水略增。方案修改后更利于堤围安全。
关键词 肇庆大旺第一水厂一级抽水泵站 防洪 一、二耦合水动力学数学模型
在线阅读 下载PDF
富水干流下游洪水风险图编制 被引量:2
16
作者 何娟 阮鹏高 +2 位作者 董苇 范杨臻 丁俊芝 《人民长江》 北大核心 2018年第S1期20-25,共6页
中小流域洪水分析资料匮乏且流域控制性水利工程较多,因此中小流域设计洪水过程推求及洪水分析方案的概化相对复杂。以富水干流下游洪水风险分析为例,通过综合分析富水下游洪灾成因及典型洪水、流域控制性水利工程调度方案、水利工程现... 中小流域洪水分析资料匮乏且流域控制性水利工程较多,因此中小流域设计洪水过程推求及洪水分析方案的概化相对复杂。以富水干流下游洪水风险分析为例,通过综合分析富水下游洪灾成因及典型洪水、流域控制性水利工程调度方案、水利工程现状防洪能力等因素确定富水下游洪水分析方案。采用水文分析计算方法推求富水下游设计洪水过程;采用一、二维耦合水动力学模型进行洪水分析计算;结合富水下游土地利用状况,进行了灾情分析和损失评估;根据淹没水深与到达时间计算结果,将风险信息添加到地图上,最终编制完成富水下游洪水风险图。 展开更多
关键词 洪水风险图 中小流域 降雨径流经验相关图法 一、二耦合水动力学模型
在线阅读 下载PDF
黄庄洼蓄滞洪区洪水演进数值模拟研究
17
作者 杨丽萍 周军 +1 位作者 周志华 刘海涛 《海河水利》 2015年第1期50-52 60,共4页
建立了潮白新河、永定新河及黄庄洼蓄滞洪区联合调度的一、二维耦合水动力学模型,模拟计算了20年一遇、50年一遇设计洪水在蓄滞洪区的演进过程,分析了蓄滞洪区淹没面积、蓄水量的变化过程。此项研究可为蓄滞洪区防洪工程建设、土地利用... 建立了潮白新河、永定新河及黄庄洼蓄滞洪区联合调度的一、二维耦合水动力学模型,模拟计算了20年一遇、50年一遇设计洪水在蓄滞洪区的演进过程,分析了蓄滞洪区淹没面积、蓄水量的变化过程。此项研究可为蓄滞洪区防洪工程建设、土地利用、分洪决策等提供有益的参考。 展开更多
关键词 蓄滞洪区 一、二水力学耦合模型 洪水演进 数值模拟
在线阅读 下载PDF
北部湾及广西近海潮流数值模拟 被引量:6
18
作者 施华斌 牛小静 余锡平 《清华大学学报(自然科学版)》 EI CAS CSCD 北大核心 2012年第6期791-797,共7页
该文采用二维浅水模型模拟分析了北部湾海域潮流系统,并在此基础上建立了广西沿海一、二维耦合水动力模型,模型中考虑了广西沿岸6条主要入海河流对近岸水动力的影响,进一步分析了广西近海海域潮致余流场及其对近海污染物输移扩散的影响... 该文采用二维浅水模型模拟分析了北部湾海域潮流系统,并在此基础上建立了广西沿海一、二维耦合水动力模型,模型中考虑了广西沿岸6条主要入海河流对近岸水动力的影响,进一步分析了广西近海海域潮致余流场及其对近海污染物输移扩散的影响。模型潮位计算值和实测值符合良好。结果显示:K1分潮在北部湾口形成一个逆时针旋转潮波系统,M2分潮在湾顶越南海防附近形成逆时针旋转潮波系统,广西近海海域K1分潮占优。广西沿海海域内,离岸较近的海域余流场相对较强,余流速度受浅水效应及湾口等地形缩窄影响而增强;离岸较远海域余流明显减弱,水体交换能力降低。 展开更多
关键词 北部湾潮流 一、二维耦合模型 潮致余流
原文传递
上一页 1 下一页 到第
使用帮助 返回顶部