方法考虑了因水排出引起的煤基质收缩的应变机理，这一机理对煤储层动态渗透率的jingque反演起到了关键作用，将因水排出引起的煤基质收缩的应变机理和渗透率与压力的经验模型相结合，实现排采过程中煤储层动态渗透率的预判，进而实现排采过程中煤地层动态渗透率变化的快速、准确监测；

②本发明所述的方法不仅可以预判煤储层动态渗透率的变化规律，还可以反演确定煤基质孔隙不同含水条件对储层动态渗透率的影响程度，为实现煤层气井的合理开发奠定了理论支撑。

上述说明仅是本申请技术方案的概述，为了能够更清楚了解本申请的技术手段，而可依照说明书的内容予以实施，并且为了让本申请的上述和其它目的、特征和优点能够更明显易懂，以下特举本申请的具体实施方式。

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

本发明实施例公开一种计算煤层气不同生产阶段气水渗透率的方法，该方法包括：

根据开采阶段的不同，将生产阶段分为以下四个阶段：单相水流动阶段、近井地带非饱和单相流动阶段、近井地带气水两相流动阶段、单相气流动阶段；如图1所示，排采初期，煤层主要产水，同时伴随有少量游离气、溶解气产出；当煤层降至临界解吸压力以下时，煤层甲烷分子迅速解吸，然后扩散到裂隙中，表现为产气量逐渐增大，产水量逐渐减小；随着采出水量的增加、生产压差的进一步增大，煤层中含水饱和度相对降低，变为以产气为主，并逐渐达到产气高峰，产水量则相对稳定在一个较低的水平上。

根据有效应力效应、因气体解吸引起的煤基质收缩效应、因水解吸引起的煤基质收缩效应、克林肯伯格效应对渗透率的制约，结合渗透率与压力的经验模型计算渗透率。在煤层气排采初期单相流阶段，煤储层所受有效应力不断增大，使得裂隙宽度变窄，渗透率降低；当储层压力降到临界解吸压力以下时，煤层气开始解吸，煤基质收缩效应逐渐加强，使得裂隙变宽，渗透率出现反弹；在排采后期，储层压力已降至较低水平，低压条件下气体的克林肯伯格效应更加明显，有利于改善煤储层渗透率。本发明结合渗透率与压力的经验模型来预测渗透率变化，综合考虑排采过程中煤储层物性受应力敏感性、因气体解吸引起的基质收缩、因水解吸引起的基质收缩和克林肯伯格效应4种效应的制约，实现排采过程中储层动态渗透率的预测。

在一个实施例中，计算煤层气不同生产阶段气水渗透率包括：

煤储层在单相水流动阶段时煤层气井气水渗透率的计算，和/或煤储层在近井地带非饱和单相流动阶段时煤层气井气水渗透率的计算，和/或煤储层在近井地带气水两相流动阶段时煤层气井气水渗透率的计算，和/或煤储层在单相气流动阶段时煤层气井气水渗透率的计算。

在一个优选实施例中，煤储层在单相水流动阶段时煤层气井气水渗透率的计算包括以下步骤：

步骤301、构建应力敏感对煤储层影响的应力表征模型：

式(1)中：σ是有效水平应力，σ0是原始条件下有效水平应力，ν是泊松比，p是储层压力，p0是原始储层压力；

步骤302、建立因气体解吸引起的煤基质收缩下的应变量模型：

式(2)中：是非平衡条件下因气体吸附引起的煤基质膨胀量，是煤基质内部气体压力与外部气体压力平衡时对应的气体吸附引起的煤基质膨胀量，t是时间，τ是扩散时间；

由于煤基质内部气体压力与外部气体压力平衡时对应的气体吸附引起的煤基质膨胀量是时间变量，无法求解，但可以作如下假设：①煤样为均质体；②煤样的外部压力恒定。则如下等式成立，进而实现模型简化：

式(3)中：εl是朗格缪尔膨胀系数，pl是压力系数，p∞是Zui终压力值；

那么，因气体解吸引起的煤基质收缩膨胀量模型可简化为：

式(4)中：是非平衡条件下因气体吸附引起的煤基质膨胀量，是原始气体压力条件下因气体吸附引起的煤基质膨胀量，是Zui终气体压力条件下因气体吸附引起的煤基质膨胀量，t是时间，τ是扩散时间；

步骤303、建立煤储层在单相水流动阶段下的煤层气井气水渗透率表征模型：

式(5)中：k是渗透率，k0是初始渗透率，cf是煤岩割理的压缩系数，ν是泊松比，p是储层压力，p0是原始储层压力，e是煤岩弹性模量，是原始气体压力条件下因气体吸附引起的煤基质膨胀量，τ是扩散时间，t是时间。

日本TOYOOKI丰兴电磁阀型号 AD-SL222-303D, AD-SL222-304D，AD-SL222-506D, AD-SL222-508D，AD-SL222-810D, AD-SL222-812D, AD-SL222-916D， AD-SL-231-101D,AD-SL231-102D,AD-SL231-303D,AD-SL231-304D,AD-SL231-506D,AD-SL231-508D AD-SL231-810D, AD-SL231-812D,AD-SL231-916D AD-SL041-303D,AD-SL041-304D,AD-SL041-506D,AD-SL041-508D,AD-SL041-810D,AD-SL041-812D AD-SL041-916D AD-SL122(E)-303B，AD-SL122(E)-304B，AD-SL122(E)-506B，AD-SL122(E)-508B，AD-SL122(E)-810B，AD-SL122(E)-812B， AD-SL122(E)-916B，AD-SL122(E)-303B，AD-SL222(E)-304B，AD-SL222(E)-506B，AD-SL222(E)-508B，AD-SL222(E)-916B，AD-SL222(E)-810B，AD-SL222(E)-812B

如图2和图3所示，在投产初期，储层各点压力均高于临界解吸压力，气体主要以吸附态赋存于煤基质块中。初始状态下的割理中仅存在少量自由气，可近似认为被水饱和。同时，由于煤层气藏属于低压气藏，导致水中的溶解气含量极少，可被忽略。气井初投产时，割理被单相水饱和，故该阶段下的流动特征为单相水流，即煤层气排水降压阶段。该阶段的主要目的是降低储层压力，促进煤岩基质块中的吸附气发生解吸并进入割理，形成产气。在该生产阶段，水相相对渗透率始终为1。煤岩易受应力敏感影响，即煤储层的渗透率会随着储层压力降低而降低，故该阶段中，全区的水相相对渗透率始终为1，近井区域的水相有效渗透率要低于远井区域。

在一个实施例中，煤储层在近井地带非饱和单相流动阶段时煤层气井气水渗透率的计算包括以下步骤：

步骤401、构建应力敏感对煤储层影响的应力表征模型：

式(6)中，σ是有效水平应力，σ0是原始条件下有效水平应力，ν是泊松比，p是储层压力，p0是原始储层压力；

步骤402、建立因气体解吸引起的煤基质收缩下的应变量模型：

式(7)中，是非平衡条件下因气体吸附引起的煤基质膨胀量，是煤基质内部气体压力与外部气体压力平衡时对应的气体吸附引起的煤基质膨胀量，t是时间，τ是扩散时间；

由于煤基质内部气体压力与外部气体压力平衡时对应的气体吸附引起的煤基质膨胀量是时间变量，无法求解，但可以作如下假设：①煤样为均质体；②煤样的外部压力恒定。则如下等式成立，进而实现模型简化：

式(8)中，εl是朗格缪尔膨胀系数，pl是压力系数，p∞是Zui终压力值；

那么，因气体解吸引起的煤基质收缩膨胀量模型可简化为：

式(9)中，是非平衡条件下因气体吸附引起的煤基质膨胀量，是原始气体压力条件下因气体吸附引起的煤基质膨胀量，是Zui终气体压力条件下因气体吸附引起的煤基质膨胀量，t是时间，τ是扩散时间；