稀土金屬改性對負載型Ni2P/M41催化劑 加氫脫硫活性的影響

姜楠 蔣博龍 王帥 宋華

摘要 ：開發高效加氫脫硫催化劑是實現油品深度脫硫的重要途徑。采用浸漬-程序升溫法制備不同稀土金屬（La、Ce、Pr、Nd、Y）改性的Ni 2P/M41催化劑。通過X射線衍射、N 2吸附-脫附、CO吸附和X射線光電子能譜技術分析催化劑的結構以及化學組成，并以二苯并噻吩為模型化合物，考察稀土金屬對Ni 2P/M41催化劑加氫脫硫性能的影響。結果表明：稀土金屬可以不同程度地影響催化劑的尺寸以及活性相Ni 2P在載體表面的分散性；稀土金屬La、Pr能夠顯著提高負載型Ni 2P/M41催化劑的HDS（加氫脫硫）活性，加入稀土金屬Nd、Ce和Y反而使催化劑活性略有降低；Pr-Ni 2P/M41表現出最優的加氫脫硫能力，反應6 h后，脫硫率達到92.6%，與未改性的Ni 2P/M41催化劑的轉化率相比，提高15.1%。

關鍵詞 ：加氫脫硫； 改性； 稀土金屬； 負載型催化劑

中圖分類號 ：O 643 ???文獻標志碼 ：A

引用格式 ：姜楠，蔣博龍，王帥，等.稀土金屬改性對負載型Ni 2P/M41催化劑加氫脫硫活性的影響［J］.中國石油大學學報（自然科學版），2024，48（1）：206-212.

JIANG Nan， JIANG Bolong， WANG Shuai， et al. Influence of rare earth metal modification on hydrodesulfurization activity of supported Ni 2P/M41 catalysts［J］.Journal of China University of Petroleum（Edition of Natural Science），2024，48（1）：206-212.

Influence of rare earth metal modification on hydrodesulfurization

activity of supported Ni 2P/M41 catalysts

JIANG Nan 1， JIANG Bolong 1， WANG Shuai 2， SONG Hua 2

（1.Innovation Institute for Sustainable Maritime Architecture Research and Technology in Qingdao University of

Technology， Qingdao 266000， China;

2.College of Chemistry and Chemical Engineering Department in Northeast Petroleum University， Daqing 163000， China）

Abstract ： It is an important way to achieve deep desulfurization of oil products for the development of efficient hydrodesulfurization （HDS） catalysts. A series of rare earth metals （Nd， La， Ce， Pr， Y） modified nickel phosphide catalysts （ X -Ni 2P/M41） were prepared by a temperature programmed reduction method. The structure and chemical composition of catalysts were analysed by X-ray diffraction （XRD）， N 2 adsorption-desorption， CO uptake and X-ray photoelectron spectroscopy （XPS）. The effect of rare earth metals on the hydrodesulfurization performance of Ni 2P / M41 catalyst was investigated by using dibenzothiophene as a model compound. The results show that the rare earth metals can affect the distribution of Ni 2P crystalline on the surface of M41 and promote the formation of the active Ni 2P crystalline phase at different levels. La and Pr can significantly increase the HDS activity of supported Ni 2P/M41 catalysts， while the addition of Nd， Ce， and Y could slightly reduces its activity. Especially， Pr-Ni 2P/M41 presents the optimal HDS preformance. After 6 h of reaction， its desulfurization conversion reaches 92.6%， which is 15.1% higher than that of the unmodified Ni 2P/M41 catalyst.

Keywords ： hydrodesulfurization; modification; rare earth metals; supported catalyst

工業上主要采用加氫脫硫（HDS）技術來改善油品質量和生產清潔燃料油 ?［1-3］ ，傳統過渡金屬硫化物HDS催化劑因其抗硫穩定性差，以及本身的層狀結構決定了活性中心的分布只能是在其邊緣處等缺點，已難以實現深度及超深度脫硫 ?［4-6］ 。過渡金屬磷化物是一類具有高加氫脫硫活性的催化加氫精制新材料 ?［7-8］ ，其中以Ni 2P催化劑的HDS性能最優，Ni 2P常規制備方法包括程序升溫還原磷酸鹽法（TPR） ?［9］ 、次/亞磷酸鹽歧化法 ?［10］ 和有機磷液相合成法 ?［11］ 等。程序升溫還原磷酸鹽法制備催化劑制備的產物大多分散性差。組成和結構決定著Ni 2P催化劑的催化性能，摻雜金屬元素可以明顯提高其催化活性，其中稀土元素具有獨特的電子結構，最外層有未充滿的4f電子殼層，因而受到廣泛關注 ?［12-15］ 。Wang等 ?［16］ 考察稀土金屬鑭（La）對以TiO 2-Al 2O 3為載體，Ni、W為活性金屬組分的加氫脫硫催化劑的性能影響，發現La可以改善活性組分與載體間的相互作用，且增加Ni-W-S活性相的數量，進而提高催化劑噻吩脫除效率。使用稀土金屬元素Y對非負載的Ni 2P催化劑進行改性，可以抑制Ni ?5 P ?4 的生成，促使生成更多的活性相Ni 2P，從而促進催化劑活性的提高。筆者在傳統程序升溫法的基礎上，制備不同稀土金屬改性的負載型磷化鎳催化劑，以二苯并噻吩（DBT）為模型化合物，結合各種表征方法，考察稀土金屬對負載型Ni 2P/M41催化劑結構及HDS活性的影響。

1 試 驗

1.1 磷化鎳催化劑前驅體的制備

將硝酸鎳和磷酸氫二銨（初始物質的量比 n （Ni）/ n （P）為1/2）溶解于一定量的蒸餾水中，待其混合均勻后，加入載體MCM-41，經過攪拌，干燥以及焙燒等步驟，最終得到初始 n （Ni）/ n （P）為1/1的催化劑前驅體，記為PNi 2P/M41（其中P代表前驅體）。

將催化劑前軀體PNi 2P/M41進行程序升溫還原，升溫速率為2.0 ℃/min，將溫度加熱至600 ℃，保持2 h，然后降至室溫。最終得到Ni負載量（質量分數）為10%，初始 n （Ni）/ n （P）為1/2的負載型磷化鎳催化劑，將其命名為為Ni 2P/M41。

1.2 稀土改性的負載型磷化鎳催化劑的制備

稱取一定量的硝酸鑭、硝酸鈰、硝酸鐠、硝酸釹、硝酸釔加入蒸餾水制成溶液，備用。將稱取的PNi 2P/M41分別加入到已配好的不同稀土金屬溶液中，攪拌、浸漬以及干燥，得到含稀土金屬元素的氧化物前驅體，記為P X -Ni 2P/M41（ X =La、Nd、Ce、Y、Pr）。進行程序升溫還原，得到稀土金屬改性的負載型磷化鎳催化劑，記為 X -Ni 2P/M41。

1.3 催化劑的表征

催化劑的X射線衍射（XRD）。利用日本理學公司D/max-2200PC型X射線衍射儀進行分析，設置管電壓為40 kV，管電流為30 mA，掃描速率為10°/min，掃描范圍10°～80°。催化劑的氮氣吸脫附曲線和比表面積采用NOVA2000e測定儀測定，利用低溫（-196 ℃）氮氣吸附法測定。樣品CO吸附表征采用美國Quantachrome公司化學吸附儀 Autosorb- 1-C儀器測定。催化劑的表面原子價態分析采用Thermofisher Scientific 公司的 K-Alpha 型 X 射線光電子能譜儀（XPS），采用外來污染碳的C 1s作為基準峰來進行校準，以測量值和參考值（284.8 eV）之差作為荷電校正值（Δ）來矯正譜中其他元素的結合能。

1.4 催化劑的活性評價

催化劑活性評價采用自制連續流動式固定床高壓微型反應裝置，以1%（質量分數，下同）二苯并噻吩、98%十氫萘、1%十四烷混合溶液為反應物原料。裝填總量為含0.7 g催化劑的前驅體。將催化劑前驅體用氫氣（200 mL/min）還原，以2.0 ℃/min的升溫速率從室溫升至500 ℃，保持2 h，然后降至設定的反應溫度開始反應。在質量空速（WHSV） 1.5 h ?-1 、氫油體積比為500、反應壓力為3.0 MPa和在340 ℃溫度條件下，液態反應物經過催化劑床層進行反應。每小時取液體樣品1次，樣品用GC-14C型氣相色譜儀（日本島津公司）進行分析。

2 結果分析

2.1 XRD結果

未改性Ni 2P/M41和不同稀土改性的 X -Ni 2P/M41催化劑的XRD譜圖見圖1。由圖1可知，在XRD譜圖中，所有催化劑在2 θ =23.0°處出現了M41的彌散衍射峰。未改性Ni 2P/M41和不同稀土改性的 X -Ni 2P/M41均在2 θ =40.5°、44.6°、47.2°、54.0°、54.7°、66.0°、72.6°和74.4°處出現了較強的Ni 2P物相的衍射峰，且均沒有檢測到磷化鎳的其他衍射峰和稀土金屬的衍射峰。表明催化劑的活性相以單一的Ni 2P形式存在，且稀土金屬以非晶態形式存在或其濃度過低，超出XRD能夠監測到的范圍。與未改性的Ni 2P/M41相比較可以看出，稀土金屬改性的 X -Ni 2P/M41催化劑的Ni 2P相衍射峰較寬，說明稀土金屬的引入促進更小尺寸的Ni 2P相的生成。通過Scherrer公式 ?［17］ 計算Ni 2P {111}面可得到的不同稀土金屬改性的 X -Ni 2P/M41催化劑的Ni 2P物相晶粒見表1。由表1可知，未改性的Ni 2P/M41催化劑的Ni 2P晶粒粒徑為27.3 nm，而Y-Ni 2P/M41和Pr-Ni 2P/M41催化劑的Ni 2P物相晶粒粒徑分別為14.8和15.7 nm。稀土金屬改性的 ?X - Ni 2P/M41催化劑的Ni 2P晶粒粒徑均低于Ni 2P/M41催化劑，表明稀土金屬的引入能夠改變Ni 2P晶相的直徑。

2.2 N 2吸附-脫附

未改性Ni 2P/M41和不同稀土改性的 X -Ni 2P/M41催化劑BET分析結果見表1。由表1可知，未改性的Ni 2P/M41催化劑的比表面積和孔容分別為 145.4 m 2/g和0.118 cm 3/g，平均孔徑為3.3 nm。不同稀土改性的 X -Ni 2P/M41催化劑的孔徑對比Ni 2P/M41，均有所降低，但是比表面積和孔容都有所增大。這說明稀土金屬已進入Ni 2P/M41催化劑的晶格內部，并且稀土金屬替代了部分Ni的位置。 X -Ni 2P/M41催化劑中，La和Pr改性的催化劑的比表面積和孔容最大，說明稀土金屬的引入能提高催化劑的比表面積。

2.3 CO吸附

未改性的Ni 2P/M41和改性的 X -Ni 2P/M41催化劑的CO的吸附容量，得到的結果見表2。未改性的Ni 2P/M41催化劑的CO吸附容量為15.1 μmol/g。引入稀土金屬后， X -Ni 2P/M41催化劑的CO化學吸附量明顯提高，均大于未改性的Ni 2P/M41催化劑，表明稀土金屬改性的催化劑表面暴露的Ni 2P活性位數目更多。稀土改性的催化劑表面金屬分散度（表2）高于未改性催化劑也證明了一點。這是由于引入稀土金屬后，促進了較小晶粒的Ni 2P相形成，有利于促進Ni 2P相的分散性。此外，由于稀土金屬可以抑制催化劑表面的磷富集現象，有利于暴露出更多的Ni 2P活性位。不同稀土金屬改性的催化劑中，Pr-Ni 2P/M41催化劑的CO吸附量最大，為46.5 μmol/g，是未改性Ni 2P/M41催化劑的3.1倍。

2.4 XPS分析

未改性的Ni 2P/M41和不同稀土金屬改性的 ?X - Ni 2P/M41催化劑的Ni 2p和P 2p的XPS譜圖見圖2，相應的XPS分析見表3。由圖2看出，在XPS譜圖中未改性的Ni 2P/M41和不同稀土金屬改性的 ?X - Ni 2P/M41催化劑中，結合能為852.3～ 852.7 ?eV（Ni 2p ?3/2 ）和856.0～856.7 eV（Ni 2p ?3/2 ）的峰分別歸屬于磷化鎳物相中的Ni ??δ + 物種和Ni 2P與P相互作用的Ni ?2+ 物種 ?［18-19］ 。值得注意的是，與未改性的Ni 2P/M41相比，不同稀土金屬改性的 X -Ni 2P/M41催化劑在結合能為852.3～852.7 eV（Ni 2p ?3/2 ）處的峰強度明顯提高，表明加入稀土金屬能有效地促進Ni 2P物相的生成。

對于P 2p，結合能為128.7～130.0 eV的峰歸屬于磷化鎳中的P ??δ - 物種，而結合能為134.2～135.2 eV的峰歸屬于氧化的磷化物。由圖2（b）可知，稀土改性的 X -Ni 2P/M41催化劑對應與P 2p的結合能為134.2～135.2 eV的峰明顯藍移，表明電子由P向稀土金屬轉移。對于所制備的未改性Ni 2P/M14和不同稀土金屬改性的 X -Ni 2P/M41催化劑，可以觀察到結合能為133.4 eV的歸屬于H 2PO - 3的特征峰和結合能為134.5 eV 的歸屬于PO ?3? 4的特征峰， 以及結合能為129.3 eV較弱的特征峰，該峰歸屬于Ni 2P中的P ??δ - 物種 ?［20-21］ 。值得注意的是歸屬于Ni 2P中的P ??δ - 物種的峰的結合能發生了藍移，電子由金屬向P轉移，稀土金屬改變了催化劑表面的電子分布。

利用XPS分析結果計算未改性的Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑中，磷化鎳催化劑表面的Ni ??δ + 物種占Ni 2p ?3/2 物種的比值（ R ??Ni ???δ + ?），其結果列于表3。未改性的Ni 2P催化劑的 R ??Ni ??δ + ?為0.332。與未改性的Ni 2P催化劑的相比，不同稀土金屬改性的 X -Ni 2P/M41催化劑的 R ????δ + 變化不大。

由XPS分析數據得到的未改性Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑表面P/Ni物質量比見表3。由表3可知，雖然各催化劑的初始 n （Ni）/ n （P）為2.0，Ni 2P物相的理論 n （Ni）/ ?n （P） 為0.5，但各催化劑表面 n （Ni）/ n （P）遠高于初始 n （Ni）/ n （P）和理論 n （Ni）/ n （P），說明各催化劑表面均有大量的磷的富集。未改性Ni 2P/M41催化劑 n （Ni）/ n （P）為5.19，與之相比，不同稀土金屬改性的 X -Ni 2P/M41催化劑表面 n （Ni）/ n （P）均提高，表明不同稀土金屬改性的 X -Ni 2P/M41催化劑中，磷富集現象比較嚴重，這可能是稀土改性 X -Ni 2P/M41活性提高不大的原因之一。

2.5 催化劑 HDS活性

2.5.1 稀土金屬對催化劑 HDS 活性的影響

未改性Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑在340 ℃， 3.0 MPa， ?V （H 2）/ V （油）=500，空速為1.5 h ?-1 反應條件下的DBT轉化率隨時間的變化關系見圖3。各催化劑的HDS活性均先隨反應時間的延長而增大，在反應6 h時，基本達到穩定。由圖3看出，反應6 h后未改性的Ni 2P/M41催化劑的HDS活性為77.5%。而所有稀土金屬改性的 X -Ni 2P/M41催化劑的DBT轉化率在71.5%～86.6%?？梢钥闯?，Pr-Ni 2P/M41和La-Ni 2P/M41催化劑的HDS活性高于未改性Ni 2P/M41催化劑，其中反應6 h后，Pr-Ni 2P/M41催化劑的HDS活性達到92.6%，與未改性的Ni 2P/M41催化劑的DBT轉化率相比，提高15.1%， 且該催化劑經過24 h測試后，活性基本保持不變 （圖4），再次證明Pr-Ni 2P/M41催化劑的穩定性以及對于工業應用的適用性。值得注意的是，Pr-Ni 2P/M41的HDS催化活性也優于當今較為成熟的HDS催化劑以及大部分最新報道的催化劑 （表4）， ?如Ni 2P-TiO 2、NiMo、Ni 2P/Zr-SBA、NiMo/MCM-41等。而Y、Nd和La改性的催化劑，其HDS活性略低于Ni 2P/M41催化劑。 對負載型磷化鎳型催化劑來說，只有稀土Pr和La能提高催化劑的HDS活性，且提高幅度不大，稀土Pr和La反而降低催化劑的HDS活性。這主要是由于載體與稀土的強相互作用導致 ?［28］ 。各催化劑活性順序為Pr-Ni 2P/M41 > La-Ni 2P/M41 > Ni 2P/M41 > Y-Ni 2P/M41 > Nd-Ni 2P/M41 > Ce-Ni 2P/M41。未改性的Ni 2P/M41催化劑的轉化頻率為8.2×10 ?-3 ?s ?-1 ，而 X -Ni 2P/M41催化劑的轉化頻率為2.8×10 ?-3 ～6.5×10 ?-3 ?s ?-1 ，均低于未改性Ni 2P/M41催化劑。這可能是由于嚴重表面P富集導致的。稀土改性負載型磷化鎳催化劑制備及加氫脫硫反應機制見圖5。稀土金屬能夠促進小顆粒催化劑的形成，提高催化劑的比表面積以及分散度，這與TEM中對催化劑顆粒分布情況的結論相吻合。

2.5.2 稀土金屬對催化劑 HDS 選擇性的影響

未改性Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑在340 ℃、3.0 Mpa、 V （H 2）/ ?V （油） =500，空速為1.5 h ?-1 反應條件下的DBT的HDS產物選擇性隨時間的變化見圖6。

由圖6看出，未改性Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑的BP的選擇性均明顯高于CHB的選擇性，表明在各催化劑在HDS反應過程中，DBT主要與催化劑的中的Ni（1）位發生作用，DDS途徑明顯占優勢。未改性Ni 2P/M41和不同稀土金屬改性的 X -Ni 2P/M41催化劑的BP選擇性隨反應時間的延長變化不大，達到穩定后各催化劑的BP選擇性相差不大，為79.1%～80.6%，表明引入稀土金屬不改變催化劑的反應路徑。

3 結 論

（1）向負載型Ni 2P/M41催化劑加入稀土金屬La、Nd、Ce、Y、Pr等，均可以提高催化劑的比表面積，有利于活性相的分散，能促進小晶粒的Ni 2P相的形成。

（2）稀土金屬La和Pr能夠顯著提高Ni 2P/M41催化劑的HDS活性，而加入稀土金屬Nd、Ce和Y，反而使催化劑活性略有降低。

參考文獻 ：

［1］ ?LIU J， WANG Z C， ZHANG D， et al. Systematic engineering on Ni-based nanocatalysts effectively promote hydrogen evolution reaction［J］. Small， 2022，18（13）：2108072.

［2］ ?JIANG N， SHI S J， CUI Y Y， et al. The effect of calcination temperature on the hydrogen evolution reaction performance of Co/NiCoP nano-heterojunction［J］.Journal of Alloys and Compounds， 2022，929：167229.

［3］ ?戈軍，石斌，楊圣闖，等.分散型鎳、鐵催化劑用于重油-煤的加氫共煉研究［J］.中國石油大學學報（自然科學版），2011，35（4）：157-163.

GE Jun， SHI Bin， YANG Shengchuang， et al. Study on dispersed Ni and Fe catalysts for coal heavy oil co processing［J］. Journal of China University of Petroleum（Edition of Natural Science）， 2011，35（4）：157-163.

［4］ ?CAROLINA V， DIEGO G， KENNETH F， et al. Effect of lanthanum doping on the reactivity of unsupported CoMoS 2 catalysts［J］. Applied Catalysis A-General， 2021，611：117891.

［5］ ?CUI T Y， ANTONY R， FAN H X， et al. Review on hydrodesulfurization over zeolite-based catalysts［J］.Industrial & Engineering Chemistry Research， 2021，60（8）：3295-3323.

［6］ ?SONG H， DAI M， SONG H L， et al. A novel synthesis of Ni 2P/MCM-41 catalysts by reducing a precursor of ammonium hypophosphite and nickel chloride at low temperature［J］. Applied Catalysis A：General， 2013，462/463：247-255.

［7］ ?張穎，蘭曉孟，孔垂旸，等.Cu/In 2O 3界面作用增強CO 2催化加氫性能［J］.中國石油大學學報（自然科學版），2023，47（2）：166-172.

ZHANG Ying， LAN Xiaomeng， KONG Chuiyang， et al. Interfacial interaction of Cu/In 2O 3 ???to enhance performance of CO 2 catalytic hydrogenation［J］. Journal of China University of Petroleum （Edition of Natural Science）， 2023，47（2）：166-172.

［8］ ?RIYAJUDDIN S， AZMI K， PAHUJA M， et al. Super-hydrophilic hierarchical Ni-foam-graphene-carbon nanotubes-Ni 2P-CuP 2 nano-architecture as efficient electrocatalyst for overall water splitting［J］. ACS Nano， 2021，15：5586-5599.

［9］ ?SONG H， ZHANG F Y， SONG H L， et al. The effect of neodymium content on dibenzothiophene HDS performance over a bulk Ni 2P catalyst［J］. Catalysis Communications， 2015，69：59-62.

［10］ ?JIANG N， SHI S J， CUI Y Y， et al. Effect of phosphorization temperature on the structure and hydrogen evolution reaction performance of nickel cobalt phosphide electrocatalysts［J］.Catalysis Communications， 2022，171：106507.

［11］ ?ZHANG X Y， GUO B Y， CHEN Q W， et al. Ultrafine and highly dispersed bimetal Ni 2P/Co 2P encapsulated by hollow N-doped carbon nanospheres for efficient hydrogen evolution［J］. International Journal of Hydrogen Energy， 2019，44（29）：14908-14917.

［12］ ?WANG G， SHI Y， MEI J L， et al. DFT insights into hydrogen activation on the doping Ni 2P surfaces under the hydrodesulfurization condition［J］. Applied Surface Science， 2021，538：148160.

［13］ ?JIANG N， JIANG B L， WANG S， et al. Efficient Ni 2P /Al 2O 3 hydrodesulfurization catalysts from surface hybridization of Al 2O 3 particles with graphite-like carbon［J］. Journal of the Taiwan Institute of Chemical Engineers， 2021，121：139-146.

［14］ ?ABU I， SMITH K J. The effect of cobalt addition to bulk MoP and Ni 2P catalysts for the hydrodesulfurization of 4，6-dimethyldibenzothiophene ［J］. Journal of Catalysis， 2006，241：356-366.

［15］ ?SONG H， XU X W， SONG H L， et al. Synthesis of an yttrium-modified bulk Ni 2P catalyst with high hydrodesulfurization activity［J］. Catalysis Communications， 2015，63：52-55.

［16］ ?WANG G， WU C，QIANG Z， SUN X， FANG W. Effects of lanthanum and ethylenediamine tetraacetic acid on structure and hydrodesulfurization performance of Ni-W/TiO 2-Al 2O 3 catalyst［J］. Industrial Catalysis， 2017，25：53-57.

［17］ ?ISLAM A， ABDULRAHMAN A A， TAWIFIK A S， et al. Carbon nanofiber-doped zeolite as support for molybdenum based catalysts for enhanced hydrodesulfurization of dibenzothiophene ［J］. Journal of Molecular Liquids， 2020，304：112376.

［18］ ?JIANG N， ZHANG F Y， SONG H. Effect of reduction temperature on the structure and hydrodesulfurization performance of Na doped Ni 2P/MCM-41 catalysts［J］. RSC Advances， 2019，9，15488-15494.

［19］ ?LIN R Y， PAN H D， XU W D， et al. Hydrodesulfurization of benzothiophene on Ni 2P surface［J］. Energy Exploration & Exploitation， 2020，38（6）：2711-2728.

［20］ ?TOPALIAN P J， CARRILLO B A， COCHRAN P M， et al. Synthesis and hydrodesulfurization properties of silica-supported nickel-ruthenium phosphide catalysts［J］.Journal of Catalysis， 2021，403：173-180.

［21］ ??VASILEIOS， PILLAI V， GIOTI C， et al. Ultrafine Ni 2P nanoparticle-decorated r-GO： a novel liquid-phase approach and dibenzothiophene hydrodesulfurization［J］. Industrial & Engineering Chemistry Research， 2021，60（11）：4300-4309.

［22］ ?LI X， SUN Z C， WANG A J， et al. Effect of TiO 2 on the hydrodesulfurization performance of bulk Ni 2P［J］. Applied Catalysis A-General， 2012，417：19-25.

［23］ ?SONG H， YU Q， JIANG N， et al. Effect of surface modification temperature on the hydrodesulfurization performance of Ni 2P/MCM-41 catalyst［J］. Research on Chemical Intermediates， 2018，44：3629-3640.

［24］ ?MA H Q， LI Q R， SHI Y， et al. Ni 2P/ZrO 2-SBA-15 dibenzothiophene hydrodesulfurization catalysts： preparation， characterization and evaluation［J］. Transactions of Tianjin University， 2018，24：340-350.

［25］ ?PUELLO-POLO E， MARQUEZ E， BRITO J L， et al. One-pot synthesis of Nb-modified Al 2O 3 support for NiMo hydrodesulfurization catalysts［J］. Journal of Sol-Gel Science and Technology， 2018，88：90-99.

［26］ ??FRANKLIN J M， OSCAR E F L， XIM B， et al. Dibenzothiophene hydrodesulfurization with NiMo and CoMo catalysts supported on niobium-modified MCM-41［J］. Applied Catalysis B-Environmental， 2017，219：479-491.

［27］ ?LIU H， YIN C L， LI X H， et al. Effect of NiMo phases on the hydrodesulfurization activities of dibenzothiophene［J］. Catalysis Today， 2017，282：222-229.

［28］ ?WANG G， WANG A C， WANG E H， et al. DFT insights into competitive adsorption and reaction mechanism of benzothiophene and naphthalene on Fe-doped Ni 2P catalyst［J］. Fuel， 2022，314：123114.

（編輯 劉為清）