地球信息科学学报 ›› 2016, Vol. 18 ›› Issue (2): 182-189.doi: 10.3724/SP.J.1047.2016.00182

• • 上一篇    下一篇

基于圆盘概化模型的岩溶裂隙渗流路径构建算法

李祥琛(), 芮小平*(), 曲小康, 宋现锋, 王慧芳   

  1. 中国科学院大学资源与环境学院,北京 100049
  • 收稿日期:2015-06-09 修回日期:2015-08-10 出版日期:2016-02-10 发布日期:2016-02-04
  • 通讯作者: 芮小平 E-mail:lixiangchen13@mails.ucas.ac.cn;ruixpsz@163.com
  • 作者简介:

    作者简介:李祥琛(1990-),男,河北枣强人,硕士生,研究方向为地理信息系统开发和地理数据挖掘。E-mail:lixiangchen13@mails.ucas.ac.cn

  • 基金资助:
    基金项目:国家环境保护公益性行业科研专项项目“环保地下水污染监控预警与事故应急技术体系研究”(201409030)

Algorithm of Karst Fracture Seepage Path Construction Based on Disc Generalized Model

LI Xiangchen(), RUI Xiaoping*(), QU Xiaokang, SONG Xianfeng, WANG Huifang   

  1. College of Resource and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
  • Received:2015-06-09 Revised:2015-08-10 Online:2016-02-10 Published:2016-02-04
  • Contact: RUI Xiaoping E-mail:lixiangchen13@mails.ucas.ac.cn;ruixpsz@163.com

摘要:

岩溶裂隙网络的识别与重构,一直是岩溶地下水资源和岩溶区地质环境保护研究的热点与难点。基于地表原始裂隙数据、经过克里格插值后的地表裂隙数据、以及通过蒙特卡洛预测得到的地下裂隙数据,本文针对三维裂隙构造中的水流通路识别模块,提出了一种基于圆盘模型的渗流路径构建方法,采用有向图的数据结构来模拟圆盘模型概化后水流沿裂隙的渗流路径,并以邻接矩阵的形式存储其渗流路径数据。同时,为满足大规模裂隙数据快速求交的需要,本文给出了三维R树索引算法缩减遍历各个采样点所用的时间,提高了整套算法运行的效率。最后,以北京市房山区张坊地区地表采样裂隙数据为基础,借助遥感地质调查、重点岩溶区精细地质测量及取样分析等手段,开展三维裂隙空间网络分布模型研究,对岩溶裂隙渗流路径构建算法进行计算机仿真实现,为张坊地区岩溶发育机理的研究及岩溶水裂隙数值模拟提供可视化的分析手段。

关键词: 岩溶区, 裂隙, 圆盘模型, 渗流路径

Abstract:

The identification and reconstruction of karst fracture network has always been a hot and difficult topic in the research of karst groundwater resources and geological environment protection. Based on the original surface fracture data obtained in Zhangfang, Fangshan district in southwest Beijing and the surface fracture records processed through Kriging interpolation method, a construction method of the flow path based on disc model is proposed, focusing on the identification module of flow passage in three dimensional fracture structure and the underground fracture data obtained by Monte Carlo method. The flow path along the fracture generalized by disc model is simulated using directed graph data structure, and is stored through the form of adjacency matrix. In order to meet the demand of fast intersection operation in large scale fracture data and to improve the operation efficiency of karst fracture simulation system, this study gives a three dimension R-tree index algorithm to reduce the time of traversing each sampling point. Finally, based on the sampling surface fracture data of Zhangfang, Fangshan district in southwest Beijing, with the assistance of remote sensing geological survey, fine geological survey of key karst areas, and sampling analysis, this pager commits to the studies of the network distribution model for three dimension fracture space. The algorithm for building the path of karst fractures has also been simulated by computer. It provides a visualized analysis method for the study of karst development mechanism and the numerical simulation of karst water fracture in Zhangfang area.

Key words: karst area, fracture, disc model, seepage path