基于ADIFDTD模拟的岩石核磁共振横向弛豫特征分析.pdfVIP

第 3O卷 第 2期 计 算 物 理 Vo1．30．N0．2 2013年 3月 CHINESE JOURNAL OF COMPUTATIONALPHYSICS Mar．，2013 文章编号 ：1001-246X(2013)02-0237-07 基于 ADI-FDTD模拟的岩石核磁共振 横 向弛豫特征分析 葛新 民 ， 范宜仁 ， 曾青冬 ， 王明方 ， 徐拥军 (1．中国石油大学地球科学与技术学院，山东 青岛 266580；2．中国石油大学 CNPC测井重点实验室 ，山东 青岛 266580； 3．中国石油大学石油工程学院，山东 青岛 266580) 摘 要 ：将复杂 的骨架 一孑L隙系统抽象成等效双孔介质 ，根据 Bloch方程构建数学模型，用交替隐式时域有限差分 (ADI—FDTD)和联合反演迭代法 (SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振 T谱的反演，定量研究扩散 系数 、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响．结果表明：横向宏观磁化矢量衰减速率与扩散系数和 微孔隙分量成正比，与孑L隙宽度成反比，与表面弛豫速率基本无关．当扩散系数较大、孔隙宽度较小时，核磁共振 T谱难 以直观反映孔隙组分及孔隙结构．应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响． 关键词 ：核磁共振 ；ADI—FDTD；扩散系数 ；弛豫速率 ；孔隙组分 ；孔隙尺寸 中图分类号 ：P631．84 文献标识码 ：A 0 引言 核磁共振是研究孑L隙结构的重要方法之一 ，通过测量氢核的横 向宏观磁化矢量并进行反演得到核磁 T： 谱 ，进而获得与骨架无关的孑L隙度、渗透率 、孔隙分布情况等 ．横 向弛豫 由体积弛豫、扩散弛豫和表面弛 豫组成 ，其中表面弛豫所 占权重最大 ，人们通常忽略体积弛豫和扩散弛豫对总信号的贡献．由于横向核磁共 振未考虑孔隙问的扩散耦合效应 ，在复杂储层 (致密砂岩 、碳酸盐岩)中应用效果较差，核磁共振 T谱与毛管 压力、铸体薄片、扫描 电镜等实验得到 的孔 隙分布不 同，给孔 隙结构研究带来 了巨大的挑战 ．本文从核 磁共振原理出发，构建等效双孔介质模型，采用交替隐式有限差分 (ADI．FDTD)对 Bloch方程进行数值模拟， 分析扩散系数、表面弛豫速率及孔隙组分 比等对横 向宏观磁化矢量及 (与)核磁 T’谱的影响． 1 双孔介质等效模型及核磁共振响应 沉积 、成岩及后生改造作用使得储层 的原生孔隙受到破坏，导致孑L隙类型、孔隙大小及分布的复杂性． 对于复杂多孔介质岩石而言，其储集空间可看成是两组 (或多组)不同性质孔隙的组合 ，为了分析方便 ，将复 杂多孔介质等效成双孔介质．图1是真实的孑L隙结构到抽象的等效双孔介质示意图，图1(a)为某岩样的背 散射扫描电镜 (BSE)(暗色部分为孔隙，可看 出不 同大小的孔 隙随机分布在骨架中)，图 1(b)为简化后的等 效双孑L介质 (大圆部分代表大孑L隙，小圆部分代表微孑L隙)． 为了计算方便 ，将双孔介质模型等效成正立方体 ，其 —l，平面如 图2所示．设总孑L隙长度为 ，大孑L隙 组分的孑L隙度为 (b… 微孔隙的孔隙度为 咖 总孑L隙度为 咖 ．将微孑L隙体积与岩石总孔 隙体积 的比值 定义为微孔 隙组分 c，其公式为 i 。 一 ‘