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

user: 2. 相关概念 2.1. 网络安全态势感知 网络安全态势感知的最早定义是1995年由Endsley [ 2 ] 提出，认为“态势感知是指在一定时空条件下，对环境因素的获取、理解和对未来的预测”，也就是说，态势感知首先需要通过一定方法获取过去与当前的各类环境要素，然后通过分析与研究从环境要素中有所“感知”，从而对未来态势的发展有所预测。2010年，Springer [ 3 ] 等人再次给出更加细化的定义，认为“网络空间(防御)态势感知涵盖态势识别、态势理解和态势预测三个阶段，至少包括7个方面的内容：态势认知、攻击影响评估、态势跟踪、对手趋势和意图分析、态势因果关系与取证分析、态势信息质量评估、态势预测”，这一定义中包含了态势感知具体工作的流程与技术七个方面，这七方面概括起来又可归纳为三个层面，即“态势识别”、“态势理解”和“态势预测” [ 4 ] 。 2.2. 威胁情报 威胁情报的定义随着时代和技术的发展逐步演进 [ 5 ] 。2013年美国IT咨询公司Gartner对威胁情报作了较为权威的定义：“威胁情报是一种基于证据的知识，它就网络资产可能存在或出现的风险、威胁，给出了相关联的场景、机制、指标、内涵及可行的建议等，可为主体响应相关威胁或风险提供决策信息 [ 6 ] 。”此外，还有国内学者 [ 7 ] 提出“威胁情报是通过大数据、分布式系统或其他特定收集方式获取，包括漏洞、威胁、特征、行为等一系列证据的知识集合和可操作性建议”。 本文选取Gartner对威胁情报的定义，并且将实际应用中的威胁情报分为三个层面：一是数据层面，例如基础的域名、IP、证书、CVE信息等，这一层面的情报一般数量巨大，但仅仅是最基本的数据采集的结果，研究和分析人员并不能从中直接得到具有指导性的信息；第二层是信息层面，是在数据层面情报的基础上，通过数据的统计、关联分析与可视化，得到事件或攻击组织层面的具有一定指导意义的情报信息，例如黑白名单、关于恶意域名/IP的地域统计分布等；第三层面是知识层面，是在信息层面情报的基础上，再进行深层的多源异构数据关联分析，得到对某一攻击事件或某一黑客组织的较为完善的分析报告，并且在一定程度上含有对未来安全态势的分析和预测，提出较具有参考价值的安全防护建议等。
system: [SYSTEM]

user: 4. 总结 在大数据、云计算、物联网等新型信息技术飞速发展的背景下，网络空间威胁也向朝泛化和复杂化的趋势发展，各类网络攻击也更加具有持续性和隐蔽性。基于威胁情报进行网络安全防御能够及时分析已发生的入侵，对未来威胁态势进行预判，并据此评估潜在的安全风险以指导用户制定有效的安全决策，系统化增强网络空间防御能力。本文通过分析威胁情报在态势感知中的基本应用环节，探讨了威胁情报对于需要安全防护的企业、风控与应急响应部门、日常运维人员和网络安全专家等不同角色人员的应用；基于特定需求，介绍了威胁情报在典型场景中的应用流程。随着威胁情报技术的不断成熟及相关规范、框架的日渐完善，在安全问题日益严峻的未来，它将进一步在社会、国家乃至全球的各个层面发挥巨大的作用。
system: [SYSTEM]

user: 摘 要 通过云南省昆明、曲靖和玉溪3个典型城市2014年1月1日~2017年12月31日空气污染状况数据，分析了PM 2.5 和O 3 的时间变化特征，讨论了冬季两个典型污染观测期间的气团后向轨迹。结果表明：2014年至2017年PM 2.5 质量浓度呈下降趋势，浓度水平受污染源排放和气象条件影响，存在明显的季节和时段变化特征；O 3 浓度呈逐年上升的趋势；冬季两个典型污染观测期间采样点的气团后向轨迹显示，云南重污染天气的出现，与外来源叠加有密切关系，最可能的是云南西南部污染和来自缅甸等国家的长距离传输叠加导致。 关键词 :PM 2.5 ，O 3 ，空气污染，后向轨迹分析 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 自2013年9月10日《大气污染防治行动计划》实施以来，我国大气质量有明显改善，统计2015年1月1日~2017年12月31日中国空气质量在线监测分析平台公布的367个城市中数据有效完整的320个城市的大气污染监测情况，PM 2.5 年平均值分别为50.5 μg·m −3 、47.0 μg·m −3 和44.9 μg·m −3 ，呈逐年下降趋势，然而现状仍不容乐观，2017年PM 2.5 年平均值是目前我国《环境空气质量标准》(GB3095-2012)中规定的二类区35 μg·m −3 限制的1.28倍，320个城市中有235个城市超出我国PM 2.5 二类区标准，超标率高达73%，若参照世界卫生组织给出的PM 2.5 标准10 μg·m −3 ，则320个城市中，达标率为0。 空气污染没有安全阈值，浓度越高，危害越大，每增加10 μg·m −3 的PM 2.5 浓度，记忆力和注意力测试中的表现水准就会有相当于衰老两岁的下降危害，而且总死亡率，心肺疾病死亡和肺癌的死亡率会增加 [ 1 ] 。地表空气中的O 3 对人体极为有害，一些易于过敏的人长时间暴露在臭氧含量超过180 μg·m −3 的环境中，会产生皮肤刺痒、眼睛刺痛、呼吸不畅、咳嗽及鼻炎等症状。空气中每增加100 μg·m −3 的O 3 ，人的呼吸功能就会减弱3%。因此，研究并尽快改善我国大气质量是目前急需解决的环境问题。 云南省地理位置特殊，同时与缅甸、越南和老挝三个国家接壤，云南省境内的空气质量既受到本地排放的影响，也受到省外和国外的传输影响。1997年，常蓉等 [ 2 ] 最先对云南省大气污染控制进行了探讨，研究表明云南省大气污染主要来源于工业生产、居民生活和汽车废气，以煤烟型为主。1997年，沈鹰等 [ 3 ] 首次对云南省大气污染扩散输送特征进行了研究，结果表明云南省大气稳定度以中性类居多，频率占50%以上，风速廓线指数值比平原地区大得多，也比一般山区的大，大气混合层高度春冬季较大，夏秋季较小。另外，罗燕等 [ 4 ] 利用2015年1月~2016年12月云南省16个城市大气污染物浓度监测结果与同期气象观测资料，分析了云南大气污染物分布的时空特征及其与气温、相对湿度、云量、降水、风速、日照时数等气象要素的关系。结果表明：云南SO 2 、NO 2 、PM 2.5 大气污染物浓度冬季最高，春秋季次之，夏季最低；PM 10 、O 3 平均浓度春季最高，冬季次之，夏季最低。CO平均浓度春秋季最高、冬季次之，夏季最低。云南CO浓度与气象要素相关关系不显著，其明显的月变化主要与污染排放源的变化有关。SO 2 、PM 10 、PM 2.5 及NO 2 浓度在不同地区对不同气象要素敏感性不一致，地区差异性较大。尚小宴等 [ 5 ] 等利用普洱市环境监测站收集的普洱市思茅城区2013~2014年度API指数，对API指数进行时间序列分析和空间分析，包括2013~2014年API的质量状况、季节变化以及比较南北市区的API值。结果表明：普洱市城区的空气状况在2013~2014年度优良，夏季的空气质量最好，春季最差；空间上存在南北差异，南市区API指数低于北市区，南市区空气质量优于北市区。Qu等 [ 6 ] 第一次分析了2000~2006年我国86个城市的PM 10 浓度特征，曲靖的PM 10 超过了80 μg·m −3 。毕丽玫等 [ 7 ] 通过采集PM 2.5 四季环境样品，分析了昆明PM 2.5 质量浓度及其4类化学组分特征。许燕婷等 [ 8 ] 基于2014~2016年中国环境保护部公布的城市空气质量指数(AQI)日报数据，综合运用空间插值法、空间自相关模型、核密度估计模型等空间数据统计模型探究中国城市空气质量时空分布及动态演进特征，结果显示：2014~2016年中国城市AQI季节均值为：冬季 > 春季 > 秋季 > 夏季，夏、秋两季空气质量空间分异格局不显著，春、冬两季的空间分异格局显著，呈北高南低、内陆高沿海低的季节性空间分布格局。程雪雁等 [ 9 ] 基于生态环境部发布的空气质量状况月报中京津冀城市群13市空气污染综合指数、主要污染物以及PM 2.5 、PM 10 、SO 2 、NO 2 、CO、O 3 等6种污染物质量浓度数据，对京津冀城市群2015~2018年的空气污染状况的时空分布进行了可视化呈现与特征总结、解释，研究发现：河北省东部的唐山、西部的保定、石家庄以及南部的邢台、邯郸等城市空气品质较差；京津冀城市群的主要污染物是PM 2.5 ，第2大主要污染物是PM 10 ，近几年来O 3 逐渐成为京津冀城市群夏季的主要污染物；此外，6大污染物中，PM 10 和SO 2 在冬季有明显的空间集聚现象。随着国家“一带一路”、长江经济带等重大发展战略和一系列重大政策的实施，云南工业化进程不断加快，2015年来，云南空气污染事件明显上升，其中最主要的是PM 2.5 和O 3 污染。本文根据2014年1月1日至2017年12月31日云南省昆明、曲靖和玉溪3个城市的空气污染监测数据，分析了三个城市的PM 2.5 和O 3 的时间变化特征，以期为科学有效地控制和治理大气污染提供科学参考。
system: [SYSTEM]

user: 2. 资料与分析方法 2.1. 数据处理 数据基于中国空气质量在线监测分析平台(https://www.aqistudy.cn/)，该平台每小时公布一次空气污染状况数据，包括PM 2.5 、PM 10 、CO、NO 2 、O 3 和SO 2 的小时平均值，公布的数据包括昆明、曲靖和玉溪等云南省内16个城市或州(图1)，由于各个点布设的时间不一致，本研究仅统计包含连续观测4年的3个典型城市(2014年1月1日至2017年12月31日)，并基于日平均值进行统计，需要说明的是，数据首先剔除了异常值0，并以日平均值统计分析导致云南空气污染最严重物质——PM 2.5 和O 3 的污染特征。 2.2. 研究方法 箱线图是一种可以简洁直观的表示数据分布特征的方法，用于反映一组或多组连续型定量数据分布的中心位置和散布范围。目前已见箱线图用于环境领域的分析研究，使复杂庞大的监测数据得以有规律的呈现出来 [ 10 ] [ 11 ] [ 12 ] 。本文用箱线图分析PM 2.5 、O 3 数据规律，并对比各城市污染状况。箱线图与描述统计中的最大值、上四分位数、中位数、下四分位数、最小值这5个统计量密切相关，本研究所作的箱线图，剔除了极值，统计的须子范围在5%~95%之间。 图1. 观测点分布图 HYSPLIT—Hybrid Single Particle Lagrange Integrated Trajectory Model (即拉格朗日混合单粒子轨道模型)，是由美国国家海洋和大气管理局(NOAA)的空气资源实验室和澳大利亚气象局在过去20年间联合研发的一种用于计算和分析大气污染物输送、扩散轨迹的专业模型。该模型具有处理多种气象要素输入场、多种物理过程和不同类型污染物排放源功能的较为完整的输送、扩散和沉降模式，已经被广泛地应用于多种污染物在各个地区的传输和扩散的研究中。采用HYSPLIT后向轨迹模式绘出云南省典型城市大气污染物输送，该模式结合了大气水平和垂直运动，在污染物输送、扩散和沉降等方面都有较完整的考虑，被广泛的运用于分析大气污染过程 [ 13 ] [ 14 ] [ 15 ] 。 气团对一个地区影响的路径主要与当地气压系统维持的位置及移动规律相关，也和地形地貌的特点及局地的环流密不可分。HYSPLIT模型的轨迹计算原理如下：假设一个被动示踪粒子随风运动，t时刻处于 P ( t ) 位置，时间步长 Δ t 后的位置 P ( t + Δ t ) 为(其中V为风速)： P ( t + Δ t ) = P ( t ) + 0.5 [ V ( P ， t ) + V ( P ′ ， t + Δ t ) ] Δ t 式中， P ′ 为一个中间猜想位置，计算过程中对时间步长 的要求是： P ( t + Δ t ) = P ( t ) + V ( P ， t ) Δ t Δ t ⋅ U max 小于0.75倍的网格距( U max 为最大风速，由气块上一时刻的最大迁移速度决定)。计算中如果轨迹点超出模型上边界则轨迹终止，如果轨迹点到达地面则沿地表继续运行 [ 16 ] 。
system: [SYSTEM]

user: 4. 结论 1) 2014~2017年云南省典型城市PM 2.5 浓度整体呈下降趋势，在我国处于较优水平，但同国外的一些国家相比，仍有一定差距。 2) PM 2.5 浓度水平受污染源排放和气象条件影响，存在明显的季节和时段变化特征，总体趋势大致表现为冬春季高，夏秋季低，浓度的变化与天气背景密切相关。 3) 2014~2016年昆明、曲靖O 3 -8h浓度呈逐年上升趋势，到了2017年情况好转，浓度明显下降；玉溪则相反，2014年玉溪O 3 -8h浓度较高，达到75.47 μg·m −3 ，2015年浓度显著下降，2015~2017年逐年增加。 4) 冬季两个典型污染观测期间采样点的气团后向轨迹显示，云南重污染天气的出现，可能与外来源叠加有密切关系，最可能的是云南西南部污染和来自缅甸等国家的长距离传输叠加导致。
system: [SYSTEM]

user: 摘 要 以典型冶炼企业的砷锑烟尘为研究对象，对砷、汞的赋存形态和分布特性进行研究。结果表明，砷锑烟尘中的重金属污染物主要包括砷、汞、镉、锑等，其中回转窑工段以砷尘为主。回转窑烟尘的砷含量高达48.00%、汞含量为0.005%，电炉烟尘的砷含量仅1.15%、汞含量为0.0064%。烟尘的分布规律表明，回转窑烟尘为粗颗粒烟尘，粒径大于10 μm的烟尘总量占34.83%，电炉烟尘主要集中分布小粒径范围内，1 μm以下的颗粒总量41.72%。砷锑烟尘中砷形态表以As(III)为主，汞含量较低，主要以单质汞和氧化汞的形态存在。 关键词 :砷锑烟尘，砷，汞，赋存形态，分布特征 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 砷锑烟尘是火法处理铅阳极泥的过程中产生的一种含砷高的烟尘，还含大量的锑、铅、镉等有价金属 [ 1 ] 。砷锑烟尘携带的砷含量较高，由此带来的大气污染和砷排放问题十分严重，已成为制约有色冶炼行业可持续发展的主要因素。对含砷物料处置不当，或在系统中循环累积，或转化成白砷产品，甚至直接堆存或排放，造成严重的安全与环境隐患，是导致我国砷污染事件频发的重要原因 [ 2 ] 。 砷是一种严重污染环境的剧毒物质，从砷锑烟尘中脱除砷并回收其他有价金属，是消除砷害，提高资源综合回收的有效途径。传统的挥发脱砷处理烟尘，砷以As 2 O 3 挥发，但产生的含砷烟气难以处理，造成严重的环境污染 [ 3 ] [ 4 ] [ 5 ] [ 6 ] 。国内冶炼厂处理高砷锑烟尘的传统工艺主要为反射炉还原挥发，但由于砷、锑物理化学性质相似，挥发的含砷烟尘中仍含有20%~50%作用的锑，烟气中特征物浓度居高不下，难以解决砷污染问题 [ 7 ] [ 8 ] 。此外，开发除砷的技术理论研究越来越受到重视 [ 9 ] [ 10 ] [ 11 ] ，但现有报道尚未对砷、锑的赋存形态及分布特征有实质性研究结论。 基于目前的研究现状，本论文开展对砷锑烟尘重点污染物砷和汞的分布特征和存在形态研究，明确砷锑烟尘的主要污染物，是解决冶炼系统砷汞去除特性、改善砷汞污染的重要理论，研究结果对我国有色金属行业发展的提供理论基础，也为改善环境质量，维护公众健康权益，促进经济、社会和环境的协调发展，具有重要的现实意义。
system: [SYSTEM]

user: 2. 原料与方法 2.1. 试验样品 本文以云南省某典型冶炼企业的传统的砷锑烟尘为研究对象，烟尘通过除尘器捕集，采集的试验样品取自电炉还原熔炼除尘系统和回转窑焙烧除尘系统。调查采样过程中，原料成分性质相同，冶炼工况稳定。 2.2. 分析方法 待测的烟尘样品在65℃下干燥12 h，其粒度分布采用激光粒度分析仪(型号：MS2000，英国马尔文公司)和振动筛分机(型号：SPB-200，柳州探矿机械厂)进行测定。采用电感耦合等离子体–原子发射光谱(ICP-AES，ICAP 6300，美国热电公司)检测分析烟尘样品主要元素的含量 [ 12 ] 。
system: [SYSTEM]

user: 3. 结果与讨论 3.1. 烟尘的重金属种类 为了研究冶炼企业所排砷锑烟尘中的重金属成分，并对其重点工段的污染源烟尘进行监测分析，研究所排烟尘的重金属成分。经对电炉还原工艺和回转窑工艺的除尘器烟尘进行采样监测，分析了烟尘中砷、汞、镉、锑重金属污染物的浓度，监测数据见表1。 Table 1 Heavy metal concentration of key pollution source in precipitator outlet/(% 取样点位 Sn Cd As Hg 回转窑烟尘 8.65 5.73 46.92 0.005 电炉烟尘 48.00 1.20 1.15 0.0064 表1. 主要砷锑烟尘的重金属浓度 对检测数据进行分析表明：各烟尘中均有砷、汞、镉、梯等重金属污染物，但其含量各有不同。其中回转窑工段以砷尘为主，电炉工段以锑尘为主，回转窑烟尘和电炉还原烟尘的砷含量分别达到48.00%、6.95%、1.15%，汞含量分别达到0.005%、0.0055%、0.0064%。 3.2. 烟尘粒径分布 实际工艺中，为了分离烟尘中的砷和锑，在烟气约310℃以上时采用电除尘回收锑尘，在温度低于100℃时采用袋式除尘回收砷尘，以利于砷、锑分离。根据湿式电除尘和布袋除尘捕集烟尘的粒径范围有明显的区别。因此，本论文为研究冶炼企业砷锑烟尘粒径分布规律，整理了相关粒径检测数据进行研究分析。 采用激光粒度分析以对三个样品的粒径进行分析，其累积分布率见图1。可以看出，不同冶炼工段下烟尘都符合累积质量分数先增大再减少的趋势，但不同工段下的烟尘粒径大小差异很大，电炉烟尘的粒径最小，同一粒径范围下的烟尘，电炉烟尘累积质量分数 > 回转窑烟尘累积质量分数。研究表明，粒径小于5 μm的颗粒能被吸入人体肺部深处，从而影响人体健康。因此，烟尘的粒径分布特征是作为判断砷锑烟尘是否进行高效回收和深度处理的必要依据。 图1. 烟尘粒径累积分布率 为进一步验证工业实际工况，采用振动筛分机对烟尘颗粒的表观粒度分布进行分析，其粒径分布统计表见表2。回转窑烟尘粒径主要分布在1~3 μm、3~5 μm、5~7 μm、7~10 μm的粒径范围内，其分布率分别为10.77%、14.95%、19.01%和30.37%。回转窑烟尘表现为粗颗粒烟尘，≥10 μm的烟尘总量为34.83%。电炉烟尘粒径主要集中分布在<0.5 μm、0.5~1 μm、1~3 μm、3~5 μm的粒径范围内，其分布率分别为24.03%、17.69%、37.04%和12.50%。电炉烟尘表现为细颗粒，3 μm内的颗粒总量78.76%，1 μm以下的颗粒总量41.72%。 Table 2 Statistical table of different particle size of dus 烟尘类型 各粒径范围内的烟尘体积分数(%) <0.5 μm 0.5~1 μm 1~3 μm 3~5 μm 5~7 μm 7~10 μm 10~20 μm >20 μm 回转窑烟尘 1.75 3.84 10.77 14.95 14.85 19.01 30.37 4.46 电炉烟尘 24.03 17.69 37.04 12.50 2.39 3.74 2.32 0.29 表2. 烟尘粒径分布统计表 3.3. 烟尘砷汞的分布特征 烟尘的粒径直接决定其逸散难易程度，采用振动筛分机将烟尘颗粒进行分级处理，按环保PM 10 、PM 2.5 的划分标准将烟尘分级为≥10 μm、2.5~10 μm、1~2.5 μm、≤1 μm。回转窑烟尘和电炉烟尘的重金属含量分布表分别见表3、表4，锑富集在粗颗粒。烟尘分级处理后，砷含量随粒径增大而减小，锑含量随粒径增大而增大，汞的含量随粒径变化的影响不明显。粒径>10 μm回转窑烟尘烟尘占34.83%，锑含量增大至18.90%；砷主要富集在细颗粒，PM 10 的砷含量明显增大。电炉烟尘主要富含锑，粗颗粒的砷含量显著降低，分离出PM 2.5 后实现粗粒中砷的去除，有利于锑的回收。 Table 3 Statistical table of different heavy metals of rotary-kiln dus 类别 粒径/μm 质量分数/% As/% Hg/% Sb/% 1 ≥10 34.83 35.20 0.0052 18.90 2 2.5~0 51.20 47.10 0.0053 7.34 3 1~.5 8.38 50.23 0.0050 4.38 4 ≤1 5.59 53.12 0.0047 2.37 合计 / 100 表3. 回转窑烟尘重金属含量分布表 Table 4 Statistical table of different heavy metals of furnace dus 类别 粒径/μm 质量分数/% As/% Hg/% Sn/% 1 ≥10 2.61 0.67 0.0052 58.9 2 2.5~10 23.27 0.92 0.0053 52.3 3 1~2.5 32.4 1.03 0.0050 49.31 4 ≤1 41.72 1.26 0.0047 43.07 合计 / 100 表4. 电炉烟尘重金属含量分布表 3.4. 烟尘砷汞的形态研究 XPS对重金属的价态监测有优良的灵敏度，能分析颗粒表面约10 nm深度的重金属赋存形态。对三种烟尘样品进行XPS扫描，分析其颗粒表面砷和汞的定性分析，砷的XPS精细扫描的拟合结果见图2。由于汞含量太低，难以通过XPS图谱的峰值进行分析，其价态分析结果见表5。 图2. 冶炼烟尘XPS图谱 由图2可知，As 3d的谱峰对应的XPS结合能较高，属于复杂As配合物中的As-O键。三个样品中砷的谱峰变化不大，表明其物相组成无明显变化，均为砷的氧化态。烟尘中汞主要以单质和化合物存在，单质汞含量共占45.9%；汞的化合物主要为氧化汞，占49.5%，还检测到少量氯化汞。 Table 5 XPS analyze of mercur 元素峰 化合物 结合能/eV 百分率/% Hg4f 7/2 HgCl 2 101.400 4.60 Hg4f 5/2 Hg 0 104.000 45.90 Hg4f 5/2 HgO 104.800 49.50 表5. 汞的XPS分析 3.5. SEM分析 对烟尘分级处理后，采用XL30型环境扫描电镜对各烟尘进行分析，SEM图谱见图3，作为冶炼烟尘，其自然颗粒很细，胶结和相互嵌布很密。对烟尘样品对比分析表明，细烟尘和粗烟尘在微观结构有明显的不同，大颗粒表面粘结的微粒脱落，晶粒均匀，与其成分相关，符合砷汞赋存形态的特征。
system: [SYSTEM]

user: 4. 结论 1) 冶炼企业排放砷锑烟尘中的重金属污染物主要包括砷、汞、镉、梯等，其中回转窑工段以砷尘为主，回转窑烟尘和电炉还原烟尘的砷含量分别达到48.00%、1.15%，汞含量分别达到0.005%、0.0064%。 图3. 不同烟尘的扫描电镜形貌 2) 烟尘的分布规律表明，回转窑烟尘表现为粗颗粒烟尘，粒径大于10 μm的烟尘总量占34.83%，电炉烟尘主要集中分布小粒径范围内，1 μm以下的颗粒总量41.72%。 3) 砷锑烟尘中砷形态表现为As(Ⅲ)为主，汞的含量较低，主要以单质汞和氧化汞的形态存在。
system: [SYSTEM]

user: 摘 要 为了扩大处理规模及提高出水水质，对湖北某白酒公司产生的白酒废水进行工程设计。采用预处理-UASB-两级A/O-芬顿工艺对白酒废水进行深度处理，并介绍了主要构筑物及相关设计参数。实践结果表明：色度 < 10，COD < 30 mg/L，SS < 25 mg/L，氨氮 < 4 mg/L，总氮氨氮 < 19 mg/L，总磷氨氮 < 0.9 mg/L，出水水质优于《发酵酒精和白酒工业水污染物排放标准》(GB27631-2011)直接外排标准，满足设计要求。 关键词 :白酒废水，深度处理，UASB，两级A/O，芬顿反应器 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 白酒酿造产生的废水是具有高浓度、高温且呈酸性的有机废水 [ 1 ]，主要污染物为淀粉、糖类、蛋白质、醇类等物质，具有COD浓度高、pH值低、SS高、氨氮和磷浓度高等特点 [ 1 ]。可生化性好，具有高污染性和潜在的资源可利用性 [ 2 ]。针对白酒行业废水的特点，一般采用“分别收集、分级削减，资源回收削减污染负荷”的治理技术路线 [ 3 ]，采用厌氧–好氧组合工艺进行处理 [ 4 ]，主要采用A2/O工艺 [ 5 ]、SBR工艺 [ 6 ]、氧化沟 [ 7 ]、MBR工艺 [ 8 ]、相关改进工艺及联合处理工艺等脱氮除磷 [ 9 ]。目前，我国大型白酒企业已具备成熟的污水处理工艺，出水水质达到《发酵酒精和白酒工业水污染物排放标准》(GB 27631-2011)一级排放标准 [ 10 ]。 湖北某白酒公司是一家大曲清香型白酒酿造企业，采用固态发酵工艺，以高粱、青稞、大麦、稻谷为原料，年产原酒5000吨。由于公司生产规模扩大后，废水量增大，污染物浓度升高，超出原污水处理系统处理能力，且国家标准《发酵酒精和白酒工业水污染物排放标准》(GB27631-2011)对直排废水作出了更严格的要求，急需对现有处理系统进行升级改造。
system: [SYSTEM]

user: 2. 工程设计 2.1. 进出水水质与水量 白酒酿造产生的废水主要为蒸馏锅底废水、发酵阶段黄水、原料清洗水、车间清洗废水及部分生活污水。生产不受时间和季节性影响，但每日排放的废水水质、水量变化较大。该废水有机物浓度高，温度高，悬浮物含量大，色度高，酸度大的特点，属于典型的高浓度有机废水 [ 11 ]。按业主要求，设计处理能力为400 m 3 /d。具体设计进、出水质指标见表1，处理出水需达《发酵酒精和白酒工业水污染物排放标准》(GB27631-2011)直接排放标准及公司内部要求。 Table 1 Designing inlet and outlet water quality indicator 项目 pH 色度 COD/(mg/L) BOD5/(mg/L) SS/(mg/L) 氨氮/(mg/L) 总氮/(mg/L) 总磷/(mg/L) 进水指标 3.5~5 ≤1000 ≤18，000 ≤10，000 ≤3000 ≤120 ≤400 ≤400 出水指标 6~9 10 30 6 25 4 19 0.9 表1. 设计进、出水水质指标 2.2. 设计要点 本工程主要包括预处理、生化处理、深度处理等阶段。出水水质要求高，常规二级处理工艺无法保证水质稳定达标，需设置深度处理工艺。在工艺设计时应考虑的因素： 1) 原水中含有稻谷壳、粉碎秸秆等粗纤维，生化处理前应设置粗格栅和机械转筛将其去除，以保证后续构筑物及设备的正常运行。 2) 设置调节池，解决进水水质及水量不均衡的问题。 3) 原水COD浓度较高，需进行厌氧处理以降低负荷，采用UASB厌氧塔一是降解90% COD及25%总氮、50%总磷，保障后续两级A/O处理工艺的稳定运行，二是可以产生沼气，回收资源利用。 4) 两级A/O具有脱氮除磷处理效果好，抗负荷冲击能力强，出水水质比较稳定的优点，用以保障废水的处理效果 [ 12 ] [ 13 ]。 5) 深度处理段采用芬顿反应–混凝工艺处理残余COD、总氮、总磷及色度，具有处理效率高，稳定性好，成本低等特点，进一步提高出水水质 [ 14 ]。 2.3. 废水处理工艺 工艺流程图如图1所示。 图1. 工艺流程图 废水首先由粗格栅去除碎秸秆及部分稻谷壳和酒糟后进入集水井，经泵提升至转筛过滤器去除剩余稻谷壳、酒糟，出水自流入调节池，匀质匀量，保证生化处理段的稳定运行。调节池出水经泵提升至UASB厌氧塔，同时设置超越管，可将调节池出水定量流入至缺氧池1与缺氧池2，废水在UASB厌氧塔内完成大部分有机污染物的厌氧降解，有机氮转化为无机氮，出水自流入缺氧池1，与超越管出水、二沉池回流污泥及硝化液回流液混合后进入反应区，采用潜水搅拌机混合搅拌，微生物在缺氧池1内利用废水中碳源反硝化回流污泥与硝化液回流中硝态氮与亚硝态氮，同时部分COD被降解，并产生一定的碱度。出水进入好氧池1，废水在好氧池1中大部分COD被降解同时完成氨氮的硝化反应，反应区末端硝化液回流至缺氧池1，回流比为200%~400%。由于池内微生物浓度较高，微生物系统存在明显的微观好氧与缺氧环境，在好氧池1内发生同时硝化与反硝化、短程硝化与反硝化反应，减少系统的反硝化的碳源需求量、耗氧量与碱度补充量。好氧池1出水自流入缺氧池2，缺氧池2根据情况补加碳源，与废水混合后，微生物利用外加碳源快速反硝化，继续降解废水残余总氮，出水进入好氧池2。废水在好氧池2内残留的有机污染物或氨氮继续被降解，保证出水水质，好氧池2出水进入二沉池，泥水分离后进入后处理单元，池内污泥回流至缺氧池1，污泥回流量100%~200%，剩余污泥排至污泥浓缩池。中间水池出水提升至高效芬顿反应器，在反应器内利用H 2 O 2 的强氧化性与硫酸亚铁的催化作用氧化废水中残余COD、色度，并去除大部分总磷，出水进入后反应池，利用未反应完全的氧化剂继续降解废水有机污染物，并投加PH调节剂与助凝剂，在混凝沉淀池内进行泥水分离后达标排放。
system: [SYSTEM]

user: 3. 主要构筑物 3.1. 调节池 调节池2座(1用1备，可以临时作为应急水池)，地下封闭式砼结构，规格12 m × 10 m × 4 m，每座有效容积420 m 3 ，HRT = 25 h。配备排污泵2套(1用1备)，Q = 20 m³/h，H = 10 m，N = 1.5 KW。 3.2. UASB厌氧塔 厌氧池为2个规格为φ8 m × 10 m的罐体，内壁防腐，总有效容积900 m 3 ，HRT = 54 h，容积负荷=8 kgCOD/(m 3 ∙d)。配备布水系统2套，pH监测仪2套，三相分离器2套，温控仪2套，沼气收集装置2套。 3.3. 两级A/O 一级A/O，通过内回流将好氧池1中产生的硝态氮输送至缺氧池1，利用原污水中碳源进行反硝化作用，强化脱氮能力 [ 15 ] [ 16 ] ；二级A/O，缺氧池2去除剩余硝态氮，解决了一级A/O工艺脱氮能力有限的问题，也使工艺具有更强的抗冲击负荷能力；好氧池2利用聚磷菌在好氧环境下吸磷的作用，强化了除磷功能 [ 16 ]。缺氧池1为规格φ6 m × 9 m的罐体，总有效容积200 m 3 ，HRT = 12 h，容积负荷=1.4 kgCOD/(m 3 ∙d)，DO控制在0.2~0.5 mg/L；好氧池1由规格6 m × 4 m × 4.5 m的池子构成，总有效容积96 m 3 ，HRT = 32 h，容积负荷=0.7 kgCOD/(m 3 ∙d)，DO控制在2~4 mg/L，回流比200%~400%；缺氧池2由1个规格3.75 m × 4 m × 4.5 m的池子构成，总有效容积60 m 3 ，HRT = 4 h，容积负荷=0.6 kgCOD/(m 3 ∙d)；好氧池2由1座规格为2.25 m × 4 m × 4.5 m的池子构成，HRT = 2.4 h，污泥负荷=0.3kgCOD/(m 3 ∙d)；配套设备：缺氧池1配备潜水搅拌机2套(1用1备)，功率为1.5 kw，叶轮直径为260 mm，转速为960 r/min；好氧池1配备硝化液回流泵2套(1用1备)，Q = 45m³/h，H = 10 m，N = 4 KW，加装微孔曝气盘300个，配备DO和pH监测仪各1套；缺氧池2配备潜水搅拌机1套，功率为0.85 kw，叶轮直径为260 mm，转速为740 r/min；好氧池2加装微孔曝气盘22个；好氧系统由2台罗茨鼓风机(1用1备) Q = 20 m³/min，H = 5 m，N = 30 KW供氧。 3.4. 二沉池 二沉池规格为φ8 m × 5 m，总有效容积220 m 3 ，HRT = 12 h，水力负荷为0.36 m 3 /(m 2 ∙h)，污泥回流比=100%~200%，每天排泥1~2次。配备2台污泥回流泵(1用1备)，规格为Q = 20 m³/h，H = 10 m，N = 1.5 KW。 3.5. 中间水池 中间水池由1座规格为5 m × 2 m × 2.4 m的池子构成，总有效容积20 m 3 ，HRT = 2 h，水力负荷=1 m 3 /(m 2 ∙h)，配套设备：规格为Q = 20 m³/h，H = 10 m，N = 1.5 KW的潜水提升泵2套(1用1备)，电磁流量计和电缆浮球液开关各1套。 3.6. 芬顿反应器 芬顿反应器的规格为φ1.5 m × 4 m，总有效容积5 m³，反应器内设置催化填料层，通过加入H 2 O 2 和硫酸亚铁，利用H 2 O 2 的强氧化性与硫酸亚铁的催化作用氧化废水中残余COD、色度，从而保证出水达到设计要求。 3.7. 混凝沉淀池 混凝沉淀池由1座规格为6 m × 4 m × 4.5 m的池子构成，总有效容积96 m³，HRT = 6 h，水力负荷=0.7 m 3 /(m 2 ∙h)，配备2套(1用1备)规格为Q = 10 m³/h，H = 10 m，N = 0.75 KW的排泥泵。 3.8. 污泥浓缩池 污泥浓缩池由1座规格为6 m × 3 m × 4.5 m的池子构成，总有效容积72 m³，污泥浓缩时间12 h，配备厢式压滤机1套，规格为Q = 10 m³/h，P = 10 Bar，N = 5.5 kw，转速258 rpm。
system: [SYSTEM]

user: 4. 运行处理效果 该污水处理系统自2018年11月建设完成，经过8个月的调试，目前已经稳定运行。出水水质已经达到设计要求。出水监测结果如表2所示。 Table 2 Effluent water quality indicator 项目 pH 色度 COD/(mg/L) BOD5/(mg/L) SS/(mg/L) 氨氮/(mg/L) 总氮/(mg/L) 总磷/(mg/L) 2019.7.15 7.7 4 17.4 — 14 1.9 18.1 0.7 2019.7.16 7.3 4 18.4 — 14 1.8 18.1 0.5 2019.7.17 7.6 4 25.5 — 11 1.3 17.8 0.4 2019.7.18 7.5 4 24.8 — 17 1.3 17.5 0.6 2019.7.19 7.4 4 24.5 — 17 1.2 17.4 0.3 平均值 7.5 4 22.1 — 14.6 1.5 17.8 0.5 设计标准 6~9 10 30 6 25 4 19 0.9 国家标准 6~9 40 100 30 50 10 20 1 表2. 出水水质指标 由表2可知，处理后色度、COD、BOD、SS、氨氮、总氮、总磷的总去除率分别为99.6%、99.9%、99.9%、99.5%、96.8%、95.6%、99.9%，出水水质全部达到了《发酵酒精和白酒工业水污染物排放标准》(GB27631-2011)直接外排标准及设计的指标要求。
system: [SYSTEM]

user: 5. 运行费用 废水处理运行费用如表3所示。 Table 3 Wastewater treatment operating cost 项目 吨水费用/(元/吨) 吨水成本 用水 0.1 6.88元/吨 用电 2.27 (按工业用电计算) PAM 0.25 硫酸亚铁 2.1 液碱 1 双氧水 1.3 表3. 废水处理运行成本 该项目总投资200万元，包括基建和设备费用.污水处理系统年运行总成本为90.8万元，吨水运行成本为6.88元/吨，年生产天数330天左右，计算水量为400 t/d。本运行费用为直接费用，不包括人工费、设备折旧费、维修费、企业管理费等。
system: [SYSTEM]

user: 6. 结论 1) 白酒废水经过预处理-UASB-两级A/O-芬顿工艺进行处理后，出水水质优于《发酵酒精和白酒工业水污染物排放标准》(GB27631-2011)直接外排标准，达到了设计的目标。 2) UASB工艺，大大降低了原水污染负荷，为后续处理及脱氮除磷创造了良好的条件。 3) 两级A/O工艺对水质的达标排放起到了关键性作用，该工艺处理效果好，抗负荷冲击能力强，出水水质稳定。 4) 芬顿–混凝沉淀工艺进一步提高水质，保障了出水达到设计要求。最终出水结果：色度 < 10，COD < 30 mg/L，SS < 25 mg/L，氨氮 < 4 mg/L，总氮氨氮 < 19 mg/L，总磷氨氮 < 0.9 mg/L。 5) 本工艺比较适合中小规模的白酒生产企业的废水处理，对没有市政排水管网及排放要求严格的地区，具有一定的借鉴意义。
system: [SYSTEM]

user: 摘 要 目的：评价经肝动脉化疗栓塞(TACE)联合调强放疗(IMRT)治疗复发性或残留性原发性肝癌(PHC)的临床疗效和毒性反应。方法：39例接受两次及以上TACE治疗复发或残留的PHC患者，19例进行调强放疗(IMRT)，放疗剂量50~60 Gy/25~30F；20例TACE维持治疗。比较两组治疗前后肿瘤大小和甲胎蛋白(AFP)水平变化、毒副反应及生存期。结果：至随访时IMRT组和TACE组的有效率(CR + PR)和中位生存时间分别为63.2% vs 25.0%和20.2个月vs 13.7个月，差异有统计学意义(P = 0.016和P = 0.021)，IMRT组AFP下降水平明显优于TACE组，不良反应I、II级肝功能损伤高于TACE组。结论：IMRT对TACE治疗后局部复发或残留的PHC患者具有良好的临床疗效，毒性反应可以耐受。 关键词 :原发性肝癌，甲胎蛋白，肝动脉栓塞，放射治疗 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 原发性肝癌(PHC)作为消化系统常见的恶性肿瘤之一，早期PHC的主要治疗手段是外科治疗，早期患者由于临床症状隐匿难以发现，确诊时往往多为中晚期，对于中晚期PHC最主要的治疗方法是肝动脉化学栓塞(TACE) [ 1 ] 。单纯的TACE可以将药物混合栓塞剂碘化油直接注入肿瘤组织，阻断肿瘤的动脉血供并增加肿瘤组织内药物浓度。由于肝脏的供血具有肝动脉、侧枝动脉及门静脉双重供血的特点，即使肿瘤肝动脉被TACE完全栓塞，侧枝动脉及门静脉仍可能为残留的肿瘤细胞提供血液 [ 2 ] ，从而为以后的复发和转移提供了血供来源。因此，本研究旨在评估IMRT治疗TACE术后复发或残留HCC的临床疗效和毒副反应。
system: [SYSTEM]

user: 2. 材料和方法 2.1. 临床资料 2014年6月至2018年7月，TACE治疗后局部复发或残留原发性肝癌患者39例，平均年龄为55岁(42~75岁)，男性24例、女性15例，19例采取IMRT治疗，20例TACE维持治疗。入选标准如下：1) 病理诊断为肝细胞癌11例，经AASLD标准 [ 3 ] 临床诊断28例；2) 接受TACE至少两次以上；3) 复发和残留诊断标准：在TACE治疗后CT/MR检查，发现肿瘤区碘化油沉积欠佳或肿瘤边缘增长现象定义为肿瘤复发或残余；4) 肝功能Child-Pugh评分A或B级；5) KPS ≥ 80分；6) 正常肝组织体积 > 700 ml；7) 无浸润或淋巴结转移；8) 无腹水征象。在治疗过程中，排除了接受肝肿瘤手术、使用靶向药物或全身化疗的患者。所有患者入组前均签署治疗知情同意书，本研究获得新沂市人民医院伦理委员会批准。 2.2. 治疗方法 放射治疗患者仰卧位，负压真空袋联合体膜固定体位，螺旋CT 1.25 mm层距连续增强扫描定位；飞利浦Pinnacle 3 治疗计划系统进行计划设计，由两名放射肿瘤科医师及物理师共同勾画靶区(GTV)和危及器官(OARs)，以GTV外扩0.5~1 cm为计划靶体积(PTV)的勾画范围，95%的等剂量曲线覆盖靶体积(PTV)，剂量体积直方图(DVH)进行评估并优化放射治疗方案，传输到Elekta公司Precise医用直线加速器按照设定的程序进行治疗。放疗处方剂量：2 Gy/次，共25~30次，5次/W，总剂量50~60 Gy。靶区邻近器官不超过其耐受剂量，平均肝脏剂量 < 23 Gy；接受>30 Gy剂量照射的正常肝脏体积(V30)达到最小化(均<30%)，治疗期间给予常规保肝治疗。 2.3. 评价标准 客观疗效评价按照RECIST1.1 [ 4 ] ，分为完全缓解(CR)、部分缓解(PR)、稳定(SD)和进展(PD)，以CR + PR为有效(RR)，CR + PR + SD为获益(DCR)。评价入组治疗前及在IMRT和TACE治疗后3个月AFP水平和肿块大小变化情况。毒副反应按照美国NCICTC标准评价，分为0~IV级。 2.4. 统计方法 采用SPSS 17.0统计学软件进行数据分析及处理，Kaplan-Meier法描绘生存曲线、Log-rank检验比较两组间总生存期(Overall survival OS)的差异，百分比表示客观疗效和毒副反应的发生率，组间比较为X 2 检验，所有P < 0.05为差异有统计学意义。
system: [SYSTEM]

user: 3. 结果 3.1. 疗效 IMRT组：有效率为63.2% (12/19)高于TACE组的25.0% (5/20)，两组差异有显著性统计学意义(P = 0.016)；IMRT组：11例患者AFP下降到正常范围及治疗前一半，明显优于单纯TACE组25.0% (5/20) (P = 0.037) (表1所示)。 Table 1 Comparison of the effect between two group 治疗方法 例数 CR PR SD PD AFP ( 下降至或近正常 ) RR P (RR) P (AFP) 研究组 19 1 11 4 3 11 63.2% (12/19) 0.016 0.037 对照组 20 0 5 7 8 5 25.0% (5/20) 表1. 两组疗效的对比(%) 3.2. 毒副反应 IMRT组消化道反应6例低于TACE组9例，IMRT组有轻度乏力9例、I和II级肝功能损伤15例多于对照组(表2所示)，两组间差异无显著统计学意义(P = 0.06)。 Table 2 Toxic and side effects in two group 毒副反应 研究组 合计(%) 对照组 合计(%) 0 I II III 0 I II III 骨髓抑制 9 7 3 0 26.3 (5/19) 14 5 1 0 35.0 (7/20) 消化道反应 13 6 0 0 21.1 (4/19) 11 6 3 0 45.0 (9/20) 肝功能损伤 4 10 5 0 78.9 (15/19) * 10 7 3 0 50 (10/20) * 肾功能损伤 19 0 0 0 0 (0/19) 20 0 0 0 0 (0/20) 乏力 10 6 3 0 47.4 (9/19) 20 0 0 0 0 (0/20) 表2. 两组的毒副反应(%) * P = 0.060. 3.3. 生存分析 至随访日研究组和对照组的中位生存时间分别为20.2个月(95% CI：16.4~24.1个月)和13.7个月(95% CI：10.6~16.7个月)，Log Rank检验示两组间P = 0.021，差异有统计学意义(图1所示)；治疗前后AFP水平下降明显的患者生存时间优于下降不显著的患者P = 0.076，差异无显著统计学意义(图2所示)。 图1. 两组间生存时间曲线(P = 0.021) 图2. AFP下降水平生存时间曲线(P = 0.076)
system: [SYSTEM]

user: 4. 讨论 肝动脉化学栓塞(TACE)是中晚期原发性肝癌(PHC)主要的治疗手段之一，有研究报道 [ 5 ] ，单纯TACE治疗PHC并不能完全控制肿瘤组织，栓塞后的病理学检查显示大约只有22%~50%的肿瘤组织被完全破坏。另外由于肝癌的血供来源于肝动脉、侧枝动脉及门静脉多重供血，肿块直径较大的肝癌患者经TACE后，肿块边缘的仍有侧枝动脉及门静脉参与供血，在一定程度上影响了TACE 的治疗效果，也是单纯TACE治疗后肿块周围继续增长的主要原因之一。对于残余癌或肝癌复发的治疗没有统一的标准，文献 [ 6 ] [ 7 ] 报道了放射治疗、射频消融、冷冻、微波、酒精等方法，取得的临床效果并不是很满意。有报道 [ 8 ] [ 9 ] [ 10 ] SBRT联合TACE对原发性小肝癌治疗优于单纯TACE治疗组，且毒副反应轻微。 我们对于位于近肝脏边缘处原发性肝癌(最大直径8.5 cm)，在两次或两次以上TACE治疗后影像学检查发现，肿瘤栓塞不完全或肿块边缘有进展的患者可以采用调强放疗(IMRT)作为补充治疗手段。在我们的研究中，TACE联合IMRT组，12、18和24个月生存率分别为78.9% (15/19)、63.2% (12/19)和31.6% (6/19)，中位总生存期为20.2个月。单纯TACE治疗，12、18和24个月生存率分别为50.0% (10/20)、25.0% (5/20)和10.0% (2/20)，中位总生存期为13.7个月。我们还发现在所有患者中治疗前后AFP水平下降显著的患者生存时间要优于下降不显著的患者，差异无显著统计学意义可能与入组的样本量少有关系。在毒副反应方面和相关研究报道 [ 11 ] [ 12 ] 类似，IMRT组中有9名患者(47.4%)出现1~2级短暂疲劳，6名患者(31.6%)在1个月内出现1~2级胃肠道反应，对症支持治疗症状得到缓解或消除后，所有患者在随访期间未发现因IMRT引起的3级或更高的急性毒性。 综上所述，IMRT对原发性肝癌是一种的姑息性无创治疗方法，对TACE术后复发或残留的肝癌疗效显著，毒副反应患者可接受。在IMRT治疗后AFP水平下降明显，可能是生存预后的预测因素。
system: [SYSTEM]

user: 摘 要 为有效解决光学临近修正中的异常图片处理效率低的问题，本文建立一种可应用于芯片制造光学临近修正过程中的异常图像识别方法，能实现自动化的异常图像识别和异常数据过滤，从而有效减少了工程师人工经验干预，既保障了数据采集的有效性。 关键词 :OPC，异常图像，数据，芯片制造 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 概述 集成电路的发展特征是逻辑尺寸不断减小，光刻机波长是制约核心技术发展的决定因素，由于光学临近效应带来的光罩版图图形畸变也是制约集成电路发展的主要障碍。这种由于光波衍射、干涉而使光刻图形与掩模图形产生偏差的现象称为光学临近效应(OPE: Optical Proximity Effect)。在光刻工艺中，现有方法是采用光学临近修正(OPC: Optical Proximity Correction)方法来尽可能地减小掩模版图形投射到硅片图形的变形与偏差，使得曝光后的图形符合设计要求 [ 1 ]。 归纳起来，OPC技术的发展经历了：基于规则的OPC修正方法、基于模型的OPC修正方法、曝光辅助图形的应用、光源优化与掩模板优化、反演光刻技术5个发展过程。基于规则的OPC方法，其实现修正的关键是制定规则库。规则库设计的要求是要能够包含原始版图上所要求的所有图形的修正方案，随着技术节点的进一步缩小，基于规则的OPC方法就显得适用性较差 [ 2 ] ；基于模型的光学邻近修正的关键是建立精确的光学模型和光刻胶模型，因为每一层掩模版都有上千万个图形，因此不仅要求模型精度高，而且要求计算速度快 [ 3 ] ；因原始版图中包含密集分布图形和稀疏图形，密集分布图形的光刻工艺窗口和稀疏图形的光刻工艺窗口是不一样的，这就会导致共同工艺窗口的变小。适用于密集图形的光照条件并不适合稀疏图形的曝光，因此当技术节点发展到90 nm时引入曝光辅助图形(sub-resolution assistant feature， SRAF)来解决这一难题，在实际光刻中，所加入的曝光辅助图形对光线只起到散射作用，不影响实际的晶圆表面光刻成像，因此，光学邻近修正发展成为SRAF + OPC修正的模式 [ 4 ]。光源优化与掩模板优化通过截取原始版图上的关键图形分别进行SMO优化，通过SMO优化的结果进行OPC修正，并用于光刻机的光照条件整理，而后将修正好的掩模版进行制造；反演光刻技术与前面所说的OPC技术的目的一样，但其实现方法完全不同。它主要是以晶圆曝光后得到的图形为目标，反演计算出修正掩模版上所需要的图形 [ 5 ]。 OPC一直是芯片制造的核心。制约OPC发展的关键是OPC模型的建立与模型的验证，而数据量测是OPC建模和模型验证的基础。目前普遍存在的问题是：量测机台有时不能准确对准、量测条件容易设置错误，目标量测图形未能准确曝光等，导致量测结果都需要工程师结合人工经验进行检查和校正，这不仅在很大程度上加重了工程师的工作负荷。 为有效解决光学临近修正中的异常数据带来的图片处理效率低的问题，本文建立了一种可应用于芯片制造光学临近修正过程中的异常图像识别方法，通过该方法能实现自动化的异常图像识别，从而有效减少了工程师人工经验干预，保障了数据采集的有效性。
system: [SYSTEM]

user: 2. 技术方法 OPC数据采集过程中的异常图像处理方法，包括：真彩色图形转换、灰度图生成和灰度图均衡化、图像的灰度直方图获取、相邻图像的图形相关系数计算等步骤 [ 6 ] [ 7 ] [ 8 ]，能实现自动化的异常图像识别和异常数据过滤，功能模块图如图1所示。 图1. 功能模块图 2.1. 真彩色图转换 机台上原始采集的量测图片为红色、绿色、蓝色为基色的真彩色图片，因此在对图片进行预处理转换时，首先要对其进行转换，将其转换级灰度图像，转换公式如下： (1) 公式1中，R、G、B分别为每个像素点对应的RGB颜色分量值，Y为对应该点的灰度值。 2.2. 灰度图生成与均衡化 ①、灰度图生成 经过原始图像->灰度图像转换以后，便可得到变化是连续的灰度图 [ 9 ]。 ②、灰度图均衡化 灰度图均衡化的目的是筛选图片本征，便于计算像素点。通过像素点的灰度计算，便于后续的灰度相关性因素计算。根据之前所生成的灰度图，故可以用下面的公式对其进行均衡： D B = f ( x ) = D M a x ⋅ ∫ H ( u ) d u / A 0   ( u = 0 ~ x ) (2) 公式2中，x为灰度级， H ( u ) 为第i级灰度的像素个数，A0为图像的面积，DMax为最大灰度值。 2.3. 灰度直方图生成 在数字图像处理中，需用灰度直方图带表征灰度图中各灰度级出现的频率。所绘制的灰度直方图如图2，图3所示 [ 8 ]。 图2. 异常图像和灰度直方图 图3. 正常图像及灰度直方图 2.4. 图像相关系数求取与统计分析 异常图像筛选主要是通过公式(3)分别计算两幅图像的相关程度，来确定图像的相似程度，来达到筛除异常图像的方法。 异常图像的图像相关系数定义为R，用于表征2幅图片中的灰度相关属性，其计算公式如下： R = f ( t ) ⋅ g ( − t ) = ∫ f ( t ) ⋅ g ( t ) d t   ( − ∞ < t < + ∞ ) (3) 其中，函数 f ( t ) 和 g ( t ) 分别是要比较的两幅图像的某种相同属性，R值越大，说明两幅图像的相似程度越高，如图4所示。 图4. 灰度直方图对比 求取的R值结果：0.44，鉴定结果是两幅图片相关度较低。 2.5. 异常图像去除 采用上述方法，对所取得的10张量测图片进行两两比较，根据公式3可得到9个相关性系数R值，如表所示，根据统计学里的四分位数原理，将数列中的小于下四分位数(前三个点对应的图片)删除，即可实现异常图像的过滤，R值计算结果如表1所示，异常图像筛选如图5所示。 Table 1 R value resul 0.2 0.24 0.22 0.79 0.82 0.84 0.89 0.94 0.95 0.96 表1. R值计算结果 图5. 被过滤的异常图像
system: [SYSTEM]

user: 3. 结束语 异常图像处理效率低，且会对OPC相关工作产生重要影响，为有效解决该问题，在图像识别理论的基础上，本文建立一种可应用于芯片制造光学临近修正过程中的异常图像识别方法，该方法通过真彩色图形转换、灰度图生成和灰度图均衡化、图像的灰度直方图获取、相邻图像的图形相关系数计算等步骤，能实现自动化的异常图像识别和异常数据过滤，从而有效减少了工程师人工经验干预，保障了数据采集的有效性。
system: [SYSTEM]

user: 摘 要 假设对一般线性随机效应模型添加新的变量项成为过参数线性随机效应模型，这时两模型相同未知参数的统计推断不一定相同。针对这一问题，本文利用带约束限制二次矩阵值函数最优化问题的解，给出过参数线性随机效应模型下未知参数函数的最佳线性无偏预测/最佳线性无偏估计的解析表达式，并利用一些代数与矩阵理论工具，得到了两模型未知参数函数最佳线性无偏预测/最佳线性无偏估计等价的条件。 关键词 :线性随机效应模型，最佳线性无偏预测，最佳线性无偏估计，协方差矩阵，等价性 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 在做统计分析时，经常需要去分析两个或多个不同线性模型间的相似性或等价性。例如当回归模型出现过度参数化或过度拟合时，尽管原模型与新模型的响应变量相同，但两模型中相同未知参数的统计推断不一定相同，因此就有兴趣去研究这两个模型间的关系。关于这类问题的已做了的研究有，文献 [ 1 ] 研究了一般线性模型下的两个转换模型参数函数最佳线性无偏预测(BLUP)/最佳线性无偏估计(BLUE)间的关系。文献 [ 2 ] 考虑原线性模型添加新的回归因子表示了成过参数模型的形式，并研究了这两个模型中所有未知参数函数BLUP/BLUE间的关系。作为一般线性模型与过参数线性模型的更进一步研究，文献 [ 3 ] 考虑对一般线性模型添加新的参数或删除已有的参数，给出了具有固定系数和相同参数限制的两个具有竞争关系的线性模型，并研究了两模型未知参数函数BLUP间的关系。无论一般线性模型是否出现过参数或欠参数情况，都可以将这类模型看作是原模型的错误指定模型，文献 [ 4 ] 对原模型与错误指定模型间的关系进行研究，并利用一些代数与矩阵理论工具来刻画两模型下未知参数函数BLUP间的关系等。 在线性随机效应模型下做了的研究有，文献 [ 5 ] 利用线性代数和矩阵理论中最新的公式，给出两个相关线性随机效应模型下统计分析和推断的一种代数方法，并得出单个模型与联合模型下相同未知参数函数BLUP间等价的条件。作为单个模型与两个联合模型间关系的延伸，文献 [ 6 ] 考虑在一组带有随机系数和相关协方差矩阵的线性模型中，建立所有未知参数函数BLUP/BLUE解析表达式，并利用一些矩阵分析工具得到单个模型与和组合模型下相同未知参数函数的BLUP/BLUE等价的充要条件。此外文献 [ 7 ] 讨论了两个具有竞争关系的线性随机效应模型下未知参数函数BLUP间的关系等。虽然在线性随机效应模型中分析两个或多个模型间关系这类问题已经做了很多研究，但仍有一些问题有待研究，如对一般线性随机效应模型添加新的变量项成为过参数线性随机效应模型，这时两模型相同未知参数函数预测量间存在什么样的关系？因此本文考虑了一般线性随机效应模型与过参数线性随机效应模型下未知参数函数预测量间的关系。