本文目录一览：
- 1、实验室测定土的渗透系数方法，从原理来讲可分为 和 两种
- 2、土工试验设备有哪些?
- 3、现场做土的压实度的试验仪器有哪些
- 4、在土石坝里测渗透压的管叫什么
- 5、土石混合体渗透性能的正交试验研究
- 6、土石混合体渗透性能的试验研究
实验室测定土的渗透系数方法，从原理来讲可分为 和 两种
从试验原理上大体可分为”常水头法“和”变水头法“两种！
另请参考：
渗透系数的测定方法主要分“实验室测定”和“野外现场测定“两大类。
1.实验室测定法
目前在实验室中测定渗透系数 k 的仪器种类和试验方法很多，但从试验原理上大体可分为”常水头法“和”变水头法“两种。
常水头试验法 就是在整个试验过程中保持水头为一常数，从而水头差也为常数。 如图：
[常水头法测渗透系数k]
试验时，在透明塑料筒中装填截面为A，长度为L的饱和试样，打开水阀，使水自上而下流经试样，并自出水口处排出。待水头差△h和渗出流量Q稳定后，量测经过一定时间 t 内流经试样的水量V，则
V = Q*t = ν*A*t
根据达西定律，v = k*i，则
V = k*（△h/L）*A*t
从而得出
k = V*L / A*△h*t
常水头试验适用于测定透水性大的沙性土的渗透参数。粘性土由于渗透系数很小，渗透水量很少，用这种试验不易准确测定，须改用变水头试验。
变水头试验法 就是试验过程中水头差一直随时间而变化，其装置如图：
[变水头法测渗透系数]
水从一根直立的带有刻度的玻璃管和U形管自下而上流经土样。试验时，将玻璃管充水至需要高度后，开动秒表，测记起始水头差△h1，经时间 t 后，再测记终了水头差△h2，通过建立瞬时达西定律，即可推出渗透系数 k 的表达式。
设试验过程中任意时刻 t 作用于两段的水头差为△h，经过时间dt后，管中水位下降dh，则dt时间内流入试样的水量为
dVe = －a dh
式中 a 为玻璃管断面积；右端的负号表示水量随△h的减少而增加。
根据达西定律，dt时间内流出试样的渗流量为：
dVo = k*i*A*dt = k*（△h/L）*A*dt
式中，A——试样断面积；L——试样长度。
根据水流连续原理， 应有dVe = dVo，即得到
k = （a*L/A*t）㏑（△h1/△h2）
或用常用对数表示，则上式可写为
k = 2.3*（a*L/A*t）㏒（△h1/△h2）
土工试验设备有哪些?
试验室设备齐全，拥有非饱和土三轴仪，CKC多功能三轴仪、X射线扫描仪、扫描电子显微镜、土水特征仪、土体动强度特性研究设备、共振柱、空心圆柱扭剪仪等先进设备和大量土工常规试验设备
一则消息
现场做土的压实度的试验仪器有哪些
有罐砂筒、环刀、核子密度仪、地基承载力测定仪、水准仪(沉降量观测)、灌水筒。这都是试验仪器
在土石坝里测渗透压的管叫什么
测压管。
1、在土石坝里测渗透压的是测压管。
2、测压管可用来监测坝体浸润线、渗压压力、地下水位及绕坝渗流等。
3、埋设渗压计除可以监测上述项目外，还可用来监测土石坝的孔隙水压力。
土石混合体渗透性能的正交试验研究
周中1 傅鹤林1 刘宝琛1 谭捍华2 龙万学2 罗强2
（1.中南大学土木建筑学院 湖南 长沙 410075
2.贵州省交通规划勘察设计研究院 贵州 贵阳 550001）
摘要 土石混合体作为土和石块的介质耦合体，具有非均质性、非连续性及试样的难以采集性等独特的性质，从而给研究带来极大的困难。土石混合体属于典型的多孔介质，其渗透特性与颗粒的大小、孔隙比及颗粒形状关系密切。本文采用室内正交实验，利用自制的常水头渗透仪，研究了砾石含量、孔隙比和颗粒形状三个因素在不同水平下对土石混合体渗透系数的影响。通过正交试验确定了三种因素对土石混合体渗透系数的影响顺序及各因素的显著性水平。提出了土石混合体渗透系数计算公式，并通过试验结果验证了计算公式的正确性，为土石混合体渗透系数的理论计算提供了一个简明有用的计算工具。
关键词 土石混合体 多孔介质 渗透性能 计算公式 正交试验
土石混合体一般由作为骨料的砾石或块石与作为充填料的粘土或砂组成，是介于土体与岩体之间的一种特殊的地质体，是土和石块的介质耦合体［1］。因为土石混合体具有物质组成的复杂性、结构分布的不规则性以及试样的难以采集性等独特的性质，从而给研究带来极大的困难，目前人们对于它的研究仍处于探索之中［2］。渗透与强度和变形特性，都是土力学中所要研究的主要力学性质，其在土木工程的各个领域中都有重要的作用［3］。土石混合体属于典型的非均质多孔介质［4］，其渗透特性与颗粒的大小、颗粒组成、孔隙比及颗粒形状关系密切。土的渗透系数可以通过室内试验由达西定理计算得出，然而土石混合体的渗透系数却难以确定，主要原因是：取样困难；难以进行常规的渗透试验；大尺度的渗透试验不仅造价高准确性差，而且试验结果离散度大，难以掌握其规律性。迄今为止，国内还没有对土石混合体渗透性能进行研究的资料，现有研究成果局限于利用物理和数值模拟试验对其变形和力学性质进行研究，而对渗透性还未涉及。因此，能够求出土石混合体渗透系数的计算公式具有重要的理论意义和工程应用价值。
本文研究土石混合体中砾石含量、孔隙比（压实度）和颗粒形状三个因素在不同水平下对土石混合体的渗透系数的影响，找出三因素与土石混合体渗透系数之间的关系，并提出土石混合体渗透系数计算公式。
1 土石混合体渗透性能的正交试验
1.1 正交试验方案设计
在室内试验中考虑砾石含量、孔隙比（压实度）和颗粒形状三个因素对土石混合体渗透系数的影响，就每种因素拟考虑3个水平。对于这种3因素3水平的试验，如果考虑每一个因素的不同水平对基材的影响，则根据组合可得有33组试验，这对人力、物力与时间来说都是一种浪费，因此采用正交试验设计来研究这一问题更为合理。本试验所选取的正交表为L9（34），考虑试验误差的影响，但不考虑各因素间的交互作用（即假定他们之间相互没有影响）。共需9组试验，每组作平行试验3次，共27次渗透试验。本试验中采用的因素与对应的水平数如表1所示，其中粗粒形状分为球形体、六面体和三棱锥3个水平，分别由卵石、强风化石块和新打碎的碎石来近似替代。
表1 正交试验的因素水平
1.2 试样的基本物理力学性质
试验所取土样为正在修建的上瑞高速公路贵州段晴隆隧道出口处典型性土石混合体，其天然状态土的物理指标及颗粒级配曲线见表2和图1。由图1可知现场取回土样的不均匀系数Cu为12.31，说明土样中包含的粒径级数较多，粗细粒径之间差别较大，颗粒级配曲线的曲率系数Cc为1.59，级配优良。
表2 天然状态土的基本物理指标
图1 天然状态土的颗粒级配曲线
1.3 大型渗透仪的研制
《土工试验规程》（SL237—1999）规定粗粒土的室内渗透系数需由常水头渗透仪测试，国内常用的常水头渗透仪是70型渗透仪。70型渗透仪的筒身内径为9.44cm，试验材料的最大粒径为2cm，规范［5］要求筒身内径应为最大粒径的8～10倍，因此70型渗透仪的筒身内径过小，有必要研制大尺寸的渗透仪。自制渗透仪的内径和试样高度至少应为最大颗粒粒径的8倍，即至少应为16cm，另外，考虑到边界效应，试样的上下两头分别增加2cm，因此，自制渗透仪的内径和试样高分别取为16cm和20cm。考虑到土石混合体的渗透性较强，选取进排水管的口径为2cm。自制的大型常水头渗透仪的如图2 和图3所示。
图2 常水头渗透仪示意图
数据单位为cm
图3 自制渗透仪
2 试验结果分析
2.1 试验结果
按正交试验表L9（34）的安排，共需作9组试验，每组试验作平行试验3次，取3次测量的平均值，并乘以温度校正系数
，即可求出每组试验20℃时的渗透系数，渗透系数的测量结果见表3。
表3 渗透试验测定结果
续表
2.2 试验分析
运用正交试验的直观分析法和方差分析法，分析各因素对土石混合体渗透系数影响的主次顺序，绘出因素水平影响趋势图，求出各因素的显著性水平。
2.2.1 直观分析
对试验所得的土石混合体的渗透系数进行正交试验的极差分析，并画出各因素的水平影响趋势图。正交试验的极差分析表见表4，3个因素与渗透系数的关系见图4。
表4 极差分析表
图4 各因素与渗透系数的关系
A—砾石含量；B—孔隙比；C—粗粒形状
由正交试验的极差分析表可以看出，对土石混合体渗透系数影响的主次顺序为A→B→C，即砾石含量→孔隙比→颗粒形状。由各因素与渗透系数的关系图可以看出砾石含量越多渗透系数越大，孔隙比越大渗透系数越大，颗粒磨圆度越大渗透系数越小。在路基工程及大坝工程中，可以通过调节粗颗粒的含量、压实度及颗粒形状以获得工程所需的渗透系数。
2.2.2 方差分析
为了确定因素各水平对应的试验结果的差异是由因素水平不同引起的，还是由试验误差引起的，并对影响土石混合体渗透系数的各因素的显著性水平给予精确的数量评估，需采用正交试验的方差分析法对试验数据进行分析，分析结果如表5所示。
表5 方差分析结果
方差分析结果表明：
（1）因素各水平对应的试验结果的差异是由因素水平不同引起的，而不是由试验误差引起的；
（2）砾石含量对土石混合体渗透系数的影响高度显著，孔隙比对土石混合体渗透系数的影响显著，颗粒形状土石混合体渗透系数的影响不显著。
3 土石混合体渗透系数
3.1 渗透系数与砾石含量之间的关系
众所周知，土石混合体的渗透系数与颗粒的大小及级配有关，本文选择等效粒径d20和曲率系数Cc来表示土的颗粒大小和颗粒级配，原因是文献［3］认为等效粒径d20比其他粒径特征系数更能准确地表示颗粒的大小，而与颗粒级配有关的系数是不均匀系数Cu和曲率系数Cc，不均匀系数Cu只反映土粒组成的离散程度，曲率系数Cc能在一定程度上反映颗粒组成曲线的特性，因而曲率系数Cc更适合于评价土的颗粒级配。不同砾石含量的颗粒级配曲线如图5所示。由图5可以求出各曲线的粒径特征系数，见表6。
图5 试样的颗粒级配曲线
表6 不同粗粒含量时的粒径特征
由图6可知，其他条件相同时，土石混合体的渗透系数k与函数f（d20，Cc）呈线性关系，其中
。
图6 k20-f（d20，Cc）关系曲线
3.2 渗透系数与密实度之间的关系
由正交试验的方差分析可知，孔隙率e对渗透系数的影响虽不如粗粒含量大，但也是很显著的。在其他条件相同时，k与
呈线性关系，如图7所示。
土石混合体
3.3 渗透系数与颗粒形状之间的关系
狄凯尔与海阿特（Tikell and Hiatt）于1938年探讨了颗粒的“棱角性”与“圆度”对渗透系数的影响，并指出颗粒的棱角性越大，渗透系数越大［6］。由正交试验分析表可知Cs1∶Cs2∶Cs3＝0.9∶1∶1.2，并且将试验数据进行回归分析，当形状系数Cs1＝0.18，Cs2＝0.2，Cs3＝0.24时与试验结果最为接近，此结论与卡门（Carmen）的研究成果［7］相近。
3.4 土石混合体的渗透系数
由以上分析可知土石混合体的渗透系数与颗粒大小、颗粒级配、颗粒形状及孔隙比有关，同时渗透流体对渗透性也有一定的影响，主要是受液体的动力粘滞度η的影响，大量研究成果表明渗透系数k 与g/η 成正比［3，4，7］。因此，土石混合体的渗透系数计算公式为
土石混合体
式中：k为土石混合体的渗透系数，cm/s；Cs为颗粒的形状系数，m-3；d20为等效粒径，小于该粒径的土重占总土重的20%，m；Cc为颗粒级配曲率系数，
；e为孔隙比；g为重力加速度，9.8 N；η 为液体的动力粘滞度，kPa · s（10-6），η20＝1.01×10-6kPa·s。
由公式（1）计算出20℃时土石混合体的渗透系数k20列于表7。与其他物理力学参数相比，土石混合体的渗透性变化范围要大得多。同时，受宏观构造和微观结构复杂性的影响，其渗透性具有高度的不均匀性［8］。为进一步验证公式（1）的正确性，将实测值与由公式（1）得出的计算值进行对比分析，见图8。由图8可知由公式（1）计算出的渗透系数值与实测值基本吻合，9组试样的平均相对误差为21%，这对于离散性很强的土石混合体的渗透系数来说已经具有足够的精确性。
表7 计算值与实测值对应关系
图8 计算值与实测值关系
4 结论
（1）通过正交试验获取了砾石含量、孔隙比和颗粒形状对土石混合体渗透系数影响的主次顺序，并得出各因素的显著性水平，工程设计中可以通过合理调整土石混合体的砾石含量、孔隙比（压实度）和颗粒形状，以达到控制其渗透能力的目的。
（2）土石混合体的渗透系数与等效粒径d20和曲率系数Cc组成的函数
成正比，并与孔隙比函数
成正比。
（3）提出了土石混合体渗透系数的计算公式，并通过试验结果验证了计算公式的正确性，为土石混合体渗透系数的定量预测提供了一个简明有用的计算工具。
参考文献
［1］油新华.土石混合体随机结构模型及其应用研究.北方交通大学博士论文，2001：1～18
［2］油新华，汤劲松.土石混合体野外水平推剪试验研究.岩石力学与工程学报，2002，21（10）：1537～1540，60～129
［3］刘杰.土的渗透稳定与渗流控制.北京：水利电力出版社，1992：1～20
［4］薛定谔A E.多孔介质中的渗流物理.北京：石油工业出版社，1984：141～173
［5］中华人民共和国水利部.土工试验规程（SL237—1999）.北京：中国水利水电出版社，1999：114～120
［6］ Tickell FG，Hiatt WN.Effect of angularity of grains on porosity and permeability of unconsolidated sands.AAPG Bulletin，1938，22（9）：1272～1274
［7］黄文熙.土的工程性质.北京：水利电力出版社，1984：60～129
［8］邱贤德，阎宗岭，刘立等.堆石体粒径特征对其渗透性的影响.岩土力学，2004，25（6）：950～954
土石混合体渗透性能的试验研究
周中1 傅鹤林1 刘宝琛1 谭捍华2 龙万学2 罗强2
（1.中南大学土木建筑学院 湖南 长沙 410075
2.贵州省交通规划勘察设计研究院 贵州 贵阳 550001）
摘要 土石混合体属于典型的多孔介质，其渗透特性与砾石的百分含量关系密切。通过自制的常水头渗透仪，测定了不同含砾量时土石混合体渗透系数值，研究发现含砾量与土石混合体渗透系数之间存在指数关系；基于幂平均法，提出了土石混合体复合渗透系数的计算公式，并通过试验结果验证了该式的正确性，为土石混合体渗透系数的理论计算提供了一个简明有用的计算工具。
关键词 土石混合体 多孔介质 渗透性能 复合渗透系数 经验公式
土石混合体一般是由作为骨料的砾石或块石与作为充填料的粘土或砂组成，它是介于土体与岩体之间的一种特殊的地质体，是土和石块的介质耦合体［1］。因为土石混合体具有物质组成的复杂性、结构分布的不规则性以及试样的难以采集性等特殊的性质，从而给研究带来极大的困难，目前人们对于它的研究仍处于探索之中［2］。渗透与强度和变形特性都是土力学中所要研究的主要力学性质，其在土木工程的各个领域都有重要的作用［3］。土石混合体属于典型的非均质多孔介质［4］，其渗透系数是由高渗透性的砾石和低渗透性的土体复合而成的。土的渗透系数可以通过室内试验由达西定理计算得出，然而土石混合体的渗透系数却难以确定，主要原因是：取样困难；难以进行常规的渗透试验；大尺度的渗透试验不仅造价高、准确性差，而且试验结果离散度大，难以掌握其规律性。因此能够求出土石混合体复合渗透系数的计算公式具有重要的理论意义和工程应用价值。
土石混合体中土与砾石粒径的界限值为5mm，即将粒径小于5mm的颗粒称为土、大于5mm的颗粒称为石，砾石含量用P5表示［1］。利用自制的常水头渗透仪，研究砾石体积百分含量P5从0%逐步过渡到100%（间隔10%）时土石混合体的渗透系数，每种配比作平行试验3次，共33次渗透试验。
1 土石混合体渗透性能试验
1.1 试样的基本物理力学性质
试验所取土样为正在修建的上瑞高速公路贵州段晴隆隧道出口处典型性土石混合体，其天然状态土的物理指标及颗粒级配曲线见表1和图1。由图1可知现场取回土样的不均匀系数Cu为12.31，说明土样中包含的粒径级数较多，粗细粒径之间差别较大，颗粒级配曲线的曲率系数Cc为1.59，级配优良。
表1 天然状态土的基本物理指标
图1 天然状态土的颗粒级配曲线
1.2 大型渗透仪的研制
《土工试验规程》（SL237—1999）规定粗粒土的室内渗透系数需由常水头渗透仪测试，国内常用的常水头渗透仪是70型渗透仪。70型渗透仪的筒身内径为9.44cm，试验材料的最大粒径为2cm，规范［5］要求筒身内径应为最大粒径的8～10倍，因此70型渗透仪的筒身内径过小，有必要研制大尺寸的渗透仪。自制渗透仪的内径和试样高度至少应为最大颗粒粒径的8倍，即至少应为16cm，另外，考虑到边界效应，试样的上下两头分别增加2cm，因此，自制渗透仪的内径和试样高分别取为16cm和20cm。考虑到土石混合体的渗透性较强，选取进排水管的口径为2cm。自制的大型常水头渗透仪如图2和图3所示。
图2 自行研制的渗透仪
图3 常水头渗透仪示意图
数据单位为cm
1.3 试验步骤
首先，将由现场取回的土样烘干、过筛，并根据粒径的大小分为0～5 mm的土和5～20mm的砾石两部分。然后，按照试验要求的砾石体积百分含量P5，以10%的初始含水量配制试样，静置24 h。试验时，将配制好的试样分层装入圆桶中，每层装料厚度30mm左右，分层压实，记录每层的击实数。按上述步骤逐层装样，至试样顶部高出测压孔约3cm为止。测出装样高度，准确至0.1cm。在试样顶部铺一层2cm厚的细砾石作缓冲层。之后，由进水管注入蒸馏水，直至出水孔有水流出，静置24 h使试样充分饱和。用量筒从渗透水出口测定渗透量，同时用温度计测量水温，用秒表测记经一定时间的渗水量，共测读6次，取其平均值，6次结果相差不得超过7%，否则需重新测定。
1.4 试验数据
按照试验设计的各种砾石体积百分含量P5共需作11组试验，每组试验作平行试验三次，取3次测量的平均值，并乘以温度校正系数
，即可求出每组试验20℃时的渗透系数，渗透系数的测量结果见表2。
表2 渗透系数测定结果
2 试验结果分析
2.1 渗透系数与砾石含量的关系
不同含砾量的颗粒级配曲线如图4所示，由图4可以求出各曲线的粒径特征系数及不均匀系数Cu和曲率系数Cc。
图4 试样的颗粒级配曲线
图5为土石混合体砾石含量P5与20℃时渗透系数的关系曲线。从图5可以看出，随着含砾量的增加，渗透系数急剧增加，可见，在设计中可以通过调节砾石的含量来控制土石混合体的宏观渗透性能。
图5 粗粒含量与渗透系数的关系
从图5还可以发现，土石混合体中砾石的含量P5与渗透系数k之间存在指数关系，与文献［6］的研究成果相似，即
土石混合体
式中：k0为P5＝0时土的初始渗透系数；n为与土石混合体本身性质相关的常数。对于文中试验值，k0与n分别为0.0006cm/s和8.82。在工程中可以通过少量试验来确定k0，n值，以此来预测不同级配土石混合体的渗透性。
2.2 土石混合体的复合渗透系数
近几十年来，许多学者在揭示影响和决定土的渗透系数内在因素及其相互关系方面进行了大量工作，并取得了有益的成果［7～12］，被认为依然有效且目前常用的确定渗透系数的半经验、半理论公式有：
（1）水利水电科学研究院公式［7］：
土石混合体
式中：k10，k20分别为温度为10℃和20℃时的渗透系数（cm/s）；η10/η20为温度为10℃和20℃的粘滞系数比；n为孔隙率；d20为等效粒径（mm）。
（2）泰勒（Taylor）［9］用毛管流的哈根-伯努力（Hange-Poiseuille）方程导出渗透系数的表达式：
土石混合体
式中：ds为当量圆球直径，可以用等效粒径d20代替；γw为液体容重；μ为液体粘滞度；e为孔隙比；C为形状系数，通常取C＝0.2。
式（2）和式（3）均是针对土体的渗透特性提出的半经验、半理论公式，然而对于非均质性更强、粒径差别更大的土石混合体来说，其适用性不是很强。土石混合体中砾石形成骨架，细颗粒充填孔隙，其渗透系数是由低渗透介质土体的渗透系数kS和高渗透性介质砾石的渗透系数kG复合而成。土石混合体复合渗透系数不是按体积百分含量的简单复合，而是高低渗透性介质的耦合。在参考相关文献［10～12］的基础上，基于幂平均法，本文提出的土石混合体复合渗透系数k复合的表达式为
土石混合体
式中：P5为砾石的体积百分含量，%；kG为砾石的渗透系数，cm/s；kS为土的渗透系数，cm/s；f为系数。
砾石的体积百分含量P5可以由筛分法求出；土的渗透系数kS和砾石的渗透系数kG可以由室内试验直接求出或参考相关资料确定；系数f可以通过少量试验回归分析确定，因此可以说（4）式是一个简明实用的土石混合体复合渗透系数计算公式。
图6 不同计算方法结果比较
为进一步验证（4）式，我们将试验测得的k值与用（2），（3），（4）式计算得到的k值进行对比分析。结果见图6，具体数值见表3。由图6和表3可知据水利水电科学研究院公式和泰勒公式计算结果均高于实测值，尤其是当P5≤30%时，（2）式计算结果和（3）式计算结果比实测值大2～3个数量级，与实测值相差较大。而用本文方法得到的土石混合体的渗透系数最接近实测值，平均相对误差仅为0.6%，能够作为土石混合体渗透系数定量预测的有效工具。在工程设计中，可以根据工程对土石混合体渗透性的要求，依据本文提供的经验公式，调整土石混合体中砾石的含量，达到控制土石混合体渗透能力的目的。
表3 土石混合体渗透系数及相关参数
3 结论
（1）利用自制的常水头渗透仪，测定了不同含砾量时土石混合体的渗透系数值，并指出含砾量与土石混合体渗透系数之间存在指数关系。在工程设计中可以通过合理调整土石混合体中砾石的含量，达到控制其渗透性能的目的。
（2）指出土石混合体的渗透系数是一种由高渗透性的砾石和低渗透性的土体复合而成的，给出了土石混合体复合渗透系数的计算公式，并通过试验结果验证了计算公式的正确性，为土石混合体渗透系数的定量预测提供了一个简明有用的计算工具。
参考文献
［1］油新华.土石混合体随机结构模型及其应用研究.北方交通大学博士论文.2001：1～18
［2］油新华，汤劲松.土石混合体野外水平推剪试验研究.岩石力学与工程学报.2002，21（10）：1537～1540，60～129
［3］黄文熙.土的工程性质.北京：水利电力出版社.1984：60～129
［4］薛定谔 A E.多孔介质中的渗流物理.北京：石油工业出版社.1984：141～173
［5］中华人民共和国水利部.土工试验规程（SL237—1999）.北京：中国水利水电出版社，1999：114～120
［6］邱贤德，阎宗岭，刘立等.堆石体粒径特征对其渗透性的影响.岩土力学，2004，25（6）：950～954
［7］刘杰.土的渗透稳定与渗流控制.北京：水利电力出版社，1992：1～20
［8］ Wen X H，Gomez-Hernandez J J.Upscaling hydraulic conductivities in heterogeneous media：An overview.Journal of Hydrology，1996，183：ix～xxxii
［9］ Taylor D W.Fundamentals of soil mechanics.John Wiley ＆ SONS，Inc.，1948
［10］ Brown W F.Solid mixture permitivities.Journal of Chemical Physic，1955，23（8）：1514～1517
［11］ Dagan G.Analysis of flow through heterogeneous random aquifers by the method of embedding matrix—1：Steady flow.Water Resources Research，1981，17（1）：107～122
［12］ Noetinger B.The effective permeability of a heterogeneous porous medium.Transport in Porous Media，1994，15：99～127