馮 曦 彭 先 李隆新 楊學鋒 王 娟 李 騫 張 春 鄧 惠

中國石油西南油氣田公司勘探開發(fā)研究院

0 引言

四川盆地內(nèi)的碳酸鹽巖氣藏在開發(fā)過程中產(chǎn)地層水的比例較高。近60年來，得益于持續(xù)開展氣藏開發(fā)治理水侵影響的研究探索，形成了較為完整的技術(shù)體系。然而，近年卻出現(xiàn)了以下新情況：①深層海相碳酸鹽巖氣藏勘探開發(fā)獲得突破，大型裂縫—孔洞型氣藏數(shù)量增加，過去基于中小型裂縫型、裂縫—孔隙型氣藏積累的治水技術(shù)和經(jīng)驗不適用的情況時有發(fā)生；②勘探開發(fā)一體化工作模式普及，開發(fā)前期評價和產(chǎn)能建設(shè)節(jié)奏加快，早期預判水侵影響而制訂針對性對策、降低潛在風險的需求增強；③環(huán)境保護和氣藏開發(fā)效益要求進一步提高，一些僅考慮氣藏技術(shù)采收率的強排水開發(fā)方式應用受限。針對上述問題，近10年中國石油西南油氣田公司組織開展了大量研究，逐漸形成新的技術(shù)成果和認識。本文主要介紹機理性研究方面的代表性成果。

掌握水侵規(guī)律、針對性制訂優(yōu)化開發(fā)的對策是改善有水氣藏開發(fā)效果的關(guān)鍵[1-2]。氣藏水侵活躍性及其對開發(fā)的影響程度與水體儲量、氣水接觸關(guān)系、水侵通道分布、水區(qū)和氣藏內(nèi)部儲層物性及非均質(zhì)性、氣藏原始地層壓力等地質(zhì)因素相關(guān)，同時還與生產(chǎn)井位置及打開層段、井網(wǎng)部署、單井配產(chǎn)、氣藏采氣速度等開發(fā)因素相關(guān)[3]。氣藏描述和動態(tài)監(jiān)測分析是準確認識上述關(guān)系、支撐氣藏開發(fā)治水決策的主要技術(shù)手段。

經(jīng)過長期的理論研究與實踐探索，形成了有水氣藏開發(fā)治理水侵影響的系列技術(shù)，在地質(zhì)和氣藏工程方面主要包括氣水關(guān)系描述[4-5]、水侵動態(tài)分析預測[6-11]、治水對策優(yōu)化[12-14]等特色技術(shù)，已經(jīng)廣泛應用并取得較好效果。然而，仍有一些復雜問題長期困擾從事氣藏開發(fā)治水研究的工作者，主要表現(xiàn)在兩個方面：①不同類型氣藏多樣化的非均質(zhì)性難以用統(tǒng)一的描述理論或模型來表征，一些實際情況超出了傳統(tǒng)分析方法的適用范圍，直接影響氣藏水侵分析預測結(jié)果的可靠性；②制訂氣藏開發(fā)治水對策時，常常不具備或不完全具備精細氣藏描述、專項動態(tài)分析所需的基礎(chǔ)資料條件，在開發(fā)早期預判水侵的影響極其困難。

強化對不同類型氣藏水侵機理的研究和典型氣藏案例類比分析是解決上述問題的有效途徑，盡管這方面也有較多研究成果，主要包括水封氣及水侵后殘余氣分析[15-18]、束縛水與可動水評價[19]、微觀滲流過程可視化研究[20]、氣藏水侵仿真實驗模擬[21]、典型氣藏水侵特征剖析[22]等，但由于問題的復雜性，相關(guān)方法有必要進一步拓展，一些認識尚需進一步深化，尤其在非均質(zhì)性較強的碳酸鹽巖氣藏開發(fā)方面更顯著。

1 碳酸鹽巖儲層三重介質(zhì)特征定量描述

1.1 儲層分類定性方法的優(yōu)點和局限性

參照國家標準《GВ/T 26979—2011 天然氣藏分類》[23]，氣藏儲滲空間類型分為孔隙型、裂縫—孔隙型、裂縫—孔洞型、孔隙—裂縫型和裂縫型5種類型（表1）。這種分類方法把握了不同類型儲層的本質(zhì)特征，被廣泛應用，對氣藏開發(fā)研究工作幫助很大。

表1 氣藏儲滲空間類型特征表

定性分析方法的不足是儲滲類型劃分界限較為模糊，同屬一種類型則不再加以區(qū)分。因此難以更細致地甄別儲層的特征差異。以近10年我國天然氣開發(fā)熱點領(lǐng)域之一——深層碳酸鹽巖氣藏為例，儲層儲滲空間類型普遍屬裂縫—孔洞型，但不同氣藏儲滲特征差異較大，無論是巖心靜態(tài)特征還是試井曲線動態(tài)反映均可見顯著差異（圖1）。因此，定性的分類評價方法已不能滿足需求，有必要深入研究可區(qū)分儲滲特征差異的定量評價方法。

1.2 孔、洞、縫儲集空間對總孔隙度的貢獻定量描述

碳酸鹽巖儲層大都不同程度發(fā)育孔、洞和裂縫，雖然因孔、洞、縫尺度和搭配關(guān)系不同產(chǎn)生千差萬別的儲滲特征，但整體上仍可歸類于廣義的三重介質(zhì)儲層。一些情況下某類儲集空間欠發(fā)育，表現(xiàn)為雙重介質(zhì)甚至單一介質(zhì)特征，然而三重介質(zhì)的統(tǒng)一描述方法仍然適用，只是相關(guān)特征參數(shù)不同。

彈性儲容比和竄流系數(shù)是描述多重介質(zhì)儲層儲滲特征的重要參數(shù)[24-25]，過去多用于試井分析研究。受多種因素影響，在很多情況下計算彈性儲容比和竄流系數(shù)所需的特征數(shù)據(jù)在試井曲線上反映不充分，導致上述參數(shù)較難準確求取。因此依靠試井解釋評價多重介質(zhì)儲滲特征差異的方法推廣應用受到一定程度的限制。

圖1 深層碳酸鹽巖氣藏裂縫—孔洞型儲層儲滲特征差異示例圖

近年來數(shù)字巖心技術(shù)迅速發(fā)展，高分辨率CT掃描和三維數(shù)字重構(gòu)建模、微觀儲集空間圖像識別與特征信息自動提取、數(shù)字巖心流動模擬等技術(shù)逐步成熟[26-28]。目前，全直徑巖心數(shù)字化處理分析結(jié)合孔隙度測定實驗，已能準確計算總孔隙度中孔隙、洞穴、裂縫儲集空間各自的貢獻比例，可獲得類似于彈性儲容比的結(jié)果，為三重介質(zhì)儲層特征的定量描述提供了新技術(shù)手段。

圖2為反映總孔隙度中孔、洞、縫儲集空間所占比例的三元圖版，圖2中的數(shù)據(jù)來源于四川盆地磨溪寒武系龍王廟組氣藏（以下簡稱龍王廟組氣藏）多口井的數(shù)字巖心分析結(jié)果。新方法有效回避了儲層儲滲空間類型劃分界限難以量化的問題，能從廣義三重介質(zhì)的全局角度定量對比不同儲層孔、洞、縫搭配關(guān)系的細節(jié)差異，具有較強實用性。

圖2 儲層總孔隙度中孔、洞、縫儲集空間占比的三元分析圖版（龍王廟組氣藏）

結(jié)合數(shù)字巖心流動模擬計算，揭示儲層縫、洞分布相對均勻情況下不同孔、洞、縫搭配對儲層滲透率的影響規(guī)律。龍王廟組氣藏平均孔隙度約為5%，氣藏中高滲透區(qū)域普遍發(fā)育毫米級溶洞和微裂縫，從數(shù)字巖心流動模擬計算結(jié)果看（圖3），當溶洞發(fā)育時，即使裂縫發(fā)育程度稍低也能形成滲透率較高的儲層，而溶洞欠發(fā)育時，儲層滲透率對裂縫發(fā)育程度依賴性增強。

圖3 儲層滲透率與總孔隙度中孔、洞、縫儲集空間占比的關(guān)系圖（龍王廟組氣藏）

碳酸鹽巖儲層儲滲特征不僅僅與總孔隙度中孔、洞、縫儲集空間的占比相關(guān)，還與孔隙連通率及喉道大小、洞和裂縫的尺度及分布密度等因素相關(guān)，需要多方面的影響因素分析相結(jié)合才能得到完整認識。

1.3 微小裂縫分布的統(tǒng)計規(guī)律定量描述

碳酸鹽巖儲層常見不同尺度的裂縫。分形理論是整體描述多尺度裂縫分布統(tǒng)計規(guī)律的一種有效方法[29]。該方法需要統(tǒng)計不同尺度范圍內(nèi)的裂縫數(shù)量，受深層巖石取樣尺寸較小的制約，更適合于描述長度為分米至微米級別的裂縫分布密度的統(tǒng)計特征。已有研究成果顯示，儲層中裂縫尺度越小其分布密度越大的情況較為常見，在分形描述雙對數(shù)圖上，表現(xiàn)為延伸長度大于觀測盒子范圍的裂縫數(shù)量與觀測盒子大小負相關(guān)，大體包括全觀測范圍裂縫分形維數(shù)穩(wěn)定（圖4-a）和局部范圍裂縫分形維數(shù)穩(wěn)定（圖4-b）的情況。當然，也有分形維數(shù)不穩(wěn)定、裂縫分布規(guī)律隨機性強的情況（圖4-c）。

2 儲層非均質(zhì)性對氣藏水侵動態(tài)的影響

2.1 沿宏觀大裂縫水侵規(guī)律

在氣藏開發(fā)實踐中，人們早已發(fā)現(xiàn)邊水或底水沿斷層突進至生產(chǎn)井的活躍水侵現(xiàn)象，極端情況下延伸數(shù)十米甚至更遠的大裂縫中地層水的流動類似于管道流，易導致氣井投產(chǎn)后快速地大量產(chǎn)水甚至水淹停產(chǎn)。這種大裂縫水竄機理描述被廣泛應用，并形成了相應的分析方法和針對性的治水對策[30-31]。

圖4 巖石中多尺度裂縫分布密度的分形特征典型示意圖

邊、底水氣藏開發(fā)過程中水侵極活躍現(xiàn)象是否僅與大裂縫相關(guān)？這是曾經(jīng)長期困擾人們的問題。后續(xù)研究陸續(xù)發(fā)現(xiàn)反例，即氣藏開發(fā)過程中水侵極活躍，按裂縫水竄模型定量分析推斷地層中存在延伸距離較遠的大裂縫，而地震解釋、試井探測及綜合地質(zhì)研究未發(fā)現(xiàn)相應尺度大裂縫的跡象，甚至可以完全否定其可能性。由此，人們開始關(guān)注更多類型的裂縫水竄機理。

2.2 沿小裂縫網(wǎng)絡(luò)形成的裂縫發(fā)育帶水侵規(guī)律

從總體情況看，網(wǎng)狀小尺度裂縫發(fā)育帶連通生產(chǎn)井和水區(qū)的情況比單一大裂縫水侵通道更普遍，氣藏開發(fā)研究應充分關(guān)注這類水侵規(guī)律。

考慮裂縫—孔隙型氣藏邊水沿裂縫發(fā)育帶侵入，在網(wǎng)狀小裂縫符合滲流力學連續(xù)介質(zhì)場假設(shè)的條件下，建立并求解相關(guān)滲流模型，計算分析不同條件下地層壓力、水侵量變化規(guī)律[32]，獲得定量化認識。當水侵通道儲層基質(zhì)物性差、竄流系數(shù)（λ）低至10－8或更小時，在一定時間段孔隙系統(tǒng)與裂縫系統(tǒng)之間壓差較為顯著：竄流系數(shù)越小，孔隙和裂縫介質(zhì)之間的壓差越大，并且較大壓差的持續(xù)時間越長（圖5），相應容易導致地層水沿壓力較低的網(wǎng)狀小裂縫向生產(chǎn)井竄進的發(fā)生。

圖5 不同竄流系數(shù)條件下孔隙系統(tǒng)、裂縫系統(tǒng)之間壓差隨時間變化圖

有關(guān)網(wǎng)狀小裂縫易發(fā)生水竄的條件——竄流系數(shù)臨界值的認識，對近年來活躍水侵氣藏數(shù)值模擬應用研究的發(fā)展起到了良好促進作用。以精細氣藏描述和巨量網(wǎng)格多重介質(zhì)模型大型并行計算數(shù)值模擬為基礎(chǔ)，不同于大裂縫水竄動態(tài)預測的新技術(shù)逐漸發(fā)展成熟，生產(chǎn)應用也取得突破性進展。

2.3 微裂縫對水侵規(guī)律的影響

深層碳酸鹽巖氣藏常見儲層微裂縫發(fā)育的情況，超壓氣藏這一特征更顯著，推測應與異常高壓的成因，以及地應力和超壓流體的長期耦合作用相關(guān)。

微裂縫改善儲層局部滲透率，但不會直接改變儲層宏觀非均質(zhì)性，這種情況下氣藏的水侵規(guī)律與大裂縫水竄或網(wǎng)狀小裂縫發(fā)育帶水侵明顯不同。微裂縫對水侵規(guī)律的影響情況視微裂縫分布、孔洞儲集空間特征不同而有所差異。在受水侵影響之前，對于微裂縫發(fā)育而大尺度縫洞不發(fā)育的氣藏而言，其開采動態(tài)往往與孔隙型氣藏的特征類似；水侵影響顯著后，微裂縫發(fā)育的儲層水相相對滲透率升高的變化特征比孔隙型儲層更顯著。

對微裂縫發(fā)育的儲層進行滲流特殊性的定量描述較為困難。逾滲理論是描述多孔介質(zhì)儲滲特征的一種獨特方法，尤其適用于孔隙度低、粒度細、微裂縫密度較高的儲層。已有相關(guān)研究建立了裂縫—孔隙型雙重介質(zhì)逾滲理論分析圖版[33]，揭示了不同裂縫分形維數(shù)情況下逾滲概率與孔隙度的關(guān)系。從圖6可看出，當裂縫不發(fā)育或不同尺度裂縫搭配不佳（裂縫分形維數(shù)較低）時，孔隙度超過20%才具有高度連通的特征；反之，裂縫發(fā)育能使低孔隙度條件下形成較大逾滲概率。依據(jù)理論圖版分析形成重要認識：如果特低孔隙度儲層具有高滲透能力，在排除局部大尺度縫、洞產(chǎn)生宏觀非均質(zhì)影響的情況下，可以推斷必然是微裂縫發(fā)育導致的滲流特征。

圖6 不同裂縫分形維數(shù)情況下孔隙度與逾滲概率關(guān)系圖

在缺乏鉆井取心的情況下，針對試井曲線呈現(xiàn)視均質(zhì)特征的低孔隙度儲層，通過對比測井解釋孔隙度和試井分析滲透率的關(guān)系，利用上述研究認識可以大致評估微裂縫對儲層滲流能力的貢獻，以及水侵后的影響程度。結(jié)合氣水相對滲透率實驗分析、試井計算、數(shù)值模擬，能夠進一步預測水侵后氣相滲透率下降對氣井產(chǎn)量和壓力的影響。

2.4 溶洞對水侵規(guī)律的影響

除地層中存在巨型洞穴或溶洞密度極高的特殊情形外，通常情況下各溶洞相對孤立，需要通過孔隙系統(tǒng)連通或裂縫溝通才能發(fā)揮增強宏觀滲透性的作用，這時儲層滲流特征往往仍受裂縫影響，但因孔隙系統(tǒng)物性和裂縫分布特征的不同而表現(xiàn)出差異。因此，前面關(guān)于裂縫影響水侵規(guī)律的描述在一定程度上仍然適用。

溶洞使儲層容納水侵量的能力增強，這是溶洞影響水侵規(guī)律的主要特殊性。在微觀滲流方面，表現(xiàn)為含水飽和度上升初期氣相相對滲透率有所下降，但水相相對滲透率增加不明顯，但當含水飽和度超過一定界限后，水相相對滲透率急劇增加（圖7），儲層滲透率較低時這一特征更明顯。在氣藏開發(fā)的宏觀顯現(xiàn)方面，表現(xiàn)為早期水侵活躍跡象在一定程度上表現(xiàn)滯后或被掩蓋，而當水侵危害明顯時其趨勢已經(jīng)較難逆轉(zhuǎn)。

圖7 裂縫—孔洞型巖心高溫高壓氣水相對滲透率曲線圖

3 氣藏水侵差異化特征早期預判

3.1 不同開發(fā)階段地層水體能量釋放強度差異

經(jīng)過長期研究，已經(jīng)確認水區(qū)儲集空間應力敏感是深層邊、底水氣藏強水侵能量的主要來源[3][8]。近年高溫高壓流固耦合巖心實驗分析技術(shù)取得進步，模擬地層水富集區(qū)上覆壓力和流體高壓條件，測定含水巖石有效壓縮系數(shù)的準確性大幅度提高。

圖8 不同儲滲類型巖樣有效壓縮系數(shù)實驗測定結(jié)果圖

圖8是同一碳酸鹽巖氣藏具有不同類型儲集特征的全直徑巖心實驗分析結(jié)果。實驗中圍壓固定，流體壓力從76 MPa逐級下降開展多點測試。實驗結(jié)果顯示流體壓力下降/有效應力（圍壓與流體壓力之差）上升初期，巖石有效壓縮系數(shù)從最高值開始下降，變化幅度較大，后期變化趨勢減緩并逐漸趨于穩(wěn)定；孔隙度較高的裂縫—孔洞型巖心上述特征最顯著，孔隙度較低的孔隙型巖心變化幅度較小，裂縫—孔隙型巖心變化幅度居中。裂縫—孔洞型巖心有效壓縮系數(shù)的變化幅度不僅僅與孔隙度相關(guān)，也受縫洞搭配關(guān)系（不同尺度縫洞分布密度）的影響。實驗所用裂縫—孔洞型、裂縫—孔隙型巖心僅含小裂縫，未獲得大裂縫發(fā)育情況下的結(jié)果，從原理方面分析大裂縫儲集空間應力敏感特征應更顯著。

實驗結(jié)果預示在邊、底水氣藏開發(fā)初期水區(qū)的單位壓降水侵強度最高，可能數(shù)倍于中后期水平，對于縫洞發(fā)育的儲層以及開發(fā)過程中有效應力變化較大的超壓氣藏該特征更突出。由此，結(jié)合實際氣藏具體情況，可在開發(fā)早期預判地層水體能量釋放高峰與水侵影響的對應關(guān)系，增強治水對策的針對性。

3.2 不同類型儲層水侵影響程度差異

在水體儲量、氣水接觸關(guān)系、水侵通道特征、氣藏開采方式等影響因素相同的情況下，僅比較儲滲空間類型對水侵動態(tài)的影響，發(fā)現(xiàn)以下規(guī)律：不同類型儲層產(chǎn)生非均勻水侵的可能性由大到小及水侵影響速度由快到慢的排序均為裂縫型儲層、裂縫—孔隙型儲層、裂縫—孔洞型儲層和孔隙型儲層，而水淹后排水消耗水侵能量及逆轉(zhuǎn)水侵影響的難度由易到難的排序也呈同樣的規(guī)律。表2展示了對不同類型儲層水侵影響差異化特征的機理分析結(jié)果，可作為預判水侵影響程度的類比性參照。

3.3 實際情況與理想化狀態(tài)的差異

根據(jù)已有認識劃分氣藏水侵類型，通過類比，歸類分析不同氣藏水侵特征的差異，是一種常見的研究模式。然而，僅僅依據(jù)對實際氣藏地質(zhì)特征和生產(chǎn)數(shù)據(jù)的統(tǒng)計分析，往往不足以形成對本質(zhì)特征的精細認識；對比氣藏開發(fā)實際情況與理論研究理想化狀態(tài)之間的差距，有利于獲得更準確的量化判斷。

表2 不同類型儲層水侵影響差異化特征分析結(jié)果表

根據(jù)滲流力學理論分析，當氣水接觸面水相流度（水相有效滲透率與黏度比值）不大于氣相流度、開發(fā)過程中氣藏各部位壓力均衡下降時，不會引起危害性水侵影響。上述理想化狀態(tài)在實際氣藏中很少出現(xiàn)，難以據(jù)此直接預判氣藏水侵特征，但通過定量對比分析實際情況與理想狀態(tài)的差距，能夠甄別不同氣藏、不同井區(qū)水侵影響的細分特征類型，間接支撐個性化治水對策的制定。例如，根據(jù)圖7的數(shù)據(jù)可繪制相應的氣、水流度與含水飽和度的關(guān)系曲線圖，氣、水流度曲線交點對應的含水飽和度為水相流度超過氣相流度的臨界值；從理論分析角度看，含水飽和度在臨界值以下時不易發(fā)生非均勻水侵，由此得知氣藏開發(fā)預防水侵危害的理想狀態(tài)——通過優(yōu)化調(diào)控始終使生產(chǎn)井周圍儲層含水飽和度保持在臨界值以下；雖然實際生產(chǎn)中幾乎不可能達到理想狀態(tài)，但分析實際含水飽和度與理論臨界值的差異以及對應的氣、水流度比值關(guān)系，可從一個側(cè)面評估氣藏開發(fā)受水侵影響的程度，為氣藏治水對策的制定和優(yōu)化提供參考依據(jù)。

4 研究成果應用情況

相關(guān)成果首先在四川盆地龍崗地區(qū)埋藏深度約6 000 m的三疊系飛仙關(guān)組鮞灘氣藏和二疊系長興組生物礁氣藏應用，其可靠性得到生產(chǎn)實踐的驗證，在水體能量評估、水侵影響預測、治水措施有效性預判等方面取得很好的應用效果[8]。隨后，全面推廣應用到在四川盆地新近發(fā)現(xiàn)的有水氣藏開發(fā)工作中。同時，應用新技術(shù)重新審視過去的經(jīng)驗性認識，梳理普適性技術(shù)理念和對策，在重點老氣田中后期開發(fā)的優(yōu)化治水工作中發(fā)揮較好作用[3]。

5 結(jié)論

針對三重介質(zhì)儲滲機理研究形成的方法和認識，豐富了碳酸鹽巖氣藏描述和水侵動態(tài)預測的技術(shù)手段，強化了對氣藏開發(fā)治水決策的技術(shù)支撐。

1）總孔隙度中孔洞縫占比三元圖分析、裂縫—孔隙型雙重介質(zhì)考慮裂縫分形維數(shù)的逾滲特征分析、高溫高壓條件下含水巖石流固耦合實驗測定儲集空間壓縮系數(shù)等方法已成熟，可用于研究氣藏儲滲特征及水侵規(guī)律。

2）微裂縫發(fā)育是特低孔視均質(zhì)儲層具備高滲透能力的必要條件、裂縫—孔隙型儲層孔隙系統(tǒng)與裂縫系統(tǒng)之間竄流系數(shù)小于10－8易產(chǎn)生網(wǎng)狀小裂縫水竄、洞發(fā)育情況下在水侵初期水相相對滲透率上升滯后及含水飽和度超過臨界值后水相相對滲透率急劇增加、異常高壓氣藏開采初期孔隙度應力敏感導致單位壓降水侵強度更高等規(guī)律性認識，對早期預判氣藏開發(fā)水侵動態(tài)及后續(xù)影響程度有較大幫助。

[ 1 ] 孫志道. 裂縫性有水氣藏開采特征和開發(fā)方式優(yōu)選[J]. 石油勘探與開發(fā), 2002, 29(4): 69-71.Sun Zhidao. Production characteristics and the optimization of development schemes of fractured gas reservoir with edge or bottom water[J]. Petroleum Exploration and Development, 2002,29(4): 69-71.

[ 2 ] 夏崇雙. 不同類型有水氣藏提高采收率的途徑和方法[J]. 天然氣工業(yè), 2002, 22(增刊1): 73-77.Xia Chongshuang. Ways and methods of enhancing recovery in various water-carrying gas reservoirs[J]. Natural Gas Industry,2002, 22(S1): 73-77.

[ 3 ] 馮曦, 鐘兵, 楊學鋒, 鄧惠. 有效治理氣藏開發(fā)過程中水侵影響的問題及認識[J]. 天然氣工業(yè), 2015, 35(2): 35-40.Feng Xi, Zhong Вing, Yang Xuefeng & Deng Hui. Effective water inf l ux control in gas reservoir development: Problems and countermeasure[J]. Natural Gas Industry, 2015, 35(2): 35-40.

[ 4 ] 閆海軍, 賈愛林, 郭建林, 程立華, 何東博. 龍崗礁灘型碳酸鹽巖氣藏氣水控制因素及分布模式[J]. 天然氣工業(yè), 2012,32(1): 67-70.Yan Haijun, Jia Ailin, Guo Jianlin, Cheng Lihua & He Dongbo.Gas-water controlling factors and distribution models of reef fl at carbonate rock gas reservoirs in the Longgang gas fi eld, Sichuan Вasin[J]. Natural Gas Industry, 2012, 32(1): 67-70.

[ 5 ] 劉勇, 楊洪志, 馬智超, 張春, 王容, 阮明龍, 等. 四川盆地東部高陡構(gòu)造帶巖性封隔不對稱氣水分布模式[J]. 天然氣勘探與開發(fā), 2017, 40(3): 1-7.Liu Yong, Yang Hongzhi, Ma Zhichao, Zhang Chun, Wang Rong,Ruan Minglong, et al. The lithologically sealed gas-water asymmetric distribution model of high-steep structural belt in east Sichuan Вasin[J]. Natural Gas Exploration and Development, 2017,40(3): 1-7.

[ 6 ] 陶詩平, 馮曦, 肖世洪. 應用不穩(wěn)定試井分析方法識別氣藏早期水侵[J]. 天然氣工業(yè), 2003, 23(4): 68-70.Tao Shiping, Feng Xi & Xiao Shihong. Using transient testing to identify early water invasion of gas reservoirs[J]. Natural Gas Industry, 2003, 23(4): 68-70.

[ 7 ] 何曉東, 鄒紹林, 盧曉敏. 邊水氣藏水侵特征識別及機理初探[J]. 天然氣工業(yè), 2006, 26(3): 87-89.He Xiaodong, Zou Shaolin & Lu Xiaomin. A preliminary discussion on mechanism and recognition of water invasion characteristics in edge water drive gas reservoirs[J]. Natural Gas Industry,2006, 26(3): 87-89.

[ 8 ] 馮曦, 楊學鋒, 鄧惠, 陳林, 朱占美. 根據(jù)水區(qū)壓力變化特征辨識高含硫邊水氣藏水侵規(guī)律[J]. 天然氣工業(yè), 2013, 33(1):75-78.Feng Xi, Yang Xuefeng, Deng Hui, Chen Lin & Zhu Zhanmei. Identif i cation of the water invasion law in high-sulfur and edge-water gas reservoirs based on the characteristics of pressure variation in the water zone[J]. Natural Gas Industry, 2013, 33(1):75-78.

[ 9 ] 康博, 張烈輝, 王健, 范坤, 王宏. 裂縫—孔洞型碳酸鹽巖凝析氣井出水特征及預測[J]. 西南石油大學學報(自然科學版),2017, 39(1): 107-113.Kang Вo, Zhang Liehui, Wang Jian, Fan Kun & Wang Hong. Features and forecast of water output in fractured vuggy carbonate condensate reservoir[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2017, 39(1): 107-113.

[10] Al-Ghanim JA, Nashawi IS & Malallah A. Prediction of water Infl ux of edge-water drive reservoirs using nonparametric optimal transformations[C]//North Africa Technical Conference and Exhibition, 20-22 February 2012, Cairo, Egypt. DOI: https://dx.doi.org/10.2118/150662-MS.

[11] Kabir CS, Parekh В & Mustafa MA. Material-balance analysis of gas reservoirs with diverse drive mechanisms[C]//SPE Annual Technical Conference and Exhibition, 28-30 September 2015,Houston, Texas, USA. DOI: https://dx.doi.org/10.2118/175005-MS.

[12] 李熙喆, 郭振華, 萬玉金, 劉曉華, 張滿郎, 謝武仁, 等. 安岳氣田龍王廟組氣藏地質(zhì)特征與開發(fā)技術(shù)政策[J]. 石油勘探與開發(fā), 2017, 44(3): 398-406.Li Xizhe, Guo Zhenhua, Wan Yujin, Liu Xiaohua, Zhang Manlang, Xie Wuren, et al. Geological characteristics and development strategies for Cambrian Longwangmiao Formation gas reservoir in Anyue gas fi eld, Sichuan Вasin, SW China[J]. Petroleum Exploration and Development, 2017, 44(3): 398-406.

[13] 周曉林. 帶油環(huán)底水凝析氣藏開發(fā)方式優(yōu)化及均衡開采對策——以東海盆地西湖坳陷CX北區(qū)塊為例[J]. 天然氣工業(yè),2017, 37(6): 55-63.Zhou Xiaolin. Development mode optimization and balanced production strategies for the gas condensate reservoirs with bottom water and oil rim: A case study of CX North Вlock in the Xihu Depression, East China Sea Вasin[J]. Natural Gas Industry, 2017,37(6): 55-63.

[14] Glumov DN, Sokolov SV & Strekalov AV. Assessment of drained gas reserves in the process of gas and gas condensate field operation in water drive[C]//SPE Russian Petroleum Technology Conference, 16-18 October 2017, Moscow, Russia. DOI: https://dx.doi.org/10.2118/187863-MS.

[15] 周克明, 李寧, 張清秀, 唐顯貴. 氣水兩相滲流及封閉氣的形成機理實驗研究[J]. 天然氣工業(yè), 2002, 22(增刊1): 122-125.Zhou Keming, Li Ning, Zhang Qingxiu & Tang Xiangui. Experimental research on gas-water two phase fl ow and conf i ned gas formation mechanism[J]. Natural Gas Industry, 2002, 22(S1):122-125.

[16] 疏壯志, 杜志敏, 劉建儀, 李穎川, 郭肖, 陳小凡. 碳酸鹽巖裂縫性水驅(qū)氣藏水鎖實驗研究[J]. 天然氣工業(yè), 2004, 24(6):89-92.Shu Zhuangzhi, Du Zhimin, Liu Jianyi, Li Yingchuan, Guo Xiao& Chen Xiaofan. Water-lock experiments of carbonate fractured gas reservoirs with water drive[J]. Natural Gas Industry, 2004,24(6): 89-92.

[17] Firoozabadi A, Olsen G & Van Golf-Racht T. Residual gas saturation in water-drive gas reservoirs[C]//SPE California Regional Meeting, 8-10 April 1987, Ventura, California, USA. DOI:https://dx.doi.org/10.2118/16355-MS.

[18] Вabadimas JM. Modelling trapped gas expansion in water-drive reservoirs[C]//SPE/IATMI Asia Pacific Oil & Gas Conference and Exhibition, 17-19 October 2017, Jakarta, Indonesia. DOI:https://dx.doi.org/10.2118/186272-MS.

[19] 郭平, 黃偉崗, 姜貽偉, 畢建霞, 陳召佑. 致密氣藏束縛與可動水研究[J]. 天然氣工業(yè), 2006, 26(10): 99-101.Guo Ping, Huang Weigang, Jiang Yiwei, Вi Jianxia & Chen Zhaoyou. Research on the irreducible and movable water of tight sandstone gas reservoir[J]. Natural Gas Industry, 2006, 26(10):99-101.

[20] 方飛飛, 李熙喆, 高樹生, 薛蕙, 朱文卿, 劉華勛, 等. 邊、底水氣藏水侵規(guī)律可視化實驗研究[J]. 天然氣地球科學, 2016,27(12): 2246-2252.Fang Feifei, Li Xizhe, Gao Shusheng, Xue Hui, Zhu Wenqing,Liu Huaxun, et al. Visual simulation experimental study on water invasion rules of gas reservoir with edge and bottom water[J].Natural Gas Geoscience, 2016, 27(12): 2246-2252.

[21] 劉華勛, 任東, 高樹生, 胡志明, 葉禮友, 劉夏. 邊、底水氣藏水侵機理與開發(fā)對策[J]. 天然氣工業(yè), 2015, 35(2): 47-53.Liu Huaxun, Ren Dong, Gao Shusheng, Hu Zhiming, Ye Liyou& Liu Xia. Water inf l ux mechanism and development strategy of gas reservoirs with edge and bottom water[J]. Natural Gas Industry, 2015, 35(2): 47-53.

[22] 楊朝強, 彭小東, 汪新光, 羅佼, 童璐一. 對水驅(qū)氣藏生產(chǎn)指示曲線的重新認識——以崖城13-1氣田北塊氣藏為[J]. 中國海上油氣, 2017, 29(1): 77-82.Yang Chaoqiang, Peng Xiaodong, Wang Xinguang, Luo Jiao &Tong Luyi. New understanding on production index curve of water drive gas reservoir: A case study of north block in YC13-1 gas fi eld[J]. China Oあshore Oil and Gas, 2017, 29(1): 77-82.

[23] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局, 中國國家標準化管理委員會. 天然氣藏分類: GВ/T 26979—2011[S]. 北京:中國標準出版社, 2011.General Administration of Quality Supervision, Inspection and Quarantine of the PRC & standardization Administration of the PRC. The classif i cation of natural gas pool: GВ/T 26979-2011[S].Вeijing: China Standards Press, 2011.

[24] Warren JE & Root PJ. The behavior of naturally fractured reservoirs[J]. SPE Journal, 1963, 3(3): 245-255. DOI: https://dx.doi.org/10.2118/426-PA.

[25] Wu Yushu, Ehlig-Economides CA, Qin Guan, Kang Zhuijiang,Zhang Wangming, Ajayi ВT, et al. A Triple-continuum pressure-transient model for a naturally fractured vuggy reservoir[C]//SPE Annual Technical Conference and Exhibition, 11-14 November 2007, Anaheim, California, USA. DOI: https://dx.doi.org/10.2118/110044-MS.

[26] Вerg CF, Lopez O & Вerland H. Industrial applications of digital rock technology[J]. Journal of Petroleum Science and Engineering, 2017, 157: 131-147. DOI: https://dx.doi.org/10.1016/j.petrol.2017.06.074.

[27] Zhao Tianyi, Zhao Huawei, Ning Zhengfu, Li Xiangfang & Wang Qing. Permeability prediction of numerical reconstructed multiscale tight porous media using the representative elementary volume scale lattice Вoltzmann method[J]. International Journal of Heat and Mass Transfer, 2018, 118: 368-377. DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2017.11.004.

[28] 鄢友軍, 李隆新, 徐偉, 常程, 鄧惠, 楊柳. 三維數(shù)字巖心流動模擬技術(shù)在四川盆地縫洞型儲層滲流研究中的應用[J]. 天然氣地球科學, 2017, 28(9): 1425-1432.Yan Youjun, Li Longxin, Xu Wei, Chang Cheng, Deng Hui &Yang Liu. Application of 3D digital core flow simulation technique in the study of gas fl ow in fractured-vuggy gas reservoirs in Sichuan Вasin[J]. Natural Gas Geoscience, 2017, 28(9): 1425-1432.

[29] Вabadagli T. Scanline method to determine the fractal nature of 2-D fracture networks[J]. Mathematical Geology, 2002, 34(6):647-670. DOI: https://dx.doi.org/10.1023/A:1019897008208.

[30] 陳立官, 李鴻智, 劉文碧, 張長盛, 趙玉隆, 楊茂康, 等. 試論在川南二疊系陽新統(tǒng)中找氣的新途徑——排水找氣[J]. 天然氣工業(yè), 1986, 6(3): 35-41.Chen Liguan, Li Hongzhi, Liu Wenbi, Zhang Changsheng, Zhao Yulong, Yang Maokang, et al. Trying to demonstrate the new way for fi nding gas in Yangxi Series of Permian in Southern Sichuan—f i nding gas by water withdrawal[J]. Natural Gas Industry,1986, 6(3): 35-41.

[31] 羅濤, 王陽. 裂縫水竄氣藏單井數(shù)值模擬研究[J]. 天然氣工業(yè),2002, 22(增刊1): 95-97.Luo Tao & Wang Yang. Research on individual well numerical simulation of fracture water breakthrough gas reservoir[J]. Natural Gas Industry, 2002, 22(S1): 95-97.

[32] 馮曦, 鐘兵, 潘寧, 王浩, 李海濤. 裂縫—孔隙型氣藏邊水沿裂縫發(fā)育帶侵入動態(tài)特征[J]. 天然氣工業(yè), 2008, 28(增刊1):79-81.Feng Xi, Zhong Вing, Pan Ning, Wang Hao & Li Haitao. Invasion dynamic characteristics of edge water along fractured zone in fractured pore gas reservoir[J]. Natural Gas Industry, 2008,28(S1): 79-81.

[33] 趙陽升. 多孔介質(zhì)多場耦合作用及其工程響應[M]. 北京: 科學出版社, 2010.Zhao Yangsheng. Multi field coupling of porous media and its engineering response[M]. Вeijing: Science Press, 2010.