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

user: 5. 分离膜在印染废水处理中的应用 国家颁布实施《纺织染整工业水污染物排放标准》(GB4287-2012)，对纺织印染废水治理提出了更高的要求，提高水的重复利用率是今后发展的必然趋势，也是纺织企业发展的重要任务 [ 12 ] 。印染废水经处理后，水质需要达到一定的要求之后才能回用，具体要求见表2。 分离膜因分离效率高、操作简单、过程易控制、无二次污染等优点，符合环境保护的可持续发展战略，在印染废水的处理中得到了广泛应用。目前比较成熟的处理印染废水的分离膜技术有微滤膜、超滤膜、纳滤膜及反渗透膜等。 5.1. 微滤膜 微滤膜因其膜孔相对较大、操作压力较低，只能用于简单的粗过滤，去除杂质和细菌病毒等微生物，对于可溶性的有机物去除率较低，通常作为废水处理的预处理阶段。 5.2. 超滤膜 超滤膜可以截留住污水中溶解性的大分子物质，使小分子溶质通过。超滤膜在印染废水处理中可以有效地去除悬浮物、胶体以及降低废水的色度和浊度，出水较为稳定。常见的超滤膜是作为反渗透的预处理阶段。 5.3. 反渗透膜 反渗透膜可以截留微生物、COD Cr 、BOD 5 、小分子以及离子等，只让水和溶剂透过，对钠离子的截留率高达99.5%。在印染废水深度处理回用中，反渗透膜可以对染料进行回收利用，形成了以反渗透为主的一系列废水处理工艺，例如，微滤–反渗透、超滤–反渗透等，其中以超滤–反渗透最为常见。 5.4. 纳滤膜 纳滤膜的孔径介于超滤膜和反渗透膜之间，操作压力也介于超滤和反渗透之间，与反渗透比，操作压力低，能耗也有所降低，而且纳滤的截留率又优于超滤，因此，纳滤在纺织废水处理回用等方面受到人们的重视。纳滤膜的过滤本质除了筛分作用还有静电作用，纳滤膜表面有荷电基团，可以与杂质离子产生静电作用，杂质离子进而吸附成团，被过滤掉 [ 13 ] 。 Table 2 Water quality indicators and limit 序号 项目 限值 1 pH值 6.5~8.5 2 化学需氧量(COD Cr ) (mg/L) ≤50 3 悬浮物(mg/L) ≤30 3 透明度(cm) ≥30 5 色度(稀释倍数) ≤25 6 铁(mg/L) ≤0.3 7 锰(mg/L) ≤0.2 8 总硬度(CaCO 3 计) (mg/L) ≤450 9 电导率(μs/cm) ≤2500 表2. 回用水水质指标及限值 5.5. 膜生物反应器 膜生物反应器(MBR)就是把分离膜与生物反应器相结合的废水处理技术，使膜的性能更加完善。膜生物反应器的特点为占地面积较小，维护比较简单，没有繁琐的反冲洗等环节 [ 14 ] 。根据机理不同可以分为好氧膜生物反应器和厌氧膜生物反应器。 好氧膜生物反应器可以有效的降低纺织废水中的COD Cr 以及色度等参数，通常使用的膜数量极少，运行成本较低。因此，好氧膜生物反应器适合作为纺织废水的预处理，用于降低废水中的COD Cr 及色度，减少污染物在后续膜处理过程中对膜的损害。而带有污泥的废水更适合厌氧膜生物反应器来进行处理，厌氧膜生物反应器能耗低，对矿化有机物的去除更加彻底 [ 15 ] [ 16 ] [ 17 ] 。同时，厌氧膜生物反应器还可以产生沼气，沼气可以燃烧提供能量，达到能量回收的作用。
system: [SYSTEM]

user: 6. 组合工艺处理法 微滤膜、超滤膜、纳滤膜以及反渗透膜处理废水的效果是逐级递增的，废水处理效果越来越明显，相对地，要求也越来越高，对于预处理技术和后处理技术的要求都十分严格。目前，对于单一的膜法水处理技术很难满足当前的产水水质要求以及回用水质要求，所以对于当下膜法水处理技术，更加适宜的是使用多种膜法水处理技术相结合，实现对纺织废水的深度处理回用 [ 18 ] [ 19 ] 。常见的组合工艺有超滤–反渗透、超滤–纳滤、微滤–反渗透等。 6.1. 超滤–反渗透 超滤去除微小杂质效果较好，在反渗透的预处理和中水回用等领域发挥了重要作用，反渗透是当下最先进最有效的膜分离技术，能够去除溶解在水中的离子、微生物、有机物等，去除率均较高。超滤–反渗透集成工艺是当下处理废水使用最多的一种技术。超滤作为反渗透的预处理阶段，能够有效地避免反渗透膜污染，高效稳定的去除污染物 [ 20 ] [ 21 ] 。常见的超滤–反渗透工艺流程如下： 6.1.1. 水解酸化–A/O–超滤–反渗透–回用 对于印染废水中色度高的问题采用水解酸化–A/O作为超滤–反渗透双膜系统的预处理阶段，可以有效地降低色度，出水水质满足后处理系统的进水水质。例如苏南某印染公司以开发销售高档棉类针织印花布，企业生产过程中主要使用的是活性染料，助剂种类繁多、复杂。使用上述工艺对印染废水进行深度处理(图1为工艺流程)，运行结果表明：COD Cr 和色度去除率分别为94.2%和98.7%，系统出水电导率55.3~83 μs/cm，硬度 ≤ 10 mg/L，出水未检出铁和锰，出水水质达到了《纺织染整工业水污染物排放标准》(GB4287-2012)的要求。水的回用率在50%以上，回用水可直接用于企业印花、水洗阶段 [ 22 ] [ 23 ] 。 但水解酸化–A/O系统只能适用于低浓度的印染废水，对于高浓度的印染废水的去除效率较低，不能满足超滤–反渗透的进水水质要求。 图1. 工艺流程图 6.1.2. 臭氧氧化–曝气生物滤池–超滤–反渗透–回用 针对纺织废水中难生物降解有机物的去除和脱盐等问题，在印染废水的处理过程中采用臭氧氧化–曝气生物滤池的预处理系统，系统中的臭氧–曝气生物滤池具有臭氧氧化和生物处理的双重作用，可以有效地提高难生物降解有机物的去除率和脱盐率，该预处理已经广泛应用到了废水处理等方面 [ 21 ] [ 24 ] [ 25 ] 。如广东溢达纺织有限公司选用该工艺流程作为纺织废水的深度处理工艺，经过预处理系统，出水COD Cr 、浊度、SS、氨氮、色度的去除率分别为72.7%、45.3%、82.6%、78.8%、89.3%，研究表明，该预处理系统有效保证了后续膜处理的进水水质，保证了膜系统的稳定运行。再经过超滤–反渗透膜系统的处理，膜系统产水pH值7.4~7.9，电导率50~200 μs/cm，总硬度2~10 mg/L，总碱度25~60 mg/L，硫酸盐22.4~90.5 mg/L，膜系统产水达到回用标准要求 [ 10 ] 。 预处理系统中的臭氧–曝气生物滤池一体化装置与先臭氧氧化、再进曝气生物滤池的分开系统相比具有处理效率高，占地面积小，能耗低，易于操作等优势。且经过该预处理系统的出水可以回用于生活杂质水、景观水体以及要求不高的印染工序 [ 26 ] 。 6.1.3. 多介质过滤器–超滤–反渗透–回用 超滤/反渗透膜系统另一常见的预处理系统为多介质过滤器，常见的为石英砂过滤器和活性炭过滤器。其中石英砂过滤器过滤掉进水中的细小颗粒、悬浮物及胶体等杂质，使出水污染指数(SDI) ≤ 5；活性碳过滤器的主要作用是有效去除水中残余的游离氯，以及吸附从石英砂过滤器泄漏过来的小分子有机物等污染性物质，进一步降低RO进水的SDI [ 27 ] [ 28 ] 。 如广东某印染公司主要染色加工真丝、仿丝、全棉等质地坯布，采用石英砂过滤器–活性炭过滤器–超滤–反渗透的工艺处理印染废水(图2为该公司的工艺流程图)，运行表明：经超滤/反渗透膜系统处理后，出水各项指标均优于纺织染整工业废水治理工程技术规范中的染色回用水水质标准，其中，出水色度为2倍，电导率为17 μs/cm，COD Cr 和BOD 5 均未检出，出水水质满足染色工序的回用 [ 28 ] 。 针对超滤反渗透体系的预处理阶段的出水水质进行了对比(见表3)，多介质过滤器对于废水的处理只是简单的吸附过滤，出水水质相对较差，主要针对低浓度印染废水；水解酸化–A/O体系与臭氧–曝气生物滤池体系的废水处理的出水水质相当，且两者对色度和COD Cr 的去除效果较为显著，但是臭氧–曝气生物滤池的运行成本较低，占地面积较小，不会产生污泥膨胀，氧传输效率高。 6.2. 膜生物反应器–纳滤 用膜生物反应器(MBR)来替代传统的生物处理–臭氧氧化–活性炭吸附工艺作为纳滤的预处理系统，得到的出水不但水质高而且水质稳定，虽然出水不能直接达到回用的标准，但是与传统的工艺相比，COD Cr 和色度去除率显著提高。再通过纳滤可以去除残余的COD Cr 、色度以及部分脱盐，出水水质满足回用标准，回用率为54.4% [ 29 ] 。 6.3. 微滤–纳滤 由于纳滤膜对进水浊度有着严格的要求，微滤作为纳滤的预处理阶段，以去除废水中的絮状沉淀和杂质降低废水浊度，有研究表明，微滤膜对浊度的去除效果较为显著，出水水质浊度 ≤ 0.1NTU [ 14 ] [ 30 ] 。 连续微滤技术 [ 31 ] (CMF)是集微絮凝技术、现代膜分离技术和控制技术于一身，即原水通过加药混合、混凝反应、微滤膜的微孔分离，从而使原水达到凝聚、过滤净化处理的目的，而设备的运行、气水双洗、反洗和排污等过程均由PLC自动控制完成。连续微滤技术相比于微滤具有出水量连续稳定的特点，且操作简单、维修方便、系统的出水满足纺织印染回用水质的要求。 例如南方某纺织工业园用上述工艺处理废水，工艺流程如图3，连续微滤技术采用的是天津膜天膜 图2. 超滤/反渗透膜技术工艺流程 图3. 工艺流程图 Table 3 Comparison of water quality in the pretreatment stag 工艺参数 色度去除率(%) COD cr 去除率(%) 出水电导率μs/cm 硬度mg/L 浊度NTU 运行成本(元/t) 水解酸化–A/O 88.5 75.7 55.3~83 ≤10 - 3.2 臭氧–曝气生物滤池 89.3 72.7 50~200 2~10 78.8 2.9 多介质过滤器 62 56.8 1346 - - 2.75 表3. 预处理阶段出水水质对比 工程技术有限公司生产的三支膜中试试验装置，产水量为6 t/h。采用外压式中空纤维膜元件和独特的气水双洗工艺技术，可以除去废水中的微生物、胶体和悬浮物等，出水浊度 ≤ 0.1 NTU。纳滤系统采用的是由海德能公司生产的单支膜系统，型号为ESNA1-4040。该膜能在超低运行压力下工作，从而达到节能、降低设备投资和操作费用的目的 [ 32 ] 。因此，该工艺可以在工业园区印染废水处理中进行推广。
system: [SYSTEM]

user: 7. 总结 纺织废水因其水量大、成分复杂、色度高、碱性大、有机污染物浓度高，处理起来较为困难。目前单一的处理工艺已经不能满足出水水质要求，将不同处理工艺进行优化组合，可以得到更好的处理效果。近年来，分离膜因其独特优异的技术特点在纺织废水处理工艺中得到广泛的应用，传统废水处理技术与分离膜技术的组合工艺成为人们研究的焦点，同时对组合工艺进行不断的优化改进，使工艺技术更加成熟稳定，这样既提高了处理效率，又降低了处理成本。
system: [SYSTEM]

user: 摘 要 本实验对4只患有子宫积脓的犬和2只正常的母犬进行了子宫切除，并在扫描电镜、透射电镜下观察两者黏膜层发生的变化。通过观察所发生的病变，了解发病机理，为该病的治疗提供相应的指导，从而提高该病的治愈率，也为国内在这方面所进行的微观结构的观察提供更多的资料。扫描电镜20 KV下50倍到300倍其上皮脱落呈现蜂窝状结构，下层组织暴露，管腔扩张，组织间隙增宽，呈现水肿状态。透射电镜3500倍到8000倍下观察细胞内膜破裂细胞基质及细胞器外渗，并且可以看到部分的细胞器，例如线粒体和内质网都会增大，线粒体空泡化，分泌的分泌物会增多。 关键词 :子宫积脓，犬，扫描电镜，透射电镜 Copyright © 2018 by authors 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 ] 。修晓茜认为，子宫内膜囊性增生，腺体增多是发生子宫蓄脓的重要过程 [ 3 ] 。钱存忠等发现，超显微结构下子宫内膜上皮细胞形态、结构、绒毛长短及细胞完整性有明显的变化 [ 4 ] 。张志强发现，发生子宫积脓的犬，子宫的组织学损伤表现为内膜的增厚或萎缩，以及子宫腺不同程度的囊性变化 [ 5 ] 。子宫积脓多发于发情期之后，引起其发病因素较多，例如不洁交配，细菌感染、生殖激素分泌、年龄等，患病的犬如果不进行治疗，常会因脓毒败血症而死亡，本课题组近年来一直在从事相关研究，本次实验是从超显微结构对比的角度，观察正常与病理情况下子宫黏膜层所发生的变化。
system: [SYSTEM]

user: 2. 样品采集与处理 2.1. 样品的采集 病理子宫样品来自于动物医院门诊接诊的病例，通过手术方法进行采集，包括卵巢、子宫角和部分子宫颈，先用生理盐水进行清洗，将污染物和分泌物清洗干净后放入戊二醛固定液中，进行固定，编号标记；正常样品为门诊接诊的做绝育手术的母犬子宫颈和卵巢，处理方法同病理样品。 2.2. 扫描电镜的样品制作 将清洗干净的子宫切成0.8 cm × 0.8 cm的样品，放入2%戊二醛固定液中固定24 h，固定好的样品用0.1 mol/L磷酸缓冲溶液冲洗3次，每次15 min，后将样品置于不同浓度的乙醇逐级脱水(30% 15 min、50% 15 min、70%过夜、80% 15 min、90% 15 min、95% 15 min、100% 10 min 3次)，之后无水乙醇/丙酮(1:1)、无水乙醇/丙酮(1:4)、纯丙酮进行置换，后喷镀金膜，扫描电镜观察。 2.3. 透射电镜的样品制作 固定：将清洗干净的子宫切成大约1 mm × 1 mm × 3 mm长块状，加入1 ml 2.5%的戊二醛固定液固定24 h。细胞固定后，弃固定液加入1 ml 0.1 mol/L磷酸盐缓冲液(pH7.2)连续冲洗5次，每次15 min。再用1%锇酸固定组织1.5 h，之后再用磷酸盐缓冲液(0.1 mol/L)冲洗5次，每次15 min。脱水：依次用乙醇逐级脱水30% 15 min、50% 15 min、70% 15 min、80% 15 min、90% 15 min、95% 15 min、100% 10 min。置换： 无水乙醇/丙酮(1:1)、无水乙醇/丙酮(1:4)、纯丙酮进行置换，EPon812包埋剂/丙酮(1:1)、EPon812包埋剂/丙酮(4:1)、纯EPon812包埋剂对细胞进行渗透，EPon812包埋剂对细胞进行包埋，分别在37℃、42℃、60℃温度下聚合72 h，100 nm厚度切片，并用铜网捞取切片，进行醋酸铀和柠檬酸铅双重染色，电镜观察 [ 6 ] 。
system: [SYSTEM]

user: 3. 观察结果 3.1. 扫描电镜 正常的子宫表面较为光滑腺体开口排列有序，开口也较为的明显，并且能看到表面的褶皱，正常的子宫纤毛和微绒毛及组织结构排列有序，管腔结构很明显(图1，图2)，发生病变后上皮脱落呈现蜂窝状 图1. 正常子宫内膜表面 图2. 正常子宫内膜表面 结构，下层组织暴露，管腔扩张，组织间隙增宽，呈现水肿状态(图3，图4)。表明：发生子宫积脓后首先黏膜脱落较为，后影响黏膜下层组织结构。 3.2. 透射电镜 正常的子宫细胞排列较为整齐(图5，图6)，发生子宫积脓后细胞内膜破裂细胞基质及细胞器外渗(图7)，并且可以看到部分的细胞器，例如线粒体和内质网都会增大(图8，图9)，分泌的分泌物会增多。我们在图10中可以看到链球菌，并且在链球菌周围看到大量的分泌颗粒，说明在细菌感染下，会引起分泌细胞进行分泌。表明：蓄脓后对子宫内黏膜表面细胞发生很大的影响，使其结构发生改变，细胞器大小改变，细胞膜发生改变。 图3. 异常子宫黏膜表面 图4. 异常子宫组织间隙 图5. 正常子宫内膜表面细胞 图6. 正常子宫细胞器 图7. 异常子宫胞质内细胞器 图8. 异常子宫线粒体 图9. 异常子宫细胞间质 图10. 异常子宫细菌感染
system: [SYSTEM]

user: 4. 讨论 在扫描电镜下张志强认为发生子宫积脓后整个内膜表面凹凸不平，内膜细胞大小形状不均匀，囊性变化明显，有的顶端出现破损，内部形成囊腔 [ 5 ] 奶牛产后子宫内膜只在某些区域脱落，并非整个内膜都脱落。在那些上皮细胞相对完整的部位，上皮细胞表面看不见纤毛，微绒毛很稀少并且不整齐。在那些内膜破损的区域，破损部分表面凸凹不平，有纤维状结构，有的地方还可见到内膜破损的断面。在产后13 d，破损区域缩小，因为在产后的这个时期，破损部位有新生细胞出现，说明子宫内膜开始了修复的过程 [ 7 ] 。透射电镜下部分病例的炎症与细菌感染有关。子宫内膜增生，上皮细胞的分泌增强，大量分泌物在子宫腔内集聚，炎性细胞浸润加之细菌的感染导致该病的发生与发展，几乎所有的子宫蓄脓病例均出现炎症 [ 4 ] 。 在扫描电镜下我们可以观察其正常生理结构表面是光滑，并且能看到排列整齐的分泌腺口发现发生子宫积脓的犬的子宫内膜表面黏膜脱落，零散排列，呈现蜂窝组织，组织间隙增宽呈现水肿状态。因此在治疗该病我们首先应该治疗其表面黏膜，制止黏膜的脱落，及修复黏膜，下层组织在发生水肿的状态，会影响其渗透压的变化， 透射电镜下正常的子宫细胞排列较为整齐，发生子宫积脓后细胞内膜破裂细胞基质及细胞器外渗，并且可以看到部分的细胞器，例如线粒体和内质网都会增大，线粒体空泡化，细胞器分泌的分泌物会增多。我们在图10中可以看到链球菌，并且在链球菌周围看到大量的分泌颗粒，说明在细菌感染下，会引起分泌细胞进行分泌。由此可以看出炎症反应后细胞内分泌会增加，治疗该病我们应阻止其分泌细胞的继续分泌，线粒体快速增大，会发生空泡化，治疗时应从防止空泡化的继续发生的角度来进行。此实验对我们保守治疗开放型子宫蓄脓，提供治疗方向，从微观的方向来继续我们的治疗，以达到更好的治疗效果。
system: [SYSTEM]

user: 5. 结论 犬子宫积脓的发生由多种因素引起，正常子宫内膜光滑，腺体开口明显，排列整齐；发生子宫积脓后子宫内膜脱落，子宫下层组织发生水肿，呈现蜂窝组织结构，组织间隙增宽，管腔腺体结构萎缩。细胞结构发生的变化：正常细胞器，细胞核周围细胞器分布均匀，细胞外膜完整；病变后细胞器增大，出现空泡化，部分细胞发生破裂，细胞器外渗，内分泌物增多。
system: [SYSTEM]

user: 摘 要 为使高分辨雷达图像更精确反映舰船目标的结构、形状特点，研究了基于后向散射场数据的舰船目标高分辨雷达成像技术。建立了舰船目标三维模型，采用三角面元对模型进行剖分，利用物理光学法计算面元的镜面反射，应用等效电流法计算模型棱边部分绕射，依据机载雷达聚束式合成孔径(SAR)成像原理，得到舰船目标在不同姿态角和频率采样点的后向散射场。通过数据变换得到目标的一维距离像以及SAR图像，在此基础上研究了入射频率、方位角和极化方式对舰船目标探测的影响。仿真结果表明，该技术应用后向散射场数据得到的雷达SAR图像能反应出舰船的基本特征，可应用于舰船目标自动识别。 关键词 :舰船目标，后向散射场，雷达成像，目标识别 Copyright © 2018 by authors 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 ] - [ 10 ] 。国内外学者在舰船目标的高分辨雷达成像和自动识别技术方面做了大量工作，其中基于高分辨一维距离像与合成孔径雷达(SAR)图像的雷达目标特征提取及识别方法是研究的热点 [ 11 ] [ 12 ] [ 13 ] [ 14 ] [ 15 ] 。 当目标雷达散射处于光学散射区，雷达发射信号带宽很大时，复杂目标占据多个雷达距离分辨单元，其雷达回波称为一维距离像，若雷达在方位角旋转一定角度，就可以得到目标合成孔径雷达SAR成像。与普通雷达相比，宽带雷达探测精度更高，能够精细反应目标的尺寸结构、形状和散射点分布等信息，更有利于雷达进行自动目标识别。 以往大多数研究高分辨雷达成像技术通常采用强散射点建模，而在实际应用中，当雷达照射目标角度变化时，目标的强散射点会发生变化，因此采用散射点模拟目标散射会存在一定的误差，不能够真实反映雷达工作情况。本文采用三角面元构建复杂目标三维模型，利用高频近似算法计算目标的后向散射场数据，通过与FEKO软件计算结果进行比较，验证了算法的正确性，研究舰船目标的散射特性，并在此基础上得到舰船目标的SAR图像，以及不同角度雷达照射情况下目标的一维距离像。分析了频率、姿态角、极化方式等各项参数对高分辨雷达成像的影响。
system: [SYSTEM]

user: 2. 机载SAR场景及成像参数 机载聚束式SAR成像场景如图1所示。飞机飞行高度为H，距离目标的距离为R，飞行速度为v，其视线相对于目标坐标系俯仰角为θ，方位角为φ。 进行电磁散射建模，假设入射方位角和俯仰角均匀采样，采样点为 { ( θ m ， φ m ) ， m = 1 ， 2 ， ⋯ ， M } ，据预先确定的中心视角 ( θ c ， φ c ) 、视角变化范围 ( Δ θ ， Δ φ ) 。照射频率变化采用频率步进的方式均匀采样，采样点为 { f n ， n = 1 ， 2 ， ⋯ ， N } ，计算舰船目标的后向散射场数据用以模拟实际探测雷达回波，在此基础上利用适当的算法对其进行成像。根据设定的中心频率f c 和带宽B，以及的角度采样点M和频率采样点N，可以得到雷达成像场景条件参 [ 3 ] ： 图1. 机载SAR成像场景 θ m = ( θ c − Δ θ 2 ) + m ⋅ Δ θ M − 1 ，     m = 0 ， 1 ， ⋯ ， M − 1 (1) φ m = ( φ c − Δ φ 2 ) + m ⋅ Δ φ N − 1 ，     m = 0 ， 1 ， ⋯ ， M − 1 (2) f n = ( f c − B 2 ) + n ⋅ B N − 1 ，     n = 0 ， 1 ， ⋯ ， N − 1 (3) 假设成像时间段内雷达视线转过了δ角度，并且在成像中心时刻，雷达处于正视状态，其方位角为φ c ，俯仰角为θ c ，与目标中心相距R 0 ，可以得到： θ c = arcsin ( H R 0 ) ， φ c = 0 (4) 入射角度的采样数M即为整个观测期间探测雷达发射的脉冲个数，它由雷达与目标之间距离R、入射角度变化范围δ、雷达脉冲重复时间(Pulse Repetition Time， PRT)和载机飞行速度v决定： M = R δ v ⋅ P R F (5) 依据实际情况下探测目标最小不模糊距离来确定频率采样点个数N。
system: [SYSTEM]

user: 3. 舰船目标电磁散射计算 考虑工作波长为厘米级的机载雷达对舰船目标的探测性能。在该情况下，舰船目标为电大尺寸目标，此时目标散射处于光学区域，各散射单元之间的耦合作用很小，通常采用高频近似算法计算目标的电磁散射特性。 3.1. 目标剖分模型 首先建立舰船目标的三维模型，将其导入Hypermesh软件，依据平面波入射频率以及计算精度进行网格剖分，图2为利用三角面元对某型舰船目标进行剖分后的三维模型，该舰船长100 m，宽13 m，桅杆高度20 m。 3.2. PO + MEC算法计算目标散射 物理光学(PO)法的基本原理，当目标表面的曲率半径远远大于照射电磁波的波长时，假设在被入射 图2. 舰船剖分模型 波直接照明的区域存在表面电磁流。可以使用目标表面的等效感应电磁流的散射来代替散射体自身的散射，此时目标的散射场通过解析公式进行求解。考虑平面波入射，物体表面远区散射场可由根据Stratton-Chu方程 [ 16 ] 进行计算： E s p o ( r ) = j k 4 π e − j k r r ∫ S s ^ × ( M s ( r ′ ) + η 0 s ^ × J s ( r ′ ) ) ⋅ exp [ j k r ′ ⋅ ( s ^ − i ^ ) ] d s ′ (6) 式中， i ^ 是入射方向单位矢量， s ^ 是散射方向的单位矢量， η 0 表示自由空间的波阻抗， S 表示目标表面被照射波照射部分。 J s ( r ′ ) 和 M s ( r ′ ) 表示 S 上的电磁流矢量幅度，其中的相位因子已包含在指数项中。 J s ( r ′ ) 和 M s ( r ′ ) 在两个基本物理光学假定的条件下，根据等效原理， M s ( r ′ ) = E ( r ′ ) × n ^ ， J s ( r ′ ) = n ^ × H ( r ′ ) (7) 式中， E ( r ′ ) 是目标表面上总电场， H ( r ′ ) 是目标表面总磁场。此时，目标表面的总场是入射场和反射场的总和，即 E ( r ′ ) = E i ( r ′ ) + E r ( r ′ ) (8) H ( r ′ ) = H i ( r ′ ) + H r ( r ′ ) (9) 式中， E i ( r ′ ) ， H i ( r ′ ) 和 E r ( r ′ ) ， H r ( r ′ ) 分别为入射场反射场矢量。 通过考虑目标棱边的绕射对PO的结果进行改进。采用文献 [ 17 ] 中的等效电流法(MEC)计算目标棱边绕射，该方法应用等效电磁流的散射等效边缘绕射场。通过对等效电磁流进行积分来计算棱边在任意方向的散射场。对于边缘为 C 的劈边，其远区边缘绕射场可通过下式进行计算： E s d = j k 4 π e − j k r r ∫ C { η 0 s ^ × [ s ^ × J ( r ′ ) ] + s ^ × M ( r ′ ) } e j k s ^ ⋅ r ′ d l (10) 式中， J ( r ′ ) = I e ( r ′ ) t ^ 表示边缘等效电流， M ( r ′ ) = I m ( r ′ ) t ^ 表示边缘等效磁流， t ^ 表示劈边 C 的切向单位矢量。I e 与I m 的定义见文献 [ 17 ] 。 目标总的散射场为PO法计算的镜面散射以及MEC法计算的棱边绕射之和，即 E t o t a l = E s p o + E s d (11) 用宽带信号照射舰船目标，利用PO + MEC算法计算舰船目标在不同频段、不同姿态角下的后向散射场数据，流程如图3所示。 例如，对某一舰船目标在n个姿态角进行照射，得到的散射场数据为： Y = [ Y ( 1 ) ， Y ( 2 ) ， ⋯ ， Y ( n ) ] (12) 图3. PO + MEC算法流程图 每一列数据对应一个姿态角下不同频率采样，每一行数据对应一个频率下不同姿态角采样。对数据做二维傅里叶变换(FFT)可以得到目标的SAR成像，对每一列直接做逆傅里叶变换(IFFT)即可得到目标高分辨距离像。
system: [SYSTEM]

user: 4. 仿真计算 采用第二节中的三维舰船模型，舰船放置于xoy平面，中心位于坐标原点，舰首对应x轴向正方向，舰船轴线垂直于y轴(见图4)，剖分精度为 λ / 6 ，总共18225个顶点，36，446个三角面元，54，701个公用棱边。 设照射频率为2 GHz，应用纯PO算法以及PO + MEC算法计算边长为0.5 m的立方体后向散射，得到水平极化(HH)雷达散射截面积(RCS)如图5所示，图中给出了应用FEKO软件计算的结果。计算结果表明，两种方法的计算结果吻合的较好，验证了PO + MEC算法的正确性。由图可知，通过加入计算棱边绕射，MEC很好地改善了PO的计算结果。 图6所示为应用本文算法计算舰船模型的后向散射，方位角φ = 0˚，俯仰角θ变化范围0~180˚，雷达照射频率2 GHz，图6(a)为HH极化结果，图6(b)为VV极化结果。计算数据表明，舰船散射在0˚ (舰首)、90˚(舰船正上方)以及180˚(舰尾)方向最大，VV极化情况下RCS略小于HH极化结果，尤其是30˚~70˚范围内HH极化结果明显较大。 取俯仰角θ分别为30˚和60˚，方位角φ变化范围0˚~360˚，相当于雷达绕舰船一周，其余参数不变，得到的散射结果如图7所示，其中图7(a)为HH极化结果，图7(b)为VV极化结果。由图可知，舰船RCS在舰首(0˚)、舰尾(180˚)、船舷(90˚和270˚)等方向最强，且随着俯仰角增大而增大，同图5中结果类似，HH极化散射强于VV极化结果，这是由于HH极化二向性散射强于VV极化，因此HH极化更能反映舰船的上层建筑结构，导致HH极化散射强于HH极化。因此，应用HH极化雷达探测舰船目标更具有优势。 图4. 电磁散射场景 图5. 算法验证 图6. 不同俯仰角舰船目标RCS；(a) HH极化；(b) VV极化 考虑典型L波段机载搜索雷达，照射频率为f c = 2.0 GHz，带宽为B = 150 MHz，频率范围：1.925 GHz~2.075 GHz，采样间隔Δf = 1 MHz，共N = 150个采样点，此时距离分辨率为ρ r = C/2B = 1.0 m，假设雷达载机飞行高度距离地面高度9 km，目标距离200 Km，此时俯仰角θ = 2.6˚。在HH极化情况下，计算后向散射场数据，进行IFFT变换并对幅值归一化的到舰船目标的一维距离像如图8所示，图8(a)为方位角φ = 0˚的结果，此时雷达从舰首入射，有图可以看出，舰船在雷达视线的径向长度为100 m左 图7. 不同方位角舰船目标RCS；(a) HH极化；(b) VV极化 图8. 不同方位角舰船目标一维距离像 右，主要反射部位集中于船身20~90 m区域，这正是舰船上层建筑最为集中的区域。图8(b)为方位角φ = 90˚的结果，相当于雷达从船舷部照射，此时雷达径向距离径向长度为14 m左右，对应的是舰船目标的船体宽度。因此，通过不同姿态角舰船高分辨一维距离像可以得到舰船长度、舰身宽度、舰桥上层建筑尺寸和数量等船体结构信息，这些信息均可应用于后续进行雷达自动目标。 其余参数不变，HH极化情况下，取方位角变化范围10˚，采样间隔为0.05˚，采样数M = 201个，得到150行和 201列大小的矩阵数据，对其进行FFT变换，即可得到舰船目标SAR成像，见图9，其中图9(a)为直接FFT变换结果，图9(b)为加窗FFT变换，窗函数为布莱克曼窗，经过对比可知窗函数能够很好的改善成像结果。从图像中可以分辨出舰船的尺寸与图7中相同，同时可以发现舰船强散射点集中于甲板上建筑反射。由于不同舰船目标的长度、宽度和长宽比都是独特的，因此经过提取这些特征参数，可以很好的应用于雷达自动目标识别。图9(c)为舰船旋转45˚的SAR图像，与图9(b)相比，此时的舰船目标强散射点增加，这是由于船舷部分的散射较强造成的。同时该结果证明了本文算法对不同姿态目标能得到很好的结果，具有很好的适应性。 图9. 舰船目标SAR图像；(a) 不加窗图像；(b) 加窗图像；(c) 旋转45˚成像
system: [SYSTEM]

user: 5. 结论 本文在建立舰船目标三维模型的基础上，研究了舰船目标L波段探测雷达照射下不同姿态角和极化方式下的电磁散射特性。结果表明，舰船目标的散射主要集中于舰首、舰尾及船舷侧，并且HH极化反射要强于VV极化。根据SAR成像原理，计算了不同方位角和频率采样点的后向散射数据，通过数据变换得到目标的一维距离像和SAR图像。计算结果表明，舰船目标高分辨雷达图像能够较好的反应舰船的尺寸、形状和强散射点数量，这些特征在雷达自动目标识别中具有重要应用，同时该结果也证明了文中三维模型及高频近似算法的准确性。研究结果对提高雷达舰船目标探测与识别具有一定的参考意义。
system: [SYSTEM]

user: 摘 要 本工作中利用高脉冲Nd:YAG脉冲激光器(波长532 nm，脉宽：4 ns)产生Ni等离子体，研究了Ar、He、N 2 气体对激光诱导Ni等离子体光谱的影响。通过实验研究发现氩气对Ni等离子体有明显增益效果。计算三种环境气体中产生Ni等离子体电子温度随气压变化关系，发现气压对Ni等离子体具有束缚作用。 关键词 :光谱分析，缓冲气体，Ni等离子体 Copyright © 2018 by authors 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 ] 。在激光诱导等离子体特性研究中等离子体电子温度是诊断等离子体特性和判断电离度的物理量 [ 3 ] 。等离子体作为物质的不稳定状态，易受到外界环境影响，不同环境条件下产生的等离子体特性不一样，因此研究不同环境条件对激光诱导等离子体产生影响一直备受关注 [ 4 ] 。郑贤锋等人 [ 5 ] 利用CO 2 分子激光器在铝靶面产生等离子体以氩气，氦气，空气为环境气体研究所产生Al等离子体特性，发现压强10 −4 Pa~2 × 10 3 Pa改变过程中电子密度从1.95 × 10 16 cm −3 增加10.5 × 10 16 cm −3 。吴涛等人 [ 6 ] 利用氩气，氦气，空气为环境气体对Al等离子体光谱进行研究发现Al等离子体光谱在氩气中强度最高，氮气中次之，空气中产生的等离子体光谱背景信号干扰较大。沈跃良等人对环境气体为氩气，氮气，空气，CO 2 气对飞灰中C，Si，Fe，Al和Mg谱线特性进行了研究发现氩气中谱线强度最高 [ 7 ] 。虽然很多研究中发现氩气对谱线强度有较大的增益效果但是解释其中原因时存在不同看法。 为了更好的解释氩气对等离子体谱线增益效果原因本文利用高脉冲Nd:YAG激光器(波长532 nm，10 Hz)为激光源，金属Ni为作为样品，氩气，氦气，氮气作为缓冲气体，观测和讨论了三种缓冲气体对激光诱导Ni等离子体光谱特性影响并通过计算等离子体温度随气压变化，讨论三种气体气压对Ni等离子体影响。最后通过碰撞来讨论和解释三种气体对等离子体影响原因。
system: [SYSTEM]

user: 2. 实验装置和测量方法 图1为实验设备装置图，其简述如下：试验中把样品放在可充缓冲气体的不锈钢腔体中，利用Nd:YAG脉冲激光器作为激光源，以7 cm的透镜聚焦在腔体里面的样品表面，45˚方向利用4.5 cm的透镜从样品表面散射回来的光聚焦在光纤头上，再利用光谱仪(Andor Solis)对光进行分光，通过增强型器件ICCD信号增强放大后其数据发送电脑软件生成Ni等离子体光谱图像。实验中激光能量80 mJ，气压范围2000 Pa~15，000 Pa，测量波长范围324~374 nm，光谱仪光栅条数为1200/mm，门宽为100 ns，在不同环境气体中采集等离子体时间分辨光谱谱线。 2.1. 三种缓冲气体对谱线的增益讨论 时间分辨光谱图是反映处于激发态原子随延迟时间变化关系。一般讨论时间分辨光谱图是选择特征谱线后，再绘出此条谱线强度随时间演化关系。本文为了比较三种不同缓冲气体对Ni等离子体的影响，利用光谱强度–波长–时间的三维时间分辨光谱图，反映出324~374 nm范围所有谱线随时间演化关系。 图2是Ni等离子体在不同缓冲气体中光谱强度随时间演化图，在氩气中产生等离子体光谱最强，在 图1. 实验设备装置图 图2. Ni等离子体在不同缓冲气体中光谱强度随时间演化图 氮气中产生等离子体光谱次之，在氦气中产生的等离子体光谱强度最低，从图中看出随着延迟时间1~10 μs变化过程中谱线强度减弱，其原因是处于高激发态的原子、电子和离子数目减少引起。氩气中衰减时间比氮气和氦气约长1 μs。 2.2. 气压对Ni等离子体电子温度影响 电子温度是诊断等离子体的一种物理参数，一般等离子体电子温度计算利用多谱线Boltzmann拟合方法 [ 8 ] [ 9 ] ，其表达形式为： ln ( I m λ m A m g m ) = E m K B T + C (1) 其中I m 是谱线强度，λ m 是谱线波长，A m 是谱线跃迁几率，g m 统计权重，E m 是激发能量，K B 是波尔兹曼常数，T e 是电子温度，C是常数。本文在329 nm~374 nm之间所观测的Ni的最强谱线处于349.30 nm，352.45 nm，356.64 nm，361.94 nm。从NIST美国谱线库网站中查到有关谱线参数 [ 10 ] ，其参数值如表1所示。 计算等离子体电子温度时，根据公式(1)绘出图，利用Origin进行线性拟合就可以得到直线斜率，从中可以求出电子温度T e 。 本文中为了更好比较氩气，氮气氦气对不同气压下对Ni等离子体的影响，环境气压范围改变为2000~15，000 Pa，测量了Ni等离子体的电子温度三种气体中随气压变化关系，其结果如下所示。 从图3中可以看出气压从2000 kp到15，000 kpa改变过程中，三种气体环境中产生的Ni等离子体电子温度随着气压增加而升高。比较三种气体同气压中电子温度得知氩气中产生的等离子体电子温度最高，氮气中次之，氦气对等离子体电子温度的增益效果比其他两种气体低。环境气压从2000 kp到15，000 kpa改变过程中氩气中电子温度从6450 K升高7778 k，氮气中电子温度从5225 K升高6970 k，在氦气中从4017 K升高到5524 K。从总体来看2000 kpa~15，000 kpa改变中，同等气压中氩气对Ni等离子体电子 Table 1 Spectral line related parameter λ m A m (×10 8 /S) g m E m (eV) 349.30 0.98 3 3.66 352.45 1.00 7 3.56 356.64 0.56 5 3.90 361.94 0.66 7 3.85 表1. 谱线相关参数 图3. 三种环境气体中 Ni等离子体电子温度随压强变化 温度有明显增益作用。从随着气压温度变化角度看随着气压变化氦气对电子温度的增益明显一些。这说明气压对等离子体的膨胀具有束缚作用。
system: [SYSTEM]

user: 3. 结论 通过研究在三种不同环境气体中产生的Ni等离子体光谱，可以发现氩气对谱线有明显的增益效果。郑贤锋等人比较氦气和氩气对Al等离子体的影响，发现氩气对等离子体光谱增强作用很明显，其原因认为氦气电离(24.59 eV)能量比氩(15.76 eV)气大，氩气比氦气更容易电离产生自由电子，使电子数目增加，碰撞加剧，这样电子温度升高。但本实验中氩气对谱线增益效果很明显，氮气次之，氦气最低。如果按电离能出发，氦气、氮气和氩气电离能分别是24.59 eV、14.53 eV、15.76 eV，氮气比氩气更容易电离。因此在氮气中观察到的谱线强度应该比氩气强，按电离能大小来解释此问题很困难。本文认为可以利用碰撞理论来解释。等离子体碰撞是属于多体碰撞，当激光照射在Ni表面时，在激光作用下Ni发生电离，产生等离子体。所产生的等离子体粒子从样品表面往外热膨胀，在热膨胀过程中等离子体粒子与环境气体原子分子发生弹性和非弹性碰撞，能量以辐射形式发出。比较三种环境气体摩尔质量得知氩气(39.5 g/mol)最大，氮气(28 g/mol)次之，氦气(4 g/mol)最小。从样品表面出来的电子和离子与环境气体发生碰撞时因为氩气的摩尔质量大，碰撞后反回到靶核附近，这个过程反复出现，使靶核附近的电子温度升高。通过观察等离子体温度随压强变化关系发现气压对等离子体具有束缚作用。激光能量不变情况下，随着气压的升高，所产生的Ni电子温度增加。其原因随着气压升高产生的Ni等离子体周围环境气体原子分子增多，随着所包围的气体原子核分子增多，碰撞越激励。等离子体与环境气体原子分子碰撞后，返回到靶核附近的可能性越大，这样使等离子体困住在靶核附近引起电子温度升高。氩气原子摩尔质量比其他两种气体较大把大量等离子体困住在靶核附近。然而氦气因为原子质量很小，在气压低时碰撞返回到靶核附近的粒子很少，所以束缚作用不太明显。但是随着气压增高氦气气体分子增多，这让碰撞反射回到靶核附近粒子数目增大，使电子温度升高，起束缚作用。
system: [SYSTEM]

user: 摘 要 房地产行业属于资金密集型的行业，房地产上市公司的规模大，业务多，固有风险比较高。因此有必要聘请会计师事务所对公司的内部控制进行审计，出具内部控制审计报告评价企业内部控制的有效性。本文以2014~2016年A股房地产87家上市公司数据为依据，利用STATA软件和面板数据模型进行实证分析，研究了A股房地产上市公司内部控制审计质量对高管薪酬的影响。研究结果表明：内部控制审计质量越高，高管薪酬越高。本文为我国房地产行业设计有效的内部控制制度，提高内部控制审计质量，促进内部控制的设计和运行的有效性，研究内部控制审计质量对公司治理结构的影响提供了经验证据。 关键词 :内部控制，审计质量，高管薪酬，面板数据 Copyright © 2018 by authors 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. 引言 现代经济活动是建立在一系列契约关系之上，公司制的两权分离导致委托代理关系的产生，继而产生代理成本、信息不对称和利益冲突等问题。高管作为理性经纪人，出于自利动机追求个人利益最大化而偏离股东目标，于是产生代理冲突等一系列问题。如何实施有效内部控制使两者目标趋于一致，使企业行为符合相关法律法规，保证企业财务报告可靠并最终实现经营活动的效率和效果的制度安排，这就是内部控制制度(COSO， 1992)所应发挥的作用。随着2016年最新内部审计准则的制定和颁布，我国更加重视上市公司内部控制设计和运行的有效性，以防范上市公司内部控制的重大缺陷从而导致财务报表的重大错报。早在2008年财政部会同证监会、审计署、银监会、保监会制定了《企业内部控制基本规范》，要求上市公司可聘请具有证券、期货业务资格的会计师事务所对内部控制的有效性进行审计。而房地产行业在我国国民经济中为支柱产业，上下游相关产业链很长，在资源、人才方面均属于密集型，同时又受到政策调控的影响很大，故其潜在风险很大。因此，聘请会计师事务所对内部控制进行审计，可以促进公司应对和防范风险，促进房地产上市公司持续稳定的发展。本文也正是基于房地产业的在我国国民经济中的特殊地位，以其为研究对象。 从目前对高管薪酬的研究现状上看，从内部控制质量、公司治理结构影响的研究较多，但是从会计师事务所对公司内部控制评价的角度研究内部控制审计质量对高管薪酬影响的较少。 制定科学的高管薪酬激励约束机制可以促进企业战略的制定以及企业内部控制制度的制定和实施。本文选取2014~2016年A股房地产上市公司为样本数据，从公司治理的角度，分析内部控制审计质量对高管薪酬的影响，从而为上市公司设计有效的内部控制制度，提高内部控制审计质量提供经验证据和相应建议。
system: [SYSTEM]

user: 2. 理论分析和研究假设 内部控制审计质量是一个抽象的概念，不可进行量化研究。孙宝厚 [ 1 ] (2008)从狭义和广义两个角度对内部控制审计质量做了阐述，其认为审计质量从狭义上来看，表现为被审计单位信息的真实情况与审计师所报告出来的信息的吻合程度，而从广义上来看，还体现在审计报告是否能满足会计信息使用者的需求。因此本文选取会计师事务所规模、内部控制审计费用和内部控制审计意见类型作为内部控制审计质量的衡量指标。 (一) 会计师事务所规模和高管薪酬 Lennox [ 2 ] (1999)等通过实证研究发现，相比于小规模事务所，大规模事务所审计后的财务报表的信息准确度更高，这也意味着大规模事务所提供了更高质量的审计。周莉，刘端 [ 3 ] (2017)研究发现，相对于“小所”聘请“大所”的公司内部控制审计对会计盈余质量的影响更显著。张向丽 [ 4 ] (2016)研究发现，外部审计确实具有减弱高管薪酬粘性的作用：经国内“十大”会计师事务所审计的上市公司的高管薪酬粘性有所减缓，高质量外部审计具有抑制高管薪酬粘性的作用。因此本文提出假设1。 假设H1：聘请“前十大”会计师事务所进行内部控制审计的公司其高管薪酬较高。即会计师事务所规模与高管薪酬呈正相关关系。 (二) 内部控制审计费用和高管薪酬 刘运国等 [ 5 ] (2006)认为，审计师提供审计服务所收取的费用的高低可以反映其提供的审计质量的高低，收取的审计费用越多，提供的审计质量就越好，用收取的审计费用的大小来衡量审计质量的高低，对审计师提供审计的质量与企业进行盈余质量之间的相关性进行研究，研究发现：审计费用在一定程度上与盈余管理具有正相关关系。肖作平(2006)、张奇峰 [ 6 ] (2007)分别以审计费用为审计质量的替代变量进行了研究，认为审计费用越高审计质量越大。内部控制审计费用越多，注册会计师对内部控制的了解越多，执行了更多的审计程序，高管越重视内部控制，内部控制的审计质量越高，高管的薪酬越高。因此本文提出假设2。 假设H2：内部控制审计费用与高管薪酬呈正相关关系。 (三) 内部控制审计意见类型和高管薪酬 Ann L.Watkins [ 7 ] (2004)认为投资者在对公司进行投资前会对公司的财务报表信息的可靠性作出判断，注册会计师的审计报告会影响投资者作出的判断，因此审计质量就体现在注册会计师对投资者的这种影响力，意味着审计质量越高，注册会计师对投资者的影响力就越大。张继勋、何亚南 [ 8 ] (2013)从心理和行为的视角出发，采用实验研究方法检验了内部控制审计意见类型对个体投资者关于标准无保留财务报告审计意见信心的影响。黄文 [ 9 ] (2012)认为高管薪酬激励会导致企业的盈余管理行为，由注册会计师出具的外部审计意见没有能够揭示出企业的盈余管理行为。因此提出假设3。 假设H3：内部控制被出具标准无保留意见的公司其高管薪酬较高。即内部控制审计意见类型与高管薪酬呈正相关关系。
system: [SYSTEM]

user: 3. 研究设计 (一) 样本选择与数据来源 本文以2014~2016年深证A股房地产上市公司作为研究样本，剔除ST、非ST和不需要出具内部控制审计报告的公司，按照国有和非国有企业进行分类，最终得到87家样本公司的数据。本文选取的数据来源于Wind数据库和CSMAR治理结构模块研究数据库查询系统，使用STATA14.0统计软件进行实证分析。 (二) 变量设计 本文采用三类研究变量来分析内部控制审计质量对高管薪酬的影响。 1) 被解释变量——高管薪酬。 高管薪酬指的是在企业中可以做出经营决策，对企业内部控制设计和运行的有效性承担最终责任的高级管理人员。本文借鉴方军雄 [ 10 ] (2016)的做法，将公司年度财务报告中所披露的“薪酬最高的前三名高级管理人员总额”定义为高管，取其薪酬总额的对数作为高管薪酬的替代变量。 2) 解释变量——内部控制审计质量。 陈丽蓉 [ 11 ] (2015)以事务所规模作为内部控制审计的替代变量，分别研究了内部控制审计质量与公司治理以及内部控制审计质量与产品市场竞争的相关性，实现了很好的效果。陈影 [ 12 ] (2014)研究董事会特征对内部控制审计质量的影响时，分别以内部控制审计收费和内部控制审计意见类型作为内部控制审计质量的衡量变量。何勰 [ 13 ] (2013)在研究内部控制审计质量与盈余质量的相关性时，分别以内部控制审计收费和会计师事务所规模来度量内部控制审计质量。本文分别以内部控制审计费用、内部控制审计意见类型以及会计师事务所规模三类指标作为内部控制审计质量的替代变量。 本文选取会计师事务所规模、内部控制审计费用和内部控制审计意见类型三类指标作为替代变量度量内部控制审计质量。具体变量描述见表1。 3) 控制变量。 张博 [ 14 ] (2014)的研究结论表明，高管薪酬与企业业绩正相关并且内部控制质量与高管薪酬业绩敏感性之间存在显著的正相关关系。财务风险和盈利能力参考刘力涵 [ 15 ] (2015)的相关研究成果。本文将企业成长性、企业性质、财务风险、董事会规模和盈利能力作为控制变量，各变量具体描述见表1。 (三) 模型设计 为了分析内部控制审计质量对高管薪酬的影响，设计如下三个回归模型： 模型一： W A G E i t = β 0 + β 1 I C A F i t + β 2 I B R i t + β 3 L E V i t + β 4 R O A i t + β 5 S O E i t + β 6 G R O W T H i t + ε i t 模型二： 模型三： W A G E i t = β 0 + β 1 B I G 10 i t + β 2 I B R i t + β 3 L E V i t + β 4 R O A i t + β 5 S O E i t + β 6 G R O W T H i t + ε i t Table 1 Definition of variable 类型 名称 符号 定义 被解释变量 高管薪酬 Wage 排名前三位高管薪酬总额的自然对数 解释变量 事务所规模 Big10 根据2016年中注协公布的会计师事务所综合排名， 当事务所排名前十时，取值为“1”，否则取值为“0” 内部控制审计费用 Icaf 采用内部控制审计费用的自然对数 内部控制审计意见类型 Opinion 内部控制审计意见为标准无保留意见时取“0”，否则取“1” 控制变量 企业成长型 Growth 用净利润增长率度量，(当期净利润 − 上期净利润)/上期净利润 企业性质 Soe 实际控制人性质，国有取值为“0”，非国有性质取值为“1” 财务风险 Lev 资产负债率度量，负债总额与资产总额的比值 董事会规模 IBR 用独立董事人数占董事会人数的比例度量 盈利能力 Roa 息税前营业利润率度量，(净利润 + 所得税 + 财务费用)/营业收入 表1. 变量定义
system: [SYSTEM]

user: 4. 实证结果分析 (一) 描述性统计分析 本文共选取2014~2016年间中国房地产行业的261家上市公司进行研究。从表2的描述性统计中可以看出，房地产行业高管薪酬(wage)的均值约为14.80，说明排名前三位的高管薪酬总额较多，且其标准差为0.73，离散程度较大，说明三年内房地产行业高管薪酬总额出现较大变化。在内部审计质量的衡量指标中，审计意见类型(opinion)的均值和标准差较小，事务所规模(big10)的均值和标准差较为适中，而内部控制审计费用(icaf)的均值和标准差较大，说明房地产上市公司近年来所花费的内部控制审计费用较大，且行业内部存在较大差异。 (二) 相关性检验 从表3中变量间的相关性分析可以看出，房地产行业的高管薪酬总额(wage)与内部控制审计费用(icaf)显著正相关，表明当房地产上市公司内部控制审计费用(icaf)较高时，其高管薪酬总额(wage)一般也比较高；而房地产行业的审计意见类型(opinion)、事务所规模(big10)则与高管薪酬总额(wage)无明显相关关系， Table 2 Descriptive statistics of major variable 变量 观测值 均值 标准差 最小值 最大值 Wage 261 14.80 0.73 12.15 16.76 opinion 261 0.01 0.09 0.00 1.00 Icaf 261 12.69 0.51 11.51 14.58 Big10 261 0.67 0.47 0.00 1.00 IBR 261 0.38 0.05 0.30 0.60 Lev 261 0.64 0.17 0.05 0.94 roa 261 0.20 1.04 −3.97 15.02 Soe 261 0.98 0.15 0.00 1.00 growth 261 0.30 5.56 −21.06 74.96 表2. 主要变量的描述性统计 Table 3 Correlation analysis of variable 变量 wage opinion Icaf Big10 IBR Lev roa Soe growth wage 1.000 opinion 0.062 1.000 Icaf 0.375*** −0.053 1.000 Big10 0.065 −0.102* 0.100 1.000 IBR −0.148** 0.041 0.052 −0.017 1.000 Lev 0.380*** −0.035 0.286*** 0.104* −0.050 1.000 roa 0.198*** 0.547*** −0.003 −0.036 −0.044 0.000 1.000 Soe −0.003 0.013 0.098 0.023 −0.095 −0.091 0.014 1.000 growth 0.110* −0.032 0.045 0.053 −0.057 0.070 0.142** −0.021 1.000 表3. 各变量间的相关性分析 注：***、**、*分别表示在1%、5%、10%的水平上显著。 Table 4 Regression analysi (1) (2) (3) wage wage wage Opinion −0.249 (−0.44) icaf 0.433*** (5.39) Big10 0.0413 (0.49) IBR −1.530** −1.881*** 1.559*** (−2.02) (−2.62) (6.50) lev 1.568*** 1.188*** 1.559*** (6.57) (5.01) (6.50) roa 0.142*** 0.132*** 0.131*** (2.95) (3.48) (3.29) soe 0.0909 −0.109 0.0854 (0.33) (−0.41) (0.31) growth 0.00636 0.00548 0.00659 (0.84) (0.77) (0.88) _cons 14.26*** 9.330*** 14.25*** (31.12) (9.21) (31.09) N 261 261 261 Adj-R 2 0.90 0.87 0.98 表4. 回归分析 表明近年来房地产行业的高管薪酬总额(wage)受其审计意见类型(opinion)、事务所规模(big10)的影响十分小。同时也可以看出，董事会规模(IBR)、资产负债率(lev)、总资产净利率(roa)与高管薪酬总额(wage)均呈现显著正相关，说明房地产行业高管薪酬总额(wage)除了受到内部控制审计费用(icaf)影响外，还受到董事会规模(IBR)、资产负债率(lev)、总资产净利率(roa)等的影响。 (三) 回归分析结果 回归结果如表4所示。从回归的统计学角度上看，各模型调整后的R 2 相应为0.90，0.87，0.98，这表明三个模型的拟合优度都很好，解释程度较高。分析结果表明，房地产行业的高管薪酬总额(wage)与内部控制审计费用(icaf)在1%的显著性水平上呈正相关，与相关性分析结果相符，说明房地产行业的内部审计质量对高管薪酬总额存在一定的影响，且在内部审计质量的衡量指标中，内部控制审计费用(icaf)尤为突出。因此，企业要想提高高管薪酬总额，应着眼于加强内部控制审计质量，比如增加内部控制审计费用。
system: [SYSTEM]

user: 5. 研究结论与政策建议 本文以2014~2016年87家A股房地产上市公司为研究样本，利用面板模型进行实证分析，研究了内部控制审计质量对高管薪酬的影响。发现会计师事务所规模和内部控制审计意见类型和高管薪酬没有明显的相关关系，内部控制审计费用与高管薪酬呈显著的正相关关系。即内部控制审计质量越高，高管薪酬越高。区分企业性质进行分样本检验时仍然支持本文提出的假设。 根据以上结论，本文提出四点建议：一是房地产上市公司要重视内部控制，建立有效的内部控制制度，定期报告内部控制运行情况及出现的问题；二是房地产上市公司要加大内部控制的信息披露，积极改进内部控制缺陷，提高高管工作的积极性，促进高管薪酬提升；三是聘请规模较大会计师事务所审计公司的内部控制，积极与其加强沟通，建立更加有效的内部控制制度和高管薪酬激励制度；四是监管机构要加强对公司内部控制审计的执行力度，对没有披露内部控制审计报告的公司采取相应的惩罚措施。