一种水平应力差异系数测量方法及装置

摘要

本发明提供了一种水平应力差异系数测量方法及装置，包括：确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值；以待测量岩芯的底面中点为原点，以过中点的两条垂直线为横轴和纵轴，以过中点的轴线为竖轴，建立坐标系；确定各个监测位置与待测量岩芯的底面平行的切线相对于坐标系的横轴和纵轴的方向余弦，根据方向余弦确定应变值与待测量岩芯的水平应变分量之间的比例关系；根据上述应变值和比例关系，确定水平应变分量；根据水平应变分量，确定水平应力差异系数。本发明中，仅通过测量圆柱状岩芯的滞弹性恢复应变值，并根据已知的滞弹性应变恢复参数，就可以确定出泥页岩储层的水平应力差异系数，测量速度快，操作简单，成本低。

权利要求

1.一种水平应力差异系数测量方法，其特征在于，所述方法包括：分别采集各个预设的监测位置在多个预设时刻对应的应变值；

统计每个监测位置对应的应变值在多个预设时刻出现的次数；

将出现次数最多的所述应变值确定为所述监测位置对应的滞弹性恢复应变值；

以待测量岩芯的底面的中点作为原点，以过所述中点的两条垂直线作为横轴和纵轴，以过所述中点且垂直于所述横轴和所述纵轴的轴线作为竖轴，建立坐标系；

确定各个所述监测位置与所述待测量岩芯的底面平行的切线相对于所述坐标系的横轴和纵轴的方向余弦，根据所述方向余弦确定所述滞弹性恢复应变值与所述待测量岩芯的水平应变分量之间的比例关系；

根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

在每个监测位置上均设置一个应变计，通过每个监测位置上的应变计确定出每个监测位置对应的滞弹性恢复应变值，根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

根据所述水平应变分量，通过公式(1)得到未知参数λ的两个根λ1和λ2，将λ1和λ2中较大的一个根确定为待测量岩芯的水平最大应变值，将λ1和λ2中较小的一个根确定为待测量岩芯的水平最小应变值；

(ε'x-λ)(ε'y-λ)-ε'xyε'yx＝0       (1)其中，ε'x为X轴方向的应变分量，ε'y为Y轴方向的应变分量，ε'xy和ε'yx均为剪应变，且ε'xy和ε'yx大小相等，但是正负号不一样；

根据所述水平最大应变值和所述水平最小应变值确定待测量岩芯的滞弹性恢复平均应变值；

根据所述水平最大应变值和所述平均应变值，确定待测量岩芯的水平最大应力值；

根据所述水平最小应变值和所述平均应变值，确定待测量岩芯的水平最小应力值；

根据所述水平最大应力值和所述水平最小应力值的比值，确定所述水平应力差异系数，其中，计算所述水平最大应力值与所述水平最小应力值的比值与1之间的差值；

将所述差值确定为所述水平应力差异系数。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：根据所述水平最大应变值、矿物失水导致的应变值以及孔隙压力逸散导致的应变值确定水平实际最大应变值；

根据所述水平最小应变值、矿物失水导致的应变值以及孔隙压力逸散导致的应变值确定水平实际最小应变值；

根据所述水平实际最大应变值和所述水平实际最小应变值确定所述水平应力差异系数。

3.根据权利要求1所述的方法，其特征在于，各个所述监测位置在圆柱状的所述待测量岩芯曲面上均匀分布。

4.根据权利要求1所述的方法，其特征在于，相邻的所述监测位置之间间隔45°。

5.根据权利要求1所述的方法，其特征在于，所述监测位置的个数大于或等于3。

6.一种水平应力差异系数测量装置，其特征在于，所述装置包括：第一确定模块，用于分别采集各个预设的监测位置在多个预设时刻对应的应变值；

统计每个监测位置对应的应变值在多个预设时刻出现的次数；

将出现次数最多的所述应变值确定为所述监测位置对应的滞弹性恢复应变值；

建立模块，用于以待测量岩芯的底面的中点作为原点，以过所述中点的两条垂直线作为横轴和纵轴，以过所述中点且垂直于所述横轴和所述纵轴的轴线作为竖轴，建立坐标系；

第二确定模块，用于确定各个所述监测位置与所述待测量岩芯的底面平行的切线相对于所述坐标系的横轴和纵轴的方向余弦，根据所述方向余弦确定所述滞弹性恢复应变值与所述待测量岩芯的水平应变分量之间的比例关系；

第三确定模块，用于根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

第四确定模块，用于在每个监测位置上均设置一个应变计，通过每个监测位置上的应变计确定出每个监测位置对应的滞弹性恢复应变值，根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

根据所述水平应变分量，通过公式(1)得到未知参数λ的两个根λ1和λ2，将λ1和λ2中较大的一个根确定为待测量岩芯的水平最大应变值，将λ1和λ2中较小的一个根确定为待测量岩芯的水平最小应变值；

(ε'x-λ)(ε'y-λ)-ε'xyε'yx＝0       (1)其中，ε'x为X轴方向的应变分量，ε'y为Y轴方向的应变分量，ε'xy和ε'yx均为剪应变，且ε'xy和ε'yx大小相等，但是正负号不一样；

第五确定模块，用于根据所述水平最大应变值和所述水平最小应变值确定待测量岩芯的滞弹性恢复平均应变值；

根据所述水平最大应变值和所述平均应变值，确定待测量岩芯的水平最大应力值；

根据所述水平最小应变值和所述平均应变值，确定待测量岩芯的水平最小应力值；

根据所述水平最大应力值和所述水平最小应力值的比值，确定所述水平应力差异系数，其中，计算所述水平最大应力值与所述水平最小应力值的比值与1之间的差值；

将所述差值确定为所述水平应力差异系数。

说明书

一种水平应力差异系数测量方法及装置

技术领域

[0001] 本发明涉及非常规油气开发技术领域，具体而言，涉及一种水平应力差异系数测量方法及装置。

背景技术

[0002] 非常规油气开发技术是近十年来全球最大的能源技术革命，以水平井多级压裂和“井工厂模式”为核心的页岩油气勘探开发技术促使了全球页岩油气的成功开发。在页岩油气压裂开发中，水平应力差异系数是决定页岩储层能否体积压裂和实现经济开采的关键参数之一，水平应力差异系数指的是两相水平主应力差值与最小水平主应力的比值。

[0003] 现有技术中，大多是通过小型测试压裂和井壁破坏相结合的方法确定水平应力差异系数，即通过小型测试压裂获取最小水平主应力，通过孔壁破坏或岩石蠕变参数估算最大水平主应力的范围，进而确定水平应力差异系数，或者直接通过水压致裂地应力测量获取水平主应力信息。

[0004] 但是，通过小型测试压裂和井壁破坏相结合的方法确定水平应力差异系数，成本较高，且无法在压裂设计和评估前获取该系数；通过水压致裂地应力测量的方法，需要利用钻井平台进行，因此影响施工作业，且需要在裸眼段进行测试，容易发生事故，测试成本高，周期长。

发明内容

[0005] 有鉴于此，本发明实施例的目的在于提供一种水平应力差异系数测量方法及装置，以试图解决或者缓解上述问题。

[0006] 第一方面，本发明实施例提供了一种水平应力差异系数测量方法，所述方法包括：

[0007] 确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值；

[0008] 以所述待测量岩芯的底面的中点作为原点，以过所述中点的两条垂直线作为横轴和纵轴，以过所述中心的轴线作为竖轴，建立坐标系；

[0009] 确定各个所述监测位置与所述待测量岩芯的底面平行的切线相对于所述坐标系的横轴和纵轴的方向余弦，根据所述方向余弦确定所述滞弹性恢复应变值与所述待测量岩芯的水平应变分量之间的比例关系；

[0010] 根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

[0011] 根据所述水平应变分量，确定所述待测量岩芯的水平最大应变值和水平最小应变值；

[0012] 根据所述水平最大应变值和所述水平最小应变值确定所述待测量岩芯的水平应力差异系数。

[0013] 结合第一方面，本发明实施例提供了上述第一方面的第一种可能的实现方式，其中，所述根据所述水平最大应变值和所述水平最小应变值确定所述待测量岩芯的水平应力差异系数，包括：

[0014] 根据所述水平最大应变值和所述水平最小应变值确定待测量岩芯的滞弹性恢复平均应变值；

[0015] 根据所述水平最大应变值和所述平均应变值，确定待测量岩芯的水平最大应力值；

[0016] 根据所述水平最小应变值和所述平均应变值，确定待测量岩芯的水平最小应力值；

[0017] 根据所述水平最大应力值和所述水平最小应力值的比值，确定所述水平应力差异系数。

[0018] 结合第一方面的第一种可能的实现方式，本发明实施例提供了上述第一方面的第二种可能的实现方式，其中，所述根据所述水平最大应力值和所述水平最小应力值的比值，确定所述水平应力差异系数，包括：

[0019] 计算所述水平最大应力值与所述水平最小应力值的比值与1之间的差值；

[0020] 将所述差值确定为所述水平应力差异系数。

[0021] 结合第一方面，本发明实施例提供了上述第一方面的第三种可能的实现方式，其中，所述方法还包括：

[0022] 根据所述水平最大应变值、矿物失水导致的应变值以及孔隙压力逸散导致的应变值确定水平实际最大应变值；

[0023] 根据所述水平最小应变值、矿物失水导致的应变值以及空隙压力逸散导致的应变值确定水平实际最小应变值；

[0024] 根据所述水平实际最大应变值和所述水平实际最小应变值确定所述水平应力差异系数。

[0025] 结合第一方面，本发明实施例提供了上述第一方面的第四种可能的实现方式，其中，各个所述监测位置在圆柱状的所述待测量岩芯曲面上同一高度处均匀分布。

[0026] 结合第一方面，本发明实施例提供了上述第一方面的第五种可能的实现方式，其中，相邻的所述监测位置之间间隔45°。

[0027] 结合第一方面，本发明实施例提供了上述第一方面的第六种可能的实现方式，其中，所述监测位置的个数大于或等于3。

[0028] 结合第一方面，本发明实施例提供了上述第一方面的第七种可能的实现方式，其中，所述确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值，包括：

[0029] 分别采集各个所述监测位置在多个预设时刻对应的应变值；

[0030] 统计每个监测位置对应的应变值在多个预设时刻出现的次数；

[0031] 将出现次数最多的所述应变值确定为所述监测位置对应的滞弹性恢复应变值。

[0032] 结合第一方面的第七种可能的实现方式，本发明实施例提供了上述第一方面的第八种可能的实现方式，其中，各个所述监测位置处设置有应变计，通过所述应变计采集各个所述监测位置在多个预设时刻对应的应变值。

[0033] 第二方面，本发明实施例提供了一种水平应力差异系数测量装置，其中，所述装置包括：

[0034] 第一确定模块，用于确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值；

[0035] 建立模块，用于以所述待测量岩芯的底面的中点作为原点，以过所述中点的两条垂直线作为横轴和纵轴，以过所述中线的轴线作为竖轴，建立坐标系；

[0036] 第二确定模块，用于确定各个所述监测位置与所述待测量岩芯的底面平行的切线相对于所述坐标系的横轴和纵轴的方向余弦，根据所述方向余弦确定所述滞弹性恢复应变值与所述待测量岩芯的水平应变分量之间的比例关系；

[0037] 第三确定模块，用于根据确定的所述滞弹性恢复应变值和所述比例关系，确定所述待测量岩芯的水平应变分量；

[0038] 第四确定模块，用于根据所述水平应变分量，确定所述待测量岩芯的水平最大应变值和水平最小应变值；

[0039] 第五确定模块，用于根据所述水平最大应变值和所述水平最小应变值确定所述待测量岩芯的水平应力差异系数。

[0040] 本发明实施例提供的水平应力差异系数测量方法及装置，仅通过测量圆柱状岩芯的滞弹性恢复应变值，并根据已知的滞弹性应变恢复参数，就可以确定出泥页岩储层的水平应力差异系数，测量速度快，操作简单，成本低。

[0041] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。

附图说明

[0042] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

[0043] 图1示出了本发明实施例1所提供的水平应力差异系数测量方法的流程图；

[0044] 图2示出了本发明实施例1所提供的水平应力差异系数测量方法中确定各个监测位置对应的滞弹性恢复应变值的流程图；

[0045] 图3示出了本发明实施例1所提供的水平应力差异系数测量方法中的待测量岩芯的示意图；

[0046] 图4示出了本发明实施例2所提供的水平应力差异系数测量装置的结构示意图。

具体实施方式

[0047] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0048] 考虑到现有技术中，大都是通过小型测试压裂和井壁破坏相结合的方法测量水平应力差异系数，或者直接通过水压致裂地应力测量获取水平主应力信息。但是，如果通过小型测试压裂和井壁破坏相结合的方法确定水平应力差异系数，成本较高，且无法在压裂设计和评估前获取该系数；通过水压致裂地应力测量方法，需要利用钻井平台进行，因此影响施工作业，且需要在裸眼段进行测试，容易发生事故，测试成本高，周期长。基于此，本发明实施例提供了一种水平应力差异系数测量方法和装置，下面通过实施例进行描述。

[0049] 实施例1

[0050] 本发明实施例提供了一种水平应力差异系数测量方法，如图1所示，该方法包括步骤S110-S160,具体如下。

[0051] S110，确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值。

[0052] 在采用本发明实施例提供的方法测量泥页岩的水平应力差异系数之前，需要现场钻取泥页岩岩芯，且钻取的泥页岩岩芯为圆柱形，然后将钻取的泥页岩岩芯的表面进行清洗，得到待测量岩芯。

[0053] 将泥页岩岩芯钻取出来后，这时，泥页岩岩芯不再受力，即处于卸荷状态，本发明实施例，正是根据待测量岩芯卸荷后的滞弹性恢复变形与水平应力之间的对应关系，来确定泥页岩储层的水平应力差异系数。

[0054] 其中，上述多个监测位置可以位于曲面上同一高度处，也可以位于曲面上不同高度处，本发明实施例并不对上述多个监测位置在曲面上的具体高度进行限制。

[0055] 其中，上述各个监测位置在圆柱状的待测量岩芯曲面上均匀分布，即每个监测位置和相邻的监测位置之间的角度间隔相等，比如说，如果，在待测量岩芯上设置三个监测位置，则可以是每个监测位置和其相邻的监测位置之间间隔120°。

[0056] 当然，上述只是以三个监测位置为例具体说明监测位置之间如何均匀分布，并没有限定设置的监测位置的具体个数，另外，上述监测位置之间也可以不均匀分布。

[0057] 为了更方便的计算待测量岩芯的水平应力差异系数，相邻的监测位置之间间隔45°。

[0058] 如果相邻的监测位置之间间隔45°，在圆柱状待测量岩芯曲面上最多可以设置八个监测位置，因此，在测量待测量岩芯的水平应力差异系数时，可以确定八个监测位置的应变值，确定的监测位置的应变值的个数也可以小于八，只要能够测量出待测量岩芯的水平应力差异系数即可，本发明实施例，并不对上述使用的监测位置的个数进行限定。

[0059] 当然，为了能够确定出待测量岩芯的水平应力差异系数，需要确定待测量岩芯在水平多个方向上的滞弹性恢复应变值，即上述监测位置的个数大于或等于3，即在待测量岩芯曲面上至少设置三个监测位置。

[0060] 由于待测量岩芯卸荷后水平方向的滞弹性恢复应变值在一段时间内会随着时间而变化，一段时间过后，应变值趋于稳定，这时将稳定的应变值确定为待测量岩芯在该方向的滞弹性恢复应变值，如图2所示，确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值，包括步骤S210-S230，具体如下。

[0061] S210，分别采集各个监测位置在多个预设时刻对应的应变值。

[0062] S220，统计每个监测位置对应的应变值在多个预设时刻出现的次数。

[0063] S230，将出现次数最多的应变值确定为监测位置对应的滞弹性恢复应变值。

[0064] 其中，上述在采集各个监测位置的应变值时，需要保持待测量岩芯的温度恒定，这样能够消除待测量岩芯温度变化带来的误差。

[0065] 其中，在本发明实施例中，可以测量一段时间内各个监测位置对应的应变值，并且测量的这一段时间应该大于待测量岩芯的滞弹性应变值趋于稳定的时间，将测量的这一段时间段预先划分为多个时刻，采集各个监测位置在预设时刻对应的应变值，在开始一段时间内，各个监测位置对应的应变值会随着时间发生变化，当一段时间过后，各个监测位置对应的应变值趋于稳定，即各个监测位置对应的应变值不再随着时间发生变化，因此，会在连续多个预设时刻测量的该监测位置的应变值是相等的，即该监测位置对应的应变值会在多个预设时刻出现，因此，可以统计每个监测位置对应的应变值在多个预设时刻出现的次数，将出现次数最多的应变值确定为该监测位置对应的滞弹性恢复应变值，也就是，出现次数最多的应变值就是该监测位置稳定的应变值。

[0066] 其中，为了采集各个监测位置的应变值，上述各个监测位置处设置有应变计，通过应变计采集各个监测位置在多个预设时刻对应的应变值。

[0067] 其中，上述应变计贴在待测量岩芯表面的监测位置处，在上述各个监测位置处可以设置一个应变计，也可以设置一组应变计，一组应变计中可以包括多个应变计，如果，在每个监测位置上均设置一组应变计，则在每个监测位置对应的应变计组中包括的应变计的个数相等。

[0068] 如果在每个监测位置上设置一组应变计，则将该组应变计内各个应变计测量的应变值的平均值作为该监测位置的应变值。

[0069] S120，以上述待测量岩芯的底面的中心点作为原点，以过该中点的两条垂直线作为横轴和纵轴，以过上述中点的轴线作为竖轴，建立坐标系。

[0070] 其中，上述可以以圆柱形待测量岩芯的上底面的中心点作为原点，也可以以圆柱形待测量岩芯的下底面的中心点作为原点，在本发明实施例中，以待测量岩芯的上底面的中心点作为原点。

[0071] 在本发明实施例中，将建立的坐标系的原点用O表示，横轴用X表示，纵轴用Y表示，竖轴用Z表示，建立待测量岩芯的O-XYZ坐标系。

[0072] S130，确定各个监测位置与待测量岩芯的底面平行的切线相对于坐标系的横轴和纵轴的方向余弦，根据该方向余弦确定上述滞弹性恢复应变值与待测量岩芯的水平应变分量之间的比例关系。

[0073] 下面将以在每个监测位置设置一个应变计、总共设置四个监测位置，且相邻每个监测位置之间间隔45°为例，详细介绍如何确定上述应变值与待测量岩芯的水平应变分量之间的比例关系。

[0074] 如图3所示，上述四个监测位置分别记为第一监测位置a1、第二监测位置a2、第三监测位置a3和第四监测位置a4，在每个监测位置上均设置一个应变计，第一监测位置a1上的两个应变计为C1，且将应变计C1测量得到的应变值确定为第一监测位置a1对应的应变值，第二监测位置a2上的两个应变计为C3，且将应变计C3测量得到的应变值确定为第二监测位置a2对应的应变值，第三监测位置a3上的两个应变计为C2，且将应变计C2测量得到的应变值确定为第三监测位置a3对应的应变值，第四监测位置a4上的两个应变计为C4，且将应变计C4测量得到的应变值确定为第四监测位置a4对应的应变值。

[0075] 分别计算第一监测位置a1、第二监测位置a2、第三监测位置a3和第四监测位置a4与待测量岩芯的底面平行的切线相对于上述坐标系的横轴和纵轴的方向余弦，该方向余弦指的是上述切线分别与横轴和纵轴的夹角的余弦值，如图3所示，第一监测位置a1与待测量岩芯的底面平行的切线与X轴平行，因此，该切线与X轴的夹角为0°，与Y轴的夹角为90°，由此，可以计算出第一监测位置a1分别与X轴和Y轴的方向余弦分别为1和0,，第二监测位置a2与待测量岩芯的底面平行的切线与Y轴平行，因此，该切线与X轴的夹角为90°，与Y轴的夹角为0°，由此，可以计算出第二监测位置分别与X轴和Y轴的方向余弦分别为0和1，第三监测位置a3与待测量岩芯的底面平行的切线与X轴的夹角为45°，与Y轴的夹角为135°，由此，可以计算出第三监测位置a3分别与X轴和Y轴的方向余弦分别为0.707和-0.707，第四监测位置a4与待测量岩芯的底面平行的切线与X轴的夹角为135°，与Y轴的夹角为45°，由此，可以计算出第四监测位置a4分别与X轴和Y轴的方向余弦分别为-0.707和0.707，上述第一监测位置a1、第二监测位置a2、第三监测位置a3和第四监测位置a4分别与X轴和Y轴的方向余弦如表1所示。

[0076] 表1

[0077]

[0078] 其中，在上表中，监测位置与待测量岩芯的底面平行的切线相对于X轴的方向余弦用li表示，监测位置待测量岩芯的底面平行的切线相对于Y轴的方向余弦用mi表示，其中，i＝1,2,3,4。当i取值为1时，l1表示第一监测位置与待测量岩芯的底面平行的切线相对于X轴的方向余弦，m1表示第一监测位置与待测量岩芯的底面平行的切线相对于Y轴的方向余弦。

[0079] 当确定出每个监测位置与待测量岩芯的底面平行的切线相对于X轴和Y轴的方向余弦后，根据如下公式确定出每个监测位置对应的应变值与待测量岩芯的水平应变分量之间比例系数，

[0080]

[0081] 通过上述公式确定监测位置对应的滞弹性恢复应变值与待测量岩芯的水平分量之间的比例系数后，得到，监测位置对应的应变值与待测量岩芯的水平分量之间的比例关系为：b＝Aε′，其中，在该公式中A为上述计算得出的比例系数，b为每个监测位置对应的滞弹性恢复应变值，ε′为待测量岩芯的水平应变分量，且ε′＝[ε′x,ε′y,ε′xy]T，其中，ε′x为X轴方向的应变分量，ε′y为Y轴方向的应变分量，ε′xy为剪应变。

[0082] S140，根据确定的上述滞弹性恢复应变值和上述比例关系，确定待测量岩芯的水平应变分量。

[0083] 通过每个监测位置上的应变计可以确定出每个监测位置对应的滞弹性恢复应变值，如果在每个监测位置上设置一个应变计，则将该应变计采集的稳定的应变值确定为该监测位置对应的应变值，如果在每个监测位置上设置有一组应变计，则将该组内所有应变计采集的稳定的应变值的平均值确定为该监测位置对应的应变值。

[0084] 第一监测位置对应的应变值用b1表示，第二监测位置对应的应变值用b2表示，第三监测位置对应的应变值用b3表示，第四监测位置对应的应变值用b4表示，则b＝[b1,b2,b3,b4]T为待测量岩芯的应变值。

[0085] 根据公式b＝Aε′以及确定待测量岩芯的应变值b＝[b1,b2,b3,b4]T，可以计算得出水平应变分量，具体计算过程如下：

[0086] 在公式b＝Aε′的两边各乘以应变值与水平应变分量之间的比例系数的转置矩阵，即：

[0087] AT Aε′＝AT b

[0088] 通过解上述公式，可以得出ε′＝(AT A)-1AT b。

[0089] S150，根据上述水平应变分量，确定待测量岩芯水平最大应变值和水平最小应变值。

[0090] 通过上述方法确定出了待测量岩芯的水平应变分量，接下来介绍如何根据水平应变分量确定待测量岩芯水平最大应变值和水平最小应变值。

[0091] 水平面上的应变状态可以通过两个主应变值来表示，这两个主应变就是水平最大应变值和水平最小应变值，通过水平应变分量计算水平最大应变值和水平最小应变值可以通过求解如下方程组得到：

[0092]

[0093] 其中，在上述公式中，λ为待求解的未知参数，li为监测位置与待测量岩芯的底面平行的切线相对于X轴的方向余弦，mi为监测位置与待测量岩芯的底面平行的切线相对于Y轴的方向余弦，ε′x为X轴方向的应变分量，ε′y为Y轴方向的应变分量，ε′xy和ε′yx均为剪应变，且ε′xy和ε′yx大小相等，但是正负号不一样。

[0094] 上述方程组为齐次线性方程组，上述方程组非零解的必要和充分条件时系数行列式为零，即：

[0095]

[0096] 将上述行列式展开后为一元二次方程式，即：

[0097] (ε′x-λ)(ε′y-λ)-ε′xyε′yx＝0

[0098] 通过求解上述公式，可以得到关于未知参数λ的两个根λ1和λ2，将λ1和λ2中较大的一个根确定为待测量岩芯的水平最大应变值，将λ1和λ2中较小的一个根确定为待测量岩芯的水平最小应变值。

[0099] S160，根据水平最大应变值和水平最小应变值确定水平应力差异系数。

[0100] 其中，上述根据水平最大应变值和水平最小应变值确定待测量岩芯的水平应力差异系数，包括：

[0101] 根据水平最大应变值和水平最小应变值确定待测量岩芯的滞弹性恢复平均应变值；根据水平最大应变值和平均应变值，确定待测量岩芯的水平最大应力值；根据水平最小应变值和平均应变值，确定待测量岩芯的水平最小应力值；根据水平最大应力值和水平最小应力值之间的比值，确定水平应力差异系数。

[0102] 在本发明实施例中，当确定了水平最大应变值和水平最小应变值之后，根据水平最大应变值和水平最小应变值，通过如下公式确定待测量岩芯的平均应变值；

[0103] ε平＝(εmax+εmin)/2

[0104] 其中，在上述公式中，ε平为待测量岩芯的滞弹性恢复平均应变值，εmax为待测量岩芯的水平最大应变值，εmin为待测量岩芯的水平最小应变值。

[0105] 根据水平应变值和平均应变值，通过如下公式，确定出水平最大应力值和水平最小应力值；

[0106] σi＝ei/Jas+em/Jav+P0

[0107] 其中，在上述公式中，i的取值为H和h，当i为H时表示水平最大应力值，当i为h时表示水平最小应力值，ei为滞弹性恢复偏应变，即水平应变值与平均应变值的差值，平均应变值为水平最大应变值和水平最小应变值的平均值，当i为H时，ei为水平最大应变值与平均应变值的差值，当i为h时，ei为水平最大应变值与平均应变值的差值，em为滞弹性恢复体应变，即滞弹性恢复平均应变值ε平，Jas为滞弹性恢复应变柔度，Jav为体积非滞弹性应变柔度，P0为孔隙压力，σi为水平应力值。

[0108] 具体的，根据水平最大应变值和平均应变值，通过如下公式，确定待测量岩芯的水平最大应力值为：

[0109] σH＝(εmax-ε平)/Jas+ε平/Jav+P0

[0110] 根据水平最小应变值和平均应变值，通过如下公式，确定待测量岩芯的水平最小应力值为：

[0111] σh＝(εmin-ε平)/Jas+ε平/Jav+P0

[0112] 通过上述公式，可以确定出水平最大应力值与水平最小应力值之间的比值为：

[0113]

[0114] 由于Jav/Jas的取值一定，在某些具体实施例中Jas/Jav的取值为1.56，因此，水平最大应力值与水平最小应力值之间的比值可以确定为：

[0115]

[0116] 将已经确定出的水平最大应变值与水平最小应变值之间的比值代入上述公式，可以确定出水平最大应力值与水平最小应力值之间的比值，将水平最大应变值与水平最小应变值之间的比值记为εmax/εmin＝K，因此，可以得出：

[0117]

[0118] 在上式中，P0为已知的，Jav及εmin均为已知的，因此，通过上述可以确定出水平最大应力值与水平最小应力值之间的比值。

[0119] 另外，从上述公式中也可以看出，水平最大应力值与水平最小应力值之间的比值随着K值的增加而增加。当确定了待测量岩芯的水平最大应力值与水平最小应力值之间比值后，根据水平最大应力值和水平最小应力值之间的比值，确定待测量岩芯的水平应力差异系数，具体包括：

[0120] 计算水平最大应力值与水平最小应力值的比值与1之间的差值；将上述差值确定为水平应力差异系数。

[0121] 在本发明实施例中，水平应力差异系数为水平最大应力值和水平最小应力值之间的差值，与水平最小应力值之间的比值，即DHSR＝(σH-σh)/σh，其中，在上述公式中，DHSR为水平应力差异系数，σH为水平最大应力值，σh为水平最小应力值，将上述公式进行变形可以得到， 因此，将水平最大应力值与水平最小应力值之间的比值与1之间的差值确定为水平应力差异系数。

[0122] 在测量泥页岩储层的水平应力差异系数时，因为矿物失水和孔隙压力逸散也会带来应变，在计算水平应力差异系数时，还需要考虑矿物失水导致的应变值以及孔隙压力逸散导致的应变值，本发明实施例提供的方法还包括：

[0123] 根据水平最大应变值、矿物失水导致的应变值以及孔隙压力逸散导致的应变值确定水平实际最大应变值；根据水平最小应变值、矿物失水导致的应变值以及孔隙压力逸散导致的应变值确定水平实际最小应变值；根据水平实际最大应变值和水平实际最小应变值确定水平应力差异系数。

[0124] 在步骤S150中，已经确定出水平最大应变值为εmax，水平最小应变值为εmin，将上述矿物失水导致的应变值记为εkw，将上述孔隙压力逸散导致的应变值记为εkx，因此，可以通过公式εmax＝εASR max+εkw+εkx计算出水平实际最大应变值，其中，上述公式中εASRmax表示水平实际最大应变值，即待测量岩芯因原地应力释放引起的最大应变值，根据公式εmin＝εASR min+εkw+εkx计算出水平实际最小应变值，其中，εASR min表示水平实际最小应变值，即待测量岩芯因原地应力释放引起的最小应变值。

[0125] 当确定出了待测量岩芯的水平实际最大应变值εASRmax和水平实际最小应变值εASR min后，计算水平实际最大应变值εASRmax和水平实际最小应变值εASR min的平均值εASR平，其中εASR平＝(εASR max+εASR min)/2。

[0126] 接下来，根据公式 计算水平实际最大应力值和水平实际最小应力值的比值。

[0127] 在上述公式中，σH实为水平实际最大应力值，σh实为水平实际最小应力值，P0为孔隙压力，Jav为体积非滞弹性应变柔度。

[0128] εASRmax、εASR min和εASR平、P0和Jav均为已知的，因此，通过上述公式可以确定出水平实际最大应力值和水平实际最小应力值之间的比值。

[0129] 当确定了水平实际最大应力值和水平实际最小应力值之间的比值后，根据确定待测量岩芯的水平应力差异系数。

[0130] 本发明实施例提供的水平应力差异系数测量方法，仅通过测量圆柱状岩芯的滞弹性恢复应变值，并根据已知的滞弹性应变恢复参数，就可以确定出泥页岩储层的水平应力差异系数，测量速度快，操作简单，成本低。

[0131] 实施例2

[0132] 本发明实施例提供了一种水平应力差异系数测量装置，如图4所示，该装置包括第一确定模块410、建立模块420、第二确定模块430、第三确定模块440、第四确定模块450和第五确定模块460；

[0133] 上述第一确定模块410，用于确定圆柱状的待测量岩芯曲面上多个监测位置的滞弹性恢复应变值；

[0134] 上述建立模块420，用于以待测量岩芯的底面的中点作为原点，以过该中点的两条垂直线作为横轴和纵轴，以过该中点的轴线作为竖轴，建立坐标系；

[0135] 上述第二确定模块430，用于确定各个监测位置与待测量岩芯的底面平行的切线相对于上述坐标系的横轴和纵轴的方向余弦，根据方向余弦确定上述滞弹性恢复应变值与待测量岩芯的水平分量之间的比例关系；

[0136] 上述第三确定模块440，用于根据确定的上述滞弹性恢复应变值和比例关系，确定待测量岩芯的水平应变分量；

[0137] 上述第四确定模块450，用于根据上述水平应变分量，确定待测量岩芯的水平最大应变值和水平最小应变值；

[0138] 上述第五确定模块460，用于根据上述水平最大应变值和水平最小应变值确定待测量岩芯的水平应力差异系数。

[0139] 本发明实施例提供的水平应力差异系数测量装置，仅通过测量圆柱状岩芯的滞弹性恢复应变值，并根据已知的滞弹性应变恢复参数，就可以确定出泥页岩储层的水平应力差异系数，测量速度快，操作简单，成本低。

[0140] 本发明实施例所提供的水平应力差异系数测量装置可以为设备上的特定硬件或者安装于设备上的软件或固件等。本发明实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，前述描述的系统、装置和单元的具体工作过程，均可以参考上述方法实施例中的对应过程，在此不再赘述。

[0141] 在本发明所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

[0142] 所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

[0143] 另外，在本发明提供的实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。

[0144] 所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

[0145] 应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释，此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

[0146] 最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。