本发明涉及一种石油天然气开采技术，尤其涉及一种用于评价流体对岩石造缝能力的实验方法。
背景技术：
水力压裂技术是提高页岩气、致密油气等非常规储层产量的重要手段，也是提高干热岩热能开采的有效方式之一。水力压裂过程中要向地层中注入大量的流体(压裂改造工作液)，流体可以显著改变岩石力学性质，因此，流体对岩石的造缝能力(或可压性)影响不可忽视。
目前评价岩石可压性的方法主要通过脆性来评价，认为脆性越强，可压性越强；也有学者发展了新的可压性评价方法，例如通过断裂韧性和脆性指数来综合评价岩可压性，但是目前方法忽略了流体对岩石造缝能力的大小。
技术实现要素：
本发明的目的是提供一种用于评价流体对岩石造缝能力的实验方法。
本发明的目的是通过以下技术方案实现的：
本发明的用于评价流体对岩石造缝能力的实验方法，包括以下步骤：
步骤1：将直径25毫米、高度50mm的标准圆柱体岩心放入恒温烘干箱中，温度恒定为100摄氏度，连续烘干24小时，以除去岩心中的残余液体；
步骤2：测试干燥岩心的高度h、直径d和质量m，计算干燥岩心的密度ρ，公式如下：
式中，ρ为岩心密度；d为岩心直径；h为岩心高度；π为圆周率，常数；m为岩心质量；
步骤3：测试干燥岩心沿轴向的纵波vp、横波vs的大小，计算其动态弹性模量ed和动态泊松比νd的大小，公式如下：
式中，ed为动态弹性模量；νd为动态泊松比；vp为纵波(或p波)波速；vs为横波(或s波)波速；
步骤4：将第3步骤中获得的动态岩石力学参数转换为静态岩石力学参数，转换公式如下：
es＝aed+b(4)
νs＝cνd+d(5)
式中，a、b、c、d为常数，通过对实验数据拟合来确定；es为静态弹性模量；νs为静态泊松比；
步骤5：计算岩心的脆性指数bi，公式如下：
式中，bi为岩心脆性指数；ebrit为归一化弹性模量；νbrit为归一化泊松比；es,max、es,min为步骤4中所有岩心静态弹性模量es构成向量中的最大元素值和最小元素值；νs,max、νs,min为步骤4中所有岩心静态泊松比νs构成向量中的最大元素值和最小元素值；
步骤6：将p波(或s波)探头置放于探头夹具内，探头端面与岩心两个端面接触，中间均匀涂上耦合剂，以一定的加载速率沿轴向加载，最大加载应力σmax为岩心单轴压缩强度的30％至50％，测试干燥岩心在不同应力下的波形和波速到时大小tp，然后根据互相关函数法，得出不同应力σ下的波速大小vp,str，同时记录相应应力σ下产生的声发射事件个数num，将单位应力下声发射事件数定义为aenum＝num/σ；
步骤7：对每块岩心，绘制不同加载应力σ与波速vp,str间的拟合直线，直线斜率kvel为岩心的波速-应力敏感性系数；
步骤8：根据获得波速数据，通过反演计算，获得每块岩心的天然裂隙密度大小ρc，反演公式如下：
式中，g0为岩石基质的剪切模量；gdry为干燥岩样的剪切模量；g0和gdry可由式(10)计算获得；ν0为岩石基质的泊松比，可由式(3)获得；ρc为岩样的天然裂隙密度；e为岩心弹性模量；ν为岩心泊松比；
步骤9：将干燥岩心浸泡于压裂改造工作液中一定时间(为压裂施工时间)，重复步骤二至步骤八，获得润湿岩心的脆性指数biwet、单位应力下声发射事件数aenum,wet、波速-应力敏感性系数kvel,wet以及天然裂隙密度ρwet,c的大小，其中润湿岩样的天然裂隙密度ρwet,c的计算公式如下(其余参数同上述步骤中的公式)：
式中，gwet为润湿岩样的剪切模量，可由公式(10)计算得出；δf为无量纲量，由公式(12)可得，其中k0和kf分别为岩石基质和流体(即压裂改造工作液)的体积模量；e0为岩石基质的弹性模量；ζ为裂隙纵横比，其值大小为10-3左右，ζ＝w/c，w和c分别为裂纹的宽度和半径；ρwet,c为润湿岩心内的天然裂隙密度；
步骤10：计算每块岩样干燥与润湿前后的脆性系数相对变化量bivar、单位应力下声发射事件数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及天然裂隙密度相对变化量ρvar大小，公式如下：
式中，bi和biwet分别是干燥、润湿岩样的脆性指数；bivar为岩样干燥与润湿前后的脆性系数相对变化量；aenum和aenum,wet分别是干燥、润湿岩样在单位应力下的声发射事件数量；aevar为岩样在单位应力下的声发射事件数量相对变化量；kvel和kvel,wet分别是干燥与润湿岩样波速-应力敏感性系数；kvar为岩样波速-应力敏感性系数相对变化量；ρc和ρwet,c为干燥、润湿岩样内的天然裂隙密度；ρvar为天然裂隙密度相对变化量。
步骤11：对于每块岩心，将步骤10中的岩心脆性指数相对变化量bivar、单位加载应力下声发射事件总数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及岩样天然裂隙密度相对变化量ρvar进行无量纲化处理，公式如下：
式中，bid、aed、kd、ρd分别为无量纲化脆性指数、无量纲化声发射数、无量纲化波速-应力敏感性系数和无量纲化天然裂隙密度；bivar,max、bivar,min为步骤10中所有岩心脆性指数相对变化量构成向量中的最大元素值和最小元素值；aevar,max、aevar,min为步骤10中所有岩心单位应力下声发射总数相对变化量构成向量中的最大元素值和最小元素值；
kvar,max、kvar,min为步骤10中所有岩心波速-应力敏感性系数相对变化量构成向量中的最大元素值和最小元素值；ρvar,max、ρvar,min为步骤10中所有岩心天然裂隙密度相对变化量构成向量中的最大元素值和最小元素值；
步骤12：计算流体作用下岩心的造缝能力指数大小fif，公式如下：
步骤13：比较步骤12中的不同流体作用下每块岩心的造缝能力指数大小，其值越大，说明该流体对岩心的造缝能力越强。
由上述本发明提供的技术方案可以看出，本发明实施例提供的用于评价流体对岩石造缝能力的实验方法，可用于优选压裂液类型，促进水力压裂过程中复杂缝网的形成。该方法可以快速既可以获得不同流体对岩石的造缝能力指数，也可以获得同一种流体对不同岩样的造缝能力指数，可以方便地指导水力压裂施工优选压裂液，以促使压裂过程中人工复杂缝网的形成。
附图说明
图1为本发明实施例提供的用于评价流体对岩石造缝能力的实验方法中的应力-波速拟合直线图。
具体实施方式
下面将对本发明实施例作进一步地详细描述。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
本发明的用于评价流体对岩石造缝能力的实验方法，其较佳的具体实施方式是：
包括以下步骤：
步骤1：将直径25毫米、高度50mm的标准圆柱体岩心放入恒温烘干箱中，温度恒定为100摄氏度，连续烘干24小时，以除去岩心中的残余液体；
步骤2：测试干燥岩心的高度h、直径d和质量m，计算干燥岩心的密度ρ，公式如下：
式中，ρ为岩心密度；d为岩心直径；h为岩心高度；π为圆周率，常数；m为岩心质量；步骤3：测试干燥岩心沿轴向的纵波vp、横波vs的大小，计算其动态弹性模量ed和动态泊松比νd的大小，公式如下：
式中，ed为动态弹性模量；νd为动态泊松比；vp为纵波(或p波)波速；vs为横波(或s波)波速；
步骤4：将第3步骤中获得的动态岩石力学参数转换为静态岩石力学参数，转换公式如下：
es＝aed+b(4)
νs＝cνd+d(5)
式中，a、b、c、d为常数，通过对实验数据拟合来确定；es为静态弹性模量；νs为静态泊松比；
步骤5：计算岩心的脆性指数bi，公式如下：
式中，bi为岩心脆性指数；ebrit为归一化弹性模量；νbrit为归一化泊松比；es,max、es,min为步骤4中所有岩心静态弹性模量es构成向量中的最大元素值和最小元素值；νs,max、νs,min为步骤4中所有岩心静态泊松比νs构成向量中的最大元素值和最小元素值；
步骤6：将p波(或s波)探头置放于探头夹具内，探头端面与岩心两个端面接触，中间均匀涂上耦合剂，以一定的加载速率沿轴向加载，最大加载应力σmax为岩心单轴压缩强度的30％至50％，测试干燥岩心在不同应力下的波形和波速到时大小tp，然后根据互相关函数法，得出不同应力σ下的波速大小vp,str，同时记录相应应力σ下产生的声发射事件个数num，将单位应力下声发射事件数定义为aenum＝num/σ；
步骤7：对每块岩心，绘制不同加载应力σ与波速vp,str间的拟合直线，直线斜率kvel为岩心的波速-应力敏感性系数；
步骤8：根据获得波速数据，通过反演计算，获得每块岩心的天然裂隙密度大小ρc，反演公式如下：
式中，g0为岩石基质的剪切模量；gdry为干燥岩样的剪切模量；g0和gdry可由式(10)计算获得；ν0为岩石基质的泊松比，可由式(3)获得；ρc为岩样的天然裂隙密度；e为岩心弹性模量；ν为岩心泊松比；
步骤9：将干燥岩心浸泡于压裂改造工作液中一定时间(为压裂施工时间)，重复步骤二至步骤八，获得润湿岩心的脆性指数biwet、单位应力下声发射事件数aenum,wet、波速-应力敏感性系数kvel,wet以及天然裂隙密度ρwet,c的大小，其中流体作用下岩样的天然裂隙密度ρwet,c的计算公式如下(其余参数同上述步骤中的公式)：
式中，gwet为润湿岩样的剪切模量，可由公式(10)计算得出；δf为无量纲量，由公式(12)可得，其中k0和kf分别为岩石基质和流体的体积模量；e0为岩石基质的弹性模量；ζ为裂隙纵横比，其值大小为10-3左右，ζ＝w/c，w和c分别为裂纹的宽度和半径；ρwet,c为润湿岩心内的天然裂隙密度；
步骤10：计算每块岩样干燥与润湿前后的脆性系数相对变化量bivar、单位应力下声发射事件数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及天然裂隙密度相对变化量ρvar大小，公式如下：
式中，bi和biwet分别是干燥、润湿岩样的脆性指数；bivar为岩样干燥与润湿前后的脆性系数相对变化量；aenum和aenum,wet分别是干燥、润湿岩样在单位应力下的声发射事件数量；aevar为岩样在单位应力下的声发射事件数量相对变化量；kvel和kvel,wet分别是干燥与润湿岩样波速-应力敏感性系数；kvar为岩样波速-应力敏感性系数相对变化量；ρc和ρwet,c为干燥、润湿岩样内的天然裂隙密度；ρvar为天然裂隙密度相对变化量；
步骤11：对于每块岩心，将步骤10中的岩心脆性指数相对变化量bivar、单位应力下声发射事件总数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及岩样天然裂隙密度相对变化量ρvar进行无量纲化处理，公式如下：
式中，bid、aed、kd、ρd分别为无量纲化脆性指数、无量纲化声发射数、无量纲化波速-应力敏感性系数和无量纲化天然裂隙密度；bivar,max、bivar,min为步骤10中所有岩心脆性指数相对变化量构成向量中的最大元素值和最小元素值；aevar,max、aevar,min为步骤10中所有岩心单位应力下声发射总数相对变化量构成向量中的最大元素值和最小元素值；
kvar,max、kvar,min为步骤10中所有岩心波速-应力敏感性系数相对变化量构成向量中的最大元素值和最小元素值；ρvar,max、ρvar,min为步骤10中所有岩心天然裂隙密度相对变化量构成向量中的最大元素值和最小元素值；
步骤12：计算流体作用下岩心的造缝能力指数大小fif，公式如下：
步骤13：比较步骤12中的不同流体作用下每块岩心造缝能力指数大小，其值越大，说明该流体对岩心的造缝能力越强。
所述非常规储层不仅包括非常规油气藏中的页岩气、煤层气或致密油气储层或干热岩储层，还包括常规油气藏中的低渗透或超低渗透储层。
本发明的用于评价流体对岩石造缝能力的实验方法，可用于优选压裂液类型，促进水力压裂过程中复杂缝网的形成。不仅适用于干热岩、页岩气、煤层气、致密油气等非常规储层，还可应用于低渗透和超低渗透的常规油气藏。该方法可以快速获得不同流体对岩石的造缝能力指数，或同一流体作用下不同岩样的造缝能力指数，指导水力压裂施工优选压裂液，以促使压裂过程中人工复杂缝网的形成。
具体实施例，如图1所示：
评价不同流体对岩石造缝能力的实验方法，岩心来自某油田致密砂岩，具体过程如下：
步骤1：将标准圆柱体岩心放入恒温烘干箱中，在100摄氏度恒温下，连续烘干24小时，除去岩心内的液体；
步骤2：测试干燥岩心的高度h、直径d和质量m，按公式(1)，计算干燥岩心的密度ρ；
步骤3：测试干燥岩心沿轴向的纵波、横波大小，按公式(2)、(3)，分别计算其动态弹性模量ed和动态泊松比νs大小；
步骤4：将第3步骤中获得的动态岩石力学参数转换为静态岩石力学参数，转换公式为：es＝0.85ed和νs＝0.95νd；
步骤5：按公式(6)(8)，先对步骤4中的静态弹性模量es和静态泊松比νs进行归一化，然后计算干燥岩心的脆性指数bi；
步骤6：将p波(或s波)探头置放于探头夹具内，探头端面与岩心两个端面接触，中间均匀涂上耦合剂，以2mpa/次的速率沿轴向加载，最大加载应力σmax为30mpa，测试干燥岩心在不同应力下的波形和波速到时，根据互相关函数法，得出不同应力σ下的波速大小vp,str，同时记录相应应力σ下产生的声发射事件个数num，计算单位应力下的声发射事件数aenum＝num/σ。
步骤7：对每块岩心，绘制不同加载应力σ与波速vp,str间的拟合直线，直线斜率kvel为岩心的波速-应力敏感性系数；
步骤8：根据获得波速数据，通过反演计算公式(9)、式(10)，获得每块岩心的天然裂隙密度大小ρc；
步骤9：将干燥岩心浸泡于蒸馏水中24小时，重复步骤二至步骤八，获得润湿岩心的脆性指数biwet、声发射事件数aenum,wet、波速-应力敏感性系数kvel,wet以及天然裂隙密度ρwet,c的大小，其中饱和流体作用下岩样的天然裂隙密度ρwet,c的计算公式按式(11)、式(12)，其余参数计算公式同步骤2至8中公式；
步骤10：按公式(13)至式(16)，计算每块岩样干燥与润湿前后的脆性系数相对变化量bivar、单位应力下声发射事件数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及天然裂隙密度相对变化量ρvar大小，公式如下；
步骤11：对于每块岩心，按公式(17)至式(20)，将步骤10中的岩心脆性指数相对变化量bivar、声发射事件总数相对变化量aevar、波速-应力敏感性系数相对变化量kvar以及岩样天然裂隙密度相对变化量ρvar进行无量纲化处理；
步骤12：按公式(21)，计算流体对岩心的造缝能力指数大小fif，见表1。
表1不同岩心的造缝能力指数
步骤13：比较步骤12中的流体作用下每块岩心造缝能力指数大小，可以看出编号为10-1的造缝能力指数最大，说明流体对此块岩心的造缝能力最强。
施工效果：
通过较高泵注排量(8～10m3/min)、小粒径支撑剂(40～70目)、滑溜水等工艺措施，对这口井进行了体积压裂施工。施工过程中，在恒定排量注入滑溜水情况下，泵压多次起伏波动，这是由于多次破裂成缝引起的。体积压裂后，此井日产天然气76000m3/天。这说明此口井滑溜水对岩石的造缝能力较好。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。