CN105068120B - 致密砂岩裂缝的声波实验方法与识别方法 - Google Patents
致密砂岩裂缝的声波实验方法与识别方法 Download PDFInfo
- Publication number
- CN105068120B CN105068120B CN201510416937.1A CN201510416937A CN105068120B CN 105068120 B CN105068120 B CN 105068120B CN 201510416937 A CN201510416937 A CN 201510416937A CN 105068120 B CN105068120 B CN 105068120B
- Authority
- CN
- China
- Prior art keywords
- wave
- coefficient
- attenuation
- rock
- tight sand
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种致密砂岩裂缝的声波实验方法,包括步骤S10,根据纵波衰减系数、横波衰减系数与微米级裂缝宽度的指数拟合关系,建立利用所述纵波衰减系数、横波衰减系数获得确定微米级裂缝宽度的公式:y=B×eA·x,式中,y为微米级裂缝宽度,x为纵波衰减系数或横波衰减系数,A为第一修正系数,范围在24~35之间,B为第二修正系数,范围在0.23~0.38之间,B值与砂岩的孔隙度呈正相关。同时,本发明还提供了一种致密砂岩裂缝的识别方法。本发明达到了在微米级上对致密砂岩裂缝进行识别的目的。
Description
技术领域
本发明涉及石油天然气勘探的测井评价技术,特别涉及一种致密砂岩裂缝的声波实验方法与识别方法。
背景技术
致密砂岩裂缝识别是测井评价的难题。过去多采用微电阻率扫描成像测井来识别、评价地层裂缝,在碳酸盐岩地层的水基泥浆井中,得到了较好地应用。但在致密砂岩微裂缝地层中,特别是在油基泥浆井中,应用效果不理想,而声波测井受井中泥浆特性影响小,因此,可用声波幅度衰减识别地层裂缝。
声波幅度衰减与裂缝特性的关系,在理论上研究较多;但是受实验条件的限制,有关声波测井的实验室内研究,国内外学者研究的相对较少。20世纪80年代,Chen在小模型井中对声波全波列进行过实验观测,其后针对模拟的硬地层和软地层在小模型井中分别开展了多极子横波测井实验研究。Winkler等在小模型井中观测了斯通利波衰减与裂缝渗透性的关系,并与理论数值分析结果作了对比。Tang等对水平及垂直裂缝与斯通利波之间的响应关系进行了理论分析和实验研究。李长文等在人造砂岩模型井中进行过全波列观测实验和数据处理分析,重点开展了储层渗透性变化及裂缝大小对声波波列及其频谱特征参数影响的研究。国内外学者关于模型井的声波测井实验研究取得了不少成果,但是实验所研究的裂缝宽度都停留在毫米级,对微米级的裂缝宽度的实验研究很少,这对致密砂岩微裂缝识别、评价应用不利。
孔隙度是指岩石中孔隙体积或岩石中未被固体物质充填的空间体积与岩石总体积的比值。
发明内容
本发明旨在至少解决现有技术中存在的问题之一。
本发明要解决的技术问题之一在于解决在微米级上对致密砂岩裂缝进行识别的问题。
为解决上述技术问题,本发明提供一种致密砂岩裂缝的声波实验方法,包括步骤S10,根据纵波衰减系数、横波衰减系数与微米级裂缝宽度的指数拟合关系,建立利用所述纵波衰减系数、横波衰减系数获得确定微米级裂缝宽度的公式:
y=B×eA·x
式中,y为微米级裂缝宽度,x为纵波衰减系数或横波衰减系数,A为根据拟合关系获得拟合曲线的斜率,范围在24~35之间,B为根据拟合关系获得拟合曲线的截距,范围在0.23~0.38之间,B值与砂岩的孔隙度呈正相关。
进一步的,包括步骤S9,根据受测量的岩心的长度和测量的声波幅值,确定纵波衰减系数或横波衰减系数,步骤中采用公式如下:
式中:x为纵波衰减系数或横波衰减系数,l为对接岩心总长度;C0为岩心零缝宽时测量的声波幅度,C为有裂缝时的声波波幅。
进一步的,还包括步骤:
S1,将一块固定长度的受测量岩心垂直锯成两块用于模拟致密砂岩裂缝,切口刨光磨平,然后抽真空,称量1次并记录;再饱和淡水,同时称量,与真空下所测重量结合利用称重法获得岩心孔隙度,之后裹上保鲜膜备用;
S2,将数字示波器与超声方波脉冲发生/接收仪连接,所述超声方波脉冲发生/接收仪连接的发送换能器和接收换能器分别装在处理后的所述两块岩石夹持器的两端;
S3,将所述两块岩心之间涂抹黄油后完全一致对合装在夹持器上,所述发送换能器和接收换能器通过耦合剂分别连接在两块岩心上;
S4,通过岩心夹持器对两块岩心加压至8~10个大气压,获取两块岩心间裂缝宽度为零对应的声波波幅;
S5,松开岩心夹持器,向所述两块岩心之间夹入塑料薄膜;
S6,固定两块岩心位置,再次对两块岩心加压至8~10个气压,获取两岩心间裂缝宽度对应的声波波幅;
S7,再向所述两块岩心之间夹入塑料薄膜,重复步骤S6;
S8,重复步骤S6、S7对两块岩心及夹入的塑料薄膜进行测量,获得夹入塑料薄膜后的声波波幅。
进一步的,所述塑料薄膜的厚度在10~30μm之间。
进一步的,受测量的所述岩心长度为5cm,直径为2.5cm。为解决上述技术问题,本发明还提供了一种致密砂岩裂缝的识别方法,所述识别方法基于上述的确定微米级裂缝宽度的公式y=B×eA·x对致密砂岩裂缝进行识别。
本发明的有益效果在于提供了一种致密砂岩裂缝的声波实验方法,本发明明确提出了利用声波全波列的纵波和横波信息评价裂缝的实验方法和实验条件,通过实验测得纵波幅度和横波幅度并通过计算得出纵波衰减系数和横波衰减系数,从而拟合出裂缝评价识别公式,达到了在微米级上对致密砂岩裂缝进行识别的目的;本发明保证了实验环境和条件的同一性,有助于测量数据的有效性,进而保证校正后声波仪器检测信号的准确性。同时,本发明还提供了一种致密砂岩裂缝的识别方法。
附图说明
图1所示为本发明实施例岩心实验仪器组装示意图;
图2所示为本发明实施例孔隙度为3.7%不同缝宽测量的全波波形图;
图3所示为本发明实施例孔隙率为4.5%时归一化后的岩心纵波幅度与裂缝宽度关系图;
图4所示为本发明实施例孔隙率为4.5%时归一化后的岩心横波幅度与裂缝宽度关系图;
图5所示为本发明实施例声波衰减系数与裂缝宽度一般关系图;
图6所示为本发明实施例归一横波化衰减系数与裂缝宽度关系图;
图7所示为本发明实施例声波测井计算裂缝宽度处理效果图。
具体实施方式
下文将结合具体实施例详细描述本发明。应当注意的是,下述实施例中描述的技术特征或者技术特征的组合不应当被认为是孤立的,它们可以被相互组合从而达到更好的技术效果。
实施例1
本实施例提供一种致密砂岩裂缝的声波实验方法,如图1~图6所示,包括步骤:
S1,将一块固定长度的受测量岩心垂直锯成两块用于模拟致密砂岩裂缝,切口刨光磨平,然后抽真空,称量1次并记录;再饱和淡水,同时称量,与真空下所测重量结合利用称重法获得岩心孔隙度,之后裹上保鲜膜备用,受测量的所述岩心长度为5cm,直径为2.5cm;
S2,如图1所示,将数字示波器与超声方波脉冲发生/接收仪连接,所述超声方波脉冲发生/接收仪连接的发送换能器和接收换能器分别装在处理后的所述两块岩石夹持器的两端,数字示波器采用DS1102E型数字示波器,所述超声方波脉冲发生/接收仪采用CTS-8077PR型超声方波脉冲发生/接收仪;
S3,将所述两块岩心之间涂抹黄油后完全一致对合装在夹持器上,所述发送换能器和接收换能器通过耦合剂分别连接在两块岩心上;
S4,通过岩心夹持器对两块岩心加压至8~10个大气压,获取两块岩心间裂缝宽度为零对应的声波波幅;
S5,松开岩心夹持器,向所述两块岩心之间夹入塑料薄膜,所述塑料薄膜的厚度在10~30μm之间;
S6,固定两块岩心位置,再次对两块岩心加压至8~10个气压,获取两岩心间裂缝宽度对应的声波波幅;
S7,再向所述两块岩心之间夹入塑料薄膜,重复步骤S6;
S8,重复步骤S6、S7对两块岩心及夹入的塑料薄膜进行测量,获得夹入塑料薄膜后的声波波幅,图2是孔隙度为3.7%不同缝宽测量的全波波形图,从图中可以看出,随着裂缝宽度的增加,波形幅度有减小趋势,特别是横波幅度的衰减非常明显,再如图3和图4所示,纵横波幅度与裂缝宽度关系明显,随着裂缝宽度的增加,声波幅度减小;
S9,根据受测量的岩心的长度和测量的声波幅值,确定纵波衰减系数或横波衰减系数,步骤中采用公式如下:
式中:x为纵波衰减系数或横波衰减系数,l为对接岩心总长度;C0为岩心零缝宽时测量的声波幅度,C为有裂缝时的声波波幅;
S10,如图5所示,根据纵波衰减系数、横波衰减系数与微米级裂缝宽度的指数拟合关系,建立利用所述纵波衰减系数、横波衰减系数获得确定微米级裂缝宽度的公式:
y=B×eA·x
式中,y为微米级裂缝宽度,x为纵波衰减系数或横波衰减系数,A为根据拟合关系获得拟合曲线的斜率,范围在24~35之间,B为根据拟合关系获得拟合曲线的截距,范围在0.23~0.38之间,B值与砂岩的孔隙度呈正相关。
在致密砂岩地层中,利用井中声波全波列测量仪器检测到纵横波幅度信息,然后根据步骤S9中的公式转化为纵横波衰减系数后,依据地层孔隙度范围选取合适的A与B后带入步骤S10中的公式计算得到裂缝宽度。
考虑到岩心孔隙度对声能量的吸收,本次实验考察了不同孔隙度情况下纵波、横波衰减系数与裂缝宽度的关系。图6所示为孔隙度4.5%、6.3%、7.5%情况下归一化横波衰减系数与裂缝宽度的关系,从图中可见两者具有很好的指数关系,且孔隙度越大,相同衰减系数对应的裂缝宽度越大。所以,利用声波衰减系数可以建立计算裂缝宽度的方法。
实施例2
本实施例提供了一种致密砂岩裂缝的识别方法,所述识别方法基于上述实施例1的确定微米级裂缝宽度的公式y=B×eA·x对致密砂岩裂缝进行识别,如图7所示,其中,第1道显示深度,第2道为岩性测井曲线,第4道为地层孔隙度计算曲线,充填部分为孔隙空间,第5到为归一化的纵波、横波幅度曲线,第6道为纵横波衰减曲线,第7道为裂缝宽度曲线,包括由纵波衰减系数计算的裂缝宽度和由横波衰减系数计算的裂缝宽度。从图中最右侧的一道可见,由纵波、横波衰减系数计算的裂缝宽度显示该井在6907-6916米、6929-6945米两段地层裂缝宽度较大,属于裂缝发育地层;同时,对比两段裂缝宽度数值大小,下段地层6929-6945米的裂缝宽度又比下段地层6907-6916米的数值大,说明下段裂缝更为发育。
所以,利用根据纵横波幅度估算裂缝宽度的方法及公式,可以对砂岩地层的裂缝进行计算,并评价地层裂缝发育状况。
本发明提供了一种致密砂岩裂缝的声波实验方法,本发明明确提出了利用声波全波列的纵波和横波信息评价裂缝的实验方法和实验条件,通过实验测得纵波幅度和横波幅度并通过计算得出纵波衰减系数和横波衰减系数,从而拟合出裂缝评价识别公式,达到了在微米级上对致密砂岩裂缝进行识别的目的;本发明保证了实验环境和条件的同一性,有助于测量数据的有效性,进而保证校正后声波仪器检测信号的准确性;同时,本发明还提供了一种致密砂岩裂缝的识别方法。
本文虽然已经给出了本发明的一些实施例,但是本领域的技术人员应当理解,在不脱离本发明精神的情况下,可以对本文的实施例进行改变。上述实施例只是示例性的,不应以本文的实施例作为本发明权利范围的限定。
Claims (4)
1.一种致密砂岩裂缝的声波实验方法,其特征在于,包括步骤为:
S1,将一块固定长度的受测量岩心垂直锯成两块用于模拟致密砂岩裂缝,切口刨光磨平,然后抽真空,称量1次并记录;再饱和淡水,同时称量,与真空下所测重量结合利用称重法获得岩心孔隙度,之后裹上保鲜膜备用;
S2,将数字示波器与超声方波脉冲发生/接收仪连接,所述超声方波脉冲发生/接收仪连接的发送换能器和接收换能器分别装在处理后的两块岩心夹持器的两端;
S3,将所述两块岩心之间涂抹黄油后完全一致对合装在夹持器上,所述发送换能器和接收换能器通过耦合剂分别连接在两块岩心上;
S4,通过岩心夹持器对两块岩心加压至8~10个大气压,获取两块岩心间裂缝宽度为零对应的声波波幅;
S5,松开岩心夹持器,向所述两块岩心之间夹入塑料薄膜;
S6,固定两块岩心位置,再次对两块岩心加压至8~10个气压,获取两岩心间裂缝宽度对应的声波波幅;
S7,再向所述两块岩心之间夹入塑料薄膜,重复步骤S6;
S8,重复步骤S6、S7对两块岩心及夹入的塑料薄膜进行测量,获得夹入塑料薄膜后的声波波幅;
S9,根据受测量的岩心的长度和测量的声波幅值,确定纵波衰减系数或横波衰减系数,步骤中采用公式如下:
式中:x为纵波衰减系数或横波衰减系数,l为对接岩心总长度;C0为岩心零缝宽时测量的声波幅度,C为有裂缝时的声波波幅;
S10,根据纵波衰减系数、横波衰减系数与微米级裂缝宽度的指数拟合关系,建立利用所述纵波衰减系数、横波衰减系数获得确定微米级裂缝宽度的公式:
y=B×eA·x
式中,y为微米级裂缝宽度,x为纵波衰减系数或横波衰减系数,A为根据拟合关系获得拟合曲线的斜率,范围在24~35之间,B为根据拟合关系获得拟合曲线的截距,范围在0.23~0.38之间,B值与砂岩的孔隙度呈正相关。
2.根据权利要求1所述的致密砂岩裂缝的声波实验方法,其特征在于,所述塑料薄膜的厚度在10~30μm之间。
3.根据权利要求1所述的致密砂岩裂缝的声波实验方法,其特征在于,受测量的所述岩心长度为5cm,直径为2.5cm。
4.一种致密砂岩裂缝的识别方法,其特征在于,基于如权利要求1至3任一项所述的确定微米级裂缝宽度的公式y=B×e A·x对致密砂岩裂缝进行识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510416937.1A CN105068120B (zh) | 2015-07-16 | 2015-07-16 | 致密砂岩裂缝的声波实验方法与识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510416937.1A CN105068120B (zh) | 2015-07-16 | 2015-07-16 | 致密砂岩裂缝的声波实验方法与识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105068120A CN105068120A (zh) | 2015-11-18 |
CN105068120B true CN105068120B (zh) | 2019-01-22 |
Family
ID=54497535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510416937.1A Active CN105068120B (zh) | 2015-07-16 | 2015-07-16 | 致密砂岩裂缝的声波实验方法与识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105068120B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106568846B (zh) * | 2016-10-25 | 2019-10-29 | 长江大学 | 岩心裂缝宽度下快慢横波衰减系数测量方法 |
CN111122572A (zh) * | 2019-12-30 | 2020-05-08 | 中国地质大学(武汉) | 一种分析石英粒内缝发育程度的方法 |
CN112065361B (zh) * | 2020-09-16 | 2021-03-12 | 成都理工大学 | 一种基于声波衰减的致密储层含气饱和度确定方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183585A (zh) * | 2011-03-09 | 2011-09-14 | 西南石油大学 | 一种岩心选样的方法 |
CN102900425A (zh) * | 2012-10-29 | 2013-01-30 | 中国石油大学(华东) | 一种新的页岩气井压裂裂缝监测方法 |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
CN104345346A (zh) * | 2014-10-30 | 2015-02-11 | 中国石油天然气股份有限公司 | 一种获取裂缝宽度的方法 |
-
2015
- 2015-07-16 CN CN201510416937.1A patent/CN105068120B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183585A (zh) * | 2011-03-09 | 2011-09-14 | 西南石油大学 | 一种岩心选样的方法 |
CN102900425A (zh) * | 2012-10-29 | 2013-01-30 | 中国石油大学(华东) | 一种新的页岩气井压裂裂缝监测方法 |
CN103256046A (zh) * | 2013-04-28 | 2013-08-21 | 北京大学 | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 |
CN104345346A (zh) * | 2014-10-30 | 2015-02-11 | 中国石油天然气股份有限公司 | 一种获取裂缝宽度的方法 |
Non-Patent Citations (4)
Title |
---|
carbon dioxide gas permeability of coal core samples and estimation of fracture aperture width;Phung Quoc Huy 等;《international journal of coal geology》;20101231(第83期);第1-10页 |
孔洞型碳酸盐岩超声波衰减特性的研究及应用;周龙涛;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20150315(第03期);第1-84页 |
孔洞结构对超声波衰减特性的影响研究;梁利喜 等;《岩石力学与工程学报》;20150531;第34卷(第增1期);摘要,第3210页第3节,第3211-3212页第4.2节 |
裂缝模型声波衰减系数的数值模拟;陈乔 等;《地球物理学报》;20120630;第55卷(第6期);第2044-2052页 |
Also Published As
Publication number | Publication date |
---|---|
CN105068120A (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105068120B (zh) | 致密砂岩裂缝的声波实验方法与识别方法 | |
CN103244104B (zh) | 利用频散校正提取套管波衰减率的方法 | |
CN104501939B (zh) | 一种利用单水听器测量非消声水池混响时间的反演方法 | |
CN110456418A (zh) | 阵列声波成像测井资料的处理和解释方法 | |
CN107907595A (zh) | 声波斜入射材料吸声系数与隔声量测量装置及其测量方法 | |
CN202832520U (zh) | 在声波测井中消除直达波干扰的***及声波测井仪 | |
CN103513270A (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
CN110348135A (zh) | 一种随钻声波测井评价地层渗透率的方法 | |
CN102169188A (zh) | 一种基于Morlet谱勘测油气的方法 | |
CN107153224A (zh) | 检波器动态性能综合测试与评价方法 | |
CN107436451A (zh) | 一种自动计算地震数据光缆耦合噪声强弱程度的振幅谱方法 | |
CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
CN106568846A (zh) | 岩心裂缝宽度下快慢横波衰减系数测量方法 | |
CN109324113B (zh) | 一种混响水池壁面平均吸声系数的测量方法 | |
Tran et al. | A comparison of shear wave velocity profiles from SASW, MASW, and ReMi techniques | |
CN105723214B (zh) | 通过回波分析对部件进行无损超声检测的方法 | |
CN105888648B (zh) | 基于音频声波电声***的油井动液面深度探测装置及方法 | |
CN201554458U (zh) | 一种数字声波和变密度综合测井仪 | |
Donald et al. | Adaptation of a triaxial cell for ultrasonic P-wave attenuation, velocity and acoustic emission measurements | |
WO2019036351A1 (en) | NETWORK OF ULTRASONIC TRANSDUCERS WITH LIQUID COUPLING FOR MEASURING ROCK ELASTIC PROPERTIES | |
CN105954170B (zh) | 一种考虑背景噪声信号量的核磁孔隙度计算方法 | |
CN108872402A (zh) | 超声波巴特沃斯、汉宁窗组合帯阻滤波方法 | |
CN110307946B (zh) | 一种燃气管道泄漏检测用模拟测量装置及其测试方法 | |
Huanran et al. | Dynamic fluid transport property of hydraulic fractures and its evaluation using acoustic logging | |
CN105651863A (zh) | 一种虫蛀孔洞玉米种子检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |