CN104948164B - The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter - Google Patents
The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter Download PDFInfo
- Publication number
- CN104948164B CN104948164B CN201510223616.XA CN201510223616A CN104948164B CN 104948164 B CN104948164 B CN 104948164B CN 201510223616 A CN201510223616 A CN 201510223616A CN 104948164 B CN104948164 B CN 104948164B
- Authority
- CN
- China
- Prior art keywords
- density
- matrix parameter
- density matrix
- fluid
- carbon dioxide
- 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
- Carbon And Carbon Compounds (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter, comprise the following steps：Carbon dioxide stream sampler body；Obtain multigroup test experiments data；Build primary calculations model：Using data fitting analyzing method, parameter fitting analysis is carried out to above-mentioned multi-group data, obtain primary calculations model；Build density matrix parameter computation model：On the basis of described primary calculations model, using density log response theory, the bulk density of CO 2 fluid is converted into apparent bulk density value, obtains density matrix parameter computation model；The temperature value and pressure value that will be measured, are updated to density matrix parameter computation model, obtain the density matrix parameter of CO 2 fluid in HTHP reservoir.The degree of accuracy for obtaining CO 2 fluid skeletal density parameter can be significantly improved using the method, and then improves the precision that HTHP evaluates porosity rich in carbon dioxide gas-bearing formation well log interpretation, with stronger versatility.
Description
Technical field
The present invention relates in a kind of oil-gas exploration technology, HTHP is rich in the well logging that carbon dioxide gas-bearing formation porosity is explained
Assessment technique, specifically refers to the acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter.
Background technology
, it is necessary to logging technique personnel are to acquisition in HTHP is rich in carbon dioxide gas-bearing formation well log interpretation evaluation procedure
Density log data in underground carries out timely well log interpretation according to stratum component volume-based model, obtains the porosity parameter of gas-bearing formation,
It is critical only that for work determines carbon dioxide stream volume density matrix parameter under the conditions of stratum.Existing carbon dioxide stream volume density bone
The acquisition methods of frame parameter are generally region empirical method, and the method mainly uses the carbon dioxide gas sample that underground is got to carry out PVT
Experiment, then obtains the density value under certain temperature and pressure, as the matrix parameter empirical value in region, but for a long time
Production practices in technical staff find the above method obtain carbon dioxide stream volume density matrix parameter accuracy it is relatively low, study carefully
Its pressure coefficient excursion of the main HTHP reservoir of its reason is larger, and HTHP is rich in carbon dioxide gas-bearing formation between well
Buried depth scope is also different on well, causes the strata pressure of gas-bearing formation, temperature value inconsistent；The subjectivity of other the method compared with
Greatly, the process of some estimations has been used；So that in carbon dioxide stream volume density matrix parameter empirical value and actual formation
Value difference is larger.
The content of the invention
It is an object of the invention to provide a kind of highly versatile and accuracy HTHP reservoir CO 2 fluid high is close
Spend the acquisition methods of matrix parameter.
To achieve the above object, the technical solution adopted in the present invention is：
A kind of acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter, comprise the following steps：
Step 1, carbon dioxide stream sampler body：Using sample apparatus CO 2 fluid is chosen from HTHP gas field well site
Sample；
Step 2, obtains multigroup test experiments data：Under different temperatures and pressure condition, to CO 2 fluid sample
The determination test of bulk density is carried out, the multi-group data being made up of bulk density parameter, temperature parameter and pressure parameter is obtained；
Step 3, builds primary calculations model：Using the data fitting analyzing method in mathematical statistics category, to above-mentioned many
Group data carry out parameter fitting analysis, obtain the bulk density of CO 2 fluid with temperature and the primary calculations mould of pressure change
Type；
Step 4, builds density matrix parameter computation model：On the basis of described primary calculations model, using density
Log response principle, apparent bulk density value is converted into by the bulk density of CO 2 fluid, obtains the two of HTHP reservoir
The density matrix parameter of carbon oxide fluid is with temperature and the density matrix parameter computation model of pressure change；
Step 5, obtains density matrix parameter：Temperature and pressure in HTHP reservoir is carried out actually detected, will surveyed
The temperature value and pressure value for obtaining, are updated in density matrix parameter computation model, obtain carbon dioxide stream in HTHP reservoir
The density matrix parameter of body.
Further, in the step 2, in determination test, equipment used is fluid volume penetron.
Further, in the step 2, multi-group data is 20 groups, wherein, the selected value of temperature is 20 DEG C, 50 DEG C, 100
DEG C and 150 DEG C, the selected value of pressure is 20MPa, 30MPa, 40MPa, 50MPa and 60MPa.
Further, in the step 3, primary calculations model is specific as follows：
Y=0.28758*ln (X)+1.0122, wherein, Y is the CO 2 fluid bulk density that experiment measurement is obtained, single
Position is (g/cm3), X be pressure value divided by temperature value, unit is (MPa/ DEG C).
Further, in the step 4, the density matrix parameter computation model is specific as follows：
ρa(CO2)=1.0697* (0.2875*ln (x)+1.0122) -0.1883,
Wherein ρa(CO2) it is HTHP reservoir carbon dioxide stream volume density matrix parameter, unit is (g/cm3), and X is pressure
Divided by temperature value, unit is (MPa/ DEG C) to force value.
Beneficial effects of the present invention：Proved by production practices, acquisition carbon dioxide can be significantly improved using the method
The degree of accuracy of fluid skeletal density parameter, and then HTHP is improved rich in carbon dioxide gas-bearing formation well log interpretation evaluation porosity
Precision, with stronger versatility.The present invention is CO 2 fluid skeletal density in the evaluation of HTHP gas field well log interpretation
Parameter is chosen and provides more preferable more efficiently acquisition methods.
Brief description of the drawings
The present invention is described in further detail with reference to the accompanying drawings and detailed description：
Fig. 1 is flow chart of the invention；
Fig. 2 is that the CO 2 fluid volume for obtaining multigroup all formation of test experiments data in step 2 of the invention is close
Spend the scatter plot of data under different temperatures, pressure condition；
Fig. 3 is present invention result and CO 2 fluid bulk density relation of the pressure value divided by temperature value in step 3
Scatter plot of data.
Specific embodiment
As shown in figure 1, a kind of acquisition methods of HTHP reservoir CO2 fluid density skeletons, it comprises the following steps：
Step 1, carbon dioxide stream sampler body：CO 2 fluid sample is chosen：Using sample apparatus from HTHP gas
Field well site collects CO 2 fluid sample；
Step 2, obtains multigroup test experiments data：Common laboratory fluids bulk density examination is carried out to CO 2 fluid sample
Test analysis, respectively temperature be 20 DEG C, 50 DEG C, 100 DEG C and 150 DEG C, pressure be 20MPa, 30MPa, 40MPa, 50MPa with
CO 2 fluid bulk density is determined under the conditions of 60MPa, measurement altogether obtains 20 groups of CO 2 fluid bulk density experiment numbers
According to, this 20 groups of CO 2 fluid bulk density experimental datas are embodied in the way of scheming, obtain CO 2 fluid body
Scatter plot of data of the product density under different temperatures, pressure condition, as shown in Figure 2.Can through the scatter plot of data obtained by experiment
Know, what CO 2 fluid bulk density and temperature, pressure were presented is dyadic correlation relation, under the conditions of same temperature, with
The increase of pressure, the increase of CO 2 fluid bulk density, and under same pressure condition, with the rising of temperature, titanium dioxide
Carbon flow body bulk density reduction, that is to say, that CO 2 fluid bulk density has different under different temperatures pressure condition
Numerical value.
Step 3, builds primary calculations model：Using the data fitting analyzing method in mathematical statistics category, to above-mentioned many
Group data carry out parameter fitting analysis, obtain the bulk density of CO 2 fluid with temperature and the primary calculations mould of pressure change
Type.
Comprise the following steps that：In order to more preferably, more accurately state between CO 2 fluid bulk density and temperature, pressure
Changing Pattern, take optimization data Fitting Analysis, Treatment Analysis are carried out to temperature, pressure parameter first, by pressure data values
Divided by temperature data value, new parameter X is constructed, obtain the scatter diagram with CO 2 fluid bulk density and X as coordinate.
In a particular embodiment, experimental data is processed, with X as abscissa, CO 2 fluid bulk density be vertical seat
Mark sets up scatter diagram, and CO 2 fluid bulk density data and X parameter are set up into logarithmic function relational expression Y=0.28758*ln
(X)+1.0122, as shown in figure 3, wherein, Y is the CO 2 fluid bulk density that experiment measurement is obtained, unit is (g/
Cm3), X be pressure value divided by temperature value, unit is (MPa/ DEG C)；In figure, R2 represents CO 2 fluid bulk density and pressure
The coefficient correlation of power/temperature value square, the value is bigger to represent CO 2 fluid bulk density and Pressure/Temperature Value Data
Correlation is better.
Step 4, builds density matrix parameter computation model：On the basis of described primary calculations model, using density
Log response principle, apparent bulk density value is converted into by the bulk density of CO 2 fluid, obtains the two of HTHP reservoir
The density matrix parameter of carbon oxide fluid is with temperature and the density matrix parameter computation model of pressure change.
Its specific implementation step is as follows：From density log response theory, what density log was measured is the electronics on stratum
Density pe, the scale equation shown with aqueous limestone graduation apparatus is：
ρa=1.07 (ρe)i-0.1883 (1)
Wherein (ρe)iIt is electron density index, is defined as：
In formula, N is Avogadro's number；
By after scale equation scale, that density logging instrument is recorded is the apparent bulk density ρ on stratuma, the volume on stratum
The relation of density and electron density：
In formula, N is Avogadro's number, and Z is atomic number, and A is atomic weight；
(3) and (2) are substituted into formula (1).For the compound that polyatom is constituted, the apparent bulk density after its scale is：
The molecular formula of carbon dioxide is CO2, and the atomic of C, O and atomic weight substituted into formula (4), and by step 103 institute
CO 2 fluid bulk density and temperature, pressure parameter computation model substitutes into formula (4) and can obtain CO2 under formation condition
Apparent bulk density mathematic(al) representation is：
ρa(CO2)=1.0697* (0.2875*ln (x)+1.0122) -0.1883 (5)
Wherein ρa(CO2) it is HTHP reservoir carbon dioxide stream volume density matrix parameter, unit is (g/cm3), and X is pressure
Divided by temperature value, unit is (MPa/ DEG C) to force value.
Step 5, obtains density matrix parameter：Temperature and pressure in HTHP reservoir is carried out actually detected, will surveyed
The temperature value and pressure value for obtaining, are updated in density matrix parameter computation model, obtain carbon dioxide stream in HTHP reservoir
The density matrix parameter of body.
Above-described specific embodiment, has been carried out further to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail, should be understood that and the foregoing is only specific embodiment of the invention, be not intended to limit the present invention
Protection domain, all any modification, equivalent substitution and improvements within the spirit and principles in the present invention, done etc. all should include
Within protection scope of the present invention.The content that this specification is not described in detail belongs to known to professional and technical personnel in the field
Prior art.
Claims (3)
1. a kind of acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter, it is characterised in that including as follows
Step：Step 1, carbon dioxide stream sampler body：Using sample apparatus CO 2 fluid sample is chosen from HTHP gas field well site
Product；Step 2, obtains multigroup test experiments data：Under different temperatures and pressure condition, body is carried out to CO 2 fluid sample
The determination test of product density, obtains the multi-group data being made up of bulk density parameter, temperature parameter and pressure parameter；Step 3, structure
Build primary calculations model：Using the data fitting analyzing method in mathematical statistics category, line parameter plan is entered to above-mentioned multi-group data
Analysis is closed, the bulk density of CO 2 fluid is obtained with temperature and the primary calculations model of pressure change, primary calculations model
It is specific as follows：Y=0.28758*ln (X)+1.0122, wherein, Y is the CO 2 fluid bulk density that experiment measurement is obtained, single
Position is g/cm3, X be pressure value divided by temperature value, unit is MPa/ DEG C；Step 4, builds density matrix parameter computation model：
On the basis of described primary calculations model, using density log response theory, the bulk density of CO 2 fluid is converted
It is apparent bulk density value, obtains the density matrix parameter of CO 2 fluid of HTHP reservoir with temperature and pressure change
Density matrix parameter computation model；Obtain density matrix parameter：Actual inspection is carried out to the temperature and pressure in HTHP reservoir
Survey, the temperature value and pressure value that will be measured are updated in density matrix parameter computation model, obtain dioxy in HTHP reservoir
Change the density matrix parameter of carbon flow body, the density matrix parameter computation model is specific as follows：
, wherein
Be HTHP reservoir carbon dioxide stream volume density matrix parameter, unit is g/cm3, X be pressure value divided by temperature value, unit is
MPa/℃。
2. acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter according to claim 1, it is special
Levy and be：In the step 2, in determination test, equipment used is fluid volume penetron.
3. acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter according to claim 1, it is special
Levy and be：In the step 2, multi-group data is 20 groups, wherein, the selected value of temperature is 20 DEG C, 50 DEG C, 100 DEG C and 150 DEG C,
The selected value of pressure is 20MPa, 30 MPa, 40MPa, 50MPa and 60MPa.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510223616.XA CN104948164B (en) | 2015-05-05 | 2015-05-05 | The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510223616.XA CN104948164B (en) | 2015-05-05 | 2015-05-05 | The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104948164A CN104948164A (en) | 2015-09-30 |
| CN104948164B true CN104948164B (en) | 2017-06-16 |
Family
ID=54163172
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510223616.XA Active CN104948164B (en) | 2015-05-05 | 2015-05-05 | The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104948164B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105952446B (en) * | 2016-04-26 | 2019-03-01 | 中国海洋石油集团有限公司 | A kind of petroleum gas gas-bearing formation component content measuring method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101963056A (en) * | 2010-08-19 | 2011-02-02 | 中国石油大学(北京) | Method for predicting carbonate formation pore pressure by using log information |
| WO2012006871A1 (en) * | 2010-07-16 | 2012-01-19 | 陈再迁 | Method of prospecting directly for free gas in stratum |
| CN102748004A (en) * | 2012-06-06 | 2012-10-24 | 中国石油化工股份有限公司 | Method for determining hydrocarbon saturation of gas logging |
| CN104213899A (en) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | Logging identifying method of stratum rock framework |
| CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101197211B1 (en) * | 2010-11-29 | 2012-11-02 | 한국지질자원연구원 | System and method for controlling temperature of fluid for improving injectivity in stratum of supercritical carbon dioxide |
-
2015
- 2015-05-05 CN CN201510223616.XA patent/CN104948164B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012006871A1 (en) * | 2010-07-16 | 2012-01-19 | 陈再迁 | Method of prospecting directly for free gas in stratum |
| CN101963056A (en) * | 2010-08-19 | 2011-02-02 | 中国石油大学(北京) | Method for predicting carbonate formation pore pressure by using log information |
| CN102748004A (en) * | 2012-06-06 | 2012-10-24 | 中国石油化工股份有限公司 | Method for determining hydrocarbon saturation of gas logging |
| CN104213899A (en) * | 2013-06-04 | 2014-12-17 | 中国石油化工股份有限公司 | Logging identifying method of stratum rock framework |
| CN104573344A (en) * | 2014-12-25 | 2015-04-29 | 中国海洋石油总公司 | Method for acquiring gas content of shale reservoir through well logging data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104948164A (en) | 2015-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103792338B (en) | A kind of determination method and apparatus of content of organic carbon of hydrocarbon source rock | |
| Lu et al. | Determination of diffusion coefficients of carbon dioxide in water between 268 and 473 K in a high-pressure capillary optical cell with in situ Raman spectroscopic measurements | |
| Dong et al. | Fracture identification and evaluation using conventional logs in tight sandstones: A case study in the Ordos Basin, China | |
| CN104298883A (en) | Establishment method for hydrocarbon source rock hydrocarbon production rate charts in petroleum resource assessment | |
| CN103344541A (en) | Method for measuring total porosity of shale | |
| CN105003258B (en) | A kind of acquisition methods of HTHP gas-bearing formation methane stream density matrix parameter | |
| CN109060928A (en) | Quantitatively judge the method and apparatus of carbonate rock cement formation and hydrocarbon charge timing | |
| CN105651912A (en) | Rock pyrologger and pyrolytic analysis method | |
| Liu et al. | Experimental investigation of pore structure and movable fluid traits in tight sandstone | |
| CN106525681A (en) | Method for determining pore diameter of shale reservoir | |
| CN109061099B (en) | Nondestructive experimental evaluation method for damage degree of heterogeneous compact rock | |
| Szlązak et al. | Estimation of gas loss in methodology for determining methane content of coal seams | |
| Mukhametdinova et al. | Reservoir properties of low-permeable carbonate rocks: experimental features | |
| CN106501149B (en) | A kind of porosity calculation method based on the analysis of coal petrography industrial components | |
| Yang et al. | The characteristics of oil migration due to water imbibition in tight oil reservoirs | |
| CN105089632B (en) | A kind of HTHP reservoir CO2The acquisition methods of fluid compressional wave time difference matrix parameter | |
| CN104948164B (en) | The acquisition methods of HTHP reservoir carbon dioxide stream volume density matrix parameter | |
| Bobrova et al. | Laboratory investigation of hydraulic fracture behavior of unconventional reservoir rocks | |
| CN105952446A (en) | Measurement method for component content of petroleum and natural gas reservoir | |
| Zhou et al. | An evaluation method of brittleness characteristics of shale based on the unloading experiment | |
| CN106248444A (en) | Soil respiration CO2the collection of gas and carbon stable isotope value method of testing thereof | |
| CN106706887A (en) | Measuring method of solvent concentration in thickened oil | |
| Guo et al. | Saturation Determination and Fluid Identification in Carbonate Rocks Based on Well Logging Data: A Middle Eastern Case Study | |
| Zhan et al. | Fractal Characteristics of Deep Shales in Southern China by Small-Angle Neutron Scattering and Low-Pressure Nitrogen Adsorption | |
| Wałowski | Gas permeability model for porous materials from underground coal gasification technology |
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 | ||
| CP03 | Change of name, title or address |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee after: CNOOC (China) Limited Zhanjiang Branch Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 China oil tower, 25 Chaoyangmen North Street, Dongcheng District, Beijing Co-patentee before: CNOOC (China) Limited Zhanjiang Branch Patentee before: China National Offshore Oil Corporation |
|
| CP03 | Change of name, title or address |