一种底水气藏水侵识别及产水阶段划分方法与流程

1.本发明涉及石油和天然气开发技术领域，特别涉及一种底水气藏水侵识别及产水阶段划 分方法。


背景技术：

2.目前，在大量底水气藏水侵机理的理论和实验研究基础上，针对底水气藏早期水侵识别， 国内外学者采用常规性的方法包括视地层压力法、视地质储量法、水侵体积系数法、产量递 减曲线判别法、水锥极限产量法、试井监测识别法，这些方法要求水侵量计算准确或是需要 气井生产较长时间，储层非均质性的存在和水气分布关系的复杂性使常规水侵识别方法在实 际应用中受到较大影响，前瞻性和适用性较差。
3.现阶段，关于气井采气阶段划分和水侵阶段识别的相关研究较少，前人多采用水气比和 产水量(或产气量)关系识别水侵阶段，这些方法相对单一，缺乏对比性，尚需丰富完善。


技术实现要素：

4.本技术的目的在于克服现有技术中所存在的上述不足，提供一种在水侵初期即能识别出 气藏水侵状态的底水气藏水侵识别方法。
5.本技术还提供了一种底水气藏产水阶段划分方法。
6.为了实现上述发明目的，本技术提供了以下技术方案：一种底水气藏水侵识别方法，包 括以下步骤：
7.采集气井产出液在气井水侵各阶段的水化学系数，并根据气井水侵阶段划分水化学系数 的分布范围；
8.根据所述水化学系数的分布范围筛选气井水侵识别系数；
9.采集气井产出液在气井水侵各阶段的不同离子的离子浓度，并根据气井水侵阶段划分不 同离子的离子浓度分布范围；
10.根据所述离子浓度的分布范围筛选气井水侵识别离子；
11.根据气井水侵识别系数和气井水侵识别离子绘制水侵识别相图；
12.根据所述水侵识别相图判断待测气井的水侵阶段。
13.进一步地，所述水化学系数包括钾钠氯系数、镁钙系数、氯镁系数、碳酸盐平衡系数和 脱硫系数。所述钠钾氯系数为气井采出液中钾离子和钠离子浓度之和与氯离子浓度的比值； 所述镁钙系数为气井采出液中镁离子浓度与该离子浓度的比值；所述氯镁系数为气井采出液 中氯离子浓度与镁离子浓度的比值；所述碳酸盐平衡系数为气井采出液中碳酸根离子浓度和 碳酸氢根离子浓度之和与钙离子浓度的比值；所述脱硫系数为气井采出液中硫酸根离子浓度 与氯离子浓度的比值*100。
14.进一步地，所述气井水侵阶段根据所述气井产出液的组成而划分为无水采气期、水侵初 期和带水采气期。所述无水采气期的产出液主要为天然气因温度和压力变化而产
出的凝析液； 所述水侵初期的产出液主要为凝析液和地层水的混合液；所述带水采气期的产出液主要为地 层水。
15.进一步地，根据各个气井水侵阶段中每一水化学系数的分布范围的重叠范围来筛选气井 水侵识别系数；当一水化学系数的分布范围在各个气井水侵阶段除端点以外无重叠部分，则 该水化学系数被筛选为一气井水侵识别系数；当一水化学系数的分布范围在各个气井水侵阶 段有端点以外的重叠部分，则该水化学系数不作为气井水侵识别系数。
16.即在气井水侵阶段，气井产出液中的每一水化学系数都随着气井产出液的变化而变化， 但每一水化学系数在不同气井水侵阶段的变化并不相同，因此对于每一水化学系数在不同气 井水侵阶段的分布范围需要单独划分。部分水化学系数在不同气井水侵阶段的分布范围可能 有较大的重叠区域，因此采用这些水化学系数难以判断气井所属的生产阶段，需要筛除这一 类水化学系数。
17.进一步地，所述离子浓度包括钾钠离子浓度、钙离子浓度、镁离子浓度、氯离子浓度和 碳酸氢根离子浓度。所述钾钠离子浓度为所述气井采出液中钠离子浓度和钾离子浓度之和。
18.进一步地，根据各个气井水侵阶段中每一离子浓度的分布范围的重叠范围来筛选气井水 侵识别离子；当一离子浓度的分布范围在各个气井水侵阶段除端点以外无重叠部分，则该离 子被筛选为一气井水侵识别系数；当一离子浓度的分布范围在各个气井水侵阶段有端点以外 的重叠部分，则该离子不作为气井水侵识别离子。
19.所述气井水侵识别离子的筛选与所述气井水侵识别系数的筛选原理基本一致，故不再一 一陈述。
20.进一步地，所述水侵识别相图以水侵识别系数和水侵识别离子分别作为x轴，以y轴确 定每一x轴的零点；再连接每一所述气井水侵阶段的水侵识别系数和水侵识别离子的端点值， 每一所述气井水侵阶段的水侵识别系数和水侵识别离子的端点值连线围成的区域即为对应的 气井水侵阶段。
21.所述水侵识别相图以每一个水侵识别系数和每一个水侵识别离子的都分别作为x轴，然 后均匀的分布在y轴的两边；选取每一个水侵识别系数和每一个水侵识别离子在每一气井水 侵阶段对应x轴上的端点值(两端端点值：最大值或最小值)；然后用直线连接每一气井水 侵阶段每一水侵识别系数的最大值连接，用直线连接每一气井水侵阶段每一水侵识别系数的 最小值连接，二者圈出的区域即为每一气井水侵阶段的水侵识别相图中的分布。
22.应用时，监测相同气藏中待测气井的采出液中水侵识别系数和水侵识别离子的数据，即 可根据水侵识别相图获得待测气井出现水侵的最早时间，提前识别出气井水侵，为气井生产 计划调整提供数据。此外，还可以根据待测气井的采出液中水侵识别系数和水侵识别离子的 数据与水侵识别相图中无水采气期和水侵初期的临界值之间的差值，提前预判待测气井是否 临近出现水侵，为气田工作人员提供提前预警，便于工作人员提前作出气井生产计划的调整 准备。
23.本技术还公开了一种底水气藏产水阶段划分方法，包括以下步骤：
24.采用上述底水气藏水侵识别方法制备待划分气井所在气藏的标准水侵识别相图；
25.采集待划分气井产出液的气井水侵识别系数和气井水侵识别离子，并绘制出气井
的实时 水侵识别线图；
26.将气井的实时水侵识别线图垂直投影在所述标准水侵识别相图中，根据所述气井的实时 水侵识别线图落在所述标准水侵识别相图中的区域，判断待划分气井所在的采气阶段。
27.本技术还公开了一种底水气藏产水阶段划分方法，包括以下步骤：采用上述底水气藏水 侵识别方法制备待划分气井所在气藏的标准水侵识别相图；
28.采集待划分气井产出液的气井水侵识别系数和气井水侵识别离子；
29.将采集到的气井水侵识别系数和气井水侵识别离子的数据代入所述标准水侵识别相图中， 判断待划分气井所在的采气阶段。
30.进一步地，上述两种底水气藏产水阶段划分方法可以联用，即先通过将采集到的待划分 气井产出液的气井水侵识别系数和气井水侵识别离子初步得到气井所在的采气阶段，然后再 通过采集到的待划分气井产出液的气井水侵识别系数和气井水侵识别离子绘制的实时水侵识 别线图进一步确定待划分气井所在的采气阶段，对气井所在采气阶段的划分结果做一次印证。
31.与现有技术相比，本发明的具有以下有益效果：本发明公开的底水气藏水侵识别方法， 筛选出可准确表征气井在不同生产阶段的水化学系数和离子来作为气井水侵识别系数和气井 水侵识别离子，并根据气井水侵识别系数和气井水侵识别离子作出水侵识别相图；通过水侵 识别相图来判断气井所在的气井水侵阶段，增加了气井水侵阶段判断的因素，提高了气井水 侵阶段识别的准确性；同时，本技术公开的水侵识别方法通过采出液中的离子变化来判断气 井所在的水侵阶段，其相较于通过气井采出液中的水气比和产水量来判断，更具有前瞻性， 更加灵敏，便于气井水侵初期的识别；此外，所采用的数据均为气井采出液中各离子的浓度， 通过实时监控气井采出液的水质即可得到，无需另外设置数据采集装备。本公开的底水气藏 产水阶段划分方法通过气井的气井水侵识别系数、气井水侵识别离子中的一个或多个在本申 请公开的水侵识别相图中的位置来进行气井产水阶段划分，其划分简单快捷，数据采集量少， 便于气井水侵阶段的初步判断、快速判断。本技术公开的另一底水气藏产水阶段划分方法通 过气井的气井水侵识别系数、气井水侵识别离子绘制的水侵识别线图在本技术公开的水侵识 别相图上垂直投影来进行气井产水阶段划分，更加准确，便于了解气井的水侵情况，为气井 的生产计划的制定赢得了更长的窗口期。
附图说明
32.图1本发明公开的底水气藏水侵识别方法的流程示意图；
33.图2本发明一些实施例中生物礁底水气藏的气井水侵识别相图；
34.图3本发明一些实施例中生物礁底水气藏的钾钠氯系数分布图；
35.图4本发明一些实施例中生物礁底水气藏的钙镁系数分布图；
36.图5本发明一些实施例中生物礁底水气藏的氯镁系数分布图；
37.图6本发明一些实施例中生物礁底水气藏的碳酸盐平衡系数分布图；
38.图7本发明一些实施例中生物礁底水气藏的脱硫系数分布图；
39.图8本发明一些实施例中生物礁底水气藏的碳酸氢根离子浓度和氯离子浓度的分布图；
40.图9本发明一些实施例中生物礁底水气藏的钾钠离子浓度和氯离子浓度的分布图；
41.图10本发明一些实施例中生物礁底水气藏的钙离子浓度和氯离子浓度的分布图；
42.图11本发明一些实施例中生物礁底水气藏的镁离子浓度和氯离子浓度的分布图；
43.图12本发明一些实施例中y1#井水侵阶段初步识别(ⅰ)；
44.图13本发明一些实施例中y1#井水侵阶段初步识别(ⅱ)；
45.图14本发明一些实施例中采用水侵识别相图进一步确认y1#井采气阶段及水侵阶段(ⅰ)；
46.图15本发明一些实施例中采用水侵识别相图进一步确认y1#井采气阶段及水侵阶段(ⅱ)。
具体实施方式
47.下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发 明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范 围。
48.传统的气井水侵阶段划分常常根据产出液中的地层水的含量范围来进行划分，在产出液 中含有地层水到产出液中地层水含量达到水侵标准之间具有较大的区间，因此，导致了气井 水侵阶段的划分具有滞后性，当气井的生产阶段被划分为水侵阶段时，气井中已经出现了明 显甚至严重的水侵，不再具有前瞻性，也不利于气井生产计划调整。
49.基于上述技术问题的考虑，本技术公开了一种底水气藏水侵识别方法，该方法根据气井 产出液中离子浓度变化规律，构建在这一气藏或者同类地层特征的气藏中的标准气井水侵识 别相图，进而判断气井所在的水侵阶段。相较于传统的气井水侵阶段识别，气井水侵阶段被 识别的时间更早，可以识别出气井水侵初期，为工作人员针对气井提前制定生产计划来延长 气井的无水采气期赢得了窗口时间；此外，通过气井产出液中的离子浓度变化规律来识别和 划分气井水侵阶段，丰富了气井水侵识别及气井水侵阶段划分的标准，避免了传统的气井水 侵识别和水侵阶段划分只能通过水气比和产水量划分而导致的划分标准单一，没有可以对比 和检验的手段。
50.本技术公开的一种底水气藏水侵识别方法，参阅图1，包括以下步骤：
51.采集气井产出液在气井水侵各阶段的水化学系数，并根据气井水侵阶段划分水化学系数 的分布范围；
52.根据所述水化学系数的分布范围筛选气井水侵识别系数；
53.采集气井产出液在气井水侵各阶段的不同离子的离子浓度，并根据气井水侵阶段划分不 同离子的离子浓度的分布范围；
54.根据所述离子浓度的分布范围筛选气井水侵识别离子；
55.根据气井水侵识别系数和气井水侵识别离子绘制水侵识别相图；
56.根据所述水侵识别相图判断待测气井的水侵阶段。
57.需要说明的是，所述气井水侵阶段根据所述气井产出液的组成而划分为无水采气期、水 侵初期和带水采气期。所述无水采气期的产出液主要为天然气因温度和压力变化而产出的凝 析液；所述水侵初期的产出液主要为凝析液和地层水的混合液；所述带水采气期
的产出液主 要为地层水。需要说明的是，所述待测气井为该气藏内任一投产气井。
58.需要说明的是，气井水侵识别系数的筛选方法如下：根据各个气井水侵阶段中每一水化 学系数的分布范围的重叠范围来筛选气井水侵识别系数；当一水化学系数的分布范围在各个 气井水侵阶段除端点以外无重叠部分，视为该水化学系数在不同的气井水侵阶段的差异明显， 则该水化学系数可作为气井水侵识别系数；当一水化学系数的分布范围在各个气井水侵阶段 有端点以外的重叠部分，视为该水化学系数在不同的气井水侵阶段的差异不明显，则该水化 学系数不能作为气井水侵识别系数。
59.需要说明的是，气井水侵识别离子的筛选方法如下：根据各个气井水侵阶段中每一离子 浓度的分布范围的重叠范围来筛选气井水侵识别离子；当一离子浓度的分布范围在各个气井 水侵阶段除端点以外无重叠部分，视为该离子的浓度在不同的气井水侵阶段的差异明显，则 该离子可以作为一气井水侵识别系数；当一离子浓度的分布范围在各个气井水侵阶段有端点 以外的重叠部分，视为该离子的浓度在不同的气井水侵阶段的差异不明显，则该离子不能作 为气井水侵识别离子。
60.在一些实施例中，所述水化学系数包括钾钠氯系数、镁钙系数、氯镁系数、碳酸盐平衡 系数和脱硫系数。所述钠钾氯系数为气井采出液中钾离子和钠离子浓度之和与氯离子浓度的 比值；所述镁钙系数为气井采出液中镁离子浓度与该离子浓度的比值；所述氯镁系数为气井 采出液中氯离子浓度与镁离子浓度的比值；所述碳酸盐平衡系数为气井采出液中碳酸根离子 浓度和碳酸氢根离子浓度之和与钙离子浓度的比值；所述脱硫系数为气井采出液中硫酸根离 子浓度与氯离子浓度的比值*100。
61.在一些实施例中，所述离子浓度包括钾钠离子浓度、钙离子浓度、镁离子浓度、氯离子 浓度和碳酸氢根离子浓度。所述钾钠离子浓度为所述气井采出液中钠离子浓度和钾离子浓度 之和。
62.需要说明的是，所述水侵识别相图以气井水侵识别系数和气井水侵识别离子分别作为x 轴，以y轴确定每一x轴的零点；再连接每一所述气井水侵阶段的水侵识别系数和水侵识别 离子的端点值，每一所述气井水侵阶段的水侵识别系数和水侵识别离子的端点值连线围成的 区域即为对应的气井水侵阶段。即所述x轴的数量与所述气井水侵识别系数和所述气井水侵 识别离子的数量一致，每一x轴与y轴的交点都是对应x轴的零点。
63.在一些实施例中，所述气井水侵识别离子为钾钠离子浓度、钙离子浓度和氯离子浓度； 所述气井水侵识别系数为钾钠氯系数、氯镁系数和碳酸盐平衡系数。
64.在一些实施例中，所述水侵识别相图包括6个x轴和一个与x轴垂直的y轴，6个x轴 对称分布在y轴的两侧，每个x轴与y轴的交点即为该x轴的零点。6个x轴分别对应钾 钠离子浓度、钙离子浓度、氯离子浓度、钾钠氯系数、氯镁系数和碳酸盐平衡系数。
65.所述水侵识别相图以每一个水侵识别系数和每一个水侵识别离子的都分别作为x轴，然 后均匀的分布在y轴的两边；选取每一个水侵识别系数和每一个水侵识别离子在每一气井水 侵阶段对应x轴上的端点值(两端端点值：最大值或最小值)；然后用直线连接每一气井水 侵阶段每一水侵识别系数的最大值连接，用直线连接每一气井水侵阶段每一水侵识别系数的 最小值连接，二者圈出的区域即为每一气井水侵阶段的水侵识别相图中的分布。
66.本技术还公开了一种底水气藏产水阶段划分方法，包括以下步骤：
67.采用上述底水气藏水侵识别方法制备待划分气井所在气藏的标准水侵识别相图；
68.采集待划分气井产出液的气井水侵识别系数和气井水侵识别离子，并绘制出气井的实时 水侵识别线图；
69.将气井的实时水侵识别线图垂直投影在所述标准水侵识别相图中，根据所述气井的实时 水侵识别线图落在所述标准水侵识别相图中的区域，判断待划分气井所在的采气阶段。
70.本技术还公开了一种底水气藏产水阶段划分方法，包括以下步骤：采用上述底水气藏水 侵识别方法制备待划分气井所在气藏的标准水侵识别相图；
71.采集待划分气井产出液的气井水侵识别系数和气井水侵识别离子；
72.将采集到的气井水侵识别系数和气井水侵识别离子的数据代入所述标准水侵识别相图中， 判断待划分气井所在的采气阶段。
73.在一些实施例中，上述两种底水气藏产水阶段划分方法可以联用，即先通过将采集到的 待划分气井产出液的气井水侵识别系数和气井水侵识别离子初步得到气井所在的采气阶段， 然后再通过采集到的待划分气井产出液的气井水侵识别系数和气井水侵识别离子绘制的实时 水侵识别线图进一步确定待划分气井所在的采气阶段，对气井所在采气阶段的划分结果做一 次印证。
74.以下结合在某生物礁气藏某一礁带为例对本技术进行说明，该气藏为生物礁高含硫气藏， 于2014年底建成投产，随着气藏开发不断深入，地层压力持续下降，元素硫溶解度逐渐从未 饱和到饱和状态，直至元素硫沉积，适用本技术所公开的方法。其具体实施过程如下：
75.该生物礁气藏为台地边缘生物礁沉积，发育4个礁带，21个礁群、90个单礁体，单个生 物礁规模小，纵向上多期发育、横向同期多个礁体叠置，为高含硫、局部存在底水、受礁滩 体控制的构造—岩性气藏。该礁带长约21.9km，礁带宽度0.8～2km，礁相储层厚30～130m， 均厚65m；气水关系十分复杂，各礁群及礁群内各礁体为独立的气水系统，无统一气水界面。 该礁带的元坝x礁群礁相储层面积约6.76km2，容积法计算地质储量29.24亿方,共发育礁盖 储层相对独立的8个小礁体。
76.该礁带的底水气藏水侵识别方法，具体包括以下步骤：
77.步骤s1：采集该礁带的气井分别在无水采气期、水侵初期和带水采气期的采出液中各离 子的浓度，包括钾离子浓度、钠离子浓度、镁离子浓度、氯离子浓度、碳酸根离子浓度、碳 酸氢根离子浓度、钙离子浓度和硫酸根离子浓度，并根据采集到的各离子浓度，获得气井分 别在无水采气期、水侵初期和带水采气期的钾钠氯系数、镁钙系数、氯镁系数、碳酸盐平衡 系数和脱硫系数的数据，如图3～7。如图3所示，钾钠氯系数在无水采气期、水侵初期和带 水采气期有区别明显的分布范围，且三者之间几乎没有交叉；如图4所示，镁钙系数在无水 采气期、水侵初期和带水采气期在0.1～0.3的范围内均匀分布，在不同的水侵阶段没有明显的 分布范围；如图5、图6所示，氯镁系数和碳酸盐平衡系数在无水采气期、水侵初期和带水 采气期均有区别明显的分布范围，且三者之间几乎没有交叉；如图7所示，脱硫系数在无水 采气期、水侵初期和带水采气期在1～10的范围内有交叉，在不同的水侵阶段没有明显的分布 范围。统计钾钠氯系数、镁钙系数、氯镁系数、碳酸盐平衡系数和脱硫系数在气井在每一个 水侵阶段的具体分布范围，其统计结果如表1所示。
[0078][0079][0080]
表1
[0081]
步骤s2：
[0082]
根据产出液的水化学系数：钾钠氯系数、镁钙系数、氯镁系数、碳酸盐平衡系数、脱硫 系数在气井的无水采气期、水侵初期、带水采气期的数值分布范围是否有除端点以外的重叠 区筛选气井水侵识别系数。从表1可以看出：钾钠氯系数、氯镁系数、碳酸盐平衡系数的数 值分布范围在气井生产的三个采气阶段无端点以外的重叠区、数值区别较明显；同时，镁钙 系数和脱硫系数的数值分布范围在气井生产的三个采气阶段具有明显的重叠区、区别不明显； 据此，筛选出钾钠氯系数、氯镁系数、碳酸盐平衡系数作为气井水侵识别系数。
[0083]
步骤s3：
[0084]
根据采集到的该礁带气井分别在无水采气期、水侵初期和带水采气期的采出液中各离子 的浓度，包括钾离子浓度、钠离子浓度、镁离子浓度、氯离子浓度、碳酸氢根离子浓度和钙 离子浓度的数据，如图8～11所示。在图8中可以看出，在不同的水侵阶段，各气井采出液的 氯离子浓度分布范围有差异明显；从图9中可以看出，在不同的水侵阶段，钾钠离子浓度分 布范围有明显差异；从图10可以看出，在不同的水侵阶段，钙离子浓度分布范围有明显差异； 从图11可以看出，镁离子在不同的水侵阶段，镁离子浓度分布范围没有明显差异。统计其在 气井在每一个水侵阶段的具体分布范围，其统计结果如表2所示。
[0085][0086]
表2
[0087]
需要说明的是，在图8～11中，不同的标识表示由不同气井采出液测得的各离子浓度的数 据。
[0088]
步骤s4：
[0089]
根据产出液中钾钠离子浓度、钙离子浓度、镁离子浓度、氯离子浓度、碳酸氢根离子浓 度在气井的无水采气期、水侵初期、带水采气期数值分布范围是否有除端点以外的重
叠区筛 选气井水侵识别离子。从表2中可以得出：钾钠离子浓度、钙离子浓度、氯离子浓度的数值 分布范围在气井生产的三个采气阶段无端点以外的重叠区、数值区别较明显，同时，镁离子 浓度、碳酸氢根离子浓度的数值分布范围在气井生产的三个采气阶段具有明显的重叠区、区 别不明显。据此，筛选出钾钠离子、钙离子、氯离子作为气井水侵识别离子。
[0090]
步骤s5：
[0091]
根据步骤s2和步骤s4得到的气井水侵识别系数和气井水侵识别离子的分布范围端点值 做气井水侵识别相图，参阅图2。
[0092]
如图2所示，所述气井水侵识别相图包括六个x轴和一个y轴，其中x轴分别为表示氯 离子浓度的x1轴、表示氯镁系数的x2轴、表示碳酸盐平衡系数的x3轴、表示钙离子浓度 的x4轴、表示钾钠氯系数的x5轴和表示钾钠离子浓度的x6轴；通过垂直x轴的y轴将 x1轴、x2轴、x3轴、x4轴、x5轴和x6轴分隔在y轴的两侧，每一个x轴与y轴的交 点即为该x轴的零点，然后在6个x轴上分别对应取每一个气井水侵识别系数在每一个气井 水侵阶段的端点值和每一个气井水侵识别离子在每一个气井水侵阶段的端点值，并将端点值 用直线连接，使气井水侵识别相图被分为四个区域，从内到外依次为：无水采气期、带水采 气期、水侵初期、水侵期。
[0093]
步骤s6
[0094]
采集待识别气井的采出液的气井水侵识别系数和气井水侵识别离子的数据，并将其带入 水侵识别相图中，即可知道待测气井所在的水侵阶段，并判断出是否出现水侵，并通过监测 该气井的气井水侵识别系数和气井水侵识别离子，即可在气井进入水侵初期即可快速判断出 气井的水侵阶段，不存在气井水侵阶段识别的滞后。应用时，还可以采集待识别气井的气井 水侵识别系数和气井水侵识别离子的数据，并根据气井水侵识别系数和气井水侵识别离子的 数据绘出待测气井实时的水侵线图，并将其垂直投影在水侵识别相图中，根据水侵线图落入 水侵识别相图的区域，即可获得该气井所处在的水侵阶段。
[0095]
以该礁带的某一气井y1#为例，通过传统的水侵识别方法：日产水量和水气比来判断该 气井的水侵阶段结果如图12、图13所示。从图12、图13中可以看出：2017年11月30日 开始该井日产水量超过20方/天(图12中时间轴的b点)，出现产地层水情况。
[0096]
采用本技术公开的水侵识别方法，y1#井在2017年4月11日的钾钠氯系数为0.52、氯镁 系数为43.28、碳酸盐平衡系数为0.82，表明y1#井在2017年4月11日开始的水侵识别系数 就已经分布在水侵初期范围内(见图13)，据此，初步判断y1#井的水侵初期时间开始在2017 年4月11日(见图12中时间轴的a点)，采用本发明方法a点时间较b点时间提前6个月。
[0097]
在上述获得的水侵识别相图的基础上，还包括用于对气井进行产水阶段划分的步骤，其 具体包括：
[0098]
步骤s7：
[0099]
采集待划分气井产出液的气井水侵识别系数和气井水侵识别离子；
[0100]
将采集到的气井水侵识别系数和气井水侵识别离子的数据代入所述标准水侵识别相图中， 初步判断待划分气井所在的采气阶段。
[0101]
根据采集到的待划分气井产出液的气井水侵识别系数和气井水侵识别离子绘制出气井的 实时水侵识别线图；
[0102]
将气井的实时水侵识别线图垂直投影在所述标准水侵识别相图中，根据所述气井的实时 水侵识别线图落在所述标准水侵识别相图中的区域，判断待划分气井所在的采气阶段。
[0103]
再以该礁带的某一气井y1#为例，采集该气井在不同日期的气井水侵识别系数和气井水 侵识别离子，将该气井水侵识别系数和气井水侵识别离子的中一个或者多个带入上述获得的 标准水侵识别相图中，初步判断出该气井的水侵阶段；根据采集到的不同日期的该气井水侵 识别系数和气井水侵识别离子的数据做出该气井在该日期的实时水侵识别线图，如图14、图 15所示。将该气井在不同日期获得的实时水侵识别相图垂直投影到图2中，即可获得该气井 在不同日期所处于的水侵阶段。
[0104]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原 则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。