沉积磷石膏的物理力学特性试验研究

摘要：磷石膏是磷酸生产过程中的副产品，目前的综合利用率尚不足40%，大部分需要堆存存放。受地形限制和经济效益考虑，中国主要为湿法堆存的山谷型堆场。依托柳树箐磷石膏堆积坝，针对沉积磷石膏首先开展了密度、含水率、渗透、土水特征和颗分等物理性质试验，然后开展了三轴cu、蠕变及动三轴等力学特性试验。试验结果表明：①沉积磷石膏的干密度与埋深没有相关关系；②沉积磷石膏不具有自然分级现象，但具有明显的各向异性；③沉积磷石膏具有较高的摩擦角和抗液化能力，但其蠕变变形较大、渗透比降较小。上述工作为分析磷石膏堆积坝的坝体稳定性提供了基础，对现行磷石膏库的运行管理以及新建工程的设计具有重要的借鉴意义。

引言

磷石膏是湿法磷酸生产过程中的副产品，2018年，全国磷石膏产生量为7800万吨，且呈逐年增长的态势。磷石膏的主要成分是CaO和SO3，但含有一定量的氟化物和其它放射性物质，在中国通常按Ⅱ类一般工业固体废物处理，鉴于无害化处理成本较高，目前综合利用率尚不足40%，故大部分磷石膏需要堆存存放。按堆存场地的不同，可分为平地型、傍山型、山谷型和截河型堆场，在中国基本上是山谷型堆场。相比较干法堆存，湿法堆存经济优势明显，因而如地质条件为非碳酸盐岩地区，一般均采用湿排湿堆方式。随着中国磷肥工业的快速发展，本世纪初中国相继建设了几座湿法堆存的大型磷石膏库，例如云天化国际化工三环分公司的柳树箐磷石膏库堆积坝和富瑞分公司的杨家箐磷石膏库堆积坝，这两座坝设计坝高均超过100m、处于8度地震区，其安全性备受关注。

据统计，110多年(1901年一2013年)来，全世界有118座尾矿坝曾发生过破坏或溃坝事故，原因主要有地震、洪水漫顶、渗透破坏和基础失稳。尾矿库失事后会造成巨大的生态灾难和人员伤亡，近几十年来，

国内外对金属尾矿的沉积规律、物理力学特性及其稳定性开展了大量的研究，但对于与金属尾矿库相近的磷石膏库，一般**于在可研阶段采用人工制备样进行物理力学性质试验，并据此进行稳定性分析，尚未有人针对己运行若干年的大型湿法堆存磷石膏堆积坝开展过磷石膏沉积规律及其物理力学特性的专项研究。

本文依托柳树箐磷石膏堆积坝，首先进行了钻探取样，采用原状样开展了密度、含水率、渗透、土水特征和颗分等物理性质试验，在此基础上有针对性的选择原状样开展了三轴CU、蠕变及动三轴等力学特性试验。通过上述试验研究，总结了磷石膏的沉积规律、渗透特性、渗透破坏特性以及静动力特性，上述研究工作对研究和评估磷石膏库堆积坝的稳定性提供了基础数据，对现行磷石膏库的运行管理以及新建工程的设计具有重要的借鉴意义。

一、依托工程概况

1.1柳树箐磷石膏堆积坝堆存设计方案

由初期坝和堆积坝组成，设计总坝高约130m。

(1)初期坝

初期坝坝高约30m，采用土料填筑。上游坡面、坝底和下游坝脚设置堆石排水体，三者相连通构成整个堆积坝的主要排渗系统。

(2)堆积坝及其辅助排渗措施

采用上游式筑坝法，共20级子坝，顶宽6～9m，高度5m，堆积高度约100m。采用排渗管网作为辅助排渗方案，目前己在5级、9级子坝和13级子坝坝前120m范围内设置了井字形排渗管网。

1.2沉积磷石膏的钻探取样

钻孔平面位置见图l。2008年6月，堆积至5级子坝时，布设了9个取样钻孔，钻孔编号K1～K9，取原状样76件；2013年5月，堆积至13级子坝时，又布设了11个取样钻孔，钻孔编号K10～K20，取原状样112件，为比较子坝加高和磷石膏堆积过程中磷石膏物理力学性质的变化，在2，4级子坝K2孔和K6孔附近各布设了一个钻孔，钻孔编号分别为K17和K18。

1.3运行概况

柳树箐磷石膏尾矿库2005年开工建设，2006年1月投入运行，截至2013年5月己堆至13级子坝，尚有7级子坝即堆存至设计高程。鉴于磷石膏库地形、地质条件较好，具备扩容改造的条件，以提高堆存库容，减少堆存占地，节约土地资源。

本文主要对沉积磷石膏的物理力学特性进行了全面总结，限于篇幅，有关现状磷石膏库堆积坝的安全性评价及其加高可行性的研究将另文发表。

二、沉积磷石膏的物理力学特性

2.1物理特性

(1)干密度分布

图2给出了14个钻孔的取样深度和试验所得干密度的关系，图中UWL表示水位线上，(系钻孔期间的初见水位线，下同)，DWL表示水位线下。图3给出了水位线上下的饱和度分布图。由于磷石膏中的主要成分为CaSO4·2H2O，不同温度和烘烤时间对测定结果有一定影响，不能照搬现行的土工试验规范，本文磷石膏的含水率测定方法为55℃温度下烘培24h。

图2沉积磷石膏干密度与埋深的关系

Fig.2Variation of dry density with depth of depositlon PG

由图3，水位线上的磷石膏饱和度平均值为50.4%，处于非饱和状态，水位线以下的磷石膏饱和度平均值为85.0%，基本处于饱和状态，由于水位下降后磷石膏来不及排水固结，故而水位线上局部试样的饱和度较高。

由图2，水位线以上的干密度在0.98～1.67g/cm3之间，均值为1.30g/cm3；水位线以下的干密度在1.15～1.73g/cm3之间，均值为1.4g/cm3。可见磷石膏与一般的尾矿有所不同，磷石膏的干密度并不随埋深的增大而明显增大，但水位线以下的磷石膏干密度从统计意义上来看仍大于水位线以上的磷石膏干密度，这主要是由于水位线随库水位的变化反复升降而使得磷石膏排水固结所致。

室内击实得到的磷石膏**干密度一般在1.36～1.46g/cm3之间，从图2可以看出，自然沉积的磷石膏**干密度可达到1.73g/cm3，原因如下：与经典的土骨架不可压缩的理论不同，石膏本身可压缩，同时由于颗粒结构不稳定，击实试验过程中磷石膏结构被破坏，受夯击处下陷，四周鼓起，出现了类似于橡皮土的现象。而在现场条件下，石膏骨架被破坏后，会导致颗粒中的结合水渗出至孔隙内，变成孔隙水，排水固结后会使得磷石膏的孔隙比减小，干密度增大。

图4给出了相邻钻孔K2和K17(位于2级子坝河床部位)以及相邻钻孔K6和K18(位于4级子坝河床部位)干密度的对比图。K2钻孔的干密度在1.12～1.47g/cm3之间，均值为1.30g/cm3，K17钻孔的干密度在1.2～1.39g/cm3之间，均值为1.32g/cm3；K6钻孔的干密度在1.13～1.57g/cm3之间，均值为1.36g/cm3，K18钻孔的干密度在1.14～1.59g/cm3之间，均值为1.37g/cm3。可见，即使从统计意义上来看，磷石膏的干密度也并未随后续磷石膏的堆积而有较为明显的增大。

(2)级配分布

颗粒分析试验采用密度计法，制备悬液时不煮沸，不加六偏磷酸钠。图5给出了试验得到的级配包线、平均粒径d50、不均匀系数(1u和曲率系数Cc的分布图。

从图5(a)可见，磷石膏的粒径主要分布在0.005～0.075mm之间，总体上属于粉土，但可能由于矿石来源或生产工艺有所不同，局部属于粉砂～中砂。粒径分布范围比Blight和张超等的试验结果要宽一些。

图4子坝加高后沉积磷石膏的干密度变化

图5沉积磷石膏的平均粒径和级配分布

由图5(b)和5(c)，无论是水平向还是垂直向，磷石膏与金属矿山尾矿的“前粗后细，上粗下细”的自然分级现象不同，也即粗颗粒并不是沿埋深逐步减小或距离放浆口越远颗粒越细，其原因如下：①磷石膏颗粒粒径组成较为集中、均匀，主要以粉粒组(0.005mm<d≤0.074mm)为主，级配较差；②相比较金属尾矿，磷石膏的比重较小，磷石膏的比重一般为2.3～2.4，远小于金属尾矿的比重，例如铁尾矿的比重可达2.9；③放浆口随子坝高度不断增加而不断变动并向库尾延伸，造成沉积磷石膏的粒径变化不明显。

从图5(d)可见，不均匀系数Cu范围值1.61～21.5，平均值为4.18，曲率系数(1c范围值0.28～9.78，平均值为1.21，在统计的100多个试样中，属于级配不良土的占93%。这种级配特性决定了磷石膏具有较高的压缩性、渗透破坏型式表现为流土破坏。

2.2渗透特性

(1)渗透系数

影响渗透系数的主要因素是粒径大小、级配和孔隙比，因而磷石膏的渗透系数与粉土较为接近。由于孔隙比e减小，使得过水通道面积减小，渗透系数k也将减小，k与e呈正相关关系。对砂土，一般认为渗透系数k与e3/(1+e)的线性关系较好，图6给出了二者间的关系曲线，由于沉积磷石膏的不均匀系数变化较大，使得沉积磷石膏的渗透系数变化范围较大(平均值为10-4cm/s数量级)，二者问的线性关系较差。

图6沉积磷石膏渗透系数与孔隙比的关系

Fig.6Relationship between hydraulic conductivity and void ratio of depositlon PG

图7给出了水平与垂直向渗透系数比值的分布。沉积磷石膏的水平向渗透系数kh一般大于垂直向的渗透系数kv,kh/kv平均值约为2.86，这一点与成层分布的金属尾矿规律一致。造成沉积磷石膏水平向渗透系数大于垂直向渗透系数的原因是由于磷石膏具有明显的晶体结构，电镜扫描显示多为菱形和棱柱状形式(见图8)，在沉积过程中，由于扁平状磷石膏颗粒多呈水平排列，使得水平方向的透水性大于垂直方向的透水性，从而使磷石膏呈现明显的各向异性。

另根据中国有色金属工业昆明勘察设计研究院在杨家箐磷石膏堆积坝开展的现场渗透试验，kh/kv的平均值约为1.9。但张超等的室内试验显示，kh/kv的平均值约为0.46，也即垂直向渗透系数大于水平向渗透系数，与本文和现场试验结果恰恰相反。从磷石膏的微观结构来看，本文试验结果更为合理。

图7沉积磷石膏干密度与水平和垂直渗透系数比值的关系

(2)渗透变形

图9为磷石膏干密度为1.40g/cm3的水力梯度J与流速V的关系曲线，试验在水平管涌仪中采用水平方向的渗流形式进行。试验得到的临界坡降Jc=0.355，破坏坡降Jp=0.375。一般来说，对1，2级工程，渗透安全系数取为2.5，则允许出逸坡降为0.355/2.5=0.142，对3级以下工程，渗透安全系数取2.0，则允许出逸坡降为0.355/2.0=0.178。其允许比降与粉土一粉砂大致相同。但上述渗透变形试验是采用自来水进行的，自来水对磷石膏具有一定的溶蚀作用，而实际上磷石膏中残余磷、硫和氟酸，库水的pH值一般小于3(称之为酸性水)，在酸性水作用下，磷石膏不会发生破坏，上述试验结果是偏于保守的，但对非酸性水条件(例如特大暴雨)下的渗透稳定判断有一定的借鉴意义。

图9水力梯度J-流速v试验过程线(ρd=140g/cm3)

(3)土水特征试验

试验在5Bar的压力板仪中进行，环刀尺寸6.18cm。干密度分1.1，1.2和1.29g/cm3共3组，吸力范围0～500kPa。表l列出了含水率特征值，试验曲线见图10所示。

试验结果表明：①干密度对进气值没有明显的影响，不同干密度的试样的进气值大致在10kPa左右；②土样残余含水率随干密度的增加而减少，残余含水率约为饱和含水率的10%。上述特性与粉土一粉砂基本一致。

2.3静力力学特性

(1)三轴CU试验

由于沉积磷石膏的密度变化较大，而进行三轴试验需要若干原状样，为使试验结果具有较好的一致性，有针对性的选择平均干密度分别为1.2(1.2～1.21g/cm3)，1.28(1.27～1.28g/cm3)，1.4(1.39～1.41g/cm3)，1.5(1.49～1.51g/cm3)和1.58g/cm3(1.57～1.6g/cm3)的若干试样，进行了5组三轴CU试验。图11是为干密度为1.2和1.58g/cm3的三轴CU试验曲线。

从图11可以看出，磷石膏的应力应变关系曲线在低围压下表现为软化型，在高围压下表现为硬化型，与一般土类相似。但与一般粉土一粉砂不同的是，即使在较为疏松的状态下，磷石膏仍表现了较为强烈的剪胀，随密实度增大，剪胀作用愈发明显。

表2给出了不同干密度下的内摩擦角，图12给出了内摩擦角与干密度的关系曲线。随干密度的增大，内摩擦角也随之增大，且可采用幂函数较好地拟合。磷石膏的干密度由1.20g/cm3增加为1.58g/cm3，增加了32%，总应力摩擦角由34.1°增加为37.3°(根据拟合曲线求得)，增幅为9.4%，有效应力摩擦角由37.6°增加为38.8°，增幅为3.2%，由于随围压的增大，孔压也明显增大，故有效摩擦角增幅较之总应力摩擦角要小。另外较之于干密度增幅，摩擦角的增幅并不显著，表明即使较为疏松的磷石膏仍具有较高的强度指标，这也表明磷石膏堆积坝的稳定性较高。

(2)蠕变(次固结)变形试验

磷石膏的蠕变变形试验在侧限压缩仪中进行，试验状态相当于K0固结。试样直径61.8mm，高度20mm，试样干密度为1.30g/cm3，对试样饱和后分别开展了上覆压力ρ为100，200，400，800kPa的试验，

试验从2012年8月27日开始，试验己进行了1年多，试验结果见图13所示。

从图13中可以看出：上覆荷载越大，磷石膏的蠕变变形也越大，荷载施加1年后，磷石膏的蠕变变形仍非常显著，尚未达到稳定状态，这也是磷石膏堆积坝后期变形较大的原因，原型观测资料表明，在5级子坝河床部位的表面沉降量已经达到3.1m，且尚未完全稳定。

按时间对数法，可求得各级荷载下的主次固结变形量如表3所示。试验结果表明，在100～400kPa上覆荷载作用下，在试验时间范围内蠕变(次固结)变形是主固结变形的1.8～3.1倍，当然由于蠕变变形尚未完成，实际的蠕变变形应更大。对土体而言，发生蠕变的原因是由于土体在主固结完成之后，土体中仍有微小的超静孔隙压力存在，驱使水在颗粒问流动，一般来讲土体的次固结远小于主固结变形；对磷石膏而言，其渗透系数在10-4cm/s数量级，远大于黏性土，但其却发生了极为显著的次固结变形，其原因在于磷石膏晶体结构发生了压缩、破坏，接触点晶格发生歪曲和变形，而破坏后晶格之间的重新排列、调整到最后趋于相对静止需要相当长的时间才能完成。

2.4动力力学特性

试验设备采用英国GDS公司进口的电机控制动三轴试验系统，试样直径39.1mm，高度80mm。

(1)动模量和阻尼比

同样由于自然沉积的磷石膏密度变化较大，为此根据物理性质试验结果，选择两种平均干密度1.34g/cm3(变化范围1.33～1.35g/cm3)和1.45g/cm3(变化范围1.44～1.46g/cm3)进行试验。

Hardin-Dmevich建议的动剪切模量G、阻尼比与剪应变幅值γd的关系式如下：

式中k1为参数，表示动剪切模量的衰减或阻尼比的增长速率；λmax为**阻尼比；Gmax，γd分别为**动剪切模量和归一化的动剪应变，表示为

式中k2,n为参数；σm为球应力；Pa为标准大气压；vd为动泊松比；εa为归一化的动应变，表达为

图14给出了动剪切模量、阻尼比与归一化动应变的关系曲线，另外图中还给出了式(1)的拟合曲线以及Seed等给出的砂样的上下边界线，图例中，σ2表示围压，Kc表示固结应力比。

从图中可以看出：①式(1)可较好地描述磷石膏的动应力-动应变试验曲线，表明采用等价黏弹性模型进行循环荷载作用下的分析是可行的；②图14(a)中干密度为1.45g/cm3的拟合线位于干密度为1.34g/cm3的拟合线上方，图14(b)中则位于下方，表明密度越大，动弹模越大、阻尼比越小；③图14(a)中，两种干密度的拟合线基本位于Seed等给出的边界线上方，而图14(b)中则基本处于边界线中间，表明相比较砂样，磷石膏动弹模较大，会导致磷石膏堆积坝的动力反应较大，但由于阻尼比也较大，这样又会削弱坝体的动力反应，二者的相互影响下，磷石膏堆积坝坝体的动力反应将不会过于强烈，这对磷石膏堆积坝的抗震稳定性是有利的。

(2)动强度

选择两种平均干密度为1.12(变化范围1.1～1.14g/cm3)和1.306g/cm3(变化范围1.29～1.3lg/cm3)进行试验，破坏标准为总应变达到10%。

图15给出了动剪应力比τd/σ0′与破坏振次Nf的关系曲线图，其中σ0为振前试样45°面上的有效法向应力，表达为：σ0=(Kc＋1)σ2′/2，Kc为固结比。从试验结果可以看出，沉积磷石膏的动强度与其它土体相似，表现为围压和固结应力比与动剪应力比呈负相关关系。

为判别沉积磷石膏的抗液化能力，假定抗震设计烈度为8度，即等效振次Ⅳ可取为30。首先由式(4)

所示的幂函数关系式得到振次为30时各个围压和固结比下的动剪应力比：

τd/σ0′=aNr-b(4)

然后可拟合求得动剪应力比与围压和固结应力比的关系式：ρa=1.12g/cm3：τd/σ′0=0.5395－0.01σ3/pa0.037Kc；ρd=1.306g/cm3：τd/σ0′=0.589－0.008σ3/ρa0.039Kc。

从上述拟合关系式可见，密实度提高后，动剪应力比提高，表明抗液化能力也提高。但即使是在低密度下，其动剪应力比较之同类的粉土或粉砂也大出许多，表明磷石膏具有较高的抗液化能力。

三、结论

依托柳树箐磷石膏堆积坝，在钻探取样工作的基础上，首先开展了物理性质试验，然后开展了静动力力学特性试验。通过上述试验研究，得出如下结论：

(1)沉积磷石膏总体上属于级配不良的粉土，局部属于粉砂一中砂，无自然分级现象。其干密度和粒径变化随埋深或距放浆口距离的变化不明显。

(2)沉积磷石膏水平方向的渗透系数大于垂直方向的渗透系数，呈现明显的各向异性。

(3)与土体颗粒不可压缩不同，磷石膏的晶体结构会发生压缩破坏，具有较大的压缩性，其次固结变形量远大于主固结变形量。

(4)沉积磷石膏的静动强度较之同等密实度下的粉土、粉砂要高。