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

user: 4. 农业科研单位人才建设对策与建议 4.1. 加强顶层设计，建立长远人才培养计划 首先，要加强单位发展定位方面的顶层设计。要按照国家发展目标和政策引导并结合生产实际需求，在科研方向、人才引进、人才培养等方面综合制定未来发展规划，并按照总体规划，合理安排布局，以“按需设岗、按岗聘人”为原则，设置合理的岗位，同一学科领域要考虑人才梯队效应，既要形成老中青搭配合理的人才梯队，又要形成科研创新人才和成果转化人才配合默契的人才团队，建立领军人才、骨干人才及科研辅助搭配合理的人才结构，实施人才分散引进并梯队培养，形成相互补充、各司其职的立体发展体系，使聘用到相应岗位上的人才能够最大限度地发挥其才能，避免人才资源的浪费。 其次，要做好科研人才源头把控。在选人方面要将热爱农业科研工作作为首选条件，改变以往只看学历、外表、笔试和面试成绩的狭隘用人观，要结合备选人员的学习经历、专业课成绩、专业作品及毕业论文质量、学校各方评价等综合考虑。要确保把正确的人选到正确的岗位上来，使其人岗合一、专业对口、乐于工作，充分发挥人的主观能动性是科研成功的关键。 第三，在人才培养方面应制定合理的人才跟踪培养机制，进行分类培养，科研人才和成果转化人才需要统筹兼顾 [ 11 ]。一是积极做好团队领军人才或学术带头人的培养、遴选和支持工作，这类人员除在学术上有所建树外，还得培养其组织管理能力；二是做好骨干人才的培养，通过制定绩效考核制度、搭建科研创新和成果转化平台等，鼓励和支持青年骨干人才在项目申报、国际交流、成果转化、科技实践活动中勇挑重担；三是做好科研辅助人员的培养，锻炼和培养他们的常规试验技能操作速度和质量。此外，还要建立必要的容错机制，鼓励青年人才大胆创新，敢于挑战的精神，如在攻坚克难、先行先试未知探索等方面虽然付出了极大努力，但因客观条件限制、技术发展不同步等因素导致科研失败，要给予工作肯定和鼓励，使其放下思想包袱，始终保持积极进取的科研态度。 4.2. 加强思想建设，解决内因问题 思想是管总的，要扭住“思想建设”这个总开关，培养踏踏实实的科研队伍，凝聚科研力量，营造良好科研氛围，农业科研单位要不定时地给科技人员做思想工作，通过教育塑造科学道德，包括学术道德、团结合作、勇于创新、乐于分享等。良好的科研道德是科技人员开展农业科研工作的重要内在素质保证 [ 10 ]。要教育他们尊重科研规律，科学研究的基本任务就是探索、认识未知。既然是探索未知，就包括能达到预期目标和达不到预期目标两个结果，科研工作者必须遵循这个规律。因此，在科研上既要有明确目标、合理规划、科学实验、久久为功的进取精神，也要准备着接受实现不了预想结果的现实，科研人员要端正心态正确面对，决不能因为达不到目标就心浮气躁、全盘否定科研全过程，要有持续探索的精神，不断总结、归纳、调整科研方向和方法，进而在后续过程中保持取得实效的勇气和决心。 4.3. 凝聚科研精神，形成科研文化 科研精神、科研文化是科研单位的内在力量，是可持续发展的根基和动力。求实、创新、协作、奉献、自律是科研人员的应有素质。科研需要实事求是的态度和坚定不移的意志，需要不断探索、积极开拓的闯劲，需要善于学习、团结协作的团队意识，需要甘于寂寞，坐得住冷板凳的献身品质，需要遵守职业道德、诚实守信的约束力量。科研单位要营造和推动这样的科研文化，将这些品质在每一个科研人员的身上得到传播和固化，成为科研人员挥之不去的精神名片。科研人员要经常出现在田间地头，用步伐去丈量数据，用深思来探究路径，用吃苦奉献来践行科研精神。 4.4. 创造优质条件，稳定人才队伍 创新是科研发展的原动力，真正创新性的科研是需要在宽松环境下，科研人员兴趣所致的一项持久、深入的研究事业。科研单位要积极采取有效措施，通过人才管理体制改革，创造优质条件，形成尊重知识、尊重人才、鼓励创新的良好氛围。首先要保证科技人员的薪酬待遇，薪酬是对人才所做贡献最直接的奖励。因此，农业科研单位要结合各自的能力与特点，通过制定相关有效的绩效考核制度进行薪酬改革。在现行事业单位岗位绩效工资制度的基础上，着重修改奖励性绩效制度，可与科研进展、取得的成绩、所获成果、效能转换等各种因素挂钩，给予科研人员一个坚实的保障。其次要不断改善科研人员的工作与生活环境，为大家排忧解难，使其能够潜心研究。一是在单位大的研究框架和发展定位内，要充分落实科技人员的科研自主权，鼓励大家从事他们兴趣内的科技创新工作，并且要积极营造和谐开放、宽容失败的科研环境 [ 8 ]。二是积极搭建各种学术交流平台，包括加强国际交流和合作，活跃学术氛围，让科技人员在“百花齐放、百家争鸣”的学术争辩中凝练科学问题，了解行业发展动向，紧盯行业前沿技术水平，促进科技创新。三要适当减轻科研人员的心理和生活压力，在其住房条件、子女入学等方面提供一定的帮助与支持，尽可能创造一个能让科研人员专心搞科研的环境。 4.5. 优化人才评价机制，激发人才创新活力 2018年2月26日，中共中央办公厅、国务院办公厅印发了《关于分类推进人才评价机制改革的指导意见》，明确指出“人才评价是人才发展体制机制的重要组成部分，是人才资源开发管理和使用的前提 [ 12 ]。农业科研单位要结合实际突破以往评价方式，由单一评价向综合评价转变，由内部评价同市场评价和社会评价结合，发挥多元评价主体作用研究切实可行的评价机制，重点要考虑五个方面的问题。一是突出“品德”评价。把品德作为科研人才评价的首要内容，倡导高尚的职业道德和从业操守，强化诚信体系建设，坚决杜绝学术造假、剽窃论文等失信行为，建议在道德和诚信方面设置“一票否决”评价，以强硬的学术手段来维护科研领域的风清气正。二是设立综合评价标准。坚持以能力、实绩、贡献评价人才，调整以往唯学历、唯资历、唯论文等“一刀切”的评价方式，构建考察各类人才的评价体系，既要突出专业性，也要兼顾创新成果和实际付出，鼓励人才在不同领域、不同岗位作出贡献、求实奋进。三是灵活确定评价周期。探索实施灵活多变的评价周期，打破年度评价的传统方式，将科研过程、科研预期结合起来，将短期评价和长期评价结合起来，针对领军研究人才、青年人才的科研项目来设置考核周期，鼓励持续研究和长期积累。四是拓展人才评价渠道。要盯住科研贡献这个关键，凡是对科研创新、科研进步做出贡献的科研人员都应被尊重，打破户籍、地域、身份、人事关系等限制条条框框，特别是针对引进的海外人才，要建立特殊评价机制，不能搞“一刀切”。五是建立能上能下的评价机制。针对取得特殊科研成果的专家，在尊重自身价值、给予其必要的荣誉和待遇的基础上，更要加强后续考核问效，以此来促进科研人员再接再厉，百尺竿头更进一步，始终保持科研热情，同时也要防止出现躺在功劳簿上睡大觉、奋斗一时、享受一世的情况发生，一旦出现名不副实的状况，要及时纠正，维护良好的科研生态。
system: [SYSTEM]

user: 5. 结语 相对于发达国家，我国的农业发展水平还不高，规模化、集约化、现代化的农业发展模式尚未形成，农业发展势在增速提质，农业科技人才是发展的关键，因此农业科研单位需要从农业发展的真实需求方面综合考虑人才建设的方方面面。相信，随着农业科研院所对自身职责定位、科研发展规划、人才队伍建设、体制机制完善等方面逐步清晰、规范，我国的农业科研发展一定会再上新台阶，为农业现代化的实现加油助力。
system: [SYSTEM]

user: 摘 要 本论文基于“推–拉”模型理论，采用文献分析法，以老年人群对冬季避寒旅游消费需求为切入点深入探讨中国老年人群避寒旅游消费需求驱动力的影响因子，构建理论模型及计算各影响因子所占比重，然后进一步分析、筛选与识别各驱动力因子。通过阅读、整理与研究大量文献资料得知：气候舒适度和自然景观资源两大因子是中国老年人群避寒旅游消费需求的核心驱动力因子；可支配收入、旅游地文化、交通条件、空气质量、基础设施、旅游动机、生态环境、自然地理位置八大因子是关键驱动力因子；高患病率和死亡率、高纬地区供暖不足的危机与环境污染、闲暇时间、网络媒体宣传与亲朋好友推荐、身心健康现状、绿色健康饮食、身心健康的保持、生活便捷八大因子是一般驱动力因子。通过构建理论模型试图提升避寒旅游资源潜力与竞争力，为建设可持续发展的避寒旅游目的地提供可操作的科学依据。 关键词 :避寒旅游，消费需求，驱动力因子，中国老年人群 Copyright © 2019 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. 引言 自2000年以来全球60岁及以上老年人口占总人口比重已达10%，这意味着全世界从此开始进入了人口老龄化时代。2014年，全球已有60多个国家或地区已步入老龄化社会，人口老龄化已成为一个全球性现象 [ 1 ] 。中国于1999年进入了“未富先老”的老龄化社会。截至2018年底，中国60岁及以上人口达2.49亿，占全国总人口的17.9%，65岁及以上人口达1.67亿，占全国总人口的11.9%，如表1。可见，中国老年人口比重大，增长快。 目前，中国已成为全球最活跃的旅游市场和旅游业发展最快的地区，当然，旅游开发离不开气候，气候资源是旅游业中不可或缺的一种宝贵的旅游资源。而近年来我国冬季严寒、雾霾、雪灾和冻雨等极端天气导致老年人呼吸道、心血管等疾病频发，入院率与死亡率倍增 [ 2 ]，据统计，2017年诊疗人数达 Table 1 Proportion of total population of different age groups in China from 1990 to 201 1990年 2000年 2010年 2018年 60岁及以上 8.56% 10.33% 13.26% 17.9% 65岁及以上 5.51% 6.96% 8.87% 11.9% 表1. 1990~2018年中国不同年龄段人口占总人口比例 资料来源：《中国人口统计年鉴》。 818，311万人次，入院人数达24，436万人，据2018年中国统计年鉴得知，我国的医疗卫生投资额为5243.3亿元，国家预算额为811.7亿元，且2013~2017年全国居民人均消费支出医疗保健费用呈上升趋势，如图1。其中，我国65岁以上老人年均医疗费用已达到1072.3元，远超于其他年龄段的年均医疗费用。因此，老年人群所需要的巨额医疗保险费用给社会家庭带来巨大的医疗、养老经济负担，老年人的身体健康与生活质量，也对中国的经济以及医疗保障面带来了极大挑战。 图1. 2013~2017年全国居民人均消费支出医疗保健费用 在这样的背景下，老年人群更加憧憬天然避寒气候环境，使得老年人对冬季避寒旅游需求剧增，避寒旅游研究成为了新时代的新课题，且赋予了新的使命。过亿老年人口的生活质量与幸福指数尤其是健康问题倍受社会的关注，其中相当部分离退休职工的体质较好，退休金高，有充足的闲暇时间，因此，这部分老年群体为发展避寒旅游提供了庞大的市场主体，避寒旅游蕴藏巨大市场机遇和消费潜能 [ 3 ] 。我国避寒旅游研究起步晚，理论研究重物轻人，因此，完善老年人群避寒旅游消费需求驱动力因子研究显得十分必要且迫切。
system: [SYSTEM]

user: 2. 相关概念界定与旅游驱动力文献回顾 2.1. 相关概念界定 2.1.1. 驱动力 驱动力一词源于物理上的概念，指驱动轮得到发动机动力时产生对地面的作用力，地面反力称为驱动力。世界上第一台蒸汽机，就是利用水沸腾产生的高压蒸汽能量转化为机械功来驱动机器运转；三轮汽车是将煤气燃烧产生蒸汽来驱动三轮汽车行驶；现代汽车是汽油在发动机气缸内和空气混合燃烧，把化学能转化成动能带动发动机做功，驱动汽车行驶；电动车即电力驱动车，是以电池作为能量来源，将电能转化为机械能来驱使电动车动起来，从而为人类服务。通过以上案例联到旅游者进行旅游活动的驱动力。 通过查阅文献得知，旅游驱动一词最早由国内学者保继刚教授提出，他认为城市旅游发展的动力是城市主题 [ 4 ] 。本文将旅游驱动力界定为旅游者群体参与旅游活动的原因与动力，着重剖析影响旅游者群体参与旅游活动的驱动力因子 [ 5 ] 。目前，驱动力研究正广泛运用于旅游研究领域，这对探索旅游驱动力研究具有十分重要的突破意义。 2.1.2. 避寒旅游 本论文将避寒旅游的概念界定为：以中高纬度带为主体的旅游者在每年冬季(北半球，12月至次年2月)以南方低纬地区，如三亚、海口、云南等为旅游活动场地，以温暖舒适的气候为主要吸引物，在旅游目的地进行至少为期一周以上且满足其休闲度假、健康疗养需求的旅游活动，同时也是一种利用冬季天然的气候资源与自然的生态环境的低碳旅游。 2.2. 旅游驱动力文献回顾 迄今，对旅游驱动力的研究主要集中于旅游发展驱动力机制、城市旅游驱动力、乡村旅游驱动力等方面。保继刚教授是较早探讨城市旅游动力机制的国内学者，他于1998年提出珠海市旅游发展驱动模型 [ 4 ] ；1999年彭华教授首次构建了旅游发展动力系统的结构模型，将旅游动力系统分为四大子系统，分别是引力系统、需求系统、支持系统、中介系统，并运用到对汕头市旅游动力系统的分析 [ 6 ] ；学者丁金胜依据旅游系统结构理论，将青岛市城市旅游动力系统分为吸引、需求、保障、决策五个子系统 [ 7 ] ；陈小春(2016)提出我国传统村落旅游发展的因素，将驱动力要素分为资源驱动力、经济驱动力、市场驱动力、社会驱动力及行政驱动力五个方面，每个驱动力要素由许多驱动力因子组成，进而提出传统村落旅游发展的驱动力结构 [ 8 ] 。
system: [SYSTEM]

user: 3. “推–拉”模型理论及其在旅游领域中运用价值 推拉理论(Push-pull theory)的起源可以追溯到19世纪末，最早由英国学者列文斯坦(E. Ravenstein，1889)在有关人口迁移行为研究中提出，他的观点为推拉理论的发展奠定了基础 [ 9 ] 。20世纪50年代末博格(D. J. Bogue， 1959)系统地阐述了推拉理论，主要研究人口学中劳动力迁移的推拉因素以及迁移行为的产生是推力和拉力相互作用的结果 [ 10 ] 。20世纪60年代美国学者Lee在前人研究的基础上，更加系统、完善的提出了推拉理论，将个人因素引入解释框架 [ 11 ] 。而后英国学者格拉汉姆M·S·戴恩(Dann G. M. S.)于1977年较为系统地提出了关于旅游活动的“推–拉”动机模型，他认为“推–拉”理论能有效解释旅游者流动机制，从而将该理论引入旅游研究领域 [ 12 ] 。至1979年，学者Crompton将游客出游动力总结为两大因素，即推动和拉动因素，随后推拉理论开始被广泛应用到旅游研究领域 [ 13 ] 。 在旅游研究领域，Dann (1977)提出的“推–拉”模型理论是研究游客动机和出游行为关系，该理论将影响游客出游行为的因素分为“推”和“拉”两种因素，并进一步解释推力和拉力的概念，确立了“推–拉”模型理论在旅游研究中的重要地位 [ 9 ] [ 14 ] 。其中，在研究游客行为、游客动机、旅游地吸引力、旅游目的地形象等方面推拉理论得到了广泛运用。之后其他学者在应用时又不断丰富和发展了“推–拉”模型理论的内涵，如美国学者Muzaffe Uysal等对“推力”因素和“拉力”因素赋予了明确的含义，即“推力”因素是指旅游者个体内在的旅游需求，主要与旅游者的需求动机及旅游者自身因素相关，如旅游者出游是为了逃离日常生活、休息放松、健康、社会交际等；“拉力”因素是指某一特定旅游目的地和有形的吸引物对旅游者的吸引使其离开居住地到目的地旅游的作用力，主要与吸引物和目的地自身属性相关，如旅游目的地的资源、环境、特色等等 [ 15 ] 。总的来说，相关研究已显示出以“推–拉”模型理论所具有的旅游运用价值，但鲜见将该理论运用于避寒旅游驱动力研究中。因此，本论文尝试将“推–拉”模型理论应用于中国老年人群避寒旅游消费需求驱动力的影响因子研究，既有其必要性，又有可行性。
system: [SYSTEM]

user: 4. 中国老年人群避寒旅游消费需求驱动力影响因子统计 老年人群作为避寒旅游活动的主力军，为了解老年人口避寒旅游消费需求，以中国知识资源总库CNKI收录的文献为对象，以“避寒旅游”为篇名，剔除不相关论文，共检索出15篇论文 [ 2 ] [ 3 ] [ 16 ] - [ 28 ]，然后采用手动统计方法对论文进行梳理、分析与总结归纳，得出以下中国老年人群避寒旅游消费需求驱动力的影响因子，见表2。 Table 2 Statistics on the influencing factors of driving forces in the winter resort tourism of consumption demand for the elderly in Chin 研究地域 作者(时间) 研究的主要内容 影响因子 四川省 攀枝花 李建(2009) 阳光避寒疗养度假模式 气候舒适度、空气清洁度、环境优美、交通条件好、医疗保障度与生活便捷度、城市基础设施好、文化资源、旅游资源、较强的旅游动机、较强的支付能力 仁和区 罗成德(2008) 分析避寒旅游资源 特殊的气候环境、空气质量、良好的交通条件、区位条件较好、观光旅游资源、基础设施较完善 贵州省 贵州省 李忠燕(2018) 气候舒适度及旅游潜力分析 气候舒适度、较好的交通、区位条件、文化旅游资源、旅游规划与宣传力度 赤水 夏兆林(2018) 分析避寒旅游气候 条件 优质的气候条件、空气质量、丰富的旅游资源、良好的生态环境、可支配收入、交通条件便利、旅游需求或动机、环境污染影响大 云南省 西双版纳 林锦屏(2013) 评价西双版纳避寒 旅游潜力 气候舒适度、自然景观、人文景观、体质好、退休金高、充足的闲暇时间、冬季供暖污染、较好的区位条件、老年人入院率与死亡率倍增、老年人避寒保健需求 元谋 李鹏坤(2016) 避寒养老旅游产业 开发策略 气候资源优势、充裕的储蓄、交通便利、区位优势、绿色饮食产业、文化旅游资源、保持身心健康 西双版纳 罗赟敏(2017) “避寒”养老旅游者 行为特征 气候温暖、空气质量、自然生态环境、经济收入、闲暇时间、旅游服务设施、口碑推荐、自然风貌、风土人情、健康养生、旅游动机、严寒气候诱发和加重老年人疾病 西双版纳 陈莹(2018) 避寒旅游地生态文化透视 天然气候舒适度、优越的地理位置、多样的生物资源、民族文化资源、旅游者动机、高入院率与死亡率、供暖能耗与污染、体质较好、退休金高、充足闲暇时间 西双版纳腾冲 保继刚(2018) 气候因素对度假地型第二居所需求的影响 宜人的气候、可支配收入、交通条件、闲暇时间、景观资源、文化、旅游服务配套设施 广西壮族自治区 梧州 冯琼兰(2010) 梧州避寒旅游资源的评价 冬季气候舒适、山清水秀、空气清新、生活便捷度、基础设施完善、良好旅游区位、较强的旅游接待能力、较丰富旅游资源、绿色农产品 广东省 湛江 罗松英(2017) 旅游气候舒适度分析与评价 气候舒适度、滨海旅游资源 阳江 吴长春(2018) 避寒旅游气候舒适度分析 气候舒适度、旅游资源丰富(包括自然、人文景观)、旅游动机 海南省 海口 李宁(2017) 居民对避寒旅游影响的感知与态度 气候舒适度、旅游资源丰富、经济高速发展、环境质量、宣传推广、基础设施、文化内涵 中国 林锦屏(2003) 中国南方十一座旅游名城避寒疗养气候旅游资源评估 优越气候舒适度、空气质量优、绿色生态环境、绿色食品、交通发达、服务设施较佳、北方冬季采暖、体质较好、低温严寒季节导致高死亡率、保持健康延长寿命 黄宇宸(2016) 中国避寒旅游目的地城市气候舒适度 气候舒适度、旅游资源吸引力、空气质量、绿化覆盖率、生活水平提高、交通发达、低温严寒季节导致多种疾病发生、旅游动机 表2. 中国老年人群避寒旅游消费需求驱动力影响因子统计 通过文献分析法得知，国内学者对避寒旅游的研究多以分析避寒旅游目的地的气候舒适度、评价避寒旅游资源、研究避寒旅游的最佳时间及避寒旅游地的发展潜力等，甚少与微观层次避寒旅游者的需求相关联，对老年人群避寒旅游消费需求驱动力研究更是罕见。老年人群作为避寒旅游活动的主力军，如何改善老年人避寒旅游条件与旅游环境、提高避寒旅游的幸福指数以及更好的促进避寒旅游地可持续发展等问题成为旅游开发经营者、政府部门和学术界普遍关注的问题。
system: [SYSTEM]

user: 6. 结论与讨论 6.1. 结论 1) 在影响中国老年人群避寒旅游消费需求的驱动力因素中，旅游资源和旅游者个体因素是中国老年人群避寒旅游消费需求的核心驱动力因素，且旅游资源对老年人群避寒旅游消费需求影响更为显著；区位条件和旅游地接待能力是关键驱动力因素，且区位条件占比高于旅游地接待能力，这两大驱动力因素对旅游者的影响较大；保持身心健康、冬季高患病率和死亡率、网络媒体宣传与口碑推荐、高纬地区供暖危机与环境污染是一般驱动力因素，对旅游者影响相对较小。 2) 在影响中国老年人群避寒旅游消费需求的驱动力因子中，气候舒适度和自然景观资源两大因子是中国老年人群避寒旅游消费需求的核心驱动力因子；可支配收入、旅游地文化、交通条件、空气质量、基础设施、旅游动机、生态环境、自然地理位置八大因子是关键驱动力因子；高患病率和死亡率、高纬地区供暖危机与环境污染、闲暇时间、网络媒体宣传与亲朋好友推荐、身心健康现状、绿色健康饮食、保持身心健康、生活便捷度八大因子是一般驱动力因子。目前，国内学者甚少把性别、客源地、受教育程度、职业等因子考虑进去，系统探讨这四个驱动力因子对老年人群避寒旅游消费需求的影响。 6.2. 讨论 1) 避寒旅游是一种生态旅游，它有利于旅游活动的主体——人的身心健康，郭来喜(2017)首次倡导增加“科教、康体、环境”新旅游三要素，把健康提到议事日程上，由此构成“食、住、行、游、购、娱、教、体、环”旅游九大要素 [ 29 ] 。当代老年人，冬季较少接受太阳光的照射，多有骨质疏松和抑郁症等多种疾病，导致其身体和心理存在多种健康隐患。而避寒旅游拥有的特殊气候环境，气候温暖适宜、阳光充足，当老年人沐浴冬季柔和的阳光，不仅能防止佝偻、骨质疏松等疾病，增强人体新陈代谢和抗病能力，还会使他们心情愉悦，心态良好，从而达到老年人追求身心健康、延年益寿的重要目标。 2) 然而，当旅游者大量涌入避寒旅游地时，也产生了其他旅游地在旅游旺季存在的问题，如产生大量的废弃物、汽车尾气，严重污染旅游地环境，威胁旅游地避寒旅游资源和生态环境，影响避寒旅游地的生命周期。同时，交通拥堵、缺乏宣传、基础设施配置不合理等问题也同样突出。因此，老年人群作为避寒旅游活动的主力军，为了满足他们避寒保健需求，躲过冬季的恶劣气候，必须保护避寒旅游资源与优化生态环境，改善旅游地交通条件，加强避寒旅游地的服务和基础设施建设，加大避寒旅游宣传力度，从而改善老年人避寒旅游条件与旅游环境步伐。 3) 笔者在大量阅读有关避寒旅游、“推–拉”理论、旅游驱动力等文献的基础上，通过梳理、分类、归纳与详细研究文献资料，尝试构建了中国老年人群避寒旅游消费需求驱动力理论模型，从理论视角探讨了中国老年人群避寒旅游消费需求驱动力因子的不同影响程度。在此构建的理论模型基础上，将在接下来的研究中进行实地调研获取相应数据来验证中国老年人群避寒旅游消费需求驱动力理论模型，通过实证理论模型，识别与鉴定影响因子是否有偏差，再进一步进行分析与调整，筛选出最显著影响中国老年人群避寒旅游消费需求的驱动力因子，从而提出具有科学性与针对性的建议与对策，促进我国避寒旅游目的地的可持续发展。
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 ] 。车联网系统往往通过在车辆仪表台安装车载终端设备，实现对车辆所有工作情况和静、动态信息的采集、存储并发送。车联网系统一般具有实时实景功能，利用移动网络实现人车交互。车联网获得的数据可以用于对车辆行驶轨迹的分析，通过其了解道路交通状况，合理规划外出行程，缓解交通压力等。但是由于地理环境影响以及仪器故障等原因使得车联网采集的数据往往会出现异常，所以我们需要对车辆行驶轨迹数据进行分析与挖掘，对异常数据进行智能校正 [ 2 ] 。
system: [SYSTEM]

user: 2. 问题分析 针对各种可能的因素会造成车联网采集的数据缺失或者数据异常，所以首先对轨迹数据做预处理。车联网在工作过程中，由于传感器所采集的数据量是巨大的，常常会造成大量数据冗余，通过预处理去掉轨迹中经纬度坐标值一直无变化的点，然后对轨迹数据中异常偏移的点进行筛选剔除 [ 3 ] 。 预处理去除错误以及冗余数据之后，由于车辆在行驶过程中经纬度坐标变化比较连续，基于运输车辆路线是走最短路径，所以根据电子实时地图按照最大概率进行路网匹配，通过编写程序代码，可得到补全之后车辆完整的行驶路线图。
system: [SYSTEM]

user: 3. 建立模型 3.1. 轨迹数据预处理 运输车辆的运动轨迹在地理空间通常可以由一系列按照时间序列分布的坐标点来表示，它们的轨迹可以由 P 1 → P 2 → P 3 → ⋯ → P n 表示，其中每个点由对应的经纬度坐标和一个时间戳组成，由 P = ( x ， y ， t ) 表示。由于各种可能的因素会造成被记录的点数据缺失或者数据异常，首先应对轨迹数据做预处理 [ 4 ] 。 3.2. 去除数据冗余 传感器采集的数据量是巨大的，但是往往也会造成大量数据冗余，例如在分析车辆轨迹时，可以去除掉经纬度坐标值一直无变化的点，在这一段时间内说明运输车辆是未运动的。以此筛选后的数据在地图上进行标注时会极大的提高效率。 3.3. 样条插值 计算运输车辆需要经纬度坐标变化比较连续，减少毛刺现象的产生，对存在缺失的异常值采用样条插值的方式。数学上将具有一定光滑性的分段多项式称为样条函数，具体的说给定区间 [ a ， b ] 的一个分划： Δ : a = x 0 < x 1 < ⋯ < x n = b 如果函数 S ( x ) 满足在每个小区间 [ x i ， x i + 1 ] ( i = 0 ， 1 ， ⋯ ， n − 1 ) 上 S ( x ) 是m次多项式且有 S ( x ) 在 [ a ， b ] 上具有 m − 1 阶连续导数，则称 S ( x ) 为关于划分 Δ 的m次样条函数，其图形为m次样条曲线。 3.4. 噪点过滤 在空间路径不是特别精确的情况下，噪点过滤算法的作用是删除或者修正轨迹中位置存在有错误的点，如图1所示，存在偏差过大的点，需要在开始数据挖掘任务前先将其过滤，主要从三种过滤算法入手：均值过滤，卡尔曼过滤和基于启发式的异常检测模型。 图1. 路径噪声示意图 1) 均值过滤。如果某一个异常点对于它前后两个观测点位置偏移量极大时，可以采用均值过滤的方法，使车辆在单位时间内的位移距离更加合理。第K个点的位置估值可以表示为： P K = P k − 1 + P k + 1 2 采用均值过滤的方法可以去除明显不合理的噪点。不过当出现连续的噪点时，会导致估值与实际值的差距增大，会使均值滤波效果变差。 2) 卡尔曼滤波过滤。卡尔曼过滤根据测量值和运动趋势得到轨迹估值，卡尔曼过滤基于轨迹符合线性模型的假设，根据运动物体转移概率矩阵的关系，对异常值的采集点预测并行校正 [ 5 ] 。 卡尔曼滤波通过系统输入输出的观测数据对系统状态进行最优的估计，能够实现实时预测的卡尔曼滤波动态轨迹预测的状态方程为： X ( k + 1 ) = A ( k ) X ( k ) + T ( k ) W (k) 其中 X ( k ) 表示为系统的状态向量， A ( k ) 表示为状态转移矩阵，可以用来描述前一时刻到当前时刻下的运动状态转移方式，其中 T ( k ) 为干扰转移矩阵， W ( k ) 表示运动模型的系统状态噪声。 观测方程为： Z ( k ) = H ( k ) X ( k ) + V (k) 其中 Z ( k ) 为观测向量，描述了在某个时刻的观测值， H ( k ) 在本文中为一维的观测矩阵， V ( k ) 为运动估计过程中的产生的观测噪声。 卡尔曼滤波算法的核心是运用递归算法来达到最优状态估计的估计模型。主要是利用前一时刻的估计值和现时刻的观测值来更新当前状态的变量的估计，基于前k个观测值得到在k时刻下的最优状态估计，计算最小方差计算的公式表示为： M = E [ ( x i − x ′ i ) ( x i − x ′ i ) T ] 对时间和状态分别进行更新。时间的更新过程是根据前一时刻最优状态估计预测出当前时刻的状态，同时也要更新预测状态的协方差，时间更新方程为： X ( k + 1 ， k ) = A ( k ) X ( k ， k ) Z ( k + 1 ， k ) = H ( k ) X ( k + 1 ， k ) 当预测出轨迹点后，根据观测值和预测值通过观测更新方程来推测出最优估计点。 3) 基于启发式的异常检测模型。 第一条：该类模型基于GPS坐标点与其后续点之间的时间间隔和距离间隔来计算出运输车在每个点的行驶速度。 如果一个路段相邻两个采集的经纬度坐标值的距离与时间间隔的比值，表示为： V K = S t 当 V K 明显大于某个阈值即高于200 km/h时，认为为异常点。 第二条：如果相邻两个GPS数据小于10米的点被认为是异常点。 第三条：GPS数据的速度低于5千米每小时和相邻两条GPS数据的时间大于1分钟为异常点。 第四条：如果一条轨迹最多有4个GPS数据则被认为是异常数据并删除。 3.5. 基于隐马尔科夫的地图匹配算法 地图匹配就是通过车辆的行驶轨迹和电子地图数据库中的道路网进行比较，并结合车辆原来的运行趋势，在地图上找出相近的路线，然后将实际定位数据映射到直观的显示地图上，如图2所示。 图2. 地图匹配原理 Step 1：确定参数： HMM的5个基本因素，由状态变量数，观测变量数，状态转移矩阵，观测概率矩阵，初始状态概率向量组成 [ 6 ] 。 1) 确定观察变量。从车联网采集的数据中得到的位置信息(包括经度，纬度)； 2) 隐常状态。此时在路网上的运输车的实际所在位置； 3) 观测概率。观测经纬度离旁边道路近，那么真实点在这个路段的概率越大 [ 7 ] ，假设服从高斯分布设有： p ( z t | r i ) = 1 2 π σ z e − 0.5 ( ‖ z t − x t ， i ‖ σ z ) 2 σ z = 1.4826 m e d i a n t ( ‖ z t − x t ， i ‖ ) 4) 状态转移概率。在两个时间点候选路段的转换概率： p ( d t ) = 1 β e − d t / β d t = | ‖ z t − z t + 1 ‖ − ‖ x t ， i − x t + 1 ， i ‖ | ， β = 1 ln ( 2 ) m e d i a n t (dt) 5) 状态转移矩阵A： a i j = P ( q t + 1 = r j | q t = r i ) = e − r i j r j 是时刻t + 1时所在路段， r i 是时刻t所在路段， r i j 有三种状态，当 r i j 为1时，相邻， r i j 为为无穷大时，不相邻， r i j 等于0时，在同一路段。 Step 2：Viterbi算法求解 原理基于动态规划利用概率最大找最优路径，其递推公式为： δ t ( i ) = max P ( i t = i ， i t − 1 ， ⋯ ， i 1 ， o t ， ⋯ ， o 1 | λ ) δ t + 1 ( i ) = max i 1 ， ⋯ ， i t P ( i t = i ， i t − 1 ， ⋯ ， i 1 ， o t ， ⋯ ， o 1 | λ ) 其中 i = 1 ， 2 ， ⋯ ， N 。那么定义在时刻t状态为i的所有单个路径中概率最大的路径的第t − 1个结点为： ψ t ( i ) = arg max [ δ t − 1 ( j ) a i j ]
system: [SYSTEM]

user: 4. 实验求解 图3为车联网系统采集到某辆运输车行驶轨迹的经纬度坐标所做出的行驶路线图，从图3中可以明显的看出行驶路线中存在大量的异常值，由于一辆车的速度等方面是有限的，所以短时间内不会在经纬度路线图上面有很大的变化(尖端毛刺)，所做出的路线图应该是平滑的才符合实际情况。 图3. 实际数据路线图 由于各种可能的因素会造成被记录的点数据缺失或者数据异常，所以首先对轨迹数据做预处理后。根据上述的异常点检测模型的检测方法，对轨迹数据进行噪点过滤，剔除异常点后路线图如图4所示。 图4. 去除异常数据图 预处理去除异常以及冗余数据之后，由于车辆在行驶过程中经纬度坐标变化比较连续，并且考虑到运输车行驶方向的趋势，再结合实时电子地图，采用基于隐马尔科夫模型的地图匹配算法可得到补全之后的路线图如图5所示。 图5. 地图匹配算法处理后 从图形中可以看出，地图匹配算法的补全效果的非常好，减少了许多毛刺，能较好的将缺失路段补全，非常贴合实际情况路线。
system: [SYSTEM]

user: 5. 总结 本文通过实验验证了异常点检测模型的检测效果比较好，能迅速找到轨迹数据存在漂移的地方，对异常点合理的进行剔除，并且基于隐马科夫模型的地图匹配算法根据车辆行驶方向的趋势并且同时结合电子地图信息，能大概率的将偏离数据匹配到路网上，得到的轨迹路线图效果较好。但是本处理方法也存在一些问题，利用地图匹配算法需要联网更新实时电子地图，效率偏低，对轨迹的匹配仍然存在一定的误差。
system: [SYSTEM]

user: 摘 要 为了获得L-苹果酸高产的谷氨酸棒状杆菌(C. glutamicum)，通过基因编辑，在草酰乙酸回补途径的关键酶丙酮酸羧化酶基因(pyc)前面，敲入超氧化物歧化酶基因启动子(Psod)，重组菌株命名为CZM03，通过强启动子Psod增强PYC酶活性，促进丙酮酸→草酰乙酸→苹果酸，利于苹果酸的生成。结果表明：与野生型C. glutamicum相比，重组菌株CZM03的PYC活性上调了约2.2倍；到平台期，两者的生长速率和最大生物量没有显著差异。这表明Psod强启动子的敲入，对菌的生长没有明显影响，但是PYC的酶活却有大幅度提高，这对于目的产物的积累是有利的。 关键词 :L-苹果酸，谷氨酸棒状杆菌，基因编辑，丙酮酸羧化酶，超氧化物歧化酶基因启动子Psod 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. 引言 L-苹果酸是一种重要的四碳化合物，已被美国能源部列为十二大生物基平台化合物之一，在食品、医药、化工等行业有着广泛的应用。目前苹果酸的合成主要有化学合成法、酶法和发酵法 [ 1 ] 。化学合成法虽然成本低廉，但合成的苹果酸为DL型，不易吸收，且有一定的毒性，不宜在食品与医药工业中应用。酶转化法生产L-苹果酸主要以富马酸酶转化为主，依赖高纯度的富马酸为原料，而富马酸价格昂贵。发酵法合成苹果酸，目前通常采用曲霉属菌株作为生产菌株 [ 2 ] [ 3 ] [ 4 ] [ 5 ] ，产量较高，但是生产周期长、产杂酸多、且难以对其进行遗传改造。 C. glutamicum是经典的食品性安全菌，遗传背景清晰，已广泛用于化学品的生物制造 [ 6 ] [ 7 ] [ 8 ] 。C. glutamicum在厌氧条件下，丙酮酸经过还原TCA (rTCA)途径合成苹果酸，理论上的糖酸转化率为200%，是目前所有已知L-苹果酸合成途径中最高的 [ 9 ] 。本研究采用系统生物学方法对C. glutamicum中L-苹果酸代谢途径进行分析，通过基因编辑，对野生型C. glutamicum进行代谢途径改造，在草酰乙酸回补途径的关键酶丙酮酸羧化酶基因(pyc)前面，敲入超氧化物歧化酶基因启动子(Psod)，增强PYC酶活性，促进丙酮酸→草酰乙酸→苹果酸，利于目标产物的积累。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 材料 本研究所用到的菌株和质粒见表1。 野生型C. glutamicum ATCC13032菌株、E. coli DH5α由本实验室保存；质粒DNA提取试剂盒和琼脂糖凝胶DNA回收试剂盒均购自Omega Bio-Tek公司；限制性内切酶购自Thermo scientific；T4 DNA ligase购自TaKaRa；fast Pfu PCR Mix DNA聚合酶购自北京佳兰生物科技有限公司；引物均由江苏金唯智生物科技有限公司合成；其它试剂均为国产分析纯。 用于扩增pyc基因上游同源臂的引物为CPyr-U1和CPyr-U2；用于扩增pyc基因下游同源臂的引物 Table 1 Lists of strains and plasmids used for this stud Name Description Sources C. glutamicum C. glutamicum ATCC 13032 Stored in lab C. glutamicum CZM03 C. glutamicum ATCC 13032 (Psod) pyc This study pK19ms-Psod-pyc pK19mobsacB derivative containing upstream and downstream homology arms of pyc gene and Psod This study 表1. 本研究所用菌株、质粒 为CPyr-D1和CPyr-D2；用于扩增Psod启动子的引物分别CPs 1和CPs 2；Psod敲入后的检测引物为OPs1和Ops2。具体引物序列组成见表2。 Table 2 Primers used in this stud 引物名称 引物序列(5’→3’) 退火温度 CPyr-U 1 CTATGACATGATTACAAGCTTGCTGTGGCAGTGACCAACCG 58℃ CPyr-U 2 CCCGGAATAATTGGCAGCTATAGAGTAATTATTCCTTTCAACAAGAGACC CPyr-D 1 TTTTACCCGTGTCGACTCACACATCTTCAACG 58℃ CPyr-D 2 GAGCTCGGTACCCGGGGATCCCAAAGATGGGGTAAGTCTGGCC CPs 1 TAGCTGCCAATTATTCCGGG 58℃ CPs 2 ATGTGTGAGTCGACACGGGTAAAAAATCCTTTC OPs 1 TGCAATATCAACAAAAAGGCCTA 60℃ OPs 2 GAGCTCGGTACCCGGGGATCCCAAAGATGGGGTAAGTCTGGCC 表2. 本研究所用的引物 注：引物CPyr-D 2与OPs 2共用。 发酵种子液培养基(g/L)：葡萄糖20，蛋白胨5，酵母提取物2.5，氯化钠5，NaAce 2.5。 发酵培养基(g/L)：葡萄糖20，KHCO 3 20，NaAce 5，KH 2 PO 4 0.75，K 2 HPO 4 0.75，(NH 4 ) 2 SO 4 5，MgSO 4 ∙7H 2 O 0.5，CaCl 2 0.2，MnSO 4 ∙4H 2 O 0.02，FeSO 4 ∙7H 2 O 0.005，V H 0.2 mg/L，V B1 0.2 mg/L。 2.2. 方法 2.2.1. 基因敲出载体构建 Psod启动子的敲入载体pK19ms-Psod-pyc构建：先扩增用于同源重组的上、下游同源臂片段，再扩增Psod启动子片段，然后利用Soeing PCR，获得上游同源臂、Psod启动子、下游同源臂三者的融合片段；将融合片段以双酶切的方式克隆至pK19mobsacB，获得Psod启动子的敲入载体pK19ms-Psod-pyc。 2.2.2. 敲入载体电转化C. glutamicum 制备C. glutamicum感受态细胞，制备方法参照文献 [ 10 ] 进行。敲入载体电转化C. glutamicum感受态细胞的参数设置为：1.8 kV，4 mS。具体操作步骤参照文献 [ 10 ] 进行。 2.2.3. 酶活测定 丙酮酸羧化酶酶活测定参照文献 [ 11 ] 进行。 2.2.4. 生长曲线绘制 每隔3 h取样，于OD 600 处测菌液浓度，并绘制生长曲线。 2.2.5. 残糖测定 取0. 5 mL培养好的菌液于2 mL的离心管中，加入0. 5 mL的DNS，在沸水中孵育5 min立即置于冰水浴冷却，然后在A 540 测OD值，再根据标准曲线计算出残糖，具体参考文献 [ 12 ] 。
system: [SYSTEM]

user: 3. 结果与分析 3.1. 敲入载体pK19ms-Psod-pyc的构建 以野生型C. glutamicum基因组DNA为模板，扩增目的片段，参照方法1.2.1，构建敲入载体pK19ms-Psod-pyc。结果见图1，图1(a)的1号为pyc基因上游同源臂(约500 bp)，2号为Psod启动子片段(约200 bp)，3号pyc基因下游同源臂(约500 bp)，4号为上游同源臂、Psod启动子、下游同源臂的融合片段(大小约为1200 bp)。 图1. 敲入载体pK19ms-Psod-pyc的构建。(a) 敲入片段的扩增：1：pyc基因上游同源臂，2：Psod启动子片段，3：pyc基因下游同源臂；(b) 敲入载体的酶切鉴定：1：重组质粒pK19ms-Psod-pyc，2：重组质粒pK19ms-Psod-pyc用Hind III和BamH I双酶切 将此融合片段克隆(酶切位点为Hind III和BamH I)至pK19mobsacB载体(大小为5719 bp)，得到重组载体pK19ms-Psod-pyc。重组质粒进行双酶切鉴定，结果见图1(b)，所得条带大小与预期相一致(1200 bp)，说明上游同源臂、Psod启动子、下游同源臂的融合片段已成功克隆到pK19mobsacB载体上，重组载体命名为pK19ms-Psod-pyc。 3.2. 重组菌株CZM03的鉴定 将所构建的敲入质粒载体，电转化谷氨酸棒状杆菌，筛选阳性重组菌株并鉴定，结果见图2。M为DNA分子量标记Marker，1号和2号分别以野生型C. glutamicum和重组菌株CZM03的基因组DNA为模板，以上游同源臂的上游引物CPyr-U1和下游同源臂的下游引物CPyr-D2，检测Psod启动子的敲入情况，Psod启动子大小为200 bp，上、下游同源臂大小约为500 bp，由图2(a)可以看出，Psod启动子已成功敲入。 同时，在Psod内部设计选取引物为OPs 1，然后分别以野生型C. glutamicum和重组菌株CZM03的基因组DNA为模板，用OPs 1和CPyr-D 2 (OPs 2)进行扩增检测，结果见图2(b)，泳道3为重组菌株CZM03的基因组DNA为模板扩增结果，有特异条带扩增出来，泳道4为野生型C. glutamicum的基因组DNA为模板扩增结果，无特异条带出现，这进一步证实Psod已成功插入到pyc基因上游。 图2. 重组菌株的鉴定 3.3. 重组菌株的生长特性 野生型C. glutamicum (WT)、重组菌株CZM03的生长情况，结果见图3。WT、CZM03在接种后5 h左右进入对数生长期，对数生长约15 h后，进入平台期，平台期维持约10 h左右，进入衰亡期。两株菌的生长速率和平台期达到的最大生物量，没有明显差异。这表明Psod强启动子的敲入，对菌的生长没有明显影响，但是PYC的酶活有显著提高(见图4)，这对于目的产物的积累是有利的。 图3. 不同菌株的生长曲线 3.4. 重组菌株PYC酶活分析 丙酮酸至苹果酸的厌氧代谢通路为：丙酮酸→草酰乙酸→苹果酸。丙酮酸至草酰乙酸的转化步骤由PYC催化，通过Psod增强PYC酶活，可以加速丙酮酸向草酰乙酸的转化，进而利于苹果酸的生成。为此，我们检测了CZM03菌株的PYC酶活性，实验结果表明，敲入Psod后，C03菌株的PYC酶活较野生型上调了约2.2倍。 图4. PYC酶活测定和代谢产物分析
system: [SYSTEM]

user: 4. 结论 1) 构建强启动子Psod敲入载体，成功地在野生型C. glutamicum的丙酮酸羧化酶基因(pyc)前面敲入Psod启动子。 2) 强启动子Psod敲入后，所获得的重组菌株，与野生型C. glutamicum相比，生长速度和最大生物量无明显差异。 3) 与野生型C. glutamicum相比，重组菌株CZM03的PYC活性上调了约2.2倍。
system: [SYSTEM]

user: 摘 要 本文研究了三次微分系统 { x ˙ = − y + δ x + m x y + l x 2 + y 3 y ˙ = x ( 1 + a x + b y ) 的异宿环分支极限环问题。作者利用Melnikov函数计算未扰系统的异宿轨经扰动破裂以后的稳定流形和不稳定流形之间的相对距离，结合奇点的类型与稳定性分析，给出了系统存在稳定极限环和不稳定极限环的条件。 关键词 :三次多项式系统，奇点，异宿环，分支，极限环 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. 引言 自D. Hilbert 16问题提出的一百多年以来，关于二维多项式微分系统的研究一直受到常微分方程专家学者们的高度重视并一直以来都是研究的热点，同时也获得了丰富的成果(见 [ 1 ] [ 2 ] )。但遗憾的是，该问题的研究距离完全解决还有非常遥远的距离，其中，即使对于二次、三次这样的低次幂多项式系统，至今也未能彻底解决。传统的方法采取通过做适当的变换将原方程化为一种特殊形式(liénard形式)之后再利用定性的方法讨论极限环的存在性。近年来，已有学者采用分支的方法研究了二次多项式微分系统的奇异环分支极限环以获得极限环的存在性条件(见 [ 3 ] [ 4 ] [ 5 ] [ 6 ] )， [ 7 ] 研究了Duffing方程的奇异环分支极限环问题， [ 8 ] 研究了软弹簧型方程在摄动下分支出极限环问题。本文作者将原方程化成一种特殊形式，再采用分支的方法对系统进行分析，来研究三次多项式系统的异宿环分支极限环的存在性问题，给出产生极限环的存在性条件和稳定性分析。
system: [SYSTEM]

user: 2. 预备知识 考虑平面自治系统 { x ˙ = X ( x ， y ) y ˙ = Y ( x ， y ) (2.1) 其中 X ， Y ∈ C 1 ， x ， y ∈ R 1 。 2.1. 环域定理 [ 9 ] [ 10 ] 设由两条不相交的简单闭曲线 L 1 与 L 2 ( L 2 ⊃ L 1 ) 所围成的环形域D内及其边界上不含系统(2.1)的奇点， L 1 与 L 2 均不是整条闭轨线，若平面自治系统(2.1)的凡与D的边界 ∂ D 相交的正半轨线均进入(均离开)环域D，则在D内至少存在此系统的一个内稳定(不稳定)的极限环和一个外稳定(不稳定)的极限环，二者可能重合为一。 L 1 与 L 2 称为环域的内外境界。 注1：环域定理的条件还可减弱，只要进入D的正(负)半轨不能跑出 D ¯ 也不能以 ∂ D 为 Ω ( A ) 极限集，则结论照样成立。这种减弱意味着：1) ∂ D 上可以有此系统的轨线段，但D的内或外境界线不能全由轨线构成；2) ∂ D 上可以出现有限个奇点，只要保证轨线一旦进入(离开) D后不再离开(进入)即可。 注2：环域D的内境界线可以缩小成一个不稳定(稳定)的奇点。 注3：若D内仅存在唯一的闭轨线 Γ ，则 Γ 必是一稳定(不稳定)的极限环；若D内仅存在有限个闭轨，则逐个分析这些闭轨线的稳定性易知，D内至少有一个稳定(不稳定)的极限环；D内还有可能存在周期或复合极限环情形，但对解析系统这种情况不可能发生。 2.2. Melnikov函数 [ 11 ] [ 12 ] 考虑平面自治系统 { x ˙ = f ( x ， y ) y ˙ = g ( x ， y ) (2.2) 及其扰动系统 { x ˙ = f ( x ， y ) + p f 0 ( x ， y ， p ， q ) y ˙ = g ( x ， y ) + p g 0 ( x ， y ， p ， q ) (2.3) 其中 f ， g ∈ C 1 ， x ， y ∈ R 1 。 f 0 ， g 0 ∈ C 1 ， p ∈ R 1 ， q ∈ R k ， k ≥ 0 。假设 1) 系统(2.2)存在异宿于鞍点 O 1 和 的两条异宿轨 Γ 1 和 Γ 2 构成的异宿环 Γ ， P 0 为 Γ i ， i = 1 ， 2 上任意一点，过 P 0 作(2.2)的横截线l与 Γ i 在 P 0 点的外法线方向 n 共线。 2) 扰动系统(2.3)在 O 1 和 O 2 点附近的鞍点分别为 O ¯ 1 和 O ¯ 2 ，过 O ¯ i ， i = 1 ， 2 的稳定流形 W O ¯ i s 和不稳定流形 W O ¯ i u 与l的交点分别为 P i s 和 P i u 。 则在小扰动下，从 P i s 到 P i − 1 u 的有向距离 d ( P i s ， P i − 1 u ) ( P i s P i − 1 u 与 n 同向时为正)为： d ( P i s ， P i − 1 u ) = p M f 2 ( P 0 ) + g 2 ( P 0 ) + o ( p ) (2.4) 其中， M = ∫ − ∞ + ∞ e − ∫ 0 t ( f x + g y ) d τ ( f g 0 − g f 0 ) p = 0 d t 为Melnikov函数。
system: [SYSTEM]

user: 摘 要 作为中国西南大省，四川省长期受到干旱灾害的影响。本文选用四川省37个气象站1977~2017年的逐日降水数据，四川省21市州的社会经济和地理信息数据，基于自然灾害风险理论，从致灾因子危险性、成灾环境敏感性、承灾体易损性和防灾减灾能力四个方面分别选取了10个评价指标，通过计算标准化降水指数SPI，利用层次分析和加权综合评价法，以ArcGIS 10.2为工具构建了四川省干旱灾害风险评估模型，并对四川省干旱灾害风险进行综合分析和区划。结果表明：四川中部及东部地区的成都、德阳、遂宁、资阳、内江、眉山、自贡、宜宾、广安、绵阳、乐山、南充市的旱灾风险偏高；四川西部地区、东南部的泸州和东北部的广元、巴中市旱灾风险偏低，多为低和较低风险区。构建的模型总体上可以反映四川省旱灾综合风险的分布情况，有利于制定和完善防灾减灾措施，从而主动规避旱灾风险，减少旱灾损失。 关键词 :干旱灾害，SPI指数，风险分析，风险区划，四川省 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. 引言 近年来，在全球气候变暖的背景条件下，干旱灾害日益严重，其发生的频率不断升高，影响范围也逐渐扩大。据世界气象组织(WMO)统计，在所有的自然灾害中，气象灾害的比例达70%，其中干旱灾害约占气象灾害的半数之多 [ 1 ]。四川省地理环境分布复杂、水土资源分配不均、人口增长迅速、社会总需水量增加，构成了形成干旱灾害的基本背景，根据四川省统计年鉴中历史旱灾灾情资料记载，近年来四川省常年干旱缺水的农田占全省耕地面积的60%以上。1952~2017年间，四川省旱灾的受灾面积为4664.1万ha，占受灾总面积的55.42%，为受灾率最高的气象灾害。 干旱灾害风险的定义为自然环境与社会经济受干旱灾害影响和危害的可能性 [ 2 ]。由自然灾害风险理论可知，干旱灾害风险(D)是致灾因子危险性(R)、成灾环境敏感性(S)、承灾体易损性(V)和防灾减灾能力(C)四个因子综合作用的结果 [ 3 ]。当前国内外许多学者都对自然灾害风险评估进行了相关研究。如，Philip等人利用气象、农业和社会经济数据绘制了加纳地区受干旱灾害影响的分布图 [ 4 ]。Shahid在GIS环境下，选取了7个自然及社会类评价指标，对孟加拉国西部地区的干旱灾害风险及综合干旱脆弱性进行了风险评估 [ 5 ]。Hungsoo提出了以降水为基础的干旱发生概率和干旱脆弱性指数(DVI)、干旱危害指数(DHI)，对韩国的干旱灾害进行了风险评估 [ 6 ]。姚玉璧等通过建立评估指标，对中国南方干旱灾害风险特征进行分析，并提出旱灾防御对策 [ 7 ]。徐玉霞等人基于GIS，从致灾因子危险性等四个方面对陕西省的干旱灾害风险进行了评估与区划分析 [ 8 ]。何娇婻等将自然灾害风险理论与云南省实际情况相结合，进行了云南省的干旱灾害风险评估和区划 [ 9 ]。王莺等依据相关的气象和社会数据等多方面因素，在GIS平台下做出了甘肃省河东地区的干旱灾害风险分析 [ 10 ]。 从研究现状来看，对于干旱灾害风险，国内外目前都尚未形成较为标准的评估指标体系，因此在研究时，需要结合研究区域的自然环境特征、社会经济特征以及旱灾发生规律，建立符合当地的评价指标模型。本文基于自然灾害风险理论，采用气象数据、社会经济数据和地理信息数据，构建针对四川省的干旱灾害风险评估模型，对四川省干旱灾害风险进行分析及区划，从而为旱灾风险管理提供科学依据，主动采取防灾减灾措施以规避旱灾风险，减少经济损失。