張海曦 李燁 邱稼儀

（1.西安長慶科技工程有限責任公司；2.西安交通大學人居學院）

長慶油田保障點建筑室外顆粒物污染預測

張海曦1李燁1邱稼儀2

（1.西安長慶科技工程有限責任公司；2.西安交通大學人居學院）

長慶油田的前線生產(chǎn)保障點建筑大多距離油氣生產(chǎn)區(qū)域較近；因此，附著有油氣田污染性廢氣的顆粒物有可能隨空氣流動進入辦公區(qū)域，對工作人員的身體健康產(chǎn)生影響。為了保證油田保障點的室內(nèi)空氣質(zhì)量，采用計算流體動力學的方法對保障點建筑的室外顆粒物擴散和分布情況進行預測。以環(huán)江作業(yè)區(qū)保障點為例，建立了其建筑物的三維模型，根據(jù)當?shù)氐臍庀髼l件，模擬計算了以PM10為代表的污染性顆粒物對保障點建筑的污染情況。計算結(jié)果顯示，室外顆粒物對保障點建筑造成了嚴重污染。通過分析計算結(jié)果，對保障點建筑的顆粒物污染防治問題提出了若干建議，能夠改善保障點的室內(nèi)外空氣質(zhì)量和建筑風環(huán)境，降低室外顆粒物在建筑內(nèi)的濃度，保證工作人員的身體健康。對新建保障點的選址和已建的建筑室外顆粒物污染較嚴重的保障點的治理具有重要的參考價值。

室外顆粒物；計算流體動力學；空氣質(zhì)量；建筑環(huán)境

引言

長慶油田地處我國鄂爾多斯盆地，大部分一線生產(chǎn)單位氣候條件惡劣，使廣大基層員工有一個良好的工作和生活環(huán)境是油田公司穩(wěn)產(chǎn)上產(chǎn)的必要保證。我國西北地區(qū)最基本的環(huán)境特征是沙塵廣布和干旱[1]，空氣中可吸入顆粒物的含量明顯偏高?？晌腩w粒物是粒徑小于或等于10 μm的顆粒物（PM10），能夠長時間懸浮于空氣中，并伴隨呼吸作用進入人體。流行病學的研究顯示，可吸入顆粒物對人體健康有著非常顯著的負面影響，可能引起多種呼吸系統(tǒng)疾病，甚至引發(fā)癌癥[2]。另外，原油的開采、處理和集輸?shù)倪^程中會發(fā)生跑冒滴漏現(xiàn)象，產(chǎn)生一些污染性氣體，包括烴類、二氧化硫、氮氧化物和一氧化碳等[3]。這些對人體有害的污染性氣體在顆粒物的吸附作用下會附著于顆粒物表面[4]。盡管這些附著有污染性氣體的顆粒物主要產(chǎn)生于油氣生產(chǎn)區(qū)域，但長慶油田前線生產(chǎn)保障點大多距離生產(chǎn)區(qū)域較近。而相關研究表明，建筑附近的大氣顆粒物能夠隨著建筑通風進入室內(nèi)，室內(nèi)環(huán)境中來源于室外大氣的顆粒物占比可以達到75%以上[5]。因此，有必要對長慶油田前線生產(chǎn)保障點附近的大氣顆粒物擴散和分布情況進行研究。

選擇長慶油田典型的前線生產(chǎn)保障點進行研究，以長慶油田第七采油廠環(huán)江作業(yè)區(qū)前線生產(chǎn)保障點為例，該保障點周圍部署有油井并且與環(huán)江第一聯(lián)合站僅一墻之隔。環(huán)江油田第一聯(lián)合站于2010年建成投運，主要由集輸、加熱、消防、原油穩(wěn)定、采出水處理五大系統(tǒng)組成，日處理原油約1400m3，日外輸凈化油700m3，日處理采出水400 m3。在采油、集輸、加熱、消防、原油穩(wěn)定、采出水處理的過程中，井口揮發(fā)、儲罐“呼吸”及管線泄漏等原因都會產(chǎn)生污染性油田廢氣。

本文以作業(yè)區(qū)內(nèi)的綜合樓和宿舍樓為主要研究對象，采用計算流體動力學的方法，通過數(shù)值計算得出建筑物周圍PM10的擴散及分布情況進行預測，并根據(jù)結(jié)果提出了改善污染情況的建議。

1 數(shù)值計算模型

1.1 建筑概況

環(huán)江作業(yè)區(qū)前線生產(chǎn)保障點位于甘肅省環(huán)縣四合塬鄉(xiāng)，根據(jù)當?shù)貧庀髼l件，分別對夏季工況和冬季工況進行了模擬計算。氣象參數(shù)如下：夏季平均溫度約為25℃，平均風速2.4 m/s，風向西南偏南；冬季平均溫度約為-5℃，平均風速2.2 m/s，風向西北偏北。適當簡化后的建筑模型如圖1所示，建筑區(qū)域的大小為100.6 m×90.6 m。主要研究對象為綜合樓（A）（長61.4 m×寬17.3 m×高15 m）和公寓樓（B）（長39.2 m×寬19.4 m×高7.5 m）。另外，參考COST Action 732[6]設置數(shù)值模擬的計算域，入口邊界、上空邊界和側(cè)面邊界的距離均取為5Hmax，出口邊界的距離取為12Hmax，其中Hmax為15 m。

圖1 建筑模型

1.2 數(shù)學模型

常溫空氣可看作不可壓縮的黏性流體，低速湍流時，符合Boussinesq假設。因此，本文選取標準k—ε紊流模型[7]，采用SIMPLE算法對問題進行模擬求解。

1）不可壓縮流體紊流時均流動的連續(xù)性方程：

式中：uˉi—— 時均速度，m/s；

xi——坐標。

2）雷諾方程：

式中：uˉj—— 時均速度，m/s；

ui——瞬時速度，m/s；

uˉ′i、uˉ′j—— 脈動速度，m/s；

xj——坐標；

t——時間，s；

ρ——密度，kg/m3；

p——壓強，Pa；

μ——黏度，N·s/m2。

3

）湍流黏性系數(shù)方程：

式中：ηt——湍流黏性系數(shù)，N·s/m2；

Cμ——模型常數(shù)，取0.09；

k——湍流中單位質(zhì)量流體脈動動能；

ε——湍流中單位質(zhì)量流體脈動動能的耗散率。4）湍動能k的輸運方程：

式中：uˉk—— 時均速度，m/s；

xk——坐標；

Ck——模型常數(shù)，取0.09～0.11。

5）湍動能耗散率ε的輸運方程：

式中：ν——流體的運動黏性系數(shù)，m2/s；

Cε——模型常數(shù)，取0.07～0.09；

Cε1——模型常數(shù)，取1.41～1.45；

Cε2——模型常數(shù)，取1.91～1.92。

6）離散相方程：

式中：up——顆粒速度，m/s；

FD——顆粒的曳力，m/s2；

u——流體速度，m/s；

gx——重力加速度，m/s2；

ρ——流體密度，kg/m3；

ρp——顆粒密度，kg/m3；

Fx——顆粒所受的其他力，m/s2。

1.3 邊界條件

2）出口邊界條件：出流面設為自由流動邊界。

3）上空及兩側(cè)邊界條件：計算區(qū)域的上空及兩側(cè)設為對稱邊界。

4）地面及建筑物壁面邊界條件：地面、墻面及屋頂設為固體壁面，用壁面函數(shù)法并選取標準壁面函數(shù)進行計算。

1.4 污染物分析

環(huán)江作業(yè)區(qū)所處環(huán)境條件下，附著有油田廢氣的沙塵隨風飛揚為當?shù)刂饕拇髿馕廴拘问?，而其中PM10對人體的危害尤甚；因此，針對來流大氣中的PM10進行計算。由于慶陽地區(qū)沒有實測數(shù)據(jù)，所以設置來流中PM10的濃度為150 μg/m3，此值為國家二級標準[9]（即國家規(guī)定空氣質(zhì)量達標的最低標準）?？諝赓|(zhì)量分指數(shù)與污染物（PM10）濃度的對應關系見表1。

表1 空氣質(zhì)量分指數(shù)與污染物濃度限值

2 計算結(jié)果與分析

2.1 夏季工況

2.1.1 夏季風場壓力分析

夏季風向為西南偏南，平均風速2.4 m/s（圖2）。由圖2可知，由于圍墻距離公寓樓過近，導致公寓樓上風向和下風向的壓差大于5 Pa，會出現(xiàn)局部漩渦和死角，對行人造成安全隱患，不利于自然通風，易引起氣態(tài)污染物的積聚。

圖2 夏季風壓云圖（Z=1.5 m）

2.1.2 夏季濃度場分析

夏季工況，PM10在不同高度的分布情況見圖3和圖4。其中，圖3為1層工作人員坐姿呼吸高度的截面圖，圖4為綜合樓2—5層工作人員坐姿呼吸高度的截面圖。圖例中的數(shù)值單位是kg/m3，其中的6E-07 kg/m3對應的空氣質(zhì)量分指數(shù)為500，也就是所謂的污染程度“爆表”，因此圖中綠色、黃色、紅色區(qū)域污染情況非常嚴重。

由圖3可知，在綜合樓的東半部分，尤其是北側(cè)區(qū)域，PM10的濃度非常高，不僅會對行人的呼吸造成影響，而且能夠通過門窗進入一樓的辦公區(qū)域。公寓樓東側(cè)污染也比較嚴重，附著有污染性氣體的顆粒物會通過側(cè)門進入公寓樓內(nèi)，危害工作人員的健康。

由圖4可知，隨著高度的增加污染情況有所減緩，但是綜合樓的北側(cè)和東側(cè)墻外PM10濃度高，2—5層均受到污染，且污染情況比較嚴重。公寓樓的東側(cè)和南側(cè)墻外PM10濃度較高，2—4層均受到污染。

圖3 夏季PM10濃度云圖（Z=1.5 m）

圖4 夏季PM10濃度云圖（2—5層）

綜上所述，綜合樓和公寓樓受到PM10的污染都比較嚴重，綜合樓尤甚。由于PM10受到重力影響，有沉降現(xiàn)象，所以污染程度隨著高度的增加有所減弱。因此，在綜合樓的北側(cè)和公寓樓東側(cè)種植一排靠近建筑的樹木，即可在很大程度上阻隔顆粒物的擴散，阻止其進入綜合樓。另外，可以看到圍墻對于污染物有明顯的阻隔作用，適當增加圍墻高度可以改善建筑區(qū)域的空氣質(zhì)量。

為了保證工作人員的身體健康和居住舒適性，當受實際情況所限無法在在污染嚴重的區(qū)域種植樹木和適當加高圍墻時，也可以考慮在建筑物內(nèi)部采取空氣凈化措施或加裝機械送風系統(tǒng)。加裝機械送風系統(tǒng)時，需將進風口設在空氣質(zhì)量較好的區(qū)域（綜合樓設在南側(cè)，公寓樓設在西側(cè)為宜），使建筑物內(nèi)保持正壓送風，防止污染物侵入室內(nèi)。

2.2 冬季工況

2.2.1 冬季風場壓力分析

冬季風向為西北偏北，平均風速2.2 m/s（圖5）。由圖5可知，綜合樓北側(cè)和南側(cè)以及宿舍樓和食堂之間的風壓差較大，影響行人的安全和舒適性，也可能會影響一樓房間的熱舒適[10]。為減小壓差，同樣可以采取夏季減弱污染的策略，即在綜合樓的北側(cè)種植一排靠近建筑的樹木。另外，公寓樓入口處與食堂南側(cè)之間風壓變化急劇，對行人造成嚴重的安全隱患，應當在兩者之間種植樹木，減小風壓差。

圖5 冬季風壓云圖（Z=1.5 m）

2.2.2 冬季濃度場分析

冬季工況，PM10在不同高度的分布情況見圖6和圖7。其中，圖6為1層工作人員坐姿呼吸高度的截面圖，圖7為綜合樓2—5層工作人員坐姿呼吸高度的截面圖。

由圖6可知，冬季工況，綜合樓的污染情況較夏季輕，但污染物同樣集中在建筑壁面附近，容易進入建筑物內(nèi)。在公寓樓的東側(cè)，PM10大范圍聚集，并有明顯的渦流現(xiàn)象，污染物不易擴散，與夏季情況相同。圍墻對于污染物的阻隔作用，在冬季依然有所體現(xiàn)。

圖6 冬季PM10濃度云圖（Z=1.5 m）

由圖7可知，綜合樓的東南角墻外PM10濃度高，2—5層均受到污染，但污染范圍較小。公寓樓的東側(cè)和南側(cè)墻外PM10濃度較高，2—4層均受到污染。

圖7 冬季PM10濃度云圖（2—5層）

與夏季相比，冬季污染區(qū)域相似且情況較輕，主要是由于冬季地面溫度較低，浮升力弱于夏季；因此，改善措施主要考慮夏季工況，兼顧冬季即可。

3 結(jié)論

通過模擬計算的結(jié)果可知，粒徑小于或等于10 μm的顆粒物（PM10）一方面受浮升力的影響能夠長時間懸浮于空氣中，另一方面也會在重力作用下有一定的沉積現(xiàn)象，同時也會隨氣流的“爬墻”效應聚集于建筑壁面附近。

基于PM10的擴散特性，油田附近附著有油田廢氣的沙塵隨氣流擴散時，遇到建筑物便有一定的聚集，導致建筑附近空氣質(zhì)量嚴重惡化；一旦氣流進入建筑內(nèi)部，就會對室內(nèi)人員產(chǎn)生影響，易引起呼吸系統(tǒng)相關的疾病。

所以，為了有效減弱PM10的污染程度，保證工作人員身體健康，提出以下建議：

1）應在以后的油田前線生產(chǎn)保障點前期選址工作中將室內(nèi)外空氣質(zhì)量和建筑周圍風環(huán)境的因素考慮在內(nèi)綜合評估，設計人員應提高對保障點建筑室外顆粒物污染的重視程度，應通過模擬計算來驗證設計方案的可行性。

2）對于已建的建筑室外顆粒物污染較嚴重的油田前線生產(chǎn)保障點，可以進行模擬計算對污染情況進行分析，通過在適當?shù)牡胤皆黾訃鷫Ω叨群头N植樹木的方法可以有效改善污染情況。增加圍墻高度能夠有效地降低進入院內(nèi)的污染物濃度；種植樹木不僅能夠阻隔污染物的擴散、降低污染物濃度，而且能夠減小風壓、改善風環(huán)境、保證行人的安全和舒適性。

3）對于某些由于生產(chǎn)需要無法兼顧室內(nèi)外空氣質(zhì)量和建筑周圍風環(huán)境的新建保障點，或已建的建筑室外顆粒物污染較嚴重的保障點治理較困難時，可以考慮購買空氣凈化裝置或加裝正壓送風系統(tǒng)，提高室內(nèi)空氣品質(zhì)，改善工作人員的呼吸環(huán)境。

[1]閆炎，趙昕奕，周力平.近50年中國西北地區(qū)干濕演變的時空特征及其可能成因探討[J].干旱區(qū)資源與環(huán)境，2010，24（4）：38-41.

[2]POPE C A.Review：Epidemiological Basis for Particulate Air Pollution Health Standards[J].Aerosol Science and Technology，2000，32（1）：4-14.

[3]吳小華，張士權，龍鳳樂，等.石油開采中總烴對大氣環(huán)境影響的研究[J].油氣田環(huán)境保護，1996（4）：32-35.

[4]蔡建程.流動條件下SiO2顆粒表面甲醛及NO吸附特性實驗研究[D].西安：西安交通大學，2012.

[5]BACHMANN J D，DAMBERG R J，CALDWELL J C，et al.Review of the National Ambient Air Quality Standards for Ozone：Policy AssessmentofScientific and Technical Information[J].Proceedings of SPIE- The International Society for Optical Engineering，1996，8557（1）：85571V-6.

[6]FRANKE J L.The COST 732 Best Practice Guideline for CFD simulation of flows in the urban environment：a summary[J].International Journal of Environment&Pollution，2011，44（1-2）：419-427.

[7]陶文銓.數(shù)值傳熱學[M].2版.西安：西安交通大學出版社，2001：347-353.

[8]中國建筑設計研究院.建筑結(jié)構荷載規(guī)范：GB 50009—2012[S].北京：中國建筑工業(yè)出版社，2012：31-33.

[9]中國環(huán)境科學研究院，中國環(huán)境監(jiān)測總站.環(huán)境空氣質(zhì)量標準：GB 3095—2012[S].北京：中國環(huán)境科學出版社，2012：3.

[10]住房和城鄉(xiāng)建設部科技發(fā)展促進中心.綠色辦公建筑評價標準：GB/T 50908—2013[S].北京：中國建筑工業(yè)出版社，2013：7-8.

2017-05-22

（編輯 李發(fā)榮）

10.3969/j.issn.2095-1493.2017.10.016

張海曦，工程師，2013年畢業(yè)于西安交通大學建筑環(huán)境與設備工程系，從事油氣田地面工程暖通設計工作，E-mail：bravolfc@126.com，地址：陜西省西安市未央?yún)^(qū)長慶興隆園小區(qū)，710018。