CN113035285A - 微气泡尺寸对油品加氢脱硫效果影响的计算方法 - Google Patents
微气泡尺寸对油品加氢脱硫效果影响的计算方法 Download PDFInfo
- Publication number
- CN113035285A CN113035285A CN202110354325.XA CN202110354325A CN113035285A CN 113035285 A CN113035285 A CN 113035285A CN 202110354325 A CN202110354325 A CN 202110354325A CN 113035285 A CN113035285 A CN 113035285A
- Authority
- CN
- China
- Prior art keywords
- reactor
- hydrodesulfurization
- reaction
- size
- concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000000694 effects Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000006243 chemical reaction Methods 0.000 claims abstract description 55
- 238000012546 transfer Methods 0.000 claims abstract description 30
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 22
- 239000001257 hydrogen Substances 0.000 claims abstract description 22
- 229910052739 hydrogen Inorganic materials 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 239000007791 liquid phase Substances 0.000 claims abstract description 17
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 239000000376 reactant Substances 0.000 claims abstract description 8
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 7
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 7
- 239000001301 oxygen Substances 0.000 claims abstract description 7
- 239000007789 gas Substances 0.000 claims description 23
- 238000006477 desulfuration reaction Methods 0.000 claims description 17
- 230000023556 desulfurization Effects 0.000 claims description 17
- 239000007788 liquid Substances 0.000 claims description 14
- 239000003054 catalyst Substances 0.000 claims description 10
- UCKMPCXJQFINFW-UHFFFAOYSA-N Sulphide Chemical compound [S-2] UCKMPCXJQFINFW-UHFFFAOYSA-N 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 6
- 238000005984 hydrogenation reaction Methods 0.000 claims description 6
- 239000012535 impurity Substances 0.000 claims description 6
- 150000001335 aliphatic alkanes Chemical class 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 4
- 150000004767 nitrides Chemical class 0.000 claims description 3
- 239000012071 phase Substances 0.000 claims description 3
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 claims description 3
- 239000011949 solid catalyst Substances 0.000 claims description 3
- 150000003568 thioethers Chemical class 0.000 claims description 2
- 239000000470 constituent Substances 0.000 claims 1
- 239000000047 product Substances 0.000 description 20
- 230000008569 process Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000005501 phase interface Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000005728 strengthening Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000005587 bubbling Effects 0.000 description 2
- 238000009903 catalytic hydrogenation reaction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000012495 reaction gas Substances 0.000 description 2
- -1 sulfide compound Chemical class 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000002283 diesel fuel Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002779 inactivation Effects 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 239000002574 poison Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
-
- C—CHEMISTRY; METALLURGY
- C10—PETROLEUM, GAS OR COKE INDUSTRIES; TECHNICAL GASES CONTAINING CARBON MONOXIDE; FUELS; LUBRICANTS; PEAT
- C10G—CRACKING HYDROCARBON OILS; PRODUCTION OF LIQUID HYDROCARBON MIXTURES, e.g. BY DESTRUCTIVE HYDROGENATION, OLIGOMERISATION, POLYMERISATION; RECOVERY OF HYDROCARBON OILS FROM OIL-SHALE, OIL-SAND, OR GASES; REFINING MIXTURES MAINLY CONSISTING OF HYDROCARBONS; REFORMING OF NAPHTHA; MINERAL WAXES
- C10G45/00—Refining of hydrocarbon oils using hydrogen or hydrogen-generating compounds
- C10G45/02—Refining of hydrocarbon oils using hydrogen or hydrogen-generating compounds to eliminate hetero atoms without changing the skeleton of the hydrocarbon involved and without cracking into lower boiling hydrocarbons; Hydrofinishing
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Computing Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Organic Chemistry (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明提供了一种微气泡尺寸对油品加氢脱硫效果影响的计算方法,该计算方法包括如下步骤:(A)构建加氢脱硫宏观反应动力学方程,(B)当不考虑氢气参与的其它反应时,依据氧气传质速率和反应消耗速率相等的关系,建立传质平衡方程,(C)构建液相反应物浓度数学模型,(D)建立所述方程式(2)中的PA与FA之间的关系式,通过联立上述方程式,在已知反应器入口处浓度CA0、CB0,计算得到微气泡大小一定时体系实际CA以及CB的大小以体现微气泡大小对加氢脱硫反应效果的影响。本发明的计算方法确定微界面强化固定床油品加氢脱硫效果的计算方法,以优化固定床油品加氢脱硫微界面强化反应器的设计和操作。
Description
技术领域
本发明涉及微界面强化领域,具体而言,涉及一种微气泡尺寸对油品加氢脱硫效果影响的计算方法。
背景技术
目前柴油脱硫使用最为普遍和有效的就是加氢脱硫技术。其中以固定床鼓泡反应器作为反应器的催化加氢技术产品质量高,应用技术成熟,目前应用最为广泛。该工艺以氢气为反应气体原料,Co和Mo等催化剂,操作温度300~420℃,脱硫率86~99.94%,液时空速1~14h-1。尽管该工艺持液量大、液固接触效率高、液体分布均匀性好,减少了催化剂的结焦和失活。但仍存在一些不足,为使原料油携带足够反应所需的氢气,常规的加氢技术需要通过加压或物料循环的方式提高氢气溶解度,由此产生了氢耗大,设备费用高昂,工艺要求严苛等问题。
微界面强化反应器气液相界面积构效调控数学模型建模方法(公开号:CN109684769A)量化了反应器气泡气液相界面积与反应器的结构参数、操作参数以及物性参数之间的关系。当通气量一定时,气泡越小,气液相界面积越大,传质性能更好。但气泡尺寸如何影响脱硫过程,如何获得最优于反应过程的气泡大小,还未有具体的研究。当脱硫要求一定时,如何优化反应器内气泡尺度以最大程度实现微界面强化技术的经济性是一个重要的现实问题。
有鉴于此,特提出本发明。
发明内容
有鉴于此,本发明提供一种微气泡尺寸对油品加氢脱硫效果影响的计算方法,该计算方法依据实际加氢脱硫反应情况为出发点,确定微界面强化固定床油品加氢脱硫效果的计算方法,以优化固定床油品加氢脱硫微界面强化反应器的设计和操作。
决定微界面强化固定床油品加氢脱硫反应体系宏观反应速率的因素包括两类,一类是传质推动力,一般为气泡内反应气体(氢气)的分压;另一类是包括传质阻力和本征反应阻力在内的两类反应阻力。对于多相反应体系,传质阻力一般是决定体系宏观反应速率大小的主因。微界面强化技术的优势在于,通过微界面机组在体系中形成大量的微气泡,使体系中气液相界面积大幅提升,从而极大减小传质阻力,最终实现反应的充分强化。
在油品的加氢脱硫过程中,由于氢气是难溶气体,在传质过程中受液膜阻力的影响,限制了加氢过程的反应速率。气泡尺度减小有利于传质阻力的减小,使得传质速率获得大幅加速或强化,从而全部或部分消除了由于宏界面体系的相间传质速率偏低而造成的传质瓶颈之故,加氢效率能够大大提升。
本发明探讨了气泡大小对固定床油品加氢脱硫效果的计算方法,当然对于微界面强化反应器气泡尺度的相关建模方法,之前的专利也多有涉及,比如专利公开号CN109684769A、CN107346378A、CN107335390A、CN107561938A等等,所以虽然本发明的整个计算方法中并没有涉及到气泡尺寸d32的计算,但是平衡方程式中以及其他诸多计算公式所涉及到的气液相界面积a、液侧传质系数kL、气侧传质系数kG、气含率QG、宏观反应速率rB等参数均与d32有一定的关联性,只不过相关公式属于现有技术已经在之前专利涉及到,所以通过本发明的计算方法仍然脱离不开与气泡尺寸d32的关联性,本发明是在之前诸多专利的基础之上进行再创新,以探讨微气泡尺寸对加氢脱硫效果的影响,比如XB、rB、CB、气体利用率FA的计算方法,这些参数可以很好的反应出微气泡尺寸对脱硫效果如何影响,并且这些计算方法均是之前专利所没有涉及过的。
本发明提供的微气泡尺寸对油品加氢脱硫效果影响的计算方法,具体包括如下步骤:
(B)当不考虑氢气参与的其它反应时,依据氧气传质速率和反应消耗速率相等的关系,可得如下传质平衡方程式:
kGa(PA/HA-CAi)=kLa(CAi-CA)=kSaS(CA-CAS)=(1-φG)(-rA) (2),A~D分别代表H2、硫化物、相应的脱硫产物烷烃和H2S,它们在液相中的浓度分别以CA、CB、Cc、CD来表示;
(C)构建液相反应物浓度数学模型:假定反应器可近似为平推流反应器,反应物A、B在固体催化剂表面进行反应,A和B的进料量分别为FB0和FA0,单位为mol/s,反应器内有效体积为V/m3,催化剂床层高度为Lb/m;
依据平推流反应器数学模型,根据硫化物B的物料平衡关系可以推导出下式:
FB0dXB=-rBdV (3)
XB为脱硫率,及dV=Adx=S0εdx,FB=FB0(1-XB),可推导得到FA~FD在催化剂床层不同高度处变化情况的数学表达式分别为:
其中,rA=-kACAS mCBS n,rB=-kBCAS mCBS n,rA=2rB,且FB=QBCB,则可推导得到CB的数学表达式为:
边界条件:
当x=0时,CA=CA0,CB=CB0,CC=0,CD=0
当x=Lb时,CA=CAf,CB=CBf,CC=CCf,CD=CDf;
(D)建立所述方程式(2)中的PA与FA之间的关系式,具体公式如下:
其中,PT为操作压力,PA为氢气分压,FA、FD分别是氢气与气相产物的气体流量,通过联立上述方程式,在已知反应器入口处浓度CA0、CB0,计算得到微气泡大小一定时体系实际CA以及CB的大小以体现微气泡大小对加氢脱硫反应效果的影响;
上述公式中,HA-亨利系数,MPa.m3/mol;
P-操作压力,atm;
a-气-液相界面积,m2/m3;
aS-液-固相界面积,m2/m3;
n,m-反应级数;
F-摩尔流量,mol/s;
Q-体积流量,m3/s;
P-压力,Pa;
V-反应器体积,m3;
r-反应速率,mol·m-3·s-1;
X-转化率,%;
C-组分浓度,mol/m3;
k-反应速率常数,s-1;
kG、kL、kS——气侧,液侧、固侧传质系数,m/s
S0-反应器横截面积,m2;
ε-床层空隙率;
在本发明的计算方法中,先通过步骤(A)构建宏观动力学方程,然后通过后续的步骤(B)建立氧气传质速率与反应消耗速率之间的等式关系,以得到对于整个反应来说最佳的微气泡尺寸,因为通过大量实践发现,只有微气泡尺寸大小正好满足氧气传质速率等于氧气在液相的消耗速率时,其反应状态是最佳的,然后根据所确定的微气泡尺寸,来计算XB、rB、CB、气体利用率FA等,如何确定微气泡尺寸的方法之前的专利已有涉及,属于现有技术,本发明的关键改进点在于如何利用微气泡尺寸来计算能够表征脱硫效果的XB、rB、CB、气体利用率FA等参数。
为了计算CB,本发明在步骤(C)中构建了液相反应物浓度数学模型,为了计算气体利用率FA,本发明在步骤(D)中构建了PA与FA之间的关系式。
优选地,作为进一步可实施的方案,当油品中除了硫化物,含还有其他杂质时,公式(4)-(7)的数学关系采用以下多组分方程式(10)-(13)代替,具体计算方式如下步骤:
多组分加氢脱除效果方程组的建立;
γH2+R=X→R+ηXHn;
多组分的加氢脱硫数学关系如下所示:
其中,Bi(i=1,2,3…)依次分别指硫化物、氮化物与多环芳烃杂质,Ci指相应的烷烃产物,Di指相应的气体产物。
因为在石油中,除了硫化物之外,氮化物、多环芳烃等杂质同样是石油中很重要的一类非烃组分,它既能使某些二次加工过程的催化剂中毒,也影响某些石油产品的安定性。同时也能够被催化加氢脱除,因此为了扩大本发明方案的适用面,本发明的上述公式适用于多组分加氢脱硫的反应情况。
本发明还涉及了一种反应器,该反应器采用上述计算方法设计得到。通过上述计算方法设计得到的反应器,更符合实际应用情况,能够将微气泡尺寸控制到最理想的状态,以达到良好的微界面反应效果。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为不同氢油比条件下气泡尺寸d32对气液相界面积a的影响;
图2为不同氢油比条件下气泡尺寸d32对液侧传质系数kL的影响;
图3为不同氢油比条件下气泡尺寸d32对气侧传质系数kG的影响;
图4为不同氢油比条件下气泡尺寸d32对气含率QG的影响;
图5为不同氢油比条件下气泡尺寸d32对宏观反应速率rB的影响;
图6为不同氢油比条件下气泡尺寸d32对脱硫率XB的影响。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
在本公开使用的术语是仅仅出于描述特定实施例的目的,而非旨在限制本公开。在本公开和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本公开可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本公开范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
下面结合具体的实施例对本发明的技术方案做进一步阐述。
本实施例基于本发明的建模方法,针对某企业柴油加氢脱硫的固定床鼓泡反应器及现有的操作工况,研究微界面机组在纯气动条件下反应器内气泡大小对硫化物宏观反应速率的影响。反应器内气泡尺度受反应器结构、物性参数、操作条件影响,其在纯气动操作条件下的调控模型可参考申请人的在先专利CN109684769A。
计算条件如下:
反应器设计参数:
反应器高度H0=8.2m;反应器横截面积S0=2.0096m2;
液相密度ρL=840kg/m3;液相的摩尔质量M1=0.1844kg/mol;
操作压力Pm=3.6MPa;
柴油流量QL=17.857t/h;氢油比Qoil/QH2=260Nm3 H2/m3 oil。
计算过程包括:
(B)当不考虑氢气参与的其它反应时,依据氧气传质速率和反应消耗速率相等的关系,可得如下传质平衡方程式:
kGa(PA/HA-CAi)=kLa(CAi-CA)=kSaS(CA-CAS)=(1-φG)(-rA) (2),A~D分别代表H2、硫化物、相应的脱硫产物烷烃和H2S,它们在液相中的浓度分别以CA、CB、Cc、CD来表示;
(C)构建液相反应物浓度数学模型:假定反应器可近似为平推流反应器,反应物A、B在固体催化剂表面进行反应,A和B的进料量分别为FB0和FA0,单位为mol/s,反应器内有效体积为V/m3,催化剂床层高度为Lb/m;
依据平推流反应器数学模型,根据硫化物B的物料平衡关系可以推导出下式:
FB0dXB=-rBdV (3)
XB为脱硫率,及dV=Adx=S0εdx,FB=FB0(1-XB),可推导得到FA~FD在催化剂床层不同高度处变化情况的数学表达式分别为:
其中,rA=-kACAS mCBS n,rB=-kBCAS mCBS n,rA=2rB,且FB=QBCB,则可推导得到CB的数学表达式为:
边界条件:
当x=0时,CA=CA0,CB=CB0,CC=0,CD=0
当x=Lb时,CA=CAf,CB=CBf,CC=CCf,CD=CDf;
(D)建立所述方程式(2)中的PA与FA之间的关系式,具体公式如下:
其中,PT为操作压力,PA为氢气分压,FA、FD分别是氢气与气相产物的气体流量,通过联立上述方程式,在已知反应器入口处浓度CA0、CB0,计算得到微气泡大小一定时体系实际CA以及CB的大小以体现微气泡大小对加氢脱硫反应效果的影响。
另外,本发明对于油品加氢脱硫工艺具体微气泡尺寸与各个关键脱硫效果参数的影响关系进行了模拟,具体结果参见图1~6所示。
图1为不同氢油比条件下气泡尺寸d32对气液相界面积a的影响,具体关系式为:该公式之前专利已有涉及,由图1可知,当氢油比条件一定时,体系气泡尺寸的减小有利于气液相界面积的增加。当然,这也与体系气含率增大有关。
图2为不同氢油比条件下气泡尺寸d32对液侧传质系数kL的影响,关系式为该公式之前专利已有涉及,图3为不同氢油比条件下气泡尺寸d32对气侧传质系数kG的影响,关系式为该公式之前专利已有涉及,由图2-3可知,体系气侧和液侧体积传质系数均随气泡直径的减小而增大;尤其是当气泡尺度由毫米级减小至微米级时,该变化趋势较为显著。因此,减小体系气泡大小有利于体系体积传质系数的增加。
图4为不同氢油比条件下气泡尺寸d32对气含率QG的影响,关系式为该公式之前专利已有涉及,图4表明,在其它条件不变时,当气泡由毫米级减小至微米级时,体系气含率将大幅增加。在实际生产过程中,反应器内气含率不应过高。一般情况下,反应器内气含率不应超过0.5,由图4可知此时反应器内气泡直径不应小于0.3mm。
图5表明,当体系气泡尺度由毫米级减小至微米级时,不同氢油比条件下的宏观反应速率增大的趋势均呈先显著后缓慢、宏观反应速率差异逐渐小的特征。因此,当将气泡尺度减小至微米级时,反应器内中间产物的累积量减少,反应选择性提高。
由图6可知,当其它条件不变时,硫化物转化率随体系中气泡尺度的减小而增大。在给定的条件下,若要达到国六标准,即产品硫含量小于10ppm,氢油比为100、150、200、250、300时d32分别不能大于4.7mm、6.3mm、7.6mm、8.7mm、9.8mm。这表明,在目前的条件下,单就提高脱硫率而言,将气泡尺度过度减小无实际意义。
虽然在本发明的方案中并没有涉及到d32,但是可见方案中出现的诸多重要参数均是与d32有一定关联性的,只是因为与d32的计算关系之前专利均涉及过,因此该专利方案中并没有过多涉及,本发明的方案重点在于如何利用d32来对加氢脱硫效果的关键参数进行计算,以达到对反应器优化设计的目的,尤其是如何针对特定的油品加氢脱硫工艺的效果进行评价的方法,是现有技术中所未涉及过的。
由图1~6的结果可知,当给定生产指标(如脱硫率、氢气利用率等)时,采用本方案所采用的计算方法可以确定固定床油品加氢脱硫微界面反应体系中最经济的气泡大小。理论计算结果表明,当其它条件一定时,减小体系气泡大小有利宏观反应速率及脱硫率、氢气利用率的提升。当气泡减小至一定数值后,继续减小气泡尺度对于指标的进一步改善作用有限。
以上所述仅为本公开的较佳实施例而已,并不用以限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开保护的范围之内。
Claims (3)
1.一种微气泡尺寸对油品加氢脱硫效果影响的计算方法,其特征在于,包括如下步骤:
(B)当不考虑氢气参与的其它反应时,依据氧气传质速率和反应消耗速率相等的关系,可得如下传质平衡方程式:
kGa(PA/HA-CAi)=kLa(CAi-CA)=ksas(CA-CAS)=(1-φG)(-rA) (2),A~D分别代表H2、硫化物、相应的脱硫产物烷烃和H2S,它们在液相中的浓度分别以CA、CB、Cc、CD来表示;
(C)构建液相反应物浓度数学模型:假定反应器可近似为平推流反应器,反应物A、B在固体催化剂表面进行反应,A和B的进料量分别为FB0和FA0,单位为mol/s,反应器内有效体积为V,单位为m3,催化剂床层高度为Lb,单位为m;
依据平推流反应器数学模型,根据硫化物B的物料平衡关系可以推导出下式:
FB0dXB=-rBdV (3)
XB为脱硫率,及dV=Adx=S0εdx,FB=FB0(1-XB),可推导得到FA~FD在催化剂床层不同高度处变化情况的数学表达式分别为:
其中,rA=-kACAS mCBS n,rB=-kBCAS mCBS n,rA=2rB,且FB=QBCB,则可推导得到CB的数学表达式为:
边界条件:
当x=0时,CA=CA0,CB=CB0,CC=0,CD=0
当x=Lb时,CA=CAf,CB=CBf,CC=CCf,CD=CDf;
(D)建立所述方程式(2)中的PA与FA之间的关系式,具体公式如下:
其中,PT为操作压力,PA为氢气分压,FA、FD分别是氢气与气相产物的气体流量;
通过联立上述(1)-(9)方程式,在已知反应器入口处浓度CA0、CB0,计算得到微气泡大小一定时体系实际CA以及CB的大小以体现微气泡大小对加氢脱硫反应效果的影响;
上述公式中,HA-亨利系数,MPa〃m3/mol;
P-操作压力,atm;
a-气-液相界面积,m2/m3;
aS-液-固相界面积,m2/m3;
n,m-反应级数;
F-摩尔流量,mol/s;
Q-体积流量,m3/s;
P-压力,Pa;
V-反应器体积,m3;
r-反应速率,mol〃m-3〃s-1;
X-转化率,%;
C-组分浓度,mol/m3;
k-反应速率常数,s-1;
kG、kL、kS——气侧,液侧、固侧传质系数,m/s
S0-反应器横截面积,m2;
ε-床层空隙率;
3.一种反应器,采用权利要求1-2任一项所述的计算方法设计得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110354325.XA CN113035285B (zh) | 2021-04-01 | 2021-04-01 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110354325.XA CN113035285B (zh) | 2021-04-01 | 2021-04-01 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113035285A true CN113035285A (zh) | 2021-06-25 |
CN113035285B CN113035285B (zh) | 2023-04-11 |
Family
ID=76453588
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110354325.XA Active CN113035285B (zh) | 2021-04-01 | 2021-04-01 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113035285B (zh) |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107335390A (zh) * | 2017-08-30 | 2017-11-10 | 南京大学 | 微界面强化反应器相界面积构效调控模型建模方法 |
CN107346378A (zh) * | 2017-08-30 | 2017-11-14 | 南京大学 | 微界面强化反应器传质速率构效调控模型建模方法 |
CN107561938A (zh) * | 2017-08-30 | 2018-01-09 | 南京大学 | 微界面强化反应器反应速率构效调控模型建模方法 |
CN107563051A (zh) * | 2017-08-30 | 2018-01-09 | 南京大学 | 微界面强化反应器气泡尺度构效调控模型建模方法 |
CN107589667A (zh) * | 2017-08-30 | 2018-01-16 | 南京大学 | 微界面强化反应器能量耗散构效调控模型建模方法 |
CN109657406A (zh) * | 2019-01-09 | 2019-04-19 | 南京大学 | Miha纯气动操作条件下相界面积调控模型建模方法 |
CN109684769A (zh) * | 2019-01-09 | 2019-04-26 | 南京大学 | Miha纯气动操作条件下气泡尺度调控模型建模方法 |
CN109887550A (zh) * | 2019-01-09 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下能量耗散调控模型建模方法 |
CN109887551A (zh) * | 2019-02-04 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下传质调控模型建模方法 |
CN110245434A (zh) * | 2019-06-18 | 2019-09-17 | 南京延长反应技术研究院有限公司 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
CN111068588A (zh) * | 2018-10-22 | 2020-04-28 | 中国石油化工股份有限公司 | 一种柴油超深度脱硫装置及柴油加氢反应*** |
CN111359542A (zh) * | 2019-03-15 | 2020-07-03 | 南京延长反应技术研究院有限公司 | 一种微界面强化柴油加氢精制反应***及方法 |
CN111482144A (zh) * | 2019-01-29 | 2020-08-04 | 南京延长反应技术研究院有限公司 | 用于渣油加氢反应的下置式微界面强化反应装置及方法 |
CN111482145A (zh) * | 2019-01-29 | 2020-08-04 | 南京延长反应技术研究院有限公司 | 用于渣油加氢反应的上置式微界面强化反应装置及方法 |
CN111523194A (zh) * | 2019-02-04 | 2020-08-11 | 南京大学 | Miha纯气动操作条件下反应速率及转化率调控模型建模方法 |
CN111686646A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种汽油加氢脱硫改质微界面强化反应***及方法 |
CN111686654A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种煤焦油加氢微界面乳化床强化反应*** |
CN111686643A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种上行式固定床油品加氢微界面强化反应*** |
CN111686644A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种下行式固定床油品加氢微界面强化反应*** |
CN111871338A (zh) * | 2020-07-16 | 2020-11-03 | 南京延长反应技术研究院有限公司 | 一种柴油加氢的智能微界面反应***及方法 |
-
2021
- 2021-04-01 CN CN202110354325.XA patent/CN113035285B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107346378A (zh) * | 2017-08-30 | 2017-11-14 | 南京大学 | 微界面强化反应器传质速率构效调控模型建模方法 |
CN107561938A (zh) * | 2017-08-30 | 2018-01-09 | 南京大学 | 微界面强化反应器反应速率构效调控模型建模方法 |
CN107563051A (zh) * | 2017-08-30 | 2018-01-09 | 南京大学 | 微界面强化反应器气泡尺度构效调控模型建模方法 |
CN107589667A (zh) * | 2017-08-30 | 2018-01-16 | 南京大学 | 微界面强化反应器能量耗散构效调控模型建模方法 |
CN107335390A (zh) * | 2017-08-30 | 2017-11-10 | 南京大学 | 微界面强化反应器相界面积构效调控模型建模方法 |
CN111068588A (zh) * | 2018-10-22 | 2020-04-28 | 中国石油化工股份有限公司 | 一种柴油超深度脱硫装置及柴油加氢反应*** |
CN109657406A (zh) * | 2019-01-09 | 2019-04-19 | 南京大学 | Miha纯气动操作条件下相界面积调控模型建模方法 |
CN109684769A (zh) * | 2019-01-09 | 2019-04-26 | 南京大学 | Miha纯气动操作条件下气泡尺度调控模型建模方法 |
CN109887550A (zh) * | 2019-01-09 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下能量耗散调控模型建模方法 |
CN111482145A (zh) * | 2019-01-29 | 2020-08-04 | 南京延长反应技术研究院有限公司 | 用于渣油加氢反应的上置式微界面强化反应装置及方法 |
CN111482144A (zh) * | 2019-01-29 | 2020-08-04 | 南京延长反应技术研究院有限公司 | 用于渣油加氢反应的下置式微界面强化反应装置及方法 |
CN109887551A (zh) * | 2019-02-04 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下传质调控模型建模方法 |
CN111523194A (zh) * | 2019-02-04 | 2020-08-11 | 南京大学 | Miha纯气动操作条件下反应速率及转化率调控模型建模方法 |
CN111359542A (zh) * | 2019-03-15 | 2020-07-03 | 南京延长反应技术研究院有限公司 | 一种微界面强化柴油加氢精制反应***及方法 |
CN111686646A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种汽油加氢脱硫改质微界面强化反应***及方法 |
CN111686654A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种煤焦油加氢微界面乳化床强化反应*** |
CN111686643A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种上行式固定床油品加氢微界面强化反应*** |
CN111686644A (zh) * | 2019-03-15 | 2020-09-22 | 南京延长反应技术研究院有限公司 | 一种下行式固定床油品加氢微界面强化反应*** |
CN110245434A (zh) * | 2019-06-18 | 2019-09-17 | 南京延长反应技术研究院有限公司 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
CN111871338A (zh) * | 2020-07-16 | 2020-11-03 | 南京延长反应技术研究院有限公司 | 一种柴油加氢的智能微界面反应***及方法 |
Non-Patent Citations (7)
Title |
---|
E.RODRI´GUEZ-CASTELLO´N等: "Nickel and cobalt promoted tungsten and molybdenum sulfide mesoporous catalysts for hydrodesulfurization", 《SCIENCEDIRECT》 * |
M.A.RODRIGUEZ等: "Modeling of Hydrodesulfurization (HDS),Hydrodenitrogenation (HDN), and the Hydrogenation of Aromatics (HDA) in a Vacuum Gas Oil Hydrotreater", 《ENERGY & FUELS》 * |
SU WANG等: "Catalyst-free biphasic oxidation of Thiophenes in continuous-flow", 《JOURNAL OF FLOW CHEMISTRY》 * |
孙国权等: "柴油加氢脱硫三集总动力学", 《石油学报(石油加工)》 * |
张志炳等: "多相反应体系的微界面强化简述", 《化工学报》 * |
张志炳等: "气液反应体系相界面传质强化研究", 《化学工程》 * |
田洪舟等: "微界面强化重油浆态床低压加氢的传质基础", 《化学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113035285B (zh) | 2023-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2261172B1 (en) | Process for producing liquid hydrocarbon from natural gas | |
US5817701A (en) | Slurry hydrocarbon synthesis with cyclic CO purge and catalyst rejuvenation | |
CN102382672B (zh) | 一种合成液态烃的方法 | |
US20060167118A1 (en) | Modification of a methanol plant for converting natural gas to liquid hydrocarbons | |
WO2005054163A1 (en) | Control of co2 emissions from a fischer-tropsch facility by use of dual functional syngas conversion | |
US8686206B2 (en) | Single loop multistage fuel production | |
SA99200195B1 (ar) | تحويل غاز باستخدام هيدروجين ناتج من غاز تخليق لإزالة الكبريت من سوائل هيدروكربونية لبئر غاز | |
CN103289738B (zh) | 一种加氢裂化尾油加氢生产高档润滑油基础油的组合方法 | |
CN103429710B (zh) | 在蒸馏塔内的三个催化反应区上的粗汽油异构化 | |
WO2014188243A1 (en) | Methods and apparatus for three phase contacting and reactions in a cross flow reactor | |
US11247897B2 (en) | Base oil production via dry reforming | |
US20080237090A1 (en) | Process and system for redcuing the olefin content of a hydrocarbon feed gas and production of a hydrogen-enriched gas therefrom | |
EA005789B1 (ru) | Процесс фишера-тропша | |
US6914082B2 (en) | Slurry bubble reactor operated in well-mixed gas flow regime | |
US2762683A (en) | Contacting of solid material with liquid phase reactant streams | |
CN105820839B (zh) | 一种轻烃异构化方法 | |
CN113035285B (zh) | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 | |
US6921778B2 (en) | Process for converting synthesis gas in reactors that are arranged in series | |
CN112143521B (zh) | 一种生产催化重整原料的加氢方法和*** | |
US7300567B2 (en) | Hydrotreating process | |
CN102165039B (zh) | 液体燃料合成方法及液体燃料合成装置 | |
CN205347269U (zh) | 一种环保发泡剂与溶剂油联产装置 | |
Eldarsi et al. | Methyl-Tert-Butyl-Ether Catalytic Distillation Column: Part II: Optimization | |
US20040235968A1 (en) | Process and apparatus for controlling flow in a multiphase reactor | |
CN101942341A (zh) | 一种石油蜡类物质加氢改质方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |