CN117310812A - Methane fluid longitudinal wave time difference skeleton parameter acquisition method - Google Patents

Methane fluid longitudinal wave time difference skeleton parameter acquisition method Download PDF

Info

Publication number
CN117310812A
CN117310812A CN202311273005.7A CN202311273005A CN117310812A CN 117310812 A CN117310812 A CN 117310812A CN 202311273005 A CN202311273005 A CN 202311273005A CN 117310812 A CN117310812 A CN 117310812A
Authority
CN
China
Prior art keywords
longitudinal wave
time difference
methane
wave time
pressure
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.)
Pending
Application number
CN202311273005.7A
Other languages
Chinese (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.)
Guangdong Ocean University
China France Bohai Geoservices Co Ltd Zhanjiang Branch
Original Assignee
Guangdong Ocean University
China France Bohai Geoservices Co Ltd Zhanjiang Branch
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 Guangdong Ocean University, China France Bohai Geoservices Co Ltd Zhanjiang Branch filed Critical Guangdong Ocean University
Priority to CN202311273005.7A priority Critical patent/CN117310812A/en
Publication of CN117310812A publication Critical patent/CN117310812A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • G01V1/305Travel times
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time

Abstract

The invention discloses a methane fluid longitudinal wave time difference skeleton parameter acquisition method, which comprises the following steps: step one, methane fluid sampling; step two, obtaining a plurality of groups of test experiment data: under different temperature and pressure conditions, performing a longitudinal wave velocity measurement test on the methane fluid sample to obtain corresponding data; step three, constructing a primary calculation model; step four, constructing a longitudinal wave time difference skeleton parameter calculation model; and fifthly, acquiring longitudinal wave time difference skeleton parameters. The method can remarkably improve the accuracy of acquiring the longitudinal wave time difference parameters of the methane fluid skeleton, further improve the accuracy of calculating the high-temperature high-pressure methane gas layer logging interpretation and evaluation porosity by utilizing the acoustic logging data, and has stronger universality.

Description

Methane fluid longitudinal wave time difference skeleton parameter acquisition method
Technical Field
The invention relates to the technical field of oil and gas exploration, in particular to a method for acquiring methane fluid longitudinal wave time difference skeleton parameters.
Background
In the process of methane gas layer logging interpretation and evaluation, logging technicians are required to timely log and interpret the obtained downhole acoustic logging data according to a stratum component volume model so as to obtain porosity parameters of the gas layer, however, the key of the work is to determine methane fluid longitudinal wave time difference skeleton parameters under different stratum temperatures and stratum pressures.
The existing method for acquiring the methane fluid longitudinal wave time difference skeleton parameters is to inquire a chart recorded by related literature, roughly acquire the related parameters from the chart, and take the related parameters as skeleton parameter values for logging interpretation. However, in long production practices, technicians find that the methane fluid longitudinal wave time difference skeleton parameters obtained by the method are low in accuracy, so that the accuracy of porosity parameters obtained by interpretation of logging processing by using acoustic logging data is low. The reason for this is mainly: (1) The actual stratum temperature and pressure parameter change range is larger; (2) Methane gas layers have different buried depth ranges on the well, so that formation pressure and temperature values of the gas layers are inconsistent. In addition, the existing method is relatively high in subjectivity, and some estimation processes are used, so that the empirical value of the methane fluid longitudinal wave time difference skeleton parameter is relatively different from the value in an actual stratum.
Disclosure of Invention
In order to solve the problems, the invention provides a method for calculating the methane longitudinal wave time difference skeleton parameters under the condition of high-resolution stratum environment parameters, so that the problem of low calculation precision of the methane longitudinal wave time difference skeleton parameters between methane gas layer wells and between methane gas layers with different burial depths is effectively solved.
In order to achieve the above purpose, the invention adopts the following technical scheme:
a methane fluid longitudinal wave time difference skeleton parameter acquisition method comprises the following steps:
step one, methane fluid sampling: selecting a methane fluid sample from a methane gas field well site by using a sampling instrument;
step two, obtaining a plurality of groups of test experiment data: under different temperature and pressure conditions, performing a longitudinal wave velocity measurement test on the methane fluid sample to obtain corresponding data;
step three, constructing a primary calculation model by adopting the test experimental data obtained in the step two;
step four, constructing a methane fluid longitudinal wave time difference skeleton parameter calculation model;
step five, acquiring longitudinal wave time difference skeleton parameters: and actually detecting the temperature and the pressure in the high-temperature high-pressure reservoir, substituting the detected temperature value and pressure value into a longitudinal wave time difference skeleton parameter calculation model, and finally obtaining the longitudinal wave time difference skeleton parameter of the methane fluid in the high-temperature high-pressure reservoir.
In the first step, the methane gas sampling instrument is a thermal-insulation pressure-maintaining in-situ fluid sampling instrument.
Further, in the second step, the device used in the measurement test is an in-situ measuring instrument of the velocity of the fluid longitudinal wave.
Further, in the second step, the plurality of sets of measurement test data are 20 sets.
Further, 20 sets of measurement test data were obtained: for the temperature conditions of 20 ℃, 50 ℃, 100 ℃, 150 ℃ and the pressure conditions of 20MPa, 30MPa, 40MPa, 50MPa and 60MPa, any one temperature and any one pressure are selected to be matched to form corresponding environmental conditions, and the longitudinal wave velocity of methane fluid is measured in the environmental conditions, so that 20 groups of test data consisting of longitudinal wave velocity parameters, temperature parameters and pressure parameters are obtained.
In the third step, a data fitting analysis method is used for carrying out parameter fitting analysis on the methane longitudinal wave velocity data, so that a primary calculation model of the longitudinal wave velocity of the methane fluid along with the temperature and pressure parameter changes is obtained.
Further, in the third step, a primary calculation model of the change of the longitudinal wave velocity of the methane fluid with the temperature and pressure parameters is as follows:
V=a*T2+b*T+c
a=-0.00018*P2+0.00142*P-0.0176
b=1.2215*ln(P)-7.1024
c=10.921*P+388.69
v is the methane fluid longitudinal wave velocity obtained by experimental measurement, and the unit is m/s; t is the temperature in degrees Celsius; p is the pressure in MPa.
In the fourth step, the longitudinal wave velocity of the methane fluid is converted into a longitudinal wave time difference value by using the acoustic logging response principle, and a longitudinal wave time difference skeleton parameter calculation model of the methane fluid, wherein the longitudinal wave time difference skeleton parameter varies with temperature and pressure, is obtained.
Further, the longitudinal wave time difference skeleton parameter calculation model is specifically as follows:
ΔT(CH4)=1000000*1/V
wherein DeltaT (CH 4) is a methane fluid longitudinal wave time difference skeleton parameter, and the unit is us/m; v is the methane fluid longitudinal wave velocity in m/s.
Advantageous effects
The accuracy of acquiring the methane fluid skeleton longitudinal wave time difference parameter can be remarkably improved, so that the accuracy of calculating the high-temperature high-pressure methane gas layer logging interpretation and evaluation porosity by utilizing the acoustic logging data is improved, and the method has strong universality. In other words, the invention provides a better and faster acquisition method for the longitudinal wave time difference parameter selection of the methane fluid skeleton in the interpretation and evaluation of high-temperature high-pressure air Tian Cejing.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions of the prior art, the drawings that are necessary for the description of the embodiments or the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention and that other drawings can be obtained according to these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a flow chart of a method of an embodiment of the present application;
FIG. 2 is a plot of the longitudinal wave velocity of methane fluid at various temperature and pressure conditions according to an embodiment of the present application.
Detailed Description
In order to better understand the above technical solutions, the following detailed description will refer to the accompanying drawings and specific embodiments.
The invention provides a methane fluid longitudinal wave time difference frame parameter acquisition method, which can remarkably improve the accuracy of acquiring the methane fluid frame longitudinal wave time difference parameter, further improve the accuracy of calculating the high-temperature high-pressure methane gas layer logging interpretation and evaluation porosity by using acoustic logging data, and has stronger universality. Therefore, the invention provides a better and faster acquisition method for the longitudinal wave time difference parameter selection of the methane fluid skeleton in the explanation and evaluation of the high-temperature high-pressure air Tian Cejing.
As shown in FIG. 1, the method for acquiring the methane fluid longitudinal wave time difference skeleton of the high-temperature high-pressure reservoir comprises the following steps:
step one, methane fluid sampling: and acquiring a methane fluid sample from the high-temperature high-pressure air Tian Jingchang by using a sampling instrument. More specifically, the methane gas sampling instrument is a heat-preserving and pressure-maintaining in-situ fluid sampling instrument.
Step two, obtaining a plurality of groups of test experiment data: laboratory fluid longitudinal wave velocity test analysis is carried out on the methane fluid sample, namely: for the temperature of 20 ℃, 50 ℃, 100 ℃ and 150 ℃ and the pressure of 20MPa, 30MPa, 40MPa, 50MPa and 60MPa, any one temperature and any one pressure are selected to be matched to form corresponding environmental conditions, the methane fluid longitudinal wave velocity is measured under the environmental conditions, 20 groups of methane fluid longitudinal wave velocity experimental data are obtained, the 20 groups of methane fluid longitudinal wave velocity experimental data are reflected in a graph mode, and a data scatter diagram of the methane fluid longitudinal wave velocity under different temperature and pressure conditions is obtained, as shown in figure 2. As can be seen from fig. 2, the methane fluid longitudinal wave velocity and the temperature and pressure show binary correlation, under the same temperature condition, the methane fluid longitudinal wave velocity increases with the increase of the pressure, and under the same pressure condition, the methane fluid longitudinal wave velocity decreases with the increase of the temperature, that is, the methane fluid longitudinal wave velocity has different values under different temperature and pressure conditions. More specifically, in this assay, the apparatus used is a fluid longitudinal wave velocity in situ meter.
Step three, constructing a primary calculation model: and carrying out parameter fitting analysis on the plurality of groups of data by using a data fitting analysis method in a mathematical statistics category, thereby obtaining a primary calculation model of the longitudinal wave velocity of the methane fluid along with the temperature and pressure parameter change. The method comprises the following steps:
in order to better and more accurately express the change rule between the methane fluid longitudinal wave speed and the temperature and pressure, adopting optimized data fitting analysis, firstly, processing and analyzing the temperature and pressure parameters, and obtaining a calculation relational expression between the methane fluid longitudinal wave speed and the temperature and pressure by using binary polynomial fitting analysis:
V=a*T2+b*T+c
a=-0.00018*P2+0.00142*P-0.0176
b=1.2215*ln(P)-7.1024
c=10.921*P+388.69
v is the methane fluid longitudinal wave velocity obtained by experimental measurement, and the unit is m/s; t is the temperature in degrees Celsius; p is the pressure in MPa.
Step four, constructing a longitudinal wave time difference skeleton parameter calculation model: based on the primary calculation model, the longitudinal wave velocity of the methane fluid is converted into a longitudinal wave time difference value by utilizing the acoustic logging response principle, and the longitudinal wave time difference skeleton parameter calculation model of the methane fluid of the high-temperature and high-pressure reservoir, wherein the longitudinal wave time difference skeleton parameter of the methane fluid changes along with the temperature and the pressure, is obtained. The specific implementation mode is as follows:
from the acoustic logging response principle, the acoustic logging measures the time taken by the sliding wave to traverse the unit length of the stratum, the instant difference is in us/m, the time difference deltat of the acoustic logging record is only related to the stratum speed, and the time required by converting the time difference into the unit distance of the acoustic wave is called the time difference through instrument scales, namely:
Δt=1/ν (1)
substituting the methane fluid longitudinal wave velocity and temperature and pressure parameter calculation model obtained in the step three into the formula (1) to obtain a longitudinal wave time difference mathematical expression of methane under stratum conditions, wherein the longitudinal wave time difference mathematical expression is as follows:
wherein DeltaT (CH 4) is a methane fluid longitudinal wave time difference skeleton parameter, and the unit is us/m; t is the temperature in degrees Celsius; p is the pressure in MPa.
Step five, acquiring longitudinal wave time difference skeleton parameters: and (3) actually detecting the temperature and the pressure of the methane gas layer, substituting the detected temperature value and pressure value into the longitudinal wave time difference skeleton parameter calculation model, and finally obtaining the longitudinal wave time difference skeleton parameter of the methane fluid in the high-temperature high-pressure reservoir.
The present invention is not limited to the preferred embodiments, but is capable of modification and variation in detail, and other modifications and variations can be made by those skilled in the art without departing from the scope of the present invention.

Claims (9)

1. The methane fluid longitudinal wave time difference skeleton parameter acquisition method is characterized by comprising the following steps of:
step one, methane fluid sampling: selecting a methane fluid sample from a methane gas field well site by using a sampling instrument;
step two, obtaining a plurality of groups of test experiment data: under different temperature and pressure conditions, performing a longitudinal wave velocity measurement test on the methane fluid sample to obtain corresponding data;
step three, constructing a primary calculation model by adopting the test experimental data obtained in the step two;
step four, constructing a methane fluid longitudinal wave time difference skeleton parameter calculation model;
step five, acquiring longitudinal wave time difference skeleton parameters: and actually detecting the temperature and the pressure in the high-temperature high-pressure reservoir, substituting the detected temperature value and pressure value into a longitudinal wave time difference skeleton parameter calculation model, and finally obtaining the longitudinal wave time difference skeleton parameter of the methane fluid in the high-temperature high-pressure reservoir.
2. The method for obtaining parameters of a methane fluid longitudinal wave time difference skeleton according to claim 1, wherein in the first step, the methane gas sampling instrument is a thermal insulation pressure-maintaining in-situ fluid sampling instrument.
3. The method for obtaining the methane fluid longitudinal wave time difference skeleton parameter according to claim 1, wherein in the second step, the equipment used in the measurement test is an in-situ measuring instrument of the fluid longitudinal wave velocity.
4. The method for obtaining methane fluid longitudinal wave time difference skeleton parameters according to any one of claims 1-3, wherein in the second step, the plurality of sets of measurement test data are 20 sets.
5. The method for acquiring methane fluid longitudinal wave time difference skeleton parameters according to claim 4, wherein the 20 sets of acquired measurement test data are: for the temperature conditions of 20 ℃, 50 ℃, 100 ℃, 150 ℃ and the pressure conditions of 20MPa, 30MPa, 40MPa, 50MPa and 60MPa, any one temperature and any one pressure are selected to be matched to form corresponding environmental conditions, and the longitudinal wave velocity of methane fluid is measured in the environmental conditions, so that 20 groups of test data consisting of longitudinal wave velocity parameters, temperature parameters and pressure parameters are obtained.
6. The method for obtaining methane fluid longitudinal wave time difference skeleton parameters according to claim 5, wherein in the third step, a data fitting analysis method is used to perform parameter fitting analysis on each methane longitudinal wave velocity data, so as to obtain a primary calculation model of the methane fluid longitudinal wave velocity along with the temperature and pressure parameter changes.
7. The method for obtaining the longitudinal wave time difference skeleton parameter of the methane fluid according to claim 6, wherein in the third step, a primary calculation model of the longitudinal wave velocity of the methane fluid according to the temperature and pressure parameters is as follows:
V=a*T 2 +b*T+c
a=-0.00018*P 2 +0.00142*P -0.0176
b=1.2215*ln(P)-7.1024
c=10.921*P+388.69
v is the methane fluid longitudinal wave velocity obtained by experimental measurement, and the unit is m/s; t is the temperature in degrees Celsius; p is the pressure in MPa.
8. The method for obtaining the longitudinal wave time difference skeleton parameter of the methane fluid according to claim 6, wherein in the fourth step, the longitudinal wave velocity of the methane fluid is converted into a longitudinal wave time difference value by using the acoustic logging response principle, and a longitudinal wave time difference skeleton parameter calculation model of the longitudinal wave time difference skeleton parameter of the methane fluid, which changes along with the temperature and the pressure, is obtained.
9. The method for obtaining the longitudinal wave time difference skeleton parameter of the methane fluid according to claim 8, wherein the longitudinal wave time difference skeleton parameter calculation model is specifically as follows:
ΔT (CH4) =1000000*1/V
wherein DeltaT (CH4) The unit is us/m for the methane fluid longitudinal wave time difference skeleton parameter; v is the methane fluid longitudinal wave velocity in m/s.
CN202311273005.7A 2023-09-27 2023-09-27 Methane fluid longitudinal wave time difference skeleton parameter acquisition method Pending CN117310812A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311273005.7A CN117310812A (en) 2023-09-27 2023-09-27 Methane fluid longitudinal wave time difference skeleton parameter acquisition method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311273005.7A CN117310812A (en) 2023-09-27 2023-09-27 Methane fluid longitudinal wave time difference skeleton parameter acquisition method

Publications (1)

Publication Number Publication Date
CN117310812A true CN117310812A (en) 2023-12-29

Family

ID=89286207

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311273005.7A Pending CN117310812A (en) 2023-09-27 2023-09-27 Methane fluid longitudinal wave time difference skeleton parameter acquisition method

Country Status (1)

Country Link
CN (1) CN117310812A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105003258A (en) * 2015-08-07 2015-10-28 中国海洋石油总公司 Method for acquiring density framework parameters of methane fluid in high temperature high pressure air layer
CN105089632A (en) * 2015-08-04 2015-11-25 中国海洋石油总公司 Method for obtaining CO2 fluid longitudinal wave time difference framework parameters of high-temperature and high-pressure reservoir
CN105114068A (en) * 2015-09-07 2015-12-02 中国地质大学(北京) Method of predicting high-water-yield area in coalbed methane area via logging information
CN105952446A (en) * 2016-04-26 2016-09-21 中国海洋石油总公司 Measurement method for component content of petroleum and natural gas reservoir
WO2017024530A1 (en) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 Method for calculating content of organic carbon in hydrocarbon source rock
JP2020106344A (en) * 2018-12-27 2020-07-09 日立Geニュークリア・エナジー株式会社 Calculation method and calculation device of gas sound speed
CN112664186A (en) * 2020-12-25 2021-04-16 中法渤海地质服务有限公司 Method for early warning and identifying other high-voltage source by using methane carbon isotope data

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089632A (en) * 2015-08-04 2015-11-25 中国海洋石油总公司 Method for obtaining CO2 fluid longitudinal wave time difference framework parameters of high-temperature and high-pressure reservoir
CN105003258A (en) * 2015-08-07 2015-10-28 中国海洋石油总公司 Method for acquiring density framework parameters of methane fluid in high temperature high pressure air layer
WO2017024530A1 (en) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 Method for calculating content of organic carbon in hydrocarbon source rock
CN105114068A (en) * 2015-09-07 2015-12-02 中国地质大学(北京) Method of predicting high-water-yield area in coalbed methane area via logging information
CN105952446A (en) * 2016-04-26 2016-09-21 中国海洋石油总公司 Measurement method for component content of petroleum and natural gas reservoir
JP2020106344A (en) * 2018-12-27 2020-07-09 日立Geニュークリア・エナジー株式会社 Calculation method and calculation device of gas sound speed
CN112664186A (en) * 2020-12-25 2021-04-16 中法渤海地质服务有限公司 Method for early warning and identifying other high-voltage source by using methane carbon isotope data

Similar Documents

Publication Publication Date Title
Sharma et al. Low-frequency distributed acoustic sensing for early gas detection in a wellbore
CN103926184B (en) Rock core gas surveys porosity detection method and detection device thereof
CN108694264B (en) Method for determining permeability of shale gas reservoir
CN102220865B (en) Method for detecting limestone formation pore pressure
CN110296931B (en) Characterization method and system for oil-water relative permeability information of tight sandstone
WO2016179593A1 (en) Full immersion pressure-pulse decay
CN109138975B (en) Novel method for solving phase-permeability characteristic curve based on time-shifting logging data
CN105863626A (en) Evaluation method for physical and chemical action of drilling fluid and shale formation
CN109061099B (en) Nondestructive experimental evaluation method for damage degree of heterogeneous compact rock
CN103321634A (en) Method for calculating stratum pressure by aid of liquid level process
WO2009102526A2 (en) Methods for the identification of bubble point pressure
CN110043254B (en) Method for obtaining stratum effective permeability based on cable stratum test data
CN105003258B (en) A kind of acquisition methods of HTHP gas-bearing formation methane stream density matrix parameter
Ling et al. Determining the permeability of tight rock with gas transient flow
CN105952446A (en) Measurement method for component content of petroleum and natural gas reservoir
CN112799140B (en) Permeability estimation method based on natural potential inversion
CN108664677B (en) Oil and gas well production data analysis method
CN105089632A (en) Method for obtaining CO2 fluid longitudinal wave time difference framework parameters of high-temperature and high-pressure reservoir
CN117310812A (en) Methane fluid longitudinal wave time difference skeleton parameter acquisition method
CN112268923B (en) Method for acquiring formation thermal conductivity based on logging curve
Ling et al. A method to determine pore compressibility based on permeability measurements
Popov et al. Evolution in information on crustal geothermal parameters due to application of advanced experimental basis
Anwar et al. Detecting and characterizing fluid leakage through wellbore flaws using fiber-optic distributed acoustic sensing
CN110965992B (en) Method for determining viscosity of stratum gas-containing crude oil
CN108647417B (en) Simple method for determining gas saturation of shale gas reservoir

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