CN108548765A - A kind of porosity calculation method for the soil skeleton that becomes sticky - Google Patents

A kind of porosity calculation method for the soil skeleton that becomes sticky Download PDF

Info

Publication number
CN108548765A
CN108548765A CN201810260747.9A CN201810260747A CN108548765A CN 108548765 A CN108548765 A CN 108548765A CN 201810260747 A CN201810260747 A CN 201810260747A CN 108548765 A CN108548765 A CN 108548765A
Authority
CN
China
Prior art keywords
density
mineral
clay
porosity
skeletal
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.)
Granted
Application number
CN201810260747.9A
Other languages
Chinese (zh)
Other versions
CN108548765B (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.)
China Petroleum and Natural Gas Co Ltd
China Petroleum Logging Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
China Petroleum Logging Co Ltd
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 China Petroleum and Natural Gas Co Ltd, China Petroleum Logging Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201810260747.9A priority Critical patent/CN108548765B/en
Publication of CN108548765A publication Critical patent/CN108548765A/en
Application granted granted Critical
Publication of CN108548765B publication Critical patent/CN108548765B/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

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention discloses a kind of porosity calculation methods for the soil skeleton that becomes sticky.This method includes:Physical property experiment is carried out to rock core first, X diffraction total rock mineral are tested, acquisition core porosity, mineral constituent, mineral content data;Then the scanning well logging of X diffraction total rock mineral experimental data scale lithology is utilized, it is accurate to calculate mineral section;Obtained core porosity, mineral section are combined inverting with density curve, mixed matrix density computation model, density porosity volume-based model and obtain clay skeletal density value;Obtained clay skeletal density value is obtained clay skeletal density computation model with spectrum curve, density curve combination inverting;Clay skeletal density computation model, mineral section, mixed matrix density are calculated into models coupling inverting mixed matrix density value again;Obtained mixed matrix density value combination density porosity volume-based model inverting is finally obtained into porosity.

Description

A kind of porosity calculation method for the soil skeleton that becomes sticky
Technical field
The invention belongs to logging technique field, more particularly to a kind of porosity calculation method for the soil skeleton that becomes sticky.
Background technology
With the continuous development of exploration engineering and geological theory, is explored in addition to conventional gas and oil and constantly obtain new discovery, it is complicated The exploration of lithology, low porosity and low permeability petroleum resources is also obtained compared with quantum jump.In recent years complex lithology oil-gas reservoir is paid attention to by various countries, I The complex lithology reservoir exploration of state also achieves progress, the complex lithology reservoirs such as Lacustrine Carbonates, the basement rock of some blocks Probing exploitation has obtained high yield commercial hydrocarbon flow.Complex lithology reservoir mineralogical composition and content are changeable, result in its rock object Characteristic is managed, if skeletal density, skeleton neutron, the variation of skeleton sound wave are complicated, porosity is difficult to calculate accurately.
The common method for seeking porosity has following three kinds:The first is nuclear magnetic resonance log method.This method is to pass through Observation to proton signal in formation pore fluid is obtained for the porosity needed for formation evaluation.Second method is area Empirical statistics equation, this method utilize core analysis porosity scale log, establish regional experience porosity calculation mould Type;The third method is that different reservoir is given when handling Reservoir Section based on well logging volume physical model porosity calculation method Rock matrix value carrys out inverting reservoir porosity.In the December, 2008 such as open non-patent literature such as Tan Fengqi exist《Current Foreign Well Logging Technology》 On deliver《Application of the Element Capture Spectroscopy Logging in volcanic rock reservoir porosity calculation》, it is bent that this method is based on ECS sensitive elements Line establishes rock core skeletal density model, obtains the matrix density curve with change in depth, then inverted porosity curve.
Porosity calculation method set forth above has respective benefit and limitation.Nuclear magnetic resonance log method can be direct Continuous porosity curve is obtained, but this method cost is higher, and in high salinity stratum and low porosity and low permeability stratum nuclear magnetic signal Signal-to-noise ratio it is low, measurement accuracy is low.Regional experience statistical formula method is simple, easy to operate, but various in complex lithology mineralogical composition Stratum, log is influenced big by rock matrix, and the relative error for calculating porosity is larger.Based on well logging volume physics mould Type hole porosity computational methods applicability is wide, but in the reservoir of mineral constituent complexity, clay mineral type is more, and different buried depth is at Tibetan item The aqueous degree difference of clay is big under part, and the selection of clay skeleton value is the difficult point that it is faced.
Invention content
The purpose of the present invention is to provide a kind of porosity calculation methods for the soil skeleton that becomes sticky, to solve the above problems.
To achieve the above object, the present invention uses following technical scheme:
A kind of porosity calculation method for the soil skeleton that becomes sticky, which is characterized in that include the following steps:
Step 1:To rock core carry out physical property experiment, X diffraction total rock mineral experiment, obtain core porosity, mineral constituent, Mineral content data;
Step 2:Well logging is scanned with the X diffraction total rock mineral experimental data scale lithology that step 1 obtains, accurately calculates mine Object section;
Step 3:Core porosity, mineral section and density curve, the mixed matrix obtained according to step 1 and step 2 Density computation model, density porosity volume-based model are combined inverting and obtain clay skeletal density value;
Step 4:The clay skeletal density value that step 3 obtains is glued with spectrum curve, density curve combination inverting Soil skeleton density computation model;
Step 5:By mineral section that clay skeletal density computation model that step 4 obtains, step 2 obtain with mix Skeletal density computation model is combined inverting and obtains mixed matrix density value;
Step 6:The mixed matrix density value combination density porosity volume-based model inverting that step 5 obtains is obtained into hole Degree.
Further, well logging is scanned with the X diffraction total rock mineral experimental data scale lithology that step 1 obtains, obtained with depth Spend the continuous mineral section of variation.
Further, the core porosity that is obtained according to step 1 and step 2, mixes bone at mineral section and density curve Frame density computation model, density porosity volume-based model are combined inverting and obtain clay skeletal density value;It specifically will be according to as follows Formula carries out:
Wherein, ρmaFor mixed matrix density, ρshIndicate clay skeletal density, vshIndicate clay mineral content, ρmaiIt indicates The density skeleton of i-th kind of mineral and include clay mineral, vmaiIt indicates i-th kind of mineral content and does not include clay mineral, Φ tables Show core porosity, ρbIndicate log measurement density value, ρfIndicate the density of water flooding, n > 1;
Known parameters ρ in formulamai、ρfIt is obtained by common skeletal density parameter list, vmai、vshWell logging determination is scanned by lithology Mineral section obtains, and Φ is obtained by rock core physical property experiment data, ρbIt is obtained by density curve, solves equation group and obtain ρshValue.
Further, the clay skeletal density value obtained by step 3 corresponds to the thorium member of the spectrum curve of depth in conjunction with it The density value of the density curve of the content and corresponding depth of element and potassium element, ρ is established by statistical analysisshWith THOR, POTA, ρb Relational expression, i.e. clay skeletal density computation model
ρsh=a* ρb+b*THOR+c*POTA+M
It is clay skeletal density in formula, THOR is the content of element thorium, POTA is Element Potassium content is that log measurement is close Angle value, a, b, c are coefficient, and M is constant term;
Coefficient a, b, c's seeks according to equation group;It is as follows:
Known parameters ρ in equation groupsh1、ρsh2..., ρshnFor clay skeletal density value, ρb1、ρb2..., ρbnFor clay skeleton The corresponding log measurement density value of density value, (THOR)1, (THOR)2..., (THOR)nIt is worth corresponding thorium for clay skeletal density Content, (POTA)1、(POTA)2..., (POTA)nBe worth corresponding potassium content for clay skeletal density, solve equation group obtain a, b, C, M values.
Further, mineral section clay skeletal density computation model that step 4 obtains, step 2 obtained with it is mixed Conjunction skeletal density computation model is combined inverting and obtains mixed matrix density value;
ρ in formulamaFor mixed matrix density value.
Further, mixed matrix density value combination density porosity volume-based model step (5) obtained calculates hole Degree:
Φ is to calculate porosity in formula.
Compared with prior art, the present invention has following technique effect:
The present invention is a kind of technology based on the soil skeleton inverting complex lithology reservoir porosity that becomes sticky, can be accurate and effective Lithology inversion mineral constituent it is complicated, clay type is more, the porosity of the big reservoir of buried depth.This technology not only avoids number The disadvantage that statistic law calculates porosity universality difference is managed, also solves the reservoir of mineral constituent complexity, clay mineral type is more, no Big with the aqueous degree difference of clay under buried depth accumulating condition, clay skeleton value chooses difficult problem, while also saving special point Analyse the expense of the experiment of clay matrix parameter and mixed matrix density experiment.The technology effectively raises complex lithology reservoir hole Inversion accuracy is spent, the evaluation demand for storing up work in oilfield stimulation is met.
The technology of the present invention provides technological borrowing for the porosity calculation of complex lithology reservoir.
Description of the drawings
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment Attached drawing is briefly described.It should be evident that drawings in the following description are only some embodiments of the invention, for For those of ordinary skill in the art, without creative efforts, other are can also be obtained according to these attached drawings Attached drawing.
Fig. 1 is that determination provided in an embodiment of the present invention becomes sticky the method flow diagram of soil skeleton porosity;
Fig. 2 is that skeletal density parameter list is commonly used in present example;
Fig. 3 is XX wells X diffraction total rock mineral scale lithology scanning well logging result map in present example;
Fig. 4 is the XX well porosity calculations well logging result map obtained using the method for the present invention in present example.
Specific implementation mode
Below in conjunction with attached drawing, the present invention is further described:
- Fig. 4 is please referred to Fig.1, a kind of porosity calculation method for the soil skeleton that becomes sticky, this method comprises the following steps:
Step 1:To rock core carry out physical property experiment, X diffraction total rock mineral experiment, obtain core porosity, mineral constituent, Mineral content data;
The serial rock core of complex lithology reservoir feature can be represented by choosing, according to《Core analysis method SY/T5336- 2006》Flow as defined in standard carries out physical property experiment, the experiment of X diffraction total rock mineral to rock core, and the rock core physical property experiment includes Core porosity is tested and core permeability experiment, and obtains core porosity and core experiment permeability value, and the X diffraction is complete The experiment of rock mineral includes rock core mineral constituent and mineral content experiment, and obtains rock core mineral constituent and mineral content.
Step 2:Well logging is scanned with the X diffraction total rock mineral experimental data scale lithology that step 1 obtains, accurately calculates mine Object section.
By X diffraction total rock mineral experimental analyses, XX well essential mineral groups are divided into Yun Yan, limestone, sandstone, clay, pyrite, Anhydrite is scanned using X diffraction total rock mineral experimental data scale lithology and is logged well, obtains cuing open with the continuous mineral of change in depth Face.The roads Tu3 2-7 be XX well X diffraction total rock mineral data scale lithology scanning well logging curve, the 8th be it is continuous Mineral section.
Step 3:Core porosity, mineral section and density curve, the mixed matrix obtained using step 1 and step 2 Density computation model, density porosity volume-based model are combined inverting and obtain clay skeletal density value.It specifically will be according to following public affairs Formula carries out:
Wherein, ρmaFor mixed matrix density, ρshIndicate clay skeletal density, vshIndicate clay mineral content, ρmaiIt indicates The density skeleton (not including clay mineral) of i-th kind of mineral, vmaiIndicate i-th kind of mineral content (not including clay mineral), Φ tables Show core porosity, ρbIndicate log measurement density value, ρfIndicate the density of water flooding, n > 1.
Known parameters ρ in formulamai、ρfIt is obtained by common skeletal density parameter list, vmai、vshWell logging determination is scanned by lithology Mineral section obtains, and Φ is obtained by rock core physical property experiment data, ρbIt is obtained by density curve, solves equation group and obtain ρshValue.
Fig. 4 is XX well porosity calculations well logging result map, and the 5th is the clay skeletal density value of inverting.
Step 4:Clay skeletal density value and spectrum curve, density curve combine and obtain clay skeletal density computation model.
The soil skeletal density value obtained by step 3 is discrete data, and the clay skeleton that cannot characterize each sampled point is close Angle value, a curve for continuously reacting each depth point clay skeletal density value, needs to establish these core samples in order to obtain The relational expression of the rock core clay skeletal density value of product and sensitive log, continuous clay skeleton can be sought by relational expression Density value is a continuous curve in logging diagram.
It can get the content of uranium in stratum, three kinds of thorium, potassium elements, clay skeleton and thorium and potassium element by spectrometry logging It is in close relations, while clay skeleton value is influenced by density curve, therefore chooses thorium element, potassium element and log measurement density Value characterizes clay skeletal density value.
The clay skeletal density value obtained using step 3 corresponds to the thorium element of the spectrum curve of depth and potassium member in conjunction with it The density value of the density curve of the content of element and corresponding depth, ρ is established by statistical analysisshWith THOR, POTA, ρbRelationship Formula, i.e. clay skeletal density computation model.It will specifically be carried out according to following invention formula:
ρsh=a* ρb+b*THOR+c*POTA+M
ρ in formulashFor clay skeletal density, THOR is the content of element thorium, POTA is the content of Element Potassium, ρbLog measurement Density value, a, b, c are coefficient, and M is constant term.
Coefficient a, b, c's seeks according to equation group.It is as follows:
Known parameters ρ in equation groupsh1、ρsh2..., ρshnFor clay skeletal density value, ρb1、ρb2..., ρbnFor clay skeleton The corresponding log measurement density value of density value, (THOR)1, (THOR)2..., (THOR)nIt is worth corresponding thorium for clay skeletal density Content, (POTA)1、(POTA)2..., (POTA)nBe worth corresponding potassium content for clay skeletal density, solve equation group obtain a, b, C, M values.
Obtained clay skeletal density value inverse model is as follows:
ρsh=3.34* ρb+0.004*THOR+0.029*POTA-6.43
Fig. 4 seconds are conventional logging density curve, and third road is thorium and potassium curve in spectrum curve, and the 6th is meter The clay skeletal density curve of calculation.
Step 5:By mineral section that clay skeletal density computation model that step 4 obtains, step 2 obtain with mix Skeletal density computation model is combined inverting and obtains mixed matrix density value.Specifically carried out according to following formula:
ρ in formulamaFor mixed matrix density value.
Fig. 4 the 7th is mixed matrix density curve.
Step 6:The mixed matrix density value combination density porosity volume-based model inverting that step 5 obtains is obtained into hole Degree.Specifically carried out according to following formula:
The 8th A-POR of Fig. 4 are experimental analysis core porosity, are indicated with rod-shaped in figure;POR is to utilize foregoing invention The porosity curve of calculating.As can be seen that porosity and experiment core analysis porosity one that the inventive technique calculates from figure Cause property is preferable, and it is high to calculate porosity precision.
The present invention is elaborated above, specific case is used in combination to carry out the principle of the present invention and embodiment Detailed introduction, the explanation of above example are intended merely to help to understand embodiments of the present invention and core concept, the present invention Implementation be not only limited to above-described embodiment, other are any based on changing made by the method for the present invention principle, replace, group It closes, is equivalent substitute mode, within the scope of the present invention.

Claims (6)

1. a kind of porosity calculation method for the soil skeleton that becomes sticky, which is characterized in that include the following steps:
Step 1:Physical property experiment, the experiment of X diffraction total rock mineral are carried out to rock core, obtain core porosity, mineral constituent, mineral Content data;
Step 2:Well logging is scanned with the X diffraction total rock mineral experimental data scale lithology that step 1 obtains, the accurate mineral that calculate cut open Face;
Step 3:The core porosity, mineral section and the density curve that are obtained according to step 1 and step 2, mixed matrix density Computation model, density porosity volume-based model are combined inverting and obtain clay skeletal density value;
Step 4:The clay skeletal density value that step 3 obtains is obtained into clay bone with spectrum curve, density curve combination inverting Frame density computation model;
Step 5:The mineral section and mixed matrix that clay skeletal density computation model that step 4 obtains, step 2 are obtained Density computation model is combined inverting and obtains mixed matrix density value;
Step 6:The mixed matrix density value combination density porosity volume-based model inverting that step 5 obtains is obtained into porosity.
2. a kind of porosity calculation method of soil skeleton that becomes sticky according to claim 1, which is characterized in that obtained with step 1 The X diffraction total rock mineral experimental data scale lithology scanning well logging arrived, obtains the continuous mineral section with change in depth.
3. a kind of porosity calculation method of soil skeleton that becomes sticky according to claim 1, which is characterized in that according to step 1 Core porosity, mineral section and density curve, mixed matrix density computation model, the density porosity body obtained with step 2 Product module type is combined inverting and obtains clay skeletal density value;It will specifically be carried out according to following formula:
Wherein, ρmaFor mixed matrix density, ρshIndicate clay skeletal density, vshIndicate clay mineral content, ρmaiIndicate i-th kind The density skeleton of mineral and include clay mineral, vmaiIt indicates i-th kind of mineral content and does not include clay mineral, Φ indicates rock Heart porosity, ρbIndicate log measurement density value, ρfIndicate the density of water flooding, n > 1;
Known parameters ρ in formulamai、ρfIt is obtained by common skeletal density parameter list, vmai、vshThe determining mineral of well logging are scanned by lithology Section obtains, and Φ is obtained by rock core physical property experiment data, ρbIt is obtained by density curve, solves equation group and obtain ρshValue.
4. a kind of porosity calculation method of soil skeleton that becomes sticky according to claim 1, which is characterized in that obtained by step 3 The clay skeletal density value arrived corresponds to the thorium element of spectrum curve of depth and the content of potassium element in conjunction with it and corresponds to depth The density value of density curve establishes ρ by statistical analysisshWith THOR, POTA, ρbRelational expression, i.e., clay skeletal density calculate Model
ρsh=a* ρb+b*THOR+c*POTA+M
It is clay skeletal density in formula, THOR is the content of element thorium, POTA is Element Potassium content is log measurement density Value, a, b, c are coefficient, and M is constant term;
Coefficient a, b, c's seeks according to equation group;It is as follows:
Known parameters ρ in equation groupsh1、ρsh2..., ρshnFor clay skeletal density value, ρb1、ρb2..., ρbnFor clay skeletal density It is worth corresponding log measurement density value, (THOR)1, (THOR)2..., (THOR)nIt is worth corresponding thorium content for clay skeletal density, (POTA)1、(POTA)2..., (POTA)nIt is worth corresponding potassium content for clay skeletal density, solves equation group and obtain a, b, c, M Value.
5. a kind of porosity calculation method of soil skeleton that becomes sticky according to claim 1, which is characterized in that obtain step 4 To clay skeletal density computation model, the obtained mineral section of step 2 be combined inverting with mixed matrix density computation model Obtain mixed matrix density value;
ρ in formulamaFor mixed matrix density value.
6. a kind of porosity calculation method of soil skeleton that becomes sticky according to claim 1, which is characterized in that by step (5) Obtained mixed matrix density value combination density porosity volume-based model calculates porosity:
Φ is to calculate porosity in formula.
CN201810260747.9A 2018-03-27 2018-03-27 Porosity calculation method for clay-changing framework Active CN108548765B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810260747.9A CN108548765B (en) 2018-03-27 2018-03-27 Porosity calculation method for clay-changing framework

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810260747.9A CN108548765B (en) 2018-03-27 2018-03-27 Porosity calculation method for clay-changing framework

Publications (2)

Publication Number Publication Date
CN108548765A true CN108548765A (en) 2018-09-18
CN108548765B CN108548765B (en) 2020-11-06

Family

ID=63517417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810260747.9A Active CN108548765B (en) 2018-03-27 2018-03-27 Porosity calculation method for clay-changing framework

Country Status (1)

Country Link
CN (1) CN108548765B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113720745A (en) * 2020-05-26 2021-11-30 中国石油化工股份有限公司 Method for calculating porosity of reservoir stratum containing carbon debris by geophysical logging

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101832133A (en) * 2010-01-28 2010-09-15 中国石油集团川庆钻探工程有限公司 Reservoir fluid type discrimination method based on difference value of density porosity and neutron porosity
CN104564037A (en) * 2013-10-10 2015-04-29 中国石油天然气股份有限公司 Shale gas reservoir brittle mineral content logging calculation method
CN105205296A (en) * 2014-06-09 2015-12-30 中国石油化工股份有限公司 Method for calculating porosity of shale gas reservoir
CN107290782A (en) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method
CN107829731A (en) * 2017-11-06 2018-03-23 陈国军 A kind of volcanic rock porosity correction method of clay alteration

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101832133A (en) * 2010-01-28 2010-09-15 中国石油集团川庆钻探工程有限公司 Reservoir fluid type discrimination method based on difference value of density porosity and neutron porosity
CN104564037A (en) * 2013-10-10 2015-04-29 中国石油天然气股份有限公司 Shale gas reservoir brittle mineral content logging calculation method
CN105205296A (en) * 2014-06-09 2015-12-30 中国石油化工股份有限公司 Method for calculating porosity of shale gas reservoir
CN107290782A (en) * 2016-03-30 2017-10-24 中国石油化工股份有限公司 Reservoir porosity, water saturation and shale content parameter Simultaneous Inversion new method
CN107829731A (en) * 2017-11-06 2018-03-23 陈国军 A kind of volcanic rock porosity correction method of clay alteration

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113720745A (en) * 2020-05-26 2021-11-30 中国石油化工股份有限公司 Method for calculating porosity of reservoir stratum containing carbon debris by geophysical logging

Also Published As

Publication number Publication date
CN108548765B (en) 2020-11-06

Similar Documents

Publication Publication Date Title
CN104516017B (en) A kind of carbonate rock petrophysical parameter seismic inversion method
CN106950347B (en) A method of evaluation mud shale each group partial volume
US9260948B2 (en) Multi-level reservoir history matching
CN101414013B (en) Method for determining underground fluid by using seismic data
CN106154351A (en) A kind of evaluation method of low porosity permeability reservoir permeability
Thiele et al. Polymer-flood modeling using streamlines
US10359529B2 (en) Singularity spectrum analysis of microseismic data
Gherabati et al. Assessment of hydrocarbon in place and recovery factors in the Eagle Ford Shale play
CN104806232B (en) A kind of method for determining porosity lower limit of fracture
CN105822302A (en) Oil-water distribution recognition method based on well-ground potential method
CN104863574B (en) A kind of Fluid Identification Method suitable for tight sandstone reservoir
CN108301825A (en) A kind of high temperature and pressure reservoir pore pressure analysis method
CN104252007A (en) Compatibility rock physics modeling method
Alpak et al. Simplified modeling of turbidite channel reservoirs
CN105840189A (en) Shale gas rock physics and logging method
CN107728205B (en) A kind of Formation pressure prediction method
CN107748399A (en) Mountain front infrastructure layer method is identified using Interface Inversion of Gravity Field
CN109425900A (en) A kind of Seismic Reservoir Prediction method
CN109509111A (en) The prediction technique and system of prospect pit strata pressure
CN107795320A (en) A kind of computational methods of horizontal well carbonate reservoir parameter
CN103946896A (en) System and method for assessing heterogeneity of a geologic volume of interest with process-based models and dynamic heterogeneity
CN106204302A (en) Method for calculating original water saturation and application
Glegola et al. Gravimetric monitoring of water influx into a gas reservoir: A numerical study based on the ensemble kalman filter
Prather et al. Hierarchy of deep-water architectural elements with reference to seismic resolution: implications for reservoir prediction and modeling
CN108548765A (en) A kind of porosity calculation method for the soil skeleton that becomes sticky

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