《盆地模拟 新》由会员分享，可在线阅读，更多相关《盆地模拟 新（80页珍藏版）》请在金锄头文库上搜索。
1、盆地模拟盆地模拟BASIN MODELINGBASIN MODELING（4）埋藏史及成熟度演化该井的门限深度为1961m；该井的核桃园组烃源岩现今的演化程度均较低，处于低 成熟至中等成熟阶段的初期。其中东9井核三段的镜质体反射率为0.620.93%，地温为 105139，核二段的镜质体反射率为0.470.62%，地温为81105。东9井模拟现今温度趋势和现今Ro趋势与实测值的吻合关系 东9井D9D9D9D9井成熟度演化井成熟度演化该井的门限深度为1948m；该井的核桃园组烃源岩现今的演化程度均较低，处于低成熟至中等成熟阶段的初期，其中核三段的镜质体反射率为0.710.73%（至完井深度处），
2、核二段的镜质体反射率为0.470.71%。南65井模拟现今温度趋势和现今Ro趋势与实测值的吻合关系 南65井N65N65N65N65井成熟度演化井成熟度演化该井的门限深度为1787m；该井的核桃园组烃源岩现今的演化程度均较低，处于低成熟至中等成熟阶段的初期，其中张15井核三段的镜质体反射率为0.730.76%（至完井深度处），核二段的镜质体反射率为0.390.73%。张15井模拟现今温度趋势和现今Ro趋势与实测值的吻合关系 张张1515井井 Zh15Zh15张张1515井成熟度演化井成熟度演化（二）油气生成与运移 D4东庄地区东庄地区L480L480剖面油气运移路径模拟结果剖面油气运移路径模拟
3、结果1、核三段的原油成熟度都只是中等成熟度，生烃时间较晚2、烃源岩层内出现油气的初次渗流，而部分原油生成后沿边 界断裂出现大量的运移，但并未运移至上部储层聚集。 L480L480H21H21H22H22 H23H23Tg-T51Tg-T51H31H31北马庄地区北马庄地区L780L780剖面的油气运移聚集模拟剖面的油气运移聚集模拟油层油层1、该剖面上烃源岩成熟度较低，没有出现的大规模油气运 移，只是在边界的扇三角洲储层聚集了少量的油层。2、油气运移目前以就近运移为主，而连通砂体主要是横向 远距离的良好通道。L780L780H21H21 H22H22H23H23Tg-T51Tg-T51H31H3
4、11、边界断裂附近的油气沿断层纵向运移至高渗透储层（同 时是输导层）后在横向运移。2、该剖面中间部分的油气运移主要以横向运移的，而边界 断裂边缘则是以纵向运移为主。黑龙庙地区黑龙庙地区L1080L1080剖面油气运移路径模拟结果剖面油气运移路径模拟结果H21H21 H22H22H23H23Tg-T51Tg-T51H31H31L1080L1080黑龙庙地区黑龙庙地区L1080L1080剖面的油气运移聚集模拟剖面的油气运移聚集模拟 储层聚集的主要为油，含少量的气H21H21 H22H22H23H23Tg-T51Tg-T51H31H31 L1080L1080 油层油层核三1段烃源岩顶面在核二1段末
5、期成熟度分布图（30Ma）核三1段烃源岩顶面在廖庄组末期 成熟度分布图（22.5Ma） 核三1段烃源岩顶面在现今成熟度 分布图（0Ma）核三1段烃源岩顶面成 熟度不同时期演化图核二3段烃源岩顶面在廖庄组末期 成熟度分布图（22.5Ma）核二3段烃源岩顶面在现今成熟度 分布图（0Ma）核二3段烃源岩顶面成 熟度不同时期演化图盆地模拟盆地模拟一、盆地模拟概述二、盆地模拟理论与方法三、典型实例分析一、盆地模拟概述地质模型数学模型成果解释盆地模拟内容5“史”模拟地史模型热史模型生烃史模型排烃史模型运移聚集史模型盆地模拟发展简史1978年原西德尤利希核能研究有限 公司石油与有机地球化学研究所建 立了一维
6、盆地模拟系统，主要是通 过对欠压实地带古地温史和埋藏史 的计算，在此基础上求有机质成熟 度，进而对生烃量和排烃量计算。盆地模拟发展简史1984年法国石油研究院建立了一个 较完整的二维盆地模拟系统，利用 回剥法求埋藏史，利用地球动力学 法求古地温史，进而求烃类成熟史 ，通过两相运移法求流体压力史和 油聚集史。TEMISPACK 盆地模拟发展简史1984年南卡罗拉那大学地质科学系 提出了用镜质体反射率确定古热流 方法；1988年又提出了用其它几种 地化资料确定古热流的方法；90年 代初，又提出盆地流体四维监控模 拟系统。GEOPET 盆地模拟发展简史1987年英国不列颠石油公司（BP石 油公司）提
7、出了一个关于油气二次 运移聚集的二维模型，即在地下流 体运移相态中考虑到油相和气相。 该模型成为目前公开发表的论述油 气二次运移与圈闭问题的较好模型 。盆地模拟发展简史1981年日本石油勘探有限公司勘探 部建立了一个简化的二维盆地模拟 系统。1988年又与南卡大学合作完 成较为完善的二维盆地模拟系统， 该系统考虑流体运动、传热（传导 和热对流）、烃类生成和运移。盆地模拟发展简史1986年美国伊利诺斯大学地质系利 用有限差分法推出了二维盆地流体 模拟软件（BASIN2）。随后又与水 岩相互作用模拟结合起来，推出了 “GWB”（Geochemists Workbench ）模拟软件。BASIN2盆
8、地模拟发展简史我国盆地模拟研究起步较晚，在80 年代末开始在国外盆地模拟软件的 基础上进行修改完善。主要有中石 油勘探规划研究院石广仁等、中海 油研究中心与加州大学伯克利分校 地质地球物理系合作完成的盆地模 拟软件。BAS2（石广仁等）PRO-BAS二、盆地模拟基本原理（一）、地史模拟（二）、热史模拟（三）、生烃史模拟（四）、排烃史模拟（五）、运移聚集史模拟（一）、地史模拟地史模拟的功能是重建沉积盆地的沉积埋藏史和构造 沉降史。沉积埋藏史是指盆地中的某一沉积单元或一系列单元 自沉积开始至今或某一地质时期的厚度和埋藏深度 的变化历史，通常由沉积埋藏史曲线表示。构造沉降史是指盆地基底的某一单元或一
9、系列单元在 盆地发育时期上下运动的历史，通常由沉降史曲线 表示。沉积埋藏史模拟沉积埋藏史通常是由以深度为纵坐标、时 间为横坐标所表示的某一沉积单元的埋 深时间曲线（即沉积埋藏史曲线）给 出。沉积埋藏史的分析方法主要是地层 孔隙度厚度恢复法，即回剥法。沉积埋藏史模拟沉积埋藏史模拟1.压实校正的基本原理压实校正最常用的方法是回剥法。回剥法所使用的关键参数是 孔隙度深度曲线，其基本原理就是沉积压实原理，即：沉 积物在压实过程中岩石骨架体积保持不变，只是因岩石孔隙 度的变小而导致岩石体积的变小。回剥法的具体做法就是从已知的单井分层参数出发，按地质年 代逐层剥去，其间考虑了沉积压实、超压、剥蚀以及沉积间
10、 断等地质事件，直至全部地层剥完为止，最终得到的结果是 该井各地层的埋深地质年代关系历史。沉积埋藏史模拟设hs为沉积物的骨架厚度，h为沉积物厚度，则：hsh（1）在一般情况下，同一沉积盆地中岩性相同的沉积物， 其孔隙度与深度Z有如下函数关 系：i=ioexp（aiZ） 其中io为i种岩性沉积物的原始孔隙度，ai为i种岩性的压缩系数。各种岩性的Z关系式通常可由孔隙度测井、地震层 速度或实验室分析得到。 去压实校正沉积埋藏史模拟2、剥蚀厚度的恢复剥蚀厚度的恢复方法很多，概括起来主要有以下三种 ：即地质学方法、地球物理方法和成熟度法。1）区域地层对比法通过与邻区沉积层序完整的地层进行对比来确定剥蚀
11、时间和剥蚀厚度沉积埋藏史模拟2）剥蚀系数法RE 2h1/（h1h1+h2h2）HE （h1h）RE 其中RE代表剥蚀系数h1代表研究区最大沉积厚度h2代表研究区最小沉积厚度h1代表研究区最小残余厚度h2 代表研究区最大残余厚度HE代表剥蚀厚度h代表研究目的层的残余厚度。沉积埋藏史模拟3) 地球物理学方法主要是应用声波测井资料，根据在正常压实条件下，泥岩孔隙 度或声波时差与深度呈指数关系的数学模型来外推剥蚀厚度 。该方法最早是由Magara(1978)提出的。其前提条件是(1)剥蚀后新沉积的地层厚度的必须小于剥蚀厚度，否则原泥 岩孔隙度(或声波时差)将被改造而失去作用；(2)必须知道未压实前泥岩
12、沉积物的原始孔隙度0或原始地表 声波时差（t0）；(3)必须能够比较系统地获得纯泥岩的孔隙度或声波时差值。 沉积埋藏史模拟4）成熟度剖面法该方法是由华莱士(W.G.Dow,1977)提出来的。在正常 地质背景下，镜质体主要受最大温度和有效加热时 间的影响，镜煤反射率Ro的对数值与埋深之间呈良 好的线性关系。利用这一关系即可恢复地层剥蚀厚 度和计算平均古地温梯度。 据Dow的研究表明，地表附近Ro的最低值约为0.18 0.2，因此只要把剥蚀面以下的深度和Ro对数值的 回归直线向上延伸至Ro=0.2处，就可求出地层的剥 蚀厚度。沉积埋藏史模拟沉积埋藏史模拟3、异常超压影响异常压力型剖面(利101井
13、)构造沉降史模拟引起沉积盆地沉降的原因有两个，一个是盆地基底下 部活跃的地球动力作用，另一个是盆地内沉积物和 水的载荷。通常认为这两种过程是一种互相独立的 事件，所以盆地基底总的沉降就等于由地球动力作 用和沉积载荷这两者引起的下沉总和。麦肯齐在他所提出的拉张型盆地模型中预测了盆地的 构造沉降，他将构造沉降分成两部分。 构造沉降史模拟引起沉积盆地沉降的原因有两个，一个是盆地基底下 部活跃的地球动力作用，另一个是盆地内沉积物和 水的载荷。通常认为这两种过程是一种互相独立的 事件，所以盆地基底总的沉降就等于由地球动力作 用和沉积载荷这两者引起的下沉总和。麦肯齐在他所提出的拉张型盆地模型中预测了盆地的
14、 构造沉降，他将构造沉降分成两部分。 构造沉降史模拟均衡原理L*w + C* c + M* m = H* s + C* c + (M-Dl)* mDl = H* (s w)/(m - w)构造沉降史模拟1、初始沉降Si： 它是由于岩石圈在拉张前后的均衡补偿作用造成的最早沉降量， 相当于盆地初期大断层的初始垂直断距。其数学表达式为：Si = hl（oc）（hc/hl）（1Tlhc/hl）hlTlo/2（1 1/）/ o（1Tl）w hl为岩石圈厚度，125km；o为地幔密度，3.33g/cm3；w为海水密度，1.00g/cm3；c为陆壳密度，2.8g/cm3；hc为原始陆壳厚度；Tl为软流圈顶界
15、温度，1333；为岩石圈热膨胀系 数，3.281051；为地壳的始拉伸量(拉伸量参数)；正演模拟Mckenzic模型构造沉降史模拟2、热沉降St： 它是由于岩石圈冷却收缩引起的缓慢下沉的沉降量。当4时， 其数学表达式为：St = 4hloTl/2（ow）（/）sin（/） 1 exp（t/） 其中，hl为岩石圈厚度，125km；o为地幔密度，3.33g/cm3；w 为海水密度，1.00g/cm3；c为陆壳密度，2.8g/cm3；hc为原始 陆壳厚度；Tl为软流圈顶界温度，1333；为岩石圈热膨胀 系 数，3.281051；k为岩石圈的热扩散率，8.0410 7m2S1；是一个常数，=h l2/2k ，它的意义为瞬时热流影 响地温的最大限度时间，在麦肯齐模型中定为62.8Ma；t为盆 地沉降时 间，Ma；为地壳的始拉伸量(拉伸量参数)；为圆 周率。构造沉降史模拟初始沉降Si和热沉降St都是与重力均衡补偿及热膨胀有关的方程 ，它们的表达式中都有拉伸 量参数。该参数是一个不易确 定但又十分重要的参数，一般可通过实测的现今基底热流值 与模拟的热流现今值拟合的方法，利用计算机自动调整值 并且进行反复迭代而求取，但也可用沉积剖面计算构造热沉 降的方法求得，其公式为：StHs（
