CN110470581B - 确定储层应力敏感程度的方法、装置及存储介质 - Google Patents

Abstract

本说明书提供了一种确定储层应力敏感程度的方法、装置及存储介质，该方法包括：获取储层的有效应力范围；在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。本说明书可以实现对储层应力敏感程度的定量评价。

Description

确定储层应力敏感程度的方法、装置及存储介质

技术领域

本说明书涉及油气藏开发技术领域，尤其是涉及一种确定储层应力敏感程度的方法、装置及存储介质。

背景技术

在油气藏钻井、完井和生产过程中，储层上覆岩层压力保持不变，储层中的流体不断地产出，孔隙压力降低，作用在岩心骨架上的有效应力不断增大，使得岩石发生变形。岩石的形变会引起岩石孔隙结构和孔隙体积的变化，例如孔喉的变形、裂缝的闭合等，导致储层渗透率和孔隙度发生变化，该现象称之为应力敏感。由于储层渗透率和孔隙度受应力状态的影响而不为常数，从而影响储层流体向井筒中的流动，给合理高效地开发油气藏增加了难度。

现有关于储层渗透率和孔隙度应力敏感性研究如下：

对储层应力敏感性测量方法的研究，主要是通过研发改进常规应力敏感测试仪器和实验流程，尽可能地模拟储层的实际条件，结合实验数据建立与储层特征(岩性、粘土矿物含量、孔隙结构、裂缝倾角、胶结程度和胶结物类型等)相匹配的应力敏感模型；或者在测试中引入声发射、核磁共振和CT扫描技术，构建与孔隙度和渗透率的相关性模型，缩短了致密储层岩心渗透率和孔隙度的测量时间，提高测量精度。但以上研究缺少具体储层应力敏感程度评价指标。

对储层应力敏感性评价方法的研究主要有行业标准和应力敏感系数法。中华人民共和国石油天然气行业标准SY/T 5358—2010《储层敏感性流动实验评价方法》给出了应力敏感性评价实验方法和评价指标，该方法通过计算最大岩心渗透率损害率和不可逆渗透率损害率并结合应力敏感性损害程度评价指标，确定损害程度，但是最大有效应力点的选择将会影响渗透率损害率最终数值。应力敏感系数法是对有效应力和渗透率进行无因次化，分别求以10为底的对数和1/3次幂，将其拟合成截距为1的直线，直线的斜率的相反数即为应力敏感系数，可以看出一块岩样只对应一个应力敏感系数，但是提出的评价指标仍然只能反映储层应力敏感性的强弱，没能给出具体的有效应力警戒值或临界值。

综上所述，目前关于测试储层岩石渗透率的应力敏感性的研究多是通过研发改进实验仪器和测试方法来提高渗透率的准确性；或是提出评价储层应力敏感性的指标，但仅能笼统地表明储层应力敏感性的强弱，并不能直接有效地指导油气生产。因此迫切需要可以定量地评价储层应力敏感程度的方法。

发明内容

本说明书实施例的目的在于提供一种确定储层应力敏感程度的方法、装置及存储介质，以实现对储层应力敏感程度的定量评价。

为达到上述目的，一方面，本说明书实施例提供了一种确定储层应力敏感程度的方法，包括：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

其中，所述根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度，包括：

对所述第一关系曲线进行求导，获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力用于表征所述岩心样品的应力加载敏感程度；

对所述第二关系曲线进行求导，获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力用于表征所述岩心样品的应力卸载敏感程度。

其中，还包括：

确定所述第一临界有效应力与所述第二临界有效应力的差值；

根据所述差值确定所述岩心样品的渗透率的滞后程度。

其中，所述第一关系曲线包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

其中，所述获取储层的有效应力范围，包括：

获取所述储层的上覆岩层压力、初始地层孔隙流体压力和实际生产压力；

根据所述上覆岩层压力和所述初始地层孔隙流体压力确定所述储层的初始有效应力；

根据所述上覆岩层压力和所述实际生产压力确定所述储层的最大有效应力；

根据所述初始有效应力和所述最大有效应力确定所述储层的有效应力范围。

其中，所述应力加载序列包括：由多个等间隔的应力加载点形成的应力递增序列；所述应力卸载序列包括：由多个等间隔的应力卸载点形成的应力递减序列；所述应力递增序列及所述应力递减序列中的最小值为所述有效应力范围的最小值，所述应力递增序列及所述应力递减序列中的最大值为所述有效应力范围的最大值。

其中，所述应力递增序列中的应力加载点与所述应力递减序列中的应力卸载点对应。

另一方面，本说明书实施例还提供了一种确定储层应力敏感程度的装置，包括：

应力范围获取模块，用于获取储层的有效应力范围；

渗透率确定模块，用于在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

关系曲线确定模块，用于确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

敏感程度确定模块，用于根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

另一方面，本说明书实施例还提供了另一种确定储层应力敏感程度的装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

另一方面，本说明书实施例还提供了一种计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

由以上本说明书实施例提供的技术方案可见，通过本说明书实施例可以定量获得岩心样品的渗透率对不同有效应力范围的敏感程度，从而可以作为开采油/气有效应力的警戒值，以用于指导实际生产。

附图说明

为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本说明书一实施例中确定储层应力敏感程度的方法的流程图；

图2为本说明书一实施例中所示的归一化渗透率随归一化有效应力加载/卸载变化的拟合曲线示意图；

图3为本说明书一实施例中所示的归一化渗透率随归一化有效应力加载/卸载变化的拟合曲线在求导后的示意图；

图4为本说明书一实施例中确定储层应力敏感程度的装置的结构框图；

图5为本说明书另一实施例中确定储层应力敏感程度的装置的结构框图。

具体实施方式

为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。

参考图1所示，本说明书一实施例的确定储层应力敏感程度的方法，可以包括以下步骤：

S101、获取储层的有效应力范围。

本说明书一实施例中，所述获取储层的有效应力范围例如可以包括：获取所述储层的上覆岩层压力、初始地层孔隙流体压力和实际生产压力；根据所述上覆岩层压力和所述初始地层孔隙流体压力确定所述储层的初始有效应力；根据所述上覆岩层压力和所述实际生产压力确定所述储层的最大有效应力；根据所述初始有效应力和所述最大有效应力确定所述储层的有效应力范围。

例如，可以根据油气藏的实际地质和生产资料，获取储层的上覆岩层压力、初始地层孔隙流体压力和实际生产压力(本文以定压力生产为例)。其中，有效应力的定义为：σ_eff＝σ-αp，其中σ_eff为有效应力，σ为上覆岩层压力，α为毕奥特(Biot)系数，p为地层孔隙流体压力。当p为初始地层孔隙流体压力时，该式求得的是初始有效应力；而将上覆岩层压力减去实际生产压力(小于初始地层孔隙流体压力)可以获得储层的最大有效应力。在其他实施例中，上述初始有效应力和最大有效应力，可以根据油气藏实际情况进行选择，也可以根据实验研究的需要和实验设备的性能来确定。

S102、在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率。

本说明书一实施例中，为获得储层的应力敏感程度，需获取可以代表储层的岩心样品。其中，储层的岩心样品制备要充分考虑油气藏的实际特征，例如岩性、粘土矿物含量、孔隙结构、裂缝倾角、胶结程度和胶结物类型等，尽可能地钻取井下岩心，若没有井下岩心，也可采用与井下岩心矿物成分相同的露头岩心或人造岩心，但获得的应力敏感临界值可能因岩心的限制而与实际油气藏有一定的偏差。

本说明书一实施例中，所述应力加载序列例如可以包括：由多个等间隔的应力加载点形成的应力递增序列；所述应力卸载序列例如可以包括：由多个等间隔的应力卸载点形成的应力递减序列；所述应力递增序列及所述应力递减序列中的最小值为所述有效应力范围的最小值，所述应力递增序列及所述应力递减序列中的最大值为所述有效应力范围的最大值。此外，为了便于后续分析岩心样品的渗透率的滞后性，所述应力递增序列中的应力加载点与所述应力递减序列中的应力卸载点对应(即应力加载点的数量、大小与应力卸载点的数量、大小相同)，一般地，应力递增序列中的应力加载点的数量以及应力递减序列中的应力卸载点的数量可以为5个或更多个。

本说明书一实施例中，岩心样品在各应力加载点及应力卸载点处的渗透率可以通过渗透率测量仪等设备测量得到。并且，需确保岩心样品在每个设定的有效应力点(包括每个应力加载点及应力卸载点)处，达到内外部应力平衡，并稳定一定时间后再进行测量。其中，达到内外部应力平衡后的稳定时间可以根据岩心样品的致密程度确定。

本说明书实施例中，上述有效应力的加载和卸载可以通过定围压变内压或者变围压定内压等方式实现。

S103、确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线。

本说明书一实施例中，为了提高数据处理的便捷性，又能凸显出渗透率与有效应力之间的相关关系，所述第一关系曲线可以包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线可以包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

例如，可以将每个有效应力点下测量的渗透率除以初始渗透率，从而可以获得对应的归一化渗透率；可以将每个有效应力点除以最大有效应力，从而可以获得对应的归一化有效应力。其中，初始渗透率为初始有效应力条件下测量的渗透率。

本说明书一实施例中，在获得确定有效应力对应的渗透率的基础上，通过拟合可以获得渗透率与其对应的有效应力之间的关系曲线。例如，以归一化有效应力作为横坐标轴、归一化渗透率作为纵坐标轴，以指数形式拟合散点图，可以得到归一化渗透率与归一化有效应力之间的函数关系。

S104、根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

本说明书一实施例中，所述根据第一关系曲线和第二关系曲线确定所述应力敏感程度，包括如下步骤：

可以对所述第一关系曲线进行一阶求导，从而获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力可以用于表征所述岩心样品的应力加载敏感程度。并可以对所述第二关系曲线进行一阶求导，从而可以获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力可以用于表征所述岩心样品的应力卸载敏感程度。

本说明书一实施例中，上述临界点例如可以是拟合曲线上斜率为-1的点。

由此可见，通过本说明书实施例可以定量获得岩心样品的渗透率对不同有效应力范围的敏感程度，从而可以作为开采油/气时有效应力的警戒值，以用于指导实际生产。

本说明书一实施例中，在获得第一临界有效应力和第二临界有效应力的基础上，还可以确定所述第一临界有效应力与所述第二临界有效应力的差值，并根据所述差值确定所述岩心样品的渗透率的滞后程度。渗透率的滞后性是指同一块岩心样品在对应有效应力加载/卸载条件下，渗透率曲线不重合的现象。临界有效应力的差值量即可反映岩心渗透率对有效应力变化的滞后程度。具体而言，差值量与滞后程度正相关，即差值量越大，滞后程度越严重，反之亦然。

在油气藏开采过程中，当经过一次采油/气后，储层的渗透率会出现较大下降，为了提高采收率，在二次或三次采油/气阶段，最好将渗透率恢复到一定值。因此，通过分析储层岩心样品的渗透率的滞后程度，可以为二次采油/气补充储层的能量，恢复地层的渗透率提供重要参考，对实际油田开发具有重要的指导意义。

为便于理解，下面以定量评价某气藏储层渗透率的应力敏感程度为例进行说明。

某致密气藏，储层平均埋深4500m，地层压力系数约为1.8，构造裂缝发育。根据油气藏的实际地质和生产资料，得到储层的初始有效应力为2610psi；为保证数据完整性，并结合实验设备及实验条件，确定实验中初始有效应力为600psi，最大有效应力6600psi，有效应力间隔为1000psi。

由于该储层构造裂缝发育，为研究因裂缝发育对储层应力敏感造成的影响，故在井下取出的岩心柱上沿储层延伸方向上钻取一块裂缝岩心，该裂缝贯穿岩心两端面；同时钻取一块基质岩心，以用于将二者的应力敏感程度进行对比。

对岩心样品逐级施加600psi、1600psi、2600psi、3600psi、4600psi、5600psi及6600psi的有效应力，再以1000psi的有效应力间隔逐级降低有效应力至600psi，每个有效应力点需稳定足够长时间，以达到内外部应力平衡，然后就可以利用渗透率测试装置测量不同有效应力下岩心样品的渗透率。

对渗透率和有效应力进行归一化处理，将所测得的渗透率除以初始渗透率，有效应力除以最大有效应力，以归一化有效应力作为横坐标，归一化渗透率作为纵坐标，作出不同岩心样品的归一化渗透率与归一化有效应力的散点图，以指数形式确定归一化处理后的无因次渗透率与无因次有效应力之间的关系曲线，例如图2所示。

对以上关系曲线进行一阶求导，可以分别计算出有效应力加载和有效应力卸载曲线上临界点对应的临界归一化有效应力，例如图3所示。从图3中可以看出，裂缝岩心有效应力加载时的临界归一化有效应力为0.39(对应的加载临界有效应力为2574psi)，有效应力卸载时的临界归一化有效应力为0.23(对应的卸载临界有效应力为1518psi)；基质岩心有效应力加载时的临界归一化有效应力为0.44(对应的加载临界有效应力为2904psi)，有效应力卸载时的临界归一化有效应力为0.31(对应的卸载临界有效应力为2046psi)。

对于裂缝岩心和基质岩心，当归一化有效应力分别小于0.39(对应的加载临界有效应力为2574psi)和0.44(对应的加载临界有效应力为2904psi)时，拟合曲线的斜率小于-1，渗透率随有效应力的加载急剧减小；当归一化有效应力分别大于0.39(对应的加载临界有效应力为2574psi)和0.44(对应的加载临界有效应力为2904psi)时，拟合曲线的斜率大于-1，渗透率随有效应力的加载缓慢下降；当有效应力卸载时，当归一化有效应力降低至0.23(对应的卸载临界效应力为1518psi)以下，裂缝的岩心渗透率有明显的恢复，有效应力小于0.31(对应的卸载临界有效应力为2046psi)时，基质岩心渗透率恢复明显。在加载过程中，由于裂缝的存在，导致裂缝岩心在较小有效应力的条件下，损失了大量的渗透率，故在一次采油/气过程中，尽量控制有效应力在裂缝岩心的临界压力以下，以保证储层有较高的整体渗透率。

由于裂缝岩心和基质岩心的加载卸载归一化有效应力的差值分别为0.16(对应的有效应力为1056psi)和0.13(对应的有效应力为726psi)，相比于基质岩心，裂缝岩心在卸载过程中需要降低更多的有效应力，渗透率才能有明显的恢复。故对于该致密储层，在生产过程中生产压差不宜过大，并且要及时补充地层能量，保证地层中的天然裂缝仍有较高的导流能力。

参考图4所示，与上述确定储层应力敏感程度的方法对应，本说明书一实施例的确定储层应力敏感程度的装置可以包括：

应力范围获取模块41，可以用于获取储层的有效应力范围；

渗透率确定模块42，可以用于在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

关系曲线确定模块43，可以用于确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

敏感程度确定模块44，可以用于根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

参考图5所示，与上述确定储层应力敏感程度的方法对应，本说明书另一实施例的确定储层应力敏感程度的装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度。

虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。

为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本说明书时可以把各单元的功能在同一个或多个软件和/或硬件中实现。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。

内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。

计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁盘式存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。

还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。

本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本说明书可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

以上所述仅为本说明书的实施例而已，并不用于限制本说明书。对于本领域技术人员来说，本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书的权利要求范围之内。

Claims (7)

1.一种确定储层应力敏感程度的方法，其特征在于，包括：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度；

其中，所述根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度，包括：

对所述第一关系曲线进行求导，获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力用于表征所述岩心样品的应力加载敏感程度；所述第一关系曲线上临界点为所述第一关系曲线上归一化渗透率变化率为-1的点；

对所述第二关系曲线进行求导，获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力用于表征所述岩心样品的应力卸载敏感程度；所述第二关系曲线上临界点为所述第二关系曲线上归一化渗透率变化率为-1的点；

确定所述第一临界有效应力与所述第二临界有效应力的差值；

根据所述差值确定所述岩心样品的渗透率的滞后程度；

所述第一关系曲线包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

2.如权利要求1所述的确定储层应力敏感程度的方法，其特征在于，所述获取储层的有效应力范围，包括：

获取所述储层的上覆岩层压力、初始地层孔隙流体压力和实际生产压力；

根据所述上覆岩层压力和所述初始地层孔隙流体压力确定所述储层的初始有效应力；

根据所述上覆岩层压力和所述实际生产压力确定所述储层的最大有效应力；

根据所述初始有效应力和所述最大有效应力确定所述储层的有效应力范围。

3.如权利要求1所述的确定储层应力敏感程度的方法，其特征在于，所述应力加载序列包括：由多个等间隔的应力加载点形成的应力递增序列；所述应力卸载序列包括：由多个等间隔的应力卸载点形成的应力递减序列；所述应力递增序列及所述应力递减序列中的最小值为所述有效应力范围的最小值，所述应力递增序列及所述应力递减序列中的最大值为所述有效应力范围的最大值。

4.如权利要求3所述的确定储层应力敏感程度的方法，其特征在于，所述应力递增序列中的应力加载点与所述应力递减序列中的应力卸载点对应。

5.一种确定储层应力敏感程度的装置，其特征在于，包括：

应力范围获取模块，用于获取储层的有效应力范围；

渗透率确定模块，用于在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

关系曲线确定模块，用于确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

敏感程度确定模块，用于根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度；

其中，所述根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度，包括：

对所述第一关系曲线进行求导，获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力用于表征所述岩心样品的应力加载敏感程度；所述第一关系曲线上临界点为所述第一关系曲线上归一化渗透率变化率为-1的点；

对所述第二关系曲线进行求导，获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力用于表征所述岩心样品的应力卸载敏感程度；所述第二关系曲线上临界点为所述第二关系曲线上归一化渗透率变化率为-1的点；

确定所述第一临界有效应力与所述第二临界有效应力的差值；

根据所述差值确定所述岩心样品的渗透率的滞后程度；

所述第一关系曲线包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

6.一种确定储层应力敏感程度的装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，其特征在于，所述计算机程序被所述处理器运行时执行如下步骤：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度；

其中，所述根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度，包括：

对所述第一关系曲线进行求导，获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力用于表征所述岩心样品的应力加载敏感程度；所述第一关系曲线上临界点为所述第一关系曲线上归一化渗透率变化率为-1的点；

对所述第二关系曲线进行求导，获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力用于表征所述岩心样品的应力卸载敏感程度；所述第二关系曲线上临界点为所述第二关系曲线上归一化渗透率变化率为-1的点；

确定所述第一临界有效应力与所述第二临界有效应力的差值；

根据所述差值确定所述岩心样品的渗透率的滞后程度；

所述第一关系曲线包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

7.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤：

获取储层的有效应力范围；

在所述有效应力范围内，确定所述储层的岩心样品在以预设的应力加载序列作用下的第一渗透率，并确定所述岩心样品在以预设的应力卸载序列作用下的第二渗透率；

确定所述第一渗透率与其对应的有效应力的第一关系曲线，并确定所述第二渗透率与其对应的有效应力的第二关系曲线；

根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度；

其中，所述根据所述第一关系曲线和所述第二关系曲线确定所述应力敏感程度，包括：

对所述第一关系曲线进行求导，获得所述第一关系曲线上临界点对应的第一临界有效应力，所述第一临界有效应力用于表征所述岩心样品的应力加载敏感程度；所述第一关系曲线上临界点为所述第一关系曲线上归一化渗透率变化率为-1的点；

对所述第二关系曲线进行求导，获得所述第二关系曲线上临界点对应的第二临界有效应力，所述第二临界有效应力用于表征所述岩心样品的应力卸载敏感程度；所述第二关系曲线上临界点为所述第二关系曲线上归一化渗透率变化率为-1的点；

确定所述第一临界有效应力与所述第二临界有效应力的差值；

根据所述差值确定所述岩心样品的渗透率的滞后程度；

所述第一关系曲线包括：归一化第一渗透率及其对应的归一化有效应力之间的关系曲线；所述第二关系曲线包括：归一化第二渗透率及其对应的归一化有效应力之间的关系曲线。

Images