水 利 学 报
2003年3月
文章编号:0559 9350(2003) 03 0033 04SHUILI XUEBAO 第3期
地下水运动数值模拟过程中边界条件问题探讨
卢文喜1
(1 吉林大学环境与资源学院, 吉林长春 130026)
摘要:本文对地下水运动数值模拟过程中边界条件的涵义和处理方法进行了分析和讨论。阐述了边界条件所包含的双重意义。指出随着人类活动影响强度的日益增大, 边界条件的处理要面临一些新的更为复杂的问题。在模型预报之前必须首先对边界条件做出预报。边界条件的预报既要考虑自然因素的作用, 同时也要考虑人类活动(人工开采和人工补给) 的影响及由于邻区水流条件变化而产生的耦合效应。之后, 给出了两个应用实例。关键词:地下水; 数值模拟; 边界条件
中图分类号:P641 2文献标识码:A
在地下水运动数值模拟的过程中, 模拟预报结果的正确与否与边界条件处理得是否恰当密切相关[1,2]。尤其是在人类活动影响强度日益增大的今天, 在处理边界条件时, 常常会面临一些新的更为
[3, 4]复杂的问题。原因在于边界处的水流状况往往不仅受到自然因素的控制, 而且还深受人类活动(如人工开采和人工补给) 的影响
[5], 同时还可能受到邻区水流条件变化的扰动, 而对于人为边界更是如此。所以必须对边界条件给予应有的重视, 深入探讨其多重的内涵并研究出切实可行的处理方法。1 边界条件涵义探讨
在地下水运动数值模拟的过程中, 一般都是在概念模型的基础上, 建立描述地下水流的数学模型, 然后再采用某种数值方法, 对模型离散并求解。对于分布参数的地下水流数学模型而言, 模型主要由两部分内容组成: 描述地下水运动规律的偏微分方程; 反映地下水模拟区域具体特征的边界条件和初始条件(若为稳定运动则没有初始条件) [6, 7]。
这里的边界条件具有两重意义:一是它与初始条件一起构成地下水流数学模型的定解条件, 用来说明具体目标系统的边界所具有的特定状态, 从而使模型的求解能够得到切合实际状况的特解。二是它描述了目标系统与其周围环境之间的相互作用关系, 即它们之间物质、能量和信息的交换。也就是说, 周围环境的变化会对目标系统产生影响, 而目标系统的变化, 也会对周围环境产生作用。
边界条件的第二重意义的重要性伴随着水资源匮乏的日益严重而逐渐凸显出来。一方面, 在过去, 由于开采规模往往比较小, 人工开采对于边界水流状态的影响还比较小, 可以忽略不计。因此, 在边界处可近似将水流视为仅受自然因素影响, 这样的边界条件处理起来比较简便。然而, 随着开采规模的扩大, 人工开采对边界处的水流已造成强烈影响, 使边界处的水流状态同时受到自然、人为双重因素的制约。另一方面, 在过去, 往往周围邻区尚不存在其它的开采水源地。而现今面临的情况是:在相邻区域, 在目标系统的上游或下游存在其它开采水源地, 相邻水源地常常各自划定边界并分别建立模型, 它们在边界处的耦合效应要通过边界条件给予描述
收稿日期:2001 11 14
作者简介:卢文喜(1956-) , 男, 吉林德惠人, 教授, 博士生导师, 主要从事生态水文和地下水系统数值模拟和优化管理方面
研究。[8,9]。
在地下水运动数值模拟过程中, 边界条件的处理可划分为两个阶段。第一个阶段即在模型识别检验阶段, 边界条件利用以往的观测资料即可确定, 处理起来较为容易。第二阶段即在模型预报阶段, 边界条件处理起来较为复杂。这里存在首先要对边界条件做出预报的问题, 然后才能进行模型预报。因为在模型预报时, 边界条件是作为已知条件出现的。边界条件的预报既要考虑本区和邻区未来自然条件的变化, 同时也要考虑本区和邻区未来人为影响(如人工开采和人工补给) 的变化, 即边界条件的预报要综合考虑各种作用(包括自然因素、人为因素及相邻系统的扰动等) 的复合效应。否则, 边界条件就不能刻画边界处水流状态的真实变化, 并将最终导致整个模型计算结果的不正确。2 处理方法实例
下面给出两种处理边界条件的方法和实例, 它们分别都刻画出了自然因素和人为作用对于边界水流的复合影响。
2 1 北京市平谷电厂水源地数值模拟一类边界水位的预报 在该水源地的集中人工开采没有上马投产之前, 边界处的水位只受降水、地质条件等自然因素控制, 通过建立回归方程
h 末= 0+ 1P + 2h 初
来表达某时段的末刻水位h 末与该时段内的降水量P 和初始水位h 初之间的数量关系。其中 0为常数, 1、 2分别为回归系数。这里h 末只受自然因素控制。在该水源地的集中人工开采投产之后, 集中人工开采的影响, 这里主要是指由于大量人工抽水而引起的水位降落漏斗, 要扩展到离抽水井较近的一类边界所在的位置, 其最近距离为r =10km, 按设计开采量Q =80000m /d, 可根据Theis 公式近似算出离开采井最近的一类边界处的水位降深变化状况如表1所示。
由此可以看出, 一类边界处的水位不仅受自然因素控制, 还要受到集中人工开采的影响, 是二者叠加复合作用的结果。在预报阶段, 在模型预报之前首先对边界条件做出预报。作者采用的方法是根据预报阶段的设计开采量, 运用Theis 公式近似算出集中人工开采在一类边界上引起的水位降深S , 然后运用下式计算某时段末刻水位h 末(自然影响) 与水位降深S (人为影响) 的代数和(二者叠加) , 从而确定H 末:
H 末=h 末-S
H 末即为在h 末的基础上叠加了人工影响后该时段的末刻水位, 该水位综合反映了自然因素和集中人工开采对一类边界水位的共同影响。
2 2 某水源地数值模拟二类边界单宽流量的预报 该水源地的含水层为我国北方中、下奥陶统的碳酸岩岩层。南部山区灰岩裸露, 地势较高, 接受降水补给, 为岩溶潜水。北部丘陵区地势较低, 灰岩含水层为石炭系和二叠系砂页岩岩层所覆盖, 向北埋深逐渐增大, 为岩溶承压水。地下水总体流向为由南向北。在划定模拟区域时, 在北部边界上, 由于承压岩溶含水层深伏地下, 根据实际情况将其确定为二类弱透水的排泄边界。同样在该水源地没有上马抽水之前, 二类边界处的单宽流量q 只受降水、地质条件等自然因素控制。而在水源地上马抽水之后, 由于人工抽水而产生的水位降落漏斗会波及到离抽水井较近的二类边界所在的位置, 其最近的距离为r =1375m, 按设计开采量Q =10000m d, 根据Theis 公式近似算出离开采井最近的二类边界处的水位降深变化状况如表2所示。
表1 距离开采井最近一类边界处的降深变化
抽水时间/d
水位降深/m600. 0051000. 0183650. 1187300. 20333表2 距离开采井最近的二类边界处的降深变化抽水时间/d 水位降深/m 300 445600 5001000 5413650 645可见, 此处边界上水位的实际变化, 以及由此诱导出来的水力梯度I 和单宽流量q 的实际变化, 也是自然因素与人为作用复合影响的结果。在模型预报之前, 也必须首先要对边界条件进行预报。对该水源地北部的二类边界而言, 就是要对单宽流量q 做出预报。考虑到在北部边界处承压含水层深
伏地下, 可用资料很少, 作者采用迭代逼近方法来预报北部二类边界的单宽流量。
由达西定律, 单宽流量q =T ∀I , 其中T 为导水系数, I 为水力梯度。由于北部二类边界位于承压区, 边界处的承压含水层厚度M 和渗透系数K 均为定值, 故T =KM 可视为常数。因此, 单宽流量仅随水力梯度I 的变化而变化。无论是天然降水还是集中人工开采状况的变化, 它们对北部二类边界单宽流量的作用都最终表现为对边界处水力梯度的改变。所以, 在预报时只要抓住水力梯度这个关键因素即可。
具体步骤是:对于第i 个区段, 先给出水力梯度的初估值I i , 由q i =T I i 算出相应的单宽流量q i , 据此运转预报模型, 并根据预报结果(水位值) 计算出北部二类边界附近的水力梯度, 记为I i , 然后再运用I i 根据公式q i =T ∀I i 算出q i , 并再运转预报模型。重复进行上述迭代过程, 直到前后相邻两次迭代中的单宽流量的差值的绝对值小于预先给定的允许误差 为止, 即要满足条件:
|q i -q i
标识迭代的次序数, i 为边界的区段标号。(k ) (k -1) (2) (2) (2) (2) (2) (1) (1) (1) (1)|# 才停止迭代。然后, 取最后一次迭代过程中的单宽流量即为所求。这里 为预先给定的允许误差, k
3 分析与结论
地下水运动数值模拟预报结果的精度与边界条件处理得是否正确密切相关。边界条件的确定不仅要考虑自然因素, 还要考虑人类活动(如水源地总采量和开采布局) 及邻区水流条件变化对边界水流状态产生的复合影响。在地下水资源数值模拟过程中, 在模型预报之前, 应首先对边界条件做出预报。边界条件的预报要综合考虑多种作用的耦合效应。当资料较少时, 迭代逼近方法亦是一种可取的方法。
参 考 文 献:
[1] Bruce hunt. Mathematical Analysis of Ground water Resources [M]. Cambrige Universi ty Press, 1983.
[2] 孙纳正. 地下水流的数学模型和数值方法[M]. 北京:地质出版社, 1991.
[3] 薛禹群, 吴吉春. 地下水数值模拟在我国! ! ! 回顾与展望[J]. 水文地质工程地质, 1997, (4) :21-24.
[4] 方淑镇, 陈崇希. 文英湖地区湖周地下水排水优化管理模型[J]. 地球科学, 1993, 18(2) :211-220.
[5] 卢文喜. 响应矩阵法在实际应用中存在的问题[J].水文地质工程地质, 1994, (2) :52-53.
[6] 卢文喜. 地下水系统的模拟预测和优化管理[M]. 北京:科学出版社, 1999.
[7] 张蔚榛. 地下水与土壤水动力学[M]. 北京:中国水利水电出版社, 1996.
[8] Vermulst J A, P H, de Lange W J. An analytic based approach for coupling models for unsaturated and saturated
ground water flow at different scales [J].Journal of Hydrology, 1999, 226:262-273.
[9] De Lange W J. A Cauchy boundary condition for the lumped interacti on between an arbi trary number of surface waters
and a reg i onal aquifer [J]. Journal of Hydrology, 1999, 226:250-261.
Approach on boundary condition in numerical simu lation of groundwater flows
LU Wen xi 1
(1 Sc hool o f Environme nt and Re sources , Jilin Unive rsit y , Changc hun 130026, China )
Abstract :The implication and disposal methods of the boundary condition in the numerical simulation of groundwater flows was discussed in this paper. The double meanings that the boundary c ondition possess
es were elaborated. The boundary condition is being faced with some ne w and more complicated situations due to the intense growth of mankind activities. The boundary c ondition must be predic ted before the model simulation. In the boundary prediction, the natural factors, mankind activities and the c oupling ef fect induced by the water flow variation in neighbor areas should all be considered. Two actual examples were studied by applying the proposed approach.
Key words :groundwater flow; numerical simulation; boundary c ondition
(上接第32页)
参 考 文 献:
[1] 谢新宇, 朱向荣, 曾国熙. 饱和土体一维大变形固结理论新进展[J].岩土工程学报, 1997, 19(4) :30
-38.
[2] 钱家欢, 殷宗泽. 土工原理与计算[M]. 北京:中国水利水电出版社, 1996.
[3] 顾晓鲁, 钱鸿缙, 等. 地基与基础[M].北京:中国建筑工业出版社, 1993.
[4] 龚晓南, 沉降浅议[J].地基处理, 1996, 7(1) :41-41.
[5] 殷宗泽, 土体沉降与固结[M]. 北京:中国电力出版社, 1998.
[6] 宰金珉, 梅国维. 全过程的沉降量预测方法研究[J].岩土力学, 2000, 21(4) :322-325.
[7] 邓聚龙. 灰色系统理论教程[M]. 武汉:华中理工大学出版社, 1990.
[8] 刘思峰, 郭天榜, 党耀国, 等. 灰色系统理论及其应用[M]. 北京:科学出版社, 1999.
[9] 宰金珉, 宰金璋. 高层建筑基础分析与设计[M]. 北京:中国建筑工业出版社, 1993.
Grey Verhulst prediction of settlement during the whole
process of construction and operation
LUO Zhan you , GONG Xiao nan , YANG Xiao jun 111
(1.Geotec hnic al Enginee ring Instit ute , Zhejiang U nive rsity , Hangzhou 310027, China )
Abstract :The authors proposed a gray Verhulst forecasting model based on the gray system theory and the settlement time relationship during the whole process of construction and operation. The applicability of this model was discussed and the model was applied to one practical project. B y analyzing the result the model was verified. The model can provide an effective way to predic t the settlement that would hap pen after the c onstruction was completed.
Key words :gray Verhulst forecasting model; settlement time relationship; gray system theory; Verhulst theory
水 利 学 报
2003年3月
文章编号:0559 9350(2003) 03 0033 04SHUILI XUEBAO 第3期
地下水运动数值模拟过程中边界条件问题探讨
卢文喜1
(1 吉林大学环境与资源学院, 吉林长春 130026)
摘要:本文对地下水运动数值模拟过程中边界条件的涵义和处理方法进行了分析和讨论。阐述了边界条件所包含的双重意义。指出随着人类活动影响强度的日益增大, 边界条件的处理要面临一些新的更为复杂的问题。在模型预报之前必须首先对边界条件做出预报。边界条件的预报既要考虑自然因素的作用, 同时也要考虑人类活动(人工开采和人工补给) 的影响及由于邻区水流条件变化而产生的耦合效应。之后, 给出了两个应用实例。关键词:地下水; 数值模拟; 边界条件
中图分类号:P641 2文献标识码:A
在地下水运动数值模拟的过程中, 模拟预报结果的正确与否与边界条件处理得是否恰当密切相关[1,2]。尤其是在人类活动影响强度日益增大的今天, 在处理边界条件时, 常常会面临一些新的更为
[3, 4]复杂的问题。原因在于边界处的水流状况往往不仅受到自然因素的控制, 而且还深受人类活动(如人工开采和人工补给) 的影响
[5], 同时还可能受到邻区水流条件变化的扰动, 而对于人为边界更是如此。所以必须对边界条件给予应有的重视, 深入探讨其多重的内涵并研究出切实可行的处理方法。1 边界条件涵义探讨
在地下水运动数值模拟的过程中, 一般都是在概念模型的基础上, 建立描述地下水流的数学模型, 然后再采用某种数值方法, 对模型离散并求解。对于分布参数的地下水流数学模型而言, 模型主要由两部分内容组成: 描述地下水运动规律的偏微分方程; 反映地下水模拟区域具体特征的边界条件和初始条件(若为稳定运动则没有初始条件) [6, 7]。
这里的边界条件具有两重意义:一是它与初始条件一起构成地下水流数学模型的定解条件, 用来说明具体目标系统的边界所具有的特定状态, 从而使模型的求解能够得到切合实际状况的特解。二是它描述了目标系统与其周围环境之间的相互作用关系, 即它们之间物质、能量和信息的交换。也就是说, 周围环境的变化会对目标系统产生影响, 而目标系统的变化, 也会对周围环境产生作用。
边界条件的第二重意义的重要性伴随着水资源匮乏的日益严重而逐渐凸显出来。一方面, 在过去, 由于开采规模往往比较小, 人工开采对于边界水流状态的影响还比较小, 可以忽略不计。因此, 在边界处可近似将水流视为仅受自然因素影响, 这样的边界条件处理起来比较简便。然而, 随着开采规模的扩大, 人工开采对边界处的水流已造成强烈影响, 使边界处的水流状态同时受到自然、人为双重因素的制约。另一方面, 在过去, 往往周围邻区尚不存在其它的开采水源地。而现今面临的情况是:在相邻区域, 在目标系统的上游或下游存在其它开采水源地, 相邻水源地常常各自划定边界并分别建立模型, 它们在边界处的耦合效应要通过边界条件给予描述
收稿日期:2001 11 14
作者简介:卢文喜(1956-) , 男, 吉林德惠人, 教授, 博士生导师, 主要从事生态水文和地下水系统数值模拟和优化管理方面
研究。[8,9]。
在地下水运动数值模拟过程中, 边界条件的处理可划分为两个阶段。第一个阶段即在模型识别检验阶段, 边界条件利用以往的观测资料即可确定, 处理起来较为容易。第二阶段即在模型预报阶段, 边界条件处理起来较为复杂。这里存在首先要对边界条件做出预报的问题, 然后才能进行模型预报。因为在模型预报时, 边界条件是作为已知条件出现的。边界条件的预报既要考虑本区和邻区未来自然条件的变化, 同时也要考虑本区和邻区未来人为影响(如人工开采和人工补给) 的变化, 即边界条件的预报要综合考虑各种作用(包括自然因素、人为因素及相邻系统的扰动等) 的复合效应。否则, 边界条件就不能刻画边界处水流状态的真实变化, 并将最终导致整个模型计算结果的不正确。2 处理方法实例
下面给出两种处理边界条件的方法和实例, 它们分别都刻画出了自然因素和人为作用对于边界水流的复合影响。
2 1 北京市平谷电厂水源地数值模拟一类边界水位的预报 在该水源地的集中人工开采没有上马投产之前, 边界处的水位只受降水、地质条件等自然因素控制, 通过建立回归方程
h 末= 0+ 1P + 2h 初
来表达某时段的末刻水位h 末与该时段内的降水量P 和初始水位h 初之间的数量关系。其中 0为常数, 1、 2分别为回归系数。这里h 末只受自然因素控制。在该水源地的集中人工开采投产之后, 集中人工开采的影响, 这里主要是指由于大量人工抽水而引起的水位降落漏斗, 要扩展到离抽水井较近的一类边界所在的位置, 其最近距离为r =10km, 按设计开采量Q =80000m /d, 可根据Theis 公式近似算出离开采井最近的一类边界处的水位降深变化状况如表1所示。
由此可以看出, 一类边界处的水位不仅受自然因素控制, 还要受到集中人工开采的影响, 是二者叠加复合作用的结果。在预报阶段, 在模型预报之前首先对边界条件做出预报。作者采用的方法是根据预报阶段的设计开采量, 运用Theis 公式近似算出集中人工开采在一类边界上引起的水位降深S , 然后运用下式计算某时段末刻水位h 末(自然影响) 与水位降深S (人为影响) 的代数和(二者叠加) , 从而确定H 末:
H 末=h 末-S
H 末即为在h 末的基础上叠加了人工影响后该时段的末刻水位, 该水位综合反映了自然因素和集中人工开采对一类边界水位的共同影响。
2 2 某水源地数值模拟二类边界单宽流量的预报 该水源地的含水层为我国北方中、下奥陶统的碳酸岩岩层。南部山区灰岩裸露, 地势较高, 接受降水补给, 为岩溶潜水。北部丘陵区地势较低, 灰岩含水层为石炭系和二叠系砂页岩岩层所覆盖, 向北埋深逐渐增大, 为岩溶承压水。地下水总体流向为由南向北。在划定模拟区域时, 在北部边界上, 由于承压岩溶含水层深伏地下, 根据实际情况将其确定为二类弱透水的排泄边界。同样在该水源地没有上马抽水之前, 二类边界处的单宽流量q 只受降水、地质条件等自然因素控制。而在水源地上马抽水之后, 由于人工抽水而产生的水位降落漏斗会波及到离抽水井较近的二类边界所在的位置, 其最近的距离为r =1375m, 按设计开采量Q =10000m d, 根据Theis 公式近似算出离开采井最近的二类边界处的水位降深变化状况如表2所示。
表1 距离开采井最近一类边界处的降深变化
抽水时间/d
水位降深/m600. 0051000. 0183650. 1187300. 20333表2 距离开采井最近的二类边界处的降深变化抽水时间/d 水位降深/m 300 445600 5001000 5413650 645可见, 此处边界上水位的实际变化, 以及由此诱导出来的水力梯度I 和单宽流量q 的实际变化, 也是自然因素与人为作用复合影响的结果。在模型预报之前, 也必须首先要对边界条件进行预报。对该水源地北部的二类边界而言, 就是要对单宽流量q 做出预报。考虑到在北部边界处承压含水层深
伏地下, 可用资料很少, 作者采用迭代逼近方法来预报北部二类边界的单宽流量。
由达西定律, 单宽流量q =T ∀I , 其中T 为导水系数, I 为水力梯度。由于北部二类边界位于承压区, 边界处的承压含水层厚度M 和渗透系数K 均为定值, 故T =KM 可视为常数。因此, 单宽流量仅随水力梯度I 的变化而变化。无论是天然降水还是集中人工开采状况的变化, 它们对北部二类边界单宽流量的作用都最终表现为对边界处水力梯度的改变。所以, 在预报时只要抓住水力梯度这个关键因素即可。
具体步骤是:对于第i 个区段, 先给出水力梯度的初估值I i , 由q i =T I i 算出相应的单宽流量q i , 据此运转预报模型, 并根据预报结果(水位值) 计算出北部二类边界附近的水力梯度, 记为I i , 然后再运用I i 根据公式q i =T ∀I i 算出q i , 并再运转预报模型。重复进行上述迭代过程, 直到前后相邻两次迭代中的单宽流量的差值的绝对值小于预先给定的允许误差 为止, 即要满足条件:
|q i -q i
标识迭代的次序数, i 为边界的区段标号。(k ) (k -1) (2) (2) (2) (2) (2) (1) (1) (1) (1)|# 才停止迭代。然后, 取最后一次迭代过程中的单宽流量即为所求。这里 为预先给定的允许误差, k
3 分析与结论
地下水运动数值模拟预报结果的精度与边界条件处理得是否正确密切相关。边界条件的确定不仅要考虑自然因素, 还要考虑人类活动(如水源地总采量和开采布局) 及邻区水流条件变化对边界水流状态产生的复合影响。在地下水资源数值模拟过程中, 在模型预报之前, 应首先对边界条件做出预报。边界条件的预报要综合考虑多种作用的耦合效应。当资料较少时, 迭代逼近方法亦是一种可取的方法。
参 考 文 献:
[1] Bruce hunt. Mathematical Analysis of Ground water Resources [M]. Cambrige Universi ty Press, 1983.
[2] 孙纳正. 地下水流的数学模型和数值方法[M]. 北京:地质出版社, 1991.
[3] 薛禹群, 吴吉春. 地下水数值模拟在我国! ! ! 回顾与展望[J]. 水文地质工程地质, 1997, (4) :21-24.
[4] 方淑镇, 陈崇希. 文英湖地区湖周地下水排水优化管理模型[J]. 地球科学, 1993, 18(2) :211-220.
[5] 卢文喜. 响应矩阵法在实际应用中存在的问题[J].水文地质工程地质, 1994, (2) :52-53.
[6] 卢文喜. 地下水系统的模拟预测和优化管理[M]. 北京:科学出版社, 1999.
[7] 张蔚榛. 地下水与土壤水动力学[M]. 北京:中国水利水电出版社, 1996.
[8] Vermulst J A, P H, de Lange W J. An analytic based approach for coupling models for unsaturated and saturated
ground water flow at different scales [J].Journal of Hydrology, 1999, 226:262-273.
[9] De Lange W J. A Cauchy boundary condition for the lumped interacti on between an arbi trary number of surface waters
and a reg i onal aquifer [J]. Journal of Hydrology, 1999, 226:250-261.
Approach on boundary condition in numerical simu lation of groundwater flows
LU Wen xi 1
(1 Sc hool o f Environme nt and Re sources , Jilin Unive rsit y , Changc hun 130026, China )
Abstract :The implication and disposal methods of the boundary condition in the numerical simulation of groundwater flows was discussed in this paper. The double meanings that the boundary c ondition possess
es were elaborated. The boundary condition is being faced with some ne w and more complicated situations due to the intense growth of mankind activities. The boundary c ondition must be predic ted before the model simulation. In the boundary prediction, the natural factors, mankind activities and the c oupling ef fect induced by the water flow variation in neighbor areas should all be considered. Two actual examples were studied by applying the proposed approach.
Key words :groundwater flow; numerical simulation; boundary c ondition
(上接第32页)
参 考 文 献:
[1] 谢新宇, 朱向荣, 曾国熙. 饱和土体一维大变形固结理论新进展[J].岩土工程学报, 1997, 19(4) :30
-38.
[2] 钱家欢, 殷宗泽. 土工原理与计算[M]. 北京:中国水利水电出版社, 1996.
[3] 顾晓鲁, 钱鸿缙, 等. 地基与基础[M].北京:中国建筑工业出版社, 1993.
[4] 龚晓南, 沉降浅议[J].地基处理, 1996, 7(1) :41-41.
[5] 殷宗泽, 土体沉降与固结[M]. 北京:中国电力出版社, 1998.
[6] 宰金珉, 梅国维. 全过程的沉降量预测方法研究[J].岩土力学, 2000, 21(4) :322-325.
[7] 邓聚龙. 灰色系统理论教程[M]. 武汉:华中理工大学出版社, 1990.
[8] 刘思峰, 郭天榜, 党耀国, 等. 灰色系统理论及其应用[M]. 北京:科学出版社, 1999.
[9] 宰金珉, 宰金璋. 高层建筑基础分析与设计[M]. 北京:中国建筑工业出版社, 1993.
Grey Verhulst prediction of settlement during the whole
process of construction and operation
LUO Zhan you , GONG Xiao nan , YANG Xiao jun 111
(1.Geotec hnic al Enginee ring Instit ute , Zhejiang U nive rsity , Hangzhou 310027, China )
Abstract :The authors proposed a gray Verhulst forecasting model based on the gray system theory and the settlement time relationship during the whole process of construction and operation. The applicability of this model was discussed and the model was applied to one practical project. B y analyzing the result the model was verified. The model can provide an effective way to predic t the settlement that would hap pen after the c onstruction was completed.
Key words :gray Verhulst forecasting model; settlement time relationship; gray system theory; Verhulst theory