线弹性含水层井水位、孔压对引潮位响应的研究及其应用.pdfVIP

第 33卷 第 1期 地 震 地 质 Vo1．33．No．1 2011年 3月 Mar．．2011 SEISM0L0GY AND GEOLOGY 线弹性含水层井水位、孑L压 对引潮位响应的研究及其应用 刘春平 唐彦东 廖 欣 万 飞 石 云 (防灾科技学院，三河 065200) 摘 要 应用线弹性介质力均衡方程 ，研究不排水条件下饱水岩体在潮汐力作用下的体应变 ，提 出 了承压含水层孔压对引潮高的线性响应方程，并给出了该方程响应系数 ()的物理意义。结合 Hsieh 等 (1987)提出的井水位对孑L压响应的振幅比(A)和位相差 (．)公式，进一步推导出了井水位一引潮高 振幅比M=EA和位相差 =d + ，公式。 和 口是基于实测井水位和理论引潮高数据计算的，在孔 压与引潮高位相差 (：)已知的条件下，由 和 n，可依次计算含水层导水系数 ()、水位一孔压振幅 比(A)，以及孔压一引潮高振幅比(E)。选择JlI18、川06井作为实例 ，计算和分析 了T、A、E和 的变 化 。 关键词 含水层 孔弹性介质 引潮高 孔压 中图分类号：P312．4 文献标识码 ：A 文章编号 ：0253—4967(2011)Ol一0133—08 0 引言 潮汐力、气压、土壤水分和雪荷载等高频荷载的作用周期从 日到季，在地质意义上是非常短 的。由于这些荷载的重复性和小量级 ，它们所引起的地球形变非常接近于线弹性 ，因此 ，这些荷 载所导致 的水压变化可以由线性孔弹性理论精确描述 (Wang，2000)。Biot(1962)构建 了孔压 和介质形变耦合孔弹性理论 ，并广泛用于解释含水层孔压与应力一应变关系 (Bodvarsson，1970； Geeta1．，1992；Neuzil，2003)。潮汐力作为一个 已知的小量级高频荷载 ，所 引起 的井水位潮汐 现象揭示 了含水层孔压与岩石形变的耦合关系，一直受到岩石力学和地震学界的关注。他们应 用 Biot孑L弹性理论建立了孔压 、井水位与应力 (或应变)方程 ，研究了井水位对潮汐力的响应函 数 (Bredehoeft，1967；vanderKampeta1．，1983；张昭栋等，1991，1995)，并通过响应(或固体潮) 系数研究含水层介质特征 ，估计水动力学参数 (Hsieheta1．，1987；张昭栋等 ，1999；2002)。 Hsieh等提出了水位和含水层孔压振幅比(A)计算公式，从理论上建立 了井水位一孔压的响应关 系。本文运用孔弹性应力平衡理论推导 出了含水层孔压与引潮位关系式 ，结合 Hsieh等 (1987) 给出的井水位与孔压之间的关系 ，分析了井水位 、含水层孔压和引潮位三者之间的振幅 比和位 相差 。 l 线弹性含水层孔压随引潮位变化方程 高频荷载周期相对于排水或水力响应发生的时间尺度很短，可应用不排水条件 的应力一应 [收稿 日期] 2010—10—10收稿 ，2011—03—04改回。 【基金项 目) 地震行业科研专项 (200808055，200808079)资助。 134 地 震 地 质 变和流体压力式。 = (A+卢N)6，1+2x／e (1) P=一 Ⅳ， (2) 式(1)、(2)中or是应力 ；P是流体压力 ；A， 为多孔介质固体骨架的拉梅参数 ；N为饱和地质 岩体的Biot模量，是与流／固组合的压缩性有关的弹性参数； 为 Biot系数，又称为有效应力系 数；6提 Kronecker函数；s是岩石应变张量；，1=s1+ 22+833是岩石体应变。 如果不考虑岩体运动加速度 ，在 (i=1，2，3表示空问3个坐标轴)方向上的应力分量总和