a
25
成果
• 编写出一条完善的计算崩塌体弹跳高度， 速度能量的源程序代码
a
26
博士阶段研究计划
一 用证据权法对国内地热进行评价和预测
二 地埋管提取地热能量理论研究及软件编写
三 中国科技部863计划，多重流体介质与岩体 应力耦合的理论研究与软件编写
a
27
一 用证据权法对国内地热进行评价和预测
研究内容：
3. 在证据权模型的基础上，对两类二进制证据 图处理，提出I/O和WofE两种方法进行评价。用 概率论两种判别方法对I/O和WofE两种方法进行 优劣性评价。最终得出合理的期望图。
a
28
技术路线
1 收集相关质料，从公共网站上收集中国地震震中数据 （震级、震源深度和经纬度坐标）然后用FORTRAN程序进 行编程处理网格划分并求解古登堡李斯特B值，处理成所需 要的证据图。同时其他三个原始数据图，包括中国断层图， 布格重力异常图，磁异常图。
在反求边坡稳定系数时，传统Sarma法 采用人工试算，这在编程计算边坡稳定性 时带来很大的局限性。在传统Sarma的基础 上，采用牛顿迭代法计算安全稳定系数， 并给出相应的推导公式。该法避免了人工 试算的繁琐，算法收敛迅速。
a
16
算 法 线 路 图
a
17
部 分 源 程 序 代 码
a
18
研究结果图示
a
11
部 分 源 程 序 代 码
a
12
研究结果图示1
0° ～ 150°临界地震系数变化
a
临界地震系数比较
13
研究结果图示2
临界地震系数和稳定性系数关系
a
临界地震系数与迭代
14
成果
• 编写出计算边坡稳定性的一套源程序代码 • 提出任意角度下临界地震系数求解公式 • 发表中文核心一篇
a
15
（3）Sarma法中使用牛顿迭代对安 全系数求解
安全系数初始预测值选取范围
不同地震系数值对应的安全系数
a
19
成果
• 提出一套简洁算法，快速有效求解安全系数 • 形成完善的计算边坡稳定的源程序代码 • 发表一篇EI(审核中)
a
20
（4）混凝土碳化区化学变化数值分 析和研究
基于混凝土碳化机理, 建立了水灰比、 水 泥用量等混凝土碳化主要影响因素与理论 模型中有效扩散系数及单位体积混凝土的 CO2 吸收量之间的定量关系, 利用推导得出 的混凝土碳化深度实用数学模型，进行编 程形成一套源程序代码。
临界地震系数值与地震力与垂直向下的夹
角呈抛物线分布， 存在极值。
地震力方向对临界地震系数有影响，但不影响安全 稳定系数
a
10
编程所用主要公式
Kc (Pn anPn1aenn1enPna2enne2enn1e n1.. ...P . 1enae1enn1e.n..e12..).es2in （PnPn1enPE n12eenneenn11.. .e1.. .EPn1e1nen1...e2） sin
TOUGHREACT 输入
FLAC3D 边界和计算控制
技 术 路 线 图
TOUGHREACT 计算输出
MESH文件
t
3D插值 T*,P*
FLAC3D
SAVE文件
, PC , K ti1 ti t
t tult
ENa D
4
技 术 难 点
： 中 心 点 与4 节 点 插 值 转
1
换
节点或角点
15
a
8
成果
• 在编写的耦合软件基础上，已发3篇EI和一 篇SCI
• 提出三维拉格朗日插值公式 • 准备申请软件专利
a
9
（2）地震力方向对Sarma法影响和研究
在原 Sarma 的基础上， 假设地震力非水平方向， 在任意角度上尝试给出推导公式，通过编程计算实 例，得出以下几点结论。
临界地震系数值关于地震力与垂直向下的夹角90° 左右对称。
注: (1)为国家自然科学基金( 40972172 )， 上海市科学技术委员会资助项目( 10dz1202300 )， 吉 林省留学人员科技创新创业项目( 3D5112204424) 和教育部留学回国人员科研启动基金 ( 3C111F284424) 资助．
a
2
（1） CO2灌注后盖层力学变化和 安全性研究
1. 选择资源评价相适应的软件，理论基础和依 据，找出信息因素所对应的参数范围。
2．选择合适的模型，本研究在成矿预测广泛采 用的证据权模型方法的基础上，通过一定的数学 计算方法确定考虑与地热相关的影响因素，进行 证据权评估。该模型以贝叶斯条件概率为基础, 通 过一定的数学计算方法确定与地热因素之间的关 系。
2四张原始数据图转化成研究所需要的四张证据图。布格 重力异常图与地热空间关系不明显，要处理成梯度图，而磁 异常图不需要任何转换直接做为证据图
3四张证据图用证据权法和平均值标准差法两种方法进行 二进制处理。证据权模型以贝叶斯条件概率为基础, 通过一 定的数学计算方法确定与地热因素之间的关系，而平均值标 准差法是去平均值，在标准差范围内的值都设为1其余为0。
主要研究CO2灌注后盖层力学变化和安全性情况， 本人主要负责编写反映溶质运移某一软件与FLAC3D耦合搭接的某
一程序代码，模拟一定时间后CO2灌注后应力温度等的进展与情 况。
主要解决的技术难题有：
FLAC3D与TOUGH2及TECPLOT三个软件的网格转换
节点和中心点转换通过三维三点拉格朗日插值
a
3
难点：涉及时间和空间，公式中含有导数， 采用一维有限差分编程。
a
21
部 分 源 程 序 代 码
a
22
成果
• 形成一套计算混凝土碳化区化学变化的源 程序代码
a
23
（5）崩塌拦石网高度的计算和研究
• 研究内容：主要编写相关参数计算的源程 序代码
• 崩塌体落体能量，速度，及崩塌体的弹跳 高度
a
24
部 分 源 程 代 码
16
18
中ห้องสมุดไป่ตู้点
6
8
12
3•
7 13
1•
3
4•
14 11
2•
5
17 10
2
9
a
5
文件关系调用图
Mesh文件数据 Flac3d网格建立
生成文件1 Flac3d调用文件1
生成Flac3d 文件3
Tough react 文件
Save文件数据 三维三点插值
生成文件2
a
6
为解决搭接而提出的插值公式
a
7
部 分 源 程 序 代 码
博士阶段研究计划
学校：吉林大学 学院：建设工程学院 姓名: 李建明 报考导师：张延军
a
1
一 硕士阶段工作总结
硕士阶段主要完成工作： （1） CO2灌注后盖层力学变化和安全性研究 （2）地震力方向对Sarma法影响和研究 （3）Sarma法中使用牛顿迭代对安全系数求解 （4）混凝土碳化区化学变化数值分析和研究 （5）崩塌拦石网高度的计算和研究