基于频散测试分析的时频电磁资料校正处理及储层评价

沈义斌，杨 俊，曹 阳，刘雪军，王财富，何展翔

1.深圳市深远海油气勘探技术重点实验室（南方科技大学），广东 深圳 518055；2.广东省地球物理高精度成像技术重点实验室（南方科技大学），广东 深圳 518055；3.南方海洋科学与工程广东省实验室，广东 广州 511458；4.东方地球物理公司，河北 涿州 072751；5.南方科技大学 地球与空间科学系，广东 深圳 518055

0 前言

电磁法是油气勘探的主要方法之一，经过多年发展取得了长足进步.随着勘探目标的不断加深，以及勘探地区的地质构造越来越复杂，常规的电磁勘探方法已不能满足勘探开发的要求，为解决常规方法勘探能力的不足，人们开始研发其他有效的电磁勘探方法.

时频电磁法（Time Frequency Electromagnetic Method，TFEM）勘探技术就是在油气勘探走向深入，难度越来越大的背景下，为适应油气勘探开发的需要，经过20 余年的探索发展起来的适用于深部资源探测的新方法［1］.TFEM 作为一项比较成熟的油气检测技术，在国内外多个盆地已得到了广泛的应用，并取得了较好的勘探效果［2-6］.经过多年的实践表明，该技术是在目标储层岩石物性研究基础上，有效地提取储层油气相关异常信息，弥补了地震在油气检测方面的不足.

岩石的导电性是电磁勘探的基础.储层岩石电阻率参数一直是区分油水界面、评价储层含油性与程度的核心参数.相较于岩石物理弹性理论研究方面，岩石物理电性研究起步稍晚［7-10］.通过岩石物性研究可充分了解储层含油气与岩石物性参数之间的相关性，为电磁油气检测提供必要基础［11-12］.但是，时频电磁油气检测一直停留在定性评价的基础上，直到2021 年赵云生等利用Archie 公式对井地电磁法的电阻率、极化率资料进行火成岩油气储层饱和度计算才突破了电磁方法定量储层评价的方法瓶颈［13］.另外，曲昕馨等基于温压与电阻率关系对电磁反演的电阻率数据进校校正研究，建立了不同类型岩石的电阻率与温度、压力定量关系，提出了一种基于温压与电阻率关系的电磁反演校正方法［14］，通过对实测数据的校正处理，证明了该方法可有效提高对薄层的精细划分能力，但并没有进一步开展油气储层的评价.

本文利用西部CXN 地区火山岩发育区时频电磁反演数据，在岩石物理实验测试数据分析的基础上，研究储层岩石电性参数随温度和压力变化的特征及校正方法，进而进行储层饱和度预测的探索.

1 岩石温压频散测试分析

岩石复电阻率频散实验的岩心均取自CXN 地区.岩心有砂岩、灰岩、白云岩3 类（表1），合计25 块岩石样本.本研究用于岩石频散实验的实验仪器为AutoLab-1000，测量的频率范围0.01～10 000 Hz，一共31 个频点.岩心压力的频散测试时，温度控制为25 ℃，测试10、20、30、40、50、60、70 MPa 下岩心电阻率频散曲线；温度的岩心频散实验，设置围压10 MPa，温度为30、45、60、75、90、100 ℃，测量岩心的电阻率频散曲线.

表1 实验岩心参数表Table 1 Parameters of cores for experiment

1.1 压力与电性参数之间的关系

储层的主要岩性多是含有孔隙的岩石，随着地层深度的增加，储层对应的地层压力增大，岩石受应力作用更加致密.储层的主要岩性是少含导电矿物的岩石，地层水会作为岩石孔隙内部的主要成分，所以岩石内部的导电机制中，电磁感应所占比重将比较小，主要为离子导电.本实验所选用的岩心包括灰色白云质灰岩、灰色灰岩、灰色含砂质硅质灰岩3 种岩石.不同压力下的电性参数频谱特征均大致相同，这里仅以137号灰色含砂质硅质灰岩为例，从10 MPa 开始，70 MPa 结束，压力递增梯度为10 MPa（图1）.

压力对电阻率数值的影响显而易见.从图1 中不难看出，在10 MPa 压力下岩心电阻率的幅值处于1 200 Ωm左右；但在70 MPa 压力下时，电阻率数值接近2 900 Ωm.同样，相位也会随压力变化逐渐变化，但整体的数值差异不大.曲线首端接近重合，整体呈现下降的趋势，且施加的压力越大，相位数值越小，对于岩石的极化模式影响很大.

图1 不同压力条件下灰色含砂质硅质灰岩（137 号岩心）的频谱特征Fig.1 Spectrum characteristics of sandy siliceous limestone（core 137）under different pressure conditions

1.2 温度与电性参数之间的关系

当外界温度增大的时候，岩石内部的导电离子活动性能也会随之变强，因此岩石内部的导电能力也会变强，而岩石的电阻率则会随之减小.图2 为114号岩心在不同温度情况下复电阻率随温度变化的图像.

图2 不同温度条件下灰色白云质灰岩（114 号岩心）的频谱特征Fig.2 Spectrum characteristics of dolomitic limestone（core 114）at different temperatures

从图2a 可见，在不同温度条件下，灰色白云质灰岩复电阻率振幅曲线首端有明显差距，而且随着温度的增大，曲线趋势没有发生太大改变，复电阻率振幅的减幅越来越小.曲线随着频率的增大逐渐靠拢，并一直保持下降的趋势，直到进入高频段以后，曲线开始重合并最终相交.而相位曲线的情况与复电阻率振幅曲线相反，呈现稳定的下降趋势（图2b）.可以看到在中低频段，曲线能完全重合在一起；随着频率增大，曲线开始下降并逐渐分开，且在同一频率时，温度越低，相位越小.

值得注意的是，电阻率的幅值差异较大，甚至能达到几倍的差距.30 ℃时超过了120 000 Ωm，而100 ℃的灰色白云质灰岩复电阻率降低到40 000 Ωm.由此可见温度对于岩石复电阻率有很大的影响.这说明在电磁资料反演过程中，忽视温度对于电阻率的影响，可能会给地球物理数据处理带来明显的误差.

1.3 温度、压力与电性参数之间的关系

对5 种岩石进行温压实验.图3 所示为116 号岩石样品灰色白云质灰岩的复电阻率振幅和相位随温度的变化曲线.从图3 可见，电阻率随温度升高而降低，相位随温度升高而增加.图4 为116 样品在30 ℃下复电阻率振幅和相位随压力的变化曲线.观察图4 可见，电阻率随压力升高而升高，相位随压力升高而降低.因此，在野外工作实际处理资料过程中，通过上述关系对电性参数进行温度和压力校正后，再进行数据处理，则可以提高储层解释精度和效果.

图3 灰色白云质灰岩（116 号岩心）复电阻率温度拟合Fig.3 Temperature fitting of complex resistivity for dolomitic limestone（core 116）

图4 灰色白云质灰岩（116 号岩心）复电阻率压力拟合Fig.4 Pressure fitting of complex resistivity for dolomitic limestone（core 116）

2 电阻率反演数据校正方法

在对目标位置的电阻率数据进行校正之前，笔者查阅了CXN 地区的地区温压系统，明确了研究区的压力梯度和温度梯度.根据前人对CXN 拗陷20 多个含气构造的主要地层压力的分析资料［15］，对研究区的地层压力分布特征有了一定的认识.按照压力系数对CXN 地区地层压力进行了分类（表2）.

表2 研究区地层压力分类表Table 2 Classification of formation pressure in the study area

刘震等［16］提出将地层压力和地温看作一个整体系统，那么含油气盆地的地温-地层压力系统就是一个独立的封闭系统，并给出了地温-地层的压力方程［16］：

其中K 和L 均为于封闭系统内部气体体积有关的常量，T 和P 分别为地温和地层压力.从西部某盆地地温、地层压力特征图（图5）可得知，随着深度增加，在800～2 000 m 时，地层压力处于2～4 MPa 的低幅超压.通过测量得出该地区的地温梯度为1.93 ℃/hm，地压梯度为2 MPa/hm，平均地表温度为27 ℃.

图5 研究区温压梯度曲线（据文献［17］修改）Fig.5 Temperature and pressure gradient curve of the study area（Modified from Reference［17］）

2.1 岩石温压关系曲线的拟合

在研究区根据重磁电震和YT-1 井测井资料，确定了测线TFEM1902 的储层范围：深度5 500～6 300 m处，主要岩性包含灰岩、页岩和白云岩3 种岩石.故以岩心编号为120 号的灰色灰岩进行温压数据拟合，获取复电阻率振幅、相位和温度、压力对应的数值关系，其中激发频率取0.01 Hz.对选取的数据拟合后得到了不同的温度和压力条件下灰色白云质灰岩的拟合图像，如图6、7 所示.

图6 灰色灰岩（120 号岩心）复电阻率温度拟合Fig.6 Temperature fitting of complex resistivity for gray limestone（core 120）

由图6a、7a 可知，随着温度和压强的增大，复电阻率振幅分别减小和增大，所以在实际的反演过程中，随着地层深度的增加，温压的综合作用存在互相抵消的情况.为了使反演结果更加精确，利用岩石物理方法获取复电阻率振幅对应温度压力的数值关系，为后面的时频电磁数据进行温压校正提供基础.

2.2 典型剖面储层目标的圈定及数据校正

1）提取TFEM1902 测线资料，绘制时频电磁反演电阻率剖面（图8）；同样，绘制TFEM1902 测线的极化率剖面（图9）.

图9 测线EM1902 极化率剖面Fig.9 Polarizability profile of survey line EM1902

2）将对应测线的极化率剖面与时频电磁反演剖面结合起来比较，结合一级有利区某些深度对应的极化率较高的特点，可以在图8 中深度约3.6 km 处大致确定测线TFEM1902 的储层范围.

图8 测线EM1902 电阻率反演剖面Fig.8 Resistivity inversion profile of survey line EM1902

为了同时考虑到温度、压力对电阻率的影响.曲昕馨等［14］提出电阻率温压校正的综合函数：

图7 灰色灰岩（120 号岩心）复电阻率压力拟合Fig.7 Pressure fitting of complex resistivity for gray limestone（core 120）

其中C 为校正系数，可通过温度、压力梯度共同定义，可表示为：

式中，ρ（H）为地层深度H 下的电阻率定义，即常规勘探测得的电阻率数据；ρ（H0）为初始深度下的地层电阻率；a 为岩石对数电阻率的地温梯度；b 为岩石对数电阻率的地压梯度；ml为岩石的胶结因数，l=1，2，3……；Φ 为地层岩石的孔隙度，可由剖面的测井曲线给出.由图5 可知CXN 地区的地温梯度和压力梯度，其中，∂P/∂H=2，∂T/∂H=1.93.根据前人研究，有a=1.05，m=1.53，b=1.01，n=1.93［18］.根据以上数据，可以通过式（2）（3）获得校正后的地层电阻率ρ（H），如图10 所示.

图10 测线EM1902 校正后反演剖面Fig.10 Corrected inversion profile of survey line EM1902

3 储层含油饱和度的估算

饱和度对于储层的评价十分的重要.测饱和度模型有威克斯曼-史密斯（W-S）［19］等类型，为便于计算一般使用Archie 公式［20］：

其中，Rt为储层地层电阻率（本研究为考虑温压影响的电阻率）；Rw为储层地层水的电阻率；Sw储层的含水饱和度.根据探区测井资料，地层水电阻率为0.65 Ωm，储层平均孔隙度为0.17，m 取值为2，饱和度指数n 设置为2.

这里公式中储层地层电阻率Rt为考虑激发极化效应的复电阻率振幅，带入前面的储层电阻率和极化率数值，同时假定时间常数和频率相关系数分别为10和0.5，利用Cole-Cole 模型计算近似的等效复电阻率.

利用公式（4）可对储层目标进行饱和度计算，选取测线有利储层对应25 km（YT-1 井）附近，纵向深度为1.6～3.9 km 处的数据绘制得到随地层深度变化，电阻率和饱和度变化的曲线，如图11 所示.可见，估算的储层饱和度达到80%以上，这与研究区火成岩储层工业油井的饱和度一致.

图11 测线TFEM1902（1.6～3.9 km）储层饱和度估算Fig.11 Estimation of reservoir saturation of survey line TFEM1902（1.6-3.9 km）

4 结论

通过岩石温压电阻率频散测试，获得研究区储层岩石随着温度和压强的增大，复电阻率分别减小和增大.因此，随着地层深度的增加，温度和压力的变化对于地层物性影响是综合作用.同时，通过岩石饱和度电阻率极化率测试研究，获得饱和度随复电阻率变化关系；利用温压校正的综合函数关系式计算地层电阻率，进而基于研究区相关储层特性开展饱和度评价.该研究无疑还存在诸多不足，但不失为时频电磁从储层含油气定性预测走向储层定量评价进行的一种探索.

致谢：本项研究受到国家自然科学基金项目“深地/深海探测中强电流激发下可控源电磁法激电效应机理研究及应用”（41874085）、深圳市深远海油气勘探技术重点实验室项目（ZDSYS20190902093007855）、南方海洋科学与工程广东省实验室（广州）人才团队引进重大专项（GML2019ZD0203）、广东省海洋经济发展（海洋六大产业）专项资金项目“海洋电子信息装备浅海试验场”（GDNRC［2021］60）、广东省地球物理高精度成像技术重点实验室项目（2022B1212010002）、深圳市科技计划项目“深海深地资源探测技术系统研发”（KQTD 20170810111725321）的联合资助，在此表示感谢.