楊宏偉，李 軍，劉金璐，柳貢慧，高 旭，趙軒剛

（1.中國(guó)石油大學(xué)(北京)石油工程學(xué)院，北京 102249；2.中國(guó)石油長(zhǎng)慶油田分公司第二采氣廠，陜西榆林 719000）

控壓放水技術(shù)通過(guò)將高壓鹽水透鏡體中的鹽水釋放，使井筒周?chē)牡貙訅毫档?，從而擴(kuò)大安全鉆井密度窗口，解決鉆進(jìn)高壓鹽水地層的難題[1]。田徑[2]提出了一種放噴泄壓技術(shù)，為控壓放水技術(shù)應(yīng)用提供了參考。程天輝等人[3]介紹了控壓放水技術(shù)在克深905 井的應(yīng)用，為克深區(qū)塊高壓鹽水層鉆進(jìn)提供了指導(dǎo)。近年來(lái)，該技術(shù)逐漸得到重視，多位學(xué)者進(jìn)行了理論和實(shí)驗(yàn)研究[4-5]。目前來(lái)看，該技術(shù)是一項(xiàng)綜合了壓井和試井理論的新技術(shù)，但現(xiàn)場(chǎng)應(yīng)用大多依靠經(jīng)驗(yàn)完成，還不清楚施工過(guò)程中關(guān)鍵參數(shù)對(duì)控壓放水效果的影響規(guī)律。

為明確關(guān)鍵參數(shù)對(duì)控壓放水效果的影響規(guī)律，筆者根據(jù)控壓放水的技術(shù)特征、工藝特點(diǎn)，總結(jié)了控壓放水工藝流程，將控壓放水施工劃分為4 個(gè)階段，并分階段建立了相關(guān)模型。在此基礎(chǔ)上，建立了參數(shù)動(dòng)態(tài)變化的數(shù)學(xué)模型，對(duì)影響控壓放水效果的關(guān)鍵參數(shù)進(jìn)行了模擬分析，得到了關(guān)井時(shí)間、節(jié)流閥承壓極限等參數(shù)與施工效果的定量關(guān)系。然后，以塔里木油田克深A(yù) 井為例，通過(guò)對(duì)比模擬計(jì)算結(jié)果與實(shí)測(cè)結(jié)果，驗(yàn)證了模型的可靠性。

1 控壓放水工藝流程

控壓放水技術(shù)是一項(xiàng)具有多個(gè)物理過(guò)程的新技術(shù)。整個(gè)控壓放水過(guò)程如圖1所示。根據(jù)工藝特點(diǎn)，可將控壓放水過(guò)程劃分為以下4 個(gè)階段：

圖1 控壓放水工藝流程Fig.1 Pressure-controlled drainage process

1）控制節(jié)流閥，循環(huán)調(diào)整鉆井液密度，每次降低0.02～0.03 kg/L。通過(guò)液面監(jiān)測(cè)，控制套壓保證壓穩(wěn)地層，循環(huán)一周，保證出入口密度均勻，關(guān)井求壓。該階段的目的是增大井底負(fù)壓差，便于后續(xù)施工中將地層中的高壓水順利排出。

2）控制套壓在0～2 MPa，每放水5～10 m3關(guān)井求壓一次，如果立壓較高且套壓沒(méi)有超過(guò)極限套壓，則繼續(xù)放水5～10 m3；如果立壓較高，且套壓已經(jīng)非常接近極限套壓，則節(jié)流循環(huán)排污，調(diào)整鉆井液性能。該階段是施工過(guò)程中最為關(guān)鍵的，也是難度最大的，既要盡可能排出更多地層水，還要保持井口套壓處于安全范圍內(nèi)。

3）重復(fù)階段2）的操作，直至放水后關(guān)井立壓接近于0。

4）重復(fù)階段1）—階段3）的操作，直至將地層壓力當(dāng)量密度降至要求值。

針對(duì)上述工藝流程，建立了對(duì)應(yīng)的數(shù)學(xué)模型。需要說(shuō)明的是，在關(guān)井期間，地層壓力會(huì)不斷恢復(fù)，這在上述工藝流程中沒(méi)有體現(xiàn)出來(lái)。

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

首先，針對(duì)控壓放水不同階段的具體情況，建立了相關(guān)參數(shù)動(dòng)態(tài)變化的數(shù)學(xué)模型；然后，基于地層鹽水滲流理論和井筒流動(dòng)理論，建立了地層壓力動(dòng)態(tài)分布模型、放水量模型、循環(huán)降密度和循環(huán)排污模型；最后，在深入研究關(guān)井期間地層壓力恢復(fù)過(guò)程的基礎(chǔ)上，建立了關(guān)井期間地層壓力恢復(fù)模型。

2.1 地層壓力動(dòng)態(tài)分布模型

控壓放水過(guò)程中，每降低一次鉆井液密度，放出一定量的地層水后，通過(guò)觀察關(guān)井立壓確定井底壓力的變化。因此，控壓放水類(lèi)似于一個(gè)定井底壓力放水過(guò)程。高壓鹽水層透鏡體中的鹽水向井筒內(nèi)流動(dòng)的過(guò)程，可以分為非穩(wěn)定滲流和擬穩(wěn)定滲流2 個(gè)階段[6]。

2.1.1 非穩(wěn)定滲流階段

井底壓力低于地層壓力時(shí)，由于巖石骨架膨脹變形，透鏡體中的鹽水受“彈性驅(qū)動(dòng)”進(jìn)入井筒。非穩(wěn)定控壓放水階段的地層壓力傳播過(guò)程：1）初始關(guān)井時(shí)，地層內(nèi)各點(diǎn)的壓力相同；2）每次鉆井液密度降低0.02～0.03 kg/L，井筒壓力下降，地層流體受巖石骨架變形擠壓流出，近井壁地層中的壓力也開(kāi)始下降，壓力傳播半徑逐漸增大；3）當(dāng)壓力傳播至透鏡體邊界re時(shí)，不穩(wěn)定滲流開(kāi)始向擬穩(wěn)定滲流階段過(guò)渡。

據(jù)此，建立了非穩(wěn)定控壓放水階段的地層鹽水滲流模型，其表達(dá)式為[7]：

式中：p為地層孔隙壓力，MPa；pi為原始地層壓力，MPa；q為高壓鹽水層出水速率，m3/s；μ為地層鹽水黏度，mPa·s；K為地層滲透率，D；h為打開(kāi)高壓鹽水層厚度，m；φ為高壓鹽水層孔隙度；Ct為綜合壓縮系數(shù)，MPa-1；r為徑向距離，m；t為時(shí)間，s；rw為井筒半徑，m。

2.1.2 擬穩(wěn)定滲流階段

當(dāng)?shù)貙訅毫鞑サ酵哥R體邊界時(shí)，因沒(méi)有能量進(jìn)行補(bǔ)充，透鏡體邊界上的壓力開(kāi)始下降。直至某一時(shí)刻放水量基本趨近于0，地層內(nèi)各點(diǎn)壓力不再變化。此后，可以繼續(xù)降低井底壓力，使地層內(nèi)各點(diǎn)的壓力繼續(xù)降低。擬穩(wěn)定滲流階段的地層壓力變化情況如圖2所示。

圖2 擬穩(wěn)定滲流階段地層壓力變化示意Fig.2 Formation pressure changes in quasi-stable seepage stage

基于此，建立了擬穩(wěn)定控壓放水階段的地層鹽水滲流模型，其表達(dá)式為[8]：

式中：pw為井底壓力，MPa。

結(jié)合初始條件和邊界條件，式（1）和式（4）可以采用有限差分法離散，即：

式中：η為導(dǎo)壓系數(shù)，cm2/s。

2.2 放水量計(jì)算模型

根據(jù)達(dá)西定律，建立了各向同性高壓鹽水層中鹽水的滲流速度模型，其表達(dá)式為[9]：

式中，vsw為地層鹽水滲流速度，m/s。

根據(jù)高壓鹽水層的孔隙壓力隨時(shí)間和空間的分布，對(duì)壓降進(jìn)行有限差分得：

將式（9）代入式（8），得到鹽水滲流速度：

2.3 循環(huán)降密度模型

在循環(huán)降密度階段，使用低密度鉆井液替換原來(lái)的高密度鉆井液；在循環(huán)排污階段，使用排污鉆井液將受鹽水污染的鉆井液循環(huán)排出。因此，降密度和循環(huán)排污過(guò)程都類(lèi)似于壓井過(guò)程，包括2 個(gè)階段[10-11]：1）注入的流體由井口進(jìn)入鉆柱，后向下流至鉆頭；2）注入的流體由環(huán)空底部上返至井口，替換出原來(lái)的流體。其中，循環(huán)排污過(guò)程如圖3所示。

圖3 循環(huán)排污過(guò)程示意Fig.3 Cyclic sewage disposal process

第一階段，套壓保持為關(guān)井套壓不變，根據(jù)以下公式計(jì)算立壓的變化情況：

式中：pt為立壓，MPa；psp為關(guān)井立壓，MPa；vpi為鉆柱內(nèi)流體流速，m/s；t1為注入流體到達(dá)鉆頭所需要的時(shí)間，s；ρk為注入流體密度，kg/m3；ρm為原井筒流體密度，kg/m3；Δpm為摩阻，MPa。

第二階段，立壓保持不變，根據(jù)以下公式計(jì)算套壓的變化情況：

式中：pa為套壓，MPa；pap為關(guān)井套壓，MPa；vo為環(huán)空流體流速，m/s；tt為循環(huán)一周的時(shí)間，s。

2.4 關(guān)井期間地層壓力恢復(fù)方程

由疊加原理可得，若在tp時(shí)刻關(guān)井，則tp至tp+Δt時(shí)間段地層壓力的恢復(fù)可用Horner 方程描述[12]，其表達(dá)式為：

式中：pws（Δt）為關(guān)井Δt時(shí)段后的地層壓力，MPa；pws（Δt=0）為關(guān)井時(shí)的地層壓力，MPa；qw為地層產(chǎn)水量，m3/d；μ為地層流體黏度，mPa·s；K為高壓鹽水層滲透率，D；h為打開(kāi)鹽水層厚度，m。

從式（13）可以看出，在半對(duì)數(shù)坐標(biāo)系中，地層壓力恢復(fù)值pws與lg[(tp+Δt)/Δt]呈線性關(guān)系，該直線即為Horner 直線，其中直線的斜率表達(dá)式為：

考慮表皮因子S，當(dāng)關(guān)井時(shí)間為1 h 時(shí)，地層壓力恢復(fù)值為：

3 計(jì)算分析

3.1 計(jì)算參數(shù)

克深地區(qū)鹽膏層地質(zhì)特征復(fù)雜，層間普遍發(fā)育高壓或超高壓鹽水層[13-19]，其主要具有如下特征：不同區(qū)塊的高壓鹽水層厚度分布不均，縱橫向上分布規(guī)律性差；各個(gè)區(qū)塊的高壓鹽水層壓力梯度變化大，且當(dāng)量密度高達(dá)2.40～2.60 kg/L；鹽膏層間含有抗壓強(qiáng)度較低的白云巖夾層，采用高密度鉆井液鉆進(jìn)時(shí)極易發(fā)生漏失或溢漏同存，井筒壓力控制難度極大[20-21]。

克深A(yù) 井鉆至井深7 200.74 m 時(shí)發(fā)現(xiàn)溢流，判斷該處存在高壓鹽水層。該段地層密度窗口窄，若采用提高鉆井液密度的方法處理溢流，極可能造成地層漏失，因此采用控壓放水技術(shù)處理該井段?？松預(yù) 井的基礎(chǔ)參數(shù)為：滲透率0.09 mD，鹽水黏度0.5 mPa·s，打開(kāi)鹽水層厚度20 m，初始地層壓力系數(shù)2.58，導(dǎo)壓系數(shù)11.4 cm2/s，鉆桿外徑127.0 mm，鉆桿內(nèi)徑89.1 mm，套管長(zhǎng)度7 086 m，關(guān)井套壓8.1 MPa，關(guān)井立壓0.5 MPa。

3.2 計(jì)算對(duì)比及誤差分析

模擬計(jì)算控壓放水15 d 的情況，地層壓力當(dāng)量密度由2.63 kg/L 降至2.48 kg/L?？貕悍潘M及修正結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果的對(duì)比如圖4所示。

圖4 模擬及修正結(jié)果與實(shí)測(cè)結(jié)果對(duì)比Fig.4 Comparison of simulation and correction results with measured results

由圖4可知，a 區(qū)域（前7 次）放水的模擬計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)結(jié)果一致性較好，但從第8 次放水開(kāi)始，模擬結(jié)果與實(shí)測(cè)結(jié)果的偏差逐漸增大，這主要是因?yàn)榈? 次控壓放水速率突然增大（見(jiàn)圖5）。根據(jù)放水速率變化情況對(duì)模擬結(jié)果進(jìn)行了修正，修正后的結(jié)果與實(shí)測(cè)結(jié)果基本吻合。

圖5 克深A(yù) 井放水速率變化曲線Fig.5 Variation curve of water drainage rates in Well Keshen A

3.3 關(guān)井時(shí)間分析

分析了不同滲透率條件下關(guān)井時(shí)間分別為1，3 和5 h 時(shí)的地層壓力，結(jié)果見(jiàn)圖6。

由圖6可以發(fā)現(xiàn)，關(guān)井時(shí)間與地層壓力之間存在2 種規(guī)律：1）關(guān)井時(shí)間越短，地層壓力下降越快，這是因?yàn)?，關(guān)井時(shí)間越短，地層壓力恢復(fù)的時(shí)間越短，地層壓力相對(duì)下降幅度越大；2）地層滲透率越低，關(guān)井時(shí)間對(duì)地層壓力的降低效果影響越顯著，這一方面是因?yàn)殛P(guān)井時(shí)間縮短，地層壓力恢復(fù)緩慢，另一方面是因?yàn)殡S著關(guān)井時(shí)間縮短，有效放水時(shí)間隨之增長(zhǎng)。

圖6 不同滲透率條件下關(guān)井時(shí)間分別為1,3 和5 h 時(shí)的地層壓力Fig.6 Formation pressure at shut-in time of 1 hr,3 hrs and 5 hrs under different permeability conditions

3.4 不同承壓極限節(jié)流閥分析

前6 次控壓放水過(guò)程中套壓和放水量的變化曲線如圖7所示。

由圖7可知，節(jié)流閥的極限套壓是限制單次放水量的關(guān)鍵因素。隨著累計(jì)放水量增加，關(guān)井套壓隨之升高。當(dāng)關(guān)井套壓接近節(jié)流閥承壓極限時(shí)，需結(jié)束該次放水操作，進(jìn)行循環(huán)排污。循環(huán)排污越頻繁，整個(gè)控壓放水周期越長(zhǎng)。

進(jìn)行了不同承壓能力節(jié)流閥單次放水量對(duì)比分析，結(jié)果見(jiàn)圖8（圖中：pc為節(jié)流閥承壓極限，MPa）。

圖8 不同承壓能力節(jié)流閥單次放水量對(duì)比Fig.8 Comparison of single drainage of throttle with different pressure-bearing capacity

由圖8可知，節(jié)流閥承壓極限越大，單次放水釋放的地層鹽水量越多，若將節(jié)流閥承壓極限從5 MPa提高至15 MPa，循環(huán)排污次數(shù)可以減少一半，這將大大縮短放水周期。因此，應(yīng)盡可能提高節(jié)流閥承壓極限，一方面可以延長(zhǎng)放水時(shí)間，減少循環(huán)排污次數(shù)，從而縮短放水周期；另一方面放水時(shí)間延長(zhǎng)，關(guān)井時(shí)間相對(duì)會(huì)縮短，地層壓力恢復(fù)也會(huì)減緩，地層壓力下降會(huì)更快，即控壓放水效果更好。

3.5 地層滲透率影響分析

考慮關(guān)井期間地層壓力的恢復(fù)過(guò)程，模擬了不同滲透率條件下地層壓力當(dāng)量密度隨放水次數(shù)的變化曲線，結(jié)果如圖9所示。

由圖9可知，當(dāng)?shù)貙訚B透率大于0.5 D 時(shí)，地層壓力當(dāng)量密度隨放水次數(shù)增加幾乎無(wú)變化；當(dāng)?shù)貙訚B透率較低（K=0.01，0.05，0.10 和0.30 D）時(shí)，地層壓力當(dāng)量密度隨放水次數(shù)的變化趨勢(shì)基本一致，前7 次放水時(shí)地層壓力當(dāng)量密度下降顯著，7～10 次的下降速度減緩，而10 次之后下降速度顯著降低；同時(shí)，地層滲透率越低，最終地層壓力當(dāng)量密度下降幅度也越大。這主要是因?yàn)椋貙訚B透率越高，滲流阻力越小，井筒周?chē)貙訅毫謴?fù)速率越快，所以在高滲透地層實(shí)施控壓放水的效果并不明顯。此外，隨著地層壓力降低，地層滲流的壓差也加大，使井底壓力得到補(bǔ)充，從而表現(xiàn)為地層壓力當(dāng)量密度降低速度減緩。實(shí)際控壓放水過(guò)程中，可以對(duì)高壓鹽水層進(jìn)行7 次試放水，若此期間地層壓力下降顯著，表明控壓放水技術(shù)在該地層的適用性較好。

圖9 不同滲透率下地層壓力當(dāng)量密度隨放水次數(shù)的變化曲線Fig.9 Variation curves of equivalent density of formation pressure with times of water drainage under different permeability conditions

4 結(jié)論與建議

1）分析了控壓放水的技術(shù)特征、工藝特點(diǎn)，總結(jié)了控壓放水的工藝流程。在此基礎(chǔ)上，建立了參數(shù)動(dòng)態(tài)變化的數(shù)學(xué)模型，來(lái)模擬控壓放水全過(guò)程。經(jīng)計(jì)算驗(yàn)證，地層壓力變化模擬結(jié)果與實(shí)測(cè)結(jié)果基本一致。

2）模擬研究可得，關(guān)井時(shí)間是影響控壓放水效果的重要因素，適當(dāng)縮短關(guān)井時(shí)間，可以增大地層壓力降低幅度；對(duì)于低滲透率地層，關(guān)井時(shí)間的影響更加顯著；提高節(jié)流閥的承壓能力，不僅可以縮短控壓放水周期，控壓放水效果也更好。

3）實(shí)際控壓放水過(guò)程中，建議對(duì)高壓鹽水層進(jìn)行7 次試放水，若此期間地層壓力下降明顯，表明控壓放水技術(shù)在該地層的適用性較好。此方法可以作為該控壓放水工藝的適應(yīng)性評(píng)價(jià)方法。