技术摘要：
本发明公开了一种基于断裂预测与动态响应的井间连通性判别方法，包括以下步骤：根据断裂预测方法确定储集体目标井间断裂构造；根据动态响应方法获取储集体目标井间连通信息；结合目标井间断裂构造和连通信息确定目标井间连通性。本发明能够将静态与动态判别方法相结合  全部
背景技术：
目前，油藏井间连通性通常从静态和动态两个方面进行研究。井间静态连通性指 应用地质和物探方法得到连通性结果，其由油藏地质特征和储层特点决定。井间静态连通 性通常采用地层对比、测井等方法来描述地层参数，并依此确定井间连通性。然而，由于缝 洞型油藏储集体具有复杂的特点，因此，传统地质和物探研究方法依然无法有效认识缝洞 体的连通性。国内外常用的井间动态连通性研究方法主要包括示踪剂测试、压力测试、试井 方法如干扰试井、脉冲试井、不稳定试井等，以及基于系统分析思想利用生产动态数据建立 多元线性回归等多种连通性识别方法。 压力资料是判断井间连通性最直接的资料。压力系统的分析是判断井间连通性的 重要依据。根据地层各处原始折算压力近似相等，开采期间，各井地层压力下降趋势、产量 递减趋势类似，以及各井原油密度、组份等一致信息来判断井间连通性是现场运用最多的 方法。 试井方法是研究油藏连通性的有效方法，其包括干扰试井、脉冲试井(如万新德 等，脉冲试井在油田开发中的应用，特种油气藏，2006。聂晶，基于脉冲试井的低渗透油藏井 间连通性分析.长江大学学报(自科版).2013(10)。翟云芳等.利用脉冲试井方法研究低渗 透油藏的连通性.石油学报.2003)和不稳定试井(如林加恩等，应用不稳定压力试井资料定 性分析注采平衡及井间连通性，油气井测试.1997(02)。廖红伟等，应用不稳定试井判断井 间连通性，石油勘探与开发.2002(04))。然而，由于干扰试井和脉冲试井需要改变井的工作 制度，且测试本身周期也长，因此，必然会影响油田的生产计划，同时，高精度的压力计成本 也高。 有学者利用“类干扰试井”方法将大量油井动态数据资料(油压、油嘴变化、产量、 含水率、原油密度等)进行井组分析，以相邻井为基本单元，利用开发过程中的井间干扰信 息，如新井投产、工作制度改变、井采取措施等，在邻井观察能否受到干扰信息，如果存在井 间干扰现象，则说明井间是连通的。如杨敏(塔河油田4区岩溶缝洞型碳酸盐岩储层井间连 通性研究.新疆地质.2004(02))、张小衡，戚明辉等(利用类干扰试井研究塔河油田井间连 通关系.内蒙古石油化工.2009) 多元线性回归分析和电容模型能够定量反映井间相对连通关系。如Albertoni采 用多元线性回归解决了井间连通问题，并取得了较好的效果。该模型回归出定量表征井间 连通程度的权重系数。Yousef等人(A  Capacitance  Model  to  Infer  Innterwell  Connectivity  from  Production  and  Injection  Rate  Fluctuations.SPE  95322.2005) 利用多元线性回归分析，考虑压力数据，以及注入信号的时滞性、衰减性特征，建立基于注 采数据、压力数据的压缩模型。Yousef等人(A  Analysis  and  Interpretation  of  5 CN 111582528 A 说　明　书 2/10 页 Interwell  Connectivity  from  Production  and  Injection  Rate  Fluctuations  Using  Capacitance  Model.SPE  99998.2006)基于电容模型建立渗透率分布和地质特征的诊断工 具，应用连通系数和时间常数比应用单个参数能够更准确地推断地层状况。但是这两种方 法无明确的地质意义。 由于碳酸盐岩缝洞型油藏储集体类型多样，缝洞组合关系复杂，连续性差，非均质 性极强，很难采用单一方法来准确判断和预测井间裂缝分布和井间连通性，因此，有必要建 立静态与动态结合的综合性判断方法。
技术实现要素：
针对上述问题，本发明的目的是提供一种基于断裂预测与动态响应的井间连通性 判别方法，其能够将静态与动态判别方法相结合，简捷高效地识别目标井间裂缝和连通性， 准确性好。 为实现上述目的，本发明采取以下技术方案：一种基于断裂预测与动态响应的井 间连通性判别方法，包括以下步骤：根据断裂预测方法确定储集体目标井间断裂构造；根据 动态响应方法获取储集体目标井间连通信息；结合目标井间断裂构造和连通信息确定目标 井间连通性。 在一个具体实施例中，断裂构造包括主断裂和次级断裂。 在一个具体实施例中，结合相干属性和振幅切片确定主断裂。 在一个具体实施例中，结合蚂蚁体和体曲率属性确定次级断裂。 在一个具体实施例中，利用注采响应法和示踪剂法获取注采井间连通信息。 在一个具体实施例中，利用类干扰法获取采油井间连通信息。 在一个具体实施例中，注采响应法通过分析生产动态数据含水率、产液量或产油 量波动特征获取注采井间连通信息；示踪剂法是指向注入井注入示踪剂段塞，在周围生产 井监测示踪剂的产出情况，确定示踪剂产出浓度随时间变化的曲线，对该曲线进行分析或 历史拟合，获取注采井间连通信息。 在一个具体实施例中，类干扰法包括采用新井类干扰、生产特征相似性、见水时间 一致和压力趋势下降方法相结合获取采油井间连通信息。 在一个具体实施例中，新井类干扰、生产特征相似性或见水时间一致方法获取的 井间连通信息分别赋权值为1，压力趋势下降方法获取的井间连通信息赋权值为2；若采油 井之间的累计连通权值大于2，则确定采油井之间连通。 在一个具体实施例中，生产特征相似性方法包括利用灰色关联度分析或动态时间 规整方法确定生产动态曲线相似度。 在一个具体实施例中，利用灰色关联度分析方法确定采油井的连通信息包括以下 步骤：(1)获取采油井的产液量数据；(2)计算采油井之间的灰色关联系数；(3)确定采油井 之间的灰色关联度；(4)根据采油井之间灰色关联度确定采油井之间是否连通。 在一个具体实施例中，确定采油井之间的灰色关联度包括以下步骤： ①选择多口采油井，确定各采油井生产数据对应序列所采用的公式为： Xi＝{xi(1)，xi(2)，......xi(n)}    (1) Xj＝{xj(1)，xj(2)，......xj(n)}    (2) 6 CN 111582528 A 说　明　书 3/10 页 式中，Xi为第i口采油井的生产数据，xi(n)为第i口采油井第n处位置的生产数据， Xj为第j口采油井的生产数据，xj(n)为第j口采油井第n处位置的生产数据，其中，每一口采 油井各位置处的生产数据为无量纲值； ②计算多口采油井对应位置生产数据的差值所采用的公式为： Δij(k)＝|xi(k)-xj(k)|，(k＝1，2，...n)    (3) Xij＝{Δij(1)，Δij(2)，...Δij(n)}    (4) 式中，Δij(k)为第i、j口采油井第k处位置生产数据对应的差值，xi(k)为第i口采 油井第k处位置的生产数据，xj(k)为第i口采油井第k处位置的生产数据，Xij为第i、j口采油 井各对应位置生产数据差值的集合； ③确定各采油井极大和极小位置处的生产数据最小绝对差所采用的公式为： 式中，Δmin为最小绝对差，x0(k)为参考井的生产序列，灰色关联进行匹配的时候 选取生产数据波动最大的井作为参考井，该参考井的连通可能性最大，其它的井和该参考 井进行相似度计算，内层的 为当前待对比的序列与参考序列间的对应位置 距离的最小值，外层的 为所有待对比的序列与参考序列之间对应位置 距离的全局最小值。 ④确定各采油井极大和极小位置处的生产数据最大绝对差所采用的公式为： 式中，Δmax为最大绝对差，内层的 为当前待对比的序列与参考序 列间的对应位置距离的最大值，外层的 为所有待对比的序 列与参考序列之间对应位置距离的全局最大值。 ⑤计算多口采油井的各项生产数据关联度系数所采用的公式为： 式中，ρ＝0.5，ρ越小表示精确度越高，ξij(k)为多口采油井的生产数据各项关联度 系数； ⑥计算多口采油井的各项生产数据灰色关联度系数所采用的公式为： 式中，rij为多口采油井的各项生产数据灰色关联度系数。 在一个具体实施例中，动态时间规整方法确定生产动态曲线相似度包括在确定灰 色关联系数的基础上采用动态时间规整算法对生产动态曲线相似度进行量化判别。 7 CN 111582528 A 说　明　书 4/10 页 在一个具体实施例中，动态时间规整方法的计算流程如下：首先，选择多口采油 井，分别确定其序列分别为： Q＝q1，q2，...，qi，...，qn    (9) C＝c1，c2，...，ci，...，cm    (10) 式中，序列Q中q1......qn是采油井Q的生产数据序列；序列C中c1 ......cm是采油 井C的生产数据序列。 接着，执行步骤如下： ①判断n是否等于m，若等于，则直接比较对应位置，否则进入下一步； ②构造一个n×m的矩阵，矩阵元素d(qi，cj)表示qi和cj之间的距离，确定序列Q的 每一点和C的每一点之间的相似度，距离越小则相似度越高； ③依据状态转移方程计算n×m矩阵总的累计距离，并记录累计次数t； ④计算动态时间规整值所采用的公式为： Dis(i，j)＝d(qi，cj) min{Dis(i-1，j-1)，Dis(i，j-1)，Dis(i-1，j)}    (11) 式中，Dis(n，m)为矩阵总的累计距离，t为累积次数，Dis(i，j)、Dis(i-1，j-1)、Dis (i，j-1)、Dis(i-1，j)分别为i×j矩阵总的累计距离，(i-1)×(j-1)矩阵总的累计距离，i× (j-1)矩阵总的累计距离，(i-1)×j矩阵总的累计距离。 在一个具体实施例中，确定矩阵元素之间的距离所采用的公式为： d(qi，cj)＝(qi-cj)    (13) 式中，d(qi，cj)为qi和cj之间的距离。 本发明由于采取以上技术方案，其具有以下优点：本发明方法结合相干属性和振 幅切片能够有效确定目标井间主断裂，利用不受反射振幅和波形影响的体曲率属性来改进 蚂蚁体追踪，能够更好地确定目标井间次级断裂，利用动态响应方法从动态上获取目标井 间连通信息，最后将动态响应结果与断裂预测相结合，确定出目标井间连通状况，本发明方 法能够克服单一方法评价井间连通性的局限性，并且能够提高复杂储集体井间连通性判断 的准确度，为注水参数优化提供基础。 本发明的其它特征和优点将在随后的说明书中阐述，并且，部分的从说明书中变 得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利 要求书以及附图中所特别指出的结构来实现和获得。 附图说明 为了更清楚的说明本发明实施例中的技术方案，下面将对实施例描述中所需要的 附图做简单的介绍： 图1是本发明基于断裂预测与动态响应的井间连通性判别方法的一个具体实施例 的流程结构示意图； 图2是本发明动态时间规整法的结构示意图； 图3是本发明 面主断裂分布的结构示意图； 图4是本发明 面次级断裂分布的结构示意图； 8 CN 111582528 A 说　明　书 5/10 页 图5是本发明示踪剂突破时间的示意图； 图6是本发明示踪剂峰值浓度的示意图； 图7是本发明注采响应图的示意图。