19期通排P13开始



Similar documents
1 79 Key words Tarim Basin Bachu -Maigati area Silurian growth fault sandy dispersion system Fig. 1 1 Overview of Bachu -Maigaiti regio

( 1) ; 6 5 [ 3 ] S1 S2 1 Fig. 1 Seismic sequence prof ile framework of Wuliyasitai Sag S1 () ( 2) S2 ( ) Ο 2 S1 Fig. 2 Sedimentary fa


1 引言

4 C C 1 2 S C C > m 6 C C Fig. 1 Geological sketch map of block-c in Yitong Basin

<4D F736F F D20C8EDCDC1B5D8BBF9CDB2BBF9CAD4B2C9C6BDCCA8B5C4CACAD3C3D0D4B7D6CEF6>

14 17.,., Fig. 1 Bougeur gravity map of research region 2 Fig. 2 Magnetic map of research region

* CUSUM EWMA PCA TS79 A DOI /j. issn X Incipient Fault Detection in Papermaking Wa

km t 2 4 Fig Tectonic divisions of Siberian platform

km BIF km BIF 6 7 BIF MT 14. MT. Fig. 1 1 Geological sketch map of Hanling -P

ACTA SEDIMENTOLOGICA SINICA Vol. 31 No. 2 Apr

13-15 Lagrange 3. 1 h t + hu + hv = 0 1 x y hu + t x hu gh 2 ( ) + y huv = - gh z 0 ( + x u u 2 2 槡 + v + W C ) 2 x + fhv + z h x 2hv u ( t x )

cm /s c d 1 /40 1 /4 1 / / / /m /Hz /kn / kn m ~

and explored deeply. Key words South Sumatra Basin petroleum geology hydrocarbon distribution exploration potential Barisan 1 1

~ 6km km a 1 Fig. 1 Location and tectonic division of Songliao Basin 19 5b 5c 5d Fig. 4 Seismic profile of channel in th

<D2BDC1C6BDA1BFB5CDB6C8DAD7CAB8DFB7E5C2DBCCB3B2CEBBE1C3FBB5A52E786C7378>

doc

(1) 1 ; 5 ; ; km ~200 :YC( )13-1 YC13-4 YC13-6 m; YC7-4 YC14-1 YC19-1 YC21-1 BD( ) V2 BD13-1 BD13-3 BD13-3S BD

增 刊 谢 小 林, 等. 上 海 中 心 裙 房 深 大 基 坑 逆 作 开 挖 设 计 及 实 践 745 类 型, 水 位 埋 深 一 般 为 地 表 下.0~.7 m 场 地 地 表 以 下 27 m 处 分 布 7 层 砂 性 土, 为 第 一 承 压 含 水 层 ; 9 层 砂 性 土

准北春晖油田油气勘探快速突破的三点启示

3 期 李 宗 发 : 贵 州 构 造 岩 土 体 分 区 及 其 与 地 质 灾 害 形 成 的 关 系 33 较 齐 全, 中 晚 元 古 宙 以 海 相 陆 源 碎 屑 沉 积 为 主, 主 要 分 布 于 黔 东 梵 净 山 雷 公 山 地 区 及 黔 中 一 带, 尤 其 在 黔 东 地

Microsoft Word 定版

Microsoft Word - 刘 慧 板.doc

在 培 养 职 前 卓 越 化 学 教 师 的 院 校, 会 编 一 本 过 去 称 作 化 学 教 学 论 实 验, 现 在 拟 为 卓 越 化 学 教 师 教 育 实 验 教 学 研 究 的 教 材 各 院 校 对 这 门 课 程 所 给 的 学 时 不 太 一 样, 但 都 是 围 绕 实 验

GLOBAL GEOLOGY Vol. 35 No. 2 Jun TOC 4. 87%

Microsoft Word tb 赵宏宇s-高校教改纵横.doc

77

T R 1 t z v 4z 2 + x 2 t = 2 槡 v t z 200 m/s x v ~

untitled

Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candi

E T 0 = γ 0 = 1 + R γ = nσ n nσ n ΔT 2 i - Σ n ΔT i T Pi - Σ n Σ n

TOPIC 专 题 45 1 加 快 农 业 大 数 据 发 展 的 现 实 意 义 农 业 大 数 据 运 用 大 数 据 的 理 论 技 术 和 方 法, 解 决 农 业 领 域 数 据 的 采 集 存 储 计 算 和 应 用 等 一 系 列 问 题, 大 数 据 技 术 是 保 障 国 家 粮

第 8 期 张 桂 荣 等 : 基 于 WEBGIS 和 实 时 降 雨 信 息 的 区 域 地 质 灾 害 预 警 预 报 系 统 1313 常 见 地 质 灾 害 具 有 发 生 频 率 高 分 布 地 域 广 的 特 点 随 着 我 国 经 济 的 发 展 和 工 程 建 设 规 模 的 扩

85% NCEP CFS 10 CFS CFS BP BP BP ~ 15 d CFS BP r - 1 r CFS 2. 1 CFS 10% 50% 3 d CFS Cli

蒋 哲等 基于 CVM 的有机农业认知程度和消费意愿研究 201 (P 0.000). On the basis of these results, enhancing publicity, pricing reasonably, focusing on high- income groups ar

GLOBAL GEOLOGY Vol. 29 No. 3 Sept P531 A doi /j. issn

Microsoft Word - 11-秦华伟.doc

% km % Nazca 15 8 cm Fig. 1 Schemat

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2

11 25 stable state. These conclusions were basically consistent with the analysis results of the multi - stage landslide in loess area with the Monte

Microsoft Word 方刚_new_.doc

2013年高等医学教育临床教学研究高峰论坛

Fig. 1 Frame calculation model 1 mm Table 1 Joints displacement mm

1 GIS 95 Y = F y + (1 F) (1) 0 0 Y0 kg/hm 2 /day F y 0 y c kg/hm 2 /day [12] y m 20 kg/hm 2 /hour Y = cl cn ch G [ F( y ) T m yo + (2) (1 F)(

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流

GLOBAL GEOLOGY Vol. 35 No. 4 Dec

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315.

海相油气地质 M A R I N E P E T R O L E U M 柯坪塔格组储层岩矿成分各分区不尽相同 塔东地区主要是砂岩 塔中 塔北及柯坪地区砂岩 中含量有所增加 以砂岩为主 夹部分砂 岩 横向上由南 塔中 向北 塔北 含量逐渐增 大 成分成熟度变高 由塔东地区 英吉苏凹陷 满 东地区 塔

3 : 505.,,,,,,,,,, 21 [1,2 ] , 21,, 21,, : [3 ]. 1. 3,, 10, 2 ( ),,, ; ; 40, [4 ]. 46, : (1),, (2) 16,,,,, (3) 17, (4) 18,, (5) 19,, (6) 20

Mnq 1 1 m ANSYS BEAM44 E0 E18 E0' Y Z E18' X Y Z ANSYS C64K C70C70H C /t /t /t /mm /mm /mm C64K

2 345 on the deposition. The high-value areas of mudstone thickness were the main developmental areas for the mudstone where were mainly semi-deep lak

2

Microsoft Word - 19王建华.doc

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

华山山前断裂中段晚第四纪活动的地貌表现及响应

5 : 7 [11-12] ( ) [13-14] Diessel [12] TPI/GI [15-16] (K3) 1 K1 50% K2 K3 K7 K2 K3 (20%~30%) K2 K10 1 Fig.1 Lithotypes of major coal seams 2.2 1

km km mm km m /s hpa 500 hpa E N 41 N 37 N 121

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang

~ 121 E 40 ~ 41 N MT Fig. 1 1 Location map of geographical exploration

Microsoft Word - A doc

Microsoft Word - d27 彭玉柱.doc

O V E R S E A E X P L O R A T I O N 66

KJWJ01 Fig. 1 Geological map of the Kajiwa landslide PDH01 60 ~ 70m 80 ~ 90m 80m 2 930m 556m 641m 598m m m 3 2 Ⅰ m 2

10 中 草 药 Chinese Traditional and Herbal Drugs 第 43 卷 第 1 期 2012 年 1 月 生 药 打 粉 入 药 的 基 本 特 点, 借 鉴 材 料 学 粉 体 学 等 学 科 的 研 究 成 果, 在 中 药 传 统 制 药 理 念 的 启 发

标题

698 39,., [6].,,,, : 1) ; 2) ,, 14,, [7].,,,,, : 1) :,. 2) :,,, 3) :,,,., [8].,. 1.,,,, ,,,. : 1) :,, 2) :,, 200, s, ) :,.

科胜伟达说明书

Microsoft Word - A _ doc

132 包 装 工 程 2016 年 5 月 网 产 品 生 命 周 期 是 否 有 与 传 统 产 品 生 命 周 期 曲 线 相 关 的 类 似 趋 势 旨 在 抛 砖 引 玉, 引 起 大 家 对 相 关 问 题 的 重 视, 并 为 进 一 步 研 究 处 于 不 同 阶 段 的 互 联 网

构象-----高分子物理的基础问题

[1] Nielsen [2]. Richardson [3] Baldock [4] 0.22 mm 0.32 mm Richardson Zaki. [5-6] mm [7] 1 mm. [8] [9] 5 mm 50 mm [10] [11] [12] -- 40% 50%

Vol. 22 No. 4 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Aug GPS,,, : km, 2. 51, , ; ; ; ; DOI: 10.

4 115,,. : p { ( x ( t), y ( t) ) x R m, y R n, t = 1,2,, p} (1),, x ( t), y ( t),,: F : R m R n.,m, n, u.,, Sigmoid. :,f Sigmoid,f ( x) = ^y k ( t) =

ph ph ph Langmuir mg /g Al 2 O 3 ph 7. 0 ~ 9. 0 ph HCO - 3 CO 2-3 PO mg /L 5 p

cm hpa hpa 2 45 N hpa ~ 12 Fig. 1 The observed rainfall distribution of Shanxi

United Nations ~ ~ % 2010

by industrial structure evolution from 1952 to 2007 and its influence effect was first acceleration and then deceleration second the effects of indust

4 731 Table TOC > 2% Log characteristics of different types of organic rich source rocks in Well FS2 /m TOC /% GR /API AC / μs /ft CNL /% DEN /

Mixtions Pin Yin Homepage

[1] Liu Hongwei,2013, Study on Comprehensive Evaluation of Iron and Steel Enterprises Production System s Basic Capacities, International Asia Confere

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN

第 02 期 1 医 疗 信 息 现 状 20% EMR Electronic Medical Record HIS HIS [1-2] 张 肖 等 : 基 于 大 数 据 的 医 疗 健 康 创 新 应 用 2 大 数 据 环 境 下 医 疗 数 据 特 征 分 析 PC [3]

Fig. 1 1 a-a b-b a-a σ ma = MPa σ a = MPa σ 0a = MPa 0. 9 σ t =135 MPa b-b σ mb = MPa τ b = MPa σ 0b =

实 践 探 讨 高 丽 : 从 少 数 民 族 大 学 生 的 阅 读 需 求 看 民 族 院 校 图 书 馆 的 资 源 建 设 有 区 域 性 和 民 族 性 很 强 的 传 统 学 科 特 色 学 科 及 优 势 学 科, 因 此 图 书 馆 的 资 源 建 设 也 要 顺 应 这 一 特 性

Scoones World Bank DFID Sussex IDS UNDP CARE DFID DFID DFID 1997 IDS

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

标题

(1) ( 1965 ),, 1952 [9] 2.1 (2) 1 53 (E i ), 2 (P i ) (G E (G P, 31 (Q i ) 3, : G E (x,y)= (E i Q(x i, y i )) E i G P (x,y)=

<B3ACBDDD>

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

STEAM STEAM STEAM ( ) STEAM STEAM ( ) 1977 [13] [10] STEM STEM 2. [11] [14] ( )STEAM [15] [16] STEAM [12] ( ) STEAM STEAM [17] STEAM STEAM STEA

Microsoft Word - 1--齐继峰_new_.doc

a 1b / 2 t t ms λ /4 < λ /4 1 Fig. 1 Model of high-speed wedge Marangakis f 100 E f = A f A* f df A f A* f d

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz


标题

34 22 f t = f 0 w t + f r t f w θ t = F cos p - ω 0 t - φ 1 2 f r θ t = F cos p - ω 0 t - φ 2 3 p ω 0 F F φ 1 φ 2 t A B s Fig. 1

m m m ~ mm


198 天 津 体 育 学 院 学 报 谁 来 教, 怎 么 教 的 问 题 长 时 间 悬 而 未 决, 或 泛 泛 而 谈, 或 知 行 不 一, 以 至 于 使 中 小 学 武 术 教 育 积 重 难 返, 停 滞 不 前 因 此, 切 实 推 进 中 小 学 武 术 教 育 的 发 展 要

Transcription:

科 技 导 报 2015,33(19) 井 震 结 合 精 细 储 层 建 模 方 法 与 实 践 张 义 1 1,2, 尹 艳 树, 秦 志 勇 3 1, 周 文 1. 长 江 大 学 地 球 科 学 学 院, 武 汉 430100 2. 山 东 科 技 大 学 ; 山 东 省 沉 积 成 矿 作 用 与 沉 积 矿 产 重 点 实 验 室, 青 岛 266590 3. 大 港 油 田 滨 港 集 团 博 弘 石 油 化 工 有 限 公 司, 天 津 300280 摘 要 融 合 地 震 资 料 与 井 资 料 分 别 在 横 向 和 纵 向 上 高 密 度 采 样 的 特 点, 以 中 国 海 上 某 气 田 B 区 块 为 例, 通 过 井 震 联 合 标 定, 进 行 精 细 构 造 解 释, 用 准 确 的 速 度 模 型 搭 建 了 时 间 域 与 深 度 域 联 系 的 桥 梁 ; 利 用 地 震 解 释 的 断 层 层 位 构 建 储 层 骨 架 模 型, 应 用 地 质 统 计 学 协 变 量 方 法, 以 井 数 据 作 主 变 量, 地 震 反 演 实 现 统 计 得 到 的 岩 性 概 率 体 作 为 第 二 变 量 约 束 井 间 插 值 外 推 和 模 拟, 建 立 砂 体 展 布 模 型 ; 依 靠 砂 体 展 布 模 型 约 束 沉 积 相 建 模, 在 相 模 型 基 础 上 构 建 物 性 模 型 实 践 表 明, 井 震 结 合 储 层 建 模 研 究 断 层 更 精 确, 能 够 弥 补 井 间 信 息 不 足 的 缺 陷, 有 效 降 低 井 间 砂 体 预 测 的 不 确 定 性, 提 高 模 型 预 测 精 度 关 键 词 井 震 结 合 ; 井 震 联 合 标 定 ; 储 层 建 模 ; 构 造 格 架 ; 岩 性 概 率 体 中 图 分 类 号 TE122 文 献 标 志 码 A doi 10.3981/j.issn.1000-7857.2015.19.015 Development of well- to- seismic integration reservoir modeling and its application ZHANG Yi 1, YIN Yanshu 1,2, QIN Zhiyong 3, ZHOU Wen 1 1. School of Geosciences, Yangtze University, Wuhan 430100, China 2. Key Laboratory of Depositional Mineralization & Sedimentary Minerals of Shandong Province; Shandong University of Science and Technology, Qingdao 266590, China 3. Bohong Oil Petrochemical Co. Ltd., Bingang Group, Dagang Oilfield Company, PetroChina, Tianjin 300280, China Abstract The characteristics of high- density sampling of seismic data and well data in the horizontal and vertical directions are integrated to construct the reservoir model for the Section B of a gas field of China. Firstly, the well-to-seismic integration calibration is used to interpret the fine structure, and the accurate velocity model is adopted to set up a bridge between the time domain of seismic data and depth domain of well data. Secondly, the reservoir structural model is built by well correlation under the constraint of the fault and horizon from seismic interpretation. Thirdly, a co- simulation method is developed to forecast the sand distribution, in which well data is the master variable and the lithology probability from seismic inversion is the second variable. The facies model is built on the basis of sand distribution model. Finally, the model of petro-physical property is built under the constraint of the facies model. Appliction has indicated that reservoir modeling with well-to-seismic integration is more accurate in terms of fault and sand body prediction between wells, thereby the improvement of accuracy of the reservoir model. Keywords integration of wells and seismic data; well- seismic integration calibration technique; reservoir modeling; tectonic framework; probability model of lithologic 收 稿 日 期 :2015-05-07; 修 回 日 期 :2015-06-25 基 金 项 目 : 国 家 自 然 科 学 基 金 项 目 (41572081); 山 东 省 沉 积 成 矿 作 用 与 沉 积 矿 产 重 点 实 验 室 开 放 基 金 项 目 (DMSM201404); 长 江 大 学 青 年 教 师 基 金 项 目 作 者 简 介 : 张 义, 硕 士 研 究 生, 研 究 方 向 为 精 细 油 藏 描 述, 电 子 信 箱 :1551269385@qq.com; 尹 艳 树 ( 通 信 作 者 ), 教 授, 研 究 方 向 为 油 气 藏 开 发 地 质, 电 子 信 箱 :yys6587@126.com 引 用 格 式 : 张 义, 尹 艳 树, 秦 志 勇, 等. 井 震 结 合 精 细 储 层 建 模 方 法 与 实 践 [J]. 科 技 导 报, 2015, 33(19): 91-95. 91

海上油田钻井资料相对较少 地震资料品质相对较好且 2 覆盖面广 常规建模算法依靠测井数据插值随机性太强 难 [1] 速度模型的建立 时深转换是井震结合的关键环节 高精度的时深转换取 以对井间区域进行控制 不确定性影响很大 应用地震资 决于速度模型的准确性 不合理的时深转换造成深度与时间 料不仅可以在空间上对构造格架的建立进行调整验证 还可 的不匹配 从而使井震结合失去意义 以利用其平面上丰富的信息提高井间估计的确定性对储层 对测井曲线进行多井一致性处理 对于扩径严重的井 砂体展布进行预测及控制 近 10 年来 综合地震与井资料建 段 由于泥浆充填 其密度信息严重失真 利用对扩径不敏感 立储层三维地质模型技术得到快速发展与应用[2,3] 利用地震 的自然伽马 自然电位等深度曲线通过多元线性拟合方法对 解释建立模型构造格架 地震反演预测井间砂体的分布 井 扩径段曲线进行重构 重构后的曲线能够有效反映储层岩性 震结合相模型 储层物性模型的建立等技术的发展 也逐渐 特征 [5] 从井旁地震道提取目的层子波 本着波形 相位 能 形成和完善了井震结合的储层三维地质建模技术 量相近的原则反复调整 得到吻合较好的合成记录 利用研 中国海上某气田 B 区块所处构造为挤压背斜构造 整体 是一个 NNE 向展布的复合大背斜 具有东缓西陡的特征 目 前该气田已进入开发期 由于河道的摆动和迁移 研究区砂 体总体表现为横向变化快 纵向砂体叠置 自下而上砂层变 究区 12 口井合成记录确定的时深关系 产生每口井的层速度 和平均速度曲线 [6] 通过分析发现工区目的层平均速度与埋 深具有较好的线性相关性 图 2 速度模型建立时利用井分 层对其做约束 确保了速度模型的精确性 薄 泥质含量增加的特征 亟需进行精细储层地质建模研 究 为开发调整提供有力依据 以往对该区的研究较多 但 对于薄互层储层建模方面研究较少 本研究针对研究区储 层砂泥叠置交互严重等问题 采取井震结合储层地质建模技 术建立研究区精细地质模型 利用地震解释的断层和层位结 合井资料建立高精度构造模型 以多个地质统计学反演实现 统计分析得到的岩性概率体控制储层砂体的空间展布 1 建模思路 井震标定是连接地震 地质 测井工作的桥梁 其目的是 利用测井资料的高分辨率精细标定各岩性界面在地震剖面 图2 上的反射位置[4] 建立时间域与深度域的对应关系 在此基础 Fig. 2 上进行精细的构造解释 有效结合地震解释的层位 断层数 and time 据和丰富的井数据 建立高精度构造格架 利用统计学反演 建立波阻抗与岩相之间的关系 通过对反演得到的多个实现 进行统计分析 得到能够充分表征三维空间岩相发育概率的 平均速度 埋深 时间交汇 Intersection of average speed, buried depth 3 地震约束下构造格架模型的建立 岩性概率体 以岩性概率体作软约束 井数据作硬约束 建 构造格架模型建立是储层建模第一步 也是非常关键的 立砂体展布模型 所得模型充分融合了井震分别具有的高 一步 其准确与否直接决定地质模型是否真实可靠 构造格 纵 横向分辨率特点 利用砂体展布约束建立沉积相模型 最 架模型的建立 实质上是在地层精细划分和对比基础上 综合 后根据 相控 原则建立属性模型 其流程如图 1 所示 地质 地球物理等信息 对研究区储层层面 包括断层部位 在 空间上分布的一种预测 在本次构造格架模型建立中 由于 井资料较少 综合地震解释成果约束构造格架模型建立就较 为关键了 构造格架模型由断层模型和层面模型组成 3.1 断层模型 断层模型主要反映空间上断层展布特征 揭示地层错动 情况以及不同断块差异性 井震结合建立断层模型的步骤 1 地震解释成果时深转换 2 利用地震解释成果初步建立断 层模型 [7] 检查调整断层倾向及连接关系 3 断点归位 进行 三维空间校正 4 以深度域方差体 蚂蚁体作为空间约束和 检查方式 在三维空间内合理直观的对断层进行修正 [8] 根 井震结合储层建模流程 据上述步骤反复调整和验证 保证断层形态既符合实际的情 Flow chart of reservoir modeling constrained by 况又能够精细刻画其在三维空间的展布特征 较好地实现了 图1 Fig. 1 well-to-seismic integration 92 断层建模 图 3

关键问题之一 通过对研究区进行岩石物理研究及测井资料 综合评价可知 波阻抗能较好的区分砂泥岩 这为以地震数 据作为建模约束打下了良好的基础 本次研究采用两步建 模策略 1 通过序贯指示同位协同模拟方法以岩性概率体作 为协变量模拟砂岩展布特征 2 用模拟的砂岩展布模型约束 沉积相模型建立 即 相约束相 建模 4.1 网格转换 地震网格到角点网格的简单重采样会导致属性窜层 层 位信息不能准确保留 对后面约束造成误差 应用网格转换 图3 Fig. 3 B 区块断层模型 Fault model of Section B 工具 将时间域反演结果无缝地转换到深度域油藏建模所需 要的角点网格属性体 保证了层位的连续性 解决了对地震 及其属性进行网格转换时的关键问题 转化的过程中保持体 3.2 层面模型 模型采用角点网格建立 综合考虑 B 区块井的分布情况 构造走向及开发需要等特征 平面网格步长为 25 m 25 m 方 案 适当添加趋势线 最大可能提高平面网格质量 垂向上 为考虑层内夹层对油水分布影响 同时兼顾模型大小 按平 均 0.5 m 大小进行划分 积和几何特征一致 具体步骤 1 创建时间域构造格架模 型 2 创建 CPG 角点网格模型 3 将反演属性体导入时间域 构造格架模型 4 应用 Jason 软件的 Rescale 模块将属性体粗 化为深度域 CPG 角点网格模型 4.2 岩性概率体 地质统计学反演通过对井点处波阻抗数据与岩相类型 针对研究区井网稀疏及个别部位地层倾角较大 单纯利 进行统计分析 确定不同波阻抗值所对应不同岩相的百分 用井分层数据建立层面模型造成井点间误差较大 影响模型 比 以此作为局部先验概率 通过合适的插值方法进行内插 精度 采用地震解释层面作为约束 以井点数据为条件数 外推 得到多个岩性体实现 由于反演结果随机性很强 常 据 可有效解决这一问题 地震解释的相对构造幅度是准确 常需要对多个实现进行统计分析来消除误差 砂岩概率体正 的 采用井震结合构建层面模型 充分保留井点数据的准确 是其统计分析的结果之一 本研究通过对反演得到的多个 性和地震的相对趋势 以此建立了 3 个主要油组顶面的层面 岩性体在每个节点上砂岩出现的概率进行统计 得到的砂岩 方法 结合手工修改方式 完成研究区构造模型建立 图 4 最大概率 [9] 模型 在此基础上 以主要油组层面为趋势 采用厚度插值 概率体即为最终结果 它在每个节点上都代表了砂岩出现的 4.3 数据分析 变差函数反映了储层在三维空间的变化规律 是两点地 质统计学建模的核心 它包括变程 块金值 拱高 3 个主要特 征值 反映了储层参数的空间变化特征 利用变差函数提供 的全部结构信息 既可以用来分析和认识所研究的地质问 题 也可以从地质角度对变差函数进行一次检验[11] 变差函 数的变程大小 不仅能反映区域化变量在空间上具有相关性 的范围 同时还能从总体上反映出区域化变量的载体 如砂 体 在空间上的展布特征 从而可利用变程来预测砂体在三 维空间某方向的延伸尺度 以实现预测砂体规模的目的 [12] 图4 Fig. 4 B 区块构造格架模型 Tectonic framework model of Section B 变差函数拟合时 根据各沉积单元的特征 选取合适模型 反 复多次进行参数实验 拟合出最佳结果 研究区物源方向为北西向 稳定的三角洲沉积环境适合 实践表明 地震解释面的约束使井间构造幅度保留了地 采用球状模型 由于井网稀疏 仅依靠井数据内插外推难以 震解释的相对趋势 有效解决了数学插值在研究区井网稀疏 预测无井控地区砂体发育情况 应用岩性概率体做约束变得 及个别部位地层倾角较大情况下造成的井间构造误差 更好 尤为重要 岩性概率体在建模过程中的约束作用主要体现在 地体现出井间微小构造的变化[10] 对储层的空间形态和砂体边界的不确定性约束方面 序贯指示同位协同模拟是在序贯指示模拟过程中增加 4 岩性概率体约束下的沉积相建模 了与主变量相关的协同变量 通过沿不同随机路径求取各节 地震资料纵向分辨率与井资料相比低很多 地震资料在 点的累积条件分布概率函数 进而从累积分布函数中提取模 储层建模中的应用是否具有可行性是井震结合储层建模的 拟值 利用序贯指示同位协同模拟方法进行储层横向预测 93

能够忠实于已知井数据 同时井间与外推区域的随机模拟又 [13] 受到协同变量 岩性概率体 的约束 具有较高的预测精度 变差函数分析时将岩性概率体作为协同变量输入 以井 上岩性数据作为主变量 岩性概率体只是代表砂 泥分布的 可能性 是岩相建模的软约束 当岩性概率体与井资料发生 冲突时 相概率服从井数据 这样建立的砂岩展布模型既能 确保在井点处与地下实际地质情况相符 又发挥了地震资料 横向信息丰富的特点 4.4 沉积相建模 分析可知 研究区发育浅湖泥 水下分流河道 席状砂 河口坝 4 种微相 由于沉积微相的分布受控于砂岩展布 沉积 相模型建立采用 相约束相 的方法 即用上一步生成的砂岩 展布模型控制沉积相分布范围 在砂岩分布边界内应用序贯 指示模拟方法生成水下分流河道 席状砂 河口坝 3 种微相 砂岩边界外直接用确定性方法赋值为浅湖泥 模拟产生 20 个实现 分析发现 20 个实现相似程度非常高 这说明采用井 震结合建模能够有效降低多解性 增加井间确定性信息 对 20 个实现进行对比分析 结合地质分析的结果 优选出实现 8 作为参数建模用的相模型 图 5 展示了实现 8 对应的各模型 对比[14] (a) 砂岩概率体 (b) 砂体展布模型 (d) 仅依靠井点插值的沉积相模型 (c) 砂体展布约束下的沉积相模型 图5 Fig. 5 B 区块 H8c2 小层各模型对比 Model contrast in H8c2 layer of Section B 图 5 展示了 B 区块 H8c2 小层砂岩概率体 砂体展布模 依据 序贯高斯建模是属性建模常用的方法 它通过对储层 型 砂体展布约束下的沉积相模型 单纯使用井插值所得的 物性进行正态变换 采用高斯函数对参数进行分析和预测 相模型 对比分析表明 井震结合储层建模结合了井震资料 储层参数的分布受到砂体成因单元的控制 只有砂体内才具 各自的优点 储层砂体的展布和尖灭更加符合实际的地质认 有有效的储层参数 不同相的储层参数分布规律不同[16] 采用 识 模型既符合井数据的地质统计学特征 又能反映出地震 相控建模方法 按沉积相模拟储层物性参数的空间展布 建 [15] 数据中观测到的大尺度结构和储层连续性 与仅依靠井数 据内插外推得到的模型相比 大大降低了井间和无井控区域 预测的随机性 5 属性模型建立 储层属性建模采用序贯高斯建模方法 分析知所建相模 型能够较好地反映沉积相空间变化 可作为相控参数建模的 94 立了研究区孔隙度 渗透率模型 采用该模型计算的储量与 地质分析结果较为吻合 误差在 10%以内 以此开展数值模 拟与生产实际能够较快地拟合上 验证了模型的准确性 6 结论 1 井震联合标定 对研究区进行了精细构造解释并得到 单井时深关系 在此基础上建立了速度模型 精确的速度模

科 技 导 报 2015,33(19) 型 是 联 系 时 间 域 与 深 度 域 的 桥 梁 2) 利 用 地 震 解 释 的 断 层 层 位 所 构 建 的 储 层 骨 架 模 型 更 为 精 确, 能 有 效 解 决 单 纯 使 用 数 据 内 插 造 成 的 井 间 构 造 误 差, 更 好 地 体 现 出 井 间 微 小 构 造 的 变 化 3) 采 用 两 步 建 模 策 略, 应 用 地 质 统 计 学 协 变 量 方 法, 以 井 数 据 作 主 变 量, 地 震 反 演 实 现 统 计 得 到 的 岩 性 概 率 体 作 为 第 二 变 量 约 束 井 间 插 值 外 推 和 模 拟, 建 立 砂 体 展 布 模 型, 利 用 相 约 束 相 建 模 方 法, 用 模 拟 的 砂 岩 展 布 模 型 约 束 沉 积 相 模 型 建 立, 由 此 建 立 的 沉 积 相 模 型 既 符 合 井 数 据 的 高 纵 向 分 辨 率 特 征, 又 真 实 地 反 映 储 层 横 向 展 布 的 连 通 性 特 征, 有 效 降 低 了 井 间 砂 体 预 测 的 不 确 定 性, 提 高 了 储 层 相 建 模 精 度 参 考 文 献 (References) [1] 崇 仁 杰, 刘 静. 以 地 震 反 演 资 料 为 基 础 的 相 控 储 层 建 模 方 法 在 BZ25-1 油 田 的 应 用 [J]. 中 国 海 上 油 气 : 地 质, 2003, 17(5): 17-21. Cong Renjie, Liu Jing. Application of a facies- controlled reservoir modeling method based on seismic inversion data in BZ25-1 Oilfield [J]. China Offshore Oil and Gas: Geology, 2003, 17(5): 17-21. [2] 陈 恭 洋, 胡 勇, 周 艳 丽, 等. 地 震 波 阻 抗 约 束 下 的 储 层 地 质 建 模 方 法 与 实 践 [J]. 地 学 前 缘, 2012, 19(2): 67-73. Chen Gongyang, Hu Yong, Zhou Yanli, et al. Methods and practice of reservoir geologic modeling constrained by seismic impedance[j]. Earth Science Frontiers, 2012, 19(2): 67-73. [3] 周 丽 清, 赵 丽 敏, 赵 国 梁, 等. 高 分 辨 率 地 震 约 束 相 建 模 [J]. 石 油 勘 探 与 开 发, 2002, 29(3): 56-58. Zhou Liqing, Zhao Limin, Zhao Guoliang, et al. Stochastic facies modeling constrained by high resolution seismic data[j]. Petroleum Exploration & Development, 2002, 29(3): 56-58. [4] 舒 克 栋, 张 才 智, 程 超, 等. 约 束 稀 疏 脉 冲 反 演 在 泌 阳 凹 陷 南 部 陡 坡 带 中 的 应 用 [J]. 新 疆 石 油 天 然 气, 2011, 7(1): 16-20. Shu Kedong, Zhang Caizhi, Cheng Chao, et al. The application of reservoir sparse spike inversion in South Slope of Biyang Depression[J]. Xinjiang Oil & Gas, 2011, 7(1): 16-20. [5] 邱 立 伟, 杨 安 平, 龙 胜 祥, 等. 孤 家 子 - 后 五 家 户 气 田 地 震 反 演 与 储 层 建 模 综 合 研 究 [J]. 勘 探 地 球 物 理 进 展, 2004, 27(1): 50-52. Qiu Liwei, Yang Anping, Long Shengxiang, et al. Geostatistical seismic inversion and reservoir modeling study in Gujiazi- Houwujiahu region [J]. Progress in Exploration Geophysics, 2004, 27(1): 50-52. [6] 李 艳 丽, 柳 成 志, 王 文 广, 等. 井 震 结 合 断 层 建 模 方 法 在 萨 北 北 二 西 的 应 用 [J]. 科 学 技 术 与 工 程, 2011, 11(18): 4182-4185. Li Yanli, Liu Chenzhi, Wang Wenguang, et al. Fault modeling method with well- to- seismic integration used in the North Ⅱ West tract of Saertu North Area[J]. Science Technology and Engineering, 2011, 11 (18): 4182-4185. [7] 曲 良 超, 卞 昌 蓉. 井 震 结 合 断 层 建 模 技 术 在 复 杂 断 块 中 的 应 用 [J]. 断 块 油 气 田, 2012, 19(4): 426-429. Qu Liangchao, Bian Changrong. Application of fault modeling combined with well data and seismic data in complicated fault block[j]. Fault- Block Oil & Gas Field, 2012, 19(4): 426-429. [8] 刘 雷, 金 利, 文 宇 飞. 开 发 地 震 构 造 建 模 技 术 的 研 究 方 法 [J]. 高 原 地 震, 2009, 21(3): 46-51. Liu Lei, Jin Li, Wen Yufei. Discussion on the ways of development of seismic structural modeling research[j]. Plateau Earthquake Research, 2009, 21(3): 46-51. [9] 白 雪. 井 震 结 合 构 造 建 模 技 术 及 在 BSX 区 块 开 发 中 的 应 用 [J]. 科 学 技 术 与 工 程, 2013, 13(11): 3064-3068. Bai Xue. Research on high efficiency price heparin sodium with enzymatic combined membrane preparation technology[j]. Science Technology and Engineering, 2013, 13(11): 3064-3068. [10] 傅 志 明. 高 台 子 油 田 井 震 结 合 构 造 建 模 技 术 研 究 [D]. 杭 州 : 浙 江 大 学, 2012. Fu Zhiming. The research on the structure modeling technology by integrating the well and seismic data in Gaotaizi Oilfield[D]. Hangzhou: Zhejiang University, 2012. [11] 刘 泽 容, 杜 庆 龙, 蔡 忠. 应 用 变 差 函 数 定 量 研 究 储 层 非 均 质 性 [J]. 地 质 论 评, 1993, 39(4): 297-301. Liu Zerong, Du Qinglong, Cai Zhong. The qualitative study of the reservoir heterogeneity with the variation function [J]. Geological Review, 1993, 39(4): 297-301. [12] 孙 磊, 孙 龙 德, 江 同 文, 等. 塔 里 木 盆 地 轮 南 油 田 TⅢ 油 组 底 水 块 状 油 藏 夹 层 分 布 规 律 研 究 [J]. 海 相 油 气 地 质, 2001, 6(4): 47-52. Sun Lei, Sun Longde, Jiang Tongwen, et al. Research on bottom water block reservoir interbed distribution in Tarim Basin Lunnan Oilfield T-Ⅲ oil group [J]. Marine Origin Petroleum Geology, 2001, 6(4): 47-52. [13] 张 淑 品, 于 兴 河. 同 位 协 同 随 机 建 模 方 法 在 储 层 预 测 中 的 应 用 [J]. 天 然 气 地 球 科 学, 2006, 17(3): 378-381. Zhang Shupin, Yu Xinghe. The application of collocated Co- Kriging stochastic modeling reservoir description [J]. Natural Gas Geoscience, 2006, 17(3):378-381. [14] 杨 勇, 聂 海 峰, 张 雅 玲, 等. 基 于 多 点 地 质 统 计 学 的 岩 性 气 藏 精 细 建 模 方 法 与 应 用 [J]. 断 块 油 气 田, 2013, 20(6): 723-726. Yang Yong, Nie Haifeng, Zhang Yaling, et al. Method and application of stochastic modeling for lithologic gas reservoir based on multiplepoint geostatistics[j]. Fault-Block Oil &Gas Field, 2013, 20(6): 723-726. [15] 刘 文 岭. 地 震 约 束 储 层 地 质 建 模 技 术 [J]. 石 油 学 报, 2008, 29(1): 64-68. Liu Wenling. Geological modeling technique for reservoir constrained by seismic data[j]. Acta Petrolei Sinica, 2008, 29(1): 64-68. [16] 段 天 向, 刘 晓 梅, 张 亚 军, 等. Petrel 建 模 中 的 几 点 认 识 [J]. 岩 性 油 气 藏, 2007, 19(2): 102-107. Duan Tianxiang, Liu Xiaomei, Zhang Yajun, et al. Discussion on geologic modeling with Petrel[J]. Lithologic Reservoirs, 2007, 19(2): 102-107. ( 责 任 编 辑 科 技 导 报 研 究 论 文 栏 目 征 稿 研 究 论 文 栏 目 专 门 发 表 自 然 科 学 工 程 技 术 领 域 具 有 创 新 性 的 研 究 论 文, 要 求 学 术 价 值 显 著 实 验 数 据 完 整 具 有 原 始 性 和 创 造 性, 同 时 应 重 点 突 出 文 字 精 炼 引 证 及 数 据 准 确 图 表 清 晰, 并 附 中 英 文 摘 要 以 及 作 者 姓 名 所 在 单 位 通 信 地 址 关 键 词 等 信 息 在 线 投 稿 : 陈 广 仁 ) 95