CO2驱正演模拟与研究.pdfVIP

中国地球科学联合学术年会 2014 ·1093 · CO 驱正演模拟研究 2 1* 1 2 2 2 李 军 张军华 谭明友 崔世凌 曲志鹏 1.中国石油大学（华东）地学院 青岛 266580; 2.胜利油田物探研究院 东营 257022 1.引言 CO 驱油作为提高采收率（EOR ）的“三次采油” 的关键技术，既能够捕获CO ，减少温室气体排放， 2 2 又能够提高原油采收率，具有应用范围广，成本较低等特点[1] 。在当今资源短缺和节能减排的背景下， CO2 驱油技术被认为有着非常广阔的前景，得到世界各国的重视，但目前关于 CO2 驱替地震响应研究的 文献还很少。本文研究了 CO2 注入的浓度变化、不同方式及不同扩散形态引起的地震响应、速度、密 度等的变化，得到了一些有价值的认识，将对下一步 CO2 驱油地震监测评价提供理论指导。 2.正演方法介绍 本文分别采用了褶积与基于波动方程的正演模拟方法，来进行不同浓度的 CO2 驱替正演模拟。褶 积方法原理简单明确，但对于研究简单的水平层状模型的振幅变化还是很有效的。无噪声时，褶积方法 的具体表达式为： s (t) w(t) * r (t) (1) 其中s (t) 表示地震记录，w(t) 代表地震子波，r (t) 代表界面的反射系数。地震子波 w(t) 可以从井旁道 获取，反射系数r (t) 可以从声波时差和密度测井转换得到，由此可以模拟出地震记录。 波动方程的正演模拟时，我们使用的是交错网格有限差分方法，该方法是从非均匀各向同性介质中 的波动方程出发，具体采用了应用比较广泛的带PML 吸收边界条件的交错网格高阶有限差分格式，能 够有效地压制数值频散，具有较高的模拟精度和稳定性。在该方法中，震源函数我们选择的是雷克子波， 而稳定性分析方法采用的是傅氏分析方法。利用该方法，可以得到非均匀各向同性介质三维声波方程交 [2] 错网格高阶有限差分格式的稳定性条件为 ： 1 1 1 1 (2) Δ tV + + ≤ 2 2 2 N Δ x Δ y Δ z ∑ c l l 1 3.模型测试及分析 模拟胜利油田高 89-1 二氧化碳驱试验区块，设计了一个 7 层的水平层状模型，其中第四层为目标 层，CO2 从层中间注入，随着 CO2 向外扩散，其饱和度则随之减小，即 CO2 浓度从中间向两侧递减。 分别向目的层注入 20%、40% 、60%、80%及 100%浓度的 CO2 [3] ，然后根据具体的岩石物理参数，由 Gassmann 方程计算得到目的层不同 CO2 浓度下其余各层的速度、密度参数。而后分别采用褶积正演模 拟方法及交错网格有限差分正演模拟方法对模型进行正演，并对得到的正演模拟记录进行了分析，得知 注入的 CO2 浓度