1、SHANDONG ELECTRIC POWER山东电力技术第50卷（总第312期）2023年第11期DOI：10.20097/ki.issn1007-9904.2023.11.001硅烷偶联剂改性TiO2掺杂对绝缘纸中水分子扩散行为的影响李健飞1，2，王鑫1，2，王洋1，2，李华强3，朱庆东4（1.西安工程大学电子信息学院，陕西西安710048；2.西安市电气设备互联感知与智能诊断重点实验室，陕西西安710048；3.西安交通大学电力设备电气绝缘国家重点实验室，陕西西安710049；4.国网山东省电力公司电力科学研究院，山东济南250003）摘要：在变压器绝缘纸的改性研究中，纳米粒子和水分含量
2、均是影响绝缘纸绝缘性能的重要因素，然而有关硅烷偶联剂改性TiO2掺杂对绝缘纸中水分子扩散行为的影响尚未可知。采用分子动力学方法，建立了3种硅烷偶联剂改性前后TiO2与绝缘纸的纳米尺度界面模型（未改性TiO2/绝缘纸、KH570改性TiO2/绝缘纸、KH792改性TiO2/绝缘纸），并对不同界面模型中水分子的扩散行为进行动力学计算。结果表明，在未改性TiO2/绝缘纸的界面模型中水分子的扩散系数最大，质心运动范围最广，且界面附近水分子分布最少；而在KH570和KH792改性TiO2/绝缘纸界面模型中水分子的扩散系数及其质心运动范围减小，界面附近水分子分布浓度增大，且KH792改性后TiO2对绝缘纸
3、中水分子扩散的抑制能力最强。此外，通过分析水分子在3种界面模型中自由体积分数、相互作用能和氢键数目，诠释在上述不同模型中水分子扩散系数差异的物理机理。研究思路与结果为纳米粒子掺杂绝缘纸的改性实验提供了理论支撑。关键词：绝缘纸；水分子扩散；TiO2掺杂；硅烷偶联剂；表面改性中图分类号：TM214；TM215；TM411文献标识码：A文章编号：1007-9904（2023）11-0001-10Effect of Silane Coupling Agent Modified TiO2on the DiffusionBehavior of Water Molecules in CelluloseLI
4、Jianfei1，2，WANG Xin1，2，WANG Yang1，2，LI Huaqiang3，ZHU Qingdong4（1.School of Electronics and Information，Xi an Polytechnic University，Xi an 710048，China；2.Xi an Key Laboratory of Interconnected Sensing and Intelligent Diagnosis for Electrical Equipment，Xi an 710048，China；3.State Key Laboratory of Elec
5、trical Insulation and Power Equipment，Xi an Jiaotong University，Xi an 710049，China；4.State Grid Shandong Electric Power Research Institute，Jinan 250003，China）Abstract：In the research on the modification of transformer insulation paper，both nanoparticles and moisture content areimportant factors affe
6、cting the insulation performance of insulation paper.However，the influence of silane coupling agentmodified TiO2on the diffusion behavior of water molecules in insulation paper is unknown.Therefore，by molecular dynamics，weestablished three types of nano-scale interface models of TiO2and insulating p
7、aper before and after the modification，such asunmodified TiO2/insulating paper，KH570 modified TiO2/insulating paper，KH792 modified TiO2/insulating paper.From this，weconducted kinetics calculation on the diffusion behavior of water molecules in different interface models.The results reveal thatthe di
8、ffusion coefficient of water molecules in the unmodified TiO2/insulating paper interfacial model is the largest，accompaniedby the widest range of centroid motion and the smallest distribution of water molecules near the interface.Conversely，in theKH570 and KH792 modified TiO2/insulating paper interf
9、ace models，the diffusion coefficient and range of centroid motion ofwater molecules decrease with the increase of water molecules concentration near the interface.Meanwhile，the KH792 modifiedTiO2demonstrate the strongest ability to inhibit the diffusion of water molecules in the insulating paper.Fur
10、thermore，the study 面向新型电力系统的电气设备绝缘技术 基金项目：陕西省自然科学基础研究计划项目（2021JM450）；中国博士后科学基金（2019M653630）；国网山东省电力公司科技项目（52062623000G）。Natural Science Basic Research Program of Shaanxi Province（2021JM450）；China Postdoctoral Science Foundation（2019M653630）；Science andTechnology Project of State Grid Shandong Electr
11、ic Power Company（52062623000G）.1山东电力技术第50卷（总第312期）2023年第11期elucidates the underlying physical mechanisms behind the variations in diffusion coefficients according to the free volumefraction，interaction energy，and number of hydrogen bonds of water molecules in the three types of interface models.The
12、researchmethod and results provide theoretical support for the experimental modification of nanoparticle-doped insulating paper.Keywords：insulating paper；diffusion of water molecules；TiO2doping；silane coupling agent；surface modification0引言油浸式变压器是电力系统的核心设备，承担着电压变换和电能传输的重要任务，其稳定可靠运行对电网安全至关重要1-2。绝缘纸作为油
13、纸绝缘系统的重要组成部分，在油浸式变压器中起着绝缘和支撑的重要作用3。然而，在变压器长期运行过程中，由于受到温度、水分和电场等环境因素的影响，导致绝缘纸逐渐劣化，严重时会发生局部放电，严重影响变压器的安全运行与使用寿命。因此，对绝缘纸改性以提升其电气特性一直是输变电设备领域的研究热点。近年来，纳米掺杂改性技术被逐步应用到变压器油纸绝缘当中4，并取得了一系列较好的实验研究成果。Choi 等5-6研究发现在变压器油中添加一定浓度的纳米粒子能够改善变压器油的散热特性，并首次提出了纳米流体的概念。研究发现，在变压器油中加入适量纳米粒子，不仅能够提高变压器油的导热性能，还可以提升其介电性能7-8和抗老化
14、特性9。在绝缘纸改性方面，张松等10和 Gao 等11发现将适量的纳米氧化铝添加至绝缘纸中，不仅可以提高其机械性能，还可以改善其电气性能。Liao等12研究发现 TiO2和 AlN 纳米粒子掺杂均可有效提高油浸绝缘纸的交流击穿强度。Zhang 等13对纳米Al2O3掺杂改性绝缘纸的陷阱分布特性进行了测量，发现纳米粒子的引入可以提高绝缘纸的陷阱能级，增大深陷阱密度，减缓其表面电位的衰减速度，进而提高绝缘纸的击穿场强。在掺杂改性实验中，为解决无机纳米材料与绝缘纸之间存在相容性差和分散不均匀的问题，Liu 等14和 Huang 等15在掺杂前采用KH550 和油酸对纳米 TiO2进行表面处理，两种处
15、理方法均可提升纳米粒子在绝缘纸中的分散性与相容性。虽然目前纳米粒子掺杂油纸绝缘的改性实验取得大量研究进展，然而其改性机理尚不明确，尤其是在物理层面与分子层面。众所周知，水分是影响绝缘纸老化的重要因素之一，同时又是绝缘纸裂解的重要产物，当变压器中水分含量增加时，会加速油纸绝缘老化，降低其绝缘性能和电气特性。基于此，部分学者采用分子模拟的方法，研究了纸中含水量对其裂解过程的影响，发现高温条件下水分子易与纤维素二糖分子形成氢键，破坏纤维素二糖分子自身稳定结构，促进纤维素二糖分子裂解，加速绝缘纸的老化16。为研究纳米掺杂对油纸绝缘中水分扩散行为的影响，Tang 等17利用分子动力学的方法，模拟了SiO
16、2掺杂绝缘纸中水分子的扩散过程，发现 SiO2纳米粒子可以有效抑制水分子在绝缘纸中的扩散行为，使水分子与纤维素之间的氢键数目减少，降低水分子对纤维素二糖自身结构的破坏作用，提高了绝缘纸的抗老化性能。虽然上述研究分析了纳米粒子引入对绝缘纸中水分子扩散行为的影响，但在实际掺杂过程中，硅烷偶联剂改性纳米粒子可以有效增大纳米粒子的比表面积，降低纳米粒子的表面能，提高纳米粒子在有机材料中的相容性，因此，硅烷偶联剂改性方法常用于纳米粒子掺杂前的预处理，以提升纳米粒子在基底材料中的分散性18。截至目前，针对不同偶联剂改性的纳米粒子掺杂对纸中水分子扩散行为的影响未见报道。因此，采用分子动力学的方法，研究两种常
17、见硅烷偶联剂改性前后 TiO2掺杂绝缘纸中水分子的扩散行为，通过分析绝缘纸中水分子的扩散系数、质心轨迹、相对分子浓度及其与 TiO2纳米粒子之间的氢键数目和相互作用能的大小，从物理层面和分子层面诠释了不同硅烷偶联剂的改性机制及其对水分子扩散行为的影响。1分子建模和模拟掺杂实验中的 TiO2纳米颗粒的粒径均为 30 nm以上，而目前针对变压器油纸绝缘电气与理化性能分子动力学仿真的分子模型尺度局限在 30 30 30 以下19。因此，在现有计算资源情况下难以构建粒径为 30 nm 的 TiO2纳米颗粒在保持高计算精度2的条件下完成分子动力学计算。与此同时，纳米球粒径越大，其表面曲率则越小，在局部范
18、围内更接近平面结构，且纳米粒子掺杂对基底材料的影响更多取决于表面积的大小，与表面曲率关联性较弱。故采用绝缘纸与 TiO2构成的平面界面来近似表征绝缘纸与 TiO2纳米颗粒之间的界面，在 Material Studio2020 软件中完成建模、优化以及水分子迁移行为的模拟计算。1.1绝缘纸层分子建模变压器绝缘纸的主要成分为纤维素，是由纤维素二糖分子聚合而成的聚合物。实际绝缘纸的聚合度通常在 1 000 以上，单体分子链较长，致使建模体系过大难以完成动力学计算，故选择聚合度为 10 的纤维素分子链进行绝缘纸建模。虽然这一聚合度显著低于实际绝缘纸正常状态下的聚合度，但其单体结构和官能团并未发生变化，
19、因此选择低聚合度的纤维素分子对绝缘纸进行分子建模是目前分子模拟的常用手段20。为减少绝缘纸与 TiO2界面处的真空层区域，采用 Confined Layer 模块构建了在 Z 轴方向具有平面界面、尺寸为 30.6 30.2 20 的绝缘纸层。变压器中的绝缘纸在投入运行初期，其含水量通常小于 0.5%，随着运行时间的加长以及油纸绝缘老化等因素的影响，后期含水量通常高于 3%。因此，在绝缘纸层中加入了质量分数为 3%的水分子。该结构含有 5 个聚合度为 10 的纤维素分子，密度为1.442 g/cm3，与实验值吻合良好21。纤维素二糖分子、水分子和绝缘纸平面层的分子结构分别如图 1（a）、图 1（
20、b）、图 1（c）所示。（a）纤维素二糖分子（b）水分子（c）绝缘纸平面层图1绝缘纸层及其组分分子结构Fig.1 Insulation paper layer and its componentmolecular structure1.2TiO2层分子建模将 TiO2的元胞利用 Cleave Surface 沿晶面指数为（1 0 1）的平面切割，并对其能量进行优化，在此基础上建立超晶包为（38）TiO2（1 0 1）表面结构，将表面未饱和 Ti 原子上进行加氢处理，即完成了未改性TiO2平面结构的构建，尺寸为 30.6 30.2 25。硅烷偶联剂对纳米粒子改性的本质是在其表面接枝含有不同官能团
21、结构的特征分子链的过程。选用的偶联剂为实验中改性效果较好的两种硅烷偶联剂，KH570 和 KH792，其改性后接枝在纳米粒子表面的特征官能团分子结构如图 2 所示，其能量优化后分子链的长度分别为 11.767 与 12.823。在 TiO2表面的 Ti 原子上接枝相应的官能团即可构建经过KH570 和 KH792 改性的 TiO2平面结构。（a）KH570（b）KH792图2KH570和KH792接枝官能团分子结构Fig.2 Molecular structure of KH570 and KH792 graftedfunctional groups1.3绝缘纸/TiO2界面结构分子建模在 B
22、uild Layers 模块中将绝缘纸层和改性前后李健飞，等：硅烷偶联剂改性TiO2掺杂对绝缘纸中水分子扩散行为的影响3山东电力技术第50卷（总第312期）2023年第11期TiO2层合并，构建了尺寸为 30.6 30.2 45 的绝缘纸/TiO2界面模型，如图 3 所示。其中上层为绝缘纸层，下层为 TiO2层。初始界面模型建立后，对该模型先后进行几何优化、100 ps 的正则系综（canonicalensemble，NVT）弛 豫 和 100 ps 的 恒 温 恒 压 系 综（constant-pressure，constant-temperature，NPT）弛豫，压强设置为 101.32
23、5 kPa22，其中几何优化采用Smart 优化算法，以 4.186 810-4kJ/mol 的能量收敛值为目标进行最大迭代次数为 10 000 次的结构寻优。弛豫过程中控制器及相关参数的设置如表 1所示。（a）未改性TiO2/绝缘纸界面（b）KH570改性TiO2/绝缘纸界面（c）KH792改性TiO2/绝缘纸界面图3绝缘纸与TiO2界面结构分子模型Fig.3 Molecular model of the interface structurebetween insulating paper and TiO2表1控制器及相关参数设置Table 1 Controller and related
24、 parameter settings条件参数温度控制器压力控制器模拟步长/fs模拟温度/K模拟压力/kPaNVTNos-Hoover1343NPTNos-HooverBerendsen1343101.325在完成几何优化和弛豫后，绝缘纸/TiO2界面模型进入平衡态，随后选取 COMPASS 力场、Nos-Hoover 温度控制器对优化后的界面模型在 NVT（温度为 343 K）系综内进行了时长为 500 ps、步长为 1 fs的动力学仿真，并以 1 帧/ps 的频率进行数据采集，用于研究不同硅烷偶联剂改性前后 TiO2表面对绝缘纸中水分子扩散行为的影响。2模拟结果与分析2.1水分子扩散系数水
25、分子在绝缘纸中的扩散行为和运动状态可以用其均方位移（mean square displacement，MSD）曲线来描述，MSD 曲线代表粒子随时间移动后的位置相对于参考位置偏差的量度，用SMD表示，如式（1）所示。SMD=|r(t)-r(0)2（1）式中：r(t)、r(0)分别为水分子在 t 时刻和 0 时刻的位置向量；为计算结果的统计平均。3 种 TiO2/绝缘纸界面模型中水分子的 MSD 曲线如图 4 所示。图43种界面模型中水分子的MSD曲线Fig.4 MSD curves of water molecules in threeinterface models从图 4 中可以看出，未改
26、性的 TiO2/绝缘纸界面模型中水分子的 MSD 曲线斜率最大，而经过 KH570和 KH792 改性后的 TiO2/绝缘纸界面模型中水分子的 MSD 曲线斜率逐渐减少，且经过 KH792 改性后的界面模型中水分子的 MSD 曲线斜率最小。4扩散系数 D 是表征物质扩散能力强弱的重要参数，扩散系数越大表明粒子受到的束缚就越小，扩散能力越强，其计算方法如式（2）所示。D=16Nlimt ddti=1N|ri(t)-ri(0)2（2）式中：N 为模型中水分子的数目；ri(t)和ri(0)分别为第 i 个原子在 t 时刻与 0 时刻的位置向量。式中的微分部分可通过 MSD 曲线拟合直线的斜率近似代替
27、。由于均方位移的计算值为 N 个水分子随时间移动后的位置相对于参考位置偏差量的统计平均，因此，式（2）可简化为D=16a（3）式中：a 为 MSD 曲线拟合直线的斜率。为进一步分析偶联剂改性 TiO2对绝缘纸中水分子扩散行为的影响，对不同界面模型中水分子的扩散系数进行了统计对比，其结果如表 2 所示。表23种界面模型中水分子MSD曲线斜率及扩散系数Table 2 Slope of the MSD curve and diffusion coefficient ofwater molecules in three interface models单位：2/ps参数aD未改性3.210.54KH5
28、70 改性1.500.25KH792 改性0.890.15由表 2 可知，未改性的 TiO2/绝缘纸界面模型中水分子的扩散系数最大；经过 KH570 改性的 TiO2/绝缘纸界面模型中水分子的扩散系数次之；而经过KH792 改性的 TiO2/绝缘纸界面模型中水分子的扩散系数最小，较未改性的 TiO2/绝缘纸界面模型中水分子的扩散系数降低了 72%。综上，经过硅烷偶联剂改性的 TiO2能够进一步抑制绝缘纸中水分子的扩散能力，且 KH792 的改性效果更优。这是因为 KH792 中的极性官能团氨基（-NH2）和亚氨基（-NH-）的极性均强于 KH570 中酯基（-COOR）的极性，因此，KH792
29、对绝缘纸中水分子的吸附效果更强，从而抑制了水分子的扩散能力，降低了水分子的扩散系数。2.2水分子质心轨迹为更直观地观察水分子在改性前后 TiO2/绝缘纸界面模型的运动情况，对 NVT 系综内 500 ps 动力学模拟中的水分子的质心轨迹进行了统计平均。3种界面模型中水分子的质心轨迹如图 5 所示，图中红色球棍代表水分子的质心运动轨迹，蓝色圆点代表水分子质心轨迹在 YZ 平面上的投影。为定量描述水分子的运动范围，对 3 种模型中水分子在 X 轴、Y轴、Z 轴 3 个方向的运动范围进行了统计计算，结果如表 3 所示。表3水分子运动范围Table 3 Range of motion of water
30、 molecules单位：类型未改性KH570 改性KH792 改性X 轴10.735.254.24Y 轴17.56.343.06Z 轴5.086.966.04从表 3 可以看出，相比于未改性的界面模型，经KH570 和 KH792 改性后的界面模型中水分子在 X轴、Y 轴、Z 轴方向的运动范围大幅减少，表明改性后的 TiO2表面对水分子有着更好的束缚效果，且KH792 的改性效果最优，与上文水分子扩散系数变化规律相吻合。2.3水分子相对分子浓度水分子相对分子浓度是描述水分子在整个界面模型中分布情况的重要参数，其物理意义为模型中沿某方向单位距离内水分子的个数。通过对比模拟前后 3 种界面模型中
31、水分子的分布情况，可以清晰地反映水分子在不同模型中的空间位置变化及分子运动的变化趋势，进而分析不同偶联剂改性 TiO2表面对水分子扩散行为的影响。基于此，计算了绝缘纸/TiO2界面法向（沿 Z 轴方向）的水分子相对分子浓度分布情况，其计算结果如图 6 所示。（a）绝缘纸/TiO2（未改性）界面李健飞，等：硅烷偶联剂改性TiO2掺杂对绝缘纸中水分子扩散行为的影响5山东电力技术第50卷（总第312期）2023年第11期（b）绝缘纸/TiO2（KH570改性）界面（c）绝缘纸/TiO2（KH792改性）界面图5水分子质心轨迹Fig.5 Centroid trajectory of water mol
32、ecules图6水分子相对分子浓度Fig.6 Relative molecular concentration of water molecules图 6 中 蓝色矩形表示模型中界面所在位置。从图 6 可知，水分子在模拟结束时均分布于 TiO2/绝缘纸界面附近，说明改性前后的 TiO2表面对水分子均存在较强的吸附作用，因此在绝缘纸中添加 TiO2纳米粒子可以促使水分子吸附于纳米粒子表面，进而限制水分子的扩散行为。通过分析界面附近水分子浓度的峰值大小可知，未改性的 TiO2表面对水分子的吸附效果较弱，虽然水分子有向 TiO2表面运动的趋势，但绝缘纸中仍存在浓度较高的水分含量，且界面附近水分子相对
33、分子浓度最大值仅有 2.22 个/。相比于未改性的 TiO2/绝缘纸界面模型，经过 KH570 和 KH792 改性后的 TiO2表面对水分子的吸附作用更强，在界面处水分子的相对分子浓度最大值分别为 3.24 个/和2.84个/，水分子向界面扩散的趋势更明显。从图 6 还可知，KH792 改性的界面模型中模拟前后水分子浓度峰值的变化幅度小于 KH570 改性的界面模型中模拟前后水分子浓度峰值的变化幅度，但其峰值附近水分子浓度却明显高于模拟前峰值附近水分子浓度，使得经 KH792 改性后的模型界面附近水分子浓度更高。这是因为经 KH792 改性后在 TiO2表面接枝的特征分子比经 KH570 改
34、性后在 TiO2表面接枝的特征分子的链长更长，拓展范围更广，如图 2 所示。故经 KH792 改性后在 TiO2表面接枝的特征分子在纤维素中分布更远更广，使其在距离界面较远的地方开始吸附水分子，进而使得水分子在 1530 的区域内整体水分含量最高，而在峰值位置变化并不明显，故在 KH792 改性的模型中，水分子含量分布更为广泛。为进一步区分 KH792 和 KH570 对 TiO2表面的改性效果，对 1535 内水分子个数进行统计计算，发现经 KH570 改性后的 TiO2界面附近的水分子数目为 20，而经 KH792 改性后界面附近的水分子数目为 23，多于 KH570 改性后界面附近的水分
35、子浓度。这是因为 KH792 特征分子的链长较长且极性较强，由于其对水分子的吸附作用，水分子会随着其特征分子链在绝缘纸中分布，因此水分子在界面附近的分布范围更广，进而使得界面附近水分子浓度更高。为研究上述现象的物理机制，对不同界面模型中的氢键数目、自由体积分数及水分子与 TiO2表面之间的相互作用能进行了计算。2.4氢键数目氢键是不同于静电力和范德华力的一种特殊的分子间作用力，其键能介于二者之间。通常采用几何判定法判断氢键的生成，即当某一氢原子与其他6电负性基团之间的距离小于 0.25 nm，且电负性基团与氢原子及氢原子供体之间形成的夹角大于 100，则判定二者之间形成了氢键。在 3 种界面模
36、型中，理论上可以形成的氢键种类有两种，第一种是纳米界面，水分子和纤维素分子两两之间形成的氢键；第二种是水分子单体之间或纤维分子单体之间形成的氢键。3种界面模型中氢键数目的统计结果如图7所示。图73种界面模型中的氢键数目Fig.7 Number of hydrogen bonds in three interface models从图 7 中可以看出，经 KH792 改性的 TiO2/绝缘纸界面模型中氢键数目最多，经 KH570 改性的 TiO2/绝缘纸界面模型中氢键数目次之，未经改性的 TiO2/绝缘纸界面模型中氢键数目最少。体系内氢键越多，代表体系内组分分子的约束越强，越不易发生扩散。为进一
37、步研究氢键对体系内水分子扩散行为的影响，对界面模型中不同类型的氢键数目进行了分类统计，结果如表 4 所示。表4各分子间氢键数目Table 4 Number of hydrogen bonds between molecules类型未改性KH570 改性KH792 改性H2O 之间（）112纤维素之间（）99113124H2O/纤维素（）151312H2O/TiO2表面（）589纤维素/TiO2表面（）81011从表 4 可以看出，氢键类型、为水分子形成的氢键，经 KH792 改性后的界面模型中与水分子相关的氢键数目最多；KH570 改性后的界面模型中与水分子相关的氢键数目次之；未改性界面模型中
38、与水分子相关的氢键数目最少。此外，相比于未改性模型，改性界面模型中水分子与 TiO2表面之间氢键数目的增多进一步说明硅烷偶联剂修饰的 TiO2表面可以增强对水分子的吸附作用，且 KH792 的改性效果优于 KH570 的改性效果。上述各模型中氢键数目的变化规律与经不同硅烷偶联剂改性后接枝在 TiO2表面分子中特征官能团的种类与数量密切相关。从图 2 可以看出经 KH792改性后 TiO2表面接枝的特征分子中含有 2 个极性基团（-NH2、-NH-），而经 KH570 改性后 TiO2表面接枝的特征分子中仅含有 1 个极性基团（-COOR）。在TiO2表面接枝相同链数特征分子的情况下，KH792
39、改性的 TiO2表面的极性基团数目是 KH570 改性的TiO2表面极性基团数目的 2 倍。因此，KH792 改性的 TiO2/绝缘纸界面模型中水分子和 TiO2表面之间的氢键数目多于 KH570 改性的 TiO2/绝缘纸界面模型中水分子和 TiO2表面之间的氢键数目。这就从氢键的角度解释了 3 种界面模型中水分子扩散系数产生差异的原因。2.5自由体积分数自由体积分数是描述分子扩散能力强弱的另一重要参量。根据 Fox 和 Flory 的自由体积理论23，材料内部的总体积可分为占有体积 V0和自由体积 Vf。自由体积分数VFF由自由体积和占有体积与自由体积之和的比值确定，即VFF=VfVf+V0
40、 100%（4）图 8 为 3 种界面模型中水分子自由体积分数柱状图，可以看出，未改性界面模型中的水分子的自由体积分数最大，经 KH570 改性的 TiO2/绝缘纸界面模型中水分子的自由体积分数次之，经 KH792 改性的TiO2/绝缘纸界面模型中水分子的自由体积分数最小。这是因为在 TiO2表面接枝的硅烷偶联剂特征分子占据了一定的自由体积，使得界面模型中可用于水分子扩散的自由体积减小，自由体积分数也随之减小。同时，由于硅烷偶联剂中极性基团的存在，使得界面模型中的纤维素分子、水分子与 TiO2表面之间的结合更加紧密，进一步减少了可用于水分子扩李健飞，等：硅烷偶联剂改性TiO2掺杂对绝缘纸中水分
41、子扩散行为的影响7山东电力技术第50卷（总第312期）2023年第11期散的自由体积。由图 2 可知，相比于 KH570，KH792接枝在 TiO2表面特征分子的链长更长，极性官能团更多，因此经 KH792 改性后的模型中分子间的吸引作用更强，结构更为紧凑，进而降低了水分子的自由体积，抑制了水分子扩散。图83种界面模型中水分子自由体积分数Fig.8 Free volume fraction of water molecules in threeinterface models2.6相互作用能相互作用能是评价分子间相互作用强弱的重要物理量，水分子与纳米粒子之间的相互作用能的强弱是影响水分子扩散系
42、数大小的重要因素。相互作用能公式为Eint=ETotal-(EW+EN)（5）式中：ETotal为界面模型的总能量；EW为水分子的势能；EN为纳米粒子表面的势能；Eint为水分子与纳米粒子表面之间的相互作用能。表 5 详细描述了 3 种界面模型中水分子与 TiO2表面之间相互作用能的大小。其中，经 KH792 改性后的 TiO2表面与水分子间的相互作用能绝对值最大，经 KH570 改性后的 TiO2表面与水分子间的相互作用能绝对值次之，未经改性的 TiO2表面与水分子间的相互作用能绝对值最小。这一规律进一步验证了相较于未改性 TiO2，硅烷偶联剂改性 TiO2与水分子之间有着更强的吸附作用，更
43、容易将绝缘纸中水分子吸附到 TiO2表面并约束其分子热运动，从而降低了水分子在绝缘纸中的扩散系数，且相较于KH570，KH792 改性效果更佳。从分子结构层面分析，这是因为在 TiO2表面接枝的硅烷偶联剂含有大量极性基团，这些极性基团的存在，增加了模型中的氢键数目，增强了二者之间的相互作用能，使 TiO2表面对水分子的吸附作用更加明显，从而抑制水分子在界面模型中的扩散能力，降低了水分子的扩散系数。由于经 KH792 改性后接枝到 TiO2表面极性官能团（-NH2、-NH-）的数目多于 KH570 引入极性官能团（-COOR）的数目，且极性更强，所以 KH792 改性效果更佳。这从相互作用能的角
44、度解释了 2.1 节中水分子扩散系数降低的原因。表5水分子与TiO2表面之间的相互作用能Table 5 Interaction energy between water molecules andTiO2surface单位：kJ/mol参数EWENEintTiO2（未改性）31 292.14186 036.27-32 724.03KH570 改性31 384.25185 835.30-32 920.81KH792 改性31 890.86186 672.66-33 356.243结论从微观的角度揭示偶联剂对水分子运动的束缚作用，为纳米粒子掺杂改性绝缘纸的实验研究提供了理论支撑。采用分子动力学方法
45、研究不同偶联剂改性 TiO2纳米粒子对绝缘纸中水分子扩散行为的影响，揭示了 3 种界面模型中水分子扩散系数降低的物理机制，具体结论如下：1）经硅烷偶联剂改性后的 TiO2/绝缘纸界面模型中水分子的扩散系数显著降低，其中，经 KH792改性后模型中的水分子扩散系数最小。2）研究水分子质心轨迹和相对分子浓度，发现相较于未改性的 TiO2/绝缘纸界面模型，经 KH570和 KH792 改性后的 TiO2/绝缘纸界面模型中水分子的质心轨迹依次减小，界面附近的水分子的浓度逐渐增大，这与水分子扩散系数的计算结果保持一致。上述结论均可通过氢键数目、自由体积分数和相互作用能的计算结果进行物理诠释。由于在TiO
46、2表面接枝的偶联剂含有极性较强的特征官能团，这些特征官能团的存在使 TiO2表面与水分子之间的相互作用能增强，氢键数目增多，进而增强了对水分子的束缚作用，降低了水分子的扩散能力；同时，接枝在 TiO2表面的硅烷偶联剂占据一定的自8由体积，使得水分子扩散的空间减少，进一步降低了水分子的扩散系数。KH792 比 KH570 特征分子的链长较长，含有更多的极性官能团，因此其对水分子的束缚作用最优，对水分子扩散行为的抑制作用更为明显。参考文献1顾朝亮，朱孟兆，朱文兵，等.油纸绝缘介电响应宽频等效建模研究 J.山东电力技术，2022，49（5）：43-49.GU Chaoliang，ZHU Mengzh
47、ao，ZHU Wenbing，et al.Broadbandmodeling of the dielectric response of oil paper insulationJ.Shandong Electric Power，2022，49（5）：43-49.2吴东乐，孙胜然，刘群华，等.变压器用绝缘纸击穿强度的多因素协同影响研究 J.绝缘材料，2020，53（10）：96-102.WU Dongle，SUN Shengran，LIU Qunhua，et al.Study on multi-factor synergy effect of breakdown strength for in
48、sulating paper J.Insulating Materials，2020，53（10）：96-102.3梁兆杰，陈硕，田杰.变压器绝缘纸材料的老化力学特性研究J.绝缘材料，2021，54（8）：52-59.LIANG Zhaojie，CHEN Shuo，TIAN Jie.Research on ageingmechanical characteristics of insulating paper for transformerJ.Insulating Materials，2021，54（8）：52-59.4FOFANA I，WASSERBERG V，BORSI H，et al.C
49、hallenge ofmixed insulating liquids for use in high-voltage transformers.1.Investigation of mixed liquidsJ.IEEE Electrical InsulationMagazine，2002，18（3）：18-31.5CHOI S U S，SINGER D A，WANG H P.Developments andapplications of non-newtonian flowsJ.ASME FED，1995，66：99-105.6DA S S K，CHOI S U S，PATEL H E.H
50、eat transfer in nanofluidsa review J.Heat Transfer Engineering，2006，27（10）：3-19.7董明，杨凯歌，马馨逸，等.纳米改性变压器油中载流子输运特性分析 J.电工技术学报，2020，35（21）：4597-4608.DONG Ming，YANG Kaige，MA Xinyi，et al.Analysis of charge-carrier transport characteristics of transformer oil-based nanofluidsJ.Transactions of China Electrot
