CN102052070B - 复杂岩性碎屑岩地层组分测井处理方法 - Google Patents
复杂岩性碎屑岩地层组分测井处理方法 Download PDFInfo
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 39
- 238000003672 processing method Methods 0.000 title claims description 10
- 239000003245 coal Substances 0.000 claims abstract description 76
- 239000010433 feldspar Substances 0.000 claims abstract description 34
- 239000010453 quartz Substances 0.000 claims abstract description 32
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims description 24
- 239000000203 mixture Substances 0.000 claims description 19
- 239000012530 fluid Substances 0.000 claims description 10
- 239000000470 constituent Substances 0.000 claims description 7
- 238000000205 computational method Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000000034 method Methods 0.000 abstract description 5
- 239000000463 material Substances 0.000 abstract 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 101100058498 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) CNL1 gene Proteins 0.000 description 2
- 101100294209 Schizosaccharomyces pombe (strain 972 / ATCC 24843) cnl2 gene Proteins 0.000 description 2
- 101100401683 Schizosaccharomyces pombe (strain 972 / ATCC 24843) mis13 gene Proteins 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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为含煤地层密度的上限值），当以上测井曲线同时满足条件时，则判断该地层含煤。
当地层含煤时，则岩石骨架由长石石英砂岩和煤组成，通过联立补偿声波、补偿中子及物质平衡方程建立方程组：
求解出岩石骨架成分及孔隙度；
方程组中，V1为长石石英砂岩的含量，V2为煤的含量，
为孔隙度，Tma1、Tma2、TF分别为长石石英砂岩、煤和流体(水)的补偿声波时差，CNL1、CNL2、NF分别为长石石英砂岩、煤和流体的补偿中子。
当地层不含煤时，则岩石骨架由长石石英砂岩、钙质、岩屑组成，与孔隙流体构成两个三角形，在补偿中子～补偿声波交会图上，从补偿中子与补偿声波的位置，作到钙质与岩屑连线的垂线，如果垂足位于钙质与岩屑的连线中间，则混合骨架成分由垂足到两种岩性点的距离进行计算；如果垂足位于钙质的左边，说明没有岩屑，混合骨架等于钙质的骨架；如果垂足位于岩屑的右边，说明没有钙质，混合骨架等于岩屑的骨架。
所述三角形中，设垂足到钙质点的距离为L1,到岩屑点的距离为L2，通过计算式
得到混合骨架的补偿中子和补偿声波值，式中，
为混合岩性的补偿中子值，
为钙质的补偿中子值，为岩屑的补偿中子值，
为混合岩性的补偿声波值， 
为钙质的补偿声波值，
为岩屑的补偿声波值；
用长石石英砂岩、混合骨架及水组成的岩性三角形可得到长石石英砂岩、混合矿物及孔隙度的含量，然后对混合矿物成分V混用
进行分配，得出钙质和岩屑的含量；
通过以上处理，得出长石石英砂岩、钙质、岩屑及孔隙度各组分的大小。
采用本发明的优点在于：
一、本发明采用多条曲线门槛值作为是否含煤判别标准，有利于提高含煤地层的识别的准确率。
二、本发明对于含煤地层，将地层简化为由长石石英、煤和孔隙组成，所以提高了含煤地层组分计算的准确性。
三、本发明对于非含煤地层，采用了混合岩性计算方法，应用两条孔隙度曲线就能得到三种岩石骨架组分，提高了岩石成分处理的准确性。
四、本发明得到了合理的岩石成分，可提高测井计算孔隙度的准确性，通过本发明可得出碎屑岩地层的主要岩石成分，从而提高地层组分及孔隙度处理的准确性，易于实现。
五、本发明可识别出含煤与非含煤地层，并可计算复杂岩性碎屑岩地层的长石石英、钙质、岩屑、煤的含量，有利于提高测井计算孔隙度等储层参数的精度。
附图说明
图1为本发明当地层不含煤时岩石骨架与孔隙流体构成两个三角形的示意图
图2为本发明在补偿中子～补偿声波交会图上的垂线示意图。
具体实施方式
一种复杂岩性碎屑岩地层组分测井处理方法，根据测井曲线特征判断地层是否含煤，根据地层是否含煤将地层分类分别进行处理；当地层含煤时，计算出长石石英、煤和孔隙度的含量；当地层不含煤时，应用混合岩性计算方法，计算出长石石英、钙质、岩屑和孔隙度的含量。
具体处理过程如下：
根据补偿声波AC、补偿中子CNL、密度DEN、深侧向电阻率RT测井曲线值是否符合含煤地层标准判别当前深度点地层是否含有煤。判别标准为：AC≥AC0（AC0为含煤地层补偿声波的下限）；RT≥RT0（RT0为含煤地层深侧向电阻率的下限）；CNL≥CNL0（CNL0为含煤地层补偿中子的下限）；DEN≤DEN0（DEN0为含煤地层密度的上限）。当以上4条测井曲线同时满足条件时，则判断该地层含煤。
当地层含煤时，则岩石骨架由长石石英砂岩和煤组成，通过联立补偿声波、补偿中子及物质平衡方程建立的方程组就可求解出岩石骨架成分及孔隙度。
上式中，V1为长石石英砂岩的含量，V2为煤的含量，
为孔隙度，Tma1、Tma2、TF分别为长石石英砂岩、煤和流体(水)的补偿声波时差，CNL1、CNL2、NF分别为长石石英砂岩、煤和流体的补偿中子。
当地层不含煤时，则岩石骨架由长石石英砂岩、钙质、岩屑组成，与孔隙流体构成两个三角形，如图1所示。
钙质与岩屑处于同一个岩性三角形，如果测井资料点落入该三角形内，则只能求出钙质与岩屑，没有长石石英组分，这与地质特征不符。另外钙质与岩屑所处的岩性三角形非常小，当岩屑的中子或声波时差落入钙质线时，则无法构成三角形。鉴于上述原因，将钙质与岩屑以一种骨架参与岩性三角形的处理。具体处理方法为：
在补偿中子～补偿声波交会图上，由补偿中子与补偿声波的位置，作到钙质与岩屑连线的垂线，如果垂足位于钙质与岩屑的连线中间，则混合骨架成分由垂足到两种岩性点的距离进行计算；如果垂足位于钙质的左边，说明没有岩屑，混合骨架等于钙质的骨架；如果垂足位于岩屑的右边，说明没有钙质，混合骨架等于岩屑的骨架。如图2所示。
设垂足到钙质点的距离为L1,到岩屑点的距离为L2，则混合骨架的补偿中子和补偿声波值可用下式计算：
上式中，
为混合岩性的补偿中子值，
为钙质的补偿中子值，
为岩屑的补偿中子值，
为混合岩性的补偿声波值， 
为钙质的补偿声波值，
为岩屑的补偿声波值；
这样，用长石石英砂岩、混合骨架及水组成的岩性三角形可求解长石石英砂岩、混合矿物及孔隙度的含量，然后对混合矿物成分V混用下式进行分配，就可得出钙质和岩屑的含量。
通过以上求解，就可得出长石石英砂岩、钙质、岩屑及孔隙度各组分的大小。
设置长石石英含量为V1，钙质含量为V2，岩屑含量为V3，煤含量为V4，孔隙度为。对测量井段进行逐点处理，就可得到连续的地质剖面。
Claims (1)
1.一种复杂岩性碎屑岩地层组分测井处理方法，其特征在于：根据测井曲线特征判断地层是否含煤，根据地层是否含煤将地层分类分别进行处理；当地层含煤时，计算出长石石英、煤和孔隙度的含量；当地层不含煤时，应用混合岩性计算方法，计算出长石石英、钙质、岩屑和孔隙度的含量；
根据补偿声波AC、补偿中子CNL、密度DEN、深侧向电阻率RT测井曲线值是否符合含煤地层标准判别当前深度点地层是否含有煤，判别标准为：AC≥AC0，AC0为含煤地层补偿声波的下限值；RT≥RT0，RT0为含煤地层深侧向电阻率的下限值；CNL≥CNL0，CNL0为含煤地层补偿中子的下限值；DEN≤DEN0，DEN0为含煤地层密度的上限值，当以上测井曲线同时满足条件时，则判断该地层含煤；
当地层含煤时，则岩石骨架由长石石英砂岩和煤组成，通过联立补偿声波、补偿中子及物质平衡方程建立方程组：
求解出岩石骨架成分及孔隙度；
方程组中，
为长石石英砂岩的含量，为煤的含量，
为孔隙度，
、
、
分别为长石石英砂岩、煤和流体的补偿声波时差，
、
、
分别为长石石英砂岩、煤和流体的补偿中子；
当地层不含煤时，则岩石骨架由长石石英砂岩、钙质、岩屑组成，与孔隙流体构成两个三角形，在补偿中子～补偿声波交会图上，从补偿中子与补偿声波的位置，作到钙质与岩屑连线的垂线，如果垂足位于钙质与岩屑的连线中间，则混合骨架成分由垂足到两种岩性点的距离进行计算；如果垂足位于钙质的左边，说明没有岩屑，混合骨架等于钙质的骨架；如果垂足位于岩屑的右边，说明没有钙质，混合骨架等于岩屑的骨架。
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)
| 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 | 西安石油大学 | 一种复杂致密储层孔隙度计算方法 |
| CN113720745B (zh) * | 2020-05-26 | 2024-07-16 | 中国石油化工股份有限公司 | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1766284A (zh) * | 2005-11-29 | 2006-05-03 | 大庆石油管理局 | 一种通过异常流体浸入井筒获取地层压力的识别方法 |
Family Cites Families (1)
| 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 |
-
2010
- 2010-12-13 CN CN 201010585174 patent/CN102052070B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1766284A (zh) * | 2005-11-29 | 2006-05-03 | 大庆石油管理局 | 一种通过异常流体浸入井筒获取地层压力的识别方法 |
Non-Patent Citations (4)
| 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) | 复杂岩性碎屑岩地层组分测井处理方法 | |
| CN104516017B (zh) | 一种碳酸盐岩岩石物理参数地震反演方法 | |
| CN104698497B (zh) | 一种砂泥岩薄互层层位的精细标定方法 | |
| CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
| Bouma et al. | Fine-Grained Turbidite Systems: AAPG Memoir 72 | |
| CN103617345A (zh) | 用元素含量计算地层岩石矿物成分的方法 | |
| CN103670383A (zh) | 一种识别泥页岩油藏有效储层的方法及设备 | |
| CN103670388A (zh) | 一种泥页岩有机碳含量的评价方法 | |
| CN104675392B (zh) | 一种基于叠前多参数降维的储层岩性识别方法 | |
| CN106600436A (zh) | 计算页岩气地层矿物组分含量及孔隙度的方法 | |
| CN103867198B (zh) | 碳酸盐岩天然气层地层密度判别法 | |
| CN106528707A (zh) | 一种煤层顶板砂岩含水层富水性评价方法 | |
| CN107576982A (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) | 碎屑岩储层埋藏过程中成岩相划分方法 | |
| CN112784404B (zh) | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 | |
| CN104698494A (zh) | 一种计算异常地层压力的方法 | |
| CN104808259A (zh) | 一种膏岩分布的检测方法及装置 | |
| Sarkar et al. | Seismic-petrophysical reservoir characterization in the northern part of the Chicontepec Basin, Mexico | |
| CN105093330B (zh) | 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法 | |
| Bazzell | Origin of brecciated intervals and petrophysical analyses; the three forks formation, williston basin, North Dakota, USA |
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 |