本發(fā)明涉及油氣田開發(fā)，尤其涉及一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法及系統(tǒng)。

背景技術(shù)：

1、整裝油藏特高含水后期受儲(chǔ)層非均質(zhì)性及長(zhǎng)期注水開發(fā)影響，油藏流線波及范圍固化，如何識(shí)別流動(dòng)單元，有效轉(zhuǎn)變流線打破原有固定的流動(dòng)體系，是特高含水整裝油藏注水開發(fā)當(dāng)前最迫切的難題，對(duì)整裝油藏持續(xù)提高水驅(qū)采收率尤其重要。

2、中國(guó)發(fā)明專利cn201410131727.3一種碳酸鹽巖水平井流動(dòng)單元?jiǎng)澐址椒?。該方法井間連通關(guān)系判斷，儲(chǔ)集體特征及泄油面積確定，并通過(guò)數(shù)值試井模型建立進(jìn)行驗(yàn)證，能夠精確劃分流動(dòng)單元。

3、中國(guó)發(fā)明專利cn201911000486.8一種基于滲流實(shí)驗(yàn)的深層非均質(zhì)灘相氣藏流動(dòng)單元?jiǎng)澐址椒?。該方法選取不同儲(chǔ)層類型的樣品開展?jié)B流特征實(shí)驗(yàn)，明確不同樣品的滲流特征；分析樣品孔隙度、滲透率、r35與滲流指數(shù)的相關(guān)性，優(yōu)選流動(dòng)單元?jiǎng)澐株P(guān)鍵參數(shù)；在滲流指數(shù)與滲透率、r35的交匯圖中，將具有相近fi值的樣品作為同一流動(dòng)單元，并讀取對(duì)應(yīng)的滲透率及r35值，完成流動(dòng)單元分類與閥門值的確定。

4、中國(guó)發(fā)明專利cn201910154632.6一種特高含水期油藏動(dòng)態(tài)流動(dòng)單元?jiǎng)澐址椒āＴ摲椒ǜ鶕?jù)油田已有的相滲曲線，回歸得到不同絕對(duì)滲透率下的相滲曲線表達(dá)式，進(jìn)而計(jì)算得到各網(wǎng)格點(diǎn)的水相相對(duì)滲透；根據(jù)滲流系數(shù)定義，計(jì)算得到模型中每個(gè)網(wǎng)格點(diǎn)的水相滲流系數(shù)；繪制水相滲流系數(shù)半對(duì)數(shù)累積曲線圖，找出兩個(gè)平行段(斜率較小，幾乎為0)以及兩個(gè)平行段間斜率較大的一段，作為動(dòng)態(tài)流動(dòng)單元的劃分界限；對(duì)各個(gè)網(wǎng)格的水相滲流系數(shù)進(jìn)行分區(qū)判定，進(jìn)而劃分流動(dòng)單元。

5、中國(guó)發(fā)明專利cn201610890913.4流動(dòng)單元的劃分方法和裝置及井間對(duì)比方法和裝置。將油氣儲(chǔ)層層段劃分為多個(gè)儲(chǔ)層細(xì)分層段，并計(jì)算各儲(chǔ)層細(xì)分層段的累計(jì)產(chǎn)能百分比和累計(jì)儲(chǔ)能百分比；并根據(jù)所述sml圖對(duì)流動(dòng)單元進(jìn)行一次劃分，在單井上獲得多個(gè)初步流動(dòng)單元；計(jì)算每相鄰的兩個(gè)所述初步流動(dòng)單元之間的流動(dòng)單元差異指數(shù)，根據(jù)所述流動(dòng)單元差異指數(shù)對(duì)所述初步流動(dòng)單元進(jìn)行二次劃分，獲得多個(gè)流動(dòng)單元。

6、目前流動(dòng)單元的研究主要針對(duì)儲(chǔ)層的靜態(tài)特征或開發(fā)動(dòng)態(tài)統(tǒng)計(jì)規(guī)律進(jìn)行流動(dòng)單元的識(shí)別與劃分，但流動(dòng)單元應(yīng)受動(dòng)靜態(tài)參數(shù)共同影響，油藏的沉積環(huán)境、儲(chǔ)層的物性特征、非均質(zhì)特征，長(zhǎng)期注水開發(fā)油藏不同區(qū)域、層位驅(qū)替程度的差異，現(xiàn)開發(fā)技術(shù)政策條件下油藏驅(qū)替強(qiáng)度的差異等均影響流動(dòng)單元的分布。本專利提出的一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法，將儲(chǔ)層沉積、地質(zhì)、開發(fā)歷史、開發(fā)現(xiàn)狀等各類動(dòng)靜態(tài)信息充分利用，提出了能夠反映儲(chǔ)層”沉積特征、非均質(zhì)特征、驅(qū)替過(guò)程和特征”的流動(dòng)單元識(shí)別劃分方法，建立了一套注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法，從而為整裝特高含水后期油藏開發(fā)調(diào)整提供了有力的支撐。

技術(shù)實(shí)現(xiàn)思路

1、鑒于上述問(wèn)題，提出了本發(fā)明以便提供克服上述問(wèn)題或者至少部分地解決上述問(wèn)題的一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法。

2、根據(jù)本發(fā)明的一個(gè)方面，提供了一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法，所述判識(shí)方法包括：

3、步驟s1：優(yōu)選流動(dòng)單元?jiǎng)屿o態(tài)表征參數(shù)集合；

4、步驟s2：采用層次分析法，確定流動(dòng)單元表征參數(shù)權(quán)重；

5、步驟s3：建立并計(jì)算研究區(qū)流動(dòng)單元綜合表征參數(shù)--流動(dòng)指數(shù)；

6、步驟s4：繪制流動(dòng)單元指數(shù)分布曲線及流動(dòng)指數(shù)概率分布圖，研究單元流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)。

7、可選的，所述步驟s1：優(yōu)選流動(dòng)單元?jiǎng)屿o態(tài)表征參數(shù)集合具體包括：

8、根據(jù)參數(shù)間遞進(jìn)關(guān)系、等效關(guān)系分析，優(yōu)選流動(dòng)單元表征參數(shù)集為巖相、滲透率、流動(dòng)帶指數(shù)、累積過(guò)水倍數(shù)、液相流速。

9、可選的，所述步驟s2：采用層次分析法，確定流動(dòng)單元表征參數(shù)權(quán)重具體包括：

10、步驟s101：構(gòu)建流動(dòng)單元流場(chǎng)表征結(jié)構(gòu)模型；

11、步驟s102：構(gòu)建判斷矩陣；

12、步驟s103：基于構(gòu)建的各判斷矩陣，對(duì)每個(gè)判斷矩陣計(jì)算權(quán)重并逐一進(jìn)行一致性檢驗(yàn)；

13、步驟s104：基于全部通過(guò)一致性檢驗(yàn)的各判斷矩陣，得到每個(gè)指標(biāo)的相對(duì)權(quán)重，根據(jù)流動(dòng)單元流場(chǎng)表征的層級(jí)結(jié)構(gòu)，對(duì)各層級(jí)的權(quán)重進(jìn)行合并計(jì)算，得到流動(dòng)單元流場(chǎng)表征中各項(xiàng)指標(biāo)的權(quán)重賦值。

14、可選的，所述步驟s101：構(gòu)建流動(dòng)單元流場(chǎng)表征結(jié)構(gòu)模型具體包括：

15、根據(jù)優(yōu)選的流動(dòng)單元表征參數(shù)構(gòu)成流動(dòng)單元流場(chǎng)表征結(jié)構(gòu)模型。

16、可選的，所述步驟s102：構(gòu)建判斷矩陣具體包括：

17、將流動(dòng)單元流場(chǎng)表征結(jié)構(gòu)模型中的指標(biāo)進(jìn)行比較，構(gòu)建判斷矩陣e。

18、可選的，所述將流動(dòng)單元流場(chǎng)表征結(jié)構(gòu)模型中的指標(biāo)進(jìn)行比較的規(guī)則為：

19、若指標(biāo)m和指標(biāo)n同等重要?jiǎng)t標(biāo)為1；

20、若指標(biāo)m比指標(biāo)n稍微重要?jiǎng)t標(biāo)為3；

21、若指標(biāo)m比指標(biāo)n明顯重要?jiǎng)t標(biāo)為5；

22、若指標(biāo)m指標(biāo)n極其重要?jiǎng)t標(biāo)為7；

23、若指標(biāo)m比指標(biāo)n強(qiáng)烈重要?jiǎng)t標(biāo)為9；2、4、6、8表示上述相鄰判斷情形的中間值；

24、若后者比前者重要?jiǎng)t取相應(yīng)的倒數(shù)；判斷矩陣e中各指標(biāo)權(quán)重的計(jì)算公式為：

25、w＝∑(aij/aj)/n

26、其中，aij為判斷矩陣e中指標(biāo)a的第i行第j列的值，aj為判斷矩陣e中指標(biāo)a的第j列所有值的和，n為該判斷矩陣e的行列階數(shù)。

27、可選的，所述步驟s103：基于構(gòu)建的各判斷矩陣，對(duì)每個(gè)判斷矩陣計(jì)算權(quán)重并逐一進(jìn)行一致性檢驗(yàn)具體包括：

28、構(gòu)建判斷矩陣e得到各層的權(quán)重，計(jì)算判斷矩陣e的一致性指標(biāo)ci：

29、

30、其中，λmax是判斷矩陣e的最大特征值，n為判斷矩陣e的行列階數(shù)；最后計(jì)算判斷矩陣e的一致性比例cr：

31、

32、其中，ri是隨機(jī)一致性指標(biāo)；當(dāng)cr<0.1時(shí)，則判斷矩陣e滿足一致性檢驗(yàn)，否則，則不滿足一致性檢驗(yàn)，重新構(gòu)建該層級(jí)的判斷矩陣e。

33、可選的，所述步驟s3：建立并計(jì)算研究區(qū)流動(dòng)單元綜合表征參數(shù)--流動(dòng)指數(shù)具體包括：

34、流動(dòng)單元計(jì)算公式如下：

35、fl＝0.052st’+0.158k’+0.211fzi’+0.316qi’+0.263vl’(1)

36、式中：st’為標(biāo)準(zhǔn)化后巖相，f；k’為標(biāo)準(zhǔn)化后滲透率，t/d；fzi’為標(biāo)準(zhǔn)化后流動(dòng)帶指數(shù)，f；qi’為標(biāo)準(zhǔn)化后累積過(guò)水倍數(shù)，f；vl為標(biāo)準(zhǔn)化后液相流速，f。

37、可選的，所述步驟s4：繪制流動(dòng)單元指數(shù)分布曲線及流動(dòng)指數(shù)概率分布圖，研究單元流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)具體包括：

38、根據(jù)計(jì)算的流動(dòng)指數(shù)，統(tǒng)計(jì)注水開發(fā)油藏流動(dòng)單元數(shù)值模擬結(jié)果中所有網(wǎng)格的流動(dòng)指數(shù)值，并按照流動(dòng)指數(shù)由小到大的順序統(tǒng)計(jì)對(duì)應(yīng)網(wǎng)格數(shù)的累加值，以累計(jì)網(wǎng)格個(gè)數(shù)為橫坐標(biāo)，流動(dòng)指數(shù)為縱坐標(biāo)繪制流動(dòng)單元指數(shù)分布曲線；

39、根據(jù)曲線形態(tài)特征，選取曲線形態(tài)變化點(diǎn)位置所對(duì)應(yīng)的流動(dòng)指數(shù)作為流動(dòng)單元分級(jí)界限，為流動(dòng)單元分級(jí)標(biāo)準(zhǔn)提供依據(jù)。

40、可選的，所述步驟s4：繪制流動(dòng)單元指數(shù)分布曲線及流動(dòng)指數(shù)概率分布圖，研究單元流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)還包括：

41、根據(jù)流動(dòng)指數(shù)值，計(jì)算不同流動(dòng)指數(shù)的概率值，以流動(dòng)指數(shù)為橫坐標(biāo)，以不同流動(dòng)指數(shù)概率值為縱坐標(biāo)繪制流動(dòng)指數(shù)概率分布圖，根據(jù)流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)單元流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)。

42、本發(fā)明還提供了一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)系統(tǒng)，應(yīng)用上述所述的一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法，所述判識(shí)系統(tǒng)包括：

43、動(dòng)靜態(tài)表征參數(shù)優(yōu)選模塊，用于優(yōu)選流動(dòng)單元?jiǎng)屿o態(tài)表征參數(shù)集合；

44、層次分析法模塊，用于采用層次分析法，確定流動(dòng)單元表征參數(shù)權(quán)重；

45、計(jì)算模塊，用于建立并計(jì)算研究區(qū)流動(dòng)單元綜合表征參數(shù)--流動(dòng)指數(shù)；

46、分類標(biāo)準(zhǔn)確定模塊，用于繪制流動(dòng)單元指數(shù)分布曲線及流動(dòng)指數(shù)概率分布圖，研究單元流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)。

47、本發(fā)明提供的一種注水開發(fā)油藏流動(dòng)單元數(shù)值模擬判識(shí)方法及系統(tǒng)，所述判識(shí)方法包括：步驟s1：優(yōu)選流動(dòng)單元?jiǎng)屿o態(tài)表征參數(shù)集合；步驟s2：采用層次分析法，確定流動(dòng)單元表征參數(shù)權(quán)重；步驟s3：建立并計(jì)算研究區(qū)流動(dòng)單元綜合表征參數(shù)--流動(dòng)指數(shù)；步驟s4：繪制流動(dòng)單元指數(shù)分布曲線及流動(dòng)指數(shù)概率分布圖，研究單元流動(dòng)指數(shù)概率分布規(guī)律，確定流動(dòng)指數(shù)分級(jí)分類標(biāo)準(zhǔn)。能夠更加準(zhǔn)確地識(shí)別儲(chǔ)層流動(dòng)單元為油藏潛力認(rèn)識(shí)、剩余油挖潛及開發(fā)調(diào)整方案編制提供技術(shù)支撐。

48、上述說(shuō)明僅是本發(fā)明技術(shù)方案的概述，為了能夠更清楚了解本發(fā)明的技術(shù)手段，而可依照說(shuō)明書的內(nèi)容予以實(shí)施，并且為了讓本發(fā)明的上述和其它目的、特征和優(yōu)點(diǎn)能夠更明顯易懂，以下特舉本發(fā)明的具體實(shí)施方式。