CN105528648B - 缝洞单元的生产动态预测方法及装置 - Google Patents
缝洞单元的生产动态预测方法及装置 Download PDFInfo
- Publication number
- CN105528648B CN105528648B CN201510855840.0A CN201510855840A CN105528648B CN 105528648 B CN105528648 B CN 105528648B CN 201510855840 A CN201510855840 A CN 201510855840A CN 105528648 B CN105528648 B CN 105528648B
- Authority
- CN
- China
- Prior art keywords
- fracture
- result information
- reservoir
- numerical
- flow behavior
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 139
- 238000000034 method Methods 0.000 title claims abstract description 100
- 238000004458 analytical method Methods 0.000 claims abstract description 104
- 238000012360 testing method Methods 0.000 claims abstract description 83
- 238000004088 simulation Methods 0.000 claims abstract description 30
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims abstract description 25
- 238000012512 characterization method Methods 0.000 claims description 94
- 230000007423 decrease Effects 0.000 claims description 65
- 238000011161 development Methods 0.000 claims description 46
- 239000000463 material Substances 0.000 claims description 17
- 238000011160 research Methods 0.000 claims description 13
- 125000005587 carbonate group Chemical group 0.000 claims description 8
- 235000013399 edible fruits Nutrition 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000000126 substance Substances 0.000 claims description 2
- 238000005303 weighing Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 208000010392 Bone Fractures Diseases 0.000 description 100
- 206010017076 Fracture Diseases 0.000 description 100
- 230000006399 behavior Effects 0.000 description 69
- 238000011156 evaluation Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 18
- 230000003247 decreasing effect Effects 0.000 description 10
- 230000008901 benefit Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000035699 permeability Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 210000002615 epidermis Anatomy 0.000 description 2
- 244000265913 Crataegus laevigata Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 244000153665 Ficus glomerata Species 0.000 description 1
- 235000012571 Ficus glomerata Nutrition 0.000 description 1
- 235000009047 Nauclea latifolia Nutrition 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 235000002639 sodium chloride Nutrition 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种缝洞单元的生产动态预测方法及装置,包括:确定碳酸盐岩油藏缝洞单元的第一储层动态描述结果信息;建立现代产量递减分析数值模型,确定缝洞单元的第二储层动态描述结果信息;建立数值试井模型,确定缝洞单元的第三储层动态描述结果信息;比较第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;若第一差值小于预设阈值,建立缝洞单元的数值模拟动态模型,确定第四储层动态描述结果信息;比较第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同预设阈值的大小,若第二差值小于预设阈值,基于现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测。
Description
技术领域
本发明涉及缝洞型碳酸盐岩油藏的动态描述及预测技术领域,特别涉及一种缝洞单元的生产动态预测方法及装置。
背景技术
国内缝洞型碳酸盐岩油藏孔、洞、缝普遍发育,储集空间具有多尺度性,非均质性极强,国外没有可借鉴的成熟开发经验及技术,缝洞型碳酸盐岩油藏的开发是亟待解决的世界级难题。其中,缝洞单元的综合动态描述是缝洞型油藏有效开发的基础。只有正确认识缝洞单元,才能制定合理的开发对策,正确预测缝洞单元的开发动态。目前,国内外对油藏的动态描述手段多采用试井分析、现代产量递减分析方法进行研究。其中,试井技术发展较早,研究文献较多,在解析求解与数值求解方面均有深入研究。现代产量递减分析法自上世纪八十年代起发展迅速,目前已经有Blasingame典型曲线法、Agarwal-Gardner典型曲线法、流动物质平衡法等多种解析求解方法,同时也形成了数值求解方法。目前现代产量递减分析法研究中多采用解析法,部分文献有采用现代产量递减分析数值方法进行描述。同时,国内外已有研究人员采用解析法试井及现代生产递减分析解析法相结合的方式进行储层动态描述,从而减小研究的不确定性、使描述结果更可靠。文献中较少提及采用数值试井和数值现代生产递减分析相结合的方式进行储层动态描述,以及如何实现解析法与数值法相结合。而对于国内复杂的缝洞型碳酸盐岩油藏来说,单纯采用试井、现代生产递减分析一种方法无法准确评价其储层参数,即使采用试井与现代生产递减分析相结合,也不能解决所有问题,尤其是对于存在边底水的缝洞单元来说。因此,直到目前为止,尚未有有效的缝洞型碳酸盐岩油藏缝洞单元综合动态描述及预测的方法,未见有综合物质平衡方法、解析试井及数值试井、解析现代生产递减分析及数值现代生产递减分析以及数值模拟 技术,综合进行缝洞型碳酸盐岩油藏缝洞单元动态描述及预测的报道。
发明内容
本发明实施例提供一种缝洞单元的生产动态预测方法及装置,以综合利用物质平衡法、解析试井及数值试井、现代生产递减分析解析法及现代生产递减分析数值法以及数值模拟技术进行储层动态描述,提升评价结果的准确性及可靠性。
为了实现上述目的,本发明实施例提供一种缝洞单元的生产动态预测方法,该生产动态预测方法包括:
步骤1:确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
步骤2:基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
步骤3:基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
步骤4:比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
步骤5:如果所述第一差值小于所述预设阈值,基于所述历史数据、所述第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
步骤6:比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测。
一实施例中,如果所述第一差值大于所述预设阈值,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致,然后执行所述步骤2及步骤3,依次迭代,直至所述第一差值小于所述预设阈值。
一实施例中,如果所述第二差值大于所述预设阈值,将所述现代产量递减分析数值模型的参数及所述数值试井模型的参数更新为与所述数值模拟动态模型的参数一致,然后执行所述步骤2至步骤5,依次迭代,直至所述第二差值小于所述预设阈值。
一实施例中,所述现代产量递减分析数值模型与所述数值试井模型采用相同的基础地质模型及网格模型。
一实施例中,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,包括:根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
为了实现上述目的,本发明实施例提供一种缝洞单元的生产动态预测装置,该生产动态预测装置包括:
储层动态描述初步确定单元,用于确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
现代产量递减分析数值模型创建单元,用于基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
数值试井模型创建单元,用于基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
比较单元,用于比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
数值模拟动态模型创建单元,用于基于所述历史数据、第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
生产动态预测预测单元,比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈 值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测。
一实施例中,所述的生产动态预测装置还包括:第一更新单元,用于当所述第一差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致。
一实施例中,所述的生产动态预测装置还包括:
第二更新单元,用于当所述第二差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数及所述数值试井模型的参数更新为与所述数值模拟动态模型的参数一致。
一实施例中,所述现代产量递减分析数值模型与所述数值试井模型采用相同的基础地质模型及网格模型。
一实施例中,所述储层动态描述初步确定单元具体用于:根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
本发明综合利用物质平衡法、解析试井及数值试井、现代生产递减分析解析法及现代生产递减分析数值法以及数值模拟技术进行储层动态描述,并基于动态描述结果基础上进行开发动态的预测,综合利用了各种动态监测资料及各方法的优点,提升了评价结果的准确性及可靠性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的生产动态预测方法流程图;
图2为国内某缝洞型碳酸盐岩油藏某生产井的物质平衡法评价结果示意图;
图3A至图3C为国内某缝洞型碳酸盐岩油藏某生产井的现代产量递减分析解析法曲线拟合结果示意图;
图4A、图4B及图4C为国内某缝洞型碳酸盐岩油藏某生产井的解析试井曲线拟合结果示意图;
图5A及图5B为国内某缝洞型碳酸盐岩油藏某区块缝洞单元划分结果及现代产量递减分析数值法初始模型示意图;
图6国内某缝洞型碳酸盐岩油藏某区块缝洞单元经综合动态描述后的缝洞单元边界示意图;
图7国内某缝洞型碳酸盐岩油藏某生产井的生产动态预测结果示意图;
图8为本发明实施例的生产动态预测装置结构框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种缝洞单元的生产动态预测方法,如图1所示,该生产动态预测方法包括:
S101:确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
S102:基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
S103:基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
S104:比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
S105:如果所述第一差值小于所述预设阈值,基于所述历史数据、所述第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
S106:比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测。
由图1所示的流程可知,本发明基于储层动态描述结果信息,综合利用物质平衡法、解析试井及数值试井、现代生产递减分析解析法及现代生产递减分析数值法以及数值模拟技术进行储层动态描述,并基于动态描述结果基础上进行开发动态的预测,综合利用了各种动态监测资料及各方法的优点,提升了评价结果的准确性及可靠性。
一实施例中,如图1所示,如果S104中比较得出第一差值大于所述预设阈值,进行S108,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致,然后循环执行所述S102及S103,直至第一差值小于预设阈值。
一实施例中,如图1所示,如果S106中比较得出第二差值大于所述预设阈值,进行S109,将现代产量递减分析数值模型的参数及数值试井模型的参数更新为与数值模拟动态模型的参数一致,然后循环执行S102至S105,直至第二差值小于预设阈值。
为了提高动态预测的准确性及一致性,具体实施时,使现代产量递减分析数值模型与数值试井模型采用相同的基础地质模型及网格模型,所有参数设置一致,采用数值试井技术完成缝洞单元内有压力恢复测试井的压恢数据的拟合。
S101中,在确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息时,可以根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
S106中,需要比较不同方法预测结果的差异,如有数值模拟结果,可采用数值模拟结果,如若没有,选择现代产量递减分析数值模型预测结果即可。
图1所示的流程介绍了本发明的具体实施步骤。在采用上述步骤进行动态描述之前,需要首先准备需要的数据,包括基本的地震、地质数据(包括缝洞雕刻结果、地震反演结果、缝洞单元边界、测井解释孔隙度、有效厚度等),缝洞单元内所有生产 井的生产动态数据(包括日产量、压力数据),所有井的动态监测数据(包括压力恢复测试资料、流压及静压测试资料)。其中,还需要将所有井的井口压力数据折算为井底流压,折算的方法包括Beggs&Brill方法、Gray方法、Petalas&Aziz方法、Hagedorn&Brown方法等。通过折算使得折算流压与测试流压一致,从而在现代产量递减分析解析法及数值法动态描述中,采用折算流压进行分析。因此,要完整按照本发明进行缝洞单元综合动态描述及预测,要求研究缝洞单元内部至少有1口井进行过压力恢复测试(进行试井解释),至少1口井有2个以上静压测试点(进行物质平衡分析),至少有1口生产井进行了生产(进行现代产量递减分析)。如果确实缺少相应数据,其相应步骤可省略,例如若没有生产井具有2个以上的静压测试数据点,则S101中的物质平衡分析法可不做。
本发明之所以综合了物质平衡法、解析试井及数值试井、现代生产递减分析解析法及现代生产递减分析数值法以及数值模拟技术进行储层动态描述,是因为这些方法均有相似的特征,均具有储层动态描述的功能,每种方法各有所长且不同方法有许多共性,但同时每种方法又有各自的优缺点。各种方法在储层参数评价的过程中可以有机结合、互相约束,从而保证解释结果的可靠性及一致性。下面重点介绍一下上述每种方法的优缺点,以便于对本发明实施步骤的理解。
不同方法所采用的基础输入数据不同。物质平衡法是基于测试的静压数据及对应时间的累积产量数据进行分析;现代生产递减分析及数值模拟方法主要采用每日测量的产量、压力数据进行分析;试井分析需要采用专门的高精度的压力不稳定测试数据进行分析。而各种方法均可进行储层参数的评价,如渗透率、动态储量等;均可进行生产动态的指标预测,如定压生产预测产量或定产量生产预测生产压力。但同时各种方法又具有各自的特点,现代生产递减分析、物质平衡方法评价动态储量较可靠,试井分析对储层渗透率、井表皮系数、缝洞单元边界响应等评价结果更可靠。各种方法的评价手段不同,现代生产递减分析法与试井分析法相似,均采用典型曲线拟合方法进行储层参数的评价,均有解析法和数值法两种不同分析方法。其中,解析法求解速度快、动态描述快,但由于只能考虑简单边界、单相流情况且不能考虑井间干扰等问题,对于复杂的碳酸盐岩油藏来说解析法求解结果可作为参考结果。而两种方法均可以通过建立模型采用数值求解方法来考虑复杂边界、多相流、多井干扰等情况,但是数值解法求解速度慢、动态描述慢。正因为如此,本发明的具体实施步骤中首先采用 解析法进行分析,数值法初始参数采用解析法,这样数值法采用解析法的结果作为初始结果后可以大大减少数值法的分析时间,大大提高了工作效率。另外,试井解释的基础数据为高精度压力计测试的大量压力恢复数据,现代生产递减分析的基础数据为生产井每天的日产量、压力数据,一般来说,现代产量递减分析的数据点或数据量要明显小于试井分析的数据量。因此,对于一口井的分析来说,现代产量递减分析所用时间一般要少于试井分析所用的时间。同时,现代生产递减分析评价动态储量可靠,可以落实缝洞单元的规模,从而确定缝洞单元的面积或大小,有助于采用数值法进行缝洞单元边界的划分,而试井分析常常受限于测试时间长短的影响,测试数据不一定能完全反映缝洞单元边界的响应,但试井分析曲线受边界的响应更加明显。因此,本发明S102中首先采用现代产量递减数值分析法落实缝洞单元动态储量及缝洞单元边界,再由S103采用数值试井进一步校正缝洞单元边界的解释。这样也是会大大提高综合动态描述的效率,大大减少模型反复修改的次数。而数值模拟法耗时最长,因此数值模拟法是在所有方法分析完成后,再基于这些方法的评价结果进行动态的描述及预测。
下面具体以实施例来介绍本发明中各种方法,以辅助对本发明的理解。
图2所示为采用物质平衡法评价的国内某缝洞型碳酸盐岩油藏某生产井的压力拟合结果及动态储量、水体大小评价结果示意图。物质平衡法需要生产井具有2个以上的静压测试点才可进行分析,且若想评价结果准确,还需要生产井具有较长的生产历史,一般采出程度达到20%以上计算结果较准确。物质平衡方法原理简单,计算方便快捷,主要用于评价单井或缝洞单元的动态储量,尤其对于含边低水的缝洞单元的动态储量及水体大小评价来说,该方法具有较强的优势。
如图3A至图3C所示为采用现代产量递减分析解析法对国内某缝洞型碳酸盐岩油藏某生产井的动态数据拟合结果示意图。由图3A至图3C可以看出,该井生产数据拟合较好。同时,现代产量递减分析解析法分析快、耗时短,可结合物质平衡法评价结果,准确落实生产井动态储量及缝洞单元规模,从而有利于数值法对缝洞边界的预判。
如图4A、图4B及图4C所示为采用解析试井法对国内某缝洞型碳酸盐岩油藏某生产井的实测压力恢复数据的拟合结果示意图。由图4A、图4B及图4C可以看出,该井生产数据拟合较好。同时,解析试井同现代产量递减分析解析法分析法类似,也 具有分析快、耗时短的特点,同时试井分析除了可准确评价储层渗透率、井表皮系数外,还可进一步落实缝洞单元边界距离井的距离,如评价该井三边具有封闭边界的响应,三个边距离井的距离分别为517米、361米、350米。
基于物质平衡分析、现代产量递减分析解析法及解析试井法评价结果,并综合地震解释与反演、地质研究对缝洞单元的认识,建立现代产量递减分析数值法求解动态模型,如图5A所示为地震及地质研究对缝洞单元边界划分结果示意图,图5B所示为建立的缝洞单元初始数值模型示意图。
按照本发明步骤2、步骤3、步骤4、步骤5、步骤6进行分析后,实现缝洞单元内所有井的生产动态数据的拟合及压力恢复资料的拟合,拟合情况如如图3A至图3C、图4A、图4B及图4C类似,这里不再赘述。通过综合分析后,可最终准确确定缝洞单元内所有单井的动态储量、渗透率、井表皮、水体大小等,同时准确确定了缝洞单元的边界。如图6所示为最终确定的国内某缝洞型碳酸盐岩油藏某区块的各个缝洞单元边界情况。对比图5A及图5B、图6,可明显发现两者之间的差异,对部分缝洞储集体的边界进行了修正。
基于综合动态描述的评价结果及最终建立的动态模型,可以进行缝洞单元内生产井的开发动态预测。如图7所示,为该缝洞单元内某口生产井的动态预测结果。井的实际生产动态证明本发明预测结果的准确性。
本发明综合了多种方法进行缝洞单元的动态描述,综合利用了各种方法的优势、避免了各种方法的劣势,各种方法应用步骤具有较强的逻辑性及关联性。通过采用该发明进行缝洞单元动态描述,各种方法相互结合、互相约束,描述结果准确性高,建立的动态模型可靠,预测结果准确,减小了油藏开发的风险,提高了油藏开发效益。本发明技术在碳酸盐岩气藏缝洞单元动态描述中也适用。
本发明实施例提供一种缝洞单元的生产动态预测装置,如图8所示,该生产动态预测装置包括:储层动态描述初步确定单元801,现代产量递减分析数值模型创建单元802,数值试井模型创建单元803,比较单元804,数值模拟动态模型创建单元805及生产动态预测预测单元806。
储层动态描述初步确定单元801用于确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
现代产量递减分析数值模型创建单元802用于基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
数值试井模型创建单元803用于基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
比较单元804用于比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
数值模拟动态模型创建单元805用于基于所述历史数据、第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
生产动态预测预测单元806比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测。
一实施例中,所述的生产动态预测装置还包括:第一更新单元,用于当所述第一差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致。
一实施例中,所述的生产动态预测装置还包括:
第二更新单元,用于当所述第二差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数及所述数值试井模型的参数更新为与所述数值模拟动态模型的参数一致。
一实施例中,所述现代产量递减分析数值模型与所述数值试井模型采用相同的基础地质模型及网格模型。
一实施例中,所述储层动态描述初步确定单元具体用于:根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
本发明综合利用了各种动态监测资料及各方法的优点,准确进行储层动态描述,并基于动态描述结果基础上进行开发动态的预测,大大提升了评价结果的准确性、可靠性。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种缝洞单元的生产动态预测方法,其特征在于,包括:
步骤1:确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
步骤2:基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
步骤3:基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
步骤4:比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
步骤5:如果所述第一差值小于所述预设阈值,基于所述历史数据、所述第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
步骤6:比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测;
其中,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,包括:根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
2.根据权利要求1所述的缝洞单元的生产动态预测方法,其特征在于,如果所述第一差值大于所述预设阈值,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致,然后循环执行所述步骤2及步骤3,直至所述第一差值小于所述预设阈值。
3.根据权利要求1所述的缝洞单元的生产动态预测方法,其特征在于,如果所述第二差值大于所述预设阈值,将所述现代产量递减分析数值模型的参数及所述数值试井模型的参数更新为与所述数值模拟动态模型的参数一致,然后循环执行所述步骤2至步骤5,直至所述第二差值小于所述预设阈值。
4.根据权利要求1至3中任一项所述的缝洞单元的生产动态预测方法,其特征在于,所述现代产量递减分析数值模型与所述数值试井模型采用相同的基础地质模型及网格模型。
5.一种缝洞单元的生产动态预测装置,其特征在于,包括:
储层动态描述初步确定单元,用于确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息,所述第一储层动态描述结果信息包括初始压力、储层物性、动态储量及水体大小;
现代产量递减分析数值模型创建单元,用于基于地震解释结果信息、地震反演结果信息、测井解释及地质研究认识结果信息、缝洞单元的生产历史数据及所述第一储层动态描述结果信息建立现代产量递减分析数值模型,确定所述缝洞单元的第二储层动态描述结果信息,所述生产历史数据包括所述缝洞单元的产量、压力信息,第二储层动态描述结果信息包括初始压力、储层物性、动态储量、水体大小及边界;
数值试井模型创建单元,用于基于所述第二储层动态描述结果信息建立数值试井模型,确定所述缝洞单元的第三储层动态描述结果信息;
比较单元,用于比较所述第二储层动态描述结果信息与第三储层动态描述结果信息的第一差值同预设阈值的大小;
数值模拟动态模型创建单元,用于基于所述历史数据、第二储层动态描述结果信息及第三储层动态描述结果信息建立缝洞单元的数值模拟动态模型,对所述缝洞单元内所有生产井生产历史进行拟合,确定所述缝洞单元的第四储层动态描述结果信息;
生产动态预测预测单元,比较所述第三储层动态描述结果信息与第四储层动态描述结果信息的第二差值同所述预设阈值的大小,如果所述第二差值小于所述预设阈值,基于所述现代产量递减分析数值模型、数值试井模型及缝洞单元的数值模拟动态模型进行缝洞单元的生产动态预测;
其中,所述储层动态描述初步确定单元具体用于:根据物质平衡方法、解析试井法、现代产量递减分析解析法对缝洞单元内各生产井进行综合描述,确定碳酸盐岩油藏的缝洞单元的第一储层动态描述结果信息。
6.根据权利要求5所述的缝洞单元的生产动态预测装置,其特征在于,还包括:
第一更新单元,用于当所述第一差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数更新为与所述数值试井模型的参数一致。
7.根据权利要求5所述的缝洞单元的生产动态预测装置,其特征在于,还包括:
第二更新单元,用于当所述第二差值大于所述预设阈值时,将所述现代产量递减分析数值模型的参数及所述数值试井模型的参数更新为与所述数值模拟动态模型的参数一致。
8.根据权利要求5至7中任一项所述的缝洞单元的生产动态预测装置,其特征在于,所述现代产量递减分析数值模型与所述数值试井模型采用相同的基础地质模型及网格模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510855840.0A CN105528648B (zh) | 2015-11-30 | 2015-11-30 | 缝洞单元的生产动态预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510855840.0A CN105528648B (zh) | 2015-11-30 | 2015-11-30 | 缝洞单元的生产动态预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105528648A CN105528648A (zh) | 2016-04-27 |
CN105528648B true CN105528648B (zh) | 2019-07-09 |
Family
ID=55770862
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510855840.0A Active CN105528648B (zh) | 2015-11-30 | 2015-11-30 | 缝洞单元的生产动态预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105528648B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106351651B (zh) * | 2016-08-26 | 2019-10-11 | 中国石油天然气股份有限公司 | 页岩气井产能的预测方法及装置 |
CN109138974A (zh) * | 2017-06-19 | 2019-01-04 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏离散数值试井分析方法及系统 |
CN109424362B (zh) * | 2017-08-29 | 2021-10-29 | 中国石油化工股份有限公司 | 计算底水油藏单井控制原油储量的方法及系统 |
CN108798634B (zh) * | 2018-07-09 | 2022-04-29 | 合肥辰工科技有限公司 | 一种缝洞型储层单洞-单层均值试井解释方法 |
CN111310338B (zh) * | 2020-02-19 | 2021-03-02 | 西南石油大学 | 一种原始地层压力修正的变储量系统典型曲线产生方法 |
CN111441767B (zh) * | 2020-05-11 | 2022-05-20 | 中国石油大学(华东) | 油藏生产动态预测方法及装置 |
CN113236207B (zh) * | 2021-07-13 | 2021-09-10 | 西南石油大学 | 一种强非均质性储层中产水气井的定产量递减预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455438A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型储层体积预测方法 |
CN103324770A (zh) * | 2012-03-21 | 2013-09-25 | 中国石油天然气股份有限公司 | 油藏的有效天然裂缝平面分布规律的预测方法 |
CN104153769A (zh) * | 2014-07-04 | 2014-11-19 | 中国石油大学(北京) | 一种缝洞型油藏流动单元的划分及评价方法 |
CN104750896A (zh) * | 2013-12-31 | 2015-07-01 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏数值模拟方法 |
-
2015
- 2015-11-30 CN CN201510855840.0A patent/CN105528648B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455438A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型储层体积预测方法 |
CN103324770A (zh) * | 2012-03-21 | 2013-09-25 | 中国石油天然气股份有限公司 | 油藏的有效天然裂缝平面分布规律的预测方法 |
CN104750896A (zh) * | 2013-12-31 | 2015-07-01 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏数值模拟方法 |
CN104153769A (zh) * | 2014-07-04 | 2014-11-19 | 中国石油大学(北京) | 一种缝洞型油藏流动单元的划分及评价方法 |
Non-Patent Citations (1)
Title |
---|
缝洞型油气藏储层类型对制定合理工作制度的影响;张学磊 等;《地质科技情报》;20110131;第30卷(第1期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105528648A (zh) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105528648B (zh) | 缝洞单元的生产动态预测方法及装置 | |
AU2017382546B2 (en) | Subsurface modeler workflow and tool | |
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN109740800B (zh) | 适用于隧道tbm掘进岩爆风险分级及预测方法与系统 | |
CN107402176B (zh) | 裂缝孔隙度的确定方法和装置 | |
CN104533400B (zh) | 一种重构测井曲线的方法 | |
CN109655903B (zh) | 页岩层横波速度预测方法及系统 | |
US20220206184A1 (en) | Integrated rock mechanics laboratory for predicting stress-strain behavior | |
CN107808032B (zh) | 页岩气储层含水饱和度计算方法及系统 | |
CN104316958B (zh) | 一种识别不同尺度地层断裂的相干处理方法 | |
CN107194180B (zh) | 油藏水驱体积波及系数的确定方法、装置及系统 | |
CN110414723A (zh) | 基于微震事件的裂缝油气藏历史拟合的方法、装置及系统 | |
CN106489068B (zh) | 测量值分析装置以及测量值分析方法 | |
CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
CN106597542B (zh) | 一种储层特征参数的预测方法及装置 | |
CN108131127B (zh) | 一种泡沫油型超重油油田的生产气油比的获取方法及装置 | |
CN111781662B (zh) | 一种储层裂缝参数获取方法、装置及设备 | |
CN106338768B (zh) | 一种生成储层预测属性数据的处理方法、装置及系统 | |
Piedipalumbo et al. | Updated f (T) gravity constraints from high-redshift cosmography | |
Dragovich et al. | FDE index for goodness‐of‐fit between measured and calculated response signals | |
CN111815769B (zh) | 逆冲推覆构造带构造的建模方法、计算设备及存储介质 | |
CN113138412B (zh) | 深层页岩气孔隙度地震预测方法及装置 | |
CN108415074A (zh) | 一种缝洞储集体地震响应特征分析方法及装置 | |
Kabir et al. | Understanding variable well performance in a chalk reservoir | |
CN114580048A (zh) | 隧道围岩压力拱计算方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |