CN110146940B - 一种气藏束缚水饱和度、气藏可动水范围确定方法和装置 - Google Patents

Abstract

本发明公开了一种气藏束缚水饱和度、气藏可动水范围确定方法和装置。气藏束缚水饱和度确定方法包括：按照预设规则将研究区划分为多个网格，目标井的位置位于研究区和多个网格的中央；根据网格与目标井的距离、目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；根据网格的有效应力和有效应力影响系数、研究区的孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；根据网格受有效应力影响的毛细管半径分布的分形维数和毛细管迂曲度的分形维数、束缚水膜厚度，确定网格的束缚水饱和度。能够确定不同生产时期、不同位置处的束缚水饱和度。

Description

一种气藏束缚水饱和度、气藏可动水范围确定方法和装置

技术领域

本发明涉及气藏开发技术领域，特别涉及一种气藏束缚水饱和度、气藏可动水范围确定方法和装置。

背景技术

虽然气藏往往具有较高的束缚水饱和度，但在生产过程中，生产压力的改变会使得气藏的部分束缚水变为可动水，而气藏大范围出现可动水会造成不同程度的水锁效应，进而导致气井减产或停产。因此，气藏生产过程中可动水范围的预测，对于降低气藏产水风险和延长稳产时间具有重要意义。

现有技术中，往往是在生产现场中利用密闭取心或测井方法确定气藏原始含水饱和度；通过室内实验，例如核磁共振法、离心称重法、加湿气驱水动态法、压汞法和半渗隔板法等实验方法进行气藏束缚水饱和度的计算；再根据原始含水饱和度与束缚水饱和度来确定可动水范围。其中，计算的气藏束缚水饱和度是在特定条件下的特定位置的束缚水饱和度，但在实际生产过程中，不同生产阶段、不用位置处的压力梯度是不同的，对束缚水饱和度变化的影响也不同，而现有技术无法确定任意生产阶段、任意位置处的束缚水饱和度。

发明内容

鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的一种气藏束缚水饱和度、气藏可动水范围确定方法和装置。

第一方面，本发明实施例提供一种气藏束缚水饱和度确定方法，包括：

按照预设规则将研究区划分为多个网格，目标井的位置位于所述研究区和所述多个网格的中央；

根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；

根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，确定网格的有效应力影响系数，所述有效应力影响系数用于根据参数的原始值确定受有效应力影响的所述参数的值；

根据网格的所述有效应力和有效应力影响系数、所述研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；

根据网格的所述压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；

根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

第二方面，本发明实施例提供一种气藏可动水范围确定方法，包括：

利用上述气藏束缚水饱和度确定方法确定每个网格的束缚水饱和度；

确定束缚水饱和度小于原始含水饱和度的网格所在的位置为可动水范围。

第三方面，本发明实施例提供一种气藏束缚水饱和度确定装置，包括：

划分模块，用于按照预设规则将研究区划分为多个网格，目标井的位置位于所述研究区和所述多个网格的中央；

第一确定模块，用于根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；

第二确定模块，用于根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，确定网格的有效应力影响系数，所述有效应力影响系数用于根据参数的原始值确定受有效应力影响的所述参数的值；

第三确定模块，用于根据网格的所述有效应力和有效应力影响系数、所述研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；

第四确定模块，用于根据网格的所述压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；

第五确定模块，用于根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

第四方面，本发明实施例提供一种计算机可读存储介质，其上存储有计算机指令，当该指令被处理器执行时实现上述气藏束缚水饱和度确定方法。

本发明实施例提供的上述技术方案的有益效果至少包括：

按照预设规则将研究区划分为多个网格，目标井的位置位于研究区和多个网格的中央；根据网格与目标井的距离、目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；根据研究区的岩石泊松比和杨氏模量，以及网格的有效应力，确定网格的有效应力影响系数，有效应力影响系数用于根据参数的原始值确定受有效应力影响的参数的值；根据网格的有效应力和有效应力影响系数、研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；根据网格的压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。能够有效的确定不同生产时期、不同位置处的束缚水饱和度，为开发井附近可动水范围的确定提供了参考，进而可以有效的指导开发。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：

图1为本发明实施例一中所述气藏束缚水饱和度确定方法的流程图；

图2为本发明实施例一中离散网格示例图；

图3为本发明实施例一中数据准备工作的流程图；

图4为本发明实施例二中所述气藏可动水范围确定方法的流程图；

图5为本发明实施例二中确定的束缚水饱和度值分布图；

图6为本发明实施例二中气藏可动水范围确定方法的示意图；

图7为本发明实施例中所述气藏束缚水饱和度确定装置的结构图。

具体实施方式

下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。

为了解决现有技术无法确定气藏任意生产阶段、任意位置处的束缚水饱和度的问题，本发明实施例提供了一种气藏束缚水饱和度确定方法，该方案能够确定气藏不同生产时期、不同位置处的可动水饱和度，为气藏开发过程中可动水范围的确定提供了参考，进而可以有效的指导气藏开发。

实施例一

本发明实施例一提供一种气藏束缚水饱和度确定方法，其流程如图1所示，包括如下步骤：

步骤S11：按照预设规则将研究区划分为多个网格，目标井的位置位于研究区和多个网格的中央。

具体的，可以是利用编程软件或者是建模软件，例如利用matlab编程软件，将研究区即目标井区平面上离散为n个网格，网格步长为l，目标井well位于研究区中央位置，故目标井也位于网格中央位置。

其中，研究区是将目标层位模型化成一个二维区域。

参照图2所示，可以是，以目标井well为中心，向四周离散网格，每一个网格为一个小方形，每个网格的长即为网格步长。每个网格内的孔隙结构特征和毛细管分布特征如图2中所示，下述步骤S12-步骤S16根据每个网格的孔隙结构特征和毛细管分布特征，利用预先建立的数学模型(后面详细介绍的每个数学方程)确定每个网格的束缚水饱和度S_wir。

在进行下述步骤S12-步骤S16之前，要先进行数据的准备工作，具体步骤参照图3所示，可以是：

步骤S31：采集目标层位的岩心样本、地层水样本。

步骤S32：对采集的岩心样本和地层水样本进行室内实验研究，通过常规实验获得研究区的预设参数的值。

具体的，预设参数包括：通过常规实验获得研究区的孔隙结构特征指数β；通过压汞实验获得研究区的毛细管半径最大值r_max、毛细管半径最小值r_min和毛细管半径平均值r_ave；通过岩心力学实验获得研究区的岩石泊松比ν和杨氏模量E；通过润湿性实验获得研究区的岩石润湿角θ；通过界面张力实验获得研究区的气水界面张力σ。

步骤S33：收集目标井基本数据。

具体的，包括目标井的当前井底压力p_w、目标井的表皮系数s、目标井井眼半径R_w、目标井的泄气半径R_e、目标井的原始地层压力p_e、原始地层含水饱和度Swi、目标井的地层上覆岩石压力p₀和地层温度T。

步骤S12：根据网格与目标井的距离、目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力。

具体的，可以是，根据公式(1)确定第i个网格的压力p_i：

根据公式(2)确定第i个网格的压差

根据网格步长和第i个网格的压差

利用公式(3)确定第i个网格的压力梯度Δp_i：

根据目标井的上覆岩石压力和第i个网格的压力p_i，利用公式(4)确定第i个网格的有效应力

上式(1)-(4)中，p_e为目标井的原始地层压力；p_w为目标井的当前井底压力；x_i为第i个网格与目标井的距离；R_e为目标井的泄气半径；l为网格步长；R_rw为目标井的折算井眼半径，R_rw＝R_we^-s，R_w为目标井井眼半径，s为目标井的表皮系数；p₀为目标井的上覆岩石压力。

由于钻井、完井及井下作业对地层的污染或改善，近井地层的渗透率将发生变化，因此产生附加阻力。可以设想井壁贴一层表皮，流体流过它时所产生的附加阻力正好等于因近井地层渗透率变化所产生的附加阻力。引入表皮后可以认为近井地层的渗透率未发生变化，从而避免了因近井地层渗透率发生变化所造成的数学处理困难。表皮所造成的阻力大小由表皮系数S表示。

步骤S13：根据研究区的岩石泊松比和杨氏模量，以及网格的有效应力，确定网格的有效应力影响系数。

上述确定的有效应力影响系数用于根据参数的原始值确定受有效应力影响的参数的值。例如，

为第i个网格受有效应力影响的毛细管半径最大值，r_max为所述研究区的毛细管半径最大值；

为第i个网格受有效应力影响的毛细管半径最小值，r_min为所述研究区的毛细管半径最小值。

具体有效应力影响系数确定方法可以是：

根据公式(5)确定第i个网格的有效应力影响系数εⁱ：

公式(5)中，ν为研究区的岩石泊松比，E为研究区的杨氏模量；

为第i个网格的有效应力；β为研究区的孔隙结构特征指数。

步骤S14：根据网格的有效应力和有效应力影响系数、研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数。

确定网格受有效应力影响的毛细管半径分布的分形维数，可以是：

根据公式(6)-公式(8)联立求解，确定第i个网格受有效应力影响的毛细管半径分布的分形维数

上式(6)-(8)中，

为第i个网格原始毛细管半径分布的分形维数；d为欧拉系数；

为第i个网格受有效应力影响的毛细管半径修正值；

为第i个网格受有效应力影响的毛细管半径最大值，

r_max为研究区的毛细管半径最大值；

为第i个网格受有效应力影响的毛细管半径最小值，

r_min为研究区的毛细管半径最小值。

确定网格受有效应力影响的毛细管迂曲度的分形维数，包括：

根据第i个网格受有效应力影响的毛细管半径修正值

利用公式(9)确定第i个网格受有效应力影响的毛细管平均迂曲度

根据公式(10)确定第i个网格受有效应力影响的毛细管表观长度

据公式(11)确定第i个网格受有效应力影响的毛细管半径平均值

上式(11)中，

r_ave为研究区的毛细管半径平均值。

根据第i个网格受有效应力影响的毛细管平均迂曲度

受有效应力影响的毛细管表观长度

和受有效应力影响的毛细管半径平均值

利用公式(12)确定第i个网格受有效应力影响的毛细管迂曲度的分形维数

步骤S15：根据网格的压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度。

具体的，可以是，根据研究区的地层温度T，利用公式(13)确定研究区的地层水粘度μ_w：

μ_w＝1.4e^-0.0176T (13)

对应的，根据公式(14)确定第i个网格中每一毛细管半径下的束缚水膜厚度δ_i：

公式(14)中，r_i为第i个网格的毛细管半径；

为第i个网格的压差。

步骤S16：根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

具体的，可以是，先根据网格的压差、研究区的气水界面张力和岩石润湿角，利用公式(13)确定第i个网格在当前压差条件下的临界毛细管半径

公式(15)中，Δp_i为第i个网格的压差；σ为研究区的气水界面张力；θ为研究区的岩石润湿角；临界毛细管半径是在某个压力下气体能进入毛细管的半径和不能进入毛细管的半径直接的界限，即在某个压力下低于临界毛细管半径气体不能进入毛细管，等于或高于临界毛细管半径气体能进入毛细管。

对应的，根据公式(16)确定第i个网格的束缚水饱和度

公式(16)中，r_i为第i个网格的毛细管半径，

为第i个网格受有效应力影响的毛细管半径，

为第i个网格受有效应力影响的毛细管半径最大值，

为第i个网格受有效应力影响的毛细管半径最小值；

为第i个网格受有效应力影响的毛细管半径分布的分形维数；

为第i个网格受有效应力影响的毛细管迂曲度的分形维数；δ_i为第i个网格中每一毛细管半径下的束缚水膜厚度；

第i个网格在当前压差条件下的临界毛细管半径，

为第i个网格在受有效应力影响的临界毛细管半径，

本发明实施例的一种气藏束缚水饱和度确定方法，按照预设规则将研究区划分为多个网格，目标井的位置位于研究区和多个网格的中央；根据网格与目标井的距离、目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；根据研究区的岩石泊松比和杨氏模量，以及网格的有效应力，确定网格的有效应力影响系数，有效应力影响系数用于根据参数的原始值确定受有效应力影响的参数的值；根据网格的有效应力和有效应力影响系数、研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；根据网格的压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。能够有效的确定不同生产时期、不同位置处的束缚水饱和度，为开发井附近可动水范围的确定提供了参考，进而可以有效的指导开发。

实施例二

本发明实施例二提供一种气藏可动水范围确定方法，其流程如图4所示，包括如下步骤：

步骤S41：按照预设规则将研究区划分为多个网格，目标井的位置位于研究区和多个网格的中央。

步骤S42：根据网格与目标井的距离、目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力。

步骤S43：根据研究区的岩石泊松比和杨氏模量，以及网格的有效应力，确定网格的有效应力影响系数。

步骤S44：根据网格的有效应力和有效应力影响系数、研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数。

步骤S45：根据网格的压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度。

步骤S46：根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

上述步骤S41-步骤S46同实施例一中的步骤S11-步骤S16，此处不做赘述。

步骤S47：确定束缚水饱和度小于原始含水饱和度的网格所在的位置为可动水范围。

图5为最终确定的每个网格的束缚水饱和度确定值，即研究区内确定的束缚水饱和度值的分布，图中的颜色代表确定的束缚水饱和度值的大小。可以看出随着到目标井距离的减小，确定的束缚水饱和度值是逐渐减小的。

通常情况下，原始含水饱和度等于束缚水饱和度加可动水饱和度，故，上述步骤S41-步骤S46确定的每个网格的束缚水饱和度，理论上应该是不高于原始含水饱和度的，但由于计算存在误差，可能计算出的束缚水饱和度存在高于原始含水饱和度的情况。例如图6所示，横坐标为网格与目标井的距离，纵坐标为网格的束缚水饱和度确定值，可以看出随着网格与目标井的距离的增大，束缚水饱和度值也是越来越大的，当束缚水饱和度值高于原始含水饱和度值时，理论上是不存在这种情况的，这是由于计算的误差所致，可以认为计算的束缚水饱和度值高于原始含水饱和度值时，真实的束缚水饱和度值约等于原始含水饱和度值，即此时不存在可动水；当计算的束缚水饱和度值低于原始含水饱和度值时，确定为存在可动水，即图6中曲线中的点为可动水范围的临界点，当网格到目标井的距离小于图6中曲线临界点到目标井的距离时，确定为可动水范围。

基于同一发明构思，本发明实施例还提供一种气藏束缚水饱和度确定装置，该装置的结构如图7所示，包括：

划分模块71，用于按照预设规则将研究区划分为多个网格，目标井的位置位于所述研究区和所述多个网格的中央；

第一确定模块72，用于根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；

第二确定模块73，用于根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，确定网格的有效应力影响系数，所述有效应力影响系数用于根据参数的原始值确定受有效应力影响的所述参数的值；

第三确定模块74，用于根据网格的所述有效应力和有效应力影响系数、所述研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；

第四确定模块75，用于根据网格的所述压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；

第五确定模块66，用于根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。

基于同一发明构思，本发明实施例还提供一种计算机可读存储介质，其上存储有计算机指令，当该指令被处理器执行时实现上述气藏束缚水饱和度确定方法。

除非另外具体陈述，术语比如处理、计算、运算、确定、显示等等可以指一个或更多个处理或者计算系统、或类似设备的动作和/或过程，所述动作和/或过程将表示为处理系统的寄存器或存储器内的物理(如电子)量的数据操作和转换成为类似地表示为处理系统的存储器、寄存器或者其他此类信息存储、发射或者显示设备内的物理量的其他数据。信息和信号可以使用多种不同的技术和方法中的任何一种来表示。例如，在贯穿上面的描述中提及的数据、指令、命令、信息、信号、比特、符号和码片可以用电压、电流、电磁波、磁场或粒子、光场或粒子或者其任意组合来表示。

应该明白，公开的过程中的步骤的特定顺序或层次是示例性方法的实例。基于设计偏好，应该理解，过程中的步骤的特定顺序或层次可以在不脱离本公开的保护范围的情况下得到重新安排。所附的方法权利要求以示例性的顺序给出了各种步骤的要素，并且不是要限于所述的特定顺序或层次。

在上述的详细描述中，各种特征一起组合在单个的实施方案中，以简化本公开。不应该将这种公开方法解释为反映了这样的意图，即，所要求保护的主题的实施方案需要清楚地在每个权利要求中所陈述的特征更多的特征。相反，如所附的权利要求书所反映的那样，本发明处于比所公开的单个实施方案的全部特征少的状态。因此，所附的权利要求书特此清楚地被并入详细描述中，其中每项权利要求独自作为本发明单独的优选实施方案。

本领域技术人员还应当理解，结合本文的实施例描述的各种说明性的逻辑框、模块、电路和算法步骤均可以实现成电子硬件、计算机软件或其组合。为了清楚地说明硬件和软件之间的可交换性，上面对各种说明性的部件、框、模块、电路和步骤均围绕其功能进行了一般地描述。至于这种功能是实现成硬件还是实现成软件，取决于特定的应用和对整个系统所施加的设计约束条件。熟练的技术人员可以针对每个特定应用，以变通的方式实现所描述的功能，但是，这种实现决策不应解释为背离本公开的保护范围。

结合本文的实施例所描述的方法或者算法的步骤可直接体现为硬件、由处理器执行的软件模块或其组合。软件模块可以位于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、移动磁盘、CD-ROM或者本领域熟知的任何其它形式的存储介质中。一种示例性的存储介质连接至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于ASIC中。该ASIC可以位于用户终端中。当然，处理器和存储介质也可以作为分立组件存在于用户终端中。

对于软件实现，本申请中描述的技术可用执行本申请所述功能的模块(例如，过程、函数等)来实现。这些软件代码可以存储在存储器单元并由处理器执行。存储器单元可以实现在处理器内，也可以实现在处理器外，在后一种情况下，它经由各种手段以通信方式耦合到处理器，这些都是本领域中所公知的。

上文的描述包括一个或多个实施例的举例。当然，为了描述上述实施例而描述部件或方法的所有可能的结合是不可能的，但是本领域普通技术人员应该认识到，各个实施例可以做进一步的组合和排列。因此，本文中描述的实施例旨在涵盖落入所附权利要求书的保护范围内的所有这样的改变、修改和变型。此外，就说明书或权利要求书中使用的术语“包含”，该词的涵盖方式类似于术语“包括”，就如同“包括，”在权利要求中用作衔接词所解释的那样。此外，使用在权利要求书的说明书中的任何一个术语“或者”是要表示“非排它性的或者”。

Claims (10)

1.一种气藏束缚水饱和度确定方法，其特征在于，包括：

按照预设规则将研究区划分为多个网格，目标井的位置位于所述研究区和所述多个网格的中央；

根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；

根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，通过下式确定网格的有效应力影响系数：

上式中，εⁱ为第i个网格的有效应力影响系数；ν为所述研究区的岩石泊松比，E为所述研究区的杨氏模量；

为第i个网格的有效应力；β为所述研究区的孔隙结构特征指数；所述有效应力影响系数用于根据参数的原始值确定受有效应力影响的所述参数的值；

根据网格的所述有效应力和有效应力影响系数、所述研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；

根据网格的所述压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；

根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

2.如权利要求1所述的方法，其特征在于，根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力，包括：

根据公式(1)确定第i个网格的压力p_i：

根据公式(2)确定第i个网格的压差

根据网格步长和第i个网格的压差

利用公式(3)确定第i个网格的压力梯度Δp_i：

根据所述目标井的上覆岩石压力和第i个网格的压力p_i，利用公式(4)确定第i个网格的有效应力

上式(1)-(4)中，p_e为所述目标井的原始地层压力；p_w为所述目标井的当前井底压力；x_i为第i个网格与所述目标井的距离；R_e为所述目标井的泄气半径；l为网格步长；R_rw为所述目标井的折算井眼半径，R_rw＝R_we^-s，R_w为所述目标井井眼半径，s为所述目标井的表皮系数；p₀为所述目标井的上覆岩石压力。

3.如权利要求1所述的方法，其特征在于，根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，确定网格的有效应力影响系数，包括：

根据公式(5)确定第i个网格的有效应力影响系数εⁱ：

公式(5)中，ν为所述研究区的岩石泊松比，E为所述研究区的杨氏模量；

为第i个网格的有效应力；β为所述研究区的孔隙结构特征指数。

4.如权利要求3所述的方法，其特征在于，确定网格受有效应力影响的毛细管半径分布的分形维数，包括：

根据公式(6)-公式(8)联立求解，确定第i个网格受有效应力影响的毛细管半径分布的分形维数

上式(6)-(8)中，

为第i个网格原始毛细管半径分布的分形维数；d为欧拉系数；

为第i个网格受有效应力影响的毛细管半径修正值；

为第i个网格受有效应力影响的毛细管半径最大值，

r_max为所述研究区的毛细管半径最大值；

为第i个网格受有效应力影响的毛细管半径最小值，

r_min为所述研究区的毛细管半径最小值。

5.如权利要求4所述的方法，其特征在于，确定网格受有效应力影响的毛细管迂曲度的分形维数，包括：

根据第i个网格受有效应力影响的毛细管半径修正值

利用公式(9)确定第i个网格受有效应力影响的毛细管平均迂曲度

根据公式(10)确定第i个网格受有效应力影响的毛细管表观长度

据公式(11)确定第i个网格受有效应力影响的毛细管半径平均值

根据第i个网格受有效应力影响的毛细管平均迂曲度

受有效应力影响的毛细管表观长度

和受有效应力影响的毛细管半径平均值

利用公式(12)确定第i个网格受有效应力影响的毛细管迂曲度的分形维数

6.如权利要求1所述的方法，其特征在于，确定网格的束缚水膜厚度前，还包括：

根据所述研究区的地层温度T，利用公式(13)确定所述研究区的地层水粘度μ_w：

μ_w＝1.4e^-0.0176T (13)

对应的，根据公式(14)确定第i个网格中每一毛细管半径下的束缚水膜厚度δ_i：

公式(14)中，r_i为第i个网格的毛细管半径；

为第i个网格的压差。

7.如权利要求3所述的方法，其特征在于，确定网格的束缚水饱和度前，还包括：

根据网格的所述压差、所述研究区的气水界面张力和岩石润湿角，利用公式(15 )确定第i个网格在当前压差条件下的临界毛细管半径

公式(15)中，Δp_i为第i个网格的压差；σ为所述研究区的气水界面张力；θ为所述研究区的岩石润湿角；

对应的，根据公式(16)确定第i个网格的束缚水饱和度

公式(16)中，r_i为第i个网格的毛细管半径，

为第i个网格受有效应力影响的毛细管半径，

为第i个网格受有效应力影响的毛细管半径最大值，

为第i个网格受有效应力影响的毛细管半径最小值；

为第i个网格受有效应力影响的毛细管半径分布的分形维数；

为第i个网格受有效应力影响的毛细管迂曲度的分形维数；δ_i为第i个网格中每一毛细管半径下的束缚水膜厚度；

第i个网格在当前压差条件下的临界毛细管半径，

为第i个网格在受有效应力影响的临界毛细管半径，

8.一种气藏可动水范围确定方法，其特征在于，包括：

利用如权利要求1～7任一所述的气藏束缚水饱和度确定方法确定每个网格的束缚水饱和度；

确定束缚水饱和度小于原始含水饱和度的网格所在的位置为可动水范围。

9.一种气藏束缚水饱和度确定装置，其特征在于，包括：

划分模块，用于按照预设规则将研究区划分为多个网格，目标井的位置位于所述研究区和所述多个网格的中央；

第一确定模块，用于根据网格与所述目标井的距离、所述目标井的当前井底压力、原始地层压力、井眼半径和泄气半径，确定网格的压力、压差、压力梯度和有效应力；

第二确定模块，用于根据所述研究区的岩石泊松比和杨氏模量，以及网格的所述有效应力，通过下式确定网格的有效应力影响系数：

上式中，εⁱ为第i个网格的有效应力影响系数；ν为所述研究区的岩石泊松比，E为所述研究区的杨氏模量；

为第i个网格的有效应力；β为所述研究区的孔隙结构特征指数；所述有效应力影响系数用于根据参数的原始值确定受有效应力影响的所述参数的值；

第三确定模块，用于根据网格的所述有效应力和有效应力影响系数、所述研究区的最大孔隙半径和最小孔隙半径，确定网格受有效应力影响的毛细管半径分布的分形维数和受有效应力影响的毛细管迂曲度的分形维数；

第四确定模块，用于根据网格的所述压力梯度、网格的毛细管半径，确定网格的束缚水膜厚度；

第五确定模块，用于根据网格受有效应力影响的毛细管半径分布的分形维数、受有效应力影响的毛细管迂曲度的分形维数和束缚水膜厚度，确定网格的束缚水饱和度。

10.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，当该指令被处理器执行时实现权利要求1～7任一所述的气藏束缚水饱和度确定方法。

Images