CN104533405B - 泥页岩储层有机碳含量的检测方法和装置 - Google Patents
泥页岩储层有机碳含量的检测方法和装置 Download PDFInfo
- Publication number
- CN104533405B CN104533405B CN201410799475.1A CN201410799475A CN104533405B CN 104533405 B CN104533405 B CN 104533405B CN 201410799475 A CN201410799475 A CN 201410799475A CN 104533405 B CN104533405 B CN 104533405B
- Authority
- CN
- China
- Prior art keywords
- hydrodynamic
- hydrodynamic zone
- carbon content
- organic carbon
- zone
- 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
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 149
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 149
- 238000001514 detection method Methods 0.000 title claims abstract description 37
- 238000012360 testing method Methods 0.000 claims abstract description 90
- 239000011435 rock Substances 0.000 claims abstract description 79
- 230000005251 gamma ray Effects 0.000 claims abstract description 73
- 238000003860 storage Methods 0.000 claims description 9
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 22
- 239000003208 petroleum Substances 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000004575 stone Substances 0.000 description 3
- 230000008021 deposition Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
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
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)
- Analysing Materials By The Use Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种泥页岩储层有机碳含量的检测方法和装置，涉及石油地质勘探技术领域。方法包括：根据待测的泥页岩储层的岩心数据中的水动力数据，将待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带；获取各水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值，以及确定各水动力带中的泥页岩岩样的有机碳含量；确定各水动力带对应的测试参数；分别计算各水动力带中的待测的泥页岩储层的有机碳含量。本发明能够解决现有技术中基于△logR方法检测有机碳含量繁琐复杂，且检测结果不准确的问题。
Description
技术领域
本发明涉及石油地质勘探技术领域，尤其涉及一种泥页岩储层有机碳含量的检测方法和装置。
背景技术
在石油开采领域中，页岩气作为一种非常规天然气资源，正在被广泛研究，而泥页岩中的有机碳含量(简称TOC)更是当前国内外研究的重点。
目前，现有技术中在通过测井资料检测有机碳含量时大多使用△logR方法或者是基于△logR方法。△logR技术是一种利用测井资料识别和计算含有机质岩层的有机碳含量的方法。当前的△logR的计算方法可以由公式(1-1)表示：
△logR＝log(R/R基线)+0.02(△t-△t基线)(1-1)
其中：△logR是测井曲线中电阻率与声波时差曲线叠合幅度差；R是实测的电阻率；R基线是基线对应的电阻率；△t是实测的声波时差；△t基线是基线对应的声波时差；0.02是对数坐标下的一个电阻率单位与算术坐标下一个声波时差周期50μs/m的比值。
由于△logR与有机碳含量呈线性相关,并且是成熟度的函数，因此当前可以得到的△logR计算有机碳含量的经验公式为如下公式(1-2)：
TOC＝△logR×10(2.297-0.l688LOM)+△TOC (1-2)
其中，TOC为计算的有机碳含量；LOM反映有机质成熟度，△TOC为有机碳含量背景值。
当前，在确定R基线和△t基线时，一般需要人为来根据测井曲线中的电阻率曲线和声波时差曲线进行绘制，这样造成△logR的不确定性，进而根据△logR所得到的有机碳含量不准确，且整个检测过程繁琐缓慢。另外，当前的有机碳含量的检测是采用一条或多条测井曲线资料对整段泥页岩进行检测，由于同一段泥页岩有机质的赋存状态不同，有机碳的含量差异很大，造成当前通过一条或多条测井曲线资料对整段泥页岩进行检测的结果不准确。
发明内容
本发明的实施例提供一种泥页岩储层有机碳含量的检测方法和装置方法和装置，以解决现有技术中基于△logR方法检测有机碳含量繁琐复杂，且检测结果不准确的问题。
为达到上述目的，本发明采用如下技术方案：
一种泥页岩储层有机碳含量的检测方法，其特征在于，包括：
根据待测的泥页岩储层的岩心数据中的水动力数据，将所述待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带；
根据待测的泥页岩储层的测井数据获取所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；
根据所述岩心数据分别确定所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量；
根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数；
根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
具体的，所述根据所述各测试参数、第一水动力带、第二水动力带以及第二水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算所述第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1。
具体的，所述根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算所述第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2。
具体的，所述根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算所述第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3。
具体的，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4。
具体的，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5。
具体的，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6。
一种泥页岩储层有机碳含量的检测装置，其特征在于，包括：
泥页岩储层分段单元，用于根据待测的泥页岩储层的岩心数据中的水动力数据，将所述待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带；
泥页岩岩样曲线值获取单元，用于根据待测的泥页岩储层的测井数据获取所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；
泥页岩岩样有机碳含量确定单元，用于根据所述岩心数据分别确定所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量；
测试参数确定单元，用于根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数；
泥页岩储层有机碳含量计算单元，用于根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
另外，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算所述第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1。
另外，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算所述第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2。
另外，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算所述第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3。
此外，所述测试参数确定单元，具体用于：
根据所述第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4。
此外，所述测试参数确定单元，具体用于：
根据所述第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5。
此外，所述测试参数确定单元，具体用于：
根据所述第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6。
本发明实施例提供的泥页岩储层有机碳含量的检测方法和装置，通过水动力数据将待测的泥页岩储层划分为多个水动力带区域，并根据测井数据获取各水动力带区域中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；并确定各水动力带区域中的泥页岩岩样的有机碳含量；在确定各水动力带区域的各测试参数之后，能够根据各水动力带区域待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，以及各测试参数，计算各水动力带区域中的待测的泥页岩储层的有机碳含量。可见本发明无需人为确定R基线和△t基线，且将一整段泥页岩根据水动力数据进行划分水动力带区域，进而检测各区域的有机碳含量，避免了现有技术中有机碳含量不准确，且整个检测过程繁琐缓慢的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的泥页岩储层有机碳含量的检测方法的流程图；
图2为本发明实施例中的计算TOC和测试TOC的相关性曲线示意图一；
图3为本发明实施例中的计算TOC和测试TOC的相关性曲线示意图二；
图4为本发明实施例中的计算TOC和测试TOC的相关性曲线示意图三；
图5为本发明实施例中的TOC曲线与△logR方法的TOC曲线对比图；
图6为本发明实施例提供的泥页岩储层有机碳含量的检测装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供的一种泥页岩储层有机碳含量的检测方法，包括：
步骤101、根据待测的泥页岩储层的岩心数据中的水动力数据，将待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带。
步骤102、根据待测的泥页岩储层的测井数据获取第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值。
步骤103、根据岩心数据分别确定第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量。
步骤104、根据第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数。
步骤105、根据各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
本发明实施例提供的泥页岩储层有机碳含量的检测方法，通过水动力数据将待测的泥页岩储层划分为多个水动力带区域，并根据测井数据获取各水动力带区域中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；并确定各水动力带区域中的泥页岩岩样的有机碳含量；在确定各水动力带区域的各测试参数之后，能够根据各水动力带区域待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，以及各测试参数，计算各水动力带区域中的待测的泥页岩储层的有机碳含量。可见本发明无需人为确定R基线和△t基线，且将一整段泥页岩根据水动力数据进行划分水动力带区域，进而检测各区域的有机碳含量，避免了现有技术中有机碳含量不准确，且整个检测过程繁琐缓慢的问题。
值得说明的是，由于不同沉积构造的泥页岩沉积时的水动力条件不同，其沉积环境存在很大差异，如水底扰动情况、含氧量的不同，导致了沉积时有机质的保存条件不同，进而一个泥页岩段的有机碳含量差异很大。因此，本发明实施例采用将一整段泥页岩根据水动力数据进行划分水动力带区域，进而检测各区域的有机碳含量，使得检测结果较为准确。
上述的第一水动力带属于强水动力带，其对应的泥页岩段为块状泥岩段，例如泥质粉砂岩；上述的第二水动力带对应的泥页岩段为弱纹层状泥岩段，例如钙质弱纹层状泥页岩、粉砂质弱纹层泥页岩；上述第三水动力带属于弱水动力带，其对应的泥页岩段为纹层状泥页岩、碳质纹层泥页岩。
在本发明实施例中，采用了泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线。因为电阻率曲线、声波时差曲线和自然伽马曲线都对有机质有较为敏感的响应。
在上述步骤105中，第一水动力带中的待测的泥页岩储层的有机碳含量可以通过如下方式计算：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1 公式一
在上述步骤105中，第二水动力带中的待测的泥页岩储层的有机碳含量可以通过如下方式计算：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2 公式二
在上述步骤105中，第三水动力带中的待测的泥页岩储层的有机碳含量可以通过如下方式计算：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3 公式三
上述步骤104中，确定第一水动力带中的各测试参数，可以通过如下方式实现：
根据第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4 公式四
上述步骤104中，确定第二水动力带中的各测试参数，可以通过如下方式实现：
根据第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5 公式五
上述步骤104中，确定第三水动力带中的各测试参数，可以通过如下方式实现：
根据第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6 公式六
值得说明的是，上述的TOC4、TOC5和TOC6是预先获取到的泥页岩岩样的有机碳含量，属于整个泥页岩段中的一个采样点的值。而上述步骤105中的TOC1、TOC2和TOC3是各水动力带对应的泥页岩段的有机碳含量。
通过上述步骤101至步骤105，应用到昭104井区，存在如下结果。
其中，对第一水动力带下块状泥岩拟合后得到如下计算公式：
TOC＝-7.65526+0.686675×logR-0.00357×AC+3.351002×logGR
如图2所示，得到的计算TOC(即通过公式计算得到的泥页岩的TOC值)和测试TOC(即泥页岩岩样的TOC值)具有很好的相关性，相关系数达到了0.763，可见该计算TOC结果可信度较高，十分接近原始储层有机碳含量。
对于第二水动力带下弱纹层状泥页岩段拟合后得到如下计算公式：
TOC＝-1.80981+2.022668×logR-0.00336×AC-0.2315×logGR
如图3所示，得到的计算TOC和测试TOC具有很好的相关性，相关系数达到了0.7348，可见该计算TOC结果可信度很高，十分接近原始储层有机碳含量。
对于第三水动力带下弱水动力带下纹层状泥页岩段拟合后得到如下计算公式：
TOC＝5.97815-4.73635×logR+0.013451×AC+2.00562×logGR
如图4所示，得到的计算TOC和测试TOC具有很好的相关性，相关系数达到了0.7638，可见该计算TOC结果可信度很高，十分接近原始储层有机碳含量。
将本发明实施例应用到一昭104井区，如图5所示，通过△logR方法和本发明方法得有机碳含量的比较结果，△logR方法的采样点的TOC与计算得到的TOC曲线较为离散，与TOC曲线重合度较差，而本发明实施例的方法的采样点的TOC与计算得到的TOC曲线重合度高，可见本发明检测的有机碳含量较为准确。
对应于图1的泥页岩储层有机碳含量的检测方法的实施例，本发明实施例提供了一种泥页岩储层有机碳含量的检测装置，如图6所示，包括：
泥页岩储层分段单元21，可以根据待测的泥页岩储层的岩心数据中的水动力数据，将待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带。
泥页岩岩样曲线值获取单元22，可以根据岩心数据获取第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值。
泥页岩岩样有机碳含量确定单元23，能够根据岩心数据分别确定第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量。
测试参数确定单元24，能够根据第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数。
泥页岩储层有机碳含量计算单元25，可以根据各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
另外，该泥页岩储层有机碳含量计算单元25，具体可以：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1。
另外，该泥页岩储层有机碳含量计算单元25，还可以：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2。
另外，该泥页岩储层有机碳含量计算单元25，还可以：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3。
此外，该测试参数确定单元24，具体可以：
根据第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4。
此外，该测试参数确定单元24，还可以：
根据第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5。
此外，该测试参数确定单元24，还可以：
根据第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6。
值得说明的是，本发明实施例提供的泥页岩储层有机碳含量的检测装置的具体实现方式可以参见图1中的方法实施例，此处不再赘述。
本发明实施例提供的泥页岩储层有机碳含量的检测装置，通过水动力数据将待测的泥页岩储层划分为多个水动力带区域，并根据测井数据获取各水动力带区域中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；并确定各水动力带区域中的泥页岩岩样的有机碳含量；在确定各水动力带区域的各测试参数之后，能够根据各水动力带区域待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，以及各测试参数，计算各水动力带区域中的待测的泥页岩储层的有机碳含量。可见本发明无需人为确定R基线和△t基线，且将一整段泥页岩根据水动力数据进行划分水动力带区域，进而检测各区域的有机碳含量，避免了现有技术中有机碳含量不准确，且整个检测过程繁琐缓慢的问题。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (14)
1.一种泥页岩储层有机碳含量的检测方法，其特征在于，包括：
将一整段泥页岩根据水动力数据进行划分水动力带区域：
根据待测的泥页岩储层的岩心数据中的水动力数据，将所述待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带；所述水动力数据包括水底扰动情况和含氧量；
根据待测的泥页岩储层的测井数据获取所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；
根据所述岩心数据分别确定所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量；
根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数；
根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
2.根据权利要求1所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述各测试参数、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算所述第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1。
3.根据权利要求1所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算所述第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2。
4.根据权利要求1所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量，包括：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算所述第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3。
5.根据权利要求2所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4。
6.根据权利要求3所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5。
7.根据权利要求4所述的泥页岩储层有机碳含量的检测方法，其特征在于，所述根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数，包括：
根据所述第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6。
8.一种泥页岩储层有机碳含量的检测装置，其特征在于，包括：
泥页岩储层分段单元，用于将一整段泥页岩根据水动力数据进行划分水动力带区域：
根据待测的泥页岩储层的岩心数据中的水动力数据，将所述待测的泥页岩储层分为第一水动力带、第二水动力带和第三水动力带；所述水动力数据包括水底扰动情况和含氧量；
泥页岩岩样曲线值获取单元，用于根据待测的泥页岩储层的测井数据获取所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值；
泥页岩岩样有机碳含量确定单元，用于根据所述岩心数据分别确定所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的有机碳含量；
测试参数确定单元，用于根据所述第一水动力带、第二水动力带以及第三水动力带中的泥页岩岩样的自然伽马曲线值、声波时差曲线值、电阻率曲线值和有机碳含量确定各测试参数；
泥页岩储层有机碳含量计算单元，用于根据所述各测试参数、第一水动力带、第二水动力带以及第三水动力带中的待测的泥页岩储层的自然伽马曲线、声波时差曲线、电阻率曲线，分别计算所述第一水动力带中的待测的泥页岩储层、第二水动力带中的待测的泥页岩储层、第三水动力带中的待测的泥页岩储层的有机碳含量。
9.根据权利要求8所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第一水动力带对应的各测试参数W1、X1、Y1、Z1、第一水动力带中的待测的泥页岩储层的自然伽马曲线值G1、声波时差曲线值A1、电阻率曲线值R1，计算所述第一水动力带中的待测的泥页岩储层的有机碳含量TOC1：
TOC1＝W1+X1·logR1-Y1·A1-Z1·logG1。
10.根据权利要求8所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第二水动力带对应的各测试参数W2、X2、Y2、Z2、第二水动力带中的待测的泥页岩储层的自然伽马曲线值G2、声波时差曲线值A2、电阻率曲线值R2，计算所述第二水动力带中的待测的泥页岩储层的有机碳含量TOC2：
TOC2＝W2+X2·logR2-Y2·A2-Z2·logG2。
11.根据权利要求8所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述泥页岩储层有机碳含量计算单元，具体用于：
根据第三水动力带对应的各测试参数W3、X3、Y3、Z3、第三水动力带中的待测的泥页岩储层的自然伽马曲线值G3、声波时差曲线值A3、电阻率曲线值R3，计算所述第三水动力带中的待测的泥页岩储层的有机碳含量TOC3：
TOC3＝W3+X3·logR3-Y3·A3-Z3·logG3。
12.根据权利要求9所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述测试参数确定单元，具体用于：
根据所述第一水动力带中的泥页岩岩样的自然伽马曲线值G4、声波时差曲线值A4、电阻率曲线值R4和有机碳含量TOC4进行多元拟合，确定第一水动力带对应的各测试参数W1、X1、Y1、Z1：
TOC4＝W1+X1·logR4-Y1·A4-Z1·logG4。
13.根据权利要求10所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述测试参数确定单元，具体用于：
根据所述第二水动力带中的泥页岩岩样的自然伽马曲线值G5、声波时差曲线值A5、电阻率曲线值R5和有机碳含量TOC5进行多元拟合，确定第二水动力带对应的各测试参数W2、X2、Y2、Z2：
TOC5＝W2+X2·logR5-Y2·A5-Z2·logG5。
14.根据权利要求11所述的泥页岩储层有机碳含量的检测装置，其特征在于，所述测试参数确定单元，具体用于：
根据所述第三水动力带中的泥页岩岩样的自然伽马曲线值G6、声波时差曲线值A6、电阻率曲线值R6和有机碳含量TOC6进行多元拟合，确定第三水动力带对应的各测试参数W3、X3、Y3、Z3：
TOC6＝W3+X3·logR6-Y3·A6-Z3·logG6。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410799475.1A CN104533405B (zh) | 2014-12-19 | 2014-12-19 | 泥页岩储层有机碳含量的检测方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410799475.1A CN104533405B (zh) | 2014-12-19 | 2014-12-19 | 泥页岩储层有机碳含量的检测方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104533405A CN104533405A (zh) | 2015-04-22 |
| CN104533405B true CN104533405B (zh) | 2018-02-02 |
Family
ID=52849005
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410799475.1A Active CN104533405B (zh) | 2014-12-19 | 2014-12-19 | 泥页岩储层有机碳含量的检测方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104533405B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104849427A (zh) * | 2015-04-27 | 2015-08-19 | 中国石油天然气股份有限公司 | 一种泥页岩层系有机碳含量的测定方法 |
| CN105653854A (zh) * | 2015-12-29 | 2016-06-08 | 中国石油天然气股份有限公司 | 计算碳酸盐岩烃源岩有机碳含量的方法和装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010118147A1 (en) * | 2009-04-07 | 2010-10-14 | Halliburton Energy Services, Inc. | Downhole carbon logging apparatus, systems, and methods |
| CN103114840A (zh) * | 2013-01-09 | 2013-05-22 | 中国石油天然气股份有限公司 | 一种高-过高成熟页岩有机碳含量计算方法及装置 |
| CN103670388A (zh) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种泥页岩有机碳含量的评价方法 |
| CN103670383A (zh) * | 2012-09-17 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种识别泥页岩油藏有效储层的方法及设备 |
| CN103792338A (zh) * | 2014-01-28 | 2014-05-14 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量的确定方法和装置 |
| CN103883322A (zh) * | 2014-04-16 | 2014-06-25 | 中国地质大学（北京） | 页岩气储层勘探方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7538319B2 (en) * | 2006-06-29 | 2009-05-26 | Baker Hughes Incorporated | Use of thorium-uranium ratio as an indicator of hydrocarbon source rock |
-
2014
- 2014-12-19 CN CN201410799475.1A patent/CN104533405B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010118147A1 (en) * | 2009-04-07 | 2010-10-14 | Halliburton Energy Services, Inc. | Downhole carbon logging apparatus, systems, and methods |
| CN103670383A (zh) * | 2012-09-17 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种识别泥页岩油藏有效储层的方法及设备 |
| CN103114840A (zh) * | 2013-01-09 | 2013-05-22 | 中国石油天然气股份有限公司 | 一种高-过高成熟页岩有机碳含量计算方法及装置 |
| CN103670388A (zh) * | 2013-12-12 | 2014-03-26 | 中国石油天然气股份有限公司 | 一种泥页岩有机碳含量的评价方法 |
| CN103792338A (zh) * | 2014-01-28 | 2014-05-14 | 中国石油天然气股份有限公司 | 一种烃源岩有机碳含量的确定方法和装置 |
| CN103883322A (zh) * | 2014-04-16 | 2014-06-25 | 中国地质大学（北京） | 页岩气储层勘探方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104533405A (zh) | 2015-04-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA3140364C (en) | Systems and methods for electromagnetic waterfront surveillance in a vicinity of an oil well | |
| CN103454678B (zh) | 一种地震切片等时性的确定方法及系统 | |
| US8645070B2 (en) | System and method for estimating fluid distribution in a subterranean reservoir | |
| US20150088424A1 (en) | Identifying geological formation depth structure using well log data | |
| US10422917B2 (en) | Sand-resistivity/conductivity and saturation evaluation in laminated formations with biaxial anisotropy | |
| US9938818B2 (en) | Hydraulic fracture permeability characterization from electromagnetic measurements | |
| CN107918156A (zh) | 检测海底节点采集地震数据极性的方法及装置 | |
| CN107144877B (zh) | 异常构造识别方法及装置 | |
| WO2016014485A1 (en) | Monitoring matrix acidizing operations | |
| US11008835B2 (en) | Horizontal reservoir description systems | |
| CN104533405B (zh) | 泥页岩储层有机碳含量的检测方法和装置 | |
| CN113550728A (zh) | 暂堵评价方法及装置 | |
| US11249218B2 (en) | Horizontal reservoir description systems | |
| Wang et al. | Dynamic monitoring of coalbed methane reservoirs using Super-Low Frequency electromagnetic prospecting | |
| CN108089228B (zh) | 一种确定地层岩石特性的解释数据方法和装置 | |
| CN110924937B (zh) | 套管井水淹层段的识别方法及装置 | |
| CN107179556B (zh) | 一种核磁共振的超前探测方法、系统及电子设备 | |
| CN111273354A (zh) | 地下开采引起的地下水变化分析方法和系统 | |
| CN105717016A (zh) | 一种煤样渗透率的计算方法 | |
| US20230060528A1 (en) | Reservoir fluid flow profiling in a wellbore environment | |
| CN109145375A (zh) | 裂缝渗透率模型校正方法及装置 | |
| CN112711076B (zh) | 提取石油钻井中泥浆侵入地层深度的方法和装置 | |
| CN118091786B (zh) | 基于井下测井资料的裂缝线密度校正方法和装置 | |
| CN105510986A (zh) | 一种确定海绿石砂岩中海绿石含量的方法和装置 | |
| Krainov et al. | PALEOCLIMATE RECORD IN THE BOTTOM SEDIMENTS OF LAKE BAIKAL FROM MAGNETIC SUSCEPTIBILITY DATA |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |