CN105093313B - A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device - Google Patents
A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device Download PDFInfo
- Publication number
- CN105093313B CN105093313B CN201510395511.2A CN201510395511A CN105093313B CN 105093313 B CN105093313 B CN 105093313B CN 201510395511 A CN201510395511 A CN 201510395511A CN 105093313 B CN105093313 B CN 105093313B
- Authority
- CN
- China
- Prior art keywords
- karst
- geologic
- attribute
- geologic parameter
- parameter
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device, this method includes：Three dimensional seismic data is gathered for Karst-type oil reservoir target area；Target area is divided into multiple geologic elements；M geologic parameter of each sample well is extracted as research unit with geologic element；Selected from M geologic parameter and the N number of geologic parameter of initial production degree of correlation highest preferably matter parameter；Set up the corresponding relation of initial production and important reservoir attribute；It is linear relationship that preferred geologic parameter is set with important reservoir attribute；Set up the non-linear relation equation between initial production data and preferred geologic parameter；The undetermined coefficient in relation equation is determined using Marquardt method, so that it is determined that going out non-linear relation equation；Value according to this relation equation and preferred geologic parameter predicts each sample well capacity.The present invention can accurately predict single well capacity in Karst-type oil reservoir, effectively facilitate the reasonable development of Karst-type oil reservoir.
Description
Technical field
The present invention relates to oil-gas exploration and development technology, in particular it relates to a kind of Karst-type oil reservoir individual well oil and gas productivity prediction
Method and device.
Background technology
About 1/3rd petroleum resources richness is stored in carbonate rock in the world, and in recent years, carbonate rock is gradually
Key areas as oil-gas exploration.With the development of geophysical exploration technology, 3-d seismic exploration technology is used as a kind of hand
Section, plays more and more important effect in oil-gas exploration, and application of the geologic parameter in reservoir prediction is also more and more extensive.
Carbonate Rocks Karst-type reservoir has very strong anisotropism, and this anisotropism leads to not use traditional product
Energy formula carries out capability forecasting, and the productivity prediction model of clastic rock cannot be applied directly in Karst-type reservoir, gives Karst-type storage
Layer prediction causes certain difficulty.
In Karst-type oil reservoir, crack and solution cavity increased the porosity and permeability of reservoir rock, cause Karst-type reservoir
Ooze and change greatly in hole.The inertia resistance of flow of fluid increased the difficulty of oil and gas productivity prediction in the metamorphosis of rock and crack
Degree.When traditional oil gas PRODUCTION FORECASTING METHODS is applied to Karst-type oil reservoir, high-precision forecast difficulty is larger.
Geologic parameter refers to by prestack or post-stack seismic data, by the geometric form derived from mathematic(al) manipulation about seismic wave
State, kinematics character, dynamic characteristic and statistics feature.Based on different application purposes, geologic parameter has different classification
Method.At present, the geologic parameter such as simple extraction amplitude, frequency, relevant, can only qualitative or sxemiquantitative Study In Reservoir information, and
The yield for predicting reservoir that can not be fine.
In sum, under current research conditions, when being predicted to Karst-type oil reservoir oil-production capacity using geologic parameter,
Precision of prediction higher can not be reached.
The content of the invention
It is an object of the invention to overcome above shortcomings in the prior art, and provide a kind of Karst-type oil reservoir individual well
The method and apparatus of initial productivity prediction, so that solve the problems, such as in the prior art cannot accurate prognostic reserves.
The invention provides a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method, described method includes：
Step 1, three dimensional seismic data is gathered for Karst-type oil reservoir target area；Institute is determined according to the three dimensional seismic data
State the interval of interest of target area；
Step 2, extracts coherent body, the dessert attribute of each interval of interest；Determined according to the coherent body and dessert attribute
Fracture river development situation in the target area, multiple is divided into according to the fracture river development situation by the target area
Geologic element；
Step 3, M geologic parameter of each sample well is extracted with the geologic element as research unit；The geologic parameter
Including seismic properties, beading area, bottom hole location away from Dominated Factors distance；
Step 4, the value of initial production and each geologic parameter to each sample well carries out data normalization pretreatment respectively,
Correlation analysis are carried out respectively to being normalized pretreated initial production and the M value of geologic parameter, from described M ground
Selected in matter parameter and the N number of geologic parameter of initial production degree of correlation highest preferably matter parameter；
Step 5, three dense mediums are set up according to the oilwell performance data in the target area, static data and crude oil property data
Matter physical model and yield Mathematical Modeling, according to the treble medium physical model and the yield Mathematical Models just
The corresponding relation of phase yield and important reservoir attribute；The important reservoir attribute includes elastic storativity ratio and interporosity flow coefficient；
Step 6, it is linear relationship to set the preferred geologic parameter with the important reservoir attribute；According to this linear relationship
And the initial production sets up the initial production data with the preferred geology with the corresponding relation of important reservoir attribute
Non-linear relation equation between parameter；
Step 7, the undetermined coefficient in the relation equation is determined using Marquardt method, so that it is determined that going out the non-thread
Sexual intercourse equation；Value according to this relation equation and preferred geologic parameter predicts each sample well capacity.
The above method can also have the characteristics that：
The value of the M is more than or equal to 4；The value of the N is 2,3 or 4.
The above method can also have the characteristics that：
The seismic properties include：Frequency decay percentage, RMS amplitude, amplitude change rate, dessert maximum, dessert
Minimum value, dessert geometrical mean.
The above method can also have the characteristics that：
Correlation analysis in the step 4 include using Pearson came correlation analysis and Si Baiman correlation analysis
Method carries out correlation analysis.
The above method can also have the characteristics that：
The treble medium physical model is the physical model with dielectric-dielectric-Fractured reservoir as standard.
The above method can also have the characteristics that：
Also include carrying out high s/n ratio, high-resolution and high-fidelity to the three dimensional seismic data for collecting in the step 1
Degree treatment.
Present invention also offers a kind of Karst-type oil reservoir individual well oil and gas productivity prediction device, including：
Seismic data acquisition unit, for gathering three dimensional seismic data for Karst-type oil reservoir target area；According to described three
Dimension seismic data determines the interval of interest of the target area；
Seismic data interpretation unit, coherent body, the dessert attribute of each interval of interest are extracted for dividing；According to the phase
Stem body and dessert attribute determine the fracture river development situation in the target area, according to the river development situation that is broken by institute
State target area and be divided into multiple geologic elements；
Geologic parameter extraction unit, the M geology ginseng for extracting each sample well as research unit with the geologic element
Number；The geologic parameter includes seismic properties, beading area, bottom hole location away from Dominated Factors distance；
Geologic parameter select unit, line number is entered for the initial production to each sample well and the value of each geologic parameter respectively
Pre-processed according to normalization, correlation is carried out respectively to being normalized pretreated initial production and the M value of geologic parameter
Analysis, selects with the N number of geologic parameter of initial production degree of correlation highest as preferred from the M geologic parameter
Geologic parameter；
Modeling unit, for being built according to the oilwell performance data in the target area, static data and crude oil property data
Vertical treble medium physical model and yield Mathematical Modeling, build according to the treble medium physical model and the yield Mathematical Modeling
Found the corresponding relation of the initial production and important reservoir attribute；The important reservoir attribute includes elastic storativity ratio and channelling system
Number；
Predicting unit, is linear relationship for setting the preferred geologic parameter with the important reservoir attribute；According to this
Linear relationship and the initial production and the corresponding relation of important reservoir attribute set up the initial production data with it is described
It is preferred that the non-linear relation equation between geologic parameter；The system undetermined in the relation equation is determined using Marquardt method
Number, so that it is determined that going out the non-linear relation equation；Value according to this relation equation and preferred geologic parameter predicts each sample
Well capacity.
Said apparatus can also have the characteristics that：
The value of the M is more than or equal to 4, and the value of the N is 2,3 or 4.
Said apparatus can also have the characteristics that：
The seismic properties include：Frequency decay percentage attribute, RMS amplitude attribute, amplitude change rate attribute, shake
Amplitude variation rate, dessert maximum, dessert minimum value, dessert geometrical mean.
Said apparatus can also have the characteristics that：
The geologic parameter select unit, is also used for Pearson came correlation analysis and Si Baiman correlation analysis
Method carries out correlation analysis.
Compared with prior art, the present invention can accurately predict single well capacity in Karst-type oil reservoir to the present invention, effectively promote
Enter the reasonable development of Karst-type oil reservoir.
Brief description of the drawings
Fig. 1 is Karst-type oil reservoir individual well RECOVERABLE RESERVE PREDICTION method flow diagram.
Specific embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
A part of embodiment of the present invention, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art
The every other embodiment obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.Need
Illustrate, in the case where not conflicting, the feature in embodiment and embodiment in the application can mutually be combined.
There is foregoing description to understand, in Karst-type oil reservoir, crack and solution cavity increased the porosity and permeability of reservoir rock,
Cause Karst-type reservoir hole to be oozed to change greatly.The inertia resistance of flow of fluid increased oil gas in the metamorphosis of rock and crack
The difficulty of capability forecasting.When traditional oil gas PRODUCTION FORECASTING METHODS is applied to Karst-type oil reservoir, high-precision forecast difficulty is larger.It is based on
This, the present invention proposes a kind of method and device of Karst-type oil reservoir individual well oil and gas productivity prediction, preferably geologic parameter, and according to
Relation equation between the initial production data and preferred geologic parameter of foundation, accurately predicts production capacity.
Below in conjunction with accompanying drawing, the present invention is described in detail.
Embodiment one
The embodiment of the present invention provides a kind of method of Karst-type oil reservoir individual well oil and gas productivity prediction, as shown in figure 1, the method
Including：
Step 1, three dimensional seismic data is gathered for Karst-type oil reservoir target area；Mesh is determined according to above-mentioned three dimensional seismic data
Mark the interval of interest in area.
Step 2, frequency dividing extracts coherent body, the dessert attribute of each interval of interest；Determined according to coherent body and dessert attribute
Fracture river development situation in target area, multiple geologic elements are divided into according to this fracture river development situation by target area.
Step 3, M geologic parameter of each sample well is extracted with geologic element as research unit；Geologic parameter includes earthquake
Attribute, beading area, bottom hole location are away from Dominated Factors distance；
Step 4, the value to the initial production of sample well and each geologic parameter carries out data normalization pretreatment respectively, right
Be normalized pretreated initial production and the M value of geologic parameter carries out correlation analysis respectively, from M geologic parameter
In select and the N number of geologic parameter of initial production degree of correlation highest preferably matter parameter.
Step 5, treble medium thing is set up according to the oilwell performance data in target area, static data and crude oil property data
Reason model and yield Mathematical Modeling, according to above-mentioned treble medium physical model and the initial production of yield Mathematical Models with it is important
The corresponding relation of reservoir attribute；Above-mentioned important reservoir attribute includes elastic storativity ratio and interporosity flow coefficient.
Step 6, it is linear relationship to set preferred geologic parameter with important reservoir attribute；According to this linear relationship and initial stage
The corresponding relation of yield and important reservoir attribute sets up the relation equation between initial production data and preferred geologic parameter.
Step 7, the undetermined coefficient in above-mentioned relation equation is determined using Marquardt method, so that it is determined that going out relation side
Journey；Value according to this relation equation and preferred geologic parameter predicts sample well capacity.
The following detailed description of each step of the above method.
In the step of above method 1, the method for gathering three dimensional seismic data includes：Explosive is buried by measurement, brill shallow well hole
(when using dynamite source), bury wave detector, arrangement cable to several procedures receptions of instrument cab and transmitting seismic signal to obtain
Take three dimensional seismic data.After collection three dimensional seismic data, high s/n ratio, high-resolution also are carried out to the three dimensional seismic data for collecting
Rate and high fidelity are processed.
Determining the method for the interval of interest of target area in step 1 according to three dimensional seismic data includes：It is soft using seismic interpretation
Part, seismic profile is shown using wiggle variable area Faxian, carries out seismic interpretation；Using well logging sound wave, density curve and seismic data
It is combined, makes composite traces, completes the demarcation to earthquake interval.
Explanation on interval includes herein below：
By browse line, road, etc. when section etc. seismic data, understand work area in earthquake reflected wave reflectance signature, mainly
Target zone earth formation, the time-domain section of local structure, plane configuration determine the concrete thought and method of Fine structural interpretation.
Continuity according to earthquake reflected wave carries out comparative interpretation respectively, the back wave poor for continuity, according to by
Principle small, from coarse to fine is arrived greatly, since the well profile excessively demarcated, it is established that the general explanation framework of the whole district, for continuous
Property preferable back wave, using three-dimensional visualization explain seed point method for tracing, carry out automatic tracing contrast.
In order to strengthen the uniformity that earthquake reflected wave follows the trail of contrast, when manual comparative interpretation is carried out, as far as possible using automatic
Tracer technique, it is ensured that the feature such as crest, trough, zero phase that seismic wave is followed the trail of is consistent in the whole district, is well-shooting joint inversion,
Reservoir parameter is extracted and reservoir prediction, there is provided accurate seismic interpretation layer position.
The achievement of comparative interpretation, is checked and is contrasted using three-dimensional visualization means.
In the step of above method 2, frequency dividing extracts coherent body, the dessert attribute of each interval of interest, observes in the plane
The distribution situation of tomography, river course.According to the distribution situation of river course, fracture in plane, the geologic element of target area is divided.
Coherent technique is a kind of method for predicting similitude, can export similitude between seismic channel, structure dip and
The attributes such as azimuth.The change of similitude is often relevant with tomography, the equal geological phenomenon of deposition, and continuous geologic body is generally corresponded to
Similar value higher, exceptional value occurs when stratum is discontinuous.3-D seismics amplitude data body is converted into coefficient correlation number
According to body, unexpected incoherent anomaly.
Dessert is obtained by instantaneous frequency and reflected intensity (also referred to as instantaneous amplitude or amplitude envelope).From number
Dessert is the ratio of reflected intensity and root mean square instantaneous frequency from the perspective of.
Reflected intensity, also known as instantaneous amplitude, amplitude envelope, is the square root of a certain moment seismic signal gross energy.Its value is total
It is positive, is the amplitude information unrelated with phase, can provide the information of acoustic impedance difference is used for the attributional analysis of amplitude anomaly, uses
Reflected to detect tomography, river course, subterranean deposit, thin layer tuning effect and differentiate thick layer from complex wave.Reflected intensity is
Identification and effective instrument to bright spot and dim spot, the cross directional variations of reflected intensity generally with the aggregation of lithology and hydrocarbon
Relevant, the acute variation of reflected intensity is also relevant with the deposition characteristicses change such as river course.Its mathematical definition is：
Wherein, A (t) is the seismic reflection intensity level of moment t；F (t) is the real signal of seismic reflection signals, and h (t) is f
The Hilbert conversion of (t), be in the computational methods of time-domain：H (t)=f (t)/π t.
Instantaneous frequency is rate of change of the instantaneous phase to the time, and instantaneous phase is empty seismic channel and real seismic trace ratio
Arc tangent.Computing formula is：
Or
Its codomain is (- f ,+f).But, most of instantaneous phases are all for just.Instantaneous frequency can provide the effective of lineups
The information of frequency absorption effect and crack influence and reservoir thickness.Identification, determination deposition for air accumulation band and low-frequency band
The jumping phenomenons such as thickness, display pinching, hydrocarbon water termination.
Seismic volume has high amplitude concurrently and low-frequency part enables to dessert property value bigger, and its of both attributes
He combines can be so that dessert property value diminishes.
In the step of above method 3, the value of M is more than or equal to 4.Wherein, seismic properties include：Frequency decay percentage,
RMS amplitude, amplitude change rate, dessert maximum, dessert minimum value, dessert geometrical mean.
Wherein, the loss of the frequency decay gross energy that refers to seismic wave in underground medium propagation, be in medium attribute
The factor for causing seismic wave frequency spectrum to be decayed, is the energy in medium between solid and solid, solid and fluid, fluid and fluid boundary
Amount consume.Theoretical research and practical application show, if pores'growth in geologic body, and when filling oil, gas, seismic reflection is inhaled
Receive and increase, high-frequency absorption decay aggravation, low frequency energy increases.
RMS amplitude is that the average value of Amplitude-squared extracts square root again.Because amplitude is before average square, therefore,
It is very sensitive to king-sized amplitude.It is suitable for the analysis of karst reservoir solution cavity, is also used for formation lithology Phase Transition Analysis, calculates thin
Sand layer thickness, identification bright spot, dim spot indicate hydro carbons display feature.Its computing formula is
Wherein, aiIt is the amplitude of ith sample point.
Amplitude change rate is similar with RMS amplitude, can highlight in work area amplitude value mutation in stratum, is capable of identify that
The size and scale of solution cavity in carbonate reservoir, are the representative attributes for recognizing karst Reservoir Body.
In the step of above method 4, the value of the value less than M of N.By many experiments and checking, the value of N be preferably 2,3 or
4.Specific correlation analysis include related to Si Baiman (Spearman) using Pearson came (Pearson) correlation analysis
Property analysis method carry out correlation analysis, obtain N number of preferred geologic parameter.
In statistics, Pearson product moment correlation coefficients be for measuring two correlations between variable X and Y,
Span is between [- 1 ,+1].Pearson product-moment correlation coefficient is widely used to measure two variables in academic research
The power of linear dependence.
Pearson correlation coefficient is a kind of situation, because necessarily assuming that data in Pearson correlation coefficient calculating process
It is to be obtained from normal distribution in couples, the requirement for data is higher, therefore we employ Spearman order phases simultaneously
Relation number is analyzed.Spearman rank correlation coefficients are the order statistical parameters of a nonparametric property (unrelated with distribution), only
There is dull functional relation in X and Y, then X is exactly complete Spearman related to Y.
In the step of above method 5, step 6 and step 7, according to the oilwell performance data in target area, static data and
Crude oil property data sets up treble medium physical model.This model is that basement rock sillar is divided into according to its porosity and permeability
Two classes：It is connective good between one class and Fracture System；It is another kind of then poor, using both media as two independent liquid
Supply source, fluid flows into crack from the two independent pore medias respectively, is flowing to shaft bottom.Can be with according to this kind of sorting technique
Fracture medium containing isolated cave is summarized as the treble medium oil pool of a class matrix-matrix-slit formation.
When setting up the Mathematical Modeling of treble medium oil pool, three kinds of media meet the respective equation of motion, state equation respectively
And continuity equation, and the channelling between medium is represented with a source sink term in continuity equation.With Ge Jiali et al.
The Mathematical Modeling of treble medium is set up as a example by the matrix-matrix-Fractured reservoir of proposition.
Because crack is the passage of oil-gas migration, only consider two kinds of matrix as oil reservoirs in prediction.For crack body
System：
Wherein, K3It is the permeability in crack, its unit is md；μ is fluid viscosity, its unit for μm2；P3For crack is original
Strata pressure, its unit is MPa；φ is porosity；C is that well stores up coefficient, and its unit is m3/MPa；Split for basement rock 1,2 is flowed into
Seam flow.
For basement rock 1：
Wherein, K1It is the permeability of basement rock 1, its unit is md；P1It is the original formation pressure of basement rock 1, its unit is MPa.
For basement rock 2：
Wherein, K2It is the permeability of basement rock 2, its unit is md；P2It is the original formation pressure of basement rock 2, its unit is MPa.
Assuming that basement rock permeability is relatively low, ignore left right-hand vector, abbreviation can be obtained
In the case of quasi-stable state channelling,WithBe given by following formula：
α1And α2It is to fix value coefficient；
Define dimensionless variable：
Wherein r is the seepage flow radius in percolation equationk, rwIt is wellbore radius；tDIt is non dimensional time.
Zero dimension is carried out to equation (4) (5) (6), after arrangement：
rDIt is non dimensional time,
Wherein,
Zero dimension boundary condition：
Zero dimension primary condition：
PDjRefer to the zero dimension pressure in j-th network.
To tDLaplace changes are done to get in return,
pDj(rD,s)|S=0=0 (j=1,2,3；1≤rD≤+∞)
S is Laplace operators.
Arrange,
Solution procedure according to the rank of empty argument 0 and 1 rank Bessel functions can be in the hope of：
And oil well output qDProduct according to the pressure solution under the yield solution and fixed output quota under Duhamel's principle, i.e. level pressure is
Yield：
Therefore
The formula in Laplace spaces is transformed into the real space using the inverting of Stehfest algorithms, inversion principle is：
V (i) refers to i-th fluid-flow rate of time step.
The relation equation q=f between yield and matrix 1, the elastic storativity ratio of matrix 2, interporosity flow coefficient is obtained by inverting
(ω1,ω2,λ1,λ2), the elastic storativity ratio of matrix 1 is ω1, interporosity flow coefficient is λ1；The elastic storativity ratio of matrix 2 is ω2, channelling
Coefficient is λ2。
It is linear relationship that preferred geologic parameter is set with important reservoir attribute：
ω1=a1x1+a2x2+…+anxn
ω2=b1x1+b2x2+…+bnxn
λ1=c1x1+c2x2+…+cnxn
λ2=d1x1+d2x2+…+dnxn
Wherein x1,x2…xnIt is preferred geologic parameter.ω1It is the elastic storativity ratio of matrix 1, λ1It is the channelling system of matrix 1
Number；ω2It is the elastic storativity ratio of matrix 2, λ2It is the interporosity flow coefficient of matrix 2.
Four equations bring q=f (ω into it will be assumed1,ω2,λ1,λ2), obtain between initial production and preferred geologic parameter
Nonlinear equation：Q=f (x1,x2…xn,a1,a2…an,b1,b2…bn,c1,c2…cn,d1,d2,…dn)；；a1, a2…an,
b1, b2…bn, c1, c2…cn, d1, d2…dnIt is the undetermined coefficient in equation.
According to the geologic parameter and corresponding single well productivity that extract, under least square meaning, can with Marquardt method
The undetermined coefficient in nonlinear equation is determined, so that it is determined that nonlinear multivariable equation of the production capacity on geologic parameter.According to this pass
It is that the value of equation and preferred geologic parameter predicts each sample well capacity.
Embodiment two
The invention provides a kind of device of Karst-type oil reservoir individual well oil and gas productivity prediction, this device includes：
Seismic data acquisition unit, for gathering three dimensional seismic data for Karst-type oil reservoir target area；According to dimensionally
Shake data determines the interval of interest of target area；
Seismic data interpretation unit, coherent body, the dessert attribute of each interval of interest are extracted for dividing；According to coherent body
The fracture river development situation in target area is determined with dessert attribute, is divided into target area according to fracture river development situation many
Individual geologic element；
Geologic parameter extraction unit, the M geologic parameter for extracting each sample well as research unit with geologic element；Ground
Matter parameter includes seismic properties, beading area, bottom hole location away from Dominated Factors distance；
Geologic parameter select unit, line number is entered for the initial production to each sample well and the value of each geologic parameter respectively
Pre-processed according to normalization, correlation is carried out respectively to being normalized pretreated initial production and the M value of geologic parameter
Analysis, selected from M geologic parameter with the N number of geologic parameter of initial production degree of correlation highest preferably matter ginseng
Number；
Modeling unit, for setting up three according to the oilwell performance data in target area, static data and crude oil property data
Dense media physical model and yield Mathematical Modeling, according to treble medium physical model and the initial production of yield Mathematical Models with
The corresponding relation of important reservoir attribute；Important reservoir attribute includes elastic storativity ratio and interporosity flow coefficient；
Predicting unit, is linear relationship for setting preferred geologic parameter with important reservoir attribute；According to this linear relationship
And the corresponding relation of initial production and important reservoir attribute set up it is non-between initial production data and preferred geologic parameter
Linear relationship equation；The undetermined coefficient in relation equation is determined using Marquardt method, so that it is determined that going out non-linear relation side
Journey；Value according to this relation equation and preferred geologic parameter predicts each sample well capacity.
Wherein, value of the value of M more than or equal to 4, N is 2,3 or 4.
Seismic properties include：Frequency decay percentage attribute, RMS amplitude attribute, amplitude change rate attribute, amplitude become
Rate, dessert maximum, dessert minimum value, dessert geometrical mean.
Geologic parameter select unit, is also used for Pearson came correlation analysis and Si Baiman correlation analysis
Carry out correlation analysis.
The function of each unit is corresponding respectively with the technical characteristic in the above method in the present apparatus, and here is omitted.
Concrete application example
Comprised the following steps in this application example：
1st, 48 mouthfuls of samples mouthful wells are had in western certain work area, western certain work area by measurement, explosive is buried in brill shallow well hole (makes
During with dynamite source), bury that wave detector, arrangement cable are received to several procedures of instrument cab and transmitting seismic signal obtains karst
The three dimensional seismic data of type oil reservoir target area.
2nd, the treatment of high s/n ratio, high-resolution and high fidelity is carried out to the three dimensional seismic data for collecting.
3rd, in seismic interpretation software, seismic profile is shown using wiggle variable area Faxian, carries out seismic interpretation.Using well logging
Sound wave, density curve are combined with seismic data, make composite traces, and it is to determine target zone to complete to the demarcation of seismic horizon
Section.
4th, frequency dividing extracts coherent body, the dessert attribute of each interval of interest, according to coherent body and dessert attribute in the plane
Observation fracture, the distribution situation in river course, geologic element is further divided according to fracture, river course distribution situation.
5th, with geologic element as research unit, 6 values of geologic parameter at 48 mouthfuls of sample well shaft bottoms, this 6 ground are extracted
Matter parameter includes：RMS amplitude attribute, interval frequency decay percentage attribute, beading area, amplitude change rate, dessert are minimum
Value, dessert maximum.
6th, using Pearson correlation analysis and Spearman correlation analysis to the geologic parameter value extracted
Correlation analysis are carried out with initial production data, the correlation coefficient charts for obtaining are as follows：
7th, according to the actual conditions and basic data of each geologic element in work area, treble medium yield model is set up.
8th, the initial productivity data and geologic parameter data to 48 mouthfuls of wells are corrected and standardization.
9th, the data after the input standardization of mode input end, record the prediction data of output after forecast model convergence
And the undetermined coefficient value in equation, the result being calculated is as follows：
| Undetermined coefficient | Solution value | Undetermined coefficient | Solution value |
| a1 | 0.26 | c1 | -0.03 |
| a2 | 0.63 | c2 | -0.10 |
| a3 | 0.64 | c3 | -0.10 |
| b1 | -1.24 | d1 | -0.23 |
| b2 | -2.74 | d2 | -15.30 |
| b3 | 1.65 | d3 | 4.88 |
By actual contrast, predicting the outcome, it is preferable to be coincide with record of production data, and consensus forecast precision reaches 86.8%.
The present invention can accurately predict single well capacity in Karst-type oil reservoir, effectively facilitate rationally opening for Karst-type oil reservoir
Hair.
Furthermore, it is necessary to explanation, the specific embodiment described in this specification, is named the shape of its parts and components
Claiming etc. can be with difference, and the above content described in this specification is only to structure example explanation of the present invention.
Descriptions above can combine implementation individually or in a variety of ways, and these variants all exist
Within protection scope of the present invention.
Herein, term " including ", "comprising" or any other variant thereof is intended to cover non-exclusive inclusion, from
And cause to include the article or equipment of a series of key elements not only including those key elements, but also its including being not expressly set out
His key element, or it is this article or the intrinsic key element of equipment also to include.In the absence of more restrictions, by language
The key element that sentence " including ... " is limited, it is not excluded that also exist in the article or equipment including the key element other identical
Key element.
The above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted, reference only to preferred embodiment to this hair
It is bright to be described in detail.It will be understood by those within the art that, technical scheme can be modified
Or equivalent, without deviating from the spirit and scope of technical solution of the present invention, all should cover in claim model of the invention
In the middle of enclosing.
Claims (10)
1. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method, it is characterised in that described method includes：
Step 1, three dimensional seismic data is gathered for Karst-type oil reservoir target area；The mesh is determined according to the three dimensional seismic data
Mark the interval of interest in area；
Step 2, extracts coherent body, the dessert attribute of each interval of interest；According to the coherent body and dessert attribute determine
Fracture river development situation in target area, multiple geology are divided into according to the fracture river development situation by the target area
Unit；
Step 3, M geologic parameter of each sample well is extracted with the geologic element as research unit；The geologic parameter includes
Seismic properties, beading area, bottom hole location are away from Dominated Factors distance；
Step 4, the value of initial production and each geologic parameter to each sample well carries out data normalization pretreatment respectively, to entering
Row normalizes pretreated initial production and the M value of geologic parameter and carries out correlation analysis respectively, from M geology ginseng
Selected in number and the N number of geologic parameter of initial production degree of correlation highest preferably matter parameter；
Step 5, treble medium thing is set up according to the oilwell performance data in the target area, static data and crude oil property data
Reason model and yield Mathematical Modeling, the initial stage according to the treble medium physical model and the yield Mathematical Models produce
The corresponding relation of amount and important reservoir attribute；The important reservoir attribute includes elastic storativity ratio and interporosity flow coefficient；
Step 6, it is linear relationship to set the preferred geologic parameter with the important reservoir attribute；According to this linear relationship and
The initial production sets up the initial production data with the preferred geologic parameter with the corresponding relation of important reservoir attribute
Between non-linear relation equation；
Step 7, the undetermined coefficient in the relation equation is determined using Marquardt method, so that it is determined that going out the nonlinear dependence
It is equation；Value according to this relation equation and preferred geologic parameter predicts each sample well capacity.
2. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method according to claim 1, it is characterised in that
The value of the M is more than or equal to 4；The value of the N is 2,3 or 4.
3. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method according to claim 1, it is characterised in that describedly
Shake attribute includes：Frequency decay percentage, RMS amplitude, amplitude change rate, dessert maximum, dessert minimum value, dessert are several
What average value.
4. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method according to claim 1, it is characterised in that the step
Correlation analysis in rapid 4 include carrying out correlation using Pearson came correlation analysis and Si Baiman correlation analysis
Analysis.
5. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method according to claim 1, it is characterised in that described three
Dense media physical model is the physical model with dielectric-dielectric-Fractured reservoir as standard.
6. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction method according to claim 1, it is characterised in that the step
Also include carrying out the three dimensional seismic data for collecting the treatment of high s/n ratio, high-resolution and high fidelity in rapid 1.
7. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction device, it is characterised in that including：
Seismic data acquisition unit, for gathering three dimensional seismic data for Karst-type oil reservoir target area；According to it is described dimensionally
Shake data determines the interval of interest of the target area；
Seismic data interpretation unit, coherent body, the dessert attribute of each interval of interest are extracted for dividing；According to the coherent body
The fracture river development situation in the target area is determined with dessert attribute, according to the river development situation that is broken by the mesh
Mark zoning is divided into multiple geologic elements；
Geologic parameter extraction unit, the M geologic parameter for extracting each sample well as research unit with the geologic element；Institute
Geologic parameter is stated including seismic properties, beading area, bottom hole location away from Dominated Factors distance；
Geologic parameter select unit, carries out data and returns respectively for the initial production to each sample well and the value of each geologic parameter
One changes pretreatment, and correlation analysis are carried out respectively to being normalized pretreated initial production and the M value of geologic parameter,
Selected from the M geologic parameter and the N number of geologic parameter of initial production degree of correlation highest preferably matter
Parameter；
Modeling unit, for setting up three according to the oilwell performance data in the target area, static data and crude oil property data
Dense media physical model and yield Mathematical Modeling, according to the treble medium physical model and the yield Mathematical Models institute
State the corresponding relation of initial production and important reservoir attribute；The important reservoir attribute includes elastic storativity ratio and interporosity flow coefficient；
Predicting unit, is linear relationship for setting the preferred geologic parameter with the important reservoir attribute；It is linear according to this
It is preferred with described that relation and the initial production set up the initial production data with the corresponding relation of important reservoir attribute
Non-linear relation equation between geologic parameter；The undetermined coefficient in the relation equation is determined using Marquardt method, from
And determine the non-linear relation equation；Value according to this relation equation and preferred geologic parameter predicts the product of each sample well
Energy.
8. as claimed in claim 7 a kind of Karst-type oil reservoir individual well oil and gas productivity prediction device, it is characterised in that the value of the M
More than or equal to 4, the value of the N is 2,3 or 4.
9. a kind of Karst-type oil reservoir individual well oil and gas productivity prediction device as claimed in claim 7, it is characterised in that earthquake category
Property includes：Frequency decay percentage attribute, RMS amplitude attribute, amplitude change rate attribute, amplitude change rate, dessert are maximum
Value, dessert minimum value, dessert geometrical mean.
10. as claimed in claim 7 a kind of Karst-type oil reservoir individual well oil and gas productivity prediction device, it is characterised in that the geology
Parameter selection unit, being also used for Pearson came correlation analysis and Si Baiman correlation analysis carries out correlation point
Analysis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510395511.2A CN105093313B (en) | 2015-07-07 | 2015-07-07 | A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510395511.2A CN105093313B (en) | 2015-07-07 | 2015-07-07 | A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105093313A CN105093313A (en) | 2015-11-25 |
| CN105093313B true CN105093313B (en) | 2017-07-07 |
Family
ID=54574167
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510395511.2A Expired - Fee Related CN105093313B (en) | 2015-07-07 | 2015-07-07 | A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105093313B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105700018A (en) * | 2016-03-31 | 2016-06-22 | 中国石油天然气集团公司 | Earthquake attribute optimizing method and earthquake attribute optimizing device |
| CN106503834A (en) * | 2016-09-30 | 2017-03-15 | 中国石油天然气股份有限公司 | Prediction method for lake-phase ultra-low pore permeability reservoir compact oil dessert area |
| CN107944620B (en) * | 2017-11-21 | 2021-11-09 | 西南石油大学 | Nonlinear prediction method for single-well steady-state production performance |
| CN109975189B (en) * | 2017-12-28 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for predicting productivity of pore type sandstone reservoir |
| CN110297264B (en) * | 2018-03-23 | 2021-01-01 | 中国石油化工股份有限公司 | Low-permeability gas reservoir thin reservoir dessert earthquake prediction method |
| CN109162693B (en) * | 2018-09-17 | 2020-06-02 | 中国地质大学（北京） | Method for rapidly testing rock mass block index by using monitoring while drilling technology without coring |
| CN110593865B (en) * | 2019-09-29 | 2022-07-29 | 中国石油集团川庆钻探工程有限公司 | Well testing interpretation method for characteristic parameters of oil reservoir fracture hole |
| CN112948513B (en) * | 2019-12-11 | 2024-03-26 | 中国石油天然气股份有限公司 | Method, device and storage medium for generating energy distribution trend graph |
| CN113627639A (en) * | 2020-05-07 | 2021-11-09 | 中国石油化工股份有限公司 | Well testing productivity prediction method and system for carbonate fracture-cave reservoir |
| CN113627069A (en) * | 2020-05-08 | 2021-11-09 | 中国石油化工股份有限公司 | Well testing dynamic yield evaluation method and system for fracture-cavity type oil reservoir oil and gas well |
| CN112346121B (en) * | 2020-11-04 | 2023-09-12 | 东北石油大学 | Reservoir stratum separation treatment method based on full waveform |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008028122A2 (en) * | 2006-09-01 | 2008-03-06 | Chevron U.S.A. Inc. | History matching and forecasting in the production of hydrocarbons |
| US8095349B2 (en) * | 2008-05-30 | 2012-01-10 | Kelkar And Associates, Inc. | Dynamic updating of simulation models |
| CN103410502B (en) * | 2013-08-05 | 2016-03-16 | 西南石油大学 | A kind of acquisition methods of three-dimensional permeability fields of network-like fracture-pore reservoir |
| CN104695950B (en) * | 2014-10-31 | 2017-10-17 | 中国石油集团西部钻探工程有限公司 | Volcanic Reservoir PRODUCTION FORECASTING METHODS |
-
2015
- 2015-07-07 CN CN201510395511.2A patent/CN105093313B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105093313A (en) | 2015-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105093313B (en) | A kind of Karst-type oil reservoir individual well oil and gas productivity prediction method and device | |
| CN101738639B (en) | Method for improving computing precision of rock fracture parameters | |
| CN104879103B (en) | Layered water injection effect analysis method | |
| CN103454685B (en) | Method and device for predicting sand body thickness by utilizing logging constrained wave impedance inversion | |
| CN103135135B (en) | Method and device for quantitative prediction of hydrocarbons based on loose sandstone model | |
| CN109613612A (en) | A kind of carbonate rock particle beach meticulous depiction and prediction technique | |
| CN109425896A (en) | Dolomite oil and gas reservoir distribution prediction method and device | |
| CN102736107B (en) | Energy constraint heterogeneous reservoir thickness identification system | |
| CN105089615B (en) | A kind of log data historical regression processing method based on reservoir model | |
| CN105386751B (en) | A kind of horizontal wellbore logging PRODUCTION FORECASTING METHODS based on reservoir model | |
| CN104965979A (en) | Tight sandstone effective reservoir identifying method | |
| CN112746837A (en) | Shale oil reservoir exploration data acquisition system and method based on distributed optical fiber sensing | |
| CN104863574B (en) | A kind of Fluid Identification Method suitable for tight sandstone reservoir | |
| CN106597543B (en) | Stratum sedimentary facies division method | |
| CN103914620B (en) | Method for computing distribution of opening spaces of fractures of fault rupture zones | |
| CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
| CN105445800A (en) | Thick sand body top differentiation lithologic reservoir identification method | |
| CN103775057A (en) | Method and device for identifying effective reservoir of tight oil and gas reservoir | |
| CN103439740B (en) | Method and device for predicting relative impedance based on dipole seismic wavelet multiple integral | |
| CN106246158B (en) | Method and device for distributing wells in ultra-deep low-hole fractured sandstone gas reservoir | |
| CN103643949A (en) | Quantitatively forecasting method and device for oil-gas possibility of reservoirs | |
| CN105863628A (en) | Fine reservoir prediction method of oilfield development phase | |
| CN105938503A (en) | Multi-layer interface recognition method of direction signals | |
| CN105116449A (en) | Method for identifying weak reflection reservoir | |
| CN113821956B (en) | Evaluation method for disturbance quantity of current geostress structure of deep shale reservoir |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170707 Termination date: 20190707 |