CN106405638A - Method and apparatus for determining gas saturation - Google Patents
Method and apparatus for determining gas saturation Download PDFInfo
- Publication number
- CN106405638A CN106405638A CN201510452356.3A CN201510452356A CN106405638A CN 106405638 A CN106405638 A CN 106405638A CN 201510452356 A CN201510452356 A CN 201510452356A CN 106405638 A CN106405638 A CN 106405638A
- Authority
- CN
- China
- Prior art keywords
- impedance
- petrofacies
- ray elastic
- elastic impedance
- angle
- 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
Links
Abstract
The disclosure provides a method and apparatus for determining the gas saturation. The method comprises the steps of performing ray elastic impedance inversion of data bodies superimposed on different angles to obtain ray elastic impedance inversion bodies of a plurality of angles; calculating an AVO impedance body according to the extracted well side data and logging lithofacies; dividing an impedance range of the AVO impedance body and conducting calculation to form a seismic lithofacies body; performing the ray elastic impedance inversion again to obtain a reservoir elastic parameter body; under the control of the seismic lithofacies body, conducting scattered point intersection analysis of the reservoir elastic parameter body to obtain an elastic parameter body sensitive to the gas saturation; extracting the well side data to conduct compaction analysis to form a dry rock skeleton body; and according to the fluid volume module and fluid density calculated by the use of the dry rock skeleton body, conducting quantitative calculation of the gas saturation. The disclosure realizes the exclusion of multiplicity of solution of an existing pre-stack inversion gas-bearing prediction method and accurately predicts the gas saturation effect in an effective reservoir.
Description
Technical field
The present invention relates to geophysical exploration seismic interpretation FIELD OF THE INVENTIONThe, it is based on rock more particularly, to a kind of
The gas saturation that phase body controls determines method and system.
Background technology
Since twentieth century nineties, China is added to natural gas exploration dynamics, and proved reserves are quickly steady every year
Step increases, but its discovery degree is still very low, and only 6.7%.With the increase of exploration and Exploitation degree,
Evaluating reservoir and gas distribution prediction become the core in work.How to obtain from geological data saturation degree,
The reservoir physical parameters such as porosity are by the key of reservoir pore space fluid prediction.At present, entered using geological data
The qualitative forecasting substantially or for fluid contained by reservoir for the row petroleum-gas prediction.
But the complexity of subsurface reservoir lithology determine can not be with an attribute or rock physics equations
Subsurface reservoir distribution, the reservoir that the fluid factor being obtained by correlation predictive method or inverting are obtained to be described
Physical parameter all can be disturbed by non-effective reservoir, and precision of prediction is low.
Content of the invention
The present disclosure proposes a kind of high gas saturation of the interference of non-effective reservoir, precision of prediction of excluding determines
Method and system.
One side according to the disclosure it is proposed that a kind of gas saturation determines method, including：
The data volume of bipartition angle degree superposition carries out the ray elasticity resistance that ray elastic impedance inverting obtains multiple angles
Anti- invertomer；
Petrofacies division is carried out to the ray elastic impedance invertomer of multiple angles according to lithology and porosity value, obtains
To log-petrofacies；
Carry out being calculated AVO impedance body according to the well lie data extracted and log-petrofacies；
Division to the impedance codomain of AVO impedance body, calculates and forms seismic facies body；
Multiple ray elastic impedance invertomers are carried out ray elastic impedance inverting and obtains reservoir elastic parameter body；
Under the control of seismic facies body, reservoir elastic parameter body is carried out with scatterplot cross analysis, obtain to containing
The sensitive elastic parameter body of gas saturation；
Well lie data is extracted to sensibility elasticity parameter body, carries out compaction analysis and form dry rock matrix body；
According to using the calculated fluid volume module of dry rock matrix body and fluid density, quantitative calculating contains
Gas saturation.
According to the disclosure in a first aspect, in the first possible implementation of the disclosure, described basis
Lithology and porosity value carry out petrofacies division to the ray elastic impedance invertomer of multiple angles, obtain rock of logging well
Mutually include：
Extract low-angle ray elastic impedance from the ray elastic impedance invertomer of multiple angles and wide-angle is penetrated
Linear elasticity impedance；
Extract the porosity information in low-angle ray elastic impedance and the lithology in wide-angle ray elastic impedance
Information；
According to the porosity information in lithology and porosity value, low-angle ray elastic impedance and wide-angle ray
Lithological information in elastic impedance carries out petrofacies division, obtains log-petrofacies.
According to the first possible implementation of disclosure first aspect, in the possible reality of the second of the disclosure
In existing mode, the described well lie data according to extraction and the log-petrofacies marking off carry out being calculated AVO
Impedance body includes：
According to the well lie data extracting and the log-petrofacies marking off, carry out scatterplot cross analysis；
Method using attribute projection is entered to low-angle ray elastic impedance body and wide-angle ray elastic impedance body
Row coordinate rotates, and obtains coordinate transformation formula；
Result according to scatterplot cross analysis and coordinate transformation formula are to low-angle ray elastic impedance and wide-angle
Ray elastic impedance carries out being calculated AVO impedance body.
According to the disclosure in a first aspect, in the third possible implementation of the disclosure, described to quick
Sense elastic parameter body extracts well lie data, carries out compaction analysis and form dry rock matrix body including：
Well lie data is extracted to sensibility elasticity parameter body；
Compaction analysis are carried out to well lie data；
Result according to compaction analysis obtains the attribute change trend of dry rock matrix；
Many attributes well interpolation is carried out according to the attribute change trend of dry rock matrix and forms dry rock matrix body.
According to the disclosure in a first aspect, the described elastic parameter body sensitive to gas saturation include density and
Bulk modulus.
According to the second aspect of the disclosure, a kind of gas saturation determines device, including：
First inverting module, the data volume for the superposition of bipartition angle degree carries out ray elastic impedance inverting, obtains
The ray elastic impedance invertomer of multiple angles；
Petrofacies division module, for the ray elastic impedance inverting to multiple angles according to lithology and porosity value
Body carries out petrofacies division, obtains log-petrofacies；
Impedance body computing module, based on carrying out according to the well lie data extracted and the log-petrofacies marking off
Calculation obtains AVO impedance body；
Petrofacies body computing module, for the division of the impedance codomain to AVO impedance body, calculates and forms earthquake rock
Xiang Ti；
Second inverting module, for multiple ray elastic impedance invertomers are carried out ray elastic impedance inverting,
Obtain reservoir elastic parameter body；
Cross analysis module, for, under the control of seismic facies body, carrying out scatterplot to reservoir elastic parameter body
Cross analysis, obtain the elastic parameter body sensitive to gas saturation；
Compaction analysis module, for extracting well lie data to sensibility elasticity parameter body, carries out compaction analysis shape
Become dry rock matrix body；
Saturation computation module, for according to using the calculated fluid volume module of dry rock matrix body and
Fluid density, quantitative calculating gas saturation.
According to the second aspect of the disclosure, in the first possible implementation in second aspect, described petrofacies
Division module includes：
First extracting sub-module, penetrates for extracting low-angle from the ray elastic impedance invertomer of multiple angles
Linear elasticity impedance and wide-angle ray elastic impedance；
Second extracting sub-module, for extracting porosity information and wide-angle in low-angle ray elastic impedance
Lithological information in ray elastic impedance；
Divide submodule, for according to the porosity in lithology and porosity value, low-angle ray elastic impedance
Lithological information in information and wide-angle ray elastic impedance carries out petrofacies division, obtains log-petrofacies.
The first possible implementation of second aspect according to the disclosure is possible in second aspect second
In implementation, described impedance body computing module includes：
First analysis submodule, for according to the well lie data extracting and the log-petrofacies marking off, carrying out
Scatterplot cross analysis；
Coordinate rotates submodule, for the method that projected using attribute to low-angle ray elastic impedance body and big
Angle ray elastic impedance body carries out coordinate rotation, obtains coordinate transformation formula；
Calculating sub module, for the result according to scatterplot cross analysis and coordinate transformation formula to low-angle ray
Elastic impedance and wide-angle ray elastic impedance carry out being calculated AVO impedance body.
According to the second aspect of the disclosure, in the third possible implementation in second aspect, described compacting
Analysis module includes：
Data extracting sub-module, for extracting well lie data to sensibility elasticity parameter body；
Compaction analysis submodule, for carrying out compaction analysis to well lie data；
Acquisition submodule, obtains the attribute change trend of dry rock matrix for the result according to compaction analysis；
Interpolation submodule, forms for carrying out many attributes well interpolation according to the attribute change trend of dry rock matrix
Dry rock matrix body.
According to the second aspect of the disclosure, in the 4th kind of possible implementation of second aspect, described to containing
The sensitive elastic parameter body of gas saturation includes density and bulk modulus.
The each side of the disclosure crosses by using ray elastic impedance and forms new AVO Impedance Inversion earthquake rock
Xiang Ti, predicting reservoir is favorably distributed, and excludes non-reservoir impact, carry out quantitative gassiness and satisfy under petrofacies body controls
With degree conversion, in conversion process, using the simplified style of Gassmann equation, porous rocks are interpreted as drying
Skeleton and pore-fluid two parts, are removing drying nest, between the fluid density obtaining and gas saturation
There is linear relationship it is achieved that excluding the multi-solution of existing prestack inversion gas distribution prediction method, accurately in advance
Survey the gas saturation effect in Effective Reservoirs.
Brief description
By combining accompanying drawing, disclosure illustrative embodiments are described in more detail, the disclosure above-mentioned
And other purpose, feature and advantage will be apparent from, wherein, in disclosure illustrative embodiments
In, identical reference number typically represents same parts.
Fig. 1 shows that a kind of gas saturation that the disclosure one exemplary embodiment provides determines the flow process of method
Schematic diagram.
Fig. 2 shows a kind of schematic flow sheet of the step 102 shown in the disclosure one exemplary embodiment.
Fig. 3 shows a kind of schematic flow sheet of the step 103 shown in the disclosure one exemplary embodiment.
Fig. 4 shows a kind of schematic flow sheet of the step 107 shown in the disclosure one exemplary embodiment.
Fig. 5 shows that a kind of gas saturation that the disclosure one exemplary embodiment provides determines the structure of device
Schematic diagram.
Fig. 6 shows that a kind of structure of the petrofacies division module 202 shown in the disclosure one exemplary embodiment is shown
It is intended to.
Fig. 7 shows a kind of structure of the impedance body computing module 203 shown in the disclosure one exemplary embodiment
Schematic diagram.
Fig. 8 shows that a kind of structure of the compaction analysis module 207 shown in the disclosure one exemplary embodiment is shown
It is intended to.
Fig. 9 show it relates to petrophysical model schematic diagram.
Specific embodiment
It is more fully described the preferred embodiment of the disclosure below with reference to accompanying drawings.Although showing in accompanying drawing
The preferred embodiment of the disclosure, however, it is to be appreciated that may be realized in various forms the disclosure and should be by
Embodiments set forth herein is limited.On the contrary, these embodiments are provided so that the disclosure is more saturating
Thorough and complete, and the scope of the present disclosure intactly can be conveyed to those skilled in the art.
Fig. 1 shows the flow chart that gas saturation according to an embodiment of the invention determines method, such as
Shown in Fig. 1, the method comprises the following steps.
Step 101, the data volume of bipartition angle degree superposition carries out ray elastic impedance inverting and obtains penetrating of multiple angles
Linear elasticity Impedance Inversion body.
Described step 101 can be implemented as in the disclosed embodiments, due to the ray elasticity of multiple angles
Impedance Inversion body includes low-angle ray elastic impedance and wide-angle ray elastic impedance, and low-angle ray
Elastic impedance contains porosity information, and wide-angle ray elastic impedance contains lithological information it is possible to pass through
Therefore rock can be distinguished by wide-angle ray elastic impedance (REI25) and low-angle ray elastic impedance (REI5)
Phase.Depending on the number range of angle typically sees the angular range of real data, for example general prestack road collection is incident
Angular region is 0～30, is equally divided into three parts, and low-angle is 5, and middle angle is 15, and wide-angle is 25.
Step 102, carries out petrofacies according to lithology and porosity value to the ray elastic impedance invertomer of multiple angles
Divide, obtain log-petrofacies.
Step 103, carries out being calculated AVO impedance body according to the well lie data extracted and log-petrofacies.
Described step 103 can be implemented as in the disclosed embodiments, the in figure that crosses scatterplot shape, permissible
Explained by log-petrofacies and obtain, such as can make triangular representation mud stone, square is dry sandstone reservoir, circular
For sandstone reservoir etc..
The number range of petrofacies body should be petrofacies classification codomain (1,2,3 ... n).In the present embodiment, two angles
The ray elastic impedance of degree passes through Coordinate Conversion, forms new impedance body (AVOImp), by AVOImp with
The density figure that crosses shows as can be seen that reservoir, dried layer and mud stone type are distinguished in AVO impedance in codomain, logical
Cross threshold value division mode, can by AVO impedance conversion become lithologic body (define 1 mud stone, 2 dried layer, 3
Reservoir).
Step 104, the division to the impedance codomain of AVO impedance body, calculate and form seismic facies body.
Described step 104 can be implemented as in the disclosed embodiments, petrofacies body LITH and AVO impedance
Body should meet equation below：
Value1 and value2 is the minimum of a value of AVOImp, the maximum distinguishing mud stone；Value1 and value2
For distinguishing minimum of a value, the maximum of the AVOImp of dried layer；Value1 and value2 is the AVOImp distinguishing reservoir
Minimum of a value, maximum.
Wherein, AVO impedance and the big low-angle ray elastic impedance body of input meet equation below：
AVOImp=A*X+B*Y+C
In formula, X is low-angle ray elastic impedance body, and Y is wide-angle ray elastic impedance body, A, B and
C is fitting parameter.
Step 105, multiple ray elastic impedance invertomers is carried out ray elastic impedance inverting and obtains reservoir elasticity
Parameter body.
Described step 105 can be implemented as in the disclosed embodiments, under the control of petrofacies body, reservoir
Elastic parameter can be described with Gassmann equation, as shown in figure 9, when reservoir is for saturated rock：
μs=μd
ρs=(1- φ) ρm+φρf
In formula：KsIt is the bulk modulus of saturated rock；KdIt is the bulk modulus of dry rock skeleton；KmIt is rock
The bulk modulus of stone solid matrix；KfIt is the bulk modulus of pore-fluid；φ is porosity；μsIt is saturated rock
The modulus of shearing of stone；μdIt is the modulus of shearing of dry rock skeleton；ρsIt is the density of saturated rock.
Here Gassmann equation pore-fluid relates only to monophasic fluid.
Step 106, under the control of seismic facies body, carries out scatterplot cross analysis to reservoir elastic parameter body,
Obtain the elastic parameter body sensitive to gas saturation.
Described step 106 can be implemented as in the disclosed embodiments, described sensitive to gas saturation
Elastic parameter body includes density and bulk modulus, if it is assumed that fluid comprises gas and water two-phase, fluid density can
Write as the weighted sum of two-phase fluid, the bulk modulus of fluid can be represented with Wood equation：
Or adopting Brie experience fluid mixed equation, e is empirical value, typically takes 2.5, value is less to approach Wood
Equation, approaches more greatly patch saturation relation fluid mixed equation：
Kf=(Kw-Kg)Sg e+Kg
Fluid density ρfCan be written as：
ρf=(1-Sg)ρw+Sgρg
By three above formula, when when containing two-phase fluid in hole, the wherein saturation degree of one phase flow body
Can be expressed with non-linear formula by fluid modulus, fluid density is expressed with linear formula.
Step 107, extracts well lie data to sensibility elasticity parameter body, carries out compaction analysis and form dry rock bone
Support body.
Described step 107 can be implemented as in the disclosed embodiments, extracts well to sensibility elasticity parameter body
Lie data, due to Kd、ρdRelevant with porosity and diagenesis, with depth, pressure trend change, permissible
By log compaction analysis, form dry rock matrix body.
Step 108, according to using the calculated fluid volume module of dry rock matrix body and fluid density, fixed
Amount calculates gas saturation.
By rearranging to Gassmann equation, fluid density and fluid modulus can be write as and be satisfied
The amount relevant with porosity with rock, skeleton, matrix：
As can be seen from the above equation it is desirable to obtain fluid density and fluid modulus, the physical quantity known is needed to have：
The bulk modulus of dry rock and density Kd、ρd, the bulk modulus of solid matrix and density Km、ρm, and hole
Degree φ；Because porosity φ is with change in depth, in the classification of specific thing phase, may be defined as constant；Km、ρmRow
Except relation with shale content, unrelated with porosity, not with change in depth, rock sample data test and tabling look-up obtains
?.
The present embodiment crosses by using ray elastic impedance and forms new AVO Impedance Inversion seismic facies body,
Predicting reservoir is favorably distributed, and excludes non-reservoir impact, carry out quantitative gas saturation and turn under petrofacies body controls
Change, in conversion process, using the simplified style of Gassmann equation, by porous rocks be interpreted as being dried skeleton with
Pore-fluid two parts, are removing drying nest, there is line between the fluid density obtaining and gas saturation
Sexual intercourse, it is achieved that excluding the multi-solution of existing prestack inversion gas distribution prediction method, is accurately predicted effectively
Gas saturation effect in reservoir.
Because low-angle ray elastic impedance contains porosity information, wide-angle ray elastic impedance contains lithology
Information, therefore can distinguish petrofacies by large and small angle ray elastic impedance, as shown in Fig. 2 described step
Rapid 102 include：
Step 1021, extracts low-angle ray elastic impedance from the ray elastic impedance invertomer of multiple angles
With wide-angle ray elastic impedance；
Step 1022, extracts the porosity information in low-angle ray elastic impedance and the elasticity resistance of wide-angle ray
Lithological information in anti-；
Step 1023, according to the porosity information in lithology and porosity value, low-angle ray elastic impedance and
Lithological information in wide-angle ray elastic impedance carries out petrofacies division, obtains log-petrofacies.
The present embodiment passes through porosity information and wide-angle ray elastic impedance in low-angle ray elastic impedance
In lithological information carry out petrofacies division, log-petrofacies can be accurately obtained.
In order to the interference of mud stone and dried layer can be excluded on AVO impedance, directly find porosity with
The relation of AVO impedance, as shown in figure 3, described step 103 includes：
Step 1031, according to the well lie data extracting and the log-petrofacies marking off, carries out scatterplot and crosses point
Analysis；
Step 1032, using attribute projection method to low-angle ray elastic impedance body and wide-angle ray bullet
Property impedance body carries out coordinate rotation, obtains coordinate transformation formula；
Step 1033, the result according to scatterplot cross analysis and coordinate transformation formula are to the elasticity resistance of low-angle ray
Anti- and wide-angle ray elastic impedance carries out being calculated AVO impedance body.
The present embodiment passes through, according to the well lie data extracting and the log-petrofacies marking off, to carry out scatterplot and cross
Analysis, using attribute projection method to low-angle ray elastic impedance body and wide-angle ray elastic impedance body
Carry out coordinate rotation, obtain coordinate transformation formula, the result according to scatterplot cross analysis and coordinate transformation formula
Carry out being calculated AVO impedance body to low-angle ray elastic impedance and wide-angle ray elastic impedance, permissible
Realize excluding the interference of mud stone and dried layer using calculated AVO impedance, directly find porosity and AVO
The relation of impedance.
In order to remove the rock matrix of compaction, obtain the gas saturation body of reflection fluid, reach
Quantitative forecast reservoir gas-bearing property purpose, as shown in figure 4, described step 107 includes：
Step 1071, extracts well lie data to sensibility elasticity parameter body；
Step 1072, carries out compaction analysis to well lie data；
Step 1073, the result according to compaction analysis obtains the attribute change trend of dry rock matrix；
Step 1074, carries out many attributes well interpolation according to the attribute change trend of dry rock matrix and forms dry rock
Skeleton body.
The present embodiment, by extracting well lie data to sensibility elasticity parameter body, is compacted to well lie data
Analysis, the result according to compaction analysis obtains the attribute change trend of dry rock matrix, according to dry rock matrix
Attribute change trend carry out many attributes well interpolation and form dry rock matrix body it is achieved that removing compaction
Rock matrix, obtains the gas saturation body of reflection fluid.
In addition, as the realization to the various embodiments described above, the embodiment of the present disclosure additionally provides a kind of gassiness saturation
Degree determines device, as shown in figure 5, including：First inverting module 201, petrofacies division module 202, impedance
Body computing module 203, petrofacies body computing module 204, the second inverting module 205, cross analysis module 206,
Compaction analysis module 207 and saturation computation module 208.
First inverting module 201, the data volume for the superposition of bipartition angle degree carries out ray elastic impedance inverting, obtains
Ray elastic impedance invertomer to multiple angles.
Petrofacies division module 202, for anti-to the ray elastic impedance of multiple angles according to lithology and porosity value
Drill body and carry out petrofacies division, obtain log-petrofacies.
Impedance body computing module 203, for carrying out according to the well lie data extracted and the log-petrofacies marking off
It is calculated AVO impedance body.
Petrofacies body computing module 204, for the division of the impedance codomain to AVO impedance body, calculates and forms ground
Shake petrofacies body.
Second inverting module 205, anti-for multiple ray elastic impedance invertomers are carried out ray elastic impedance
Drill, obtain reservoir elastic parameter body.
Cross analysis module 206, for, under the control of seismic facies body, dissipating to reservoir elastic parameter body
Point cross analysis, obtain the elastic parameter body sensitive to gas saturation, described sensitive to gas saturation
Elastic parameter body includes density and bulk modulus.
Compaction analysis module 207, for extracting well lie data to sensibility elasticity parameter body, carries out compaction analysis
Form dry rock matrix body；
Saturation computation module 208, for according to using the calculated fluid volume module of dry rock matrix body
And fluid density, quantitative calculating gas saturation.
The present embodiment crosses by using ray elastic impedance and forms new AVO Impedance Inversion seismic facies body,
Predicting reservoir is favorably distributed, and excludes non-reservoir impact, carry out quantitative gas saturation and turn under petrofacies body controls
Change, in conversion process, using the simplified style of Gassmann equation, by porous rocks be interpreted as being dried skeleton with
Pore-fluid two parts, are removing drying nest, there is line between the fluid density obtaining and gas saturation
Sexual intercourse, it is achieved that excluding the multi-solution of existing prestack inversion gas distribution prediction method, is accurately predicted effectively
Gas saturation effect in reservoir.
Because low-angle ray elastic impedance contains porosity information, wide-angle ray elastic impedance contains lithology
Information, therefore can distinguish petrofacies by large and small angle ray elastic impedance, as shown in fig. 6, described rock
Phase division module 202 includes：First extracting sub-module 2021, the second extracting sub-module 2022 and division submodule
Block 2023.
First extracting sub-module 2021, for extracting little angle from the ray elastic impedance invertomer of multiple angles
Degree ray elastic impedance and wide-angle ray elastic impedance.
Second extracting sub-module 2022, for extracting the porosity information in low-angle ray elastic impedance and big
Lithological information in angle ray elastic impedance.
Divide submodule 2023, for according to the hole in lithology and porosity value, low-angle ray elastic impedance
Lithological information in porosity information and wide-angle ray elastic impedance carries out petrofacies division, obtains log-petrofacies.
The present embodiment passes through porosity information and wide-angle ray elastic impedance in low-angle ray elastic impedance
In lithological information carry out petrofacies division, log-petrofacies can be accurately obtained.
In order to the interference of mud stone and dried layer can be excluded on AVO impedance, directly find porosity with
The relation of AVO impedance, as shown in fig. 7, described impedance body computing module 203 includes：First analysis submodule
Block 2031, coordinate rotation submodule 2032 and calculating sub module 2033.
First analysis submodule 2031, for the well lie data according to extraction and the log-petrofacies marking off,
Carry out scatterplot cross analysis.
Coordinate rotates submodule 2032, for the method that projected using attribute to low-angle ray elastic impedance body
Carry out coordinate rotation with wide-angle ray elastic impedance body, obtain coordinate transformation formula.
Calculating sub module 2033, for the result according to scatterplot cross analysis and coordinate transformation formula to low-angle
Ray elastic impedance and wide-angle ray elastic impedance carry out being calculated AVO impedance body.
The present embodiment passes through, according to the well lie data extracting and the log-petrofacies marking off, to carry out scatterplot and cross
Analysis, using attribute projection method to low-angle ray elastic impedance body and wide-angle ray elastic impedance body
Carry out coordinate rotation, obtain coordinate transformation formula, the result according to scatterplot cross analysis and coordinate transformation formula
Carry out being calculated AVO impedance body to low-angle ray elastic impedance and wide-angle ray elastic impedance, permissible
Realize excluding the interference of mud stone and dried layer using calculated AVO impedance, directly find porosity and AVO
The relation of impedance.
In order to remove the rock matrix of compaction, obtain the gas saturation body of reflection fluid, reach
Quantitative forecast reservoir gas-bearing property purpose, as shown in figure 8, described compaction analysis module 207 includes：Data carries
Take submodule 2071, compaction analysis submodule 2072, acquisition submodule 2073 and interpolation submodule 2074.
Data extracting sub-module 2071, for extracting well lie data to sensibility elasticity parameter body.
Compaction analysis submodule 2072, for carrying out compaction analysis to well lie data.
Acquisition submodule 2073, the attribute change obtaining dry rock matrix for the result according to compaction analysis becomes
Gesture.
Interpolation submodule 2074, for carrying out many attributes well interpolation according to the attribute change trend of dry rock matrix
Form dry rock matrix body.
The present embodiment, by extracting well lie data to sensibility elasticity parameter body, is compacted to well lie data
Analysis, the result according to compaction analysis obtains the attribute change trend of dry rock matrix, according to dry rock matrix
Attribute change trend carry out many attributes well interpolation and form dry rock matrix body it is achieved that removing compaction
Rock matrix, obtains the gas saturation body of reflection fluid.
Application example
For ease of understanding scheme and its effect of the embodiment of the present invention, a concrete application example given below.
It will be understood by those skilled in the art that this example is only for the purposes of understanding the present invention, its any detail is not
It is intended to limit by any way the present invention.
Chuan Xi somewhere Penglaizhen Formation group develops many set gas-bearing formations, illustrates good Exploration Potential.But developing
Cheng Zhong, finds that this low porosity and low permeability Gas Reservoirs water-gas relation is complicated, the ground of different gassiness abundance gas-bearing formations rings
Answer feature difference inconspicuous, lead to conventional gas and oil detection method to be difficult to solve local area air water forecasting problem, for
Local area geology characteristic and seismic data present situation, control lower gas saturation prediction to solve using petrofacies current
Exploitation difficult point.
Research area definition seismic facies have considered deposition and diagenesis, porosity, and combine
Research area petrofacies are marked off three class petrofacies (table 1), obtain ground by inverting by the standard of logging Reservoir Evaluation
Shake petrofacies body.Divided by petrofacies body, I, II class thing phase porosity is defined as
The table 1 research area petrofacies criteria for classifying
| Petrofacies | Reservoir levels | The criteria for classifying |
| I class | 1 class reservoir, 2 class reservoirs | 7<Por<15 |
| II class | Dried layer | Por<7 |
| Group III | Mud stone | Por<7 |
On the basis of petrofacies calculate, crossed Rock physical analysis using the four-dimension, i.e. density (abscissa), in length and breadth
Wave velocity ratio (ordinate), porosity (color) and four kinds of parameters of seismic facies (symbol) are crossed
Analysis, is analyzed by intersection it is recognised that in Favorable Reservoir area, controlling lower density in certain porosity and contain
Gas saturation is linear.
The major parameter of this prediction gassiness abundance by the use of density parameter as us.On this basis by fluid
Density is converted into gas saturation by Gassman equation deformation type.In figure is red, yellow represents gas-bearing formation, green
Color represents air water area, and blue is water layer, and white is reservoir agensis section, predicts the outcome and each testing well result
Substantially identical, effectively demonstrate the applicability of method.
The disclosure can be system, method and/or computer program.Computer program can include
Computer-readable recording medium, containing for making processor realize the computer of various aspects of the disclosure
Readable program instructions.
Computer-readable recording medium can be to keep and store the instruction being used by instruction execution equipment
Tangible device.Computer-readable recording medium for example may be-but not limited to-storage device electric,
Magnetic storage apparatus, light storage device, electromagnetism storage device, semiconductor memory apparatus or above-mentioned any conjunction
Suitable combination.The more specifically example (non exhaustive list) of computer-readable recording medium includes：Portable
Formula computer disks, hard disk, random access memory (RAM), read-only storage (ROM), erasable type can
Program read-only memory (EPROM or flash memory), static RAM (SRAM), portable pressure
Contracting disk read-only storage (CD-ROM), digital versatile disc (DVD), memory stick, floppy disk, machinery are compiled
Decoding apparatus, the punch card being for example stored thereon with instruction or groove internal projection structure and above-mentioned any conjunction
Suitable combination.Computer-readable recording medium used herein above is not construed as instantaneous signal itself, such as
The electromagnetic wave of radio wave or other Free propagations, the electromagnetic wave propagated by waveguide or other transmission mediums
(for example, by the light pulse of fiber optic cables) or by the electric signal of wire transfer.
Computer-readable program instructions as described herein can download to each from computer-readable recording medium
Calculating/processing equipment, or downloaded by network, such as internet, LAN, wide area network and/or wireless network
To outer computer or External memory equipment.Network can include copper transmission cable, Optical Fiber Transmission, wireless biography
Defeated, router, fire wall, switch, gateway computer and/or Edge Server.Each calculates/processes and sets
Adapter in standby or network interface receive computer-readable program instructions from network, and forward this meter
Calculation machine readable program instructions, for being stored in the computer-readable recording medium in each calculating/processing equipment.
Computer program instructions for executing disclosure operation can be assembly instruction, instruction set architecture
(ISA) instruction, machine instruction, machine-dependent instructions, microcode, firmware instructions, condition setup data,
Or the source code write with any combination of one or more programming language or object code, described programming language
Speech includes OO programming language Smalltalk, C++ etc., and the process type programming language of routine
Speech such as " C " language or similar programming language.Computer-readable program instructions can fully be counted in user
Execution on calculation machine, partly on the user computer execution, as independent software kit execution, a part
Part executes or on the remote computer completely on remote computer or server on the user computer
Execution.In the situation being related to remote computer, remote computer can pass through the network bag of any kind
Include LAN (LAN) or wide area network (WAN) is connected to subscriber computer, or it may be connected to outside meter
Calculation machine (for example using ISP come by Internet connection).In certain embodiments, pass through
Using computer-readable program instructions status information come personalized customization electronic circuit, such as FPGA
Circuit, field programmable gate array (FPGA) or programmable logic array (PLA), this electronic circuit can
To execute computer-readable program instructions, thus realizing various aspects of the disclosure.
Stream referring herein to method, device (system) and computer program according to the embodiment of the present disclosure
Journey figure and/or block diagram describe various aspects of the disclosure.It should be appreciated that each side of flow chart and/or block diagram
In frame and flow chart and/or block diagram, the combination of each square frame, can be realized by computer-readable program instructions.
These computer-readable program instructions can be supplied to all-purpose computer, special-purpose computer or other and can compile
The processor of journey data processing equipment, thus produce a kind of machine so that these instructions are by computer
Or other programmable data processing unit computing device when, create in flowchart and/or block diagram
The device of function/action specified in one or more square frames.Can also be these computer-readable program instructions
Storage in a computer-readable storage medium, these instruction make computer, programmable data processing unit and/
Or other equipment works in a specific way, thus, the computer-readable medium of the instruction that is stored with then includes one
Manufacture, it includes function/action specified in one or more of flowchart and/or block diagram square frame
The instruction of various aspects.
Can also computer-readable program instructions be loaded into computer, other programmable data processing unit,
Or so that executing one on computer, other programmable data processing unit or miscellaneous equipment on miscellaneous equipment
Series of operative steps, to produce computer implemented process, so that in computer, other programmable number
According to one or more of the instruction flowchart of execution and/or block diagram side in processing meanss or miscellaneous equipment
Function/action specified in frame.
Flow chart in accompanying drawing and block diagram show the system of multiple embodiments according to the disclosure, method and meter
The architectural framework in the cards of calculation machine program product, function and operation.At this point, flow chart or block diagram
In each square frame can represent the part of a module, program segment or instruction, described module, program segment
Or a part for instruction comprises the executable instruction of one or more logic functions for realizing regulation.Having
In a little realizations as replacement, the function of being marked in square frame can also be with suitable different from marked in accompanying drawing
Sequence occurs.For example, two continuous square frames can essentially execute substantially in parallel, and they sometimes can also
Execute in the opposite order, this is depending on involved function.It is also noted that block diagram and/or flow chart
In each square frame and square frame in block diagram and/or flow chart combination, can be with the function of execution regulation
Or the special hardware based system of action is realizing, or can be with specialized hardware and computer instruction
Combine and to realize.
It is described above the presently disclosed embodiments, described above is exemplary, and non-exclusive,
And it is also not necessarily limited to disclosed each embodiment.In the scope and spirit without departing from illustrated each embodiment
In the case of, many modifications and changes will be apparent from for those skilled in the art.
The selecting it is intended to best explain the principle of each embodiment, practical application or to market of term used herein
In technology technological improvement, or so that other those of ordinary skill of the art is understood that and discloses herein
Each embodiment.
Claims (10)
1. a kind of gas saturation determines method it is characterised in that including：
The data volume of bipartition angle degree superposition carries out the ray elasticity resistance that ray elastic impedance inverting obtains multiple angles
Anti- invertomer；
Petrofacies division is carried out to the ray elastic impedance invertomer of multiple angles according to lithology and porosity value, obtains
To log-petrofacies；
Carry out being calculated AVO impedance body according to the well lie data extracted and log-petrofacies；
Division to the impedance codomain of AVO impedance body, calculates and forms seismic facies body；
Multiple ray elastic impedance invertomers are carried out ray elastic impedance inverting and obtains reservoir elastic parameter body；
Under the control of seismic facies body, reservoir elastic parameter body is carried out with scatterplot cross analysis, obtain to containing
The sensitive elastic parameter body of gas saturation；
Well lie data is extracted to sensibility elasticity parameter body, carries out compaction analysis and form dry rock matrix body；
According to using the calculated fluid volume module of dry rock matrix body and fluid density, quantitative calculating contains
Gas saturation.
2. gas saturation according to claim 1 determine method it is characterised in that described according to rock
Property and porosity value carry out petrofacies division to the ray elastic impedance invertomer of multiple angles, obtain log-petrofacies
Including：
Extract low-angle ray elastic impedance from the ray elastic impedance invertomer of multiple angles and wide-angle is penetrated
Linear elasticity impedance；
Extract the porosity information in low-angle ray elastic impedance and the lithology in wide-angle ray elastic impedance
Information；
According to the porosity information in lithology and porosity value, low-angle ray elastic impedance and wide-angle ray
Lithological information in elastic impedance carries out petrofacies division, obtains log-petrofacies.
3. gas saturation according to claim 2 determines method it is characterised in that described basis carries
The well lie data taking and the log-petrofacies marking off carry out being calculated AVO impedance body and include：
According to the well lie data extracting and the log-petrofacies marking off, carry out scatterplot cross analysis；
Method using attribute projection is entered to low-angle ray elastic impedance body and wide-angle ray elastic impedance body
Row coordinate rotates, and obtains coordinate transformation formula；
Result according to scatterplot cross analysis and coordinate transformation formula are to low-angle ray elastic impedance and wide-angle
Ray elastic impedance carries out being calculated AVO impedance body.
4. gas saturation according to claim 1 determine method it is characterised in that described to sensitivity
Elastic parameter body extracts well lie data, carries out compaction analysis and form dry rock matrix body including：
Well lie data is extracted to sensibility elasticity parameter body；
Compaction analysis are carried out to well lie data；
Result according to compaction analysis obtains the attribute change trend of dry rock matrix；
Many attributes well interpolation is carried out according to the attribute change trend of dry rock matrix and forms dry rock matrix body.
5. gas saturation according to claim 1 determine method it is characterised in that described to gassiness
The sensitive elastic parameter body of saturation degree includes density and bulk modulus.
6. a kind of gas saturation determines device it is characterised in that including：
First inverting module, the data volume for the superposition of bipartition angle degree carries out ray elastic impedance inverting, obtains
The ray elastic impedance invertomer of multiple angles；
Petrofacies division module, for the ray elastic impedance inverting to multiple angles according to lithology and porosity value
Body carries out petrofacies division, obtains log-petrofacies；
Impedance body computing module, based on carrying out according to the well lie data extracted and the log-petrofacies marking off
Calculation obtains AVO impedance body；
Petrofacies body computing module, for the division of the impedance codomain to AVO impedance body, calculates and forms earthquake rock
Xiang Ti；
Second inverting module, for multiple ray elastic impedance invertomers are carried out ray elastic impedance inverting,
Obtain reservoir elastic parameter body；
Cross analysis module, for, under the control of seismic facies body, carrying out scatterplot to reservoir elastic parameter body
Cross analysis, obtain the elastic parameter body sensitive to gas saturation；
Compaction analysis module, for extracting well lie data to sensibility elasticity parameter body, carries out compaction analysis shape
Become dry rock matrix body；
Saturation computation module, for according to using the calculated fluid volume module of dry rock matrix body and
Fluid density, quantitative calculating gas saturation.
7. gas saturation according to claim 6 determines device it is characterised in that described petrofacies are drawn
Sub-module includes：
First extracting sub-module, penetrates for extracting low-angle from the ray elastic impedance invertomer of multiple angles
Linear elasticity impedance and wide-angle ray elastic impedance；
Second extracting sub-module, for extracting porosity information and wide-angle in low-angle ray elastic impedance
Lithological information in ray elastic impedance；
Divide submodule, for according to the porosity in lithology and porosity value, low-angle ray elastic impedance
Lithological information in information and wide-angle ray elastic impedance carries out petrofacies division, obtains log-petrofacies.
8. gas saturation according to claim 7 determines device it is characterised in that described impedance body
Computing module includes：
First analysis submodule, for according to the well lie data extracting and the log-petrofacies marking off, carrying out
Scatterplot cross analysis；
Coordinate rotates submodule, for the method that projected using attribute to low-angle ray elastic impedance body and big
Angle ray elastic impedance body carries out coordinate rotation, obtains coordinate transformation formula；
Calculating sub module, for the result according to scatterplot cross analysis and coordinate transformation formula to low-angle ray
Elastic impedance and wide-angle ray elastic impedance carry out being calculated AVO impedance body.
9. gas saturation according to claim 6 determines device it is characterised in that described compacting divides
Analysis module includes：
Data extracting sub-module, for extracting well lie data to sensibility elasticity parameter body；
Compaction analysis submodule, for carrying out compaction analysis to well lie data；
Acquisition submodule, obtains the attribute change trend of dry rock matrix for the result according to compaction analysis；
Interpolation submodule, forms for carrying out many attributes well interpolation according to the attribute change trend of dry rock matrix
Dry rock matrix body.
10. gas saturation according to claim 6 determine device it is characterised in that described to containing
The sensitive elastic parameter body of gas saturation includes density and bulk modulus.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510452356.3A CN106405638B (en) | 2015-07-28 | 2015-07-28 | A kind of gas saturation determines method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510452356.3A CN106405638B (en) | 2015-07-28 | 2015-07-28 | A kind of gas saturation determines method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106405638A true CN106405638A (en) | 2017-02-15 |
| CN106405638B CN106405638B (en) | 2018-09-18 |
Family
ID=58009104
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510452356.3A Active CN106405638B (en) | 2015-07-28 | 2015-07-28 | A kind of gas saturation determines method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106405638B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957534A (en) * | 2018-06-12 | 2018-12-07 | 中国石油天然气股份有限公司 | Gas saturation prediction technique and device |
| CN111812709A (en) * | 2020-07-08 | 2020-10-23 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103558629A (en) * | 2013-11-18 | 2014-02-05 | 中国石油天然气集团公司 | Physical simulation method for different air-containing saturation sandstone reservoirs |
| CN103760600A (en) * | 2014-01-07 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas saturation inversion method |
| US20140149042A1 (en) * | 2012-11-23 | 2014-05-29 | Fugro Geoconsulting, Inc. | Method and System for identification of gas hydrates and free gas in geologic beds |
-
2015
- 2015-07-28 CN CN201510452356.3A patent/CN106405638B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140149042A1 (en) * | 2012-11-23 | 2014-05-29 | Fugro Geoconsulting, Inc. | Method and System for identification of gas hydrates and free gas in geologic beds |
| CN103558629A (en) * | 2013-11-18 | 2014-02-05 | 中国石油天然气集团公司 | Physical simulation method for different air-containing saturation sandstone reservoirs |
| CN103760600A (en) * | 2014-01-07 | 2014-04-30 | 中国石油天然气股份有限公司 | Gas saturation inversion method |
Non-Patent Citations (5)
| Title |
|---|
| SARMADI DOOST ,M ET AL.: "AVO analysis for Gas Hydrate Bearing Zone of Oman Sea", 《ISTANBUL INTERNATIONAL GEOPHYSICAL CONFERENCE AND OIL & GAS EXHIBITION》 * |
| 王才志等: "基于地层弹性模量建立含气饱和度模型", 《石油勘探与开发》 * |
| 王磊等: "利用弹性模量计算含气饱和度方法研究", 《西南石油大学学报（自然科学版）》 * |
| 王钧等: "基于斑块饱和模型利用地震波频散特征分析含气饱和度", 《石油地球物理勘探》 * |
| 贺保卫等: "用AVO方法从定性到半定量检测砂岩含气性", 《断块油气田》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957534A (en) * | 2018-06-12 | 2018-12-07 | 中国石油天然气股份有限公司 | Gas saturation prediction technique and device |
| CN108957534B (en) * | 2018-06-12 | 2020-02-14 | 中国石油天然气股份有限公司 | Method and device for predicting gas saturation |
| CN111812709A (en) * | 2020-07-08 | 2020-10-23 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
| CN111812709B (en) * | 2020-07-08 | 2022-03-04 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106405638B (en) | 2018-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200183047A1 (en) | Automated Reservoir Modeling Using Deep Generative Networks | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| Soleimani et al. | Integrated petrophysical modeling for a strongly heterogeneous and fractured reservoir, Sarvak Formation, SW Iran | |
| CN109388817A (en) | A kind of Reservoir Fracture three-dimensional modeling method | |
| Tolstukhin et al. | Ekofisk 4D seismic-seismic history matching workflow | |
| CN105425292A (en) | Oil and gas prediction method and oil and gas prediction device | |
| EP3077618B1 (en) | Tuning digital core analysis to laboratory results | |
| CN106062310A (en) | Facies definition using unsupervised classification procedures | |
| Bokulich | Towards a taxonomy of the model-ladenness of data | |
| CN106168676A (en) | Formation lithology based on seismic data and Fluid Identification Method and device | |
| CN109522578A (en) | Complex Lithofacies prediction technique and system based on Bayes's fuzzy discrimination | |
| CN108663714A (en) | A kind of sedimentary micro depicting method | |
| EP4067944A1 (en) | Systems and methods for identifying geostructural properties as a function of position in a subsurface region of interest | |
| Gilbert et al. | Reservoir modeling: Integrating various data at appropriate scales | |
| WO2016187048A1 (en) | Prospect assessment using structural frameworks | |
| Alameedy et al. | Predicting dynamic shear wave slowness from well logs using machine learning methods in the Mishrif Reservoir, Iraq | |
| CN106405638A (en) | Method and apparatus for determining gas saturation | |
| US9880321B2 (en) | Defining non-linear petrofacies for a reservoir simulation model | |
| Pyrcz et al. | Representative input parameters for geostatistical simulation | |
| US11754736B2 (en) | System and method for classifying seismic data by integrating petrophysical data | |
| Lefranc et al. | Reservoir Characterization for Fracture Distribution Prediction and New Well Target Identification | |
| Aghchelou et al. | Lithofacies Estimation by Multi-Resolution Graph-Based Clustering of Petrophysical Well Logs: Case Study of South Pars Gas Field of Iran | |
| Stephen | Seismic history matching with saturation indicators combined with multiple objective function optimization | |
| Qi et al. | Correlation of seismic attributes and geomechanical properties to the rate of penetration in the Mississippian Limestone, Oklahoma | |
| Macé et al. | Improved Resources and Reserves Estimation with the Combined Use of Scenario-Based Geological Models and Flow-Based Connectivity |
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 |