绿地生态功能导向下的CO2扩散效应研究——以沈北新区为例

绿地生态功能导向下的CO2扩散效应研究*

——以沈北新区为例

时泳，付士磊，屠兵

摘要:目前国内对绿地的研究大多限制在“量”上，还没有从空间格局上来进行研究的。二氧化碳浓度与植物的固碳能力之间有着紧密的联系，本文以沈北新区为研究对象通过运用流体力学模拟软件对二氧化碳扩散效应进行研究进而来确定绿地在空间上的布局。通过宏观、中观、微观三个角度采用实测-模拟-验证的方法对二氧化碳在空间中的分布进行分析研究。结果表明：二氧化碳主要集中在主导风向的下风向；开敞空间通风良好的地区比较利于二氧化碳的扩散；上风向比下风向扩散效果好；较高处比较低处扩散效果好；有通风廊道且与主导风向一致时（廊道与城市风向之间的夹角在30°-60°的范围内）扩散效果更好；二氧化碳会在建筑背风向形成集聚，在较高建筑上风向脚下形成湍流风降低该空间的二氧化碳浓度。

关键词:生态功能,二氧化碳,扩散效应,绿地

气候变化已然成为人类关注的重要环境问题之一，并成为影响人类生存及发展的根本问题[1]。二氧化碳含量的上升成为人类共同面临的重大挑战，要解决目前人类生活和自然的矛盾就需要构建基于目前情况下的绿地生态优化体系来满足人类生活对自然生态系统的需求[2]。然而有限的城市面积不可能无限的增加绿地的供给，我们不能一味地通过增加绿地面积来增加碳汇量，因此如何在有限绿地面积的情况下，通过大气动态模拟技术获取城市二氧化碳的空间分布格局，揭示二氧化碳空间分布特征，界定城市植物固碳效应最佳边界，优化城市绿地系统的空间布局，从空间格局的变化上增加碳汇量，从而在城市有限绿地面积下发挥植被最大生态效应就显得尤为重要。

二氧化碳浓度与植物的固碳能力之间有着紧密的联系，因此本文以沈北新区为例开展基于植物固碳释氧效应的城市绿地系统优化布局[3]，采用生态规划中的定量分析与空间技术分析相结合的研究方法，以实地外业调查数据为基础，以模型模拟为手段，以绿地生态功能的静态和动态空间分析为主要内容，构建基于定量分析的城市绿地适宜空间布局。

1 材料与方法

1.1 区域概况分析

沈北新区地处沈阳四环北郊，位于沈阳、大连、长春、哈尔滨―东北城市走廊‖中部，南靠三环，北隔辽河、万泉河与铁岭、法库县相望，东与抚顺市、铁岭县毗邻，西接辽西走

*课题“基于二氧化碳扩散模拟的城市绿地系统空间规划研究（编号：51308350）”

廊，与新民市、于洪区相连。是连接吉林、黑龙江和内蒙古三省区的黄金通道和"东北城市走廊"的枢纽重地。坐标介于东经123°16'至48'，北纬41°54'至42°11'之间。沈北新区总面积819平方千米，总人口319380人。

沈阳市沈北新区现有的农业生态绿地大约占规划区总面积的68% ，这些绿地构成了沈北新区绿地的主要基质；另外山体绿地、景区的景观绿地水域河道周围的自然景观绿地等部分作为规划区内的绿地斑块，而河流水域及其沿岸周围的绿地一起构成了沈北新区整个绿地系统的生态廊道[4]。沈北新区现状绿地布局存在的问题：（1）农田面积过大，碳汇能力低;（2）冬季西北部风沙侵袭;（3）绿地系统布局不合理，各功能类型绿地发展不均衡;（4）建成区的绿地植被类型单一，生态系统的同质化现象比较严重，生态稳定性较差。

1.2 数据来源

本文研究是基于植物的固碳释氧功能，根据植物的一般生理活性，在沈阳地区绿地植被在秋季和冬季叶片枯萎不具有固碳功能[5]，因此本文选取沈阳市春季和夏季的气候参数来对二氧化碳的扩散进行研究。二氧化碳浓度的测量采用手持二氧化碳测试仪，从2015年1月的冬天开始，到10月的秋天截止，选择早（9:00）、中（12:00）、晚（18:00）三个时段进行监测获得。同期监测点的气温、湿度及风速等参数通过中国气象网进行获取。将所获数据的平均值作为Fluent模拟初始条件。Gambit物理模型的高度确立通过阴影长度法来求得[6]。1.3 研究方法

本文采用实测-模拟-验证的方法，根据实测获取的数据运用Fluent模拟软件从宏观的域尺度、中观的建成区尺度以及微观的小区尺度对沈北新区内部的二氧化碳进行动态模拟，并通过实地采样监测的方法来验证模型精度，从而对沈北新区内二氧化碳的扩散效应进行对比分析。

1.3.1 三维模型的建立

使用Gambit进行三维模型的建立，为保证来流充分扩散，分别将宏观尺度、中观尺度、微观尺度的计算域设置为50000m×50000m×1500m、10000m×10000m×500m、1000m×1000m×100m的箱体空间。

1.3.2 网格划分与边界条件的设定

本研究运用模糊理论对模型进行简化处理，并采用结构化网格对Gambit模型进行划分。为了提高模拟的精确度，保证每个计算面都能够正确的划分网格，本文进行多次尝试根据不同的模块尺度大小，形成5m、10m、15m、20m不等的分析网格系统[7]。

在数值模拟过程中，不考虑温度对二氧化碳气体流场的影响，也不考虑空气受温度影

响的上升浮力，城市流场近地面气流速度数量级不大（20m/s 以下），因此可以将模型内部的流场计算按照不可压缩流动来处理[8]。研究区春季的主导风向为西南风，夏季主导风向为南风，模拟中风速设置为沈阳的常年平均风速3.0m/s 。设定流场的岀流面上的流体流动已得到了充分扩散，流动形态已经恢复为无建筑影响时的流动特征，因此将出口边界的压力设置为0，建筑物表面设为有摩擦的平滑墙壁[9]，重力加速度值定为9.81 m/s 2。

1.3.3 控制方程的选择

因二氧化碳密度高于室内的空气密度，将Gambit 建成的模型导入Fluent 模拟软件中设定初始模拟条件后，我们选择以下四组控制方程来进行5000次的迭代运算得出最终的模拟结果[10] 。

（1）连续性方程：

0i i U x ∂=∂ 1,2i =

（2）动量守恒方程：

()21i i j i j j i j U U p U u u v U t x x x ρ∂∂∂∂+=--+∇∂∂∂∂ 1,2j =

式中 ,,,,,,i j i j i j i j p v u u i j u u x x i j U U i j ρ—流体时均压强

—流体密度

—流体运动黏性系数

—流体脉动速度在方向上的分量

—流体脉动切应力

—方向上的坐标

—流体时均速度在方向上的分量

（3）标准的k ε-方程组：

12j t i i j t j j

k j j i j j t i i j j j j j i j U v U U k k k U v t x x x x x x U v U U U C v C t x x x k x x x εεεεσεεεεεσ⎛⎫⎛⎫∂∂∂∂∂∂∂+=++- ⎪ ⎪ ⎪ ⎪∂∂∂∂∂∂∂⎝⎭⎝⎭⎡⎤⎛⎫⎛⎫∂∂∂∂∂∂∂+=++-⎢⎥ ⎪ ⎪ ⎪ ⎪∂∂∂∂∂∂∂⎢⎥⎝⎭⎝⎭⎣

⎦ 式中 ε—湍动耗散率

（4）组分方程：