[发明专利]宾汉姆流体在多孔介质中的有效渗透率的预测方法

权利要求书

1.一种宾汉姆流体在多孔介质中的有效渗透率的预测方法，所述宾汉姆流体的本构方程为：

τ=τ0+μγ·;]]>

其中，式中τ为切应力，τ₀为所述宾汉姆流体的屈服应力，μ为所述宾汉姆流体的粘度，是剪切速率；

其特征在于，所述预测方法包括如下步骤：

步骤a，提供所述多孔介质的样本，测量所述多孔介质的孔隙度（φ）和颗粒半径（R）；

步骤b，根据所述多孔介质的孔隙度（φ）和颗粒半径（R）计算所述多孔介质在分形模型中的结构参数，所述结构参数包括最小孔隙半径（r_min）、最大孔隙半径（r_max）、毛细管的直线长度（L₀）、流体路径的弯曲度（Γ）、毛细管的迂曲度分形维数（D_T）、孔隙分形维数（D_f），其分别由如下公式计算得出：

rmax=R4[2(1-0342φ1-φ-1)+2π(1-0.342φ)3(1-φ)-2];]]>

rminrmax=2241-φ1-0.342φ;]]>

L0=R2π(1-0.342φ)3(1-φ);]]>

Γ=12[1+121-φ+1-φ(11-φ-1)2+141-1-φ];]]>

DT=1+lnΓlnL02rav;]]>

rav=DfrminDf-1[1-(rminrmax)Df-1];]]>

Df=2-lnφlnrminrmax;]]>

步骤c，测量所述宾汉姆流体的特性参数，所述特性参数包括所述宾汉姆流体的屈服应力（τ₀）、所述宾汉姆流体的粘度（μ）、所述宾汉姆流体的表面张力（T）、所述宾汉姆流体与所述多孔介质的孔壁之间的接触角（θ）；

步骤d，根据所述宾汉姆流体的表面张力（T）、所述宾汉姆流体与所述多孔介质的孔壁之间的接触角（θ）、所述多孔介质的孔隙度（φ）按照如下公式计算所述宾汉姆流体在所述多孔介质中的平均毛细压差

Δpc‾=2TcosθDf(1-φ)rminφ(1+Df);]]>

步骤e，根据所述多孔介质的结构参数、所述宾汉姆流体的特性参数及所述宾汉姆流体在所述多孔介质中的平均毛细压差及动力压差（Δp_m）按照如下公式计算所述宾汉姆流体的有效渗透率（K_e）：

K_e＝K₁+K₂-K₃；

其中，K1=rmax1+DTφ(2-Df)τΔpm24-DTL0DT-1(3+DT-Df)(1-φ)(τ-τ0)(Δpm+Δpc‾);]]>

K2=rmaxDT(2-Df)τcosθT23-DTL0DT-1(2+DT-Df)(τ-τ0)(Δpm+Δpc‾);]]>

K3=rmaxL0τφ(2-Df)τ04(3-Df)(1-φ)(τ-τ0)(Δpm+Δp‾c).]]>