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

user: 5. 适用方法 若要将测谎意见作为证据种类采用，那么其在证明方面的适用也必然要接受法律的规制。而确定具体适用的原则与规则、制定完善的机制不仅可以保证测谎过程的有序进行和测谎意见的正确适用，也可以保护被测试人的合法权利。鉴于测谎技术本身的特殊性和证据适用的原则，归纳测谎意见的适用原则如下： 5.1. 自愿原则 由于强制测谎、直接窥探他人内心的意思不仅会对被测人的内心产生压力，使得测谎结果不真实，极易出现测谎结果假阳性 10 的现象，导致测谎意见的适用反而出现消极作用；也在很大程度上侵犯了被测人的人格自由，属于非法取证。因此，在测谎之前，测试者应当与被测人达成合意之后签订代表被测人同意的测谎协议之后，再进行测谎。才能达到将测谎意见列为法定证据的立法原意，使得测谎技术实现其真正价值。 5.2. 非法测谎排除原则 非法证据排除原则旨在维护程序正义、保障被害人的合法权益，在我国刑事诉讼法中也有明确规定。该原则对所有法庭证据的可采性均有影响，故其理所应当对测谎意见也具有强制限制。 5.3. 严格限制原则 由于测谎技术是通过将人的脉搏、心脏电位等内征变化通过具体数值外化，以判断其所述真伪。对数值变化的规律采取经验法则，即适用测试员所掌握的一般规律做出判断。但诸多因素也可对数值变化的结果造成影响，比如被测人本身的身体状况(例有无心脏疾病)、被测人接受测谎的次数、被测人在接受测试时的心理状态等等。故在适用对象和适用范围上的限制，需要制定出具体的实施细则后结合个案进行分析和考量。 5.4. 有限采用规则 证据的“有限采用规则”(limited admissibility rule)源于英美法系国家的规定，是指可以为某个限定的目的而被采纳的证据，证据的适用以该限定目的的存在为前提 [ 26 ]。测谎技术作为侦查中的特殊手段，测谎意见为特定情况下为了辅证有关言词证据、查明案件事实的鉴定意见，其应属于“有限采用”的证据。其适用限制可以体现在以下几点：其一，其可以用来审查言词证据的真实性，但不可直接论证被测试人有罪之事实。其对言词证据的作用正如笔迹鉴定对书证之作用，特别是用来证明被测试人陈述是否真实可靠，以解决伪证问题，但不能直接用来证明案件事实，进而确定被测试人的法律责任；其二，考虑到其属于情态证据，具有“论心定罪”之危险，故其在适用时应当遵循刑事诉讼中有利被告人原则，即测谎意见所指向的诚实论述可以作为无辜者的认定，而审判者认为的说谎论述不得作为重刑证据，亦不可作为直接定罪之理由。 5.5. 对鉴定人员之专业素质应当严格把关 目前如美国、日本等将测谎意见采纳为法庭证据的国家均对测谎人员的素质以及能力有着硬性要求，具体要求上文已有涉及。而我国目前测谎行业并未有相应的机构以及具体的实施规则，测谎人员的能力素质也良莠不齐，这也是我国测谎意见信任度不高的原因之一。测谎意见要在我国提高支持率与可信度，测谎行业必须进行行业的内部建设以及相应的规制，对能上业的测谎人员也必有素质要求，才能保证测谎意见的正确率与可信度。 5.6. 对测谎意见有异议的救济措施应当参照我国鉴定意见的有关规定 世界上不存在正确率为100%的鉴定意见，故须有对鉴定意见出错后救济渠道之规定。上文已经论述，测谎意见的证据属性为鉴定意见，故其在证明力存在瑕疵时的救济途径也适用于我国法律中关于鉴定意见的有关规定。
system: [SYSTEM]

user: 6. 结语 测谎证据在特定条件下基本符合证据的三个属性特征，也具有应然层面上的证据能力与证明力。故可以将其归为鉴定结论的证据种类，在证据理论上属于间接证据范畴。测谎技术有利于案件的侦办，对待其正确的态度应为在刑事诉讼中事实认定确有困难时辅以测谎证据来认定事实。鉴于测谎证据的应用实践先于理论与立法，适时推动测谎证据立法化、制度化不仅有利于解决司法实践中的问题，还能对测谎意见的适用进行正确的规制，更好地惩治犯罪、保障人权。
system: [SYSTEM]

user: 摘 要 目的：探讨小切口锁定钢板内固定治疗肩胛骨骨折的手术效果。方法：选取2010年3月至2018年5月重庆市第五人民医院收治的30例肩胛骨骨折的患者，采用微创切口锁定钢板内固定。使用Hardegger评分系统评价疗效。结果：术后门诊随访6~12个月，平均4.5个月骨折愈合。评分结果显示：功能优：26例，良：4例。结论：小切口锁定钢板固定肩胛骨骨折创伤小，切口感染率低，术后能早期功能锻炼，肩关节功能恢复好，疗效满意。 关键词 :肩胛骨，小切口，锁定钢板，骨折 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 肩胛骨为三角形扁骨，表面覆盖有冈上肌、冈下肌、大圆肌、小圆肌、斜方肌等，肩胛骨的运动可上提，下抑，外旋，内旋，外展及内收等，具有重要的功能。随着国民经济的发展，在生产生活中，意外事故的增多导致肩胛骨骨折的发病率逐渐上升，肩胛骨骨折多由于直接暴力如高处坠落伤、重物砸伤等所致，多为粉碎性骨折，既往对于肩胛骨体部骨折多采取保守治疗，保守治疗采用三角巾或者绷带悬吊患肢制动，保守治疗时间长，患者肩关节功能恢复差。肩胛骨体部及颈部的骨折的手术入路可采取经典的后方入路如Judet入路、改良Judet入路、后方直切口入路等，使用钢板内固定。但存在创伤大，出血多，术后易出现手术切口愈合不良，影响治疗效果。我院自2010年3月至2018年8月采用后方微创切口治疗肩胛骨体部骨折，治疗效果比较满意，创伤小，现报道如下：
system: [SYSTEM]

user: 2. 资料与方法 该研究遵从世界医学协会赫尔辛基宣言，获得相应的伦理许可，所有患者均签署知情同意书。 选取2010年3月至2018年5月重庆市第五人民医院收治的30例肩胛骨体部骨折的患者，采用后路微创切口重建锁定钢板内固定手术治疗。其中：男性26例，女性4例。年龄：20~50岁，平均40.6岁。左侧肩胛骨骨折13例，右侧肩胛骨骨折14例，双侧3例。重物砸伤5例，高处坠落伤15例，车祸伤7例，电击伤3例。纳入标准：1) 骨折时间小于3周；2) 随访时间大于3个月；3) 所有患者均完善肩关节正侧位片、肩胛骨CT + 三维重建，患者均知情理解，要求手术。排除标准：1) 肩胛骨病理性骨折、老年骨质疏松患者及儿童患者；2) 合并严重系统性疾病不能耐受全麻及手术者；3) 要求保守治疗的患者；4) 合并颅脑、胸腹等全身重要器官损伤的患者及臂丛神经损伤的患者。手术日均在伤后6 d以后，避开水肿高发期，手术时皮肤肿胀已明显消退，出现皮肤皱褶。 手术方法：术前注射氨甲环酸以减少出血，使用气管插管进行全身麻醉，受治疗者，俯卧位，患肢与身体成90度角，进行消毒，并根据骨折部位的具体情况进行切口。背部偏外侧的切口与肩胛骨外侧边缘保持平行，沿身体纵轴方向切开皮肤、皮下的组织、浅筋膜、深筋膜至肌肉层，切开三角肌的一部分并将其向外拉，向外侧拉开可暴露其下的冈下肌和小圆肌，从冈下肌与小圆肌间隙进入即可显露外侧骨折部位，清理骨折末端的血块，并把肩胛骨外侧缘的骨折块手法复位后，为了牢固的固定骨折断端，将钢板放在肩胛骨的外侧边缘，成型后，使用螺丝钉固定。如合并肩胛岗处的骨折，可向内侧牵开暴露肩胛岗，必要时可辅助增加肩胛岗处皮肤的横切口暴露肩胛岗的骨折。如合并内侧脊柱边缘的骨折，沿肩胛骨内侧脊柱边缘切开皮肤，皮下，显露肩胛骨内侧脊柱边缘，把肩胛骨内侧缘的骨折块手动复位后，将钢板沿肩胛骨的内侧脊柱沿放置，固定螺丝钉。电外科止血，用大量盐溶液冲洗伤口后，在伤口上安装真空引流软管以保持引流畅通。将伤口逐层缝合，并应用无菌敷料覆盖。手术后，可以静脉注射氨甲环酸以止血。术后处理：术后常规伤口换药，第二天拔除伤口引流管，术后第一天即开始肘关节、腕关节、掌指关节的屈伸，以减少手臂肿胀和保持关节运动。疼痛缓解后可练习外展肩关节，练习摆动运动，练习被动抬高患肢，外旋、内旋、外展等。使用可调节外展支具固定，伤口拆线后，术后2到3周加强上肢被动功能锻炼，同时复查X片，发现骨折是否移位。术后3到6周继续加强患肢主动功能锻炼，术后1、3、6和12个月定期随访。 结果：30例患者全部获得骨性愈合，肩关节功能好，平均手术时间45~120 min，出血量100 ml，术后无切口感染发生，无静脉血栓形成。术后随访情况，术后X片提示复位满意，术后平均骨折愈合时间：3.5月，术后肩关节评分采用Hardegger评分 [ 1 ]，优(肩关节前屈、后伸、上举、外展低于健侧10~15度，活动范围无明显受限)：26例，良(肩关节外展、上举及内外侧旋转低于健侧16~30度)：4例。如图1~4所示。 图1. Judet入路手术切口 图2. 术中使用钢板 图3. 后方微创入路手术切口 图4. 术中使用钢板
system: [SYSTEM]

user: 3. 讨论 肩胛骨浮于胸廓后面，是人体最大的扁骨，分为两面、三缘及三个角，肩胛岗延伸为肩峰 [ 2 ]。肩胛骨的位置较深，肩胛骨的腹侧面有肩胛下窝，背侧面由肩胛岗分为冈上窝及冈下窝，肩胛骨表面覆盖肌肉众多，包括三角肌、大圆肌、小圆肌、冈上肌、冈下肌、肩胛舌骨肌、肩胛下肌等，肩胛骨上缘的肩胛上切迹有肩胛上神经及肩胛上血管通过 [ 3 ]，支配肩胛骨上的肌肉，对肩关节的外展、外旋具有重要的功能。肩胛骨的血管主要由肩胛上动脉、旋肩胛动脉及肩胛背动脉组成，肩胛上动脉由肩胛上切迹进入岗上窝，旋肩胛动脉由三边孔进入肩胛骨后侧，肩胛背动脉临近肩胛骨内侧缘。肩胛骨骨折的受伤原因多为严重暴力所致，多合并肋骨骨折、颅脑外伤、脊柱骨折等，属于严重复合伤的一部分，合并伤易危及生命 [ 4 ] [ 5 ] [ 6 ]。肩胛骨骨折的治疗主要是在保障生命安全的基础上进行的，CT及三维成像检查在肩胛骨骨折的诊断上有重要的作用 [ 7 ]，可以指导肩胛骨骨折的分型及治疗。肩胛骨骨折分型主要有根据肩胛骨形态特点分型的Miller分型及根据肩胛骨解剖部位分型的AO分型以及根据肩胛骨骨折部位分类的Hardegger分型 [ 8 ]，肩胛盂骨折常采用Ideberge分型 [ 9 ]。其中肩胛骨体部骨折在肩胛骨骨折中较常见，肩胛骨骨折后，患者常见的临床表现为患侧肩胛部疼痛剧烈，患侧肩部不能抬举，可合并其他损伤，如不及时治疗，容易遗留肩关节功能障碍，对患者的生活质量有严重的影响。对肩胛骨体部骨折，移位不明显者可采取保守治疗，保守治疗采用三角巾悬吊患肢制动，疼痛缓解后进行肩关节的功能锻炼，无明显移位的肩胛骨体部骨折采用保守治疗效果较好，能够恢复较好的肩关节功能，但肩胛骨骨折多为粉碎性骨折，对于移位明显的骨折，保守治疗不能有效的恢复盂极角，不能恢复肩胛骨的解剖结构，导致骨折畸形愈合，同时患者疼痛缓解不明显，固定时间长，康复时间长，不能早期功能锻炼 [ 10 ]，可导致肩胛骨周围肌肉粘连，遗留肩关节的功能障碍及慢性疼痛，降低生活质量，影响青壮年的劳动力，故目前对肩胛骨粉碎性骨折多采取手术治疗。对累及肩胛骨体部及颈部骨折的手术治疗既往常采用后侧入路，如Judet入路 [ 11 ]、改良Judet入路及外侧缘直切口入路。Judet入路采用倒L切口，切口起于肩胛骨内侧缘及肩胛岗，剥离皮瓣，切口长，暴露范围广，出血多，容易破坏肌肉的血供，导致冈下肌等肌肉坏死 [ 12 ]。改良Judet入路皮肤切口与Judet入路皮肤切口相似，从冈下肌及小圆肌的间隙进入，对肌肉的损伤较小，但存在皮下剥离范围大，如果术后引流不充分或皮下出血，易导致皮下积液及皮瓣坏死，同时于肩胛骨颈部骨折，对外侧缘的暴露不够充分，过度牵拉皮瓣及肌肉可能损伤血管及神经 [ 13 ]。单纯使用外侧缘直切口治疗肩胛骨体部及颈部骨折，从冈下肌及小圆肌的间隙进入，可暴露肩胛骨的外侧缘及颈部，但对肩胛骨的内侧缘暴露不佳，如合并肩胛骨体部内侧缘及肩胛岗内侧部分的骨折，单纯使用外侧缘直切口对骨折断端的暴露不够充分，不利于骨折的复位。采用后方微创切口治疗肩胛骨体部及颈部骨折，沿肩胛骨外侧缘做切口，从冈下肌和小圆肌之间暴露骨折断端，内侧缘切口沿肩胛骨内侧缘进入，肩胛骨的内外侧缘及肩胛岗处骨质较厚，螺钉固定较稳固，可对肩胛骨骨折内侧缘及外侧缘进行固定 [ 14 ]，不需要显露肩胛骨体部如冈上窝及冈下窝等处的多个骨折线，创伤小，便于复位。累及肩胛骨内外侧缘及肩胛岗处的严重的粉碎性骨折对骨折复位后可使用重建钢板进行三角形框架固定 [ 15 ]，钢板根据肩胛骨的生理弧度进行塑形后比较贴附，有研究表明，对累及肩胛骨内外侧缘及肩胛岗处的严重粉碎性骨折，单纯固定肩胛骨骨折部位的内侧缘或外侧缘，不能达到良好的复位及固定。采用后方微创切口治疗肩胛骨体部及颈部骨折具有创伤小，出血少，便于复位及固定骨折，术后具有患者疼痛缓解明显，可以早期功能锻炼的优点 [ 16 ]，是一种较好的手术方法，值得推广。
system: [SYSTEM]

user: 摘 要 森林是地球之“肺”，具有保持水土，调节气候，防止水、旱、风、沙等灾害的作用。文章利用绥宁1963年~2017年森林火灾资料和全市2001~2017年的火灾资料，以绥宁作为森林火险研究的主要对象，结合气象因子对森林火险的影响，根据国家和省的标准，确定森林火灾等级，根据气象要素来确定森林火险天气指数，结果表明，1) 森林火灾具有明显的季节性，春季出现的频率最高，达60.1%，冬季次之，占21.8%，夏季最少，只有4.2%。2) 森林火灾可划分为4个等级，即一般、重大、比较重大、特别重大。3) 森林火险天气指数划分为6级。4) 森林火险等级可划分为5级，当森林火险等级达到4级或以上时，就要发布火险预警信息，做好防范工作。 关键词 :森林火险等级，预报技术，火灾特点，气象因子，火险天气指数 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 森林是地球之“肺”，它具有保持水土，调节气候，防止水、旱、风、沙等灾害的作用 [ 1 ]。 近年来由于森林火灾的不断发生，给国家经济、人民财产和生命安全带来很大危害。森林火灾的危害相当严重，它不仅会烧毁林木，降低木材资源的数量和质量，造成大量动植物死亡，破坏大自然的生态平衡，而且还会引起森林小气候的变化，使林区局部地区温差增大，水分蒸发增加，温度降低，无霜期缩短。森林火灾发生的次数和危害程度以春季最多，最严重，秋季次之。由于每年气象条件变化的明显差异，火险期的早晚、长短也会出现不同的情况，例如冬季降雪量偏少，春季气温回升快，大风日数多，火险期会相应提前；若雨季开始晚，春旱持续时间长，火险期则会相应延长。如果降水量比常年偏多，冷空气活动频繁，气温降低，火险期则会相应推迟或缩短。防火的重要意义就在于此，森林火灾是破坏森林资源安全，威胁人类环境最为严重的灾害之一，森林火灾的发生与否，发生后的蔓延和扩大，除与林木种类、疏密程度、树木年龄等条件有关外，还与气象条件有着密切的关系。气象条件不仅综合影响森林的可燃性，而且影响火的蔓延和林火行为的特点。因此，掌握气象条件与森林着火和火灾蔓延的关系，做好森林火灾发生可能性的预报也就是森林火险预报，对于及时扑火和消灭森林火灾，防止造成灾害，具有重要意义。一直以来，森林火险没有一个定量的预报，而定性预报往往存在着不确定性。因此，制定定量的森林火险等级预报，是完全必要的。我们根据《中华人民共和国林业行业标准》中，全国森林火险天气等级的规定来制定绥宁森林火险等级预报，提高我们对森林火灾的防御能力，保护人民的生命财产安全，减少森林火灾所造成的损失，是防灾减灾的迫切要求。对社会经济发展、社会和谐、人民生命财产安全等方面都会产生深远的政治影响。 在全球性持续变暖的气候背景下，一些地区林火发生频率有显上升趋势。据统计我国平均每年发生森林火灾1.43万次，受害森林面积8.212 ×10 5 hm 2 ，占世界火灾次数的14%，平均森林受灾面积为世界的20%。其中大、小兴安岭平均每年发生森林火灾近30次，为我国森林火灾高发区之一，必需做好林火监测预报工作。把卫星遥感技术应用于森林火灾后植被恢复情况的监测，早在上世纪70年代美国就开始利用地形、植被和气象因素的差异而导致森林火灾中林地火烧强度的不同，而用卫星红外遥感技术取得的数据来反映火烧程度的差别，并可得到Landsat影像图，植被和火烧程度的变化都可由卫星数据表现出来。目前欧美一些国家普遍应用雷击火自动定位系统、遥感自动林火天气站网，以及GIS和GPS与林火管理相结合等。芬兰用NOAA、AVHRR数据而设的自动化林火管理工作站已投入使用；1992年波兰发生林火，用法国的SPOT XS和XP图象以及ERS-1，及时地监视了林火燃烧状况。林火数据库(FFD)是提供信息的基础；波兰的FFD储存32种关于林火燃烧的信息：气象要素、燃烧类型、面积数、林火特性及火灾控制信息等，加上GIS对分析林火危险程度非常有用 [ 2 ] [ 3 ]。国内，林火监测大多利用卫星遥感、航空、望塔等手段，闪电定位系统才刚刚起步，至于GPS、GIS系统的结合处于初期试验研究阶段，与先进国家相比差距较大。1987年大兴安岭地区发生森林大火后，国家卫星气象中心、黑龙江省气象台首先使用极轨卫星NOAAPAVHRR对林火进行监测，至今黑龙江省森林大火无一漏报。
system: [SYSTEM]

user: 2. 资料来源与方法 本文所用资料为绥宁县林业局提供的1963年~2017年森林火灾资料和邵阳市林业公安局提供的邵阳市2001~2017年的火灾资料(期间因种种原因资料有部分中断)，以及1963~2017年绥宁气象局的地面气象观测资料，利用统计分析方法，统计出林火面积100亩以上的情况，找出火灾与气象因子的关系，划分出森林火险等级，用于日常森林火险预报，发布森林火险预警。
system: [SYSTEM]

user: 3. 邵阳森林火灾发生特点 邵阳市以绥宁森林面积为全市之最，因此，我们以绥宁作为森林火险研究的主要对象，结合全市2001~2017年的火灾情况，提炼森林火险等级指标。据不完全资料统计，1963~2002年、2008年~2011年根据1973~1975年、1978年、1981年、1983~1988年、2008~2017年较为详细的资料统计，森林火灾具有明显的季节性，春季出现的频率最高，达60.1%，冬季次之，占21.8%，夏季最少，只有4.2%，见表1。林火出现频率的月分布规律，以4月最多，占24.%，3月次之，占21.1%，6月、8月和9月最少占1.2%，见表2。历年最大森林火灾面积出现在3月，1月次之，见表3。 Table 1 Frequency table of forest fire in each seaso 名称 春季(3~5月) 夏季(6~8月) 秋季(9~11月) 林火(含火警) 60.1 4.4 13.7 ≥100亩的 46.5 4.7 16.3 ≥500亩的 36.4 0 15.1 ≥1000亩的 27.6 0 11.2 表1. 各季林火出现频率 Table 2 Forest fire frequency in each mont 名称 1 2 3 4 5 6 7 8 9 10 11 12 林火(含火警) 6.7 8.3 21.1 24.9 14.0 1.2 2.0 1.2 1.2 4.9 7.6 6.9 ≥100亩的 9.3 11.6 33.7 9.3 3.5 0 3.5 1.2 0 7.0 9.3 11.6 ≥500亩的 12.1 18.2 30.3 6.1 0 0 0 0 0 3.0 12.1 18.2 ≥1000亩的 16.7 16.7 27.6 0 0 0 0 0 0 5.6 5.6 27.8 表2. 各月林火出现频率 Table 3 Area affected by the largest forest fire in each month of the calendar year (Unit: hm2 月份 1 2 3 4 5 6 7 8 9 10 11 12 面积 400.2 166.75 433.55 62.90 10.67 10.01 13.34 6.67 4.67 106.72 106.72 500.25 表3. 历年各月最大一次林火受害面积(单位：公顷)
system: [SYSTEM]

user: 4. 森林火险与气象因子关系 大量研究表明，森林火险与气象条件是非常密切相关的。为进一步研究森林火险与气象条件的密切关系，我们根据日常业务情况和预报的可操作性，对几个主要气象因子与森林火灾的关系进行研究，即湿度、降水量、温度、风、雷暴。 4.1. 空气相对湿度与森林火灾的关系 空气相对湿度能影响森林火灾的发生一般情况下，当相对湿度大于70%时，不易发生火灾；但如果长期未下雨，即使短期内相对湿度达80%，也可能发生火灾；相对湿度在50%~70%时，可能发生火灾；小于50%时，可发生大火灾，若小于35%时，发生特大火灾的可能性就更大。根据班玛县气象局多年观测的气象数据，县区内年平均空气湿度为60%，林区内湿度大于60%，但在冬春两季空气干燥、湿度较小森林火险等级较高，冬春季节是林区防火重要时期。故加强两季巡山与护林防火宣传是重点工。 4.2. 森林火灾发生的次数与降水量和降水日数的关系 森林火灾发生的次数与降水量和降水日数有很好的对应关系降水量减少，无雨日较长，森林可燃物的含水量将不断下降，森林火灾发生的可能性和严重性也随之增大。一般大雨过后，3~4 d之内不会出现火险。在降水偏少的月份发生的次数也较多；连续无降水达3 d以上，发生森林火灾的可能性较大，日降水量在2~5 mm时，因为有了降水的滋润，降低了林区的燃烧性，仍有可能发生火灾，但机率不是很高；大于5 mm时，林区地面可燃物和树枝等吸收了较多的水分，有的甚至达到饱和状态，杜绝了星星之火的后患。 4.3. 蒸发量和气温日较差与森林火灾也有密切的关系 气温较高，日较差7℃~13℃，日蒸发量3~5 mm，能引发火灾；气温很高，日较差达13℃以上，日蒸发量超过5 mm，极易发生火灾。因此，气温越高，日较差越大，蒸发作用越强，越容易引发森林火灾。空气温度对森林火灾的影响是多方面的。温度越高，可燃物中水分蒸发和变干的速度越快，火灾发生的可能性越大。气温影响可燃物的着燃性，高温还会促使火势更加猛烈。所以高温干燥的天气情况下更应加强防范，在有必要和有条件的情况下最好和气象部门取得联系进行人工降雨或飞机洒水等措施增加空气湿度减少火灾发生的几率。 4.4. 风与森林火灾的关系 风是森林火灾扩大蔓延的重要因素风力的大小与森林火灾的发生关系尤为密切。风不仅能把植被吹干，有助于燃烧，同时还使小火扩大，并能使死灰复燃。而且在火灾发生后，还能使火源得到充分的氧气供应，加速燃烧，同时使火花飞溅，影响火灾的形状，延伸火灾的外形，扩大火灾面积，使地面火变为树冠火。看看各起火灾的规模与风的关系，则可知风越强越易发生大火，大规模的火灾也都是因风蔓延的。所以及时了解火灾发生的风向与风速是灭火的重要因素。 4.5. 雷暴与森林火灾的关系 雷暴与闪电可带来森林火源。众所周知，大自然给人类带来了火，人类利用火走向了文明。但是地球上大部分的森林火灾因闪电而引发的，所以建立森林地区雷电预警系统，并且应有一定的防雷措施。可以在气象部门的配合下在林区内建造独立避雷塔，接地装置的接地电阻值为4 欧姆。这在一定程度上降低了因雷电引发的火灾。
system: [SYSTEM]

user: 5. 森林火险等级划分及其预报方法 5.1. 森林火灾等级划分标准 我国1995年12月1日颁布实施了《全国森林火险天气等级》标准 [ 4 ]，1992年12月1日颁布实施了《全国森林火险区划等级》，2008年进行重新修订 [ 5 ]，1988年经国务院制定发布《森林防火条例》，2008年11月19日经国务院第36次常务会议重新修订通过，根据新修订后的《森林防火条例》进行森林火灾划分，确定等级。 按照国家和省市标准，森林火灾划分为四种类型，见表4： Table 4 Classification of forest fire type 火灾类型 受害森林面积(公顷) 伤亡人数 一般 1公顷以下 死亡1人以上3人以下或重伤1人以上10人以下 较大 1公顷亩以上~100公顷以下 死亡3人以上10人以下或重伤10人以上50人以下 重大 100公顷上~1000公顷以下 死亡10人以上30人以下或重伤50人以上100人以下 特别重大 1000公顷以上 死亡30人以上或重伤100人以上 表4. 森林火灾类型划分 注：本表所称“以上”包括本数，“以下”不包括本数。 5.2. 森林火险等级划分及其预报方法 为了保护森林资源，掌握防火的主动权，有关部门根据气象条件变化的季节特点，规定了防火期，一般分春季防火期和秋季防火期。许多地方还依据气象与森林火灾的关系，划分了森林火灾等级，并且在发布未来天气预报时，也发布森林防火预报。我们根据季节和各月气象要素的长期预报，估计相应季节和月份内林火发生的可能趋势，在防火期间每天分析气象要素的变化，预测未来林火发生的可能性大小，确定多个气象因子与森林火灾的相关关系，制作出森林火险天气等级预报。森林火险气象等级的划分主要是考虑天气和气象要素的变化对森林可燃物易燃性的影响。冬天森林可燃物较多，如果有火种，有合适的气象条件(湿度、风速等)就可以助火越烧越大；“天干物燥”说的就是这个意思。目前中央气象台业务上预报的森林火险气象等级有五级。如果预报的森林火险气象等级超过了(含)四级，表示预报区域“气温高，湿度小，风力大，林区干燥，有利起火，容易蔓延。属高级火险，易引起火灾且蔓延较快”。如果到了五级就比四级火险等级还要高。我们根据气象要素来确定森林火险天气指数(我们用SHZ表示)，SHZ的值分为5级，其计算公式为：SHZ = A + B + C + D + E − F，式中A、B、C、D、E、F分别如下表表示(见表5~11)。各级别森林火险及危险程度见表12。 Table 5 Forest fire risk index A value of daily maximum temperature in forest fire prevention perio 空气温度等级 最高气温(℃) 森林火险指数A 一 ≤5.0 0 二 1~10.0 4 三 10.0~15.0 8 四 15.1~20.0 12 五 20.1~25.0 16 六 ≥25.0 20 表5. 森林防火期每日最高气温的森林火险指数A值 Table 6 Relative humidity B value at 14 hours corresponding to forest fire prevention fire risk inde 14时湿度等级 14时相对湿度% 森林火险指数B 一 ≥65 0 二 55~64 4 三 45~54 8 四 35~44 12 五 25~34 16 六 ≤24 20 表6. 森林防火期火险指数对应的14时相对湿度B值 Table 7 Number of consecutive days without precipitation corresponding to C value of fire risk inde C值 0 5 10 15 20 25 30 35 40 45 50 连续无降水日/d 当天 1 2 3 4 5 6 7 8 9 >10 表7. 火险指数C值对应的连续无降水日数 注：其中日降水量为0.1 mm视为无降水日。 Table 8 Wind force at 14:00 of the day corresponding to fire risk index D valu D值 5 10 15 20 25 30 35 40 风力 ≤1 2 3 4 5 6 7 ≥8 表8. 火险指数D值对应的当日14时风力 Table 9 Forest fire risk weather index E value of daily maximum wind power rating during forest fire prevention perio 风力等级 距地面10 m高处风速(m/s) 范围中数 森林火险天气指数E 0 0.0~0.2 0 0 1 0.3~1.5 1 5 2 1.6~3.3 2 10 3 3.4~5.4 4 15 4 5.5~7.9 7 20 5 8.0~10.7 9 25 6 10.8~13.8 12 30 7 13.9~17.1 16 35 8 17.2~20.7 19 40 表9. 森林防火期每日的最大风力等级的森林火险天气指数E值 Table 1 0 Daily precipitation corresponding to F value of fire risk weather inde 雨量等级 雨量 森林火险天气指数F 一 0.1~9.9 25 二 10.0~16.9 30 三 17.0~24.9 35 四 25.0~37.9 40 五 38.0~74.9 45 六 ≥75 50 表10. 火险天气指数F值对应的当日降水量 Table 1 1 Fire risk rating corresponding to SHZ value of fire risk prediction inde 火险等级 SHZ值 火险程度 一 ≤25 无火险 二 26~50 低度火险 三 51~72 中度火险 四 73~90 高度火险 五 ≥91 极度火险 表11. 火险预报指数SHZ值对应的火险等级 Table 1 2 All levels of forest fire risk and hazardous degre 等级 危险程度 易燃性 蔓延性 森林火险天气指数 一 没有危险 不能燃烧 不能蔓延 ≤25 二 低度危险 难以燃烧 难以蔓延 26~50 三 中度危险 较易燃烧 较易蔓延 51~72 四 高度危险 容易燃烧 容易蔓延 73~90 五 极度危险 极易燃烧 极易蔓延 ≥91 表12. 各级别森林火险及危险程度 5.3. 森林火险预报的发布 根据每天的森林火险天气指数最后计算出火险等级，再做出森林火险预报服务，发布森林火险等级，当森林火险等级达到4级或以上时，就要发布火险预警信息，并及时向林业局和政府部门进行书面和口头汇报，做好防范准备，并抓住时机实施人工增雨作业，预防森林火灾的发生。
system: [SYSTEM]

user: 6. 小结 总之，从林火的监测技术、林火的预报来说，我国与先进国家相比还存在较大的差距。本文就邵阳市绥宁县的地面气象要素与森林火灾发生关系进行分析探讨，确定森林火灾等级，根据气象要素来确定森林火险天气指数，发布森林火险等级预报，但未考虑大气环流背景、林区下垫面等其他因素。今后应建立基于气象要素、利用卫星遥感技术，雷击火自动定位系统、遥感自动林火天气站网以及GIS和GPS与林火管理相结合等对林火进行实时定位和追踪，建立完善邵阳森林防火预报服务系统，提高森林火险预报水平。
system: [SYSTEM]

user: 摘 要 应用2014~2018年MODIS产品、气象观测数据、DEM等数据资料，采用SEBAL模型估算分析吐鲁番市的蒸散特征。结果表明：2014~2018年吐鲁番市平均年总蒸散量为3.25 × 10 9 m 3 ，且蒸散量时空差异显著，97%的蒸散量集中在4~10月份，年际间呈波动下降趋势；对于不同土地利用类型，单位面积年蒸散量耕地 > 城乡、工矿及居民用地 > 林地 > 草地 > 水域 > 其他用地，其中草地及其他用地所占面积较大，年总蒸散量最高，草地、耕地及其他用地三者蒸散量占总蒸散90%以上，水域多为积雪区、季节性水域且面积较小，年总蒸散量及单位面积年均蒸散量均不高。 关键词 :遥感，蒸散发，SEBAL模型，吐鲁番市 Copyright © 2020 by author(s) and Wuhan University. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 蒸散发(Evapotranspiration， ET)指区域内土壤蒸发、水面蒸发和植物蒸腾的总称，作为唯一同时出现在水量平衡和能量平衡方程中的参数，是水循环中重要的水文过程 [ 1 ]。目前随着遥感技术的飞速发展，利用多时相、多光谱、多分辨率的遥感数据定量反演蒸散发的技术逐渐成熟 [ 2 ]，多年来，国内外学者建立了诸如SEBI、SEBAL、SEBS、METRIC等遥感ET计算模型并进行了广泛应用和研究 [ 3 ] [ 4 ] [ 5 ] [ 6 ]。其中SEBAL模型是由Bastiaanssen在1998年提出的基于地表能量平衡的遥感反演模型 [ 4 ]，Wang J.、李宝富、杜嘉、周妍妍等 [ 7 ] [ 8 ] [ 9 ] [ 10 ] 在我国西北干旱区、塔里木河流域、三江平原、疏勒河等地进行了应用研究，结果显示该模型具有较好的适用性，能够很好地反映研究区地表蒸散发的实际情况。 吐鲁番市属于干旱内流区艾丁湖流域，近30年来，吐鲁番市水资源环境急剧恶化，地下水多年超采，地下水水位不断下降，坎儿井80%以上断流，艾丁湖湖面严重萎缩，已成为季节性湖泊 [ 11 ]，水资源短缺导致的植被退化、土壤盐碱化、荒漠化等水环境恶化问题突出，对此一些学者开展了吐鲁番地下水水位变化、水资源优化配置等方面的研究 [ 12 ] [ 13 ] [ 14 ]，但仍缺乏对吐鲁番水资源耗散过程的探讨。 研究吐鲁番市水资源耗散变化过程，对吐鲁番市水资源合理开发利用，促进水生态恢复，保障社会、经济和生态可持续发展具有重要意义。本文以MODIS数据产品、气象观测数据、土地利用数据、DEM数据为基础，利用SEBAL模型对吐鲁番市2014~2018年蒸散量数据进行计算验证，分析研究吐鲁番市蒸散空间分布规律，以期为吐鲁番市水资源决策管理和生态保护提供支持。
system: [SYSTEM]

user: 2. 研究区概况及数据材料 2.1. 研究区概况 吐鲁番市位于北纬41˚12'~43˚40'，东经87˚16'~91˚55'之间，属天山东部山间盆地，西临天山和觉罗塔格山，东以了墩隆起和哈密盆地分界，东西长约245千米；北部为天山东段博格达山南麓，南抵库鲁克塔格山，南北宽约75千米，总面积6.9 × 10 4 平方千米，其中平原区占77%，山地占23%。盆地整体北高南低，海拔在−155~4372 m，平均海拔1013 m (图1)。吐鲁番市属于典型的暖温带大陆性荒漠气候，降水主要集中在北部山区，盆地内全年日照充足，降雨稀少，极其干燥，多年平均降水16.4毫米，年蒸发潜力3000 mm以上，径流依靠山区降雨及融雪，河流汛期在春夏季节。人类活动集中在平原中部绿洲区域，主要农业作物是葡萄、哈密瓜、棉花、果树，由于近30年来绿洲持续扩张，水资源短缺、地下水超采、艾丁湖季节性干涸等问题凸显。 图1. 吐鲁番市DEM和气象站位置图 2.2. 数据资料 本文使用的数据包括：1) 遥感数据：2014~2018年MODIS地表温度产品(MYD11A1)、地表反照率产品(MYD09GA)、BRDF/Albedo数据(MCD43A1)，上述数据下载自NASA的LADDS数据网站(https://ladsweb.modaps.eosdis.nasa.gov)；GDEMV2 30M分辨率数字高程数据，下载自地理空间数据云(http://www.gscloud.cn/)。2) 研究区内5个国家基本气象站点数据(图1)，主要包括日最高气温、日最低气温、平均温度、气压、风速、相对湿度等，获取自吐鲁番市气象局。3) 2016吐鲁番市土地利用类型数据采用吐鲁番市世界银行贷款新疆坎儿井保护及节水灌溉项目成果，数据由Landsat 8 OLI数据解译得到，分辨率30 m。 上述所有数据资料，经在ArcGIS及ERDAS软件中格式转换、投影转换、插值、裁剪、重采样等，处理为WGS-1984坐标系、通用墨卡托投影、250 m分辨率的TIFF格式，并剔除研究区云覆盖量大于5%的影像。
system: [SYSTEM]

user: 3. SEBAL模型原理 SEBAL模型以地表能量平衡方程为基础，忽略了水平能量交换和光合作用 [ 4 ]，如式(1)。 λ E T = R n − H − G (1) 式中： λ 为汽化潜热(J∙kg −1 )，ET为蒸散量(mm)， λ E T 即为潜热通量， R n 为地表净辐射通量，H为显热通量，G为土壤热通量(单位均为W∙m −2 )。 地表净辐射通量指地表获得的太阳净辐射和大气逆辐射之和减去地表辐射能量后所得的净能量，是地表能量的主要来源，地表净辐射平衡方程见式(2)。 R n = ( 1 − α ) R s + R l ↓ − R l ↑ (2) 式中： R s 指入射至地表的太阳辐射， R l ↓ 指大气逆辐射， R l ↑ 指地表辐射(单位均为W∙m −2 )， α 指地表反照率， R s 、 R l ↓ 、 R l ↑ 计算方法如式(3)至式(5)。 R S = G s c ⋅ cos θ ⋅ d r ⋅ τ s w (3) 式中： G s c 是太阳常数，取1367 W∙m −2 ； cos θ 是太阳天顶角余弦，在平原地区 cos θ = cos ( π / 2 − φ ) ， φ 为纬度，在山区则根据坡度s、坡向 γ 、太阳赤纬 δ 、太阳时角 ω 计算(均为弧度)， cos θ = sin ( δ ) sin ( φ ) cos ( s ) − sin ( δ ) sin ( s ) cos ( γ ) + cos ( δ ) cos ( φ ) cos ( s ) cos ( ω )     + cos ( δ ) sin ( φ ) sin ( s ) cos ( γ ) cos ( ω ) + cos ( δ ) sin ( γ ) sin ( s ) sin ( ω ) ； d r 是日地距离因子， d r = 1 + 0.33 cos ( D O Y ⋅ 2 π / 365 ) ，DOY为儒略日； τ s w 是大气透射率， τ s w = 0.75 + 2 × 10 − 5 × Z ，Z为海拔高度(m)。 R l ↓   = ε a ⋅ σ ⋅ T a 4 (4) R l ↑   = ε s ⋅ σ ⋅ T s 4 (5) 式中： ε a 为大气比辐射率， ε a = 1.08 ( − ln τ s ω ) 0.265 ； τ s ω 为大气透射率，计算方式同前式； ε s 为地表比辐射率， ε s = 1.009 + 0.047 ln ( N D V I ) ，式中 N D V I > 0 ，否则为 ε s = 1 ； σ 指斯蒂芬–玻尔兹曼常数，取5.67 × 10 −8 W∙m −2 ； T a 指2 m处气温(℃)； T s 指地表温度(℃)。 土壤热通量是指土壤内部的热交换，受植被盖度影像较大。Bastaanssen等在SEBAL模型中提出利用净辐射通量、地表反照率、NDVI、地表温度的计算植被覆盖区域土壤热通量的方法 [ 15 ]，如式(6)；对于水体 G = 0.41 R n − 0.51 ；对于城镇用地 G = 0.3 R n 。 G = R n × T s / α ( 0.0038 α + 0.0074 α 2 ) ( 1 − 0.98 N D V I 4 ) (6) 显热通量指单位时间内地面和大气垂直方向间单位面积通过的热量，其产生的唯一结果是温度变化，计算方法如式(7)。 H = ρ ⋅ C p ⋅ d T / r a h (7) 式中： ρ 为空气密度(kg∙m −3 )， C p 为空气定压比热容，取1004 J∙kg −1 ∙K −1 ， d T 为 z 1 、 z 2 高度温度差值(K)， r a h 为热量传输空气动力学阻抗。 对于中性层结大气，空气动力学阻抗计算公式如下： r a h = ln ( z 2 / z 1 ) / ( k u ∗ ) (8) u ∗ = k u x / ln ( z x / z o m ) (9) z o m = exp ( a ⋅ N D V I / α + b ) (10) 式中： z 1 、 z 2 为参考植被零平面位移高度，经验值一般取0.1 m和2 m；k为卡曼常数，取0.41； u ∗ 为摩擦速度(m∙s −1 )； z x 为距地表某一高度； u x 为此高度风速； z o m 为地表动量粗糙度长度(m)。 由经验公式 z o m = 0.123 h [ 16 ]，h为植被高度，建立样本空间，据式(10)，求出常数a，b，从而计算出每个像元的 z o m 。 随着高度增加，大气受下垫面影响逐渐减弱，到达一定高度时，湍流速度不再受地表影响，此高度称为混合高度，经验值取200 m，并认为此风速作用于整个影像区域，由气象站点数据根据式(11)计算出 u 200 ，带入式(9)得出每个像元 u ∗ 。 u 200 = u ∗ k ln ( z 200 z o m ) (11) SEBAL假定 d T 和地表温度 T s 呈线性关系，即 d T = b + a T s ，并选取图像中冷热极值点，冷点为温度较低近似于气温且水分充足，显热通量为零，一般选取水面区域；热点为温度较高且极为干旱区域，ET为零，一般选取裸地或荒漠区域。根据实际冷点选取水库区域，热点选择沙漠区域，据冷热点计算出a、b系数。 a = ( R h o t − G h o t ) ⋅ r a h h o t C p ⋅ ρ a i r ⋅ ( T s h o t − T s c o l d ) (12) b = − a ⋅ T s c o l d (13) 考虑浮力效应对大气状态影响，SEBAL模型引入莫宁-奥布霍夫长度，对大气稳定状态及不稳定状态下的 r a h 、 u ∗ 按照式(14)和式(15)进行修正，并通过多次选取冷热点迭代计算a、b、dT、H，直至热点 d T h o t 、 r a h 趋于稳定。 r a h = ln ( z 2 z 1 ) − ψ h ( z 2 ) + ψ h ( z 1 ) k u ∗ (14) u ∗ = k ⋅ u 200 ln ( z 200 z o m ) − ψ m ( z 200 ) (15) 式中： ψ h ( z 1 ) 、 ψ h ( z 2 ) 是 z 1 、 z 2 高度热量传输稳定度修正因子， ψ m ( z 200 ) 为200 m处动量传输稳定度修正因子，按照Paulson [ 17 ]、Webb [ 18 ] 提出公式计算。 由于上述热通量均为瞬时状态，计算出的蒸散量是卫星影像获取时的瞬时蒸散量，SEBAL模型采用蒸发比法将瞬时蒸散量扩展为日尺度，假定每日潜热通量 λ E T 与净辐射减去土壤热通量( R n − G )的比值不变，如式(16)，则每日 E T 24 可用式(17)计算得到。 E F = λ E T / ( R n − G ) (16) E T 24 = 86400 ⋅ E F ⋅ ( R n 24 − G 24 ) / λ (17) 由于获取的部分影像数据受云覆盖影像无法使用，为计算月、年累计蒸散量，这些时段蒸散量采用气象站点蒸散量数据比值和晴好日反演的蒸散量来推算，如式(18)。 E T j = E T s t j E T s t i E T r e i (18) 式中： 为第j天日蒸散量， 、 为第i天和第j天气象站点蒸散量， 为第i天的反演蒸散量。
system: [SYSTEM]

user: 4. 精度检验 选取研究区范围内5个国家气象站的气象观测数据，以FAO 56推荐的Penman-Monteith公式 [ 19 ] 计算每个站点的地表日蒸散量作为校验依据，对SEBAL模型结果进行精度检验。 (19) 式中： 为作物系数； 为饱和水气压(kPa)； 为实际水汽压( )； 为饱和水汽压曲线的斜率； 为湿度计常数(kPa∙℃ −1 )； 为2米处风速(m∙s − 1 )。 检验结果显示，SEBAL模型计算结果的平均相对误差为11.99% (表1)，从二者相关性来看，拟合优度R 2 值为0.87 (图2)，说明蒸散量反演结果与实测值整体上有较好的正相关，认为结果精度是可以适用的。 Table 1 Relative errors calculated by SEBAL mode 站点 经度 纬度 平均相对误差(%) 吐鲁番 89.23 42.95 13.16 吐鲁番农业站 89.25 42.83 14.61 托克逊站 90.23 42.85 9.49 鄯善站 88.60 42.77 13.78 库米什站 88.22 42.23 8.91 平均相对误差(%) 11.99 表1. SEBAL模型计算的相对误差值 图2. SEBAL值与P-M值散点图
system: [SYSTEM]

user: 5. 结果分析 5.1. 日蒸散量 吐鲁番市冬季寒冷干燥，蒸散发比较微弱，从多年蒸散平均来看，12月、1月、2月日蒸散量在0~2 mm/d；春季气温回升较快，3月份月均温即从2月份的零下回升至15℃左右，植被开始复苏，且春季多大风天气，5月份日蒸散量最大值已达6.4 mm/d；夏季蒸散量在0~9.0 mm/d，秋季蒸散量在0~3.7 mm/d。 从不同土地覆被类型来看，蒸散量差异性较大，日蒸散量极大值水体 > 耕地 > 林地 > 草地 > 城乡、工矿及居民用地 > 其他用地。水域主要为冰川积雪、河流、水库、季节性湖泊等，日蒸散量极大值出现在7月中下旬，在8.5~9.0 mm/d左右。耕地主要种植作物为葡萄、棉花、瓜，冬季棉花、瓜等一年生作物收获，葡萄掩埋过冬，蒸散发量主要集中在4~10月份作物生长季，占全年蒸散量99%，日蒸散量极大值出现在7月初，在8.3~8.8 mm/d之间。林地包括人工林地和自然林地，蒸散量极大值出现在8月初，在4.3~4.6 mm/d之间。草地蒸散量极大值同样出现在8月初，在3.5~3.9 mm/d之间。城乡、工矿及居民用地呈零星分布，多与耕地交错，受混合像元影响较大，日蒸散量极大值在3.4~4.5 mm/d之间。其他用地包括裸土、裸岩、沙漠、戈壁等，这些地方降水稀少，全年日蒸散量基本在0~2 mm/d。 5.2. 月蒸散量 月蒸散量时间变化规律与日蒸散量类似，2014年~2018年全年蒸散数据显示，蒸散量在各月份差异显著，主要集中在4~10月份，占全年总蒸散量的97.04%，从季节上看，春季占比22.06%，夏季占比62.25%，秋季占比14.83%，冬季占比0.86%。冬季月均蒸散量在3 mm以下，春季开始增大，3月、4月、5月平均蒸散量分别为5.46 mm、19.89 mm、43.12 mm，夏季达到峰值，5月、6月、7月平均蒸散量分别为67.14 mm、68.29 mm、67.41 mm，秋季逐渐降低，9月、10月、11月平均蒸散量分别为43.34 mm、17.26 mm、1.47 mm。 各季节典型月份蒸散量分布如图3，蒸散量较大区域集中在中部绿洲区及北部山区草地，其他平原区域多为荒漠、沙漠、戈壁，南部山区裸岩、裸土分布较广，干旱少雨，单位面积蒸散量小。 对于不同土地利用类型，月均蒸散量有较大差异，受冰川积雪占比、及混合像元影响，水域月均蒸散量在4~10月份低于耕地，最大值在63.9 mm。耕地受灌溉影响，在作物生长季节蒸散发强烈，6月份月均蒸散量达179.6 mm，由于吐鲁番种植两季瓜，每年在6月底7月初收获，6月份蒸散量稍高于7月，秋季作物收获后，蒸散量迅速下降。林地、草地蒸散量趋势相近，林地蒸散量8月份最大为67.1 mm，草地为59.4 mm。城乡、工矿及居民用地受混合像元中耕地影响，蒸散量6月最大为49.2 mm。其他用地蒸散量较低，月均蒸散量最大为10.9 mm。 图3. 2014~2018年各季节典型月份平均月蒸散量空间分布 5.3. 总蒸散量 2014~2018年吐鲁番全市范围蒸散情况如表2，平均总蒸散量为3.25 × 10 9 m 3 ，从蒸散量年际变化来看，总蒸散量有波动下降趋势，对于不同土地利用类型，耕地、草地、城乡、工矿及居民用地总蒸散量总体呈下降趋势，林地蒸散量则有上升趋势，水域及其他用地耗水变化不大，上述变化的原因主要包括：一是吐鲁番推行退地减水政策，逐步压减了耕地面积，同时通过改进灌溉方式提高了用水效率；二是地下水水位不断下降，草场面积有所缩减；三是当地节水政策促使工矿企业提高了中水利用率和居民节水意识；四是近年来吐鲁番人工造林面积不断扩大，林地用水及蒸散量均有所增加。 2014年~2018年全市单位面积平均年蒸散量46.56 mm，从不同土地类型来看，草地平均年蒸散量1.23 × 10 9 m 3 、其他用地1.06 × 10 9 m 3 、耕地7.54 × 10 8 m 3 ，分别占全年蒸散量的37.94%、32.66%和23.21%，林地1.11 × 10 8 m 3 ，城乡、工矿及居民用地8.01 × 10 7 m 3 ，水域9.74 × 10 6 m 3 ，而单位面积平均年总蒸散量耕地 > 城乡、工矿及居民用地 > 林地 > 草地 > 水域 > 其他用地。耕地由于人工灌溉单位面积年蒸散量最高，但由于面积只占全市范围1.26%，年总蒸散量低于草地和其他用地蒸散量；水域单位面积蒸散量较低的原因在于大部分水域为季节性水域或积雪区域，草地和其他用地单位面积年蒸散量不高，但由于面积分别占全市面积82.98%和14.06%，总蒸散量反而最高。 Table 2 Total evapotranspiration and different land use at Turpan city in 2014-2018 (m3 年份 总蒸散量 耕地 林地 草地 城乡、工矿 及居民用地 水域 其他用地 2014年 3.40 × 10 9 8.13 × 10 8 1.12 × 10 8 1.32 × 10 9 9.51 × 10 7 1.08 × 10 7 1.05 × 10 9 2015年 3.18 × 10 9 7.54 × 10 8 1.06 × 10 8 1.19 × 10 9 8.14 × 10 7 9.85 × 10 6 1.04 × 10 9 2016年 3.28 × 10 9 7.45 × 10 8 1.08 × 10 8 1.25 × 10 9 7.87 × 10 7 9.84 × 10 6 1.09 × 10 9 2017年 3.20 × 10 9 7.41 × 10 8 1.14 × 10 8 1.22 × 10 9 7.17 × 10 7 8.97 × 10 6 1.05 × 10 9 2018年 3.18 × 10 9 7.17 × 10 8 1.14 × 10 8 1.18 × 10 9 7.44 × 10 7 9.22 × 10 6 1.08 × 10 9 平均 3.25 × 10 9 7.54 × 10 8 1.11 × 10 8 1.23 × 10 9 8.01 × 10 7 9.74 × 10 6 1.06 × 10 9 百分比 100% 23.21% 3.42% 37.94% 2.47% 0.30% 32.66% 面积占比 100% 1.26% 1.02% 14.06% 0.51% 0.17% 82.98% 单位面积平均年蒸散量 (mm) 46.56 859.87 156.43 125.55 223.19 83.07 18.33 表2. 2014~2018年全市及不同土地类型蒸散量/m 3
system: [SYSTEM]

user: 6. 结论 本文以MODIS数据、气象观测数据、DEM等为基础，通过SEBAL模型对吐鲁番市2014~2018年蒸散量进行估算，分析蒸散量时空分布特征，结论如下： 1) 吐鲁番市日蒸散量时空差异较大，夏季日蒸散量最大，春季、秋季次之，冬季蒸散量最小；日蒸散量极大值：水域 > 耕地 > 林地 > 草地 > 城乡、工矿及居民用地 > 其他用地。 2) 不同土地利用类型的月平均蒸散量极大值：耕地 > 水域 > 林地 > 草地 > 城乡、工矿及居民用地 > 其他用地。 3) 年总蒸散量以草地、未利用地和耕地蒸散量为主，分别占全区总蒸散量的37.94%、32.66%和23.21%，林地、城乡、工矿及居民用地和水域总蒸散量占比较少，分别占3.42%、2.47%和0.17%，单位面积年总蒸散量耕地 > 城乡、工矿及居民用地 > 林地 > 草地 > 水域 > 其他用地。 4) 年总蒸散量97.04%集中在4~10月份，从季节上看，春季占比22.06%，夏季占比62.25%，秋季占比14.83%，冬季占比0.86% 利用FAO 56推荐的P-M公式和气象站观测数据对结果进行精度检验，平均误差为11.99%，拟合优度R 2 值为0.87。
system: [SYSTEM]

user: 摘 要 目的：通过对芳香化酶抑制剂联合基因重组人生长激素治疗大骨龄矮身材儿童疗效的分析，为大骨龄矮身材儿童的治疗提供临床依据。方法：选取2017年10月至2019年10月在青岛市妇女儿童医院内分泌代谢科诊断为矮身材、骨龄 ≥ 13岁的男童，同时应用来曲唑和基因重组人生长激素治疗，观察基线、治疗6个月、12个月疗程结束后的身高、体重、生长速率、骨龄、预测成年身高、睾酮、胰岛素样生长因子-1的变化及副作用。结果：芳香化酶抑制剂联合基因重组人生长激素治疗后患儿睾酮、胰岛素样生长因子-1水平明显升高，骨龄进展缓慢，生长速率、预测成年身高改善，差异有统计学意义(P < 0.01)。治疗过程中监测患儿均有正常男性化进程，治疗前后体重指数无统计学差异(P > 0.05)，且安全性及耐受性良好。结论：芳香化酶抑制剂联合基因重组人生长激素治疗大骨龄矮身材男童可延缓骨龄进展，改善预测成年身高，近期安全性及耐受性良好。 关键词 :来曲唑，基因重组人生长激素，大骨龄，矮身材，男童 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 矮身材是指在相似生活环境下，身高低于同种族、同性别和同年龄健康人群第3百分位数或平均身高减2个标准差者 [ 1 ]。基因重组人生长激素(recombinant human growth hormone， rhGH)的应用可有效提高矮身材患儿的生长速率、最终成年身高的疗效已得到公认。但对于骨龄偏大的矮身材儿童，由于生长潜能不足，治疗空间有限，导致成年身高不理想。近年有研究表明，芳香化酶抑制剂(aromatase inhibition， AI)联合基因重组人生长激素治疗大骨龄矮身材儿童可有效延缓骨龄进展，改善预测成年身高，但是目前国内外相关研究样本量较少，随访时间较短，对其疗效及副作用存在一定的争议。目前国内外对大骨龄的定义尚未明确，结合文献报道 [ 2 ]，本研究以骨龄13岁以上为大骨龄，分析AI联合rhGH治疗大骨龄矮身材儿童的疗效及安全性。