基于Boltzmann方法粗糙岩石裂隙渗流-传热-传质耦合机制研究

    Investigation of coupled seepage-heat-mass processes in rough rock fracture based on Lattice Boltzmann Method

    • 摘要: 为研究粗糙岩石裂隙渗流-传热-传质耦合作用机制,基于格子Boltzmann方法,采用三个分布函数分别模拟速度场、温度场和浓度场的演化过程,考虑流体温度对其运动黏度、热扩散系数以及溶质扩散系数的影响,建立了模拟三维粗糙岩石裂隙渗流-传热-传质耦合过程的数值计算模型,并通过两个经典算例证明其准确性。最后,讨论了岩石温度、压力梯度和裂隙粗糙度等因素对裂隙内溶质运移过程的影响,并基于参数敏感性综合评价方法,探讨了裂隙开度、压力梯度、岩石温度、裂隙粗糙度、溶质初始扩散系数5个影响因素的综合敏感度。研究结果表明:随着岩石温度的升高,裂隙内流体的运动黏度降低,溶质扩散系数变大,使得岩石裂隙溶质运移显著加快。当岩石温度从20℃升高至70℃时,流体运动黏度降低至初始值的48%,溶质扩散系数增大为初始值的5.17倍。当渗流流速增大时,在流速和温度综合影响下加快了溶质运移,但流速的增大在一定程度上削弱了温度对溶质运移的加速作用。随着裂隙粗糙程度的增大,裂隙内流体流速变缓,裂隙内温度逐渐升高,流速降低对溶质运移的影响明显高于温度升高诱导的溶质扩散作用,最终呈现裂隙粗糙度越大,流体流速越低,溶质运移越慢。粗糙岩石裂隙溶质运移过程中,裂隙开度表现最为敏感,而溶质初始扩散系数则敏感性最低。

       

      Abstract: To elucidate the coupled mechanisms of seepage, heat transfer, and mass transfer in rough rock fractures, a three-dimensional numerical model was developed based on the lattice Boltzmann method (LBM). The model employs three distinct distribution functions to capture the evolution of the velocity, temperature, and concentration fields, explicitly accounting for the dependence of fluid kinematic viscosity, thermal diffusivity, and solute diffusivity on temperature. Model accuracy was verified against two benchmark cases. The model was subsequently applied to investigate the effects of rock temperature, pressure gradient, and fracture roughness on solute transport within fractures. A comprehensive sensitivity analysis was performed to evaluate the combined influence of five key parameters: fracture aperture, pressure gradient, rock temperature, fracture roughness, and initial solute diffusion coefficient. The results reveal that increasing rock temperature reduces the kinematic viscosity of the fluid while enhancing the solute diffusion coefficient, thereby markedly accelerating solute transport. Specifically, when rock temperature rises from 20°C to 70°C, kinematic viscosity decreases to 48% of its initial value, whereas the solute diffusion coefficient increases by a factor of 5.17 relative to its initial value. Higher seepage velocity further promotes solute transport through the combined effects of advection and thermal enhancement; however, increased flow velocity partially offsets the temperature-induced acceleration. Increasing fracture roughness leads to reduced fluid velocity and elevated temperature within the fracture. The deceleration of flow exerts a substantially greater inhibitory effect on solute transport than the enhancement of diffusion arising from temperature elevation. Consequently, greater fracture roughness results in slower solute transport. Among all investigated factors, fracture aperture exhibits the highest sensitivity, whereas the initial solute diffusion coefficient shows the lowest sensitivity.

       

    /

    返回文章
    返回