CN114086933B - 一种试井解释方法、装置、存储介质以及计算机设备 - Google Patents
一种试井解释方法、装置、存储介质以及计算机设备 Download PDFInfo
- Publication number
- CN114086933B CN114086933B CN202010767808.8A CN202010767808A CN114086933B CN 114086933 B CN114086933 B CN 114086933B CN 202010767808 A CN202010767808 A CN 202010767808A CN 114086933 B CN114086933 B CN 114086933B
- Authority
- CN
- China
- Prior art keywords
- fracture
- karst
- stratum
- model
- well
- 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.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 259
- 238000000034 method Methods 0.000 title claims abstract description 90
- 239000002131 composite material Substances 0.000 claims abstract description 81
- 239000011148 porous material Substances 0.000 claims description 94
- 239000012530 fluid Substances 0.000 claims description 73
- 230000015572 biosynthetic process Effects 0.000 claims description 28
- 230000035699 permeability Effects 0.000 claims description 26
- 238000011084 recovery Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 17
- 230000006835 compression Effects 0.000 claims description 16
- 238000007906 compression Methods 0.000 claims description 16
- 230000014509 gene expression Effects 0.000 claims description 14
- 239000003129 oil well Substances 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000002500 effect on skin Effects 0.000 claims description 5
- 230000035515 penetration Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims 1
- 206010017076 Fracture Diseases 0.000 description 160
- 208000010392 Bone Fractures Diseases 0.000 description 134
- 238000005259 measurement Methods 0.000 description 17
- 238000004364 calculation method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 239000011324 bead Substances 0.000 description 2
- 230000005465 channeling Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005325 percolation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 208000013201 Stress fracture Diseases 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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 DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明公开了一种试井解释方法、装置、存储介质以及计算机设备，通过基于油藏参数针对包括溶孔和溶洞与裂缝组成的串珠状缝洞模式构建裂缝‑溶孔型复合地层的试井解释模型，根据该试井解释模型求取试井的实空间井底压力解并利用实空间井底压力解绘制实空间试井理论曲线，通过将实空间试井理论曲线与试井实测数据进行拟合，最后可以利用得到的拟合曲线进行试井解释。该方法通过结合地层性质，综合考虑地层中的溶孔，提高对地层中串珠状缝洞模式描述的精度，可以有利于确定出大尺寸溶洞、裂缝的分布状况，并能够直接获取到溶洞、裂缝和溶孔的相关物性参数，有利于提高试井解释的精度，方法简单、便捷且更贴近实际的地层情况。
Description
技术领域
本发明涉及油气开发技术领域，尤其涉及一种试井解释方法、装置、存储介质以及计算机设备。
背景技术
在现有技术中，针对缝洞型碳酸盐岩油藏的试井方法主要有：
一种是基于等效连续介质渗流理论的连续介质模型。该模型把整个溶洞看成是一个假想的连续体，重点研究介质整体所表现出来的宏观流动特征，基于流量等效原理，将裂缝、溶洞的渗透性等效平均到整个介质中，利用现有的渗流理论来研究流体流动特征。适用于只有溶蚀小孔洞、均匀网络状分布破碎裂缝、表征单元体存在且尺度较小的情况，多看成是双重介质或三重介质地层，不能用于描述大尺度溶洞和裂缝存在的缝洞型油藏储层特征。
另一种是基于非连续介质理论所建立的离散缝洞模型。相对连续介质来说，这类模型考虑了储集介质的具体形态和尺度等特征，在一定程度上较真实的描述了裂缝的地质分布特征及其对流体流动的影响。由于其模型的建立过程复杂，计算过程中网格需要单独生成，计算量大，是一种对油藏描述技术依赖性很强的油藏数值模拟方法，也难以与当前地质建模结果相结合。当有多个大尺度裂缝或者溶洞存在时，需要付出的计算代价很大，考虑到油田实际应用情况，离散模型不适用于缝洞型油藏试井模型建立。
目前的试井解释方法多是基于等效连续介质渗流理论，没有考虑地层中除大尺度缝、洞之外的溶孔等储层情况，对于存在大溶洞和大裂缝组成的串珠状缝洞型油藏是不适用的，其理论曲线难以与的实测曲线相互拟合，试井解释结果误差大。
发明内容
本发明要解决的技术问题是：现有技术中试井解释结果误差较大，如何提高试井解释精度的问题。
为解决上述技术问题，本发明提供了一种试井解释方法、装置、存储介质以及计算机设备。
本发明的第一个方面，提供了一种试井解释方法，其包括：
基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
根据所述试井解释模型，求取试井的实空间井底压力解；
利用所述实空间井底压力解绘制实空间试井理论曲线；
将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
通过所述拟合曲线进行试井解释。
优选的，所述基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，还包括：
基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型；
根据所述地层渗流模型和所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
优选的，所述方法还包括：根据所述裂缝-溶孔型复合地层的渗透强度进行地层区域划分；
所述基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型，包括：针对划分后的每个地层区域，分别建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型，其中，每个地层区域中裂缝和溶孔的地层渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，pf表示裂缝压力，r表示径向距离，Q表示地层区域且Q为正整数，α表示形状因子，Km表示基质渗透率，pm表示基质压力，φf表示裂缝孔隙度，Ctf表示裂缝压缩系数，t表示时间，φm表示基质孔隙度，Ctm表示基质压缩系数。
优选的，所述基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力，建立所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，包括：
根据钻遇溶洞周围的其余溶洞通过裂缝向所述钻遇溶洞窜流的过程、连接条件以及所述钻遇溶洞周围其余各个溶洞的渗流模型，建立所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，其中，所述串珠状缝洞模式的流体渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，W表示裂缝宽度，pv表示裂缝压力，pf表示裂缝压力，r表示径向距离，L表示裂缝长度，R表示溶洞半径，h表示有效厚度，Ctv表示溶洞压缩系数，φv表示溶洞孔隙度，t表示时间，m表示裂缝的总数且m＝1，2，……，j……。
优选的，根据所述地层渗流模型和所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型，包括：
根据各个地层区域中的所述地层渗流模型，结合各个地层区域的连接条件、约束边界条件、初始条件以及所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
优选的，所述方法还包括：对所述裂缝-溶孔型复合地层的试井解释模型进行无因次化，以获取无因次试井解释模型。
优选的，根据所述试井解释模型，求取试井的实空间井底压力解，包括：
对所述无因次试井解释模型进行拉普拉斯变换，以得到拉式空间的井底压力解；
利用Stehfest数值反演法将拉式空间的井底压力解反演到实空间，获取所述试井的实空间井底压力解。
优选的，对所述无因次试井解释模型进行拉普拉斯变换，以得到拉式空间的井底压力解之后，还包括：
叠加井筒储集系数和表皮系数，以得到综合井筒储集和表皮效应后的拉式空间井底压力解。
优选的，所述实空间试井理论曲线包括所述实空间井底压力解的双对数理论曲线；所述试井实测数据包括压力恢复数据。
本发明的第二个方面，提供了一种试井解释装置，其包括：
模型构建模块，用于基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
数据分析模块，用于根据所述试井解释模型，求取试井的实空间井底压力解；以及利用所述实空间井底压力解绘制实空间试井理论曲线；
曲线拟合模块，用于将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
试井解释模块，用于通过所述拟合曲线进行试井解释。
本发明的第三个方面提供了一种存储介质，该存储介质中存储有计算机程序，该计算机程序被处理器执行时能够实现如上所述的试井解释方法。
本发明的第四个方面提供了一种计算机设备，其包括：处理器以及存储介质，该存储介质中存储有计算机程序，该计算机程序被所述处理器执行时能够实现如上所述的试井解释方法。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
应用本发明的试井解释方法、装置、存储介质以及计算机设备，通过基于油藏参数针对包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，构建裂缝-溶孔型复合地层的试井解释模型，根据该试井解释模型求取试井的实空间井底压力解并利用实空间井底压力解绘制实空间试井理论曲线，通过将实空间试井理论曲线与试井实测数据进行拟合，最后可以利用得到的拟合曲线进行试井解释。该方法通过结合地层性质，综合考虑地层中的溶孔，提高对地层中串珠状缝洞模式描述的精度，可以有利于确定出大尺寸溶洞、裂缝的分布状况，并能够直接获取到溶洞、裂缝和溶孔的相关物性参数，有利于提高试井解释的精度，方法简单、便捷且更贴近实际的地层情况。
附图说明
通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
图1示出了本申请实施例提供的一种试井解释方法的流程示意图；
图2示出了本申请实施例提供的构建裂缝-溶孔型复合地层的试井解释模型的方法流程示意图；
图3示出了本申请实施例提供的一种试井解释方法的流程示意图；
图4(1)示出了本申请实施例提供的井钻遇单个溶洞的缝洞模式示意图；图4(2)示出了本申请实施例提供的串珠状缝洞模式的结构示意图；
图5示出了本申请实施例提供的求取实空间的井底压力解的方法流程示意图；
图6示出了本申请实施例提供的实空间井底压力解的双对数理论曲线示意图；
图7示出了本申请实施例提供的实空间井底压力解的双对数理论曲线与试井实测压力恢复数据的双对数曲线的拟合曲线示意图；
表1为本申请实施例提供的实测井基础参数数据表；
表2为本申请实施例提供的实测井解释结果数据表。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
在现有技术中，在现有技术中，针对缝洞型碳酸盐岩油藏的试井方法主要有：
一种是基于等效连续介质渗流理论的连续介质模型。该模型把整个溶洞看成是一个假想的连续体，重点研究介质整体所表现出来的宏观流动特征，基于流量等效原理，将裂缝、溶洞的渗透性等效平均到整个介质中，利用现有的渗流理论来研究流体流动特征。适用于只有溶蚀小孔洞、均匀网络状分布破碎裂缝、表征单元体存在且尺度较小的情况，多看成是双重介质或三重介质地层，不能用于描述大尺度溶洞和裂缝存在的缝洞型油藏储层特征。
另一种是基于非连续介质理论所建立的离散缝洞模型。相对连续介质来说，这类模型考虑了储集介质的具体形态和尺度等特征，在一定程度上较真实的描述了裂缝的地质分布特征及其对流体流动的影响。由于其模型的建立过程复杂，计算过程中网格需要单独生成，计算量大，是一种对油藏描述技术依赖性很强的油藏数值模拟方法，也难以与当前地质建模结果相结合。当有多个大尺度裂缝或者溶洞存在时，需要付出的计算代价很大，考虑到油田实际应用情况，离散模型不适用于缝洞型油藏试井模型建立。
目前的试井解释方法多是基于等效连续介质渗流理论，没有考虑地层中除大尺度缝、洞之外的溶孔等储层情况，对于存在大溶洞和大裂缝组成的串珠状缝洞型油藏是不适用的，其理论曲线难以与的实测曲线相互拟合，试井解释结果误差大。因此，需要基于非连续介质理论建立离散缝洞模型，考虑流体的流动机理，描述井周围储容性及尺寸各不相同的溶洞和不同导流能力、不同长度的裂缝与溶洞组成的缝洞模式，建立更符合实际的缝洞型碳酸盐岩油藏特点的试井解释方法。
有鉴于此，本申请提供了一种试井解释方法、装置、存储介质以及计算机设备，通过基于油藏参数针对包括溶孔和溶洞与裂缝组成的串珠状缝洞模式构建裂缝-溶孔型复合地层的试井解释模型，根据该试井解释模型求取试井的实空间井底压力解并利用实空间井底压力解绘制实空间试井理论曲线，通过将实空间试井理论曲线与试井实测数据进行拟合，最后可以利用得到的拟合曲线进行试井解释。该方法通过结合地层性质，综合考虑地层中的溶孔，提高对地层中串珠状缝洞模式描述的精度，可以有利于确定出大尺寸溶洞、裂缝的分布状况，并能够直接获取到溶洞、裂缝和溶孔的相关物性参数，有利于提高试井解释的精度，方法简单、便捷且更贴近实际的地层情况。
实施例一
参见图1所示，图1示出了本申请实施例提供的一种试井解释方法，其包括步骤S101至步骤S105。
在步骤S101中：基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型。
其中，本领域技术人员可知的是溶孔和溶洞的尺度差别很大，在本申请中综合考虑了地层中溶孔和溶洞与裂缝组成的串珠状缝洞模式。
在本申请实施例中，参见图2所示，图2示出了本申请实施例提供的构建裂缝-溶孔型复合地层的试井解释模型的方法流程示意图，具体可以通过步骤S1011至步骤S1013实现：
步骤S1011：基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型。
步骤S1012：基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力，建立裂缝-溶孔型复合地层中串珠状缝洞模式的流体渗流模型。
步骤S1013：根据地层渗流模型和流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
在步骤S1011中，地层渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，pf表示裂缝压力，r表示径向距离，Q表示地层区域且Q为正整数，α表示形状因子，Km表示基质渗透率，pm表示基质压力，φf表示裂缝孔隙度，Ctf表示裂缝压缩系数，t表示时间，φm表示基质孔隙度，Ctm表示基质压缩系数。
作为一优选示例，在构建裂缝-溶孔型复合地层的试井解释模型之前还可以根据裂缝-溶孔型复合地层的渗透强度进行地层区域划分，具体的，可以根据储层性质，将地层划分为渗透性较强的1区和渗透性较弱的2区两个地层区域，再针对划分后的每个地层区域，分别建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型，则1区中裂缝和溶孔的地层渗流模型可以表示为：
2区中裂缝和溶孔的地层渗流模型可以表示为：
井钻遇溶洞的质量守恒方程：
其中，pw表示井底压力，C表示井筒储集系数，
溶洞壁与地层之间的衔接条件：
其中，rw表示井筒半径，S表示表皮系数，
在界面流量相等条件下的连接条件：
在界面压力相等条件下的连接条件：
pf1(r1,t)＝pf2(r1,t) (1-8)
初始条件：
pj(r,t＝0)＝pi(j＝1,2,f,v,m) (1-9)
外边界条件：
pf2(r→∞,t＝0)＝pi(无限大边界) (1-10)
作为一示例，在步骤S1012中，构建串珠状缝洞模式的流体渗流模型可以具体为：
根据钻遇溶洞周围的其余溶洞通过裂缝向所述钻遇溶洞窜流的过程、连接条件以及钻遇溶洞周围其余各个溶洞的渗流模型，建立裂缝-溶孔型复合地层中串珠状缝洞模式的流体渗流模型，其中，串珠状缝洞模式的流体渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，W表示裂缝宽度，pv表示裂缝压力，pf表示裂缝压力，L表示裂缝长度，R表示溶洞半径，h表示有效厚度，Ctv表示溶洞压缩系数，φv表示溶洞孔隙度，t表示时间，m表示裂缝的总数且m＝1，2，……，j……。
具体的，当井筒钻遇溶洞1时，即钻遇单洞的缝洞模式，该钻遇溶洞的控制方程为：
其中，q表示产量，
当钻遇溶洞1通过裂缝1连接溶洞2时，溶洞1、裂缝1和溶洞2组成的串珠状缝洞模式的流体渗流方程为：
当钻遇溶洞1通过m个裂缝连接m个溶洞时，则可得该串珠状缝洞模式的流体渗流模型为：
其中，F表示连通溶洞之间的大裂缝，f表示双重介质储层中天然微裂缝。
根据以上表达式(1-1)至(1-11)可以得到裂缝-溶孔型复合地层的试井解释模型。
作为一优选示例，在本申请实施例中还可以对裂缝-溶孔型复合地层的试井模型进行无因次化，以获取无因次试井解释模型，有利于后续进行计算。
通过结合储层性质，综合考虑溶洞、溶洞与裂缝组成的串珠状缝洞模式的渗流模型，建立了更加符合实际的缝洞型油藏试井的试井解释模型，可以有效提高试井解释模型的精确度。
在步骤S102中，根据试井解释模型，求取试井的实空间井底压力解。
在步骤S103中，利用实空间井底压力解绘制实空间试井理论曲线。
在本申请实施例中，对基于裂缝和溶孔的地层渗流模型以及串珠状缝洞模式的流体渗流模型构建的试井解释模型进行拉普拉斯变换，求解获得拉普拉斯空间下的井底压力解。再将求解得到的拉普拉斯空间井底压力解通过Stehfest数值反演法，计算得到实空间的井底压力解，并利用实空间井底压力解绘制实空间试井理论曲线。其中，实空间试井理论曲线可以为实空间井底压力解的双对数理论曲线。
作为一优选示例，为了更全面的考虑储层性质，提高对试井解释的精确度，该步骤还可以具体为，对基于裂缝和溶孔的地层渗流模型以及串珠状缝洞模式的流体渗流模型构建的试井解释模型进行拉普拉斯变换后，叠加井筒储集系数和表皮系数，以得到综合井筒储集和表皮效应后的拉式空间井底压力解，从而可以综合考虑井筒中流体的存储释放以及井底周围受污染的情况。
在步骤S104中，将实空间试井理论曲线与试井实测数据进行拟合，以得到实空间试井理论曲线与试井实测数据的拟合曲线。
在本申请实施例中，试井实测数据可以为实测试井的压力恢复数据，将实空间试井理论曲线与试井实测数据进行拟合可以为，将实空间井底压力解的双对数理论曲线与实测试井的压力恢复数据进行拟合，具体的，在实空间井底压力解的双对数坐标系中，绘制实测试井的压力恢复及压力恢复倒数的双对数曲线，利用基础参数，通过不断调整油藏参数及拟合参数，反复调用正演过程，使实空间井底压力解的双对数理论曲线与实测试井的双对数曲线相互拟合，以得到最终的拟合曲线。
在步骤S105中，通过拟合曲线进行试井解释。
具体的，根据拟合曲线可以计算裂缝渗透率Kf1、溶洞体积V、窜流系数λ和溶洞储容比w等油藏参数，从而可以得到试井解释结果数据，提供了一种综合考虑储层性质的试井解释方法。
以上为本申请实施例一提供的一种试井解释方法，基于油藏参数构建裂缝-溶孔型复合地层的试井解释模型，通过结合储层性质，考虑了缝洞型储层中溶孔的特性以及溶洞与裂缝构成的串珠状模式，从而为确定大尺度溶洞、裂缝的分布状况及溶孔参数提供了完整的技术手段，能够直接获得溶洞大小、溶洞体积、溶洞储容比等参数，所得的结果与地质认识相结合，能够提高地层描述精度，有利于获取更多的缝洞结构参数信息，使试井解释结果数据更加贴近实际。
实施例二
参见图3所示，图3示出了本申请实施例提供的一种试井解释方法的流程示意图，其包括步骤S301至步骤S308。
在步骤S301中，根据复合地层的渗透强度进行地层区域划分。
作为示例，根据图4(1)所示，图4(1)示出了井钻遇单个溶洞的缝洞模式示意图，可以根据复合地层的渗透强度将地层区域划分，具体的，可以根据形成溶洞的过程中，溶洞周围的储层性质，将地层区域划分为两个区域，其中，溶洞周围容易坍塌形成渗透性较强的储层对应Q＝1区，1区中包含溶控、溶洞与裂缝组成的串珠状缝洞模式，远离溶洞中心渗透性较差的储层对应Q＝2区。通过结合储层性质对地层区域进行划分，有利于提高地层描述的精度以及试井解释的准确性。
在本申请实施例中，均以将地层区域划分为1区和2区两个区域为例进行说明。
步骤S302，基于油藏参数，针对划分后的每个地层区域，分别建立裂缝和溶孔的地层渗流模型。
步骤S303，建立裂缝-溶孔型复合地层中串珠状缝洞模式的流体渗流模型。
步骤S304，对裂缝-溶孔型复合地层的试井解释模型进行无因次化，以获取无因次试井解释模型。
在本申请实施例中，在步骤S302中，可以通过定义无因此量得到无因次化后1区中裂缝和溶孔的地层渗流模型：
无因次化后得到的2区中裂缝和溶孔的地层渗流模型：
其次，分别约束1区和2区的连接条件、边界条件及初始条件。
作为一示例，可以在假设界面流量相等条件下，使用如下表达式作为连接条件：
作为另一示例，可以在假设界面压力相等条件下，使用如下表达式作为连接条件：
pf1D(r1D,tD)＝pf2D(r1D,tD)(2-6)
边界条件可以为无限大边界条件：pf2D(rD→∞,tD＝0)＝0 (2-7)
初始条件可以为：pjD(rD,tD＝0)＝0(j＝1,2,f,v,w,m)(2-8)
其中，无因次压力：
无因次时间：
无因次井筒储集系数；
无因次距离：其中reD表示有效边界距离；
1区裂缝储容比：2区裂缝储容比：/>
溶洞储容比：1区窜流系数：/>α代表性质因子；
2区窜流系数：1区和2区储容比：/>
1区和2区流度比：无因次溶洞体积：/>
在步骤S303中，可以通过考虑储层性质，假设溶洞中各点压力均相等，利用物质平衡方程来描述溶洞中流体的流动。
作为一示例，当井筒为一个，井钻遇溶洞1，溶洞1周围的裂缝-溶孔型复合地层中存在m个溶洞，且该m个溶洞通过m条裂缝相互连接组成串珠状的缝洞模式。
针对该串珠状缝洞模式，可以按照各个溶洞向着溶洞1窜流的过程、连接条件以及各个溶洞的控制方程来建立串珠状缝洞模式的流体渗流模型，具体的建立过程可以为：
当井筒钻遇溶洞1时，该单溶洞的渗流表达式可以为：
当钻遇溶洞1通过裂缝1与溶洞2连接时：
……
当钻遇溶洞1通过m个裂缝连接m个溶洞时，参见图4(2)所示，图4(2)示出了本申请实施例提供的串珠状缝洞模式的结构示意图，该串珠状缝洞模式的渗流方程可以表示为：
步骤S305：根据无因次试井解释模型，求取实空间井底压力解。
在本申请实施例中，参见图5所示，图5示出了本申请实施例提供的求取实空间的井底压力解的方法流程示意图，步骤S305可以通过以下步骤实现：
在步骤S3051中，对无因次试井解释模型进行拉普拉斯变换，以得到拉式空间的井底压力解。
作为一优选示例，在对无因次试井解释模型进行拉普拉斯变换之后，解得无限大地层边界条件下拉普拉斯空间井底压力解，在所求压力解得基础上，叠加无因次井筒储集系数CD和表皮系数S，得到综合井筒储集和表皮效应后的拉式空间井底压力解：
其中，u代表拉普拉斯变量，表达式(3-1)对应单洞模式下的拉式空间井底压力解，表达式(3-2)对应串珠状模式下的拉式空间井底压力解，其中，
其中，
第一类虚宗量0阶贝塞尔函数；
第一类虚宗量1阶贝塞尔函数；
第二类虚宗量0阶贝塞尔函数；/>
第二类虚宗量1阶贝塞尔函数；
其中，
在步骤S3052中，利用Stehfest数值反演法将拉式空间的井底压力解反演到实空间，获取实空间井底压力解。
在本申请实施例中，利用Stehfest数值反演法对拉式空间的井底压力解进行反演后，可以得到井底压力解的表达式如下：
其中，其中，Z为Laplace因子；j为整数变量。
Vj是常数，取决于N值，用下面方程确定：
其中，N为偶数，一般取值在8—16之间。
需要说明的是，在本申请实施例中只是示例性的给出了无限大边界条件下的拉普拉斯空间井底压力解，本申请对边界条件并不做特别的限定，也可以通过其他边界条件获取拉普拉斯空间井底压力解。
在步骤S306中，利用实空间井底压力解绘制实空间试井理论曲线。
在本申请实施例中，利用实空间井底压力解绘制实空间试井理论曲线，其中，实空间试井理论曲线可以为实空间井底压力解的双对数理论曲线，作为一示例，串珠状模式下的时空间井底压力解的双对数理论曲线可以如图6所示，图6示出了本申请实施例提供的实空间井底压力解的双对数理论曲线示意图。
步骤S307：将实空间试井理论曲线与试井实测数据进行拟合，以获取空间试井理论曲线与试井实测数据的拟合曲线。
其中，试井实测数据可以为试井实测的压力恢复数据。
在本申请实施例中，该步骤可以具体为，在获取到实空间井底压力解的双对数理论曲线后，在同一个双对数坐标系中，利用试井实测的压力恢复数据，绘制试井实测的压力及压力导数的双对数曲线。利用基础参数数据，通过不断调整油藏参数及拟合参数，反复调用正演过程，使实空间井底压力解的双对数理论曲线与试井实测的双对数曲线相互拟合，以得到实空间井底压力解的双对数理论曲线与试井实测的双对数曲线的最终拟合图，参见图7所示，图7示出了本申请实施例提供的实空间井底压力解的双对数理论曲线与试井实测压力恢复数据的双对数曲线的拟合曲线示意图。
步骤S308：通过实空间试井理论曲线与试井实测数据的拟合曲线进行试井解释。
在本申请实施例中，可以根据实空间井底压力解的双对数理论曲线与试井实测的双对数曲线的最终拟合图获取到如下多个拟合参数，其中下标M代表拟合得到的数据：
压力拟合值：
时间拟合值：
无因次井筒储集系数：表皮系数：(S)M；
无因次溶洞半径：
1区窜流系数：2区窜流系数：/>
12区储容比：12区流度比：/>
1区裂缝储容比：
2区裂缝储容比：
溶洞储容比：无因次1区半径：/>
进一步地，根据获取到的拟合参数可以相应得到多个油藏参数，从而得到试井解释结果数据。例如：
可以根据压力拟合值，计算出1区裂缝渗透率Kf1：
可以利用时间拟合值和无因次井筒储集系数，计算出井筒储集系数C：
可以利用无因次溶洞体积，计算出溶洞体积V：
可以利用无因次1区半径，计算出1区半径r1：r1＝rw(r1D)M；
可以利用表皮系数、溶洞储容比、1区裂缝储容比、2区裂缝储容比、1区窜流系数、2区窜流系数：
1区窜流系数：λf1＝(λf1)M；2区窜流系数：λf2＝(λf2)M；
1区裂缝储容比：ωf1＝(ωf1)M；2区裂缝储容比：ωf2＝(ωf2)M；
溶洞储容比：ω1＝(ω1)M；表皮系数：S＝(S)M；
可以利用12区流度比，计算出2区流度M2：
可以得到的12区储容比，计算出2区储容：
以上为本申请实施例提供的一种试井解释方法，该方法中通过基于油藏参数构建裂缝-溶孔型复合地层的试井解释模型，通过结合储层性质，考虑了缝洞型储层中溶孔的特性以及大尺度溶洞与裂缝构成的串珠状缝洞模式，并根据渗透强度对地层进行区域划分，以针对不同地层区域建立渗流模型，从而为确定大尺度溶洞的分布状况及洞外溶孔参数提供了完整的技术手段，能够直接获得溶洞大小、溶洞体积、溶洞储容比等参数，所得的结果与地质认识相结合，能够提高地层描述精度并获取更多的缝洞结构参数信息，使试井解释结果数据更加贴近实际，有利于有效开发缝洞型碳酸盐油藏。
另外，本发明中提供的试井解释模型建立简单，根据该试井解释模型求取在拉普拉斯空间的解析解，不会涉及到复杂函数的计算，方法简便、有利于提高计算速度。
实施例三
以下为本申请提供的一应用试井解释方法的具体示例，实测井的设计井深5730.0m，完钻井深5754.0m，人工井底5520.0m，完钻层位奥陶系，裸眼酸压完井，压恢测试时间112.88h。钻井过程中发生漏失，漏失量33.3m3，无水采油期154d，无水采油量3.89万吨。实例井基本参数数据见表1。
表1实测井基础参数数据
利用实例井的基础参数数据，通过不断调整油藏参数及拟合参数，反复调用正演过程，使实空间井底压力解的双对数理论曲线与实测井的双对数曲线相互拟合，曲线拟合之后得到实空间双对数理论曲线与实测井双对数曲线的最终拟合结果，如图7所示。
利用实空间井底压力解的双对数理论曲线与实测压力恢复数据的最终拟合结果，通过计算油藏参数，获得裂缝导流能力、裂缝长度、溶洞体积及溶洞储容比等试井解释结果数据，见表2。
表2实测井解释结果数据
以上为应用本申请实施例提供的试井解释方法的具体示例，基于构建的裂缝-溶孔型复合地层的试井解释模型，通过实空间井底压力解的双对数理论曲线与实测压力恢复数据进行拟合后的拟合曲线，获得溶洞体积和溶洞储容比等试井解释结果数据，通过应用本申请实施例提供的试井解释方法，能够直接解释出地层和溶洞的相关参数，求解方法简单快捷。
实施例四
本申请的另一个方面还提供了一种试井解释装置，其包括：
模型构建模块，用于基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
数据分析模块，用于根据所述试井解释模型，求取试井的实空间井底压力解；以及利用所述实空间井底压力解绘制实空间试井理论曲线；
曲线拟合模块，用于将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
试井解释模块，用于通过所述拟合曲线进行试井解释。其中，本领域技术人员可知的是溶孔和溶洞的尺度差别很大，在本申请中综合考虑了地层中溶孔和溶洞与裂缝组成的串珠状缝洞模式。
在本申请实施例中，参见图2所示，图2示出了本申请实施例提供的构建裂缝-溶孔型复合地层的试井解释模型的方法流程示意图，具体可以通过步骤S1011至步骤S1013实现：
步骤S1011：基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型。
步骤S1012：基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力，建立裂缝-溶孔型复合地层中串珠状缝洞模式的流体渗流模型。
步骤S1013：根据所述地层渗流模型和所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
作为一优选示例，在构建裂缝-溶孔型复合地层的试井解释模型之前还可以根据裂缝-溶孔型复合地层的渗透强度进行地层区域划分，具体的，可以根据储层性质，将地层划分为渗透性较强的1区和渗透性较弱的2区两个地层区域，再针对划分后的每个地层区域，分别建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型。
根据钻遇溶洞周围的其余溶洞通过裂缝向钻遇溶洞窜流的过程、连接条件以及钻遇溶洞周围其余各个溶洞的渗流模型，建立裂缝-溶孔型复合地层中述串珠状缝洞模式的流体渗流模型。
作为一优选示例，在本申请实施例中还可以对裂缝-溶孔型复合地层的试井模型进行无因次化，以获取无因次试井解释模型，有利于后续进行计算。
通过结合储层性质，综合考虑溶洞、溶洞与裂缝组成的串珠状缝洞模式的渗流模型，建立了更加符合实际的缝洞型油藏试井的试井解释模型，可以有效提高试井解释模型的精确度。
在本申请实施例中，对基于裂缝和溶孔的地层渗流模型以及串珠状缝洞模式的流体渗流模型构建的试井解释模型进行拉普拉斯变换，求解获得拉普拉斯空间下的井底压力解。再将求解得到的拉普拉斯空间井底压力解通过Stehfest数值反演法，计算得到实空间的井底压力解，并利用实空间井底压力解绘制实空间试井理论曲线。其中，实空间试井理论曲线可以为实空间井底压力解的双对数理论曲线。
作为一优选示例，为了更全面的考虑储层性质，提高对试井解释的精确度，该步骤还可以具体为，对基于裂缝和溶孔的地层渗流模型以及串珠状缝洞模式的流体渗流模型构建的试井解释模型进行拉普拉斯变换后，叠加井筒储集系数和表皮系数，以得到综合井筒储集和表皮效应后的拉式空间井底压力解，从而可以综合考虑井筒中流体的存储释放以及井底周围受污染的情况。
在本申请实施例中，试井实测数据可以为实测试井的压力恢复数据，将实空间试井理论曲线与试井实测数据进行拟合可以为，将实空间井底压力解的双对数理论曲线与实测试井的压力恢复数据进行拟合，具体的，在实空间井底压力解的双对数坐标系中，绘制实测试井的压力恢复及压力恢复倒数的双对数曲线，利用基础参数，通过不断调整油藏参数及拟合参数，反复调用正演过程，使实空间井底压力解的双对数理论曲线与实测试井的双对数曲线相互拟合，以得到最终的拟合曲线。
具体的，根据拟合曲线可以计算裂缝渗透率Kf1、溶洞体积V、窜流系数λ和溶洞储容比w等油藏参数，从而可以得到试井解释结果数据，提供了一种综合考虑储层性质的试井解释方法。
实施例五
本申请的另一个方面还提供了一种存储介质，该存储介质中存储有计算机程序，该计算机程序被处理器执行时能够实现如下所述的试井解释方法：
基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
根据所述试井解释模型，求取试井的实空间井底压力解；
利用所述实空间井底压力解绘制实空间试井理论曲线；
将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
通过所述拟合曲线进行试井解释。
本申请的另一方面，还提供了一种计算机设备，其包括：处理器以及存储介质，该存储介质中存储有计算机程序，该计算机程序被处理器执行时能够实现如下所述的试井解释方法：
基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
根据所述试井解释模型，求取试井的实空间井底压力解；
利用所述实空间井底压力解绘制实空间试井理论曲线；
将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
通过所述拟合曲线进行试井解释。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (11)
1.一种试井解释方法，其特征在于，包括：
基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
根据所述试井解释模型，求取试井的实空间井底压力解；
利用所述实空间井底压力解绘制实空间试井理论曲线；
将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
通过所述拟合曲线进行试井解释；
所述基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，包括：
根据钻遇溶洞周围的其余溶洞通过裂缝向所述钻遇溶洞窜流的过程、连接条件以及所述钻遇溶洞周围其余各个溶洞的渗流模型，建立所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，其中，所述串珠状缝洞模式的流体渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，W表示裂缝宽度，pv表示裂缝压力，pf表示裂缝压力，r表示径向距离，L表示裂缝长度，R表示溶洞半径，h表示有效厚度，Ctv表示溶洞压缩系数，φv表示溶洞孔隙度，t表示时间，m表示裂缝的总数且m＝1，2，……，j……。
2.根据权利要求1所述的方法，其特征在于，所述基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，还包括：
基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型；
根据所述地层渗流模型和所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
3.根据权利要求2所述的方法，其特征在于，所述方法还包括：根据所述裂缝-溶孔型复合地层的渗透强度进行地层区域划分；
所述基于流体粘度、径向距离、时间、形状因子以及裂缝和基质的渗透率、压力、孔隙度、压缩系数，建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型，包括：针对划分后的每个地层区域，分别建立裂缝-溶孔型复合地层中裂缝和溶孔的地层渗流模型，其中，每个地层区域中裂缝和溶孔的地层渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，pf表示裂缝压力，r表示径向距离，Q表示地层区域且Q为正整数，α表示形状因子，Km表示基质渗透率，pm表示基质压力，φf表示裂缝孔隙度，Ctf表示裂缝压缩系数，t表示时间，φm表示基质孔隙度，Ctm表示基质压缩系数。
4.根据权利要求3所述的方法，其特征在于，根据所述地层渗流模型和所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型，包括：
根据各个地层区域中的所述地层渗流模型，结合各个地层区域的连接条件、约束边界条件、初始条件以及所述流体渗流模型，构建裂缝-溶孔型复合地层的试井解释模型。
5.根据权利要求1所述的方法，其特征在于，所述方法还包括：对所述裂缝-溶孔型复合地层的试井解释模型进行无因次化，以获取无因次试井解释模型。
6.根据权利要求5所述的方法，其特征在于，根据所述试井解释模型，求取试井的实空间井底压力解，包括：
对所述无因次试井解释模型进行拉普拉斯变换，以得到拉式空间的井底压力解；
利用Stehfest数值反演法将拉式空间的井底压力解反演到实空间，获取所述试井的实空间井底压力解。
7.根据权利要求6所述的方法，其特征在于，对所述无因次试井解释模型进行拉普拉斯变换，以得到拉式空间的井底压力解之后，还包括：
叠加井筒储集系数和表皮系数，以得到综合井筒储集和表皮效应后的拉式空间井底压力解。
8.根据权利要求1所述的方法，其特征在于，所述实空间试井理论曲线包括所述实空间井底压力解的双对数理论曲线；所述试井实测数据包括压力恢复数据。
9.一种试井解释装置，其特征在于，包括：
模型构建模块，用于基于油藏参数，构建裂缝-溶孔型复合地层的试井解释模型，其中，所述裂缝-溶孔型复合地层中包括溶孔和溶洞与裂缝组成的串珠状缝洞模式，所述试井解释模型包括基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型；
数据分析模块，用于根据所述试井解释模型，求取试井的实空间井底压力解；以及利用所述实空间井底压力解绘制实空间试井理论曲线；
曲线拟合模块，用于将所述实空间试井理论曲线与试井实测数据进行拟合，以得到所述实空间试井理论曲线与所述试井实测数据的拟合曲线；
试井解释模块，用于通过所述拟合曲线进行试井解释；
所述基于裂缝长度、裂缝宽度、有效厚度、流体粘度、径向距离、时间、孔隙度以及裂缝、溶洞和油井的压力建立的所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，包括：
根据钻遇溶洞周围的其余溶洞通过裂缝向所述钻遇溶洞窜流的过程、连接条件以及所述钻遇溶洞周围其余各个溶洞的渗流模型，建立所述裂缝-溶孔型复合地层中所述串珠状缝洞模式的流体渗流模型，其中，所述串珠状缝洞模式的流体渗流模型的表达式为：
其中，μ表示流体粘度，Kf表示裂缝渗透率，W表示裂缝宽度，pv表示裂缝压力，pf表示裂缝压力，r表示径向距离，L表示裂缝长度，R表示溶洞半径，h表示有效厚度，Ctv表示溶洞压缩系数，φv表示溶洞孔隙度，t表示时间，m表示裂缝的总数且m＝1，2，……，j……。
10.一种存储介质，其特征在于，该存储介质中存储有计算机程序，所述计算机程序被处理器执行时能够实现如权利要求1至8中任一项所述的试井解释方法。
11.一种计算机设备，其特征在于，包括：处理器以及存储介质，该存储介质中存储有计算机程序，该计算机程序被处理器执行时能够实现如权利要求1至8中任一项所述的试井解释方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010767808.8A CN114086933B (zh) | 2020-08-03 | 2020-08-03 | 一种试井解释方法、装置、存储介质以及计算机设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010767808.8A CN114086933B (zh) | 2020-08-03 | 2020-08-03 | 一种试井解释方法、装置、存储介质以及计算机设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114086933A CN114086933A (zh) | 2022-02-25 |
| CN114086933B true CN114086933B (zh) | 2024-03-19 |
Family
ID=80295226
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010767808.8A Active CN114086933B (zh) | 2020-08-03 | 2020-08-03 | 一种试井解释方法、装置、存储介质以及计算机设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114086933B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106150477A (zh) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | 一种确定单井控制储量的方法 |
| WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110096718A (zh) * | 2018-01-30 | 2019-08-06 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中溶洞的体积的方法 |
| CN110532592A (zh) * | 2019-07-15 | 2019-12-03 | 西南石油大学 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
| CN110593865A (zh) * | 2019-09-29 | 2019-12-20 | 中国石油集团川庆钻探工程有限公司 | 一种油藏缝洞特征参数试井解释方法 |
| CN110984973A (zh) * | 2019-11-29 | 2020-04-10 | 中国石油集团川庆钻探工程有限公司 | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| MX346226B (es) * | 2012-03-30 | 2017-03-07 | Inst Mexicano Del Petróleo | Método de análisis integral de pruebas de trazadores entre pozos petroleros. |
| US20180179891A1 (en) * | 2016-01-04 | 2018-06-28 | Saudi Arabian Oil Company | Measuring effective fracture half-length and quantifying flux distribution in and around fractures in petroleum reservoirs |
-
2020
- 2020-08-03 CN CN202010767808.8A patent/CN114086933B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106150477A (zh) * | 2015-04-23 | 2016-11-23 | 中国石油化工股份有限公司 | 一种确定单井控制储量的方法 |
| WO2016192077A1 (zh) * | 2015-06-04 | 2016-12-08 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种致密气压裂水平井数值试井模型建立求解方法 |
| CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
| CN110096718A (zh) * | 2018-01-30 | 2019-08-06 | 中国石油化工股份有限公司 | 一种获取碳酸盐岩油藏中溶洞的体积的方法 |
| CN110532592A (zh) * | 2019-07-15 | 2019-12-03 | 西南石油大学 | 一种缝洞油气藏压裂井大溶洞试井解释方法 |
| CN110593865A (zh) * | 2019-09-29 | 2019-12-20 | 中国石油集团川庆钻探工程有限公司 | 一种油藏缝洞特征参数试井解释方法 |
| CN110984973A (zh) * | 2019-11-29 | 2020-04-10 | 中国石油集团川庆钻探工程有限公司 | 缝洞型碳酸盐岩气藏单井控制储量的确定方法 |
Non-Patent Citations (3)
| Title |
|---|
| "串珠"状地震反射与试井分析"串珠"模型关系探讨;胥洪俊;孙贺东;张万能;王峰;王明;吕波;骆蕾;高上钧;;油气井测试(05);第6-9页 * |
| 串珠状缝洞型碳酸盐岩油藏数值试井分析;朱波宁;张磊;李丁;;辽宁化工;46(02);第164-166页 * |
| 大尺度溶洞发育的缝洞型油藏试井解释模型研究;尹洪军 等;特种油气藏;第25卷(第5期);第84-88页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114086933A (zh) | 2022-02-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104533370B (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
| CN106599449B (zh) | 一种溶洞体积计算的试井解释方法 | |
| WO2016192077A1 (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
| CN109241588B (zh) | 一种基于拟连续地质力学模型的单裂缝扩展的模拟方法 | |
| Yang et al. | A comprehensive model for real gas transport in shale formations with complex non-planar fracture networks | |
| CN111581786A (zh) | 用于分析缝洞串联模式双孔复合储层的试井解释模型 | |
| CN110096718B (zh) | 一种获取碳酸盐岩油藏中溶洞的体积的方法 | |
| CN109162701B (zh) | 一种煤层裸眼井破裂压力预测方法 | |
| CN110107277B (zh) | 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法 | |
| CN106886046B (zh) | 确定缝洞型气藏未投产区块可动用储量的方法 | |
| CN106555575A (zh) | 一种深层页岩的地层起裂压力的计算方法 | |
| CN112096362B (zh) | 非常规储层多簇射孔竞争起裂与扩展模拟方法及装置 | |
| CN116306385B (zh) | 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质 | |
| CN113255092A (zh) | 体积压裂水平井特征点拟合反演方法、电子设备及介质 | |
| CN110096669B (zh) | 一种获取碳酸盐岩油藏中裂缝体积的方法 | |
| CN112377182B (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| Sarvaramini et al. | Characterizing the stimulated reservoir volume during hydraulic fracturing-connecting the pressure fall-off phase to the geomechanics of fracturing | |
| CN114086933B (zh) | 一种试井解释方法、装置、存储介质以及计算机设备 | |
| CN112377184B (zh) | 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置 | |
| Mazo et al. | Interpretation of indicator studies of multistage fracturing | |
| CN114153007B (zh) | 水锁效应损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 | |
| CN114086934B (zh) | 一种试井解释方法、装置、存储介质以及计算机设备 | |
| CN113887110A (zh) | 油井产能预测方法、装置、电子设备及计算机存储介质 | |
| CN106930759A (zh) | 一种低渗透砂岩油藏产能的预测方法 | |
| CN114325845A (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 |