Aspen高级班石油炼制过程的模拟解读.ppt

2020年5月23日,AspenPlus高级班石油炼制过程的模拟,2020年5月23日,高级精馏在化工领域的应用,RADFRAC单元操作模型的应用、收敛方法和参数的选取；三相蒸馏、萃取蒸馏、共沸蒸馏和反应蒸馏等特殊蒸馏过程的模拟；电解质过程的模拟等。,2020年5月23日,石油炼制过程的模拟方法,以一套典型的常减压装置为主线，讲授下列5部分内容：1)油品馏程分析数据的处理；2)初馏塔的模拟计算；3)常压塔的模拟计算；4)减压塔的模拟计算；5)计算原油中所含杂质硫在各种产品中的分布。,2020年5月23日,掌握如下新的内容：,1)石油炼制过程蒸馏单元的模拟步骤；2)流程模拟软件中石油馏份的特殊处理方法；3)掌握PETROFRAC单元操作模型的使用；4)产品中杂质分布的计算。,2020年5月23日,复习：建立流程模拟模型的基本步骤,确定模拟范围，绘制模拟流程图；将已知条件输入相应表格：SETUPCOMPONENTSPROPERTIESSTREAMSBLOKS。运行模拟计算检查结果,油品分析和虚拟组分,目标:学习输入石油分析数据,创建调和组分(油品混合),以及表征石油馏分AspenPlus参考资料:UserGuide,Chapter6,SpecifyingComponentsUserGuide,Chapter32,PetroleumAssaysandPseudocomponents,2020年5月23日,本章概况,表征石油流体的挑战油品馏程数据和性质数据的输入蒸馏曲线的转换生成虚拟组分油品调和用户自定义的虚拟组分,2020年5月23日,石油的表征,石油是烃类化合物和杂质组成的宽沸程的混合物.石油或石油馏分的性质取决于它们的组成，但是其组成难以定量分析石油用虚拟组分(pseudocomponents)来表征,虚拟组分是蒸馏曲线上较小的片断.每一个虚拟组分的基本性质由下列关联式来计算:分子量(MW)正常沸点(TBP)重度或者WatsonUOPK,2020年5月23日,Assays和Blends组分,AssayData(TBP,D86),Pseudocomponents,Assay-1Assay-2,Blend,Pseudocomponents,输入输出,2020年5月23日,输入油品馏程分析数据,1.在DataBrowser下面,选择Components,然后选Assay-Blend.2.在Assay-Blend目标管理器下,创建一个assay并且提供一个名称(ID).3.在BasicDataDistCurve表中，输入蒸馏曲线数据.4.在BasicDataDistCurve表中，输入油品重度，或者在BasicDataGravity/UOPK表中，输入重度曲线或UOPK曲线数据.虚拟组分将在模拟计算之初自动生成.,2020年5月23日,ASSAY组分的表征,输入蒸馏曲线蒸馏曲线分割为很多个片段(TBPCutpoints)每一个片段生成一个虚拟组分每一个虚拟组分的性质以下列参数的平均值为基础进行估算：分子量(MW)正常沸点(NBP)密度或比重(APIorSG)虚拟组分的流率在流股结果表中报告.,2020年5月23日,Assay组分BasicDataDistCurve表,2020年5月23日,重度的规定,油品整体重度比重API重度重度曲线比重API重度API=141.5/SG-131.5WatsonUOPK曲线NBP=(UOPK*SG)3,2020年5月23日,Assay组分BasicDataGravity/UOPK表,2020年5月23日,蒸馏曲线,蒸馏曲线类型:ASTMD86,D1160,D2887实沸点蒸馏曲线(TBP)(重量或体积基准)真空蒸馏曲线(重量或体积基准)每一个蒸馏曲线至少需要输入4个点.报告中的蒸馏曲线为:ASTMD86ASTMD1160Vacuum(liquidvolume),2020年5月23日,油品分析数据计算程序Version9或更新的版本Version8或更老的版本蒸馏曲线转换方法ASTMD86至实沸点(Trueboilingpoint)EdmisterEdmister-OkamotaPMLD2887至TrueboilingpointAPI87API94TBPWT,蒸馏曲线计算方法,2020年5月23日,曲线处理选项初馏点(Initialboilingpoint)0.5%(缺省)终馏点(Finalboilingpoint)99%(default)外延方法Probability(缺省)Quadratic锯齿段拟合方法Harwell(缺省)Hermite曲线处理选项在ComponentsPetroCharacterizationAnalysisOptionsAssayProcedure表上.,蒸馏曲线计算方法(续),2020年5月23日,虚拟组分从油品分析数据生成.你可以规定切割温度或切割片段(窄馏份)的宽度.ComponentsPetroCharacterizationGenerationCuts表上，规定切割温度等.缺省值是40个切割点:,切割温度,2020年5月23日,蒸馏曲线的切割,2020年5月23日,虚拟组分性质,用Watson/UOP“K”和API关联式估算：典型地10UOPK13此处10.0表示芳烃含量很高13.0表示烷烃含量很高,2020年5月23日,所生成虚拟组分的命名,ComponentsPetroCharacterizationGeneration规定命名选项自动生成的名称可以是:实际计算出来的NBP(如:PC864F)沸程下限温度(如:PC800F)沸程上限温度(如:PC900F)用户选择的顺序号(如:PC1,PC2)用户定义的名称,2020年5月23日,附加的ASSAY组分数据,可选的性质曲线:轻端组分分析(light-endsmatching)比重或者API重度曲线WatsonUOPK曲线分子量曲线性质曲线(多于30种性质)绝对&动力粘度曲线如果输入油品整体性质，单独的曲线数值被归一到整体性质.可以应用用户性质曲线,2020年5月23日,Assays组分的用户性质曲线,在表格PropertiesAdvancedUserProperties上，建立用户性质曲线.在表格中激活AssayCurveproperties选项在表格ComponentAssay/BlendPropertyCurve上，输入用户性质曲线的数据。对于高度非线性调和现象，如马达法辛烷值，建议使用用户调和子程序。,2020年5月23日,性质曲线的流股计算,性质曲线被切割为与虚拟组分同样多的片段。为每一个片段计算出平均的性质数值流股的性质数值由虚拟组分以及它们在流股中的浓度做调和计算得出计算虚拟组分混合物的性质，有各种调和计算选项。,2020年5月23日,石油性质的调和,石油性质调和方法:标准体积(StdVol)摩尔(Mole)质量(Mass)用户(User)用于性质计算的用户调和子程序：用户子程序名称为BLDPPU.“调和选项”可以规定.性质可以外延与否.在ComponentsPetroCharacterizationAnalysisOptionsBlendOptions上，调和规则可以修改。,2020年5月23日,Assay(石油)数据库,软件内置的Assay数据库，包含194种石油的分析数据,其中有:-来自Phillips石油数据库的10种石油-来自文献的184种石油在ComponentSpecificationsPetroleum表，查找Assay库文件。和PhillipsPetroleum石油分析数据库的接口，其中收集了超过250种石油的数据*需要与PhillipsPetroleum签订License协议。,2020年5月23日,Assay(石油)数据库,每一种石油至少提供如下信息:馏程分析重度和/或重度曲线某些石油有附加的性质曲线:苯胺点闪点凝固点硫含量金属含量硫醇研究法/马达法辛烷值,烷烃环烷烃芳烃雷诺蒸汽压,2020年5月23日,访问Assay数据库,2020年5月23日,创建Blends组分,通过规定blend中每个assay的百分数，由assay创建blend.调和可以是Mole,Mass,或StdVol基准.要创建一个blend，至少规定2个assay组分.一个blend中的assay的数量没有上限.报告中的蒸馏曲线为:ASTMD86ASTMD1160Vacuum(liquidvolume),2020年5月23日,输入Blend组分,1.在数据浏览器DataBrowser,选择组分Components,然后选择Assay-Blend.2.IntheAssay-Blend目标管理器,创建一个blend并且提供名称(ID).3.在MixtureSpecifications表上，输入各assay组分的调和分数.虚拟组分将在模拟计算之初自动生成.,2020年5月23日,Blend组分之MixtureSpecifications表,2020年5月23日,Assay/Blend结果,Assay/blend的计算结果可以在表格ComponentsAssay-BlendResults查看.每一个assay或blend组分有自己单独的结果表格(resultsform).,2020年5月23日,GenerationCuts表格,2020年5月23日,虚拟组分的生成,缺省情况,一套完整的虚拟组分由输入的所有assay组分生成.在表格ComponentsPetroCharacterizationGenerationSpecifications上，用户可以规定生成虚拟组分时，包括/不包括特别的assay.只有在表格GenerationSpecifications上选出的assays和blends，才用于生成虚拟组分。在表格GenerationSpecifications上未选用的Assaysorblends，将不被使用，它们的IDs不在流股输入表格上显示.,2020年5月23日,虚拟组分命名,虚拟组分的命名约定计算出的沸点(例如：PC125,PC136,)规定起始编号切割温度上限(例如：将225-250F窄馏份命名为PC250)切割温度下限(例如：将225-250F窄馏份命名为PC225)用户定义的名称列表在表格ComponentsPetroCharacterizationGenerationNamingOptions上，规定命名约定.如果你要将AspenPlus生成的虚拟组分与反应器或者分离模块关联，应使用用户定义的名称列表。,2020年5月23日,GenerationNamingOptions表格,2020年5月23日,1.在表格ComponentsSpecificationsSelection上，输入组分IDs并且选择Pseudocomponent类型。2.在表格ComponentsPseudocomponentsSpecifications上,至少提供下列参数中的两个：正常沸点NBP,重度Gravity(Specificgravity,APIgravity,orLiquiddensity),分子量MolecularweightAspenPlus可以估算第3性质.3.可选地,输入与温度相关的性质数据:液相粘度Liquidviscosity水溶解度Watersolubility饱和蒸汽压Vaporpressure,定义用户虚拟组分,2020年5月23日,虚拟组分规定表格,2020年5月23日,虚拟组分物性方法,方法名称,描述,API-METH,COAL-LIQ,LK,API-TWU,APIprocedure,Correlationsforcoal-liquidfractions,LeeKeslercorrelations,APIprocedurewithTwucorrelationsforcriticalproperties,ASPEN(default),APIprocedurewithAspenmodifications,注意:不要将虚拟组分物性方法与流程模拟物性方法混淆，前者只是用于虚拟组分石油性质参数的计算。,2020年5月23日,虚拟组分计算结果,生成的虚拟组分在表格ComponentsPetroCharacterizationResults上浏览检查。,2020年5月23日,表格AnalysisOptionsAssayProcedure,2020年5月23日,石油性质模型,基本的虚拟组分性质模型分子量MolecularWeight临界温度CriticalTemperature临界压力CriticalPressure临界体积CriticalVolume偏心因子AccentricFactor饱和蒸汽压VaporPressure液相摩尔体积LiquidMolarVolume水溶解度WaterSolubility粘度Viscosity基本的虚拟组分性质模型可以在表格ComponentsPetroCharacterizationPropertiesMethodsBasic上修改。,2020年5月23日,石油性质模型(续),虚拟组分热力学性质模型理想气体热容IdealGasHeatCapacity汽化焓EnthalpyofVaporization标准生成焓StandardEnthalpyofFormation标准生成自由能StandardFreeEnergyofFormation基本的虚拟组分性质模型可以在表格ComponentsPetroCharacterizationPropertiesMethodsThermodynamic上修改.,2020年5月23日,石油性质模型(续),虚拟组分状态方程(EOS)性质模型RKSInteractionParametersBWROrientationParameterBWRCriticalTemperatureBWRCriticalVolume基本的虚拟组分性质模型可以在表格ComponentsPetroCharacterizationPropertiesMethodsEOS上修改.,2020年5月23日,石油交互作用参数,虚拟组分与下列常规组分的交互作用参数Kijs:H2SN2CO2COH2CH4,2020年5月23日,物性方法的修改,用户可以建立自己的虚拟组分石油性质方法.各种性质的计算模型可以改变.大多数虚拟组分性质可以应用子程序.在线帮助中有每一种性质方法的简要说明,2020年5月23日,流程模拟计算,在流程模拟计算期间,不区分虚拟组分和真实组分.产品流股中的虚拟组分调和在一起，并且生成蒸馏曲线产品终馏点(即D8695%point)或者初馏点(D865%point)可以规定为单元操作的性能指标：塔内部的设计规定用于其它单元操作的外部设计规定,2020年5月23日,检查蒸馏塔各产品之间的分离效果(重叠度),PetroleumPropertiesandPropertySets,Objective:BecomefamiliarwithpetroleumpropertiesandlearntheuseofstreamanalysisAspenPlusReferences:UserGuide,Chapter32,PetroleumAssaysandPseudocomponents,2020年5月23日,PetroleumProperties,ReidvaporpressureAPIandspecificgravityWatsonUOPKfactorAnilinepointPourpointFlashpointCetanenumberCarbontohydrogenratioGrossheatingvalue,2020年5月23日,PetroleumProperties(Continued),MeanandvolumeaverageboilingpointHydrateformationpressureandtemperatureRefractiveindexLiquidviscosityindexD86temperaturesat(0,5,10,95%)D1160temperaturesat(0,5,10,95%)D2887temperaturesat(0,5,10,95%)TBPtemperaturesVacuumtemperatures,2020年5月23日,PetroleumProperties(Continued),Properties(curves)fromAssayAnalysisPointsAniline,cloud,flash,freeze,smokeContentsAromatic,olefin,paraffin,naphthalene,mercaptanC,H,O,NImpurities:S,V,Fe,NiIndices:Antiknock,vaporknock,refractive,warm-upNumbers:MotorandresearchoctaneLuminometerViscosityGenerallyinputcurvesforassaystogeneratethesecurvesforproductstreamsisneeded.,2020年5月23日,DistillationCurves,VolumeandweightbasisASTMD86ASTMD1160ASTMD2887(weightpercentonly)VacuumTBPAPIGravityGravityMW,2020年5月23日,PropertySets,Apropertyset(Prop-Set)isawayofaccessingacollection,orset,ofpropertiesasanobjectwithauser-givenname.Onlythenameofthepropertysetisreferencedwhenusingthepropertiesinanapplication.Usepropertysetstoreportthermodynamic,transport,andotherpropertyvalues.Currentpropertysetapplicationsinclude:Designspecifications,Fortranblocks,sensitivityStreamreportsPhysicalpropertytables(PropertyAnalysis)Trayproperties(RadFrac,PetroFrac,etc.)Heating/coolingcurves(Flash2,MHeatX,etc.),2020年5月23日,PropertiesincludedinProp-Sets,Propertiescommonlyincludedinpropertysetsinclude:API-APIgravitySGSTD-Standardspecificgravity(60F,1atm)REIDVP-ReidvaporpressureD86T-D86temperatureMUMX-ViscosityAvailablepropertiesinclude:ThermodynamicpropertiesofcomponentsinamixturePurecomponentthermodynamicpropertiesTransportpropertiesElectrolytepropertiesPetroleum-relatedpropertiesReference:AspenPlusPhysicalPropertyDataReferenceManual,Chapter4,PropertySetshasacompletelistofpropertiesthatcanbeincludedinapropertyset.,2020年5月23日,DefiningPropertySets,Specificationsforapropertysetinclude:PropertynameandunitsTemperatureandpressureatwhichtocalculateIncludeorexcludewaterphasequalifier(V,L,L1,L2,T,S),2020年5月23日,SpecifyingPropertySets,UsethePropertiesProp-Setsformtospecifypropertiesinapropertyset.TheSearchbuttoncanbeusedtosearchforaproperty.Allspecifiedqualifiersapplytoeachpropertyspecified,whereapplicable.UserscandefinenewpropertiesonthePropertiesAdvancedUser-PropertiesformbyprovidingaFortransubroutine.,2020年5月23日,PredefinedPropertySets,SomesimulationTemplatescontainpredefinedpropertysets.ThefollowingtablelistspredefinedpropertysetsandthetypesofpropertiestheycontainfortheGeneralTemplate:,2020年5月23日,PredefinedPetroleumProp-Sets(Continued),PredefinedProp-SetswiththePetroleumApplicationtemplate:,PREDEFINED,PROPERTYSET,TYPESOFPROPERTIES,CUTS-EVolumetricflowrateofpetroleumcutsat100FintervalsD86-5ASTMD86temperatureat5%liquidvolumedistilledD86-95ASTMD86temperatureat95%liquidvolumedistilledGAS-PROPSCompressibility,volumetricflowrate,heatcapacityratioKINVISCLiquidkinematicviscosityat100Fand212FLIGHTPetroleumcharacteristicsforlightdistillatesMIDDLEPetroleumcharacteristicsformiddledistillatesPETROVolumetricflowrate,gravity,WatsonK,distillationcurvesTBP-5Trueboilingpointtemp.at5%liquidvolumedistilledTBP-95Trueboilingpointtemp.at95%liquidvolumedistilled,2020年5月23日,PetroleumProp-Sets,Mostpetroleumpropertiesarecalculatedfrominputpropertycurves.Somepetroleumpropertiesarecalculatedfromcorrelations(FLSH-API,POUR-API)Propertiesareaccessedandusedinexactlythesamewayasotherprop-setproperties,2020年5月23日,StreamResultsOptions,OntheSetupReportOptionsStreamsheet,use:FlowBasisandFractionBasischeck-boxestospecifyhowstreamcompositionisreportedPropertySetsbuttontospecifynamesofpropertysetscontainingadditionalpropertiestobereportedforeachstream,2020年5月23日,Prop-setsDiscussion,1.Givetworeasonsthatapropertysetisneeded.2.Howdoyoudefineapropertyset?3.Auserdefinedaprop-setforviscosityandthenranthesimulation.Afterthesimulationwascomplete,heopenedthestreamresults.Tohissurprise,hedidnotseeviscosity.Whathappened?,PetroFrac,Objective:IntroducethefeaturesofthePetroFracmodelanditsuseinsimulatingacrudetowerAspenPlusReferences:UnitOperationModelsReferenceManual,Chapter4,Columns,2020年5月23日,CharacteristicsofRefiningProcess,WideboilingrangefluidsLargetemperaturegradientsalongcolumnMostfluidcompositionsareunknownbuttheirdistillationanalysisisknownThreephasewithfreewaterSidestrippersandpumparoundsPetroleumspecificphysicalpropertiesHighdegreeofheatintegrationMetalandsulfurcontents,2020年5月23日,PetroleumRefiningApplications,PetroFracisadistillationmodelspecificallydesignedforpetroleumrefiningapplications:preflashtowercrudeunitvacuumunitvacuumlubefractionatorFCCmainfractionatordelayedcokermainfractionatorquenchtower/primaryfractionator,2020年5月23日,PetroFracFeatures,Integrationoffurnace,column,sidestripperandpumparoundADA/PCSforassayinput,crudeblendandbreakdownConversionbetweendifferentASTMcurvesUserspecificationofcutsPetroleumorientedspecificationsTrayandpackingsizingandratingOtherfeaturesInitializationandconvergence,2020年5月23日,ExampleCrudeTower,GAS,DECANTWATER,NAPHTHA,JETFUEL,KEROSENE,DIESEL,GASOIL,RESID,CRUDEFEED,STEAM,2020年5月23日,PetroFrac-Connectivity,VaporDistillate,Bottoms,Feeds,Feeds,Overflash,Pumparound,Pumparound,StripperBottoms,LiqReturn,SideProducts,Overhead,LiquidDistillate,WaterDistillate,Feed,Heat,BottomStageorReboilerHeatDuty,TopStageorCond.HeatDuty,HeatorReboilerHeatDuty,SteamFeed,2020年5月23日,PetroFrac-FlowsheetConnectivity,AnynumberofpumparoundsandsidestrippersSteamorreboilerstrippersStripperhotliquidreturnMultiplefeedstostrippersTotalpumparoundandproductdrawoffsInterstagewaterdecantersHeatstreamsforpumparoundcoolersPseudostreams,2020年5月23日,PetroFrac-ModelSpecifications,Setupform(required)NumberofstagesCondenserandReboilerDistillaterate,Bottomsrate,Condenserduty,Reboilerduty,orRefluxRatioFeed/Product(required)Feedconvention:furnace,on-stage,above-stagePres-Profile(required)PumparoundsStrippers,2020年5月23日,PetroFrac-PumparoundSpecifications,Drawandreturnstagelocations(required)Fortotaldrawoffs,specifyoneofthefollowing:returntemperaturetemperaturechange(negativeifcooled)duty(negativeifcooled)Forpartialdrawoffs,specifyliquidflowrateandoneoftheabove.,2020年5月23日,PetroFrac-SideStripperSpecifications,StageNumberForeachstripper,thetraysarenumberedfromtoptobottomstartingwithstage1.SteamStripperSpecifyingasteam-to-productratiowilloverridethesteamflowrate.ReboilerStripperSpecifyareboilerduty.,2020年5月23日,PetroFrac-SideStripper,ConventionsLiquiddrawalwaysgoestothestrippertopstage.VaporalwaysreturnstothemaincolumnOtherSpecificationsAfractionofthestripperbottomscanbereturnedtothemaincolumnMultipledrawoffsfromdifferentstagesasfeedstothesamestripperareallowed,2020年5月23日,PetroFrac-OtherSpecifications,Steamtoproductflowratio(mass)FeedoverflashFurnaceoutlettemperatureLiquidrunbackfromastageASTMtempsandgapsPetroleumproperties(RVP,flashpoint,pourpoint,etc)MurphreeefficiencyMultiplepropertyoptionspecificationsFreewatercalcsinthecondenseroronalltraysBlocksPetroFracPropertiesformtospecifyinternalstagescheckedforfreewater-condenserischeckedautomaticallyTray/packing,2020年5月23日,PetroFrac-DesignSpecifications,Meetproductpropertybyvarying:strippingsteamflowrate/reboilerdutyproductflowratepumparoundflowrateorreturntemperaturefurnacetemperatureMeetaparticularstagetemperaturebyvarying:furnaceoutlettemperaturepumparoundflowrateorreturntemperatureMeetproductflowratebyvarying:furnaceoutlettemperature,2020年5月23日,PetroFrac-Definitions,RunbackLiquidflowratefromastagetoalowerstage,equivalenttostageliquidflowratelesssidedrawsTheRunbackformisusedtospecifytherunbackfromastage.ThemanipulatedvariableisthepumparounddutyorflowrateLiquidFlowInPetroFrac,stageliquidflowreportedinthetrayprofileisthesumofrunback,sidedraw,andpumparound.InRadFrac,stageliquidflowexcludesanysideproductdrawandpumparound.,2020年5月23日,OutsideLoopPhysicalProperties(ai,KB,HV,HL),InsideLoop,cmi=cmioutyi=KixiSxi=1.0,Syi=1.0,SHiin=SHiout,l,v,x,y,L=Sli,V=Svi,T,InsideloopsolvedbyNewtonmethodwithorwithoutlog-leg,PetroFrac-Convergence,Sum-Ratesvariantoftheinside-outalgorithm,2020年5月23日,PetroFrac-Convergence,Algorithm:sum-ratesvariantoftheinside-outInitialization:designedforrefiningapplicationToenhanceconvergence:providecolumn/strippertemperatureestimatesvaryconvergencedampingfactorstoenhancestabilitychangedesignspecifications,2020年5月23日,PetroFrac-Summary,1.PetroFraccombinesFurnace,Column,Pumparound,andSidestripperintoonemodelwhichmakesthesimulationeasytosetupandconverge.2.Feedcharacteristicscanbeuserinputorfromtheassaylibrary.3.Multiplefeedscanbeblendedtocreateafeedmixture.4.Additionalpropertiesofaparticularproductcanbereportedbyaddingpropertysets.5.Runbackfromastageistheliquidflowrateleavingthestageandenteringthestagebelow.Runbackandliquidflowcanbedifferent.6.Mostcrudetowersdonothaveareboileratthebottom.Eitherthebottomsteamrateorthefurnacedutyisvariedtomeettheenergyrequirement(energybalance).,2020年5月23日,PetroFrac-Discussion,1.RadFracandMultiFraccanalsobeu