3. 流变模型与参数解读
固井水泥浆的流变性能决定了它在环空中的流动压降、顶替效率以及对井壁的悬浮能力。现场测得的旋转黏度计读数(如600、300、200、100、6、3 r/min)本身只是一组转速—刻度盘读数,必须借助流变模型把它们转换成可计算的参数,才能用于压降和顶替设计。本节围绕三种最常用的本构模型展开:宾汉塑性模型、幂律模型和赫巴模型,重点讲清每个模型的方程形式、参数物理意义,以及如何从实测读数反算参数。
一、三种本构方程的形式与物理意义
宾汉塑性模型是最经典的描述,其本构方程为 τ = τ₀ + μ_p·γ̇,其中 τ 为剪切应力(Pa 或 lbf/100ft²),γ̇ 为剪切速率(s⁻¹),τ₀ 为动切力(Yield Point,YP),μ_p 为塑性粘度(PV)。它的物理图像是:浆体存在一个屈服门槛,只有当外加剪切应力超过 τ₀ 后才开始流动,流动后应力与剪切速率呈线性关系。塑性粘度反映的是浆体内部固相颗粒与液相之间摩擦造成的黏性阻力,主要受固相含量、颗粒级配和温度影响;动切力则反映浆体在流动状态下的结构强度,与黏土颗粒表面的电化学作用、絮凝结构有关。宾汉模型对常规密度水泥浆在中高剪切速率段拟合较好,但在低剪切速率段往往高估剪切应力。
幂律模型的本构方程为 τ = K·γ̇ⁿ,其中 K 为稠度系数(Pa·sⁿ),n 为流性指数(无量纲)。幂律模型没有屈服应力项,认为浆体从极低剪切速率起就服从幂函数关系。n 反映浆体的剪切稀释能力:n < 1 时为假塑性流体,n 越小剪切稀释越强;n = 1 时退化为牛顿流体;n > 1 时为膨胀性流体。K 则综合反映浆体的稠度水平,K 越大浆体越稠。幂律模型在描述低剪切速率下的剪切稀释行为方面优于宾汉模型,适合低密度、高水灰比或聚合物类水泥浆体系。
赫巴模型(Herschel-Bulkley)是本构方程为 τ = τ₀ + K·γ̇ⁿ 的三参数模型,它同时包含屈服应力 τ₀、稠度系数 K 和流性指数 n。赫巴模型可以看作宾汉模型与幂律模型的结合:当 n = 1 时退化为宾汉模型,当 τ₀ = 0 时退化为幂律模型。由于它既能描述屈服门槛,又能描述剪切稀释,对水泥浆在全剪切速率范围内的拟合精度通常最高,尤其适用于含降失水剂、分散剂等外加剂的复杂体系。
二、从黏度计读数反算模型参数
现场最常用的是范氏旋转黏度计,其读数 θ(刻度盘偏转角)与剪切应力的换算关系为 τ = 0.511·θ(单位 Pa),剪切速率与转速的关系为 γ̇ = 1.703·N(N 为 r/min)。基于 600 r/min 和 300 r/min 两个读数,宾汉模型参数可直接得到:
PV = θ₆₀₀ − θ₃₀₀ (单位 cP,即 mPa·s)
YP = θ₃₀₀ − PV (单位 lbf/100ft²)
= 0.511 × (θ₃₀₀ − PV) (换算为 Pa)
AV = θ₆₀₀ / 2 (单位 cP)
举例说明(以下为教学示例数据,非实测结果):某常规密度水泥浆测得 θ₆₀₀ = 120,θ₃₀₀ = 70。则 PV = 120 − 70 = 50 cP,YP = 70 − 50 = 20 lbf/100ft² ≈ 10.2 Pa,AV = 60 cP。可以看出 AV 介于 PV 与 PV+YP 之间,它把屈服应力的贡献也折算进了“等效黏度”,因此 AV 只能用于横向比较,不能直接代入宾汉方程计算压降。
幂律模型的 n 和 K 通常用 600/300 两读数计算:
n = 3.32 × lg(θ₆₀₀ / θ₃₀₀)
K = 0.511 × θ₃₀₀ / (511)^n (单位 Pa·sⁿ,511 为 300 r/min 对应剪切速率)
代入上例:n = 3.32 × lg(120/70) = 3.32 × 0.234 ≈ 0.78,说明该浆体具有明显剪切稀释性;K = 0.511 × 70 / 511^0.78 ≈ 0.28 Pa·sⁿ。n 值偏低时,低剪切速率下黏度上升快,有利于停泵时悬浮加重剂,但循环压耗也会增大,需要结合环空压耗核算。
赫巴模型为三参数,需要至少三组读数(常用 600、300、6 或 600、300、3 r/min)进行非线性回归。工程上可先由低转速读数估算 τ₀(如取 θ₃ 换算),再用 600/300 读数回归 K 与 n,最后用中间转速点校核拟合偏差。若三个模型对同一组读数的拟合残差差异明显,应优先选用残差最小的模型;对含外加剂较多的体系,赫巴模型通常更贴合。
三、参数在固井设计中的判读要点
解读参数时不能只看单个数值,而要结合工况判断。塑性粘度偏高(如超过 80 cP)往往意味着固相含量过高或分散不良,会推高循环压耗,在窄环空井段可能压漏地层;动切力偏低则悬浮能力不足,停泵后加重剂易下沉,影响浆柱稳定性。流性指数 n 过低(如小于 0.5)时,浆体在低剪切速率下黏度急剧上升,虽然悬浮性好,但顶替时界面推进不均,容易形成滞留区。实际设计中常以“PV 适中、YP 满足悬浮、n 在 0.6~0.9 之间”作为初步判据,再通过环空压耗和顶替效率模拟最终确认。
此外,温度对参数影响显著。井底循环温度下水泥浆的 PV 通常下降、n 略升,而地面测得的数据偏保守。工程上可依据 API 推荐的温度校正关系或实测高温流变曲线进行修正,避免用地面参数直接计算井底压耗而造成误判。掌握三种模型的方程形式与参数含义,并能从黏度计读数正确反算,是后续密度—流变配方联算的基础。