发明名称 一种街区污染物浓度分布模型建立方法
摘要 本发明提供一种街区污染物浓度分布模型建立方法,利用CFD模型获取风场矩阵,以此为基础建立街区内风场多元线性回归方程,实现街区内风场的计算,基于高斯烟团模型和GPU并行技术获取街区内污染物浓度分布,基于高斯烟羽模型和背景风获取街区外污染物浓度分布。该方法简化并加快了风场的计算过程,提高了街区内外部污染物浓度的计算效率。
申请公布号 CN104008229A 申请公布日期 2014.08.27
申请号 CN201410181396.4 申请日期 2014.04.30
申请人 北京大学 发明人 刘峻峰;黄欣;陶玮;孟靖;陈卓
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京三聚阳光知识产权代理有限公司 11250 代理人 寇海侠
主权项 一种街区污染物浓度分布模型建立方法,其特征在于,包括如下步骤:S1:建立多元线性回归方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><msub><mi>a</mi><mn>0</mn></msub><mo>+</mo><msub><mi>a</mi><mn>1</mn></msub><mo>&times;</mo><mfrac><msub><mi>H</mi><mn>1</mn></msub><msub><mi>H</mi><mn>2</mn></msub></mfrac><mo>+</mo><msub><mi>a</mi><mn>2</mn></msub><mo>&times;</mo><mfrac><msub><mi>H</mi><mn>3</mn></msub><msub><mi>H</mi><mn>4</mn></msub></mfrac><mo>+</mo><msub><mi>a</mi><mn>3</mn></msub><mo>&times;</mo><msub><mi>V</mi><mi>y</mi></msub><mo>+</mo><msub><mi>a</mi><mn>4</mn></msub><mo>&times;</mo><msub><mi>V</mi><mi>x</mi></msub></mrow>]]></math><img file="FDA0000499155470000011.GIF" wi="999" he="160" /></maths>以上各参数表示的含义为:u是风速;H<sub>1</sub>、H<sub>2</sub>、H<sub>3</sub>、H<sub>4</sub>为待分析街区的高度或者宽度或者长度,且H<sub>1</sub>≠H<sub>2</sub>,H<sub>3</sub>≠H<sub>4</sub>;V<sub>x</sub>和V<sub>y</sub>是背景风平行于街区方向的分量或者是背景风垂直于街区方向的分量,且V<sub>x</sub>≠V<sub>y</sub>;a<sub>0</sub>‑a<sub>4</sub>是待定系数,基于CFD模型获得的风场矩阵中的风速进行多元回归计算得到;S2:将H<sub>1</sub>、H<sub>2</sub>、H<sub>3</sub>、H<sub>4</sub>,V<sub>x</sub>和V<sub>y</sub>代入到上述多元线性回归方程中后得到风速u;S3:基于高斯烟团模型获得街区内污染物浓度分布C<sub>1</sub>(x,y,z,T):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>C</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mo>&infin;</mo></msubsup><mfrac><mi>QdT</mi><mrow><msup><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><msub><mi>&sigma;</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><msub><mi>&sigma;</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><msub><mi>&sigma;</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mo>&times;</mo><msup><mi>e</mi><mrow><mo>[</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>uT</mi><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><msup><mi>y</mi><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>y</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>z</mi><mo>-</mo><mi>H</mi><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac></mrow></msup></mrow>]]></math><img file="FDA0000499155470000012.GIF" wi="1563" he="236" /></maths>以上各参数表示的含义为:x,y,z是坐标值,T是时间,Q是点源排放速率,σ<sub>x</sub>,σ<sub>y</sub>,σ<sub>z</sub>是三个方向上的扩散系数,u是步骤S2中得到的风速,H是点源高度;S4:基于高斯烟羽模型获得街区外污染物浓度分布C<sub>2</sub>(x,y,z):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>C</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>Q</mi><mrow><mn>2</mn><mi>&pi;V</mi><msub><mi>&sigma;</mi><mi>y</mi></msub><msub><mi>&sigma;</mi><mi>z</mi></msub></mrow></mfrac><mo>&times;</mo><msup><mi>e</mi><mrow><mo>[</mo><mo>-</mo><mfrac><msup><mi>y</mi><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>y</mi><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mi>z</mi><mo>-</mo><mi>H</mi><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>z</mi><mn>2</mn></msubsup></mrow></mfrac><mo>]</mo></mrow></msup></mrow>]]></math><img file="FDA0000499155470000013.GIF" wi="1069" he="262" /></maths>以上各参数表示的含义为:V是背景风风速,Q是点源排放速率,H为点源高度,σ<sub>y</sub>和σ<sub>z</sub>是两个方向上的扩散系数。
地址 100871 北京市海淀区颐和园路5号