基于随机建模技术水文地质岩相模拟初探.pdfVIP

基于随机建模技术的水文地质岩相的模拟初探 杜慧丽 河海大学土木工程学院，江苏南京 (210098) E-mail ：dhl1110@163.com 摘 要：对于溶质运移等小尺度的水文地质模型，水文地质参数在土壤孔隙介质空间分布的 变异性会严重影响地下水流动与溶质运移，但实际资料通常无法提供模型所需要的详细的水 文参数说明。因此如何适当地估测地下水溶质传输模拟所需的水文参数在空间的分布成为水 文地质领域的重要研究课题。针对此种情况，本文运用地质统计学方法中的序贯高斯和模拟 退火算法分别进行模拟，可以根据钻孔资料，建立小尺度的水文地质相模型，并比较分析了 两种算法的模拟效果，探讨两者对地质统计模型构建的影响，为进一步进行小尺度的水文地 质模型的计算研究提供一种新的思路方法。 关键词：地质统计学，空间变异性，模拟退火，序贯高斯，随机建模 1. 引言 通常水文地质计算中是将含水层视为均匀水平分布，水文地质参数取为常数。这主要是 因为在当前技术条件下，受限于经费与现实条件难以掌握任意尺度的水文地质参数的真实特 征，所以采取的简化方法。 实际上由于地质作用的影响，水文地质参数在土壤孔隙介质的空间分布具有空间变异 性。此变异性严重影响地下水流动与溶质运移的行为，因此计算地下水流动与溶质运移时， 地层差异所造成水文地质参数的空间变异性必须加以重视。如何适当地估测地下水污染传输 模拟所需的水文参数在空间的分布成为水文地质领域的重要研究课题。 本文使用地质统计学中的模拟退火算法和高斯序贯算法进行水文地质相的模拟，根据钻 孔资料分析得出水文地质相的分布规律，随机模拟满足分布规律的水文地质相，为分析探讨 水文地质参数空间变异性对水流和污染物运移方面的影响提供了新的思路。 2. 地质统计学方法简介 2.1 区域变量理论 地质统计学是以空间随机函数Z(X)表示任意和空间坐标有关的参数，称为区域性变量， 其中 X=(x,y,z)代表点所在的空间位置。一般而言，区域性变量具有两类主要特征，第一为 随机性，即所在分析区域中任意点的观测值，皆具有不确定性。第二为相关性，即区域性变 量在不同位置的观测值可能具有某种程度的相关性。在地质统计学方法中，最常用的是以变 异函数作为相关性的量化表示[1] 。 2.2 半变异函数与半变异图 半变异函数指区域化变量的增量平方的数学期望之半, 即区域化变量的方差之半，其定 义如下[2] ： N (h ) 1 ⎧ 2 ⎫ =+ − γ h z x h z x ( ) 2N (h ) ⎨∑[ ( i ) ( i ) ] ⎬ (1) ⎩i 1 ⎭ z (x ) z (x +h) h i xi i 式中 为在间距范围内配对点间的距离， 是位于 的区域变量值， 是位 x +h N (h ) x x +h h i i i 于 的区域变量值， 为 和 之间距离 的配对数。以半变异函数与距