CN106093350A - The method determining heterogeneous carbonate reservoir saturation exponent - Google Patents

The method determining heterogeneous carbonate reservoir saturation exponent Download PDF

Info

Publication number
CN106093350A
CN106093350A CN201610649337.4A CN201610649337A CN106093350A CN 106093350 A CN106093350 A CN 106093350A CN 201610649337 A CN201610649337 A CN 201610649337A CN 106093350 A CN106093350 A CN 106093350A
Authority
CN
China
Prior art keywords
rock core
rock
saturation
reservoir
core sample
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
CN201610649337.4A
Other languages
Chinese (zh)
Other versions
CN106093350B (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
Original Assignee
China Petroleum and Natural Gas 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 filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201610649337.4A priority Critical patent/CN106093350B/en
Publication of CN106093350A publication Critical patent/CN106093350A/en
Priority to PCT/CN2017/084304 priority patent/WO2018028258A1/en
Priority to US16/095,432 priority patent/US11112527B2/en
Application granted granted Critical
Publication of CN106093350B publication Critical patent/CN106093350B/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
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK 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
    • G01V20/00Geomodelling in general
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Geophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Mining & Mineral Resources (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Fluid Mechanics (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

A kind of method determining heterogeneous carbonate reservoir saturation exponent of disclosure, including: according to pre-defined rule, target rock reservoir is divided at least two Reservoir type;Obtain the corresponding relation between saturation exponent and irreducible water saturation in every kind of Reservoir type;Reservoir type belonging to rock core to be measured is determined according to described pre-defined rule;Obtain the irreducible water saturation of described rock core to be measured;Irreducible water saturation according to described rock core to be measured described corresponding relation based on described affiliated Reservoir type calculates the saturation exponent of described rock core to be measured.

Description

The method determining heterogeneous carbonate reservoir saturation exponent
Technical field
The present invention relates to oil-gas reservoir Logging Evaluation of Fractured Reservoir technical field, determine that heterogeneous carbonate stores up particularly to one The method of layer saturation exponent.
Background technology
The electric conductivity of reservoir rock depends primarily on the fluid properties in interstitial space, saturation and spatial distribution thereof, For heterogeneous carbonate reservoir, pore structure affects highly significant to resistivity, sometimes the seam hole shadow to resistivity Ringing the impact of remote super oil-gas possibility, Many researchers attempts to break away from simple dependent resistor rate curve over more than 70 year and calculates saturation Traditional method, attempts utilizing various non-electrical log to calculate carbonate rock saturation, but effect is unsatisfactory.Therefore, it is so far Only, the saturation computation based on electrical log is still that the most practical, feasible method.
Archie (1942) proposes oil-bearing reservoir resistivity Magnification water saturation, formation factor hole the earliest Relational expression between porosity.Li Ning (1989) is based on Inhomogeneous Anisotropic stratigraphic model, by complete mathematical derivation, Give the relational expression between resistivity Magnification water saturation, formation factor porosity.Universal relation formula is from theory On solve the accurate quantification computational problem of heterogeneous body complicated reservoirs saturation, but under prior art conditions, be applied Also needing to solve two key technical problems in actual production: first, universal relation formula is a general solution equation, is not easy to straight Connect for program calculation, i.e. under the conditions of existing logging technique, determine that all parameters in universal relation are extremely difficult, therefore real Border application needs to select to meet according to actual reservoir characteristic the short-form (also referred to as optimised form) of required precision;Secondly, Need to select a kind of method reliably accurately to determine each undetermined parameter in optimised form so that the computation model energy finally determined True rule between reflection reservoir resistivity and water saturation the most to greatest extent.
During hydrocarbon saturation with Archie formula as core calculates, relate to cementation factor m, saturation exponent this Two important litho-electric parameters.The most accurately determine saturation exponent, be petrophysicist and log analyst all the time The emphasis of research.At present, determine that the method for saturation exponent substantially can be divided into two big classes: one is to utilize rock-electric test, logical Relation between over-fitting resistance Magnification water saturation determines;Two is to utilize the saturation exponent and hole set up Empirical relation between the reservoir parameters such as degree, permeability determines.First kind method belongs to direct method, is current saturation exponent The conventional method determined, but there are two shortcomings in actual application: first, in order to determine saturation exponent, need every block of rock The heart carries out change saturation rock-electric test (such as gas drive method), and the time that this experiment spends is the longest;Secondly, such method is utilized It is difficult to be dynamically determined according to reservoir characteristics the numerical value of saturation exponent by well-log information.Equations of The Second Kind method belongs to indirect Method, is the emphasis of research the most both at home and abroad, utilizes such method easily to realize being dynamically determined of saturation exponent.By research, Many researchers it is also proposed and gives high porosity sandstone in SHENGLI PETROLEUM AREA according to reservoir physical parameter, such as Zhu Jiajun (2010) Reservoir saturation exponent and porosity and the relation of formation water salinity, Fu Aibing (2007) passes through piecewise regression and dependency Analyze and research the method for saturation exponent.These researchs have focused largely on sandstone reservoir, carbonate rock saturation exponent It is dynamically determined technique study relatively fewer.Additionally, existing research is more based on saturation exponent and reservoir porosity, infiltration Rate relation analysis, for heterogeneous body carbonate complex rock, research shows the change of saturation exponent under same holes porosity, permeability Change scope is relatively big, the most accurately determines that the numerical value of carbonate rock saturation exponent is faced in current logging evaluation Big challenge.
Summary of the invention
For above-mentioned technical problem, the present invention proposes a kind of side determining heterogeneous carbonate reservoir saturation exponent Method, can accurately determine heterogeneous carbonate saturation exponent, thus it is dynamic to realize carbonate rock saturation model parameter Determine, improve the computational accuracy of hydrocarbon saturation.
For reaching above-mentioned purpose, the application provides a kind of method determining heterogeneous carbonate reservoir saturation exponent, Including:
According to pre-defined rule, target rock reservoir is divided at least two Reservoir type;
Obtain the corresponding relation between saturation exponent and irreducible water saturation in every kind of Reservoir type;
Reservoir type belonging to rock core to be measured is determined according to described pre-defined rule;
Obtain the irreducible water saturation of described rock core to be measured;
Irreducible water saturation according to described rock core to be measured described corresponding relation based on described affiliated Reservoir type calculates The saturation exponent of described rock core to be measured.
As one preferred embodiment, described according to pre-defined rule, target rock reservoir is divided at least two reservoir class Type includes:
Choose the rock core sample of polylith target rock reservoir;
Originally it is divided at least two class rock cores to represent at least two Reservoir type core sample described in polylith according to pre-defined rule.
As one preferred embodiment, the rock core sample of more than 10 pieces target rock reservoirs is chosen.
As one preferred embodiment, the rock core sample choosing polylith target rock reservoir described in includes:
Determine the position of target rock reservoir, effective thickness;
The position of coring of target rock reservoir is determined according to Image Logging Data;
Multiple rock core sample is drilled through in described position of coring.
As one preferred embodiment, described according to pre-defined rule, core sample described in polylith is originally divided at least two classes Rock core includes to represent at least two Reservoir type:
According to pre-defined rule, core sample described in polylith is originally divided into two class rock cores;Described two class rock cores include: a class is hole Hole is grown but the rock core sample of hole poor connectivity;Another kind of is containing crack or the rock core sample based on intracrystalline pore.
As one preferred embodiment, described pre-defined rule includes:
If it can be seen that Porous Characteristic but free from flaw feature in the Image Logging Data that rock core is corresponding, or oozing intersection in hole In figure, the data point of this rock core is positioned at lower right, or during core observation naked eyes it can be seen that corrosion hole but deposit without microcrack , then this rock core can be classified as hole grow but the rock core of hole poor connectivity, otherwise rock core can range containing crack or with Intracrystalline pore is main rock core.
As one preferred embodiment, saturation exponent and irreducible water saturation in every kind of Reservoir type of described acquisition Between corresponding relation include:
Obtain the irreducible water saturation of every piece of described rock core sample;
Obtain the saturation exponent of every piece of described rock core sample;
By the saturation exponent of rock core sample described in every kind of Reservoir type with irreducible water saturation linear fit to obtain Saturation exponent and the corresponding relation of irreducible water saturation in every kind of Reservoir type;Described corresponding relation is as follows:
N=aSwir+b;
N is saturation exponent;Swir is irreducible water saturation.
As one preferred embodiment, saturation exponent and irreducible water saturation in every kind of Reservoir type of described acquisition Between corresponding relation also include:
The described rock core sample of partial amt in every kind of Reservoir type is chosen according to described irreducible water saturation;
Corresponding, the saturation exponent of described acquisition every piece described rock core sample includes:
Obtain in every kind of Reservoir type the saturation of every piece of described rock core sample in the described rock core sample of partial amt to refer to Number.
As one preferred embodiment, the irreducible water saturation of described acquisition every block described rock core sample includes:
Rock core sample every piece described is carried out hole and oozes parameter measurement and centrefuge experiment to determine irreducible water saturation.
As one preferred embodiment, the saturation exponent of described acquisition every block described rock core sample includes:
The described rock core sample of partial amt in every kind of Reservoir type is carried out rock-electric test to obtain saturation exponent.
As one preferred embodiment, the irreducible water saturation of the described rock core to be measured of described acquisition includes:
Described rock core to be measured is carried out hole and oozes parameter measurement and centrefuge experiment to determine irreducible water saturation.
By above description it can be seen that compared with the determination method of tradition reservoir saturation exponent, the present invention proposes Method has following significant advantage: (1) determines that method is easy: once this area's saturation n and irreducible water saturation Swir Between relation set up, utilize the method, be not required to again interval of interest rock core be carried out reservoir conditions rock-electric test, can by bundle Tie up water saturation and determine the numerical value of other layer of reservoir saturation exponent in this area, thus it is multiple to overcome reservoir conditions rock-electric test Miscellaneous difficulty;(2) parameter is the most accurate: utilize the method can determine the saturated of correspondence according to target zone position different reservoir feature Degree index, embodies the different pore structures difference on electrically impact;(3) it is easy to utilize well-log information to be dynamically determined: utilizing should Method, utilizes well-log information can conveniently realize the dynamic calculation of saturation exponent, thus improves carbonate complex rock oil The computational accuracy of gas saturation.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments of invention, for those skilled in the art, on the premise of not paying creative work, it is also possible to root Other accompanying drawing is obtained according to these accompanying drawings.
Fig. 1 is the method for the determination heterogeneous carbonate reservoir saturation exponent that a kind of embodiment of the application is provided Flow chart of steps;
Fig. 2 is that Changqing oilfields hole is grown but the reservoir saturation exponent of hole poor connectivity and irreducible water saturation Corresponding relation figure between Swir;
Fig. 3 is that Changqing oilfields contains crack or the reservoir saturation exponent based on intracrystalline pore and irreducible water saturation Swir Between corresponding relation figure.
Detailed description of the invention
For the technical scheme making those skilled in the art be more fully understood that in the application, real below in conjunction with the application Execute the accompanying drawing in example, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described enforcement Example is only some embodiments of the present application rather than whole embodiments.Based on the embodiment in the application, this area is common The every other embodiment that technical staff is obtained on the premise of not making creative work, all should belong to the present invention and protect The scope protected.
As it is shown in figure 1, the application one embodiment provides one to determine heterogeneous carbonate reservoir saturation exponent Method, the method comprises the following steps:
S100: target rock reservoir is divided at least two Reservoir type according to pre-defined rule;
This step is to consider the pore structure characteristic of heterogeneous carbonate reservoir, and target rock reservoir is divided at least two Plant Reservoir type.There is difference in the pore structure of different reservoir type.
In one embodiment, described according to pre-defined rule, target rock reservoir is divided at least two Reservoir type (step Rapid S100) include following sub-step:
S110, choose the rock core sample of polylith target rock reservoir;
Reflection Reservoir levels preferably, in this step, preferably chooses the core sample of more than 10 pieces target rock reservoirs This.
This step S110 can be chosen survey region representational rock core sample, refer to conventional, imaging logging On Information integration analysis foundation, determine the position of target zone position, effective thickness, and according to data such as routine, imaging loggings Variation characteristic in the vertical, determines and cores position and drill through representational rock core sample.
Specifically, choose the rock core sample (step S110) of polylith target rock reservoir described in and can include following sub-step: S111, determine the position of target rock reservoir, effective thickness;S112, determine coring of target rock reservoir according to Image Logging Data Position;S113, drill through multiple rock core sample in described position of coring.
S120, originally it is divided at least two class rock cores to represent at least two reservoir core sample described in polylith according to pre-defined rule Type;
In this step, represented the Porous Characteristic of its place reservoir by the Porous Characteristic of rock core sample, on this is also The purpose of one many more than the 10 pieces rock core samples of step S110 choosing, represented place rock storage by core sample exactly originally with finer The feature of layer.
It is also preferred that the left described, core sample described in polylith is originally divided at least two class rock cores to represent at least two according to pre-defined rule Kind of Reservoir type (step S120) may include that S121, core sample described in polylith be originally divided into according to pre-defined rule two class rock cores; Described two class rock cores include: a class is that hole is grown but the rock core sample (first kind) of hole poor connectivity;Another kind of is containing splitting Seam or the rock core sample (Equations of The Second Kind) based on intracrystalline pore.
Wherein, if described pre-defined rule may include that in the Image Logging Data that rock core is corresponding it can be seen that Porous Characteristic But free from flaw feature, or in cross plot (porosity-permeability cross plot) is oozed in hole, the data point of this rock core is positioned at lower right, Or naked eyes during core observation it can be seen that corrosion hole but without microcrack exist, then this rock core can be classified as hole grow but hole The rock core of hole poor connectivity, otherwise rock core can range containing crack or the rock core based on intracrystalline pore.
S200: obtain the corresponding relation between saturation exponent and irreducible water saturation in every kind of Reservoir type;
Wherein, this step S200 can be performed by following steps: S210, obtains the irreducible water of every piece of described rock core sample Saturation;S220, obtain the saturation exponent of every piece of described rock core sample;S230, by core sample described in every kind of Reservoir type This saturation exponent and irreducible water saturation linear fit are satisfied with irreducible water to obtain saturation exponent in every kind of Reservoir type Corresponding relation with degree.Described corresponding relation is: n=aSwir+b;Wherein, n is saturation exponent;Swir is that irreducible water is saturated Degree.
In step s 200, first to selected rock core sample, carry out porosity, Permeability Parameters measurement, the most right Selected rock core is centrifuged experiment.In experiment, the selection of centrifugal force should be with reference to the hole of reservoir and pressure characteristic.Centrifugal real In testing, should first measure the full wet weight of rock core sample, under certain rotating speed (centrifugal force in other words), then be centrifuged experiment, Experiment measures the rock core sample weight after being centrifuged after terminating, and calculates the irreducible water saturation of rock core.That is, described every piece of institute of acquisition The irreducible water saturation (step S210) stating rock core sample may include that rock core sample every piece described is carried out hole oozes parameter measurement And centrefuge experiment is to determine irreducible water saturation.
In experiment, type and the salinity of formation water determine according to survey region formation water analysis of data, experimental temperature and Confined pressure determines according to reservoir temperature and the pressure of interval of interest.On the basis of (displacement) rock-electric test, draw the electricity of rock core sample Graph of a relation between resistance Magnification and water saturation, and utilize Archie formula (Archie formula) to determine that every block of rock core is saturated The numerical value of degree index n.That is, the saturation exponent (S220) of described every piece of described rock core sample of acquisition may include that every kind of storage In layer type, the described rock core sample of partial amt carries out rock-electric test to obtain saturation exponent.
Further, for obtaining rationally result of calculation the most accurately, in every kind of Reservoir type of described acquisition saturation exponent with Corresponding relation (step S200) between irreducible water saturation can also include step: S225, according to described irreducible water saturation Choose the described rock core sample of partial amt in every kind of Reservoir type.
Corresponding, saturation exponent step S230 of described acquisition every piece described rock core sample then include (being): S231, Obtain in every kind of Reservoir type the saturation exponent of every piece of described rock core sample in the described rock core sample of partial amt.
Irrational data can be rejected, i.e. irreducible water saturation does not substantially meet this rock core sample by step S225 Place Reservoir levels, thus improve the precision of calculating.Research shows, first kind rock core irreducible water saturation is the most relatively low, second The irreducible water saturation of class rock core is the highest.In each class rock core, select there is representative according to the numerical value of irreducible water saturation The rock core of property carries out rock-electric test.Such as: if survey region, first kind rock core irreducible water saturation is 50% to the maximum, then need to be First kind rock core select irreducible water saturation carry out rock-electric test at the representative rock core of 0-50%, the number more than 50% According to rejecting.After rejecting irrational rock core sample, every piece of remaining rock core sample is all carried out rock-electric test (i.e. step S231)。
In step S230, according to selected irreducible water saturation Swir of representative rock core, saturation exponent Result of calculation, determined between this area's different reservoir type saturation exponent and irreducible water saturation Swir by matching Relation, as a example by above-mentioned two class Reservoir types:
1st class: hole is grown but the saturation exponent of the reservoir of hole poor connectivity pass corresponding with irreducible water saturation System is:
N=a1Swir+b1 (1)
2nd class: the saturation exponent containing crack or the reservoir based on intracrystalline pore and the corresponding relation of irreducible water saturation For
N=a2Swir+b2 (2)
Parameter a in above-mentioned formula1、b1And a2、b2It is constant for the certain layer position of particular locality, but different regions, The numerical value of different layers position above-mentioned parameter there are differences, parameter a1、b1And a2、b2Numerical value can be determined by core experiment.
S300: determine Reservoir type belonging to rock core to be measured according to described pre-defined rule;
In this step, S300 can confirm Reservoir type belonging to rock core to be measured, rock core to be measured by above-mentioned pre-defined rule Do not carry out rock-electric test.
Wherein, if described pre-defined rule may include that in the Image Logging Data that rock core is corresponding it can be seen that Porous Characteristic But free from flaw feature, or the data point of this rock core is positioned at lower right, or naked eyes energy during core observation in cross plot is oozed in hole Enough see corrosion hole but exist without microcrack, then this rock core can be classified as hole growth but the rock core of hole poor connectivity, no Then rock core can range containing crack or the rock core based on intracrystalline pore.
S400: obtain the irreducible water saturation of described rock core to be measured;
In this step, rock core to be measured carrying out porosity, Permeability Parameters measurement, rock core the most to be measured is centrifuged reality Test.In experiment, the selection of centrifugal force should be with reference to the hole of reservoir and pressure characteristic.In centrefuge experiment, should first measure rock core to be measured Full wet weight, is then centrifuged experiment under certain rotating speed (centrifugal force in other words), and experiment is measured after terminating after being centrifuged Rock core weight to be measured, and calculate the irreducible water saturation of rock core to be measured.
S500: according to the irreducible water saturation of described rock core to be measured described corresponding relation based on described affiliated Reservoir type Calculate the saturation exponent of described rock core to be measured.
Such as, determine that Reservoir type belonging to rock core to be measured is the 1st class by step S300, determine that it is restrainted by step S400 Tiing up water saturation is SwirTreat, it is Swir in step S500TreatThat brings in step S230 in the 1st class Reservoir type is saturated Degree index and the corresponding relation of irreducible water saturation, thus obtain the saturation exponent of rock core to be measuredTreat=a1SwirTreat+b1
By above description it can be seen that compared with the determination method of tradition reservoir saturation exponent, the present invention proposes Method has following significant advantage: (1) determines that method is easy: once this area's saturation n and irreducible water saturation Swir Between relation set up, utilize the method, be not required to again interval of interest rock core be carried out reservoir conditions rock-electric test, can by bundle Tie up water saturation and determine the numerical value of other layer of reservoir saturation exponent in this area, thus it is multiple to overcome reservoir conditions rock-electric test Miscellaneous difficulty;(2) parameter is the most accurate: utilize the method can determine the saturated of correspondence according to target zone position different reservoir feature Degree index, embodies the different pore structures difference on electrically impact;(3) it is easy to utilize well-log information to be dynamically determined: utilizing should Method, utilizes well-log information can conveniently realize the dynamic calculation of saturation exponent, thus improves carbonate complex rock oil The computational accuracy of gas saturation.
Below in conjunction with the accompanying drawings, the detailed description of the invention of the present invention is described in further detail.
On to Information integration analysis foundations such as Changqing oilfields interval routine, imaging loggings, it is determined that target zone position Position, effective thickness, and according to data variation characteristics in the vertical such as routine, imaging loggings, determine the position boring of coring Take 19 pieces of representational plunger rock core (this is rock core sample).
First to 19 blocks of selected plunger rock cores, porosity, Permeability Parameters measurement are carried out, then to selected rock The heart is centrifuged experiment.In centrefuge experiment, first measure every block of rock core and satisfy wet weight, then under rotating speed is 3000 revolutions per seconds Being centrifuged experiment, the weight after every block of rock core is centrifuged is measured in experiment after terminating.
Image Logging Data, properties of pore and permeability and core observation according to the selected rock core correspondence degree of depth is by above-mentioned 19 blocks of rock cores Being divided into two classes that pore structure is different: a class is that hole is grown but the rock core of hole poor connectivity, totally 10 pieces, another kind of is containing splitting Seam or the rock core based on intracrystalline pore, totally 9 pieces.Carry out classifying according to rock core pore structure method particularly includes: if rock core is corresponding Image Logging Data on it can be seen that Porous Characteristic but free from flaw feature, or the data point of this rock core in cross plot is oozed in hole Be positioned at lower right, or during core observation naked eyes it can be seen that corrosion hole but exist without microcrack, then this rock core can be classified as Hole is grown but the rock core of hole poor connectivity, and otherwise rock core can range Equations of The Second Kind, i.e. containing crack or based on intracrystalline pore Rock core.If rock core has carried out three dimensional CT test, then according to CT data, rock core is divided into hole and grows but hole poor connectivity Rock core, containing crack or the rock core based on intracrystalline pore, the easiest, accurately.
Rock core to above-mentioned two class different pore structures, selects representative rock core to carry out according to irreducible water saturation Reservoir conditions displacement rock-electric test, and calculate the saturation exponent of rock core.According to irreducible water saturation, select Changqing oilfields The layer position first kind (i.e. hole is grown but hole poor connectivity) 7 pieces of rock core, Equations of The Second Kind (i.e. containing crack or based on intracrystalline pore) rock 5 pieces of the heart has carried out reservoir conditions semi-permeable plate gas drive rock-electric test.In experiment, the salinity of saturated brine be 100000ppm, Water type is NaCl type, and formation water resistivity Rw is 0.032 ohm meter.According to purpose region depth of reservoirs, in experiment, confined pressure is 15MPa。
According to the resistivity of rock core under the different water cut saturation obtained in experiment, depict every piece of rock core resistance Magnification And the graph of a relation between water saturation, and Archie formula is utilized to determine the numerical value of every piece of rock core saturation exponent.
According to selected irreducible water saturation Swir of representative rock core, the result of calculation of saturation exponent, logical Over-fitting determines the relation between this area's different reservoir type saturation exponent and irreducible water saturation Swir:
Changqing oilfields layer position first kind reservoir, i.e. hole are grown but the reservoir irreducible water saturation of hole poor connectivity Graph of a relation between Swir and saturation exponent is shown in accompanying drawing 2, and quantitative fit correlation formula is:
N=2.342Swir+0.857 (1)
Changqing oilfields layer position Equations of The Second Kind reservoir, i.e. containing crack or reservoir irreducible water saturation Swir based on intracrystalline pore And the graph of a relation between saturation exponent is shown in accompanying drawing 3, quantitative fit correlation formula is:
N=-1.957Swir+2.919 (2)
Changqing oilfields interval is not carried out store up rock-electric test rock core A, B (this is rock core to be measured), first analyze this two The pore structure type of block rock core, according to the above analysis method (pre-defined rule), determine that rock core A is the first kind, and rock core B is second Class, irreducible water saturation Swir calculating this rock core A according to centrifuge results is 0.38, and irreducible water saturation Swir of rock core B is 0.71, utilize the relation between aforementioned two class reservoir core irreducible water saturations Swir and the saturation exponent having built up, can With the quick saturation exponent of rock core A that is calculated preparatively for 1.75, the saturation exponent of rock core B is 1.53.Similarly, may be used In order to determining other rock core and the saturation exponent of corresponding interval reservoir in aforementioned manners.
The ultimate principle of the present invention, principal character and advantages of the present invention have more than been shown and described.The technology of the industry Personnel, it should be appreciated that the present invention is not restricted to the described embodiments, simply illustrating this described in above-described embodiment and description The principle of invention, without departing from the spirit and scope of the present invention, the present invention also has various changes and modifications, and these become Change and improvement both falls within the scope of protection of present invention.Claimed scope by appending claims and Equivalent defines.

Claims (11)

1. the method determining heterogeneous carbonate reservoir saturation exponent, it is characterised in that including:
According to pre-defined rule, target rock reservoir is divided at least two Reservoir type;
Obtain the corresponding relation between saturation exponent and irreducible water saturation in every kind of Reservoir type;
Reservoir type belonging to rock core to be measured is determined according to described pre-defined rule;
Obtain the irreducible water saturation of described rock core to be measured;
Irreducible water saturation according to described rock core to be measured described corresponding relation based on described affiliated Reservoir type calculates described The saturation exponent of rock core to be measured.
2. the method for claim 1, it is characterised in that described according to pre-defined rule, target rock reservoir is divided at least two Plant Reservoir type to include:
Choose the rock core sample of polylith target rock reservoir;
Originally it is divided at least two class rock cores to represent at least two Reservoir type core sample described in polylith according to pre-defined rule.
3. method as claimed in claim 2, it is characterised in that choose the rock core sample of more than 10 pieces target rock reservoirs.
4. method as claimed in claim 2, it is characterised in that described in choose the rock core sample of polylith target rock reservoir and include:
Determine the position of target rock reservoir, effective thickness;
The position of coring of target rock reservoir is determined according to Image Logging Data;
Multiple rock core sample is drilled through in described position of coring.
5. method as claimed in claim 2, it is characterised in that described according to pre-defined rule, core sample described in polylith is originally divided into At least two class rock cores include to represent at least two Reservoir type:
According to pre-defined rule, core sample described in polylith is originally divided into two class rock cores;Described two class rock cores include: a class is that hole is sent out Educate but the rock core sample of hole poor connectivity;Another kind of is containing crack or the rock core sample based on intracrystalline pore.
6. method as claimed in claim 5, it is characterised in that described pre-defined rule includes:
If it can be seen that Porous Characteristic but free from flaw feature in the Image Logging Data that rock core is corresponding, or in cross plot is oozed in hole The data point of this rock core is positioned at lower right, or during core observation naked eyes it can be seen that corrosion hole but exist without microcrack, then This rock core can be classified as hole and grow but the rock core of hole poor connectivity, and otherwise rock core can range containing crack or with intracrystalline pore It it is main rock core.
7. the method for claim 1, it is characterised in that saturation exponent and constraint in every kind of Reservoir type of described acquisition Corresponding relation between water saturation includes:
Obtain the irreducible water saturation of every piece of described rock core sample;
Obtain the saturation exponent of every piece of described rock core sample;
By the saturation exponent of rock core sample described in every kind of Reservoir type with irreducible water saturation linear fit to obtain every kind Saturation exponent and the corresponding relation of irreducible water saturation in Reservoir type;Described corresponding relation is as follows:
N=aSwir+b;
N is saturation exponent;Swir is irreducible water saturation.
8. method as claimed in claim 7, it is characterised in that saturation exponent and constraint in every kind of Reservoir type of described acquisition Corresponding relation between water saturation also includes:
The described rock core sample of partial amt in every kind of Reservoir type is chosen according to described irreducible water saturation;
Corresponding, the saturation exponent of described acquisition every piece described rock core sample includes:
Obtain in every kind of Reservoir type the saturation exponent of every piece of described rock core sample in the described rock core sample of partial amt.
9. method as claimed in claim 7, it is characterised in that the irreducible water saturation of described acquisition every piece described rock core sample Including:
Rock core sample every piece described is carried out hole and oozes parameter measurement and centrefuge experiment to determine irreducible water saturation.
10. method as claimed in claim 7, it is characterised in that the saturation exponent of described acquisition every piece described rock core sample Including:
The described rock core sample of partial amt in every kind of Reservoir type is carried out rock-electric test to obtain saturation exponent.
11. the method for claim 1, it is characterised in that the irreducible water saturation bag of the described rock core to be measured of described acquisition Include:
Described rock core to be measured is carried out hole and oozes parameter measurement and centrefuge experiment to determine irreducible water saturation.
CN201610649337.4A 2016-08-10 2016-08-10 The method for determining heterogeneous carbonate reservoir saturation exponent Active CN106093350B (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201610649337.4A CN106093350B (en) 2016-08-10 2016-08-10 The method for determining heterogeneous carbonate reservoir saturation exponent
PCT/CN2017/084304 WO2018028258A1 (en) 2016-08-10 2017-05-15 Method and system for determining heterogeneous carbonate rock reservoir saturation index
US16/095,432 US11112527B2 (en) 2016-08-10 2017-05-15 Method and system for determining heterogeneous carbonate reservoir saturation exponent

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610649337.4A CN106093350B (en) 2016-08-10 2016-08-10 The method for determining heterogeneous carbonate reservoir saturation exponent

Publications (2)

Publication Number Publication Date
CN106093350A true CN106093350A (en) 2016-11-09
CN106093350B CN106093350B (en) 2018-05-04

Family

ID=57456241

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610649337.4A Active CN106093350B (en) 2016-08-10 2016-08-10 The method for determining heterogeneous carbonate reservoir saturation exponent

Country Status (3)

Country Link
US (1) US11112527B2 (en)
CN (1) CN106093350B (en)
WO (1) WO2018028258A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018028258A1 (en) * 2016-08-10 2018-02-15 中国石油天然气股份有限公司 Method and system for determining heterogeneous carbonate rock reservoir saturation index
CN107807078A (en) * 2017-09-13 2018-03-16 中国石油天然气股份有限公司 Rock core mobile water saturation online test method under formation condition
CN111079313A (en) * 2019-12-31 2020-04-28 成都理工大学 Carbonate reservoir information classification processing method and information data processing terminal
CN113640456A (en) * 2020-05-11 2021-11-12 中国石油天然气股份有限公司 Method and device for determining gas pipeline sample

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108647417B (en) * 2018-04-28 2022-11-01 中石化石油工程技术服务有限公司 Simple method for determining gas saturation of shale gas reservoir
CN111441765B (en) * 2019-01-16 2023-12-26 中国石油天然气股份有限公司 Experimental method and device for evaluating air gravity displacement potential of fractured tight oil reservoir
CN112049623B (en) * 2019-06-05 2024-05-28 中国石油天然气股份有限公司 Method and device for judging water outlet of medium-base volcanic bound water
CN111425189B (en) * 2019-06-13 2024-01-30 中国石油天然气股份有限公司 Quality evaluation method and device for ultra-deep fracture-cavity carbonate reservoir and storage medium
CN112627810B (en) * 2019-09-24 2023-08-22 中国石油天然气股份有限公司 Method for detecting karst cave ratio in carbonate reservoir
CN113031064B (en) * 2019-12-25 2022-11-04 中国石油大学(北京) Method, device and equipment for identifying oil layer with low oil saturation
CN113029892B (en) * 2020-03-17 2022-12-13 中国海洋石油集团有限公司 Method for evaluating reasonability of oil-water relative permeability curve based on regional statistical rule
CN111894568A (en) * 2020-08-04 2020-11-06 中国地质大学(北京) Digital core analysis method for fractured carbonate reservoir saturation model
CN114086938B (en) * 2020-08-05 2024-03-08 中国石油化工股份有限公司 Method for predicting gas saturation of heterogeneous sandstone reservoir
CN114428089B (en) * 2020-09-15 2023-11-28 中国石油化工股份有限公司 Clay hydrocarbon source rock classification rock type evaluation method and device, electronic equipment and medium
CN114280686B (en) * 2020-09-27 2024-05-28 中国石油天然气股份有限公司 Rock core physical property analysis method and equipment for carbonate reservoir
CN112304849A (en) * 2020-11-13 2021-02-02 北京拓普莱博油气田开发技术研究院 Method for testing gas permeability of irregular-shaped loose rock sample
CN112595566B (en) * 2020-12-03 2023-05-16 四川蜀科仪器有限公司 Centrifugal dewatering method for compact rock core
CN112924360B (en) * 2021-02-06 2022-05-13 西安石油大学 Method for quantitatively evaluating homogeneity degree of rock core
CN113586043B (en) * 2021-07-19 2023-10-27 中国石油天然气股份有限公司 Method and device for determining saturation parameters of bound water and computer equipment
CN113969777A (en) * 2021-11-17 2022-01-25 中国石油天然气股份有限公司 Reservoir oil saturation determining method, storage medium and electronic equipment
CN114112841B (en) * 2021-11-19 2023-06-09 中国石油大学(华东) Calculation method for irreducible water saturation of tight sandstone reservoir
CN114089421B (en) * 2021-12-06 2023-03-21 中国矿业大学 Oil and gas reservoir heterogeneity analysis method
CN114487341A (en) * 2021-12-17 2022-05-13 成都理工大学 Rock electricity water saturation determination system and method suitable for shale and application
CN114544459B (en) * 2022-02-09 2024-04-19 武汉左晟检测技术有限公司 Rock core oil-water permeability curve construction method based on reservoir parameters
CN116774279B (en) * 2023-06-14 2024-04-16 西南石油大学 Shale water saturation prediction method based on gas-water balance lithology combination unit

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806215A (en) * 2010-03-05 2010-08-18 中国石油集团川庆钻探工程有限公司 Method for judging reservoir fluid type with irreducible water saturation data
CN102243196A (en) * 2011-04-26 2011-11-16 中国石油大学(华东) Method for quickly detecting saturation exponent n in laboratory
CN102434152A (en) * 2011-12-05 2012-05-02 中国石油天然气股份有限公司 Method for calculating oil saturation of reservoir
CN103422853A (en) * 2012-05-24 2013-12-04 中国石油天然气集团公司 Method and device for determining water saturation of stratum
WO2013180593A1 (en) * 2012-06-01 2013-12-05 Общество С Ограниченной Ответственностью "Сплит" Method for determining an oil-and-gas saturation index by a geophysical well-investigation system on the basis of pulsed neutron logging methods
CN104278989A (en) * 2013-07-02 2015-01-14 中国石油天然气股份有限公司 Method for obtaining saturability index of low porosity and low permeability reservoir
CN105114064A (en) * 2015-08-04 2015-12-02 中国石油天然气股份有限公司 Method for determining saturation of compact sandstone reservoir

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4790180A (en) * 1988-02-16 1988-12-13 Mobil Oil Corporation Method for determining fluid characteristics of subterranean formations
FR2844360B1 (en) * 2002-09-11 2004-10-15 Inst Francais Du Petrole METHOD FOR DETERMINING THE RESISTIVITY INDEX BASED ON THE WATER SATURATION OF CERTAIN ROCKS OF COMPLEX POROSITY
US7133779B2 (en) * 2003-05-08 2006-11-07 Schlumberger Technology Corporation Automated borehole geology and petrophysics interpretation using image logs
EP1672392B1 (en) * 2004-12-17 2009-10-28 Services Petroliers Schlumberger Method for determining the water saturation of an underground formation
US8805616B2 (en) * 2010-12-21 2014-08-12 Schlumberger Technology Corporation Method to characterize underground formation
CN102175832B (en) 2011-01-10 2013-10-16 中国石油天然气股份有限公司 Method for determining optimal saturation computing model for typical reservoir
CA2862951C (en) * 2011-10-21 2018-10-30 Saudi Arabian Oil Company Methods, computer readable medium, and apparatus for determining well characteristics and pore architecture utilizing conventional well logs
CN103543474A (en) * 2012-07-16 2014-01-29 中国石油化工股份有限公司 Hydrocarbon saturation evaluation method based on non-conductive pore model
MX361132B (en) * 2013-03-01 2018-11-28 Halliburton Energy Services Inc Downhole differentiation of light oil and oil-based filtrates by nmr with oleophilic nanoparticles.
CN105445441B (en) * 2015-12-03 2017-06-13 中国石油天然气股份有限公司 Reservoir gas-bearing saturation degree determines method and device
CN106093350B (en) * 2016-08-10 2018-05-04 中国石油天然气股份有限公司 The method for determining heterogeneous carbonate reservoir saturation exponent

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806215A (en) * 2010-03-05 2010-08-18 中国石油集团川庆钻探工程有限公司 Method for judging reservoir fluid type with irreducible water saturation data
CN102243196A (en) * 2011-04-26 2011-11-16 中国石油大学(华东) Method for quickly detecting saturation exponent n in laboratory
CN102434152A (en) * 2011-12-05 2012-05-02 中国石油天然气股份有限公司 Method for calculating oil saturation of reservoir
CN103422853A (en) * 2012-05-24 2013-12-04 中国石油天然气集团公司 Method and device for determining water saturation of stratum
WO2013180593A1 (en) * 2012-06-01 2013-12-05 Общество С Ограниченной Ответственностью "Сплит" Method for determining an oil-and-gas saturation index by a geophysical well-investigation system on the basis of pulsed neutron logging methods
CN104278989A (en) * 2013-07-02 2015-01-14 中国石油天然气股份有限公司 Method for obtaining saturability index of low porosity and low permeability reservoir
CN105114064A (en) * 2015-08-04 2015-12-02 中国石油天然气股份有限公司 Method for determining saturation of compact sandstone reservoir

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗少成 等: "致密砂岩储层饱和度指数n 计算方法研究", 《西南石油大学学报(自然科学版)》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018028258A1 (en) * 2016-08-10 2018-02-15 中国石油天然气股份有限公司 Method and system for determining heterogeneous carbonate rock reservoir saturation index
US11112527B2 (en) 2016-08-10 2021-09-07 Petrochina Company Limited Method and system for determining heterogeneous carbonate reservoir saturation exponent
CN107807078A (en) * 2017-09-13 2018-03-16 中国石油天然气股份有限公司 Rock core mobile water saturation online test method under formation condition
CN107807078B (en) * 2017-09-13 2019-11-08 中国石油天然气股份有限公司 Rock core mobile water saturation online test method under formation condition
CN111079313A (en) * 2019-12-31 2020-04-28 成都理工大学 Carbonate reservoir information classification processing method and information data processing terminal
CN113640456A (en) * 2020-05-11 2021-11-12 中国石油天然气股份有限公司 Method and device for determining gas pipeline sample
CN113640456B (en) * 2020-05-11 2024-05-28 中国石油天然气股份有限公司 Method and device for determining gas pipeline sample

Also Published As

Publication number Publication date
US20190137648A1 (en) 2019-05-09
US11112527B2 (en) 2021-09-07
CN106093350B (en) 2018-05-04
WO2018028258A1 (en) 2018-02-15

Similar Documents

Publication Publication Date Title
CN106093350B (en) The method for determining heterogeneous carbonate reservoir saturation exponent
CN109838230B (en) Quantitative evaluation method for oil reservoir water flooded layer
CN102434152B (en) Method for calculating oil saturation of reservoir
CN104564041B (en) Hyposmosis clastic reservoir rock efficiency evaluation method based on exploitation permeability limits
CN109653725B (en) Mixed reservoir flooding degree logging interpretation method based on sedimentary microfacies and rock facies
CN103352693B (en) A kind of output prediction system based on oily content and method thereof
CN104847340B (en) A kind of logging in water flooded layer quantitative evaluation method
CN106323835B (en) The method for determining heterogeneous carbonate reservoir cementation factor
CN102565858B (en) Method for calculating water saturation of porous medium
CN104750884B (en) Shale rich accumulation of oil and gas index quantitative evaluation method based on multifactor nonlinear regression
CN103198363B (en) Reservoir stratum gas production amount prediction method and device based on CT pore analysis
CN102175832A (en) Method for determining optimal saturation computing model for typical reservoir
CN105445441B (en) Reservoir gas-bearing saturation degree determines method and device
Zifei et al. A study on remaining oil distribution in a carbonate oil reservoir based on reservoir flow units
CN104594888B (en) A kind of Identification Method of Oil based on compact reservoir conduction Factor Experiment
CN107102377B (en) The method of quantitative forecast tight sand favorable oil/gas exploration area
CN111963159A (en) Method for identifying fluid properties of conglomerate compact oil reservoir
CN109031436A (en) A kind of hypotonic loose sand reservoir water saturation evaluation method of low-resistance based under Reservoir Classification
CN105401937B (en) A kind of saturation exponent prediction technique based on pore structure
RU2009143585A (en) METHOD FOR DEVELOPING AN INHOMOGENEOUS MASSIVE OR MULTIPLAST GAS-OIL OR OIL AND GAS-CONDENSATE DEPOSIT
CN105240003A (en) Ideal formation factor determination method, hydrocarbon saturation determination method and porosity effectiveness evaluation method
Li et al. In situ estimation of relative permeability from resistivity measurements
CN105604548B (en) A kind of formation oil based on oil base drilling fluid sentences knowledge method
CN113969781B (en) Sea-phase carbonate microporous pore structure cause low-resistance oil layer identification method
KR101405298B1 (en) Evaluation of high or low permeability of fractured rock using well head losses from step-drawdown test

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
GR01 Patent grant