基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法
【專(zhuān)利摘要】本發(fā)明涉及基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法，獲取工區(qū)已試井頁(yè)巖氣層參數(shù)，獲取已試井頁(yè)巖氣層壓裂參數(shù)，通過(guò)測(cè)試產(chǎn)能資料獲取已試井頁(yè)巖氣層無(wú)阻流量；由獲取的參數(shù)依次計(jì)算已試井氣層含氣指數(shù)、壓裂改造指數(shù)、產(chǎn)氣指數(shù)，與對(duì)應(yīng)已試井計(jì)算的產(chǎn)氣指數(shù)IPG利用最小二乘法確定模型AOFg＝A×eB×IPG的模型系數(shù)；獲取工區(qū)待預(yù)測(cè)井頁(yè)巖氣層參數(shù)，獲取待預(yù)測(cè)井頁(yè)巖氣層壓裂參數(shù)，獲取的參數(shù)根據(jù)公式依次計(jì)算待預(yù)測(cè)井氣層壓裂改造指數(shù)、含氣指數(shù)、產(chǎn)氣指數(shù)；將產(chǎn)氣指數(shù)帶入模型AOFg1＝A×eB×IPG1計(jì)算待預(yù)測(cè)井頁(yè)巖氣層無(wú)阻流量，輸出計(jì)算結(jié)果。本發(fā)明已在某頁(yè)巖氣田預(yù)測(cè)30口井采用，平均誤差不超過(guò)20％。
【專(zhuān)利說(shuō)明】
基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法
技術(shù)領(lǐng)域
[0001]本發(fā)明涉及一種基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法，主要為頁(yè)巖儲(chǔ)層測(cè)錄井解釋及壓裂改造提供關(guān)鍵評(píng)價(jià)參數(shù)。
【背景技術(shù)】
[0002]無(wú)阻流量預(yù)測(cè)是天然氣勘探開(kāi)發(fā)中一個(gè)十分重要的環(huán)節(jié)，是對(duì)儲(chǔ)層品質(zhì)、儲(chǔ)層流體含油氣性質(zhì)和其產(chǎn)能的綜合評(píng)價(jià)指標(biāo)，同時(shí)也是最重要的指標(biāo)之一。
[0003]頁(yè)巖氣是指主體位于暗色泥頁(yè)巖或高碳泥頁(yè)巖中，以吸附或游離狀態(tài)為主要存在方式的天然氣聚集。頁(yè)巖氣藏具有超低孔隙度、滲透率的特點(diǎn);且部分氣體吸附于儲(chǔ)層巖石顆粒表面，以吸附態(tài)存在于儲(chǔ)層中，這給頁(yè)巖氣解釋評(píng)價(jià)帶來(lái)困難。頁(yè)巖氣井投產(chǎn)前都需要進(jìn)行大型壓裂，從而導(dǎo)致頁(yè)巖氣藏產(chǎn)能與常規(guī)氣藏相比更有其特殊復(fù)雜的影響因素。
[0004]傳統(tǒng)的絕對(duì)無(wú)阻流量獲取方法主要為一點(diǎn)法和擬壓力法等，這些都是確定原始?xì)獠貕毫ο陆o(wú)阻流量方程。對(duì)于頁(yè)巖氣藏來(lái)說(shuō)，由于滲流機(jī)理復(fù)雜，通過(guò)這些方法難以建立適用的無(wú)阻流量模型。

【發(fā)明內(nèi)容】

[0005]本發(fā)明的目的是針對(duì)上述技術(shù)現(xiàn)狀，旨在提供一種方法簡(jiǎn)便、適用范圍廣，并且能為頁(yè)巖氣層提供可靠資料，提高預(yù)測(cè)率的基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法。
[0006]本發(fā)明目的的實(shí)現(xiàn)方式為，基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法，具體步驟為:
[0007]I)通過(guò)測(cè)錄井資料獲取工區(qū)已試井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H、水平井段長(zhǎng)度L、壓力系數(shù)Kf、孔隙度P0R、含氣飽和度Sg、有機(jī)碳含量TOC和脆性指數(shù)Brit;
[0008]2)獲取已試井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf、壓裂支持劑體積Vs;
[0009]3)通過(guò)測(cè)試產(chǎn)能資料獲取已試井頁(yè)巖氣層無(wú)阻流量AOFg;
[0010]4)運(yùn)用步驟I)和步驟2)獲取的參數(shù)依次計(jì)算已試井氣層含氣指數(shù)Ig、壓裂改造指數(shù)IFR、產(chǎn)氣指數(shù)IPG，具體公式如下:
[0011]IFR=(Vf+Vs)/10000,
[0012]Ig = HXKf XPORX Sg X TOC,
[0013]IPG = LX (Ig+IFR) XBrit；
[0014]式中:壓裂液體積Vf量綱為m3，
[0015]壓裂支持劑體積Vs量綱為m3，
[0016]壓裂改造指數(shù)IFR量綱為10km3，
[0017]厚度H量綱為m，
[0018]水平井段長(zhǎng)度L量綱為hm，
[0019]壓力系數(shù)Kf無(wú)量綱，
[0020]孔隙度P0R、含氣飽和度Sg、有機(jī)碳含量TOC量綱為小數(shù)，
[0021 ]含脆性指數(shù)Brit量綱為小數(shù)；
[0022]5)將獲取的已試井頁(yè)巖氣層無(wú)阻流量AOFg與對(duì)應(yīng)已試井計(jì)算的產(chǎn)氣指數(shù)IPG利用最小二乘法確定模型A0Fg=AXeBXIPG的模型系數(shù)A、B;
[0023]式中無(wú)阻流量AOFg量綱為104m3/d;
[0024]6)通過(guò)測(cè)錄井資料獲取工區(qū)待預(yù)測(cè)井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H 1、水平井段長(zhǎng)度L 1、壓力系數(shù)Kf 1、孔隙度POR 1、含氣飽和度Sg 1、有機(jī)碳含量TOC I和脆性指數(shù)Brit I;
[0025]7)獲取待預(yù)測(cè)井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf 1、壓裂支持劑體積 V s I ；
[0026]8)運(yùn)用步驟6)和步驟7)獲取的參數(shù)根據(jù)公式
[0027]IFR l = (Vf I+Vs I)/10000，
[0028]Ig I=H IXKf 1XP0R IXSg 1XT0C 1，
[0029]IPG I=L I X (Ig 1+IFR I) XBrit I；
[0030]依次計(jì)算待預(yù)測(cè)井氣層壓裂改造指數(shù)IFR1、含氣指數(shù)Ig 1、產(chǎn)氣指數(shù)IPGl;
[0031]9)將步驟5)計(jì)算出的模型系數(shù)A、B帶入AOFg I=AXe6xipc 1中計(jì)算出的待預(yù)測(cè)井頁(yè)巖氣層無(wú)阻流量AOFg I；
[0032]10)輸出計(jì)算結(jié)果。
[0033]本發(fā)明解決了傳統(tǒng)無(wú)阻流量獲取方法在預(yù)測(cè)頁(yè)巖氣層無(wú)阻流量精確度不高的問(wèn)題，并且比傳統(tǒng)絕對(duì)無(wú)阻流量獲取方法更簡(jiǎn)便，適用范圍更廣。
[0034]本發(fā)明已在某頁(yè)巖氣田預(yù)測(cè)30口井絕對(duì)無(wú)阻流量，本發(fā)明預(yù)測(cè)的頁(yè)巖氣層絕對(duì)無(wú)阻流量與現(xiàn)場(chǎng)生產(chǎn)開(kāi)發(fā)的絕對(duì)無(wú)阻流量接近，平均誤差不超過(guò)20%，符合現(xiàn)場(chǎng)絕對(duì)無(wú)阻流量預(yù)測(cè)需要。
【附圖說(shuō)明】
[0035]圖1為本發(fā)明工作流程框圖，
[0036]圖2為本發(fā)明J工區(qū)9口井絕對(duì)無(wú)阻流量與產(chǎn)氣指數(shù)交會(huì)圖。
【具體實(shí)施方式】
[0037]參照?qǐng)D1，本發(fā)明的具體步驟為:
[0038]I)通過(guò)測(cè)錄井資料獲取工區(qū)已試井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H、水平井段長(zhǎng)度L、壓力系數(shù)Kf、孔隙度P0R、含氣飽和度Sg、有機(jī)碳含量TOC以及脆性指數(shù)Brit;
[0039]2)獲取已試井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf、壓裂支持劑體積Vs;
[0040]3)獲取已試井頁(yè)巖氣層無(wú)阻流量AOFg;
[0041]3)通過(guò)測(cè)試產(chǎn)能資料獲取已試井頁(yè)巖氣層無(wú)阻流量AOFg;
[0042]4)運(yùn)用步驟I)和步驟2)獲取的參數(shù)依次計(jì)算計(jì)算已試井氣層含氣指數(shù)Ig、壓裂改造指數(shù)IFR、產(chǎn)氣指數(shù)IPG，具體公式如下:
[0043]IFR=(Vf+Vs)/10000,
[0044]Ig = HXKf XPORX Sg X TOC,
[0045]IPG = LX (Ig+IFR) XBrit；
[0046]5)將獲取的已試井頁(yè)巖氣層無(wú)阻流量AOFg與對(duì)應(yīng)已試井計(jì)算的產(chǎn)氣指數(shù)IPG利用最小二乘法確定模型A0Fg=AXeBXIPG的模型系數(shù)A、B;
[0047]6)通過(guò)測(cè)錄井資料獲取工區(qū)待預(yù)測(cè)井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H 1、水平井段長(zhǎng)度L 1、壓力系數(shù)Kf 1、孔隙度POR 1、含氣飽和度Sg 1、有機(jī)碳含量TOC I以及脆性指數(shù)Br it I ；
[0048]8)運(yùn)用步驟6)和步驟7)獲取的參數(shù)依次計(jì)算待預(yù)測(cè)井氣層壓裂改造指數(shù)IFR、含氣指數(shù)Ig、產(chǎn)氣指數(shù)IPG;
[0049]7)獲取待預(yù)測(cè)井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf 1、壓裂支持劑體積 V s I ；
[0050]8)運(yùn)用步驟6)和步驟7)獲取的參數(shù)依次根據(jù)公式[0051 ] IFR l = (Vf I+Vs I)/10000，
[0052]Ig I=H IXKf 1XP0R IXSg 1XT0C 1，
[0053]IPG I=L I X (Ig 1+IFR I) XBrit I；
[0054]依次計(jì)算待預(yù)測(cè)井氣層壓裂改造指數(shù)IFR1、含氣指數(shù)Ig 1、產(chǎn)氣指數(shù)IPG I;
[0055]9)將步驟5)計(jì)算出的模型系數(shù)A、B帶入AOFg I=AXe6xipc 1中計(jì)算出的待預(yù)測(cè)井頁(yè)巖氣層無(wú)阻流量AOFg I；
[0056]10)輸出計(jì)算結(jié)果。
[0057]下面用具體實(shí)施例詳述本發(fā)明。
[0058]實(shí)例一:某頁(yè)巖氣田J工區(qū)R2HF井
[0059]I)通過(guò)測(cè)錄井資料獲取J工區(qū)9口已試井(WlH井、W1-2H井、W1-3HF井、W6-2HF井、W7-2HF井、W8-2HF井、W10_2HF井、W11-2HF井、W12-3HF井)頁(yè)巖氣層的氣層厚度H、水平井段長(zhǎng)度L、壓力系數(shù)Kf、孔隙度P0R、含氣飽和度Sg、有機(jī)碳含量TOC以及脆性指數(shù)Brit;
[0060]2)獲取9 口已試井頁(yè)巖氣層壓裂所用的壓裂液體積Vf、壓裂支持劑體積Vs;
[0061]3)通過(guò)測(cè)試產(chǎn)能資料獲取9 口已試井頁(yè)巖氣層無(wú)阻流量AOFg;
[0062]4)根據(jù)步驟I)和步驟2)獲取的參數(shù)依次計(jì)算9 口已試井氣層壓裂改造指數(shù)IFR、含氣指數(shù)Ig、產(chǎn)氣指數(shù)IPG;
[0063]5)將9 口已試井的頁(yè)巖氣層無(wú)阻流量AOFg與步驟4)計(jì)算的產(chǎn)氣指數(shù)IPG利用最小二乘法確定模型A0Fg=AXeBXIP(%模型系數(shù)A=9.3173、B = 0.0238(見(jiàn)圖2);
[0064]6)通過(guò)測(cè)錄井資料獲取R2HF井頁(yè)巖氣層參數(shù)，R2HF井氣層厚度H為40m、水平井段長(zhǎng)度L I為14.58m、壓力系數(shù)Kf I平均為1.45、孔隙度POR I平均為0.059、含氣飽和度Sg I平均為0.671、有機(jī)碳含量TOC I平均為3.6%、脆性指數(shù)Brit I平均為0.65;
[0065]7)獲取R2HF井頁(yè)巖氣層壓裂參數(shù)，壓裂液體積Vf I為30963.2m3、壓裂支持劑體積Vs I為778.4m3;
[0066]8)運(yùn)用步驟6)和步驟7)獲取的參數(shù)計(jì)算的R2HF井氣層壓裂改造指數(shù)IFR I為3.174、含氣指數(shù)Ig I為8.266、產(chǎn)氣指數(shù)IPG I為116.825;
[0067]9)利用模型AOFg 1=9.3173父6()'()2381[5(；1計(jì)算1?2冊(cè)井頁(yè)巖氣層絕對(duì)無(wú)阻流量；
[0068]10)輸出計(jì)算結(jié)果，R2HF井計(jì)算絕對(duì)無(wú)阻流量為150.25X 104m3/d，頁(yè)巖氣開(kāi)發(fā)測(cè)試絕對(duì)無(wú)阻流量為144.89 X 104m3/d，誤差為3.7%,小于20.0% ,符合現(xiàn)場(chǎng)絕對(duì)無(wú)阻流量預(yù)測(cè)需要。
[0069]實(shí)例二:某頁(yè)巖氣田J工區(qū)R3HF井
[0070]I )R3HF井和與實(shí)例一中的R2HF井都屬于J工區(qū)的生產(chǎn)井，可以沿用實(shí)例一中9 口已試井計(jì)算確定的模型系數(shù)A、B，模型AOFg=AX eBXIP(%模型系數(shù)A=9.3173、B = 0.0238;
[0071]2)通過(guò)測(cè)錄井資料獲取R3HF井頁(yè)巖氣層參數(shù)，R3HF井氣層厚度H為40m、水平井段長(zhǎng)度L為13.98m、壓力系數(shù)Kf平均為1.45、孔隙度POR平均為0.055、含氣飽和度Sg平均為0.552、有機(jī)碳含量TOC平均為2.36%、脆性指數(shù)Brit平均為0.55 ；
[0072]3)獲取R3HF井頁(yè)巖氣層壓裂參數(shù)，壓裂液體積Vf為28627.5m3、壓裂支持劑體積Vs為1015.4m3;
[0073]4)運(yùn)用步驟2)和步驟3)獲取的參數(shù)計(jì)算的R3HF井氣層壓裂改造指數(shù)IFR為2.964、含氣指數(shù)Ig為4.156、產(chǎn)氣指數(shù)IPG為56.359;
[0074]5)利用模型40?8 = 9.3173\6—[5(；計(jì)算1?冊(cè)井頁(yè)巖氣層絕對(duì)無(wú)阻流量；
[0075]6)輸出計(jì)算結(jié)果，R3HF井計(jì)算絕對(duì)無(wú)阻流量為35.63X 104m3/d，頁(yè)巖氣開(kāi)發(fā)測(cè)試絕對(duì)無(wú)阻流量為40.40X104m3/d，誤差為11.8%，小于20.0%，符合現(xiàn)場(chǎng)絕對(duì)無(wú)阻流量預(yù)測(cè)需要。
【主權(quán)項(xiàng)】
1.基于產(chǎn)氣指數(shù)的頁(yè)巖氣層絕對(duì)無(wú)阻流量的獲取方法，其特征在于:具體步驟為: 1)通過(guò)測(cè)錄井資料獲取工區(qū)已試井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H、水平井段長(zhǎng)度L、壓力系數(shù)Kf、孔隙度POR、含氣飽和度Sg、有機(jī)碳含量TOC和脆性指數(shù)Br i t; 2)獲取已試井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf、壓裂支持劑體積Vs; 3)通過(guò)測(cè)試產(chǎn)能資料獲取已試井頁(yè)巖氣層無(wú)阻流量AOFg; 4)運(yùn)用步驟I)和步驟2)獲取的參數(shù)依次計(jì)算已試井氣層含氣指數(shù)Ig、壓裂改造指數(shù)IFR、產(chǎn)氣指數(shù)IPG，具體公式如下:IFR= (Vf+Vs )/10000，Ig = H X Kf X POR X Sg X TOC， IPG = LX (Ig+IFR) XBrit； 式中:壓裂液體積Vf量綱為m3， 壓裂支持劑體積Vs量綱為m3， 壓裂改造指數(shù)IFR量綱為10km3， 厚度H量綱為m， 水平井段長(zhǎng)度L量綱為hm， 壓力系數(shù)Kf無(wú)量綱， 孔隙度P0R、含氣飽和度Sg、有機(jī)碳含量TOC量綱為小數(shù)， 含脆性指數(shù)Brit量綱為小數(shù)； 5)將獲取的已試井頁(yè)巖氣層無(wú)阻流量AOFg與對(duì)應(yīng)已試井計(jì)算的產(chǎn)氣指數(shù)IPG利用最小二乘法確定模型AOFg=A X eBXIPG^模型系數(shù)A、B; 式中無(wú)阻流量AOFg量綱為104m3/d; 6)通過(guò)測(cè)錄井資料獲取工區(qū)待預(yù)測(cè)井頁(yè)巖氣層參數(shù)，具體參數(shù)為氣層厚度H1、水平井段長(zhǎng)度L 1、壓力系數(shù)Kf 1、孔隙度POR 1、含氣飽和度Sg 1、有機(jī)碳含量TOC I和脆性指數(shù)Brit I; 7)獲取待預(yù)測(cè)井頁(yè)巖氣層壓裂參數(shù)，具體參數(shù)為壓裂液體積Vf1、壓裂支持劑體積VsI; 8)運(yùn)用步驟6)和步驟7)獲取的參數(shù)根據(jù)公式 IFR l = (Vf I+Vs I )/10000， Ig I=H IXKf IXPOR IXSg IXTOC 1， IPG I = L I X (Ig 1+IFR I) XBrit I； 依次計(jì)算待預(yù)測(cè)井氣層壓裂改造指數(shù)IFR 1、含氣指數(shù)Ig 1、產(chǎn)氣指數(shù)IPG1; 9)將步驟5)計(jì)算出的模型系數(shù)A、B帶入AOFgI=AXe6xipc 1中計(jì)算出的待預(yù)測(cè)井頁(yè)巖氣層無(wú)阻流量AOFg I; 10)輸出計(jì)算結(jié)果。
【文檔編號(hào)】G06Q50/02GK105822298SQ201610262650
【公開(kāi)日】2016年8月3日
【申請(qǐng)日】2016年4月25日
【發(fā)明人】馮愛(ài)國(guó), 廖勇, 張建平, 石文睿, 馮詩(shī)旸, 彭超, 趙紅燕, 張占松, 張超謨, 任元, 石元會(huì), 饒海濤, 陳志蓉, 魏煒
【申請(qǐng)人】中石化石油工程技術(shù)服務(wù)有限公司, 中石化江漢石油工程有限公司測(cè)錄井公司