本發(fā)明屬于co2驅(qū)油藏開(kāi)發(fā)
技術(shù)領(lǐng)域：
，具體涉及一種co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法。
背景技術(shù)：
：隨著溫室氣體減排和地質(zhì)埋存日益受到全世界的關(guān)注，co2驅(qū)油技術(shù)在國(guó)內(nèi)越來(lái)越受到重視。制約co2驅(qū)開(kāi)發(fā)效果的主要因素為氣竄通道導(dǎo)致的co2氣竄。由于co2的黏度低，氣竄通道和不利的流度比使注入的co2普遍地繞過(guò)被驅(qū)替的油相而發(fā)生竄流。當(dāng)壓力梯度相同時(shí)，co2比油和水的移動(dòng)速度快。因此，氣竄通道的存在直接降低了co2波及效率，嚴(yán)重影響co2驅(qū)油藏開(kāi)發(fā)效果。目前co2驅(qū)油藏氣竄通道識(shí)別多從靜態(tài)入手，通過(guò)分析儲(chǔ)層發(fā)育特征，尋找井間高滲層。但國(guó)內(nèi)大量礦場(chǎng)試驗(yàn)已經(jīng)證明co2驅(qū)油藏，特別是低滲透油藏，其氣竄受人工裂縫溝通隱裂縫形成的縫網(wǎng)控制。縫網(wǎng)的形成受人工儲(chǔ)層改造、高壓注入改造作用雙重控制，單純從孔隙度、滲透率等靜態(tài)參數(shù)分析，只能定性描述氣竄通道，且描述結(jié)果與動(dòng)態(tài)數(shù)據(jù)吻合度低，無(wú)法準(zhǔn)確指導(dǎo)co2驅(qū)油藏的開(kāi)發(fā)管理。技術(shù)實(shí)現(xiàn)要素：為解決上述問(wèn)題，本發(fā)明的目的是提供一種co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法。為達(dá)到上述目的，本發(fā)明提供了一種co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，其中，該方法包括以下步驟：獲取實(shí)際氣油比曲線特征參數(shù)，所述氣油比曲線特征參數(shù)包括氣竄對(duì)應(yīng)注氣時(shí)間、最大氣油比和氣油比曲線半坡值時(shí)間比例；通過(guò)數(shù)值模擬，獲得理論氣油比曲線，建立理論氣油比曲線和氣竄通道特征參數(shù)的關(guān)聯(lián)關(guān)系；根據(jù)上述關(guān)聯(lián)關(guān)系，建立反演指標(biāo)體系；利用模擬得到的實(shí)際氣油比曲線特征參數(shù)和反演指標(biāo)體系對(duì)應(yīng)關(guān)系，實(shí)現(xiàn)實(shí)際氣竄通道特征參數(shù)的動(dòng)態(tài)反演。本發(fā)明提供co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，在co2驅(qū)油藏生產(chǎn)動(dòng)態(tài)數(shù)據(jù)抽象方法的基礎(chǔ)上，結(jié)合數(shù)值模擬技術(shù)，利用分級(jí)劃分抽象后的生產(chǎn)動(dòng)態(tài)數(shù)據(jù)，實(shí)現(xiàn)了co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演，對(duì)co2驅(qū)油藏，特別是低滲透油藏的科學(xué)開(kāi)發(fā)和有效調(diào)整具有重要意義。該方法的創(chuàng)新之處在于提供的是一種從動(dòng)態(tài)數(shù)據(jù)(即氣油比曲線)反演氣竄通道的方法，有別于原先通過(guò)靜態(tài)數(shù)據(jù)(即地層信息，例如滲透率)或測(cè)試手段(微地震、示蹤劑)分析氣竄通道。該方法的最大優(yōu)勢(shì)在于，氣油比曲線非常容易獲取，且數(shù)據(jù)密度大(一般為1個(gè)/天)，比靜態(tài)數(shù)據(jù)分析精度高、比測(cè)試手段成本低。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，在所述獲取實(shí)際氣油比曲線特征參數(shù)步驟中，可以采用以下3個(gè)過(guò)程獲得待反演氣竄通道的實(shí)際氣油比曲線特征參數(shù)：過(guò)程1，利用公式1，根據(jù)采油井氣竄日期dgor和注氣井注co2日期dinj，計(jì)算氣竄對(duì)應(yīng)注氣時(shí)間t：t＝dgor-dinj公式1；該過(guò)程中，所述實(shí)際氣油比曲線指co2驅(qū)油藏采油井實(shí)測(cè)氣油比曲線；所述氣油比曲線特征參數(shù)指氣竄對(duì)應(yīng)注氣時(shí)間t、最大氣油比gormax和氣油比曲線半坡值時(shí)間比例r；所述采油井氣竄日期dgor指注co2后，采油井氣油比開(kāi)始連續(xù)3個(gè)月大于原始?xì)庥捅鹊娜掌?；所述注氣井注co2日期dinj指采油井來(lái)氣方向注氣井開(kāi)始注入co2的日期；所述原始?xì)庥捅戎缸o2前采油井氣油比平均值。過(guò)程2，計(jì)算最大氣油比gormax：該過(guò)程中，所述最大氣油比gormax指采油井氣竄日期dgor后有效氣油比的最大值；所述有效氣油比指除去異常氣油比數(shù)值后的采油井氣油比曲線；所述異常氣油比指由于含水率突變或產(chǎn)氣數(shù)據(jù)計(jì)量不準(zhǔn)導(dǎo)致的非連續(xù)變化所對(duì)應(yīng)的氣油比。過(guò)程3，利用公式2，分別根據(jù)氣油比半坡值gormid對(duì)應(yīng)注氣時(shí)間tgormid、最大氣油比對(duì)應(yīng)注氣時(shí)間tgormax計(jì)算半坡值時(shí)間比例r：該過(guò)程中，所述氣油比半坡值gormid指最大氣油比gormax與原始?xì)庥捅萭orii的算術(shù)平均值；所述氣油比半坡值gormid對(duì)應(yīng)注氣時(shí)間tgormid，指氣油比半坡值日期dgormid與注氣井注co2時(shí)間dinj之差，如公式3所示：tgormid＝dgormid-dinj公式3。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，所述通過(guò)數(shù)值模擬，獲得理論氣油比曲線，建立理論氣油比曲線和氣竄通道特征參數(shù)的關(guān)聯(lián)關(guān)系的具體方式為：根據(jù)收集到的實(shí)際co2驅(qū)油藏的靜態(tài)參數(shù)、流體參數(shù)和開(kāi)發(fā)歷程，利用數(shù)值模擬軟件計(jì)算不同氣竄通道特征參數(shù)對(duì)應(yīng)的理論氣油比曲線特征參數(shù)。所述理論氣油比曲線為數(shù)值模擬軟件計(jì)算得到的氣油比曲線。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，所述靜態(tài)參數(shù)包括油藏埋深、地層厚度、孔隙度場(chǎng)、滲透率場(chǎng)、凈毛比場(chǎng)、井網(wǎng)類型、單井射孔厚度、原始地層壓力和原始地層溫度；優(yōu)選地，所述流體參數(shù)包括地下原油組分組成、油氣水三相相滲曲線和毛管力曲線；優(yōu)選地，所述開(kāi)發(fā)歷程為實(shí)際co2驅(qū)油藏注co2前注入井和采油井工作制度。另外，數(shù)值模擬可采用一般的商業(yè)油藏?cái)?shù)值模擬軟件(例如eclipse等)，當(dāng)然也可以使用自行開(kāi)發(fā)的油藏?cái)?shù)值模擬軟件。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，所述氣竄通道特征參數(shù)包括氣竄通道類型、厚度和滲透率級(jí)差；進(jìn)一步優(yōu)選地，所述氣竄通道類型包括高滲層和高滲通道。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，所述根據(jù)上述關(guān)聯(lián)關(guān)系，建立反演指標(biāo)體系的具體方式為：根據(jù)不同氣竄通道特征參數(shù)對(duì)應(yīng)的理論氣油比曲線特征參數(shù)的分布，按照等概率分布原則對(duì)氣竄對(duì)應(yīng)的注氣時(shí)間、最大氣油比和氣油比曲線半坡值時(shí)間比例進(jìn)行分類，建立起反演指標(biāo)體系。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，在所述等概率分布原則下，可以按照實(shí)際需求對(duì)氣油比曲線特征參數(shù)進(jìn)行分類。在發(fā)明的一種優(yōu)選實(shí)施方式中，在所述等概率分布原則下，氣竄對(duì)應(yīng)注氣時(shí)間被分為四類；最大氣油比被分為四類；氣油比曲線半坡值時(shí)間被分為小于40％、介于40％-60％和大于60％三類。在上述co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法中，優(yōu)選地，所述油藏為低滲透油藏。本發(fā)明提供的co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，通過(guò)建立基于采油井生產(chǎn)動(dòng)態(tài)數(shù)據(jù)的co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，解決了現(xiàn)有的靜態(tài)分析方法只能定性描述氣竄通道，且與描述結(jié)果動(dòng)態(tài)數(shù)據(jù)吻合度低，無(wú)法準(zhǔn)確指導(dǎo)co2驅(qū)油藏開(kāi)發(fā)管理的問(wèn)題。附圖說(shuō)明圖1為實(shí)施例2中的中國(guó)東部某co2驅(qū)油藏w1井氣油比曲線圖。具體實(shí)施方式為了對(duì)本發(fā)明的技術(shù)特征、目的和有益效果有更加清楚的理解，現(xiàn)對(duì)本發(fā)明的技術(shù)方案進(jìn)行以下詳細(xì)說(shuō)明，但不能理解為對(duì)本發(fā)明的可實(shí)施范圍的限定。實(shí)施例1本實(shí)施例提供了一種co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，該方法包括以下步驟：步驟一：獲取實(shí)際氣油比曲線特征參數(shù)，所述氣油比曲線特征參數(shù)包括氣竄對(duì)應(yīng)注氣時(shí)間、最大氣油比和氣油比曲線半坡值時(shí)間比例。在本實(shí)施例中，實(shí)際氣油比曲線特征參數(shù)通過(guò)以下過(guò)程獲得：過(guò)程1，利用公式1，根據(jù)采油井氣竄日期dgor和注氣井注co2日期dinj，計(jì)算氣竄對(duì)應(yīng)注氣時(shí)間t：t＝dgor-dinj公式1；該過(guò)程中，實(shí)際氣油比曲線指co2驅(qū)油藏采油井實(shí)測(cè)氣油比曲線；氣油比曲線特征參數(shù)指氣竄對(duì)應(yīng)注氣時(shí)間t、最大氣油比gormax和氣油比曲線半坡值時(shí)間比例r；采油井氣竄日期dgor指注co2后，采油井氣油比開(kāi)始連續(xù)3個(gè)月大于原始?xì)庥捅鹊娜掌?；注氣井注co2日期dinj指采油井來(lái)氣方向注氣井開(kāi)始注入co2的日期；原始?xì)庥捅戎缸o2前采油井氣油比平均值。過(guò)程2，計(jì)算最大氣油比gormax：該過(guò)程中，最大氣油比gormax指采油井氣竄日期dgor后有效氣油比的最大值；有效氣油比指除去異常氣油比數(shù)值后的采油井氣油比曲線；異常氣油比指由于含水率突變或產(chǎn)氣數(shù)據(jù)計(jì)量不準(zhǔn)導(dǎo)致的非連續(xù)變化所對(duì)應(yīng)的氣油比。過(guò)程3，利用公式2，分別根據(jù)氣油比半坡值gormid對(duì)應(yīng)注氣時(shí)間tgormid、最大氣油比對(duì)應(yīng)注氣時(shí)間tgormax計(jì)算半坡值時(shí)間比例r：該過(guò)程中，氣油比半坡值gormid指最大氣油比gormax與原始?xì)庥捅萭orii的算術(shù)平均值；氣油比半坡值gormid對(duì)應(yīng)注氣時(shí)間tgormid，指氣油比半坡值日期dgormid與注氣井注co2時(shí)間dinj之差，如公式3所示：tgormid＝dgormid-dinj公式3。步驟二：通過(guò)數(shù)值模擬，獲得理論氣油比曲線，建立理論氣油比曲線和氣竄通道特征參數(shù)的關(guān)聯(lián)關(guān)系。該步驟具體為：根據(jù)收集到的實(shí)際co2驅(qū)油藏的靜態(tài)參數(shù)、流體參數(shù)和開(kāi)發(fā)歷程，利用數(shù)值模擬軟件計(jì)算不同氣竄通道特征參數(shù)對(duì)應(yīng)的理論氣油比曲線特征參數(shù)；所述理論氣油比曲線為數(shù)值模擬軟件計(jì)算得到的氣油比曲線。其中，氣竄通道特征參數(shù)包括氣竄通道類型、厚度和滲透率級(jí)差；氣竄通道類型包括高滲層和高滲通道。靜態(tài)參數(shù)包括油藏埋深、地層厚度、孔隙度場(chǎng)、滲透率場(chǎng)、凈毛比場(chǎng)、井網(wǎng)類型、單井射孔厚度、原始地層壓力和原始地層溫度；流體參數(shù)包括地下原油組分組成、油氣水三相相滲曲線和毛管力曲線；開(kāi)發(fā)歷程為實(shí)際co2驅(qū)油藏注co2前注入井和采油井工作制度；使用的一般的商業(yè)油藏?cái)?shù)值模擬軟件即可。步驟三：根據(jù)步驟二獲得的關(guān)聯(lián)關(guān)系，建立反演指標(biāo)體系。該步驟具體為：根據(jù)不同氣竄通道特征參數(shù)對(duì)應(yīng)的理論氣油比曲線特征參數(shù)的分布，按照等概率分布原則對(duì)氣竄對(duì)應(yīng)的注氣時(shí)間、最大氣油比和氣油比曲線半坡值時(shí)間比例進(jìn)行分類，建立起反演指標(biāo)體系。其中，在等概率分布原則下，氣竄對(duì)應(yīng)注氣時(shí)間被分為四類；最大氣油比被分為四類；氣油比曲線半坡值時(shí)間被分為小于40％、介于40％-60％和大于60％三類。步驟四：利用模擬得到的實(shí)際氣油比曲線特征參數(shù)和反演指標(biāo)體系對(duì)應(yīng)關(guān)系，實(shí)現(xiàn)實(shí)際氣竄通道特征參數(shù)的動(dòng)態(tài)反演。本實(shí)施例提供的co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法，可以解決了現(xiàn)有的靜態(tài)分析方法只能定性描述氣竄通道，且與動(dòng)態(tài)數(shù)據(jù)吻合度低，無(wú)法準(zhǔn)確指導(dǎo)co2驅(qū)油藏開(kāi)發(fā)管理的問(wèn)題。實(shí)施例2本實(shí)施例提供了采用實(shí)施例1中co2驅(qū)油藏氣竄通道動(dòng)態(tài)反演方法的應(yīng)用實(shí)驗(yàn)，具體為：圖1為中國(guó)東部某co2驅(qū)油藏w1井氣油比曲線。w1井氣竄日期為2011年4月，對(duì)應(yīng)注氣井注co2日期為2008年5月，計(jì)算得到氣竄對(duì)應(yīng)注氣時(shí)間為1065天。最大氣油比為2012年8月對(duì)應(yīng)的263.2m3/t。原始?xì)庥捅燃s42.3m3/t，對(duì)應(yīng)半坡值為152.8m3/t，氣油比半坡值日期為2012年3月15日，計(jì)算半坡值時(shí)間比例為71.5％。針對(duì)該東部co2驅(qū)油藏，根據(jù)其靜態(tài)參數(shù)、流體參數(shù)和開(kāi)發(fā)歷程，利用數(shù)值模擬軟件(eclipse)計(jì)算得到了不同氣竄通道特征參數(shù)對(duì)應(yīng)的理論氣油比曲線特征參數(shù)。靜態(tài)參數(shù)如表1所示。流體參數(shù)中原油組分組成如表2所示，油氣水三相相滲曲線和毛管力曲線參考《油藏物理》(秦積舜著)中低滲油藏典型曲線取值。開(kāi)發(fā)歷程為4口采油井定液6rm3/d生產(chǎn)，衰竭開(kāi)采至地層壓力降至15mpa后注入井以35rm3/d的速度注入co2。氣竄通道特征參數(shù)設(shè)計(jì)如表3所示。利用商業(yè)軟件eclipse的組分模型模擬，建立的反演指標(biāo)體系如表4所示。表1東部co2驅(qū)油藏靜態(tài)參數(shù)參數(shù)取值參數(shù)取值油藏埋深/m3000凈毛比場(chǎng)均質(zhì)，1地層厚度/m11井網(wǎng)類型五點(diǎn)井網(wǎng)孔隙度場(chǎng)/f均質(zhì)，0.12單井射孔厚度/m11滲透率場(chǎng)/10-3μm2均質(zhì)，1.5原始地層壓力/mpa20.5原始地層溫度/℃98.9表2東部co2驅(qū)油藏流體參數(shù)(原油組分組成)擬組分co2c1+n2c2c3c4c5+c6c7-c10c11-c28c29+組分組成0.003430.18710.0590.038430.016960.03950.171850.377790.10594表3氣竄通道特征參數(shù)設(shè)計(jì)氣竄通道類型氣竄通道厚度/m滲透率級(jí)差高滲層、高滲通道0.25，1，52.5，10，40表4建立的反演指標(biāo)體系根據(jù)反演指標(biāo)體系，w1井對(duì)應(yīng)的歸類為“3-4-3”。查詢模擬得到的實(shí)際氣油比曲線特征參數(shù)和反演指標(biāo)體系對(duì)應(yīng)關(guān)系，得到氣竄通道類型為高滲層，厚度為5m，滲透率級(jí)差為40。當(dāng)前第1頁(yè)12