多孔介质非均质性如何控制搅拌与弥散：二维稳态达西流中流体元拉伸与弥散的升尺度随机模型

《Journal of Fluid Mechanics》：How porous medium heterogeneity controls stirring and dispersion

　　

研究人员研究二维非均质多孔介质中稳态流动下的搅拌（stirring）与弥散（dispersion）。多孔介质结构通过引起空间流速变异和流体元变形来搅拌流体。弥散度量溶质分布的扩展，而拉伸（stretching）支撑溶质混合、化学反应和生物活动。尽管如此，关于多孔介质中流体拉伸的认识仍少得多。其一方面与弥散过程的关系，另一方面与流动和介质属性的关系，仍是开放问题。为弥补这一空白，研究人员探究导致流体元拉伸和弥散的介质非均质性机制，并量化支配其动力学的升尺度规律（upscaled laws）。因此，研究人员对以空间变化水力传导率为特征的二维多孔介质中的达西流（Darcy flow）进行详细数值模拟。弥散通过流体元的首达时间分布（first-passage time distributions）及其位移均值和方差的时间演化分析，变形则通过流体元的概率密度函数（probability density function, PDF）和平均伸长率（average elongation）研究。研究人员发现弥散和拉伸行为密切相关。支配首达时间分布衰减和位移方差演化的指数，也决定平均伸长率的演化。该动力学被量化于一个单一升尺度理论框架中，该框架连接欧拉流动结构（Eulerian flow structure）、大尺度弥散和流体元拉伸。

论文解读

研究背景与问题。非均质多孔与裂隙介质中的流动和溶质输运受水力传导率空间变异控制，该变异可跨越多个数量级。稳态流沿介质特征尺度组织，并产生宽流速分布及宽停留时间分布，导致非Fickian大尺度输运，表现为宽首达时间分布（first-passage time distribution，粒子到达控制面时间的分布）和非线性弥散（dispersion，溶质分布扩展）演化。弥散、搅拌（stirring，流体元变形与拉伸）和混合是不同过程：孔隙尺度上介质通过复杂孔隙空间流速变异搅拌流体，三维孔隙空间对流通常混沌、拉伸指数增长；达西尺度上，各向同性多孔介质稳态流受拓扑约束不能混沌，拉伸必为亚指数。纵向弥散仅由沿流线速度波动控制，横向弥散由流线几何控制，而搅拌由垂直于流线的速度梯度（剪切和涡量）控制。二者由同一水力传导率场控制，预期有强联系。二维稳态流中Lyapunov指数为零，纯对流横向弥散为零，但其与各向同性非均质多孔介质中稳态三维流类似，因为该类流螺旋度密度（helicity density）为零、流线限于二维流面。已有连续时间随机游走（continuous time random walk, CTRW，升尺度随机输运模型）和耦合模型描述稳态随机流变形，但非均质达西流中拉伸率与介质属性关系未建立；Le Borgne et al. 的乘法随机过程需由拉伸观测校准，未揭示与流场/介质结构及弥散动力学关系。因此，研究人员同时分析二维非均质多孔介质中流体元的弥散和变形，识别非均质性机制，并将大尺度拉伸与弥散联系到欧拉流组织和介质结构；其意义涉及远大于深度的含水层、粗糙裂隙中流与对流主导输运及三维各向同性多孔介质。

技术方法（≤250字）。研究人员以多高斯随机场生成水力传导率场，包括对数正态分布（σ_f²=0.25,1,3,6）和伽马分布（α=0.25,0.5,0.75）；在100?_Y×200?_Y、1000×2000节点域上生成200个实现。采用有限体积两点通量近似求解达西流，三阶样条插值获取速度梯度；用Pollock方法追踪粒子，并以Protean坐标（随流线方向旋转的局部坐标）数值积分变形方程及双粒子法计算线元演化。观测量包括首达时间分布、位移均值与方差、平均伸长率和伸长率PDF；并建立速度与剪切率的CTRW/Ornstein–Uhlenbeck过程（OU过程，平稳高斯马尔可夫过程）。

研究结果。3.1 Continuous time random walk for particle dispersion：通过将沿流线等距采样的速度序列建模为平稳马尔可夫过程，以正态分数变换和OU过程表述CTRW，给出位移均值、方差和首达时间分布表达式；结论是该CTRW可作为稳态随机流中粒子弥散的升尺度模型。3.2 Stochastic dynamics of deformation：通过分析剪切率σ(s)与速度v(s)统计，发现条件平均绝对剪切率与速度线性相关?|σ||v?=k_cv，ξ(s)近似高斯、相关函数近似指数并在与速度相同?_c尺度去相关，从而将剪切率建模为OU过程，得到变形闭合随机模型。3.3 Summary：将粒子位置、速度、时间、线元分量、ζ和ξ的演化写成差分CTRW方程组，形成同时描述弥散和拉伸的单一框架；结论是弥散与变形可在同一CTRW中统一。3.4 Discussion of model parameters：通过速度PDF和参数分析，CTRW输入为p_e(v)、χ、?_c、k_c；伽马场小速度p_e(v)∝v^β?1（β=0.3,0.55,0.85对应α=0.25,0.5,0.75）并在v_m截止；弱非均质有?_c=8/3?_Y，强非均质用经验式；结论是模型参数可由欧拉速度统计和介质相关尺度确定。4.1 Dispersion：通过比较数值模拟与CTRW的首达时间分布和位移方差，CTRW再现非Fickian特征；伽马场f(t,x₁)∝t^?1?β，位移方差κ(t)∝t^2?β，在τ_c?t?τ_m；对数正态场小σ_f²≤1无重尾，大σ_f²有显著拖尾和幂律交叉，长期方差渐近Fickian但首达时间仍可重尾，逆高斯拟合不能捕捉尾部；结论是弥散非Fickian可由速度分布经CTRW解释。4.2 Stretching：通过平均伸长率λ(t)及p_?(l,t)比较数值与CTRW，早期线性，随后幂律λ(t)∝t^γ；伽马场0<β<1时γ=2?β，对数正态场落入β>3类γ=1/2；伸长率PDF短时尖峰、随时间向大值偏斜，CTRW与数值一致；结论是介质非均质性经速度PDF控制拉伸，平均伸长率和分布可由CTRW描述。4.3 Discussion：通过Protean坐标和CTRW分析，指出多孔介质搅拌不同于混沌/湍流：稳态达西流拓扑约束使拉伸亚指数、由变持续时间剪切事件主导，持续时间随局部流速变化；该框架从第一原理推导随机拉伸模型，避免Le Borgne et al.需拟合几何布朗运动的局限，并将拉伸与弥散联系到欧拉速度分布。

讨论与结论翻译。讨论表明，多孔介质非均质性通过欧拉速度分布控制搅拌；拉伸导致垂直于材料条带的浓度梯度，加速混合和反应，强拉伸区与反应热点有关，可解释溶解/沉淀、岩溶和矿床分布。研究结论：研究人员研究非均质二维多孔介质中流体元拉伸与弥散的机制和随机动力学，定量联系拉伸与弥散以及拉伸与介质/流动属性，建立升尺度建模方法并识别预测大尺度流体变形和弥散的关键参数，这些是混合模型的关键输入。研究人员在Protean坐标中表述变形问题，揭示粒子速度和剪切率拉格朗日序列对材料条带伸长演化的关键作用；沿流线等距采样速度与剪切率，以平稳马尔可夫过程量化，均由单一相关长度表征；速度用正态分数变换的OU过程表示，剪切率满足与速度线性相关、条件分布高斯、与速度同相关尺度演化，故也建模为OU过程，从而形成统一CTRW框架。该框架用于伽马分布和对数正态水力传导率场直接数值模拟的流体粒子变形分析，确认首达时间分布和位移矩的弥散结果，观察到宽首达时间分布和非线性位移方差等非Fickian特征；这些特征也反映在平均伸长率中，初始线性演化到由欧拉速度PDF决定的渐近幂律标度，并由CTRW拉伸模型完全预测；粗粒化变形过程为耦合CTRW，在幂律转移时间分布下为Lévy walk（Lévy行走，耦合CTRW的幂律转移时间极限），从而将标度与速度PDF标度解析联系。总之，所导出的CTRW模型以欧拉速度PDF预测拉伸和弥散行为，阐明介质非均质性如何控制搅拌。理解非均质性对拉伸的控制对于量化和预测多孔介质中混合与化学反应至关重要。