基于声发射分形参数分析的储层岩石裂纹演化特征研究

唐 杰,王 浩,温 雷,孙成禹,关 键

(1.中国石油大学(华东)地球科学与技术学院,山东青岛266580;2.中国石油化工股份有限公司胜利油田分公司物探研究院,山东东营257022)

唐杰,王浩,温雷,等.基于声发射分形参数分析的储层岩石裂纹演化特征研究[J].石油物探,2017,56(6):-781

TANG Jie,WANG Hao,WEN Lei,et al.Evolution characteristics of reservoir rock cracks based on acoustic emission fractal parameter analysis[J].Geophysical Prospecting for Petroleum,2017,56(6):-781

基于声发射分形参数分析的储层岩石裂纹演化特征研究

唐 杰1,王 浩1,温 雷1,孙成禹1,关 键2

(1.中国石油大学(华东)地球科学与技术学院,山东青岛266580;2.中国石油化工股份有限公司胜利油田分公司物探研究院,山东东营257022)

研究储层岩石压裂过程中所产生的信号特征,能够分析岩石压裂过程中的裂纹演化特征。对不同类型的储层岩石进行单轴受压破坏全过程的声发射实验,结合声发射信号振幅、上升时间及持续时间的累积频度和离散频度分形b值分析,研究了声发射信号在岩石破坏全过程中的变化规律。研究结果表明:①不同类型的岩石受其内部结构特征的影响,应力、应变状态有显著差异,声发射参数分形b值在不同的应力加载阶段分布特征有明显的差异;②振幅、上升时间以及持续时间的分形b值在储层岩石加载的不同阶段表现出明显的规律性。室内试验结果和现场微地震研究结果都表明,分形b值的研究对于致密储层水压致裂过程中的裂纹演化研究具有一定的指导意义。

裂纹演化;声发射;分形;b值;Kaiser效应

致密储层需要通过压裂改造才能形成裂缝网络,水力压裂可以激活天然裂缝和岩石层理,进而实现储层改造[1]。通过对储层岩石压裂过程中所产生的信号开展特征研究,能够分析压裂过程中的震源特征和裂纹演化特征,这对于致密储层的石油勘探开发具有重要意义[2-3]。针对震源特征,朱海波等[4]研究了水力压裂微地震监测的震源机制反演方法,分析了利用矩张量描述震源属性的效果。野外水压致裂和实验室压裂试验虽然尺度不同,但都会产生破裂信号(微地震信号和声发射信号)。要研究微地震信号的演化特征,可以首先分析小尺度的实验室岩石声发射实验测试结果,为分析微地震演化特征提供参考[5-6]。岩石在应力加载过程中以及产生宏观破裂时都会产生声发射,针对岩石在不同加载条件和过程中的声发射特征,WUESTEFELD等[7]做了测试分析,结果表明,在岩石破裂的不同过程中声发射参数会发生明显的变化。尹贤刚等[8]分析了岩石在加载破坏过程中声发射分形特征的变化情况;李元辉等[9]分析了岩石破裂过程中声发射b值的变化特征,结果表明,分形b值与岩石内部微裂纹的演化有着密切的联系。以往的研究大多是针对花岗岩、玄武岩等,而致密储层大多是致密砂岩、泥岩或白云岩等,因此有必要对其破裂过程中的声发射分形演化特征展开研究[10]。

本文开展了泥岩、白云岩、砂岩等不同岩性的岩石单轴压裂试验,采用声发射仪器获得了储层岩石破裂过程中的声发射参数,利用分形b值研究了储层岩石应力加载以及宏观破裂过程中的裂纹演化特征。

1 岩石声发射实验

声发射事件(AE)是岩石内部产生微裂纹时释放能量所导致的,因此声发射与岩石内部裂纹的演化发展密切相关[11-12]。本文选取泥岩、白云岩、砂岩3种类型的岩石,利用电子伺服压机进行单轴压裂试验,采用PCI-II型声发射记录仪采集了储层岩石应力加载和宏观破裂过程中的声发射参数。

1.1 实验设备及装置

首先选择测试岩心,包括白云岩、砂岩以及泥岩。岩石加工成直径25mm,长50mm的岩心,岩心两端采用打磨和抛光设备加工平整,以满足实验测试需要。在圆柱形岩样的四周安置4个微声发射探头。应力加载仪器采用MTS刚性压机,加载方法为单轴多级循环加载,加载仪器采用轴向变形位移控制。使用PCI-II型声发射记录系统采集应力加载和宏观破裂过程中产生的声发射波形及其它参数信息(如图1所示)。PCI-II型声发射采集系统具有高速采集能力,其采样率高达40MHz,此外还具有记录连续波形的能力。

图1 岩石破裂声发射测试系统

1.2 实验结果分析

图2给出了白云岩、泥岩和砂岩在应力加卸载过程中应力与应变的关系曲线。由图2可知,不同类型的岩石应力、应变状态有显著差异。白云岩较为致密,强度大,弹性特征明显,加载过程中应力、应变曲线接近线性;泥岩内部含有较多的粘土质矿物,塑性相对较强,加载过程中应力、应变曲线的非线性特征较为明显,在应力卸载变为0时,应变不能完全恢复之前的状态;砂岩内部相对而言孔隙较多,孔隙度较大,加载过程中应力、应变曲线的非线性特征也较为明显,当应力加载到一定强度,原本存在的一些较薄的裂隙会闭合,在应力减小的时候也不能完全恢复。

图3给出了白云岩、泥岩和砂岩的声发射(AE)振幅与应力关系曲线;图4给出了实验过程中白云岩、泥岩和砂岩的声发射(AE)持续时间与应力变化关系曲线;图5给出了白云岩、泥岩和砂岩的声发射(AE)上升时间与应力变化关系曲线。由图3到图5可知,不同类型的岩石声发射参数演化特征各异,这主要是由不同类型岩石内部颗粒特征和弱结构面特征不同所导致的。砂岩颗粒相对较大,其声发射数目明显多于白云岩和泥岩;白云岩所含的矿物一般为白云石、长石、石英等脆性矿物,其脆性特征较为明显,当岩石应力加载到后期,岩体内部出现较大裂缝,记录到的声发射事件数量急剧增加;泥岩是由黏土矿物固化而成,具有较强的塑性,其声发射信号的持续时间和上升时间大于白云岩和砂岩。从图3到图5还可以看出,储层岩石在应力加载过程中会表现出Kaiser效应[13],在岩石重复加载的过程中,如果再次加载的应力不超过前次加载的最大应力,则岩石不会有明显的声发射现象,当再次加载的应力超过前次加载的最大应力时,岩石声发射现象又会再次出现。

图2 白云岩(a)、泥岩(b)、砂岩(c)的应力与应变关系曲线

图3 白云岩(a)、泥岩(b)、砂岩(c)的AE振幅与应力关系曲线

图4 白云岩(a)、泥岩(b)、砂岩(c)的AE持续时间与应力关系曲线

图5 白云岩(a)、泥岩(b)、砂岩(c)的AE上升时间与应力关系曲线

2 实验数据的声发射分形特征分析

2.1 累积频度与离散频度分形b值

岩石破坏过程中内部表现出微裂纹产生、扩展和形成大破裂等演化过程,在裂纹演化过程中伴随着声发射信号的产生,GUTERBERG等[14]给出了声发射信号大小与发生频度之间的统计关系式:

(1)

式中:M=AdB/20,为声发射信号的大小,AdB是采用分贝表示的声发射信号的最大振幅;N是信号大小大于或等于M时的累积声发射数目;a和b是常数,a代表声发射的活动性,主要取决于序列中的最大声发射信号大小,b值反映了岩石承受应力和接近强度极限的程度。研究表明,b值具有分形特征,分形b值随着岩石内部微裂纹的演化发生明显的变化,b值的增大或减小反映了不同尺寸微裂纹所占比例的变化情况。采用最大似然估计可以计算b值:

(2)

式中:∑(nilgAi)为声发射信号的振幅对数和,Ai为单个声发射信号的振幅大小,ni为不同阶段声发射有效信号的数量;Amin为声发射信号的最小振幅;n是声发射信号总数。将(2)式扩展,可以分析上升时间T以及持续时间t的分形特征参数b1和b2:

式中:Ti为单个声发射信号的上升时间;ti为单个声发射信号的持续时间;Tmin为计算对象中的最小上升时间;tmin为计算对象中的最小持续时间。

2.2 不同应力阶段分形b值计算结果

针对声发射测试,可以采用离散频度方法和累积频度方法求取b值。采用离散频度方法求取b值时,采取了滑动平均的方式,选取150个有效声发射信号作为一组数据,采用的滑动窗口大小为40个声发射事件。图6给出了不同类型储层岩石在破裂过程中采用两种方法得到的声发射振幅分形b值变化特征。由图6可知,在开始进行应力加载的阶段,计算所得到的累积频度和离散频度分形b值相对较大,这主要是由于这一阶段振幅较小的信号所占比例较多,也就是说这一阶段所产生的破裂尺度很小;随着加载应力的增加,累积频度振幅分形b值的规律性开始增强,振幅分形b值逐渐减小,这表明储层岩石内部相对较大的微裂纹开始产生,而且较大裂纹所占的比例也逐步增大;随着加载应力的持续增大,累积频度声发射信号振幅分形b值在岩心破坏前下降到最低值,此时岩石产生明显的宏观破裂面。结合LEI等[15-16]的研究结果发现,岩石破坏前b值下降的物理机制一是微破裂从以张破裂为主转为以剪破裂为主;二是裂纹相互作用明显增强[17]。比较不同类型的储层岩石分形b值的演化特征可知,白云岩在破坏前累积频度声发射振幅分形b值急剧下降,离散频度声发射振幅分形b值在形成破裂面后增加;砂岩其累积频度声发射振幅分形b值随着加载逐步下降,离散频度声发射振幅分形b值在加载过程中扰动明显,这是由于砂岩在加载过程中容易产生微破裂与宏观破裂面共同作用的结果。

图7和图8分别给出了不同类型岩石在破裂过程中声发射上升时间与持续时间分形b值的变化特征。由图7和图8可以发现,对于不同岩性的岩石而言,累积频度方法获得的上升时间和持续时间b值变化趋势基本一致,离散频度法b值变化趋势有所不同。对于泥岩和砂岩来说,总体趋势呈现缓慢的下降趋势,在破裂发生的时间点b值最小;而对于白云岩来说,声发射上升时间b值先呈下降趋势,在破裂时间点达到极小值,然后增大再减小,声发射信号持续时间有相同的变化趋势。研究表明,声发射参数的分形b值是评价岩石破坏过程的有效指标。

图6 岩石破裂过程中AE振幅的分形b值变化情况 a 白云岩; b 砂岩; c 泥岩

图7 岩石破裂过程中AE上升时间的分形b值变化情况 a 白云岩; b 砂岩; c 泥岩

图8 岩石破裂过程中AE持续时间的分形b值变化情况 a 白云岩; b 砂岩; c 泥岩

3 结论

对泥岩、白云岩和砂岩进行单轴加载压裂声发射试验,分析讨论了岩石破裂过程中的声发射特性,得出以下结论。

不同类型的岩石受内部结构特征的影响,其应力、应变状态有显著差异,声发射振幅参数的分形b值在不同的应力加载阶段分布特征有明显的差异。应力初始加载阶段,分形b值相对较大;随着加载应力的增加,分形b值开始规律性增强,之后逐渐减小;随着加载应力的持续增大,分形b值在岩心破坏前下降到最低值。

在储层岩石加载的不同阶段,上升时间和持续时间的分形b值表现出明显的规律性,总体呈现缓慢的下降趋势,在破裂发生的时间点b值最小。

考虑到水压致裂的特点,未来可以开展实验室水压致裂过程中的声发射特征分析。此外。近年来国内外对于水压致裂过程中微地震信号的分形b值做了一些研究和分析,室内试验结果和现场微地震研究结果都表明,分形b值的研究对于致密储层水压致裂过程中的裂纹演化特征研究具有一定的指导意义。下一步将结合水力压裂诱发的微地震数据进行声发射特征分析。此外,可以开展其它的分形、分维参数的研究,进一步认识裂纹的演化过程。

[1] 朱海波,杨心超,廖如刚,等.基于微地震裂缝参数反演的解释与应用研究[J].石油物探,2017,56(1):150-157

ZHU H B,YANG X C,LIAO R G,et al.Microseismic fracture interpretation and application based on parameters inversion[J].Geophysical Prospecting for Petroleum,2017,56(1):150-157

[2] 杨心超,朱海波,李宏,等.基于P波辐射花样的压裂微地震震源机制反演方法研究及应用[J].石油物探,2016,55(5):640-648

YANG X C,ZHU H B,LI H,et al.Microseismic focal mechanism inversion based on P-wave radiation pattern and its application[J].Geophysical Prospecting for Petroleum,2016,55(5):640-648

[3] 刘勇,方伍宝,李振春,等.基于叠前地震的脆性预测方法及应用研究[J].石油物探,2016,55(3):425-432

LIU Y,FANG W B,LI Z C,et al.Brittleness prediction and application based on pre-stack seismic inversion[J].Geophysical Prospecting for Petroleum,2016,55(3):425-432

[4] 朱海波,杨心超,王瑜,等.水力压裂微地震监测的震源机制反演方法应用研究[J].石油物探,2014,53(5):556-561

ZHU H B,YANG X C,WANG Y,et al.The application of microseismic source mechanism inversion in hydraulic fracturing monitoring[J].Geophysical Prospecting for Petroleum,2014,53(5):556-561

[5] 梁天成,葛洪魁,郭志伟.利用声发射和波速变化判定岩石损伤状态[J].中国地震,2012,28(2):154-166

LIANG T C,GE H K,GUO Z W.Evaluation of rock damage state with acoustic emission and velocity variation[J].Earthquake Research in China,2012,28(2):154-166

[6] 刘玉海,尹成,潘树林,等.基于互相关函数法的地面微地震信号检测研究[J].石油物探,2012,51(6):633-637

LIU Y H,YIN C,PAN S L,et al.Ground microseismic signal detection based on cross-correlation function[J].Geophysical Prospecting for Petroleum,2012,51(6):633-637

[7] WUESTEFELD A,KENDALL J M,VERDON J P,et al.In situ monitoring of rock fracturing using shear wave splitting analysis:an example from a mining setting[J].Geophysical Journal International,2011,187(2):848-860

[8] 尹贤刚,李庶林,唐海燕.岩石破坏声发射强度分形特征研究[J].岩石力学与工程学报,2005,24(19):3512-3516

YIN X G,LI S L,TANG H Y.Study on strength fractal features of acoustic emission in process of rock failure[J].China Journal of Rock Mechanics and Engineering,2005,24(19):3512-3516

[9] 李元辉,刘建坡,赵兴东,等.岩石破裂过程中的声发射b值及分形特征研究[J].岩土力学,2009,30(9):2559-2563

LI Y H,LIU J P,ZHAO X D,et al.Study onb-value and fractal dimension of acoustic emission during rock failure process[J].Rock and Soil Mechanics,2009,30(9):2559-2563

[10] 宋维琪,喻志超,杨勤勇,等.低信噪比微地震事件初至拾取方法研究[J].石油物探,2013,52(6):596-608

SONG W Q,YU Z C,YANG Q Y,et al.The first arrival picking method of microseismic events with low S/N[J].Geophysical Prospecting for Petroleum,2013,52(6):596-608

[11] MAIN I G.Damage mechanics with long-range interactions:correlation between the seismicb-value and the fractal two-point correlation dimension[J].Geophysical Journal International,1992,111(3):531-541[12] HATTON C G,MAIN I G,MEREDITH P G.A comparison of seismic and structural measurements of scaling exponents during tensile subcritical crack growth[J].Journal of Structural Geology,1993,15(12):1485-1495

[13] KAISER E J.A study of acoustic phenomena in tensile test[D].Munich Germany:Technische Hochschule Munchen,1953

[14] GUTERBERG B,RICJHTER C.Seismicity of the earth and associated phenomena[M].USA:Princeton University Press,1954:1-284

[15] LEI X L,KUSUNOSE K,NISHIZAWA O.Quasi-static fault growth and cracking in homogeneous brittle rock under triaxial compression using acoustic emission monitoring[J].Journal of Geophysical Research,2000,105(B3):6127-6139

[16] LEI X L,SATOH T S.Indicators of critical point behavior prior to rock failure inferred from pre-failure damage[J].Tectonophysics,2007,431(1/2/3/4):97-111

[17] 唐杰,方兵,孙成禹,等.压裂诱发的微地震震源机制及信号传播特性研究[J].石油地球物理勘探,2015,50(4):643-649

TANG J,FANG B,SUN C Y,et al.Study on focal mechanism of micro-seismic induced by hydrofracture and signal propagation characteristics[J].Oil Geophysical Prospecting,2015,50(4):643-649

(编辑：陈 杰)

Evolutioncharacteristicsofreservoirrockcracksbasedonacousticemissionfractalparameteranalysis

TANG Jie1,WANG Hao1,WEN Lei1,SUN Chengyu1,GUAN Jian2

(1.SchoolofGeosciences,ChinaUniversityofPetroleum,Qingdao266580,China;2.GeophysicalResearchInstitute,ShengliOilfieldBranchCo.,SINOPEC,Dongying257022,China)

Acoustic emission signal characteristics can be used to analyze the characteristics of crack evolution during reservoir rock fracturing process.Using acoustic emission tests under uniaxial compression on different types of reservoir rocks,the change law of the acoustic emission signal with stress load variation in the whole process of rock failure was studied by analyzing the fractal b values of amplitude,rising time,and duration of the acoustic emission signal obtained using the cumulative frequency and discrete frequency methods.The test results show that,during uniaxial compression process,the stress and strain state of different types of rocks differ significantly,because they have different internal structures,and the distribution characteristics of the fractal b value in different stress loading stages differ obviously.Moreover,the fractal b values of amplitude,rising time,and duration show obvious regularity at different stress loading stages.Both the laboratory experiment and the field microseismic test results show that the analysis of fractal b value can be used to clarify crack evolution characteristics in tight reservoir hydraulic fracturing.

crack evolution,acoustic emission,fractal,b value,Kaiser effect

2017-02-13;改回日期2017-05-23。

唐杰(1980—),男,博士,副教授,主要从事地震岩石物理学研究工作。

国家自然科学基金项目(41504097,41374123)、国家科技重大专项(2016ZX05006-002)与中央高校基本科研业务费(15CX08002A)联合资助。

This research is financially supported by the Natural Science Foundation of China (Grant Nos.41504097,41374123),the National Science and Technology Major Project of China (Grant No.2016ZX05006-002) and the Fundamental Research Funds for the Central Universities (Grant No.15CX08002A).

P631

A

1000-1441(2017)06-0775-07

10.3969/j.issn.1000-1441.2017.06.001