本發(fā)明涉及一種復(fù)雜斷塊油藏高含水后期剩余油二次富集的快速識別方法，屬于復(fù)雜斷塊油田開發(fā)領(lǐng)域。
背景技術(shù)：
：復(fù)雜斷塊油藏具有“小、散、貧、碎”的地質(zhì)特點，開發(fā)井網(wǎng)及儲量控制程度受儲層平面展布形狀的影響很大，多數(shù)斷塊不能形成的完善的注采井網(wǎng)。由此導(dǎo)致復(fù)雜斷塊油藏普遍存在含水上升快、采出程度低的問題，如何有效提高水驅(qū)開發(fā)后剩余油的二次富集程度和采出程度，是提高復(fù)雜斷塊油田水驅(qū)開發(fā)效果和提高油藏最終采收率的關(guān)鍵。常規(guī)油藏高含水后期剩余油富集的研究，大多沒有考慮注采井網(wǎng)不完善這一因素帶來的影響，富集規(guī)律和計算方法相對簡單。然而，在實際情況中，儲層的地質(zhì)特征和注采井網(wǎng)不完善雙重因素均會對復(fù)雜斷塊油藏產(chǎn)生影響，使得高含水后期儲層中依然會存在大量剩余油，且主要分布在非均質(zhì)性強的層段以及井網(wǎng)控制不到的區(qū)域。因此，常規(guī)油藏高含水后期剩余油富集的研究并不能滿足復(fù)雜斷塊油藏實際開發(fā)的需要，其在實際應(yīng)用過程中存在無法有效進(jìn)行剩余油二次富集的快速計算，不能準(zhǔn)確預(yù)測驅(qū)替前緣動態(tài)參數(shù)變化等一系列問題。綜上所述，提供一種新型的復(fù)雜斷塊油藏高含水后期剩余油二次富集的計算方法成為本領(lǐng)域亟待解決的技術(shù)問題。技術(shù)實現(xiàn)要素：為解決上述技術(shù)問題，本發(fā)明提供了一種復(fù)雜斷塊油藏高含水后期剩余油二次富集的快速識別方法，其能夠準(zhǔn)確預(yù)測驅(qū)替前緣參數(shù)的變化情況，使得復(fù)雜斷塊油藏高含水后期剩余油二次富集的快速計算得以實現(xiàn)。為達(dá)到上述目的，本發(fā)明提供了一種復(fù)雜斷塊油藏高含水后期剩余油二次富集的快速識別方法，該方法包括以下步驟：步驟S10，對目標(biāo)油藏進(jìn)行測量，以獲得目標(biāo)油藏的地質(zhì)參數(shù)和井網(wǎng)參數(shù)；步驟S11，根據(jù)目標(biāo)油藏的地質(zhì)參數(shù)和井網(wǎng)參數(shù)，建立目標(biāo)油藏的初步物理模型；步驟S12，基于目標(biāo)油藏的初步物理模型，根據(jù)流線流管法進(jìn)行擬合計算，以獲得采油井的動態(tài)擬合特征和剩余油二次富集前的飽和度場，并對目標(biāo)油藏的初步物理模型進(jìn)行修正，以得到修正后的物理模型；其中，所述采油井的動態(tài)擬合特征可以表示采油井的實際生產(chǎn)動態(tài)；步驟S13，基于修正后的物理模型，對剩余油富集過程中的各節(jié)點分別進(jìn)行垂直富集和水平富集的計算，以得到各節(jié)點的儲層飽和度和含水率；完成目標(biāo)油藏高含水后期剩余油二次富集的識別。與現(xiàn)有技術(shù)相比，本發(fā)明提供的技術(shù)方案同時對復(fù)雜斷塊油藏儲層的地質(zhì)特征和注采井網(wǎng)特征這兩個因素進(jìn)行了研究，并在此基礎(chǔ)上，通過建立目標(biāo)油藏的物理模型(該物理模型具有目標(biāo)油藏的典型地質(zhì)特征)，并結(jié)合本領(lǐng)域中的流線流管法理論和公式，最終得到的計算方法能夠有效應(yīng)用于復(fù)雜斷塊油藏，實現(xiàn)剩余油二次富集的快速計算，及軀替前緣動態(tài)參數(shù)變化的準(zhǔn)確預(yù)測。在本發(fā)明提供的技術(shù)方案中，所述流線流管法是指流線法和流管法，其中，流線法是對地下流體平面運移規(guī)律進(jìn)行數(shù)學(xué)描述的一種數(shù)學(xué)表征方法；流管法是對地下流體所處的空間結(jié)構(gòu)和空間內(nèi)各節(jié)點運移關(guān)系進(jìn)行數(shù)學(xué)描述的一種物理模型構(gòu)建方法；這兩種方法都是本領(lǐng)域已知的計算方法，是本領(lǐng)域的專業(yè)術(shù)語。在上述計算方法中，優(yōu)選地，在步驟S12中，所述采油井的動態(tài)擬合特征包括以下參數(shù)：產(chǎn)油量、產(chǎn)水量和含水率。在上述計算方法中，優(yōu)選地，在步驟S11中，根據(jù)目標(biāo)油藏的地質(zhì)參數(shù)和井網(wǎng)參數(shù)，建立目標(biāo)油藏的物理模型包括以下步驟：根據(jù)獲得的井網(wǎng)參數(shù)，設(shè)置初步物理模型的底面類型；根據(jù)獲得的地質(zhì)參數(shù)，設(shè)置初步物理模型的地質(zhì)參數(shù)。在上述計算方法中，優(yōu)選地，根據(jù)獲得的井網(wǎng)參數(shù)，設(shè)置初步物理模型的底面類型包括以下步驟：根據(jù)獲得的井網(wǎng)參數(shù)，確定目標(biāo)油藏的井網(wǎng)類型、油水井?dāng)?shù)比、井距和排距；根據(jù)所述井網(wǎng)類型和油水井?dāng)?shù)比，設(shè)置初步物理模型的底面類型(設(shè)置流程如圖2所示)：所述井網(wǎng)類型為交錯時，設(shè)置初步物理模型的底面類型為正三角形；所述井網(wǎng)類型為正對，且油水井?dāng)?shù)比＝1時，設(shè)置初步物理模型的底面類型為矩形；所述井網(wǎng)類型為正對，且油水井?dāng)?shù)比≠1時，設(shè)置初步物理模型的底面類型為正三角形；初步物理模型的底面類型為矩形時，設(shè)置矩形的長等于井距，矩形的寬等于排距；初步物理模型的底面類型為正三角形時，設(shè)置正三角形的邊長的等于井距，正三角形的高等于排距。在本發(fā)明提供的技術(shù)方案中，初步物理模型是一個具有三維空間結(jié)構(gòu)的模型，該三維空間結(jié)構(gòu)包括水平X方向、水平Y(jié)方向和垂直Z方向(水平X方向與水平Y(jié)方向相互垂直，垂直Z方向分別垂直于X方向和Y方向)，初步物理模型的底面在所述水平X方向和水平Y(jié)方向共同構(gòu)成的水平面上，按上述方法確定好底面類型后，將底面沿著垂直Z方向平移一段距離(該段距離即為初步物理模型的厚度，其等于實際儲層的厚度)后得到的三維空間結(jié)構(gòu)即為初步物理模型的空間結(jié)構(gòu)，例如按照上述方法確定底面類型為矩形，將矩形沿著豎直方向平移一段距離后得到立方體結(jié)構(gòu)，該立方體結(jié)構(gòu)即為初步物理模型的空間結(jié)構(gòu)。在該空間結(jié)構(gòu)中，將每個單油層按照水平X、水平Y(jié)、垂直Z劃分三個方向，在水平面上按照流體在地下滲流時的主流線方向，將平面劃分成多條流管(如圖8所示)；每條流管的長度根據(jù)油藏形狀的不同而不同，數(shù)值上等于單根流管所在平面對應(yīng)位置的油藏長度。對于規(guī)則的矩形油藏，同一平面上的每條流管的長度相同；對于三角形油藏，同一平面上的每條流管的長度不同。每條流管的寬度(所述寬度為流管的直徑)和長度(所述長度為主流線方向上流管的長度)在數(shù)值上相同，劃分尺度根據(jù)油藏大小靈活決定，取值范圍為1-20m。平面流管劃分完成后，在與主流線垂直的方向，依據(jù)與平面設(shè)置相同的流管的寬度和長度繼續(xù)劃分，三個方向上的流管會在空間中相交形成網(wǎng)格，相交點就是計算節(jié)點(如圖9所示，圖中圓點表示節(jié)點，圓管表示流管)；主流線方向上的每條流管都對應(yīng)一個產(chǎn)出端。由于地層是非均質(zhì)的，所以流體在地下滲流過程中通常具有各向異性，即流動方向不是單一的，可能四面八方都有流動，但一定會存在一個主要流動方向，這與開采有關(guān)。因為不論是采油井還是注水井，在井點附近都會造成地層壓力的變化(采油井點附近的地層壓力會降低，注水井點附近的地層壓力會升高)，從而在開采區(qū)域內(nèi)形成一個壓力差，該壓力差就是驅(qū)動地下流體流動的主要動力來源，這個方向就是本發(fā)明所述的主流線方向。在上述計算方法中，優(yōu)選地，所述步驟S12包括以下步驟：步驟S121，對初步物理模型的儲層特征參數(shù)進(jìn)行預(yù)設(shè)，得到預(yù)設(shè)的儲層特征；利用采油井的實際動態(tài)歷史特征對預(yù)設(shè)的儲層特征參數(shù)進(jìn)行擬合，以確定最終的采油井的動態(tài)擬合特征和初步物理模型的儲層特征；其中，所述采油井的實際動態(tài)歷史特征包括產(chǎn)油量、產(chǎn)水量和含水率這些參數(shù)；步驟S122，根據(jù)流線流管法，初步擬合得到剩余油二次富集前的含油飽和度場；利用飽和度修正關(guān)系式對初步擬合得到的剩余油二次富集前的含油飽和度場進(jìn)行擬合修正，以確定最終的剩余油二次富集前的含油飽和度場；步驟S123，對目標(biāo)油藏的初步物理模型進(jìn)行修正，得到修正后的物理模型。在上述計算方法中，優(yōu)選地，在步驟S121中，所述儲層的特征參數(shù)包括初始含油飽和度、儲層厚度、儲層傾角、注采速度和平均滲透率。在上述計算方法中，優(yōu)選地，在步驟S121中，利用采油井的實際動態(tài)歷史特征對預(yù)設(shè)的儲層特征參數(shù)進(jìn)行擬合，以確定最終的采油井的動態(tài)擬合特征和初步物理模型的儲層特征，包括以下步驟：基于目標(biāo)油藏的初步物理模型，根據(jù)流線流管法進(jìn)行擬合計算，獲得初步擬合得到的采油井的動態(tài)擬合特征；初步擬合得到的采油井的動態(tài)擬合特征與目標(biāo)油藏的采油井的實際動態(tài)歷史特征相符合時，則判斷初步擬合得到的采油井的動態(tài)擬合特征是符合要求的；初步擬合得到的采油井的動態(tài)擬合特征與目標(biāo)油藏的采油井的實際動態(tài)歷史特征不相符合時，則判斷初步擬合得到的采油井的動態(tài)擬合特征是不符合要求的；此時，需要對初步物理模型的儲層特征參數(shù)進(jìn)行修改，并根據(jù)流線流管法對采油井的動態(tài)擬合特征進(jìn)行重新擬合，直至采油井的動態(tài)擬合特征與目標(biāo)油藏上采油井的實際動態(tài)歷史特征相符合；初步擬合得到的采油井的動態(tài)擬合特征與目標(biāo)油藏的采油井的實際動態(tài)歷史特征相符合時，確定其對應(yīng)的儲層特征即為初步物理模型的儲層特征。在上述計算方法中，優(yōu)選的，在步驟S122中，根據(jù)流線流管法，初步擬合得到剩余油二次富集前的含油飽和度場包括以下步驟：計算主流線方向上各個節(jié)點的含水率和含水上升率；計算主流線方向上的含水前緣飽和度及含水前緣位置；計算主流線方向上產(chǎn)出端的總含水率。在上述計算方法中，優(yōu)選地，在步驟S122中，利用飽和度修正關(guān)系式對初步擬合得到的剩余油二次富集前的含油飽和度場進(jìn)行擬合修正，包括以下步驟：設(shè)置飽和度修正關(guān)系式，對初步擬合得到的剩余油二次富集前的含油飽和度場進(jìn)行修正，得到修正后的含油飽和度場；所述修正后的含油飽和度場與目標(biāo)油藏的實際含油飽和度場相符合時，則判斷修正后的含油飽和度場是符合要求的；所述修正后的含油飽和度場與目標(biāo)油藏的實際飽和度場不相符合時，則判斷修正后的含油飽和度場是不符合要求的；此時，重新設(shè)置飽和度修正關(guān)系式，對剩余油二次富集前的含油飽和度場進(jìn)行重新修正，直至修正后的含油飽和度場與目標(biāo)油藏的實際含油飽和度場相符合；修正后的含油飽和度場與目標(biāo)油藏的實際含油飽和度場相符合時，確定其為最終的剩余油二次富集前的含油飽和度場。在本發(fā)明提供的技術(shù)方案中，地下實際的飽和度場是未知的，通過單井測試獲得的信息，僅能代表地下一個點，即該井所在位置當(dāng)時的飽和度大小。對于整體油藏而言，通過若干個井點獲得的飽和度數(shù)據(jù)，對于整個油藏范圍而言仍然是非常有限的，因此需要通過這些已知數(shù)據(jù)對油藏其它位置的飽和度進(jìn)行預(yù)測，預(yù)測數(shù)據(jù)通常是根據(jù)相鄰已知點的數(shù)據(jù)，通過插值方法得到相鄰兩點間的未知數(shù)據(jù)，但這樣插值得到的預(yù)測結(jié)果，并不能代表真實的油藏飽和度，因此需要檢驗和修正。而檢驗的方法就是通過不斷修正插值方法預(yù)測得到的飽和度，來擬合單井動態(tài)(包括產(chǎn)油、產(chǎn)水、含水率)，因為單井動態(tài)是實際能測得的數(shù)據(jù)。如果單井?dāng)M合結(jié)果與實際單井動態(tài)是吻合的，那么說明預(yù)測飽和度或者修正后的飽和度，與實際飽和度場吻合；否則，需要繼續(xù)修正。在上述計算方法中，優(yōu)選地，在步驟S11中，該計算方法還包括按照目標(biāo)油藏的實際地質(zhì)結(jié)構(gòu)特征，設(shè)置初步物理模型的沉積韻律、地質(zhì)傾角和滲透率各向異性；其中：沉積韻律是指按顆粒大小、比重的順序先后分層沉積而成成巖的規(guī)律，分為均質(zhì)韻律、正韻律、反韻律三種；均質(zhì)韻律砂層：水動力條件相對穩(wěn)定，層內(nèi)巖性，物性均質(zhì)；正韻律砂層：下部粒度粗，上部粒度細(xì)，反映了沉積環(huán)境水動力條件由下到上變?nèi)?；反韻律砂層：下部粒度?xì)，上部粒度粗，沉積環(huán)境水動力條件由下到上變強；地層傾角是指油水井之間油層走向與水平面之間的夾角；滲透率各向異性是指油層在不同方向上都存在滲透率變化的差異性。在上述計算方法中，優(yōu)選地，根據(jù)獲得的地質(zhì)參數(shù)，設(shè)置初步物理模型的地質(zhì)參數(shù)，主要是指按照獲得的目標(biāo)油藏的實際地質(zhì)參數(shù)，設(shè)置初步物理模型的平均滲透率、初始含油飽和度、孔隙度、束縛水飽和度、殘余油飽和度、滲透率級差、地下原油粘度、儲層厚度、儲層劃分層數(shù)、注采速度等油藏測量中需要獲取的常規(guī)參數(shù)；其中：平均滲透率是指在一定壓差下，巖石允許流體通過的能力大小；初始含油飽和度是油層有效孔隙中含油體積和巖石有效孔隙體積之比，以百分?jǐn)?shù)表示；孔隙度是指巖樣中所有孔隙空間體積之和與該巖樣體積的比值；束縛水飽和度是指由于巖石表面潤濕性導(dǎo)致殘留在巖石孔隙中的最小水體所占儲集層孔隙體積的百分比；殘余油飽和度是指殘余油在巖石孔隙中所占體積的百分?jǐn)?shù)；滲透率級差是指最大滲透率與最小滲透率的比值；地下原油粘度是指在地層條件下原油內(nèi)部某一部分相對于另一部分流動時摩擦阻力的度量。在上述計算方法中，優(yōu)選地，在步驟S11結(jié)束后，步驟S12開始前，該計算方法包括基于目標(biāo)油藏的初步物理模型，對初步物理模型的縱向劃分層數(shù)和縱向非均質(zhì)性進(jìn)行預(yù)設(shè)的步驟，其包括以下過程：初步物理模型的單層厚度＝初步物理模型的總厚度/預(yù)設(shè)的縱向劃分層數(shù)；初步物理模型的初始非均質(zhì)性＝預(yù)設(shè)的縱向非均質(zhì)性。在上述計算方法中，優(yōu)選地，在步驟S122中，基于目標(biāo)油藏的初步物理模型，根據(jù)流線流管法，初步擬合得到剩余油二次富集前的飽和度場包括以下步驟：①計算主流線方向上各個節(jié)點的含水率和含水上升率；此處的含水率指的是主流線方向上利用貝克萊-列維爾特公式(簡稱B-L公式)計算得到的含水率，物理模型中主流線上的各個節(jié)點通過該公式計算后都能獲得一個對應(yīng)的數(shù)據(jù)，因此得到的是一系列數(shù)據(jù)點；②計算主流線方向上的含水前緣飽和度及含水前緣位置；③計算主流線方向上產(chǎn)出端的總含水率；④對主流線方向上產(chǎn)出端的總含水率進(jìn)行判斷，根據(jù)采油井的實際含水率對儲層特征參數(shù)進(jìn)行修正，使采油井的擬合含水率與實際含水率相符；⑤基于步驟①至步驟④的計算方法，得到采油井的擬合含水率達(dá)到剩余油二次富集要求時對應(yīng)的含油飽和度場，該含油飽和度場為預(yù)測值；⑥對邊界至主流線進(jìn)行含油飽和度插值，以修正步驟⑤得到的含油飽和度場，從而得到剩余油二次富集前的飽和度場。在上述計算方法中，優(yōu)選地，在步驟①中：所述主流線方向上各個節(jié)點的含水率的計算公式如式1所示，式1可以用于計算物理模型中任意一個點對應(yīng)的含水率在式1中，fw為含水率，小數(shù)；krw、kro為水相、油相的相對滲透率，無因次；μw、μo為水相、油相的粘度，mPa·s；更優(yōu)選地，所述含水上升率的計算公式如式2所示在式2中，f'w為含水上升率(含水率的倒數(shù))，小數(shù)；Sw為含水飽和度，小數(shù)；i為第i個節(jié)點位置處，i-1為i節(jié)點的前一個節(jié)點位置處。在上述計算方法中，優(yōu)選的，在步驟②中，計算含水前緣飽和度及含水前緣位置包括：根據(jù)Sw—fw關(guān)系曲線，通過Swi點與每個Sw—fw關(guān)系曲線節(jié)點連接成線，求該線的導(dǎo)數(shù)，導(dǎo)數(shù)值最大的節(jié)點，則為前緣含水飽和度Swf；再根據(jù)Sw—f'w關(guān)系，得到f'w(Swf)值；最后，求出含水前緣位置xf。在上述計算方法中，優(yōu)選地，在步驟②中，所述主流線上的含水前緣位置的計算公式如式3所示其中，在式3中，xf為含水前緣位置，m；x0為含水初始位置，m；f'w(Swf)為含水前緣飽和度對應(yīng)的含水上升率，小數(shù)；Φ為孔隙度，小數(shù)；A為截面積，m2；Q為流量，m3；t為驅(qū)替時間，天；fw(Swf)為含水前緣飽和度對應(yīng)的含水率，小數(shù)；Swf為含水前緣飽和度，小數(shù)；Swc為束縛水飽和度，小數(shù)。在上述計算方法中，優(yōu)選地，在步驟③中，計算主流線方向上產(chǎn)出端的總含水率包括：根據(jù)各層的注水量和每條流管的滲透率的比值，對每條流管的流量進(jìn)行劈分，計算各層的見水時間和各層產(chǎn)出端的含水率；根據(jù)每條流管的流量對各層產(chǎn)出端的含水率進(jìn)行加權(quán)平均，以得到主流線上產(chǎn)出端的總含水率。在上述計算方法中，優(yōu)選地，在步驟③中，優(yōu)選地，所述各層的見水時間的計算公式如式4所示在式4中，xf為含水前緣位置，m；x0為含水初始位置，m；f'w(Swf)為含水前緣飽和度對應(yīng)的含水上升率，小數(shù)；A為截面積，m2；Q為流量，m3。在上述計算方法中，優(yōu)選地，在步驟③中，各層產(chǎn)出端的含水率的計算包括：在計算的時候先判斷各流管的見水時間，然后采用流量劈分法，根據(jù)式5所示的公式，計算各層產(chǎn)出端的含水飽和度對應(yīng)的含水上升率，然后根據(jù)式2所示的公式，計算各層產(chǎn)出端的含水率fw；在式5中，t為生產(chǎn)時間，天；T為見水時間，天；L為油水井距，m；f'w(SwL)為各層產(chǎn)出端的含水飽和度對應(yīng)的含水上升率，小數(shù)。在上述方法中，優(yōu)選地，在步驟④中，對主流線方向上產(chǎn)出端的總含水率進(jìn)行判斷包括以下步驟：從初始時刻開始，將計算得到的主流線方向上產(chǎn)出端的總含水率與采油井的實際含水率進(jìn)行數(shù)值對比；當(dāng)兩者數(shù)值不符時，對儲層特征參數(shù)進(jìn)行修正，若修正后的模型得到的計算含水率與實際含水率的數(shù)值差異較大，則需繼續(xù)修正儲層特征參數(shù)；若修正后的模型得到的計算含水率與實際含水率的數(shù)值一致或差異很小，則該時刻點對應(yīng)的含水率的擬合過程完成，開始判斷下一個時刻點對應(yīng)的含水率；直至實際含水率測量的終止時刻，則判斷過程結(jié)束，含水率擬合完成。其中，對于稀油油藏而言，數(shù)值差異在-5％-5％以內(nèi)時，認(rèn)為差異較小，否則認(rèn)為差異較大；對于稠油油藏而言，數(shù)值差異在-15％-15％以內(nèi)時，認(rèn)為差異較小，否則認(rèn)為差異較大。在上述計算方法中，優(yōu)選地，在步驟S13中，所述垂直富集的計算過程包括：①計算每條流管各節(jié)點的垂向高度差，然后分別按照從上往下的方式和從下往上的方式對每條流管各節(jié)點的垂向高度差進(jìn)行累計求和，以獲得垂向上每條流管各節(jié)點的總垂向高度差(該總垂向高度差包括垂向上從上往下每條流管各節(jié)點的總垂向高度差，和垂向上從下往上每條流管給階段的總垂向高度差)；②計算垂向上每條流管各節(jié)點完全平衡時的平均含水飽和度；其中，所述完全平衡指的是剩余油經(jīng)過富集后再不能富集的狀態(tài)，此時垂直方向與水平方向的含油飽和度不再發(fā)生變化；將每條流管所在層位的厚度與步驟①獲得的垂向上每條流管各節(jié)點的總垂向高度差進(jìn)行大小對比，從而得到垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最大飽和度平均值和最小飽和度平均值；③將每條流管各節(jié)點完全平衡時的平均含水飽和度與飽和度剖面上各節(jié)點的最大含水飽和度平均值和最小含水飽和度平均值進(jìn)行大小對比，以獲得剖面平衡態(tài)類型A、剖面平衡態(tài)次級類型B、剖面平衡態(tài)類型的節(jié)點間比例系數(shù)a、剖面平衡態(tài)次級類型的節(jié)點間比例系數(shù)b、剖面平衡態(tài)類型A的平衡位置和剖面平衡態(tài)次級類型B的平衡位置；④計算次生油水界面的剖面含水飽和度和剖面平衡態(tài)等效毛管力；⑤計算垂向上每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率。在上述計算方法中，所述剖面平衡態(tài)類型A是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度大于最大含水飽和度的平均值的節(jié)點個數(shù)；所述剖面平衡態(tài)次級類型B是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度小于最小含水飽的平均值的節(jié)點個數(shù)；所述節(jié)點間比例系數(shù)包含剖面平衡態(tài)類型節(jié)點間比例系數(shù)a和剖面平衡態(tài)次級類型節(jié)點間比例系數(shù)b；其中，剖面平衡態(tài)類型節(jié)點間比例系數(shù)a是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度小于最大含水飽和度的平均值的節(jié)點個數(shù)占該條流管節(jié)點總數(shù)的比值；剖面平衡態(tài)次級類型節(jié)點間比例系數(shù)b是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度大于最大含水飽和度的平均值的節(jié)點個數(shù)占該條流管節(jié)點總數(shù)的比值；所述平衡位置是指在完全平衡條件下平均含水飽和度值所處的油層剖面的位置，其包括剖面平衡態(tài)類型A的平衡位置和剖面平衡態(tài)次級類型B的平衡位置；所述次生油水界面的剖面含水飽和度是指由于重力影響下產(chǎn)生次生油水界面下計算出剖面各節(jié)點的含水飽和度值；所述剖面平衡態(tài)等效毛管力是指根據(jù)次生油水界面的剖面含水飽和度計算出考慮重力影響的剖面各節(jié)點處的毛管力值。在上述計算方法中，優(yōu)選地，在步驟①中，所述每條流管各節(jié)點的垂向高度差的計算公式如式6所示在式6中，Δh為由毛管力引起的高度差，m；Pc為每個節(jié)點處的毛管力，atm；γw為水的重度，×104N/m3；γo為油的重度，×104N/m3；i為下標(biāo)表示第i個節(jié)點位置處，i+1表示i節(jié)點的下一節(jié)點位置處，例如Pc(i)表示第i個節(jié)點位置處的毛管力，Pc(i+1)表示第i+1個節(jié)點位置處的毛管力。優(yōu)選地，在步驟①中，按照從上往下的方式對每條流管各節(jié)點的垂向高度差進(jìn)行累計求和，得到垂向上各節(jié)點的總垂向高度差時，計算公式如式7所示在式7中，hud(i)為垂向上從上往下每條流管各節(jié)點的總垂向高度差；i為第i個節(jié)點位置處；j為第j個節(jié)點位置處；優(yōu)選地，在步驟①中，按照從下往上的方式對每條流管各節(jié)點的垂向高度差進(jìn)行累計求和，得到垂向上各節(jié)點的總垂向高度差時，計算公式如式8所示在式8中，hdu(i)為垂向上從下往上每條流管各節(jié)點的總垂向高度差；i為第i個節(jié)點位置處；j為第j個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)。在上述計算方法中，優(yōu)選地，在步驟②中，垂向上每條流管各節(jié)點完全平衡時的平均含水飽和度的計算公式如9所示在式9中，為垂向上每條流管各節(jié)點完全平衡時的平均含水飽和度；Swz為垂向上每條流管完全平衡時各節(jié)點的含水飽和度；i為第i個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)；優(yōu)選地，在步驟②中，垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最大含水飽和度的平均值的計算公式如式10所示在式10中，為垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最大含水飽和度的平均值；Sw為各節(jié)點含水飽和度；hud為垂向上從上往下每條流管各節(jié)點的總垂向高度差；h為小層厚度；Δh為垂向上每條流管各節(jié)點的垂向高度差；i為第i個節(jié)點位置處；j為第j個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)；優(yōu)選地，在步驟②中，垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最小含水飽和度的平均值的計算公式如式11所示在式11中，為垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最小含水飽和度的平均值；Sw為各節(jié)點含水飽和度；hdu為垂向上從下往上每條流管各節(jié)點的總垂向高度差；h為小層厚度；Δh為垂向上每條流管各節(jié)點的垂向高度差；i為第i個節(jié)點位置處；j為第j個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)。在上述計算方法中，優(yōu)選地，在步驟③中，所述剖面平衡態(tài)類型的節(jié)點間比例系數(shù)的計算公式如式12所示在式12中，a為剖面平衡態(tài)類型節(jié)點間比例系數(shù)；為垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最大含水飽和度的平均值；為垂向上每條流管各節(jié)點的平均含水飽和度；A為剖面平衡態(tài)類型；i為第i個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)。在上述計算方法中，優(yōu)選地，在步驟③中，所述剖面平衡態(tài)次級類型的節(jié)點間比例系數(shù)的計算公式如式13所示在式13中，b為剖面平衡態(tài)次級類型節(jié)點間比例系數(shù)；為垂向上每條流管各節(jié)點完全平衡時，飽和度剖面上各節(jié)點的最小含水飽和度的平均值；為垂向上每條流管各節(jié)點的平均含水飽和度；B為剖面平衡態(tài)次級類型；i為第i個節(jié)點位置處；n為垂向上每條流管節(jié)點總數(shù)。在上述計算方法中，優(yōu)選地，在步驟③中，剖面平衡態(tài)類型的平衡位置的計算公式如式14所示在式14中，xa為剖面平衡態(tài)類型平衡位置；h為小層厚度，m；hud為垂向上從上往下每條流管各節(jié)點的總垂向高度差，m；a為剖面平衡態(tài)類型節(jié)點間比例系數(shù)；i為第i個節(jié)點位置處；A為剖面平衡態(tài)類型。在上述計算方法中，優(yōu)選地，在步驟③中，剖面平衡態(tài)次級類型的平衡位置的計算公式如式15所示在式15中，xb為剖面平衡態(tài)次級類型平衡位置；h為小層厚度，m；hdu為垂向上從下往上每條流管各節(jié)點的總垂向高度差，m；b為剖面平衡態(tài)次級類型節(jié)點間比例系數(shù)；i為第i個節(jié)點位置處；B為剖面平衡態(tài)次級類型。在上述計算方法中，優(yōu)選地，在步驟④中，所述次生油水界面的剖面含水飽和度的計算公式如式16所示在式16中，Swpm為次生WOC界面位置每個節(jié)點的剖面含水飽和度，小數(shù)；Sw為相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Δh為由毛管力引起的油水高度差，m；Hpm為次生WOC界面位置每個節(jié)點的剖面高度，m；i為第i個節(jié)點位置處。在上述計算方法中，優(yōu)選地，在步驟④中，所述次生油水界面的剖面平衡態(tài)等效毛管力的計算公式如式17所示在式17中，Swpm為次生WOC界面位置每個節(jié)點的剖面含水飽和度，小數(shù)；Sw為相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Pcpm為次生WOC界面位置每個節(jié)點的剖面毛管力，atm；Pc為毛管壓力曲線中每個節(jié)點對應(yīng)的毛管力，atm；i為下標(biāo)表示第i個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處。在上述計算方法中，優(yōu)選地，在步驟⑤中，垂向上每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率的計算過程包括：步驟1，計算在富集初始時刻，垂向上每條流管各節(jié)點對應(yīng)的油相、水相的相對滲透率，計算公式如式18所示在式18中，Kr(t0)為富集初始時刻油水每個節(jié)點相對滲透率，無因次；Kr為油水相滲中每個節(jié)點對應(yīng)的相對滲透率，無因次；Sw為油水相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Sw(t0)為富集初始時刻每個節(jié)點的含水飽和度，小數(shù)；i為下標(biāo)表示第i個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處；步驟2，計算在富集過程中，垂向上每條流管各節(jié)點對應(yīng)的含油飽和度，計算公式如式19所示在式19中，So為每個節(jié)點對應(yīng)的含油飽和度，小數(shù)；KZ為Z方向上每個節(jié)點對應(yīng)的滲透率，mD；Kro為每個節(jié)點的油相對滲透率，無因次；Krw為每個節(jié)點的水相對滲透率，無因次；μw為水的粘度，mPa·s；μo為油的粘度，mPa·s；Pc為每個節(jié)點的毛管力，atm；Pc∞為次生WOC界面位置每個節(jié)點的剖面毛管力，atm；φ為孔隙度，小數(shù)；t為富集時間，天；n為劃分的時間段次數(shù)；i為第i個節(jié)點位置處；步驟3，計算在富集完成時刻，垂向上每條流管各節(jié)點對應(yīng)的油相、水相的相對滲透率，計算公式如式20所示在式20中，Kr(tj)為計算富集完成時刻每個節(jié)點的相對滲透率，無因次；Kr為油水相滲中每個節(jié)點對應(yīng)的相對滲透率，無因次；Sw為油水相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Sw(tj)為計算富集完成時刻每個節(jié)點的含水飽和度，小數(shù)；i為下標(biāo)表示第i個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處。在上述計算方法中，優(yōu)選地，在步驟S13中，所述水平富集的計算過程包括：步驟1，計算在富集初始時刻，平面上每條流管各節(jié)點對應(yīng)的毛管力，計算公式如式21所示在式21中，Pc(t0)為富集初始時刻每個節(jié)點的毛管力，atm；Pc為毛管壓力曲線中每個節(jié)點對應(yīng)的毛管力，atm；Sw為油水相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Sw(t0)為富集初始時刻每個節(jié)點的含水飽和度，小數(shù)；i為下標(biāo)表示第i個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處；步驟2，計算在富集過程中，平面上每條流管各節(jié)點對應(yīng)的含油飽和度，計算公式如式22所示在式22中，So為每個節(jié)點對應(yīng)的含油飽和度，小數(shù)；KZ為X方向上每個節(jié)點對應(yīng)的滲透率，mD；Kro為每個節(jié)點的油相對滲透率，無因次；Krw為每個節(jié)點的水相對滲透率，無因次；μw為水的粘度，mPa·s；μo為油的粘度，mPa·s；Pc為每個節(jié)點的毛管力，atm；φ為孔隙度，小數(shù)；t為時間，天；n為劃分的時間段次數(shù)；i為下標(biāo)表示第i個節(jié)點位置處，i+1表示i節(jié)點的后一個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處；步驟3，計算在富集完成時刻，平面上每條流管各節(jié)點對應(yīng)的毛管力，計算公式如式23所示在式23中，Pc(tj)為富集完成時刻每個節(jié)點毛管力，atm；Pc為毛管壓力曲線中每個節(jié)點對應(yīng)的毛管力，atm；Sw為油水相滲中每個節(jié)點對應(yīng)的含水飽和度，小數(shù)；Sw(tj)為富集完成時刻每個節(jié)點的含水飽和度，小數(shù)；i為下標(biāo)表示第i個節(jié)點位置處，i-1表示i節(jié)點的前一個節(jié)點位置處；步驟4，根據(jù)計算得到的富集過程中平面上每條流管各節(jié)點對應(yīng)的含油飽和度，以及公式20，得到平面上每條流管各節(jié)點的油相、水相的相對滲透率。在一個具體實施方式中，步驟S12包含以下過程(如圖3和圖4所示)：(1)對初步物理模型的特征參數(shù)進(jìn)行預(yù)設(shè)，得到預(yù)設(shè)后的初步物理模型；其中，設(shè)置初步物理模型的單層厚度＝初步物理模型的總厚度/預(yù)設(shè)縱向劃分層數(shù)；設(shè)置初步物理模型的初始非均質(zhì)性＝預(yù)設(shè)縱向非均質(zhì)性；(2)基于預(yù)設(shè)后的初步物理模型，根據(jù)流線流管法計算得到采油井的擬合含水率；(3)對計算得到的采油井的擬合含水率進(jìn)行判斷：所述擬合含水率與目標(biāo)油藏上采油井的實際含水率相符合，判斷所述擬合含水率是符合要求的；所述擬合含水率與目標(biāo)油藏上采油井的實際含水率不相符合，判斷所述擬合含水率是不符合要求的；此時，需要對初步物理模型的縱向劃分層數(shù)和縱向非均質(zhì)性進(jìn)行修改，并重復(fù)步驟(2)-(3)的操作過程，直至所述擬合含水率與目標(biāo)油藏上采油井的實際含水率相符合；(4)基于預(yù)設(shè)后的初步物理模型，設(shè)置目標(biāo)含水率，根據(jù)流線流管法計算采油井的含水率達(dá)到目標(biāo)含水率時的含油飽和度場；其中，所述目標(biāo)含水率即為采油井的含水率達(dá)到關(guān)井富集時機時對應(yīng)的含水率；(5)設(shè)置飽和度修正關(guān)系式，并利用其對步驟(4)計算得到的含油飽和度場進(jìn)行修正，得到修正后的含油飽和度場；(6)對修正后的含油飽和度場進(jìn)行判斷：所述修正后的含油飽和度場與目標(biāo)油藏的實際飽和度場相符合，判斷修正后的含油飽和度場是符合要求的；所述修正后的含油飽和度場與目標(biāo)油藏的實際飽和度場不相符合，判斷修正后的含油飽和度場是不符合要求的；此時，重新設(shè)置飽和度修正關(guān)系式，重復(fù)步驟(5)至步驟(6)的操作過程，直至修正后的含油飽和度場與目標(biāo)油藏的實際飽和度場相符合；(7)對目標(biāo)油藏的初步物理模型進(jìn)行修正，得到修正后的物理模型。在一個具體實施方式中，步驟S13包含以下過程：(1)計算每條流管各節(jié)點的垂向高度差，分別按照從上往下的方式和從下往上的方式進(jìn)行累計求和，得到垂向各節(jié)點位置對應(yīng)的總垂向高度差；(2)計算垂向上每條流管各節(jié)點完全平衡時的平均含水飽和度，根據(jù)每條流管所在層位的厚度與獲得的垂向上各節(jié)點位置對應(yīng)的總垂向高度差進(jìn)行大小對比，得到所述飽和度剖面上各節(jié)點的最大含水飽和度平均值和最小含水飽和度平均值；(3)將垂向上的平均含水飽和度，與飽和度剖面上各節(jié)點的最大含水飽和度平均值和最小含水飽和度平均值進(jìn)行大小對比，得到剖面平衡態(tài)類型、剖面平衡態(tài)次級類型、剖面平衡態(tài)類型的節(jié)點間比例系數(shù)、剖面平衡態(tài)次級類型的節(jié)點間比例系數(shù)、剖面平衡態(tài)類型的平衡位置和剖面平衡態(tài)次級類型的平衡位置；(4)計算次生油水界面的剖面含水飽和度和剖面平衡態(tài)等效毛管力；(5)預(yù)設(shè)富集時間和富集次數(shù)，在每次循環(huán)計算中，分別進(jìn)行垂向和水平向富集計算，計算出每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率(如圖5所示)；(6)每循環(huán)完一次垂向和水平向富集計算，進(jìn)行判斷是否達(dá)到預(yù)設(shè)的富集次數(shù)，若未達(dá)到，再進(jìn)行(5)步驟計算，若達(dá)到富集次數(shù)，則輸出保存最后一次富集計算完的每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率數(shù)據(jù)；(7)根據(jù)(6)步所得的數(shù)據(jù)，計算出采油井端含水率，即是采油井開井含水率，輸出采油井開井含水率，計算結(jié)束。本發(fā)明的有益效果：現(xiàn)有方法不適用與復(fù)雜斷塊油藏的地質(zhì)和井網(wǎng)特點，無法有效進(jìn)行剩余油二次富集的快速計算，本發(fā)明提供的技術(shù)方案則有效解決了現(xiàn)有方法存在的這一弊端，其能夠準(zhǔn)確預(yù)測驅(qū)替前緣參數(shù)的變化情況，使得復(fù)雜斷塊油藏高含水后期剩余油二次富集的快速計算得以實現(xiàn)。附圖說明圖1為復(fù)雜斷塊油藏高含水后期剩余油二次富集快速識別方法的流程示意圖；圖2為實施例1中復(fù)雜斷塊油藏的初步物理模型的建立流程示意圖；圖3為實施例1中采油井的動態(tài)歷史特征的擬合方法的流程示意圖；圖4為實施例1中剩余油二次富集前飽和度場的計算流程示意圖；圖5為實施例1中剩余油二次富集過程中各節(jié)點參數(shù)的計算流程示意圖；圖6為實施例1中復(fù)雜斷塊油藏物理模型的空間構(gòu)造示意圖；圖7為實施例1中剩余油富集前后的平面含油飽和度場與垂向含油飽和度場對比圖，圖中A為富集前水平方向的含油飽和度場，C為富集前的垂直方向的含油飽和度場，B為富集后水平方向的含油飽和度場，D為富集后的垂直方向的含油飽和度場；圖8為水平面沿主流線方向的流管示意圖；圖9為流管的空間結(jié)構(gòu)示意圖。具體實施方式為了對本發(fā)明的技術(shù)特征、目的和有益效果有更加清楚的理解，現(xiàn)對本發(fā)明的技術(shù)方案進(jìn)行以下詳細(xì)說明，但不能理解為對本發(fā)明的可實施范圍的限定。實施例本實施例以國內(nèi)某一復(fù)雜斷塊油藏為研究對象，提供一種復(fù)雜斷塊油藏高含水后期剩余油二次富集快速識別方法，如圖1所示，該方法包括：步驟S10，對目標(biāo)油藏進(jìn)行測量，以獲得目標(biāo)油藏的地質(zhì)參數(shù)和井網(wǎng)參數(shù)(如表1所示)。步驟S11，根據(jù)目標(biāo)油藏的地質(zhì)參數(shù)和井網(wǎng)參數(shù)，建立目標(biāo)油藏的初步物理模型，具體包括以下步驟：①根據(jù)獲得的井網(wǎng)參數(shù)，設(shè)置初步物理模型的底面類型，具體包括以下過程：a、根據(jù)獲得的井網(wǎng)參數(shù)，確定目標(biāo)油藏的井網(wǎng)類型、油水井?dāng)?shù)比、井距和排距；b、根據(jù)所述井網(wǎng)類型和油水井?dāng)?shù)比，設(shè)置初步物理模型的底面類型(設(shè)置流程如圖2所示)：當(dāng)目標(biāo)油藏的井網(wǎng)類型為交錯時，將初步物理模型的底面類型設(shè)置為正三角形；當(dāng)目標(biāo)油藏的井網(wǎng)類型為正對，且油水井?dāng)?shù)比＝1時，將初步物理模型的底面類型設(shè)置為矩形；當(dāng)目標(biāo)油藏的井網(wǎng)類型為正對，且油水井?dāng)?shù)≠1時，將初步物理模型的底面類型設(shè)置為正三角形；在本實施例中，初步物理模型的底面類型為正三角形；c、根據(jù)上述確定的井距和排距，設(shè)置正三角形的邊長等于井距，正三角形的高等于排距；d、底面類型確定為正三角形后，將底面沿著垂直方向平移一段距離，該距離位實際儲層的厚度，得到的三維空間結(jié)構(gòu)即為初步物理模型的空間結(jié)構(gòu)，②按照目標(biāo)油藏的實際地質(zhì)結(jié)構(gòu)特征，設(shè)置初步物理模型的沉積韻律、地質(zhì)傾角和滲透率各向異性；其中，沉積韻律是指按顆粒大小、比重的順序先后分層沉積而成成巖的規(guī)律，分為均質(zhì)韻律、正韻律、反韻律三種；均質(zhì)韻律砂層：水動力條件相對穩(wěn)定，層內(nèi)巖性，物性均質(zhì)；正韻律砂層：下部粒度粗，上部粒度細(xì)，反映了沉積環(huán)境水動力條件由下到上變?nèi)?；反韻律砂層：下部粒度?xì)，上部粒度粗，沉積環(huán)境水動力條件由下到上變強；地層傾角是指油水井之間油層走向與水平面之間的夾角；滲透率各向異性是指油層在不同方向上都存在滲透率變化的差異性。③按照獲得的目標(biāo)油藏的實際地質(zhì)參數(shù)，設(shè)置初步物理模型的平均滲透率、初始含油飽和度、孔隙度、束縛水飽和度、殘余油飽和度、滲透率級差、地下原油粘度、儲層厚度、儲層劃分層數(shù)、注采速度等油藏測量中需要獲取的常規(guī)參數(shù)；其中，平均滲透率是指在一定壓差下，巖石允許流體通過的能力大小；初始含油飽和度是油層有效孔隙中含油體積和巖石有效孔隙體積之比，以百分?jǐn)?shù)表示；孔隙度是指巖樣中所有孔隙空間體積之和與該巖樣體積的比值；束縛水飽和度是指由于巖石表面潤濕性導(dǎo)致殘留在巖石孔隙中的最小水體所占儲集層孔隙體積的百分比；殘余油飽和度是指殘余油在巖石孔隙中所占體積的百分?jǐn)?shù)；滲透率級差是指最大滲透率與最小滲透率的比值；地下原油粘度是指在地層條件下原油內(nèi)部某一部分相對于另一部分流動時摩擦阻力的度量。本實施例中目標(biāo)油藏的初步物理模型如圖6所示，模型的底面為正三角形，該模型在X、Y、Z方向每個網(wǎng)格的大小分別為：10m、10m、2m，網(wǎng)格總數(shù)為35×30×25＝26250，油藏低部有邊水。生產(chǎn)井(PRO)位于油藏的高部位，注水井(INC)位于油藏的低部位，油水井正對排列。具體參數(shù)如表1所示表1目標(biāo)油藏的初步物理模型的對應(yīng)參數(shù)參數(shù)名稱參數(shù)取值參數(shù)名稱參數(shù)取值井網(wǎng)類型正三角形孔隙度/％0.28井距/m300地下水粘度/mPa.s0.6排距/m300地下原油粘度/mPa.s7地層傾角/°10地下水比重0.96儲層韻律正韻律地下油比重0.78儲層厚度/m50初始含油飽和度/％0.55油水最大毛管力/bar0.5束縛水飽和度/％0.25平均滲透率/10-3μm2300殘余油飽和度/％0.20垂向平面滲透率比值0.3關(guān)井前含水率/％92滲透率級差5單層注水量/m3500步驟S12，基于目標(biāo)油藏的初步物理模型，根據(jù)流線流管法進(jìn)行擬合計算，以獲得采油井的動態(tài)擬合特征和剩余油二次富集前的飽和度場，并對目標(biāo)油藏的初步物理模型進(jìn)行修正，以得到修正后的物理模型，具體包括以下步驟：①對初步物理模型的儲層參數(shù)特征進(jìn)行預(yù)設(shè)(如圖3所示)，得到預(yù)設(shè)后的初步物理模型：模型的單層厚度＝模型的總厚度/縱向劃分層數(shù)，模型的初始非均質(zhì)性＝縱向非均質(zhì)性；②利用式1和式2分別計算主流線方向上各個節(jié)點的含水率和含水上升率。③計算主流線方向上的含水前緣飽和度及含水前緣位置，具體過程包括：根據(jù)式1和式2可以建立含水飽和度(Sw)與含水率(fw)的關(guān)系曲線(記為Sw—fw關(guān)系曲線)，將Sw(i)點與每個Sw—fw關(guān)系曲線節(jié)點連成線，求該線的導(dǎo)數(shù)，其中最大的導(dǎo)數(shù)值所對應(yīng)的節(jié)點，即為前緣含水飽和度(Swf)；根據(jù)式1和式2可以建立含水飽和度(Sw)與含水上升率(f’w)的關(guān)系(記為Sw—f’w)，由此關(guān)系可以得到f’w(Swf)值，并根據(jù)式3進(jìn)一步計算得到含水前緣位置(xf)。④計算主流線方向上產(chǎn)出端的總含水率，具體過程包括：a、根據(jù)各層的注水量和每條流管的滲透率的比值(各層的注水量和每條流管的滲透率數(shù)據(jù)均可以直接測量獲得)，對每條流管的流量進(jìn)行劈分，以計算得到各層的見水時間(計算公式如式4所示)和各層產(chǎn)出端的含水率(計算公式如式5和式2所示)；b、根據(jù)每條流管的流量對各層產(chǎn)出端的含水率進(jìn)行加權(quán)平均，以得到主流線方向上產(chǎn)出端的總含水率。⑤對上述主流線方向上產(chǎn)出端的總含水率進(jìn)行判斷和修正，直至產(chǎn)出端的總含水率與實際含水率相符合(此時對應(yīng)的縱向劃分層數(shù)和縱向非均質(zhì)性數(shù)據(jù)符合要求，即為最終接受值)，從而擬合得到采油井的擬合含水率。其中，修正過程包括對初步物理模型的縱向劃分層數(shù)和縱向非均質(zhì)性進(jìn)行修改，擬合儲層非均質(zhì)性(如圖3所示)。⑥設(shè)置目標(biāo)含水率(所述目標(biāo)含水率是指滿足關(guān)井富集時機時對應(yīng)的含水率)，根據(jù)上述步驟②至步驟⑤的計算方法，得到采油井的含水率達(dá)到剩余油二次富集要求時對應(yīng)的含油飽和度場，該含油飽和度場為預(yù)測值。⑦對邊界至主流線進(jìn)行含油飽和度插值(即建立飽和度修正關(guān)系式)，以修正上述步驟⑥得到的含油飽和度場，從而得到剩余油二次富集前的飽和度場；其中，上述步驟①-⑤的計算流程如圖3所示，步驟⑥-⑦的計算流程如圖4所示。⑧根據(jù)步驟⑤中確定的縱向劃分層數(shù)、縱向非均質(zhì)性的最終接受值，對目標(biāo)油藏的初步物理模型進(jìn)行修正，得到修正后的物理模型(如圖3所示)。步驟13，基于修正后的物理模型，對剩余油富集過程中的各節(jié)點分別進(jìn)行垂直富集和水平富集的計算，以得到各節(jié)點的儲層飽和度和含水率完成目標(biāo)油藏高含水后期剩余油二次富集的識別(如圖5所示)，具體包括以下步驟：1)垂直富集過程的計算流程如下所述：①按照式6，計算得到每條流管各節(jié)點的垂向高度差，然后分別按照從上往下的方式(計算公式如式7所示)和從下往上(計算公式如式8所示)的方式對每條流管各節(jié)點的垂向高度差進(jìn)行累計求和，得到垂向各節(jié)點位置對應(yīng)的總垂向高度差。②計算垂向上每條流管各節(jié)點完全平衡時的平均含水飽和度(計算公式如式9所示)，將每條流管所在層位的厚度與上一步獲得的垂向上每條流管各節(jié)點的總垂向高度差進(jìn)行大小比對，得到所述飽和度剖面上各節(jié)點的最大含水飽和度平均值(計算公式如式10所示)和最小含水飽和度平均值(計算公式如式11所示)。③將每條流管各節(jié)點完全平衡時的平均含水飽和度，與得到的飽和度剖面上各節(jié)點的最大飽和度平均值和最小飽和度平均值進(jìn)行大小對比，得到剖面平衡態(tài)類型A、剖面平衡態(tài)次級類型B、剖面平衡態(tài)類型的節(jié)點間比例系數(shù)a、剖面平衡態(tài)次級類型的節(jié)點間比例系數(shù)b、剖面平衡態(tài)類型A的平衡位置和剖面平衡態(tài)次級類型B的平衡位置；其中，所述剖面平衡態(tài)類型A是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度大于最大含水飽和度的平均值的節(jié)點個數(shù)；所述剖面平衡態(tài)次級類型B是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度小于最小含水飽的平均值的節(jié)點個數(shù)；剖面平衡態(tài)類型節(jié)點間比例系數(shù)a是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度小于最大含水飽和度的平均值的節(jié)點個數(shù)占該條流管節(jié)點總數(shù)的比值(計算公式如式12所示)；剖面平衡態(tài)次級類型節(jié)點間比例系數(shù)b是指垂向上每條流管各節(jié)點完全平衡時，各節(jié)點的平均含水飽和度大于最大含水飽和度的平均值的節(jié)點個數(shù)占該條流管節(jié)點總數(shù)的比值(計算公式如式13所示)；剖面平衡態(tài)類型A的平衡位置(計算公式如式14所示)和剖面平衡態(tài)次級類型B的平衡位置(計算公式如式15所示)；④計算次生油水界面的剖面含水飽和度(計算公式如式16所示)和剖面平衡態(tài)等效毛管力(計算公式如式17所示)。⑤計算垂向上每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率；具體包括以下過程：a、根據(jù)公式18，計算在富集初始時刻，垂向上每條流管各節(jié)點對應(yīng)的油相、水相的相對滲透率；b、根據(jù)公式19，計算在富集過程中，垂向上每條流管各節(jié)點對應(yīng)的含油飽和度；c、根據(jù)公式20，計算在富集完成時刻，垂向上每條流管各節(jié)點對應(yīng)的油相、水相的相對滲透率。2)水平富集過程的計算流程如下所述①計算在富集初始時刻，平面上每條流管各節(jié)點對應(yīng)的毛管力(計算公式如式21所示)；②計算在富集過程中，平面上每條流管各節(jié)點對應(yīng)的含油飽和度(計算公式如式22所示)；③計算在富集完成時刻，平面上每條流管各節(jié)點對應(yīng)的毛管力(計算公式如式23所示)。3)按照圖5所示的流程：a、預(yù)設(shè)富集時間和富集次數(shù)，在每次循環(huán)計算中，分別進(jìn)行垂向和水平向富集計算，計算出每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率；b、每循環(huán)完一次垂向和水平向富集計算，進(jìn)行判斷是否達(dá)到預(yù)設(shè)的富集次數(shù)，若未達(dá)到，再進(jìn)行步驟a的計算，若達(dá)到富集次數(shù)，則輸出保存最后一次富集計算完的每條流管各節(jié)點的含水飽和度和油、水相的相對滲透率數(shù)據(jù)；c、根據(jù)步驟b所得數(shù)據(jù)，計算出采油井端含水率(計算公式如式1所示)，即是采油井開井含水率，輸出采油井開井含水率，計算結(jié)束。本實施例中，按照上述復(fù)雜斷塊油藏高含水后期剩余油二次富集快速識別方法得到了剩余油富集前后的平面含油飽和度場與垂向含油飽和度場對比圖(如圖7所示，圖7中A為富集前水平方向的含油飽和度場，C為富集前的垂直方向的含油飽和度場，B為富集后水平方向的含油飽和度場，D為富集后的垂直方向的含油飽和度場)。當(dāng)前第1頁1 2 3