CN115496014A - 一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 - Google Patents
一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 Download PDFInfo
- Publication number
- CN115496014A CN115496014A CN202211432710.2A CN202211432710A CN115496014A CN 115496014 A CN115496014 A CN 115496014A CN 202211432710 A CN202211432710 A CN 202211432710A CN 115496014 A CN115496014 A CN 115496014A
- Authority
- CN
- China
- Prior art keywords
- crack
- well
- fracture
- water
- data
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 124
- 238000004088 simulation Methods 0.000 title claims abstract description 79
- 239000007924 injection Substances 0.000 title claims abstract description 67
- 238000002347 injection Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000009826 distribution Methods 0.000 claims abstract description 42
- 238000012360 testing method Methods 0.000 claims abstract description 27
- 238000005094 computer simulation Methods 0.000 claims abstract description 20
- 238000010586 diagram Methods 0.000 claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000003384 imaging method Methods 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 30
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 239000002245 particle Substances 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 238000010276 construction Methods 0.000 abstract description 3
- 206010017076 Fracture Diseases 0.000 description 116
- 208000010392 Bone Fractures Diseases 0.000 description 78
- 239000003129 oil well Substances 0.000 description 12
- 239000012530 fluid Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 239000012634 fragment Substances 0.000 description 4
- 239000013505 freshwater Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000005012 migration Effects 0.000 description 4
- 238000013508 migration Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005429 filling process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 235000020681 well water Nutrition 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
Images
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- E21B33/00—Sealing or packing boreholes or wells
- E21B33/10—Sealing or packing boreholes or wells in the borehole
- E21B33/13—Methods or devices for cementing, for plugging holes, crevices or the like
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及裂缝型油藏分析技术领域,更具体地,涉及一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备。该方案包括通过截取获得单井的截面图,设置单井模型;选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;准备充填动态模拟数据,判断裂缝分布反映井周流动能力;判断含水率数据的模拟结果与实际的吻合程度;计算裂缝开度大小和裂缝长度的更新值。该方案通过对水注入后井周围裂缝分布规律的动态建立,以裂缝扩张和延伸模型为基础,综合分析实际充填井的成像测井数据、清水实验数据、充填施工数据、措施后生产数据,得到了靶区裂缝充填主控因素。
Description
技术领域
本发明涉及裂缝型油藏分析技术领域,更具体地,涉及一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备。
背景技术
在本发明技术之前,裂缝型油藏出水的规律的研究不足,导致目前还存在以下问题:第一方面,对于裂缝型油藏,由于其储集空间的复杂性,在油井水淹后往往关注储集空间对油井出水的影响,而在实际的堵水措施中,有效地封堵出水通道才能控制含水上升。但是,目前对于裂缝型油藏中裂缝出水规律的认识还十分有限,当同时存在宽缝和窄缝时,宽缝的存在会对窄缝中的流体运移形成一定的干扰。
第二方面,由于对流体在不同大小的裂缝中的运移规律认识不足,难以确定宽缝对窄缝的流动干扰程度及其界限。同时对于宽裂缝出水、窄裂缝出水或者宽窄裂缝都出水的情况下油井的生产动态响应特征还不清楚。因此,难以通过油井的动态资料去认识和诊断地下裂缝出水情况,还需建立典型的数值模拟模型,得到裂缝干扰出水规律并得到动态特征去指导实际油井出水规律的认识。
第三方面,由于裂缝型油藏储集空间中流体物性复杂,数值模拟表征上还存在很多问题。目前的裂缝型油藏数值模拟技术主要有等效模拟法和离散模型。等效模拟法虽然能够适应矿场尺度的研究,但是无法精确表征裂缝参数,在研究裂缝出水模拟中难以得到有效的认识规律;而离散模型虽然可以表征裂缝参数并且模拟出水规律,但是对于堵水模拟等难以与商业数值模拟软件相比。此外,对于裂缝出水规律的认识在矿场和动态上的识别和诊断目前还没有成熟的方法,需要进行油井出水规律的识别才能针对性地进行分类堵水。
发明内容
鉴于上述问题,本发明提出了一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备,通过提供一种井封隔体注入控堵水模拟方法,实现对于水注入后的井周围裂缝分布规律的动态建立,以裂缝扩张和延伸模型为基础,综合分析实际充填井的成像测井数据、清水实验数据、充填施工数据、措施后生产数据,得到了靶区裂缝充填主控因素。
根据本发明实施例第一方面,提供一种裂缝型油藏封隔体注入控堵水模拟方法。
在一个或多个实施例中,优选地,所述一种裂缝型油藏封隔体注入控堵水模拟方法包括:
通过截取获得单井的截面图,设置单井模型;
选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
在一个或多个实施例中,优选地,所述通过截取获得单井的截面图,设置单井模型,具体包括:
通过截取获得单井的截面图;
根据所述单井的截面图以井边距离为基础设置单井模型;
对所述单井模型进行水平方向的加密。
在一个或多个实施例中,优选地,所述选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片,具体包括:
选取时间较长的清水测试数据,所述时间较长的清水测试数据具体为超过预设的时间长度的测试数据;
以净注入量为控制条件,生成注水制度文件;
根据管流模拟生成注入VFP表,其中,所述VFP表为关系型数据表;
根据所述关系型数据表计算比较井口压力;
充填模拟需要离散化的裂缝;
通过成像测井得到的裂缝统计表,生成过井的裂缝片。
在一个或多个实施例中,优选地,所述在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度,具体包括:
在所述过井的裂缝片区域之外,生成随机裂缝,形成全部的裂缝基础数据;
根据所述裂缝基础数据,设置所有的所述裂缝的初始开度;
拟合清水测试中的井口压力,使得压力的拟合值与井周流动能力一致的情况下,更新调整所述裂缝基础数据中的所述裂缝分布和所述裂缝的初始开度。
在一个或多个实施例中,优选地,所述准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力,具体包括:
准备充填动态模拟数据,根据开度、导流能力、渗透率、颗粒浓度、应力的场数据变化开展充填动态模拟;
在三维视窗中查看力饱和度;
在二维曲线中查看充填压力和排量;
通过调整A井井上裂缝开度渗透率、井周裂缝纵向联通程度和密度拟合清水实验压力曲线,判断所述裂缝分布反映井周流动能力。
在一个或多个实施例中,优选地,所述获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度,具体包括:
获取当前的模拟结果,判断充填是否主要发生在井上裂缝,且距离井更远的裂缝较少波及;
判断井上裂缝开度随充填进行均逐渐增大,增大幅度从井眼处向周围逐渐减小;
判断裂缝导流能力由于开度增大一开始普遍增大,随后颗粒进入封堵,从中心向外围逐渐变小;
判断半年内含水率数据的模拟结果与实际动态趋势和幅度的吻合程度。
在一个或多个实施例中,优选地,所述获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值,具体包括:
获取所述裂缝开度大小和所述裂缝长度,利用第一计算公式计算裂缝开度变化率和裂缝长度变化率;
利用第二计算公式计算所有的相邻的裂缝之间的关联度,不相邻的裂缝之间关联度为零;
利用第三计算公式计算每个裂缝的总体关联度;
利用第四计算公式计算裂缝开度大小的更新值;
利用第五计算公式计算裂缝长度的更新值;
所述第一计算公式为:
其中,D 1_i_T 为裂缝开度变化率,D 2_i_T 为裂缝长度变化率,K 1_i_T 为裂缝开度大小,K 2_i_T 为裂缝长度;
所述第二计算公式为:
所述第三计算公式为:
其中,G i 为第i裂缝的总体关联度,n为裂缝总数,g ii_T 为0;
所述第四计算公式为:
其中,TK 1_i_T 为裂缝开度大小的更新值;
所述第五计算公式为:
其中, TK 2_i_T 为裂缝长度的更新值。
根据本发明实施例第二方面,提供一种裂缝型油藏封隔体注入控堵水模拟系统。
在一个或多个实施例中,优选地,所述一种裂缝型油藏封隔体注入控堵水模拟系统包括:
单井模型生成模块,用于通过截取获得单井的截面图,设置单井模型;
井周裂缝生成模块,用于选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
裂缝拟合模块,用于在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
充填动态模拟模块,用于准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
指标输出模块,用于获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
规则分析模块,用于获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
根据本发明实施例第三方面,提供一种计算机可读存储介质,其上存储计算机程序指令,所述计算机程序指令在被处理器执行时实现如本发明实施例第一方面中任一项所述的方法。
根据本发明实施例第四方面,提供一种电子设备,包括存储器和处理器,所述存储器用于存储一条或多条计算机程序指令,其中,所述一条或多条计算机程序指令被所述处理器执行以实现本发明实施例第一方面中任一项所述的方法。
本发明的实施例提供的技术方案可以包括以下有益效果:
本发明方案中,建立了非结构化网格双重介质油藏模型及完成历史拟合,分析了靶区裂缝发育特征描述和分布规律。
本发明方案中,通过对裂缝充填动态过程进行模拟,实现了靶区裂缝充填控堵水机理及主控因素的分析。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法的流程图。
图2是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的通过截取获得单井的截面图,设置单井模型的流程图。
图3是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片的流程图。
图4是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度的流程图。
图5是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力的流程图。
图6是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度的流程图。
图7是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值的流程图。
图8是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟系统的结构图。
图9是本发明一个实施例中一种电子设备的结构图。
具体实施方式
在本发明的说明书和权利要求书及上述附图中的描述的一些流程中,包含了按照特定顺序出现的多个操作,但是应该清楚了解,这些操作可以不按照其在本文中出现的顺序来执行或并行执行,操作的序号如101、102等,仅仅是用于区分开各个不同的操作,序号本身不代表任何的执行顺序。另外,这些流程可以包括更多或更少的操作,并且这些操作可以按顺序执行或并行执行。需要说明的是,本文中的“第一”、“第二”等描述,是用于区分不同的消息、设备、模块等,不代表先后顺序,也不限定“第一”和“第二”是不同的类型。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明技术之前,裂缝型油藏出水的规律的研究不足,导致目前还存在以下问题:第一方面,对于裂缝型油藏,由于其储集空间的复杂性,在油井水淹后往往关注储集空间对油井出水的影响,而在实际的堵水措施中,有效地封堵出水通道才能控制含水上升。但是,目前对于裂缝型油藏中裂缝出水规律的认识还十分有限,当同时存在宽缝和窄缝时,宽缝的存在会对窄缝中的流体运移形成一定的干扰。
第二方面,由于对流体在不同大小的裂缝中的运移规律认识不足,难以确定宽缝对窄缝的流动干扰程度及其界限。同时对于宽裂缝出水、窄裂缝出水或者宽窄裂缝都出水的情况下油井的生产动态响应特征还不清楚。因此,难以通过油井的动态资料去认识和诊断地下裂缝出水情况,还需建立典型的数值模拟模型,得到裂缝干扰出水规律并得到动态特征去指导实际油井出水规律的认识。
第三方面,由于裂缝型油藏储集空间中流体物性复杂,数值模拟表征上还存在很多问题。目前的裂缝型油藏数值模拟技术主要有等效模拟法和离散模型。等效模拟法虽然能够适应矿场尺度的研究,但是无法精确表征裂缝参数,在研究裂缝出水模拟中难以得到有效的认识规律;而离散模型虽然可以表征裂缝参数并且模拟出水规律,但是对于堵水模拟等难以与商业数值模拟软件相比。此外,对于裂缝出水规律的认识在矿场和动态上的识别和诊断目前还没有成熟的方法,需要进行油井出水规律的识别才能针对性地进行分类堵水。
本发明实施例中,提供了一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备。该方案通过提供一种井封隔体注入控堵水模拟方法,实现对于水注入后的井周围裂缝分布规律的动态建立,以裂缝扩张和延伸模型为基础,综合分析实际充填井的成像测井数据、清水实验数据、充填施工数据、措施后生产数据,得到了靶区裂缝充填主控因素。
根据本发明实施例第一方面,提供一种裂缝型油藏封隔体注入控堵水模拟方法。
图1是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法的流程图。
在一个或多个实施例中,优选地,所述一种裂缝型油藏封隔体注入控堵水模拟方法包括:
S101、通过截取获得单井的截面图,设置单井模型;
S102、选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
S103、在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
S104、准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
S105、获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
S106、获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
在本发明实施例中,为了能够进行对实际井隔体注入控堵水的模拟,确定了具体的思路如下,对井周裂缝分布直接影响充填效果,首先建立并校准井周天然裂缝,结合建立井周天然裂缝时,以成像测井为约束随机生成,对长度分布、初始开度分布等作为待校准量,在充填前的清水注入测试数据对井周裂缝分布及流动能力的最直接反应,以此作为校准依据。
图2是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的通过截取获得单井的截面图,设置单井模型的流程图。
如图2所示,在一个或多个实施例中,优选地,所述通过截取获得单井的截面图,设置单井模型,具体包括:
S201、通过截取获得单井的截面图;
S202、根据所述单井的截面图以井边距离为基础设置单井模型;
S203、对所述单井模型进行水平方向的加密。
在本发明实施例中,优选地,截取单井模型,模型边界距井最近处200m为基础设置为单井模型,并对单井模型水平方向加密至15m,纵向不变。
图3是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片的流程图。
如图3所示,在一个或多个实施例中,优选地,所述选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片,具体包括:
S301、选取时间较长的清水测试数据,所述时间较长的清水测试数据具体为超过预设的时间长度的测试数据;
S302、以净注入量为控制条件,生成注水制度文件;
S303、根据管流模拟生成注入VFP表,其中,所述VFP表为关系型数据表;
S304、根据所述关系型数据表计算比较井口压力;
S305、充填模拟需要离散化的裂缝;
S306、通过成像测井得到的裂缝统计表,生成过井的裂缝片。
在本发明实施例中,选取A井第一次清水测试,保证较多数据点,并以净注入量为控制条件,生成注水制度文件,此外要根据管流模拟生成关系型数据表,意义在于将测试数据生成为计算机能够识别使用的对应数据表格式的一种标准化数据处理,在此基础上,充填模拟需要离散化的裂缝,不能使用双孔模型,双孔模型为模型包括基质和裂缝两种通道,基质用于储油,裂缝用于油气流通;若直接使用双孔模型作为离散裂缝,也过于复杂。因此通过成像测井得到的裂缝统计表,生成过井的裂缝片建立的过井裂缝片为待充填的初始模型,后续需对裂缝分布进行拟合。
图4是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度的流程图。
如图4所示,在一个或多个实施例中,优选地,所述在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度,具体包括:
S401、在所述过井的裂缝片区域之外,生成随机裂缝,形成全部的裂缝基础数据;
S402、根据所述裂缝基础数据,设置所有的所述裂缝的初始开度;
S403、拟合清水测试中的井口压力,使得压力的拟合值与井周流动能力一致的情况下,更新调整所述裂缝基础数据中的所述裂缝分布和所述裂缝的初始开度。
在本发明实施例中,生成随机裂缝时的纵向沟通性、密度、开度产生均作为拟合变量,可认为是真实裂缝系统的简化,但保证关键流动指标一致,在此基础上,调整裂缝分布,以及所有裂缝的初始开度,拟合清水测试井口压力。压力拟合可保证模型井周流动能力与真实情况一致。
图5是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力的流程图。
如图5所示,在一个或多个实施例中,优选地,所述准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力,具体包括:
S501、准备充填动态模拟数据,根据开度、导流能力、渗透率、颗粒浓度、应力的场数据变化开展充填动态模拟;
S502、在三维视窗中查看力饱和度;
S503、在二维曲线中查看充填压力和排量;
S504、通过调整A井井上裂缝开度渗透率、井周裂缝纵向联通程度和密度拟合清水实验压力曲线,判断所述裂缝分布反映井周流动能力。
在本发明实施例中,通过调整A井井上裂缝开度渗透率,及井周裂缝纵向联通程度、密度等,拟合清水实验压力曲线,反应当前裂缝分布井周流动能力。
图6是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度的流程图。
如图6所示,在一个或多个实施例中,优选地,所述获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度,具体包括:
S601、获取当前的模拟结果,判断充填是否主要发生在井上裂缝,且距离井更远的裂缝较少波及;
S602、判断井上裂缝开度随充填进行均逐渐增大,增大幅度从井眼处向周围逐渐减小;
S603、判断裂缝导流能力由于开度增大一开始普遍增大,随后颗粒进入封堵,从中心向外围逐渐变小;
S604、判断半年内含水率数据的模拟结果与实际动态趋势和幅度的吻合程度。
在本发明实施例中,充填模拟拟合后,充填井口压相对一致时,正确反映充填流动特征,对模拟中的整片裂缝进行充填程度,根据导流能力变化关系统计,基本满足前期概念模型下的理论规律,通过上述方式为进行填充,对比不充填的情况,初期含水率降低接近10个百分点,3个月后控水效果逐渐减弱。
图7是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟方法中的获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值的流程图。
如图7所示,在一个或多个实施例中,优选地,所述获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值,具体包括:
S701、获取所述裂缝开度大小和所述裂缝长度,利用第一计算公式计算裂缝开度变化率和裂缝长度变化率;
S702、利用第二计算公式计算所有的相邻的裂缝之间的关联度,不相邻的裂缝之间关联度为零;
S703、利用第三计算公式计算每个裂缝的总体关联度;
S704、利用第四计算公式计算裂缝开度大小的更新值;
S705、利用第五计算公式计算裂缝长度的更新值;
所述第一计算公式为:
其中,D 1_i_T 为裂缝开度变化率,D 2_i_T 为裂缝长度变化率,K 1_i_T 为裂缝开度大小,K 2_i_T 为裂缝长度;
所述第二计算公式为:
其中,g ij_T 为第i裂缝与第j裂缝的关联度,P为分辨率系数,
为第i裂缝与第j裂缝开度的二级最大差,为第i裂缝与
第j裂缝开度的二级最小差,为第i裂缝与第j裂缝长度
的二级最大差,为第i裂缝与第j裂缝长度的二级最小差;
所述第三计算公式为:
其中,G i 为第i裂缝的总体关联度,n为裂缝总数,g ii_T 为0;
所述第四计算公式为:
其中,TK 1_i_T 为裂缝开度大小的更新值;
所述第五计算公式为:
其中, TK 2_i_T 为裂缝长度的更新值,e为纳皮尔常数。
在本发明实施例中,根据实时的井周裂缝的纵向关联程度和横向关联程度进行在线的分析,实现井周的裂缝的开度进行根据不同裂缝的大小,而根据时间尺度的分级进行动态的调整和修正,使得井周裂缝模拟过程更加能够呈现裂缝之间的相互影响,并能够适应不同裂缝的大小的发展与扩展。
根据本发明实施例第二方面,提供一种裂缝型油藏封隔体注入控堵水模拟系统。
图8是本发明一个实施例的一种裂缝型油藏封隔体注入控堵水模拟系统的结构图。
在一个或多个实施例中,优选地,所述一种裂缝型油藏封隔体注入控堵水模拟系统包括:
单井模型生成模块801,用于通过截取获得单井的截面图,设置单井模型;
井周裂缝生成模块802,用于选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
裂缝拟合模块803,用于在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
充填动态模拟模块804,用于准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
指标输出模块805,用于获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
规则分析模块806,用于获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
本发明实施例中,提供了一种系统结构,这种系统结构下,各个模块分别进行模型生成、井周裂缝生成、裂缝拟合等操作,进而实现对于综合分析和数据扩展。
根据本发明实施例第三方面,提供一种计算机可读存储介质,其上存储计算机程序指令,所述计算机程序指令在被处理器执行时实现如本发明实施例第一方面中任一项所述的方法。
根据本发明实施例第四方面,提供一种电子设备。图9是本发明一个实施例中一种电子设备的结构图。图9所示的电子设备为通用裂缝型油藏封隔体注入控堵水模拟装置,其包括通用的计算机硬件结构,其至少包括处理器901和存储器902。处理器901和存储器902通过总线903连接。存储器902适于存储处理器901可执行的指令或程序。处理器901可以是独立的微处理器,也可以是一个或者多个微处理器集合。由此,处理器901通过执行存储器902所存储的指令,从而执行如上所述的本发明实施例的方法流程实现对于数据的处理和对于其它装置的控制。总线903将上述多个组件连接在一起,同时将上述组件连接到显示控制器904和显示装置以及输入/输出(I/O)装置905。输入/输出(I/O)装置905可以是鼠标、键盘、调制解调器、网络接口、触控输入装置、体感输入装置、打印机以及本领域公知的其他装置。典型地,输入/输出装置905通过输入/输出(I/O)控制器906与系统相连。
本发明的实施例提供的技术方案可以包括以下有益效果:
本发明方案中,建立了非结构化网格双重介质油藏模型及完成历史拟合,分析了靶区裂缝发育特征描述和分布规律。
本发明方案中,通过对裂缝充填动态过程进行模拟,实现了靶区裂缝充填控堵水机理及主控因素的分析。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (10)
1.一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,该方法包括:
通过截取获得单井的截面图,设置单井模型;
选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
2.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述通过截取获得单井的截面图,设置单井模型,具体包括:
通过截取获得单井的截面图;
根据所述单井的截面图以井边距离为基础设置单井模型;
对所述单井模型进行水平方向的加密。
3.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片,具体包括:
选取时间较长的清水测试数据,所述时间较长的清水测试数据具体为超过预设的时间长度的测试数据;
以净注入量为控制条件,生成注水制度文件;
根据管流模拟生成注入VFP表,其中,所述VFP表为关系型数据表;
根据所述关系型数据表计算比较井口压力;
充填模拟需要离散化的裂缝;
通过成像测井得到的裂缝统计表,生成过井的裂缝片。
4.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度,具体包括:
在所述过井的裂缝片区域之外,生成随机裂缝,形成全部的裂缝基础数据;
根据所述裂缝基础数据,设置所有的所述裂缝的初始开度;
拟合清水测试中的井口压力,使得压力的拟合值与井周流动能力一致的情况下,更新调整所述裂缝基础数据中的所述裂缝分布和所述裂缝的初始开度。
5.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力,具体包括:
准备充填动态模拟数据,根据开度、导流能力、渗透率、颗粒浓度和应力的场数据变化开展充填动态模拟;
在三维视窗中查看力饱和度;
在二维曲线中查看充填压力和排量;
通过调整A井井上裂缝开度渗透率、井周裂缝纵向联通程度和密度拟合清水实验压力曲线,判断所述裂缝分布反映井周流动能力。
6.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度,具体包括:
获取当前的模拟结果,判断充填是否主要发生在井上裂缝,且距离井更远的裂缝较少波及;
判断井上裂缝开度随充填进行均逐渐增大,增大幅度从井眼处向周围逐渐减小;
判断裂缝导流能力由于开度增大一开始普遍增大,随后颗粒进入封堵,从中心向外围逐渐变小;
判断半年内含水率数据的模拟结果与实际动态趋势和幅度的吻合程度。
7.如权利要求1所述的一种裂缝型油藏封隔体注入控堵水模拟方法,其特征在于,所述获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值,具体包括:
获取所述裂缝开度大小和所述裂缝长度,利用第一计算公式计算裂缝开度变化率和裂缝长度变化率;
利用第二计算公式计算所有的相邻的裂缝之间的关联度,不相邻的裂缝之间关联度为零;
利用第三计算公式计算每个裂缝的总体关联度;
利用第四计算公式计算裂缝开度大小的更新值;
利用第五计算公式计算裂缝长度的更新值;
所述第一计算公式为:
其中,D 1_i_T 为裂缝开度变化率,D 2_i_T 为裂缝长度变化率,K 1_i_T 为裂缝开度大小,K 2_i_T 为裂缝长度;
所述第二计算公式为:
其中,g ij_T 为第i裂缝与第j裂缝的关联度,P为分辨率系数,
为第i裂缝与第j裂缝开度的二级最大差,为第i裂缝与第j
裂缝开度的二级最小差,为第i裂缝与第j裂缝长度的二
级最大差,为第i裂缝与第j裂缝长度的二级最小差;
所述第三计算公式为:
其中,G i 为第i裂缝的总体关联度,n为裂缝总数,g ii_T 为0;
所述第四计算公式为:
其中,TK 1_i_T 为裂缝开度大小的更新值;
所述第五计算公式为:
其中, TK 2_i_T 为裂缝长度的更新值。
8.一种裂缝型油藏封隔体注入控堵水模拟系统,其特征在于,该系统用于实施如权利要求1-7中任一项所述的方法,该系统包括:
单井模型生成模块,用于通过截取获得单井的截面图,设置单井模型;
井周裂缝生成模块,用于选取时间较长的清水测试数据,生成注水制度文件,生成过井的裂缝片;
裂缝拟合模块,用于在过井的裂缝片区域之外,生成随机裂缝,更新裂缝分布和裂缝的初始开度;
充填动态模拟模块,用于准备充填动态模拟数据,判断所述裂缝分布反映井周流动能力;
指标输出模块,用于获取当前的模拟结果,判断含水率数据的模拟结果与实际动态趋势和幅度的吻合程度;
规则分析模块,用于获取裂缝开度大小和裂缝长度,计算裂缝开度大小的更新值和裂缝长度的更新值。
9.一种计算机可读存储介质,其上存储计算机程序指令,其特征在于,所述计算机程序指令在被处理器执行时实现如权利要求1-7中任一项所述的方法。
10.一种电子设备,包括存储器和处理器,其特征在于,所述存储器用于存储一条或多条计算机程序指令,其中,所述一条或多条计算机程序指令被所述处理器执行以实现如权利要求1-7任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211432710.2A CN115496014B (zh) | 2022-11-16 | 2022-11-16 | 一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211432710.2A CN115496014B (zh) | 2022-11-16 | 2022-11-16 | 一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115496014A true CN115496014A (zh) | 2022-12-20 |
CN115496014B CN115496014B (zh) | 2023-03-10 |
Family
ID=85115872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211432710.2A Active CN115496014B (zh) | 2022-11-16 | 2022-11-16 | 一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115496014B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN105606775A (zh) * | 2016-02-29 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种新型裂缝模拟评价实验方法 |
CN105756651A (zh) * | 2016-04-18 | 2016-07-13 | 中国海洋石油总公司 | 一种调整井正压钻井堵水与分段压裂增产工艺方法 |
CN110210144A (zh) * | 2019-06-05 | 2019-09-06 | 西南石油大学 | 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法 |
CN112012712A (zh) * | 2020-08-27 | 2020-12-01 | 西安石油大学 | 一种嵌入式离散裂缝的注水生长缝数值模拟方法和装置 |
CN114482913A (zh) * | 2020-10-27 | 2022-05-13 | 中国石油化工股份有限公司 | 裂缝性见水油井水平井封堵长度的确定方法、封堵方法 |
-
2022
- 2022-11-16 CN CN202211432710.2A patent/CN115496014B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102852516A (zh) * | 2012-04-19 | 2013-01-02 | 北京大学 | 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置 |
CN105606775A (zh) * | 2016-02-29 | 2016-05-25 | 中国石油天然气股份有限公司 | 一种新型裂缝模拟评价实验方法 |
CN105756651A (zh) * | 2016-04-18 | 2016-07-13 | 中国海洋石油总公司 | 一种调整井正压钻井堵水与分段压裂增产工艺方法 |
CN110210144A (zh) * | 2019-06-05 | 2019-09-06 | 西南石油大学 | 一种暂堵剂促进水平井压裂裂缝均匀扩展的优化设计方法 |
CN112012712A (zh) * | 2020-08-27 | 2020-12-01 | 西安石油大学 | 一种嵌入式离散裂缝的注水生长缝数值模拟方法和装置 |
CN114482913A (zh) * | 2020-10-27 | 2022-05-13 | 中国石油化工股份有限公司 | 裂缝性见水油井水平井封堵长度的确定方法、封堵方法 |
Non-Patent Citations (4)
Title |
---|
张轶栋等: "裂缝性底水油藏单井水锥数值模拟研究", 《石油化工应用》 * |
杜殿发等: "底水潜山油藏出水规律研究", 《石油钻采工艺》 * |
杨庆军等: "含天然闭合缝开启的低渗透油藏精细数值模拟", 《江汉石油学院学报》 * |
龙旭等: "塔河缝洞型油藏单井含水变化类型定量评价", 《西南石油大学学报(自然科学版)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN115496014B (zh) | 2023-03-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sheng et al. | A new fractal approach for describing induced-fracture porosity/permeability/compressibility in stimulated unconventional reservoirs | |
US11892592B2 (en) | Pressure transient testing with sensitivity analysis | |
US10909281B2 (en) | History matching of hydrocarbon production from heterogenous reservoirs | |
RU2605192C2 (ru) | Обновление микросейсмических гистограммных данных | |
US10294759B2 (en) | Method and system for hydraulic fracturing based on skin factor analysis | |
US10309216B2 (en) | Method of upscaling a discrete fracture network model | |
Elmo et al. | Guidelines for the quantitative description of discontinuities for use in discrete fracture network modelling | |
Yang et al. | Assessment of production interference level due to fracture hits using diagnostic charts | |
Dong et al. | Quantitative evaluation and classification method of the cataclastic texture rock mass based on the structural plane network simulation | |
CN115496014B (zh) | 一种裂缝型油藏封隔体注入控堵水模拟方法、系统及设备 | |
WO2012114146A2 (en) | Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region | |
Xue et al. | Estimating the permeability of fractured rocks using topological characteristics of fracture network | |
US20220245300A1 (en) | Method and system for constructing statistical emulators for reservoir simulation models with uncertain inputs | |
CN111155980B (zh) | 一种水流优势通道识别方法及装置 | |
CN114236624B (zh) | 基于电磁法估算压裂改造空间体积的方法和系统 | |
Liu et al. | Classification and evaluation for stimulated reservoir volume (SRV) estimation models using microseismic events based on three typical grid structures | |
Kluckner et al. | Estimation of the in situ block size in jointed rock masses using three-dimensional block simulations and discontinuity measurements | |
Morelli | Empirical assessment of the mean block volume of rock masses intersected by four joint sets | |
Palleske | A study of biases, assumptions and practical consideration for the use of discrete fracture networks in geomechanical practice | |
Torres et al. | A New Discrete Fracture Network Calibration Workflow using EDFM Method | |
US20180306939A1 (en) | Correcting Biases In Microseismic-Event Data | |
US11680465B2 (en) | Systems and methods for multiscale sector hydrocarbon reservoir simulation | |
Bellal et al. | Prediction of discrete and continuous wellbore logs values using artificial neural network, decision trees and support vectors machine algorithms | |
Karimi Sharif et al. | DFNAnalyzer: A web-based application for discrete fracture network analysis | |
Yan et al. | A Novel Workflow from StimPlan to EDFM for Complex Hydraulic Fracture Modeling and Production Simulation |
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 |