数值模拟方法进展
油藏数值模拟
名词解释 油藏模拟油藏数值模拟数学模拟物理模型数值模型质量守恒定律适定问题初始条件黑油模型组分模型网格节点块中心网格点中心网格离散化有限差分法显示差分 隐式差分前差分后差分中心差分点交替排列格式交替对角排列格式标准排列格式 对角排列格式隐式差分格式差分方程稳定性截断误差松弛法IMPES方法历史拟合 动态预测灵敏度实验 选择题 由于油藏各点的渗透率不同，束缚水饱和度不同，因而需要对相对渗透率曲线进行归一化处理 以X方向为例，传导系数为 块中心网格是用（）来表示小块坐标的 A网格块中心B节点C网格块边缘D网格块夹角 下述表达式表示定产量内边界条件的是 认识油田的主要方法有直接观察法和模拟法 相对渗透率取值一般取上游权的处理方法 IMPES方法是（）的求解方法 A隐式压力B隐式饱和度C全隐式 历史拟合在含水拟合时主要是对（）的修改 A孔隙度B相对渗透率曲线C渗透率D地层厚度 在隐式差分格式中，有多个未知数，当已知第n时刻的值P i n时，为了求出第n+1时刻的P i n+1，需要（） A解n个方程B解一个线性代数方程组C直接求解D解一个方程 根据每一组分的质量守恒建立的渗流数学模型称为（）模型 A热采B化学驱C黑油D组分 一维径向模拟时r=10cm，r=40cm，那么可以推断r s的大小是 A120 B200 C400D 640 下列哪一种方法不属于迭代求解方法 A雅克比法B超松弛法CLU分解法D交替方向隐式法 对于二位6*4网络系统，如果按行标准排列，气半带宽W= A6 B4 C12 D8 克兰克?尼克森差分格式的截断误差为（） 块中心网格和点中心网格的差分方程相比较，结果（） A一样的B有半个网格的误差C相差流动项系数D维数不同 三．判断题2分*10 1.黑油模型中水相与其他两相不发生质量转移，气可以从油中出入，但不能汽化液相 2.离散化的核心是把整体分为若干单元来处理，它是油藏对象的空间离散 3.显式差分格式是有条件收敛的 4.差分方程组的直接解法的特点是计算工作量小，精确度较高，计算程序简单 5.差分方程组的迭代解法主要用于处理系数矩阵阶数较高的问题 6.相对渗透率取值一般取上游权的处理方法 7.油藏模拟的基础在于油藏描述和生产动态，若油层参数和生产数据不准确，通过数值模 拟的算法也可以消除 8.显示差分格式的稳定条件是△t/△x2≤0.5 9.有限差分法就是用差商来代替微商
西南石油大学石油工程专业开题报告模板
本科毕业设计（论文）附件 题目：多层油藏不稳定渗流理论研究 学生姓名学号 教学院系石油工程学院 专业年级石油工程2###级 指导教师职称 单位西南石油大学
一、设计（论文）的选题意义及国内外研究现状 1、选题的意义 随着中国经济的高速发展，对石油的需求量也不断增加。2010年上半年中国原油需求量达875万桶/天，同比增长19.8%。从全球来看，中国石油需求增速领先于其他新兴市场以及世界其他地区，全球石油消费有一半来自中国。 正确评价已开发油藏的动态参数，尤其是剩余油分布是科学、合理制定提高采收率措施方案的基础，利用试井测试资料确定储层参数和剩余油分布是一种简便、经济、可靠的实用方法。我国东部地区是我国石油生产的主要集中地，但东部大多为老油田（如大庆、胜利等），都已进入高含水时期，石油产量逐年下降，实现稳产已是十分困难。为了保证老油田稳产、挖掘老油田的潜力，必须实施更精细化的开发。这其中就包括对油藏进行更精细化的描述，原来作为单层处理的油层需要根据它的非均质性划分为多层，并获得各层的渗透率等地层参数，为油田的后续开发提供依据。另外，水驱采油是目前国内的主要采油方式，为提高油井产量，需要改变原来的传统注水方式，根据剩余油分布给各油层分配不同的注水量。 综上可知，多层油藏不稳定渗流理论的研究对提高老油田的采收率、稳定老油田的产量、实施更精细化的开发都有重要的意义。 2、国内外研究现状 20世纪50年代以来，试井已成为人们认识油藏的重要手段。随着现代科技的发展，特别是进入20世纪80年代以来，计算机技术的发展以及试井分析新方法的提出极大的丰富了试井解释功能，有效的提高了分析结果的准确性。长期以来，建立在对渗流方程进行解析求解基础上的常规试井方法一直在试井领域占据着主导地位，特别是在20世纪80年代Bourdet图版问世以来，结合Bourdet典型曲线与Gringarten典型曲线对诸如渗透率、原始地层压力、地层表皮系数及
16种常用数据分析方法
一、描述统计描述性统计是指运用制表和分类，图形以及计筠概括性数据来描述数据的集中趋势、离散趋势、偏度、峰度。 1、缺失值填充：常用方法：剔除法、均值法、最小邻居法、比率回归法、决策 树法。 2、正态性检验：很多统计方法都要求数值服从或近似服从正态分布，所以之前需要进行正态性检验。常用方法：非参数检验的K-量检验、P-P图、Q-Q图、W 检验、动差法。 二、假设检验 1、参数检验 参数检验是在已知总体分布的条件下（一股要求总体服从正态分布）对一些主要的参数（如均值、百分数、方差、相关系数等）进行的检验。 1）U验使用条件：当样本含量n较大时，样本值符合正态分布 2）T检验使用条件：当样本含量n较小时，样本值符合正态分布 A 单样本t检验：推断该样本来自的总体均数卩与已知的某一总体均数卩0 （常为理论值或标准值）有无差别； B 配对样本t 检验：当总体均数未知时，且两个样本可以配对，同对中的两者在可能会影响处理效果的各种条件方面扱为相似； C 两独立样本t 检验：无法找到在各方面极为相似的两样本作配对比较时使用。 2、非参数检验 非参数检验则不考虑总体分布是否已知，常常也不是针对总体参数，而是针对总体的某些一股性假设（如总体分布的位罝是否相同，总体分布是否正态）进行检验。 适用情况：顺序类型的数据资料，这类数据的分布形态一般是未知的。 A 虽然是连续数据，但总体分布形态未知或者非正态； B 体分布虽然正态，数据也是连续类型，但样本容量极小，如10 以下； 主要方法包括：卡方检验、秩和检验、二项检验、游程检验、K-量检验等。 三、信度分析检査测量的可信度，例如调查问卷的真实性。 分类： 1、外在信度：不同时间测量时量表的一致性程度，常用方法重测信度 2、内在信度；每个量表是否测量到单一的概念，同时组成两表的内在体项一致性如何，常用方法分半信度。 四、列联表分析用于分析离散变量或定型变量之间是否存在相关。对于二维表，可进行卡 方检验，对于三维表，可作Mentel-Hanszel 分层分析列联表分析还包括配对计数资料的卡方检验、行列均为顺序变量的相关检验。 五、相关分析 研究现象之间是否存在某种依存关系，对具体有依存关系的现象探讨相关方向及相关程度。 1、单相关：两个因素之间的相关关系叫单相关，即研究时只涉及一个自变量和一个因变量； 2、复相关：三个或三个以上因素的相关关系叫复相关，即研究时涉及两个或两个以
CO2驱油数值模拟研究现状与发展趋势
注气驱油数值模拟方法研究现状与发展趋 势 姬泽敏，秦积舜，李实，廉黎明 （提高石油采收率国家重点实验室—中国石油勘探开发研究院，北京市 100083） 摘要：注气驱油技术是提高石油采收率的重要方法之一，应用和发展前景广阔。气驱油过程伴随着油气体系间的组分传质和系统压力变化，进而引起油气体系的相态转化，使得注气驱油过程的物理化学现象的表征和数学描述变得十分复杂，至今尚未形成统一和精确的油气体系相态表征和描述方法。通过考察国内外已有的相关数学模型和计算模拟方法，本文较为系统的梳理了注气驱油数值模拟方法的发展历程，评价了现有方法的优缺点，并结合我国油藏储层及流体特征，提出了适合中国油藏特点的注气驱油数值模拟方法的发展方向。 关键词：注气驱油技术；油气体系；数值模拟；组分传质；相态 Research Status of Gas Flooding Numerical Simulation and Its Development Trend Ji Zemin1，Qin Jishun，Li Shi，Lian Liming (The State Key Laboratory of Enhanced Oil Recovery—RIPED, Beijing 100083,China) Abstract: Gas flooding technology is one of the pivotal EOR methods, which has a broad prospect of application and development. However, mass transfer and change of system pressure during the process of gas flooding lead to the change of phase behavior, which draws great difficulties to the mathematical description and characterization of physical and chemical phenomenon during the process, so far there has not been a set of uniform and accurate methods to describe and characterize the phase behavior of oil and gas system. According to the mentioned above, based on the investigation of several gas flooding numerical simulation methods at home and abroad, this paper hackled the development process of the gas flooding numerical simulation methods, evaluated the advantages and disadvantages of these methods. Finally, combined with the characters of reservoirs and fluid, development direction of gas flooding numerical simulation with Chinese characteristics was proposed. Key words：gas flooding，oil-gas system，numerical simulation，compositional transfer，phase 1收稿日期： 第一作者简介：姬泽敏（1985—），男，博士研究生，主要从事注气提高采收率技术及数值模拟研究。 基金项目:国家973项目04课题“孔隙介质中相态实验与理论研究”(No. 2011CB707304)，国家科技重大专项(No. 2011ZX05016-001)。
第三章-数值模拟理论与方法
第三章 数值模拟理论与方法 §3.1 流体力学的基本方程 流体运动所遵循的规律是由物理学三大守恒定律规定的，即质量守恒定律，动量守恒定律和能量守恒定律[44]。 （一）连续方程 0)(=?+??v t ρρ （3.1） 式中 ρ－流体密度 u －流体速度分量 （二）动量方程（x 方向） 对于不可压流体（即0=?v ） x p f v u v x u x ??-+??=??+??ργρρρ)()()( （3.2） 式中 γ－运动粘性系数 p －压力 对于可压缩流体 ()()()()()x p f v x u u v x u x ??-+???+??=????ργργρρρ 31 （3.3） 式中等号后前两项是粘性力 y ，z 方向上的动量方程可类似推出。 （三）能量方程 ()()()v q T k e v e t ερρ++???=??+?? （3.4） 其中 T C e v = 式中等号左边第一项是瞬变项，第二项是对流项，等号右边第一项是扩散项，第二、三项是源项。 所以，流体力学基本方程组为： ()0=?+??v t ρρ
()x p f u u v f t u x ??-+??=??+??ργρ)( ()()y p f v v v f t v y ??-+??=??+??ργρ （3.5） ()()w p f w w v f t w w ??-+??=??+??ρλρ ()()v q e c k e v f e t v ερ++??? ? ????=??+?? §3.2 紊流模式理论概况 §3.2.1 基本方程 在自然界中，真实的流体都具有粘性。粘性流体存在两种不同的运动方式和流态，即层流和紊流。而在自然界和工农业生产中所遇见的流体流动大部分都是紊流。 三维的N-S 方程是目前描述粘性流体运动较为理想的模型，其优点一是应用范围广，在空气、水流、传热等方面均用N-S 方程描述；二是对于有分离、旋涡等情况的复杂三维流动更为适用。 三维直角坐标下的N-S 方程[45]，[46]，即不可压缩粘性流体的动量方程式为： ?????????????+??+??+??-=??+??+??+??-=??+??+??+??-=)()()(222222222222222222z w y w x w z p F Dt Dw z v y v x v y p F Dt Dv z u y u x u x p F Dt Du z y x μρρμρρ μρρ （3.6） 不可压缩流体的连续性方程为： （3.7） 式（3.6）和（3.7）共有四个未知数（u 、v 、w 、p ）和四个方程，加上边界条件，从理论上来讲其解是存在的。但是，要直接求解复杂而详细的粘性流体运动是十分复杂和困难的。其原因是：直接求解N-S 方程要求求解从反映消散运动的最小涡漩尺度到反映大尺度涡体的所有流动尺度，因而只有对简单情况下才有理论解。 0=??+??+??z w y v x u
油藏数值模拟方法
第一章油藏数值模拟方法分析 油藏数值模拟 油藏数值模拟简述 油藏数值模拟是根据油气藏地质及开发实际情况，通过建立描述油气藏中流体渗流规律 的数学模型，并利用计算机求得数值解来研究其运动变化规律。其实质就是利用数学、地质、物理、计算机等理论方法技术对实际油藏的复制。其基础理论是基于达西渗流定律。 油藏数值模拟就是利用建立起的数学模型来展现真实油藏动态，同时采用流体力学来模 拟实际的油田开采的一个过程。基本原理是把生产或注人动态作为确定值，通过调整模型的不确定因素使计算的确定值（生产动态）与实际吻合。其数学模型，是通过一组方程组，在一定假设条件下，描述油藏真实的物理过程。充分考虑了油藏构造形态、断层位置、油砂体分布、油藏孔隙度、渗透率、饱和度和流体PVT性质的变化等因素。这组流动方程组由运动方程、状态方程和连续方程所组成。油藏数值模拟是以应用数学模型为基础的用来再现油田实际生产动态的过程。具体是综合运用地震，地质、油藏工程、测井等方法，通过渗流力学，借助大型计算机为介质条件建立三维底层模型参数场中，对数学方程求解重现油田生产历史，解决实际问题。 油藏数值模拟技术从50年代的提出到90年代间历经40年的发展，日益成熟。现在进入另外一个发展周期。近十年油藏数值模拟为油田开发研究和解决实际决策问题提供强有力的支持。在油田开发好坏的衡量、投资预测及油田开发方案的优选、评价采收指标等应用 非常广泛。 油藏数值模拟功能包括两大部分：①复杂渗流力学研究，②实际油气藏开发过程整体模 拟研究，且可重复、周期短、费用低。
图1油藏数值模拟流程图 油藏数值模拟的类型 油藏数值模拟类型的划分方法有多种， 划分时最常用的标准是油藏类型、 需要模拟的油 藏流体类型和目标油藏中发生的开采过程， 也可以根据油气藏特性及开发时需要处理的各种 各样的复杂问题而设定， 油气藏特性和油气性质不同， 选择的模型也不同， 还可以根据油藏 数值模拟模型所使用的坐标系、空间维数和相态数来划分。 以油藏和流体类型来划分，其模型有：气体模型、黑油模型和组分模型； 以开采过程来 划分，其模型包括：常规油藏、化学驱、热采和混合驱模型。 以油藏和流体描述为基础的油藏模型分为两类：黑油模型和组分模型。 （1） 黑油模型，是常规油田开发应用的油藏数值模型，用于开采过程中，对油藏 流体组分变化不敏感的情况， 是最完善、最成熟的。黑油模型假设质量转移完全取决于 压力变化，适应于油质比较重的油藏类型，在这些模型中，流体性质 E O 、B g 、R S 决定PVT 的 变化，如普通稠油及中质油的油气藏。 （2） 组分模型，应用于开采过程中对组分变化敏感的情况。这些情况包括：挥发性油 藏和凝析气藏的一次衰竭采油阶段， 用组分模型进行模拟。在组分模型中，适用于油质比较轻、气体组分比较高的油气藏， 使用 数据化 流体的PVT 数据、相 渗曲线、岩石数据 建立地质模型 建立网格 参数场 表格数据 油水井产量、井史 数据 T 动态模拟 含油边界拟合 非井点地质静态参数拟合 区块、单井压力拟合 生产指数拟合 以及压力保持阶段。同时，多次接触混相过程通常也采
低渗透油藏数值模拟技术研究
低渗透油藏数值模拟技术研究 摘要 随着现代石油工业的发展，低渗透油田的开发愈来愈为人们所重视。与中、高渗油藏相比较，低渗透油藏具有以下明显渗流特征：流体渗流不遵循达西定律，渗流中存在启动压力梯度，并且应力敏感性影响不可忽略。 本文在广泛调研的基础上，综合考虑启动压力梯度和应力敏感性，对低渗透油田的数值模拟方法进行研究。主要考虑油水两相具有相同的启动压力梯度，并且启动压力梯度的大小与地层渗透率和含水饱和度有关，同时考虑地层弹塑性变形引起的渗透率改变。在此基础上，建立了符合实际情况的二维油水两相数学模型，并通过IMPES方法进行数值求解，对低渗透油藏注水开发的生产特征，产能影响因素及不同注采井网方式进行了研究，研究表明启动压力梯度、应力敏感性、油水粘度差和毛管力对低渗透油藏产能影响明显。 关键词：低渗透油藏，启动压力梯度，应力敏感性，IMPES方法，数值模拟
Abstract With the development of modern oil industry, the exploitation of low permeability reservoir has been paid more and more attention in recent years. Compared with medium and high permeability reservoirs, the low permeability reservoir has different seepage flow characteristics,which are as follows: The flow of the fluid does not obey Darcy`s Law; The starting pressure gradient exists in transfusion; The stress sensitivity influence can not neglect. This paper, based on widespread investigation and study, considered starting pressure gradient and stress sensitivity influence, makes some research on numerical simulation of low permeability reservoir. Specifically, it is considered that oil and water have the same starting pressure gradient, which has the relationship with absolute permeability and water saturation.At the same time, it is also considered that the absolute permeability will vary, because of stratum elastoplasticity distortion. On this basis,a two-dimensional mathematical model of oil-water two-phase is set up, which is solved with the IMPES numerical method. The production characteristics, the factors affecting the deliverability and the different injection pattern were studied. The result indicates that starting pressure gradient、stress sensitivity、oil/water viscosity ratio and capillary pressure have obvious effect to the deliverability of low permeability reservoir. Key words: Low permeability reservoir, Starting pressure gradient, Stress sensitivity,The method of IMPES, Numerical simulation
湍流的数值模拟方法进展
3 大涡模拟（LES ） 湍流大涡数值模拟（LES ）是有别于直接数值模拟和雷诺平均模式的一种数值模拟手段。利用次网格尺度模型模拟小尺度湍流运动对大尺度湍流运动的影响即直接数值模拟大尺度湍流运动， 将N-S 方程在一个小空间域内进行平均（或称之为滤波），以使从流场中去掉小尺度涡，导出大涡所满足的方程。 3.1 基本思想 很多尺度不同的旋涡一起组成了湍流运动平均流动主要取决于大漩涡的流动，大尺度运动则受到小旋涡的影响。流动中的大涡实现了动量、能量质量、热量的交换，耗散主要是由于小涡作用的。大旋涡中受到流场形状、阻碍物的影响，，使大漩涡的各向异性更加明显。然而小漩涡之间各项同性，相互没有太大的区别，所以建立统一的模型比较容易一些。综上所述，大涡模拟将湍流瞬时运动量通过滤波将运动分成小尺度和大尺度。大尺度的运动受到小尺度的运动的影响可以通过应力项（类似于雷诺应力项）来表示，即为亚格子雷诺应力，以建立这种模型的方法来模拟。而大尺度则是求解运动微分方程而计算出来的，也就是说大涡模拟，要先过滤掉小尺度的脉动，然后再推出小尺度的运动封闭方程以及大尺度的运动控制方程。 3.2 滤波函数 正如上面提到，大涡模拟要先将流动变量分解成小尺度量和大尺度量，我们把这个作用叫做滤波。滤波运算就是在一区域内按照一定的条件对函数进行加权平均，作用是将高波数滤掉，使低波数保留，滤波函数的特征尺度决定了截断波数的最大波长，下面三种滤波函数是最为常用的主要有以下三种：盒式、富氏截断以及高斯滤波函数。 不可压常粘性系数的湍流运动控制方程为N-S 方程: j ij i j j i i x S x P x u u t u ???+??-=??+??)2(1γρ 式中：S 拉伸率张量，表达式为：2/)//(i j j i ij x u x u S ??+??=；γ分子粘性系数；ρ流体密度。设将变量i u 分解为方程(11)中i u 和次网格变量(模化变量)'i u ，
中石油2012年秋《油气田开发方案设计》在线作业答案
中石油2012年秋《油气田开发方案设计》在线作业答案第1—3阶段 第一阶段在线作业答案 第1题生产试验区是在详探程度高的地区，任意划出一块面积，用正规井网正式投入开发，并进行各项开发试验。 您的答案：错误 题目分数：0.5 此题得分：0.5 批注：考查生产试验区的概念 第2题试油是在油井完成后（固井、射孔），把某一层的油气水从地层中诱到地面上来，并经过专门测试取得各种资料的工作，一般时间较短。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查试油的概念 第3题区域勘探(预探)阶段主要任务是寻找和查明油气田，计算探明储量，为油气田开发做好准备。 您的答案：错误 题目分数：0.5 此题得分：0.5 批注：考查勘探各阶段的划分 第4题油藏模拟就是模仿和再现油藏内流体的渗流现象和渗流规律。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查油藏模拟的概念 第5题油气田开发方案是集合勘探开发各工程环节的数据信息，对油气田开发给出最佳的综合技术对策。 您的答案：正确
题目分数：0.5 此题得分：0.5 批注：考查对油气田开发方案的理解 第6题黑油模型主要用于模拟因粘滞力、重力和毛管力作用而引起的油、气、水三相的等温流动的模型。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查黑油模型的概念 第7题随机建模指以已知的信息为基础，以随机函数为理论，应用随机模拟方法，产生可选的、等概率的储层模型的方法，即对井间未知区应用随机模拟方法给出多种可能的预测结果。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查随即建模的概念 第8题油田开发方案编制的原则是确保油田开发取得好的经济效益和较高的采收率。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查对油田开发方案的理解 第9题我国油田开发的方针是在一个较长的时期内实现稳产高产。 您的答案：正确 题目分数：0.5 此题得分：0.5 批注：考查对方针的理解 第10题凡是能够储集和渗滤流体的岩层都叫储集层. 您的答案：正确 题目分数：0.5 此题得分：0.5
数据分析软件和工具
以下是我在近三年做各类计量和统计分析过程中感受最深的东西，或能对大家有所帮助。当然，它不是ABC的教程，也不是细致的数据分析方法介绍，它只是“总结”和“体会”。由于我所学所做均甚杂，我也不是学统计、数学出身的，故本文没有主线，只有碎片，且文中内容仅为个人观点，许多论断没有数学证明，望统计、计量大牛轻拍。 于我个人而言，所用的数据分析软件包括EXCEL、SPSS、STATA、EVIEWS。在分析前期可以使用EXCEL进行数据清洗、数据结构调整、复杂的新变量计算（包括逻辑计算）；在后期呈现美观的图表时，它的制图制表功能更是无可取代的利器；但需要说明的是，EXCEL毕竟只是办公软件，它的作用大多局限在对数据本身进行的操作，而非复杂的统计和计量分析，而且，当样本量达到“万”以上级别时，EXCEL的运行速度有时会让人抓狂。 SPSS是擅长于处理截面数据的傻瓜统计软件。首先，它是专业的统计软件，对“万”甚至“十万”样本量级别的数据集都能应付自如；其次，它是统计软件而非专业的计量软件，因此它的强项在于数据清洗、描述统计、假设检验（T、F、卡方、方差齐性、正态性、信效度等检验）、多元统计分析（因子、聚类、判别、偏相关等）和一些常用的计量分析（初、中级计量教科书里提到的计量分析基本都能实现），对于复杂的、前沿的计量分析无能为力；第三，SPSS主要用于分析截面数据，在时序和面板数据处理方面功能了了；最后，SPSS兼容菜单化和编程化操作，是名副其实的傻瓜软件。 STATA与EVIEWS都是我偏好的计量软件。前者完全编程化操作，后者兼容菜单化和编程化操作；虽然两款软件都能做简单的描述统计，但是较之 SPSS差了许多；STATA与EVIEWS都是计量软件，高级的计量分析能够在这两个软件里得到实现；STATA的扩展性较好，我们可以上网找自己需要的命令文件（.ado文件），不断扩展其应用，但EVIEWS 就只能等着软件升级了；另外，对于时序数据的处理，EVIEWS较强。 综上，各款软件有自己的强项和弱项，用什么软件取决于数据本身的属性及分析方法。EXCEL适用于处理小样本数据，SPSS、 STATA、EVIEWS可以处理较大的样本；EXCEL、SPSS适合做数据清洗、新变量计算等分析前准备性工作，而STATA、EVIEWS在这方面较差；制图制表用EXCEL；对截面数据进行统计分析用SPSS，简单的计量分析SPSS、STATA、EVIEWS可以实现，高级的计量分析用 STATA、EVIEWS，时序分析用EVIEWS。 关于因果性 做统计或计量，我认为最难也最头疼的就是进行因果性判断。假如你有A、B两个变量的数据，你怎么知道哪个变量是因（自变量），哪个变量是果（因变量）？ 早期，人们通过观察原因和结果之间的表面联系进行因果推论，比如恒常会合、时间顺序。但是，人们渐渐认识到多次的共同出现和共同缺失可能是因果关系，也可能是由共同的原因或其他因素造成的。从归纳法的角度来说，如果在有A的情形下出现B，没有A的情形下就没有B，那么A很可能是B的原因，但也可能是其他未能预料到的因素在起作用，所以，在进行因果判断时应对大量的事例进行比较，以便提高判断的可靠性。 有两种解决因果问题的方案：统计的解决方案和科学的解决方案。统计的解决方案主要指运用统计和计量回归的方法对微观数据进行分析，比较受干预样本与未接受干预样本在效果指标（因变量）上的差异。需要强调的是，利用截面数据进行统计分析，不论是进行均值比较、频数分析，还是方差分析、相关分析，其结果只是干预与影响效果之间因果关系成立的必要条件而非充分条件。类似的，利用截面数据进行计量回归，所能得到的最多也只是变量间的数量关系；计量模型中哪个变量为因变量哪个变量为自变量，完全出于分析者根据其他考虑进行的预设，与计量分析结果没有关系。总之，回归并不意味着因果关系的成立，因果关系的判定或推断必须依据经过实践检验的相关理论。虽然利用截面数据进行因果判断显得勉强，但如果研究者掌握了时间序列数据，因果判断仍有可为，其
激光加工数值模拟技术研究现状
激光焊接数值模拟技术研究现状 摘要：介绍了激光焊接数值模拟技术在激光焊接温度场分析、激光焊接应力应变分析、激光焊接熔池流动场分析、激光焊接接头微观组织分析方面的研究现状, 并对激光焊接数值模拟技术在这几方面的模拟方法、原理及模型的建立进行了较为详细的介绍。最后, 对我国焊接数值模拟技术的发展进行了展望。 关键词：激光焊接数值模拟温度场应力应变熔池模拟接头 1.引言 激光焊接是利用高能量的激光光束作为热源照射到材料表面从而使材料汽化、熔化并冷却结晶形成焊缝的一种先进焊接方法。由于具有高能量密度，高效率，高精度，柔性好等优点，激光焊接受到了广泛的重视，并且已经应用到了航天航空，汽车制造等材料加工的领域。 随着激光焊接应用的增加，人们对激光焊接过程的研究也更加重视。在计算机仿真技术应用于焊接学科之前，人们为了某些材料制定合适的激光焊接工艺，往往需要进行大量的实验，耗费大量的物力人力财力。因此，建立激光焊接的数学模型并对激光焊接进行全过程的模拟仿真，对于预测焊接结果，实现激光焊接工艺参数预选和优化，减少工艺试验次数，甚至控制激光焊接过程，防止出现焊接缺陷都具有十分重要的意义。 2.焊接数值模拟发展历史 焊接过程的数值模拟研究由来已久。70年代，有限元法逐渐在焊接温度场分析计算中使用。1975年，加拿大的Poley和Hibbert提出利用有限元法研究焊接温度场，并编制了简单的温度场计算程序。1976年，Krutz在博士论文中专门研究了利用焊接温度场预测接头强度的问题[1]。随着80年代末90年代初，热弹塑性计算理论的逐步完善，焊接应力应变的数值模拟也逐渐发展起来并日益成熟。同时，计算机技术的发展也为焊接数值模拟提供了更有力的支持，使人们能够进一步对激光焊接的接头组织，熔池流动等进行更深入的数值模拟研究。
现代数值模拟方法及其应用
现代数值模拟方法及其应用 这是一门什么样的课？ 研究生的全校公选课。 （怎么讲，有待实践和探讨） 假设应当具有的基本知识 高等数学 如微积分、级数展开、微分方程 线形代数、概率统计 问题：关于级数展开及其应用 21 ()(0)'(0)''(0)...2! f x f f x f x =+++ 答： * 当x 较小时，可取前面几项作为函数的近似 * 当函数形式未知时，可用级数逐项逼近 计算机编程 包括 Linux 系统、画图和数据分析软件， 例如 xmgrace ，mitlab 问题：A=0.0D+00 DO 10 I=1,10 A=A+1.0D+00*I 10 CONTINUE 代表什么含义 物理学 （50％内容或多或少与物理学有关） 最理想是学习过普通物理学 或者中学的物理学，能理解基本的物理问题 比如，物理是研究物质的结构和运动的学科 物质有各种形态，如气态、液态和固态等
物质的运动遵从一定的运动规律 如运动方程，分布函数等 问题：力学、统计物理和量子力学的基本知识 化学、生物学和经济学 简单的基本知识 基本的英文阅读和书写能力 不打算非常系统地讲授种种数值模拟方法 因为时间有限、精力有限 重点讲两种方法 Monte Carlo 模拟 和 分子动力学 简单介绍一些重要的基本方法 一定程度上给出数值模拟方法的概况 目的是学习应用计算机模拟方法研究科学问题 至少了解如何用计算机模拟方法研究科学问题 包括 方法本身 科学问题的表述，模型化 Ising 模型的种种应用 {} 1 1 1 i j i j i H k T i S i H K S S h S kT Z S e -- =+==±∑∑∑
油藏数值模拟入门指南
[转]【推荐】油藏数值模拟入门指南 尝试写一写油藏数值模拟入门指南，希望对那些刚刚开始进入油藏数值模拟领域的工作者有所帮助。 第一：从掌握一套商业软件入手。 我给所有预从事油藏数值模拟领域工作的人员第一个建议是先从学一套商业数值模拟软件开始。起点越高越好，也就是说软件功能越强越庞大越好。现在在市场上流通的ECLIPSE,VIP 和CMG都可以。如果先学小软件容易走弯路。有时候掌握一套小软件后再学商业软件会有心里障碍。 对于软件的学习，当然如果能参加软件培训最好。如果没有机会参加培训，这时候你就需要从软件安装时附带的练习做起。油藏数值模拟软件通常分为主模型，数模前处理和数模后处理。主模型是数模的模拟器，即计算部分。这部分是最重要的部分也是最难掌握的部分。它可以细分为黑油模拟器，组分模拟气，热采模拟器，流线法模拟器等。数模前处理是一些为主模拟器做数据准备的模块。比如准备油田的构造模型，属性模型，流体的PVT参数，岩石的相渗曲线和毛管压力参数，油田的生产数据等。数模后处理是显示模拟计算结果以及进行结果分析。 以ECLIPSE软件为例，ECLIPSE100,ECLIPSE300和FrontSim是主模拟器。ECLISPE100是对黑油模型进行计算，ECLISPE300是对组分模型和热采模拟进行计算，FrontSim是流线法模拟器。前处理模块有Flogrid,PVTi,SCAL,Schedule,VFPi等。Flogrid用于为数值模拟建立模拟模型，包括油田构造模型和属性模型；PVTi用于为模拟准备流体的PVT参数,对于黑油模型，主要是流体的属性随地层压力的变化关系表，对于组分模型是状态方程；SCAL为模型准备岩石的相渗曲线和毛管压力输入参数；Schedule处理油田的生产数据，输出ECLIPSE 需要的数据格式（关键字）；VFPi是生成井的垂直管流曲线表，用于模拟井筒管流。ECLIPSE OFFICE和FLOVIZ是后处理模块，进行计算曲线和三维场数据显示和分析，ECLIPSE OFFICE同时也是ECLIPSE的集成平台。 对于初学者，不但要学主模型，也需要学前后处理。对于ECLISPE的初学者，应该先从ECLISPE OFFICE学起，把ECLISPE OFFICE的安装练习做完。然后再去学Flogrid，Schedule 和SCAL。PVTi主要用于组分模型，做黑油模型可以不用。 第二：做油藏数值模拟都需要准备什么参数 在照着软件提供的安装例子做练习时经常遇到的问题是：虽然一步一步按照手册的说明做，但做的时候不明白每一步在做什么，为什么要这么做。这时候的重点在于你要知道你一开始做的工作都是为数值模拟计算提供满足软件格式要求的基础参数。有了这些基础参数你才能开始进行模拟计算。这些基础参数包括以下几个部分： 1。模拟工作的基本信息：设定是进行黑油模拟，还是热采或组分模拟；模拟采用的单位制（米制或英制）；模拟模型大小（你的模型在X,Y,Z三方向的网格数）；模拟模型网格类型（角点网格，矩形网格，径向网格或非结构性网格）；模拟油藏的流体信息（是油，气，水三相还是油水或气水两相，还可以是油或气或水单相，有没有溶解气和挥发油等）;模拟油田投入开发的时间；模拟有没有应用到一些特殊功能（局部网格加密，三次采油，端点标定，多段井等）；模拟计算的解法（全隐式，隐压显饱或自适应）。 2。油藏模型：模型在X,Y,Z三方向的网格尺寸大小，每个网格的顶面深度，厚度，孔隙度，渗透率，净厚度（或净毛比）。网格是死网格还是活网格。断层走向和断层传导率。
地下水数值模拟研究进展和发展趋势
地下水数值模拟研究进展与发展趋势 摘要：地下水数值模拟的应用研究进展国外对地下水数值模拟的研究和应用较早，且理论、技术等各方面相对成熟，目前已经从“水量问题”的应用研究逐步过渡到“水质问题”的应用研究上，以解决各种更复杂的地下水问题。国内相关研究起步较晚、同国外存在一定的差距，主要应用研究在地下水位预测、地下水资源开发利用、地下水循环机制研究、地下水资源预报评价等水量、水位问题方面，但在加油站渗漏场、石油渗漏场、垃圾填埋场、工业废料填埋场、矿区、核废料处置场等污染场地污染物的迁移问题方面的应用研究逐渐增多，并已取得了一定的成果。 关键词：数值模拟、进展、发展趋势 随着计算机技术的快速发展，科学有效的数值计算方法在处理地下水污染、分析地下水资源评估等问题中的应用越来越广泛; 利用数值模拟软件对地下水流等问题进行模拟，以其有效性、灵活性和相对廉价性逐渐成为地下水研究领域的一种不可缺少的重要方法［1］。尤其针对加油站渗漏场、石油渗漏场、垃圾填埋场、工业废料填埋场、矿区、核废料处置场等污染场地污染物的迁移问题，建立准确的数值模型进行预测是查明污染物污染潜水范围、程度及其分布特征最有效最直观的方法之一，同时还可以为污染区实施污染防治与修复等优化配置提供科学技术支持［2］。 地下水数值模拟的应用研究进展国外对地下水数值模拟的研究和应用较早，且理论、技术等各方面相对成熟，目前已经从“水量问题”的应用研究逐步过渡到“水质问题”的应用研究上，以解决各种更复杂的地下水问题。国内相关研究起步较晚、同国外存在一定的差距，主要应用研究在地下水位预测、地下水资源开发利用、地下水循环机制研究、地下水资源预报评价等水量、水位问题方面，但在加油站渗漏场、石油渗漏场、垃圾填埋场、工业废料填埋场、矿区、核废料处置场等污染场地污染物的迁移问题方面的应用研究逐渐增多，并已取得了一定的成果［4］。 近几十年来，随着地下水科学和计算机科学的发展，地下水数值模拟也得到了快速发展，主要体现在：加拿大Borden基地、美国Cape Cod基地与Columbus基地开展的大型野外试验场研究，大大丰富了地下水溶质运移的理论和方法，取得不少新的认识，并为发展和检验溶质运移理论和相应数学模型提供了大量数据（MacKay et al，1986； LeBlanc et al，1991；Bogga et al，1992；Zheng and Gorelick，2003）；随机方法在非均质介质渗流和溶质运移的模拟中得到比较多的应用，从而加深、甚至改变了人们对此类介质中流体运动和溶质运移的认识（Dagan and Neuman，1997； Zhang D，2002）；通过多孔介质中水流运动、溶质运移和化学反应，甚至生物过程的耦合建立模型来集成地研究这些过程也取得很多进展（van Genuchten and Sudicky，1999； Yeh and Tripathi，1989； Barry et al，2002）。此外，计算方法也取得不少进展，但溶质运移模拟中数值弥散和振荡问题的解决和地下水模拟逆问题的求解进展比较缓慢（Sun and Yeh，2007）。 由于种种原因，国内地下水数值模拟开展得比较晚，始于20世纪70年代初，当时文化大革命还没有结束，所以从事这项工作困难重重，而且人也不多，主要来自高等学校和研究部门，以后才逐步扩展到产业部门。为了加快我国地下水数值模拟的发展，深切感到有必要
湍流大涡数值模拟进展
第22卷第2期空气动力学学报Vol.22,No.2 2004年06月ACTA AERODYNAMICA SINICA Jun.,2004 文章编号:0258-1825(2004)02-0121-09 湍流大涡数值模拟进展 崔桂香,许春晓,张兆顺 (清华大学工程力学系,北京100084) 摘要:本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题。文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径。文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向。 关键词:湍流;大涡数值模拟;亚格子模型;近壁模型;标量湍流 中图分类号:V211.3文献标识码:A * 0引言 复杂流动的准确数值预测是当前航空、航天器研究和设计中迫切需要解决的空气动力学前沿问题之一。随着计算空气动力学方法的不断完善,计算机品质的不断提高,湍流的数值模拟方法成为提高数值预测航天器空气动力特性的瓶颈。 目前,数值预测湍流的方法有:直接数值模拟(DNS),大涡数值模拟(LE S)和雷诺平均模拟(RANS)。 直接数值模拟(DNS)是精确数值模拟湍流的方法,它的主要困难在于湍流是一种不规则多尺度运动,无论在空间上或者时间上湍流都有十分宽广谱。准确数值模拟湍流既要精确计算大尺度流动;又要足够准确地计算小尺度运动。在最简单的各向同性湍流中湍流的空间尺度有以下估计:L ma x/l mi n~Re3/4K,同样最大和最小时间尺度之比T max/t min~Re3/4K,它们都和流动的泰勒雷诺数Re K有关。按照上述估计,空间网格数至少应有:N=N x@N y@N z~Re9/4K,运算量超过Re3K,航空航天器复杂绕流计算的网格数和运算时间远远超过上述估计。因此,目前不具备直接数值模拟复杂工程湍流需要的计算机,湍流直接数值模拟只能作为低雷诺数简单湍流的研究工具。 工程中常用的复杂湍流数值模拟方法是求解雷诺平均的控制方程,这种方法只计算大尺度平均流动,而所有湍流脉动对平均流动的作用,即雷诺应力,用模型假设封闭。由于雷诺应力主要由大尺度脉动贡献,而大尺度脉动和流动的几何特性密切相关,因此雷诺平均模式不是普适的,而是和流动有关,就是说,不存在对一切流动都适用的统一模式;对于不同类型的流动,模式的形式或系数需要修正,而这种修正常常带有经验性。所以,雷诺平均模式不是理想的封闭方法。 湍流大涡数值模拟是有别于直接数值模拟和雷诺平均模式的一种数值预测湍流的方法。这种方法是基于对各种尺度湍流脉动在输运和耗散中作用的认识:大尺度湍流脉动具有主要的能量和动量并支配湍流脉动的动量和能量输运;而湍动能的耗散主要发生在小尺度脉动中[1];根据这一认识产生了湍流大涡数值模拟。它的具体实施方法如下:首先,用滤波方法将小尺度脉动从湍流脉动中去掉,假设空间任意一点的滤波函数为G(x-x0),最简单的滤波器是盒式滤波: G(G)=1,|G i|[$/2, G(G)=0,|G i|>$/2(1)利用滤波器对湍流速度场过滤,过滤后的速度脉动中不存在过滤尺度$以下的脉动成分,称为可解湍流: u i(x,t)= 1 $3 m D u(y,t)G(x-y)d y(2) *收稿日期:2003-03-26;修订日期:2003-06-02. 基金项目:国家自然科学基金资助项目(批准号:10272065,10232020). 作者简介:崔桂香(1950-),女,清华大学工程力学系教授,从事湍流大涡数值模拟和标量湍流的研究.