主权项 |
一种街区污染物浓度分布模型建立方法,其特征在于,包括如下步骤: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>×</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>×</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>×</mo><msub><mi>V</mi><mi>y</mi></msub><mo>+</mo><msub><mi>a</mi><mn>4</mn></msub><mo>×</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>∫</mo><mn>0</mn><mo>∞</mo></msubsup><mfrac><mi>QdT</mi><mrow><msup><mrow><mo>(</mo><mn>2</mn><mi>π</mi><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup><msub><mi>σ</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><msub><mi>σ</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><msub><mi>σ</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mo>×</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>σ</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>σ</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>σ</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>πV</mi><msub><mi>σ</mi><mi>y</mi></msub><msub><mi>σ</mi><mi>z</mi></msub></mrow></mfrac><mo>×</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>σ</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>σ</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>是两个方向上的扩散系数。 |