CN111444462A - 根据不稳定试井测算串珠体数据的方法及设备 - Google Patents
根据不稳定试井测算串珠体数据的方法及设备 Download PDFInfo
- Publication number
- CN111444462A CN111444462A CN202010258793.2A CN202010258793A CN111444462A CN 111444462 A CN111444462 A CN 111444462A CN 202010258793 A CN202010258793 A CN 202010258793A CN 111444462 A CN111444462 A CN 111444462A
- Authority
- CN
- China
- Prior art keywords
- bead body
- well
- data
- bottom hole
- bead
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011324 bead Substances 0.000 title claims abstract description 174
- 238000012360 testing method Methods 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 79
- 238000005553 drilling Methods 0.000 claims abstract description 47
- 238000011084 recovery Methods 0.000 claims abstract description 36
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 230000000704 physical effect Effects 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 22
- 230000002093 peripheral effect Effects 0.000 claims description 19
- 239000012530 fluid Substances 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 8
- 230000035699 permeability Effects 0.000 claims description 7
- 238000007906 compression Methods 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 9
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000001788 irregular Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 241001146209 Curio rowleyanus Species 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- Mathematical Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Data Mining & Analysis (AREA)
- Geophysics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明实施例提供了一种根据不稳定试井测算串珠体数据的方法及设备。所述方法包括：获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。本发明实施例提供的根据不稳定试井测算串珠体数据的方法及设备，可以基于不稳定测试井数据快速测算串珠体物性及体积。
Description
技术领域
本发明实施例涉及碳酸盐岩开采技术领域，尤其涉及一种根据不稳定试井测算串珠体数据的方法及设备。
背景技术
我国塔里木盆地海相碳酸盐岩储层非均质极强，存在的孔-缝-洞系统类型丰富多样：溶蚀孔隙发育的孔隙型碳酸盐岩储层、天然构造裂缝和微裂缝广泛存在的裂缝型碳酸盐岩储层，空间尺寸较大(填充)溶洞体存在的缝洞型碳酸盐岩储层。由于(填充)溶洞体内流体与外围储层的波阻抗值不同，地震反射波在洞体与储层之间会形成明显的波阻抗界面。因此，该类储层在地震响应特征和地震属性上具有明显的串珠状反射现象(简称为“串珠体”)。储层串珠体较外围储层物性极好，是油气储集与流动的主要空间；且串珠体空间大小是评价储层物性、钻遇串珠生产井产能的重要指标。目前，对于串珠体大小的刻画还未有较为成熟的技术。因此，开发一种根据不稳定试井测算串珠体数据的方法，可以有效填补上述相关技术缺陷，就成为业界广泛关注的技术问题。
发明内容
针对现有技术存在的上述问题，本发明实施例提供了一种根据不稳定试井测算串珠体数据的方法及设备。
第一方面，本发明的实施例提供了一种根据不稳定试井测算串珠体数据的方法，包括：获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述变化数据获取井底压力随关井时长的测试点关系数据，包括：将关井后每个测试点的时刻数值及井底压力数值分别减去关井时刻对应的时刻数值和井底压力数值，得到各测试点的井底恢复压力与关井时长差的关系数据。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线，相应地，所述井底恢复压力导数包括：
其中，j为第j个测试点；Δpj为第j个测试点的井底恢复压力；Δtj为第j个测试点的关井时长差。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述双对数曲线获取若干特征值，包括：获取所述双对数曲线的第一水平段拐点对应的横坐标和纵坐标数值为第一特征值，第二水平段拐点对应的横坐标和纵坐标数值为第二特征值，第三水平段拐点对应的横坐标和纵坐标数值作为第三特征值。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体物性，包括：
其中，χ为串珠体对外围储层的流动系数比；DL1-2为所述第一特征值和第二特征值在双对数曲线上的纵向距离；DL2-3为所述第二特征值和第三特征值在双对数曲线上的纵向距离；UL为双对数曲线上纵坐标主刻度的步长值；ω为串珠体对外围储层的储容比；a和b为系数；k为串珠体外围储层渗透率；φ为串珠体外围储层孔隙度；ct为储层综合压缩系数；μ为流体粘度；η为串珠体储层波导系数；η′为串珠体外围储层波导系数。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体体积，包括：
其中，Rmin为串珠体垂直轴心距离井筒中轴线的距离；rw为井筒半径；Δt1为所述第一特征值的关井时刻值；Requ为串珠体等效渗流体垂直边界距离井筒中轴线的距离；Δt2为所述第二特征值的关井时刻值；V为串珠体体积。
在上述方法实施例内容的基础上，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体物性，还包括：
其中，β为串珠体的垂向尺寸极大值；α为串珠体的径向尺寸极大值。
第二方面，本发明的实施例提供了一种根据不稳定试井测算串珠体数据的装置，包括：
双对数曲线获取模块，用于获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；
分析图版获取模块，用于根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；
串珠体数据测算模块，用于采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
第三方面，本发明的实施例提供了一种电子设备，包括：
至少一个处理器；以及
与处理器通信连接的至少一个存储器，其中：
存储器存储有可被处理器执行的程序指令，处理器调用程序指令能够执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的根据不稳定试井测算串珠体数据的方法。
第四方面，本发明的实施例提供了一种非暂态计算机可读存储介质，非暂态计算机可读存储介质存储计算机指令，计算机指令使计算机执行第一方面的各种可能的实现方式中任一种可能的实现方式所提供的根据不稳定试井测算串珠体数据的方法。
本发明实施例提供的根据不稳定试井测算串珠体数据的方法及设备，通过采用若干测试点关系数据绘制出井底恢复压力及其导数的双对数曲线，在双对数曲线上选取特征值并进一步构建串珠体储层井底压力分析图版，可以基于不稳定测试井数据快速测算串珠体物性及体积。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单的介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的根据不稳定试井测算串珠体数据的方法流程图；
图2为本发明实施例提供的不稳定试井对应的井底压力测试数据示意图；
图3为本发明实施例提供的双对数坐标下的关井井底恢复压力与井底恢复压力导数示意图；
图4为本发明实施例提供的串珠体等效渗透流体积示意图；
图5为本发明实施例提供的单串珠空间信息示意图；
图6为本发明实施例提供的根据不稳定试井测算串珠体数据的装置结构示意图；
图7为本发明实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。另外，本发明提供的各个实施例或单个实施例中的技术特征可以相互任意结合，以形成可行的技术方案，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时，应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
为了快速测算得到缝洞型碳酸盐岩储层中钻遇串珠体的物性参数和空间信息，需要一种利用钻遇井的不稳定试井双对数曲线来直接确定串珠体物性参数及空间信息的方法，其目的在于：通过对钻遇串珠体的垂直井进行不稳定试井测试，依据测试数据在双对数图版的位置信息以实现对钻遇串珠物性参数及空间体积的确定，为生产井单井产能的评价、区块井网加密调整方案的制定提供理论依据。具体包括以下内容：考虑不同串珠体的孔隙度、渗透率等物性参数以及串珠体积、位置、个数等空间信息对不稳定试井数据的影响，提出“串珠等效渗流体积”原理；基于“串珠等效渗流体积”原理，根据不同钻遇井的测试数据在时间-压力双对数坐标下的位置信息，可以表征串珠物性参数和空间信息的特征点以及特征值，并定性地比较不同井钻遇的串珠物性好坏程度、空间体积大小；利用建立的串珠储层井底压力分析图版以及特征点数值，可以定量地确定串珠孔隙度、渗透率等物性参数以及串珠体积、位置、个数等空间信息。基于这种思想，本发明实施例提供了一种根据不稳定试井测算串珠体数据的方法，参见图1，该方法包括：
101、获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；
102、根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；
103、采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述变化数据获取井底压力随关井时长的测试点关系数据，包括：将关井后每个测试点的时刻数值及井底压力数值分别减去关井时刻对应的时刻数值和井底压力数值，得到各测试点的井底恢复压力与关井时长差的关系数据。具体地可以参见图2，通过设置井底精密压力计，测试关井后目标井的井底压力变化数据(时间数值tj、压力数值pj)；将关井对应时刻后的每个测试的时间数值tj、压力数值pj分别减去关井时刻对应的时间数值t0、压力数值p0，得到各个测试点关井时长差△tj和恢复压力△pj数据。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线，相应地，所述井底恢复压力导数包括：
其中，j为第j个测试点；Δpj为第j个测试点的井底恢复压力；Δtj为第j个测试点的关井时长差。具体可以参见图3，采用(1)式得到各个测试点的关井时长差△tj和井底恢复压力导数△p’·△t数据；在双对数坐标系中绘制关井时长差△t和井底恢复压力△p的曲线△p(△t)，以及关井时长差△t和井底恢复压力导数△p’·△t的曲线△p’·△t(△t)。在图3中，△p’1为第一特征值；△p’2为第二特征值；△p’3为第三特征值；△t1为关井时间的第一拐点值；△t2为关井时间的第二拐点值；DL1-2为△p’1与△p’2在双对数曲线上的纵向距离；DL2-3为△p’2与△p’3在双对数曲线上的纵向距离；UL为双对数曲线上纵坐标主刻度的步长值。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述根据所述双对数曲线获取若干特征值，包括：获取所述双对数曲线的第一水平段拐点对应的横坐标和纵坐标数值为第一特征值，第二水平段拐点对应的横坐标和纵坐标数值为第二特征值，第三水平段拐点对应的横坐标和纵坐标数值作为第三特征值。具体可以参见图3，从双对数坐标中取恢复压力导数曲线上第一水平段拐点对应的横坐标和纵坐标数值(△t1,△p’1)为第一特征值，第二水平段拐点对应的横坐标和纵坐标数值(△t2,△p’2)为第二特征值，第三水平段对应的纵坐标数值△p’3为第三特征值。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体物性，包括：
其中，χ为串珠体对外围储层的流动系数比；DL1-2为所述第一特征值和第二特征值在双对数曲线上的纵向距离；DL2-3为所述第二特征值和第三特征值在双对数曲线上的纵向距离；UL为双对数曲线上纵坐标主刻度的步长值；ω为串珠体对外围储层的储容比；a和b为系数；k为串珠体外围储层渗透率；φ为串珠体外围储层孔隙度；ct为储层综合压缩系数；μ为流体粘度；η为串珠体储层波导系数；η′为串珠体外围储层波导系数。系数a、b由实际目标区岩石取芯的串珠体外围储层孔隙度φ和串珠体外围储层渗透率k交汇图或者目标井所在油田矿场经验公式得到。具体可以参见图4，串珠等效渗流体体积原理为任何一个、多个串珠体或具有非规则外边界形状的串珠体等复杂的储层区域可以用一个规则的储层区域替换，两者虽然在不稳定试井过程中具有不相同的井底压力变化曲线，但两者对应的井底压力变化曲线的特征值是相同的。基于此现象，直接利用曲线特征值以及一个规则的储层区域体积来计算复杂不规则的串珠体积。图4中Rmin为串珠体垂直轴心距离井筒中轴线的距离；Requ为串珠体等效渗流体垂直边界距离井筒中轴线的距离。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体体积，包括：
其中，Rmin为串珠体垂直轴心距离井筒中轴线的距离；rw为井筒半径；Δt1为所述第一特征值的关井时刻值；Requ为串珠体等效渗流体垂直边界距离井筒中轴线的距离；Δt2为所述第二特征值的关井时刻值；V为串珠体体积。
基于上述方法实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的方法，所述通过所述特征值比较不同钻遇井对应的串珠体物性，还包括：
其中，β为串珠体的垂向尺寸极大值；α为串珠体的径向尺寸极大值。
具体可以参见图5，α为串珠体的径向尺寸极大值；β为串珠体的垂向尺寸极大值；Rj为串珠体离散后的第j层串珠外边界距离井筒中轴线的距离；Rmin为串珠体垂直轴心距离井筒中轴线的距离。
本发明实施例提供的根据不稳定试井测算串珠体数据的方法，通过采用若干测试点关系数据绘制出井底恢复压力及其导数的双对数曲线，在双对数曲线上选取特征值并进一步构建串珠体储层井底压力分析图版，可以基于不稳定测试井数据快速测算串珠体物性及体积。
由于考虑了串珠体物性参数以及串珠空间信息，所以进行串珠钻遇井的不稳性试井解释提高了压力监测解释的精准度，且更符合实际情况。另外，提出了“串珠等效渗流体积”原理能根据不稳定试井的双对数曲线，直接、快速地测算钻遇串珠的物性参数及空间信息。本发明各个实施例提出的方法为钻遇井的单井产能、储量的评价以及后续井网密度、生产制度的调整具有指导作用。具体体现在：提出了一种通过钻遇井不稳定试井数据测算串珠物性参数、空间信息的方法和系统，与现存的通过数值模拟方法刻画钻遇串珠体的方法相比，本发明提出方法可以避免大量的微分方程求解过程以及理论模型与实际数据繁琐拟合的过程；考虑了串珠体孔隙度、渗透率等物性参数以及串珠大小、位置、数量等空间信息，克服了不得已假设“串珠体为一个圆柱体内区储层”而强行使用复合储层模型解释的不合理性，提高了压力监测解释的精准度，且更符合实际情况；本专利提出的“串珠等效渗流体积”原理简单易懂、特征参数均具有明确物理意义，丰富了径向复合类储层不稳定渗流理论；本专利提出的“串珠等效渗流体积”原理，仅通过不稳定试井数据在双对数曲线特征值的位置信息便可以快速定量地判断串珠体的物性参数、定性的判断串珠体的空间信息，对测试数据的质量要求低，可操作性强。
本发明各个实施例的技术方案属于油气田开发领域，涉及确定缝洞型碳酸盐岩生产井钻遇串珠体的物性参数和空间信息的油藏工程方法，具体涉及通过利用井底测试压力分析储层缝洞体信息的不稳定试井方法。本发明尤其针对具有明显串珠状反射(强串珠反射、串珠反射、弱串珠反射)的缝洞型碳酸盐岩垂直钻遇井，对于钻遇圆柱状(填充)溶洞、缝洞单元体的垂直井仍具有使用价值，但由于与现存复合储层的试井解释方法得到结果相同而无明显的优势。
本发明各个实施例的实现基础是通过具有处理器功能的设备进行程序化的处理实现的。因此在工程实际中，可以将本发明各个实施例的技术方案及其功能封装成各种模块。基于这种现实情况，在上述各实施例的基础上，本发明的实施例提供了一种根据不稳定试井测算串珠体数据的装置，该装置用于执行上述方法实施例中的根据不稳定试井测算串珠体数据的方法。参见图6，该装置包括：
双对数曲线获取模块601，用于获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；
分析图版获取模块602，用于根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；
串珠体数据测算模块603，用于采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
本发明实施例提供的根据不稳定试井测算串珠体数据的装置，采用双对数曲线获取模块、分析图版获取模块和串珠体数据测算模块，通过采用若干测试点关系数据绘制出井底恢复压力及其导数的双对数曲线，在双对数曲线上选取特征值并进一步构建串珠体储层井底压力分析图版，可以基于不稳定测试井数据快速测算串珠体物性及体积。
需要说明的是，本发明提供的装置实施例中的装置，除了可以用于实现上述方法实施例中的方法外，还可以用于实现本发明提供的其他方法实施例中的方法，区别仅仅在于设置相应的功能模块，其原理与本发明提供的上述装置实施例的原理基本相同，只要本领域技术人员在上述装置实施例的基础上，参考其他方法实施例中的具体技术方案，通过组合技术特征获得相应的技术手段，以及由这些技术手段构成的技术方案，在保证技术方案具备实用性的前提下，就可以对上述装置实施例中的装置进行改进，从而得到相应的装置类实施例，用于实现其他方法类实施例中的方法。例如：
基于上述装置实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的装置，还包括：减差模块，用于将关井后每个测试点的时刻数值及井底压力数值分别减去关井时刻对应的时刻数值和井底压力数值，得到各测试点的井底恢复压力与关井时长差的关系数据。
基于上述装置实施例的内容，作为一种可选的实施例，本发明实施例中提供的根据不稳定试井测算串珠体数据的装置，还包括：特征值模块，用于获取所述双对数曲线的第一水平段拐点对应的横坐标和纵坐标数值为第一特征值，第二水平段拐点对应的横坐标和纵坐标数值为第二特征值，第三水平段拐点对应的横坐标和纵坐标数值作为第三特征值。
本发明实施例的方法是依托电子设备实现的，因此对相关的电子设备有必要做一下介绍。基于此目的，本发明的实施例提供了一种电子设备，如图7所示，该电子设备包括：至少一个处理器(processor)701、通信接口(Communications Interface)704、至少一个存储器(memory)702和通信总线703，其中，至少一个处理器701，通信接口704，至少一个存储器702通过通信总线703完成相互间的通信。至少一个处理器701可以调用至少一个存储器702中的逻辑指令，以执行如下方法：获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
此外，上述的至少一个存储器702中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。例如包括：获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。基于这种认识，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
在本专利中，术语"包括"、"包含"或者其任何其它变体意在涵盖非排它性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其它要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句"包括……"限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种根据不稳定试井测算串珠体数据的方法，其特征在于，包括：
获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；
根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；
采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
2.根据权利要求1所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述根据所述变化数据获取井底压力随关井时长的测试点关系数据，包括：将关井后每个测试点的时刻数值及井底压力数值分别减去关井时刻对应的时刻数值和井底压力数值，得到各测试点的井底恢复压力与关井时长差的关系数据。
3.根据权利要求2所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线，相应地，所述井底恢复压力导数包括：
其中，j为第j个测试点；Δpj为第j个测试点的井底恢复压力；Δtj为第j个测试点的关井时长差。
4.根据权利要求3所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述根据所述双对数曲线获取若干特征值，包括：获取所述双对数曲线的第一水平段拐点对应的横坐标和纵坐标数值为第一特征值，第二水平段拐点对应的横坐标和纵坐标数值为第二特征值，第三水平段拐点对应的横坐标和纵坐标数值作为第三特征值。
5.根据权利要求4所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述通过所述特征值比较不同钻遇井对应的串珠体物性，包括：
其中，χ为串珠体对外围储层的流动系数比；DL1-2为所述第一特征值和第二特征值在双对数曲线上的纵向距离；DL2-3为所述第二特征值和第三特征值在双对数曲线上的纵向距离；UL为双对数曲线上纵坐标主刻度的步长值；ω为串珠体对外围储层的储容比；a和b为系数；k为串珠体外围储层渗透率；φ为串珠体外围储层孔隙度；ct为储层综合压缩系数；μ为流体粘度；η为串珠体储层波导系数；η′为串珠体外围储层波导系数。
6.根据权利要求5所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述通过所述特征值比较不同钻遇井对应的串珠体体积，包括：
其中，Rmin为串珠体垂直轴心距离井筒中轴线的距离；rw为井筒半径；Δt1为所述第一特征值的关井时刻值；Requ为串珠体等效渗流体垂直边界距离井筒中轴线的距离；Δt2为所述第二特征值的关井时刻值；V为串珠体体积。
7.根据权利要求6所述的根据不稳定试井测算串珠体数据的方法，其特征在于，所述通过所述特征值比较不同钻遇井对应的串珠体物性，还包括：
其中，β为串珠体的垂向尺寸极大值；α为串珠体的径向尺寸极大值。
8.一种根据不稳定试井测算串珠体数据的装置，其特征在于，包括：
双对数曲线获取模块，用于获取钻遇井在关井阶段的井底压力随关井时长的变化数据，根据所述变化数据获取井底压力随关井时长的测试点关系数据，根据所述测试点关系数据获取井底恢复压力及井底恢复压力导数相对所述关井时长的双对数曲线；
分析图版获取模块，用于根据所述双对数曲线获取若干特征值，通过所述特征值比较不同钻遇井对应的串珠体物性和体积，并构建串珠体储层井底压力分析图版；
串珠体数据测算模块，用于采用所述串珠体储层井底压力分析图版，对钻井串珠体数据进行测算。
9.一种电子设备，其特征在于，包括：
至少一个处理器、至少一个存储器和通信接口；其中，
所述处理器、存储器和通信接口相互间进行通信；
所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令，以执行权利要求1至7任一项权利要求所述的方法。
10.一种非暂态计算机可读存储介质，其特征在于，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行权利要求1至7中任一项权利要求所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010258793.2A CN111444462B (zh) | 2020-04-03 | 2020-04-03 | 根据不稳定试井测算串珠体数据的方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010258793.2A CN111444462B (zh) | 2020-04-03 | 2020-04-03 | 根据不稳定试井测算串珠体数据的方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111444462A true CN111444462A (zh) | 2020-07-24 |
| CN111444462B CN111444462B (zh) | 2022-02-08 |
Family
ID=71652785
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010258793.2A Expired - Fee Related CN111444462B (zh) | 2020-04-03 | 2020-04-03 | 根据不稳定试井测算串珠体数据的方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111444462B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113673101A (zh) * | 2021-08-18 | 2021-11-19 | 中国石油大学(北京) | 串珠状缝洞体溶洞体积的计算方法和处理器 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030102163A1 (en) * | 2001-11-29 | 2003-06-05 | Hemphill Alan Terry | Method for determining sweep efficiency for removing cuttings from a borehole |
| US20170053046A1 (en) * | 2015-08-17 | 2017-02-23 | Saudi Arabian Oil Company | Capillary pressure analysis for petrophysical statistical modeling |
| CN107237626A (zh) * | 2016-03-29 | 2017-10-10 | 中国石油化工股份有限公司 | 一种井钻遇溶洞体积的确定方法 |
| CN107577831A (zh) * | 2017-03-17 | 2018-01-12 | 西南石油大学 | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110096718A (zh) * | 2018-01-30 | 2019-08-06 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中溶洞的体积的方法 |
| CN110107277A (zh) * | 2018-01-30 | 2019-08-09 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 |
-
2020
- 2020-04-03 CN CN202010258793.2A patent/CN111444462B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030102163A1 (en) * | 2001-11-29 | 2003-06-05 | Hemphill Alan Terry | Method for determining sweep efficiency for removing cuttings from a borehole |
| US20170053046A1 (en) * | 2015-08-17 | 2017-02-23 | Saudi Arabian Oil Company | Capillary pressure analysis for petrophysical statistical modeling |
| CN107237626A (zh) * | 2016-03-29 | 2017-10-10 | 中国石油化工股份有限公司 | 一种井钻遇溶洞体积的确定方法 |
| CN107577831A (zh) * | 2017-03-17 | 2018-01-12 | 西南石油大学 | 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110096718A (zh) * | 2018-01-30 | 2019-08-06 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中溶洞的体积的方法 |
| CN110107277A (zh) * | 2018-01-30 | 2019-08-09 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 |
Non-Patent Citations (4)
| Title |
|---|
| WENYANG SHI等: "Pressure transient behavior of layered commingled reservoir with vertical inhomogeneous closed boundary", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
| YANG XIAO等: "Dynamic and Static Combination Method for Fracture-Vug Unit Division of Fractured-Vuggy Reservoirs", 《ARABIAN JOURNAL FOR SCIENCE AND ENGINEERING》 * |
| 刘洪等: "缝洞型油藏钻遇溶洞油井的压力曲线特征", 《岩性油气藏》 * |
| 熊万里: "缝洞型碳酸盐岩试井实验研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113673101A (zh) * | 2021-08-18 | 2021-11-19 | 中国石油大学(北京) | 串珠状缝洞体溶洞体积的计算方法和处理器 |
| CN113673101B (zh) * | 2021-08-18 | 2024-03-29 | 中国石油大学(北京) | 串珠状缝洞体溶洞体积的计算方法和处理器 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111444462B (zh) | 2022-02-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11782741B2 (en) | Modeling of fluid introduction and/or fluid extraction elements in simulation of coreflood experiment | |
| CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN104533400B (zh) | 一种重构测井曲线的方法 | |
| EP3108099B1 (en) | Measuring behind casing hydraulic conductivity between reservoir layers | |
| CN110656915B (zh) | 一种页岩气多段压裂水平井多工作制度产能预测方法 | |
| CN109033504B (zh) | 一种油水井套管损坏预测方法 | |
| WO2017120099A1 (en) | Modeling to characterize fractures network in homogeneous petroleum reservoirs | |
| CN108572401B (zh) | 缝洞组合模型的构建方法及探测储层缝洞变形的方法 | |
| CN107038290A (zh) | 一种考虑增压效应力的井壁坍塌压力的计算方法 | |
| EP3004542B1 (en) | A system for computing the radius of investigation in a radial, composite reservoir | |
| CN107237626B (zh) | 一种井钻遇溶洞体积的确定方法 | |
| CN111444462B (zh) | 根据不稳定试井测算串珠体数据的方法及设备 | |
| CN109798108B (zh) | 深层低渗-致密储层复杂流体性质的识别方法和装置 | |
| CN111963149B (zh) | 一种考虑滞地液量增压的压裂后地层压力求取方法 | |
| CN114922614A (zh) | 一种控压钻井工况下的地层压力监测方法 | |
| Johnston et al. | Predicting success in the Haynesville Shale: a geologic, completion, and production analysis | |
| CN114565116A (zh) | 一种低渗储层物性参数的预测方法和装置 | |
| US10753203B2 (en) | Systems and methods to identify and inhibit spider web borehole failure in hydrocarbon wells | |
| CN107291968B (zh) | 压裂选段方法与系统 | |
| US20230259670A1 (en) | Permeability modeling in a reservoir simulation model using dynamic pressure transient analysis | |
| CN113792932B (zh) | 一种利用微震-损伤-渗流关系的页岩气产量预测方法 | |
| CN117307142A (zh) | 一种确定填充溶洞储层渗透率敏感系数的方法及系统 | |
| US20220229950A1 (en) | Hybrid procedure for evaluating stress magnitude and distribution on a liner | |
| CN115271367A (zh) | 高原隧道围岩分级方法、装置、设备以及存储介质 | |
| CN115344821A (zh) | 断溶体油藏静态储量门槛值确定方法、电子设备及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220208 |