3乳状液油水分层预测模型

在高含水率下,乳状液油水分层过程受多个因素的影响,包括液滴聚并、界面不稳定性、黏度变化等。传统Stokes定律主要适用于低含水率乳状液,但在高含水率下,液滴运动行为复杂,需对Stokes模型和经验模型进行修正。

高含水率下乳状液油水分层的影响机制:当含水率高时,水滴间的相互作用强,导致小粒径液滴容易发生聚并,形成较大粒径液滴,沉降速率与粒径的二次方成正比(Stokes沉降理论),因此会加速油水分层;当含水率增加时,油水两相界面张力下降,导致破乳剂的破乳效果减弱,相界面趋于不稳定,界面不稳定性的加剧促使水滴更快聚并,进而加快沉降速率;当含水率增加时,乳状液从W/O型乳状液向O/W型乳状液转变,该转变直接影响液滴的沉降行为;油相黏度随含水率的增加呈非线性变化,且高含水率下分层过程的流动形态可能从Stokes层流转变为湍流，因此需对流动阻力系数进行修正以适应不同流动状态。

基于上述高含水率对乳状液油水分层的影响,引入液滴聚并效应修正项fc、界面不稳定修正项fi、含水率W、经验常数k4、比例系数A,对Stokes模型和经验模型进行修正,模型中界面与含水率修正项的设置与已有稳定性表征方法的结论相符。

1)液滴聚并效应修正。聚并效应导致有效液滴粒径增加,可引入修正因子de,其表达式为:

de=d0(1+k1W)

式中:d0为初始液滴粒径,由实验测定,μm;W取值为0~60%;k1为经验常数,通过实验数据拟合得到。

2)针对高含水率下界面不稳定性对沉降速率的强化作用,引入界面不稳定修正项fi,其表达式为:

fi=1+k2W

式中:k2为经验常数,通过实验数据拟合得到。

3)高含水率下乳状液黏度呈非线性变化,引入有效黏度ηeff,其表达式为:

ηeff=η0(1+k3W)

式中:η0为初始黏度,mPa·s;k3为经验常数,通过实验数据拟合得到,用于表征含水率对乳状液黏度的影响。

综合上述修正项,最终建立高含水率下的析水速率修正模型:

v=A·{[g(ρd-ρc)d0²(1+k1W)²(1+k2W)]/[18η0(1+k3W)]·(1-e^(-k4W))}

式中:v为油水分层速率,m/s;ρd为水相密度,kg/m³;ρc为油相密度,kg/m³;g为重力加速度,m/s²。

采用非线性最小二乘法对40组实验数据进行拟合,得到最优参数:A=0.92,k1=0.15,k2=0.08,k3=0.12,模型的拟合决定系数R²=0.956>0.900,表明修正模型的预测值与实验数据吻合良好。为进一步验证模型的有效性,选取乳状液含水率为40%~60%的10组独立实验数据进行了验证,结果如图8所示。由图8可知,预测值与实测值的平均相对误差为8.7%,最大误差不超过12.3%,总体预测误差可控制在10.0%以内,表明该修正模型可用于高含水率下的油水分层速率计算。

目前所建模型是假设乳状液粒径已知,且是基于静态条件下的实验数据建立的。然而,在实际生产过程中乳状液的粒径分布可能受温度、压力、CO2饱和浓度等因素的影响发生变化。此外,本模型未考虑流动条件下的影响因素。为了使模型更广泛适用于动态流动条件,后续研究将进一步考虑流动阻力系数的修正和粒径变化对模型预测结果的影响。

4结论

1)在低搅拌转速(400~500r/min)下，剪切力较弱，液滴破碎受限，液滴粒径较大，易发生聚并和沉降，乳状液稳定性较差;在中等搅拌转速(600~700r/min)下，剪切力可满足液滴均匀分散需求，不仅可有效降低液滴粒径，还可增强界面膜稳定性，乳化效果最佳;在高搅拌转速(800r/min)下,剪切力虽可进一步降低液滴粒径，但过高的剪切力会破坏液滴界面，诱发液滴聚并，导致乳状液稳定性下降。经综合分析,600~700r/min的搅拌转速为优化搅拌转速区间，在此区间可形成粒径较小且稳定性优良的乳状液。

2)当含水率升高至40%~50%时，液滴间相互作用显著增强，导致液滴沉降行为偏离传统Stokes定律，需采用修正的经验模型对其沉降行为进行描述。模型预测结果表明，高含水率下析水速率整体较低，但液滴聚并和界面不稳定性等因素会导致乳状液局部析水现象趋缓，有利于维持短期内的油水分层稳定性。

3)CO2过饱和处理可显著降低油水界面张力，优化乳化剂在界面的吸附，从而使液滴形态趋于规则、粒径减小，提升乳状液的稳定性，但过高的CO2浓度可能诱发界面失稳，加剧液滴聚并。

4)通过对低含水率的Stokes方程和高含水率的经验模型进行修正，拓展了模型的含水率适用范围，使其可对低、中、高含水率乳状液的析水速率进行预测。验证结果表明，该模型的总体预测误差可控制在10%以内，能够有效表征乳状液在静态条件下的稳定性。

5)当前模型存在两项假设:一是乳状液液滴的粒径已知，二是模型建立基于静态实验数据。在实际应用中，随流动条件与温度的变化，模型需进一步修正流动阻力系数以匹配实际流态，并考虑粒径的动态变化。后续研究将围绕流动条件下的模型适用性、流动阻力系数的修正、粒径变化的影响，对模型进行进一步的优化和验证。