发明名称 一种冰封期水功能区纳污能力计算方法
摘要 本发明的目的在于提供一种冰封期水功能区纳污能力计算方法,其特征在于包括构建水功能区纳污能力一维水质模型等步骤,目的在于提出了一种适合寒区河流冰封期水功能区纳污能力计算的方法,对于针对寒区冰封期水环境风险进行分析和调控具有重要意义。
申请公布号 CN104615871A 申请公布日期 2015.05.13
申请号 CN201510036674.1 申请日期 2015.01.26
申请人 中国水利水电科学研究院 发明人 胡鹏;周祖昊;贾仰文;周娜;朱乾德;贺华翔;李海红;刘盈斐;王鹏
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 四川君士达律师事务所 51216 代理人 芶忠义
主权项 一种冰封期水功能区纳污能力计算方法,其特征在于包括以下步骤:步骤一、构建水功能区纳污能力一维水质模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>C</mi><mi>x</mi></msub><mo>=</mo><msub><mi>C</mi><mn>0</mn></msub><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>KL</mi><mo>/</mo><mi>u</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mi>p</mi><mi>Q</mi></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>Kx</mi><mo>/</mo><mi>u</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>-</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000661373600000011.GIF" wi="1286" he="131" /></maths>C<sub>x</sub>——下断面污染物浓度,mg/L;C<sub>0</sub>——初始断面的污染物浓度,单位mg/L;L——计算河段长度,单位m;u——设计流量下河道断面的平均流速,单位m/s;K——污染物综合衰减系数,单位1/s;p——污染物入河速率,单位g/s;Q——水功能区设计流量,单位m<sup>3</sup>/s;x——河段排污口距下断面的纵向距离,单位m;步骤二、计算河段最大纳污能力:令C<sub>x</sub>=C<sub>s</sub>,解得计算河段最大纳污能力为:M=[C<sub>s</sub>‑C<sub>0</sub>exp(‑KL/u)]·Q·exp(Kx/u)    (5‑3)式中,C<sub>s</sub>——下断面污染物目标浓度,M——水功能区纳污能力,单位g/s;步骤三、污染物综合衰减系数K的计算:不同水温条件下K值估算关系式如下:K<sub>t</sub>=K<sub>20</sub>·θ<sup>(T‑20)</sup>    (5‑4)式中:K<sub>t</sub>——T℃时的K值,单位1/d;T——水温,单位℃,冰封期取值为0℃;K<sub>20</sub>——20℃时的K值,1/d;θ——温度系数,工业废水一般在1.03‑1.1的范围内;步骤四、设计流量下河道断面的平均流速u:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><msup><mi>R</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><msup><mi>i</mi><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>-</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000661373600000012.GIF" wi="953" he="131" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mi>A</mi><mrow><mi>b</mi><mo>+</mo><mn>2</mn><mi>h</mi><msqrt><mn>1</mn><mo>+</mo><msup><mi>m</mi><mn>2</mn></msup></msqrt></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>-</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000661373600000013.GIF" wi="1141" he="146" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>i</mi><mo>=</mo><mfrac><mi>&Delta;z</mi><mi>L</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>-</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000661373600000014.GIF" wi="1031" he="127" /></maths>其中,A为过流断面面积,b为过流断面概化为等腰梯形后的底宽,m为概化断面的边角余切值,h为过流断面水深,R为水力半径,以上各项参数均指能代表整个水功能区的平均值;n为曼宁系数;i为计算水功能区的平均坡降,Δz为水功能区上下断面高程差,根据各水功能区上、下断面所在子流域的高程数据作差求得;步骤五,水功能区设计流量Q的计算方法:根据步骤四中所得出的u值,得出,Q=uA=u·(b+h·m)·h    (5‑5)。
地址 100038 北京市海淀区复兴路甲一号