【正文】 to decreasing in the same differential pressure. As is shown, the influence of kickoff pressure gradient on the oil well productivity will be 20% or more, which depends on the size of the kickoff pressure gradient. (4) Multilayers will influence each other. When bottomhole wellbore pressure declines under the minimum of formation pressure, each layer will produce fluid。 油氣田 的產(chǎn)能是油氣田開發(fā)過程中最重要的的問題之一， 是油氣田生產(chǎn)系統(tǒng)分析以及生產(chǎn)動態(tài)預(yù)測的基礎(chǔ)， 是油氣 開發(fā)建設(shè)和合理開發(fā)油氣資源的重要依據(jù)，必須在油氣藏評價階段認(rèn)真研究并加以確定。它為油 井 的合理配產(chǎn)、油田開發(fā)規(guī)劃、開發(fā)方案的制定和生產(chǎn)制度的調(diào)整等重大決策提供科學(xué)的依據(jù)。油田的產(chǎn)能大小是油田開發(fā)的核心問題，如何高效上產(chǎn)和穩(wěn)產(chǎn)是油田開發(fā)的中心環(huán)節(jié)。 通俗的說產(chǎn)能試井就是通過若干次改變氣井的工作制度，并測量在各個工作制度下氣井的穩(wěn)定產(chǎn)量及與之對應(yīng)的井底流壓，根據(jù)產(chǎn)量及對應(yīng)的井底流壓確定氣井的產(chǎn)能方程以及無阻流量。本文主要 對 某一油田的產(chǎn)能 進(jìn)行分析和預(yù)測，然后提供合理的開發(fā)建議。 其中試井產(chǎn)能評價方法最為常用。同登科等 人 [7]建立了應(yīng)力敏感地層 的 雙重介質(zhì)模型， 分析了只有裂縫發(fā)生形變以及裂縫與巖塊同時發(fā)生形變的情況 ,并用 DouglasJones 的方法 預(yù)估－校正法求解了這兩種模型，并 進(jìn)一步 探討了變形參數(shù)和雙重 介質(zhì)參數(shù)變化時壓力的變化規(guī)律。何國良 [11]、同登科等 [12]在建立變形雙重介質(zhì)非線性滲流數(shù)學(xué)模型時，在 WarrenRoot 模型基礎(chǔ)上，引入分形參數(shù)和壓縮系數(shù)，考慮壓力 對具有分形特征的滲透率和孔隙度的影響，并用預(yù)估校正方法進(jìn)行求解； 2020 年，韓洪寶等 [13]通過室內(nèi)基礎(chǔ)實驗和基礎(chǔ)理論研究，以及物理模擬和數(shù)值模擬方法， 分析了由于啟動壓力梯度的存在而使得低滲透儲層產(chǎn)生非達(dá)西滲 流的機理和規(guī)律，建立了描述低滲透油氣非達(dá)西滲流模型。模型既考慮基質(zhì)系統(tǒng)和裂縫系統(tǒng)不同的滲流特性，又考慮了裂縫變形的影響；而且在模型中同時考慮了驅(qū)替、重力和滲吸作用。西南石油大學(xué)的 霍進(jìn)、賈永祿、王海濤等人研究 了具有竄流的有界三層油藏井底壓力動態(tài)，并成功地計算出了典型曲線 ； 2020 年，張烈輝，王海濤，賈永祿等 [16]首次建立了每一層均為 雙孔介質(zhì)的層狀油藏的滲流模型，并成功地繪制出了該模型的典型曲線 ； 馮國慶，劉啟國等 [17]建立了 考慮啟動壓力梯度的低滲透氣藏不穩(wěn)定滲流模型準(zhǔn)確描述存在啟動壓力梯度的低滲透氣藏的不穩(wěn)定滲流問題 。而在 國內(nèi) ， 2020 年，陳元千 [22]教授 基于 Joshi、 Giger 對水平井產(chǎn)量公式的研究思路，利用面積等值的擬圓形驅(qū)動半徑和產(chǎn)量等值的擬圓形生產(chǎn)坑道的物理概念，采用水電相似原理的 等值滲流阻力法，得到了一個新的水平井產(chǎn)量計算公式。故本文研究的主要內(nèi)容 包括 非均質(zhì)條件下產(chǎn)能分析 、 應(yīng)力敏感條件下產(chǎn)能分析 、 啟動壓力條件下產(chǎn)能分析 、 不同井型對產(chǎn)能的影響 、 復(fù)雜儲層下的產(chǎn)能 預(yù)測 ，并進(jìn)行產(chǎn)能評價和預(yù)測。區(qū)域構(gòu)造位置上，油田位于北部灣盆地北部 拗 陷潿西南凹陷 ② 號斷層上升盤，是一個被斷層 復(fù)雜化 的半背斜構(gòu)造。每年 5～ 11 月份為臺風(fēng)季節(jié)，夏季通常 3 級～ 4 級風(fēng)，最大陣風(fēng) 6級～ 7 級，風(fēng)向西南；冬季一般 6 級～ 7 級，最大陣風(fēng) 9 級～ 10 級，風(fēng)向東北。海區(qū)最大潮差達(dá) 。流三段砂巖分選性以中等為主，顆粒磨圓度以次棱～次圓狀為主， 支撐類型以顆粒支撐為主，顆粒間多為點 線狀接觸，孔隙式膠結(jié)。 從開發(fā)井儲層測井解釋物性統(tǒng)計結(jié)果可以看出，整體上 L3Ⅰ 、 L3Ⅱ 油組儲層物性好于 L3Ⅲ 油組。平均意義上， L3Ⅰ 、 L3Ⅱ和 L3Ⅲ A 油組為中孔高滲儲層， L3Ⅲ B油組為中孔中滲儲層 。 油 田分別于 2020 年 5 月 8 日、 2020 年 8 月 3 日分兩批投產(chǎn) 8 口井， 2020 年 8 月實施 2 口調(diào)整井。油田 ODP 設(shè)西南 石油大學(xué)工程碩士 研究生 學(xué)位 論文 5 計建立注采的層位分為兩個層系： L3Ⅰ +Ⅱ ， L3Ⅲ 油組。 2020 年 7月以后， AA7 井正常生產(chǎn)，氣油比趨于平穩(wěn)， 穩(wěn)定在 300m3/m3上下。由此 推測 ，在關(guān)井保護生產(chǎn)期間，原油飽和壓力發(fā)生了變化，地層溶解氣進(jìn)行了重新分配。 潿洲 111 油田生產(chǎn)曲線0300600900120015000705 0711 0804 0810 0904 0910date(yymm)qo(m3/d),ql(m3/d),GOR(m3/m3)020406080100fw(%)ql OPD GOR fw 圖 21 該 油田生產(chǎn)曲線 全油田探明地質(zhì)儲量為 104m3，其中東塊 2 井區(qū) 104m3，西塊 Wan9井區(qū) 104m3，西塊 1 井區(qū) 104m3。截至 2020 年 9 月30 日，全油田 累積采油 104m3。 針對 該油田所在 儲層 的實際 情況， 開展了如下四種模型 的研究： 徑向非均質(zhì) 模型 、線性非均質(zhì) 模型 、塊狀地層部分打開井 模型 、薄互層 模型 。 圖 31 徑向復(fù)合油藏模型 對于外區(qū)，在穩(wěn)態(tài)滲流條件下，如果邊界壓力為 ep ，則地層壓力解可以寫成： e22e2 ln2rqBpp k h r???? (31) 令 frr? 處，得到內(nèi)外區(qū)銜接面處的地層壓力： e2fe2fln2rqBpp k h r???? (32) 對于內(nèi)區(qū)，在穩(wěn)態(tài)滲流條件下的地層壓力解可以寫成： 1f1f1 ln2q B rpp k h r???? (33) 令 wrr? 處，則得到井底流壓： 1fwf1wln2q B rpp k h r???? (34) 內(nèi)內(nèi)區(qū)區(qū) 外外區(qū)區(qū) 111 ??、k222 ??、kerfrwr西南 石油大學(xué)工程碩士 研究生 學(xué)位 論文 7 將（ 32）式代入（ 34）式，得 e2 1 fwe2 f 1 wl n l n22rq B q B rpp k h r k h r????? ? ? (35) 將（ 35）式變形，得 ? ?ewe1 f 21 w 2 f2ln lnppqrB r Bk h r k h r?????? (36) 如果考慮表皮系數(shù)，則上式可以改為： ? ?ewe1 f 21 w 2 f2l n l nppqrB r BSk h r k h r????????????? (37) 式（ 37）就是徑向復(fù)合油藏的穩(wěn)態(tài)產(chǎn)能公式。 q — 油井產(chǎn)量， cm3/s； K — 地層滲透率， um2； ； h — 油層厚度， cm； pe — 供給壓力（地層壓力）， atm； pwf— 井底壓力， atm； ? — 液體的粘度， mPa 參考該油田的參數(shù)范圍，假定 m250?er ， ? ， ? 。從圖上可以看出，當(dāng)復(fù)合油藏的內(nèi)區(qū)與外區(qū)的流度比12M 逐漸增大時，采油指數(shù)之比 Jb 逐漸減小，且隨著流度比 12M 的增大，采油指數(shù)之比Jb 減小的速度越來越慢。 00 5 10 15 20流度比 M 12Jb=Jo復(fù)合/Jorf=25mrf=50mrf=100mrf=150mr e = 2 5 0 m 、 r w = 0 . 0 7 8 5 m 、 S =3 圖 32 流度比 M12 對采油指數(shù)之比 Jb 的影響圖 圖 33 是內(nèi)區(qū)半徑 fr 對采油指數(shù)之比 Jb 的影響圖，從圖上可以看出，當(dāng) M121 時，隨著 fr 的增大， 采油指數(shù)之比 Jb 逐漸減小，且減小的速度越來越慢。當(dāng) M121 時，所有 Jb 值均 1? ，當(dāng) M121，所有 Jb 值均 1。 不穩(wěn)態(tài)產(chǎn)能分析 油氣井的產(chǎn)能不僅與油藏和流體的性質(zhì)有關(guān)，也與壓力波在地層中的傳播范圍有關(guān)，即與測試時間有關(guān) 。 q —— 油井產(chǎn)量， cm3/s； K —— 地層滲透率， um2； ； h —— 油層厚度， cm； pi —— 供給壓力（地層壓力）， atm； pwf—— 井底壓力， atm； ? —— 液體的粘度， mPa 利用 Stehfest 數(shù)值反演編程計算，將拉氏空間的解反演到實空間 得到 圖 34，它 以該油田某 井 資料為基礎(chǔ) 計算 的不 穩(wěn)態(tài) 擬采油指數(shù) Jpo 隨測試時間變化關(guān)系圖，從該 圖 可以看出， 擬 采油指數(shù) Jpo 的大小隨生產(chǎn)測試的時間而下降，在早期 由于壓力波傳播的范圍有限，擬 采油指數(shù) Jpo 下降快 ，在達(dá)到一定時間后，下降速度變緩，變得相對穩(wěn)定。 0204060801000 1000 2020 3000 4000 5000時間（h ）Jpo（m3/(D = 0 . 9 4 5 6 c P 、 B = 1 . 1 6 2 2 、 C t = 0 . 0 0 1 5 6 7 8 M P a1 圖 34 滲透率對不穩(wěn)態(tài)采油指數(shù)的影響圖 0123450 1000 2020 3000 4000 5000時間（h ）擬采油指數(shù)日變化率（%）K=179mDK=100mDK=50mDK=25mDh = 4 . 6 m 、 ? = 0 . 2 9 、 r w = 0 . 1 1 0 2 5 m 、 S = 0 . 7 4181。對于物性較差的儲層，采用 DST測試 [23~25]時，測試時間短，壓力波波及范圍非常有限，初期測試時的產(chǎn)能常常不能反應(yīng)油氣井在生產(chǎn)較長時間下的產(chǎn)能 。 圖 36 是 復(fù)合油藏 模型流度比對 非穩(wěn)態(tài) 比擬采油指數(shù) Jpb 的影響 關(guān)系圖，從該 圖 可以看出， 在內(nèi)區(qū)半徑一定的情況下，比擬 采油指數(shù) Jpb 的大小隨生產(chǎn)測試的時間而下降，同時受內(nèi)外區(qū)流度比的顯著影響。 圖 37 是 復(fù) 合油藏 模型內(nèi)區(qū)半徑對 非穩(wěn)態(tài) 比擬采油指數(shù) Jpb 的影響 關(guān)系圖，從該 圖可以看出， 在流度比一定的情況下，比擬 采油指數(shù) Jpb的大小隨生產(chǎn)測試的時間而下降，同時受內(nèi)區(qū)半徑的影響，相對來說內(nèi)區(qū)半徑影響較小。 02468100 1000 2020 3000 4000 5000時間（h ）Jpb（m3/(Dm)）M12=1M12=2M12=5M12=10k / 181。MPa = 1 0 0 m D / c P 、 ? C t = 0 . 0 0 0 3 M P a1M 12 =5 、 r w = 0 0 . 1 1 0 2 5 m 、 S =0 、 B = 1 . 1 6 2 2 圖 37 內(nèi)區(qū)半徑對比擬采油指數(shù)的影響圖 WZ 油田產(chǎn)能分析與預(yù)測 12 線性非均質(zhì)條件下測試井產(chǎn)能分析 由于 該油田所在區(qū)域 斷塊油藏復(fù) 雜的地質(zhì)特征，強烈的非均質(zhì)性，本次研究還建立了線性非均質(zhì)復(fù)合油藏模型（條帶狀復(fù)合油藏模型）。 圖 38 條帶狀復(fù)合油藏示意圖 對于 Ⅰ 區(qū)的滲流微分方程為： 221 1 122 2 ( ) ( )D D DD w D D w DD D Dp p px x y yx y t? ? ?? ? ?? ? ? ? ?? ? ? ， xD≥0 (316) 對于 Ⅱ 區(qū)的滲流微分方程為： 222 2 222 1D D DD D Dp p px y t?? ? ???? ? ? ， xD≤0 (317) 初始條件： 1 ( , , 0) 0D D Dp x y ?， 2 ( ,0) 0DDpx ? (318) 邊界條件： 1100D D e DDDy y yppyy???? ， 2200D D e DDDy y yppyy???? 1 ( , , ) 0D D Dp y t??， 2 ( , , ) 0D D Dp y t?? ? (319) 界面 銜接條件： 1200DDDDxxppMxx????? ， 1200