基于变异函数的岩体各向异性研究.docxVIP

基于变异函数的岩体各向异性研究 水利水电工程的渗透、地下水资源的开发和利用以及岩石边坡的稳定性评价经常与岩石裂缝的渗透有关。因此,在地下水渗透参数的模拟中,对岩体渗透性参数空间分布规律准确的认识和描述具有重要的实际意义。裂隙岩体渗透性的空间分布具有结构性和随机性特征,地质统计学方法中区域化变量能够较好地描述变量的结构性和随机性。非参数地质统计学主要依赖于变异函数模型,因而变异函数理论模型的选取极为重要。目前,在裂隙岩体渗透性研究领域,主要通过球状变异函数模型来研究岩体渗透性。本文引入孔穴效应模型,运用该模型对不同方向上的渗透性进行结构套合,并以锦屏水电站坝址右岸的岩体渗透性为例运用指示Kriging方法进行了估值和分析。结果表明了指示Kriging方法在确定水文地质参数方面具有良好的效果。 1 khig理论 1.1 变异函数理论 Ι(X,Ζc)={1,当Ζc-1Ζ(X)≤Ζc0?否则I(X,Zc)={1,当Zc?1Z(X)≤Zc0?否则 式中:Zc称为边界值,I(X,Zc)为指示函数。 变异函数既能描述区域化变量的空间结构性,也能描述其随机性,是描述地质参数方面的好工具。设空间某点位置为X(x,y,z),则区域化变量Z(X)在X,X+h两点处的值之差的方差之半定义为Z(X)在空间上的变异函数,记为γ(X,h)。在地质统计学研究中,常作本征假设(intrinsic assumption),此时变异函数与空间位置X无关,而仅与h(滞后距)有关,则指示变异函数定义为: γΙ(h,Ζc)=12E(Ι(X+h,Ζc)-Ι(X,Ζc))2γI(h,Zc)=12E(I(X+h,Zc)?I(X,Zc))2 实验指示变异函数为 γ*Ι(h,Ζc)=1Ν(h)Ν(h)∑i=1(Ι(Xi+h,Ζc)-Ι(Xi,Ζc))2γ?I(h,Zc)=1N(h)∑i=1N(h)(I(Xi+h,Zc)?I(Xi,Zc))2 式中:N(h)是相距|h|的数据点的数目。 变异函数的理论模型可分为有基台值和无基台值模型两类。有基台值的模型包括球状模型和指数函数模型等,无基台值的模型包括纯块金效应模型和孔穴效应模型等。如果变异函数值不是单调的增大而是呈现上下波动时,则该变异函数就表现出孔穴效应。常用的孔穴效应模型的公式为 γ(h)=C0+C1(1-e-|h|/acos(h/b))(1) 式中:e-|h|/a起控制振幅的作用,h/b控制着孔穴效应的频率。 对于三维各向异性可以通过不同方向上的结构套合为各向同性。首先,确定各向异性比ai/aj=k(常数),使任意两方向的变程相同。把总的套合结构γL(h)看成是在一个各向同性结构γ1(h)基础上叠加了一个在i方向上多出来的附加结构γ2(hi),即 γL(h)=γ1(h)+γ2(hi)=γ1(h)+γ(hi)-γ1(hi)(2) 套合结构的一般表达式为 γL(h)=γ1(h)+∑iγ(hi)-γ1(hi)(3) 1.2 基于约束信息的非线性动力系统的方程 以x0为中心的待估区的区域化变量I*(x,Z),即由定义在点承载xi上的n个指示函数值I(xα,Z)来确定,其概率估计量为 Ι*c(x,Ζ)=n∑α=1λα?Ιc(xα,Ζ)(4) 式中的权系数λα可由指示Kriging方程组求出,其方程组及方差表示为 {n∑β=1λβγΙ(xα,xβ;Ζ)+μ=ˉγΙ(xα,V;Ζ)n∑α=1λα=1(α?β=1,2,??n)(5)σ2ΚΙ=n∑α=1λαˉγΙ(xα,V;Ζ)-ˉγΙ(V?V;Ζ)+μ(6) 式中:μ为拉格朗日乘子;γI(xα,xβ;Z)、ˉγI(xα,V;Z)分别代表在给定边界信息Z条件下,各信息点xα与xβ之间指示变异函数值,以及待估区域V与信息点xα之间的平均指示变异函数值。 最后,把式(4)带入下式,待估区域的估计值可求出 Ζ*(X)=C∑c=1[Ι*c(X)][Ζ(X)|Ζc1Ζ(X)≤Ζc]*=C∑c=1[Ι*c(X)][∑bλbΖ(X)](7) 式中的权系数λb可用普通的Kriging方法求解得到。 2 实例分析 2.1 地层及构造环境 锦屏一级水电站位于雅砻江中游锦屏大河湾之西侧,地处青藏高原向四川盆地过渡的斜坡地带,地势总的趋势是西北高东南低,海拔从5000～4000m降低至约2000m,呈典型高山峡谷地貌。坝址区出露地层主要是三叠统杂谷脑组二段大理岩,厚度约600m。第一段绿片岩,在地表未出露,深埋于河床190m以下及右岸350m以里,为主要隔水层。坝址区地质构造复杂,褶皱断裂十分发育。褶皱主要分为NNE向和NWW向两组。其中NNE向三滩向斜通过坝区,其核部位于坝段左侧,为以紧闭的倒转向斜,轴面走向N20E,NW∠43°,核部地层为杂谷脑组三段六组的板岩,两翼对称出现杂谷脑三段1～5层砂板岩和二段大理岩