地学数值模拟技术01

1、2009年资源与地球科学学院地学数值模拟技术课程提纲 参考文献资源与地球科学学院地学数值模拟技术数值模拟的核心技术 潘钟祥潘钟祥. 石油地质学石油地质学. 北京：地质出版社，北京：地质出版社，1986 陈荣书陈荣书. 天然气地质学天然气地质学. 北京：中国地质大学出版社，北京：中国地质大学出版社，1986 李思田，林畅松李思田，林畅松. 沉积盆地分析沉积盆地分析. in：肖庆辉主编：肖庆辉主编. 当代地质科学前沿当代地质科学前沿. 北京：中国地质大学北京：中国地质大学出版社，出版社，1993 郝石生，陈明章，高耀斌郝石生，陈明章，高耀斌. 天然气藏的形成和保存天然气藏的形成和保存. 北京：石油

2、工业出版社，北京：石油工业出版社，1995 庞雄奇，陈章明，陈发景庞雄奇，陈章明，陈发景. 含油气盆地地史、热史、生留排烃史数值模拟研究与烃源岩定量含油气盆地地史、热史、生留排烃史数值模拟研究与烃源岩定量评价评价. 北京：地质出版社，北京：地质出版社，1993 李明诚，李明诚， 石油与天然气运移石油与天然气运移(第二版第二版). 北京：石油工业出版社，北京：石油工业出版社，1994 j.w.哈博，哈博，c.博纳姆博纳姆-卡特卡特. 地质过程的计算机模拟地质过程的计算机模拟. 北京：地质出版社，北京：地质出版社，1986 范士芝，姜鹏，陶一川等范士芝，姜鹏，陶一川等. 塔里木盆地北部油气运移二维

3、二相流数值模拟分析塔里木盆地北部油气运移二维二相流数值模拟分析. 地球科学地球科学中国地质大学学报，中国地质大学学报，1995，20(3)：321327 韦重韬，煤层甲烷地质演化史数值模拟，徐州：中国矿业大学出版社，韦重韬，煤层甲烷地质演化史数值模拟，徐州：中国矿业大学出版社，1998 陆金甫，顾丽珍，陈景良陆金甫，顾丽珍，陈景良. 偏微分方程差分方法偏微分方程差分方法. 北京：高等教育出版社，北京：高等教育出版社，1988 plays h w, flaneily b p, tokolsky s a, et al. 数值方法大全数值方法大全. 王璞，何玉江，苗天住译王璞，何玉江，苗天住译. 兰

4、州：兰州：兰州大学出版社，兰州大学出版社，1991 作业要求资源与地球科学学院地学数值模拟技术数值模拟的核心技术针对某一讲座内容，查阅参考文献，编写图书报告。针对某一讲座内容，查阅参考文献，编写图书报告。要求：要求： 反映本学科领域的发展现状和发展趋势；反映本学科领域的发展现状和发展趋势； 字数字数3000-5000。编写计算机程序编写计算机程序天然气扩散方程的结算天然气扩散方程的结算条件：条件： 烃源层不生烃，初始含气量烃源层不生烃，初始含气量30m3/m3; 扩散通道和储层的初始甲烷浓度均为扩散通道和储层的初始甲烷浓度均为0； 扩散通道厚度为扩散通道厚度为30m，扩散系数为，扩散系数为2.

5、5*10-8； 所有岩石的密度均为所有岩石的密度均为2.5g/cm3； 模拟模拟100ma之后，储层的甲烷含量。之后，储层的甲烷含量。 模拟技术-什么是模拟？/simulation/modeling 模拟建立模拟现实的模型，并以此模型做实验。其目的是某一过程的理解、预测和控制。 在地质学中，模拟的主要作用是理解地质过程和控制、预测地质过程的结果。资源与地球科学学院地学数值模拟技术概述水平羽状煤层气井 平均单井71120m3/d，产水75m3 模拟技术-模拟技术的种类 数值模拟建立地质/物理模型、数学模型、编制计算机程序，利用程序的设定不同条件运行程序，考察地质过程的性状，预测地质过程的结果。资

6、源与地球科学学院地学数值模拟技术概述水平羽状煤层气井 平均单井71120m3/d，产水75m3 模拟技术-模拟技术的种类资源与地球科学学院地学数值模拟技术概述 物理模拟/实验室模拟构建地质/物理模型，根据模型，采用适当的材料建立缩小比例的实体模型以及所需的外界物理、化学条件，由此模仿实际的地质过程，并考察地质过程的性状，预测地质过程的结果。煤体三轴应力-应变模拟测试装置模拟地层条件下吸附甲烷的煤体在应力作用下的变形特征昆士兰大学物理、化学和建筑系 模拟技术-模拟工作的步骤资源与地球科学学院地学数值模拟技术概述 定义系统系统是由若干成分组成的，当系统的任何一部分发生改变时，将会影响到整个系统。建

7、立系统时需要考虑的内容： 系统的元素：即系统的组成部分 系统的界限：即确定系统的外边界，因为考虑的问题肯定是一个大系统的一小部分，所以必须做这个工作 系统的内在联系：通常要建立组成元素内在联系结构图 物质和能量的迁移和记载：运送和储存的过程、物质的交换等等 时间流程：动态系统的变化，很大程度上是以时间为基础来进行的，在描述一个系统时，没有时间因素几乎是不可能的。所以，时间是系统中的一个最重要的因素。此外，在定义系统时需要非常注意的一个问题是系统的简化。因为地质过程是非常复杂的，如果不简化，对下一步的建模、计算非常不利，有时甚至是不可解。 模拟技术-模拟工作的步骤资源与地球科学学院地学数值模拟技

8、术概述 建立模型地质模型和数学模型模型确性检校模型地质模型/物理模型数学模型一种假说，定性地描述整个地质过程对地质模型的定量表述采用对已知的过程进行模拟，观察模拟结果是否与实际结果相符合来进行opening faultground surfacecoal organic mattergenerates gasdiffusion dissipation inverticle directioncap outburst dissipationpermeation dissipationalong cleats and faults general description of the proces

9、s by mass conservationgsd=qqqqg quantity of gas generated by coal organic material, m3/t;qs gas storing in coal reservoir, m3/t;qd gas dissipated from the coal reservoir, m3/t; gas generation qg by thermal simulation test of coal organic matterggo,maxn()qf rmro, max vitrinite reflectance of coal, %;

10、mn pure organic material in coal (dry, ash-free) , t/t-coal; and fg(ro, max)432go,maxo,maxo,maxo,maxo,max()1.691() -23.216()90.950() -23.928()-2.831f rrrrrin qinshui basin:maturation history of coal organic material is a key toaccess gas generation history 模拟技术-模拟工作的步骤资源与地球科学学院地学数值模拟技术概述观察模型的性态研制计算机

11、程序观察模型的性态编程、调试、检校、调整模型、完成程序研制在模拟软件的入口输入不同的条件，观察和研究模拟的结果，包括过程的时间变化特征、最终结果特征等等，从而达到认识自然过程、根据模拟结果进行控制、预测的目的 数值模拟技术在地质学中的应用资源与地球科学学院地学数值模拟技术概述 油气地质学 盆地模拟/储层模拟 构造地质学 沉积学 水文地质学 工程地质学 学科领域 模拟内容 相关软件petromod (ies)/temispack (ifp)/basinmod (platte river inc.)coal gas (holdich inc.)/comet (ari) ies（德国有机地化研究所）

12、ifp（法国石油研究院）wasy（德国水资源规划与系统研究所）adina(automatic dynamic incremental nonlinear analysis)feflow (finite element subsurface flow system) 构造应力场模拟 沉积体形成过程模拟 地下水动力学过程 岩体应力-应变feflow (wsy) adina (adina r&d co.) adina (adina r&d co.) matlab (mathworks co.) matlab (mathworks co.) 数值模拟技术的核心偏微分方程资源与地球科学学

13、院地学数值模拟技术数值模拟的核心技术 偏微分方程（partial differential equations/pde方程） 含有未知函数及其偏导数的方程。拉普拉斯方程： 偏微分方程的依据 质量守恒定律：参加反应前各物质的质量总和等于反应后生成各物质的质量总和 能量守恒定律：能量既不会凭空产生，也不会凭空消失，它只能从一种形式转化为 别的形式，或者从一个物体转移到别的物体，在转化或转移的过程 中其总量不变0222222zuyuxu u为一个未知的函数 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 根据质量守恒定律推导的天然气扩散运移方程 地质模型：扩散是一种

14、自然界中常见的物质传递过程，特点： 动力是性状的差异，如浓度差、温度差、密度差； 速度非常缓慢，天然气的扩散运移是一种浓度扩散； 遵循费克定律，即扩散的量与扩散系数成正比，与浓度梯度成正比，假定有天然气（甲烷）从烃源层已扩散的方式向储层运移烃 源 层扩散通道储 层在天然气运移的通道上选择的任意单元控制体 zxyz扩散入扩散出 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 根据质量守恒定律推导的天然气扩散运移方程根据质量守恒定律，在某一时刻t，有控制体内部的甲烷质量变化率扩散进入控制体的甲烷量扩散出控制体的甲烷量甲烷在控制体内质量的变化率为：tzyxtc zx

15、yz扩散入扩散出c为t时间内控制体中的甲烷浓度由费克定律，在t时间内沿z轴扩散进入控制体的甲烷质量为：tzyxnec. ：煤层甲烷的通量矢量ecn.tzyxzntzyxnecec.同理，扩散出控制体的甲烷质量为：根据质量守恒和控制体 的任意性：zyx0.zntcec 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 根据质量守恒定律推导的天然气扩散运移方程表示甲烷浓度变化的pde方程 zxyz扩散入扩散出cdnec.d：煤层甲烷在封盖层中的扩散系数 ：哈密尔顿算子222222zcycxc所以：0d22tczc烃 源 层扩散通道储 层进一步：cz表明在t时刻，从烃

16、源层到扩散层，甲烷浓度c在空间上的变化量；表明，甲烷浓度是空间和时间的函数；扩散运移量可以通过费克定律来计算),(tzfc 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 pde方程的解算数值法/numerical solution目前常见的解算方法有： 经验公式法 有限差分法 有限元法 边界元法 0d22tczc本质上是想办法把方程转化为加、减、乘、除的计算按照偏导数的原始定义，把 、 和 等离散化：22zczctc所以，在t=j 时刻：将扩散通道等分为n份，则每份的高度为n+1ni+1ii -1210z同时假定，总的扩散时间为t，将其等分为m份，则扩散时间

17、可以表示为：zcmt nhh/h 和h分别称为时间步长和空间步长 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 pde方程的解算数值法/numerical solution-离散化的过程i-1ii+1h有：z在t=j 时刻：thcchcchcczcjijijijijijiji21111hhcchcczcjijijijiji2112221122hcccjijijii-1ii+1hz在t=（j+1） 时刻：tjicjic1jic11jic11jic11jicjijijicctc1以上为偏微分的差分格式 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技

18、术数值模拟的核心技术 pde方程的解算数值法/numerical solution-离散化的过程各种差分格式hcczcjijiji21111显式蛙跳格式hcczcjijiji211隐式蛙跳格式hcccczcjijijijiji4111111隐式二层格式2112222hccczcjijijiji隐式格式2111111122422hcccccczcjijijijijijiji隐式二层格式1jijijicctc隐式格式i-1ii+1hzjicjic1jic10d22tczc 数值模拟技术的核心偏微分方程资源与地球科学学院地学数值模拟技术数值模拟的核心技术 pde方程的解算数值法/numerical solution-离散化的过程令：022d1211jijijijijicchccc通式：22dh有：111)2(jijijijijiccccc由此得线性方程组jmjmjmjmjmjmjmjmjmjmjijijijijijjjjjjjijjjccccccccccccccccccccccccc11111111111112123201121)2()2()2()2()2(转化为：jmjmjmjjjjmjmjjrccccrccccccrrrrrrrrrr11111201112121212121三对角线性方程组，很容易求解但在求解之前，需要设定初始条件和边界条件pde方程差分后，可能