流-固-化耦合條件下硬脆性泥頁(yè)巖井壁漸進(jìn)破壞效應(yīng)探討

肖志強(qiáng)，賈善坡,2，亓憲寅,，戴永浩，呂方，賈陸鋒，溫曹軒

(1.長(zhǎng)江大學(xué)巖土力學(xué)與工程研究中心，湖北荊州，434023；2.東北石油大學(xué)非常規(guī)油氣研究院，黑龍江大慶，163318；3.中國(guó)科學(xué)院武漢巖土力學(xué)研究所巖土力學(xué)與工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北武漢，430071；4.湖北省油田化學(xué)產(chǎn)業(yè)技術(shù)研究院，湖北荊州，434000)

中深層勘探開(kāi)發(fā)是東部油田近年勘探開(kāi)發(fā)的重點(diǎn)，因其埋藏深，巖性水敏性強(qiáng)，裸眼浸泡時(shí)間長(zhǎng)，同時(shí)存在裂縫，加劇了井壁失穩(wěn)，使中深層硬脆性泥頁(yè)巖的井壁穩(wěn)定成為制約鉆速的瓶頸[1-3]。因此，從降低鉆井成本以及加快油氣田勘探開(kāi)發(fā)的角度出發(fā)，對(duì)泥頁(yè)巖井壁穩(wěn)定性問(wèn)題進(jìn)行深入的研究具有非常重要的意義。NP 凹陷屬于典型的陸相淡水湖泊沉積盆地，形成了半深湖-深湖相暗色泥頁(yè)巖和油頁(yè)巖等優(yōu)質(zhì)烴源巖，主力烴源巖為沙河街組巖性圈閉，埋藏深度超過(guò)4 000 m，發(fā)育著幾百米的硬脆性泥頁(yè)巖。沙河街組硬脆性泥頁(yè)巖黏土礦物含量較高、微裂縫發(fā)育，在實(shí)鉆過(guò)程中，壓差滲透、裂縫自吸水及長(zhǎng)時(shí)間浸泡引發(fā)的泥頁(yè)巖水化效應(yīng)顯著，剝落掉塊，導(dǎo)致井徑變化較大，經(jīng)常造成鉆具阻卡、電測(cè)遇阻等井下復(fù)雜情況，嚴(yán)重影響鉆井效率。地層原有的力學(xué)平衡及化學(xué)平衡，在井眼鉆開(kāi)后便被打破，由于液壓梯度和化學(xué)梯度的共同作用，導(dǎo)致了井壁圍巖浸泡區(qū)的出現(xiàn)，滲流擴(kuò)散效應(yīng)和鉆井液對(duì)泥頁(yè)巖的水化作用，使得浸泡區(qū)圍巖成為了隨著時(shí)間而變化的非均質(zhì)體，這一系列的影響導(dǎo)致了泥頁(yè)巖井壁經(jīng)歷了一個(gè)非常復(fù)雜的過(guò)程。針對(duì)泥頁(yè)巖井壁穩(wěn)定性的研究，國(guó)內(nèi)外學(xué)者已經(jīng)開(kāi)展過(guò)很多工作，但由于研究方向的不一致以及問(wèn)題的高度復(fù)雜性，導(dǎo)致了研究成果的不一致，并且每一種理論都有其一定的局限性。YEW 等[4]將泥頁(yè)巖水化膨脹應(yīng)力比擬為熱膨脹變溫應(yīng)力，將水向頁(yè)巖中的運(yùn)動(dòng)比擬成熱擴(kuò)散。根據(jù)熱彈性力學(xué)理論，對(duì)泥頁(yè)巖水化產(chǎn)生的力學(xué)效應(yīng)建立了定量化模型。該理論模型考慮了含水量的變化對(duì)泥頁(yè)巖井壁穩(wěn)定的影響效應(yīng)，但是忽略了泥頁(yè)巖與鉆井液之間化學(xué)作用的本質(zhì)；張樂(lè)文等[5-6]應(yīng)用半透膜等效孔隙壓力理論研究了泥頁(yè)巖與水基鉆井液之間的相互作用，建立了泥頁(yè)巖地層井壁穩(wěn)定的力化耦合模型，但是該研究未考慮離子擴(kuò)散和離子交換對(duì)泥頁(yè)巖水化的影響；HEIDUG 等[7-9]建立的泥頁(yè)巖吸水膨脹耦合模型，該模型較全面地研究了井壁穩(wěn)定的力學(xué)-化學(xué)耦合效應(yīng)，但沒(méi)有考慮到流體孔隙壓力與巖石變形之間的聯(lián)合作用；CHEE 等[10-13]在綜合考慮泥頁(yè)巖-鉆井液體系中的電化勢(shì)、離子運(yùn)移和流體流動(dòng)、溶質(zhì)擴(kuò)散過(guò)程的非線(xiàn)性以及流體流動(dòng)和離子運(yùn)移對(duì)固體變形影響的基礎(chǔ)上，建立了泥頁(yè)巖井壁穩(wěn)定流-固-化耦合模型，但是該模型中存在著很多難以確定的參數(shù)，不方便實(shí)際應(yīng)用，同時(shí)也未考慮到地層鉆井卸載以及塑性變形。泥頁(yè)巖地層一經(jīng)鉆開(kāi)，便打破了其原有的平衡狀態(tài)，井壁圍巖參數(shù)會(huì)隨著與鉆井液接觸的時(shí)間變化而發(fā)生變化。因此，一個(gè)有效的井眼預(yù)測(cè)模型應(yīng)該使發(fā)生在泥頁(yè)巖地層中的力學(xué)以及物理化學(xué)動(dòng)態(tài)過(guò)程得到體現(xiàn)。NP凹陷W號(hào)構(gòu)造為巖性圈閉，其埋藏深度為3 800～4 500 m。東二段至沙河街井段中鉆遇到層理和微裂隙都比較發(fā)育的大段硬脆性泥頁(yè)巖，在鉆井過(guò)程中表現(xiàn)為易塌、易漏等復(fù)雜情況，導(dǎo)致該地區(qū)的鉆井效率被嚴(yán)重制約。本文作者總結(jié)前人的研究成果，并對(duì)硬脆性泥頁(yè)巖理化性能、微觀結(jié)構(gòu)以及力學(xué)特性進(jìn)行分析，解釋其失穩(wěn)機(jī)理。為了探索井眼破壞的演化規(guī)律，通過(guò)考慮鉆井過(guò)程中巖石的實(shí)際卸載過(guò)程、強(qiáng)度弱化以及塑性變形的特性，建立了泥頁(yè)巖井眼漸進(jìn)破壞耦合模型，從機(jī)理上提高了保持硬脆性泥頁(yè)巖井壁穩(wěn)定的認(rèn)識(shí)，為鉆井液的優(yōu)選提供參考依據(jù)。

1 泥頁(yè)巖組成及理化性能試驗(yàn)

在實(shí)鉆過(guò)程中，沙河街組泥頁(yè)巖井徑變化較大，井壁失穩(wěn)問(wèn)題突出。因此，需要對(duì)泥頁(yè)巖的巖性、礦物組成、微觀結(jié)構(gòu)及水化等方面開(kāi)展測(cè)試研究。

1.1 地層礦物組成

利用X線(xiàn)衍射儀，對(duì)W號(hào)構(gòu)造3口井沙河街組泥頁(yè)巖段巖芯進(jìn)行全巖礦物和黏土礦物分析，全巖礦物和黏土礦物組成與質(zhì)量分?jǐn)?shù)如表1和表2所示。

從表1和表2可以看出：沙河街組泥頁(yè)巖的礦物主要為黏土和石英，另外不同程度地發(fā)育著長(zhǎng)石、方解石、白云石等，脆性礦物(石英、長(zhǎng)石和方解石)較為發(fā)育，其中石英質(zhì)量分?jǐn)?shù)為6.71%～39.94%；黏土含量整體較高，其質(zhì)量分?jǐn)?shù)為19.09%～43.06%，黏土礦物中較為發(fā)育的主要為伊利石以及伊/蒙混層，沙河街組礦物中均未見(jiàn)蒙脫石的發(fā)育，其中伊利石質(zhì)量分?jǐn)?shù)為34.47%～56.64%，伊/蒙混層質(zhì)量分?jǐn)?shù)為6.72%～39.74%。

1.2 泥頁(yè)巖細(xì)觀構(gòu)造

探究黏土礦物晶體的定向排列、膠結(jié)結(jié)構(gòu)以及微裂隙的發(fā)育和分布狀況要從泥頁(yè)巖的微觀結(jié)構(gòu)入手。目前研究巖石內(nèi)部微裂縫等微觀結(jié)構(gòu)普遍所采用的有效手段是利用環(huán)境掃描電鏡進(jìn)行觀察。利用掃描電鏡分析W 號(hào)構(gòu)造沙河街組泥頁(yè)巖微觀構(gòu)造和黏土礦物賦存形態(tài)，掃描結(jié)果如圖1所示。

從圖1可以看出：W號(hào)構(gòu)造沙河街組泥頁(yè)巖內(nèi)部微裂縫、微孔洞等發(fā)育程度較高，微裂縫以及微孔洞的發(fā)育一方面破壞了泥頁(yè)巖的結(jié)構(gòu)完整性，弱化原巖的力學(xué)性能，另一方面也為鉆井流體侵入地層提供了流動(dòng)通道。鉆井流體通過(guò)這樣的通道侵入地層不僅可能誘發(fā)水力劈裂作用促使井眼加劇破壞，更重要的是與地層中的黏土礦物以及有機(jī)質(zhì)等發(fā)生作用，造成井壁坍塌等事故。

1.3 鉆井液體作用對(duì)泥頁(yè)巖結(jié)構(gòu)的影響

泥頁(yè)巖中的黏土礦物在遇到水的情況下，會(huì)發(fā)生水化膨脹的現(xiàn)象，沙河街組泥頁(yè)巖黏土礦物中以伊利石以及伊/蒙混層較為發(fā)育，蒙脫石晶層間的連接力很弱，因此，能在其晶層間以及顆粒外表面進(jìn)行水化作用，水化膨脹性很強(qiáng)，而伊利石晶層間的連接較緊密，水化膨脹僅能發(fā)生在顆粒外表面，水化膨脹性相對(duì)較弱，各種黏土礦物的水化膨脹能力由大到小排序?yàn)椋好擅撌?，?蒙混層，伊利石，高嶺石，綠泥石?？梢?jiàn)：不同黏土礦物的水化膨脹能力不同，以至于不同的黏土礦物遇水后吸水膨脹速率相差很大，進(jìn)而產(chǎn)生相差很大的膨脹壓力；再者鉆井流體會(huì)使微裂縫尖端產(chǎn)生應(yīng)力集中，這些因素的存在很可能導(dǎo)致泥頁(yè)巖的最終破壞。

圖1 W-82井泥頁(yè)巖巖芯掃描電鏡圖像Fig.1 SEM images of shale core in well W-82

采用高倍偏光顯微鏡觀察泥頁(yè)巖在流體作用后細(xì)觀結(jié)構(gòu)的變化。圖2和圖3所示分別為清水及水基鉆井液作用下，高倍偏光顯微鏡下觀察到的巖樣結(jié)構(gòu)變化及微裂縫的形態(tài)特征。可見(jiàn)：清水、鉆井液侵入前的巖樣結(jié)構(gòu)較為致密、完整，而經(jīng)過(guò)清水、鉆井液一段時(shí)間的作用后，泥頁(yè)巖表面發(fā)生了變化，這是泥頁(yè)巖內(nèi)部微裂紋的萌生、擴(kuò)展以及分叉所導(dǎo)致的結(jié)果，從而進(jìn)一步破壞了原巖的完整性以及弱化原巖的力學(xué)性能，并且清水、鉆井液還將沿微裂縫侵入到泥頁(yè)巖的內(nèi)部，與泥頁(yè)巖發(fā)生水化作用進(jìn)一步破壞泥頁(yè)巖的結(jié)構(gòu)，加劇井壁的破壞。

1.4 陽(yáng)離子交換容量

陽(yáng)離子交換即陽(yáng)離子交換性吸附，是泥質(zhì)巖重要特征之一，陽(yáng)離子交換容量是指黏土礦物的表面可用于交換陽(yáng)離子的數(shù)量。泥頁(yè)巖的陽(yáng)離子交換容量是全部黏土礦物陽(yáng)離子交換的總和。陽(yáng)離子的交換容量越大，表示黏土能夠攜帶的負(fù)電荷越多，越能從周?chē)慕橘|(zhì)中吸附更多的陽(yáng)離子來(lái)中和多余的負(fù)電荷，即陽(yáng)離子交換容量越大越有利于泥頁(yè)巖表面水化。因此，可以利用陽(yáng)離子的交換容量，來(lái)預(yù)測(cè)地層潛在的水敏性。黏土在水中分散時(shí)，吸附的陽(yáng)離子要從黏土表面擴(kuò)散，泥頁(yè)巖地層中黏土礦物的含量較大，通過(guò)測(cè)定泥頁(yè)巖的陽(yáng)離子交換容量(CEC)，可以反映出泥頁(yè)巖的水化、膨脹和分散能力，是評(píng)價(jià)地層水化膨脹能力的重要依據(jù)。

對(duì)泥頁(yè)巖巖樣處理后測(cè)定陽(yáng)離子交換容量，W 號(hào)構(gòu)造沙河街組泥頁(yè)巖的陽(yáng)離子交換容量分布范圍為90～235 mmol/kg，平均為146.25 mmol/kg，表明沙河街組泥頁(yè)巖具有一定的水化能力，水基鉆井液作用下易發(fā)生井下失穩(wěn)。

圖2 清水沿微裂縫侵入Fig.2 Clear water intrusion along micro-cracks in shale

圖3 鉆井濾液作用下的微裂縫形態(tài)變化Fig.3 Morphological changes of micro-cracks under action of drilling filtrate

2 鉆井液對(duì)泥頁(yè)巖力學(xué)性質(zhì)影響試驗(yàn)

泥頁(yè)巖的力學(xué)特性是影響井眼穩(wěn)定性的重要因素。由于硬脆性泥頁(yè)巖的水化特性顯著，在對(duì)其進(jìn)行鉆井時(shí)，對(duì)鉆井液具有很高的要求，鉆井液的優(yōu)選對(duì)保持硬脆性泥頁(yè)巖地層的井壁穩(wěn)定起著決定性的作用

2.1 原巖力學(xué)特性

考慮到泥頁(yè)巖的水敏性，針對(duì)鉆井巖芯進(jìn)行干冰冷凍，采用液氮冷卻鉆頭，鉆取試驗(yàn)所用巖樣。取自W-82井沙河街組泥頁(yè)巖巖芯進(jìn)行原巖三軸壓縮力學(xué)特性測(cè)試，三軸應(yīng)力應(yīng)變曲線(xiàn)如圖4所示，泥頁(yè)巖彈性模量為18.27～23.59 GPa，泊松比為0.12～0.23，內(nèi)聚力為24.14 MPa，內(nèi)摩擦角為21.7°。另外，對(duì)5個(gè)沙河街組巖樣進(jìn)行了巴西劈裂試驗(yàn)，抗張強(qiáng)度為1.13～6.09 MPa，平均值為4.41 MPa。

脆性是巖石的一種固有力學(xué)特性，并受巖石組分、結(jié)構(gòu)以及圍壓、溫度等因素的影響、控制。目前，巖石礦物組分及彈性參數(shù)是石油工程領(lǐng)域應(yīng)用較為廣泛的2種巖石脆性評(píng)價(jià)指標(biāo)：基于巖石礦物組分的脆性評(píng)價(jià)和基于巖石彈性參數(shù)的脆性評(píng)價(jià)。參考賈善坡等[14-15]的研究成果，基于脆性礦物的泥頁(yè)巖脆性評(píng)價(jià)結(jié)果為：W-82 井沙河街組泥頁(yè)巖脆性指數(shù)為0.664，W-96井沙河街組泥頁(yè)巖為0.809，W-81井沙河街組泥頁(yè)巖為0.569。

圖4 W-82井泥頁(yè)巖原巖三軸試驗(yàn)結(jié)果(深度為4 150 m)Fig.4 Triaxial test results of original rock in W-82 well(depth of 4 150 m)

RICKMAN等[16]在對(duì)北美地區(qū)泥頁(yè)巖的彈性模量、泊松比進(jìn)行統(tǒng)計(jì)的基礎(chǔ)上，考慮彈性模量越大，泊松比越小，則脆性越大，提出了利用彈性模量和泊松比2個(gè)參數(shù)表征巖石的脆性強(qiáng)弱。利用上述泥頁(yè)巖的三軸壓縮測(cè)試分析結(jié)果：巖芯1 號(hào)(圍壓為0 MPa)對(duì)應(yīng)的脆性指數(shù)為0.424，巖芯2號(hào)(圍壓為25 MPa)為0.367，巖芯3 號(hào)(圍壓為50 MPa)為0.326。表明該區(qū)的泥頁(yè)巖呈現(xiàn)出較高的脆性。

2.2 鉆井液體作用對(duì)泥頁(yè)巖力學(xué)性質(zhì)的影響

為探究泥頁(yè)巖力學(xué)性質(zhì)在鉆井液作用下的變化規(guī)律，為鉆井液保持巖石強(qiáng)度性能的評(píng)價(jià)提供依據(jù)，開(kāi)展了工區(qū)在用鉆井液作用下巖石力學(xué)特性測(cè)試。

由于泥頁(yè)巖鉆井巖芯獲取困難，且用于開(kāi)展三軸試驗(yàn)的有效巖芯非常有限，采用史氏壓入硬度法測(cè)試鉆井液對(duì)泥頁(yè)巖壓入硬度影響。硬度反映了巖石抵抗工具侵入破壞的能力，壓入硬度試驗(yàn)對(duì)巖樣形狀的要求較直剪試驗(yàn)和三軸壓縮試驗(yàn)低，便于大量進(jìn)行測(cè)試分析。

將泥頁(yè)巖浸泡在鉆井液中，通過(guò)測(cè)試不同浸泡時(shí)間下泥頁(yè)巖的壓入硬度，描述泥頁(yè)巖強(qiáng)度的變化特征。選取同一位置泥頁(yè)巖分別開(kāi)展原巖、浸泡鉆井液6 h 和12 h 的壓入硬度試驗(yàn)，共開(kāi)展4組，試驗(yàn)結(jié)果如表3所示。從表3可以看出：隨著浸泡時(shí)間的增加巖樣的壓入硬度逐漸降低，并且浸泡時(shí)間越長(zhǎng)壓入硬度降低越明顯，這些規(guī)律表明了鉆井液對(duì)其強(qiáng)度的影響不可忽略。

表3 壓入硬度隨著鉆井液浸泡時(shí)間的變化Table 3 Change of indentation hardness with drilling fluid immersion time

由于該工區(qū)巖芯短缺，對(duì)鄰區(qū)同層位泥頁(yè)巖進(jìn)行了測(cè)試分析，結(jié)果表明：在對(duì)巖樣進(jìn)行浸泡的過(guò)程中，沙河街組樣品的完整性被破壞，發(fā)現(xiàn)有宏觀裂縫產(chǎn)生，部分巖樣破碎成碎塊，水化程度較為嚴(yán)重。分別進(jìn)行泥頁(yè)巖在鉆井液中浸泡0 h(原巖)，6 h 和12 h 后的單軸和三軸壓縮試驗(yàn)，試驗(yàn)結(jié)果如表4所示。從表4可知：在同一圍壓下，隨著硬脆性泥頁(yè)巖在鉆井液中浸泡時(shí)間的增加，其力學(xué)強(qiáng)度逐漸降低，浸泡6 h 后強(qiáng)度平均降低16.32%，浸泡12 h后強(qiáng)度平均降低23.10%。

表4 鉆井液浸泡作用下泥頁(yè)巖強(qiáng)度變化Table 4 Shale strength changes under drilling fluid immersion

3 泥頁(yè)巖井眼破壞流-固-化耦合模型

綜合上述對(duì)泥頁(yè)巖理化性能和力學(xué)特性試驗(yàn)分析結(jié)果可知，導(dǎo)致該區(qū)鉆井過(guò)程中井壁失穩(wěn)的主要因素有以下幾點(diǎn)：一是泥頁(yè)巖黏土含量較高且黏土礦物中伊/蒙混層較為發(fā)育，作為膨脹性與非膨脹性黏土相間的礦物，較易吸水致使非均勻性水化膨脹，類(lèi)似的規(guī)律通過(guò)陽(yáng)離子交換試驗(yàn)再次證明，該區(qū)的泥頁(yè)巖易與外來(lái)流體發(fā)生水化作用，減弱泥頁(yè)巖結(jié)構(gòu)強(qiáng)度。二是沙河街組泥頁(yè)巖內(nèi)部微裂隙和微孔洞較為發(fā)育，提供了外來(lái)流體侵入泥頁(yè)巖內(nèi)部的通道以及水化作用的空間，造成巖石強(qiáng)度降低，增加井壁坍塌風(fēng)險(xiǎn)。三是W 號(hào)構(gòu)造沙河街組泥頁(yè)巖脆性相對(duì)較強(qiáng)，在鉆井過(guò)程中的應(yīng)力釋放(卸荷)、井筒壓力波動(dòng)、機(jī)械鉆速低或浸泡時(shí)間長(zhǎng)更易導(dǎo)致井壁巖石破裂，為鉆井液侵入提供通道，誘發(fā)、加劇井壁垮塌。因此，研究泥頁(yè)巖井眼漸進(jìn)破壞動(dòng)態(tài)演化時(shí)須考慮鉆井液在地層中的擴(kuò)散與滲透以及巖石的強(qiáng)度弱化行為[17]。

3.1 吸水?dāng)U散模型

泥頁(yè)巖在井眼鉆開(kāi)后便與鉆井液接觸，受到水力梯度和化學(xué)勢(shì)梯度的驅(qū)動(dòng)，引起了水和離子的相互傳遞。泥頁(yè)巖的吸水過(guò)程十分復(fù)雜，這里令ω表示吸附水量的質(zhì)量分?jǐn)?shù)，它隨著時(shí)間和空間而發(fā)生變化，依據(jù)質(zhì)量守恒方程建立了水分的擴(kuò)散方程[18]：

式中：Cf為泥頁(yè)巖吸水?dāng)U散系數(shù)。吸水?dāng)U散系數(shù)表征了地層水變化傳遞的速度，可通過(guò)泥頁(yè)巖吸水?dāng)U散試驗(yàn)測(cè)得。相應(yīng)的邊界和初始條件按照遠(yuǎn)場(chǎng)以及井壁的含水量來(lái)確定。

3.2 壓力傳遞模型

考慮到泥頁(yè)巖的半透膜特性，化學(xué)勢(shì)差的存在會(huì)促使鉆井液中的一些離子發(fā)生運(yùn)移，導(dǎo)致了地層的水活度發(fā)生改變。活度指的是鹽溶液和純水的逸度比，表征了溶液中化學(xué)勢(shì)的強(qiáng)弱，鉆井過(guò)程中地層流體的活度計(jì)算模型為[19]

式中：D為活度擴(kuò)散系數(shù)。活度擴(kuò)散系數(shù)表征了地層水活度傳遞的快慢，可以通過(guò)頁(yè)巖的壓力傳遞試驗(yàn)數(shù)據(jù)擬合得到。

地層的孔隙壓力在鉆井液進(jìn)入地層后，受到地層流體與鉆井液之間的化學(xué)勢(shì)差和水力壓差的共同作用而重新分布。在綜合考慮了化學(xué)勢(shì)差和流體的流動(dòng)之后，建立了孔隙壓力的計(jì)算模型[20-21]：

式中：k為泥頁(yè)巖的滲透率；υ為流體黏度，mPa ?s；C為流體的壓縮系數(shù)，Pa-1；φ為孔隙度；p為地層的流體壓力，MPa；Im為膜效率；R為氣體常數(shù)，8.314 m3?Pa/(K ?mol)；T為熱力學(xué)溫度，K；V為水的偏摩爾體積，1.80×10-5m3/mol；θ為鉆井過(guò)程中地層流體活度；θshale為泥頁(yè)巖中地層水活度。

3.3 井眼破壞模型

有效應(yīng)力可根據(jù)Biot 有效應(yīng)力原理，表示如下：

式 中：σij為總應(yīng)力張量；α為Biot 系數(shù)；δij為Kronecker符號(hào)。

根據(jù)已建立的壓力傳遞模型并結(jié)合有效應(yīng)力原理以及彈塑性力學(xué)理論[22]，建立綜合考慮化學(xué)勢(shì)變化、流體流動(dòng)與骨架變形耦合的本構(gòu)方程，表示如下：

式中：De為彈性剛度張量；ε為總應(yīng)變張量；εp為塑性應(yīng)變張量；εij=(ui,j+uj,i)/2，dεij= dui,j；ui為巖石骨架的位移分量。

考慮到泥頁(yè)巖的實(shí)際抗拉能力，井眼的破壞準(zhǔn)則采用的是修正的Drucker-Prager 準(zhǔn)則(簡(jiǎn)稱(chēng)D-P準(zhǔn)則)。針對(duì)于拉剪型的應(yīng)力狀態(tài)，常規(guī)的D-P 準(zhǔn)則不適用，它只是考慮了壓剪型的應(yīng)力狀態(tài)，而在拉剪區(qū)將會(huì)失去實(shí)驗(yàn)基礎(chǔ)，高估了巖石的抗拉能力。當(dāng)采用過(guò)平衡鉆井方式時(shí)，由于井底的壓力高于地層壓力，在井眼處可能會(huì)有拉剪型應(yīng)力狀態(tài)的出現(xiàn)，故而有必要對(duì)常規(guī)D-P 準(zhǔn)則進(jìn)行修正，修正后的D-P準(zhǔn)則示意圖如圖5所示?；谫Z善坡等[23]的修正思想，采用Rankine 最大拉應(yīng)力準(zhǔn)則對(duì)常規(guī)的D-P準(zhǔn)則進(jìn)行修正，屈服函數(shù)表達(dá)式為

式中：q為Mises 應(yīng)力；σm=(σ′1+σ′2+σ′3)/3 為有效圍壓應(yīng)力；c為巖石黏聚力；?為巖石內(nèi)摩擦角；l0=c0-σt，c0為原巖黏聚力，σt為巖石抗拉強(qiáng)度。

圖5 修正的Drucker-Prager準(zhǔn)則Fig.5 Modified Drucker-Prager criteria

從式(6)可知：若l0= 0，修正過(guò)后的D-P 準(zhǔn)則就變?yōu)榱嗽璂-P準(zhǔn)則。修正過(guò)后的D-P準(zhǔn)則處處都是光滑的，是一個(gè)正則屈服準(zhǔn)則且便于數(shù)值的計(jì)算。

泥頁(yè)巖的強(qiáng)度參數(shù)隨著地層含水量而發(fā)生變化，結(jié)合相關(guān)試驗(yàn)成果[21]可知，強(qiáng)度參數(shù)的變化趨勢(shì)近似是隨著含水量動(dòng)態(tài)線(xiàn)性衰減，將式(1)和式(6)進(jìn)行耦合，可建立水化泥頁(yè)巖強(qiáng)度弱化動(dòng)態(tài)模型，即：

式中：cb為當(dāng)含水量為ωb時(shí)的黏聚力；?b為當(dāng)含水量為ωb時(shí)的內(nèi)摩擦角；Ks為黏聚力系數(shù)；Ls為內(nèi)摩擦角系數(shù)。

3.4 模型的數(shù)值實(shí)現(xiàn)

基于上述分析，本文通過(guò)采用解耦的數(shù)值處理方法，首先針對(duì)相對(duì)獨(dú)立的水分?jǐn)U散方程進(jìn)行求解，然后將壓力傳遞模型與固體變形進(jìn)行間接協(xié)同耦合進(jìn)行求解，整個(gè)過(guò)程涉及到ABAQUS 軟件內(nèi)嵌的2個(gè)計(jì)算模塊，也就是巖土固結(jié)模塊和質(zhì)量擴(kuò)散模塊。

雖然巖土介質(zhì)的流-固耦合場(chǎng)與質(zhì)量擴(kuò)散場(chǎng)在數(shù)值計(jì)算上表現(xiàn)為較大的差異性，但究其本質(zhì)兩者都是包含了線(xiàn)性化與時(shí)步離散(或載荷增量)2 個(gè)基本內(nèi)容，也就是將流-固體場(chǎng)和擴(kuò)散場(chǎng)按照2 個(gè)獨(dú)立的系統(tǒng)分別進(jìn)行設(shè)計(jì)、計(jì)算，每一時(shí)步上的參數(shù)耦合是通過(guò)數(shù)據(jù)通訊的方式進(jìn)行實(shí)現(xiàn)，不斷地在每一時(shí)步上修正相關(guān)的系數(shù)，并且這種相互修正是在一系列的時(shí)步上交叉進(jìn)行的。在賈善坡等[24-25]的研究基礎(chǔ)上，以MATLAB 為平臺(tái)，ABAQUS 為求解器，開(kāi)發(fā)了計(jì)算泥頁(yè)巖井眼動(dòng)態(tài)破壞的軟件，通過(guò)編制ABQMAIN 子程序?qū)崿F(xiàn)各個(gè)計(jì)算模塊之間數(shù)據(jù)的存儲(chǔ)和通訊，強(qiáng)度參數(shù)的弱化則是通過(guò)USDFLD 子程序來(lái)進(jìn)行實(shí)現(xiàn)，程序的流程見(jiàn)圖6。

圖6 多場(chǎng)耦合程序系統(tǒng)程序結(jié)構(gòu)Fig.6 Flowchart of multifield coupling code

4 實(shí)例分析

本文以欠平衡鉆井方式下鉆開(kāi)地層為例，分析是否考慮水化作用對(duì)井眼漸進(jìn)破壞的影響規(guī)律。

4.1 工程概況

W號(hào)構(gòu)造位于渤海灣盆地黃驊坳陷NP凹陷的西部，是重要的勘探區(qū)帶，已鉆遇地層自下而上依次為古近系沙河街組、東營(yíng)組，新近系館陶組、明化鎮(zhèn)組和第四系平原組。W 號(hào)沙河街組巖性圈閉，面積為356 km2，最大厚度為840 m，閉合幅度為1 000 m，上覆沙一段大套泥頁(yè)巖可作為良好蓋層，有利于形成自生自?xún)?chǔ)的大型巖性油氣藏。

依據(jù)鉆井井史、鉆井液總結(jié)報(bào)告及完井卡片，對(duì)該區(qū)內(nèi)已鉆的20 余口井目標(biāo)層井下復(fù)雜情況進(jìn)行統(tǒng)計(jì)，分析表明：井下復(fù)雜情況以井漏現(xiàn)象最為突出，其次為井眼掉塊、坍塌，井漏發(fā)生次數(shù)約為井壁失穩(wěn)的4倍，目標(biāo)層位地層破裂壓力系數(shù)大都在1.80 以上(表5)，井下復(fù)雜井段所用鉆井液當(dāng)量密度最大為1.6，結(jié)合巖芯觀察和成像測(cè)井資料，井漏段地層井周裂縫較為發(fā)育，提高鉆井液當(dāng)量密度，對(duì)井壁垮塌、擴(kuò)徑的抑制效果不理想，部分井的鉆井液當(dāng)量密度由1.23 提高至1.39 時(shí)，井徑反而增大。

沙河街組硬脆性泥頁(yè)巖發(fā)育微裂縫，同時(shí)存在部分微孔洞，這些孔縫為鉆井濾液侵入地層內(nèi)部提供通道，加之地層黏土含量較高，常規(guī)鉆井過(guò)程中濾液的侵入或自吸水的發(fā)生容易引起水化膨脹，再加上高滲流能力裂縫的存在，使常規(guī)鉆井過(guò)程中容易出現(xiàn)井漏等鉆井問(wèn)題。為避免或減少在鉆井過(guò)程中出現(xiàn)相關(guān)的問(wèn)題，實(shí)施欠平衡鉆井方式是解決上述有關(guān)問(wèn)題的有效途徑。從表5可知：沙河街組地層壓力系數(shù)高于1.3，故欠平衡鉆井循環(huán)鉆井液選用水基鉆井液體系。本文以W-82井為研究對(duì)象，該井上覆無(wú)海水，井深為4 745 m，鉆遇地層主要為明化鎮(zhèn)組、館陶組、東營(yíng)組和沙河街組地層，用于勘探W 號(hào)構(gòu)造沙河街組巖性圈閉含油氣情況，四開(kāi)采用液相欠平衡鉆井技術(shù)，欠平衡鉆進(jìn)于4 146.93～4 745.00 m，欠平衡進(jìn)尺為598.07 m。

表5 W號(hào)構(gòu)造地層壓力測(cè)試結(jié)果Table 5 Pressure test results of W structural formation

4.2 計(jì)算模型及參數(shù)

為模擬沙河街組4 250 m附近井眼穩(wěn)定性，對(duì)模型進(jìn)行離散采用的是平面應(yīng)變四邊形單元，由于井眼的對(duì)稱(chēng)性，有限元計(jì)算模型取為井眼的1/4，模型長(zhǎng)為15 m，長(zhǎng)寬比為1:1，模型的中央是井眼，井眼半徑為0.108 m。利用漸變網(wǎng)格技術(shù)實(shí)現(xiàn)對(duì)井眼及模型中間線(xiàn)的局部加密，劃分的網(wǎng)格如圖7所示，網(wǎng)格總數(shù)為1 993，節(jié)點(diǎn)總數(shù)為2 087。

圖7 網(wǎng)格劃分Fig.7 Meshes of analysis model

定義模型邊界條件為：在BC邊施加最大水平地應(yīng)力，CD邊施加最小水平地應(yīng)力；OB邊為水平滑移邊界，OD邊為垂向滑移邊界；在模型的內(nèi)部定義初始孔隙壓力、初始應(yīng)力、初始地層含水量、初始地層水活度；地層鉆開(kāi)之后，AE邊施加井眼液柱壓力、含水量條件、鉆井液活度條件，并根據(jù)鉆井液活動(dòng)設(shè)置滲流邊界條件(定義滲透壓力)。

通過(guò)有限元軟件模擬井眼漸進(jìn)破壞過(guò)程的計(jì)算步驟如下：第一階段為地應(yīng)力平衡階段，該階段是為了還原初始應(yīng)力場(chǎng)也就是地層鉆開(kāi)之前的應(yīng)力狀態(tài)，通過(guò)定義模型內(nèi)部與外部的應(yīng)力場(chǎng)，使得還原的初始地應(yīng)力場(chǎng)符合計(jì)算的要求；第二階段為鉆井階段，本階段是為了模擬井眼的鉆開(kāi)，通過(guò)單元生死，并施加泥漿壓力模擬鉆井過(guò)程(地層被鉆開(kāi))，由于該步的時(shí)間比較短，并且考慮到泥頁(yè)巖的滲透性較低，該步不考慮井壁的滲透性；第三階段為滲流擴(kuò)散階段，本階段是為了模擬井筒內(nèi)、外發(fā)生的滲流擴(kuò)散效應(yīng)，井眼圍巖強(qiáng)度因水化而被減弱，模擬了井眼鉆開(kāi)完成后53 d 內(nèi)井壁的破壞過(guò)程。

根據(jù)室內(nèi)試驗(yàn)、現(xiàn)場(chǎng)資料、測(cè)井資料以及鄰區(qū)地質(zhì)資料獲得計(jì)算所需基本參數(shù)：上覆巖層壓力為91.63 MPa，水平最大地應(yīng)力為81.22 MPa，水平最小地應(yīng)力為68.72 MPa，地層壓力系數(shù)為1.3，地層彈性模量為20.2 GPa，泊松比為0.16，黏聚力為24.1 MPa，內(nèi)摩擦角為21.7°；根據(jù)測(cè)井解釋?zhuān)貙訚B透率為(0.13～12.56)×10-3μm2，孔隙度為8.00%。實(shí)際鉆井液當(dāng)量密度為1.2，活度擴(kuò)散系數(shù)為5×10-9m2/s，地層吸水?dāng)U散系數(shù)為9.5×10-9m2/s；初始地層含水量為2%，飽和含水量為10%。

4.3 不考慮水化情況下的計(jì)算結(jié)果

為了區(qū)分泥頁(yè)巖滲流和水化擴(kuò)散對(duì)井眼破壞的影響機(jī)制，暫不考慮兩者的耦合，這里僅探討滲流對(duì)井眼穩(wěn)定性的影響。

采用欠平衡鉆井方式打開(kāi)地層，以鉆井液當(dāng)量密度為1.0，1.1，1.2 和1.3 為例，對(duì)上述計(jì)算模型進(jìn)行數(shù)值模擬，探討井眼坍塌破壞過(guò)程，地層鉆開(kāi)后井眼破壞區(qū)分布如圖8所示。

從圖8可知：井眼在平行最小地應(yīng)力方向發(fā)生剪切破壞，在這一方向尺寸變大，破壞區(qū)為近似的橢圓形狀，與井下成像觀測(cè)和室內(nèi)試驗(yàn)結(jié)果一致。鉆井液當(dāng)量密度越大，塑性區(qū)越小，適當(dāng)提高鉆井液當(dāng)量密度，有利于井壁保持穩(wěn)定。地層鉆開(kāi)數(shù)天后，井眼破壞程度并沒(méi)有發(fā)生任何變化。

不考慮泥頁(yè)巖水化影響，井眼擴(kuò)大率隨時(shí)間的變化曲線(xiàn)如圖9所示。從圖9可知：井眼擴(kuò)大率隨著鉆井液當(dāng)量密度的增大而減小，當(dāng)鉆井液當(dāng)量密度為1.2 時(shí)，鉆開(kāi)瞬間井眼擴(kuò)大率在20%以?xún)?nèi)，隨著時(shí)間的增加井眼并沒(méi)有繼續(xù)擴(kuò)大。由此可以說(shuō)明，不考慮水化效應(yīng)時(shí)，井眼鉆開(kāi)后便趨于穩(wěn)定，鉆開(kāi)初期，滲透壓力的梯度較大，滲流場(chǎng)的變化會(huì)引起應(yīng)力場(chǎng)的變化，盡管滲流場(chǎng)的變化會(huì)引起應(yīng)力場(chǎng)的變化，但其變化程度較小，并不足以引起井眼的繼續(xù)破壞。

圖10所示為鉆井液當(dāng)量密度為1.2時(shí)，地層鉆開(kāi)后井壁孔隙壓力分布圖。從圖10可知：孔隙壓力隨著距井壁距離的增大而逐漸增大并趨于初始孔壓，其變化速率隨鉆開(kāi)時(shí)間的增大而減小，鉆開(kāi)井眼1 d后，滲流場(chǎng)的擾動(dòng)范圍約為10倍的井眼半徑，鉆開(kāi)井眼53 d 后的擾動(dòng)范圍約為50 倍的井眼半徑。

圖8 鉆開(kāi)后井眼破壞云圖分布Fig.8 Distribution of cloud damage in wellbore after excavation

圖9 不考慮水化作用下的井眼擴(kuò)大率隨時(shí)間變化曲線(xiàn)Fig.9 Curve of wellbore enlargement rate with time under without considering hydration

根據(jù)上述分析可知，對(duì)硬脆性泥頁(yè)巖地層實(shí)施欠平衡鉆井來(lái)說(shuō)，鉆井初期出現(xiàn)井壁坍塌等鉆井問(wèn)題的主要原因是鉆井液所能提供的液柱壓力不能達(dá)到地層支撐的壓力。為此在欠平衡鉆井過(guò)程中，應(yīng)減少循環(huán)和壓力波動(dòng)，保證欠平衡鉆井的順利進(jìn)行。

4.4 考慮水化情況下的計(jì)算結(jié)果

圖11所示為鉆井液當(dāng)量密度為1.2時(shí)，地層鉆開(kāi)后井壁孔隙壓力分布圖。從圖11可知：孔隙壓力隨著距井壁距離的增大而逐漸增大并趨于初始孔壓，其變化速率隨鉆開(kāi)時(shí)間的增大而減小，這與不考慮水化影響所得出的規(guī)律一致，不同的是考慮水化的情況下，擾動(dòng)區(qū)域的范圍更大，增加的這一部分便是由于水化擴(kuò)散的影響。

圖12所示為鉆井液當(dāng)量密度為1.2時(shí)，地層鉆開(kāi)后井眼圍巖含水量分布圖。從圖12可知：含水量隨著距井壁距離的增大而逐漸減小并趨于最小值2%，地層鉆開(kāi)后早期，含水量變化較為劇烈，存在一定的波動(dòng)，此后含水量變化速率隨浸泡時(shí)間的增大而減小。

圖10 不考慮水化作用下地層鉆開(kāi)后井壁孔隙壓力分布Fig.10 Pore pressure distribution of wellbore after drilling without considering hydration

圖11 考慮水化作用下地層鉆開(kāi)后井壁孔隙壓力分布Fig.11 Pore pressure distribution of wellbore after drilling considering hydration

圖12 地層鉆開(kāi)后井壁含水量分布Fig.12 Water content distribution of wellbore after drilling

圖13 不同距離下地層黏聚力的變化Fig.13 Variation of formation cohesion at different distances

圖14 不同距離下地層內(nèi)摩擦角的變化Fig.14 Variation of internal friction angle of formation at different distances

圖13和圖14所示分別為鉆井液當(dāng)量密度為1.2時(shí)，地層巖石黏聚力和內(nèi)摩擦角分布圖。由于二者都受到含水量的控制以至于兩者的分布規(guī)律與井壁圍巖含水量的分布規(guī)律相類(lèi)似，在地層鉆開(kāi)早期變化較為劇烈，此后變化減弱，然后逐漸地趨于一個(gè)定值。

鉆井液當(dāng)量密度為1.0 和1.3 時(shí)井眼漸進(jìn)破壞過(guò)程分別如圖15和16所示?？芍捍瓜?y方向)的破壞程度和范圍明顯比水平向(x方向)的大，破壞區(qū)基本呈橢圓形狀，當(dāng)鉆井液當(dāng)量密度為1.0 時(shí)，浸泡1，7，14 和24 d 后的最大塑性應(yīng)變分別為0.80%，1.03%，1.17%和1.30%，而當(dāng)鉆井液當(dāng)量密度為1.3時(shí)，浸泡1，7，14和24 d后的最大塑性應(yīng)變分別為0.50%，0.66%，0.73%和0.79%，可見(jiàn)，較小的鉆井液當(dāng)量密度加劇了井眼的水化破壞。就井眼破壞深度而言，由于水化效應(yīng)的影響，在井眼浸泡早期，井眼破壞區(qū)域增加較為迅速，此后趨于穩(wěn)定并逐漸呈線(xiàn)性增加趨勢(shì)。

圖15 鉆井液當(dāng)量密度為1.0時(shí)井眼漸進(jìn)破壞過(guò)程Fig.15 Progressive failure process of wellbore with a drilling fluid equivalent density of 1.0

圖16 鉆井液當(dāng)量密度為1.3時(shí)井眼漸進(jìn)破壞過(guò)程Fig.16 Progressive failure process of wellbore with a drilling fluid equivalent density of 1.3

圖17所示為不同鉆井液當(dāng)量密度下的井眼擴(kuò)大率隨時(shí)間變化曲線(xiàn)。從圖17可以看出：隨著浸泡時(shí)間增加井眼擴(kuò)大率逐漸增大，不同鉆井液當(dāng)量密度下的井眼擴(kuò)大率隨著鉆井液當(dāng)量密度增大而減小。對(duì)比圖9可知：以鉆井液當(dāng)量密度為1.2為例，在不考慮水化效應(yīng)時(shí)，雖然有滲流場(chǎng)的影響，但其影響程度還不足以使得井眼擴(kuò)大率隨著時(shí)間發(fā)生改變，而考慮水化效應(yīng)后，鉆井后對(duì)應(yīng)的擴(kuò)眼率為18.07%，7 d 后對(duì)應(yīng)的擴(kuò)眼率為58.22%，井眼擴(kuò)大率增加約40%，可見(jiàn)水化效應(yīng)對(duì)井眼漸進(jìn)破壞的影響要遠(yuǎn)大于單純滲流場(chǎng)的影響。

圖17 考慮水化作用下的井眼擴(kuò)大率隨時(shí)間變化曲線(xiàn)Fig.17 Curve of wellbore enlargement rate with time under considering hydration

根據(jù)上述分析可知：在欠平衡鉆井施工條件下，鉆井?dāng)_動(dòng)致使原有裂隙延伸變寬、新的微裂隙產(chǎn)生，微裂隙發(fā)生自吸水效應(yīng)致使泥頁(yè)巖水化膨脹，裂縫間聯(lián)結(jié)強(qiáng)度降低，圍巖整體強(qiáng)度逐漸降低，發(fā)生垮塌掉塊現(xiàn)象，是導(dǎo)致井壁漸進(jìn)失穩(wěn)的內(nèi)在因素。為此，在欠平衡鉆井過(guò)程中，鉆井速度要快，嚴(yán)格控制鉆井液對(duì)泥頁(yè)巖的抑制性能，并盡量減少循環(huán)時(shí)間，才能有效保證欠平衡鉆井的順利實(shí)施。

4.5 與現(xiàn)場(chǎng)結(jié)果對(duì)比

為了驗(yàn)證力-化耦合作用下井眼漸進(jìn)破壞模型的可靠性與準(zhǔn)確性，將數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)資料進(jìn)行了對(duì)比分析。

在實(shí)際的井眼鉆開(kāi)過(guò)程中，井眼擴(kuò)大率只要控制在一定的范圍之內(nèi)，并且及時(shí)地沖洗干凈井底巖屑，不造成井底巖屑的堆積，一般就不會(huì)引起卡鉆的事故。工程上一般以平均井徑擴(kuò)大率不超過(guò)15%，產(chǎn)層段最大井徑擴(kuò)大率不超過(guò)30%且最大井徑段占整個(gè)產(chǎn)層段比例不超過(guò)30%為合格。該井在四開(kāi)井段采用密度不大于1.2 g/cm3的水基鉆井液體系鉆進(jìn)沙河街組地層，鉆進(jìn)過(guò)程中井壁基本穩(wěn)定，僅有幾次輕微井塌，井下未發(fā)生復(fù)雜情況，井壁穩(wěn)定條件滿(mǎn)足欠平衡鉆井施工要求，平均機(jī)械鉆速為1.78 m/h，根據(jù)圖17所示的計(jì)算結(jié)果，當(dāng)井眼鉆開(kāi)后，井眼擴(kuò)大率計(jì)算值約為18%，與現(xiàn)場(chǎng)結(jié)果相符。

在采用水基鉆井液體系實(shí)施欠平衡鉆井過(guò)程中，井筒和地層之間產(chǎn)生負(fù)壓差，導(dǎo)致地層中的水向井筒中流動(dòng)，從而減少了井漏的風(fēng)險(xiǎn)，并且降低了鉆井液對(duì)地層的損害，但由于毛細(xì)管效應(yīng)的存在會(huì)造成地層巖石發(fā)生自吸水化作用，特別是發(fā)育的微裂隙，造成泥頁(yè)巖的吸水膨脹，水化效應(yīng)更加顯著，加上井眼缺乏了鉆井液的壓力平衡，加劇了井壁垮塌。中途起鉆和測(cè)井工程中曾返出大量掉塊，接立柱困難，短起下鉆必須使用倒劃眼才能進(jìn)行，從井徑曲線(xiàn)可以看到裸眼段存在多處大井眼井段(圖18)，其中4 265～4 275 m 段井眼擴(kuò)大率在55%～85%之間，井眼擴(kuò)大率計(jì)算結(jié)果約為84%。隨后進(jìn)行中途測(cè)試，至中途測(cè)試下鉆通井泥漿浸泡地層長(zhǎng)達(dá)53 d，掉塊多(圖19)，造成地層失穩(wěn)，井壁大段坍塌，而本文預(yù)測(cè)的井眼擴(kuò)大率已超過(guò)100%，垮塌嚴(yán)重，與實(shí)際鉆井基本吻合，說(shuō)明本文預(yù)測(cè)模型可以有效反映泥頁(yè)巖井眼漸進(jìn)動(dòng)態(tài)破壞過(guò)程。

圖18 W-82井徑擴(kuò)大率Fig.18 Diameter expansion rate of W-82 well

圖19 W-82井返出掉塊情況Fig.19 W-82 well returns block

4.6 討論

W 號(hào)構(gòu)造沙河街組泥頁(yè)巖破壞的內(nèi)在原因是微孔洞、微裂隙的發(fā)育程度較高，主要包括3個(gè)方面：1)原生裂縫較為發(fā)育；2)鉆井卸載作用使圍巖產(chǎn)生新的裂縫；3)硬脆性泥頁(yè)巖毛細(xì)管效應(yīng)較為顯著，而毛細(xì)管效應(yīng)產(chǎn)生的根本原因則是微裂縫的是否發(fā)育以及開(kāi)度大小，次生裂縫伴隨著毛細(xì)管自吸水化作用的發(fā)生而產(chǎn)生。過(guò)平衡鉆井方式下鉆井液密度較高，使得井筒內(nèi)的壓力高于地層壓力，產(chǎn)生井筒與地層之間的正壓差，促使原有的孔縫擴(kuò)張和新的孔縫的產(chǎn)生，使得鉆井流體侵入到地層內(nèi)部產(chǎn)生水化作用；若采用欠平衡鉆井方式，孔縫的存在會(huì)導(dǎo)致產(chǎn)生毛細(xì)管效應(yīng)造成自吸水化作用促使泥頁(yè)巖宏觀破壞。無(wú)論是以哪種方式侵入地層，產(chǎn)生的水化作用將會(huì)促使原有的裂縫擴(kuò)張和新的裂縫產(chǎn)生，進(jìn)一步加劇了水化作用?？梢?jiàn)，泥頁(yè)巖微裂隙發(fā)育和逐步擴(kuò)展是井壁漸進(jìn)失穩(wěn)的重要原因。

在欠平衡鉆井方式下(圖20(a))，地層的流體有控制地進(jìn)入井筒，由于泥頁(yè)巖的親水特性致使鉆井流體侵入地層產(chǎn)生倒吸現(xiàn)象，造成水化作用?；瘜W(xué)勢(shì)差作用下的水流方向取決于鉆井液活度的影響，當(dāng)鉆井液活度較大時(shí)，化學(xué)勢(shì)差作用下的水流方向指向地層，當(dāng)鉆井液活度較小時(shí)，水流方向則指向井筒?；瘜W(xué)勢(shì)差和負(fù)壓差的綜合效應(yīng)是否能抵制毛細(xì)管力下的水流方向，主要取決于鉆井液的活度和密度，但過(guò)低的鉆井液密度會(huì)導(dǎo)致地層的支撐壓力減小，從而增大了井眼圍巖擾動(dòng)區(qū)，加大了裂縫的擴(kuò)展，也會(huì)加劇自吸現(xiàn)象。針對(duì)這種問(wèn)題，有效的解決方法是提高鉆井液的抑制能力，控制鉆井液的活度，降低自吸水效應(yīng)。

圖20 鉆井過(guò)程中水相流動(dòng)示意圖Fig.20 Schematic diagram of water phase flow during drilling

過(guò)平衡鉆井方式下(圖20(b))，井筒壓力大于地層壓力，正壓差驅(qū)動(dòng)井筒中的流體向地層侵入，使得鉆井流體更多地進(jìn)入地層，加劇了泥頁(yè)巖的水化作用，弱化了圍巖強(qiáng)度。同欠平衡鉆井方式一樣，化學(xué)勢(shì)差作用下的水流方向取決于鉆井液活度的影響，但其影響的程度不足以抵抗正壓差和毛細(xì)管力的共同作用，總的水流方向指向地層，從而造成泥頁(yè)巖的水化作用。總體而言，過(guò)平衡鉆井方式下的井眼擾動(dòng)程度較欠平衡鉆井方式弱，但過(guò)平衡鉆井方式下的水化作用比欠平衡鉆井方式的強(qiáng)。針對(duì)這種問(wèn)題，有效的解決方法是在提高鉆井液抑制性的同時(shí)優(yōu)選出有效的封堵材料加入鉆井液中以強(qiáng)化鉆井液的封堵性。

水化作用是影響井壁坍塌的重要因素，無(wú)論是欠平衡鉆井方式還是過(guò)平衡鉆井方式，都會(huì)產(chǎn)生水化作用，過(guò)平衡鉆井方式較欠平衡鉆井方式下水化更劇烈，更加不利于井壁的穩(wěn)定，而影響水化作用的根本原因是微裂縫的發(fā)育與擴(kuò)展。綜上，有效的解決對(duì)策為：1) 優(yōu)選有效的抑制劑，加入鉆井液中以抑制硬脆性泥頁(yè)巖中的黏土礦物；2)優(yōu)選有效的封堵材料，加入鉆井液中以封堵微裂縫、微孔洞，降低鉆井液侵入地層的能力；3)優(yōu)選潤(rùn)濕反轉(zhuǎn)劑，加入鉆井液中以改變巖石表面潤(rùn)濕性，防止毛細(xì)管自吸效應(yīng)。

5 結(jié)論

1)綜合有效應(yīng)力原理并結(jié)合吸水?dāng)U散和壓力傳遞模型，建立了綜合考慮化學(xué)勢(shì)變化、流體流動(dòng)與骨架變形耦合的泥頁(yè)巖井眼漸進(jìn)破壞分析模型，為了利用該模型模擬硬脆性泥頁(yè)巖井眼動(dòng)態(tài)破壞的整個(gè)過(guò)程，以ABAQUS 軟件為求解器，編制了計(jì)算程序，計(jì)算得到的結(jié)果與實(shí)際鉆井基本吻合，驗(yàn)證了提出模型的有效性。

2)考慮水化影響條件下，井眼擴(kuò)大率隨著鉆井液浸泡時(shí)間的增大而增大，水化效應(yīng)對(duì)井眼漸進(jìn)破壞的影響要大于滲流場(chǎng)的影響。以鉆井液當(dāng)量密度1.2為例，不考慮水化效應(yīng)時(shí)，隨著時(shí)間的增加井眼并沒(méi)有繼續(xù)擴(kuò)大，而考慮水化效應(yīng)后，7 d 后的井眼擴(kuò)大率增加約40%，水化起著非常重要的作用。

3)中深層硬脆性泥頁(yè)巖微裂隙較為發(fā)育，地層中伊/蒙混層礦物和伊利石含量較高，泥頁(yè)巖易自吸水水化膨脹，為鉆井過(guò)程中易發(fā)生井壁失穩(wěn)的根本原因，在欠平衡鉆井初期井眼并未出現(xiàn)失穩(wěn)問(wèn)題，隨著浸泡時(shí)間的增加，井壁逐漸失穩(wěn)。因此，在硬脆性泥頁(yè)巖地層中實(shí)施鉆井技術(shù)，僅僅依賴(lài)于提高鉆井液密度來(lái)維護(hù)井壁穩(wěn)定的方法不太可取，有效的解決方法則是強(qiáng)化鉆井液的封堵性以及提高鉆井液的抑制能力。