岩溶多重介质地下水模拟技术及应用进展

由于岩溶发育的极不均匀性和岩溶含水介质的多重介质性等特征, 使得岩溶水系统的地下水运动规律呈现出复杂多变的特点, 在孔隙或微小裂隙含水介质中岩溶水满足达西流, 而在岩溶管道中可能出现非达西运动, 这给岩溶区地下水模拟带来了极大挑战。系统地梳理了岩溶水流和溶质运移模拟的主要方法, 总结了复杂岩溶多重介质系统地下水模拟技术的现状, 并结合我国南方岩溶和北方岩溶各自的特点, 提出岩溶水模型在实际应用中需要关注的要点。主要包括: 需加强野外调查和观测, 提高模型中对岩溶含水介质结构非均质性表征的精度; 进一步深化对岩溶水运动机理的研究, 刻画岩溶多重介质的水流特征以及水交换机理; 针对北方岩溶水系统模拟, 可选择等效连续多孔介质模型, 重在耦合区域分布的岩溶溶孔－裂隙介质中慢速渗流与脉状分布的强径流带快速流; 对于南方岩溶区, 建议考虑建立分布式流域水文模型与能刻画集中管道流多流态变化的多重介质模型的岩溶地下水耦合模型。      

基于岩溶演化模型的隧道突水危险性评价

岩溶隧道是否发生涌突水,主要取决于洞身岩溶发育程度和岩溶系统结构类型,而这与岩溶系统的演化过程密切相关。岩溶系统的发育演化始于具有侵蚀性的水对可溶性岩体中裂隙的溶蚀扩展。基于此,本研究耦合了渗流模型和溶蚀扩展模型,并结合鸡公岭隧道区水文地质背景资料,建立了隧道区岩溶系统岩溶演化模拟模型,模拟再现了隧道区岩溶含水介质的发育情况。结果显示:随着岩溶系统的发育演化,整个系统中裂隙不断被溶蚀扩展,差异性溶蚀越来越显著,岩溶含水介质的非均质性越来越强;同时随着埋深的增加,岩溶发育程度及非均质性越来越弱;浅部和深部裂隙的平均溶蚀扩展率相差近1 500倍;鸡公岭隧道洞身标高处岩溶发育微弱,裂隙溶蚀扩展较小,仍为裂隙岩体,其等效渗透系数为0.51m/d,岩溶突水概率极低,隧道发生涌水时的涌水量为127.9m3/d,对隧道施工影响较小。      

岩溶洼地内涝蓄水量与不同级次裂隙对溶质迁移影响的室内实验与模拟

岩溶洼地在西南地区分布广泛,岩溶地下水为当地居民生产和生活的主要水源,研究岩溶洼地内涝条件下裂隙系统中的溶质迁移及响应规律可为峰丛洼地区岩溶地下水的污染防控提供科学依据。基于相似原理构建了岩溶洼地-裂隙系统物理模型,采用对流-弥散模型和脉冲响应模型对溶质迁移过程进行了模拟分析,探讨了不同级次裂隙和内涝蓄水量大小对溶质迁移的影响。结果表明:大裂隙作为优先水流通道,对溶质迁移起主要的控制作用;小裂隙对溶质迁移起暂时储存和缓慢释放的调蓄作用,其溶质迁移的叠加过程使总出口溶质迁移的"拖尾"效应更加明显;随着内涝蓄水量的增加,稀释作用增强,导致溶质峰值质量浓度降低、溶质质量浓度回归至本底值的时间延长、延缓释放溶质的质量增加。该物理模型是对复杂岩溶洼地系统内污染物迁移研究的一种新探索,其模拟结果进一步深化了对洼地内涝条件下面源污染物在岩溶水系统中迁移规律的认识。      

澜沧江上游河谷卸荷松弛区地下水位分布成因

澜沧江上游两岸山体普遍存在强卸荷松弛和地下水位低缓现象,深入研究岸坡卸荷松弛区水文地质特性可为开展水电工程运行期渗流特征分析和渗控设计提供依据。澜沧江深切河谷区地下水位低缓成因复杂,选取上游某水电站坝址区为研究对象,从区域条件和局部条件综合分析地下水位低缓的原因,构建研究区考虑卸荷松弛分带的水文地质结构和地下水流模型,模拟与验证自然条件下地下水渗流规律。结果表明：研究区地下水位低缓既受控于区域地下水流与澜沧江的共同作用,又是河谷岸坡坡度、表层滑带、岩层结构、卸荷裂隙等局部水文地质条件阻碍降雨入渗补给地下水的内部表现;降雨径流受地形控制,主要经地表及强卸荷裂隙排泄至澜沧江,降雨入渗补给仅占地下水总补给量的5.50%。     

3D GIS下地下水流有限元数值模拟方法与应用

目前缺乏针对较强空间特性水文地质计算参数的高效采集方法与管理机制,而且现有孔隙地下水流有限元数值模拟在三维可视化与空间分析等方面仍存在很多不足。因此,本文以有限元数值计算方法和三维地理信息平台为基础,结合GIS空间分析算法和计算机图形学理论,针对孔隙承压地下水有限元分析过程的关键步骤（概念模型构建、空间离散、水文地质参数提取与赋值等）,提出了孔隙地下水有限元数值模拟过程在3D GIS下的实现方法和技术框架。基于水文地质钻孔数据和GTP（Generalized Tri-Prism）空间数据模型,以盐城市滨海平原水文地质区为研究区域,实现了地下水系统空间结构特征的三维表达,同时与地下水流有限元数值模型相耦合,提出了3D GIS模型与地下水非稳定流有限元数值模型耦合的核心处理算法,并在此基础上开发了3D GIS地下水流数值模拟系统。应用实例表明,该技术框架可有效简化有限元分析流程、优化模型计算效率、实现地下水有限元数值模拟过程及计算结果的三维可视化。     

基于离散裂隙网络法和水流数值模拟技术的地下水封洞库水封性研究

水封性是保证地下水封洞库安全运营的关键，目前研究地下水封洞库水封性的主要难点在于裂隙岩体含水层强烈的非均质性和各向异性。以某地下水封洞库工程为背景，通过实测裂隙几何参数分析裂隙发育特征，建立了研究区离散裂隙网络模型。之后对研究区进行网格化划分计算了各单元体渗透系数张量，基于对裂隙渗流基本规律的理解和假设建立了研究区高精度的非均质各向异性渗流模型，并对研究区渗透性特征进行分析。通过建立非均质各向异性渗流模型对水幕系统水封性能进行了探讨并与传统均质各向同性渗流模型进行对比分析，利用该模型可以弥补传统均质各向同性渗流模型的不足，更加准确地评价由于裂隙岩体非均质各向异性造成的储品泄漏风险。      

基于无人机影像和面向对象随机森林算法的岩溶湿地植被识别方法研究

目前对岩溶湿地的重视程度远低于其他湿地类型,缺乏利用遥感技术进行岩溶湿地植被高精度识别的研究,但岩溶湿地同其他湿地类型一样,湿地面积退化严重,亟待需要解决。因此,本文选取受人类活动影响较大、湿地退化较为严重的广西桂林会仙喀斯特国家湿地公园的部分核心区域作为研究区,以DJI大疆御Mavic Pro无人机航摄影像为数据源,利用泛化能力强、分类精度高的面向对象随机森林算法实现了会仙岩溶湿地植被的高精度分类,探究无人机可见光影像和面向对象随机森林算法在岩溶湿地植被识别中的适用性,为无人机遥感技术应用于岩溶湿地的研究和保护提供技术参考。首先,在eCognition Developer9.0中利用多尺度迭代分割算法对影像图层进行分割;然后,基于以往在进行面向对象分类研究的经验来指导我们进行特征选择,充分考虑了影像的光谱和纹理特征、植被指数、无人机遥感数据派生的研究区数字地表模型（DSM）和几何特征;最后,在RStudio中实现了随机森林算法参数的调优、模型的构建以及分类。结果显示,面向对象随机森林算法对会仙湿地植被具有较高的识别能力,在95%置信区间内总体精度为86.75%,Kappa系数为0.83。在单一典型岩溶湿地植被识别精度中,狗牙根-白茅-水龙植被群丛的用户精度在90%以上,生产者精度高于80%,竹子-马甲子-桂花生产者精度高于80%,但是用户精度较低,仅为70.59%。     

平阴县前寨—凌庄水源地地下水数值模拟

根据平阴县前寨一凌庄水源地的水文地质条件，建立了裂隙岩溶水的平面二维非稳定地下水流模型，并应用Modflow软件对该区地下水流进行数值模拟。在进行计算区剖分、空间离散、参数分区、边界条件及源汇项处理与模型识别、验证的基础上，通过预测分析得出：该水源地按设计开采量开采后，地下水位不会出现持续下降趋势，最大降深在允许范围之内，因而也不会诱发环境地质问题。     

黄陵穹隆周缘岩溶水流系统特征及成因

岩溶水流系统特征研究有助于地下水资源的合理评价和开发利用。借鉴水文学的研究方法，统计并量化了典型岩溶水流系统的空间特征以及其水文动态响应、温度场和电导率特征。划分了扇状、树枝状、平行状、梳状4种地下水系来综合反映岩溶水流系统的地表-地下岩溶特征，前两者主管道垂直于地层走向，构造裂隙起汇水作用，后两者主管道平行于地层走向，层面裂隙起汇水作用。黄陵穹隆西北翼、西翼和南翼以平行状和树枝状为主，东翼和北翼则以扇状水系和平行状水系为主。不同地下水系结构的形成及区域差异与含水系统和水系的空间关系和级次性密切相关，并表现出不同的动态特征。扇状和平行状岩溶水流系统对降雨响应最为敏感，而梳状水系岩溶水流系统响应和衰减过程最慢；基于岩溶地下水温度与出露高程和循环深度显著相关的关系建立了鄂西山区地下水温度线。这一基础性研究可为岩溶地下水流系统研究和当地工程实践提供一定的理论支持。      

涿鹿盆地三维多裂隙地质模型地温场数值模拟

地温场预测是地热资源储量评价的重要前提,深部裂隙岩体渗流-传热耦合温度场研究正受到越来越多的关注。涿鹿盆地是新生代地堑式断陷盆地,地表温泉出露点多,地热资源丰富。依据实际地质调查结果和3条地球物理勘探测线解译的地质构造剖面,建立了矿场尺度的涿鹿盆地三维不平行多裂隙地质模型和地热成因模型,赋值调整参数,得到考虑地下水流动情况下多裂隙岩体渗流-传热数学模型,模型的合理性得到验证后,模拟计算了盆地内部温度场,指出盆地地热异常区集中在盆地中心蓟县系（70℃左右）和盆地北东长城系部位（90~98℃）,最后研究分析了水力梯度的变化对盆地内部温度场的影响规律。研究结果为涿鹿盆地或同类型地区地热资源勘探、评价和开发利用提供一定的依据。      

基于权重线性加和模型的城市污泥土地利用适用性评估

城市污泥土地利用的环境问题是限制其大规模利用的主要因素。本文建立基于权重线性加和模型的城市污泥土地利用环境风险控制方法,提出北京市城市污泥土地利用环境风险控制方案。不同因素对城市污泥施用的环境影响依次为：土壤重金属含量（0.22）>土地利用类型（0.17）≈土壤类型（0.17）≈自然降雨（0.17）>坡度（0.13）>与自然水体距离（0.09）>与城镇居民区距离（0.05）。北京市城市污泥低风险适宜施用区域主要集中在平谷和顺义交界处、昌平-延庆中部以及房山的东部山间林地,施用面积为2033 km²。中风险施用区分布在西南部、东南部的林地和旱地的混合区域,面积为5079 km²。高风险施用区面积为380 km²,分布在石景山以及门头沟东北部、房山西南部以及平谷北部区域。禁止施用区面积达8916 km²,主要分布在城区及城区周边的郊县、延庆、怀柔、密云等部分区域。     

济南市长清区地质灾害易发程度评价及防治分区划分

济南长清区地质灾害发育较强烈,包括崩塌、滑坡、泥石流和岩溶塌陷等4种类型。 该文在现状调查的基础上,采用“地质灾害综合危险性指数法”,以地质、地形地貌、气候植被、地质灾害隐患点、地质灾害规模、分布密度、活动频次和险情等因素为评价因子,将长清区地质灾害易发程度划分为中易发区、低易发区及不易发区3个区,并在此基础上进行了防治分区划分 ,为地质灾害的预防和治理提供了科学依据。     

康定北部高原构造岩溶发育特征与地下水径流带识别

受自然环境和技术方法制约，青藏高原岩溶发育演化和岩溶地下水循环特征研究相对薄弱，制约了青藏高原碳酸盐岩区的经济发展、民生设施建设和地质灾害防治。通过野外地质测量，岩溶地下水、地表水和大气降水水化学和同位素特征分析，泉水流量动态，水均衡计算和物探等技术方法，系统分析了四川省康定市北郊碳酸盐岩分布区的岩溶发育特征，识别了岩溶径流通道和岩溶大泉主要补给来源。结果表明:康定市北郊碳酸盐岩分布于高山峡谷地貌类型区，可溶岩地层分布、岩溶发育程度和岩溶水补给、径流、排泄均受构造控制，可溶岩与非可溶岩接触带和活动断裂附近的岩溶发育程度较强。岩溶水呈管道流径流，主要以岩溶大泉形式集中排泄，泉流量约1.5×104 m3/d且动态较为稳定。通过水文地质条件分析，识别出研究区存在通化组岩溶水径流带和雅拉河断裂岩溶水径流带。水化学-同位素数据、岩溶泉流量动态和水均衡计算结果显示，雅拉河河水是岩溶大泉的主要补给源，岩溶地下水主要沿雅拉河断裂岩溶水径流带径流并集中排泄。      

山东省地质灾害分区与防治对策研究

根据山东省地质环境条件、地质灾害分布与发育规律,结合灾害所造成的经济损失,将山东省分为7个主要地质灾害区(地面沉降、地裂缝、采空塌陷、岩溶塌陷、崩塌滑坡泥石流、崩塌泥石流、海咸水入侵)、26个地质灾害地段。通过调查及综合研究,对各区段地质灾害特征和发生规律等进行了概略评价;提出了地质灾害防治的行政性对策建议,并针对不同类型地质灾害提出了技术性对策。     

沂蒙山区典型断陷盆地岩溶地下水系统特征: 以莱芜盆地为例

2011年以来，中国地质调查局在沂蒙山区组织实施了1:5万标准图幅水文地质调查4万余km2，并在严重缺水村镇开展了大量的找水打井示范工作，获得了较为丰富的地质数据，对断陷盆地地下水流系统取得了新认识。受中新生代构造影响，沂蒙山区发生断裂褶皱、伸展滑脱及岩浆侵入活动，形成一系列"南超覆北断陷"的地堑-半地堑盆地，并最终形成现今典型的"盆-山"岩溶水文地质结构，及以盆地为单元的相对独立的岩溶地下水流系统。为研究沂蒙山区地下水流系统发育特征，选择莱芜盆地为典型研究区，基于野外地质调查，通过综合分析盆地南北两侧地貌单元、含水岩组立体空间分布及地下水水位、水化学及同位素结果等，探讨了沂蒙山区典型"盆-山"结构塑造的多级岩溶地下水流系统特征。结果表明，莱芜盆地岩溶地下水由盆地外围向盆地中心呈"向心式"径流；受人类活动影响、地质构造控制和含水层分布制约，盆地南北两侧地下水流系统特征存在差异:盆地南部发育中间和局部两级地下水流系统；盆地北部则仅发育单一的局部地下水流系统，但占有已勘查论证的近1/2的地下水水源地。此外，研究发现在大汶河最低侵蚀基准面，区域滑脱构造及热液混合作用拆离的空隙，与岩层面、层间裂隙、顺层溶蚀空隙等共同构成立体的岩溶地下水网络，影响着岩溶水循环途径及深度。以此认识为指导，实施的探采结合井成井率达到86%，强化了地下水流系统理论在北方基岩山区水文地质工作中理论与实践的结合。      

碳酸盐岩古溶洞层级约束地质建模方法探讨:以塔河油田奥陶系某缝洞单元为例

碳酸盐岩地下古溶洞及其内部充填结构一直以来是地质建模的难点。综合地震、地质、测井等资料,提出了溶洞复合体、单一溶洞及溶洞内部充填结构3个层次逐级约束建模的思路,探索了基于多点地质统计学的溶洞内部充填结构建模方法。其中溶洞复合体属于较大尺度地质体单元,地震分辨率基本可识别,井震联合反演约束下采用截断高斯模拟方法建立溶洞复合体模型;单一溶洞是溶洞复合体的重要组成单元,属于中等尺度地质体单元,地震分辨率难以直接识别,基于大量露头实测数据,在溶洞复合体约束下采用基于目标的方法建立单一溶洞分布模型;溶洞内部充填结构属于小尺度地质体单元,地震分辨率很难识别,在单一溶洞约束下,尝试采用多点地质统计学方法建立溶洞内部充填结构模型。以上层级约束建模思路有助于提升地下古溶洞及其内部储集体地质建模的准确度,在塔河油田某缝洞单元测试中取得了初步效果,同时讨论了存在问题及攻关方向。该研究对于精细表征缝洞型碳酸盐岩油气藏具有指导和启示作用。      

乌蒙山地区岩溶地下水流系统结构及其找水应用

为推动乌蒙山贫困缺水区生态环境建设、地下水资源合理利用,以昭通市作为重点调查区,基于地质调查、泉流量统计、水质检测,展开岩溶地下水富集规律及物探找水方法的研究。结果表明:①研究区岩溶地下水系统以条带岭谷型、埋藏型为主。条带岭谷型岩溶地下水系统以现代岩溶为主,目标含水层多、发育深度有限、岩溶发育及富水程度差异大,富水块段为背斜核部及两翼、向斜核部、断层影响带;埋藏型岩溶地下水系统以古岩溶为主,目标含水层单一、发育深度较深且极不均匀,富水块段为断陷谷盆埋藏的古岩溶。②地下水化学类型以HCO₃型、HCO₃·SO₄型为主,条带岭谷型、埋藏型岩溶水分别占96.73%,92.93%,水质总体较好,综合水质评价Ⅰ~Ⅲ类水占比分别为80.84%,64.41%。③建议对>50 L/s大泉进行提引、丰储冬用,对富水块段进行综合物探探测和钻探验证的方法找水。④综合物探找水方法:先通过高密度电法、联合剖面法查明岩溶破碎带及断层,再激电测深确定极化率高的含水层,最后综合测井和钻孔揭露确定具体出水段及涌水量,找水成功率为86.67%,适用于条带岭谷型及浅埋岩溶地下水。      

临沂市柳航头地区地热资源特征研究

区域地质调查和地热地质勘探资料显示,临沂市柳航头地区区域构造发育较强烈,大气降水通过构造破碎带向地下深部运动,深部寒武-奥陶纪灰岩岩溶地层厚度可达1 500余米。区域地下水在深部渗透性较好的岩溶地层内径流,向柳航头地区运动过程中,不断受到地球深部热源的地温加热而形成地热水,寒武-奥陶纪岩溶地层成为热储层,其上覆地层导热率较低,是较好的盖层。结合研究区内的地热井深部测温数据和地热水化学特征数据,建立柳航头地区地热概念模型。     

Simulation of land-use scenarios for Beijing using CLUE-S and Markov composite models

This study investigated and simulated land use patterns in Beijing for the year 2000 and the year 2005 from the actual land use data for the year 1995 and the year 2000,respectively,by combining spatial land allocation simulation using the CLUE-S model,and numerical land demand prediction using the Markov model.The simulations for 2000 and 2005 were confirmed to be generally accurate using Kappa indices.Then the land-use scenarios for Beijing in 2015 were simulated assuming two modes of development:1) urban development following existing trends;and 2) under a strict farmland control.The simulations suggested that under either mode,urbanized areas would expand at the expense of land for other uses.This expansion was predicted to dominate the land-use conversions between 2005 and 2015,and was expected to be accompanied by an extensive loss of farmland.The key susceptible to land-use changes were found to be located at the central urban Beijing and the surrounding regions including Yanqing County,Changping District and Fangshan District.Also,the simulations predicted a considerable expansion of urban/suburban areas in the mountainous regions of Beijing,suggesting a need for priority monitoring and protection.