CN112147051B - 一种基于渗透率分布形态的测压流度标准化方法 - Google Patents
一种基于渗透率分布形态的测压流度标准化方法 Download PDFInfo
- Publication number
- CN112147051B CN112147051B CN201911355166.4A CN201911355166A CN112147051B CN 112147051 B CN112147051 B CN 112147051B CN 201911355166 A CN201911355166 A CN 201911355166A CN 112147051 B CN112147051 B CN 112147051B
- Authority
- CN
- China
- Prior art keywords
- permeability
- data
- well
- fluidity
- standard
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
本发明涉及一种基于渗透率分布形态的测压流度标准化方法,包括:1)选择某种地层测试仪器在某地区的测试资料作为标准,建立不同井、不同深度的流度及其对应渗透率数据体;2)建立拟标准化井测试流度与储层渗透率对应关系数据,并制作渗透率数据的频率分布图;3)按照与拟标准化井相同的渗透率区间,并制作每口井的渗透率数据频率分布图,分别计算拟标准化井的渗透率分布形态与每口标准井渗透率分布形态的相似性系数;4)找到拟标准化井数据与标准测试结果的渗透率分布形态相似的井或子数据体,对拟标准化井数据进行标准化处理。
Description
技术领域
本发明涉及石油勘探开发领域,具体的说,涉及一种基于渗透率分布形态的测压流度标准化方法。
背景技术
储层流体流度信息大都依靠电缆地层测试方法来获取。电缆地层测试技术已日益成为石油勘探开发中地质家和工程师们所重视的一种测试手段。目前使用的电缆地层测试器包括模块式电缆地层动态测试器(MDT)、多次地层测式器(FMT)、油藏特性测试器(RCI)等,它们具有连续测压、实时光学流体性质分析以及常规取样和PVT取样等功能,并可根据测试资料计算地层渗透率和目的层流体密度等物理参数,从而能直接、快速、准确判别地层流体性质和帮助进行油藏综合评价。
由于不同公司仪器设备的原理和测试工艺存在差异,井眼和地层条件的影响程度不同,测试结果反演参数的算法不同,都会使相同地层和流体条件下获取的最终流度数据存在一定的偏差,这些差异会影响油藏产能的综合评价。因此,要统一不同仪器测得的数据就需要以一种测试仪器测试结果为标准对其它类型仪器进行标准化处理。
发明内容
针对上述问题,本发明的目的是提供一种基于渗透率分布形态的测压流度标准化方法,该方法能够基于常规岩心实验测试渗透率或测井计算渗透率资料,以一种地层测试仪器测试结果为标准对其他类型测试仪器获取的流度进行标准化处理。
为实现上述目的,本发明采取以下技术方案:一种基于渗透率分布形态的测压流度标准化方法,包括以下内容:
步骤1):建立标准地层测试流度及其对应渗透率数据体。选择某种地层测试仪器在某地区的测试资料作为标准,建立不同井、不同深度的流度(mob)及其对应渗透率(K)数据体(HBwi,MobBwi,KBwi),其中H表示数据点深度,下标B表示标准地层测试仪器类型,下标w代表井序号,下标i代表数据点序号。数据体中的流度Mob由地层测试获得,且为质量合格的数据点,K为与流度深度一致的储层渗透率值,它可以是岩心实验测试分析的结果,也可以是由测井资料计算的结果。
步骤2):建立拟标准化井储层渗透率分布形态。
(1)建立拟标准化井测试流度与储层渗透率对应关系数据(HSrj,MobSrj,KSrj),其中H表示数据点深度,下标S表示不同于标准地层测试仪器类型,下标r代表拟标准化井序号,下标j代表数据点序号,流度Mob由地层测试获得,且为质量合格的数据点,K为与流度深度一致的储层渗透率值,它可以是岩心实验测试分析的结果,也可以是由测井资料计算的结果。
(2)把拟标准化井的储层渗透率数据划分为n个小区间,分别统计在每个渗透率区间内的渗透率个数PSrk,k=1,…,n,并制作渗透率数据的频率分布图。
步骤3):建立拟标准化井与标准数据体渗透率分布形态的相似性关系。
(1)按照与拟标准化井相同的渗透率区间,分别统计标准数据体中每口井在每个渗透率区间内的渗透率个数PBwk,k=1,…,n,并制作每口井的渗透率数据频率分布图。
(2)分别计算拟标准化井的渗透率分布形态与每口标准井渗透率分布形态的最大相似性系数γw,w=1,…,mB。mB为标准井个数。
在mB个相似性系数中寻找最大值,且该值大于0.75,则认为该标准井渗透率数据的频数分布图与拟标准化井渗透率数据频率分布图相似。
(3)如果(2)中不存在和拟标准化井渗透率频率分布相似的井,此时可采用多井组合或由多井组合数据通过抽取选择的方式,构造一个子数据体,使这个组合或子数据体的渗透率数据频数分布图与拟标准化井的渗透率数据频数分布图形态相似,则该组合数据或子数据体即可作为虚拟的标准井数据。
步骤4):地层测试流度的标准化。假设通过分析,找到拟标准化E井数据(HSEj,MobSEj,KSEj)与标准井F测试结果(HBFi,MobBFi,KBFi)的渗透率分布形态相似。对这两个渗透率分布相似的井的流度数据进行对比,选取两口井的流度中值作为代表,由于中位数是所有数据的代表值,不受最大值或最小值的影响,则标准化后流度
Mob* SEj=MobSEj×BF/BE,j=1,…,mr
式中Mob* SEj是标准化后的流度,BF是标准井(或数据体)的流度数据MobBFi中值,BE是拟标准化井的流度数据MobSEj中值,mr是拟标准化井数据点个数。
本发明的有益效果在于:本发明基于地层测试仪器探测深度浅,在一个工区泥浆滤液的粘度差别不大,所以有相似渗透率分布井的流度分布也应该相似,通过分析储层渗透率分布的相似性,对不同地层测试仪器测量的差异进行标准化处理,对储层产能的动态分析具有重要的意义,有着非常重要的应用价值。本发明基于岩心或测井资料计算渗透率,数据容易获得,适用范围广、精确度高,可进行推广使用,为不同测试仪器测量的储层流度进行区域性的定量分析提供新的有效手段。避免了由于不同公司仪器设备的原理和测试工艺存的差异、井眼和地层条件的影响程度差异、测试结果反演参数的算法差异而造成的最终流度数据的偏差,减少了对油藏产能的综合评价的影响。
附图说明
图1为本发明实施例拟标准化井数据频数分布直方图,其中,(a)为拟标准化井渗透率频数分布直方图,(b)为拟标准化井流度频数分布直方图;
图2为本发明实施例标准井数据频数分布直方图,其中,(a)为标准井渗透率频数分布直方图,(b)为标准井流度频数分布直方图;
图3为本发明拟标准化井标准化后流度频数分布直方图。
具体实施方式
以下结合附图来对本发明进行详细的描绘。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
实施例1
通过本发明提供的基于渗透率分布形态相似的地层测试流度资料的标准化方法进行标准化流度的测算,包括以下步骤:
1、建立标准地层测试流度及其对应渗透率数据体。选择某种地层测试仪器(例如MDT)在某地区的测试资料作为标准,针对工区该测试仪器测得的每一个合格流度数据MobBwi,获取对应井某深度点HBwi的岩心实验测试分析渗透率或测井解释渗透率值KBwi,按不同井分别建立地层测试流度及其对应渗透率的数据体(HBwi,MobBwi,KBwi),其中H表示数据点深度,下标B表示标准地层测试仪器类型,下标w代表井序号,下标i代表数据点序号。数据体中的流度Mob由地层测试获得,且为质量合格的数据点,K为与流度深度一致的储层渗透率值,这一数据的获得可以是岩心实验测试分析的结果,也可以是由测井资料计算的结果。
2、建立拟标准化井储层渗透率分布形态。
(1)建立拟标准化井测试流度与储层渗透率对应关系数据(HSrj,MobSrj,KSrj),其中H表示数据点深度,下标S表示不同于标准地层测试仪器类型(例如RCI),下标r代表拟标准化井序号,下标j代表数据点序号,流度Mob由地层测试获得,且为质量合格的数据点,K为与流度深度一致的储层渗透率值,这一数据的获得可以是岩心实验测试分析的结果,也可以是由测井资料计算的结果。
(2)把拟标准化井的储层渗透率数据划分为n个小区间,分别统计在每个渗透率区间内的渗透率个数PSrk,k=1,…,n,并制作渗透率数据的频率分布图(图1(a))。图1(a)共划分出11个渗透率区间,分别是0~100、100~200、200~300、300~400、400~500、500~600、600~700、700~800、800~900、900~1000和>1000。
3、建立拟标准化井与标准数据体渗透率分布形态的相似性关系。
(1)按照与拟标准化井相同的渗透率区间,分别统计标准数据体中每口井在每个渗透率区间内的渗透率个数PBwk,k=1,…,n,并制作每口井的渗透率数据频率分布图,图2(a)是其中一口井的渗透率数据频率分布图实例,同样划分11个渗透率区间,分别是0~100、100~200、200~300、300~400、400~500、500~600、600~700、700~800、800~900、900~1000和>1000。
(2)分别计算拟标准化井的渗透率分布形态与每口标准井渗透率分布形态的最大相似性系数γw,w=1,…,mB。
在mb个相似性系数中寻找最大的值,且该值大于0.75,则认为该标准井渗透率数据的频数分布图与拟标准化井渗透率数据频率分布图相似,图1(a)和图2(a)具有明显的渗透率分布的相似性,相似性系数为0.87。
4、地层测试流度的标准化。通过分析,找到拟标准化E井数据(HSEj,MobSEj,KSEj)与标准井F测试结果(HBFi,MobBFi,KBFi)的渗透率分布形态相似(图1a与图2a)。对这两个渗透率分布相似的井的流度数据进行对比(如图1(b)和图2(b)所示),发现二者的分布相似,但其值有一定的差别。选取两口井的流度中值作为代表,由于中位数是所有数据的代表值,不受最大值或最小值的影响,图1(b)流度最小值2.3,最大值2909.4,平均值400.1,中值98.0;图2(b)流度最小值1.3,最大值4141.9,平均值580.0,中值173.7;标准化系数1.77;则标准化后流度
Mob* SEi=MobSEj×1.77,j=1,…,mr
式中Mob* SEj是标准化后的流度,MobSEj是标准化前的流度,mr是拟标准化井数据点个数。
图3是对图1(b)标准化后流度分布数据,标准化后的结果在渗透率分布相似的情况下,不同仪器的测试结果具有一致性。从图3中的结果可以明显看出,基于渗透率分布形态相似的地层测试流度资料的标准化方法是有效的,而且易于实现,是可以满足生产需求。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (5)
1.一种基于渗透率分布形态的测压流度标准化方法,其特征在于,包括:
步骤1):建立标准地层测试流度及其对应渗透率数据体;
步骤2):建立拟标准化井储层渗透率分布形态;
步骤3):建立拟标准化井与标准数据体渗透率分布形态的相似性关系;
步骤4):拟标准化井的测试流度的标准化计算;所述步骤1)包括:
1.1)标准的确立:选定地层测试仪器在某地区的测试资料作为标准地层测试仪器;
1.2)采集待测区域内不同井、不同深度的流度及其对应的渗透率,建立不同井、不同深度的流度及其对应渗透率数据体(HBwi,MobBwi,KBwi),其中H表示数据点深度,Mob表示流度,K表示渗透率,下标B表示标准地层测试仪器类型,下标w代表井序号,下标i代表数据点序号;
所述步骤2)包括:
2.1)建立拟标准化井测试流度与储层渗透率对应关系数据(HSrj,MobSrj,KSrj),其中H表示数据点深度,下标S表示不同于标准地层测试仪器类型,下标r代表拟标准化井号,下标j代表数据点序号;
2.2)把渗透率数据划分为n个小区间,分别统计在每个渗透率区间内的渗透率个数PSrk,k=1,…,n,并制作渗透率数据的频率分布图;
所述步骤4)具体包括,对这两个渗透率分布相似的井的流度数据进行对比,选取两口井的流度中值作为代表,由于中位数是所有数据的代表值,不受最大值或最小值的影响,则标准化后流度
Mob* SEj=MobSEj×BF/BE,j=1,…,mr
式中Mob* SEj是标准化后的流度,BF是标准井的流度数据MobBFi中值,BE是拟标准化井的流度数据MobSEj中值,mr是拟标准化井数据点个数。
2.根据权利要求1所述的基于渗透率分布形态的测压流度标准化方法,其特征在于,所述数据体中的流度MobBwi通过标准地层测试仪器采集,渗透率KBwi为与流度MobBwi深度一致的储层的渗透率值,所述渗透率KBwi通过岩心实验测试分或者测井资料计算获取。
3.根据权利要求2所述的基于渗透率分布形态的测压流度标准化方法,其特征在于,所述流度MobSrj由地层测试获得,渗透率KSrj为与流度MobSrj深度一致的储层的渗透率值,所述渗透率KSrj通过岩心实验测试分或者测井资料计算获取。
5.根据权利要求4所述的基于渗透率分布形态的测压流度标准化方法,其特征在于,如果在mB个相似性系数中寻找最大的值,且该值大于0.75,则认为该标准井渗透率数据的频率分布图与拟标准化井渗透率数据频率分布图相似,以步骤1中的数据体作为标准井数据进行步骤4)的计算;如果步骤3.2)中不存在和拟标准化井渗透率分布相似的井,此时采用多标准井组合数据通过抽取选择方式构造的一个子数据体,当所述多标准井组合数据或者构造的子数据体的渗透率数据频率分布图与拟标准化井的渗透率数据频率分布图形态相似,则所述多标准井组合数据或子数据体即可作为虚拟的标准井数据,进行步骤4)的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911355166.4A CN112147051B (zh) | 2019-12-25 | 2019-12-25 | 一种基于渗透率分布形态的测压流度标准化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911355166.4A CN112147051B (zh) | 2019-12-25 | 2019-12-25 | 一种基于渗透率分布形态的测压流度标准化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147051A CN112147051A (zh) | 2020-12-29 |
CN112147051B true CN112147051B (zh) | 2023-04-25 |
Family
ID=73892124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911355166.4A Active CN112147051B (zh) | 2019-12-25 | 2019-12-25 | 一种基于渗透率分布形态的测压流度标准化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147051B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN102619503A (zh) * | 2012-04-11 | 2012-08-01 | 上海大学 | 电缆地层测试器 |
CN104379870A (zh) * | 2012-02-13 | 2015-02-25 | 普拉德研究及开发股份有限公司 | 用于确定储层的径向流响应的渗透率或迁移率的方法 |
CN109710891A (zh) * | 2018-12-24 | 2019-05-03 | 核工业北京地质研究院 | 一种基于流动单元分类计算砂岩型铀矿砂体渗透率的方法 |
CN110043254A (zh) * | 2019-05-11 | 2019-07-23 | 中国海洋石油集团有限公司 | 一种基于电缆地层测试资料地层有效渗透率的获取方法 |
CN110162851A (zh) * | 2019-05-09 | 2019-08-23 | 中国海洋石油集团有限公司 | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 |
-
2019
- 2019-12-25 CN CN201911355166.4A patent/CN112147051B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN104379870A (zh) * | 2012-02-13 | 2015-02-25 | 普拉德研究及开发股份有限公司 | 用于确定储层的径向流响应的渗透率或迁移率的方法 |
CN102619503A (zh) * | 2012-04-11 | 2012-08-01 | 上海大学 | 电缆地层测试器 |
CN109710891A (zh) * | 2018-12-24 | 2019-05-03 | 核工业北京地质研究院 | 一种基于流动单元分类计算砂岩型铀矿砂体渗透率的方法 |
CN110162851A (zh) * | 2019-05-09 | 2019-08-23 | 中国海洋石油集团有限公司 | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 |
CN110043254A (zh) * | 2019-05-11 | 2019-07-23 | 中国海洋石油集团有限公司 | 一种基于电缆地层测试资料地层有效渗透率的获取方法 |
Non-Patent Citations (2)
Title |
---|
低渗透率储层流度计算改进方法探讨;张国栋等;《测井技术》;20160228;第40卷(第01期);全文 * |
电缆地层测试技术的发展及其在地层和油藏评价中的角色演变;孙华峰等;《测井技术》;20100831;第34卷(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147051A (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
CN109932297B (zh) | 一种致密砂岩储层渗透率的计算方法 | |
CN107292074B (zh) | 一种井间连通性判断方法 | |
NO339596B1 (no) | Kalibrering av en permeabilitetsmodell ved bruk av målinger av reservoarfluidstrømning | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN108194077B (zh) | 一种气测全烃校正方法 | |
CN110162851B (zh) | 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法 | |
CN110043254B (zh) | 一种基于电缆地层测试资料地层有效渗透率的获取方法 | |
CN105604546B (zh) | 双重介质碳酸盐岩储层的定量分类方法 | |
CN105445441A (zh) | 储层含气饱和度确定方法及装置 | |
CN108894778B (zh) | 一种利用气测录井资料识别油气层流体性质的方法 | |
CN105447298A (zh) | 一种储层分析的方法及其应用 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN102562051B (zh) | 处理气测图版原始数据的方法 | |
CN111520126A (zh) | 密度测井扩径校正方法及装置 | |
CN109583113B (zh) | 一种岩石地层压实系数和有效孔隙体积压缩系数计算方法 | |
CN114370269A (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN112147051B (zh) | 一种基于渗透率分布形态的测压流度标准化方法 | |
CN110322363A (zh) | 页岩气储层改造体积计算方法及系统 | |
CN110017136B (zh) | 一种基于视水层电阻率的水淹层识别与产水率预测方法 | |
CN111241652B (zh) | 一种确定地层原油粘度的方法及装置 | |
CN105089632A (zh) | 一种高温高压储层co2流体纵波时差骨架参数的获取方法 | |
CN105298477A (zh) | 一种基于流动单元的地层孔隙结构解释方法 | |
CN114060015B (zh) | 一种致密砂岩含气性的评价方法及装置 | |
CN110792425B (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 |