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

張海曦 李燁 邱稼儀

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

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

張海曦1李燁1邱稼儀2

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

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

室外顆粒物；計算流體動力學；空氣質量；建筑環境

引言

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

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

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

1 數值計算模型

1.1 建筑概況

環江作業區前線生產保障點位于甘肅省環縣四合塬鄉，根據當地氣象條件，分別對夏季工況和冬季工況進行了模擬計算。氣象參數如下：夏季平均溫度約為25℃，平均風速2.4 m/s，風向西南偏南；冬季平均溫度約為-5℃，平均風速2.2 m/s，風向西北偏北。適當簡化后的建筑模型如圖1所示，建筑區域的大小為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]設置數值模擬的計算域，入口邊界、上空邊界和側面邊界的距離均取為5Hmax，出口邊界的距離取為12Hmax，其中Hmax為15 m。

圖1 建筑模型

1.2 數學模型

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

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

式中：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

）湍流黏性系數方程：

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

Cμ——模型常數，取0.09；

k——湍流中單位質量流體脈動動能；

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

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

xk——坐標；

Ck——模型常數，取0.09～0.11。

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

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

Cε——模型常數，取0.07～0.09；

Cε1——模型常數，取1.41～1.45；

Cε2——模型常數，取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）上空及兩側邊界條件：計算區域的上空及兩側設為對稱邊界。

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

1.4 污染物分析

環江作業區所處環境條件下，附著有油田廢氣的沙塵隨風飛揚為當地主要的大氣污染形式，而其中PM10對人體的危害尤甚；因此，針對來流大氣中的PM10進行計算。由于慶陽地區沒有實測數據，所以設置來流中PM10的濃度為150 μg/m3，此值為國家二級標準[9]（即國家規定空氣質量達標的最低標準）。空氣質量分指數與污染物（PM10）濃度的對應關系見表1。

表1 空氣質量分指數與污染物濃度限值

2 計算結果與分析

2.1 夏季工況

2.1.1 夏季風場壓力分析

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

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

2.1.2 夏季濃度場分析

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

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

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

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

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

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

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

2.2 冬季工況

2.2.1 冬季風場壓力分析

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

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

2.2.2 冬季濃度場分析

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

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

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

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

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

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

3 結論

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

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

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

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

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

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

[1]閆炎，趙昕奕，周力平.近50年中國西北地區干濕演變的時空特征及其可能成因探討[J].干旱區資源與環境，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]吳小華，張士權，龍鳳樂，等.石油開采中總烴對大氣環境影響的研究[J].油氣田環境保護，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]陶文銓.數值傳熱學[M].2版.西安：西安交通大學出版社，2001：347-353.

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

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

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

2017-05-22

（編輯 李發榮）

10.3969/j.issn.2095-1493.2017.10.016

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