Numerical Study on Evolution Law of the Mechanical Properties of Granite under Ultra-High Temperature Conditions
-
摘要: 为提高深层干热岩破碎效率,以干热岩地层常见岩性-花岗岩为研究对象,提出超高温条件下花岗岩力学性质演化规律研究。鉴于室内实验方法的局限性,主要采用数值模拟和理论计算等研究方法,建立了高温花岗岩巴西圆盘劈裂数值模型,探究了硬质花岗岩在温度-压力联合作用下岩石应力场分布特征与扰动机制。研究发现:在单轴径向压缩载荷下,高温花岗岩各向应力分布受到显著扰动,岩石沿横轴拉伸应力、两端载荷施加位置附近压缩应力和剪切应力均显著减小,导致裂纹越容易扩展。温度越高,岩石损伤越明显,其岩石抗拉强度越低。研究结果揭示了温度作用下花岗岩力学性质演化规律,可为深层干热岩资源的高效开发提供理论依据。Abstract: To improve the rock-breaking efficiency of deep dry hot rock, taking granite, a common lithology of dry hot rock formation, as the research object, the changes in the mechanical properties of granite under ultra-high temperature was proposed. In view of the limitations of laboratory experimental methods, the Brazilian disc splitting numerical model of high-temperature granite was established by using numerical simulation and theoretical calculation methods, and the distribution characteristics and disturbance mechanism of stress field of hard granite under temperature-pressure conditions were investigated. It is found that under uniaxial and radial compression load, the stress distributions of high-temperature granite get significantly disturbed, leading to lower tensile stress along the horizontal direction, smaller compression and shear effect near the loading points, and easier fracture conduction. The higher the temperature, the more obvious the rock damage and the lower the tensile strength of the rock. The results reveal the evolution law of mechanical properties of granite under the effect of temperature, which can provide theoretical basis for efficient development of deep dry hot rock resources.
-
Keywords:
- hot dry rock /
- granite /
- temperature /
- Brazilian splitting /
- tension /
- stress field
-
随着工业化的迅猛发展,全球正面临着日益严峻的能源危机和环境污染问题。氢能是一种高效、清洁的新能源,可被广泛用于氢能发电、氢动力汽车和氢燃料电池等领域[1]。电催化分解水作为一种可持续地利用丰富水资源制氢的方法备受关注[2-4],通常包含2个半反应:阴极的析氢反应(HER)和阳极的析氧反应(OER),其中OER涉及四电子转移,从而制约了电催化分解水制氢的效率[5-6]。Ir和Ru化合物已被证明具有优秀的OER性能,然而这些贵金属催化剂储量稀少,价格昂贵,无法大规模制备和应用[7-8]。因此,利用非贵金属催化剂取代贵金属催化剂是研究者追求的目标。金属硫化物具有较高的电导率、优良的电化学稳定性和丰富的活性位点,已被广泛用作OER催化剂,例如Fe-Ni3S2/NF[9]。金属硫化物一般通过水热或溶剂热法、湿化学法和化学气相沉积法等方法合成[10-11],不可避免地需要高温或高压等苛刻条件,同时伴有H2S等有害气体的释放。
本研究采用电沉积技术在镍网(NF)上原位制备氧化镍-硫化铁(NiO-FeS)复合材料,研究电催化分解水的产氢效率[12-14]。电沉积法具有操作简单、可控性强等特点[15-17]。采用两步电沉积法制备的NiO-FeS@NF呈现独特的二维纳米片状结构,极大提升了电催化OER性能。
1. 实验部分
1.1 试剂及仪器
试剂:硫脲(CH4N2S)、硫酸镍(NiSO4·6H2O)、硫酸亚铁(FeSO4·7H2O)、盐酸(HCl)和氢氧化钾(KOH),所有试剂均为分析纯。
仪器:扫描电镜(SEM,Sigma 300,德国Zeiss)、透射电镜(TEM,Talos F200X,美国FEI)、能量色散X射线谱(EDS,X-Max 50,英国Oxford)、X射线光电子能谱(XPS,Lab220i-XL,美国Thermo Fisher)、电化学工作站(CHI660E,上海辰华)。
1.2 NiO-FeS@NF的合成
用盐酸、乙醇和去离子水依次超声清洗镍网(NF,1 cm×0.5 cm)约20 min,干燥待用;以NF为工作电极、Ag/AgCl电极为参比电极、铂线为对电极进行电沉积,电解质为0.1 mol/L FeSO4·7H2O和0.1 mol/L CH4N2S的混合溶液,沉积电位为-0.85 V vs Ag/AgCl,沉积时间为15 min,沉积完成后用去离子水清洗并干燥;采用FeS@NF为工作电极,电解质是0.1 mol/L NiSO4·6H2O水溶液,沉积电位为-1 V vs Ag/AgCl,沉积时间为15 min;将沉积产物用去离子水清洗并干燥,得到NiO-FeS@NF复合材料。
1.3 电化学测试方法
采用电沉积法使NiO-FeS在NF表面原位生长,保证了复合材料和基底之间的电子传递,因此NiO-FeS@NF可以直接用作电催化水解的工作电极。在1 mol/L KOH电解质中测试NiO-FeS@NF复合材料的线性扫描伏安(LSV)曲线。泡沫镍在电催化过程中易被氧化,出现较强的氧化还原峰。为了准确测试材料的电催化OER性能,避免镍的氧化还原反应,LSV的扫描电位采用从高电位到低电位的方式,即从0.8 V到0 V,扫描速率为5 mV/s。
在10-2~105 Hz的频率、0.45 V vs Ag/AgCl下测试材料的电化学阻抗谱(EIS),测试电解质为1 mol/L的KOH溶液。
依据LSV曲线,利用塔菲尔公式分析过电位。塔菲尔公式:
$$ \eta = b\lg j + a, $$ (1) 其中,b是Tafel斜率,j是电流密度,η是过电位。
在25 ℃下,可逆氢电极(RHE)与Ag/AgCl电极电位的换算关系:
$$ {V_{{\rm{RHE}}}} = {V_{{\rm{Ag/AgCl}}}} + 0.059{\rm{pH}} + 0.197。 $$ (2) 2. 结果与讨论
2.1 NiO-FeS@NF复合材料的表征结果
2.1.1 NiO-FeS@NF的表面形貌
分别对FeS@NF、NiO@NF以及NiO-FeS@NF复合材料的形貌采用SEM形貌观察(图 1),FeS@NF为错落交织的二维片状结构,NiO@NF为堆积的片状结构,而NiO-FeS@NF复合材料兼具这两种形貌,主要由二维纳米片形成的团簇组成,表明FeS与NiO已成功复合。
复合材料的TEM和HRTEM表征结果(图 1D)表明:NiO-FeS@NF呈现明显的二维纳米片结构,且高度透明,说明所得材料具有超薄的特点。超薄二维片状材料具有大的比表面积和丰富的暴露活性位点[18],有利于催化效率的提高[19-20]。从图 1E中可以清晰观察到晶格间距为0.210、0.217 nm的2种晶格条纹,分别对应于NiO的(200)晶面和FeS的(202)晶面[21],证实NiO-FeS@NF具有良好的结晶度,同时也表明FeS与NiO形成了较好的异质结[22]。一方面这种异质结具有特殊的电荷转移特性,另一方面可以保证相邻组分的良好接触,进而实现有效的界面电荷传输,提升电催化性能[23-25]。
2.1.2 NiO-FeS@NF的化学成分
对NiO-FeS@NF复合材料进行XPS分析(图 2),在Ni 2p谱(图 2A)中,结合能为855.6 eV和873.7 eV处的主峰分别代表Ni 2p3/2和Ni 2p1/2,对应于Ni—O键中的Ni2+[26-27]。类似地,Fe 2p谱(图 2B)中结合能为711.5 eV和725.4 eV处的主峰分别对应于Fe 2p1/2和Fe 2p3/2,表明复合材料中存在Fe2+和Fe3+,其中Fe3+可能来自于材料表面的氧化[5, 28]。O 1s谱(图 2C)中结合能为530.6、531.6、533.1 eV处的谱峰分别代表Ni—O、O—H以及吸附的水或氧,从而证实了NiO的存在。S 2p谱(图 2D)中结合能为160.0、167.9 eV处的峰证明存在金属硫化物(S2-),材料表面少部分已被氧化成SOx[16, 28-30]。采用两步电沉积法成功制备出具有二维片状结构的NiO-FeS@NF复合材料。
2.1.3 NiO-FeS@NF的生长机理
根据前述对NiO-FeS@NF表面形貌和化学成分的XPS分析,提出NiO-FeS@NF复合材料在两步电沉积法制备过程的生长机理(图 3)。电沉积法是一种“自下而上”的自组装过程[7]。在第一步沉积过程中,电解液中的Fe2+和S2-在电场作用下发生迁移并在NF基底表面富集,2种离子发生反应生成FeS核。随着沉积时间的延长,FeS核逐渐定向生长,最后形成二维片状结构。在第二步沉积过程中,电解液中的Ni2+也迁移到FeS表面并逐渐富集,最终形成NiO-FeS二维复合材料。
2.2 NiO-FeS@NF电化学性能
2.2.1 不同沉积条件的优化结果
电催化性能与电催化材料的物理、化学和电子特性等相关,往往可以通过制备过程中的参数进行调控。因此,考察电沉积过程所用不同电解液中硫酸镍的浓度、沉积电位和沉淀时间等对NiO-FeS@NF复合材料OER性能的影响。
采用3种电解液硫酸镍浓度(0.05、0.10、0.20 mol/L),沉积15 min,沉积电位-1.0 V vs Ag/AgCl,制备3种NiO-FeS@NF复合材料,并对材料进行线性扫描伏安测试(图 4A)。当c(Ni2+)为0.05、0.10、0.20 mol/L时,OER过电位分别为263.8、270.2、271.0 mV。一般认为,金属硫化物中的金属为催化活性中心,本研究在FeS表面沉积NiO后成功在两者界面处构筑了异质结,可加速电子转移,同时还引入镍铁双金属协同效应,提升电催化OER性能。然而,当c(Ni2+)过低(0.05 mol/L)时,协同效应不明显,当c(Ni2+)过高(0.20 mol/L)时,过量的NiO会影响Fe活性中心,造成OER性能降低。因此,本研究最佳的c(Ni2+)选用0.10 mol/L。
改变操作参数考察沉积电位和沉积时间对OER性能的影响(图 4B)。采用3种沉积电位(-0.9、-1.0、-1.1 V vs Ag/AgCl),沉积15 min,c(Ni2+)为0.10 mol/L。当沉积电位为-1.0 V时,NiO-FeS@NF的OER过电位最低(263.8 mV)。一般外加电位的大小能够改变离子的迁移速度,过大或过小都会影响材料的均匀生长,因此本研究最佳的沉积电位选用-1.0 V vs Ag/AgCl。
考察3种沉积时间(10、15、20 min)对材料OER性能的影响(图 4C),沉积电位为-1.0 V vs Ag/AgCl,c(Ni2+)为0.1 mol/L。当沉积15 min时,材料的OER过电位最低。我们认为沉积时间较短时,NiO不足以形成良好的超薄片状结构,而沉积时间较长时,所得纳米片结构可能太厚、太密集,进而影响活性位点的数量。因此,采用的最佳沉积时间为15 min。
综合以上实验可确定最佳沉积条件:电解液中硫酸镍的浓度为0.1 mol/L,沉积电位为-1.0 V vs Ag/AgCl,沉积时间为15 min。在该条件下,FeS和NiO界面的电子传输性能最好,并且镍铁双金属的协同效应最佳,使得复合材料的OER性能最优。
2.2.2 不同材料的OER性能
对比NiO-FeS@NF复合材料和单独FeS@NF、单独NiO@NF的OER性能(图 5A)结果表明:NiO的加入使FeS@NF的氧化还原峰的峰电流变大,反应进行速率加快; NiO-FeS@NF仅需263.8 mV的过电位即可达到电流密度10 mA/cm2,过电位低于FeS@NF的282.2 mV和NiO@NF的319.2 mV; NiO-FeS@NF的过电位降低说明NiO-FeS@NF复合材料的OER性能得到了提升。
通过塔菲尔曲线(图 5B)进一步考察电催化剂的OER动力学,塔菲尔曲线的线性部分能够反映较简单的电子传递过程。通过计算得到NiO-FeS@NF的塔菲尔斜率仅为23.7 mV/dec,远远低于FeS@NF的55.1 mV/dec和NiO@NF的114.7 mV/dec[8, 27, 31-32]。塔菲尔斜率越小,表明在相同动力学电流密度下该催化过程的过电位越低,进一步说明NiO-FeS@NF复合材料的OER动力学性能得到了大幅提升。
2.2.3 复合材料OER稳定性测试
对NiO-FeS@NF电催化稳定性进行研究,采用i-t法测试材料在1 mol/L KOH溶液中的电流曲线(图 6)。在1.494 V vs RHE的恒电位下,NiO-FeS@NF的电流密度在10 h内一直保持在10 mA/cm2,未发现明显的衰竭现象,表明电极优异的稳定性。首先,NiO-FeS@NF复合材料独特的二维片状结构能够快速释放产生的O2,材料在连续OER反应过程中未出现严重变形;其次,NiO-FeS在NF基底表面原位生长,不需要使用任何粘合剂,因此材料能够与基底紧密连接,不易脱落。
2.3 NiO-FeS@NF的OER机理探讨
电催化剂的催化性能往往与电极表面的电荷转移电阻有关。Nyquist图可以反映表面电荷的转移过程,通过拟合电化学阻抗谱(EIS)可以评价电极的电导率,进而揭示OER反应机理。在相同电位(0.45 V vs Ag/AgCl)下,在10-2~105 Hz频率范围测试了几种材料的EIS谱Nyquist图(图 7)。单独NiO@NF和单独FeS@NF电极具有较大的半圆直径,说明它们的电阻较大,而NiO-FeS@NF复合材料的电阻仅为7 Ω,表明其导电性最好。在复合材料中,FeS和NiO形成了良好的异质结,同时存在镍铁双金属协同效应,使两者界面处的电子转移速率加快,从而OER电催化性能最好[12, 16, 33]。
一般认为,在碱性电解液中OER反应主要分为4个步骤[34-35]:
$$ {\rm{M}} + {\rm{O}}{{\rm{H}}^ - } \to {\rm{M — OH + }}{{\rm{e}}^ - }, $$ (3) $$ {\rm{M — OH + O}}{{\rm{H}}^ - } \to {\rm{M — O + }}{{\rm{H}}_2}{\rm{O}} + {{\rm{e}}^ - }, $$ (4) $$ {\rm{M — O + O}}{{\rm{H}}^ - } \to {\rm{M — OOH + }}{{\rm{e}}^ - }, $$ (5) $$ {\rm{M — OHH + O}}{{\rm{H}}^ - } \to {\rm{M}} + {{\rm{O}}_2} + {{\rm{H}}_2}{\rm{O}} + {{\rm{e}}^ - }。 $$ (6) OER是一种四电子-质子耦合反应,涉及M—OH,M—O和M—OOH中间体的形成。这种间接途径导致电催化反应需要更高的能量和过电位,使OER动力学过程变得极为缓慢,主要取决于第三步(式(5)),因此调控M—O的结合能对OER性能的提升至关重要。制备的NiO-FeS@NF复合材料具有二维超薄结构,能够为电催化反应提供较大的比表面积和更多暴露的活性位点,同时异质结的存在促进了界面电子转移过程,降低了M—O的结合能,进而降低了OER反应中的动力学势能垒,使NiO-FeS@NF复合材料的OER性能得到提升。
3. 结论
通过简易的电沉积法制备了二维复合材料NiO-FeS@NF,将其应用于电解水OER反应。SEM、TEM和XPS等表征结果表明:NiO-FeS@NF复合材料具有二维超薄片状结构,FeS和NiO能够形成良好的异质结。由于电沉积法的原位合成,NiO-FeS材料和NF基底之间存在较好的电子转移,因此可直接用作电催化OER电极。通过调节电沉积参数确定了最佳沉积条件:电解液中硫酸镍的浓度为0.1 mol/L,沉积电位为-1.0 V vs Ag/AgCl,沉积时间为15 min。在1 mol/L KOH碱性电解液中,NiO-FeS@NF电极展现出优异的OER性能,当电流密度为10 mA/cm2时,OER过电位为263.8 mV,塔菲尔斜率为23.7 mV/dec,优于单独的NiO@NF和FeS@NF电极。通过EIS谱分析发现NiO-FeS@NF电阻仅为7 Ω,在所有样品中最小,说明其电子转移速度最快。一方面,二维超薄结构为电催化OER提供更多的活性位点;另一方面,FeS和NiO形成的异质结以及镍铁双金属协同效应能够增强界面电子传递,促进OER动力学过程。这种简单的电沉积法有利于制备高效的催化剂材料,有望在能源和环境等领域得到广泛应用。
-
表 1 巴西劈裂圆盘材料主要物性参数
Table 1 The main material parameters of Brazilian test modelling
材料 温度/℃ 密度/(g·mm-3) 弹性模量/GPa 泊松比 单轴抗压强度/MPa 承压板 — 7.8 210.00 0.077 — 花岗岩 25 2.7 39.30 0.200 90.42 150 2.7 32.26 0.200 87.73 200 2.7 29.22 0.200 86.35 300 2.7 27.43 0.200 79.36 400 2.7 18.97 0.200 83.77 500 2.7 20.99 0.200 68.20 600 2.7 19.87 0.200 69.56 650 2.7 16.16 0.200 59.83 表 2 花岗岩的巴西劈裂圆盘的抗拉强度模拟结果
Table 2 The numerical results of Brazilian tensile strength of granite discs
温度/℃ 弹性模量/GPa 单轴抗压强度/MPa 最大载荷/N 巴西劈裂强度/MPa 25 39.30 90.42 621.77 13.19 150 32.26 87.73 532.13 11.29 200 29.22 86.35 494.11 10.49 300 27.43 79.36 474.93 10.08 400 18.97 83.77 374.06 7.94 500 20.99 68.20 382.35 8.11 600 19.87 69.56 140.66 2.99 650 16.16 59.83 122.71 2.60 -
[1] 林伯强. 中国能源发展报告2018[M]. 北京: 北京大学出版社, 2019. [2] 蔺文静, 刘志明, 马峰, 等. 我国陆区干热岩资源潜力估算[J]. 地球学报, 2012, 33(5): 807-811. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201205018.htm LIN W J, LIU Z M, MA F, et al. An estimation of HDR resources in China's mainland[J]. Acta Geoscientica Sinica, 2012, 33(5): 807-811. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201205018.htm
[3] 王贵玲, 张薇, 梁继运, 等. 中国地热资源潜力评价[J]. 地球学报, 2017, 38(4): 449-450. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201704002.htm WANG G L, ZHANG W, LIANG J Y, et al. Evaluation of geothermal resources potential in China[J]. Acta Geoscientica Sinica, 2017, 38(4): 449-450. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201704002.htm
[4] 董晋鹏. 高温后不同粒径花岗岩三轴力学特性及颗粒流模拟研究[D]. 北京: 中国矿业大学, 2020. [5] 李天斌, 高美奔, 陈国庆, 等. 硬脆性岩石热-力-损伤本构模型及其初步运用[J]. 岩土工程学报, 2017, 39(8): 1477-1484. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708019.htm LI T B, GAO M B, CHEN G Q, et al. A thermal-mechanical-damage constitutive model for hard brittle rocks and its preliminary application[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(8): 1477-1484. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201708019.htm
[6] 杜守继, 刘华, 职洪涛, 等. 高温后花岗岩力学性能的试验研究[J]. 岩石力学与工程学报, 2004, 23(14): 2359-2364. doi: 10.3321/j.issn:1000-6915.2004.14.010 DU S J, LIU H, ZHI H T, et al. Testing study on mechanical properties of post-high-temperature granite[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(14): 2359-2364. doi: 10.3321/j.issn:1000-6915.2004.14.010
[7] ZHANG W, SUN Q, ZHU S, et al. Experimental study on mechanical and porous characteristics of limestone affected by high temperature[J]. Applied Thermal Engineering, 2017, 110: 356-362. doi: 10.1016/j.applthermaleng.2016.08.194
[8] YIN T, LI X, CAO W, et al. Effects of thermal treatment on tensile strength of Laurentian Granite using Brazilian test[J]. Rock Mechanics & Rock Engineering, 2015, 48(6): 2213-2223.
[9] CHEN Y L, WANG S R, NI J, et al. An experimental study of the mechanical properties of granite after high temperature exposure based on mineral characteristics[J]. Engineering Geology, 2017, 220: 234-242. doi: 10.1016/j.enggeo.2017.02.010
[10] KANG F, LI Y, TANG C. Grain size heterogeneity controls strengthening to weakening of granite over high-temperature treatment[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 145: 104848/1-17.
[11] 张宁. 高温三轴应力下花岗岩蠕变-渗透-热破裂规律与地热开采研究[D]. 太原: 太原理工大学, 2013. [12] 程绍桐. 加热方式对花岗岩热损伤影响的研究[D]. 北京: 北京工业大学, 2020. [13] 李春, 胡耀青, 张纯旺, 等. 不同温度循环冷却作用后花岗岩巴西劈裂特征及其物理力学特性演化规律研究[J]. 岩石力学与工程学报, 2020, 39(9): 1797-1807. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202009007.htm LI CHUN, HU Y Q, ZHANG C W, et al. Brazilian split characteristics and mechanical property evolution of granite after cyclic cooling at different temperatures[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(9): 1797-1807. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202009007.htm
[14] 余莉, 彭海旺, 李国伟, 等. 花岗岩高温-水冷循环作用下的试验研究[J]. 岩土力学, 2021, 42(4): 1025-1035. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104015.htm YU L, PENG H W, LI G W, et al. Experimental study on granite under high temperature-water cooling cycle[J]. Rock and Soil Mechanics, 2021, 42(4): 1025-1035. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202104015.htm
[15] LIU S, XU J. An experimental study on the physico-mechanical properties of two post-high-temperature rocks[J]. Engineering Geology, 2015, 185: 63-70.
[16] MIAO S, PAN P Z, ZHAO X, et al. Experimental study on damage and fracture characteristics of Beishan Granite subjected to high-temperature treatment with DIC and AE techniques[J]. Rock Mechanics and Rock Engineering, 2021, 54: 721-743.
[17] QIN Y, TIAN H, XU N X, et al. Physical and mechanical properties of Granite after high-temperature treatment[J]. Rock Mechanics and Rock Engineering, 2020, 53(1): 305-322.
[18] 郤保平, 吴阳春, 王帅, 等. 青海共和盆地花岗岩高温热损伤力学特性试验研究[J]. 岩石力学与工程学报, 2020, 39(1): 69-83. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202001007.htm XI B P, WU Y C, WANG S, et al. Experimental study on mechanical properties of granite taken from Gonghe basin, Qinghai province after high temperature thermal damage[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(1): 69-83. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202001007.htm