CN108979629B - 一种基于密度和中子和核磁共振测井的气层密度计算方法 - Google Patents

一种基于密度和中子和核磁共振测井的气层密度计算方法 Download PDF

Info

Publication number
CN108979629B
CN108979629B CN201810751517.2A CN201810751517A CN108979629B CN 108979629 B CN108979629 B CN 108979629B CN 201810751517 A CN201810751517 A CN 201810751517A CN 108979629 B CN108979629 B CN 108979629B
Authority
CN
China
Prior art keywords
density
gas
neutron
expressed
logging
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
Application number
CN201810751517.2A
Other languages
English (en)
Other versions
CN108979629A (zh
Inventor
张恒荣
胡向阳
吴一雄
王丽
陈嵘
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC China Ltd Zhanjiang Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201810751517.2A priority Critical patent/CN108979629B/zh
Publication of CN108979629A publication Critical patent/CN108979629A/zh
Application granted granted Critical
Publication of CN108979629B publication Critical patent/CN108979629B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Fluid Mechanics (AREA)
  • Operations Research (AREA)
  • Environmental & Geological Engineering (AREA)
  • Algebra (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于密度和中子和核磁共振测井的气层密度计算方法，具体步骤为：S1、三种测井曲线获取：选用随钻中子、密度、核磁同一趟测量最佳，或采用电缆中子、密度、核磁同一趟测量，保证测量的对象地层是时间统一、侵入深度接近；S2、密度测井测量的信息可用公式；S3、中子测井测量的信息可用公式；S4、核磁共振测井测量的信息可用公式；S5、通过联立中子、密度、核磁共振的三个方程可以计算出三个未知数，气层的真实总孔隙度
Figure DDA0001725666150000011
含气饱和度(1‑Sxo)和气层密度ρh；S6、输出计算结果。本发明解决了含气储层测井评价提供精确可靠的总孔隙度、冲洗带含气饱和度和气体密度，提高地质评价准确度，模型计算结果准确、适用范围广。

Description

一种基于密度和中子和核磁共振测井的气层密度计算方法
技术领域
本发明涉及含气储层的流体识别及孔隙度、流体密度测算技术领域，具体为一种基于密度和中子和核磁共振测井的气层密度计算方法。
背景技术
在气田的测井储层评价中，准确计算气层的真实总孔隙度、气体密度和冲洗带含气饱和度一直是个难题；比较传统的方法是利用中子密度曲线联合求解，一般沉积岩中的含氢指数和电子密度指数决定了中子测井和密度测井响应，它们之间是相关的，以氢元素的存在而建立联系，且中子密度测井响应都与岩石体积密度相关，因此可利用这两种测井信息计算孔隙度和识别流体；但是利用两个方程求解三个未知数可靠性很差，因此常常代入先验的气体密度参数，气体密度值一般通过地层测试资料获得，成本比较昂贵，且缺乏每口井气体密度可能发生变化，因此中子密度联合计算方法在定量得到冲洗带含气饱和度和气体密度参数的应用上一直有缺陷。随着核磁共振测井的发展，核磁共振具有只对流体敏感的优点，测量的对象也是含氢指数，与中子测井不同的是，核磁共振基本不受骨架中含氢的影响，且测量的地层深度与中子密度测井基本一致，因此可以将核磁共振测井与中子密度测井联合应用，针对以上所述，在这里我们提出一种基于密度和中子和核磁共振测井的气层密度计算方法。
发明内容
本发明的目的是针对上述技术现状，旨在为含气储层测井评价提供精确可靠的总孔隙度、冲洗带含气饱和度和气体密度，提高地质评价准确度。
本发明目的的实现方式为，一种基于密度和中子和核磁共振测井的气层密度计算方法，具体步骤为：
S1、三种测井曲线获取：选用随钻中子、密度、核磁同一趟测量最佳，或采用电缆中子、密度、核磁同一趟测量，保证测量的对象地层是时间统一、侵入深度接近；
S2、密度测井测量的信息可用公式表示为：
Figure GDA0003175672210000021
式中：
Figure GDA0003175672210000022
Figure GDA0003175672210000023
Figure GDA0003175672210000025
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA0003175672210000026
为密度测量视孔隙度，用％表示；
ρma为岩石骨架密度，用g/cm3表示，ρb为密度测量曲线值，用g/cm3表示，ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
S3、中子测井测量的信息可用公式表示为：
Figure GDA0003175672210000024
式中：
Figure GDA0003175672210000031
Figure GDA0003175672210000032
Figure GDA0003175672210000033
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA0003175672210000034
为中子测量视孔隙度，用％表示；
Figure GDA0003175672210000035
为岩石骨架中子，用pu表示，
Figure GDA0003175672210000036
为中子测量曲线值，用pu表示，
Figure GDA0003175672210000037
为孔隙中水的流体中子，用pu表示；
ρh为气层的流体密度，用g/cm3表示；
S4、核磁共振测井测量的信息可用公式表示为：
Figure GDA0003175672210000038
式中：
P＝1-e-Tw/T1
Figure GDA0003175672210000039
Figure GDA00031756722100000310
Figure GDA00031756722100000311
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA00031756722100000312
为核磁共振测量视孔隙度，用％表示；
ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
Mw为流体的相对分子质量，无量纲；
Tw为核磁测量等待时间，用S表示；
T1为纵向核磁弛豫时间，用S表示；
S5、通过联立中子、密度、核磁共振的三个方程可以计算出三个未知数，气层的真实总孔隙度
Figure GDA0003175672210000041
含气饱和度(1-Sxo)和气层密度ρh；
S6、输出计算结果。
作为本发明的进一步优选方式，还包括，基于三种孔隙度测井响应模型计算出气层真实的总孔隙度、含气饱和度和气体密度,确定模型中
Figure GDA0003175672210000049
Sxo、ρh参数，确定方法包括：
定义密度测井响应方程：
Figure GDA0003175672210000042
定义中子测井响应方程：
Figure GDA0003175672210000043
Figure GDA0003175672210000044
定义核磁共振测井响应方程：
Figure GDA0003175672210000045
Figure GDA0003175672210000046
上述三式中：
Figure GDA0003175672210000047
Sxo、ρh为未知数，其他参数均为测量值或固定参数，如三种孔隙度测井响应的骨架参数、水的流体参数和核磁共振测井的测量模式参数等待时间；
联合上述三个式子，通过迭代方法，选取气体密度在0-1之间，进行0.01步长的循环，即可求解出
Figure GDA0003175672210000048
参数。
与现有技术相比，本发明的有益效果如下：
本发明解决了含气储层测井评价提供精确可靠的总孔隙度、冲洗带含气饱和度和气体密度，提高地质评价准确度，模型计算结果准确、适用范围广，本发明适合在南海西部海域计算含气储层总孔隙度、冲洗带含气饱和度和气体密度曲线，流体识别符合率较高，能较好满足地质研究需要，值得推广。
附图说明
图1为本发明实际工作流程详图；
图2为发明应用成果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1-2，本实用发明提供一种技术方案：一种基于密度和中子和核磁共振测井的气层密度计算方法，具体步骤为：
S1、三种测井曲线获取：选用随钻中子、密度、核磁同一趟测量最佳，或采用电缆中子、密度、核磁同一趟测量，保证测量的对象地层是时间统一、侵入深度接近；
S2、密度测井测量的信息可用公式表示为：
Figure GDA0003175672210000051
式中：
Figure GDA0003175672210000061
Figure GDA0003175672210000062
Figure GDA0003175672210000063
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA0003175672210000064
为密度测量视孔隙度，用％表示；
ρma为岩石骨架密度，用g/cm3表示，ρb为密度测量曲线值，用g/cm3表示，ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
S3、中子测井测量的信息可用公式表示为：
Figure GDA0003175672210000065
式中：
Figure GDA0003175672210000066
Figure GDA0003175672210000067
Figure GDA0003175672210000068
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA0003175672210000069
为中子测量视孔隙度，用％表示；
Figure GDA00031756722100000610
为岩石骨架中子，用pu表示，
Figure GDA00031756722100000611
为中子测量曲线值，用pu表示，
Figure GDA00031756722100000612
为孔隙中水的流体中子，用pu表示；
ρh为气层的流体密度，用g/cm3表示；
S4、核磁共振测井测量的信息可用公式表示为：
Figure GDA0003175672210000071
式中：
P＝1-e-Tw/T1
Figure GDA0003175672210000072
Figure GDA0003175672210000073
Figure GDA0003175672210000074
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure GDA0003175672210000075
为核磁共振测量视孔隙度，用％表示；
ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
Mw为流体的相对分子质量，无量纲；
Tw为核磁测量等待时间，用S表示；
T1为纵向核磁弛豫时间，用S表示；
S5、通过联立中子、密度、核磁共振的三个方程可以计算出三个未知数，气层的真实总孔隙度
Figure GDA0003175672210000076
含气饱和度(1-Sxo)和气层密度ρh；
S6、输出计算结果。
还包括基于三种孔隙度测井响应模型计算出气层真实的总孔隙度、含气饱和度和气体密度,确定模型中
Figure GDA0003175672210000077
Sxo、ρh参数，确定方法包括：
定义密度测井响应方程：
Figure GDA0003175672210000081
定义中子测井响应方程：
Figure GDA0003175672210000082
Figure GDA0003175672210000083
定义核磁共振测井响应方程：
Figure GDA0003175672210000084
Figure GDA0003175672210000085
上述三式中：
Figure GDA0003175672210000086
Sxo、ρh为未知数，其他参数均为测量值或固定参数，如三种孔隙度测井响应的骨架参数、水的流体参数和核磁共振测井的测量模式参数等待时间；
联合上述三个式子，通过迭代方法，选取气体密度在0-1之间，进行0.01步长的循环，即可求解出
Figure GDA0003175672210000087
Sxo、ρh参数。
综上述，本发明解决了含气储层测井评价提供精确可靠的总孔隙度、冲洗带含气饱和度和气体密度，提高地质评价准确度，模型计算结果准确、适用范围广，本发明适合在南海西部海域计算含气储层总孔隙度、冲洗带含气饱和度和气体密度曲线，流体识别符合率较高，能较好满足地质研究需要，值得推广。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点,对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (2)

1.一种基于密度和中子和核磁共振测井的气层密度计算方法，其特征在于，具体步骤为：
S1、三种测井曲线获取：选用随钻中子、密度、核磁同一趟测量最佳，或采用电缆中子、密度、核磁同一趟测量，保证测量的对象地层是时间统一、侵入深度接近；
S2、密度测井测量的信息可用公式表示为：
Figure FDA0003241880350000011
式中：
Figure FDA0003241880350000012
Figure FDA0003241880350000013
Figure FDA0003241880350000014
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure FDA0003241880350000015
为密度测量视孔隙度，用％表示；
ρma为岩石骨架密度，用g/cm3表示，ρb为密度测量曲线值，用g/cm3表示，ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
S3、中子测井测量的信息可用公式表示为：
Figure FDA0003241880350000016
式中：
Figure FDA0003241880350000017
Figure FDA0003241880350000021
Figure FDA0003241880350000022
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure FDA0003241880350000023
为中子测量视孔隙度，用％表示；
Figure FDA0003241880350000024
为岩石骨架中子，用pu表示，
Figure FDA0003241880350000025
为中子测量曲线值，用pu表示，
Figure FDA0003241880350000026
为孔隙中水的流体中子，用pu表示；
ρh为气层的流体密度，用g/cm3表示；
S4、核磁共振测井测量的信息可用公式表示为：
Figure FDA0003241880350000027
式中：
P＝1-e-Tw/T1
Figure FDA0003241880350000028
Figure FDA0003241880350000029
Figure FDA00032418803500000210
为地层真实总孔隙度，用％表示；
Sxo为近井带气层含水饱和度，用小数表示；
Figure FDA00032418803500000211
为核磁共振测量视孔隙度，用％表示；
ρmf为孔隙中水的流体密度，用g/cm3表示；
ρh为气层的流体密度，用g/cm3表示；
Mw为流体的相对分子质量，无量纲；
Tw为核磁测量等待时间，用S表示；
T1为纵向核磁弛豫时间，用S表示；
S5、通过联立中子、密度、核磁共振的三个方程可以计算出三个未知数，气层的真实总孔隙度
Figure FDA0003241880350000031
含气饱和度(1-Sxo)和气层密度ρh；
S6、输出计算结果。
2.根据权利要求1所述的一种基于密度和中子和核磁共振测井的气层密度计算方法，其特征在于，还包括，基于三种孔隙度测井响应模型计算出气层真实的总孔隙度、含气饱和度和气体密度,确定模型中
Figure FDA0003241880350000032
Sxo、ρh参数，确定方法包括：
定义密度测井响应方程：
Figure FDA0003241880350000033
定义中子测井响应方程：
Figure FDA0003241880350000034
Figure FDA0003241880350000035
定义核磁共振测井响应方程：
Figure FDA0003241880350000036
Figure FDA0003241880350000037
上述三式中：
Figure FDA0003241880350000038
Sxo、ρh为未知数，其他参数均为测量值或固定参数；
联合上述三个式子，通过迭代方法，选取气体密度在0-1之间，进行0.01步长的循环，即可求解出
Figure FDA0003241880350000039
Sxo、ρh参数。
CN201810751517.2A 2018-07-10 2018-07-10 一种基于密度和中子和核磁共振测井的气层密度计算方法 Active CN108979629B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810751517.2A CN108979629B (zh) 2018-07-10 2018-07-10 一种基于密度和中子和核磁共振测井的气层密度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810751517.2A CN108979629B (zh) 2018-07-10 2018-07-10 一种基于密度和中子和核磁共振测井的气层密度计算方法

Publications (2)

Publication Number Publication Date
CN108979629A CN108979629A (zh) 2018-12-11
CN108979629B true CN108979629B (zh) 2021-10-26

Family

ID=64537691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810751517.2A Active CN108979629B (zh) 2018-07-10 2018-07-10 一种基于密度和中子和核磁共振测井的气层密度计算方法

Country Status (1)

Country Link
CN (1) CN108979629B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111812736B (zh) * 2020-07-22 2021-03-16 西南石油大学 一种致密砂岩无水气藏含气性评价方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1325494A (zh) * 1998-12-30 2001-12-05 贝克休斯公司 由钻孔电阻率图像仪、横向感应测井纪录以及含水饱和度张量模型确定水饱和度及砂岩分数的方法
CN105426650A (zh) * 2015-06-08 2016-03-23 中国石油集团川庆钻探工程有限公司 利用测井资料优选页岩气井压裂改造井段的方法
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN108049866A (zh) * 2017-12-11 2018-05-18 中石化石油工程技术服务有限公司 二维核磁共振测井致密气藏定量评价方法
CN108138555A (zh) * 2015-02-23 2018-06-08 奈克森能量无限责任公司 预测储层性质的方法、系统及设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2016005569A (es) * 2013-12-19 2016-10-26 Halliburton Energy Services Inc Integracion de porosidad de resonancia magnetica nuclear (rmn) en una metodologia de interpretacion probabilistica de registros multiples.

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1325494A (zh) * 1998-12-30 2001-12-05 贝克休斯公司 由钻孔电阻率图像仪、横向感应测井纪录以及含水饱和度张量模型确定水饱和度及砂岩分数的方法
CN108138555A (zh) * 2015-02-23 2018-06-08 奈克森能量无限责任公司 预测储层性质的方法、系统及设备
CN105426650A (zh) * 2015-06-08 2016-03-23 中国石油集团川庆钻探工程有限公司 利用测井资料优选页岩气井压裂改造井段的方法
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN108049866A (zh) * 2017-12-11 2018-05-18 中石化石油工程技术服务有限公司 二维核磁共振测井致密气藏定量评价方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
综合常规与核磁共振测井资料评价储层产液性质;张恒荣等;《测井技术》;20141031;第38卷(第5期);第574-580页 *

Also Published As

Publication number Publication date
CN108979629A (zh) 2018-12-11

Similar Documents

Publication Publication Date Title
EP2965126B1 (en) Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data
CN109753755B (zh) 一种确定储层含水饱和度的方法
Mao et al. Estimation of permeability by integrating nuclear magnetic resonance (NMR) logs with mercury injection capillary pressure (MICP) data in tight gas sands
CN103422853B (zh) 一种确定地层含水饱和度的方法及装置
CN104278989B (zh) 一种获取低孔低渗储层饱和度指数的方法
US8538700B2 (en) Method of determining subterranean formation parameters
CN108694264B (zh) 一种确定页岩气储层渗透率的方法
CN110231272A (zh) 致密砂岩孔径与核磁共振t2值转换关系的确定方法及系统
CN109100812B (zh) 基于核磁共振的岩石孔隙分形维数评价方法及装置
CN107688037B (zh) 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法
US8005619B2 (en) Method of determining reservoir parameters
CN106897531B (zh) 一种低渗透石灰岩储层渗透率的定量评价方法
CN103615230B (zh) 一种双泥质指示因子含水饱和度模型的建立方法
CN102052074A (zh) 一种岩心数字化确定岩电参数的方法
Garcia et al. Integrated workflow to estimate permeability through quantification of rock fabric using joint interpretation of nuclear magnetic resonance and electric measurements
CN108979629B (zh) 一种基于密度和中子和核磁共振测井的气层密度计算方法
Brun et al. Use of prior information in gradient-based history matching
CN105844011A (zh) 一种基于毛管模型的渗透率计算方法
CN113775326A (zh) 可动水饱和度的评价方法、装置、电子设备及介质
CN103675945B (zh) 一种测定孔洞型储层的饱和度的方法及设备
CN114060015B (zh) 一种致密砂岩含气性的评价方法及装置
CN114086938B (zh) 一种非均质性砂岩储层的含气饱和度预测方法
CN111475982B (zh) 一种岩石内部磁场梯度加权几何平均值的三因素估算方法
CN114035236A (zh) 基于新三水模型的核磁和电阻率联合反演含水饱和度方法
CN115901563A (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