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

user: 3. 讨论 腹膜假性粘液瘤(Pseudomyxoma peritonei PMP)是一种以粘液外分泌性细胞在腹膜或网膜种植而导致腹腔内大量胶冻状粘液腹水特征的疾病[ 2 ] 。多数患者缺乏特异性临床表现，多表现为大量腹腔积液和全身消耗等症状，如腹痛、腹胀、腹部包块、恶心、呕吐、贫血、体重下降、疲乏无力、发热等。有的会呈肠梗阻，但以腹痛、腹胀、腹部包块最为常见。1884年Werth首次报道一例腹膜假性粘液瘤，患者是来源于阑尾囊腺瘤，1901年，FrankeI报道了一例来源于阑尾的囊肿。目前认为PMP虽可有不同的来源，但主要是来源阑尾粘液腺瘤[ 2 ] ，其次是卵巢，少部分可来源于胰、肠、胃、结肠等。阑尾粘液腺瘤是一种病理上分化很好的腺癌，破溃后可以种植播散，分布于全腹，一般无淋巴结转移，亦不出现远处脏器的转移，故称之为假性粘液腺瘤。 PMP发病率较低，由于以前人们对此病认识较少，故该病误诊率较高，且术前常难以确诊．大多数病例是在术中发现或剖腹探查确诊的。目前B超、CT、MRI对PMP的诊断有一定的价值，尤其对于肿块型及囊腺瘤具有较高的诊断价值，有助于发现病变及术前对肿瘤的定位，但无法做定性诊断。对于部分表现为糜烂或结节样的阑尾粘液腺癌上述检查阳性率较低。结肠镜检查对于来源于阑尾的腹膜假性黏液瘤具有较高的诊断价值。6例患者中通过结肠镜检查及活检做病理后5例得以明确诊断。1例因考虑为阑尾囊腺瘤而未取活检，后通过手术确诊。通过分析该6例患者镜下特点，主要表现为阑尾囊腺瘤、阑尾内口处糜烂或结节样增生以及肿块样表现。“果冻样”腹水为本病的特征之一，另外如果有大量腹水但不易抽出，可换用大号穿刺针或注射器在超声引导下反复用力抽取。如抽出果冻样腹水或仍不能抽出腹水则首先考虑本病，结合腹部CT等检查如除外来源于卵巢等部位病变后可行结肠镜检查，以尽快明确诊断。 另外一些研究报告指出，在腹膜假性粘液瘤病人中一部分患者CEA、CA125等肿瘤标志升高[ 3 ] [ 4 ] 。本报告中6例患者有4例CEA，CA125升高，故上述肿瘤标志物或许可以为诊断腹膜假性粘液瘤提供一定的线索。
system: [SYSTEM]

user: 摘 要 采用ANSYS程序建立反应堆吊篮在空气中的有限元模型，计算其在空气中的振动特性。重点研究压紧弹簧刚度对反应堆吊篮的振动特性影响。计算结果表明压紧弹簧的弯曲刚度和扭转刚度对吊篮结构的一阶梁式和轴向振动频率有明显的影响，刚度的存在使得结构的这几阶频率明显增大，对结构振型出现的先后顺序没有变化。在一定区域范围内随着压紧弹簧弯曲刚度的增加，吊篮的一阶梁式和轴向振动频率明显增加而吊篮的壳式振动频率几乎不受压紧弹簧弯曲刚度的影响。 关键词 :反应堆吊篮，压紧弹簧，振动特性
system: [SYSTEM]

user: 1. 引言 反应堆吊篮是装载和支承反应堆内部各种部件及堆芯组件的重要结构，研究表明，吊篮结构在正常工作时不仅会出现壳式振动，还会引发整体的梁式振动，这种振动有时会对反应堆安全运行造成严重的危害性 [ 1 ] 。对反应堆吊篮的振动特性研究具有重要的工程应用价值。目前已发表的工作中，常对结构作一定简化或分别采用吊篮梁式或者壳式来计算。文献 [ 2 ] 研究吊篮在空气中的振动频率，结合有限元建立了两类压紧弹簧边界简化模型，并建立吊篮的振动模型。文献 [ 3 ] 采用吊篮的梁式和壳式模型分别计算吊篮结构的梁式模态和壳式模态，并利用FRIZ理论 [ 4 ] 和ANSYS程序声流–结构耦合考虑了下降腔流体对吊篮的影响。文献 [ 5 ] 采用流固耦合方法对吊篮在空气和静水中的振动特性进行分析，计算结果表明吊篮在空气和静水中前几阶固有频率都包含梁式振动和壳式振动。文献 [ 6 ] 对压紧弹簧刚度进行了多种方法的计算研究，表明压紧弹簧分析接触问题时能采用二维模型计算。本文拟采用数值模拟的方法，合理简化压紧弹簧处的约束，研究压紧弹簧约束对吊篮筒体的振动特性的影响。进一步探讨压紧弹簧的工程应用价值，调整压紧弹簧参数，分析吊篮振动特性对压紧弹簧参数的敏感性。
system: [SYSTEM]

user: 2. 反应堆吊篮结构及模型简化 图1为吊篮结构的示意简图。吊篮通过吊篮法兰固定在压力容器的内壁凸台(如图2所示)，吊篮法兰的上部一般由预紧螺栓、上凸缘和压紧弹簧等组成，下部由压力容器的下凸缘支承。 实际建模时吊篮筒体采用Solid-shell单元，吊篮下支承板结构较厚采用SOLID185单元。在吊篮法兰处，下支承简化为间支，在上支承处由于压紧弹簧形变，给吊篮上支承提供弯曲刚度和扭转刚度(图3)。为准确模拟吊篮法兰处压紧弹簧的约束，有必要先计算出压紧弹簧的压缩刚度和扭转刚度。 本文计算中，压紧弹簧材料参数：泊松比 ，密度 ，杨氏模量 ；吊篮筒体材料参数：泊松比 ，密度 ，杨氏模量 。 基于ANSYS程序建立压紧弹簧有限元模型(图4)，根据预紧位移确定出吊篮法兰和压紧弹簧上下表面的接触区域，并在模型底面相应区域加固支约束，顶面对应区域施加垂直于顶面的面力F，进行静力学分析，得到节点的转角及位移(图5)。根据如下公式求得弯曲刚度值： (1) 图1. 吊篮结构示意图 图2. 法兰处压紧弹簧示意图 图3. 吊篮法兰约束示意图 图4. 压紧弹簧有限元模型图 图5. 压紧弹簧位移云图 (2) 式中， ——力； ——位移； ——压紧弹簧上节点的弯矩； ——节点转角。
system: [SYSTEM]

user: 3. 吊篮在空气中的振动特性计算 为考察压紧弹簧对吊篮振动特性的影响，考虑两种情况：(1) 考虑吊篮法兰上支承不受压紧弹簧约束；(2) 考虑吊篮法兰上支承受压紧弹簧约束。模型1为考虑法兰上支承无压紧弹簧的弯曲和扭转刚度；模型2为考虑吊篮法兰上支承受压紧弹簧的弯曲和扭转刚度。采用三维实体建模(模型及网格图如图6~8所示)，分别计算吊篮结构在空气中的固有振动频率以及梁式和壳式振型，计算结果见表1。 由表1梁式固有振动频率可见，模型2的一阶梁式振动频率和轴向振动频率明显比模型1的频率高；整体而言，模型2的各阶固有频率都较模型1高，除一阶梁式振动频率和轴向振动频率，其它的振动频率相差不大。由于模型1与模型2的区别就在于是否考虑压紧弹簧的约束，计算结果表明：压紧弹簧对吊篮筒体的一阶梁式振动和轴向振动有明显的影响；压紧弹簧刚度的存在使得结构的这几阶振型振动频率增大；压紧弹簧的刚度不会影响吊篮结构各阶振型出现的先后顺序(图9~11)。究其原因，相比较模型1， 图6. 吊篮结构实体模型 图7. 吊篮结构有限元模型(模型1) Table 1 Natural frequencies of core barrel vibrating in air (beam mode， shell mode 编号 振动频率(Hz) 振动属性 振动频率(Hz) 振动属性 模型1 模型2 模型1 模型2 轴向半波数m 环向半波数n 1 24.294 25.920 一阶梁式振动 69.651 70.557 1 2 2 81.559 81.562 环向扭转 59.737 60.066 1 3 3 87.724 92.798 轴向振动 81.107 81.176 1 4 4 130.132 130.26 二阶梁式振动 164.248 164.420 2 2 表1. 吊篮筒体在空气中的固有振动频率(梁式振型与壳式振型) 图8. 吊篮弹簧有限元模型(模型2) 图9. 一阶梁式振型图 图10. 轴向振动振型图 图11. 壳式振型图(m = 1， n = 2) 模型2仅在吊篮法兰处增加压紧弹簧的作用，体现在边界条件上，即增加了弯曲刚度和扭转刚度，此刚度对吊篮筒体结构的梁式频率影响较大，而整体的壳式刚度变化不大，因此，对壳式振动频率的影响相对较小。
system: [SYSTEM]

user: 4. 压紧弹簧刚度对吊篮振动特性的影响 由压紧弹簧的对称性分析，压紧弹簧的扭转刚度对结构整体振动影响较弱。为进一步考察吊篮振动特性对压紧弹簧参数的敏感性，本节通过改变压紧弹簧的弯曲刚度计算吊篮的振动频率，讨论吊篮振动频率随压紧弹簧弯曲刚度的变化规律。图12，图13给出了反应堆吊篮的固有振动频率随压紧弹簧弯曲刚度的变化图。 由图12和图13可见，当压紧弹簧弯曲刚度取值小于 以及大于 的区域，吊篮的振动频率变化很小；在 到 之间的区域，随着压紧弹簧弯曲刚度的增加，一阶梁式振动、轴向振动和壳式振动(m = 1， n = 2)振型的振动频率都有明显升高，其它振型振动频率变化不明显；整体而言，压紧弹簧的弯曲刚度的改变对吊篮壳式振型频率影响较小。研究结果表明：在本文的参数下，在压紧弹簧的弯曲刚度取值为 ~ 区间时，压紧弹簧的约束对吊篮筒体的固有振动特性有较明显的影响。当压紧弹簧弯曲刚度较小时，意味着相对整个吊篮筒体的结构以及约束而言，其刚度不足以影响吊篮筒体的动力学特性；而当压紧弹簧的弯曲刚度取值足够大时，其约束性质和固支的约束非常接近，这时如进一步增加其弯曲刚度，也不会对结构的动力学特性有很大的影响。 图12. 梁式振动随弹簧单元弯曲刚度变化图 图13. 壳式振动随弹簧单元弯曲刚度变化图
system: [SYSTEM]

user: 5. 结论 本文首先采用有限元方法获得压紧弹簧作用在吊篮法兰的弯曲刚度和扭转刚度，计算分析了在其约束下反应堆吊篮在空气中的振动特性。进一步调整了压紧弹簧的弯曲刚度，讨论了吊篮的振动频率随其弯曲刚度的参数变化敏感性。计算结果表明： (1) 压紧弹簧的刚度对吊篮筒体一阶梁式振型和轴向振动有较大的影响。 (2) 压紧弹簧的约束使吊篮结构整体的振动频率升高，但是不改变吊篮结构各阶振型的出现顺序。 (3) 在一定区域时，吊篮结构的一阶梁式振动、轴向振动和壳式振动(m = 1， n = 2)对压紧弹簧的刚度变化较为敏感，随着压紧弹簧弯曲刚度的增加，其振动频率增加。吊篮结果的壳式振型对压紧弹簧弯曲 刚度的改变不敏感。
system: [SYSTEM]

user: 沙漠、森林和河流的孤独
system: [SYSTEM]

user: 摘 要 《绿房子》是巴尔加斯略萨的代表作。小说通过其人物的设定阐释了深刻的文化内涵。本文着重分析了天主教文化在拉丁美洲的传播带有的殖民主义的强权性质，其“逻各斯中心主义”的特点使其与拉丁美洲传统文化和新兴资本主义文化产生尖锐的对立冲突。揭示出谋求主体间性，使不同文化在拉丁美洲这块土地上形成对话机制，是打破人与人、民族与民族、文化与文化之间冷漠隔膜，杂糅而不融合的孤独状态的可行途径。 关键词 :绿房子，天主教文化，资本主义文化，主体间性
system: [SYSTEM]

user: 1. 引言 小说《绿房子》是秘鲁著名小说家、诗人马里奥巴尔加斯略萨的代表作。《绿房子》发表于1965年，叙述从上个世纪20年代到60年代，时间横跨40年，展现了一幅广阔的秘鲁社会图景。故事是在沿海沙漠地带的皮乌拉城、森林小镇圣玛丽亚德涅瓦镇和马腊尼昂河各支流上发生的。小说以其结构的新颖著称。打乱的时序，时空的分割，细节的混溶，都显示出作者的一种有意为之的“混乱”。略萨也因其在小说形式技巧方面的创新被称为“结构主义大师”。对作品叙事时空机制、结构形态的分析历来是研究的热点。然而，《绿房子》并非仅是一部形式新颖的先锋之作，其破碎断裂的叙述、混乱杂糅的时空结构正是作品反映的文化内涵的一种外在形态，这种形式背后的“意味”也是值得探究的。 本文旨在通过分析《绿房子》中天主教文化影响下的人物命运，探究秘鲁乃至整个拉丁美洲天主教文化的一些消极特点，进而说明这些天主教文化的消极因素正是造成人与人之间冷漠隔膜、文化杂糅性强而融合性不够，并且阻碍拉丁美洲发展的文化根源。
system: [SYSTEM]

user: 2. 殖民者的传教所与殖民地的森林 自哥伦布发现美洲以后，来自西班牙和葡萄牙的老牌殖民者对拉丁美洲的征服就是武力征服和精神控制是齐头并进的。中世纪晚期，在新兴资产阶级的文艺复兴运动的大潮下，加上宗教改革对天主教会的冲击，天主教在拉丁美洲殖民地的大肆传播可以说是天主教会凭借其殖民者的强力，对精神文化影响控制权争夺上的一次全力反扑。正如何塞卡洛斯马里亚特吉所说：“对美洲的征服是最后一次‘十字军东征’，而西班牙的威力也随着征服者而衰落了。它的十字军的性质，决定了征服在本质上只是一项军事和宗教事业。它是由士兵和传教士共同完成的”[ 1 ] 。但是不得不承认，这次反扑他们成功了，而且可谓影响深远。至今，中外学界大多认为，拉丁美洲文化是西方天主教文化的一个分支，是天主教文化的一支“亚文化”。这种观点强调的就是，天主教因素在拉美文化中的压倒性优势地位[ 2 ] 。 天主教文化在拉丁美洲这块土地上是作为一种强大的殖民文化而存在的。然而，早在“西班牙人用剑和十字架统治了美洲”[ 3 ] 之前，在秘鲁这块土地上，古代印加人就已经创造了光辉灿烂的文明，并且这种文明也是将宗教和社会和政治制度结合在一起的。因此，可以说原著民所共同拥有着一些他们所独有的模式化的思维方式与价值观念是根深蒂固地埋藏在他们血液中的。这一点就不仅仅是靠着殖民者的铁蹄就足以踏平的，这也就是其强大的文化惯性。印第安原著民文化和天主教文化正是拉丁美洲“前殖民时代文化”和“殖民时代文化”这两端，联系着拉丁美洲几百年受殖民统治的现实和本民族遥远而深厚的历史。 小说中鲍妮法西娅是贯穿全书的中心人物，她的命运正是由于受到两种文化的影响而跌宕起伏，她在天主教文化和原著民文化的价值判断中左右摇摆。可以说鲍妮法西娅是文化碰撞的一个典型。她出生于森林地区的一个土著部落，后来被圣玛丽亚德聂瓦镇传教所的嬷嬷收养，在传教所受到天主教教育，并且学会了读书写字。但是当孤儿院又抓来两个土著女孩时候，她却在自己都不知道为什么的情况下私自放走了她们。 “‘也可能是魔鬼在诱惑我，嬷嬷’鲍妮法西娅说道，‘但我并没有察觉，我只是感到可怜她们，请您相信我。’”[ 4 ] 鲍妮法西娅说不出来放走奇凯斯土著女孩们的原因，只是说可怜她们，是魔鬼诱惑了自己，其实说明了来自土著部落的身世使她与土著女孩有种天然地联系，土著女孩的遭遇对她造成了心灵上的触动。这种触动是她身上埋藏的印第安文化的血脉的一种唤醒。尽管鲍妮法西娅是试图调和这两者的，她对传教所的嬷嬷充满了感恩，甚至遵从天主教习俗在安塞尔莫去世时去找加西亚神父为安塞尔莫做临终祈祷，但是她始终在二者之间徘徊。拒绝穿高跟鞋，在“绿房子”中化名塞尔瓦蒂卡(“丛林里的女人”)是她身上联系这印第安原著民文化这一端的再次显现。在鲍妮法西娅身上，两种文化同时并存并激烈碰撞，彼此都不能将另一方完全消灭。 鲍妮法西娅的命运悲剧也是由两种文化的对立引起。天主教文化是殖民宗主国的文化。“西班牙美洲的物质征服开始之后，精神的征服几乎立刻接踵而至。神甫们紧步士兵的后尘，又是还走在他们的前面，企图使沉沦于异教的印第安群众改变信仰”[ 5 ] 。小说中胡姆对殖民当局的抗议和失败就交代了天主教文化话语霸权的根基——征服者在政治经济领域的绝对优势地位。和殖民运动相联系，使得天主教文化在拉丁美洲的传播不是按照自觉的“接触–认同–接受”的步骤进行的，而是一种伴随着武力征服施行的文化强权，因此宗主国的殖民文化带有一种唯我独尊的话语霸权。鲍妮法西娅放走两个土著女孩的行为是背叛行为，不能被容忍，只能被逐出传教所。 “‘是你自己把自己赶走的，’住持说道，‘我们对你比对任何人都好，鲍妮法西娅，你本来一辈子都可以留在传教所里，可现在，等孩子们一回来，你就不能留在这里了。我也很遗憾，可是你的行为太坏了。我也知道安赫利卡嬷嬷舍不得你，但是为了传教所，你必须离开此地。’”[ 6 ] 与两个土著女孩被“抓”来一样，鲍妮法西娅来到传教所和离开都不是她自己的选择，而是由天主教“一元神论”所主导。鲍妮法西娅被驱逐辗转沦落是因为她背叛了天主教文化，在印第安文化面前发生动摇的结果。这也充分说明了天主教在其“一元神论”思想和占有殖民者绝对优势地位情况下面对印第安原著民文化显示出的自我中心主义和排他性。
system: [SYSTEM]

user: 3. 天主教的十字架与沙漠中的绿房子 长达350年的殖民历史已经在拉丁美洲的思想意识上形成了一种固定的结构形态，而资本主义以一种“新”的姿态对其产生冲击。小说中另一个主要人物是绿房子的建设者——安塞尔莫，他建造的“绿房子”和他传奇的一生可以看做是一场具有资本主义精神个人主义冒险，反映出了天主教文化对资本主义的排斥。 20世纪20年代的皮拉乌还是个外地人眼中孤独凄凉的城市，道路缺乏，沙漠使其全国各地隔绝，处在孤立状态，毫无商业气息可言。就在这个时候，安塞尔莫出现了。安塞尔莫似乎是从天而降，他的身世、绿房子的建立都被看做是奇迹，带有神秘主义的色彩。其实，安塞尔莫对绿房子的经营并没有什么特殊：建房子、招募乐师、雇妓女、烹饪美食、提供各类酒和饮料，灯火通明，夜夜笙歌，他的经营方式俨然就是资本主义的经营方式而已。但是，由于皮乌拉人闭塞的生活状态和心理状态，在他们看来这种资本主义的经营方式仍然是新鲜事物，是神秘的存在。 这是由于在殖民统治下时日已久的皮乌拉，固然其本身自然环境不利于发展，但更重要的是天主教文化的“正统”观念有维护传统制度的倾向，对这种离经叛道的“新”是持否定态度的。“天主教对人们从事世俗活动设置了种种规定，不能对无所不能的上帝的意志提出任何挑战，试图以满足上帝的要求来压抑人们许多生活需要的本能，在一种宗教炫耀中以精神的安慰来弥补物质上的匮乏。只有这样，人们才能获得拯救，死后才能进入基督徒在尘世梦寐以求的‘天堂’”[ 7 ] 。时代发展到20世纪，与伊比利亚封建贵族结合在一起的天主教文化轻视物质财富和劳动因素，已经不能成为拉丁美洲发展的促进因素。“从总的趋势来看，天主教尽管也在努力适应时代变迁所带来的挑战，但以它为核心的伦理观在根本上有敌视或蔑视新生事物的传统，至少在相当长的历史时期内是这样”[ 7 ] 。并且，正如前文提到的，西班牙、葡萄牙殖民者在拉丁美洲的天主教文化传播是一次封建统治的成功的反扑，从根本上就是与崇尚资本主义的新教划清界限的。因此可以说天主教文化是一种反资本主义的文化。 四个二流子的故事也是小说中的一条线索，他们代表了皮乌拉城人原本生活的真实状态，也是一般状态——懒散、颓废、终日闲荡、无所事事，最后甚至靠塞尔瓦蒂卡做妓女赚得的钱来养活。正如马里亚特吉所说：“作为人生观和精神戒律，西班牙的天主教却没有能力在其殖民地创造劳动因素和财富因素”[ 8 ] 。天主教文化中对物质财富追求的蔑视消磨了人对事业的进取心，在物质条件匮乏的情况下形成了一种畸形的空虚心理，表现在精神文化方面就是怠惰和拖延。有人甚至把“懒惰”作为拉丁美洲“种族精神”的三个基本特征 1 之一。四个二流子没有合法的工作，生活就是混日子，似乎每日只在饮酒聊天。由于在天主教文化中贫穷被看做是在世应受的苦难，为拉美人安于现状提供了借口，加之又没有资本主义的重视时间价值的观念，排斥资本主义实干精神，形成了整个社会像患上瘟疫一般疲软、怠惰，发展十分缓慢。 加西亚神父称“绿房子”为地狱，称皮乌拉城为“蛾摩拉”。他完全站在天主教文化卫道士的立场上，只看到“绿房子”作为妓院在伦理上的恶，并不能注意到它给全城带来的商业气息。可是资本主义取代封建制度这一来势汹汹的大潮不是靠着天主教文化就能够阻挡的。所以，就算加西亚神父率众烧毁了“绿房子”，更多的妓院建起来了，沙漠中的“绿房子”必然享有存在于世的第二次机会。 从另一方面来看，更难能可贵的是，作者对这种资本主义力量并没有盲目的赞扬，而是清醒的看到这种现代化力量有其进步的一面——使皮乌拉城增添了活力，“城市的面貌发生了变化。那安静的内地街道上充满了外地来的人......”[ 9 ] ——当然它也有明显的弊端——“绿房子”毕竟是座妓院，是堕落的场所。安塞尔莫建造“绿房子”和他与安东尼娅的爱情其实是这种“新”的力量的两个方面，前者是对财富、事业的追求，后者是带有反叛色彩的对个人主义的爱情的追求。“绿房子”是引人堕落的妓院，是以资本主义的方式肆无忌惮地宣泄欲望的场所(表现了资本原始积累的罪恶)；而安塞尔莫对安东尼娅的诱拐却是人类爱欲的表现，由同情开始，是真挚美好的。而以加西亚神父为代表的天主教文化却以其禁欲主义的伦理将二者混同，甚至安东尼娅的死成了“绿房子”被烧毁的直接原因，可谓本末倒置，这是对天主教伦理教条伦理的莫大讽刺。在天主教文化占主导的社会，这种资本主义的“新”的力量毕竟是微弱的，这也就是安塞尔莫悲剧命运的原因。 沙漠中的绿房子的兴衰史可以看做是资本主义在天主教文化统治下对皮乌拉城的一次渗透，联系着殖民统治下的拉丁美洲孤独的现实与未来要走的资本主义发展道路。一端是被延宕的“现实”(后殖民现状)，一端是迟缓的“未来”(现代化进程)，天主教文化固有的保守性，和其对资本主义的排斥显然成了阻碍拉美社会进步的一种力量。
system: [SYSTEM]

user: 4. 亚马孙河上的孤独 小说的中伏屋是以一个殖民者的形象出现的。他的大部分故事是是在亚马孙河上漂流中通过与阿基里诺的对话中回忆出来的。他勾结官府，占岛为王，杀人越货，对印第安原著民族进行无耻的掠夺剥削和残酷的压迫。这些行为中可以看到当初殖民者初登南美大陆那些血腥历史的影子。城市已经渗透着资本主义的商业气息，殖民者的幽灵仍然在河流上游荡，拉丁美洲创伤的历史、闭塞的现在和飘渺的未来同时在《绿房子》中上演。作者巧妙的运用时空断裂、破碎的结构将情节杂糅并置，正是表现了秘鲁这块土地上“历史”、“现在”和“未来”杂糅并存的社会现实。而究其原因，天主教文化缺乏主体间性是造成拉丁美洲发展不平衡，发展阶段共时性并置，彼此隔离、不相融合的原因。 主体间性问题就是在“自我与他者”的关系之内逐步消减“我思”的优越性，承认他人存在的意义，寻找“自我与他者”之间理想的关系样态[ 10 ] 。也就是寻求各不同民族文化间的可交流性。350年的殖民统治，天主教文化毫无疑问已经成为拉美社会的主流意识形态，其他文化尽管犹如彩虹般多姿多彩，仍然只能处在他者的边缘地位。作为文化主导的天主教文化缺乏主动与他种文化对话的自觉，只企图用暴力的手段同化和压抑它种文化，如小说开头警察和嬷嬷用暴力的手段从奇凯斯抓土著女孩到传教所。 “两个阿瓜鲁纳男人在步枪前面显得很驯服。悬崖近了，讨厌鬼：要不要把他们干掉？开枪吗？安赫利卡嬷嬷：野蛮，不许杀人。黄毛在岸边野兽抓起一个女孩的胳膊，另一手抓住另一个女孩的的脖子：警长，她们怎么不愿意下去？两个女孩企图挣脱黄毛，她们不喊不叫，只是用头和肩，用脚和退乱撞乱踢，左右扭动地挣扎。”[ 11 ] 可以看到天主教文化带有其文化殖民的情势，本身就是以“逻各斯中心主义”的对立态度来看待他种文化的。站在其自我中心的角度上，他种文化无疑是外围、是依附、是边缘化的。因此，在向印第安原著民传播天主教文化时，为了达到自以为正确的目的他们不惜使用武力手段。海德格尔把这种自我中心主义描述为：“人把自身建立为一切尺度的尺度，即人们据以测度和测量(计算)什么能被看做确定的——也即真实的或存在着的——东西的那一切尺度的尺度”[ 12 ] 。这种自我中心主义实质上是一种话语霸权，使天主教文化成为了排他的暴力主体。消解了在拉丁美洲这块土地上，与其固有的印第安传统文化和新兴资本主义文化沟通共识、互识互补以促进文化乃至社会发展进步的可能性，造成多种文化以对立的形式存在的局面。更为严重的是，天主教文化这种缺乏主体间性的传播割裂了其与印第安文化传统和资本主义经济发展的联系，本质上是反历史的。 值得欣慰的是，作者在小说的最后让加西亚神父说出了一句感人至深的“我来”，象征着不同文化和谐共存的一抹曙光。“从历史文化的角度来说，人类正是通过人与人之间、人与自然之间、不同文明之间的对话才能达到平等相处与和谐共存，从而使人类自身能够生活在一个没有危险的没有纷争的伊甸园”[ 13 ] 。天主教文化唯有探求主体间性，抛弃“我”而谋求“我们”的“共在”和“对话”，才是对其自身也是拉丁美洲发展有益的途径。 综上所述，天主教文化在拉丁美洲的传播带有殖民主义的强权性质，其“逻各斯中心主义”的特点使他与拉丁美洲印第安文化历史传统和新兴资本主义文化对立冲突，造成了《绿房子》中鲍妮法西娅和安塞尔莫悲剧的命运。扩大来看，这也是引起秘鲁乃至整个拉丁美洲社会不发达的一个重要原因。谋求主体间性，使不同文化在拉丁美洲这块土地上形成对话机制，是打破人与人、民族与民族、文化与文化之间冷漠隔膜，杂糅而不融合的孤独状态的可行途径。
system: [SYSTEM]

user: 摘 要 纸是由纤维素层状无序排布构成的多孔介质。近年来，检测试纸和纸基微流控芯片在疾病诊断、水污染和食品安全的快速检测中的应用受到人们广泛关注。纸内的渗流规律对于解释检测结果和提高检测的准确度和灵敏度都具有重要意义。本文概述了纸内渗流理论的研究进展，包括宏观毛细渗流及孔隙结构内微观渗流，并介绍了纸内流动的控制方法及其对检测试纸和纸基微流控芯片检测效果的影响。 关键词 :纸基微流控芯片，检测试纸，渗流，流动控制
system: [SYSTEM]

user: 1. 引言 纸是由纤维素和纤维素之间的孔隙所构成的多孔结构。纤维素通过机械或者化学制浆的方法从木料中分离出来，纤维素与水的悬浮物经过搅拌、过滤、丝网沉降、干燥等单元操作便形成了纸张。纸的孔隙结构和形貌具有各向异性和无序的特点(图1)。在沉降过程中纤维素逐层堆积，一层内的纤维一般不会交织到邻近两层以外的层内，不同层间纤维素通过氢键连接，从而维持结构的整体性[ 1 ] 。各层内纤维分布的位置、方向和形貌都是随机的。单根纤维素的截面中空(图2) [ 2 ] ，其截面直径为几十微米(~20 μm)，纤维长度为毫米量级(~2 mm)远大于截面尺寸。一般而言，纸在厚度方向具有10层左右纤维，在面内方向的特征孔隙尺寸为几十微米[ 1 ] 。 纸在许多应用中的性能都与其孔隙微结构内的流动特性相关。例如，办公用纸为保持表面光洁且不透光需要在表面喷涂微细颗粒进行修饰，纸的孔隙结构对喷涂料的渗流过程具有重要影响；墨料在纸孔隙内的扩散对印刷质量具有很大的影响；包装纸中气体的穿透性能也与孔结构密切相关。近年来，纸基微流控芯片被广泛应用于即时检测，具有快速、成本低和操作简单等优点，在疾病、环境和食品安全检测和监控等领域具有重大应用前景。这其中涉及到毛细力驱动下检测液体在试纸孔隙结构内的反应渗流过程，因此研究纸内渗流特性及其与微结构的关系对于解释检测结果和提高检测的灵敏度和准确度都有重要意义。本文主要概述了纸内毛细渗流理论和流动控制及其在即时检测应用中的研究进展。 图1. 纸表面微结构SEM图片[ 1 ] 图2. 单根纤维中空结构[ 2 ]
system: [SYSTEM]

user: 2. 宏观毛细渗流 2.1. 单向渗流 润湿液体在试纸内的自吸渗流由三相界面处的毛细力驱动，毛细力的大小由拉普拉斯方程确定： (1) 其中σ为液体的表面张力，θ为接触角，r eff 为多孔介质的等效半径。 饱和液体在试纸内的粘性流动阻力通过达西定律计算： (2) 其中，μ为液体的动力粘度；ε和K分别为孔隙率和渗透率，两者都是多孔介质的结构参数，v为流动速度，L为已吸入试纸内液体的长度。将流速写成吸入试纸内液体的长度对时间的倒数并代入式(2)可得： (3) 吸入液体的流动加速度仅在试纸与液体接触的初始阶段具有影响[ 3 ] ，如果忽略加速度引起的惯性力，吸入液体流动受毛细力与粘性力平衡控制： (4) 解 (5) 式(5)即为描述润湿液体在多孔介质中毛细渗流过程的Washburn方程，最早由Washburn考虑液体在毛细管内的上升过程推导得到[ 4 ] 。式(5)表明吸入液体长度与时间的0.5次方成正比。 Washburn方程适用于纸片水平放置下的吮吸过程，当纸竖直放置时还需要考虑重力的影响，考虑毛细力、粘性力、重力平衡的控制方程为 [ 5 ] ： (6) 液体沿纸面上升的过程中会从表面蒸发，蒸发导致纸内液体的流量(流速)增大，附加流量会产生附加粘性流动阻力，这部分附加阻力可采用达西定律描述，考虑毛细力、粘性力、重力、蒸发效应的控制方程为[ 5 ] ： (7) 式(7)中最后一项为蒸发导致液体在纸内流动产生的达西粘性流阻，P为纸的润湿周长，A为横截面积。 为湿纸面的蒸发速率——单位面积单位时间蒸发质量，假定为常量，即不随纸面位置和时间改变，其值取决于液体种类、温度、以及空气饱和度[ 5 ] [ 6 ] 。 Washburn方程及其改进模型假定液面所经过部分的孔隙内的液体完全饱和，孔隙内液体饱和度在吸入前缘位置存在跳跃变化，但实验观察发现在液体上升到一定高度后液体饱和度在吸入前缘附近存为梯度变化，通过求解Richard方程考虑非饱和液体在多孔介质中毛细吮吸过程可参考文献[ 7 ] 。 2.2. 径向渗流 径向渗流考虑在多孔介质中心半径为R 0 的区域不受限制供液，液体在毛细力驱动下在多孔介质内的径向扩散过程。为推导液体扩散半径与时间的关系，我们从达西定律出发： (8) 其中Q为体积流量，A (= 2πRH)为横截面积，H为纸的厚度。式(8)积分得到： (9) ΔP c 为毛细压差，R为湿部多孔介质的半径。将体积流量(流速)写成湿部半径对时间的倒数： (10) 注意，式(10)及式(2)中的速度v为多孔介质孔隙平均速度。将式(10)代入式(9)可得： (11) 式(11)积分并将毛细压差 代入湿部半径与时间的解析关系[ 8 ] [ 9 ] ： (12) 2.3. 液滴渗流 包括即时检测在内的很多应用中都涉及到液滴在纸内的渗流，目前有关液滴在多孔介质内毛细渗流的理论研究较少，主要是实验研究。实验结果表明液滴掉落在纸上后的渗流过程可分为两个阶段[ 10 ] ：初始侵入纸内过程和后期慢速扩散过程。液滴在纸面上扩散开的半径随时间关系为[ 10 ] ： (13) 其中V为液滴体积，K为经验参数，时间项指数n为0.3左右，比径向渗流情况下n = 1小得多。Danini和Marmur [ 11 ] 进行了如图3所示的四种情况纸吸水试验，(a) 水平单相渗流；(b) 液滴渗流；(c) 无限供液径向渗流；(d) 受限供液径向渗流。其中图3(d)中受限供液试验是将图3(c)中液池与纸之间连接的毛细导体在一段时间后移除，以类比图3(b)中液滴渗流时受限供液的情形。结果表明，液滴渗流的扩散速度比径向渗流慢，而与受限供液径向渗流的动力学接近，因此揭示了液滴渗流比径向渗流扩散速度慢的原因是由于受限和无限供液所导致的。
system: [SYSTEM]

user: 3. 微观孔隙结构渗流 宏观模型中将纸视作为均质材料，不能反映孔隙结构对流动的影响。纸的无序孔隙结构中孔隙率、渗透率、接触角发生局部变化，从而可能导致宏观模型所不能解释的现象。例如，宏观模型中假定吸入前缘是均匀推进，但这仅是在粘性力占主导时成立，当毛细力占主导时，吸入前缘呈指状突进，形成不均匀推进的前缘，如图4给出的纸巾浸入墨水所形成的锯齿状吸入前缘 [ 12 ] 。Horvath和Stanley [ 13 ] 的实验结果发现水面上升高度与时间的幂函数关系L ~ t0.38，这与Washburn方程预测结果L ~ t0.5不符合，当考虑动态渗透率后就能较好的解释该实验现象。宏观模型假定吸入前缘经过后的所有孔隙内液体完全饱和，但实际上部分孔隙的润湿需要一个过程甚至是不能被润湿。此外，宏观模型中一些与纸孔隙结构相关的参数，例如渗透率、毛细压差等也需要通过微观模拟或者实验的方法来确定。 纸孔隙结构内渗流模拟分为两个步骤，首先需要获得代表纸微结构的几何模型，然后在该几何模型的基础上求解流动控制方程。生成纸维结构的方法大体有根据纸在制造过程中纤维沉降动力学生成、随机生成法、三维重构建模、拓扑网络模型等。流动控制方程的求解方法可以分为采用有限元或有限容积法计算方法求解基于连续介质假说的控制方程如NS方程，以及格子玻尔兹曼介观模拟方法。下面对以上方法的一些代表性工作进行介绍。 Koponen等[ 14 ] 通过模拟纸制备过程中纤维在纸厚度方向的沉降过程生成了纸微结构模型(图5(a))，模型生成时考虑了纤维受力导致的弯曲变形。采用格子玻尔兹曼方法模拟了沿纸厚度方向流动的渗透率，模拟得到渗透率与孔隙率之间的关系(图5(b))。并基于数值模拟结果给出了预测渗透率与孔隙率和纤维半径之间的经验关系式。 Becker和Wiegmann等[ 15 ] [ 16 ] 采用随机方法生成了碳纤维纸的微结构(图6)，所生成的随机模型与实际纸的微结构类似：纤维在纸平面任意分布，在厚度方向层叠。纤维直径和孔隙率通过实验如SEM图 图3. 四种情况纸吸水试验：(a) 水平单向渗流；(b) 液滴渗流；(c) 无限供液径向渗流；(d) 受限供液径向渗流[ 11 ] 图4. 墨水吸入纸巾形成的锯齿状吸入前缘[ 12 ] 图5. (a) 基于纤维沉降动力学生成的纸微结构模型；(b) 模拟得到渗透率与孔隙率之间的关系[ 14 ] 图6. 随机法生成的纸的微结构：(a) Toray090纸；(b) SGL10BA纸[ 15 ] 片确定后输入到模型中。基于生成的随机模型，求解了纸内的单相流动，获得了纸在平面和厚度方向的渗透率；同时通过求解毛细两相流动，成功获得了毛细力–饱和度曲线以及两相界面分布。 三维重构建模一般是先利用CT或SEM对纸进行二维断层扫描，得到图形数据库，然后对二维断层扫描图样进行三维重建，获得固体骨架模型。目前空间精度可达到1 μm，模型还原度高。Mark等基于SEM图片得到的纸微结构统计信息生成了纸的三维微结构模型(图7)，利用VOF方法通过求解NS方程模拟了此三维微结构的两相流动，获得了纸的微结构流动参数(如相对渗透率)，然后将微观模拟得打的参数代入到宏观两相流模型中进行多尺度模拟[ 17 ] 。 在随机生成或者三维重构几何模型上采用有限元或者有限容积方法求解NS方程，需要大量的网格来捕捉不规则区域的流动和热边界层，对计算资源要求高且耗时。格子Boltzmann方法(LBM)方法是一种介于分子动力学与连续介质模拟之间介观模拟方法，处理复杂边界非常方便，而且可以自动追踪界面，在多孔介质流动模拟方面具有突出优势，因此许多文献采用该方法模拟了纸微结构内的单相流动和多相毛细渗流[ 9 ] [ 14 ] [ 18 ] [ 19 ] 。基于纸三维微结构重构模型，Hyvaluoma等采用LBM模拟了润湿液体在毛细渗流过程，模拟得到液滴在纸内渗流过程如图8所示[ 9 ] 。Wiklund和Uesaka采用LBM模拟了随机生成的二维纤维多孔介质内的毛细渗流过程，数值模拟结果展示了孔隙结构内的非饱和流动现象(图9)，即吸入前缘经过的区域有些孔隙没有被液体润湿，而是仍然由气相占据[ 18 ] 。 随机生成和三维重构方法都是尽量逼真多孔介质真实不规则的孔隙微结构。拓扑网络模型则是通过建立转化原则，将实际多孔介质中的孔隙部分简化为由规则管道和节点相连的网络构成，流体在管道和节点内流动，不参与流动的固体骨架部分则被排除到网络以为而忽略。为保证模拟结果与真实多孔结构一致，通过CT扫描技术等对实际多孔介质进行扫描，获得孔隙分布的信息，并据此确定拓扑网络模型中管道和节点形状、尺寸、连接特定等参数，使其流动阻力性质与真实孔隙内流动一致。拓扑网络模型虽然不能完全重现多孔介质的内局部流动分布，但是在很多情况下我们并不关心每个孔隙内的具体流动状态，而是多孔介质内的整体流场分布。拓扑网络模型可以给出多孔介质内的大致气液分布界面，压力分布，流量，渗透概率等等，能够满足大部分细观渗流问题的需要[ 20 ] 。 Bijeljic等[ 21 ] 采用拓扑网络模型模拟了重量场下润湿液体在堆积颗粒多孔介质内的毛细爬升动力过程。拓扑网络模型由球体和圆管按立方体方式排布而构成(图10)，其中球体和圆管半径都按实际多孔介质孔隙结构的统计信息给定。即在液体被吸入的初始阶段，毛细力大于重力，液体爬升速度很快，流动受粘性力主导，这时气液界面分明(即液体饱和度跳跃变化)；当液面爬升达到一定高度以后，毛细力与重力 图7. 基于断层扫描重构的三维微结构[ 17 ] 图8. 格子Boltzmann方法模拟得到的液滴在纸内渗透过程[ 9 ] 图9. 格子Boltzmann方法模拟得到液体在二维随机多孔介质中的吸入过程；图中灰色表示固体，黑色为液相，白色为气相[ 18 ] 图10. 多孔介质立方体网络拓扑模型；图中黑色和灰色分别表示由气相和液相占据[ 21 ] 相当，液面上升速度逐渐减慢，这时气液界面为多相流区域(即液体饱和度存在梯度变化)。拓扑网络模型成功预测了气液界面区域饱和度变化的现象。Pan等[ 22 ] 采用三维拓扑网络模型模拟了纸表面喷涂液体后的干燥过程，但该模型的建立并没有考虑纸的实际微结构。Ghassemzadeh和Sahimi首先通过SEM图片获得纸的孔隙分布，并在此基础之上建立了立方体拓扑网络模型描述纸的微结构，模拟了纸的单相[ 23 ] 和多相流动[ 24 ] 。
system: [SYSTEM]

user: 4. 流动控制及在即时检测中的应用 近年来世界各地艾滋病、登革热、心脏病、癌症等疾病的病例逐渐增加。比如，艾滋病导致每年大约两亿人死亡[ 25 ] ；食源性疾病导致每年大约48亿美国人染上疾病，128，000人入院，3000人死亡[ 26 ] ；全球每年大约有两百万人因使用污浊之水、环境和个人卫生原因死亡，并导致超过六十万霍乱病例发生[ 27 ] 。 病毒和水质的检测对于疾病诊断和降低水质污染对健康的危害具有重要意义。传统上一般采用台式检测仪器(例如聚合酶链反应分析仪和酶标仪)来检测及分析定量多种病毒、细菌或食品毒素。但这些分析方法非常昂贵，耗力和耗时，预处理的步骤繁琐[ 28 ] 。随着微流体技术的发展，低成本、小型化、便捷化检测手段日益受到重视。玻璃基微流体芯片需要复杂的制造工艺和昂贵的设备进行操作[ 29 ] 。相比而言，检测试纸(侧流及微流体试纸)设计简单、成本低廉、易于携带、操作简便，适用于发展中国家、农村、家庭的临床检验[ 30 ] ，作为体外快速检测装置，被广泛应用于医疗诊断[ 31 ] -[ 33 ] 、环境监测[ 34 ] 和食品安全及营养监测[ 35 ] 等领域。 市场上供应的侧流试纸(LFA) (例如艾滋病检测试纸)可检测血清或血浆里的抗原(Ag)或抗体(Ab)。一般试纸条采用双抗原夹心法免疫测定Ag，将制备的抗体(Ab)预先固定在纸上。当检测样本含有Ag时，Ag先与胶体金(Au)标记Ab结合，通过毛细作用前往硝酸纤维膜，与检测线上的Ab结合，在线上形成Ab-Ag-Ab-Au复合物，使检测线显为红色，而剩余的金标记Ab将与对照线上的二抗体(Ab’)结合，形成红色对照线(图11) [ 36 ] 。因此，检测结果显示两条红线判为阳性，一条红色对照线则判为阴性。检测结果根据累积到的量显现出检测线颜色的深浅，从而可以做出临床判断[ 37 ] 。除了检测Ab或Ag，LFA也可用于核酸检测，采用合成的DNA直击病毒或其它检测物的RNA或DNA结构，通过分子杂交或碱基互补配对原理，检测血液中是否含有病毒核酸，以诊断有无病毒感染(图11) [ 36 ] [ 38 ] -[ 40 ] 。 除了LFA，近年来，多项研究已采用纸(如色谱纸、滤纸、硝酸纤维素膜等)作为芯片主要制作材料，形成纸基微流控芯片(paper-based microfluidics)，简称纸芯片。纸芯片中一般使用疏水材料把亲水性的纸 图11. 侧流试纸：(a) 结构示意图[ 38 ] ；(b) 病毒抗原检测[ 14 ] ；(c) 病毒核酸检测[ 40 ] 材料分割成亲水和疏水相间的区域，形成纸通道。由于毛细力驱动，液体便可在纸通道内自发运动[ 41 ] 。通过控制纸芯片的流动，可以简化复杂的多步骤操作(如多个样品前处理步骤和有序化学反应)，实现样品或检测试剂的有效混合，减短检测时间及增强检测灵敏度[ 42 ] 。文献报道了采用物理、化学方法及阀门等各种流动控制方法，下面一一分别介绍。 物理流动控制方法包括石蜡打印式通道[ 32 ] 、中空通道[ 43 ] 及剃刀制作开放式通道[ 34 ] 等。一般而言，为改进检测的灵敏度，流体流速应该稍微降低，以允许样品和试剂在检测前达到充分的反应。例如，采用石蜡打印物理方法为LFA制造疏水区域，延缓流体流动。与一般的LFA相比之下，该模型可提高灵敏度高达三倍[ 32 ] 。此外，采用中空通道可提高七倍流速，减短了检测时间，已成功地在纸芯片上实现了牛血清白蛋白(BSA)与葡萄糖的并行快速检测。剃刀制作开放式通道也用同样原理实现了环境检测，如重金属离子(total Fe， Fe 3+ 与Ni 3+ )检测(图12) [ 34 ] 。 化学流动控制方法包括采用干蔗糖通道[ 44 ] ，可溶性聚合物桥梁[ 45 ] 等，这些方法大多是为了降低流体流速，简化多步骤操作，提高检测效率。例如，在纸上采用干蔗糖通道控制流体流速大小。将纸芯片设计成四通道式，实现疟疾检测的四大步骤：分别将样本-Au结合物、洗涤缓冲液I、信号增强试剂、洗涤缓冲液II，依照次序灌入四条含有不同浓度的蔗糖通道，以控制各通道的流体流动和反应时间，简化多步操作(图13) [ 44 ] 。 阀门装置被通常设置在流体通道中控制流动开闭，文献设计了包括按钮激活阀门[ 33 ] 、磁驱动阀门[ 46 ] 等各种方法。下面重点介绍按钮激活阀门工作原理。在激活按钮阀门前，亲水通道预先被小空间隔开，使用者可用圆珠笔按压并结合亲水通道，激活阀门，使流体流通。在生物检测应用方面(图14)，该装置可具有四个阀门，让使用者有选择性的检测尿液里的蛋白质、葡萄糖、酮或亚硝酸盐，适合用于低量样品下的检测[ 33 ] 。 图12. 物理流动控制方法：(a) 剃刀制作的开放式通道；(b) 环境重金属离子检测应用[ 34 ] 图13. 化学流动控制方法：(a) 干蔗糖通道；(b) 疟疾检测运用[ 44 ] 图14. 阀门装置：(a) 按钮激活阀门；(b) 尿液检验应用[ 33 ]
system: [SYSTEM]

user: 5. 结论 本文概述了纸内宏观和细观渗流的理论研究进展，以及检测试纸和纸基微流控芯片中纸内渗流的流动控制方法及其对检测效果的影响。纸内宏观毛细渗流主要采用Washburn方程及其改进模型描述，包括单相渗流、径向渗流、液滴渗流三种情形。基于纸微结构的细观渗流模拟可以考虑孔隙结构的不规则分布导致的变化渗透率和动态接触角效应，以及非饱和流动现象。细观模拟中生成纸维结构方法有纤维沉降动力学、随机生成法、三维重构建模、拓扑网络模型，流动方程求解方法包括连续介质的方法以及介观格子Boltzmann方法。纸内渗流的流动控制方法主要有物理和化学方法，结合在流道中设计阀门装置，可以实现对流动速度和反应时间的控制，从而提高检测速度和灵敏度。
system: [SYSTEM]

user: 摘 要 黔西县是贵州主要的产煤县之一，煤炭资源丰富，煤层主要产于二叠系上统龙潭组(P 3 1)地层中，属于织纳煤田。通过对主要矿区主采煤层系统取样分析，研究了不同煤层、不同矿区微量元素空间分布特征；与相邻的大方、赫章等区域煤层中微量元素含量对比发现，在陆源供给多的三角洲平原环境形成的煤层中微量元素含量相对高，在受海水影响较大的潮坪环境煤层中微量元素相对含量低，揭示了研究区陆源物质供给母岩——玄武岩具有较高的微量元素地球化学背景是其主要制约因素。个别煤层中有害元素Pb、Tl等含量高，是后期热液活动和断裂构造影响的结果。 关键词 :微量元素，沉积环境，主采煤层，二叠系龙潭组，贵州省黔西县
system: [SYSTEM]

user: 1. 引言 贵州素有“江南煤海”之誉，已探明煤炭资源储量约542亿t [ 1 ] ，列全国第五，是中国重要的煤炭生产基地。黔西县位于贵州西部的织纳煤田内，煤炭资源储量大，质量好，煤炭开发是县域经济发展的重要支柱，但煤层中元素地球化学分布规律研究程度低，成煤环境与微量元素的富集、贫化机理不清，是否能合理利用煤炭中稀有、有益元素无科学资料支撑，相关方面的研究严重滞后。在实施“西电东送”的煤转电工程中，大量的粉煤灰形成，为煤层中稀有元素综合利用研究创造了条件，同时，对整个贵州全省煤炭资源的合理利用、环境保护、煤炭中潜在的共伴生资源开发具有重要的理论及示范指导意义。
system: [SYSTEM]

user: 2. 黔西地区二叠系含煤地层特征 贵州主要含煤地层为上二叠统，形成时代为晚二叠世，自东向西因沉积环境不同分为吴家坪组(海相沉积)、龙潭组(海陆交互相沉积)和宣威组(陆相沉积)。黔西地区大地构造单元位于上扬子地台南部被动边缘褶冲带的黔中隆起区，含煤地层为海陆交互相和三角洲相沉积龙潭组(P 3 l )，由泥岩、炭质泥岩、粉砂岩、岩屑砂岩及煤层等组成 [ 2 ] ，厚度67.4~264 m，与上覆下三叠统飞仙关组呈整合接触，与下伏峨眉山玄武岩组呈假整合接触。根据其含煤特征、岩性组合，以标三顶为界龙潭组可划分上下两段。上段(P 2 l 2 )以细砂岩、粉砂岩为主，含煤6~19层，可采煤层较多(M2、M4、M5、M6、M7、M9)，且稳定性较好，煤层结构较单一，一般不夹矸，厚度37~131 m。下段(P 2 l 1 )以深灰色粉砂岩、泥岩为主夹煤线和煤层，含大量黄铁矿，由上而下逐渐增多。其含煤性较上段差，可采煤层数较少(M12、M13、M14、M15局部可采)，且稳定性较差，煤层结构一般较复杂，厚度31~128 m (图1)。
system: [SYSTEM]

user: 5. 结论 通过对贵州黔西地区主采煤层中微量元素分布的研究，结合威宁、赫章、织金、纳雍、大方、龙里、 图11. 贵州西部二叠纪龙潭期岩相古地理图 [ 19 ] 天柱等地煤层中微量元素分布情况，分析贵州二叠系煤层中微量元素分布的主要控制因素，并利用煤层中微量元素分布探讨成煤环境，获得以下认识： 1) 系统分析了黔西主采煤层的微量元素，认为煤层中微量元素含量主要受控于陆源碎屑供给情况，当时成煤环境陆源碎屑主要是峨眉山玄武岩，其具有高稀土和富集Cr、Ni、Co、Sb、Se等元素，直接影响煤层中微量元素的含量，陆源碎屑供给大的三角洲成煤环境，形成的煤层微量元素含量高，陆源碎屑供给少的潮坪成煤环境，形成的煤层微量元素含量低。 2) 贵州二叠系煤层具有较高的F含量，一般为193~341 μg/g，平均251 μg/g (N = 31)，其他微量元素均低。低温热液流体是影响黔西晚二叠世煤中微量元素局部异常富集和赋存状态的主控地质因素。煤中有害元素的富集主要受大断裂及其衍生的断裂所控制，多期次的低温热液流体成为有害元素的主要载体，通过断裂带进入煤层。因此，断裂带及附近的煤层由于后期热液改造，Tl、Pb、Sb、Zn含量很高，再利用断裂带附近煤层时应该进行处理，以减少污染。
system: [SYSTEM]

user: 摘 要 水电站建设及建成发电运行以后，伴随着水库蓄水，往往不同程度地触发岸坡岩土体失稳，不仅影响工程建设的实施，而且对拦河大坝和周围人民群众的生命财产安全造成严重威胁。本文以某水电站近坝库区的荒田滑坡为例，分析其滑动成因和机制；在启动速度和下滑速度计算的基础上，对其滑动入水后形成的涌浪进行分析和量化计算。根据山区河道型水库的特点提出弯道系数的概念并修正潘家铮计算公式，该系数主要由大坝附近涌浪监测结果反分析获得。基于数值计算确定不同水位情况下的失稳范围，采用上述修正方法对库区一处潜在失稳区岩土体进行滑动入水后的涌浪预测，从而获得较为可靠的预测结果，为水电站的安全运行提供重要依据。 关键词 :近坝库岸滑坡，水库蓄水，涌浪预测，弯道系数法，安全运行
system: [SYSTEM]

user: 1. 引言 水库库岸大型潜在滑坡的失稳会导致岩土体高速入水，不仅会堵塞河道、威胁航行船只及沿岸居民生命财产的安全，更严重的是还会激起巨大涌浪，造成坝顶过水甚至冲毁水工建构筑物。特别是正在蓄水或已经蓄水不久的库区，库水的变化往往对库岸潜在滑坡的稳定性造成巨大影响，一旦受到暴雨、库水位突然变化等多重外界影响因素的联合作用，一些处于潜在稳定状态或基本不稳定状态的库岸滑坡就可能失稳，如果失稳的滑坡体积较大，滑块速度较快，不可避免会激起涌浪。例如，1961年3月6日，正值施工期间的柘溪水库大坝上游右岸1.5 km处的塘岩光发生滑坡，体积超过140万m 3 的滑体以高达25 m/s的速度滑入深50余m的山区水库，激起的最大涌浪高度为21 m，至坝前达2.5 m，并漫过尚未完建的大坝顶部泄向下游，造成了巨大损失，当场死亡40余人[ 1 ] [ 2 ] 。1963年10月9日，意大利瓦依昂水库大坝上游峡谷区左岸山体突然滑下体积为2.4亿m 3 的巨型滑坡体，滑体以15~30 m/s的运动速度进入水深达235 m的水库，形成的涌浪高度达到175 m，翻坝的库水给下游数公里范围内的市镇造成巨大灾难，死亡人数就达到3000人 [ 3 ] 。1985年6月12日，三峡库区秭归县新滩镇发生体积超过2000万m 3 的特大型滑坡，滑坡毁掉大半个新滩古镇的同时，滑坡体前缘的土石堵塞江面达1/3，激起的涌浪波及上下游共40多km的江段，最高可达80多m，将上游2 km处的渔船打翻四艘，到达小镇上游秭归老县城时还有1 m高 [ 4 ] 。2003年7月13日，长江支流青干河边的千将坪村发生总方量约2400万m 3 的特大型滑坡，掀起大于20 m高的涌浪，打翻船舶22艘，造成14人死亡，10人失踪的重大事故 [ 5 ] [ 6 ] 。因此，对库区滑坡涌浪灾害预测研究是非常重要的。 影响水库库岸滑坡涌浪的因素是十分复杂的，而且许多因素不能明确确定而仅为估计值，另外，计算时边界条件和初始条件也较为复杂，还没有一种通用的计算方法。在计算水库库岸滑坡引起的涌浪时，主要使用试验计算方法和数值模型方法 [ 7 ] ，经验公式如E. Noda [ 8 ] 垂直下落模式和水平壁移动模式，无量纲组合模式，瑞士方法，潘家铮算法等 [ 9 ] [ 10 ] ；试验计算方法主要为具体滑坡的水工模型试验及多组试验后得到的统计公式 [ 11 ] ；数值模型方法如有限差分法、有限单元法以及耦合的欧拉与拉格朗日(CEL)算法、无网格SPH法等 [ 12 ] - [ 14 ] 。本文通过对荒田滑坡的现场调查分析其滑动破坏的成因机制，结合室内岩土体试验数据和实际监测结果，对其滑坡入水后的涌浪情况进行量化分析，提出适合该该滑坡涌浪计算方法。在此基础之上，对近坝附近一处重点潜在不稳定岸坡进行涌浪预测，讨论其对库岸周围和大坝的影响。
system: [SYSTEM]

user: 2. 荒田滑坡概况及成因机制分析 2.1. 基本概况 荒田滑坡位于水电站大坝上游约9.5 km的库区右岸，其剖面特征及外观形态全貌分别如图1和图2 图1. 荒田滑坡示意图 图2. 荒田滑坡全貌 所示。2009年7月20日凌晨3:36左右滑坡所在山体突发滑坡，形成一宽约500余m，高约350 m，总方量超过500万m 3 的大型滑坡，其中滑入江中山体约100万m 3 ，滑坡形成的巨大冲击波，掀起平均30~40 m高的巨浪。形成的滑坡后缘陡壁高约100 m，坡度约60˚，滑坡发生后在堆积体的后缘及中部形成了两个缓坡平台，其中后缘平台剖面方向宽约300 m，中部平台剖面方向宽约150 m。 2.2. 成因机制 据现场调查和已有的资料分析，荒田滑坡形成的原因主要有以下几点：1) 岸坡岩土体物质组成及结构。岸坡位于河床右侧凹岸出口部位，由于长期的凹岸水流冲刷作用，加之岸坡属坡向与岩层倾向小角度相交的顺向坡，为表层破碎岩体各类重力地质现象的产生提供了有利条件，从而在岸坡中下部形成较厚的第四系覆盖层，主要为崩坡积的块碎石土，广泛分布于区内台地及缓坡地带，为黄褐色夹块碎石，一般厚度10~20 m，个别在30 m左右。岩基为无量山群上统深灰色片岩、板岩、变质砂岩，在滑坡右侧前缘可见板岩，薄层状，岩层产状为N80˚W/NE∠50˚~60˚。区域地质资料揭示在岸坡附近发育有两条规模较大的断裂，其中八字耳朵断裂从岸坡前缘澜沧江河谷通过，另外一条老君山断裂则从岸坡下游一侧的垭口通过，由于断裂形成过程中构造挤压作用的影响，岸坡岩体结构完整性差；伴随后期的风化营力作用，在岸坡浅表层形成了较厚的结构破碎的强全风化岩体，根据滑坡后缘陡壁揭露情况，基岩强、全风化厚度可达20~30 m。2) 降雨和水库蓄水影响。根据降雨资料，滑坡发生前数天区内普将暴雨，造成的影响是：浅表层岩体饱水对岸坡加载的同时，降雨下渗造成岸坡全强风化岩体和第四系覆盖层的强度急剧降低。滑坡发生前水库已经开始蓄水，岸坡前缘坡脚处于蓄水位以下，由于水的软化弱化作用，前缘坡脚岩土体的强度被进一步降低，同时也导致其有效应力的减小。3) 岸坡前缘临空条件。从滑坡发生前的地形地貌来看，岸坡前缘临江高程约990~1100 m段地形较陡，坡度约50˚~55˚，这为滑坡的产生提供了良好的临空条件。 其滑动机制为：1) 岸坡表层较厚的第四系覆盖层和强全风化岩土体结构松散破碎，易于饱水，滑坡发生前的强降雨过程，使岸坡岩土体一直处于加载过程中，造成岸坡岩土体下滑力增大；此外，降雨下渗形成的软化弱化作用导致岸坡岩土体内部抗滑能力降低。2) 水库蓄水过程中，岸坡前缘坡脚岩土体处于水下并且岸坡地下水位线被抬升，一方面蓄水浸泡降低了岸坡岩土体的强度，另一方面地下水位线抬升产生的浮托力减小了岸坡前缘岩土体的有效应力导致抗滑能力降低。简言之，强降雨导致岸坡岩土体下滑力增大以及抗剪强度降低，蓄水过程导致岸坡前缘抗滑段抗滑能力降低，前缘良好的临空条件最终造成了滑坡的整体下滑。
system: [SYSTEM]

user: 3. 荒田滑坡启动和下滑速度计算 3.1. 滑坡启动速度计算 对于推移式滑坡而言，假定滑坡滑动前由主动滑动段(质量为 )和被动滑动段(质量为 )组成，一旦滑坡锁固段被剪断破坏，则滑坡主动滑动段内的岩土体颗粒结构变性能 会瞬间转化为主动滑动段内岩土体的动能，从而立即获得初始速度 ，滑坡中前部的被动滑动段内的岩土体就会被其推动一起运动，滑坡开始大规模滑动，可以说此时滑坡主动段和被动段内岩土体的速度一致，启动速度 被滑坡整体获得。计算公式如下： (1) (2) (3) 根据现场取样并室内试验的岩土体大剪试验结果，计算时参数选取如下，认为滑带土在运动情况下的内聚力 ，内摩擦角为 ，滑坡体物质的密度为 ，整个滑坡体积为 ，主动滑动段滑动的岩土体体积为被动滑动端岩土体体积的2倍，粘滞阻力系数取0.18，滑坡岩土体颗粒结构的应变为 ，岩土体颗粒体的弹性模量 。带入上述公式获得荒田滑坡的启动速度为 。 3.2. 滑体下滑速度计算 滑体下滑速度的预测和计算，有多种方法。本文主要侧重理论计算法和滑坡特征参数推算法。 1) 理论计算法 目前具代表性的有能量法、谢德格尔法和科内尔法。 ① 能量法根据能量守恒定理，导出滑块沿滑面下滑 距离(水平距 )后的滑速 为： (4) 式中： ——滑面倾角(沿滑动方向)； ——滑体单宽重量； 、 ——滑动时滑面抗剪强度参数； ——滑体质心落差； ——滑块与滑面接触面长(沿滑动方向)。 ② 谢德格尔法(A. E. Scheidegger， 1978) 考虑了滑坡体积效应。他根据33个滑坡的调查资料发现，滑坡的体积和等价摩擦系数 在对数坐标上呈直线关系。两者的关系的经验式为： (5) 式中： 为滑体体积； ； 。据 求出 后，按下式计算滑速： (6) ③ 科内尔法(H. J. Korner， 1976) 将大型滑坡视为流体，按流体力学方法计算滑速。设定阻力有两个部分，一部分由滑动摩擦形成，另一部分决定于滑动速度的平方，即： (7) 式中： 为法向作用力； ( 、 为流体密度和厚度)： 为流体紊流系数。经推算， (8) 式中： ； 为初速度。 这三种方法在应用上都有一定局限性，能量法基于已知滑面的纯理论假设，计算参数容易获取，但不能反应滑坡滑动的实际情况；谢德格尔法根据已有调查资料推导的经验公式，考虑了滑坡体积效应，在一定条件下能够接近真实结果；科内尔法对颗粒流滑坡更有适用性，但采用的参数难于确定，有待进一步研究。 2) 滑坡特征参数推算法 根据已发生的滑坡自身提供的表征运动特征的参数，可按下式(王兰生、董孝璧，1988)推算滑体最大滑速： (9) 式中： ——滑坡后缘落差，它近似等效于滑体质心的落差； ——滑面等效坡度，平面滑面的 ； ——滑面综合动摩擦角，以块状滑坡残体下伏滑面平均坡度表示。 对荒田滑坡采用不同方法进行计算的结果见图3。综合分析后确定其下滑速度为28 m/s，作为后续计算的参考。