梁海波 宋 洋 于志剛 辛小軍 楊 海

(1.西南石油大學(xué)機(jī)電工程學(xué)院 2.中海石油(中國)有限公司湛江分公司 3.中法渤海地質(zhì)服務(wù)有限公司湛江分公司)

0 引 言

在石油工程中，一口油氣井鉆井成功與否，在很大程度上取決于鉆井液的性質(zhì)和功能。鉆井液具有清洗井底、攜帶巖屑、冷卻潤滑鉆頭和鉆柱、傳遞水功率、控制和平衡地層壓力以及提高鉆井速度等作用[1]。隨著石油鉆井行業(yè)的發(fā)展，鉆井液需要滿足鉆井和完井等越來越多的功能需求，比如實時監(jiān)測安全密度窗口等[2]，及時準(zhǔn)確地監(jiān)測到鉆井液的工藝性能是很關(guān)鍵的一步。隨著鉆井自動化逐步實現(xiàn)，鉆井液性能測量受到業(yè)內(nèi)人士的關(guān)注，目前鉆井液性能數(shù)據(jù)主要來自于人工測量，鉆井液樣品要經(jīng)過收集、運輸、處理及分析，然后才能得到性能報告。人工測量不利于分析流體特性對預(yù)期問題的影響，而且也不能正確表征溫度和壓力對鉆井液性能的影響。

目前，針對鉆井液流變性測量，在鉆井現(xiàn)場使用比較廣泛的是馬氏漏斗黏度計和旋轉(zhuǎn)黏度計[3]。馬氏漏斗是API標(biāo)準(zhǔn)規(guī)定的測量鉆井液黏度的儀器，作為目前井場測量鉆井液黏度較為常見的手動測量儀器，它能夠簡單快速得到結(jié)果。但是其所得的流變參數(shù)只是粗略值，且無法實現(xiàn)測量的實時性，存在人工操作的誤差。旋轉(zhuǎn)黏度計是目前井場用于鉆井液黏度測量最常用的儀器，但其為手動操作儀器，缺乏測量的及時性，存在操作的誤差及需要一定的人工參與，不具備自動化測量的條件[4]。2008年，A.SAASEN等[5]為解決目前鉆井液無法自動測量的問題，在旋轉(zhuǎn)黏度計的基礎(chǔ)上，通過增加調(diào)速機(jī)構(gòu)進(jìn)行自動化測量，但是由于在旋轉(zhuǎn)黏度計中內(nèi)外筒的流體置換并讓流體處于穩(wěn)定狀態(tài)需要較長時間，無法實現(xiàn)實時在線測量。哈里伯頓公司開發(fā)出一款全自動鉆井液流變性測量設(shè)備——BaraLogix DRU，但其流變性測量周期為15 min，不具備測量的及時性與實時性。創(chuàng)升益世公司提出的DFS系統(tǒng)具備溫度跟隨系統(tǒng)，適用范圍廣，但仍舊無法解決測量周期過長的問題。

針對當(dāng)前鉆井液測量中無法實現(xiàn)在線測量、工作量大及存在操作誤差等問題，本文提出一種鉆井液流變性在線測量裝置，旨在實時準(zhǔn)確地監(jiān)測鉆井液流變性能，其測量參數(shù)包括密度、溫度、質(zhì)量流量、塑性黏度、動切力、稠度系數(shù)及流性指數(shù)等。通過測量鉆井液在不同管徑下的壓力損失[6]，結(jié)合質(zhì)量流量計所測得的數(shù)據(jù)，上傳到上位機(jī)進(jìn)行處理，從而得到鉆井液的各項流變性能參數(shù)。該鉆井液流變性在線測量裝置對鉆井完井自動化進(jìn)程有較好的推動作用，具有顯著的社會效益和經(jīng)濟(jì)效益。

1 理論方案模型

1.1 設(shè)計原理

用電動隔膜泵從鉆井液罐抽取一定體積鉆井液到測量裝置中，在通過濾網(wǎng)時濾去鉆井液內(nèi)雜質(zhì)后，經(jīng)阻尼器通過質(zhì)量流量計，將測量的數(shù)據(jù)值傳輸給上位機(jī)，判斷是否為穩(wěn)定流量后，對泵的流量進(jìn)行反饋調(diào)節(jié)[7]，保證輸出鉆井液流量恒定，同時質(zhì)量流量計測得溫度、密度及流量信號傳輸?shù)缴衔粰C(jī)。鉆井液經(jīng)過兩根不同直徑的測量管[8]，用一套平面膜差壓傳感器測量通過測量段的壓差損耗，將壓差信號傳輸?shù)缴衔粰C(jī)，輸出傳統(tǒng)六速測量結(jié)果，得到塑性黏度、動切力、靜切力、稠度系數(shù)以及流性指數(shù)等值，最后回到鉆井液罐。根據(jù)傳輸?shù)男盘柡蜏y量管尺寸，連續(xù)計算并記錄鉆井液流變性。測量裝置工作原理示意圖如圖1所示。

圖1 工作原理示意圖

1.2 數(shù)據(jù)處理算法

在試驗中，由于流量、溫度、密度及固相含量等都對壓差測量數(shù)據(jù)存在不可避免的波動干擾[9]，故使用數(shù)據(jù)處理算法對測量參數(shù)進(jìn)行修正，以提高流變性參數(shù)測量精度。本文采用偏最小二乘回歸方法對自變量數(shù)據(jù)進(jìn)行處理[10]。

偏最小二乘回歸方法是一種多元統(tǒng)計分析方法，能夠廣泛應(yīng)用在各個領(lǐng)域，其主要在變量之間存在的多重相關(guān)性上回歸建模，同時也將信息綜合和篩選應(yīng)用在了建模過程中，也就是一種將主成分分析、多元線性回歸方法融為一體的分析方法，其建模過程如下。

在自變量X和因變量Y中分別提取出第一個成分a1、b1，a1和b1包含著原數(shù)據(jù)中最大的變異信息，即有：

Var(a1)→max

(1)

Var(b1)→max

(2)

為了更好地解釋自變量與因變量之間的關(guān)系，a1和b1間還需要滿足相關(guān)程度達(dá)到最大，即有：

r(a1,b1)→max

(3)

這兩個成分不僅能夠最優(yōu)地代表各自的原數(shù)據(jù)，同時a1也對b1有最佳的解釋能力。

在偏最小二乘回歸方法中，可令成分間滿足

cov(a1,u1)=

(4)

分別求出X0對a1的回歸方程和Y0對b1的回歸方程：

(5)

針對第二成分a2的提取，則將X0用X1取代，將Y0用Y1取代，使用上述相同的方法進(jìn)行回歸。

以此作為基礎(chǔ)，當(dāng)X的秩為M時，就有：

X0=a1P1T+……+aMPMT

(6)

Y0=a1R1T+……+aMRMT+YM

(7)

由于a1，……，aM能夠以X01，……，X0P的線性組合方式來表達(dá)，所以將上式還原為：

YK=ck1X1+……+ckPXP+YMK(k=1,……,P)

(8)

式中：YMK為殘差矩陣YM的第k列。

1.3 數(shù)學(xué)模型

由差壓傳感器測得的測量管兩端壓差可以得到管道壁面的剪切應(yīng)力[11]：

(9)

式中：τw為管壁處的剪切應(yīng)力，Pa；D為測量管的內(nèi)徑，m；p為測量管測量段的壓差，Pa；l為測量管測量段的長度，m。

對于管道，管壁處的剪切速率為：

(10)

令：

(11)

則有：

(12)

式中：γw為管壁處的剪切速率，s-1；v為流體經(jīng)測量管的速度，m/s；N為無因次參數(shù)。

因流體流變學(xué)復(fù)雜性和紊流渦漩等因素導(dǎo)致非牛頓流體湍流的數(shù)學(xué)建模十分復(fù)雜，所以，通常引入摩擦因子項f。則用摩擦因數(shù)表示的壓力損失方程為[13]：

(13)

式中：ρ為鉆井液密度，kg/m3。

管道中流體為層流時，f的計算式為：

(14)

式中：Re為雷諾數(shù)，無量綱。

對于紊流，相關(guān)學(xué)者建立了經(jīng)驗關(guān)聯(lián)式來獲得摩擦因數(shù)。為了獲得牛頓流體在管道中湍流的摩擦因數(shù)，提出了幾種關(guān)聯(lián)式(如Colebrook 1939，Blasius 1913)，但對非牛頓流體進(jìn)行了一些試驗研究的Dodge和Metzner提出了非牛頓流體在光滑管道中湍流與時間無關(guān)的經(jīng)驗?zāi)Σ烈驍?shù)關(guān)聯(lián)式[12]：

(15)

式中：n為流動行為指數(shù)，無量綱。

雖然這種相關(guān)性提高了某些類型鉆井液流變性測量的精度，但這種相關(guān)性不適用于聚合物基鉆井液流變性的測量。因此，對于特定的鉆井液系統(tǒng)，更傾向于實時測量和獲取摩擦因數(shù)，而不是依賴于經(jīng)驗相關(guān)性，經(jīng)驗相關(guān)性僅在某些條件下適用。絕大多數(shù)鉆井液都具有觸變性，在很多基礎(chǔ)性的鉆井液試驗研究中，我們能夠觀察到鉆井液經(jīng)過晃動然后靜止并形成凝膠的形態(tài)，再經(jīng)過第二次晃動則恢復(fù)到初始形態(tài)。所以，如果從表面現(xiàn)象來講，鉆井液的觸變性是指在恒溫恒壓下，其搖動后變稀，靜止后變稠的特征。嚴(yán)格意義的觸變性定義為：在固定的剪切速率梯度下，流體的剪切應(yīng)力會呈現(xiàn)隨著時間的延長而逐漸減小的特性。從膠體化學(xué)的層面上講，觸變性是指在恒溫狀態(tài)下，流體能夠在溶膠與凝膠之間相互轉(zhuǎn)變的特性。

牛頓流體的流變方程為：

τ=μγ

(16)

式中：τ為剪切應(yīng)力，Pa；μ為表觀黏度，Pa·s；γ為剪切速率，s-1。

根據(jù)鉆井液的流體力學(xué)性質(zhì)，鉆井液大多屬于純黏性非牛頓流體，描述其流體力學(xué)性質(zhì)的一般本構(gòu)方程式為[14]：

(17)

式中：η為有效黏度，Pa·s。

目前井場鉆井液的流體模式有賓漢流體模式、假塑性流體模式、膨脹性流體模式、卡森流體模式[15]、赫-巴流體模式以及考慮溫度和壓力的流變模式，其對應(yīng)的流變方程如下。

一維簡單剪切流動情況下的賓漢流體模式為：

τ=τ0+μpγ

(18)

有效黏度函數(shù)為：

(19)

式中：τ0為動切力(屈服應(yīng)力)，Pa；μp為塑性黏度，Pa·s。

τ=Kγn(n<1)

(20)

式中：K為稠度系數(shù)，Pa·sn。

有效黏度函數(shù)為：

(21)

一維簡單剪切流動情況下膨脹性流體模式為：

τ=Kγn(n>1)

(22)

有效黏度函數(shù)同式(2)。

一維簡單剪切流動情況下的卡森流體模式為：

(23)

有效黏度函數(shù)為：

(24)

式中：η∞為卡森極限高剪黏度，mPa·s；τc為卡森動切應(yīng)力，Pa。

卡森流體模式的另一特征參數(shù)是剪切稀釋指數(shù)Im，Im值用于反映鉆井液剪切稀釋性的強(qiáng)弱，可表示為：

式中，i為第 i個產(chǎn)業(yè)；j為第 j個地區(qū)，Xij表示第j個地區(qū)的第 i產(chǎn)業(yè)的銷售收入或產(chǎn)值指標(biāo)；產(chǎn)業(yè)規(guī)模區(qū)位商也稱區(qū)域規(guī)模優(yōu)勢指數(shù)。在本文中，產(chǎn)業(yè)規(guī)模區(qū)位商表示該地區(qū)工業(yè)總產(chǎn)值占全省三大產(chǎn)業(yè)總產(chǎn)值的比重。當(dāng)產(chǎn)業(yè)規(guī)模區(qū)位商大于1時，表明該地區(qū)工業(yè)比較集中，一定程度上顯示出工業(yè)具有較強(qiáng)的競爭力，產(chǎn)業(yè)規(guī)模區(qū)位商越大，表示該地區(qū)工業(yè)集聚越顯著；產(chǎn)業(yè)規(guī)模區(qū)位商等于1時，表示該地區(qū)工業(yè)處于均勢，該產(chǎn)業(yè)的集中還不明顯；產(chǎn)業(yè)規(guī)模區(qū)位商小于1時，表明該產(chǎn)業(yè)在該地區(qū)沒有形成集中。

Im=[1+(100τc/η∞)0.5]2

(25)

一維簡單剪切流動情況下的赫-巴流體模式表達(dá)式為：

(26)

具有3個參數(shù)的赫-巴流體模式[16]在水力學(xué)計算上顯得更為復(fù)雜，但同時也可以明顯看出其具有相對于其他流變模式更高的精確度，故其使用廣泛性正在逐步提升。

W.C.MCMORDIE，JR等[17]在1974年提出了同時考慮溫度和壓力的油基鉆井液流變模式：

(27)

式中：T為溫度，A、B為常數(shù)。

對于具有屈服應(yīng)力的鉆井液，其流變模式方程可以改寫為：

(28)

針對不同的鉆井液選擇不同的流變模式，目前廣泛使用的多數(shù)鉆井液為賓漢流體和假塑性流體。根據(jù)流變曲線或數(shù)值擬合分析來決定采取哪種流變模式，選擇相關(guān)性好的流體模式來分析和表示實際鉆井液的流變行為[18]。

2 鉆井液流變性實時測量裝置

測量裝置整體框架采用不銹鋼矩形管搭建。測量管采用?16.0 mm×1.6 m和?20.0 mm×2.0 m，壁厚均為2 mm的不銹鋼圓管，管線之間采用法蘭連接方式，便于管線的拆裝和分段清理。電動隔膜泵出口端與質(zhì)量流量計間用不銹鋼圓管連接，電動隔膜泵入口端與鉆井液罐、測量出口端均用軟管連接。

整個裝置要求減振，泵與框架均不連接。橇裝時，在底座與泵體底座間加減振橡膠或其他性能良好的減振器。由于裝置距離井場振動源較近，框架存在振動，故在每個傳感器的固定點下方增加了減振設(shè)計。

鉆井液黏度受溫度影響較大，為排除此干擾，現(xiàn)將整個測量裝置置于箱體內(nèi)。溫度浮動較大時，增加必要的保溫設(shè)備，以保證裝置測量的準(zhǔn)確性。裝置的現(xiàn)場試驗圖如圖2所示。

1—電動隔膜泵；2—空氣壓縮機(jī)；3—測量管；4—差壓傳感器；5—截止閥；6—質(zhì)量流量計；7—離心泵；8—電動三通球；9—脈沖阻尼器。

在測量過程中，每次流量的測量需要保持一定周期，等流量達(dá)到穩(wěn)定之后再進(jìn)行，剪切速率511和1 022 s-1可在一定流量下同時讀出，其他剪切速率需要通過程序控制電動隔膜泵調(diào)節(jié)流量來得到。輸出傳統(tǒng)六速測量結(jié)果，得到塑性黏度、動切力、靜切力、稠度系數(shù)以及流性指數(shù)等值的工藝流程圖，如圖3所示。

圖3 工藝流程圖

3 試驗測試

3.1 測量數(shù)據(jù)精確度試驗

配制質(zhì)量分?jǐn)?shù)5%的膨潤土漿作為鉆井液樣品[19]，將其導(dǎo)入本鉆井液流變性實時測量裝置進(jìn)行試驗測試。測試時，保證在每個排量下的壓差保持穩(wěn)定以減小瞬態(tài)效應(yīng)[20]，然后開始記錄數(shù)據(jù)，可得鉆井液樣品的壓差測量值，如圖4所示。

從圖4可以看出，壓差值存在波動，但均在一定范圍內(nèi)，且符合鉆井液流變性規(guī)律，利用所述科學(xué)方法進(jìn)行壓差測量值數(shù)據(jù)處理，依照處理后壓差數(shù)據(jù)輸出部分黏度值曲線圖。其中流變曲線反映的是非牛頓流體從低剪切速率到高剪切速率的寬切范圍內(nèi)的特性[21]，流動曲線是鉆井液流變性實時測量計算所需的曲線，結(jié)果如圖5、圖6、圖7及圖8所示。

圖4 流量及壓差變化曲線

圖5 表觀黏度20～46 mPa·s時lnγ-lnτw流動曲線

圖6 表觀黏度52～78 mPa·s時lnγ-lnτw流動曲線

由圖7可以看到，表觀黏度為20.0、24.0及38.5 mPa·s時，試驗測量輸出的流變特性曲線和旋轉(zhuǎn)黏度計所測量的曲線幾乎相同，僅在46.0 mPa·s高剪切速率時有些許偏差，這可能是隨機(jī)系統(tǒng)及環(huán)境誤差所導(dǎo)致的。

由圖7與圖8可以看到，表觀黏度為52、62、71及78 mPa·s時，試驗測量輸出的流變特性曲線和旋轉(zhuǎn)黏度計所測量的曲線幾乎相同，能夠驗證試驗裝置測量值的準(zhǔn)確性。

總之，從流變特性曲線可以得出，試驗測量值與范式六速旋轉(zhuǎn)黏度計測量值非常接近，其變化趨勢幾乎相同，且測量誤差在很小的范圍內(nèi)，在有系統(tǒng)及環(huán)境影響的情況下，其誤差值也在4%以內(nèi)，能夠滿足現(xiàn)場自動測量要求，即根據(jù)各個黏度的lnγ-lnτw流動曲線可以實時得出計算所需廣義流動行為指數(shù)n值。從圖5和圖7可以看到，在試驗過程中，由于系統(tǒng)流道的優(yōu)化設(shè)計，n值一直處于比較穩(wěn)定的狀態(tài)，且在本系統(tǒng)中所采取的擬合算法使得其曲線擬合程度極高。

圖7 表觀黏度20～46 mPa·s時流變特性曲線

依據(jù)各黏度的流變曲線，在輸出廣義流動行為指數(shù)n后，便可實時計算出鉆井液的黏度。在試驗中，還使用旋轉(zhuǎn)黏度計進(jìn)行了對比，獲取作為參照的流變曲線及數(shù)據(jù)，將兩組試驗數(shù)據(jù)進(jìn)行對比，結(jié)果如表1所示。

由表1可得，本試驗裝置所測得數(shù)據(jù)與旋轉(zhuǎn)黏度計測量數(shù)據(jù)相差較小，在合理的誤差范圍內(nèi)，能夠替代六速旋轉(zhuǎn)黏度計。

表1 試驗結(jié)果對比

3.2 鉆井液流變性實時測量試驗

配制質(zhì)量分?jǐn)?shù)5%的膨潤土漿，并通過純堿與燒堿的混合物將其pH值調(diào)節(jié)至8～9，在一定時刻開始使用鉆井液流變性實時測量裝置進(jìn)行試驗。

測量過程如下：在原有樣品下自動化測量40 min后，添加重晶石粉等物質(zhì)增加鉆井液密度[22]，繼續(xù)自動化測量40 min后，通過添加PAC-LV、XC和KPAM等聚合物增加鉆井液樣品的黏度[23]，繼續(xù)自動化測量40 min后，再次添加重晶石粉等物質(zhì)增加鉆井液密度，并測量40 min，模擬井下鉆井液變化情況，實時記錄黏度變化趨勢。

在測量開始20 min后，取一定量鉆井液樣品進(jìn)行六速旋轉(zhuǎn)黏度計標(biāo)定測量，之后，每隔40 min進(jìn)行一次取樣測量，結(jié)果如圖9所示。圖9中AV、PV和YP分別表示表觀黏度、塑性黏度和動切力。從圖9可以看到，在添加重晶石粉提高鉆井液密度后，鉆井液流變性自動化測量在初始時受一定程度影響，而后快速趨于平穩(wěn)，且與六速旋轉(zhuǎn)黏度計測量結(jié)果高度重合，其測量結(jié)果誤差為±2 mPa·s。在添加PAC-LV、XC和KPAM等聚合物后，鉆井液流變性有大幅度的變化，本試驗裝置自動化測量結(jié)果依然能夠快速反映出鉆井液流變性的變化差異，并與六速旋轉(zhuǎn)黏度計測量結(jié)果高度重合。當(dāng)表觀黏度為25 mPa·s以下時，測量誤差為±2 mPa·s；當(dāng)表觀黏度為25～150 mPa·s時，測量誤差為±3%。最后再次添加重晶石粉，本試驗裝置依然可以快速準(zhǔn)確測量當(dāng)前鉆井液流變性，誤差在±4%，證明本試驗裝置的流變性測量在鉆井液密度變化的情況下仍然可以正常進(jìn)行。

圖9 實時測量結(jié)果輸出曲線

4 結(jié) 論

(1)針對目前井場所用的鉆井液流變性測量存在因測量時間過長、手工操作誤差及耗費人力等導(dǎo)致鉆井安全風(fēng)險的問題，研制了一套鉆井液流變性實時測量裝置。利用雙測量管、質(zhì)量流量計及電動隔膜泵等器件搭建了鉆井液流變性測量試驗平臺，并通過調(diào)節(jié)不同參數(shù)的鉆井液進(jìn)行了試驗驗證。

(2)對試驗平臺進(jìn)行了數(shù)據(jù)測量精度對比試驗，試驗結(jié)果與現(xiàn)階段的API 2003的鉆井液流變性測量結(jié)果相差無幾，誤差在3%以內(nèi)。

(3)在對鉆井液密度及黏度變化情況下的試驗平臺實時測量中，所得結(jié)果與不間斷采樣旋轉(zhuǎn)黏度計測量數(shù)據(jù)幾乎相同，誤差在4%以內(nèi)。