分类:全 站  |  新闻  |  通知公告  |  重大计划  |  工作动态  |  队伍建设  |  关于我们  |  互动  |  图 片  |  视 频 查询到17条相关结果,系统用时0.022秒。

查询时间:时间不限

  • 时间范围


查询范围:标题+正文

排序方式:按相关度排序

检索范围:

    近日,2024年中国地球物理学会工程技术创新成果奖评审结束,由中国地质调查局成都地质调查中心(以下简称“成都地调中心”)牵头申报的“上扬子西缘盆山结合带区域地球物理调查关键技术及应用”成果荣获银奖。

    该成果依托成都地调中心承担的“四川会理-盐源地区地球物理调查”二级项目以及“重力场定量校正方法研究”国家自然科学基金,针对西南地区典型代表性的盆山结合带地形复杂条件下导致地球物理场畸变和各向异性问题,研究开发了盆山带重磁资料精细处理技术、盆山带大地电磁数据处理技术,建立了盆山带区域地球物理方法探测技术体系,解决了利用传统的地球物理处理技术难以“干净”地恢复地下介质的真实情况的问题,提高了对地球深部结构解释的精度。该技术在西南地区资源能源地球物理勘查中广泛应用。其中,有效支撑云南大关县木杆区块页岩气钻井定位,经过近年来的勘探、压裂,该区探获页岩气储量32亿立方,直接经济价值60亿元,实现了理论创新和应用效益,为乌蒙山区脱贫攻坚提供了技术支撑。

    下一步,成都地调中心将持续开展勘查技术方法创新,将理论创新和应用创新相结合,在“新一轮找矿突破战略行动中”做好科技支撑,发挥更大更有效的作用,保障国家资源能源安全。

    获奖单位

    应用证明

    成都地调中心牵头完成项目获2024年中国地球物理学会...

    近日,地调局地学文献中心(中国地质图书馆)举办了“地质专业知识服务系统推广及应用情况交流研讨会”,本次会议基于地质专业知识服务系统,发布了包括服务网站、微信订阅号、知识服务应用、院士团队及咨询课题的专题服务等一系列成果。

    中国地质调查局地学文献中心此次发布的“地质专业知识服务系统”是由地学文献中心承担的中国工程院项目,是作为中国工程科技知识中心地质专业分中心为公众提供服务的平台,共建单位还包括中国地质调查局实物地质资料中心、中国地质调查局实验测试中心、中国地质科学院岩溶地质研究所。目前该系统已收集整合了地质领域各类数据资源70余万条,涵盖文献、科学数据、统计数据、术语百科、专家机构名录、地质图件等,数据量达800G。

    当前已开发提供6个专题知识服务应用,其中术语服务被列为“地质云”第一批特色权威服务产品,“地质云”已将地质术语应用作为工具型服务进行集成。地质术语知识服务应用将地质领域权威词典资源进行提取,整合发布了7万余条地质术语。用户可对地质术语词条进行搜索,或按照学科分类进行浏览,了解术语名称及各不同来源解释。其他知识服务应用还包括地质专家学术圈,知识立方,地质图件GIS查询,地质灾害事实数据GIS图等,还专门开设了微信订阅号。

    地质分中心旨在构建地质专业领域的知识服务平台,聚合地质行业公开资源数据,建立资源知识关联,为地质工程科技用户提供知识搜索和知识服务应用,为院士团队提供咨询服务。在系统试运行期间,平台已经为众多科研人员提供了大量科研情报,获得了郑绵平,袁道先,赵文津和多吉等院士的好评和肯定。

    下一步,地质分中心将积极拓展资源收集范围,拓展资源共建单位,不断迭代升级服务系统,提升用户使用体验,根据专业需求开发新的知识服务应用,为院士团队及咨询课题提供更为有效的推送和专题知识服务。将以更大力度全面开展知识服务系统的运营和推广工作,未来有望为更多专家学者乃至全社会提供优秀的知识产品及服务。

    来自中国地质调查局局属单位、天津市国土资源和房屋管理局、天津市地质调查研究院、天津市地质环境监测总站、中国地震台网中心、国家基础地理信息中心等单位的专家对该系统的前期成果表示了肯定,并就知识数据挖掘与更新、系统服务对象定位、提升系统影响力等方面进行了研讨。

    地学文献中心副主任致开幕词

    项目负责人介绍项目进展情况

    与会代表观看成果展

    地质专业知识服务系统推广研讨会在津召开
    土壤是人类的衣食之源和生存之本,是最基本的生产要素。土壤的健康质量直接影响动植物和人类健康。为了更好地保护我们的生存之本,今天我们就来认识了解土壤环境问题中较突出的重金属污染。

    土壤重金属污染指的是什么 

    重金属通常是指密度大于5克/立方厘米的所有金属元素,包括汞、镉、砷、铅、铬、镍、铜、锌、钒、锰、锑等,其中前5种元素因其毒性大被称为“五毒元素”。

    顾名思义,土壤重金属污染就是重金属或其化合物通过各种途径进入土壤造成的污染。土壤遭受重金属污染的典型事例最早可追溯到19世纪发生在日本足尾铜矿山的公害事件,由于铜矿山废水排入农田,使土壤中铜含量高达200毫克/千克,不仅造成水稻严重减产,而且使矿山周围农田变为不毛之地。进入20世纪五六十年代,相继发生了举世瞩目的“八大公害事件”,其中发生在日本的“痛痛病”和“水俣病”公害事件就是土壤受到重金属镉和汞污染的两个典型。

    土壤重金属污染的危害 

    重金属可以污染水体、大气、土壤、作物等,但重金属不会像有机污染物那样被降解,因此通过食物链被生物体吸收后,会在体内积累,对人类健康造成巨大的威胁。有毒有害的重金属元素,例如砷、镉、铬、汞和铅,会对人体造成严重的危害,可能导致高血压、语言障碍、疲劳、睡眠障碍、提高攻击性倾向、注意力不集中、易怒、过敏反应、自身免疫疾病、血管闭塞以及记忆力下降等疾病和症状。重金属元素还会对人体细胞酶产生毒害作用。

    土壤重金属污染来源 

    土壤中重金属的来源可分为地质过程内源和人为活动外源两部分。地质过程内源又可分为继承型和次生富集型两类。继承型是指母岩中镉、汞、铅等有害重金属含量本底高,在后期的风化成土过程中,这些有害重金属继续保留在土壤中。资料显示,我国土壤大面积的重金属高异常主要是由成土母岩引起的,这些成土母岩多是富含铜、铅、锌、砷、镉等有害元素的硫化物矿床、黑色岩系、煤系地层等地质体,以及含锰、铬、镍的基性岩等。

    次生富集型是指成土母质中重金属元素含量并不高,但是在母岩风化成土过程中,化学性质活跃的元素,如钾、钠、钙、镁等易进入水体流失。而化学性质不活跃的元素,如汞、铅、砷等有害元素在原地的风化残留物中反而富集了。

    人为活动外源主要是指大量重金属通过人为活动进入到土壤环境中,其中主要是现代化工业,例如电镀、电池、化肥、矿业、造纸、杀虫剂、制革、塑料制品、冶金、采矿、化石燃料等制造、使用、活动过程中产生的含重金属的废水、废渣和废气。

    土壤中重金属的活性 

    土壤中重金属的含量和存在形态,很大程度上决定了其对环境、人体的风险高低。目前,土壤重金属的形态分级可分为离子态(水溶态)、可交换态、碳酸盐结合态、铁锰氧化物结合态、有机物结合态和残渣态。不同形态的重金属,其毒性、迁移性和生物有效性均有不同差异。

    一般来说,离子态的重金属移动性强,易被植物吸收,多存在于土壤溶液中或土壤黏粒表面,有着较高的生物有效性;而某些重金属离子会与土壤中的盐类(如磷酸盐、碳酸盐、硫化物、铁锰氧化物等)、有机质形成沉淀物、复合物或螯合物,移动性有所降低,生物有效性也随之下降;最为稳定的则是残渣态,一般存在于硅酸盐、原生和次生矿物等晶格中。

    影响土壤重金属形态分布的因素有很多,归纳起来可分为两大类:一类是土壤内因,即土壤理化性质,如pH值、土壤有机质、土壤质地、胶体含量、离子含量、Eh值、营养元素等;另一类是人类活动,如输入到城市土壤中的重金属的数量、种类的影响。

    土壤重金属污染修复的方法 

    根据修复方式以及处理后土壤位置是否改变,土壤重金属污染治理方法分为原位治理和异位治理。异位治理环境风险低,见效快,成本高,环境扰动大,如客土法、换土法、土壤淋洗法等。原位治理中主要包括物理修复、化学修复、生物修复以及农业生态修复。

    物理修复主要包括电动修复、电热修复等。电动修复主要通过电流的作用,土壤中重金属离子和无机离子以电渗透和电迁移方式向电极运输,然后进行集中收集处理。

    化学修复就是向土壤中投入改良剂,通过对重金属的吸附、氧化还原等作用,降低重金属的生物有效性,常用的土壤改良剂有石灰、沸石、碳酸钙、磷酸盐、有机质等。

    生物修复利用生物削减、净化土壤中重金属或降低重金属毒性。1983年美国科学家Chaney提出了利用超富集植物清除土壤中重金属污染的思想,即利用植物对土壤中的污染元素具有特殊的吸收富集能力,将植物收获并进行妥善处理后可将该重金属移出土壤,达到污染治理与生态修复的目的。

    农业生态修复主要包括两个方面:一是农艺修复措施,包括改变耕作制度,调整作物品种,种植不进入食物链的植物,选择能降低土壤重金属污染的化肥,或增施能够固定重金属的有机肥等措施,来降低土壤重金属污染;二是生态修复,通过调节诸如土壤水分、土壤养分、土壤pH值和土壤氧化还原状况及气温、湿度等生态因子,实现对污染物所处环境介质的调控。

    钨尾矿资源,待挖掘的宝藏 

    □郝小非

    钨矿是重要的矿产资源,已发现钨矿物和含钨矿物20余种,最具有开采价值的是白钨矿和黑钨矿,主要分布于江西和湖南。钨矿品位一般在0.1%~0.5%,对钨矿选别后,91%以上的固体废弃物被作为尾矿丢弃。我国每年约有1000多万吨钨尾矿被排放到尾矿库中,未被有效开发利用,不仅占用大量土地,还存在安全隐患。但尾矿不是没有任何利用价值的废物,在技术经济条件达到的情况下也是待挖掘的宝藏。

    随着人们对尾矿资源综合利用认识的提高,人们也越来越注意到尾矿也是一种资源,可以被开发利用。钨尾矿中一般含铜、钼、铋等重要金属矿物及石英、萤石、绿柱石和石榴子石等非金属矿物,随着选矿技术的提高,我们不但可以再次回收利用钨,也可以回收铜、钼、铋、铷、锂和钾等有价元素,还可以回收萤石、石英、石榴子石等非金属矿物。另外,我们还可以利用钨尾矿制备地聚物、微晶玻璃、陶瓷原料、免烧砖、透水砖等环保建筑材料。

    总之,只要我们结合钨尾矿资源特点,因地制宜地寻求钨尾矿资源化利用途径,一定可使钨尾矿变废为宝,实现经济效益、生态效益、社会效益的有效统一。

    “锆”诉 

    □雷晴宇

    对普通大众来说,单纯提到化学元素锆,可能相对比较陌生,但它却与我们的生活密切相关。比如,我们最常用到的卫生洁具、瓷砖等陶瓷产品,就是因为其中含有一定量的锆才让它的外表变得那么艳丽和富于光泽,也让它具有了极强的抗腐蚀和耐磨等性能。

    在工业实践中,由于提炼和加工困难、产量不多,锆被列为稀有金属。

    锆的应用领域非常广泛。其中,63%以上的锆以硅酸锆、氧化锆的形式应用于陶瓷、耐火材料领域,约有13%用在锆化学品领域,仅有3%~4%的锆矿石被加工成金属锆。因其具有惊人的抗腐蚀性能、极高的熔点、超高的硬度和强度等特性,锆还被广泛应用在航空航天等领域。

    锆的元素符号为Zr,锆元素在地壳中的含量仅约为0.025%。

    地壳中大部分锆呈分散状态存在于许多矿物中,已知含锆的独立矿物有38种,锆英石(ZrSiO4)和斜锆石(ZrO2)是主要的具有工业价值的含锆矿物。锆英石主要赋存于海滨砂矿中,是世界冶炼金属锆的主要来源。斜锆石主要产于碱性火成岩中,与霞石、霓石、磷灰石、萤石、钙钛矿、锆石、烧绿石等共生。

    含锆的天然硅酸盐ZrSiO4被称为锆石或风信子石,广泛分布于自然界中,具有从橙到红的各种美丽的颜色,自古以来被认为是宝石,印度洋中的岛国斯里兰卡盛产锆石。

    1789年,德国人M.H.Klaproth对锆石进行研究时发现,将它与氢氧化钠共熔,用盐酸溶解冷却物,在溶液中添加碳酸钾,沉淀、过滤并清洗沉淀物,再将沉淀物与硫酸共煮,然后滤去硅的氧化物,在滤液中检查钙、镁、铝的氧化物均未发现,在溶液中添加碳酸钾后出现沉淀,这个沉淀物不像氧化铝那样溶于碱液,也不像镁的氧化物那样和酸作用,Klaproth认为这个沉淀物和以前所知的氧化物都不一样,是由Zirkonerde(锆土,德文)构成的。不久,法国化学家de Morueau和Vauquelin两人都证实M.H.Klaproth的分析是正确的,该元素拉丁名为Zirconium,符号认为Zr,中国译成锆。

    1808年,英国的H.Davy利用电流分解锆的化合物,没有成功;1824年,瑞典的J.J.Berzelius用钾还原K2ZrF6时制得金属锆,但不够纯。直到1914年,荷兰一家金属白热电灯制造厂的两位研究人员Lely和Hambruger用无水四氯化锆和过量金属钠同盛入一空球中,利用电流加热500℃,取得了纯金属锆。

    锆矿资源是稀有金属矿产资源之一。世界各大洲均发现有锆资源,主要分布在大洋洲和非洲,美洲、亚洲、欧洲也有发现。锆矿床按其成因可分为脉型岩矿和砂矿两种类型,但由于岩矿结构形态复杂,分离共生矿物成本高及开采难度较大等原因,导致目前全球工业开采多以砂矿型为主 。砂矿主要包括滨海砂矿、冲积砂矿以及残积砂矿,其中滨海砂矿最具工业开采价值,规模和产量远大于冲积砂矿及残积砂矿。

    据美国地质调查局数据显示,2012年~2018年,全球锆资源储量维持在7500万吨左右,澳大利亚、南非、肯尼亚及莫桑比克4个国家锆矿储量合计6140 万吨,占全球的84.11%,矿床类型多以滨海砂矿为主; 印度、马达加斯加、巴西、中国、美国、乌克兰、印度尼西亚及俄罗斯等国家锆矿储量1160万吨,占全球的15.89%。

    我国的锆储量和美国基本一致,约有50万吨,仅占全球储量的0.68%。相较于其他国家,我国锆资源非常缺乏,所以,我国的锆资源主要靠进口获得。随着需求量不断增大,近年来进口以每年6%的速度增长,国内每年锆进口需求量达到90%以上。

    锆英砂主要用于生产化学锆、电熔锆、硅酸锆、金属锆等。2019年,我国锆英砂市场需求量为62.02万吨,而我国自有资源产量不足1万吨,近3年的进口量均超过100万吨。

    中国和欧洲是锆的主要消费市场,中国对锆的需求在全球占比高达52%。

    锆矿按照主要用途分为金属锆和工业锆两类。金属核级锆处于锆产业链最顶端,工业锆主要用在化工耐酸碱设备、电子行业等领域。中国是世界陶瓷工业生产和出口大国,硅酸锆则是陶瓷行业的直接和主要原料,陶瓷制品离不开装饰,好的装饰使制品身价百倍,装饰材料是装饰的物质基础,陶瓷色料是最重要的陶瓷装饰材料。由此可见,陶瓷色料在陶瓷装饰中的地位,也可知氧化锆在陶瓷装饰中的地位。同时,随着中国陶瓷产业的迅速发展,锆需求也随之猛增。

    矿山废水变废为宝的秘诀 

    □胡四春

    在矿山开采、矿物富集分离过程中,会产生大量的矿山废水,其中包括矿坑水、露采厂废水、选厂废水、尾矿库和废石场的淋滤水,这些废水不仅被白白浪费掉,而且还污染了地表水和地下水,危害环境。

    根据产生的途径不同,矿山废水性质相差很大。例如,矿坑废水pH值要么是强酸性,要么是碱性;选厂废水可能含有大量的重金属离子和有机药剂,这些都给废水处理及回用造成了巨大的麻烦。因此,根据废水产生的途径和废水处理后的性质进行分类处理和分质利用就成了把矿山废水变废为宝的关键。

    矿山采选废水常见处理方法 

    一般来说,矿山采选废水常见处理方法主要包括七方面:

    自然净化法。自然净化法作为最廉价、最简单的废水治理方法,被我国的选矿厂普遍采用。自然净化法常以尾矿库为构筑物,废水通过管道运输至尾矿库,在库内发生沉淀、水解、氧化、挥发、光照降解甚至生物分解等作用,使悬浮颗粒和残余药剂浓度降低,甚至基本去除。

    自然净化的效果与曝晒时间、光照强度、水体温度、初始pH值、溶解氧等因素有关。通常曝晒时间越长、光照强度越强、温度越高,自然净化效果越好。

    特点:自然净化法具有成本低、管理方便、无二次污染等特点,但存在净化不彻底、耗时长、气候等自然因素干扰大等问题,特别在高寒地区,往往会因为净化效率低下而影响废水的回用。因此,自然净化法通常可作为选矿废水的预处理方法,或用于成分相对简单的重、磁选废水的处理。

    酸碱中和法。酸碱中和法是一种传统的废水治理方法,因简单实用而被广泛采用。这其中既包括酸性废水中的H+(或碱性废水的OH-)与中和剂中的OH-(或H+)发生反应,生成中性水分子,同时矿浆的合适碱度也有利于重金属离子与氢氧根离子反应生成难溶的氢氧化物沉淀,从而消除重金属污染。

    生产实践中,常用的中和剂有石灰、消石灰、硫酸、碱性废水废渣(电石渣等)、酸性废水废气等。在选择中和剂时,应优先考虑厂区周边的废料,以达到“以废治废”的目的。理论上各重金属在一定pH范围内均能沉淀,因此控制好pH值是中和法的关键。

    特点:酸碱中和法具有管理方便、费用较低、操作简便、处理量大、适应性强和运行稳定等优点,但也存在一些问题,如在用石灰中和时,设备及管壁结垢严重、污泥增量较大、易产生二次污染等。

    混凝沉淀法。混凝沉淀法是目前治理选矿废水较成熟的一种方法,常与活性炭吸附或氧化法组成混凝沉淀——活性炭吸附法和混凝沉淀——氧化法。

    混凝沉淀法使用的药剂主要包括凝聚剂和絮凝剂两大类。凝聚剂主要有氯化铁、硫酸铁、硫酸铝、氯化铝、聚合氯化铝(PAC)、聚合氯化铁(PFC)、聚合硫酸铁(PFS)等,使用最普遍的絮凝剂是聚丙烯酰胺(PAM)。混凝剂的选择至关重要,它直接关系到净化效果的好坏。

    特点:混凝沉淀法可以有效去除废水中的悬浮颗粒和一些重金属离子,是一种成熟、稳定、高效的废水治理方法,但也存在对有机化学药剂净化不彻底,因药剂用量过大易产生二次污染等问题。

    化学氧化法。化学氧化法是深度治理废水中残留浮选药剂的有效方法,特别是近年发展起来的高级氧化技术(AOP)能彻底去除废水中持久性难降解有机污染物。

    化学氧化法的实质是,氧化剂通过夺取废水中有机污染物中的H原子等途径,将有机污染物氧化成无毒或低毒的小分子物质,或转化为容易从水中分离的物质,从而降低废水的COD、BOD。常见的氧化剂有臭氧、Fenton试剂、双氧水、次氯酸钠等。

    特点:化学氧化法治理废水具有操作稳定、反应彻底、处理效率高并能提高废水的可生化性等特点,特别对于处理高COD的有机废水具有显著优势,但也存在运行费用较高等问题。

    人工湿地法。人工湿地是仿照自然湿地人工修建并参与监督控制的具有流动或静止水体的浅水水域,是以基质-植物-微生物为核心的综合生态系统,可通过基质截留、过滤、吸附,植物吸收、拦截,微生物摄食、分解等途径去除废水中的污染物,充分发挥了物理、化学和生物的协同作用。

    特点:人工湿地法为治理废水提供了一条绿色化、生态化的技术路线,但也存在基质易堵塞、占地面积大、受气候等因素干扰大等局限性。

    微生物处理法。微生物处理法对于矿山酸性废水具有显著的优势。其净化原理是利用微生物的新陈代谢作用降解水体中的污染物,从而达到净化废水的目的。

    特点:微生物法治理废水拥有巨大的发展潜力,具有环境友好、选择性好、二次污染少等特点,甚至还可以回收某些重金属原料,但如何筛选出适应性强的菌种是个难题。

    矿山废水的分步处理和分质利用 

    一般来说,铅锌矿矿山废水的分步处理和分质利用分4步来进行:

    一是将铅精矿和锌精矿的浓密溢流水直接回用到各自的选别流程。

    二是向尾矿水加入一定量的钾明矾和阴离子PAM进行混凝反应和絮凝沉降,将尾矿废水中影响选矿指标的铜、铅、锌、镉、铬等重金属离子去除掉。然后,再加入一定量的椰壳型粉末活性炭,并通过纤维球过滤塔来去除掉影响选矿指标的部分有机残留药剂,适度处理后的废水大部分回用到选矿流程。

    三是采场废水一部分用于厂区绿化用水和尾矿干堆库区降尘喷淋用水。

    四是加入一定量的纯碱来降低水的硬度,通过砂滤和膜滤工艺来降低水的浊度等,然后将深度处理后的水分别用于陶瓷过滤机的槽洗水、酸洗水及浮选药剂的配制溶解用水。

    根据其性质及成分存在差异,选用合适的废水处理技术及回用方法,可提高选矿废水循环利用率,实现废水的清洁排放,真正做到变废为宝。

    生态画卷 资源综合利用有新突破

    资源、环境与生态问题已成为事关人类发展前景的全球性问题。近几十年来,随着人口急剧增长与经济快速发展,世界工业化、城市化进程不断加快,人类活动已成为全球变化的重要驱动力。在经济全球化、区域一体化不断深化的推动下,各国经济发展对相互之间资源、环境与生态的影响不断加大,人类进入了生态全球化时代。面对前所未有的重大而紧迫的全球性环境问题,世界各国在持续努力探索解决之道。党的十八大从新的历史起点出发,做出“大力推进生态文明建设”的战略决策;习近平总书记从新时代基本方略的高度提出要树立“两个共同体”理念——“人类命运共同体”理念与“山水林田湖草生命共同体”理念,为推进全球经济社会发展指明了方向,地质调查工作迎来了新的转型发展。地质调查工作如何适应与服务全球与国内生态文明建设并推动全球与区域问题的解决,亟待深入思考。

     

    11990~2015年不同国家矿产资源人均开采量与消费量变化

    地球系统问题的全球性与区域性

    20世纪50年代以来,人类活动对地球系统影响的程度和频度发生了急剧变化,人类施加于地球系统的各种压力进入“大加速”时期,地球从全新世跨入了新的地质年代——人类世。人类活动对地球系统的影响已经接近或超过自然因素引发的环境变化,并正在继续加剧,有可能产生不可逆转的后果。在第23届联合国气候大会上,来自世界各国的科学家发出警告:地球系统越来越抵近危险的“临界点”。

    1. 全球自然资源开发从线性增长转变为指数增长,发展中国家增长尤为突出

    过去的100多年,矿产、水、土地等自然资源开发经历了从线性增长到指数增长的转变。

    (1)矿产资源:全球开采总量快速增长,发达国家主导矿产消费,发展中国家开采快速增加

    1901年以来,全球矿产开采总量经历了缓慢增长、快速增长、稳定增长与急剧增长的变化。与1901年比较,2015年全球矿产开采总量增长了32.0倍,其中化石能源增长14.6倍,金属矿石增长41倍,非金属矿石增长49.3倍。根据开采量增长情况,矿产资源开发可划分为4个阶段:1945年以前,矿产开采量缓慢增长,年均增长0.59亿吨,人均开采量1.73吨;1946~1973年,矿产开采量快速增长,年均增长6.40亿吨,人均开采量增长到5.78吨,年均增长4.0%;1974~1997年,矿产开采增速减缓,年均增长6.15亿吨,人均开采量增至6.34吨,年均增长0.4%;1998~2015年,矿产开采量急剧增长,年均增长16.05亿吨,人均开采量增至9.01吨,年均增长2%。

    近几十年来,全球矿产开采与消费格局发生了重大变化。从开采来看,20世纪90年代中期之前,OECD国家主导全球,开采量占全球的41.8%,之后开采量占全球比例不断降低,到2015年降至23.0%,并且自2007年开始由增长转变为下降趋势;金砖国家开采量快速增长,在1995年超过OECD,占全球比例由1995年的37.9%升至2015年的51.6%。从消费来看,直到2007年,OECD国家消费量呈不断增长趋势,1990~2007年平均占全球总量的52.1%,2007年之后消费量降中趋稳,近年来稳定在295.42亿吨左右,占全球比例降至2015年的36.4%;金砖国家消费量在2000年之后快速增长,年均增长6.3%,在2010年超过OECD国家,到2015年增至360.57亿吨,占全球总量的44.0%;其余国家矿产消费量保持稳定增长趋势,年均增长3.1%。

    全球资源治理体系变革滞后于全球矿产开采消费格局的变化。1990~2015年,OECD国家人均矿产消费量大大高于其人均开采量,平均高出42.2%,且这一比例有增大的趋势。这表明,发达国家所开发的矿产根本满足不了其消费需求,通过进口越来越多的原矿石、矿产品与各种制成品来补充。金砖国家、其余国家人均开采量一直大于其消费量,说明发展中国家所开采的矿产在满足本国需求之外,有相当比例以原矿石、矿产品、各种制成品等形式出口。以金砖国家为例,2015年矿产开采量14.6吨/人,消费量11.7吨/人,在满足本国需求的同时,每人平均为其他国家贡献了2.9吨的矿产。目前的全球资源治理体系与发展中国家的贡献不相适应,亟需变革,以促进全球资源优化配置。

    (2)水资源:开采总量保持增长态势下呈现出显著的区域分化

    全球水资源开采在总量持续增长态势下呈现出显著的区域性差异。1901年~1950年,全球水资源开采量缓慢增长,由6713亿立方米增至12265亿立方米,年均增长1.3%;1951年~1980年,水资源开采量快速增长,年均增长3.2%;1981年以来,水资源开采量增速趋缓,年均增长0.8%。OECD国家水资源开采量在1980年由快速增长转变为稳定波动趋势,近年来稳定在9200亿立方米,占全球总量的23%。金砖国家水资源开采量自20世纪60年代以来保持快速增长的趋势,1960年~2000年年均增长2.4%以上,2000年以后增速有所减缓,到2015年增至17500亿立方米,占全球总量的43.7%。全球水资源开采量增长的主要原因是灌溉农业的快速发展与农业经济的持续增长。中国、印度等新兴经济体农业快速发展,加上持续的工业化和城市化,用水量有较大幅度的增长;欧盟、美国等发达经济体由于越来越多地进口工业制造产品与粮食,同时技术进步促使工业与城市用水下降,用水量自以前的增长转变为稳定或下降。

    地下水开采量快速增加,部分发展中国家含水层疏干问题严重。全球地下水开采量自20世纪60年代的3120亿立方米增至2010年的9820亿立方米,增长了3倍多。与水资源类似,地下水开采亦呈现出显著的区域差异。发达国家地下水开采在经历了一段时期的快速增长后已趋于稳定或缓慢下降。例如,美国地下水开采1950年~1980年保持了30年的增长,之后趋于稳定。发展中国家地下水开采自20世纪七八十年代以来处于快速增加的态势。例如,埃及1972年~2000年地下水开采量增长了6倍。地下水开采主要集中在亚洲国家,印度、中国、巴基斯坦、伊朗、孟加拉国等5个国家地下水开采量占全球总量的53.2%。地下水开采量的快速增加导致部分地区地下水位持续下降,引发了严重的生态环境问题,如泉水消失、湿地萎缩、地面沉降、海水入侵等。

    (3)土地资源:城市与农业用地持续扩展,生态空间不断萎缩

    1901年~2015年,全球土地利用变化的趋势是拓荒草原与森林来扩展农业用地,开发农业用地来扩展城市和基础设施建设用地,森林、草原、湿地等生态空间不断萎缩。农业用地面积扩展趋势趋于减缓。1901年~1955年,全球农业用地面积快速增长,年均增长0.88%,占全球土地面积的比例由20.6%增至33%;1955年~2015年,农业用地面积增速趋缓,年均增长0.23%,约占全球土地面积的38.0%。从区域上看,欧盟、东欧和北美的耕地面积有所下降,而南美、非洲和亚洲的耕地面积呈扩大态势。全球森林面积不断减少。1901年~1960年,森林面积平均以每年减少0.18%的速度逐年缩小,1960年以后森林面积缩小速度减缓,年均减少0.1%。

    城市化以前所未有的速度在扩张。遥感图像分析表明,全球城市面积6587.6万公顷,占全球土地面积的0.51%。城市用地占土地面积比例最高的地区是西欧(2.11%),其次是东亚(0.97%)、北美(0.72%)、东南亚(0.63%)。据统计,1950年~2015年人口大于1000万的城市群数量由2个增加到29个,人口500万~1000万的城市群数量由5个增加到45个。联合国粮农组织(FAO)估计,目前城市面积以每年200万公顷的速度扩展,80%的土地来自于农业用地。虽然城市占用土地面积比例很小,但是由于城市集聚了全球一半以上的人口,城市发展对生态环境的影响是巨大而深远的。

    2. 全球生态环境恶化趋势加剧,区域分化明显

    在不断加快的世界工业化、城市化进程作用下,气候变暖、自然灾害、水土污染等日益成为影响全球发展的重大生态环境问题。

    (1)二氧化碳等温室气体浓度不断攀升,全球气候变化加剧

    根据观测数据,大气中二氧化碳等温室气体浓度上升呈加剧趋势。1901年~1960年,大气二氧化碳浓度由296ppm增至316ppm,年均增长0.11%;1960年之后,增长速度逐渐加快,1961年~1997年均增长0.36%,1997年~2015年均增长0.55%,2015年大气二氧化碳浓度增至399.57ppm。大气二氧化碳浓度升高的主要原因是化石燃料燃烧和水泥生产排放了大量的二氧化碳。2015年化石燃料燃烧与水泥生产排放了360.2亿吨二氧化碳,是1990年的1.6倍。

    发展中国家开采了越来越多的化石能源,来满足发达国家的能源消费需求。在世界经济发展竞争加剧的背景下,很多发展中国家为了获得竞争优势,降低或放松了环境标准要求,推动高耗能、高污染、高碳产业发展;而发达国家对环境标准要求不断提高,以提高本国环境质量和生活舒适度。受此影响,高碳产业可能从环境标准高的发达国家向环境标准宽松的发展中国家转移,从而导致碳排放转移。全球碳计划(GCP)对1990年~2015年二氧化碳排放量估算表明:OECD国家因消费造成的碳排放大于其生产造成的碳排放,且差值越来越大;相反,金砖国家生产造成的碳排放大于其消费造成的碳排放,差值亦越来越大。这说明,发展中国家开发了本国越来越多的化石能源,加工、制造成各种产品出口到发达国家,承担了碳排放量上升与环境污染的代价。

    (2)重大突发性地质灾害呈上升趋势,经济损失快速增加

    全球重大地质灾害发生频次不断上升。联合国国际减灾战略机构EM-DAT灾害数据库收集了各国发生的重大自然灾害。入库灾害至少满足下列条件之一:造成10人以上死亡;100人以上受到灾害影响;政府宣布应对灾害紧急状态;政府在救灾过程中呼吁国际援助。1940年~2015年,全球发生重大崩塌、滑坡、泥石流地质灾害697次,造成6.5万人死亡,有记录的经济损失约89.4亿美元。上世纪40年代到80年代初重大地质灾害增长较慢,80年代以后发生频率快速增加,从80年代初的年均不足10次增加到近10年的年均18次。虽然发生频次增加,但是因灾死亡人数没有明显增长,单次地质灾害造成的死亡人数总体上是下降的,从1970年~1979年的136人/次下降到近5年的38人/次,说明各国地质灾害防治取得了一定成效。然而,地质灾害造成的经济损失自80年代以来快速增加,从70年代的平均每年0.14亿美元增加到近10年的平均每年1.76亿美元。

    不同国家地质灾害发生与防治情况存在显著差异。美国1960年~2009年地质灾害共造成336人死亡,直接经济损失12.4亿美元(按1960年折算)。1970年以后,美国地质灾害造成的死亡人数保持在很低的水平,平均年死亡人数在4人以下;1985年以前直接经济损失呈快速增加趋势,之后直接经济损失则呈减少的趋势。墨西哥1997年以前地质灾害发生在低水平波动,平均每年发生10次左右,平均每年导致近14人死亡;1998年以来,地质灾害显著增加,平均每年发生的地质灾害增加至86次,平均每年导致50人以上死亡。尼泊尔1971年~1992年发生地质灾害频次保持稳定,多在19次上下波动;1993年以后发生频次明显增加并呈周期性波动,平均每年发生120次以上,在高发年可达380次以上。

    (3)全球水土污染处于上升态势

    已有数据研究表明,全球水土污染呈上升趋势,随着部分工业企业(特别是高污染企业)由发达国家向新兴市场国家转移,新兴市场国家水体和土壤面临着越来越大的污染压力。

    地表水和地下水污染日趋严重。据联合国估计,全球每天大约有200万吨工农业和生活废弃物排入地表水体中,全球每年污水产生量高达1500立方千米。在发展中国家,80%的污水未经处理直接排放到河流、湖泊和海洋中。世界卫生组织统计显示,全球有8.84亿人缺乏安全饮用水,全球88%的腹泻与不安全饮用水、缺乏卫生条件有关,大部分分布在发展中国家。在快速城市化和农业种植区,地下水中的氮浓度不断上升,地下水质趋于恶化。在人类活动的作用下,孟加拉国、缅甸、阿富汗、柬埔寨、印度、中国等地区发生了地下水砷污染,影响了3500万~7500万人口的饮水安全。土壤污染问题在发达国家和发展中国家普遍存在。由于长达200年的工业化过程和现代工农业的发展,欧洲土壤污染严重。据欧盟调查,38个欧洲国家发现大约有250万个场地存在污染风险,其中有34.2万个已被确认为污染场地,需要进行修复。由于土壤污染的隐蔽性和复杂性,土壤污染问题在很多国家尚没有引起足够重视。

    地球系统问题解决的理论框架 

    不断加速的工业化、城市化与全球化耦合在一起对地球系统产生了前所未有的影响,促使人们必须从全球尺度去认识地球系统的变化机理;同时,不同区域或国家自然资源与生态环境变化出现了明显分化,与人类相互联系最为密切的近地表圈层资源、环境与生态问题呈现显著的区域性特征,促使人们必须从近地表圈层去认识地球系统的变化机理。在问题驱动下,随着全球观测、信息等技术进步,地球科学形成了一门新的分支——地球系统科学;在地球系统科学理论指导下,聚焦近地表圈层形成了一个新兴领域——地球关键带。

    近年来,我国从生态文明建设实践出发,提出了“构建人类命运共同体”和“山水林田湖草生命共同体”的理念。“人类命运共同体”的内涵是从生态、经济、政治、合作等方面构建全球治理体系,推动形成新型国际关系和国际新秩序;在生态方面强调生态环境问题无边界,保护地球系统是全人类的共同责任。“山水林田湖草生命共同体”的内涵是按照生态系统的整体性、系统性及其内在规律,统筹考虑自然生态各要素、山上山下、地上地下、陆地海洋以及流域上下游,进行整体保护、系统修复和综合治理。由此,学术界与政界在应对人类面临的地球系统问题方面高度契合,共同构成了完整的理论框架。

    1. 地球系统科学:服务构建人类命运共同体

    地球系统科学把地球看成一个由相互作用的岩石圈、水圈、大气圈、生物圈等圈层构成的统一系统,重点研究各组成部分之间的相互作用,了解整个地球系统的过去、现今及未来的行为,为全球生态环境问题的解决提供理论基础与对策方案。上世纪80年代以来,地球系统科学以全球气候变化研究为重点,技术方法不断发展,研究内容不断丰富,研究体系日趋完善与成熟。

     

    地球系统问题解决的理论框架

    (1)以观测、机理、建模与解决方案为重点,地球系统科学研究取得重大进展

    地球系统观测网不断扩展与升级,地球系统监测能力不断增强。美国NASA于1991年建立地球观测系统(EOS),利用卫星与其他手段对全球陆地表面、生物圈、地球空间、大气以及海洋进行长期观测;EOS之后,启动了地球系统任务(ESM),加深对气候系统与气候变化的认识;2017年,启动了下一代联合极轨卫星系统,用于天气预报和环境监测。美国地质调查局自1972年起陆续发射LandSat系列卫星,用于探测地球资源与环境,包括调查地下矿藏、海洋资源和地下水资源,监视农、林、畜牧业和水利资源利用,监测自然灾害和环境污染等。法国国家空间研究中心自1986年开始研发SPOT系列卫星,进行土地利用/覆盖变化、植被监测、自然灾害评估等。欧盟与欧洲航天局自2005年资助地球观测计划——全球环境与安全监测系统(GMES),由遥感卫星与陆地、海洋、大气等监测传感器组成,2013年更名为“哥白尼计划”,以扩大地球观测计划在公众中的影响力。

    地球系统变化与过程机理研究不断深化,揭示了地球系统要素不同时空尺度下的变化规律与影响。地球系统变化包括大气过程、海洋过程、陆地过程、冰冻圈过程等,这些过程相互影响、相互作用。由于碳循环是地球系统物质和能量循环的核心,全球碳循环及其对全球变化的响应研究一直是被广泛关注的前沿问题。人们对岩石圈、陆地生态系统、海洋、大气以及人类社会等碳库的储量、在全球碳循环中的地位及其作用机制有了深入的认识。人们认识到土地利用、覆盖变化是造成全球变化的重要原因,很多学者对土地利用变化引起的区域气候、土壤、水文、地质等因子变化及其对生态系统影响进行了大量研究。针对全球变化的生态系统影响,学者从植物群落、植物生理生态、地下生态、水生态系统、生物入侵、生物多样性等方面开展了深入研究。

    先后建立了多个地球系统模拟模型,地球系统变化预测能力大幅度提升。上世纪80年代以来,很多研究机构陆续开展了大气模式、海洋模式、陆面模式、海冰模式等地球系统模拟模型的研发和应用。2000年美国NASA提出构建地球系统建模框架ESMF,包括核心框架、天气及气候建模、数据同化应用等,为地球系统建模提供了一个标准的开放资源的软件平台。ESMF发展至今,已经拥有40多个模型,包含大气圈模型、大气动力学/物理学相关模型、海洋模型、陆地和陆表模型、水文学/分水岭模型等。欧洲提出了欧洲地球系统模拟网络(ENES)计划,包括地球系统模拟集成和气候资料存储与分发两个计划,目标是建立一个高效的欧洲地球系统模拟和气候预测系统进行集成模拟研究。日本在上世纪90年代启动了“地球模拟器”计划,于2002年研制成功,并在国际上率先开展了超高分辨率的全球气候系统模式的发展和模拟研究。中国科学院开发了地球系统模式CAS-ESM,集成了大气、陆面、陆冰、海洋、海冰等分量模式。

    应对全球变化提出了系列减缓、适应方案,服务制定政策、编制规划和措施决策。基于地球系统观测、机理研究与模型模拟预测,开展全球变化的适应与可持续发展研究是地球系统科学研究的重点之一。2015年,《联合国气候变化框架公约》近200个缔约方在巴黎气候变化大会上达成《巴黎协定》,将所有国家都纳入了呵护地球生态确保人类发展的命运共同体当中,目标是把全球平均气温较工业化前水平升高控制在2℃之内,并为把升温控制在1.5℃之内努力。越来越多的研究强调通过人类自身行为的改变,主动适应地球系统变化;通过土地系统和景观的重新设计,协调生态系统服务和人类福祉之间的相互关系;通过社会-经济-环境可持续性的综合协同,降低地球系统变化的风险。

    (2)促进自然科学与人文科学融合和推进更加平衡的多学科集成,成为地球系统科学发展的未来趋势

    国际科学理事会(ICSU)于2010年提出了面向全球可持续发展地球系统科学面临的5大挑战:一是如何提高对未来环境条件及其影响预测的实用性;二是如何发展、增强和集成必要的观测系统用以管理全球和区域环境变化;三是如何预见、识别、避免与管理破坏性全球环境变化;四是采取什么样的制度、经济和行为变化以迈入全球可持续发展路径;五是如何在技术研发、政策制定与社会响应中鼓励创新来实现全球可持续性。

    面临这些重大挑战,地球系统科学将会从自然科学主导的研究转变为有广泛的科学和人文领域参与的研究,从单学科主导的研究转为更加平衡的多学科集成研究。“未来地球计划”未来10年将集中在3个方面:动态行星地球——观测、解释、了解和预测地球、环境和社会系统趋势、驱动力和过程及其相互作用;全球发展——获得管理食物、水、能源、材料、生物多样性和其他生态系统功能和服务所需要的知识;可持续性转型——了解转型过程与选择,评估跨部门和跨尺度的全球环境治理与管理战略。

    中国所提出的构建人类命运共同体理念,得到了国际社会的高度认可。这一理念被联合国纳入相关决议,与“未来地球计划”等一起共同引导与推进全球生态文明建设。

    2. 地球关键带理论:服务构建山水林田湖草生命共同体

    地球关键带是指异质的近地表环境,岩石、土壤、水、空气和生物在其中发生着复杂的相互作用,在调控着自然生境的同时,决定着维持经济社会发展所需的资源供应。地球关键带科学为近地表圈层地球系统研究提供了一个整体框架,在此框架内开展全面、系统、持续、深入的跨学科研究。可以说,地球关键带科学是地球系统科学在近地表圈层的具体实现,为地球系统科学提供区域理论基础并服务于区域与全球可持续发展。

    (1)融合地质、水文、土壤、生态等学科,地球关键带科学快速发展

    通过探索,地球关键带科学形成了一条整合研究的技术框架:循环上升的调查-监测-研究体系。通过调查、监测和研究的循环进行,不断深化对关键带及其过程时空变化规律的认识;在此基础上,通过对图件、数据和成果集成分析,针对管理者、科学家、社会公众等服务对象生产各种产品,将关键带研究成果最大程度地传递给社会。

    调查是了解地球关键带组成与结构的基础,也是部署监测和开展建模的基础。2012年,美国地质调查局发布了其核心科学体系科学战略(2013~2023),明确将地球关键带作为其研究的核心靶区,提出针对关键带的结构和过程进行调查,建立关键带3D/4D地质框架模型。针对土壤侵蚀、盐渍化、有机质减少和滑坡等土壤环境问题,欧盟委员会发布了土壤保护主题战略,将传统的1~2m深的土壤层扩展到地表至基岩之间的未固结土层进行调查和研究。关键带调查的主要目标之一是回答“关键带如何形成与演化”这一基本科学问题。欧盟资助的欧洲流域土壤变化项目选择了代表土壤形成不同阶段的4个地区进行调查研究,分析确定关键带形成演化的影响因素和关键带生态服务的可持续性。

    监测是了解地球关键带随时间变化的基础,为建模提供所需的输入数据和校正数据。美国国家科学基金会于2007年启动了关键带观测计划,先后建立了10个关键带观测站,以流域为单元,对关键带各种要素进行长期观测。德国亥姆霍兹联合会于2008年启动了陆地环境观测建设项目,先后建成了4个陆地环境观测站,为区域尺度气候变化研究提供地下水、包气带水、地表水、生物和大气的基础观测数据。法国则通过提升现有的“河流盆地网络”所属的观测站,建设关键带观测设施,以流域为单元对关键带要素进行观测。欧盟委员会于2009年启动了“欧洲流域土壤变化”项目,选择4个典型地点建立了地球关键带观测站,将土壤监测作为长期观测的重点。

    建模对于深化对关键带形成、运行与演化的科学认识具有重要的作用,始终是关键带科学研究的重要领域之一。例如,美国关键带观测计划的重要目标之一是建立能够描述关键带生态过程、生物地球化学过程和水文过程的系统模型,定量预测气候变化、地质作用和人类活动下关键带结构和功能的响应。关键带过程模型大致可分为两类:一类是描述单个过程的数学模型,一类是描述多个过程叠加的耦合过程的数学模型。对于前者,目前已建立了较为成熟的模拟模型;而对于后者,是关键带建模的重点和难点,尽管近年来做了很多探索工作,耦合模型还远不成熟,仍在不断发展中。

    (2)随着地球关键带科学的形成与发展,或将促使地球表层研究发生科学变革

    地球关键带将与经济社会最密切的近地表环境作为独立的开放系统,为区域资源、环境和生态问题研究提供一个完整的系统框架。地球关键带科学研究尚处于探索阶段,近年的进展表明地球关键带科学有潜力促使地球表层研究发生科学变革,为经济社会面临的气候变化、生态系统管护、水资源安全、自然灾害防治等重大问题的解决展示了一种新的图景。未来地球关键带科学研究发展方向包括4个方面:开发一个统一的地球关键带演化理论框架;开发耦合的系统模型来探究地球关键带服务;开发一个集成的数据和测量框架并进行验证;建立多学科集成的地球关键带观测站。

    从国内生态文明建设的实践中,我国提出了“山水林田湖草是一个生命共同体”的理念。在内涵上,地球关键带与山水林田湖草异曲同工,前者侧重理论,后者侧重实践,目标均是推进区域生态环境治理。地球关键带科学是山水林田湖草系统治理的理论基础,后者则是前者与实践相结合的应用体现。地球关键带科学与山水林田湖草生命共同体理念共同构成了区域生态环境治理的理论框架,共同推进区域可持续发展。

    对地质调查工作的思考

    地球系统问题得到了政府与学术界的高度关注。在社会治理层面,围绕人类社会持续发展需求形成了“两个共同体”理念——人类命运共同体与山水林田湖草生命共同体。在学术层面,随着全球观测、信息等技术的进步,以问题为导向,地球科学形成了新的分支——地球系统科学,聚焦近地表圈层衍生了“地球关键带”新领域。由此,政府与学界在应对地球系统问题方面高度契合,共同构成了完整的理论框架。地质调查工作应树立人类命运共同体与山水林田湖草生命共同体理念,以地球系统科学理论为指导,以地球关键带为重点,加强调查、监测与机理研究,加强综合评价,服务和支撑生态文明建设。

    一是以地球关键带为重点加强综合调查评价。将地球关键带作为地质调查工作的重点靶区。按照统一的技术规范和标准,开展不同尺度的专业性基础性地质调查,充分反映地质框架的成土条件、成矿条件、水文条件等多种属性,建立地球表层三维地质框架模型。充分利用现代信息、网络、大数据等技术,加强区域问题综合评价,形成基础扎实、数据可靠、形式多样的综合评价产品,服务区域生态治理与自然资源综合管理。

    二是以服务生态保护修复为目标加强生态地质调查。根据自然资源管理与生态保护修复需要,选择典型地区探索开展生态地质调查,形成生态地质调查技术规范。根据自然资源勘查开发的源头保护、利用节约与破坏修复全过程需要,推进不同尺度生态地质调查,提出生态保护修复地质解决方案。

    三是以服务全球资源治理为重点加强全球问题合作研究。以“一带一路”倡议为抓手,加快推进矿产资源勘查开发国际合作,加强产能合作,促进全球资源优化配置。立足我国优势,在前沿与关键领域,策划实施地学大科学计划,以全球岩溶动力系统资源环境、地球化学调查、青藏高原特提斯演化与资源-环境效应等为重点,推进国际地学大计划合作。

    四是以资源环境要素为重点加强地球系统探测与监测。采用卫星遥感、航空遥感等对地观测技术,定期采集全球与区域资源环境要素数据。协调、整合、新建观测站点,形成地球关键带综合监测网。开展区域自然资源数量、质量与生态综合监测,及时提出预警。围绕深部资源勘查开发与灾害防治需要,加强地壳深部探测。

    五是以提升自然资源管理决策支撑能力为重点加强地质大数据建设。整合现有地质、资源、环境、生态等调查数据,构建地质大数据核心数据库体系。建立资源环境要素数据动态更新机制,实现地质大数据与自然资源管理需求在时空上的契合。与经济、管理、社会等相关基础数据无缝链接,为自然资源管理与资源环境治理提供全方位支撑。

     

    地球关键带研究的调查-监测-研究循环体系框架

    六是以过程机理研究为基础加强综合评价。基于三维地质框架模型,加强地球系统物理过程、化学过程、生物过程的机理研究,建立地球系统或地球关键带模拟模型。基于机理模型,考虑不同社会经济发展情景,对所面临的问题进行综合评价,有针对性地提出地质解决方案。

    (作者单位:自然资源部中国地质调查局发展研究中心)

    解决地球系统问题需要新站位

    8月6日凌晨,因连续强降雨,云南省富宁县阿用乡-那来村发生山体滑坡,阻断河流形成堰塞湖,蓄水约15万立方米,威胁下游20户91名村民生命安全。

    自然资源部地质灾害防治技术指导中心迅速组建5人专家组于10日17时,到达灾害现场。与当地自然资源主管部门组成联合专家组,利用无人机三维实景建模技术,对滑坡体方量、边界及成因进行分析,并制定初步防治方案。经调查,滑坡体积约20万立方。

    11日,专家组重点对滑坡后缘稳定性及可能影响范围进行调查研究,针对滑坡体处置及坝体加固等工作,建议及时填埋后缘及周界裂缝,并设专人进行实时监测。目前,堰塞湖已成功泄洪,山体滑坡趋于稳定,险情解除。

    现场调查

    灾害现场 
     
    自然资源部地质灾害防治技术指导中心派专家组赴云南...

    1 前言

    近年由于常规天然气资源量和产量的下降,特别是在北美洲,非常规天然气得到了高度的重视。一些估计表明,全球非常规天然气资源量(不含水合物)超过30000万亿立方英尺,大约有50%的资源来自页岩气。Julander能源公司的首席执行官Fred Julander认为页岩气(SG)是“自发现石油以来最重要的能源进展”。

    水平钻井技术的进步、水力压裂、相对高的天然气价格(相比2009年之前)和近来在巴内特页岩(Barnett Shale)和美国其他几个页岩气藏的商业成功都使页岩气在美国成为了热门能源,而且页岩气的勘探开发已开始蔓延到加拿大和世界其他几个地区。

    由于页岩气远景的复杂性和广泛性,针对页岩气的应用不能采用普遍用于常规气和煤层气的应用技术,而需专门设计开发工具和方法。多名学者包括Gray等人(2007)和Harding(2008)认为基于确定性解决方案的决议不适用于页岩气开发,因其没有考虑与复杂成藏有关的风险和不确定性,且经常导致过于乐观的结果。

    到目前为止,尽管在北美和欧洲的勘查活动活跃以及近期商品价格下降,页岩气远景分析工作也只完成了极少的部分。商品价格的下降使最高质量远景区的开发至关重要,这些区域的开发不仅最符合公司的利益,并且赋予公司与国外的低成本常规气田(即卡塔尔和沙特阿拉伯相关的天然气)竞争的最佳潜力。Williams-Kovacs和Clarkson(2011)提供了与非常规的远景分析有关的现有工作的回顾,并提供了一种专为页岩气应用而设计的综合的六阶段远景分析及开发评价方法(PADEM)。本文中,作者还展示了一个专门开发用以筛查页岩气远景区并且选择最适合详细分析远景的工具。本文以Williams-Kovacs和Clarkson的工作为基础,致力于远景评价并选择进行更深入分析的远景区的试点位置。

    当前工作的目标是:①开发一种协助页岩气勘探开发阶段的方法和配套的分析工具;②演示已开发技术在加拿大西部致密砂岩/页岩远景区的应用。这项工作的主要贡献是开发与示范一种针对页岩气远景区的严格分析方法。当考虑共存关系时,基于先导试验井输入变量的不确定性,该方法能生成其预测的分布。以前所有的工作一直专注于全域开发方案,然而无法利用勘探开发早期阶段可获取的少量数据快速形成这种全域开发方案。

    2 工具开发

    在这项工作中,开发了一种用于分析页岩气远景的工具。该工具选择使用(以Williams-Kovacs和Clarkson提出的方法(2011)为例的)预筛选的方法。本文将重点放在该工具的开发和应用,分析某一远景区的不同区域,以确定它们是否是适合的试点项目,并描述了图1所示的PADEM工作流程的勘探阶段。勘探阶段的目的是对从更多的详细资料中筛选的远景进行调查,以增加对油藏流动性和碳氢化合物生成能力的了解。在这项工作中,我们对个别类型油井采用概率范围经济学(probabilistic scoping economics)作为勘探标准,以确定该远景区是否适合实行试点项目。表1中完整提供了Williams-Kovacs和Clarkson(2011)详细讨论整体勘探开发方法的总结。

    表1  勘探开发方法概况

    发展阶段

    概述

    靶区筛选

    评估所有潜在的远景区,并选择能提供最好的商业成功机会的远景区

    勘探

    对远景区进行更详细地调查,提高对油藏流体特性和相应碳氢化合物生产能力认识。确定有代表性的试点项目适合的地区

    试采

    继续提高对远景区的认识,集中验证试采区单井的供给能力,评估完井方法

    商业示范

    在项目提交全部资金预算之前,完成开发部分(30%)针对错误的试验结果的测试

    全域开发

    完成全域开发计划,开始制定退出战略

    新的远景/退出

    完成项目详细回顾,评估区域及具体化开发过程中新的远景相关区域。调整和实施退出战略以及任何所需的补充措施

    在这项应用中解析模型比数值模拟更适用,其原因在于应用程序自设置和初始化的时间很短,整合的蒙特卡罗模拟法简单易行,并且在勘探早期阶段不容易获得形成精准的数值模拟所需的详细数据。尽管数值模拟技术已得到改进,但解析方法在工业和文献中依然被大量使用。下文给出了开发工具的关键部分的概要。

    2.1 属性图

    勘查方法最关键的组成部分可能是关键储层、地质力学、岩石物理和地球化学特性的精确属性图的开发。从地质模型、产量不稳定分析(RTA)、压力不稳定分析(PTA)、岩石物理调查等组合中可以推导出这些属性图。这些属性图用于远景的可视化、区块选区以及单一区块的分析。天然气原始地质储量图(OGIP)、Km-h图、压裂脆性图等有助于选择代表性区块以及具备更大开发潜力的区块,甚至高度非均质性区块。区块作为一种评价不同区块远景生产特性的方法,基于地质和岩石物理的观察,比较简单易于操作。采用区块方法不需要针对每个勘探网区块开发一种标准井进行分析,然而通过应用蒙特卡罗法依然解释了其变化性和不确定性。Clarkson和McGovern(2005)采用区块方法评价了煤层气(CBM)远景。通过输入X-Y坐标值以及PetrelTM软件的储层属性Z值可以在Excel中创建储层属性图。随后,数据透视表程序被用于对数据排序,并利用二维绘图应用软件创建属性图。由于早期的岩石物理模型通常利用有限的数据集开发,单一区块在蒙特卡罗模拟中选择不确定的输入数据和参数范围可以解释模型参数的不确定性。这种解释不确定性的方法将在本文所示实例中进行演示。

    2.2 水力压裂模型

    该项工作中,水力压裂裂缝的半长采用Valko(2001)提出的在常规和致密气中应用的简单双翼压裂模型来预测。该模型采用基质渗透率、剪切模量(杨氏模量与泊松比的函数)以及其他储层参数作为输入数据,且如果建模的输入参数不确定,则都必须重新计算每次蒙特卡罗迭代。采用简单的关联(Acm=4xfh)可将裂缝半长转换为与压裂有关的面积。这个压裂模型可能无法代表部分更复杂的页岩气裂缝。为了更好的表示引入到大部分页岩气储层的复杂压裂网,Xu(2009,2010)等人建立了一个更具有代表性的水力压裂模型,该模型将被结合到本次工作中所演示的更新版本的方法中。该区的微地震观测表明,在本文预测的远景区横向双翼压裂的假设是合理的。

    作为所应用的速率预测模型中的关键组成部分必须估算裂缝半长,这一问题将在下面部分开展讨论。水力压裂裂缝半长在随机分析中作为不确定的输入量,其分布主要根据该地区的微地震事件或者其他方法来确定。

     

     

    图1  非常规天然气勘探阶段的勘探/开发方法工作流程

    2.3 速率预测

    Clarkson(2013)提供了关于页岩气井生产分析和速率预测综合全面的概述。在该工作中,我们将页岩气井理想化为一个矩形双孔介质系统,气体从基质岩块流入到裂缝且储层不随着裂缝延展(如图2的概念模型)。该模型忽略了包括体积压裂(SRV)在内的影响,其他作者认为大部分低渗页岩气井在合理的时间内不会发生体积压裂。此外,图2所示的概念模型假设了一个均质的完井——Amborse等(2011)和Nobakht等(2011a)讨论了非均质储层完井的预测。

    在本次工作中,该模型的解决方案首先由EI-Banbi(1998)提出来。人们普遍认为在页岩气藏中占主导地位的瞬时流动状态是从基质到裂缝的线性流。同时,也可能出现一个与水力压裂线性流动相关的线性流动周期,但是通常认为这个阶段持续时间很短,或者被水力压裂清理以及表皮效应所掩盖,而很少可用于分析。本项工作中,我们假设瞬时线性流(从基质到裂缝)之后是边界控制流,该流态与受表皮效应(见等式7)影响的线性流体模型存在早期偏差。压裂段之间的不渗透边界结构导致了边界控制流产生。由Wattenbarger等(1998)首先将早期线性到边界控制流体的假设引入到致密气的应用中,并且该假设被广泛应用于文献和页岩气行业的解析模型。

     

     

    图2  从线性流到边界流的解的概念模型

    2.3.1 瞬时线性流的速率预测

    EI-Banbi(1998)提出通过恒定速率和恒定流体压力来描述瞬时线性流的公式。本项工作中采用恒定流体压力的条件,这也是本文其他部分的重点——该边界条件最接近大部分产生达到最大水位降低值的页岩气井的流动条件。Samandarli等人(2011)采用不同的流体压力迭代方法,对页岩气生产进行分析建模,但是他们表明在大部分情况下采用恒定流体压力的假设就可以了。

    与常用于表征简单横向双翼压裂的裂缝半长(Xf)相比,相关储层面积(Acm)能更好的表示完井措施和增产措施效果以及生成复杂裂缝的能力。因此,在这一分析中,采用相关的储层(气藏)面积(Acm)取代裂缝半长(Xf)。许多业内专家相信由于页岩气藏超低的基质渗透率,复杂压裂对于页岩气的商业生产至关重要。

    无因次时间,tD,Acm,相关储层面积(Acm)依据公式1在恒定压力条件下定义。

                               (1)

    无因次速率,qD,Acm,由无因次时间定义:

                                           (2)

    基于储层特性的无因次速率表达式,如果可获得关于KmAcm估算值,通过公式(3)可确定气体流速。采用不稳定产量分析或者其他的模拟技术可估算KmAcmKm也可以通过实验室技术单独确定。

                                 (3)

    Ibrahim和Wattenbarger(2006)认为线性流的性能受水位下降程度的影响,同时提出水位下降量修正因子(fcp)。此次工作中采用的修正因子(fcp)由公式4给出。

                                (4)

    此处,

     

    Nobakht等人2011a和Nobakht等人(2011b)通过分析中采用校正时间(本次工作未采用)提出一种更严格的校正水位下降量的方法。

    将水位下降量修正因子应用到公式3得出公式5:

                           (5)

    除了水位下降量的修正,这些公式经过进一步修改可直接应用于页岩气井。与致密气井相比,大部分页岩气井在时间曲线的平方根中表现出的较大截距(在致密气井中曲线通常穿过原点),而在流量和时间双对数曲线上页岩气井则呈现出的一半斜率的偏差。多名作者最初认为是裂缝的有限导流能力造成了这种偏差,但是Bello(2009)和Bello和Wattenbarger(2009,2010)认为这种偏差可以通过采用表面效应来更好的解释。Bello(2009)、Bello和Wattenbarger(2009)在恒定流量和恒定流体压力条件下完成了大量的受表皮效应(skin effect)影响的线性流分析,且推导出了恒定流体压力条件下的解析解。在他们的分析中,将表皮效应作为一个常量。Bello(2009)和Bello和Wattenbarger(2009)证明恒定流量情况下表皮是附加量,而恒定流体压力情况下表皮的作用是非线性的。由Bello和Wattenbarger(2009)提出的解析式可以使用下面的近似代数方程:

                    (6)

    从方程(6)可以看出,当tD(t)值大时,包含表皮的项就会变小。

    Nobakht等人(2012)研究了巴内特、马塞勒斯和蒙特利的大量页岩气井(这些气井在相对恒定的流压下产量不断降低),同时得出结论:通常这些页岩气井更多表现出恒定流量的情况而不是恒定流压的情况。作者假设这种意想不到的表现可能是由于Bello(2009)以及Bello和Wattenbarger(2009)提出的表皮模型太过理想化,因此无法代表野外条件。通过假设恒定的表皮效应,模型不能说明由压裂清理、压力敏感地层、变化的压裂导流能力、变化的井底流压、压力相关的流体性质、变化的井筒流体梯度、液体加载等导致的表皮改变。作为这项工作的结果,作者提出了一个可应用于公式(2)的替代表皮修正项:

                       (7)

    包括水位最低量和表皮的影响,公式(1)、(5)、(7)能够利用预测的气体流量,作为时间的函数,在线性流区域可对KmAcm给出独立的估测。

    2.3.2 边界控制流的流量预测

    上面描述的方法适用于有效的储层边界相互接触,边界控制流形成之前。基于图2所示的几何图形,边界控制流紧随着瞬时线性流的末期出现。当外部SRV的影响较为显著时,这一观点较为保守。Clarkson和Beierle(2011)认为如果遇到了其他的瞬时流区,则应采用多重分区的方法,此外,如果多级压裂井需要进行非均质性储层的完井(heterogeneous completion),早期线性流之后不会立刻发生真实边界控制流,且需要更复杂“混合”预测技术。如同下面叙述的,我们选择采用更为保守预测程序,假设线性流之后紧随边界控制流。

    利用公式8计算达到线性流的拟稳态时间(或者是瞬时线性流的结束时间):

                           (8)

    正如图2中看到Ye是压裂到储层边界的距离,计算公式如下:

                           (9)

    多名作者已经提出了页岩气井拟稳态线性流的预测方法。包括Fraim和Wattenbarger(1987),Palacio和Blasingame(1993),Doublet等(1994),Agarwal等(1999)和Mattar和Anderson(2005)认为可采用物质平衡类模拟程序预测边界控制流。Clarkson和Pedersen(2010)将这种方法应用于致密油研究,同时本文也将采用这种方法。公式(10)给出采用物质平衡方法预测边界控制流的生产速度:

                  (10)

    此处qpssi-Linear是边界控制流初始的页岩气流体速度,Pri)pss是边界控制流初始的平均储层压力,且Pwfi)pss边界控制流体初始时井筒流体压力。通过物质平衡计算平均储层实际气体拟压力。对于含有大量吸附气的页岩气开采(application),一般使用Clarkson和McGovern(2005)提出的MBE方法。而在以游离气为主的情况下,则使用定容气藏的常规MBE方法。物质平衡计算需要地质储量和气体特性(比如天然气压缩因子),这两者都是由关键PVT输入量和状态公式(EOS)确定的。

    (a)

    收入总额

    (b)

    收入总额

    扣减

    使用费

    扣减

    使用费

    扣减

    运营成本

    扣减

    运营成本

    得出

    税前运营现金收入(OCIBT)

    扣减

    资金成本补助(CCA)

    扣减

    收入税

    扣减

    加拿大开发费用(CDE)

    得出

    税后运营现金收入(OCIAT)

    扣减

    加拿大勘查费用(CEE)

    扣减

    资本支出

    扣减

    加拿大油气物业费(COGPE)

    得出

    税后现金流(CFAT)

    得出

    生产应税所得

    贴现

    税后贴现现金流(DCFAT)

    生产税率

       

    得出

    应付税款

       

    扣减

    免税额度

       

    得出

    应付净所得税

    图3  现金流分析:(a)现金流;(b)收入税(加拿大税制)

    结合El-Banbi(1998)改进的瞬时线性流的无因次公式和边界控制流的物质平衡模拟方法,可以开发一种综合的预测方法:

    1)        获取Acm(或者Xf)和Km(来源于微地震和/或RTA模拟/已有生产数据或者其他估计)的独立估算值。

    2)        使用公式(1)和(7)作为时间函数计算tD,AcmqD,Acm

    3)        线性流部分的数据利用公式(5)作为时间函数计算qg

    4)        指定排放区(来源FMB模拟/已有的生产数据或者其他估算)。

    5)        使用公式(8)和(9)计算tPSS-LinearYe

    6)        确定

    7)        采用公式(10)通过废弃量(边界控制流)从tPSS-Linear预测产量。

    上面描述的解析模型是假设模型(最小变化)区块内的体积平均值参数是恒量,并从认为是不确定的参数的概率分布中选择一个值。每一次蒙特卡罗迭代将选择不同的值,导致不同的流量预测和不同的主要经济指标值。在许多参数高异质性水平的情况下,存在明显的不确定性,这种不确定性反映在关键输出参数的显著变化。

    2.4 经济模块

    将经济模块与速率预测集成来计算与生产相关的现金流。因为通常行业采用名义美元计算实际(通常的)现金流和名义(现行的)现金流,虽然采用实际的盈利指数计算项目的最低预期资本回收率,且通过不同的通货膨胀率来比较项目。采用图3中的业务流程计算现金流和收入税(加拿大税收制度)。

    该模块中的天然气价格的确定实行了价格操纵,而非价格预测。采用价格操纵表明了项目十分稳定(不论是单独而言还是相较于其他项目),并且不再需要预测极不稳定的天然气价格,该模块中也设置了以价格预测为基础引导经济的选项。

    方法中建立了多个实际盈利能力的指标,包括净现值(NPV)、内部收益率(IRR)和投资收益率(ROI),用来比较项目和公司设定的最低预期资本回收率,同时可给项目进行排序。

    2.5 蒙特卡罗模拟的一体化

    本次工作将蒙特卡罗模拟整合到方法开发中。采用@RISKTM(Palisade Corporation,2010)对关键PVT和储层属性(原始参数)进行概率分布和模拟操作。概率分布的输入变量根据不同项目的数据数量和质量而变化。Clarkson和McGovern(2005),Haskett和Brown(2005)和Harding(2008)认为对数正态分布最能代表PVT、储层和经济特性,因此本文使用了这种分布类型。这些概率分布拟合按P10(低)、P50(中)和P90(高)不同的值输入各个不确定变量。这些输入值可能来自勘探/远景数据、个人经验、模拟数据等。缩减所有输入变量的分布保证每个实现只选择合理的数值(缩减分布将选择少量接近无穷大的数值,从而影响输出变量)。

    上面讨论了@RISKTM输出变量定义的关键经济参数,以及气体速率和累积天然气产量。由于每个输出变量允许量化与项目相关的不确定性,可对其生成一个概率分布,以便做出与远景选取和开发有关的明智决策。

    通过在x轴上找到相应的最低预期资本回收率时的位置,向上垂直移动至曲线处,然后再水平投影到y轴,这样可以从累积概率分布计算出超过设定最低预期资本回收率的概率。用1减去y轴上求出的值,得出超过最低预期资本回收率的概率。这个方法在本文中将作为范例进行演示。

    在这一应用中(如在孔隙度和渗透率之间),采用了拉丁超立方体抽样,如果有必要的话,还可合并相关性(如孔隙度与渗透率)。典型的多相(气+水)页岩气/致密气应用的主要参数如表2.3所示。在某些情况下,参数的依赖关系可使用行业普遍接受的经验模型进行解释,而在其他情况下会使用来自现场数据或者估算得到的基于方向的相关性(如较高的正相关关系)。例如,与压力有关的渗透率(绝对的渗透率比值)使用Yilmaz等人(1991)的方法可与储层压力和岩石力学特性关联。相反,束缚水饱和度与孔隙度密切正相关。可能的参数关系如表2所示。

    蒙特卡罗模拟运用了一个类似于Clarkson和McGovern(2005)使用的煤层气气藏远景分析的方法。

    表2  基本参数、可能的相关性和参数关系

    基本参数

    可能的相关性

    关系

    有效厚度/英尺

    孔隙度/%

    粒径,有机质

    适用于某些情况下和正相关情况的实证模型

    初始含水饱和度/%

    孔隙度

    高度正相关

    束缚水饱和度/%

    孔隙度

    高度正相关

    基质渗透率/毫达西,初始状态

    孔隙度,有机物

    适用于某些情况下和正相关情况的实证模型

    基质渗透率/毫达西,初始比

    储层压力,力学性能

    野外/岩心数据经验曲线

    相对渗透率

    含水饱和度,束缚水饱和度

    野外/岩心数据经验曲线

    初始储层压力/磅/平方英寸

    深度,渗透率(超压)

    气压梯度

    储层温度/℉

    深度

    温度梯度

    天然气比重

    朗缪尔体积/标准立方英尺/吨

    容积密度

    来自岩心/岩屑的线性关系

    朗缪尔压力/磅/平方英寸

    体积密度/克/立方厘米

    流泄区/英亩

    含气量/标准立方英尺/吨

    TOC

    正相关

    井眼半径/英尺

    表面

    增产效果

    高度正相关

    压裂总半径/英尺

    剪切模量(+),渗透率(-),有效厚度(-),井眼半径/英尺

    变化—见括号中相关方向

    井底流压

    井眼长度

    高度正相关

    3 该方法应用于远景勘探

    本文中开发的方法广泛应用于SG远景将其分成区块进行分析的目的,以确定是否适合作为一个试点项目。由于SG试点和开发项目成本高,且其详细分析需要大量数据,页岩气远景勘探至关重要。

    对于远景勘探应用而言,其方法的选择以当前远景数据和模拟数据相结合为基础。理想情况下,对于关键PVT和储层参数情况良好的估计,作为空间坐标的函数可用于远景勘探。如果事实并非如此,可以对模拟气藏或者其他数据源进行估算以获取数据,同时分析该方法带来的不确定性。

    假定整个远景区PVT和其他储层特性不变,输入数据可用于生成主要储层特性图。关键生产指标图如OGIP和基质渗透率乘以可以开发的净投入(千米/小时),可用于区块的选择。区块的选择基于区域类似的关键生产指标的值。对页岩气储层而言,压裂的指标,如压裂指数或脆性也可能用于区块选择,同时许多作者表明建立复杂裂缝网的能力对于页岩气商业开采至关重要。

    选择区块后,开始进行蒙特卡罗模拟,按照P10、P50、P90的概率预测和可以开发累积产气的区块,且结合使用关键经济指标的分析来确定区块能否适合一个试点项目。其他因素比如公司的经验,企业和商业策略,可用的资源和基础设施等都将纳入评估,以便为公司以及股东们确定哪些区域可以作为最佳试点选项作出明智的决策。

    远景勘探方法工作流程见图4所示。

    4 采用两段页岩开发模型的样本示例

    为了进一步说明该方法的应用,对加拿大西部的某处致密砂岩/页岩(假定没有吸附气体)远景区的两段进行了分析。在之前的研究中,PetrelTM开发的远景地质模型采用可用的岩石物性、储层和生产数据。图5所示研究区域内4口井的三维孔隙度模型和孔隙度相关的钻/录/测井记录。在该区域,存在两处可获益的产气水平井段(井段3和井段4)。

     

    输入数据

    关键储层属性的填图属性

    PVT,其他储层和水力压裂属性

    生产数据

    经济投入

    区块选择

    根据OGIP或者其他关键属性确定区块

    蒙特卡罗模拟

    模拟输出

    P10、P50、P90的概率预测和累积产气量

    水力压裂运行情况

    经济参数

    可行的商业区块标志

    其他

     

     

    图4  远景勘探方法的工作流程

     

     

    图5  三维孔隙度模型和孔隙度相关的测井

    模型开发期间这个开发区拥有11口垂直井,2口倾斜井和4口水平井。最初钻完成垂直井,紧随其后的是开始于2008年的水平井。Clarkson和Beierle(2011)在该区选择一系列井进行不稳定产量试井(RTA)。模型开发中使用的水平井的总结显示在下面表3中,同时在图6中(在下面描述)该区域的天然气原始地质储量(OGIP)图上显示了井的近似轨迹。

    表3  研究区水平井概况

    井名

    井向

    进入层位

    完井方式

    1号井

    水平

    井段4

    尾管注水泥

    2号井

    水平

    井段3

    自膨胀封隔器

    3号井

    水平

    井段4

    自膨胀封隔器

    4号井

    水平

    井段4

    自膨胀封隔器

    所做的分析主要集中在大部分是水平井的井段4。为了简化分析,采用孔隙度下限为4%,通过Excel加权平均井段4层位,将PetrelTM多层模型转换成一个单层模型。这一平均化过程是为了完成对基质的孔隙度、初始含水饱和度和渗透率的处理。利用孔隙度下限值还可以计算总有效收益和毛净收益(有效收益假设包括所有孔隙度下限值以上的层)。图7a和图8a显示了OGIP和Km-h属性图。

    模型采用的网格大小如表4所示。在整个开发过程中假设为常量的PVT、储层和生产参数如表5所示。

    表4  网格属性

    网格属性

    数值

    网格尺寸

    135×129

    区块长度,X/英尺

    49.76

    区块长度,Y/英尺

    49.76

    网格区块面积/Ac

    0.057

    对于这种情况,人们认为井筒流动压力(pwf)为常量1750磅/平方英寸,接近开发区水平井最初的井筒流动压力。随着时间的推移井筒流动压力降低,后期模型中压力驱动力低于开发井,模拟气率并不乐观。这种情况下,在可获取日常生产和流动压力期间内,平均两个收益井的流动压力大约是1550磅/平方英寸,因此到开发后期之前,这种假设的影响并不很明显。在实际勘探中,该地区还没有投入生产,由于我们不需要将可用的生产数据与模型匹配,而是采用实际的流动压力估计值尝试得到一个准确的潜在生产能力估计值,所以这种假设的影响不是一个值得关注的问题。

    表5  PVT常数、储层和生产投入参数

    参数

    PVT参数

     

    气体比重

    0.648

    N2/%

    0.46

    CO2/%

    0.2

    H2S/%

    0.0

    温度/℉

    166.5

    Cw/磅/平方英寸-1

    2.9×10-6

    Cr/磅/平方英寸-1

    5.6×10-6

    VL/标准立方英尺/吨

    N/A

    PL/磅/平方英寸

    N/A

    储层参数

     

    Pi/磅/平方英寸

    3500

    排放面积/Ac

    80

    生产参数

     

    Pwf/磅/平方英寸

    1750

    rw/英尺

    0.3

    3个区块中假设关键属性的变化情况如表6所示。各属性的数值是每个区块的各个网格值的算术平均数。由于基质渗透率是蒙特卡罗输入量,且利用基质渗透率值可计算总压裂半径(虽然也可使用压裂分析模型在每次迭代时作为基质渗透率函数计算总压裂半径),故给出了一个基质渗透率值以显示区块之间总值的变化情况。

    表6  储层变量和水力压裂输入参数

    参数

    区块1

    区块2

    区块3

    储层参数

         

    有效厚度/英尺

    102

    74

    58

    孔隙度/%

    7.1

    6.5

    6.0

    Sw/%

    18

    15

    16

    Km/毫达西

    0.0084

    0.0079

    0.0077

    水力压裂参数

         

    剪切模量/磅/平方英寸

    2×106

    2×106

    2×106

    总压裂半径/英尺

    1432

    1477

    1489

     

     

    图6  研究区地质储量图呈现近似水平井轨迹

    4.1 区块选择

    利用从PetrelTM多层模型开发的单层模型,其单层等量地质储量如图7a所示。根据类似颜色为代表的区域具有类似地质特征和岩石物理性质,通过视觉观察可选择区块。虽然已知气藏具有高度的横向非均质性,可以看到关键的地质和岩石物理性质明显凸出部分。该图形显示了更复杂的异质性模式的情况,需要更多的区块并且可能有必要用区块代表具有相似属性的不连续块段。图7b显示基于天然气原始地质储量选择的区块远景区。在计算天然气原始地质储量时,虽然该远景区吸附气体量很容易被包含其中,但还是假设其可以忽略不计。

     

     

    图7  地质储量图:(a)地质储量;(b)选区

    从图7b可以看出选取的三个区块中,区块1具有最高的天然气原始地质储量(红色和橙色),区块2具有的地质储量(光和暗绿色)次之,区块3具有的地质储量(紫色和蓝色)最低。从这幅图中可以推断出区块1将有最理想的属性,因此可能具有最高的产量,而区块3产气物性最不理想,因此可能具有最不理想产气量。如同气藏地质储量图(图7)一样,如果绘制Km-h图我们也可以分辨出三个相似的区块。此次应用区块选区采用的天然气原始地质储量图和Km-h图作为代表资源的程度/密度和储层特性的两个要素,这是工业上常用的评估致密砂岩和页岩远景好坏的关键因素。区块选区的属性根据不同项目而变化,取决于驱动特定资源类型远景的关键要素。

    对于这种情况,假设简单的水平双翼压裂(如所使用的压裂模型所假定的)就足够了,因为微地震数据对同一区域的补充水平压裂井的解译说明复杂程度较低,如果不是水平情况,则进行压裂(图9)。采用水平和垂直观察井用以观察,同时采用双阵列处理会产生一个好的数据集。一般情况下,各个阶段仅出现一个水力压裂裂缝。水力压裂裂缝通常选择北东-南西方向,与加拿大西部沉积盆地(WCSB)部分最大水平应力方向一致。

     

     

    图8  Km-h图:(a)Km-h;(b)选区

    通过比较图6与图7b和8b可以看出在开发区所有水平井部分或全部在区块1范围内。因为这个原因,剩余的分析还将在区块1中开展。对区块1区域的水平井的预测情况而言稍微乐观,因为这些水平井水平延伸超出区块1区域进入地质储量和Km-h更低的区域(该区水平井采用恒定的流体压力与(Pwf)i相比将获得相反的影响)。

     

     

    图9  根据微地震数据解译的研究区内水平井水力压裂裂缝几何图形

    4.2 经济分析

    分析假设只有天然气价格是变量,而所有其他经济参数都保持常量。表7列出了其他主要经济参数的值(基于Magyar和Jordan的估算(2009))和表8介绍了主要的专利权使用费、税和贴现参数。

    在本文的分析中,净现值(NPV)作为重要的收益经济指标且最低资本回报率为0。

    分析远景的工作流程图如图4。

    表7  资本和运营成本参数

    参数

    土地成本

     

    租金/美元/亩

    2500

    代理费/美元/亩

    50

    单井成本

     

    钻井/百万美元

    1.5

    完井/模拟/百万美元

    2

    配套设施/管道/百万美元

    0.35

    储层表征

     

    地震/百万美元

    0

    测井/百万美元

    0

    提取岩心/百万美元

    0

    其他/百万美元

    0

    运营成本

     

    固定成本/美元/月

    5800

    可变成本/美元/千标准立方英尺

    1.25

    表8  使用费、税收和折现率

    经济参数

    费率

    使用费率

    20%

    税率

    30%

    实际贴现率

    15%

    名义贴现率

    18.45%

    通货膨胀率

    3%

    4.3 蒙特卡罗模拟

    在区块选择之后,本文进行了蒙特卡罗模拟研究。蒙特卡罗模拟中,基质渗透率(km)和页岩气价格不断变化,而所有其他的PVT、储层参数和经济参数保持不变。为了更好地进行说明,我们选择了将“不确定”的输入变量的数量显著限制在基本控制远景的油藏性能(储层渗透率)和经济情况(天然气价格)。基于P10、P50和P90值按照对数正态分布模拟参数。在大多数的勘探情况下,许多参数都是不确定的,可以通过这些参数的概率分布(见表2)来定义。对于需要使用概率分布进行定义的一些关键参数,可通过评估给定区块内重大变化的属性图来直接确定,或用更严格的统计技术,如采用区块内部数值计算变异系数(Cv)。由于基质渗透率是基质流动的主要控制要素,以及未来商品价格造成的天然气价格的高度不确定,针对这种情况,我们选择基于视觉观察的基质渗透率。

    基质渗透率按照P10、P50和P90的值计算如下。通常情况下,可以通过岩石物理模型中的参数值拟合分布来生成概率分布,但是因为我们处理的是远景的早期评估,因此我们采用了替代的方法,即最大限度提高模型获取的不确定性来解释其他早期参数估算无法获取的变化性。如果需要,对其他不确定参数也可以使用相似的方法。

    P10——区块1中比第十百分位值的基质渗透率低20%

    P50——区块1中的基质渗透率值居中间数

    P90——区块1中比基质渗透率的九十百分位值高20%

    表9中定义了2个输入变量的分布。将模型内部不确定参数合并关联(见表2)也很重要。虽然孔隙度和渗透率之间的相关性被加入到原始岩石物理模型(幂律相关),并且压裂半径与剪切系数(正相关)、基质渗透率(负相关),净收益(负相关)和压裂模型井眼半径(负相关)相关,但是出于演示的目的,本文对这一方法进行了简化,使蒙特卡罗模拟中的主要变量之间没有相关性。由于压裂半径取决于基质渗透率,压裂模型必须在每次迭代时重新计算。气体流量,累积产气量和净现值被定义为@RISKTM输出变量。

    本文进行了5000次蒙特卡罗迭代,以确保蒙特卡罗输入变量充分覆盖样本空间。要求覆盖足够的样品空间,是为了确保每个模拟输入相同参数运行时,能得出同样的结论。出于演示的目的,用上述方法获得的迭代数并不是最优化。但是,通过将无限大(非常大)的样本输出分布与减少样本数量的输出分布比较,同时寻找要求充分重复“已知”输出分布的最小值,可以获得优化的迭代数。当进行多个模拟时,优化处理可用于减少处理时间和容量。

    4.4 结果

    图10显示了区块1中单口气井的确定产气量和累积产气量预测。这个“确定性”的基质渗透率的值来自于表9所示输入分布的斯旺森平均值(SM),假设这个值代表区块收益的平均水平(静态平均Km=0.0095毫达西)。虽然Bickel等人(2011)指出了斯旺森平均值(SM)的缺点,但它仍然被广泛地用于工业,因此在这种情况下还将使用。此外,斯旺森平均值在输入分布的平均值的5%范围内(使用@RISKTM计算),因此认为在这个例子中的平均值是准确的。另外,可以使用另一个估计的平均值(即分布平均值、区块值的算术平均值等)。图10a显示的产气速率与时间半对数图以及累积气体的产生与时间的笛卡尔曲线,而图10b显示了产气速率和时间的对数分布图。

     

     

    图10  开发模型情况下的确定速率预测:(a)产气速率和时间、累积产气量和时间的半对数;(b)产气速率和时间的对数关系

    图11显示了产气速率与时间的半对数图,图11b显示一个产气速率与时间的对数图和图11c显示预测(约14年)最初5000天累积产气与时间的笛卡尔曲线。

    通过比较图10和图11,可以再次看到确定性预测与P50概率预测相比,具有更大的IP,持续的生产速度和累积产气量,表明确定性预测是比中位数情况稍微乐观,并且明显远超过P10的情况。这些结果再次支持使用概率分析取代非常规应用的确定性分析。

     

      

    图11  开发模型情况中概率速率预测:(a)产气速度和时间的半对数关系;(b)产气速率和时间的对数关系;(c)累积产气量和时间

    随后,P10、P50和P90产量预测与区块1内水平井可获取的生产数据进行对比,以测试开发方法的稳健性和准确性。在这个比较中,由于完井的复杂性,只有井3和井4可用,而井1表现不佳,且井2在此次分析区块外部。井3的产量被缩减了30天,以便使该井产量自然下降的初始时间与概率预测的一致(指修正井3)。生产的前430天的对比曲线如图12所示。

    如图12所示,两口井的生产数据(修正井3和井4)普遍落在P10和P90之间(使用@RISKTM生成的预测)。除了生产的前20天和第300天左右时的大约20天两个时间段(模型没有指出的操作问题导致的结果)外,约80%的数据点如预期处在P10和P90预测之间。初步预测产量可能更高,因为它不考虑压裂清理干扰、启动效应等,该模型增加了表皮效应来提高与IP的匹配程度。但是,在真正的勘探情况下表皮效应的大小无从得知,这是因为无法获取产气远景区域的数据且需要将其作为不确定的输入量以最大限度地提高模型的准确性。

     

     

    图12  3号井和4号井生产数据和概率速率预测的对比:(a)产气速率和时间的半对数关系;(b)产气速率和时间的对数关系;(c)累积产气量和时间

    虽然这不是一个令人满意的统计样本,只有一个关键属性(Km)被认为是不确定的,但结果令人鼓舞。图13显示了净现值的增加的累积概率分布,直方图和回归系数托那多图。图13a再次显示超过最低预期资本回收率概率计算的累积概率分布图。

    从图13a可以看出这个模拟平均净现值为53万美元,可能超过最低预期资本回收率的50%。然后,可将平均净现值和超过最低预期资本回收率的概率与相同远景的其他区块,以及与其他潜在远景的区块进行比较,从而确定哪些远景区域可提供最好的经济成功机会。这一分析显示了积极的NPV平均值和超过最低预期资本回收率的适度概率。基于这样的分析,可以得出结论:区块1的样品远景对于试点项目是极好的备选。这一分析支持了该地区的开发,但是这一测试中所采用的天然气价格网格假设对其结果影响极大。图13C中托那多图表明天然气价格对净现值带来的影响最大,基质渗透率给净现值带来的影响其次(区块1中最小的基质渗透率变化的结果)。这表明假设较高的气体价格(比如该区水平井钻探时期的气体价格)将提高远景的可取性。从图13b直方图可以看出模拟中大部分的净现值在300万美元和350万美元之间,众数等于-1.5万美元,相当于平均数53万美元左右。

     

     

    图13  开发模型应用NPV法得出的经济结果:(a)累积概率分布;(b)柱状图;(c)回归系数的龙卷风图

    此分析程序可在在开发区的其他2个区块内完成,以协助选择最适合公司的试点项目的位置。2号和3号区块的填图属性的直观观察(图7b和8b)表明,这些地区情况没有区块1理想,因此在本次分析所使用的气体价格假设中可能不适合作为试点项目。

    5 结论

    在本文中,开发了一种方法理论和基于excel的方法以协助页岩气和致密砂岩气藏的勘探。这个方法包含了来自不同来源的映射属性、一个用于估算水力压裂半径的简单的压裂模型、目前应用于页岩气井开采的速率预测技术、计算关键盈利能力指标的经济模块以及解释非常规资源中内在的风险和不确定性的蒙特卡罗模拟。本文所描述的方法和工具可被工业界用于评估远景区域内的各个区块和选择适合试点项目的地区。该方法较为严谨,以岩石物理、地质和现在产业应用的分析储层模型为基础,且通过重建现有实例的油藏动态来证明其准确性。由于不需要建立复杂的数值模型和详细的开发方案(所需数据是在开发早期通常无法获取),这种方法既简单又高效。

    感谢代金友副教授对本文提出的宝贵意见。本文受中国地质调查“地学情报综合研究与产品研发”(121201015000150002)项目支持。

    资料来源:Williams-Kovacs J. D., Clarkson C. R. A new tool for prospect evaluation in shale gas reservoirs. Journal of Natural Gas Science and Engineering,2014,18(5):90-103.

    一种用于页岩气藏远景评价的新方法

        在迷你系列“页岩气神话和现实”三部曲中的第二部分里,我想论述的问题是围绕页岩气这一新型资源的一系列核心问题中的一个:页岩气究竟是便宜还是贵?上周,在安得烈蒙特福德的“Bishop Hill”博客上的激烈讨论中,我声称页岩气是贵的,而非便宜。这一观点引发了一系列长篇大论的怀疑和反对。
        相反观点认为,页岩气便宜而且很便宜。因为使用页岩气还可以节省许多钻探和输送页岩气所需的钢管的费用,以及制造钢管所需的人工费。
        根据前期对页岩气成本所做的一些调查,我认为保本价格可能在4美元到8美元/千立方英尺之间(每一千立方英尺的价格大致相当于每百万英制热量单位的价格,即美元/百万btu)。一个应注意到的极其重要的情况是:如果是在欧洲,目前4美元到8美元的价格是便宜的;如果是在日本则是非常便宜的;但是如果是在北美,这个价格很贵(图1)。因此,对全生命周期生产成本简单地做出贵或便宜的结论不是一个明智的方式。


    图1  全球天然气年平均价格的变化(据英国石油公司)

    北美页岩气开发、福岛核事故和液化天然气的短缺是推动价格变化的几个主要因素。右边的灰色栏显示美国页岩气保本价格的最小值、最大值和中值(据贝克研究所Ken Medlock)。如果是在日本,这3种价格都非常便宜,但在美国却会被认为是很高的价格,基本上仍然高于美国当前的现货价格(每百万英热量单位3.85美元)。
       

        在对已公布的页岩气全生命周期盈亏成本估计数字进行更详细的分析前,这里将首先非常简要地概述美国页岩气生产的历史和页岩气产量惊人的增长速度。因为了解过去,有助于理解现在,并有可能预见未来。

       一、美国页岩气产量

        图2显示了David Hughes编制的美国页岩气产量增长的历史。可以看出,美国页岩气主要产自Haynesville、Barnett、Marcellus和Fayetteville这4个地区(图2和图3)。页岩气产量已经从2007年每天约40亿立方英尺增长到2012年每天约260亿立方英尺!页岩气产量的增长对美国天然气总产量有实质的影响,美国天然气总产量从1995年~2002年每月2万亿立方英尺的稳定水平增长到2012年每月2.5万亿立方英尺(图4)。一些人认为这是美国向能源独立的迈进,然而也有人认为受成本和递减率的影响,这可能是一个能源泡沫。


    图2  2000~2012年美国各地区页岩气产量累计
    (产量从21世纪初的停滞状态到现在惊人的每天250亿立方英尺)

    图3  美国和加拿大主要页岩气和石油产地(资料来源:美国能源信息署)

    图4  自1993年以来美国不同类型天然气产量的变化(基于美国能源信息署数据制作)
    2006年数据系列的异常是因为页岩气数据的引入。在此之前,页岩气产量包括在传统的干气产量中。2011年末,页岩气占美国天然气总产量的32%。据David Hughes的调查结果,页岩气产量到2012年末已占美国天然气总产量的40%(图2)。

        美国页岩气产量的超常增长并不是偶然发生的。从2003年开始美国钻探工业进入超速发展期,到2008年夏天钻探设备总数是2003年的两倍多,达到2000台(图5)。大多数钻探设备是开采天然气的,把图5、图2和图4结合起来,就会明显地看到越来越多的钻机把开采页岩气作为目标,这是以常规天然气产量的减少为代价的,其维持多年的稳定水平从2007年开始下降(图5)。2008年8月金融危机爆发后,钻探活动急剧减少,因为当时所有人都在评估今后全球金融体系的格局。美国的钻探活动在2009和2010年迅速恢复。随后又开始了从页岩气开采向页岩油开采的大转移(图5)。金融危机以前,大约有400台石油钻机,1600台天然气钻机。现在是大约1400台石油钻机,400台天然气钻机。


    图5  美国贝克休斯油田服务公司提供的美国钻机数与美国能源信息署提供的天然气总产量对比


        1995年以来,美国钻探队伍的规模已经从大约700台钻探设备增长到2000台。在很长一段时间里,美国增长的钻探工作量仅仅能够弥补减少的工作量,从而维持总体工作量不下降的局面,直至钻探行业把目标转向页岩。与开采非常成熟的常规天然气相比,新的页岩气开采项目显然给公司带来了更高的产量回报。是否页岩气开采也能给公司带来更高的经济回报尚待分晓。有证据显示天然气总产量在过去的24个月里始终维持在一个稳定水平上。
        油气公司从开采天然气转向开采石油的原因很简单。页岩气的生产过剩使价格从2008年6月的超过13美元/千立方英尺下降到2012年4月的不到2美元/千立方英尺(图6)。在相同的时间段里,石油价格要高得多,从而激励了页岩气开采向页岩油开采的大转移。不到2美元/千立方英尺的天然气价格对于美国消费者和经济增长来说是一个异想天开的数字,但却导致了埃克森美孚公司首席执行官Rex Tillerson的著名宣告:“今天我们的财产全部化为乌有”。他在纽约外交关系理事会的一次发言里说,“我们没赚到钱,账面上都是红色赤字。”


    图6  亨利中心天然气现货价格(资料来源:美国能源信息署)


        欧洲和亚洲各国政府对美国的低价格羡慕不已,因为福岛核事故和2012年液化天然气供应量的下降导致了国际市场上天然气短缺,使价格飙升到历史最高水平,同时也导致了不同地区天然气价格的巨大差距(图1)。美国天然气价格下降到小于2美元/千立方英尺的情况制造了一个错觉,即页岩气是便宜的。但实际上是生产过剩使得美国天然气价格大跌至公司可以获益的水平以下。这种情况对社会是不利的,因为社会的生存有赖于可以赢利的能源工业,但社会没有意识到这一点。

       二、页岩气的实际成本是多少?

        在北美,页岩气的全生命周期实际生产成本是多少?没有单一的和简单的回答,有大量的变量需要考虑:
         1. 气井的最终回收率(EUR)
         2. 钻井的沉没成本
         3. 土地成本
         4. 管道、加工设备和运输至市场的成本
         5. 税和权利金
         6. 付出利息和利率
         7. 公司日常管理费
        不同州之间在某些变量上的差异是很大的,例如税制,北美和欧洲之间税制差异更大。在以下给出的各种各样的评估结果里,在某些情况下还不清楚是否使用了全生命周期经济评估法或“远期点”(point-forward)经济评估法。下面是莱斯大学贝克研究所能源研究中心资深总监Ken Medlock的评估结果。
        “一些井每千立方英尺只投入2.65美元成本就可以盈利,其它的需要投入8.10美元……中值是4.85美元。”


    图7  Ruud Weijermars在其学术论文中指出,美国不同产地的页岩气最低价格在每千立方英尺4到6美元之间。

    左图:2007~2011年美国主要的页岩气产区的产量;右图:美国主要页岩气产区的盈亏平衡点边际价格

        图1中已标示了上述3个价格。图7引自一项学术研究的结果,显示最低成本为每千立方英尺4到6美元。
        页岩气分析和评论家Arthur Berman在其2012年的论文中指出,按8%的贴现率计算,美国Barnett、Fayetteville和Haynesville三地生产的页岩气的全生命周期成本超过了8美元。Berman认为公司一直在高估最终回收率,因此相对成本实际上比公司公布的要高。按天然气的现价计,美国所有页岩气产区都是非盈利的。
       一家主要石油公司的前首席执行官说:一切都取决于页岩气的价格。如果价格在每千立方英尺4~5美元,则页岩气产量会很高。
        总之,情况似乎很明显,如果价格低于每千立方英尺2美元,则整个页岩气行业就无利可图。现有的证据表明,贝克研究所Ken Medlock对页岩气价格的估计是比较合理的。各个公司的利润率取决于他们总体业务的质量。只有时间能够证明Weijermars和Berman估计的高价格是否会因页岩气产量低于当前的预测量而实现。但这是一个极端不稳定和活跃的行业,效率一直在提升,因此对成本,尤其是今后成本的预测极其困难。

       三、美国页岩气的未来

        有迹象显示,美国页岩气产量增长速度的急剧放缓已经开始导致美国天然气总产量增速放缓。在过去24个月内,美国天然气总产量始终稳定在一个水平上可能就是证据(图5)。令人惊奇的是开采速率下降之后,产气量继续上升。部分是由于开采与输送系统之间的时间滞差。气井需要等管道接通,任何时候都有储备的新井闲置在那里等待连接管道。而且钻井效率也提高了,水力压裂和完井技术也正在改善,类似宾夕法尼亚州Marcellus的一些地区,生产的热点地区正在被发现和钻探。以下这段话引自“傻瓜投资指南”网:
        根据贝克休斯的数据,相当多公司都在采用垫块钻井方法。因此北美地区2013年第三季度的钻井工作量比2012年第一季度高出15%不应该感到惊讶。
    根    据美国能源信息署的数据,在美国北达科他州的巴肯页岩区,今年新井每台钻机的日产油量预计将提高将近一倍。去年的日产油量略高于250桶,而今年12月的日产油量估计可达到496桶。同样,在德克萨斯州的老鹰堡页岩区,新井每台钻机的日产量预计可从去年的约200桶增加至今年的413桶。
        这样一来,油气的供求似乎有可能实现平衡。400台钻机能够生产出足够多的石油以弥补天然气产量的下降。美国对廉价天然气的需求已经提高了,还需要从加拿大进口少量天然气。在现在市场条件下,很难想象还会有钻机从开采石油转向开采天然气。因此,合理的预测似乎是:石油产量的增长可能会停下来,甚至在中期可能逆转,这种情况将持续到天然气价格回升到某个水平上,足以使上述在巴肯和老鹰堡页岩区开采石油的钻机又转向开采页岩气。
        页岩气开采效率的惊人增长速度肯定会在达到某个水平后放缓。从长远看,产量将取决于能够发现多少高效和高产的气井。只要新井的质量继续高于老井,那么页岩气的奇迹似乎将永远继续下去。
        根据上面的分析,美国天然气价格可能将继续上涨,并在短期内升至每千立方英尺5美元。市民可能会因此而抱怨,但应该认识到只有把价格维持在这个水平,页岩气开采业才有可能至少在一段时期内繁荣发展,美国经济才会受益。按照5美元/立方英尺的价格,那些拥有低成本、高效率气井的公司可以赚许多钱,而那些只有高成本、低效率气井的公司则可能无利可图。

    美国页岩气的真实成本是多少?

     

    广西防城港地区地层多为碎屑岩,地下水资源匮乏。农村地区则更多饮用3-5米深的敞口井,在干旱的冬季许多井都会干涸,季节性缺水严重。2015年以来,中国地质调查局武汉地质调查中心承担的“防城港地区水文地质工程地质调查评价”子项目开展了江平幅1:5万水文地质调查,在面上调查工作基础上,重新优化布设了水文地质钻探,于8月27日全面进入施工阶段。截止9月26日,完成水文地质钻探总进尺860.79米,成井13口,在已完成的水文地质钻孔中,多数出水量在1.5-3.5立方每小时,水量匮乏。

    江平13号(JP-13)钻孔位于侏罗系碎屑岩地层内,地质构造较为特殊,在钻进至20米处即有水自流而出,在打穿73米一处断裂破碎带后,水量急剧增大,经三角堰测量达55.3立方每小时,水质清澈,直至150米终孔时,水量未见衰减。这在整个防城港地区也是极为罕见,对周围两个村庄的安全供水具有较大的现实意义。“这么大的水量、这么清澈甘甜的井水,在防城港地区绝对是绝无仅有!”,每一个路过江平13号钻探现场的人无不发出如此的赞叹,连在本地区从事水文地质工作40余年、参加过1:20万水文地质区调的广西水文队黄显坤老人也如是说。

    目前,项目组对该井进行了止水并取样,待抽水试验、水样分析测试及水源地综合评价完成之后,将与地方政府合作,解决当地实际饮水困难,实现地质调查服务社会经济发展的目标。

     

    江平13号(JP-13)施工现场

    江平13号(JP-13)施工现场

    广西防城港地区探明重要水源地

    中国地质调查局青岛海洋地质研究所日前组织专家对南海三维多道地震调查获取的野外资料进行验收,专家组对“海洋地质九号”船首次自主开展的双源双缆地震作业取得的高品质资料给予了充分肯定,一致通过验收并被评为“优秀”级。

    外业项目负责人介绍,项目组历时1个月时间,在南海开展了三维多道地震调查连续作业。本次三维多道地震共完成测线120条,累计工作量2833千米(满覆盖2626千米),覆盖面积220平方千米。目前,三维地震数据处理工作正在有序进行,从现场预处理完成的地震剖面来看,此次采集的地震数据信噪比高,浅部地层和断裂显示清晰,在峡谷的两侧能够清晰的识别出BSR,且BSR沉积层的厚度可达0.5s,具有较高的调查目标勘探价值。完全满足项目组后续开展研究工作的需要。

    历经4年多的发展,青岛海洋所多道地震调查队伍从无到有,技术能力日臻完善,继2021年实现8000立方英寸大容量气枪阵列/12千米缆长二维多道地震作业能力后,再次升级为双源双缆三维作业。本项目的成功实施,开创了青岛海洋所自主实施三维多道地震的先河,标志着青岛海洋所多道地震调查能力再上台阶,为后续开展斜缆、垂直缆等作业方法研究积累了宝贵经验。 

     

     
    青岛海洋所多道地震调查能力迈入三维时代

    在第59个学雷锋纪念日到来之际,为传承雷锋精神,大力弘扬奉献、友爱、互助、进步的志愿精神,3月4日,中国地质调查局地学文献中心团委组织开展以“志愿服务行 一起向未来”为主题的冬奥会志愿者分享活动。地学文献中心两名借调至奥组委服务团队的团员青年袁晓诺、童铃青分享了冬奥会志愿服务期间的见闻感受。

    分享会上,袁晓诺讲述了从2021年11月中旬借调到冬奥组委国家游泳中心场馆运行团队的工作经历。水(冰)立方作为开赛最早、赛程最多、耗时最长的冬奥场馆,袁晓诺及其所在团队, 1月起停休,且在冬奥会期间平均每日工作时间超过16小时,周末和节假日均未休息,负责赛前讲解、接待和赞助商接待空间的布置工作。曾给多个到访团队进行水冰转换、冰壶相关内容讲解,并获得充分肯定。她表示,在志愿服务期间,让她增长了见识、开阔了眼界、提升了技能、积累了经验,未来将带着奥运精神投入到工作当中。

    童铃青从2021年11月底借调至冬奥组委机场运行团队从事注册业务领域工作。从1月3日进入闭环,到3月冬残奥会保障结束后移出闭环,童铃青身穿二级防护服,在闭环内全力以赴保障冬奥会、冬残奥会。首都国际机场是涉奥人员入境第一站,服务保障质量将直接影响涉奥人员对北京冬奥会的第一印象。“国门第一岗”“北京冬奥会初印象”“奥运会利益相关方的必经之路”这是激活中心的特殊标签,也是童铃青下定决心守护的责任与荣誉。在过去的一个多月时间里,她站好了“国门第一岗”,为入境的涉奥旅客提供了热情、专业的注册激活服务,她表示,“作为冬奥志愿的‘大管家’,我们一起出色地完成了志愿服务任务,用自己的实际行动,奉献冬奥,兑现了庄严的承诺!”

    她们用笑容和奉献温暖了这个冬天,用专业素养和细心耐心为冬奥会顺利举办贡献着力量。她们在为冬奥会奉献的同时,个人也得到了成长,收获了满满的成就感。团员青年纷纷表示,要以实际行动践行志愿精神,充分发挥学雷锋志愿服务站的作用,依托业务工作积极行动,推动地学文化传播,做好地学文献资源的守护者和地学科普知识的传播人。

     

    袁晓诺在为志愿者团队做培训

     

    1月21日0点 童铃青激活中心第一班岗

     
    地学文献中心团委开展冬奥会志愿者分享会

    核心提示

    搞清楚“隐患在哪里”是我国地质灾害防治工作的核心问题之一。2019年,自然资源部部署开展了基于综合遥感技术的地质灾害隐患识别示范工作,力争3~5年内掌握全国地灾隐患“底数”。

    去年,首轮“国家级”大规模地灾隐患综合遥感识别将全国地灾高中易发区8450处疑似隐患“揪”了出来,其成果得到充分认可。但专家也指出,地灾隐患综合遥感识别技术不能“包治百病”,需要科学认识,因地制宜。要通过多层次多技术手段的综合运用,人防﹢技防,方能最大限度保障人民群众生命财产安全。

    “天”——在卫星层次,用InSAR和高分辨率光学遥感技术开展广域普查
    “空”——在航空平台层面,用无人机机载LiDAR和三维摄影测量技术开展局部详查
    “地”——地面专业人员有重点的进行调查复核,以及地面勘查和监测
    名词解释

    ◆InSAR:合成孔径雷达干涉测量,是指利用雷达卫星获取的相位信息,通过干涉测量的方式实现地表形态与形变测量的雷达遥感应用技术。特点是观测范围大,精度高,快速高效,可发现正在“活动”的隐患,特别是一些微小缓慢,人不易觉察的地表变形。

    ◆高精度DEM:即数字高程模型,指分辨率和精度在米级的三维地形数据。特点是精细刻画地表三维形态变化,与光学遥感、InSAR等联合使用,有助于识别灾害体的位置、范围、朝向、坡向等形态信息。

    ◆高分辨率光学遥感:即一般意义上的“遥感”,主要通过可见光谱段对地观测成像,是目前遥感应用最广泛的影像数据。高分辨率遥感影像可直观地反映地表覆盖、地物类型等要素,揭示与隐患相关的地表特征,如裂缝、鼓丘、滑塌等现象。

    ◆LiDAR:一般指三维激光扫描,通过激光雷达对地观测,获取精细三维地形数据。特点是可穿透植被,清晰看到斜坡结构、断层裂缝等地表信息。

    浩瀚宇宙,卫星凝视着地球,以“上帝视角”记录着地表的动态变化;山峦之上,无人机悄然掠过,精准“透视”出高山峡谷的旧伤新痕;边坡危岩,专业人员“按图索骥”,结合孕灾规律和地质条件靶向核查。

    上述天空地作业,均指向同一个目标——识别地质灾害隐患。

    小试牛刀

    全国地灾高中易发区迎来“体检”

    京城五月,繁花似锦,绿意盎然。

    位于北四环边上的奥运大厦6层,此刻却充斥着紧张的气氛:中国地质调查局自然资源航空物探遥感中心的技术团队争分夺秒解译着卫星遥感数据,赶在汛前,将全国121万平方千米地质灾害高中易发集中连片区疑似隐患“揪”出来,为各地防灾减灾提供参考。这是他们第二年开展“国家级”大规模地质灾害隐患识别。

    我国地质灾害点多面广,滑坡、崩塌、泥石流等致灾因素众多。一个无法回避的现实是,全国每年新发生的地灾70%都在已查明的隐患点范围之外。其中,滑坡最为常见,占全国灾害发生次数和影响程度的92%以上。

    消除“盲区”,解决“隐患在哪里”始终是“防”的难题。

    “2017年6月24日四川茂县新磨村高位远程滑坡和2018年金沙江白格滑坡的发生让业内清晰地认识到,仅靠传统的调查排查和群测群防手段很难对高位、高隐蔽性的滑坡隐患进行早期识别,必须借助现代对地观测技术,尽可能全面的从广域尺度上探测和发现隐患点。”中国地质调查局自然资源航空物探遥感中心遥感技术研究所负责人葛大庆说。

    为攻克防灾“卡脖子”环节,2019年,自然资源部部署开展了基于综合遥感技术的地质灾害隐患识别示范工作,目标两个:一是论证规模化开展隐患识别的技术可行性与应用效果;二是推动构建国家级和省级隐患识别中心,形成“1﹢N”的分级组织模式,由国家级隐患识别中心统筹全国尺度上的普查识别和技术引领,省级中心负责本省详细调查与现场核查,力求在3~5年内掌握全国地灾隐患的“底数”。其中,国家级隐患识别中心主要依托航遥中心推进隐患识别业务。

    经过充分论证与应用实践,确立了综合应用空天地多源遥感观测技术,以“形态、形变、形势”为识别内容的重大隐蔽性地灾隐患早期识别技术思路,形成了包括专题信息提取、隐患特征识别、野外核查验证的业务流程,丰富了地灾调查评价业务工作方式。

    葛大庆介绍,地灾隐患综合遥感识别技术就是通过高分辨率光学遥感影像与高精度DEM进行潜在灾害体“形态”调查,揭示地质背景、几何形态、地表覆被变化等成灾要素;以不同分辨率lnSAR监测获取斜坡体地表“形变”信息,表征活动状态,判别变形部位、滑移规模和活动阶段;以长时序InSAR形变数据、地质资料、威胁对象等综合分析变形活动的发展趋势以及致灾“形势”,判别疑似隐患,并最终经现场核查验证完成识别确认。

    “综合遥感识别从技术上体现了多源遥感手段的互补,降低了单一技术有效识别不足的概率。同时充分利用孕灾背景、基础地质条件以及已有隐患状况,从信息解译、判断、现场核查与最终确认上体现综合性,避免信息不足带来的认识偏颇。”葛大庆强调。

    2020年该项工作有了实质性进展。

    在自然资源部部署下,航遥中心组织相关力量对全国地灾高中易发区10省市开展地灾隐患综合遥感识别,覆盖黄河上游、四川强震区、藏东南地区、滇西北地区、三峡库区五类典型地灾分布区,面积118万平方千米。

    “此次识别聚焦大型活动性灾害隐患,以滑坡、不稳定斜坡为主,兼顾泥石流、崩塌等其他类型。同时,以人为本,突出对村庄、集镇等人口聚集地区和重大基础设施威胁严重的隐患。”葛大庆告诉记者。

    在识别层次上,一是发挥摸底的作用,对于“人迹罕至、人不能至”的地区,以及其他工作薄弱区或者空白区,综合遥感识别作为先行手段,为地面调查指明“靶区”;二是体现动态监测作用,通过技术手段揭示与灾害发育相关的“变化”特征,对于地面工作程度较高的地区,与大比例尺调查并行,互为补充。

    在工作模式上,实行分级协同。即由航遥中心牵头组织相关技术单位进行隐患识别,为省市提供疑似隐患识别成果。各地核查时首先将疑似隐患与已有隐患库对比,现场核查主要针对新发现的隐患,并根据危险性和危害程度开展相应的防灾工作。

    经过4个多月作业,全国地灾高中易发区221个县区8450处疑似隐患被“揪”了出来。经地方核查,识别正确率在30%以上。其中,青海、四川、湖北分别达到了76.5%、69.9%和62.86%,特别是四川省,核查确认的隐患近半为新增隐患,很好地体现了综合遥感技术发现未知隐患的作用。

    统计显示,2020年在罕见极端天气等多种因素影响下,全国地灾发生数比2019年增加26.8%,地灾造成的死亡(失踪)人数减少37.9%。这份守护人民群众生命和财产安全的答卷上,留下了综合遥感识别技术的印记。

    “地灾隐患综合遥感识别是一项非常有意义的探索性工作。”成都理工大学副校长、地质灾害防治与地质环境保护国家重点实验室常务副主任许强如是评价:

    其一,通过示范区域研究,基本构建了一套基于综合遥感技术的地灾隐患识别技术方法体系和流程。其二,从单一技术手段到多种手段的综合运用,形成了统一认识并对各类技术的优点和局限性有了比较清楚的认知。其三,地勘队伍在技术推广应用中,思想观念和工作方法有了明显改变。其四,识别了一批新的地灾隐患,尤其是滑坡,这是实实在在的成果。

    甘肃省地矿局第一地质矿产勘查院灾防院院长王岩对此深有体会。去年“8·13”陇东南特大群发地灾后,航遥中心下发了一批隐患分布图,灾防院对张家川县113个疑似隐患点进行了野外调查验证,确定地灾隐患80处,全部为滑坡,其中新增地灾隐患70处。

    “综合遥感识别技术能准确判断出斜坡是否变形,为技术人员判定滑坡灾害提供了可靠依据,同时加强了工作区地灾隐患点调查程度,补充完善了隐患点数据库。希望这类遥感核查项目能辐射更多的县区,相关技术能普及到地勘单位,提高技术人员对灾害的辨识率。”王岩说。

    广泛实践

    地方“自选动作”有声有色

    目前,地灾隐患综合遥感识别技术已被广泛接受,四川、贵州、陕西、重庆、甘肃等省市陆续开展识别任务,落实《地质灾害防治三年行动实施纲要》规定任务。技术力量雄厚的省份积极尝试“自选动作”,探索形成了具有地方特色的防灾减灾方式。

    “多灾多难”的四川,便是其中之一。

    四川省国土空间生态修复与地质灾害防治研究院副院长马志刚介绍,去年四川在全国开展示范工作的基础上,部署实施了省级地灾隐患综合遥感识别与监测项目,分片区对川北、川西南、川西、川东川南开展了InSAR监测和光学遥感筛查,发现和圈定了一批高位地灾隐患点;查明了重要城镇及交通沿线、地灾高易发区等重点地段的地灾情况。同时,选取了重大地灾隐患点或处于持续变形的重大工程治理点开展遥感追踪监测评价。

    记者了解到,和国家示范工作相比,四川省级地灾隐患综合遥感识别有4方面特点:

    一是分层级部署。在国家部署“面”的基础上,增设了“重点地段”和“重大地灾隐患点”两个层次,实现了综合遥感识别工作的粗细结合。二是多手段应用。为克服高密度植被覆盖区InSAR识别效果较差的技术短板,增加了“机载LiDAR”和无人机航空摄影,发挥其可穿透植被或高精度的优势。三是全灾种解译。将泥石流纳入隐患识别范围,实现了突发地灾滑坡、崩塌、泥石流等全灾种的解译分析。四是双路径分析。将找出存在形变且未登记在册的隐患点和查找正在变形的登记在册的隐患点作为识别的两个路径,一体化考虑,为地灾风险评价、汛前隐患排查提供了重要靶区。

    “目前四川已将地灾隐患识别纳入常态化防灾工作,并将解译成果作为地灾风险调查、汛前排查、重要强降雨天气后隐患排查的工作依据,实现了从被动应急救灾到主动防灾的转变。”马志刚说。

    “天无三日晴”的贵州,2016年在全国率先引进了综合遥感识别技术中的InSAR 开展地灾隐患的早期识别。5年来,全省共监测发现疑似滑坡形变区2000余处,经核查确认新发现地灾隐患600余处。

    “InSAR也抓住了一些地表变形不明显的在册地灾隐患点,帮助我们加强管理,采取综合防治措施,该治理的治理,该搬迁的搬迁。把‘隐患在哪里’问题又向前推了一小步。”贵州省地质灾害应急技术指导中心研究员刘秀伟告诉记者。

    近年来,贵州积极将技术逻辑与行政逻辑结合。去年9月,水城县发耳镇发生大滑坡,方量达80余万立方,InSAR和地面监测同时捕捉到地表动态,适时发布了预警信息,政府及时将受地灾隐患威胁的群众及重要财产转移到安全区域,把损失降到了最低。

    “事实上,2017年我们就监测到了发耳滑坡隐患。专家调查分析认为滑坡形变在可承受范围内,暂时不用大规模搬迁,通过持续监测和防控既能保证群众安全,又可保障当地生产经营活动有序开展,最大限度降低隐患对社会的负面影响。”刘秀伟介绍,3年来当地经济社会发展有序,老百姓生产生活如常,技术保障为行政决策提供了强有力地支撑。

    记者了解到,目前贵州地灾防治管理模式从以隐患点为主要管理对象,逐渐转变为以隐患点和风险斜坡为对象的“点面双控”风险管理。今年年底,全省以县为单元的精度优于1∶5万的地灾风险斜坡详细调查将全面完成。

    “斜坡与隐患点不同,‘病’不算重,常规设备监测投入大,管控效果不理想,但风险又不得不防。”刘秀伟说。对此,贵州正积极探索将InSAR运用到风险斜坡管控中,力争使该技术从地灾隐患识别迈向监测预警,提高风险斜坡的综合管控能力,为“灾害什么时候发生”探路。

    科学认识

    地灾隐患综合遥感识别技术不能“包治百病”

    多地实践证明,地灾隐患综合遥感识别成果可“指导排查,辅助调查,引导勘察”,提高了调查的针对性,减少了地面工作量,部分解决了“隐患在哪里”的难题。

    但专家同时强调,地灾隐患识别是多种技术手段、多学科知识的综合运用。综合遥感识别技术有其局限性,不可能“包治百病”,现有的技术手段仍难以识别前期没有明显变形或变化很快的突发性灾害隐患,这一点务必要科学理性地认识。

    “每种技术都有优势和短板,实操中须要因地制宜。”葛大庆说,比如,光学遥感容易受云雾天气和分辨率影响,大规模应用中往往难以获取亚米级的遥感数据,解译对专家经验和知识依赖度较高。

    “有的观点认为,综合遥感识别就是应用InSAR技术,这也是不对的。”葛大庆进一步阐释,InSAR能敏锐捕捉地表的动态变化,与光学遥感影像综合运用,可在植被覆盖度较低、高陡地区大范围普查,识别高位、隐蔽性的灾害隐患。比如在地层内部有损伤,地表持续活动的四川强震区,效果就很好。

    但是,InSAR会受到观测角度、植被状况、地形结构、水汽以及数据处理技术的影响和制约,普适性存在一定局限。同时,对于事前没有明显变形迹象,或者变形过程极短的突发性地灾,因其超出了雷达卫星有效观测的能力范围,InSAR很难发挥作用。比如2019年7月23日贵州水城县突发的特大山体滑坡,事前未出现明显变形迹象,InSAR未探测到变形信息。此外,在广东、浙江、福建等植被覆盖率较高的地区,光学遥感、InSAR等技术适用性受限,漏识率也较高。

    “天空地一体化防灾,以广域普查﹢局部详查﹢重点核查的方式分层次开展,是识别和防范地灾隐患最科学有效的方式。”许强认为。

    “天”指在卫星层次,用InSAR和高分辨率光学遥感技术开展广域普查;“空”指在航空平台层面,用无人机机载LiDAR和三维摄影测量技术开展局部详查;“地”即地面专业人员有重点的进行调查复核,以及地面勘查和监测。

    “通常,我们会用两种或多种技术手段背靠背地做,然后相互校验和补充。如果两种方法的结果能对上和相互检验,地灾隐患就可以确认了。”许强说。

    2015年,许强团队通过InSAR发现川西某村寨的山后有一处正在变形的大型滑坡区,由于最高海拔达到了4600米,人工没法上去调查复核,一直没能确认。3年后,在四川省支持下,团队运用LiDAR和三维摄影测量对该区域进行了航飞,经多种技术手段比对分析最终确认了蠢蠢欲动的滑坡隐患。

    多名业内人士指出,专业知识决定了技术应用的深度。在地灾防治领域,地质专家好比“医生”,InSAR、LiDAR等技术手段应用,如同为疑似隐患点做“CT”“B超”。是否患“病”,“病”征如何,最终还得依靠“医生”根据各种测试结果综合判定。如果不掌握足够的地质知识,对灾害形成条件、发育特征、成因机理等不了解,很难用好先进技术。

    扬长补短

    人防技防一起上,天上地下一起看

    甘肃立节北山老滑坡“苏醒”、四川洪雅突发高位远程崩滑、湖南安化切坡建房引发山体滑坡……今年以来,全国已发生200余起地灾,随着汛期的到来,地灾防治工作迎来大考。

    如何牢固树立以人民为中心的思想,进一步提升防灾减灾能力?——“多层次多技术手段综合运用,人防﹢技防”是专家们的一致观点。

    许强指出,近年来我国地灾防治工作取得显著进步,一是构建了一套具有中国特色的地灾群测群防体系,解决了我国地灾点多面广的问题;二是现代技术的有效运用,包括综合遥感技术、监测预警、计算机信息技术等,成效显著;三是引入了风险防控的理念,从以群测群防为代表的自发性防灾,发展到理论化防灾、技术化防灾、体系化防灾逐渐成熟。但是,地灾防治工作也面临瓶颈,最突出的就是对第三类隐患——具有突发性失稳破坏的潜在不稳定斜坡,目前还很难依靠技术手段提前发现。

    “地灾隐患可分为三类,第一类是正在变形的区域。第二类是历史上曾出现过变形破坏、有损伤的区域。第三类就是潜在不稳定斜坡,即历史上从未出现变形破坏,目前也无明显变形迹象,但其稳定性相对较差,暂时处于基本稳定或欠稳定状态,一旦外部条件发生变化,比如一场强降雨或人工开挖坡角就有可能突发地灾。”许强阐释。

    记者了解到,对于第一类隐患,InSAR和高分辨率光学影像能很好地进行识别。第二类隐患,运用LiDAR 最为有效,并可以三维摄影测量辅助。而第三类隐患,靠遥感技术已无能为力,必须动用物探、钻探等手段,把斜坡地下结构和含水状况查清楚才能判断其稳定性,但潜在不稳定斜坡点多面广,防不胜防。

    “现阶段,第三类隐患还主要依靠群测群防来防范,当然我们也正在研发航空物探技术,希望能快速查明地下结构,评价其稳定性。”许强说。

    今后一段时期,地灾防治可以从哪些方面发展突破?

    许强建议,一是要构建航空遥感联盟,搭建信息共享平台,整合全国航空遥感力量,尽快实现航空遥感数据的全国全覆盖。二是要发展航空物探,尽可能多地将潜在不稳定斜坡的地下结构和含水状态搞清楚,把第三类隐患找出来。

    葛大庆认为,提高综合遥感识别能力,找到更多地灾隐患,一方面要着力发展专业目标的遥感卫星,强化源头的有效观测与数据保障,进一步提升卫星监测的频率、精度和针对性。“地灾隐患识别目前主要用欧空局哨兵-1号卫星C波段数据,以及日本ALOS-2卫星L波段数据,但受卫星观测周期和范围限制,一些地方还不能有效覆盖,分辨率也不够,限制了可识别的范围、规模和精度。”葛大庆说。

    另一方面,大规模业务化识别工作要进一步提高遥感识别信息处理的自动化和智能化水平,尽可能将专家知识转为模型或者方法,提升智能识别“去伪识真”的本领,降低人员劳动强度。此外,还要加快培养具有地灾专业知识,又掌握遥感、InSAR以及地理信息技术的人才。

    马志刚指出,现阶段地灾源头管控理念还有待提升,“‘源头管控’虽然是旧词新用,但在地灾防治上是一个全新的认识,是推动地灾‘早期防’到‘源头控’观念转变的重要变革。”

    近年来,四川按照以防为主、防治结合的方式开展了地灾易发县县域地灾风险评价、隐患识别分析等工作,初步摸清了隐患和风险的家底,获得了较为丰富的基础性成果资料。但是,扎实推进基础工作并不代表地灾“防”的工作就做好了,如何依托现有成果和技术把地灾隐患“管起来”,最大限度地降低灾害风险是当前的难题。

    对此,马志刚提出两点建议,一是加快构建地灾隐患点、风险区双控新格局。充分运用地灾易发县县域地灾风险评价、遥感识别监测等成果,选择重点区域开展风险管控试点,形成“识别、监测、管控”结合的风险管控体系,不断提高地灾隐患识别和风险管控能力。

    二是努力做好科技防灾。在加大群专结合监测预警体系建设的同时,重点建立一个综合监测预警平台,构建一套地灾监测预警管理体系,培养一支技术团队,形成“风险提前管控,空天地隐患识别,气象风险动态预报、专业设备持续监测、信息平台实时预警,现场人员动态核查”的监测预警新局面。

    (绘图 陈舒)

    地灾隐患识别:寻找大地上的“潜伏者”

    ▶中国地质调查局岩溶地质研究所实施的《乌蒙山水文地质调查》项目,助力四川省昭觉县利用温泉发展生态渔业。 邹胜章 供图

    ▶汾河的源头—雷鸣寺泉。申豪勇 供图

    ▶江西省宁都县小源村安全饮水示范井,可保障全村3500余人饮水。 邵长生 供图

    3月22日是第28届“世界水日”,3月22日~28日是第33届“中国水周”。我国纪念2020年“世界水日”和“中国水周”活动的宣传主题为“坚持节水优先,建设幸福河湖”。

    水是人类赖以生存和发展的珍贵资源。地球上虽然“三分陆地七分水”,储水量有14.5亿立方千米之多,但其中淡水储量仅占2.5%,而人类可以直接利用的湖泊、河流和浅层地下的淡水资源仅占地球总水量的0.26%左右。随着人口的增长和经济的发展,人类对水的需求不断增长,而与之对应的是,在工业化、气候变化等影响下,淡水资源日渐短缺,许多国家陷入缺水困境。

    我国人多水少、水资源时空分布不均,水安全问题事关经济社会发展稳定和人民健康福祉。近年来,自然资源部中国地质调查局坚持山水林田湖草是一个生命共同体,从地下水调查转向地表水与地下水一体化调查,从注重水的资源功能转向资源、环境、生态功能并重,从水资源阶段性评价转向周期性评价,并着力开展流域水文地质与水资源调查,在服务生态文明建设、精准扶贫和区域经济社会发展中发挥了重要作用,为保护和合理开发利用水资源提供了支撑与服务。

    今天,我们邀请水文地质调查领域的专家和有关项目负责人,谈谈如何合理开发利用水资源。

    黄河流域: 量水而行须摸清水资源家底

    韩双宝

    黄河是中华民族的母亲河。在自然环境变化与人类活动的双重影响下,黄河流域面临严重的水资源与水平衡问题。流域水资源总量贫乏,开发利用程度高;水资源时空分布不均,产水区、用水区不匹配,有效利用难度大;耗水量持续增长,供需矛盾突出;生态环境需水量研究不深、配置不足,造成湿地萎缩退化等一系列水资源与水平衡问题。

    党的十八大以来,习近平总书记多次实地考察黄河流域生态保护和发展并作出重要指示,明确黄河流域生态保护和高质量发展是重大国家战略,发展要坚持生态优先,把水资源作为约束,以水而定,量水而行,合理规划黄河流域人口、城市和产业发展。

    为摸清流域水文地质条件,以及水资源数量、质量、开发利用及其影响,在自然资源部中国地质调查局的统一部署下,中国地质调查局水文地质环境地质调查中心联合地调局有关单位,开展了黄河流域水文地质调查工程。

    目前,该工程已完成流域内34万平方千米地下水主要分布区的水位统测工作,掌握了黄河流域地下水流场及水位变化特征。在湟水河流域、渭河流域、汾河流域、黄河下游干流区等典型流域开展水文地质与水资源调查,实施了一系列地表水—地下水转化断面监测、入渗与蒸发试验,为流域水资源评价提供关键水文地质参数。此外,还开展了流域水平衡分析研究,明确了目前黄河流域上中下游重点地区的水平衡问题。

    为支撑黄河流域生态保护与高质量发展,水文地质环境地质调查中心将继续组织开展水文地质与水资源调查、地下水统测、地下水资源评价、信息化与数据库建设、综合研究等工作,以干支流河谷区、银川平原、河套平原、关中平原、鄂尔多斯盆地、太原盆地、运城盆地、临汾盆地和黄河下游干流区为重点,开展年度统测,掌握地下水位状况、盈亏量及其相关的生态环境问题,分析区域地下水位变化趋势、成因及应对对策;科学划分水文地质单元,开展地下水资源数量、质量评价,分析生态环境问题,提交国情数据;针对典型子流域重点生态环境问题,查明水文地质条件与水资源状况,揭示水资源变化与生态环境的关系,提出地下水资源合理利用、生态保护建议。

    围绕黄河上游局部地区生态系统退化,中上游主要平原(盆地)地下水开发利用程度高、生态环境变化大,下游悬河段地表水—地下水转化关系变化、三角洲湿地生态变化等问题,项目将开展重点调查和水平衡分析,服务国土空间规划与生态修复。

    (作者系黄河流域水文地质调查工程首席科学家)

    汾河流域: 保护岩溶大泉迫在眉睫

    申豪勇 谢 浩

    汾河发源于山西省宁武县东寨镇管涔山脚下的雷鸣寺泉,在山西省万荣县荣河镇汇入黄河,全长716千米,控制流域面积39741平方千米,是黄河的第二大支流。汾河横贯整个山西省中南部,流经包括省会太原市在内的6个地市,流域面积约占山西省总面积的1/4。

    汾河流域以发育有岩溶大泉为主要特征,受地形地貌和地质构造的控制,流域内共发育了8个规模不等、流量不等的岩溶大泉,从上游到下游依次为雷鸣寺泉、晋祠泉、兰村泉、洪山泉、郭庄泉、霍泉、龙子祠泉和古堆泉。岩溶大泉的泉水流量稳定、水质良好,且便于开发利用,是周围城市生产和生活的主要供水水源。但随着经济的快速发展,流域内水资源供需矛盾日益突出,水资源过度开发引发的生态环境问题日益严重,加强岩溶大泉的保护迫在眉睫。

    随着泉域内自然条件的改变和人类活动强度的加剧,岩溶泉域的地质结构、水循环要素以及输入输出条件都发生了根本性变化,并诱发了一系列水资源环境问题。首先,岩溶大泉的泉水流量大幅度衰减甚至断流:20世纪60年代汾河流域8个岩溶大泉的总流量达每秒27.09 立方米,到2017年已衰减至每秒8.56立方米,其中兰村泉、晋祠泉、古堆泉和洪山泉已经断流。其次,区域岩溶地下水位持续下降:20世纪80年代到2010年,年平均下降速度为1米~2米;2010年以后,由于当地政府实行地下水压采政策,下降速度有所减缓。再次,岩溶地下水的水质呈恶化趋势:受地质结构影响,汾河流域多数岩溶水系统为“水煤共存”系统。煤炭资源开发,特别是近年来闭坑煤矿形成的“老窑水”对岩溶地下水资源的污染日益严重。

    中国地质调查局岩溶地质研究所组织实施的《汾河流域晋中南大型岩溶泉域1:5万水文地质调查》项目,对汾河流域中游的晋祠泉域和兰村泉域开展了1:5万水文地质调查,掌握了泉域内岩溶地下水资源的开发利用现状,并提出了晋祠泉复流措施。利用项目建立的地下水渗流数值模型,对各种措施的效果进行定量评价预测,晋祠泉有望在未来2~3年恢复出流。

    针对汾河流域的水资源供需矛盾问题及生态环境问题,加强汾河流域地下水资源保护、开展岩溶大泉的生态环境修复需要重点关注以下几个方面的工作:一、加强岩溶大泉基本水文地质条件的调查工作,尤其是掌握其随自然或人类活动产生的变化情况,这是开展地下水资源合理开发利用与保护的前提;二、查明各岩溶大泉泉水流量衰减或断流的成因,为泉水的生态修复提供基础依据;三、开展“老窑水”对岩溶地下水的环境效应分析,揭示岩溶大泉水质恶化的原因。

    (作者分别系中国地质调查局岩溶地质研究所《汾河流域水文地质调查》项目负责人、成员)

    华北平原: 加强地下水合理开发是关键

    石建省 张翼龙

    地下水长期过度开采会引发一系列资源、生态、地质环境问题。目前,我国华北平原已经成为世界上地下水降落漏斗面积最大、地面沉降最严重的地区。2014年,国家启动河北省地下水超采综合治理试点;2019年,国务院印发实施《华北地区地下水超采综合治理行动方案》。通过城镇自备井关停、最大限度利用外调水、置换农业农村生活水源、探索旱作雨养种植模式、开展休耕轮作试点、大力发展节水灌溉等一系列措施,目前,华北地区用水总量增长得到一定程度遏制,城市集中供水区地下水位下降态势初步扭转,但区域地下水降落漏斗发展态势依然严峻。

    地下水超采治理是一个长期过程。现阶段地下水超采治理以控采为主,从长远来看,超采治理不能仅是禁采和限采地下水,地下水位的过度回升也会导致一系列的生态环境问题。例如,浅层地下水水位过高会导致地下水蒸发增加,形成土壤盐渍化;过高的地下水水位可能会造成地下空间和建筑地基失效。因此,从水平衡角度出发,加强地下水合理开发是地下水超采治理需要研究解决的重要问题。

    自然资源部中国地质调查局部署开展的“海河—淮河流域水文地质调查工程”对华北平原水平衡变化进行了重点研究。基于此,就科学开展地下水超采治理提出以下建议:

    一是注重水平衡变化分析。加强超采治理的水文变化过程研究,进行水平衡动态分析评价,根据影响水平衡的关键要素,优化调整治理措施,制订科学合理的地下水开发方案,充分发挥地下水资源的作用。

    二是科学确定超采治理目标。针对不同区域地下水功能,根据稳定地下水位、防止水资源枯竭、减缓地面沉降、保证城市地下空间工程安全、防止土壤盐渍化、保护湿地资源等综合指标,分层、分区划定科学化、精准化超采治理“红线”。

    三是遵从和维护水循环自然规律。建议尽量恢复水循环自然功能,水资源利用要因“地”制宜,因“质”制宜,宜地表则地表,宜地下则地下。要保证优质地下水资源首先用于人民饮用,下决心全面置换工业用水的地下水开采,根据水源置换能力和地下水可供性综合调整农业产业结构,切实做到“以水定产”。

    四是推进非常规水源开发利用。华北平原中东部地区的中浅部分布有大范围微咸水、半咸水,资源量约36亿立方米/年,尚未得以有效开发利用。应进一步创新工作思路,转变研究理念,加强咸水、微咸水改善与开发利用技术研究,最大限度地实现咸水、微咸水资源化,以咸水、微咸水的区域性改善与规模化开发,促进华北平原地下水可持续利用。

    五是加大地下水安全回补。充分利用地下水超采形成的地下库容,发挥地下水的调蓄功能,开展多水源联合调度,实施地下水安全回补与调蓄工程,构建集地下水回补及调蓄—风险评价—监控预警—含水层管理于一体的地下水涵养体系,以实现地下水安全回补以及含水层综合调蓄与管理,进而提高区域水资源及水环境承载能力。

    六是优化地下水监测体系。现有的国家和省级地下水监测点以区域性控制为主,针对不同水文地质单元布设,尚不能有效支撑地下水超采治理目标划定和效果评估工作。建议针对地下水严重超采区、地下水降落漏斗分布区、地面沉降严重区等区域,优化和加密地下水监测网,建立精准地下水超采治理监测网,为全国地下水超采治理监测起到示范作用。

    (作者分别系中国地质调查局水文地质环境地质研究所党委书记、副所长,水文地质与水资源调查计划首席科学家;中国地质调查局水文地质环境地质研究所科学技术处处长)

    赣州四县: 让贫困区缺水百姓喝上干净水

    邵长生

    “终于能喝上干净水了!”今年春节前夕,江西省宁都县青塘中学收到了一份特殊的新年礼物——由中国地质调查局武汉地质调查中心建设的安全饮水示范井。该井是武汉地调中心落实自然资源部和中国地质调查局党组工作部署,在江西宁都县青塘中学组织实施的一项惠学惠民工程,解决了困扰该校20余年的安全饮水难题。

    让缺水地区更多的群众喝上干净、卫生的地下水,一直是中国地质调查局支撑脱贫攻坚的目标之一。江西省赣州市兴国县、赣县区、于都县、宁都县(以下简称“赣州四县”)是自然资源部定点帮扶县。这里不仅季节性缺水问题突出,而且由于地处红层区、浅变质岩区和花岗岩区,地下水基本水文条件不佳,找水打井困难。自然资源部加强对定点帮扶的赣州四县支持力度,要求中国地质调查局发挥专业优势,做好“两不愁三保障”中的饮水安全保障工作。

    针对赣州市贫困地区农村饮水困难及旱情导致的季节性严重缺水等问题,2019年,地调局组织武汉地质调查中心、南京地质调查中心和水文地质环境地质调查中心3家直属单位,在赣州四县开展找水打井专项水文地质工作,共实施了183口水井,涉及106个村组、10个千吨万人水厂,直接解决了7.74万人的安全饮水问题,并为12.1万群众提供补充水源。

    2019年,赣州地区遭遇近百年一遇干旱灾害,持续时间长达6个月,河塘干涸、溪水断流,群众饮用水困难问题越发突出,部分水厂限时供水,甚至无水可供。随着脱贫摘帽时间接近,该类村组的安全饮水问题更成为当地政府部门“两不愁三保障”工作的头等大事。中国地质调查局急群众之所急,及时开展找水打井工作,提供应急备用或补充水源,让当地政府和群众吃下定心丸。如赣县杨西村水厂地表水源干枯,武汉地调中心及时为该村建设了一处日供水量可达256吨饮水示范井,大埠乡正规划将该处作为乡政府驻地近1万人的补充水源。

    下一步,自然资源部中国地质调查局及直属单位将继续落实党中央和自然资源部党组关于脱贫攻坚的指示精神,坚持脱贫不脱政策、脱贫不脱帮扶、脱贫不脱责任,继续充分发挥专业优势进行帮扶。一是查漏补缺,在安全饮水保障程度不足、高铁锰氟等水源质量差的问题地区,继续开展找水打井和改善治理示范;二是聚焦产业扶贫,针对赣州四县特色的天然矿泉水资源、地热资源和富硒蔬菜产业灌溉水资源等需求,开展勘查与评估,为当地特色水产业可持续发展提供服务;三是提出赣州四县地下水资源开发利用建议,指导当地地下水资源开发利用,为实现脱贫攻坚和乡村振兴有效衔接提供地质支撑。

    (作者系中国地质调查局武汉地质调查中心《赣南地区饮水安全水文地质调查》项目负责人)

    幸福之水哪里来?