白丝美女被狂躁免费视频网站,500av导航大全精品,yw.193.cnc爆乳尤物未满,97se亚洲综合色区,аⅴ天堂中文在线网官网

首頁 / 專利庫 / 地球科學(xué) / 巖石 / 沉積巖 / 石灰?guī)r / 深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法

深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法

閱讀:562發(fā)布:2020-05-16

專利匯可以提供深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法專利檢索,專利查詢,專利分析的服務(wù)。并且本 發(fā)明 公開了深層 石灰?guī)r 儲層酸壓酸液有效作用距離的數(shù)值計算方法,包括:(1)計算施工結(jié)束時的井底 溫度 ;(2)計算施工結(jié)束時的裂縫縫長,在縫長方向劃分網(wǎng)格單元并計算各單元的平均縫寬,得到各單元酸液流速;(3)以井底溫度為初始條件,計算各單元不考慮酸巖反應(yīng)熱的酸液溫度以及酸液濃度;(4)計算各單元的酸巖摩爾反應(yīng)熱;(5)計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,建立酸液濃度分布曲線,定義酸液初始注入濃度的10%作為殘酸極限濃度,以殘酸極限濃度對應(yīng)的縫長為酸液有效作用距離。本發(fā)明考慮了酸壓過程中酸液受到的熱效應(yīng),提高了酸液有效作用距離計算的準(zhǔn)確性,對深層石灰?guī)r儲層酸壓設(shè)計具有指導(dǎo)意義。,下面是深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法專利的具體信息內(nèi)容。

1.深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,依次包括以下步驟:
(1)采用Eickmeier井筒溫度場模型模擬井筒溫度場并計算施工結(jié)束時的井底溫度;
(2)采用PKN公式計算施工結(jié)束時的裂縫縫長,在縫長方向劃分網(wǎng)格單元并計算各單元的平均縫寬,再計算各單元酸液流速;
(3)以步驟(1)計算的井底溫度為初始條件,通過步驟(2)計算的酸液流速計算各單元不考慮酸巖反應(yīng)熱的酸液溫度以及酸液濃度;
(4)通過步驟(3)中計算的酸液溫度和濃度計算各單元的酸巖摩爾反應(yīng)熱;
(5)將步驟(4)計算的酸巖摩爾反應(yīng)熱代入步驟(3)中,計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,建立酸液濃度分布曲線,定義酸液初始注入濃度的10%作為殘酸極限濃度,根據(jù)分布曲線以殘酸極限濃度對應(yīng)的縫長為酸液有效作用距離。
2.如權(quán)利要求1所述的深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,其特征在于,所述步驟(1)過程如下:
建立目標(biāo)井井筒溫度場模型并劃分網(wǎng)格單元,以油管中心為軸:徑向上,定義ri為網(wǎng)格單元半徑,劃分N個網(wǎng)格單元(i=0,1,2,…,N),其中模型邊界點rN點的選擇應(yīng)滿足該點處溫度始終等于地層溫度;縱向上,設(shè)目的層深度為H,單元深度為ΔH,則縱向上可劃分為M個單元(j=0,1,2…,M),邊界點M為(H/ΔH)取整;定義施工總時間為tT,時間步長為Δt,時間步為n(n=0,1,2…,tT),利用下式計算n時刻單元(i,j)的溫度 最終得到施工結(jié)束時的井底溫度
式中:Q——注酸排量,m3/min;
ρi——網(wǎng)格單元i的密度,kg/m3;
Ci——網(wǎng)格單元i的比熱容,J/(kg·℃);
ΔHj——網(wǎng)格單元j的高度,m;
ri——模型徑向網(wǎng)格半徑,m;
λi——網(wǎng)格單元i的熱傳導(dǎo)系數(shù),W/(m·℃);
Δt——時間步長,min;
tT——施工總時間,min;
n——時間步,min;
——n時刻網(wǎng)格單元(i,j)的溫度,℃。
3.如權(quán)利要求1所述的深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,其特征在于,所述步驟(2)過程如下:
計算施工結(jié)束即n=tT時的裂縫縫長,定義x為縫長方向并劃分網(wǎng)格單元,計算各單元對應(yīng)的平均縫寬w:
式中:L——施工結(jié)束時的縫長,m;
c——酸液濾失系數(shù),m/min0.5;
μf——酸液粘度,Pa·min;
h——裂縫高度,m;
E——儲層巖石楊氏模量,Pa
ε——儲層巖石泊松比,無因次;
x——縫長方向上單元坐標(biāo),m;
w——施工結(jié)束時各單元的平均縫寬,m;
從而得到各單元酸液流速vx:
式中:vl——各單元濾失速度,m/s;
vx——各單元x方向流速,m/s。
4.如權(quán)利要求1所述的深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,其特征在于,所述步驟(3)過程如下:
假設(shè)裂縫入口溫度等于井底溫度,設(shè)ΔrQm(Tw,p)為0,求出各單元不考慮酸巖反應(yīng)熱的酸液溫度與裂縫壁面溫度:
hT=NNuKhf/w
hT(Tw-T)=q(tT)-kRCmΔrQm(Tw,p)
式中:T——各單元酸液平均溫度,℃;
Tw——各單元裂縫壁面溫度,℃;
2
hT——對流換熱系數(shù),J/(m·min·℃);
NNu——努謝爾特數(shù),無因次;
Khf——酸液的導(dǎo)熱系數(shù),J/(m·min·℃);
ρf——儲層巖石密度,kg/m3;
Cf——儲層巖石比熱容,J/(kg·℃);
kR——酸巖反應(yīng)速率常數(shù),(kg/m3)-m·mol/(m2·s);
C——各單元酸液濃度,kg/m3;
ΔrQm(Tw,p)——酸巖摩爾反應(yīng)熱,kJ/mol;
m——反應(yīng)級數(shù),無因次;
q(tT)——施工結(jié)束地層傳熱熱量,kJ;
Mma——地層巖石體積熱容,J/m3·℃;
Khr——地層巖石導(dǎo)熱系數(shù),W/(m·℃);
——孔隙度,%;
通過下式得到各單元酸液濃度與裂縫壁面濃度:
kg(C-Cw)=kRCwm
kg=DeNSh/w
式中:kg——對流傳質(zhì)系數(shù),m/min;
Cw——各單元裂縫壁面酸液濃度,kg/m3;
De——酸液氫離子有效傳質(zhì)系數(shù),m2/s;
Nsh——舍伍德數(shù),無因次;
k0——頻率因子,(kg/m3)-m·mol/(m2·s);
Ea——反應(yīng)活化能,J/mol;
R——氣體常數(shù),8.314J/(mol·℃)。
5.如權(quán)利要求1所述的深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,其特征在于,所述步驟(4)過程如下:
首先計算高溫高壓條件酸巖摩爾反應(yīng)熱ΔrHm(Tw,p):
式中:ΔrHm(Tw)——溫度Tw下的摩爾反應(yīng)熱,kJ/mol;
p——地層壓,MPa;
Cp,m(i,T)——組分i在溫度T下的常壓摩爾熱容,J/(℃·mol);
fg——CO2摩爾分?jǐn)?shù),無因次;
Vco2——CO2摩爾體積,m3/mol;
Vacid——殘酸體積,m3;
Sco2——殘酸中CO2溶解度,m3/m3;
ΔrHm(Tw,p)——溫度Tw壓力p下的摩爾反應(yīng)熱,kJ/mol;
從而得到地層條件下酸巖摩爾反應(yīng)熱ΔrQm(Tw,p):
6.如權(quán)利要求1所述的深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,其特征在于,所述步驟(5)過程如下:
1)通過步驟(4)計算的酸巖摩爾反應(yīng)熱,由步驟(3)計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,并作縫長方向各單元酸液濃度分布曲線;
2)定義酸液初始注入濃度的10%作為殘酸極限濃度,取殘酸極限濃度對應(yīng)的裂縫縫長為酸液有效作用距離。

說明書全文

深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法

技術(shù)領(lǐng)域

[0001] 本發(fā)明涉及石油工程領(lǐng)域,尤其是深層石灰?guī)r儲層酸壓過程中酸液有效作用距離的數(shù)值計算方法。采用數(shù)值計算方法計算高溫條件下的酸液有效作用距離,提高深層石灰?guī)r儲層酸壓設(shè)計的準(zhǔn)確性,實現(xiàn)深層石灰?guī)r儲層高效酸壓改造。

背景技術(shù)

[0002] 川西地區(qū)二疊統(tǒng)深層酸鹽巖儲層展現(xiàn)出良好的勘探潛,已成為四川盆地現(xiàn)階段最現(xiàn)實的接替層系之一。儲層以石灰?guī)r、白巖為主,具有埋藏深(>5000m)、溫度高(>145℃)等特征。該地區(qū)開發(fā)難以自然獲得工業(yè)油氣流,酸壓裂技術(shù)是區(qū)儲層主要的增產(chǎn)改造手段之一(齊寶權(quán),趙佐安,賀洪舉等.應(yīng)用測井儲層品質(zhì)綜合評價指數(shù)來判別儲層的有效性——以四川盆地西部地區(qū)二疊系棲霞組為例[J].天然氣工業(yè),2018,38(02):25-32)。
[0003] 酸化壓裂(簡稱酸壓)是指在高于儲層破裂壓力或天然裂縫的閉合壓力下,將酸液擠入儲層,在儲層中形成裂縫,同時酸液與裂縫壁面巖石發(fā)生化學(xué)反應(yīng),非均勻刻蝕裂縫壁巖石,形成溝槽狀或凹凸不平的刻蝕裂縫,施工結(jié)束后裂縫不完全閉合,最終形成具有一定幾何尺寸和導(dǎo)流能力的酸蝕裂縫,實現(xiàn)油氣井增產(chǎn)。酸液有效作用距離是評價酸壓效果的重要指標(biāo)之一,它主要受溫度、酸液濃度等因素影響。其中,溫度是影響石灰?guī)r儲層酸壓效果的重要因素(李穎川.采油工程[M].石油工業(yè)出版社,2009)。
[0004] 目前常用的酸液有效作用距離計算方法主要關(guān)注殘酸極限濃度的計算(李建輝,慕立俊,趙振峰,等.一種基于殘酸極限評價酸液有效作用距離的方法,CN105044288A[P].2015),很少綜合考慮壓裂過程中酸液的溫度變化對于酸壓改造效果的影響。因此,針對川西地區(qū)儲層溫度高的特點,研究酸壓過程酸液溫度變化規(guī)律,明確溫度對于酸壓效果(酸液有效作用距離)的影響,對于川西地區(qū)儲層酸壓改造設(shè)計至關(guān)重要。

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

[0005] 本發(fā)明的目的在于提供深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,該方法重點研究酸壓過程中酸液受到的熱效應(yīng)(井筒傳熱、酸液流動換熱、酸巖反應(yīng)放熱),采用數(shù)值方法計算考慮以上熱效應(yīng)的酸液溫度與濃度,最終基于酸液濃度確定酸液有效作用距離。該方法原理可靠,有利于計算高溫地層的酸液溫度及濃度,提高酸液有效作用距離計算的準(zhǔn)確性,對深層石灰?guī)r儲層酸壓設(shè)計具有指導(dǎo)意義。
[0006] 為達(dá)到以上技術(shù)目的,本發(fā)明提供以下技術(shù)方案。
[0007] 首先,建立井筒溫度場模型模擬井筒傳熱并計算井底溫度;其次,采用PKN公式計算酸壓裂縫的縫長和縫寬,根據(jù)裂縫縫寬分布求解酸液流速;再者,以井底溫度為初始條件,根據(jù)流速求解酸液溫度和濃度;然后,通過酸液溫度和濃度計算酸巖反應(yīng)熱;最后,計算考慮酸巖反應(yīng)熱后的酸液溫度和濃度,建立酸液濃度與裂縫長度的關(guān)系曲線,確定酸液有效作用距離。
[0008] 深層石灰?guī)r儲層酸壓酸液有效作用距離的數(shù)值計算方法,依次包括以下步驟:
[0009] (1)采用Eickmeier井筒溫度場模型模擬井筒溫度場并計算施工結(jié)束時的井底溫度(王鴻勛,張士誠.力壓裂設(shè)計數(shù)值計算方法[M].北京:石油工業(yè)出版社,1998);
[0010] (2)采用PKN公式計算施工結(jié)束時的裂縫縫長(李穎川.采油工程[M].北京:石油工業(yè)出版社,2001),在縫長方向劃分網(wǎng)格單元并計算各單元的平均縫寬,再計算各單元酸液流速;
[0011] (3)以步驟(1)計算的井底溫度為初始條件,通過步驟(2)計算的酸液流速計算各單元不考慮酸巖反應(yīng)熱的酸液溫度以及酸液濃度;
[0012] (4)通過步驟(3)中計算的酸液溫度和濃度計算各單元的酸巖摩爾反應(yīng)熱;
[0013] (5)將步驟(4)計算的酸巖摩爾反應(yīng)熱代入步驟(3)中,計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,建立酸液濃度分布曲線。定義酸液初始注入濃度的10%作為殘酸極限濃度,根據(jù)分布曲線以殘酸極限濃度對應(yīng)的縫長為酸液有效作用距離。
[0014] 在本發(fā)明中,所述步驟(1)中,采用Eickmeier井筒溫度場模型模擬井筒溫度場并計算施工結(jié)束時的井底溫度,過程如下:
[0015] 采用Eickmeier井筒溫度場模型建模方法建立目標(biāo)井井筒溫度場模型并劃分網(wǎng)格單元。以油管中心為軸:徑向上,定義ri為網(wǎng)格單元半徑,劃分N個網(wǎng)格單元(i=0,1,2,…,N)。其中模型邊界點rN點的選擇應(yīng)滿足該點處溫度始終等于地層溫度;縱向上,設(shè)目的層深度為H,單元深度為ΔH,則縱向上可劃分為M個單元(j=0,1,2…,M),邊界點M為(H/ΔH)取整。定義施工總時間為tT,時間步長為Δt,時間步為n(n=0,1,2…,tT)。結(jié)合目標(biāo)井地質(zhì)工程參數(shù),利用下式計算n時刻單元(i,j)的溫度 最終得到施工結(jié)束時的井底溫度[0016]
[0017] 式中:Q——注酸排量,m3/min;
[0018] ρi——網(wǎng)格單元i的密度,kg/m3;
[0019] Ci——網(wǎng)格單元i的比熱容,J/(kg·℃);
[0020] ΔHj——網(wǎng)格單元j的高度,m;
[0021] ri——模型徑向網(wǎng)格半徑,m;
[0022] λi——網(wǎng)格單元i的熱傳導(dǎo)系數(shù),W/(m·℃);
[0023] Δt——時間步長,min;
[0024] tT——施工總時間,min;
[0025] n——時間步,min;
[0026] ——n時刻網(wǎng)格單元(i,j)的溫度( 為施工結(jié)束時的井底溫度),℃。
[0027] 在本發(fā)明中,所述步驟(2)中,采用PKN公式計算施工結(jié)束時的裂縫縫長,在縫長方向劃分網(wǎng)格單元并計算各單元的平均縫寬,再計算各單元酸液流速,過程如下:
[0028] 1)通過計算施工結(jié)束即n=tT時的裂縫縫長,定義x為縫長方向并劃分網(wǎng)格單元,計算各單元對應(yīng)的平均縫寬w:
[0029]
[0030] 式中:L——施工結(jié)束時的縫長,m;
[0031] c——酸液濾失系數(shù),m/min0.5;
[0032] μf——酸液粘度,Pa·min;
[0033] h——裂縫高度,m;
[0034] E——儲層巖石楊氏模量,Pa
[0035] ε——儲層巖石泊松比,無因次;
[0036] x——縫長方向上單元坐標(biāo),m;
[0037] w——施工結(jié)束時各單元的平均縫寬,m。
[0038] 2)通過下式得到各單元酸液流速vx:
[0039]
[0040] 式中:vl——各單元濾失速度,m/s;
[0041] vx——各單元x方向流速,m/s。
[0042] 公式(3)的推導(dǎo)過程如下:
[0043] 酸液連續(xù)性方程為:
[0044]
[0045] 假設(shè)裂縫內(nèi)酸液流速考慮為平均流速,縫內(nèi)壓力均勻分布,在裂縫軸線附近(即y=0處)的壓力梯度為0,所以有vy|y=0=0。另外,裂縫壁面上酸液流速等于濾失速率,即因此將式(4)由裂縫軸線(y=0)至裂縫壁面(y=w/2)積分,即可推導(dǎo)公式(3)。
[0046]
[0047] 在本發(fā)明中,所述步驟(3)中,以步驟(1)計算的井底溫度為初始條件,通過步驟(2)計算的酸液流速計算各單元不考慮酸巖反應(yīng)熱的酸液溫度以及酸液濃度,過程如下:
[0048] 1)假設(shè)裂縫入口溫度等于步驟(1)計算的井底溫度,縫長方向各單元溫度考慮為平均溫度,將步驟(2)求解的酸液流速代入下式求解各單元酸液溫度與裂縫壁面溫度:
[0049]
[0050] 其中,對流換熱系數(shù)可由努謝爾特準(zhǔn)則求解:
[0051] hT=NNuKhf/w???(7)
[0052] 式中:T——各單元酸液平均溫度,℃;
[0053] Tw——各單元裂縫壁面溫度,℃;
[0054] hT——對流換熱系數(shù),J/(m2·min·℃);
[0055] NNu——努謝爾特數(shù),無因次;
[0056] Khf——酸液的導(dǎo)熱系數(shù),J/(m·min·℃);
[0057] ρf——儲層巖石密度,kg/m3;
[0058] Cf——儲層巖石比熱容,J/(kg·℃)。
[0059] 2)假設(shè)酸液對流換熱熱量hT(Tw-T)主要由地層傳熱熱量以及酸巖反應(yīng)熱熱量提供:
[0060] hT(Tw-T)=q(tT)-kRCmΔrQm(Tw,p)???(8)
[0061] 式中:kR——酸巖反應(yīng)速率常數(shù),(kg/m3)-m·mol/(m2·s);
[0062] C——各單元酸液濃度,kg/m3;
[0063] ΔrQm(Tw,p)——酸巖摩爾反應(yīng)熱,kJ/mol;
[0064] m——反應(yīng)級數(shù),無因次;
[0065] q(tT)——施工結(jié)束地層傳熱熱量,kJ。
[0066] 其中,熱流函數(shù)公式為(Whitsitt,Dysart.The?effect?of?temperature?on?stimulation?design[J].SPE?2497,1970):
[0067]
[0068] 式中:Mma——地層巖石體積熱容,J/m3·℃;
[0069] Khr——地層巖石導(dǎo)熱系數(shù),W/(m·℃);
[0070] ——孔隙度,%。
[0071] 由于酸巖反應(yīng)熱是關(guān)于酸液溫度和濃度的函數(shù),因此在步驟(3)中先不考慮酸巖反應(yīng)熱(設(shè)其值為0),將式(9)、式(8)代入式(6)即可求出各單元不考慮酸巖反應(yīng)熱的酸液溫度。
[0072] 式(6)的推導(dǎo)過程如下:
[0073] 假設(shè)酸液溫度場為穩(wěn)態(tài)。各單元可建立能量守恒方程:
[0074]
[0075] 同公式(5)的處理方法,對公式(10)由裂縫軸線(y=0)至裂縫壁面(y=w/2)積分:
[0076]
[0077] 由于裂縫軸線(y=0)速度梯度為0,因此溫度梯度也為0。各單元溫度考慮為平均溫度,將式(11)處理為常微分方程:
[0078]
[0079] 對流換熱系數(shù)為表征流體與固體表面之間換熱能力的參數(shù),定義為單位面積交換熱量與流體固體溫度差的比值(何潮洪,馮霄.化工原理[M].北京:科學(xué)出版社,2001:190)。因此,定義hT為對流換熱系數(shù):
[0080]
[0081] 需要說明的是,式(7)是基于努謝爾數(shù)定義推導(dǎo)的求解公式,是hT另一種表現(xiàn)形式。將式(3)、式(13)代入式(12),即可推導(dǎo)式(6)。
[0082] 3)假設(shè)各單元酸液濃度為平均濃度,將步驟(2)求解的酸液流速代入下式得到各單元酸液濃度與裂縫壁面濃度:
[0083]
[0084] 式中:kg——對流傳質(zhì)系數(shù),m/min;
[0085] Cw——各單元裂縫壁面酸液濃度,kg/m3;
[0086] 假設(shè)裂縫壁面的酸液與裂縫壁面巖石完全反應(yīng),因此有:
[0087] kg(C-Cw)=kRCwm???(15)
[0088] 式中,由舍伍德數(shù)定義,對流傳質(zhì)系數(shù)可以表示為:
[0089] kg=DeNSh/w???(16)
[0090] 酸巖反應(yīng)速率常數(shù)KR由阿倫尼烏斯公式求解(李穎川.采油工程[M].北京:石油工業(yè)出版社,2001):
[0091]
[0092] 式中:De——酸液氫離子有效傳質(zhì)系數(shù),m2/s;
[0093] Nsh——舍伍德數(shù),無因次;
[0094] k0——頻率因子,(kg/m3)-m·mol/(m2·s);
[0095] Ea——反應(yīng)活化能,J/mol;
[0096] R——氣體常數(shù),8.314J/(mol·℃)。
[0097] 因此,將步驟2)求解的各單元酸液溫度代入式(17)求解得到KR,而后聯(lián)立求解式(14)、式(15)即可求得各單元酸液濃度。
[0098] 式(14)的推導(dǎo)過程如下:
[0099] 各單元建立酸液濃度物質(zhì)平衡方程:
[0100]
[0101] 同公式(5)的處理方法,對公式(18)由裂縫軸線(y=0)至裂縫壁面(y=w/2)積分:
[0102]
[0103] 由于裂縫軸線(y=0)速度梯度為0,因此濃度梯度也為0。各單元濃度考慮為平均濃度,將式(19)處理為常微分方程:
[0104]
[0105] 由對流傳質(zhì)系數(shù)定義(劉瑾,王豪,王治紅.傳熱與傳質(zhì)[M].西南石油大學(xué),2007):
[0106]
[0107] 需要說明的是,式(16)是基于舍伍德數(shù)定義推導(dǎo)的求解公式,是kg另一種表現(xiàn)形式。將式(3)、式(21)代入式(20),即可推導(dǎo)式(14)。綜上,由步驟(3)可求得各單元不考慮酸巖反應(yīng)熱的酸液溫度和濃度。
[0108] 在本發(fā)明中,所述步驟(4)中,結(jié)合步驟(3)中計算的酸液溫度和濃度計算各單元的酸巖摩爾反應(yīng)熱,過程如下:
[0109] 1)考慮目標(biāo)工區(qū)儲層巖性為石灰?guī)r,通過查閱物理化學(xué)手冊(孫艷輝,何廣平,國正等.物理化學(xué)實用手冊[M].北京:化學(xué)工業(yè)出版社.2016)獲得反應(yīng)體系中各組分的標(biāo)準(zhǔn)摩爾反應(yīng)熱,并計算反應(yīng)體系的摩爾反應(yīng)熱
[0110]
[0111] 式中:υi——反應(yīng)體系中i組分的化學(xué)計量系數(shù),無因次;
[0112] ——i組分的標(biāo)準(zhǔn)摩爾反應(yīng)熱,kJ/mol;
[0113] ——反應(yīng)體系標(biāo)準(zhǔn)摩爾反應(yīng)熱,kJ/mol。
[0114] 2)結(jié)合步驟1)的計算結(jié)果,計算高溫高壓條件酸巖摩爾反應(yīng)熱ΔrHm(Tw,p)(Guo?J,Liu?H,Zhu?Y,et?al.Effects?of?acid–rock?reaction?heat?on?fluid?temperature?profile?in?fracture?during?acid?fracturing?in?carbonate?reservoirs[J].Journal?of?petroleum?science&engineering,2014,122:31-37):
[0115]
[0116] 其中,fg為游離態(tài)CO2摩爾分?jǐn)?shù):
[0117]
[0118] 式中:ΔrHm(Tw)——溫度Tw下的摩爾反應(yīng)熱,kJ/mol;
[0119] p——地層壓力,MPa;
[0120] Cp,m(i,T)——組分i在溫度T下的常壓摩爾熱容,J/(℃·mol);
[0121] fg——CO2摩爾分?jǐn)?shù),無因次;
[0122] Vco2——CO2摩爾體積,m3/mol;
[0123] Vacid——殘酸體積,m3;
[0124] Sco2——殘酸中CO2溶解度,m3/m3。
[0125] ΔrHm(Tw,p)——溫度Tw壓力p下的摩爾反應(yīng)熱,kJ/mol。
[0126] 3)酸巖反應(yīng)過程中,反應(yīng)體系對CO2做體積功。因此,綜合考慮溫度、壓力和CO2的體積功的影響,得到地層條件下酸巖摩爾反應(yīng)熱ΔrQm(Tw,p):
[0127]
[0128] 在本發(fā)明中,所述步驟(5)中,將步驟(4)計算的酸巖摩爾反應(yīng)熱代入步驟(3)中,計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,建立酸液濃度分布曲線。定義酸液初始注入濃度的10%作為殘酸極限濃度,根據(jù)分布曲線以殘酸極限濃度對應(yīng)的縫長為酸液有效作用距離,過程如下:
[0129] 1)將步驟(4)計算的酸巖摩爾反應(yīng)熱代入公式(8)中計算各單元酸巖反應(yīng)熱,而后由步驟(3)計算考慮酸巖反應(yīng)熱的各單元酸液溫度以及酸液濃度,并作縫長方向各單元酸液濃度分布曲線;
[0130] 2)定義酸液初始注入濃度的10%作為殘酸極限濃度,參考酸液濃度分布曲線,取殘酸極限濃度對應(yīng)的裂縫縫長為酸液有效作用距離。
[0131] 與現(xiàn)有技術(shù)相比,本發(fā)明提供了一種針對深層石灰?guī)r儲層的酸液有效作用距離的數(shù)值計算方法。本方法在計算過程中耦合了井筒傳熱、酸液對流換熱、酸巖反應(yīng)熱三種熱效應(yīng):①首先通過井筒溫度場模型計算施工結(jié)束井底溫度;②利用PKN公式計算的縫寬求解縫長方向各單元酸液流速,以井底溫度為初始條件,通過酸液流速求解各單元酸液溫度以及濃度,并實現(xiàn)井筒傳熱與對流換熱的耦合;③通過各單元酸液溫度、濃度求解酸巖反應(yīng)熱,耦合酸巖反應(yīng)熱后計算各單元酸液濃度得到酸液濃度分布曲線,確定酸液有效作用距離。本方法克服了現(xiàn)有方法未能綜合考慮壓裂過程中多種熱效應(yīng)對于酸壓改造效果影響的局限。
附圖說明
[0132] 圖1是本發(fā)明中井筒溫度場模型網(wǎng)格劃分示意圖。
[0133] 圖2是本發(fā)明中井底溫度變化模擬結(jié)果圖。
[0134] 圖3是本發(fā)明中酸巖反應(yīng)熱計算結(jié)果。
[0135] 圖4是本發(fā)明中酸液有效作用距離計算結(jié)果圖。

具體實施方式

[0136] 下面結(jié)合附圖和應(yīng)用實例進(jìn)一步說明本發(fā)明,旨在對本發(fā)明做示意性說明和解釋,并不限定本發(fā)明的范圍。具體如下:
[0137] 1.基于目標(biāo)工區(qū)某探井X井酸壓設(shè)計基礎(chǔ)參數(shù)(表1),建立X井井筒溫度場物理模型并劃分網(wǎng)格(見圖1),采用步驟(1)建立的井筒溫度場模型計算該井施工過程中動態(tài)井底溫度,得到施工結(jié)束時的井底溫度。對比實測井底溫度曲線,可以證明該模型可靠(見圖2)。
[0138] 表1?X井酸壓設(shè)計部分基礎(chǔ)參數(shù)
[0139]
[0140] 2.以施工結(jié)束時的井底溫度為初始條件,采用數(shù)值方法計算步驟(2)和步驟(3)中的裂縫酸液溫度和濃度方程,得到施工結(jié)束時各單元酸液溫度和濃度。
[0141] 3.基于具體實施方式中的步驟2計算得到的各單元的酸液溫度和濃度計算各單元的酸巖反應(yīng)熱,如圖3所示。
[0142] 4.將具體實施方式中的步驟3計算的各單元酸巖反應(yīng)熱代入步驟2中,計算考慮酸巖反應(yīng)熱的酸液溫度以及酸液濃度,得到縫長方向各單元酸液濃度分布曲線,取殘酸極限濃度對應(yīng)的裂縫縫長為酸液有效作用距離,如圖4所示。
高效檢索全球?qū)@?/div>

專利匯是專利免費檢索,專利查詢,專利分析-國家發(fā)明專利查詢檢索分析平臺,是提供專利分析,專利查詢,專利檢索等數(shù)據(jù)服務(wù)功能的知識產(chǎn)權(quán)數(shù)據(jù)服務(wù)商。

我們的產(chǎn)品包含105個國家的1.26億組數(shù)據(jù),免費查、免費專利分析。

申請試用

分析報告

專利匯分析報告產(chǎn)品可以對行業(yè)情報數(shù)據(jù)進(jìn)行梳理分析,涉及維度包括行業(yè)專利基本狀況分析、地域分析、技術(shù)分析、發(fā)明人分析、申請人分析、專利權(quán)人分析、失效分析、核心專利分析、法律分析、研發(fā)重點分析、企業(yè)專利處境分析、技術(shù)處境分析、專利壽命分析、企業(yè)定位分析、引證分析等超過60個分析角度,系統(tǒng)通過AI智能系統(tǒng)對圖表進(jìn)行解讀,只需1分鐘,一鍵生成行業(yè)專利分析報告。

申請試用

QQ群二維碼
意見反饋