油田含油污水的污泥沉降及油液上浮模型

王 妍 曹柏寒 Sam REIFSNYDER 高 勝 Diego ROSSO

（1.東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院；2.加利福尼亞大學(xué)歐文分校土木與環(huán)境工程系）

注水注聚采油是我國油田開發(fā)的主要方式，這使得采出液中含有大量的顆粒物、油液及化學(xué)藥劑等［1］。油田污水處理就是將轉(zhuǎn)油站、脫水站和聯(lián)合站處理后的含油污水通過沉降、過濾等手段凈化處理，以達(dá)到水質(zhì)指標(biāo)的要求，回注地下再利用［2］。顯然，其處理效果直接影響到油田的開采和環(huán)境保護(hù)［3］。開展含油污水沉降機(jī)理研究，建立污水全系統(tǒng)沉降模型并進(jìn)行模擬分析，將揭示污水處理的過程和規(guī)律，對(duì)提高油田污水處理效果意義重大。

目前，國內(nèi)對(duì)于沉降系統(tǒng)的研究，主要是通過流場(chǎng)分析的方法優(yōu)化沉降系統(tǒng)的結(jié)構(gòu)工藝來達(dá)到一定的凈化效果。 李杏杏利用數(shù)值模擬、流場(chǎng)實(shí)驗(yàn)等方式的結(jié)合，優(yōu)化改進(jìn)了混凝沉降罐的內(nèi)部結(jié)構(gòu)［4］。韓海倉等利用Fluent對(duì)沉降罐內(nèi)配液管不同開孔方式的流場(chǎng)進(jìn)行研究，分析沉降罐的油水分離效果［5］。 雖然沉降系統(tǒng)的研究取得了一定的進(jìn)展和效果，但學(xué)者們大多將城市污水的沉降理論和分析方法轉(zhuǎn)用于含油污水，沒有針對(duì)含油的特殊性研究油液的上浮機(jī)理，更鮮有建立系統(tǒng)整體的數(shù)學(xué)模型進(jìn)行模擬，整個(gè)含油污水的精確建模和模擬研究仍遠(yuǎn)落后于其他領(lǐng)域。 筆者通過與美國加利福尼亞大學(xué)歐文分校城市水研究中心合作， 首次將Takacs沉降模型改進(jìn)用于油液上浮的模擬， 力求建立完善的系統(tǒng)整體沉降模型，并結(jié)合現(xiàn)場(chǎng)和實(shí)驗(yàn)數(shù)據(jù)，深入分析仿真結(jié)果，驗(yàn)證模型對(duì)污水凈化的描述性。

1 模型依據(jù)

1.1 模型應(yīng)用依據(jù)

沉降是油田含油污水處理的首要環(huán)節(jié)，核心為重力沉降［6］。 Takacs模型是重力沉降的經(jīng)典模型，可用于模擬在穩(wěn)態(tài)或動(dòng)態(tài)下污水中各物質(zhì)的分布和沉降速度［7，8］。含油污水的沉降過程總體上為基本穩(wěn)態(tài)， 故Takacs模型對(duì)此具有適應(yīng)性［9，10］。但由于油田污水含油的特殊性，懸浮固體沉降與油液上浮同時(shí)發(fā)生， 兩者獨(dú)立進(jìn)行又彼此干涉，而油田中除油率和除懸率是評(píng)價(jià)凈化效果的兩個(gè)重要指標(biāo)。 顯然，建立油液上浮的數(shù)學(xué)模型，并考慮物質(zhì)間的彼此耦合和層間耦合特性，從整體上描述含油污水沉降過程和效果很有必要。

1.2 沉降系統(tǒng)的離散化

如圖1所示，按照經(jīng)驗(yàn)和實(shí)際沉降系統(tǒng)工況，將沉降罐劃分為10層，第3層為進(jìn)水口層，第1、7、10層分別為收油層、出水口層和污泥層，水流速度可根據(jù)各流量和沉降罐截面積計(jì)算得出。 污水流入后，分別向上和向下流動(dòng)。 對(duì)于任一層，接受著上一層下沉的懸浮固體和下一層上升的油液的同時(shí)，本層中的懸浮固體也向下一層沉降且油液向上一層上升。 層與層間存在緊密的耦合關(guān)系（圖2），在各層閾值濃度Xt干涉下，可以通過各層的通量J來描述，J是各物質(zhì)的濃度和速度的函數(shù)。

圖1 沉降罐離散化示意圖

圖2 模型層間耦合關(guān)系

2 模型建立

2.1 污泥沉降模型

2.2 油液上浮模型

油田含油污水有毒有害，不易降解且含油率高，將影響水質(zhì)穩(wěn)定劑的正常使用，阻止空氣中的氧溶解于水，因此除油率亦是評(píng)價(jià)污水處理效果的重要參數(shù)。 油液的凈化過程與懸浮固體的重力沉降剛好相反，一直處于上浮狀態(tài)且具有聚合性——油液在上浮過程中，不同粒徑的油滴不斷發(fā)生碰撞和融合，使其粒徑不斷改變，致使上浮速度和擴(kuò)散通量發(fā)生變化。 因此，必須綜合考慮以上因素，才能建立準(zhǔn)確的數(shù)學(xué)模型，描述油液的運(yùn)動(dòng)規(guī)律和耦合特性。

由于上浮和沉降的原理相似， 筆者仍將Takacs模型作為油液上浮模擬的基礎(chǔ)， 在此前提下，改進(jìn)Takacs模型，需引入油液聚集速度vc和油液聚集通量Jc：

在vc，i的計(jì)算中， 等號(hào)右側(cè)的第1項(xiàng)反映顆粒中徑較大且聚集快的油滴聚集速度， 第2項(xiàng)反映顆粒中徑較小且聚集慢的油滴聚集速度。

最終，建立油液上浮模型如下：

由式（14）～（20）可見，各層油液上浮模型中不只包含油液上浮通量Js，亦包含油液聚集通量Jc，此時(shí)Js仍用來描述各層間由于重力作用油液的耦合關(guān)系，Jc用來描述油液的聚集性帶來的層間耦合關(guān)系。

含油污水中還包含如鹽、硫化物及揮發(fā)酚等可溶性有機(jī)物。 油田只要求處理懸浮固體與油液的污水全部回注地下，確保污水不外排亦不做他用， 對(duì)油和懸浮物以外的物質(zhì)并無嚴(yán)格要求，筆者對(duì)此進(jìn)行簡化處理， 不做詳細(xì)建模和仿真分析。

3 仿真分析

3.1 參數(shù)設(shè)置

采用Simulink建立仿真模型， 并進(jìn)行GUI設(shè)計(jì)。 模型中最小沉降速度和最大上浮速度的取值直接決定懸浮固體、 油液的濃度和凈化效果，與研究對(duì)象的理化性質(zhì)密切相關(guān)，故理化數(shù)據(jù)基于大慶油田某污水處理廠實(shí)際技術(shù)參數(shù)和實(shí)驗(yàn)室樣本實(shí)測(cè)數(shù)據(jù)：

污油含量 194g/m3

懸浮物含量 216g/m3

其他有機(jī)物含量 406g/m3

最大固體沉降速度 321m/d

最大油液上升速度 298m/d

最大Vesilind沉降速度 505m/d

最大Vesilind上升速度 411m/d

受阻沉降區(qū)固體沉降參數(shù) 5.76×10-4m3/g

受阻沉降區(qū)油液上升參數(shù) 5.76×10-4m3/g

懸浮固體閾值濃度 4 000g/m3

油液閾值濃度 5 000g/m3

不可沉降參數(shù) 0.002 3

最大油液濃度 840 000g/m3

最大油液聚集速度 6m/d

離散區(qū)沉降參數(shù) 0.002 9m3/g

離散區(qū)聚集參數(shù) 8×10-5g/m3

集合區(qū)聚集參數(shù) 9×10-5g/m3

3.2 仿真結(jié)果

3.2.1 實(shí)驗(yàn)設(shè)計(jì)

實(shí)驗(yàn)裝置主要由量筒、燒杯、攪拌棍、分液漏斗、漏斗、電子天平、加熱盤、真空抽濾瓶、濾紙、烘干爐、電磁攪拌器、激光粒度分析儀及高速攝影器等組成。 為了反映現(xiàn)場(chǎng)水質(zhì)，保證實(shí)驗(yàn)數(shù)據(jù)真實(shí)有效，所有實(shí)驗(yàn)均在采樣當(dāng)天進(jìn)行，具體過程為：

a. 測(cè)量各物質(zhì)濃度。 取200mL污水原液進(jìn)行實(shí)驗(yàn)。 初始稱重之后，按時(shí)間梯度靜置，并分別測(cè)量漂浮至頂部的油液濃度；利用真空抽濾管加濾紙，對(duì)污水原液進(jìn)行過濾實(shí)驗(yàn)，測(cè)量懸浮固體濃度。

b. 測(cè)量油液沉降參數(shù)。 分別取20、80、120、160mL油液與凈化后的污水樣品混合， 充分?jǐn)嚢韬螅?進(jìn)行靜置實(shí)驗(yàn)， 并在1、2、3、4、5、10、15、20、30、45、60min時(shí)用高速攝像機(jī)進(jìn)行拍攝記錄，測(cè)量油液在污水中的沉降參數(shù)。

c. 進(jìn)行粒徑分析。 利用激光粒度分析儀進(jìn)行粒徑分析，包括未經(jīng)處理的污水原液和經(jīng)過沉降后的樣品，來獲取懸浮固體與油液的部分粒徑參數(shù)。

3.2.2 仿真結(jié)果

所有測(cè)量的參數(shù)作為數(shù)學(xué)模型的基礎(chǔ)數(shù)據(jù)，建立懸浮固體和油液數(shù)學(xué)模型并進(jìn)行仿真模擬。仿真結(jié)果見表1，可見，采用持續(xù)4h沉降過程，懸浮固體含量從第1～10層逐層遞增，油液含量逐層遞減，第3～9層中懸浮固體含量在68～70g/m3，第8層上/下的油液含量均達(dá)到20g/m3左右，各層模擬計(jì)算值與實(shí)測(cè)值的相對(duì)誤差在6%以內(nèi)。

表1 懸浮固體和油液的濃度仿真與實(shí)測(cè)值 g/m3

由表2可知， 沉降系統(tǒng)的油液去除率和懸浮固體去除率與實(shí)測(cè)值的偏差率的絕對(duì)值控制在1.50%以內(nèi)。 因此，依據(jù)Takacs模型建立的油液上浮模型和污泥沉降模型能夠反映油田含油污水的真實(shí)沉降狀況，揭示沉降機(jī)理，而且懸浮固體的模擬效果稍優(yōu)于油液的模擬。

表2 油液去除率與懸浮固體去除率 %

3.3 仿真分析

3.3.1 仿真曲線

圖3為各層懸浮固體濃度變化曲線。 由圖3可知，懸浮固體除第10層污泥層外，總體都呈現(xiàn)出下降，最終達(dá)到恒定的趨勢(shì)。 第3層之上由于入水口影響了物質(zhì)交換， 導(dǎo)致濃度的初期波動(dòng)較大，其中第1、2層受到下方各層的油液上浮影響最大，懸浮固體含量逐漸降低；第3、4層先小幅升高隨后下降，并趨于穩(wěn)定，原因在于初期第1、2層下沉至第3層的懸浮固體含量多于第3層下沉至第4層的含量；第5～9層，由于位置比較靠下，受到上方各層的累積沉積效應(yīng)有所延后，初期出現(xiàn)穩(wěn)定期，層數(shù)越靠下，穩(wěn)定期越長；作為污泥層的第10層， 其懸浮固體的含量隨著時(shí)間一直在增加。

圖3 第1～10層中懸浮固體濃度變換曲線

圖4為各層油液濃度變化曲線。 由圖4可知，油液濃度總體呈現(xiàn)下降并恒定的規(guī)律。 第3層及之上出現(xiàn)了濃度的初期波動(dòng)，油液含量先小幅升高，隨后趨于穩(wěn)定，原因是第2、3層初期接受下層油液上浮的作用大于向第1層的上浮；第4～10層油液含量迅速下降，并在一定時(shí)間后趨于穩(wěn)定。

圖4 第1～10層中油液濃度變換曲線

對(duì)比可見，同樣的沉降時(shí)間，由于油液上浮速度較懸浮物下沉速度慢，使得趨于穩(wěn)態(tài)的油液凈化周期較懸浮固體凈化周期長。

3.3.2 非穩(wěn)態(tài)層仿真分析

非穩(wěn)態(tài)層由于與外界進(jìn)行物質(zhì)交換，物質(zhì)上浮或沉降態(tài)勢(shì)具有其特殊的變化規(guī)律。 分別針對(duì)4個(gè)非穩(wěn)態(tài)層進(jìn)行仿真分析。

第1層收油層， 主要受到下層油液的上浮作用，并將本層懸浮固體沉降至下層。 如圖5所示，懸浮固體濃度總體呈現(xiàn)一直下降的趨勢(shì)，因?yàn)闆]有上層污水的影響，短期內(nèi)以指數(shù)級(jí)迅速達(dá)到穩(wěn)態(tài)，初期的效果明顯優(yōu)于后期；該層的油液濃度呈現(xiàn)先上升后逐步下降趨于穩(wěn)定的趨勢(shì)，由于下面層數(shù)較多，初期接受下層油液上浮的累積效果起決定作用，含油率上升達(dá)到頂峰后油液溢流并逐步下降，最終達(dá)到穩(wěn)定，且達(dá)到穩(wěn)態(tài)的過程比較緩慢，這也是由于下面層數(shù)多和油液上浮速度慢所致。

圖5 第1層各物質(zhì)變換曲線

第3層進(jìn)水層， 與上下鄰層都存在懸浮固體和油液的雙交換。 如圖6所示，懸浮固體濃度曲線顯示出急速上升隨后下降的趨勢(shì)，該層初期接受上兩層懸浮固體的沉降量快于向下層的沉降，懸浮固體含量快速上升，上兩層影響減弱后達(dá)到穩(wěn)態(tài)，這個(gè)過程非常短，源于沉降速度較大和上面的層數(shù)較少，也可從圖3可見，出水口層之上均具有類似趨勢(shì)；油液濃度呈現(xiàn)先上升后逐步下降至穩(wěn)定的趨勢(shì)，初期下層油液的上浮作用大于向上兩層上浮的作用，此趨勢(shì)與第1層類似，但明顯凈化周期縮短很多，因?yàn)橄聦訉訑?shù)減少，而且本層油液亦可上浮至其他層。 另外，油液總體的凈化時(shí)間較懸浮固體長，是油液的上浮速度比懸浮固體沉降速度小所致。

圖6 第3層各物質(zhì)變換曲線

第7層出水層， 同樣與上下鄰層都受到懸浮固體和油液的交換作用。 如圖7所示，由于位置比較靠下， 受到上方各層的累積沉積效應(yīng)有所延后。 此外，隨時(shí)的出液和向下層沉降也會(huì)對(duì)懸浮固體濃度有影響，因此導(dǎo)致懸浮固體濃度在初期一段時(shí)間內(nèi)，處于緩慢上升后趨于基本穩(wěn)態(tài)的式樣，從圖3中亦可看出，越靠下的層，其達(dá)到穩(wěn)態(tài)的處理周期將越長；油液濃度與第1層、第3層及其以上各層一樣，均呈現(xiàn)先上升后逐步下降至穩(wěn)定的趨勢(shì)，但顯然上升的幅度和達(dá)到穩(wěn)態(tài)的時(shí)間都在逐步縮短，這是因?yàn)樵酵聦邮芨聦拥挠绊憣p少。

圖7 第7層各物質(zhì)濃度變換曲線

第10層污泥層，主要負(fù)責(zé)污泥沉降和收集處理，只與上層存在耦合交換，受上層懸浮固體的下沉作用，并將本層油液上浮至上層。 圖8為該層的懸浮固體和油液濃度變化都不再波動(dòng)，懸浮固體濃度呈持續(xù)增大趨勢(shì)且增大速率穩(wěn)定；油液濃度呈指數(shù)級(jí)下降趨勢(shì)，相對(duì)其他各層都快。

圖8 第10層各物質(zhì)濃度變換曲線

3.4 綜合評(píng)價(jià)

相對(duì)來看，油液上浮速度較懸浮固體沉降速度慢，穩(wěn)態(tài)和非穩(wěn)態(tài)層都體現(xiàn)出油液的凈化周期長的特點(diǎn)。 此外，從各模擬曲線均可見，懸浮固體和油液濃度最終都將趨于穩(wěn)態(tài)， 達(dá)到凈化效果。沉降過程中，沉降罐結(jié)構(gòu)、污水水質(zhì)和凈化時(shí)間都是影響凈化效果的核心因素，但要保證凈化效果，其中的時(shí)間因素仍然是影響凈化效果的關(guān)鍵參數(shù)。

4 結(jié)束語

油田含油污水處理效果直接關(guān)系到回注地下的水質(zhì)，影響油田的開采和環(huán)境保護(hù)。 在凈化水系統(tǒng)模擬中，提高沉降系統(tǒng)油液和懸浮固體建模和計(jì)算的準(zhǔn)確性，具有很大的實(shí)用價(jià)值。 通過分析研究， 筆者將Takacs模型應(yīng)用于油液上浮的模擬，增加了油液凝聚速度參數(shù)，建立了油液上浮模型和污泥沉降模型。 計(jì)算結(jié)果表明，兩種模型能夠反映油田含油污水的真實(shí)狀況，可揭示上浮和沉降機(jī)理。 并通過仿真模擬借助現(xiàn)場(chǎng)數(shù)據(jù)進(jìn)行深入分析，準(zhǔn)確預(yù)測(cè)污水凈化效果，為尋求高效、綠色且智能化的油田污水處理系統(tǒng)提供一定的理論依據(jù)。