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

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

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