CN104101649B - 利用声速判定页岩中的有机质含量的方法 - Google Patents

利用声速判定页岩中的有机质含量的方法 Download PDF

Info

Publication number
CN104101649B
CN104101649B CN201410270856.0A CN201410270856A CN104101649B CN 104101649 B CN104101649 B CN 104101649B CN 201410270856 A CN201410270856 A CN 201410270856A CN 104101649 B CN104101649 B CN 104101649B
Authority
CN
China
Prior art keywords
shale
content
velocity
sound
organic
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
CN201410270856.0A
Other languages
English (en)
Other versions
CN104101649A (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.)
Nanjing University
Sinopec Geophysical Research Institute
Original Assignee
Nanjing University
Sinopec Geophysical Research Institute
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 Nanjing University, Sinopec Geophysical Research Institute filed Critical Nanjing University
Priority to CN201410270856.0A priority Critical patent/CN104101649B/zh
Publication of CN104101649A publication Critical patent/CN104101649A/zh
Application granted granted Critical
Publication of CN104101649B publication Critical patent/CN104101649B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

利用声速判定页岩中的有机质含量的方法,利用含有有机质的页岩中的声速与有机质含量的变化的特性判断页岩中有机质的含量,页岩在有机质含量变化后纵波和横波声速发生明显的变化;采用超声的方法无损检测页岩的有机质的含量,即采用测量页岩的声速随有机质含量的变化来检测页岩的有机质的含量;通过发射和接收超声换能器,测量声波在页岩样品的传播时间和样品的长度,得到页岩的声速,根据声速与页岩中有机质含量的关系来判定页岩中有机质的含量。

Description

利用声速判定页岩中的有机质含量的方法
一、技术领域
本发明涉及利用超声波进行页岩中有机质含量的判定的方法,尤其是对有机质含量较高的情况。
二、背景技术
页岩是一种沉积岩,成分复杂,但都具有薄页状或薄片层状的节理,主要是由黏土沉积经压力和温度形成的岩石,但其中混杂有石英、长石的碎屑以及其他化学物质。页岩通常会表现出各向异性,这种性质来源于沉积之前、过程中、和之后的地质过程的作用。页岩各向异性的物理模型需要多尺度的结构信息,包含粘土域中晶体和细小孔隙的结构,这些粘土域的队列去向,以及包含泥沙和富含有机质的粘土的细小分层。针对富含油母岩质的页岩,如何有效地判定其有机质的含量目前还没有很好的方法。
声速是声学的一个基本参量,它可以反映材料的弹性特性,有望反映页岩的有机质含量的变化。
三、发明内容
本发明的目的是:采用超声的方法来研究声速与有机质(油母质)含量的变化关系,即采用测量页岩的声速来判定页岩中的有机质含量的变化。
本发明的技术方案:利用声速判定页岩中的有机质含量的方法,利用含有有机质的页岩中的声速与有机质含量的变化的特性判断页岩中有机质的含量,页岩在有机质含量变化后纵波和横波声速发生明显的变化;采用超声的方法无损检测页岩的有机质的含量,即采用测量页岩的声速随有机质含量的变化来检测页岩的有机质的含量;通过发射和接收超声换能器,测量声波在页岩样品的传播时间和样品的长度,得到页岩中的料的声速,根据声速与页岩中有机质含量的关系来判定页岩中有机质的含量:
1)由不同矿物组成无油的自由“页岩”即干页岩,不同矿物被认为是各向同性的或各向异性的。不同矿物和它们各自的体积百分比,从X-射线衍射(XRD)分析得到的;粘土,石英,方解石等页岩的无机成分矿物存在于岩石的弹性常数使用Reuss-Voigt-Hill平均计算;
2)考虑油母质为背景材料,使用带自洽的各向异性微分有效介质模型添加“干页岩”包裹到油母质背景材料中,形成油母质和干页岩复合材料;
3)添加气孔进入油母质和干页岩复合材料,使用带自洽的各向异性微分有效介质模型重新估计带气孔油母质和干页岩的弹性模量;
4)添加液体进入气孔中,使用Brown-Korringa模型,获得流体饱和油母质和干页岩的弹性模量;根据上述模型计算得到饱和流体页岩的弹性模量随油母质含量的变化图;
5)根据弹性模量和声速的关系得到流体饱和油母质和干页岩中的纵波和横波的声速;得到饱和流体页岩的声速随油母质含量的变化曲线。
本发明尤其是对于Bazhenov页岩,可由图2根据测量到的声速来判定页岩中有机质的含量。本发明方法可推广到其他类型的页岩。
我们利用带自洽的各向异性微分有效介质模型得到页岩的弹性模量随有机质含量的变化关系,根据弹性模量与声速的关系可得到页岩中的声速,从而获得声速和有机质(油母质)含量的变化曲线,利用此变化曲线,可以通过测量声速来获得页岩中有机质(油母质)的含量。
本发明的有益效果:与现有技术相比具有的特点是:通过声学方法来测量页岩中声速随有机质含量的变化,寻求页岩声速随有机质含量之间的关系。此方法的特点是无损和实时,方法简单易行。这里主要针对有机质含量较高的页岩(有机质的体积比大于10%),因为有机质含量低,此模型可能失效。
发展页岩的有机质含量的检测,是一项具有重大现实意义的课题。因此目前页岩的研究已成为石油工业的热点课题,由于页岩材料的有机质含量的高低直接影响到页岩的开采利用的价值,因此准确和实时判定页岩有机质的含量显得特别的重要。利用声速判定其相变的方法。本发明以Bazhenov页岩为列来说明采用声速来判定有机质含量的方法。利用页岩的声速与有机质含量变化的敏感性,页岩在不同有机质含量后声速发生明显的变化;采用超声的方法无损检测页岩的声速,即根据页岩的声速和有机质的含量的关系来检测页岩的有机质的含量。
四、附图说明
图1饱和流体页岩的弹性模量随油母质含量的变化;
图2饱和流体页岩的声速随油母质含量的变化。
五、具体实现方案
本发明方法包括测量和计算,过程分为以下几步:
(1)矿物(粘土,石英,方解石等页岩的无机成分)存在于岩石中,干页岩体积模量和剪切模量使用Reuss-Voigt-Hill平均计算。矿物组成无油母质的自由“页岩”即干页岩,它可以被认为是各向同性的或各向异性的。矿物质和它们各自的体积百分比,通常可从X-射线衍射(XRD)分析得到的。
(2)考虑油母质为背景材料,使用带自洽的各向异性微分有效介质模型添加“干页岩”包裹到背景材料中,形成油母质和干页岩复合材料,
(3)添加气孔进入油母质和干页岩复合材料,使用带自洽的各向异性微分有效介质模型重新估计带气孔的油母质和干页岩的弹性模量。
(4)添加液体进入气孔中,使用Brown-Korringa模型,以获得流体饱和油母质和干页岩的弹性模量。根据上述模型计算得到饱和流体页岩的弹性模量随油母质含量的变化图;
(5)根据弹性模量和声速的关系得到流体饱和油母质和干页岩中的纵波和横波的声速;
得到饱和油母质页岩中的声速随油母质含量的变化曲线。
例如:根据Bazhenov页岩的参数:
由Reuss-Voigt-Hill模型知,
M v = Σ i = 1 N φ i M i
其中Mi是Voigt有效弹性模量,它可以是体积模量Kv或者剪切模量μV,φi和Mi是分别是第i种成分(粘土,石英,方解石等页岩的无机成分)的含量和弹性模量。
1 M R = Σ i = 1 N φ i M i
M vrh = M v + M r 2
Mv,MR和Mvrh的定义与Mi相同,可以是体积模量或剪切模量。
代入数据计算得:
干页岩体积模量和剪切模量为:K=32.08GPa,μ=23.92GPa.
相应地,干页岩的各向异性微分有效介质模型中的有关参数(弹性模量):
C33=K+4μ/3=63.97GPa,C44=μ=23.92GPa,C12=K-2μ/3=16.13GPa.
油母质的各向异性微分有效介质模型中的有关参数(弹性模量):
C11=C33=6.5GPa,C44=C66=2.7GPa,C13=C12=C33-2C44=1.10GPa.
取油母质的体积分数为16.8%,干页岩的纵横比(纵波比横波声速a)为0.1,密度为g=2.34g/cm3,油母质-干页岩的混合物的弹性模量为:
C11=52.74GPa,C33=46.61GPa,C44=21.30GPa,C66=23.61GPa,C13=11.72GPa,
页岩可以看成是横向各向同性介质,横向各向同性介质的弹性张量由5个独立的弹性常数描述,用矩阵表示为:
C = C 11 C 12 C 13 0 0 0 C 12 C 11 C 13 0 0 0 C 13 C 13 C 33 0 0 0 0 0 0 C 44 0 0 0 0 0 0 C 44 0 0 0 0 0 0 C 66
其中 C 66 = 1 2 ( C 11 - C 12 )
取孔隙率4.12%,得到带气孔的页岩的弹性模量:
C11=50.03GPa,C33=42.17GPa,C44=20.89GPa,C66=23.58GPa,C13=10.74GPa
用Brown-Korringa模型,得到饱和流体的页岩的弹性模量:
C11=50.45GPa,C33=42.60GPa,C44=21.30GPa,C66=23.61GPa,C13=11.16GPa
可根据上述模型计算得到饱和流体页岩的弹性模量随油母质含量的变化图(图1)根据(现有的公式):参见横向各向同性介质声波传播特性研究,《应用声学》1984年第04期作者:蔡忠理,刘克,梁治安,p31-35
V ph = C 11 ρ , V pv = C 33 ρ ,
V sh = C 66 ρ , V pv = C 44 ρ ,
可得到饱和流体页岩的声速随油母质含量的变化曲线(图2)。

Claims (1)

1.利用声速判定页岩中的有机质含量的方法,其特征是利用含有有机质的页岩中的声速与有机质含量的变化的特性判断页岩中有机质的含量,页岩在有机质含量变化后纵波和横波声速发生明显的变化;采用超声的方法无损检测页岩的有机质的含量,即采用测量页岩的声速随有机质含量的变化来检测页岩的有机质的含量;通过发射和接收超声换能器,测量声波在页岩样品的传播时间和页岩样品的长度,得到页岩的声速,根据声速与页岩中有机质含量的关系来判定页岩中有机质的含量:
1)由不同矿物组成无油的自由“页岩”即干页岩,不同矿物被认为是各向同性的或各向异性的;不同矿物和它们各自的体积百分比,从X-射线衍射分析得到的;包括粘土,石英,方解石的无机成分矿物存在于岩石的弹性常数使用Reuss-Voigt-Hill平均计算;2)考虑油母质为背景材料,使用自洽的各向异性微分有效介质模型添加“干页岩”包裹到油母质背景材料中,形成油母质和干页岩复合材料;3)添加气孔进入油母质和干页岩复合材料,使用自洽的各向异性微分有效介质模型重新估计带气孔油母质和干页岩的弹性模量;4)添加液体进入气孔中,使用Brown-Korringa模型,获得流体饱和油母质和干页岩的弹性模量;根据上述模型计算得到饱和流体页岩的弹性模量随油母质含量的变化图;
5)根据弹性模量和声速的关系得到流体饱和油母质和干页岩中的纵波和横波的声速;得到饱和流体页岩的声速随油母质含量的变化曲线;针对Bazhenov页岩根据测量到的声速来判定页岩中有机质的含量。
CN201410270856.0A 2014-06-17 2014-06-17 利用声速判定页岩中的有机质含量的方法 Expired - Fee Related CN104101649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410270856.0A CN104101649B (zh) 2014-06-17 2014-06-17 利用声速判定页岩中的有机质含量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410270856.0A CN104101649B (zh) 2014-06-17 2014-06-17 利用声速判定页岩中的有机质含量的方法

Publications (2)

Publication Number Publication Date
CN104101649A CN104101649A (zh) 2014-10-15
CN104101649B true CN104101649B (zh) 2016-08-17

Family

ID=51669972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410270856.0A Expired - Fee Related CN104101649B (zh) 2014-06-17 2014-06-17 利用声速判定页岩中的有机质含量的方法

Country Status (1)

Country Link
CN (1) CN104101649B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104374623A (zh) * 2014-11-17 2015-02-25 中国石油大学(北京) 非规则页岩样品、制备方法及超声波测试装置
CN105606703B (zh) * 2016-01-28 2018-07-31 中国石油大学(华东) 页岩吸附气和游离气的计算方法及其测量装置
CN106226400B (zh) * 2016-06-29 2019-04-30 中国石油大学(北京) 页岩各向异性测量装置及测量方法
CN106370578B (zh) * 2016-08-30 2019-05-07 中国石油天然气股份有限公司 一种确定页岩有机质孔隙度的方法及系统
CN107575218A (zh) * 2017-09-01 2018-01-12 中国地质大学(北京) 一种快速判断富有机质成熟页岩的方法及应用
CN110836836B (zh) * 2018-08-17 2022-08-16 中国石油化工股份有限公司 干酪根弹性模量测试方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO301736B1 (no) * 1994-09-14 1997-12-01 Nyfotek As Fremgangsmåte for måling av lydhastigheten samt pröveholder
IT1313693B1 (it) * 1999-11-26 2002-09-09 Eni Spa Procedimento per determinare in continuo l'interazione tra fluidi diperforazione e formazione argillose.
IT1313691B1 (it) * 1999-11-26 2002-09-09 Eni Spa Procedimento per verificare l'efficacia dei fluidi di perforazionenello stabilizzare le pareti di pozzi petroliferi.
US8906690B2 (en) * 2012-04-13 2014-12-09 Schlumberger Technology Corporation Methods for simultaneous estimation of quantitative minerology, kerogen content and maturity in gas shale and oil-bearing shale

Also Published As

Publication number Publication date
CN104101649A (zh) 2014-10-15

Similar Documents

Publication Publication Date Title
CN104101649B (zh) 利用声速判定页岩中的有机质含量的方法
Madhubabu et al. Prediction of compressive strength and elastic modulus of carbonate rocks
Mondol et al. Experimental mechanical compaction of clay mineral aggregates—Changes in physical properties of mudstones during burial
CN103760081B (zh) 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统
Chawre Correlations between ultrasonic pulse wave velocities and rock properties of quartz-mica schist
Mondol et al. Elastic properties of clay minerals
Last Mineralogical analysis of lake sediments
CN102508296B (zh) 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置
CN105221141B (zh) 一种泥页岩脆性指数预测方法
Uyanık et al. Prediction of mechanical and physical properties of some sedimentary rocks from ultrasonic velocities
CN103984027B (zh) 基于椭球体双重孔隙模型的岩石纵波速度预测方法
Croizé et al. Experimental mechanical and chemical compaction of carbonate sand
Delle Piane et al. Frequency-dependent seismic attenuation in shales: experimental results and theoretical analysis
Knuth et al. Evolution of ultrasonic velocity and dynamic elastic moduli with shear strain in granular layers
Hampton et al. Microcrack damage observations near coalesced fractures using acoustic emission
Soroush et al. Evaluation of rock properties using ultrasonic pulse technique and correlating static to dynamic elastic constants
Knez et al. Discrepancy between measured dynamic poroelastic parameters and predicted values from Wyllie’s equation for water-saturated Istebna sandstone
Ding et al. P-and S-wave velocity and anisotropy in saturated rocks with aligned cracks
Arman et al. Strength estimation of evaporitic rocks using different testing methods
Benson et al. Role of void space geometry in permeability evolution in crustal rocks at elevated pressure
Mondol Seal quality prediction using E-Poisson's ratio rock physics templet-A case study from the Norwegian Barents Sea
Hefny et al. A laboratory approach for the calibration of seismic data in the western part of the Swiss Molasse Basin: the case history of well Humilly-2 (France) in the Geneva area
Abdlmutalib et al. The relationship between lithological and geomechanical properties of tight carbonate rocks from Upper Jubaila and Arab-D Member outcrop analog, Central Saudi Arabia
Bayuk et al. Multiscale rock-physics modeling of effective elastic properties of fractured reservoir rocks
Cholach et al. Flin Flon Belt seismic anisotropy: elastic symmetry, heterogeneity, and shear-wave splitting

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160817

Termination date: 20170617