研究方法二

1. 研究方法及流程

作者采用如下方法和步骤（图1-7），以保证顺利完成主要研究内容，达到发展煤层气地质基础理论，服务于煤层气勘探开发实践的最终目的。

第一，查阅前人有关文献、区域地质调查资料和地质勘探成果，分析现有成果的研究程度以及有待于深化的主要问题，制定本文研究计划，明确研究思路和主攻方向，提出解决方法和研究流程。

图1-6山西南部二叠系山西组上主煤层含气量等值线图

第二，开展野外地质和矿井地质调查，收集现场资料，进行样品采集，同时对资料进行初步的分析整理，以进一步增进对区内煤层、煤级、瓦斯、煤层物性等特征的了解，并为后续的室内研究工作奠定坚实基础。

第三，进行样品的分析测试，包括煤岩学分析、煤层含气性分析、煤孔隙性分析、煤的顺磁共振实验、岩样X射线衍射分析、脉体流体包裹体矿物学和地球化学测试等，目的是为深化研究提供可靠的数据资料。

第四，进行地质调研、实验数据的整理分析及数值模拟软件的编制工作，进而查明研究区地质发展的四方面的历史，即煤层沉积埋藏史、盆地热历史、煤的热演化史以及煤层气生成保存史，特别是与煤层气资源评价密切相关的二次生烃史及其空间分布规律。

图1-7研究方法及流程示意图

第五，综合分析上述研究成果，从煤中有机质生烃演化角度，对本区煤层气资源前景作出评价和预测，为煤层气地质选区和勘探开发部署提供理论依据。

2. 论文的研究方法一栏该怎么写

论文的研究方法主要有以下几种：

1、调查法

它是有目的、有计划、有系统地搜集有关研究对象现实状况或历史状况的材料的方法。调查方法是科学研究中常用的基本研究方法，它综合运用历史法、观察法等方法以及谈话、问卷、个案研究、测验等科学方式，对教育现象进行有计划的、周密的和系统的了解。

2、观察法

观察法是指研究者根据一定的研究目的、研究提纲或观察表，用自己的感官和辅助工具去直接观察被研究对象，从而获得资料的一种方法。

3、实验法

实验法是通过主支变革、控制研究对象来发现与确认事物间的因果联系的一种科研方法。其主要特点是：第一、主动变革性和控制性。

4、文献研究法

文献研究法是根据一定的研究目的或课题，通过调查文献来获得资料，从而全面地、正确地了解掌握所要研究问题的一种方法。

5、实证研究法

在科学研究中，通过定量分析法可以使人们对研究对象的认识进一步精确化，以便更加科学地揭示规律，把握本质，理清关系，预测事物的发展趋势。

3. 研究方法有哪些

研究方法，哲学术语，是指在研究中发现新现象、新事物，或提出新理论、新观点，揭示事物内在规律的工具和手段。这是运用智慧进行科学思维的技巧，一般包括文献调查法、观察法、思辨法、行为研究法、历史研究法、概念分析法、比较研究法等。

研究方法是人们在从事科学研究过程中不断总结、提炼出来的。由于人们认识问题的角度、研究对象的复杂性等因素，而且研究方法本身处于一个在不断地相互影响、相互结合、相互转化的动态发展过程中，所以对于研究方法的分类很难有一个完全统一的认识。

4. 研究方法是什么

过点A作直线a和b的垂直线

5. 课题研究的基本方法有哪些

课题研究的基本方法有行动研究法、资料收集法、学生带动法、教育实验法、个案研究法。

1、行动研究法：制定个性研究方案，通过学生实践情况进行分析，再研究调整重新进行实践。并将经验总结、记录，形成有价值的文字。

2、资料收集法：深入班级，深入学生个体，对学生现状进行调查，利用不同的资源进行收集，找准问题所在，明确研究对象。

3、学生带动法：通过一小部分学生先学、先走，在带动、感染他周围的学生也来学习。

如要详细、全面拍摄一堂课，一部摄像机是不够的。观察者应准备几部摄像机，并事先作好分工。即使是作观察记录，也需要事先作好设计。在记录纸上印好以一定的格式排列的必须记录的项目，还可以约定一些记录符号，以尽量减少现场记录时书写文字的时间。

6. 研究方法

地下水系统演化:即在自然条件与人类活动影响叠加状态下地下水动力场及地下水化学场的演化。本书主要研究内容是地下水系统演化和地下水水环境保护。因此，本书综合采用水文地质、水文地球化学、环境地质学、数学等研究方法。以演化的观点为指导，动态地研究地下水系统的变化特征;从整体的观点出发，研究地下水系统中地下水动力场、化学场的相互关系;以环境的观点为依托，系统地研究地下水主要超标因子的来源与形成机理，并进行地下水系统的脆弱性分区，为地下水环境的改善和保护提供可靠的信息和科学的依据。具体方法如下:

图1.1 研究技术路线框图

1)用MapGIS和SUFER8.0刻画了研究区不同时期地下水流场、地下水位降落漏斗、地下水水化学场，以及主要因子的动态变化和空间分布特点，阐明了区域地下水系统的演化。

2)用水文地球化学模拟软件PHREEQC对地下水Fe离子升高的机理进行了探讨。

3)用DRASTIC模型和GIS平台的空间分析功能对研究区进行地下水系统内在脆弱性分区，为保护水资源提供理论依据。

首先，在查阅国内外大量专业文献，了解本领域研究现状的基础上，根据存在的环境水文地质问题，结合目前环境水文地质研究特点和发展方向，紧紧围绕地下水系统演化和地下水环境保护两条主线，充分收集资料，然后进行环境水文地质条件综合分析;最后在分析研究地下水动力场、地下水化学场演化规律及其影响因素的基础上，充分分析影响地下水水质的因素，探讨了总硬度升高、Fe离子含量升高的机理，并对研究区进行了地下水脆弱性分区。为保护水环境以及改善水质提供理论依据。研究方法与研究思路见图1.1。

7. 　主要研究方法

研究金属矿床成矿时代的常用方法有三种，一是矿石铅同位素年代学方法，二是蚀变矿物的同位素测年方法，三是据赋矿围岩、控矿构造及与矿化有关岩脉的时代间接推断矿脉形成时代。本书主要应用这三种不同的年代学方法确定矿床成矿时代，同时注意不同方法所得年龄的对比分析与相互验证。近年来发展起来的铼-锇同位素年代学方法能直接测定辉钼矿等矿石矿物的形成时代，然而这种方法在我国目前尚处在试用阶段，在燕山地区尚未全面展开该项测年工作。

一、普通铅同位素的演化模式与年龄计算公式

矿石铅同位素年代学方法是直接测定成矿时代的重要研究方法，被广泛用于世界各地的金属矿床。目前常用的铅同位素演化模式包括单阶段模式如Holms-Houtermans模式，二阶段模式如正常铅混合模式、瞬间增长模式与连续增长模式，多阶段模式如简单的三阶段铅混合模式等。但这些模式都存在严格的应用条件。单阶段模式只适合于封闭体系、无后期铅混染的少数几个整合矿床；简单的二、三阶段模式要求体系相对封闭，各阶段异常铅只能来自于单一的且铀、钍、铅同位素比值均一的源区，还要求体系在各阶段的铅同位素均匀分布。这些模式在一般的造山带与地盾、地台区，都能有效地用于确定矿床成矿时代。然而，燕山陆内造山带具有十分复杂的地质过程，矿质具有两种以上的复杂来源；成矿体系多属开放体系，铀-钍-铅同位素混合过程也颇为复杂，存在多种不同的情况；上述几个特殊的铅同位素模式不足以概括本区常见的开放体系铅的混合过程，以至于使本区已积累的近百组铅同位素资料长期以来得不到充分利用，求不出有地质意义的成矿时代。为此，笔者首先从理论上分析常见开放体系铅同位素混合过程，建立开放体系铅同位素演化模式，推导其年龄计算公式。这些模式在燕山地区成岩成矿时期的研究中，取得了良好的应用效果。

1.基本假设

（1）同一来源的²⁰⁶Pb、²⁰⁷Pb、²⁰⁸Pb、²⁰⁴Pb以相同的概率进入同一样品。不同铅同位素化学性质的相似性，使这一假设在各种地质过程中都能成立。

（2）同一时代地质体的N（²³⁸U）/N（²⁰⁴Pb）（即μ值）与N（²³⁵U）/N（²⁰⁴Pb）（v值）可以变化；铀的丢失与加入常造成这种结果。

（3）当铅混合时，铅同位素可来源于两种以上不同的铅源，包括正常铅铅源与放射成因异常铅铅源；同一铅源对不同样品的贡献可以不一样，即同一体系不同样品的铅同位素来自于任一源区的概率可以不一样。

（4）体系中的铅可以来自于一个至数个放射性成因铅源，将N（²³⁸U）/N（²⁰⁴Pb）＝μ_i的源区叫μ_i源。

（5）铀、铅及其同位素在地幔中均匀分布。

（6）铅在最后一阶段混合后，保持其同位素比值，直至现代。

2.二阶段铅混合的系统模式

设样品来自于t₁时形成正常铅的概率为α₁，来自于T至t₂时期形成的放射成因铅的概率为α₂。t₂混合时，设有m个μ_i源，样品中混合铅来自于μ_i源的概率为β_i。t₂混合之后，样品铅同位素组成可表示为：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

式中：

为第二阶段（t₂）体系的铅同位素组成；

，

为第一阶段（t₁）体系铅同位素组成，由H-H模式确定：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

a₀、b₀为T=4550Ma时地球的初始铅同位素组成；α₁＋α₂=1，

；T为地球年龄。

模式Ⅰ当α₁=1，α₂=0时，由（3.1.1）、（3.1.2）式知，二阶段铅退化为单阶段铅。这时为正常铅，样品点在N（²⁰⁷Pb）/N（²⁰⁴Pb）—N（²⁰⁶Pb）/N（²⁰⁴Pb）坐标图中分布于一点。据（3.1.3）、（3.1.4）式得：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

由（3.1.5）式与（3.1.3）、（3.1.4）式可计算成岩或成矿年龄t₁与源区μ，v值。该模式相当于H-H模式（G.福尔，1983）。

模式Ⅱ0＜α_i＜1,i=1,2；β₁=1，β_j=0（2≤j≤m），μ₁＝μ；这时（3.1.1）,（3.1.2）式可简化为：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

令

由（3.1.6）、（3.1.7）式得：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

当α_i对不同样品取值不一样时，样品点呈线性分布，直线斜率为R，如图3-1所示。样品点分布于增长曲线的弦上，等时线与增长曲线的两交点对应时代t₁与t₂相当于两次普通铅的形成时代。该模式相当于前述已有的正常铅与正常铅混合二阶段模式。当已知t₁与t₂之一时，可据R求出另一时代。

模式Ⅲ－1当i=1,2时，0＜α_i＜1,0＜β_j＜1（1≤j≤m），β_j及

（β_j·μ_j）对不同样品不取恒定值，但α_i对所有样品取恒定值。这时，由（3.1.1）、（3.1.2）式导出：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

（3.1.8）式中，

为混合铅同位素比值。这时样品点分布在一条直线上，直线斜率较大（图3-2），据（3.1.8）式能求出t₂。当其它条件相同，而β_j对所有样品取定值（1≤j≤m）时，由（3.1.1）、（3.1.2）式可知，样品点的铅同位素组成均匀分布，在坐标图中分布于一点；在这种情况下，难以求出t₁或t₂值。

模式Ⅲ－20＜α_i＜1,i=1,2;0≤β_j＜1,1≤j≤m；α₁对不同样品皆非恒定值，β_j对不同样品非定值；这时，若

（β_j·μ_j）趋于μ，则由（3.1.1）、（3.1.2）式导出：

图3-1模式Ⅱ图解

Fig.3-1Lead-isotope evolution of modelⅡ

图3-2模式Ⅲ-1图解

Fig.3-2Lead-isotope evlution of modelⅢ-1

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

令

由于

（β_j·μj）趋于定值μ，所以X′_t1与Y′_t1近为定值。代入（3.1.9）、（3.1.10）式，得：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

这时，样品点呈线性分布，据直线斜率能求出t₁与t₂之一。

令

，则由（3.1.9）、（3.1.10）式可以看出，当μ′＜0时，样品点靠近t₁分布，甚至会落在t₁左侧；当μ′≥0时，样品点靠近t₂点分布，部分样品点会落在t₂右侧。增长曲线如图3-3。当t₁与t₂相差较大时，该模式相当于连续增长模式；当t₁与t₂近似相同时，则等时线由弦而渐趋于切线，这时相当于瞬间增长模式。

模式Ⅲ-3当0＜α_i＜1,0≤β_j＜1（i=1,2,1≤j≤m），β_j、α_j对不同样品皆非常数时，若样品的α₁值仅取几个定值之一，当样品点足够多时，样品点呈图3-4所示分布状态，即分布于一组平行直线上。据直线斜率能求出t₁或t₂，斜率R可表示为：

图3-3模式Ⅲ-2图解

Fig.3-3The first lead-isotope evolution of model Ⅲ-2

图3-4模式Ⅲ-3图解之一

Fig.3-4The first lead-isotope evolution of model Ⅲ-3

当

（β_j·μ_j）为定值时，

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

当

（β_j·μ_j）不为定值时，

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

若α_i对不同样品都不一样时，样品点呈星散状分布（图3-5），这时无法求出t₁或t₂的真实值。

图3-5模式Ⅲ－3图解之二

Fig.3-5The second lead-isotope evolution of model Ⅲ－3

3.三阶段铅混合的系统模式

设一阶段铅的分离时代为t₁，二阶段铅的混合时代为t₂，三阶段铅的混合时代为

、

为二阶段普通铅源i的同位素比值，

、

为t₃体系中铅同位素比值；设有m个放射成因铅源μ_i，有n个普通铅源；t₃时刻混合时，体系铅来自于普通铅i源的概率为ε_i，来自于放射成因铅的概率为ε_n+1；当ε_n+1＞0时，μ_j源铅进入样品的概率为β_j，则

=1，且

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

（3.1.12）、（3.1.13）式为一般情况下三阶段铅混合时的定量关系式。不同条件下，三阶段混合铅具有不同特征，对应于不同的铅演化图，下面分别予以讨论。

（1）ε₁=1，ε_i=0,2≤i≤n+1，这时三阶段铅退化为二阶段铅。

（2）0＜ε₁＜1；ε_i=0,2≤i≤n;0＜ε_n+1＜1，这时（3.1.12）、（3.1.13）式可写成：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

模式Ⅳ当β₁=1，β_j=0,2≤j≤m时，（3.1.14）与（3.1.15）式可写成：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

若放射性铅与普通铅在T到t₂期间有相同的演化过程和成分，即

则

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

这时相当于G.福尔提出的简单三阶段模式；且样品点或呈线性分布（图3-6），或分布于一点。据等时线斜率R能求出t₂与t₃之一：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

模式Ⅴ－1当所有样品的ε₁、X_t2、Y_t2取相同值时，则ε₁·X_t2、ε₁·Y_t2为常量。若β_j对所有样品取相同值，0≤β_j≤1,1≤j≤m；这时三阶段样品铅同位素构成一点。据（3.1.14）、（3.1.15）式，有

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

只有当ε₁及X_t2、Y_t2都已知时，才能求出t₃；一般情况下，若上述三参数未知，则无法计算出真实年龄t₃。

模式Ⅴ-2当ε₁及X_t2、Y_t2为常量，而不同样品β_j不同时，1≤j≤m，若

不为恒定值，则据（3.1.14）、（3.1.15）式，可推导出：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

这时，样品点呈线性分布（图3-7），直线斜率一般较大。据R能求出t₃。

图3-6混合铅模式Ⅳ图解

Fig.3-6Lead-isotope evolution of model Ⅳ

图3-7模式Ⅴ－2图解

Fig.3-7Lead-isotope evolution of model V-2

模式Ⅵ当所有样品点的X_t2、Y_t2恒定时，若0≤β_j≤1,1≤j≤m，β_j非常数；而X_t2

，则据（3.1.14）、（3.1.15）式，可导出：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

令

则

、

近为常量。据（3.1.16）、（3.1.17）式可导出：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

这时，样品点呈线性分布，分布特征类似于模式Ⅲ-2，如图3-8所示。

模式Ⅶ当X_t2、Y_t2恒定，ε₁、β，对不同样品取不同值时，若

不恒定，且ε₁仅有几个可能的值，则混合铅样品点分布于几条平行直线上，直线斜率

据之能求出t₃，否则，样品点呈星散状分布。混合铅演化如图3-9所示。

图3-8模式Ⅵ图解

Fig.3-8Lead-isotope evolution of model Ⅵ

（3）当不同样品的X_t2、Y_t2不同，0≤ε_i＜1,1≤i≤n+1时，有下列模式：

模式Ⅷ若X_t2、Y_t2呈线性分布，不同样品点ε_i相同（1≤i≤n），0≤β_j＜1（1≤j≤m）；则有几种可能性：

模式Ⅷ－1若β_j恒定，1≤j≤m，则（3.1.12）、（3.1.13）式可写成：

图3-9模式Ⅶ图解

Fig.3-9Lead-isotope evolution of model Ⅶ

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

令

由于β_i为常量，对所有1≤j≤m都成立，所以

c与d皆为常数；样品点仍呈线性分布，其斜率与二阶段等时线相同，如图3-10示。据R能求出t₁与t₂之一，而求不出t₃。

图3-10模式Ⅷ－1图解

Fig.3-10Lead-isotope evolution of modelⅧ－1

模式Ⅷ－2若

相同，但β_j不同时，则任一（X_t2，Y_t2）点都对应一条三阶段等时线，所有样品点沿两组平行直线分布（图3-11），r₁一般大于r₂；

，据之能求出t₂与t₃之一；r₁为二阶段等时线斜率，据之能求出t₁与t₂之一。只有当样品点足够多时，才有可能据该模式求出t₁、t₂或t₃，否则，r₁与r₂难以确定，无法计算年龄。

图3-11模式Ⅷ－2图解

Fig.3-11Lead-isotope evolution of modelⅧ－2

模式Ⅸ若（X_t2,Y_t2）呈线性分布，不同样品ε_j值相同，β_j值不同，

亦因样品不同而不同，这时（3.1.12）、（3.1.13）式中

为常数，由（3.1.12）、（3.1.13）式可导出：

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

这时样品点沿两组斜率较大的平行直线分布。当样品点足够多而能求出r₁，与r₂时，则可据此求出t₁、t₂或t₃。

模式X若（X_t2,Y_t2）呈线性分布，但ε_i，β_j对不同样品不取恒定值时，则据（3.1.12）、（3.1.13）式，样品点呈星散状分布，或呈线性沿两组平行直线分布。后一种分布状状只有当ε_i对不同样品点仅取几组确定值时才能出现，据平行直线的斜率能求出t₃，斜率r₂为：

当

非定值时

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

当

恒定值时

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

模式Ⅺ当（X_t2,Y_t2）不呈线性分布，而呈星散状分布时，则三阶段铅样品点仍呈星散状分布，这时无法求出t₃与t₂的真实值。

模式Ⅻ当（X_t2,Y_t2）分布于数条平行直线上，而β_j、ε_j恒定时，由（3.1.12）、（3.1.13）式可得出样品点的（X_t3,Y_t3）仍呈线性分布，斜率与二阶段等时线相同（图3-12）；据斜率r₁可求出t₁或t₂，详见模式Ⅲ－3，但无法求出t₃。

模式ⅩⅢ当（X_t2,Y_t2）呈线性分布于数条平行直线上（其斜率为r₁），若ε₁恒定，β_j对不同样品取值不尽相同，则由（3.1.12）、（3.1.13）式可导出：当样品点足够多时，样品点分布于一个菱形区域内，类似于图3-11所示的三阶段铅样品点的分布状态；据两组直线斜率r₁、r₂能求出t₁、t₂或t₃,r₂表达式为：

当

非常数时

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

当

为常数时

燕山陆内造山带金-多金属成矿作用与构造-成矿关系

以上从理论上分析了开放体系多种情况下铅同位素的演化模式。可以看出，混合铅样品点呈同一或类似分布状态时，可对应一个至数个不同的地质过程。因此在应用铅同位素研究地质问题时，应尽量取足够多的样品；在样品点足够多的前提下，结合其它地质与地球化学资料进行综合分析，以便合理地解释铅同位素的演化，求出成岩、成矿时代。这些模式在燕山地区成岩成矿时期研究中，取得了较好的应用效果。

图3-12模式Ⅻ图解

Fig.3-12Lead-isotope evolution of modelⅫ

二、其它研究方法简介

1.据矿石蚀变矿物的K-Ar法、Rb-Sr等时线法、裂变径迹法确定成矿时代

上一章已述，燕山地区大部分类型的矿化都伴有强烈的蚀变，蚀变阶段性与矿化阶段性存在良好对应关系，两者形成时间相近。因此，蚀变矿物的同位素年龄能代表成矿时代。

蚀变矿物绢云母、白云母、钾长石等适合于K-Ar法年龄测定，白云母、绢云母的K-Ar法年龄能较好地反映同期矿化时代。

近矿蚀变矿物绢云母、白云母等的单矿物Rb-Sr等时线法年龄也能准确地反映成矿时代，是确定矿床形成时代的良好方法。

蚀变矿物的裂变径迹法年龄常较实际成矿时代偏小，其上限能大致代表成矿时间（杨应平，1985，硕士论文）。

2.据赋矿围岩时代与矿区岩脉时代间接推断成矿时代

当有充分的资料说明矿化与围岩成岩作用存在成因联系时，围岩时代能代表成矿时代下限。表3-1说明燕山地区中生代赋矿岩体时代与矿化时代的一致性。

当矿区内存在大量岩脉时，根据岩脉时代及岩脉与矿体相互穿切关系，也能较好地推断成矿时代。

表3-1岩体与其中金矿时代对比表

3.据同成矿期控矿构造的成生、活动时间推断成矿时代

任何控矿构造都属于某一个或某些构造体系，皆有一定的形成与活动时期；因此据同成矿期控矿构造的时代能定性推断部分矿床的成矿时代。古构造筛分有助于这方面的研究工作。

8. 研究方法

利用注水井吸水剖面、小层沉积微相和数值模拟三种方法综合研究南区沙二下_1-5层系剩余油分布规律。

1.注水井吸水剖面法

注水井吸水剖面法是利用历年来注水井吸水剖面资料，将注水井累积注水量分配到小层，再根据室内岩心水驱油试验结果，注入体积倍数与采收率、含水率之间的关系，来确定小层剩余油分布规律。

（1）建立静态数据库，统计小层渗透率分布规律

系统建立南区沙二下_1-5层系油、水井静态参数数据库。利用算术平均法和有效厚度加权平均法，分别计算出各小层渗透率平均值。利用概率统计的方法，求出各小层渗透率分布变异系数。

（2）建立吸水剖面数据库，计算小层累积注水量

在静态数据的基础上，建立注水井吸水剖面数据库。利用吸水剖面数据库可以统计出历年单井、小层吸水厚度变化趋势和吸水强度分布规律。利用吸水剖面数据库和注水井单井累积注水量，可以计算出历年小层累积注水量。

（3）建立注入体积倍数与采收率、含水之间关系，计算小层采出程度

根据濮城油田南区濮检1井非稳定流油水相对渗透率、水驱油试验报告和沙二下第446号岩心试验结果，由小层累积注水量计算出小层注入体积倍数，再根据以上关系内插求出各小层的采出程度和含水率。

（4）确定小层驱油效率

根据利用中原油田开发室内试验数据统计出来的驱油效率E_D试验公式：

高含水油田剩余油分布研究：以辽河油田欢26断块为例

驱油效率E_D可以做为小层在均质条件下的最终值，驱油效率E_D1可以做为小层在非均质条件下油田开发的最终值，或称测算采收率。在油田开发中，驱油效率还受注采井网及工艺技术条件的限制。

（5）计算小层剩余油量

根据小层驱油效率计算出可采储量，再由小层采出程度计算出剩余油量。

2.小层沉积相法

通过对濮城油田沙二下段沉积相的研究，认为濮城沙二下段沉积环境为浅水湖泊相和浅水三角洲相，其特点是水下分支河道异常发育，水下河道亚相是沙二下段沉积主体和骨架，河道层序具有对称性，底部粗粒段和顶部细粒段较薄、中间段厚度大且粒度均匀，河道砂体是本区沙二下段主要储集层；南区沙二下长期处于水下河道沉积区，砂层多，分选好，是濮城油田沙二下中的最好储集层。

针对沙二下_1-5油层目前开发现状，结合沉积相研究和油水生产剖面的初步分析，得到以下认识：

（1）河道砂是主要的吸水层，也是目前的主要产出层

在油田开发初期，河道砂（包括水下河道主水流线上的SH型砂体，居非主水流线上的H型砂体和居水下河道中的相对高台上的T型砂体）是主要的吸水层，也是主要的产油层。到油田开发中后期，由于油田含水的升高，主产层逐步过渡到主产水层。

根据1987年至1991年注水井吸水状况分类统计，河道砂是注水井的主要吸水层，统计48口注水井的吸水剖面，河道砂的射孔厚度204.5m，占总射孔厚度的45.7%，河道砂的绝对吸水量2692.2m³/d，占总吸水量的66.3%。其中1988年至1990年，河道砂射孔厚度占总射孔厚度的53%左右，绝对吸水量的百分数却高达80%以上。1987年至1990年，在射开河道砂厚度相对稳定的情况下，注水井中河道砂体的吸水能力有增大趋势，相对吸水百分数由57%增大到90%。

根据9口生产井产出剖面统计资料（表4-14），河道砂也是目前主要的产出层。统计沙二下_1-5层系河道砂射孔厚度45.1m，占总射孔厚度的40.1%，河道砂产液量122.3m³/d，占总产液量的64.8%。

（2）河道砂在注水井和生产井之间已经形成地下水道，是主要的产水层

根据濮3-284井环空测井资料分析，射开16层，产出层5个，产出层占31.3%；射开厚度33.5m，产出厚度16.4m，产出厚度占49.0%。其中主要产水层3²小层，2层5.0m，日产油1.7m³，日产水19.7m³，含水92.1%。

濮3-284井的一线注水井是3-282井，由于濮3-28井处于河流的边滩部位，油层物性差，吸水状况差。根据历次吸水剖面资料解释，射开有效厚度1.4m，日吸水量只有5m³左右，分析结果一线注水井不是主要的来水方向。

濮3-278井是濮3-284井的二线注水井，根据吸水剖面资料分析，是其主要的来水方向。濮3-278井沙二下3²小层，射开吸水厚度3.2m，日吸水量66.3m³。根据沉积相分析，濮3-278井和濮3-284井的沙二下3²小层处于同一河道砂体，它们之间连通性好、渗透性好，在油田注水开发中已经形成了地下水道。

（3）前缘砂和滨湖砂是目前主要的产油层

前缘砂分布在水道的两侧，滨湖砂距河道砂较远。前缘砂属中渗透砂体，滨湖砂属于低渗透性砂体。

统计沙二下_1-5层系主要处于前缘砂和滨湖砂部位的21口生产井，1992年9月份日产油水平289t，井数占全层系开井数的34.4%，日产油水平占56.1%。21口生产井平均单井日产水平13.8t，平均含水37.0%。其中处于前缘砂亚相的濮3-41井，生产沙二下_3-5，射开5层13.4m，其中有效厚度3层7.6m,9月份平均日产油16t，含水61%，累积产油7.09×10⁴t。

统计沙二下3²和沙二下5²两个典型含油小层，前缘砂2.32km²，滨湖砂3.02km²，分别占两小层含油面积的30.1%和39.0%。前缘砂和滨湖砂在平面上分布面积比较大，由于油层物性差、渗透率低，目前水驱动用状况差，剩余油量比较大，是今后挖潜的主要方向。

综合以上分析，河道砂是主要的吸水部位，同时也是主要的产出部位，过去是主要的产油层，目前是主要的产水层。含水一般均在80%以上，局部含水达到90%以上。目前剩余油很少，已到水洗油的阶段。大庆的河流过渡相和河漫相部位（濮城的前缘相与滨湖相）是目前主要的剩余油聚集带，也是目前主要的产油层，因此下步调整挖潜的方向应为河床过渡相和河漫相。

3.数值模拟法

（1）建立模型

①网格的划分

该模拟区块共有25小层，模型建立纵向上以主力层单独模拟层为原则划分为13个模拟层；平面上选取不等间距的矩形网格系统。整个模型网格总数为13×18×13=7254，其中有效节点4873个，死节点为2381个。

②油藏参数的选取

油藏流体物性参数。

相对渗透率数据：由于没有本区块油藏的相对渗透率数据借用邻近区濮检1井的数据进行了修正。沙二下_1-5共选用七条相对渗透率曲线。

PVT数据：南区沙二下_1-5层系没有取得PVT数据，故借用与其相近的东区文35井的数据进行了处理修正。

网格节点参数：网格节点数据除网格步长外，其他地质参数均来自每口井的电测解释结果，在工作站上用插值法算得每个网格的数据。

初始化计算结果：濮53块沙二下_1-5油藏由于未对每一小层储量进行标定，利用每小层体积百分数来计算每一小层储量。利用三维三相模拟各小层储量结果。

（2）历史拟合

根据生产历史对单井，全油田的压力、含水进行了拟合，均得到了较满意的结果。

9. 研究方法及技术思路

1) 充分收集和利用前人已有研究成果和资料，避免重复。周口店岩体和北京西山中生代区域热变质岩研究程度较高，已发表大量分析测试和实验数据，这是本书充分利用的研究起点。

2) 周口店岩体形成深度和源区深度的研究，采用: ①详细、精确测定周口店岩体各期次的化学成分、造岩矿物化学成分。特别注意分析数据的配套，造岩矿物化学成分主要用电子探针测定钾长石、斜长石、角闪石、黑云母等矿物的成分，强调矿物环带研究。②详细、精确测定周口店岩体造岩矿物和所含包体的化学成分。③源区地球化学示踪是本书的主要研究手段，中生代地壳加厚如果是由岩浆底侵作用带来的新生地壳物质加入造成的，那么物质在地壳中居留的时间应该较短，Nd 同位素的模式年龄应较年轻，且 ε_Nd值较高，而初始 Sr 同位素比值较小。

3) 区域和接触热变质作用深度、变质变形事件时代研究，采用: ①综合构造地层剖面法。在工作区选择地层发育齐全，变质带完整的地质剖面，从岩石的岩相学、构造变形等入手，用构造解析和地层分析的方法建立区域热变质作用发生的综合构造地层剖面。②矿物平衡的热力学分析法。在工作区选择变质程度最高的变质带，从岩石的岩相学观察入手，确定变质矿物的形成世代、期次以及矿物相平衡关系，用热力学计算分析方法计算矿物相平衡的 P － T 条件。根据北京西山地区区域热变质岩的岩石类型，矿物共生组合特征，重点对变泥质岩的石榴子石、黑云母、十字石、白云母、硬绿泥石的化学成分进行电子探针测定。以此为基础，选择矿物对、矿物组合地质温压计，求出区域热变质岩形成的温度、压力条件。③选择北京西山中窖剖面热动力变质岩中的白云母为研究对象，对其进行⁴⁰Ar /³⁹Ar 定年，限定该期变质变形事件 ( 岩体侵入作用) 发生的时间，进而探讨华北中生代构造体制转折的关键时限和岩石圈减薄的可能时间。

4) 南大岭玄武岩形成深度和源区深度的研究，主要是通过对前人资料的再利用，利用新的理论约束岩石圈地幔和下地壳。

5) 在地质事实和理论依托基础上，建立华北东部岩石圈减薄模式。