碾压混凝土坝等效黏弹性参数反演.pdf

第 3 5卷第 8期 2 0 1 3年 8月 人民黄河 YELLOW RI VER Vo 1 3 5 No 8 Au g ， 2 0 1 3 【 水利水电工程 】 碾压混凝土坝等效黏弹性参数反演 李姝昱 ， 白家泽 ， 苏 琼 ( 1 黄河水利科学研究院， 河南 郑州4 5 0 0 0 3 ； 2 河南东龙控股有限公司， 河南 郑州 4 5 0 0 0 0 ； 3 河海大学 水利水电学院， 江苏 南京 2 1 0 0 9 8 ) 摘要： 根据碾压混凝土坝的成层结构特性， 将层面影响概化到整个坝体中， 将碾压混凝土坝视为横观各向同性体， 建立 了碾压混凝土坝等效黏弹性模型。根据位移实测资料， 采用偏最小二乘法建立位移统计模型， 分 离水压分量和时效分 量， 建立反演碾压混凝土坝双向异性黏弹性参数的目标函数 ， 并利用混沌遗传算法对黏弹性参数进行反演。工程实例 中， 将反演得到的黏弹性参数作为已知值， 计算测点在相应时间段变形的水压和时效分量之和， 并与分离出的实测值进 行对比， 结果表明： 计算值与实测值比较接近 ， 说明反演方法是正确、 可行的。 关键词： 碾压混凝土坝；黏弹性参数；混沌遗传算法；偏最小二乘法 中图分类号 ：T V 6 4 2 2 文献标 志码 ： A d o i ： 1 0 3 9 6 9 j i s s n 1 0 0 0 1 3 7 9 2 0 1 3 0 8 0 3 2 I n v e r s i o n o f t h e Eq u i v a l e nt Vi s c o e l a s t i c Pa r a me t e r s o f t he Ro l l e r Co mpa c t e d Co n c r e t e Da m L I S h u y u ，B a i J i a z e ，S U Q i o n g 3 ( 1 Y e l l o w Riv e r I n s t i t u t e o f H y d r a u l i c R e s e a r c h ，Z h e n g z h o u 4 5 0 0 0 3， C h i n a ；2 H e n a n D o n g l o n g Ho l d i n g s C o L t d， Z h e n g z h o u 4 5 0 0 0 0， C h i n a ； 3 C o l l e g e o f Wa t e r C o n s e r v a n c y a n d H y d rop o w e r ， H o h a i U n i v e r s i t y ，N a n j i n g 2 1 0 0 9 8 ，C h i n a ) Ab s t r a c t ：A c c o r d i n g t o t h e l a y e r f e a t u r e o f r o l l e r c o mp a c t e d c o n c r e t e d a m ( R C C D) ， t h i s p a p e r g e n e r a l i z e d t h e l a y e r i n fl u e n c e t o t h e e n t i re d a m， a n d t o o k rol l e r c o mp a c t e d c o n c r e t e d a m a s t r a n s v e r s e i s o t ropy ，t h e n b u i l t t h e e q u i v ale n t v i s c o e l a s t i c mo d e 1 Ac c o r d i n g t o t h e me asu r e d d a t a ，t h e d a m d i s p l a c e me n t s t a t i s t i c a l mo d e l wa s e s t a b l i s h e d t o s e p a r a t e t h e wa t e r p r e s s u re c o mp o n e n t a n d t h e a g i n g c o mp o n e n t u s i n g t he p a r t i al l e a s t s q u are s r e g r e s s i o n me t h odA n d t h e o b j e c t i v e f u n c t i o n w a s e s t a b l i s h e d t o i n v e rs e t h e b i a x i al a n i s o t rop i c v i s c o e l a s t i c p aram e t e r s o f rol l e r c o mp a c t e d c o n c r e t e d a m u t i l i z i n g c h a o s g e n e t i c alg o rit h mT h rou g h a n e n g i n e e rin g e x a mp l e ，t h e i n v e r s e d v i s c oel ast i c p aram e t e r s we r e s e t as k n o wn v alu e s ；t h e s u m o f t h e w a t e r p r e s s u re c o mp o n e n t a n d t h e a g i n g c o mp o n e n t wa s c alc u l a t e d a n d wa s c o mp a r e d wi t h t h e me asu red v alu e T h e res u l t s s h o w t h a t t h e c alc u - l a t e d v alu e a p p rox i ma t e s t o t h e o b j e c t i v e v a l u e a n d p rov e t h a t t h i s i n v e rsi o n m e tho d i s f e asib l e Ke y wo r d s ：mi l e r c o mp a c t e d c o n c r e t e d a m ；v i s c o e l a s t i c p a r am e t e rs；c h a o s g e n e t i c a l g o rithm ；p a r t i al l e as t s q u a r e s me tho d 碾压混凝土筑坝技术始于2 O世纪 7 O年代末， 因其具有工 艺简单 、 工期短、 造价低、 适应性强等特点而在大坝建设中得 到广泛应用。不同于常规混凝土坝， 碾压混凝土坝有几十个甚 至几百个施工层面， 其层面的力学性能表明： 层面是碾压混凝 土坝的薄弱环节， 各项强度指标都小于本体， 是影响碾压混凝 土坝强度、 稳定和渗流的关键部位 J 。因此， 要对碾压混凝土 坝运行时的工作状态进行分析就需要正确模拟层面结构。 对碾压混凝土坝层面结构进行模拟有两种思路： 一种是直 接模拟层面， 如采用节理单元 j 、 裂缝单元 J 、 薄层单元 J 等， 但碾压混凝土坝施工层面多， 直接模拟层面计算量大， 也给 反演层面力学参数带来了困难 ； 另一种是将层面影响概化到整 个坝体中， 在建立碾压混凝土坝整体连续模型的同时模拟本体 和层面的力学特性。笔者采取第二种思路 ， 将碾压混凝土坝视 为横观各向同性体 ， 利用实测资料反演出垂直和平行于层面方 向的等效黏弹性参数。 1 0 2 1 碾压混凝土坝等效黏弹性本构模型 1 1 等效弹性本构模型 在逐层碾压过程中， 碾压混凝土坝的层与层之间形成过渡 区， 即层面影响带。假设碾压混凝土沿高程方向每层厚度为 h ， 本体厚度为 h ， 层面厚度为 h ， 则 h=h +h 。碾压混凝土平 行于层面方向的等效弹性模量、 泊松比和等效切变模量分别为 巨 、 G ， 垂直层面方向的等效弹性模量、 泊松比和等效切变模 量分别为 E , t 、 G ， 其中E 。 l 、 E 、 G 为 5 个相互独立的参 数。等效弹性模型对应的应力应变关系为 收稿日 期 ： 加1 2 1 2 2 4 基金项 目 ： 水利部科技推 广计 划项 目( T G 1 1 1 6 ) ； 水利部 公益性 行 业科研 专项 ( 2 0 1 3 0 1 0 6 1 ) ； 水利标准化工作资助项 目( 1 2 6 1 1 2 0 1 6 2 5 0 9 ， 1 2 6 1 1 2 0 1 6 2 5 1 0 ) ； 水利 重大政策研究和制度建设专项资助项目( 1 2 6 1 3 2 0 1 6 2 5 1 4 ) 。 作者简介： 李妹昱( 1 9 8 8 一) ， 女， 河南濮阳人， 助理工程师， 硕士， 主要从事大坝 安全监测工作 E - ma i l ：h h u l s y 1 6 3 c o m 人 民 黄 河2 0 1 3 年第 8期 =D e ( 1 ) 式中： 为应力矩阵； 为应变矩阵； D 为横观各向同性弹性 矩 阵。 D = E。 E E 等 告 E， 2 ( 1+ 1 ) G ( 2 ) 式中： g= E 1 E ； p：1一 l 一 2 q l x 。 1 2 等效黏弹性本构模型 采用三元件模型模拟碾压混凝土坝的黏弹性特性， 碾压混 凝土总应变由弹性应变和复合黏性应变两部分组成 ， 即其总应 变 为 s = + ( 3 ) 式中： 为总应变 ； 8 为弹性应变 ； e 为复合黏性应变。 弹性应变满足 H o o k e 定律， 弹性应变为 E = 1 ( 4 ) 式中： 为应力矩阵； D 为横观各向同性弹性矩阵， 与碾压混凝 土瞬时弹性模量 。 ( q= 口 ， z ) 有关。 复合黏性应变由三元件固体模型确定： ， 7 。 s +E 2 8 =o ( 5 ) 式中： 为复合黏性应变对时间的导数； 。 ( q= ， z ) 为碾压混 凝土的延迟弹性模量； 仉( q= ， Z ) 为碾压混凝土黏性系数。 对式 ( 5 ) 积分得 s = + C 爰 e x p ( 一 d c 【 e x p 一 ( ( 6 ) 式中： 为 t 。时刻结构的初始黏性应变 ； D：为弹性矩阵， 与碾 压混凝士延迟弹性模量 ( g= ， f ) 有关。 考虑 为常应力的情况得 s =和 一 ( 0 ) 1 ) 2 1 一 x p 一 ( 7 ) 2 黏弹性参数优化反演分析 将碾压混凝土视为横观各向同性体进行三元件的黏弹性 非线性有限元分析， 需要反演出平行于层面方向及垂直于层面 方向的瞬时弹性模量和延迟弹性模量， 共 6个参数 ， 分别为瞬 时弹性模量 E l 1 、 延迟弹性模量 。 、 和黏性系数 町 。 、 。 碾压混凝土坝的位移是各种因素影响下大坝稳定性的综合反 映， 在建立碾压混凝土坝等效黏弹性模型的基础上 ， 利用坝体 的位移值对坝段的等效力学参数进行反演。碾压混凝土坝的 位移主要 由水压分量、 温度分量和时效分量组成 j ， 位移为 6=6 H+6 T+6 0 ( 8 ) 式中： 艿为碾压混凝土坝位移； 占 为位移水压分量； 6 为位移温 度分量； 6 。 为位移时效分量。 反演瞬时弹性模量 E E 时， 需要通过数学模型从总位移 中分离出水压分量 ， 即6 = 6一 研 一 氐， 并将不同时刻的6 作为 目标进行反分析。同理 ， 通过建立位移数学模型分离出水压和 时效分量之和 6 + 。 = 6一 砩， 将不同时刻的6 + 。 作为 目标进行 反演。在之前反演出的E E。 基础上， 进一步反演得到碾压混 凝土坝的延迟弹性模量 。 、 和黏性系数 叼 ， 、 7 。 使用优化反分析法反演坝体的黏弹性参数时还需要确定 目标函数 ， 通常将有限元计算出的各位移分量与实测各位移分 量的残差加权平方和作为目标函数， 即目标函数为 m l 1 F ( ) = ( 一 6 ) X c n ( 9 ) 式中： m 。 为测点数； n 为测点 i 的监测次数； 伪 权重， 其值取 决于测点的重要性和测量精度， 精度高的点相应的权重值大； 氏 为在力学参数 =( 。 ， ： ， ， ) 下， 通过等效黏弹性模型 计算出的位移分量 ； 6 为通过建立位移统计模型分离出的实测 值分量； 力为力学参数 的取值范围， 通常依据工程经验及现 场试验给出， 即 口 b 。 ( 1 i n ) ， 0 、 b 分别为第 i 个力学 参数 的最小值 、 最大值。 确定了目标函数后 ， 参数的反演问题实际上就转变为求解 最优化问题， 就是要在满足约束条件的前提下， 求出能够使优 化反演 目标 函数 达到极 小值 的参数组合 ， 即求 使得 F= F ( x) 一mi n满足约束条件 n 6 。通常优化反分析法不能 够直接求解， 而是要通过正分析法进行反复试算 ， 最终得到参 数范围内的最优解。求解过程中， 要选取合适的方法对参数进 行优化， 笔者采用改进的混沌遗传算法进行优化求解。 2 1 混沌遗传算法 遗传算法是一类可用于解决复杂非线性优化问题的随机 搜索方法， 但在优化求解过程中也存在一些不足 ： 缺乏产生最 优个体的强大能力 ， 搜索在接近全局最优解时速度变慢， 不易 收敛， 易陷入局部最优 ； 存在振荡、 随机性大等问题。为弥补遗 传算法的不足， 通常将遗传算法与其他算法集成起来。混沌足 自然界普遍存在的一种非线性现象 ， 充分体现了系统的复杂 性， 它具有随机性 、 规律性和遍历性 。根据这些特性， 将混沌 优化算法和遗传算法集成起来， 可充分发挥两种算法的优点： 可以利用混沌算法在较小参数空间范围内较强的搜索能力和 遗传算法的全局搜索能力； 同时， 弥补了混沌优化算法搜索范 围小和遗传算法收敛速度慢或未成熟收敛的不足。 2 1 1混 沌 优 化 方 法 混沌优化过程主要分为两个阶段： 首先， 在变量取值范围 内按 自身规律不重复地遍历各点， 根据适应度值的大小挑选出 最优点； 然后， 以当前确定的最优点为第二阶段的搜索起始点， 通过附加小幅度的扰动， 继续开展局部区域内的细搜索 ， 当满 足算法终止准则时停止搜索 J 。具体步骤如下。 ( 1 ) 初始化， 采用 L o s t i c映射， 赋予式( 1 0 ) n个不同的初 值( 不能取方程的不动点 0 、 0 7 5 ) ， 其中 n为待优化参数个数， 得到 n个不同轨迹的混沌变量 t ( k ) ( k：1 ， 2 ， 其中 s 为混 】 0 3 人 民 黄 河2 0 1 3年第 8期 沌迭代的次数) 。 t ( 后+1 )=g t ( ) 1一t ( Jj ) 0t ( ) 1 ( 1 O ) 式中 为控制参数， = 4时， 为典型的 L o g i s t i c映射， 此时系统 达到完全混沌 。 ( 2 ) 通过式( 1 1 ) ， 用载波方法将混沌变量的取值范围变换 到待反演参数的取值范围。 ( k )=a +( b 一0 ) t ( k ) ( 1 1 ) ( 3 ) 置 k =1 ， 由 ( k ) 计算适应度值， ( k ) 。令 = ( 0 ) ， f = ， ( 0 ) 当， ( k ) ， 时 = ， ( k ) ， = ( k ) ； 当f ( k ) ， 时， 放弃 ( k ) 。 ( 4 ) 置 k =k+1 ， 重复步骤( 3 ) ， 经s 步搜索， 直到l厂 保持不 变， 即搜索到当前最优解 = ， ( ) 。 ( 5 ) 按式( 1 2 ) 进行第二次载波。 ( k )= + t ( k ) ( 1 2 ) 式中 t ( k ) 为遍历区间很小的混沌变量， 为调节常数， 可以 小 于 1 。 ( 6 ) 置 k 1 ， 对第二次载波后的变量进行迭代， 由 ( k ) 计算适应度值， ( k ) 。当， ( k ) 厂 时 = ， ( k ) ， = ( k ) ； 当， ( k ) ， 时， 放弃 ( k ) 。 ( 7 ) 置 k =k +1 ， 当满足条件时停止搜索， 得到最优解 ，否则返回步骤( 6 ) 继续搜索。 2 1 2 混沌遗传算法步骤 将遗传算法和混沌算法智能集合起来， 通常分为 3个步 骤： 第一步， 将遗传算法生成的初始种群映射到混沌区间进行 混沌迭代， 并通过逆映射将迭代个体的取值范围变换到待反演 参数的取值范围； 第二步， 根据适应度大小， 选择可作为父代的 个体进行交叉、 变异基本的遗传操作； 第三步， 对挑选出的适应 度最大的个体进行混沌优化。循环进行这 3个步骤直到满足 终止条件， 得到最优解。 2 2 偏最小二乘法 在反演碾压混凝土坝黏弹性参数时， 需要建立实测位移的 统计模型分离出实际的水压分量 6 和时效分量 占 。 ， 因此要选 取合理的建立统计模型的方法。偏最小二乘法是一种新型的 多元统计方法， 在一个算法下可以同时实现回归建模、 数据结 构简化 ( 主成分分析) 以及变量间的相关分析( 典型相关分 析) 。 2 2 1 偏最小二乘法算法流程 偏最小二乘法不直接建立因变量 y与 自变量 x的回归模 型， 而是采用循环式的信息提取和分解方法， 逐次筛选组合 自 变量， 挑选 出对因变量具有最佳解释能力的新综合变量( 成 分) ， 进而建立回归模型。 ( 1 ) 对因变量 l ， 和自变量 x分别进行标准化处理， 得到数 据矩阵 F 0 和 E n 。 ( 2 ) 从 F 0中提取一个成分 l f ， = F0 C 1 ，Il c l II=1 ； 从 E 0 中提 取一个 成 分 t ， ， t ，要 尽 可能 携带 的变 异信 息， 即 V a r ( t 1 ) 一m a 】 【 ， 且 t 与 y的相关程度最大， 即 r ( t 1 ， Y ) 一m a ) 【 。 t 。 =E o ， 。 ， 1 ， 为 的第一主轴， I ， ， ： F o l l E F o l l 。分别实 施 和 F n 在 t 上的 回归 ： 】 0 4 f E o -p E j ( 1 3 ) I Fo t 1 ， + F 式中： E。 、 F。 为残差矩阵； p 。 、 r I 为回归系数向量 = ， f 一 ll t l lI 。 ( 3 ) 用 E 、 F。 代替 E 。 、 F 。 ， 按照上述步骤， 求得 “ ： 、 t 、 P 、 r 2 ， 如此往复类推得到 t 、 t 2 、 、 t ( mA， A=秩( ) ) ， 实施 F。 关 于 t 的回归 ， 得到 F o=， 1 t l+r 2 f 2+ +r m t ( 1 4 ) 式中： t 。 、 t 2 、 、 t 均是 E 0的线性 组合； F 。可写作 E 。的表 达式 。 ( 4 ) 最后进行变换处理， 将所得方程还原为因变量 y对自 变量 X的回归方程。 2 2 2交叉有 效性 采用偏最小二乘法建立模型时， 需要确定成分个数 ， 通常 采用交叉有效性判别方法。记 Y 为原始数据； 为对所有样本 点取成分 t 。 、 t ： 、 、 t 建立回归方程后第 i 个样本点的拟合值； )为删去第 i 个样本点， 实施对 t 。 、 t ： 、 、 t 的回归后计算得 到的 Y 拟合值， 则有判别式 l一 ( 1 5 ) 式中 ： P S S =( 一 ( ) ； S S 为 全 部样本 拟合的 具 有 h 一 1 个 成 分的 回 归 方程的 拟 合 误差 平方 和， S S =( 一 夕 川) 。 当Q 0 0 9 7 5 ( 即 0 9 5 ) 时， 新引进的 成分 t 会对模型的预测能力有显著改善作用， 应选择成分 t ， 否则认为t 对模型的贡献不显著， 不应选择 t 。 3 工程实例 某大坝为碾压混凝土重力坝， 正常蓄水位为3 7 5 O 0 m， 校 核洪水位为 3 8 1 8 4 m， 死水位为3 3 0 O 0 m， 坝顶高程为 3 8 2 O 0 m。该大坝共分为3 2个坝段， 其中， 1 一 4 、 6 1 1 为右岸重力 坝段， 5 为通航坝段 ， l 2 一l 9 为溢流坝段， 2 2 一3 0 为厂房坝 段， 2 0 2 1 、 3 1 3 2 为左岸重力坝段。以右岸挡水坝段中的 1 1 坝段为例 ， 沿坝基上下游方向和坝基深度方向分别取约 1 5 倍坝高， 建立三维有限元模型， 坝体和坝基物理力学参数见表 1 。模型采用六面体 8节点等参单元和五面体 6节点单元 ， 单 元总数为1 6 1 0 个， 节点总数为 2 6 2 2个， 其中坝体 3 2 2个单元 ， 坝基 1 2 8 8个单元。 为了监测大坝的位移 ， 在该坝段布设倒垂线 I P一1 1 ， 正垂 线 P L 1 1 3 、 P L 1 122 、 P L 1 121 、 P L 1 11 ， 共 5个测点 ， 观测高程2 2 2 7 5 、 2 7 0 O 0 、 3 1 0 O 0 、 3 4 2 O 0 、 3 7 9 2 0 m相对高程 1 5 1 5 0 m的位移。结合该工程的实际情况， 反演资料系列取 P L 1 1 2 2测点 2 0 0 61 0 0 3 -2 0 0 8 0 4 0 2的位移测值。 根据 P L 1 1 2 2水平位移测点测值 ， 采用偏最小二乘法建立 水平位移统计模型分离各分量。根据工程经验及相关试验资 人 民 黄 河2 0 1 3年第 8期 料给出待反演参数的取值范围， 见表 2 。首先根据水压分量的 实测值及待反演参数的取值范围， 采用混沌遗传算法反演出 1 1 坝段的瞬时弹性模量； 再将瞬时弹性模量作为已知值， 根据 水压分量和时效分量的和， 反演得到 1 1 坝段的延迟弹性模量 和黏性系数， 结果见表 2 。 表 1 坝体和坝基物理力学参数 表 2 待反演参数取值范围及黏弹性参数反演结果 为了验证结果的正确性 ， 将反演得到的黏弹性参数作为已 知值 ， 计算 P L 1 1 2 2测点在相应时间段变形的水压和时效 分量之和， 并与分离出的实际测值进行对比， 见表 3 。由表3可 以看出， 计算值与实测值比较接近， 说明计算结果是合理的， 求 解方法是可行的。 表 3 水平位移计 算值与实测值对 照 mm 实测值 黏弹性模型计算值 最大值 最小值 最大值 最小值 l 1 6 O 8 4 4 l 1 8 4 8 6 7 4 结语 根据碾压混凝土坝的成层结构特性， 将层面影响概化到整 个坝体中， 将碾压混凝土坝视为横观各向同性体， 研究了采用 双向异性的三元件模型模拟碾压混凝土坝黏性效应的方法， 建 ( 上接第 1 0 1页) 为了避免指标之间的相互关联和重叠， 使设计的指标能够 客观 、 全面地反映移民可持续生计水平， 可以采取专家咨询法 等手段对指标进行归属度分析， 剔除归属度低的指标 ， 从而得 到最能反映移民可持续发展水平的指标体系。 3 结语 走出水利水电工程项目移民“ 发展悖论” ， 形成各方利益均 衡局面是促进移民社会稳定、 生活水平提高的关键。在新时 期， 移民安置需要坚持“ 有发展的安置” 这一思路， “ 有发展的 安置” 是社会发展的内在属性和要求 ， 坚持“ 有发展的安置” 有 利于促进移 民“ 落地生根” 和“ 安居乐业” ， 否则可能诱发各种 社会冲突。 立碾压混凝土坝的等效黏弹性模型。相 比于采用裂缝单元 、 薄 层单元等直接模拟层面， 该等效模型原理简单、 便于计算 ， 并从 宏观上有效地模拟了碾压混凝土坝的各向异性。根据位移实 测资料 ， 采用偏最小二乘法建立位移统计模型， 分离水压分量 和时效分量 ， 建立了反演碾压混凝土坝双向异性黏弹性参数的 目标函数 ， 利用混沌遗传算法对黏弹性参数进行 了反演， 并研 制了相应的反演分析程序。工程实例中， 将反演得到的黏弹性 参数作为已知值 ， 计算测点在相应时间段变形的水压和时效分 量之和， 并与分离出的实际测值进行对比， 结果表明： 计算值与 实测值比较接近， 说明反演方法是正确、 可行的。 参考文献 ： 1 贾金生， 改新， 锋玲， 等 碾压混凝土坝发展水平和工程实例 M 北京： 中 国水利水 电出版社 ， 2 0 0 6 2 彭友文， 郑东健， 吴中如 碾压混凝土坝层面弹性模量和厚度反演 J 水力 发电， 2 0 0 5 ， 3 1 ( 8 ) ： 2 7 2 9 3 G o o d m a n R E， T a y l o r R L ， B e r k k e T L A M o d e l fo r t h e M e c h a n i c s o f J o i n t e d R o c k J J A S C E ， 1 9 6 8 ， 9 4 ( 3 ) ： 6 3 7 6 5 9 4 Z i e n k ie w ic z 0 c ，B e s t B ，D u l l a g e C ， e t a 1 A n a ly s i s o f N o M i n e P r o b l e m s wit h P a r t i c u l a r Re fe I e n c e t o J o i n t e d R o c k S y s te m s J I n t e rna t i o n a l S o c i e t y f o r Rock Me c h ani c s ，1 9 7 0( 3) ：5 0 15 0 9 5 杨学堂， 哈秋龄， 张永兴， 等 裂隙岩体宏观力学参数数值仿真模拟研究 J 水力发电， 2 0 0 4 ， 3 0 ( 7 ) ： 1 4 1 6 6 顾冲时， 程乐群， 李婷婷 探讨碾压混凝土坝薄层单元有限元分析法 J 计 算力学学报， 2 0 0 4 ， 2 1 ( 6 ) ： 7 1 8 7 2 1 7 吴中如 水工建筑物安全监控理论及其应用 M 北京： 高等教育出版社， 2 0 0 3 8 易尚潭 基于混沌遗传算法的大地电磁测深反演 D 长沙： 中南大学， 2 0 0 4 9 姚俊峰， 梅炽， 彭小奇， 等 改进的混沌遗传算法及其在炼铜转炉操作优化 中的应用 J 中国有色金属学报， 2 0 0 1 ， 1 1 ( 5 ) ： 9 2 0 9 2 4 1 0 王惠文 偏最d x - - 乘 回归方法及其应用 M 北京： 国防工业出版社， 1 9 9 9 【 责任编辑昌艳梅】 参考文献 ： 1 余文学 水库移民问题社会经济分析 J 河海大学学报： 哲学社会科学版， 2 0 0 0 ( 4 ) ： 1 - 5 2 迈克尔 M塞尼 移民 重建 发展： 世界银行移民政策与经验研究( -) M 南京 ： 河海 大学出版社 ， 1 9 9 8 3 斯卡德 大坝的未来 M 北京： 科学出版社， 2 0 0 8 4 王茂福 移民安置目标基于福利经济学的分析 J 人口与经济， 2 0 0 3 ( 4 ) ： 7 3- 7 6 5 张宝欣 三峡开发性移民政策与实践 M 北京： 中国三峡出版社 ， 1 9 9 9 6 于谨凯， 高磊 基于 P S R模型的海洋生物资源可持续开发政府诱导研究 J 经济问题探索， 2 0 0 9 ( 9 ) ： 1 1 一 l 5 7 施国庆 水库移民社会经济影响分析 G 水库移民理论与实践 北京： 中 国水利水电出版社， 2 0 0 5 ： 1 2 7 1 3 0 【 责任编辑张华岩】 1 0 5