尿素节能改造及优化控制小结（精选合集）

第一篇：尿素节能改造及优化控制小结

尿素装置节能改造与优化操作小结

刘肇庆

（山东阳煤恒通化工股份有限公司

山东郯城

276100）

一、基本概况

山东阳煤恒通化工有限公司现有二套化四院设计年产四万吨、六万吨水溶液全循环法尿素装置各一套，近年来采用同行业先进的尿素系统化改造技术，分步骤对尿素装置进行了技术改造，已具备日产800吨的生产能力，各项消耗指标均有所下降。

二、系统节能改造情况 1.高压系统

1.1在原设计为40kt/a和60kt/a生产能力的基础上，先后对两套尿素进行双塔并联流程改造，由于产能的扩大，合成塔生产强度的提高，相应合成塔转化率有所下降，通过应用新型高效塔板，在塔板上增加气泡式泡罩，使气液相接触反应的几率增加，并加快反应热的扩散速度，从而提高第一反应式的完成程度。同时增加了塔板的安装数量，增加了塔内的反应区段，从而提高了气液传质效果。改造后，在高生产强度下，合成转化率可提高到67%，生产能力提高的同时，中压系统负荷也大大降低。

1.2入塔立式止回阀全部更换成卧式止回阀，不仅方便了检修，降低了维修费用，而且使用效果较好。

1.3两套尿素各增加一台14m3/h的3JA-14/21-TB一甲泵和一台20m3/h的3A-20/21-T的液氨泵。2.中压系统

2.1中压分解系统改造。中压分解系统的改造重点通过新技术的应用改造一段分解系统，降低一段分解蒸汽耗。原一段分解采用预蒸馏流程，2000年5月，2#尿素新上一台φ800mm预分离器，采用预分离—预蒸馏流程，即合成塔出来的反应熔融物，经减压绝热膨胀，进入预分离器，分解所需的热量由合成塔出来的熔融物自身降温提供，使大部分游离氨得到分离。提高了CO2的吸收率，大大降低了一吸塔的热负荷，减少了一吸塔塔顶和塔底回流氨用量，同时一分塔蒸汽消耗也明显下降。但是随着负荷的增加，由于预分离器容积较小（￠800×2600，V=0.8m3），安装位差达不到要求，在运行中预分离器气相带液严重，达不到理想的预分离效果。2006年11月，2#尿素新上一台φ1200mm F=333m2的自气提式一分塔，分上下两段，上段为五层GC型高效翅片式塔板，下段为降膜式换热器。改造开车后得到了较好的效果：一分塔气相温度由原来的130℃下降为123℃，一吸塔运行稳定；吨尿素蒸汽消耗下降约50Kg。由于效果显著，在2009年6月1#尿素进行了同样改造。

2.2中压吸收系统改造。中压吸收系统的改造重点在于中压吸收塔的优化，扩大一吸冷却器面积，尽可能的降低一吸塔精洗段的负荷，并相应增加氨冷器面积，满足扩能的需要。一吸塔是水溶液全循环法工艺的“心脏”，其塔盘泡罩高度，鼓泡段容积、分布器开孔情况对热平衡和物料平衡影响很大。1#尿素原φ1000mm的一吸塔更换为φ1200mm的一吸塔，2#尿素在原φ1100mm一吸塔下部加长了2m，对 分布器重新开孔；对一吸冷却器设备流程进行优化改造，1#尿素增设一台φ600mm、F=71m2的第二一吸冷却器，与原71m2串联；2#尿素增设一台φ800mm、F=90m2的第二一吸冷却器及一台φ1000mm、F=150m2的第三一吸冷却器，与原71m2串联，使一吸塔的吸收负荷外移。从而提高了一段吸收系统的生产能力，使一吸塔内的吸收负荷下降，减少了顶底回流氨用量，增加了一吸塔调节裕量，进一步提高了一吸塔的操作稳定性和生产能力。

2.3对氨冷器进行了改造，1#尿素增设一台φ800mm、F=220m2氨冷器；2#尿素增设一台φ800mm、F=110m2氨冷器，同时把氨冷器A由原来F109m2更换为F350m2的氨冷器，增加一根φ57mm气相平衡管，平衡了氨冷器的负荷，提高了氨冷器的冷却效率，从而提高了整个吸收系统的生产能力。解决了以前因水温高造成氨冷器积液，引起工况波动，为稳产高产提供了保证。

2.4两套尿素的二甲泵由原来2台增加到3台；改变了原来二循二冷溶液中CO2含量超标的现象，吨尿素氨耗进一步下降。

2.5对两套尿素的脱盐水冷却系统进行改造。1#尿素增设一台换热面积为139m2卧式脱盐水换热器。2#尿素增设一台换热面积为200m2立式脱盐水换热器，提高了换热效率，吨尿素脱盐水消耗由3.0m3左右降到2.0 m3以下。

2.6更换了两套尿素框架内的原料氨管线，对两套系统液氨缓冲槽液位调节阀、中压压力调节阀进行了改造，由六楼移到三楼，降低了操作人员的劳动强度，提高了工作效率，减少了工况波动。3.低压及解析系统

为降低改造的投资费用，低压分解系统基本利用一段分解系统设备进行改造，扩大低压分解系统的能力，并利用新技术，降低分解消耗，降低低压吸收系统的负荷。

3.1把原一分塔改为二分塔使用。并把一分加热器液位槽的膨胀蒸汽送给二分加热器使用，不足部分由中压蒸汽管网提供。1#尿素二循一冷更换为F=220m2，原F=168m2的二循一冷作为二循二冷使用。利用废旧的闲置设备，在二循一冷前增设预冷器，减轻了二循一冷和二循二冷的吸收负荷。改造后，不仅完全满足了高负荷生产的要求，而且吨尿素氨耗和蒸汽消耗进一步降低。

3.2对尾吸系统进行了改造，采用双尾吸流程。在常压尾吸塔前，增设0.4Mpa带压尾吸，提高了尾吸吸收效率，减少了尾吸塔排放气中的氨含量，由原来设计的8%降至1%以下，尿素氨耗进一步降低。

3.3对尿素深度水解装置的回流冷凝器流程进行改造，改变了以往因回流液温度高，解析超压问题。深度水解回流液温度控制在45~50℃送二循一冷。处理后含氨≤20ppm，尿素≤5ppm的废水送往水解废液和甲醇残液回处理系统，经处理合格后给造气夹套汽包使用，既节约了脱盐水，又提高了经济效益和环保效益。

3.4对高压柱塞泵填料密封水进行改造。正常运行时，氨泵填料密封水、一甲泵填料密封水、熔融泵密封水、二甲泵、氨水泵的填料密封水收集后送往碳铵液槽，经解析泵泵送入水解系统，既控制了环保排放，又降低了氨耗。同时增加了φ5000mm、V=100 m3和V=50 m3 的碳铵液事故槽，供生产不正常时使用。

4.蒸发造粒及包装系统

4.1二蒸分离器进行改造，改为旋风分离器，二蒸加热器的换热面积由原来7.5m2增加为F=12.85 m2，使用钛材材质，有效的延长了加热器的使用寿命，由原6个月左右的使用周期到目前使用10年。

4.2一尿的闪蒸由原来￠500mm改为￠700mm。并增设了闪蒸冷凝器，有利于蒸发真空度的控制。

4.3对表冷器流程进行优化改造。一表冷进行更换由F=112m2改为F=207m2，新增了二表水冷，与原二表水冷串联，改造后，保证系统的冷却效果，满足了生产能力提高的需要。

4.4一尿二表氨冷改用冷却水冷却，因原二表氨冷换热面积偏小，在其后又串联一个第二表冷器，以满足二段真空度的要求。氨冷改用水冷后，吨尿素节约冷量2.63~3.05GJ，可节约合成氨冰机吨氨电耗1度。但近年来由于一次水资源比较紧张，夏季循环水温度高，氨冷改用水冷后，二段真空度难以控制。

4.5对二蒸下液管进行改造，对易出现冲刷腐蚀的二蒸出口管道进行加粗，由￠159mm改为￠219mm有效的减少了设备的腐蚀情况，同时提高了管道的使用材质，使用316材质的管道，延长了管道的使用寿命，由一年到现在的3年左右。

4.6增设水力喷射装置。原有的蒸发系统喷射器采用蒸汽作为动力，既浪费了蒸汽，又造成氨损失。经过考察论证，两套尿素各增1套水力双吸喷射装置，取代了原一段蒸汽喷射器和二段蒸汽喷射器。水力喷射循环槽用水由熔融泵密封水供给，循环提浓后送水解系统，这样不仅节约了密封用脱盐水，而且减轻了环保压力。改造后，吨尿素蒸汽消耗下降了50Kg，日回收氨0.5t。

4.7对蒸发系统的蒸汽冷凝液系统进行改造，一、二蒸加热器不使用疏水器疏水，使用了液封自流式液位槽，降低了蒸汽消耗。

4.8对造粒系统进行改造，有效的改善了造粒塔粘塔的现象，同时尿素质量也有所提高，保证了系统扩能的要求。

5.空压站循环水系统

5.1在2009年1月份对循环水流程进行改造，增设4台12SA-13A扬程为26米的低压循环水泵，压缩机及二循一冷、二循二冷改用低压水，保证了六楼氨冷器水量，避免了夏季氨冷积液现象，同时节约了水量，减少了开泵台数，节约了电耗。

5.2在2009年2月份循环水系统增设旁滤装置，循环水水质得到改善，保证了冷却器冷却效果。同时减少了排水量，节约了一次水用水量。

5.3在2007年12月份仪表空气增设了冷干机，增加一台100m3的空气储罐，解决了仪表空气带水问题，保证了调节阀门的稳定运行。

三、工况的优化控制

设备具备了降耗的基础，同时再以正确的调整方法，即可达到节能降耗的目的，优化工艺控制操作如下: 1.用水平衡的概念指导操作

水溶液全循环法工艺是指用水吸收未反应物，分解出来的氨、CO2，以甲铵液返回合成塔。循环的内容有：①过量的氨和未反应物中分解的氨循环；②未反应物分解后CO2的循环；③水的循环。水溶液全循环法工艺的核心是保持系统水平衡，水的来源有三个：①用于吸收分解气的水；②解析气相中带的水；③一、二分气相带来的水。如果操作不稳定，带进尿塔的水增多，则转化率下降，进入循环系统的未反应物增多，分解气中水蒸气也增多；一甲液浓度下降速度增加。致使操作难以控制，发展到恶性循环而不能控制的地步。因此，维持系统水平衡的关键在于控制吸收用水和分解气、解析气中的水含量。

根据系统工艺工况，系统水平衡有三个标志：①入塔H2O/CO2比在0.65~0.85；②一甲液中CO2含量在28~34%；③二分后尿液浓度为约67%。由工艺物料衡算可知，吨尿素一段进入合成塔的水330Kg，合成塔生成的水300Kg，最后由解析废液排出。

2.提高CO2转化率

二氧化碳转化率是尿素生产中的重要工艺指标，直接影响蒸汽消耗量和尿素产量，转化率每上升1%，吨尿素在一分加热器的蒸汽耗下降约40Kg。因此，提高尿素合成塔的转化率，是降低消耗的重要因素之一，首先维持好全系统水平衡，控制入塔NH3/CO2比在4.0，H2O/CO2比在0.65（最大限度不能超过0.85），塔温188温度，尿塔顶、底温差控制在6~8度。提高转化率，减少循环负荷，减少损失量。其次是控制好CO2气的纯度，我公司NHD脱碳工艺CO2纯度在98.5%以上，CO2气的氧含量能控制的是加空气量，正常生产时控制在指标的中下限，以减少带入系统的惰性气，减少尾气排放量。同时要稳定 合成塔负荷，避免合成塔压力温度的波动。入塔甲铵液温度严格控制在95℃以内，保证一甲泵长周期稳定运行。

3.控制好一、二段的分解温度

操作中控制好一、二段的分解温度在指标内，保证使一、二段的甲铵分解率和总氨蒸出率达到指标要求。即一段甲铵分解率达89%；二段甲铵分解率达98%，总氨蒸出率达99%以上。若一段分解不完全，则负荷带至二段，一是吸收用水增多，影响系统水平衡，二是尾气中氨含量升高，氨耗高。因为二段分解负荷增加，去蒸发尿液中氨含量上升。

4.吸收系统操作

4.1一吸塔操作是关键，一吸塔工况稳定与否是体现系统水平衡程度的主要标志。稳定一吸塔的操作，控制进水量，使一甲液浓度稳定，维持液位稳定，不排放。一吸塔尽量不用底部回流氨，少用顶部回流氨，塔底部温度控制不能太底，一般在90~95度；对于水溶液全循环法，按吨尿素为基准：氨泵2.2~2.4m3，一甲液1.1 m3，氨水0.16 m3，二甲液0.32 m3，一冷加水50kg，二冷加水100kg。在实际操作中，为提高热能利用段的效率和提高一吸塔下部鼓泡段吸收CO2 的能力，应尽量加大二甲液量，但总的补水量是不变的，但要注意在改变补水量的分配后，不要忘记一吸塔的最少喷淋量，防止过低的减氨水造成精洗段超温。从一吸塔上部进入精洗段的氨水量减少后，提高了该段氨水的浓度，所以该段的温度较好控制，顶氨的补入量也减少了，能够降低一定的氨耗和顶部温度易于控制；由于进入一吸塔下部 鼓泡段的二甲液增加，吸收一分气相中CO2 的能力增强，提高了一吸塔的生产能力，同时到精洗段气相中的CO2 相减少，顶部温度更加容易控制，所以能够达到降低氨耗、增强产能和稳定一吸塔操作的目的。脱盐水温度控制在90~95℃，使浓甲铵在一吸冷却器中不产生结晶的危险。水温调至95℃是上限，如在95~100℃下循环，因汽化而使泵的循环量下降，一吸塔底部温度会上升。

4.2操作中一、二冷液位要保持在下视镜中，防止一、二冷气相带液，同时要注意不能任意排放。

4.3控制好尾吸塔的补水量和塔温，减少放空气中的氨含量。4.4开好解吸，稳定好水解操作。

4.4.1解吸操作的依据：①解析塔气相中水含量尽量少；②解析废液中氨含量，尿素含量达标且尽量低；③达到前两个条件下尽量节约蒸汽。

4.4.2解吸目的：将氨水中的二氧化碳、氨分离出来，得到二氧化碳和氨的混合气体和水的过程，解吸是吸收的逆过程，解吸需降压和提温。为此，解吸压力越低越利于解吸，但解吸气要送至水解回流冷凝故需具有一定压力，克服阻力的影响，解析塔压力控制在0.30~0.43mpa为宜。

4.4.3解析塔上、下塔温度的控制。只要解吸压力确定了，就可以查出对应压力下水的沸点。在此室温下，即可达到解吸效果。解析塔上塔Ｔ≥135℃，解析塔出液Ｔ为145℃~153℃，增大蒸汽量不会改变塔中部和底部温度，而会使顶部温度升高，故应使底部蒸汽刚能使 中底部温度达到沸点为宜。

4.4.4解析塔出气温度的控制。为了保证解析塔气相中水含量尽量少应降低解吸气相温度，通过水解自身回流来控制，但同时要兼顾回流液浓度和解析压力，故解析气相温度不能降的过低，一般控制Ｔ≤121℃为宜。

4.4.5尽量提高碳铵液槽氨水浓度。进解析塔的氨水来自氨水槽、提氢等压回收塔氨水以及铜洗再生气回收氨水，其浓度受系统影响，不易选择，碳铵液槽氨水槽氨含量一般为5％~7％。浓度升高，会使进塔溶液中的氨分压增大，从而影响解析塔的正常操作，当氨水槽浓度升高时，操作要小心。当进解析塔的量及其浓度一定时，适当提高氨水温度，可以降低进解析塔的蒸汽用量，进液温度低，耗蒸汽多，进液温度高，易发生闪蒸，使解吸气相含水量增加。进解析塔氨水经解析换热器温度Ｔ＞90℃，一般为115℃左右。

4.4.6水解目的是回收氨水中的尿素。通过直接加热，产生水解反应。得到氨和二氧化碳。水解反应是在液相中产生的，故为使塔内液体不气化，故压力必须大于该温度下的饱和蒸汽压力。水解压力一般控制为1.5~1.7mpa。水解塔温度越高，水解越充分，但由于水解塔材质允许温度的限度。我公司使用的是316衬里水解塔，故水解塔下部壁温Ｔ≥199℃，出液温度Ｔ≥200℃。水解塔上部温度≥190℃，水解气相Ｔ≥180℃。水解塔液位控制为50％~90％。

5.蒸发系统操作

5.1闪蒸真空度维持在0.02~0.03Mpa，使尿液温度达90~95℃，能稳定真空送料至一蒸加热器底部。若闪蒸真空度过高，则下液不畅，蒸发操作不稳定。

5.2维持一、二段操作在工艺指标内，一、二段真空差保持在0.03 Mpa，否则通过二段加热器时间长，缩二尿增高。

一、二段分离器中不能存液。

5.3二表氨冷的操作，对提高二段真空度很重要，必须使液位维持在30%以下，气化后的气氨能释放出去，二表液温度在15~25℃左右。

5.4减少造粒过程中尿素的损失，根据负荷大小及时调整造粒喷头转速，减少造粒粉尘的产生和尿素粘塔现象。

6.压缩机操作

重点是提高压缩机的打气量，特别是在夏季，水温高，传热效果差，能耗高是制约生产的主要问题，因此如何有效的降低一段进气温度是提高压缩机打气量关键问题所在。

6.1提高一入气体压力，增设了一入压力调节阀，便于提高和稳定一入气体压力。

6.2加强循环冷却水管理，降低水温，保证水质；定期清理各段水冷器，提高冷却效果，降低压缩机功耗。

6.3定期清理汽缸夹套；定期清理活门，提高压缩机出力率。6.4进口气阀选择弹性适宜的弹簧，降低吸入阻力 6.5坚持每班按时检查各排油阀，防止泄漏。

6.6保证气质，提高CO2纯度，提高一入气体中的有效成分，减 少压缩机虚功。

四、结束语

通过以上的资金投入和技术改造后，经生产实践不断调整优化操作思路，装置运行效果较为理想，氨耗≤585kg/吨尿素，蒸汽≤1250kg/吨尿素，脱盐水消耗≤2.0m3/吨尿素，产品质量一级品率大于90%。实现了系统高产低耗、安全稳定长周期运行，达到了理想的经济效益。

第二篇：尿素改造的经验教训

2010年大修和尿素装置技改工作中吸取的教训

根据2010年技术改造项目规划，公司于2010.10.26安排了大修和550吨/年尿素技术改造项目。合成氨系统检修历时9天，尿素系统由于受节能减排的影响于12月4日原始开车。本次大修和技改合成氨系统共计实施了20项大修和小改小革项目，尿素系统技改涉及的内容较多，仅设备即更新了17台（套），除高压合成、压缩机及动力设备立足现有设备未作改动外，其余都做了改动。

大修技改结束投入生产运行后，合成氨系统虽然消除了安全隐患、实现了安全自动化，生产运行也比较平稳，但是原料煤消耗不但未达到预想指标，反而高于大修前。尿素系统先后经过8次大大小小的整改，目前仍存在尿素成品质量不稳定缩二脲超标、消耗高的问题。静下心来认真地进行思考和反思，我觉得主要有以下几点教训： 一是大修期间部分检修项目质量较差，开车后未达到预期效果。造气工段检修投入资金260多万元，新上了联合废锅、D500风机、彻底检修了余热锅炉，并且搞了一些小改小革，重新调整了煤气管线减少了阻力，但是煤气炉的发气量不但没有增加反而减少，原料煤消耗和蒸汽消耗不但没有降低反而生高。二次脱硫实施的一些改造项目总体来讲效果也不理想，被迫两次局部停车整改，至今脱硫效率仍未恢复到大修以前。脱碳岗位检修后同样也没有检修前效果好。之所以出现这种现象，暴露出我们部分管理干部心浮气躁、急于求成、对一些成功的技术采用后没有很好的消化吸收。原料煤消耗始终是我公司的

一大心病，为此我们在大修前多次外出考察学习、请专家会诊，并聘周生贤教授为我们的技术顾问。在这次大修中，对造气岗位实施了小阀加氮和内防流改造，由于在技术和操作上没有真正掌握要领，并且有些技术脱离我们公司的实际状况，导致原料煤消耗上升。净化车间在大修过程中把大量精力和时间放在小改小革上，忽视了正常的设备检修，小改小革没有达到预想的效果，检修质量也没有得到保证。二是尿素技改由于缺乏经验和受能力水平的限制，反复整改，造成重大损失。日产550吨尿素改造项目是公司2010年两大技术改造项目之一，从考察论证、方案确定、设备制作到安装开车，经历了近一年时间，并且整个过程还不断得到集团公司的技术支持，就是在这种情况下，在整个的试车过程中，先后8次停车整改，大大小小改动了21部位，影响生产长达两个月，使企业蒙受了重大损失。暴露出我们部分领导和干部本领偏弱、平平庸庸，技术人员盲目自信、缺乏严谨细致的科学态度。本次尿素改造原计划是由宁波远东实施交钥匙工程。但我们盲目自信，觉得没什么问题，在方案制定、流程及设备图纸设计等方面完全由公司技术部门自主完成，未通过专业设计部门进行计算和核算，在方案设计上欠科学、不系统，特别是具体细节的考虑上不够周全，造成反复的停车整改。例如：开车准备工作做的不好，吹除、试压没有严格按规程去做，导致开车后多次堵塞阀门管道：蒸发系统部分管线采用直角弯头，从配管的角度达到了横平竖直的美观但影响了尿液流速及过料；对新技术的应用缺乏考察论证，一吸塔设置多孔板不成功；一分塔塔板安装前调整泡罩高度不尽合理，两次进

行拆检；；一分加热器与预蒸馏塔之间的位差不够，造成预蒸馏塔内尿液不能全部入一分加热器，造成缩二脲超标等等。

三是作风飘浮、基础不扎实、不牢固，业务素质不强、缺少发现问题和解决问题的能力。生产系统部分管理干部，平时不爱学习、不下功夫钻研业务，不注重出主意想办法，提高解决问题的能力，关键时刻冲不上去。在造气的降耗攻关过程中，对问题吃不透、看不准采取措施不得力，职能处室与车间步调不统一，工艺调整多次出现反复造成炉矿不稳、消耗上生。

四是高素质的工程技术人员和管理干部匮乏、一线生产骨干年龄偏大。由于前几年公司生产经营困难、效益较差，员工的工资福利较低，一些有技术、懂管理、会管理的人才大量流失。企业人员虽然不少，但很大一部分是七十、八十年代的老工人，这部分人年龄大、学历低，接受新知识的意识差，而专业对口的学生和年轻人又不愿来公司，造成目前生产一线骨干青黄不接，出现了严重的断层。企业维持现状人员上已感到不能满足，如要加快发展新上项目和技术改造、实施新技术、新工艺、新设备、新材料则更是力不从心。这也是这次大修和尿素技改失误较多的重要原因。

整改措施

一是统一思想，克服畏难发愁、悲观失望的情绪，坚定信心、攻坚克难，对生产中出现的这些问题找准原因，对症下药。

二是推进“人才强企”战略，对内完善动态管理，对外积极引进人才，全力推进学习型企业进程。

三是完善机制、理顺职能、健全制度、严格考核。四是多措并举进一步节能降耗和降本增效。2011年企业的总体发展思路及目标

2011年是“十二五”开局之年，也是公司在2010年打下的良好基础上应该加快发展的关键一年。总体思路是：以国家产业政策为指导，以市场为导向，以提高经济效益为中心，突出“转方式、调结钢”的主线，通过改造提升现有合成氨、尿素装置，研制开发橡胶助剂产品，整合资源、拉伸产业链条发展下游产品，实现结构优化调整，真抓实干，务实创新，全力推进企业又好又快发展。奋斗目标：实现销售收入4.3亿元，利润680万元,氨醇11万吨，尿素15万吨，碳酸氢铵1.8万吨,液体CO2 1.84万吨，工业氨水1.1万吨；助剂产品RT 2500吨、4020 3700吨、FAO-7 300吨、其它100吨。为此，将主要抓好以下工作：

（1）、立足主业做优化肥，做到化肥“保吃饭”。依托原料煤的优势从节能降耗，降低成本上下功夫。

化肥是企业做了40多年的产品，其他后续产品都是依托化肥而做的，这是我们的根，是我们的立足之本。作为尿素这个产品，质量没有多大的差别，品牌效应也不是很明显，拼的就是价格，看谁的成本低。2011年我们将依托原料煤的优势从节能降耗，降低成本上下功夫。对现有装置填平补齐，挖潜改造，消化吸收行业内最先进的技术，不断完善合成氨和尿素产品的低成本竞争的各项措施，将现有的装置能力发挥到极致，吨合成氨煤耗1.16吨、综合电耗1480KWH；尿素

单耗液氨580Kg、电140KWH、蒸汽1200Kg，尿素成本下降100元，实现真正意义的“做优化肥”目标。

（2）、借助我们的优势和经验，调整产品结构，完成2万吨/年橡胶助剂产品。

我们的橡胶助剂产品品种比较多，但规模比较小，主要是受中间体的制约，通过与南京理工大学合作，依托我们的原料优势，开发橡胶助剂的中间体，将橡胶助剂系列产品做到20000吨/年。

（3）、整合资源，优势互补，拉伸产业链条，形成上下游结合的循环经济。

循环经济是今后经济社会发展的方向，是化工行业发展的必然趋势，而与公司毗邻的几个厂家基本上形成了典型的“减量化、再利用和资源化”的循环经济产业链条模式。硝酸是我们合成氨产品的下游产品，在一个厂区内合成氨装置和硝酸装置属于两个企业，资源浪费巨大。液氨和硝酸都属危化品，给运输和公路交通增大压力，又浪费大量能源。通过整合全部由管道输送，靠两装置之间工艺存在的压力差和反应热使资源达到最优化的配置，大大降低硝酸成本，提高竞争力。降解塑料是二氧化碳的下游产品，将降解塑料装置，合成氨装置和二氧化碳装置三套装置整合在一起，相互之间形成完整的产业链，装置之间相互依托，一套装置的排放废气是另一套装置的原料，使资源达到最充分的利用，既节约了资源减少了排放，又大幅度的降低了成本。一是收购华阳迪尔现有4万吨/年浓硝酸装置，新增销售收入8000万元，利润500万元。二是探讨整合金明热电，新增280吨供

热，4.5万Kw的发电能力，实现互补双赢。三是进一步考察论证10000吨/年降解塑料项目。

（4）、依托基础，走边缘化战略。

在氮肥行业我们的装置能力小，工艺、设备落后，怎么赶都赶不上那些单套能力上百万吨的新装置。在橡胶助剂行业，我们也没法与圣奥等大企业相比，但是他们这些大企业都是一门心思的发展单一产品，做成行业老大，全国老大。有很多相关的小产品他们根本不屑一顾，这些产品被他们边缘化。我们下一步就是搞那些被大企业边缘化的产品。

张大涛 2011年1月30日

第三篇：空调节能优化改造项目施工方案

安徽移动西园枢纽楼节能优化改造项目施

工方案

1、工程概况：

该项目位于合肥市长江西路99号安徽移动西园枢纽楼，由于空调安装较早，使用年限较久，系统老化，加之机房设备增加，现有的水冷螺杆机系统及风冷热泵系统已经不能满足该机房楼正常的生产需要；且现有系统故障率升高，可靠性降低，机组能耗增加。严重影响通信设备的安全运行，需要通过对空调系统进行节能及优化改造，以降低机房空调能耗，提高系统运行安全性及稳定性。

由于本次为改造项目施工，机房内设备已投入使用，因此在施工过程中，需在做好安全防护措施，保证机房内设备正常运行，降低施工过程中的风险。

2、工程整改内容：

1)拆除现有的风冷热泵系统及；

2)完成对梅州第二机楼二、三楼气流组织进行改造，加强局部回风工作；

3)完成对梅州第二机楼冷却水系统水管更换、冷却水泵更换、配电系统的整改工作； 4)配合设计单位完成节能优化改造后的效果分析及机房PUE值评估报告。

主要工作内容对原有空调系统的管路、设备及配电进行改造：新建环状供回水主干管，新增2台热回收离心式冷水机组，与原有螺杆机一并为现有3套系统提供冷冻水；原有热泵系统、集分水器、泵组、无关管路及无关配电设备等均拆除；新增一套主干立管，与原有立管同时为2-11层机房内末端空调提供冷冻水，形成双立管供水，并设置末端切换机制；12-15层新建机房区域，由新建立管和原有高区热泵立管同时为末端空调提供冷冻水，形成双立管供水；16-19层舒适性空调区域，夏季由本系统直接供冷，冬季由热回收离心机提供35/40℃-40/45℃热水供热，紧急情况时，该立管可通过阀门切换回制冷状态，为12-15层机房供冷。

3、施工风险评估分析及影响：

由于机房为在用机房，消防、空调、监控、通信等全部设备已投入运行使用，在施工过程中需要做好安全保护措施。现三楼有3套全年供冷的3组360RT的风冷热泵系统（2用1备），此次需拆除，由于设备老化，吊装时应做好防护加固措施；新增一套主干立管与原有螺杆机系统立管并联形成双立管供水，并设置末端切换装置，需要在各楼层支管进行割接工作；本次的主要风险影响有：机房空调冷源、机房内防尘、防水、防火安全等影响。1）、本次施工需要对新增一套主干立管与原有螺杆机系统形成双立管供水，需要在各楼层支管进行管道割接，并安装电动二通阀进行管路切换系统内冷水需放空，机房缺少后备空调冷源保障，故将管道割接施工安排在冬季最冷的一月份进行，并增加施工班组在各楼层同时进行割接工作，在保证施工质量的前提下，尽可能减少施工时间，尽快恢复机房空调供冷。施工期间需要打开机房内窗户进行通风，并加设临时工业风扇，增强室内外热冷空气的对流交换，以保障机房设备的安全运行；必要时，在机房内安装分体空调，以保障冷源。

2）、各楼层墙面开孔工作，采用专用墙体水钻并在施工现场布置防水帆布、集水盆和吸尘器，以免施工时产生的水滴、尘埃飘落到机房设备上，造成设备损坏。

3）、供水管道的割接工作：焊工应持证上岗，作业时需取得安全部门开具的动火证，并严格检查操作现场和作业点下方10m内不得存有易燃易爆和杂物。必须将焊接点周围10m内的易燃易爆物排除或采取防火卷帘隔离后，方可进行操作。焊接现场同时必须配备足够数量的灭火器。

4、施工步骤：

第一阶段：施工前准备工作：合同签订、设备订货、施工员人员机房施工安全知识培训、组织图纸会审；

第二阶段：施工前做好机房温度测试记录，并安装100A新装空调配电箱及进行调试，安装分体空调并进行调试试运行，开启新装及原有分体空调；

第三阶段：停二楼水冷式柜机空调系统，观察分体空调是否满足供冷需求，如满足则开始改造施工二楼水冷式柜机空调系统及一楼天面冷却水系统。如不满足供冷需求，则增加安装分体空调；

第四阶段：停三楼水冷式柜机空调系统，观察分体空调是否满足供冷需求，如满足则开始改造施工三楼水冷式柜机空调系统。如不满足供冷需求，则增加安装分体空调；

第五阶段：更换天面配电柜、敷设不锈钢线槽及电缆，并割接安装至冷却水塔及水泵中； 第六阶段：

二、三楼新装水冷式柜机空调系统投入试运行阶段；

第七阶段：测试改造后效果并进行分析总结，并进行资料归档及竣工验收。

5、具体施工步骤：

1)前期施工准备：

2)冷冻站施工前整理：拆除隔墙、平整地坪、做离心机、水泵基础等； 3)新增离心机、水泵搬运至地下室；

4)地下室DN400环网施工，新建DN400立管施工，新立管与各层末端沟通管道施工； 5)2、8层AHU转入高区热泵系统工作；

低区热泵系统（含主机、水泵、无关管道阀门）拆除，以下设备就位：

软化水箱搬迁

冷却塔就位

冷水机组就位，冷却水泵、冷冻水泵；

定压补水设备；

6)低区热泵系统立管与环网沟通；

7)新建管道系统试压（支管逐层试压，做好安全防护措施）；

自控系统同步实施；

8)各层割接至新立管，试运行1周时间。新系统自控系统调试。

第一阶段系统水力平衡；

9)关停螺杆机组；螺杆机与冷却塔、冷却水泵深度保养；拆除原有螺杆机冷冻水泵组； 10)螺杆机冷冻水接入环网；

11)高区热泵出水管接入环网；高区立管接入环网；热回收管道接入高区立管； 12)螺杆机接入自控系统；第二阶段系统水力平衡； 13)逐层拆除原有AHU；

14)各膨胀水箱拆除；无关管道拆除封堵； 15)自控系统调试；工程验收； 16)

 施工进场前期准备工作：合同签订、安装材料、辅材采购、施工人员通信机房安全知识培训、机房注意事项培训、技术交底、与业主单位、设计单位进行方案会审等；  施工人员现场勘察并办理人员出入机楼证件、动火证件等；

 在各层机房内做好温度测试记录，在加设的工业风扇或应急用分体空调开启前，做好机房的温度测试记录，在机房内选定多处测试点，对机房进行温度测试并记录数据，在系统停机，管路冷水放空后，做好机房温度对比，保障机房供冷安全； 17)安装临时工业风扇

 在各楼层新增安装可移动工业风扇10台，促使机房内气流组织流动顺畅，以增强机房内热空气与室外冷空气的对流，达到降低机房温度的效果，防止出现局部高温现象；电源采用机房内现有电源，电源线的摆放布设要求紧贴墙角，并加装线管等保护措施，在醒目位置放置指示牌，避免人员误操作或不小心而发生触电等意外伤害。18)停2-11楼螺杆机空调系统：

 为减少螺杆机系统停用的时间，保障机房设备安全运行，我司将采用各楼层管路同时进行割接施工的方法进行该系统的改造施工，减少施工风险。先将该系统停机，管路中的冷水放空，然后在各楼层布置一个班组进行管路的割接工作。最大限度的减少施工时间，减少系统停机时间，保障机房冷源及时供应。

 冷水管割接完成后进行试压，保压测试合格后，进行下一步工作；  墙体开孔防护措施：

在进行墙体开回风口的过程中，在此施工过程中需要特别做好安全防护措施，在墙体开孔过程中，通过制作防尘罩的方法把开孔处遮盖住（防尘罩规格：1200×1800mm），防尘罩周边粘贴橡胶密封条，这样防尘罩能够严实的压在开孔处，避免灰尘流进机房设备，开孔时使用墙体专用水钻开孔，减少开孔的灰尘，在施工过程中在机房内使用吸尘器把少量灰尘及时吸走，同时使用集水盆将水钻流出的水收集起来，避免流入机房设备，确保不影响设备运行使用。 管道割接防护措施：

原有风管上加装风管及风口，为防止安装过种中有铁屑掉落到通信设备上，需在施工位置放置接铁屑盆以及吸尘器，同时使用专用的镀锌板开孔工具，防止铁屑飞溅，控制铁屑掉落范围，在施工区域周围遮盖防火布。

 对一、二楼新装水冷式柜机空调系统进行调试，调试合格后对空调系统进行试运行，在确保一、二楼新装水冷式柜机空调系统运行正常，满足空调机房冷量要求后，准备开始三楼水冷式柜机空调系统的改造工作，同时需密切观察动环监控中心； 19)停三楼水冷式柜机空调系统

 停三楼水冷式柜机空调系统，拆除天面冷却水管及水泵，并对天面水平冷却水管、冷却水泵及相应的阀门进行更换，对冷却水管加装保护层，水管用不锈钢进行包裹防止腐蚀；  进行三楼空调系统冷却塔风机轴承、皮带、防护网及接水盘补漏翻新工作；

 拆除三楼空调机房的冷却水管及阀门，逐台拆除空调风柜并更换新的空调风柜，在空调机房新增回风口、风阀、风口软连接以及防火阀，同时在三楼加装新风管及开风口（三楼机房墙体回风口开孔及风管开风口施工及防护方法与二楼一致）；  敷设新装从空调电箱到水冷式空调机的电源线路，电缆规格：ZR-VV1KV3*16+1*10mm2-BVV500V-10mm2  冷水管更换完成后进行试压，保压测试合格后，进行下一步工作；  安装吊顶式离心通风机、风管以及电源线路等；

 对三楼新装水冷式柜机空调系统进行调试，调试合格后对空调系统进行试运行，在确保三楼新装水冷式柜机空调系统运行正常，满足空调机房冷量要求后，开始天面空调配电柜割接更换工作。

 为提高机房的空调冷源保障，我司在更换天面配电柜时以不长时间停二、三楼新装风柜的形式，进行更换天配电柜，把原有5AC1、5AC2空调配电柜冷却水塔及水泵负载割接至临时配电柜中，由临时配电柜进行供电，确认临时配电柜供电稳定及新装水冷式柜机空调运行稳定后，拆除原有配电柜并安装新的空调配电柜；

 对天面新装的空调配电柜进行调试，敷设新装不锈钢线槽及相应电缆，并把新敷设的电缆逐台割接至冷却塔及相应的冷却水泵中；

 停二、三楼水冷式柜机空调系统（此时分体空调进行供冷），把二、三楼空调冷却水系统在天面阀门处进行连通，此时5台冷却塔可以相互互为备用。

 二、三楼空调系统改造完成，重新启动水冷式柜机空调系统正式投入运行阶段。 停用二、三楼分体空调，并对二、三楼机房做温度测试记录，测试改造后效果并进行分析总结。

 施工完成，资料归档及竣工验收。

6、应急保障方案：

本次在用机房改造施工过程中，需做好安全保障措施后方可进行施工，通过各程方法降低施工风险，把施工风险降到最低限度。方案一：

在更换空调风柜施工过程中，如遇到部分保障分体空调出现故障，及时起动大风量工业风扇，马上把分体空调故障消除解决（在施工过程中，有专人空调维修工程师驻现场，及时处理故障）； 方案二：

在更换空调风柜施工过程中，如出现大面积分体空调停机不能工作，马上起动大风量工业风扇，在走廊放置桶装冰块，能够起到降温作用，打开所有机房的门，起动消防排风机，利用大楼自带的消防排风机把机房内的热量排走，同时机房内形式负压，机房门打开后，室外的空气流到机房内可起到降温作用（由于现在处于冬季，室外的温度约为16~20℃，温度较低，把室外空气导流进机房，能够起到很好的降温作用。）马上增补分体空调，工程师维修空调消除故障，待空调消除故障后，重新起动空调，把机房关好以及停止消防排风机工作。方案三：

在更换空调风柜施工过程中，准备移动油机车到现场，确保供电保障。

7、施工难点、重点技术措施：

本次施工的为在用机房的改造项目，需要考虑到机房通信设备的使用安全，施工过程存在危险通信设备安全运行的因素，需要分析难点、重点把施工风险因素降到最低。难点、重点因素一：

在通信设备上方安装风管及开风口，安装过程中有铁屑掉落的危险源存在，因此必须要危险源去除掉，做好安全保护措施；

1)在通信设备上方施工时，旁边必须要有监护人员，必须提醒施工人员注意安全，以及要做好防护措施；

2)安全防护措施做好后，需要由监护人员进行检查并报监理、业主人员审批后方可进行施工，做到严格的操作流程，防止意外发生。

3)在通信设备上方安装风管及开风口时，通知通信业务对口人员做好数据备份工作； 4)在通信设备上方敷设防火布（注意：防火布在敷设时不能阻挡通信设备的散热），根据现场情况专门制作一个可以接尘、接铁屑的盒子，在开孔时掉落的少量灰尘及铁屑不会掉到设备上；

5)使用镀锌板专用开孔工具，该工具开孔时不会产生火花，掉落的铁屑也不会到处飞溅，能够提高施工的安全系统数，做到零风险！难点、重点因素二：

需要在投入使用的机房内进行土建开孔，必须保证机房内环境卫生、无尘，不能让灰尘影响设备的使用安全。

1)在机房内进行土建开孔时，我司预先制作好一个防尘罩（防尘罩规格：1200×1800mm），把防尘罩安装在开孔的位置，防尘罩四周紧贴着墙体，防尘罩安装在通信机房一侧，施工有空调房内进行开孔；

2)施工人员通过水钻专用开孔工具在空调机房一侧进行开孔作业，在水钻头处安装一个接水罩，并在下方放置一个接水盘，把水钻少量的水接住，这样可以做到无尘化施工作业。保证机房通信设备使用安全。如图所示：

密闭防尘罩通信机房空调机房施工人员要空调机房内开孔

难点、重点因素三：

由于本次工程需要更换水管，因此需要空调机房内进行动火作业，为减少动火次时及时间，以减少危险因素源，我司采用的施工方法是：

1)风柜就位安装后，施工人员在空调机房内测量水管的尺寸长度，然后把水管拿到室外进行焊接加工，做到每台风柜水管只预留一个焊接口，以达到减少动火作业次数目的，降施工风险源；

2)在施工焊接过程中，施工人员必须在旁边放置手提式灭火筒，有专人进行监护施工，提高安全施工性；

3)在焊接过程中产生的烟雾通过小型抽风机向空调机房的新风口排走，同时通知消防值班人员把附近的烟感进行屏蔽，以免消防产生误报警； 难点、重点因素四：

在通信设备上方楼板打冲击钻时，打钻产生的灰尘及沙子会掉落到设备上，会影响设备正常使用。

1)在进行打楼板膨胀螺丝孔时使用专用的自带吸尘式冲击钻，该工具在打孔过程中，自带的吸尘会工作起来，起到很好的吸尘效果；

2)在打孔时，再放置一个外置的吸尘器，把吸尘器对准打孔的位置及时把灰尘吸走，同时设备上方放置一块防火布进行遮盖，每打一个孔清理一次防火布，做无灰化作业施工。

第四篇：优化小结

优化第一步：

确定分子构型，可以根据对分子的了解通过GVIEW和CHEM3D等软件来构建，但更多是通过实验数据来构建（如根据晶体软件获得高斯直角坐标输入文件，软件可在大话西游上下载，用GVIEW可生成Z-矩阵高斯输入文件），需要注意的是分子的原子的序号是由输入原子的顺序或构建原子的顺序决定来实现的，所以为实现对称性输入，一定要保证第一个输入的原子是对称中心，这样可以提高运算速度。我算的分子比较大，一直未曾尝试过，希望作过这方面工作的朋友能补全它。以下是从本论坛，大话西游及宏剑公司上下载的帖子。将键长相近的，如

B12 1.08589 B13 1.08581 B14 1.08544 键角相近的，如

A6 119.66589 A7 120.46585 A8 119.36016 二面角相近的 如

D10-179.82816 D11-179.71092 都改为一致，听说这样可以减少变量，提高计算效率，是吗？ 在第一步和在以后取某些键长键角相等，感觉是一样的。

只是在第一步就设为相等，除非有实验上的证据，不然就是纯粹的凭经验了。在前面计算的基础上，如果你比较信赖前面的计算，那么设为相等，倒还有些依据。但是，设为相等，总是冒些风险的。对于没有对称性的体系，应该是没有绝对的相等的。

或许可以这么试试： 先PM3，再B3LYP/6-31G.（其中的某些键长键角设为相等），再B3LYP/6-31G（放开人为设定的那些键长键角相等的约束）。

比如键长，键角，还有是否成键的问题，Gview看起来就是不精确，不过基本上没问题，要是限制它们也许就有很大的问题，能量上一般会有差异，有时还比较大

如果要减少优化参数，不是仅仅将相似的参数改为一致，而是要根据对称性，采用相同 的参数。例如对苯分子分子指定部分如下：

C

C 1 B1

C 2 B2 1 A1

C 3 B3 2 A2 1 D1

C 4 B4 3 A3 2 D2

C 1 B5 2 A4 3 D3

H 1 B6 2 A5 3 D4

H 2 B7 1 A6 6 D5

H 3 B8 2 A7 1 D6

H 4 B9 3 A8 2 D7

H 5 B10 4 A9 3 D8

H 6 B11 1 A10 2 D9

B1 1.395160

B2 1.394712

B3 1.395427

B4 1.394825

B5 1.394829

B6 1.099610

B7 1.099655

B8 1.099680

B9 1.099680

B10 1.099761

B11 1.099604

A1 120.008632

A2 119.994165

A3 119.993992

A4 119.998457

A5 119.997223

A6 119.980770

A7 120.012795

A8 119.981142

A9 120.011343

A10 120.007997

D1-0.056843

D2 0.034114

D3 0.032348

D4-179.972926

D5 179.953248

D6 179.961852

D7-179.996436

D8-179.999514

D9 179.989175

参数很多，但是通过对称性原则，并且采用亚原子可以将参数减少为：

X

X 1 B0

C 1 B1 2 A1

C 1 B1 2 A1 3 D1

C 1 B1 2 A1 4 D1

C 1 B1 2 A1 5 D1

C 1 B1 2 A1 6 D1

C 1 B1 2 A1 7 D1

H 1 B2 2 A1 8 D1

H 1 B2 2 A1 3 D1

H 1 B2 2 A1 4 D1

H 1 B2 2 A1 5 D1

H 1 B2 2 A1 6 D1

H 1 B2 2 A1 7 D1

B0 1.0

B1 1.2

B2 2.2

A1 90.0

D1 60.0

对于这两个工作，所用的时间为57s和36s，对称性为C01和D6H，明显后者要远远优于前者。

好，我真的很想知道这句话“你手工把坐标写好输入，输入为C2对称性”是什么意思？老是见到版面上有人说对输入施加对称性限制，就是不明白是怎么回事。您能详细解释一下吗？有例子就更好了。麻烦了谢谢！

比如，你要输入CO2，那么你这么输入 C 0.0 0.0 0.0 O 0.0 0.0 1.203212313123 O 0.0 0.0 1.203212313123 写得很精确，那就是对称性输入了。这样g98可自动判定为D∞h 如果你输入的差不多，比如： C 0.0 0.0 0.0 O 0.0 0.0 1.203212313123 O 0.0 0.0 1.204212313123 精度不够，那么g98就只能判定为 C∞h了。

再给你举个例子，比如甲烷的Td点群的输入坐标，考虑对称性，把C放在直角坐标的原点，四个氢放在四个角的顶点。C 0.0.0.H 1.1.1.H-1.-1.1.H 1.-1.-1.H-1.-1 1.这样出来的结果点群判断就会是正确的，当然这里的1我是随便取的一个值，只是来说明这个问题而矣。

谢谢二位指教，也就是说第一个输入的原子应该是对称中心，其它输入的原子应保持数值上的对应性。那对于用Gview或3D生成的分子构型，也可以根据这个原则将处于对称位置的原子的坐标改为一致了？如果分子很大（如我这次算的C14H14N4S20）有三个环，整体性的对称输入就有些麻烦，是不是可以先生成分子构型，然后利用对称性对它的笛卡儿坐标或矩阵坐标进行修改，将相近的坐标的数据改为一致？

“第一个输入的原子应该是对称中心，其它输入的原子应保持数值上的对应性”没错,做到这点就很方便了,“对于用Gview或3D生成的分子构型，也可以根据这个原则将处于对称位置的原子的坐标改为一致”那还不如用手动输入,大体系很难办了,如果体系太大,不用保持严格的对称性吧.有机分子，原子多了以后构型就没有什么对称性了。

但对于有些系统，比如 C60, 你当然得按照它的实际的对称性来输入。如果你希望你的计算按照某种对称性来做，那就一定要手工输入严格的 满足对称性的坐标.用软件生成的坐标，不是严格满足对称性，即使 用symm=loose,有时也不一定能判断出来。

各位大虾是用什么软件建模的呀，我是用SGI里sybyl软件进行构像搜索后取最低能量构像导入高斯程序的；有的是用chem3D直接建模，这样子有没有什么问题啊，应该怎么做，谢谢

通常也就是这些方法,建议还是认真的核实分子结构的是否合理,然后再试着增大迭代圈数,降低收敛判据.用什么软件建模并不是问题，对于我们量化计算，如vliant所说，核实分子的合理构型才是最重要的，包 括对个原子间的成键分析。就算不同的软件，也不可能保证给出的构型合理，因此我们一般都提倡手动输入坐标，就是为了加强对构型的理解，只有在这个基础上使用软件建模才能尽可能保证不出错。一个合理的计算结果不能寄希望于软件。

我们都知道常用的分子坐标输入有： 内坐标 直角坐标 混合坐标

分子构型输入方法的不同会对计算过程产生影响么？或者说，对于哪些体系，用哪种方法更好优化？ “混合坐标”怎么讲？

我通常用内坐标，它能很好的控制对称性，且时间要短。

“混合坐标”是直角坐标和内坐标的混合使用，高斯手册上有讲这部分内容。有看过讲哪种体系使用哪种方法更好优化，好像这都是凭经验的结果。

在“免费资源”论坛中有本《计算化学》，其中第八章（67－72页）好像提到你所关心的问题，你可以下来看一下，或许有所帮助。

如果你使用缺省的优化方案且初始几何构型完全相同，那么，不论你使用何种坐标输入方式，优化结果应该是一样的。简言之，优化结果只取决于优化方案和初始几何构型，与坐标录入方式无甚关联。

于对称性限制的优化来说，显然以内坐标输入，且选用opt=z-matrix可以减少变量，使得优化效率提高。免费资源中《gauss入门》有介绍。

您冗余坐标是指有过剩的内坐标。冗余内坐标算法可以用于结构优化。

优化第二步：

对于大分子，一般要先预优化（如用PM3或AM1），有时是用部分优化，然后才用能达到实验精度方法和基组进行优化，就我所了解大多数情况下，6-31G（D）基组就差不多。

不好意思，我是新手，刚刚接触guassian，最近要算两个较大有机分子中原子的静电荷急相关能量和几何性质，分子中已经超过了80个原子（包括氢），用chemdraw画出图形，在chem3d中生成的input文件，可是，软件报错——2070~~ 好象说是连接（last link）错误，怎么回事，哪里出错了，各位指点一下~ 怎么，先用pm3预优化就可以了？？ 我的意思是，初试构象对计算量有这么大的影响么？？

我以前也遇到过相似情况，同一个分子，相同的计算条件，用gview画出的模型进行计算报2070的错，用chemdraw画的模型，chem3d生成的input文件进行计算就OK，完成可以计算完成。我还以为你画分子式的时候出什么问题了呢~~

只想优化H原子，其他重原子不做优化。怎样进行分子坐标输入，下面的坐标输入方法对吗？就是仅H原子的参数设置成Vairable.%chk=sample #p HF/6-31G* opt

Title line 0 1 cu 7.577854 5.926989 4.935050 cu 9.091644 8.489386 4.644861 cu 10.605857 5.766733 4.988313 p 10.432217 3.600341 5.598077 p 12.551342 6.888289 4.775263 p 7.164457 9.584280 4.209578 p 11.125158 9.457819 4.661391 p 7.503791 3.715690 5.278502 p 5.670551 7.057133 4.652208 s 9.100712 6.979012 6.564150 s 9.088921 6.492596 3.113104 n 8.884097 3.048614 5.950712 n 5.827205 8.602518 4.003874 n 12.406788 8.511755 5.225239 H 1 B14 2 A13 7 D12 H 1 B15 2 A14 7 D13 H 1 B16 2 A15 7 D14 H 1 B17 2 A16 7 D15 H 2 B18 1 A17 3 D16 H 2 B19 1 A18 3 D17 H 2 B20 1 A19 3 D18 H 2 B21 1 A20 3 D19 H 3 B22 1 A21 2 D20 H 3 B23 1 A22 2 D21 H 3 B24 1 A23 2 D22 H 3 B25 1 A24 2 D23 H 1 B26 2 A25 7 D24 H 2 B27 1 A26 9 D25 B14 2.990841 B15 2.990841 B16 2.985833 B17 2.985833 B18 3.009872 B19 3.009872 B20 2.998051 B21 2.998051 B22 3.015712 B23 3.015712 B24 3.002042 B25 3.002042 B26 3.489189 B27 3.477113 A13 161.290094 A14 146.871331 A15 108.517371 A16 96.114237 A17 106.549047 A18 100.057030 A19 155.415250 A20 145.871488 A21 163.425916 A22 138.991013 A23 105.076579 A24 99.387878 A25 66.120736 A26 65.475681 D12 64.540783 D13-53.098055 D14-163.347541 D15 152.573174 D16 152.164746 D17-162.671553 D18-56.766147 D19 42.781046 D20 46.122424 D21-46.534404 D22 150.210575 D23-164.670929 D24-78.866142 D25 109.837455 你如果要这么做的话，必须加opt＝Z-matrix,这样的话它就只优化后面的H，前面的直角坐标就不会优化了。你只用了opt的话，那就是全优化了 你也可以通过molredunt来实现部分优化 固定原子的方法有许多,最简单的是这样的： #rhf/sto-3g opt nosymm test the two C atoms are frozen 0,1 8 0-1.000000 0.000000 0.000000 6-1 0.000000 1.000000 0.000000 6-1 0.000000-1.000000 0.000000 8 0 1.000000 0.000000 0.000000 我不知道你的那种方法行不行，但我这个肯定可以。前面我都贴出来了：）这可是我的心血噢。＝＝＝＝＝＝＝＝＝ opt=z-matrix的部分优化。%Chk=reac-hf3 #p RHF/STO-3G opt(z-matrix,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 x2 y2 z2 H 0 x3 y3 z3 H 0 x4 y4 z4 O-1 2.16296 0.18089-0.15503 H 0 x6 y6 z6 H 0 x7 y7 z7 H 0 x8 y8 z8 x2=-2.2088 y2=0.39705 z2=0.46255 x3=-2.14589 y3=-0.9279 z3=-0.4899 x4=-1.88451 y4=-1.06386 z4=1.11629 x6=1.36256 y6=0.28882 z6=-0.67401 x7=2.32973 y7=0.98241 z7=0.34634 x8=2.90378 y8=0.01009 z8=-0.74115 *********************************************** opt=modredundant的部分优化。

分子坐标部分可以为具体的数值，而不用为参数形式。%Chk=reac-hf3 #p RHF/STO-3G opt(modredundant,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 H 0 H 0 O-1 2.16296 0.18089-0.15503 H 0 H 0 H 0 1 5 F

感谢各位大侠的帮助！

看了之后有个不明白的地方，还请大侠不吝赐教！在opt=modredundant的部分优化时最后一行的 1 5 F 是什么意思？

就是将第一和第五号原子的距离固定

问一个菜鸟的问题：N 后面的-1和H后面的0是什么意思？是像手册上所说的把坐标增加所需要的值吗？ 见笑了！

－1指的是固定，0指的是优化，前面helpme贴的输入方法也是这样的

部分优化输入格式之二

opt=z-matrix的部分优化。%Chk=reac-hf3 #p RHF/STO-3G opt(z-matrix,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 x2 y2 z2 H 0 x3 y3 z3 H 0 x4 y4 z4 O-1 2.16296 0.18089-0.15503 H 0 x6 y6 z6 H 0 x7 y7 z7 H 0 x8 y8 z8 x2=-2.2088 y2=0.39705 z2=0.46255 x3=-2.14589 y3=-0.9279 z3=-0.4899 x4=-1.88451 y4=-1.06386 z4=1.11629 x6=1.36256 y6=0.28882 z6=-0.67401 x7=2.32973 y7=0.98241 z7=0.34634 x8=2.90378 y8=0.01009 z8=-0.74115 *********************************************** opt=modredundant的部分优化。

分子坐标部分可以为具体的数值，而不用为参数形式。%Chk=reac-hf3 #p RHF/STO-3G opt(modredundant,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 H 0 H 0 O-1 2.16296 0.18089-0.15503 H 0 H 0 H 0 1 5 F

冻结坐标的目的是什么了?我虽然会做但不太清楚Gaussian为什么会有这个功能.还有最近做了一个体系.有一个键如果不冻结的话.就优化不出那个反应中间体.那个键会断.我这样冻结键优化出来的结构计算频率也无虚频.它能看作为一个minima吗?这个结构可以作为我反应机理的中间体吗?我用溶剂化优化结构可以在不冻结坐标的情况得到一个和这个结构很接近的结构.冻结坐标有时是有用的，举几个例子：

1，如优化溶质分子和nH2O的supermolecule，周围的水分子，可能结构变化不大，这时就可以固定水的键长或键角，减少工作量。

2，如用cluster模拟大块固体，这时也许就要按照晶格参数固定原子键长和键角。3，又如不固定NH3－H－H2O中N－O键长是得不到质子传递势垒。

4，还有在优化稳定构型或过渡态得时候也有时用到部分优化得到接近得结果，再放开优化，这样逐步得到你想要得结果。

冻结坐标对一些大的体系,超分子体系好象是很有用的.好象做生物分子的反应时氨基酸残基必须先固定.但我这个体系如果不固定,只要一放开,结构就跑了会优化到其它构型上去.就算stepsize设得很小也没有用.优化第三步：

优化结果分析，一般要进行频率计算，无虚频，表示是极小值；有虚频，表示是过渡态。如果出现虚频，可根据虚频的简正模式进行调整以获得极小值，我还没有碰到这种情况，希望各位朋友能给些详细的资料。

如前信所说，键长键角有差别是很正常的，因为气相分子和晶体内分子的结构不可能完 全吻合，只可能近似相同。

另外，对不同元素的限制主要看基组是否能计算该类元素，例如你所采用的方法中用了

6－31G(D)基组，该基组支持H-Kr的元素，即可以计算Br原子但是不能够计算I原子。如果你所

用的基组不支持你所计算的某个原子可以采用扩展基组来进行计算，或者降低基组的精度。

我用Gaussian优化一个分子结构，由于目前实验上没有办法得到纯净的样品来实验验证，那么我怎么判断优化得到的结构是正确的？

一些主要的键长和键角的实验值或者别人计算的值是否和与你计算的分子类似？ 这个判断不好做。

所谓正确不正确，如果计算出来的频率没有虚频，就表示你找到了一个局域稳定点。要确定是不是全局稳定点很难，只有通过别的方式比如HOMO-LUMO能级差等信息看你算出来态是不是稳定。

也觉得高级的方法和大基组的结果更可信,谢谢各位高手的指点！再问一个问题,如果出现虚频如何修改得到分子稳定的结构

用GV看一下，根椐产生虚频的振动模式调整结构

消去虚频有什么技巧？我有一个计算结果，虚频有三个，都只有－30~－20左右，如果结构改动过大，一不小心，虚频就会更大。头痛ing～～～

那你就把结构朝振动的方向改小一点试一下 你

确实据我了解，消除虚频的主要方法是改变构型。

其次在计算上还可以尝试：nosymm；加大循环次数；提高收敛度；iop(1/8=1)等。其实我自己也是有这个问题，对于势能面很平，较小的虚频很难消除。首先，1在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果 还有虚频，参见下一步。

2.对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。

3.如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中 产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到 虚频没有。实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的 东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做 一般的电子结构计算，这样，应当没有问题。

不要讲你还没有得到最稳定的结构，那么，是你的分子有问题，要么计算错了，要么就是游离出现代计算的范畴 zixia上大虾的回复: 将内坐标改成直角坐标优化”这一方法没有多大作用，我试过，还是一样的结果 另外如果是根据振动频率的模式来调整分子结构的话，还要注意调整完后的分子结构 对应的能量是不是比未调整之前有所降低，如果没有，说明还是不行。

alwens兄所讲的第三点，经过我的猜测和与alwens的探讨，终于知道应该怎么做： 分子结构输入用直角坐标，将这些产生虚频的振动坐标直接与原始坐标相加。alwens兄进一步提到： 这个方法有时候不太好。

不一定将坐标全部加上。有时候可以选择加1/2等。主要的目的是去掉由分子对称所形成的虚频。

象上面的那个计算，完全可以再利用scf(tight)的选项消掉1－2个虚频。强度很小。

又多学了一招，自我感觉很幸福中„„^_^ 怎么把内坐标改成直角坐标啊 看你的输出文件

Full point group CS NOp 2 Largest Abelian subgroup CS NOp 2 Largest concise Abelian subgroup CS NOp 2 Standard orientation:-------Center Atomic Atomic Coordinates(Angstroms)Number Number Type X Y Z-------1 46 0-0.242506 1.028743 2.204311 2 46 0-0.476137 2.619863 0.000000 3 46 0-0.242506 1.028743-2.204311 4 46 0-0.242506-1.618028-1.382208 5 46 0-0.242506-1.618028 1.382208 6 46 0-1.725240 0.349678 0.000000

这一部分就是，如果没有，看 Input orientation:也行。

在opt的选项中加入“tight”或“saddle＝0”，若是计算过度态，应是“saddle＝1”。但是，也不一定能消除虚频！最好，将结构重新做一下，注意分子的构型，如CH3的三个H与相邻基团是交叉，还是重叠！请教虚频

请问如何知道虚频的结构是什么，是不是频率计算负值都可以认为是虚频，那我用gaussianview 看有三个负值，接下来我想做过渡态，用IRC，那我应该怎么弄呢？谢谢指导 你的输出文件中这部分就是看是否有虚频的

A“ A' A” Frequencies---65.9781 56.7053 59.7893 Red.masses--105.9032 105.9032 105.9032 Frc consts--0.2716 0.2006 0.2231 IR Inten--4.2436 0.0151 1.4107 Raman Activ--0.0000 0.0000 0.0000 Depolar--0.0000 0.0000 0.0000 Atom AN X Y Z X Y Z X Y Z 1 46-0.49 0.13-0.15-0.41-0.16 0.08-0.10 0.15 0.22 2 46 0.00 0.00 0.16 0.64-0.01 0.00 0.00 0.00-0.18 3 46 0.49-0.13-0.15-0.41-0.16-0.08 0.10-0.15 0.22 4 46-0.22 0.05 0.08 0.05-0.06 0.09-0.26-0.19 0.00 5 46 0.22-0.05 0.08 0.05-0.06-0.09 0.26 0.19 0.00 6 46 0.00 0.00 0.44-0.06 0.23 0.00 0.00 0.00-0.31 7 46 0.00 0.00-0.23-0.11 0.01 0.00 0.00 0.00 0.52 8 46 0.00 0.00-0.24 0.23 0.20 0.00 0.00 0.00-0.48 Frequencies 那一行数值为负，即有虚频。有三个负值，说明找的不是过渡态。但是如果另外两个虚频很小，而另外一个虚频很大，而且它的振动模式对应于你所计算的反应模式（比如反就前后某一键的拉长，缩短等），在计算所用的方法充许校正范围内，也还可以认为是找到过渡态。但是一般要求能够把另外两个虚频去掉。方法是要据振动模式，调整结构，再优化。

般认为，过渡态仅有一个虚频，中间体没有虚频！建议仔细阅读《Exploring Chemistry with Electronic Structure Methos》和《G98W程序使用入门》，后者稍微简单一点，厦大的免费资源里有！

有一个虚频是过渡态的必要条件；有几个虚频肯定不行了，有时候有一个虚频也不一定就找准了，做IRC成功就可靠了。另：复合物也会有小的虚频，只要不大于100，也可认为是一极小值。

没错,有些体系可能存在several or many local minima,那么找到的具有一个虚频的过渡态可能连接两个loccal minima,而非连接反应物和产物,因此必须做IRC来判断.关于虚频的物理意义

自恨理论功底浅薄，所以特地在此向诸位高手求教：

在寻找过渡态计算中opt=ts, 需要计算频率，如果有一个虚频就说明优化所得的构型是个 过渡态。有时在做部分限制优化opt=minimal时，也能得到一个虚频（一般发生在frozen 部分的化学键上）。我想问一下，虚频的数值大小本身具有什么物理意义。后一种情况是否也可以认定为过渡态。

我的理解：

虚频就是在这个振动方向上力常数是负的

过渡态之所以有一个虚频，是因为它是一个鞍点，鞍点上有一个负的梯度的方向 分子沿这个方向振动时将转化为反应物或产物。就象从山上掉下来，受到的力可以 认为是负的。

但并不是有一个虚频就是过渡态，还要在这个虚频的振动方向上分别指向反应物和 产物

有一个虚频并不一定是你想要的过渡态，你要看虚频的振动方向是不是对应着反应的方向。从概念上说，虚频就是化学反应方向上的简正坐标的二次偏导小于零。

另外，根据我的经验，这个虚频的负值应该大一些，如果这个值太小，那也不足信。不知对不对，呵呵！

优化输出信息：

这说明优化失败了，没有收敛，你可以尝试：

1。降低收敛标准（算一般的有机分子不要紧，如果是过渡金属，那就不要了）2。增大循环次数 3。换基组和算法 4。改变初始猜测

是有机分子阿，呵呵，用的是ＡＭ１，还是太大吗ｈｓｉ

这种不收敛的情况可能是你最初用软件建模的时候构型不合理，仔细考虑一下你的分子的构型是否合理，是否还存在其它的构型，单纯地降低收敛标准意义不大，不过你可以先尝试一下增大循环次数，最多scf(maxcyc=300),再不收敛的话就不是循环的问题，是构型不对。

S能不能说说，π轨道和σ轨道的组成上有什么差别？这方面我比较欠缺。ure.You can directly analyze the components of the molecule orbitals from the population in the *.log file.我在版面上说过了，当初的讨论还真是热烈：）你可以翻看一下

我这里简述几句：

π轨道主要是由npx，npy，nd－1，nd+1组成σ轨道主要是由ns，npz（键轴方向），nd0

还有，nd＋2，nd－2主要是组成δ轨道。

opt 是可以restart，前面的中间文件是少不了的。1, %chk=...#..opt Stopped 2, %chk=...#..opt(restart)guess=read geom=allcheck-----------------对弈freq的restart，是要采用数值计算方法才可以restart。（以下我省略chk和方法的输入，自己加上）1, freq(numerical)2, freq(restart,numerical)guess=read

OPT=restart不能进行，why 因为死机OPT被中断，我采用OPT=restart继续进行，但是出现以下情况：

Restoring state from the checkpoint file “Mn7O14V11O2-PCOPT1-2.chk”.FileIO operation on non-existent file.FileIO: IOper= 2 IFilNo(1)=-997 Len= 120004 IPos= 0 Q= 8710320

dumping /fiocom/, unit = 1 NFiles = 16 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 667136 FType=2 FMxFil=10000 Number 0 0 0 0 0 508 522 598 Base 144896 149504 142848 22016 444416 143360 143872 172032 End 146944 172032 143360 142336 667136 143375 144572 172034 End1 146944 172032 143360 142336 667136 143872 144896 172544 Wr Pntr 143360 149504 142848 22016 199168 143360 143872 172032 Rd Pntr 143360 149504 142948 142020 199168 143360 143872 172032 Number 634 991 992 994 995 996 998 999 Base 199168 147456 146944 20480 142336 21504 20992 172544 End 444168 149097 146949 20510 142346 21554 21042 198796 End1 444416 149504 147456 20992 142848 22016 21504 199168 Wr Pntr 199168 147456 146944 20480 142336 21504 20992 172544 Rd Pntr 199168 149097 146949 20510 142346 21554 21042 172545

dumping /fiocom/, unit = 2 NFiles = 5 SizExt = 0 WInBlk = 512 defal = F LstWrd = 267264 FType=2 FMxFil=10000 Number 0 508 522 634 998 Base 266245 20480 20545 21245 20495 End 267264 20495 21245 266245 20545 End1 267264 20495 21245 266245 20545 Wr Pntr 266245 20480 20545 21245 20495 Rd Pntr 266245 20480 20545 21245 20495

dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 66048 FType=2 FMxFil=10000 Number 0 Base 20480 End 66048 End1 66048 Wr Pntr 20480 Rd Pntr 20480 Error termination in NtrErr: NtrErr Called from FileIO.为什么我的file IO会不存在？是什么参数设漏了吗？请各位大虾支持，谢谢！

在你前一次的计算中，第一单点的计算未完成，所以，checkpoint文档中尚无优化过程的信息，无法restart.过渡态优化过程中遇到一个这样的问题

我在过渡态优化过程中用这样的方法 # hf/6-31g(d)opt=(ts,calcfc)freq当运行到L999时Link died请问可能是什么原因？是不是我的过渡态输入有问题？另外，寻找过渡态有哪些技巧，请指教！

你的计算在设定的步数里面没有找到过渡态。这有两个可能：一是初始输入的构型不好。二是循环次数不够(可以maxcycle=大一点的数)。具体还是等高手来指点吧。

对于决定收敛结构的标准，做了一些很小但很重要的变动。当力比截断值小两个数量级时（即，极限值的1/100），结构即被认为是收敛的，即便是位移比截断值还大。这个变动有利于大分子的优化，因为它们的最小值周围可能有非常平缓的势能曲面。(G98manual p106)所以你的opt有2个yes 频率计算完成后还有一次优化，可能这次就优化到3个了

！做频率分析的时候只需要用Freq就可以，不需要加Raman选项，对吗？ freq=noraman不作拉曼强度计算,可以提高10%-20%的计算速度.算一个分子的电离能，里面有I，CL 在优化起正离子的时候老是没有没有结果 体系的能量趋势不固定，有时大有时小，MAXCYCLE=80后出现2个YES，然后DIE LINK 999 但是COPY。OUT最后的构型再去优化后 结构4个全是NO 现在还在算，但是还是NO 怎么办？ 要OPT=CALCALL吗？有什么用？ 另外，还有一个题外的问题 这样写对吗？

OPT（MAXCYCLE=80，STEPSIZE=10）可以吗？

增大优化次数，optcyc=N;N给大一些。opt=(maxcycle=N)也一样

前面的老帖子已经讲到这些，好象在免费资源里的Gaussian98 tips.一部分引用如下：

First, whenever you encounter SCF convergence failures you need to look at the progress of the SCF by using #P in place of # which will print out the status at each iteration.It is useful to know if the SCF is converging slowly, oscillating, taking large jumps etc.Slow convergence can be a sign that a higher order converger like SCF=QC would be helpful.Oscillations can be a sign that a small HOMO/LUMO gap is causing a flip in the occupied orbitals and SCF=VShift=n is useful.Large jumps can indicate that a poor initial guess or a change in the electronic state from the initial guess has occured and normal convergence will eventually achieved but after a larger number of cycles and a better initial guess would help.Second, whenever possible it is helpful to start with a smaller basis set, i.e.STO-3G, LANL1MB, etc for single points.This is cheaper and so you can often generate a better initial guess which can be used with GUESS=READ and the larger basis set and without the more expensive options.In your case the SCF converges initially quite smoothly and then bounces around without making progress.This makes SCF=QC a likely approach to clean up this problem, at least at the initial point.You might want to look at adding IOP(5/13=1)which allows the SCF to print out the MO's even in the event that convergence was not reached.This should only be used with single points, i.e.remove OPT.As to optimization constraints.To freeze atomic coordinates you should turn on NoSymm.This method of freezing structures causes problems for redundant internals which are being looked into but NoSymm fixes the most serious ones.下面是血红素的化学合成的模型，由于分子有点大，又加上有个Fe，所以优化非常困难，请求斑竹支 高招：

（1）这是一个高度对称的分子，如何利用加虚原子，提高计算速度？

（2）优化Fe，我采用了LANL2DZ，请问有没有更好的基组？怎么手工加基组？（3）由于变量太多，您认为我通过怎样的办法，减少变量来提高计算速度？ 非常感谢！！做出文章了，我一定挂你的大名，呵呵。

S O 1 B1 Fe 2 B2 1 A1 N 3 B3 2 A2 1 D1 N 3 B4 2 A3 1 D2 N 3 B5 2 A4 1 D3 N 3 B6 2 A5 1 D4 C 4 B7 3 A6 2 D5 C 8 B8 4 A7 3 D6 C 9 B9 8 A8 4 D7 C 4 B10 3 A9 2 D8 C 5 B11 3 A10 2 D9 C 12 B12 5 A11 3 D10 C 13 B13 12 A12 5 D11 C 5 B14 3 A13 2 D12 C 6 B15 3 A14 2 D13 C 16 B16 6 A15 3 D14 C 17 B17 16 A16 6 D15 C 6 B18 3 A17 2 D16 C 7 B19 3 A18 2 D17 C 20 B20 7 A19 3 D18 C 21 B21 20 A20 7 D19 C 7 B22 3 A21 2 D20 C 23 B23 7 A22 3 D21 H 24 B24 23 A23 7 D22 C 11 B25 4 A24 3 D23 H 26 B26 11 A25 4 D24 C 16 B27 6 A26 3 D25 H 28 B28 16 A27 6 D26 C 20 B29 7 A28 3 D27 H 30 B30 20 A29 7 D28 H 1 B31 3 A30 2 D29 H 9 B32 8 A31 4 D30 H 10 B33 9 A32 8 D31 H 13 B34 12 A33 5 D32 H 14 B35 13 A34 12 D33 H 17 B36 16 A35 6 D34 H 18 B37 17 A36 16 D35 H 21 B38 20 A37 7 D36 H 22 B39 21 A38 20 D37 B1 4.21561119 B2 1.64053835 B3 2.01170682 B4 2.00223534 B5 2.01140630 B6 2.01948030 B7 1.38502525 B8 1.45027869 B9 1.36206816 B10 1.38792662 B11 1.38546948 B12 1.45027336 B13 1.36191526 B14 1.38377173 B15 1.38899517 B16 1.45046091 B17 1.36116892 B18 1.38527197 B19 1.39574494 B20 1.44813170 B21 1.36193008 B22 1.39696571 B23 1.37779525 B24 1.08374598 B25 1.38063627 B26 1.08377153 B27 1.37911143 B28 1.08372124 B29 1.37753366 B30 1.08361344 B31 1.37864549 B32 1.07802759 B33 1.07801256 B34 1.07797961 B35 1.07795952 B36 1.07798544 B37 1.07794319 B38 1.07803016 B39 1.07806743 A1 6.31961440 A2 94.07195309 A3 95.30512753 A4 94.38669042 A5 94.58713162 A6 127.10303272 A7 109.89556522 A8 107.10994453 A9 126.54307980 A10 126.81263629 A11 109.85724161 A12 107.06302128 A13 126.91593763 A14 126.68048694 A15 109.77898998 A16 107.13028054 A17 126.98352930 A18 126.67933314 A19 109.65642922 A20 107.43086078 A21 126.85650291 A22 125.36309437 A23 117.60829677 A24 125.56378588 A25 117.60334742 A26 125.52442366 A27 117.63863540 A28 125.46104554 A29 117.64536447 A30 98.48331793 A31 124.49018481 A32 128.37464062 A33 124.51087771 A34 128.41640662 A35 124.47708155 A36 128.39197501 A37 124.32863598 A38 128.24305473 D1-108.51711460 D2 161.17665430 D3 70.83221148 D4-18.82295908 D5 85.79136868 D6-174.21655652 D7 0.09638172 D8-87.07481057 D9 86.40576040 D10-175.60372334 D11-0.05103700 D12-88.28171080 D13 88.94769251 D14-174.64817315 D15 0.35121394 D16-84.34232126 D17 84.15970765 D18-172.91552513 D19 0.72935930 D20-85.63442533 D21-5.66240036 D22 179.49688145 D23-6.19734957 D24-179.54995066 D25 4.93855843 D26 179.77279426 D27 6.47288831 D28-179.56429342 D29-89.13791474 D30-179.79660459 D31-179.88612542 D32-179.83142492 D33-179.77466949 D34-179.60732447 D35-179.84766166 D36-179.39904616 D37 179.72163849

我也在做大分子的计算，希望能交流交流。我的电子邮件：chsun@imr.ac.cn 我的经验是先做PM3优化，然后再用精度高的方法和基组。

你的前三个原子，不是一条直线吗？如果是直线，那对称性就可以是C4v了。如果不是直线，那至少要做成cs对称性吧？把前三个分子所在的平面作为对称面。你现在的对称性是没有——虽然看上去很对称。Fe用赝势基组不知道会怎么样，你可以看看别人的文章。

做优化的时候，可以固定一些外层的原子，只让和中间的原子相连的动。

大约有3M 另外，Fe、O、S分子是在一条直线上么？如果不是，你可以将吡啶环上的参数先减少了，然后再加Fe、O、S，肯定能够减少不少参数，但是不知道由没有实际意义了

1）我谈一下我做过的体系加虚原子的经验，希望对你有帮助：一般情况下，虚原子加在分子的对称中心，或键轴Cn轴上，如要加第二个，可放在垂直于Cn轴的方面上。（2）你可以到基组库去寻找基组，可能要比单单的Lanl2Dz效果好

（3）你的体系高度对称，肯定有不少的CH键长键角是一样的（或者差别很小），一律设为一个参数，让它同步优化。

PS：你的输入文件的参数也太多了点

我看了一下你的结构，你说是高度对称？我想问你一下这是你手建的吗？从你的坐标看好象没有对称性。32号H怎么位置这么怪！我认为你的结果最高对称性是 C(4V)，最好手建一下为好（如不考虑外场作用）。

这种体系的计算现在非常热门也非常有难度，上次我们听的以色列教授介绍的也是这方面的内容，这种体系以我们的水平可能很难处理好。你可能要找到做这方面的专家指导才会有很好的效果，否则很难分析清楚计算出的结果。

上次以色列教授沙松过来我们这边就是介绍的有关P450的计算，虽然我没有听懂什么，不过以我自己计算含金属的小分子经验来看，关于这种体系很难处理的部分就在于如果控制中心金属原子铁和周围配体的配位情况，以及了解中心金属原子的电子占据，如果你仅是靠软件产生体系的输入坐标，可能计算很难达到你想要的结果。这些坐标的产生可以要靠自己的经验来步步设置，所以就像你所说的太细节的东西是很难说清楚的，因为他们在计算过程中拥有比你更多的经验。只能给你这些参考性建议了，没有什么实质性的帮助。

我的输入文件是：

%chk=zux3

#p hf/lanl2dz opt symm=loose scf=(maxcycle=200)test

zux3

0 1

0.000000 0.000000 1.463625 6-2.059783 0.151849 1.234523 8-3.189332 0.235187 0.995024 6 0.153438 2.073727 1.558845 8 0.238770 3.225066 1.644723 6 0.000000 0.000000 3.544137 8 0.000000 0.000000 4.703699 6 0.000000 0.000000-0.664692 6 0.000000 0.000000-1.891052 6 0.000000 0.000000-3.251669 7-1.186271-0.125670-3.989741 6-2.427060-0.359632-3.218024 6-1.381218 0.710276-5.201852 1-2.263789-1.158724-2.497766 1-2.749787 0.537294-2.683916 1-3.209142-0.660105-3.913517 1-0.432645 1.124763-5.529600 1-1.810550 0.110340-6.005013 1-2.062473 1.534158-4.980011 6 2.059783-0.151849 1.234523 8 3.189332-0.235187 0.995024 6-0.153438-2.073727 1.558845 8-0.238770-3.225066 1.644723 7 1.186271 0.125670-3.989741 6 2.427060 0.359632-3.218024 1 2.263789 1.158724-2.497766 1 2.749787-0.537294-2.683916 1 3.209142 0.660105-3.913517 6 1.381218-0.710276-5.201852 1 0.432645-1.124763-5.529600 1 1.810550-0.110340-6.005013 1 2.062473-1.534158-4.980011

分子的点群是c2，但是优化第二轮是就出错了，Omega: Change in point group or standard orientation.Error termination via Lnk1e in d:g98wl202.exe.如果我不用关键词symm=loose,优化时分子点群就变成了c1。

请高手指教，这是何原因呢？若我想保持分子仍是c2，改如何处理？

3x!换内坐标吧

opt=Z-matrix 你用内坐标可能会好一点吧？而且也能保证你的对称性的。[em12][em12] Angle between quadratic step and forces= 29.72 degrees.Linear search not attempted--first point.Iteration 1 RMS(Cart)= 0.00036483 RMS(Int)= 0.00000015 Iteration 2 RMS(Cart)= 0.00000015 RMS(Int)= 0.00000006 Variable Old X-DE/DX Delta X Delta X Delta X New X(Linear)(Quad)(Total)R1 2.05001-0.00006 0.00000-0.00017-0.00017 2.04984 R2 2.83834 0.00000 0.00000 0.00009 0.00009 2.83844 R3 2.05240-0.00007 0.00000-0.00023-0.00023 2.05218 R4 2.05240-0.00007 0.00000-0.00023-0.00023 2.05218 R5 2.03466-0.00003 0.00000-0.00007-0.00007 2.03459 R6 2.47511 0.00010 0.00000 0.00015 0.00015 2.47526 R7 2.02743-0.00009 0.00000-0.00024-0.00024 2.02719 R8 2.51854 0.00015 0.00000 0.00038 0.00038 2.51892 A1 1.94852-0.00009 0.00000-0.00069-0.00069 1.94783 A2 1.88382 0.00007 0.00000 0.00042 0.00042 1.88424 A3 1.88382 0.00007 0.00000 0.00042 0.00042 1.88424 A4 1.93636-0.00006 0.00000-0.00034-0.00034 1.93602 A5 1.93636-0.00006 0.00000-0.00034-0.00034 1.93602 A6 1.87204 0.00008 0.00000 0.00063 0.00063 1.87267 A7 2.06277 0.00007 0.00000 0.00058 0.00058 2.06335 A8 2.15661 0.00001 0.00000 0.00012 0.00012 2.15674 A9 2.06380-0.00008 0.00000-0.00071-0.00071 2.06310 A10 2.19767 0.00004 0.00000 0.00041 0.00041 2.19808 A11 2.13077-0.00008 0.00000-0.00050-0.00050 2.13027 A12 1.95474 0.00004 0.00000 0.00009 0.00009 1.95483 D1 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D2 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 D3 1.03955 0.00001 0.00000 0.00017 0.00017 1.03972 D4-2.10204 0.00001 0.00000 0.00017 0.00017-2.10187 D5-1.03955-0.00001 0.00000-0.00017-0.00017-1.03972 D6 2.10204-0.00001 0.00000-0.00017-0.00017 2.10187 D7 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 D8 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D9 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D10 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 Item value Threshold Converged? Maximum Force 0.000148 0.000450 YES RMS Force 0.000059 0.000300 YES Maximum Displacement 0.000912 0.001800 YES RMS Displacement 0.000365 0.001200 YES Optimization completed.--Stationary point found.详细解释一下好吗，这一段几乎全看不懂

Iteration 1 RMS(Cart)= 0.00036483 RMS(Int)= 0.00000015 Iteration 2 RMS(Cart)= 0.00000015 RMS(Int)= 0.00000006 这应该是按笛卡尔坐标和内坐标计算的均方误差，比较的标准不清楚！

Variable列是输入的变量；Old X列是上一次的变量值；-DE/DX是能量对坐标的一阶导数，是判断驻点的依据；Delta X是与上次相比变量的变化；New X是新的变量数值！

SCF Done: E(RHF)=-113.866331170 A.U.after 12 cycles Convg = 0.1529D-08-V/T = 2.0020 S**2 = 0.0000 NROrb= 34 NOA= 8 NOB= 8 NVA= 26 NVB= 26 最后一行是什么东东？

总轨道数34，Alpha占居轨道8个，Beta占居轨道8个，Alpha虚轨道26个，Beta虚轨道26个！NROrb＝number of orb NOA= numberof alpha occupied orb NVA=number of virtual alpha orb

结构优化的时候有没有可能出现这种情况：

即构型在两个或几个结构循环变化，总也不能达到终点 可能出现么？如果这样的话该如何处理？非线形搜索？ 多谢！

“构型在两个或几个结构循环变化，总也不能达到终点，可能出现么？”这个问题是你实际遇到的，还是凭空想出来的？几个结构在势能面上对应着局部最小，自洽场的结果是会到达一个局部最小值，而不会出现在几个个结构循环变化。我觉得你想说的是收敛振荡的问题吧？

他指的可能是结构优化振荡，我曾经用过iop（1/8)里的一些选项，对控制结构优化有一定帮助

我现在遇到的问题是这样的，优化了很多圈，但每次的能量总是在很小的范围内起伏

只能在Force上达到两个YES，而Displacement上总是NO,我觉得就象是在几个结构中振荡一样，这是所谓的优化振荡么？

iop（1/8)里有设定优化步长的功能：

0 DXMAXT = 0.1 BOHR OR RADIAN(L103, Estm or UnitFC).= 0.3 Bohr or Radian(L103, Read or CalcFC).= 0.2 Bohr or Radian(L105).= 0.3 Bohr or Radian(L113, L114).N DXMAXT = 0.01 * N 哪里有关于iop（1/8)的详细内容呢？？

在免费资源版有IOP的pdf文件

谢supi!使用了iop(1/8=1)结构终于优化完了，但这说明什么呢？体系的能量对于结构的微小变化很敏感？还是说势能面比较平坦？

版上的今天刚贴的Guassian tips这篇贴子中最后有提到收敛遇到的一种情况就是“taking large jumps”.对于优化默认采用的步长为0.3Bohr对你优化的结构来说太大,所以你看到结果能量总是在很小的范围内起伏.因此你使用IOP(1/8=1)设置步长为0.01Bohr,而达到收敛的最终结果,说明这个结构的势能面是起伏不平的,对结构的变化非常敏感.不过,这种情况也像是oscillating,建议用SCF(Vshift=300)试一下,看能不能收敛到同一结果?

试了好几次,总觉得使用二次收敛方法scf=qc,对收敛没有什么太大的效果,反而不如以前.而且对于oscillating采用SCF(Vshift=N)来说,由于DFT已经使用N=200,即使增大N的值对收敛的结果也无明显效果,比不上减小优化步长.但是对于减少步长来说,机时消耗太大,虽然在能量上比前者来得低,但是仍可能达不到收敛的预期效果,因此这些都说明了一个问题:The initial guess isn't very good.关于如果处理这些问题,《计算化学》上给出了一些更好的建议。

看了关于学术腐败的讨论，很惭愧，学习高斯已经有好几个月，还停留在Opt阶段，而且还有不少问题，先问两个吧，请莫见笑。

1.我发现构型优化时，Opt(tight)和只有Opt时，机时差别若干倍，但结果好像差的不是很多，为什么？ 2.对于大分子，可以先用小基组优化，而后再用更高级基组优化，具体写输入文件时，小基组时的分子说明部分用初始分子构型的数据，而更高基组优化时的分子说明部分应该用小基组的计算结果，怎么写呢？

。opt=tight采用了更为严格的收敛标准，使用选项Opt=Tight的优化计算比使用默认截止点的优化计算要多算几步。所用的时间当然也大得多了。要注意我们使用的能量的原子单位是个很大的数值1au＝627.51kcal/mol，也许用两种方法算出来的能量只差了0。05个单位，可换算一下就有30多的kcal/mol！2。你可以先用小基组算，然后用前面的chk文件来作大基组的计算。在控制行加入geom＝check，guess＝read，就不用分子构型输入了

结构优化不收敛，一般有什么OPT=ReadFc, opt=CalcFc, opt=calcAll,或者设置maxcycle,opt=restart等等，但有没有人知道那种最为有效？具体方法怎样？当然有示例是最好的了！呵呵：）谈谈我的想法: 我感觉结构不收敛可能与你的初始猜测有关,也许先固定某些键长键角再优化,事后再放开优化会比较容易收敛.采用低级算法先计算只是其中的方法之一,你为何一直要用这种方法?有什么依据? 如果你顺利解决了你的计算,那就谈谈你的经验如何?

我觉得采用OPT=CALCFC并不一定得到好的结果,有时候收不收敛实际和我们处理的化学问题相关,而不仅仅是用计算上的技巧进行解决.比如两个原子间电子占据没有sigma成键,但是存在强的sigma非键作用,那么它们之间的斥力很强,计算就很难收敛形成稳定的体系.因此在处理收敛的问题时,我们应该多分析不收敛的问题实质,是电子的什么运动状态引起的,而非看到不收敛就想利用一些key来解决,我自己也试过一些key,包括opt=calcfc,allcfc,NoDIIS,QC等等(包括一些强制性的手段),并不是很有效,因为它们都是从技术的角度来解决问题.如果一个体系能形成稳定的体系,那么它不依靠这些key也是能解决的.反而如果一体系不收敛只是利用了这些技术就收敛了,那么这种结果如何可信?前面的贴子中,有些技巧诸如调整原子间的距离和键角,这些手段才是比较有效的,因为它直接影响了原子间的电子运动状态,即波函数,我们有时候初始猜测不正确就是因为程序给出的电子的占据状态错误而直接影响了原子间的成键,这时就需要对电子占据的轨道进行调整,Guess=alter,这对于过渡金属的计算来说特别要注意,因为G98总不能给出正确的初始猜测.谢谢诸位意见，收益非浅。我起初一直认为我的初始模型是正确的，但后来我试过4个模型，才发现初始结构非常重要。开始三个是MnO2理想的晶体结构，通过改变键长希望结构优化收敛，但是失败。我又重新定义采用MnO2的真实晶胞参数，现在结果有些眉目，能收敛只是没有正常结束，可能是我吸附分子的坐标出了点问题，出现

Error on Z-matrix card number 23 angle Alpha is outside the valid range of 0 to 180.所以我深切体会到楼上师兄的话：调整原子间的距离和键角,这些手段才是比较有效的。我现在才发觉，G98真是一个费内存，费时间的东西，很多计算要的是经验，你计算的垃圾多了，你才知道你需要什么！这个过程真够悲惨。

第五篇：节能改造服务合同

附件

编号：

上海市节能改造服务合示范文本（2014版）

上海市工商行政管理局

制定

同上海市节能工程技术协会

使用说明

一、本合同文本是根据《中华人民共和国合同法》、《中华人民共和国节约能源法》等有关法律法规制定的示范文本，供双方当事人约定采用，合同签订前请仔细阅读。

二、节能服务行业专业性强，涉及诊断、设计、改造、系统集成等内容，为更好地维护双方当事人的权益，双方签订合同时应尽到审慎义务，力求合同条款约定具体、严密。

三、节能服务公司应具有相应的资质，用能单位签订合同前应验看节能服务公司的《企业法人营业执照》和企业资质证书。

四、本合同中名词的含义：

用能单位：指在合同中约定的具有项目发包主体资格，承担支付项目价款义务的当事人。

节能服务公司：指在合同中约定的被用能单位接受的具有项目承包主体资格的当事人。

项目：指用能单位与节能服务公司在合同中约定的承包范围内的节能改造服务内容。

项目价款：指用能单位与节能服务公司在合同中约定，用能单位用以支付节能服务公司按照合同约定完成承包范围内全部项目的款项。

节能率：报告期节能量与基准期的能源消耗量之比，即采取节能措施之后节约的能源量与未采取节能措施之前能源消费量的比值。

五、经备案的节能服务公司在实施合同能源管理模式的节能改造项目时，应按GB/T24915的要求执行。

六、本合同文本自发布之日起使用。今后凡未制定新的版本前，本版本延续使用。

合同编号：

上海市节能改造服务合同

（2014版）

甲方（用 能 单 位）： 乙方（节能服务公司）：

根据《中华人民共和国合同法》、《中华人民共和国节约能源法》等有关法律法规的规定，结合本节能改造服务项目（以下简称“项目”）的具体情况，甲乙双方在自愿、平等、诚实守信的基础上，经双方协商一致，签订本合同。

第一条 项目内容

1.1项目名称：。1.2项目地址：。1.3节能改造服务范围：。1.4节能要求：□□节能率 □节能量 □节能效益。第二条 项目工期

2.1施工日期：。2.2竣工日期：。第三条 项目标准 3.1能耗基准：。3.2测量与验证方法

□ 方案A：隔离改造部分：测量关键参数 □ 方案B：隔离改造部分：测量所有参数 □ 方案C：整个耗能设施 □ 方案D：经校准的模拟

□。3.3执行技术标准：。第四条 项目总价

本合同总价款为人民币 元，（大写）元（价款明细另列附件）。

第五条 付款方式。

第六条 项目设计

6.1本合同签订后 日内，双方根据能耗统计结果及节能潜力协商确认节能改造设计方案。

6.2双方根据方案确定待采购的主要设备和原材料的清单作为本合同附件，乙方于双方确认后开始项目的深化设计和设备的采购。

6.3在节能改造项目实施过程中，甲方提出变更设计的，应当签订项目设计变更单。

第七条 项目实施

7.1乙方的节能设备和原材料到达现场时间 年 月 日，双方应对设备到达现场的事实予以确认。

7.2节能改造过程中，乙方发现甲方指定或提供的原材料、设备有质量问题、工程缺陷或规格差异的，应及时向甲方提出。甲方认为可以继续使用的，应向乙方作出书面确认。

7.3乙方应按照双方确认的设计文件安装到位，设备稳定运行后 个工作日内，双方应进行项目竣工验收。

7.4对竣工验收或测量与验证结果存在争议的，双方可共同委托具有相应法定资质的第三方检测机构对争议事项进行检测确认。

7.5双方约定维保期限自竣工验收之日起 月，维保日期为（每月/季/半年）。

7.6甲乙双方竣工验收后，乙方应自确认签字之日起承担质保责任，质保期自竣工验收合格之日起 月，质保范围为。对实施节能量分享型的节能改造项目，乙方应实施分享期内的质保服务。

第八条 甲方权利义务

8.1甲方应提供节能改造项目的设计、施工所需的交通、环卫和防治施工噪音管理等必备手续。

8.2甲方应根据合同约定提供项目的相关资料，并确保其真实、准确、完整。

8.3乙方运送至甲方的设备、原材料、施工工具等物品，甲方在安装前负有无偿保管义务，并应提供存放场地或仓库。

8.4在项目施工阶段，甲方应为乙方提供项目实施的工 作条件：

（1）进出甲方项目现场的通行证；

（2）项目场地在甲方专属区域的，甲方应提供办公和生活场所等，并负责现场保安工作；

（3）提供节能项目实施所需要的现场条件，如提供动力源、管道图纸、清理施工现场等；

（4）为利于节能改造项目实施而合理地调整生产经营；（5）其他必要的条件。

8.5对乙方提交的设计、施工方案在收到之日起__个工作日内以书面形式予以确认；如甲方认为乙方提交的设计、施工方案与合同不符，应在该期限内提出书面意见。

8.6甲方应根据项目操作规程和保养要求，对设备进行操作、维保，以保证系统与设备正常运行和满足设计要求，并对乙方服务质量进行评价。

8.7如设备发生故障、损坏，甲方应在获悉情况后 个工作日内书面通知乙方，并对乙方检测和维保工作予以配合。

第九条 乙方权利义务

9.1乙方应根据本合同约定的节能设计方案、采购和供应相关设备，进行施工、安装和调试。

9.2项目开工前，应将设计、施工等资料及项目计划表提交甲方。

9.3在收到甲方确认开工之日起 个工作日内，书面答复甲方对设计、施工方案提出的要求或意见。9.4项目竣工合格验收后，乙方应对甲方指派的操作人员进行系统的工作原理、操作规程、常见故障及处理措施等的免费培训。

9.5乙方应做好分项验收记录、竣工验收记录、项目结算记录、项目竣工图，设备和系统调试记录及报告，设备和原材料合格证等的归档工作。项目结算后应移交本项目档案资料及继续运行所必需的相关技术资料。

第十条 违约责任

10.1由于乙方原因逾期竣工的，每逾期一日，乙方支付甲方本项目价款 ‰的违约金。

10.2由于甲方原因导致延期开工或中途停工的，甲方应补偿乙方因停工所造成的损失。需乙方后续进场施工的，甲方承担后续进场费用。

10.3如甲方未按合同约定付款的，每延误一日，按未支付部分价款的 ‰向乙方支付违约金，逾期付款金额达到合同总价款的 %，乙方有权解除本合同。

10.4如乙方提供的节能改造设备存在质量问题或节能改造方案存在瑕疵，致使节能设备不能稳定运行，则甲方可要求乙方修理更换，并赔偿相应损失。

10.5一方违反保密义务，致使另一方损失的，应当予以赔偿。

第十一条 知识产权

本合同涉及的专利实施许可和技术秘密许可，双方约定如下：。

第十二条 其他约定。

第十三条 争议解决

双方发生争议的，可协商解决，或向有关部门申请调解：也可提请上海仲裁委员会仲裁（不愿意仲裁而选择向法院提．．．．．．．．．．．．起诉讼的，请双方在签署合同时将此仲裁条款划去）。．．．．．．．．．．．．．．．．．．．．．第十四条 附则

本合同自双方签字或盖章之日起生效。本合同一式两份，具有同等效力，甲乙双方各执一份。

甲方（盖章）： 乙方（盖章）： 住

所： 住

所： 法定代表人： 法定代表人： 授权代表签字： 授权代表签字： 电

话： 电

话： 传

真： 传

真： 开户银行： 开户银行： 帐

号： 帐

号：

签约日期： 签约地点：