绿地生态功能导向下的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 jk 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）组分方程：()()()()s s s s s c div uc div D grad c s tρρρ∂+=+∂式中 s s s c s c s D s ρ—组分的体积浓度—组分的质量浓度—组分的扩散系数1.3.4 模拟结果验证本文选取典型的定位测量和十字样线两种采样布点方法来对模拟结果进行监测验证，共选取了8组监测点，24个测量点进行测量，监测点位布置如图1所示：图1 监测点在研究区域内的空间分布图将监测点位测量所得的二氧化碳浓度与对应点位模拟结果所得浓度两组数据运用spss 软件进行样本T 检验可得Sig 值为0.390＞0.05，表明模拟结果与实测二氧化碳数值之间差异不显著，满足模拟的精度，证明CFD 模型对二氧化碳扩散效应场的模拟是有效可行的[11]。

2 结果与分析2.1 春夏季二氧化碳模拟结果分析图2 沈北新区春、夏季二氧化碳扩散模拟结果宏观域尺度的模拟主要体现了大范围内的二氧化碳走势，通过对沈北新区春、夏季二氧化碳的模拟结果分析我们可以总结出以下几个方面的结论：第一，二氧化碳的浓度较大的区域主要分布在沈北新区的下风向；第二，通过不同主导风向下的二氧化碳浓度分布可以看出，在有通风廊道的地区二氧化碳比较容易扩散，相应的二氧化碳浓度较低。

第三，工业用地地块较长边垂直于主导风向时，其下风向上的二氧化碳扩散范围较其较短边垂直于主导风向时的扩散范围广。

图3 新城子新城春、夏季二氧化碳扩散模拟结果中观尺度的模拟主要选取了周围环境对其影响较小的新城子新城作为模拟对象，通过对模拟结果的分析我们可以总结出以下几点结论：第一，新城子新城通风状况较为良好，二氧化碳的扩散较为适宜，其下风向上的二氧化碳没有形成集聚。

第二，在滨河的通风廊道处通风效果好，二氧化碳的浓度较其他地区更低。

第三，对于没有障碍物的西部工业用地部分在通风良好的春季其所产生的二氧化碳可以得到及时扩散。

在夏季，主导风向为南风，其下风向上的二氧化碳扩散距离较远，浓度稀释较慢扩散效果不及春季主导风向为西南风时的好，这主要说明了两个问题，其一，说明了通风廊道对于城区二氧化碳扩散的重要性。

其二，说明了在城市规划中道路、水系、绿道等形成的通风廊道对城市二氧化碳扩散的影响；图4 天泰小区春、夏季二氧化碳扩散模拟结果小区尺度的模拟比较精确直观，通过对小区内春、夏季二氧化碳的模拟结果的对比分析，可总结出以下几方面的结论：第一，二氧化碳占空气体积比较大区域为整个小区的下风向；第二，二氧化碳浓度最高的区域一是该区域缺乏通风道，上风向建筑体块比较长，或者上风向前一排的建筑挡住了当前建筑的风道。

二是该区域的下风向存在高度较高的建筑从而形成了湍流风，湍流风带走了较高建筑脚下的二氧化碳，然而却阻挡了上风向建筑二氧化碳的扩散从而使上风向二氧化碳形成集聚，提高了该区域的二氧化碳浓度；2.2 典型水平剖面二氧化碳分布规律高度H=0.1米时为地面高度；H=1.5米是人类呼吸的敏感范围；H=9米是地被植物及灌木植物生长的极限点；H=18米是一般植物生长的极值，文章选取这四个典型高度的模拟结果来进行分析，并进行相应的植物配置。

2.2.1 典型水平剖面模拟结果通过模拟分析可得宏观、中观、微观尺度的二氧化碳空间分布结果，以下为0.1m，1.5m，9m，18m不同高度水平截面的颗粒物浓度分布，可见不同平面的大气污染物的分布均匀程度、浓度降低趋势各不相同。

宏观模拟结果：H = 0.1米H = 1.5米H = 9米H =18米图5 沈北新区春季二氧化碳扩散平剖面分布结果H = 0.1米H = 1.5米H = 9米H =18米图6 沈北新区夏季二氧化碳扩散平剖面分布结果中观模拟结果：H = 0.1米H = 1.5米H = 9米H =18米图7 新城子新城春季二氧化碳平剖面分布结果H = 0.1米H = 1.5米H = 9米H =18米图8 新城子新城夏季二氧化碳平剖面分布结果微观模拟结果：H = 0.1米H = 1.5米H = 9米H =18米图9 小区春季二氧化碳平剖面分布结果H = 0.1米H = 1.5米H = 9米H =18米图10 小区夏季二氧化碳平剖面分布结果2.2.2 不同截面二氧化碳模拟结果分析表1 宏观水平截面模拟结果分析表2 中观水平截面模拟结果分析表3 微观水平截面模拟结果分析2.3 典型垂直剖面二氧化碳分布规律2.3.1 垂直剖面的选取为进一步揭示和验证二氧化碳在不同空间形态的分布特征，选取小区春夏二季小区同一位置的垂直剖面对以下12个不同空间形态点位的二氧化碳空间分布来进行分析比较，剖面位置见图11 ：图11 轴剖面位置示意图春夏两季模拟结果的垂直剖面二氧化碳扩散效应见图12，通过城市冠层内轴剖面的二氧化碳分布云图可以看出，小区内释放的二氧化碳在风环境的影响下空间格局扩散分布与静态分布有所不同。