CN106483056B - 一种基于纵波波速的页岩孔隙度测量方法及测量装置 - Google Patents
一种基于纵波波速的页岩孔隙度测量方法及测量装置 Download PDFInfo
- Publication number
- CN106483056B CN106483056B CN201610859521.1A CN201610859521A CN106483056B CN 106483056 B CN106483056 B CN 106483056B CN 201610859521 A CN201610859521 A CN 201610859521A CN 106483056 B CN106483056 B CN 106483056B
- Authority
- CN
- China
- Prior art keywords
- shale
- porosity
- effective stress
- formula
- longitudinal wave
- 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
- 238000000691 measurement method Methods 0.000 title claims description 7
- 239000011148 porous material Substances 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000002245 particle Substances 0.000 claims description 5
- 238000005056 compaction Methods 0.000 claims description 4
- 230000000644 propagated effect Effects 0.000 claims description 4
- 239000013049 sediment Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 14
- 239000011435 rock Substances 0.000 abstract description 13
- 238000005259 measurement Methods 0.000 description 8
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 6
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 239000002734 clay mineral Substances 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000008398 formation water Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000010720 hydraulic oil Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000003350 kerosene Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Dispersion Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于纵波波速的页岩孔隙度测量方法及测量装置，旨在克服常规岩石孔隙度测量方法需要测量介质且测量介质无法充满页岩所有孔隙的问题，该方法基于“地层声波因素公式”建立孔隙度与纵波波速之间的联系，通过沉积物压实理论得到页岩的孔隙度；该装置由加压系统、声波发射系统、声波采集系统和数据处理系统组成，该装置的加压系统由计算机进行控制，实现自动加压，采用本发明测量页岩的孔隙度，不需要测量介质，不损害岩样性质；增压系统设定压力后增压泵能实现自动增压，自动化程度高，降低了操作者的劳动强度。
Description
技术领域
本发明属于石油、地质、矿业领域室内测量页岩孔隙度的方法和装置技术领域，尤其涉及基于纵波波速测量页岩孔隙度的方法及装置。
背景技术
页岩孔隙度是认识和评价油气藏的重要基础参数。页岩中粘土矿物含量高，常成条状分布在骨架颗粒的表面，粘土矿物与有机质赋存关系复杂，常规的孔隙度测量方法无法直接运用到页岩孔隙度测量中。
在CN202267616U中提到一种测量岩石孔隙度的方法和装置，其采用乙醇作为测量介质。若将此方法用在测量页岩孔隙度测试中会存在以下不足：富有机质页岩中有机质孔中存在吸附气，乙醇无法完全进入到纳米级的有机质孔中，且页岩层理发育，乙醇会使得页岩中的裂纹沿着弱结构面扩展。在常规的岩石孔隙度测量方法中常采用煤油和氮气作为测量介质，而页岩中粘土矿物是强水湿，且有机质对氮气具有吸附性，会使得页岩的孔隙体积测量存在一定的误差。在CN105651805A中提出运用核磁共振测量岩石的孔隙度，运用核磁共振测量岩石的孔隙度成本相对较高，且存在测量介质无法全部进入页岩孔隙中的问题。
发明内容
本发明的主要目的在于克服常规岩石孔隙度测量方法需要测量介质且测量介质无法充满页岩孔隙的问题，提出一种在不损害岩样本身性质的情况下快速准确的测量页岩孔隙度的方法。
本发明是采用以下技术方案实现：一种基于纵波波速的页岩孔隙度测量方法，它主要包括以下步骤：
S1、将页岩加工成直径为25mm，长度为50mm的柱状，在烘箱中烘24小时备用；
S2、测量页岩在有效应力为5MPa，7MPa，10MPa，15MPa，20MPa，40MPa，45MPa，50MPa，55MPa，60MPa下的纵波波速Vp；
S3、采用式1建立Vp与有效应力σ之间的联系，计算有效应力为300MPa时纵波波速Vp300：
Vp＝V0+Dσ-Be-kσ (1)
式中，Vp-纵波波速，km/s；σ-有效应力,MPa；V0、D、B、k为拟合参数,对应量纲分别为，km/s、10-4km/s/MPa、km/s、10-2MPa4；
S4、运用“声波地层因素公式”式2
式中，Vma-声波在骨架颗粒中的传播速度，km/s，采用页岩在300MPa有效应力时传播的波速Vp300代替Vma；φ-孔隙度，无量纲；X-岩性指数，无量纲，页岩孔隙结构复杂，X一般取3，得到有效应力σ下的孔隙度；
S5、基于泥质沉积物压实理论，建立孔隙度与有效应力之间的关系式3，
φ＝φ0e-bσ (3)
式中，φ0-有效应力为零时页岩的孔隙度，无量纲，b-拟合参数，无量纲。
一种基于纵波波速的页岩孔隙度测量装置，它包括加压系统、声波发射系统、声波采集系统和数据处理系统组成。加压系统可以实现对页岩施加0～90MPa的压力，采用全自动增压泵自动进行加压；声波发射系统由声波发射器和声波发射端组成，声波发射器发射的纵波频率为170KHz；声波采集系统由数据线和示波器、声波接收端等组成，声波采集系统可以记录纵波首波到达的时间，从而得到声波在页岩中的传播时间；数据处理系统运用original软件得到纵波波速Vp与有效应力σ之间拟合公式中各参数的值，同时建立孔隙度φ与有效应力σ之间的关系，拟合公式中φ0为有效应力为零时对应的页岩孔隙度。
本发明具有以下优点：(1)一种页岩孔隙度测试方法，运用类比法建立声波在页岩中传播的波速、声波在骨架中的传播速度与孔隙度三者之间的联系，同时用页岩在300MPa时的波速代替声波在骨架中的传播速度，无需对页岩组分进行分析测试同时不需要测量介质，保证精度的同时节约了时间和成本，且不损害岩样性质；(2)一种基于纵波波速的页岩孔隙度测试装置可以快速测定不同有效应力下的纵波波速，采用高精度的示波器进行声波波形采集，采集精度高。增压系统设定压力后增压泵能实现自动增压，自动化程度高，降低了操作者的劳动强度。
附图说明
图1为该方法的具体实施步骤图。
图2为一种基于纵波波速的页岩孔隙度测量装置结构示意图；
图中，1-夹持器外壳，2-液压油，3-橡胶胶套，4-页岩样品，5-声波发射器，6-计算机，7-声波接收端，8-柱塞，9-金属管线，10-阀门，11-自动增压泵，12-压力传感器，13-示波器，14-声波发射端；图中虚线部分表示导线。
图3为有效应力σ与纵波波速Vp之间拟合图。
图4为有效应力σ与孔隙度φ拟合图。
图5为本发明测量的孔隙度和He测量的孔隙度关系图。
具体实施方式
下面通过具体实施案例进一步详细说明本发明的方法的特点及使用效果，但本发明并不因此而受到任何限制。
本发明的理论基于“声波地层因素公式”，
式中，Vma-声波在骨架颗粒中的传播速度，km/s，采用页岩在300MPa有效应力时传播的波速Vp300代替Vma；φ-孔隙度，无量纲；X-岩性指数，无量纲，页岩孔隙结构复杂一般取3，得到有效应力σ下的孔隙度；该公式模仿电阻率地层因素公式(5)
Rw-的是岩样的电阻率，Ω˙m，R0-地层水电阻率，Ω˙m，b-比例系数，无量纲，φ-孔隙度，无量纲，m-地层因素，无量纲。
岩样中电流是依靠孔隙中的地层水传播的，传播路径是弯曲的，弯曲程度由胶结指数m决定。同样，声波在页岩中传播主要依靠固体传播，(所占的百分比为1-φ，φ为孔隙度)其路径也是弯曲的，弯曲程度由X决定，X表示岩性指数，X值由岩性决定。
纵波波速与有效应力之间的关系符合公式6
Vp＝V0+Dσ-Be-kσ (6)
式中，Vp-纵波波速，km/s；σ-有效应力,MPa；V0、D、B、k为拟合参数,对应量纲分别为，km/s、10-4km/s/MPa、km/s、10-2MPa4；
根据公式6可预测高有效应力下的纵波波速。本发明中采用300MPa下的纵波波速Vp300代替纵波在骨架中的传播速度Vma，具有节省了时间的优点，且得到的纵波波速更能代表页岩中固体颗粒的纵波波速；
通过公式1、公式2得到页岩不同有效应力下的孔隙度，基于压实函数式7得到页岩的孔隙度。
φ＝φ0e-bσ (7)
式中，φ0-有效应力为零时页岩的孔隙度，无量纲，b-拟合参数，无量纲。
具体实施例如下：
S1、选取四川盆地龙马溪组页岩岩样5块，将页岩加工成直径约为25mm，长度约为50mm的柱状，编号分别为A-1、A-2、A-3、A-4、A-5，在烘箱中烘24小时备用；
S2、将页岩样品放入加压系统中的岩心夹持器中，在计算机6上设定压力分别为5MPa，7MPa，10MPa，15MPa，20MPa，40MPa，45MPa，50MPa，55MPa，60MPa，当自动增压泵11压力达到设定的为5MPa时，声波发射器5将声波信号传到声波发射端14，声波接收端7接到信号后将信号传送到示波器13，计算机6采集示波器13的信号，记录有效应力为5MPa时首波到达的时间t5，μs。而后，自动增压泵11将压力自动增加至7MPa，计算机6记录有效应力7MPa时的首波到达时间t7，μs。以此类推，直至增压泵压力达到60MPa时结束测量。
S4、依据公式8得到在有效应力为5MPa，7MPa，10MPa，15MPa，20MPa，40MPa，45MPa，50MPa，55MPa，60MPa时纵波的波速(表1)
其中，L-岩样长度/cm，t-声波在岩样中传播的时间，μs，V-纵波波速，km/s；
S5、运用original软件采用公式Vp＝V0+Dσ-Be-kσ对有效应力与孔隙度进行拟合，拟合结果见(图3)，计算有效应力为300MPa时纵波波速Vp300(表2)；
S6、运用“声波地层因素公式”采用页岩在300MPa有效应力时传播的波速Vp300代替Vma；φ表示孔隙度；X表示岩性指数，页岩中X取3，得到有效应力σ下的孔隙度；
S7、基于泥质沉积物压实理论(式7)，建立孔隙度与有效应力之间的关系，得到页岩在有效应力为零时的孔隙度(图4)。
通过本发明测量得到的页岩孔隙度与He孔隙度进行对比(图5)，可以看出本发明提供的方法测量结果与He孔隙度具有一致性。
Claims (1)
1.一种基于纵波波速的页岩孔隙度测量方法，该方法包括:
S1、将页岩加工成直径为25mm，长度为50mm的柱状，在烘箱中烘24小时备用；
S2、测量页岩在有效应力为5MPa，7MPa，10MPa，15MPa，20MPa，40MPa，45MPa，50MPa，55MPa，60MPa下的纵波波速Vp；
S3、采用式1建立Vp与有效应力为σ之间的联系，计算有效应力为300MPa时纵波波速Vp300：
Vp＝V0+Dσ-Be-kσ (1)
式中，Vp-纵波波速，km/s；σ-有效应力,MPa；V0、D、B、k为拟合参数,对应量纲分别为，km/s、10-4km/s/MPa、km/s、10-2MPa4；
S4、运用“声波地层因素公式”式2
式中，Vma-声波在骨架颗粒中的传播速度，km/s，采用页岩在300MPa有效应力时传播的波速Vp300代替Vma；φ-孔隙度，无量纲；X-岩性指数，无量纲，页岩孔隙结构复杂一般取3，得到有效应力σ下的孔隙度；
S5、基于泥质沉积物压实理论，建立孔隙度与有效应力之间的关系式3，
φ＝φ0e-bσ (3)
式中，φ0-有效应力为零时页岩的孔隙度，无量纲，b-拟合参数，无量纲。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610859521.1A CN106483056B (zh) | 2016-09-28 | 2016-09-28 | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610859521.1A CN106483056B (zh) | 2016-09-28 | 2016-09-28 | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106483056A CN106483056A (zh) | 2017-03-08 |
| CN106483056B true CN106483056B (zh) | 2019-12-10 |
Family
ID=58268117
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610859521.1A Active CN106483056B (zh) | 2016-09-28 | 2016-09-28 | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106483056B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109655528B (zh) * | 2017-10-11 | 2021-09-17 | 中国石油化工股份有限公司 | 原煤吸附状态下弹性波速度分析方法 |
| CN111443025B (zh) * | 2020-04-12 | 2022-09-23 | 北京工业大学 | 一种基于纵波波速预估高温后花岗岩围岩孔隙度的方法 |
| CN113848158B (zh) * | 2020-06-28 | 2023-09-26 | 中国石油天然气股份有限公司 | 二维大岩石模型孔隙度分布测试方法及装置 |
| CN112630120B (zh) * | 2020-11-27 | 2024-04-02 | 中国科学院深圳先进技术研究院 | 一种测量涂层孔隙率模型的建立方法及该模型的使用方法 |
| CN112525803A (zh) * | 2021-02-10 | 2021-03-19 | 西南石油大学 | 基于黏土壳状模型的砂岩孔隙度有效应力系数计算方法 |
| CN112800632B (zh) * | 2021-04-01 | 2021-11-30 | 西南石油大学 | 基于粒状模型的砂岩储层渗透率有效应力系数计算方法 |
| CN113588416B (zh) * | 2021-05-13 | 2024-02-02 | 中铁大桥勘测设计院集团有限公司 | 一种礁灰岩岩体基本质量分级方法 |
| CN113776914B (zh) * | 2021-09-24 | 2023-09-08 | 广西大学 | 应力环境中结构面上下盘试件一体化制作装置和制作方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651805A (zh) * | 2016-03-29 | 2016-06-08 | 西南石油大学 | 利用核磁共振测量岩石孔隙度的方法 |
| CN105866003A (zh) * | 2016-05-16 | 2016-08-17 | 中国地质大学（武汉） | 一种页岩气储层总孔隙度获取方法及装置 |
-
2016
- 2016-09-28 CN CN201610859521.1A patent/CN106483056B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651805A (zh) * | 2016-03-29 | 2016-06-08 | 西南石油大学 | 利用核磁共振测量岩石孔隙度的方法 |
| CN105866003A (zh) * | 2016-05-16 | 2016-08-17 | 中国地质大学（武汉） | 一种页岩气储层总孔隙度获取方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| P wave velocities, anisotropy and hysteresis in ultrahigh-pressure metamorphic rocks as a function of confining pressure;Shaocheng Ji等;《Journal of Geophysical Research: Solid Earth》;20070930;第1-24页 * |
| 基于声波测井资料确定泥质砂岩储层孔隙度的方法研究;章钰等;《内蒙古石油化工》;20081231;第82-84页 * |
| 王森等.碳酸盐岩储层孔隙度超声波评价数值模拟.《地球物理学进展》.2015, * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106483056A (zh) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106483056B (zh) | 一种基于纵波波速的页岩孔隙度测量方法及测量装置 | |
| US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
| CN103278389B (zh) | 一种岩石动、静态弹性参数同步测量的方法 | |
| CN110006738B (zh) | 一种基于应力应变曲线和划痕测试的岩石脆性评价方法 | |
| CN102220865B (zh) | 一种检测灰岩地层孔隙压力的方法 | |
| CN201242522Y (zh) | 适用于低渗特低渗岩心孔隙度的测量装置 | |
| CN103528934A (zh) | 一种测量超低渗岩石渗透率应力敏感性的互相关技术 | |
| CN104101647A (zh) | 一种模拟储层条件的岩石超声波速度测试系统及方法 | |
| WO2020228386A1 (zh) | 利用声发射技术识别岩石裂缝起裂应力的方法 | |
| CN106019403B (zh) | 自生自储油气层孔隙度测量方法 | |
| CN105317435A (zh) | 一种水平井裂缝识别方法 | |
| CN110020785A (zh) | 基于深度学习的稠油热采井固井评价系统 | |
| CN108918396A (zh) | 一种用三轴应力测试致密储层岩石渗透率的方法 | |
| CN105092815A (zh) | 一种模拟储层条件的岩石声学和电学参数联合测试装置 | |
| CN109594971A (zh) | 基于核磁共振测井增强扩散气层识别因子的流体性质识别方法 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN112432996A (zh) | 瞬变电磁时域多通道多层管柱损伤检测系统及方法 | |
| CN104484573B (zh) | 一种确定地层刚性系数的方法 | |
| CN102383784A (zh) | 存储式漏层位置综合测量仪 | |
| CN106679871A (zh) | 一种利用内摩擦角获得岩石单轴声发射条件下真实地应力值的方法 | |
| CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
| CN112147050A (zh) | 油气储层岩石动态液测渗透率的方法和装置 | |
| CN108875115B (zh) | 一种确定岩石强度的方法 | |
| Duchkov et al. | Equipment for the studies of the acoustic properties of hydrate-containing samples in laboratory conditions | |
| CN111426617B (zh) | 一种地质条件下陆相湖盆致密储层孔隙度的定量预测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |