白云岩声电各向异性实验测量及分析
时间：2023-06-05 22:25:07 来源：雅意学习网 本文已影响 人 
唐军，刘沁园，赖强，吴煜宇，许巍
(1.长江大学 地球物理与石油资源学院，湖北 武汉 430100；
2中国石油西南油气田分公司 勘探开发研究院，四川 成都 610041)
随着勘探技术的不断进步，中深层碳酸盐岩油气藏的开发为中国油气增储上产做出巨大贡献，其中白云岩储层较多，展现出了较好的开发前景[1-2]。但是中深层白云岩储层受到高温高压的影响，具有较强的各向异性。这不仅影响了水平井技术在白云岩储层中的应用效果，同时也大大降低了测井解释的正确率[3]。
声波、电阻率的差异是白云岩各向异性的直观反映，用声波和电阻率表征白云岩储层的各向异性效果较好[4]。经过广泛调研发现：专家和学者利用大量理论和物理实验研究证明了各向同性模型有较大的局限性后，各向异性模型成为了油气勘探的研究热点。Klein[5]指出，电各向异性主要出现在不同岩性的薄互层中或在电阻率差异较大的层状地层中；
徐亦鸣等[6]结合理论模型和实际地震资料，提取了纵波各向异性参数；
Moran等[7]提出假设,在均匀的无限厚各向异性地层中，电法测井仪器所测得电阻率值都会受到井斜和地层倾角的干扰；
伍文杰[8]利用声波测量系统对圆球形岩样进行声波测量和各向异性分析，发现该方法比传统测量方法效果好；
孙晶波[9]利用交互式速度分析方法，证明了当介质存在各项异性时，利用非双曲时差速度分析方法，能进一步改善声波波形的信噪比；
但是并没有表征声波时差各向异性的大小，也没有讨论声波时差各向异性对速度分析方法影响的强弱；
蔡琳等[10]在储层测井资料评价中，发现碳酸盐岩由于上覆地层压力与围压的不同，导致碳酸盐岩电各向异性较强；
但是并没有考虑该储层的声各向异性。综上所述，以往的专家和学者仅从电学或声学方面评价了储层的各向异性，其结果往往是不全面的。
本文以四川盆地高石—磨溪地区野外露头白云岩为实验对象，采用物理实验的方法，主要分析在高温高压下岩心声各向异性系数与电各向异性系数的实验关系和不同含水饱和度下声电各向异性系数的大小变化；
通过以上研究以期为白云岩储层声波、电阻率差异校正提供帮助，为白云岩各向异性评价提供更多方法。
岩心声电各向异性系数可以定量表征出岩心不同方向上声波和电阻率的差异大小，该表征方法是白云岩声电各向异性评价的基础。
1.1 岩心尺寸的选取
以往电各向异性的研究多是分析在不同倾角下，砂岩岩心电阻率差异的问题。其实，岩石电阻率的差异不仅仅与地层倾角有关，与储层温度、压力和方向等因素也密切相关。本次实验为了保证岩心数据与测井资料有较好的可比性，同时为了满足岩心夹持器的尺寸要求，选用了大尺寸方形岩心为实验对象，具体尺寸大小为50 mm×50 mm×50 mm。
1.2 电各向异性表征方法
图1是双侧向电极系在水平井中的测量示意。其中，A0是主电极，M1和M2是监督电极，A1和A2是屏蔽电极。如图1所示，设与储层层理平行的方向为x方向，该方向电阻率为Rx，单位为Ω·m；
与储层层理垂直的方向为z方向，该方向电阻率为Rz，单位为Ω·m。如图2所示，实验所用岩心较致密且没有裂缝，y方向电阻率与x、z方向电阻率差异不大，所以不考虑y方向电阻率的大小。本实验电各向异性系数λ[11]的表征方法为：
图1 水平井各向异性储层测井示意
。
(1)
1.3 声各向异性表征方法
声各向异性系数表征方法一般分为两类，第一类是当岩心无裂缝时，通常利用岩心水平方向纵波速度和垂直方向纵波速度的比值表征声各向异性系数；
第二类是当岩心存在明显的裂缝和孔洞时，用快慢横波表征岩心的声各向异性系数。本次实验所用岩心较致密且无裂缝孔洞，选用第一类方法表征声各向异性较为合适。如图2所示，设与岩心纹层平行的方向为x方向，该方向纵波波速为Vx，单位为m/s; 与纹层垂直的方向为z方向，该方向纵波波速为Vz，单位为m/s，则声各向异性系数α的表征方法为：
图2 岩心x方向和z方向示意
。
(2)
本节将介绍实验仪器的实验原理、实验方案设计和误差分析，科学严谨的实验流程将提高实验数据的真实性。
2.1 声波实验测量原理
图3是声波实验测量原理，如图所示，左边的激发探头将高压脉冲信号转换成高频机械振动，机械振动沿水平轴穿过岩心，由接收探头转换为脉冲信号输出，最后由计算机和示波器记录和读取声波波形。实验前选用声波探头对接法测量声波的迟滞时间T0，这种方法操作简便，也可以测试实验仪器的测量效果。经过测量，本次实验纵波探头T0值为7.12 μs，横波探头T0值为22.12 μs。
图3 声波脉冲透射法测量原理
2.2 电阻率实验测量原理
实验所用的电阻率测量仪并没有采用古典电桥，而是采用微处理器构成的新式电桥测量电阻。该电桥避免了古典电桥操作复杂、程序步骤多等缺点，程序简便，可自动测量被测元件的电阻值，有效提高电阻率测量仪的测量效率。
如图4所示，正弦波发射器S输出电流Ix流过被测岩心和电阻为rS的标准电阻，由具有相同增益K的差分放大器分别输出电压e1和e2，通过微机自动输出被测元件的电阻值rx。具体表达式为：
图4 电阻率测量仪基本测量原理
rx=rS[e1/e2]
。
(3)
2.3 方形岩心夹持器
如图5所示，方形岩心夹持器采用三轴加压系统控制夹持器的压力，采用温度控制器控制夹持器的温度。声波探头的测量面制作成和岩心测量面大小一样的方形形状，即方便测量该方向岩心的声波数据也方便装载和卸载岩心；
在测量声波的同时电极片测量岩心另一个方向的电阻，实现不同方向的声电同步测量。
图5 实验仪器连接示意
2.4 实验操作过程及误差分析
该实验步骤较多较复杂，操作过程中应尽量减少每个环节的操作误差，保持数据的准确性。
2.4.1 实验操作过程
1)测量前对所有岩心进行洗油、洗盐处理，清理掉岩心中残余的杂质；
2)将所有岩心置于80℃的恒温干燥箱中烘干24 h，测量岩心的干重和尺寸；
3)本次实验配制10 000 ppm的NaCl溶液饱和岩心。测量岩心饱和后的重量并计算孔隙度和密度。由于高温高压下实验极易导致夹持器中胶套破损，影响实验安全，最终确定实验温度和压力分别为60 ℃和60 MPa；
4)将第一块岩心沿x方向放入夹持器中密封，在60 ℃和60 MPa下测量x方向的声波和z方向的电阻；
5)沿z方向将岩心放入夹持器中密封，在60 ℃和60 MPa下测量z方向的声波和x方向电阻；
6) 完成所有岩心声波和电阻率的测量后，根据式(1)、(2)、(4)和(5)计算声电各向异性系数。
声波速度计算公式为：
，
(4)
式中L为岩心测量方向的长度，单位为m；
T1为声波时差，单位为s；
T0为声波迟滞时差，单位为s。
电阻率计算公式为：
R=(r×S)/L
，
(5)
式中r是测量电阻值，单位为Ω；
S是岩心样品测量面的横截面积，单位为m2；
L是岩心测量方向上的长度，单位为m；
R是电阻率值，单位为Ω·m。
2.4.2 实验误差分析
1)测量误差：岩心的尺寸、干重、湿重等数据测量不准确造成测量误差。
2)岩心制备误差：岩心在切割、打磨以及运送过程中的磕碰导致岩心端面不平整造成制备误差。
3)岩心饱和误差：全直径方形岩心尺寸较大难以达到完全饱和而造成饱和误差。
本节将根据所采集的8块岩心声波波形和电阻率，总结归纳声各向异性系数与电各向异性系数实验关系和含水饱和度对声电各向异性系数的影响规律[12-15]。
3.1 岩心物性测量结果
表1是8块岩心测量参数值，8块岩心平均孔隙度为2.724%，平均密度为2.78 g/cm3。其中2号岩心孔隙度较小，密度较大；
原因是岩心制备误差和测量误差所导致的。
表1 岩心物性参数测量结果统计
3.2 声波波形各向异性分析
本节将根据采集的声波波形分析不同方向白云岩岩心声波幅度和声波速度差异[16-22]。
由图6和图7可知2号岩心x方向的纵波波速为7 280.8 m/s，横波波速为4 203.6 m/s；
z方向的纵波波速为7 062.1 m/s，横波波速为4077.4 m/s。图8是2号岩心x方向和z方向的声波波形对比，分析该图可知，同一块岩心x方向和z方向声波速度和声波幅度大小都不一样，均表现出各向异性，归纳为以下几个原因：
图6 2号岩心x方向声波波形
图7 2号岩心z方向声波波形
图8 2号岩心声波波形对比
1)内摩擦消耗
沿z方向穿过岩心的声波与岩石骨架、硅质条带和孔隙流体产生内摩擦消耗了一部分声波的能量；
沿x方向穿过岩心的声波，产生的内摩擦较小，减少了声波能量的消耗；
所以z方向的声波波形相比于x方向的声波波形幅度小、声波速度慢。
2)散射消耗
声波从z方向穿过岩心时，由于岩石骨架和硅质条带中分子的非均质性，造成z方向的声波一部分由于散射消耗掉了；
x方向声波穿过岩心时，散射衰减的声波能量较少。所以，x方向的声波幅度较大，声波速度较快；
而z方向的声波幅度较小，声波速度较慢。
表2是8块岩心样品x方向和z方向纵波速度、横波速度以及声各向异性系数统计表。如表所示，8块岩心x方向和z方向的纵波速度和横波速度差异较小，纵波各向异性系数和横波各向异性系数在1～1.05之间。
表2 白云岩岩心声各向异性系数统计
3.3 电阻率各向异性分析
表3是8块岩心x、z方向电阻率以及电各向异性系数统计结果，如表所示，x方向和z方向电阻率差异较大，电各向异性系数在1.2～1.75之间。
表3 白云岩岩心电各向异性系数统计
3.4 声电各向异性系数实验关系分析
对8块岩心的声、电各向异性系数进行整理分析，得到了60 ℃、60 MPa下声电各向异性系数实验关系图(图9)和表达式：
图9 声各向异性系数与电各向异性系数实验关系
y=53.561x-53.871,R2=0.9534
，
(6)
式中，y为电各向异性系数，无单位；
x为声各向异性系数，无单位。
由图可知声、电各向异性系数呈正相关线性关系，且相关系数R2=0.9534，声电各向异性系数拟合关系较好。该关系表明利用声波和电阻率评价各向异性是存在联系的，白云岩声电各向异性系数可以相互转化，实现声电多参数各向异性评价方法。
3.5 不同含水饱和度下岩心声电各向异性系数实验关系分析
研究1号岩心、3号岩心和6号岩心在不同含水饱和度下声、电各向异性系数规律关系。在改变岩心含水饱和度的过程中，因为岩心孔隙度较小，较致密；
且实验仪器水平放置，采用氮气驱替改变岩心的含水饱和度效果不理想。采用自然风干法改变岩心的含水饱和度效果较好。风干后重新测量岩心的重量，计算岩心的含水饱和度S，测量岩心不同含水饱和度下的声波波形和电阻，按照式(1)、(2)、(4)、(5)得到不同含水饱和度下岩心声电各向异性系数[23-31]。
经过实验多次测量，得到了不同含水饱和度下声波波形和不同含水饱和度下声电各向异性系数，如图10和表4所示。总结得到以下认识：
图10 1号岩心不同含水饱和度x方向声波波形对比
表4 3块岩心不同含水饱和度下声、电各向异性系数统计
1)岩心随着含水饱和度的升高，电阻率下降，且x方向和z方向的电阻率差异增大，电各向异性系数增大。随着含水饱和度的上升，孔隙中游离的Na+离子和Cl-离子增多，电流更容易通过岩心，使岩心电阻率下降；
z方向的电流垂直经过纹层方向，需要穿过电阻较高的硅质条带，使z方向电阻率高于x方向电阻率，且z方向电阻率与x方向电阻率差异增大，电各向异性系数增大；
2)含水饱和度越高的岩心，声波速度上升，x方向和z方向的声波速度差异增大，声各向异性系数增大。干燥岩心x方向和z方向声波速度较为接近，声各向异性系数较小；
随着含水饱和度增加，x方向和z方向声波速度上升，但是z方向声波垂直经过岩石骨架和孔隙流体时摩擦消耗和散射消耗较大，导致z方向声速与x方向声速差异增大，声各向异性系数变大；
3)3块岩心中3号岩心的孔隙度最高，完全饱和水后电各向异性系数最大。说明除去含水饱和度对电各向异性系数的影响外，岩心的孔隙度大小对电各向异性系数也有影响。若要获得更加详细的孔隙度与电各向异性系数的实验关系，则需要更多的岩心和更深入的实验研究来探索。
1)通过物理实验测量，得到8块岩心孔隙度、声波波形和电阻率等数据。研究发现，孔隙度最小的2号岩心，声、电各向异性系数最小；
孔隙度最大的8号岩心，声、电各向异性系数最大；
同时求得8块岩心平均声各向异性系数α=1.033，平均电各向异性系数λ=1.47。
2)通过实验发现，岩心在60 ℃、60 MPa下声、电各向异性系数实验关系为正相关线性关系，关系式为λ=53.561α-53.871，相关程度为0.953，且该关系式在100% 含水饱和度时误差最小。该关系式说明白云岩声各向异性与电各向异性是内在相关的，声、电各向异性强弱程度是一致的。该关系式为水平井和直井声波、电阻率差异建立联系，实现声电多参数各向异性评价方法。
3)对比1号岩心、3号岩心和6号岩心不同含水饱和度下声、电各向异性系数后发现，随着含水饱和度的上升，声各向异性系数变大，电各向异性系数变大。声各向异性系数随含水饱和度平均变化率为0.023，电各向异性系数随含水饱和度平均变化率为0.216。该实验结果说明含水饱和度是声电各向异性评价方法中的重要参数，且电各向异性系数对含水饱和度的变化较为敏感。
猜你喜欢 含水岩心声波 保压取心工具连续割心系统设计中国石油大学学报（自然科学版）(2022年4期)2022-09-05镇北油田某油藏延长低含水采油期技术研究化工管理(2021年7期)2021-05-13基于含水上升模式的海上油田水平井稳产调控策略复杂油气藏(2021年4期)2021-03-08——以双河油田Eh3Ⅳ5-11岩心为例">交联聚合物在岩心孔隙中长期滞留性能研究——以双河油田Eh3Ⅳ5-11岩心为例石油地质与工程(2020年2期)2020-04-17基于声波检测的地下防盗终端电子制作(2018年23期)2018-12-26声波杀手小猕猴智力画刊(2017年6期)2017-07-03声波实验小溪流（画刊）(2017年5期)2017-06-15岩心对复配型驱油剂采油效率的影响当代化工(2015年7期)2015-10-21十几年后的真相小小说月刊(2015年4期)2015-04-15浅议地质岩心实物档案管理办公室业务(2015年23期)2015-03-19
