張麗華，潘保芝，莊華，郭立新，李慶峰，趙小青

（1.吉林大學地球探測科學與技術(shù)學院，吉林長春130026；2.大慶鉆探工程公司測井公司，黑龍江大慶163412）

低孔隙度低滲透率儲層壓裂后產(chǎn)能測井預測方法研究

張麗華1，潘保芝1，莊華1，郭立新2，李慶峰2，趙小青2

（1.吉林大學地球探測科學與技術(shù)學院，吉林長春130026；2.大慶鉆探工程公司測井公司，黑龍江大慶163412）

通過研究影響壓裂后產(chǎn)能的各種參數(shù)，對低孔隙度低滲透率的朝陽溝油田進行了壓裂后產(chǎn)能預測。對壓裂后產(chǎn)能影響較大的參數(shù)有地層有效滲透率、有效厚度、井底流壓等。確定了該區(qū)產(chǎn)能級別劃分標準，提出2種預測方法：方程法和神經(jīng)網(wǎng)絡法。建立了朝陽溝油田預測壓裂后產(chǎn)能的無限導流能力垂直裂縫穩(wěn)態(tài)流方程，并給出了各參數(shù)的計算方法。應用方程法對11口井的壓裂后產(chǎn)能進行預測，判別單層測試產(chǎn)油級別符合率為73%。應用神經(jīng)網(wǎng)絡法，選取研究區(qū)13口井進行訓練和驗證，單層測試產(chǎn)油級別符合率為100%。

測井評價；產(chǎn)能預測；壓裂；低孔隙度；低滲透率；產(chǎn)能級別；方程法；神經(jīng)網(wǎng)絡法

0 引 言

油氣田勘探與開發(fā)的一項基本任務是對儲層產(chǎn)能進行定性或定量評價。Russell和Truitt［1］對壓裂井滲流方式的分析比較清晰，指出在無限導流垂直裂縫附近地層將出現(xiàn)線性流，而遠方地層出現(xiàn)徑向流，并給出了校正用普通直井壓力恢復半對數(shù)直線段解釋的kh值方法。鞠江慧［2］建立了由次生孔隙度、孔隙結(jié)構(gòu)系數(shù)、總孔隙度等參數(shù)回歸的壓裂后產(chǎn)能預測方程。張松揚［3］將測井特征參數(shù)比值加權(quán)合成法用于產(chǎn)能預測。此外，袁云福［4］、何巖峰［5］等許多研究者［6－17］都做了產(chǎn)能預測方面的研究工作。但多數(shù)文獻都是對自然產(chǎn)能進行預測，對壓裂后產(chǎn)能預測的研究較少。利用常規(guī)測井資料主要以建立經(jīng)驗關系的方式實現(xiàn)產(chǎn)能預測，其精度、適用性都不強。隨著算法的不斷改進，產(chǎn)能預測的精度有所提高。但到目前為止，還沒有一種綜合方法，既從產(chǎn)能公式的各項參數(shù)入手，又綜合利用現(xiàn)代測井的各種資料進行儲層的壓裂后產(chǎn)能預測。大慶外圍朝陽溝油田扶楊油層孔隙度范圍為14.77%～17.77%，屬于中孔隙度油層。由于原生孔隙小，沉積物中鈣積現(xiàn)象比較普遍，成巖再生作用比較強，致使孔隙度變低，也使?jié)B透率變得很低，空氣滲透率為（2.6～22.2）×10－3μm2。儲層自然產(chǎn)能較低，要對儲層進行改造才能見到較好的效果。

本文從影響壓裂后產(chǎn)能的各種參數(shù)入手，建立預測朝陽溝油田壓裂后產(chǎn)能的方程和方法，對產(chǎn)能進行預測。

1 壓裂后產(chǎn)能影響因素及產(chǎn)能級別劃分

低孔隙度低滲透率儲層一般沒有自然產(chǎn)能，需要經(jīng)過壓裂改造才能投產(chǎn)。現(xiàn)代水力壓裂過程是通過井筒向目標地層泵注壓裂液以高壓形成裂縫，而后注入混有支撐劑的攜砂液繼續(xù)延展裂縫，同時攜帶支撐劑深入裂縫以備支撐，再后注入破膠劑使壓裂液降解為低黏度的流體，最后井口放噴卸載，注入液體流向井底反排而出。這樣，就能夠在地層中留下1條水力裂縫——高導流能力的油氣通道。未壓裂的水力完善井的滲流方式是平面徑向流，對于低孔隙度低滲透率儲層，進行水力壓裂后進行開采，如果壓裂后能夠產(chǎn)生有一定傳導能力的垂直裂縫，那么原來的平面徑向滲流模式將改變?yōu)槠矫婢€性流（高導流能力裂縫）或雙線性滲流模式（低導流能力裂縫），其線性流或雙線性流的流場特點是流線和等勢線平行分布，和徑向流模式相比，線性流或雙線性流模式的滲流阻力要小得多，是油氣開采中的有利滲流模式，故經(jīng)壓裂后可產(chǎn)生增產(chǎn)效果（見圖1）。

圖1 壓裂前后滲流方向變化圖

影響產(chǎn)能的主要因素包含幾何因素、油藏類型、鉆井液與完井、強化增產(chǎn)等4個方面。其中幾何因素包含儲層厚度、滲透率各向異性、井眼尺寸等；油藏類型包含井在油藏中的位置、井網(wǎng)等，鉆井液與完井包含表皮效應對產(chǎn)能的影響、完井方式等；強化增產(chǎn)包含壓裂、酸化等。

從油層壓裂機理的研究得知，油層壓裂后的產(chǎn)能增長與所形成的裂縫填砂面積大小有直接關系。填砂面積大，則效果好。而加砂量又取決于排量、壓裂液用量與壓裂液類型等一系列工藝參數(shù)的選擇。只有在大排量下或優(yōu)質(zhì)壓裂液（高黏度，低濾失量）和足夠體積的壓裂液用量下才能獲得大的裂縫填砂面積和較好的壓效。

表1 產(chǎn)能級別劃分標準表

2 壓裂后產(chǎn)能預測方程

2.1 建立壓裂后產(chǎn)能預測方程

無限導流能力垂直裂縫井穩(wěn)態(tài)滲流產(chǎn)能公式［18］為

式中，q為儲層的產(chǎn)量，t／d；K為儲層內(nèi)流體的有效滲透率，mD＊＊非法定計量單位，1mD＝9.87×10－4μm2，下同；h為儲層有效厚度，m；pavg為泄流區(qū)的平均壓力，MPa；pw為流壓，MPa；μ為流體黏度，mPa·s；B為體積系數(shù)；re為有效供油半徑，m；ws、Ks分別為裂縫污染部分的寬度和滲透率；xf為水力垂直裂縫半長，m；Sf為表皮因子。

則式（1）轉(zhuǎn)化為

2.2 方程中基本參數(shù)的確定

由式（1）、式（2）可見，與油層的壓裂后產(chǎn)能有關的參數(shù)不僅為測井綜合參數(shù)，還有儲層中流體的黏度、流體的有效滲透率、流體的體積系數(shù)、供油半徑以及壓裂施工參數(shù)等。這些參數(shù)是利用無限導流產(chǎn)能公式解決油層壓裂后產(chǎn)能預測問題的關鍵。

2.2.1 地層原油黏度

根據(jù)式（1），儲層中原油產(chǎn)量與原油的黏度密切相關。對于同樣物性、埋藏深度、地溫及含油豐度條件下，原油黏度低的儲層比原油黏度高的儲層獲得的產(chǎn)能要高。原油的黏度主要受地層的溫度、構(gòu)成原油的組分和天然氣在原油中溶解度的影響。壓力低于或等于泡點壓力時

壓力大于泡點壓力時

浮游動物豐度及季節(jié)變化如圖2，四次調(diào)查汾河二庫浮游動物平均豐度是8 550個/L，變化范圍是3 441～16 811個/L，浮游動物總豐度在2017年春季最高，為16 811個/L，2016年春季及夏季，2017年夏季浮游動物總豐度較低，分為是9 951個/L，3 441個/L，3996個/L。原生動物和輪蟲的豐度在2016年夏季和2017年夏季較低，枝角類和橈足類在2017年春季最高，汾河二庫浮游動物豐度的主體為原生動物和輪蟲的。

式中，pb為飽和壓力，MPa；μ0為pwf壓力下的地層原油黏度，×10－3Pa·s；μob為pb壓力下的地層原油黏度，×10－3Pa·s。

2.2.2 飽和壓力

飽和壓力又稱泡點壓力，表示在地層條件下，原油中的溶解氣開始分離出來時的壓力。

式中，RS為溶解氣油比，m3／m3；Twf為地層溫度，℃；γg為天然氣的相對密度。

2.2.3 流壓和滲透率

通過測試資料中壓裂井的流壓資料，確定研究區(qū)的平均流壓為1.2MPa。

通過測試資料中壓裂井油層的壓裂后滲透率和地層滲透率的關系求取壓裂后滲透率Ky。確定研究區(qū)的壓裂后滲透率計算公式為

2.2.4 X的確定

通過研究區(qū)13個壓裂后單層測試井的X與混砂量和壓入地層總液量之間的關系看出，在X小于5時，X與二者之間有很好的相關性。X小于5的數(shù)據(jù)點測試的深度都大于500m。當深度大于500 m時，X與混砂量Vsand和壓入地層總液量Vfluid之間的關系如圖2、圖3所示。得到的關系式為

相關系數(shù)為0.898。

2.3 應用效果分析

應用式（13）對11口井進行計算，結(jié)果見表2。從表2中可知有3個層的產(chǎn)能級別與試油的不符合。產(chǎn)能級別符合率為73%。

3 BP神經(jīng)網(wǎng)絡法預測壓裂后產(chǎn)能

BP（Back Propagation）神經(jīng)網(wǎng)絡［19］能學習和存貯大量的輸入－輸出模式映射關系，而無需事前揭示描述這種映射關系的數(shù)學方程。它的學習規(guī)則是使用最速下降法，通過反向傳播不斷調(diào)整網(wǎng)絡的權(quán)值和閾值，使網(wǎng)絡的誤差平方和最小。利用神經(jīng)網(wǎng)絡預測產(chǎn)能時，選擇不同的測井曲線和壓裂施工參數(shù)作為輸入層，單位厚度產(chǎn)油量作為輸出，從而對比產(chǎn)能與各輸入?yún)?shù)的相關程度及預測效果。統(tǒng)計分析了朝陽溝油田13口井13個單層壓裂改造油層段的資料，根據(jù)該地區(qū)測井響應特征，優(yōu)選自然伽馬、聲波時差、電阻率、中子、密度等5條測井曲線的特征參數(shù)，以及黏度、壓力、壓入地層總液量和混砂液量4個動態(tài)參數(shù)。依據(jù)該區(qū)經(jīng)驗規(guī)律，將這9個參數(shù)提取出來作為BP神經(jīng)網(wǎng)絡模型的輸入，單位厚度產(chǎn)油量作為輸出進行訓練，然后利用該網(wǎng)絡對壓裂后儲層單位厚度產(chǎn)油量進行預測。

在研究區(qū)13口單層壓裂井中選取其中10口井作為訓練用數(shù)據(jù)，3口作為驗證井。訓練和驗證的結(jié)果見表3。經(jīng)計算，產(chǎn)能級別符合率為100%。

表2 壓裂后實際產(chǎn)能級別與計算產(chǎn)能級別對比表

表3 訓練和驗證結(jié)果表

圖4 ×井方程法和神經(jīng)網(wǎng)絡法計算產(chǎn)能成果圖

4 應用實例

圖4是×井計算壓裂產(chǎn)能的實例。1 182～1 183.8m井段方程法計算的產(chǎn)油量為1.258 5t／d，據(jù)此判斷產(chǎn)能級別為低產(chǎn)工業(yè)油層。用神經(jīng)網(wǎng)絡法預測的產(chǎn)油量為1.408t／d。該井段試油結(jié)論是低產(chǎn)工業(yè)油層，日產(chǎn)油1.419t?？梢钥闯?，計算和預測的產(chǎn)能級別與試油的級別是符合的。

5 結(jié) 論

（1）采用無限導流能力垂直裂縫穩(wěn)態(tài)流方程對壓裂后產(chǎn)能預測，建立了方程中重要參數(shù)X的計算方法，利用該方程對壓裂產(chǎn)能進行預測，判別的單層測試產(chǎn)油級別符合率為73%。雖然利用該方法預測產(chǎn)能級別的符合率達不到100%，但該方法有理論基礎，適用于各個地區(qū)，只要方程中各參數(shù)能確定，就可直接應用該公式計算產(chǎn)量。

（2）采用神經(jīng)網(wǎng)絡方法預測壓裂后產(chǎn)能，對于單層測試樣本，輸入電阻率Rt、聲波時差AC等測井特征參數(shù)和地層壓力、原油黏度、壓入地層總液量和混砂液量，判別的單層測試產(chǎn)油級別符合率為100%。利用該方法預測產(chǎn)能級別的符合率很高，但該方法適用性較差，只能適用于研究區(qū)，如果在別的地區(qū)應用，要重新訓練網(wǎng)絡才能應用該網(wǎng)絡對產(chǎn)能進行預測。

［1］ Russell D G，Truitt N E.Transient Pressure Behavior in Vertically Fractured Reservoir［C］∥SPE 967，1964.

［2］ 鞠江慧，王建功，高瑞琴，等.二連油田低孔隙度低滲透率儲層壓裂后產(chǎn)能預測［J］.測井技術(shù)，2005，29（4）：379－382.

［3］ 張松揚，范宜仁，黃國騫，等.常規(guī)測井特征比值法在大牛地氣田產(chǎn)能評價研究中的應用［J］.測井技術(shù)，2006，30（5）：420－424.

［4］ 袁云福，高楚橋.利用MDT和NMR資料確定儲層產(chǎn)能預測參數(shù)［J］.江漢石油學院學報，2004，26（增刊）：68－71.

［5］ 何巖峰，吳曉東，韓增軍，等.低滲透油藏油井產(chǎn)能預測新方法［J］.中國石油大學學報：自然科學版，2007，31（5）：69－73.

［6］ 楊雄文，樊洪海，廖新維，等.基于電纜地層測試的儲層產(chǎn)能預測方法研究［J］.海洋石油，2008，28（2）：105－109.

［7］ 湯述安，鄭澤忠，朱學波.多元回歸分析在油氣產(chǎn)能預測中的應用［J］.科技創(chuàng)新導報，2008，21：168－171.

［8］ 劉英，吳剛，燕金梅，等.M油田儲層流體性質(zhì)識別與產(chǎn)能預測［J］.測井技術(shù)，2004，28（5）：458－461.

［9］ 譚成仟，宋子齊，吳少波.利用灰色關聯(lián)分析法綜合評價油氣儲層產(chǎn)能［J］.河南石油，2001，15（2）：20－24.

［10］譚成仟，馬娜蕊，蘇超.儲層油氣產(chǎn)能的預測模型和方法［J］.地球科學與環(huán)境學報，2004，26（2）：43－48.

［11］高建華.測井儲層產(chǎn)能預測方法研究［D］.成都：成都理工大學，2002.

［12］李瑞，楊光惠，胡奇凱.鄂爾多斯盆地碳酸鹽巖儲層測井產(chǎn)能預測研究［J］.勘探地球物理進展，2003，26（2）：109－113.

［13］田冷，何順利，顧岱鴻，等.改進BP神經(jīng)網(wǎng)絡模型在長慶氣田產(chǎn)能預測中的應用［J］.石油天然氣學報，2008，30（5）：106－111.

［14］李昌彪，宋建平，夏克文.利用RBF神經(jīng)網(wǎng)絡自適應調(diào)整算法預測儲層產(chǎn)能［J］.石油地球物理勘探，2006，41（1）：53－57.

［15］蔡左花，匡建超，曾劍毅，等.新場氣田沙溪廟組氣藏單井產(chǎn)能預測［J］.鉆采工藝，2009，32（3）：34－37.

［16］張鋒，張星，張樂，等.利用支持向量機方法預測儲層產(chǎn)能［J］.西南石油大學學報，2007，29（3）：24－29.

［17］童凱軍，單鈺銘，李海鵬，等.支持向量回歸機在氣井產(chǎn)能預測中的應用［J］.新疆石油地質(zhì)，2008，29（3）：382－385.

［18］蔡長宇.水力壓裂井產(chǎn)能研究［D］.北京：中國地質(zhì)大學，2005.

［19］楊銘震，王燕霞.人工神經(jīng)網(wǎng)絡及其在石油勘探中的應用［M］.北京：兵器工業(yè)出版社，1993.

Productivity Log Forecasting Method for Post－frac Reservoir with Low Porosity and Low Permeability

ZHANG Lihua1，PAN Baozhi1，ZHUANG Hua1，GUO Lixin2，LI Qingfeng2，ZHAO Xiaoqing2
（1.College of Geo－exploration Science and Technology，Jilin University，Changchun，Jilin 130026，China；2.Well Logging Company of Daqing Petroleum Administrative Bureau，Daqing，Heilongjiang 163412，China）

Post－frac productivity forecasting in Chaoyanggou oilfield with low porosity and low permeability is made by means of researches on relative factors.The major factors which influence the post－frac productivity include：formation effective permeability，effective thickness，downhole fluid pressure，etc.Determined is the productivity level standard，and given are 2prediction methods：equation method and neural network method.Built is the stable state flow equation method of infinite flow conductivity vertical fracture，and worked out is computing method for each parameter.Equation method is used to forecast the post－frac productivity in 11wells，and productivity grade coincidence rate of single layer is up to 73%.When using neural network method to do training and prediction in 13wells，the productivity grade coincidence rate reaches 100%.

log evaluation，productivity forecast，fracturing，low porosity，low permeability，productivity grade，equation method，neural network method

P631.84 文獻標識碼：A

2011－05－25 本文編輯 王小寧）

國家重大專項課題（No.2008ZX05010－002）

陳亦寒，女，1974年生，在站博士后，主要從事儲層地質(zhì)研究。