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

user: 3. 数据准备 3.1. 实验数据 本文主要以国内某公司自动售卖机的销售数据为研究对象，涉及销售设备154台，覆盖该城市大部分区域，时间主要为201801~201812销售数据，该数据集覆盖地区较广，设备数目较多，数据量相对完整。原始数据集包含主要数据信息如表1所示。 由于自动售卖机的销量受到许多因素的影响，本次考虑的外部因素主要时气象数据，节假日等，因此本文爬取分析了可获得的外部环境数据，主要数据信息如表2所示。 Table 1 Main fields of vending machine sales dat 序号 字段 数据描述 1 付款时间 例：2019/9/17 19:44 2 商品id 例：2914 3 商品类别 例：饮料，水等 4 设备id 例：73c5fd40940e4e0eb25625a1092d1717 5 实际价格 例：5.5 6 经纬度 例：31.091778，121.520671 7 地铁线 例：1号线，2号线 表1. 自动售货机销售数据主要字段 Table 2 Main fields of external data of vending machin 序号 字段 数据描述 1 节假日 例：休息日，节假日，工作日 2 气象数据 例：温度，降水量，湿度等 表2. 自动售货机外部数据主要字段 3.2. 数据预处理 由于目前该企业正处于上升阶段，自动贩卖机设备的数量处于不断上升的过程，所产生的数据由于各种原因会存在部分、整体缺失或者不合理的情况，无法直接用于分析与预测研究，因此需要对数据进行预处理操作，主要包含以下步骤： 1) 数据归约： 由于销量数据太过庞大，并且每条销量数据记录之间的间隔并不确定，因此需要对销量数据进行聚合处理。提取出销售数据中的关键信息，按照合适的时间粒度对初始数据进行整理，完成数据归约操作，得到售卖机销量数据表，本次选取的预测时间粒度为一天，反映了一天内用户购买的饮品数量。 2) 数据清洗： 对于接入的数据，很多都包含了缺失、异常或者类型不一致的字段数据，此时需要对数据进行数据清洗，数据清洗就是发现这些类型的数据并纠正这些类型的数据，其中包含数据一致性、无效值以及缺失值的处理，处理方法有补全、删除等。本次由于自动售卖机机器损坏或数据为接入等原因导致销量存在较多的缺失值与部分异常值；对于该部分缺失值与异常值军进行异常处理，采用离散系数和小波分析相结合的方式进行识别，识别后采用对应方法处理。 3.3. 数据分析 依据主观分析，温度对饮品的销量会产生较大的影响，此次采用散点图的形式判断温度与销量之间是否存在相关关系。 图1 为2018年销量与温度对应的关系图表，其中纵坐标表示不同温度对应的平均日销量，从图中可以发现饮品的销量与温度之间存在明显的相关关系，尤其是当温度越高时，相关关系越明显。因此可以将温度作为销量预测的一个影响变量。 图1. 销量与温度关系图
system: [SYSTEM]

user: 5. 结论 对于时间序列的预测而言，不同时间长度的销量序列与不同粒度的预测分析对预测结果都会产生一定的影响。对于本次的自动售卖机的销量预测而言，采用时间长度为一年的销量序列作为基础时间序列进行时，进行总体预测的预测精度是最优的，但是对于单机预测而言，需要单台机器进行具体分析，对于不同的机器可以适当的调整基础时间序列的长度来进行短期预测，但就单台机器的整体预测效果而言，当时间序列长度为4个月或者6个月时，预测精度最优。
system: [SYSTEM]

user: 摘要 目的：探讨卡氏肺孢子菌肺炎发生在艾滋病与非艾滋病患者中的影像学特点。方法：对确诊的26例患者，以入院首次肺部影像作为研究资料进行回顾性分析。结果：26例患者中，HIV患者呈现磨玻璃阴影、肺气肿影及月弓征居多，分别占比为21例(87.5%)，8例(88.9%)和5例(83.3%)，双轨征及铺路石征较少；非HIV患者，双轨征及铺路石征居多，分别为5例(71.4%)和3例(100%)。结论：PCP患者的胸部CT影像学表现复杂多样，多呈双肺弥漫对称性分布，磨玻璃影、肺囊肿和月弓征在HIV/ARDS患者占比较多，而双轨征及铺路石征在非HIV患者中占比较多。
system: [SYSTEM]

user: 1. 引言 肺孢子菌肺炎(Pneumocystis pneumonia， PCP)，又称卡氏肺孢子虫肺炎，卡氏肺囊虫肺炎，是由肺孢子菌引起的间质性浆细胞性肺炎，为条件性肺部感染性疾病。常常归属于真菌机会性肺部感染的一类。过去曾经认为肺孢子菌具有原虫的形态学特征而被学界认为是一类新型的锥体型原虫。直到1988年，对肺孢子菌rRNA亚型充分分析确定了该病原体属真菌类，此后相继众多学者行基因检测也证实此病原体确属于真菌中的子囊菌类 [ 1 ]。肺孢子菌本身具备特定的生物遗传特性和宿主特异性，学术界将引起人类肺孢子菌肺炎的菌属称为P. jirovecii [ 1 ]。 自身获得性免疫缺陷综合征(acquired immunodeficiency syndrome， AIDS)合并PCP的病例报道逐年增多，越来越多的被临床重视。对AIDS患者的急性发病率及死亡率分析，合并PCP的患者占比较高。但近年来随着科学技术水平的突飞猛进，尤其是器官移植的广泛实施，人类免疫缺陷病毒(human immunodeficiency virus， HIV)阴性的普通患者及免疫抑制患者病例层出不穷，如结缔组织疾病(尤其是胶原血管病)、皮质类固醇长期使用、糖尿病、恶性肿瘤(尤其实体及血液系统恶性肿瘤)等等。综合分析，上述两类群体不仅造成PCP的原因有异，而且在两类人群的发病率和病死率方面亦有所不同，参考报道的相关统计数字，AIDS并发PCP的患者急性加重住院比例明显升高，重症患者死亡率进行性升高，可达29%~62%；比较而言，非AIDS并发PCP患者的病死率更加严峻，可达48%~67% [ 2 ]。PCP一经诊断，病情迅速进展，早期恰当治疗可显著降低病死率。此外肺孢子菌在免疫功能正常人群中常常存在定植现象，可能合并呼吸道传播的可能 [ 3 ]。就呼吸科常见病而言，存在免疫功能缺陷患者肺孢子菌定植的比例较高，另一方面最常见的肺孢子菌病传染源是具有免疫功能正常、无症状的健康人群。
system: [SYSTEM]

user: 2. 资料与方法 2.1. 入选标准 选取研究对象为青岛大学第一附属医院从2010年3月至2020年3月确诊的PCP (肺孢子菌肺炎)的26例门诊及住院患者，PCP诊断符合下列标准第1项或2~7中任何4项。① 支气管镜活检或支气管肺泡灌洗检测PC阳性；② ARDS诊断标准成立 [ 3 ]；③ 血清淋巴细胞亚群检测CD4 + T淋巴细胞计数 < 200个/µl [ 3 ]；④ 伴有咳嗽，以干咳为主、可伴呼吸困难、发热、胸痛及体重下降等症状，但查体无阳性胸部体征；⑤ 胸部影像学可见典型征象；⑥ 给予经验性抗PCP治疗临床症状有效；⑦ 反复3次血清乳酸脱氢酶检测均升高 [ 4 ]。要求ARDS诊断以符合2011年中华医学会感染病学会所修订的诊断为准 [ 3 ]。 2.2. 一般资料 入组患者共26例均诊断为PCP，男性21例，女性5例，男女之比为4.2:1。年龄22~68岁，平均42.3岁。HIV感染途径包括：性传播16例(61.5%)，其中异性传播12例，同性性传播4例；血液传播2例(7.7%)；感染途径不确定8例(30.8%)。 临床表现：发热首诊患者23例(88.5%)，活动性呼吸困难19例(73.1%)，静息憋喘4例(15.4%)，干咳21例(80.8%)，紫绀6例(23.1%)，肺喘鸣音9例(34.6%)。 2.3. 方法 对过去10年间符合PCP诊断的患者初始确诊时的胸部CT影像作为研究资料进行病例对照研究分析，内容包括影像征象及病变具体部位。在分析时以肺叶为单位，双肺均按照上叶/上中叶、下叶、全肺来统计描述。
system: [SYSTEM]

user: 3. 结果 3.1. 全部PCP患者胸部CT影像分布特征 在入组的26例病例中，弥漫性磨玻璃影占比最大，共计24例(92.3%)，其他影像依次为网格影12例(46.2%)，肺气肿影9例(24.6%)，双轨征7例(26.9%)，斑片影8例(30.8%)，弓月征6例(23.1%)，铺路石征3例(11.5%) (如图1~6)。具体分布中肺气肿以上肺叶分布为主，磨玻璃阴影以全肺叶居多，月弓征及斑片阴影以肺下叶为主，未见合并胸腔积液。与文献报道基本一致，见表1。 图1. 双肺呈磨玻璃密度影，伴对称、弥漫性分布 图2. 双肺磨玻璃影散在分 图3. 双肺磨玻璃密度影伴条索状影、网格影 图4. 双肺磨玻璃密度影，左肺上叶见斑片实变影 图5. 多发小叶间隔增厚，铺路石征改变 图6. 双肺磨玻璃影，伴右上肺气肿影 Table 1 The distribution of chest CT imaging in 26 PCP patients [cases (%) 形态 例数 上/中叶 下叶 全肺 磨玻璃影 24 (92.3%) 2 (8.3%) — 22 (91.7%) 网格 12 (46.2%) — 3 (25%) 8 (75%) 肺气肿 9 (24.6%) 7 (77.8%) — 2 (22.2%) 双轨征 7 (26.9%) 3 (42.9%) — 4 (57.1%) 斑片影 8 (30.8%) — 3 (37.5%) 5 (62，5%) 月弓征 6 (23.1%) — 5 (83.3%) 1 (16.7%) 铺路石征 3 (11.5%) — 1 (33.3%) 2 (66.7%) 表1. 26例PCP患者胸部CT影像学分布[例(%)] 注：“—”表示此肺叶无相关影像学征象。 3.2. PCP患者原发疾病谱特征分析 对26例诊断为PCP的患者通过病史询问及既往病历统计，确诊为HIV/ARDS的患者为19例(73.1%)，非HIV患者7例(26.9%)，其中肾移植患者1例(3.8%)，淋巴瘤患者1例(3.8%)，干燥综合征患者1例(3.8%)，嗜酸细胞性肉芽肿性多血管炎(EGPA) 1例(3.8%)，2型糖尿病患者1例(3.8%)，2例非HIV患者未明确原发病因。见表2。 Table 2 Etiology distribution of 26 PCP patients [cases (%) 疾病 例数 HIV 19 (73.1%) 非HIV (7; 26.9%) 器官移植(肾移植) 1 (3.8%) 恶性肿瘤(淋巴瘤) 1 (3.8%) 自身免疫性疾病 2 (7.7%) (干燥综合征 + 血管炎) 0 (0%) 2型糖尿病 1 (3.8%) 未知 2 (7.7%) 表2. 26例PCP患者病因分布[例(%)] 3.3. HIV及非HIV中PCP患者胸部CT影像分布特征 26例患者中，HIV患者呈磨玻璃影及肺气肿居多，分别占比为18例(94.7%)和8例(42.1%)，月弓征及斑片影各5例(26.3%)，双轨征及铺路石征较少，分别占比为2例(28.6%)和0例(0%)。对非HIV患者影像分析中双轨征及铺路石征居多，分别为5例(71.4%)和3例(100%)，而磨玻璃阴影、肺气肿影及月弓征占比较少，分别为3例(12.5%)、1例(11.1%)和1例(16.7%)。对于网格影及斑片影无明显特征性分布人群。其中，磨玻璃影、网格、双轨征、铺路石征差异有统计学意义(表3)。 Table 3 The distribution of chest CT images of PCP patients in HIV and non-HIV [cases (%) 类型 HIV (19) 非HIV (7) P值 磨玻璃影(24) 18 (94.7%) 3 (42.9%) 0.01 网格(12) 5 (26.3%) 7 (100%) 0.005 肺气肿(9) 8 (42.1%) 1 (14.3%) 0.357 双轨征(7) 2 (10.5%) 5 (71.4%) 0.006 斑片影(8) 5 (26.3%) 3 (42.9%) 0.635 月弓征(6) 5 (26.3%) 1 (14.3%) 1.0 铺路石征(3) 0 (0%) 3 (42.9%) 0.013 表3. HIV及非HIV中PCP患者胸部CT影像分布[例(%)]
system: [SYSTEM]

user: 4. 讨论 1. PCP在临床属相对少见病，是免疫抑制患者严重的感染性疾病之一，常常伴有很高的病死率，据统计HIV继发PCP的病死率为10%~20%，而非HIV患者继发PCP的病死率可达30%~60% [ 1 ]。过去PCP的报道主要源于HIV/AIDS患者，分布于各年龄阶段。但随着临床技能和科研水平的进步，越来越多的PCP患者来自于其他免疫功能抑制患者(如器官移植者、恶性肿瘤、长期应用激素、糖尿病及血液病等)，曾有文献报道诊断为PCP的2587例患者，其中HIV/AIDS相关PCP 2061例(79.7%)，器官移植术后相关PCP 316例(12.3%)，恶性肿瘤相关PCP 78例(3.0%)，自身免疫性疾病相关PCP 11例(0.4%) [ 5 ]。 2. 鉴别诊断：1) 弥漫性肺泡出血：胸部X线检查表现为斑片影，实变影及磨玻璃影，中下肺野多见，HRCT影像可见腺泡填充影弥漫均匀分布，初期主要表现为磨玻璃阴影，中后期表现为结节影，当出血量较大时可表现为大片浸润影或实变影。在DAH (弥漫性肺泡出血)复发病例中，多数以间质性病变为主，最常见的是肺纤维化 [ 6 ]。2) 肺泡蛋白沉积症：常规胸部X线平片表现为多种多样，最常见的是磨玻璃样密度影，从模糊结节影到大片实变影，病变呈中央或周边分布均可见，上下肺分布未见明显差异。高分辨CT影像表现为磨玻璃密度影合并网格状影或斑片状影，病变与周围肺组织边界清晰呈地图样改变。可显示肺间质改变，部分肺泡实变，小叶间隔增厚呈铺路石征或铺路石征。典型的支气管肺泡灌洗呈“牛奶样”或“泥浆样”改变，静置后出现沉淀 [ 7 ]。3) 特发性肺纤维化：典型的IPF双肺病灶多位于胸膜下或基底部，伴有异常的网格状影，同时可见蜂窝样改变，伴或不伴有牵拉性支气管扩张，最后应当除外不符合UIP型改变① 双肺病灶主要出现于中，上肺野；② 病变主要延支气管血管束分布；③ 广泛磨玻璃样影；④ 大量微结节(双侧，上肺野分布为主)；⑤ 散在囊状病变；⑥ 弥漫性马赛克征/气体陷闭；⑦ 支气管肺段/肺叶实变影 [ 8 ]。 3. 查阅相关文献HIV/AIDS合并PCP的胸部HRCT影像包括① 两肺弥漫磨玻璃阴影：早期以弥漫浸润的实质病变为主，随着病情加重逐步出现多发磨玻璃阴影，磨玻璃影作为在其最具特征性的CT影像。此时双肺出现透亮度降低，病变以肺门周围及两肺中、下部分布为主，以间质改变为主 [ 9 ]；② 肺气肿：大多早期见于肺上叶，后期中叶及下叶也可见，肺气肿影源于肺泡与肺间质炎性改变和成纤维细胞增殖导致气道重塑，当肺气肿不断扩大扩大甚至发生破裂，破裂时可伴有自发性气胸形成；肺实质多发的囊性改变，提示肺间质的终末期样改变，常常是一种不可逆的表现 [ 10 ] [ 11 ]；③ “月弓征”影像：常常出现在发病早期或进展期，主要为胸膜下可见的正常肺组织形成的类似新月形样或弓形样豁免区域，表示该区域为尚未因病变受累的肺组织，目前学界部分专家认为“月弓征”与HIV继发PCP有着较高的相关性 [ 9 ]。 同时结合文献及本次回顾性研究数据分析非HIV合并PCP的胸部螺旋CT表现主要有① 双轨征：该征象可见于双上肺，或全肺受累，本研究患者中出现双轨征者共7例(26.9%)，非HIV占有5例(71.4%)，根据患者基础病史及部分病理提示双轨征即支气管扩张所致，与患者局部起到慢性炎症导致官腔渗出物增多及气道壁水肿相关。治疗中除了积极治疗PCP外还应覆盖结构性肺病的基础治疗。因此在一些文献中提出非HIV患者继发PCP临床症状及死亡率均高于HIV相关性PCP患者。② 铺路石征：本研究铺路石征有3例(11.5%)，铺路石征在HRCT上表现为地图状分布，重叠有网状的光滑的细线磨玻璃影，细线影多为小叶间隔水肿及慢性炎症增厚，该征象可见于很多疾病。 综上所述，无论是HIV/ARDS还是非HIV患者，若在随诊或近期出现肺部感染性病变时，HRCT扫描可见前述双肺间质性病变、如弥漫性磨玻璃影、新发肺气肿，月弓征，双轨征及铺路石征等表现应高度警惕PCP感染的可能，尽可能做到早期诊断PCP，早期治疗，从而改善预后；同时可以通过不同的螺旋CT影像来提示HIV及非HIV疾病的早期筛查 [ 12 ]。
system: [SYSTEM]

user: 摘要 CO 2 近混相驱对地层压力、原油物性和注入气纯度等要求不高，降低其操作成本，更具有应用可行性。本文基于低渗油藏E区块油藏和非纯CO 2 气性质，采用细管模拟及数值模拟等手段，研究近混相驱界限、CO 2 纯度下限要求以及应用可行性。研究结果表明，目标油藏最小近混相驱压力(MNMP)为(0.80~0.84) MMP、CO 2 纯度下限为57.6%；CO 2 纯度80%下采出程度比CO 2 纯度40%下提高了12.53%，仅比纯CO 2 混相驱采出程度降低了3.88%，进一步提高CO 2 纯度对采收率的提高幅度很有限，说明CO 2 近混相驱开发效果接近于混相驱。考虑到提纯CO 2 成本，推荐近混相注入气的合理纯度为80%。CO 2 近混相驱10年开发期采出程度高达58.28%，比水驱提高了17.67%，目标区块实施非纯CO 2 近混相驱开发可行性高。研究成果对我国低渗透油藏应用非纯CO 2 近混相驱开发技术具有重要指导意义。
system: [SYSTEM]

user: 1. 引言 国内外大量研究和矿场应用 [ 1 ] [ 2 ] 表明，混相驱其驱油效率最高，对混相驱技术的研究与应用比较成功。对于中国大多数油藏属于陆相沉积，其油品和物性较差，还缺乏纯CO 2 气源，导致最小混相压力(MMP)较高，难以满足混相条件。与混相驱相比，CO 2 近混相驱对地层压力、原油物性和注入气纯度等要求不高，大大提高其应用可行性 [ 3 ] [ 4 ]。如果将非纯CO 2 注入油藏即能提高原油采收率，还能减轻温室气体排放，因此，探索非纯CO 2 近混相驱更有技术和环保意义。Dong等人对CO 2 近混相驱进行了实验研究 [ 5 ] [ 6 ]，其驱油效率比混相驱并未明显下降。Parvazdavani采用室内物模到油藏规模模拟手段，研究了近混相驱效果 [ 7 ]。尚宝兵、杨胜来等人 [ 8 ] [ 9 ] [ 10 ] 研究了杂质气体对二氧化碳驱最小混相压力和原油物性的影响，得出了回注CO 2 伴生气可行性和注入方式，证明其应用。张贤松等 [ 11 ] 针对某海上低渗油藏E区块富含CO 2 (24%~90 mol%)，最小混相压力(MMP)高于地层压力，探讨了非纯CO 2 近混相驱的控制条件。与混相驱判别不同，近混相驱是介于混相驱和非混相驱之间的区域，如何界定该区域以及对CO 2 纯度要求具有应用价值。本文基于E区块5井油气样和油藏性质以及细管实验结果，建立细管模型和典型油藏模型，采用细管模拟及数值模拟等手段，研究了非纯CO 2 最小近混相压力及对CO 2 纯度的要求，以及应用CO 2 近混相驱开发可行性。结果表明，目标油藏近混相驱CO 2 纯度区间为下限为57.6%~99%，其下限为57.6%；CO 2 纯度80%下近混相驱的采出程度比非混相CO 2 (纯度40%下)提高了12.53%，仅比混相驱(CO 2 纯度100%)采出程度下降了3.88%，说明CO 2 近混相驱开发效果接近于混相驱。进一步CO 2 纯度对采出程度的提高幅度较小，考虑进一步提纯CO 2 的成本更高，推荐目标区块近混相驱注入气的CO 2 纯度为80%。研究成果对海上低渗油藏气驱高效开发提供技术依据，对陆地低渗油田CO 2 提高采收率技术应用也具有指导意义。
system: [SYSTEM]

user: 2. 最小近混相压力与CO 2 纯度下限 E区块油藏富含CO 2 伴生气，可作为CO 2 驱的气源，回注油藏实施CO 2 驱。其油品及气体组成分析见表1。为使非纯CO 2 注入气与地层原油达到近混相程度，首先要研究不同纯度CO 2 对近混相区间及纯度下限所要求。本文基于细管实验结果 [ 11 ]，建立长细管模型，在模型始端和末尾端分别设置一口注入井和油井，细管模型图如图1所示。 Table 1 Main property parameters of crude oil and well flui 气油比(闪蒸气/闪蒸油) 107.5 m 3 /m 3 地层油体积系数 1.3357 地层油密度 0.6986 g/cm 3 地层油粘度 0.29 mPa.s C 1 +N 2 42.92 mol% CO 2 +C 2 ~C 10 27.33 mol% C 11+ 29.75 mol% 表1. 原油和井流物主要性质参数 图1. 细管模拟模型图 在低于MMMP下，即在近混相及非混相段加密压力点，在混相和非混直线段之间，细管模拟得到一条新的直线段即近混相段，其与混相和非混直线段交叉得到两个交点，定义这两点之间的压力区间为近混相区间。以注纯CO 2 为例，设定不同压力，注入1.2 PV的纯CO 2 ，细管模拟计算得到不同压力下对应的驱油效率。作不同压力与其驱油效率关系图，得到明显的三段线性关系，如图2。将中间段线性关系做直线延长，与其前后的直线段交点所对应压力点为近混相压力区间，其中，与非混相直线段交点所对应的压力为近混相压力下限，定义为最小近混相压力(MNMP)，得到E区块纯CO 2 最小近混相驱压力为27.05 MPa，近混相压力区间为27.05 MPa~33.64 MPa。 如同注纯CO 2 驱一样，通过细管模拟得到某一CO 2 纯度下的驱油效率与注入压力关系，获得不同CO 2 纯度下对应的MMP和MNMP，见表2。根据不同CO 2 纯度下的MNMP和MMP，进一步得到MNMP为(0.8~0.84) MMP。 根据表2数据，绘制MMP和MNMP与 P C O 2 (CO 2 纯度)之间关系曲线，发现MMP和MNMP均与CO 2 纯度呈线性关系： M M P = − 0.1703 P C O 2 + 49.86 (1) M N M P = − 0.1516 P C O 2 + 41.73 (2) 对纵坐标压力平行于横坐标做直线，与MMP和MNMP线相交，交点所对应的CO 2 纯度即为该压力下近混相驱的CO 2 纯度上下限，如图3所示。如目标油藏平均地层压力为33 MPa，近混相驱的CO 2 纯度区间范围为57.6%~99.0%，其中，CO 2 纯度下限57.6%。目标油藏产出气CO 2 含量为24%~90%，CO 2 近混相驱纯度窗口区间较宽，因此，注入伴生非纯CO 2 近混相驱的可行性和潜力较大。 图2. 驱油效率与注入压力关系图(纯CO 2 ) 图3. 近混相驱CO 2 纯度图版 Table 2 MNMP and MMP with different CO2 purit CO 2 纯度 % MNMP MPa MMP MPa MNMP /MMP 混相程度 100 27.05 33.64 0.80 混相驱 90 28.04 34.78 0.81 近混相 80 29.27 35.63 0.82 70 30.74 37.01 0.83 60 32.73 39.86 0.82 50 34.30 41.03 0.84 40 35.42 42.95 0.82 非混相 24 38.24 46.28 0.83 表2. 不同CO 2 纯度下的MNMP和MMP 根据E区块5井和4井四个井段的资料，得到不同井段近混相CO 2 纯度下限稍有差异，分析其原因主要与各个井段所处的油藏条件有关。随井深增加，油藏温度和压力增加，CO 2 近混相驱纯度下限也随之增加，见表3。 Table 3 Comparison of reservoir conditions and oil propertie 井号 5井1段 4井2段 5井2段 4井3段 纯度下限，% 64 52 43 34 油藏中深，m 3319 3324 3489 3503 油藏温度，℃ 112.1 114.2 116.7 122.5 油藏压力，MPa 31.96 33.85 34.04 36.53 表3. 油藏条件及原油物性对比
system: [SYSTEM]

user: 3. 不同条件和方式的开发效果 采用Eclipse数值模拟软件，建立典型油藏模型，网格尺寸为10 m × 10 m × 20 m，初始地层压力33 Mpa，初始含油饱和度为60%。纵向上设置三层，模型参数见表4。建立一注一采的布井方式，如图4所示。 Table 4 Permeability and porosity parameters of typical model 层名 孔隙度 渗透率mD 层厚m I-1 0.25 10 20 I-2 0.22 5 20 I-3 0.15 0.8 20 表4. 典型模型渗透率和孔隙度参数 图4. 典型油藏模型垂向渗透率示意图 3.1. 不同混相程度 设计不同CO 2 纯度，数值模拟得到CO 2 气驱过程中的含油饱和度变化。图5为3种CO 2 纯度下(不同混相状态)驱替中同一时刻的含油饱和度变化图，CO 2 纯度为40%代表非混相驱，CO 2 纯度为80%代表近混相驱，CO 2 纯度为100%代表混相驱。可以看出，CO 2 近混相驱其驱油效率和波及范围上居中，更接近于混相驱。 图5. 不同CO 2 纯度下气驱含油饱和度图 表5为不同CO 2 纯度下气驱采出程度效果对比，当CO 2 纯度为40%时，累计采出程度为42.27%，而当CO 2 纯度80%下采出程度达到54.8%，是前者的1.29倍；纯CO 2 混相驱采出程度为58.68%，比CO 2 纯度80%下近混相驱仅提高了3.88%，提高幅度有限，纯度高于80%后，进一步提高纯度对采出程度的提高幅度较小，这说明CO 2 近混相驱开发效果接近于混相驱。考虑到CO 2 提纯成本，推荐E区块CO 2 驱近混相注入气的理想纯度为80%。 Table 5 Simulation results of different CO2 purity flooding scheme CO 2 纯度% 混相方式 累计采油量m 3 采出程度% 40 非混相 121737.6 42.27 60 近混相 142300.8 49.41 80 近混相 157824.0 54.80 90 近混相 163958.4 56.93 100 混相 168998.4 58.68 表5. 不同CO 2 纯度气驱方案模拟计算结果 3.2. 不同开发方式 根据油藏特征与上述研究结果，采用如下注气参数 [ 12 ] [ 13 ]：CO 2 纯度为80%，注气速度为0.08 HCPV/a，数值模拟计算CO 2 近混相驱和水驱两种开发方式10年开发期内的采出程度，计算结果见表6。CO 2 近混相驱开发10年期采出程度高达58.28%，比水驱提高了17.67%，表明非纯CO 2 近混相驱可获得更好的开发效果。 Table 6 Comparison of effects with different development pattern 方式 累计采油量m 3 累计增油量m 3 采出程度% 近混相驱 167846.4 94924.8 58.28 水驱 116956.8 44035.2 40.61 表6. 不同开发方式开发效果对比
system: [SYSTEM]

user: 4. 结论与建议 (1) 近混相驱最小近混相压力为(0.8~0.84) MMP，CO 2 纯度在57.6%以上即可达到近混相，目标区块满足CO 2 近混相驱潜力空间大。 (2) 高于80%后，继续增大CO 2 纯度，其采出程度提高幅度较小，CO 2 近混相驱注入气的合理纯度为80%。CO 2 纯度80%下近混相驱采出程度比纯度40%下非混相驱提高了12.53%，比纯CO 2 混相驱仅下降了3.88%，CO 2 近混相驱开发效果接近于混相驱。 (3) 目标区块CO 2 近混相驱开发10年期采出程度高达58.28%，比水驱提高了17.67%，非纯CO 2 近混相驱具有很高的应用可行性。
system: [SYSTEM]

user: 摘要 为研究降雨条件下淤积层的降雨入渗规律及评价淤积水量，基于内蒙古达拉特旗哈他土沟砒砂岩区淤积层为研究对象，开发一套模拟降雨条件下一维土柱的入渗装置，对淤积土柱做了4组不同降雨强度下的一维垂直入渗试验，得到不同降雨强度下垂直土柱的入渗率时程曲线、浸润峰深度时程曲线及监测点体积含水率的变化规律。结果表明：1) 降雨强度对垂直土柱的入渗影响较大，当降雨强度小于土柱最小入渗能力时，入渗率等于降雨强度；当降雨强度大于土柱最小入渗能力时，入渗率时程曲线呈无压入渗、有压入渗和饱和入渗三阶段变化。2) 不同降雨强度下，土柱出现积水点和饱和点的历时不同，降雨强度越大出现积水点和饱和点的时间越短，有压入渗阶段越长。3) 在同一降雨强度下，监测点距土柱上表面越远，其体积含水率时程曲线越密集；而同一监测点，降雨强度越大，其体积含水率时程曲线越稀疏。
system: [SYSTEM]

user: 1. 引言 砒砂岩区是黄河粗沙的重要的来源之一，也成为黄土高原水土保持治理难度最大的区域之一，被国内外专家称为“世界水土流失之最”。建造在黄土高原砒砂岩区的淤地坝，是拦蓄降水、减少进入黄河的泥沙的重要水土保持措施之一，起到调节径流、滞洪拦沙 [ 1 ]、增加土壤入渗的作用。降水通过类似于海绵体的淤地坝淤积层，全部或部分地表向下流动，并且在淤积层中运动和储存 [ 2 ] 从而形成土壤水库，在防洪保安全的前提下，充分高效利用淤地坝淤积层拦蓄雨洪资源的潜力，提高水资源利用率，对缓解半干旱生态脆弱区水资源供需矛盾具有十分重要的作用，因此，开展降雨条件下砒砂岩区淤地坝淤积层入渗规律研究具有重要的现实意义。 G. A. Ampt [ 3 ] 最早提出了入渗理论，并由此假定建立了Green-Ampt入渗模型。J. R. Phihp [ 4 ] 认为在入渗过程中任意时刻的入渗率与时间呈现幂级数关系。王全九等 [ 5 ] 对比分析了Philip入渗模型和Green-Ampt入渗模型，并通过维垂直入渗试验描述了2个模型中参数之间的关系。毛丽丽等 [ 6 ] [ 7 ] 根据水量/物质平衡原理和水平土柱中土壤剖面含水率分布所遵循的活塞假定，计算土壤入渗性能，并通过室内水平土柱试验计算得到了土壤的入渗性能曲线。J. Almedeij和I. I. Esen [ 8 ] 对Green-Ampt入渗模型进行了修改与完善，并将修改后的Green-Ampt入渗模型应用于案例分析，最后与Mein-Larson模型进行了比较。戴经梁等 [ 9 ] [ 10 ] 通过室内入渗试验，对黄土的积水入渗和水平入渗进行了研究，分别得到了不同压实度下黄土的积水入渗规律和水平黄土的渗水特性。陈洪凯等 [ 11 ] 对重庆境内的侏罗系强风化泥岩进行了6种降雨强度下的入渗试验研究，提出了可表征该类泥岩入渗过程的降雨入渗公式。朱伟等 [ 12 ] 通过室内降雨入渗土柱试验，探讨了反映降雨入渗量的有限元计算方法，并通过数值模拟揭示了室内试验未能明确的降雨入渗规律。王春颖等 [ 13 ] 通过室内层状夹砂土柱一维薄层积水入渗试验和相应情况下均质土柱的对照试验，研究了夹砂层对入渗强度、湿润锋行进和沿程土壤含水率变化的影响。苗强强等 [ 14 ] [ 15 ] 对非饱和含黏砂土分别做了毛细上升试验和一维入渗试验，分别得到了含黏砂土含水率和吸力随时间的变化规律及渗透曲线的变化规律。上述研究对非饱和土体入渗规律的分析具有重要的现实意义。但是，目前大多数研究均是针对水平土柱的稳定入渗过程，没有考虑重力作用的影响，而对于非饱和淤积土在降雨条件下的垂直入渗规律的研究较少。本次自主设计加工了一套一维垂直土柱入渗模型试验装置，基于该试验装置研究了降雨条件一维淤积土柱的入渗规律。
system: [SYSTEM]

user: 2. 土壤水分传感器的标定与数据采集 2.1. EC-5土壤水分传感器原理 试验所用水分传感器为EC-5土壤水分传感器，EC-5是由传感器壳体以及探针组成，通过测量介质的介电常数来确定体积含水量(VWC)。其测量范围：矿物土±3%VWC，大多数矿物土，高达8 dS/m ± 1%~2% VWC，土壤特定校准：岩棉 ± 3% VWC，0.5~8 dS/m；灌封土：±3% VWC，3~14 dS/m：0.1% VWC (矿物土)，0.25% VWC(岩棉)；测量时间10 ms；工作温度−40℃至+50℃。 EC-5土壤水分传感器是通过使用电容/频域技术测量介质的介电常数来确定体积含水量(VWC)，其内部是由方波信号发生电路、RC充放电电路以及时间电压转换电路等三部分组成，它他们之中的传感器部件电容C是用耐腐蚀的PCB为模板、待测土壤为电解质、土壤介电常数随着水分的变化表现为电容量的一种变化。当测试电极几何因数、电阻R、电压输出V是定值的时候，土壤介电常数可以由充电时间来决定，然后通过时间电压转化电路得到输出电压。所以说输出电压的变化就是所反映的土壤水分的变化。具有结构小巧、简单，防水能力强，耐腐性强，测量精度高，性能可靠，受土壤含盐量影响较小等优势。同时EC-5土壤水分传感器的温度补偿功能是其他传感器所无法比拟的。因此，该传感器在土壤含水率的测量中得到广泛应用。 2.2. EC-5土壤水分传感器的标定 根据EC-5土壤水分传感器的原理可知，对同一土体，传感器的输出信号值与被测土壤的体积含水率之间具有对应的映射关系；因此，试验前对所用的5支传感器分别编号为S1，S2，S3，S4和S5，然后将试验所用淤积土加水分别配制体积含水率为5.6%，11.3%，15.1%，34.9%和44.9%的试样5个，最后对不同编号的土壤水分传感器进行标定，5支传感器的输出信号值与体积含水率如图1所示。 图1. 水分传感器标定结果 图2. 降雨入渗试验设备图 由图2可知，S1，S3，S4，S5传感器的输出信号值与体积含水率具有较好的线性关系，S3传感器因标定过程中的失误，在体积含水率为34.9%时的读数具有较大出入值，线性拟合时舍去该含水率下数据输出值。 2.3. CR1000数据采集器 CR1000数据采集器是美国Campbell公司设计和制造的一款智能数据采集仪，扫描速率能够达到100 Hz，拥有模拟输入、脉冲计数、电压激发转换、数字等多个端口。可以在恶劣环境、电池供电环境下进行精确测量，并以表格形式存储的数据带有时间标记和记录标号，有可时实查看数据，画图并分析数据软件。当CR1000数据采集器从主电源断电后，电池支持SRAM内存和时钟，可以确保数据、程序和准确时间等信息不会丢失。CR1000所具有的高精度性、高适应性、高可靠性以及价格合理等特点，使其广泛应用于气象观测、农业研究、土壤水分研究、通量观测等领域。
system: [SYSTEM]

user: 3. 试验装置与试验方案 3.1. 试验装置 模拟降雨条件下淤积体垂直入渗的试验装置包括垂直土柱制作系统、数据采集系统及模拟降雨系统，如图2所示。 垂向淤积土柱的制作采用外径300 mm，内径280 mm，管长2000 mm的有机玻璃管，控制某一干密度分层进行击实；在试验前将有机玻璃管按每隔50 mm进行均分，每层土体根据试验控制干密度进行称量和夯实，考虑到会出现积水情况，制备好后土柱的高度为1800 mm；土柱制作完毕后在有机玻璃管外用胶带黏接刻度清晰的卷尺，用于读取浸润峰和积水的深度。试验的数据采集系统为EC-5土壤水分传感器和CR1000智能数据采集仪。EC-5土壤水分传感器其探针部分全部插入土柱，连接线路，接入CR1000智能数据采集仪，然后数据采集仪与电脑连接进行数据采集。人工降雨系统根据试验前设置的降雨强度，按试验前设定的时间间隔计算出向土柱上表面所加水的质量，并采用精度为0.01 g的电子天平进行称量后均匀缓慢加至土柱上表面。 3.2. 供试土样 试验所用淤积层淤积土经过0.1 cm的筛子过筛后，加水配制成初始质量含水率为7.2%的土样(初始体积含水率的理论计算值为11.5%)，然后击实为干密度为1.350 g/cm 3 的一维垂直土柱。通过室内常水头渗流试验，测得重塑淤积土在干密度为1.350 g/cm 3 下的饱和渗透系数均值为4.03 × 10 −6 m/s (约为0.24 mm/min)。 3.3. 试验方案 为模拟不同降雨强度下一维淤积土垂向入渗过程，分别设置了降雨强度为10、25、35和45 mm/h的4组模拟试验。当设置降雨强度为10 mm/h时，因为该时段降雨强度小于饱和渗透系数，每隔10 min向土柱上表面加一次水，理论计算可得每次的加水质量为85.48 g，累积降雨历时和降雨量分别为20 h和300 mm。当设置降雨强度为25、35、45 mm/h时，每隔5 min向土柱上表面加一次水，不同雨强下加水质量分别为64.11、85.48、128.22 g，降雨历时分别为20、15、10 h，累积降雨量分别为400、525和450 mm。EC-5土壤水分传感器插入土柱内，连接好CR1000智能数据采集仪，并设置每隔10 min采集一次数据。并将制备好的土柱，在玻璃管黏接刻度清晰的卷尺，并在试验过程中每次加水前读取浸润峰和积水的深度。
system: [SYSTEM]

user: 5. 结论 降雨垂直入渗补给过程是淤地坝淤积层水循环的重要过程，也是降雨径流挖潜的重要组成部分。通过自行设计加工的试验装置对初始体积含水率为11.50%的淤积土柱进行了4组不同雨强下的垂向入渗特性试验，主要得到以下结论： 1) 降雨垂向入渗结果表明，当降雨强度小于淤积层的最小入渗能力时，淤积层垂直入渗率时程曲线为一平行于横轴的直线；当降雨强度大于最小入渗能力时，垂向入渗率时程曲线呈无压入渗、有压入渗和饱和入渗三个阶段变化。 2) 不同降雨强度下积水点和饱和点出现的时间点不同，降雨强度越大时，积水点和饱和点出现的时间越短，而有压入渗经历的阶段越长。 3) 不同降雨条件下浸润峰深度时程变化呈现三个阶段变化趋势；且表现出平均传导率在降雨初期最大，随着降雨历时的)增加逐渐衰减，最后趋于稳定值。 4) 从土壤水分传感器测试的含水率分布可知，淤积土层表面的体积含水率不会迅速达到饱和状态，而是经历一个由非饱态和到暂态饱和态，最后达到饱和态的变化过程。
system: [SYSTEM]

user: 摘要 自2015年以来，P2P网络借贷平台发展出现滑铁卢，在短短几年间，市场上正常运营的P2P网络借贷平台迅速从3000多家锐减为300多家。一方面是因为外部需求在不断发生变化，平台传统的经营模式已不适应发展需求；另一方面平台自身的经营绩效也存在一定的问题。本文运用数据包络分析方法研究P2P网络借贷平台的DEA有效性，同时对比经过ICP认证与未经ICP认证的P2P网络借贷平台的DEA有效性差异，并对非DEA有效的P2P网络借贷平台进行DEA有效性指标数值调整。
system: [SYSTEM]

user: 1. 引言 P2P网络借贷起源于欧洲，但在我国发展较快。2007年，中国首个P2P网络借贷平台“拍拍贷”在上海成立。一些具有创业精神的P2P投资者开始尝试建立网上借贷平台。在随后的5年里，我国P2P网络借贷平台进入了起步阶段，拍拍贷、人人贷、红岭创投等知名的网络借贷平台也在这一阶段出现。2012年11月，国际互联网集团董事长兼首席执行官首次在第五届移动互联网博览会上推出“互联网+”概念，一时间各行各业都加入到“互联网+”的行列之中。在2012年到2015年，P2P网络借贷平台数量呈现爆发式增长，至2015年，正常运营的P2P网络借贷平台高达3576个，创下历史新高，相比2012年的涨幅超过25倍。但从2016年开始整个行业进入下坡，至2019年底，市场上正常运营的平台数量下降至343家。近10年P2P网络借贷平台数量变化如图1所示。 图1. 各年P2P网络借贷平台数量变化 根据国家《互联网管理办法规定》，如果经营性网站不办理ICP (Internet Content Provider)经营许可证即视为非法经营。平台运营的状态决定了平台上全体用户的效益，平台的绩效也直接影响到平台上的每一位直接投资人和潜在投资人。P2P网络借贷平台存在的目的是盈利，平台的经营绩效同时受到国家宏观调控和行业自身交易指标影响，因此平台的绩效研究对平台本身和整个行业的发展都尤为重要。经营绩效作为一个多方面指标共同作用的综合结果，能很好的反映出平台本身的经营能力，这也正是决定平台发展和保护投资人合法权益的首要因素。因此本文以2020年6月的P2P网络借贷平台为研究对象，通过前人的文献研究选取合适的输入输出指标，运用DEA方法对我国的P2P网络借贷平台的综合效率进行评定，对比分析经过ICP认证和未经过ICP认证的P2P网络借贷平台的DEA有效性差异，并以此为基础为非DEA有效平台提供指标性数值调整，最后根据研究结果向P2P网络借贷平台提出合理化的建议。
system: [SYSTEM]

user: 2. 文献综述 2.1. P2P网络借贷平台研究综述 自2007年国内首个P2P网络借贷平台成立以来，经过十余年的发展，P2P网络借贷平台已从最初的几家发展成如今的6000多家，而正常运营的P2P网络借贷平台数量不过300余家，国家对于P2P平台的监管也越来越严格。 近年来，越来越多的学者对P2P网络借贷平台进行了不同方面的研究，主要包括平台的运营绩效、经营风险、法律(监管)和对策等相关问题。郑思海等运用不同的研究方式分别从量化的角度、基于改进TOPSIS法、基于DEA模型和绩效影响因素等方面对P2P网络借贷平台的运行效率进行了不同层面的研究 [ 1 ] [ 2 ] [ 3 ] [ 4 ]。朱宗元 [ 5 ] 运用数据包络分析(DEA)和层次分析法(AHP)研究了22个P2P网络借贷平台的运营效率，研究结果表明，各平台的运行效率普遍较低，不同平台之间存在显著差异。王重润 [ 6 ] 运用DEA-Tobit方法对P2P平台的融资效率进行了测度，并研究了平台效率的影响因素。研究结果显示平台的注册资本、借贷款人数、运营时间、信用等级、正常还款次数和借款用途与平台的融资效率正相关，逾期还款次数和借款期限等因素与平台的融资效率负相关。 孟靖贻 [ 7 ] 将P2P网络借贷平台的运营模式分为四类：1) 机构担保交易模式；2) P2P平台下的债权合同转让模式；3) 大型金融集团创建的互联网服务平台；4) 以交易参数为基点，结合O2O的综合交易模式。并总结了现有P2P网络借贷平台面临的四大风险：信用风险、技术风险、市场风险和流动性风险、法律及政策环境风险。同时针对现有问题提出四个对策建议：1) 甄别和监管网络平台及借款人信用风险；2) 控制信贷规模；3) 培养专业人才，提高技术水平；4) 加强P2P网贷立法。 2.2. 数据包络分析研究综述 在实际应用中，数据包络分析(DEA， Data Envelopment Analysis)多用于对能源、资源和管理等效率的研究。Castañeda [ 8 ] 提出了一个基于数据包络分析的软件工厂及其项目相对效率评价模型，通过软件工厂不同工作单元的活动来衡量软件工厂软件生产构件的生产率。Salas-Velasco [ 9 ] 以西班牙公立大学为例，基于两个阶段的程序来评估他们在提供教学和研究方面的内部效率。在第一阶段，估计了每个大学的技术效率得分，西班牙大学的平均效率约为92%。在第二阶段，将效率得分与一组协变量进行回归，以研究它们与大学效率水平的关系。研究发现，受赠者比例较高的大学往往效率较低，而拥有终身教职的学者比例较高，会提高西班牙高等教育部门的生产效率。最后，计算了Spearman在SCImago和Shanghai大学排名中的DEA效率得分与西班牙大学分类之间的等级相关性，结果显示DEA分数给西班牙大学的排名与他们在公认的排名中的位置相匹配。现有的研究中也有不少学者对传统DEA方法进行改进以更好地解决实际问题。Feng Li [ 10 ] 提出了一个基于DEA方法的扩展生态低效率模型，以解决银行活动及其可能的相对绩效。由于生态效率是通过一组最优权值得到的，因此从数据驱动的角度考虑所有可行权值集，进一步构建了一个性能优势结构。然后，发展了两个成对生态无效率优势概念，并提出了无效率优势概率模型。最后，以2014~2018年32家中国上市银行为例，通过对生态无效率优势概率法的实证研究，验证了该方法的有效性和实用性。 国内对P2P网络借贷平台的经营效率研究多用DEA或DEA与其他方法的结合进行研究，但大多数学者只是判断出P2P平台经营效率的DEA有效性，然后根据研究结果提出理论建议，并未给出具体的数值调整使平台的经营效率达到DEA有效。本文拟用DEA中的CCR模型对P2P网络借贷平台进行运行效率的DEA有效性评价研究，并对非DEA有效的平台进行具体的指标性数值调整至DEA有效。
system: [SYSTEM]

user: 3. 模型介绍 1978年，美国著名运筹学家、克萨斯大学教授A. Charnes、W.W. Cooper和E. Rhodes提出数据包络分析。(DEA， Data Envelopment Analysis)。DEA是一种基于数学和运筹学的线性规划模型。它表示为产出与投入的比率。通过比较一个特定单元的效率与一组提供相同服务的相似单元的性能，DEA试图使单元的效率最大化。在这个过程中，一些效率达到100%的单元称为相对效率单元，而效率得分低于100%的单元称为无效单元。DEA作为一种效率评价方法，具有以下优点：1) 不涉及中间函数；2) 在操作过程中不涉及指标权重；3) 评价决策单元(DMU)的有效性相对客观。 CCR是由Cooper、Charnes和Rhodes三人提出并由三人姓名的首字母命名组成，是评价生产效率的一种非常重要的非参数方法。其主要原理是保持决策单元的投入或产出不变，通过数学规划确定相对有效的生产前沿，然后将每个决策单元投影到DEA生产前沿，比较决策单元偏离DEA前沿的程度，评价其相对有效性。 设有m个DMU，每个DMU有a种生产要素投入和b种生产要素产出，记输入指标 X j = ( x 1 j ， x 2 j ， ⋯ ， x a j ) T 和输出指标 Y j = ( y 1 j ， y 2 j ， ⋯ ， y b j ) T ， j = 1 ， 2 ， ⋯ ， m 。输入指标权重 v = ( v 1 ， v 2 ， ⋯ ， v a ) T ，输出指标权重 u = ( u 1 ， u 2 ， ⋯ ， u b ) T 。在研究过程中，为方便记忆，将被研究单位输入指标记为 X 0 = X j ，输出指标记为 Y 0 = Y j ， j = 1 ， 2 ， ⋯ ， m 。第j个决策单元的效率评价指数表示为： h j = u T y v T x = ∑ i = 1 b u i y i j ∑ i = 1 a v i x i j ， j = 1 ， 2 ， ⋯ ， m (1) 在公式(1)中，得到的h j 值越大则说明对应的决策单元DMU j 效率越高。若将第j个决策单元的效率指标作为评价目标，其余的决策单元作为约束条件，则可以得到CCR模型表示为： ( C 2 R ) T { max u T Y j v T X j = max ， u T Y j v T X j ≤ 1 ， j = 1 ， 2 ， ⋯ ， m ， u ≥ 0 ， v ≥ 0. (2) 将公式(2)的分式规划问题线性化，则第j个DMU j 的生产效率模型可表示为： ( P c ) { max u T Y j = Z v T X j − u T Y j ≥ 0 ， j = 1 ， 2 ， ⋯ ， m v T X j = 1 u ≥ 0 ， v ≥ 0 (3) 设 ε > 0 为非阿基米德无穷小量( ε 小于任何大于0的正数)，则其对偶模型为： ( D c ) { min [ θ − ε ( e ^ T S − + e T S + ) ] = z ∑ j = 1 m X j λ j + S − = θ X 0 ∑ j = 1 m Y j λ j − S + = Y 0 λ j ≥ 0 ， j = 1 ， 2 ， ⋯ ， m ; S + ≥ 0 ， S − ≥ 0 (4) 其中， e ^ T = ( 1 ， ⋯ ， 1 ) ∈ E m ， e T = ( 1 ， ⋯ ， 1 ) ∈ E s ， S + 为松弛变量， S − 为剩余变量，且 S + 和 S − 的值越小越有效。在实际应用中，利用该模型的最优值 θ 0 ， S − 0 ， S + 0 ， λ 0 来判断对应DMU 0 的DEA有效性： 1) 若 θ 0 = 1 ， S − 0 和 S + 0 存在非零值，则对应决策单元为弱DEA有效。若 S − 0 的子向量存在非零值，则说明该子向量对应的投入指标未得到充分利用；若 S + 0 的子向量存在非零值，则说明该子向量对应的产出指标离最优产出指标相比尚存不足。 2) 若 θ 0 = 1 ， S − 0 = 0 ， S + 0 = 0 ，则对应决策单元为DEA有效。此时对应的决策单元同时达到技术有效和规模有效。 3) 若 θ 0 < 1 ， S − 0 > 0 ， S + 0 > 0 ，则对应决策单元为非DEA有效。 θ 0 → 1 表明该决策单元有效性越大； S + ， S − 的值越大，则表明对应投入或产出离有效前沿面越远。现实原因可能是因为生产投入过大，或生产能力未能得到充分发挥，而造成产出不足。 根据对偶模型分析结果，还可以了解到决策单元的规模收益情况： 1) 若 1 θ 0 ⋅ ∑ j = 1 m λ j = 1 ，则决策单元规模收益不变，且已达到最佳。 2) 若 1 θ 0 ⋅ ∑ j = 1 m λ j < 1 ，则决策单元规模收益递增，此时增加投入可加快产出水平增长速度。 3) 若 1 θ 0 ⋅ ∑ j = 1 m λ j > 1 ，则决策单元规模收益递减，此时增加投入反而会降低产出的增长速度。 在CCR模型中， x 0 * = θ 0 x 0 − s − 0 ， y 0 * = y 0 + s + 0 表示DMU 0 在生产有效前沿面上的投影， x 0 * 表示DMU 0 调整后的投入量， y 0 * 表示DMU 0 调整后的输出量。
system: [SYSTEM]

user: 5. 结论与展望 从DEA有效性方面看，有超过半数的P2P网络借贷平台为非DEA有效，可见目前市场上的P2P网络借贷平台的运营绩效还存在一定不足，在未来还有待于进一步提高各大P2P网络借贷平台的投入产出效率。从规模收益方面看，在选取的17个样本P2P网络借贷平台中，DEA有效的8个P2P网络借贷平台的规模收益均处于不变状态；非DEA有效的9个P2P网络借贷平台中有8个平台的规模收益呈现递减状态，且此8个平台的“注册资金”和“平均利率”两个指标的DEA有效性调整值与平台目前的指标值相差较大，说明现有部分P2P网络借贷平台存在盲目投入现象，同时其平台项目“水分”较重，经营易出现裂口。在激烈的竞技场中，部分P2P网络借贷平台为赢得市场，不断增加投入，导致吸收的资本超出产品发展需求，造成资源浪费；同时为获取利润，疯狂贷款项目，不断寻找借款人，导致平台的项目质量愈加参差不齐，平台的规模收益呈现递减状态。从ICP认证与否方面看，总体来说经过ICP认证之后的平台的DEA有效性相较于未经过ICP认证平台的DEA有效性更高，且经过ICP认证的平台中DEA有效的平台所占比例也更高。 DEA模型使用的前提是指标之间具备相互独立性，本文的研究对出现的指标之间未进行相关性检测而直接使用DEA模型，研究结果缺少一定的说服力。希望未来的研究能在检测指标之间的关联性的基础上进行DEA有效性分析，以便得出更加准确有效的结论。
system: [SYSTEM]

user: 摘要 光谱成像技术因其高准确率、高分辨率、连续性、遥感性等多方面优势，使其在农业、矿物识别、植物生态学、食品安全、医学诊断、航天等领域有着无可替代的作用，并且因其广泛的应用前景激发了众多研究者的研究兴趣。本文阐述了光谱相机的原理及遥感技术用于水质检测现状，为探究光谱相机用于探测水体污染物的可行性以及探测结果质量，设计并进行了模拟实验。采集图像分析后发现模拟的污水在绿光波段反射率最高，在红光波段和蓝光波段反射率不是特别高，相比之下红光波段的反射优于蓝光波段，同时利用灰度均值的计算结果找到了最佳观测点。此研究可作为实际中光谱相机监测水污染的一种参考。
system: [SYSTEM]

user: 1. 引言 多光谱成像技术是应用于很多不同的研究领域的一种新兴技术 [ 1 ]，多光谱成像技术可以得到目标各像素点的高分辨率光谱。 水的质量分类指标其中包括：色度、气味、浑浊度、PH值、COD (化学需氧量)、BOD (生化需氧量)等等 [ 2 ]。现阶段被广泛使用的水质监测技术有原子光谱技术、生物传感技术、化学分析技术、电化学分析技术、色谱分离技术和分子光谱技术等 [ 3 ] [ 4 ]。由于操作复杂、检测时间长、成本高、寿命短等不足，目前被使用最多的方法是分子光谱技术 [ 5 ]。目前水质监测技术的方法大多是利用传感器投入水中或者取样后进行检测，对面积较大的水域进行监测时，有工作量较大、周期长等缺点，很难做到宏观、准确地监测愈发严重的水体污染环境。遥感技术是通过传感器采集物体辐射、反射的光波，然后进行处理得到影像资料，以此对物体镜像进行探测、识别 [ 6 ]，这样就可以实时、迅速的得到水域环境的质量情况，将其应用于水体污染物的检测中可以发挥其测量精度高、测量范围广等优势，同时我们考虑将光谱相机获取的水样图像通过一种全新的分析方法BRDF (双向反射分布函数法)进行分析，由此来评测水样被污染的情况。遥感技术从诞生开始就获得了极高的关注度，在对地观测、空间探测方面研究人员想要同时得到目标物的光谱信息和影像信息，这就导致了光谱技术与成像技术互相融合，从而衍生出一种全新的技术即光谱成像技术 [ 7 ] [ 8 ] [ 9 ]。近些年由于无人机的兴起与蓬勃发展，多光谱技术和高光谱技术有了新的需求与发展，现如今人们正在将多光谱技术、高光谱技术应用于无人机上进行低空作业，包括农作物监测、水污染监测、土壤检测等等 [ 10 ] [ 11 ] [ 12 ]。
system: [SYSTEM]

user: 2. 光谱相机及光谱分析的技术原理 2.1. 光谱相机 20世纪80年代，第一台成像光谱仪在美国被研制出来。光谱相机拥有多光谱相机、高光谱相机、超光谱相机等多种类型，其中多光谱相机可以在红外、近红外、紫外光等范围进行拍摄，利用滤光片、分光镜等设备对光线进行过滤、分解，可以在感光胶片上形成目标在不同波长的光波下的光谱图片，而且图片分辨率较高。光谱相机的基本结构由图像采集系统、光学系统、控制与信息处理器、探测器等装置组成。 2.2. 双向反射分布函数 BRDF最早是在光辐射角度进行定义，是用来描述光线如何在物体表面进行反射，表现其空间反射分布特性和光谱特性，其理论认为一个表面被一个方向光照射后会发生反射，但是在表面上方整个半球的每个方向都会有反射光出射，各个方向上反射的能量是不同的，BRDF表示的某个确定反射方向的光与指定方向入射光的关系： f ( l ， v ) = d L 0 ( v ) d E l (2-1) 其中f就是BRDF函数，l、v分别为入射光和反射光的方向向量。 d L 0 ( v ) 反射光方向上的微分辐射率， L 0 ( v ) 是指表面上全部入射光对v这一个反射方向的辐射率，微分辐射率是指仅从l这一个方向入射的光对v方向的反射辐射率， d E ( l ) 是指仅从v方向上接收到l方向入射光的辐照度，E为v方向上接收到所有入射光的辐照度，而 d E ( l ) 同样是指特定为l入射方向的辐照度。 辐射功率是指单位时间内单位面积上的总辐射能，用Φ表示，单位为J/s，定义公式为： Φ = d Q d t (2-2) E表示的是辐射照度是指受照面单位面积上的辐射通量，定义如下： E = d Φ d A (2-3) L为辐射率，指的是某单位面积内的单位立体角的辐射通量密度，其中 d A ⊥ 是面积微分dA在与光线垂直方向上的投影，如下图1。 L = d Φ d ω d A ⊥ (2-4) 图1. 微分面积投影图 BRDF也可以用表面上半球建立球坐标系来表示，如下图2。θ为极角是指入射或反射光线与法线的夹角，同时在平面上建立一个坐标轴，入射反射光线在平面上投影与该坐标轴的夹角为ϕ，称为方位角。由此BRDF也可以用极角和方位角表示： f ( θ i ， ϕ i ， θ o ， ϕ o ) ，对于各向同性的物质，可以用两光线投影的夹角ϕ来代替 ϕ i 和 ϕ o ，即 f ( θ i ， θ o ， ϕ ) 。 图2. 球坐标系