DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 3. 飞行员跳槽违约金适用的现状 自2004年起，为规范飞行员稳定有序的流动，民用航空总局等相继出台了《关于规范飞行人员流动管理保证民航飞行队伍稳定的意见》、《关于规范通用航空飞行人员流动管理有关问题的通知》、《民航华东地区飞行员流动管理办法》，对飞行员跳槽进行严格限定并对违约金的数额进行具体的限定，但是这些规定与《劳动法》、《劳动合同法》保护劳动者的自主择业权的基本精神相违背。由于在立法上对于违约金适用的分歧，使得飞行员跳槽的处理更为复杂。飞行员方面认为，作为一名劳动者，理应受到《劳动法》及《劳动合同法》的保护，享有单方面解除权以及自主择业权；而航空公司方面则认为，飞行员作为特殊的劳动者应适用民航局为其颁布的特殊规定，其流动应受到合理限制。 在司法实践中，飞行员劳动合同纠纷几乎都是由飞行员主动向劳动仲裁会提起，而后双方对于劳动仲裁庭的决议不服，从而向人民法院提起诉讼，人民法院对于飞行员解除劳动合同的诉讼请求都予以支持。对航空公司而言，在飞行员提出解除劳动合同时，其也相应的提出了高额的违约金对飞行员的流动进行了限制。常常见诸媒体的飞行员跳槽“天价”违约金主要包括培训费补偿金以及飞行员提前解除劳动合同的违约金。对于航空公司提出的培训费用补偿金，大多数法院对此请求予以支持；对于航空公司提出的飞行员提前解除劳动合同应予以支付违约金的请求仅有一部分法院支持，其支持率远低于对培训费用的请求。对于飞行员而言，其主要引用《劳动合同法》第31条 3 以及第38条 4 ，以期在免除高额违约金的前提下与航空公司解除劳动合同，但法院对此请求的支持率甚低。通过对劳动合同违约金适用的理论研究以及飞行员跳槽违约金的现状，将违约金适用所存在的问题归纳如下： 1) 服务期违约金的适用范围过于狭窄 我国《劳动合同法》采纳了限制违约金的立法模式，其基础为劳动者分层理论 [ 5 ] ，此理论涉及三个群体：一是强势员工，此类劳动者具备了成熟的技术才能，为市场所需的高级人才；二是具有巨大潜力成为强势劳动者的员工，此类劳动者需要通过用人单位提供的专项技术培训提升其劳动素质成为所需人才；三是弱势员工，此类劳动者仅需从事较为简单的劳动，本身不需要用人单位对其进行专项技术培训。《劳动合同法》的规定仅对第二类劳动者进行了限制，而本身较为强势，用人单位为吸引其到企业工作提供了例如购房福利、为劳动者解决户口问题、为子女出国提供便利等优惠的第一类劳动者却不在其调整之列，对其他劳动者来说有失公平。 民营航空公司为提升自身的营业能力，通常会“挖取”国有航空公司中拥有成熟的飞行技能的机长。为了吸引机长的加入，民营航空公司提供了非常丰厚的条件，而这些优惠的条件不属于《劳动合同法》规定的违约金调整之列。当飞行员再次跳槽时，违约成本仅为航空公司提供的专项培训费用，飞行员违约成本降低，跳槽越来越频繁，不利于劳动市场的稳定和谐。 2) 劳动合同服务期限缺乏合理规范 劳动合同服务期是指用人单位与劳动者约定的，在用人单位为劳动者提供专业技术培训的前提下，劳动者应当为用人单位服务的期限。用人单位在给劳动者提供专项培训服务之后才可与劳动者约定服务期，其本质为限制劳动者辞职权的适用，一旦在服务期内辞职便构成违约 [ 6 ] ，是对劳动者的劳动自由权与用人单位的营业自由权与财产权的一种平衡。《劳动合同法》仅规定了服务期的前提条件，即用人单位需提供专项培训，但没有明确服务期确定的标准，如服务期的期限长短，服务期开始和终止时间的合理限定等。许多用人单位便利用此不足，通过为劳动者提供培训服务，与劳动者签订与培训成本不相符的服务期劳动合同，部分用人单位甚至于劳动者签订了终身劳动合同，以极小的成本获得长期的服务期条款约定。无论劳动者为公司工作了几十年，一旦中途辞职，都被视为违约，需向用人单位交纳违约金。而劳动者为了提升自身的劳动素质，增强综合竞争力，在接受专项培训服务后通常都会被迫签订此类劳动合同。 在飞行员跳槽案中，用人单位为飞行员提供了数百万的高昂培训费用，并以此为条件与飞行员签订了无固定期限合同，这类合同在司法实践中通常被解释为终身服务合同，飞行员需要为公司提供服务直至退休，即使是自费参加飞行员培训的劳动者也需要与公司签订服务期长达15年的劳动合同。如果不对劳动合同中的服务期进行合理的限制，其会严重影响到劳动者的利益。 3) 服务期违约金未体现惩罚性 我国对于违约金数额的规定主要体现在《劳动合同法》第22条以及《劳动合同法实施条例》第16条。违约金的主要目的在于对违约当事人实行制裁，这不仅有助于当事人正确履行合同，而且对于减少违约行为的发生，维护交易的秩序。《合同法》实施以后，违约金的补偿性质得到加强，但在我国违约金仍然具有惩罚性。如果否认违约金的惩罚性，违约金与损害赔偿金无实质性的差异 [ 7 ] ，而如果违约金单纯具有补偿性，则违约当事人就有可能在违约造成的损失不超过预定的违约金数额时，选择违约而追求更高的利益，其相当于变相的鼓励违约，使违约金丧失保障合同履行的功能。一方面，过高的违约金会破坏劳资双方利益的平衡，侵害劳动者的劳动自由权；另一方面，过低的违约金会减少企业培养人才的动力，因此，需要以其损失为限，将违约金的惩罚性控制在合理的范围内。 在飞行员跳槽案中，根据《劳动合同法》的规定，违约金以为航空公司提供给飞行员的培训费用为限，但在司法实践中甚至出现了飞行员跳槽“零赔付”的案例 [ 8 ] 。倘若不计原航空公司的利益，片面地保护飞行员的劳动自由权，“挖人”的航空公司便能以极小的代价获得一名成熟的飞行员，如果不能有效地保护现行投入的一方，可能我们就在保护不“投”而获的一方，从而使得用人单位忌惮高层次技术人才的跳槽而减弱投入培养人才的资本，转而以廉价的方式挖去其他企业的墙角，社会将会缺乏对高级技术人才培养的动力，降低航空领域劳动力市场的运作效率。
system: [SYSTEM]

user: 4. 完善我国飞行员跳槽违约金适用的思考与建议 1) 扩大服务期违约金的适用范围 人才是当前市场竞争中制胜的关键，企业为了提升自身的核心竞争力，为吸引优秀的劳动者为其提供了优厚的福利待遇，如“安家费”、解决城市户口、优惠的住房条件等，有的企业甚至还为就业者家属提供就业机会、上学优待、大额商业保险等条件。而根据《劳动合同法》的规定，此类优惠的福利待遇是不属于违约金所调整的范围，在劳动者与用人单位解除劳动合同时，劳动者仅需赔偿用人单位专业培训的费用即可，用人单位为吸引劳动者而提供的条件则得不到法律的保护。如果不对讲服务期违约金的适用范围进行扩大，用人单位的合法利益得不到保障，其吸引优秀劳动者的投资成本减少，怠于培养优秀劳动者的热情，不利于促进整个市场经济的发展。因此，需要扩大服务期违约金的范围，将用人单位吸引优秀人才就业的特殊福利待遇也算进去，保障用人单位为培养优秀人才而投资的合法利益，稳定劳动者与用人单位之间的利益平衡。 2) 合理规范劳动合同的服务期限 用人单位在向劳动者提供了专项培训服务后，可要求劳动者履行服务期约定义务，该规定有一定的正当性，但这种正当性也应建立在合理的基础之上，一旦服务期明显不合理，那么其正当性也不复存在。服务期限是否合理的判断因素应有：① 训练期间；② 训练成本；③ 补偿措施；④ 劳动力的替代可能性。一名劳动者所接受的训练时间越长、训练成本越高、补偿措施越多、劳动力的替代性越弱，则此劳动者的服务期约定就越长。服务期的长短关乎劳动者的劳动自由权，对服务期长度加以合理的限制显得十分重要。但在我国法律规范中缺乏对服务期限的具体标准的规定，因此可借鉴其他国家对服务期的相关规定，如《德国民法典》中将服务期约定为5年 5 。对于飞行员劳动合同等用人单位需要巨大的资本投入、极长的训练时间并具有稀缺性和不可替代性的技术人才，可允许其与劳动者签订超过5年的服务期，例如可以签订10年的服务期。在合理的服务期限规定前提下，飞行员为公司工作满服务期约定后即可寻找其他的工作机会，将改变飞行员必须为一家航空公司工作到退休的局面，从而有利于规范飞行员劳动市场，缓解劳资矛盾。 3) 加大劳动合同违约金惩罚性 如前文所述，违约金是兼具补偿性与惩罚性的，而我国的立法中主要强调了其补偿性而忽略了其惩罚性，使一般辞职和服务期辞职发生冲突，不利于用人单位利益的保护，从而使其怠于培养高端人才，不利于市场经济的发展。因此，加大违约金设置的惩罚力度则显得十分必要，具体而言，就是在合同双方违反其义务时，使其承担更高的违约责任，提高违约金保障劳动合同稳定性的作用。当前劳动法合同法中所规定的违约金，不仅没有惩罚性，而且在赔偿性上也打折，违约金与实际赔偿两者从低适用。在司法实践中，为了平衡飞行员和航空公司的利益关系，法院大多采取扩大专项培训费用来增加飞行员跳槽的违约成本。此做法虽然能够从一定程度上增加违约成本，但是其不能对当事人的违约性行为进行制裁，减少违约行为的发生，帮助当事人正确履行合同。因此，适当增加对飞行员劳动合同的惩罚性，增加飞行员跳槽的违约金成本，相应地减少飞行员守约的成本，使得其他航空公司在“挖人”时能够有所忌惮，保障航空公司享有的合法利益，才能从根本上提高用人单位培养技术人才的活力。
system: [SYSTEM]

user: 摘 要 本文分析网厂协同调度业务类型(流程交互、非流程交互)，遵循OS2标准同时兼顾发电厂侧集中式访问和分布式应用需要，采用应用功能整合和调度管理流程接入两种方式开发网厂协同调度运行管理技术支撑平台，通过对业务数据的共享以及流程的交互，实现涉厂调度运行业务的统一管理，即简化用户访问复杂度又具备良好的兼容性和可扩展性，进一步加强了对网厂交互业务的统一监管，提供有力技术支撑。 关键词 :协同调度，流程对接，功能整合，OS2
system: [SYSTEM]

user: 1. 引言 目前国内多数省级电网公司在涉厂调度运行管理方面均遇到了电网同电厂业务及数据交互方式复杂、业务流程不畅，电厂运行数据上报方式不规范等问题。有的省调也曾进行了电厂接入网调系统的研究工作，例如：通过在二区建立数据上报通道实现了业务信息在网、厂之间流转 [ 1 ] - [ 3 ] 。然而这类系统架构在业务流程流转过程中仍存在一些瓶颈，存在着业务终端设备难以管控、系统可访问范围小，业务使用困难等问题 [ 4 ] 。广东省调同样进行了电厂接入网调系统的研究工作，通过二区传送调度指令，一定程度上实现了调度生产业务的自动化与信息化，并在三区建立DMZ区(数据交换区)来实现数据上报业务，然而仍未彻底解决上述的问题。
system: [SYSTEM]

user: 2. 设计思想 在广东中调与电厂之间建设协同调度运行业务服务总线，在遵循二次系统安全防护体系架构的前提下，构建跨电网公司和发电厂不同企业之间安全可靠的统一业务交互平台。 通过整合优化业务信息和流程、消除“延伸终端”业务访问模式，实现各类涉厂运行管理业务统一的功能整合、流程接入及交互方式，简化用户访问复杂度，集中维护管理等平台建设的优良特性。 根据涉厂业务类型(流程交互、非流程交互)，同时兼顾发电厂侧集中式访问和分布式应用系统建设，采用应用功能整合和调度管理流程接入两种方式开发网厂协同调度运行业务综合管理平台，通过对业务数据的共享以及流程的交互，实现涉厂调度运行业务统一管理的数据接入、信息发布、网厂在线互动及流程运转优化，进而提高对网厂交互业务需求的支撑能力 [ 5 ] - [ 7 ] 。 1) 实现网厂协同调度运行管理业务一体化管理。 2) 网厂调度运行管理业务系统应用功能的模块化建设，支持未来功能的灵活扩展，实现涉业务运行信息的“灵活互动”和功能的“即插即用”。 3) 提高涉厂业务系统对电网调度生产工作的支撑能力，最终实现保障电力生产活动安全、稳定、高效开展。
system: [SYSTEM]

user: 4. 特色与创新 1) 将OS2“横向协同，纵向贯通”的系统设计方法延伸至网厂间调度运行管理业务，为调度运行管理业务系统提供统一的业务数据交互和信息发布通道，实现网厂调度运行管理业务的一体化。 图9. 技术实现方式二 图10. 网厂协同调度运行业务综合管理平台主要功能 2) 构建网厂协同调度运行管理业务全局架构，实现网厂之间调度运行管理信息的统一交互，避免了多业务系统独立部署终端设备的竖井式系统架构。 3) 业务层面整合各个网厂间调度运行管理业务，规范业务交互模式，实现各业务数据的共享及业务流程无缝衔接；运维层面收集电厂各业务系统和设备运行情况信息，实现电厂涉网业务系统运行可靠性、安全性的综合监管。 4) 避免了生产控制区内大量使用延伸用户终端的部署方式，消除了系统安全隐患。 5) 网厂调度业务分布式系统间采用流程对接模式，实现异构系统流程对接的灵活性和装配的适应性。
system: [SYSTEM]

user: 5. 结论 本文从设计思想、总体框架、业务整合、流程接入、应用功能、关键技术、创新点等方面对广东网厂调度运行业务综合管理平台进行全面的介绍。该平台通过将各类电网与电厂间调度运行管理业务和数据进行整合，统一了广东网厂协同调度运行业务及流程的接入和交互方式。并提供全面、灵活、易访问的业务服务，同时平台定义了一套灵活的接口规范，采用表单文件式的业务交互方式，实现电网与电厂、平台与电厂其他系统的对接，为新增业务需求和流程优化提供良好的支持和扩展能力。通过对业务数据的共享以及异构系统间流程的对接规范，实现网厂协同调度运行业务的统一管理和流程的一体化管理，进而提高调度运行效率。
system: [SYSTEM]

user: 摘 要 宁夏电网考虑大用户直购电约束的中长期调度计划优化研究与应用在优化模式上兼顾了节能发电调度、经济调度、兼顾节能的经济调度等模式，在时间维度上跨越了年度、季度、月度、周、日前等多个时间维度，立足于电力系统基础理论和数学优化理论，充分考虑实际调度经验和实际运行中的实际需求，以及大用户直购电交易约束，采用先进的设计思想，并对应研究开发相应的关键性技术和模型，全面实现年度电量分解、机组组合制定、年度电量分解与机组组合相协调的中长期调度计划一体优化决策，从调度计划的制定精度、制定效率、计划结果合理性等多个方面全面提升调度计划的质量。 关键词 :宁夏电网，大用户直购电，中长期调度计划，年度电量分解
system: [SYSTEM]

user: 1. 引言 在我国的电力体制改革实施过程中，电网调度机构的地位日益凸显，电网调度运行管理已成为我国电网企业众多业务中最为核心的业务之一。受限于大规模优化计算技术的制约，长期以来，发电机组月度电量的分解与机组组合是分离决策的，没有考虑机组启停费用对电量分解的影响，没有考虑大用户直购电对机组组合计划的影响；如有问题将采用人工干预的方法影响机组月度电量的分解。这种分离决策的模式难以保证机组组合与月度电量分解相协调，难以保证电网安全与合同电量完成相配合，难以保证机组组合方案与电网安全相配合。考虑大用户直购电交易约束的中长期调度计划系统将研究突破这种传统分离计算的模式，全面实现月度电量分解、机组组合制定、月度电量分解与机组组合相协调的月度调度计划一体优化决策，从而获得高效、高精度、高可靠性的月度调度计划，确保这一决策过程能够在15分钟内计算完成。将三个独立的问题统一决策必将产生巨大的经济效益、公平效益、节能效益、安全效益，提升月度电量分解、机组组合的协调效率 [ 1 ] - [ 3 ] 。
system: [SYSTEM]

user: 2. 系统总体功能 如图1，宁夏电网考虑大用户直购电交易约束的中长期调度计划系统主要功能包括： 基础支撑模块：系统管理、基础信息管理、数据采集模块等； 负荷预测模块：数据准备、负荷预测、结果管理、计划发布等； 年度电量分解模块：计算数据准备、约束条件设置、年度电量分解、结果分析等模块； 月度机组组合模块：计算数据准备、约束条件设置、月度机组组合、结果分析等模块； 分析评估模块：统计报表及后评估模块。
system: [SYSTEM]

user: 3. 系统编制流程 如图2，考考虑大用户直购电交易约束的中长期调度计划系统编制流程如下： 图1. 中长期调度计划主要功能 首先根据历史负荷数据及节假日等信息预测未来月份的发电量、最大、最小负荷以及未来各天的最大和最小负荷、工作日和节假日的负荷曲线； 然后导入EMS系统最新的实际电网模型，导入审批后的电网设备检修计划； 根据实时电网模型和计划检修生成计划月的未来态电网模型； 年度剩余电量分解，考虑年度电量计划、已发生月份电量、大用户直购电约束等信息； 月度机组组合优化决策，考虑安全校核及校正。
system: [SYSTEM]

user: 4. 约束设置 考虑大用户直购电交易约束的中长期调度计划编制过程中考虑各种约束功能，主要包括： 1) 负荷预测：各个节点母线有功负荷是考虑电网安全、系统负荷平衡、系统备用预留等约束的基础； 2) 机组运行功能：机组煤耗参数、机组容量、出力上下限、机组爬坡能力和降出力能力、最小启停时间、机组电量约束条件、机组检修计划等主要机组参数； 3) 电厂参数功能：包括电厂的容量约束参数、电厂的电量约束参数、电厂最小开机台数、电厂运行负荷率约束等，如图3； 4) 大用户直购电：电厂年度和月度大用户直购电约束； 5) 电网拓扑参数：电网在发电计划决策周期内的拓扑信息、拓扑变化信息以及线路实际参数、发输变电设备的检修数据； 6) 电网的实际运行参数设置：线路、断面、变压器的传输容量等； 7) 机组指定状态：机组的实际运行中，部分机组的状态由于供热、检修等原因需要指定开停机状态，不参与优化求解，其他为自由优化状态，如图4； 8) 系统受限约束：包括并网减发，电网原因造成的窝电容量、电煤原因造成的缺煤停机容量、煤质差原因造成的降出力容量。
system: [SYSTEM]

user: 5. 大用户直购电与计划编制模型的整合 考虑大用户直购电交易约束的中长期调度计划是一项综合性的工作，需要考虑电网模型、电网潮流 图2. 中长期调度计划编制流程 图3. 电厂运行参数设置 图4. 机组指定状态 断面、检修申报数据、母线负荷预测、机组发电计划、大用户直购电计划、外网交换计划等诸多数据。近几年宁夏电网规模发展很快，电网结构日新月异，特别是大用户直购电的快速增长，为了保证中长期计划的正确性、科学性和合理性，需要及时更新EMS电网模型。通过对各类数据的研究我们发现，电网模型(来自EMS)数据最庞大、结构更规范、信息更全面，以电网模型为源数据，实现大用户直购电与中长期计划编制系统的规范整合，是切实可行的方案。 1) 考虑线路、主变、母线检修方式对电网拓扑的影响，对直购电电量进行安全校核。 2) 考虑大用户直购电对机组开停机方式的影响，机组最大可发电量必须大于大用户直购电量。 3) 年度电量分解结果必须满足大用户直购电电量约束。
system: [SYSTEM]

user: 6. 年度电量分解 根据所制定的年度电量计划，需要将其按照一定的原则滚动分解至月度(或季度)，以实现年度计划到月度(季度)的最佳分解。 年度电量分解至各月的目标是在电厂年度合同确定的情况下，根据月度负荷预测、机组检修，按电厂将合同电量滚动分解到每个月，保证各月分配的总电量等于该月预测电量，并尽可能保证各电厂每月的合同完成进度接近 [ 4 ] 。 年度电量分解至各月的数学模型可描述如下： 1) 目标函数： 各电厂月平均负荷率与其年平均负荷率的偏差大致相近。 2) 求解方法： 为确保年度电量分解的灵活性，适应现阶段实际应用的需求，在建立标准线性规划模型的同时，也为用户提供了启发式的年度电量分解方法。 年度电量分解的流程图如图5所示。 步骤1：计算各电厂年度实际需要分配的电量 由于月度电量预测的不断更新，各电厂年度剩余合同电量的总和并不与月度电量预测的总和相等， 所以首先计算各电厂应分配的电量( )，对超出合同外的电量(或合同的缺额)按各电厂容量大小进行分 配(也就是利用小时相等原则)。公式如下： 图5. 年度电量分解基本流程 步骤2：计算各电厂各月需要分配的电量。 优化目标是以电厂月平均负荷率与年平均负荷率偏差大致相近，故可将电厂月平均负荷率表示为年平均负荷率加上一个负荷率的偏差，相同月内各电厂的负荷率偏差相同。第j月的电量可表示为： 由于 ，可得电厂i在第j月的发电量为： 由此式计算的月度分解电量，对于各个电厂是满足年度需要完成的电量的。 步骤3：修正电厂的分解电量。 上一步骤所得的电量没有考虑最小电量和最大电量。本步骤进行修正使其满足最大和最小电量的约束，其流程如图6所示。 上述模型中，下标i表示电厂i，j表示月份； G为电厂的集合，M为月集合。 为电厂i年度应分电量； 为电厂i年度合同； 为系统年度电量预测； 为第j月预测的电量； 为电厂i年发电能力，等于电厂i的容量乘以电厂i年度非检修时间。 为电厂i在j月的发电能力，等于电厂i的容量乘以电厂i在第j月的非检修时间。 图6. 修正电厂的分解电量 为电厂i年平均负荷率； 为j月各电厂与其年平均负荷率的偏差； 为电厂i在第j月的平均负荷率。
system: [SYSTEM]

user: 7. 月度机组组合 月度机组组合是衔接年度发电计划和周、日前发电计划的关键环节，其核心是在保证未来月份电力供需平衡的基础上，优化决策火电机组次月的开停机方案，实现同类型机组之间交错开机，有序配合，在保证机组完成月度分解电量的前提下，提高机组的运行负荷率，并确保整个电网发电计划由长期至短期的有序性、安全性和平稳性 [ 5 ] - [ 7 ] 。 月度机组组合的数学模型可描述如下： 1) 目标函数： 系统发电煤耗最低、三公调度(计划跟踪)、系统购电费用最低、电力市场等多种优化目标 [ 8 ] 。 2) 求解方法： 以传统的分枝割平面算法为框架，通过对节点问题进行“松弛解优化趋整”处理，取得问题的次优解，进而快速地获得最小化问题的上界，促进分枝问题尽早尽多地裁剪，提高传统分枝割平面算法的计算速度。 月度机组组合问题最显著的计算特点有两条：一是它具有关于整数变量的时段耦合约束，二是其启停成本参数与运行成本参数具有基本相当的优化效益。这两者都导致了月度机组组合问题的计算困难。据此，可以依据下式设置惩罚向量 的各分量 ： 其中， 表示t的时段邻域，它的中心为t，半径为m； 表示在t的邻域内机组i所有时段松弛最优解分量 的平均值。从上式中可以看到， 随 的增加而减少；当 时， ，说明 对于已经在t的邻域内取1的各时段的 都不会受到惩罚；当 时， 无意义，为了避免这种 情况，应追加附加条件：当 时， ，这里 和M分别代表一极小正数与一极大正数。 这里，首先用 代替 来说明上式的“趋整”效果。在机组组合模型中，受系统备用约束的限制以及目标最小化的引导，同一时段的机组启停变量 的加和几乎为一个定值，它反映的是该时段内系统总的开机容量。从这一角度看，整个机组组合问题的实质是如何将同一时段内的系统总开机容量分配给各个机组。对于同一时段t，如果 ，则根据上式罚因子的设置，在第k + 1次迭代时， 会趋向于取更小的值，而 会趋向于取更大的值，这样就通过多次的惩罚迭代形成了一种将较大的非整数变量逼向整数1，将较小的非整数变量逼向0的驱动力，并且变量一旦取得1或者0，则其数值将不会在迭代过程再发生变化，因为其罚因子已经在同时段中取得最小或最大。这种惩罚向量的取值方式异化了问题部分变量的价值系数，打破了机组运行费用参数与启停费用参数的优化效益均衡模式，有利于问题的快速收敛。按此方式取得的整数解其优化性能也将是较好的。 然而，在月度机组组合模型中，还有最小启停时间约束这一关于启停状态变量 的时段耦合约束，其收敛效果必然会受到 时段间耦合的限制。因此，在对每个时段t的 进行惩罚时，都使用时段t邻域内 的均值，充分考虑了时段耦合特性，使其收敛效果更好。 另外“松弛解优化趋整”的迭代过程需要其他辅助方法将松弛解中的剩余非整数分量取整，可以利用邻域搜索的办法将剩余的非整数分量取整。 图7. 大用户直购电比例分析 图8. 月度机组组合结果
system: [SYSTEM]

user: 8. 结果分析评估 1) 对大用户直购电在电厂分月计划电量所占的比例进行分析，如果所占比例过高，则进行相应的手工调整，如图7。 2) 对月度机组组合结果进行分析，根据组合结果计算电厂月度运行负荷率以及系统电力电量平衡，如图8。
system: [SYSTEM]

user: 9. 结论 考虑大用户直购电交易约束的中长期计划编制系统是电网智能调度的重要创新，全面实现了中长期负荷预测、考虑大用户直购电约束的年度电量分解、月度机组组合中长期计划编制全过程，可以为宁夏电网调度工作提供科学高效的技术支持，有利于提高和改善宁夏电网调度工作的管理水平与实施效果，实现发电计划编制的最优化、协调化、规范化、流程化，并创造显著的社会效益与经济效益，实现多方共赢。
system: [SYSTEM]

user: 摘 要 为进一步提高ICEAS工艺污水处理厂的总氮去除率，本文研究了如何通过对ICEAS工艺的改良提高其脱氮效率，并最终应用于城镇污水处理厂ICEAS工艺高效脱氮升级改造。研究结果表明，连续进水方式限制了总氮去除率的提高，将连续进水方式变为间歇进水方式，系统的脱氮效率明显提高。运行模式为搅拌56 min、曝气112 min、沉淀50 min，滗水70 min，运行周期为4.8 h时，可明显提高总氮的去除率。 关键词 :ICEAS工艺，低碳氮比，实际生活污水，总氮去除率，间歇进水
system: [SYSTEM]

user: 1. 引言 滇池是云南省面积最大的高原湖泊，有着“高原明珠”之称。1999年国务院就批准实施了《昆明市城市总体规划》，用于切实保护和改善滇池的生态环境，《云南省滇池保护条例》也于2013年1月1日起正式开始施行，然而滇池的富营养化问题始终是滇池面临的最重要问题。根据2012年监测数据表明，草海和外海(滇池分为草海和外海)水质类别均为劣五类，草海综合营养状态指数为69.8，外海综合营养状态指数为68.4，营养状态均为中度富营养。 我国的污水排放标准越来越严格，采用脱氮除磷工艺的城镇污水处理厂有95%都无法达到《城镇污水处理厂污染物排放标准(GB18918-2002)》中的一级A标准，即使城镇污水处理厂出水达到了一级A的标准，也与地表水四类标准有着较大的差距，当总磷浓度超过0.1 mg/L (如果磷是限制因素)或总氮浓度超过0.3 mg/L (如果氮是限制因素)时，水体就可能出现富营养化 [ 1 ] [ 2 ] ，进而引起藻类及其它浮游生物迅速繁殖，水体溶氧量下降，鱼类及其它生物大量死亡。因此尽可能的降低进入水体中的氮、磷浓度是缓解水体富营养化十分有效和必要的措施之一，其中磷可以通过混凝沉淀等物理化学方法去除 [ 3 ] - [ 5 ] ，而要想通过物理方法去除污水中的氮素只能通过反渗透膜法 [ 6 ] ，然而极高的处理成本使得这种物理去除氮素的方法并不具有可行性，而经济可行的方法只有通过生物法去除。目前针对实际污水高效脱氮特别低碳氮比污水的研究取得了巨大的进展，主要技术有短程硝化–反硝化脱氮 [ 7 ] 、同步硝化反硝化脱氮 [ 8 ] [ 9 ] 、厌氧氨氧化(ANAMMOX)脱氮 [ 10 ] [ 11 ] 等，然而这些技术研究基本还处在实验室小试或者中试阶段 [ 12 ] ，距离实际应用还很远，所以在原有传统生物脱氮处理工艺基础上研究高效生物脱氮的策略与方法，是我国污水处理厂升级改造面临的迫切问题。 ICEAS全称为间歇式循环延时曝气活性污泥法，是SBR工艺中最具代表性的变种工艺之一，1976年美国建成世界上第一座ICEAS工艺污水处理厂，随后在日本、美国、加拿大、澳大利亚等地得到广发应用 [ 13 ] 。我国第一座ICEAS工艺污水处理厂昆明市第三水质净化厂于1994年建成，处理规模为15万吨/天，2007年8月昆明市第三水质净化厂二期建成，设计处理能力6万吨/天，至此处理能力达到21万吨/天，本文围绕昆明第三水质净化厂运行现状开展高效脱氮的策略的研究，为污水处理厂的高效脱氮升级改造提供理论支持。
system: [SYSTEM]

user: 2. ICEAS工艺特征 图1为ICEAS反应池构造简图。反应池分为预反应区和主反应区。预反应区具有缓冲连续进水和生物选择的功能：污水从预反应区以很低的流速(一般为3~5 cm/s)通过隔墙底部的过水孔进入主反应区，缓冲了进水对主反应区内泥水分离效果的影响；污水进入预反应区中与污泥接触，形成较高的污泥负荷，此时菌胶团细菌为优势菌，从而有效抑制丝状菌的生长进而防止污泥发生丝状膨胀 [ 14 ] 。 ICEAS反应池内不存在严格的厌氧环境，聚磷菌是一类特殊的兼性细菌，在厌氧条件下，聚磷菌可以利用体内储存的能量，将外部的易生物降解有机碳源吸收进体内，储存碳源及能源，同时消耗体内的ATP，将磷酸根释放到水体内；在好氧条件下，聚磷菌利用在厌氧阶段储存的碳源或能源进行新陈代谢，通过吸收外界磷酸根合成能源。在厌氧阶段吸收的碳源越多，在好氧阶段氧化分解产生的能量越多，从外界吸收的磷越多。所以ICEAS工艺的除磷功能特别弱，仅通过生物二级处理很难使磷达到较高的去除率，但磷可以通过深度处理系统添加化学药剂(PAM和PAC)经过混凝、沉淀去除，总之，磷的去除对污水处理厂来说并非难事。 ICEAS反应池通过缺氧时段和好氧时段的交替运行具有脱氮的功能。以昆明第三水质净化厂为例，其运行模式为连续进水、间歇排水两段A/O式运行模式，每天运行5周期，每周期4.8小时，具体运行模式为：第一次搅拌36 min、第一次曝气56 min、第二次搅拌20 min、第二次曝气56 min、沉淀50 min、滗水70 min，如图2。 通过以上运行模式，ICEAS反应池可实现进水中总氮的去除，然而在城镇污水厂实际运行中发现出水总氮难以达到一级A标准，以昆明第三水质净化厂为例，主要有以下几个因素限制了总氮的去除。 1) 进水碳氮比较低，反硝化阶段碳源不足。进水碳氮比较低是我国城镇污水处理厂面临的共同问题，进水碳氮比(COD/TN)为5左右，进水COD范围为250~350 mg/L，其中溶解性COD的浓度只有100 mg/L左右，污水进入ICEAS反应池后COD浓度由于稀释作用进一步降低，第一次搅拌阶段和第二次搅拌阶段的缺氧反硝化反应由于碳源的缺失经行的并不充分； 图1. ICEAS反应池构造简图 图2. 昆明第三水质净化厂运行示意图 2) 两段A/O式运行模式没有发挥其应有作用。ICEAS反应池反硝化主要发生在第一次搅拌阶段，经过第一次曝气阶段好氧曝气后有机碳源COD浓度降到40 mg/L左右，在第二次搅拌阶段缺氧搅拌时水中的碳源严重不足，致使其反硝化功能基本丧失； 3) 连续进水方式使得进水中的碳源没有被高效利用。第一次曝气阶段和第二次曝气阶段好氧曝气时依然在进水，此时进水中的有机碳源仅仅通过好氧曝气方式去除，即原水中约40%的碳源在曝气阶段被氧化降解，而并没有用作反硝化碳源，这种连续进水的方式不仅使进水中的碳源利用率较低，同时还增加曝气阶段系统的有机物负荷，造成了不必要的能源浪费。 为解决以上存在的ICEAS工艺总氮去除率不高的问题，对其展开SBR工艺小试装置实验研究，重点研究高效脱氮策略，将连续进水模式改为间歇进水模式，同时将两段搅拌方式合为一段，并且增加A/O比，最终确定最佳运行模式，为污水处理厂的高效脱氮升级改造提供理论支持。
system: [SYSTEM]

user: 3. 实验材料和方法 3.1. 小试污泥和污水 SBR小试用水取自昆明第三水质净化厂一期(老厂)一级配水渠道，其水质特征为COD (42.5~181.9) mg/L，BOD 5 (20.3~92.4) mg/L， -N (25.6~38.6) mg/L，NO 3 -N (0~0.52) mg/L，TN (27.4~39.6) mg/L，COD/TN (1.3~5.7)。该生活污水平均COD/TN为2.9，碳氮比较低；SBR小试所用污泥直接取自昆明第三水质净化厂ICEAS池。 3.2. SBR小试装置 实验所用SBR小试试验装置由有机玻璃制成，圆柱形，高36 cm，直径20 cm，总有效容积12 L，如图3所示，进水通过蠕动泵控制，根据需要由加热棒来调节温度，采用机械搅拌，装置底部用粘砂块作为微孔曝气头，由空气泵进行曝气。反应装置内设有DO、pH和ORP探头，可以实时监测反应过程中的DO、pH和ORP值及其变化规律。 3.3. 检验分析项目 试验检测项目及方法见表1所示，测样前水样都经过滤纸过滤。 3.4. 实验方案 SBR小试实验装置分为A、B两组，且在相同的污泥浓度(5000 mg/L左右)和进水条件下同时运行。 A组模拟目前昆明第三水质净化厂一期的运行模式，具体如图2所示，进水方式为连续进水，每天运行5个周期，每周期4.8小时； B组模拟间歇进水运行模式，具体如图4所示，进水方式为间歇进水(沉淀、滗水、搅拌阶段进水而曝气阶段不进水)，每天运行5个周期，每周期4.8小时，具体运行周期分为二个阶段，具体如表2所示。
system: [SYSTEM]

user: 4. 结果与讨论 4.1. COD去除效果分析 SBR小试A、B两组运行模式的进水条件相同，进水的COD浓度范围为92.5~187.1 mg/L，平均浓度为142.7 mg/L。 A组连续进水运行模式的出水COD平均浓度为24.8 mg/L，COD的平均去除率为82.7%。如图5所示； 图3. SBR小试试验装置图 图4. (a) B组阶段一运行模式；(b) B组阶段二运行模式 图5. A组COD去除情况 Table 1 Analysis of items and method 检测项目 分析方法 检测项目 分析方法 COD 快速消解法(兰州连华) DO WTW Multi340i -N 纳氏试剂光度法 ORP WTW Multi340i TN 过硫酸钾氧化-紫外分光光度法 MLSS 滤纸称重法 表1. 检测项目与分析方法 Table 2 The operating mode of stage two of group 项目 搅拌 (min) 曝气 (min) 沉淀 (min) 滗水 (min) 阶段一 48 120 50 70 阶段二 56 112 50 70 表2. B组二阶段具体运行模式 B组间歇进水运行模式阶段一的出水COD平均浓度为25.0 mg/L，COD的平均去除率应为81.7%；阶段二的出水COD平均浓度为13.9 mg/L，COD的平均去除率为83.9%。如图6所示。 由上可以看出，将连续进水模式改为间歇进水模式后，对COD的去除效果并没有产生负面的影响，其出水COD的平均浓度仍然小于50 mg/L，出水水质有机物指标可以稳定达到城镇污水排放一级A标准。 4.2. 氨氮去除效果的分析 SBR小试A、B两组运行模式进水氨氮的浓度范围为20.3~40.3 mg/L，氨氮平均浓度为31.0 mg/L。 A组连续进水运行模式的硝化效果较好，出水氨氮的平均浓度为0.5 mg/L，氨氮的平均去除率达到98.6% (如图7所示)。 B组间歇进水运行模式阶段一出水氨氮平均浓度为0.2 mg/L，氨氮的平均去除率为99.7%；阶段二的出水氨氮平均浓度为0.3 mg/L，氨氮的平均去除率为98.8% (如图8所示)。 通过分析，将连续进水模式改为间歇进水模式后，系统的硝化能力并没有减弱，间歇进水模式A/O比分别为1:2.5和1:2时，依然都能够实现很高的氨氮去除率，其出水氨氮的平均浓度均小于1.0 mg/L，出水水质氨氮指标都可以稳定达到城镇污水排放一级A标准。 4.3. 总氮去除效果的分析 SBR小试A、B两组运行模式进水总氮的浓度范围为31.8~41.9 mg/L，总氮平均浓度为31.8 mg/L，为消除不同除碳氮比对两种不同运行模式总氮去除率的影响，本实验严格控制碳氮比在4.5左右。 A组连续进水运行模式出水总氮的平均浓度为13.8 mg/L，总氮的平均去除率为56.1% (如图9所示)。 B组间歇进水运行模式阶段一A/O比为1:2.5时：出水总氮平均浓度为12.0 mg/L，总氮平均去除率为62.9%，此时A组连续进水运行模式出水总氮平均浓度为13.3 mg/L，总氮平均去除率为58.0%； B组间歇进水运行模式阶段一A/O比为1:2时：出水总氮平均浓度为11.0 mg/L，氨氮的平均去除率为65.0%，此时A组连续进水运行模式出水总氮平均浓度为14.4 mg/L，总氮平均去除率为56.0%。 通过以上可以看出，在碳氮比为4.5左右时，间歇进水模式的总氮去除率始终高于连续进水模式。 间歇进水模式在A/O比为1:2.5时，总氮去除率比连续进水模式高4.9%；间歇进水模式在A/O比为1:2时，总氮去除率比连续进水模式高9%。 间歇进水模式A/O比为1:2时比A/O比1:2.5时的总氮去除率高2.1%。 图6. B组COD去除情况 图7. A组氨氮去除情况 图8. B组氨氮去除情况 图9. A、B两组总氮去除情况
system: [SYSTEM]

user: 5. 结论 SBR工艺起初是为去除碳源有机物开发的，后来随着处理要求的提高，增加了脱氮功能，ICEAS工艺虽然在传统SBR的基础上做出了改进，进水方式改为连续进水后简化了配水的设施，也省去了对配水的控制，从而使得ICEAS工艺可用于较大型的污水处理厂，然而其脱氮的功能并没有得到强化，为进一步提高ICEAS工艺污水处理厂的总氮去除率，还需要对其进行进一步的优化改进。 1) 连续进水方式不利于ICEAS工艺总氮的去除，改变进水模式，即搅拌前集中进水，系统内的微生物可最大程度的利用进水中的碳源进行反硝化，在提高系统反硝化效率的同时也为后继的曝气阶段减轻有机物负荷，适当降低能耗。曝气阶段停止进水，系统内留有大量的氨氮和少量的有机物。此时，硝化菌处于竞争优势，可快速利用水中溶解氧气对氨氮进行彻底的氧化降解，系统的硝化效果得到加强。将连续进水方式改为间歇进水方式后，总氮去除率平均可提高至少5%； 2) 适当提高ICEAS工艺的A/O比，将A/O比从1:2.5提高到1:2后，总氮去除率得到进一步的提高，但当进一步提高A/O比出水氨氮浓度有增加的趋势，所以最佳A/O比为1:2； 3) 连续进水模式改为间歇进水模式后，出水COD、氨氮浓度均可以稳定达到城镇污水排放一级A标准，而磷由于后期的深度处理，出水浓度也可以实现城镇污水排放一级A标准。
system: [SYSTEM]

user: 摘 要 欧姆接触是碳化硅(Silicon Carbide， SiC)新一代电力电子器件研究中的技术难点之一。除了金属选择外，相对于n型掺杂，SiC材料中p型杂质的离化能比n型杂质的离化能高，优质的p型SiC欧姆接触更难于形成。该文对近十几年来极具代表性的传统Al基金属体系和非传统Al基金属体系在p型SiC材料上形成欧姆接触的研究进行了总结，并对其发展前景进行了展望。 关键词 :p型SiC，电力电子器件，欧姆接触
system: [SYSTEM]

user: 1. 引言 SiC电力电子器件主要包括功率整流器(SBD、PiN和JBS等)、单极型功率晶体管(MOSFET、JFET和SIT等)和双极型载流子功率晶体管(BJT、IGBT和GTO等)。其中极具代表的SiC PiN、SiC MOSFET和SiC BJT等器件，均需要形成良好的p型欧姆接触。经过大量的研究工作，SiC欧姆接触问题目前主要集中在两个技术点上：一是大剂量的离子注入，这能显著地增加SiC的载流子密度，并减少它的耗尽层宽度，形成隧穿；这种方法的关键问题是：在离子注入过程中容易形成晶格缺陷或非晶状态，这些缺陷是非常顽固的，需经过大约2000K的极高温退火来复原，从而使制造SiC器件复杂化 [ 1 ] ；二是通过淀积金属和退火技术，产生一种具有较窄带隙或高载流子密度的中间层。为形成这种层，大范围的材料可供选择，包括金属化物、硅化物、碳化物、氮化物和石墨等等 [ 2 ] 。实际欧姆工艺制备是在综合以上两个技术点的情况下完成的。
system: [SYSTEM]

user: 3. 总结和展望 p型SiC欧姆接触的研究是在缺乏高质量的p型SiC欧姆接触的背景下被逐渐放到研究前沿的，该文总结了近十几年来极具代表性的关于传统Al基金属体系和非传统Al基金属体系在p型SiC材料上形成欧姆接触的研究。目前p型SiC欧姆接触的研究主要是通过掺杂工艺结合沉积退火工艺来实现。掺杂浓度均高于10 18 cm −3 ；传统的TiAl接触退火温度均高达1000℃左右，而高的退火温度会引起了极粗糙的表明形态，且传统Al基金属体系是不可靠的，在传统TiAl两层金属体系中加入Ni、Ge、Au等金属后，有利于退火温度的降低及稳定性的改善。截止目前，大量p型SiC欧姆接触研究报道出的比接触电阻值在10 −3 ~10 −5 Ωcm 2 范围之内，和相对成熟的n型SiC欧姆接触比接触电阻值依旧存在着差距。在可预见的将来，p型SiC欧姆接触研究的方向主要还在TiAl基金属体系的选择以及掺杂激活工艺的改进上。相信随着工艺技术的不断改进和完善，p型SiC欧姆接触的难题一定能够成功解决。
system: [SYSTEM]

user: 摘 要 本文研究第一类多环芳香烃的化学拓扑指数。对第一类多环芳香烃的分子结构进行分析，计算出每一对顶点之间的距离，根据广义度距离和维纳相关指数的计算公式得到相应的结果。 关键词 :多环芳香烃，分子结构，广义度距离
system: [SYSTEM]

user: 1. 引言 从化学实验的角度研究一种新的物质，需要有良好的实验环境支撑，有充足的资金来购买仪器设备和试剂，同时还将花费大量人力物力。而用图模型研究分子结构性质的方法，就能避免这些问题，在无硬件设备的条件下获得关于化学物的化学、药学特征。在具体计算模型中，将化学分子的结构用图模型来表示，化学分子的特性就可以用一些参数来衡量，比如PI指数、维纳指数、舒尔茨指数等等。文献 [ 1 ] 讲述了拓扑指数在生物学和化学的使用开始于1947年，当时化学家哈罗德.维纳证明有机物化合分子的理化性质之间的化学图表，维纳指数所有无序对简单图G的顶点之间的距离之和，得出维纳指数。文献 [ 2 ] 根据1988年对于拓扑指数Hosoya引入Hosoya多项式，定义了维纳指数和相关多项式。Schultz在文献 [ 3 ] 中定义了基于结构描述符“分子拓扑指数”(即舒尔茨指数)。Klavžar和Gutman在文献 [ 4 ] 中定义了修改舒尔茨指数。文献 [ 5 ] - [ 10 ] 研究和计算了舒尔茨指数和修改的舒尔茨指数。其他关于拓扑指数计算以及在化学工程中的应用可参考文献 [ 11 ] - [ 20 ] 。 多环芳香烃(Polycyclic Aromatic Hydrocarbons，简称PAH或PAHs)又称多环性芳香化合物或多环芳香族碳氢化合物，是芳香族碳氢化合物的一种特例。由不包含杂环或取代基的芳香环所组成。第一类多环芳香烃简称PAH1，是多环芳香烃家族的第一个成员，俗称苯环，如图1所示。 苯环是由6个sp2杂化碳原子通过 键和 键构成平面正六边形的碳环，分子式为C 6 H 6 ，苯环是最简单的芳香烃环，由六个碳原子构成一个六元环，每个碳原子接一个基团，苯的6个基团都是氢原子。 图1. 第一类多环芳香烃 苯分子中6个碳原子各以3个sp2杂化轨道分别跟相邻的两个碳原子的sp2杂化轨道和氢原子的1s轨道重叠，形成6个碳碳 键和6个碳氢 键。两个sp2杂化轨道的夹角是120˚，正适合6个碳原子处于一个平面上，形成一个正六边形的苯环。苯环上6个碳原子各有一个未杂化的2p轨道，6个2p轨道的对称轴都垂直于环的平面，并从侧面相互重叠，形成一个闭合的 键。它均匀地对称分布在环平面的上方和下方。通常把苯的这种键型叫做大 键。苯分子中 键电子云完全平均化，使苯环中每个碳碳键的键长和键能都是相等的。这就说明苯分子的对称性和稳定性。苯环的主要化学特性是环平面上下的 键电子容易受到亲电试剂的进攻，结果通常发生环上的取代反应。由于苯环较稳定，较难发生环上的加成反应。 下面介绍本文关注的一些重要指数 一个分子图的顶点带权维纳数定义为： 。 倒数度距离定义为： 。 倒数乘法度距离定义为： 。 单顶点带权维纳多项式定义为： 。 双顶点带权维纳多项式定义为： 。 广义度距离定义为： 。 其中 为实数。 另外，修改的广义度距离则定义为： ， 其中 是实数 它们对应的多项式为： ， 。 修改的维纳指数是原维纳指数的扩展，定义为 ， 其中 是非零实数。可见当 时，修改的维纳指数即为维纳指数。 修改的超维纳指数是超维纳指数的扩展，定义为 。 乘法维纳指数是维纳指数的另一种表现形式，定义为 。 此外，对数乘法维纳指数定义为 。 Harary指数是一类和维纳指数相关的指标，定义为 。 对应的Harary多项式则定义为 。 上述Harary指数通常称为第一类Harary指数，而第二类和第三类Harary指数则定义为 ， 。 更一般地，广义Harary指数则定义为 ， 其中 是一个非负整数。 本文将给出第一类多环芳香烃的上述化学指数计算公式。
system: [SYSTEM]

user: 2. 主要结论及证明 设 为第一类多环芳香烃，它由6个碳原子(C)和6个氢原子(H)组成。假设 和 是分子图上的两个顶点， 和 分别表示顶点 和 的度。设 是分子图 中顶点 和 之间的最短距离。 表示分子图中距离为 的顶点对的集合。顶点原子 和 之间的最短路径的原子集合数为66，顶点之间的最短距离有从一到五一共五种。即 ，显然， 。下面，我们将把它分为五部分来一一列举。 2.1. 顶点u和v之间最短距离为1 6个碳原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，36，54。 6个氢原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，24，18。 从整个分子图来看，第一项 ， ，所以，维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为12，60，72。 2.2. 顶点u和v之间最短距离为2 6个碳原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为12，72，108。 6个氢原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，24，18。 从整个分子图来看，第二项 ， ，所以，维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为18，96，180。 2.3. 顶点u和v之间最短距离为3 6个碳原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为3，18，27。 6个氢原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，12，6。 所有碳原子与氢原子之间：最短距离为3的顶点集合数为12。 从整个分子图来看，第三项 ， ，所以，维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为21，78，69。 2.4. 顶点u和v之间最短距离为4 6个碳原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，12，6。 6个氢原子之间： ， 则维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为6，24，18。 从整个分子图来看，第四项 ， ，所以，维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为12，36，24。 2.5. 顶点u和v之间最短距离为5 从整个分子图来看，第五项 ， ， 所以，维纳多项式、广义舒尔茨多项式和广义修改的舒尔茨多项式中第一项的系数分别为3，6，3。 通过列举PAHs所有原子不同顶点之间的最短距离，可利用指数计算公式得到如下关于第一类多环芳香烃的指数计算公式。 顶点带权维纳数为： 。 倒数度距离： 倒数乘法度距离为： 单顶点带权维纳多项式为： 双顶点带权维纳多项式为： 广义维纳指数为： 其中 是实数。 超维纳指数为： 。 修改的超维纳指数为 。 乘法维纳指数为 。 对数乘法维纳指数为 。 Harary指数是一类和维纳指数相关的指标为 对应的Harary多项式为： 第二类Harary指数为： 。 第三类Harary指数为： 。 广义Harary指数为 。 其中 为非负整数。 广义度距离和修改的广义度距离对应的多项式为： ， 。 广义度距离为： 。 修改的广义度距离为： 其中 是实数。
system: [SYSTEM]

user: 3. 总结 通过对第一类多环芳香烃的分子结构分析，可以计算出每一对顶点之间的距离，再根据指数计算公式得到相应的结果。本文理论结果可以为其它化合物的分析计算提供参考。
system: [SYSTEM]

user: 摘 要 拓扑指数是分子结构数值化的一种方式，它能充分反映分子图的连接信息和化学环境，能有效地表达化合物的结构与性质。本文利用边集划分技术，得到Pent-Heptagonal纳米管若干基于度的拓扑指数。 关键词 :Pent-Heptagonal纳米管，分子图，广义Randić指数，广义Zagreb指数
system: [SYSTEM]

user: 1. 引言 在化学图论中，我们用一个图来描述分子的拓扑结构，用顶点表示原子，边表示化学键，该图称为分子图。拓扑指数是从化合物的结构图衍生出来的不变量。它不但可以定量的描述分子的结构，而且可以分析相关分子的结构与性能之间的关系。相关内容可参考 [ 1 ] - [ 12 ] 。 用 表示顶点/原子 的度。分子图 的拓扑指数是分子结构数值化的一种方式。1972年Trinajstić和Gutman引入了的第一Zagreb指数和第二Zagreb指数： ， 。 1975年，美国化学家Milan Randić [ 2 ] 提出了最古老的连通性指数，定义为 1988年，Bollobas和Erdos在 [ 6 ] 中推广了Randi指数，用任意实数k代替了，从而定义了广义的Randić指数： ， 其中k是非零实数。以下若没有特别说，k均表示非零实数。 当k = 1时，广义的Randić指数就是我们通常说的第二Zagreb指数。当k = −1时，广义的Randić指数就是修改的第二Zagreb指数。 Azari和Iranmanesh [ 13 ] 引入广义Zagreb指数，定义为 ， 其中 和 为任意的非负整数。 和连通指数则定义为： 推广形式的和连通指数为 。 而调和指数则定义为： 它的广义形式为： 。 Eliasi和Iranmanesh [ 14 ] 定义广义几何代数指数如下： 。 当k = 1时，即为一般GA指数如下： 。 Estrada等人 [ 15 ] 定义原子键连通指数(atom-bound connectivity index，简称ABC指数)如下： 。 两类Zagreb指数对应的多项式分别为 ， 。 此外，第三类Zagreb指数和多项式定义为 ， 。 Gutman [ 16 ] 给出分子图的乘法Zagreb指数定义如下 ， 。 其中 为第一类乘法Zagreb指数， 为第二类乘法Zagreb指数。 Usha等人 [ 17 ] 重新定义了三类Zagreb指数如下： ， ， 。 Furtula和Gutman [ 18 ] 定义了F指数如下： 。 设 为满足 和 的边e = uv的个数，则M-多项式定义为 。 广义ABC指数 和and增强的Zagreb指数 则定义为 [ 19 ] ： ， 。 最后，减弱的倒数Randić指数定义为 。 第二减弱的倒数Randić指数定义为 。 Pent-Heptagonal纳米管 是由周期长度为5的 (五边形)和周期长度为7的 (七边形)相互交替构成的网状化合物。设 纳米管的2维拉丁网格中，第一行和列的7边形个数分别为m和n (见图2)。将 纳米管记为 。 纳米管的圆柱和2维拉丁网络可参考图1和图2。 图1. Pent-Heptagonal纳米管的圆筒晶格 图2. Pent-Heptagonal纳米管的二维晶格 本文利用边集合划分方法得到Pent-Heptagonal纳米管若干基于度的拓扑指数计算公式。
system: [SYSTEM]

user: 2. 主要结果 将Pent-Heptagonal纳米管的顶点集和边集合作如下划分： ， ， ， ， 。 从图2来看，在Pent-Heptagonal纳米管 中分别用m和n来表示第一行和第一列的七边形数量，则有纳米管的顶点数 ，并且有： 并且，我们在图2中用黄色标记 、 的边，用红色标记 、 的边， 用白色标记。因此，得到 的大小等于 ， 的大小为 。换句话说，以上边的数量划分为 和 ，我们有以下的表格： 于是，利用上表以及各种拓扑指数的定义，我们得到： 广义的Randić指数为： 修改的第二Zagreb指数为： 广义Zagreb指数为： 和连通指数为： 推广形式的和连通指数为： 调和指数为： 调和指数的广义形式为： 广义几何代数指数为： 。 一般GA指数为： 。 原子键连通指数为： 。 两类Zagreb多项式分别为 ， 。 第三类Zagreb指数和多项式分别为： ， 。 两类乘法Zagreb指数分别为： 。 重新定义的三类Zagreb指数分别为： ， ， 。 F指数为： 。 M-多项式为 。 广义ABC指数 和and增强的Zagreb指数 分别为： ， 。 两类减弱的倒数Randić指数分别为 。 第二减弱的倒数Randić指数定义为 。
system: [SYSTEM]

user: 3. 总结 随着化学实验条件的改善和新化学生物技术的发现及应用，每年有大量复杂的化学分子被发现。在化学领域，化学家们会通过大量实验，来研究新物质的性质，比如：熔点、沸点、毒性、氧化性等等。因而对每年发现的大量新物质的研究需要大量仪器设备和药品试剂，造成了研究的高成本。但是，如果把化学分子的结构用图模型来表示的话，化学分子的特性就可以用一些参数来衡量，比如Zagreb指数、维纳指数、Randić指数等。本文利用边集合划分方法得到的Pent-Heptagonal纳米管若干基于度的拓扑指数计算公式，计算结果可运用到新化学分子的性质研究上，进而可以大大节约研究成本以及研究时间。
system: [SYSTEM]

user: 摘 要 化学图论其模型是用图结构来表示分子：原子用顶点来表示，原子之间的化学健用边来表示。本文利用图分析的方法计算Harary图H 2m，n 每对顶点的距离，并根据拓扑指数的定义给出维纳相关指数的表达式。 关键词 :Harary图，修改的维纳指数，Harary指数，乘法维纳指数
system: [SYSTEM]

user: 1. 引言 在化学图论和分子拓扑学中，对于不同的拓扑指数虽然可以用不同的方法计算得到，但其基本的思想及方法基础是一样的，那就是要确定原子的属性或分子图节点的化学特征(例如采用原子的分子连接度和原子共价半径等)、原子间的连接关系(键属性)以及利用各种分子图不变量(例如顶点数目、边数目、等长度路径数目、连接矩阵顶点度、分子图的特征多项式和图谱等) [ 1 ] 。 随着化学实验条件的改善和新化学生物技术的革新，大量的新化合物和药物涌现出来。近年来，有太多新发现的化合物、材料和药物急剧增长。这些新化学结构，都要鉴定它们的物理、化学和生物特征，这使得化学和制药领域的研究者的工作量大大增加。怎样才能精确、有效、快捷地知道新物质的特性成为现在科研的一个新的挑战。更多详细信息可以参看文献 [ 2 ] - [ 6 ] 。 研究实验表明：化合物和药物的部分特性离不开化学分子结构。所以，部分学者在分子结构图上定义了一些指标，如维纳指数，超维纳指数等，通过计算这些指数来得知它们的化学性质。这种方法可以大大减轻了鉴定药物特性的工作量，可以直接通过计算得出这些化学物质的特性，免去了化学实验，节约成本。关于指数信息可以参看文献 [ 7 ] - [ 14 ] 。 在理论化学中，其计算模型是用图结构来表示分子结构：其中原子用顶点来表示，原子之间的化学健用边来表示。一般地，设 是一个图，用顶点集合 即为原子集合，边集合 则为化学键集合。 修改的维纳指数是原维纳指数的扩展，定义为 ， 其中 是非零实数。可见当 时，修改的维纳指数即为维纳指数。 修改的超维纳指数是超维纳指数的扩展，定义为 。 乘法维纳指数是维纳指数的另一种表现形式，定义为 。 此外，对数乘法维纳指数定义为 Harary指数是一类和维纳指数相关的指标，定义为 = 。 对应的Harary多项式则定义为 。 上述Harary指数通常称为第一类Harary指数，而第二类和第三类Harary指数则定义为 ， 。 更一般地，广义Harary指数则定义为 ， 其中 是一个非负整数。 用 来表示Harary图，其中 表示Harary图的顶点数。若将其顶点表示成 ，则两顶点用 之间存在边当且仅当 。图1给出当 时的Harary图。 本文将Harary图 的顶点对距离进行分析，并给出修改的维纳指数，Harary指数和乘法维纳指数的计算公式。
system: [SYSTEM]

user: 2. 主要结论及证明 根据Harary图 的定义，可知 。图中取两个顶点 且满足 ，它们之间的距离为 ，因此这个正则图边的数量为 。另一方面从定义可知所有 图的顶点都有相似的几何特征和拓扑结构，距离i点所有的路径总数就是顶点的个数。例如对 ，有2m条边，即有2m条距离为1的路径。从图的定义中可以看到 和 ， 因为 和 。显然 。因为 ， 满足 。因此 。显然 。由此可知，对于Harary图上的任意顶点 ，图的直径 就是 。同样顶点 和 之间的距离也等于 。对于所有的 。容易看出如果 是奇数就由 ， ，当n为偶数时 。所以当n为奇数时对应的Hosoya多项式为： 当n为偶数时，设n = 2q，对应的Hosoya多项式为： 图1. Harary图H 6，10 综合以上分析和计算，当 为奇数时， ， ， ， ， ， ， 。 当 为偶数时， ， ， ， ， 。