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

user: 6. 结论 1) 英西地区E 3 2 储层为咸化湖盆背景下的深湖–半深湖混积岩沉积，储层矿物有方解石、(铁)白云石、石英、长石、(硬)石膏、黄铁矿等。成岩作用类型有：压实作用、胶结作用、破裂作用和溶蚀作用，其中储集空间类型以晶间孔隙为主，伴有粒间孔隙、溶蚀孔隙和微裂缝的发育。孔隙度主要介于4%~12%之间，平均值为7.85%，渗透率普遍小于1 mD，平均值为0.469 mD。孔渗关系复杂，储层具备低孔、低渗的特点。 2) 储层白云石含量增加有利于孔隙度的增大，含砂或砂质灰云岩有利于微裂缝的发育，溶蚀作用多发育在碳酸盐矿物中，且多与溶蚀缝伴生，广泛发育的泥粉晶级别晶粒是导致储层致密的根本原因。根据矿物含量及成岩作用类型，将成岩相划分为4种类型：成岩微裂缝相、溶蚀相、晶间孔相和压实致密相。四种类型成岩相中压实致密相和晶间孔相广泛发育，成岩微裂缝相和溶蚀相为有利成岩相。 图16. 英西地区狮41-2井下干柴沟组上段成岩相测井连续识别 3) 基于测井对储层原始沉积矿物组分及成岩作用类型和强度的响应，通过挑选成岩相响应敏感的测井参数，运用主成分分析法对储层成岩相进行了定量表征。其中，成岩微裂缝相具有低密度值、中自然伽马值、高碳酸盐矿物含量的特征，溶蚀相具有中密度值、低自然伽马值、高碳酸盐矿物含量的特征，晶间孔相具有中高密度值、中自然伽马值、中碎屑矿物含量的特征，压实致密相具有高密度值、高自然伽马值、高黏土矿物含量的特征。此外，英西地区复杂的成岩过程导致了单井纵向成岩相变化频率快，成岩微裂缝相常与构造裂缝伴生，因此是主要的裂缝型“甜点”发育带，而在无构造裂缝的叠加作用时，溶蚀相物性最好，压实致密成岩相物性最差。
system: [SYSTEM]

user: 摘 要 为探究庐山森林土壤中N、C、S元素迁移转化的规律，本文在野外观察的基础上，选取了6个受人类活动干扰较小的采样点进行样品采集，并对其元素含量进行了测定与分析，结果表明：1) 庐山森林土壤中TN、TC、TS含量分别为0.03%~0.48%、0.18%~7.00%、0.02%~0.13%，TN、TC值变化幅度大，TS值变化幅度小；2) 除RZF采样点外，土壤剖面不同发生层TN、TC值随着深度的增加有显著的降低；RZF点分层不明显，故各层TN、TC值差异不显著；3) LL采样点因地形坡度大，土壤元素易被淋失，因此元素含量没有表现出与其他5处采样点相同的变化特征，即TN、TC值随着海拔上升而增加；4) 庐山森林土壤C/N值远小于适宜值25，也低于全国生态系统平均值，变化趋势与TN、TC一致；5) 土壤C/S值均小于200，有机硫矿化作用强。进一步讨论发现，庐山森林土壤生态系统中N、C、S元素的迁移转化主要受温度、地形、地表枯枝落叶量及植被类型影响。 关键词 :庐山，森林土壤，碳氮硫，有机质，环境因子 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. 引言 作为地球上最大的陆地碳库，森林生态系统碳库的碳储量约占陆地总碳库的46% [ 1 ] ，而森林土壤碳库又是森林生态系统碳库的主体，因此其对全球碳循环和碳平衡的重要性不言而喻，近年来对森林土壤碳含量研究也已经成为气候变化研究的主要关注点之一 [ 2 ] [ 3 ] [ 4 ] [ 5 ]。研究显示，土壤碳储量约为植被碳储量的2.5~3倍，大气碳储量的2~3倍。森林土壤碳库由有机碳库与无机碳库两大部分组成，有机碳库是指存在于土壤中的各种含碳有机物质，其在好氧条件下可以通过微生物矿化分解，最终以CO 2 形式进入大气；无机碳库即全球干旱半干旱地区土壤中的碳酸盐类，包括原生性碳酸盐和次生性碳酸盐 [ 6 ]。森林土壤无机碳储存量相对较少，并且变化幅度不大，有机碳储量则十分庞大，约占全球土壤有机碳储量的39%，森林生态系统有机碳库的2/3 [ 7 ]。同样是以CO 2 形式向大气中排放碳，但因其巨大的碳储量，土壤每年排放的碳通量约为化石燃烧排放量的12~16倍 [ 8 ]。除此之外，氮和硫元素也是土壤元素的主要组成。氮元素是土壤有机质的重要组成部分，因此土壤氮含量高低是评价土壤质量的基本标准之一；同时土壤氮含量变化还会改变微生物呼吸作用，进而影响土壤CO 2 排放 [ 9 ]。硫的营养功能和循环机制与氮相似，也是植物生长发育所必需的矿质营养元素，植物所需的硫主要来源于根系对土壤中无机硫的吸收 [ 10 ] [ 11 ] ；硫化合物以多样的形态在陆地生态系统中发挥着重要的生态功能，是参与生物地球化学循环的主要元素 [ 12 ]。 土壤圈是地球表层与大气圈、生物圈、水圈、岩石圈相交的界面，并进行着物质循环和能量转换，常因成土环境的变化导致土壤性质的改变，进而影响土体内部物质的淋溶与淀积，导致土壤在发生发展的过程中产生质和量的变化 [ 13 ]。土壤圈的元素转换会影响大气中温室气体成分的动态平衡，进而导致全球气候变化；不同的土壤类型，其元素迁移转化过程也不同。近些年来，通过对美国北部大平原、我国松嫩平原等地土壤中碳氮硫元素迁移转化过程的研究，学者们发现土壤中碳氮硫元素变化主要受到降水量、年均温、有机质含量、植被类型及植物根系发展程度等因素影响 [ 14 ] [ 15 ] [ 16 ] [ 17 ] [ 18 ]。庐山因其特殊的地理位置和条件，形成了一套典型的山地森林土壤垂直带谱；不同森林土壤类型以及周围环境因子的变化，会直接或间接影响土壤中碳氮硫元素的累积与消耗，使得对土壤元素组成特征的研究变得困难。 本文依据庐山山地森林土壤垂直带谱的演变规律，选择具有代表性的采样点，试图找出森林土壤元素含量变化与森林土壤类型及其周围环境因子之间的相关性，以期探讨出庐山山地森林土壤生态系统演变对气候变化的影响。
system: [SYSTEM]

user: 2. 研究区概况 庐山位于江西省九江市庐山市境内，西北濒临长江，东南濒临鄱阳湖，范围介于115˚52'E~116˚8'E，29˚26'N~29˚41'N之间。山体南北长约29 km，东西宽约15 km，海拔25~1474 m，主峰汉阳峰为最高峰。气候属典型的亚热带季风气候，全年温和湿润，年均气温为11.4℃，年均降水量为1917 mm。另外，作为典型的地垒式断块山体，山地小气候特征显著，年均相对湿度78%，年均雾日191天。 此外，在不同的海拔高度下，庐山山体表面水热组合状况略有不同，进而影响气候、生物、土壤等环境条件，导致自然带呈现出垂直方向上的明显变化。庐山土壤类型的分布状况为：400 m以下山麓及山麓以外的丘陵和沉积阶地为红壤和黄壤分布区域；400 m~1200 m山坡地带，为黄壤和棕壤的分布区域；1000 m以上的山地，为山地棕壤和亚高山草甸土分布区域。植被由山麓到山顶垂直分布主要呈现出常绿阔叶林 → 常绿、落叶阔叶混交林 → 落叶阔叶林的变化规律 [ 19 ]。
system: [SYSTEM]

user: 3. 材料与方法 为详细了解庐山不同海拔高度土壤碳氮硫元素含量的变化特征，通过野外考察，根据“典型性、代表性、一致性”的采样原则，选取了6个受人类活动影响小，且垂直分层明显的采样点。采样点地理位置与详细信息如图1、表1所示。采样时将地面的落叶与杂物清理干净，用土铲挖出垂直的土壤剖面， 图1. 庐山各采样点地理位置示意图 根据土壤的颜色、干湿度、松紧度、新生体等指标将土层分为腐殖质层、淋溶层和淀积层三个发生层。在每一土层取土样约30 g装入样品袋中，写上编号、采集地点、采集时间、地表植被类型、土壤名称、采集深度、土壤性质的描述等信息。样品采集完成后将土坑填平。 Table 1 The detailed information of geographic location and the soil profile of the sampling sites in Mt. Lusha 采样点 经纬度 海拔 /m 剖面 深度 /cm 采样点代号 采样点 深度 /cm 分层 界限 /cm 地表植被类型 分层概况 白鹿洞书院 116˚22'11''E，29˚31'9''N 190 160 BLD-6 6 1~9 常绿 阔叶林 腐殖质层A：浅黄色、紧实、较干 BLD-23 23 9~35 淋溶层E：黄色、紧实、干 BLD-160 160 35~160 淀积层B：红棕色、紧实、干 黄龙寺 115˚57'47''E，29˚33'11''N 900 55 HLS-10 10 2~16 常绿落叶阔叶 混交林 腐殖质层A：黑褐色、疏松、稍润 HLS-25 25 16~29 淋溶层E：褐色、疏松、稍润 HLS-40 40 29~55 淀积层B：黄色、紧实、润 大天池 115˚57'4''E，29˚33'37''N 950 48 DTC-6 6 2~11 黄山松 针叶林 腐殖质层A：暗棕色、疏松、稍润 DTC-17 17 11~22 淋溶层E：棕色、疏松、稍润 DTC-35 35 22~48 淀积层B：黄棕色、紧实、润 芦林 大坝 115˚57'19''E，29˚33'16''N 980 66 LL-4 4 1~6 柳杉 针叶林 腐殖质层A：黑灰色、松散、稍润 LL-20 20 6~33 淋溶层E：灰黄色、疏松、润 LL-47 47 33~66 淀积层B：棕黄色、紧实、润 日照峰 115˚59'15''E，29˚34'37''N 1240 90 RZF-6 6 1~9 芒草 腐殖质层A：深褐色、松散、稍润 RZF-30 30 9~45 淋溶层E：棕色、疏松、稍润 RZF-60 60 45~90 淀积层B：棕黄色、紧实、润 大月山 115˚59'56''E，29˚34'5''N 1300 113 DYS-30 30 2~40 落叶阔叶林，黄山松针叶林 腐殖质层A：暗棕色、疏松、稍润 DYS-45 45 40~55 淋溶层E：灰黄棕色、疏松、润 DYS-80 80 55~113 淀积层B：黄棕色、紧实、润 表1. 庐山各采样点地理位置及土壤剖面信息 样品带回实验室后，剔除土壤中的杂质(石子、枯枝落叶等)，1) 将样品置于室内自然晾干，研磨过10目筛；2) 取一定量的样品，加入5%的稀盐酸多次搅拌，不断加入稀盐酸至反应完全，浸泡二十四小时；3) 用中性去离子水水洗至中性(pH = 7)，烘干后将样品用玛瑙研钵研磨至粉末状，过200目筛；4) 根据样品中总氮、总碳和总磷的含量(预先根据不同发生层土壤样品性状测定了若干不同的控制样)称取一定量(约200 mg)的被测样品在锡纸紧密包裹下送入德国Elementar公司生产的Vario EL cube元素分析仪直接测定出TN、TC、TS和C/N、C/S值。使用Microsoft Office Excel 2007软件对数据进行处理分析，OriginPro 8.0和ArcGIS 10.1绘制图表(表2)。 Table 2 The contents of N， C， S and C/N， C/S values of the sampling sites’ soils in Mt. Lusha 分层 采样点 元素含量 N (%) C (%) S (%) C/N C/S 腐殖质层 BLD-6 0.13 1.11 0.02 8.33 55.50 HLS-10 0.23 3.12 0.04 13.55 78.00 DTC-6 0.24 3.49 0.04 14.56 87.25 LL-4 0.17 1.65 0.04 9.92 41.25 RZF-6 0.26 2.92 0.06 11.38 48.67 DYS-30 0.48 7.00 0.07 14.68 100.00 淋溶层 BLD-23 0.08 0.54 0.04 6.79 13.50 HLS-25 0.15 1.46 0.02 9.76 73.00 DTC-17 0.15 1.71 0.04 11.68 42.75 LL-20 0.08 0.74 0.03 8.84 21.76 RZF-30 0.21 2.10 0.13 10.18 16.15 DYS-45 0.32 4.33 0.07 13.38 61.86 淀积层 BLD-150 0.06 0.18 0.02 3.12 9.00 HLS-40 0.05 0.39 0.12 7.80 3.25 DTC-35 0.05 0.37 0.03 7.33 12.33 LL-47 0.03 0.13 0.02 4.88 6.50 RZF-60 0.16 1.67 0.06 10.68 27.83 DYS-80 0.20 1.98 0.06 10.05 33.00 表2. 庐山各采样点土壤N、C、S含量以及C/N、C/S值
system: [SYSTEM]

user: 5. 结论 1) 随着剖面深度的增加，庐山森林土壤中TN、TC含量不断下降；采样点海拔越高，土壤中TN、TC含量则越高。 2) 庐山森林土壤中TN、TC含量受外界环境因子的影响显著，影响因素包括土壤温度、地表枯枝落叶量、地形以及地表植被类型等。 3) 庐山森林土壤C/N值普遍低于全国生态系统的平均值，土壤氮的矿化率高；C/N值的变化趋势以及导致其变化的原因与TN、TC的一致；C/S值均低于200，土壤有机硫矿化作用较强，且底层土壤与高海拔地区土壤中有机硫的矿化度分别高于表层土壤与低海拔地区土壤；土壤C/N、C/S值的变化主要取决于土壤中有机碳含量的变化。
system: [SYSTEM]

user: 摘 要 通过分析基于模板匹配的含线特征混合像元亚像元定位算法存在的问题，提出改进的模板匹配新方法。在相关系数选择模板的基础上提出两种新的模板选择方法：1) 对含线特征混合像元八邻域拟合直线进一步确定模板，在一定程度上缩小了模板选择的范围；2) 为避免模板选择的不确定性，对含线特征混合像元八邻域及模板分别进行直线拟合，计算直线相关性，得出最佳模板。进一步地，结合像元引力解决包含三种及以上地物混合的含线特征像元的亚像元定位问题。实验结果表明所提出的方法能提高亚像元定位精度。 关键词 :混合像元，亚像元定位，模板匹配，相关系数，直线拟合 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 ] 。软分类技术，能够获取像元中每一个地物类别对应的丰度，获得与类别个数相等的丰度影像，有效的解决混合像元问题，但无法确定各种地物在像元空间中的具体位置 [ 3 ] [ 4 ] 。Atkinson [ 5 ] 首次提出了亚像元定位(subpixel mapping， SPM)的概念。 亚像元定位主要思想是将混合像元按需求放大比例进行分割，确定每个亚像元的所属类别即确定混合像元中不同地物类型具体空间分布。目前主要的亚像元定位算法主要有像元交换 [ 6 ]，亚像元/像元空间引力 [ 7 ]，神经网络 [ 8 ] [ 9 ]，软硬属性转换 [ 10 ] [ 11 ]，线性优化 [ 12 ] [ 13 ] [ 14 ] 等模型。上述方法都忽略了混合像元内部不同的结构特征。针对该问题，国内学者许雄提出了顾忌地物空间分布的亚像元定位方法 [ 15 ]，将混合像元分为边界混合像元和线特征混合像元，利用特征探测、定位和决策等多个智能体分别实现不同特征的获取、不同特征的定位和最后的优化处理；葛勇等通过形状密度指数将地物目标分为点、线、面三类，针对不同的目标使用相应的定位方法 [ 16 ] 。对含线特征的混合像元内的亚像元定位用了模板匹配方法，但是该方法无法解决多个模板中如何智能选取最优模板问题，为此本文提出两种模板选择方法，实验证明所提出的模板选择方案能有效解决模板选择时多个模板最优问题。通过分析基于模板匹配的亚像元定位算法存在的问题，提出了基于改进模板匹配的线特征亚像元定位算法。
system: [SYSTEM]

user: 2. 含线特征混合像元亚像元定位的模板匹配法 2.1. 含线特征混合像元 由于传感器空间分辨率的限制和地物的复杂多样性，遥感图像中混合像元普遍存在。Fisher等 [ 17 ] 将遥感图像中混合像元划分为四种类型：亚像元目标、边界混合像元、致密型混合像元和线性目标，如图1所示，其中每个小方格表示1个像元，不同的颜色表示混合像元内不同的地物及其不同类型的空间分布。如图1(d)所示的线性目标对应的混合像元称为含线特征混合像元，该类混合像元一般为某类线性地物，如道路、河流等，由于其宽度和长度分别小于和大于影像空间分辨率，从而造成其贯穿多个混合像元造成的。 每种地物在一个混合像元内占的面积百比分称为该地物在该混合像元内的丰度。 图1. 混和像元的四种不同类型。(a) 亚像元目标；(b) 致密型混合像元；(c) 边界混合像元；(d) 线性目标 2.2. 亚像元定位的模板匹配法 基于模板匹配的线状地物亚像元定位方法基本原理包括以下3个步骤 [ 16 ] ： 1) 对每个混合像元，确定像元尺度下最佳匹配的线状模板。 图2. 3 × 3模板 以图2所示的3 × 3的线状模板为例，每个模板可用二进制矩阵表示为： T k = [ t ( − 1 ， − 1 ) t ( − 1 ， 0 ) t ( − 1 ， 1 ) t ( 0 ， − 1 ) t ( 0 ， 0 ) t ( 0 ， 1 ) t ( 1 ， − 1 ) t ( 1 ， 0 ) t ( 1 ， 1 ) ] (1) 其中，矩阵各元素 t ( m ， n ) ∈ { 0 ， 1 } ， m ， n = − 1 ， 0 ， 1 ， t ( m ， n ) = 1 表示对应图1中的灰色小方格， t ( m ， n ) = 0 对应图1中的白色方格。同理，可以构建重建尺度为 的二进制模板矩阵。 根据各模板和以混合像元P ij 为中心的 s × s 邻域各像元中线状分布地物的丰度之间的匹配程度确定P ij 的最佳匹配模板。匹配程度的度量方法有多种 [ 17 ] 。以标准相关匹配(Standard correlation， SC)方法为例，对混合像元中呈线性分布的第c类地物进行亚像元定位。设 w = ( s − 1 ) / 2 ，计算图2中各模板和以混合像元P ij 为中心的 s × s 邻域各像元第c类地物丰度之间的相关系数 r i ， j ， k = ∑ m = − w w ∑ n = − w w   T k ( m ， n ) ⋅ S c ( i + m ， j + n ) ∑ m = − w w ∑ n = − w w     S c ( i + m ， j + n ) 2 ⋅ ∑ m = − w w ∑ n = − w w   T k ( m ， n ) 2 (2) 其中， T k ( m ， n ) 表示模板矩阵中第k个模板 ( m ， n ) 位置的元素值， S c 表示图像中第c类地物的丰度值。 由式(1)可知，r i ，j，k 的取值范围在0~1之间，值越接近1模板的匹配效果越好，反之效果越差。对图2中的所有模板依次计算相关系数，取值最大的模板为最佳模板。 2) 利用最优模板 T k 对混合像元P ij 中第c类地物进行定位。 根据丰度求得第c类地物在混合像元P ij 内所占的亚像元个数 N c = S c ( i ， j ) ⋅ s 2 ，建立像元和亚像元的统一坐标系 [ 17 ]，在统一坐标系中，计算亚像元v pq 与最优模板 T k 中各元素之间的欧式距离，离模板距离最小的 N c 个亚像元类别为c。 3) 像元内其他亚像元类别确定 分别计算像元内其他地物的亚像元个数，随机确定像元内剩余亚像元的类别。 2.3. 存在的问题 上述模板匹配亚像元定位方法存在如下问题： 1) 直接根据相关系数r i ，j，k 进行模板选择时最大r i ，j，k 对应的模板可能不唯一，如何从多个匹配的模板中选出唯一的最优模板仍是个问题； 2) 当混合像元内包含三种及以上地物时，用随机方法确定像元内其他亚像元类别不够准确。
system: [SYSTEM]

user: 3. 两种最佳匹配模板确定的新方法 针对问题1)提出标准相关匹配结合像元直线拟合(Pixel Line Fitting after Standard Correlation， SC-PLF)的方法，标准相关匹配结合直线拟合相关性(Line Fitting Correlation after Standard Correlation， SC-LFC)的方法。 3.1. 基于拟合直线的最佳匹配模板确定 根据混合像元P ij 及其八邻域中地物 的丰度值进行二值化，即8邻域中如果含有第c类地物则将该混合像元标记为1，否则标记为0。图3(a)为第c类地物8邻域丰度示意图，图3(b)中0~1值表示二值化操作后的结果。求出二值化后八邻域中值为1的像元中心坐标，基于线性最小二乘法拟合直线，如图3(b)所示。 图3. 直线拟合示意图。(a) 混合像元第c类地物八邻域丰度图；(b) 二值化后直线拟合 以重建尺度 s = 4 为例说明像元中心坐标的表示。将混合像元P ij 及其8邻域放在同一坐标系并以左上角为坐标原点，如图3(b)所示，以左上角为统一坐标系原点，将每一个小方格长度记为1，因此在该同一坐标系中，左上角的亚像元 p 0 ， 0 坐标为(0.5， 0.5)，左上角像元 p 0 ， 0 坐标为(2， 2)，像元 P i ′ ， j ′ 的坐标为 ( s ⋅ ( i ′ + 0.5 ) ， s ⋅ ( j ′ + 0.5 ) ) ，亚像元 p a ， b 的坐标为 ( a + i ′ ⋅ s + 0.5 ， b + j ′ ⋅ s + 0.5 ) 。 图4. 4种直线拟合结果 对于SC确定的若干模板，依次判断模板矩阵中值为1的像元是否被拟合后的直线所经过，并统计其个数N T 。设拟合直线L a 的方程为 A x + B y + C = 0 ，模板矩阵中值为1的像元 P m n 中心坐标为 ( x m n ， y m n ) ，则像元P mn 到直线L a 的距离为 d ( P m n ， L a ) = | A x m n + B y m n + C A 2 + B 2 | (3) 如果 d ( P m n ， L a ) ≤ 2 / 2 ，则像元P mn 被直线L a 经过。 N T 值最大的模板为最优模板。图4给出了图2中部分模板的直线拟合结果，其中各模板中直线所经过的灰色区域个数N T 分别为1，3，2，1，图4(b)所示模板为最优模板。 3.2. 标准相关匹配结合直线拟合相关性选择最佳模板 SC-PLF方法虽然进一步缩小了模板选择的范围，但是对于模板唯一性仍无法保证。为此本节提出了标准相关匹配结合直线拟合相关性(Line Fitting Correlation after Standard Correlation， SC-LFC)的方法，在基于标准相关性的模板匹配选择方法选择出的若干模板基础上，分别对所选的模板和含线特征混合像元邻域进行直线拟合，并计算各模板和像元拟合直线间的相关性，可以保证模板选择的唯一性，算法主要思想是： 1) 对混合像元P i ，j 及其邻域像元根据地物 的丰度进行二值化，拟合直线 L a = k 1 ⋅ x + b 1 ； 2) 对经过SC算法选择出的各模板，根据 t ( m ， n ) = 1 的位置拟合直线 L i = k 2 ⋅ x + b 2 ， i = 1 ， 2 ， ⋯ ， k ； 3) 利用式(4)计算L a 和直线 L 1 ， L 2 ， ⋯ ， L k 的相关性： L C i = w 1 ( k 1 − k i ) + w 2 ( b 1 − b i ) ，     i = 1 ， 2 ， ⋯ ， k (4) 其中，w1和w2为权重，k1和 k i ， i = 1 ， 2 ， ⋯ ， k 分别为直线L a 和直线L i 的斜率，同理，b1和bi则为直线L a 和直线L i 的截距。实验证明 w 1 = 0.8 ， w 2 = 0.2 时，所选模板能够保证唯一且最优。 图5给出了部分模板的拟合直线结果，LC i 的值越大说明相关性越大。因此通过选择LC i 最大值对应的模板得出最优结果。 图5. 直线拟合结果
system: [SYSTEM]

user: 4. 改进的含线特征混合像元亚像元定位流程 针对传统模板匹配方法存在的第2个问题，提出以下改进策略：如果混合像元P ij 内仅包含两种地物，则P ij 内其他亚像元P ab 类别可以直接确定；如果混合像元P ij 内地物类别大于两类，则提出基于像元引力确定其类别。采用环绕关系确定亚像元P ab 的邻域像元集 [ 6 ] ： N ( P a b ) = { P i j | ( a s ≠ i ∨ b s ≠ j ) ∧ ( | a s − i | ≤ 1 ∧ | b s − j | ≤ 1 ) } (5) 第c类地物对亚像元P ab 的引力为： P a b ( c ) = avg { S c ( x ， y ) d ( P i j ， P a b ) | P i j ∈ N ( P a b ) } (6) d ( P a b ， P i j ) = [ a + 0.5 − s × ( i + 0.5 ) ] 2 + [ b + 0.5 − s × ( j + 0.5 ) ] 2 (7) 其中 max c P a b ( c ) 为亚像元P ab 的类别。 图6给出了改进的含线特征混合像元亚像元定位(Modified subpixel mapping of linear feature mixed pixels， MSPM L )流程图。 图6. 含线特征混合像元亚像元定位过程示意图
system: [SYSTEM]

user: 6. 结论 模板匹配和直线拟合是目前常用的线特征定位算法。由于模板匹配算法提高亚像元空间分辨率主要依赖于模板选择，仅采用标准相关系数(SC)对模板进行选择并不能保证所选模板最优，即存在最优模板存在多个的情况。本文针对基于模板选择部分提出两种模板选择方案，分别为SC-PLF和SC-LFC。其次，为进一步对混合像元内其余地物更好地定位，可采用模板匹配结合像元引力的方法，实验证明所提出的基于模板匹配的含线特征混合像元定位方法不仅能保证线特征的连通性而且使得混合像元内其余地物在定位时更加符合空间相关性。
system: [SYSTEM]

user: 摘 要 连州汛塘地区紧靠汛塘岩株，区内断裂构造发育。根据1:5万水系沉积物测量显示Sn、Cu、Bi、Ag、As、Au、Zn、Pb元素套合好，1:1万土壤测量显示锡铜等元素在浓集中心含量高，变异系数大。本文重点分析该区元素异常特征以进一步圈定找矿远景区。 关键词 :汛塘地区，锡铜多金属异常，找矿潜力 Copyright © 2019 by author 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/3以上，国内个旧、马关和广西大厂三个超大型矿床储量占了全国70%以上 [ 1 ] 。但经过多年开采，很多矿山进入“老年期”，资源接续成为重大问题，严重制约矿山可持续发展，因此，通过地质调查增加该矿种储量，具有重大的社会经济意义。 连州地区以往地球化学综合研究程度较低，2014年广东省地质调查院在此地开展了1:5万水系沉积物测量工作，并进一步开展了1:1万土壤测量等查证工作，取得丰硕成果 [ 2 ] [ 3 ] [ 4 ] 。汛塘位于连州大路边镇，根据前人地质资料，该地区已发现矿点汛塘锡铜锌矿，本文通过该区水系沉积物地球化学特征结合研究区土壤剖面测量工作发现汛塘锡铜多金属化探异常区有良好的找矿潜力，可为下一步地质找矿工作作铺垫。
system: [SYSTEM]

user: 2. 地质矿产特征 汛塘地区出露地层主要为石炭系地层，主要包括石磴子组(C 1 s)、测水组(C 1 c)、梓门桥组(C 1 z)、船山组(C 2 c)，偶见大埔组(C 2 dp)出露，总体岩性以灰岩及白云岩为主。南西部出露二叠系地层，为栖霞组(P 1-2 q)、茅口组(P 2 m)、童子岩组(P 2 t)，岩性为灰岩及砂岩。西部地层走向为北东向展布，东部走向以近东西向展布，总体上地层东老西新。 区内紧靠汛塘岩株，局部见侵入岩脉，燕山期岩浆侵入活动影响明显，与该区成矿作用关系密切。 汛塘异常区断裂发育，受北东向的汛塘断裂与山塘断裂、北西向的河佳汉断裂影响明显。北西向断裂老，北东向断裂新，总体上受北东向断裂影响范围更大。 本区见一处锡铜多金属矿点，汛塘村民采砂锡矿历史久远。见大理岩化、褐铁矿化、铜绿、铜蓝、硅化等现象。
system: [SYSTEM]

user: 3. 地球化学特征 该异常以Sn、Cu、Bi、Ag、As、Sb、Zn、Pb单元素异常范围圈定(图1)，其中异常最强为Sn、Bi，其次为Cu、Pb、Zn与Ag、As、Sb元素组，属中高温成矿元素为主的综合异常，近等轴状，长轴北东向、长5 km。异常较集中，呈面状分布，具南北向延伸形态，Sn、Bi、Cu等元素具有明显的三级浓度分带，且套合较好。异常地球化学参数表见表1。
system: [SYSTEM]

user: 4. 异常检查 在汛塘锡浓集中心，部署了1:1万土壤剖面测量，图2表明：1:1万土壤异常浓集中心较明显，大多浓集中心集中在北部，且锡铜等元素含量极高，为成矿元素最集中区域。其中，锡、锌异常呈北东向走向，未封闭，可进一步向北东向开展工作；钨铜铋异常存在北西向高值区，且未封闭，可在浓集中心及北部开展工作。 土壤测量结果(表2)表明：土壤中锡含量高，均值达129.95 × 10 −6 ，最大值达2189 × 10 −6 ，变异系数为2.24；铜、铋、锌含量均较高，变异系数大。表明该区易于富集锡铜铋等元素 [ 5 ] 。 图1. 汛塘锡铜多金属异常剖析图 Table 1 1:50，000 chemical geochemical abnormality parameter 元素 异常面积(km 2 ) 平均值 最高值 异常下限 衬度(CD) 规模值(NAP) Sn 7.83 333 1858 10 33.3 144.3 Cu 5.77 148.37 603 50 2.97 17.1 Ag 5.00 0.54 3.2 0.2 2.70 13.5 As 5.62 217 554 60 3.62 20.3 Bi 3.30 9.47 29.1 2 4.74 15.6 Pb 2.96 233 346 75 3.11 9.2 Zn 5.30 394.7 982 250 1.58 8.4 Sb 3.17 30.1 75.4 10 3.01 9.5 表1. 1:5万化探异常特质参数 含量单位：Au × 10 −9 ，其它×10 −6 。 图2. 汛塘工区土壤异常剖析图 Table 2 Table of parameters of 1:1 million soil mineralization elements in Qiangtan 元素 观测数 最大值 最小值 平均 标准差 变异系数 Cu 104 831 14.6 107.67 164.09 1.52 Pb 104 448 10 78.16 83.74 1.07 Zn 104 2250 19.1 267.5 325.47 1.22 Bi 104 62.3 0.3 5.11 11.56 2.26 As 104 966 6.48 142.94 186.19 1.30 Sn 104 2189 2.55 129.95 291.55 2.24 Ag 104 0.86 0.03 0.16 0.16 0.96 Au 104 47.3 0.37 5.21 6.71 1.29 表2. 汛塘1:1万土壤成矿元素参数表 含量单位：Au × 10 −9 ，其它×10 −6 。 对土壤数据进行R型聚类，谱系图(图3)表明：该区大致有三组关系密切的元素组，第一组为Zn、Pb、Au、Ag；第二组为W、Cu、Bi；第三组为As、Sn、Hg，基本的高温-中低温元素组合特征，揭示了该区成矿要素单一，推测为岩浆接触带热液成矿作用导致 [ 6 ] 。 通过地质踏勘，发现异常区北东部有废弃民窿，经民窿编录，见富锡铜多金属破碎带，刻槽样连续4米原生矿石中铜、锡均达工业品位，铜含量最大值4.5%，锡含量最大值1.66%，锌含量最大值1.41%。 对民窿矿石光片鉴定结果显示，矿石为块状锡石闪锌矿矿石，矿石中主要金属矿物为锡石和闪锌矿，其次为黄铜矿。呈粒状结构，团块状、浸染状构造。锡石形成较早，闪锌矿和少量黄铜矿交代锡石，充填于锡石孔隙，有的锡石呈孤岛状分布于闪锌矿之中。闪锌矿中可见乳滴状黄铜矿分离物。矿床为高–中温热液型，闪锌矿、黄铜矿、方解石组合为后期叠加热液形成。矿石为块状锡石闪锌矿矿石。锡石(Cst)约20%，灰白微褐色，高硬度，非均质性，难磨光，有自形晶。闪锌矿(Sp)约30%，灰色，中硬度，均质性，常有乳滴状黄铜矿。 图3. 汛塘1:1万土壤成矿元素R型聚类谱系图
system: [SYSTEM]

user: 5. 找矿潜力分析 本异常区成矿条件好，找矿潜力大。异常区存在锡、铜异常强度大，且见小型锡铜矿点；南西部位汛塘岩株，燕山期岩浆岩，为成矿提供物质来源 [ 7 ] [ 8 ] ；区域上北东、北西向断裂发育，利于矿液输送与储藏；出露地层主要为石炭系地层，分别为石磴子组、测水组、梓门桥组、船山组，偶见大埔组出露，岩性以灰岩及白云岩为主，为有利成矿的地球化学障，利于矿液的沉淀成矿 [ 9 ] [ 10 ] 。 在锡高值点附近已有民采矿坑、民窿及选矿设备。汛塘矿山地表捡块样铜达4.92%，在马占新发现4 m宽锡铜矿破碎带，表明该区含矿性较好，矿体具北东向延伸特征，为中高温热液充填破碎带型矿床。
system: [SYSTEM]

user: 6. 结论 1) 汛塘地区成矿条件有利，靠近汛塘岩株，区内断裂发育，地层主要为石炭系地层，这一区域的Sn、Cu、Bi、Ag、As、Au、Zn、Pb等化探异常已反映了控矿成矿信息。 2) Sn、Cu为本区主要成矿元素，浓集系数高，变异系数大。主体沿北东向展布，与断裂方向一致。北东向、北西向断裂的发育方向可能指示了成矿和容矿的矿化信息。 3) 汛塘地区有一定的找矿潜力，有进一步开展地质工作的必要。下一步建议加强汛塘原有矿点的外围找矿工作，在新发现矿体附近开展预查工作。
system: [SYSTEM]

user: 摘 要 马泰壕煤矿距离红碱淖最近短距离2 km，正确评价其开采对红碱淖补给量的影响，成为合理开发利用煤炭资源过程中急需解决的问题。在区域及井田地质水文地质调查基础上，运用概率积分模型进行地表变形预测，并根据“三下”规程进行导水裂缝带发育高度计算结果。分析煤矿开采对红碱淖地表、地下产汇流及其含水系统的影响。研究表明：全井田开采后，大的地形趋势没有发生变化，大部分地表水及潜水仍向红碱淖方向汇流，局部因采煤沉陷导致河道水力坡度变化和潜水位埋深变浅，甚至出露。最大减少向红碱淖补给的地表水量只占其天然地表径流量的1.47%。在直罗组稳定隔水层阻隔下，煤矿开采不会对红碱淖含水层中的水造成直接影响。 关键词 :马泰壕煤矿，红碱淖补给量，采煤沉陷，产汇流 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. 引言 20世纪90年代，尤其是1992年以后，我国最大的沙漠淡水湖泊——红碱淖湖面面积不断萎缩。国内外众多学者对其萎缩原因 [ 1 ] - [ 6 ] 进行研究。刘萍萍等 [ 1 ] 通过流域水量平衡计算，指出自然因素是影响红碱淖流域水量失衡的主要因素。李登科等 [ 2 ] 认为20世纪90年代以来红碱淖水面面积在大幅度萎缩是气候变化和人为影响因素综合作用的结果，上游河流截流是主要原因。王亚俊等 [ 3 ] 指出干旱区湖泊对人类活动的干扰最敏感。 煤炭开采作为人类活动的典型代表，国内外诸多学者分别从影响机理 [ 7 ] [ 8 ] [ 9 ] [ 10 ] [ 11 ] 、因素识别 [ 12 ] 、趋势预测 [ 12 ] [ 13 ] 等方面阐述其对水资源的影响。主要包括： 曾庆铭等 [ 7 ] 指出煤矿开采对水资源影响机理主要是煤矿开采产生导水裂隙带，使地表水与地下水连通，从而造成地表水大量渗漏，河川径流量明显减少。张发旺等 [ 8 ] [ 9 ] 认为煤炭开发对水资源的影响机制主要是煤炭开发引起覆岩和地表移动，从而破坏含水层结构和改变水循环规律变化、引起补排关系发生变化，结果造成河川径流量减少，井泉干涸，地下水位下降等。并结合神东矿区采煤塌陷对地表水体的影响分析指出，补给源的影响居于主要地位。张思锋等 [ 11 ] 应用有限分布滞后模型改进了多元回归模型，并建立了煤炭开采与乌兰木伦河流径流量相关关系的代数表达式，研究指出在影响乌兰木伦河径流量变化的诸多因素中，采煤活动占人类活动的77.3%。徐良骥等 [ 13 ] [ 14 ] 根据地面沉陷程度预计煤矿开采对水系破坏趋势，为水系治理保护提供依据。 在众多研究中缺乏煤矿开采对地表水、地下水产汇流，尤其是产汇流条件的改变，进而造成补给源变化的定量研究更少。因此，本文以马泰壕煤矿为例，定量化研究煤矿开采对红碱淖补给资源量可能造成的影响，为马泰壕煤矿合理取、用水提供依据，也为黄河流域水资源规划与合理利用提供技术支撑。
system: [SYSTEM]

user: 2. 研究区概况及研究思路 2.1. 红碱淖流域及马泰壕井田概况 红碱淖流域面积1493.3 km 2 ，多年平均年降水深为356.1 mm，多年平均水面蒸发量为1296.5 mm。多年平均地表水资源量为5447.8万m [ 15 ] ，不重复计算的地下水资源量为4829.7万m 3 。多年平均水资源总量为10，277.5万m [ 15 ] 。 红碱淖流域是一个以红碱淖为中心的四周高、中心低洼的内陆封闭状湖盆盆地，这种特定的地形条件形成了地表径流向红碱淖汇集。除地表水体外，红碱淖主要接受周边地下水排泄补给。因此，无论是煤矿开采对地表水和地下水补给量的影响都关系到湖面面积的消长。 马泰壕井田位于红碱淖东北部，占红碱淖流域面积的4.28%。其井田边界与红碱淖最近距离为2.65 km，开采边界与红碱淖最近距离为3.2 km，地形整体上北高南低，地面最大标高为1398.7 m，最低标高为1223.6 m，最大地形相对高差为175.10 m；一般地形相对高差为90 m左右。 井田设计开采面积63.9 km 2 ，地质资源量140，107万t，主要可采煤层为3-1、4-1、5-1、5-2、6-1、6-2。根据地质填图及钻探成果对比分析，区内地层由新至老发育有第四系(Q)、白垩系下统志丹群(K 1zh )、安定组(J 2a )、直罗组(J 2z )，侏罗系中统延安组(J 2y )，三叠系上统延长组(T 3y )。主要含水层为第四系松散岩类孔隙潜水含水岩组和白垩系下统志丹群、侏罗系直罗组下段、延安组孔隙、裂隙承压含水层。其中厚度67.10~130.00 m侏罗系中统安定组及直罗组上段(140~330 m)岩层组成全区层位稳定、隔水性能良好隔水层。 井田西南界外有札萨克河，由西北流向东南，最后汇入红碱淖，在距札萨克源头23 km处修建了一座总库容为5117万m 3 的札萨克水库，并于2006年5月建成运行。井田内有木独石犁河和哈拉木河穿过，木独石犁河由北向南流入红碱淖自然保护区，哈拉木河是札萨克河的一条支流，在札萨克水库下游汇入红碱淖。 2.2. 研究思路 因目前马泰壕井田没有建成投产，缺乏沉陷量和导水裂隙带发育高度的实测资料。根据东胜矿区多年开采实测资料表明，马泰壕煤矿符合概率积分法 [ 16 ] ，故采用概率积模型进行沉陷量预测。导水裂隙带、冒落带高度计算主要依据“三下压煤开采规程”中“缓倾斜和倾斜煤层开采时导水裂缝带高度计算”和“厚煤层分层开采的垮落带高度计算”经验公式。 H l i = 100 ∑ M 1.6 ∑ M + 3.6 ± 5.6 (1) 式中：H li 为导水裂缝带高度， ∑ M 为累计采厚，公式应用范围为单层采厚1~3 m，累计采厚不超过15 m。 H k = 100 ∑ M 4.7 ∑ M + 19 ± 2.2 (2) 式中：H k 为垮落带高度， 为累计采厚，公式应用范围为单层采厚1~3 m，累计采厚不超过15 m。 根据沉陷预测结果从对井田内地表水产流、河道径流、潜水汇水等方面分析煤矿开采对红碱淖地表、地下产汇流和补给量的影响；根据导水裂隙带发育高度及影响半径计算结果判断煤矿开采对红碱淖湖区含水层的影响。
system: [SYSTEM]

user: 3. 预测结果 3.1. 采煤沉陷预测 上世纪50年代开始，“三下”采煤技术作为一种随着煤炭生产活动发展而诞生的一门新型实用技术，在我国受到关注和重视。在经过了几十年不断的探索、实验与实践后，该项技术从应用范围、核心技术、理论水平等多个方面看，已经成为一个实用性强、应用范围广并兼具相应理论和技术的学科。在多年的水体下采煤实践中，研究人员归纳出了不同覆岩岩性、倾角、采煤方法的上覆岩石破碎规律。根据总结出的上覆岩石破坏规律，相关学者不仅构建出了一套适用不同岩性、倾角、开采方法的冒落带和导水裂缝带高度的计算公式，而且形成了一套地表变形预测理论。该理论被广泛的应用于采煤引起的沉陷预测中。 3.1.1. 预测参数的确定 马泰壕井田内各煤层直接顶板岩性主要为砂质泥岩、细粒砂岩、粉砂岩，次为中粗粒砂岩，抗压强度较低，平均为26.2 MPa。根据“三下规程”中分层岩性评价系数选取表，并结合煤层顶板岩性及开拓开采情况，进行计算参数的选取，详见表1。 Table 1 The estimation of the Surface deformation paramete 参数 符号 参数值 备注 下沉系数 q 0.8 第一次，第二次重复采动分别取0.88，0.95 主要影响正切 tgβ 2.3 重复采动取2.5 水平移动系数 b 0.3 拐点偏移距 S 0.15～0.25H H为平均采深 影响传播角 θ 90－ 0.68a 表1. 地表移动变形预计参数 3.1.2. 预测方案及结果 根据矿区的盘区划分和接续计划，按照“远粗近细”的原则，同时考虑红碱淖作为研究中最重要的保护目标，为说明各煤层开采后对红碱淖的影响，本次将各煤层距红碱淖最近的盘区开采完成后作为一个阶段。因此将首采区(3-1煤一盘区)作为沉陷预测的第一阶段。全井田开完最为最后一个阶段。 根据上述预测方案计算得出：3-1、4-1、5-1、5-2、6-1、6-2煤依次开采后地表沉降量为4.96 m、6.62 m、7.67 m、8.90 m、10.80 m、12.99 m。首采区和全井开采结束后不同地表沉陷深度影响面积，详见表2。 Table 2 Along with the exploitation of the first mining area and the mine field， the land surface has subsided for different degrees， which may affect different field 沉陷深度范围 首采区 全井田开采结束 影响面积(km 2 ) 占总影响面积比例 影响面积(km 2 ) 占总影响面积比例 0.010~2 m 6.39 17.74% 7.52 12.92% 2~4 m 4.05 11.24% 3.31 5.69% 4~6 m 25.59 71.02% 3.24 5.57% 6~8 m — — 3.25 5.58% 8~10 m — — 2.65 4.55% 10~12 mm — — 14.23 24.45% 12 m以上 — — 24 41.24% 合计 36.03 100% 58.2 100.00% 表2. 首采区和全井开采结束后不同地表沉陷深度影响面积统计表 由表2可知：3-1煤所有盘区开采后大部分沉陷区下沉值均在4~6 m之间，占到总沉陷面积的71.02%，不会改变井田内的地形地貌，下沉状况不明显。随着各煤层的开采，地表下沉值和沉陷面积逐渐增大，全井田开采结束后地表下沉值达到最大的12.99 m，41.24%的沉陷范围下沉值在12 m以上，下沉深度在10 m以上的区域占总沉陷面积的65.69%，因此全井田开采结束后地表会形成明显的下沉盆地。马泰壕现有地下水位埋深1.0~50.0 m左右，一般在5.0~10.0 m，煤矿开采将导致沉陷区内浅层地下水位埋深变浅甚至出露于地表形成积水盆地，盆地的大小与开采煤层数量有关。 3.2. 导水裂隙带高度预测 井田主采煤层为3-1煤层，煤厚4.76~7.63 m，平均6.32 m。顶板岩性为砂质泥岩，底板岩性以砂质泥岩、泥岩为主，其次为粉砂岩。根据钻孔岩石物理、力学性试验结果表明，岩石的抗压强度多数均小于30 MPa，个别在30~60 MPa之间。井田内煤层顶底板岩石以软弱~半坚硬岩石为主，导水裂缝带高度选用覆岩岩性为“中硬”的计算公式。详见表3。同时根据井田内第1、3号走向线和2、5、6、13、14号勘探线的各钻孔煤层分别进行导水裂隙带高度计算结果进行验证。 Table 3 Predication results of the height of water-flowing fractured zone after 3-1seam minin 煤层 最大采厚(m) 导水裂隙带高度(m) 冒落带高度(m) 3-1 6.2 65.24 16.11 表3. 3-1煤层开采后导水裂缝带发育高度 由表3可知：导水裂隙额发育最高高度为3-1煤层之上的65.24 m处，只穿透延安组含水层，并发育至侏罗系中统直罗组(J 2Z )下段含水层，不会导通侏罗系中统直罗组(J 2Z )上段隔水层，也不会发育到安定组隔水层，顶点距白垩系下统志丹群底界面最近距离为83 m，距第四系含水层底部最近距离为275 m，白垩系与第四系地层只是随着下部地层的下沉而产生整体沉降。
system: [SYSTEM]

user: 4. 采煤沉陷对红碱淖地表水及地下水产汇流影响 4.1. 采煤沉陷对井田河流产流的影响 4.1.1. 沉陷对札萨克河产流的影响 札萨克水库坝址以上河道长23 km，占札萨克河道总长的67.9%；流域面积234 km 2 ，占札萨克全流域面积的65.7%，札萨克水库坝址以上3.4 km 2 的集水面积在井田内，占井田面积的5.3%，自产天然径流量约为11.9万m 3 。由于札萨克水库全部拦截了上游的来水量，该区域由于地表沉陷全部不产流情况下，只减少札萨克水库的供水量，不会减少札萨克河向红碱淖的补水量。 井田位于札萨克水库下游，开采范围内有46.1 km 2 的札萨克河集水面积(包括哈拉木河)，约占札萨克河全流域面积的12.9%，该区域自产天然径流量约为162.1万m 3 。2000年以后经济社会的用水消耗、水库工程的建设影响、植被用水增加，以及地下水开采的消耗导致的减水量为81.2万m 3 ，也即该区域现状实际向札萨克河补水量为80.9万m 3 。即使该区域由于地表沉陷全部不产流，最大减少的红碱淖补给量80.9万m 3 。 4.1.2. 沉陷对木独石犁河产流的影响 木独石犁河全流域面积35.9 km 2 。在流入马泰壕井界上游260 m处由东胜区宝恒电厂修建1个水库，该水库控制流域面积15.2 km 2 ，自产天然径流量为57.9万m 3 ，主要解决电厂生活及工业用水，基本将其范围内木独石犁河径流量全部拦截。因此，马泰壕井田开采造成的地表沉陷仅对木独十里河在井田范围内地表径流产生影响。 木独石犁河在井田内河段长度7 km。且河谷切割较深，在矿井范围内河宽约为3~5 m，深约2~3 m，井田内的积水面积自产天然径流量为55.1万m 3 。采煤沉陷将其全部截流，减少的向红碱淖补给量最大为其自产天然流量55.1万m 3 。由于当地企业及群众在木独石梨河修建宝恒电厂水库并沿河挖渗坑取水，2000年后已经不向红碱淖补水。 4.2. 采煤沉陷对河道径流的影响 木独石犁河和哈拉木河，河道比降分别为5.98‰和7.64‰，木独十里河进出井田原始河道标高高差为64 m [ 17 ] ，哈拉木河进出井田原始河道标高高差为89 m [ 17 ] ，详见图1、图2。远大于河道下沉深度0~12.99 m，整个河道径流流向不会发生改变。煤矿开采过程中河道局部沉陷深度不同，局部河段水力坡度发生变化，一方面在局部河段将会出现河道水面扩大的现象，从而增加河道的水面蒸发，以及向潜水的补给。另一方面由于采煤沉陷，地下水位下降，增加河道渗漏量，减少地表径流向红碱淖的补给量。 图1. 开采前后哈拉木河河道变化图 图2. 开采前后木独石犁河河道变化图 4.3. 采煤沉陷对地下水汇水的影响 当地表发生大面积的沉陷时，安定组及直罗组隔水层也随之产生沉陷，该隔水层虽然隔水性能良好，能够隔断2个含水层向下部采空区渗漏，但是由于地表沉陷，必然在采动影响区和非影响区产生一定落差，这样就会使非采动区潜水以下沉区域为排泄区，源源不断的流入到下沉区域内，以抬高下沉区域下降的潜水面。随着重复采动的不断进行，地表下沉幅度逐渐增大，必然会进一步导致开采区域及其周边的潜水水位整体下降。但在更远范围地下水共同补给下，这种影响将会逐渐消除。因此，在马泰壕矿区地势北高南低、东高西低、最大地形标高差167.7 m、一般标高差90 m的地形条件下，开采引起的地表沉陷只会在局部地段将之前马泰壕井田内潜水直接朝红碱淖方向汇流，转变为先向沉陷区汇流，仅减慢了马泰壕井田内潜水向红碱淖的径流速度，并不会改变潜水向红碱淖的汇流方向。
system: [SYSTEM]

user: 5. 采煤对红碱淖含水层的影响 5.1. 采煤对红碱淖浅部含水层的影响 红碱淖区域与马泰壕井田含水层如表1所示。由表4可知：红碱淖区域与马泰壕井田具有相同的含水层是第四系松散层潜水，下白垩系砂岩裂隙、孔隙潜水。相异的含水层均属于侏罗系地层，其中红碱淖区域内的安定组含水层位于马泰壕井田内直罗组底至3煤层之间含水层之上，两者之间存在直罗组隔水层。 Table 4 The contrast of the comparability and differences between the Hongjiannao region and Ma Tai Hao well-field aquifer 红碱淖区域 马泰壕井田 备注 第四系松散层潜水含水层 第四系松散层潜水含水层 不受开采扰动影响 下白垩系洛河组砂岩裂隙、孔隙潜水(马泰壕井田内缺失) 白垩系下统志丹群孔隙潜水～承压水含水层 不受开采扰动影响 侏罗系安定组砂泥岩裂隙潜水(可视为隔水层) 侏罗系中统碎屑岩类承压水含水层 受开采扰动影响 直罗组下段裂隙承压水含水层 受开采扰动影响 表4. 红碱淖区域与马泰壕井田含水层对比表 在侏罗系地层中，红碱淖区域内的安定组潜水含水岩组虽然单独划分为1组含水层，但实质上该含水层只有上部的风化带含有裂隙水，这部分裂隙水赋存在风化砂泥岩中，其含水量是极其微弱的，可视为隔水层。因此在侏罗系地层中，红碱淖区域内实际上没有含水层，也不存在与马泰壕井田具有相同层位的含水层。 根据导水裂隙带发育对上覆地层影响结果可知：煤矿开采时上部2个含水层基本没有影响。如这两个区域内的潜水含水层为等厚均质连续体，且存在水力联系，在稳定隔水层阻隔作用下，从含水层受威胁程度讲，煤矿开采不会直接影响到红碱淖区域内的潜水含水层。 5.2. 采煤对红碱淖深部承压含水层的影响 煤层开采后，井田范围内的侏罗系直罗组下部含水层和煤系地层中的延安组含水层的水将全部被疏干，同时，周边未开采区域内的同层含水层水位也随之产生下降。井田开采对周边同层含水层影响半径采用下式计算，结果详见表5。 R = 10 S K (3) Table 5 Influnce radius of the coal mining on perimeter aquifer 含水层名称 渗透系数(m/d) 水位标高(m) 水位降深值(m) 影响半径(m) 直罗组下部含水层与延安组含水层 0.00362748~0.0265 1256.832~1270.516 316.832~330.516 190.8~538.0 表5. 煤矿开采对周边含水层影响半径 由表5可知：直罗组下部含水层与延安组含水层在被完全疏干的情况下，对周边的最大影响半径约为538.0 m。
system: [SYSTEM]

user: 6. 结论及建议 1) 井田开采形成的采煤沉陷，导致该区域地表水完全不产流情况下，最大减少向红碱淖补给的地表径流量为80.1万m 3 ，占红碱淖天然地表径流量5447.8万m 3 的1.47%。 2) 木独石犁河和哈拉木河进出井田界河床高程差远大于河道下沉深度0~12.99 m，河道径流流向不会发生改变，洪水期及时采取疏排措施，基本不会对河道汇流产生影响。 3) 由于地表存在较大的地势差，开采不会改变马泰壕潜水向红碱淖的汇流方向，仅减慢了径流速度。若地表不形成积水盆地，开采也不会影响马泰壕井田潜水补给红碱淖的径流量。因此，矿井开采过程中，应避免在地表形成积水盆地，最大限度地减少潜流转明流而带来的水面蒸发。 4) 在直罗组隔水层隔水作用下，煤矿开采不会直接影响到红碱淖区域内的潜水含水层。马泰壕采煤对周边的最大影响半径约为538.0 m，而井田开采边界与红碱淖边界之间的最近距离为3.2 m。因此，煤矿开采也不会波及到红碱淖湖域范围内的直罗组下部含水层与延安组含水层。
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/