[go: up one dir, main page]

CN114813509B - Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference - Google Patents

Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference Download PDF

Info

Publication number
CN114813509B
CN114813509B CN202210438404.3A CN202210438404A CN114813509B CN 114813509 B CN114813509 B CN 114813509B CN 202210438404 A CN202210438404 A CN 202210438404A CN 114813509 B CN114813509 B CN 114813509B
Authority
CN
China
Prior art keywords
time difference
porosity
rock sample
wave time
formula
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
Application number
CN202210438404.3A
Other languages
Chinese (zh)
Other versions
CN114813509A (en
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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202210438404.3A priority Critical patent/CN114813509B/en
Publication of CN114813509A publication Critical patent/CN114813509A/en
Application granted granted Critical
Publication of CN114813509B publication Critical patent/CN114813509B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/011Velocity or travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0232Glass, 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 Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种利用声波时差计算岩石孔隙度的压实校正系数确定方法,其特征在于:基于怀利公式声波计算孔隙度的方法,以干燥后岩样纵波时差利用怀利公式计算的孔隙度为未压实的声波计算孔隙度,以抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度为压实过后的声波计算孔隙度,根据岩样在此两种状态下利用怀利公式计算的孔隙度与压实校正系数的物理关系,建立了一种利用声波时差计算岩石孔隙度的压实校正系数确定方法,对于石油勘探开发油气储量计算具有重要的意义。

The invention discloses a method for determining a compaction correction coefficient for calculating rock porosity using acoustic wave time difference. The method is characterized in that: based on the method for calculating porosity using acoustic wave time difference of the Wylie formula, the porosity calculated using the Wylie formula by using the longitudinal wave time difference of a rock sample after drying is the uncompacted acoustic wave calculated porosity, and the porosity calculated using the Wylie formula by using the longitudinal wave time difference of a rock sample saturated with pores filled with vacuum distilled water is the compacted acoustic wave calculated porosity. According to the physical relationship between the porosity calculated using the Wylie formula and the compaction correction coefficient of the rock sample in these two states, a method for determining the compaction correction coefficient for calculating rock porosity using acoustic wave time difference is established, which is of great significance for the calculation of oil and gas reserves in oil exploration and development.

Description

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法A method for determining compaction correction coefficients for calculating rock porosity using acoustic wave transit time

技术领域Technical Field

本发明属于石油天然气勘探开发领域,具体涉及一种利用声波时差计算岩石孔隙度的压实校正系数确定方法。The invention belongs to the field of oil and gas exploration and development, and in particular relates to a method for determining a compaction correction coefficient for calculating rock porosity using acoustic wave time difference.

背景技术Background technique

岩石孔隙度是石油天然气勘探开发领域关乎地下油气储量计算的一项重要参数,从怀利(Wyllie)发明时间平均公式以来,声波计算孔隙度被学术界广泛的研究使用,声波计算孔隙度的怀利公式如下:Rock porosity is an important parameter related to the calculation of underground oil and gas reserves in the field of oil and gas exploration and development. Since Wyllie invented the time-averaged formula, the porosity calculated by acoustic waves has been widely studied and used in academia. The Wyllie formula for calculating porosity by acoustic waves is as follows:

式中,φs为利用怀利公式声波计算的孔隙度,无量纲;Δt为实测岩样的纵波时差,μs/ft;Δtma为岩样骨架的纵波时差,μs/ft;Δtf为流体的纵波时差,μs/ft。Wherein, φs is the porosity calculated by acoustic wave using the Wylie formula, dimensionless; Δt is the longitudinal wave time difference of the measured rock sample, μs/ft; Δtma is the longitudinal wave time difference of the rock sample skeleton, μs/ft; Δtf is the longitudinal wave time difference of the fluid, μs/ft.

在实际使用时,研究者们发现,使用怀利公式声波计算的岩石孔隙度往往比实测值大得多,所以研究者引入了声波时差计算岩石孔隙度的压实校正系数Cp,后经校正关系如下:In actual use, researchers found that the rock porosity calculated by using the Wylie formula is often much larger than the measured value, so the researchers introduced the compaction correction coefficient Cp for calculating the rock porosity using the acoustic time difference. After correction, the relationship is as follows:

式中,Cp为声波时差计算岩石孔隙度的压实校正系数,无量纲。Where Cp is the compaction correction coefficient for calculating rock porosity using acoustic time difference, which is dimensionless.

对怀利公式声波计算孔隙度的公式进行压实校正后,被学术界普遍认可,在科学研究和项目工程开发使用一直沿用至今,其中最重要的就是声波时差计算岩石孔隙度的压实校正系数Cp的确定。After the compaction correction was performed on the Wylie formula for calculating porosity by acoustic waves, it was generally recognized by the academic community and has been used in scientific research and project engineering development to this day. The most important of these is the determination of the compaction correction coefficient Cp for calculating rock porosity by acoustic wave time difference.

发明内容Summary of the invention

本发明的目的在于:解决声波计算孔隙度压实校正的问题,提供一种利用声波时差计算岩石孔隙度的压实校正系数确定方法。The purpose of the present invention is to solve the problem of compaction correction of porosity calculated by acoustic waves and to provide a method for determining the compaction correction coefficient for calculating rock porosity using acoustic wave time difference.

本发明采用的技术方案如下:The technical solution adopted by the present invention is as follows:

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法,其确定步骤如下:A method for determining a compaction correction coefficient for calculating rock porosity using acoustic wave time difference, wherein the determination steps are as follows:

步骤1.1、将岩样在105℃温度环境干燥72小时,去除孔隙内的流体和孔隙表面吸附的挥发性物质,然后对干燥岩样进行声波测量,得到干燥后岩样的纵波时差Δtcd,利用怀利公式计算干燥后岩样的声波孔隙度,得到干燥后岩样的纵波时差利用怀利公式计算的孔隙度φsdStep 1.1, drying the rock sample at 105°C for 72 hours to remove the fluid in the pores and the volatile substances adsorbed on the pore surface, then performing acoustic wave measurement on the dried rock sample to obtain the longitudinal wave time difference Δtcd of the dried rock sample, and calculating the acoustic wave porosity of the dried rock sample using the Wylie formula to obtain the porosity φsd calculated by the longitudinal wave time difference of the dried rock sample using the Wylie formula;

步骤1.2、将经过步骤1.1的岩样抽真空进行蒸馏水饱和,将岩样孔隙用蒸馏水填实,直至吸水后的岩样重量不再增加,然后对岩样进行声波测量,得到抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs,利用怀利公式计算孔隙被蒸馏水填实饱和岩样的声波孔隙度,得到抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度φssStep 1.2, vacuum the rock sample obtained in step 1.1 and saturate it with distilled water, fill the pores of the rock sample with distilled water until the weight of the rock sample after absorbing water no longer increases, then perform acoustic wave measurement on the rock sample to obtain the longitudinal wave time difference Δ tcs of the rock sample filled with pores by vacuum saturated distilled water, calculate the acoustic wave porosity of the rock sample filled with pores by distilled water using the Wylie formula, and obtain the porosity φ ss calculated by the longitudinal wave time difference of the rock sample filled with pores by vacuum saturated distilled water using the Wylie formula;

步骤1.3、利用干燥后岩样的纵波时差Δtcd和抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs分别利用怀里公式计算的孔隙度,φsd和φss的物理关系,得到声波时差计算岩石孔隙度的压实校正系数CpStep 1.3: Use the longitudinal wave time difference Δtcd of the dried rock sample and the longitudinal wave time difference Δtcs of the vacuum saturated distilled water filled rock sample to calculate the porosity using the Huyley formula, and the physical relationship between φsd and φss to obtain the compaction correction coefficient Cp for calculating rock porosity using the acoustic wave time difference:

进一步,所述步骤1.3.1的具体步骤如下:Further, the specific steps of step 1.3.1 are as follows:

步骤1.3.1以干燥后岩样的纵波时差利用怀利公式计算的孔隙度φsd为未压实的声波计算孔隙度,以抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度φss为压实过后的声波计算孔隙度,根据声波压实校正系数的物理定义,得到声波计算孔隙度的压实校正系数CpStep 1.3.1 The porosity φ sd calculated by the longitudinal wave time difference of the dried rock sample using the Wylie formula is the uncompacted acoustic wave calculated porosity, and the porosity φ ss calculated by the longitudinal wave time difference of the rock sample saturated with vacuum distilled water and filled with pores using the Wylie formula is the compacted acoustic wave calculated porosity. According to the physical definition of the acoustic wave compaction correction coefficient, the compaction correction coefficient C p of the acoustic wave calculated porosity is obtained:

式中,Cp为声波时差计算岩石孔隙度的压实校正系数,无量纲;φss为抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度,无量纲;φsd为干燥后岩样的纵波时差利用怀利公式计算的孔隙度,无量纲;Δtcs为抽真空饱和蒸馏水填实孔隙岩样的纵波时差,μs/ft;Δtcd为干燥后岩样的纵波时差,μs/ft;Δtma为岩样骨架的纵波时差,μs/ft;Δtf为流体的纵波时差,μs/ft。Wherein, Cp is the compaction correction coefficient for calculating rock porosity by acoustic time difference, dimensionless; φss is the porosity calculated by using the Wylie formula using the longitudinal wave time difference of the rock sample saturated with vacuum distilled water to fill the pores, dimensionless; φsd is the porosity calculated by using the Wylie formula using the longitudinal wave time difference of the rock sample after drying, dimensionless; Δtcs is the longitudinal wave time difference of the rock sample saturated with vacuum distilled water to fill the pores, μs/ft; Δtcd is the longitudinal wave time difference of the rock sample after drying, μs/ft; Δtma is the longitudinal wave time difference of the rock sample skeleton, μs/ft; Δtf is the longitudinal wave time difference of the fluid, μs/ft.

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法,包括如下步骤:A method for determining a compaction correction coefficient for calculating rock porosity using acoustic wave time difference comprises the following steps:

步骤(1)、将岩样在105℃温度环境干燥72小时,去除孔隙内的流体和孔隙表面吸附的挥发性物质,然后对干燥岩样进行声波测量,得到干燥后岩样的纵波时差Δtcd,利用怀利公式计算干燥后岩样的声波孔隙度,得到干燥后岩样的纵波时差利用怀利公式计算的孔隙度φsdStep (1), drying the rock sample at 105° C. for 72 hours to remove the fluid in the pores and the volatile substances adsorbed on the pore surface, then performing acoustic wave measurement on the dried rock sample to obtain the longitudinal wave time difference Δ tcd of the dried rock sample, and calculating the acoustic wave porosity of the dried rock sample using the Wylie formula to obtain the porosity φ sd of the dried rock sample calculated using the longitudinal wave time difference using the Wylie formula;

步骤(2)、将经过步骤(1)的岩样抽真空进行蒸馏水饱和,将岩样孔隙用蒸馏水填实,直至吸水后的岩样重量不再增加,然后对岩样进行声波测量,得到抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs,利用怀利公式计算孔隙被蒸馏水填实饱和岩样的声波孔隙度,得到抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度φssStep (2), vacuuming the rock sample after step (1) to saturate it with distilled water, filling the pores of the rock sample with distilled water until the weight of the rock sample after absorbing water no longer increases, then performing acoustic wave measurement on the rock sample to obtain the longitudinal wave time difference Δtcs of the rock sample filled with pores by vacuuming and saturating it with distilled water, using the Wylie formula to calculate the acoustic wave porosity of the rock sample filled with pores by distilled water, and obtaining the porosity φss calculated by the longitudinal wave time difference of the rock sample filled with pores by vacuuming and saturating it with distilled water using the Wylie formula;

步骤(3)、利用干燥后岩样的纵波时差Δtcd和抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs分别利用怀里公式计算的孔隙度,φsd和φss的物理关系,得到声波时差计算岩石孔隙度的压实校正系数CpStep (3), using the longitudinal wave time difference Δtcd of the dried rock sample and the longitudinal wave time difference Δtcs of the vacuum saturated distilled water filled rock sample to calculate the porosity, the physical relationship between φsd and φss using the Huyley formula, and obtain the compaction correction coefficient Cp for calculating the rock porosity by the acoustic wave time difference.

综上所述,由于采用了上述技术方案,本发明的有益效果是:In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are:

本发明中,将干燥后的岩样视作无有效介质填充孔隙,将抽真空饱和蒸馏水后的岩样视作孔隙被填实,采用两种状态下岩样的纵波时差怀利公式计算的孔隙度与压实校正系数的物理关系,确定了声波时差计算岩石孔隙度的压实校正系数CpIn the present invention, the dried rock sample is regarded as having no effective medium filling the pores, and the rock sample saturated with distilled water after vacuuming is regarded as having filled the pores. The physical relationship between the porosity calculated by the longitudinal wave time difference Wylie formula of the rock sample in the two states and the compaction correction coefficient is adopted to determine the compaction correction coefficient Cp of the rock porosity calculated by the acoustic wave time difference.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为岩样的基础参数Figure 1 shows the basic parameters of the rock sample

图2为不同岩石骨架及流体纵波时差参数Figure 2 shows the longitudinal wave time difference parameters of different rock skeletons and fluids

图3为本发明确定Cp计算的孔隙度φs与实测孔隙度φt的1:1对比图FIG. 3 is a 1:1 comparison diagram of the porosity φs calculated by the present invention to determine Cp and the measured porosity φt

具体实施方式Detailed ways

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明,应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法:A method for determining the compaction correction coefficient for calculating rock porosity using acoustic wave time difference:

式中,φsd为干燥后岩样的纵波时差利用怀利公式计算的孔隙度,无量纲;Δtcs为抽真空饱和蒸馏水填实孔隙岩样的纵波时差,μs/ft;φss为抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度,无量纲;Δtcd为干燥后岩样的纵波时差,μs/ft。Wherein, φ sd is the porosity calculated by the Wylie formula using the longitudinal wave time difference of the rock sample after drying, dimensionless; Δ tcs is the longitudinal wave time difference of the rock sample filled with pores by vacuum saturated distilled water, μs/ft; φ ss is the porosity calculated by the Wylie formula using the longitudinal wave time difference of the rock sample filled with pores by vacuum saturated distilled water, dimensionless; Δ tcd is the longitudinal wave time difference of the rock sample after drying, μs/ft.

本发明的确定过程如下:The determination process of the present invention is as follows:

以干燥后的岩样视作无有效介质填充孔隙,将抽真空饱和蒸馏水后的岩样视作孔隙被填实,采用两种状态下岩样的纵波时差怀利公式计算的孔隙度与压实校正系数的物理关系,确定了声波时差计算岩石孔隙度的压实校正系数CpThe dried rock samples were regarded as having no effective medium filling the pores, and the rock samples saturated with distilled water after vacuuming were regarded as having filled the pores. The physical relationship between the porosity and the compaction correction coefficient calculated by the longitudinal wave transit time Wylie formula of the rock samples in the two states was adopted to determine the compaction correction coefficient C p of the rock porosity calculated by acoustic wave transit time.

干燥后岩样的纵波时差利用怀利公式计算的孔隙度,计算公式如下:The porosity of the dried rock sample is calculated using the Wylie formula using the longitudinal wave time difference. The calculation formula is as follows:

抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度,计算公式如下:The porosity is calculated by the longitudinal wave time difference of the rock sample filled with vacuum distilled water saturated with pores using the Wylie formula. The calculation formula is as follows:

根据声波计算孔隙度的压实校正系数的物理意义,则可得声波时差计算岩石孔隙度的压实校正系数,计算公式:According to the physical meaning of the compaction correction coefficient of porosity calculated by acoustic wave, the compaction correction coefficient of rock porosity calculated by acoustic wave time difference can be obtained. The calculation formula is:

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法,包括步骤如下:A method for determining a compaction correction coefficient for calculating rock porosity using acoustic wave time difference comprises the following steps:

步骤1:将岩样在105℃温度环境干燥72小时,去除孔隙内的流体和孔隙表面吸附的挥发性物质,然后对干燥岩样进行声波测量,得到干燥后岩样的纵波时差Δtcd,利用怀利公式计算干燥后岩样的声波孔隙度,得到干燥后岩样的纵波时差利用怀利公式计算的孔隙度φsdStep 1: Dry the rock sample at 105°C for 72 hours to remove the fluid in the pores and the volatile substances adsorbed on the pore surface, then perform acoustic wave measurement on the dried rock sample to obtain the longitudinal wave time difference Δtcd of the dried rock sample, calculate the acoustic wave porosity of the dried rock sample using the Wylie formula, and obtain the porosity φsd calculated by the longitudinal wave time difference of the dried rock sample using the Wylie formula;

步骤2:将经过步骤1的岩样抽真空进行蒸馏水饱和,将岩样孔隙用蒸馏水填实,直至吸水后的岩样重量不再增加,然后对岩样进行声波测量,得到抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs,利用怀利公式计算孔隙被蒸馏水填实饱和岩样的声波孔隙度,得到抽真空蒸馏水饱和填实孔隙岩样的纵波时差利用怀利公式计算的孔隙度φssStep 2: vacuum the rock sample after step 1 and saturate it with distilled water, fill the pores of the rock sample with distilled water until the weight of the rock sample after absorbing water no longer increases, then perform acoustic wave measurement on the rock sample to obtain the longitudinal wave time difference Δ tcs of the rock sample filled with pores by vacuum saturated distilled water, calculate the acoustic wave porosity of the rock sample filled with pores by distilled water using the Wylie formula, and obtain the porosity φ ss calculated by the longitudinal wave time difference of the rock sample filled with pores by vacuum saturated distilled water using the Wylie formula;

步骤3:利用干燥后岩样的纵波时差Δtcd和抽真空饱和蒸馏水填实孔隙岩样的纵波时差Δtcs分别利用怀里公式计算的孔隙度,φsd和φss的物理关系,得到声波时差计算岩石孔隙度的压实校正系数CpStep 3: Use the longitudinal wave time difference Δtcd of the dried rock sample and the longitudinal wave time difference Δtcs of the vacuum saturated distilled water filled rock sample to calculate the porosity using the Huyley formula, and the physical relationship between φsd and φss to obtain the compaction correction coefficient Cp for calculating rock porosity using the acoustic wave time difference.

实施例Example

一种利用声波时差计算岩石孔隙度的压实校正系数确定方法:A method for determining the compaction correction coefficient for calculating rock porosity using acoustic wave time difference:

该实施例采用不同比例砂岩骨架的成岩岩样采用本发明提出的方法进行计算校验,岩样的基础参数如图1所示,不同岩石骨架及流体纵波时差参数如图2所示。本发明确定Cp计算的孔隙度φs与实测孔隙度φt的1:1对比图如图3所示,较大的R2(0.873)与较小的AAREP(1.88%),说明采用本发明提出的声波时差计算岩石孔隙度的压实校正系数Cp的计算方法能够很好应用到声波计算孔隙度当中,预测精度高。This embodiment uses diagenetic rock samples with different proportions of sandstone skeletons to perform calculation verification using the method proposed in the present invention. The basic parameters of the rock samples are shown in Figure 1, and the longitudinal wave time difference parameters of different rock skeletons and fluids are shown in Figure 2. The 1:1 comparison chart of the porosity φs calculated by the present invention and the measured porosity φt is shown in Figure 3. The larger R2 (0.873) and the smaller AAREP (1.88%) indicate that the method for calculating the compaction correction coefficient Cp of rock porosity calculated by the acoustic wave time difference proposed in the present invention can be well applied to the acoustic wave porosity calculation, and the prediction accuracy is high.

Claims (1)

1. A compaction correction factor determination method for calculating rock porosity by using acoustic time difference comprises the following steps:
Step (1), drying a rock sample in a temperature environment of 105 ℃ for 72 hours, removing fluid in pores and volatile substances adsorbed on the surfaces of the pores, then carrying out acoustic wave measurement on the dried rock sample to obtain a longitudinal wave time difference delta tcd of the dried rock sample, calculating the acoustic wave porosity of the dried rock sample by using a Huai equation, and obtaining a porosity phi sd calculated by using the Huai equation of the longitudinal wave time difference of the dried rock sample, wherein the calculation formula is as follows:
Wherein phi sd is the porosity calculated by using a Chlamydia formula and is dimensionless; Δ tcd is the longitudinal wave time difference of the dried rock sample, μs/ft; Δ tma is the longitudinal wave time difference of the rock sample skeleton, μs/ft; Δ tf is the longitudinal wave time difference of the fluid, μs/ft;
Step (2), vacuumizing the rock sample subjected to the step (1) to saturate distilled water, filling the pores of the rock sample with distilled water until the weight of the rock sample after water absorption is not increased, then performing acoustic wave measurement on the rock sample to obtain a longitudinal wave time difference delta tcs of the rock sample filled with the vacuumized saturated distilled water, calculating the acoustic wave porosity of the rock sample filled with distilled water by using a Chler formula, and obtaining a porosity phi ss calculated by using the Chler formula of the longitudinal wave time difference of the rock sample filled with the vacuumized distilled water, wherein the calculation formula is as follows:
Wherein phi ss is the porosity calculated by using a Chlamydia formula and is dimensionless, wherein the longitudinal wave time difference of the vacuumized distilled water saturated filled pore rock sample is calculated by using the Chlamydia formula; delta tcs is the longitudinal wave time difference of the vacuumized saturated distilled water filled pore rock sample;
And (3) using the longitudinal wave time difference delta tcd of the dried rock sample and the longitudinal wave time difference delta tcs of the vacuumized saturated distilled water to fill up the porosity calculated by using the Win formula respectively, and obtaining a compaction correction coefficient C p of the rock porosity calculated by using the acoustic wave time difference according to the physical relationship between phi sd and phi ss, wherein the calculation formula is as follows:
wherein C p is a compaction correction coefficient for calculating the porosity of the rock by acoustic time difference, and the compaction correction coefficient is dimensionless.
CN202210438404.3A 2022-04-21 2022-04-21 Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference Active CN114813509B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210438404.3A CN114813509B (en) 2022-04-21 2022-04-21 Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210438404.3A CN114813509B (en) 2022-04-21 2022-04-21 Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference

Publications (2)

Publication Number Publication Date
CN114813509A CN114813509A (en) 2022-07-29
CN114813509B true CN114813509B (en) 2024-06-11

Family

ID=82506630

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210438404.3A Active CN114813509B (en) 2022-04-21 2022-04-21 Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference

Country Status (1)

Country Link
CN (1) CN114813509B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001255186A (en) * 2000-03-09 2001-09-21 Matsushita Electric Ind Co Ltd Flow measurement device
CN102129084A (en) * 2010-12-17 2011-07-20 中国石油天然气股份有限公司 Method and device for acquiring seismic thin reservoir speed through well control
JP2011185706A (en) * 2010-03-08 2011-09-22 Central Res Inst Of Electric Power Ind Method and device for monitoring surface defect using laser ultrasonic waves
CN103513270A (en) * 2012-06-21 2014-01-15 中国石油天然气集团公司 Gas reservoir identification evaluating method based on acoustic characteristic of rock and device thereof
CN104316978A (en) * 2014-10-29 2015-01-28 中国石油天然气股份有限公司 Near-surface three-dimensional velocity field research method and device for geophysical
FR3026313A1 (en) * 2014-09-29 2016-04-01 Francois Parmentier METHOD OF CHROMATOGRAPHY ON A MULTICAPILLARY TRIM
CN105804737A (en) * 2016-05-17 2016-07-27 西南石油大学 Method for solving formation porosity on basis of iterative algorithm
WO2019049442A1 (en) * 2017-09-08 2019-03-14 ソニー株式会社 Microparticle measurement device, information processing device, and information processing method
CN109989746A (en) * 2017-12-29 2019-07-09 中国石油天然气股份有限公司 Method and device for evaluating carbonate reservoir
JP2020007730A (en) * 2018-07-04 2020-01-16 株式会社安藤・間 Movable measuring device of water content, and rolling compaction method
CN113376709A (en) * 2021-06-21 2021-09-10 西南石油大学 Method for predicting reservoir natural gas hydrate saturation by using logging data
CN113625364A (en) * 2021-08-18 2021-11-09 西南石油大学 A calculation method of shale formation pore pressure based on double correction
CN113820744A (en) * 2020-06-18 2021-12-21 中国石油化工股份有限公司 Reservoir pore pressure prediction method and device, electronic equipment and medium

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001255186A (en) * 2000-03-09 2001-09-21 Matsushita Electric Ind Co Ltd Flow measurement device
JP2011185706A (en) * 2010-03-08 2011-09-22 Central Res Inst Of Electric Power Ind Method and device for monitoring surface defect using laser ultrasonic waves
CN102129084A (en) * 2010-12-17 2011-07-20 中国石油天然气股份有限公司 Method and device for acquiring seismic thin reservoir speed through well control
CN103513270A (en) * 2012-06-21 2014-01-15 中国石油天然气集团公司 Gas reservoir identification evaluating method based on acoustic characteristic of rock and device thereof
FR3026313A1 (en) * 2014-09-29 2016-04-01 Francois Parmentier METHOD OF CHROMATOGRAPHY ON A MULTICAPILLARY TRIM
CN104316978A (en) * 2014-10-29 2015-01-28 中国石油天然气股份有限公司 Near-surface three-dimensional velocity field research method and device for geophysical
CN105804737A (en) * 2016-05-17 2016-07-27 西南石油大学 Method for solving formation porosity on basis of iterative algorithm
WO2019049442A1 (en) * 2017-09-08 2019-03-14 ソニー株式会社 Microparticle measurement device, information processing device, and information processing method
CN109989746A (en) * 2017-12-29 2019-07-09 中国石油天然气股份有限公司 Method and device for evaluating carbonate reservoir
JP2020007730A (en) * 2018-07-04 2020-01-16 株式会社安藤・間 Movable measuring device of water content, and rolling compaction method
CN113820744A (en) * 2020-06-18 2021-12-21 中国石油化工股份有限公司 Reservoir pore pressure prediction method and device, electronic equipment and medium
CN113376709A (en) * 2021-06-21 2021-09-10 西南石油大学 Method for predicting reservoir natural gas hydrate saturation by using logging data
CN113625364A (en) * 2021-08-18 2021-11-09 西南石油大学 A calculation method of shale formation pore pressure based on double correction

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Damage constitutive model of brittle rock considering the compaction of crack;Gu Qingheng;《GEOMECHANICS AND ENGINEERING》;20181228;第15卷(第5期);1081-1089 *
Mathematical methods for evaluating a reservoir based on gas dynamic monitoring during underbalanced drilling;Fan Xiangyu;《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING 》;20160106;第26卷;1068-1079 *
基于体积模型法的孔隙度计算方法研究;王惠宁;《石化技术》;20201228;第27卷(第12期);78-79 *
声速测井标准刻度井设计与建造进展;任晓荣;;《工业计量》;20200917(第05期);46-50 *
平衡深度法估算地层压力-以辽东湾为例;宋海强;《内蒙古石油化工》;20080828(第10期);198-200 *
碳酸盐岩不同含水饱和度对声波传播特性影响;张明明;;《测井技术》;20170220(第01期);18-23 *
致密油储层岩石各向异性的动态声学实验研究;石玉江;《测井技术》;20190620;第43卷(第3期);228-234 *

Also Published As

Publication number Publication date
CN114813509A (en) 2022-07-29

Similar Documents

Publication Publication Date Title
CN105927218B (en) A method and device for predicting gas content in continental mud shale reservoirs
CN105910971B (en) The simultaneous measuring method of rich organic matter compact rock core gas permeability and diffusion coefficient
CN108982817B (en) A method for evaluating gas content in shale based on methane carbon isotope
US7453266B2 (en) Electronic humidity chamber for vapor desorption to determine high capillary pressures
CN105866002B (en) An accurate NMR porosity test method for oil-bearing shale
CN105866009A (en) Method and device for calculating effective porosity of compact oil storage layer
CN104569350B (en) Test method for original water saturation of irregular full-diameter rock samples from sealed coring
CN108446526B (en) Prediction method of shale adsorption capacity based on adsorption characteristic curve and characteristic equation
CN105259080A (en) Shale gas reservoir gas diffusion coefficient experiment test method
CN104182647B (en) The method for obtaining characteristic curve of soil moisture Van Genuchten model parameters
Meng et al. Experimental study on sorption induced strain and permeability evolutions and their implications in the anthracite coalbed methane production
CN110672494B (en) Method for rapidly measuring different porosities of porous concrete
CN107462936A (en) Utilize the method for pressure monitoring Data Inversion low permeability reservoir non-Darcy percolation law
CN104849178A (en) Water vapor diffusion coefficient unsteady state measurement device and method thereof
CN104697914A (en) Predication method of organic holes of different types of clay shale
CN117113884B (en) A method for determining the remaining water saturation around gas wells in water gas reservoirs
CN106019403A (en) Self-generation self-storage hydrocarbon reservoir porosity measurement method
CN114813509B (en) Compaction correction coefficient determination method for calculating rock porosity by using acoustic wave time difference
CN107063935A (en) One kind measurement CO2Diffusion process and the method for diffusion coefficient between water-oil phase
CN110095584B (en) Reservoir oil-water saturation correction method
CN111650108A (en) Method and device for measuring effective porosity of shale rock
CN110702586A (en) Method and device for measuring permeability of rock core
Gross et al. Dynamic pressurization: novel method for measuring fluid permeability
CN108268712A (en) The method and apparatus that nuclear magnetic resonance determines pore media capillary pressure
CN116165116A (en) Prediction method based on compact sandstone elasto-electric property joint inversion pore structure

Legal Events

Date Code Title Description
PB01 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