含水层非均质结构的马尔可夫链地质统计方法及应用含水层非均质结构的马尔可夫链地质统计方法及应用.pdf

2 2 2 2 第 卷 　第 期 地 质 科 技 情 报 　 25 5 　 　　　 　 　　　Vol. 25 No. 5 2006年 　　9 月 Geological Science and Technology Information Sep. 　2006 含水层非均质结构的马尔可夫链 地质统计方法及应用 李 　杨a ,唐仲华b ,方 　琼a (中国地质大学 研究生院 环境学院 武汉 ) 2 2 a. ;b. , 430074 摘 　要 :针对含水层非均质性对地下水流动和溶质运移模拟的准确性有着重要的影响 ,但非均质性的刻画十分困难的问题 ,介 绍了马尔可夫转移概率理论方法及其在含水层非均质结构研究中的应用 ,利用 Walther 定律建立了三维马尔可夫链 ,并将该方法 用于研究位于山西省介休市龙凤河冲积扇岩相的空间分布规律。结果表明该地区细颗粒具有在垂直方向上向上沉积、水平方向 上向外沉积的趋势 ,细颗粒岩相与粗颗粒岩相之间的关联性强而粗颗粒岩相之间的关联性差。该方法的缺点在于 Walther 定律的 应用可能会增加不确定性。 关键词 :转移概率 ;马尔可夫链 ;含水层 ;非均质性 中图分类号 :P641. 2 　 　　　　文献标识码 :A 　 　　文章编号 :1000 7849(2006) 05 0092 05 　① 　　多孔介质体的形成经历了复杂的物理化学、成 的调整缺少客观标准 ,从而增加了模型建立的不确 岩作用过程 ,因此含水层在空间上常表现出非均质 定性。 ②基于交叉变差函数的地质统计方法主要以 性 ,准确掌握含水层结构的非均质性对地下水流动 钻孔岩芯数据为基础而很少考虑地质信息。③交叉 和溶质运移问题的研究十分重要 。 变差函数中恒有 rij ( h ) = rij ( - h ) , 所以它不能反 目前研究含水层非均质性的方法较多 , Kolter 映地质空间分布上的不对称性及其单调分布趋势。 mann 等[ 1 ] 对这些方法进行了详细的介绍和归纳 。 为弥补传统的基于交叉变差函数的地质统计方 其中传统的空间统计方法是一种基于交叉变差函数 法存在的上述不足 ,将基于马尔可夫链的地质统计 的地质统计方法 ,它以空间相关性和连续性模型为 方法引入到含水层非均质性的研究中。 基础 ,从已知的钻孔资料 、地质信息推测未知区域含 对于基于马尔可夫链的地质统计方法的理论研 水层水力特征的分布状况。基于交叉变差函数的地 究及其在刻画含水层非均质性方面的应用国外早有 质统计方法是以协同区域化理论为基础 ,以交叉变 报道。如 Carle 等[5 ] 对该方法进行了详细的论述 ,