[發明專利]基于含水率的稠油邊底水油藏多輪CO2吞吐的評價方法有效
| 申請號: | 201710003679.3 | 申請日: | 2017-01-03 |
| 公開(公告)號: | CN106802985B | 公開(公告)日: | 2019-09-24 |
| 發明(設計)人: | 周建;王佳;成杰;江禮武;劉同敬;第五鵬祥;劉睿 | 申請(專利權)人: | 北京金士力源科技有限公司;中國石油大學(北京) |
| 主分類號: | G06F17/50 | 分類號: | G06F17/50;G06Q50/02;E21B43/16 |
| 代理公司: | 北京三友知識產權代理有限公司 11127 | 代理人: | 姚亮;沈金輝 |
| 地址: | 100083 北京市海淀*** | 國省代碼: | 北京;11 |
| 權利要求書: | 查看更多 | 說明書: | 查看更多 |
| 摘要: | |||
| 搜索關鍵詞: | 基于 含水率 稠油邊 底水 油藏 co2 吞吐 評價 方法 | ||
1.一種基于吞吐后開井含水率的稠油邊底水多輪CO2吞吐的評價方法,其包括以下步驟:
步驟一,測定實際油藏各因素的實際參數,以此為依據建立油藏模型,模擬得到吞吐后開井含水率;
步驟二,以吞吐后開井含水率為評價關聯對象,計算各因素的修正權重;
步驟三,借助油藏數值模擬技術建立各因素與隸屬度的函數關系式,分別得到模擬方案和實際方案的評價矩陣;
步驟四,根據模擬方案評價矩陣和各因素的修正權重,計算各模擬方案所有因素的綜合得分的平均值
步驟五,根據實際方案評價矩陣和各因素的修正權重,計算第j項實際方案所有因素的綜合得分Cj并排序;
步驟六,對比和Cj的大小,給出評價結論;若則評價結果為較好;若則評價結果為較差;
借助油藏數值模擬技術建立各因素與隸屬度的函數關系式,分別得到模擬方案和實際方案的評價矩陣的步驟包括:
根據測得的實際油藏各因素的實際參數,對各因素設計參數變化范圍;
借助油藏數值模擬技術,預測各因素參數變化時對應的吞吐后開井含水率;
依據吞吐后開井含水率的大小,建立各因素與隸屬度的函數關系式;
依據各因素對應的隸屬度,得到模擬方案評價矩陣;
依據建立的各因素與隸屬度的函數關系式,計算實際方案各因素取值對應的隸屬度,得到實際方案評價矩陣;
所述各因素包括隔夾層位置、沉積韻律、出水段長度、出水段位置、水平段位置、平行井距離、吞吐輪次、地層傾角、油層有效厚度、地層非均質性、投產含水、出水段非均質程度、吞吐時機、周期注氣量、注氣速度、燜井時間、開井后采液速度和平行井產液速度中的多種;
所述各因素與隸屬度的函數關系式分別如下:
利用式(1)得到隔夾層位置對應的隸屬度y1a、y1b和y1c:
利用式(2)得到沉積韻律對應的隸屬度y2a、y2b和y2c:
利用式(3)得到出水段長度對應的隸屬度y3a、y3b和y3c:
利用式(4)得到出水段位置對應的隸屬度y4a、y4b和y4c:
其中,A點為水平井的趾端、B點為水平井的跟端;
利用式(5)得到水平段位置對應的隸屬度y5a、y5b和y5c:
利用式(6)得到平行井距離對應的隸屬度y6a、y6b和y6c:
利用式(7)得到吞吐輪次對應的隸屬度y7:
利用式(8)計算地層傾角對應的隸屬度y8a、y8b和y8c:
第1輪:y8a=-0.0065θ2-0.0582θ+1.233
第2輪:y8b=-0.0076θ2-0.242θ+1.657
第3輪:y8c=-0.0065θ2-0.0582θ+1.233 式(8)
其中,θ為地層傾角,°;
利用式(9)計算油層有效厚度對應隸屬度y9a、y9b和y9c:
第1輪:y9a=-0.0826h+1.091
第2輪:y9b=-0.0826h+1.091
第3輪:y9c=-0.0826h+1.091 式(9)
其中,h為油層有效厚度,m;
利用式(10)計算地層非均質性對應的隸屬度y10a、y10b和y10c:
其中,Vk為洛倫茲系數;
利用式(11)計算投產含水對應的隸屬度y11a、y11b和y11c:
第1輪:y11a=-fw0+1
第2輪:y11b=-fw0+1
第3輪:y11c=-fw0+1 式(11)
其中,fw0為投產含水;
利用式(12)計算出水段非均質程度對應的隸屬度y12a、y12b和y12c:
其中,Kv為滲透率極差;
利用式(13)計算吞吐時機對應的隸屬度y13a、y13b和y13c:
其中,f為吞吐時機;
利用式(14)計算周期注氣量對應的隸屬度y14a、y14b和y14c:
其中,Qig為周期注氣量,sm3;
利用式(15)計算注氣速度對應的隸屬度y15a、y15b和y15c:
其中,vig為注氣速度,sm3/d;
利用式(16)計算燜井時間對應的隸屬度y16a、y16b和y16c:
其中,t為燜井時間,d;
利用式(17)計算開井后采液速度對應的隸屬度y17a、y17b和y17c:
其中,vl為開井后采液速度,rm3/d;
利用式(18)計算平行井采液速度對應的隸屬度y18a、y18b和y18c:
其中,ql為平行井采液速度,rm3/d;
根據模擬方案評價矩陣和各因素的修正權重,計算各模擬方案所有因素的綜合得分的平均值的步驟包括:
首先將計算得到的各因素的修正權重組成修正權重向量,計算修正權重向量與模擬方案評價矩陣的乘積,得到模擬方案的得分向量;利用式(II)計算第j項模擬方案所有因素的綜合得分cj:
其中,Ai為第i項因素的修正權重;i=1,2,……,n(n≤N);N為方案評價所需因素總數;n為現場資料實際提供的因素總數;bji為第j項模擬方案的第i項因素對應的隸屬度;j=1,2,……,m;m為模擬方案總數;cj為第j項模擬方案所有因素的綜合得分;
然后利用式(III)計算各模擬方案所有因素的綜合得分的平均值作為評價界限的參考值:
其中,為各模擬方案所有因素的綜合得分的平均值;cj為第j項模擬方案所有因素的綜合得分;j=1,2,……,m;m為模擬方案總數;
根據實際方案評價矩陣和各因素的修正權重,計算實際方案所有因素的綜合得分并排序的步驟包括:
首先將計算得到的各因素的修正權重組成修正權重向量,計算修正權重向量與實際方案評價矩陣的乘積,得到實際方案的得分向量;利用式(IV)計算第j項實際方案所有因素的綜合得分Cj:
其中,Ai為第i項因素的修正權重;i=1,2,……,n(n≤N);N為方案評價所需因素總數;n為現場資料實際提供的因素總數;Bji為第j項實際方案的第i項因素對應的隸屬度;j=1,2,……,M;M為實際方案總數;Cj為第j項實際方案所有因素的綜合得分;
然后對第j項實際方案所有因素的綜合得分按Cj大小進行排序。
2.根據權利要求1所述的評價方法,其特征在于,各因素的修正權重是采用式(I)計算得到的:
其中,ai為第i項因素的經驗權重;Ai為第i項因素的修正權重;i=1,2,……,n(n≤N);N為方案評價所需的因素總數;n為現場資料實際提供的因素總數。
該專利技術資料僅供研究查看技術是否侵權等信息,商用須獲得專利權人授權。該專利全部權利屬于北京金士力源科技有限公司;中國石油大學(北京),未經北京金士力源科技有限公司;中國石油大學(北京)許可,擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作,請聯系【客服】
本文鏈接:http://www.szxzyx.cn/pat/books/201710003679.3/1.html,轉載請聲明來源鉆瓜專利網。





