毕业论文输油管线参数优化设计软件的编制

输油管线参数优化设计软件的编制摘要：管道优化设计是一项复杂的系统工程，目标是寻求既满足工艺要求，同时经济效益较好的方案，但工作量和计算量庞大，通过编制计算机软件既提高工了作效率，节省时间，而且能更直观更准确地对管线参数进行操作，从而达到优化的效果。论文对优化设计的数学模型类型、寻优方法及发展趋势，以及油气管道工程的优化技术研究进行了全面综述，根据流体力学及传热学相关理论建立了热含蜡原油管道系统优化设计的数学模型,并编制了输油管线参数优化软件。本软件采用VisualBasic6.0编制，以MicrosoftOfficeAccess为数据库，通过ADO数据控件绑定链接，在窗体上采用Dategrid控件显示数据，具有网格式大量数据读入、温度压力计算及分布曲线的绘制、停输时间计算、推荐运行参数和数据输出等几大功能。通过对白豹到华池的输油管线的优化设计的实例计算,表明计算结果与管线实际运行参数相符，优化结果令人满意，证实了本软件的可靠性与实用性。该优化问题的解决对于提高输油管系统的设计、管理水平具有重要的实际意义。关键词：输油管道；含蜡原油；优化设计；加热输送；软件编制ThedesignofsoftwareforoptimezingtheparameterofpipelineAbstract:Pipelineoptimizationdesignisacomplexsystemengineering.Itaimsatfindingtheprogramthateconomicallymeetsthetechnologyrequirement.However,thesystemengineeringneedsenormouscalculationsandworks.Theuseofsoftwarecannotonlyhighlyraisetheefficiency,butalsocalculatetheparameteraccuratelyinintuitionisticway,consequentlyachievetheoptimizationofpipeline.Thispaperfirstgivesaoverallsummaryofmathematicalmodelsofoptimaldesign,optimum-seekingtechniquesandtheirdevelopmenttrendsandtheoptimizationresearchesinpipelineproject.Accordingtothetheoryofhydromechanicsandheattransfor,theauthorbuildthemathematicalmodelfortheoptimaldesignofpipelinesystemtransportinghotwaxycrudes,thendesignthesoftwareforoptimezationofpipelinebaseonthemathematicalmodelinvisualbasic6.0language.ThesoftwareuseMicrosoftOfficeAccessasdatebase,whichwaslinkedbyADOcontrol,thendateswereshowedintheformbydategridcontrol.Themainfunctionsofthesoftwarearedateinputing,calculatingtemperatureandpressureanddrawingtheirgraphs,calculatingtimeofhaltandoptimalparameter.ThecalculationsofthepipelinewhichisfromBaibaotoHuacheindicatethattheresultsaccordwithpracticalparameter,consequentlyachievetheoptimezationofpipeline.Thereby,thereliabilityandpracticabilityofthesoftwareisconfirmed.Itisespeciallysignificantforimprovingthedesignandmanagementartsofpetroleumpipelinetothesolutionoftheproblem.KeyWord:Oilpipeline;Waxycrudeoil;Optimizationdesign;Hottransportation;software目录1前言 11.1问题的提出 11.2输油管道优化运行研究及其软件编制的发展状况11.2.1国外输油管道优化运行及其软件编制研究的发展状况 21.2.2国内输油管道优化运行研究的发展状况 31.3管道优化及其软件编制的意义 51.3.1节约能源 51.3.2降低油品的生产成本 51.3.3有利于环境保护 51.4主要研究内容62基本算法的确定 62.1热油管道沿程温降计算 72.2热油管道的摩阻计算82.3保证管线正常输运所允许的停输时间 82.4原油输送能耗的计算123软件的设计153.1软件的设计环境及功能说明153.2软件的结构图163.3软件的操作界面及使用说明173.3.1启动界面及主操作界面173.3.2数据输入193.3.3温度及压力计算213.3.4停输时间233.3.5推荐运行参数244计算实例254.1管道温度及压力分布254.2停输时间的计算274.2推荐运行参数295结论及建议31毕业设计小结32参考文献33致谢34附录1软件源代码351前言1.1问题的提出石油管道运输是近三四十年发展起来的一种运输方式，特别是近二十年来发展尤为迅速。随着国民经济发展和石油工业的壮大，石油管道运输已经成为自己的运输体系，并与铁路、公路、水路、航空运输一起成为五大运输体系。目前的管道运输已不仅仅运送石油，它已向多品种管道运输发展，已经使用和正在研究使用管道运输的物品有汽体、成品油、水、粮食、煤桨、矿浆及LPG等不仅陆地广泛使用管道运输，海底油气管道也得到大量运用，可以说管道运输已在国民经济发展中起着越来越重要的作用。在现代管道工业中，优化技术在管道规划、设计、运行管理、控制等方面得到了广泛的应用。优化技术是在现代计算机技术广泛应用的基础上发展起来的一项新技术，是根据现代数学最优化原理和方法，综合各方面因素，在现有的工程条件下，从问题的众多可行方案中选出最经济的方案[1]。输油管道优化运行就是在管线系统本身的结构、流程及其外部环境条件为给定的前提下，合理制定管线的输油计划及运行方案，使之在规定期内完成给定输油任务的同时运行费用最低[2]。输油局是输油企业，也是能源消耗大户，在原油输送中每年就要消耗26.5万吨燃料油，12.4亿度动力电。因此，如何才能实现管道的安全运行、低耗、节能、提高系统的经济效益，这是科研攻关和运行管理的大目标[3]。输油管道优化运行解决的也是这一问题，在确保完成输油任务的前提下尽可能降低输油能耗，寻求最经济的输送条件。输油管道在稳态(包括准稳态)工况下，给定管道流量及其他不可控因素，在满足管道运行的约束条件下，要求选择管道的最优运行方式、设备组合以及运行参数使得单位时间内的能耗费用最低[4]。对正在运行中的输油管道，通过降低热能损失来降低总能耗对降低输油成本、提高经济效益影响重大，输油管道优化运行在现实生活中具有重大的意义。但是，输油管道优化运行的研究是一项浩大的工程，需要耗费大量的人力物力，在科技飞速发展的今天，计算机技术日益发达，通过编制软件不仅可以大大节省时间和资源，而且更能准确直观的实现对管线参数的操作。目前，国内尚无成熟的输油管线优化软件可供使用，所以该软件的编制具有降低输油成本、提高经济效益等重大意义。1.2输油管道优化运行研究及其软件编制的发展状况最优化技术是研究和解决最优化问题的一门学科，属数学规划的范畴，第二次世界大战以前，处理最优化问题的数学方法主要是古典的微分法和变分法，二次大战中，由于军事上的需要产生了运筹学，提出了大量不能用上述经典方法解决的最优化问题，随之产生了线性规划、动态规划、决策论、对策论、图论等新的方法。随着近代科技与生产的发展，工程与技术的复杂化，使得一个决策的好坏，对经济效益有着重大的影响，这就要求人们必须寻求最优化决策，以获得最好的经济效益，但由于最优化的任务是要在一切可能的方案中寻求最佳方案，常常需要大量的计算，若仅依靠手工计算，则不仅工作量大，计算时间长，而且往往难于实现，令人欣慰的是，由于电子计算机的飞速发展，为最优化技术提供了强有力的计算工具，使其迅速发展成为一门新兴的学科。自60年代以来，最优化理论和方法不断得到丰富和发展，其应用几乎已深入到各个生产科研领域，成为促进国民经济多快好省发展的有效方法[7]。国外输油管道优化运行及其软件编制研究的发展状况[9][8]油气管道系统最优化技术国外研究较早，水平相对较高。早在50年代，苏联学者切尔尼金就提出了热油管道最优加热温度这一概念。他针对当时采用的往复泵开式的流程，流型为牛顿流的热油管线，用微分法导出了确定最优出站各参数的关系式。六十年代初，美国壳牌石油公司的杰斐逊将动态规划法用于等温密闭输油管道的站间压力分配，随后又提出了隐枚举法确定站内最优泵组合。随后，由于计算机技术和应用数学的飞速发展，国外油气管道最优化技术迅猛发展，涉及油气管道系统设计、运行、扩建、控制等领域。尤其在最优运行与最优设计领域硕果累累，且达到了相当的深度。如Cheesernan用坐标轮换法优化管道。程序有两条主线:站位和压缩比优化。站位程序的核心是压降计算。费用程序包括两个主要参数:动力费和材料投资。假设对每一个搜索的参数只有一个最小费用。从不同起点出发看有无相同局部优化值，如果有就认为它是全局优化值。虽然这对于有约束的非线性问题并不是一个高效的算法，仍使设计时间减少70%，设计费用减少30%，设计质量得到提高。Grelli和Gilmou:设计的管线优化程序用于太平洋输气公司和阿尔伯达天然气公司的一条输气管线。程序采用稳态分析和动态规划。1985年1月和8月进行了试验，节约燃料8%-12%。在国外油气动力系统中，有不少是采用热机(如蒸汽机，燃气轮机等)作为原动力的。为使油气管道中的泵和压缩机运转，每年用于燃料的费用是很可观的。因此，不少管道公司开始对压缩机站和泵站的燃料消耗进行了最优化研究。1980年，美国的德克萨斯气体运输公司的Copal提出了一种管道泵站优化运行的方法。目标是根据每一原动机的燃料消耗率，决定要开动哪些泵机组，在保证所需流量的前提下使总燃料费用最少;用整数规划算法来选择最优的泵机组搭配，仅需对很少几种组合进行计算。此外，还用动态规划的方法对出站压力进行优化。Singh(美，1981年)则对此问题提出了一种简便的算法。苏联学者在这方面的研究主要集中在如何选择油管的最佳工况上。国内输油管道优化运行研究的发展状况对输油管道优化运行的研究是输油工作者研究的重点，建国以来，我国的石油工业和管道事业有了很大的发展，在管道优化运行研究方面所做的工作也很多，取得了较令人满意的成果。(1)优化的数学模型[10][4][8]把一个复杂的工程系统抽象成数学模型时，设计者常常面临多种可供选择的数学模型，恰当的选择数学模型是一个至关重要的问题。较理想的数学模型应尽可能准确的反映工程实际问题，同时又便于处理。已投产运行的输油管道，其设备、管径、管长、站间距等己确定，我们只能对运行参数进行优化，如温度、泵组合等。这样影响输油能耗的只有热、动力费用，一般都以二者费用之和最小来建立模型。最初建立的是线性稳态模型，后来，石油大学的严大凡、吴长春考虑了油品的温度沿管线真实分布及热机和泵效率的影响，建立了包含全线泵组合与管路的最优匹配和输油温度的最优化两级子问题的动态规划数学模型，该模型适合于我国大多数输油管道。1998年，张增强等建立了秦京输油管道优化运行的动态规划数学模型，把全线总压头和全线总热负荷在各站间进行了最优分配，提出了一种二维动态规划求解方法，在求解时用逐次逼近法降维，这种方法既充分考虑了输油管道的特点又兼顾了数学上的严密性。(2)优化方法用于优化的方法很多。线性规划问题:单纯形法可以用来解决所有线性规划问题，线性离散优化问题可以使用割平面法、分枝定界法及0-1规划的布尔运算法等。非线性规划问题:计算方法有解析法、搜索法、数值法、线性规划或二次规划逐步逼近非线性规划的方法等。解析法亦称间接法，包括古典微分法、古典变分法、极值原理和库恩一塔克定理。搜索法亦称直接法，包括消去法(格点法、黄金分割法、Fibonacci法)、爬山法(变量轮换法、复合形法)、随机试验法等。数值法包括牛顿法、变尺度法、可行方向法、鞍点迭代法、梯度法、梯度投影法、简约梯度法、制约函数法等，线性规划或二次规划逐步逼近非线性规划的方法有序列线性规划法(SLP)和(移动限制)近似规划法(MAP)等。动态规划问题:动态规划法是一种能很好的解决所谓多阶段决策过程最优化问题的方法。输油管道优化运行研究所要解决的是如何确定各站最佳的进出站油温，各种加热、加压设备的最优组合，可以十分有效的利用这种方法来解决。尽管优化的方法很多，而应用到输油管道优化运行上的方法还不是很多。在我国输油管道优化运行方面使用过的方法有:1991年，严大凡、吴长春使用了动态规划法:1992年，吴长春又提出了所谓的DFS(深度优先搜索法)，提高了搜索过程效率;1994年，吴长春又提出了先以随机搜索产Powell法的初始点，采用改进0.618法进行一维搜索的方法:1996年，抚顺石油学院的吴明使用了二分法，相继牛顿-拉普森法也得到了应用；1997年，张其敏使用了MDOD(混合离散变量直接搜索法)，该方法优化效果好，速度快:1991年，吴立峰考虑到数学模型中既有连续变量(温度)，又有离散变f(泵组合)，将MDCP法(混合离散变A组合型法)应用到了输油管道的优化运行上，该方法优化速度快、优化结果显著;后来有人采用几种功能不同的算法进行寻优，以期获得最佳运行方案。(3)泵组合和温度优化方面1980年，吴长春等人对热含蜡原油管道的经济运行方案的求解方法进行了探讨，郑国华等人也讨论过低输量管道的安全经济运行，金德馨则应用动态规化算法阐述了输油管道最佳运行方案的选择方法训。1986年起，东北输油管理局与华东石油学院对铁秦线输油优化运行进行了协作研究。同期，华东石油管理局也开始了对鲁宁线(山东临邑至南京仪征)低耗节能(即开式流程优化运行)的工业试验研究。由于铁秦线和督宁线自投产后一直处于不满负荷状态下运行，而且实际输量波动较大，因而在输量安排、输油温度选择、管线与泵组合的合理匹配等方面，均为优化运行提供了较宽的决策面，一但实现优化运行，其效益就相当可观。这两个研究项目都获得了可喜的成果。1991年，西南石油学院吴立峰的硕士论文对热含蜡原油管道进行了优化设计，采用的方法是混合离散变量复合形法，有一定的研究价值。1992年，蒲家宁运用动态规划原理，针对密闭输送工艺，简明、扼要、具体的阐述了长输管道优化运行的分析与计算方法，有助于对输油管道优化运行的研究。1992年，金岗、张心等对热油管道优化运行技术进行了探讨，指出输油管道的能耗及运行费用除受油泵压力和油品温度变化影响外，还必须考虑动力和热力设备在不同负荷下系统效率变化的影响。通过计算机系统的反馈验证，找出了偏差并制定了能控制最小运行费用的新优化运行方案。1993年，重庆大学的王勇勤等将线性规划法用于油田输油泵站原油外输机组的优化运行，减少了原油外输中油泵工况搭配和搭配机组运行时间的盲目性，有利于优化运行。1995年，王鉴光根据泵的机械特性，提出了泵站节能运行方案和泵速最优节能控制算法，节约了电能，带来了经济效益。目前辽河油田考虑到运行温度决定了能耗费用，而运行费用中可变的参数就是能耗费用，所以对沱一鞍输油管道的运行温度进行了优化计算，将中间加热站改成了热泵站，同时，也提高了沱--鞍输油管道输油能力，设计输量为500X104吨，优化前输量仅有350X104吨，优化后输量达到7450X104吨。通过实践证明，优化了运行温度后，每年可节省能耗费用30.12万元，如果管道较长，输量较大，经济效果会更加明显。1998年，吴明等用最优化方法计算了热油输送的经济温度，对节能有重要意义。2001年，魏安河等以庆铁线新庙和牧羊泵站为例，提出通过运用降低进站温度的方法，对降低油电消耗总费用的可行性进行了具体分析，通过计算证明在保证低输量的安全运行的前提下，降低进站温度的方法对减少节流损失，降低输油成本是完全可行的。我国虽然在输油管道优化运行方面的研究起步较晚，但其进步很快，储运科技人员在不断的朝着这个方向努力，建立了许多数学经济模型，寻求最优运行乡数和最优运行方案，为节能降耗做出了重大贡献。目前在这方面的研究还不是很完善，建立的模型都比较简单，求解的方法也不是很先进，还有待于我们后来之人继续研究和探讨。(4)软件编制方面1999年，清华大学的张信荣等以东辛含硫管道为例，对孤岛原油和清河原油相混后的流变性进行了研究，提出了一种计算东辛管道输送混合原油粘度的新模型。将管道运行参数优选和节能技术方案的研究相结合，建立了节能优化的改造模型，开发了东辛管道优化运行软件，既考虑了管道运行能耗费用又考虑了投资改造费用，实现了优化与改造方案相结合，该模型既可用于寻求最优运行参数又可指导最优改造方案。对我们的优化运行研究很有帮助。2004年，长江大学石油工程学院张顶学等受长庆油田分公司第二输油处委托，编制了集输管线优化软件。该软件界面友好、操作方便，功能全面，但其不足之处是在停输时间的计算上还不够准确，需要进一步完善。1.3管道优化及其软件编制的意义1.3.1节约能源节约能源及提高经济效益是当前国家对工交部门的两项重要要求。石油企业既是能源生产企业，又是能源消耗大户，仅输油的能源消费就占生产成本的50%,能源消耗费用是企业生产成本中重要的可控部分，降低单位产品的能源消耗是企业能否有效的控制生产成本的重要因素。能源的开发和合理利用是社会发展的源泉和战略依据，并标志和决定着一个国家的竞争实力和综合实力[1]。例如管径为720m，长435km的铁岭至大连的管道，年流量为200万吨，如六个热泵站每站都将出站油温提高1摄氏度，全年将多烧原油8500吨。如果每个泵站节流l00kPa，则全年将损失380万度电。若能将输油成本降低0.1分/吨.公里，则一年将节约870万元，可见节省下来的运行费用是惊人的[5]。1.3.2降低油品的生产成本[1]随着国内市场的不断规范以及同国际市场的逐步接轨，石油企业在激烈的市场竞争中通过国家政策和地区保护而获得较好的经济效益的机会越来越少，只有按照市场的需求对运行的管道进行优化，降低产品的单位成本，才能提高企业在市场经济中的竞争能力。1.3.3有利于环境保护[1]大量燃料不能合理和有效的利用，增加了污染物的排放，加重了对环境的污染。对输油管道的优化运行，不但达到了节能降耗的目的，而且可以降低企业治理工业污染的费用，对企业自身也有着直接的经济效益。由于管道建设的大型化和现代化，需要巨额的投资，如在运行中引进优化技术，对流量、温度、泵机组特性等参数进行很小的优化改变，便可节省出巨额资金。根据技术经济学原理，通过对各种优化措施和可行方案的技术分析、经济比较和效益评价、寻优技术与经济的最优组合，从而确定出经济合理的运行参数，使系统工作在最佳状态下，达到节约能源、少投入、多产出、提高经济效益的目的。可见，输油管道优化运行方案研究是一个很有潜力的发展方向，对输油管道优化运行进行研究具有重大的现实意义[6]。综上所述，油气管道系统最优化技术研究是向复杂化、多目标方向发展。随着计算机技术向高速、多功能、智能化方向发展，油气管道系统最优化技术研究的领域更广，设计、运行、控制、模拟技术会相互融合形成综合最优化技术。1.4主要研究内容本论文主要研究集输管线参数优化软件的编制。管道优化设计是一项复杂的系统工程，目标是寻求既满足工艺要求，同时经济效益较好的方案，但工作量和计算量庞大，通过编制计算机软件既提高工了作效率，节省时间，而且能更直观更准确地对管线参数进行操作，从而达到优化的效果。本论文首先对优化设计的数学模型类型、寻优方法及发展趋势，以及油气管道工程的优化技术研究进行了全面综述，建立了热含蜡原油管道系统优化设计的数学模型,然后编制了输油管线参数优化软件。本软件采用VisualBasic6.0编制，以MicrosoftOfficeAccess为数据库，通过ADO数据控件绑定链接，在窗体上采用Dategrid控件显示数据，具有网格式大量数据读入、温度压力计算及分布曲线的绘制、停输时间计算、推荐运行参数和数据输出等几大功能。最后，对白豹到华池的输油管线的优化设计进行实例计算，验证计算结果与管线实际运行参数是否相符，从而证实本软件的可靠性与实用性。2基本算法的确定2.1热油管道沿程温降计算易凝、高粘的原油当其凝固点高于管道周围环境温度，或在环境温度下油流的粘度很高时，工程上常采用加热、加添加剂的办法输送，以保证原油输送能安全、经济地运行。热油在输送过程中存在着两个方面的能量损失，即摩阻损失和散热损失。因此，也必须从两个方面给油流供给能量，即由加热站供给热能，由泵站该给压能。在热油管道的设计和管理中，要正确处理这两种能量的供求平衡关系。这两种能量损失的多少又是互相影响的，其中散热损失往往是起取得决定作用的因素。因为摩阻损失的大小决定于原油的粘度，而粘度的大小则决定于输送温度的高低。提高油温，使原油在较高温度下输送，散热损失将增大，但摩阻损失则可减少。故对于某一输送任务，存在着能量最小消耗的最优输送条件，这也就是热油输送管道工艺计算所要确定的目标。为此，首先要确定热油管道沿程的温度分布。油流在加热站加热到一定的温度后进入热油管道。油流在沿管道流动的过程中不断向周围散热，使其温度下降。散热量和原油温度的沿程分布受到诸如输油量、加热温度、环境条件、管道散热条件等许多因素的影响。严格地讲，这些因素是随时间变化的，故热油管道经常处于热力不稳定状态。但在实际工程中，将正常运行工况近似为热力、水力稳定状态，在此前提下进行轴向温降的计算。设管道周围介质的温度为T0，dl微元管段上的油温为T，管道的质量输油量为G，管段的水力坡降为i。流经dl段后散热油流产生的温降为dT。在稳定工况下，dl微元管段上的能量平衡式为：上式中左边为dl管段单位时间内向周围介质的散热量，右边的第一项为管内油流温降dT的散热量；右边第二项为dl段上油流摩阻损失转化的热量。由于dl与dT的方向相反，故取负号。假定在dl管段内总传热系数K为常数，且忽略水力坡降沿dl的变化，对上式分离变量积分，得到沿程温度分布公式，即列宾宗公式：记，；则有：(2—1)(2—2)或(2—3)上式中，G——原油的质量流量，kg/s；c——输油平均温度下原油的比热容，J/(kg.℃)；D——管道外径，m；L——管道加热输送的长度，m；K——管道的总传热系数，W/（m2.℃）；TR——管道起点油温，℃；TL——距起点为L处的油温，℃；T0——周围介质温度，埋地管道则为管中心埋深处的自然地温，℃；I——油流的水力坡降，m/m；g——重力加速度，m/s2。若油流出加热站的温度TR为定值，则管道内原油沿程的温度分布可表示为：(2—4)在推导上式时，假定水力坡降i为定值。实际上，水力坡降要由管道的摩阻计算得出，在进行温度分布计算时，它还未精确给出。在本研究里我们采用迭代算法确定。从温度分布的算式可以看出，两加热站间的温度梯度在各处是不同的。在加热站的出口处，油流与环境温度的温差大，温降快；而在进入下一加热站前，由于油温低，温降慢。故加热温度愈高，散热愈多，温降就愈快。因此，试图靠提高加热站出口的油温来提高管段末端的油温，往往收效不大。常常是首端温度提高10℃，而末端的油温仅提高2～3℃。2.2热油管道的摩阻计算由于热油管道油流的温度沿程是不断下降的，故油流的粘度不断增大，其水力坡降也是不断增大的。因此，在计算热油管道的摩阻时，必须考虑管道沿线的温降情况和油品的粘温特性。在本研究里，我们采用分段耦合计算热油管道的压力和温度分布。在本研究中，我们采用平均油温的计算法计算油流的压力沿程分布。其计算步骤为：a．按甲方的要求，将两加热站间的热油管道分划成需要小的分段，在计算程序里，我们大致按0.1千米分划，记下每个计算点的管道里程和管道中心的高程；b．由反推计算确定的总传热系数K、原油的比热容和假定的管段内的水力坡度，从加热站开始，逐段计算该段的终点温度T2；c．按下式计算该段原油的平均温度：(2－5)d．由实测的原油粘温曲线，内插得到平均温度下的原油粘度，计算油流的雷诺数，并判别流态；e.由科尔布鲁克公式计算油流的沿程阻力系数，进而求得该段的水力摩阻损失和计算的水力坡度；f．比较计算的水力坡度和原假定的水力坡度的差别是否精度要求，若不满足，则返回（b），重新计算；g．若计算的水力坡度和原假定的水力坡度的差别在给定的误差范围内，则由伯努利方程求的该管段末端的压力P2；h．将计算得到的末端温度和压力作为下一段的起点温度和压力，重复上述计算，直到下一加热站为止。2.3保证管线正常输运所允许的停输时间本节所讨论的问题实质是要弄清楚热油管线的停输时间和再启动压力的计算问题。热油管路的计划检修、电源中断、输油量过低、输油温度过低、管道泄漏都可能造成热油管线的停输。在现场，通常在输量不足的情况下，为节省动力并为了保持一定的热力条件，也常采用间歇输油的方式来代替反输。停输后，由于管路中所输油品的温度不断下降，油的粘度不断增加，给管路输油的再启动造成困难，甚至造成油品在管路中凝结，使管路“冻结”。为了避免冻结事故，确保管线安全平稳生产，必须了解管路在各种条件下的允许最长停输时间和停输后的温度变化情况，确定出再启动时所需的压力和输量，以确保管线的安全生产。（1）埋地热油管线停输后温降的计算埋地热油管线停输后，由于经长期稳定输油，管线周围土壤中所积蓄的热量要比管路中原油所含热量大很多。故埋地热油管线停输后的温度下降情况与架空管道不同，主要决定于管线周围土壤的冷却情况。实践表明，由于刚停输时管内壁处的油温与管外壁土壤间的温差较大。在地温为12～15℃，油温为45～50℃时，大约相差3～8℃。当地温为2℃，油温为35～40℃时，约相差8～10℃。故埋地管线停输后地温降可分为两个阶段：●管内油温较快地冷却到略高于管外壁土温，其中靠近管壁处的油温下降更快；●管内存油和管外土壤作为一个整体而缓慢地冷却。ⅰ)开始阶段的温降计算对于这一阶段的温降计算，我们作如下假设：①忽略管内横断面上温度梯度和轴向温降；②管道总传热系数K保持不变；③忽略油品特性随温度的变化；④外界温度不变；⑤不计保温层的热容量。在时间内的热平衡方程为：(2—6)将上式积分，并整理后得到：(2—7)式中，K—停输后油管内原油至油管外壁处的总传热系数，W/(m2·℃)；D,D1,D2—管道平均直径、油管内径和油管外径，m；—原油的比热和密度，J/(㎏·℃),㎏/m3；—钢材的比热和密度，J/(㎏·℃),㎏/m3；T0———管道埋深处的地温，℃；———停输后小时的油温，℃；———停输开始时的油温，℃；———停输时间，h(小时)。若考虑到钢管的壁厚不大，可以不计钢管的热容量，则上式可以更进一步化简为：(2—8)ⅱ)管内存油和管外土壤的共同缓慢冷却若输油管道正常运行时管道周围的土壤温度已接近稳定，则停输后土壤温度场的衰减过程，也就是启动预热的反过程。可按恒热流法计算停输不同时间后的各点土壤温度，得到在冷却过程中管壁处的土壤温度随时间的变化规律。我们选用恒热流法计算停输不同时间后的各点土壤温度。这种方法是把埋地管道当作半无限大均匀介质中连续作用的线热源，即认为土壤是各向均匀同性的，管道传往各方向的热流强度是相等的。即：(2—9)式中， —开始停输时管壁处的土壤温度； —停输小时后管壁处的土壤温度；——土壤的导热系数，W/(m·℃)；——土壤的导温系数，；——管道埋深，m；——停输时间，小时，h；——稳定工况的管道散热量，W/m；若将上式也可以(2—10)如果近似取停输后的恒定散热损失等于稳定输送时的热损失，则(2—11)代入（3－5）式，并考虑当(3～4)时，<<,忽略前一项，并将展开为级数，取前两项，整理后得到第二阶段的停输时间：(2—12)由（3－10）式，可以求得温降第一阶段的停输时间为：式中，以秒计。上述恒热流法系一种二维不稳定传热计算。它只能求得某一断面或某一管段处的停输温降。下面为一种计算管道任意断面上停输后油温的近似计算方法。（2）短期停输后管内油温计算由于土壤的蓄热量和热阻都很大，埋地管道在停输后，温度的下降是比较缓慢的。因此，可以将管内原油的冷却过程看成是一系列准稳定状态，列出微元时段内的热平衡方程，进行近似求解，得出停输后管内的油温。假定在时段内，管内所存原油及钢管温降的放热量等于管道向环境的散热，且总传热系数对于稳定状态时的总传热系数K，类似可导得：(2—13)令，(2—14)则上式可写成：对上式积分，当。此时，距起点l处，油温Tl可按苏霍夫公式计算：(2—15)于是，可得距起点任意距离l处，停输温降时间后管内油温的计算式为：(2—16)式中，T——距管道起点l处，停输时间后的油温，℃；T0——周围环境温度，℃；TR——管道起点处，开始停输时的油温，℃；K——管道在稳定状态工况下的总传热系数，W/(m2.℃)；L——计算点距管道起点的距离，m；b——由（3—14）式定义的系数；——停输时间，s。（3—16）式可以计算停输后管道任意断面上平均油温随时间的变化。各种突然事故都可能造成输油管线停输。因此，停输可能发生在不同的季节（不同的地温）、不同的排量、不同的起点温度、不同的起点压力下，停输的时间长短也不一。在停输温降计算的软件设计中，应考虑上述各种可能的情况的停输事件的发生，及时和准确地计算出停输管线沿线各断面处的油温随时间变化的情况。限于篇幅，在本论文中的各加热站的停输计算中，只对冬季最寒冷的月份发生停输事件的温降进行了计算。其他的各种可能的停输温降计算可以很方便地利用本研究提供的软件进行。2.4原油输送能耗的计算原油在加热输送的过程中，能量消耗主要为提高原油位能和克服管道摩阻损失的机械能损失和加热输送的管道散热损失。a）机械能消耗在无翻越点的原油输送管道中，机械能消耗仅为提高原油位能和克服原油输送时的摩阻损失。泵提供给单位重量原油的有效能头为：（2-17）而在输送管道存在翻越点的管线中，若不允许管道出现负压，则泵提供给单位重量原油的有效能头为：（2-18）上两式中，泵的有效功率为：（2-19）上式中，若以单位能耗，即输送1吨原油每1千米所消耗的能量表示，再计入泵和电机的效率，则有：（2-20）式中：b）热能消耗单位热能的消耗以下述公式计算：（2-21）式中，3软件的设计3.1软件的设计环境及功能说明软件的开发环境采用WindowsXP操作系统，因为Windows是一个多任务的操作系统，它是集多任务、多用户、多功能和多媒体为一体的软件集成环境。本软件采用的开发工具为VisualBasic6.0，它有以下特点:VB继承了Basic简单易用的特点，适用于Windows环境下快速编程；采用可视化技术,操作方便；采用面向对象技术,没有复杂的程序流程；编程模块化、事件化、程序由许多小程序组成。本软件采用MicrosoftOfficeAccess数据库，以ADO数据控件链接，在窗体上用Dategrid控件显示数据。该软件由数据输入、模拟计算、结果输出等三大模块组成。(1)数据输入模块该模块包含基本管线基础参数输入和管线高程参数输入及原油粘温数据三个子模块。(2)模拟计算模块该模块由管道温度、压力的分布计算模块、保证任务数量下可将能耗将为最低的优化运行参数计算模块及停输降温计算模块组合而成，用户可根据屏幕提示和计算要求进行选择，以便节省计算时间。(3)结果输出模块该模块包含文件保存、数据显示、粘温曲线、温度、压力分布曲线等子模块。3.2软件的结构图集输管线参数优化软件启动界面集输管线参数优化软件启动界面主控菜单输入模块管道基础数据管道高程数据计算模块温度计算压力计算推荐参数计算停输时间计算输出模块粘温曲线温度、压力曲线原油粘温数据管道纵断面图3.3软件的操作界面及使用说明3.3.1启动界面及主操作界面本软件具有友好的界面设计，在主操作界面上，用菜单编辑器制作了主菜单和其子菜单，为了用户能够方便的使用该软件，使用工具栏（ToolBar）控件，在主操作界面中制作了工具按钮，然后链接Imagelist控件为按钮添加了彩色图标。主菜单包括有管道及原油数据，温度压力计算，停输时间，推荐运行参数，帮助及退出选项。在管道及原油数据中，包括有子菜单管道基础数据、管道高程数据和粘温数据。工具栏中的工具按钮包括有管道基础数据、管道高程数据、粘温数据、温度压力计算、停输时间、推荐运行参数及退出。在状态栏中，有“欢迎使用输油管线参数优化软件”字样和当前的时间显示。其源代码见附录1.1。3.3.2数据输入本软件采用MicrosoftOfficeAccess数据库，以ADO数据控件链接，在窗体上用Dategrid控件显示数据，从而实现了网格式的大量数据的输入。（1）管线基础数据（2）管线高程数据以白豹至华池的输油管线为例，管线全长27.072km，以0.1km为步长，共277个计算点。在界面左侧为Dategrid控件，可以输入不同管线的各点高程数据，右侧为管线纵断面图。在图上方，可实时显示鼠标位置处的座标，从而可以直观的查询到管道任意位置处的高程。此功能是通过Picturebox(图片框)控件的line方法实现的，以管线里程为横坐标，管线高程为纵坐标，两点连线，再加以循环即可。具体源代码见附录1.3。（3）粘温数据此窗口用来显示原油的粘温数据，并绘制原油的粘温曲线。以原油温度为横坐标，原油粘度为纵坐标，具体实现方法与管道高程数据窗口相同。实现此功能的源代码见附录1.4。3.3.3温度及压力计算此窗口功能为在给定的管线起点温度与起点压力条件下，来计算管线沿程的温度与压力的分布，绘制温度及压力分布曲线，并将计算结果保存到文本文件中。通过Textbox（文本框）来实现管道起点温度与管道起点压力的输入，根据输入的数据，来计算管道沿线的温度与压力的分布。具体源代码见附录1.5。3.3.4停输时间此窗口用来实现管道停输时，在已知地温的条件下，根据输入的停输时间来计算管道沿线的温度分布，并将计算结果保存到文本文件，从而可以判断管道停输多久可能发生原油凝管事故。源代码见1.6。3.3.5推荐运行参数此窗口用来实现对管道运行参数的优化。根据不同的输送温度下，管道的运行费用不同，以费用最低时的输送温度为最优温度，并计算此起点温度下，管线沿程温度和压力的分布，绘制温度、压力的分布曲线，然后将优化结果保存到文本文件。具体源代码见附录1.7。4计算实例以白豹至华池的输油管线为例：管线全长27.072km，起点温度为32℃，起点压力为3.0MPa，以0.1km为步长，进行计算。由于计算点较多，这里仅列出整数km的计算点数据。4.1管道温度及压力分布由软件计算得出的温度压力数据如下表：表4-1管道温度压力数据表里程（km）高程（m）温度（℃）压力（MPa）01358323.01.001136331.302.962.001136730.632.933.094137629.912.854.029139129.322.735.087139928.662.666.049154428.091.477.069153327.501.568.094158426.921.149.084156626.381.2910.094164425.840.6511.044167225.350.4212.00156524.871.3013.056146524.362.1214.028140223.902.6415.028137823.442.8416.128136322.962.9617.034136122.572.9718.045135122.143.0619.134134321.703.1220.034133221.353.2121.034132720.973.2522.034131220.603.3823.034130620.233.4324.034130819.893.4125.034129119.553.5526.034129019.223.5627.022129118.913.55图4-1管线纵断面图图4-2管线温度分布曲线图图4-3管线压力分布曲线图4.2停输时间的计算停输时间的计算结果如下：表4-2停输后温度分布数据表停停输时间油温里程1h5h10h20h030.6726.0621.5615.591.00130.01225.5321.1615.362.00129.3825.0220.7815.143.09428.7024.4820.3614.904.02928.1424.0320.0214.715.08727.5223.5319.6514.496.04926.9723.1019.3214.317.06926.4122.6518.9814.118.09425.8722.2118.6513.929.08425.3521.8018.3413.7510.09424.8521.4018.0313.5711.04424.3821.0317.7513.4112.023.9320.6617.4813.2513.05623.4520.2717.1913.0814.02823.0119.9316.9312.9415.02822.5819.5816.6612.7916.12822.1219.2116.3812.6317.03421.7518.9216.1612.5018.04521.3518.6015.9112.3619.13420.9318.2615.6412.2120.03420.6018.0015.4612.1021.03420.2417.7115.2511.9722.03419.8917.4315.0311.8523.03419.5517.1614.8311.7424.03419.2216.8914.6311.6225.03418.9016.6414.4411.5126.03418.5916.3914.2511.4027.02218.2916.1514.0711.30图4-4管线不同停输时间后的温度分布B--停输1h后的温度分布C--停输5h后的温度分布D--停输10h后的温度分布E--停输20h后的温度分布由上图可看出，随停输时间的延长，温度曲线趋于平缓。停输时间达到一定长后，油温将趋于地温。4.2推荐运行参数经优化后，管线运行参数如下表：表4-3推荐运行参数数据表里程（km）高程（m）温度（℃）压力（MPa）0135831.031.001136330.332.962.001136729.682.933.094137629.002.854.029139128.432.735.087139927.802.666.049154427.251.477.069153326.681.568.094158426.131.149.084156625.611.2910.094164425.100.6411.044167224.630.4212.0156524.171.3013.056146523.682.1214.028140223.242.6415.028137822.802.8316.128136322.342.9617.034136121.962.9718.045135121.563.0619.134134321.143.1220.034133220.803.2121.034132720.433.2522.034131220.083.3823.034130619.833.4324.034130819.403.4125.034129119.083.5526.034129018.763.5627.022129118.463.55图4-5管线的推荐运行温度E--优化前管线的温度分布C--优化后管线的温度分布图4-6管线的推荐运行压力5结论及建议通过编制输油管线参数优化软件和对实例的计算分析，得到以下结论：1.本论文对以往的输油管道优化运行所做的工作进行了分析和研究，考虑了运行费用的影响因素，在进站油温、出站油温、热负荷、水力约束、管线强度、流量等综合约束条件下，建立了输油管道优化运行的数学模型。2.由于计入了摩擦生热对沿线温降的影响，所以本文推出了各种流态下的水力坡降i的计算公式，并将沿线的水力坡降i取作变量，但由于i随油流温度T的变化而改变，同时计算温度T时，T又受i的影响，故使用迭代法，较精确的计算了水利坡降i与油流温度T沿管线的分布，从而大大减少了由于摩擦生热而产生的误差。3.软件具有网格式的数据输入界面，计算管线沿程温度、压力的分布，停输时间及推荐运行参数，并绘制了管线纵断面图、原油粘温曲线、温度压力分布曲线及停输时的温度分布曲线等功能。4.以白豹至华池的输油管线为例，进行了计算，与实际测量值相符，从而验证了本软件的正确性与实用性，优化后取得了令人满意的效果。几点建议：1．计算时没有考虑传热系数K的变化，将其设为定值，但随着温度的降低，必将出现结蜡现象，势必会影响到管道的传热系数，应将数学模型更加完善，分段计算传热系数K的值并且考虑结蜡的影响，使计算结果与工程实际更接近。2．只考虑了稳态的密闭型输油管道，为了尽可能反映工程实际问题，应当将非稳态、旁接油罐等输油方式考虑进去，使输油管道优化运行的数学模型进一步完善化。毕业设计小结三个多月的毕业设计已经结束，而我也将告别我的大学生活。这三个多月可以说是我在大学里最重要的三个月，这次毕业设计不但加深了我对基础知识的理解，更增强了我的实际动手能力，对我以前所学知识是一次系统的、全面的总结，对我以后的学习和工作都有着重要的意义和深远的影响。这次毕业设计是对我以前所学课程的一次全面总结。因为题目涉及数学、流体力学、传热学、专业课程及计算机等多方面的知识。因此，又是对我综合能力的考察。在毕业设计过程中，曾遇到很多困难，我都一一克服。有很多老师和同学给了我很多的帮助和启示。在这里向他们表示真诚的感谢！由于水平有限，时间仓促，不足之处在所难免，还请各位老师批评指正。参考文献[1]俞伯炎、吴照云等，石油工业节能技术。第一版.北京：石油工业出版社.2000.397-401页[2]吴长春：热油管道的优化运行.石油规划设计.1994.5(2):48-49[3]茹旅灵、闰宝东：物油管道节能技术概论.第一版.北京：石油工业出版.2000.110-118页[4]张思俊译，高粘原油的优化输送.国外油气储运，Vo1.11,No1,1993.2[5]吴长春、严大凡：热油管道运行管理最优化.江汉石油学院学报.1992.83-87[6]王为民：国内外石油管道输送技术发展综述.管道技术与设备，1997.4[7]杨筱茜、张国忠：抽油管道设计与管理.第1版.北京:石油大学出版社.1996.[8]谢英：输气管网系统研究.西南石油学院硕士论文.1997[9]吴立峰：热含蜡原油管道系统的优化设计.西南石油学院硕士论文.1991[10]黄志潜等.石油与天姗气管道工程最优化技术的综述.油气储运.1987,6(2);1-5[11]张增强等.秦京物油管道优化运行研究.油气储运.1998，17(10):9-10[12]徐严波、徐严义等，MDCP法在输油管道优化运行中的应用，油气储运，2002,21(12)19--220[13]高松竹、汪玉春等，混合遗传算法在输油管道优化运行中的应用,油气储运,2004,23(7)34～37.[14]杨筱蘅、张国忠等，《输油管道设计与管理》，石油大学出版社[15]张其敏：输气干线模糊优化设计.西南石油学院硕士论文.1997[16]侯铜瑞：贯彻管运输技术政策提高管道运输综合能力.油气储运,Vo1.18,No7,1999[17]MicrosoftCorporation著，微软（中国）有限公司译，VisualBasic6.0中文版程序员指南，1998[18]林慕新等，VisualBasic6.0实例教程，电子工业出版社，1999[19]EricBrierley、AnthonyPrintce、DavidRinaldi著，王建华、陈一飞等译，VisualBasic6开发人员指南,机械工业出版社.致谢向蒋华义老师和油气储运教研室的各位老师致以衷心的谢意！这次毕业设计是在蒋华义老师悉心指导下完成的。蒋老师渊博的知识和严谨的治学态度使我感触颇深。在毕业设计过程中，每当遇到困难，蒋老师都给予我充分的鼓励和帮助，使我能够顺利完成毕业设计。在此，对蒋老师表示由衷的感谢，并感谢油气储运教研室的各位老师，感谢他们给我们的帮助，从他们那里我学到了很多全新的知识，大大开拓了我的知识面，他们丰富的经验，孜孜不倦的教诲，给我留下了深刻的印象。在软件编制的过程中，鲁明俊老师给与了细心的指导和热心的帮助，在此予以诚挚的谢意。除此之外，很多同学都给了我极大的帮助，在这里向他们表示真诚的谢意！附录1软件源代码1启动窗体PrivateSubMain()Form7.Left=3000Form7.Top=3000Form7.ShowForm7.RefreshFori=1To300000000NextiUnloadForm7MDIForm1.ShowEndSub1.2程序主操作界面PrivateSubDatebest_Click()Form6.ShowEndSubPrivateSubDateHight_Click()Form4.ShowEndSubPrivateSubDateofoil_Click()Form5.ShowEndSubPrivateSubexit_Click()UnloadMeEndSubPrivateSubMDIForm_Load()MDIForm1.Width=18000:MDIForm1.Height=11000DimL0()AsDouble,Z0()AsDoubleDimnAsInteger,TAsIntegern=1OpenApp.path+”\管线里程数据表.txt"ForInputAs#1DoWhileNotEOF(1)ReDimPreserveL0(n)Input#1,L0(n-1)n=n+1Loopn=n-1Close#1ReDimZ0(n)T=0OpenApp.path+”\管线高程数据表.txt"ForInputAs#1DoWhileNotEOF(1)Input#1,Z0(T)T=T+1LoopClose#1EndSubPrivateSubP_Click()Form3.ShowEndSubPrivateSubT_Click()Form2.ShowEndSubPrivateSubTandP_Click()Form2.ShowEndSubPrivateSubToolbar1_ButtonClick(ByValButtonAsMSComctlLib.Button)SelectCaseButton.KeyCase"DateBasic"Form1.ShowCase"DateHight"Form4.ShowCase"Dateoil"Form5.ShowCase"TandP"Form2.ShowCase"Tstop"Form3.ShowCase"Datebest"Form6.ShowCase"exit"UnloadMeEndSelectEndSubPrivateSubTstop_Click()Form3.ShowEndSub1.3管道高程数据PrivateSubForm_Load()DimL0()AsDouble,Z0()AsDoubleDimnAsInteger,mAsInteger’读入数据n=1OpenApp.Path+"\管线里程数据表.txt"ForInputAs#1DoWhileNotEOF(1)ReDimPreserveL0(n)Input#1,L0(n-1)n=n+1Loopn=n-1Close#1ReDimZ0(n)m=0OpenApp.Path+"\管线高程数据表.txt"ForInputAs#1DoWhileNotEOF(1)Input#1,Z0(m)m=m+1LoopClose#1’绘制管道纵断面图Picture1.Scale(0,Picture1.Height)-(L0(n-1),0)Picture1.DrawWidth=1.5m=0DoWhile(m<n-1)Picture1.Line(L0(m),Z0(m)*15-18100)-(L0(m+1),Z0(m+1)*15-18100),vbBluem=m+1LoopEndSub’显示里程、高程坐标PrivateSubPicture1_MouseMove(ButtonAsInteger,ShiftAsInteger,XAsSingle,YAsSingle)Label3.Caption=Str((Y+18100)/15):Label2.Caption=Str(X)EndSub1.4粘温数据PrivateSubForm_Load()DimT0()AsDouble,N0()AsDouble,aAsDoubleDimnAsInteger,mAsInteger’读入数据n=1OpenApp.Path+"\粘温数据表.txt"ForInputAs#1DoWhileNotEOF(1)ReDimPreserveT0(n),N0(n)Input#1,T0(n-1),N0(n-1)'Picture1.PrintT0(n-1),N0(n-1)n=n+1Loopn=n-1Close#1’绘制粘温曲线Picture1.Scale(0,Picture1.Height)-(T0(0)+5,0)Picture1.DrawWidth=1.5m=n-2DoWhile(m>=0)Picture1.Line(T0(m),N0(m)*70+400)-(T0(m+1),N0(m+1)*70+400),vbRedm=m-1LoopEndSub’显示温度、粘度坐标PrivateSubPicture1_MouseMove(ButtonAsInteger,ShiftAsInteger,XAsSingle,YAsSingle)Label3.Caption=Str((Y-400)/70):Label2.Caption=Str(X)EndSub1.5温度压力计算PrivateSubCommand1_Click()DimkAsDouble,dAsDouble,cAsDouble,GAsDouble,aAsDouble,mdAsDoubleDimiAsDouble,i0AsDouble,bAsDouble,lAsDouble,L0()AsDouble,Z0()AsDouble,P0()AsDoubleDimTAsDouble,T0()AsDouble,Tp()AsDouble,V0()AsDouble,Re()AsDouble,Re1AsDouble,Re2AsDouble,EAsDoubleDimnAsInteger,mAsInteger,fAsDoubleConstpi=3.1415926,G0=9.8'原始数据的读入n=1OpenApp.Path+"\管线里程数据表.txt"ForInputAs#1DoWhileNotEOF(1)ReDimPreserveL0(n)Input#1,L0(n-1)n=n+1Loopn=n-1Close#1ReDimZ0(n)m=0OpenApp.Path+"\管线高程数据表.txt"ForInputAs#1DoWhileNotEOF(1)Input#1,Z0(m)m=m+1LoopClose#1'开始进行温度分布的计算ReDimT0(m),Re(m),V0(m),Tp(m),P0(m)T0(0)=32:T=7.6:P0(0)=3*10^6G=37/3600:k=1:d=0.1583:md=838.6:i=0.001:E=2*0.0008/dc=(1.687+3.39/1000*30)/Sqr(md/1000)*1000Re1=59.5/E^(8/7):Re2=(665-765*Log(E)/Log(10))/Ef=0Forn=1Tom-1Doi0=ia=k*pi*d/G/md/c:b=G0*i/c/aT0(n)=(T+b)+(T0(n-1)-(T+b))*Exp(-a*(L0(n)-L0(n-1))*1000)Tp(n)=T0(n-1)/3+T0(n)*2/3V0(n)=3.759*10^-4*Exp(-0.126*Tp(n)):Picture1.PrintV0(n)Re(n)=4*G/pi/0.1583/V0(n)If(Re(n)<=2000)Theni=64*(4*G/pi/d^2)^2/2/9.8/dElseIf(Re(n)<=Re1)Theni=0.0246*G^1.75*V0(n)^0.25/d^4.75ElseIf(Re(n)<=Re2)Theni=0.0802*10^((0.127*Log(E/2)/Log(10))-0.627)*G^1.877*V0(n)^0.123/d^4.877Elsei=0.0826*0.11*(E/2)^0.25*G^2*V0(n)/d^5EndIff=f+1LoopWhileAbs(i-i0)>0.000001P0(n)=P0(n-1)+md*G0*(Z0(n-1)-Z0(n))-i*md*G0*(L0(n)-L0(n-1))'Picture1.Print"L=",L0(n-1),"T=",T0(n-1),"P=",P0(n-1)/1000000,fNextn'将计算结果输出为文本OpenApp.Path+"\温度分布.txt"ForOutputAs#1Forn=0Tom-1Write#1,L0(n),T0(n)NextnClose#1OpenApp.Path+"\压力分布.txt"ForOutputAs#1Forn=0Tom-1Write#1,L0(n),P0(n)/1000000NextnClose#1'绘制温度分布曲线图Picture1.Scale(0,Picture1.Height)-(L0(n-1),0)Picture1.DrawWidth=1Picture1.Line(0,0)-(0,Picture1.Height),vbRedPicture1.Line(0,20)-(Picture1.Width,20),vbReda=Picture1.Width/L0(n-1)n=0DoWhile(n<m-1)Picture1.Line(L0(n),T0(n)*100+300)-(L0(n+1),T0(n+1)*100+300),vbRedPicture1.Line(L0(n),P0(n)/900+300)-(L0(n+1),P0(n+1)/900+300),vbBluen=n+1LoopEndSubPrivateSubCommand3_Click()UnloadForm2EndSub’显示里程、温度及压力的坐标PrivateSubPicture1_MouseMove(ButtonAsInteger,ShiftAsInteger,XAsSingle,YAsSingle)Label2.Caption=Str(X):Label3.Caption=Str((Y-300)/100):Label4.Caption=Str((Y-300)*900/1000000)EndSub1.6停输时间PrivateSubCommand1_Click()DimL0()AsDouble,T0()AsDouble,TT()AsDouble,T2()AsDoubleDimnAsInteger,mAsInteger,iAsInteger,TtimeAsDoubleDimdAsDouble,timeAsInteger,kAsDouble,cAsDouble,TAsDouble,mdAsDouble,bAsDoubleConstpi=3.1415926’从文本读入数据n=1OpenApp.Path+"\温度分布.txt"ForInputAs#1DoWhileNotEOF(1)ReDimPreserveL0(n),T0(n)Input#1,L0(n-1),T0(n-1)n=n+1Loopn=n-1Close#1ReDimTT(10,n)'进行停输温降的计算Ttime=Text1k=1:d=0.1583:md=838.6:T=7.6:time=0:G=37/3600*mdc=(1.687+3.39/1000*30)/Sqr(md/1000)*1000c1=300:md1=