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

查询时间:时间不限

  • 时间范围


查询范围:标题+正文

排序方式:按相关度排序

检索范围:

    近期,中国地质科学院岩溶地质研究所在西南岩溶石漠化区生态恢复领域取得新进展。

    岩溶石漠化生态系统植被恢复普遍面临土壤氮素匮乏的制约。豆科植物通过根瘤菌共生体系实现生物固氮作用,将惰性大气氮转化为可利用的铵态氮,因此在岩溶区生态恢复演替初期,豆科植物常作为先锋物种占据优势生态位。然而,不同豆科物种在促进土壤氮有效性方面存在显著种间异质性。除共生固氮外,豆科植物可能通过根际调控作用影响土壤氮循环关键过程,进而影响土壤氮有效性。这种植物-土壤互作对氮素有效性的调控路径及其影响机制仍未明晰,制约着石漠化治理中豆科植物的科学配置。

    岩溶所科研人员通过系统研究云南岩溶断陷盆地石漠化区豆科植物和非豆科植物生长的氮限制程度和土壤氮转化过程,发现豆科植物生长的氮限制程度显著低于非豆科植物。豆科植物经过有效改善土壤环境,可提高土壤无机氮供应能力,进而增加氮有效性,能够缓解植物生长氮限制程度。当然,这种能力因物种类型不同有较大差异,土壤无机氮供应能力越大,豆科植物氮限制程度越低。因此,豆科植物在土壤氮循环和有效性中发挥着重要作用,对于石漠化区生态恢复初期植物群落物种的配置具有重要理论指导意义。

    岩溶石漠化区豆科植物和非豆科植物氮限制及其影响机制概念图

    岩溶所在西南岩溶石漠化区生态恢复领域取得新进展

    近日,中国地质调查局海口海洋地质调查中心(以下简称“海口中心”)海岸带地质调查团队在粤东近岸海域表层沉积物碳、氮、磷的分布规律及污染评价研究方面取得新进展。

    粤东是广东省经济发展的重要动力引擎,其近海海鲜养殖和旅游业尤受欢迎。但随着海岸带经济的快速发展,各滨海城市均面临着巨大的生态环境保护与治理的压力。为了解粤东河口及近岸海域表层沉积物中总有机碳(TOC)、总氮(TN)、总磷(TP)的分布特征及污染状况,海口中心海岸带地质调查团队以粤东近岸海域为研究对象,全面分析表层沉积物中TOC、TN、TP的含量分布情况,并对其环境污染状况进行评价。此外,该研究还对沉积物中有机质的来源进行了探讨,并对比研究了该区域近15年TOC、TN、TP含量及来源的变化。

    研究结果显示:义丰溪河口、榕江河口、濠江河口、柘林湾和南澳岛养殖区为TOC、TN、TP主要富集区。TOC、TP含量分布具有空间异质性,表现为湾内及河口区域含量显著高于湾外及离岸海域。TN含量分布则呈现离岸海域高于近岸趋势。研究区表层沉积物中有机质主要以海源输入为主,要素之间具有显著正相关性,同时要素与粒径间相关性较好。研究区表层沉积物大部分为轻度污染状态,TOC、TP受危害等级为Ⅰ类,基本未构成污染,TN为Ⅱ类污染,处于受污染状况。本研究对了解粤东近岸海域表层沉积物现状与污染生态风险具有重要意义,可为当地近海海洋环境保护和治理提供科学依据,同时为海洋碳汇研究提供基础数据支撑。

    该研究由海口中心海岸带地质调查团队林聪、王洋等合作完成,相关研究成果以“粤东近岸海域表层沉积物碳、氮、磷(TOC-TN-TP)分布规律及环境评价分析”为题发表于《水生态学杂志》期刊,得到中国地质调查项目(潮汕海岸带综合地质调查:DD20208013,惠州-汕尾海岸带自然资源综合调查DD20230415)的资助。

     

     

     
     
    海口中心在粤东近岸海域表层沉积物分布和污染评价方...

    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.

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

    近日,自然资源部、广西壮族自治区岩溶动力学重点实验室在岩溶区水体水气界面碳循环研究方面取得新进展。 

    岩溶区水体水气界面碳循环是受到广泛关注的全球内陆水体水气界面碳循环的重要组成部分,也是岩溶作用碳循环研究的重要内容,是定量准确评价流域岩溶碳汇效应的关键。但由于水体水气界面碳交换存在高度的时空异质性,使得准确量化水体水气界面碳通量存在着极大的不确定性和面临巨大的挑战,也给岩溶碳汇的准确评估带来极大的困难。针对上述科学问题,在国家自然科学基金,地质调查项目,广西自然科学基金,基本科研业务费项目的支持下,岩溶所蒲俊兵研究员团队利用自主设计的浮游静态箱和薄边界模型法两种方法,分别选择岩溶地下水补给的溪流、水库和流经岩溶区的大型河流作为研究区,进行不同时空尺度的野外监测研究,获得了新的科学认识。

    对流经岩溶区的大型河流(桂江)的25个断面的监测研究表明,桂江水气界面二氧化碳交换表现出显著的季节和空间变化。桂江流域全年总体上为大气二氧化碳的源,但是在冬季,由于水生生物代谢过程和碳酸盐缓冲体系的影响,二氧化碳脱气受到极大的限制,甚至表现为大气二氧化碳的汇。由于在冬季二氧化碳脱气受到极大的限制,使得二氧化碳脱气占到流域无机碳通量的20.1%,形成了每年12.8万吨碳的碳汇通量。对以岩溶地下河为主要补给的(补给量占90%以上)的大龙洞水库的研究表明,大龙洞水库水气界面二氧化碳交换量表现出显著的季节变化和昼夜差异。大龙洞水库类似湖泊的“热分层效应”控制了二氧化碳脱气量的季节变化,在“热分层效应”显著的夏季,二氧化碳交换量受到抑制且出现显著的吸收大气二氧化碳的现象。相反,在冬季水库发生混合时,水库水气界面发生强烈的二氧化碳的脱气,其脱气量最高。此外,在“热分层效应”显著的夏季,二氧化碳脱气量昼夜变化差异较小,但是在翻库期间,二氧化碳脱气量昼夜变化差异较大。对以岩溶地下河为主要补给的(补给量占98%以上)的官村地表溪流的4个断面的研究表明,官村溪流水气界面二氧化碳交换表现出强烈的空间差异性,在1.3公里的距离内,二氧化碳交换量从源头的423.4毫克每平方米每小时急剧下降到最下游点的257.0毫克每平方米每小时,说明在岩溶地下河-溪流的体系中,当岩溶地下水转移到地表后会发生强烈且迅速的碳排放。官村溪流每年向大气释放21.75吨二氧化碳,但仅占流域无机碳损失通量的1.71-5.62%。同时,在官村溪流中降雨期的高分辨率监测研究表明,降雨期间受到水生光合生物新陈代谢过程的影响,溪流仍旧存在溶解无机碳的昼夜变化现象。计算表明,在降雨事件期间,溪流内沉水植物每天通过新陈代谢过程会形成5.6公斤的碳汇量。研究结果对于岩溶区溪流中沉水植物可以稳定碳酸盐岩风化作用的碳具有重要意义。

    上述研究成果,强调了水生光合作用和碳酸盐缓冲体系在抑制岩溶水体中二氧化碳脱气方面的重要性,定量评估了具有高度时空异质性的水体二氧化碳交换通量,强调了要准确评估岩溶水体水气界面二氧化碳交换通量和碳汇量,必须根据其高度的时空异质性特征,改进监测方法,提高通量评估的准确性。该系列研究成果于近期发表在《Journal of Hydrology》(国际TOP期刊)、《Environmental Science and Pollution Research》、《Environmental Earth Sciences》等国际期刊上。

    水生光合生物新陈代谢过程和碳酸盐缓冲体系对岩溶河流水气界面碳循环过程的影响示意图

    岩溶区水体水气界面碳循环研究取得新进展

      

    石漠化地区

      

    石漠化地区三七种植基地

     

    近日,由中国地质调查局岩溶地质研究所牵头的国家重点研发计划项目《喀斯特断陷盆地石漠化演变及综合治理技术与示范》正式启动。本项目将以国家石漠化治理工程区的蒙自、建水、泸西三个盆地为研究区,通过研发和集成喀斯特断陷盆地石漠化演变机理、水土漏失阻控、植被恢复、特色资源开发及生态衍生产业构建,建立喀斯特断陷盆地生态环境基础数据库、技术模式与方法库及信息共享平台等,形成断陷盆地生态安全评估技术体系,建立喀斯特断陷盆地综合治理理论与技术体系。

     

    1

    中国在喀斯特生态及石漠化治理研究中居领先地位,石漠化治理仍需加强

     

    项目首席科学家、岩溶所曹建华研究员介绍,喀斯特生态环境的脆弱性导致石漠化问题在世界不同喀斯特地区普遍存在。石漠化在中国西南、中南半岛、中东地区、地中海沿岸均有分布。除中国西南、中南半岛外,其他国家更多关注的是石漠化区的资源禀赋及其导致的干旱缺水问题,而对生态退化及修复关注较少。

    中国在喀斯特生态及石漠化治理研究中居领先地位。我国从“十五”时期就开始了石漠化的治理工作,中国科学家连续组织实施5个IGCP项目,针对石漠化治理进行研究。在国家石漠化综合治理工程的努力下,我国石漠化面积由逐年增加变为逐年下降,总体趋势好转,西南地区石漠化土地总面积减少7.4%,植被盖度增加4.4%;生态产业得到大力培植与发展,实现200万人脱贫;岩溶地下水开发工程解决了部分人畜饮用水问题。

    2015年4月,中共中央国务院发布《关于加快推进生态文明建设的意见》,明确提出继续推进石漠化综合治理。2016年3月,国务院发布的“十三五”规划纲要明确指出,荒漠化、石漠化、水土流失综合治理是推进国家重点区域生态修复的主要内容,标志石漠化综合治理工作又进入了新的阶段。

    目前,我国石漠化治理在以下几方面仍需加强:石漠化演变机理研究需加强,断陷盆地类型区尤其薄弱;缺水问题没有得到根本解决,水资源利用效率有待提高;水土漏失的阻控技术尚未形成;石漠化区植被退化修复生态学机制有待阐明;特色生态产业物种配置,引进物种适应性及功能缺乏综合效益评价;断陷盆地石漠化治理模式和技术集成欠缺。

     

    2

    我国断陷盆地石漠化治理起步晚,治理模式和技术相对单一,亟须实现理论突破和技术创新

     

    尽管我国石漠化治理取得一定成效,然而,目前已开展的系列国家“973”计划、科技支撑计划和重点自然科学基金等国家重大项目主要集中在喀斯特高原和峰丛洼地区,断陷盆地鲜有涉及。

    断陷盆地区是我国石漠化最为典型、程度最严重的地区之一,研究显示,我国断陷盆地石漠化面积比例达1.51万平方千米,占喀斯特面积高达32%,集中分布在我国的滇东和攀西一带。由于断陷盆地盆—山地形变化剧烈,气候反差大,降水偏少,蒸发量大,地下水埋藏深、水土资源分离、植被立地条件差等一系列特点,造成断陷盆地石漠化治理起步晚、投入少,导致一系列资源环境问题。

    一是水土漏失加剧。调查结果显示,喀斯特地区土壤侵蚀以地下漏失为主,地下漏失模数占年均总土壤侵蚀的70%~80%。蒙自南洞地下河是水土漏失过程的典型代表,该地下河出口泥沙含量达0.77kg/m3,是南盘江的2.33倍。整个南洞流域水土流失总量为43.6万t/a,高原面土壤垂向漏失是造成地下河泥沙含量提高的主要原因。断陷盆地水土漏失严重威胁长江、珠江中下游经济区的生态安全。

    二是水资源开发与利用不匹配,缺水困难没有根本解决。据调查,西南喀斯特区仍有6637万人存在饮水不安全问题,缺水人口达1947万人。目前完成的西南喀斯特区1:5万水文地质调查面积23万平方千米,可开发利用水资源量615亿立方米/年,但开采程度仅16%。断陷盆地地下水深埋,受制于技术手段落后,地下水勘探、开采、利用难度大,缺水问题仍十分严重。

    断陷盆地地区隶属国家级连片特困地区——乌蒙山区和滇黔桂石漠化区,其中分布在滇东的国家级贫困县有24个,22个是喀斯特县。生态环境恶劣的石漠化区人口密度高,贫困人口近百万,贫困面大且程度深,是西南乃至全国扶贫攻坚的重点地区。

    最新遥感数据显示,滇东—攀西断陷盆地重度石漠化面积呈增加的趋势。滇东—攀西地区隶属国家“两屏三带”生态安全屏障的黄土高原—川滇生态屏障,国家14个连片特困地区的乌蒙山区、滇黔桂石漠化区,影响着长江、珠江、澜沧江等国内、国际重要河流的生态安全,是我国重要的生态保护功能区。

    相比其他喀斯特地貌类型区,我国断陷盆地石漠化治理在基础理论研究和治理关键技术研发方面相对薄弱,治理模式和技术单一,效果较差,存在着“地头水柜水不见,山高水深饮水难,水土流漏无措施,地下“红”河愁不清,山路崎岖生活艰,居民贫根难拔除”的现象,亟须实现理论突破和技术创新,形成科技支撑服务的、产学研用为一体的生态治理—生态产业协同发展的断陷盆地石漠化综合治理模式。

     

    3

    选择蒙自、建水和泸西3个典型研究区,拟解决3个重大科学问题

     

    据曹建华介绍,《喀斯特断陷盆地石漠化演变及综合治理技术与示范》重点研发计划项目将由岩溶所牵头、联合中科院亚热带农业生态研究所、中国水利水电科学研究院、北京林业大学、云南省地质调查局等单位产学研用攻关,针对断陷盆地盆—山共存的环境地质结构、水土资源不匹配和石漠化严重等问题,选择断陷盆地典型、石漠化程度严重的蒙自、建水和泸西3个研究区开展研究工作。

    三个研究区各有特点。其中,蒙自研究区为断陷盆地喀斯特发育强烈、水土流/漏失严重的典型代表,拟在该地区进行石漠化演变机理研究、水土流/漏失阻控技术和生态产业研发和示范。建水研究区气候干热,缺水严重,具有多级盆地结构特点,拟在该区进行地表—地下水联合调控、植被恢复及功能提升研究。泸西示范区,以小江流域断陷盆地为核心区,流域上游区为喀斯特槽谷与喀斯特峰丛分布区,中游为断陷盆地,下游为峡谷区,将重点研究示范区断陷盆地生态环境地质分异特征及其石漠化演变机理,对断陷盆地“五水”循环过程及水资源均衡进行分析,研发地表—地下水联合调控及水资源高效利用及优化配置技术,构建立体生态产业模式,集成泸西断陷盆地水土资源高效利用、石漠化生态治理—生态富民耦合调控技术,进行综合效益、适用性评价,在县域(流域)范围进行应用示范。

    项目将重点解决3个重大科学问题:

    一是地质—气候制约的断陷盆地生态系统运行规律及石漠化驱动机制。

    针对断陷盆地急剧变化的地形、气候特点以及“盆—山”共存的环境地质结构,加之断陷盆地生态环境地质具有高度异质性,形成特殊的生态水文及其耦联的碳、氮、钙物质循环过程,影响和控制了区域喀斯特生态系统运行和石漠化的发生、演变过程,重点揭示地质—气候制约的断陷盆地生态系统运行规律及石漠化驱动机制,有效服务于断陷盆地水土资源高效利用,进行土壤侵蚀阻控,退化生态修复,生态产业模式构建。

    二是断陷盆地地表—地下水、土过程与资源高效利用调控机理。

    断陷盆地特殊的地质结构和水文地质条件,导致地下水埋深大,“土在楼上,水在楼下”,空间分布极不均匀又不配套,水资源以地下水为主,土壤资源集中分布在盆地平坝区。将对断陷盆地不同喀斯特形态影响水土过程、制约了水资源和土壤资源的开发利用机理等进行研究。

    三是生态服务功能提升与生态产业培育协同机制。

    基于喀斯特断陷盆地富钙偏碱的地球化学背景,导致喀斯特区植被具有喜钙性、旱生性、石生性。缺水少土导致立地条件差,生态恢复的难度大。生态服务功能的提升需要系统阐明植被退化与恢复生态学机制;生态产业涉及的物种与断陷盆地区域地质生态环境协同的生态学原理需进一步揭示;拟建立断陷盆地石漠化治理评价指标体系与评价方法。

     

    4

    设置6个课题,系统研究石漠化演变机理及综合治理技术与示范

     

    根据项目的重大科学问题、构建技术、生态产业模式等目标,本重点研发项目设置6个课题,拟系统研究石漠化演变机理及综合治理技术与示范。

    《断陷盆地生态环境地质分异及石漠化演变机理》课题,将选择代表性的喀斯特断陷盆地石漠化区,定量刻画典型流域植被与水文过程交互作用,确立流域生态需水关键期及需水量,明确植被生态水文耦合过程对碳、氮、钙、水等物质传输的影响。在此基础上,解译喀斯特断陷盆地生态环境地质分异规律及其内在联系,揭示喀斯特断陷盆地生态环境地质分异与石漠化协同演变机制、生态系统演替趋势,为喀斯特断陷盆地石漠化区面向生态的水资源合理配置和生态功能恢复提供理论依据。

    《断陷盆地地表、地下水资源高效利用与优化调控》课题,针对断陷盆地水资源以地下水为主、且埋深大、雨水集中、水资源时空分异大、水资源与土壤资源空间上分隔、时间上存在错位等突出问题,以泸西、建水、蒙自三个盆地为解剖对象,开展断陷盆地地表—地下水过程研究、资源高效利用调控技术研发,阐明断陷盆地流域尺度水文地质结构与水文过程,分析“五水(雨水、地表水、喀斯特表层带水、地下水、土壤水)”转换过程,研究喀斯特系统水均衡问题,研发地下水开发、地表水调蓄、雨水收储、水污染风险评估及水质修复、节水灌溉、土壤保墒等综合高效利用技术,探索断陷盆地水资源可持续利用管理技术途径提供科技支撑。

    《断陷盆地土壤流/漏失阻控与质量提升》课题,选择蒙自、建水断陷盆地,针对喀斯特断陷盆地周边山地(陡坡带、高原面)土层浅薄不连续、土壤流/漏失严重、土地持续利用难等问题,查明陡坡带、高原面岩土组构与土壤生态特点,揭示土壤流/漏失驱动因素、土壤退化对石漠化演变的响应关系,及特色作物连作障碍机理,研发深根系乔—灌—草生物篱等阻控技术,消减土壤连作等障碍因子,提高土壤质量,为该区土壤资源抢救、保护与高效利用的有机结合提供理论和技术支撑,促进区域生态改善与生态富民。

    《断陷盆地石漠化区植被恢复与功能提升》课题,针对断陷盆地生态环境地质分异、植被退化和石漠化程度严重等问题,以建水断陷盆地为主要对象,开展植被退化与恢复生态学机制研究,阐明断陷盆地高原面、坡面和盆地植被退化与恢复生态学机制,以生态系统服务功能提升为核心,研发喜钙、耐旱植物材料筛选与快繁技术、稳定群落构建与可持续性经营管理技术,构建仿自然植物群落与林—灌—草优化配置模式并示范,为喀斯特断陷盆地林草植被生态恢复提供理论和技术支撑。

    《断陷盆地特色生态产业的培育、技术开发与示范》课题,针对断陷盆地光热资源丰富,但石漠化严重、水资源短缺、植被退化严重、产业单一、经济效益低下、生态产业欠发达、农民贫困等问题,以提升生态衍生产业发展为主要目标,选取蒙自市和泸西县作为课题研究区域,研发断陷盆地特色经济价值的林、灌、草植物品种种植技术,构建具生态经济效益的复合经济林—草(药)—畜(禽)等系统和复合水/干果技术模式,形成生态衍生产业模式并示范;研究市场潜力大的生物资源(塔拉等)与断陷盆地地质生态环境的适应机制,研发绿色产业提质增效加工利用技术,培育特色生物资源加工企业,为喀斯特断陷盆地的石漠化综合治理和农民脱贫致富提供科技支撑。

    《断陷盆地石漠化综合治理模式与技术集成》课题,针对断陷盆地石漠化区生态治理—生态产业协同石漠化综合治理模式欠缺、石漠化综合治理效益评估缺乏等关键科学问题,开展断陷盆地石漠化治理技术集成与规范,断陷盆地石漠化治理模式集成与范式,断陷盆地石漠化治理县域示范与效益评价,断陷盆地石漠化区生态环境保护、资源综合开发利用建议与对策研究。集成水土资源高效利用、石漠化生态治理—生态富民耦合调控技术,形成基于生态服务功能提升的石漠化生态富民综合治理模式;对石漠化治理工程进行综合效益、适用性评价,并在县域(流域)范围进行应用示范,为断陷盆地退化生态修复和精准扶贫发挥技术支撑与示范引领作用,同时利用联合国教科文组织国际岩溶研究中心平台作用,积极开展国际科技合作,促进技术共享转化。

     

    5

    重点研发6项关键技术,致力于改善我国老、少、边、穷地区区域生态环境质量

     

    项目将重点研发6项关键技术:

    一是大功率充电联合跨孔CT 成像准确定位地下河技术。

    在断陷盆地区,地下水空间分布及运移规律复杂,难以定位,开采难度大。采用大功率充电联合跨孔彩色CT成像技术,可以降低噪音,提高识别分辨率,对地下水的空间分布进行准确定位。

    二是流域尺度水土资源高效利用与优化调控技术。

    从水资源承载力出发,集成地下水开发、地表水调蓄、雨水收储、水污染风险评估及水质修复、节水灌溉、土壤保墒等技术,达到地表地下水联合调配高效利用。从断陷盆地土壤属性出发,集成有机培肥、保护性耕作、特色作物连作障碍消减等土壤质量提升技术,达到土壤资源高效利用。

    三是深根系乔—灌—草生物篱土壤流/漏失阻控及水土保持生态工程技术。

    以具有丰富溶蚀裂隙喀斯特表层带和落水洞为单元,开展深根系乔—灌—草配置阻控土壤漏失试验,研发相关生物篱技术;引进传统土壤流失防治工程和生物措施,集成断陷盆地水土保持生态工程技术。

    四是喜钙耐旱植物材料筛选与快繁及仿自然植物群落优化配置技术。

    筛选喜钙、抗旱、抗风、耐瘠薄、保水保土固氮能力强等适宜优良植物材料,进行栽培试验;通过种、营养体、最适合环境营造等多技术途径,研发不同植物材料的快繁技术体系。

    五是特色经济林—草(药)—畜(禽)复合系统构建技术。

    重点突破断陷盆地石漠化修复的林下种草养牛(羊)与圈养技术、林—果—药复合经营技术、特色林(果、药)的种苗繁育与优质栽培技术。

    六是生物资源利用提升与生态衍生产业培育技术。

    从塔拉引种、栽培、适应性试验,到生物资源获取、生物产品研制,延长生态产业链,进行喀斯特特色资源利用加工企业的培育。

    曹建华表示,项目参与单位囊括国土资源部、国家林业局、教育部、水利部、中国科学院等部门,覆盖了地学、生态学、林学、土壤学和水土保持学等学科,项目各研究团队均具有深厚的学术研究基础和很强的技术优势,为本项目的顺利执行奠定了坚实基础。

    项目的实施地点与示范区位于我国老、少、边、穷地区,生态环境恶劣,项目将技术研发、试验示范与推广相结合、新技术与已有技术集成和组装相结合、生态效益与社会经济效益相结合,研发技术的应用推广将极大改善区域生态环境质量,使示范区土壤侵蚀模数降低20%~30%,示范区生态服务功能提高30%以上,推动相关生态衍生产业的发展,为生活在生态环境脆弱的断陷盆地区贫困农民提供新的发展机会,对促进少数民族地区的社会稳定具有巨大的社会效益。

    地质扶贫,瞄准断陷盆地精准发力

     

    岩溶地质专家现场研讨

    由贵阳市出发,沿着国道G60上省道S45,一路向南,车行近4小时,是贵州省平塘县。

    贵州省平塘县最近声名鹊起。因为这里建了迄今为止世界最大的单口径射电望远镜FAST工程。贵州大窝凼洼地是典型的喀斯特洼地,一大片漏斗天坑群,它就像一个天然的“巨碗”,可以架设望远镜,而且喀斯特地质条件可以保障雨水向地下渗透,而不在表面淤积,腐蚀和损坏望远镜。

    然而,却很少有人知道,这一全球聚焦的伟大的天文工程却是集中体现了岩溶地质地貌景观综合利用的最典型的范例。更少有人知道,岩溶地质地貌对所在区域经济社会发展产生的系列影响。

    在我国,岩溶地质地貌多分布于南方。包括云南、贵州、广西、湖南、湖北、重庆、四川和广东8个省(区、市)的部分地区,总面积达到78万平方千米。区内多为亚热带气候,湿热多雨。径流时空分布不均,丰水年和枯水年水量之比高达7倍之差。区域内极易发生洪涝或旱灾。

    从区位角度看,西南岩溶区为珠江和长江流域的补给和径流区。珠江流域和西江主干流贯穿该区域南部,主要支流有南盘江、北盘江、红水河、柳江、右江、左郁江等;而长江中上游位于该区域的北部,主要有大渡河、岷江、嘉陵江、乌江、沅江和湘江等。而这些区域恰恰是珠三角、长三角和西部开发的几大战略区域的汇集区,区域位置非常重要。

    相关专家表示,受地形、地貌和岩溶条件的控制,西南地区地表水和地下水在空间上分布极其不均匀,地表水具有支流少、切割深、坡降大等典型特征,对区域性供水十分不利;地下岩溶含水介质则以管道、大型溶蚀裂隙和溶蚀孔隙为主,具有极高的非均质性和各向异性特征。西南岩溶水系统范围较小,且分布分散,地表和地下水的分布、动态和转化特征完全不同于北方地下水系统,这也决定了开发利用条件和方式相较于北方更加复杂多样,应因地制宜,分类开发利用。

    (一)

    长期以来,西南岩溶地区亟待解决的水资源环境问题主要在4个方面:受全球气候变化和特殊岩溶地质条件影响,岩溶区干旱缺水严重,急需开展水文地质调查工作,查明水资源开发利用条件,为当地居民生活和经济发展提供必要的水资源保障。西南岩溶系统对水资源的调蓄能力较差,水土流失严重,导致岩溶石漠化,其综合治理需要水文地质和环境地质综合调查作为技术支撑。岩溶地区地下水质恶化,有部分地下河水质变差,给供水安全和当地居民健康带来安全隐患,急需解决地下水污染防治和水质保护问题。岩溶管道堵塞导致洼地排洪能力不足,岩溶内涝成为当地农村经济发展的制约因素。

    但是,随着城镇化及交通工程建设的快速推进,环境受到的影响也在日趋明显。

    中国地质调查局“岩溶地区水文地质环境地质综合调查”工程首席专家蒋忠诚表示,目前,西南岩溶面临石漠化、岩溶塌陷、植被退化、水土漏失、干旱、内涝、水污染等一系列环境问题。石漠化地区已经无土可流。且气候变化导致的地质灾害、生物入侵、地下水药物污染、作物连作障碍等新的生态环境问题也应高度重视。

    中国地质调查局地质科学院岩溶地质研究所资源研究室主任夏日元也表示,受全球气候变化和特殊岩溶地质条件影响,岩溶地区干旱缺水严重。

    数据显示,2009年秋至2012年,西南地区持续发生严重干旱,仅云南、贵州和广西有1900多万人饮水困难,1亿多亩耕地受旱,直接经济损失 250多亿元。截止2010年,贵州全省仍有1299.8万农村人口饮水不安全,其中水质不达标人口9.4369万人;全省耕地2358万亩,尚有1461.43万亩旱地无水浇灌。而在云南,2010年旱灾造成2702万人受灾,965万人、619万头大牲畜饮水困难,744条中小河流断流、564座小型水库和7599个小坝塘干涸,农作物受灾4740多万亩。在广西,饮水不安全人口1779.64万人,无灌溉设施耕地面积4330万亩;万亩以上岩溶旱区共有83片,总面积643.14055万亩,其中干旱片57片共464.45705万亩,既旱又涝26片共178.6835万亩。

    夏日元表示,岩溶地质的特点是既旱又涝,洪涝灾害频发,给人民生命财产安全带来了威胁,急需查明洪涝分布、规模、产生原因,提出洪涝灾害治理措施或方案。岩溶管道堵塞,洼地排洪能力不足,岩溶内涝成为农村经济发展的制约因素。同时,水土流失与石漠化问题严重。江河上游生态功能区石漠化治理、生态恢复与保护规划需要水文地质和环境地质综合调查技术支撑。目前,石漠化整体上呈改善趋势,已经由1999年11.4万平方千米减少到2015年9.2万平方千米,平均每年减少1.2%。但局部地区石漠化恶化现象严重。岩溶石山地区经济发展缓慢,居民生活困难,为西南地区主要的贫困分布地区。以水资源有效开发利用为龙头,以生态经济综合发展为模式,可以促进扶贫攻坚工作的进程。需要综合集成岩溶地质、水文地质、环境地质和农业地质等地质背景资料,建立水资源开发利用、经济发展和生态环境保护可持续模式。另外需要特别关注的问题是,国家重点生态功能区生态环境恶化。岩溶地下水污染越来越严重,地下河水质变差,给供水安全和人民生命健康带来了安全隐患。急需解决地下水污染防治和水质保护问题。近20年来,部分地下河水质变差,且岩溶地下水深埋于地下,污染具有隐蔽性,处理难度大。

    (二)

    在距离平塘县克度镇不远的地方,早在2005年,由中国地质调查局地质科学院岩溶地质研究所主导的巨木地下河开发示范工程就已经建成应用。

    巨木地下河为大小井地下河系统中游,发源于惠水县抵季乡和羡塘乡,于平塘县塘边镇巨木寨集中排泄。示范工程由地下水开发工程、洪涝洼地治理工程及农业种植3部分组成。其中地下河开发工程又包括了蓄水、提水和输水3个部分。

    平塘县国土资源局和平塘县水利局于2008年联合发布的关于巨木地下河示范工程成果应用评估报告中这样评价道:“地下河出口一下6000亩农田在有了稳定的灌溉水源后,可增收粮食90万千克/年。按照当地市场价格测算,仅农业生产一项就可以实现经济增收171万/年,其中还不包括地下河开发带动的其他经济收益。在巨木地下河流域中,利用石漠化治理地学模式对环境综合治理、石漠化综合治理取得了良好的社会效益、经济效益和环境效益,石漠化得到了有效的控制。”

    事实上,在超过10年的时间里,中国地质调查局地质科学院岩溶地质研究所一直为解决上述问题而努力。

    2003年~2015年,中国地质调查局地质科学院岩溶地质研究所联合西南8省地勘单位及有关科研院所、高校等30多家单位,开展“西南岩溶地区地下水与环境地质调查综合研究”、“西南岩溶石山地区地下水及环境地质调查”、“西南岩溶地区水文地质环境地质调查”、“南方岩溶水文地质环境地质综合调查与整治示范”计划项目,共完成1∶5万水文地质调查面积25万平方千米,综合地球物理探测90000点,岩溶洞穴探测60000米,水文地质钻探68000米。查明了工作区水资源开发利用条件,制订了岩溶水资源开发利用方案和地质环境综合整治区划。选择典型地段进行了岩溶地下水开发利用和生态环境综合治理示范,建立了地下水开发利用和生态环境综合治理模式,解决了620万人饮水困难,为西南岩溶干旱缺水区水资源合理利用、石漠化治理和生态环境保护提供了技术支撑。

    据蒋忠诚介绍,目前,正在实施的国家岩溶环境治理项目《岩溶地区石漠化综合治理规划大纲》(2006-2015) 是以县为单位,共有451石漠化县纳入规划,涉及面积105.1万平方千米,主要实行封山育林育草、人工造林、改良草地、人工种草、发展草食畜牧业、基本农田建设和配套水利水保设施“六大工程”。这一项目从2008年100个试点县算起,到2014年已经扩大到314个县,石漠化治理区面积7万平方千米。中央财政共投入77亿元。云南、贵州、广西、四川、湖南、湖北、重庆、广东石漠化分布的8省区自2000年到2015年石漠化共减少21502平方千米。50%的区域植被覆盖度增加,植被覆盖率提高了6.22%,森林覆盖率增加10.29%。金银花、花椒、火龙果、文山三七等生态种植业得到培植与发展,也带动了很多地区生态旅游产业的兴盛。岩溶水开发工程解决了1000多万人的饮水和一些耕地的灌溉问题。促进了新农村建设的开展;形成了一批综合治理有效模式和经验。

    另据了解,而在其他生态环境治理方面,国家也投入了相当大的力度。21世纪以来,在西南岩溶区还实施了退耕还林工程,长江、珠江防护林工程,“长防”和“长治”工程、“珠治”试点工程,以及国土资源大调查、国家科技计划、抗旱打井工程,各省的土地整理工程等。

    当前,中国地质调查局地质科学院岩溶地质研究所已经形成了6大特色业务,并且成果迭现。一是创立了岩溶动力学。将全球气候变化引入现代岩溶学研究,认识到岩溶作用具有短时间尺度特点,在全球碳循环中具有重要意义。石笋气候变化记录研究达年际精度,解释了岩溶碳循环的主要过程与碳汇稳定性影响因子,提出了人工固碳增汇新途径,可为应对气候变化及实现碳减排做出潜在贡献。二是岩溶生态理论与技术支撑石漠化治理。查明了中国西南石漠化分布、成因、演变及危害,明确了显著石漠化的县域和八大典型区域。建立了广西马山弄拉、平果果化等综合治理示范基地,率先阐明岩溶生态系统、石漠化、水土流失的概念和评价方法,创建了石漠化综合治理模式和技术体系,为岩溶生态修复和精准扶贫提供依据。三是岩溶地下水资源调查为饮水安全提供保障。针对岩溶地区干旱缺水问题开展水文地质综合调查研究,查明了岩溶地下水分布规律和开发利用条件,建立了岩溶地下河系统水循环理论,提出岩溶地下水有效开发利用模式。为岩溶石漠化、地下水污染等环境地质问题治理提出有效措施,探索了水资源利用与生态恢复、经济发展协同途径。通过示范推广直接服务于脱贫攻坚,为饮水安全提供保障。四是岩溶塌陷调查服务岩溶区城镇化和重大工程建设。针对具有隐蔽性和突发性特点的岩溶塌陷地质灾害问题,开展岩溶塌陷调查评价、监测预警和早期识别技术研究,形成中国岩溶塌陷调查评价技术方法体系,为岩溶地区城市规划建设、新型城镇化建设和重大工程建设的地质环境安全提供有力保障。五是岩溶景观与洞穴调查研究支撑世界遗产保护。建立完善了岩溶景观与洞穴发育演化理论和评价方法体系,完成20多处世界遗产、地质公园申报,为岩溶地质遗迹保护和开发利用、生态文明建设和脱贫攻坚提供科技支撑。六是油气藏岩溶储层研究助力油气勘探与开发。以塔里木盆地、四川盆地、鄂尔多斯盆地及渤海湾盆地为研究对象,建立不同地貌单元岩溶发育模式与地震形影模式,研究古岩溶发育分布对油气岩溶储层的控制,建立了碳酸盐岩岩溶储层形成与演化地质模式,实现了对油气聚集有力的岩溶储层预测,对碳酸盐岩油气勘探与开发提供地质理论依据。

    (三)

    中国地质调查局地质科学研究院岩溶地质研究所所长刘同良表示,岩溶地区资源环境问题突出,制约着经济社会的发展。科技创新不仅有助于解决和应对岩溶地区的资源环境问题,也将推动岩溶科学的发展,服务我国生态文明建设。当前,岩溶研究面临一系列科技创新方面的挑战。比如,如何落实地球系统科学在岩溶学中的应用,发挥我国岩溶研究的地域优势,探索我国岩溶关键带的特征和重要过程;如何加强岩溶作用应对全球变化、岩溶碳汇速率和稳定性的研究,建立应对极端气候的长效机制;怎样考虑古纬度和古气候对古岩溶形成的影响;如何深入探索微生物对深部碳酸盐岩岩溶形成的作用;如何系统梳理和总结我国石漠化治理工程的经验和存在的问题,更好地指导下一步的治理工作,思考从南北方岩溶分界线的角度开展岩溶自然遗产地的申报工作;如何将现代大数据等技术运用到岩溶资源环境、水文地质研究中,做好地质灾害的预警预报和应对资源短缺问题的研究,服务国家需求等等。

    夏日元表示,事实上,从解决问题的具体层面,岩溶环境治理也取得了部分重要的科技创新成果。一是建立了岩溶环境监测研究技术体系:自动化监测仪器,自动化监测站,便携式监测仪器。岩溶环境条件、水化学、碳循环观测和试验;岩溶塌陷自动化监测;查明了21世纪以来石漠化环境演变规律。石漠化面积由逐年增加演变为下降,总体趋势好转;区域发展不平衡,局部地区石漠化进一步恶化。二是查明岩溶作用驱动元素迁移,形成特殊岩溶环境。从岩溶动力系统驱动元素迁移的角度,揭示了岩溶环境富钙、偏碱的特点及其制约土壤营养元素的有效性、影响植物和农作物生长的规律,成为因地制宜选择适生树种治理石漠化的理论基础。三是揭示了岩溶环境的空间异质性。但从更高层面来看,岩溶环境研究科技创新仍然面临诸多尚未解决的问题。比如区域生态功能没有形成;多数生态产业缺乏持续性;水土漏失没有得到治理;缺水困难没有根本解决;对新的环境问题缺乏对策等。仅就西南岩溶生态脆弱区地下水而言,水文地质具有岩溶发育与控制因素的复杂性、岩溶分布的非均一性、岩溶含水介质多重性、岩溶水流的多相性、动态变化的剧变性、水资源开发的调节性等特点。由此,解决岩溶地区地下水资源环境问题,必须首先解决天然水资源量大,开发利用率低的问题。目前,西南岩溶地区地下水可开采资源量534亿立方米/年,现开采量仅66亿立方米/年,开发利用潜力巨大。

    目前,中国地质调查局地质科学院岩溶地质研究所针对不同类型区开发条件,提出了堵洞成库、建柜蓄水、抽水调节和束流壅水4种岩溶地下水有效开发利用模式,通过示范推广,直接服务于脱贫攻坚。拦截地下河形成地下调节水库5座,地表溶洼成库2处,钻探成井300多眼,表层岩溶水调蓄20余处,开挖浅井50多口,应急抗旱打井2800多口,解决了620多万人饮用水、100万多亩耕地的灌溉用水问题。例如,岩溶峰林平原区开发利用浅层地下水与节水生态农业模式。利用峰林平原浅层岩溶化强、储水性和透水性好的特点,在农作物需水的干旱季节,开采地下岩溶径流带水资源,形成调节空库容,雨季来临后通过天然降水入渗而恢复地下库容;通过喷灌、滴灌、移动式灌溉等节水灌溉技术的应用,提高水资源利用率。宾阳县黎塘镇谢村桥美示范区,通过开发表层岩溶系统地下水,建设以蔬菜为主的高效农业基地,年产值达1300多万元,增收20%~30%。再如溶丘洼地区堵洞引水与庭院经济模式。在地下水补给区的高位洼地,选择岩溶管道集中发育的部位,堵洞形成地表——地下联合水库;径流区开凿隧洞,拦截地下河,引水灌溉和发电。贵州省巨木地下河岩溶水开发示范工程。在地下河出口处筑坝拦蓄地下水成库,库容239.2万立方米地下水库、300立方米的高位蓄水池、水轮泵站。解决了区内5000余人和10000余头大牲畜的饮水、缓解了6000余亩农田缺水灌溉。

    (四)

    不可否认,所有环境地质问题都与人有关。所以,一切地质环境研究都必须“以人为本”。

    国家发展研究院瑞意高等研究所城市与环境学院土地科学中心蔡运龙教授就表示,岩溶地区的地质环境当前面对诸如岩溶区地下水循环转化、岩溶地下水污染、水土流(漏)失与石漠化、石漠化综合治理与生态修复、岩溶生态脆弱区水资源开发利用等多个问题,但可以更广义地归纳为两大问题:即石漠化治理模式的可推广性和石漠化治理模式的可持续性;退化土地生态重建的社会工程途径。但是,必须意识到,生态系统退化表面上是生态环境问题,归根结底是社会发展问题,要从社会经济方面寻求根本的解决途径。只有社会经济发展了,才能改变“土里刨食”那样对土地(生态系统)的严重依赖,进而解除人类活动对生态系统的干扰。

    他认为,欲使农村人口尤其是岩溶分布区的农村人口摆脱掠夺土地(生态系统)的状况,需要提高人口素质,掌握其他的生存和发展手段。教育与培训的作用至关重要。生态修复工程项目的投资不要只用在“地”上,而需要把更多的财力和科技投入到“人”的能力建设上。

    当岩溶地区遭遇科技克星