CN102052070B - 复杂岩性碎屑岩地层组分测井处理方法 - Google Patents

复杂岩性碎屑岩地层组分测井处理方法 Download PDF

Info

Publication number
CN102052070B
CN102052070B CN 201010585174 CN201010585174A CN102052070B CN 102052070 B CN102052070 B CN 102052070B CN 201010585174 CN201010585174 CN 201010585174 CN 201010585174 A CN201010585174 A CN 201010585174A CN 102052070 B CN102052070 B CN 102052070B
Authority
CN
China
Prior art keywords
coal
stratum
landwaste
calcareous
porosity
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.)
Expired - Fee Related
Application number
CN 201010585174
Other languages
English (en)
Other versions
CN102052070A (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.)
CNPC Chuanqing Drilling Engineering Co Ltd
Original Assignee
CNPC Chuanqing Drilling Engineering Co Ltd
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 CNPC Chuanqing Drilling Engineering Co Ltd filed Critical CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN 201010585174 priority Critical patent/CN102052070B/zh
Publication of CN102052070A publication Critical patent/CN102052070A/zh
Application granted granted Critical
Publication of CN102052070B publication Critical patent/CN102052070B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种复杂岩性碎屑岩地层组分测井处理方法,根据测井曲线特征判断地层是否含煤,根据地层是否含煤将地层分类分别进行处理;当地层含煤时,计算出长石石英、煤和孔隙度的含量;当地层不含煤时,应用混合岩性计算方法,计算出长石石英、钙质、岩屑和孔隙度的含量。本发明解决了含煤地层识别及含煤地层组分计算问题,解决了应用两条孔隙度测井曲线计算不含煤地层三种岩石组分及孔隙度的问题。

Description

复杂岩性碎屑岩地层组分测井处理方法
技术领域
本发明涉及一种复杂岩性碎屑岩地层组分测井处理方法,适用于岩石含有煤、长石石英、岩屑、钙质的复杂岩性碎屑岩地层的测井资料处理。
背景技术
目前碎屑岩地层组分处理方法,其骨架部分均认为是由单一的砂岩组成,没有计算岩石骨架具体成分的相对含量。例如,王宇文等于1994年在《中国海上油气(地质)》第1期第8卷中公开的“沉积盆地硅质碎屑岩地层模拟及应用实例”。
当地层存在有多种岩性的情况下,将导致储层参数计算结果不准确,同时影响与岩石含量有关的一系列评价的精度。出现上述问题的主要原因在于:①现有碎屑岩地层组分测井计算方法将地层骨架看成单一的,只能适用于岩性简单的地层;②当碎屑岩地层含有多种岩石成分时,现有算法无法计算其相对含量;③地层含有多种岩石成分时,计算的孔隙度精度差。
发明内容
本发明的目的在于克服现有碎屑岩地层组分处理方法存在的上述不足,提供一种复杂岩性碎屑岩地层组分测井处理方法,本发明解决了含煤地层识别及含煤地层组分计算问题,解决了应用两条孔隙度测井曲线计算不含煤地层三种岩石组分及孔隙度的问题。
为实现上述目的,本发明采用的技术方案如下:
一种复杂岩性碎屑岩地层组分测井处理方法,其特征在于:根据测井曲线特征判断地层是否含煤,根据地层是否含煤将地层分类分别进行处理;当地层含煤时,计算出长石石英、煤和孔隙度的含量;当地层不含煤时,应用混合岩性计算方法,计算出长石石英、钙质、岩屑和孔隙度的含量。
根据补偿声波AC、补偿中子CNL、密度DEN、深侧向电阻率RT测井曲线值是否符合含煤地层标准判别当前深度点地层是否含有煤,判别标准为:AC≥AC0(AC0为含煤地层补偿声波的下限值);RT≥RT0(RT0为含煤地层深侧向电阻率的下限值);CNL≥CNL0(CNL0为含煤地层补偿中子的下限值);DEN≤DEN0(DEN0为含煤地层密度的上限值),当以上测井曲线同时满足条件时,则判断该地层含煤。
当地层含煤时,则岩石骨架由长石石英砂岩和煤组成,通过联立补偿声波、补偿中子及物质平衡方程建立方程组:
                                                 
Figure 469616DEST_PATH_IMAGE001
求解出岩石骨架成分及孔隙度;
方程组中,V1为长石石英砂岩的含量,V2为煤的含量,
Figure 386757DEST_PATH_IMAGE002
为孔隙度,Tma1、Tma2、TF分别为长石石英砂岩、煤和流体(水)的补偿声波时差,CNL1、CNL2、NF分别为长石石英砂岩、煤和流体的补偿中子。
当地层不含煤时,则岩石骨架由长石石英砂岩、钙质、岩屑组成,与孔隙流体构成两个三角形,在补偿中子~补偿声波交会图上,从补偿中子与补偿声波的位置,作到钙质与岩屑连线的垂线,如果垂足位于钙质与岩屑的连线中间,则混合骨架成分由垂足到两种岩性点的距离进行计算;如果垂足位于钙质的左边,说明没有岩屑,混合骨架等于钙质的骨架;如果垂足位于岩屑的右边,说明没有钙质,混合骨架等于岩屑的骨架。
所述三角形中,设垂足到钙质点的距离为L1,到岩屑点的距离为L2,通过计算式
Figure 228811DEST_PATH_IMAGE003
得到混合骨架的补偿中子和补偿声波值,式中,
Figure 624020DEST_PATH_IMAGE004
为混合岩性的补偿中子值,
Figure 251441DEST_PATH_IMAGE005
为钙质的补偿中子值,为岩屑的补偿中子值,
Figure 403254DEST_PATH_IMAGE007
为混合岩性的补偿声波值, 
Figure 602154DEST_PATH_IMAGE008
为钙质的补偿声波值,
Figure 271033DEST_PATH_IMAGE009
为岩屑的补偿声波值;
用长石石英砂岩、混合骨架及水组成的岩性三角形可得到长石石英砂岩、混合矿物及孔隙度的含量,然后对混合矿物成分V
Figure 77446DEST_PATH_IMAGE010
进行分配,得出钙质和岩屑的含量;
通过以上处理,得出长石石英砂岩、钙质、岩屑及孔隙度各组分的大小。
采用本发明的优点在于:
一、本发明采用多条曲线门槛值作为是否含煤判别标准,有利于提高含煤地层的识别的准确率。
二、本发明对于含煤地层,将地层简化为由长石石英、煤和孔隙组成,所以提高了含煤地层组分计算的准确性。
三、本发明对于非含煤地层,采用了混合岩性计算方法,应用两条孔隙度曲线就能得到三种岩石骨架组分,提高了岩石成分处理的准确性。
四、本发明得到了合理的岩石成分,可提高测井计算孔隙度的准确性,通过本发明可得出碎屑岩地层的主要岩石成分,从而提高地层组分及孔隙度处理的准确性,易于实现。
五、本发明可识别出含煤与非含煤地层,并可计算复杂岩性碎屑岩地层的长石石英、钙质、岩屑、煤的含量,有利于提高测井计算孔隙度等储层参数的精度。
附图说明
图1为本发明当地层不含煤时岩石骨架与孔隙流体构成两个三角形的示意图
图2为本发明在补偿中子~补偿声波交会图上的垂线示意图。
具体实施方式
一种复杂岩性碎屑岩地层组分测井处理方法,根据测井曲线特征判断地层是否含煤,根据地层是否含煤将地层分类分别进行处理;当地层含煤时,计算出长石石英、煤和孔隙度的含量;当地层不含煤时,应用混合岩性计算方法,计算出长石石英、钙质、岩屑和孔隙度的含量。
具体处理过程如下:
根据补偿声波AC、补偿中子CNL、密度DEN、深侧向电阻率RT测井曲线值是否符合含煤地层标准判别当前深度点地层是否含有煤。判别标准为:AC≥AC0(AC0为含煤地层补偿声波的下限);RT≥RT0(RT0为含煤地层深侧向电阻率的下限);CNL≥CNL0(CNL0为含煤地层补偿中子的下限);DEN≤DEN0(DEN0为含煤地层密度的上限)。当以上4条测井曲线同时满足条件时,则判断该地层含煤。
当地层含煤时,则岩石骨架由长石石英砂岩和煤组成,通过联立补偿声波、补偿中子及物质平衡方程建立的方程组就可求解出岩石骨架成分及孔隙度。
 
Figure 831775DEST_PATH_IMAGE001
上式中,V1为长石石英砂岩的含量,V2为煤的含量,
Figure 631104DEST_PATH_IMAGE002
为孔隙度,Tma1、Tma2、TF分别为长石石英砂岩、煤和流体(水)的补偿声波时差,CNL1、CNL2、NF分别为长石石英砂岩、煤和流体的补偿中子。
当地层不含煤时,则岩石骨架由长石石英砂岩、钙质、岩屑组成,与孔隙流体构成两个三角形,如图1所示。
钙质与岩屑处于同一个岩性三角形,如果测井资料点落入该三角形内,则只能求出钙质与岩屑,没有长石石英组分,这与地质特征不符。另外钙质与岩屑所处的岩性三角形非常小,当岩屑的中子或声波时差落入钙质线时,则无法构成三角形。鉴于上述原因,将钙质与岩屑以一种骨架参与岩性三角形的处理。具体处理方法为:
在补偿中子~补偿声波交会图上,由补偿中子与补偿声波的位置,作到钙质与岩屑连线的垂线,如果垂足位于钙质与岩屑的连线中间,则混合骨架成分由垂足到两种岩性点的距离进行计算;如果垂足位于钙质的左边,说明没有岩屑,混合骨架等于钙质的骨架;如果垂足位于岩屑的右边,说明没有钙质,混合骨架等于岩屑的骨架。如图2所示。
设垂足到钙质点的距离为L1,到岩屑点的距离为L2,则混合骨架的补偿中子和补偿声波值可用下式计算:
 
Figure 154489DEST_PATH_IMAGE003
上式中,
Figure 584334DEST_PATH_IMAGE004
为混合岩性的补偿中子值,
Figure 639009DEST_PATH_IMAGE005
为钙质的补偿中子值,
Figure 914132DEST_PATH_IMAGE006
为岩屑的补偿中子值,
Figure 619920DEST_PATH_IMAGE007
为混合岩性的补偿声波值, 
Figure 955086DEST_PATH_IMAGE008
为钙质的补偿声波值,
Figure 497057DEST_PATH_IMAGE009
为岩屑的补偿声波值;
这样,用长石石英砂岩、混合骨架及水组成的岩性三角形可求解长石石英砂岩、混合矿物及孔隙度的含量,然后对混合矿物成分V用下式进行分配,就可得出钙质和岩屑的含量。
 
Figure 575872DEST_PATH_IMAGE010
通过以上求解,就可得出长石石英砂岩、钙质、岩屑及孔隙度各组分的大小。
设置长石石英含量为V1,钙质含量为V2,岩屑含量为V3,煤含量为V4,孔隙度为。对测量井段进行逐点处理,就可得到连续的地质剖面。

Claims (1)

1.一种复杂岩性碎屑岩地层组分测井处理方法,其特征在于:根据测井曲线特征判断地层是否含煤,根据地层是否含煤将地层分类分别进行处理;当地层含煤时,计算出长石石英、煤和孔隙度的含量;当地层不含煤时,应用混合岩性计算方法,计算出长石石英、钙质、岩屑和孔隙度的含量;
根据补偿声波AC、补偿中子CNL、密度DEN、深侧向电阻率RT测井曲线值是否符合含煤地层标准判别当前深度点地层是否含有煤,判别标准为:AC≥AC0,AC0为含煤地层补偿声波的下限值;RT≥RT0,RT0为含煤地层深侧向电阻率的下限值;CNL≥CNL0,CNL0为含煤地层补偿中子的下限值;DEN≤DEN0,DEN0为含煤地层密度的上限值,当以上测井曲线同时满足条件时,则判断该地层含煤;
当地层含煤时,则岩石骨架由长石石英砂岩和煤组成,通过联立补偿声波、补偿中子及物质平衡方程建立方程组:
Figure DEST_PATH_IMAGE002
求解出岩石骨架成分及孔隙度;
方程组中,
Figure DEST_PATH_IMAGE004
为长石石英砂岩的含量,为煤的含量,
Figure DEST_PATH_IMAGE008
为孔隙度,
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE014
分别为长石石英砂岩、煤和流体的补偿声波时差,
Figure DEST_PATH_IMAGE016
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE020
分别为长石石英砂岩、煤和流体的补偿中子;
当地层不含煤时,则岩石骨架由长石石英砂岩、钙质、岩屑组成,与孔隙流体构成两个三角形,在补偿中子~补偿声波交会图上,从补偿中子与补偿声波的位置,作到钙质与岩屑连线的垂线,如果垂足位于钙质与岩屑的连线中间,则混合骨架成分由垂足到两种岩性点的距离进行计算;如果垂足位于钙质的左边,说明没有岩屑,混合骨架等于钙质的骨架;如果垂足位于岩屑的右边,说明没有钙质,混合骨架等于岩屑的骨架。
CN 201010585174 2010-12-13 2010-12-13 复杂岩性碎屑岩地层组分测井处理方法 Expired - Fee Related CN102052070B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010585174 CN102052070B (zh) 2010-12-13 2010-12-13 复杂岩性碎屑岩地层组分测井处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010585174 CN102052070B (zh) 2010-12-13 2010-12-13 复杂岩性碎屑岩地层组分测井处理方法

Publications (2)

Publication Number Publication Date
CN102052070A CN102052070A (zh) 2011-05-11
CN102052070B true CN102052070B (zh) 2013-06-05

Family

ID=43956828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010585174 Expired - Fee Related CN102052070B (zh) 2010-12-13 2010-12-13 复杂岩性碎屑岩地层组分测井处理方法

Country Status (1)

Country Link
CN (1) CN102052070B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102508310B (zh) * 2011-10-31 2013-08-28 中国矿业大学 一种煤田火区上部地层孔隙度分布的探测方法
CN103122762B (zh) * 2013-01-30 2015-07-08 中国石油天然气股份有限公司 一种非常规泥页岩油气藏有效压裂层段的检测方法及装置
CN103775072A (zh) * 2014-01-16 2014-05-07 燕山大学 基于测井资料的煤岩类型确定方法
CN103867197B (zh) * 2014-04-04 2016-07-20 中国石油集团川庆钻探工程有限公司 复杂岩性天然气层声波时差判别法
CN105298481B (zh) * 2015-10-15 2016-06-15 山东科技大学 一种深部煤系及薄煤层多级叠加测井识别方法
CN105756673B (zh) * 2016-02-26 2017-04-05 山东科技大学 一种基于测井曲线的煤与天然焦快速识别方法
CN107795320B (zh) * 2016-09-05 2021-01-05 中国石油化工股份有限公司 一种水平井碳酸盐岩储层参数的计算方法
CN108562950B (zh) * 2017-12-11 2020-01-03 中国石油天然气集团公司 一种基于测井信息对地层层位进行智能化划分的方法
CN109061752B (zh) * 2018-06-26 2020-01-17 西南石油大学 一种含灰质地层的电阻率曲线校正方法
CN109339773A (zh) * 2018-10-09 2019-02-15 中国地质大学(北京) 基于砾岩基质含量测井孔隙度模型及其构建方法和应用
CN111241460B (zh) * 2020-01-19 2021-08-24 西安石油大学 一种复杂致密储层孔隙度计算方法
CN113720745A (zh) * 2020-05-26 2021-11-30 中国石油化工股份有限公司 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766284A (zh) * 2005-11-29 2006-05-03 大庆石油管理局 一种通过异常流体浸入井筒获取地层压力的识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2442014B (en) * 2006-09-19 2011-08-03 Reeves Wireline Tech Ltd A method of acquiring and processing neutron log data, and related apparatuses

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766284A (zh) * 2005-11-29 2006-05-03 大庆石油管理局 一种通过异常流体浸入井筒获取地层压力的识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘芳.浅谈数字测井技术在靖宇县煤田地质勘探中的应用.《中国水运》.2010,第10卷(第5期),
江苏地区复杂岩性储层测井评价方法研究;马静等;《内蒙古石油化工》;20050425(第2期);98 *
浅谈数字测井技术在靖宇县煤田地质勘探中的应用;刘芳;《中国水运》;20100515;第10卷(第5期);149 *
马静等.江苏地区复杂岩性储层测井评价方法研究.《内蒙古石油化工》.2005,(第2期),

Also Published As

Publication number Publication date
CN102052070A (zh) 2011-05-11

Similar Documents

Publication Publication Date Title
CN102052070B (zh) 复杂岩性碎屑岩地层组分测井处理方法
CN104453873B (zh) 页岩油气经济有效层段的评价方法
CN103122762B (zh) 一种非常规泥页岩油气藏有效压裂层段的检测方法及装置
CN104636588B (zh) 计算烃源岩中有机碳含量的方法及装置
CN104698497B (zh) 一种砂泥岩薄互层层位的精细标定方法
CN103867197B (zh) 复杂岩性天然气层声波时差判别法
CN103670383B (zh) 一种识别泥页岩油藏有效储层的方法及设备
CN104806232B (zh) 一种确定碳酸盐岩储层孔隙度下限的方法
CN103617345A (zh) 用元素含量计算地层岩石矿物成分的方法
CN103670388A (zh) 一种泥页岩有机碳含量的评价方法
CN106528707A (zh) 一种煤层顶板砂岩含水层富水性评价方法
CN103867198B (zh) 碳酸盐岩天然气层地层密度判别法
CN104632201B (zh) 一种页岩油气“甜点区”有机碳含量下限值的确定方法
CN109870720B (zh) 一种页岩气微裂缝测井识别方法
CN107576982A (zh) 一种砂岩型铀矿地震综合预测方法
CN105676294B (zh) 一种碳酸盐岩缝洞型储集体定量描述方法
CN105840189A (zh) 一种页岩气岩石物理与测井方法
CN102116871A (zh) 一种提高地震反演数据预测储层能力的方法
CN101315428B (zh) 基于横波速度的高孔隙碎屑岩的高分辨率勘探方法
Sun et al. Investigation of pore-type heterogeneity and its inherent genetic mechanisms in deeply buried carbonate reservoirs based on some analytical methods of rock physics
CN105093304A (zh) 一种地球物理勘探利用测井曲线自动求取岩性曲线的方法
CN102536223B (zh) 煤层裂隙孔隙度快速计算方法
CN109458176A (zh) 碳酸盐岩储层压力的预测方法及其应用
CN108345046B (zh) 碎屑岩储层埋藏过程中成岩相划分方法
KR20180014294A (ko) 물리검층 자료를 이용한 s파 속도 산출방법

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130605

Termination date: 20171213

CF01 Termination of patent right due to non-payment of annual fee