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

user: 4. 掌子面及掌子面超前岩体挤出变形跟踪监测 和新奥法一样，新意法也要求在隧道开挖和支护施作期间对隧道受力变形进行跟踪监测，根据监测结果对设计施工参数进行动态调整，确保隧道安全和工程质量。除了新奥法要求的掌子面后方隧道收敛变形观测及初始衬砌的受力变形量测外，新意法特别重视掌子面的挤出变形量测，以及掌子面超前岩体的挤出变形量测。 本工程对于掌子面挤出变形量测采用全站仪完成，对于掌子面超前岩体挤出变形量测采用了本工程专利技术——基于磁环感应原理的隧道掌子面前方深部岩体挤出变形的量测技术。 4.1. 掌子面挤出变形实测 掌子面挤出变形量测要求在每次爆破清渣后即布设监测点，按每4小时间隔频率进行量测。本工程试验段掌子面挤出变形测点布置如图11所示，该测量断面位于桩号K54+246处。 量测得到出渣后8小时及12小时时刻各点挤出变形，如表2所示。将表2数据绘于掌子面上，即得挤出变形在掌子面的分布，如图12所示。实线和虚线分别为出渣后12小时和8小时数值。 从表2和图12可以看出，掌子面挤出变形的分布规律是：中央数值最大，边缘数值最小；挤出变形随时间增大，但增速不大。 图11. 掌子面挤出变形测点布置图(图中尺寸单位：cm) Table 2 Extrusions of the monitoring points at workface of tunnel (unit: mm 测点 1 2 3 4 5 出渣后8小时 14.1 15.5 14.3 14.8 14.8 出渣后12小时 15.7 17.8 17.0 18.1 17.7 变形增量 1.6 2.3 2.7 3.3 2.9 表2. 掌子面挤出变形一览表(单位：mm) 图12. 掌子面挤出变形(图中尺寸单位：mm) 4.2. 掌子面超前岩体挤出变形实测 本工程对于掌子面超前岩体挤出变形量测采用基于磁环感应原理的隧道掌子面前方深部岩体挤出变形的量测技术。其原理近似于岩土工程的分层沉降量测技术，只是将垂直沉降量测变为水平位移量测而已。 从掌子面水平钻孔进入掌子面超前岩体中，在掌子面前方不同位置埋设磁环，磁环套接于PVC管上，PVC管内部中空。PVC管外注水泥浆，确保磁环和岩体共同位移。另有可以探测磁环位置的电磁感应仪安装在有刻度的测杆上。电磁感应仪电源开关开启后，电磁感应探头通过磁环时会发出蜂鸣声。根据蜂鸣时测杆读数，即可以获得蜂鸣时磁环的位置。测试示意图如图13所示。 图13. 基于磁环感应原理的掌子面超前岩体挤出变形量测技术示意图 量测步骤如下：先通过全站仪确定孔口的水平位移，然后利用磁环量测计量测得到磁环与孔口的相对位移，最后两者求和，得到磁环的绝对水平位移。 本工程试验段在桩号K54+246处布设了磁环感应计，如图14所示。 图14. 掌子面超前岩体挤出变形测点布置图 量测得到爆破出渣后8小时及12小时时刻，各测点的挤出变形值，如表3所示。如图15所示。图中虚线和实线分别表示爆破出渣8小时和12小时的数值。 Table 3 Extrusions of the monitoring points in front rock (unit: mm 测点 0 1 2 3 4 5 6 与掌子面距离(m) 4 6 8 10 12 15 测孔 SA 15.7 4.7 2 0 0 0 0 SB 17.8 6.3 2.5 0 0 0 0 SC 17.0 4.0 1.5 0 0 0 0 表3. 掌子面超前岩体挤出变形一览表(单位：mm) 图15. 掌子面超前岩体挤出变形 从表3和图15可以看出，掌子面中央挤出变形最大。距离掌子面6米处的挤出变形小于5mm，至8米处已无挤出变形，说明开挖影响半径在6米至8米之间。从图中还可以看出，掌子面挤出变形在前后两次的增量不大，超前岩体的挤出变形的增量衰减更快。可见玻璃纤维锚杆加固效果是明显的，掌子面是稳定的和安全的。
system: [SYSTEM]

user: 5. 三维有限元数值分析 5.1. 三维有限元计算模型 本工程试验段在确定玻璃纤维密度、搭接长度和初始衬砌参数的优化，均通过三维有限元数值计算来实现，所用计算分析软件为Z-Soil软件。 Z-Soil软件是瑞士联邦理工学院(Swiss Federal Institute of Technology in Lausanne)开发的真三维岩土工程有限元软件，可以提供涉及土力学和岩石力学、隧道工程、基坑开挖、土–结构相互作用、地下水和温度分析的统一方法，在隧道工程界得到广泛应用。 图16为三维有限元计算模型。计算中围岩本构模型采用摩尔一库仑(M-C)模型，玻璃纤维锚杆采用锚杆单元模拟。隧道的初期支护按弹性材料考虑，采用壳单元模拟。计算不考虑二次衬砌的影响。计算所用的玻璃纤维锚杆力学参数如表4所示。计算所用的材料力学参数如表5所示。 图16. 三维有限元计算模型 Table 4 Properties of GFRP bol 弹性模量E(GPa) 锚杆直径(mm) 水泥浆粘聚力(kPa) 内摩擦角φ(˚) 40 22 20 × 10 4 25 表4. 玻璃纤维锚杆力学参数 Table 5 Properties of material for 3D-FEM analysi 材料 容重γ(kN/m 3 ) 弹性模量(E/GPa) 泊松比μ 内聚力c(kPa) 内摩擦角(˚) IV围岩 21.5 3 0.3 100 25 V围岩 19 1 0.4 70 20 超前核心岩体(注浆加固) 19 1.1 0.35 75 21 初始支护 25 23 0.20 表5. 三维有限元数值分析的材料力学参数 5.2. 三维有限元数值计算结果 5.2.1. 掌子面超前岩体的挤出变形云图 计算得到掌子面超前岩体挤出变形云图，如图17所示。掌子面最大挤出变形在隧道中央，为17.5 mm。从图中可以看出，掌子面前方岩体的挤出变形衰减很快，在距离掌子面4 m处挤出变形已下降至6 mm以下，距离掌子面6 m处，几乎无挤出变形。说明，玻璃纤维锚杆锚杆加固后，掌子面稳定，采用全断面开挖是可行的。 图17. 掌子面超前岩体挤出变形云图(单位：mm) 5.2.2. 有限元计算结果和实测结果比较 针对掌子面超前岩体的挤出变形，三维有限元计算结果和实测结果对比分析如图18所示。从图中可以看到，计算结果和实测结果变化规律是一致的，量值也比较接近。说明采用的计算模型和力学参数是合理的。依此制定的施工参数是合理的。 图18. 掌子面超前岩体挤出变形计算值和实测值的比较
system: [SYSTEM]

user: 6. 新意法与新奥法施工效果的对比分析 新奥法采用上下台阶法进行隧道开挖，其具体工艺是：开挖上台阶，对拱部施工系统锚杆和钢拱架及挂网喷射混凝土；开挖下台阶，完成侧墙钢拱架及挂网喷射混凝土，初始支护完成封闭；最后开挖仰拱，施工仰拱结构。新意法采用全断面进行隧道开挖，其具体工艺是：掌子面打设超前玻璃纤维锚杆，全断面开挖，施工系统锚杆和钢拱架及挂网喷射混凝土，初始支护一次封闭。前者属于分部开挖，后者属于全断面开挖。为了比较二者工艺的效果，分别对有限元计算结果和实测洞周位移进行了比较。 6.1. 新意法和新奥法有限元计算结果比较 通过三维有限元计算，分别得到两种方法下的隧道变形情况。新意法计算结果见本文前节。 图17和图19分别为新意法和新奥法施工时隧道轴向变形云图。 图19. 新奥法时隧道的轴向挤出变形云图，单位：m 从图19中可以看出，采用新奥法当上台阶开挖后，掌子面的轴向变形为33.61 mm，下台阶开挖后掌子面的轴向变形继续增大，达到54.21 mm。从图17中得到的新意法开挖后掌子面的轴向挤出变形为17.55 mm。计算结果表明，新意法对于控制掌子面变形是显著的。 6.2. 实测结果比较 根据现场监测结果，对采用新奥法和新意法施工的两段隧道的收敛变形进行了比较。表6为新意法和新奥法施工时拱顶沉降和收敛位移实测值。 Table 6 Comparison of measurement between tunnels by NATM and by NIT 施工方法 新奥法 新意法 掌子面位置 ZK54+310 ZK54+250 收敛位移(mm) 上测线 27.9 24.9 下测线 24.4 19.5 拱顶沉降(mm) 18.2 15.9 表6. 新奥法和新意法实测值结果比较 从表6可以看出，新意法施工时隧道的收敛位移和拱顶沉降均比新奥法施工时来得小，说明新意法对于掌子面后方隧道洞周变形的控制比新奥法有效。
system: [SYSTEM]

user: 7. 结论 1) 对于通过强度较差的IV围岩和V围岩的大断面隧道，采用新意法可以实现隧道的全断面开挖； 2) 新意法施工实践表明：全断面开挖、掌子面超前岩体加固、掌子面超前岩体挤出变形实时监控是控制隧道变形及掌子面稳定的关键因素；计算和实测结果均表明新意法比新奥法具有优势。 3) 玻璃纤维锚杆应用于掌子面前方岩体的超前加固方面效果可靠，方便开挖，是一种值得推荐的超前加固方法。
system: [SYSTEM]

user: 摘要 裂隙在数理方法层面可表达为渗流力学边界条件。对于裂隙环境下的弹性波传播，在裂隙处设置法向位移为零的渗流力学边界条件，与实际渗流现象更加符合。对这一数理表达进行求解，不仅能够生动描述弹性波穿透裂隙传播的效应，而且能够有效反映裂隙长度、位置、端点等因素对弹性波传播的共同影响。
system: [SYSTEM]

user: 1. 引言 渗流是指流体渗入物质裂隙过程中形成的微观流动作用 [ 1 ]。裂隙在自然界普遍存在，大到可能发生于受到力学作用的岩土层，小到可以产生于增材制造等工艺下制备出的材料或元件。源于裂隙的渗流现象会改变物质原有的物理特性，对资源勘探开发、物质无损检测等技术成果的可解释性 [ 2 ] 具有广泛影响。 弹性波在物质中的传播机制是以上技术领域常常需要面对和解决的力学基本问题。当弹性波传播到物质上的裂隙时，受到渗流力学作用的影响，传播方向会发生偏离，波动性态也会有所改变 [ 3 ] [ 4 ]。为了准确描述上述现象，有必要聚焦于弹性波传播问题，构建一套行之有效的渗流力学数理方法。本文假设裂隙是作为渗流力学边界条件而存在的，首先在弹性波方程的基础上给出边界条件的若干具有直观意义的数理表达，然后绘制几种典型的渗流模型，采用频域有限元方法模拟弹性波传播到裂隙处时发生的奇特现象，揭示各种渗流力学作用对弹性波传播的影响。
system: [SYSTEM]

user: 2. 理论方法与基本假设 2.1. 带有完全匹配层的弹性波方程的等效积分弱形式 不失一般性，考虑二维直角坐标系(x， y)下无限大空间内的弹性波传播问题。在频域下引入完全匹配层(PML)，以有限空间区域表达无限大空间内的弹性波传播。推导弹性波方程的等效积分弱形式(EIWF)，以适用于频域有限元求解 [ 5 ]。EIWF的数理表达为 { ∫ L ( s y τ x x n ⇀ δ u + s x τ x y n ⇀ δ v ) d L − ∬ S ( s y τ x x ∂ δ u ∂ x + s x τ x y ∂ δ v ∂ y ) d S + ∬ S s x s y ρ ω 2 u δ u d S = 0 ∫ L ( s y τ x y n ⇀ δ u + s x τ y y n ⇀ δ v ) d L − ∬ S ( s y τ x y ∂ δ v ∂ x + s x τ y y ∂ δ u ∂ y ) d S + ∬ S s x s y ρ ω 2 v δ v d S = 0 (1) 其中，S为一有限空间区域，L泛指区域内外边界， n ⇀ 为边界外法向，ρ为物质密度，ω为弹性波角频率，u、v分别为弹性波在x、y分量的位移，τ xx 、τ yy 分别为弹性波在x、y分量的法向应力，τ yy 为切向应力。另，s x 、s y 与PML有关，分别表征弹性波在x、y分量的吸收，具体设置方式详见文献 [ 6 ]。 2.2. 对裂隙的基本假设 如果物质表现为纯粹的固态，那么弹性波在物质上任一位置处的传播均满足法向位移连续、切向位移连续、法向应力连续、切向应力连续这四个边界条件；如果物质表现出固态和流态的接触作用，那么弹性波在接触处的传播满足法向位移连续、法向应力连续、切向应力为零这三个边界条件 [ 7 ]。据此，本文对物质上的裂隙性质作出如下假设： ① 裂隙宽度趋于零； ② 裂隙存在与否对弹性波传播具有实质影响，而裂隙内填充的物质对弹性波传播的影响不大； ③ 弹性波传播到裂隙处时，部分或全部法向连续性边界条件丧失； ④ 为确保弹性波穿透裂隙，切向连续性边界条件不作处理。 2.3. 边界条件的数学表达 根据章节2.2提出的假设，尤其是③号假设，给出渗流力学边界条件的三种具有直观意义的数理表达。 2.3.1. 弱渗流裂隙 弱渗流裂隙处满足的边界条件有两种情形： ① 法向位移为零、法向应力连续。此时有 n ⇀ ⋅ [ u L 0 v L 0 ] T = 0 (2) 其中，L 0 为裂隙所在位置。 ② 法向位移连续、法向应力为零。此时有 n ⇀ ⋅ τ L 0 = 0 (3) 其中， τ 为应力矩阵， n ⇀ ⋅ τ 代表法向应力。 2.3.2. 强渗流裂隙 强渗流裂隙处同时满足法向位移为零、法向应力为零这两个边界条件。此时式(2)和式(3)共同成立。
system: [SYSTEM]

user: 3. 应用 3.1. 裂隙对弹性波传播的影响 3.1.1. 模型与参数 如图1所示绘制1号模型。物质长1.2 m，宽0.8 m，外围施加一厚度为0.2 m的PML。一条无限长的竖向裂隙将物质分割为左右对称的两部分。在1号模型下，渗流力学边界条件简化为： ① 弱渗流裂隙，法向位移为零、法向应力连续。此时根据式(2)，有 u L 0 = 0 (4) ② 弱渗流裂隙，法向位移连续、法向应力为零。此时根据式(3)，有 τ x x L 0 = 0 (5) 在有限元求解框架下，式(5)化为 ρ v p 2 ∂ v x L 0 ∂ x + ρ ( v p 2 − 2 v s 2 ) ∂ v z L 0 ∂ z = 0 (6) 其中，v p 、v s 分别为弹性波在物质中传播的纵波、横波速度。 ③ 强渗流裂隙，同时满足法向位移为零、法向应力为零。此时式(4)和式(6)共同成立。 在模型1中，点源和接收端处于同一水平面。点源到裂隙的距离为0.5 m。接收端一点在裂隙上，另一点在点源另一侧，到裂隙的距离也为0.5 m。 图1. 1号模型 物质的弹性参数如表1所示。 Table 1 Elastic parameters of a materia 参数 数值 密度(kg/m 3 ) 2000 纵波速度(m/s) 3000 横波速度(m/s) 1500 表1. 物质的弹性参数 3.1.2. 结果与讨论 在点源处激发一个频率为10 kHz的u分量冲击函数，考察三种边界条件数理表达下1号模型内的u分量弹性波传播响应(以实部为例)，分别如图2(b)、图2(c)、图2(d)所示。为便于比较，图2(a)展示了裂隙不发生作用的情形。同时，表2给出了以上四种情形下接收端处的u分量响应值。分析发现，法向应力为零的边界条件对弹性波传播响应的影响不大，图2(c)的波环相比图2(a)只是多了一些抖动，弹性波几乎毫无障碍地穿透到了裂隙另一侧，无法合理反映渗流力学作用；但是，法向位移为零的边界条件，会使得弹性波传播到边界条件处时产生的特性变化尤为明显，即图2(b)的波环形态(如幅度、波数等)相比图2(a)发生了更为明显的变化，这与渗流力学实际现象更加符合。当然值得注意的是，将裂隙处的法向位移强制取零，并不必然代表弹性波法向位移无法穿透裂隙传播，这可能是因为部分切向位移也在裂隙处发生了传播方向的偏转，“衍生”了法向位移，并在裂隙另一侧传播；此外，根据式(1)，各个位移、应力分量的相互耦合作用也促进了这一穿透效应的产生。 图2. 1裂隙对10 kHz弹性波传播响应的影响(模型内)。(a) 无裂隙；(b) 弱渗流裂隙①；(c) 弱渗流裂隙②；(d) 强渗流裂隙 Table 2 Effect of a crack on the elastic propagation response with a frequency of 10 kHz (At the receiver 裂隙类型 裂隙处接收端/(10 −12 ) 物质内接收端/(10 −12 ) 无裂隙 −1.2711 − 2.7712i 2.5190 + 0.4845i 弱渗流裂隙① 0.0000 + 0.0000i 0.2910 + 0.0973i 弱渗流裂隙② −1.2188 − 3.1541i 2.3614 + 0.6462i 强渗流裂隙 0.0000 + 0.0000i 0.2582 + 0.2264i 表2. 裂隙对10 kHz弹性波传播响应的影响(接收端处) 3.2. 裂隙端点对弹性波传播的影响 3.2.1. 模型与参数 考虑一个更复杂些的模型。如图3所示绘制2号模型。物质长、宽、外围施加的PML厚度不变。模型内部有一条有限长的斜向45˚裂隙。鉴于章节3.1得到的重要结论，在2号模型下，仅考虑弱渗流裂隙①的边界条件，此时根据式(2)，有 u L 0 cos θ + v L 0 sin θ = 0 (7) 当 θ = 45 ∘ 时，式(7)简化为 u L 0 + v L 0 = 0 (8) 若执意要考虑法向应力为零的边界条件，则根据式(3)，有 τ 45 ∘ L 0 = τ x x L 0 cos 2 θ + τ y y L 0 sin 2 θ − 2 τ x y L 0 cos θ sin θ = 0 (9) 当 θ = 45 ∘ 时，式(9)简化为 τ x x L 0 + τ y y L 0 − 2 τ x y L 0 = 0 (10) 然后将各应力分量转换为位移即可，在此不再展开。 在模型2中，点源依旧与两个接收端处于同一水平面。裂隙一端在点源正上方，点源到该端的距离为0.5 m；裂隙另一端在最右侧接收端正下方，距离也为0.5 m。为了考察裂隙端点对弹性波传播的影响，以下将裂隙两端也视作接收端。 物质的弹性参数同如表1所示。 图3. 2号模型 3.2.2. 结果与讨论 在点源处激发一个频率为10 kHz的u分量冲击函数，考察弱渗流裂隙①的数理表达下2号模型内的弹性波传播响应(以实部为例)，分别如图4(a)~(b)所示。为便于比较，读者也可将图2(a)展示的裂隙不发生作用的情形拿来作为对比。分析发现，裂隙端点附近的波环发生了一定程度的扭曲，这就如同文献 [ 4 ] 所提到的，裂隙端点“似乎成为了一个新的‘波源’”，会“源源不断地在模型中传播某种波”。此外，当裂隙不发生作用时，根据点源处激发u分量冲击函数时的并矢格林函数解 [ 8 ]，点源左右两侧、上下两侧的v分量本应呈现反对称特征，然而，倾斜裂隙的存在、裂隙端点的波源效应，共同使这一反对称特征受到了破坏。 图4. 裂隙端点对10 kHz弹性波传播响应的影响(模型内)。(a) u分量；(b) v分量 表3给出了上述情形下各接收端处的u、v分量响应值。其中，最右侧接收端的u分量由于受到裂隙长度、位置、端点等因素的共同影响，其响应值−0.2483 − 0.3743i与1号模型中相应的0.2910 + 0.0973i 相差较大，且发生了极性翻转。顺带地，表3还给出了叠加分量 1 2 ( u + v ) 的响应值。其中，裂隙及其端 点处的叠加分量均为零，验证了模型求解的自圆性。 Table 3 Effect of crack end point on the elastic propagation response with a frequency of 10 kHz (At the receiver and the end point 裂隙类型 裂隙处接收端/(10 −12 ) 物质内接收端/(10 −12 ) 裂隙上端点/(10 −12 ) 裂隙下端点/(10 −12 ) 无裂隙 −1.2711 − 2.7712i 2.5190 + 0.4845i − − 弱渗流裂隙① u分量 −0.6367 − 1.3855i −0.2483 − 0.3743i 4.8803 + 1.0155i −1.3954 − 0.6770i 弱渗流裂隙① v分量 0.6367 + 1.3855i 1.2998 + 0.4108i −4.8803 − 1.0155i 1.3954 + 0.6770i 弱渗流裂隙① 1 2 ( u + v ) 分量 0.0000 + 0.0000i 0.7435 + 0.0258i 0.0000 + 0.0000i 0.0000 + 0.0000i 表3. 裂隙端点对10 kHz弹性波传播响应的影响(接收端、裂隙端点处)
system: [SYSTEM]

user: 4. 结论 本文聚焦于弹性波传播问题，对渗流力学边界条件的数理表达进行了探索，主要认识有： (1) 将裂隙作为边界条件进行处理，不失为探究其渗流力学作用的有效策略； (2) 对于弹性波传播问题，法向位移为零的渗流力学边界条件与实际现象更加符合； (3) 裂隙处的法向位移为零，并不必然代表弹性波法向位移无法穿透裂隙传播。 同时，结合弹性波传播响应在不同渗流模型中的差异可知，本文提出的渗流力学边界条件的数理表达，能够合理反映裂隙长度、位置、端点等因素的共同影响。
system: [SYSTEM]

user: 摘 要 民勤盆地水文与湿地演化对当地生态与沙漠化具有重要作用。通过文献、遗迹及荒漠湿地样方调查和分析得出：历史上，民勤盆地水文与湿地演化分为自然水系、半自然水系和人工水系演化过程，从大面积尾闾湖泊、河流、沼泽、湿地到干涸与沙漠化发展，并随水文演化，植被也经历了从水生、沼泽、湿地逐步向中生、沙旱生、强旱生、超旱生植物演化，沙漠化扩展至全境，是长期气候干旱化和人类活动的结果。近年来，随着全流域水资源统筹，民勤盆地来水量增加，在局部地下水位上升与人工补给下，荒漠区形成了一些微域性湿地，表现在湿地区植被覆盖度大，生态正向发展，随着离湿地距离增加，地下水位下降，植被盖度渐低，土地沙漠化仍很严重，生态逆向发展。 关键词 :水文湿地，沙漠化，演化，民勤盆地 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. 引言 甘肃河西地区深居内陆，南部祁连山高山、亚高山的冰雪冻土带与植被带降水量大于蒸发量，是内陆河流域径流源区，平原盆地绿洲带和荒漠带蒸发量远大于降水量，是径流利用区，各级径流从山区到平原汇集成河流、尾闾湖和湿地，调节气候和生态，也是我国重要的绿洲农业生产基地。然而，随着历史上气候干旱加强、人口增加与水资源不合理利用，一些河流、湖泊和湿地的消失，形成了河西内陆河流域水资源短缺、生态环境恶化与沙漠化不断发展的局面 [ 1 ]。石羊河流域是河西三大内陆河流域中人口最多，产流量最小的内陆河，从历史上水丰草美的草原，绿洲农业大发展，到基本生态用水散失，植被衰退，沙漠化蔓延的过程，已演变成为我国沙漠化最为严重的典型区域 [ 2 ] [ 3 ]。近年来，通过流域统筹用水，生态用水补给，局部湿地恢复，湿地及周边生态有所改善，而外围土地沙漠化仍在发展的情况 [ 4 ]。本文通过石羊河流域中下游民勤盆地水文与湿地历史演化和当代进展研究，以便对流域发展有所借鉴。
system: [SYSTEM]

user: 2. 研究区自然概况与研究方法 自然概况：民勤县地处河西走廊东段石羊河流域中下游，地理位置在N38˚3'45''~39˚27'37''，E101˚49'41''之间，东西长206 km，南北宽156 km，面积1.58 × 10 4 km 2 ，平均海拔1400米，由沙漠、低山丘陵和荒漠平原三种基本地貌组成，绿洲面积10%，其余为荒漠；属温带大陆型干旱荒漠气候，东西北三面被腾格里、巴丹吉林沙漠与河西走廊流沙包围，年均降水量为115.8 mm，年均蒸发量2623 mm，年均气温8.3℃，日照时数为3073.5 h，无霜期162天，年均风速2.4 m/s，最大风速25 m/s；土壤为风沙土、灰棕漠土、棕漠土、绿洲灌淤土为主；自然植被稀疏，以沙旱生、超旱生荒漠植物为主，分布有白刺、柽柳、沙拐枣、梭梭、沙蒿、绵刺、胡杨、沙米、盐生草、碱蓬、盐爪爪等典型天然荒漠植被。 研究方法：通过文献、遗迹调查及居民访问等掌握水系与湿地演化情况；对民勤县境内现有河流、河漫滩，河岸湿地及湖泊、沼泽等湿地进行实地调查与绿洲边缘沙漠化遥感分析，并设置样方(1 × 1 m或5 × 5 m)进行湿地植物种类、株数、植被覆盖度、生长状况、土地类型调查、地上(下)水位监测，进行历史上民勤盆地水文、湿地及其植被演化分析。
system: [SYSTEM]

user: 4. 讨论与结论 民勤盆地水文与湿地演化分为自然、半自然和人工水系的演化过程，从占据多数水域、沼泽、湿地的盆地，变为现在90%的沙漠、戈壁和荒漠化土地与仅10%左右的绿洲，由于各时期气候和人地关系不同，演化历史各具特点 [ 24 ] - [ 29 ]。 民勤盆地从历史上大面积水域、湿地向沙漠化土地不断发展的过程，植被也从水生、沼泽、湿地草甸植被逐步向中生、沙生、沙旱生、强、超旱生植被演变的过程。总体上，湿地区植被覆盖度大，生态正向发展，沙漠化被遏制，也随着离湿地距离增加，地下水位下降，植被盖度渐低，土地沙漠化发展，生态逆向发展。具体表现为，在水域、沼泽区多为茂密的芦苇、香蒲优势植被群落，植被盖度大；地表常年湿润区湿地随土壤盐分变化为芦苇群落、盐爪爪群落、柽柳 + 草本群落、草本 + 灌木群落等；在地表干燥、地下水位1~3 m区柽柳、白刺沙包发育良好，沙旱生、盐生植被发育；地下水位3 m以下，芦苇、盐爪爪等退化，6~7 m以下柽柳、沙蒿、中麻黄等大部分植被退化，白刺、红砂、膜果麻黄、霸王、泡泡刺等优势度增加，并随着地下水位继续下降，植被雨养性增强，植被覆盖度进一步下降，风沙活动加强，沙漠化重度发展 [ 30 ] [ 31 ] [ 32 ]。 近年来，随着石羊河流域综合治理工程实施，全流域开展了水资源统筹，关井压田、退耕还林(草)、节水灌溉、生态用水供给等措施，民勤盆地来水量显著增加，在局部地下水位上升和人工补给下，当前民勤盆地形成了4块较大的湿地，并在低地的低洼处形成了一些小面积湿地，在荒漠区形成了一些以湿地为中心的微域性生态屏障，对当地气候调节、改善生态、阻止沙漠扩展与沙漠化防治起到了重要作用。
system: [SYSTEM]

user: 摘 要 本文结合青衣江流域雅安段2019年7个断面水质监测数据，运用基于SPSS20的主成分分析法和Origin2018实现了青衣江雅安段水质时空分布特征的可视化分析，选取DO、COD Mn 、COD Cr 、BOD 5 、NH 3 -N、TN、TP为水质评价指标，运用综合水质标识指数法对水质进行分水期评价。结果表明：青衣江流域雅安段TN指标超标严重，全年多月属于地表水IV~V类；水质的时空分布差异显著，结合研究区水文气象和地理位置的特点，从时间分布看各断面丰水期水质优于平水期和枯水期，从空间分布看青衣江雅安段的支流濆江、陇西河和周公河监测断面水质最差，以连续性氮磷污染为主，其余断面水质能达到地表水III类标准。 关键词 :水质，时空分布，水质评价，主成分分析，水质标识指数法 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. 引言 河流是水进行迁移转化的主要途径，也是人类可利用淡水资源的主要载体 [ 1 ]。随着社会经济的快速发展和城镇化进程的不断推进，流域水环境污染问题变得越发突出，特别是在经济发展速度快和人口密度大的区域，河流水环境日益恶化已成为影响流域水安全和区域水资源可持续利用的主要因素 [ 2 ]。因此，对河流水质进行监测评价，研究其时空分布特征及演变规律，解析污染成因，能够为流域水环境污染治理和水资源的利用与保护提供动态理论依据。 雅安市地处长江经济带上游四川盆地西部边缘，是川西高原与成都平原的过渡带。青衣江流域贯穿整个雅安，地理位置在东经102˚16'至103˚43'、北纬29˚24'至30°56'之间，是大渡河左岸最大的一级支流，主源为宝兴河，流域区内面积10，053 km 2 ，境内流长184 km [ 3 ]，属亚热带湿润气候区，春季少雨干旱，盛夏暴雨洪涝，秋天阴雨连绵，冬季雨雪霜少。流域由于受特殊地理位置、地形作用的影响(南有东西走向的大相岭、峨眉山，北有邛崃山脉环绕，西有南北走向的夹金山，形成马蹄形地形)，构成了四川省著名的青衣江暴雨区。该区雨量非常充沛，流域多年平均降雨量达1776.7 mm。但在地区上变化较大，大致由西北向东南递增。流域内的气温与降雨在地区上的分布趋势较为一致，由西北向东南递增。域内多年平均气温大致介于14℃~18℃之间。最高气温出现在7、8月份，其多年平均气温约22℃~26℃，最低的1、2月份约4.6℃~8℃。极端最高气温36.6℃，极端最低气温−3.3℃。无霜期约300天。 雅安市随着国际化区域性生态城市建设和城镇化、工业化的不断推进，用水规模不断增加，对水资源质量的要求越来越高，水环境生态保护压力不断增大，近年来，青衣江雅安段支流水少、水浑、水脏问题逐渐明显。 本文以青衣江流域雅安段2019年7个监测断面的7项水质指标为研究对象，运用基于SPSS的主成分分析法(PCA)和综合水质标识指数法实现了研究区水质时空分布特征可视化分析以及水质综合评价，厘清了该区域水环境污染状况、主要污染因子时空演变规律，为该区制定合理有效的水资源开发利用规划和水环境污染防治措施提供了科学依据 [ 4 ] - [ 9 ]。
system: [SYSTEM]

user: 2. 数据与方法 2.1. 数据来源 本文水质监测数据由雅安市环境中心站提供，选取了DO、COD Mn 、COD Cr 、BOD 5 、NH 3 -N、TN、TP为水质评价指标。由于青衣江水期特征明显，分水期进行了水质评价，结合流域水文气象特征将水期划分为：丰水期为6、7、8、9月，平水期为3、4、5、10月，枯水期为1、2、11、12月。 青衣江流域雅安段共布设了7个监测断面，分别是青衣江雨城区上游的飞仙关镇(D 1 )、多营镇(D 2 )断面，雨城区青衣江3条支流汇入点的濆江(D 3 )、陇西河(D 4 )和周公河(D 5 )断面，水中坝(D 6 )和出境处的龟都府(D 7 )两个国控监测断面，青衣江流域图及雅安段监测断面布设如图1、图2所示。 2.2. 研究方法 2.2.1. 主成分分析 利用主成分分析(PCA)降低数据维度的原理，运用基于SPSS的PCA对流域水质监测数据进行主成分提取，根据水质数据的特征向量计算出主成分得分情况，结合Origin软件实现青衣江流域雅安段水质指标时空分布可视化分析。 主成分分析数学模型如下：设n维原始水质监测样本数据矩阵A 图1. 青衣江流域图 图2. 青衣江流域雅安段监测断面布置图 A = [ A 11 A 12 L A 1 n A 21 A 22 L A 2 n L L A n 1 A n 2 L A n n ] (1) 矩阵A经标准化处理后，即 A i j = ( X i j − X j ¯ ) S j (2) 变量间的相关矩阵为R，得到特征方程：Rb = λb。解特征方程将其特征根进行大小排列，λ 1 > λ 2 > L > λ n 。 ∑ j = 1 m λ j / ∑ j = 1 n λ j ≥ 0.85 (3) 确定一个m值，使信息的利用率达85%以上，选择前m个特征值对应的单位特征向量得到主成分计算式。 2.2.2. 综合水质标识指数法 单因子水质标识指数P i 由一位整数和小数点后两位有效数字组成，P i 值计算式由水环境质量标准确定 [ 10 ]。其中，整数代表了水质标识的水质类别，小数代表了监测数据在此类水质变化区间中所处的位置，单因子水质标识指数计算表达式为： P i = K i + ( C i − P i ， K i 上 ) ( P i ， K i 上 − P i ， K i 下 ) (4) 式中：K i 表示第i项水质指标所处的水质类别，取值为1，2，…，6；C i 表示第i项指标的实测质量浓度，mg/L；P i，Ki 上 表示第i项水质指标第K i 类水区间质量浓度的上限值，mg/L；P i ， Ki 下 表示第i项水质指标第K i 类水区间质量浓度的下限值，mg/L。 综合水质标识指数法是在单因子水质标识指数法的基础上对严重超标的水质指标进行加权处理，本文取权重为0.5，其表达式为： P = 1 2 ( n − 1 ) ∑ P i + P max ⋅ ( n − 2 ) 2 ( n − 1 ) (5) 式中：P max 表示n项单项指标水质标识指数中的最大值。水质标识指数评价标准如表1所示。 Table 1 Evaluation criteria for water quality labeling inde P值 1 ≤ P ≤ 2 2 < P ≤ 3 3 < P ≤ 4 4 < P ≤ 5 5 < P ≤ 6 P > 6 水质类别 I II III IV V 劣V 表1. 水质标识指数评价标准
system: [SYSTEM]

user: 4. 结论 依据2019年青衣江雅安段7个断面的水质监测数据，运用基于SPSS20和Origin2018对水质时空分布特征进行分析，以综合水质标识指数法为水质评价模型，分丰、平、枯水期进行水质综合评价。具体结论包括以下几个方面： 1) 2019年青衣江流域雅安段TN指标超标严重，TN指标浓度全年多月属于地表水IV~V类，雨城区支流的濆江、陇西河监测断面TN指标超标最为严重。 2) 青衣江流域雅安段水质时空差异性显著，干流水质总体较好，城区濆江、陇西河和周公河水质较差，丰水期水质优于平水期和枯水期。各监测断面水质时间变化趋势从优到差的顺序为：丰水期 > 平水期 > 枯水期；空间结构上，飞仙关镇、濆江和陇西河TN指标浓度较其他断面偏高，濆江和陇西河断面枯水期时段TP指标浓度值较高，其余指标变化不大。各监测断面水质从优到差的顺序为：多营镇 > 龟都府 > 水中坝 > 飞仙关镇 > 周公河 > 陇西河 > 濆江。 3) 2019年青衣江流域雅安段综合水质标识指数评价结果，丰水期水质优于平水期和枯水期，濆江、陇西河和周公河断面全年多月综合水质标识指数P值为4.01~5.04，水质类别为地表水IV~V类，TN指标P值为4.93~7.18之间，属劣V类，严重污染。飞仙关镇水质类别为地表水III类，多营镇为II类，水中坝和龟都府为II类。 4) 青衣江流域雅安段主要以连续性氮污染为主，濆江和陇西河断面枯水期磷污染较为严重，加之各支流上游水电站开发和季节性降雨影响，导致雨城区支流水少、水浑和水脏的现象。建议加快完善城区污水管网设施建设和提高污水处理效率，建立健全流域水环境污染防治法规。
system: [SYSTEM]

user: 摘要 目的：在新冠肺炎疫情下，探究个体防护病毒的公共危机能力对焦虑情绪的影响。方法：于2020年2月20日~2月29日运用焦虑自评量表(SAS)与个体公共危机能力问卷，采用网上发放问卷的方式对西南地区民众在疫情中的焦虑水平与公共危机能力进行匿名调查，内容包括性别、年龄和所属群体类别。结果：青年的焦虑情绪显著高于老年(t = 3.73， p < 0.05)；医护人员的焦虑情绪显著高于普通群体(t = 10.4， p < 0.01)和密切接触人群(t = 14.7， p < 0.05)。焦虑情绪与个体公共危机能力(−0.342 ** )显著负相关，表明个体公共危机能力得分越高，则个体的焦虑情绪就越低。进一步的回归分析发现，个体公共危机能力对焦虑情绪得分的回归系数Beta值达到了极其显著的水平(F = 47.405， P < 0.001)，并且预测系数解释了焦虑得分30%的变异，表明在新冠疫情下个体的公共危机能力是焦虑情绪的显著预测因素。结论：个体的公共危机能力对焦虑情绪具有缓冲作用，在疫情危机下青年人、医护人员的焦虑情绪水平显著高于其它群体，提高个体公共危机能力是改善疫情下民众焦虑的有效途径。
system: [SYSTEM]

user: 1. 引言 新冠肺炎疫情属于突发公共卫生事件，具有很强的突发性和不确定性，传染性强、范围广、传播速度快，使个体依靠理性和获取知识排除患病的能力大大降低了，不确定性、紧急性和危害性导致的安全感都降低了，大众对COVID-19的调查显示，98.54%的受访者出现过度恐惧、担忧和紧张(齐晔，陈刘欢，张栗等，2020)，而针对新冠肺炎患者的调查显示，76例病毒感染患者SAS评分为(54.39 ± 7.51)，存在轻度焦虑33例，中度12例，重度5例，由此可见感染患者中产生焦虑情绪的比例达到66%，对疾病本身不确定的焦虑感更为强烈(程丽，郑丽平，晏苏玉等，2020)。获取信息不确定性越多，信息辨识度越低、可信度差，越容易导致个体的焦虑水平增高，焦虑感也更为强烈。焦虑情绪是人们在日常生活中常常体验到的一种重要情绪。焦虑(Anxiety)是个人对即将到来的、可能会造成某种不良后果或模糊性威胁时所持有的警戒态度，同时伴随有一种紧张、不安、忧虑、烦恼、无助等不愉快的心理状态(黄希庭，2003)。通常也可以把每一种情绪都从3个方面进行定义：生理体验、心理感受、以及外在行为表现。在生理体验上，焦虑会导致肌肉紧张、出汗、呼吸急促和心悸，慢性焦虑则伴随着失眠、注意力涣散、耐心减退、兴奋易怒等问题；在主观心理感受上，高焦虑者会出现强烈的不确定感和不可控制感，对环境中的威胁相关信息投入更多的注意资源，同时倾向于将模糊信息解释为危险信息，在外在行为表现上，焦虑会促使个体回避潜在的威胁相关刺激，或者选择拖延而不采取能引发焦虑的行动(古若雷，徐鹏飞，徐蕊等，2016)。 焦虑不完全是一种消极情绪，适度的焦虑对动物和人类的生存有着重要的意义，但过度的、不当的焦虑则会严重干扰个人的社会生活。在重大疫情期间，适度的焦虑有助于引起市民对疫情的重视，积极采取措施防控疫情。但过度的焦虑会严重影响人们的生活，焦虑和恐慌可以瞬间在人群中蔓延，导致群体的盲目性和盲从性。处于焦虑状态下的人们往往会混淆主观感受和客观事实，更容易去相信谣言、阴谋论、夸大的数字和非黑即白的论断(邱稚博，2020)，也会过度地去解读一些不确定的信息。由此引起更多恐慌、焦虑，形成一个恶性循环，难以顺利度过疫情。 在面对重大危机时想要顺利地度过危机，需要具备良好的心态和聚焦问题的解决方式和行动方案。习近平总书记在统筹推进新冠肺炎疫情防控和经济社会发展工作部署会议上的重要讲话中，强调要总结经验、吸取教训，深入研究如何强化公共卫生法治保障、改革完善疾病预防控制体系、改革完善重大疫情防控救治体系、健全重大疾病医疗保险和救助制度、健全统一的应急物资保障体系等重大问题，抓紧补短板、堵漏洞、强弱项，提高应对突发重大公共卫生事件的能力和水平(本报评论员，2020)。可见危机应对和危机管理是有效控制危机情况、防止危机蔓延和扩大，并且最终战胜和攻克危机的必备素质。众多有关疫情危机的研究中，都提出要把个体应对危机的意识和能力纳入到公共危机治理中，成为公共危机管理的重要组成内容。 对“危机”的定义有很多，其中一般倾向于使用美国著名学者罗森塔尔的观点：危机是指“对一个社会系统的基本价值和行为准则构架产生严重威胁，并在时间压力和不确定性极高的情况下必须对其做出关键决策的事件。(杨晓燕，周懿瑾，2005)”国内学者对“公共危机”提出了自己的界定，李煜明提出，公共危机是突然发生并危及公众生命财产、社会秩序和公共安全，对社会价值和行为准则构成威胁，在时间压力和不确定性极高的情况下有待于政府和社会采取措施共同应对的公共性事件(李煜明，2005)。李燕凌、陈冬林和周长青认为：公共危机是指对社会公众具有巨大现实或潜在危险(危害或风险)的事件(苌瑞月，2007)。傅瀚霄认为公共危机即公共利益意义上的危机，是指一个突然发生并对正常的社会秩序、社会价值以至公共的生命财产构成威胁的事件(傅瀚霄，2007)。张伟指出公共危机是在社会运行过程中，由于自然灾害、社会运行机制失灵而引发的，可能危及公共安全和正常秩序的危机事件(张伟，2010)。可见公共危机具有公共性、突发性、不确定性、破坏性，危机侵害的对象可能是每个公民，公共危机能力不仅仅涉及到管理部门和智能部门，而且也是个体在危机情境下根据信息搜集、认知判断和自我整合而做出的应对方式，同时能够合理地选择和采取行动来应对危机带来的困难和障碍，是在公共卫生事件中特别针对公共危机而对个体所提出的要求和规则，也是公民基本素养的一种体现。基于上述，本研究将公共危机能力定义为，个体在发生危机的情景下，通过认知调整、行为适应和情绪调节等方式来适应危机，改善和降低危机的影响，从而有效处理危机的能力。 袁慧等人在研究中表明应对方式是影响个体心理健康的重要因素，在抗击疫情时医务人员产生的负性情绪反应与消极应对是呈显著的正相关，其产生的负面情绪越多，所采用的应对方式越消极，同时也指出消极的应对方式同样也可能对个体的心理健康状态造成不良的影响(袁慧，罗琳，吴建元，胡迪，雷凯，黄建英，2020)。有研究表明，在校大学生面对非典病毒的考验时，使用不成熟或消极的应对方式，可能会给心理健康带来损害，采取积极的应对方式可以缓解心理压力，维持心理平衡(高延，许明璋，杨玉凤，姚凯南，2004)。综上所述，在面对公共危机时，个体的应对方式、应对能力影响着个体的心理健康，情绪反应。因此本研究在次此新冠肺炎疫情下对公共危机能力对焦虑情绪的影响进行考察，同时提出假设：公共危机能力与焦虑情绪负相关，公共危机能力负向预测焦虑水平。
system: [SYSTEM]

user: 2. 对象与方法 2.1. 对象 2020年2月20日~2月29日期间，集中在云南、贵州、四川和重庆等地区选取18岁及以上民众作为研究对象，问卷回收后将答卷时间小于100 s的问卷删除，在问卷中有80%以上均选择同一个条目的视为无效问卷，删除无效问卷41份，最后有效问卷为359份，总体有效问卷率为89.75%。根据主要观测指标性别、年龄段(青年(18~30)、中年(31~50)、老年(51及以上))、群体类别(医护人员、普通人群、密切接触人群)，如表1所示。 Table 1 Distribution table of basic information of survey sample 性别 年龄段 群体类别 男 女 青年 中年 老年 医护人员 普通群众 密切接触人群 77 282 190 148 21 50 252 52 表1. 调查样本基本情况分布表 2.2. 调查方法 调查主要采用网上自填式匿名问卷的形式，基于微信用滚雪球抽样法邀请研究对象，被邀请者在线填写调查问卷，并请研究对象将问卷转发给其5名同伴。 2.3. 调查工具 2.3.1. 防疫公共危机能力问卷 采用自编防疫公共危机能力自评问卷，分4个维度，40个初始题目，5级量表评分。经过探索性因素分析和验证性因素分析，最终正式问卷共32题，分4个维度，这4个维度分别为：疫情熟知度、疫情适应、责任意识和自我约束。疫情熟知度是指个体对疫情相关的信息储备、危险识别、信息评估和风险管理等能力；环境适应是对现实环境变化进行有目的、灵活的心理和行为调节，减少不良情绪反应及恐慌、惊恐发作等急性应激反应；责任意识是主动和自觉地把防疫、控疫作为自己的一份职责，将自身力量与防疫事件联系起来，做好社会一份子应有的参与、组织、宣传和支持行为；自我约束是克服困难并加强自我监督、计划和执行防控措施的能力。该量表采用5点计分，1为非常不符合，2为比较符合，3为一般，4为比较符合，5为非常符合。其中1、2、6、14、17、21、23、25、26、24、27这11条为反向计分，分数越高，代表公共危机能力越强。该量表的内部一致性信度系数为0.807，各维度的内部一致性信度系数在0.720~0.759之间；量表的四因素假设模型满足拟合优度模型的条件，χ 2 /df的值小于5，RMSEA的值小于0.08，CFI及NNFI的值均在0.85以上，模型拟合较好，问卷的结构效度较好。 2.3.2. 焦虑自评量表 采用Zung于1971年编制的焦虑自评量表，采用4级评分，主要评定项目为所定义症状出现的频度，1表示没有或很少有时间有，2是小部分时间有，3代表相当多时间有，4是绝大部分或者全部时间都有。评定的时间范围是“现在或过去一周”。其中5、9、13、17、19这5个条目反向计分，20项总分均值为29.78 ± 0.46作为常模总分均值之上限。 2.4. 统计方法 应用SPSS25.0统计软件进行描述统计、方差分析、相关分析、线性回归。
system: [SYSTEM]

user: 4. 讨论 4.1. 不同性别、年龄段、所属群体对公共危机能力的影响分析 进一步分析不同性别、年龄段、所属群体这三个变量对公共危机能力的影响，对公共危机能力各维度及总分的得分情况及其差异进行比较分析，结果显示年龄段会显著影响个体的公共危机能力，而性别和所属群体对公共危机能力造成影响较小。 4.2. 不同性别、年龄段、所属群体对公共危机能力焦虑的影响分析 进一步分析不同性别、年龄段、所属群体这三个变量对焦虑的影响，对焦虑得分及其差异进行比较分析，结果显示年龄段、所属群体显著影响焦虑，说明不同的年龄阶段、所属群体在疫情下所感受的焦虑程度有很大的不同，而性别上焦虑水平不存在显著差异。 性别对焦虑影响不显著，这与其他有关疫情研究所表明的女性在疫情中比男性更容易感知到不良情绪，焦虑水平要高于男性的焦虑水平结果不符合，这可能是由于取样不均匀导致的。 年龄阶段是影响焦虑、公共危机能力的一个显著变量。对不同年龄段青年、中年和老年的公共危机能力的比较中，可以看到青年的危机能力是低于中年和老年人的，而在焦虑得分中则是青年人的焦虑高于中年和老年人。这与多数研究是一致的，例如，朱越等人对新冠肺炎疫情下负面情绪影响心理健康的条件过程模型的研究中显示：青年的焦虑水平高于其他年龄阶段。因此青年群体是值得关注的群体，其对于危机应对的能力不足，同时容易产生焦虑和抑郁的不良情绪，心理承受力较弱，所以社会、学校、家庭应该帮助青年群体提高公共危机能力，从而改善青年群体的焦虑情绪。 所属群体是影响焦虑的另一个显著变量，对不同群体医护人员、普通群众和密切接触人群的焦虑情况进行比较，可以看出医护人员的焦虑水平显著高于普通群体和密切接触人群。盛小燕等在COVID-19感染爆发期护理实习生的心理状况及睡眠质量调查中，发现护理实习生群体SAS的得分显著高于全国常模焦虑，表现出明显的焦虑；袁慧等人就疫情中医务人员心理状态进行了分析，结果也显示抗击疫情期间医务人员突发公共卫生事件心理状态各维度(包括焦虑)得分较普通人群(袁慧，罗琳，吴建元，胡迪，雷凯，黄建英，2020)。本研究结果与之一致。这是因为与新冠肺炎的战斗中，医护人员是主力军，尤其是一线医护人员战斗在最前线，与新冠肺炎密切接触，承担起高强度、高风险的工作，他们忙于照顾患者，无暇顾及自己的身体情况，有些医护人员感染上新冠肺炎，有些甚至离世，这些医护人员承受着巨大的心理和工作压力，因此我们需要关注医护人员，特别是一线医护人员的心理健康，给他们进行及时的进行心理疏导和干预，帮助他们提高公共危机能力，改善其焦虑情绪。 4.3. 公共危机能力与焦虑的相关关系及回归分析 4.3.1. 公共危机能力与焦虑的相关关系分析 对公共危机能力总分及四因素与焦虑总分之间进行相关分析，结果显示焦虑与公共危机能力的四因素疫情熟悉度、环境适应、责任意识、自我约束存在显著的相关性。 个体对疫情的熟知度与焦虑情绪呈负相关，这与袁静等人研究所指出的疑似患者的疾病认知程度与焦虑水平有关的结果是符合(袁静，杜云海，徐巍，陈亚梅，2020)。也与苏茂玲的研究所表明的对脑肿瘤患者及家属疾病认知干预，可以显著降低患者和家属焦虑的水平符合(苏茂玲，2020)。也有研究指出在非典期间获取相关信息是有利于降低民众的内心恐慌。因此，个体疫情熟知度提高，即个体了解更多与疫情有关的信息，对新冠肺炎病毒的认知程度提高，掌握更多自我防疫技能，可以降低个体自身的焦虑情绪。环境适应是指面对公共危机时个体对危机的风险感知和应对策略，从以往的研究中可以看到不同的应对方式会对个体的心理健康有着不同的影响，采取积极的应对方式有利于缓解负面情绪；危机意识的不足，应付突发性危机性事件的心理准备不足，会使个体在面对突发事件时容易产生心理紧张，使个人在认知和行为上的应对能力降低(朱蕴丽，苗元江，2005)。因此提高个体的环境适应能力，增强个体的危机意识，选择积极有效的应对方式可以改善个体在面对公共危机时的焦虑情绪。责任意识是指个体以一个公民身份来要求自己，主动和自觉地把防疫、控疫作为自己的一份职责，以理性的态度看待疫情危机，关注积极正性的抗疫力量，将自身力量与防疫事件联系起来，做好社会一份子应有的参与、组织、宣传和支持行为；《当代大学生公共危机防范和应对能力的透视与反思》一文，在阐述对大学生危机教育中，强调要突出社会责任意识的培养，特别是对公民意识教育，要认同公民身份，为他人和社会的整体利益服务，积极承担应尽的责任(周庆行，高定伦，韦丽，2006)。自我意识的增强，可以帮助个体在抗疫期间积极承担起社会责任，在以理性的态度看待疫情危机时，也能够引导和帮助身边的人正确科学的防控疫情，帮助他人和自己去有效调控情绪。自我约束是指个体有目的、有意识地执行疫情防护措施，克服困难并加强自我监督、计划和执行防控措施的能力。个体自我约束力提高，是疫情得以控制和管理工作得以顺利进行的必要保障。综上所述，要提高个体公共危机能力，缓解个体面对公共危机时的焦虑情绪，可以从信息辨识度、责任意识、环境适应和自我约束的几个维度入手。 4.3.2. 公共危机能力与焦虑的回归分析 根据公共危机能力与焦虑的回归分析显示，公共危机能力是焦虑的显著预测因素(F = 47.405)，公共危机能力越高则个体感知到的焦虑情绪就越低。 焦虑情绪与公共危机能力存在显著负相关，俞晓静等探讨的先前SARS疫情焦虑水平与应对方式的关系中，个体越是采用消极应对方式应对疫情，焦虑水平就越高，消极应对正向预测非典的焦虑水平(俞晓静，张巍，沈晓红，2005)。另一项研究也进一步证明了个体在疫情下的情绪反应取决于个体采取的态度和行为，积极应对方式负向预测了个体的焦虑、抑郁情绪(何宏灵，杨玉凤，2005)。综上所述，个体在疫情之下体会到的情绪反应与个体采取的应对方式和应对能力有关。但对此还未从公共危机能力的角度进行考查，公共危机能力是特别针对公共危机事件，有效应对危机的认知、态度和行为上的综合体现，是针对特定公共事件的应对方式和应对能力的综合体现。研究结果与假设相符。通过分析可知，焦虑与公共危机能力存在显著的相关性，且公共危机能力显著负向预测个体的焦虑水平，公共危机能力有助于提高个体的主动性和适应性，让每个公民更好的去应对疫情所带来的各种变化，从而缓冲焦虑情绪反应；低危机能力的个体在应对一个突发危机事件时，其在主动学习防疫知识、疫情适应、约束自我行为方面都存在较大不足，难以忍受长期的不确定性，对于由于疫情发展而产生的各种生活变化无法适应，更容易产生焦虑情绪。
system: [SYSTEM]

user: 5. 结论 1) 不同年龄阶段和所属群体，在疫情下其焦虑水平具有显著差异，青年人和医护人员的焦虑情绪值得特别关注，开展相应的心理教育和心理援助； 2) 公共危机能力与个体焦虑情绪显著相关，公共危机能力越强的个体，体验到的焦虑情绪可能越低下； 3) 公共危机能力负向预测焦虑，有助于缓冲焦虑情绪反应，因此提高公共危机能力是改善疫情下民众焦虑的有效途径，应从信息辨识度、责任意识、环境适应和自我约束的几个维度开展相关宣传教育。
system: [SYSTEM]

user: 摘 要 黄河流域水沙变化剧烈，为更好地利用和开发水资源，定量的分析出引起河川径流变化的原因。采用Mann-Kendall、Pettitt突变点检验和累积距平曲线法，对1971~2016年窟野河流域内气象数据进行处理，基于budyko水热平衡理论对窟野河径流变化归因识别，得出气象因素和认为活动因素对径流变化的贡献率。结果表明：降雨和潜在蒸发均为非显著性增加，径流显著性下降，人为因素的下垫面变化对径流减少的贡献率达到98%是最主要的影响因素。人为活动中，通过土地利用类型图的变化，看出大量的煤矿开采，退耕还林还草，以及城乡居民用水导致了径流的急剧减少。 关键词 :窟野河流域，Budyko理论，水热平衡，归因分析，土地利用，气候变化 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. 概述 由于气候和认为的因素导致近50年来黄河流域径流发生着剧烈的变化，引起人们的广泛关注，黄河流的气候条件却没有出现显著性的变化，人类活动就起到了很关键的影响因素。窟野河流域有大量的煤矿和村庄，所以探究各种因素对流域径流变化的贡献率就非常重要，可以有效地对水资源进行合理可持续的利用。 流域径流的变化在年尺度上受到通过气候变化和人类因素的影响，Bukyko假设构建了流域水量与能量耦合平衡方程 [ 1 ]。近年来很多研究人员基于这一假设提出了一系列反应下垫面参数的经验公式并验证了假设合理性 [ 2 ]。孙福宝等 [ 3 ] 通过对黄河中游多个子流域气候变化和下垫面变化进行分析验证了Budyko水热耦合平衡方程在黄土高原流域的应用合理性。郭巧玲等 [ 4 ] 认为窟野河流域1997~2010年，该阶段人类活动对径流减少的影响所占比例高达93.62%。刘二佳等 [ 5 ] 认为窟野河流域径流变化在1980年~2005人类活动占78.25%其中煤炭开发的影响大于水保措施导致的径流减少。白乐等 [ 6 ] 采用SCRAQ法估算认为在1997~2010年的变化期人类活动对降水的影响为87.99%。本文在Budyko水热耦合理论的基础上计算降水、潜在蒸散发人类活动的弹性系数，最终定量化的得出三者对径流变化的贡献率，为窟野河水资源利用提供理论支持。
system: [SYSTEM]

user: 2. 研究区概况与数据获取 2.1. 研究区概况 窟野河流域，黄河一级支流，坐标范围为108˚28'E~110˚45'E，38˚22'N~39˚50'N (地理位置见图1)，地处黄河中游支流部分，发源于内蒙古自治区巴定沟，经伊金霍洛旗和陕西省府谷县境内，于神木县注入黄河。两条主要支流乌兰木伦河和悖牛川河自神木境内汇合，由温家川站流出，全长242 km，流域面积达8706 km 2 ，窟野河流域因属于干旱半干旱的大陆性气候，径流量随季节降雨影响大，因此径流量具有明显的季节性变化的特点。降水时空分布极不均匀，从时间上看，汛期6~9月降水量占年降水量的75%~81%；7、8两个月的降水量占年降水量的50%~60%，暴雨为主；从空间上来看，降水量南多北少、东多西少。窟野河流域多年平降雨量约为399.13 mm (1971年~2016年)，多年平均潜在蒸散发量为1273.66 mm (1971年~2016年)。 图1. 窟野河流域地理概况图 2.2. 数据来源与处理 本文选取窟野河流域的把口水文站温家川站(见表1)，年径流量数据来源于控制面积最大的温家川水文站(1971~2016年)的日流量数据。本文考虑到站点空间分布以及气象站资料的时间序列和完整性选择了窟野河流域周边的神木、兴县、河曲、东胜4个国家基本气象站的1971~2016年逐日气象资料。 Table 1 Basic information of hydrological statio 流域 站点 经度 纬度 年均径流量/亿m 3 水文年份 窟野河 温家川 110.75 38.43 4.30 1971年~2017年 表1. 水文站基本情况