[0048] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0049] 本发明的目的是提供一种多层合采气井产量劈分方法及系统，能够准确有效的确定气井产量的劈分结果。
[0050] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0051] 如图1所示，本发明提供的一种多层合采气井产量劈分方法，包括以下步骤：
[0052] 步骤101：根据单井气层组分年产出剖面测试资料与测井解释参数，建立气层组分年分层产量贡献率与测井解释参数的综合关系。
[0053] 气层组分年分层产量贡献率与测井解释参数的综合关系建立过程如下：
[0054] (1)筛选测井解释参数。
[0055] 测井解释参数包括有利参数和不利参数。
[0056] 所选取的有利参数为反映储层渗流能力的地层系数(Kh)、反映气层含气能力的含气饱和度(Sg)，不利参数为岩性参数‑泥质含量(Vsh)。
[0057] (2)根据所述测井解释参数，计算每层含气小层参数。
[0058] 对于地层系数、含气饱和度而言，其值越大，含气小层性能越好，分层产量贡献越大。用本项分层参数值除以层组分层参数值之和，即：
[0059]
[0060] 式中
[0061] Kh—地层系数，毫达西米mD.m；
[0062] X—参数，无因次，小数；
[0063] 下标i—第i层；
[0064] n—合采总层数；
[0065]
[0066] 对于泥质含量(Vsh)而言，其值越小，含气小层性能越好，分层产量贡献越大。用本项参数值的倒数除以层组分层参数值倒数之和，即：
[0067]
[0068] 式中
[0069] Vsh—泥质含量，无因次，小数；
[0070] Y—参数，无因次，小数。
[0071] (3)根据每层所述含气小层参数，构建气层组分年分层产量贡献率与测井解释参数的综合关系。
[0072] 给上述三项测井解释参数分别设置不同权重系数x、y、z，得到综合系数F：
[0073] Fi＝Xi·x+Yi·y+Zi·z(4)；
[0074] 其中：Fi为第i层的综合系数；x+y+z＝1。
[0075] 步骤102：根据单井气层组分年产出剖面测试资料，计算气层组分年分层历史平均产量贡献率。计算公式如下：
[0076]
[0077] 式中
[0078] —气层组分年分层历史平均产量贡献率；
[0079] Vg—气层组分年分层历史产量贡献率，无因次，小数；
[0080] Z—参数，无因次，小数。
[0081] 步骤103：根据所述气层组分年分层历史平均产量贡献率，以及所述气层组分年分层产量贡献率与测井解释参数的综合关系，采用多元线性回归算法，拟合确定所述综合关系中各个测井解释参数的权重系数，进而确定最终的气层组分年分层产量贡献率与测井解释参数的综合关系。具体为：
[0082] 采用多元线性回归算法，将公式(4)与公式(5)进行多参数线性拟合，观察不同权重系数的拟合效果，然后不断调整权重系数，选取选取拟合结果为综合系数最大时所设置的权重系数，得到最终的气层组分年分层产量贡献率与测井解释参数的综合关系，其公式：
[0083] 式中：a，b—线性拟合系数。
[0084] 步骤104：根据所述最终的气层组分年分层产量贡献率与测井解释参数的综合关系，反推确定气层组分年分层产量贡献率计算关系式。
[0085] 其关系式为：
[0086] 式中： —气层组分年分层产量贡献率。
[0087] 步骤105：根据单井气层组分年产出剖面测试资料和反推确定的所述气层组分年分层产量贡献率计算关系式，对单井产量进行劈分。
[0088] 根据单井的具体测试和生产情况进行劈分，如图2所示，具体包括：
[0089] (1)当单井投产后一直未测产出剖面测试资料时，按所建立的气层组分年分层产量贡献率计算关系式计算气层组分年分层产量贡献率，然后按照气层组分年分层产量贡献率对单井产量进行劈分。
[0090] 特殊情况，如果综合关系无法建立，那么气层组分年分层产量贡献率计算关系式也无法建立，此时就按统计的气层组分年分层历史平均产量贡献率对单井产量进行劈分。如果计算或统计的气层组分年分层历史平均产量贡献率之和不等于1时，首先进行归一化处理后，再进行劈分，详细见图2。
[0091] (2)当单井投产后只有一次产出剖面测试资料时，按所述气层组分年分层产量贡献率计算关系式和产出剖面测试日期之前每年单井产量计算气层组分年分层产量贡献率，然后按照所述气层组分年分层产量贡献率对单井产量进行劈分。
[0092] 特殊情况，如果综合关系无法建立，那么气层组分年分层产量贡献率计算关系式也无法建立，就按统计的气层组分年分层历史平均产量贡献率进行劈分，以后每年用产出剖面测试资料进行劈分。
[0093] (3)当单井投产后有多次产出剖面测试资料时，按如下三种测试情况进行劈分：
[0094] 第一种测试情况，当单井在投产开始到第一次产出剖面测试之前时，按所述气层组分年分层产量贡献率计算关系式和第一次产出剖面测试资料计算气层组分年分层产量贡献率，然后按照所述气层组分年分层产量贡献率对单井产量进行劈分。
[0095] 第二种实际情况，当单井在多次产出剖面测试之间时，根据邻近两次产出剖面测试资料，采用线性插值方式得到过程产出剖面测试资料，然后按照所述气层组分年分层产量贡献率计算关系式和所述过程产出剖面测试资料计算气层组分年分层产量贡献率，最后按照所述气层组分年分层产量贡献率对单井产量进行劈分。
[0096] 第三种实际情况，当单井在最后一次产出剖面测试至今时，按照所述气层组分年分层产量贡献率计算关系式和最后一次产出剖面测试资料计算气层组分年分层产量贡献率，然后按照所述气层组分年分层产量贡献率对单井产量进行劈分。
[0097] 为实现上述目的，本发明还提供了一种多层合采气井产量劈分系统，如图3所示，包括：
[0098] 综合关系建立模块201，用于根据单井气层组分年产出剖面测试资料与测井解释参数，建立气层组分年分层产量贡献率与测井解释参数的综合关系。
[0099] 平均产量贡献率计算模块202，用于根据单井气层组分年产出剖面测试资料，计算气层组分年分层历史平均产量贡献率。
[0100] 综合关系确定模块203，用于根据所述气层组分年分层历史平均产量贡献率，以及所述气层组分年分层产量贡献率与测井解释参数的综合关系，采用多元线性回归算法，拟合确定所述综合关系中各个测井解释参数的权重系数，进而确定最终的气层组分年分层产量贡献率与测井解释参数的综合关系。
[0101] 反推模块204，用于根据所述最终的气层组分年分层产量贡献率与测井解释参数的综合关系，反推确定气层组分年分层产量贡献率计算关系式。
[0102] 第一劈分模块205，用于根据单井气层组分年产出剖面测试资料和反推确定的所述气层组分年分层产量贡献率计算关系式，对单井产量进行劈分。
[0103] 其中，所述综合关系建立模块201，具体包括：
[0104] 筛选单元，用于筛选测井解释参数；所述测井解释参数包括泥质含量、反映储层渗流能力的地层系数以及反映气层含气能力的含气饱和度。
[0105] 计算单元，用于根据所述测井解释参数，计算每层含气小层参数。
[0106] 综合关系建立单元，用于根据每层所述含气小层参数，构建气层组分年分层产量贡献率与测井解释参数的综合关系。
[0107] 所述综合关系确定模块203，具体包括：
[0108] 所述综合关系确定单元，用于根据所述气层组分年分层历史平均产量贡献率，以及所述气层组分年分层产量贡献率与测井解释参数的综合关系，采用多元线性回归算法，拟合所述综合关系中各个测井解释参数的权重系数，并选取拟合结果为综合系数最大时所设置的权重系数，进而确定最终的气层组分年分层产量贡献率与测井解释参数的综合关系。
[0109] 优选地，所述多层合采气井产量劈分系统还包括：第二劈分模块，用于当所述最终的气层组分年分层产量贡献率与测井解释参数的综合关系无法建立时，按照所述气层组分年分层历史平均产量贡献率对单井产量进行劈分。
[0110] 本发明提供的一种多层合采气井产量劈分方法及系统，用于确定多层合采气藏各个小层的产量贡献率。该产量劈分方法准确度明显高于常规方法，特别是对于泥质含量和含水饱和度较高的气藏更具有适应性；此外，不需要关井，不影响气田正常生产，在有无产气剖面测试的基础上都可以进行，节约了测试成本，能适应和满足复杂多层气藏开发的需要。
[0111] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0112] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
