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

    Investigation of coupled seepage - heat transfer - mass transfer mechanisms in rough rock fractures based on lattice Boltzmann method

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

       

      Abstract: To investigate the coupled mechanisms of seepage, heat transfer, and mass transfer in a rough rock fracture, a three-dimensional numerical model is proposed based on the lattice Boltzmann method (LBM). The model employs three distribution functions to capture the evolution of the velocity, velocity temperature, and concentration fields, respectively accounting for the dependence of fluid kinematic viscosity, thermal diffusivity, and solute diffusivity on temperature. Model accuracy is verified against two benchmark cases. The model is subsequently applied to investigate the effects of rock temperature, pressure gradient, and fracture roughness on solute transport within fractures. A comprehensive sensitivity evaluation analysis is performed to evaluate the overall sensitivity of five key influencing factors: 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 significantly accelerating solute transport in the rock fracture. Specifically, when rock temperature rises from 20℃ to 70℃, kinematic viscosity decreases to 48% of its initial value, while the solute diffusion coefficient increases to 5.17 times its initial value. Higher seepage velocity further promotes solute transport through the coupling effect between flow velocity and temperature. However, increased flow velocity partially offsets the temperature-induced acceleration. In the near-wall regions of fractures, solute transport exhibits higher sensitivity to temperature variations, with temperature exerting a more pronounced control than the pressure gradient. 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.

       

    /

    返回文章
    返回