请在此输入您的回答，每一次专业解答都将打造您的权威形象数据源：（是什么）研究区域描述：（如果你研究的是区域的话，要写出研究区域你要研究的那一方面的发展概况）数据处理方法：你用了什么方法，仔细描绘，比如怎么选取变量，有无修正参数或部分数据啦等等，怎么检验你处理的方法是否恰当啦
第一节论文写作的八个环环相联的步骤严格地说，论文写作并不是从提笔写（或在电脑上打字）开始的。此前的许多步骤都属于论文写作的必要环节，一定程度上比实际动手写重要得多。许多过来的人都有体会，完成一篇较大的论文，准备时间少则数月，多则数年，一旦真正准备好了，动手写作的时间不过数天到数周。第一步，确定论文的选题。从广义上说，选任何本专业范围内的题目都能够写出东西来，只要你有新观点、新发现、新角度、新研究方法、新材料等等。但是这后面的“五新”大大限制了硕士论文的选题。这是由于作者多数是第一次写这么长的学术论文，缺乏经验，也缺乏深厚的知识积累，难以把握；同时，二三万字这个条件也对选题有很大的制约，如果题目过大，无法在这个相对狭小的范围内展开。所以，选题是否得当，对于论文的成功，影响很大，甚至有人说，一个好的选题等于成功了一半。根据许多硕士论文的选题经验，这一级论文的选题可从以下几方面考虑：本专业的研究空白、发生争议的话题（自己的观点感到较为充分）、对比性的话题、从其他专业角度研究本专业的话题（这是一种选题的边际效应）、有新的插入角度的老话题、刚刚冒出来的本专业的新问题。第二步，围绕已经确定的论文选题，回顾相关的理论和研究，或者叫“文献检索”。这一步的工作是较为艰苦的，需要有思想准备。在我国，多数中文学术资料目前没有上网，需要手工查找，因而这个步骤中查找中文资料花费的时间和精力可能很大；拉丁文资料，特别是英文资料由于网络传播的方便条件，相对好查询。但是不少资料即使找到了目录，真正能够阅读到，仍需要作者不懈的努力。这一步是必要的，如果没有这一步，你的论文内容很可能重复了别人已经做过的工作，等于白做；查找的过程，也是启发思路、产生观点火花的过程，不走这一步，等于掐掉了自己新观点、新视角、新材料的来源。这也是为下一步做观点、角度、材料上的准备。第三步，提出你自己关于选题的理论假设，或要研究的具体问题。选题是指准备写的论文的大体方向和范围，真要动手写作，就会遇到两类具体的问题。第一类属于观点方面的：我的具体观点是什么？你可以设想出一个或几个观点，但它们仅仅是一种假设，通过许多证据、材料，通过严密的论证和适当的论证框架结构，证明你的假设是成立的，这才能形成论文的主体。第二类属于实用方面的：我要具体论证什么问题？你可以提出许多原因、各种环境条件的影响，它们是不是与所论证的问题相关，相关到什么程度，这需要通过科学的调查和分析。不论哪一种情况，这涉及论文的中心思想或论证主题，一定要明确，并且贯穿论文的始终。由于硕士论文字数相对长，常见的问题之一，在于作者把握不住全文，写着写着，无形中脱离了自己原来确定的假设或具体问题，说了许多无关中心思想或论证主题的内容。第四步，决定采用哪些研究方法。人文－社会科学的研究方法，大体可以归为两大范畴，思辨研究和实证研究，后者又可分为定性研究、定量研究两种具体的研究方法。人们为探究社会事实或社会现象，而采用不同的研究取向，不同的研究取向又有不同的研究方法，不同的研究假设、收集资料的方式和对结果的判断标准。但是各种研究方法在现在的论文写作中，已经越来越多地呈现相容和内在的连接。一般地说，根据自己的选题和讨论的具体问题，可以以一种研究方法为主，辅以其他的方法。例如研究“人”作为大众媒体信息的接受者其接受信息时的状况，这种研究取向就决定了研究本身要以定量分析为主，但同时也需要一些历史的、文化的、政治经济学的思辨研究。在文科硕士论文中，作者直接为论文进行的定量分析，规模一般较小，适应的范围也是有限的，较多地采用别人而不是自己直接的调查结果。这是由于论文的规模较孝给予作者的研究经费有限、作者个人进行社会调查的能力有限等原因造成的。以逻辑分析为主的论文，适当采用一些定量分析的数据，有时会给文章增添一些分量。但是，一定要根据实际需要，而不要为了显示研究方法的多样而有意去做。例如一个很宏观的话题本来适于思辨研究，硬要加进一项微观的量化调查结果证明什么，反而会弄巧成拙。这是现在写好论文要把握的一个具体问题。第五步，设计论文的框架结构。一般文章的写作也需要有这一步，但对硕士论文来说，更为必要，其要求也更细一些。一般情况下，一篇硕士论文要有绪章、入题的第一章、主体章节，以及结束语。章节的设置在写前要有个大体的布局逻辑，使之结构合理；章和章之间有一种逻辑联系，防止盲目写下去，淹没主题，不知所云。这一步很少有一次完成的，往往会根据收集材料的情况、调查访问中遇到的新情况，经常变动。但是就像建筑师在盖房子前必须有图纸一样，到了写硕士论文这个层次上，大体的文章框架不能仅仅存于脑子中，一般要形成文字，相对细致一些，具体到“节”更好（但“节”的层次开始时不要固定化），便于写作时心中有数。到了设计论文框架这一步，因为有了文字化的章节设计，除了请导师指导外，这是在正式动笔写前较广泛地征求其他专家意见的一个好机会。框架还不是厚厚的论文，看时花费的时间不多，又可以大体看出文章的价值或存在的问题。这时修改论文结构比写完后修改要轻松、容易得多，时间也较为宽余，不要错过这个机会。第六步，对已经取得的文献资料、调查材料和各种论据进行分析、归类，分别充实到各章节中，再进行解释、论证。这实际就是论文写作本身，所以这样描述，意在让作者理解论文写作的过程。各种材料和论据，不是天生就可以证明论点或说明具体问题的，需要通过作者对材料的组织和论证，才能使其变得富有生命力，极其自然、有力地为自己所论的题目服务。在这一步，需要温习一下学过的逻辑学或社会调查统计的知识，用正确的逻辑思维和严谨的数据组织方式，紧紧围绕已经确定的理论假设或具体问题，调动自己所学的各种知识，通过正论（这是主要的）、反论、设论、驳论、喻论等等手法，论证观点或问题，得出结论，完成论文。论证中肯定会出现种种材料使用或缺乏的问题、逻辑推理的问题、论据与论题不相配的问题等等，需要停下来再找材料和访问专家，充实或削减原来论文框架中的内容，必要时对框架结构进行局部调整。这种情形是正常的、经常发生的。在时间的安排上，对此要作出计划。如果时间安排不当，有时论文功亏一篑的原因就在于写作时间安排过紧，来不及调整论文结构，这很不值得。第七步，必要时重新估量选题，修正论证对象的范围。这是与第六步同时出现的另一种情形，即通过较为广泛地征求意见和本人的思考，感到原来的选题对自己不适当，或难以完成，那么就要及时调整整个论文写作的计划，改变选题。这种情形也是正常的，关键在于不要长时期犹豫不决，必须较快地作出决定，以便有时间重起炉灶。由于前面已经对本专业的学术研究有过较多的思考和文献检索，即使改变选题，重新做起，花费的时间不会很多，对此过多的担心是不必要的。选题不当、难以完成的另一种情况不在于选题本身，而在于选题论证的范围过大。解决这个问题并不难，把论证对象的范围缩小就是了。这里最大的障碍在于作者舍不得“割爱”，花费了许多功夫准备论文，一旦许多材料用不上，难以割舍。这种情况当然会涉及到重新设计论文框架结构的问题。不过，将较大的论证对象的范围缩小，总比相反的情形要容易得多。硕士论文写作中，论证对象范围过小的情况很少见，因为二三万字的论文，本来大多适宜开口较小的选题。第八步，对论文从技术上进行规范化的检查和调整。章节设计的技术问题（含目录）、文中的引证标示、注释及编号、文后的参考文献编排，以及不属于论文本身的内容提要（包括英文提要）、关键词等等，都要按照规范化的要求进行检查和调整。这些虽然属于技术性问题，但也反映出作者的治学态度。特别是引证，凡是使用了别人观点的地方，都必须注明材料来源，不能含糊不清，更不能将别人的研究成果变成自己的。标明的材料来源也要十分清楚，论著名称、作者或编者、出版社或发表的刊物名称、出版或发表时间等等，一应具全。有时，一篇较有水平的论文，答辩时提出批评的主要问题是引证的不规范，由于这个原因使论文的评价低一档，这很不值得。第二节论文的选题：一个适当的选题等于论文成功了一半由于十年“文化大革命”这场浩劫，我国的人文－社会科学研究遭到了毁灭性的破坏，林彪、“四人帮”在意识形态领域的种种谬论充斥人文－社会科学的各个学科。1978年我国恢复研究生招生时，人文－社会科学研究几乎一切从头开始，因而研究生的论文选题还没有显现为一个问题，因为写什么题目，即使不是空白，也是需要重新认识和论证的。经过改革开放后二十多年几代人文－社会科学研究者的努力，我国的人文－社会科学研究呈现一派繁荣的景象，研究生的招收数量大幅度提高。年年招生，年年写新的论文，而各学科的研究领域是相对稳定的，空白越来越少，新问题的研究又有一定难度，于是现在文科研究生的论文选题，形成一种独特的竞争局面，如何选择好论证的题目，凸现为一个新问题。有鉴于此，需要重点谈谈硕士论文的选题问题。
数据分析可以分成两部分，一部分是对分析过程及分析结果的描述，另一部分是结合专业知识对结果进一步分析，为什么会出现这样的结果。
数据分析是指用适当的统计分析方法对收集来的大量数据进行分析，将它们加以汇总和理解并消化，以求最大化地开发数据的功能，发挥数据的作用。数据分析是为了提取有用信息和形成结论而对数据加以详细研究和概括总结的过程。
论文数据分析分析目的：
数据分析的目的是把隐藏在一大批看来杂乱无章的数据中的信息集中和提炼出来，从而找出所研究对象的内在规律。在实际应用中，数据分析可帮助人们做出判断，以便采取适当行动。数据分析是有组织有目的地收集数据、分析数据，使之成为信息的过程。
这一过程是质量管理体系的支持过程。在产品的整个寿命周期，包括从市场调研到售后服务和最终处置的各个过程都需要适当运用数据分析过程，以提升有效性。
论文的数据分析怎么写如下：
首先，针对实证性论文而言，在开始撰写论文之前，必须要提前确定好数据研究方法。而数据研究方法的确定与选择需要根据大家毕业论文的研究课题来确定。
另外，大家也可以跟自己的的论文指导老师多多交流，尽可能多的了解更多关于研究方法的知识，以供自己选择。除此之外，大家还需要大量查找文献资料，见多识广有大量输入之后才能有所输出，本环节需要大家跟导师沟通商议后决定。
接下来一个比较重要的步骤是搜集和整理实验数据。在这一部分，很多同学朋友都会遇到各种各样的问题，比如，不知道去哪里找数据，找到的数据可靠性无法保障，需要的数据总是无法搜集全面等等各种问题。
那么在这里需要跟大家强调一下，推荐大家使用国家统计局、中国统计年鉴、国泰安、万方等等这些比较权威的网站去搜集数据资料。
在此需要注意的是，国泰安和万方等这些网站是需要收费的，上去看了一下，价格不是很亲民。
给大家分享一下，如果有些数据在国家官方网站确实找不到或者毕业论文所需的最新数据还没及时发布，推荐大家可以上某宝，因为某宝上电子版数据往往都很全面，而且价格大都可以接受。
在此提醒大家搜集到数据之后，一定要按照自己的习惯整理保存好，避免后期使用数据时出现差错。
少年锦时，应是山间不羁的一首诗。 骑最俊的马，喝最烈的酒，看最美的夕阳，爱着心头的姑娘
研究生论文可以参考的数据库可以有：国内的知网、万方、维普、百度学术等，国外的Google Scholar、SCI-Hub、Web of Science等，这些数据库是每个专业写论文必备的数据库，当然，你还需要与你专业相关的专业数据库，一般在学校图书馆里的资源利用中会有。
1、收录情况不同：从收录情况来看，维普收录最久，万方收录时间最短，但是现刊万方收录最好，维普最近两年的很多期刊都没有了，特别是医学类的。
2、清晰度不同：文档的清晰度：万方最高，维普最差。文本清晰度关系到用户的满意度。主要原因是CNKI的使用范围、收入论文数量、质量、用户人数较之于万方数据库要优秀，基于CNKI指标更能反映出论文在学术市场的接受度。
3、使用方便性不同：知网和万方功能一样，维普使用不方便。万方的期刊浏览功能做的很好。《中国期刊全文数据库（CJFD）》截至2007年3月共收录了国内8127种期刊，全文文献总量已达2272万多篇，收录起始时间一般为1994年，个别追溯到1979年或创刊年，按学科分126个专题，内容每日累增。产品形式有网络形式、光盘形式、《中国期刊专题全文数据库光盘版》。
相同点都是知名文档搜集网站。提供CNKI源数据库、外文类、工业类、农业类、医药卫生类、经济类和教育类多种数据库。其中综合性数据库为中国期刊全文数据库、中国博士学位论文数据库、中国优秀硕士学位论文全文数据库、中国重要报纸全文数据库和中国重要会议文论全文数据库。
参考资料来源：百度百科-万方数据库
参考资料来源：百度百科-中国知网
参考资料来源：百度百科-维普网
1、知网 知网是国内最权威、应用最广泛的数据库,想必大家都很熟悉。在国内的晋升考核和毕业考核中,知网查重的要求最为普遍。知网包括了国内绝大多数的期刊和文献,而知网是国内学术论文出版的最终要求,比国内其他数据库使用和认可也更广泛。 Sci是世界上学术价值最高、认可度最高的数据库,可以说是衡量国内外科研水平的最高标尺。国内高水平职称评审有sci 论文发表的相关要求,除此之外,研究生和博士生也有sci 论文出版的要求。 Ei是工程技术的索引,是国际范围内仅次于sci的数据库。具有很强的专业针对性
随着自由贸易的发展,我国所遭受的倾销越来越严重,运用反倾销 措施 来维护公平竞争的市场环境、保护国内企业的合法利益、保证产业安全已刻不容缓。下面是我为大家推荐的化工论文，供大家参考。
化工论文 范文 一：能源化学工程专业无机化学教学改革
能源化学工程专业[1]是利用化学、化工的理论与技术来解决能量的转换、储存及传输等问题，通过生产清洁、高效的新能源服务于人类生活的一门学科。无机化学是本专业所开设的第一门专业基础课，其教学质量直接影响到培养的应用创新型人才的质量。而目前无机化学的教学中面临着很多问题，如大一新生刚从高中迈入大学，面临如此信息量大的课程感到迷茫;教师面对课时量日趋减少的趋势，而传递的信息量大的困扰，不知如何把握日常教学;另外，加上教师科研压力等方面的因素，使得其未能全身心地投入教学中。因此，无机化学教学的改革与探讨在本专业教学过程、人才培养模式中的地位尤为重要。例如:
(1)武汉工程大学化工与制药学院从优化课程内容入手，对无机化学的 教学 方法 进行了改革[2];
(2)钦州学院化学化工学院从无机化学的重要地位出发，结合无机化学的教学目的，对无机化学多媒体课件进行了构建和探讨[3]。菏泽学院是一个应用型的地方性教学型本科院校，于2012年成功申请了与国家战略性新兴产业密切相关的能源化工专业。我系主要从教学目标、教学内容、现代化的教学手段等方面对无机化学的教学进行了改革与探索。
1明确合理的教学目标
根据能源化学工程专业的培养目标及培养模式，结合无机化学课程特点，菏泽学院化学化工系于2012年制定了能源化工无机化学教学目标。通过该课程的理论基础及实验实践的学习，能够使学生掌握无机化学基本知识和技能，为培养成高素质劳动者和化工专业技能人才做好准备;同时，也为今后学习专业知识和职业技能打下坚实的基础。此目标主要分为以下几个方面的目标。
知识目标
主要分为了解、理解、掌握三个层次方面目标。通过该课程的教学，应使学生了解:气体的扩散定律，气体分子的速率分布和能量分布;反应速率的概念及反应速率理论;强电解质解离、离子氛、活度系数的概念;微观粒子运动的特殊性;路易斯结构式，等电子体原理，分子轨道理论;化学电源与电解;卤素单质的物理性质，金属卤化物、拟卤素和拟卤化物、互卤化物和多卤化物;硫和硫化物、单质硫、硫化氢和氢硫酸的物理性质;硅的单质、硅烷、硅的卤化物、硅的含氧化合物。通过该课程的教学，应使学生理解和掌握:气体的状态方程及混合气体的分压定律;热力学第一定律，化学反应的热效应、热化学方程式、盖斯定律、生成热的概念及应用，化学反应进行方向的判断方法;浓度对反应速率的影响;缓冲溶液的原理及应用;沉淀溶解平衡及移动;核外电子运动的描述，核外电子排布和元素周期律及基本性质的周期性;价键理论，价层电子互斥理论及杂化轨道理论;基本概念:原电池、电极电势和电动势及能斯特方程;卤素单质的化学性质，卤化氢和氢卤酸的化学性质;氧、氧化物、臭氧、过氧化氢的物理化学性质，硫的含氧化合物的化学性质。掌握氮的氢化物、氮的含氧化合物的化学性质。
专业能力与素质目标
能力目标方面主要是培养学生谦虚的品格、勤奋好学的习惯以及知识迁移的能力;培养学生勤于动手创作、做事严谨的良好作风;培养学生学会运用唯物主义辨证的思维分析问题及解决问题的能力;培养学生工程质量意识和规范意识以及严谨、认真的工作态度。专业能力目标方面使学生能够掌握重要元素及其化合物的主要性质、结构、存在、制法、用途等基本知识;培养学生独立进行化学计算和利用参考资料等方面的能力;具有通过对实验数据的分析，绘制出特性曲线，能够写出规范实验 报告 并加以 总结 概括的能力。素质目标方面主要是培养学生具备良好的职业道德;培养学生勤苦奋斗、勇于创新、敬业乐业的工作作风。
2丰富合理的教学内容
科研成果与课堂教学相结合，保持教学内容的前沿性
科研成果与课堂教学相结合包含两部分内容:一是在教学过程，教师能将自己的科研成果带入教学内容之中。这就要求教师教学的同时展开科研，而科研课题也要紧紧围绕教学内容展开，这样会更能了解学科的前沿动态并能深入把握，有利于增强教学的深度、广度，有效地提高教学质量[4]。另外教师将科研成果带入课堂分析中，将科研成果与教学有机地结合起来，将最新知识与信息传递给学生，科研推动教学，教学促进科研。二是在教学过程中结合学科发展情况，充分利用别人的研究成果，及时补充教学内容，进行教材建设。另外，在教学实践中可采用“案例教学”，对具体科研案例进行讨论、分析，比较各种方案的优缺点及产生原因，选择合理方案。在项目设计过程中，通过教师的引导作用，学生可以自主查阅资料并开展项目的研究性学习。
建设开放的无机化学实验教学环境，理论与实验相结合
充分利用我系基础实验室和化学工程实验中心的仪器设备和师资力量，结合我系化学能源工程专业及无机化学教学内容的特点，试图探索出一套完善的开放式无机化学实验教学模式，注重实验与课堂教学相结合、开展系内实验技能竞赛及无机化学创新实验设计竞赛等项目，激励学生的学习积极性及培养今后创新实践的能力。开展大学生创新研究计划，引导学生在大三下学期进入教师的科研室进行锻炼，参与课题的研究，培养学生的创新意识和实践能力;鼓励大二学生参加无机化学实验技能竞赛，鼓励学生进行科技创新;另外聘请国内外无机化学研究领域的专家学者来我系作学术报告，增加学生的科研兴趣及全面了解无机化学的前沿动态，为今后的科研之路做好准备。
3多媒体与板书相结合的现代化教学手段
针对目前无机化学课时缩减而传递信息量大的情况，传统的板书教学手段已不能满足时代的需要，因此多媒体技术已广泛使用在课堂教学中。这样一方面将节省下的板书的时间能够用于重点难点的讲解，另一方面多媒体中引入一些无机化学演示实验、实物图像，将枯燥的理论教学表现的更加生动直观，提高了学生的学习积极性。然而仅利用多媒体也有一定的缺陷，如对一些公式的推导，仅利用多媒体会受到一定的限制，因此多媒体跟板书结合会更加有利于公式的推导。另外，还会避免仅利用多媒体的教学进度过快，学生不能融会贯通的缺点。总之，鼓励学生 课前预习 ，采用板书与多媒体技术相结合既能考虑教师的教学进度与学生的掌握程度，又能兼顾教学的广度与深度的问题，取得了较好的教学效果。
4结束语
无机化学是能源化学工程专业学生迈入大学的第一门专业基础课，其教学效果直接影响着学生学习本专业的积极性及掌握本专业基础知识的扎实程度。本系以上结合能源化学工程专业特点对无机化学的教学目标、教学内容及教学手段的初探具有一定的意义。今后会继续探索无机化学其他方面的改革。
化工论文范文二：油藏化学工程研究发展趋势
推动我国油藏化学工程研究与我国社会进步有着密不可分的联系。为了赶上发达国家对油藏化学工程研究的脚步，我国必须大幅度提升在这一方面的开发技术，更好地促进化学工程研究大步向前发展。
1油藏化学工程研究的发展背景
人类面临的最大危机之一就是能源问题，世界各国都在担忧石油问题。迄今为止，人类只开采了大约总储藏量1/3的原油，因此，油藏开发及提高效率是每一个科技工作人员的头等任务。半世纪以前，世界对石油的总需求量日益增长，工人们利用油藏工程的原理提高采收率来满足市场需求，同时也促进了油藏工程原理的发展。作为石油工程的重要组成部分，油藏工程主要负责各类研究，在掌握动态规律与原理的同时，也辅助了钻井与采油工程的开展。
2三次采油技术
自改革开放以来，世界各国石油界的精英们一直努力提高石油的采收率。一次和二次采油主要是靠自身压力和注气注水等方法，三次采油是采用之前的任何工业技术[2]。因而提高油藏采收率并没有局限在某一阶段或手段，它主要是靠原来油藏中没有的物料开采。它的定义与分类是不矛盾的。油藏化学工程是在三次采油的背景下发展起来的，它和化学工程学科共同发展。随着现代科技的迅猛发展，人们不断引进新技术，取得新成就。这一阶段也让人们认识到发展的多样性，开始探究多方面技术，涉及各种学科，主要有胶体与界面科学、化学工程学、化学反应动力学、渗流力学、热力学、计算数学等多种高等学科。
3化学复合驱技术
我国油田多数是陆相沉积，分布相当不均匀，原油中的蜡含量和芳烃含量比例较大，且黏度大，导致水驱采收率只在33%左右。三次采油的研究技术表明，化学复合驱能够有效提高采收率，它是在单一化学剂驱的基础上组合两种不同的化学剂，形成多种复合体系。通过实验证明，复合驱的相互作用比单一化学驱剂效果显著的多。随着各方面技术的发展和完善，复合驱逐渐成为我国提高原油采收率的主导技术。复合驱配方体系主要是由高浓度小段塞和低浓度大段塞2种体系组成。高浓度小段塞是利用表面活性剂和助剂，使油水形成中相微乳液体系，增强原油的乳化。典型的代表有胶束.聚合物驱体系，它的表面活性剂浓度在～，段塞小于，若形成微乳液，效率更大，能达到80%以上。低浓度大段塞是后期才引进的策略，它的驱油原理主要是毛管准数理论，利用碱和表面活性剂降低油水界面张力。这种体系应用相对广泛，高酸值和低酸值都适用。近年来，随着研究力度加强，新型产品不断出现，如梳形聚合物KYPAM，星形聚合物STARPAM，疏水缔合聚合物。这些新型耐温抗盐聚合物，有利于节约淡水资源，保护环境。也扩展了油藏水的矿化度和文档范围。
4油田堵水调剖技术
开发油田主要采用水驱开发在在这一过程中，因储存分布不均，导致注水过程中出现沿高渗透带窜流，水波效果差，油井含水快速上升，尤其当进入高含水阶段，会出现水短路的现象，加深开采工作难度。为改变这一现状，专家们提出采用“堵水调剖”这一方法。堵水调剖具有颇多优势，操作简洁、规模较小、周期短、效果显著，能有效提高注水开发效果。油田堵水调剖技术历经磨难，从单井油井堵水油井堵水到单井水井调剖，目前主要发展到调整深部调驱。直到2006年底，才开始着手整体堵水调剖示范工程，在采油研究院的带领下，全面开展现工作，有条理的分析堵水调剖工艺技术，给予独特的评价以及实地示范。为改善注水开发的现状，应做如下调整目标：将单井措施向区块整体转变;将近井剖面转向深部液流;阶段上实施一体化转变;评价上从单井向整体转变;应用上改用多种复杂油藏，不再局限在常规水驱油藏据调查，仍有多个区块可以进行整体调堵，由此看来，堵水调剖技术发展趋势将奋力往前。
5评价与改进
综上所述，虽然油藏采收率明显提高，技术也不断突破，但仍然要看清形势。在取得成果的同时，也要擅于总结 经验 ，找出不足，精心解析。例如耐温抗盐聚合物产品的溶解性和长期热稳定性都还不是很乐观，在现场实施过程中，不能有效地达到施工要求，高效率的完成任务。同样地，化学驱技术需要改进解决的问题也是各方面的，需要研究者在过程中分层次去进行。只有抱着永不止步的态度去钻研，去创新，去探索，才能攻克这些技术上遇到的“疑难杂症”，才能进一步将化学驱油技术往特色道路上发展，不断为油藏化学工程研究的发展做贡献。
6结束语
为推动我国油藏化学工程持续发展，还需加强工作。不停探索实验技术，顺应环境变化。掌握化学驱技术，在实际工作中解决问题。还要继续研究物理化学模型，对敏感参数进行验证。油藏化学工程研究的全方位发展，有利于解决能源紧缺问题，有利于稳定我国石油市场，有利于世界和平。
解国军1，2金之钧1
（1.中国石化石油勘探开发研究院，北京100083；2.中国石油大学（北京）资源与信息学院，北京102249）
摘要 为了对东营凹陷营11 透镜状砂岩油藏的成藏机理进行深入研究，以掌握该类油藏成藏的主要影响因素，本文利用可压缩多孔介质油水两相渗流的基本原理，对其成藏过程进行了二维数值模拟。模拟过程中考虑了砂体区地层沉积（剥蚀）、地层厚度变化、岩石孔渗条件变化、流体物性的变化、毛管压力、相对渗透率和油气生成等一系列过程和参数。模拟再现了石油在砂体中聚集的过程，模拟的砂体的最终含油饱和度及分布与现实情况基本吻合。通过模拟和分析表明，围岩和砂体的毛管压力差异才是驱动石油在类似营11这样的透镜状砂岩油藏中聚集的根本动力，而这一驱动力是由于围岩和砂体物性上的差异以及油气的生成两种因素结合在一起而形成的。
关键词 透镜状砂岩油藏 成藏过程 数值模拟 两相流动 毛管压力 东营凹陷
Two Dimensional Numerical Simulation of Oil-trapping Process of Ying11 Lentoid Sand Reservoir of Dongying Depression
XIE Guo-jun1，2，JIN Zhi-jun1
（ ＆ Production Research Institute，SINOPEC，Beijing100083； and Communication Collage，China University of Petroleum，Beijing102249）
Abstract The oil-trapping process of Ying11 lentoid sand reservoir of Dongying Depression is numerically simulated on two dimension condition based on the theory of two-phase fluid flow in compactable porous media in order to thoroughly study the oil accumulation mechanism and master the dominating influential factors of this kind of various related processes and parameters considered in the simulating process are sedimentation/denudation，the thickness of strata，the porosity and permeability of rock，the physical properties of fluids，capillary pressure，relative permeability，and oil oil accumulation process in the reservoir reappears in the simulation，and the oil saturation and distribution accord with the real is indicated that the fundamental driven force for oil accumulation in lentoid sand reservoir as Ying11 is the difference of the capillary pressures built up between the source rock and reservoir，and the driven force forms from the combination of the difference of the physical properties between source rock and reservoir and the oil generation.
Key words lentoid sand reservoir oil-trapping process numerical simulation two phase fluid flowcapillary pressure Dongying Depression
原生透镜状砂体油藏是一类典型的砂岩岩性油藏，这类油藏一般是由浊积岩砂体被低渗透性泥页岩包围所形成的，砂体内油气来源于周围的源岩，东营凹陷的营11砂体油藏就是这类油藏的典型代表。由于完全被泥岩所包围，因此对于这种油藏形成的油水交替机理在人们看来具有不同于构造或地层油藏的特殊性。关于这类砂体油藏的成藏机理和影响因素，前人有过多种有益的实验研究和理论探索。陈章明等［1］、李丕龙等［2］通过成藏物理模拟试验对原生岩性砂体的成藏过程和影响因素进行了分析。王宁等在岩性油藏成藏过程中考虑了成藏的动力和阻力两种因素［3］；庞雄奇等则从“成藏门限”的角度对砂岩透镜体的成藏控制条件进行了分析［4］。李丕龙等提出了“相”、“势”控油理论，对包括透镜状砂体油藏在内的隐蔽油藏的形成机制进行了分析［5］。隋风贵对浊积砂体油气成藏的主控因素进行了定量分析［6］。
然而无论是上述的实验研究还是理论分析，基本上都是从定性或半定量的角度对该类油藏的成藏过程机理进行的讨论，或只是对这类油藏的含油性相关影响因素进行了分析，而没有涉及其成藏机理，因此都无法更详细地了解原生岩性油藏成藏的整个过程及控制机理。由于透镜体油藏的成藏过程是与其围岩紧密相关的，因此，理解砂体的成藏过程必须将砂岩体的演化过程与围岩的演化过程结合起来统一考虑。本文则是从演化的观点，利用定量的方法，综合考虑地层沉降、温压变化、砂体和围岩物性变化、孔隙流体物性变化、石油生成等各种相关过程，模拟处于围岩包围中的营11透镜体油藏成藏的整个过程，并分析其成藏机理和含油性的主要控制因素。通过对该油藏成藏过程的二维数值模拟，可以更深刻地了解这类油藏成藏过程中的油水交替过程及其力学机制，为探讨这类油藏的成藏机理及其影响因素提供了很好的例证。
1 模拟模型的建立
由于营11透镜状砂岩油藏在成藏过程中涉及由于压实作用而导致围岩和砂体的变形以及油水两相流体在其中的流动过程，并且岩石的变形和流体流动是相互影响的，因此这是一个可变形多孔介质两相流动的流固耦合问题。
与二次运移相比，油气从低渗透源岩中的排出（初次运移）一直是比较难以理解的现象。从油气自源岩中排出的相态来看，现在普遍被接受的观点是大多数油气是通过独立相态排出的［7］，而油气排出的主要动力则来源于压实及生烃等作用产生的过剩地层压力［7～9］。描述流体在多孔介质中低速流动的通用方法是依据达西定律给出的，虽然对于在低渗透性泥页岩地层中达西定律是否适用还存在疑问，但其作为一种描述孔隙流体流动速度和压力关系的有效手段还是被广泛应用于各种排烃模拟中［10～13］。为了模拟石油从源岩中排出并进入被其所包围的砂岩中聚集这一过程，本次模拟也采用了基于达西定律的油水两相渗流模型。模型中油相和水相的压力差即为毛管压力。
由均匀介质弹性力学的广义胡克定律可以推出其应变和应力之间的关系。但对于地质过程的模拟，地层压实作用不同于弹性力学所描述的微小变形过程，从长时间看是一种非弹性的大变形过程，而对于这一过程的地质描述一般采用一种近似的简化关系，即将这种变形转化为岩石孔隙度与其所受到的垂向有效应力之间的指数关系［13～15］。根据Terzaghi方程，垂向有效应力可用岩石总负载与孔隙流体压力之差来表示［10，16］。
生油泥岩可视为由干酪根、无机杂基和孔隙3个部分组成，其中干酪根与无机杂基构成生油岩的骨架。为了处理问题简单，可将干酪根划分为有效干酪根（具有生油潜力，可全部转化为石油）和无效干酪根（不具有生油潜力）。因此，可将生油岩重新划分为以下3个部分，即有效干酪根、不可压缩骨架（包括无效干酪根和无机杂基）和孔隙。模型假设有效干酪根降解将产生同质量的烃并使泥岩骨架厚度减小。而岩石的厚度变化可根据不可压缩骨架体积不变的原理得到。对于砂岩储层可不考虑有效干酪根降解所导致的骨架厚度的变化。生油岩中烃类是其中包含的干酪根热降解的结果，而干酪根的热降解采用化学反应动力学中的一级反应定律来近似描述［17］。根据一级反应定律，干酪根的转化率与剩余的干酪根量成正比，可表示成多个平行的一级反应。而反应常数是由反应活化能、频率因子和反映温度决定的。设同质量的有效干酪根降解可产生相同质量的石油，因此石油生成的速率也就是干酪根的降解速率。
2 相关参数变化
水和油的密度是温度和压力的函数，可采用指数型状态方程来描述［13］。水和油的黏度是影响水和油渗流的参数，水的黏度一般采用与温度相关的函数［13，18］，而本次模拟油的黏度采用了考虑了油的重度和温度的Beggs ＆ Robinson公式［19］。
沉积岩的渗透率对地层流体的流动和异常压力的形成都起着至关重要的作用，一般受沉积岩类型和埋藏深度等因素的影响，其大小有时存在多个数量级上的变化。对于碎屑岩地层，一般情况下渗透率的变化可表示为孔隙度的函数，如Kozeny-Carman方程［10，18］。在本次模拟中采用渗透率与孔隙度为幂函数关系的公式［13，20］。
在包含两相或两相以上非混相流体的渗流系统中需要考虑岩石的毛细管压力特征。由于模拟中处理的基本上是石油生排及聚集的过程，因此只需考虑岩石的驱替毛管压力曲线特征。本次模拟研究采用驱替毛管压力与含水饱和度呈幂律关系的公式［21］：
油气成藏理论与勘探开发技术
式中：Pcb为毛管突破压力；γ为孔隙大小分布指数；Sw为含水饱和度。对应于突破压力的毛管半径可用其与孔隙度和渗透率的经验关系来表示［22］。由Laplace方程可知毛管压力是界面张力、润湿角和毛细管半径的函数。水烃体系界面张力可一般表达为体系温度和油水密度的函数［19］。另外，本次模拟假设岩石完全水湿，可得润湿相接触角为0。因此，将可求得岩石毛管突破压力Pcb。如果要求得驱替毛管压力曲线，还需要确定孔隙大小分布指数。对东营凹陷的28块砂岩样的压汞曲线的拟合分析表明，孔隙大小分布指数基本上是与岩石的孔隙度和绝对渗透率等物性参数无关的参数，本次模拟取其均值。本次成藏模拟对于泥岩也采用相同的突破毛管压力公式和孔隙大小分布指数值。
油和水的相对渗透率采用Brooks-Corey经验关系式表示［13，21，22］，其中油和水的相对渗透率与含水饱和度和孔隙大小分布指数有关。
3 营11砂岩油藏成藏过程模拟
营11砂岩油藏概述
营11砂岩油藏位于东营凹陷的东辛油田西南部，西邻郝家油田，南靠现河庄油田。构造上处于东营凹陷中央隆起带西部，东辛、郝家、现河庄构造断裂带之间的洼陷中央。本次模拟的是营11砂体沙河街组三段中下油藏，探明石油地质储量1248×104t，是东营凹陷迄今为止发现的最大的独立砂体油藏。营11沙河街组三段中下砂体的构造图及模拟剖面线位置见图1。
图1 营11沙河街组三段中下砂体顶面构造图及模拟剖面线位置
营11砂岩油藏模拟的前期准备
模拟的前期准备工作由剖面网格化、原始沉积剖面恢复、上覆地层沉积过程反演和模拟演化过程参数确定等几部分组成。
剖面网格化
选取的剖面长度以营75井为分界点，向砂体上倾方向延伸5600m，向砂体下倾方向延伸2400m，剖面总长度为8000m。剖面体垂直方向深度从2700m（大致为沙河街组三段上亚段的底界面）至3600m（大致为沙河街组四段上亚段底界面）。从沙河街组三段中亚段向沙河街组三段上亚段，砂岩沉积逐渐占据主导地位，由于砂岩较好的导流性，不易形成显著的异常压力，因此在剖面体顶部位置的压力边界条件以常压来考虑。由沙河街组四段上亚段地层向下膏泥岩居主导地位，因此可以沙河街组四段上亚段地层底界为剖面体的封闭边界。由此可见剖面体长8000m，高900m。在网格划分时既要考虑精度，又要考虑计算工作量的大小，因此，在砂体所对应的长度和高度方向进行网格细化，而在其他地方，尽量将网格粗化以减小计算工作量。
原始沉积剖面恢复
由于剖面显示的是现今的沉积厚度和孔隙度特征，要进行砂体成藏过程的正演模拟，需将剖面恢复到模拟零时刻的状态。本次模拟的零时刻设定为沙河街组三段上亚段沉积期末，因此，需将模拟剖面从顶部的2700m恢复到0 时的剖面状态。恢复是按地层压缩时骨架体积不变的原则进行的。地层孔隙度采用随深度按指数递减规律变化的公式，其中相关参数是根据东营凹陷实际探井的地层数据回归得到的。
上覆地层沉积过程反演
由于成藏过程为一正演过程，因此需知道模拟剖面上覆地层在不同沉积期的沉积速率以及地层的砂泥岩含量。为此，首先要了解沉积地层现今的厚度及砂泥岩含量。表1给出了营11砂体区域有代表性井的地层厚度和地层砂质含量以及地层平均沉积速率。其中的地层砂质含量由自然电位或自然伽马测井数据计算得出；地层沉积速率是指沉积物处于沉积表面时的沉积速率，根据地层的砂泥岩含量、地层厚度和深度以及沉积持续时间给出。而东营期末的沉积间断按剥蚀200m的东营组计算，并依据沉积间断的时间得到平均剥蚀速率。
表1 营11砂体上覆地层模拟参数
模拟演化过程参数确定
营11砂体区的古地温梯度采用东营凹陷的古地温梯度，距今时间为43Ma，38Ma，36Ma，，，，2Ma和0时的古地温分别是℃/100m，℃/100m，℃/100m，℃/100m，℃/100m，4℃/100m，℃/100m和℃/100m［23］。
与砂岩岩石压缩有关的参数值由东营凹陷砂岩孔隙度与深度及有效应力的关系回归得到，而与泥岩压缩相关的参数值来自Mudford等［24］。砂岩渗透率与孔隙度关系式中的参数值来自东营凹陷的数据回归，而泥岩参数值来自Luo 和 Vasseur［13］。
岩石的生烃潜力可定义为生油岩有效干酪根（可转化为烃类）占岩石骨架总量的质量比，而原始生烃潜力是指烃源岩在演化的初始时刻的生烃潜力。一般将在岩石热解分析中的S2值视为岩石的生烃潜力值，因此若想得到网格体岩石的生烃潜力值，需要本区大量的有机岩热解分析资料，而现实的情况是这种分析资料在本区非常有限，无法满足网格体的生烃潜力值的数值化。因此，本次模拟网格体的生烃潜力利用营11砂体区的测井数据进行计算。采用Passey等［25］提出的基于孔隙度和电阻率测井数据的ΔLgR方法，经过改进可以对烃源岩在演化初期的原始生烃潜力进行预测。进行网格体原始生烃潜力赋值应用了钻遇营11砂体和其附近的营76井、营101井、营102井、新营69井、营75井、营70井、营67井、营68井、营78井等的测井数据。由于上述井均未钻遇沙河街组四段上亚段地层，因此，模拟剖面沙河街组四段上亚段地层的原始生烃潜力采用河88和郝科1的计算值。
考虑到东营凹陷沙河街组四段上亚段、沙河街组三段下亚段以及沙河街组三段中亚段的烃源岩以I型干酪根为主，在生油模拟中烃源岩的干酪根依反应活化能划分的各组分初始含量和频率因子等参数采用Schenk等［26］提供的I型干酪根数据。
模拟过程及结果分析
营11砂体的成藏模拟从距今开始，即模拟的0时间点，而后每1Ma记录一次网格体各相关参数的变化情况。
含油饱和度
图2为模拟10Ma，20Ma，30Ma和 4个时刻的含油饱和度在网格体空间的分布情况。
图2 营11砂体模拟剖面4个模拟时刻的含油饱和度
营11砂体有显著的油气聚集大约从模拟的5～10Ma就已经开始。在地层演化过程中，石油在砂体中一直处于聚集状态，含油饱和度不断升高，这可以从更细致的含油饱和度随时间变化趋势上得以验证。到模拟结束，整个砂体都饱含石油，平均含油饱和度在73%左右，这与砂体实际的含油饱和度平均值（69%）很接近。
油相压力和水相压力
图3给出了在模拟30Ma时间点上油相压力和水相压力在网格体空间的分布情况，而这一时间点呈现的油、水相压力的分布特点基本上代表了整个模拟过程每一时刻的压力分布特点，只是在压力的绝对大小上有差别。网格体油相压力总体变化趋势是由地层的深部向浅部压力逐渐降低，而在这总体背景上，于砂体处存在油相压力的相对低值区。水相压力由地层深部向浅部的变化趋势是逐渐降低的，并且随着网格体埋深总体的压力是增加的。
对网格体毛管压力分布的分析表明，相对低毛管压力区存在于砂岩部位。根据多孔介质中同一点的油相压力和水相压力之差值等于毛管压力可知，油相压力和水相压力分布规律上的差异是由毛管压力的差异引起的。
油势梯度和水势梯度
图4给出了模拟30Ma时间点上油势梯度和水势梯度在网格体空间的分布情况。其中势梯度的正值表明流体流动的方向为轴的负向，而梯度负值表明流体流动方向为轴的正向。
图3 营11砂体模拟剖面在30Ma时油相压力（左图）和水相压力（右图）分布
图4 营11砂体模拟剖面在30Ma时油（上图）和水势梯度（下图）分布
位于左边的两图为水平方向势梯度，位于右边的两图为垂直方向势梯度
4 成藏过程机理分析
营11砂体是处于生油岩包围中的典型透镜状砂岩油藏，其油气来源于围岩生成的烃类。对于这类油藏成藏过程中的油水运移机理和油气聚集过程的认识还存在不足。一般的观点认为异常高压是油气初次运移的主要动力，因此，有些人也笼统地认为异常压力是驱使油气进入砂体的动力。然而，被源岩所包围的砂体内的流体同源岩内流体一样处于封闭环境，而且，在地层沉降压实的过程中，砂体的孔隙也是减小的，因此，从总体上看，砂岩体也是向外排出流体的。因此，如何理解油气自源岩中向砂体运移并聚集，在实际理解上存在一定的困难。
现在普遍的油气运移理论认为，石油是以独立相进行运移的，油水在运移中有着各自独立的流动途径和压力系统，而在同一点的油水压力之差由油水间的毛管压力来平衡。因此，在理解这类透镜状砂体成藏时，不应从单一的流体相来考虑源岩和砂体间的压力差异，而应该像本次模拟一样，将其作为两相流来考虑。
从营11砂体模拟剖面油相压力分布以及油势梯度在水平和垂直方向的变化特点可知，在砂体区存在油的相对于围岩的低势区。油势梯度的正负代表了石油的流动方向，因此砂体区油相低势的特点决定了其必然会成为石油的聚集区。而通过水相压力分布和模拟区水势梯度的变化特点可知在砂体部位不存在水的低势区，砂体对水的流向只起到了一些扰动作用，但水的总体的流动方向是由下向上排出的。
由此可见，超压是推动流体整体运移的动力，而对处于生油围岩包围中的透镜状岩性砂体，围岩和砂体间毛管压力的差异才是驱动油气在其中聚集的根本动力。而这一驱动力是由于围岩和砂体物性上的差异以及油气的生成两种因素结合在一起形成的。
5 结论
（1）通过可压缩多孔介质油水两相渗流的基本原理，并结合与油气的生成、运移和聚集相关的各种因素和作用，可以模拟类似营11砂体的透镜状砂体油藏的成藏过程。
（2）通过对成藏过程中围岩和砂体的油、水相压力及油、水相势梯度的分布特点可知，在成藏过程中砂体区相对于围岩成为油相的低势区，因此石油得以在砂体中进行聚集，而水在砂体中没有聚集的趋势，其总体的运移方向是向着上方的低势区。
（3）石油在类似于营11砂岩油藏中聚集的根本动力是围岩和砂体之间的毛管压力差，而这一差异是围岩与砂体的物性差异以及围岩中石油的生成相结合的必然结果。
参考文献
［1］陈章明，张云峰，韩有信等.凸镜状砂体聚油模拟实验及其机理分析［J］.石油实验地质，1998，20（2）：166～170.
［2］李丕龙，庞雄奇，陈冬霞等.济阳坳陷砂岩透镜体油藏成因机理与模式［J］.中国科学D辑，2004，34（增刊1）：143～151.
［3］王宁，陈宝宁，翟建飞.岩性油气藏形成的成藏指数［J］.石油勘探与开发.2001，27（6）：4～5，8.
［4］庞雄奇，陈冬霞，李丕龙等.砂岩透镜体成藏门限及控油气作用机理［J］.石油学报，2003，24（3）：38～41.
［5］李丕龙，张善文，宋国奇等.断陷盆地隐蔽油气藏形成机制——以渤海湾盆地济阳坳陷为例［J］.石油实验地质，2004，26（1）：3～10.
［6］隋风贵.浊积砂体油气成藏主控因素的定量研究［J］.石油学报，2005，26（1）：55～59.
［7］Palciauskas V migration of petroleum［C］//Merrill R and migration processes and evaluation techniques（AAPG Treatise of Petroleum Geology），Tulas，：13～22.
［8］Mann U，Hantschel T，Schaefer R G，et migration：mechanisms，pathways，efficiencies and numerical simulations［C］//Welte D H，Horsfield B，Baker D and Basin Evolution，SpringerVerlag Berlin Heidelberg，1997：403～520.
［9］罗晓容.油气初次运移的动力学背景与条件［J］.石油学报，2001，22（6）：24～29.
［10］Ungerer P，Burrus J，Dollgez B，et evaluation by integrated two-dimensional modeling of heat transfer，fluid flow，hydrocarbon generation，and migration［J］.AAPG Bulletin，1990，74（3）：309～335.
［11］Düppenbecker S J，Dohmen L，Welte D modeling of petroleum expulsion in two areas of the Lower Saxony Basin，Northern Germany［C］//England W A，Fleet A Migration-Geological Society Special Publication No ：Geological Society，1991：47～64.
［12］Okui A，Siebert R M，Matsubayashi of oil expulsion by 1-D and 2-D basin modelling saturation threshold and relative permeabilities of source rocks［C］//Dppenbecker S J，Iliffe J Modelling：Practice and Progress-Geological Society Special Publications ：Geological Society，1998：45～72.
［13］Luo X R and Guy mechanism of organic matter cracking：Numerical modeling［J］.AAPG Bulletin，（6）：856～874.
［14］Ungerer P，Besis F，Chenet P Y，et and geochemical models in oil exploration：principles and practical examples［C］//Demaison G，Murris R geochemistry and basin evaluation-AAPG Memoir ：AAPG，～57.
［15］Hermanrud C，Wensass L，Teige G M G.，et porosities from well logs on Haltenbanken（offshore mid-Norway）show no influence of overpressuring［C］//Law B E，Ulmishek G F，Slavin V pressures in hydrocarbon environments-AAPG Memoir ：AAPG，～85.
［16］Slavin V I，Smirnova E high formation pressures：origin，prediction，hydrocarbon field development，and ecological problems［C］//Law B E，Ulmishek G F，Slavin V pressures in hydrocarbon environments-AAPG Memoir ：AAPG，～114.
［17］Tissot B P，Welte D Formation and Occurrence［M］.Berlin，Heidelberg，New York，Tokyo：Springer-Verlag，1984.
［18］石广仁.油气盆地数值模拟方法［M］.北京：石油工业出版社，～47.
［19］Danesh A.沈平平，韩冬译.油藏流体的PVT与相态［M］.北京：石油工业出版社，～56.
［20］Tokunaga T，Hosoya S，Tooaka H，et estimation of the intrinsic permeability of argillaceous rocks and the effects on long-term fluid migration［C］//Düppenbecker S J，Iliffe J Modelling：Practice and Progress-Geological Society Special Publications ：Geological Society，～94.
［21］Brooks R H，Corey of porous media affecting fluid flow［J］.Journal of the Irrigation and Drainage Division，1966，92（2）：61～88.
［22］Dullien F A L.范玉平，赵东伟等译.现代渗流物理学［M］.北京：石油工业出版社，2001.
［23］翁庆萍，庞雄奇，Leonard J E.东营凹陷牛庄洼陷沙三段中亚段岩性油藏数值模拟研究［C］.见：李丕龙，庞雄奇.隐闭油气藏形成机理与勘探实践——第三届隐闭油气藏国际学术研讨会论文集.北京：石油工业出版社.2004：215～218.
［24］Mudford B S，Gradstein F M，Katsube TJ，et 1D compaction-driven flow in sedimentary basins：a comparison of the Scotian Shelf，North Sea and Gulf Coast［C］//England W A，Fleet A Migration-Geological Society Special Publication No ：Geological Society，1991：65～85.
［25］Passey Q R，Creaney S，Kulla J B，et practical model for organic richness from porosity and resistivity logs［J］.AAPG Bulletin，1990，74（12）：1777～1794.
［26］Schenk H J，Horsfield B，Krooss B，et of petroleum formation and cracking［C］//Welte D H，Horsfield B，Baker D and Basin Heidelberg，Springer-Verlag，1997：233～269.
康志江 张 允
（中国石油化工份有限公司石油勘探开发研究院，北京 100083）
摘 要：缝洞碳酸型盐岩油藏具有储集空间变化尺度大、介质复杂、流体流动形态多样等特点，无法利 用比较成熟的砂岩油藏数值模拟理论与技术，因此缝洞油藏数值模拟成了当前世界面临的难点和重点，其制 约着这类油藏的合理高效开发。为此，在缝洞油藏尺度上，依据连续性介质的思想框架，发展了双重介质，形成了等效多重介质理论，即将缝洞型油藏中的多相流动问题等效成为若干个连续介质中的多相流动问题，建立了包含溶洞、裂缝、溶孔的三重介质连续性模型，研究了表征单元体理论，提出了模型的建立准则；同 时针对缝洞型油藏大型溶洞中流体流动需要精细刻画的问题提出了耦合型数值模拟技术。主要包括建立了缝 洞型油藏数值模拟多孔介质区、洞穴区及其交界面的数学模型，实现了溶洞中Navier-Stokes流和基质中Darcy 流的耦合，解决了油水两相界面处理问题，形成了洞穴与多孔介质区的交界面条件，然后分别研究了等效多 重介质模型和耦合型数值模拟的数值算法。最后根据形成的缝洞型油藏数值模拟技术编制了的三维三相流体 数值模拟器，通过物理模拟实验和数值模拟实验模拟了一注水驱油过程，结果的一致性验证了方法的正确性。
关键词：缝洞型油藏；数值模拟；多重介质；流渗耦合
Study on Numerical Simulation Technology of Fractured-vuggy Carbonate Reservoir
Kang Zhijiang，Zhang Yun
（Exploration ＆ Production Research Institute，SINOPEC，Beijing 100083，China）
Abstract：Fractured-vuggy carbonate reservoir is characterized by different scales of reservoir space，medium complex，many fluid flow patterns， it can not make use of more mature sandstone reservoir simulation theory and the numerical simulation method of naturally fractured-vuggy carbonate reservoir is the world difficulty and emphasis，and it restricts efficient development of such according to the fractured-vuggy reservoir characteristics，the equivalent multi-media numerical simulation technology was formed based on dual media is，multiphase flow problems are equivalent to a number of multiphase flow problems in continuous medium in fractured- vuggy continuity medium of triple-medium model was established including caves，fractures and so then Representative Elementary Volume was studied，and the model rules were put in order to fine description fluid flow in the large cave of fractured-vuggy reservoir，coupled numerical simulation technology was article established the mathematical model that included porous media area，cave area，and their interface，achieved the coupling of the Navier-Stokes flow in cave and Darcy flow in matrix，and solved the oil-water two- phase interface problem，and formed caves and porous area of the interface the numerical algorithm of the numerical simulation of multiple media model and coupled model was ，the fractured-vuggy reservoir numerical simulator was physical simulation and numerical simulation of a simulation process of water flooding was it was used to verify the correctness of the numerical simulation method.
Key words：fractured-vuggy reservoir；numerical simulation；multi-media；coupled flow
引言
世界上已发现的油气储量有一半以上来自碳酸盐岩油气储集层［1］，而缝洞型碳酸盐岩油藏作为 其中的一种特殊类型，也在我国乃至世界的油气资源中也占有很大的比重。缝洞型碳酸盐岩油藏属 于非常规油气藏类型，其储量规模大，可以形成大型油气藏，也是世界碳酸盐岩油藏生产的重要组 成部分。
近十多年来，研究对象为碎屑岩的油藏数值模拟，其相关的理论与技术研究均基于多孔介质理论，已经取得了巨大的发展，形成了工业化技术应用。但对于储层空间变化尺度大，介质复杂的碳酸盐岩缝 洞型储层，目前的理论与技术方法在很多方面都不能适用，为此开展了缝洞型碳酸盐岩油藏数值模拟研 究，主要为等效多重介质数值模拟技术［2～16］和耦合型数值模拟技术［17～20］。
1 数学模型的建立
模型建立准则
多重介质理论本质上是一种连续介质理论，而连续介质理论成立的前提是其表征单元体存在。目前 在单重介质表征单元体研究方面已有很多成果［21］，对于多重介质表征单元体理论方面的研究国内外还 很少，这是由于复杂介质中不同空隙类型的空间尺度差异很大、空隙中多相流体的流动形态也是多种多 样，因此在我们研究的尺度范围内复杂介质的表征单元体往往并不存在。为解决这一问题，我们提出了 复杂介质多重表征单元体的概念。
对于复杂介质油藏，设ΩK（x0）为复杂介质区域中的一个体积，x0是体积ΩK（x0）的质心，E为 该复杂介质的外延量（质量、空隙空间、单位时间通过的流体质量等）、e为该外延量对应的内涵量（密度、孔隙度、质量流量等）。E（ΩK（x0））表示体积ΩK（x0）内的外延量，eK（x）表示点x处的内涵 量，M、F、V为基质、裂缝、溶洞，如果满足：
国际非常规油气勘探开发（青岛）大会论文集
则，外延量EK相应的内涵量eK的表征单元体存在，连续介质方法可用，采用多重介质方法。否则 就要单独处理，即用离散方法（耦合方法）处理，式（1）和式（2）即是复杂介质的多重介质模型的 建立准则。
多孔介质中的控制方程
洞缝型油藏考虑为等温条件，并且包含油、水两相流体。复杂介质区域流体流动的方程为：
水相：
国际非常规油气勘探开发（青岛）大会论文集
油相：
国际非常规油气勘探开发（青岛）大会论文集
其中，当b相流体（w为水；o为油）为Darcy流动时，其速度根据Darcy定理如下定义：
国际非常规油气勘探开发（青岛）大会论文集
式中，ρβ是β相在油藏条件下的密度； 是在油藏条件下脱去溶解气的油相密度；φ是油层的有效孔 隙度；μβ是β相的黏度；Sβ是β相的饱和度；Pβ是β相的压力；qβ是地层β组分每单位体积汇点/源 点项；g是重力加速度；k是油层的绝对渗透率；krβ是β相的相对渗透率；D是深度。
洞穴中的控制方程
洞穴自由流动区控制方程采用Navier-Stokes方程。在自由流动区域油水不可混溶形成双流体。孔洞 内油水两种流体间有明显的界面，且可以明确表示出来。控制方程包括油区域的控制方程、水区域的控 制方程以及油水界面运动方程。
首先分别对油、水存在的区域给出质量守恒、动量守恒方程。
油相方程：
国际非常规油气勘探开发（青岛）大会论文集
水相方程：
国际非常规油气勘探开发（青岛）大会论文集
其中，fσ表示表面张力。
这就是微可压缩流体的两相流动方程。
作为特例，假设油水不可压缩，则密度为常数，此时在上述方程中消去密度常数可得：
国际非常规油气勘探开发（青岛）大会论文集
以及（8）式和（9）式是两组标准的N-S运动方程组。区别在于在交界面上，两相流体性质如密 度、黏性不同。另外还有交界面的运动方程。
下面考虑油水间界面的运动方程。关于界面有两种表达形式，针对不同的算法可以选取不同的 形式。
(1)界面用点集描述，
国际非常规油气勘探开发（青岛）大会论文集
这种情况下，界面上的点以流体速度按如下规律运动。
国际非常规油气勘探开发（青岛）大会论文集
(2)界面上的点用F（x，t）＝0方程确定
国际非常规油气勘探开发（青岛）大会论文集
此时F（x，t）满足
国际非常规油气勘探开发（青岛）大会论文集
其中u表示流体运动速度。
在油水两种流体的分界面上压力、速度等物理量都是连续的。而密度、黏性等表示流体特性的物理 量则不同。
洞穴与多孔介质区的交界面条件
界面条件包括浓度连续性、压力平衡、流通量平衡等。并考虑油藏实际情况，可以对交界面条件进 行简化。由于不论是多孔介质的压力pd还是洞穴的压力ps，都很大。相对于地层压力，速度和黏性都 很小，因此可以忽略。同时可以假设洞穴流动区域和多孔介质区域在边界切线方向上没有滑移。在这种 假设下，交界面条件可以表示为
国际非常规油气勘探开发（青岛）大会论文集
这组条件实际上表示在交界面上浓度、压力和速度的连续性。
实际计算过程中，交界面条件（12）一般要比较容易使用，特别是在使用有限差分和有限体积进 行离散时。但在使用有限元方法求解推导弱形式过程中，可以直接应用。
2 数值算法研究
等效多重介质模型数值模拟技术
油藏模型考虑为等温条件，并且包含油、气、水三相流体。水和油这两个液体组分分别存在于水相 中和油相中，而气体不仅存在于气相，而且可以溶解于油中。每一相的流体在压力、重力和毛细管力的 作用下按照Darcy定理流动；溶洞内和溶洞之间的流动为非Darcy流或管流。
采用有限体积法进行空间离散后，采用向后一阶差分进行时间离散，可得离散化后单元i内方 程为：
国际非常规油气勘探开发（青岛）大会论文集
其中，M是β相的质量；上标n表示是前一时刻的量；上标n＋1表示是当前时刻的量；Vi是单元 i（基质、裂缝或溶洞）的体积；△t是时间步长；ηi是同单元i相连接的单元j的集合；Fβ，ij是单元i同 单元j之间β相的质量流动项；Qβi是单元i内β相的源汇项。多重介质单元i、j之间的流动项Fβ，ij可表 示如下：
国际非常规油气勘探开发（青岛）大会论文集
其中，λβ，ij＋1/2是β相的流度，
国际非常规油气勘探开发（青岛）大会论文集
为了描述复杂介质中的多种流动形态，在复杂介质的多重介质模型中，单元间的流动分为渗流（达西流、非达西低速流或非达西高速流）、一维管流和裂缝面上的二维流动（非达西高速流）、无充填 溶洞内的三维“洞穴流”。
耦合型数值模拟技术
针对缝洞型油藏大型洞穴内流体流动问题，在Navies-Stokes方程的理论基础上，考虑动量守衡，创 建了油水两相不混溶、微可压缩流动的数学模型，实现了复杂介质油藏Navies-Stokes流和渗流耦合的数 值模拟技术。其数值模型的建立包含两个步骤：求解区域的离散和方程的离散。求解区域的离散产生求 解区域的数值描述，包括求解点的位置和边界描述。空间被分为有限的离散区域，称为控制体积或体 网格。而对瞬态问题，时间区间也被分为有限的时间步长。方程离散则将控制方程的项转化为代数 表达。
对任意的物理量φ，其传递方程可写为：
国际非常规油气勘探开发（青岛）大会论文集
有限体积方法要求满足以P为基础的控制体VP中控制方程的积分形式：
国际非常规油气勘探开发（青岛）大会论文集
由于扩散项是φ的二阶导数，为保证一致性，有限体积中离散的阶数必须等于或大于方程的阶数。
离散方法的精确性取决于在点P附近的时空位置上假定的变化函数φ＝φ（x，t）。
要获得传递方程的离散形式，关键在于交界面f上的值及其上的垂直梯度，即φf和S·（▽φ）f。对 位于区域边界上的面，其值由边界条件计算得到。
3 油藏数值模拟方法的验证
根据对油藏数值模拟方法研究结果，编制了相应的数值模拟软件，为了验证该数值模拟方法的正确 性，开展了注水驱替油物理实验，实验中充填物为右半部为5mm白色大理石，左半部为3mm白色大理 石，注入清水速度为，模型内部充满油，从左向右驱替油，实验结果如图1所示，采用相同 的参数进行数值模拟，结果如图2所示，通过比较可以得出数值模拟实验与物理实验趋势一致，从而验 证了方法的正确性。
图1 水驱油物理实验现象
图2 油藏数值模拟实验含油饱和度图
4 结论
（1）双重介质理论在裂缝型油藏广泛应用，较好地解决了基质与裂缝中流体流动差异性大的问题，对于小型溶蚀洞，也有的专家开展了三重介质数值模拟研究，对于具有洞穴的缝洞型油藏没有相关报 道，通过研究形成了一套能处理洞穴的基于多重介质缝洞油藏自适应隐式数值模拟方法；
（2）针对洞穴内两相流界面计算和洞穴多孔介质耦合计算两个关键问题开展了研究，结合洞穴内 油水两相流物理实验，解决了油水两相界面技术问题，形成了耦合油藏数值模拟方法，确定了数值 解法；
（3）根据形成的数值模拟方法编制了相应的数值模拟软件，并通过实验验证了方法的正确性。
参考文献
［1］李珂，李允，刘明.缝洞型碳酸盐岩油藏储量计算方法研究［J］.石油钻采工艺，2007，29（2）：103～104.
［2］，， simulation of water-oil flow in naturally fractured of Reservoir Performance［C］.Loa Angeies，1976.
［3］ of naturally fractured reservoirs with semi-implicit source SPE-AIME Fo ～h SymPOsiVm on Numerical Simulation of Reservoir Performance［C］.Los Angeles，1976.
［4］ of naturally fractured Reservoir Simulatlon Symposium［C］.San Francleco，CA，1983.
［5］， new approach for simulating complex fractured SPE 1985 Middle East Oil Technical Conference and Exhibition［C］.Bahrain：SPE 13537，1985.
［6］，， simulation of naturally fractured 1986 SPE Annual Technical Conference and Exhibition［C］.New Orlean：SPE 15627，1986.
［7］ matrix，fracture fluid transfer function in dual porosity SPE International Meeting on Petroleum Engineering［C］.Tianjin，China：Paper SPE 17599，1988.
［8］Abbas Firoozabadi， spe comparative solution project：dual-porosity 1989 SPE Symposium on Reservoir Simulation［C］.Houston：SPE 18741，1989.
［9］Karimi-Fard，M.，Gong，B.，Durlofsky， of coarse-scale continuum flow models from detailed fracture .［J］.42，2006.
［10］Yu-Shu Wu，Guan Qin，Richard et Multiple-Continuum Approach for Modeling Multiphase Flow in Naturally Fractured Vuggy Petroleum 2006 SPE International Oil ＆Gas Conference and Exhibition［C］.Beijing，China：SPE 104173，2006.
［11］，， et multiphase flow in naturally fractured vuggy petroleum 2006 SPE Annual Technical Conference and Exhibition［C］.Texas，：SPE 102356，2006.
［12］σmez，， et of an Evolutionary Algorithm in Well Test Characterization of Naturally Fractured Vuggy First International Oil Conference and Exhibition［C］.Mexico：SPE 103931 ，2006.
［13］Peter Popov，Linfeng Bi，Yalchin Efendiev et and Multiscale Methods for Modeling Fluid Flow Through Naturally Fractured Vuggy Carbonate 15th SPE Middle East Oil ＆ Gas Show and Conference［C］.Kingdom of Bahrain：SPE 105378，2007.
［14］Wu，.， transient flow in naturally fractured reservoirs with three-porosity systems［J］.Acta，Mechanica Sinica，1983，15（1）：81-85.
［15］Abdassah， Ershaghis I..Triple-porosity system for representing naturally fractured reservoirs［C］.SPE .，1986，1，113-127.
［16］Wu，.，， triplecontinuum approach for modeling flow and transport processes in fractured rock［J］.Journal Of Contaminant Hydrology，2004a，73：145-179.
［17］Kang，Zhijiang，Yu-Shu Wu，Jianglong Li，Yongchao Wu，Jie Zhang，Guangfu multiphase flow in naturally fractured vuggy petroleum 2006 SPE Annual Technical Conference and Exhibition［C］.San Antonio，Texas，2006，24-27.
［18］林加恩，李亮，杨慧珠.管流与渗流耦合流动理论研究初探［J］.西安石油大学学报（自然科学版），2007，22（2）：11～12.
［19］耿克勤，吴永平.拱坝和坝肩岩体的力学和渗流的耦合分析实例［J］.岩石力学与工程学报，1997，16（2）：125～131.
［20］耿克勤.复杂岩基的渗流、力学及其耦合分析研究以及工程应用［D］.清华大学工学博士论文.1994.
［21］郭自强.碳酸盐岩缝洞型油藏水锥预测及控制研究［D］.北京：中国石油大学，2009.
［22］周创兵，陈益峰，姜清辉.岩体表征单元体与岩体力学参数［J］.岩土工程学报，2007，29（8）：1135～1142.
毕业论文数据分析的做法如下：
首先，针对实证性论文而言，在开始撰写论文之前，必须要提前确定好数据研究方法。而数据研究方法的确定与选择需要根据大家毕业论文的研究课题来确定。
另外，大家也可以跟自己的的论文指导老师多多交流，尽可能多的了解更多关于研究方法的知识，以供自己选择。除此之外，大家还需要大量查找文献资料，见多识广有大量输入之后才能有所输出，本环节需要大家跟导师沟通商议后决定。
接下来一个比较重要的步骤是搜集和整理实验数据。在这一部分，很多同学朋友都会遇到各种各样的问题，比如，不知道去哪里找数据，找到的数据可靠性无法保障，需要的数据总是无法搜集全面等等各种问题。
那么在这里需要跟大家强调一下，推荐大家使用国家统计局、中国统计年鉴、国泰安、万方等等这些比较权威的网站去搜集数据资料。
在此需要注意的是，国泰安和万方等这些网站是需要收费的，上去看了一下，价格不是很亲民。
给大家分享一下，如果有些数据在国家官方网站确实找不到或者毕业论文所需的最新数据还没及时发布，推荐大家可以上某宝，因为某宝上电子版数据往往都很全面，而且价格大都可以接受。
在此提醒大家搜集到数据之后，一定要按照自己的习惯整理保存好，避免后期使用数据时出现差错。
通过数据进行分析的论文用数据是数学方法。
数据分析方法：将数据按一定规律用列表方式表达出来，是记录和处理最常用的方法。表格的设计要求对应关系清楚，简单明了，有利于发现相关量之间的相关关系。
此外还要求在标题栏中注明各个量的名称、符号、数量级和单位等：根据需要还可以列出除原始数据以外的计算栏目和统计栏目等。
数据分析目的：
数据分析的目的是把隐藏在一大批看来杂乱无章的数据中的信息集中和提炼出来，从而找出所研究对象的内在规律。在实际应用中，数据分析可帮助人们做出判断，以便采取适当行动。数据分析是有组织有目的地收集数据、分析数据，使之成为信息的过程。
这一过程是质量管理体系的支持过程。在产品的整个寿命周期，包括从市场调研到售后服务和最终处置的各个过程都需要适当运用数据分析过程，以提升有效性。
例如设计人员在开始一个新的设计以前，要通过广泛的设计调查，分析所得数据以判定设计方向，因此数据分析在工业设计中具有极其重要的地位。
论文常用数据分析方法
论文常用数据分析方法，对好的论文分析研究方法应该从哪些方面展开，如何表达才能显得自己对该论文真的有所理解，应该看哪些书呢？下面我整理了论文常用数据分析方法，一起了解看看吧!
论文常用数据分析方法分类总结
1、 基本描述统计
频数分析是用于分析定类数据的选择频数和百分比分布。
描述分析用于描述定量数据的集中趋势、波动程度和分布形状。如要计算数据的平均值、中位数等，可使用描述分析。
分类汇总用于交叉研究，展示两个或更多变量的交叉信息，可将不同组别下的`数据进行汇总统计。
2、 信度分析
信度分析的方法主要有以下三种：Cronbach α信度系数法、折半信度法、重测信度法。
Cronbach α信度系数法为最常使用的方法，即通过Cronbach α信度系数测量测验或量表的信度是否达标。
折半信度是将所有量表题项分为两半，计算两部分各自的信度以及相关系数，进而估计整个量表的信度的测量方法。可在信度分析中选择使用折半系数或是Cronbach α系数。
重测信度是指同一批样本，在不同时间点做了两次相同的问题，然后计算两次回答的相关系数，通过相关系数去研究信度水平。
3、 效度分析
效度有很多种，可分为四种类型：内容效度、结构效度、区分效度、聚合效度。具体区别如下表所示：
4、 差异关系研究
T检验可分析X为定类数据，Y为定量数据之间的关系情况，针对T检验，X只能为2个类别。
当组别多于2组，且数据类型为X为定类数据，Y为定量数据，可使用方差分析。
如果要分析定类数据和定类数据之间的关系情况，可使用交叉卡方分析。
如果研究定类数据与定量数据关系情况，且数据不正态或者方差不齐时，可使用非参数检验。
5、 影响关系研究
相关分析用于研究定量数据之间的关系情况，可以分析包括是否有关系,以及关系紧密程度等。分析时可以不区分XY，但分析数据均要为定量数据。
回归分析通常指的是线性回归分析，一般可在相关分析后进行，用于研究影响关系情况，其中X通常为定量数据（也可以是定类数据，需要设置成哑变量），Y一定为定量数据。
回归分析通常分析Y只有一个，如果想研究多个自变量与多个因变量的影响关系情况，可选择路径分析。
、首先我们对上表数据进行细化，找到每组内受访者的具体满意度打分数值，而不是这个汇总后的得分值。2、SPSS方差分析：image分析：比较均值，单因素方差分析因变量列表：品类满意度因子：收入选项：方差同质性检验3、数据是否适合做方差分析image方差分析之前，需要进行可行性检验，原假设，各分组方差无差异。根据同质性检验可知，sig值，为大概率，原假设成立，即不同分组之间同质，没有显著差异，可进行方差分析。4、方差分析结果image原假设，各分组之间无差异。方差分析sig值，大于小概率值，为大概率，原假设成立，即不同收入水平分组之间在品类满意度上并不没有不同。不存在显著差异。5、用可视化图来揭示原因image我们可以看到，每类收入者的满意度得分都围绕平均值上下波动，这表明不同收入者对品类的态度存在明显差异，例如，同是高收入者，有的非常满意，有的却十分的不满意。同组内的差异甚至高出不同收入者之间的差异，这一点可以通过方差分析中方差得以判断。因此说，收入水平并不是导致用户对A卖场品类满意度的关键因素。可见，数据的表象往往迷惑人，尤其是综合汇总后的平均值，通过对底层数据进行分组及方差分析则可以让我们拨开云雾，看到数据的本质。同时，这个案例也告诉我们，在常规的报表分析当中，经常性的工作是对底层数据进行汇总分析，然后拿汇总数据用于决策，此时，非常容易就数字大小的对比而做出判断，报表工作人员需要注意，需要养成用统计的理念和逻辑上报数据的结果。打开CSDN，阅读体验更佳Stata：多个变量组间均值\中位数差异检验 作者：韩少真(西北大学) || 刘婉青(西北大学) Stata 连享会： 知乎 | 简书 | 码云 | CSDN 2019暑期Stata现场班，日，北京，连玉君+刘瑞明 主讲 Stata连享会 精品专题 || 精彩推文 文章目录1. 问题背景 期刊论文示例一 期刊论文示例二 期刊论文示例三2. Stata实现组间均值或中位数差异检验的常见...浏览器打开 均值已知检验方差_SPSS篇—方差分析昨天跟大家分享了如何用SPSS进行回归分析，知道了回归分析的用途以及使用的场景。今天跟大家分享的就是之前文章里面出现很多次的一个分析—方差分析。方差分析又被称作“F检验”或者“变异数分析”，主要是用于两个及两个以上样本均值差异的显著性检验。方差分析和回归分析一样，也有很多个分支。对于方差分析，一般我们是用来研究不同来源的变异对总变异的贡献大小，从而确定可控因素对因变量的影响大小。我们今天通过一个例...浏览器打开相关推荐双因素方差分析_科研常用显著性分析方法汇总及选择（T检验，Mann-Whitney U test检验，方差分析等）...科研常见的差异性分析方法汇总根据数据是否符合正态分布，分为：参数检验非参数检验非参数检验是在总体方差未知或知道甚少的情况下，利用样本数据对总体分布形态等进行推断的方法。由于非参数检验方法在推断过程中不涉及有关总体分布的参数，因而得名为“非参数”检验。 参数检验（parameter test）全称参数假设检验，是指对参数平均值、方差进行的统计检验。先由测得的样本数据计算检验统计量，若计算的统计量值落...浏览器打开如何判断组之间是否有显著性差异？怎么知道组之间是否有显著性差异？ 方法：单因素方差分析；双尾检验；K-S检验；x²检验；蒙特卡罗检验 1 K-S检验法介绍： 有人首先想到单因素方差分析或双尾检验（2 tailed TEST）。其实这些是不准确的，最好采用Kolmogorov-Smirnov test（柯尔莫诺夫-斯米尔诺夫检验）来分析变量是否符合某种分布或比较两组之间有无显著性差异。（...浏览器打开均值已知检验方差_方差分析与R1.什么是方差分析？假设有多个总体（三个及以上），都是服从正态分布且方差相同。方差分析就是检验多个总体均值是否相等的统计方法。比如用三种鸡饲料喂小鸡，三个月后小鸡的重量是随机的，假设服从正态分布。我们自然就问，这三种鸡饲料喂的小鸡三个月以后重量的均值是否相同？从这个例子中我们可以看出，在假设其它条件相同的情况下，造成小鸡三个月后平均重量不同的因素就是鸡饲料。若三种鸡饲料对小鸡重量的影响效果相同，那...浏览器打开均值已知检验方差_方差分析不显著就一定无差异吗？方差分析的零假设是：各组均值相等。这个“各组均值相等”如何理解？正确理解是：各组和所有组总均值相等，并不是真的“各组均值相等”。方差分析认为：各组和总均值无差异，那么各组均值等于总均值，意味着各组均值相等。单因素方差分析大家应该都理解的比较好，我们可以看看单因素方差分析F检验统计量的分子核心部分：（各组均值-总均值）的平方。看到没，减的是“总均值”。一般来说，如果各组和总均值无差异，
是数据处理时的错误。如果是硕士毕业论文，这个只是数据处理时的错误。除非你原文是有颠覆性的结果且改正后就没了，也没有人会在意你数据处理时出现的错误。
1、方差分析前提：不同水平下，各总体均值服从方差相同的正态分布。
2、方差齐性检验：采用方差同质性检验方法（Homogeneity of variance），在spss中打开你要处理的数据，在菜单栏上执行：analyse-compare means--one-wayanova，打开单因素方差分析对话框 在这个对话框中，将因变量放到dependent list中，将自变量放到factor中，点击post hoc，选择snk和lsd，返回确认ok。
SPSS是一个应用非常广泛的统计分析软件，它已有40余年的成长历史，已经成为全球专业统计分析软件的佼佼者。对于统计行业的从业人员来说，学会SPSS是其必备的基本技能之一。本书将带领读者进入SPSS，让SPSS帮助读者轻松完成统计分析工作。
操作设备：戴尔电脑
操作系统：win10
操作软件：SPSS 版
1、首先打开SPSS 版软件，找到要编辑的数据，可以从下图中找到方框。
2、在接下来的过程中需要在上方菜单栏中找到分析菜单，将鼠标移动到一般线性模型，然后选择单个变量，单击鼠标左键选择。
3、可以看到界面中的红色框。在单变量对话框中，将变量分别移动到因变量和协变量。在这里，将高度移动到因变量，将药物移动到协变量。
4、单击右侧菜单中的选项，将鼠标移动到单变量选项，选择参数估计值，并将参数估计值标记为勾号。
5、选择完成后，点击选项中的继续选项，然后可以选择在单变量对话框中点击确定，即可查看编辑后的操作。
6、最后可以看到界面上的方框显示在SPSS查看器中可以看到药物对身高影响的显着性分析，红框内的显着性为0<，为显着。
t检验
适用于计量资料、正态分布、方差具有齐性的两组间小样本比较，检验两个处理平均数的差异是否显著。
spss提供的T检验有3种形式，分别是单样本T检验（One-Sample T Test），独立样本T检验(Independent-Sample T Teat)和成对样本T检验(Paired-Sample T Test）。