本发明公开了一种煤储层可压裂性评价方法，包含以下步骤：步骤1：求取待评价煤储层脆性指数BI；步骤2：求取在围压下煤储层断裂韧性KC；步骤3：求取待评价煤层的煤层厚度Eh；步骤4：求取煤层与围岩弹性模量差异Ek’；步骤5：求取待评价煤层的净压力系数σn；步骤6：计算煤岩的水分含量Mad；步骤7：计算煤储层可压裂性评价分数Frac；步骤8：将煤层的可压裂性进行分级。本发明的主要优点在于采用的参数容易获取，对于煤层压裂评价实施难度较低，同时综合考虑了影响煤层压裂的各方面影响因素，能够较为准确的评价可压性，这对于后期的煤层气开发具有重要的指示作用，可以很好地避免无效井位的施工，节约煤层气开发的成本。
1.一种煤储层可压裂性评价方法，其特征在于，包含以下步骤：步骤1：求取待评价煤储层脆性指数BI；
在上述步骤2中，步骤2.1，求取煤储层抗压强度σc，σc＝(0.0045+0.0035Vsh)EKc (4)其中，Vsh是煤储层泥质含量，
Vsh＝(GR-GRmin)/(GRmax-GRmin) (5)EKc是煤储层动态弹性模量，
(4)式中，Vsh为待评价煤层泥质含量；(5)式中，GR为待评价煤层的伽马测井实测值，单位为API，GRmin为测得的井段砂岩层的自然伽马值，单位为API，GRmax为测得的井段泥岩层的自然伽马值，单位为API；(6)式中，EKc为待评价煤层的动态弹性模量，单位为GPa，vp为待评价煤层的纵波速度，单位为km/s，vs为待评价煤层的横波速度，单位为km/s；ρ为岩石密度，单位为g/cm3；(7)式中，A为常数；
KC＝0.0956pw+0.1383St-0.0820 (8)为简化计算，(8)式中围压pw用最小水平主应力代替；
步骤2.4，对待评价煤层的断裂韧性进行正向化、归一化处理，对KC进行正向化、归一化：
Kc’＝(Kc-Kcmin)/(Kcmax-Kcmin) (9)式中，Kc’是正向化的断裂韧性，Kc是待评价煤层的断裂韧性，Kcmin是所有评价煤层的断裂韧性最小值，Kcmax是所有待评价煤层断裂韧性的最大值；
式中，Kc”为正向化、归一化的断裂韧性；Kc’为正向化的断裂韧性；μKc为待评价煤层正向化后断裂韧性的均值，σKc为待评价煤层正向化后断裂韧性的标准差。
2.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤1中，获取煤储层弹性模量和泊松比，对弹性模量和泊松比进行归一化处理，求取待评价煤储层脆性指数BI：(1)式中，EBI’为待评价煤层归一化的弹性模量，EBImax和EBImin分别为煤储层最大和最小弹性模量，单位为GPa；(2)式中，μBI’为煤储层归一化的泊松比，μBImax和μBImin分别为煤储层最大和最小泊松比，无量纲。
3.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤3中，步骤
3.1，利用已有的煤田钻孔以及煤储层气参数井资料，采用内插法线性拟合约束计算研究区域的煤储层厚度，煤储层厚度为Eh，单位m，步骤3.2，对煤储层厚度进行归一化：
式中，Eh’为归一化的煤储层厚度，μEh为待评价煤层煤厚的均值，σEh为待评价煤层煤厚的标准差。
4.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤4中，步骤
4.1，利用围岩与煤储层弹性模量的比值Ek表征煤储层与围岩的弹性模量差异：式中，Ec为煤储层弹性模量，单位为GPa；Et为顶板弹性模量，单位为GPa；Eb为底板弹性模量，单位为GPa，步骤4.2，对弹性模量差异进行归一化：
式中，Ek’为归一化的弹性模量差异，μEk为待评价煤层弹性模量差异的均值，σEk为待评价煤层弹性模量差异标准差。
5.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤5中，步骤
5.1，采用水压致裂法获得煤储层最小水平主应力Pc，单位为MPa；煤储层最大水平主应力σH：σH＝3Pc-Pf-Po+T (14)
式中，Pf为煤储层破裂压力，单位为MPa；Po为煤储层压力，单位为MPa；T为煤储层抗拉强度，单位为MPa；
式中，σn’为归一化的弹性模量差异，μσn为待评价煤层净压力系数的均值，σσn为待评价煤层净压力系数的标准差。
6.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤6中，步骤
6.1，通过对煤岩含水量与测井参数相关关系的分析，选择相关性最为密切的三个参数进行多元回归分析，依此建立了煤岩水分含量的预测方程：Mad＝1.4655-0.5827×DEN-2.1115×GR+0.2319×ρs (17)式中：Mad为待评价煤层煤岩水分含量，单位为％；DEN为待评价煤层密度值，单位为g/cm3；GR为待评价煤层的自然伽马值，单位为API；ρs为待评价煤层的视电阻率值，单位为Ω·m；
步骤6.2，对煤岩水分含量进行正向化、归一化：Mad’＝(Mad-Madmin)/(Madmax-Madmin) (18)式中，Mad’是正向化的煤岩水分含量，Mad是待评价煤层的煤岩水分含量，Madmin是所有评价煤层的煤岩水分含量最小值，Madmax是所有待评价煤层煤岩水分含量的最大值，式中，Mad”为正向化、归一化的煤岩水分含量；Mad’为正向化的煤岩水分含量；μMad为待评价煤层正向化后煤岩水分含量的均值，σMad为待评价煤层正向化后煤岩水分含量的标准差。
7.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤7中，计算煤储层可压裂性评价分数Frac：Frac＝0.3BI+0.25Kc”+0.05Ed’+0.15EK+0.2σn’+0.05Mad” (20)式中，BI为归一化的脆性指数；Kc”为正向化、归一化的断裂韧性；Ed’为归一化的煤储层厚度；Ek’为归一化的围岩与煤储层弹性模量比值；σn’为归一化煤储层净压力系数；Mad”为正向化、归一化的水分含量；根据煤储层可压裂性评价分数Frac评价煤储层可压裂性。
8.根据权利要求1所述的煤储层可压裂性评价方法，其特征在于：在上述步骤8中，对可压裂性评价分数Frac进行分级。
一种煤储层可压裂性评价方法
技术领域
[0001] 本发明涉及煤层气开发技术领域，特别是涉及一种煤储层可压裂性评价方法。
背景技术
[0002] 对于煤层气开采过程中重要的一个环节就是水力压裂煤储层，因此在进行煤储层气开发有利区和甜点区优选时，煤储层可压裂性是一项重要评价参数。
[0003] 但目前在对于煤储层气的选区评价选定井位的过程中，评价方法主要集中在煤储层气的资源因素以及渗透率等煤储层的物性参数的优选，暂无合适方法或技术对煤储层可压裂性优劣进行评价。这不仅使得煤储层气选区及井位部署出现盲区，而且造成评价结果可信度降低。在工程技术施工的过程中，由于煤储层本身的物理性质的特殊性，对于压裂技术的要求也是十分苛刻的。
发明内容
[0004] 为了克服上述现有技术的不足，本发明提供了一种煤储层可压裂性评价方法。
[0005] 本发明所采用的技术方案是：一种煤储层可压裂性评价方法，包含以下步骤：步骤1：求取待评价煤储层脆性指数BI；
[0006] 步骤2：求取在围压下煤储层断裂韧性KC；
[0007] 步骤3：求取待评价煤层的煤层厚度Eh；
[0008] 步骤4：求取煤层与围岩弹性模量差异Ek’；
[0009] 步骤5：求取待评价煤层的净压力系数σn；
[0010] 步骤6：计算煤岩的水分含量Mad；
[0011] 步骤7：计算煤储层可压裂性评价分数Frac；
[0012] 步骤8：将煤层的可压裂性进行分级。
[0013] 进一步地，在上述步骤1中，获取煤储层弹性模量和泊松比，对弹性模量和泊松比进行归一化处理，求取待评价煤储层脆性指数BI：
[0014]
[0015]
[0016]
[0017] (1)式中，EBI′为待评价煤层归一化的弹性模量，EBImax和EBImin分别为煤储层最大和最小弹性模量，单位为GPa；(2)式中，μBI’为煤储层归一化的泊松比，μBImax和μBImin分别为煤储层最大和最小泊松比，无量纲；
[0018] 在上述步骤2中，步骤2.1，求取煤储层抗压强度σc，
[0019] σc＝(0.0045+0.0035Vsh)EKc (4)
[0020] 其中，Vsh是煤储层泥质含量，
[0021] Vsh＝(GR-GRmin)/(GRmax-GRmin) (5)
[0022] EKc是煤储层动态弹性模量，
[0023]
[0024] 步骤2.2，再求取煤储层单轴抗拉强度St，
[0025] St＝σc/A (7)
[0026] (4)式中，Vsh为待评价煤层泥质含量；(5)式中，GR为待评价煤层的伽马测井实测值，单位为API，GRmin为测得的井段砂岩层的自然伽马值，单位为API，GRmax为测得的井段泥岩层的自然伽马值，单位为API；(6)式中，EKc为待评价煤层的动态弹性模量，单位为GPa，vp为待评价煤层的纵波速度，单位为km/s，vs为待评价煤层的横波速度，单位为km/s；ρ为岩石密度，单位为g/cm3；(7)式中，A为常数；
[0027] 步骤2.3，在围压下岩石的断裂韧性KC，
[0028] KC＝0.0956pw+0.1383St-0.0820 (8)
[0029] 为简化计算，(8)式中围压pw用最小水平主应力代替；
[0030] 步骤2.4，对待评价煤层的断裂韧性进行正向化、归一化处理：
[0031] 对KC进行正向化、归一化：
[0032] Kc’＝(Kcij-Kcjmin)/(Kcjmax-Kcjmin) (9)
[0033] 式中，Kc’是正向化的断裂韧性，Kcij是待评价煤层的断裂韧性，Kcjmin是所有评价煤层的断裂韧性最小值，Kcjmax是所有待评价煤层断裂韧性的最大值；
[0034]
[0035] 式中，Kc”为正向化、归一化的断裂韧性；Kc’为正向化的断裂韧性；μKc为待评价煤层正向化后断裂韧性的均值，σKc为待评价煤层正向化后断裂韧性的标准差。
[0036] 进一步地，在上述步骤3中，步骤3.1，利用已有的煤田钻孔以及煤储层气参数井资料，采用内插法线性拟合约束计算研究区域的煤储层厚度，煤储层厚度为Eh，单位m。
[0037] 步骤3.2，对煤储层厚度进行归一化：
[0038]
[0039] 式中，Eh’为归一化的煤储层厚度，μEh为待评价煤层煤厚的均值，σEh为待评价煤层煤厚的标准差。
[0040] 进一步地，在上述步骤4中，步骤4.1，利用围岩与煤储层弹性模量的比值Ek表征煤储层与围岩的弹性模量差异：
[0041]
[0042] 式中，Ec为煤储层弹性模量，单位为GPa；Et为顶板弹性模量，单位为GPa；Eb为底板弹性模量，单位为GPa；
[0043] 步骤4.2，对弹性模量差异进行归一化：
[0044]
[0045] 式中，Ek’为归一化的弹性模量差异，μEk为待评价煤层弹性模量差异的均值，σEk为待评价煤层弹性模量差异标准差。
[0046] 进一步地，在上述步骤5中，步骤5.1，采用水压致裂法获得煤储层最小水平主应力Pc，单位为MPa。
[0047] 煤储层最大水平主应力σH：
[0048] σH＝3Pc-Pf-Po+T (14)
[0049] 式中，Pf为煤储层破裂压力，单位为MPa；Po为煤储层压力，单位为MPa；T为煤储层抗拉强度，单位为MPa。
[0050] 则煤储层净压力系数σn：
[0051]
[0052] 步骤5.2，对净压力系数进行归一化：
[0053]
[0054] 式中，σn’为归一化的弹性模量差异，μσn为待评价煤层净压力系数的均值，[0055] σσn为待评价煤层净压力系数的标准差。
[0056] 进一步地，在上述步骤6中，步骤6.1，通过对煤岩含水量与测井参数相关关系的分析，选择相关性最为密切的三个参数进行多元回归分析，依此建立了煤岩水分含量的预测方程：
[0057] Mad＝1.4655-0.5827×DEN-2.1115×GR+0.2319×ρs (17)
[0058] 式中：Mad为待评价煤层煤岩水分含量，单位为％；DEN为待评价煤层密度值，单位为g/cm3；GR为待评价煤层的自然伽马值，单位为API；ρs为待评价煤层的视电阻率值，单位为Ω·m。
[0059] 步骤6.2，对煤岩水分含量进行正向化、归一化：
[0060] Mad’＝(Madij-Madjmin)/(Madjmax-Madjmin) (18)
[0061] 式中，Mad’是正向化的煤岩水分含量，Madij是待评价煤层的煤岩水分含量，Madjmin是所有评价煤层的煤岩水分含量最小值，Madjmax是所有待评价煤层煤岩水分含量的最大值。
[0062]
[0063] 式中，Mad”为正向化、归一化的煤岩水分含量；Mad’为正向化的煤岩水分含量；μMad为待评价煤层正向化后煤岩水分含量的均值，σMad为待评价煤层正向化后煤岩水分含量的标准差。
[0064] 进一步地，在上述步骤6中，计算煤储层可压裂性评价分数Frac：
[0065] Frac＝0.3BI+0.25Kc”+0.05Ed’+0.15EK+0.2σn’+0.05Mad” (20)[0066] 式中，BI为归一化的脆性指数；Kc”为正向化、归一化的断裂韧性；Ed’为归一化的煤储层厚度；Ek’为归一化的围岩与煤储层弹性模量比值；σn’为归一化煤储层净压力系数；Mad”为正向化、归一化的水分含量；根据煤储层可压裂性评价分数Frac评价煤储层可压裂性。
[0067] 进一步地，在上述步骤8中，对可压裂性评价分数Frac进行分级：
[0068]
[0069]
[0070] 与现有技术相比，本发明的有益效果是：方法利用了研究区域的试井结果、测井实验参数、煤样的室内实验数据。对这些实验数据整合处理，筛选出煤储层的脆性指数、断裂韧性、煤储层厚度、煤储层与围岩的弹性模量差异、净压力系数、水分含量六个对煤储层可压裂性具有重大影响的参数，并对上述参数进行权重赋值，最后对研究区域的煤储层可压裂性进行评价优选。
具体实施方式
[0071] 为了加深对本发明的理解，下面结合实施例对本发明进一步说明，该实施例仅用于解释本发明，并不对本发明的保护范围构成限定。
[0072] 本发明所采用的技术方案是：一种煤储层可压裂性评价方法，包含以下步骤：
[0073] 步骤1：求取待评价煤储层脆性指数BI；
[0074] 步骤2：求取在围压下煤储层断裂韧性KC；
[0075] 步骤3：求取待评价煤层的煤层厚度Eh；
[0076] 步骤4：求取煤层与围岩弹性模量差异Ek’；
[0077] 步骤5：求取待评价煤层的净压力系数σn；
[0078] 步骤6：计算煤岩的水分含量Mad；
[0079] 步骤7：计算煤储层可压裂性评价分数Frac；
[0080] 步骤8：将煤层的可压裂性进行分级。
[0081] 在上述步骤1中，获取煤储层弹性模量和泊松比，对弹性模量和泊松比进行归一化处理，求取待评价煤储层脆性指数BI：
[0082]
[0083]
[0084]
[0085] (1)式中，EBI’为待评价煤层归一化的弹性模量，EBImax和EBImin分别为煤储层最大和最小弹性模量，单位为GPa；(2)式中，μBI’为煤储层归一化的泊松比，μBImax和μBImin分别为煤储层最大和最小泊松比，无量纲。
[0086] 进一步地，在上述步骤2中，步骤2.1，求取煤储层抗压强度σc，
[0087] σc＝(0.0045+0.0035Vsh)EKc (4)
[0088] 其中，Vsh是煤储层泥质含量，
[0089] Vsh＝(GR-GRmin)/(GRmax-GRmin) (5)
[0090] EKc是煤储层动态弹性模量，
[0091]
[0092] 步骤2.2，再求取煤储层单轴抗拉强度St，
[0093] St＝σc/A (7)
[0094] (4)式中，Vsh为待评价煤层泥质含量；(5)式中，GR为待评价煤层的伽马测井实测值，单位为API，GRmin为测得的井段砂岩层的自然伽马值，单位为API，GRmax为测得的井段泥岩层的自然伽马值，单位为API；(6)式中，EKc为待评价煤层的动态弹性模量，单位为GPa，vp为待评价煤层的纵波速度，单位为km/s，vs为待评价煤层的横波速度，单位为km/s；ρ为岩石密度，单位为g/cm3；(7)式中，A为常数30。
[0095] 步骤2.3，在围压下岩石的断裂韧性KC，
[0096] KC＝0.0956pw+0.1383St-0.0820 (8)
[0097] 为简化计算，(8)式中围压pw用最小水平主应力代替。
[0098] 步骤2.4，对待评价煤层的断裂韧性进行正向化、归一化处理：
[0099] 对KC进行正向化、归一化：
[0100] Kc’＝(Kcij-Kcjmin)/(Kcjmax-Kcjmin) (9)
[0101] 式中，Kc’是正向化的断裂韧性，Kcij是待评价煤层的断裂韧性，Kcjmin是所有评价煤层的断裂韧性最小值，Kcjmax是所有待评价煤层断裂韧性的最大值。
[0102]
[0103] 式中，Kc”为正向化、归一化的断裂韧性；Kc’为正向化的断裂韧性；μKc为待评价煤层正向化后断裂韧性的均值，σKc为待评价煤层正向化后断裂韧性的标准差。
[0104] 在上述步骤3中，步骤3.1，利用已有的煤田钻孔以及煤储层气参数井资料，采用内插法线性拟合约束计算研究区域的煤储层厚度，煤储层厚度为Eh，单位m。
[0105] 步骤3.2，对煤储层厚度进行归一化：
[0106]
[0107] 式中，Eh’为归一化的煤储层厚度，μEh为待评价煤层煤厚的均值，σEh为待评价煤层煤厚的标准差。
[0108] 在上述步骤4中，步骤4.1，利用围岩与煤储层弹性模量的比值Ek表征煤储层与围岩的弹性模量差异：
[0109]
[0110] 式中，Ec为煤储层弹性模量，单位为GPa；Et为顶板弹性模量，单位为GPa；Eb为底板弹性模量，单位为GPa；
[0111] 步骤4.2，对弹性模量差异进行归一化：
[0112]
[0113] 式中，Ek’为归一化的弹性模量差异，μEk为待评价煤层弹性模量差异的均值，σEk为待评价煤层弹性模量差异标准差。
[0114] 在上述步骤5中，步骤5.1，采用水压致裂法获得煤储层最小水平主应力Pc，单位为MPa。
[0115] 煤储层最大水平主应力σH：
[0116] σH＝3Pc-Pf-Po+T (14)
[0117] 式中，Pf为煤储层破裂压力，单位为MPa；Po为煤储层压力，单位为MPa；T为煤储层抗拉强度，单位为MPa。
[0118] 则煤储层净压力系数σn：
[0119]
[0120] 步骤5.2，对净压力系数进行归一化：
[0121]
[0122] 式中，σn’为归一化的弹性模量差异，μσn为待评价煤层净压力系数的均值，σσn为待评价煤层净压力系数的标准差。
[0123] 在上述步骤6中，步骤6.1，通过对煤岩含水量与测井参数相关关系的分析，选择相关性最为密切的三个参数进行多元回归分析，依此建立了煤岩水分含量的预测方程：
[0124] Mad＝1.4655-0.5827×DEN-2.1115×GR+0.2319×ρs (17)
[0125] 式中：Mad为待评价煤层煤岩水分含量，单位为％；DEN为待评价煤层密度值，单位为g/cm3；GR为待评价煤层的自然伽马值，单位为API；ρs为待评价煤层的视电阻率值，单位为Ω·m。
[0126] 步骤6.2，对煤岩水分含量进行正向化、归一化：
[0127] Mad’＝(Madij-Madjmin)/(Madjmax-Madjmin) (18)
[0128] 式中，Mad’是正向化的煤岩水分含量，Madij是待评价煤层的煤岩水分含量，Madjmin是所有评价煤层的煤岩水分含量最小值，Madjmax是所有待评价煤层煤岩水分含量的最大值。
[0129]
[0130] 式中，Mad”为正向化、归一化的煤岩水分含量；Mad’为正向化的煤岩水分含量；μMad为待评价煤层正向化后煤岩水分含量的均值，σMad为待评价煤层正向化后煤岩水分含量的标准差。
[0131] 在上述步骤7中，计算煤储层可压裂性评价分数Frac：
[0132] Frac＝0.3BI+0.25Kc”+0.05Ed’+0.15EK+0.2σn’+0.05Mad” (20)[0133] 式中，BI为归一化的脆性指数；Kc”为正向化、归一化的断裂韧性；Ed’为归一化的煤储层厚度；Ek’为归一化的围岩与煤储层弹性模量比值；σn’为归一化煤储层净压力系数；Mad”为正向化、归一化的水分含量；根据煤储层可压裂性评价分数Frac评价煤储层可压裂性。
[0134] 在上述步骤8中，对可压裂性评价分数Frac进行分级：
[0135]
[0136] 本发明的实施例公布的是较佳的实施例，但并不局限于此，本领域的普通技术人员，极易根据上述实施例，领会本发明的精神，并做出不同的引申和变化，但只要不脱离本发明的精神，都在本发明的保护范围内。
