袁新梅，周思柱，黃天成, ，華 劍，劉 奔

（1.長(zhǎng)江大學(xué) 機(jī)械工程學(xué)院，荊州 434023；2.中國(guó)石油大學(xué) 機(jī)械與儲(chǔ)運(yùn)工程學(xué)院，北京 102249；3.中鐵重工有限公司，武漢 430063）

0 引言

混砂車是實(shí)現(xiàn)按比例混砂，并能按壓裂工藝要求有效地向壓裂車供應(yīng)不同要求壓裂液的專用設(shè)備[1]，其主要作用是將液體和支撐劑按一定比例混合后，向壓裂車輸送，經(jīng)壓裂泵加壓后擠入井底巖層。混砂車的現(xiàn)場(chǎng)工況是要求大排量、大功率的[2]，這也對(duì)攪拌罐的整體結(jié)構(gòu)提出了新的要求。

本文的研究目的是通過(guò)攪拌罐的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化改進(jìn)，提高攪拌罐的混砂效率，降低攪拌功率，以得到攪拌效果較好（攪拌均勻度較高）的攪拌槳參數(shù)及操作參數(shù)的組合。

1 數(shù)值模擬過(guò)程

1.1 模型建立

混砂車使用的攪拌罐為雙層罐體，如圖1所示，水通過(guò)中間的入口進(jìn)入攪拌罐的外腔空間，然后分別由上中下三層的進(jìn)水口進(jìn)入攪拌罐內(nèi)腔。中上層部分由彎頭入水口進(jìn)入攪拌罐內(nèi)腔，其中中層入水口與重力方向成一定角度，底層入水口由罐底若干個(gè)孔進(jìn)入攪拌罐，通過(guò)這些孔進(jìn)水有以下兩個(gè)作用：一是從底部進(jìn)水沖刷沉積在罐底的砂子，二是在不用進(jìn)行混砂作業(yè)的時(shí)候，沖洗攪拌罐，經(jīng)過(guò)攪拌后由出口排出罐外，攪拌罐的罐底為球面形狀。

采用不同尺寸非結(jié)構(gòu)化網(wǎng)格對(duì)攪拌罐進(jìn)行網(wǎng)格劃分（圖2），槳葉區(qū)采用較細(xì)密的四面體網(wǎng)格，罐內(nèi)其他區(qū)域采用尺寸稍大的四面體網(wǎng)格。不同結(jié)構(gòu)尺寸的模型網(wǎng)格數(shù)量稍有不同，但其網(wǎng)格總數(shù)都在90萬(wàn)左右。

圖2 攪拌罐網(wǎng)格模型

1.2 數(shù)值求解

用于模擬計(jì)算的工作介質(zhì)為水和石英砂的混合物，其中石英砂的堆積密度為1640kg/m3。采用RNG ek- 湍流模型。旋轉(zhuǎn)葉輪和靜止罐體之間的耦合采用了滑移網(wǎng)格模型。入水口和入砂口都為速度入口，混合出口設(shè)為壓力出口。采用PISO算法，時(shí)間步長(zhǎng)0.002s，監(jiān)控出口流體的密度，達(dá)到穩(wěn)定后可以判斷解的收斂性[3～6]。

2 試驗(yàn)方案設(shè)計(jì)及試驗(yàn)結(jié)果

在多因素方差分析中，對(duì)各因素的每一種水平組合都要進(jìn)行試驗(yàn)，然而全面試驗(yàn)的試驗(yàn)次數(shù)又太多，為了解決這一問(wèn)題正交試驗(yàn)設(shè)計(jì)便應(yīng)運(yùn)而生。正交試驗(yàn)設(shè)計(jì)是利用按一定規(guī)律設(shè)計(jì)的正交表來(lái)安排試驗(yàn)方案進(jìn)行試驗(yàn)，并對(duì)試驗(yàn)結(jié)果進(jìn)行有效的分析計(jì)算，得出最優(yōu)方案條件的一種科學(xué)設(shè)計(jì)方法[7]。它只需選用一部分組合去做試驗(yàn)，不僅使試驗(yàn)次數(shù)盡可能少而且還能得到所需要的結(jié)論，正交試驗(yàn)設(shè)計(jì)已在實(shí)際中得到廣泛的使用，并證明其是安排多因素試驗(yàn)、尋求最優(yōu)水平組合的十分高效的試驗(yàn)設(shè)計(jì)方法。

2.1 試探性正交試驗(yàn)方案設(shè)計(jì)

由于攪拌罐的各因素對(duì)攪拌功率和攪拌時(shí)間的影響程度并不清楚，而且對(duì)各因素間是否有交互作用也不明確，因此，根據(jù)前人研究的相關(guān)情況，結(jié)合本課題的攪拌槽的結(jié)構(gòu)型式，攪拌罐的設(shè)計(jì)專業(yè)知識(shí)，確定考慮的攪拌葉輪和混合液體的物性參數(shù)有以下6個(gè)：A-上輪直徑，B-下輪直徑，C-上下輪距，D-葉輪罐底距，E-葉輪轉(zhuǎn)速，F(xiàn)-含砂比。

該試驗(yàn)除考察以上6個(gè)因素外，還要研究15個(gè)交互作用A×B，A×C，A×D，A×E，A×F，B×C，B×D，B×E，B×F，C×D，C×E，C×F，D×E，D×F，E×F對(duì)試驗(yàn)指標(biāo)的影響。選用L32(231)正交表，并根據(jù)L32(231)相應(yīng)的交互作用列表確定試驗(yàn)表頭并安排試驗(yàn)，試驗(yàn)表頭如表1所示。

表1 試驗(yàn)表頭設(shè)計(jì)

續(xù)（表1）

表2 試驗(yàn)結(jié)果數(shù)據(jù)

2.2 試探性正交試驗(yàn)結(jié)果

表2為按照正交試驗(yàn)方案模擬取得的試驗(yàn)數(shù)據(jù)，攪拌時(shí)間可以直接從軟件的模擬計(jì)算時(shí)間獲取，對(duì)于攪拌功率的計(jì)算，因?yàn)閺能浖@取的直接數(shù)據(jù)為攪拌槳的扭矩，且每次試驗(yàn)的攪拌槳轉(zhuǎn)速已知，因此攪拌功率可由下式計(jì)算得到：

其中：

P為功率，單位W；

M為扭矩，單位Nm；

n為轉(zhuǎn)速，單位rpm。

數(shù)據(jù)的處理利用正交試驗(yàn)設(shè)計(jì)的直觀分析和方差分析，試驗(yàn)指標(biāo)有顯著性影響因素且按主次順序排列為：

攪拌扭矩：

E ＞ B ＞ F ＞ BE ＞ A ＞ EF ＞ AE ＞ BF

攪拌時(shí)間：

EF ＞ BF ＞ E

對(duì)攪拌功率的方差分析R方=0.999（調(diào)整R方=0.996），試驗(yàn)結(jié)果的擬合相當(dāng)準(zhǔn)確，說(shuō)明在攪拌罐的結(jié)構(gòu)及攪拌操作的各參數(shù)對(duì)攪拌功率的影響已經(jīng)全部考慮，不存在遺漏的狀況，開始假設(shè)的參數(shù)包含了決定攪拌功率的主要參數(shù)。

從攪拌時(shí)間的方差分析R方=0.829（調(diào)整R方=0.469），可知試驗(yàn)結(jié)果的擬合是不太令人滿意的，說(shuō)明對(duì)攪拌時(shí)間有影響的數(shù)并沒有全部包含在開始時(shí)初選的參數(shù)中，各組試驗(yàn)中還有除變化因素外的其他因素在變動(dòng)。

為了尋求被忽略的因素，將多個(gè)試驗(yàn)的流場(chǎng)及混合密度圖對(duì)比觀察分析，發(fā)現(xiàn)有部分試驗(yàn)密度云圖顯示從入砂口進(jìn)入的砂粒只有少部分被吸入上槳的圓環(huán)之中，使攪拌槳的中間圓環(huán)失去了其導(dǎo)砂的作用，降低混合效率是顯然的，分析其原因，一是因攪拌槳中間圓環(huán)直徑較小，加上轉(zhuǎn)速較低（有的試驗(yàn)轉(zhuǎn)速為100rpm），使圓環(huán)中的壓降較小，沒有足夠的“吸”砂能力，二是由于上下攪拌槳的距離過(guò)大，上下槳的流場(chǎng)沒有能夠很好的連接起來(lái)，使大部分砂粒在上下槳之間被擴(kuò)散，降低了攪拌混砂的整體效率。而且混合罐的進(jìn)水口高度與攪拌槳葉高度之間的關(guān)系對(duì)罐內(nèi)液體的攪拌混合效果也有有較大影響。

2.3 攪拌槳優(yōu)化試驗(yàn)方案

攪拌轉(zhuǎn)速對(duì)攪拌功率的影響比攪拌槳結(jié)構(gòu)要大許多（重要程度按方差分析結(jié)果，大一個(gè)數(shù)量級(jí)）。對(duì)于混砂比，它與攪拌槳結(jié)構(gòu)參數(shù)間的交互作用對(duì)攪拌功率影響較?。ㄖ匾潭劝捶讲罘治鼋Y(jié)果，小一個(gè)數(shù)量級(jí)），為了詳細(xì)探究攪拌槳結(jié)構(gòu)參數(shù)對(duì)攪拌時(shí)間及攪拌功率的影響，故將攪拌轉(zhuǎn)速和混砂比在這里取常用工況值。

為了彌補(bǔ)在試探性分析中考慮因素遺漏的情況，著重細(xì)化分析攪拌槳的結(jié)構(gòu)參數(shù)對(duì)攪拌時(shí)間和攪拌功率的影響，將攪拌槳葉和槳葉中間的圓柱筒體的結(jié)構(gòu)尺寸分為兩因素考慮，如圖3所示。確定每個(gè)因素均取四個(gè)水平。擬考慮的攪拌槳結(jié)構(gòu)參數(shù)為：A上槳葉片長(zhǎng)度，B上槳圓環(huán)直徑，C下槳葉片長(zhǎng)度，D下槳圓環(huán)直徑，E上下槳距離，F(xiàn)下槳罐底距。

圖3 攪拌槳三維結(jié)構(gòu)圖

由試探性正交試驗(yàn)結(jié)果的結(jié)論，上下槳距離和下槳罐底距對(duì)攪拌功率和攪拌時(shí)間都沒有顯著影響，因此上下槳間的距離應(yīng)根據(jù)罐內(nèi)的流場(chǎng)分布確定，使上下槳脈動(dòng)速度高的區(qū)域能連接起時(shí)，攪拌罐的混合效率較高。對(duì)于下槳罐底距仍需要考慮其以下兩點(diǎn)影響：1）下層葉輪所產(chǎn)生的流場(chǎng)對(duì)罐底液體的混合狀態(tài)影響，2）從攪拌槳中間圓環(huán)導(dǎo)出的砂粒需要在下槳與罐底之間得到一定程度的分散混合后，從罐底的出口流出，才會(huì)使出口的攜砂液更為均勻。

同樣由試探性正交試驗(yàn)結(jié)果的結(jié)論知，攪拌結(jié)構(gòu)參數(shù)間的交互作用對(duì)攪拌時(shí)間和攪拌功率的影響并不顯著，因此對(duì)本正交試驗(yàn)不再考慮各因素間的交互作用。選用L16(45)正交表。

表3 正交試驗(yàn)結(jié)果

2.4 攪拌槳優(yōu)化試驗(yàn)結(jié)果

表3所示為攪拌槳正交試驗(yàn)的試驗(yàn)數(shù)據(jù)，對(duì)正交試驗(yàn)的結(jié)果進(jìn)行方差分析：攪拌功率R方=0.981（調(diào)整R方=0.906），對(duì)于攪拌時(shí)間R方=0.994（調(diào)整R方=0.970），可知試驗(yàn)結(jié)果的擬合是相當(dāng)準(zhǔn)確的，試驗(yàn)參數(shù)對(duì)攪拌功率都有顯著影響。

正交試驗(yàn)設(shè)計(jì)的直觀分析和方差分析的結(jié)果，可以看出對(duì)試驗(yàn)指標(biāo)具有顯著性影響因素且按主次順序排列為：

攪拌扭矩：C B A D

攪拌時(shí)間：B A C D

并按試驗(yàn)指標(biāo)值越小越好的要求優(yōu)選因素水平值，得到優(yōu)方案。將得到的優(yōu)方案結(jié)果再進(jìn)行模擬試驗(yàn)，其他條件不變，得到的試驗(yàn)結(jié)果為：攪拌扭矩：1298NM，攪拌功率：27.2kW，攪拌時(shí)間9.0s，由于優(yōu)化的目標(biāo)是以攪拌時(shí)間為標(biāo)準(zhǔn)，對(duì)比正交表中各試驗(yàn)的結(jié)果，優(yōu)化后的攪拌時(shí)間是最小的，因此本優(yōu)化過(guò)程可信。

3 結(jié)論

1）上下槳葉間的距離與下槳罐底距在選定的取值范圍內(nèi)對(duì)攪拌時(shí)間和攪拌功率無(wú)顯著影響，其取值可以根據(jù)罐內(nèi)的流場(chǎng)分布確定；

2）對(duì)于槳葉的結(jié)構(gòu)參數(shù)（上下槳直徑、上下槳距、下槳罐底距）來(lái)說(shuō)，它們之間的交互作用對(duì)攪拌時(shí)間和攪拌功率都無(wú)顯著影響。對(duì)于攪拌的操作條件（轉(zhuǎn)速與混合液濃度），它們之間的交互作用對(duì)攪拌時(shí)間和攪拌功率有都顯著影響；

3）對(duì)比攪拌罐優(yōu)化前后的結(jié)構(gòu)和操作參數(shù)及其對(duì)應(yīng)的模擬結(jié)果，主要的改進(jìn)為上攪拌槳、攪拌槳在罐內(nèi)的安裝高度等，對(duì)比改進(jìn)前后的結(jié)果，當(dāng)達(dá)到相似的混合效果，混砂濃度為40%時(shí)，攪拌功率降低了46.4%，攪拌時(shí)間降低了11.2%，性能提升明顯。

[1] 胡博仲,蘭中孝.大慶油田水力壓裂工程[M].北京:石油工業(yè)出版社,2008.

[2] 鄧斌奇,陳浩.混砂車的發(fā)展趨勢(shì)研究[J].中國(guó)石油和化工標(biāo)準(zhǔn)與質(zhì)量,2011,31(7):193.

[3] 王福軍.計(jì)算流體動(dòng)力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2004.

[4] 周思柱,陳翔,吳漢川,等.混砂車攪拌罐流場(chǎng)分析方法[J].石油機(jī)械,2011,39（增刊）:23-24.

[5] 周思柱,劉奔,華劍,等.基于均勻設(shè)計(jì)的混砂車攪拌槳結(jié)構(gòu)改進(jìn)[J].機(jī)械設(shè)計(jì)與制造,2013,5:120-122.

[6] 永斌.攪拌槽內(nèi)固液懸浮特性的試驗(yàn)研究[J].機(jī)械工程學(xué)報(bào),2011,47(6):186-192.

[7] 邱軼兵.試驗(yàn)設(shè)計(jì)與數(shù)據(jù)處理[M].北京:中國(guó)科學(xué)技術(shù)大學(xué)出版社,2008.