一种边水气藏见水时间预测方法
【技术领域】
[0001] 本发明设及石油和天然气开发技术领域，具体一种边水气藏见水时间的预测方 法。
【背景技术】
[0002] 在水驱气藏的开发中，水侵的早期识别是充分利用无水期进行生产与主动地实施 治水措施的前提。在开发全过程中，必须不断加强地层水监测及水侵动态跟踪分析预测，根 据地层水推进状况做好气井调配，进一步做好地层水防治，有利于提高气藏稳产期和气藏 最终开发效益。
[0003] 目前常用的边水气藏见水预测方法主要有W下两种：（1)气藏工程法，假设地层均 质同性，利用气水两相渗流的基本规律，通过数学推导，得到边水气藏水舌进突破时间的预 测方法，如《边水气藏气井见水时间预测方法》(王会强等.特种油气藏.2008/04)，该方法考 虑了气水流度比、残余气饱和度、储层厚度、孔隙度、距边水距离、气井产量等对气井见水时 间的影响，在相对均质的储层中，实践证明具有一定的实用性。但由于该方法考虑不了储层 的非均质性，对于非均质较严重的储层，见水时间误差很大。（2)数值模拟法，直接建立区块 数值模拟模型来研究边底水的推进是一种普遍手段，但由于其要求参数多，工作量大，特别 是在气藏开发早期的现场应用中存在一定的局限性。
【发明内容】
[0004] 本发明的目的是针对现有技术预测的边水气藏见水时间精度低，常规数值模拟 方法要求参数多、工作量大而难W满足边水气藏气井水侵预警的需要等缺陷，提供了一种 边水气藏见水时间的预测方法。
[0005] 本发明包括W下步骤：
[0006] 1、获取边水气藏地质资料和分析化验数据。
[0007] 收集边水气藏地质资料和分析化验数据。其中，地质资料包括地震、钻井、测井、录 井数据;分析化验数据包括流体高压物性、岩屯、相渗实验数据。
[000引2、依据边水气藏非均质特征，建立不同非均质储层组合形式的Ξ维地质模型，划 定边水气藏的气区和边水区，在边水气藏的气区中建立虚拟井。
[0009] 2.1结合边水气藏地质资料，利用储层预测软件，开展边水气藏储层展布研究，结 合生产测井资料，明确气藏纵向上储层非均质特征，划定边水气藏不同非均质储层的组合 形式。
[0010] 2.2利用步骤1的地质资料和2.1确定的边水气藏不同非均质储层的组合形式，建 立反映边水气藏非均质特征的四种Ξ维地质模型；
[0011] 2.3结合步骤1收集的边水气藏的分析化验数据，对步骤2.2建立的反映边水气藏 非均质特征的Ξ维地质模型进行初始化，划定边水气藏的气区和边水区。
[0012] 2.4依据上述步骤2.2建立的反映边水气藏非均质特征的Ξ维地质模型，在步骤 2.3划定的边水气藏的气区中建立虚拟井。
[0013] 3、依据边水气藏的气区中建立的虚拟井的见水时间与见水响应特征参数之间的 关系，预测边水气藏不同非均质储层的组合形式下的见水时间。
[0014] 3.1优选影响边水气藏气井见水时间的特征参数，分别模拟预测步骤2.4建立的虚 拟井的见水时间，确立不同非均质储层的组合形式下的单一特征参数与虚拟井见水时间 的关系如下：
[0020] 公式（1)、（2)、（3)、（4)、（5)、（6)、（7)分别表示当非均质储层呈现上差下好的组合 形式时，层间渗透率比值、气井产量、层间厚度比值、气井距边水距离、气层总厚度、层间孔 隙度比值、垂向渗透率与水平渗透率比值等7个单一特征参数与虚拟井见水时间的关系式。 公式(8)、（9)、（10)、（11)、（12)、（13)、（14)分别表示当非均质储层呈现上好下差的组合形 式时，上述7个单一特征参数与虚拟井见水时间的关系式。其中，t为虚拟井见水时间，ai~ ai4、bi~bi4为常数，X为单一特征参数。
[0021] 3.2利用步骤3.1得到的全部单一特征参数与虚拟井见水时间的关系进行多元回 归组合，建立非均质储层的组合形式下的多因素特征参数与虚拟井见水时间的关系。
[0022] 3.3依据步骤3.2建立非均质储层的组合形式下的多因素特征参数与虚拟井见水 时间的关系，预测边水气藏不同非均质储层的组合形式下的虚拟井的见水时间。
[0023] 本发明通过建立四种非均质储层的组合形式下Ξ维地质模型，划定边水气藏的气 区和边水区，在边水气藏的气区中建立虚拟井，利用气井见水时间与储层非均质特征参数 的响应关系，可准确预测边水气藏不同非均质储层的组合形式下的气井见水时间，用于指 导气藏开发工作制度的调整，及时下调即将出水井产量和控制生产压差，有效减缓边底水 的推进速度，从而为边水气藏长期稳产提供技术支撑。
【附图说明】
[0024] 图1是本发明技术方案流程框图。
[0025] 图2是普光气藏平均相渗曲线。
[00%] 图3、图4、图5、图6是利用现有技术编制的普光103-1、105-2、103-4、203-1四口单 井测井解释组合图。其中红色（如C所示色块)代表物性最好的一类储层、黄色（如b所示色 块)代表二类储层、绿色(如a所示色块)代表物性最差的Ξ类储层。
[0027] 图7是利用现有技术建立的普光气藏Ξ维地质模型。
[0028] 图8是利用现有技术编制的普光气藏主体边部储层展布连通图。其中红色（如C所 示色块)代表物性最好的一类储层、黄色(如b所示色块)代表二类储层、绿色(如a所示色块） 代表物性最差的Ξ类储层，蓝色(如d所示色块)代表非储层。
[0029] 图9~22是利用现有技术生成的两种非均质储层的组合形式下的单一特征参数与 模拟井见水时间的关系。其中，图9~图15针对下好上差储层，图16~图22针对下差上好储 层。
【具体实施方式】
[0030] 为使本发明目的、技术方案和优点更加清楚，W中原油田普光气藏主体气藏为例， 结合附图对本发明作进一步描述。由图1可知，本发明包含W下步骤。
[0031] 1、收集普光气藏边水气藏地质资料和分析化验数据。
[0032] 地质资料包括150km2S维地震，45 口井的钻井、测井、录井数据;分析化验数据包 括普光2、普光6、普光302-1、普光104-1等4 口井所取气样的体积系数、偏差系数、粘度、密度 等高压物性资料W及120个岩屯、相渗实验数据。
[0033] 对岩屯、相渗实验数据做归一化处理，得到反映普光气藏相渗规律的平均相渗曲 线，见图2。
[0034] 2、依据普光边水气藏非均质特征，建立不同非均质储层组合形式的Ξ维地质模 型，划定边水气藏的气区和边水区，在边水气藏的气区中建立虚拟井。
[0035] 2.1结合测井解释及储层预测手段，确定边水气藏单井储层的非均质性特征，明确 整个边水气藏纵向上储层非均质特征;划定边水气藏不同非均质储层的组合形式。
[0036] 利用步骤1的测井解释数据，对普光气藏边水区气井进行测井解释，绘制如图3、图 4、图5、图6及表1所示的测井解释组合图，确定单井储层纵向非均质特征。由图3、表1可知， 普光北部普光105-1H井纵向上储层段由底部的TSQ1-虹层序往上TSQ1-VI层序储层渗透率 逐渐减小，总体呈现呈上差下好形式，普光104-1、104-3、103-1等3口气井储层特征与普光 105-1H井相似；由图4、表1可知，普光中北部的普光105-2井纵向上中间V层序渗透率较小， 往下Π 层序及往上VI层序渗透率变大，总体呈现两端好中间差，普光103-2储层特征与普光 105-2井类似；由图5、表1可知，普光中南部的普光103-4井纵向上中间V层序渗透率较高， 往下虹层序及往上VI层序渗透率变小，总体呈现两端差中间好的特征；由图6、表1可知，普 光南部的普光203-1井纵向上储层段纵向上储层段由底部的m层序往上VI层序储层渗透率 逐渐变大，总体呈现呈上好下差形式。
[0037] 表1普光气藏8 口单井非均质特征统计表
[00；3 引
[0039] 利用步骤1的地震及测井解释数据，采用波阻抗反演技术进行储层预测，确定由图 7显示的普光气藏边部的储层展布。由图7可知，普光气藏边部储层呈现北部储层上差下好、 中部储层两端好中间差、南部储层上部好下部差的非均质特征，优质储层由北向南逐渐向 上变好的趋势。
[0040] 2.2建立上述步骤2.1所反映的四个普光边水气藏非均质特征的Ξ维地质模型。
[0041] 利用步骤1中的地质资料，采用基于震控的随机建模技术，建立如图7所示的普光 气藏Ξ维地质模型。在此模型上切割出四个反映不同非均质特征的Ξ维模型1区、2区、3区、 4区。依次代表优质储层由北向南逐渐向上变好的趋势，普光103-1、105-2、103-4、203-1井 四口井分别位于1区、2区、3区、4区，切割后单个模型网格数都为mX25X30 = 83250，可W 满足精细模拟研究的需要。
[0042] 2.3依据步骤2.2建立的四个反映普光气藏不同非均质特征的Ξ维地质模型，结合 步骤1的普光2、普光6、普光302-1、普光104-1等4 口气井的流体高压物性资料及图2所示的 反映普光气藏相渗规律的平均相渗曲线，利用静力学平衡方程原理，对四个反映不同非均 质特征的Ξ维地质模型进行初始化，得到如图8所示的四个不同非均质特征的Ξ维地质模 型的气区和边水区。
[0043] 2.4依据上述步骤2.2建立的四种Ξ维地质模型和步骤2.3划定的边水气藏的气区 中建立虚拟井。
[0044] 依据上述步骤2.2建立的反映边水气藏非均质特征的四种Ξ维地质模型，在步骤 2.3划定的边水气藏的气区中建立虚拟井。按照储层上好下差、两端好中间差、中间好两端 差、上差下好四种非均质特征，对4 口虚拟井的渗流特征参数进行赋值