CN103257089B - 压力脉冲测量装置测量基质和裂缝渗透率的方法 - Google Patents

Abstract

本发明涉及一种压力脉冲测量装置，包括通过密封管线依次连接的气瓶、压力调节阀、气动阀Ⅰ、气动阀Ⅱ、压力传感器Ⅱ、气动阀Ⅲ、压力传感器Ⅲ，还包括连接在压力传感器Ⅱ和压力传感器Ⅲ之间的岩心夹持器、环压阀、环压泵，所述环压阀控制环压泵和岩心夹持器的连通，所述压力传感器Ⅱ和压力传感器Ⅲ测量岩心夹持器两端压力。还涉及一种利用压力脉冲测量装置同时测量基质和裂缝渗透率的方法，通过绘制压力-时间关系的实测曲线和理论曲线得出基质渗透率k_m和裂缝渗透率k_f。本发明只需对一块具有裂缝的岩心进行一次测试就能分别测出裂缝、基质的渗透率，提高了测量效率。

Description

压力脉冲测量装置测量基质和裂缝渗透率的方法

技术领域

本发明涉及油气田开发实验测定常规物性参数的特殊实验方法领域，具体涉及一种压力脉冲测量装置及利用其测量基质和裂缝渗透率的方法。

背景技术

双重介质模型是油气藏开发工程长期关注的一个重要模型，但是针对同一岩样的渗透率中基质渗透率、裂缝渗透率各占多少，却没有现成的方法可以实现。原有的常规渗透率测试装置以及压力脉冲衰减方法，并不能排除裂缝对渗透率的干扰，测试结果往往体现出裂缝的特征，而并非基质的渗透率，因此目前急需一种可以同时测试同一块双重介质样品的基质渗透率、裂缝渗透率的装置和方法。

发明内容

为了解决以上技术问题，本发明提供一种压力脉冲测量装置及利用其测量基质和裂缝渗透率的方法，对一块具有裂缝的岩心进行一次测试就能分别测出裂缝、基质的渗透率，提高了测量效率。

本发明是通过以下技术方案实现的：

一种压力脉冲测量装置，包括通过密封管线依次连接的气瓶、压力调节阀、气动阀Ⅰ、气动阀Ⅱ、压力传感器Ⅱ、气动阀Ⅲ、压力传感器Ⅲ，还包括连接在压力传感器Ⅱ和压力传感器Ⅲ之间的岩心夹持器、环压阀、环压泵，所述环压阀控制环压泵和岩心夹持器的连通，所述压力传感器Ⅱ和压力传感器Ⅲ测量岩心夹持器两端压力。

在上述技术方案中，还包括连接在气动阀Ⅰ和气动阀Ⅱ之间的压力传感器Ⅰ、手动阀和标准室，所述手动阀控制压力传感器Ⅰ和标准室的连通。

在上述技术方案中，还包括连接在岩心夹持器和环压阀之间的温度传感器。

本发明还提供了一种利用压力脉冲测量装置测量基质和裂缝渗透率的方法方法，包括如下步骤：

在恒温条件下通过所述环压泵提供压力使所述岩心夹持器对岩心柱塞样加环压p_h；

开启所述气动阀Ⅰ、气动阀Ⅱ、气动阀Ⅲ，调节所述压力调节阀为p₁稳定；

关闭所述气动阀Ⅱ、气动阀Ⅲ；

调节所述压力调节阀为p₁’稳定，关闭所述气动阀Ⅰ；

瞬间打开所述气动阀Ⅱ，施加脉冲压力，通过所述压力传感器Ⅱ、压力传感器Ⅲ采集所述岩心夹持器进口端和出口端的压力，直至两端压力平衡；

两端压力平衡后，所述压力传感器Ⅱ、压力传感器Ⅲ继续采集压力值直到无压降发生；

将所述压力传感器Ⅱ、压力传感器Ⅲ记录的压力、时间无因次化并绘制压力-时间关系的实测曲线，根据理论计算的所述岩心夹持器进口端和出口端的无因次压力和无因次时间绘制理论曲线，将实测曲线与理论曲线做对比，拟合好的渗透率就是该岩心样品的对应的基质渗透率k_m和裂缝渗透率k_f。

在上述技术方案中，所述实测曲线通过如下方法绘制：

用拟压力形式表示其中μ是气体粘度，单位是mPa·s，Z为气体偏差因子，单位是1，p为压力传感器Ⅱ或Ⅲ记录的实测气体压力，单位为MPa，p₀为压力传感器Ⅱ或Ⅲ记录的起始压力，单位为MPa；

通过公式将所述拟压力无因次化为无因次压力；

通过公式将传感器记录的时间t无因次化为无因次时间，式中，所述代表t时刻岩心夹持器出口端或进口端的无因次压力值，所述代表0时刻岩心夹持器出口端或进口端的拟压力值，所述是指在没加脉冲压力前岩心夹持器夹持的岩心柱塞样的拟压力值，所述t_D是无因次时间，所述k_m是岩心柱塞样基质渗透率，单位是mD，所述φ_m代表岩心柱塞样基质的孔隙度，单位是1，所述μ_m是岩心柱塞样基质的气体粘度，单位是mPa·s，所述C_g为气体体积压缩系数，数值上等于压力传感器Ⅱ记录的压力平均值和压力传感器Ⅲ记录的压力平均值之和的倒数，单位是MPa^-1，所述l为岩心柱塞样的高，角标上D代表无因次量；所述岩心夹持器进口端即岩心夹持器靠近脉冲压力施加的一端，所述岩心夹持器出口端即岩心夹持器远离脉冲压力施加的一端。

根据所述无因次拟压力和无因次时间t_D绘制实测曲线。

在上述技术方案中，所述理论曲线通过如下方法绘制：

①通过公式分别计算无因次基质压力无因次裂缝压力

公式一：

公式二：

上式中 $L_1(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k+n}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]\mathrm{\λ}{\mathrm{\π}}^2{t}_D}{\mathrm{\ω}{\mathrm{\ξ}}_m}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$

$L_2(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]\mathrm{\λ}{\mathrm{\π}}^2{t}_D}{\mathrm{\ω}{\mathrm{\ξ}}_m}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$ 所述x_D代表岩心柱塞样轴向上的无因次坐标，所述l为岩心柱塞样的高，所述z_D代表岩心柱塞样径向上的无因次坐标，所述所述t_D代表无因次时间，k＝1,2,3…∞，n＝1,2,3…∞，所述角标f代表裂缝渗透率相关数值，所述角标m代表基质渗透率相关数值，所述角标u代表岩心夹持器进口端，角标d代表岩心夹持器出口端，*为卷积运算符，代表岩心夹持器进口端无因次拟压力，代表岩心夹持器出口端无因次拟压力；

②再将所述代入下面的积分-微分方程组Ⅰ

以上式中， ${\mathrm{\ξ}}_f=\frac{\stackrel{\‾}{{\mathrm{\μ}}_f{C}_{\mathrm{tf}}}}{\stackrel{\‾}{{\mathrm{\μ}}_0{C}_{t0}}},{\mathrm{\ξ}}_m=\frac{{\mathrm{\μ}}_m{p}_0}{{\mathrm{\μ}}_0{p}_m},\mathrm{\λ}=\frac{k_m{h}_m}{k_f{h}_f},\mathrm{\ω}=\frac{{\mathrm{\φ}}_m{\stackrel{\‾}{p}}_f{h}_m}{{\mathrm{\φ}}_f{\stackrel{\‾}{p}}_m{h}_f},{\mathrm{\ξ}}_d=\frac{\stackrel{\‾}{{\mathrm{\μ}}_d}{p}_0}{{\mathrm{\μ}}_0\stackrel{\‾}{p_d}},$ 其中所述φ_m代表岩心柱塞样基质的孔隙度，单位为1，所述φ_f代表岩心柱塞样的裂缝孔隙度，单位为1，所述C_tf代表岩心柱塞样的裂缝中气体的体积压缩系数，单位为MPa^-1，所述C_t0代表未加脉冲压力前岩心柱塞样基质中气体的体积压缩系数，单位为MPa^-1，所述p₀是未加脉冲压力前岩心夹持器出口端的气体压力，单位为MPa，所述p_f是指岩心柱塞样的裂缝中气体的压力，单位为MPa，所述p_m是岩心柱塞样基质中气体的压力，单位为MPa，所述p_d是岩心夹持器出口端气体压力，单位为MPa，所述h_f是岩心柱塞样的裂缝的半长，所述h_m是岩心柱塞样的半径减去裂缝的半长，所述μ₀是未加脉冲压力前岩心柱塞样中的气体的粘度，单位是mPa·s，所述μ_d是岩心夹持器出口端气体粘度，单位是mPa·s，所述μ_m是岩心柱塞样基质内气体粘度，单位是mPa·s，所述μ_f是岩心柱塞样的裂缝内的气体粘度，单位是mPa·s，所述x是圆柱体样品岩心高度方向上的坐标值，所述z代表样品岩心直径方向上的坐标值，l为岩心柱塞样的高，所述V_f为岩心柱塞样的裂缝的体积，所述V_u岩心夹持器进口端体积，所述V_d岩心夹持器出口端体积；所述岩心夹持器进口端体积即气动阀Ⅱ到气动阀Ⅲ之间管线的体积加上压力传感器Ⅱ到岩心夹持器进口端之间管线的体积，所述岩心夹持器出口端体积即所述气动阀Ⅲ到所述压力传感器Ⅲ之间管线的体积加上所述压力传感器Ⅲ到所述岩心夹持器出口端之间管线的体积；

将带有所述为变量的积分-微分方程组Ⅰ离散化，得到以下方程组：

其中，

k＝1，2，3…M，j＝1，2，3…L，

把实验总共所需时间t代入公式所得到的无因次时间t_D值就是T，式中所述k_m是岩心柱塞样基质渗透率，单位是mD，所述k_f为岩心柱塞样的裂缝渗透率，单位是mD，所述h_f是岩心柱塞样的裂缝的半长，所述h_m是岩心柱塞样的半径减去裂缝的半长，所述φ_f代表岩心柱塞样的裂缝孔隙度，单位为1，所述μ_f为岩心柱塞样的裂缝中的气体粘度，单位为mPa·s，所述C_gf为岩心柱塞样的裂缝中的体积压缩系数，单位为MPa^-1，所述T根据计算精度需要分成N个等份，t_D ⁱ代表第i个时间，而和分别代表第i个时间时和的值，所述M为根据计算精度需要将径向坐标轴上的岩心柱塞样基质的无因次半径等分的份数，所述L为根据计算精度需要将轴向坐标上的岩心柱塞样的无因次总长度等分的份数；

根据离散化的方程组求出对应的不同等份时间的理论的岩心夹持器进口端压力岩心夹持器出口端压力

根据所述t_D ⁱ、绘制理论曲线。

该方法通过记录岩心夹持器进口端和出口端的压力变化，针对单块岩心建立了考虑裂缝、基质渗流特征的数学模型，对一块具有裂缝的岩心进行一次测试就能分别测出裂缝、基质的渗透率各是多少，节约了时间，提高了测量效率。

附图说明

图1为本发明实施例提供的压力脉冲测量装置示意图。

图2为本发明实施例提供的绘制的实测曲线与理论曲线对比图。

图3为本发明实施例提供的岩心柱塞样径向和轴向上的坐标示意图。

具体实施方式

下面结合附图和实施例对本发明的技术方案进行详细描述。

参见附图1，一种压力脉冲测量装置，包括通过密封管线依次连接的气瓶1、压力调节阀2、气动阀Ⅰ3、气动阀Ⅱ5、压力传感器Ⅱ6、气动阀Ⅲ7、压力传感器Ⅲ8，还包括连接在压力传感器Ⅱ6和压力传感器Ⅲ8之间的岩心夹持器9、环压阀12、环压泵11，所述环压阀12控制环压泵11和岩心夹持器9的连通，所述压力传感器Ⅱ6和压力传感器Ⅲ8测量岩心夹持器9两端压力。还包括连接在气动阀Ⅰ3和气动阀Ⅱ5之间的压力传感器Ⅰ4、手动阀13和标准室14，所述手动阀13控制压力传感器Ⅰ4和标准室14的连通。所述连接在气动阀Ⅰ3和气动阀Ⅱ5之间的压力传感器Ⅰ4、手动阀13和标准室14用于测量岩心夹持器(9)进口端和出口端体积；另外调整脉冲压力时可以先关闭气动阀Ⅱ5，打开气动阀Ⅰ3，压力平衡时关闭气动阀Ⅰ3，再打开气动阀Ⅱ5，这样可以起到使压力缓冲平衡的作用，在这期间，通过压力传感器Ⅰ4可以看气动阀Ⅰ3到气动阀Ⅱ5之间的压力是否趋于稳定。还包括连接在岩心夹持器9和环压阀12之间的温度传感器10。

一种利用压力脉冲测量装置测量基质和裂缝渗透率的方法，包括如下步骤：在恒温条件下通过所述环压泵11提供压力使所述岩心夹持器9对岩心柱塞样12加环压p_h；

开启所述气动阀Ⅰ3、气动阀Ⅱ5、气动阀Ⅲ7，调节所述压力调节阀2为p₁稳定；

关闭所述气动阀Ⅱ5、气动阀Ⅲ7；

调节所述压力调节阀2为p₁’稳定，关闭所述气动阀Ⅰ3；

瞬间打开所述气动阀Ⅱ5，施加脉冲压力，通过所述压力传感器Ⅱ6、压力传感器Ⅲ8采集所述岩心夹持器9进口端和出口端的压力，直至两端压力平衡；

两端压力平衡后，所述压力传感器Ⅱ6、压力传感器Ⅲ8继续采集压力值直到无压降发生；

记录压力传感器Ⅱ6(记录岩心夹持器9进口端压力值)、压力传感器Ⅲ8(记录岩心夹持器9出口端压力值)记录的压力、时间，时间单位为秒，压力单位为MPa，见表1。

表1

将实测数据无因次化，用拟压力形式表示其中μ是气体的粘度，Z为气体偏差因子，p为压力传感器Ⅱ6或Ⅲ8记录的气体压力，p₀为压力传感器Ⅱ6或Ⅲ8记录的起始压力；Z值计算时取1；

通过公式将所述拟压力无因次化为无因次拟压力；

通过公式将传感器记录的时间t无因次化为无因次时间，所述代表t时刻岩心夹持9出口端或进口端的无因次压力值，所述代表0时刻岩心夹持器9出口端或进口端的拟压力值，所述是指在没加脉冲压力前岩心夹持器9夹持的岩心柱塞样12的拟压力值，所述t_D是无因次时间，所述k_m是岩心柱塞样12的基质渗透率，单位是mD，所述φ_m代表岩心柱塞样12的基质的孔隙度，单位是1，所述μ_m是岩心柱塞样12的基质的气体粘度，单位是mPa·s，所述C_g为气体体积压缩系数，数值上等于压力传感器Ⅱ6记录的压力平均值和压力传感器Ⅲ8记录的压力平均值之和的倒数，单位是MPa^-1，l为岩心柱塞样12的高，角标上D代表无因次量；所述岩心夹持器9进口端即岩心夹持器9靠近脉冲压力施加的一端，所述岩心夹持器9出口端即岩心夹持器9远离脉冲压力施加的一端。

将实测数据无因次化后的数据见表2。

表2

根据无因次压力和无因次时间t_D绘制实测曲线。

所述理论曲线通过如下方法绘制：

①通过公式计算无因次基质压力无因次裂缝压力

公式一：

公式二：

上式中 $L_1(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k+n}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]\mathrm{\λ}{\mathrm{\π}}^2{t}_D}{\mathrm{\ω}{\mathrm{\ξ}}_m}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$ $L_2(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]\mathrm{\λ}{\mathrm{\π}}^2{t}_D}{\mathrm{\ω}{\mathrm{\ξ}}_m}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$

所述x_D代表岩心柱塞样12轴向上的无因次坐标，所述l为岩心柱塞样的高，所述z_D代表岩心柱塞样12径向上的无因次坐标，所述t_D代表无因次时间，k＝1,2,3…∞，n＝1,2,3…∞，所述角标f代表裂缝渗透率相关数值，所述角标m代表基质渗透率相关数值，所述角标u代表岩心夹持器9进口端，所述角标d代表岩心夹持器9出口端，所述*为卷积运算符，所述代表岩心夹持器9进口端无因次拟压力，所述代表岩心夹持器9出口端无因次拟压力；

②再将代入下面的积分-微分方程组Ⅰ

以上式中， ${\mathrm{\ξ}}_f=\frac{\stackrel{\‾}{{\mathrm{\μ}}_f{C}_{\mathrm{tf}}}}{\stackrel{\‾}{{\mathrm{\μ}}_0{C}_{t0}}},{\mathrm{\ξ}}_m=\frac{{\mathrm{\μ}}_m{p}_0}{{\mathrm{\μ}}_0{p}_m},\mathrm{\λ}=\frac{k_m{h}_m}{k_f{h}_f},\mathrm{\ω}=\frac{{\mathrm{\φ}}_m{\stackrel{\‾}{p}}_f{h}_m}{{\mathrm{\φ}}_f{\stackrel{\‾}{p}}_m{h}_f},{\mathrm{\ξ}}_d=\frac{\stackrel{\‾}{{\mathrm{\μ}}_d}{p}_0}{{\mathrm{\μ}}_0\stackrel{\‾}{p_d}},$ 其中所述φ_m代表岩心柱塞样12的基质的孔隙度，单位为1，所述φ_f代表岩心柱塞样12的裂缝孔隙度，单位为1，所述C_tf代表岩心柱塞样12的裂缝13中气体的体积压缩系数，单位为MPa^-1，所述C_t0代表未加脉冲压力前岩心柱塞样12的基质中气体的体积压缩系数，单位为MPa^-1，所述p₀是未加脉冲压力前岩心夹持器9出口端的气体压力，单位为MPa，所述p_f是指岩心柱塞样12的裂缝13中气体的压力，单位为MPa，所述p_m是岩心柱塞样12的基质中气体的压力，单位为MPa，所述p_d是岩心夹持器9出口端气体压力，单位为MPa，所述h_f是岩心柱塞样12的裂缝13的半长，所述h_m是岩心柱塞样12的半径减去裂缝13的半长，所述μ₀是未加脉冲压力前岩心柱塞样12中的气体的粘度，单位是mPa·s，所述μ_d是岩心夹持器9出口端气体粘度，单位是mPa·s，所述μ_m是岩心柱塞样12的基质内气体粘度，单位是mPa·s，所述μ_f是岩心柱塞样12的裂缝13内的气体粘度，单位是mPa·s，所述x是圆柱体样品岩心高度方向上的坐标值，所述z代表样品岩心直径方向上的坐标值，所述l为岩心柱塞样12的高，所述V_f为岩心柱塞样12的裂缝13的体积，所述V_u岩心夹持器9进口端体积，所述V_d岩心夹持器9出口端体积；所述岩心夹持器9进口端体积即气动阀Ⅱ5到气动阀Ⅲ7之间管线的体积加上压力传感器Ⅱ6到岩心夹持器9进口端之间管线的体积，所述岩心夹持器出口端体积即气动阀Ⅲ7到压力传感器Ⅲ8之间管线的体积加上压力传感器Ⅲ8到岩心夹持器9出口端之间管线的体积；

将带有所述为变量的积分-微分方程组Ⅰ离散化，得到以下方程组：

其中，

k＝1，2，3…M，j＝1，2，3…L，

把实验总共所需时间t代入公式所得到的无因次时间t_D值就是T，式中所述k_m是岩心柱塞样的基质渗透率，单位是mD，所述k_f为岩心柱塞样裂缝渗透率，单位是mD，所述h_f是岩心柱塞样12的裂缝13的半长，所述h_m是岩心柱塞样的半径减去裂缝13的半长，所述φ_f代表岩心柱塞样12的裂缝孔隙度，单位为1，所述μ_f为岩心柱塞样12的裂缝13的气体粘度，单位为mPa·s，所述C_gf为岩心柱塞样12的裂缝13中的体积压缩系数，单位为MPa^-1，所述T根据计算精度需要分成N个等份，t_D ⁱ代表第i个时间，i取值到为0为止，而和分别代表第i个时间时和的值，所述M为根据计算精度需要将径向坐标轴上的柱塞样12的基质的无因次半径等分的份数，所述L为根据计算精度需要将轴向坐标上的柱塞样12的无因次总长度等分的份数；根据方程组求出对应的不同等份时间的理论的岩心夹持器9进口端压力岩心夹持器9出口端压力根据t_D ⁱ、绘制理论曲线。

实施例中选取k_m值为0.0001mD、0.00001mD和0.000001mD，k_f值为0.00312mD。绘制理论曲线所需值见表3、表4：

表3

表4

将实测曲线和理论曲线作对比，参见附图2，拟合好的曲线为以三角形点和叉形点表示的曲线，因此可以求得岩心裂缝渗透率为0.00312mD，基质渗透率为0.00001mD。

本发明只需对一块具有裂缝的岩心柱塞样进行一次测试就能分别测出裂缝、基质的渗透率，提高了测量效率。

最后所应说明的是，以上实施例仅用以说明本材料的技术实施方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (3)

1.一种压力脉冲测量装置测量基质和裂缝渗透率的方法，所述压力脉冲测量装置包括通过密封管线依次连接的气瓶(1)、压力调节阀(2)、气动阀Ⅰ(3)、气动阀Ⅱ(5)、压力传感器Ⅱ(6)、气动阀Ⅲ(7)、压力传感器Ⅲ(8)，还包括连接在压力传感器Ⅱ(6)和压力传感器Ⅲ(8)之间的岩心夹持器(9)、环压阀(12)、环压泵(11)，所述环压阀(12)控制环压泵(11)和岩心夹持器(9)的连通，所述压力传感器Ⅱ(6)和压力传感器Ⅲ(8)测量岩心夹持器(9)两端压力；还包括连接在气动阀Ⅰ(3)和气动阀Ⅱ(5)之间的压力传感器Ⅰ(4)、手动阀(13)和标准室(14)，所述手动阀(13)控制压力传感器Ⅰ(4)和标准室(14)的连通；还包括连接在岩心夹持器(9)和环压阀(12)之间的温度传感器(10)；

其特征在于，所述测量基质和裂缝渗透率的方法，包括如下步骤：

在恒温条件下通过所述环压泵(11)提供压力使所述岩心夹持器(9)对岩心柱塞样加环压p_h；

开启所述气动阀Ⅰ(3)、气动阀Ⅱ(5)、气动阀Ⅲ(7)，调节所述压力调节阀(2)为p₁稳定；

关闭所述气动阀Ⅱ(5)、气动阀Ⅲ(7)；

调节所述压力调节阀(2)为p₁’稳定，关闭所述气动阀Ⅰ(3)；

瞬间打开所述气动阀Ⅱ(5)，施加脉冲压力，通过所述压力传感器Ⅱ(6)、压力传感器Ⅲ(8)采集所述岩心夹持器(9)进口端和出口端的压力，直至两端压力平衡；

两端压力平衡后，所述压力传感器Ⅱ(6)、压力传感器Ⅲ(8)继续采集压力值直到无压降发生；

将所述压力传感器Ⅱ(6)、压力传感器Ⅲ(8)记录的压力、时间无因次化并绘制压力-时间关系的实测曲线，根据理论计算的所述岩心夹持器(9)进口端和出口端的无因次压力和无因次时间绘制理论曲线，将实测曲线与理论曲线做对比，拟合好的渗透率就是该岩心样品的对应的基质渗透率k_m和裂缝渗透率k_f。

2.如权利要求1所述的测量基质和裂缝渗透率的方法，其特征在于：所述实测曲线通过如下方法绘制：

用拟压力形式表示其中μ是气体粘度，单位是mPa·s，Z为气体偏差因子，单位是1，p为压力传感器Ⅱ(6)或Ⅲ(8)记录的实测气体压力，单位为MPa，p₀为压力传感器Ⅱ(6)或Ⅲ(8)记录的起始压力，单位为MPa；

通过公式将所述拟压力无因次化为无因次压力；

通过公式将传感器记录的时间t无因次化为无因次时间，式中，所述代表t时刻岩心夹持器(9)出口端或进口端的无因次压力值，所述代表0时刻岩心夹持器(9)出口端或进口端的拟压力值，所述是指在没加脉冲压力前岩心夹持器(9)夹持的岩心柱塞样的拟压力值，所述t_D是无因次时间，所述k_m是岩心柱塞样基质渗透率，单位是mD，所述φ_m代表岩心柱塞样基质的孔隙度，单位是1，所述μ_m是岩心柱塞样基质的气体粘度，单位是mPa·s，所述C_g为气体体积压缩系数，数值上等于压力传感器Ⅱ(6)记录的压力平均值和压力传感器Ⅲ(8)记录的压力平均值之和的倒数，单位是MPa^-1，所述l为岩心柱塞样的高，角标上D代表无因次量；所述岩心夹持器(9)进口端即岩心夹持器(9)靠近脉冲压力施加的一端，所述岩心夹持器(9)出口端即岩心夹持器(9)远离脉冲压力施加的一端；

根据所述无因次拟压力和无因次时间t_D绘制实测曲线。

3.如权利要求2所述的测量基质和裂缝渗透率的方法，其特征在于：所述理论曲线通过如下方法绘制：

①通过公式分别计算无因次基质压力无因次裂缝压力

公式一：

公式二：

上式中 $L_1(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k+n}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]{\mathrm{\λ\π}}^2{t}_D}{\mathrm{\ω\ξm}}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$ $L_2(x_D,z_D,t_D)=\underset{n=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{k=1}{\overset{\∞}{\mathrm{\Σ}}}\frac{{(-1)}^{k}n}{2k+1}{e}^{-\frac{[{(k+1)}^2+n^2]{\mathrm{\λ\π}}^2{t}_D}{\mathrm{\ω\ξm}}}\cos (k+\frac{1}{2})\mathrm{\π}{z}_D\sin \left(\mathrm{n\π}{x}_D\right),$ 所述x_D代表岩心柱塞样轴向上的无因次坐标，所述l为岩心柱塞样的高，所述z_D代表岩心柱塞样径向上的无因次坐标，所述所述t_D代表无因次时间，k＝1,2,3…∞，n＝1,2,3…∞，所述角标f代表裂缝渗透率相关数值，所述角标m代表基质渗透率相关数值，所述角标u代表岩心夹持器(9)进口端，角标d代表岩心夹持器(9)出口端，*为卷积运算符，代表岩心夹持器(9)进口端无因次拟压力，代表岩心夹持器(9)出口端无因次拟压力；

②再将所述代入下面的积分-微分方程组Ⅰ

以上式中， ${\mathrm{\ξ}}_f=\frac{\stackrel{\‾}{{\mathrm{\μ}}_f{C}_{\mathrm{tf}}}}{\stackrel{\‾}{{\mathrm{\μ}}_0{C}_{t0}}},{\mathrm{\ξ}}_m=\frac{{\mathrm{\μ}}_m{p}_0}{{\mathrm{\μ}}_0{p}_m},\mathrm{\λ}=\frac{k_m{h}_m}{k_f{h}_f},\mathrm{\ω}=\frac{{\mathrm{\φ}}_m{\stackrel{\‾}{p}}_f{h}_m}{{\mathrm{\φ}}_f{\stackrel{\‾}{p}}_m{h}_f},{\mathrm{\ξ}}_d=\frac{\stackrel{\‾}{{\mathrm{\μ}}_d}{p}_0}{{\mathrm{\μ}}_0\stackrel{\‾}{p_d}},$ 其中所述φ_m代表岩心柱塞样基质的孔隙度，单位为1，所述φ_f代表岩心柱塞样的裂缝孔隙度，单位为1，所述C_tf代表岩心柱塞样的裂缝(13)中气体的体积压缩系数，单位为MPa^-1，所述C_t0代表未加脉冲压力前岩心柱塞样基质中气体的体积压缩系数，单位为MPa^-1，所述p₀是未加脉冲压力前岩心夹持器(9)出口端的气体压力，单位为MPa，所述p_f是指岩心柱塞样的裂缝(13)中气体的压力，单位为MPa，所述p_m是岩心柱塞样基质中气体的压力，单位为MPa，所述p_d是岩心夹持器(9)出口端气体压力，单位为MPa，所述h_f是岩心柱塞样的裂缝(13)的半长，所述h_m是岩心柱塞样的半径减去裂缝(13)的半长，所述μ₀是未加脉冲压力前岩心柱塞样中的气体的粘度，单位是mPa·s，所述μ_d是岩心夹持器(9)出口端气体粘度，单位是mPa·s，所述μ_m是岩心柱塞样基质内气体粘度，单位是mPa·s，所述μ_f是岩心柱塞样的裂缝(13)内的气体粘度，单位是mPa·s，所述x是圆柱体样品岩心高度方向上的坐标值，所述z代表样品岩心直径方向上的坐标值，l为岩心柱塞样的高，所述V_f为岩心柱塞样的裂缝(13)的体积，所述V_u岩心夹持器(9)进口端体积，所述V_d岩心夹持器(9)出口端体积；所述岩心夹持器(9)进口端体积即气动阀Ⅱ(5)到气动阀Ⅲ(7)之间管线的体积加上压力传感器Ⅱ(6)到岩心夹持器(9)进口端之间管线的体积，所述岩心夹持器出口端体积即所述气动阀Ⅲ(7)到所述压力传感器Ⅲ(8)之间管线的体积加上所述压力传感器Ⅲ(8)到所述岩心夹持器(9)出口端之间管线的体积；

将带有所述为变量的积分-微分方程组Ⅰ离散化，得到以下方程组：

其中，

k＝1，2，3…M，j＝1，2，3…L，

把实验总共所需时间t代入公式所得到的无因次时间t_D值就是T，式中所述k_m是岩心柱塞样基质渗透率，单位是mD，所述k_f为岩心柱塞样的裂缝渗透率，单位是mD，所述h_f是岩心柱塞样的裂缝(13)的半长，所述h_m是岩心柱塞样的半径减去裂缝(13)的半长，所述φ_f代表岩心柱塞样的裂缝孔隙度，单位为1，所述μ_f为岩心柱塞样的裂缝(13)中的气体粘度，单位为mPa·s，所述C_gf为岩心柱塞样的裂缝(13)中的体积压缩系数，单位为MPa^-1，所述T根据计算精度需要分成N个等份，t_D ⁱ代表第i个时间，而和分别代表第i个时间时和的值，所述M为根据计算精度需要将径向坐标轴上的岩心柱塞样基质的无因次半径等分的份数，所述L为根据计算精度需要将轴向坐标上的岩心柱塞样的无因次总长度等分的份数；

根据离散化的方程组求出对应的不同等份时间的理论的岩心夹持器(9)进口端压力岩心夹持器(9)出口端压力

根据所述t_D ⁱ、绘制理论曲线。