实验报告一
（四学时）
一【实验目的】
1.了解使用VC++6.0或CodeBlocks的开发环境，掌握在编译器上如何编辑、编译、连接
和运行一个C程序；
2.通过运行简单的C程序，初步了解C源程序的特点；
3.学习变量常量的使用；
4.掌握顺序结构程序设计方法；
二【实验内容】
第一部分：编译环境熟悉及常量变量的使用（2学时）
【实验1—基础题】
要求：完成实验指导书P208页2.2.1第2章实验题：熟悉上机实验环境和不同类型变量的定义与赋值，并把实验结果（包括编写的程序和运行结果）截图粘贴到“实验结果”下方（其他试验项目要求相同，不再重复说明）
目的：熟悉VC++6.0编译环境，了解变量类型和赋值方法
实验结果：
【实验2—验证题】
要求：完成实验指导书P208页2.2.2第3章实验题：编写一个程序从键盘输入球的半径r，计算球的体积和表面积并输出，把实验结果截图到相应位置（使用宏定义或const常量定义圆周率）
目的：熟悉简单的算术运算
实验结果：
第二部分：顺序结构程序设计（2学时）
【实验—3】数学函数及运算表达式的使用
编写程序求函数的值，要求输入一个x，程序输出y的值。
程序分析：
（1）声明两个实型变量x、y
（2）用scanf函数从键盘输入一个数给x
（3）带入表达式求出y 的值并用printf函数输出
注意事项：
（1）头文件还应包含数学头文件：#include
sin(67.0/180*3.14159)，数学公式具体参照实验指导书附录
实验结果：
【实验—4】顺序结构程序设计应用
美国人习惯使用英制计量单位，他们用几尺几寸（英尺英寸）来报自己的身高。如果遇到?
个美国人告诉你他的身高是5英尺7英寸，他的身高应该是多少米呢？
编程实现一个英尺转换器，输入英尺英寸，输出对应米数。
提示：换算关系如下：
1 英尺= 1
2 英寸=30.48厘米
实验结果：
【实验5—延伸训练】算数运算符的使用
要求：编写一个程序，其功能是：将两个两位数的正整数a、b合并成一个整数放在c中。合并的方式是：将a中的十位和个位数一次放在c数的千位和十位上，b数的十位和个位数依次放在c数的百位和个位上。例如：a=45,b=12时，执行程序后，c=4152.
实验结果：
选做题：
考试目标及考试大纲 本题库的编纂目的旨在给出多套试题，每套试题的考查范围及难度配置均基于“水平测试”原则，按照教学大纲和教学内容的要求，通过对每套试题的解答，可以客观公正的评定出学生对本课程理论体系和应用方法等主要内容的掌握水平。通过它可以有效鉴别和分离不同层次的学习水平，从而可以对学生的学习成绩给出客观的综合评定结果。 本题库力求作到能够较为全面的覆盖教学内容，同时突显对重点概念、重点内容和重要方法的考查。考试内容包括以下部分： 绪论与误差：绝对误差与相对误差、有效数字、误差传播分析的全微分法、相对误差估计的条件数方法、数值运算的若干原则、数值稳定的算法、常用数值稳定技术。 非线性方程求解：方程的近似解之二分法、迭代法全局收敛性和局部收敛定理、迭代法误差的事前估计法和事后估计法、迭代过程的收敛速度、r 阶收敛定理、Aitken加速法、Ne w to n法与弦截法、牛顿局部收敛性、Ne w to n收敛的充分条件、单双点割线法（弦截法）、重根加速收敛法。 解线性方程组的直接法：高斯消元法极其充分条件、全主元消去法、列主元消去法、高斯-若当消元法、求逆阵、各种消元运算的数量级估计与比较、矩阵三角分解法、Doolittle 和Crout三角分解的充分条件、分解法的手工操作、平方根法、Cholesky分解、改进的平方根法(免去开方)、可追赶的充分条件及适用范围、计算复杂性比较、严格对角占优阵。 解线性方程组迭代法：向量和矩阵的范数、常用向量范数的计算、范数的等价性、矩阵的相容范数、诱导范数、常用范数的计算；方程组的性态和条件数、基于条件数误差估计与迭代精度改善方法；雅可比（Jacobi）迭代法、Gauss-Seidel迭代法、迭代收敛与谱半径的关系、谱判别法、基于范数的迭代判敛法和误差估计、迭代法误差的事前估计法和事后估计法；严格对角占优阵迭代收敛的有关结论；松弛法及其迭代判敛法。 插值法：插值问题和插值法概念、插值多项式的存在性和唯一性、插值余项定理；Lagrange插值多项式；差商的概念和性质、差商与导数之间的关系、差商表的计算、牛顿（Newton）插值多项式；差分、差分表、等距节点插值公式；Hermite插值及其插值基函数、误差估计、插值龙格（Runge）现象；分段线性插值、分段抛物插值、分段插值的余项及收敛性和稳定性；样条曲线与样条函数、三次样条插值函数的三转角法和三弯矩法。 曲线拟合和函数逼近：最小二乘法原理和多项式拟合、函数线性无关概念、法方程有唯一解的条件、一般最小二乘法问题、最小二乘拟合函数定理、可化为线性拟合问题的常见函数类；正交多项式曲线拟合、离散正交多项式的三项递推法。最佳一致逼近问题、最佳一致逼近多项式、切比雪夫多项式、切比雪夫最小偏差定理、切比雪夫多项式的应用(插值余项近似极小化、多项式降幂)。本段加黑斜体内容理论推导可以淡化，但概念需要理解。 数值积分与微分：求积公式代数精度、代数精度的简单判法、插值型求积公式、插值型求积公式的代数精度；牛顿一柯特斯(Newton-Cotes)公式、辛卜生（Simpson）公式、几种低价牛顿一柯特斯求积公式的余项；牛顿一柯特斯公式的和收敛性、复化梯形公式及其截断误差、复化Simpson公式及其截断误差、龙贝格（Romberg）求积法、外推加速法、高斯型求积公式、插值型求积公式的最高代数精度、高斯点的充分必要条件。正交多项式的构造方法、高斯公式权系数的建立、Gauss-Legendre公式的节点和系数。本段加黑斜体内容理论推导可以淡化，但概念需要理解。 常微分方程数值解：常微分方程初值问题数值解法之欧拉及其改进法、龙格—库塔法、阿当姆斯方法。
西南石油大学油层物理考研《油层物理CAI课件》练习检测题 目 9年第号 1。2.3.4.5.6.7. 储层流体的物理性质 天然气有哪些分类？它是如何分类的？代表天然气成分的一般方法是什么？压缩因子z的物理意义是什么？ 如何确定多组分物质的表观临界压力和表观临界温度？天然气的体积系数是多少？天然气的可压缩性是多少？什么是泡点和露点？地层油的饱和压力是多少？如何定义 地层油的溶解油气比？试分析它与天然气在原油中溶解度的区别和联系。 8。一个带刻度的活塞气缸配有45000cm3(在10325帕，00摄氏度时)。当温度和压力从 03 变化到地层条件(17.8兆帕，71C)时，测量体积为265厘米。气体压缩系数是多少？9.当天然气的相对密度为0.743，局部地层压力为13.6兆帕，地层温度为93.30℃时，计算天然气的压缩系数。 10。天然气成分分析结果见表1-4。地层压力为8.3兆帕，地层温度为320℃表1-4由CH4 C2H6 C3H8 IC4H10Mol组成，分为0.902 0.045 0.031 0.021 (1)，并得到天然气的压缩系数。(2)计算天然气的体积系数；(3)尝试将地下10000m3天然气所占的体积转换；
11。画出油层的烃相状态，在图上标出纯油层、饱和油层、凝析气藏和气藏的位置，注意压力和温度的范围 12。什么是一次脱气和多次脱气？一次脱气和多次脱气的区别和联系是什么？13.第一油层温度为750℃，饱和压力Pb=18MPa，饱和压力下溶解汽油比Rob=120 m3/m3当压力降至15兆帕时，气体已经分离，气体的相对密度为0.7，r为115。m3/m3，Pb=1.25，两相的体积系数是多少？40.地层水的分析结果如下。试着计算它的水类型。1-10钠+镁+钙+氯-硫酸-HCO 3-4952 838 620 10402 961 187 14。定性绘制天然气的体积系数Bg、压缩系数Cg和粘度系数Pg如下图所示，以及温度的变化规律 15。试画出原油体积系数B0、两相体积系数U、压缩系数Co和粘度μ0随压力和温度的定性变化规律如下(如果有饱和压力，应该注意) 16。定性图显示了不同温度下天然气在原油中溶解度随压力的变化规律。17.利用物质平衡原理，推导出以下条件下的储层储量计算公式。 油藏的原始状态是溶解气驱油藏，没有气顶和边水。经过一段时间的开发，地层压力降至饱和压力以下，形成二次气顶，但储层的孔隙体积不变，如图7所示。 图7开发前后储层流体状态变化示意图 集:原始储油量为Ni (S？M3) 累积石油产量为Np (S？M3)溶解气油比？M3/米)累计平均生产气
一、选择题（共10小题，每小题3分，共30分，将答案填入题后方括号内） 1、（本题3分） 质量分别为m 1和m 2的两滑块A 和B 通过一轻弹簧水平连结后置于水平桌面上，滑块与桌面间的摩擦系数均为μ，系统在水平拉力F 作用下匀速运动，如图所示．如突然撤消拉力，则刚撤消后瞬间，二者的加速度a A 和a B 分别为 (A) a A =0 , a B =0. (B) a A >0 , a B <0. (C) a A <0 , a B >0. (D) a A <0 , a B =0. ［ ］ 2、（本题3分） 一质量为m 的质点，在半径为R 的半球形容器中，由静止开始自边缘上的A 点滑下，到达最低点B 时，它对容器的正压力为N ．则质 点自A 滑到B 的过程中，摩擦力对其作的功为 (A) )3(21 mg N R -． (B) )3(21 N mg R -． (C) )(2 1 mg N R -． (D) )2(2 1 mg N R -． ［ ］ 3、（本题3分） 一质量为M 的弹簧振子，水平放置且静止在平衡位置，如图所示．一质量为m 的子弹以水平速度v 射入振子中，并随之一起运动．如果水平面光滑，此后弹簧的最大势能为 (A) 2 2 1v m ． (B) )(222m M m +v ． (C) 22 22)(v M m m M +． (D) 2 22v M m ． ［ ］ 4、（本题3分） 已知氢气与氧气的温度相同，请判断下列说法哪个正确？ x A B
(A) 氧分子的质量比氢分子大，所以氧气的压强一定大于氢气的压强． (B) 氧分子的质量比氢分子大，所以氧气的密度一定大于氢气的密度． (C) 氧分子的质量比氢分子大，所以氢分子的速率一定比氧分子的速率大. (D) 氧分子的质量比氢分子大，所以氢分子的方均根速率一定比氧分子的方均根速率大． ［ ］ 5、（本题3分） 关于热功转换和热量传递过程，有下面一些叙述： (1) 功可以完全变为热量，而热量不能完全变为功； (2) 一切热机的效率都只能够小于1； (3) 热量不能从低温物体向高温物体传递； (4) 热量从高温物体向低温物体传递是不可逆的． 以上这些叙述 (A) 只有(2)、(4)正确． (B) 只有(2)、(3) 、(4)正确． (C) 只有(1)、(3) 、(4)正确． (D) 全部正确． ［ ］ 6、（本题3分） 已知某简谐振动的振动曲线如图所示，位移的单位为厘米，时间单位为秒．则此简谐振动的振动方程为： (A) )3 232cos(2π+π=t x ． (B) )3 232cos(2π-π=t x ． (C) )3 234cos(2π+π=t x ． (D) )3234cos(2π-π=t x ． (E) )4 134cos(2π-π=t x ． ［ ］
西南石油大学本科生毕业设计（论文）工作管理规定 西南石大教【2015】115号 第一章总则 第一条毕业设计和毕业论文（以下称毕业设计（论文））是高等学校实现培养目标的重要教学环节，是对学生学习成果的全面检验。做好学生的毕业设计（论文）工作，对全面衡量和提高教学质量具有重要的意义。为进一步规范毕业设计（论文）工作，提高毕业设计（论文）质量，特制定本规定。 第二条本科生毕业设计（论文）的基本目标是培养学生综合运用所学基础理论和专业知识分析问题、解决问题和初步进行科学研究的能力，并培养其优良的思想品质和探求真理的科学精神，提高其综合素质。其培养和训练过程应包括综合运用知识能力、文献资料应用能力、设计（实验）能力、计算能力、外文应用能力、计算机应用能力、技术经济与社会文化分析能力、创新创业能力、团队协作能力等内容，并根据专业性质有所侧重。 第三条学生用于毕业设计（论文）的时间（包括评阅和答辩时间）不得少于专业人才培养方案规定的周数。毕业设计（论文）工作的主要环节包括：申报题目——选题——开题——完成毕业设计（论文）——审阅——评阅——答辩——成绩评定——总结归档。 第二章组织与管理 第四条毕业设计（论文）工作在学校校长的统一领导下实行校、院两级管理。学校校长和学院院长是毕业设计（论文）教学工作与教学质量的第一责任人，分管教学工作的副校长负责领导全校毕业设计（论文）工作，教务处代表学校进行统一管理与协调，学院具体组织，专业教研室负责实施。 1．学校成立毕业设计（论文）工作指导委员会，负责全校毕业设计（论文）工作的指导、协调与宏观管理。指导委员会由分管教学工作的副校长、分管学生工作的副书记、教务处处长及分管副处长，以及学生工作部（处）、规划与评估处、实验室与设备管理处、财务处、现代教育技术中心和图书馆等部门负责人组成。 2．学院成立毕业设计（论文）工作领导小组，负责组织和领导本学院毕业设计（论文）工作。领导小组由院长、分管教学工作的副院长、学院党总支副书记、学院教学办公室主任、教研室主任、学生工作辅导员组成。各学院院长对本学院毕业设计（论文）教学工作负责，分管教学工作的副院长负责具体组织领导。
西南石油大学本科生毕业设计（论文）撰写规范 毕业设计（论文）是学生综合运用所学知识和技能，学习科学研究或工程设计基本方法，培养实践能力、创新能力和科学精神的重要实践教学环节。为规范我校本科生毕业设计(论文)的撰写工作，根据《中华人民共和国国家标准科学技术报告、学位论文和学术论文的编写格式》（国家标准GB7713-87）的规定，制定“西南石油大学本科生毕业设计(论文)撰写规范”。 一、毕业设计（论文）结构要求 毕业设计（论文）应包括： 1．论文题目 题目应该简短、明确、有概括性。通过题目，能大致了解论文内容、专业特点和学科范畴。字数要适当，中文题目一般不宜超过20字，英文题目一般不宜超过10个实词。必要时可加副标题。 2．中英文摘要与关键词 中英文摘要：摘要应概括反映出毕业设计(论文)的内容、方法、成果和结论。摘要中不宜使用公式、图表，不标注引用文献编号。中文摘要一般500字左右，英文摘要内容应与中文摘要相对应。为了便于文献检索，应在摘要下方另起一行注明论文的关键词。 关键词：关键词是供检索用的主题词条，应采用能覆盖设计（论文）主要内容的通用技术词条（参照相应的技术术语标准）。关键词一般为3～5个，按词条的外延层次排列（外延大的排在前面）。 3．目录 目录按三级标题编写（例如：1……、1.1……、1.1.1……），标题应层次清晰。目录中的标题要与正文中的标题一致。 4．正文 正文是毕业设计(论文)的主体和核心部分，一般应包括绪论、论文主体及结论等部分。 绪论 绪论一般是毕业设计(论文)主体的开端篇。包括：毕业设计（论文）的背景及目的；国内外研究状况和相关领域中已有的成果；设计和研究方法；设计过程及研究内容等。
一. 填空 2．Gauss型求积公式不是插值型求积公式。（限填“是”或“不是”） 3. 设l k(x)是关于互异节点x0, x1,…, x n, 的Lagrange 插值基函数，则 0 m=1,2,…,n 5．用个不同节点作不超过次的多项式插值，分别采用Lagrange插值方法与Newton插值方法所得多项式相等(相等, 不相等)。 。 7. n个不同节点的插值型求积公式的代数精度一定会超过n-1次 8.f(x)=ax7+x4+3x+1,f[20, 21,…,27]= a，f [20, 21,…,28]= 0 10设 (i=0,1,…,n)，则＝ _x_ , 这里（x i x j,ij, n2）11.设称为柯特斯系数 则=______1____ 12采用正交多项式拟合可避免最小二乘或最佳平方逼近中常见的_法方程组病态___问题。 13辛卜生（Simpson）公式具有___3____次代数精度。 14 牛顿插商与导数之间的关系式为： 15试确定[0,1]区间上2x3的不超过二次的最佳一致逼近多项式p(x), 该多项式唯一否？答：p(x)=(3/2)x, ; 唯一。 17.给定方程组记此方程组的Jacobi迭代矩阵为B J=(a ij)33，则a23= -1； ,且相应的Jacobi迭代序列是__发散_____的。 18.欧拉预报--校正公式求解初值问题的迭代格式(步长为h) ，此方法是阶方法。 ，此方法是 2阶方法。 19. 2n阶Newton-Cotes公式至少具有2n+1次代数精度。 20.设，则关于的 ||f|| =1 21矩阵的LU分解中L是一个 _为单位下三角阵，而U是一个上三角阵____。 22.设y=f (x1,x2) 若x1,x2,的近似值分别为x1*, x2*,令y*=f(x1*,x2*)作为y的近似值,其绝对误差限的估计式为: ||f(x1*,x2*)|x1-x*1|+ |f(x1*,x2*)|x2- x*2| 23设迭代函数(x)在x*邻近有r（1）阶连续导数，且x* = (x*)，并且有(k) (x*)=0 (k=1,…,r-1)，但(r) (x*)0，则x n+1=(x n)产生的序列{ x n }的收敛阶数为___r___ 24设公式为插值型求积公式，则, 且=b-a 25称微分方程的某种数值解法为p阶方法指的是其局部截断误差 为O(h p+1)。 26.设x0, x1,x2是区间[a, b]上的互异节点，f(x)在[a, b]上具有各阶导数，过
1-1 何谓采油指数的物理意义？如何获取？影响单相渗流和油气两相渗流采油指数的主要因 素有何异同？ 1-2 已知A井位于面积4.5×104m2的正方形泄流区域中心，井眼半径r w =0.1m,根据高压物性资料 B O =1.15，μ O 为4mPa.s；由压力恢复试井资料获得S=3。试根据下表中测试资料绘制IPR曲线， 并求采油指数J O 及油层参数Kh。 1-3 0， 1-4 1-5 ， 1-6 v ）。 1-7 80m3/d，1-8 1-10 p =6.4mm, 、20 1-15 、0.65和℃/100m。 4-1 ％；产液指数 0.6；油相对密度0.86；气体偏差系数0.86。设计和选择完整的潜油电泵系统。 4-2 设计潜油电泵时为什么要进行粘度校正? 4-3 潜油电泵井中，为什么井液必须通过电机? 4-4 潜油电泵井中，为什么采用高效率的井下分离器更加优越? 4-5 设计和选择螺杆泵系统，已知条件与习题1相同。 4-6 试述螺杆泵的采油装置、采油原理和特点。 4-7 设计和选择完整的水力射流泵系统，已知条件与习题1相同。
下列计算题已知基本数据如下： 油井数据：垂直井，井距300m×300m，井深2000m；油管内、外径分别为Ф62mm和Ф73mm，套管内、外径分别为Ф121.36mm和Ф139.7mm；射孔密度16孔/m，孔眼直径10mm。 地层数据：上覆层岩石平均密度2300kg/m3，地层压力系数1.0，岩石弹性模量20000MPa、泊松比0.20，岩石抗张强度2.5MPa，地层温度70°C，孔隙弹性常数取1.0；产层有效 厚度15m，储层渗透率2×10-3μm2，孔隙度12%；原油饱和度60%，压缩系数1.7×10-2压裂液性能：牛顿流体密度1020kg/m3，粘度80mPa.s，减阻率60%；初滤失系数5×10-4m3/m2， 压力校正后造壁性滤失系数8×10-4m/m in。 支撑剂性能：φ0.45~0.9mm宜兴陶粒，颗粒密度2800 kg/m3，砂堆孔隙度35%。 生产限制条件：生产时井底流动压力10MPa，最大排量3.0m3/min，最大砂浓度720 kg/m3。6-1 根据低渗透储层中水力压裂泵压变化典型示意曲线说明地层破裂、裂缝延伸和闭合特征。 6-2 6-3 6-10 6-11 6-12 6-15 取
用复化梯形公式计算积分 1 ()f x dx ?,要把区间[0,1]一般要等分 41 份才能保 证满足误差小于0.00005的要求（这里(2) () 1f x ∞ ≤） ；如果知道(2) ()0f x >，则 用复化梯形公式计算积分1 ()f x dx ? 此实际值 大 (大，小)。 在以1 0((),())()(),(),()[0,1]g x f x xf x g x dx f x g x C = ∈?为内积的空间C[0,1] 中，与非零常数正交的最高项系数为1的一次多项式是 2 3 x - 3. (15分)导出用Euler 法求解 (0)1y y y λ'=??=? 的公式, 并证明它收敛于初值问题的精确解 解 Euler 公式 11,1,,,k k k x y y h y k n h n λ--=+== L -----------(5分) ()()1011k k k y h y h y λλ-=+==+L ------------------- (10分) 若用复化梯形求积公式计算积分1 x I e dx = ? 区间[0,1]应分 2129 等分，即要 计算个 2130 点的函数值才能使截断误差不超过 71 102 -?；若改用复化Simpson 公式，要达到同样精度区间[0,1]应分12 等分，即要计算个 25 点的函数值 1．用Romberg 法计算积分 2 3 2 x e dx -? 解 []02()()2b a T f a f b -= += 9.6410430E-003 10221()222 b a a b T T f -+=+= 5.1319070E-003 10 022243 T T S -= = 4.6288616E-003 22T = 4.4998E-003 21 122243 T T S -= = 4.E-003 10 02221615 S S C -= = 4.6588636E-003 32T = 4.7817699E-003 32 222243 T T S -= = 4.1067038E-003
第一章 储层岩石的物理特性 24、下图1-1为两岩样的粒度组成累积分布曲线，请画出与之对应的粒度组成分布曲线，标明坐标并对曲线加以定性分析。 ∑Log d i W Wi 图1-1 两岩样的粒度组成累积分布曲线 答：粒度组成分布曲线表示了各种粒径的颗粒所占的百分数，可用它来确定任一粒级在岩石中的含量。曲线尖峰越高，说明该岩石以某一粒径颗粒为主，即岩石粒度组成越均匀；曲线尖峰越靠右，说明岩石颗粒越粗。一般储油砂岩颗粒的大小均在1～0.01mm 之间。 粒度组成累积分布曲线也能较直观地表示出岩石粒度组成的均匀程度。上升段直线越陡，则说明岩石越均匀。该曲线最大的用处是可以根据曲线上的一些特征点来求得不同粒度属性的粒度参数，进而可定量描述岩石粒度组成的均匀性。 曲线A 基本成直线型，说明每种直径的颗粒相互持平，岩石颗粒分布不均匀；曲线B 上升段直线叫陡，则可看出曲线B 所代表的岩石颗粒分布较均匀。 30、度的一般变化范围是多少，Φa 、Φe 、Φf 的关系怎样？常用测定孔隙度的方 法有哪些？影响孔隙度大小的因素有哪些？ 答：1）根据我国各油气田的统计资料，实际储油气层储集岩的孔隙度范围大致为：致密砂岩孔隙度自＜1%～10%；致密碳酸盐岩孔隙度自＜1%～5%；中等砂岩孔隙度自10%～20%；中等碳酸盐岩孔隙度自5%～10%；好的砂岩孔隙度自20%～35%；好的碳酸盐岩孔隙度自10%～20%。 2）由绝对孔隙度a φ、有效孔隙度e φ及流动孔隙度ff φ的定义可知：它们之间的关系应该是a φ>e φ>ff φ。 3）岩石孔隙度的测定方法有实验室内直接测定法和以各种测井方法为基础的间接测定法两类。间接测定法影响因素多，误差较大。实验室内通过常规岩心
数值分析上机报告 班级：20级学隧2班 姓名：000000000 学号：00000000000
目录 1 序言 (6) 2 题目 (7) 2.1 题2 (7) 2.1.1 题目内容 (7) 2.1.2 MATLAB程序 (8) 2.1.3 计算结果 (8) 2.1.4 图形 (9) 2.1.5 分析 (14) 2.2 题3 (14) 2.2.1 题目内容 (14) 2.2.2 程序 (14) 2.2.3 计算结果 (14) 2.2.4 图形 (15) 2.2.5 分析 (16) 2.3 选做题5 (16) 2.3.1方法介绍 (17) 2.3.2计算结果及分析 (17) 3总结 (18) 4．附录 (19) 4.1 题1程序代码 (19) 4.2 题2程序代码 (22) 4.3 题3程序代码 (26)
数值分析2015上机实习报告要求 1．应提交一份完整的实习报告。具体要求如下： （1）报告要排版，美观漂亮（若是纸质要有封面，封面上）要标明姓名、学号、专业和联系电话； （2）要有序言，说明所用语言及简要优、特点，说明选用的考量； （3）要有目录，指明题目、程序、计算结果，图标和分析等内容所在位置，作到信息简明而完全； （4）要有总结，全方位总结机编程计算的心得体会； （5）尽量使报告清晰明了，一般可将计算结果、图表及对比分析放在前面，程序清单作为附录放在后面，程序中关键部分要有中文说明或标注， 指明该部分的功能和作用。 2．程序需完好保存到期末考试后的一个星期，以便老师索取用于验证、询问或质疑部分内容。 3．认真完成实验内容，可以达到既学习计算方法又提高计算能力的目的，还可以切身体会书本内容之精妙所在，期间可以得到很多乐趣。 4．拷贝或抄袭他人结果是不良行为，将视为不合格。 5．请按任课老师要求的时间和载体（电子或纸质）提交给任课老师。
习题16 16-6在均匀密绕的螺绕环导线内通有电流20A ，环上线圈 400匝，细环的平均周长是40cm ，测得环内磁感应强度是1.0T 。求： (1)磁场强度； (2)磁化强度； (3)磁化率； (4)磁化面电流的大小和相对磁导率。 [解] (1) 螺绕环内磁场强度 由nI d L =??l H 得 1 -42 m 100.2104020400??=??== -A L nI H (2) 螺绕环内介质的磁化强度 由M B H -= μ得 1-547 m 1076.710210 40 .1??=?-?= -= --A H B M πμ (3) 磁介质的磁化率 由H M m χ=得 8.381021076.74 5 m =??==H M χ (4)环状磁介质表面磁化面电流密度 -15m 1076.7??==A M j 总磁化面电流 A L j dL M I L 55101.34.01076.7?=??=?=?='? 相对磁导率 8.398.3811m 0r =+=+== χμμH B
16-7．一绝对磁导率为μ1的无限长圆柱形直导线，半径为R 1，其中均匀地通有电流I 。导线外包一层绝对磁导率为μ2的圆筒形不导电磁介质，外半径为R 2，如习题16-7图所示。试求磁场强度和磁感应强度的分布，并画出H -r ，B-r 曲线。 [解] 将安培环路定理∑?=?I d L l H 应用于半径为r 的同心圆周 当0≤r ≤1R 时，有 2 2 1 12r R I r H πππ?= ? 所以 2 112R Ir H π= 2111 112R Ir H B πμμ== 当r ≥1R 时，有I r H =?π22 所以r I H π22= 在磁介质内部1R ≤r ≤2R 时，r I H B πμμ22222== 在磁介质外部r ≥2R 时，r I H B πμμ20202 ==' 各区域中磁场强度与磁感应强度的方向均与导体圆柱中电流的方向成右手螺旋关系。 H -r 曲线 B-r 曲线 习题16-7图 R 1 R 2 本图中假设 B 2 12 1μμ>r r 1
油层物理复习重点 一、名词解释：7个，21分， 二、按题意完成：5个，42分， 三、计算题：3个，37分，4-5分8-9分20几分（多步完成，按步给分） 第一章 1.粒度组成概念，主要分析方法，粒度曲线的用途 2.比面概念，物理意义 3.空隙分类（大小，连通性，有效性；毛细管空隙，超毛细管空隙，微毛细管空隙）， 孔隙度概念（绝对孔隙度，有效孔隙度，流动孔隙度，连通孔隙度的概念与区别）， 孔隙度的测定（给定参数会计算，不要求测定的具体步骤） 4.岩石压缩系数及其含义，地层综合弹性压缩系数，弹性驱油量的计算 5.流体饱和度的概念（落实到具体的物质，油、水、气；初始含油、水、气饱和度，残余流体饱和度的概念，束缚水饱和度） 饱和度测定（各种饱和度，会根据给定参数计算） 7.达西定律，及达西公式的物理意义，岩石绝对渗透率感念，液测、气测渗透率的计算方法，液测气测渗透率与岩石绝对渗透率的关系，根据达西定律测定岩石渗透率要满足的三个测定条件，气体滑脱效应对气测渗透率的影响，及影响滑脱效应的因素。 8.胶结概念与类型，
粘土矿物：水敏，酸敏，速敏等，会判断具体的矿物如蒙脱石，高岭石，绿泥石 第二章 1.烃类体系P-T相图，划分相区，临界点，临界凝析温度，临界凝析压力，露点线，泡点线，等液量线，等温反凝析区等术语，露点，泡点，露点压力和泡点压力概念，等温反凝析概念，反凝析作用，对凝析气藏开发的影响，用相图判断油气藏类型。 （露点概念：气相体系生出第一滴液滴时的温度压力点；露点压力：气相体系生出第一滴液滴时的压力） 2.油气分离的两种方式，特点及其结果的差异，以及产生差异的原因，天然气分子量概念，天然气在原油中的溶解规律 3.油气高压物性参数的概念，高压物性参数随压力的变化关系，（肯定会考曲线；不考随温度的变化） 4.平衡常数概念（哪两个之间的平衡关系，） 相平衡中的一些平衡关系（物质平衡，相平衡） 第三章 1.界面张力的概念，界面吸附的两种类型 2.润湿接触角概念，润湿程度判定参数、方法（常用接触角），润湿滞后概念，前进角，后退角概念，润湿滞后对水驱油得影响。 3.油藏润湿性类型，油藏润湿性的影响因素 4毛细管压力概念，毛细管中液体上升高度计算，毛细管滞后，吸入和驱替过程等概念（毛管力是动力，阻力），毛细管压力曲线的测定方法（3种），
数值分析2018-2019第1学期上机实习题 f x，隔根第1题.给出牛顿法求函数零点的程序。调用条件：输入函数表达式() a b，输出结果：零点的值x和精度e，试取函数 区间[,] ，用牛顿法计算附近的根，判断相应的收敛速度，并给出数学解释。 1.1程序代码： f=input('输入函数表达式：y=','s'); a=input('输入迭代初始值：a='); delta=input('输入截止误差：delta='); f=sym(f); f_=diff(f); %求导 f=inline(f); f_=inline(f_); c0=a; c=c0-f(c0)/f_(c0); n=1; while abs(c-c0)>delta c0=c; c=c0-f(c0)/f_(c0); n=n+1; end err=abs(c-c0); yc=f(c); disp(strcat('用牛顿法求得零点为',num2str(c))); disp(strcat('迭代次数为',num2str(n))); disp(strcat('精度为',num2str(err))); 1.2运行结果： run('H:\Adocument\matlab\1牛顿迭代法求零点\newtondiedai.m') 输入函数表达式：y=x^4-1.4*x^3-0.48*x^2+1.408*x-0.512 输入迭代初始值：a=1 输入截止误差：delta=0.0005 用牛顿法求得零点为0.80072 迭代次数为14 精度为0.00036062 牛顿迭代法通过一系列的迭代操作使得到的结果不断逼近方程的实根，给定一个初值，每经过一次牛顿迭代，曲线上一点的切线与x轴交点就会在区间[a,b]上逐步逼近于根。上述例子中，通过给定初值x=1，经过14次迭代后，得到根为0.80072，精度为0.00036062。
习题 8 8-1．选择题 1．一定量的理想气体，分别经历习题8-1(1)(a) 图所示的abc 过程(图中虚线ac 为等温线)和习题8-1(1)(b) 图所示的def 过程(图中虚线df 为绝热线)，试判断这两过程是吸热还是放热（ ） (A) abc 过程吸热，def 过程放热 (B) abc 过程放热，def 过程吸热 (C) abc 过程def 过程都吸热 (D) abc 过程def 过程都放热 2．如习题8-1(2) 图所示，一定量的理想气体从体积V 1膨胀到体积V 2分别经历的过程是：A-B 等压过程；A-C 等温过程； A-D 绝热过程。其中,吸热最多的过程（ ） (A) A-B (B) A-C (C) A-D (D) 既是A-B ，也是A-C ，两者一样多 3．用公式E =νC V ，m T (式中C V ，m 为定容摩尔热容量，ν为气体的物质的量)计算理想气体内能增量时，此式( ) (A) 只适用于准静态的等容过程 (B) 只适用于一切等容过程 (C) 只适用于一切准静态过程 (D) 适用于一切始末态为平衡态的过程 4．要使高温热源的温度T 1升高ΔT ，或使低温热源的温度T 2降低同样的ΔT 值，这两种方法分别可使卡诺循环的效率升高Δ1和Δ2。两者相比有（ ） (A) Δ1>Δ2 (B) Δ1<Δ2 (C) Δ1= Δ2 (D) 无法确定哪个大 5. 理想气体卡诺循环过程的两条绝热线下的面积大小(如习题8-1(5)图中阴影所示)分别为S 1和S 2，则两者的大小关系是（ ） (A) S 1 > S 2 (B) S 1 = S 2 (C) S 1 < S 2 (D) 无法确定 6. 热力学第一定律表明（ ） (A) 系统对外做的功不可能大于系统从外界吸收的热量 (B) 系统内能的增量等于系统从外界吸收的热量 (C) 不可能存在这样的循环过程，在此循环过程中，外界对系统做的功不等于系统传给外界的热量 (D) 热机的效率不可能等于1 7. 根据热力学第二定律可知（ ） (A) 功可以全部转换为热，但热不能全部转换为功 (B) 热可以从高温物体传到低温物体，但不能从低温物体传到高温物体 (C) 不可逆过程就是不能向相反方向进行的过程 (D) 一切宏观的自发过程都是不可逆的 8.不可逆过程是（ ） (A) 不能反向进行的过程 (B) 系统不能回复到初始状态的过程 (C) 有摩擦存在的过程或者非准静态过程 (D) 外界有变化的过程 习题8-1(1)图 习题8-1(2)图 习题8-1(5)图
西南石油大学油层物理复习资料1.txt我的优点是：我很帅；但是我的缺点是：我帅的不明显。什么是幸福？幸福就是猫吃鱼，狗吃肉，奥特曼打小怪兽！令堂可是令尊表姐？我是胖人，不是粗人。本文由梓悟青柠贡献 doc文档可能在WAP端浏览体验不佳。建议您优先选择TXT，或下载源文件到本机查看。 第一章油层岩石的物理特性 1. 什么是油藏？油藏的沉积特点及其与岩石特性之间的关系是什么？ 2. 沉积岩有几大类？各自有些什么特点？ 3. 油藏物性参数有些什么特点？通常的测定方法是什么？ 4. 什么是粒度组成？ 5. 粒度的分析方法有哪些？其基本原理是什么？ 6. 粒度分析的结果是如何表示的？各自有些什么特点？ 7. 如何计算岩石颗粒的直径，粒度组成，不均匀系数和分选系数？ 8. 岩石中一般有哪些胶结物？它们各自有些什么特点？对油田开发过程会产生什么影响，如何克服或降低其影响程度？ 9. 通常的岩类学分析方法有哪些？ 10.如何评价储层的敏感性（具体化，包括评价地层伤害的程度）？ 11.如何划分胶结类型，其依据是什么？它与岩石物性的关系怎样？ 12.什么是岩石的比面？通常的测试方法有哪些？其原理是什么？ 13.推导岩石的比面与粒度组成之间的关系？ 14.粒度及比面有何用途？ 15.什么是岩石的孔隙度，其一般的变化规律是什么？ 16.按孔隙体积的大小可把孔隙度分为几类？各自有些什么特点及用途？ 17.孔隙度的测定方法有哪些？各自有什么特点？ 18.孔隙度有些什么影响因素，如何影响的？ 19.岩石的压缩系数反映了岩石的什么性质？是如何定义的？ 20.综合弹性系数的意义是什么？其计算式为： C * = C f + C Lφ 式中各物理量的含义是什么？ 21.当油藏中同时含有油，气、水三相时，试推导： C= C f + φ (S o C o + S w C w + S f C f ) 22.试推导分别以岩石体积，岩石骨架体积和岩石孔隙体积为基准的比面之间的关系 S = S s (1 ? ? ) = φ ? S p S―以岩石体积为基准的比面， S p ―以岩石空隙体积为基准的比面， S s ―以岩石骨架体积为基准的比面。 23.什么是岩石的渗透性？什么是渗透率？岩石渗透率的“1 达西”的物理意义是什么？ 24.什么是岩石的绝对渗透率？测定岩石绝对渗透率的限制条件是什么？如何实现这些条件？ 25.达西定律及其适用范围是什么？ 26.试从理论及实验两方面证明渗透率的物理意义。 27.渗透率可分为几大类，其依据是什么？ 28.水测，油测及气测渗透率在哪些方面表现出它们之间的差别？ 29.从分子运动论的观点说明在什么条件下滑脱效应对渗透率无影响，这一结论在理论和实验工作中有什么用途？ 30.影响渗透率的因素有哪些？是如何影响的？ 31.什么是束缚水饱和度，原始含油饱和度及残余油饱和度，在地层中它们以什么方式存在？ 32.流体饱和度是如何定义的？ 33.对低渗岩芯，能用常压下的气测渗透率方法来测其绝对渗透率吗？ 34.测定饱和度的方法有哪些？它们各自有何优劣点？ 35.什么是等效渗流阻力原理？利用等效渗流阻力原理推导出岩石的渗透率，孔隙度及孔道半径之间的关系。 36.推导引入迂回度后，孔隙度，渗透率，比面及孔道半径之间的关系。 37.推导泊稷叶方程。 38.在测定岩石的比面时，分析产生误差的原因。 39.矿场上是怎样用岩石的孔隙度和渗透率指标划分储油气岩层好坏的？ 40.有人说：“岩石的孔隙度越大，其渗透率越大”。这种说法对吗？为什么？ 41.试述实验室测定岩石孔隙度的基本原理和数据处理方法，并画出实验仪器的流程示意图。 42.试述实验室测定岩石渗透率的基本原理和数据处理方法，并画出实验仪器的流程示意图。 43.试述实验室测定岩石比面的基本原理和数据处理方法，并画出实验仪器的流程示意图。
15章习题参考答案 15-3求各图中点P 处磁感应强度的大小和方向。 [解] (a) 因为长直导线对空间任一点产生的磁感应强度为： 对于导线1：01=θ，2 2π θ= ，因此a I B πμ401= 对于导线2：πθθ==21，因此02=B 方向垂直纸面向外。 (b) 因为长直导线对空间任一点产生的磁感应强度为： 对于导线1：01=θ，2 2π θ= ，因此r I a I B πμπμ44001= = ，方向垂直纸面向内。 对于导线2：21π θ=，πθ=2，因此r I a I B πμπμ44002==，方向垂直纸面向内。 半圆形导线在P 点产生的磁场方向也是垂直纸面向内，大小为半径相同、电流相同的 圆形导线在圆心处产生的磁感应强度的一半，即 r I r I B 4221003μμ= = ，方向垂直纸面向内。 所以，r I r I r I r I r I B B B B 4244400000321p μπμμπμπμ+=++=++= (c) P 点到三角形每条边的距离都是 o 301=θ，o 1502=θ 每条边上的电流在P 点产生的磁感应强度的方向都是垂直纸面向内，大小都是 故P 点总的磁感应强度大小为 方向垂直纸面向内。 15-4在半径为R 和r 的两圆周之间，有一总匝数为N 的均匀密绕平面线圈，通有电流I ，方向如图所示。求中心O 处的磁感应强度。 [解] 由题意知，均匀密绕平面线圈等效于通以 I NI 圆盘，设单位长度线圈匝数为n 建立如图坐标，取一半径为x 厚度为dx 的 圆环,其等效电流为: 方向垂直纸面向外. 15-5电流均匀地流过一无限长薄壁半圆筒，设电流I =5.0A ，圆筒半径 R =m 100.12?如图所示。求轴线上一点的磁感应强度。 [解] 把无限长薄壁半圆筒分割成无数细条，每一细条可看作一无限长直导线，取一微元d l 则I R l I πd d = 则l d 在O 点所产生的磁场为 又因，θd d R l = 所以，R I R I B 2002d 2d d πθ μπμ== θcos d d x B B =，θsin d d y B B = 半圆筒对O 点产生的磁场为：
本科毕业设计（论文） 题目光固化石墨烯/聚吡咯复合膜 对电极的制备 学生姓名学号 教学院系化学化工学院 专业年级化学工程与工艺2009级 指导教师职称副教授 单位西南石油大学 辅导教师职称 单位 完成日期2013 年 6 月12 日
Southwest Petroleum University Graduation Thesis Uv-curable graphene/polypyrrole Composite film electrode preparation Grade: 2009 Name: Yang Yin Long Specialty: Chemical Engineering and Technology Instructor: Ke Qiang
摘要 对电极是染料敏化太阳能电池的重要组成部分，改进对电极是提高其能量转换效率及降低成本的有效手段之一,本文利用乳液聚合法制备了石墨烯/聚吡咯复合电极材料，用SEM，FT-IR，EIS等分析手段对复合材料的结构进行表征，最后在环氧丙烯酸酯的包裹下制备了光固化复合膜作为对电极。结果表明石墨烯和聚吡咯均匀的聚合在一起，光固化后的复合材料在EIS测试中具有良好的导电性能，这可能得益于良好分散的石墨烯片层增大了电极材料的比表面积所致。 关键词：石墨烯；聚吡咯；光固化；乳液聚合。
Abstract counter electrode is an important part of the dye sensitized solar cells, the improvement of electrode is to improve the energy conversion efficiency, and one of the effective means to reduce costs, the paper was prepared by emulsion polymerization method of graphene/polypyrrole composite electrode materials, SEM, FT - IR, the EIS analysis means such as the structure of the composites were characterized, the last under the parcel of epoxy acrylate light curing composite membrane was prepared as the counter electrode. Results showed that graphene and polypyrrole together evenly and light cured composite materials with good electric conductivity in the EIS tests, this may benefit from the good dispersion of graphene layers caused by increasing the specific surface area of the electrode materials. Key words: graphene; Polypyrrole; Light curing; Emulsion polymerization. II
第一章油层岩石的物理特性 1.什么是油藏？油藏的沉积特点及其与岩石特性之间的关系是什 么？ 2.沉积岩有几大类？各自有些什么特点？ 3.油藏物性参数有些什么特点？通常的测定方法是什么？ 4.什么是粒度组成？ 5.粒度的分析方法有哪些？其基本原理是什么？ 6.粒度分析的结果是如何表示的？各自有些什么特点？ 7.如何计算岩石颗粒的直径，粒度组成，不均匀系数和分选系数？ 8.岩石中一般有哪些胶结物？它们各自有些什么特点？对油田开发 过程会产生什么影响，如何克服或降低其影响程度？ 9.通常的岩类学分析方法有哪些？ 10.如何评价储层的敏感性（具体化，包括评价地层伤害的程度）？ 11.如何划分胶结类型，其依据是什么？它与岩石物性的关系怎样？ 12.什么是岩石的比面？通常的测试方法有哪些？其原理是什么？ 13.推导岩石的比面与粒度组成之间的关系？ 14.粒度及比面有何用途？ 15.什么是岩石的孔隙度，其一般的变化规律是什么？ 16.按孔隙体积的大小可把孔隙度分为几类？各自有些什么特点及用 途？ 17.孔隙度的测定方法有哪些？各自有什么特点？ 18.孔隙度有些什么影响因素，如何影响的？
19. 岩石的压缩系数反映了岩石的什么性质？是如何定义的？ 20. 综合弹性系数的意义是什么？其计算式为： φL f C C C +=* 式中各物理量的含义是什么？ 21. 当油藏中同时含有油，气、水三相时，试推导： C= )(C S C S C S C f f w w o o f +++φ 22. 试推导分别以岩石体积，岩石骨架体积和岩石孔隙体积为基准的比面之间的关系 S S p s S ?=-=φ?)1( S ―以岩石体积为基准的比面，S p ―以岩石空隙体积为基准的比面，S s ―以岩石骨架体积为基准的比面。 23. 什么是岩石的渗透性？什么是渗透率？岩石渗透率的“1达西”的物理意义是什么？ 24. 什么是岩石的绝对渗透率？测定岩石绝对渗透率的限制条件是什么？如何实现这些条件？ 25. 达西定律及其适用范围是什么？ 26. 试从理论及实验两方面证明渗透率的物理意义。 27. 渗透率可分为几大类，其依据是什么？ 28. 水测，油测及气测渗透率在哪些方面表现出它们之间的差别？ 29. 从分子运动论的观点说明在什么条件下滑脱效应对渗透率无影响，这一结论在理论和实验工作中有什么用途？ 30. 影响渗透率的因素有哪些？是如何影响的？