Molecular dynamics simulation of adsorption of polyacrylic acid on surface of dicalcium silicate
-
摘要: 采用分子动力学(MD)方法,模拟研究了聚合物抑制剂聚丙烯酸与β-硅酸二钙(110)晶面的相互作用行为。结果表明:聚丙烯酸与β-硅酸二钙(110)晶面间结合能为负值,在模拟范围内结合能随聚丙烯酸聚合度的增加而增大,聚丙烯酸分子能较好地吸附在β-硅酸二钙(110)表面。通过体系中各种相互作用、径向分布函数以及吸附构象图分析,发现聚丙烯酸主要是通过库仑静电作用吸附在β-硅酸二钙(110)晶面上。聚丙烯酸中不同位置羧基的分子动力学模拟结果差别也很大,链端羧基波动扭转比链中端羧基强度大得多,表明中端羧基与β-2CaO·SiO2晶体的吸附比链端更牢固因而更能有效地吸附在晶面上,达到抑制熟料二次反应的作用。Abstract: The interactions between polyacrylic acid (PAA) and β-dicalcium silicate (110) crystal face were simulated by molecular dynamics(MD).The results show that the binding energy for the polyacrylic acid polymer with β-dicalcium silicate (110) crystal is negative and increases with the increasing in the degree of polymerization of polyacrylic acid within simulation range, which indicate that the polyacrylic acid can be adsorbed well on the β-dicalcium silicate crystal face. The analysis of various interactions and pair correlation functions of all systems and adsorption conformation chart indicate that binding energies are mainly determined by coulomb interaction. The results of molecular dynamics simulation of carboxyl groups at different positions of polyacrylic acid also differs greatly. The strength of torsion fluctuation of the chain end carboxyl is larger than that of the mid-chain. Both of those indicate that the mid-chain carboxyl’s adsorption with β-2CaO·SiO2 crystals is firmer than the chain end carboxyl’s adsorption with β-2CaO·SiO2 crystals, and more effective. Thus, clinker secondary reaction can be inhibited.
-
坚持绿色、循环、可持续的发展路线,社会才能长期稳定的发展,经济发展必须重视自然的保护问题[1-2].然而,人类的活动还是有很多方面影响着自然平衡.因矿山地下开采空区未进行处理造成巷道变形、地表沉降、房屋开裂等对自然环境造成破坏的例子比比皆是,国内外对此问题进行的研究也不尽其数[3-4].开采环境不同,影响模式不同,研究价值及意义也就不同.
宁都硫铁矿是当地一家非常重要的矿山企业,该矿山一直采用地下开采,开采方法为浅孔房柱采矿法,矿房长60 m,开采已形成+130 m、+160 m、+190 m 3个中段,最上中段+190 m离地表标高+240 m的高程差仅为50 m.井下开采空区对应地表为居民楼、农田及公路,矿山的开采必须以地表建筑物不破坏为前提、得到百姓的默认.矿山开采至今已形成较大采空区,且采空区仅留房柱对上部围岩进行支撑.矿山开采是否造成空区失稳是该矿山工作人员一直担忧的问题.基于此,通过理论分析及现场监测、预测研究,两者对比分析、综合评判研究矿山开采的稳定性[5-6],为矿山今后的开采提供指导意见及建议.
1 上覆岩层移动变形分析
1.1 上覆岩层组合及关键层
空区上覆岩层中对岩层移动起主要控制作用的岩层为关键层[7-8],根据其定义与变形特征,假定有n层同步协调变形,则最下部的岩层为关键层.同时,根据关键层变形特征及支承特征,可得:当(qn+1)1<(qn)1时,则此时的n层为关键层.由此可知n=1,即每一个岩层组合中最下面的那个岩层为关键层.
宁都硫铁矿上覆岩层较薄,顶板岩层由下而上为较硬的石英砂岩、黏土砂土.由此可知,较硬的石英砂岩即为关键层,对上覆岩层的破坏起主要控制作用,+190 m中段石英砂岩岩层厚度约15~20 m.
1.2 岩层破坏极限跨距理论
该矿山岩层可看作为水平矩形岩层,采用托板理论进行分析[9-10],假设共有n层岩层,取其中任一i层进行研究,岩层长为L,宽度为S,厚度为hi,hi比L及S小.宁都硫铁矿特点为薄矿层开采,根据弹性薄板的假设原理可以将岩层近似为薄板[11-12],进而可以推导出第i层的物理与几何方程,物理方程:
$$ \left. \begin{array}{l} {\sigma _x} = - {E_i}z/1 - {\mu _i}^2\left( {\partial _{wi}^2/{\partial _x}^2 + {\mu _i}\partial _{wi}^2/{\partial _y}^2} \right)\\ {\sigma _y} = - {E_i}z/1 - {\mu _i}^2\partial _{wi}^2/{\partial _y}^2 + {\mu _i}\partial _{wi}^2/\partial \\ {\tau _{xy}} = - {E_i}z\partial _{wi}^2/\left( {1 + {\mu _i}} \right){\partial _x}{\partial _y} \end{array} \right\} $$ (1) 几何方程:
$$ \left. \begin{array}{l} {\varepsilon _x} = - z\partial _{wi}^2/{\partial _x}^2\\ {\varepsilon _y} = - z\partial _{wi}^2/{\partial _y}^{2}\\ {\gamma _{xy}} = - 2z\partial _{wi}^2/{\partial _x}{\partial _y} \end{array} \right\} $$ (2) x与y方向弯矩:
$$ \left. \begin{array}{l} {M_x} = - {D_i}\left( {{\partial ^2}/\partial _x^2 + {\mu _i}\partial _w^2/\partial _y^2} \right)\\ {M_y} = - {D_i}\left( {{\partial ^2}/\partial _y^2 + {\mu _i}\partial _w^2/\partial _x^2} \right) \end{array} \right\} $$ (3) 其中:Di为弯曲刚度;Ei为弹性模量;μi为泊松比;Wi(x,y)为第i层岩层扰度.n为托板的弯矩:
$$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{M_n}\left( {x, y} \right) = A{q_n}\\ A = \frac{{7(3{x^2} - {L^2}/4){{({y^2} - {S^2}/4)}^2} - 7{\mathit{\mu }_\mathit{n}}{{{\rm{(}}{\mathit{x}^{\rm{2}}}{\rm{ - }}{\mathit{L}^{\rm{2}}}{\rm{/4)}}}^{\rm{2}}}{\rm{(3}}{\mathit{y}^{\rm{2}}}{\rm{ - }}{\mathit{S}^{\rm{2}}}{\rm{/4)}}}}{{{\rm{2(}}{\mathit{L}^{\rm{4}}} + {\mathit{S}^{\rm{4}}} + {\rm{4}}{\mathit{L}^{\rm{2}}}{\mathit{S}^{\rm{2}}}{\rm{/7)}}}} \end{array} $$ 进而求得:
$$ {M_{{\rm{max}}}} = 7{q_n}({L^2}{S^4} + {\mu _n}{L^4}{S^2})/128({L_4} + {S_4} + 4{L^2}{S^2}/7) $$ 最后,求出n托板的等效跨距Ln:
$$ {L_n}\sqrt {\frac{{{\rm{ - }}{F_\mathit{n}} + \sqrt {F_n^2 + 4{G_\mathit{n}}{K_\mathit{n}}} }}{{2{G_\mathit{n}}}}} $$ (4) 其中:${F_n} = 21{q_n}{S^4} - \frac{{-256{S^2}{\mathit{h}_\mathit{n}}{\mathit{\sigma }_\mathit{n}}}}{{\rm{7}}}$;Gn=21qnμnS2-64h2nσn;Kn=64hn2S4σn
依据式(4)Ln的计算结果判定托板的稳定性,当Ln计算值小于采空区跨度时,托板受到破坏,相反,Ln计算值大于采空区跨度则稳定.依据以上岩层破坏极限跨距理论,前人总结得出了关键层的厚度与关键层极限跨度的关系,根据研究成果[13-14],矿房长度不一,其所对应值不一,当矿房长度在50 m与60 m时关系如曲线图 1所示.
综合分析,矿山矿房长度60 m,关键层为石英砂岩,厚度为15~20 m.由图 1可以看出当关键层厚度在15~20 m时的极限跨度为25~40 m,即空区跨度小于25 m才能保证稳定.宁都硫铁矿开采区域内空区的跨度主要在8~10 m区间,且区域内矿房条件及结构基本相似,由此可见开采区域空区岩层稳定.
2 地表沉降灰色系统预测分析
灰色理论中最为常用的灰色预测模型是等时距GM(1,1)模型[15-16],先假设x(0)={x(0)(t1),x(0)(t2),…,x(0)(tn)}为观测得到的原始沉降量数据序列,时间步长k=ti+1-ti,r次累加后,得到新数列:x(r)={x(r)(1),x(r)(2),…,x(r)(n)},
累加一次得到等时距预测模型微分方程:
$$ \frac{{{\rm{d}}{\mathit{x}^{\left( 1 \right)}}}}{{{\rm{d}}\mathit{t}}} + a_x^{\left( 1 \right)} = \mu \left( {a, \mu \;为常数} \right) $$ 得出时间响应函数:
$$ {\mathit{x}^{\left( 1 \right)}}\left( t \right) = \left[{{\mathit{x}^{\left( 1 \right)}}\left( 0 \right)-\frac{\mu }{a}} \right]{e^{ - a\left( {t - 1} \right)}} + \frac{\mu }{a} $$ (5) 令$Y = {\left[{{x^{\left( 0 \right)}}\left( 2 \right), {x^{\left( 0 \right)}}\left( 3 \right), \cdots, {x^{\left( 0 \right)}}\left( n \right)} \right]^{\rm{T}}}$
$$ \begin{array}{l} B = \left[\begin{array}{l} - \frac{1}{2}{\left[{x_1^{\left( 1 \right)}\left( 1 \right) + x_1^{\left( 1 \right)}\left( 2 \right)} \right]^1}\\ - \frac{1}{2}{\left[{x_1^{\left( 1 \right)}\left( 2 \right) + x_1^{\left( 1 \right)}\left( 3 \right)} \right]^1}\\ \;\;\;\;\;\;\;\; \cdots \;\;\;\;\;\;\;\;\;\;\; \cdots \\ - \frac{1}{2}{\left[{x_1^{\left( 1 \right)}\left( n \right) + x_1^{\left( 1 \right)}\left( {n-1} \right)} \right]^1} \end{array} \right] \end{array} $$ 由最小二乘法[17-18]得到:[aμ]T=[BTB]-1BTY,求出a、μ值,代入式(5)中进行灰色预测[19-20].
在后续3.2小节中具体分析了监测方法,对矿区地表进行了9次(15个月周期)的沉降监测,都取得了有效的监测数据.监测结果显示监测点MB5为最大沉降点,并归纳得到MB5沉降点7次的实测沉降数据(原始序列).同时,利用灰色预测模型进行预测(图 2),预测值与实测值进行对比,结果如表 1所列.可以看出,通过灰色系统预测得到的沉降值与实际监测沉降量基本一致,灰色预测模型预测地表沉降具有一定的可靠性.
表 1 实测与预测累积沉降值Table 1. Measured and predicted cumulative settlement value时间/(年-月) 序号 沉降差/mm 实测值/mm 预测值/mm 2013-10 1 0 -8.4 / 2014-01 2 -2.4 -10.8 -14.4 2014-03 3 -7.8 -18.6 -18.8 2014-05 4 -6.7 -25.3 -24.6 2014-07 5 -7.3 -32.6 -32.2 2014-10 6 -12.3 -44.9 -42.0 2015-01 7 -7.0 -51.9 -54.8 注:“/”指未进行预测. 3 现场监测预测研究
3.1 地表沉降监测方案
地表沉降能够反映地压的变化规律,地表沉降值可以通过长期的监测获取.依据矿山井上井下对照图及空区图,预估沉降范围,最后选取沿公路及公路两边辐射的监测线路;基准点布置于预估沉降区域范围外,根据矿区实际环境,设计了合理的监测方案,方案及工作量如表 2、图 3中描述.观测采用精密水准仪,配套水准尺.
表 2 地表沉降监测方案Table 2. Surface settlement monitoring program内容 监测点数量/个 监测点位置 周期/月 工作量/次 沉降监测 32 5个基准点(BM1~BM5),
8个一级沉降观测点(MB1~MB8),
19个二级沉降观测点(MB1-1~MB9-3)15 256 3.2 结果及分析
参考有关研究者及类似监测案例的处理分析方法[21-22],对监测周期内有效的监测数据进行处理分析,绘制沉降量—时间曲线图,同时,也得到了监测周期内各监测点的最终沉降量, 如图 4、表 3.
表 3 各监测点最终沉降量Table 3. Final settlement of monitoring points监测点 沉降量/mm 监测点 沉降量/mm 监测点 沉降量/mm MB1 0.00 MB1-2 -0.88 MB6-1 -31.05 MB2 -7.10 MB2-1 -0.65 MB6-2 -32.25 MB3 -13.55 MB2-2 -3.90 MB7-1 -40.55 MB4 -29.25 MB3-1 -19.90 MB7-2 -43.20 MB5 -51.88 MB3-2 -18.50 MB8-1 -8.05 MB6 -44.00 MB4-1 -19.89 MB8-2 -2.85 MB7 -29.75 MB4-2 -25.00 MB9-1 -18.10 MB8 -17.85 MB5-1 -30.85 MB9-2 -13.45 MB1-1 -1.80 MB5-2 -32.25 MB9-3 -15.60 由图 4与表 3显示,在监测周期内,监测点基本上都发生了下沉现象,下沉级别在毫米至厘米范围内,由表 3可以看出沉降量最小为0,最大为51.88 mm;一级沉降监测点的沉降值与预想一致,普遍较大;同时还可以看出MB4、MB4-1、MB4-2、MB5、MB5-1、MB5-2、MB6、MB6-1、MB6-2、MB7、MB7-1、MB7-2的沉降量属于大值,这些监测点刚好位于公路两侧附近,也是井下采空区的正上方,产生较大沉降量的原因必然与井下矿体的开挖有关.最大沉降处为农田,在可控范围之内,但在今后的开采中矿山必须重视.对于沉降较为严重的区域,不应进行开采,并及时采取支撑、充填等措施进行防护.
4 结论
1)采用理论分析宁都硫铁矿空区岩层的稳定性,判定得出宁都硫铁矿上覆岩层现开采阶段稳定.
2)利用灰色预测模型,结合地表监测数据预测沉降量,结果显示预测沉降量与实际监测值具有一致性,说明灰色预测模型能应用于地表沉降的预测.
3)根据矿山实际情况,进行了长期的地表沉降监测,结果显示大部分监测点都具有沉降数据,下沉级别控制在毫米至厘米级别,最小为0.00 mm,最大为51.88 mm,最大沉降处为农田,沉降值在可接受范围内,但矿山在今后必须跟踪监测,引起重视.
4)矿山开采对地表及矿山本身具有一定的影响,根据分析及监测现今暂未造成失稳.但在今后的开采中开采稳定性问题必须引起重视,有必要时采取充填等防范处理措施.
-
表 1 聚合物分子与β-2CaO·SiO2(110)晶面相互作用的结合能和非键作用能(eV)
Table 1 Binding interaction energies and non—bond interaction energies between polymers and(110)surface of β-2CaO·SiO2(eV)
能量形式 不同链长的PAA 8 10 12 14 20 Etotal -6691.299 -7922.779 -6691.311 -7924.557 -7932.705 Esurface -6680.930 -7913.624 -6680.969 -7907.592 -7906.279 Epolymer -4.442 -2.341 -3.090 -2.682 -4.020 Einter -5.927 -6.813 -7.252 -14.283 -22.406 Ebinding 5.927 6.813 7.252 14.283 22.406 Enon-bond -6.882 -5.302 -10.145 -13.740 -21.473 Evdw 2.075 5.367 -2.601 -0.923 1.323 Ecoulomb -8.944 -10.660 -7.532 -12.808 -22.777 -
[1] 陈滨,李小斌,徐华军,等. 氧化铝熟料溶出过程二次反应的热力学讨论[J].北京化工学报(自然科学版),2007,34(2):189-192. http://www.cnki.com.cn/Article/CJFDTOTAL-BJHY200702019.htm [2] 袁艺,项阳,黄芳,等. 铝酸盐烧结熟料溶出二次反应浅析[J].轻金属,1999(12):18-21. http://www.cnki.com.cn/Article/CJFDTOTAL-QJSS199912005.htm [3] LI X B, ZHAO Z, LIU G H, et.al .Behavior of calcium silicate hydratein aluminate solution [J]. Transactions of Nonferrous Metals Society of China,2005,15 (5):1145-1149 http://cn.bing.com/academic/profile?id=2349545729&encoded=0&v=paper_preview&mkt=zh-cn
[4] 陈红武,周宗科.烧结法熟料溶出条件对二次反应影响分析[J].轻金属,2001(8):21-25. http://www.cnki.com.cn/Article/CJFDTOTAL-QJSS200108003.htm [5] 陈滨,张淑英.烧结法熟料高浓度溶出过程二次反应的抑制[J].湖南工业大学学报,2011,25(3),8-12. http://www.cnki.com.cn/Article/CJFDTOTAL-ZZGX201103005.htm [6] 张建.碱石灰铝土矿熟料高浓度溶出过程中二次反应的研究[D].长沙:中南大学,2008 http://cdmd.cnki.com.cn/Article/CDMD-10533-2008164674.htm [7] 任根宽.二次反应对熟料中氧化铝溶出率的影响分析[J].轻金属,2007(4):16-18. http://www.cnki.com.cn/Article/CJFDTOTAL-QJSS200704003.htm [8] 李太昌.二次反应抑制剂及其添加工艺技术研究[J].有色金属(冶炼部分),2002(1):26-28. http://www.cnki.com.cn/Article/CJFDTOTAL-METE200201007.htm [9] 任根宽.二次反应抑制剂及其添加工艺技术研究[J].轻金属,2008(5):16-18. http://www.cnki.com.cn/Article/CJFDTOTAL-QJSS200805006.htm [10] 张程忠,于海燕,张立强,等.腐殖酸钠抑制氧化铝熟料溶出二次反应机理初探[J].化工学报,2008,59(2):526-530. http://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ200802042.htm [11] 张程忠.抑制熟料溶出二次反应添加剂的研究[D].沈阳:东北大学,2007. http://cn.bing.com/academic/profile?id=2609144788&encoded=0&v=paper_preview&mkt=zh-cn [12] YU H Y,PAN X L,DING T T,et.al.Electrochemical study on adsorption behavior of surfactants at β-2CaO?SiO2/NaAlO2 interface[J].Transactions of Nonferrous Metals Society of China,2013,8(23):214-219. http://cn.bing.com/academic/profile?id=2044104895&encoded=0&v=paper_preview&mkt=zh-cn
[13] Peng Xuan.Molecular simulations of the purification of toxic benzene gas on single-walled carbon nanotubes[J]. Acta Physico-Chimica Sinica,2014,30(11):2000-2008.
[14] 倪哲明,夏明玉,施炜,等.糠醛在Pt(111)表面的吸附和脱碳反应[J].物理化学学报,2013, 29 (9):1916-1922. http://www.cnki.com.cn/Article/CJFDTOTAL-WLHX201309010.htm [15] MEI Q Q,HOU M Q,NING H, et.al.Microstructure and intermolecular interactions of [Bmim][PF6]+Water+ Alcohol Systems: A Molecular Dynamics Simulation Study[J].Acta Physico-Chimica Sinica,2014, 30 (12):2210-2215. http://cn.bing.com/academic/profile?id=2369552101&encoded=0&v=paper_preview&mkt=zh-cn
[16] 于海燕,王波,潘晓林,等.聚乙二醇在硅酸二钙表面吸附的分子动力学模拟[J].化工学报,2013,64(3):943-948. http://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ201303023.htm [17] Sun H,Ren P,Fried J.R.The COMPASS Force Field:Parameterization and Validation for polyphosphazenes.Computational and Theoretical Polymer Science, 1998,8(1/2):229 http://cn.bing.com/academic/profile?id=2024955711&encoded=0&v=paper_preview&mkt=zh-cn
[18] K.H.Jost. B.Ziemer R.Seydel. Redetermination ofthe Structure of β-Dicaleium Silicate[J],Acta Cryst. 1997.B33:1696-1700.
[19] 黄玉成,胡应杰,肖继军,等.TATB基PBX结合能的分子动力学模拟[J].物理化学学报,2005,21(4):425-429. http://www.cnki.com.cn/Article/CJFDTOTAL-WLHX200504015.htm [20] 张曙光,王风云,雷武,等.水溶性聚合物与硬石膏晶体相互作用的分子动力学模拟[J].化学学报,2007,65(20):2249-2256. http://www.cnki.com.cn/Article/CJFDTOTAL-HXXB200720007.htm [21] 杨小震,分子模拟与高分子材料.北京:科学出版社,2002. -
期刊类型引用(0)
其他类型引用(5)