低滲透油氣藏的參數(shù)解釋方法及解釋系統(tǒng)的制作方法
【專利摘要】本發(fā)明公開一種低滲透油氣藏的參數(shù)解釋方法及系統(tǒng)，首先建立流體連續(xù)性方程，設(shè)定包含擬啟動壓力梯度的流動方程，利用流體連續(xù)性方程和流動方程建立滲流方程，之后確定低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流體壓縮系數(shù)，設(shè)定低滲透油氣藏的參數(shù)集合，利用低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和參數(shù)集合對滲流方程進行求解，得到井底計算壓力，在井底計算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精度時，當前設(shè)定的參數(shù)集合作為低滲透油氣藏的實際參數(shù)。本發(fā)明公開的解釋方法及系統(tǒng)，能夠更加準確地解釋低滲透油氣藏的參數(shù)。
【專利說明】
低滲透油氣藏的參數(shù)解釋方法及解釋系統(tǒng)
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于油氣藏開發(fā)技術(shù)領(lǐng)域，尤其涉及低滲透油氣藏的參數(shù)解釋方法及解釋 系統(tǒng)。
【背景技術(shù)】
[0002] 流體在低滲透油氣藏中滲流時，必須有一個附加的壓力梯度克服巖石表面吸附膜 或水化膜引起的阻力才能流動，該附加的壓力梯度稱為啟動壓力梯度。對于低滲透油氣藏， 很多學者認為存在絕對啟動壓力梯度，當壓力梯度達到該絕對啟動壓力梯度時，低滲透油 氣藏中的流體開始流動。
[0003] 低滲透油氣藏在試采階段的瞬態(tài)壓力響應(yīng)，普遍具有的特征是:關(guān)井所測的壓力 及壓力的導數(shù)在后期呈平行上翹趨勢。眾多學者對低滲透油氣藏的瞬態(tài)壓力響應(yīng)進行研 究，認為低滲透油氣藏的絕對啟動壓力梯度是導致壓力及壓力的導數(shù)在后期平行上翹的原 因。據(jù)此，本領(lǐng)域技術(shù)人員基于疊加原理，對包含絕對啟動壓力梯度的滲流方程進行求解， 以確定低滲透油氣藏的參數(shù)。
[0004] 但是，
【申請人】發(fā)現(xiàn)基于上述方式解釋出的低滲透油氣藏的參數(shù)存在很大的偏差。

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

[0005] 有鑒于此，本發(fā)明的目的在于提供一種準確度更高的低滲透油氣藏的參數(shù)解釋方 法及解釋系統(tǒng)。
[0006] 為實現(xiàn)上述目的，本發(fā)明提供如下技術(shù)方案：
[0007] 本發(fā)明提供一種低滲透油氣藏的參數(shù)解釋方法，包括：
[0008] 根據(jù)所述低滲透油氣藏的流體特性，建立流體連續(xù)性方程；
[0009] 設(shè)定流動方程，所述流動方程中的待擬合參數(shù)至少包括擬啟動壓力梯度、最大啟 動壓力梯度、最小啟動壓力梯度和絕對滲透率；
[0010] 利用所述流體連續(xù)性方程和所述流動方程建立滲流方程；
[0011] 確定所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流體壓縮系數(shù)；
[0012] 設(shè)定所述低滲透油氣藏的參數(shù)集合，所述參數(shù)集合包括所述低滲透油氣藏的孔隙 度、地層的有效厚度、原始地層壓力、表皮系數(shù)、井儲系數(shù)、油藏邊界、以及所述流動方程包 含的待擬合參數(shù)；
[0013] 利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和當 前設(shè)定的參數(shù)集合對所述滲流方程進行求解，獲得井底計算壓力；
[0014] 對所述井底計算壓力與所述低滲透油氣藏的井底實測壓力進行擬合，得到擬合精 度；
[0015] 如果所述擬合精度達到預(yù)設(shè)精度，則確定當前設(shè)定的參數(shù)集合為所述低滲透率油 氣藏的實際參數(shù)；
[0016] 如果所述擬合精度未達到所述預(yù)設(shè)精度，則調(diào)整所述低滲透油氣藏的參數(shù)集合， 執(zhí)行利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和當前設(shè) 定的參數(shù)集合對所述滲流方程進行求解的步驟及后續(xù)步驟。
[0017]優(yōu)選的，上述方法中，所述流動方程為：
[0019] 其中，u為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓 力梯度，X min為最小啟動壓力梯度，Ape3sud。為擬啟動壓力梯度，η為大于1的整數(shù)。
[0020] 優(yōu)選的，上述方法中，所述調(diào)整所述低滲透油氣藏的參數(shù)集合，包括：
[0021] 步驟S91:判斷井儲階段的井底計算壓力與井底實測壓力的擬合精度是否達到所 述預(yù)設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S92，若達到所述預(yù)設(shè)精度，則執(zhí)行步驟 S93；
[0022]步驟S92:調(diào)整井儲系數(shù)；
[0023]步驟S93:判斷井儲向地層的過渡階段的井底計算壓力與井底實測壓力的擬合精 度是否達到所述預(yù)設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S94，若達到所述預(yù)設(shè)精度， 則執(zhí)行步驟S96;
[0024]步驟S94:調(diào)整第一參數(shù)子集，所述第一參數(shù)子集包括最小啟動壓力梯度、最大啟 動壓力梯度、擬啟動壓力梯度、絕對滲透率和表皮系數(shù)；
[0025]步驟S95:如果調(diào)整第一參數(shù)子集的次數(shù)滿足第一預(yù)定條件，仍不能使得井儲向地 層的過渡階段的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則微調(diào)井儲 系數(shù)；
[0026]步驟S96:判斷井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬合 精度是否達到所述預(yù)設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S97，若達到所述預(yù)設(shè)精 度，則轉(zhuǎn)至步驟S99;
[0027]步驟S97:調(diào)整第二參數(shù)子集，所述第二參數(shù)子集包括最小啟動壓力梯度、最大啟 動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界和孔隙度；
[0028]步驟S98:如果調(diào)整第二參數(shù)子集的次數(shù)滿足第二預(yù)定條件，仍不能使得井儲向地 層過渡階段之后的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則微調(diào)井 儲系數(shù)和表皮系數(shù)；
[0029]步驟S99:判斷整體的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù) 設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S910;
[0030]步驟S910:微調(diào)第三參數(shù)子集，所述第三參數(shù)子集包括最小啟動壓力梯度、最大啟 動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界、孔隙度、表皮系數(shù)和井儲系數(shù)；
[0031] 步驟S911:如果微調(diào)第三參數(shù)子集的次數(shù)滿足第三預(yù)定條件，仍不能使得整體的 井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則調(diào)整所述流動方程中的指 數(shù)η。
[0032] 優(yōu)選的，上述方法中，所述對所述井底計算壓力與所述低滲透油氣藏的井底實測 壓力進行擬合，包括：
[0033]繪制井底計算壓力隨時間變化的第一曲線；
[0034]繪制井底計算壓力的導數(shù)隨時間變化的第二曲線；
[0035] 對所述第一曲線與第三曲線進行擬合，對所述第二曲線與第四曲線進行擬合，得 到擬合精度，其中所述第三曲線為井底實測壓力隨時間變化的曲線，所述第四曲線為所述 井底實測壓力的導數(shù)隨時間變化的曲線。
[0036] 本發(fā)明還提供一種低滲透油氣藏的參數(shù)解釋系統(tǒng)，包括：
[0037] 流體連續(xù)性方程建立單元，用于根據(jù)所述低滲透油氣藏的流體特性，建立流體連 續(xù)性方程；
[0038] 流動方程設(shè)定單元，用于設(shè)定流動方程，所述流動方程中的待擬合參數(shù)至少包括 擬啟動壓力梯度、最大啟動壓力梯度、最小啟動壓力梯度和絕對滲透率；
[0039] 滲流方程建立單元，用于利用所述流體連續(xù)性方程和所述流動方程建立滲流方 程；
[0040] 參數(shù)確定單元，用于確定所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù) 和流體壓縮系數(shù)；
[0041] 參數(shù)設(shè)定單元，用于設(shè)定所述低滲透油氣藏的參數(shù)集合，所述參數(shù)集合包括所述 低滲透油氣藏的孔隙度、地層的有效厚度、原始地層壓力、表皮系數(shù)、井儲系數(shù)、油藏邊界、 以及所述流動方程包含的待擬合參數(shù)；
[0042] 計算單元，用于利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流 體壓縮系數(shù)和當前設(shè)定的參數(shù)集合對所述滲流方程進行求解，獲得井底計算壓力；
[0043]擬合單元，用于對所述井底計算壓力與所述低滲透油氣藏的井底實測壓力進行擬 合，得到擬合精度；
[0044] 第一處理單元，用于在所述擬合單元得到的擬合精度達到預(yù)設(shè)精度的情況下，確 定當前設(shè)定的參數(shù)集合為所述低滲透率油氣藏的實際參數(shù)；
[0045] 第二處理單元，用于在所述擬合單元得到的擬合精度未達到所述預(yù)設(shè)精度的情況 下，調(diào)整所述低滲透油氣藏的參數(shù)集合，觸發(fā)所述計算單元執(zhí)行操作。
[0046] 優(yōu)選的，上述系統(tǒng)中，所述流動方程設(shè)定單元設(shè)定的流動方程為：
[0048] 其中，u為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓 力梯度，X min為最小啟動壓力梯度，Ape3sud。為擬啟動壓力梯度，η為大于1的整數(shù)。
[0049] 優(yōu)選的，上述系統(tǒng)中，所述第二處理單元包括：
[0050] 第一判斷子單元，用于判斷井儲階段的井底計算壓力與井底實測壓力的擬合精度 是否達到所述預(yù)設(shè)精度；
[0051] 第一處理子單元，用于在所述第一判斷子單元確定井儲階段的井底計算壓力與井 底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整井儲系數(shù)；
[0052]第二判斷子單元，用于在所述第一判斷子單元確定井儲階段的井底計算壓力與井 底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層的過渡階段的井底計 算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0053]第二處理子單元，用于在所述第二判斷子單元確定井儲向地層的過渡階段的井底 計算壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第一參數(shù)子集， 所述第一參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透 率和表皮系數(shù)；
[0054]第三處理子單元，用于在所述第二處理子單元執(zhí)行調(diào)整第一參數(shù)子集的次數(shù)滿足 第一預(yù)定條件，仍不能使得井儲向地層的過渡階段的井底計算壓力與井底實測壓力的擬合 精度達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)；
[0055]第三判斷子單元，用于在所述第二判斷子單元確定井儲向地層的過渡階段的井底 計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層過渡階 段之后的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0056]第四處理子單元，用于在所述第三判斷子單元確定井儲向地層過渡階段之后的井 底計算壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第二參數(shù)子 集，所述第二參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對 滲透率、油藏邊界和孔隙度；
[0057]第五處理子單元，用于在所述第四處理子單元執(zhí)行調(diào)整第二參數(shù)子集的次數(shù)滿足 第二預(yù)定條件，仍不能使得井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬 合精度達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)和表皮系數(shù)；
[0058]第四判斷子單元，用于在所述第三判斷子單元確定井儲向地層過渡階段之后的井 底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷整體的井底計算 壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0059] 第六處理子單元，用于在所述第四判斷子單元確定整體的井底計算壓力與井底實 測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，微調(diào)第三參數(shù)子集，所述第三參數(shù)子集 包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油臧邊界、孔隙 度、表皮系數(shù)和井儲系數(shù)；
[0060] 第七處理子單元，用于在所述第六處理子單元執(zhí)行調(diào)整第三參數(shù)子集的次數(shù)滿足 第三預(yù)定條件，仍不能使得整體的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè) 精度的情況下，調(diào)整所述流動方程中的指數(shù)η。
[0061] 優(yōu)選的，上述系統(tǒng)中，所述擬合單元具體用于：
[0062] 繪制井底計算壓力隨時間變化的第一曲線；
[0063]繪制井底計算壓力的導數(shù)隨時間變化的第二曲線；
[0064]對所述第一曲線與第三曲線進行擬合，對所述第二曲線與第四曲線進行擬合，得 到擬合精度，其中所述第三曲線為井底實測壓力隨時間變化的曲線，所述第四曲線為所述 井底實測壓力的導數(shù)隨時間變化的曲線。
[0065]由此可見，本發(fā)明的有益效果為：
[0066] 本發(fā)明公開的低滲透油氣藏的參數(shù)解釋方法和解釋系統(tǒng)，設(shè)定的流動方程包含有 擬啟動壓力梯度，由于擬啟動壓力梯度是導致低滲透油氣藏的壓力及壓力的導數(shù)在后期平 行上翹的原因，因此包含擬啟動壓力梯度的流動方程能夠更加準確地反映低滲透油氣藏的 流動規(guī)律;相應(yīng)的，對利用流體連續(xù)性方程和包含有擬啟動壓力梯度的流動方程建立的滲 流方程進行求解，在低滲透油氣藏的井底計算壓力和井底實測壓力的擬合精度達到預(yù)設(shè)精 度情況下，所設(shè)定的低滲透油氣藏的參數(shù)集合，更加接近低滲透油氣藏的實際參數(shù)，具有更 高的準確度。
【附圖說明】
[0067] 為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動的前提下，還可以根據(jù) 提供的附圖獲得其他的附圖。
[0068] 圖1為置加原理的不意圖；
[0069] 圖2為本發(fā)明公開的一種低滲透油氣藏的參數(shù)解釋方法的流程圖；
[0070] 圖3為低滲透油氣藏中流體的速度與壓力梯度之間的關(guān)系圖；
[0071] 圖4為本發(fā)明公開的一種調(diào)整低滲透油氣藏的參數(shù)集合的流程圖；
[0072] 圖5為本發(fā)明公開的一種低滲透油氣藏的參數(shù)解釋系統(tǒng)的結(jié)構(gòu)框圖；
[0073]圖6為圖5所示參數(shù)解釋系統(tǒng)中第二處理單元的一個結(jié)構(gòu)框圖。
【具體實施方式】
[0074] 目前，本領(lǐng)域技術(shù)人員基于疊加原理，對包含絕對啟動壓力梯度的滲流方程進行 求解，以確定低滲透油氣藏的參數(shù)。
【申請人】經(jīng)過大量研究，發(fā)現(xiàn)目前的處理方式是錯誤的。
[0075] 對于齊次的滲流方程，可用疊加原理進行變流量的壓力響應(yīng)計算。
[0076]這里舉例進行說明，以流量q從to時刻到4時刻生產(chǎn)、從4時刻關(guān)井且關(guān)井時間為 A t的流量史等效為：以流量q從to時刻到(tP+A t)時刻生產(chǎn)的流量史(如圖1中的曲線1)，與 以流量_q從tP時刻開始注入直到(tP+ △ t)時刻的流量史(如圖1中的曲線2)的疊加。圖1中的 曲線3為以流量q從to時刻到tP時刻生產(chǎn)、從t P時刻關(guān)井且關(guān)井時間為△ t的流量史。
[0077]但是，包含絕對啟動壓力梯度的滲流方程是非齊次的。
[0078] 公式1為低滲透油氣藏的包含絕對啟動壓力梯度的流動方程：
[0079]
[0080]其中，U為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，λ為絕對啟動壓力 梯度。
[0081 ]將公式1帶入流體連續(xù)性方程，就得到包含絕對啟動壓力梯度的滲流方程。例如， 將公式1帶入單相流體的流體連續(xù)性方程(公式2)，就可以得到包含絕對啟動壓力梯度的滲 流方程(公式3)。
[0082]
[0083] 其中，Ρ為密度，單位為Kg/m3，u為速度，單位為m/s，Φ為孔隙度，f為井的質(zhì)量流 量，單位為Kg/(m 3s)。
[0084]
[0085] 可以看到，在公式3中，由于存在λ項，因此該滲流方程不是齊次方程，如果基于疊 加原理對包含絕對啟動壓力梯度的滲流方程進行求解，就會產(chǎn)生錯誤的結(jié)果。
[0086] 本發(fā)明公開一種準確度更高的低滲透油氣藏的參數(shù)解釋方法及解釋系統(tǒng)。
[0087] 下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進行清楚、完 整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例?；?本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他 實施例，都屬于本發(fā)明保護的范圍。
[0088] 參見圖2,圖2為本發(fā)明公開的一種低滲透油氣藏的參數(shù)解釋方法的流程圖。該方 法包括：
[0089] 步驟S1:根據(jù)低滲透油氣藏的流體特性，建立流體連續(xù)性方程。
[0090]流體連續(xù)性方程是質(zhì)量守恒定律在流體力學中的具體表述形式。上文中的公式2 即為一種單相流體的流體連續(xù)性方程。
[0091] 步驟S2:設(shè)定流動方程，流動方程中的待擬合參數(shù)至少包括擬啟動壓力梯度、最大 啟動壓力梯度、最小啟動壓力梯度和絕對滲透率。
[0092]
【申請人】經(jīng)過大量研究發(fā)現(xiàn)，低滲透油氣藏的擬啟動壓力梯度是導致低滲透油氣藏 的壓力及壓力的導數(shù)在后期平行上翹的原因。
[0093]低滲透油氣藏的滲流是由邊界層、細微的喉道、非牛頓流體等因素共同形成的。儲 集層的小喉道所對應(yīng)的啟動壓力梯度較大，而大喉道所對應(yīng)的啟動壓力梯度較小。對于真 實的油氣藏，由于儲集層中喉道的大小分布在一個區(qū)間內(nèi)，因此油氣藏的啟動壓力梯度也 分布在一個區(qū)間內(nèi)。儲集層中的最小喉道對應(yīng)于最大啟動壓力梯度，即A max;儲集層中的最 大喉道對應(yīng)于最小啟動壓力梯度，即xmin。若λΜη=0,則意味著最大喉道不存在啟動壓力梯 度。
[0094]對于低滲透油氣藏，當壓力梯度從最小啟動壓力梯度λΜη逐漸增大至最大啟動壓 力梯度λΜΧ時，將有更多的小喉道的流體參與流動，流體的速度與壓力梯度呈非線性關(guān)系， 不滿足達西定律；當壓力梯度大于最大啟動壓力梯度A max時，流體的速度與壓力梯度呈線性 關(guān)系，即滿足達西定律。
[0095]圖3示出了低滲透油氣藏中流體的速度與壓力梯度的關(guān)系。可以看到，當壓力梯度 大于最大啟動壓力梯度λ_χ時，流體的速度與壓力梯度之間是線性關(guān)系。當壓力梯度介于最 大啟動壓力梯度λ_ χ與最小啟動壓力梯度Amax之間時，流體的速度與壓力梯度之間是非線性 關(guān)系。
[0096]也就是說，低滲透油氣藏的儲集層的喉道微細，在低速滲流時不符合線性滲流規(guī) 律，滲流的速度(也就是流體的速度)和壓力梯度的關(guān)系是一條曲線，而不是直線。當滲流的 速度增加到一定程度時，滲流的速度和壓力梯度的關(guān)系才能變成一條直線，但是該直線不 再通過原點，將該直線延長與壓力梯度坐標軸相交，該直線的延長線在壓力梯度坐標軸的 截距稱為擬啟動壓力梯度。相對于絕對啟動壓力梯度理論，人們稱此為擬啟動壓力梯度理 論。擬啟動壓力梯度是儲集層的孔隙結(jié)構(gòu)和固液作用的綜合體現(xiàn)，表征了儲集層滲流的非 線性程度和滲流能力，而不是滲流的門檻壓力梯度。
[0097]由于導致低滲透油氣藏的壓力及壓力的導數(shù)在后期平行上翹的原因是:低滲透油 氣藏的擬啟動壓力梯度，因此，包含擬啟動壓力梯度的流動方程，才能更加準確地反映低滲 透油氣藏的流動規(guī)律。
[0098]步驟S3:利用流體連續(xù)性方程和流動方程建立滲流方程。
[0099]流動方程是關(guān)于求解流體的速度u的公式，將流動方程帶入流體連續(xù)性方程，就可 以獲得滲流方程。
[0100]步驟S4:確定低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流體壓縮系 數(shù)。
[0101]其中，低滲透油氣藏的井筒類型包括:垂直井、垂直裂縫井、水平井和水平壓裂井 等。流體性質(zhì)是指:地下油、地下水等流體的粘度和地層體積系數(shù)隨壓力變化情況。低滲透 油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流體壓縮系數(shù)為已知量。
[0102] 步驟S5:設(shè)定低滲透油氣藏的參數(shù)集合，該參數(shù)集合包括低滲透油氣藏的孔隙度、 地層的有效厚度、原始地層壓力、表皮系數(shù)、井儲系數(shù)、油藏邊界、以及流動方程包含的待擬 合參數(shù)。
[0103] 其中，低滲透油氣藏的原始地層壓力是指:壓裂前的地層壓力分布。表皮系數(shù)是 指:用于描述近井地層的滲透率變化的系數(shù)。井儲系數(shù)是指:用于描述近井筒續(xù)流效應(yīng)的系 數(shù)。油藏邊界包括邊界類型(如定壓邊界和封閉邊界)和邊界大小(也就是油氣藏大?。?。
[0104] 實施中，如果低滲透油氣藏包含多相流體，該參數(shù)集合還包括多相流體中各相流 體的飽和度。
[0105] 步驟S6:利用低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù) 和當前設(shè)定的參數(shù)集合對滲流方程進行求解，獲得井底計算壓力。
[0106] 利用低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和當前設(shè) 定的參數(shù)集合對滲流方程進行求解得到井底壓力，該井底壓力是計算得到的，將其記為井 底計算壓力，以便與實測得到的井底壓力(記為井底實測壓力)進行區(qū)分。
[0107] 實施中，對滲流方程進行求解可以為數(shù)值求解、解析求解或者半數(shù)值半解析求解。
[0108] 步驟S7:對井底計算壓力與低滲透油氣藏的井底實測壓力進行擬合，得到擬合精 度。
[0109] 通過執(zhí)行步驟S6得到了低滲透油氣藏的井底計算壓力，同時獲取該低滲透油氣藏 的井底實測壓力，對井底計算壓力和井底實測壓力進行擬合，之后根據(jù)擬合精度執(zhí)行后續(xù) 的步驟S8或者步驟S9。
[0110] 步驟S8:如果擬合精度達到預(yù)設(shè)精度，則確定當前設(shè)定的參數(shù)集合為低滲透率油 氣藏的實際參數(shù)。
[0111] 步驟S9:如果擬合精度未達到預(yù)設(shè)精度，則重新設(shè)定低滲透油氣藏的參數(shù)集合，返 回步驟S6。
[0112]如果井底計算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精度，表明當前設(shè)定的低 滲透油氣藏的參數(shù)集合與低滲透油氣藏的真實參數(shù)非常接近。因此，當前設(shè)定的參數(shù)集合 即為低滲透油氣藏的實際參數(shù)。
[0113]如果井底計算壓力與井底實測壓力的擬合精度未達到預(yù)設(shè)精度，表明當前設(shè)定的 參數(shù)集合中的至少一個參數(shù)存在偏差，因此需要調(diào)整低滲透油氣藏的參數(shù)集合，返回執(zhí)行 步驟S6,直至對滲流方程求解得到的井底計算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精 度。
[0114] 本發(fā)明公開的低滲透油氣藏的參數(shù)解釋方法，設(shè)定的流動方程包含有擬啟動壓力 梯度，由于擬啟動壓力梯度是導致低滲透油氣藏的壓力及壓力的導數(shù)在后期平行上翹的原 因，因此包含擬啟動壓力梯度的流動方程能夠更加準確地反映低滲透油氣藏的流動規(guī)律； 相應(yīng)的，對利用流體連續(xù)性方程和包含有擬啟動壓力梯度的流動方程建立的滲流方程進行 求解，在低滲透油氣藏的井底計算壓力和井底實測壓力的擬合精度達到預(yù)設(shè)精度情況下， 所設(shè)定的低滲透油氣藏的參數(shù)集合，更加接近低滲透油氣藏的實際參數(shù)，具有更高的準確 度。
[0115] 在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋方法中，流動方程可以為：
[0116]
[0117]其中，U為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓 力梯度，Xmin為最小啟動壓力梯度，Ape3sud。為擬啟動壓力梯度，a、b和c為系數(shù)。
[0118] 作為優(yōu)選方案，在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋方法中，流動方程 為：
[0119]
[0120] 其中，U為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓 力梯度，X min為最小啟動壓力梯度，Ape3sud。為擬啟動壓力梯度，η為大于1的整數(shù)，η越大表示流 體的速度越慢。
[0121] 本發(fā)明公開的如公式5所示的流動方程，具有以下優(yōu)點：1、該流動方程中壓力梯度 的取值是連續(xù)的；2、該流動方程中除了最大啟動壓力梯度A max、最小啟動壓力梯度Amin、擬啟 動壓力梯度λΡ_<!。之外，僅引入變量η，使得應(yīng)用更加方便，并且能夠減小整個參數(shù)解釋過程 的運算量。
[0122] 在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋方法中，調(diào)整低滲透油氣藏的參數(shù) 集合，可以采用如圖4所示的方式，包括：
[0123] 步驟S91:判斷井儲階段的井底計算壓力與井底實測壓力的擬合精度是否達到預(yù) 設(shè)精度，若未達到預(yù)設(shè)精度，則執(zhí)行步驟S92，若達到預(yù)設(shè)精度，則執(zhí)行步驟S93。
[0124] 步驟S92:調(diào)整井儲系數(shù)，執(zhí)行步驟S6。
[0125] 步驟S93:判斷井儲向地層的過渡階段的井底計算壓力與井底實測壓力的擬合精 度是否達到預(yù)設(shè)精度，若未達到預(yù)設(shè)精度，則執(zhí)行步驟S94，若達到預(yù)設(shè)精度，則執(zhí)行步驟 S96〇
[0126] 步驟S94:調(diào)整第一參數(shù)子集，執(zhí)行步驟S6。其中，第一參數(shù)子集包括最小啟動壓力 梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率和表皮系數(shù)。
[0127] 步驟S95:如果調(diào)整第一參數(shù)子集的次數(shù)滿足第一預(yù)定條件，仍不能使得井儲向地 層的過渡階段的井底計算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精度，則微調(diào)井儲系 數(shù)，執(zhí)行步驟S6。
[0128] 其中，該第一預(yù)定條件可以為:調(diào)整第一參數(shù)子集的次數(shù)達到第一預(yù)定次數(shù)。該第 一預(yù)定條件也可以為:調(diào)整第一參數(shù)子集的次數(shù)已使得第一參數(shù)子集中各參數(shù)的取值遍歷 各參數(shù)的理論取值區(qū)間。
[0129] 步驟S96:判斷井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬合 精度是否達到預(yù)設(shè)精度，若未達到預(yù)設(shè)精度，則執(zhí)行步驟S97，若達到預(yù)設(shè)精度，則轉(zhuǎn)至步驟 S99〇
[0130] 步驟S97:調(diào)整第二參數(shù)子集，執(zhí)行步驟S6。其中，第二參數(shù)子集包括最小啟動壓力 梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油臧邊界和孔隙度。
[0131] 步驟S98:如果調(diào)整第二參數(shù)子集的次數(shù)滿足第二預(yù)定條件，仍不能使得井儲向地 層過渡階段之后的井底計算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精度，則微調(diào)井儲系 數(shù)和表皮系數(shù)，執(zhí)行步驟S6。
[0132] 其中，該第二預(yù)定條件可以為:調(diào)整第二參數(shù)子集的次數(shù)達到第二預(yù)定次數(shù)。該第 二預(yù)定條件也可以為:調(diào)整第二參數(shù)子集的次數(shù)已使得第二參數(shù)子集中各參數(shù)的取值遍歷 各參數(shù)的理論取值區(qū)間。
[0133] 步驟S99:判斷整體的井底計算壓力與井底實測壓力的擬合精度是否達到預(yù)設(shè)精 度，若達到預(yù)設(shè)精度，則確定當前設(shè)定的參數(shù)集合為低滲透率油氣藏的實際參數(shù)，若未達到 預(yù)設(shè)精度，則執(zhí)行步驟S910。
[0134] 步驟S910:微調(diào)第三參數(shù)子集，執(zhí)行步驟S6。其中，第三參數(shù)子集包括最小啟動壓 力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界、孔隙度、表皮系數(shù)和 井儲系數(shù)。
[0135] 步驟S911:如果微調(diào)第三參數(shù)子集的次數(shù)滿足第三預(yù)定條件，仍不能使得井底計 算壓力與井底實測壓力的擬合精度達到預(yù)設(shè)精度，則調(diào)整流動方程中的指數(shù)n，執(zhí)行步驟 S6〇
[0136] 其中，該第三預(yù)定條件可以為:調(diào)整第三參數(shù)子集的次數(shù)達到第三預(yù)定次數(shù)。該第 三預(yù)定條件也可以為：調(diào)整第三參數(shù)子集的次數(shù)已使得第三參數(shù)子集中各參數(shù)的取值遍歷 各參數(shù)的理論取值區(qū)間。
[0137] 低滲透油氣藏的開采過程通常包括以下階段:井儲階段-井儲向地層過渡階段-徑 向流階段-邊界主導的流動階段。在一些情況下，徑向流階段和邊界主導的流動階段可能不 出現(xiàn)。
[0138] 在井儲階段，低滲透油氣藏的壓力主要是受井儲系數(shù)的影響。因此，在對低滲透油 氣藏的井底計算壓力和井底實測壓力進行擬合之后，如果井儲階段的井底計算壓力與井底 實測壓力的擬合精度達到預(yù)設(shè)精度，表明當前設(shè)定的井儲系數(shù)非常接近低滲透油氣藏的實 際井儲系數(shù);如果井儲階段的井底計算壓力與井底實測壓力的擬合精度未達到預(yù)設(shè)精度， 則調(diào)整井儲系數(shù)。
[0139] 在井儲向地層過渡階段，低滲透油氣藏的壓力主要是受最小啟動壓力梯度、最大 啟動壓力梯度、擬啟動壓力梯度、絕對滲透率和表皮系數(shù)的影響。因此，在對低滲透油氣藏 的井底計算壓力和井底實測壓力進行擬合之后，如果井儲向地層過渡階段的井底計算壓力 與井底實測壓力的擬合精度未達到預(yù)設(shè)精度，則調(diào)整最小啟動壓力梯度、最大啟動壓力梯 度、擬啟動壓力梯度、絕對滲透率和表皮系數(shù)。
[0140] 這里需要說明的是，如果低滲透油氣藏包含多相流體，那么還需要調(diào)整多相流體 中各相流體的飽和度。
[0141]另外，如果多次調(diào)整最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕 對滲透率和表皮系數(shù)之后，井儲向地層過渡階段的井底計算壓力與井底實測壓力的擬合精 度未達到預(yù)設(shè)精度，可以微調(diào)井儲系數(shù)。
[0142]在井儲向地層過渡階段之后（如徑向流階段、邊界主導的流動階段），低滲透油氣 藏的壓力主要是受最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、 油藏邊界和孔隙度的影響。因此，在對低滲透油氣藏的井底計算壓力和井底實測壓力進行 擬合之后，如果井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬合精度未達 到預(yù)設(shè)精度，則調(diào)整最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、 油藏邊界和孔隙度。
[0143]實施中，如果多次調(diào)整最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、 絕對滲透率、油藏邊界和孔隙度之后，井儲向地層過渡階段之后的井底計算壓力與井底實 測壓力的擬合精度仍未達到預(yù)設(shè)精度，可以微調(diào)井儲系數(shù)和表皮系數(shù)。
[0144] 如果井儲階段、井儲向地層過渡階段、以及井儲向地層過渡階段之后的井底計算 壓力與井底實測壓力的擬合精度均達到預(yù)設(shè)精度，需要進一步判斷整體的井底計算壓力與 井底實測壓力的擬合精度是否達到預(yù)設(shè)精度，如果未達到預(yù)設(shè)精度，則微調(diào)最小啟動壓力 梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界、孔隙度、表皮系數(shù)和井 儲系數(shù)。如果多次微調(diào)最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透 率、油藏邊界、孔隙度、表皮系數(shù)和井儲系數(shù)之后，整體的井底計算壓力與井底實測壓力的 擬合精度仍未達到預(yù)設(shè)精度，需要調(diào)整流動方程中的其他待擬合參數(shù):針對公式5所示的流 動方程，需要調(diào)整該流動方程中的指數(shù)n，針對公式4所示的流動方程，需要調(diào)整系數(shù)a、b和 c〇
[0145] 本發(fā)明中所稱的微調(diào)參數(shù)是指：以較小的步長調(diào)整該參數(shù)的取值。
[0146] 本發(fā)明上述公開的調(diào)整低滲透油氣藏的參數(shù)集合的方法，根據(jù)井底計算壓力和井 底實測壓力在各階段的擬合精度，有針對性地調(diào)整特定的參數(shù)，能夠盡快完成整個參數(shù)解 釋過程。
[0147] 另外，本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋方法中，對井底計算壓力與低 滲透油氣藏的井底實測壓力進行擬合可以包括壓力史擬合、壓力降落或恢復擬合、以及壓 力導數(shù)擬合中的一項或多項。
[0148] 作為優(yōu)選方案，對井底計算壓力與低滲透油氣藏的井底實測壓力進行擬合，采用 如下方案：
[0149] 繪制井底計算壓力隨時間變化的第一曲線；
[0150]繪制井底計算壓力的導數(shù)隨時間變化的第二曲線；
[0151] 對第一曲線與第三曲線進行擬合，對第二曲線與第四曲線進行擬合，得到擬合精 度。其中，第三曲線為井底實測壓力隨時間變化的曲線，第四曲線為井底實測壓力的導數(shù)隨 時間變化的曲線。
[0152] 本發(fā)明上述公開了低滲透油氣藏的參數(shù)解釋方法，相應(yīng)的，本發(fā)明公開低滲透油 氣藏的參數(shù)解釋系統(tǒng)。下文關(guān)于參數(shù)解釋系統(tǒng)的描述，與上文關(guān)于參數(shù)解釋方法的描述，可 以相互參見。
[0153] 參見圖5,圖5為本發(fā)明公開的一種低滲透油氣藏的參數(shù)解釋系統(tǒng)的結(jié)構(gòu)框圖。該 參數(shù)解釋系統(tǒng)包括流體連續(xù)性方程建立單元10、流動方程設(shè)定單元20、滲流方程建立單元 30、參數(shù)確定單元40、參數(shù)設(shè)定單元50、計算單元60、擬合單元70、第一處理單元80和第二處 理單元90。
[0154] 其中：
[0155] 流體連續(xù)性方程建立單元10,用于根據(jù)低滲透油氣藏的流體特性，建立流體連續(xù) 性方程。
[0156] 流動方程設(shè)定單元20,用于設(shè)定流動方程，流動方程中的待擬合參數(shù)至少包括擬 啟動壓力梯度、最大啟動壓力梯度、最小啟動壓力梯度和絕對滲透率。
[0157] 滲流方程建立單元30,用于利用流體連續(xù)性方程和流動方程建立滲流方程。
[0158] 參數(shù)確定單元40,用于確定低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和 流體壓縮系數(shù)。
[0159] 參數(shù)設(shè)定單元50,用于設(shè)定低滲透油氣藏的參數(shù)集合。
[0160] 其中，參數(shù)集合包括低滲透油氣藏的孔隙度、地層的有效厚度、原始地層壓力、表 皮系數(shù)、井儲系數(shù)、油藏邊界、以及流動方程包含的待擬合參數(shù)。在低滲透油氣藏包含多相 流體時，該參數(shù)集合還包括多相流體中各相流體的飽和度。
[0161] 計算單元60,用于利用低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體 壓縮系數(shù)和當前設(shè)定的參數(shù)集合對滲流方程進行求解，獲得井底計算壓力。
[0162] 擬合單元70,用于對井底計算壓力與低滲透油氣藏的井底實測壓力進行擬合，得 到擬合精度。
[0163] 第一處理單元80,用于在擬合單元70得到的擬合精度達到預(yù)設(shè)精度的情況下，確 定當前設(shè)定的參數(shù)集合為低滲透率油氣藏的實際參數(shù)。
[0164] 第二處理單元90,用于在擬合單元70得到的擬合精度未達到預(yù)設(shè)精度的情況下， 調(diào)整低滲透油氣藏的參數(shù)集合，觸發(fā)計算單元60執(zhí)行操作。
[0165] 本發(fā)明公開的低滲透油氣藏的參數(shù)解釋系統(tǒng)，設(shè)定的流動方程包含有擬啟動壓力 梯度，由于擬啟動壓力梯度是導致低滲透油氣藏的壓力及壓力的導數(shù)在后期平行上翹的原 因，因此包含擬啟動壓力梯度的流動方程能夠更加準確地反映低滲透油氣藏的流動規(guī)律； 相應(yīng)的，對利用流體連續(xù)性方程和包含有擬啟動壓力梯度的流動方程建立的滲流方程進行 求解，在低滲透油氣藏的井底計算壓力和井底實測壓力的擬合精度達到預(yù)設(shè)精度情況下， 所設(shè)定的低滲透油氣藏的參數(shù)集合，更加接近低滲透油氣藏的實際參數(shù)，具有更高的準確 度。
[0166] 在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋系統(tǒng)中，作為一種優(yōu)選方案，流動 方程設(shè)定單元20設(shè)定的流動方程為上文中的公式5。當然，流動方程設(shè)定單元20設(shè)定的流動 方程并不限定為公式5。
[0167] 在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋系統(tǒng)中，作為一種優(yōu)選方案，第二 處理單元90的結(jié)構(gòu)如圖6所示，包括：
[0168] 第一判斷子單元901，用于判斷井儲階段的井底計算壓力與井底實測壓力的擬合 精度是否達到所述預(yù)設(shè)精度；
[0169] 第一處理子單元902,用于在所述第一判斷子單元確定井儲階段的井底計算壓力 與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整井儲系數(shù)；
[0170]第二判斷子單元903,用于在所述第一判斷子單元901確定井儲階段的井底計算壓 力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層的過渡階段的 井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0171]第二處理子單元904,用于在所述第二判斷子單元903確定井儲向地層的過渡階段 的井底計算壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第一參數(shù) 子集，所述第一參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕 對滲透率和表皮系數(shù)；
[0172] 第三處理子單元905,用于在所述第二處理子單元904執(zhí)行調(diào)整第一參數(shù)子集的次 數(shù)滿足第一預(yù)定條件，仍不能使得井儲向地層的過渡階段的井底計算壓力與井底實測壓力 的擬合精度達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)；
[0173]第三判斷子單元906,用于在所述第二判斷子單元903確定井儲向地層的過渡階段 的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層 過渡階段之后的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0174]第四處理子單元907,用于在所述第三判斷子單元906確定井儲向地層過渡階段之 后的井底計算壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第二參 數(shù)子集，所述第二參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、 絕對滲透率、油藏邊界和孔隙度；
[0175] 第五處理子單元908,用于在所述第四處理子單元907執(zhí)行調(diào)整第二參數(shù)子集的次 數(shù)滿足第二預(yù)定條件，仍不能使得井儲向地層過渡階段之后的井底計算壓力與井底實測壓 力的擬合精度達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)和表皮系數(shù)；
[0176]第四判斷子單元909,用于在所述第三判斷子單元906確定井儲向地層過渡階段之 后的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷整體的井 底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度；
[0177] 第六處理子單元910,用于在所述第四判斷子單元909確定整體的井底計算壓力與 井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，微調(diào)第三參數(shù)子集，所述第三參 數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊 界、孔隙度、表皮系數(shù)和井儲系數(shù)；
[0178] 第七處理子單元911，用于在所述第六處理子單元910執(zhí)行調(diào)整第三參數(shù)子集的次 數(shù)滿足第三預(yù)定條件，仍不能使得整體的井底計算壓力與井底實測壓力的擬合精度達到所 述預(yù)設(shè)精度的情況下，調(diào)整所述流動方程中的指數(shù)η。
[0179] 本發(fā)明圖6所示的第二處理單元90,根據(jù)井底計算壓力和井底實測壓力在各階段 的擬合精度，有針對性地調(diào)整特定的參數(shù)，能夠盡快完成整個參數(shù)解釋過程。
[0180] 另外，在本發(fā)明上述公開的低滲透油氣藏的參數(shù)解釋系統(tǒng)中，擬合單元70對井底 計算壓力與低滲透油氣藏的井底實測壓力進行擬合可以包括壓力史擬合、壓力降落或恢復 擬合、以及壓力導數(shù)擬合中的一項或多項，
[0181 ]作為一種優(yōu)選方案，擬合單元70具體用于：
[0182] 繪制井底計算壓力隨時間變化的第一曲線;繪制井底計算壓力的導數(shù)隨時間變化 的第二曲線;對第一曲線與第三曲線進行擬合，對第二曲線與第四曲線進行擬合，得到擬合 精度，其中第三曲線為井底實測壓力隨時間變化的曲線，第四曲線為井底實測壓力的導數(shù) 隨時間變化的曲線。
[0183] 最后，還需要說明的是，在本文中，諸如第一和第二等之類的關(guān)系術(shù)語僅僅用來將 一個實體或者操作與另一個實體或操作區(qū)分開來，而不一定要求或者暗示這些實體或操作 之間存在任何這種實際的關(guān)系或者順序。而且，術(shù)語"包括"、"包含"或者其任何其他變體意 在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、物品或者設(shè)備不僅包括那 些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、物品或者 設(shè)備所固有的要素。在沒有更多限制的情況下，由語句"包括一個……"限定的要素，并不排 除在包括所述要素的過程、方法、物品或者設(shè)備中還存在另外的相同要素。
[0184] 本說明書中各個實施例采用遞進的方式描述，每個實施例重點說明的都是與其他 實施例的不同之處，各個實施例之間相同相似部分互相參見即可。對于實施例公開的裝置 而言，由于其與實施例公開的方法相對應(yīng)，所以描述的比較簡單，相關(guān)之處參見方法部分說 明即可。
[0185] 對所公開的實施例的上述說明，使本領(lǐng)域?qū)I(yè)技術(shù)人員能夠?qū)崿F(xiàn)或使用本發(fā)明。 對這些實施例的多種修改對本領(lǐng)域的專業(yè)技術(shù)人員來說將是顯而易見的，本文中所定義的 一般原理可以在不脫離本發(fā)明的精神或范圍的情況下，在其它實施例中實現(xiàn)。因此，本發(fā)明 將不會被限制于本文所示的這些實施例，而是要符合與本文所公開的原理和新穎特點相一 致的最寬的范圍。
【主權(quán)項】
1. 一種低滲透油氣藏的參數(shù)解釋方法，其特征在于，包括： 根據(jù)所述低滲透油氣藏的流體特性，建立流體連續(xù)性方程； 設(shè)定流動方程，所述流動方程中的待擬合參數(shù)至少包括擬啟動壓力梯度、最大啟動壓 力梯度、最小啟動壓力梯度和絕對滲透率； 利用所述流體連續(xù)性方程和所述流動方程建立滲流方程； 確定所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流體壓縮系數(shù)； 設(shè)定所述低滲透油氣藏的參數(shù)集合，所述參數(shù)集合包括所述低滲透油氣藏的孔隙度、 地層的有效厚度、原始地層壓力、表皮系數(shù)、井儲系數(shù)、油藏邊界、W及所述流動方程包含的 待擬合參數(shù)； 利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和當前設(shè) 定的參數(shù)集合對所述滲流方程進行求解，獲得井底計算壓力； 對所述井底計算壓力與所述低滲透油氣藏的井底實測壓力進行擬合，得到擬合精度； 如果所述擬合精度達到預(yù)設(shè)精度，則確定當前設(shè)定的參數(shù)集合為所述低滲透率油氣藏 的實際參數(shù)； 如果所述擬合精度未達到所述預(yù)設(shè)精度，則調(diào)整所述低滲透油氣藏的參數(shù)集合，執(zhí)行 利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓縮系數(shù)和當前設(shè)定的 參數(shù)集合對所述滲流方程進行求解的步驟及后續(xù)步驟。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述流動方程為：其中，U為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓力梯 度，Amin為最小啟動壓力梯度，Apesud。為擬啟動壓力梯度，η為大于1的整數(shù)。3. 根據(jù)權(quán)利要求2所述的方法，其特征在于，所述調(diào)整所述低滲透油氣藏的參數(shù)集合， 包括： 步驟S91:判斷井儲階段的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù) 設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S92，若達到所述預(yù)設(shè)精度，則執(zhí)行步驟S93; 步驟S92:調(diào)整井儲系數(shù)； 步驟S93:判斷井儲向地層的過渡階段的井底計算壓力與井底實測壓力的擬合精度是 否達到所述預(yù)設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S94，若達到所述預(yù)設(shè)精度，則執(zhí) 行步驟S96; 步驟S94:調(diào)整第一參數(shù)子集，所述第一參數(shù)子集包括最小啟動壓力梯度、最大啟動壓 力梯度、擬啟動壓力梯度、絕對滲透率和表皮系數(shù)； 步驟S95:如果調(diào)整第一參數(shù)子集的次數(shù)滿足第一預(yù)定條件，仍不能使得井儲向地層的 過渡階段的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則微調(diào)井儲系 數(shù)； 步驟S96:判斷井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬合精度 是否達到所述預(yù)設(shè)精度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S97，若達到所述預(yù)設(shè)精度，貝U 轉(zhuǎn)至步驟S99; 步驟S97:調(diào)整第二參數(shù)子集，所述第二參數(shù)子集包括最小啟動壓力梯度、最大啟動壓 力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界和孔隙度； 步驟S98:如果調(diào)整第二參數(shù)子集的次數(shù)滿足第二預(yù)定條件，仍不能使得井儲向地層過 渡階段之后的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則微調(diào)井儲系 數(shù)和表皮系數(shù)； 步驟S99:判斷整體的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精 度，若未達到所述預(yù)設(shè)精度，則執(zhí)行步驟S910; 步驟S910:微調(diào)第Ξ參數(shù)子集，所述第Ξ參數(shù)子集包括最小啟動壓力梯度、最大啟動壓 力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界、孔隙度、表皮系數(shù)和井儲系數(shù)； 步驟S911:如果微調(diào)第Ξ參數(shù)子集的次數(shù)滿足第Ξ預(yù)定條件，仍不能使得整體的井底 計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度，則調(diào)整所述流動方程中的指數(shù)η。4. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述對所述井底計算壓力與所述低滲 透油氣藏的井底實測壓力進行擬合，包括： 繪制井底計算壓力隨時間變化的第一曲線； 繪制井底計算壓力的導數(shù)隨時間變化的第二曲線； 對所述第一曲線與第Ξ曲線進行擬合，對所述第二曲線與第四曲線進行擬合，得到擬 合精度，其中所述第Ξ曲線為井底實測壓力隨時間變化的曲線，所述第四曲線為所述井底 實測壓力的導數(shù)隨時間變化的曲線。5. -種低滲透油氣藏的參數(shù)解釋系統(tǒng)，其特征在于，包括： 流體連續(xù)性方程建立單元，用于根據(jù)所述低滲透油氣藏的流體特性，建立流體連續(xù)性 方程； 流動方程設(shè)定單元，用于設(shè)定流動方程，所述流動方程中的待擬合參數(shù)至少包括擬啟 動壓力梯度、最大啟動壓力梯度、最小啟動壓力梯度和絕對滲透率； 滲流方程建立單元，用于利用所述流體連續(xù)性方程和所述流動方程建立滲流方程； 參數(shù)確定單元，用于確定所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)和流 體壓縮系數(shù)； 參數(shù)設(shè)定單元，用于設(shè)定所述低滲透油氣藏的參數(shù)集合，所述參數(shù)集合包括所述低滲 透油氣藏的孔隙度、地層的有效厚度、原始地層壓力、表皮系數(shù)、井儲系數(shù)、油藏邊界、W及 所述流動方程包含的待擬合參數(shù)； 計算單元，用于利用所述低滲透油氣藏的井筒類型、流體性質(zhì)、巖石壓縮系數(shù)、流體壓 縮系數(shù)和當前設(shè)定的參數(shù)集合對所述滲流方程進行求解，獲得井底計算壓力； 擬合單元，用于對所述井底計算壓力與所述低滲透油氣藏的井底實測壓力進行擬合， 得到擬合精度； 第一處理單元，用于在所述擬合單元得到的擬合精度達到預(yù)設(shè)精度的情況下，確定當 前設(shè)定的參數(shù)集合為所述低滲透率油氣藏的實際參數(shù)； 第二處理單元，用于在所述擬合單元得到的擬合精度未達到所述預(yù)設(shè)精度的情況下， 調(diào)整所述低滲透油氣藏的參數(shù)集合，觸發(fā)所述計算單元執(zhí)行操作。6. 根據(jù)權(quán)利要求5所述的系統(tǒng)，其特征在于，所述流動方程設(shè)定單元設(shè)定的流動方程 為：其中，U為速度，K為絕對滲透率，μ為流體粘度，Vp為壓力梯度，Amax為最大啟動壓力梯 度,^min為最小啟動壓力梯度,^pesud。為擬啟動壓力梯度,η為大于1的整數(shù)。7. 根據(jù)權(quán)利要求5或6所述的系統(tǒng)，其特征在于，所述第二處理單元包括： 第一判斷子單元，用于判斷井儲階段的井底計算壓力與井底實測壓力的擬合精度是否 達到所述預(yù)設(shè)精度； 第一處理子單元，用于在所述第一判斷子單元確定井儲階段的井底計算壓力與井底實 測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整井儲系數(shù)； 第二判斷子單元，用于在所述第一判斷子單元確定井儲階段的井底計算壓力與井底實 測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層的過渡階段的井底計算壓 力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度； 第二處理子單元，用于在所述第二判斷子單元確定井儲向地層的過渡階段的井底計算 壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第一參數(shù)子集，所述 第一參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率和 表皮系數(shù)； 第Ξ處理子單元，用于在所述第二處理子單元執(zhí)行調(diào)整第一參數(shù)子集的次數(shù)滿足第一 預(yù)定條件，仍不能使得井儲向地層的過渡階段的井底計算壓力與井底實測壓力的擬合精度 達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)； 第Ξ判斷子單元，用于在所述第二判斷子單元確定井儲向地層的過渡階段的井底計算 壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷井儲向地層過渡階段之 后的井底計算壓力與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度； 第四處理子單元，用于在所述第Ξ判斷子單元確定井儲向地層過渡階段之后的井底計 算壓力與井底實測壓力的擬合精度未達到所述預(yù)設(shè)精度的情況下，調(diào)整第二參數(shù)子集，所 述第二參數(shù)子集包括最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透 率、油藏邊界和孔隙度； 第五處理子單元，用于在所述第四處理子單元執(zhí)行調(diào)整第二參數(shù)子集的次數(shù)滿足第二 預(yù)定條件，仍不能使得井儲向地層過渡階段之后的井底計算壓力與井底實測壓力的擬合精 度達到所述預(yù)設(shè)精度的情況下，微調(diào)井儲系數(shù)和表皮系數(shù)； 第四判斷子單元，用于在所述第Ξ判斷子單元確定井儲向地層過渡階段之后的井底計 算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度的情況下，判斷整體的井底計算壓力 與井底實測壓力的擬合精度是否達到所述預(yù)設(shè)精度； 第六處理子單元，用于在所述第四判斷子單元確定整體的井底計算壓力與井底實測壓 力的擬合精度未達到所述預(yù)設(shè)精度的情況下，微調(diào)第Ξ參數(shù)子集，所述第Ξ參數(shù)子集包括 最小啟動壓力梯度、最大啟動壓力梯度、擬啟動壓力梯度、絕對滲透率、油藏邊界、孔隙度、 表皮系數(shù)和井儲系數(shù)； 第屯處理子單元，用于在所述第六處理子單元執(zhí)行調(diào)整第Ξ參數(shù)子集的次數(shù)滿足第Ξ 預(yù)定條件，仍不能使得整體的井底計算壓力與井底實測壓力的擬合精度達到所述預(yù)設(shè)精度 的情況下，調(diào)整所述流動方程中的指數(shù)η。8.根據(jù)權(quán)利要求5或6所述的系統(tǒng)，其特征在于，所述擬合單元具體用于： 繪制井底計算壓力隨時間變化的第一曲線； 繪制井底計算壓力的導數(shù)隨時間變化的第二曲線； 對所述第一曲線與第Ξ曲線進行擬合，對所述第二曲線與第四曲線進行擬合，得到擬 合精度，其中所述第Ξ曲線為井底實測壓力隨時間變化的曲線，所述第四曲線為所述井底 實測壓力的導數(shù)隨時間變化的曲線。
【文檔編號】G06Q50/02GK106096300SQ201610457442
【公開日】2016年11月9日
【申請日】2016年6月20日
【發(fā)明人】李道倫, 查文舒, 鄭德溫, 盧德唐
【申請人】合肥工業(yè)大學