低滲透油田壓裂液返排規(guī)律研究報告論文樣版

1、 - 總結(jié)中國石油大學(xué)（）現(xiàn)代遠程教育中國石油大學(xué)（）現(xiàn)代遠程教育畢畢 業(yè)業(yè) 設(shè)設(shè) 計（論文）計（論文）低滲透油田壓裂液返排規(guī)律研究低滲透油田壓裂液返排規(guī)律研究姓 名：學(xué) 號：性 別：專 業(yè): 石油工程批 次：學(xué)習(xí)中心：指導(dǎo)教師： 年 月. v摘摘 要要水力壓裂是低滲透油氣藏開發(fā)評價和增產(chǎn)增注必不可少的技術(shù)措施，而油氣井壓后的壓裂液返排又是水力壓裂作業(yè)的重要環(huán)節(jié)。目前，對壓裂液返排的控制，大多采用經(jīng)驗方法，沒有可靠的理論依據(jù)。本文對壓裂液的返排過程和壓后井底壓力的確定進行了較為系統(tǒng)的研究，旨在為壓裂液返排控制提供理論依據(jù)。本文在以壓裂液的濾失量計算的基礎(chǔ)上，運用流體力學(xué)和數(shù)值模擬的相關(guān)理論以及

2、物質(zhì)平衡原理，針對返排期間裂縫閉合的過程，考慮了啟動壓力梯度的影響，建立了壓裂液返排的數(shù)學(xué)模型，給出了壓裂液返排數(shù)學(xué)模型的數(shù)值解法。研究表明，為了減少壓裂液對儲層的傷害，低滲透儲層中的壓裂井應(yīng)采用停泵后立即返排的方式，使裂縫強制閉合。實測井口壓降曲線與計算值的比較結(jié)果表明，建立的模型能夠比較準(zhǔn)確地預(yù)測裂縫閉合過程和壓裂液返排過程。最后，對返排的過程進行了室內(nèi)模擬實驗研究，通過巖心實驗，發(fā)現(xiàn)了返排過程中的一些特定規(guī)律。然后以濾失機理研究為基礎(chǔ)，通過了建立裂縫壁面上的滲流模型，編制了返排參數(shù)預(yù)測程序，可通過對壓裂返排效果的預(yù)測來指導(dǎo)壓裂液返排作業(yè)。關(guān)鍵詞：關(guān)鍵詞：水力壓裂；裂縫閉合；壓裂液返排；數(shù)

3、學(xué)模型；井底壓力. v目目 錄錄第 1 章緒論 11.1 壓裂液返排規(guī)律研究的目的和意義 11.2 目前關(guān)于壓裂液返排規(guī)律研究存在的不足 2第 2 章低滲透油田特點及壓裂液返排規(guī)律研究現(xiàn)狀 32.1 國內(nèi)外低滲透油田儲量分布及特點 32.1.1 國外低滲透油田儲量分布 32.1.2 國外低滲透油田的主要特點 32.1.3 國內(nèi)低滲透油田儲量分布 42.1.4 國內(nèi)低滲透油田的主要特點：42.2 壓裂液返排規(guī)律研究現(xiàn)狀 52.2.1 國外壓裂液返排的推薦做法 52.2.2 國內(nèi)壓裂液返排的研究現(xiàn)狀 72.3 裂縫形態(tài)的數(shù)學(xué)模型 8第 3 章裂縫閉合期間壓裂液返排模型 93.1 裂縫閉合過程中模型

4、的假設(shè)條件 93.2 壓裂液返排的二維數(shù)學(xué)模型 93.2.1 壓裂液從地層返排的數(shù)學(xué)模型 93.2.2 初始條件及邊界條件 133.3 模型的數(shù)值解法 143.3.1 返排模型的離散 143.3.2 求解方法 17. v3.4 裂縫閉合時間及壓裂液返排量的確定 183.4.1 裂縫閉合時間確定 183.4.2 壓裂液返排量的計算 183.4.3 停泵后裂縫體積變化量的計算 193.5 實例計算與分析 203.6 室內(nèi)實驗?zāi)M研究 223.6.1 實驗方法 223.6.2 實驗數(shù)據(jù)及處理 23第 4 章壓裂液返排的實驗研究 264.1 實驗儀器材料 264.2 實驗步驟 264.3 實驗數(shù)據(jù)處理

5、與結(jié)果分析 264.3.1 采用瓜膠壓裂液進行壓裂實驗 264.3.2 采用田菁膠壓裂液進行實驗的結(jié)果 29圖 4-6 累計流量與滲透率恢復(fù)值 304.4 結(jié)論與建議：30第 5 章壓裂過程中的濾失與返排效果預(yù)測 315.1 壓裂液濾失理論 315.1.1 受壓裂液黏度控制的濾失系數(shù) C1325.1.2 受地層流體壓縮性控制的濾失系數(shù) C2325.1.3 壓裂液造壁性控制的濾失系數(shù) C3345.1.4 綜合濾失系數(shù) C34. v5.2 一維總濾失體積計算 355.3 壓裂液返排數(shù)學(xué)模型 365.3.1 飽和度呈線性分布 365.4 實例計算 375.5 影響壓裂液造壁性濾失系數(shù)實驗 38結(jié)論

6、39參考文獻 40. v第第 1 1 章章 緒論緒論1.1 壓裂液返排規(guī)律研究的目的和意義壓裂工藝是油、氣藏增產(chǎn)和提高采收率的最有效的措施之一1-2。隨著水力壓裂技術(shù)在低滲透油氣田勘探、開發(fā)及其它工業(yè)生產(chǎn)領(lǐng)域中的廣泛使用，其在理論方法、工藝、設(shè)備及工具方面都得到了迅速的發(fā)展。水力壓裂具有多學(xué)科性，它是與巖石力學(xué)（控制著裂縫幾何形態(tài)） 、流體力學(xué)（控制著液體流動與支撐劑的鋪置和回流）以及化學(xué)（控制著施工的材料性能，如壓裂液）密切地聯(lián)系著的。由于這種多學(xué)科交叉的復(fù)雜性，人們對水力壓裂工藝整個過程的研究還有許多不完善的地方，有些水力壓裂的指導(dǎo)方法或控制程序到目前為止還只是建立在現(xiàn)場工程師的經(jīng)驗之上。

7、所以，結(jié)合這些學(xué)科的知識對指導(dǎo)水力壓裂工藝的理論基礎(chǔ)進行深入的探討，具有很大的現(xiàn)實意義，即可引導(dǎo)工程技術(shù)人員設(shè)計出更合理的水力壓裂工藝控制程序。壓裂是重要的油氣井增產(chǎn)措施，而壓裂液的返排是壓裂施工中的一個主要的步驟。壓裂液能否順利、及時地排出對施工效果影響非常明顯，特別是對于低滲低壓地層，由于返排困難，殘液可能造成地層再次傷害，影響增產(chǎn)效果。研究表明：施工結(jié)束后，為了防止殘液長期滯留地層造成二次傷害，一般說來應(yīng)將殘液盡可能地快速排出。在實際施工過程中常常由于對地層客觀認(rèn)識不足，對工作液、添加劑選擇不當(dāng)，或施工工藝不合理，使返排率低，造成對儲層的傷害。可見，返排在壓裂、酸化工藝中占有非常重要的位

8、置，如果不能及時排完、排凈會給地層帶來再次傷害，返排效果的好壞直接影響措施效果。由于這一過程是在裂縫閉合期間進行的，因而，適當(dāng)?shù)姆蹬懦绦蛲ǔＪ潜３至芽p良好導(dǎo)流能力的關(guān)鍵所在，壓后油氣田的生產(chǎn)能力在很大程度上取決于裂縫導(dǎo)流能力。在壓裂液返排過程中，工程技術(shù)人員往往希望通過對返排流量的控制使支撐劑在產(chǎn)層區(qū)獲得較好的鋪置，繼而使裂縫具有較高的導(dǎo)流能力。但是，由于對壓后關(guān)井期間或返排過程中裂縫閉合情況、支撐劑運移情況、壓裂液濾失及流體性質(zhì)的變化情況不能很好地把握，所以對返排流量的控制經(jīng)常顯得無據(jù)可依。到目前為止，工程現(xiàn)場反饋了很多問題，其中的一些問題導(dǎo)致了嚴(yán)重的后果。對于油氣層壓力低的井，返排困難的問

9、題十分突出；而當(dāng)油氣層壓力較高時，油氣井產(chǎn)量大時，對于油井，過大的放噴速度會使支撐劑回流；對于氣井，過大的放噴速度會產(chǎn)生氣竄，使氣體將殘液分割在地層中不能排出，因此仍然存在排液速度控制的問題。這些問題主要有兩大方面3-4：一是沒有選擇合理的時機對支撐劑的運移進行控制，導(dǎo)致大量的支撐劑回流到井筒（吐砂） ，或在近井筒帶破碎，從而使產(chǎn)層區(qū)的支撐劑很少或. v分布不合理。這在很大程度上降低了裂縫的導(dǎo)流能力，嚴(yán)重的會導(dǎo)致壓裂施工的失??；二是返排流量控制不好，使回流的支撐劑沖出井口，刺壞放噴油嘴以及破壞其他設(shè)備。而且，往往發(fā)現(xiàn)這些問題時，施工基本完成，想進行挽救為時已晚。在低滲透儲層中，一般都要進行大型

10、水力壓裂作業(yè)，壓后返排出現(xiàn)的這些問題就更加突出了。存在這些問題的一個很重要的原因就是由于壓裂液返排控制中經(jīng)驗成分過多，缺少可靠的理論依據(jù)。從合理控制返排的目標(biāo)出發(fā)，必須對壓裂液返排過程的機理進行深入分析，了解裂縫的閉合過程，認(rèn)識壓裂液返排的規(guī)律。1.2 目前關(guān)于壓裂液返排規(guī)律研究存在的不足綜上所述，雖然在低滲油氣藏壓裂液返排規(guī)律研究上已取得了很大的進步，但還存在以下幾方面不足5-10：（1）以往在壓裂液返排工藝研究上，缺乏量化的操作流程，返排時的井口壓力完全依賴經(jīng)驗。壓裂液在儲層中的返排過程類似于油(氣)驅(qū)水過程，井口壓力大小直接關(guān)系到油(氣)驅(qū)水流動過程中的壓力梯度，會最終影響排驅(qū)效率。因此

11、有必要建立井口不同返排壓力與返排效果的關(guān)系，根據(jù)油藏實際情況選擇合適的返排井口壓力。（2）以往建立的壓裂液返排模型雖然在理論上對壓裂施工具有一定的指導(dǎo)意義，但這些模型一般適合于中、高滲透油氣藏的開發(fā)，沒有考慮到低滲透油氣藏的實際情況，忽略了低滲透油氣藏中啟動壓力梯度對返排效果的影響，因此有必要建立適合于低滲透油氣藏壓裂液返排的數(shù)學(xué)模型和物理模型，分析影響因素，指導(dǎo)壓裂作業(yè)。（3）影響返排效果的因素是多方面的，對于不同油氣藏，其儲層物性和流體性質(zhì)也不同，各影響因素(返排壓差、滲透率、流體的黏度、界面張力、潤濕性和孔隙度等)對返排效果影響程度也不盡相同，因此有必要了解不同油氣藏的主要因素，從而為現(xiàn)

12、場壓裂作業(yè)提供合適可行的方案。油氣藏壓裂后壓裂液返排過程中雖然是油、水相流動，但在壓裂過程中，其有效孔隙度和滲透率發(fā)生了變化，因此特別有必要建立殘液返排的滲流模型，這樣更有利于對返排率的預(yù)測和壓裂后生產(chǎn)井的效果預(yù)測。目前國內(nèi)現(xiàn)場放噴排液基本上還處于靠經(jīng)驗操作的階段，沒有一套科學(xué)的排液理論來加以指導(dǎo)和量化，導(dǎo)致排液措施隨意性大，往往對施工效果造成非常不利的影響，但是這種影響又經(jīng)常被忽視。因此，目前急需從機理上研究殘液的返排過程，在機理研究的基礎(chǔ)上提出具有針對性的、可量化操作的排液措施及排液參數(shù)，對于提高壓裂施工成功率、提高施工效果和油氣井產(chǎn)量是非常必要的。. v第第 2 2 章章 低滲透油田特點

13、及壓裂液返排規(guī)律研究現(xiàn)狀低滲透油田特點及壓裂液返排規(guī)律研究現(xiàn)狀2.1 國內(nèi)外低滲透油田儲量分布及特點.1 國外低滲透油田儲量分布國外低滲透油田儲量分布世界上低滲透油田資源十分豐富，分布范圍廣泛，各產(chǎn)油國基本上都有該類型的油田11-13。美國中部、南部、北部和東部，前蘇聯(lián)的前喀爾巴阡山、克拉斯諾達爾、烏拉爾伏爾加、西西伯利亞油區(qū)和加拿大西部的阿爾伯達省都有廣泛的分布。近些年來，小而復(fù)雜的低滲透油田比例越來越大。例如，俄羅斯近年來在西西伯利亞地區(qū)新發(fā)現(xiàn)的低滲透、薄差層儲量已占探明儲量的 50%以上。有的地區(qū)，低滲透油田連片分布，成為低滲油區(qū)。1999 年俄羅斯各油氣區(qū)中難以開采的石

14、油儲量占剩余可采儲量的 40%以上。低滲透儲層的儲量達數(shù)百億噸，其中，滲透率低于 5010-3的低滲透層儲量約有 150108t，占2m低滲透儲量的 90%以上，占俄羅斯可采儲量的 30%以上。其中 60%在西西伯利亞。這些低滲透層大部分已投入開發(fā)。目前，全俄羅斯從低滲透儲集層中采出的石油占全部采出量的20%左右。低滲透儲集層中石油黏度在 10mPas 以下，其有效厚度多數(shù)在 210m，埋藏深度大多為 12001400m，采出程度不高。1998 年美國低滲透油氣田可采儲量占全國總儲量的 10%15%。據(jù)北美 172 個低滲透油藏統(tǒng)計，滲透率一般在幾個毫達西到幾十個毫達西。其中，2010010-

15、3的油田占2m這些低滲透油藏總數(shù)的 60%，20110-3的占 30%，少數(shù)低于 110-3，約占2m2m5%。.2 國外低滲透油田的主要特點國外低滲透油田的主要特點從國外報道的情況看，對低滲透油田大體上可以歸納出以下幾個特點：（1）儲層物性差，滲透率低。由于顆粒細、分選差、膠結(jié)物含量高，經(jīng)壓實和后生成巖作用使儲層變得十分致密，滲透率一般小于 0.1 ，少數(shù)低于 0.001 (統(tǒng)計北美2m2m172 個低滲透砂巖油藏的數(shù)據(jù))。（2）儲層孔隙度一般偏低，變化幅度大。大部分為 7%20%，個別高達 28%。（3）原始含水飽和度較高，原油物性較好。一般含水飽和度 30%40%，個別高

16、達60%（美國東堪頓油田） ，原油比重多數(shù)小于 0.85，地層油黏度多數(shù)小于 3mPas。. v（4）油層砂泥交互，非均質(zhì)性嚴(yán)重。由于沉積壞境不穩(wěn)定，砂層的厚薄變化大，層間滲透率變化大，有的砂巖泥質(zhì)含量高，地層水電阻率低，給油水層的劃分帶來很大困難。（5）天然裂縫相對發(fā)育。由于巖性堅硬致密，多存在不同程度的天然裂縫系統(tǒng)，一般受區(qū)域性地應(yīng)力的控制，具有一定的方向性，對油田開發(fā)的效果影響較大，裂縫是油氣滲透的通道，也是注水竄流的條件，且人工裂縫又多與天然裂縫的方向一致，因此，天然裂縫是低滲透砂巖油田開發(fā)必須認(rèn)真對待的問題。（6）油層受巖性控制，水動力聯(lián)系差，邊底水驅(qū)動不明顯，自然能量補給差，多數(shù)靠

17、彈性和溶解氣驅(qū)采油，油層產(chǎn)能遞減快，一次采收率低，只能達到 8%12%，采用注水保持能量后，二次采收率可提高到 25%30%。（7）由于滲透率低，孔隙度低，必須通過酸化壓裂投產(chǎn)，才能獲得經(jīng)濟價值或必須通過壓裂增產(chǎn)，才能提高經(jīng)濟效益。（8）由于孔隙結(jié)構(gòu)復(fù)雜，喉道小，泥質(zhì)含量高，以及各種水敏性礦物的存在，導(dǎo)致開采過程中易受傷害，損失產(chǎn)量可達 30%50%，因此，在整個采油工藝系列中，保護油層是至關(guān)重要的環(huán)節(jié)。.3 國內(nèi)低滲透油田儲量分布國內(nèi)低滲透油田儲量分布低滲透油田儲量就占我國陸上已探明未動用總儲量的 60%以上。大慶外圍油田均屬這類油藏。就目前石油工業(yè)的發(fā)展?fàn)顩r來看，我國大多數(shù)

18、油田已經(jīng)進入高含水和特高含水期，原油的開采難度逐漸加大，勘探的形勢是新近探明儲量中低豐度、低滲透、低產(chǎn)能（俗稱“三低” ）儲量占據(jù)的比例很大。90 年代以來在大慶、 、遼河、勝利、長慶等主要油田陸續(xù)發(fā)現(xiàn)了許多低滲透油藏。據(jù)統(tǒng)計，在近幾年探明的未動用石油地質(zhì)儲量中，低滲透油層儲量占 58%以上，而在已經(jīng)探. v明的石油地質(zhì)儲量中，低滲透油藏的石油地質(zhì)儲量所占比例高達 6070%，甚至更高。可見低滲透油藏是我國今后相當(dāng)長一段時間內(nèi)增儲上產(chǎn)的主要資源基礎(chǔ)。我國陸地低滲透油田廣泛分布于全國 20 多個油區(qū)，它們分布在不同的巖性地層中，物性參數(shù)相差很大。而且，在我國陸上原油探明儲量中低滲透油田占有非常大

19、的比例，目前探明儲量約40108t。.4 國內(nèi)低滲透油田的主要特點：國內(nèi)低滲透油田的主要特點：（1）目前發(fā)現(xiàn)的低滲透油田儲層以中深埋藏深度為主由各油區(qū)的低滲透儲層埋藏深度統(tǒng)計表明，目前發(fā)現(xiàn)的油藏以中深層為主，埋藏深度小于 1000m 約占 5.2%，10002000m 約占 43.1%，20003000m 約占 36.2%，大于 3000m約占 15.5%。（2）低滲透儲層中特低滲透及超低滲透層儲量占較大的比例根據(jù)滲透率大小，低滲透油藏可分為 3 類：類是滲透率為 1050，其儲-3210 m量占 53%，類滲透率為 11010-3，其儲量占 38.6%，類滲透率為 0.111

20、0-32m，占儲量 8.4%，二、三類低滲透儲層的儲量占到 47%。2m（3）國內(nèi)低滲透油藏巖性以砂巖為主從目前探明的低滲透油藏統(tǒng)計，砂巖油藏占 70%左右，礫巖油藏占 10%左右，其余存在于變質(zhì)巖和灰?guī)r等特殊巖性油藏中。低滲透問題是一個十分復(fù)雜的課題，屬于非線性問題。目前，國內(nèi)外關(guān)于這方面的研究處于探索階段，有許多問題尚不清楚。就已有的研究成果來看，主要存在以下問題：低速滲流時，有沒有啟動壓力梯度還存在很多爭議；如何測量低滲透巖石的啟動壓力梯度和啟動壓力產(chǎn)生的原因值得進一步研究。對于低速非達西滲流，沒有判斷低速非達西準(zhǔn)則，往往僅以地層滲透率作為界限。對低滲透非達西滲流的滲流機理的認(rèn)識還處于探

21、索階段。. v低滲透介質(zhì)中的滲流規(guī)律甚為復(fù)雜，至今還沒有一個清楚的令人滿意的表達方法。因此，合理開發(fā)低滲透油田是非常必要的。首先，必須正確認(rèn)識其儲層特征和滲流規(guī)律，準(zhǔn)確的進行滲流計算，確定合理的開發(fā)方案。達西定律一直作為一個基本的規(guī)律被廣泛地應(yīng)用于油氣田開發(fā)的滲流計算中。然而，傳統(tǒng)的達西定律面對和高滲透油田有著諸多不同的低滲透油田開發(fā)計算問題，顯得不盡準(zhǔn)確，所以，研究人員必須打破傳統(tǒng)的達西公式，尋找更合理的計算方法，來解決工程問題。2.2 壓裂液返排規(guī)律研究現(xiàn)狀.1 國外壓裂液返排的推薦做法國外壓裂液返排的推薦做法近十多年來，國外學(xué)者在壓后壓裂液返排的問題上形成了多種認(rèn)識，具有

22、代表性的觀點主要有小排量返排、強化返排和反向脫砂三種14-24。 小排量返排1988 年，Robinson 等人討論了采用小油嘴排液以減小裂縫閉合應(yīng)力的優(yōu)點，提出了“小排量早期返排” （Early Flowback）的觀點。實際上， “小排量早期返排”是一種有控制的返排。Robinson 等人所做的研究突出了裂縫閉合應(yīng)力對支撐劑破碎以及裂縫閉合時間對支撐劑沉降的影響，認(rèn)為低滲透儲層壓裂后通常需要較長的閉合時間。在裂縫閉合之前壓裂液已完全破膠，支撐劑已大量沉降，排液初期通過控制返排速度的辦法，盡可能減小地層閉合應(yīng)力，讓支撐劑停留在裂縫內(nèi)，從而減少支撐劑的破碎和倒流。實現(xiàn)有效的支撐縫

23、長。基于以上認(rèn)識，Robinson 等人提出了控制返排的推薦做法：（1）水力壓裂前，應(yīng)有地層閉合應(yīng)力預(yù)測值或測定值，這是選擇支撐劑、確定排液程序的基礎(chǔ)。（2）獲取本地區(qū)裂縫閉合所需的時間。如果壓裂液破膠時間大于裂縫閉合時間，應(yīng)當(dāng)用 0.81.6mm 的小油嘴，以 1938（L/min）的小排量返排，使裂縫降壓。待裂縫閉合后. v立即關(guān)井等待壓裂液破膠，即使這樣仍會有部分支撐劑倒流入井筒。（3）如果施工井作業(yè)后有自流能力，應(yīng)當(dāng)使用 2.4mm3.2mm 的小油嘴返排，并控制回壓至最低（即降低近井筒帶的壓降） 。（4）關(guān)井、生產(chǎn)井決不能用大油嘴瞬時開井，推薦以每次 0.8mm 的放大量逐步放大油嘴

24、開關(guān)。（5）油氣井生產(chǎn)期間，應(yīng)定期測定或計算井底流壓。當(dāng)井底流壓持續(xù)回落，或者當(dāng)?shù)貙娱]合應(yīng)力接近所用支撐劑的最大允許應(yīng)力時，就不能再放大油嘴，除非萬不得已。Robinson 等人的“小排量早期返排”實際上是一種有控制的返排，他們所做的研究突出了裂縫閉合應(yīng)力對支撐劑破碎、裂縫閉合時間對支撐劑沉降的影響，認(rèn)為低滲透儲層壓裂后通常需要較長的閉合時間。在此之前壓裂液己完全破膠，支撐劑己大量沉降，排液初期通過控制返排速度的辦法，盡可能減小地層閉合應(yīng)力，讓支撐劑留在裂縫內(nèi)，從而減少支撐劑的破碎和倒流。 強化返排強化返排與 Robinson 等人相反，1990 年 E1y 等人提

25、出用強制裂縫閉合（Forced Closure）工藝，配以較高的支撐劑濃度（高砂比）和嚴(yán)格的壓裂液質(zhì)量控制措施，能極大地改善低滲透油氣井支撐縫的導(dǎo)流能力。E1y 等人推薦的排液做法，是在頂替壓裂液的 30s 內(nèi)就完成裂縫的閉合，當(dāng)從地面壓力檢測到近井筒帶裂縫己經(jīng)閉合后，以小于 3857（L/min）的速度返排 30min，然后放大返排量至 160320（L/min） ，只要不出砂。這種“強制裂縫閉合”實際上是一種強化返排方式。強化返排減少了壓裂液在地層里的停留時間，從而減少了液體傷害，有助于改善裂縫導(dǎo)流能力。后來有學(xué)者認(rèn)為，Ely 等. v人提出的這種返排程序，非常適合特低滲透地層，能極大地改

26、善低滲透油氣井的返排效果，但不具普遍性。而且這種強制裂縫閉合工藝通常會使支撐劑形成嚴(yán)重的“砂堤” （Proppant Banking） ，并且由于大量支撐劑的運移，液體濾失加重，在近井筒帶形成“裂縫尖端”（Pinch Point） 。另外，這種返排方式不適合易出水地層。 反向脫砂反向脫砂一般而言，低滲透地層水力壓裂需要大液量的壓裂液，大量的支撐劑是通過大液量的壓裂液以較小的填砂強度泵入地層的，這樣延長了裂縫閉合時間，此時就有必要采取“反向脫砂”方式，使支撐劑在井筒附近脫砂，形成橋堵，而通過尾追支撐劑的方法即能加速這一過程。基于這種認(rèn)識，1995 年 Barree 和

27、 Mukherjee 提出了“反向脫砂” （Reverse Screenout）工藝。Barre 和 Mukherjee 使用全二維裂縫幾何模擬器，系統(tǒng)研究了裂縫閉合期間支撐劑的沉降規(guī)律，討論了返排速度、射孔段位置、最終鋪砂濃度、裂縫幾何形態(tài)對保持閉合后裂縫內(nèi)鋪砂濃度的影響，解釋了“強制裂縫閉合”與“反向脫砂”在返排程序設(shè)計上的區(qū)別，并且深入到多個射孔段之間的交叉流動對支撐劑運移規(guī)律的影響，至今仍是水力壓裂作業(yè)返排設(shè)計的指導(dǎo)性原則。Barre 和 Mukherjee 認(rèn)為，井筒內(nèi)液體的膨脹性、液體濾失、支撐劑的對流、停泵期間水擊效應(yīng)引起的壓力波動、關(guān)井后裂縫的延伸和泵送結(jié)束時剩余壓力梯度的消解

28、、不同閉合應(yīng)力的層段間的交叉流動，是引起支撐劑運移與沉降的六個主要因素。他們針對不同的裂縫形態(tài)，分別就對稱裂縫、向上延伸裂縫、向下延伸裂縫、層內(nèi)延伸裂縫、高濾失下的對稱裂縫、端部脫砂壓裂及多層限流壓裂做了模擬研究，給出了以下推薦做法：（1）以高填砂強度、小液量脫砂設(shè)計，比依賴壓后返排保持裂縫導(dǎo)流能力更可行。. v（2）只有當(dāng)層內(nèi)裂縫保持良好延伸，或者有較高應(yīng)力差的上下層遮擋時，才有通過適當(dāng)?shù)姆蹬懦绦蚋纳谱罱K鋪砂濃度的可能，而必須的返排程序是施工結(jié)束后立即開始返排，并且返排速度要高于裂縫內(nèi)液體的濾失速度。如果是多層裂縫，返排速度還必須大大超過層間交叉流動的速度。（3）在向上延伸的裂縫或裂縫高度過

29、快發(fā)育的情況下，從加強膠體殘渣返排的角度講，有控制的返排或許有益。此時裂縫閉合期間聚合物膠體聚積于濾失帶，很小的返排速度避免了裂縫過快閉合，使得膠體殘渣排出的同時，支撐劑通過沉降或濾失的方式向下運移，有助于減小有效縫高（即支撐縫高） 。（4）裂縫向下延伸時，無論泵送期間還是裂縫閉合期間，支撐劑都是沉降的，此時再大的返排速度也不足以影響裂縫體積，支撐劑可能填充在射孔孔眼處，但不可能覆蓋整個裂縫段。（5）在向上延伸的裂縫里，裂縫閉合期間支撐劑的沉降將大大增加裂縫內(nèi)鋪砂濃度，強制裂縫閉合方式會使裂縫過早閉合。此時不宜用強化排液法。（6）在層內(nèi)裂縫延伸良好的情況下，大于地層濾失速度的返排速度有助于支撐

30、劑向射孔孔眼方向運移，這是有控制的返排方式的最佳應(yīng)用。（7）高濾失地層厚度比整個裂縫高度相對薄時，利用地層的天然濾失性能比返排更有效地使支撐劑向高濾失帶運移。此時也不宜采用強化返排法。（8）有較大應(yīng)力差的多層水力壓裂作業(yè)，壓后返排速度必須快于層間交叉流動的速度，以免低應(yīng)力層支撐劑過剩而高應(yīng)力層支撐劑不足，而且必須在停泵后立即開始返排。Barre 和 Mukherjee 的“反向脫砂”工藝，實際上是一種快速返排加井筒脫砂方式，這種方式至少有三點好處：首先，減輕了壓裂液對地層的傷害；其次，井筒脫砂顯著改善了支撐劑在近井筒帶的填充，提高了裂縫的無因次導(dǎo)流能力；最后，快速返排使得支撐劑未能大量沉降到裂

31、縫端部前，裂縫己經(jīng)閉合，從而形成較長的支撐縫。為了減少快速返排和. v井筒脫砂可能對循環(huán)系統(tǒng)造成的磨損，國外一般還要采取控制支撐劑倒流的措施，使用尾追樹脂包裹支撐劑（樹脂砂）就是其中的一種。在克羅拉多州的 Dodell 低滲透（0.0030.05）致密砂巖氣層的水力壓裂作業(yè)中，就采用了支撐劑中混入纖維材23m10料的辦法，返排速度高達 320L/min 時也無支撐劑返出。.2 國內(nèi)壓裂液返排的研究現(xiàn)狀國內(nèi)壓裂液返排的研究現(xiàn)狀目前國內(nèi)學(xué)者認(rèn)為，在殘液返排上主要集中在 3 個方面：工程上提高排液效率；集中在影響排液因素上的研究，以理論結(jié)合實際提高排液效率；排液預(yù)測模型的研究上。劉川

32、生8等提出對于氣井酸化，利用間隙放噴排液的方法來增加返排率，當(dāng)氣井酸化后經(jīng)自噴和誘噴轉(zhuǎn)入放噴階段。在這一階段對每次放噴時間間隔和放噴強度進行控制，由于酸化規(guī)模、井身結(jié)構(gòu)和氣產(chǎn)量的不同，作者采用定性的方法總結(jié)了一套氣井間隙放噴的方法，應(yīng)用效果較好。劉應(yīng)學(xué)25等提出酸化解堵增能返排技術(shù)，其技術(shù)是在擠酸液之前或之后將增能液注入地層，通過控制其反應(yīng)速度，在地層中發(fā)生化學(xué)反應(yīng)，產(chǎn)生大量熱量和氣體，使射孔段地層“瞬間”升溫。利用這種熱化學(xué)反應(yīng)，配合酸化處理，可以溶解近井地帶及井筒內(nèi)復(fù)雜的有機（無機）沉積垢，疏通油流通道。反應(yīng)放出的高溫高壓氣體，借助起泡劑的作用，在地層產(chǎn)生強大的推動力，造成強烈的壓力波動，

33、沖擊喉道堵塞物，不僅有助于解堵，還起到助排作用。王尤富26論述了入井液表面張力與儲層損害關(guān)系的實驗室研究。降低入井液的表面張力，可以減小毛細管阻力和提高入井液的返排率，從而達到保護儲集層的目的。實際油氣儲集層具有不同的滲透率和不同的潤濕性，毛細管阻力對不同滲透率和不同潤濕性的儲集層造成的損害不同，因此，對入井液的表面張力應(yīng)有不同要求，以得出它們的變化規(guī)律，為現(xiàn)場保護油氣層和確定入井液的合理表面張力提供可靠依據(jù)。賀承祖27論述了水鎖效應(yīng)與儲層傷害的關(guān)系。未開發(fā)的油氣層處于殘余水飽和狀態(tài)，可以認(rèn)為油、氣驅(qū)動壓力與毛細管力處于平衡狀態(tài)。當(dāng)儲層鉆開后在地層未受到傷害時，會出現(xiàn)水鎖效應(yīng)，一般認(rèn)為外. v

34、來流體在地層中的毛細管力越高，水鎖效應(yīng)越強，油氣產(chǎn)量越低。研究表明，并不是所有能降低表面張力的物質(zhì)都能降低水鎖效應(yīng)和儲層傷害，水鎖程度的大小是受毛管力控制的，表面張力只是影響毛管力的一個因素。2.3 裂縫形態(tài)的數(shù)學(xué)模型在進行返排程序設(shè)計的時候，研究人員必須知道停泵時刻裂縫的相關(guān)參數(shù)，進而就要用到分析裂縫形態(tài)的數(shù)學(xué)模型28-32。到目前為止，確定裂縫形態(tài)的數(shù)學(xué)模型已經(jīng)從二維模型、擬三維模型發(fā)展到全三維模型。對于全三維模型，國外，Clifton 與 Abou-saved 及 Cleary 等人提出兩種具有代表性的全三維裂縫延伸模型；國內(nèi)的趙金洲、吳迪祥等人也在裂縫三維延伸方面作了大量的研究工作，并

35、取得了一些成果。在全三維模型中，縫寬方程是奇異積分方程，對于這類方程，當(dāng)源點和場點重合時，被積函數(shù)無窮大，僅在柯西主值的意義上收斂。因此，這類方程的直接數(shù)值求解是繁瑣、困難的?；谶@種情況，目前對全三維模型的求解依舊是一個很大的研究課題。在成熟的軟件里或在現(xiàn)場施工的實際應(yīng)用中，還是以二維模型和擬三維模型為主。. v第第 3 3 章章 裂縫閉合期間壓裂液返排模型裂縫閉合期間壓裂液返排模型本章將主要運用流體力學(xué)和數(shù)值模擬的相關(guān)理論以及物質(zhì)平衡原理對壓裂液返排期間裂縫閉合的過程建立壓裂液返排的數(shù)學(xué)模型。給出壓裂液返排數(shù)學(xué)模型的數(shù)值解法，根據(jù)實例計算及方巖心流動實驗考察裂縫閉合過程中井底或井口壓力的變

36、化規(guī)律，對所選模型進行驗證和改進。3.1 裂縫閉合過程中模型的假設(shè)條件壓裂施工停泵后井底或井口壓力一般都會隨時間下降。這是因為地面停泵后，壓裂液的注入量為 0，裂縫內(nèi)的壓裂液在內(nèi)外壓差的作用下繼續(xù)濾失到地層，從而導(dǎo)致井底（井口）壓力下降。壓力下降是與濾失量以及裂縫寬度的變化緊密地聯(lián)系在一起的。所以，就可以根據(jù)閉合期間裂縫內(nèi)的流體體積平衡方程并結(jié)合巖石力學(xué)的理論來得到壓力的變化情況。在壓裂液返排過程中，裂縫參數(shù)的一些基本假設(shè)為：（1）停泵后縫中壓力短時間內(nèi)平衡，裂縫立即停止延伸；（2）設(shè)地層為線彈性體，層間無滑動，停泵后裂縫形狀呈矩形；（3）裂縫在閉合期間，井底裂縫的縫高和縫長不變，僅縫寬發(fā)生變

37、化；（4）支撐劑不影響裂縫的自由閉合。則根據(jù)注入階段和閉合期間裂縫內(nèi)的流體體積平衡原理，就可以建立分析停泵后壓力遞減規(guī)律。3.2 壓裂液返排的二維數(shù)學(xué)模型無論是一維濾失系數(shù)的計算方法還是修正的濾失系數(shù)的計算方法都只考慮了一維單相流動。本節(jié)在參考前人方法的基礎(chǔ)上，考慮返排壓裂液在地層中作二維流動和壓裂液為非牛頓流體的實際情況，建立壓裂液返排的二維模型，該模型綜合考慮了地層條件、油藏邊. v界條件和壓裂液性質(zhì)對壓裂液返排的影響，目的就是讓計算的結(jié)果對現(xiàn)場施工具有指導(dǎo)意義。.1 壓裂液從地層返排的數(shù)學(xué)模型壓裂液從地層返排的數(shù)學(xué)模型在壓裂施工過程及裂縫閉合的過程中，只要裂縫是張開著，在

38、裂縫壁面與油藏之間就會形成一定的壓差，而壓裂液的返排速度的大小又受裂縫與儲層之間的壓力梯度所控制36-39。在建立數(shù)學(xué)模型時分別考慮濾液在濾餅區(qū)和侵入?yún)^(qū)的輸運過程以及地層流體在儲層區(qū)的滲流過程，而且，假定地層流體驅(qū)替濾液采取的是活塞式，侵入?yún)^(qū)與儲層區(qū)交界處的流速連續(xù)。 壓裂液在侵入?yún)^(qū)的滲流壓裂液在侵入?yún)^(qū)的滲流考慮壓裂液非牛頓性質(zhì)對返排的影響，結(jié)合壓裂液或濾液的特性，將其假定為冪律流體。本文采用由 Teeuw 和 Hesselink 提出的適合于冪律流體的已修訂的達西定律： （3-1）1112n81312nnnnnKpnKL 式中，v為返排速度，m/s；為地層孔隙度，

39、無因次；為壓裂液流態(tài)指數(shù)，無因次；n地層滲透率，；為壓裂液稠度系數(shù)，Pasn；為壓降，Pa；L為多孔介質(zhì)的長K2mnKp度（侵入?yún)^(qū)長度） ，m。結(jié)合式（3-1） ，壓裂液在侵入?yún)^(qū)的運動方程如下： （3-2）deKVpG （）其中為啟動壓力梯度。將式（3-2）寫成分量形式xy,)G ( （3-3）dxedye()()xyKpVxKpVy . v和與和同號，說明由于啟動壓力梯度的作用使得滲流速度降低。 pxpyxy （3-4） 1112dend31288nnnnnnKpnK式中，e為濾液的有效黏度，Pas；為侵入?yún)^(qū)的滲透率，m2；為侵入?yún)^(qū)的孔隙度，dKd無因次。值得注意的是，由于壓裂液濾失造成了地

40、層損害，致使侵入?yún)^(qū)的孔隙度和滲透d率與儲層區(qū)的孔隙度和滲透率是不同的。dKK根據(jù)質(zhì)量守恒定律，可得到連續(xù)方程： （3-5） 0Vt考慮流體及多孔介質(zhì)的可壓縮性，流體密度變化及地層孔隙度變化的狀態(tài)方程如下： （3-6）0f01()cpp （3-7）001()cpp式中，和分別為流體及孔隙的壓縮系數(shù)，Pa-1；令，為綜合壓縮系fcccccfttc數(shù)，Pa-1。結(jié)合式（3-6）和（3-7） ，將式（3-3）代入式（3-5） ，則可得到如下的偏微分方程： （3-8）220dxf0fxf2e22y0df0fyf2e()1()()1()KVpppcppccxxxxVKpppcppccyyyy 將式（3-8

41、）中兩式相加，略去項，這樣可以得到：2fipcx （3-9）2etf0fxyd1()cpppcpppcxyKt令，由于不考慮油藏的非均質(zhì)性，則可令fff01()cccppttf01()cccpp. v，因此式（2-9）可以寫成：xy （3-10）2etfdcppppcxyKt方程（3-10）描述了濾液在侵入?yún)^(qū)的非穩(wěn)態(tài)滲流。 地層流體在儲層區(qū)的滲流地層流體在儲層區(qū)的滲流地層流體在儲層區(qū)的滲流一般被看作牛頓流體的流動。通?？蓪⑴ｎD流體的流動視為冪律流體流動的一種特殊情況。因此，在上述方程中設(shè)n=1 并采用原始地層的孔隙度和滲透率，按照相同的推導(dǎo)步驟，則可得到地層流體在儲

42、層區(qū)的滲流方程： （3-11）2tfcppppcxyKt式中，為地層流體的黏度，Pas；為原始地層的孔隙度，無因次；為原始地層K的滲透率，m2。 濾餅區(qū)的滲流濾餅區(qū)的滲流如果將濾餅視為滲透率，孔隙度的多孔介質(zhì)，并考慮壓裂液的非牛頓特性（假cKc定為冪律流體） ，結(jié)合方程（3-1） ，可得到壓裂液在濾餅區(qū)的運動方程： （3-12）12ccfwnc81231nnnnKppKnL式中，為濾餅區(qū)的滲透率，m2；為濾餅區(qū)的孔隙度，無因次；為裂縫中的壓cKcfP力，Pa；為濾餅與侵入?yún)^(qū)交界面的壓力，Pa；L為濾餅的厚度，m。wP在濾餅的形成過程中，可以認(rèn)為濾餅的體積與濾過的液體

43、體積成正比，則有： （3-13）dlLV 式中的比例系數(shù)可以寫成下式，是形成濾餅物質(zhì)的體積濃度，dsc （3-14）sdsc1(1) (1)CC. v將式（3-10）及式（3-11）代入式（3-9）中，整理后得到： （3-15）12sccfwnlscc312118nnnnCnppKVCnK從上式看到，濾餅兩端的壓力差與濾失速度及濾失體積所形成的濾餅厚度有關(guān)。上式可以寫成更簡單的形式： （3-16）fwclnppV其中，是冪律流體的濾餅因數(shù)，它的物理意義是指單位面積的單位濾失體積與達c到n冪的濾失速度所產(chǎn)生的濾餅壓降。它是取決于流體和濾餅性質(zhì)的一個綜合參數(shù)。但是，濾餅形成過程的復(fù)雜性常使濾餅壓降

44、與式（3-16）的計算結(jié)果有差異。 （3-17）12scccnscc312118nnCnKCnK實際上，在形成濾餅以前，有些更為微小的聚合物及地層微粒已經(jīng)進入到地層的內(nèi)部，稱之為內(nèi)濾餅。隨著內(nèi)濾餅的形成，外濾餅也就在時刻開始在縫面上沉積。因此，在許0t多地層中，內(nèi)濾餅的形成常常是外濾餅形成的基礎(chǔ)。可用式（3-16）來計算由于內(nèi)濾餅的存在而發(fā)生的壓力降，此時濾餅因數(shù)逐漸增長至某一上限值。隨著壓裂液向地層的濾失及濾餅的增厚，順著裂縫流過的壓裂液對己形成的外濾餅還有一種磨蝕的作用，阻止濾餅厚度的增加。此外，濾餅中的某些物質(zhì)與從它表面流過的壓裂液之間，由于濃度的差別，甚至還會產(chǎn)生從濾餅中擴散到壓裂液中

45、去的現(xiàn)象。由上述的對濾餅形成過程的分析可以看出，應(yīng)適當(dāng)修正濾餅因數(shù)，使之適應(yīng)濾餅形c成過程中的既有增厚又有削薄的情況。通過分析動濾失的實驗數(shù)據(jù)，用下式對進行修正：c （3-18）c1c2c00cc00tttttttt 式中，為修正后的濾餅因數(shù)，無因次；為外濾餅開始形成的時間，s； 為從開始c0tt濾失算起的時間，s；c1內(nèi)濾餅參數(shù)，無因次；c2為受侵蝕濾餅參數(shù)，無因次。此外，. v、c1、c2，4 個參數(shù)都是用實驗數(shù)據(jù)確定的。分析式（3-18） ，可以定性地看出，由0tca內(nèi)濾餅的影響使c1為正值，而剪切和擴散效應(yīng)使c2為負(fù)值。式（3-16）可以改寫成： （3-19）fwclnppV3.2.2

46、3.2.2 初始條件及邊界條件初始條件及邊界條件考慮濾液沿垂直裂縫壁面方向和裂縫方向的二維濾失（設(shè)沿裂縫方向為x方向，垂直裂縫方向為y方向） ，其初始條件及邊界條件如下： 初始條件初始條件返排過程開始前，可以根據(jù)有關(guān)的壓裂資料確定整個儲層的壓力，即停泵后儲層靜壓。則有， （3-20）itPP0對于裂縫閉合過程中的返排而言，返排的初始條件即為停泵時刻的壓力分布。此壓力分布情況由計算施工過程中的返排得到，也可以由停泵時的壓裂資料估算得到。 邊界條件邊界條件四條邊界與返排方向的三個區(qū)域保持一致，裂縫內(nèi)部的濾餅面、濾餅與地層間的裂縫面、侵入?yún)^(qū)與儲

47、層區(qū)間的交界面和儲層邊界。裂縫內(nèi)部的濾餅面處， （3-21）ffx LPP式中，為停泵時刻裂縫的長度，m。濾餅與地層間的裂縫面處，縫面壓力可依據(jù)式fL（3-19）算出。在返排過程中，侵入?yún)^(qū)與儲層區(qū)間的界面隨返排量的增加而移動。在活塞式驅(qū)動中，此邊界的移動是連續(xù)的。交界處的一邊是冪律流體在多孔介質(zhì)中的滲流，另一邊是牛頓液在多孔介質(zhì)中的滲流，故有：. v （3-22）dede()()()()KpKpxxKpKpyy動邊界動邊界動邊界動邊界儲層的邊界條件可根據(jù)實際情況設(shè)定，大致可分為兩類：封閉外邊界： （3-23）00eeyyxxypxp式中，和分別為矩形儲層的長度和寬度。eXeY定壓外邊界： （3

48、-24）eeppppeeyyxx此外，對于侵入?yún)^(qū)及儲層區(qū)的二維流動有如下的條件： （3-25）e0e0)(000XxtLypYyxpyx以上即為二維返排數(shù)學(xué)模型的控制方程、邊界條件和初始條件。3.3 模型的數(shù)值解法本部分將對上面所建立的壓裂液返排二維模型中的偏微分方程尋求合適的差分格式，并離散這些方程，給出數(shù)值求解的步驟，最后求得沿裂縫壁面的壓裂液返排速度及井底壓力變化40-45。.1 返排模型的返排模型的離散離散 網(wǎng)格系統(tǒng)的劃分網(wǎng)格系統(tǒng)的劃分. v以侵入?yún)^(qū)和儲層區(qū)為劃分對象（濾餅可以當(dāng)作邊界來處理，下面將作說明） ，由于裂縫的對稱性，沿半縫長取 1

49、/4 油藏為研究單元，網(wǎng)格劃分如圖 3-1 所示。在壓裂液返排過程中裂縫附近的壓力梯度大，而遠離裂縫的壓力梯度小，因此采用不均勻網(wǎng)格。圖 3-1 網(wǎng)格系統(tǒng)示意圖 侵入?yún)^(qū)滲流方程的離散由方程（3-7）可以看出，壓裂液的有效黏度是壓力梯度的函數(shù)，所以，對壓力進e行求解的時候構(gòu)成一個非線性系統(tǒng)。為了求解的方便，對進行線性化處理，并采用隱e式格式將方程（3-7）離散為如下形式： （3-26）ji,11ji,ji,1j1,iji,1ji,ji,1j1,iji,11ji,ji,gpepdpcpbpannnnn式（3-26）中，；1ii21ixxxi1i21ixxx；1jj21jyyyj1j2

50、1jyyy；21i21iixxx21j21jjyyy此外，t為時間步長；n表示計算過程中的第n時間步。附加邊界條件和初始條件，式（3-26）構(gòu)成五對角方程組。 儲層區(qū)滲流方程的離散儲層區(qū)滲流方程的離散方程（3-11）的類型與方程（3-10）的類型一樣，都為拋物型方程。但在方程（3-11）中，和對于同一儲層而言可看成與壓力梯度無關(guān)的物性參數(shù)，則方程的求解大為簡化。K同樣采用隱式格式將方程（3-11）離散為： （3-27）ji,11ji,ji,1j1,iji,1ji,ji,1j1,iji,11ji,ji,gpepdpcpbpannnnn式（3-27）中，. v其他的表示

51、項同式（3-26） 。同樣，附加邊界條件和初始條件后，式（3-27）構(gòu)成一個五對角方程組。 侵入?yún)^(qū)與儲層區(qū)之間動邊界的處理侵入?yún)^(qū)與儲層區(qū)之間動邊界的處理由于此邊界隨時間不斷向前推移，假設(shè)在第n時間計算步時，邊界上各點在網(wǎng)格系統(tǒng)中的坐標(biāo)為（，）則坐標(biāo)計算式為：nxMnyM （3-28）tvyytuxxnnnnnnMM1MMM1M（，）分別表示邊界上第M個坐標(biāo)點在n時刻，在x和y方向上的推移速度，nMunMv它們可以由動邊界的連續(xù)條件求得，由式（3-19）有， （3-29）MMMMMMMMMM11dMMMe(1 2)11dMMMe(1 2)nnnnnnnnnnnnnnx

52、xxxnxnnnnyxyxnyppppKKuxxppppKKvyy 裂縫壁面處的返排速度裂縫壁面處的返排速度 （3-30）wii,1dbie(1 2)1()nnnppKy.2 求解方法求解方法一般先求出三個區(qū)域的壓力分布然后再求返排速度。為了得到停泵時刻的壓力分布，需從壓裂液開始濾失時計算，在此期間不需要考慮裂縫的擴展，即縫長Lf在停泵以后是不再變化的。求解過程如下：第一步，計算的啟動。在停泵后，得到裂縫的縫長為Lf不再變化。在此縫長范圍內(nèi)用作為定壓邊界，以儲層原始壓力分布作為初始條件并結(jié)合儲層的邊界條件，用強隱式fp方法求解式（3-26）所構(gòu)成的五

53、對角方程組，從而得到各網(wǎng)格點的壓力。然后，在裂ji,p縫壁面處求出壓裂液的返排速度和，并計算出其推移距離，進而形成動邊界。MuMv第二步，隨著動邊界的推移，在計算的n時刻可以將計算區(qū)域化為兩大部分：第一部. v分為濾餅區(qū)，其壓力分布由式（3-19）計算；第二部分為侵入?yún)^(qū)和儲層區(qū)，通過動邊界的推移距離來判斷網(wǎng)格系統(tǒng)中的各點是處于侵入?yún)^(qū)還是儲層區(qū)，然后在這些點分別采用式（3-27）和式（3-26）提供的格式，結(jié)果形成兩組五對角矩陣方程。連同濾餅?zāi)Ｐ秃瓦吔鐥l件，可以形成關(guān)于壓力的耦合方程組。但是，濾餅的壓降取決于返排速度和返排量，而返排速度方程（3-30）又是地層壓力分布的函數(shù)，因此，就壓力而言整個

54、藕合方程組構(gòu)成一個非線性系統(tǒng)，線化后可采用逐次超松弛法來求解，即得到的值。npji,第三步，得到后，再根據(jù)方程（3-28）和（3-29）求出動邊界的推移距離，在給出njip,裂縫壓力和縫長Lf的條件下，重復(fù)第二步的計算方法，得到n+1 時間步的壓力分布fnp。以此重復(fù)，直至返排結(jié)束為止。1ji,np第四步，停泵后，縫長Lf不再變化，以停泵時刻的壓力分布為初始值，給定的條fp件下，重復(fù)第二步和第三步，則可得到不同時刻的壓力分布。根據(jù)壓力分布，可求解出各時刻的返排速度和返排量，由式（3-19）可以看出，濾失量跟的變化有關(guān)，可寫成fp的形式。f()Vf p根據(jù)上面的步驟便可以編制計算機程序進行計算。

55、后面便以實驗和實例計算的方式驗證模型的準(zhǔn)確程度，并加以修正。3.4 裂縫閉合時間及壓裂液返排量的確定.1 裂縫閉合時間確定裂縫閉合時間確定由于關(guān)井時間一般比較長，縫中平均壓力與井底壓力基本處于平衡狀態(tài)。對于裂縫而言，當(dāng)裂縫中的平均壓力值達到裂縫閉合壓力時，可認(rèn)為裂縫已基本上閉合。則裂縫閉合時有， （3-31）fcsc( )p tp式中，s為關(guān)井后裂縫中的平均壓力與井底壓力之比，無因次；為閉合壓力，cpMPa。前面己經(jīng)求得了關(guān)井期間井底壓力變化。根據(jù)上式，就可以在求出閉合壓力。在均質(zhì)、單層內(nèi)進行壓裂時，閉合壓力就等于壓裂層的最小主應(yīng)力。.2 壓裂液返排量的計算壓裂

56、液返排量的計算對于井筒及放噴油嘴中的一維流動，模型假設(shè)井筒中的返排液為牛頓流體，并忽略井. v筒的摩阻，排液油嘴的出口壓力約等于一個標(biāo)準(zhǔn)大氣壓。對于恒定流，由伯努利方程則可得到放噴油嘴中不同時刻返排液的速度44： （3-32）5 . 05 . 0225 . 03 1 . 0)()1 (101tpRr由式（3-29）又可導(dǎo)出不同時刻的返排流量Q和累積返排量，計算式分別如下：outV （3-33）5 . 05 . 0225 . 023 1 . 0)()1 (101tpRrrQ （3-34）2320 50.50.5201 10(1) ( )0.1dt.outrVrp ttR上面的式子中，為返排液的密

57、度，kg/m3；為局部阻力系數(shù)，無量綱，在這里取，0.5； ，為放噴油嘴和井筒的半徑，m；為返排過程中井口的壓力，Pa；其與井底rR)(tp壓力的關(guān)系為，其中，為井筒靜液柱壓力，Pa。fh( )( )p tp tphp由式（3-30）可知，積分式中關(guān)于的函數(shù)為非線性函數(shù)，為了求解返排過程中不)(tp同時刻的累積返排量，需采用復(fù)合梯形公式的形式將積分式進行離散，則有：outV （3-35）out0fefh1 ( )0.1(0.1)2(0.1)2nniitVp tppppCC 其中： （3-36）5 . 0225 . 023)1 (105 . 0RrrCC式（3-32）中，t為

58、時間步長，與前面計算井口壓降和濾失量時的取法一致；上標(biāo)i表示時刻；表示剛返排時井口壓力，如果返排是從停泵后立即進行的，那么0t)(0tp等于瞬時停泵壓力 ISIP。)(0tp從式（3-33）可以看出，在返排過程中，是隨井底壓力而變化的。在時刻，outVfpnt至?xí)r刻的井底壓力值已求出，則僅是的函數(shù)，公式表示如下：1t1ntoutVfp （3-37）outf()nnVf p則根據(jù)注入階段和閉合期間裂縫內(nèi)的流體體積平衡原理，就可以建立分析停泵后壓力遞減規(guī)律的擬三維模型。根據(jù)體積平衡原理，壓裂液返排時，裂縫體積的變化量應(yīng)該等于. v返排時刻起的壓裂液總濾失量與返排量之和，公式表示如下： （3-38）

59、outlf0ff)()()(VtVtVtVV式中，和分別為剛返排和返排過程中某一時刻裂縫的體積，m3。)(0ftV)(ftV.3 停泵后裂縫體積變化量的計算停泵后裂縫體積變化量的計算在裂縫閉合期間，由于縫長和縫高不變，僅有縫寬變化。設(shè)剛停泵時井口壓力為ISIP，則停泵后某時刻裂縫體積的變化量為： （3-39）22ffhf2(1)ISIP( )VLpp t HE式中，ISIP 為停泵時刻井口壓力，Pa ；為井筒中靜液柱產(chǎn)生的壓力，Pa；hp在裂縫參數(shù)和地層參數(shù)己經(jīng)給定的情況下，上式中的為的函數(shù)，表示為：fV)(ftp （3-40）)()(f2ftpftV由以上的分析可以看出，和都

60、是的函數(shù)，結(jié)合第二章對濾失量的計算，outVfV)(ftptV在每一計算步內(nèi)對式（3-35）中的壓力項迭代求解，就可以得到停泵后不同時刻的井底壓力，進而得到井口的壓力變化，即。hf)()(ptptp3.5 實例計算與分析根據(jù)前文介紹的關(guān)井期間井底或井口壓力變化、裂縫閉合壓力以及裂縫閉合時間的計算方法，對實例進行計算，以得出返排過程中的井口壓力變化，并對計算結(jié)果進行分析?；居嬎銋?shù)如表 3-1，計算結(jié)果如圖 3-2 所示。表 3-1 壓裂液返排模擬計算中需要輸入的基本參數(shù)儲層參數(shù)參數(shù)值壓裂參數(shù)參數(shù)值有效厚度(m)19裂縫高度(mm)40泊松比(無因次)0.16瞬時停泵壓力(MPa)13.5楊氏