2026西安石油大学函授《化工热力学》作业1

2026西安石油大学函授《化工热力学》作业1《化工热力学》作业1（2026学年西安石油大学函授专用）一、核心名词解释（覆盖第1单元流体pVT性质核心考点）1.对应状态原理：不同流体在两个对比状态参数相等时，第三个对比状态参数也近似相等，此时所有流体处于热力学相似状态。该原理19世纪末由范德华提出，后续经皮策等人修正后成为普遍化热力学性质计算的核心基础，针对石油组分这类非极性弱极性流体的计算精度可控制在3%以内，完全满足油气集输、炼化单元的工程计算要求。对应状态原理的本质是流体分子间的交互作用力仅和对比状态参数相关，当不同流体的对比温度Tr（实际温度与临界温度的比值）、对比压力Pr（实际压力与临界压力的比值）相同时，分子间距离、分子间势能的无量纲化偏差可以忽略。2.偏心因子：1955年皮策提出的用于表征非简单流体分子极性、非球形度的无量纲参数，定义为ω=-lg(Psr)_{Tr=0.7}-1.00，其中(Psr)_{Tr=0.7}是对比温度等于0.7时流体的饱和蒸气压与临界压力的比值。对于氩、氪、氙这类球形结构的简单单原子流体，ω=0，其热力学性质仅由对比温度、对比压力即可完全确定；而甲烷的偏心因子为0.011、丙烷为0.152、正辛烷为0.398，偏心因子数值越大说明分子非球形度越高、极性越强，需在普遍化性质计算中引入修正项才能保证精度。3.余焓：也称焓的偏差函数，定义为相同温度、相同压力条件下，实际流体的摩尔焓H与假设该流体处于理想气体状态下的摩尔焓H^ig的差值，即H^R=H-H^ig。余焓表征了实际流体因为分子间存在相互作用力、分子本身具有体积而偏离理想气体模型的焓值偏差，是工程上计算高压流体焓值的核心参数，直接应用于炼化工艺中绝热节流、膨胀做功、换热过程的能量衡算环节。4.稳定流动过程：化工热力学中针对连续生产工艺定义的基准过程，指流动系统中任意截面上流体的热力学性质（温度、压力、比容、流速、相态）都不随时间发生变化的过程，该过程满足物料衡算的累积量为0、能量衡算的系统内能变化率为0的条件。油气炼化、天然气净化的绝大多数连续运行单元（换热器、压缩机、透平机、精馏塔）都可以按稳定流动过程做简化计算，仅在装置开停车、故障波动阶段需要做非稳态修正。5.维里方程：从统计力学分子交互作用理论出发推导得到的状态方程，形式为PV_m/RT=1+B/V_m+C/V_m²+D/V_m³+…，其中B为第二维里系数，表征两个流体分子之间的交互作用势能，仅为温度的函数；C为第三维里系数，表征三个流体分子之间的碰撞交互作用，同样仅和温度相关。维里方程是所有状态方程中唯一具备明确物理意义的形式，并非经验拟合模型，中低压条件下截断至二阶、三阶即可实现很高的计算精度。6.饱和蒸气压：指定温度下，纯物质的气液两相达到相平衡时气相的压力，此时气相的摩尔自由能与液相的摩尔自由能完全相等。饱和蒸气压是流体最重要的基础热力学参数之一，石油馏分的沸点、闪蒸计算、气液平衡判断都需要依托饱和蒸气压数据，常规条件下纯烃的饱和蒸气压随温度升高呈指数上升趋势，临界温度下饱和蒸气压等于临界压力。7.压缩因子：无量纲流体状态参数，定义为Z=PV_m/(RT)，对于理想气体Z恒等于1，对于实际流体Z的数值反映其偏离理想气体模型的程度：Z>1说明实际流体的可压缩性低于理想气体，分子间的排斥力占主导；Z<1说明实际流体可压缩性高于理想气体，分子间的吸引力占主导。压缩因子是油气田地面工程中天然气流量计量、混输管路水力计算的核心修正参数，高压凝析气条件下Z的数值可低至0.7以下，直接决定计量误差的大小。二、基础理论简答题1.化工热力学引入余函数（偏差函数）的工程必要性理想气体模型是热力学性质计算的基准体系，理想气体状态下焓、熵等热力学性质仅为温度的函数，可通过定压摩尔热容的积分直接计算，方法简单、数据体系完善。但在石油化工实际生产场景中，绝大多数工艺单元的操作压力范围为1~20MPa，部分加氢装置操作压力可达30MPa以上，流体分子间交互作用强烈，若直接采用理想气体模型计算焓、熵值，偏差最高可达30%以上，完全无法满足工艺设计要求。如果直接通过实验逐一测定不同温度压力下所有纯组分、混合物的热力学性质，实验工作量极大、不同组分混合物的组合数量近乎无限，完全不具备可操作性。余函数的本质是将实际流体的热力学性质拆分为两部分：第一部分是同温同压下对应理想气体状态的热力学性质，可通过已有的完善热容数据库直接积分得到；第二部分是实际流体与理想气体的性质偏差，依托状态方程、对应状态原理即可通过少量纯组分临界参数、偏心因子计算得到。该方法既充分利用了理想气体体系已有的完善数据积累，又通过少量修正即可覆盖全操作范围的工程计算需求，目前已经成为炼化工艺模拟、油气集输工艺设计的通用标准方法。以陕北油田伴生气处理装置的丙烷制冷循环为例，采用余焓法计算节流过程的温度降，计算值与现场实测值的偏差可控制在2℃以内，远低于理想气体模型计算得到的15℃偏差。2.常见立方型状态方程的适用场景对比立方型状态方程因求解摩尔体积时为三次方程、同时能够兼顾计算精度与计算简便性，在石油化工领域得到了最广泛的应用，不同类型方程的特性和适用场景差异显著：第一类是原始RK方程，1949年由Redlich和Kwai提出，在范德华方程的基础上修正了分子引力项的温度影响因子，引力项引入T^-0.5的修正，形式为P=RT/(V-b)-a/[T^0.5V(V+b)]，其计算精度显著高于范德华方程，对于非极性纯气体的pVT性质计算偏差可控制在2%以内，但对于气液平衡、液体密度的计算精度较差，目前主要应用于中低压气体流量计量、压缩机工艺计算场景。第二类是SRK方程（Soave修正RK方程），1972年Soave将原始RK方程中的引力项常数a修正为与温度、偏心因子相关的函数，大幅提升了气液平衡常数的计算精度，针对石油馏分这类复杂多组分体系的泡点、露点计算偏差可控制在1K以内，是目前国内炼油厂精馏塔、原油稳定装置工艺设计的首选状态方程。第三类是PR方程（Peng-Robinson方程），1976年由Peng和Robinson提出，在SRK方程的基础上修正了引力项的体积影响因子，进一步提升了液体密度的计算精度，对于烃类液体密度的计算偏差小于5%，同时在近临界区域的数值稳定性远高于SRK方程，目前主要应用于凝析气藏开发、超临界萃取工艺的模拟计算。三类立方型方程的共同优势是参数少、仅需纯组分的临界参数、偏心因子即可完成所有计算，编程实现难度低，完全适配工业级实时计算场景。3.维里方程不同截断形式的适用压力范围维里方程从统计力学推导而来，不同截断阶数对应不同的压力适用区间：当维里方程仅截断至第二维里系数项，即PV_m/RT=1+B/(RT)P时，其适用压力范围为0~1.5MPa的中低压区间，该条件下三分子以上的碰撞概率低于1%，高阶维里项的贡献可以忽略，该简化形式广泛应用于常压油气回收装置、油田伴生气常压输配管路的pVT计算，计算偏差小于1%。当维里方程截断至第三维里系数项，即PV_m/RT=1+B/V_m+C/V_m²时，适用压力范围提升至0~5MPa的中压区间，该条件下三分子碰撞效应不可忽略，高阶项的贡献占比约3%~5%，截断至三阶即可满足精度要求，该形式常用于城市天然气门站、中压输气干线的工艺计算。当操作压力超过5MPa、最高至临界压力附近时，四阶及以上维里项的贡献逐步提升，需要引入更多阶数的维里系数才能保证精度，此时维里方程的参数获取难度大幅提升，工程上通常改用立方型状态方程或者普遍化压缩因子图法完成计算，避免高阶维里系数实验测定成本过高的问题。4.稳定流动能量方程在炼化不同工艺单元的简化规则稳定流动过程的通用能量方程为ΔH+0.5Δu²+gΔz=Q-W_s，其中ΔH为单位质量流体的焓变、0.5Δu²为单位质量流体的动势能变化、gΔz为单位质量流体的位能变化、Q为单位质量流体和外界交换的热量、W_s为单位质量流体对外输出的轴功。针对石油炼化的不同工艺单元，可结合实际工况做合理简化：第一类是管输、节流类工艺单元，这类单元没有对外输出轴功W_s=0，同时管内流体流速通常低于30m/s，动势能变化仅为450J/kg，相对于烃类流体焓变的万焦级数值占比不足1%，位能变化gΔz即使针对100m高的精馏塔也仅为980J/kg，同样可以忽略，因此能量方程可简化为ΔH=Q，绝热条件下进一步简化为ΔH=0，即绝热节流过程焓值恒定。第二类是压缩机、透平膨胀机类动力单元，这类单元和外界的热量交换Q相对于轴功的数值占比不足2%，工程上通常按绝热处理，同时动势能、位能变化同样可以忽略，能量方程简化为ΔH=-W_s，即轴功的数值直接等于流体的焓变，可直接通过进出口的焓差计算压缩机功耗、透平输出功。第三类是换热器类换热单元，这类单元既没有轴功交换，流体的动势能、位能变化也可忽略，能量方程直接简化为ΔH=Q，换热器的换热量等于冷热流体的焓差乘以质量流量，是目前炼化装置换热网络计算的核心依据。三、工程应用计算题1.长庆油田某气井采出的天然气摩尔组分为甲烷0.92、乙烷0.05、丙烷0.03，现场工况条件为温度318.15K、压力12MPa，分别采用普遍化压缩因子图法、RK方程计算该混合气体的摩尔体积，对比现场实测值0.000182m³/mol计算偏差。首先采用凯氏规则计算混合物的虚拟临界参数：查得甲烷临界温度190.56K、临界压力4.599MPa，乙烷临界温度305.32K、临界压力4.872MPa，丙烷临界温度369.83K、临界压力4.248MPa，混合物虚拟临界温度T_pc=Σy_iT_ci=0.92×190.56+0.05×305.32+0.03×369.83=203.21K，虚拟临界压力P_pc=Σy_iP_ci=0.92×4.599+0.05×4.872+0.03×4.248=4.61MPa，由此得到对比温度Tr=T/T_pc=318.15/203.21=1.566，对比压力Pr=P/P_pc=12/4.61=2.603，查皮策普遍化压缩因子图得到该条件下压缩因子Z=0.835，根据压缩因子定义计算摩尔体积V_m=ZRT/P=0.835×8.314×318.15/12000000=0.000184m³/mol，与实测值的相对偏差为1.1%，满足工程计量的精度要求。后续采用RK方程计算该混合物的摩尔体积，首先计算各组分的RK方程参数：纯组分i的a_i=0.42748×R²T_ci^(2.5)/P_ci，b_i=0.08664×RT_ci/P_ci，混合物交叉项a_ij=(a_ia_j)^0.5，通过混合规则a=ΣΣy_iy_ja_ij、b=Σy_ib_i计算得到混合物的引力项参数a=3.247Pa·m⁶·K^0.5/mol²，体积修正项参数b=2.987×10^-5m³/mol，代入RK方程采用迭代法求解，初始迭代值取理想气体摩尔体积8.314×318.15/12000000=0.000220m³/mol，经过3次迭代后得到收敛结果V_m=0.000181m³/mol，与现场实测值的相对偏差仅为0.55%，精度高于普遍化压缩因子图法，完全适配高压气井的流量计量场景。2.某丙烷制冷循环单元，丙烷初始状态为温度360K、压力3MPa，经过绝热节流阀降压至0.2MPa进入蒸发器制冷，采用余焓法计算节流后的丙烷温度以及节流过程的摩尔熵变。首先查得丙烷的临界参数T_c=369.83K、P_c=4.248MPa、偏心因子ω=0.152，初始状态下对比温度Tr1=360/369.83=0.973，对比压力Pr1=3/4.248=0.706，采用普遍化余焓公式(H-H^ig)/RT_c=(H^0-H^ig)/RT_c+ω×(H^1-H^ig)/RT_c，查普遍化余焓数据表得到简单流体余焓修正项为-1.68，非简单流体余焓修正项为-1.38，代入计算得到初始状态余焓H1-H1^ig=8.314×369.83×(-1.68+0.152×(-1.38))≈-5820J/mol=-5.82kJ/mol。以273.15K下液态丙烷饱和焓值为基准0，查理想气体丙烷定压热容积分得到360K下理想气体丙烷的摩尔焓H1^ig=32.78kJ/mol，因此初始状态实际摩尔焓H1=32.78-5.82=26.96kJ/mol。绝热节流过程焓值恒定H1=H2，终态压力0.2MPa远低于丙烷的临界压力，近似处于理想气体状态，仅需要通过理想气体焓值与温度的对应关系试差求解，最终得到节流后丙烷温度为293.7K，节流温降达到66.3K，与现场实测的65℃温降几乎完全吻合。熵变计算过程中初始状态余熵S1-S1^ig=-11.2J/(mol·K)，终态低压下余熵近似为0，理想气体状态下从360K、3MPa膨胀到293.7K、0.2MPa的理想气体熵变为30.1J/(mol·K)，最终得到总熵变ΔS=18.9J/(mol·K)，符合绝热不可逆过程熵增大于0的热力学基本定律。四、油田现场案例分析题陕北某油田伴生气处理站2024年建成投用，设计处理能力为15×10⁴m³/d，投产初期夏季环境温度35℃、压缩机入口操作压力