本发明属于离子型稀土溶浸采矿技术领域，尤其涉及一种定水头条件下离子型稀土单井注液影响范围的测定方法。
背景技术：
目前，业内常用的现有技术是这样的：离子吸附型稀土(简称离子型稀土)是我国特有、世界关注的重要矿产资源，其中含有丰富的中重稀土，主要分布在江西、福建、广东、湖南、广西、云南、浙江等七省区。在该类矿床中，稀土以水合阳离子或羟基水合阳离子形式吸附在黏土矿物上，具有配分齐全、高附加值元素含量高、放射性比度低、综合利用价值大等特点。离子型稀土开采先后经历了池浸、堆浸和原地浸矿，目前主要采用原地浸矿的采矿方法，即直接在山体表面布置注液井网，注入盐类溶浸液渗透到黏土矿物中，通过离子交换过程将稀土离子交换于溶液中，浸出液通过山脚处的积液沟或收液巷道汇集到母液池，再通过沉淀母液中的稀土，实现资源回收。原地浸矿开采设计的关键因素是计算井网参数，如果注液井网布置太密，不仅增加注液井的钻孔工程量，且容易造成山体过度饱和，引起山体滑坡；如果注液井网布置太疏，则会形成许多浸矿死角与盲区，从而影响稀土资源充分回收。设置井网参数的主要依据是单井注液影响范围，即入渗过程的湿润体特征及运移规律，因此，确定单井注液影响范围对于降低浸矿剂消耗、提高稀土资源回收率，以及控制土体变形破坏具有重要意义。国内外学者研究了关于土壤入渗的湿润体特征及渗流规律，但是由于离子型稀土开采长期处于技术保密阶段，对于离子型稀土溶浸采矿过程入渗的湿润体特征及渗流规律很少，导致在开采过程中，稀土资源浸取率不确定性因素大，诸如“稀土去哪儿了”、“浸矿液流去哪里了”等问题悬而未决，如果能建立注液时离子型稀土单井注液影响范围的计算模型，可以有效解决上述问题，提高生产效率。基于上述研究，取江西龙南足洞稀土矿区稀土进行三维入渗试验，研究了不同入渗水头下湿润体形状及湿润锋运移规律，推求了横向湿润锋距离、垂向距离与湿润体体积计算模型，并与试验值进行比较，模型精度较高，此模型可以进一步提高稀土开采效率。
综上所述，现有技术存在的问题是：
(1)原地浸矿过程中溶浸液通过单个注液井向稀土矿体扩散的过程是一个动态过程，诸多因素影响溶浸液入渗过程，如何通过注液水头高度等因素确定单井注液影响范围是预测稀土开采全过程的关键问题之一。
(2)原地浸的注液井有一定的注水高度，即存在一定的水头压力，有别于其他没有水头压力的入渗过程，把水头高度作为一个重要参数建立离子型稀土浸矿的三维入渗模型，可以大幅度提高入渗模型的精准性和适用性。
(3)现有一维入渗规律及入渗模型研究较多，三维入渗规律及入渗模型较少，一维入渗和三维入渗的研究方式和研究方法具有差别，一维入渗规律及入渗模型是简化后的，具有一定的局限性。
解决上述技术问题的难度和意义：原地浸矿开采设计的关键因素是计算井网参数，如果注液井网布置太密，不仅增加注液井的钻孔工程量，且容易造成山体过度饱和，引起山体滑坡；如果注液井网布置太疏，则会形成许多浸矿死角与盲区，从而影响稀土资源充分回收。设置井网参数的主要依据是单井注液影响范围，即入渗过程的湿润体特征及运移规律，因此，通过确定单井注液影响范围计算模型可以有效预测和适时调控离子型稀土开采中溶浸液流动的全过程，对于降低浸矿剂消耗、提高稀土资源回收率，以及防治山体滑坡和控制环境污染具有重要意义。
技术实现要素：
针对现有技术存在的问题，本发明提供了一种定水头条件下离子型稀土单井注液影响范围的测定方法。
本发明是这样实现的，一种湿润体特征值计算模型，所述湿润体特征值计算模型包括：湿润体半径模型和湿润体体积模型；通过此模型可以有效预测离子型稀土开采中浸矿液的影响半径和浸矿体积，对浸矿全过程实施模拟计算。
所述湿润体半径模型为：
当sinθ＝1时，r为湿润体垂向最大湿润深度，用rb表示，得：
得到时间t与横向最大湿润距离、垂向最大湿润深度之间的关系。
式中，θ为水平面与微元体的夹角，无量纲；k为饱和导水率，量纲l/t；r0为注液井半径，量纲l；h为水头高度，量纲l。θi为土体初始含水率，量纲l3/l3；θs为土体饱和含水率，量纲1；当sinθ＝0时，r为湿润体横向最大湿润距离，用ra表示；
所述湿润体体积模型为：
上式中，gv是一个关于v的隐函数，式中的三维变量被用来计算v和t之间的几何关系。
本发明的另一目的在于提供一种基于所述湿润体特征值计算模型的定水头条件下离子型稀土单井注液影响范围的测定方法，所述定水头条件下离子型稀土单井注液影响范围的测定方法包括：
步骤一，将土样自然风干，捣碎，过筛，测得土样的基本物理参数；
步骤二，按预定容重1.56g/cm3将土样分层装入土箱内，每装入一层土样，用木锤轻轻将其均匀捣实，将捣实后的层面抓毛再装入下一层；
步骤三，入渗水头分别设置为5cm、10cm、15cm、20cm、25cm，入渗时间为8～10h。
进一步，所述定水头条件下离子型稀土单井注液影响范围的测定方法入渗的完整湿润体近似为半个椭球体；
湿润体体积v和累计入渗量w之间存在显著的线性关系，拟合方程为：
w＝0.326v，r＝0.9774；
表明土体湿润体内平均含水率增量基本上保持为定值，δθ＝0.326；干土层的含水率为初始含水率θi，湿润体内的稳定含水率为饱和含水率θs，则δθ＝θs－θi＝0.326；含水率是体积含水率；
横向入渗距离a(t)和垂向入渗距离b(t)与时间t的关系，用幂函数、指数函数和对数函数进行回归分析；幂函数形式为：a(t)＝atb；b(t)＝ctd，参数a、b、c、d为湿润锋距离和入渗时间关系的拟合参数。
本发明的另一目的在于提供一种所述湿润体特征值计算模型在土层或岩石的钻进中的应用。
综上所述，本发明的优点及积极效果为：本发明通过不同的恒定水头条件下稀土入渗试验可知，其湿润体形状近似为半椭球体，入渗初期横向湿润距离较垂向湿润距离大，随着入渗的进行，垂向湿润锋较横向湿润锋的发展速度增大，湿润体形状由半椭球体向半球体发展；湿润体体积和累计入渗量之间存在显著的线性关系，说明土体湿润体内平均含水率增量基本上保持为一定值；横向和垂向最大湿润距离与入渗时间存在极显著的幂函数关系；基于green-ampt模型和水量平衡原理，建立了湿润体特征值计算模型，模型结果与试验结果吻合度较高，具有应用价值。根据三维green-ampt模型入渗能力曲线，说明三维入渗率随着时间的增加而增加，是三维入渗模型与一维模型最显著的差异。
附图说明
图1是本发明实施例提供的定水头条件下离子型稀土单井注液影响范围的测定方法流程图。
图2是本发明实施例提供的湿润体形状观测结果示意图；
图中：(a)h＝5cm；(b)h＝10cm；(c)h＝15cm；(d)h＝20cm；(e)h＝25cm。
图3是本发明实施例提供的湿润体扁率示意图。
图4是本发明实施例提供的湿润体体积与累计入渗量的关系示意图。
图5是本发明实施例提供的入渗水头对水平湿润锋的影响示意图。
图6是本发明实施例提供的入渗水头对垂向湿润锋的影响示意图。
图7是本发明实施例提供的湿润体特征参数分析图。
图8是本发明实施例提供的湿润锋距离与时间的关系曲线示意图。
图9是本发明实施例提供的横向湿润锋计算值与实测值示意图。
图10是本发明实施例提供的三维green-ampt模型入渗能力曲线示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明在离子型稀土原地浸矿中，合理设置注液井网有利于提高资源回收率，计算井网参数的关键环节是确定单井注液影响范围，即湿润体特征及其运移规律。
如图1所示，本发明实施例提供的定水头条件下离子型稀土单井注液影响范围的测定方法包括以下步骤：
s101：将土样自然风干，捣碎，过筛，测得土样的基本物理参数；
s102：按预定容重(1.56g/cm3)将土样分层(5cm)装入土箱内，每装入一层土样，用木锤轻轻将其均匀捣实，将捣实后的层面抓毛再装入下一层；
s103：入渗水头分别设置为5cm、10cm、15cm、20cm、25cm，入渗时间要尽可能的长，一般为8～10h。
下面结合试验对本发明的应用原理作进一步的描述。
1试验材料与方法
1.1试验装置
整个试验系统由供水瓶、试验土箱、注液井、溢流管、电子秤等组成，根据单井注液的特点，试验土箱夹角为30°扇柱体土箱，土箱采用有机玻璃制作，垂直高度50cm，径向长度50cm，垂向和径向外壁上分别贴有刻度，注液井半径为1cm。
1.2试验设计
试验土样取自江西龙南足洞矿区稀土。将土样自然风干，捣碎，过筛，测得土样的基本物理参数见表1所示。按预定容重(1.56g/cm3)将土样分层(5cm)装入土箱内，每装入一层土样，用木锤轻轻将其均匀捣实，将捣实后的层面抓毛再装入下一层。入渗水头分别设置为5cm、10cm、15cm、20cm、25cm，入渗时间要尽可能的长，一般为8～10h。
表1稀土试样的基本物理参数
1.3观测内容
注液开始后，用秒表计时，每隔10min记录1次数据，按照先密后疏的原则观测土体湿润体水平方向和垂直方向上的入渗距离随时间的变化。由于所用风干土颜色较浅，水分湿润后的区域颜色明显清晰，因此可以在土箱上描绘出二维坐标系内不同入渗时刻所对应的湿润锋位置及形状，测量横向湿润距离、垂向湿润距离和15°、30°、45°、60°、75°方向的最大湿润距离。为了保证试验数据的准确性，每个水头入渗重复试验3次，取其平均值作为试验结果。通过观测供水瓶以及收液量筒中水量，便可求出土样中水的累计入渗量。
2试验结果分析
2.1湿润体形状
单井注液的湿润体形状主要受土的质地和容重、水头高度、注液强度、注液时间和初始含水率等诸多因素的影响。湿润体形状是这些因素对土入渗能力和单孔注液强度的综合表现。为了方便观察湿润体，将实测的湿润体形状表示在二维坐标系内，如图2所示。由对称性可以看出，入渗的完整湿润体近似为半个椭球体。随入渗时间的增加，湿润体对称中心逐渐下移。分析认为，在入渗初期，土水势中主要是基质势作用在驱动湿润锋运移，水平方向上水分扩散速率快，土样湿润体的形状为长轴在横向的椭球体形状；随着入渗时间的增加，垂向的重力作用变大，重力势作用逐渐大于基质势，垂向的水分入渗速率相对增加，土样湿润体的形状开始向半球体发展。
考虑到一般情况，下面对湿润体形状是椭球体进行验证。假设土体湿润体为半个椭球体，且入渗点源在坐标系的中心，则过土体湿润体对称轴的任意剖面为一椭圆，该椭圆方程式为：
式中a代表湿润体的横向最大湿润距离，cm；b代表湿润体的垂向最大湿润距离，cm。表2为某次试验(h＝10cm)中不同入渗时间不同观测角度上(湿润锋和入渗圆点与水平方向之间的角度)湿润锋与入渗圆点之间的径向距离观测结果。
将10min实测的a＝7.5cm，b＝5.5cm代入式(1)，即：
设此时刻湿润体过对称轴的纵剖面上15°角对应的湿润锋横向坐标为x，则垂向坐标y可表示为：
y＝tg15°·x＝0.268x(3)
表2不同入渗时刻湿润锋距离点源的径向距离
将式(3)代入式(2)，得：
由式(4)解得x＝7.0cm，将x代入(2)可得y＝1.9cm，从而求得入渗10min时湿润体剖面15°角对应的湿润锋距点源的径向距离为：
同理可求得30°、45°、60°、75°角湿润锋距离点源的距离分别为6.81cm，6.22cm，5.87cm，5.58cm。依次计算30min、60min、180min、420min、600min湿润锋距离点源的距离，结果见表3。从表3可以看出计算的距离和实测的距离有一定的差异，但其相对误差都保持在10％以内。依次分析入渗水头为5cm、15cm、20cm、25cm的试验数据，与10cm结果一致。综合分析，可以认为具有入渗水头的单井注液所形成的湿润体形状是半个椭球体。
表3湿润锋实测值与计算值的误差分析
以方程α＝(a-b)/b计算椭球体扁率α，其中a和b为湿润体长半轴和短半轴的测量值，其值介于0和1之间，计算结果如图3所示。扁率反映了椭球体的扁平程度，α值愈大，椭球愈扁。分析可知，在入渗初始，湿润体扁率在0.3～0.4之间，随着入渗过程进行，湿润体扁率的逐渐减小，然后稳定在0.1左右，湿润体扁率的变化趋势也验证了随着入渗时间的增加，湿润体形状从半个椭球体逐渐向半球体发展。
2.2湿润体体积和入渗时间、注液量的关系
通过上述湿润体形状分析可知，单井注液离子型稀土湿润体的形状接近于半个椭球体，因此，土壤湿润体的体积v(t)可以用椭球体体积公式计算：
式中：a(t)为横向湿润距离，b(t)为垂向湿润距离。
掌握不同水头作用下，湿润体体积和入渗时间、累计入渗量之间的关系，无论是对离子型稀土原地浸矿的液井网设计还是注液量控制管理而言均具有重要的意义。用式(6)和试验数据计算不同入渗时刻的湿润体体积，其与累计入渗量的关系曲线如图4所示。由图4可以看出，在入渗量相同的条件下注液井内水头高度对于湿润体体积有一定影响，即供水量相同时，湿润体的体积随着注液水头的增加而减小。这也说明水头高度大时，湿润体的体积质量更大，含水率更高，饱和程度也更高，在一定程度上有利于离子型稀土原地浸矿的离子交换作用的发生。然而，相同累计入渗量的湿润体体积的差异不大，保持在10％以内。分析表明，湿润体体积v和累计入渗量w之间存在显著的线性关系，拟合方程为：
w＝0.326v，r＝0.9774(7)
式(7)表明土体湿润体内平均含水率增量基本上保持为一定值，即δθ＝0.326。设干土层的含水率为初始含水率θi，湿润体内的稳定含水率为饱和含水率θs，则δθ＝θs－θi＝0.326。此处含水率是体积含水率。
2.3湿润锋距离与入渗时间的经验关系
掌握不同水头入渗过程中土的湿润体特征值与入渗时间关系，是确定工程实践中注液井网参数的重要依据。恒定水头入渗的水平湿润距离a(t)和垂向入渗深度b(t)是湿润体二个重要的特征值，其入渗距离与时间关系曲线如图5和图6所示。
对横向入渗距离a(t)和垂向入渗距离b(t)与时间t的关系，用幂函数、指数函数和对数函数进行回归分析，结果表明幂函数能够比较精确地描述横向和垂向入渗距离与时间的关系，相关系数的平方值r2均在0.99以上。幂函数形式为：a(t)＝atb；b(t)＝ctd，参数a、b、c、d为湿润锋距离和入渗时间关系的拟合参数，二者通过湿润锋的观测易于获得，拟合结果见表4和表5，利用该函数能够显著的表示湿润锋的运移规律。
表4水平湿润锋和入渗时间关系拟合结果
表5垂向湿润锋和入渗时间关系拟合结果
3计算模型
3.1湿润体半径
通过试验可知，离子型稀土单井注液的影响范围是一个椭球体。基本假设：(1)注液井在椭球的中心，且初始时刻(t＝0)供水源半径等于注液井半径。(2)注液过程中，湿润体内部是饱和区，土中含水率为饱和含水率，渗透系数为饱和渗透系数，湿润体外部是干燥区，土中含水率为初始含水率。(3)土中毛管势忽略不计。(4)水分运动和湿润锋运移是从椭球中心沿着径向往四周扩散。在湿润峰上取楔形微元da，如图7所示，则有：
da＝2πrcosθ(rdθ)＝(2πcosθdθ)r2(8)
式中，da表示微元，量纲l2；r为湿润峰到椭球中心的径向距离，量纲l；θ为水平面与微元体的夹角，无量纲。
由式(8)可知，所取微元da是湿润锋到椭球中心的径向距离平方r2的比例项，由质量守恒定理可知，通过微元da的水流通量与微元面积(湿润锋到椭球中心的径向距离平方的比例项)成反比。即：
式中i为入渗率(水流通量)，量纲l/t；n为比例系数，l3/t；r表示微元上任意点到椭球中心的径向距离，量纲l。
根据达西定律和green-ampt模型，水流通量可以表示为：
式中k为饱和导水率，量纲l/t；h为总水头，量纲l。
联立(9)和(10)，可得：
整理后，得：
nr-2dr＝-kdh(12)
两边积分，得：
式中r0为注液井半径，量纲l；h为水头高度，量纲l。
求解方程(13)，得：
把(14)代入(9)，可得：
当r＝r时，湿润锋边界处的水流通量为：
由水量平衡原理和green-ampt入渗模型中的假定得到t时间内入渗量为水量i和土体湿润锋推进距离之间r的关系：
i＝(θs-θi)×r(17)
式中，i为累计入渗量，量纲l；θi为土体初始含水率，量纲l3/l3；θs为土体饱和含水率，量纲1。
入渗率和累计入渗量之间的转换关系如下：
上式表明，任意时刻的入渗率等于时段内累计入渗量对时间的导数，且等于土体湿润前后含水率变化值与湿润锋推进速度的乘积。因此通过试验记录土体湿润锋推进过程就可以计算得到土的入渗性能曲线。
由(16)和(18)，得：
重新整理上式，可得：
对上式进行积分，得：
求解方程(21)，得：
当sinθ＝0时，r为湿润体横向最大湿润距离，用ra表示，由(20)得：
对上式积分，可得：
当sinθ＝1时，r为湿润体垂向最大湿润深度，用rb表示，可得：
由(24)和(25)可以得到时间t与横向最大湿润距离、垂向最大湿润深度之间的关系。
若已知注液井半径、注液水头高度、初始含水率、饱和含水率、饱和渗透系数等参数，则可得到湿润锋距离r与时间t的关系曲线，如图8所示。由图看出，曲线的斜率在入渗初期比较大，随后趋于平缓。这是因为，在三维入渗中，水在土中渗透是向四周扩散，入渗初期湿润体体积小，渗透阻力小，随着湿润体体积越来越大，湿润体周围的渗透阻力逐渐增大，且小于基质吸力的增加量，故湿润锋推进的速度逐渐减小。
代入相应数据，利用matlab软件编程可以得到模型的计算值，湿润距离计算值与实测值的关系曲线如图9所示。由图可知，当用直线方程拟合不同入渗水头的湿润锋距离计算值与实测值曲线时，相关系数为r2＝0.9660，直线斜率为0.9906。说明二者具有较好的吻合度，即该计算模型具有很好的适用性。当实验室测得r，t时，此公式也可以反演饱和渗透系数，亦可以作为渗透系数计算与验证的一种方法。
3.2湿润体体积
湿润体体积是另外一个重要的特征值，图7中所取楔形微元体的体积可以由棱柱体公式求得，即
湿润体的体积可以表示为
式(25)是一个以时间为变量的隐函数，在实际中不方便应用与计算，可以运用数值反演的方法，把v和t的关系用图形呈现，如图10所示。
为了减少(26)中变量的数量，先对(22)进行无量纲化处理，两边除以r02可得：
由此可知，r/r0是h/r0、θ和kt/r0(θs－θi)的函数：
fr是一个关于r的隐函数，把(29)代入(27)，进行无量纲处理，可得到：
因为θ是这个积分的一个虚拟变量(哑变量)，故上式可以变成：
上式中，gv是一个关于v的隐函数，式(31)中的三维变量被用来计算v和t之间的几何关系。
图10表示了三维green-ampt模型入渗能力曲线。由图可以看出，当入渗时间较长的情况下，即当kt/r0(θs-θi)>10，直线的斜率均大于1，说明入渗率随着时间的增加而增加，也是三维入渗模型与一维模型最显著的差异。显然，充足供水条件下，随着湿润锋推移，湿润体的体积增量对入渗率的影响比水力梯度减少量对入渗率的影响要更大。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。