范云鵬 周志軍 李長(zhǎng)春

摘 ?????要： 在經(jīng)過幾十年的注水開發(fā)后，國(guó)內(nèi)大多數(shù)油藏都進(jìn)入了特高含水期，這使得油藏的進(jìn)一步開采更加困難，尤其是經(jīng)歷了高壓注水和循環(huán)注水等過程后，儲(chǔ)層壓力、地層孔隙結(jié)構(gòu)和流體參數(shù)等不斷的發(fā)生變化，這導(dǎo)致儲(chǔ)層多孔介質(zhì)發(fā)生不規(guī)律的變形，而這種變形同時(shí)也會(huì)影響地層流體的滲流作用，致使特高含水儲(chǔ)層的流固耦合作用愈發(fā)明顯。結(jié)合流固耦合理論在特高含水期油藏開發(fā)中的研究現(xiàn)狀和發(fā)展歷程進(jìn)行綜述和分析，并對(duì)該理論未來的發(fā)展趨勢(shì)進(jìn)行展望。

關(guān) ?鍵 ?詞：特高含水期;流固耦合;油藏?cái)?shù)值模擬

中圖分類號(hào)：TE341 ??????文獻(xiàn)標(biāo)識(shí)碼： A ??????文章編號(hào)： ?1671-0460（2019）03-0575-04

Abstract： With the water injection development of oil field for several decades， most of the domestic oil reservoirs have entered an ultra-high water cut stage， which makes further exploitation of the oil reservoirs more difficult. Especially after the process of high pressure water injection and circulating water injection， the reservoir pressure， stratigraphic pore structures and fluid parameters constantly change， which leads to irregular deformation of reservoir porous medium. At the same time， this kind of deformation will also affect the seepage effect of the formation fluid， which makes the fluid-solid coupling of ultra-high water cut reservoir more obvious. In this article， the research status and development history of the fluid-solid coupling theory in the development of oil reservoirs with ultra-high water cut stage were summarized and analyzed， and the future development trend of the theory was forecasted.

Key words： Ultra-high water cut stage; Fluid-solid coupling; Oil reservoir numerical simulation

當(dāng)油藏進(jìn)入高含水期以后，由于可采儲(chǔ)量下降，剩余油的分布分散，使得采油效率急劇下降，并且耗水量增大，注采井組的井況日漸惡化，開發(fā)效益越來越差。但是，從儲(chǔ)層豐度以及油田總體產(chǎn)量來看，特高含水期的老油田依舊是整個(gè)區(qū)塊的開采主體，所以在特高含水油藏的開采方面不容小覷[1]。因此，對(duì)處于特高含水期油藏的挖潛，是關(guān)系到國(guó)家能源發(fā)展的重要因素，只有不斷的嘗試提高采收率，才能更充分的發(fā)揮老油田的經(jīng)濟(jì)優(yōu)勢(shì)。

油田的常規(guī)開采中一般都會(huì)經(jīng)歷注水開發(fā)的過程，凡是涉及到油藏?cái)?shù)值模擬的過程就不得不考慮儲(chǔ)層中流體的滲流規(guī)律。從經(jīng)典的滲流理論，到現(xiàn)在的流固耦合理論，都繞不開流體與巖石顆粒之間的相互作用關(guān)系。流體在儲(chǔ)層孔隙中的滲流會(huì)涉及到不同方面的理論，首先必須考慮的是油藏注水過程中流體壓力的變化，進(jìn)而會(huì)引起巖石顆粒的運(yùn)移以及孔隙雜質(zhì)的流動(dòng)。在原有地層壓力平衡被打破以后，巖石骨架會(huì)趨向于下一個(gè)穩(wěn)態(tài)，這一過程中就會(huì)發(fā)生多孔介質(zhì)的變形，孔隙結(jié)構(gòu)的重新分布直接影響著儲(chǔ)層物性和流體的滲流規(guī)律;其次，由于儲(chǔ)層孔隙結(jié)構(gòu)的重新分布，儲(chǔ)層空間發(fā)生變化引起流體壓力的重新分布，這反過來用于巖體的孔隙和吼道上，使得巖石骨架在受力條件下趨于新的平衡，產(chǎn)生了又一次的巖石骨架變形，直到它們之間達(dá)到一個(gè)動(dòng)態(tài)平衡狀態(tài)，這就是流固耦合作用。因此，只要涉及到流體在巖石體儲(chǔ)層中的滲流問題，就不可避免的考慮流固耦合作用，這也稱之為儲(chǔ)層開采過程中應(yīng)力場(chǎng)和滲流場(chǎng)的耦合效應(yīng)。本文將從流固耦合理論在特高含水期油藏開發(fā)中的研究現(xiàn)狀和發(fā)展趨勢(shì)展開論述與分析。

1 ?流固耦合理論在特高含水期油藏開發(fā)中的國(guó)內(nèi)外研究現(xiàn)狀

在石油開采的相關(guān)工藝中，伴隨著油氣的開采和運(yùn)移，都會(huì)涉及到流體滲流和巖石骨架或者地層裂縫之間的強(qiáng)相互耦合作用，這種耦合作用會(huì)明顯的影響油藏的開采以及相關(guān)施工（鉆井、固井、完井、增產(chǎn)壓裂改造等）難度。尤其在注水開發(fā)的油藏中，當(dāng)油藏到達(dá)特高含水期時(shí)，油藏內(nèi)的流體流動(dòng)頻繁，對(duì)巖石骨架和孔隙結(jié)構(gòu)產(chǎn)生難以估量的影響。因此針對(duì)特高含水期的油藏開發(fā)需要著重對(duì)考慮流固耦合問題。國(guó)內(nèi)外學(xué)者對(duì)流固耦合問題一直都十分看重，并且已經(jīng)成果斐然。

1.1 ?國(guó)外流固耦合研究現(xiàn)狀

在石油天然氣開發(fā)領(lǐng)域，國(guó)外學(xué)者在前人的研究成果之上做了大量工作，關(guān)于流固耦合的理論不斷豐富和完善。

M.Mainguy[2] 研究了油藏儲(chǔ)層中流體流動(dòng)與巖石之間的耦合方程，通過巖石本構(gòu)方程、儲(chǔ)層應(yīng)力方程等做了進(jìn)一步計(jì)算整合，編輯了新的完全耦合模擬器可以實(shí)現(xiàn)與傳統(tǒng)儲(chǔ)層模擬器之間的信息互換。N.Khalili[3] 提出了一種控制雙孔隙彈性介質(zhì)中非等溫流動(dòng)和變形的偏微分方程，用以分析雙孔彈性介質(zhì)中熱流體的注入所產(chǎn)生的全耦合相關(guān)規(guī)律。D.ANDREW S.REES[4]對(duì)熱流體注入多孔介質(zhì)中產(chǎn)生的局部非平衡熱效應(yīng)做出研究，提出溫度對(duì)巖石多孔介質(zhì)中流體的影響局部非平衡方程。Y.W.Kwon Mech[5]將Bltzmann法應(yīng)用于不規(guī)則形狀流體域，對(duì)流體流動(dòng)方程進(jìn)行求解，開發(fā)出新的流固耦合作用分析程序，能夠解決每個(gè)時(shí)間步長(zhǎng)的耦合瞬間流體流動(dòng)和結(jié)構(gòu)力學(xué)參數(shù)的精確計(jì)算，并給出2維和3維條件下流固耦合模擬實(shí)例。R.Gelet[6]針對(duì)雙重孔隙結(jié)構(gòu)的裂縫性油藏采用局部分析非平衡熱流體注入條件下的熱-流-固耦合模型，運(yùn)用有限元近似位移矢量描述巖石骨架變形，并用軟件模擬參數(shù)變化。

Pan Peng-Zhi[7]應(yīng)用模擬多組分流體流動(dòng)和熱傳遞的程序TOUGH2結(jié)合模擬巖石非線性和不連續(xù)地質(zhì)力學(xué)特性的程序RDCA，集成為TOUGH-

RDCA模擬器用于多相流體存在的流固耦合應(yīng)力場(chǎng)和滲流場(chǎng)的模擬分析。N.Guy[8]就蒸汽輔助重力泄油（SAGD）過程中由應(yīng)力變化導(dǎo)致儲(chǔ)層內(nèi)滲透率和孔隙度的變化，建立了熱流體流動(dòng)與地層巖石骨架之間的耦合方程，編寫有效的熱-流-固耦合模擬數(shù)值程序，能更好的實(shí)現(xiàn)油藏模擬過程中自適應(yīng)網(wǎng)格細(xì)化過程，從而對(duì)蒸汽驅(qū)動(dòng)前沿做出更精確的描述。Zhong Ruizhi[9]研究了水力壓裂過程中高壓流體和巖石骨架之間的耦合作用，通過建立一個(gè)全耦合有限元模型來模擬多個(gè)水力壓裂過程，并且針對(duì)性的模擬了高壓流體耦合作用下裂縫幾何形狀與應(yīng)力變化之間的關(guān)系。

1.2 ?國(guó)內(nèi)流固耦合發(fā)展?fàn)顩r

國(guó)內(nèi)在流固耦合滲流理論研究方面起步較晚，但近年來發(fā)展較快，并且由原來的單一靠現(xiàn)有的成熟商業(yè)化軟件模擬到現(xiàn)階段的自主研發(fā)軟件模擬。

2001年黎水泉等在已有變形介質(zhì)滲流模型基礎(chǔ)上考慮了有效應(yīng)力對(duì)滲流參數(shù)的影響，推導(dǎo)了非線性滲流模型的有限元計(jì)算格式，并編寫了流固耦合計(jì)算程序[10]。

2002年劉建軍等針對(duì)低滲透儲(chǔ)層的滲流規(guī)律，充分考慮啟動(dòng)壓力梯度，推導(dǎo)計(jì)算了低滲條件下的耦合數(shù)學(xué)模型，在已有黑油模型基礎(chǔ)上編寫了適合低滲油藏流固耦合計(jì)算軟件[11]。

2003年周志軍做了低滲條件下油水兩相滲流的綜合分析研究，建立相關(guān)數(shù)學(xué)方程并做了數(shù)值求解，在GALERKIN有限元理論的基礎(chǔ)上設(shè)計(jì)數(shù)值模擬計(jì)算軟件并對(duì)低滲透油藏有關(guān)滲流參數(shù)進(jìn)行模擬計(jì)算[12]。

2004年李勇通過建立孔隙度和滲透率模型研究巖石骨架變形對(duì)其地層參數(shù)的影響，進(jìn)一步分析壓裂體系中流固耦合作用，編制整體壓裂動(dòng)態(tài)開發(fā)模擬軟件[13]。

2005年文成楊采用Drucker-Prager原則對(duì)氣藏儲(chǔ)層的彈塑性做了描述，通過分析應(yīng)力應(yīng)變與巖石骨架之間的關(guān)系建立流固耦合全耦合模型，基于FEPG程序開發(fā)了適合氣藏的流固耦合計(jì)算程序[14]。

2007年孫輝、李玉坤等人研究了稠油熱采過程中注蒸汽對(duì)地層力學(xué)特性的影響，進(jìn)一步完善了熱-流-固耦合模型，并通過GALERKIN有限元法中的全隱式解法做了求解和模擬，得到單項(xiàng)流中流體的壓力分布規(guī)律[15，16]。

2009年孫峰等結(jié)合油藏生產(chǎn)過程中流體對(duì)巖石骨架的影響，建立了與地層應(yīng)力應(yīng)變相關(guān)的穩(wěn)定性評(píng)價(jià)模型，定量分析井底附近儲(chǔ)層變形及應(yīng)力變化情況[17]。

2010年賈善坡等人以油氣藏中基本滲流規(guī)律為基礎(chǔ)，建立了相關(guān)地層參數(shù)的動(dòng)態(tài)演變方程，推導(dǎo)油氣儲(chǔ)層中流固耦合全耦合動(dòng)態(tài)模型，并在一組注采井網(wǎng)中加以實(shí)驗(yàn)驗(yàn)證[18]。

2011年楊軍征等人提出了控制體積流量法的計(jì)算思路，有效的減少了流固耦合程序編寫及運(yùn)算的難度。姚軍、單嫻等人針對(duì)裂縫性油藏開發(fā)過程中存在的流固耦合問題，建立了油水兩相滲透率張量和滲流壓力等參數(shù)的計(jì)算模型，運(yùn)用GALERKIN加權(quán)殘量法計(jì)算后模擬驗(yàn)證[19-21]。

2012年張?jiān)实热酸槍?duì)目前油藏?cái)?shù)值模擬技術(shù)精度差、計(jì)算速度慢等問題，提出分區(qū)域單個(gè)計(jì)算后并行模擬的方法，實(shí)現(xiàn)了數(shù)值模擬高效、精確和可移植性的并行編程計(jì)算[22]。

2013年盛茂等人對(duì)頁巖氣藏的流固耦合方程做出改進(jìn)，通過考慮氣體滲流場(chǎng)對(duì)巖石骨架形變影響以及各種滲流機(jī)理的作用，綜合建立了滲流模型并用有限元法和隱式解法對(duì)模型做出求解，編輯了數(shù)值模擬軟件[23]。

2015年張宏學(xué)博士對(duì)煤系頁巖儲(chǔ)層注氣條件下的流固耦合規(guī)律做了研究，分別建立了基于雙重介質(zhì)的孔隙度和滲透率模型并用物理實(shí)驗(yàn)加以驗(yàn)證，得出頁巖儲(chǔ)層流固耦合控制方程，利用COMSOL Multiphysics軟件模擬了基質(zhì)和裂縫系統(tǒng)中的耦合作用多頁巖氣采收率的影響[24]。

2016年蘇玉亮、康永尚等從多重介質(zhì)耦合方面考慮頁巖氣的滲流規(guī)律，建立動(dòng)態(tài)耦合模型后推導(dǎo)出方程的半解析解，并通過物理模擬實(shí)驗(yàn)中的CT成像技術(shù)加以分析驗(yàn)證[25，26]。

我國(guó)的油田進(jìn)過多年的注水開發(fā)，油藏大都進(jìn)入高含水甚至特高含水期，通過以往學(xué)者在流固耦合理論和實(shí)踐方面的研究，可以為特高含水期油藏開發(fā)借鑒思路，豐富研究手段，為開發(fā)高含水時(shí)期的剩余油提供理論基礎(chǔ)和技術(shù)手段。

2 ?流固耦合理論在特高含水期油藏開發(fā)中的發(fā)展趨勢(shì)

未來流固耦合理論在特高含水期油藏開發(fā)中的應(yīng)用發(fā)展趨勢(shì)主要表現(xiàn)在以下幾個(gè)方面。

2.1 ?由簡(jiǎn)單發(fā)展到越來越復(fù)雜

以往在受計(jì)算機(jī)技術(shù)的限制，數(shù)值模擬涉及該方面的運(yùn)算一直受到制約，隨著近年來試驗(yàn)測(cè)試和計(jì)算能力的提升，耦合問題的研究從理論到實(shí)踐都有了一個(gè)長(zhǎng)足的發(fā)展。從原來的最簡(jiǎn)單了單孔單滲模型轉(zhuǎn)向更為負(fù)雜的雙孔連續(xù)性介質(zhì)模擬以及后來連續(xù)性介質(zhì)或非連續(xù)性孔隙網(wǎng)格模型的發(fā)展。從之前單一固相的介質(zhì)微小彈性變形轉(zhuǎn)為更貼合實(shí)際生產(chǎn)需求的非線性有限變形問題，其中涉及到了巖體結(jié)構(gòu)的彈塑性以及流變和蠕變理論。不但從滲流理論方面尋找突破，還可以從巖石力學(xué)性能性質(zhì)方面考慮[27]。

2.2 ?由定性研究向精確化的定量研究發(fā)展

定量科學(xué)研究是對(duì)事物由感性（定性）認(rèn)識(shí)上升到理性（定量）的一個(gè)重要手段。定量科學(xué)研究是一個(gè)典型范式，也就是基于數(shù)學(xué)模型的研究方法，常見的流固耦合問題亦是如此，再針對(duì)模型的數(shù)值求解過程中有兩種方法：有限差分法和有限元法 [28]。人們對(duì)于流固耦合理論最開始的認(rèn)知只是把一些影響因素考慮在內(nèi)，對(duì)特高含水期油藏的內(nèi)部變化有一個(gè)定性的了解，對(duì)于變化的大小、幅度并沒有概念，而現(xiàn)在隨著計(jì)算機(jī)技術(shù)的日趨進(jìn)步，相關(guān)理論日益完善，人們對(duì)于特高含水時(shí)期油藏開發(fā)中流固耦合問題，也由定性認(rèn)知提升到日益追求精確化的定量研究，而且在未來精確度也會(huì)逐漸的提升。

2.3 ?模擬軟件逐漸成熟化和個(gè)性化

對(duì)于特高含水期油藏中的流固耦合理論的研究，對(duì)于該理論的內(nèi)涵及驗(yàn)證，最開始只是停留在方程推導(dǎo)，建立耦合方程，這只是純理論階段。后來人們利用成熟的商業(yè)化軟件對(duì)流固耦合的過程進(jìn)行模擬，從而驗(yàn)證耦合項(xiàng)和耦合方程的合理性。而未來的發(fā)展趨勢(shì)則是由原來的單一靠現(xiàn)有的成熟商業(yè)化軟件模擬到現(xiàn)階段的自主研發(fā)軟件模擬。

3 ?結(jié)束語

在資源日益匱乏的大環(huán)境下，充分開發(fā)特高含水期的剩余油顯得尤為重要，而在該過程中考慮流固耦合問題更是有利于加深對(duì)特高含水期油藏內(nèi)流體運(yùn)移和巖石骨架變形的認(rèn)識(shí)。而對(duì)于長(zhǎng)期注水開發(fā)的老油田，尤其是受到高強(qiáng)度注水沖刷以后，儲(chǔ)層地質(zhì)條件變得十分復(fù)雜，以往的開采理論對(duì)現(xiàn)有狀況的指導(dǎo)效果收效甚微，只有從流固耦合滲流機(jī)理方面取得成效，才能精確的預(yù)測(cè)地層耦合規(guī)律、剩余油分布、地層物性參數(shù)變化等，提出更好的優(yōu)化方案、增產(chǎn)措施。

參考文獻(xiàn)：

[1]任剛. 基于相控模型的精細(xì)數(shù)值模擬技術(shù)研究[D]. 東北石油大學(xué)， 2011.

[2]Mainguy M， Longuemare P. Coupling Fluid Flow and Rock Mechanics： Formulations of the Partial Coupling Between Reservoir and Geomechanical Simulators[J]. Oil & Gas Science & Technology， 2002， 57（4）：355-367.

[3]Khalili N， Selvadurai A P S. A fully coupled constitutive model for thermo‐hydro‐mechanical analysis in elastic media with double porosity[J]. Geophysical Research Letters， 2003， 30（24）：665-678.

[4]Rees D A S， Bassom A P， Siddheshwar P G. Local thermal non-equilibrium effects arising from the injection of a hot fluid into a porous medium[J]. Journal of Fluid Mechanics， 2008， 594 （594）： 379-398.

[5]Kwon Y W. Multi-scale modeling of mechanical behavior of polycrystalline materials[J]. Journal of Computer-Aided Materials Design， 2004， 11（1）：43-57.

[6]Gelet R， Loret B， Khalili N. A thermo‐hydro‐mechanical coupled model in local thermal non‐equilibrium for fractured HDR reservoir with double porosity[J]. Journal of Geophysical Research Solid Earth， 2013， 117（B7）.

[7]Feng X T， Pan P Z， Zhou H. Simulation of the rock microfracturing process under uniaxial compression using an elasto-plastic cellular automaton[J]. International Journal of Rock Mechanics & Mining Sciences， 2006， 43（7）：1091-1108.

[8]Guy N， Enchery G， Renard G. Numerical Modeling of Thermal EOR： Comprehensive Coupling of an AMR-Based Model of Thermal Fluid Flow and Geomechanics[J]. Oil & Gas Science & Technology， 2013， 67（6）：1019-1027.

[9]Zhong R， Bao J， Fathi E. Fully coupled finite element model to study fault reactivation during multiple hydraulic fracturing in heterogeneous tight formations[J]. Society of Petroleum Engineers， 2014：418-429.

[10]黎水泉， 徐秉業(yè)， 段永剛. 裂縫性油藏流固耦合滲流[J]. 計(jì)算力學(xué)學(xué)報(bào)， 2001， 18（2）：133-137.

[11]劉建軍， 劉先貴， 胡雅礽，等. 低滲透儲(chǔ)層流-固耦合滲流規(guī)律的研究[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2002， 21（1）：37-41.

[12]周志軍. 低滲透儲(chǔ)層流固耦合滲流理論及應(yīng)用研究[D]. 大慶石油學(xué)院， 2003.

[13]李勇. 考慮流固耦合效應(yīng)的整體壓裂數(shù)值模擬研究[D]. 西南石油學(xué)院， 2004.

[14]文成楊. 雙重介質(zhì)氣藏流固耦合數(shù)值模擬研究[D]. 西南石油學(xué)院， 2005.

[15]孫輝， 李兆敏， 焦玉勇. 稠油油藏?zé)?流體-力學(xué)耦合模型研究及應(yīng)用[J]. 巖土力學(xué)， 2007， 28（12）：2560-2564.

[16]李玉坤， 曲曉建， 毛勇. 復(fù)雜邊界油藏滲流壓力分布的有限元法模擬[J]. 中國(guó)石油大學(xué)勝利學(xué)院學(xué)報(bào)， 2007， 21（4）：1-3.

[17]孫峰， 薛世峰， 葛洪魁，等. 基于流固耦合理論的疏松砂巖地層穩(wěn)定性評(píng)價(jià)[J]. 西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2009， 24（2）：9-12.

[18]賈善坡， 王強(qiáng)， 姚華彥. 可變形儲(chǔ)層注采過程中滲流場(chǎng)與應(yīng)力場(chǎng)動(dòng)態(tài)耦合分析[J]. 長(zhǎng)江大學(xué)學(xué)報(bào)（自科版）， 2010， 7（1）：104-107.

[19]楊軍征， 汪緒剛， 王瑞和，等. 基于有限元法的油藏開發(fā)數(shù)值模擬[J]. 新疆石油地質(zhì)， 2011， 32（1）：54-56.

[20]單嫻， 姚軍. 基于滲透率張量的各向異性油藏兩相滲流數(shù)值模擬[J]. 中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2011， 35（2）：101-106.

[21]單嫻， 姚軍. 基于離散裂縫模型的裂縫井滲流壓力場(chǎng)分析[J]. 油氣地質(zhì)與采收率， 2011， 18（3）：67-69.

[22]張?jiān)剩?袁向春. 油藏?cái)?shù)值模擬有限元并行計(jì)算方法研究[J]. 微計(jì)算機(jī)信息， 2012（1）：39-41.

[23]盛茂， 李根生， 黃中偉，等. 頁巖氣藏流固耦合滲流模型及有限元求解[J]. 巖石力學(xué)與工程學(xué)報(bào)， 2013， 32（9）：1894-1900.

[24]張宏學(xué). 頁巖儲(chǔ)層滲流—應(yīng)力耦合模型及應(yīng)用[D]. 中國(guó)礦業(yè)大學(xué)， 2015.

[25]蘇玉亮， 盛廣龍， 王文東，等. 頁巖氣藏多重介質(zhì)耦合流動(dòng)模型[J]. 天然氣工業(yè)， 2016， 36（2）.

[26]康永尚， 鄧澤， 王紅巖，等. 流-固耦合物理模擬實(shí)驗(yàn)及其對(duì)頁巖壓裂改造的啟示[J]. 地球科學(xué)-中國(guó)地質(zhì)大學(xué)學(xué)報(bào)， 2016， 41（8）：1376-1383.

[27]吳林強(qiáng)， 劉成林， 李冰，等. 應(yīng)力場(chǎng)數(shù)值模擬與油藏有利區(qū)預(yù)測(cè)--以松遼盆地乾安地區(qū)歸字井青三段為例[J]. 地質(zhì)力學(xué)學(xué)報(bào)， 2014（4）：339-351.

[28]董平川， 韓德金， 牛彥良，等. 油藏多相滲流的面向?qū)ο笥邢拊绦蛟O(shè)計(jì)[J]. 巖土力學(xué)， 2009， 30（4）：1115-1121.