技術(shù)特征：

1.二氧化碳驅(qū)替前緣的試井確定方法，其特征在于，該二氧化碳驅(qū)替前緣的試井確定方法包括：

步驟1，建試井解釋模型，包括油藏模型和滲流數(shù)學(xué)模型；

步驟2，通過(guò)壓力方程線性化，數(shù)值求解試井解釋模型；

步驟3，根據(jù)基本的油藏物性測(cè)試資料，通過(guò)基礎(chǔ)試井解釋參數(shù)，運(yùn)用模擬方法進(jìn)行求解，獲得壓力降落/恢復(fù)試井過(guò)程中的壓力數(shù)據(jù)，并進(jìn)行處理，得到壓力及壓力半對(duì)數(shù)或雙對(duì)數(shù)曲線；

步驟4，利用得到的結(jié)果，擬合實(shí)際壓力及壓力半對(duì)數(shù)、雙對(duì)數(shù)曲線，并計(jì)算目標(biāo)函數(shù)值；

步驟5，根據(jù)步驟4獲得的參數(shù)和求解二氧化碳驅(qū)替前緣關(guān)系式，確定二氧化碳驅(qū)替前緣。

2.根據(jù)權(quán)利要求1所述的二氧化碳驅(qū)替前緣的試井確定方法，其特征在于，在步驟1中，建立的油藏模型和滲流數(shù)學(xué)模型為：

$\frac{\∂}{\∂t}\[\mathrm{\φ}(x_u{\mathrm{\ρ}}_o{S}_o+y_u{\mathrm{\ρ}}_g{S}_g)\]+\▿\[x_u{\mathrm{\ρ}}_o\stackrel{\‾}{V_o}+y_u{\mathrm{\ρ}}_g\stackrel{\‾}{V_g}\]-\stackrel{\‾}{q_u}=0---\left(1\right)$

$\stackrel{\‾}{V_o}=\left\{\begin{array}{c}-\mathrm{\α}\frac{{\mathrm{kk}}_{ro}}{{\mathrm{\μ}}_o}(\▿p_o-\stackrel{\‾}{{\mathrm{\ρ}}_o}g\▿D-G_o),|\▿p_o-\stackrel{\‾}{{\mathrm{\ρ}}_{o}g}\▿D|\>G_o\\ 0,|\▿p_o-\stackrel{\‾}{{\mathrm{\ρ}}_o}g\▿D|\≤G_o\end{array}\right.---\left(2\right)$

$\stackrel{\‾}{V_g}=\left\{\begin{array}{c}-\mathrm{\α}\frac{{\mathrm{kk}}_{rg}}{{\mathrm{\μ}}_g}(\▿p_g-\stackrel{\‾}{{\mathrm{\ρ}}_g}g\▿D-G_g),|\▿p_g-\stackrel{\‾}{{\mathrm{\ρ}}_g}g\▿D|\>G_g\\ 0,|\▿p_g-\stackrel{\‾}{{\mathrm{\ρ}}_g}g\▿D|\≤G_g\end{array}\right.---\left(3\right)$

$k=k^{*}{e}^{-\mathrm{\ϵ}(p_i-p)}---\left(4\right)$

其中，φ孔隙度；x_u為液相(油相)中組分u摩爾分?jǐn)?shù)；y_u為氣相中組分u摩爾分?jǐn)?shù)；z_u為組分u總摩爾分?jǐn)?shù)；s_o為含油飽和度；s_g為含氣飽和度；為u組分的單位體積摩爾質(zhì)量流量；p_l為相壓力；p_i為原始地層壓力；k為地層絕對(duì)滲透率；k_rl為相對(duì)滲透率；k^*為原始滲透率；μ_l為粘度；ε為壓敏系數(shù)；G_l為各相的啟動(dòng)壓力梯度；為質(zhì)量密度；ρ為摩爾密度；下標(biāo)o，g分別代表為油、氣相，且l＝o,g；α為單位換算系數(shù)，當(dāng)上述變量均采用水力學(xué)單位時(shí)，即達(dá)西單位，α＝1。為重力在x、y和z方向上的梯度。

3.根據(jù)權(quán)利要求1所述的二氧化碳驅(qū)替前緣的試井確定方法，其特征在于，在步驟2中，數(shù)值求解試井解釋模型的公式為：

$\begin{array}{c}\mathrm{\Δ}\[{T_o}^n\left({{\mathrm{\Δp}}_o}^{n+1}-{\stackrel{\‾}{{\mathrm{\ρ}}_o}}^n{g}^n{\mathrm{\ΔD}}^n-{G_o}^n\right)+{T_g}^n\left({\mathrm{\Δp}}_o+{\mathrm{\Δp}}_{cog}-{\stackrel{\‾}{{\mathrm{\ρ}}_g}}^n{g}^n{\mathrm{\ΔD}}^n-{G_g}^n\right)\\ +q_o^n+q_g^n=\frac{V}{\mathrm{\Δ}t}\[{\mathrm{\φ}}^{n+1}{\left({\mathrm{\ρ}}_o{S}_o+{\mathrm{\ρ}}_g{S}_g\right)}^{n+1}-{\mathrm{\φ}}^n{\left({\mathrm{\ρ}}_o{S}_o+{\mathrm{\ρ}}_g{S}_g\right)}^n\]\end{array}---\left(5\right)$

式中，Δp_oⁿ⁺¹為n+1時(shí)刻油相壓力差，分別為油、氣質(zhì)量密度；gⁿ為n時(shí)刻重力常數(shù)，ΔDⁿ為n時(shí)刻高度差；G_lⁿ為各相的啟動(dòng)壓力梯度，l＝o,g；S_o,S_g為油、氣飽和度；為油組分及氣組分的單位體積摩爾質(zhì)量流量；Δp_cog為油氣毛管力；V為單元體積；T_oⁿ為差分表達(dá)式中與時(shí)間無(wú)關(guān)的變量；Δt為時(shí)間間隔；Φ為孔隙度。

4.根據(jù)權(quán)利要求1所述的二氧化碳驅(qū)替前緣的試井確定方法，其特征在于，在步驟4中，計(jì)算目標(biāo)函數(shù)值的公式為：

$D\left(\stackrel{\→}{x}\right)=\underset{i=1}{\overset{n}{\Σ}}{\[p\left(\stackrel{\→}{x},t_{Di}\right)-p\left(t_{Di}\right)\]}^2+\underset{i=1}{\overset{n}{\Σ}}{\[dp/d\ln t\left(\stackrel{\→}{x},t_{Di}\right)-dp/d\ln t\left(t_{Di}\right)\]}^2---\left(6\right)$

為最小二乘值，p(t_Di)分別為t_Di時(shí)刻沿搜索方向的壓力和初始?jí)毫Γ?img id="icf0015" file="FDA0000856041180000027.GIF" wi="189" he="85" img-content="drawing" img-format="GIF" orientation="portrait" inline="no" />t(t_Di)分別為沿搜索方向的時(shí)間步長(zhǎng)和初始時(shí)間步長(zhǎng)；p為壓力；n為曲線上的離散點(diǎn)數(shù)目。

用SPSA算法在變量定義域中進(jìn)行搜索，改變相關(guān)參數(shù)，包括二氧化碳區(qū)半徑R₁、表皮系數(shù)S、井筒存儲(chǔ)系數(shù)C、流度比M、導(dǎo)壓系數(shù)比η、滲透率k及壓縮因子，進(jìn)行模擬計(jì)算，使得目標(biāo)函數(shù)值最小，得到最終的擬合參數(shù)解釋結(jié)果。

5.根據(jù)權(quán)利要求1所述的二氧化碳驅(qū)替前緣的試井確定方法，其特征在于，在步驟5中，根據(jù)壓力降落測(cè)試壓力及壓力半對(duì)數(shù)曲線，由圖解法得到波及區(qū)徑向流直線段斜率m_c、波及區(qū)徑向流結(jié)束時(shí)刻t_end，同時(shí)結(jié)合試井解釋結(jié)果，構(gòu)建關(guān)系式，求取二氧化碳驅(qū)替前緣：

$X=\frac{XY+Z-1}{Z}---\left(7\right)$

X＝(R₂/R₁)² (8)

$XY=7\×{10}^5\frac{{\mathrm{qB}}_g{\mathrm{\μ}}_g}{{\mathrm{\πkhm}}_c{t}_{end}}---\left(9\right)$

Z＝η/M (10)

其中，R₂為二氧化碳驅(qū)替前緣半徑，m；R₁為純二氧化碳區(qū)半徑，m；q為二氧化碳注入量，m³/d；B_g為二氧化碳體積系數(shù)；μ_g為二氧化碳粘度，mPa·s；k為試井解釋得到的油藏滲透率，10^-3μm²；h為儲(chǔ)層厚度，m；η為導(dǎo)壓系數(shù)比，M為流度比；X為過(guò)渡區(qū)與CO₂區(qū)半徑比值的平方，Y為復(fù)合參數(shù)表達(dá)式，Z為儲(chǔ)容系數(shù)。