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

user: 4. 结论与建议 本文运用博弈理论，以中央政府进行干预作为影响因素对地方政府与排污企业关于环境治理的博弈策略进行对照探究，并得出下列结论： 第一，在中央政府不发挥作用时，地方政府基于发展当地经济的内在动力与排污企业的追求利益最大化达成共识，为二者牺牲生态效益换取经济增长的“共谋”行为提供了动机，在这种动机诱导下地方政府对企业的排污治污行为消极监管，监管方消极的管理促使企业想要牟利的动机增加，最终会导致中央政府关于环境治理的政策落实不到位，极大地破坏了社会生态环境。第二，在中央政府对地方政府与排污企业的博弈施加影响后，地方政府受到上级政府的监督，迫于问责的行政压力，其寻求经济与社会效益最大化的目标转变为提高社会生态效益，因此地方政府会加大企业污染环境行为的监管力度和惩罚力度。地方政府的选择使得企业倾向选择对污染进行“治理”的策略，这个博弈策略的采取会反向推动政府降低监管的概率，减少监管成本的可能性支出。同时污染被治理使得地方政府与排污企业获得一定程度上的奖励，中央政府虽然付出了绩效奖励，但也达到了其追求社会效益最大化的目的，由此三者都达到了自己的最优策略。 基于上述分析和结论关于环境治理对中央政府提出以下几点建议： 1) 实行环境管理中央集权，构建“中央主导–地方监管–企业执行”的环境治理体系，完善“胡萝卜+大棒”奖罚制度。通过研究发现早期的分权式管理结合大力发展经济的阶段性特征使得地方政府在促进当地经济发展的基础上与排污企业达成“共谋”，使得社会环境效益受损。在中央政府干预后地方政府和排污企业混合策略博弈达到纳什均衡表明中央政府采取“大棒加胡萝卜”的激励机制对打破地方政府与用排污企业“共谋”倾向有显著作用。因此，中央政府应当加强完善对环境污染行为追责惩处的“大棒”和对治理污染进行鼓励补贴的“胡萝卜”的激励体系。 2) 细化传统的政绩评价体制，将更多元因素纳入地方政府绩效考核体系。晋升考评体制应当降低传统GDP的权重，同时纳入产业创新、社会生态、绿色GDP等环保指标 [ 12 ]，促使地方政府产生严格监管的外部压力和激发提高生态效益的内生动力。 3) 完善关于对地方履职和企业治污的监察体系，引导社会进行第三方监督。现有制度下地方政府是否严格按照中央政府的要求进行监管和排污企业是否严格依据排污标准进行治理等信息都是未公开的，由此产生的信息不对称会弱化中央政府所获取信息的真实性，因此要加快健全生态环境的监管体系进程，同时积极引导公众、社会媒体和环保部门等第三方充分发挥监督作用，为中央政府监察破坏环境的不法行为提供有力保障。 4) 提供污染治理补贴，鼓励排污企业进行治污技术创新。中央政府可以设立治污专项资金用以企业进行资源技术创新，提高资源生产力，实现污染治理技术突破。资源利用率的提高和治污成本的降低能极大提高企业进行污染治理的概率，使得污染排放问题得到有效改善。
system: [SYSTEM]

user: 摘要 随着当前高等教育的发展，高校大学生就业形势日趋严峻，为增加就业机会，大学生实习人数节节攀高，“大学生实习”一度成为话题热点，受社会广泛关注。大学生实习是一种角色转换的过程，可以将理论知识转换到实际工作中，以便夯实自身专业知识，积累丰富的工作经验，增加自身未来就业竞争力和就业机会。为此，本研究以川内大学生实习满意度为切入口，详细剖析了其实习满意度影响因素如自身、学校、企业、社会四个层面，并给出相应的建议，旨在提升大学生实习满意度，增强其综合素质，最终促进大学生就业。
system: [SYSTEM]

user: 1. 研究背景 据统计，2021届毕业的高校大学生总人数预计909万人，同比增加35万，再创历史新高。在日趋严峻的就业形势下，大学生实习人数节节攀高，但“免费劳动力”、“廉价大学生”、“实习权益得不到保障”等实习问题接踵而来。为此，展开关于川内大学生实习满意度的研究，探讨影响大学生实习满意度的因素，以期提高大学生实习满意度，促进就业。 通过文献分析法发现：李世辉、龙思远(2017年)认为我国现行的大学生实习机制存在高校、产业、政府和学生四大主体自身作为不到位且缺乏联动的问题，提出应该引入第三方实习平台，建立高效的“五体联动”机制以保障和提高大学生实习实践能力 [ 1 ]。王志梅、龚青、李骏婷、骆南峰、程延园、葛晓蕾(2017年)认为大学生实习是体验式学习，并对目前关于大学生实习的概念、测量、影响因素、实施效果、作用机制等五方面的相关研究进行了系统梳理和评析，同时表示未来更多的研究应当关注实习领域，探索其特殊的作用机制以及影响实习效果的关键因素，从而为实践中有效发挥作用提供一定的指导 [ 2 ]。丁小浩、马世妹、朱菲菲(2018年)基于目前自主实习和学校统一安排实习的混合方式，对其影响因素和就业关系进行探讨，认为学校统一安排的实习更有利于就业对口，而大学生自主实习更能提高就业起薪，但两种实习方式本身并不影响整体就业满意度 [ 3 ]。丁小浩、魏哈、乐康、约翰·摩根认为，在大学毕业生供过于求的劳动力市场背景下，大学生被迫加紧提前学习，而在实习上花费大量的时间和精力对大学学习有害，因此通过研究表明，提高大学教学质量为学生过度实习提供了可行的选择 [ 4 ]。侯义安认为，成功的实习可以起到学校、学生与行业之间“双赢”的三角伙伴关系作用，尤其有利于学生未来职业规划和发展，而大学生实习期望和认知差距影响大学生实习，因此从性别和多元智力的角度探讨学生对实习的期望，以期事先避免潜在的差距 [ 5 ]。 目前国内外有很多关于大学生实习问题的调查研究，但是大部分侧重于对实习问题的归纳或研究，缺乏对大学生这个实习主体的实地调研。为此，本研究综合运用文献分析、访谈和问卷调查的方法，侧重从大学生的角度去探讨其实习满意度的影响因素，旨在提高大学生实习满意度，从而促进高校大学生充分就业。
system: [SYSTEM]

user: 2. 研究对象与方法 2.1. 研究对象 本文以川内高校大学生为研究对象，采用线上填写问卷的调查方式，总共收集了349份问卷。共回收有效问卷349份，有效回收率为100%。其中性别比例男生为136人(39.24%)，女生为213人(60.76%)；年级比例为大一65人(17.44%)，大二150人(43.60%)，大三73人(21.22%)，大四61人(17.74%)；学校类型比例为普通本科317人(90.70%)，“985”、“211”高校32人(9.30%)，调查对象具体数据如表1所示： Table 1 Basic information of college students under investigatio 项目 类别 样本数 百分比 性别 男 136 39.24% 女 213 60.76% 年级 大一 65 17.44% 大二 150 43.60% 大三 73 21.22% 大四 61 17.74% 学校类型 普通本科 317 90.70% 985、211高校 32 9.30% 表1. 被调查大学生基本情况 2.2. 研究方法 1) 文献分析法：通过分析国内外学者关于大学生实习问题的文献，结合川内当下大学生实习的现状，探究目前川内大学生实习存在的问题。 2) 访谈法：采用线下的方式对不同高校的大学生进行访谈，总计访谈44人，通过实习意愿以及实习满意度等问题，了解了当前大学生实习的一些情况，并加以总结分析，为编制封闭式问卷提供依据。 3) 问卷调查法：基于文献分析与访谈的基础上，自编大学生实习满意度的问卷调查，采用打分制对大学生实习满意度影响因素进行相关性调查。首先从大学生实习动机入手，然后从个人、企业、学校、社会四个维度研究大学生实习满意度的影响因素，其中个人层面包括个人学历、专业对口性、个人兴趣爱好、个人未来发展等方面；企业层面包括公司发展前景、员工福利及实习薪资、企业知名度、实习工作含金量、就业地区、工作氛围等方面；学校层面包括学校知名度及特色专业、学校对于实习的要求、对于学生实习的支持力度三个方面；社会层面包括国家实习政策、实习相关的法律法规、大学生实习平台三个方面。最后收集一定的建议，以便为提高大学生实习满意度提供参考意见。 2.3. 信效度分析 根据表2信度分析结果可以看出，在影响实习满意度的因素上，总体的标准化信度系数为0.914，根据删除项后的克隆巴赫系数可以看出都是小于总体的。因此，影响实习满意度因素的维度不需要调整。而因为总体的标准化信度系数为0.914，接近1，可靠性很高，所以，实习满意度的影响因素信度非常好。 本次问卷的KMO检验的系数结果为0.945，靠近1，说明问卷的效度好。根据球形度检验的显著性也可以看出，本次检验的显著性无限接近于0，所以问卷的效度较好。 Table 2 Four-dimensional confidence test on internship satisfactio 影响实习满意度的四维度信度检验 影响实习满意度的因素 删除项后的标度平均值 删除项后的标度方差 修正后的项与总计相关性 平方多重相关性 删除项后的克隆巴赫系数 标准化后的α 影响实习满意的个人因素 64.09 169.402 0.786 0.636 0.775 0.914 影响实习满意的学校因素 84.74 278.314 0.802 0.643 0.805 影响实习满意的企业因素 61.05 162.216 0.819 0.718 0.76 影响实习满意的社会因素 88.13 308.32 0.791 0.643 0.843 表2. 影响实习满意度的四维度信度检验
system: [SYSTEM]

user: 3. 调研结果分析 3.1. 大学生实习总体满意度分析 实习作为大学生走出校园、接触职场的一环，对于大学生能否适应职场生活、找准未来职业发展定位发挥着巨大的作用。如表3所示，对于有实习规划并且已经实习过的201位大学生，53.62%的大学生对于自己最近一次的实习是比较满意，14.01%的大学生非常满意，29.47%的大学生一般满意，仅有2.9%的大学生不满意。所以，大部分实习过的大学生实习满意度较高。 Table 3 Overall satisfaction of college internshi 排序 选项 小计 比例 1 非常满意 29 14.01% 2 比较满意 111 53.62% 3 一般满意 61 29.47% 4 不满意 6 2.90% 5 非常不满意 0 总计 201 表3. 大学生实习总体满意度 综上所述，大学期间，很大一部分大学生有实习意愿和实习规划，且付诸行动的大学生的实习动机和实习满意度存在正相关，实习动机越强烈，影响实习满意度的因素越重要，实习满意度越高。 3.2. 影响大学生实习满意度的因素分析 3.2.1. 影响大学生实习满意度的总体因素分析 根据表4可知，影响实习满意度的四大因素的平均大于3.5，说明各影响因素与大学生实习满意度的相关性较强。其中，影响实习满意度的企业因素平均值高达3.83，说明在四大因素中企业因素与大学生实习满意度的相关性较高。 3.2.2. 影响大学生实习满意度的自身因素分析 根据表5与收集数据时设定的打分原则可看出，越靠近5意味着与大学生实习满意度的相关性越强，越靠近1意味着相关性越弱。影响实习满意度的个人因素平均值大于3.1，说明以上因素与大学生实习满意度有相关性。其中，个人未来发展、专业对口性、个人学历、个人经验与实习满意度的相关性在个人维度里较强，而个人政治面貌和个人家庭社会关系的相关性较弱。因此，大学生应该规划好个人未来发展计划，以选择合适的实习职位。 Table 4 Analysis of the overall factors affecting college students’ internship satisfactio 影响实习满意度的因素 均值 影响实习满意度的个人因素 3.52 影响实习满意度的学校因素 3.65 影响实习满意度的企业因素 3.83 影响实习满意度的社会因素 3.74 表4. 影响大学生实习满意度的总体因素分析 Table 5 The influence of the individual level on internship satisfactio 影响实习满意度的个人层面 均值 标准偏差 专业对口性 3.73 1.105 个人学历 3.70 1.043 个人获得的资格证书 3.50 1.068 个人兴趣爱好 3.59 1.062 个人经验 3.67 0.999 个人未来发展 3.75 0.978 个人政治面貌 3.17 1.132 个人实习年级及实习期间 3.45 0.989 个人的家庭社会关系 3.24 1.105 个人实习成本 3.43 1.063 样本量 349 表5. 影响大学生实习满意度的自身因素分析 3.2.3. 影响大学生实习满意度的学校因素分析 根据表6可以看出，影响大学生实习满意度的学校因素平均分大于3.5，与大学生实习满意度有一定相关性。其中，学校的知名度以及特色专业、学校对于学生实习支持力度在该层面与大学生实习满意度的相关性较强，因此，学校应重点打造自身特色专业，加强社会、学生、企业三方的交流，完善大学生就业指导机制，培养与社会需求接轨的人才。 Table 6 Analysis of school factors affecting college students’ internship satisfactio 影响实习满意度的学校层面 均值 标准偏差 学校知名度以及特色专业 3.79 0.995 学校对于实习的要求 3.62 1.018 学校对于学生实习支持力度 3.70 0.993 校友资源 3.50 1.079 样本量 349 表6. 影响大学生实习满意度的学校因素分析 3.2.4. 影响大学生实习满意度的企业因素分析 由表7可知，影响大学生实习满意度的企业因素平均分也大于3.7，是4个维度里面平均分最高的，也就是说企业层面的因素与大学生实习满意度具有较强相关性。其中，实习工作含金量、实习薪酬及员工福利、公司发展前景、实习氛围在企业层面相关性很强，均值高于3.9。在前面实习动机的问卷调查中，发现了实习薪酬与福利在实习动机中占比较为靠后，而到实习满意度的影响因素时，实习薪酬及福利与实习满意度的相关性极强，两个结论相悖，通过分析得到大学生实习过程中对于自我发展、基本生活需求较为强烈，更关注实习薪酬问题。可见大学生在实习过程中渴望得到企业的回报、尊重和认可，并且得到较为满意的实习薪酬也能让他们对公司更有归属感、对工作更有责任感。因此，企业应优化实习管理制度，最大程度满足大学生实习诉求，提高大学生实习满意度。 Table 7 Analysis of enterprise factors affecting college students’ internship satisfactio 影响实习满意度的企业层面 均值 标准偏差 公司发展前景 3.92 0.929 实习薪酬及员工福利 3.92 0.907 企业知名度 3.83 0.942 实习能否转正 3.77 0.960 实习职位 3.69 0.995 就业地区和就业地点 3.84 0.985 交通便利性 3.78 0.964 工作氛围 3.90 0.956 实习工作内容含金量 3.93 0.926 工作时长 3.72 0.898 样本量 349 表7. 影响大学生实习满意度的企业因素分析 3.2.5. 影响大学生实习满意度的社会因素分析 据表8可知，影响大学生实习满意度的社会因素平均分也大于3.7，说明该层面的影响因素与大学生实习满意度具有相关性。根据访谈和文献分析结果表明目前大学生实习权益保障机制不够健全，大学生权益自我保护的法律意识也较为薄弱，因此，应健全大学生实习平台、完善实习政策和法律法规、提高大学生权益保护法律意识，从而打造良好的实习环境，提高大学生实习满意度，促进社会就业。 Table 8 Analysis of social factors affecting college students’ internship satisfactio 影响实习满意度的社会层面 均值 标准偏差 国家实习政策 3.74 0.960 实习相关的法律法规 3.73 0.904 大学生实习平台 3.74 0.950 样本量 349 表8. 影响大学生实习满意度的社会因素分析
system: [SYSTEM]

user: 4. 对策建议 4.1. 企业应完善实习管理体系，优化实习过程 企业是大学生实习的主要客体，当企业给大学生提供实习机会的同时，实习生对于企业的长远发展也有着重要价值。首先，对于企业来讲，如今人才市场竞争激烈，招纳实习生可以为企业的发展储备人才资源，占领市场先机。其次，实习生可以为企业创造更高的价值。根据影响大学生实习满意度的企业因素的分析结果可以得到，大学生实习所看重的并非完全是薪酬福利，也侧重于关注实习工作的含金量。如果企业能够将大学生实习的工作安排好、管理好、指导好，完善企业内部的实习生管理制度，让大学生在实习过程中有更高的收获感与满意度，这样更能发挥大学生在实习过程中的作用，为企业创造更高的价值。 4.1.1. 增加实习薪酬，提高实习待遇 薪酬是对于员工提供的工作和服务所做出的各种直接和间接的回报，对于保障员工个人基本生存需要和自身发展需要、激励员工具有重要作用。企业作为大学生实习的客体，在大学生实习过程扮演着不可或缺的角色，面对目前实习薪资和福利较低、实习薪酬不公平的问题，应完善实习薪酬管理体系，调整实习薪酬结构，在企业盈利的情况下，合理地增加实习薪资和福利，为大学生实习提供比较充分的物质基础，满足其基本生活需要和自身发展需要，从而更好地留住实习人才。实习生作为企业的一员，面对自己实习期间的所付出的劳动，应该得到一定的、公平的回报，以满足大学生基本生活、尊重、自身发展的需求，从而提高大学生实习满意度。 4.1.2. 增加实习工作含金量，提升实习生实习的获得感和满意度 面对目前“实习生打杂”、“实习跑腿”等实习工作内容含金量不高的问题，企业作为实习客体，应该做好人力资源管理以及实习生工作规划安排。首先，企业应据职位空缺填补计划选择最为合适的实习生；其次，注重实习生岗前培训的质量，尤其提升岗位培训和技术能力培训，以便实习生掌握更具实践性和专业性的能力，更快地适应实习环境和实习岗位，进入实习角色；再者，在大学生实习期间，为其选择最为合适的实习老师，指导实习工作，解决实习难题，从而有效地提高大学生实习工作的参与感和获得感；最后，安排合理地实习工作内容，最大程度地增加其实习工作内容的含金量，让实习生有较强的实习获得感，以便提升大学生实习满意度，也能为企业培育人才。 4.1.3. 签订合理的劳动协议，注重保障实习生个人权益 在法制社会下，面对“实习权益受到侵犯”、“实习权益得不到保障”等问题，企业作为与保护大学生实习权益息息相关的一方，应该秉持高度的社会责任和法律意识，与实习生签订合理的劳动协议，尽可能保护实习生个人权益，以树立良好的企业形象，留住优秀的实习生，从而为企业未来发展提供优秀的人才资源库。 4.2. 高校、学生与企业加强联系，实现三方联动 4.2.1. 学校应加强理论教学与社会实践的联系，以便与社会实际需要人才接轨 学校作为人才培育中心，在解决实习问题、缓解就业压力上也发挥重要作用，应根据社会发展需求不断调整不同专业人才的培养方案，优化专业课程理论和实践教学设计，以有效迎合社会实际人才需求。据问卷调查可知，学校知名度和特色专业与大学生实习满意度有一定相关性，所以，高校还需要不断发展自身，打造国家特色专业，增加学校专业多样性和特色性，以提高学校知名度，获得更多更优质的实习资源，从而提高大学生实习满意度。 4.2.2. 学生自身注重夯实专业知识基础，并提升个人实习权益保护意识 面对目前存在的实习问题和日益增大的就业压力，大学生自己也应该具备忧患意识，改进自身能力的不足，增加自身就业竞争力。在学习上，注重夯实个人专业知识基础，提升个人专业力；在工作上，积极参与校内外实习、实践等，积累更为丰富的实习、实践经验，提升自身实践和工作能力；在生活上，保持积极健康的生活态度，多了解实习政策及法律法规，提升个人权益保护意识，以保护自身实习权益。 4.2.3. 加强高校、学生、企业的交流，做到三方联动 高校、企业、学生都在大学生实习过程中扮演着不可或缺的角色，为此，提高大学生实习满意度，促进大学生就业，必然需要加强高校、学生与企业之间的联动。高校在三方联动中必然扮演着沟通桥梁的作用，对于学生，它可以通过人才培养方案和实习指导制度帮助学生更好地参与实习过程；对于企业，它可以加强与合作企业的沟通，最大程度地为合作企业提供合适的实习人才。总之，相信通过加强高校、学生、企业之间的联动，可以提高大学生实习参与度、收获感、满意度。 4.3. 国家健全实习政策，提高实习权益法律保障 国家作为中心，应发挥其领导、规范作用。面对目前“实习生权益得不到法律保障”、“免费劳动力”等实习问题，更应尽快健全实习的法律法规，从法律法规层面对大学生实习权益进行强制性的立法保护，并且完善实习政策，尽可能地给高校实习生一定的实习福利，缓解大学生实习经济压力，解决大学生实习难题，促进大学生就业。同时，国家应加强和规范普通本科高校实习管理工作，特别要明晰学生、高校以及企业三方的权利义务，做好学生权益保障工作，以提高实习质量，切实维护学生、高校和企业三方的合法实习权益。
system: [SYSTEM]

user: 摘要 水泥土搅拌桩是地基处理的常用方法。在大直径钻孔灌注桩成孔施工时，为了防止坍孔，以确保邻近地层稳定，经常采用三轴水泥土搅拌桩对钻孔桩周围地层进行加固处理。但三轴水泥土搅拌桩沉桩施工因为其高压注浆和快速钻进原因，会对邻近土层产生挤压作用，这种挤压作用将导致邻近土体的位移。为此，需要根据邻近地铁区间隧道的运营维护要求，适当控制搅拌桩的注浆压力和沉桩速度，将搅拌桩的挤土不利影响降至最低程度。本文基于Mindlin理论，采用理论解析的方法，就三轴水泥土搅拌桩施工对邻近土体及邻近隧道的影响进行了公式推导和计算分析，并将计算结果和试验监测结果进行了对比。本案例表明，采用理论解析方法来预测搅拌桩的挤土效应具有快速方便的特点，在搅拌桩施工参数选择决策中，不失为一种有益的方法。
system: [SYSTEM]

user: 1. 引言 在邻近地铁区间隧道区域进行大直径钻孔灌注桩施工时，常常采用三轴深层水泥土搅拌桩预先加固钻孔灌注桩桩周地层，以防止孔壁坍塌。但是，三轴水泥土搅拌桩施工时喷浆压力大、下沉速度大，挤土效应显著，往往导致邻近地下管线和构筑物的损坏。因此，针对三轴搅拌桩施工，需要首先对其挤土效应进行预测，为选择合适的喷浆压力及沉桩速度提供依据。 对于桩基施工及其它岩土施工活动如盾构(顶管)隧道掘进等的挤土效应，预测分析的方法有很多，其中基于Mindlin基本解 [ 1 ] 的理论解析方法具有快速易学的优点，是最为常用的理论解析方法之一。至今，采用此法的工程实践积累了较为丰富的研究成果。 魏刚，徐日庆(2005)基于Mindlin基本解推导了盾构隧道掌子面荷载作用下土体的应力和位移计算公式，提出了软土隧道盾构法施工引起的纵向地面变形预测的方法 [ 2 ]；许有俊等人(2019)针对土压平衡矩形顶管正面附加推力作用，采用基于Mindlin基本解导出的土体位移公式，对地表隆起变形影响进行了分析研究，并提出了预测顶管施工过程地表隆起变形的方法 [ 3 ]；秦世伟等人(2013)从Mindlin基本解导出了静压桩沉桩施工时小孔扩张应力下半无限体内任一点处的应力和位移计算公式，并据此就静压桩沉桩施工对临近隧道的影响进行了分析研究，提出了保护隧道的措施 [ 4 ]；谷冠思(2019)针对钻孔灌注桩基桩施工全过程对临近地铁隧道土体扰动行为，采用Mindlin基本解进行了研究，在此基础上提出了控制孔壁坍塌的建议 [ 5 ]。上述学者的研究方法为三轴深层水泥土搅拌桩施工的挤土效应预测分析提供了有益的借鉴。 本文依托上海浦东创新中路附近的某地铁二号线跨线桥工程，采用基于Mindlin基本解的理论解析方法就三轴水泥土搅拌桩施工的挤土效应进行研究，推导了三轴搅拌桩施工时邻近土体的位移计算公式，并以本工程作为案例进行计算分析。
system: [SYSTEM]

user: 3. 案例分析 3.1. 工程概况 计算案例依托于上海市浦东新区创新中路某跨线桥工程 [ 6 ] [ 7 ] [ 8 ]。该桥跨越正在运营的地铁2号线区间隧道，跨线桥的桥墩基桩采用钻孔灌注桩。为防止钻孔灌注桩成孔施工时发生坍孔，先采用三轴水泥土搅拌桩对钻孔灌注桩区域的地层进行地基加固，三轴搅拌桩直径850 mm，加固深度22 m。搅拌桩和区间隧道最小距离为2.83米。桩基与盾构隧道位置关系如图5所示。 图5. 桥梁桩基和地铁隧道的位置关系 3.2. 地质参数 地基土土层自地面而下分布如下：① 填土，② 粉质黏土，③ 淤泥质粉质黏土，③ 夹层粉质黏土，④ 淤泥质粘土，⑤ 1粘土，⑤ 3-1粉质黏土，⑤ 3-2粘质粉土，⑤ 3-3粉质黏土，⑤ 3-4粉质黏土夹粉性土，⑦ 2粉砂，⑧ 1粉质黏土，⑧ 2-1粉质黏土粉砂互层。 各层土的计算参数如表1所示。 Table 1 Properties of soils for Mindlin equatio 土层 厚度(m) 重度(kN/m 3 ) 弹性模量(E/MPa) 泊松比μ ② 层粉质黏土 1.7 18.50 30 0.4 ③ 层淤泥质粉质黏土 5.1 17.8 10 0.45 ③ 夹层粉质黏土 4.1 18.3 35 0.35 ④ 层淤泥质粘土 6.7 16.8 7 0.5 ⑤ 1层粘土 9.0 18.9 35 0.30 ⑤ 3-1层粉质黏土 23.5 18.6 40 0.30 ⑦ 2层粉砂 7.1 20.1 45 0.25 表1. 土层的主要计算参数一览表 3.3. 现场量测 为了了解三轴搅拌桩施工对邻近土体的影响，施工期间在邻近隧道边缘布置了测斜管，以便获得深层土体的水平位移。测点布置如图6所示。测斜管深度为22米，距离隧道外缘0.5 m。 图6. 测斜管布置示意图 3.4. 计算结果及其分析 根据荷载情况、土层情况，按式(13)和(14)采用Mathcad软件或Matlab软件编程计算，可得测斜管水平位置沿深度的水平位移和竖向位移。计算时土层弹性模量及泊松比取值如表1所示。 3.4.1. 搅拌桩施工导致邻近土体的水平位移 计算得到搅拌桩施工对邻近土体的水平位移位移和现场实测水平位移，如图7所示。结果表明，计算值在4 m埋深处达到最大值3 mm，实测值也在大致深度处达到最大值4.7 mm。位移变化规律是：理论分析与现场实验相比，位移沿深度变化趋势相同，位移峰值点出现在较浅层土体(0~5 m)处。理论解得出的峰值位移偏小，可能是理论计算时选取的喷浆压力偏小所致。 图7. 搅拌桩施工导致邻近土体的水平位移 3.4.2. 搅拌桩施工导致邻近土体的竖向位移 计算得到搅拌桩施工对邻近土体的竖向位移，如图8所示。结果表明，在浅层土体表现为土体隆起，随深度增加，位移减小，直至反向达到峰值，随后沉降又逐渐减弱。 图8. 搅拌桩施工对土体竖向位移影响 3.4.3. 计算结果讨论 Mindlin基本解是在假定土体为各向同性的线弹性体基础上获得的，其解答适合于小变形情况。实际上土体是各向异性的弹塑性体，因此采用Mindlin公式得到的计算结果是近似的。 由于三轴搅拌桩施工对土体的扰动是在深层土体中发生的，总体上符合小变形假定。对于变形相对微小(即实际变形远小于结构如隧道的尺寸)情况，只要土层参数准确，计算结果的误差在工程允许范围内，其结论对于施工参数的初步选择仍具有指导意义，有助于快速了解注浆压力及沉桩速度对邻近土体位移的影响程度。
system: [SYSTEM]

user: 4. 结论 1) 搅拌桩施工时的挤土效应本质上就是施工过程中在大地半无限体内某一区域作用了分布荷载，因此，可以采用Mindlin解来预测邻近土体的位移； 2) 三轴搅拌桩的荷载效应计算时假定在其外包等效圆圆周上作用了径向分布力p，然后依此对Mindlin公式进行积分得到搅拌桩施工全过程中邻近土体的位移； 3) 尽管Mindlin解是在各向同性线弹性体假定基础上得到的，它仍可以真实反映搅拌桩施工时注浆压力及沉桩速度与挤土效应(及邻近土体的位移)的关系，可以帮助现场工程师快速获得施工参数控制的关键信息，这是采用复杂的三维有限元数值计算软件难以做到的。
system: [SYSTEM]

user: 摘要 秦岭地区是我国石墨矿资源丰富的主要分布区之一，目前该地区石墨资源分布及成矿机理的研究尚较薄弱。以秦岭造山带已发现的数十处石墨矿床/点/远景区的基础资料为依据，通过典型矿床野外地质调查及测试分析，运用地质综合分析方法，研究结果表明：秦岭地区石墨资源分为晶质石墨矿床和隐晶质石墨矿床两大类，呈三大矿带展布。在秦岭地区多期构造运动叠加下，区域变质作用控制形成了鳞片晶质石墨矿床；在岩浆热力与构造应力双重作用下，煤、碳质泥页岩变质形成隐晶质石墨矿床。
system: [SYSTEM]

user: 1. 引言 石墨是一种新兴矿产资源 [ 1 ] [ 2 ] [ 3 ]。天然石墨是工业化制备新型碳材料石墨烯的首选原料，全球已发现石墨矿床主要分布于北纬30˚~50˚，集中在中国、朝鲜、墨西哥、澳大利亚、俄罗斯和欧洲等地区。中国石墨储量全球第一 [ 4 ]，近年全球石墨产量前三甲为中国(67%)、印度(14%)和巴西(8%) [ 5 ]，石墨是中国出口量最大的非金属矿产品之一。 秦岭地区矿产资源丰富 [ 6 ]，是我国石墨矿资源的主要分布区之一 [ 4 ]。王力等(2017) [ 7 ]、Sun et al. (2018) [ 8 ] 对中国石墨矿床分布进行过划分，认为秦岭地区是一个重要的石墨矿床聚集区；Feng et al. (2015) [ 9 ] 研究了秦岭地区石墨矿床的基本特征，尝试对石墨矿床进行初步划分。目前秦岭造山带非金属石墨矿产的研究程度较浅，研究工作的开展大都在上世纪60~90年代。已发现石墨矿床/点/远景区数十处，类型涉及晶质鳞片状石墨和隐晶质石墨，大多数石墨矿仅开展过初步调查，有的甚至仅仅是群众举报的矿点，研究范围小且呈孤立点状研究，亟需开展区域性矿床分布特征、石墨矿床成矿机理与主控因素等方面的研究工作，服务于秦岭地区石墨资源摸家底和企业整装勘查与开发。
system: [SYSTEM]

user: 2. 秦岭地区石墨资源分布 秦岭地区含石墨矿层系主要有新太古界太华群板石山组/水底沟组、元古界秦岭群雁岭沟组与宽坪群、韩舞–奥陶系、下志留统、中泥盆统和上石炭统。秦岭地区石墨矿床类型有区域变质型晶质石墨矿床和接触变质型隐晶质石墨矿床，大鳞片晶质石墨矿床主要分布在调查区的东部，分布范围较广，以豫西南朱阳关–柳泉铺大型石墨矿为其典型代表；隐晶质石墨矿床主要分布在调查区西部，以陕西省凤县煤沟–贯沟石墨矿为其典型代表。受秦岭东西向区域性构造的控制，秦岭地区石墨矿产资源划分为三大矿带，大致呈近东西向展布(图1)。 2.1. 华北板块南缘基底型晶质石墨成矿带 主要赋矿岩石主体为新太古界–元古界华北板块的结晶基底。 ① 新太古宙太华群板石山组/水地沟组石墨矿群西起陕西蓝田，经临潼历山、华县太华山、潼关，东入河南，经灵宝、萧山、熊耳山、鲁山至舞阳。一般构成东西西北弧形古变质岩带，长约450千米，宽约50千米。包括小秦岭区(包括太华山、老牛山)板石山组和鲁山区水底沟组两套岩石组合，主要赋矿岩石为石墨片岩、石墨混合片麻岩、含石墨透辉斜(二)长片麻岩、含石墨黑云斜(二)长片麻岩等(表1)。 图1. 秦岭造山带构造格架与主要石墨矿分布 Table 1 Geological characteristics of key graphite deposits in Qinling area (The deposit location can be seen in Figure 1 矿带 序号 矿床 名称 地理 位置 矿床 类型 含矿 地层 矿层 时代 矿石类型 矿体形态 华北 板块 南缘 基底 型晶 质石 墨成 矿带 1 豫西南朱阳关–柳泉铺大型石墨矿 河南省内乡县 晶质 石墨矿 秦岭群雁岭沟组、 宽坪群 古元 古代 宽坪群为钙质晶质石墨 片岩、黝帘石晶质石墨片岩，雁岭沟组含石墨铁质 石英岩、二云石英片岩 矿体呈透镜状、层状、似层状，矿体沿走向及倾向变化大，连续性差产状与围岩一致，随地层褶皱而弯曲。平均品位4%左右，厚度1.33~38 m，一般在10~34 m左右。 2 潼关善车峪石墨矿 陕西省潼关县 晶质 石墨矿 太华群 新太 古代 含石墨斜长片麻岩、 含石墨长石石英岩 矿体呈似层状，有5个矿体， 长35~200 m，厚20~87 m。品位3%~5%。 3 西安骊山石墨矿 陕西省西安市 晶质 石墨矿 太华群板石山组 新太 古代 石墨片麻岩、 石墨片岩 矿体呈似层状，矿体呈不规则长条状，出露长1500米，宽50~300米；矿体为一透镜体，延深150~300 m；矿体总体产状倾向310˚~340˚，倾角40˚~85˚，品位3%~15%。 南秦 岭中 部晶 质石 墨成 矿带 4 十亩地–秧田坝石墨矿 陕西省佛坪县–洋县 晶质 石墨矿 上流水店组 早志 留世 含石墨黑云母石英片岩、 含石墨二云石英片岩、 石墨黑云母石英片岩 (石墨矿) 矿化带长18.3 Km，宽100~400 m， 属变成矿床。共圈定矿体129个， 呈凸镜体。矿体平均品位3.5%~4.5%。 5 留坝县青桥河石墨矿 陕西省留坝县 晶质 石墨矿 三河口组 泥盆纪 含石墨大理岩 赋存于大理岩与侵入岩接触带，矿体及矿化体17个，矿体规模小，呈条带状，单个矿体长25米，厚0.2~1米，平均固定碳20.34%。 西秦 岭山 间盆 地群 隐晶 质石 墨成 矿带 6 户县涝浴银洞沟石墨矿 陕西省西安市户县 隐晶质 石墨矿 草凉驿组 石炭纪 石墨片岩、 含石墨泥灰质绢 云母片岩 矿体成透镜状，沿北东东–南西西方向继续出露，四个矿体。长51~253 m，厚0.74~2.4 m。固定炭含量66.5~80%，含硫0.022~0.918%。 7 眉县齐镇石墨矿 陕西省眉县 隐晶质 石墨矿 草凉驿组 石炭纪 黑色石墨片岩 呈透镜体状，厚约2米，长度不明，固定炭含量约80%左右。 8 凤县煤沟–贯沟石墨矿 陕西省凤县 隐晶质 石墨矿 草凉驿组 石炭纪 石墨板岩、 黑色碳质板岩 夹石墨板岩 石墨矿呈层状、似层状、透镜状等，矿体沿走向、倾向变化较大，厚度较不稳定变化较大，石墨矿厚度0.4~2.3 m。品位76.88%~89.31%，平均为84.64%。 表1. 秦岭地区重点石墨矿床地质特征表(矿床位置见图1) ② 古元古代秦岭群雁岭沟组石墨矿群分布于天水、宝鸡唐藏张家庄、太白县沪家垣、西峡县、内乡青山、马山口及桐柏县老湾等地区。雁岭沟组以富含石墨为其主要标志，岩性主要为石墨片麻岩、含石墨铁质石英岩、二云石英片岩夹钙质晶质石墨片岩、含石墨条带大理岩、含石墨厚层状大理岩，锤击具有极明显的银灰色痕迹，手搓表面光滑，污手严重。以朱阳关–柳泉铺大型石墨矿为典型代表。 ③ 古元古代陡岭岩群大沟组石墨矿目前仅见于河南省淅川县五里梁，矿石以含石墨片麻岩、石墨片岩为主，矿体(层)形态多为层状、似层状，产状与围岩片理及区域构造线一致。 ④ 新元古代宽坪(岩)群石墨矿群石墨矿群分布在河南省西峡、内乡、镇平等县，矿石以钙质晶质石墨片岩、黝帘石晶质石墨片岩为主，矿体呈层状、似层状、透镜状，矿体沿走向及倾向变化大。 2.2. 西秦岭山间盆地群隐晶质石墨成矿带 包括西秦岭区域一系列沿商丹缝合带近东西向断续分布的上石炭统煤系隐晶质石墨矿床，重要的煤系石墨矿床有西秦岭凤县煤沟–贯沟石墨矿、眉县铜峪石墨矿、户县涝峪银洞沟石墨矿、太白县靖口油房沟石墨矿、武山县石墨矿等。含矿地层主要为石墨板岩、石墨片岩等。为隐晶质石墨，矿体呈似层状，透镜状，含固定碳一般65%~90%，属变成石墨矿床(表1)。 2.3. 南秦岭中部晶质石墨成矿带 ① 寒武–奥陶系石墨矿群主要分布在陕西省佛坪一带，矿石主要为含石墨透闪长英变粒岩、含石墨黑云母二长片麻岩，矿体呈似层状产出，平均固定碳含量4.41%。 ② 下志留统上流水店组石墨矿群主要分布在佛坪县十亩地–洋县秧田坝一带，矿石主要以含石墨黑云母石英片岩、含石墨二云石英片岩、石墨黑云母石英片岩(石墨矿)为主。矿体呈凸镜体，平均固定碳含量3.5%~4.5% (表1)。 ③ 中泥盆世三河口组石墨矿群主要分布在留坝县青桥河，矿石主要为片状含石墨石英岩、含石墨黑云斜长片麻岩、含石墨大理岩，矿体赋存于大理岩与侵入岩接触带，矿体呈条带状，矿体规模较小，平均固定碳含量20.34%。
system: [SYSTEM]

user: 3. 秦岭地区石墨矿产成矿主控因素 3.1. 多期构造运动叠加区域变质作用控制大鳞片晶质石墨矿床形成 自太古代秦岭地区经历了多期构造运动的叠加改造，多次区域变质作用下发育了晶质石墨矿床，经受改造的次数越多，石墨晶体发育相对较好。赋矿岩石来自华北板块结晶基底的往往发育基底型大鳞片晶质石墨矿，赋矿岩石来自古生代盖层的常形成细鳞片晶质石墨矿。 秦岭地区区域变质晶质石墨矿石原岩为含碳质沉积建造 [ 10 ] [ 11 ]。无机来源的碳质以为来自原岩中的碳酸盐矿物主，碳酸盐经漫长区域变质作用使其中碳酸盐主要发生CaCO 3 (方解石) → CaO (硅灰石) + CO 2 、CaMg (CO 3 ) (白云石) + 2SiO 2 (石英) → CaMgSi 2 O 6 (透辉石) + CO 2 而脱碳形成碳气。生成的以CO 2 为主的碳气，被H + 还原，与其结合，主要经CO 2 + 4H + → C + 2H 2 O反应转变成单质的碳。Kreulen (1987) [ 12 ] 指出流体包裹体中CO 2 和CH 4 的混合物并不稳定，会发生CO 2 + CH 4 → 2C + 2H 2 O化学反应形成石墨，直到两种成分中的一种被完全消耗。Ortega等(2010) [ 13 ] 通过研究英国Borrowdale石墨矿床矿化通道中与石墨结晶相关的石英晶体内的流体包裹体，发现石墨形成于富含蒸汽的CH 4 + H 2 O流体。 研究结果表明，在涉及CO 2 和CH 4 的反应中，石墨从C-O-H流体中形成，可能有几种方式：① CO 2 + CH 4 → 2C + 2H 2 O；② CO 2 + 4H + → C + 2H 2 O；③ CO 2 → C + O 2 ；④ O 2 + CH 4 → C + 2H 2 O。经历上述水化反应后，流体从500℃冷却到400℃，流体中的碳过饱和而析出 [ 13 ] (图2)。后期岩浆侵入引起围岩温度增高，富含钾、钠、硅质的岩浆溶液交代含矿层，它受到以钾盐和硅化作用为主的混合岩化作用的影响。分散细晶石墨聚集再结晶，形成大鳞片状晶质石墨。 图2. 英国Borrowdale石墨矿床的C-O-H流体系统 [ 13 ] 3.2. 岩浆热力与构造应力双重作用控制煤变质型石墨矿床形成 3.2.1. 秦岭造山带多次岩浆热侵位为煤层受热变成石墨提供了充足的热动力条件 西秦岭煤系石墨矿石具有它形鳞片变晶结构，呈似层状、透镜状构造，矿层强烈揉皱，内部小构造发育。I号和II号矿层固定碳含量与视密度分别为67.38%与1.87和84.64%与1.90，热接触变质引起下部矿层固定碳含量升高，视密度增大。岩浆热侵位对煤层受热变质形成石墨的程度具有控制，距离越近石墨的晶体自形程度越高，晶体也越大。距离岩浆岩体较近的贯沟石墨矿石结晶程度较好，接近正六边形；吴家沟石墨矿距离岩体较远，石墨为密集细小鳞片，自形程度较差，少见六边形。 电子探针实验在热接触变质带角岩化岩石中发现大量的代表热接触变质作用的矿物及组合，出现红柱石、硬绿泥石、绿泥石、黑云母、石榴子石等，红柱石具有重要的角岩化指相意义。钻孔资料揭示地下400 m左右岩石重结晶现象显著，黑灰色石英砂岩和石英砾岩大多变为灰白色石英岩。西秦岭凤县石墨矿床中石墨的分布和红柱石密切相关，且距离岩体有一定距离。老硐中煤系石墨矿层中热液石英脉发育，石英形成于热液背景自形较差。 初步研究表明煤系经受岩浆热侵位发生接触变质作用具有层次性：a. 在岩浆岩体侵位阶段，在岩体一侧发生电气石化、云英岩化以及碱交代活动。在被接触的一侧发生局部熔融并有硼、氟、钾、钠、硅的迁移，产生硅化和电气石化；b. 在紧随其后的岩体释热、煤系围岩受热阶段，煤系围岩中泥质物变成红柱石或堇青石等，煤层则经历构造有序化过程后变成隐晶质石墨。 3.2.2. 构造应力“催化”石墨化进程，并控制石墨富集部位 秦岭造山带内纬向构造及其与其他构造体系复合处，是中生代特别是燕山期岩浆强烈而频繁活动的地方 [ 14 ] [ 15 ] [ 16 ]，秦岭受到强大的挤压和多期岩浆冲压作用 [ 17 ]。陕西凤县草滩沟–煤沟一带构造十分发育 [ 18 ]，受到华力西运动、印支运动和燕山运动以及喜山运动构造应力影响 [ 19 ] [ 20 ]，产生大量级别不等的断层，内部小断裂和揉皱十分发育，由于矿层较软容易变形且石墨是良好的润滑剂，矿层作为应力的释放点而被强烈改造且往往与断层伴生，石墨矿层在空间上呈现透镜状、鸡窝状或者豆荚状，复式褶皱系内次级背斜与次级向斜的核部成为石墨矿层富集带(图3)。 图3. 西秦岭凤县煤沟–贯沟煤系石墨矿龙王沟构造控矿剖面图
system: [SYSTEM]

user: 4. 结论 (1) 秦岭地区发育晶质石墨和隐晶质石墨两大类矿床，研究首次揭示了秦岭地区石墨矿资源分布特征，大致呈近东西向的三大矿带。 (2) 秦岭地区晶质石墨成矿主控因素：多期构造运动叠加，区域变质作用控制形成了鳞片晶质石墨矿床。 (3) 创新性的提出了秦岭地区煤系石墨的形成主要受温度与压应力力的双重作用控制。秦岭造山带多次岩浆热侵位为煤层受热变成石墨提供了充足的热条件，构造应力“催化”石墨化进程，并控制石墨富集部位。
system: [SYSTEM]

user: 摘要 研究企业的贷款违约风险不仅对金融机构解决“惜贷”问题和防范信用风险具有重要的现实意义，而且能为企业规范自身经营和改善财务状况提出有针对性的建议及措施。本文根据某机构的企业贷款违约数据对贷款违约风险进行研究，首先对原始数据进行缺失值处理、特征选择和不平衡数据处理，然后利用逻辑回归、随机森林、XGBoost和LightGBM四种机器学习方法对数据进行建模和分析并比较模型优劣，最后利用GBDT模型计算特征重要性。结果表明：1) 三种集成模型的预测效果显著优于单一模型，2) 在集成模型中LightGBM模型表现出了最好的预测性能，3) 企业的纳税情况和曾经获得的授信情况可以作为判断该企业是否会发生贷款逾期现象的重要参考。
system: [SYSTEM]

user: 1. 引言 长期以来，贷款难的问题一直阻碍着企业的发展。虽然国家出台了一系列的支持政策，但企业自身存在的经营风险以及贷款违约现象的频繁发生使得银行等金融机构越来越“惜贷”。这种情况产生的根本原因在于金融机构缺乏具体有效的信用风险评估方法和防控措施，特别是随着大数据时代的来临，基于统计学方法的传统违约风险预测模型在数据处理方面的不足日益明显，且不具备自学习能力，已经无法适应如今的要求。因此，有关机构能否获得更加高效又切实可行的信用风险评估模型成为决定其未来发展的关键。
system: [SYSTEM]

user: 2. 文献综述 随着人工智能理论的深入发展，众多学者在基于机器学习算法的信用风险评估方法上取得了丰硕的研究成果。例如，朱景文以上市公司的债券违约风险为研究对象，证明了经过遗传算法优化的BP神经网络模型在公司债违约风险评估应用中的高准确性 [ 1 ]。王晓菲等通过将AI模型与传统信用风险度量模型相结合的方法，从多角度对房地产企业的信用违约风险进行研究 [ 2 ]。李占玉将SMOTE算法与随机森林相结合构建了互联网金融公司的财务风险评估模型，并取得了良好的预测效果 [ 3 ]。潘永明等利用引入信息增益的SVM算法构建供应链中小企业信用风险分类预测模型，结果显示该模型比单一SVM模型精度提高8.97% [ 4 ]。郑建国等将PCA、SMOTE与通过网格搜索法进行参数寻优的SVM模型相结合对企业信用风险进行研究，结果表明新构建的模型具有更高的稳定性和预测能力 [ 5 ]。胡贤德等针对传统BP神经网络在实际应用中学习速度慢、易陷入局部解以及运算结果误差较大等缺陷，提出一种基于改进离散型萤火虫(IDGSO)算法的BP神经网络集成算法，并以此建立了小微企业信用风险评估模型，该模型有效提高了预测准确性 [ 6 ]。
system: [SYSTEM]

user: 3. 基本方法 本文选择使用逻辑回归、随机森林、XGBoost (Extreme Gradient Boosting)和LightGBM (Light Gradient Boosting Machine)四种模型对企业贷款违约数据进行建模研究，并根据研究结果选出预测性能最优的模型，为有关机构进行风控管理提供参考。其中，逻辑回归模型是一种单一模型，具有更加简单的算法原理和内部结构，且该模型的输出形式易于理解、建模效率较高，有十分广泛的应用；随机森林是一种由大量分类回归树组成的集成算法，其应用十分灵活，在特征较多、数据量较大的数据集上表现良好；XGBoost是一种基于GBDT (Gradient Boosting Decision Tree)的集成算法，它不仅对损失函数进行二阶泰勒展开，而且将正则化项加入到损失函数中，有效控制了模型的复杂程度并大大提高了模型的预测性能；LightGBM也是一种基于GBDT的改进集成算法，主要用于处理样本容量大且特征维度高的数据，其具有运算速度快、内存占用率低、准确率高以及支持并行计算等诸多优点。
system: [SYSTEM]

user: 5. 结论 针对当今社会众多的企业贷款违约问题，本文利用Python软件基于机器学习方法构造违约风险评估模型，并以某机构的企业贷款违约数据为例进行实证分析。首先将原始数据进行缺失值处理、特征选择和不平衡数据处理，然后将数据集划分为训练集和测试集，并分别代入逻辑回归、随机森林、XGBoost和LightGBM四种模型进行训练和预测。结果表明：综合混淆矩阵和ROC曲线的多个指标，可以发现后三种集成模型(随机森林、XGBoost和LightGBM)对分类问题的预测效果明显优于单一模型(逻辑回归模型)，并且在集成模型中LightGBM模型拥有最好的表现，其两项五折交叉验证均值均大于另外两个集成模型，这与其运算效率高、内存消耗低、支持并行计算等优点密切相关。因此，使用基于LightGBM算法的评估模型对实际贷款违约问题进行研究和预测是一个明智的选择，它能帮助有关机构更好地规避信用违约风险、加强风控措施。最后，本文利用GBDT模型对数据集中的特征重要性进行评估，发现企业的纳税情况和曾经获得的授信情况的重要性最大，分别为0.414和0.403，这两个特征的相关情况应该作为判断该企业是否会发生贷款逾期现象的重要参考。 本文的不足之处在于有效数据集较小，数据集涵盖的特征种类较少，无法充分体现出不同集成算法间的差别。同时，LightGBM模型的预测效果还有进步的空间，可以考虑将它和其他模型再次融合，进一步提升预测的准确率和稳定性。
system: [SYSTEM]

user: 摘要 以某SUV 6座纯电动汽车为研究对象，用仿真软件对车内制冷及舒适性进行研究，考虑人体模型，采用PMV-PPD和空气龄(MAA)对乘客舱内热环境的舒适性情况进行分析和评价。结果显示，车内冷空气流速、体表温度分布在舒适的范围内。通过分析评价指标的分布情况，可以清晰地了解车内舒适性分布情况，较实验具有更节约成本、缩短时效的优势。
system: [SYSTEM]

user: 1. 引言 目前，电动汽车已经越来越受群众的喜爱，但其内部乘坐热舒适性也越来为人们所关注。现在汽车为了保证驾驶乘客驾驶安全及审美疲劳，而增大汽车车窗来开阔驾乘人员的视野，但却增大了太阳辐射，影响乘坐人员的热舒适性，甚至影响驾驶安全，所以对乘客舱舒适性研究极为重要 [ 1 ]。汽车室内的舒适性主要与空气的温度、速度、质量及太阳辐射有关 [ 2 ]。热舒适性在炎热和寒冷的气候表现的最为明显，车内温度过高、过低都会影响人的舒适性 [ 3 ]。近年，有关汽车乘客舱及空调的研究比较多，大多在建立三维汽车模型的基础上对送风结构即参数进行研究，如文献 [ 4 ] [ 5 ] [ 6 ] [ 7 ] 通过建立模型，研究不同送风条件对车内的热舒适性的影响。文献 [ 8 ] 通过CFD软件建立单个通风口的地上、地下两个相同有限空间内的气流运动状态，采用不同网格长度对比分析有限空间内的气流运动，得出外界的风量和风速越大，内部空气循环的越彻底。车内热舒适不能仅考虑送风结构变化对其的影响，还应综合考虑太阳辐射和人体热生理情况，文献 [ 9 ] [ 10 ] 考虑了人体生理热，引入PMV-PPD等多种热舒适性指标对有限空间的热环境进行评估。文献 [ 11 ] [ 12 ] [ 13 ] 综合考虑了太阳辐射和人体生理热对室内热环境的影响，通过建立仿真模型以及引入评价指标对有限空间内的空气质量和热舒适性进行评估，得出室内环境的舒适性情况。 综上所述，这些研究大多侧重于风道的改造，对于综合考虑太阳辐射和人体自身因素对热环境影响的研究很少，不能精确地反映出热环境的舒适性。因此，本文对SUV 6座汽车内部的热舒适性进行模拟研究，并综合引入各种评价指标对车内的舒适性进行评估，具有一定的实用意义。
system: [SYSTEM]

user: 2. 模型的建立 2.1. 几何模型及网格划分 本文使用CATIA软件对SUV汽车和假人模型进行简化，简化过程中保留制冷吹面的风道，且风道密封良好，如图1为简化后的SUV汽车模型。然后对整车和假人模型网格划分，STAR-CCM+中，整个计算域内体网格数为850万。图2为假人及汽车网格模型。 2.2. 数学模型 STAR-CCM+中有自带的数值求解器，可以精确模拟车内的热舒适情况。车内空气流速较小，使用k-ε湍流模型，内部流体可看作三维不可压缩流 [ 14 ]；乘客舱内的空气为辐射透明介质，符合灰体辐射。求解方程如下： 连续性方程： ∂ ρ ∂ t + ∂ ρ u i ∂ x i = 0 (1) 图1. SUV汽车三维几何模型 图2. 假人及车舱网格划分示意图 动量方程： ∂ ρ u i ∂ t + ∂ ρ u i u j ∂ x i = − ∂ P ∂ x i + ∂ ∂ x j ( μ ( ∂ u ∂ x j + ∂ u j ∂ x i ) ) (2) 能量方程： ∂ ρ T ∂ t + ∂ ρ T u j ∂ x j = ∂ c p ∂ x i ( k ∂ T ∂ X i ) (3) 式中， ρ 是气流密度； μ 是动力粘度；P为压力，P；T为温度，℃；t为时间，s； x i 为x、y、z方向的张量；k为热传导； c p 是比热容。
system: [SYSTEM]

user: 3. 边界参数及评价指标 3.1. 边界参数 本文将进口设为质量流量入口，出口设为分散流出口，采用瞬时计算3900 s。在STAR-CCM+仿真软件中有自带的参数设置框。 3.2. 评价指标 3.2.1. PMV-PPD PMV表示预期平均热感觉，其将舒适状态从冷到热划分为7个等级，见表1。由于环境的舒适性因人体的生理情况不同而异。所以不适人员比例PPD被提出来表达人们对热环境的不满意程度，当PMV在−0.5~0.5时，且PPD < 15%时，人体处于舒适状态。PMV和PPD的计算如式(4)和(5)。 PMV = [ 0.028 + 0.303 exp ( − 0.036 M ) ] ⋅ [ M − W − E c − E r e s − C r e s − H ] (4) PPD = 100 − 95 ⋅ exp ( − 0.03353 ⋅ PMV 4 − 0.2179 ⋅ PMV 2 ) (5) 式中，M——人能量代谢率，W/m 2 ； W——人做的机械功，W/m 2 ； E c ——表皮及出汗的散热量，W/m 2 ； E r e s ——呼吸潜热损失量，W/m 2 ； C r e s ——呼吸显热损失量，W/m 2 ； H——体表对流和辐射热损失，W/m 2 。 Table 1 Comfort level of PMV evaluation inde PMV指标 −3 −2 −1 0 +1 +2 +3 热感觉 冷 偏冷 凉 舒适 温暖 偏热 热 表1. PMV评价指标舒适等级 3.2.2. 空气龄(MAA) 空气龄为新鲜空气在有限空间内停留的时间，可通过示踪气体实验估算其值 [ 15 ]。有限空间内某点的平均空气龄计算公式，单位：秒： τ = ∫ 0 ∞ C ( τ ) d τ C 0 (6) 式中， C 0 ——有限空间内示踪气体均匀分布时的初始浓度； C ( τ ) ——当以某种特定送风方式送入新鲜空气时测得的示踪气体浓度随时间的变化。
system: [SYSTEM]

user: 4. 模拟结果分析 4.1. 乘客舱空气流速分布情况 图3为模拟过程SUV车内空气流动状态。可以看出，风向从风口出来后直接吹响乘客，离风口较近的部位，气流密度强，对于乘客来说吹风感较猛。图中，乘客从头部到小腿出现涡旋状回流，表现为头暖脚冷，且风口直吹的地方风速较大达到3 m/s，人体表面最大风速达到2.5 m/s。总体而言，乘客头部的风速相对腿脚处较高，但因热胀冷缩原理，车内温差会尽快降到合理范围；图中可以看出人体感到舒适的风速为0.5~2 m/s，前两排乘客体表风速基本能够满足，后排乘客与C柱出风口的距离太近，使得乘客表面及肩部风速稍高。 4.2. 乘客舱成员体表温度分布情况 图4为模拟过程乘客表面的温度分布状态。图中，前排乘客腿部温度较其他乘客高，基本在38℃左右。一方面是该部分新鲜空气流速较低，换热量较少；另一方面是因为前排挡风玻璃较大，腿部吸收的辐射热较多。二排分布乘客体表温度相对均匀。一方面是因为侧窗玻璃的透射率较低，吸收的辐射热较少；另一方面是因为空气流速较均匀，能够迅速降低温度。第三排乘客除后侧头部区域温度较高，其他区域的温度均匀分布。这是因为三排乘客表面风速虽然分布较均匀，但后脑勺部位吸收的辐射热较大。可见乘客体表温度不仅与出风口温度、吹风位置有关，还与乘客舱周围玻璃透射率有关。总体而言，乘客舱处于人体舒适的温度区间内。 图3. 车内空气流线示意图 图4. 乘客表面温度 4.3. 乘客舱平均温度分布 图5为车舱截面温度分布，可以看出乘客舱温度总体处在26℃附近，但车舱轮廓处温度要比车舱内部温度稍高。这是因为边界处直接吸收大量的太阳辐射，且又与外界高温环境进行对流换热，冷空气到达车舱边界时的流速很低，对流换热系数较小，不能及时给边界降温，故使轮廓温度高于内部温度。图中前挡风玻璃与仪表盘的夹角处温度过高，这是由于三维仿真中的网格尺寸大于该处的实际几何尺寸，因而产生计算死角，所得的温度结果失真。 图5. 车内横截面温度 4.4. PMV-PPD评价指标分析 图6中，PMV值基本分布在−0.6~0.5之间，乘客处于热舒适状态。因为后两排车窗玻璃的透光率比前排玻璃低，乘客所吸收的太阳辐射热少；且对流换热较好，冷空气均匀分布在体表上，所以后两排乘客基本处于热舒适状态。前排乘客除大腿和腹部位置的PMV值较高，其余处的PMV值处在[−0.5， 0.5]的舒适区间内，这是因为前挡风玻璃的透光率较高，大腿及腹部吸收的辐射热最多，并且该处的冷空气流速较低，换热量较少，热量难被带走。 图6. 乘客表面PMV分布 图7. 乘客表面PPD分布 图7，乘客表面PPD分布基本低于15%，可认为乘客的满意程度比较高。前排乘客除大腿及腹部的PPD值高于15%，其余位置处于满意的状态，这是因为乘客的大腿和腹部吸收大量的太阳辐射，且该部位的空气流速低，换热量较少；因后两排乘客的肩颈直对出风口，风速感较强，且从车窗处吸收了部分的辐射热，使人感到极不舒适。 4.5. 空气龄(MAA) 有限空间内，空气龄在40 s内被视为可接受的舒适范围。如图8，乘客表面的MAA分布在25 s内，说明有限空间内的新鲜空气循环的比较顺畅，乘客感觉比较舒服。这里是乘客表面的平均MAA已经达到了舒适的状态，但是车内局部位置MAA较高，如图9，车舱前部的MAA较高，这是因为车舱前部吸收过多的太阳辐射，且风速循环较慢，换热量低，导致温度太高，空气干燥，MAA较高。总的来说，后两排乘客的MAA要比前排乘客的稍低。这是因为后排出乘客接收的新鲜冷空气较多、空气流速较均匀，并且吸收的辐射热较少，因此MAA要比前排的稍低，感觉比较舒适。 图8. 乘客表面MAA分布 图9. 乘客表面MAA分布
system: [SYSTEM]

user: 5. 总语 本文通过STAR-CCM+建立了三维SUV 6座车舱模型，引入人体模型和太阳辐射模型模拟车内热环境在高温夏季静置后分别以不同的运行方式启动时车内的热舒适性情况。总得来说，车内及乘客表面温湿度分布较均匀且舒适，稳定时平均温度在26℃附近，但前排乘客的大腿及腹部位置吸收较多的辐射热，使得平均辐射温度值偏高、空气比较干燥，这是因为前挡风玻璃的透射率较大，使PMV-PPD的值超出了舒适的范围。同时，后排出风口的风速偏高，后排乘客的肩颈部吹风感较强，该区域的乘客PPD值偏高。因此，出风口温度过低以及空气流速过快都会影响乘客的乘坐热舒适性，为了满足热舒适性要求，汽车空调系统不仅需要合理控制出风口风速，还需要控制阳光进入乘员舱的辐射量。