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

user: 2. 探针探测数据集解析 系统中，每个探针周期性的将探测到的WIFI终端信息发送到数据库服务器，数据库服务器中，探针探测到的一个WIFI终端的信息被保存到表格TZD中。表格TZD定义了ID、tID、tData、tTime四个字段，其中：ID字段为增量为1的标识字段，用于标识每一条被探测到的WIFI终端的信息；tID为探测到该条WIFI终端的信息的探针的名字；tData为探针探测到的WIFI终端相关信息；tTime为该条探测信息被保存到数据库的时间。依据探针数据协议中规定的格式，需要将存储在tData字段中的内容解析为WIFI终端MAC地址、AP点MAC地址、信号大类类别、信号小类类别、传输信道、信号强度等信息。最终，表格TZD的每一行被保存为表格Analyze_LocTZdata的一行。数据集解析的流程如图1示意。 表格Analyze_LocTZdata的字段包括timeFrameCounter、ID、tID、tData、tTime、WIFI终端MAC地址、AP点MAC地址、信号大类类别、信号小类类别、传输信道、信号强度、解析时间、终端位置locClass、终端水平位置locX、终端水平位置locY。其中，ID、tID、tData、tTime字段的取值为表格TZD对应字段的取值；WIFI终端MAC地址、AP点MAC地址、信号大类类别、信号小类类别、传输信道、信号强 图1. 数据解析流程 度等字段的取值从表格TZD中tData字段取值中解析获得；解析时间字段的取值为解析tData字段完成后的时间；终端位置locClass取值为WIFI终端所在区域的标识，取值为0时表示在区域外，1表示为在区域内；终端水平位置locX、终端水平位置locY字段的取值为WIFI终端所在位置的横坐标和纵坐标；timeFrameCounter为数据采集时间tTime的时间帧编号。 探针探测数据存储表格TZD中的数据被解析到Analyze_LocTZdata表格时，timeFrameCounter、终端位置locClass、终端水平位置locX、终端水平位置locY取值均为空。其中，终端位置locClass、终端水平位置locX、终端垂直位置locY的取值需通过人工标注，而timeFrameCounter由专门的编码过程赋值。
system: [SYSTEM]

user: 3. 面向室内人员定位的探针数据构造 考虑到WIFI终端发送probe帧的随机性和突发性，本文相关的研究在使用WIFI探针计量建筑物内指定区域内人员数量时将一段时间内被WIFI探针感知到的全部WIFI终端信息作为一个时间帧内被探针系统感知到的WIFI终端信息。全部时间帧自1开始顺序编号。因此，对表格Analyze_LocTZdata，需要根据每行tTime的取值确定timeFrameCounter的值。图2给出了对Analyze_LocTZdata的每一行的timeFrameCounter字段的赋值流程。图2示意的流程中，起始时间timeStart、时间帧长度timeInterval是 图2. 探针数据的时间帧计数流程 必须的输入；timeFrameCounter的值逐1增加，时间格式为YYYY-MM-DD hh:mm:ss，精确到毫秒。 为准确的判定WIFI探针感知到的WIFI终端是否在指定区域内，研究在指定区域内部署了4个探针，只有当4个探针在同一时间帧内都感知到某一WIFI终端时才对该终端是否在指定区域内进行判定。所以需要从解析后的WIFI探针数据集中提取具timeFrameCounter取值相同、被4个探针分别感知的WIFI终端的数据。虽然分析的目的以及所使用的分析模型不进相同，但目标数据的基本格式可以表述为四元组rssiV = (RSSI 1 ， RSSI 2 ， RSSI 3 ， RSSI 4 )，其中RSSIi表示第i个探针感知到WIFI终端的probe帧时测得的接收信号的强度指示(Received Signal Strength Indication)值。为适合判断是否在指定区域内的室内定位模型构建，研究使用的数据是四元组rssiV形式基本格式的扩展，rssiV’ = (RSSI1， RSSI2， RSSI3， RSSI4， locStatus)，其中，locClass为WIFI终端所在区域的标识，取值为0时表示在区域外，1表示为在区域内。从解析后的数据表格Analyze_LocTZdata中构造目标数据集的流程由图3给出。 图3. 探针定位模型数据集构造流程 图3中：所谓构造原始矩阵所需的列是指从Analyze_LocTZdata中构造具有timeFrameCounter、WIFI终端MAC地址、信号强度、终端位置locClass等列的子视图；而最终数据集矩阵Dataset是指为便于不同研究平台下数据共享而以excel文件格式保存好的目标数据。
system: [SYSTEM]

user: 4. 实验结果与分析 为快速验证基于上述探针探测数据集解析、面向室内人员定位的探针数据构造等步骤以及流程所构造的数据集的有效性，研究在Matlab环境中使用BP神经网络快速实现了一个辨识WIFI终端是否在指定区域判别器：1) 该判别器以四个WIFI探针在同一时间帧内嗅探到的同一WIFI终端的RSSI值为输入，每个WIFI探针对应着一个输入神经元，即该判别器使用的BP神经网络有四个输入神经元；2) 该判别器使用的BP神经网络的输出层包含一个输出神经元，该神经元输出值为1时表明被感知的WIFI终端在指定的区域内，而输出值为0时表明被感知的WIFI终端在指定的区域外；3) 该判别器使用的BP神经网络有1个隐层，包含5个隐含神经元；4) 该判别器使用的BP神经网络隐层神经元的激活函数为 f ( x ) = tan s i g ( x ) = 2 1 + e − 2 x − 1 ，输出层神经元的激活函数为 f ( x ) = purelin ( x ) = x 。试验中，对训练完成 的判别器，若其使用的BP神经网络的输出值x时，若 x ≥ ε ，则认为判别器的输出为1，否则认为判别器的输出为1，其中ε为实现指定的阈值。 实验使用的数据为我们部署在智能建筑与建筑节能安徽省重点实验室物联网应用研究室的WIFI探针系统采集到的WIFI终端数据；关注的WIFI终端共16部，7部放在室内的不同位置，9部放在室外走廊的不同位置。 实验中，全部数据分成样本数据和测试数据两部分。实验中，对每一个采样率，判别器被重复1000次构造并计算判别的准确率。图4显示了当采样率从0.005到0.995变化时，使用基于BP神经网络的WIFI终端是否在指定区域判别器判别的准确率随着采样率的变化而变化的情况。图4中，横坐标表示采样率，纵坐标表示识别的准确率。图4中，蓝色实线表示对样本数据进行判别时判别准确率的平均值随 图4. 采样率对识别准确率的影响 着采样率的变化情况，红色虚线表示对测试数据进行判别时判别准确率的平均值随着采样率的变化情况。显然，在采样率不需要很大时，基于样本数据，判别器可以以较高的准确率(90%)辨识WIFI终端是否在指定的区域中。类似的结论对测试数据依然成立。这表明，依据多探针同时感知到的同一WIFI终端的RSSI值判别WIFI终端的位置是可行的。 表1出来了采样率为0.1、0.2、∙∙∙、0.9时，使用判别器判别WIFI终端是否在指定区域内的误判情况。表1中，错误率1表示被样本数据的判别错误情况，错误率2表示被测试数据的判别错误情况。显然，对不同的采样率，使用基于BP神经网络的判别器辨识WIFI终端是否在指定区域内不仅保持较低的错误率，且错误率的变化较稳定，这进一步表明，依据多探针同时感知到的同一WIFI终端的RSSI值判别WIFI终端的位置是可行的。 从表1中可以发现，虽然在采用率比较小的时候无论是错误率1还是错误率2的平均值都保持较低的水平，但是由于方差都比较大(与错误率在同一数量级)。从图5和图6可以更直观的发现同样的现象： Table 1 Sample rate VS identification error rat 采样率 错误率1 错误率2 平均值 方差 平均值 方差 0.1 0.0409 0.0243 0.0785 0.0162 0.2 0.0565 0.0179 0.0715 0.0135 0.3 0.0599 0.0173 0.0672 0.0154 0.4 0.0562 0.0086 0.0651 0.0065 0.5 0.0510 0.0067 0.0658 0.0053 0.6 0.0515 0.0067 0.0677 0.0059 0.7 0.0532 0.0054 0.0645 0.0058 0.8 0.0555 0.0048 0.0594 0.0071 0.9 0.0566 0.0049 0.0469 0.0096 表1. 采样率与识别错误率 图5. 采样率对样本数据识别错误率(error1)的影响(带标准差) 图6. 采样率对样本数据识别错误率(error2)的影响(带标准差) 图5和图6分别以及给出的是采样率对样本数据识别错误率(error1)的影响采样率对测试数据识别错误率(error2)，显然，无论从图5还是图6都可以发现，当采样率较小时，判别器识别的误差率波动较大。因此故实际应用中需要选取合适的采样率，在辨别器的训练时间、错误率以及错误率的稳定性三个性能指标之间进行平衡。
system: [SYSTEM]

user: 5. 结论与研究展望 围绕WIFI探针系统采集到的WIFI终端数据的预处理，本文将WIFI终端数据预处理流程划分为探针探测数据集解析、探针探测数据时间帧编号、面向室内人员定位的探针数据构造等三个阶段。并在对每个阶段的任务目标厘清的基础上设计了相关预处理任务的流程并进行了实现。实验结果表明，所获得的预处理后的WIFI终端数据被当作基于BP神经网络的判别器的输入时，该判别器可以以很高的准确率判别WIFI终端是否在指定的区域内。基于BP神经网络的判别器的性能表明依据多探针同时感知到的同一WIFI终端的RSSI值辨识WIFI终端是否在指定的区域内是可行的。 为更好的提升基于BP神经网络的判别器在辨识WIFI终端是否在指定的区域内时的准备率，基于BP神经网络的判别器的优化，包括BP神经网络的隐层数、隐层神经元的数量、判别器阈值都是值得优化的参数。进一步，由于BP神经网络中要求所有非输入神经元的激活函数一阶可导，而是否在指定区域内的0/1辨识结果是离散的的，因此，结合具体的应用需求，面向离散判别的基于神经网络、辨识WIFI终端是否在指定区域内的判别器的前向神经网络架构以及训练方法是值得探索的楼宇智能化实现相关的研究。
system: [SYSTEM]

user: 摘 要 调水工程是解决地区水资源不足或时空分布不均衡问题的必要工程措施，不同水库间通过管道和河流调水会对调蓄水库水温产生影响。采用二维水流水温数学模型模拟了吉林省引松供水工程调水对调节水库石头口门水库水温的影响。计算结果表明，调水对调蓄水库水温影响的范围、时间主要和调水量、调水温差、水库水深、入水口位置等有密切关系。引松供水工程调水量和调蓄水库的水量相比很小，调水对水库水温分布的影响有限。调水入水口位于库尾，受影响的范围基本局限于水库库尾局部范围内，水温影响在1个月内即达到平衡，对整个库区和坝前影响很小。调水工程入水口远离坝址可以减小调水对调蓄水库及下游河流的影响程度。 关键词 :调水工程，调蓄水库，水温影响，数学模型 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 我国许多地区由于特定的地理条件，水资源时空分布很不均衡，导致地区工程性缺水、资源性缺水或水质性缺水。往往需要通过调水工程来有效利用和调配水资源。如南水北调、引江济淮、引黄入晋工程等都属于大型的跨流域调水工程。在调水工程中，必须具有足够水量的调蓄水库来保证调度运行的灵活性和便捷性。调水工程的实施必然会对调蓄水库水环境及其下游河道产生影响，许多学者都对调水工程的生态环境影响进行了研究 [ 1 ] [ 2 ] [ 3 ] [ 4 ] 。 随着工农业用水对水质要求的提高，调蓄水库除了调节水资源量外，还要兼顾调水对灌溉、渔业的影响，这些都和水库水温关系密切。一旦调水水量过大，调水水温过高或过低达不到生态用水温度要求，将会引起灌溉农作物大幅度减产，同时对下游河道中的水生生物特别是鱼类产生毁灭性打击。如刘建康、蒋固政的研究表明，春夏季节水库下泄水温低于建库前天然河道水温，将直接导致下游鱼类繁殖季节推迟、当年幼鱼生长期缩短、生长速度减缓 [ 5 ] [ 6 ] 。韩彩霞、吴佳鹏研究了水库下泄水温对农作物的影响，如果下泄水温未达到最适宜温度范围，将使得灌溉难以发挥最佳效益 [ 7 ] [ 8 ] 。因此，准确预测调水对调蓄水库水温分布的影响十分必要。 20世纪60年代初美国最先开始研究水库水温的分布和变化规律。此后，许多学者采用不同种方法和手段对水库库区水温分布进行大量研究工作，取得了许多有意义的成果。如刘平等 [ 9 ] 通过经验公式法对迪那河五一水库水温及下泄水温进行了预测分析。脱友才、杨红秀分别对丰满水库和汾河水库不同季节的水温变化分布规律进行了原型观测和资料分析 [ 10 ] [ 11 ] 。赵一慧、甘衍军等 [ 12 ] [ 13 ] 通过EFDC数学模型模拟了小浪底、二滩水库的水温分布变化规律。江春波、戴凌全等 [ 14 ] [ 15 ] 对水库水温数学模型研究进展进行了综述，指出了不同模型的优缺点以及需要进一步研究的方向。李林等 [ 16 ] 通过三维水温模型研究了支流影响下的水库水温变化特征，表明支流入流越靠近坝址对坝前水温分布影响越大。高学平、柳海涛等 [ 17 ] [ 18 ] 通过水温模型试验模拟了糯扎渡水库等大型水库的水温分布规律。 本文通过水流水温数学模型，研究了吉林省引松供水二期工程调水对石头口门调蓄水库水温的影响，为工程实施提供技术支撑，对类似调水工程设计具有借鉴意义。
system: [SYSTEM]

user: 2. 工程概况 吉林省中部城市引松供水工程位于吉林省中部，主要工程措施是从第二松花江丰满水库库区引水至吉林省中部地区，增加向长春、四平、辽源等城市群的生产生活供水，可为退还被城市用水挤占的农业用水和河道生态环境用水、限采地下水创造条件，为中部城市群的振兴和发展提供水资源保证，并改变部分地区的饮用水水质，对保障吉林省百亿斤粮食生产和饮水安全、促进地区乃至全省经济社会可持续发展、改善生态环境具有重要作用。 依据总体规划，考虑工程的实施进度安排，现状基准年采用2013年，近期水平年为2020年，设计水平年为2030年。设计水平年多年平均引水量为8.98亿m 3 ，丰满水库取水口设计引水流量为38.0 m 3 /s。目前吉林引松工程一期工程顺利实施，二期工程正在方案论证阶段，共包括12条输水支线和1座调蓄水库为石头口门水库，通过在总干线设立饮马河分水口向石头口门水库调水。 饮马河是第二松花江较大的一条支流，发源于磐石市驿马乡呼兰岭，流经磐石、双阳、永吉、九台、德惠、农安六县(市)，至农安县靠山屯与伊通河汇合，北流近20公里汇入第二松花江。饮马河分水口距离石头口门水库约40 km，调水所经河道弯曲，河槽较浅，河底为细沙和淤泥。汛期流量较大，非汛期流量较小。 石头口门水库供水对象为长春市玉米工业园、九台市、德惠市，水库控制流域面积 4944 km 2 ，原设计是以防洪除涝、农业灌溉为主，结合发电和养鱼等综合利用的大型水库。兴利增容工程完成后，水库的任务调整为防洪、城市供水、灌溉为主，兼顾发电。正常蓄水位为189.0 m，死水位为182.5 m，兴利库容38，409万m 3 。图1为饮马河流域及石头口门水库示意图。 图1. 饮马河流域及石头口门水库示意图
system: [SYSTEM]

user: 3. 数学模型的建立 3.1. 水流基本方程 1) 连续方程 ∂ Z ∂ t + ∂ H u ∂ x ∂ H v ∂ y = q (1) 2) 动量方程 ∂ H u ∂ t + ∂ H u u ∂ x + ∂ H u v ∂ y = − g H ∂ Z ∂ x − ∂ H u ′ u ′ ¯ ∂ x − ∂ H u ′ v ′ ¯ ∂ y + f H v + τ s x − τ b x ρ + q u * (2) ∂ H v ∂ t + ∂ H u v ∂ x + ∂ H v v ∂ y = − g H ∂ Z ∂ y − ∂ H v ′ u ′ ¯ ∂ x − ∂ H v ′ v ′ ¯ ∂ y − f H u + τ s y − τ b y ρ + q v * (3) 式中： u 、 v 分别为 x 、 y 方向流速分量； Z 为水位； H 为水深； f = 2 ω sin φ 为科氏力系数(ω为地球自转角速度)； q 为单位面积上的源汇强度(m 3 /s∙m 2 )， u * 和 v * 为源汇节点周边 x 、 y 方向流速分量。 ( τ s x ， τ s y ) 为表面风应力项， ( τ b x ， τ b y ) 为底部床面阻力项。 3.2. 二维温水扩散运动方程 ∂ Δ T H ∂ t + ∂ ( H u Δ T ) ∂ x + ∂ ( H v Δ T ) ∂ y = D H x ∂ 2 Δ T ∂ x 2 + D y H ∂ 2 Δ T ∂ y 2 − k Δ T ρ C P + q Δ T i (4) 式中： Δ T 为水体的温升(℃)； Δ T i 为源项水体超温值(℃)； D x 、 D y 为热扩散系数； k 为水面综合散热系数； C P 为水的比热； ρ 为水体密度(kg/m 3 )； q 为源项单位面积流量(m 3 /s∙m 2 )。 3.3. 求解方法及定解条件 1) 求解方法 采用《海岸与河口潮流泥沙模拟技术规程》(JTS/T 231-2-2010)推荐的三角形单元法进行方程的离散和求解，采用Fortran计算机语言进行实际编程计算。 2) 初始条件 初始条件的影响将会在经过几次迭代后逐渐消失，因此本文采用冷启动来处理。即给定初始条件时刻计算域内所有求解变量为定值。 3) 边界条件 开边界水动力计算给定潮位过程。 温升计算边界出流时给定： ∂ ( D Δ T ) ∂ t + 1 C ξ C η [ ∂ ∂ ξ ( D Δ T u C η ) + ∂ ∂ η ( D Δ T v C ξ ) ] = 0 (5) 温升计算边界入流时： Δ T | Γ 0 = Δ T a ( t ， ξ ， η ) (6) 式中： Δ T a 为入流外海边界热回归温升过程，当选取范围足够大，即开边界在温排水扩散范围外时，可认为入流温升为零。 闭边界水动力计算时，流速满足不可入边界条件，即其法向流速为0。温升计算时闭边界采用绝热条件，即温排水扩散法向浓度梯度为0。 4) 动边界 工程水域有大片浅滩，高潮位时淹没，低潮位时出露，计算过程中采用冻结法处理动边界。 3.4. 计算条件 根据调水方案和饮马河基本概况，数学模型计算了1月份、5月份和7月份调水情况下的水库温度变化。其中，1月份调水流量为6.31 m 3 /s，调水温度为1.7℃，水库温度为0.3℃。5月份调水流量为7.91 m 3 /s，调水水温为7.5℃，水库水温为12.95℃。7月份调水流量为3.78 m 3 /s，调水温度为11℃，水库温度为26.5℃。 3.5. 计算参数 根据工程水域实测水温和风况统计，确定计算所需的平均水温和风力条件。在数学模型计算中，影响温度扩散的主要因素是热扩散系数和水面综合散热系数。热扩散系数的取值和温排水的扩散和输移与水流运动有关，一般取热扩散系数为常值，也有学者按沿水流方向和垂直水流方向分别取值。根据已有研究成果，本次计算采用下式计算：取D x = 0.5 m 2 /s和D y = 0.5 m 2 /s。水面蒸发系数和水面综合散热系数采用《工业循环水冷却设计规范(GBT 50102-2014)》中公式计算确定。图2为水库计算范围和计算网格。
system: [SYSTEM]

user: 4. 调水影响分析 为了分析调水对水库水温的影响，分别在水库库尾、水库中央和水库坝前(1#、2#和3#点)布置监测点，进行数据提取和分析，监测点位置见图3。 4.1. 监测点水温变化 表1为不同月份监测点水温随时间的变化特征值，图4为1月份监测点水温变化过程。可知，三个月份不同的调水量、调水温度和水库水温情况下，调水只对水库入口局部区域内水温产生影响，水库中部和坝前水温基本不受影响，水温影响在1个月内即达到平衡。这主要是因为水库调水量相对于水库水量非常小，不足以影响整个水库水温。 4.2. 调水水温影响范围 图5为1月份调水后水库在30天内水温的变化分布情况，可以看出由于调水口在库尾，调水量相对于整个 图2. 水库计算网格图 图3. 水库观测点布置图 图4. 1月份监测点水温变化 图5. 1月份调水后水库30天内水温变化分布 Table 1 Water temperature change at observation poin 月份 监测点 水温(℃) 变化最大值 0 d 5 d 10 d 15 d 20 d 25 d 30 d 1月 1 0.3 0.79 1.64 1.51 1.46 1.43 1.42 1.21 2 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0 3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0 5月 1 12.90 11.65 7.55 7.5 7.5 7.5 7.5 5.4 2 12.90 12.90 12.90 12.90 12.90 12.90 12.90 0 3 12.90 12.90 12.90 12.90 12.90 12.90 12.90 0 7月 1 26.5 26.49 20.49 14.36 14.82 15.44 16.63 12.14 2 26.5 26.5 26.5 26.5 26.5 26.5 26.5 0 3 26.5 26.5 26.5 26.5 26.5 26.5 26.5 0 表1. 水库计算监测点水温变化表 Table 2 Influence distance of water temperature chang 月份 影响距离(m) 5 d 10 d 15 d 20 d 25 d 30 d 1月 724 1112 1190 1480 1742 1945 5月 711 963 1187 1409 1411 1431 7月 675 857 1119 1189 1433 1452 表2. 石头口门水库水温影响距离 调节水库的水流非常小，调水对水库水温的影响范围仅限于库尾局部区域。 4.3. 调水水温影响距离 表2为调水对石头口门水库水温的纵向距离，其中1月份调入水的影响距离最远为1945 m，5月份调入水的影响距离最远为1431 m，7月份调入水的影响距离最远为1452 m。不同调水方式下，调水对水库水温的影响主要集中在库尾局部范围内。
system: [SYSTEM]

user: 5. 结论 本文针对吉林省引松供水二期工程，建立了石头口门水库平面二维水流水温数学模型，预测了水库调水后库区水温的变化情况、影响范围和影响时间等。主要结论如下： 1) 调水对水库水温影响的范围、时间主要和调水量、调水温差、水库水深、入水口位置等有密切关系。 2) 由于调水量和水库水量相比并不大，调水对水库水温分布的影响有限，受影响的范围基本局限于库末局部范围内，水温影响在1个月内即达到平衡。 3) 调水工程入水口位置尽量设置在水库库尾，远离坝址可以防止调水对水库坝前以及下游河流产生影响。 4) 平面二维水流水温数学模型可以对水库调水影响进行合理的模拟预报。但水库水深相对较大时，水温和水库水体的掺混具有典型的三维特征，更精细的预报需要借助于三维数学模型来模拟。
system: [SYSTEM]

user: 摘 要 本论文以实验室自行筛选并鉴定的一株撕裂蜡孔菌为材料，以多糖得率为评价指标，在单因素试验的基础上，利用响应曲面法的中心组合设计，对菌丝体多糖的微波辅助提取工艺进行优化，并对多糖的体外抗氧化活性进行了研究。结果表明：以水为浸提溶剂，微波辅助提取多糖的最优条件为微波时间30 s、微波功率480 W、液料比1:20、浸提时间30 min、浸提温度35℃。在此条件下菌丝体多糖得率预测值为6.523%，实测值为6.437%。体外抗氧化试验结果证实菌丝体多糖具有一定的抗氧化活性，抗氧化能力与浓度间有很好的量效关系，浓度越大，抗氧化能力越强。其中菌丝体多糖的抗脂质过氧化能力和清除超氧阴离子自由基能力较强，对羟自由基和过氧化氢的清除效果较好，但还原能力较弱。撕裂蜡孔菌菌丝体多糖具有天然的抗氧化能力，在开发多糖类天然抗氧化剂用品方面具有一定的现实意义。 关键词 :撕裂蜡孔菌，多糖，响应曲面法，微波，抗氧化活性 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 多糖(Polysaccharide)是一类由单糖脱水缩合而成的天然高分子化合物，分子量从几千到几百万不等 [ 1 ] ，是生物体内重要的能量储备和结构物质，也是除蛋白质和核酸以外的重要生命物质 [ 2 ] [ 3 ] ，功能涉及多细胞生命体的生长发育、成熟衰老及免疫感染等多个生理和病理过程 [ 4 ] 。大量的药理学实验表明，糖类化合物具有增强和调节机体免疫、抗病毒、抗肿瘤、抗辐射、抗衰老、抗氧化等作用 [ 5 ] [ 6 ] ，现已成为营养保健和功能性食品的重要来源。 多糖种类繁多，广泛存在于植物、动物、微生物(细菌和真菌)和藻类组织中。微生物多糖是多糖研究中较为详细的一类多糖，目前国内外学者从真菌(主要是担子菌纲和子囊菌纲)中已得到真菌多糖达数百种，香菇多糖、云芝多糖、灰树花多糖、裂褶多糖和槐树多糖等已经临床用于免疫、抗肿瘤治疗 [ 4 ] ；同时真菌多糖纯天然的抗氧化功能也成为科研工作者的研究热点 [ 7 ] 。 撕裂蜡孔菌(Ceriporia lacerata)，属于真菌界、担子菌门、多孔菌目、蜡孔菌属的一种真菌，多孔菌如灵芝、云芝等很多属于传统中药或民间药物。撕裂蜡孔菌是近几年新发现的蜡孔菌属的一个新品种，目前对其研究的还很少，但是专家估测其有很高的利用价值，研究可扩大到生物制浆造纸、饲料蛋白、抗癌药物和环境净化等领域。 微波辅助提取技术将溶剂浸提与微波提取相结合，可使样品瞬时产生大量热能，细胞破裂，细胞内物质易渗出，缩短实验时间，提高物质提取率 [ 8 ] 。响应曲面法(Response Surface Methodology， RSM)因其合理的试验设计，经济的试验方式，快速有效地确定多因素系统最佳条件等优势被广泛应用于众多过程的优化控制 [ 9 ] 。本论文以实验室自行筛选和鉴定的一株撕裂蜡孔菌(Ceriporia lacerata)菌丝体为原料，以水为提取溶剂，采用单因素试验对微波提取菌丝体多糖进行研究，并运用响应曲面法对影响多糖提取的几个因素进行了分析，优化其提取工艺；同时对菌丝体多糖的抗氧化性能进行分析，为撕裂蜡孔菌的有效开发利用提供理论依据。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 试验菌株和主要仪器 撕裂蜡孔菌(Ceriporia lacerata)：由江苏科技大学生物技术学院食品生物技术实验室自行从桑枝上分离和鉴定(NCBI登录号为：KP677607)。 主要仪器：Beckman ALLEGRA 64R高速冷冻离心机，DHZ-052D双层恒温摇床，UV9600紫外可见分光光度计，SW-CJ-1F型超净工作台，MCR-3微波化学反应器，日本三洋SANYO高压蒸气灭菌器。 2.2. 菌株发酵与菌丝体收集 将撕裂蜡孔菌接种在灭菌PD培养基中(去皮马铃薯200 g，葡萄糖20 g，水1000 mL，pH自然)，28℃培养箱(120 r/min)中震荡培养8 d。四层纱布过滤，将收集的菌丝体在不高于50℃烘箱中烘干，研磨后过80目筛，粉末收集备用。 2.3. 单因素及响应曲面优化试验设计 选取对菌丝体多糖提取有显著影响的料液比(mg/mL，1:10、1:20、1:30、1:40、1:50和1:60)、微波时间(s，20、30、40、50和60)、微波强度(W，160、320、480、640和800)、浸提时间(min，20、30、40、50和60)和浸提温度(℃，15、25、35、45和55)为考察因素，以多糖得率(%)为评价指标，进行单因素试验。在单因素试验结果的基础上，应用Design Expert 7.0软件进行试验设计，响应曲面法优化菌丝体多糖的提取工艺。 2.4. 多糖得率的测定 多糖含量的测定采用硫酸–苯酚法 [ 10 ] 。将多糖提取液离心(5000 r/min，10 min，25℃)，上清液用蒸馏水稀释后进行多糖测定。取2 mL待测液，加入1 mL苯酚溶液(6%，临时配制)，摇匀后，缓慢加入5 mL浓硫酸(加入时悬空，不断反复摇动试管，使反应充分)，摇匀，室温静置30 min，在490 nm波长下测定吸光度值。以葡萄糖为标准品绘制标准曲线。 多 糖 得 率 = 提 取 液 中 多 糖 质 量 ( g ) / 菌 丝 体 原 料 质 量 ( g ) × 100 % 2.5. 菌丝体粗多糖的体外抗氧化活性 将菌丝体多糖提取液离心(5000 r/min，10 min，25℃)，上清液经旋转蒸发仪浓缩，冷冻干燥。用蒸馏水配成0.2 mg/mL、0.4 mg/mL、0.6 mg/mL、0.8 mg/mL和1 mg/mL系列浓度梯度用于体外抗氧化活性的测定。 2.5. 1 . 菌丝体多糖的羟自由基清除能力 取菌丝体多糖溶液2 mL于试管中，依次加入2 mL 6 mmol/L的FeSO 4 溶液和2 mL 6 mmol/L的H 2 O 2 溶液，混匀后室温静置10 min。混合液中再加入2 mL 6 mmol/L的水杨酸，混匀后室温静置30 min。510 nm波长处吸光值记为A 1 ，双蒸水代替水杨酸时吸光值记为A 2 ，以双蒸水代替多糖溶液时吸光值记为A 0 [ 11 ] 。羟基自由基清除率按照公式E/% = [1 − (A 1 − A 2 )/A 0 ] × 100%计算。 2.5.2 . 菌丝体多糖的过氧化氢清除能力 取菌丝体多糖溶液5 mL于试管中，加入由50 mmol/L磷酸缓冲液(pH 7.4)配制的10 mmol/L H 2 O 2 溶液5 mL，缓慢摇匀，10 min后记下在230 nm处时的吸光度值，记为A 1 。A 0 为不加样品液的H 2 O 2 溶液吸光度值，A 2 为不加H 2 O 2 溶液的样品液吸光值。过氧化氢的清除率按照公式E/% = [A 0 − (A 1 − A 2 )]/A 0 × 100%计算 [ 12 ] 。 2.5.3 . 菌丝体多糖的还原能力 取菌丝体多糖溶液1 mL，加入2.5 mL磷酸缓冲液(pH 6.6)和2.5 mL 1%铁氰化钾溶液，混匀后50℃恒温水浴20 min，再加入2.5 mL 10%三氯乙酸，4000 r/min离心10 min。取上清液2 mL加入2.5 mL蒸馏水，再加入0.5 mL 0.1% FeCl 3 ，充分反应后于700 nm处测定其吸光度。吸光值越高，还原性越强 [ 13 ] 。 2.5.4 . 菌丝体多糖的抗脂质过氧化能力 吸取1:25卵黄悬液0.4 mL，加入菌丝体多糖溶液0.8 mL，再加入25 mmol/L的FeSO 4 溶液0.4 mL，用pH7.4的磷酸缓冲液补充至4 mL，37℃振荡15 min。加入20%三氯乙酸1 mL，3500 r/min离心10 min。取上清液4 mL，加入质量分数为0.8%的TBA(硫代巴比妥酸)溶液1 mL，封口，沸水浴15 min，532 nm处下测定吸光度值。不加样品管的吸光度为A 0 ，样品加磷酸缓冲液的吸光度为A 1 。待测液对卵黄脂质过氧化的抑制率(%) = [(A 0 − A 1 )/A 0 ] × 100% [ 14 ] 。 1:25卵黄悬液的配制：新鲜鸡蛋去鸡蛋清，蛋黄用pH 7.4的0.1 mol/L磷酸盐缓冲溶液配成1:1悬液，磁力搅拌10 min。取1 mL该混合液于25 mL容量瓶中，pH 7.4磷酸盐缓冲液定容至25 mL，4℃冷藏备用。 2.5.5 . 菌丝体多糖的超氧阴离子自由基清除能力 采用邻苯三酚自氧化法测定菌丝体多糖对超氧阴离子自由基的清除能力 [ 15 ] 。取0.05 mol/L Tris-HCl缓冲液4.5 mL(pH 8.2)，25℃水浴中预热20 min后再加入l mL菌丝体多糖溶液和0.4 mL邻苯三酚溶液(25 mmol/L)，混合均匀后置于25℃水浴中反应4 min，加入1 mL 8 mol/L HCl终止反应，在299 nm处测定吸光度A 1 ，在空白对照组中以相同体积的蒸馏水代替溶液样品，吸光度记为A 0 。超氧离子自由基清除率按照公式E/% = (A 0 − A 1 )/A 0 × 100%计算。 2.6. 数据处理 每一试验均重复三次，所示数据为三次结果的平均值。
system: [SYSTEM]

user: 4. 结论 本实验应用响应曲面法对撕裂蜡孔菌菌丝体多糖提取的工艺条件进行优化，对真菌资源的综合利用具有较好的现实意义。优化后最佳水提取多糖的工艺条件为：微波时间32.27 s、料液比1:24.17、浸提时间36.56 min、微波功率480.49 W、浸提温度37.89℃。根据实际操作情况，修正最佳提取工艺条件为：微波时间30 s、微波功率480 W、液料比1:20、浸提时间30 min、浸提温度35℃。在此修正条件下菌丝体多糖实际得率为6.437%，比理论预测值6.523%低0.086%，表明所得模型可以较好的预测撕裂蜡孔菌多糖的提取情况。 撕裂蜡孔菌菌丝体多糖表现出一定的抗氧化能力，且抗氧化效果与浓度表现出良好的相关性。多糖对羟自由基和过氧化氢的清除作用较好，抗脂质过氧化能力和清除超氧阴离子自由基能力最强，但还原能力比较弱。 综上所述，撕裂蜡孔菌多糖具有较强的抗氧化能力，在开发具有清除自由基、抗脂质过氧化、延缓机体衰老的天然抗氧化剂用品方面具有较好的前景。
system: [SYSTEM]

user: 摘 要 通过对脆性煤层工作面煤壁片帮的机理和影响因素分析，并运用PFC软件对回采过程进行离散元数值模拟分析，对超大采高回采过程中可能形成的片帮情况以及片帮规模有了定量的预测。明确了大海则煤矿大采高工作面片帮属于典型的拉张型板裂破坏，与煤岩层的受力、结构的变形、煤岩物理力学性质有关。制定了采用提高围岩支护强度，提升煤体强度，优化液压支架结构设计，改进回采工艺以及保证煤层节理与煤壁倾向的综合治理技术方案，来减少煤壁片帮。 关键词 :煤矿，大采高工作面，PFC数值模拟，片帮治理 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 近年来，国内外对大采高综采工作面开采越来越重视，大采高综采技术也日臻完善，已成为厚煤层开采工艺的重要发展方向之一 [ 1 ] [ 2 ] [ 3 ] 。在该项技术方法中，采煤工作面过高，容易发生片帮。片帮煤过多、过大时会压死溜槽链子；可能砸坏电缆槽，损伤采煤机链带或电缆；还可能阻塞采煤机通过，卡在转载机进料口堵塞煤流，影响设备正常运转。煤壁片帮也给人员、设备带来一定的危害，易诱发冒顶事故，严重影响了大采高综采工作面的产能和效能，对其有效地控制成为大采高综采技术工作面现场管理的重点和难点。目前，对于大采高综采的研究国内外以前都仅仅局限在对支架支撑能力的研究和现场工作阻力与支架适应性的分析，而对于大采高条件下煤体的破坏范围及片帮机理却研究较少 [ 4 ] [ 5 ] 。
system: [SYSTEM]

user: 2. 工程概况 大海则煤矿主采2号煤层，煤层倾角0.6˚，埋深547.95~654.95 m，煤层厚度4.41~10.97 m，平均煤厚6.32 m。直接顶板一般以泥岩、粉砂质泥岩、粉砂岩为主，砂泥岩互层厚0.55~10.07 m，一般为1~3 m，强度较低。底板以泥质粉砂岩、粉砂质泥岩为主，粉砂岩主要分布北西部，泥岩主要分布在东南部，厚1~3 m左右，局部地段有0.1~0.4 m炭质泥岩，其下多为细粒砂岩。 首采的201及206盘区煤层赋存均稳定，其中201盘区2号煤层厚度4.82~6.52 m，平均厚度5.6 m；206盘区2号煤层厚度6.27~8.20 m，平均厚度6.9 m。2号煤层的天然容重平均1.29875 g/cm 3 ；天然状态抗压强度9.89~14.2 MPa，饱和状态抗压强度8.06~14.2 MPa。煤层的硬度天然状态下仅有1.0~1.4左右，属于脆性煤层。
system: [SYSTEM]

user: 3. 煤层工作面片帮原因分析 3.1. 脆性煤层工作面片帮机理 脆性煤层工作面煤壁的破坏属于典型的拉张型板裂破坏，拉张型板裂化煤壁失稳的形成一般经历三个阶段：首先工作面煤壁在切向应力作用下产生与工作面相平行的张性破裂面并将煤壁分离为互相平行的煤板；其次煤板在切向应力的进一步作用下发生压缩变形并积聚应变能；最后煤板发生屈曲失稳将积聚的应变能突然释放形成失稳。 脆性煤层工作面煤壁内张性破裂面的产生可从断裂力学角度予以解释。煤岩内部存在随机分布的微裂纹，煤岩在无围压承受一定大的轴向压力时，微裂纹发生摩擦滑动和自相似性扩张后发生弯折扩展，由于不存在围压，弯折裂纹沿最大主应力方向即轴向力方向继续发展，煤层受压发生劈裂破坏形成平行于轴向力的破裂面。 工作面每推进一段距离，新暴露的煤壁围岩应力状态由初次状态转变到二次状态，切向应力增大，径向应力降低为零。煤壁内部微裂纹发生摩擦滑动和自相似性扩张后发生弯折扩展，由于径向应力降低为零，弯折裂纹会沿切向应力方向继续发展，微裂纹贯通后形成平行于开挖面的破裂面，即张性破裂面。 拉张型破裂面产生后，形成平行于工作面的煤板，形成的煤板若没有遭到进一步扰动破坏，随着工作面围岩变形增加，煤板积聚能量，当储存能量达到一定程度，煤板遭到扰动后就会失稳发生断裂，产生的煤块被抛出或大面积垮落，造成工作面煤壁灾变和失稳。 3.2. 脆性煤层工作面片帮影响因素 煤壁大规模片帮失稳是一个较为复杂的力学过程，它除了与地质构造有关外，还与煤岩层的受力、结构的变形、煤岩物理力学性质有关。对于浅埋煤层，由于风化、裂隙发育、完整性差，煤岩体内能量耗散殆尽，所以工作面开采时往往出现零星的片帮而几乎没有大面积煤壁失稳发生。 煤壁内层对临空面煤壁层的作用可视为弹性约束，受纵向力如图1所示，在重力场 σ y 和构造应力场 图1. 煤壁薄板受力模型。(a) 工作面割煤高度；(b) 工作面割煤宽度；t．工作面板裂厚度 σ x 作用下，工作面煤壁产生向临空面的变形，同时受到内煤层的制约。层间受拉是脆性，当临空面煤壁与内层之间的拉力或相对位移达到一定值时出现离层。但是由于应力的重新分布和内岩层约束扼制了离层的迅速扩散。煤板上的无弹性支承段x随变形增长，裂隙扩大，剩余弹性支承迅速断裂，造成煤壁大面积片帮失稳，如图2所示。从断裂力学角度看，煤岩体内存在着微裂纹，当受纵向力作用时，裂纹将会平行或偏向最大主应力的方向扩张，形成劈裂破坏。所以对于脆性煤层，工作面板裂片帮并非完全呈薄板状，有可能呈现“洋葱片”状、中间厚两边薄的透镜状或楔形状形状。
system: [SYSTEM]

user: 4. 脆性煤层工作面片帮数值模拟分析 为了进一步分析超大采高回采过程中可能形成的片帮情况以及片帮规模，采用PFC软件对回采过程进行离散元数值模拟分析。 基础模型为宽度80 m，各模拟煤岩层的厚度分别为：煤层7 m、6 m、5 m、4 m，粉砂岩顶板8 m，中粒砂岩顶板8 m，泥质粉砂岩顶板9 m。其上载荷全部施加到上方载荷层。开挖过程考虑到边界效应的影响，左侧预留20 m煤柱，右侧预留20 m煤柱。 4.1. 煤壁片帮的规模分析 煤壁片帮的规模如图3所示。7 m割煤高度时煤壁相对稳定高度在4.5 m左右，该范围内发生轻微片帮，片帮深度在0.5 m以内；上方2.5 m范围内发生较严重的片帮，片帮深度达1 m，呈现三角形区域整体滑落状态。6 m割煤高度情况下片帮情况与7 m类似，中上部约2 m范围内发生深度达1 m的片帮；5 m、4 m割煤高度下仅在煤壁发生轻微片帮，片帮规模较小。 4.2. 煤壁片帮前方支撑压力分析 工作面前方的支承压力分布情况如图4示。在煤壁前方均形成明显的应力集中现象，煤壁处的煤体应力发生明显降低。7 m割煤高度时超前支承压力的峰值区域在工作面煤壁前方5 m处，割煤高度6 m、5 m、4 m时，支承压力峰值区域在工作面煤壁前方2 m处。
system: [SYSTEM]

user: 5. 预防片帮治理技术 根据脆性煤层工作面片帮发生原因及PFC数值模拟分析，决定从提高围岩支护强度，提升煤体强度，优化液压支架结构设计，改进回采工艺以及保证煤层节理与煤壁倾向方面采取技术措施，来降低煤壁损伤，减少煤壁片帮，控制片帮事故。 图2. 煤壁板裂化片帮模型。x．煤壁弹无性支承段；q 1 ．煤壁弹性支承段；l．割煤高度 图3. 煤壁片帮规模分析图。(a) 割煤高度7 m；(b) 割煤高度6 m；(c) 割煤高度5 m；(d) 割煤高度4 m 图4. 煤壁片帮前方支承压力分布情况图。(a) 割煤高度7 m；(b) 割煤高度6 m；(c) 割煤高度5 m；(d) 割煤高度4 m 5.1. 治理方案 5.1.1. 提高支架初撑力和前端支顶力 采高增加后工作面片帮加剧的一个重要原因是前方支承压力增加，范围增大，导致工作面前方煤体，尤其是煤壁附近煤体的承载能力严重降低。提高液压支架强度和提高顶梁前端支顶力使液压支架更好地适应大采高开采的围岩关系，是控制工作面片帮冒顶的最有效措施之一。 5.1.2. 超前加固，提升煤体强度 对煤壁节理裂隙发育区域，尽可能增加护帮长度可以有效改善前方煤体内微裂隙的发育，减小能够发展成为裂隙的弱面数量，从而控制工作面片帮。 对煤体松散、破碎区域，特别是煤壁片帮特别严重区域，采用加强支护或煤壁固化的方法提前加固，如打设锚杆或注射玛丽散、博特威等固化剂。 5.1.3. 优化液压支架结构设计 进行液压支架结构优化设计，提高液压支架强度，使液压支架更好地适应大采高开采的围岩关系。采高7 m的液压支架，顶梁设计为整体顶梁加伸缩梁结构，伸缩梁前端铰接护帮结构为三级机构，一级和三级护帮采用小四连杆结构，可翻转180˚，护帮收回时，二级护帮与一级护板保持水平，三级护帮折叠收回。 5.1.4. 改进回采工艺，降低煤壁损伤 加快工作面推进速度，减少煤壁暴露时间，可以减少支承压力对煤体的作用时间，降低煤壁损伤程度，从而减少煤壁片帮。当采煤机割煤通过后，采用及时移架的方式，及时地把伸缩梁伸出，对顶板进行支护。若工作面出现局部片帮和大采高支架歪倒、陷底、挤架等现象，应及时调整。 5.1.5. 保证煤层节理倾向煤壁 设计工作面时避开上覆煤层遗留煤柱对工作面造成集中压力，工作面推进方向应保证煤层节理倾向煤壁，工作面应俯斜开采，避免仰斜开采。 5.2. 应用效果分析 通过在首采的20101工作面和20601工作面采用上述片帮综合治理措施情况来看，相比相似煤层大采高工作面，工作面总体片帮现象比较轻微(补连塔22307工作面，采高6.8 m，煤壁片帮深度2.6 m；大柳塔52302工作面，采高6.6 m，煤壁片帮深度1.8 m)。周期来压期间，工作面煤壁片帮最大长度3.7 m，最大深度0.3 m，平均片帮率8%。总之，大海则煤矿首采的20101工作面和20601工作面采取一系列片面帮措施后，煤壁片帮程度被控制在一定范围内，保证了工作面的正常生产。
system: [SYSTEM]

user: 6. 结论 1) 通过对脆性煤层工作面煤壁片帮的机理和影响因素分析，认为片帮属于典型的拉张型板裂破坏。煤壁大规模片帮失稳是一个较为复杂的力学过程，它除了与地质构造有关外，还与煤岩层的受力、结构的变形、煤岩物理力学性质有关。 2) 运用PFC软件对回采过程进行离散元数值模拟分析，对超大采高回采过程中可能形成的片帮情况以及片帮规模有了定量的预测。 3) 确定了采用提高围岩支护强度，提升煤体强度，优化液压支架结构设计，改进回采工艺以及保证煤层节理与煤壁倾向的综合治理技术方案，来减少煤壁片帮。
system: [SYSTEM]

user: 摘 要 突泉县位于中国十四个集中连片特困地区——大兴安岭南麓山区，鉴于突泉县特殊的地域类型，分为北部山区，中部丘陵地区，以及南部平原地区，因此突泉县贫困问题具有特殊性。以突泉县三个地区根据实地调研抽选三个行政村为主要研究对象，基于人地关系地域系统理论，结合空间自相关方法，分析贫困村分布规律和贫困地域特征，探究突泉县农村贫困的主导因素，旨在为突泉县精准扶贫工作提供依据。研究结果表明：突泉县农村贫困具有高度空间聚类的特点，在村域空间单元水平上，Moran’s I指数的值为0.38 > 0，说明贫困村空间分布呈现出聚类分布的特征，即贫困程度相似的村域总体上倾向于聚集分布。其中Z得分为12.54，则说明随机产生此聚类模式的可能性小于1%。P检验也为0，说明通过检验。且分为高值聚类和低值聚类。 关键词 :贫困，地域差异，比较研究，精准扶贫 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着人类社会科技经济的快速发展，不断增加的人口与自然环境的矛盾日益凸显，由此带来的贫困一直是困扰发展中国家的问题之一。中国作为人口大国，长期以来一直在尝试用各种方法减少贫困，提高人民生活水平。改革开放以来，我国政府主导下实施大规模扶贫开发，特别是《国家八七扶贫攻坚计划(1994~2000年)》《中国农村扶贫开发纲要(2001~2010年)》及新《纲要》的相继实施，有力地促进了贫困地区县域经济和农村转型发展，扶贫开发的成效举世瞩目 [ 1 ] [ 2 ] [ 3 ] 。1978~2014年，全国贫困人口由2.5亿降至7017万人，贫困发生率也从30.7%降至7.2%。目前中国的农村贫困县是每人年收入2300元(按2010年不变价)，即每人每天约6.3元。据国务院扶贫办统计，全国有14个集中连片特困地区，片区县及国家扶贫开发重点县共832个，贫困村12.8万个，贫困户29，485万个，贫困人口7017万人。贫困人口不仅规模大，分布广，而且致贫原因复杂，脱贫难度大，贫困地区的转型发展、贫困人口的脱贫事关2020年全面建成小康社会目标的成效。贫困问题一直是近些年来中国进步发展所面临的问题。而内蒙古作为北方农牧交错和生态环境脆弱地带，人地关系整体上趋于紧张。海山指出内蒙古农牧交错带是我国自然生态环境破坏最严重、经济最贫困、文化科技教育最落后的地区之一，也是我国人地矛盾十分尖锐的地区之一。农牧交错带是自治区和全国经济最贫困地带之一 [ 4 ] 。修长柏提出交错带的贫困问题是交错带区别于其他区域的一个显著特征，有必要对这个问题进行深入研究 [ 5 ] 。突泉县作为国家级重点贫困县，根据扶贫办2017年最新统计数据显示，突泉县仍有20，312人的贫困人口，占兴安盟总贫困人口的25%，在2020年前，脱贫攻坚任务仍然十分艰巨。造成贫困的原因复杂多样，学者进行了大量相关研究。但是国内农村贫困的相关研究仍以定性研究居多，比如农村贫困特征，致贫原因，脱贫政策以及制度研究等 [ 6 ] [ 7 ] [ 8 ] ，而对农村贫困化进行定量研究的较少，现有文献主要包括，如李双成运用BP神经网络图研究区域贫困 [ 9 ] ，杨国涛基于GQLorenz曲线的参数方法估计贫困现状 [ 10 ] ，韩林芝、邓强等运用灰色关联分析致贫原因 [ 11 ] ，王小林、Alkire利用中国健康与营养调查数据对中国城市和农村家庭多维贫困进行了测量 [ 12 ] ，王艳慧构建基于“双临界值”的“维度加总分解”算法进行了“县级–村级”的贫困人口多维贫困量算和分析；借助Kriging法对村级多维贫困测算结果进行空间插值处理，系统分析研究区多维贫困状况空间分布格局 [ 13 ] 。 还有一些文献是基于地域差异的角度研究贫困，分析了中国贫困地域类型的划分 [ 14 ] [ 15 ] ，还有是进行实证研究 [ 16 ] [ 17 ] ，也有运用地理探测器多元线性等模型方法，诊断出县域农村贫困化分异的主导因素，揭示了农村贫困化分异的动力机制，提出了不同贫困化地域类型的扶贫政策与模式 [ 18 ] ，但是对农村贫困村域地域差异的定量研究较少。 根据前期调研测算突泉县贫困率发现，基于对突泉县实地调研的了解，因此本研究基于人地关系地域系统理论，选取突泉县的三个行政村研究其贫困成因，提出扶贫策略，旨在为突泉县精准扶贫工作提供思路与参考。
system: [SYSTEM]

user: 2. 研究区概况 2.1. 地理位置 突泉县是位于内蒙古自治区东北部、兴安盟中南部，北纬45˚11′25″~46˚5′12″，东经120˚48′45″~ 122˚10′20″。总面积4889.50平方公里，总耕地面积为1776.5676平方公里。A村位于西北部山区，B村位于中部浅丘陵区，C村位于南部平原。见图1。 2.2. 地貌特征 根据三个贫困村的高程，分层设色，可以直观的看出三个村分别位于山地，丘陵，平原。其中A村 图1. 研究区区位图 平均高程约560米；B村平均高程约400米；C村平均高程约240米左右。研究区海拔自西北向东南呈明显的递减趋势。见图2。 2.3. 自然环境 突泉县处于农牧交错带范围。地处北温带，西依大兴安岭，东连松嫩平原。县境在温带季风气候区内，由于远离海洋，有明显的大陆性特点。由南向北，降水逐渐增多：由北向南，气温明显上升。县境地势西北高，向东南逐渐降低。土壤类型主要为栗钙土和暗棕壤。 2.4. 社会经济 据统计，2015年全县地区生产总值完成734，116万元，其中第一产业生产总值占31.85%，第二产业占44.05%，第三产业占24.10%；2015年总人口有30.38万人，其中蒙古族占20.54%，其他少数民族仅占总人口数的9.91%。
system: [SYSTEM]

user: 3. 空间自相关分析与样本村选择依据 3.1. 全局空间自相关 通常情况下，地理现象会显示出聚类分布的特征。事实上，Moran’s I的值的大小仅仅是一种机会，这些值需要和其期望值相比较，并用其标准化的Z值来解释 [ 19 ] [ 20 ] 。 见表1结果显示，在村域空间单元水平上，Moran’s I指数的值为0.38 > 0，说明贫困村空间分布呈现出聚类分布的特征，即贫困程度相似的村域总体上倾向于聚集分布。其中Z得分为12.54，则说明随机产生此聚类模式的可能性小于1%。P检验也为0，说明通过检验，突泉县贫困村分布空间上呈现聚类特征。 3.2. 局部空间自相关 运用General G探测局部聚类的情况，见表2结果显示：根据Z得分为3.09，则说明随机产生此高聚类模式的可能性小于1%。说明突泉县贫困村为高聚类模式。 3.3. 样本村选择依据 选村依据：根据测算贫困发生率，运用空间自相关分析法得出突泉县的贫困村在空间分布上具有空间集聚的特征，且均为正相关。鉴于突泉县的特殊地形，分为北部山区，中部丘陵地区，以及南部平原地区 [ 21 ] ，所以在这三个不同的地形区抽选三个有代表性的行政村进行比较研究。预探究贫困的地域差异。