[K=φ3CS2] （1）
式中：[K]为多孔介质的渗透率；[φ]为多孔介质的孔隙率；[C]为KC常数；[S]是比表面积［8］。KC常数是Kozeny在1927年首次提出的，推导出KC常数为5。KC常数是一个经验参数，它被证明不是一个常数，而且对于不同介质不同环境中KC常数都不一样。这引起了研究人员的广泛关注，使得研究者们不断修正KC方程来提高计算精度。Xu等［9］回顾并总结了KC方程和KC常数相关的各种模型，推导出了不含经验参数的KC常数解析表达式。但是研究者只关注于充分发展的牛顿流体层流流动和各向同性松散多孔介质的影响，研究结果表明KC常数与孔隙率有关。Xiao等［10］ 根据分形理论，推导出了多孔纳米纤维材料的KC常数解析表达式，发现其与孔隙率、面积分形维数有关。在低孔隙率时，KC常数随着孔隙率的增加而增大，但与实验结果比较误差较大。Wei等［11］基于水力孔径的分形分布特性，分析了毛细管吸胀过程中的KC常数，但是受吸胀过程中各种因素的影响，无法准确的预测任意条件下多孔介质的KC常数。Xiao等［12］导出了纤维多孔介质中KC常数的分形模型，结果表明其与纤维多孔介质本身的微观结构参数密切相关，但是该研究没有考虑多孔介质表面粗糙度对渗流过程的影响。Xiao等［13］推导出了粗糙表面纤维多孔介质的KC常数，但是忽略了其他非线性因素的影响。结果发现KC常数随着粗糙度、孔隙率、孔隙面积分形维数和弯曲度分形维数的增加而增大。前人的大量研究表明了KC常数与多孔介质孔隙率密切相关，但对KC常数与多孔介质结构参数间的关系研究还不够深入，并且树状分叉输运网络中的KC常数研究目前还未涉及。为了分析树状分叉网络KC常数，本文采用了一个理想化的模型（嵌入在多孔介质中的分形树状网络）来模拟真实的多孔介质［14］。最后推导出KC常数的解析表达式并讨论了孔隙率和网络的微观结构参数对KC常数的影响。
1 嵌入在多孔介质中的分形树状网络复合材料的结构参数
复合材料由分形树状网络和基质多孔介质两部分组成［15］。
1.1 分形树状网络结构参数
树状分叉网络如图1所示。为了建立一个复合材料中流体渗流模型，忽略树状网络每层的壁的厚度，调整每层管子的分支角度和体积，使其不相互作用，每个通道在下一层被划分为[n]个分支（这里[n]=2），分支的最大分叉级数为[m]。为了描述分形树状网络的结构，对于其中任意一级分叉而言，第[k]（[k]=0，1，2，…，[m]）级分叉管道的长度和直径分别为[lk]和[dk]，分叉角为[θ]。用来描述分叉网络几何结构的长度比（[α]）和直径比（[β]）可以定义为［16］：
[α=lk+1lk] （2）
[β=dk+1dk] （3）
推出：
[lk=l0αk] （4）
[dk=d0βk] （5）
其中[l0]和[d0]分别表示第0分支管道的长度和直径。
<G:\武汉工程大学\2022\2022-06工程\Image\李子豪-1.tif>[dk][lk][lk+1][dk+1][2θ]
图1 树状分叉示意图［长度单位：μm，角度单位：（°）］
Fig. 1 Tree bifurcation diagram（Length unit：μm，
Angle unit：degree）
1.2 复合材料结构参数
嵌入在多孔介质中分形树状网络复合材料如图2所示，复合材料的宽度[b]、长度[L0]和分形树状网络的实际长度[Lt]的表达式可写出：
[b=k=1m2lksinθ=2l0sinθα-αm+11-α] （6）
[L0=l0+k=1mlkcosθ=l0(1+α-αm+11-αcosθ)] （7）
[Lt=k=0mlk=l01-αm+11-α] （8）
根据体积孔隙度的定义，可以得到复合网络的体积孔隙度［17］：
[φ=VfV=k=0mnkπd2klk4bd0L0=]
[πd08l0(1-α)[1-(nαβ2)m+1]sinθ[1+α-αm+11-αcosθ](α-αm+1)(1-nαβ2)] （9）
其中[Vf]表示树状分叉网络的体积，[V]是复合网络总体积。
2 树状分叉网络的渗透率
根据 Hagen-Poiseuille方程，第[k]级单根通道中流量可写为：
[Qk=πd4k128μΔpklk] （10）
其中，[μ]是流体的黏度，[Δpk]是第[k]级通道的压降［18］。根据 Darcy定律：
[Qk=Kkμ(dk2)2πΔpklk=πd2kKk4μΔpklk] （11）
其中[Kk]是第[k]级单根通道的渗透率。联立方程（10）、（11）可得：
[Kk=d2k32] （12）
根据渗流并联模型，第[k]级[nk]根管道并联起来的有效渗透率为：
[Kek=i=0nkAkiKkiAek] （13）
其中，[Aki=πd2k/4]是第[k]级单根分叉管道的横截面积，[Kki=Kk=d2k/32]是第[k]级单根分叉管道的渗透率，[Aek=i=0nkAki=nkπd2k/4]是第[k]级[nk]根并联管道的等效横截面积。由此可得到第[k]级[nk]根管道并联后的有效渗透率：
[Kek=d2k32] （14）
可以发现第[k]级[nk]根管道并联起来的有效渗透率与第[k]级单根管道的渗透率相等。可以将该网络视为由[m]个长度为[lek]，横截面积为[Aek]，渗透率为[Kek]的单通道串联而成，整个网络可以等效为一个渗透率为[Ke]的单通道。[Ke]也就是整个网络的有效渗透率。对于整个网络的等效单管，由广义Darcy定律：
[Q=KeμAeΔple] （15）
其中[Ae]和[le]分别是整个网络等效单管的横截面积和长度，[Δp=k=0mΔpk]是网络的总压降。由于质量守恒，第[k]级管道的流量与网络总流量相等
[Q=KekμAekΔpklek] （16）
[Δp=k=0mQμlekKekAek] （17）
[le=k=0mlk=l01-αm+11-α] （18）
根据Hagen-Poiseuille方程：
[Q=πd4e128μΔple] （19）
其中[de]为等效直径。
[Ae=πd4e/4] （20）
联立方程（14~20）可以得到整个分叉网络的有效渗透率：
[Ke=d20321-αm+11-α1-α/nβ41-(α/nβ4)m+112] （21）
在渗流过程中，流体的流动不可能是沿直线流动，而是弯曲曲折地流动。因此还需要考虑到迂曲度[T]对渗流过程的影响，迂曲度定义为渗流通道的实际长度[Lt]与直线长度[L0]的比值［19-20］。
[T=LtL0=1-αm+1(1-α)1+α-αm+1(1-α)cosθ] （22）
考虑迂曲度的影响后，网络的有效渗透率为：
[K=Ke1T=]
[d20321-α1-αm+11-α/nβ41-(α/nβ4)m+1121+α(1-αm)1-αcosθ] （23）
3 嵌入在多孔介质中树状分叉网络的KC常数
树状分叉网络中各级管道的总内表面积[At]可写为：
[At=k=0mnkπdklk=πd0l01-(nαβ)m+11-nαβ] （24）
于是可以得到比面[S]的表达式：
[S=AtV=πd0l01-(nαβ)m+11-nαβbd0L0=]
[π1-(nαβ)m+11-nαβ2l0sinθα-αm+11-α(1+α-αm+11-αcosθ)] （25）
联立方程（1）、（9）、（23）、（25）可以得到嵌入在多孔介质中树状分叉网络的KC常数表达式：
[C=πd04l0×]
[ (1-α)[1-(nαβ2)m+1]3sinθ1+α-αm+11-αcosθ2(α-αm+1)(1-nαβ2)3×][ (1-nαβ)2[1-(nαβ)m+1]21-α1-αm+11-α/nβ41-(α/nβ4)m+1-12]
（26）
4 分析与讨论
方程（26）给出了树状分叉网络复合材料的KC常数解析式，很明显它与分叉网络的微观结构有着直接的关系。从图3中可以看到KC常数随着孔隙率[φ]的增加而增大，这与前人的研究结果相吻合。图3显示KC常数随着分叉网络的长度比[α]的增加而增大，这是因为随着长度比的增加，下一级分叉管道就越长，流体流动阻力就越大，渗透率降低，根据方程（1），KC常数增大。图3显示KC常数随着直径比[β]的增加而减小，这是由于直径比的增加，下一级分叉管道就越粗，流体流量增加，渗透率升高，根据方程（1），KC常数减小。图4列举了一些前人研究成果中KC常数的预测值，并与本文推算出的预测值作比较，发现在孔隙率较小（[0<φ<0.2]）时，本文预测结果和Sullivan［21］及徐鹏等［22］的结果较吻合；在[0.2<φ<0.3]时，本文预测结果和Devies［23］等人及Sparrow［24］等的结果十分接近；当孔隙率在0.3附近时，本文预测结果符合Happel［25］等的结果，此时分形模型的预测值接近经典的KC常数值；在[0.4<φ<0.5]时，本文预测结果和Kyan等［26］及Sahraoui等［27］的预测值较接近。综上发现模型预测结果与实验结果吻合较好，这很好地验证了本模型的正确性。KC常数和孔隙率间的定量关系因为不同学者的研究方法和手段的差异，在学界还未有定论，故不断有学者提出讨论。本文研究的是树状分叉网络的KC常数，而参考文献中模型多是单根管道，所以数据对比会有差异性，但是数值量纲与前人研究成果相符合。
5 结 论
前人的大量研究表明了KC常数与多孔介质孔隙率密切相关，但是对于树状分叉输运网络中的KC常数研究目前还未涉及，而本文创新的推导出了嵌入在多孔介质中的树状分叉网络复合材料的KC常数的解析表达式，发现KC常数与孔隙率和网络结构参数有关系。结果表明树状分叉网络KC常数随着复合网络孔隙率的增加而增大，随着分叉网络长度比的增加而增大，随着分叉网络直径比的增加而减小。当孔隙率在0.3附近时，本文分形模型的预测值最接近经典的KC常数值（5）。尽管本文推导出了嵌入在多孔介质中树状分叉网络复合材料KC常数的解析表达式，但是该结果忽略了树状分叉管道壁表面的粗糙元素，为了进一步深入的研究树状分叉网络中的KC常数，需要将粗糙度对渗流过程的影响考虑进去。