CN109138975B - 一种基于时移测井数据的求解相渗特征曲线的新方法 - Google Patents
一种基于时移测井数据的求解相渗特征曲线的新方法 Download PDFInfo
- Publication number
- CN109138975B CN109138975B CN201810750441.1A CN201810750441A CN109138975B CN 109138975 B CN109138975 B CN 109138975B CN 201810750441 A CN201810750441 A CN 201810750441A CN 109138975 B CN109138975 B CN 109138975B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- calculating
- time
- phase
- formation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000035699 permeability Effects 0.000 claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 32
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 31
- 238000005259 measurement Methods 0.000 claims abstract description 22
- 239000012530 fluid Substances 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 238000005553 drilling Methods 0.000 claims description 25
- 230000009545 invasion Effects 0.000 claims description 13
- 239000008398 formation water Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000000704 physical effect Effects 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 3
- 239000011435 rock Substances 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 abstract description 5
- 238000012360 testing method Methods 0.000 abstract description 4
- 208000035126 Facies Diseases 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000010606 normalization Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000706 filtrate Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000005325 percolation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
Images
Classifications
-
- 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
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
- E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于时移测井数据的求解相渗特征曲线的新方法,步骤1,根据不同时刻电阻率测量数据,基于联合反演理论反演地层电阻率剖面信息;步骤2,格式化输出电阻率剖面信息;步骤3,根据电阻率和饱和度转换公式计算出不同时刻地层饱和度剖面;步骤4,计算任意两时刻间地层含水饱和度的改变量;步骤5,计算单位时间内地层流体的流量;步骤6,计算储层中各相流体的有效渗透率特征点;步骤7,拟合输出两相流体的相渗曲线,本发明以往获取有效渗透率的最佳途径主要由测试获得,但其成本极高,基于时移测井数据求解相渗特征曲线的新方法,将为储层的有效性和产能评价提供关键的基础参数,是勘探开发降本增效的一项重要举措。
Description
技术领域
本发明涉及油田开发技术领域,具体为一种基于时移测井数据的求解相渗特征曲线的新方法。
背景技术
储层的渗流特征是控制完井、油田开发及油藏管理决策的一个关键参数。对储层渗流特征评价的最重要参数是有效渗透率,获取有效渗透率的最佳途径主要由测试获得,但其成本极高,因此一般由实验室的绝对渗透率进行转换,比如岩心分析、测井对比以及井测试等。岩心分析提供的是相当于几立方英寸的样品的平均渗透率。这种样品的尺寸太小,无法概括地描述储层特征。标准的井测试所研究的样品的尺寸约是几千立方英尺。这种样品的尺寸太大,不能对储层的非均质性做出任何详细的描述。
发明内容
本发明的目的在于提供一种基于时移测井数据的求解相渗特征曲线的新方法,解决了背景技术中所提出的问题。
为实现上述目的,本发明提供如下技术方案:一种基于时移测井数据的相渗特征曲线求解新方法,其主要特征在于:
步骤1,根据不同时刻电阻率测量数据,基于联合反演理论反演地层电阻率剖面信息;
步骤2,格式化输出电阻率剖面信息;
步骤3,根据电阻率和饱和度转换公式计算出不同时刻地层饱和度剖面;
步骤4,计算任意两时刻间地层含水饱和度的改变量;
步骤5,计算单位时间内地层流体的流量;
步骤6,计算储层中各相流体的有效渗透率特征点;
步骤7,拟合输出两相流体的相渗曲线。
所述的根据时移电阻率测量数据联合反演地层电阻率剖面信息,还包括,涉及一种时移电阻率测量数据联合反演方法,主要步骤包括:
步骤1,测井数据收集及整理,包括不同时刻电阻率测量数据以及对应的测井时间、钻井日志、测井解释成果表、储层物性参数等;
步骤2,数据预处理,包括曲线有效性判别、深度对齐、自动分层等;
步骤3,时移电阻率联合反演,根据不同仪器的测量数据调用相应的正演仿真程序,采用LM迭代反演方法反演获得不同时刻电阻率剖面;
步骤4,输出反演的地层电阻率剖面,包括不同时刻侵入深度、侵入带电阻率和地层真电阻率。
所述的时移测井电阻率联合反演算法:所述联合反演包括同一仪器不同时刻的随钻实时和上提复测测量数据或不同仪器不同时刻的随钻和电缆测量数据。
所述的格式化输出反演的地层电阻率剖面,其特征在于:所述输出的电阻率剖面中地层真电阻率为同一值,不同时刻测井的原因导致泥浆侵入深度和侵入带电阻率不同,所述反演结果包括不同时刻的泥浆侵入深度、不同时刻侵入带电阻率、原状地层真电阻率等信息。
所述的计算不同时刻地层饱和度剖面,其特征在于:所述根据电阻率和饱和度转换关系,计算不同时刻侵入带含水饱和度和原状地层含水饱和度;
所述的计算单位时间内泥浆滤液的侵入量,其特征在于:
所述的由钻井日志中给出的泥浆比重,计算出钻井柱和地层压力,再根据达西公式计算出地层水、油或气的有效渗透率特征点,其特征在于:
步骤1,根据钻井日志提取目标层段的泥浆密度(比重)和地层水密度,由下面压力计算公式计算钻井柱和原状地层得压力:
所述的根据相渗透率的定义,获取某含水饱和度对应的相对渗透率特征点,拟合输出两相流体的相渗曲线,其特征在于:
步骤2,根据目标层计算的含水饱和度和对应的相流体相对渗透率特征散点值,拟合获得含水饱和度和各相流体相对渗透率函数关系。
与现有技术相比,本发明的有益效果如下:
该一种基于时移测井数据的求解相渗特征曲线的新方法,本发明提出一种基于不同时间下的电阻率测井资料、多种动静态测井数据联合表征储层的有效渗透率的新方法,该方法涉及到不同时刻测量数据的有效归一化,需要测井、储层、油藏等多学科的联合攻关,是一项全新的求解相渗特征曲线的新方法。基于时移测井数据求解相渗特征曲线的新方法,将为储层有效性和产能评价提供关键基础参数,能为测井的储层评价技术实现由静至动的飞跃,并推动测井资料的有机集成和时间归一化,是勘探开发降本增效的一项重要举措。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明的基于时移电阻率测量数据的建立相渗曲线的一具体实例流程图;
图2为本发明一实例中不同时刻电阻率测量数据示意图;
图3为本发明一实例中时移电阻率联合反演的电阻率剖面图;
图4为本发明一实例中计算的各目标储层饱和度与相对渗透率对应的特征散点图;
图5为本发明一实例中根据特征散点所拟合的饱和度和相对渗透率关系曲线。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
如图1所示,图1为本发明所涉及的基于时移电阻率测量数据建立相渗特征曲线的流程图。
步骤101,首先基于时移电阻率测量数据开展电阻率联合反演研究,其主要步骤包括:
1),测井数据收集及整理,包括不同时刻电阻率测量数据以及对应的测井时间、钻井日志、测井解释成果表、储层物性参数等,不同时刻电阻率测量数据如图2所示;
2),数据预处理,包括曲线有效性判别、时移电阻率曲线深度对齐、地层结构的自动分层处理等;
3),时移电阻率多核并行联合反演处理,根据不同仪器的测量数据调用相应的正演仿真程序,采用LM并行迭代反演方法反演获得不同时刻电阻率剖面,流程转入步骤102;
步骤102,格式化输出反演的地层电阻率剖面,包括不同时刻侵入深度、侵入带电阻率和地层真电阻率,如图3所示,其中输出的电阻率剖面中地层真电阻率为同一值,不同时刻测井的原因导致泥浆侵入深度和侵入带电阻率不同,转入步骤103;
步骤103,根据电阻率和饱和度转换公式计算出不同时刻地层饱和度剖面,计算不同时刻侵入带含水饱和度和原状地层含水饱和度,流程转入步骤104;
步骤106,根据达西公式计算出地层水、油或气的有效渗透率特征点,计算结果如图4所示,其有效渗透率特征点的具体计算过程如下:
1),根据钻井日志提取目标层段的泥浆密度和地层水密度,由下面压力计算公式计算钻井柱和原状地层得压力:
步骤107,根据相渗透率的定义,获取某含水饱和度对应的相对渗透率特征点,并拟合输出两相流体的相渗曲线,结果如图5所示,其具体计算过程如下:
2),根据步骤103中计算的各目标层含水饱和度及其对应的步骤106计算的水、油或气的相对渗透率特征散点值,拟合获得含水饱和度和各相流体相对渗透率函数关系。
本发明提出一种基于不同时间下的电阻率测井资料、多种动静态测井数据联合表征储层的有效渗透率的新方法,该方法涉及到不同时刻测量数据的有效归一化,需要测井、储层、油藏等多学科的联合攻关,是一项全新的求解相渗特征曲线的新方法。基于时移测井数据求解相渗特征曲线的新方法,将为储层有效性和产能评价提供关键基础参数,能为测井的储层评价技术实现由静至动的飞跃,并推动测井资料的有机集成和时间归一化,是勘探开发降本增效的一项重要举措。
工作原理:测井数据收集及整理,包括不同时刻电阻率测量数据以及对应的测井时间、钻井日志、测井解释成果表、储层物性参数等,不同时刻电阻率测量数据,数据预处理,包括曲线有效性判别、时移电阻率曲线深度对齐、地层结构的自动分层处理等,时移电阻率多核并行联合反演处理,根据不同仪器的测量数据调用相应的正演仿真程序,采用LM并行迭代反演方法反演获得不同时刻电阻率剖面。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点,对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (6)
1.一种基于时移测井数据的相渗特征曲线求解新方法,其特征在于:
步骤1,根据不同时刻电阻率测量数据,基于联合反演理论反演地层电阻率剖面信息,主要步骤包括:
步骤a,测井数据收集及整理,包括不同时刻电阻率测量数据以及对应的测井时间、钻井日志、测井解释成果表、储层物性参数;
步骤b,数据预处理,包括曲线有效性判别、深度对齐、自动分层;
步骤c,时移电阻率联合反演,根据不同仪器的测量数据调用相应的正演仿真程序,采用LM迭代反演方法反演获得不同时刻电阻率剖面;
步骤2,格式化输出电阻率剖面信息;
步骤3,根据电阻率和饱和度转换公式计算出不同时刻地层饱和度剖面;
步骤4,计算任意两时刻间地层含水饱和度的改变量;
步骤5,计算单位时间内地层流体的流量;
步骤6,计算储层中各相流体的有效渗透率特征点,主要步骤包括:
步骤a,根据钻井日志提取目标层段的泥浆密度和地层水密度,由下面压力计算公式计算钻井柱和原状地层的压力:
步骤7,拟合输出两相流体的相渗曲线。
2.根据权利要求1所述的一种基于时移测井数据的相渗特征曲线求解新方法,其特征在于:时移测井电阻率联合反演算法,联合反演包括同一仪器不同时刻的随钻实时和上提复测测量数据或不同仪器不同时刻的随钻和电缆测量数据。
3.根据权利要求1所述的一种基于时移测井数据的相渗特征曲线求解新方法,其特征在于:步骤2格式化输出反演的地层电阻率剖面,包括不同时刻侵入深度、侵入带电阻率和地层真电阻率,格式化输出反演的地层电阻率剖面方法,所述输出的电阻率剖面中地层真电阻率为同一值,不同时刻测井的原因导致泥浆侵入深度和侵入带电阻率不同,所述反演结果包括不同时刻的泥浆侵入深度、不同时刻侵入带电阻率、原状地层真电阻率信息。
4.根据权利要求1所述的一种基于时移测井数据的相渗特征曲线求解新方法,其特征在于:步骤3计算不同时刻地层饱和度剖面方法,所述根据电阻率和饱和度转换关系,计算不同时刻侵入带含水饱和度和原状地层含水饱和度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810750441.1A CN109138975B (zh) | 2018-07-10 | 2018-07-10 | 一种基于时移测井数据的求解相渗特征曲线的新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810750441.1A CN109138975B (zh) | 2018-07-10 | 2018-07-10 | 一种基于时移测井数据的求解相渗特征曲线的新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109138975A CN109138975A (zh) | 2019-01-04 |
CN109138975B true CN109138975B (zh) | 2022-08-02 |
Family
ID=64800231
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810750441.1A Active CN109138975B (zh) | 2018-07-10 | 2018-07-10 | 一种基于时移测井数据的求解相渗特征曲线的新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109138975B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109989743B (zh) * | 2019-05-06 | 2022-09-06 | 中海油田服务股份有限公司 | 一种确定泥浆滤液侵入深度的方法及系统 |
CN110043254B (zh) * | 2019-05-11 | 2022-11-08 | 中国海洋石油集团有限公司 | 一种基于电缆地层测试资料地层有效渗透率的获取方法 |
CN111965720B (zh) * | 2020-08-19 | 2023-05-23 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种基于地-井联合获取水力传导系数的方法 |
CN112199852A (zh) * | 2020-10-22 | 2021-01-08 | 中海油田服务股份有限公司 | 一种定向射孔的变孔密方法、存储介质和电子装置 |
CN116398126B (zh) * | 2023-06-07 | 2023-08-25 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101487390A (zh) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | 一种确定油层原始含油饱和度的阿尔奇模式方法 |
CN101649738A (zh) * | 2008-08-13 | 2010-02-17 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法 |
CN103645126A (zh) * | 2013-12-02 | 2014-03-19 | 西南石油大学 | 地层高温高压气水相渗曲线测定方法 |
CN104018829A (zh) * | 2014-05-23 | 2014-09-03 | 中国地质大学(北京) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 |
CN104196524A (zh) * | 2014-07-14 | 2014-12-10 | 中国地质大学(北京) | 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法 |
CN104895561A (zh) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | 基于并行计算技术的随钻电阻率测井联合反演方法 |
CN104912550A (zh) * | 2015-05-05 | 2015-09-16 | 中国海洋石油总公司 | 一种核磁共振测井资料定量计算储层产液剖面方法 |
-
2018
- 2018-07-10 CN CN201810750441.1A patent/CN109138975B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101649738A (zh) * | 2008-08-13 | 2010-02-17 | 中国石油天然气集团公司 | 一种确定地层含水饱和度的方法 |
CN101487390A (zh) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | 一种确定油层原始含油饱和度的阿尔奇模式方法 |
CN103645126A (zh) * | 2013-12-02 | 2014-03-19 | 西南石油大学 | 地层高温高压气水相渗曲线测定方法 |
CN104018829A (zh) * | 2014-05-23 | 2014-09-03 | 中国地质大学(北京) | 一种利用煤层气井生产数据测量气水相渗曲线的方法 |
CN104196524A (zh) * | 2014-07-14 | 2014-12-10 | 中国地质大学(北京) | 一种欠饱和煤储层开发的气水产出动态相渗曲线测量方法 |
CN104895561A (zh) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | 基于并行计算技术的随钻电阻率测井联合反演方法 |
CN104912550A (zh) * | 2015-05-05 | 2015-09-16 | 中国海洋石油总公司 | 一种核磁共振测井资料定量计算储层产液剖面方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109138975A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CN105468886B (zh) | 基于岩石物理学参数计算地层压力的方法 | |
Quinn et al. | Using constant head step tests to determine hydraulic apertures in fractured rock | |
CN101487390B (zh) | 一种确定油层原始含油饱和度的阿尔奇模式方法 | |
Wu et al. | Numerical simulation of mud-filtrate invasion in deviated wells | |
CN108694264B (zh) | 一种确定页岩气储层渗透率的方法 | |
CN104213899B (zh) | 一种地层岩石骨架的测井识别方法 | |
CA2867580C (en) | Fracking method at fracking intervals of a horizontal drilling zone in a sweet spot range determined based on measurements of resistivity and density in the horizontal drilling zone | |
CN104899411B (zh) | 一种储层产能预测模型建立方法和系统 | |
CN103628871A (zh) | 一种基于阿尔奇公式的电阻率侵入校正的新方法 | |
CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
CN106223941B (zh) | 基于测井资料的页岩气储层有机孔隙度确定方法 | |
CN103698811A (zh) | 一种碳酸盐岩岩石结构组分测井定量识别方法及其用途 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN102621586B (zh) | 一种识别地层属性的地层数据处理方法 | |
Odling et al. | Fractional flow in fractured chalk; a flow and tracer test revisited | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
CN109386285A (zh) | 一种特低渗透油层泥浆侵入程度及其影响的评价方法 | |
Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
Shchipanov et al. | A new approach to deformable fractured reservoir characterization: case study of the Ekofisk field | |
George et al. | Estimation of aquifer hydraulic parameters via complementing surfacial geophysical measurement by laboratory measurements on the aquifer core samples | |
Pontiggia et al. | New integrated approach for diagenesis characterization and simulation | |
CN107808055B (zh) | 一种基于双重影响的页岩气饱和度测井计算方法 | |
CN115099014A (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 |